Anxiety-related Endophenotypes and Hazardous Alcohol Use in Young Adults are Associated with a Functional Polymorphism in the SLC6A4 Gene

Karen M. Jiménez1, Angela J. Pereira-Morales1, Ana Adan2, 3, Sandra Lopez-Leon4, Diego A. Forero1, 5, 6, *
1 Laboratory of NeuroPsychiatric Genetics, Biomedical Sciences Research Group, School of Medicine, Universidad Antonio Nariño, Bogotá, Colombia
2 Department of Clinical Psychology and Psychobiology, University of Barcelona, Barcelona, Spain
3 Institute of Neurosciences, University of Barcelona, Barcelona, Spain
4 Novartis Pharmaceuticals Corporation, One Health Plaza, East Hanover, NJ, 07936-1080, USA
5 School of Medicine, Universidad Antonio Nariño, Bogotá, Colombia
6 Laboratory of NeuroPsychiatric Genetics, School of Medicine, Universidad Antonio Nariño, 110231, Bogotá, Colombia

Article Metrics

CrossRef Citations:
Total Statistics:

Full-Text HTML Views: 145
Abstract HTML Views: 118
PDF Downloads: 46
ePub Downloads: 36
Total Views/Downloads: 345
Unique Statistics:

Full-Text HTML Views: 98
Abstract HTML Views: 74
PDF Downloads: 32
ePub Downloads: 26
Total Views/Downloads: 230

© 2019 Jiménez et al.

open-access license: This is an open access article distributed under the terms of the Creative Commons Attribution 4.0 International Public License (CC-BY 4.0), a copy of which is available at: This license permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

* Address correspondence to this author at the Laboratory of NeuroPsychiatric Genetics, School of Medicine, Universidad Antonio Nariño, 110231, Bogotá, Colombia; Tel: +57 313 610427; Email:



A functional polymorphism (5-HTTLPR, rs4795541) in the serotonin transporter (SLC6A4) gene has been shown as an important candidate for several psychiatric and behavioral traits.


The objective of this study was to examine the possible interaction of this polymorphism with physical neglect in childhood on the presentation of anxiety traits and hazardous alcohol consumption in young Colombian subjects.


272 young adults (mean age: 21.3 years) were evaluated with the Childhood Trauma Questionnaire, the Zung Self-rating Anxiety Scale, the Big Five Inventory, the Cohen’s Perceived Stress Scale, the Alcohol, Smoking, Substance Involvement Screening Test and the Alcohol Use Disorders Identification Test. Genotyping for the 5-HTTLPR polymorphism was carried out using conventional PCR. A linear regression model, corrected by age and gender, was used.


We found that individuals with the L/L genotype showed higher scores on physical neglect (p=0.0047), anxiety symptoms (p=0.028), neuroticism (p=0.019) and perceived stress (p=0.035). L/L genotype was a risk factor for hazardous alcohol use in young adults (OR=3.06, p=0.0003). No GxE interactions were observed in our data.


Our results provide novel evidence for the role of a functional polymorphism in the SLC6A4 gene on the relationship of childhood trauma, anxiety-related traits and risky consumption of alcohol.

Keywords: Alcohol use, Anxiety, Candidate genes, Childhood trauma, Neuroticism, Psychiatric genetics.


Anxiety Disorders (AD), which include panic disorder, generalized anxiety disorder, post-traumatic stress disorder, phobias, and separation anxiety disorder, are among the most common mental disorders around the world [1, 2]. These disorders usually have an early onset followed by a chronic course, being the sixth leading cause of disability worldwide [3]. The current global prevalence of anxiety disorders is 7.3%, suggesting that one in 14 people has an anxiety disorder at any given time [1]. The prevalence is higher in women than in men, with women to men ratio of 1.9:1 and it has been remained relatively constant across time [4]. Prevalence estimates for AD were significantly higher in people exposed to an armed conflict [1] and, in this context, the Colombian population has suffered from one of the longest internal armed conflicts in the world [5]. According to the DSM-V, substance use disorders are characterized by at least two of the diagnostic criteria (such as social and legal problems related to use, among others) [6] and are important contributors to the global burden of disease [7]. Hazardous alcohol use has been defined as a pattern of consumption that leads to a higher risk of harmful effects for the individual [8].

AD and hazardous alcohol use are complex traits with a complex etiology. They both have an environmental and a genetic component. It is expected that the combination of several genetic polymorphisms interacts with environmental variables to increase the risk [9-11]. Even though several genes are associated with both disorders, the SLC6A4 gene, which encodes the serotonin transporter, became a major candidate gene for AD and hazardous alcohol use [10, 12]. It is interesting in the context that several drugs of psychiatric use, such as paroxetine and citalopram, are inhibitors of this transporter [12].

The most commonly studied polymorphism in this gene is the non-coding variant 5-HTTLPR (rs4795541), which is located in the promoter region [13]. This polymorphism is a Variable Number Tandem Repeat (VNTR), composed of 24 bp repeats which generate different alleles. The two main alleles are the short (S) with 14 repeats and the long (L) with 16 repeats [13].

