Abstract

Genetic variants of the human serotonin transporter (SERT) may contribute to HPA axis dysregulation. SERT has two promoter region polymorphisms (5-HTTLPR: VNTR and SNP: rs25531), which may alter levels of SERT protein and its function. Combining these polymorphisms creates a functional polymorphism (FN) which may modulate mRNA expression. This study examines the relationship between these genetic variants and morning and evening salivary samples of both cortisol and dehydroepiandrosterone sulfate (DHEAS) concentrations in 269 African American (AA) adults. Resultant allele frequencies for the VNTR, SNP, and FN genotypes were 70% L (2% XLL), 84% A, and 54% LA (2% XLLA), respectively. The XLL genotype was associated with significantly higher concentrations of cortisol (~3X) and DHEAS (~2X) for both VNTR and FN polymorphisms. No significant differences were found for SNP genotypes. Conclusions were that persons with VNTR and FN XLL polymorphisms had significantly higher cortisol and DHEAS concentrations than other genotypes. AAs also appear to have a higher frequency of the rare XL allele than Caucasians. Whether the XLL genotype predisposes AAs to greater health challenges will require further research to determine the implications of these findings.

1. Introduction

The serotonergic system is a component of the central nervous system associated with differential hypothalamic-pituitary-adrenal (HPA) axis responses to stress [13]. The serotonergic system can indirectly modulate the secretion of cortisol and alter feedback regulation of cortisol and the HPA axis [47]. The HPA axis has a well-characterized response to stress and strong diurnal rhythm [8, 9]. Under normal circumstances, cortisol levels are high upon waking and reach a peak about 30 to 40 minutes after waking and steadily decline throughout the day to reach a nadir around midnight [8, 10]. Diurnal pattern deviations have been correlated with negative health and psychological outcomes [1113]. Similar diurnal patterns have also been reported for dehydroepiandrosterone sulfate (DHEAS) [14, 15]. Moreover, abnormal concentrations of cortisol and DHEAS have been linked to HPA axis dysfunction [1619].

The SLC6A4 gene codes for the human serotonin transporter (SERT; NCBI Gene ID: 6532). SERT controls the rate of serotonin reuptake from the synaptic cleft [20]. This gene has specific polymorphic changes within the promoter region, with one being the insertion of a variable number tandem repeat (VNTR) within its promoter region (5-HTTLPR) [21]. The alleles of this polymorphism are classified as short (S) for those missing the insertion, long (L) for those with the insertion, and extralong (XL) for those with an additional insertion. This particular allelic categorization appears to influence the amount of SERT mRNA expression and serotonin neurotransmission [22]. Moreover, the LL diplotype is associated with higher transcriptional activity, binding efficiency, and serotonin reuptake than the LS and SS diplotypes, resulting in more efficient clearing of serotonin from the synaptic cleft [23]. However, investigations of these polymorphisms with physiological outcomes have mixed results [24]. The inconsistent results may in part be explained by the identification of a single nucleotide polymorphism, an adenine to guanine substitution (SNP; rs25531), within the previously mentioned polymorphic region [25]. Recently these two genetic variants have been combined to create new haplotypes designated LA, LG, SA, and SG. Molecular characterization of LA, LG, SA, and SG variants in cell culture, by Hu et al. [24], clearly demonstrated that the LA haplotype is associated with approximately twice the amount of mRNA expression compared to SA. Further, work by this team identified that the A to G substitution resulted in an AP2-binding site on the L allele which acts as a transcriptional suppressor. This group also reported that the mRNA expression for SA, LG, and SG haplotypes appears to be nearly equivalent [24], and they therefore designated all three variants as S+. Together these haplotypes create functional diplotypes (FN): , +, and S+S+ which are based on transcriptional regulation.

The objective of this study was to evaluate how SERT 5-HTTLPR variants affect HPA axis function as measured by cortisol levels. In particular, we hypothesized that the S+S+ diplotype would be associated with higher morning cortisol levels as compared to the other diplotypes. Specifically, we sought to determine whether levels of morning cortisol (AM-CORT), evening cortisol (PM-CORT), and mean cortisol (M-CORT), as well as diurnal cortisol variation (DI-CORT), differed by genotype. Also, we explored differences in morning (AM-DHEAS) and evening (PM-DHEAS) dehydroepiandrosterone sulfate and mean DHEAS (M-DHEAS) by genotype. Further, as the XL allele was found in this population, we documented the association between the XL variant, cortisol, and DHEAS.

