Search tips
Search criteria 


Logo of nihpaAbout Author manuscriptsSubmit a manuscriptHHS Public Access; Author Manuscript; Accepted for publication in peer reviewed journal;
Am J Med Genet B Neuropsychiatr Genet. Author manuscript; available in PMC 2012 June 1.
Published in final edited form as:
PMCID: PMC3108453

Variation in NGFB is Associated with Primary Affective Disorders in Women


Affective disorders (AFDs) are highly comorbid with substance dependence (SD) and both are genetically influenced. However, the specific etiology of the comorbidity is not well understood. We genotyped an array of 1,350 single nucleotide polymorphisms (SNPs) in or near 130 genes in 868 European-Americans (EAs), including 182 individuals with primary AFDs (PAFDs), 214 with SD comorbid with AFD (CAFD), and 472 screened controls. NGFB, which encodes nerve growth factor β and was represented in the array by 15 SNPs, showed the strongest evidence of association, but only among women with PAFDs. Six of the SNPs showed a nominally significant association with PAFDs in women (Ps = 0.0007–0.01); three (rs2856813, rs4332358, and rs10776799) were empirically significant based on 1,000,000 permutations (Ps = 0.008–0.015). Seven haplotypes were significantly associated with PAFDs in women (Ps = 0.0014–0.01), of which six were significant based on empirical permutation analysis (minimal P = 0.0045). Four diplotypes were significantly associated with PAFDs in women (global Ps = 0.001–0.01). The specific diplotype GG-TC, reconstructed from rs2856813 and rs6678788, showed the strongest evidence of association with PAFDs in women (OR = 4.07, P = 4.2E-05). No SNPs or haplotypes were associated with PAFDs in men or with CAFDs in either sex. We conclude that variation in NGFB is a risk factor for PAFDs in women, but not for CAFD.

Keywords: nerve growth factor β, NGFB, affective disorder, substance dependence, association study, sex-specific


Affective disorders (AFDs), characterized by dramatic changes or extremes, or other dysregulation of mood, are highly comorbid with substance dependence (SD). In the U.S. general population, the 12-month prevalence of AFDs and SD were 9.3% and 4.1%, respectively. The observed comorbidity was substantial: the 12-month prevalence of any AFD among survey respondents with diagnosed SD was 29.2%, and that of any SD among survey respondents who were also diagnosed with any AFD was 12.9% [Grant et al., 2004].

Although there appears to be a strong relationship between AFDs and SD, the biological nature of this relationship is not very clear. SD may increase risk for AFDs, AFDs may increase risk for SD, or they may share a common etiology. Compared to primary AFDs (PAFDs, defined here to mean those without comorbid SD), AFDs comorbid with SD (CAFDs) have an earlier age of onset, greater severity of symptoms, worse prognosis, poorer response to treatment, greater functional impairment, and greater risk of suicidal ideation [Davis et al., 2008]. Although on this basis it has been hypothesized both that PAFDs and CAFDs result from entirely different processes [Schuckit, 2006], the issue remains controversial [Markou et al., 1998]. In fact, PAFDs and CAFDs may each have unique causative risk factors as well as factors that overlap with one another. Genetic epidemiologic studies have demonstrated that AFDs and SD are individually heritable; the estimated heritability of bipolar disorders (BPs) ranged from 60–80% [Lichtenstein et al., 2009], and the estimated heritability of major depression (MD) and SD were both ~40% [Fu et al., 2002; Uhl and Grow, 2004] (although substantially higher heritability has been observed for some SD traits). However, it is difficult to identify genes that influence risk for AFDs or the comorbidity of AFDs with SD, in part because both AFDs and SD are complex disorders (i.e., they are influenced by environmental factors and numerous genes that interact with one another).

Our study examined 130 genes for their association to PAFDs and CAFDs by case-control association analysis using a candidate-gene array. Based on the results of that data mining effort, we focused on the gene encoding nerve growth factor beta (NGFB).

