The American Psychiatric Association (APA) has updated its Privacy Policy and Terms of Use, including with new information specifically addressed to individuals in the European Economic Area. As described in the Privacy Policy and Terms of Use, this website utilizes cookies, including for the purpose of offering an optimal online experience and services tailored to your preferences.

Please read the entire Privacy Policy and Terms of Use. By closing this message, browsing this website, continuing the navigation, or otherwise continuing to use the APA's websites, you confirm that you understand and accept the terms of the Privacy Policy and Terms of Use, including the utilization of cookies.

×
Published Online:https://doi.org/10.1176/foc.3.1.156

Abstract

In a prospective-longitudinal study of a representative birth cohort, we tested why stressful experiences lead to depression in some people but not in others. A functional polymorphism in the promoter region of the serotonin transporter (5-HTT) gene was found to moderate the influence of stressful life events on depression. Individuals with one or two copies of the short allele of the 5-HTT promoter polymorphism exhibited more depressive symptoms, diagnosable depression, and suicidality in relation to stressful life events than individuals homozygous for the long allele. This epidemiological study thus provides evidence of a gene-by-environment interaction, in which an individual’s response to environmental insults is moderated by his or her genetic makeup.

Depression is among the top five leading causes of disability and disease burden throughout the world (1). Across the life span, stressful life events that involve threat, loss, humiliation, or defeat influence the onset and course of depression (25). However, not all people who encounter a stressful life experience succumb to its depressogenic effect. Diathesis-stress theories of depression predict that individuals’ sensitivity to stressful events depends on their genetic makeup (6, 7). Behavioral genetics research supports this prediction, documenting that the risk of depression after a stressful event is elevated among people who are at high genetic risk and diminished among those at low genetic risk (8). However, whether specific genes exacerbate or buffer the effect of stressful life events on depression is unknown. In this study, a functional polymorphism in the promoter region of the serotonin transporter gene (SLC6A4) was used to characterize genetic vulnerability to depression and to test whether 5-HTT gene variation moderates the influence of life stress on depression.

The serotonin system provides a logical source of candidate genes for depression, because this system is the target of selective serotonin reuptake–inhibitor drugs that are effective in treating depression (9). The serotonin transporter has received particular attention because it is involved in the reuptake of serotonin at brain synapses (10). The promoter activity of the 5-HTT gene, located on 17q11.2, is modified by sequence elements within the proximal 5’ regulatory region, designated the 5-HTT gene-linked polymorphic region (5-HTTLPR). The short (“s”) allele in the 5-HTTLPR is associated with lower transcriptional efficiency of the promoter compared with the long (“l”) allele (11).

Evidence for an association between the short promoter variant and depression is inconclusive (12). Although the 5-HTT gene may not be directly associated with depression, it could moderate the serotonergic response to stress. Three lines of experimental research suggest this hypothesis of a gene-by-environment (G × E) interaction. First, in mice with disrupted 5-HTT, homozygous and heterozygous (5-HTT −/− and +/−) strains exhibited more fearful behavior and greater increases in the stress hormone (plasma) adrenocorticotropin in response to stress compared to homozygous (5-HTT +/+) controls, but in the absence of stress no differences related to genotype were observed (13). Second, in rhesus macaques, whose length variation of the 5-HTTLPR is analogous to that of humans, the short allele is associated with decreased seroto-nergic function [lower cerebrospinal fluid (CSF) 5-hydroxyindoleacetic acid concentrations] among monkeys reared in stressful conditions but not among normally reared monkeys (14). Third, human neuroimaging research suggests that the stress response is mediated by variations in the 5-HTTLPR. Humans with one or two copies of the s allele exhibit greater amygdala neuronal activity to fearful stimuli compared to individuals homozygous for the l allele (15). Taken together, these findings suggest the hypothesis that variations in the 5-HTT gene moderate psychopathological reactions to stressful experiences.

