Serotonin Transporter Gene, Depressive Symptoms, and Interleukin-6CLINICAL PERSPECTIVE
Background— We explored the relationship of genetic variants of the serotonin transporter gene SLC6A4, a key regulator of the serotonergic neurotransmission, with both depressive symptoms and plasma interleukin-6 (IL-6) levels.
Methods and Results— We genotyped 20 polymorphisms in 360 male twins (mean age, 54 years) from the Vietnam Era Twin Registry. Current depressive symptoms were measured with the Beck Depression Inventory II. IL-6 was assessed using a commercially available ELISA kit. Genotype associations were analyzed using generalized estimating equations. To study how SLC6A4 genetic vulnerability influences the relationship between depressive symptoms and IL-6, bivariate models were constructed using structural equation modeling. Of the 20 polymorphisms examined, the effective number of independent tests was 6, and the threshold of significance after Bonferroni correction was 0.008. There were 6 single-nucleotide polymorphisms significantly associated with Beck Depression Inventory (P≤0.008), including rs8071667, rs2020936, rs25528, rs6354, rs11080122, and rs8076005, and 1 single-nucleotide polymorphism was borderline associated (rs12150214, P=0.017). Of these 7 single-nucleotide polymorphisms, 3 were also significantly associated with IL-6 (P<0.008), including rs25528, rs6354, and rs8076005, and the other 4 were borderline associated (P=0.009 to 0.025). The subjects with 1 copy of the minor allele of these 7 single-nucleotide polymorphisms had higher Beck Depression Inventory scores and IL-6 levels. Further bivariate modeling revealed that ≈10% of the correlation between Beck Depression Inventory and IL-6 could be explained by the SLC6A4 gene.
Conclusions— Genetic vulnerability involving the SLC6A4 gene is significantly associated with both increased depressive symptoms and elevated IL-6 plasma levels. Common pathophysiological processes may link depression and inflammation, and implicate the serotonin pathway in neural-immune interactions.
Received April 2, 2009; accepted September 28, 2009.
The relationship between depression and the incidence of coronary artery disease (CAD) is well established.1 Adverse lifestyle behaviors, lower heart rate variability, and enhanced platelet activation have long been considered potential explanations for this association.2 Of more recent interest is the role of inflammation in the development and progression of atherosclerosis3 and its potential association with depression.4 Major depression and depressive symptoms are associated with higher levels of inflammatory biomarkers such as interleukin-6 (IL-6).5 However, the causal direction of this association remains unclear and may be bidirectional, such that neurobiological correlates of depression, or physical illness, may result in enhanced inflammation, and the latter, in turn, may increase the risk of depression.6
Clinical Perspective on p 614
It is also plausible that a common genetic vulnerability accounts for the observed association between depression and inflammation, because both phenotypes are heritable.7,8 We have recently reported evidence for a common genetic pathway linking depressive symptoms and IL-6 using twins; specifically, we found that more than two thirds of the covariance between these 2 phenotypes could be explained by the genes they have in common.8 Potential candidates for this shared genetic substrate are genetic variants involved in neurobiological pathways of the stress response, which are associated with both depression and inflammation. Among these, the serotoninergic pathway is of particular interest.
Serotonin (5-HT) is released by human mast cells and platelets after antigen challenge.9 This neurotransmitter is well known for its role in depression,10 but there is also evidence that serotonin modulates immune responses11 and proinflammatory cytokine synthesis.12 The serotonin transporter (5-HTT) is a key protein of the serotonergic system, which regulates serotonin concentration in the synaptic cleft and extrasynaptic sites.13 The human serotonin transporter gene (SLC6A4) is mapped to chromosome 17q11.1 to 17q12, and several polymorphisms within the SLC6A4 gene have been described.14 One of the most frequently studied functional polymorphisms of the SLC6A4 gene is the serotonergic transporter-linked polymorphic region (5-HTTLPR) located in the promoter region, which is a 44-bp insertion or deletion polymorphism producing 2 common alleles, long (L) and short (S) alleles.15 The homozygous or heterozygous carriage of the S allele is associated with lower transcription of 5-HTT and therefore with increased serotonin binding and reuptake.16 Some studies have found an association between 5-HTTLPR and depression,17,18 whereas others have not.19,20 Only a few studies investigated the relationship of other polymorphisms within the SLC6A4 gene and depression.21 Furthermore, no study has investigated the association of variants of the SLC6A4 gene with plasma IL-6 levels.
In this study, we investigated the association between polymorphisms of the SLC6A4 gene on depressive symptoms and IL-6 levels in a well-characterized sample of middle-aged male twins. The study of twins allowed us to quantify the role of the specific genotype of interest in explaining shared genetic effects using bivariate structural equation models (SEM).22
The Twins Heart Study (THS) is an investigation of psychological, behavioral, and biological risk factors for subclinical cardiovascular disease using twins. Twins were selected from the Vietnam Era Twin Registry,23 which includes 7369 middle-aged male-male twin pairs both of whom served in the US military during the time of the Vietnam war.
The THS included 180 twin pairs (102 monozygotic and 78 dizygotic) who were born between 1946 and 1956. The methods of construction of this sample were described in the online-only Data Supplement and were also described previously.8 Twin pairs were examined together at the Emory University General Clinical Research Center between March 2002 and March 2006 and maintained an identical schedule during examination. All twins had a comprehensive physical examination and were queried about previous diagnoses of cardiovascular diseases. We also administered the Structured Clinical Interview for DSM-IV24 to classify subjects based on a lifetime history of major depressive disorder (MDD). The final data included 105 twin pairs where neither twin had a history of MDD, 68 twin pairs discordant for history of MDD, and 7 twin pairs where both had a history of MDD. Compared with the normal twin pairs, the discordant twin pairs and the concordant-depressed twin pairs were younger, had higher prevalence of diabetes, more depressive symptoms, and more likely to use antidepressants. For other study variables, there were no significant differences among these twin pairs. Because few twins (n=8) had a current major depressive episode, we focused on current level of depressive symptoms measured with the continuous Beck Depression Inventory II (BDI-II) score. Therefore, we used all the twins in the analysis rather than study the twins based on discordant and concordant for history of MDD or current major depressive episode. This protocol was approved by the Institutional Review Board at Emory University, and informed consent was obtained from all subjects.
Assessment of Depression
Lifetime history of MDD was assessed by means of the Structured Clinical Interview for DSM-IV, which yields a clinical diagnosis of major depression based on a lifetime history of major depressive episodes.24 We also collected information on current level of depressive symptoms using the BDI-II, a standardized scale providing a continuous measure of depressive symptoms.25 It has been used extensively in community samples and has satisfactory test-retest and internal consistency reliabilities.26
Assessment of IL-6
Plasma level of IL-6 was assessed using commercially available high-sensitivity ELISA kits obtained from R&D Systems (Minneapolis, Minn). All samples were run in duplicate. Inter- and intraassay variability were reliably <10%. The values of IL-6 were log transformed to improve distribution. We were also able to evaluate the multiplicative changes due to genetic variants, ie, percent change of IL-6 for twins with the risk allele compared with those without the risk allele.
We measured weight, height, blood pressure, high- and low-density lipoprotein (HDL and LDL) cholesterol, and fasting glucose, as reported previously.8 Physical activity was assessed with a modified version of the Baecke Questionnaire of habitual physical activity; this is a 16-question instrument documenting level of physical activity at work, during sports, and nonsports activities. The global physical activity score was used in the analysis. Cigarette smoking was classified into current versus never or past smoker. A history of coronary heart disease occurring since 1990 was defined as a previous diagnosis of myocardial infarction or angina pectoris, or previous coronary revascularization procedures. Diabetes mellitus was defined as having a fasting glucose level >126 mg/dL or being treated with antidiabetic medications. Framingham risk score, a commonly used summary index of CAD risk, was also calculated to incorporate information on presence and severity of the following risk factors: age, LDL cholesterol, HDL cholesterol, blood pressure, diabetes mellitus, and current smoking.27 Information on current use of antidepressants was also collected.
Genotyping the SLC6A4 Polymorphisms and Multiple Testing
5-HTTLPR Polymorphism Genotyping
Genotyping of the 5-HTTLPR length polymorphism was carried out using polymerase chain reaction (PCR). PCR was performed in 384-well plates in a 10-μL volume with 10 ng of genomic DNA. The PCR products were mixed with a ROX 500 sizing ladder (Applied Biosystems, Foster City, Calif) then separated using an Applied Biosystems 3100 genetic analyzer and analyzed with Applied Biosystems GeneMapper 4.0 software. Fragment lengths for the L allele and S allele were 291 and 247 bp, respectively. The primers and PCR reaction conditions can be obtained by request.
Single-Nucleotide Polymorphisms Genotyping
We used the HapMap Caucasian panel (CEPH, release 24, build 36) to select SLC6A4 single-nucleotide polymorphisms (SNPs). Twenty-three SNPs with minor allele frequency (MAF) >0.05 were selected. An additional SNP (rs6355) with MAF of 0.03 was also selected because it was located in the codon region. A total of 24 SNPs were genotyped in our sample using the Beckman GenomeLab SNPstream Genotyping System (Beckman Coulter, Inc, Fullerton, Calif), whereas 5 SNPs failed to perform well in genotyping assays and were therefore discarded from further analysis. However, these 5 SNPs were in strong linkage disequilibrium (LD) with other 19 SNPs (pairwise r2 range from 0.85 to 1). Genotyping accuracy for all 19 SNPs was 99% as assessed by inclusion of duplicates (pairs of monozygotic twins) in the arrays, and negative controls (water blanks) were included on each plate.
We used the method developed by Gao et al28 for multiple testing correction because there were 20 polymorphisms being studied. Briefly, this method uses composite LD among SNPs to capture the correlation and derives the effective number of independent tests using a principal component approach. Then the Bonferroni correction can be implemented based on the number of independent tests. The method of Gao et al28 is accurate and efficient and is comparable with the permutation test. By using this method, we performed 6 independent tests from 20 correlated polymorphisms. To achieve a nominal significance of 0.05, the threshold after Bonferroni correction was 0.05/6=0.008.
The main purpose of our analyses was to test the associations between each single polymorphism and haplotype variations in the SLC6A4 gene with depressive symptoms and IL-6 plasma levels. We further investigated to what extent the covariance between depressive symptoms and IL-6 could be explained by variations in the SLC6A4 gene.
We used Haploview 4.0 software for computing Hardy-Weinberg equilibrium, MAF, and pair-wise LD between the studied polymorphisms, where 1 twin was randomly selected from each pair for the calculation.29 Association analyses were performed using generalized estimating equations (GEE), which corrects for the relatedness between twins. The exchangeable correlation structure was used for the GEE. Analyses were first performed separately for each of the studied polymorphisms and followed up by haplotype analysis. For individual polymorphism analysis, we tested the additive model of polymorphisms on depressive symptoms and IL-6 levels, respectively, where the polymorphisms were coded as 0, 1, or 2 depending on the carrier status of the minor allele. In the presence of a significant association, we further adjusted for a priori specified covariates, including body mass index, physical activity, history of coronary heart disease, use of antidepressants, and Framingham risk score. Using GEE, we also examined the association between SLC6A4 polymorphisms and MDD as a binary outcome variable.
The haplotype trend regression approach was adapted to test for associations of statistically inferred haplotypes with depressive symptoms and IL-6, respectively, as outlined by Zaykin et al,30 where linear regression was replaced with the GEE procedure to take into account the relatedness between twins. Briefly, haplotype trend regression is based on the regression of a trait on a design matrix that includes the expected proportions of haplotypes. The contributions of haplotypes are weighted with the design matrix such that unambiguous pairs of haplotypes are coded 1 for the haplotypes of homozygotes, 0.5 for each of the haplotypes of a heterozygote, and 0 for all other haplotypes. However, the contributions of ambiguous pairs of haplotypes in the design matrix are based on the probabilities of haplotype pairs (divided by 2) as estimated by PHASE 2.0.31 The most frequent haplotype was used as the reference with which the other haplotypes were contrasted. Haplotypes with estimated frequencies <5% were pooled together. The analyses were repeated after adjusting for the CAD risk factors and use of antidepressants.
We further constructed bivariate SEM, which models and takes account of residual twin covariance, to study the interplay of depressive symptoms and IL-6 with the SLC5A4 genetic vulnerability. Details of the approach for incorporating the measured genotypes in SEM have been described previously.22 Compared with univariate analysis, bivariate modeling has the advantage of allowing genetic effects on the means, variances, and relations between depressive symptoms and IL-6. The models were fitted with Mx software.32 Most of the twins from the THS are whites (94%). Further adjustment for race did not change the results. To avoid the potential influences of race/ethnicity, we repeated all the analyses after exclusion of blacks. The results were virtually identical and are not reported.
Of the 360 THS twins, 82 twins had a lifetime history of MDD. The mean duration of major depressive episodes was 8.3 months. However, the majority of twins with MDD were in remission, with only 8 subjects meeting DSM-IV criteria for a current major depressive episode and 77% having had the last depressive episode >1 year before examination. Of the 82 MDD twins, 27 (33%) were on antidepressants. There was no association between SCL6A4 polymorphisms and lifetime history of MDD. Further adjustment for the covariates did not change the results (data not shown). Table 1 shows the demographic characteristics and CAD risk factors of the THS twins. The mean age (±SD) was 54 years (±2.89), with a range of 47 to 60 years. Twenty percent of the participants were current smokers and ≈10% had a history of coronary heart disease. Fifty-two subjects (14.4%) were using antidepressants.
A total of 20 polymorphisms in the SLC6A4 gene were genotyped, including an ins/del polymorphism located in the promoter region (5-HTTLPR), a missense mutation in the codon region (rs6355), 2 SNPs in untranslated regions (UTR, rs6354 and rs1042173), and 16 SNPs in introns. The location of these polymorphisms is shown in supplemental Figure I. Mutation types on all polymorphisms and MAF are presented in Table 2, along with the Hardy-Weinberg equilibrium tests. Sixteen polymorphisms were common with MAF >10%, 2 SNPs had MAF <10% but >5%, and 2 SNPs had MAF <5%. No significant deviation from Hardy-Weinberg equilibrium was found for any polymorphism after the multiple testing corrections.
Pairwise LD coefficients of the 20 polymorphisms are displayed in supplemental Figure I. There was strong LD among rs8071667, rs12150214, rs2020936, rs25528, rs6354, rs11080122, and rs8076005 (pairwise r2 of 0.72 to 0.95). Two pairs of SNPs, rs2066713 and rs2020942 (r2=0.96) as well as rs4325622 and rs1042173 (r2=0.98), showed almost perfect LD. Strong LD was also found among rs2020939, rs140701, and rs2054847 (r2=0.85 to 0.9).
The association of each SLC6A4 polymorphism with depressive symptoms and IL-6 is shown in Table 2. No significant association was observed between 5-HTTLPR and depressive symptoms. However, 6 SNPs in other locations were significantly associated with BDI-II scores (P<0.008), including rs8071667, rs2020936, rs25528, rs6354, rs11080122, and rs8076005, and 1 SNP was borderline associated (rs12150214, P=0.017). Of these 7 SNPs, 3 were also significantly associated with IL-6 (P<0.008), including rs25528, rs6354, and rs8076005, and the other 4 were borderline associated (P=0.009 to 0.025). For these 7 SNPs, the minor allele was associated with higher BDI-II scores and IL-6 levels. Adjustment for covariates did not change the results.
The haplotype analysis using these 7 SNPs reflected the single SNP analysis because they were highly correlated. Approximately 94% of the genetic variability of the haplotypes consisting of these 7 SNPs was explained by 2 common haplotypes, including G-G-T-A-A-C-A (Hap1) and A-C-C-C-C-T-G (Hap2) with frequencies of 81% and 13%, respectively (supplemental Table I). Haplotype-based association tests are shown in Table 3. The β coefficient represents the effect size for haplotype homozygote. Compared with the subjects homozygous for Hap1, those homozygous for Hap2 had, on average, a 98% higher BDI-II score (P=0.006) and a 56% higher IL-6 level (P=0.01). Adjustment for covariates did not change the results. In addition, we performed a full haplotype analysis using all 20 polymorphisms. No significant associations were found (data not shown).
To avoid fitting multiple SNP models, we used the haplotypes consisting of these 7 SNPs for the bivariate SEM construction instead of using each single SNP. The bivariate SEM effects of the SLC6A4 risk haplotype (Hap2) on BDI-II scores (β=4.299, P=0.02) and IL-6 levels (β=0.469, P=0.01) were comparable with the β values obtained from GEE as shown in Table 3. Moreover, the bivariate SEM provided the estimations of SLC6A4 haplotypes contributing to the correlation between BDI-II and IL-6. As shown in the Figure, the explained proportions of the phenotypic variance by SLC6A4 haplotypes were 2.1% and 2.4% for BDI-II and IL-6, respectively. The overall phenotypic correlation between BDI-II and IL-6 was 0.19 (P<0.001). After taking into account the effect of SLC6A4 haplotypes, the residual correlation was estimated as 0.17 (P<0.001; see Figure), indicating that ≈10% of this correlation could be explained by the SLC6A4 haplotypes. As we have previously reported,8 more than two thirds of this correlation is due to the genetic factors. This implies that the SLC6A4 gene contributes ≈15% of the shared genetic correlation between these 2 phenotypes.
We report for the first time that genetic variations in the serotonin transporter gene SLC6A4 are associated with both depression and inflammation in middle-aged men. Participants who carried a specific haplotype consisting of 7 SNPs in the SLC6A4 gene had both significantly increased depressive symptoms and elevated IL-6 plasma levels. This risk haplotype explained ≈10% of correlation between these 2 phenotypes. Our results indicate that the SLC6A4 gene underlies a shared pathophysiological pathway linking depression and inflammation. Our data also suggest that other genetic pathways, yet to be discovered, are involved in the link between depression and inflammation.
Growing evidence has implicated the serotonergic neurotransmission pathway, and its related genes, in the pathophysiology of depression.33 The gene encoding the serotonin transporter is a fundamental regulator of serotonin signaling and therefore has been the focus of many previous studies.13 Several studies have investigated the relationship of a common functional polymorphism in the promoter region of this gene (5-HTTLPR) and depression, but the results are inconsistent.17–20 In this study, we found no association between 5-HTTLPR and depressive symptoms. It is possible that the previously reported inconsistent results may be due to an interaction between stressful life events and 5-HTTLPR, whereby carriers of the S allele are more likely to suffer from MDD but only if they had experienced multiple traumatic events.34
Only a few studies examined the relationship of other SNPs within the SLC6A4 gene and depression. A recent study examined the 5-HTTLPR and additional 4 SNPs in this gene in relation to depressive symptoms. For these 4 SNPS, no significant main effects were observed, although an interaction with stressful life events was found.21 Three of the 4 SNPs examined in this previous study (rs2020942, rs140700, and rs1042173) were also genotyped in our study. The SNP rs3794808 was not genotyped, but 3 SNPs (rs2020939, rs140701, and rs2054847) in strong LD with it (based on HapMap data) were studied in this study. No significant associations of these 6 SNPs with depressive symptoms were observed in our study.
In contrast, we found that 6 additional SNPs in the SLC6A4 gene were significantly associated with depressive symptoms and 1 SNP was borderline, where each minor allele was associated with increased BDI-II scores. These 7 SNPs were in strong LD with each other. Small to moderate LD was observed between these 7 SNPs and other polymorphisms. Six of the 7 SNPs were located in introns, and 1 SNP was in 5′ UTR (rs6354). Our study is the first one to report an association between these SNPs and depressive symptoms. Further study is necessary to explore whether these SNPs affect SLC6A4 gene expression.
We found that these same 7 SNPs were also significantly or borderline associated with IL-6 plasma levels. Substantial previous data suggest that serotonin plays a role in immune function and inflammation.11 An in vivo study showed that autologous serotonin is required for optimal T-cell activation and that the activation of suboptimally stimulated T cells can be augmented with low doses of exogenously added serotonin.35 In addition, serotonin induces NF-kB activation in human vascular smooth muscle cells and enhances IL-6 synthesis via the 5-HT2A receptor.12 Despite this evidence, no previous study examined the association between genetic variants in the serotoninergic system and IL-6 levels.
In a previous twin study using SEM, we reported that depression and inflammation share a common genetic substrate. More than two thirds of the correlation between depressive symptoms and IL-6 plasma levels was explained by shared genetic effects.8 This study extends our previous work and shows that the SLC6A4 gene plays a fundamental role in the shared genetic liability to both depressive symptoms and inflammation. We found that variants in this gene explained ≈10% of the shared phenotypic correlation of depressive symptoms and inflammation and 15% of their shared genetic correlation. These results suggest that the serotonin system is pathophysiologically implicated in both depression and inflammation, and that these 2 phenotypes may be the expression of a common pathophysiological mechanism involving neural-immune dysregulation. Our data also indicate that the remaining portion (85%) of the shared genetic correlation between depressive symptoms and IL-6 is still unexplained. Future studies should evaluate other genetic pathways that are involved in the link between depression and inflammation.
Our study is cross-sectional, thus limited in the ability to discern the temporal order between depression and inflammation. However, based on our results, the covariation of these 2 phenotypes is due in large part to a common genetic precursor rather than being a cause-effect relationship. Another limitation is that few twins met the criteria for a current major depressive episode, preventing us from examining current clinical depression. Because of this low prevalence, our analysis focused on current depressive symptoms, which are also related to CAD risk1 and which showed a graded effect with increasing depressive symptoms in this sample as previously described.8 Finally, because our twins were all middle-aged male military veterans, caution should be used in generalizing our results to women or older individuals. Our study provides the foundation for future investigations including diverse sociodemographic groups.
In conclusion, we have discovered a specific genetic vulnerability involving the serotonin transporter gene that is shared between depressive symptoms and elevated IL-6 plasma levels. Our results are consistent with the hypothesis that common pathophysiological processes link depression and inflammation and provide important insight into the role of serotonin in the pathophysiology of depression and the modulation of neural-immune interactions. Our data suggest that combining genetic and psychiatric profiling may ultimately prove useful in identifying depressed patients at highest risk for cardiovascular consequences and potentially aid in the design of novel drug targets, which may benefit clinical treatments.
The United States Department of Veterans Affairs has provided financial support for the development and maintenance of the Vietnam Era Twin Registry. Numerous organizations have provided invaluable assistance in the conduct of this study, including the Department of Defense; the National Personnel Records Center, National Archives, and Records Administration; the Internal Revenue Service; the National Institutes of Health; the National Opinion Research Center; the National Research Council, National Academy of Sciences; and the Institute for Survey Research, Temple University. Most importantly, we acknowledge the continued cooperation and participation of the members of the Vietnam Era Twin Registry and their families; without their contribution, this research would not have been possible.
Sources of Funding
This study was supported by grants K24HL077506, R01 HL68630, and R01 AG026255 from the National Institutes of Health; by grants 0245115N, 0725513B, 0730100N, and 09SDG2140117 from the American Heart Association; by grant KL2 RR025009 from the Atlanta Clinical and Translational Science Institute; by grant NIH NINR P20NR007798 from the Center for the Study of Symptoms, Symptom Interactions and Health Outcomes; and by the Emory University General Clinical Research Center MO1-RR00039.
Wulsin LR, Singal BM. Do depressive symptoms increase the risk for the onset of coronary disease? A systematic quantitative review. Psychosom Med. 2003; 65: 201–210.
Lett HS, Blumenthal JA, Babyak MA, Sherwood A, Strauman T, Robins C, Newman MF. Depression as a risk factor for coronary artery disease: evidence, mechanisms, and treatment. Psychosom Med. 2004; 66: 305–315.
McCaffery JM, Frasure-Smith N, Dube MP, Theroux P, Rouleau GA, Duan Q, Lesperance F. Common genetic vulnerability to depressive symptoms and coronary artery disease: a review and development of candidate genes related to inflammation and serotonin. Psychosom Med. 2006; 68: 187–200.
Empana JP, Sykes DH, Luc G, Juhan-Vague I, Arveiler D, Ferrieres J, Amouyel P, Bingham A, Montaye M, Ruidavets JB, Haas B, Evans A, Jouven X, Ducimetiere P. Contributions of depressive mood and circulating inflammatory markers to coronary heart disease in healthy european men: the Prospective Epidemiological Study of Myocardial Infarction (PRIME). Circulation. 2005; 111: 2299–2305.
Su S, Miller AH, Snieder H, Bremner JD, Ritchie J, Maisano C, Jones L, Murrah NV, Goldberg J, Vaccarino V. Common genetic contributions to depressive symptoms and inflammatory markers in middle-aged men: the Twins Heart Study. Psychosom Med. 2009; 71: 152–158.
Deakin JF. The role of serotonin in panic, anxiety and depression. Int Clin Psychopharmacol. 1998; 13: S1–S5.
Ito T, Ikeda U, Shimpo M, Yamamoto K, Shimada K. Serotonin increases interleukin-6 synthesis in human vascular smooth muscle cells. Circulation. 2000; 102: 2522–2527.
Lesch KP, Bengel D, Heils A, Sabol SZ, Greenberg BD, Petri S, Benjamin J, Muller CR, Hamer DH, Murphy DL. Association of anxiety-related traits with a polymorphism in the serotonin transporter gene regulatory region. Science. 1996; 274: 1527–1531.
Cervilla JA, Rivera M, Molina E, Torres-Gonzalez F, Bellon JA, Moreno B, de Dios Luna J, Lorente JA, de Diego-Otero Y, King M, Nazareth I, Gutierrez B; PREDICT Study Core Group. The 5-HTTLPR s/s genotype at the serotonin transporter gene (SLC6A4) increases the risk for depression in a large cohort of primary care attendees: the PREDICT-gene study. Am J Med Genet B Neuropsychiatr Genet. 2006; 141B: 912–917.
Mendlewicz J, Massat I, Souery D, Del-Favero J, Oruc L, Nothen MM, Blackwood D, Muir W, Battersby S, Lerer B, Segman RH, Kaneva R, Serretti A, Lilli R, Lorenzi C, Jakovljevic M, Ivezic S, Rietschel M, Milanova V, Van Broeckhoven C. Serotonin transporter 5HTTLPR polymorphism and affective disorders: no evidence of association in a large European multicenter study. Eur J Hum Genet. 2004; 12: 377–382.
First MB, Spitzer RL, Williams JBW, Gibbon M. Structured Clinical Interview for DSM IV, Patient Edition (SCID-P). Washington, DC: American Psychiatric Press; 1995.
Beck AT, Steer RA, Brown GK. Beck Depression Inventory (BDI-II). 2nd ed. San Antonio, Tex: Psychological Corporation; 1996.
Wilson PW, D'Agostino RB, Levy D, Belanger AM, Silbershatz H, Kannel WB. Prediction of coronary heart disease using risk factor categories. Circulation. 1998; 97: 1837–1847.
Gao X, Becker LC, Becker DM, Starmer JD, Province MA. Avoiding the high Bonferroni penalty in genome-wide association studies. Genet Epidemiol. 2009 May 11 [Epub ahead of print].
Barrett JC, Fry B, Maller J, Daly MJ. Haploview: analysis and visualization of LD and haplotype maps. Bioinformatics. 2005; 21: 263–265.
Neale MC, Boker SM, Xie G, Maes HH. Mx: Statistical Modeling. 6th ed. Richmond, Va: Virginia Commonwealth University Medical School; 2003.
Caspi A, Sugden K, Moffitt TE, Taylor A, Craig IW, Harrington H, McClay J, Mill J, Martin J, Braithwaite A, Poulton R. Influence of life stress on depression: moderation by a polymorphism in the 5-HTT gene. Science. 2003; 301: 386–389.
There is growing evidence that depression is associated with increased levels of inflammatory markers, but whether this link is casual remains speculative. We sought to determine whether a common genetic vulnerability may be an alternative explanation for the observed association between depression and inflammation. Consistent with this hypothesis, in a study of middle-age male twins, we found that the serotonin transporter gene (SLC6A4) is a common precursor of both depressive symptoms and circulating interleukin-6 levels. Approximately 10% of the correlation between these 2 phenotypes can be explained by the SLC6A4 genetic variants. Our results provide important insight into the role of serotonin in the pathophysiology of both depression and inflammation.
The online-only Data Supplement is available at http://circgenetics.ahajournals.org/cgi/content/full/CIRCGENETICS.109.870386.