Our results provide evidence that CSA-associated risk for alcohol-related outcomes is moderated by the H2 haplotype. We first confirmed that a history of CSA is associated with significant risk for alcohol consumption and lifetime DSM-IV alcohol dependence in our sample; the H2 haplotype was not associated with risk for either of these outcomes. We found that, for both alcohol consumption and dependence, a G × E interaction involving CSA and the H2 haplotype is associated with significant protective effects. In those with the H2 haplotype, we observed no significant CSA-associated risk for either higher alcohol consumption or alcohol dependence.
Our results extend and clarify a prior report (
Blomeyer et al. 2008) in which a significant G × E interaction involving severe stressors over the prior three years and the minor allele of rs1876831, a
CRHR1 SNP, led to protection against alcohol consumption in adolescents. We found evidence of a similar interaction and demonstrated that it involves the ~1.5 Mb H2 haplotype spanning several genes in this region of chromosome 17. rs1876831 is one of many informative markers that can be used to tag this haplotype. Another important feature that distinguishes our study from the prior report is that our larger, substantially older sample has already passed through the period of greatest risk for problematic alcohol use enabling an examination lifetime measures, peak ACFS and DSM-IV alcohol dependence diagnosis, that characterize respondents’ mature drinking patterns. Both of these measures have excellent psychometric properties and have been shown to be at least moderately heritable (
Heath et al. 1997;
Bucholz et al. 1994;
Agrawal et al. 2009;
Grant et al. 2009). Our findings demonstrate that the protection against CSA-associated risk for problematic alcohol use associated with the H2 haplotype is persistent into adulthood. Overall, our results suggest that one or more of the 5 adjacent genes within the H2 haplotype play a major role in the risk for alcohol-related outcomes associated with severe life stressors.
Additional research will be necessary to clarify which gene or genes are responsible for these protective effects. Animal studies (
Le et al. 2000;
Shaham et al. 2000;
Funk et al. 2003;
Hansson et al. 2006;
Sommer et al. 2008;
Heilig & Koob, 2007;
Pastor et al. 2008;
Funk et al. 2003) published to date provide evidence that alterations in
CRHR1 expression are involved in important facets of alcohol consumption, dependence, and relapse. Marchigian-Sardinian Preferring (msP) rats, selectively bred for high alcohol preference, have a
CRHR1 promoter polymorphism that results in increased
CRHR1 expression in limbic regions (
Hansson et al. 2006). Alcohol self-administration by non-dependent msP rats is suppressed by CRF
1 antagonists; similar effects are not seen in unselected Wistar (control) rats. msP rats also display greater sensitivity to inhibition of foot-shock-induced ethanol reinstatement by a CRF
1 antagonist than do control rats.(
Hansson et al. 2006) Rodents (not limited to those bred for ethanol preference) experience a prolonged period of increased anxiety and stress responsivity when withdrawn from ethanol after having been made dependent, symptoms that can be blocked by administration of a CRF
1 antagonist (
Sommer et al. 2008). In these animals, upregulation of
CRHR1 expression has been found in the basolateral and medial amygdalar nuclei (
Sommer et al. 2008) as have elevated CRF levels in their central amygdalar nucleus. In rats previously dependent on ethanol (similar findings have been reported with heroin and cocaine) (
Shaham et al. 2000), relapse can be induced by foot-shock stress. For each drug, this reinstatement can be blocked by intracerebroventricular infusion of CRF
1 antagonists (
Le et al. 2000;
Shaham et al. 2000). Direct CRF administration into involved brain regions also induces relapse in these animals (
Le et al. 2000) and has been shown to block the increase in
c-fos expression in the central nucleus of the amygdala that follows footshock (
Funk et al. 2003).
A recent report (
Barr et al. 2009) provides evidence that findings from the rodent literature are also applicable to non-human primates. The authors describe a functional
CRH promoter polymorphism in rhesus macaques that is associated with greater stress responsivity and report interactions between genotype and rearing condition with significantly greater ACTH and cortisol responses to social stress observed in peer-reared carriers of this polymorphism. Most importantly, they report a genotype-environment interaction in which significantly greater alcohol consumption was found only for carriers of this polymorphism who had experienced early life stress (peer-rearing). This interaction, involving a
CRH promoter polymorphism and a severe early life stressor associated with greater alcohol consumption in the absence of a significant main effect, is strikingly similar to the findings of the current report.
Other studies in rodents have demonstrated that exposure to severe, early life stressors results in epigenetically-mediated alterations in gene expression (
Weaver et al. 2004;
Weaver et al. 2005;
Champagne et al. 2006;
Weaver et al. 2007). Although admittedly speculative at present, it is possible that the protective effect which we observed could result from one or more functional polymorphisms limiting CSA-related epigenetically-mediated changes. Epigenetic modifications of gene expression due to ethanol are known to occur through a variety of different mechanisms (
Pandey et al. 2008;
Pietrzykowski et al. 2008;
Shukla et al. 2008) and thus might not be susceptible to a similar protective effect (consistent with the lack of a significant main effect for H2 haplotype in our sample).
Our failure to confirm the earlier report (
Treutlein J et al. 2006) of a significant main effect for rs1876831 on lifetime prevalence of (any) binge drinking and drunkenness in their adolescent sample and on alcohol intake in clinically-ascertained alcoholic adults may have resulted from disparity in the two studies’ outcome measures. The ACFS (
Agrawal et al. 2009;
Grant et al. 2009) used in this study is an estimate of consumption that is primarily based on the period of heaviest lifetime use. The main effects reported in the adolescent sample (
Treutlein et al. 2006) were for two very early drinking career milestones, any episode of binge drinking or drunkenness. In their clinical sample, the significant main effect was for consumption of >250 g ethanol daily prior to admission, a binary measure of daily consumption much later in the course of the disorder that could be affected by various covariates (eg, gender,
ADH genotype, alcoholic liver disease). Another study (
Dahl et al. 2005) examined whether
CRHR1 polymorphisms were associated with alcohol dependence risk in European American alcoholics found no association for any SNPs including those informative for the H2 haplotype.
Several additional issues should be considered when interpreting our results. We substituted CSA for the measure of multiple severe stressors in adolescence used in the prior report (
Blomeyer et al. 2008). Previous reports of G × E interactions have included either multiple classes of stressors (
Caspi et al. 2003) or substituted stressors for which data were available in replication studies (
Kaufman et al. 2004;
Gillespie et al. 2005;
Kendler et al. 2005;
Surtees et al. 2006). Additional study will be necessary to delineate whether our findings will generalize to different types of stressors or those occurring during other developmental periods. CSA is also associated (
Nelson et al. 2002) with risk for subsequent stressors further complicating definitive attribution. Although use of binary outcome variables for examination of G × E interactions may predispose to spurious positive findings (
Eaves, 2006), the association we observed with ACFS, a continuous measure, is less susceptible to this source of error. Given concerns about the sensitivity of inferences from G × E interactions to scaling, we ran an additional analysis using Huber robust regression which down-weights outlier observations: the G × E interaction term remained significant (p=0.026) with a modest reduction in effect size (beta= -0.37 versus beta= -0.41 in the linear regression analyses). The association with alcohol-related phenotypes could be bidirectional, however, the mean age (11.0y) of first CSA occurrence suggests that CSA typically precedes alcohol problems. The use of a binary variable for CSA combines diverse abuse experiences encompassing a wide range of severity; a stronger association might have been observed if a quantitative covariate incorporating information on abuse severity and duration were available. Retrospective assessment of CSA in adults may raise concerns regarding the introduction of bias. The CSA measure used in this study has been found to have reasonable concordance within female like-sex twin pairs and significant association with psychiatric sequelae and parental risk factors (
Dinwiddie et al. 2000;
McLaughlin et al. 2000). Examination of discordant pairs found no evidence for retrospective bias in the association of CSA with parental rejection (
McLaughlin et al. 2000). The prevalence of CSA in our sample is higher than is typical for community samples, a likely result of enrichment for heavy smoking. Because of this enrichment, additional studies will be necessary to determine whether our results will generalize to samples representative of the general population. Post-hoc analyses of genetically informative data on ancestry from a partially overlapping GWAS revealed that two families (N=4 individuals in total) are outliers on the basis of ancestry. The ACFS of these families does not differ significantly from the remainder of the sample; dropping these 4 individuals has very minimal effect on either interaction term point estimates or p values.
Our data suggest that the H2 haplotype is protective against CSA-associated risk for higher lifetime alcohol consumption and alcohol dependence. The extent of these protective effects suggests that one or more of the genes within the H2 haplotype are playing an important role in stress-associated risk for alcohol consumption and dependence. Although the evidence (
Le et al. 2000;
Hansson et al. 2006;
Sommer et al. 2008;
Heilig & Koob, 2007;
Pastor et al. 2008) from animal research is quite strong, it is premature to conclude that these protective effects are due to a
CRHR1 polymorphism. Additional research is needed to determine definitively the gene or genes responsible for these protective effects. Gene expression studies have provided strong evidence in favor of
MAPT (
Caffrey et al. 2006;
Caffrey et al. 2008) versus
CRHR1 (
Campdelacreu et al. 2006) involvement in PSP risk, although other research (
Cruchaga et al. 2009) suggests that multiple genes may be involved. However, since our finding involves a G × E interaction rather than a main effect, it will be considerably more difficult to conduct this type of investigation. If the effects that we observe are a consequence of a
CRHR1 polymorphism, our results may have immediate clinical relevance. Researchers have recently developed (
Gehlert et al. 2007) improved CRF
1 antagonist drugs which are already scheduled for pharmaceutical company-sponsored clinical trials of alcohol dependence treatment. Our findings would predict that H1 homozygotes with a history of severe early trauma exposure will preferentially respond to these agents. More generally, our findings emphasize the potential utility of screening for severe early trauma exposure in individuals presenting for alcohol dependence treatment.