We tested this G × E hypothesis among members of the Dunedin Multidisciplinary Health and Development Study (16). This representative birth cohort of 1037 children (52% male) has been assessed at ages 3, 5, 7, 9, 11, 13, 15, 18, and 21 and was virtually intact (96%) at the age of 26 years. A total of 847 Caucasian non-Maori study members, without stratification confounds, were divided into three groups on the basis of their 5-HTTLPR genotype (11): those with two copies of the s allele (s/s homozygotes; n=147; 17%), those with one copy of the s allele (s/l heterozygotes; n=435; 51%), and those with two copies of the l allele (l/l homozygotes; n=265; 31%). There was no difference in genotype frequencies between the sexes [χ2(2<)=0.02, P=0.99]. Stressful life events occurring after the 21st birthday and before the 26th birthday were assessed with the aid of a life-history calendar (17), a highly reliable method for ascertaining life-event histories (18). The 14 events included employment, financial, housing, health, and relationship stressors. Thirty percent of the study members experienced no stressful life events; 25% experienced one event; 20%, two events; 11%, three events; and 15%, four or more events. There were no significant differences between the three genotype groups in the number of life events they experienced, F(2,846)=0.56, P=0.59, suggesting that 5-HTTLPR genotype did not influence exposure to stressful life events.

Study members were assessed for past-year depression at age 26 with the use of the Diagnostic Interview Schedule (19), which yields a quantitative measure of depressive symptoms and a categorical diagnosis of a major depressive episode according to Diagnostic and Statistical Manual of Mental Disorders (DSM-IV) criteria (20). 17% of study members (58% female versus 42% male; odds ratio=1.6; 95% confidence interval from 1.1 to 2.2) met criteria for a past-year major depressive episode, which is comparable to age and sex prevalence rates observed in U.S. epidemiological studies (21). In addition, 3% of the study members reported past-year suicide attempts or recurrent thoughts about suicide in the context of a depressive episode. We also collected informant reports about symptoms of depression for 96% of study members at age 26 by mailing a brief questionnaire to persons nominated by each study member as “someone who knows you well.”

We used a moderated regression framework (22), with sex as a covariate, to test the association between depression and (i) 5-HTTLPR genotype, (ii) stressful life events, and (iii) their interaction (table S1). The interaction between 5-HTTLPR and life events showed that the effect of life events on self-reports of depression symptoms at age 26 was significantly stronger (P=0.02) among individuals carrying an s allele than among l/l homozygotes (Fig. 1A). We further tested whether life events could predict within-individual increases in depression symptoms over time among individuals with an s allele by statistically controlling for the baseline number of depressive symptoms they had before the life events occurred (table S1). The significant interaction (P=0.05) showed that individuals carrying an s allele whose life events occurred after their 21st birthday experienced increases in depressive symptoms from the age of 21 to 26 years (b=1.55, SE=0.66, t=2.35, P=0.02 among s/s homozygotes and b=1.25, SE=0.34, t=3.66, P<0.001 among s/l heterozygotes) whereas l/l homozygotes did not (b=0.17, SE=0.41, t=0.41, P=0.68).

The G × E interaction also showed that stressful life events predicted a diagnosis of major depression among carriers of an s allele but not among l/l homozygotes (P=0.056, Fig. 1B). We further tested whether life events could predict the onset of new diagnosed depression among carriers of an s allele (table S1). We excluded from analysis study members who were diagnosed with depression before age 21. The significant interaction (P=0.02) showed that life events occurring after their 21st birthdays predicted depression at age 26 among carriers of an s allele who did not have a prior history of depression (b=0.79, SE=0.25, z=3.16, P=0.002 among s/s homozygotes and b=0.41, SE=0.12, z=3.29, P=0.001 among s/l heterozygotes) but did not predict onset of new depression among l/l homozygotes (b=0.08, SE=0.20, z=0.42, P=0.67). Further analyses showed that stressful life events predicted suicide ideation or attempt among individuals carrying an s allele but not among l/l homozygotes (P=0.05, Fig. 1C). The hypothesized G × E interaction was also significant when we predicted informant reports of age-26 depression (P<0.01), an analysis that ruled out the possibility of self-report bias (Fig. 1D). The interaction showed that the effect of life events on informant reports of depression was stronger among individuals carrying an s allele than among l/l homozygotes. These analyses attest that the 5-HTT gene interacts with life events to predict depression symptoms, an increase in symptoms, depression diagnoses, new-onset diagnoses, suicidality, and an informant’s report of depressed behavior.