Functional studies have reported that the L allele is associated with higher levels of transcriptional activity and with a higher rate of serotonin uptake [14]. In addition, several works have shown that the S allele, in general, is associated with a less favorable response to pharmacological treatment in affective and anxiety disorders [15, 16].

Lesch et al conducted the first study that found an association between the SLC6A4 promoter polymorphism and anxiety-related traits, including the neurotic personality trait [17], and the authors reported that the 5-HTTLPR poly-morphism accounted for 3-4% of population variability in neuroticism [17]. In other studies, it has been reported that individuals with the S/S genotype scored higher for neuroticism on the big five personality model. In addition, in patients with depression and anxiety disorders it has been observed a higher frequency of the S allele [18, 19]. The 5-HTTLPR polymorphism has been investigated in the context of stressful life events, including childhood maltreatment, showing association with the increase in the susceptibility to depressive and anxious symptoms [20]. Subjects with both the S/S genotype and with high levels of childhood maltreatment have been shown to have a higher sensitivity to anxiety [20].

There is a growing evidence for the association of genes involved in the serotonergic system with hazardous alcohol use [10]. A recent meta-analysis showed that the S allele is associated with an increased risk for alcohol dependence [21]. In addition, other studies have found that subjects with the S allele reported more frequent drinking [22-24]. In contrast, other works have found an association with hazardous drinking in subjects with the L allele [25].

The purpose of this study was to examine the hypothesis that the 5-HTTLPR polymorphism might interact with physical negligence in childhood on the risk for the presence of anxiety symptoms and hazardous alcohol consumption in young adulthood.


2.1. Participants

Two hundred seventy-two Colombian young adults were included in this study (21.3 years of age, SD=3.8, 75% women) [26, 27]. Participants were enrolled in medical and nursing schools and were invited to participate in the study [28, 29]. All participants signed a written informed consent and this study was approved by the Institutional Ethics Committee of the Antonio Nariño University. The procedures performed were in accordance with the ethical standards of the insti-tutional research committee and with the 1964 Helsinki declaration and its later amendments.

A self-report instrument including general demographic and medical information was applied and it served to identify subjects with personal history of neuropsychiatric disorders. Ages over 18 years and absence of personal history of neuropsychiatric disorders were the inclusion criteria for the study (12 subjects were excluded due to it). All the phenotypic data were collected through self-report instruments, which were applied by a psychologist. The included individuals were mainly from low and medium socioeconomic status (34 and 45%, respectively), had mainly secondary as their education level (79%) and single was the most frequent marital status (93%).

2.2. Psychological Measurements

To evaluate physical negligence, the Childhood Trauma Questionnaire (CTQ)-brief version was used [30]. This questionnaire is a 28-item self-administered inventory that provides an assessment of child abuse and neglect through five dimensions: physical, sexual and emotional abuse and physical and emotional neglect. Scores of 9 (for both men and women) are indicative of physical abuse and physical neglect and of 5 and 7 for sexual abuse (in men and women, respectively) [31]. The CTQ has adequate psychometric properties in different languages, included the Spanish [31]. In the current study, the Cronbach´s alpha was 0.89.

Anxiety-related endophenotypes were measured using the Big Five Inventory (BFI-S) [32], in the short version (15 items), the Zung Self-Rating Anxiety Scale (ZSAS) [33] and the Cohen’s Perceived Stress Scale (PSS) [34].

Fig. (1). Overview of the associations found in the current study.

The BFI-S inventory measures 5 dimensions of personality: Factor O (Openness), Factor C (Conscien-tiousness), Factor E (Extraversion), Factor A (Agreeableness) and Factor N (Neuroticism) on a 7-point Likert scale: from 1 (strongly disagree) to 7 (strongly agree) [32]. In our analysis, we only included the factor N (Neuroticism) because it is one of the most established risk factors for anxiety and mood disorders [35]. The Cronbach’s alpha was 0.62 for Neuroticism in the current study.

ZSAS is a self-report questionnaire designed to measure anxiety symptomatology in the general population. It is composed of 20 items with 4 possible responses: (1) never, (2) rarely/sometimes, (3) frequently and (4) always [33]. Each item is scored from 1 to 4 according to the severity of anxiety symptoms and total scores range between 20 to 80 points. ZSAS has shown to have a good internal consistency in Colombia [36]. In the current study the Cronbach’s alpha was 0.85.

The PSS-10 items [34] is a self-report scale used to assess an individual’s experience of general perceived stress over the past month. The scores of this scale range from 0 to 40, higher scores mean higher levels of stress. It has been validated in Colombia population [37]. In the current study, the Cronbach’s alpha was 0.85.

Hazardous alcohol use was measured with two different scales, searching for a better characterization of the con-sumption:

The Alcohol, Smoking and Substance Involvement Screening Test (ASSIST) [38] is an instrument created by the World Health Organization (WHO) for the detection of risky or abusive consumption of substances. This instrument assesses the frequency, patterns of substance use and risky use of tobacco, alcohol, amphetamine-type substances, cannabis, cocaine, hallucinogens, inhalants, opioids, sedatives and other substances through eight questions. Total scores for each substance are categorized as ‘low-risk use’ (0-3), ‘moderate-risk use’ (4-26) or ‘high-risk use’ (>26) [38]; it has been validated in different languages, included Spanish [39]. In this study, alcohol was the substance with higher consumption (42.3%) and for that reason we selected hazardous alcohol use as one of the outcome variables and their Cronbach’s alpha was of 0.75.

The Alcohol use disorders identification test (AUDIT) [8] is a brief written screening method to identify current harmful and hazardous drinking in primary care settings, in both developing and developed countries. The 10-item scale was designed to assess three conceptual domains: drinking, behavior/dependence and alcohol-related problems. It has been validated in a Colombian sample [40]. The cut-off point used in the current study was ≥8 for risk of hazardous alcohol use, based on the recommendations by the WHO [8]. In the current sample, the Cronbach´s alpha was 0.78.

2.3. Genotyping

The extraction of DNA was performed by a salting out method, using 400 µl of peripheral blood [41]. For the genotyping of the 5-HTTLPR polymorphism (rs4795541), a conventional PCR was carried out in a Labnet MultiGene 96-well thermal cycler (Labnet International Inc, Edison, NJ, USA). The PCR reaction was performed in a total volume of 20 µl, with 20 ng of DNA, 0.8 U of Taq polymerase (Bioline, London, United Kingdom), 10 µM of each primer (F: 5´ATG CCA GCA CCT AAC CCC TAA TGT 3´, R: 5´GGA CCG CAA GGT GGG CGG GA 3´), 0.25 mM of each dNTP, 10X reaction buffer and 1.5 mM MgCl2. The PCR program consisted in one cycle of 94 oC for 30 s, 35 cycles of 94oC for 30 s, 66 oC for 30 s, 72oC for 40 s, and a final step of 1 cycle of 72 oC for 5 min. The PCR products were allowed to run for 40 min at 140 V, in a 2% agarose gel. The genotypes were determined as follows: L/L: 419/419 pb, L/S: 419/375 bp, S/S: 375/375 bp. With the aim of verifying the genetic results, 10% of the samples were randomly selected and re-analyzed.

All mandatory health and safety procedures were accomp-lished in all the experiments carried out in this study.

2.4. Statistical Analysis

The psychometric analysis of the scales and questionnaires was performed using standard approaches: Internal consistency through Cronbach’s alpha and factorial structure of the tests using exploratory factorial analysis.

Normal distributions from the psychological scales scores were examined through skewness and kurtosis and bivariate correlations with Pearson Coefficient were used. These statistical analyses were performed with the Statistical Package for the Social Sciences (SPSS v. 18)

A linear regression model corrected by age and gender, using SNPStats program, was used to identify the possible association of the SNP with anxiety symptoms, perceived stress, physical neglect, hazardous alcohol use and neuroticism [42]. SNPStats program was also used to calculate genotype and allelic frequencies and Hardy-Weinberg equilibrium. The best genetic model was selected based on the Akaike Information Criterion [42]. For the current study, the Bonferroni corrected p value (0.05 divided by 6), was 0.008.

In addition, multiple regression models were carried out to examine GxE Genotype (G) and Environment (E) interactions:

(A) Physical neglect in childhood x 5-HTTLPR genotypes to predict anxiety symptoms

(B) Physical neglect x 5-HTTLPR genotypes to predict neurotic trait

(C) Physical neglect x 5-HTTLPR genotypes to predict perceived stress

(D) Physical neglect x 5-HTTLPR genotypes to predict hazardous alcohol use, for both AUDIT and ASSIST alcohol subscale.

To know whether 5-HTTLPR genotype played a role as mediator, the mediation-in-serial models were employed, using multiple regressions, following the procedures recommended by Hayes [43]. The mediation effect occurs if the association between the Independent Variable (IV) and the Outcome Variable (DV) is eliminated (complete mediation) or reduced (partial mediation) when the Mediator Variable (MV) is introduced in the analysis. The following mediation models were carried out:

  1. Physical neglect in childhood was inserted as an IV, anxiety-related phenotypes were inserted as DVs separately and 5-HTTLPR genotype was inserted as the MV.
  2. Physical neglect was inserted as an IV, alcohol use was inserted as DV for AUDIT and ASSIST alcohol subscale, separately and 5-HTTLPR genotype was inserted as the MV.

Direct, indirect and total effects, standard errors and 95% confidence intervals (biased-corrected bootstrapped; 1000= repetitions) for each model, were calculated using speci-fication model 4 of the PROCESS plug in (V.2014) [43] in SPSS (v. 18).


In our sample, the factorial structures found were the expected: six-factor solution for CTQ (64.1% of variance), five-factor solution for BFI-S and ZSAS (62.1% and 56% of variance respectively), two-factor solution for PSS (56% of variance), one-factor solution (46.3% of variance) for alcohol use (ASSIST) and three-factor solution for AUDIT (59% of variance).

