Mental disorders are characterised by a complex-genetic aetiology, where – assuming the traditional ‘vulnerability–stress model’ – biological factors interact with negative environmental influences to shape the risk towards the respective disorder.Reference Zubin and Spring1 Translating this concept to current methodology, research focusing on gene–environment (G × E) interactions has paid much attention to the interactive detrimental effects of genetic risk factors and environmental adversity (for review see, for example, Musci et alReference Musci, Augustinavicius and Volk2). The ability to successfully cope in the face of adversity, i.e. to ‘successfully evoke the capacities of the organism to contend with problems’,Reference Zubin and Spring1 constitutes another key factor in shaping the risk-resilience spectrum towards mental disorders.Reference Rutten, Hammels, Geschwind, Menne-Lothmann, Pishva and Schruers3 However, apart from a few and mostly theoretical publications in schizophrenia research,Reference Zubin, Steinhauer and Condray4 coping has not been included in the general conception of the vulnerability–stress model and has not yet been systematically studied in synopsis with genetic and adverse environmental factors in mental disorders. Consequently, there is a need for probing a novel conceptual framework integrating both adverse events and coping factors in an extended ‘vulnerability–stress–coping model’ of mental disorders. In the present study, a three-dimensional G × E × C model expanded by the dimension of coping ability (‘C’) was therefore exemplarily applied to the phenotype of anxiety based on the following state of knowledge.
The neuropeptide S system
On a genetic level (‘G’), the neuropeptide S (NPS) system constitutes one of the most promising biological candidates related to arousal and anxiety. The gene coding for the NPS receptor (NPSR1) is located on chromosome 7p14 and contains an A/T (Asn107Ile) single nucleotide polymorphism (rs324981), with the T allele resulting in a tenfold increase in NPS potency at the NPS receptor.Reference Reinscheid, Xu, Okamura, Zeng, Chung and Pai5 This more active T allele has been linked to panic disorder per se Reference Domschke, Reif, Weber, Richter, Hohoff and Ohrmann6 as well as to panic disorder endophenotypes such as elevated levels of anxiety sensitivity,Reference Domschke, Reif, Weber, Richter, Hohoff and Ohrmann6,Reference Beste, Konrad, Uhlmann, Arolt, Zwanzger and Domschke7 enhanced response inhibition and increased error monitoring,Reference Beste, Konrad, Uhlmann, Arolt, Zwanzger and Domschke7 increased heart rate and higher symptom reports during a behavioural avoidance test,Reference Domschke, Reif, Weber, Richter, Hohoff and Ohrmann6 and enhanced neuroendocrine and subjective responses to acute stress.Reference Kumsta, Chen, Pape and Heinrichs8 On a neural level, T allele carriers exhibited increased amygdala activation in response to fearful and angry faces in healthy probands,Reference Dannlowski, Kugel, Franke, Stuhrmann, Hohoff and Zwanzger9 whereas in patients with panic disorder the T allele appears to be related to attenuated activation of the dorsolateral prefrontal, lateral orbitofrontal and anterior cingulate cortex during the processing of fearful faces,Reference Domschke, Reif, Weber, Richter, Hohoff and Ohrmann6 potentially based on a delayed maturation of corticolimbic connectivity during adolescence.Reference Domschke, Akhrif, Romanos, Bajer, Mainusch and Winkelmann10 Taken together, the NPSR1 rs324981 variant has been shown to play a significant role in anxiety and anxiety disorders, particularly panic disorder and is involved in the modulation of a variety of intermediate anxiety phenotypes.
Environment and coping
On an environmental level (‘E’), experiences of threat, interpersonal conflicts or stressful life events in adulthood as well as traumatic experiences in the sensitive period of childhood have been reported to constitute risk factors for the later development of anxiety disorders.Reference Scott, McLaughlin, Smith and Ellis11
Regarding the dimension of coping ability (‘C’), self-efficacy – referring to an individual's belief in his or her own ability to successfully cope with challenge and measured as a general construct termed general self-efficacyReference Bandura12 – constitutes a promising dimension related to coping with adversity. High general self-efficacy has been linked to lower levels of trait anxiety in healthy adolescents and adults,Reference Endler, Speer, Johnson and Flett13,Reference Muris14 less severe symptoms of depression,Reference Muris14,Reference Maciejewski, Prigerson and Mazure15 decreased risk for symptoms of social anxiety in childhoodReference Rudy, Davis and Matthews16 and resilience against the deleterious effects of daily stressors,Reference Schönfeld, Brailovskaia, Bieda, Zhang and Margraf17 thus highlighting its function as an important link between environmental adversity and anxiety-related mental health outcomes.
Gene–environment interactions
In first G × E studies, the interactive effects of NPSR1 rs324981 genotype and adverse environmental aspects on anxiety-related phenotypes were described for childhood maltreatment,Reference Klauke, Deckert, Zwanzger, Baumann, Arolt and Pauli18 family adversityReference Laas, Reif, Akkermann, Kiive, Domschke and Lesch19 and urban upbringing.Reference Streit, Haddad, Paul, Frank, Schafer and Nikitopoulos20 However, results differ regarding the allelic direction of this interaction: the TT genotype has been shown to drive anxiety risk dependent on experiences of childhood traumaReference Klauke, Deckert, Zwanzger, Baumann, Arolt and Pauli18 and to confer enhanced amygdala activation during stress exposure in the context of high early urbanicity.Reference Streit, Haddad, Paul, Frank, Schafer and Nikitopoulos20 By contrast, Laas et alReference Laas, Reif, Akkermann, Kiive, Domschke and Lesch19 observed an interaction of the AA genotype with a history of stressful life events conferring increased trait anxiety as well as with family adversity resulting in a higher risk for anxiety disorders. No study so far has looked into beneficial factors such as coping ability possibly further shaping the spectrum of risk and resilience towards anxiety and thus potentially reconciling previous diverging findings.
Against this background, the present study investigated the interactive effect of the functional NPSR1 rs324981 A/T variant, childhood trauma and general self-efficacy on trait anxiety in a large discovery and an independent replication sample of healthy volunteers applying a G × E × C approach in order to probe an extended ‘vulnerability–stress–coping model’ of mental disorders.
Method
Samples and procedures
The discovery sample included a total of 1403 healthy participants (women n = 826, mean age 25.59 years, s.d. = 5.71) and was recruited in the context of project Z02 within the Collaborative Research Centre SFB-TRR58 ‘Fear, Anxiety, Anxiety Disorders’ during the project's second funding period at the Universities of Würzburg, Hamburg and Münster, Germany, between 2013 and 2016 (see Schiele et alReference Schiele, Ziegler, Holitschke, Schartner, Schmidt and Weber21).
The replication sample, consisting of 630 healthy adults (women n = 445, mean age 25.03 years, s.d. = 5.57), was independently recruited during the above-mentioned project's third funding period at the Universities of Würzburg, Hamburg and Münster, Germany, between 2016 and 2018.
For both samples, inclusion criteria were defined as ethnic White descent (self-report up to third generation), age at inclusion between 18 and 50 years, right-handedness (see below) and fluency in German. Probands with a past or current DSM-IV axis I disorder as ascertained by experienced psychologists (Mini-International Psychiatric Interview),Reference Sheehan, Lecrubier, Sheehan, Amorim, Janavs and Weiller22 past or current severe neurological or somatic disorders, currently taking centrally active medication, excessive alcohol (≥15 units/week), nicotine (≥20 cigarettes/day) and caffeine (≥4 cups/day) consumption, utilisation of illegal drugs, or pregnancy were excluded. Participants completed a set of questionnaires (see below). For genetic analyses, a venous blood sample (~18 mL) was taken. The project has been designed to also provide a recruitment pipeline channelling comprehensively pheno- and genotyped participants into several subprojects within the SFB-TRR58 for specific research questions including neuroimaging studies (hence left-handedness was defined as an exclusion criterion).
Participation was remunerated with 50€. Written informed consent was obtained from all participants. The study was reviewed and approved by the ethical committees of the Universities of Würzburg, Hamburg and Münster and was conducted in compliance with the Declaration of Helsinki.
Self-report measures
The short form of the Childhood Trauma Questionnaire (CTQReference Bernstein and Fink23), containing 28 items, retrospectively assesses the frequency of childhood maltreatment. The total score is obtained as the sum of all items, resulting in possible sum scores between 25 and 128.
The General Self-Efficacy Scale (GSEReference Schwarzer, Jerusalem, Weinman, Wright and Johnston24) comprises ten items to assess perceived self-efficacy, i.e. a person's belief in their own ability to cope with difficulties, which can be summed up, resulting in total scores between 10 and 40 points.
The trait scale of the State-Trait Anxiety Inventory (STAI-TReference Spielberger, Gorsuch and Lushene25) consists of 20 items capturing anxiety as a stable, trans-situational disposition. Seven items are reverse-coded to minimise errors because of arbitrary responding. After re-coding of the respective items, the total STAI-T score is calculated as the sum of all items, resulting in possible scores between 20 and 80.
Genotyping
For genotyping of NPSR1 rs324981, DNA extracted from EDTA (ethylenediamine tetraacetic acid)-blood was amplified by polymerase chain reaction (PCR, 45 s at 95°C, 45 s at 58°C, 45 s at 72°C for 35 cycles) using oligonucleotide primers F: 5ʹ-TGCTTTGCATTTCCTCAGTG-3ʹ and R: 5ʹ-TTGTCTCATCACATTTGGAAGG-3ʹ. PCR products were digested with AseI at 37°C overnight, separated on 3% agarose gel containing ethidium bromide, and visualised by ultraviolet light.
Hardy–Weinberg criteria as determined by the online program DeFinetti (http://ihg.gsf.de/cgi-bin/hw/hwa1.pl) were fulfilled for both samples (discovery sample: TT, 306; AT, 710; AA, 387; P = 0.592; replication sample: TT, 131; AT, 328; AA 171; P = 0.264). For further analyses, NPSR1 rs324981 genotypes were grouped into TT genotype versus A allele (AA/AT) carriers in accordance with previous studies in healthy probands.Reference Klauke, Deckert, Zwanzger, Baumann, Arolt and Pauli18
Statistical analyses
Group differences regarding continuous variables were analysed by means of t-tests. Differences regarding categorical variables were tested by means of χ2-tests. Possibly confounding gene–environment correlations were evaluated using Pearson's correlations. The influence of NPSR1 genotype, CTQ and GSE scores, as well as their interactions on STAI-T scores were tested via hierarchical multiple regression analyses (cf. Schiele et alReference Schiele, Ziegler, Holitschke, Schartner, Schmidt and Weber21). Genotype groups were centredReference Kraemer and Blasey26 and coded as 0.5 (TT) and −0.5 (AA/AT). CTQ and GSE sum scores were centred (mean 0, s.d. = 1). Regression analyses were performed in three steps: in the first step, main effects were entered into the model, i.e. genotype, centred CTQ sum score and centred GSE sum score. Gender (centred; 0.5, women; −0.5, men) was entered as an additional factor in the first step. In a second step, all two-way interaction terms – genotype × CTQ, genotype × GSE and CTQ × GSE – were included. Finally, in a third step, the three-way interaction term comprising genotype × CTQ × GSE was added. All statistical tests were performed using SPSS V24 (SPSS Inc, Chicago, Illinois, USA). The significance level was set at 0.05.
Results
Descriptives
Descriptive characteristics of the discovery and the replication sample are given in Table 1.
Table 1 Descriptive sample characteristics
CTQ, Childhood Trauma Questionnaire, GSE, General Self-Efficacy Scale; STAI-T: Trait scale of the State-Trait Anxiety Inventory; NPSR1: neuropeptide S receptor gene.
Genotype groups did not differ with regard to age (discovery sample: t = 1.13, d.f. = 1401, P = 0.258; replication sample: t = 1.72, d.f. = 628, P = 0.125), gender (discovery sample: χ2 = 0.08, d.f. = 1, P = 0.777; replication sample: χ2 = 0.93, d.f. = 1, P = 0.335), CTQ score (discovery sample: t = −0.33, d.f. = 1401, P = 0.745; replication sample: t = −0.79, d.f. = 628, P = 0.431), GSE score (discovery sample: t = −1.45, d.f. = 1401, P = 0.149; replication sample: t = −0.50, d.f. = 628, P = 0.619), or STAI-T score (discovery sample: t = −0.15, d.f. = 1401, P = 0.878; replication sample: t = 0.85, d.f. = 628, P = 0.395).
No significant gene–environment correlations with either predictor variable were observed (discovery sample: all P ≥ 0.149; replication sample: all P ≥ 0.431).
Discovery sample
Step 1 yielded significant main effects of CTQ score (β = 0.214, t = 9.960, P = 1.265 × 10−22) and GSE score (β = −0.543, t = −25.148, P = 1.942 × 10−115). In step 2, significant main effects emerged for CTQ score (β = 0.214, t = 8.213, P = 4.976 × 10−16) and GSE score (β = −0.566, t = −21.136, P = 3.048 × 10−86). Step 3 returned significant main effects of CTQ score (β = 0.250, t = 9.405, P = 2.058 × 10−20) and GSE score (β = −0.583, t = −21.849, P = 3.031 × 10−91), significant two-way interactions of GSE × CTQ scores (β = 0.086, t = 3.669, P = 2.250 × 10−4 and genotype × GSE score (β = −0.056, t = −2.128, P = 0.034) as well as a significant three-way interaction term of genotype × CTQ × GSE score (β = 0.129, t = 5.524, P = 3.938 × 10−8). The addition of the three-way interaction term in step 3 accounted for a significant increment in explained variance (R 2 = 0.380, ΔR 2 = 0.014, ΔF = 30.520, P = 3.938 × 10−8).
The relationship between NPSR1 genotype, CTQ and GSE on STAI-T scores is depicted in Fig. 1(a); in A allele carriers, high general self-efficacy exerted a buffering effect on trait anxiety despite a history of childhood trauma, whereas low general self-efficacy led to increased anxiety scores in individuals with an otherwise equal genetic and environmental risk constellation. In turn, increases in trait anxiety as a function of childhood maltreatment were observed in TT homozygotes largely unaffected by general self-efficacy.
Fig. 1 Effect of childhood trauma on anxiety as a function of grouped NPSR1 rs324981 genotype and general self-efficacy.
(a) In the discovery sample of 1403 healthy probands, and (b) in the replication sample of 630 healthy probands (illustration of linear regression analyses).
CTQ, Childhood Trauma Questionnaire; GSE: General Self-Efficacy Scale; STAI-T, trait scale of the State-Trait Anxiety Inventory; neuropeptide S receptor (NPSR1) rs324981 genotype groups (TT genotype versus A allele carriers).
Replication sample
Significant main effects of CTQ (β = 0.251, t = 7.661, P = 7.041 × 10−14) and GSE score (β = −0.511, t = −15.467, P = 6.235 × 10−46) emerged in step 1. Step 2 returned significant main effects of CTQ (β = 0.224, t = 5.177, P = 3.043 × 10−7) and GSE score (β = −0.520, t = −11.552, P = 4.230 × 10−28). In step 3, significant main effects of CTQ (β = 0.260, t = 5.677, P = 2.109 × 10−8) and GSE score (β = −0.518, t = −11.546, P = 4.519 × 10−28) were obtained, as well as a significant three-way interaction term of genotype × CTQ × GSE (β = 0.102, t = 2.328, P = 0.020). The addition of the three-way interaction term in step 3 accounted for a significant increment in explained variance (R 2 = 0.341, ΔR 2 = 0.006, ΔF = 5.421, P = 0.020) (see Fig. 1(b)). Again, low or high GSE score, respectively, went along with increased or decreased trait anxiety depending on CTQ scores in A allele carriers, but not in TT homozygotes.
Discussion
Main findings
The present G × E × C study demonstrates a moderating influence of the coping factor ‘general self-efficacy’ on the deleterious effects of childhood maltreatment on trait anxiety in a genotype-dependent fashion. Specifically, the highest trait anxiety scores were observed in carriers of the less active NPSR1 rs324981 A allele with a history of childhood maltreatment and low general self-efficacy. However, when general self-efficacy was high, this pattern was reversed, with A allele carriers then showing the lowest trait anxiety scores despite the experience of childhood adversity. In individuals homozygous for the more active NPSR1 rs324981 T risk allele, increases in trait anxiety as a function of childhood adversity were observed to be largely unaffected by general self-efficacy. These findings suggest that in interaction with childhood adversity the TT risk genotype – irrespective of the influence of general self-efficacy – increases the risk for anxiety, whereas the A allele does not, unless accompanied by low general self-efficacy, i.e. insufficient coping ability. Therefore, it can be proposed that high general self-efficacy is able to alleviate trait anxiety in A allele carriers despite a history of maltreatment, whereas low general self-efficacy can lead to increased anxiety scores in individuals with an otherwise equal genetic and environmental risk constellation.
Contextual embedding
The present results are in line with the existing G × E literature linking the NPSR1 rs324981 TT genotype to increased anxiety sensitivityReference Klauke, Deckert, Zwanzger, Baumann, Arolt and Pauli18 or heightened amygdala activationReference Streit, Haddad, Paul, Frank, Schafer and Nikitopoulos20 via environmental variation. Furthermore, they are in accordance with previous findings of the T allele mediating panic disorder as a nosological entityReference Domschke, Reif, Weber, Richter, Hohoff and Ohrmann6 as well as intermediate phenotypes of panic disorder such as dimensional anxiety, psychophysiological responding,Reference Domschke, Reif, Weber, Richter, Hohoff and Ohrmann6,Reference Beste, Konrad, Uhlmann, Arolt, Zwanzger and Domschke7 altered neural activation patternsReference Domschke, Reif, Weber, Richter, Hohoff and Ohrmann6,Reference Dannlowski, Kugel, Franke, Stuhrmann, Hohoff and Zwanzger9,Reference Domschke, Akhrif, Romanos, Bajer, Mainusch and Winkelmann10 or neuroendocrine stress responsiveness,Reference Kumsta, Chen, Pape and Heinrichs8 whereas the A allele was generally not found to predispose to anxiety disorder risk per se or to influence anxiety traits on an intermediate level. The present G × E × C finding might furthermore reconcile an apparently discrepant previous observation of the A allele being associated with a higher frequency of anxiety disorders in women reporting a negative family environment and with higher trait anxiety dependent upon the experience of past stressful life events:Reference Laas, Reif, Akkermann, Kiive, Domschke and Lesch19 coping characteristics could function as an additional and previously unconsidered dimension able to buffer or, if maladaptive, further increase disease risk. Thus, in line with the ‘differential susceptibility hypothesis’,Reference Belsky, Jonassaint, Pluess, Stanton, Brummett and Williams27 rather than representing a ‘risk gene’, NPSR1 appears to constitute a ‘plasticity gene’ insofar as its rs324981 A allele appears to drive increased susceptibility to environmental conditions, positive as well as negative, whereas the TT genotype may be less responsive to the resilience-increasing effects of adaptive coping capabilities.
Implications
Conceptually, the present study suggests revisiting the original concept proposed by Zubin & SpringReference Zubin and Spring1 and to officially extend the general understanding of the ‘vulnerability–stress model’ in mental disorders by the dimension of coping towards a three-dimensional ‘vulnerability–stress–coping’ or G × E × C model, where in the face of environmental adversity, impaired coping ability lowers the threshold towards illness despite a low genetic vulnerability, whereas high coping ability can buffer the impact of a high genetic and environmental factor constellation and thus increase the threshold towards illness (see Fig. 2).
Fig. 2 An extended three-dimensional ‘vulnerability–stress–coping model’ of mental disorders based on the vulnerability–stress model by Zubin & Spring.1
The light blue area indicates illness, the dark blue area indicates wellness. Black lines indicate thresholds between illness and wellness. Coping ability can either decrease the threshold towards illness (solid black line) if low, or increase it (dotted black line) if high. The dark blue/light blue hatched areas indicate wellness or illness depending on an individual's constellation of vulnerability factors, stress experiences and coping ability.
From a clinical point of view, in contrast to an individual's genetic makeup and past adverse experiences, coping ability, i.e. general self-efficacy, is dynamic and malleable, and may therefore constitute a target for indicated preventive interventions in at-risk individuals increasing the threshold between wellness und illness. Indeed, general self-efficacy has been shown to increase following stress management training in university studentsReference Molla Jafar, Salabifard, Mousavi and Sobhani28 and self-management interventions in severe mental disordersReference Lean, Fornells-Ambrojo, Milton, Lloyd-Evans, Harrison-Stewart and Yesufu-Udechuku29 as well as over the course of cognitive–behavioural therapy.Reference Gallagher, Payne, White, Shear, Woods and Gorman30 General self-efficacy was furthermore shown to confer clinical symptom improvement in panic disorderReference Gallagher, Payne, White, Shear, Woods and Gorman30 and social anxiety disorder.Reference Bouchard, Gauthier, Nouwen, Ivers, Vallieres and Simard31,Reference Gaudiano and Herbert32
Limitations and future directions
The present study holds several strengths, such as investigating a functional genetic variant in large, sufficiently powered and independent samples of comparable demographic composition, under application of identical inclusion/exclusion criteria and using the same psychometric instruments and analysis methods under identical laboratory conditions, thus allowing for robust replication of the main finding. Nonetheless, some limitations should be taken into account while interpreting the present results and when designing future studies. The presently assessed samples comprised a relatively young, mostly student population of healthy volunteers with a mean age of 25.6 and 25.0 years, and may thus not be representative of the general population. Therefore, replication of the presently reported results is warranted interrogating more heterogeneous samples of healthy probands as well as extending the presently applied G × E × C approach to a clinical context by including high-risk individuals and patients with manifest anxiety disorders. In addition, in the present study childhood maltreatment was assessed retrospectively, which may be sensitive to recall bias.
The reported rates of experiences of childhood maltreatment were relatively low as indicated by low CTQ sum scores.
Extending the present cross-sectional design, longitudinal studies could provide promising insights into the developmental trajectories of anxiety from susceptibility to disorder manifestation. This approach could additionally contribute to the understanding of resilient functioning given evidence for the differential importance of environmental contributions during childhood, adolescence and adulthood, with the influence of environmental factors on the phenotypic stability of symptoms of anxiety and depression increasing with advancing age.Reference Nivard, Dolan, Kendler, Kan, Willemsen and van Beijsterveldt33 Finally, on a genetic level, beyond the single candidate-gene approach future research may want to address the G × E × C model under consideration of haplotype or epistatic genetic effects as well as in the context of genome-wide association studies in sufficiently powered samples. This is because, in particular, recent genome-wide studies have reported several loci to significantly contribute to coping and resilience phenotypes.Reference Maul, Giegling, Fabbri, Corponi, Serretti and Rujescu34,Reference Stein, Choi, Jain, Campbell-Sills, Chen and Gelernter35 Along these lines, future studies will have to disentangle the individual and interactive genetic effects on risk and resilience, respectively, given that resilience towards stress-related mental disorders has been reported to show a heritability of up to 40%Reference Waaktaar and Torgersen36 and to display distinct biological underpinnings which do not simply constitute a reversal of mechanisms leading to disease risk.Reference Choi, Stein, Dunn, Koenen and Smoller37,Reference Feder, Fred-Torres, Southwick and Charney38 In addition, given increasing evidence for a role of temporally dynamic epigenetic processes in the modulation of gene function partly mediating adaptation to environmental influences in the context of anxiety disorders and anxiety-related phenotypes,Reference Schiele and Domschke39 future G × E × C models would benefit from the inclusion of epigenetic markers such as DNA methylation towards the formulation of a comprehensive multidimensional model of anxiety risk.
Summary and conclusion
In conclusion, the present G × E × C study exemplarily suggests that adaptive coping ability – i.e. high levels of general self-efficacy – may act as a broker between an otherwise detrimental gene–environment vulnerability profile and the resulting phenotype by exerting a protective effect in a resilience-enhancing way. In detail, a history of childhood maltreatment may result in increased trait anxiety and thus in an elevated risk for the development of anxiety disorders in NSPR1 A allele carriers, but only in the absence of the ability to successfully cope with adversity, whereas NSPR1 TT homozygotes carry an increased risk of anxiety traits as a function of maltreatment experiences irrespective of coping ability. The present results provide a framework for an extended ‘vulnerability–stress–coping model’ of anxiety and possibly mental disorders in general and – along the notion of a ‘precision medicine’ approach – might allow for targeted preventive interventions for at-risk individuals carrying a G × E risk factor constellation by strengthening their adaptive coping strategies.
Data availability
The data associated with the manuscript is not currently publicly available.
Acknowledgements
We gratefully acknowledge the skilful technical support by C. Gagel as well as the support in proband recruitment (discovery sample) of D. Schümann, L. Kürbitz and K. Holitschke supervised by T.B.L. and P.Z. K.D., M.A.S. and P.Z. are members of the Anxiety Disorders Research Network (ADRN), European College of Neuropsychopharmacology (ECNP).
Author contributions
M.A.S. performed the statistical analyses and wrote the first version of the manuscript. K.D. designed the study and drafted the manuscript together with M.A.S. The discovery sample was coordinated and in part recruited by M.A.S. under the supervision of K.D., J.D., A.R. and P.P. The replication sample was recruited by K.H., E.J.L., J.R., J.B. and K.R. under the supervision of K.D., J.D., U.D., T.B.L., P.P. and M.A.S. Genotyping analyses were performed by C.S. (discovery sample) and L.K. (replication sample) supervised by K.D., J.D. and A.R. Part of this work was included in the doctoral thesis of M.A.S. Authors M.A.S. and K.D. have full and ongoing access to all the data in the study and take responsibility for the integrity of the data and the accuracy of the data analysis. All authors have read and approved the final article.
Funding
This work was funded by the German Research Foundation (DFG) – project number 44541416 – TRR 58, subprojects Z02 (to K.D., A.R., P.P., J.D., T.B.L., U.D.) and C02 (to K.D. and J.D.) and CRC 1193, subproject Z03 (to A.R.).
Declaration of interest
K.D., A.R., P.P., J.D., T.B.L. and U.D. report funding by the German Research Foundation (DFG) in relation to this study. J.D. reports grants from the European Union and the Bavarian Secretary of Commerce outside the submitted work.
An ICMJE form is in the supplementary material, available online at https://doi.org/10.1192/bjp.2020.73.
Open access
eLetters
No eLetters have been published for this article.