This evidence that 5-HTTLPR variation moderates the effect of life events on depression does not constitute unambiguous evidence of a G × E interaction, because exposure to life events may be influenced by genetic factors; if individuals have a heritable tendency to enter situations where they encounter stressful life events, these events may simply be a genetically saturated marker (23, 24). Thus, what we have identified as a gene × environment interaction predicting depression could actually reflect a gene × “gene” interaction between the 5-HTTLPR and other genes we did not measure. We reasoned that, if our measure of life events represents merely genetic risk, then life events would interact with 5-HTTLPR even if they occurred after the depression episode. However, if our measure of life events represents environmental stress, then the timing of life events relative to depression must follow cause-effect order and life events that occur after depression should not interact with 5-HTTLPR to postdict depression. We tested this -hypothesis by substituting the age-26 measure of depression with depression assessed in this longitudinal study when study members were 21 and 18 years old, before the occurrence of the measured life events between the ages of 21 and 26 years. Whereas the 5-HTTLPR × life events interaction predicted depression at the age of 26 years, this same interaction did not postdict depression reported at age 21 nor at the age of 18 years (table S2), indicating our finding is a true G × E interaction.

If 5-HTT genotype moderates the depressogenic influence of stressful life events, it should moderate the effect of life events that occurred not just in adulthood but also of stressful experiences that occurred in earlier developmental periods. Based on this hypothesis, we tested whether adult depression was predicted by the interaction between 5-HTTLPR and childhood maltreatment that occurred during the first decade of life (16, 25). Consistent with the G × E hypothesis, the longitudinal prediction from childhood maltreatment to adult depression was significantly moderated by 5-HTTLPR (table S3). The interaction showed (P=0.05) that childhood maltreatment predicted adult depression only among individuals carrying an s allele but not among l/l homozygotes (Fig. 2).

We previously showed that variations in the gene encoding the neurotransmitter-metabolizing enzyme monoamine oxidase A (MAOA) moderate children’s sensitivity to maltreatment (25). MAOA has high affinity for 5-HTT, raising the possibility that the protective effect of the l/l allele on psychiatric morbidity is further augmented by the presence of a genotype conferring high MAOA activity (13, 26). However, we found that the moderation of life stress on depression was specific to a polymorphism in the 5-HTT gene, because this effect was observed regardless of the individual’s MAOA gene status (tables S4 and S5).

Until this study’s findings are replicated, speculation about clinical implications is premature. Nonetheless, although carriers of an s 5-HTTLPR allele who experienced four or more life events constituted only 10% of the birth cohort, they accounted for almost one-quarter (23%) of the 133 cases of diagnosed depression. Moreover, among cohort members suffering four or more stressful life events, 33% of individuals with an s allele became depressed, whereas only 17% of the l/l homozygotes developed depression (Fig. 3). Thus, the G × E’s attributable risk and predictive sensitivity indicate that more knowledge about the functional properties of the 5-HTT gene may lead to better pharmacological treatments for those already depressed. Although the short 5-HTTLPR variant is too prevalent for discriminatory screening (over half of the Caucasian population has an s allele), a microarray of genes might eventually identify those needing prophylaxis against life’s stressful events (27).

Evidence of a direct relation between the 5-HTTLPR and depression has been inconsistent (12), perhaps because prior studies have not considered participants’ stress histories. In this study, no direct association between the 5-HTT gene and depression was observed. Previous experimental paradigms, including 5-HTT knockout mice (13), stress-reared rhesus macaques (14), and human functional neuroimaging (15), have shown that the 5-HTT gene can interact with environmental conditions, although these experiments did not address depression. Our study demonstrates that this G × E interaction extends to the natural development of depression in a representative sample of humans. However, we could not test hypotheses about brain endophenotypes (28) intermediate between the 5-HTT gene and depression because of the difficulty of taking CSF or functional magnetic resonance imaging measurements in an epidemiological cohort.