2. Methods

2.1. Participants

Participants included men and women between 18 and 60 years old of African descent (AA) who were recruited through public advertisements, announcements, and local churches in the Greater Washington DC area. Criteria for exclusion were pregnancy or the use of steroidal medications including prednisone. Subjects arrived at approximately 8 AM in the fasted condition at the Human Performance Laboratory at Uniformed Services University of the Health Sciences (HPL) for sample and data collection. The study was approved by the Institutional Review Board of the Uniformed Services University of the Health Sciences. Informed consent was obtained from all participants prior to the beginning of the study.

2.2. DNA Extraction and Genotyping

Blood samples for DNA were collected for analysis of selected SNPs as part of a larger study at the HPL research facility. Genomic DNA was successfully extracted from 254 participants as previously described [26]. Portions of DNA stock solutions were subsequently diluted to 25 ng/μL to create standardized working solutions.

Genotyping for the 5-HTTLPR VNTR and rs25531 SNP within the amplicon was performed by polymerase chain reaction (PCR) using 25 ng genomic DNA, 300 nmoles of each primer, forward 5′TCCTCCGCTTTGGCGCCTCTTCC3′ and reverse 5′TGGGGGTTGCAGGGGAGATCCTG3′, and PCR grade water (Fisher BioReagents, Fairlawn, NJ, USA) and QIAGEN multiplex PCR master mix kit (Valencia CA, USA) as per package instructions was used because it corrects for short allele amplification bias. Thermocycler (iCycler iQ4, BioRad, Hercules, CA, USA) conditions were initial denaturation at 95°C (15 min), followed by 40 cycles of 95°C (30 sec), 60°C (90 sec), and 72°C (30 sec), ending with a final elongation step at 72°C (4 min). All experiments were accomplished with positive controls and “no template” negative controls. This assay resulted in VNTR amplicons of 469 bp for S allele, 512 bp for the L allele, and ~600 bp for the XL allele as demonstrated in Figure 1.

Restriction fragment length polymorphisms (RFLPs) were generated by overnight digestions with HpaII (New England Biolabs, Ipswich, MA, USA) as per package instructions. Amplicons and RFLPs were visualized by electrophoresis in a 2% agarose gel and documented on a Gel Doc (BioRad). Digestion indicated the presence of the A or G allele for the rs25531 SNP with bands of 469 bp for the SA, 402 and 67 bp for the SG, 512 bp for the LA, and 402 and 110 bp for the LG; all XL amplicons carried the A allele and were ~600 bp in size as demonstrated in Figure 1.

Genetic assays were repeated by a separate lab (Children’s National Medical Center (CNMC)) blinded to previous results using previously described primers [23] with similar PCR conditions followed by a two-hour digestion for samples with L alleles with MspI (New England Biolabs) as per package instructions. The substitution of the G for A in the SNP produces an additional MspI recognition site (CCGG) on the long allele of the 5-HTTLPR PCR product. Samples with an A allele showed bands at 62, 127, and 340 bps. Samples with the G allele had bands at 62, 127, 166, and 174. Genotypes were then determined by running the digested PCR products out on a 2% agarose gel stained with ethidium bromide.

2.3. Hormone Analysis

Saliva samples for measuring cortisol and DHEAS were collected by the participants, one upon waking or by 8 AM (AM-CORT and AM-DHEAS) and another in the evening (PM-CORT and PM-DHEAS) between 5 and 6 PM. The samples were stored at −20°C until they were ready for shipment. Participants were instructed to refrain from eating or drinking 45 to 60 minutes before collection. Saliva samples were mailed to the investigators in prestamped envelopes and upon arrival stored at −80°C for no longer than 1 month before analysis. Vials were processed by Diagnos-Tech Inc. (Tukwila, WA, USA) for analyses of the bioactive (free fraction) portion of saliva analyte using an enzyme immunoassay procedure; values are expressed in nmol/L.

