Main and interaction effects of CRHR1 single-nucleotide polymorphisms (SNPs) with child abuse on adult depression. A, CRHR1 linkage disequilibrium (LD) map demonstrating the physical location (arrows) and LD pattern of the selected CRHR1 SNPs. For initial genotyping, we selected 15 SNPs covering the gene from the promoter region to the 3′ end of the gene, with an average intermarker distance of 4.1 kilobases (kb) over a total region of 57 kb. The x-axis shows the position on chromosome 17. Ten of these SNPs had a minor allele frequency greater than 5.0% and were thus included in the analysis. The lower part shows an LD plot generated with Haploview using r2 as the measure of LD, which ranges from 1 (or complete LD, indicated by black squares) to 0 (or absence of LD, indicated by white squares). In our population, we observed moderate LD in this gene, especially in the part containing the coding regions with 2 separate blocks of LD according to the confidence interval method. B, Main genetic effects. The P values (log scale) for the main genetic effect (for a model containing SNP genotype, Childhood Trauma Questionnaire scores, age, and sex as independent variables and Beck Depression Inventory scores as outcome) for all of the genotyped subjects (n = 422), subjects with no to mild child abuse (n = 262), and subjects with moderate to severe child abuse (n = 160). Concordant with a gene × environment interaction effect, the main genetic effect is carried by individuals with moderate to severe child abuse and is stronger than in the whole sample, whereas no SNPs significant in the whole sample remain significant in the group with no to mild child abuse. C, Interaction effects. The y-axis shows the P value (log scale) for the interaction effect of the 10 SNPs that were present with a minor allele frequency greater than 5.0%. Horizontal lines indicate P = .05 (nominally significant) and P = .0094 (significance level after correction for multiple testing as described in the text).
Child abuse interacts with the CRHR1 genotype to enhance risk for depression in adults. A, Mean Beck Depression Inventory (BDI) scores (n = 476) are predicted by continuous scores on the Childhood Trauma Questionnaire (CTQ) physical, sexual, and emotional abuse scales divided into quartiles (quartile 1: CTQ score = 15.0-17.5 [n = 117]; quartile 2: CTQ score = 17.5-22.0 [n = 112]; quartile 3: CTQ score = 22.0-31.0 [n = 123]; quartile 4: CTQ score > 31 [n = 124]). The error bars indicate SEM. The overall analysis of variance comparing the 4 groups on the BDI was significant (F3,472 = 37.12; P < .001). Post hoc analysis indicated significant differences (P < .001) between all pairs of groups except for between groups 2 and 3. Effect of rs7209436 (B), rs110402 (C), and rs242924 (D) genotypes and child abuse on adult depression. Error bars indicate SEM. The sample sizes for no to mild child abuse and moderate to severe child abuse are as follows: rs7209436: CC, 110 and 84, respectively; CT, 115 and 61, respectively; TT, 23 and 9, respectively (B); rs110402: GG, 108 and 84, respectively; AG, 114 and 61, respectively; AA, 23 and 9, respectively (C); and rs242924: GG, 112 and 87, respectively; GT, 111 and 60, respectively; TT, 24 and 10, respectively (D). P values of significance (based on conservative linear permutation analyses) for the main effect of genotype in the moderate to severe abuse group are shown.
Effect of rs242940 (A), rs173365 (B), rs4792887 (C), and rs242948 (D) genotypes and child abuse on adult depression. Error bars indicate SEM. The sample sizes for no to mild child abuse and moderate to severe child abuse are as follows: rs242940: AA, 98 and 74, respectively; AG, 113 and 69, respectively; GG, 31 and 13 respectively (A); rs173365: GG, 79 and 56, respectively; AG, 114 and 78, respectively; AA, 51 and 21, respectively (B); rs4792887: CC, 107 and 58, respectively; CT, 105 and 70, respectively; TT, 29 and 26, respectively (C); and rs242948: GG, 84 and 57, respectively; GT, 107 and 79, respectively; TT, 64 and 22, respectively (D). P values of significance for the main effect of genotype in the moderate to severe abuse group are shown.
Effect of estimated CRHR1 haplotypes and child abuse on adult depressive symptoms. The frequencies of the estimated individual haplotypes within the first linkage disequilibrium block (A) and of the 3 most significant single-nucleotide polymorphisms (SNPs) (B). C, The effects of the estimated copy numbers of specific haplotypes on adult depression in the presence or absence of moderate to severe child abuse. Error bars indicate SEM. Only haplotypes with an estimation likelihood greater than 95.0% were included in the analysis. The sample sizes for no to mild child abuse and moderate to severe child abuse were as follows: TCA haplotype: 0 copies, 107 and 82, respectively; 1 copy, 109 and 60, respectively; 2 copies, 21 and 9, respectively; TAT haplotype: 0 copies, 118 and 91, respectively; 1 copy, 106 and 54, respectively; 2 copies, 19 and 9, respectively. P values of significance for the interaction effect of haplotype in the moderate to severe abuse group are shown. D, A graphical representation of the number of copies of the TAT haplotype in a smaller but ethnically and socioeconomically separate sample. In this separate sample, the best SNP haplotype frequencies were 41.6% for TAT and 56.0% for CGG. The frequencies of TAT haplotypes are shown within child abuse groups for those with no history of major depressive disorder (MDD) and those with a lifetime history of MDD. The protective alleles appear to be present at a higher rate in those with moderate to severe child abuse who have no lifetime MDD compared with those with a lifetime history of MDD.
Customize your JAMA Network experience by selecting one or more topics from the list below.
Bradley RG, Binder EB, Epstein MP, et al. Influence of Child Abuse on Adult Depression: Moderation by the Corticotropin-Releasing Hormone Receptor Gene. Arch Gen Psychiatry. 2008;65(2):190–200. doi:10.1001/archgenpsychiatry.2007.26
Genetic inheritance and developmental life stress both contribute to major depressive disorder in adults. Child abuse and trauma alter the endogenous stress response, principally corticotropin-releasing hormone and its downstream effectors, suggesting that a gene × environment interaction at this locus may be important in depression.
To examine whether the effects of child abuse on adult depressive symptoms are moderated by genetic polymorphisms within the corticotropin-releasing hormone type 1 receptor (CRHR1) gene.
Association study examining gene × environment interactions between genetic polymorphisms at the CRHR1 locus and measures of child abuse on adult depressive symptoms.
General medical clinics of a large, public, urban hospital and Emory University, Atlanta, Georgia.
The primary participant population was 97.4% African American, of low socioeconomic status, and with high rates of lifetime trauma (n = 422). A supportive independent sample (n = 199) was distinct both ethnically (87.7% Caucasian) and socioeconomically (less impoverished).
Main Outcome Measures
Beck Depression Inventory scores and history of major depressive disorder by the Structured Clinical Interview for DSM-IV Axis I Disorders.
Fifteen single-nucleotide polymorphisms spanning 57 kilobases of the CRHR1 gene were examined. We found significant gene × environment interactions with multiple individual single-nucleotide polymorphisms (eg, rs110402, P = .008) as well as with a common haplotype spanning intron 1 (P < .001). Specific CRHR1 polymorphisms appeared to moderate the effect of child abuse on the risk for adult depressive symptoms. These protective effects were supported with similar findings in a second independent sample (n = 199).
These data support the corticotropin-releasing hormone hypothesis of depression and suggest that a gene × environment interaction is important for the expression of depressive symptoms in adults with CRHR1 risk or protective alleles who have a history of child abuse.
Although numerous studies suggest that both genetic and environmental influences contribute substantively to vulnerability to major depressive disorder (MDD), the search for direct genetic influences contributing to the risk for depression has been inconclusive so far. Recent studies1,2 examining genetic variants predisposing to MDD suggest that they may only do so via interaction with environmental variables. Several reports have now replicated the finding that the short allele of a functional repeat polymorphism located in the promoter region of the serotonin transporter locus SLC6A43 increases vulnerability to depression only in the presence of significant adverse life events.1,4-8 Although the serotonin transporter is an attractive genetic candidate for moderating the adverse effects of life events, epidemiological studies clearly indicate that the vulnerability to depression is influenced by several genes.9
Genes regulating the effects of stress via the hypothalamic-pituitary-adrenal (HPA) axis have been implicated in the physiological and pathological regulation of stress reactivity.10,11 The physiological stress response is primarily mediated by the release of hypothalamic corticotropin-releasing factor, also known as corticotropin-releasing hormone (CRH), from nerve terminals in the median eminence arising in the paraventricular nucleus, which stimulates adrenocorticotropin release from the anterior pituitary that in turn induces the release of cortisol from the adrenal cortex. Overactivity of the HPA axis is a hallmark of MDD and is at least partly due to a hyperactivity of CRH neurons.12,13
In addition to its effects on the HPA axis, CRH activity at the CRH type 1 receptor (CRHR1) in extrahypothalamic regions is also thought to produce symptoms of anxiety and depression. Preclinical animal models of depression14,15 and clinical studies of patients with depression16,17 have found increased CRH concentrations in the cerebrospinal fluid as well as altered CRHR1 messenger RNA expression and CRH binding in limbic brain regions (eg, hypothalamus, amygdala, and other areas mediating emotion response), lending support to this theory. Furthermore, a nonblinded study demonstrated clinical efficacy of a CRHR1 antagonist in patients with depression.18 Finally, single-nucleotide polymorphisms (SNPs) in the CRHR1 gene (Entrez Gene ID 1394 or RefSeq ID NM_004382) have been associated with MDD in Han Chinese individuals,19 and response to antidepressant medications appears to vary with CRHR1 haplotypes in Mexican American persons.20
Studies suggest that HPA axis hyperactivity related to MDD may be a function of early life stress (ELS).21,22 Given the relationships between ELS, MDD, and the dysregulation of the CRH and CRHR1 system, we hypothesized that genetic polymorphisms altering the function of CRHR1 may moderate the effects of child abuse on adult depression.
The data from this study were collected as part of a larger study investigating the roles of genetic and environmental factors in predicting response to stressful life events in a predominantly African American urban population of low socioeconomic status. The sample for the replication study is discussed later. Research participants were approached, either while waiting for their medical appointments or while waiting with others who were scheduled for medical appointments, in the waiting rooms of the primary care clinic or obstetrical-gynecological clinic of a large, urban, public hospital. Subjects who indicated willingness to participate provided written informed consent, participated in a verbal interview, and provided a salivary sample for DNA extraction (described later). The data presented in this article are from the first 560 subjects screened (we have incomplete genotype or phenotype data for some subjects, so the total number of subjects in the analysis will be fewer than 560, with exact sample sizes stated in the “Results” section). The final genetic analysis was performed on a subset of 422 of these subjects. All of the procedures in this study were approved by the institutional review boards of Emory University School of Medicine and Grady Memorial Hospital, Atlanta, Georgia. Demographic data on sex, self-identified race/ethnicity, education, employment, disability status, and monthly household income are presented in Table 1. The mean (SD) age in the sample was 38.4 (13.3) years (range, 18-81 years).
Participants verbally completed a battery of self-report measures, including the Beck Depression Inventory (BDI) and the brief Childhood Trauma Questionnaire (CTQ), which took 45 to 75 minutes to complete (dependent in large part on the extent of the participant's trauma history and symptoms). We read instruments to participants to guard against relatively high rates of impaired literacy. Each person was paid $15.00 for participation in this phase of the study.
Depressed mood was assessed with the 21-item BDI,23 a well-validated, commonly used, continuous measure of the level of depressive symptoms. In this sample, the BDI had a standardized α coefficient of .99 (mean [SD] BDI score, 14.43 [13.11]). In addition, we examined validity of the BDI using data collected in a later phase of the study for which the Structured Clinical Interview for DSM-IV Axis I Disorders (SCID-I),24 a validated interview assessment of DSM-IV mood disorders, was administered. In this subsample of participants (n = 202), 19.3% met the criteria for current MDD and 40.1% met criteria for lifetime MDD. In addition, 1.1% met criteria for current bipolar I disorder and 2.0% met lifetime criteria for bipolar I disorder; 2.2% met criteria for current bipolar II disorder and 3.0% met lifetime criteria for bipolar II disorder. Based on MDD diagnosis with the SCID-I, the BDI total scores for those with current MDD (mean [SD]
BDI total score, 28.26 [15.01]) were significantly different from the BDI total scores for those without a current MDD diagnosis (mean [SD] BDI total score, 13.48 [12.73]) (F1,173 = 32.9; P < .001). Additionally, because these analyses were performed in subjects identified in a general medical clinic, it is possible that the BDI score is a proxy for medical-related depressive symptoms separate from MDD. To examine this, we used history of prescribed medications as a proxy variable for the medical illness category. Although we have only limited diagnostic information, we have medical pharmaceutical data for approximately 120 subjects from the full genotype and phenotype data set. Similar to the linear regression analyses with permutations for SNP analyses (described later), we used linear regression with permutation procedures (using 10 000 permutations with BDI as the dependent variable and each of the 13 medical variables as predictors) to establish whether the proxy variables for different types of medical illness predict the BDI score. These permutation procedures randomly assigned the sample BDI scores to subjects (sampled without replacement) while holding each subject's medical prescription variables fixed. We conducted these analyses using appropriate components of SAS version 9.1 statistical software (SAS Institute, Inc, Cary, North Carolina). There were no significant effects found for any of the 13 available classes of medications examined (all P > .10). This demonstrates that the moderation effects of child abuse and CRHR1 genotype were unlikely affecting medical condition as a possible contributor to depressive symptoms.
The CTQ25 is a self-report inventory assessing 3 types of childhood abuse: sexual, physical, and emotional. Studies have established the internal consistency, stability over time, and criterion validity of both the original 70-item CTQ and the current brief version.26 The CTQ yields a total score and subscale scores for each of the types of child abuse. Our CTQ data demonstrated good internal consistency reliability (α = .99 for physical abuse; α = .94 for sexual abuse; α = .93 for emotional abuse; and α = .98 for the total of these 3 scales). Bernstein and Fink25 established scores for none, mild, moderate, and severe for each type of abuse. The data from the CTQ were used to classify participants into 2 categories for each type of abuse (physical, sexual, and emotional): (1) those with CTQ scale scores in the none to mild range, and (2) those with CTQ scores in the moderate to severe range. We then created a composite variable across all of the 3 types of abuse. Using this composite, we divided participants into 2 groups with respect to the numbers of types of abuse that fell into the moderate to severe range: (1) those with no type of abuse in the moderate to severe range, and (2) those with at least 1 type of abuse in the moderate to severe range.
A sample of 204 women recruited at the Women's Mental Health Center at Emory University had been assessed for lifetime psychiatric diagnoses using the SCID-I24 and for the presence of child abuse using the CTQ. The severity of child abuse was then scored as none to mild vs moderate to severe. In this sample, the mean (SD) age was 31.9 (5.0) years, and 87.7% were Caucasian or white, 7.0% were African American or black, 3.7% were Native American, and 1.6% were Asian. As compared with the previous sample, subjects in this sample were primarily middle-income Americans. The mean (SD) number of years of education in this sample was 16.0 (2.0) years. Among the subjects, 23.6% reported moderate to severe child abuse. Most subjects (96.4%) fulfilled criteria for at least 1 Axis I DSM-IV diagnosis, with 62.6% of the patients presenting with a principal lifetime diagnosis of MDD. The other reported lifetime diagnoses included bipolar I disorder, bipolar II disorder, panic disorder, generalized anxiety disorder, anxiety disorder not otherwise specified, obsessive-compulsive disorder, specific phobias, dysthymia, depressive disorder not otherwise specified, brief psychotic disorder, schizophrenia, mood disorder due to the general medical condition, cannabis abuse, opioid abuse, posttraumatic stress disorder, adjustment disorder, and binge eating. Women with these latter diagnoses did not fulfill criteria for MDD at any time and represented the group with no MDD.
We extracted DNA from saliva collected into Scope mouthwash (Procter and Gamble, Cincinnati, Ohio) (n = 46) or into Oragene saliva kits (DNA Genotek Inc, Ottawa, Ontario, Canada) using the Qiagen M48 system (Qiagen, Hilden, Germany). High-quality DNA was available for 502 individuals.
Fifteen SNPs validated in African populations were selected from public databases (dbSNP, http://www.ncbi.nlm.nih.gov/projects/SNP/; HapMap, http://www.hapmap.org) to span the CRHR1 gene (NM_004382) from the promoter region to the 3′ end of the gene, with an average intermarker distance of 4.1 kilobases (kb) over a total region of 57 kb. Ten of these SNPs had a minor allele frequency greater than 5.0% and were thus included in the analysis (full SNP information is available from us, and identification numbers of 10 selected SNPs are listed in Figure 1A). In our population, we observed moderate linkage disequilibrium (LD) in this gene, especially in the part containing the coding regions with 2 separate blocks of LD (Figure 1A). All of the SNPs were genotyped using a TaqMan allelic discrimination assay27 developed for use on the 7900HT instrument (Applied Biosystems, Foster City, California) in 502 individuals using predesigned and validated TaqMan assay reagent kits containing 1 pair of polymerase chain reaction primers and 1 pair of fluorescently labeled probes (Applied Biosystems). Assay identification numbers for all of the 15 assays are available on request. Polymerase chain reactions were performed in 5-μL reaction volumes in 384-well plates and contained 5 ng of DNA. The standard protocol provided with the kit was followed. Thermal cycler conditions were 95°C for 10 minutes and then 40 cycles of 95°C for 15 seconds and 60°C for 1 minute. The SDS 2.2 software (Applied Biosystems) was used for allelic discrimination. For quality control, 9.0% of the samples were genotyped as duplicates across and within a 384-well plate. Only 2 discordances were recorded with the 6255 unique genotypes (0.03%) and excluded from the analysis. All of the SNP identification numbers, their locations, Hardy-Weinberg equilibrium, and minor allele frequencies are available on request.
Genotyping of ancestry-informative markers (AIMs) and data analysis for estimates of ancestry were performed by DNA Print Genomics, Sarasota, Florida (http://www.dnaprint.com). One hundred thirty-four AIMs were genotyped in 280 individuals from our sample using single-base primer extension and a Beckman SNPstream instrument (Beckman-Coulter, Fullerton, California). Maximum likelihood estimates of individual biogeographical ancestry admixture were determined using methods described previously by others28 and 134 AIMs29-31 chosen from the genome based on their information for a continental, 4-population model. The choice of a 4-population model was based on previously published hypothesis-free cluster analyses of worldwide populations,31,32 and this model lent itself to use of the terms European (genetic ancestry shared among European, Middle Eastern, and, to a lower extent, South Asian populations), sub-Saharan African, East Asian, and indigenous American (genetic ancestry shared among American Indian, Latin and South American Indian, and certain Central Asian populations). The names of these parental populations were chosen to describe extant elements of genetic structure and are arbitrary in that they are reflective of modern population distributions—not necessarily the distributions of the original parental populations 15 000 to 50 000 years ago. Individuals with more than 35 failed AIM genotypes were excluded from the analysis (7.1%) because their biogeographical ancestry admixture could not be estimated with sufficient certainty.
We used Haploview33 to determine the LD structure of the SNPs within the CRHR1 gene and test for Hardy-Weinberg equilibrium (Figure 1A). We also compared the LD structure of a subgroup of CRHR1 SNPs across our study population and 2 populations genotyped in the HapMap project (Yorubans [a sub-Saharan African population] and Centre d’Étude du Polymorphisme Humain [Caucasian population]). We compared a tag SNP selection for the CRHR1 locus based on these HapMap populations covering the same area on chromosome 17 using an r2 cutoff of 0.75 and a minor allele frequency greater than 0.05. We found that across HapMap populations, our SNP sample covers at least 60.0% of the genetic variation described in HapMap. In addition, we used 3 SNPs that have not yet been typed in the HapMap project that cover additional, likely overlapping genetic variation of the gene.
We used linear regression to assess whether SNPs in the CRHR1 gene interacted with child abuse to influence BDI scores. For these analyses, a maximum sample size of 422 individuals was informative, but the total sample sizes for the single-SNP analyses vary by a maximum of 5.0% depending on the call rate of each individual SNP. We first considered single-SNP analyses that regressed the BDI outcome (continuous variable) on genotype (coded under an additive model), with child abuse coded into the 2 groups of none to mild and moderate to severe. We further adjusted the regression models for potential confounders including age, sex, and ancestry. We established the significance of genotype–child abuse interaction effects using permutation-based procedures34,35 that randomly assigned the sample BDI scores to subjects (sampled without replacement) while holding each subject's genotype and environmental variables fixed. This permutation method is robust against nonnormal distribution of the outcome variable as we observed with the BDI scores in this population. For each analysis, the empirical P value was based on 10 000 permutations. We conducted these analyses using appropriate components of SAS version 9.1 statistical software.
After completing the single-SNP interaction analyses, we next considered whether SNP-based haplotypes interacted with child abuse to influence BDI scores. Using linear regression, we modeled the BDI score as a function of the various haplotypes (coded under an additive model), child abuse, and the haplotype–child abuse interaction effects. We pooled all of the haplotypes with less than 10.0% estimated frequency into a rare category. We further adjusted the regression models for potential confounders, such as age and sex. We assessed the significance of haplotype-environment interaction effects using permutation procedures that randomly shuffled BDI scores within the sample while holding genotype and environmental data for each subject fixed. For each interaction analysis, we based inference on 100 000 permutations.
We estimated the haplotype frequencies for the 3 SNP haplotypes of LD blocks 1 and 2 (Figure 1A) as well as for the 3 SNP haplotypes consisting of the 3 SNPs with the strongest gene × environment interaction effects (rs7209436 and rs110402 located in block 1 and rs242924 between the 2 blocks). The haplotype block structure was determined using the confidence interval method according to Gabriel et al36 as implemented in the Haploview software. Given that we generated only genotype data, the phase of the underlying haplotypes of a particular subject may be unknown. Therefore, we estimated sample haplotype frequencies using an expectation maximization algorithm37 and used these frequencies to calculate each subject's posterior haplotype-pair probability given his or her genotype data. This was done using the SNPHAP program (http://slack.ser.man.ac.uk/progs/snphap.html). We then applied a conservative haplotype test that used only those subjects whose posterior haplotype-pair probability exceeded 95.0%. However, because the region has a low degree of haplotype complexity (high LD with a few common haplotypes), this excluded only 5.2% of subjects from interaction analysis.
We applied the program SNPSpD (http://genepi.qimr.edu.au/general/daleN/SNPSpD/) to adjust the nominal significance level for the testing of multiple SNP variants within the CRHR1 gene. The program, based on the work of Nyholt,38 applies spectral decomposition techniques to determine the effective number of independent SNPs in regions where LD exists. The software then implements a Bonferroni correction that divides the nominal significance level of .05 by this effective number of independent SNPs. Using this method, we determined that α = .0094 provides an appropriate threshold to declare significance within the CRHR1 gene after correcting for multiple testing.
The second supportive sample was clinical in nature, so we tested the protective effects of the protective alleles and haplotypes against development of MDD as opposed to other types of psychiatric disorders in the presence of moderate to severe child abuse (Table 2). Haplotypes with a posterior haplotype probability greater than 95.0% could only be estimated in 199 women for the block 1 haplotype and 198 women for the 3 best SNP haplotypes. Differences in the SNP genotypes of rs7209436, rs4792887, rs110402, rs242940, and rs242924 and in the copy number frequency of the TCA haplotype, formed by rs7209436, rs4792887, and rs110402, respectively, and of the TAT haplotype, formed by rs7209436, rs110402, and rs242924, respectively, between the groups with no to mild vs moderate to severe child abuse were tested separately for patients with a lifetime diagnosis of MDD and patients with no lifetime diagnosis of MDD using contingency tables and an additive genetic model (linear by linear association). Our specific directional hypotheses were that the protective alleles of these SNPs and both the TCA and TAT haplotypes (which were protective in the African American sample) would also be protective, and thus overrepresented, in women with moderate to severe child abuse but no lifetime diagnosis of MDD. Thus, we report 1-tailed P values for the analyses of these directional hypotheses. The frequencies of the SNPs and for the common haplotypes formed by rs7209436, rs4792887, and rs110402 and by rs7209436, rs110402, and rs242924 in the replication sample are reported in Table 2.
To examine the interaction of the CRHR1 genotype and ELS, we initially analyzed data from our first sample group of 560 low-income, primarily (97.4%) African American adults living in an urban area. Of the 476 subjects with complete CTQ and BDI phenotypes, we first examined the effects of child abuse on adult depression. As expected based on prior literature, higher CTQ abuse scores predicted higher adult BDI scores (F3,472 = 37; P < .001) (Figure 2A). Additionally, when child abuse is divided into the none to mild and moderate to severe groups, we see a significant difference in BDI scores (none to mild: n = 300, mean BDI score, 11.58; moderate to severe: n = 176, mean BDI score, 19.74; P < .001).
The gene encoding CRHR1 is located on chromosome 17q21.31 and contains 13 exons spanning 51 kb. Of the 502 participants in the study from whom we obtained DNA, 422 were genotyped for 15 SNPs located in a 57-kb region spanning the 5′ promoter region to the 3′ end of the CRHR1 gene, with an average intermarker distance of 4.1 kb (Figure 1A). We identified 2 blocks of LD: one within intron 1 spanning at least 8 kb (block 1), and the second from introns 2 through 9 spanning at least 18 kb (block 2). Ten of the SNPs showed a minor allele frequency greater than 5.0% and were thus included in the gene × environment regression analysis described later.
To test the primary hypothesis that genetic polymorphisms in the CRHR1 gene moderate the effects of child abuse on adult depression, we used linear regression–based methods. We first performed single-SNP analyses in which we regressed, using permutation analyses, continuous BDI scores on genotype (coded under an additive model), child abuse (none to mild vs moderate to severe), and the interaction between genotype and child abuse, adjusting for potential confounders such as age and sex. Statistical significance of main effects and interactions were evaluated using a permutation-based analysis (see the “Methods” section). Prior to adjusting for multiple comparisons, 7 of the 10 tested SNPs showed a significant interaction (lowest P = .008 for rs110402) with child abuse for the prediction of adult depression (Figure 1B and C).
To correct for multiple comparisons, we used spectral decomposition techniques to determine the effective number of independent SNPs in regions where LD exists. Using this method, we determined that α = .0094 provides an appropriate threshold to declare significance within the CRHR1 gene. This α level is derived from 5.3 independent tests accounting for the moderate extent of LD across the 10 tested SNPs.
The interactions between child abuse and 2 of the CRHR1 SNPs—rs110402 and rs7209436—remained significant after correction for multiple testing (Figure 2B and C), and 5 others showed significant interaction effects prior to correction (Figure 2D and Figure 3). The rs110402 and rs7209436 SNPs both reside in the first LD block of CRHR1 that spans intron 1 of the gene and are in high LD (r2 = 0.93). For both SNPs, no significant genotype effect on depression was detected in the no to mild abuse group, with the mean BDI scores from these subjects ranging from 10.91 to 12.67 across genotypes (Figure 2B and C). Among individuals in the moderate to severe abuse groups, the rare allele of both SNPs had a protective effect on the severity of adult depressive symptoms. While individuals homozygous for the common alleles (CC for rs7209436 and GG for rs110402) had average BDI scores of 22.33 to 22.49 (reflective of moderate depressive symptoms), heterozygous individuals had significantly lower BDI scores (mean, 16.31-16.39, reflective of mild depressive symptoms) and individuals homozygous for the rare alleles (TT or AA, respectively) had a mean BDI score of 10.22 (reflective of low to absent levels of depressive symptoms) despite the presence of at least 1 type of moderate to severe child abuse. The data suggest an additive protective genetic effect of the rare allele of these 2 SNPs in those exposed to child abuse. We also note that similar interactive gene dosage patterns were seen across the other SNPs that were significant prior to correction, including rs242924 (Figure 2D), rs4792887, rs242940, rs173365, and rs242948 (Figure 3).
A similar permutation-based linear regression approach was applied to investigate gene × environment interactions involving common CRHR1 haplotypes (haplotype frequency > 10.0%) in haplotype blocks 1 and 2 (Figure 1A) as well as a haplotype defined by the 3 SNPs carrying the strongest gene × environment interactions in the individual SNP analysis. In block 1, we observed 3 common haplotypes of approximately equal frequency (30.4%-34.1%) that summed to account for approximately 98.5% of the observed haplotypes (Figure 4). Haplotype interaction analysis showed a significant interaction effect of child abuse and haplotypes within block 1. This association was accounted for by the haplotype formed of alleles T (rs7209436), C (rs4792887), and A (rs110402), which appears to be protective in a dose-dependent manner (P < .001) (Figure 4A and C). The other 2 common haplotypes within block 1 and those within block 2 did not show a significant interaction effect.
We then investigated the haplotypes formed by the 3 most significant SNPs in the single-SNP analysis, of which 2 are located in block 1 (rs7209436 and rs110402) and 1 is between blocks 1 and 2 (rs242924). These 3 SNPs form 2 common “yin and yang” haplotypes in our sample accounting for more than 95.0% of the observed haplotypes (CGG = 66.5% and TAT = 28.8%) (Figure 4B and C). The haplotypes formed by these 3 SNPs showed a significant interaction with child abuse on adult BDI scores (P = .003). The presence of the rarer haplotype appears to decrease the risk of adult depressive symptoms in an additive manner in those with a history of child abuse.
Varying degrees of genetic admixture of sub-Saharan African, European, and Native American ancestry have been reported in African American individuals31,32 and could lead to population stratification and thus spurious association results. To control for potential differences in genetic ancestry, we genotyped 134 AIMs.31,39 Differences in the extent of admixture in our sample do not account for the observed results because genetic ancestry was correlated with neither the extent of child abuse (Pearson correlations of genetic ancestry with CTQ total abuse scores were nonsignificant, ranging from −0.033 to 0.048 for the 4 different tested populations) nor adult BDI scores (Pearson correlations of genetic ancestry with BDI total scores were nonsignificant, ranging from −0.094 to 0.099). To confirm that admixture was not an important confounder in our analysis, we reinvestigated the interaction effect of rs110402, the SNP with the lowest P value, by using linear regression and adding the percentages of sub-Saharan and European ancestry as independent variables. In this model, the interaction effect between the rs110402 genotype and child abuse on BDI scores was less than P = .05 (P = .04 in a smaller sample of 210 individuals). Percentages of sub-Saharan and European ancestry were not significant predictors in this model, with P values of .38 and .80, respectively. Because the β values for the interaction term were not significantly different with and without ancestry in the model (β [SE], 6.56 [2.70] and 5.59 [2.01], respectively), we also concluded that ancestry is not a confounder in our analysis. Overall, the interaction term remained at P ≤ .05 for 5 of the 10 CRHR1 SNPs when adjusting for admixture.
These data suggest that certain alleles of the CRHR1 gene may provide a protective effect on the risk for depression in adults with a history of child abuse. To confirm this hypothesis, we genotyped rs7209436, rs110402, rs4792887, rs242924, and rs242940 in an independent sample that was ethnically (87.7% Caucasian) and socioeconomically (less impoverished) distinct from the original population described earlier. We analyzed the allele distribution of the previously identified protective SNP alleles and haplotypes TCA (rs7209436, rs110402, and rs4792887) (Figure 4A and C) and TAT (rs7209436, rs110402, and rs242924) (Figure 4B and C) in women with and without a lifetime diagnosis of MDD and with and without child abuse (n = 199). For haplotype block 1, we found the following frequencies: TCA, 42.1%; CGG, 43.0%; and CTG, 13.1%. For the best SNP haplotype, we found frequencies of 41.6% for TAT and 56.0% for CGG.
We hypothesized that we would find an overrepresentation of 2 copies of the protective haplotypes in women with a history of moderate to severe child abuse but who had not developed MDD within their lifetimes. Indeed, we found a significant overrepresentation of the protective alleles of the tested SNPs (most significant in replication rs7209436 and rs242940) and of 2 protective copies of both the TCA and TAT haplotypes in this group (Table 2). In the group with child abuse but no MDD, 35.0% carried 2 copies of these haplotypes; however, in the group with no child abuse and no MDD, only 17.0% had 2 copies (Table 2). In contrast, in the group with a lifetime diagnosis of MDD, there were no significant differences in the distribution of haplotype copy numbers regardless of child abuse history.
The group with a history of child abuse but no MDD appears to be protected from depression despite the increased psychosocial risk. Notably, this group also has a relative increase in the copy numbers of protective CRHR1 alleles compared with those with a history of lifetime depression (Figure 4D), underscoring the potential protective effects of the identified genotypes in 2 independent and ethnically different samples.
Our results demonstrate that genetic variants in CRHR1 are moderators of the effects of child abuse on adult depressive symptoms in 2 independent, ethnically different populations. A haplotype formed by 3 SNPs located in intron 1 of the CRHR1 gene is most significantly associated with diminished effects of child abuse on adult depressive symptoms. Similar effects were seen on a single-SNP basis and with haplotypes formed by the 3 most significant SNPs. These genotypes and haplotypes potentially serve as predictors of both risk and resilience for adult depression in men and women with a history of child abuse. The protective effect was seen in a predominantly African American sample as well as in a sample that was predominantly Caucasian, and the effect was observed with both current depressive symptoms and lifetime diagnosis of MDD.
It is unlikely that any of these SNPs are actually functional variants in the CRHR1 gene; instead, they are presumably in LD with an as-yet untyped, potentially functional variant. Given the haplotype block structure of the gene, the pattern of association (Figure 1), and the fact that a haplotype in the first LD block showed the most significant interaction effect (Figure 4), it is likely that the causal variant is located in the 5′ region of the gene (ranging from the first 2 kb of the promoter to intron 2). Increasing evidence suggests that much of the regulation of gene function is within the noncoding regions.40,41 Intron 1 of CRHR1 contains 3 highly conserved regions that may have regulatory functions (according to the UCSC Genome Browser database, http://genome.ucsc.edu/). Functional intronic regulatory elements have been reported for several genes42,43; therefore, these CRHR1 intronic regions could affect transcriptional modulation of gene function. Denser fine mapping of this region combined with resequencing followed by in vitro and in vivo studies may allow for the identification of potentially functional variants.
One potential limitation to this study is the use of BDI scores as a continuous variable to measure current depression symptoms, as our genotype effects measure the interactions with child abuse on adult depression symptoms rather than MDD per se. However, this limitation could also be seen as a strength of the study. As pointed out in the “Methods” section, a significant percentage of our sample received a SCID-I–based DSM diagnosis of MDD, and these data strongly validate our BDI measures of depression. Additionally, current research on the assessment of psychopathology points to a number of advantages of a continuous or dimensional assessment of psychopathology, including depressive symptoms. Other potential limitations are the limited size and different phenotype measure (presence or absence of lifetime MDD) in the second supportive study sample. Despite the limitations of this additional sample, we believe that both having a supporting result within a different socioeconomically based sample and finding the apparent protective effects of the CRHR1 polymorphisms using the SCID-I diagnoses support the main findings from the larger primary sample. In the future, we would also like to know the differential effects of age, sex, and type of abuse among other factors that may mediate these effects.
Several studies suggest that depression-related HPA axis hyperactivity may be related to ELS, which may form the physiological basis for the finding that abuse during early life substantially elevates adult risk for depression.21,22 For example, plasma adrenocorticotropin and cerebrospinal fluid CRH concentrations appear to correlate with perceived ELS more strongly than with current depression severity.44,45 Preclinical studies indicate that the persistent hyperactivity of the HPA axis associated with ELS is mediated by hyperactive CRH neurons, with chronic activation of CRHR1 in limbic brain regions.46,47 It is reasonable to hypothesize that alterations in CRHR1 responsiveness during these early emotional critical periods could alter later risk for HPA axis overactivity and depression.
Our genetic results support the hypothesis that the CRH and CRHR1 system moderates the effects of ELS on adult psychopathology. Our data also point out that genetic association and linkage studies that fail to take the environment into account may miss important genetic variants involved in the etiology of complex diseases. Given the clinical and preclinical data12,13,16,21,45,48 that overactivity of CRH-ergic neurotransmission is associated with depressive and anxiety-like symptoms, it is reasonable to expect that the protective polymorphisms identified here associate with either decreased sensitivity of the CRHR1 or increased negative feedback regulation of its functioning. These results suggest the possibility that heritable differences in CRH-mediated neurotransmission exacerbate or dampen the effects of child abuse on the stress hormone system, potentially modulating HPA axis sensitivity, extrahypothalamic CRH-ergic circuits, and risk for depressive symptoms in adulthood.
Correspondence: Kerry J. Ressler, MD, PhD, Department of Psychiatry and Behavioral Sciences, Yerkes Research Center, Emory University, 954 Gatewood Dr, Atlanta, GA 30329 (email@example.com).
Submitted for Publication: July 26, 2007; final revision received October 5, 2007; accepted October 16, 2007.
Author Contributions: Drs Bradley and Binder contributed equally to this article. Drs Cubells and Ressler are co-senior authors.
Financial Disclosure: Dr Binder receives grant support or has received awards from Pfizer, GlaxoSmithKline, the National Alliance for Research on Schizophrenia and Depression, and the National Institute of Mental Health and has received speaker's honoraria from AstraZeneca. Dr Gillespie has received funding from the American Psychiatric Institute for Research and Education/Wyeth and the National Institute on Drug Abuse. Dr Newport has received research support from Eli Lilly, GlaxoSmithKline, Janssen, the National Alliance for Research on Schizophrenia and Depression, the National Institutes of Health, and Wyeth and has served on speakers' bureaus and/or received honoraria from AstraZeneca, Eli Lilly, GlaxoSmithKline, Pfizer, and Wyeth; he has not served on advisory boards or as a consultant to any biomedical or pharmaceutical corporations, and neither he nor his family members hold equity positions in biomedical or pharmaceutical corporations. Dr Stowe has received grants for clinical research and educational activities, has served as an advisor or consultant to GlaxoSmithKline, Wyeth, and Pfizer, has served as an advisor or consultant for Bristol-Myers Squibb, has served on speakers' bureaus for GlaxoSmithKline, Wyeth, Pfizer, and Eli Lilly, and receives grants from the National Institutes of Health for his clinical research. In the past, Dr Nemeroff consulted to, served on the speakers' bureau and/or board of directors of, has been a grant recipient from, and/or owned equity in 1 or more of the following: Abbott Laboratories, Acadia Pharmaceuticals, the American Foundation for Suicide Prevention, the American Psychiatric Institute for Research and Education, AstraZeneca, BMC-JR LLC, Bristol-Myers-Squibb, CeNeRx, Corcept, Cypress Biosciences, Cyberonics, Eli Lilly, Entrepreneur's Fund, Forest Laboratories, George West Mental Health Foundation, GlaxoSmithKline, i3 DLN, Janssen Pharmaceutica, Lundbeck, the National Alliance for Research on Schizophrenia and Depression, Neuronetics, the National Institute of Mental Health, the National Foundation for Mental Health, NovaDel Pharma, Otsuka, Pfizer Pharmaceuticals, Quintiles, Reevax, UCB Pharma, and Wyeth-Ayerst. Currently, Dr Nemeroff serves on the scientific advisory boards for AstraZeneca, Johnson & Johnson, Pharma Neuroboost, Forest Laboratories, Quintiles, and the National Alliance for Research on Schizophrenia and Depression, is a grant recipient from the National Institutes of Health, the National Alliance for Research on Schizophrenia and Depression, and the American Foundation for Suicide Prevention, serves on the board of directors of the American Foundation for Suicide Prevention, the American Psychiatric Institute for Research and Education, NovaDel Pharmaceuticals, Mount Cook Pharma, Inc, and the George West Mental Health Foundation, owns equity in CeNeRx and Reevax, and owns stock or stock options in Corcept and NovaDel. Dr Ressler has received awards and/or funding support from Lundbeck, the Burroughs Wellcome Fund, Pfizer, the National Alliance for Research on Schizophrenia and Depression, the National Institute of Mental Health, and the National Institute on Drug Abuse and has a consulting agreement with Tikvah Therapeutics for N-methyl-D-aspartate–based therapeutics.
Funding/Support: This work was primarily supported by grant MH071537 from the National Institute of Mental Health. This work was also supported by grant P50-MH 68036 from the Women's Mental Health Program at Emory University; the Emory and Grady Memorial Hospital General Clinical Research Center; grant M01 RR00039 from the National Centers for Research Resources, National Institutes of Health; grants MH-42088 (Dr Nemeroff), MH-58922 (Dr Nemeroff), and MH069884 (Dr Ressler) from the National Institute of Mental Health; grant DA015766 from the National Institute of Drug Abuse (Dr Cubells); and the Burroughs Wellcome Fund (Dr Ressler).
Additional Contributions: Jason Mercer, MS, Taira Everett, BA, Rafael Outland, BA, Abby Powers, BA, Allen Graham, BA, Mark Gapen, BA, and Rachel Hershenberg, BA, provided technical support.