Much genetic research has been guided by the assumption that genes cause diseases, but the expectation that direct paths will be found from gene to disease has not proven fruitful for complex psychiatric disorders (29). Our findings of G × E interaction for the 5-HTT gene and another candidate gene, MAOA (25), point to a different, evolutionary model. This model assumes that genetic variants maintained at high prevalence in the population probably act to promote organisms’ resistance to environmental pathogens (30). We extend the concept of environmental pathogens to include traumatic, stressful life experiences and propose that the effects of genes may be uncovered when such pathogens are measured (in naturalistic studies) or manipulated (in experimental studies). To date, few linkage studies detect genes, many candidate gene studies fail consistent replication, and genes that replicate account for little variation in the phenotype (29). If replicated, our G × E findings will have implications for improving research in psychiatric genetics. Incomplete gene penetrance, a major source of error in linkage pedigrees, can be explained if a gene’s effects are expressed only among family members exposed to environmental risk. If risk exposure differs between samples, candidate genes may fail replication. If risk exposure differs among participants within a sample, genes may account for little variation in the phenotype. We speculate that some multifactorial disorders, instead of resulting from variations in many genes of small effect, may result from variations in fewer genes whose effects are conditional on exposure to environmental risks.

Figure 1.

Figure 1. Results of multiple regression analyses estimating the association between number of stressful life events (between ages 21 and 26 years) and depression outcomes at age 26 as a function of 5-HTT genotype. Among the 146 s/s homozygotes, 43 (29%), 37(25%), 28 (19%), 15 (10%), and 23 (16%) study members experienced zero, one, two, three, and four or more stressful events, respectively. Among the 435 s/l heterozygotes, 141 (32%), 101 (23%), 76 (17%), 49 (11%), and 68 (16%) experienced zero, one, two, three, and four or more stressful events. Among the 264 l/l homozygotes, 79 (29%), 73 (28%), 57 (21%), 26 (10%), and 29 (11%) experienced zero, one, two, three, and four or more stressful events. (A) Self-reports of depression symptoms. The main effect of 5-HTTLPR (i.e., an effect not conditional on other variables) was marginally significant (b=−0.96, SE=0.52, t=1.86, P=0.06), the main effect of stressful life events was significant (b=1.75, SE=0.23, t=7.45, P<0.001), and the interaction between 5-HTTLPR and life events was in the predicted direction (b=−0.89, SE=0.37, t=2.39, P=0.02). The interaction showed that the effect of life events on self-reports of depression symptoms was stronger among individuals carrying an s allele (b=2.52, SE=0.66, t=3.82, P<0.001 among s/s homozygotes, and b=1.71, SE=0.34, t=5.02, P<0.001 among s/l heterozygotes) than among l/l homozygotes (b=0.77, SE=0.43, t=1.79, P=0.08). (B) Probability of major depressive episode. The main effect of 5-HTTLPR was not significant (b=−0.15, SE=0.14, z=1.07, P=0.29), the main effect of life events was significant (b=0.37, SE=0.06, z=5.99, P<0.001), and the G × E was in the predicted direction (b=−0.19, SE=0.10, z=1.91, P=0.056). Life events predicted a diagnosis of major depression among s carriers (b=0.52, SE=0.16, z=3.28, P=0.001 among s/s homozygotes, and b=0.39, SE=0.09, z=4.24, P<0.001 among s/l heterozygotes) but not among l/l homozygotes (b=0.16, SE=0.13, z=1.18, P=0.24). (C) Probability of suicide ideation or attempt. The main effect of 5-HTTLPR was not significant (b=−0.01, SE=0.28, z=0.01, P=0.99), the main effect of life events was significant (b=0.51, SE=0.13, z=3.96, P<0.001), and the G × E interaction was in the predicted direction (b=−0.39, SE=0.20, t=1.95, P=0.051). Life events predicted suicide ideation or attempt among s carriers (b=0.48, SE=0.29, z=1.67, P=0.09 among s/s homozygotes, and b=0.91, SE=0.25, z=3.58, P<0.001 among s/l heterozygotes) but not among l/l homozygotes (b=0.13, SE=0.26, z=0.49, P=0.62). (D) Informant reports of depression. The main effect of 5-HTTLPR was not significant (b=−0.06, SE=0.06, t=0.98, P=0.33), the main effect of life events was significant (b=0.23, SE=0.03, t=8.47, P<0.001), and the G × E was in the predicted direction (b=−0.11, SE=0.04, t=2.54, P<0.01). The effect of life events on depression was stronger among s carriers (b=0.39, SE=0.07, t=5.23, P<0.001 among s/s homozygotes, and b=0.17, SE=0.04, t=4.51, P<0.001 among s/l heterozygotes) than among l/l homozygotes (b=0.14, SE=0.05, t=2.69, P<0.01).