2.4. Statistics

The Hardy-Weinberg Equilibrium (HWE) equation followed by the Chi Square test was used to analyze genotypes/alleles for the normal distribution and validate HWE distribution adherence. Fisher’s Exact Test and were used to compare the allele distribution with other published studies. The categorical influences of genotype on hormone concentrations were determined by ANOVA and/or independent samples -test when applicable with 95% confidence intervals. Hormone values were reported as mean ± standard error (SE). Diurnal cortisol (DI-CORT) was calculated as the difference between AM-CORT and PM-CORT, while the mean cortisol (M-CORT) was calculated as the sum of AM-CORT and PM-CORT divided by 2 for each subject. Mean DHEAS (M-DHEAS) was calculated as the sum of AM-DHEAS and PM-DHEAS divided by 2 for each subject. The level of significance was set at and trends were recognized at values between 0.05 and 0.10. All calculations were performed using SPSS v20 software (SPSS, Inc., Chicago, IL, USA).

3. Results

3.1. Characteristics of the Study Population

The cohort included 269 individuals, 100 male and 169 female AA subjects. The mean weight was  Kg (range = 42.5 to 185 Kg) and mean age was years (range = 18 to 60 years old).

3.2. Genetic 5-HTTLPR Distribution

Genotyping was completed for 254 subjects. The distribution of the VNTR alleles was 70% L, 28% S, and 2% XL. The distribution of the SNP alleles was 84% for A and 16% for G. Both polymorphisms met Hardy-Weinberg Equilibrium (HWE) conditions (VNTR: ; ; rs25531 SNP: value of 0.01; ). Likewise, the FN genetic distribution was 54% for LA, 44% for S+, and 2% for the XL alleles and was within HWE (; ). Genotypic distributions for the study are presented in Table 1.

3.3. Salivary Cortisol and DHEAS by VNTR, SNP, and FN

Cortisol concentrations by VNTR showed that AM-CORT and PM-CORT concentrations differed significantly by VNTR genotype [ and 2.62, , resp.]. Further, the XLL AM-CORT and PM-CORT genotype had cortisol levels almost threefold that of others. This resulted in a significant difference for M-CORT levels [,   ]. DI-CORT values demonstrated a trend for a greater diurnal difference [, ]. Those with the XLL genotype were over 2 times the total study average while LS was nearly half. In the absence of the XLL genotype, the mean DI-CORT levels increased in significance [, ]. Those with the LS genotype had a significantly smaller difference between AM-CORT and PM-CORT concentrations than the other genotypes. Table 2 shows mean salivary cortisol concentrations by genotype. DHEAS analysis demonstrated that the mean AM-DHEAS differed significantly by VNTR [, ], and the XLL genotype had levels over two times higher than other genotypes. The M-DHEAS concentration was also significantly different [, ] by VNTR. No significant difference was found for other genotypes in the absence of the XLL genotype (). Mean DHEAS concentrations by genotype are presented in Table 3. Hormone concentration by SNP analysis showed no significant differences ().

Analysis by FN genotype demonstrated significant differences for AM-CORT, PM-CORT, M-CORT, AM-DHEAS, and the M-DHEAS (). The genotype had significantly higher hormone concentrations than the other three genotypes.

4. Discussion

The serotonin transporter (SERT), one of the primary modulators of serotonin function, has been studied with regard to anxiety, depression, obsessive-compulsive disorder, and stress reactivity [27]. In addition to the common VNTR alleles, we discovered a small subpopulation carrying the XL allele in our cohort. The XL allele appears to be more common in those of African descent than Caucasian [28]. In this study, those with the XLL genotype had threefold higher mean cortisol levels and twofold mean higher AM DHEAS levels than those with the LL, LS, or SS genotypes. These findings were also significant for the FN XLLA genotype, but not for the rs25531 SNP genotypes.

