- Research article
- Open Access
- Open Peer Review
Serotonin transporter gene polymorphism (5-HTTLPR) L allele interacts with stress to increase anxiety symptoms in Chinese adolescents: a multiwave longitudinal study
BMC Psychiatryvolume 15, Article number: 248 (2015)
Previous studies of the interaction between a functional polymorphism in the serotonin transporter gene-linked promoter region (5-HTTLPR) and stress in anxiety-related phenotypes have produced inconsistent results. The aim of the current study was to examine the effect of the 5-HTTLPR × stress interaction on anxiety symptoms in Chinese adolescents.
A total of 651 healthy adolescents [323 females and 328 males; age 14–17 (mean = 16.27, standard deviation = 0.77)] participated in this study. At the initial assessment, participants completed self-report measures assessing anxiety symptoms, depressive symptoms and stressful life events. Additionally, anxiety symptoms and stressful life events were assessed once every 3 months for the subsequent 9 months. A hierarchical linear model (HLM) was used to investigate the 5-HTTLPR × stress interaction.
The HLM indicated no main effect of 5-HTTLPR on anxiety symptoms. Significant 5-HTTLPR × stress interaction effect in predicting anxiety symptoms was found. Specifically, individuals with the 5-HTTLPR L allele exhibited more anxiety symptoms related to stressful life events.
The association between stress and anxiety symptoms is moderated by 5-HTTLPR. The 5-HTTLPR L allele increases individuals’ vulnerability to anxiety under stress situations.
Anxiety disorders are complex mental disorders that place an enormous burden on individuals and society . Environmental stress, often indexed by stressful life events (SLEs), has been consistently linked to anxiety. Nevertheless, the effects of SLEs on anxiety differ among individuals, and not all people who experience SLEs develop anxiety disorders , suggesting the existence of variation in individuals’ susceptibility to stress.
Susceptibility to stress may have biological roots, especially in the serotonergic system. Serotonin [5-hydroxytryptamine (5-HT)] is an important monoamine neurotransmitter involved in the regulation of several psychological activities, such as emotion, cognition, circadian and neuroendocrine rhythms (food intake, sleep and sexual activity) . Among genes of the 5-HT system, the serotonin transporter (5-HTT) gene has captured particular attention because 5-HTT is involved in the reuptake of 5-HT at brain synapses. The 5-HTT gene has several polymorphisms, including a functional polymorphism consisting of a 44–base pair insertion/deletion in the 5′ promoter region, known as the serotonin transporter–linked promoter region polymorphism (5-HTTLPR). 5-HTTLPR includes long (L) and short (S) alleles that influence the rate of 5-HT transcription. Specifically, the S allele has lower transcriptional efficiency compared with the L allele .
Focusing on 5-HTT, Lesch et al.  first linked anxiety-related traits with the S allele of 5-HTTLPR. However, this association was observed inconsistently in subsequent studies of anxiety-related phenotypes and was difficult to replicate [5–7]. A possible explanation for the mixed results is that genetic vulnerability might manifest in subjective markers of anxiety, but only when activated by environmental factors, such as stress. In recent years, researchers have paid attention to the interplay between genetic and environmental (G × E) factors in anxiety. Several studies have investigated the 5-HTTLPR × stress interaction in anxiety-related phenotypes. Unlike risk role of 5-HTTLPR S allele in most depression studies , there was no consolidated conclusion about 5-HTTLPR × stress interaction in anxiety to date. Gunthert et al.  and Petersen et al.  observed that individuals with at least one S allele who experienced more stressors reported more anxious mood. Conversely, Grabe et al.  found an interactive effect of the higher active L allele and traumatic events on anxiety disorder. However, a few studies indicated no interaction effect between 5-HTTLPR and SLEs on anxious/depressed symptoms or anxiety sensitivity [12, 13]. These conflicting results for the 5-HTTLPR × stress interaction may be due to differences in sample age, ethnicity (e.g., 5-HTTLPR allele frequency varies considerably according to ethnic background), as well as the methods used to rate predictors and outcomes.
Beyond the inconsistent findings, previous studies have several major limitations. Most studies conducted to test G × E factor interaction in anxiety used cross-sectional designs. Without longitudinal measures of predictors and outcomes, researchers cannot determine whether this interaction predicts the development of anxiety over time. Furthermore, the majority of previous studies focused predominantly on adult samples, few studies have thoroughly examine these effects and interactions in adolescence, one of the most important periods for the onset of anxiety. Adolescence is marked by major changes in physiology, interpersonal relationships and social interests. Because of unique stressors related to school work and socialization, anxiety symptoms and the prevalence of anxiety disorders increase substantially in adolescence . Thus, adolescence is a crucial period to study the pathogenesis of genetic and environmental factors’ effects on anxiety. Moreover, few studies have considered the effects of depressive symptoms on anxiety symptoms. As depression and anxiety are closely related in children and adolescents , it is necessary to control for depressive symptoms when testing G × E interaction in anxiety. In addition, 5-HTTLPR × stress studies have yielded inconsistent results in different ethnicities [16, 17]. Due to ethnical differences, it is necessary to test whether the conclusions of 5-HTTLPR × E studies among other ethnicities may apply to Chinese populations.
To the best of our knowledge, no study has investigated how 5-HTTLPR and SLEs may interact to predict changes in anxiety symptoms over time among Chinese adolescents. Increases in anxiety symptoms over time are expected to be related to both changing levels of stress and individuals’ genetic predispositions. We have reported interaction effects of 5-HTTLPR × stress on depression in Chinese adolescent girls in previous study . In the present study, we thus used the same longitudinal design to investigate whether the interaction between 5-HTTLPR and SLEs can predict anxiety symptoms in Chinese adolescents.
All subjects were volunteers recruited from two public senior high schools in Hunan province of China by posters and media advertisements. Participants with neurological diseases, and/or past or current episodes of anxiety disorders, major depression disorder, manic disorder, bipolar disorder, schizoaffective disorder or schizophrenia were excluded. Of the 692 adolescents administered the clinical interviews, 41 were excluded for meeting the criteria for lifetime major depressive disorder (9, 1.3 %), generalized anxiety disorder (11, 1.6 %), compulsive disorder (16, 2.3 %) or specific phobia (5, 0.7 %). Finally, the study sample consisted of 651 healthy students (323 females and 328 males) aged 14–17 (mean = 16.27, standard deviation = 0.77). All subjects were of Han ancestry, the predominant ethnic group in China.
This study complied with the Code of Ethics of the World Medical Association and was approved by the Ethics Committee of Second Xiangya Hospital, Central South University. All participants and their parents received detailed information and provided written informed consents. Trained researchers who were graduate students at Second Xiangya Hospital performed clinical assessments and administered questionnaires. Neurological physical examination and the clinical interviews were conducted one-on-one with each participant outside of class time. Interviews consisted of two parts. The Anxiety Disorders Interview Schedule for Children for DSM-IV was used to assess anxiety disorders . Affective disorder, schizoaffective disorder and schizophrenia were diagnosed using the Schedule for Affective disorder and Schizophrenia for School-Age Children . At initial assessment, each subject completed the Chinese version of Multidimensional Anxiety Scale for Children (MASC-C) , Center for Epidemiological Studies Depression Scale (CES-D) , and Adolescent Life Events Questionnaire (ALEQ) . In addition, an oral swab was collected from each participant for 5-HTTLPR genotyping during the initial assessment. The researchers returned to the schools to conduct assessments every 3 months for the subsequent 9 months (i.e., at 3, 6, 9 months). Participants completed the MASC-C and ALEQ at each follow-up assessment.
Multidimensional anxiety scale for children
The MASC-C is a 39-item measure that assesses the severity of anxious symptoms in the past week . Each item consists of a statement that youth rate on a four-point Likert scale ranging from 0 (never applies to me) to 3 (often applies to me). Total scores range from 0 to 117, with higher scores indicating higher levels of anxiety symptoms. The MASC-C has been demonstrated to be reliable and valid . In the current study, Cronbach’s alpha values for the scale ranged from 0.90 to 0.95 across administrations, indicating strong internal consistency.
Center for epidemiological studies depression scale
The CES-D is a 20-item measure designed to assess the current level of depressive symptoms, with emphasis on the affective component, and depressive mood in general populations . Each item consists of one symptom. Participants rated the frequency of each symptom within the past week on a four-point scale (0, <1 day; 1, 1–2 days; 2, 3–4 days; 3, 5–7days). Total scores range from 0 to 60, with higher scores indicating higher elevations in depressive symptoms. The Chinese version of the CES-D has shown high degrees of reliability and validity . The Cronbach’s alpha value for the scale was 0.89 at the initial assessment, indicating high internal consistency.
Adolescent life events questionnaire
The ALEQ is a self-reported questionnaire developed by Hankin and Abramson  to assess a broad range of negative life events (e.g., school/achievement problems, friendship and romantic problems, and family problems) that occur during adolescence. Participants rated the frequency of negative life events within the past month on a five-point scale (1, never; 2, rarely; 3, sometimes; 4, usually; 5, always). Total scores range from 70 to 350, with higher scores reflecting a greater number of negative life events. Past research has found the ALEQ to be reliable and valid when used in Chinese adolescents . In the current study, Cronbach’s alpha values for this scale ranged from 0.95 to 0.97 across administrations, indicating high internal consistency.
Genomic DNA was extracted from exfoliated buccal cells using the TIANamp Swab DNA Kit (TIANGEN Biotech, Beijing, China) according to standard procedures. 5-HTTLPR genotyping was performed using the primers described by Heils et al.  (forward: 5′-GGCGTTGCCGCTCTGAATTGC-3′; reverse: 5′-GAGGGACTGAGCTGGACAACCCAC-3′). Polymerase chain reaction (PCR) amplification was conducted using a Perkin-Elmer GeneAmp PCR System 2400 (Applied Biosystems, USA). Amplification system was performed in a volume of 25 μL containing 50ngDNA template, 9.5 μL nuclease-free water, 0.4 μM each primer, 0.2 mM dNTPs, 10 mM Tris–HCl (pH = 8.3), 50 mM KCl, 1.5 mM MgCl2, and 1U GoTaq DNA polymerase (Promega, USA). The cycling conditions were: (1) initial denaturation at 94 °C for 3 min, (2) 35 cycles of amplification (denaturation at 95 °C for 30s, annealing at 62 °C for 30s, and synthesis at 72 °C for 45 s), and (3) final extension at 72 °C for 7 min. The amplification products were resolved on a 1.5 % agarose gel by electrophoresis and visualized by Du Red staining (Biosharp, USA). Fragment sizes were determined on a Bio-Rad Gel Doc XR+ system (Bio-Rad, USA) by comparison with molecular length standards (50 bp ladder, TIANGEN Biotech). The results were fully validated in approximate 5 % of genotyped individuals who were randomly selected for retesting through the same procedure.
A hierarchical linear model (HLM) was used to investigate whether interaction between 5-HTTLPR and stress can predict the level of anxiety symptoms. Analyses were carried out using the SAS (version 9.0, SAS Institute) MIXED procedure and maximum likelihood estimation. A two-level structure of model was handled. Level 1 was 4-time follow-up measurements within each subject (fluctuation in anxiety and stress over times within each subject), and level 2 was all subjects (different 5-HTTLPR genotypes between subjects). The dependent variable was within-subject fluctuation in MASC-C scores during the follow-up period (ANXIETY). The primary predictors of ANXIETY were 5-HTTLPR and fluctuations in ALEQ scores during the follow-up period (STRESS). As STRESS was a within-subject predictor, ALEQ scores were centered at each individual’s mean prior to analyses, such that STRESS reflected the upward or downward fluctuation in an individual’s level of stress compared with his/her mean level of stress. 5-HTTLPR was treated as a three-classification variable by triallelic genotyping (LL = −1, SL = 0, SS = 1). To control for individual differences in age, gender, and baseline anxiety symptoms, these variables were included in this model. CES-D scores obtained at timepoint 1 were also included in the model to control for depressive symptoms. The two-level model for subject i at timepoint t was:
Level 1 (within-subject)
ANXIETY ti = β0i + β1i(STRESS) ti + e ti
Level 2 (between-subject)
β0i = γ00 + γ01(timepoint 1 MASC-C) i + γ02(timepoint 1 CES-D) i + γ03(Age) i + γ04(Gender) i + γ05(5-HTTLPR) i + u0i
β1i = γ10 + γ11(5-HTTLPR) i + u1i.
The mixed model:
ANXIETY ti = γ00 + γ10(STRESS) ti + γ05(5-HTTLPR) i + γ01(timepoint 1 MASC-C) i + γ02(timepoint 1 CES-D) i + γ03(Age) i + γ04(Gender) i + γ11(STRESS) ti (5-HTTLPR) i + [u0i + u1i(STRESS) ti + e ti ]
Frequency of 5-HTTLPR genotypes
The frequency distributions of the 5-HTTLPR genotypes are as following: SS: 351 (173 boys, 178 girls), SL: 252 (130 boys, 122 girls), LL: 48 (25 boys, 23 girls). The genotype frequencies were consistent with Hardy-Weinberg equilibrium (χ2 = 0.034, p >0.05). No gender difference in 5-HTTLPR frequency distributions was observed.
Descriptive anxiety, depression and stress data
One sample Kolmogorov-Smirnov tests showed that MASC-C scores in 4 measurements all accorded with normal distribution (all p >0.05). The Pearson’s correlation coefficients between CES-D scores and MASC-C scores in 4 measurements ranged from 0.36 to 0.58 (all p <0.001). The means and standard deviations of all assessments and their gender differences are presented in Table 1. Overall, MASC-C and ALEQ scores decreased. Females reported higher levels of anxiety symptoms than did males at all assessments. Gender differences in MASC-C scores were significant at the initial assessment (t = −2.34, p <0.05) and at 6 months (t = −2.24, p <0.05). Males reported more SLEs than did females at 9 months (t = 2.87, p <0.05). Differences of scores in all assessments between 5-HTTLRP genotypes are presented in Table 2. There were significant differences between 3 genotypes on ALEQ scores at 9 months (F = 3.41, p <0.05).
Statistical analyses of interaction between 5-HTTLPR and stress
To select a covariance structure for our analyses, we fitted the models utilizing each structure and chose the best fit based on Akaike information criterion (AIC and AICC) and Schwarz Bayesian criterion (BIC). The best fit was a heterogeneous autoregressive structure (ARH). With respect to random effects, the ARH parameter (p <0.001), random slope (p <0.01) and random intercept (p <0.001) were all retained in the model. Table 3 showed the estimates of covariance parameter for the final model. Preliminary analyses indicated no gender differences in models with 5-HTTLPR × stress interaction, and thus fixed-effects component of the model analyses are presented for the sample as a whole. Analyses of HLM showed a significant main effect of stress on anxiety symptoms (B = 0.21, p <0.001, Table 4). No significant main effect of 5-HTTLPR on anxiety symptoms was found (B = 0.32, p >0.05). After controlling for age, gender, initial anxiety and depressive symptoms, a significant two-way, cross-level interaction between 5-HTTLPR and stress were detected (B = −0.08, p <0.01). As stress levels increase, the anxiety level among SS carriers appears to increase at a slower rate compared to LL carriers. To present the form of this interaction, the model summarized in Table 4 was used to calculate predicted anxiety symptom scores for participants with SS, SL and LL genotypes who experienced a low or high level of stressful life events (plus or minus 1.5 × mean within-subject standard deviation) without controlling for age, gender or depressive symptoms (see Fig. 1).
As found in many other studies [26–28], females reported higher levels of anxiety symptoms than did males at each assessment in the current study, indicating that females are more likely to experience anxiety symptoms in adolescence. Stress and anxiety symptom levels were highest at the initial assessment and decreased during the follow-up period. A possible explanation for this phenomenon is that the students had just entered senior high school when the first assessment was carried out. This major transition may have involved increased competition and academic pressure, elevating stress and anxiety symptoms. The decreases in stress and anxiety levels may reflect acclimation to senior high school.
No evidence of a 5-HTTLPR main effect was found in the present study, in accordance with the findings of previous researches on the association between 5-HTTLPR and anxiety symptoms , generalized anxiety disorder  or anxiety sensitivity , suggesting that 5-HTTLPR does not directly affect anxiety symptoms.
Regarding gene-environment interaction, this is the first longitudinal study to clarify the interaction between 5-HTTLPR and stress in the prediction of anxiety symptoms in Chinese populations. In consideration of the limitation of sample size in our previous depression study , along with high correlation between depression and anxiety in adolescents , we performed the current study in a larger sample and controlled for the effects of depressive symptoms on fluctuation of anxiety symptoms to improve reliability of the results. A two-way interaction of 5-HTTLPR and SLEs was detected in the current study, LL genotypes increase genetic vulnerability to the experience of anxiety symptoms slightly in response to daily stressors. Consistent with our finding, an interaction effect of the 5-HTTLPR LL genotypes and family adversity , traumatic events , or childhood maltreatment , respectively, has been observed for anxiety disorders or anxiety sensitivity. In contrast, some studies indicated significant interactions between S allele and SLEs or childhood maltreatment on anxiety-related phenotypes [9, 10, 31]. Other studies even suggested no interaction of 5-HTTLPR × stress on anxiety [12, 13].
Several reasons may account for inconsistencies in 5-HTTLPR-related association findings in anxiety. Firstly, ethnic background may be one of the major reasons for the inconsistent results. For instance, most 5-HTTLPR × E studies in depression point to a general tendency for detrimental effect of the 5-HTTLPR S allele. But one of the few study in Chinese population indicated that individuals carrying L alleles could be susceptible to major depression when exposed to negative life events . There were essential uncertainties in ethnic differences of 5-HTTLPR allele frequency, 5-HTT availabilities, 5-HTT uptake and the central serotonergic activity, which might be responsible for the inconsistency. Frequency of the L allele is much lower than that of S allele in Chinese population  and Japan population  but higher in Caucasians or African-Americans [4, 34]. Meta-analyses have demonstrated the different effects of 5-HTTLPR on selective serotonin reuptake inhibitors (SSRIs) efficacy between Caucasians and Asians [35, 36]. Moreover, functional magnetic resonance imaging (fMRI) studies have shown a link between S allele and higher amygdala activation in response to emotional stimuli in Caucasians [37–39], whereas in Asians L allele was associated with amygdala hyperactivation . Therefore, these fascinating questions of different 5-HTTLPR effects between Caucasian and Asian populations require further investigation. It was notable that results of the current study were inconsistent with our previous depression study which indicated risk effect of S allele in girls . The reasons were various (e.g. sample size, other gene variants, and different pathological mechanisms of 5-HTTLPR in anxiety and depression) and should be explored in future.
Second, sample age spanning may be another possible reason contributing to inconsistencies. The majority of studies using subjects whose ages spanned a wide range might neglect developmental issues. Given the dynamics of genetic influences across the lifespan, the impact of genetic factors is likely to depend on developmental stages . Indeed, findings on the neurobiological function of the 5-HTTLPR were not stable between adolescents and adults. A fMRI study assessing amygdala function demonstrated that adolescents with current anxiety or major depressive disorder who were carries of 5-HTTLPR LL genotype exhibited higher amygdale activation to fear faces . This finding was contrary to those reported from affected adults, indicating greater amygdala response in S allele carries .
Furthermore, other genetic variants are variables to be considered in light of the inconsistent findings. It is conceivable that a single gene variation and a certain environmental factor cannot completely reflect the complexity of gene-environment interactions. Thus, the exploration of multiple genetic variations interacting with multiple environmental factors may be needed. Regarding further genetic factors, there are already studies reporting interactions between 5-HTTLPR, other polymorphisms and stress on depression [17, 44]. Therefore, future research should consider gene-gene-environment interactions to gain a more comprehensive understanding of the pathological mechanism of anxiety.
Strengths of the present study included its longitudinal design, which allowed us to control for pre-existing anxiety and depression symptoms to ascertain the effects of stress and 5-HTTLPR on fluctuation of anxiety symptoms in follow-up intervals. We examined 5-HTTLPR, in conjunction with within-subject fluctuations in stress level, to predict within-subject fluctuations in the level of anxiety symptoms. This idiographic approach and relatively reliable estimate of each adolescent’s degree of stress reaction minimized the effects of individual differences in variables. The study was also strengthened by a sample of homogenous population that reduced biases.
However, several limitations of the current study should be noticed. First, self-report measures were used to assess SLEs, which may be less accurate than the use of contextual stress interviews. Further research would benefit from using contextual stress interviews in a multiwave framework. Second, the present study examined the G × E interaction only in adolescents aged 14–17 years. Additional research is needed to test whether findings from this study can be generalized to younger adolescents or children. Third, the sample of this study was non-clinical sample which limited the generalizability of the findings to patients with anxiety disorders. Additionally, the effect size of 5-HTTLPR × stress interaction was modest. It is important to replicate these results in future so that we can be more confident in this interaction. Finally, a single nucleotide polymorphism (SNP) rs25531 within L allele has been described . Due to this SNP, the L allele can be further categorized into LA and LG allele, with LG allele functionally equivalent to S allele . In the current study, genotype analysis did not differentiate between LA and LG alleles. However, considering the frequency of L allele frequency in the present study (0.267) and the minor allele frequency of rs25531 [G = 0.138, according to the NCBI Variation Database (http://www.ncbi.nlm.nih.gov/SNP/)], the LG allele should occur infrequently in our study sample. Moreover, present findings remain unclear in studies reclassified LG allele with lower transcription efficiency. And a study by Martin et al.  questioned the functional interpretation of the LG allele.
In conclusion, results of the current study indicate an interaction between 5-HTTLPR and stress in the prediction of anxiety symptoms among Chinese adolescents. Individuals with the L allele exhibited moderately heightened anxiety response to stress. If the results can be replicated, the current study will provide new evidence for exploring the roles of genetic and environmental factors in the pathological mechanism of anxiety.
Hoffman DL, Dukes EM, Wittchen HU. Human and economic burden of generalized anxiety disorder. Depress Anxiety. 2008;25:72–90.
McLaughlin KA, Conron KJ, Koenen KC, Gilman SE. Childhood adversity, adult stressful life events, and risk of past-year psychiatric disorder: a test of the stress sensitization hypothesis in a population-based sample of adults. Psychol Med. 2010;40:1647–58.
Heils A, Teufel A, Petri S, Stober G, Riederer P, Bengel D, et al. Allelic variation of human serotonin transporter gene expression. J Neurochem. 1996;66:2621–4.
Lesch KP, Bengel D, Heils A, Sabol SZ, Greenberg BD, Petri S, et al. Association of anxiety-related traits with a polymorphism in the serotonin transporter gene regulatory region. Science. 1996;274:1527–31.
Osher Y, Hamer D, Benjamin J. Association and linkage of anxiety-related traits with a functional polymorphism of the serotonin transporter gene regulatory region in Israeli sibling pairs. Mol Psychiatry. 2000;5:216–9.
Middeldorp CM, de Geus EJC, Beem AL, Lakenberg N, Hottenga JJ, Slagboom PE, et al. Family based association analyses between the serotonin transporter gene polymorphism (5-HTTLPR) and neuroticism, anxiety and depression. Behav Genet. 2007;37:294–301.
Narasimhan S, Hodge R, Doyle GA, Kraemer DJM, Prabhakaran R, Rickels K, et al. Association analysis between the 5-HTTLPR polymorphism in the SLC6A4 gene and generalized anxiety disorder. Psychiatr Genet. 2011;21:267–8.
Karg K, Burmeister M, Shedden K, Sen S. The serotonin transporter promoter variant (5-HTTLPR), stress, and depression meta-analysis revisited: evidence of genetic moderation. Arch Gen Psychiatry. 2011;68:444–54.
Gunthert KC, Conner TS, Armeli S, Tennen H, Covault J, Kranzler HR. Serotonin transporter gene polymorphism (5-HTTLPR) and anxiety reactivity in daily life: a daily process approach to gene-environment interaction. Psychosom Med. 2007;69:762–8.
Petersen IT, Bates JE, Goodnight JA, Dodge KA, Lansford JE, Pettit GS, et al. Interaction between serotonin transporter polymorphism (5-HTTLPR) and stressful life events in adolescents’ trajectories of anxious/depressed symptoms. Dev Psychol. 2012;48:1463–75.
Grabe HJ, Spitzer C, Schwahn C, Marcinek A, Frahnow A, Barnow S, et al. Serotonin transporter gene (SLC6A4) promoter polymorphisms and the susceptibility to posttraumatic stress disorder in the general population. Am J Psychiatry. 2009;166:926–33.
Middeldorp CM, de Geus EJC, Willemsen G, Hottenga JJ, Slagboom PE, Boomsma DI. The serotonin transporter gene length polymorphism (5-HTTLPR) and life events: no evidence for an interaction effect on neuroticism and anxious depressive symptoms. Twin Res Hum Genet. 2010;13:544–9.
Zavos HM, Wong CC, Barclay NL, Keers R, Mill J, Rijsdijk FV, et al. Anxiety sensitivity in adolescence and young adulthood: the role of stressful life events, 5HTTLPR and their interaction. Depress Anxiety. 2012;29:400–8.
Irons C, Gilbert P. Evolved mechanisms in adolescent anxiety and depression symptoms: the role of the attachment and social rank systems. J Adolesc. 2005;28:325–41.
Axelson DA, Birmaher B. Relation between anxiety and depressive disorders in childhood and adolescence. Depress Anxiety. 2001;14:67–78.
Zalsman G, Huang YY, Oquendo MA, Burke AK, Hu XZ, Brent DA, et al. Association of a triallelic serotonin transporter gene promoter region (5-HTTLPR) polymorphism with stressful life events and severity of depression. Am J Psychiatry. 2006;163:1588–93.
Zhang K, Xu Q, Xu Y, Yang H, Luo J, Sun Y, et al. The combined effects of the 5-HTTLPR and 5-HTR1A genes modulates the relationship between negative life events and major depressive disorder in a Chinese population. J Affect Disord. 2009;114:224–31.
Ming QS, Zhang Y, Chai QL, Chen HY, Hou CJ, Wang MC, et al. Interaction between a serotonin transporter gene promoter region polymorphism and stress predicts depressive symptoms in Chinese adolescents: a multi-wave longitudinal study. BMC Psychiatry. 2013;13:142.
Silverman WK, Albano AM. The Anxiety Disorders Interview Schedule for Children for DSM-IV: child and parent versions. San Antonio: Psychological Corporation-Harcourt, Brace, & Company; 1996.
Kaufman J, Birmaher B, Brent D, Rao U, Flynn C, Moreci P, et al. Schedule for affective disorders and schizophrenia for school-age children-present and lifetime version (K-SADS-PL): initial reliability and validity data. J Am Acad Child Adolesc Psychiatry. 1997;36:980–8.
Yao S, Zou T, Zhu X, Abela JRZ, Auerbach RP, Tong X. Reliability and validity of the Chinese version of the multidimensional anxiety scale for children among Chinese secondary school students. Child Psychiat Hum D. 2007;38:1–16.
Radloff LS. The CES-D scale: a self-report depression scale for research in the general population. Appl Psychol Meas. 1997;1:385–401.
Hankin BL, Abramson LY. Measuring cognitive vulnerability to depression in adolescence: reliability, validity, and gender differences. J Clin Child Psychol. 2002;31:491–504.
Wang M, Armour C, Wu Y, Ren F, Zhu X, Yao S. Factor structure of the CES-D and measurement invariance across gender in Mainland Chinese adolescents. J Clin Psychol. 2013;69:966–79.
Abela JRZ, Stolow D, Mineka S, Yao S, Zhu X, Hankin BL. Cognitive vulnerability to depressive symptoms in adolescents in urban and rural Hunan, China: a multiwave longitudinal study. J Abnorm Psychol. 2011;120:765–78.
Casper RC, Belanoff J, Offer D. Gender differences, but no racial group differences, in self-reported psychiatric symptoms in adolescents. J Am Acad Child Adolesc Psychiatry. 1996;35:500–8.
Lewinsohn PM, Gotlib IH, Lewinsohn M, Seeley JR, Allen NB. Gender differences in anxiety disorders and anxiety symptoms in adolescents. J Abnorm Psychol. 1998;107:109–17.
Angelopoulos N, Economou M. Prevalence of anxiety and depressive symptoms in a high-school students population. Eur Psychiat. 1994;9:19–26.
Laucht M, Treutlein J, Blomeyer D, Buchmann AF, Schmid B, Becker K, et al. Interaction between the 5-HTTLPR serotonin transporter polymorphism and environmental adversity for mood and anxiety psychopathology: evidence from a high-risk community sample of young adults. Int J Neuropsychoph. 2009;12:737–47.
Klauke B, Deckert J, Reif A, Pauli P, Zwanzger P, Baumann C, et al. Serotonin transporter gene and childhood trauma - a G*E effect on anxiety sensitivity. Depress Anxiety. 2011;28:1048–57.
Stein MB, Schork NJ, Gelernter J. Gene-by-environment (serotonin transporter and childhood maltreatment) interaction for anxiety sensitivity, an intermediate phenotype for anxiety disorders. Neuropsychopharmacology. 2008;33:312–9.
He Q, Xue G, Chen C, Lu Z, Dong Q, Lei X, et al. Serotonin transporter gene-linked polymorphic region (5-HTTLPR) influences decision making under ambiguity and risk in a large Chinese sample. Neuropharmacology. 2010;59:518–26.
Kunugi H, Hattori M, Kato T, Tatsumi M, Sakai T, Sasaki T, et al. Serotonin transporter gene polymorphisms: ethnic difference and possible association with bipolar affective disorder. Mol Psychiatry. 1997;2:457–62.
Gelernter J, Kranzler H, Coccaro EF, Siever LJ, New AS. Serotonin transporter protein gene polymorphism and personality measures in African American and European American subjects. Am J Psychiatry. 1998;155:1332–8.
Serretti A, Kato M, De Ronchi D, Kinoshita T. Meta-analysis of serotonin transporter gene promoter polymorphism (5-HTTLPR) association with selective serotonin reuptake inhibitor efficacy in depressed patients. Mol Psychiatry. 2007;12:247–57.
Porcelli S, Fabbri C, Serretti A. Meta-analysis of serotonin transporter gene promoter polymorphism (5-HTTLPR) association with antidepressant efficacy. Eur Neuropsychopharmacol. 2012;22:239–58.
Hariri AR, Mattay VS, Tessitore A, Kolachana B, Fera F, Goldman D, et al. Serotonin transporter genetic variation and the response of the human amygdala. Science. 2002;297:400–3.
Hariri AR, Drabant EM, Munoz KE, Kolachana BS, Mattay VS, Egan MF, et al. A susceptibility gene for affective disorders and the response of the human amygdala. Arch Gen Psychiatry. 2005;62:146–52.
Murphy SE, Norbury R, Godlewska BR, Cowen PJ, Mannie ZM, Harmer CJ, et al. The effect of the serotonin transporter polymorphism (5-HTTLPR) on amygdala function: a meta-analysis. Mol Psychiatry. 2013;18:512–20.
Li S, Zou Q, Li J, Wang D, Yan C, Dong Q, et al. 5-HTTLPR polymorphism impacts task-evoked and resting-state activities of the amygdala in Han Chinese. PLoS One. 2012;7:e36513.
Reiss D, Neiderhiser JM. The interplay of genetic influences and social processes in developmental theory: specific mechanisms are coming into view. Dev Psychopathol. 2000;12:357–74.
Lau JY, Goldman D, Buzas B, Fromm SJ, Guyer AE, Hodgkinson C, et al. Amygdala function and 5-HTT gene variants in adolescent anxiety and major depressive disorder. Biol Psychiatry. 2009;65:349–55.
Munafo MR, Brown SM, Hariri AR. Serotonin transporter (5-HTTLPR) genotype and amygdala activation: a meta-analysis. Biol Psychiatry. 2008;63:852–7.
Priess-Groben HA, Hyde JS. 5-HTTLPR X stress in adolescent depression: moderation by MAOA and gender. J Abnorm Child Psychol. 2013;41:281–94.
Hu XZ, Lipsky RH, Zhu GS, Akhtar LA, Taubman J, Greenberg BD, et al. Serotonin transporter promoter gain-of-function genotypes are linked to obsessive-compulsive disorder. Am J Hum Genet. 2006;78:815–26.
Martin J, Cleak J, Willis-Owen SA, Flint J, Shifman S. Mapping regulatory variants for the serotonin transporter gene based on allelic expression imbalance. Mol Psychiatry. 2007;12:421–2.
Funding of this study was provided by grants from the National Natural Science Foundation of China (no.81471384) and the construct program of the key discipline in Hunan Province. The funder had no role in the study design, data collection and analysis, decision to publish, or preparation of the manuscript. We also thank all participants for their participation and support.
The authors declare that they have no competing interests.
All authors have contributed significantly to this manuscript and agree with its content. Author YS, ZX and YJ designed the study and wrote the protocol. Author MQ, ZY, YJ and WX took part in data collecting. Author MQ and YJ managed statistical analysis. Author MQ completed the literature searches and the first draft of the manuscript. All authors read and approved the final manuscript.