The allele frequencies in our sample were 0.51 and 0.49 for the S and L alleles, respectively. The S/S genotype was observed in 26% (n = 70), the S/L genotype in 50% (n = 136) and the L/L genotype in 24% (n = 65) of the subjects in the current sample. The genotype frequencies were in Hardy-Weinberg equilibrium (p = 1).

Correlations were significant only between physical neglect and anxiety symptoms (r = 0.278, p<0.01) and physical neglect and perceived stress (r = 0.221, p<0.01), no association was found between physical neglect and neurotic trait (r = 0.113, p>0.05), or for physical neglect and ASSIST alcohol subscale (r = 0.100, p>0.05) or AUDIT scores (r = 0.103, p>0.05).

We found significant associations between behavioral variables and 5-HTTLPR: anxiety symptoms (p = 0.029), physical neglect (p = 0.0047), alcohol use (p = 0.011), perceived stress (p = 0.035) and neuroticism (p = 0.019) (Table 1). The results for association with physical neglect were significant after a Bonferroni correction for multiple testing.

The risk of alcohol use based on recommended cut-off points in AUDIT (81 subjects for moderate risk) and ASSIST (150 subjects for moderate risk) was associated significantly with 5-HTTLPR polymorphism. The odds ratio was higher for carriers of the long variant of the 5-HTTLPR polymorphism (AUDIT, OR= 2.01, CI: 1.04 – 3.87; ASSIST alcohol subscale, OR= 3.06, 1.62 – 5.79) (Table 2). Results for the ASSIST alcohol subscale were significant after a Bonferroni correction for multiple testing. An overview of the associations found in our work is shown in Fig. (1).

No GxE interactions were observed in the data, for the multiple regression models, as well as for neither of the models of mediation analysis (Table 3). 5-HTTLPR genotype was excluded of all the regression models for anxiety symptoms, neurotic trait and perceived stress (as outcomes variables) and 5-HTTLPR genotypes and physical neglect in childhood (as predictors) (p = 0.467, p = 0.088, p = 0.809 respectively). For hazardous alcohol use models, physical negligence was excluded from both models (p = 0.397 for AUDIT and p = 0.871 for ASSIST alcohol subscale).

Table 1. Association between anxiety symptoms, physical negligence, hazardous alcohol use, neuroticism, perceived stress and 5–HTTLPR polymorphism.
Measurements Genotypes Mean Scores (SE) p value
Anxiety symptoms S/S –S/L 37.8 (0.6) 0.028
L/L 40.5 (1.3)
Physical negligence S/S – S/L 6.3 (0.1) 0.0047*
L/L 7.2 (0.3)
Hazardous alcohol use ASSIST scores S/S – S/L 6.6 (0.5) 0.011
L/L 9.2 (0.9)
AUDIT scores S/S – S/L 4.6 (0.3) 0.25
L/L 5.31 (0.6)
Neuroticism S/S – S/L 4.1 (0.1) 0.019
L/L 4.47 (0.1)
Perceived stress S/S – L/L 18.81 (0.6) 0.035
S/L 17.32 (0.5)
*Significant after Bonferroni corrections; Abbreviations: ASSIST: Alcohol, Smoking and Substance Involvement Screening Test; AUDIT: Alcohol Use Disorders Identification Test.
Table 2. Association between 5 –HTTLPR and risk of hazardous alcohol use.
Scale Genotype Risk use OR (95% CI) p value
Low Moderate
AUDIT L/S – L/L 134 (70.4%) 67 (82.7%) 2.01 (1.04 – 3.87) 0.03
S/S 56 (29.5%) 14 (17.3%) 1.0
ASSIST S/S – L/S 104 (86%) 102 (68%) 1.0
L/L 17 (14.1%) 48 (32%) 3.06 (1.62 – 5.79) 0.00033*
Note: Low and moderate risks are based on cut-off points for both scales. Cut-off point for risk use in AUDIT: 8 men, 6 women. Cut-off point for moderate-risk of consumption in the ASSIST alcohol subscale: 4.; *Significant after the Bonferroni correction.; Abbreviations: ASSIST: Alcohol, Smoking and Substance Involvement Screening Test; AUDIT: Alcohol Use Disorders Identification Test
Table 3. Interaction between physical neglect and 5-HTTLPR genotypes to predict outcome variables.
Multiple Regression Mediation Analysis
Physical neglect in childhood 5-HTTLPR interaction 5-HTTLPR interactiona
B 95% CI p B 95% CI p Effect 95% CI p
Neurotic trait 0.155 0.01-0.11 0.001 0.103 -0.02-0.37 0.088 0.002 -0.002-0.012 >0.05
Anxiety symptoms 0.281 0.52-1.25 0.001 0.043 -0.97-2.11 0.467 0.007 0.009-0.084 >0.05
Perceived stress 0.23 0.26-0.77 0.001 0.014 -0.94-1.20 0.809 0.001 -0.010-0.043 >0.05
Hazardous alcohol use (AUDIT) 0.05 -0.01-0.02 0.397 0.093 -0.16-1.32 0.128 0.007 -0.006-0.040 >0.05
Hazardous alcohol use (ASSIST) 0.01 -0.01-0.02 0.871 0.176 0.04-0.20 0.004 0.006 -0.016-0.079 >0.05
Note: a Indirect effect of physical neglect on neurotic trait, anxiety symptoms, perceived stress and hazardous alcohol use via 5-HTTLPR genotypes.