We hypothesized that the SS or S+S+ variants would have lower mRNA expression and produce fewer SERT, which should increase serotonin bioavailability in comparison to the other genotypes. Thus, we expected the SS/S+S+ genotypes to be associated with elevated cortisol levels. Our hypothesis was not supported by the results. However, the XLL genotype was associated with a threefold increase in cortisol and a twofold increase in DHEAS levels relative to the other genotypes. Further, the LS genotype had significantly smaller differences in AM-CORT and PM-CORT levels than the other genotypes. This flattened response was not dose dependent by genotype, as the LS and SS homozygotes show pronounced variation. A reduction in the DI-CORT variation suggests an impaired feedback regulation of cortisol and possibly the HPA axis. Further research is needed to discern this mechanism.

To our knowledge, this is the first study in which the FN was examined in humans with the XLLA allele. In work examining the , +, and S+S+ FN variants, Alexander et al. found that men with the S+S+ genotype and a significant history of stressful life events exhibited significantly greater elevations in cortisol levels in response to a stressor than those with the other genotypes [29]. Similarly, Wankerl and colleagues supported these findings but also found a significant sex by FN interaction. Higher cortisol levels were noted in men with the S+S+ genotype but not in women, relative to the other genotypes [3].

It is of interest that if we excluded the XLL genotype from our analyses, the SS group had slightly, elevated morning cortisol concentrations relative to the LS and LL genotypes, but it was not significant. Thus the association between the 5-HTTLPR FN and cortisol will require further research to determine the nature of these discrepancies.

The nonsignificant findings for the rs25531 SNP genotypes are similar to other studies [30] that support the suggestion that this SNP, when analyzed in the absence of the VNTR variable, may be a nonfunctional variant [24] with respect to cortisol level modulation. However our interpretation is limited since only seven people in our study were homozygous for the G allele. Further investigation with a larger sample size would be required to make that determination. Our FN analysis also demonstrated that those with the XLLA genotype had significantly elevated morning and evening cortisol levels. However, the same people were in the VNTR XLL and the FN XLLA groupings, and therefore the results were expected.

Limitations to our study include the small sample size for the XLL VNTR, GG SNP, and XLLA FN genotypes. Thus, our results need to be confirmed with a larger sample size. In addition, we did not control for age, menstrual status, mental health status, and use of medication/drugs including contraceptives and tobacco. All of these can alter cortisol or DHEA levels. We looked exclusively at genotype and steroid concentrations, with no noted gender effects. Due to the limited amount of sample DNA, we were not able to sequence the XL alleles to determine the genetic architecture of the VNTRs.

Our allelic distribution was within HWE and verified by our collaborators at CNMC with consistent results. We utilized reagents and methods to compensate for S allele overexpression as previously described [31]. In addition to the common VNTR alleles, a small subpopulation of this AA cohort carried the XL allele. For this study, participants with the XLL genotype had threefold higher mean cortisol and twofold higher DHEAS levels than those with the LL, LS, and SS genotypes. This was also true when classified by functional genotypes. Despite the small sample size of the XLL genotype for both the VNTR and FN groupings, this finding is of interest. Further research is needed to determine if these genotypes are predisposed to maladies associated with dysregulation of the HPA axis and conveying a causal factor for health disparities.

5. Conclusion

In summary, whereas the rs25531 SNP does not appear to be associated with either salivary cortisol or DHEAS levels in AA, people with the XLL genotype had significantly higher levels of salivary cortisol and DHEAS when analyzed by VNTR than the other genotypes. This relationship was maintained when the VNTR and rs25531 SNP polymorphisms were combined to create the FN polymorphism. Based on gene expression, those with the XLLA genotype () again had significantly higher cortisol and DHEAS concentrations than those with the , +, or S+S+ genotypes.

Disclaimer

The views expressed are those of the authors and do not necessarily reflect the official policy or position of the Uniformed Services University, Naval Medical Research Center, Department of the Navy, Department of the Army, or the United States Department of Defense.

Conflict of Interests

The authors report no conflict of interests.

Acknowledgments

This research was performed at USUHS. The authors acknowledge that this project was funded by NCMHD, National Institutes of Health, MD, Establishing Exploratory NCMHD Research Centers of Excellence (P20), and RFA-MD-07-00.