NGFB is of interest in this regard for a number of reasons. First, NGFB showed the strongest evidence of association (albeit only in female patients with PAFDs) among the 130 candidate genes in the initial data mining analysis. Second, the functions of the NGFB protein are related to the neurotrophic hypotheses of AFD development [Bersani et al., 2000]. Nerve growth factor (NGF), a member of the neurotrophin family, consists of 3 subunits (alpha, beta, and gamma), apparently with the beta unit (NGFB) being solely responsible for the biological activity of NGF [Barde, 1990]. Binding of NGFB to its receptor TrkA leads to activation of three signaling pathways, PI3K/Akt, Ras/raf/ERK, and PLC/PKC [Sofroniew et al., 2001]. Activation of these pathways is essential for the development, patterning, and maintenance of the nervous system. These pathways are also targets of mood stabilizers [Chen and Manji, 2006; Wada, 2009]. In addition, NGFB is regulated by testosterone [Katoh-Semba et al., 1990; Wright et al., 1987], which is relevant in light of the fact that women have a significantly higher lifetime risk of AFDs than men. The odds ratio of AFDs in females compared with males ranges from 1.6–2.2 [Seedat et al., 2009]. Finally, NGFB is closely associated to stress-related AFDs in animal models [Alleva and Santucci, 2001; Aloe et al., 2002; Marais et al., 2008].



European-American (EA) subjects were recruited at the University of Connecticut Health Center, Yale University School of Medicine, Brown University School of Medicine, and the Medical University of South Carolina. The institutional review board at each of the participating institutions approved the study protocol and consent processes. After subjects received a complete description of the study, they provided written informed consent to participate. Eight hundred sixty-eight unrelated subjects were included, of whom 396 had a lifetime AFD diagnosis (included 182 PAFDs and 214 CAFDs) and 472 were healthy controls. In the PAFD sample set, 167 of 182 subjects were diagnosed with MD, 3 with dysthymia, and 12 with seasonal depression; 11 have MD with anxiety disorders, five have MD with posttraumatic stress disorder, two have MD with obsessive compulsive disorder, and two have MD with social phobia. In the CAFD sample set, the subjects with affective disorders included 46 subjects with BP and 168 with MD. Various substances were abused by individuals in the CAFD group. The characteristics of the sample are shown in Table I.

Characteristics of the European American Samples

Cases met lifetime DSM-III-R or DSM-IV diagnostic criteria for AFDs with or without comorbid SD (including dependence on drugs or alcohol), based on a structured interview, as described below. DSM-IV defines independent mood disorder as occurring in the absence of current substance use, medical illness, or bereavement. PAFD subjects were assessed for lifetime SD and met diagnostic criteria for independent AFDs without SD. All CAFDs were differentiated from PAFDs by virtue of the former having occurred in the context of substance dependence (rather than being accounted for by a medical illness or bereavement). Healthy controls were recruited by advertisement, and screened to exclude individuals with major Axis I psychiatric disorders, including SD, mood disorders, major anxiety disorders, and psychotic disorders. Diagnoses were made using the Structured Clinical Interview for DSM-III-R (SCID-III-R) [Spitzer et al., 1992] or DSM-IV (SCID-IV) [First et al.,1997], the Computerized Diagnostic Interview Schedule for DSM-III-R (CDIS-R) [Ross et al., 1995], or the Semi-structured Assessment for Drug Dependence and Alcoholism [SSADDA][Pierucci-Lagha et al., 2005].

The self-identified ethnicity of case and control subjects was confirmed using a set of ancestry-informative markers (AIMs).

Markers and Genotyping

DNA was extracted from immortalized cell lines, blood, or saliva. A total of 1,350 SNPs within or near 130 candidate genes in 868 samples were genotyped using the Illumina GoldenGate Assay manufacturer’s methods (Illumina Inc., San Diego, CA, USA), as described previously [Hodgkinson et al., 2008]. In addition, 113 AIMs were included in the array.

Markers were excluded if they deviated significantly from Hardy-Weinberg equilibrium (HWE) expectations (P<0.001) in controls or if they had a minor allele frequency (MAF) <0.025, or a SNP call rate <95%.

Statistical Analysis

Data Mining: 130-Candidate-Gene Association Analysis

The chi-square trend model in JMP Genomics 3.1 software (SAS Institute Inc. Cary NC) was used to analyze the association of the 130 candidate genes (1,350 SNPs) to AFDs, PAFDs, and CAFDs. We performed 100,000 permutations to estimate empirical P-values and to address the issue of multiple testing based on the full set of 1,350 SNPs. Statistical significance was defined by an empirical P < 0.05.

NGFB Association Analysis

Single-Marker Association Analysis

Fifteen SNPs, spanning 53,747 bp on chromosome 1 and covering the full length of NGFB (Table II), were analyzed using Haploview 4.1 software in the same case and control samples. The LD relationships of these 15 SNPs are displayed in Figure 1. Genotype-phenotype association analyses were performed by comparing the allele frequency distributions of individual SNPs between cases and controls with χ2 tests. χ2 tests were used to test for deviation from HWE with a cut-off of P = 0.05 (for the single-marker analysis).