In the mediation models, 5-HTTLPR genotype was not significant mediator for any relationship between variables (Table 3).


The findings of the present study show an association between physical neglect in childhood, anxiety-related phenotypes, hazardous alcohol use and 5-HTTLPR poly-morphism (Fig. 1). Our results suggest that 5-HTTLPR polymorphism could be an important risk factor for hazardous alcohol consumption and anxiety-related endophenotypes, such as the neurotic trait, anxiety symptoms and perceived stress in young adults. However, we did not find GXE interactions between 5-HTTLPR polymorphism and physical neglect to predict anxiety-related phenotypes or hazardous alcohol use.

Childhood maltreatment is an important risk factor for affective disorders and a major public health problem [44]. Exposure to emotional and physical neglect is associated with increased risk for a wide range of behavioral problems, including anxiety disorders, alcohol abuse, depression and suicidal behavior [45].

Previous findings have shown that many potential mechanisms may explain the mentioned association between neglect in childhood and increased risk for mental health problems and alcohol consumption. One of the most plausible explanations is that the neurobiological development can be altered by the maltreatment during childhood, which in turn can negatively affect a child's cognitive, emotional, and social growth, leading to behavioral and learning problems that can persist throughout life [44]. Moreover, it is widely accepted that in many cases some victims of childhood maltreatment may try to manage the subsequent distress through the use of alcohol, tobacco, medicaments and other drugs [46]. In addition, epigenetic studies have proposed that early life events affect DNA methylation of the SLC6A4 gene [47, 48].

Previous findings have shown that the risk for affective disorders is elevated among carriers of the S allele when they were exposed to stressful life events [49]. However, our data suggest that the L/L genotype is likely to be important on psychosocial development, particularly in those settings that increase the risk of anxiety symptoms, neurotic trait and hazardous alcohol use. These findings are consistent with some studies that revealed an interaction between the L allele and the increased risk of depression in subjects exposed to stress [25]; however, few studies have replicated these findings in healthy young adults with measures of anxiety symptoms [25].

To date, there are few studies that have investigated the association of the functional 5-HTTLPR polymorphism with hazardous alcohol use in healthy subjects. While some studies showed evidence of an association between 5-HTTLPR and psychological variables related to hazardous alcohol use [50], others found a significant result with higher alcohol intake [22-25]. Our findings of a dominant model for the association of the S allele and hazardous alcohol use is supported by previous articles [23, 51]. In addition, our finding of heterosis in the association between HTTLPR and perceived stress is concordant with a previous study [52].

Multiple structural neuroimaging studies have shown that young adults with history of heavy drinking have lower volumes in prefrontal cortex [53], with some of these changes being partially reversible [54]. In addition, several functional neuroimaging studies have shown that these young adults also have elevated levels of brain activity in fronto-parietal areas during cognitive tasks and increased levels of neural responses in several mesocorticolimbic areas, such as the anterior cingulate cortex, the amygdala and the striatum [53]. Animal models have also found that alcohol exposure affects neurochemical processes in mesocorticolimbic brain regions [55].

In accordance with the above findings, the association of the 5-HTTLPR polymorphism with childhood maltreatment and anxiety-related phenotypes is not clearly understood and the reports of association have been inconsistent and difficult to replicate [56]. Our results have further shown a protective role for the S allele in the risk for alcohol consumption and the development of anxiety symptoms and neuroticism. Our findings are supported in part by previous results where S/S genotype of the 5-HTTLPR polymorphism appeared to be protective against depressive symptomatology in a supportive early context [57]. However, several studies have reported the opposite effect, with the S/S genotype being a risk factor for affective disorders [49].

A possible contributory factor of this inconsistency may be heterogeneity in the measurements and characteristics of environmental adversity investigated in the different studies, particularly, the nature of the stressors studied. For example, the duration of exposure to adversity and the type of stressors studied [58]. Recall bias in the psychological assessment might be another cause of inconsistent findings; it is widely accepted that retrospective reports on traits or experiences are often discrepant (in comparison with prospective studies). In many cases, individuals may attempt to mentally reconstruct information about past events, personality trait or common behaviors, by generating detailed reports based on their general beliefs [59]. In patients with anxiety disorders, such as panic disorder and obsessive-compulsive disorder, it has been documented a negative recall bias in retrospective self-reports, resulting in a retrospective overestimation of the disorder-specific symptoms [60].


In the current study, two types of stress were found as associated with the SLC6A4 gene: chronic stress and acute stress, through physical neglect in childhood and perceived stress over the past month respectively; but only results for physical neglect were significant after a Bonferroni correction for multiple testing. This is consistent with the results of a meta-analysis [18] that found that the effect of 5-HTTLPR polymorphism (S allele) on neurotic trait is only apparent when those carrying the short allele are exposed to the stressful environment and with the results of another meta-analysis [58] that reported robust evidence for 5-HTTLPR moderating the relationship between both childhood maltreatment and specific stressors and depressive symptoms.