Figure 2.

Figure 2. Results of regression analysis estimating the association between childhood maltreatment (between the ages of 3 and 11 years) and adult depression (ages 18 to 26), as a function of 5-HTT genotype. Among the 147s/s homozygotes, 92 (63%), 39 (27%), and 16 (11%) study members were in the no maltreatment, probable maltreatment, and severe maltreatment groups, respectively. Among the 435 s/l heterozygotes, 286 (66%), 116 (27%), and 33 (8%) were in the no, probable, and severe maltreatment groups. Among the 265 l/l homozygotes, 172 (65%), 69 (26%), and 24 (9%) were in the no, probable, and severe maltreatment groups. The main effect of 5-HTTLPR was not significant (b=−0.14, SE=0.11, z=1.33, P=0.19), the main effect of childhood maltreatment was significant (b=0.30, SE=0.10, z=3.04, P=0.002), and the G × E interaction was in the predicted direction (b=−0.33, SE=0.16, z=2.01, P=0.05). The interaction showed that childhood stress predicted adult depression only among individuals carrying an s allele (b=0.60, SE=0.26, z=2.31, P=0.02 among s/s homozygotes, and b=0.45, SE=0.16, z=2.83, P=0.01 among s/l heterozyotes) and not among l/l homozygotes (b=−0.01, SE=0.21, z=0.01, P=0.99).

Figure 3.

Figure 3. The percentage of individuals meeting diagnostic criteria for depression at age 26, as a function of 5-HTT genotype and number of stressful life events between the ages of 21 and 26. The figure shows individuals with either one or two copies of the short allele (left) and individuals homozygous for the long allele (right). In a hierarchical logistic regression model, the main effect of genotype (coded as s group=0 and l group=1) was not significant, b=−0.15, SE=0.21, z=0.72, P=0.47; the main effect of number of life events was significant, b=0.34, SE=0.06, z=5.70, P<0.001; and the interaction between genotype and number of life events was significant, b=−0.30, SE=0.15, z=1.97, P=0.05.

(Reprinted with permission from Science 2003 301(5631):386–389. Copyright 2003 American Association for the Advancement of Science)

Supporting online materialwww.sciencemag.org/cgi/content/full/301/5631/386/DC1 Materials and Methods Tables S1 to S5

References and Notes

1 C. J. Tang, A. D. Lopez, Lancet 349, 1498 (1997).CrossrefGoogle Scholar

2 G. W. Brown, Soc. Psychiatry Psychiatr. Epidemiol. 33, 363 (1998).Google Scholar

3 K. S. Kendler, L. M. Karkowski, C. A. Prescott, Am. J. Psychiatry 156, 837 (1999).CrossrefGoogle Scholar

4 R. C. Kessler, Annu. Rev. Psychol. 48, 191 (1997).CrossrefGoogle Scholar

5 D. S. Pine, P. Cohen, J. G. Johnson, J. S. Brook, J. Affect. Disorders 68, 49 (2002).CrossrefGoogle Scholar