FIG. 1
LD Structures of 15 SNPs Spanning the NGFB Gene in the EA Sample Sets. White color represents D′<1 and LOD <2; blue color represents D′=1 and LOD<2; shades of pink/red represent D′<1 and LOD≥2; ...
Characteristics of SNPs in NGFB

Haplotype Association Analysis

Linkage disequilibrium (LD) values between SNP pairs and haplotype association were determined with Haploview version 4.1. Permutation tests with 1,000,000 permutations were used to evaluate significance in this analysis based on 15 single markers and haplotypes.

Diplotype Association Analysis

Diplotypes, also known as multi-locus genotype patterns, were constructed from SNP 9 (rs2856813), SNP10 (rs6678788) SNP 11 (rs4529705), and SNP12 (rs6537860). These SNPs were selected for inclusion in diplotype analyses based on their haplotype block structure (Figure 1) and significant results when taken individually. Fisher’s exact test and likelihood ratio χ2 tests, in the R environment (R Project for Statistical Computing,, were performed to determine individual and global P-values, respectively. Bonferroni correction was used to avoid inflation of type 1 error, accounting for 12 comparisons resulting from the analysis of four genotype patterns for the whole, the sex-specific, and the phenotype-specific samples. This resulted in a statistical significance threshold of p<0.004 (0.05/12).

Population Structure Analysis

We genotyped an additional 348 unrelated African-American subjects with the same Illumina candidate-gene array used above for inclusion in Bayesian cluster analysis to serve as a proxy ancestral population. Based on 113 AIMs, STRUCTURE 2.2 was used to assess population stratification in the case and control subjects by estimating the ancestry proportion for each EA and AA subject. This also established whether the self-reported ethnicity was consistent with the genetically identified ethnicity. Self-identified EA and AA subjects (N=1,216, total) were included in each STRUCTURE run. Parameters were set to 100,000 burn-ins followed by 100,000 iterations; the number of genetic clusters, k, was set from 1 to 5 and the best k was inferred by the rate of change in the log probability of data, LnP(D), between successive k values [Evanno et al., 2005].


Data Mining: Association Analysis for 130 Candidate Genes

Analysis of the 1,350 markers showed that several candidate genes were significantly associated with AFDs (data not shown). NGFB, the most strongly associated gene with PAFDs in women, was of particular interest because it showed differential association in PAFD and CAFD subjects and in women and men. After 100,000 permutations based on 1,350 SNPs, one SNP (rs2856813) in NGFB was significantly associated with AFDs (empirical P = 0.005). However, after stratification by sex, six SNPs in NGFB showed significant association with AFDs in women (empirical Ps = 0.001–0.03), but none in men. After stratification by AFD type, six markers (including rs2856813) were significantly associated with PAFDs (empirical Ps = 0.004–0.04). The observed association for these six SNPs with PAFDs increased in magnitude when only women were included in the analysis (empirical Ps = 0.001–0.02), but no SNP showed significant association with PAFDs in men only. Except for rs2856813, which showed a modest, but significant association with CAFDs in women (empirical P = 0.03), no other markers were associated with CAFDs in either sex (Table III).

The Results of Single-Marker Association Analyses of NGFB after Permutations Based on the Full Set of Candidate Genes

Association Analysis for NGFB

Single-Marker Association Analysis

In this step, we focused on the NGFB gene, which showed the strongest association with PAFDs in women, for the reasons we mentioned in the introduction. We excluded three markers (SNPs 1, 2, and 5) with MAF<0.025 (Table II), leaving 12 SNPs for the association analysis. Considering all of the AFD samples, four markers showed nominally significant association with AFDs (P = 0.004–0.048), and SNP 9 (rs2856813) remained significant after 1,000,000 permutations (empirical P = 0.038). After stratification by sex, seven markers (P’s = 0.0007–0.04), showed stronger significant association with AFDs in women; two of these remained significant after permutations (empirical P’s = 0.008 and 0.018). No markers showed significant association with AFDs in men.

We also examined the associations separately for PAFDs and CAFDs. Six markers displayed stronger nominal association with PAFDs (P’s = 0.004–0.037); these were the same markers previously found in the data mining to be associated with AFDs in women (Table III). After permutation tests, two remained significantly associated with PAFDs (empirical P’s = 0.029 and 0.04). No markers showed significant association with CAFDs.

We also stratified PAFDs and CAFDs by sex. We found that the effects of all six positive markers identified in the PAFD group were attributable to the association in women, which yielded greater statistical significance (P’s = 0.0007–0.01). Three of the SNPs (SNPs 9, 13, and 14) remained significant in the permutation test (empirical P’s = 0.008–0.015). No markers showed association with PAFDs in men. Also, no markers showed significant association with CAFDs in either sex after permutations (Table IV and Figure 2).

FIG. 2
Single-Marker, Haplotype, and Diplotype Association. SM, single-marker; Hap, Haplotype; Dip, diplotype. A, B, and C: Genetic association between NGFB and AFDs, PAFDs, and CAFDs after stratification by AFD type, a, b, and c: Genetic association between ...
Single-Marker Association Analysis

Using JMP Genomics and Haploview software yielded similar results for the single marker association analysis.

Haplotype Association Analysis

The LD structure from the present haplotype analysis was consistent with that calculated using information from the HapMap database for this gene region. We observed 3 blocks, which consisted of SNPs 3, 4, and 6 (block 1), SNPs 9 and 10 (block 2), and SNPs 11 and 12 (block 3) (Figure 1). Taking into account the LD structure and the positive markers in a markerwise analysis, SNPs 9, 11, and 12, identified as Cluster 1, were selected for additional haplotype association analysis.

Haplotype association analysis showed that T-G-G in block 1, A-C and G-T in Block 2, and A-G-G, G-A-A in Cluster 1 were nominally significantly associated with AFDs (Ps = 0.004–0.03), and A-C and A-G-G retained significance (empirical Ps = 0.02 and 0.03, respectively) after 1,000,000 permutations (Table V). The association between haplotypes and AFDs was more extensive and stronger after subgroup and sex stratification, as for single-SNP analysis. Seven haplotypes, including the five listed above, showed greater statistical evidence of association to AFDs in women before multiple testing (Ps = 0.0008–0.04) and four of these remained significant after multiple testing was accounted for (empirical Ps = 0.005–0.0037). However, none of the haplotypes showed significant association with AFDs in men. In the PAFD samples, before multiple test correction, six haplotypes (Ps = 0.0009–0.04) showed significant association with PAFDs, with five remaining significant after permutations (empirical Ps = 0.002–0.04). Seven haplotypes showed significant association with PAFDs in women before (Ps = 0.001–0.01) and six after permutations (empirical Ps = 0.006–0.049). There were no haplotype associations in men. In the CAFD samples, only one haplotype showed nominally significant association with CAFDs in women (P = 0.03), but this was not significant after permutations (Table V and Figure 2).

Haplotype Association Analysis

Since the CAFD group included168 subjects with MD and 46 with BP, we also considered the MD and BP subgroups individually. After multiple test correction, no SNP or haplotype showed association with MD or BP (Supplementary Table I and II).

In this study, the controls were younger (mean age = 30.2) than the CAFD cases (mean age =39.7). As some controls may become affected later in life and thus could be misclassified, we repeated the association analysis excluding control subjects younger than 26 years (resulting mean age for controls, 37.4, N=259). Two single markers (rs2856813 and rs4332358) and two haplotypes (A–C and G-C in block 2) remained significantly association with PAFD in women in the permutation test (empirical Ps = 0.0497–0.0326) (Supplementary Table III), but no markers or haplotypes showed association with PAFD in men. We did not observe any single marker or haplotype associated with CAFD in men women.

Diplotype Analysis

Diplotype association analysis showed strong evidence of association with PAFDs in women. Four genotype patterns (SNPs 9-10, SNPs 11-12, SNPs 9-11-12, and SNPs 9-10-11-12) showed significant association (global Ps = 0.001–0.02) (Table VIa and Figure 2). After Bonferroni adjustment, the association with SNPs 9-10 remained statistically significant (Table VIa). The specific diplotype GG-TC from the SNP 9-10 pattern showed strongly significant association with PAFDs in women (OR = 4.07, P = 0.000042) (Table VIb). There was no evidence of association with PAFDs or CAFDs in men.

Diplotype Association Analysis

Population Structure Analysis

According to the distribution of the rate of change in the log probability of data between successive k values [Evanno et al., 2005], as shown in Figure 3, the number of clusters that best fit the data set was two, corresponding to two ancestries: AAs (additional samples were included to set AA reference ancestry for improved clustering accuracy) and EAs (our current study subjects) (Figure 4). Based on inferred ancestry proportions, the admixed samples were allocated into two populations: “genetic” EAs with European ancestry proportion more than 50%, and “genetic” AAs with African ancestry proportion more than 50%. The concordance between the genetic ancestry and the self-reported ancestry are 99.8% (866/868) in EAs and 98.8% (344/348) in AAs. The average European ancestry proportion was 86% in the 868 genetic EA subjects, and 18.7% in the 348 genetic AA subjects. In the EA subjects, the average European ancestry proportion was 86.3% in cases and 86.4% in controls. There were no statistically significant differences in admixture between case and control groups (t = 0.04, df = 866, P = 0.97). Thus, stratification was ruled out as a source of spurious association in our sample sets.

FIG. 3
Probability Estimates for the Number of Clusters. The ordinate shows the Ln probability corresponding to the number of clusters (k) shown on the abscissa.
FIG. 4
Ancestry Structure of Individuals, x axis = individuals; y axis = ancestry proportions, a = EAs, and b = AAs in the x-axis; European ancestry proportions: red, and African ancestry proportions: green.


We examined the role of variation in NGFB in risk of PAFDs and CAFDs. Our findings showed that NGFB is strongly associated with PAFDs, but not with CAFDs. The association between NGFB and PAFDs was observed only in women.

We excluded population stratification as a source of the observed associations through the analysis of data on 113 SNP AIMs from 868 EAs and 348 AAs. To correct for multiple comparisons in the single-marker and haplotype association analyses, we used the permutation test with 1,000,000 permutations to estimate empirical P-values. For the diplotype association analysis, because we could not identify suitable software to conduct permutation-based analysis, we used the more conservative Bonferroni adjustment. The Bonferroni adjustment is based on the assumption that the tests are independent, so it is conservative because markers in NGFB are in linkage disequilibrium. We observed strong evidence of association between NGFB and PAFDs in women after stringent correction of multiple testing. Although the positive markers within and near NGFB were selected from 1,350 SNPs from the 130-candidate-gene association analysis, we evaluated statistical significance by means of a permutation test based on the full set of 1,350 markers in a data mining step that avoided the accumulation of type I errors. Thus, there is significant evidence of association of NGFB with PAFDs in women, shown in single-marker, haplotype, and diplotype analyses.

The association between NGFB and PAFDs that we observed supports a role of NGFB in the development of PAFDs, consistent with results from animal studies [Aloe et al., 2002; Marais et al., 2008]. NGFB has been shown to play a key role in the development and function of the CNS [Levi-Montalcini, 1987]. Withdrawal of NGFB can lead to impaired neuronal function and neuronal apoptosis [Davies, 2000], which could contribute to the development of AFDs [Drevets et al., 2008; Phillips et al., 2008]. Further, it has been shown in animal models and humans that stressful life events alter endogenous NGF synthesis and/or utilization, which in humans lead to the development of AFDs [Alleva and Santucci, 2001; Aloe et al., 2002]. When the first parachute jumps of soldiers were used as a model of human stress, NGF levels in peripheral lymphocytes increased 84% in pre-jump and 107% in post-jump sampling [Alleva and Santucci, 2001]. Endogenous release of NGF during exposure to stressors has been shown to remodel damaged tissues, as a compensatory mechanism [Aloe et al., 2002; Marais et al., 2008]. In animal models, treatment with lithium and electroconvulsive stimuli increased NGF concentrations in the frontal cortex, hippocampus, amygdala and limbic forebrain [Angelucci et al., 2003; Hellweg et al., 2002].

Our results also indicated that the role of NGFB in PAFDs might be sex specific. A striking sex difference in the prevalence of AFDs can be seen in different countries and cultures, suggesting that there is a biological basis for it [Altemus, 2006]. However, efforts to identify the biological origin of these sex differences have remained largely unsuccessful. We found that variants in NGFB showed significant association in women with PAFDs, but not in men. Our finding is similar to the recent observation that variation in NGFB confers susceptibility to anxiety-related personality traits in a sex-dependent way [Lang et al., 2008]. Findings from animal studies show that NGF can be regulated by sex hormones and that females produce less NGF than males [Levi-Montalcini, 1987; Levi-Montalcini and Angeletti, 1968; Wright et al., 1987], which may help to explain the sex-dependent effect observed in our study. Castration in males drastically reduces NGF concentrations in brain, spinal cord and salivary gland, whereas exogenous administration of testosterone to either young or adult females or castrated males increases the concentration of NGF in these tissues [Katoh-Semba et al., 1990; Levi-Montalcini and Angeletti, 1968]. In rhesus macaques, the plasma concentration of NGF increased in response to stress only in males [Cirulli et al., 2009]. Animal studies have shown that stress early in pregnancy may produce a significant sex-dependent effect on placental gene expression, which could alter fetal transport of growth factors and nutrients. These alterations could then exert downstream effects resulting in sex differences in stress sensitivity and psychiatric disease risk (Goel et al., 2009). Males had higher concentrations of NGF in peripheral and brain tissues than females, possibly representing a better adaptation to stressful life events. The sex-dependent effects of NGFB in PAFDs might help to explain why AFDs are more prevalent in women than in men.

Our findings also suggest that PAFDs and CAFDs may have different underlying mechanisms, since we found no evidence for association between NGFB and CAFDs. PAFDs might result in part from inadequate amounts of endogenous NGF associated with NGFB variation, hampering neuronal development and remodeling of damaged neurons following stressful events. On the other hand, risk for CAFDs might be regulated by the effects of neurotoxic substances such as drugs of abuse, which might secondarily lead to inadequate NGFB expression, an environmental effect that does not depend on the presence of variation in NGFB per se. Consistent with this interpretation, recent studies revealed that chronic heroin, cocaine, cannabis, and alcohol abuse reduced serum concentrations of NGFB [Angelucci et al., 2007; Angelucci et al., 2008]. In view of the neurotrophic hypothesis of AFDs and the findings from these studies, reduced concentration of NGFB could increase the risk of developing AFDs.

Limitations of the current study include the fact that we did not differentiate independent from substance-induced AFDs in CAFD samples because of the limited sample and the difficulties in identifying the exact chronological order of onset of AFDs and SD in individuals with both disorders. Although the CAFD group may include some individuals with primary AFDs comorbid with SD, it is unlikely that the inclusion of such individuals in this subgroup confounded the analysis, because such misclassification would be expected to result in an association between NGFB and CAFD, which we did not observe. On the other hand, individuals with both independent and substance-induced major depressive episodes (MDEs) had more severe psychopathology than those with either independent or substance-induced MDE only [Niciu et al., 2009]. This provides support for our interpretation that our findings reflect different pathophysiologic effects in PAFDs and CAFDs, with individuals suffering from both types of AFDs having “a double dose” of pathogenic risk and a concomitantly more severe clinical picture. In addition, our findings derived from EAs, and thus they are applicable to this population only It is necessary that these findings be confirmed in other populations,

Despite these limitations, our study provides a number of possible genetic and pathophysiological insights into AFDs. In addition to showing that PAFDs and CAFDs, which represent common disorders, might have different genetic backgrounds and pathophysiologic processes, we found that the association of NGFB with PAFDs was sex dependent. This could help to explain the higher prevalence of AFDs in women. The findings may also help to explain why PAFDs respond well to medications, which may serve to alleviate genetically influenced symptoms in PAFDs. In contrast, CAFDs appear to respond effectively to abstinence from substances of abuse [Brown and Schuckit, 1988; Nakamura et al., 1983]. If these findings are confirmed in EAs and extended to other populations, it may be possible to use variation in NGFB to determine the best treatment for individuals with affective disorders.

Supplementary Material

Supp Table S1-S3

Supplementary Table I Single-Marker Association Analysis of MD and BP subgroups in CAFD

Supplementary Table II Haplotype Association Analysis of MD and BP subgroups in CAFD

Supplementary Table III Single-Marker and Haplotype Association Analysis in PAFD after excluding those controls with age <26.


This work was supported by the VA CT REAP Center, and NIH grants AA11330, DA12849, DA018432, DA12690, RR06192 (University of Connecticut General Clinical Research Center), and K05 AA017435. We thank the families and individuals who volunteered to participate in this research. Ann Marie Lacobelle, Michelle Streckenbach, Gregory Dalton-Kay, and Christa Robinson provided excellent technical assistance. We thank Shizhong Han, Ph.D. for his assistance with graphics. Colin A. Hodgkinson, Ph.D. at the National Institute on Alcohol Abuse and Alcoholism helped with genotyping by the Illumina GoldenGate Assay methodology.


Some of these data were presented at the 58th Annual Meeting of the American Society of Human Genetics, Philadelphia PA, November 11–15, 2008


  • Alleva E, Santucci D. Psychosocial vs. “physical” stress situations in rodents and humans: role of neurotrophins. Physiol Behav. 2001;73(3):313–20. [PubMed]
  • Aloe L, Alleva E, Fiore M. Stress and nerve growth factor: findings in animal models and humans. Pharmacol Biochem Behav. 2002;73(1):159–66. [PubMed]
  • Altemus M. Sex differences in depression and anxiety disorders: potential biological determinants. Horm Behav. 2006;50(4):534–8. [PubMed]
  • Angelucci F, Aloe L, Jimenez-Vasquez P, Mathe AA. Electroconvulsive stimuli alter nerve growth factor but not brain-derived neurotrophic factor concentrations in brains of a rat model of depression. Neuropeptides. 2003;37(1):51–6. [PubMed]
  • Angelucci F, Ricci V, Pomponi M, Conte G, Mathe AA, Attilio Tonali P, Bria P. Chronic heroin and cocaine abuse is associated with decreased serum concentrations of the nerve growth factor and brain-derived neurotrophic factor. J Psychopharmacol. 2007;21(8):820–5. [PubMed]
  • Angelucci F, Ricci V, Spalletta G, Pomponi M, Tonioni F, Caltagirone C, Bria P. Reduced serum concentrations of nerve growth factor, but not brain-derived neurotrophic factor, in chronic cannabis abusers. Eur Neuropsychopharmacol. 2008;18(12):882–7. [PubMed]
  • Barde YA. The nerve growth factor family. Prog Growth Factor Res. 1990;2(4):237–48. [PubMed]
  • Bersani G, Iannitelli A, Fiore M, Angelucci F, Aloe L. Data and hypotheses on the role of nerve growth factor and other neurotrophins in psychiatric disorders. Med Hypotheses. 2000;55(3):199–207. [PubMed]
  • Brown SA, Schuckit MA. Changes in depression among abstinent alcoholics. J Stud Alcohol. 1988;49(5):412–7. [PubMed]
  • Chen G, Manji HK. The extracellular signal-regulated kinase pathway: an emerging promising target for mood stabilizers. Curr Opin Psychiatry. 2006;19(3):313–23. [PubMed]
  • Cirulli F, Francia N, Branchi I, Antonucci MT, Aloe L, Suomi SJ, Alleva E. Changes in plasma levels of BDNF and NGF reveal a gender-selective vulnerability to early adversity in rhesus macaques. Psychoneuroendocrinology. 2009;34(2):172–80. [PMC free article] [PubMed]
  • Davies AM. Neurotrophins: more to NGF than just survival. Curr Biol. 2000;10(10):R374–6. [PubMed]
  • Davis L, Uezato A, Newell JM, Frazier E. Major depression and comorbid substance use disorders. Curr Opin Psychiatry. 2008;21(1):14–8. [PubMed]
  • Drevets WC, Savitz J, Trimble M. The subgenual anterior cingulate cortex in mood disorders. CNS Spectr. 2008;13(8):663–81. [PMC free article] [PubMed]
  • Evanno G, Regnaut S, Goudet J. Detecting the number of clusters of individuals using the software STRUCTURE: a simulation study. Mol Ecol. 2005;14(8):2611–20. [PubMed]
  • First M, Spitzer RL, Gibbon M, Williams JB. Structured Clinical Interview for DSM-IV AXIS I Disorders (Clinician Version) SCID-I Administration Booklet. American Psychiatric Publishing. Inc; 1997.
  • Fu Q, Heath AC, Bucholz KK, Nelson E, Goldberg J, Lyons MJ, True WR, Jacob T, Tsuang MT, Eisen SA. Shared genetic risk of major depression, alcohol dependence, and marijuana dependence: contribution of antisocial personality disorder in men. Arch Gen Psychiatry. 2002;59(12):1125–32. [PubMed]
  • Grant BF, Stinson FS, Dawson DA, Chou SP, Dufour MC, Compton W, Pickering RP, Kaplan K. Prevalence and co-occurrence of substance use disorders and independent mood and anxiety disorders: results from the National Epidemiologic Survey on Alcohol and Related Conditions. Arch Gen Psychiatry. 2004;61(8):807–16. [PubMed]
  • Goel N, Bale TL. Examining the Intersection of Sex and Stress in Modeling Neuropsychiatric Disorders. Journal of Neuroendocrinology. 2009;21:415–420. [PMC free article] [PubMed]
  • Hellweg R, Lang UE, Nagel M, Baumgartner A. Subchronic treatment with lithium increases nerve growth factor content in distinct brain regions of adult rats. Mol Psychiatry. 2002;7(6):604–8. [PubMed]
  • Hodgkinson CA, Yuan Q, Xu K, Shen PH, Heinz E, Lobos EA, Binder EB, Cubells J, Ehlers CL, Gelernter J, et al. Addictions biology: haplotype-based analysis for 130 candidate genes on a single array. Alcohol Alcohol. 2008;43(5):505–15. [PMC free article] [PubMed]
  • Katoh-Semba R, Semba R, Kashiwamata S, Kato K. Influences of neonatal and adult exposures to testosterone on the levels of the beta-subunit of nerve growth factor in the neural tissues of mice. Brain Res. 1990;522(1):112–7. [PubMed]
  • Lang UE, Hellweg R, Bajbouj M, Gaus V, Sander T, Gallinat J. Gender-dependent association of a functional NGF polymorphism with anxiety-related personality traits. Pharmacopsychiatry. 2008;41(5):196–9. [PubMed]
  • Levi-Montalcini R. The nerve growth factor 35 years later. Science. 1987;237(4819):1154–62. [PubMed]
  • Levi-Montalcini R, Angeletti PU. Nerve growth factor. Physiol Rev. 1968;48(3):534–69. [PubMed]
  • Lichtenstein P, Yip BH, Bjork C, Pawitan Y, Cannon TD, Sullivan PF, Hultman CM. Common genetic determinants of schizophrenia and bipolar disorder in Swedish families: a population-based study. Lancet. 2009;373(9659):234–9. [PubMed]
  • Marais L, van Rensburg SJ, van Zyl JM, Stein DJ, Daniels WM. Maternal separation of rat pups increases the risk of developing depressive-like behavior after subsequent chronic stress by altering corticosterone and neurotrophin levels in the hippocampus. Neurosci Res. 2008;61(1):106–12. [PubMed]
  • Markou A, Kosten TR, Koob GF. Neurobiological similarities in depression and drug dependence: a self-medication hypothesis. Neuropsychopharmacology. 1998;18(3):135–74. [PubMed]
  • Nakamura MM, Overall JE, Hollister LE, Radcliffe E. Factors affecting outcome of depressive symptoms in alcoholics. Alcohol Clin Exp Res. 1983;7(2):188–93. [PubMed]
  • Niciu MJ, Chan G, Gelernter J, Arias AJ, Douglas K, Weiss R, Anton RF, Farrer L, Cubells JF, Kranzler HR. Subtypes of major depression in substance dependence. Addiction. 2009;104:1700–1709. [PMC free article] [PubMed]
  • Phillips ML, Ladouceur CD, Drevets WC. A neural model of voluntary and automatic emotion regulation: implications for understanding the pathophysiology and neurodevelopment of bipolar disorder. Mol Psychiatry. 2008;13(9):829, 833–57. [PMC free article] [PubMed]
  • Pierucci-Lagha A, Gelernter J, Feinn R, Cubells JF, Pearson D, Pollastri A, Farrer L, Kranzler HR. Diagnostic reliability of the Semi-structured Assessment for Drug Dependence and Alcoholism (SSADDA) Drug Alcohol Depend. 2005;80(3):303–12. [PubMed]
  • Ross HE, Swinson R, Doumani S, Larkin EJ. Diagnosing comorbidity in substance abusers: a comparison of the test-retest reliability of two interviews. Am J Drug Alcohol Abuse. 1995;21(2):167–85. [PubMed]
  • Schuckit comorbidity between substance use disorders and psychiatric conditions. research report. 2006;101:76–88. [PubMed]
  • Seedat S, Scott KM, Angermeyer MC, Berglund P, Bromet EJ, Brugha TS, Demyttenaere K, de Girolamo G, Haro JM, Jin R, et al. Cross-national associations between gender and mental disorders in the World Health Organization World Mental Health Surveys. Arch Gen Psychiatry. 2009;66(7):785–95. [PMC free article] [PubMed]
  • Sofroniew MV, Howe CL, Mobley WC. Nerve growth factor signaling, neuroprotection, and neural repair. Annu Rev Neurosci. 2001;24:1217–81. [PubMed]
  • Spitzer RL, Williams JB, Gibbon M, First MB. The Structured Clinical Interview for DSM-III-R (SCID). I: History, rationale, and description. Arch Gen Psychiatry. 1992;49(8):624–9. [PubMed]
  • Uhl GR, Grow RW. The burden of complex genetics in brain disorders. Arch Gen Psychiatry. 2004;61(3):223–9. [PubMed]
  • Wada A. Lithium and neuropsychiatric therapeutics: neuroplasticity via glycogen synthase kinase-3beta, beta-catenin, and neurotrophin cascades. J Pharmacol Sci. 2009;110(1):14–28. [PubMed]
  • Wright LL, Beck C, Perez-Polo JR. Sex differences in nerve growth factor levels in superior cervical ganglia and pineals. Int J Dev Neurosci. 1987;5(5–6):383–90. [PubMed]