Association results of NTRK2 with suicide attempt among depressed patients. A, The negative common logarithm of the best genetic model P value of the 65 tagging single-nucleotide polymorphisms (SNPs) obtained in the discovery sample was plotted against the SNP's chromosome position throughout the NTRK2 locus on chromosome 9q22. The lower horizontal line indicates the nominal significance level of .05; the upper horizontal line, the significance level after correction for multiple comparisons over all tested SNPs and genetic models; asterisks, the 5 SNPs for which association with suicide attempt could be replicated in a larger independent German sample; and Mb, megabases. B, The 5′ to 3′ structures of the most prominent isoforms of NTRK2 are depicted and aligned with part A. The boxes and lines show the exon-intron structure of NTRK2. TrkB indicates receptor tyrosine kinase B. C, The r2-based linkage disequilibrium structure over all genotyped tagging SNPs in control subjects of the discovery sample is shown. The genetic information content of NTRK2 spanning across 355 kilobases (kb) can be clustered into 16 D′-based haplotype blocks. D, The linkage disequilibrium (r2) between the 5 SNPs associated with suicide attempt among depressed patients in both German samples is displayed.
Multilocus regression analysis revealed a single-nucleotide polymorphism (SNP)–SNP interaction in suicide attempt (SA). Only the 3 SNPs most strongly associated with SA among depressed patients in the combined sample were included in multilocus modeling to conserve power. In addition to the SNPs' main effects on risk for SA (Table 3), the independent SNPs (Figure 1D) rs10868235 and rs1147198 also showed an interaction effect on SA. The relative risk for SA is plotted against all genotype combinations of the 2 markers.
Replication of suicide attempt (SA) association in an African American cross-sectional study of nonpsychiatric clinic patients. The 2 tested single-nucleotide polymorphisms, rs1147198 (A) and rs1867283 (B), showed no association with SA in the complete sample (all) but did show association with SA among the subgroup of subjects with an ascertained lifetime diagnosis of major depressive disorder (MDD) for rs1147198 (A) and with trend significance also for rs1867283 (B). The risk alleles were the same as in the German samples (Table 3). OR indicates odds ratio; CI, confidence interval.
Kohli MA, Salyakina D, Pfennig A, Lucae S, Horstmann S, Menke A, Kloiber S, Hennings J, Bradley BB, Ressler KJ, Uhr M, Müller-Myhsok B, Holsboer F, Binder EB. Association of Genetic Variants in the Neurotrophic Receptor–Encoding Gene NTRK2 and a Lifetime History of Suicide Attempts in Depressed Patients. Arch Gen Psychiatry. 2010;67(4):348-359. doi:10.1001/archgenpsychiatry.2009.201
Copyright 2010 American Medical Association. All Rights Reserved. Applicable FARS/DFARS Restrictions Apply to Government Use.2010
A consistent body of evidence supports a role of reduced neurotrophic signaling in the pathophysiology of major depressive disorder (MDD) and suicidal behavior. Especially in suicide victims, lower postmortem brain messenger RNA and protein levels of neurotrophins and their receptors have been reported.
To determine whether the brain-derived neurotrophic factor (BDNF) gene or its high-affinity receptor gene, receptor tyrosine kinase 2 (NTRK2), confer risk for suicide attempt (SA) and MDD by investigating common genetic variants in these loci.
Eighty-three tagging single-nucleotide polymorphisms (SNPs) covering the genetic variability of these loci in European populations were assessed in a case-control association design.
Inpatients and screened control subjects.
The discovery sample consisted of 394 depressed patients, of whom 113 had SA, and 366 matched healthy control subjects. The replication studies comprised 744 German patients with MDD and 921 African American nonpsychiatric clinic patients, of whom 152 and 119 were positive for SA, respectively.
Blood or saliva samples were collected from each participant for DNA extraction and genotyping.
Main Outcome Measures
Associations of SNPs in BDNF and NTRK2 with SA and MDD.
Independent SNPs within NTRK2 were associated with SA among depressed patients of the discovery sample that could be confirmed in both the German and African American replication samples. Multilocus interaction analysis revealed that single SNP associations within this locus contribute to the risk of SA in a multiplicative and interactive fashion (P = 4.7 × 10−7 for a 3-SNP model in the combined German sample). The effect size was 4.5 (95% confidence interval, 2.1-9.8) when patients carrying risk genotypes in all 3 markers were compared with those without any of the 3 risk genotypes.
Our results suggest that a combination of several independent risk alleles within the NTRK2 locus is associated with SA in depressed patients, further supporting a role of neurotrophins in the pathophysiology of suicide.Published online February 1, 2010 (doi:10.1001/archgenpsychiatry.2009.201).
Suicide represents a major public health problem. Each year, an estimated 10 to 20 million suicides are attempted and 1 million suicides are completed worldwide.1 The overall lifetime prevalence of suicide attempts (SAs) ranges between 3% and 5%, with similar numbers across different populations.2,3 Patients with psychiatric disorders are at higher risk for SA, with patients with mood disorder showing the highest relative risk for lifetime SA (odds ratio [OR] = 7.8-29.9), followed by patients with other psychiatric diagnoses (OR = 2.1-6.5).3 This increased risk is also reflected by reported estimates for lifetime prevalence of completed suicide being higher in patients with mood disorder than in the general population (2%-4% vs 0.5%, respectively).4 Twin and family studies suggest that suicide and SA are heritable traits5,6 and likely part of the same phenotype, with completed suicide and SA clustering in the same families.7- 11 The genetic risk factors for suicide appear to be independent from the underlying psychiatric disorder.7,9,12- 14
Regarding the neurobiological foundation of suicidal behavior, a compelling body of evidence points to alterations in the monoamine, neurotrophin, and stress response systems.15- 19 Monoaminergic and stress hormone system dysregulations have already received support from human genetic association studies.20 On the other hand, the genetic effect of the neurotrophic system has not yet been extensively investigated, even though a series of postmortem studies now convincingly shows a downregulation of neurotrophic signaling in suicide victims.21- 24 The neurotrophin family is composed of nerve growth factor, brain-derived neurotrophic factor (BDNF), neurotrophin 3, neurotrophin 4/5, their high-affinity receptors (receptor tyrosine kinase A [NTRK1], receptor tyrosine kinase B [NTRK2], and receptor tyrosine kinase C [NTRK3]), and the low-affinity nerve growth factor receptor (p75NTR).25 Postmortem brain studies uniformly report reduced BDNF and NTRK2 messenger RNA and protein expression in the prefrontal cortex and hippocampus of suicide victims with depression or other psychiatric disorders relative to nonpsychiatric and nonsuicidal control subjects. Moreover, suicidal subjects showed reduced expression of nerve growth factor, neurotrophin 3, neurotrophin 4/5, NTRK1, and NTRK3 and reduced phosphorylation of all of these receptor tyrosine kinases in the hippocampus.21- 24,26 This suggests insufficient neurotrophic signaling and thus synaptic plasticity in suicidal subjects, which is consistent with results of other studies showing altered volumes of the hippocampus and prefrontal cortex27- 30 and reduced plasma31 and serum BDNF levels in suicidal subjects.32
Genetic studies investigating the genes encoding the neurotrophic system in suicidal behavior and major depressive disorder (MDD) have mostly focused on a functional BDNF Val66Met (rs6265) polymorphism with inconsistent results33- 36 or polymorphisms within the low-affinity receptor gene NGFR37,38 (for review, see the “Comment” section). To our knowledge, there are no published case-control studies in SA or MDD exploring the BDNF gene and the gene encoding its high-affinity receptor, NTRK2, with comprehensive genetic coverage. Here we report the results of a study investigating the association of 83 tagging single-nucleotide polymorphisms (SNPs) for these 2 genes with SA and MDD in a German sample with replication in a larger independent German case-control sample and an African American cohort. To avoid confounding by the underlying psychiatric disease, analyses testing association with SA were performed within the patient group.
A total of 405 inpatients with mood disorder (178 men, 227 women) were recruited for the Munich Antidepressant Response Signature study (http://www.mars-depression.de) of the Max Planck Institute of Psychiatry, Munich, Germany. The study has been described previously in detail.39- 41 Briefly, diagnosis was ascertained by trained psychiatrists according to DSM-IV criteria. Patients classified as having an MDD by 2 independent psychiatrists and fulfilling the criteria for at least a moderate depressive episode (score ≥14 on the 21-item Hamilton Scale for Depression) were included. Table 1 and Table 2 show demographic and illness-related data. A lifetime history of SA was ascertained using a semistructured clinical interview established for the Munich Antidepressant Response Signature study39 and scoring points in the Hamilton Scale for Depression rating scale (score on suicide item = 4). Information on lifetime history of SA was available in 394 patients (97.3%). As controls, 366 healthy subjects were randomly selected from a Munich-based community sample and matched to patients for age, sex, and ethnicity (Table 1). The healthy control subjects were screened for the absence of any DSM-IV Axis I psychiatric disorder using the face-to-face computer-assisted interview based on the Munich version of the Composite International Diagnostic Interview.42
A total of 921 patients (302 men, 619 women) with recurrent major depression were recruited at the Max Planck Institute of Psychiatry and psychiatric hospitals in Augsburg and Ingolstadt (both in Bavaria, Germany). For details on this sample, see the article by Lucae et al.43 Patients were diagnosed by World Health Organization–certified raters according to DSM-IV criteria using the World Health Organization Schedules for Clinical Assessment in Neuropsychiatry, version 2.1. Patients of European descent with at least 2 moderately severe depressive episodes were included. A lifetime history of SA was ascertained using the same semistructured clinical interview established for the Munich Antidepressant Response Signature study.39 Information on lifetime history of SA was available for 744 patients (80.8%). A total of 1029 healthy control subjects matched to the patient sample for age, sex, and ethnicity were randomly selected from a Munich-based community sample. They were free of lifetime anxiety and mood disorders according to the Munich version of the Composite International Diagnostic Interview.42
Both studies have been approved by the ethics committee of Ludwig-Maximilians University of Munich, and written informed consent was obtained from all subjects.
From a cross-sectional study of nonpsychiatric clinic patients with significant levels of trauma, 921 African American patients (390 men, 531 women) with information on lifetime history of SA were selected (119 with SA and 802 without SA). Study design, ascertainment, and rating protocols have been described elsewhere in more detail.44 The patients' mean (SD) age was 40.8 (14.1) years (men, 44.3 [12.7] years; women, 38.3 [14.5] years). Participants were seeking general medical or obstetric/gynecological care at an urban public hospital in Atlanta, Georgia. The Structured Clinical Interview for DSM-IV Axis I Disorders (nonpatient edition)45 was performed for 335 of these probands, of whom 80 fulfilled criteria for a lifetime history of MDD (32 with SA and 48 without SA; 32 men and 48 women; mean [SD] age for all 80 subjects, 45.5 [10.5] years; mean [SD] age for men, 49.3 [8.0] years; mean [SD] age for women, 42.9 [11.3] years). All procedures in this study were approved by the institutional review boards of Emory University School of Medicine and Grady Memorial Hospital.
On enrollment in the study, up to 40 mL of blood was drawn from each patient into anticoagulant EDTA tubes. The DNA was extracted from fresh blood using standard DNA extraction procedures with the Puregene whole-blood DNA extraction kit (Gentra Systems Inc, Minneapolis, Minnesota).
An SNP marker panel was used to tag the full-length NTRK2 (RefSeq NM_006180, 355.04 kilobases [kb]) and the longest BDNF isoform (RefSeq NM_170731, 66.86 kb) comprising 20 kb of flanking sequences according to Human HapMap Project phase I and II data for the central European population.46 The tagger software integrated in the HapMap project home page47 was used for tagging SNP selection. Pairwise r2 was set to ≥0.8 with a minor allele frequency (MAF) cutoff of ≥0.1 for NTRK2 and ≥0.01 for BDNF. The suggested number of tagging SNPs was 69 for NTRK2 and 18 for BDNF (total of 87 tagging SNPs). The SNPs were genotyped on a matrix-assisted laser desorption/ionization–time of flight mass spectrometer (MassARRAY system; Sequenom, Inc, San Diego, California) using the Assay Designer software (Sequenom, Inc) for primer selection, multiplexing, and assay design and using the homogeneous mass extension process for producing primer extension products. Matrix-assisted laser desorption/ionization–time of flight SNP genotyping was performed at the facilities of the Helmholtz Zentrum Munich.
For power calculations, the Genetic Power Calculator was used48 (http://pngu.mgh.harvard.edu/~purcell/gpc). Given a unipolar depression prevalence of 16%,49 a common marker in perfect linkage disequilibrium (LD) (r2 = 1) with the causative variation under an allelic, dominant, or recessive model, and 80% statistical power, we were able to detect an effect size of 1.7 or higher in our discovery sample (α = .05 divided by 83 SNPs equals 6.02 × 10−4). The prevalence of a lifetime history of SA among depressed patients was estimated to be 23% from our combined German samples, which is consistent with the reported prevalence.3 Under the earlier-mentioned settings, we had 80% statistical power to detect effect sizes of 2.1 or higher on SA among patients with depression in our discovery sample.
Case-control association screening in the discovery sample was conducted using WG-Permer software (http://www.mpipsykl.mpg.de/wg-permer/). The SNPs were tested for association with MDD in a case-control design using 2-tailed χ2 test statistics under allelic, recessive, and dominant modes of inheritance. Genetic associations with lifetime history of SA among depressed patients were tested by 2-tailed analysis of variance with residual variables that were generated by logistic regression analysis correcting for age and sex. For this analysis, patients with and without history of SA were compared to avoid confounding by the underlying psychiatric disease. Obtained association P values were verified by calculating permutation-based empirical P values under 105 permutations (WG-Permer software) and by logistic regression analysis including sex and age as covariates using SPSS statistical software version 16 for Windows (SPSS Inc, Chicago, Illinois). Empirical P values and P values from logistic regression analysis with covariates did not differ from those of the residual analysis of variance test statistics of the screening analysis. The level of significance was set to .05 for P values corrected for multiple comparisons by the permutation-based minimum P method50,51 under 105 permutations over all performed tests (3 genetic models, 83 SNPs) per phenotype (MDD, SA). In addition, the permutation-based Fisher product method (FPM)52 was performed, which allows association testing over all tested markers of a candidate gene with the phenotype of interest. The FPM was exerted under 105 permutations over all tested SNPs and genetic models for both phenotypes (MDD, SA) and genes (BDNF, NTRK2). The 4 resulting FPM P values were Bonferroni corrected for multiple comparisons (α = .05 divided by 4 equals .0125). The same statistical procedures were applied for association testing in the German replication sample but under the SNP-wise best genetic model and with the same risk allele/genotype as observed in the discovery sample, which allowed 1-sided testing. In a second replication sample (African American probands), an allelic model with the same risk allele as determined in the German samples was tested (1-sided). Multilocus interaction testing on SA among depressed patients of the combined German sample was performed using stepwise logistic regression analysis conducted in R statistical software version 2.5.0 (http://cran.r-project.org). We included the 3 most significant SNPs (under their best genetic models), their pairwise interaction term, age, and sex in the model. Meta-analyses were computed by applying Comprehensive Meta-analysis statistical software version 2 (Biostat, Englewood, New Jersey).
The LD pattern and haplotype block delineation were determined using Haploview software version 4.0 (http://www.broad.mit.edu/mpg/haploview) in control subjects of the discovery sample.53 Blocks were defined using the method described by Gabriel et al.54
In the discovery sample, a total of 87 tagging SNPs covering NTRK2 and BDNF were genotyped. Data of 4 single r2-bin SNPs of NTRK2 were excluded from association testing because of insufficient call rates. The remaining 83 SNPs were included in association analysis. Their mean (SD) call rate was 99.7% (0.82%), and no SNP displayed deviation from Hardy-Weinberg equilibrium after Bonferroni correction for multiple testing (eTable 1 and eTable 2). Overall, we obtained a 100% SNP marker coverage of BDNF with 20 kb of flanking sequences according to HapMap phase I and II data (central European population, r2 ≥ 0.8, MAF ≥ 0.01). For NTRK2, SNP marker coverage was 94.2% (central European population, r2 ≥ 0.8, MAF ≥ 0.1).
We did not detect any case-control association with MDD in NTRK2 or in BDNF, which withstood correction for multiple comparison. The smallest P value (nominal P = 9.3 × 10−3, corrected P = .38) was obtained for rs9969765, an SNP located in the NTRK2 gene. There were 7 SNPs in NTRK2 and 4 SNPs in BDNF, which showed nominal P values less than .05 (eTable 1 and eTable 2). For instance, the often reported functional BDNF polymorphism (rs6265) showed nominal significant association with MDD (P = .035). Neither NTRK2 or BDNF were associated with MDD by applying the FPM overall-tested variants within a gene.
There was 1 SNP in the promoter region of NTRK2 that was associated with SA among depressed patients after correction for multiple comparison over all tested SNPs and genetic models (rs11140714: nominal P = 2.6 × 10−4, corrected P = .043). A total of 12 SNPs located in NTRK2 were associated with SA on a nominal significance level of .05 (Figure 1 and eTable 2), whereas in BDNF there were none (eTable 1). The FPM suggested association of NTRK2 with SA among depressed patients over all tested SNPs and models with P = 4.1 × 10−3. After correction for all 4 FPM tests performed in this study (2 phenotypes × 2 genes), NTRK2 was still significantly associated with SA among depressed patients (Bonferroni-corrected P = .016). In contrast to NTRK2, BDNF was not significantly associated with SA by the FPM.
To follow-up the association between NTRK2 and SA among depressed patients in the discovery sample, we selected the 12 nominally significant NTRK2 SNPs to be genotyped in the independent German replication sample (eTable 3). We performed the FPM on SA with the best genetic model and the same risk allele/genotype per SNP as have been observed in the discovery sample. The FPM P value was significant in the replication sample (P = 3.3 × 10−3), indicating real association of SNPs with SA among depressed patients. Five single SNP associations (rs10868235, rs1147198, rs1867283, rs1187286, and rs11140800) among the 12 tested SNPs showed nominal significance with the same risk allele/genotype and genetic model as in the discovery sample (Table 3). This number of nominally replicated single SNP associations is about 16.7 times higher than would be expected by chance.
When depressed patients with SA from both studies were tested vs healthy control subjects for the 5 associated SNPs mentioned earlier, significant associations were again observed but with effect sizes smaller than when compared with depressed patients with no SA (Table 3). These 5 SNPs did not show case-control association with MDD in the combined sample, indicating a specific genetic effect on SA. Of the replicated SNPs, only rs2277193 became nominally significantly associated with MDD in the combined sample (n = 1250/1375, allelic P = .033), but this association did not withstand correction for multiple testing.
Age at onset of depressive disorder was 5.6 years lower in those with SA compared with those without SA in the combined German sample (P = 5.1 × 10−8) (Table 2). This result is consistent with a previous study.55 In addition, a higher number of previous depressive episodes in patients with a lifetime history of SA compared with patients without SA has been reported.55 In our samples, we observed the same direction of this effect, approaching significance in the combined sample (P = .061) (Table 2). We thus controlled for the age at onset of depressive disorder and the number of previous depressive episodes for the described genetic associations with SA among depressed patients, in addition to age and sex, in post hoc logistic regression analysis. In these models, all 5 SNPs remained significant predictors for SA, with sex and age at onset as significant covariates. We repeated these analyses with other potentially confounding variables showing at least trend significance between SA-positive and SA-negative patients in at least 1 sample (alcohol and nicotine consumption, remission) (Table 2). Only alcohol consumption became another significant covariate, but neither alcohol consumption nor nicotine consumption affected the reported genetic associations with SA.
The LD structure of tested NTRK2 SNPs is shown in Figure 1C and D. By definition, tagging SNPs are not in high LD with each other. Among the 5 SNPs associated with SA, rs10868235 and rs11140800 correlated most (r2 = 0.63; D′ = 0.94), while rs1147198 was not correlated with any of the other 4 associated SNPs. Its pairwise r2 value with each other SNP was below 0.01 (D′max = 0.14 with rs1867283). rs1147198 maps 8 kb 5′ of the transcription start site into the putative promoter region of NTRK2. The other 4 SNPs are located in introns 14 and 16 of full-length NTRK2 (Figure 1A and B). Because 4 of the 5 associated SNPs show low or no LD among each other, they might reflect more than 1 causal variant of the NTRK2 gene conferring risk for SA. Therefore, we performed multilocus regression analysis allowing 2-way and 3-way interactions between the 3 associated SNPs (rs10868235, rs1147198, and rs1867283) with the lowest P values in the combined sample. The model showed a multiplicative effect of the combined SNPs on SA and became highly significant (P = 4.7 × 10−7). Main effects for each single SNP, age, sex, and the interaction between rs10868235 and rs1147198 were significant predictors (P < .05). Figure 2 illustrates the interaction between genotypes of rs10868235I and rs1147198II. The OR of the association with SA was 3.3 (95% confidence interval [CI], 1.7-6.1; P = 1.5 × 10−4) when patients carrying risk genotypes in both interacting loci (CCI or CTI + CCII) were compared with patients without risk genotypes (TTI + AAII or ACII). In analogy, the effect size became even larger when risk genotype carriers in all 3 markers (rs10868235I, rs1147198II, and rs1867283III: CCI or CTI + CCII + AGIII or GGIII) were compared with those without any risk genotype in any marker (TTI + CCII or ACII + AAIII) (OR = 4.5; 95% CI, 2.1-9.8; P = 5.1 × 10−5).
The 3 SNPs (rs10868235, rs1147198, and rs1867283) showing the most robust associations with SA in both German samples were genotyped in African American, nonpsychiatric clinic patients from a cross-sectional study. rs10868235 was excluded from association analysis owing to an insufficient genotyping call rate. Both other SNPs did not show allelic association with a lifetime history of SA in the complete sample (802 subjects without SA vs 119 with SA; P > .05) (Figure 3). However, in the subgroup of probands with a lifetime diagnosis of MDD as ascertained by the Structured Clinical Interview for DSM-IV, the rs1147198 risk allele (C) for SA was significantly more common in the SA-positive MDD group than in the SA-negative MDD group (n = 78; P = .038; OR = 1.80; 95% CI, 1.04-3.11) (Figure 3). Thus, the association is in the same direction as for the German samples. The rs1867283 risk allele (G) for SA was also more frequent in the SA-positive African American MDD group compared with the SA-negative African American MDD group, with trend significance (n = 80; P = .091; OR = 1.88; 95% CI, 0.86-4.13) (Figure 3).
Meta-analyses of these associations with SA among depressed patients in the German and African American samples gave best-model P = 4 × 10−3 for rs1147198 (CC vs AC + AA) and P = 5 × 10−4 for rs1867283 (G vs A).
In this study, we show association of 5 tagging SNPs (Table 3) located within the NTRK2 locus with a lifetime history of SA within depressed patients in 2 independent German samples. This association appears to be carried by several independent risk loci within this gene, and carriers of the risk genotypes of the 3 most significant markers had a 4.5-fold higher risk for SA than depressed patients carrying the nonrisk genotypes. Moreover, we report supportive evidence for these associations in African American patients with MDD. In contrast, no polymorphisms within BDNF were found to be associated with SA. Case-control associations with MDD did not withstand correction for multiple testing in either of the genes. Our results thus support the large body of evidence implicating the neurotrophic system in the pathophysiology of suicide.
The facts that the genetic associations with SA were stronger when comparing depressed patients with SA vs depressed patients without SA than with healthy control subjects and that these SNPs were not associated with MDD suggest that these associations are specific to SA. Moreover, the discovery sample consisted of patients with different mood disorders, whereas the German replication patient sample comprised only patients with recurrent unipolar depression (Table 2). The associations with SA showed very similar effect sizes in these 2 diagnostically different samples and were unaffected by introducing type of mood disorder as a covariate, which further points toward a specificity of these genetic associations to SA.
We could replicate previous studies showing an earlier age at onset of depression and a higher number of previous depressive episodes in those with SA who have mood disorders.55 However, covarying for these variables did not change the strength of the genetic association, suggesting that they did not confound our analysis. Our results are thus in concordance with the hypothesis that psychiatric disease, social risk factors, and other environmental risk factors are additive or interactive factors with genetic susceptibility specific to suicidal behavior (for review, see the articles by Zubenko et al56 and Mann57). While we have sufficient statistical power (74%-99%) to detect each of the 5 described associations with SA among depressed patients in our combined German sample, independent positive replications are needed to further validate these results and answer the question of whether these associations are specific to SA within patients with mood disorder or are also seen in other psychiatric disorders.
Besides testing for association with SA, we also tested association with MDD with our complete marker panel in the discovery sample. We analyzed 353 patients with unipolar depression and 366 nonpsychiatric control subjects. Neither BDNF nor NTRK2 were associated over all tested markers, and no single SNP association withstood correction for multiple testing. While large effects of polymorphisms in these genes are not likely, our study lacked power to detect effects of smaller magnitude (OR < 1.6) as are often reported for complex disease.33,58,59
In contrast to our study, most reported studies have only typed the functional Val66Met polymorphism (rs6265) within the BDNF locus. Interestingly, this polymorphism showed nominal significant association with MDD but not with SA (eTable 1 and eTable 2) in our discovery sample. No association with this variant was seen in our larger German replication sample. Conflicting association results with this SNP and MDD have been reported by several groups. One study35 in Mexican American patients reported a positive association and 4 other studies34,60- 62 failed to show significant single-marker case-control association with this or any other BDNF polymorphism. Nevertheless, Schumacher et al33 showed significant association of a 3-marker haplotype (rs988748-(GT)n-rs6265) that withstood correction for multiple testing in the original German sample and replication in our independent German sample, and Kaufman et al59 reported a significant 3-way interaction between the BDNF Val66Met and serotonin transporter length (5-HTTLPR) polymorphisms and maltreatment history in children with MDD. However, our overall negative results for markers in BDNF and NTRK2 associated with MDD are concordant with results from a large genome-wide case-control association study on MDD (Genetic Association Information Network study),63 where no SNP within BDNF and NTRK2 reached a nominal association P value below .01. Our negative association results with polymorphisms in BDNF and NTRK2 and MDD are thus congruent with most published studies, including recent genome-wide association studies.
Population stratification is a concern in genetic association studies.64 We recorded ethnicity of parents and all 4 grandparents of each proband. Affiliation with an ethnicity other than one of European descent in 1 or more relatives led to the exclusion of a proband from statistical analysis in this study. The proportion of depressed patients with SA having all 4 grandparents of German origin over both studies was 89.6%, which was not different from the proportion in depressed patients without SA (90.2%) (Table 2). We repeated our analysis for the 5 associated SNPs after exclusion of all probands with 1 or more grandparent with a nationality other than German. We could not detect any significant change in the strength of association with SA among depressed patients in this restricted analysis. Furthermore, previously reported data indicate little if any detectable population stratification within individuals of German ethnicity.64 Both German samples explored in this study were recruited in Bavaria, thus exclusively in southern Germany, and therefore it is not likely that our results are substantially biased by population stratification. The fact that similar associations are also seen in an African American sample may point to the fact that these variants are close to functional polymorphisms relevant for several different populations.
In our study, we achieved an almost complete tagging SNP coverage according to HapMap phase I and II data over both genes under investigation. This allowed us to identify several independent common markers associated with SA among depressed patients in NTRK2. Multilocus interaction analysis in the combined German patient sample revealed that the 3 most significant SNPs (Table 3) contribute to the risk for SA in a multiplicative fashion. The model showed P = 4.7 × 10−7 with a significant interaction between 2 of these SNPs (Figure 2). When patients carrying risk genotypes in all 3 markers were compared with those without any risk genotype in any marker, the effect size was 4.5 (95% CI, 2.1-9.8). Four of the 5 SNPs found to be associated with SA were in low LD (Figure 1D). Our results thus suggest that a combination of several independent risk alleles within the NTRK2 locus is associated with SA in depressed patients. While the BDNF locus was tagged down to an MAF of 1%, our study was not powered to detect associations with variants below an MAF of 10%; therefore, we cannot exclude associations with the tested rare variants in this gene.
Besides the herein reported association of polymorphisms with NTRK2 and SA among depressed patients, most recent large meta-analyses support the association of monoaminergic genes with suicidal behavior, including polymorphisms in the tryptophan hydroxylase (TPH), serotonin transporter (SLC6A4),65- 67 and catechol O-methyltransferase (COMT) genes.68 Several single studies report association of corticotrophin-releasing hormone system–related genes (CRH, CRHR1, CRHR2, and CRHBP) and the angiotensin-converting enzyme (ACE) gene, all genes related to the stress hormone system, with suicidal behavior among psychiatric patients or compared with control subjects.69- 72 Animal models suggest that the monoaminergic genes and stress hormone system closely interact with the neurotrophic system. In fact, neurotrophic signaling is decreased following chronic stress, which can be prevented by selective serotonin reuptake–inhibiting antidepressants.73- 75 Individuals with a genetic liability for dysfunctional NTRK2 signaling might more easily pass the threshold for insufficient neurotrophic signaling in the additional presence of a current mood disorder. This could lead to an enhanced impairment of neuronal survival and synaptic plasticity and might be related to the development of suicidal behavior in patients with MDD. In addition, neurobiological factors related to impulsive and aggressive behavior have been shown to be of importance in the development of suicidal behavior.20,76- 78 The combined presence of risk polymorphisms in these pathways and environmental factors might be necessary to trigger an SA. Unfortunately, the sample size of our study does not allow for such analysis.
In summary, this study shows association of several independent common variants of NTRK2 with SA among depressed patients in independent samples. This supports the large body of evidence that dysfunctional neurotrophic signaling might be involved in the pathophysiology of suicidal behavior.
Correspondence: Martin A. Kohli, PhD, Max Planck Institute of Psychiatry, Kraepelinstraße 2-10, D-80804 Munich, Germany (email@example.com).
Submitted for Publication: January 27, 2009; final revision received June 26, 2009; accepted July 29, 2009.
Published Online: February 1, 2010 (doi:10.1001/archgenpsychiatry.2009.201).
Financial Disclosure: Within the last 3 years, Dr Ressler has received research funding support from the National Institute of Mental Health, the National Institute on Drug Abuse, Lundbeck, the Burroughs Wellcome Foundation, and the National Alliance for Research on Schizophrenia and Depression, and he has an unrelated consulting agreement with Tikvah Therapeutics for N-methyl-D-aspartate–based therapeutics. Dr Holsboer has received grant support from Bristol-Myers Squibb, is a founder and shareholder of Affectis Pharmaceuticals, and is a shareholder of Corcept Therapeutics and Neurocrine Biosciences. Dr Binder currently receives grant support from the National Institute of Mental Health and the Doris Duke Charitable Foundation and in the past 2 years has received grant support from Pfizer Pharmaceuticals. Drs Menke, Holsboer, and Binder are inventors of “Means and Methods for Diagnosing Predisposition for Treatment Emergent Suicidal Ideation (TESI),” European patent application 08016477.5; Drs Holsboer and Binder are inventors of “FKBP5: A Novel Target for Antidepressant Therapy,” international patent publication WO 2005/054500; and Drs Holsboer and Binder are inventors of “Polymorphisms in ABCB1 Associated With a Lack of Clinical Response to Medicaments,” international patent application PCT/EP2005/005194.
Funding/Support: The study is supported by a grant from the Exzellenz-Stiftung of the Max Planck Society and by grant FKZ 01GS0481 from the National Genome Research Network, Federal Ministry of Education and Research.
Previous Presentation: This work was presented in part as a poster at the XXIV CINP Congress; July 14, 2008; Munich, Germany.
Additional Contributions: Sabine Damast, Maik Koedel, Susann Sauer, and Alina Tontsch provided technical assistance and Benno Pütz, PhD, Hildegard Pfister, MS, and Macus Ising, PhD, provided database management.