6 E. J. Costello et al., Biol. Psychiatry 52, 529 (2002).CrossrefGoogle Scholar

7 S. M. Monroe, A. D. Simons, Psychol. Bull. 110, 406 (1991).CrossrefGoogle Scholar

8 K. S. Kendler et al., Am. J. Psychiatry 152, 833 (1995).CrossrefGoogle Scholar

9 C. A. Tamminga et al., Biol. Psychiatry 52, 589 (2002).CrossrefGoogle Scholar

10 K. P. Lesch, M. D. Greenberg, J. D. Higley, A. Bennett, D. L. Murphy, in Molecular Genetics and the Human Personality, J. Benjamin, R. P. Ebstein, R. H. Belmaker, Eds. [American Psychiatric Association (APA), Washington, DC, 2002], pp. 109–136.Google Scholar

11 K. P. Lesch et al., Science 274, 1527 (1996).CrossrefGoogle Scholar

12 K. P. Lesch, in Behavioral Genetics in the Postgenomics Era, R. Plomin, J. C. DeFries, I. W. Craig, P. McGuffin, Eds. (APA, Washington, DC, 2003), pp. 389–424.Google Scholar

13 D. L. Murphy et al., Brain Res. Bull. 56, 487 (2001).CrossrefGoogle Scholar

14 A. J. Bennett et al., Mol. Psychiatry 7, 188 (2002).Google Scholar

15 A. R. Hariri et al., Science 297, 400 (2002).CrossrefGoogle Scholar

16 Materials and methods are available as supporting material on Science Online.Google Scholar

17 A. Caspi et al., Int. J. Methods Psychiatr. Res. 6, 101 (1996).CrossrefGoogle Scholar

18 R. F. Belli, W. L. Shay, F. P. Stafford, Public Opin. Q. 65, 45 (2001).CrossrefGoogle Scholar

19 L. N. Robins, L. Cottler, K. Bucholtz, W. Compton, Diagnostic Interview Schedule for DSM-IV (Washington University, St. Louis, MO, 1995).Google Scholar

20 APA, Diagnostic and Statistical Manual of Mental Disorders (APA, Washington, DC, ed. 4, 1994).Google Scholar

21 R. C. Kessler, K. A. McGonagle, M. Swartz, D. G. Blazer, C. B. Nelson, J. Affect. Disorders 29, 85 (1993).CrossrefGoogle Scholar

22 L. S. Aiken, S. G. West, Multiple Regression: Testing and Interpreting Interactions (Sage, Thousand Oaks, CA, 1991).Google Scholar

23 K. S. Kendler, L. Karkowski-Shuman, Psychol. Med. 27, 539 (1997).CrossrefGoogle Scholar

24 R. Plomin, C. S. Bergeman, Behav. Brain Sci. 14, 373 (1991).CrossrefGoogle Scholar

25 A. Caspi et al., Science 297, 851 (2002).CrossrefGoogle Scholar

26 N. Salichon et al., J. Neurosci. 21, 884 (2001).CrossrefGoogle Scholar

27 W. E. Evans, M. V. Relling, Science 286, 487 (1999).CrossrefGoogle Scholar

28 I. I. Gottesman, T. D. Gould, Am. J. Psychiatry 160, 636 (2003).CrossrefGoogle Scholar

29 D. Hamer, Science 298, 71 (2002).CrossrefGoogle Scholar

30 A. V. S. Hill, Br. Med. Bull. 55, 401 (1999).CrossrefGoogle Scholar

31 We thank P. Silva, founder of the Dunedin Multidisciplinary Health and Development Study, Air New Zealand, and the study members, their families, and their friends. Supported by the Health Research Council of New Zealand and the University of Wisconsin Graduate School and by grants from the U.K. Medical Research Council, the William T. Grant Foundation, and the U.S. National Institute of Mental Health (MH49414 and MH45070). T.E.M. is a Royal Society–Wolfson Research Merit Award holder. The study protocol was approved by the institutional review boards of the participating universities.Google Scholar