Current understanding of the specific characteristics of the stressors that are moderated by 5-HTTLPR is insufficient [58]. However, considering previous data and our current results, specific and chronic stressors (such as childhood maltreatment) might be more important than acute stressors in the relationship of 5-HTTLPR polymorphism and anxiety related-traits.

Limitations of the current study include the analysis of a single polymorphism in the SLC6A4 gene and the use of a convenience sample design, which includes a higher proportion of females.

To our knowledge, the current study is the first work that explores the relationship between 5-HTTLPR polymorphism and anxiety-related phenotypes and hazardous alcohol consumption using validated scales for several psychosocial variables (Table S1) [18, 58]. Furthermore, this is the first study of this type in a Latin American sample.


In futures studies, it could be important to examine the relationship between 5-HTTLPR polymorphism and exposure to environmental adversity, using different phenotypic characterizations as well as different types of adversity. In this context, future longitudinal studies, involving functional neuroimaging methods, will provide an interesting approach for analyzing the cumulative effect of genetic factors, early alcohol use and childhood trauma on brain function.


KMJ participated in acquisition and analysis of genetic data and drafting and critical revision of the manuscript. AJP participated in study design, acquisition and analysis of psychological data and drafting and critical revision of the manuscript. SL-L participated in analysis of genetic data and drafting and critical revision of the manuscript. AA participated in analysis of psychological data and drafting and critical revision of the manuscript. D.A.F participated in study design, acquisition and analysis of genetic and psychological data, drafting and critical revision of the manuscript.


The study was approved by the Institutional Ethics Committee of the Antonio Nariño University.


No animals were used in this research. All human research procedures followed were in accordance with the ethical standards of the committee responsible for human experi-mentation (institutional and national), and with the 1964 Helsinki declaration and its later amendments.


Written informed consent has been obtained from all the participants.


The datasets used and/or analyzed during the current study are available from the corresponding author on reasonable request.


This study was supported by a research grant from Colciencias (grant # 823-2015) and by a grant from the Spanish Ministry of Economy, Industry and Competitiveness PSI2015-65026 (MINECO/FEDER/UE).


The authors declare no conflict of interest, financial or otherwise.


The authors thank to Andres Camargo, who assisted with recruitment and evaluations of subjects.


Supplementary material is available on the publishers Website along with the published article.

Download File


[1] Baxter AJ, Scott KM, Vos T, Whiteford HA. Global prevalence of anxiety disorders: A systematic review and meta-regression. Psychol Med 2013; 43(5): 897-910.
[2] Kessler RC, Aguilar-Gaxiola S, Alonso J, et al. The global burden of mental disorders: An update from the WHO World Mental Health (WMH) surveys. Epidemiol Psichiatr Soc 2009; 18(1): 23-33.
[3] Baxter AJ, Vos T, Scott KM, Ferrari AJ, Whiteford HA. The global burden of anxiety disorders in 2010. Psychol Med 2014; 44(11): 2363-74.
[4] Baxter AJ, Scott KM, Ferrari AJ, Norman RE, Vos T, Whiteford HA. Challenging the myth of an “epidemic” of common mental disorders: trends in the global prevalence of anxiety and depression between 1990 and 2010. Depress Anxiety 2014; 31(6): 506-16.
[5] Rodriguez C, Sanchez F. Armed conflict exposure, human capital investments, and child labor: Evidence from Colombia. Defence Peace Econ 2012; 23(2): 161-84.
[6] Bhad R, Lal R, Balhara YP. Disorders related to use of psychoactive substances in DSM-5: Changes and challenges. Indian J Psychol Med 2015; 37(4): 470-2.
[7] Whiteford HA, Degenhardt L, Rehm J, et al. Global burden of disease attributable to mental and substance use disorders: Findings from the Global Burden of Disease Study 2010. Lancet 2013; 382(9904): 1575-86.
[8] Saunders JB, Aasland OG, Babor TF, de la Fuente JR, Grant M. Development of the alcohol use disorders identification test (AUDIT): WHO collaborative project on early detection of persons with harmful alcohol consumption--II. Addiction 1993; 88(6): 791-804.
[9] Forero DA, Vélez-van-Meerbeke A, Deshpande SN, Nicolini H, Perry G. Neuropsychiatric genetics in developing countries: Current challenges. World J Psychiatry 2014; 4(4): 69-71.
[10] Kim J, Park A. A systematic review: Candidate gene and environment interaction on alcohol use and misuse among adolescents and young adults. Am J Addict 2018.
[11] Shimada-Sugimoto M, Otowa T, Hettema JM. Genetics of anxiety disorders: Genetic epidemiological and molecular studies in humans. Psychiatry Clin Neurosci 2015; 69(7): 388-401.
[12] Gelernter J. SLC6A4 polymorphism, population genetics, and psychiatric traits. Hum Genet 2014; 133(4): 459-61.
[13] Iurescia S, Seripa D, Rinaldi M. Role of the 5-HTTLPR and SNP promoter polymorphisms on serotonin transporter gene expression: A closer look at genetic architecture and In Vitro functional studies of common and uncommon allelic variants. Mol Neurobiol 2016; 53(8): 5510-26.
[14] Hu XZ, Lipsky RH, Zhu G, et al. Serotonin transporter promoter gain-of-function genotypes are linked to obsessive-compulsive disorder. Am J Hum Genet 2006; 78(5): 815-26.
[15] Kato M, Serretti A. Review and meta-analysis of antidepressant pharmacogenetic findings in major depressive disorder. Mol Psychiatry 2010; 15(5): 473-500.
[16] Bauer IE, Graham DP, Soares JC, Nielsen DA. Serotonergic gene variation in substance use pharmacotherapy: A systematic review. Pharmacogenomics 2015; 16(11): 1307-14.
[17] Lesch KP, Bengel D, Heils A, et al. Association of anxiety-related traits with a polymorphism in the serotonin transporter gene regulatory region. Science 1996; 274(5292): 1527-31.
[18] Munafò MR, Freimer NB, Ng W, et al. 5-HTTLPR genotype and anxiety-related personality traits: A meta-analysis and new data. Am J Med Genet B Neuropsychiatr Genet 2009; 150B(2): 271-81.
[19] Balestri M, Calati R, Serretti A, De Ronchi D. Genetic modulation of personality traits: A systematic review of the literature. Int Clin Psychopharmacol 2014; 29(1): 1-15.
[20] 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(2): 312-9.
[21] Cao J, Hudziak JJ, Li D. Multi-cultural association of the serotonin transporter gene (SLC6A4) with substance use disorder. Neuropsychopharmacology: Ameri College Neuropsychopharmacol 2013; 38(9): 1737-47.
[22] Covault J, Tennen H, Armeli S, et al. Interactive effects of the serotonin transporter 5-HTTLPR polymorphism and stressful life events on college student drinking and drug use. Biol Psychiatry 2007; 61(5): 609-16.
[23] Daw J, Shanahan M, Harris KM, Smolen A, Haberstick B, Boardman JD. Genetic sensitivity to peer behaviors: 5HTTLPR, smoking, and alcohol consumption. J Health Soc Behav 2013; 54(1): 92-108.
[24] Kranzler HR, Scott D, Tennen H, et al. The 5-HTTLPR polymorphism moderates the effect of stressful life events on drinking behavior in college students of African descent. Am J Med Genet B Neuropsychiatr Genet 2012; 159B(5): 484-90.
[25] Laucht M, Treutlein J, Blomeyer D, 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 Neuropsychopharmacol 2009; 12(6): 737-47.
[26] Pereira-Morales AJ, Adan A, Forero DA. Network analysis of multiple risk factors for mental health in young Colombian adults. J Ment Health 2019; 28(2): 153-60.
[27] Pereira-Morales AJ, Adan A, Lopez-Leon S, Forero DA. Personality traits and health-related quality of life: The mediator role of coping strategies and psychological distress. Ann Gen Psychiatry 2018; 17: 25.
[28] Jimenez KM, Pereira-Morales AJ, Adan A, Forero DA. Telomere length and childhood trauma in Colombians with depressive symptoms. Rev Bras Psiquiatr 2018.
[29] Jiménez KM, Pereira-Morales AJ, Adan A, Lopez-Leon S, Forero DA. Depressive symptoms are associated with a functional polymorphism in a miR-433 binding site in the FGF20 gene. Mol Brain 2018; 11(1): 53.
[30] Bernstein DP, Stein JA, Newcomb MD, et al. Development and validation of a brief screening version of the Childhood Trauma Questionnaire. Child Abuse Negl 2003; 27(2): 169-90.
[31] Hernandez A, Gallardo-Pujol D, Pereda N, et al. Initial validation of the Spanish childhood trauma questionnaire-short form: Factor structure, reliability and association with parenting. J Interpers Violence 2013; 28(7): 1498-518.
[32] Lang FR, John D, Lüdtke O, Schupp J, Wagner GG. Short assessment of the Big Five: robust across survey methods except telephone interviewing. Behav Res Methods 2011; 43(2): 548-67.
[33] Zung WW. A rating instrument for anxiety disorders. Psychosomatics 1971; 12(6): 371-9.
[34] Cohen S, Kamarck T, Mermelstein R. A global measure of perceived stress. J Health Soc Behav 1983; 24(4): 385-96.
[35] Middeldorp CM, Cath DC, Van Dyck R, Boomsma DI. The co-morbidity of anxiety and depression in the perspective of genetic epidemiology. A review of twin and family studies. Psychol Med 2005; 35(5): 611-24.
[36] De La Ossa S, Martinez Y, Herazo E, Campo A. Study of internal consistency and factor structure of three versions of the Zung’s rating instrument for anxiety disorders. Colomb Med 2009; 40(1): 71-7.
[37] Campo-Arias A, Oviedo HC, Herazo E. Escala de Estrés Percibido-10: Desempeño psicométrico en estudiantes de medicina de Bucaramanga, Colombia. Rev Fac Med (Caracas) 2014; 62(3): 1-24.
[38] Group WAW. The Alcohol, Smoking and Substance Involvement Screening Test (ASSIST): Development, reliability and feasibility. Addiction 2002; 97(9): 1183-94.
[39] Soto-Brandt G, Portilla Huidobro R, Huepe Artigas D, et al. Validity evidence of the Alcohol, Smoking and Substance Involvement Screening Test (ASSIST) in Chile. Adicciones 2014; 26(4): 291-302. [Validity evidence of the Alcohol, Smoking and Substance Involvement Screening Test (ASSIST) in Chile].
[40] Campo-Arias A, Villamil-Vargas M, Herazo E. Confiabilidad y dimensionalidad del AUDIT en estudiantes de medicina. Psicol Caribe 2013; 30(1): 21-35.
[41] Morales LC, Arboleda G, Rodríguez Y, et al. Absence of Lamin A/C gene mutations in four Wiedemann-Rautenstrauch syndrome patients. Am J Med Genet A 2009; 149A(12): 2695-9.
[42] Solé X, Guinó E, Valls J, Iniesta R, Moreno V. SNPStats: A web tool for the analysis of association studies. Bioinformatics 2006; 22(15): 1928-9.
[43] Hayes AF. Introduction to mediation, moderation, and conditional process analysis: A regression-based approach 2013.
[44] Norman RE, Byambaa M, De R, Butchart A, Scott J, Vos T. The long-term health consequences of child physical abuse, emotional abuse, and neglect: a systematic review and meta-analysis. PLoS Med 2012; 9(11)e1001349
[45] Organization WH. World report on violence and health 2002; 186.
[46] Shevlin M, Dorahy MJ, Adamson G. Trauma and psychosis: an analysis of the National Comorbidity Survey. Am J Psychiatry 2007; 164(1): 166-9.
[47] Duman EA, Canli T. Influence of life stress, 5-HTTLPR genotype, and SLC6A4 methylation on gene expression and stress response in healthy Caucasian males. Biol Mood Anxiety Disord 2015; 5: 2.
[48] Kang HJ, Kim JM, Stewart R, et al. Association of SLC6A4 methylation with early adversity, characteristics and outcomes in depression. Prog Neuropsychopharmacol Biol Psychiatry 2013; 44: 23-8.
[49] Caspi A, Sugden K, Moffitt TE, et al. Influence of life stress on depression: Moderation by a polymorphism in the 5-HTT gene. Science 2003; 301(5631): 386-9.
[50] Lovallo WR, Enoch MA, Yechiam E, et al. Differential impact of serotonin transporter activity on temperament and behavior in persons with a family history of alcoholism in the Oklahoma Family Health Patterns Project. Alcohol Clin Exp Res 2014; 38(6): 1575-81.
[51] McHugh RK, Hofmann SG, Asnaani A, Sawyer AT, Otto MW. The serotonin transporter gene and risk for alcohol dependence: A meta-analytic review. Drug Alcohol Depend 2010; 108(1-2): 1-6.
[52] Nilsson KW, Sjöberg RL, Damberg M, et al. Role of the serotonin transporter gene and family function in adolescent alcohol consumption. Alcohol Clin Exp Res 2005; 29(4): 564-70.
[53] Cservenka A, Brumback T. The burden of binge and heavy drinking on the brain: Effects on adolescent and young adult neural structure and function. Front Psychol 2017; 8: 1111.
[54] Bühler M, Mann K. Alcohol and the human brain: A systematic review of different neuroimaging methods. Alcohol Clin Exp Res 2011; 35(10): 1771-93.
[55] Alasmari F, Goodwani S, McCullumsmith RE, Sari Y. Role of glutamatergic system and mesocorticolimbic circuits in alcohol dependence. Prog Neurobiol 2018; 171: 32-49.
[56] Olsson CA, Byrnes GB, Lotfi-Miri M, et al. Association between 5-HTTLPR genotypes and persisting patterns of anxiety and alcohol use: Results from a 10-year longitudinal study of adolescent mental health. Mol Psychiatry 2005; 10(9): 868-76.
[57] Taylor SE, Way BM, Welch WT, Hilmert CJ, Lehman BJ, Eisenberger NI. Early family environment, current adversity, the serotonin transporter promoter polymorphism, and depressive symptomatology. Biol Psychiatry 2006; 60(7): 671-6.
[58] 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(5): 444-54.
[59] Solhan MB, Trull TJ, Jahng S, Wood PK. Clinical assessment of affective instability: Comparing EMA indices, questionnaire reports, and retrospective recall. Psychol Assess 2009; 21(3): 425-36.
[60] Ebner-Priemer UW, Kuo J, Welch SS, et al. A valence-dependent group-specific recall bias of retrospective self-reports: A study of borderline personality disorder in everyday life. J Nerv Ment Dis 2006; 194(10): 774-9.