The cardiovascular and hypothalamus-pituitary-adrenal axis response to stress is controlled by glucocorticoid receptor sequence variants and promoter methylation
© Li-Tempel et al. 2016
Received: 7 October 2015
Accepted: 20 January 2016
Published: 28 January 2016
Gender, genetic makeup, and prior experience interact to determine physiological responses to an external perceived stressor. Here, we investigated the contribution of both genetic variants and promoter methylation of the NR3C1 (glucocorticoid receptor) gene to the cardiovascular and hypothalamus-pituitary-adrenal (HPA) axis response to the socially evaluated cold pressor test (seCPT).
Two hundred thirty-two healthy participants were recruited and underwent the experiment. They were randomly assigned to either the seCPT group (cold water) or a control group (warm water). The seCPT group had a clear stress reaction; salivary cortisol levels and peak systolic and diastolic blood pressure all increased significantly compared to the control group. GR genotype (TthIIII, NR3C1-I, 1H, E22E, R23K, BclI and 9beta) and methylation data were obtained from 218 participants. Haplotypes were built from the GR genotypes, and haplotype 2 (minor allele of BclI) carriers had a higher cortisol response to the seCPT in comparison to non-carriers (20.77 ± 13.22; 14.99 ± 8.42; p = 0.034), as well as independently of the experimental manipulation, higher baseline heart rate (72.44 ± 10.99; 68.74 ± 9.79; p = 0.022) and blood pressure (115.81 ± 10.47; 111.61 ± 10.74; p = 0.048). Average methylation levels throughout promoter 1F and 1H were low (2.76 and 1.69 %, respectively), but there was a strong correlation between individual CpGs and the distance separating them (Pearson’s correlation r = 0.725, p = 3.03 × 10−26). Higher promoter-wide methylation levels were associated with decreased baseline blood pressure, and when incorporated into a linear mixed effect model significantly predicted lower systolic and diastolic blood pressure evolution over time in response to the experimental manipulation. The underlying genotype significantly predicted methylation levels; particularly, the homozygous BclI minor allele was associated with higher methylation in promoter 1H (p = 0.042).
This is one of the first studies linking epigenetic modifications of the GR promoter, receptor genotype and physiological measures of the stress response. At baseline, there were clear genetic and epigenetic effects on blood pressure. The seCPT induced a strong cardiovascular and HPA axis response, and both systems were affected by the functional genetic variants, although methylation also predicted blood pressure reactivity. The return to baseline was predominantly influenced by the genomic sequence. Overall, the physiological response to the seCPT is controlled by an exquisite mix of genetic and epigenetic factors.
KeywordsGlucocorticoid receptor Alternative promoter Single nucleotide polymorphism Methylation
External challenges trigger hypothalamus-pituitary-adrenal (HPA) axis activation and cortisol secretion, maintaining homeostasis and permitting adaptation . Glucocorticoid receptor (GR, gene: NR3C1, OMIM +138040) protein isoforms and levels throughout all HPA axis tissues control glucocorticoid (GC) feedback, setting individual levels of stress reactivity and responsivity. A complex interplay of genetic and epigenetic mechanisms control GR levels, protein isoforms, and potentially the end phenotype. Twin research has suggested that part of the inter-individual differences in the stress response may be explained by genetic factors , and both rodent models and human studies show an environmental influence via epigenetic mechanisms [3–5]. Both epigenetic and genetic factors influence the transcriptional control of the GR through the series of tissue-specific promoters found upstream of the 11 alternative GR first exons [6–8].
It is well established that individual genetic variants of the glucocorticoid receptor affect both the basic cellular phenotypes i.e. GR expression levels  and the overall HPA axis stress response (reviewed in ) through either an altered GC response or sensitivity. Numerous GR SNPs are in a high linkage disequilibrium resulting in commonly accepted haplotypes  (Additional file 1: Table S1). Three haplotypes, BclI alone, TthIIII + BclI and N363S alone are all associated with an increased sensitivity to GCs [11, 12]. The N363S polymorphism was associated with increased BMI, raised cholesterol levels and an increased risk for coronary artery disease . Inversely, two haplotypes TthIIII + 9β and TthIIII + 9β + ER22/23EK have been associated with GC resistance . Importantly, there are a total of 12 known genetic variants throughout the 8 confirmed promoter regions controlling the expression of the 11 alternative first exons in the variable 5′ untranslated region (UTR) of the GR [9, 14]. This 5′UTR is responsible for controlling tissue-specific alternative first exon expression, overall GR levels and isoforms [6, 7, 15, 16].
Epigenetic modifications such as DNA methylation, post translational chromatin remodelling and small RNA-based mechanisms have more recently been shown to contribute both independently or together with genetic variation to gene regulation. DNA methylation is unique, as it is the only epigenetic mechanism that may regulate gene expression, is clearly propagated through mitosis, whilst retaining its function . Although the associations of GR promoter methylation with diseases such as posttraumatic stress disorder  and depression [5, 19–21] are well studied, there is very little evidence on how it influences HPA axis (re)activity or any other aspect of the stress response such as cardiovascular reactivity. The available studies provide inconclusive data. High methylation levels were associated with an increased salivary cortisol response in infants  and with a female-specific increased cortisol secretion after stress . Conversely, increased methylation levels were also associated with a decreased response to pharmacological HPA axis stimulation  or could be explained by differences in both education and lifestyle . The effect of DNA methylation on stress-related cardiovascular reactivity remains unexplored. These studies were limited to the proximal GR promoter regions thought to control tissue and stimuli specific GR levels [6, 16]. The mechanisms underlying the effects of DNA methylation on gene expression are not, however, particularly well understood . The longstanding association of DNA methylation with gene silencing (reviewed in ) does not reflect its functional complexity, orchestrating tissue-specific regulatory elements and expression patterns , marking alternative intra-genic promoters , controlling alternative splicing [25, 29, 30] and even promoting gene transcription [27, 31, 32]. The evidence currently available suggests that methylation in these regions of the GR do not only control the relative promoter activity, and levels of individual first exon transcripts, but also the final protein isoform and its cellular localisation [7, 16].
Genetic and epigenetic factors work together to produce the overall response, reflected in the cortisol secretion and cardiovascular system activation to an external stressor; however, neither factor acts unilaterally. Whilst epigenetic factors, particularly DNA methylation, integrate the environment experienced with the genotype , the underlying DNA sequence also has a large influence on methylation levels. In both genome-wide family-based genetic studies and HapMap cell lines, genetic variants affected DNA methylation without necessarily introducing new CpG methylation sites [34, 35]. However, the only GR data available to our knowledge suggests that there is no association between the BclI variant and methylation of the GR promoter 1F in the human placenta, although none of the other promoters or exons were investigated. However, there was a potential association between methylation, genotype and infant neurobehavioural outcomes .
In this study, we investigated the relative contribution of genetic variants and promoter methylation of the GR to both cardiovascular and HPA axis stress reactivity. In a cohort of healthy adults, stress was induced using the socially evaluated cold pressor test (seCPT), and cardiovascular reactivity was assessed from heart rate (HR) and systolic and diastolic blood pressure (SBP/DBP) changes. HPA axis reactivity was assessed from salivary cortisol. We identified associations between promoter methylation and genetic variants and analysed how these impact cardiovascular as well as HPA axis activity in response to the seCPT.
Study population and randomisation
Summary of participant characteristics and group gender repartition
22.89 ± 2.7
23.27 ± 3.0
23.01 ± 2.8
22.1 ± 2.3
22.5 ± 2.4
22.2 ± 2.3
6.1 ± 5.3
5.5 ± 4.5
6.0 ± 5.1
20.1 ± 19.3
23.2 ± 20.1
23.8 ± 19.5
25.2 ± 21.8
22.4 ± 21.3
24.3 ± 21.6
12.1 ± 15.6
7.6 ± 9.5
10.6 ± 14.1
17.8 ± 4.8
17.6 ± 5.2
17.7 ± 4.9
26.3 ± 6.5
27.6 ± 6.1
26.7 ± 6.4
103 (47.2 %)
115 (52.8 %)
218 (100 %)
70 (66.66 %)
77 (68.14 %)
147 (66.97 %)
35 (33.33 %)
36 (31.86 %)
72 (33.03 %)
seCPT induced a physiological and subjective psychological stress response
As previously reported , salivary cortisol levels were significantly increased in the seCPT group compared to the control group (p AUCg = 0.001; p AUCi < 0.001, t test). Peak SBP and DBP levels were significantly increased from baseline in the seCPT group compared to the control group. However, both SBP and DBP were comparable between the seCPT and control groups in both the baseline and recovery periods (p = 0.299 and p = 0.712, ANOVA between the seCPT and the control group). Bivariate analysis by sex showed a significantly greater increase in SBP levels in men (p = 7.02 × 10−8) and a trend towards a greater increase from baseline to peak levels (p = 0.093). There was no effect of sex on diastolic blood pressure or heart rate increase, decrease or peak levels (p > 0.1). Participants rated the seCPT significantly more stressful and had significantly higher levels of arousal, anxiety, activity and tension (all p < 0.01 paired t tests) compared to the control group. Subjective ratings were not dependent on gender (p > 0.05), although female participants tended towards increased anxiety in the seCPT group compared to control (p = 0.057).
Genotype and haplotype analysis in the cohort
Descriptive data of NR3C1 single markers in 218 participants using Haploview
Haplotype associations with HPA axis reactivity
AUCg and AUCi were used as the dependent variable in separate between-participants ANOVAs with the factors seCPT group and genotype (with each SNP and each haplotype as the genotype factor). AUCg was influenced by a significant interaction of haplotype 2 (BclI alone) and seCPT group (20.77 ± 13.22; 14.99 ± 8.42; p = 0.034). Scrutinising the structure of this interaction effect, effects analyses showed that carriers of haplotype 2 had a significantly higher AUCg than non-carriers in the seCPT group (20.77 ± 13.22; 14.99 ± 8.42; p = 0.003), whereas carriers and non-carriers did not differ significantly in the control group (Additional file 1: Table S3).
Association of haplotypes with HR and SBP
To evaluate the influences of genotypes and haplotypes on SBP and HR, a series of between-participants ANOVAs were performed. The results are summarised in Additional file 1: Table S3. There were significant main effects on the baseline, peak and recovery with higher HR for carriers of haplotype 2 (72.44 ± 10.99; p = 0.022, 74.05 ± 12.25; p = 0.023, 69.50 ± 9.48; p = 0.027, respectively) as compared to non-carriers (68.74 ± 9.79; 70.65 ± 11.74; 66.20 ± 9.28). We observed, independent of the experimental group, a significantly higher decrease of HR after the water task in carriers of haplotype 3 compared to non-carriers (6.93 ± 9.80; 3.55 ± 7.03; p = 0.016).
There was an interaction between haplotype 2 and group regarding SBP baseline (115.81 ± 10.47; 111.61 ± 10.74; p = 0.048). Only in the seCPT group, carriers of haplotype 2 had a significantly higher “SBP baseline” than non-carriers (115.81 ± 10.47; 111.61 ± 10.74; p = 0.025). This effect of haplotype may have emerged only in the seCPT group because of the higher participant number in there. There was also an interaction between haplotype 3 and group regarding the recovery SBP (112.54 ± 11.22; 116.41 ± 10.68; p = 0.019). Non-carriers in the seCPT group had higher recovery SBP than the non-carriers in the control group (116.41 ± 10.67; 112.75 ± 11.88; p = 0.049), whereas carriers in both groups did not differ.
GR promoter methylation levels and distribution in the cohort
As methylation levels correlated in clusters, promoter-wide sum methylation levels were investigated. Promoter 1H sum methylation levels were significantly higher in women than men (Mann-Whitney rank sum test, p < 0.01; Fig. 4), although there was no difference for promoter 1F sum methylation levels (Mann-Whitney rank sum test, p = 0.91; Fig. 4). This difference was maintained for methylation summed throughout the two promoters (p = 0.038, Mann-Whitney rank sum test). Although sum methylation levels for promoters 1F and 1H were not normally distributed (p < 0.001 Shapiro and Kolmogorov-Smirnov tests), there was a weak but significant Pearson’s correlation between the two promoters (r = 0.287, p = 0.65 × 10−4; Fig. 4), suggesting the clusters may also cover complete promoters.
Gender repartition after median split on methylation level
Males, n (%)
Females, n (%)
Methylation level predicts SBP and DBP
Linear mixed effects models for systolic and diastolic blood pressure
High methylation group
Test group (seCPT vs control)
High methylation group
Test group (seCPT vs control)
In the statistical model for both SBP and DBP, methylation data remained a valid predictor, despite the loss of power after dichotomisation. As a separate confirmation that the sum methylation of GR promoters 1F and 1H was significantly associated with peak SBP levels, methylation data was analysed as a continuous variable. Spearman’s correlations were performed, confirming this link (rho = −0.243, p = 0.00045; Fig. 5). However, DBP only had a trend towards associating with methylation (rho = −0.122, p = 0.095; Fig. 5).
Association of haplotype and methylation levels
As the BclI genotype is part of haplotype 4 (BclI + TthIIII + 1H), the association between methylation and haplotype 4 was analysed. Haplotype 4 tended to associate with both promoter 1F and 1H methylation levels (linear association, p = 0.067 and 0.066) although the combined methylation grouping was not linked to the haplotype (p = 0.102, chi-squared test). Similarly, haplotype 5 (TthIIII + NR3C1-I + 9beta + ER22/23EK) showed an association trend to promoter 1F methylation levels (p = 0.064, linear association), but did not associate with promoter 1H methylation levels (p = 0.921, linear association).
As the BclI genotype has been previously reported to be in LD with variants in promoter 1H (rs10482614) . Sanger sequencing of this promoter was performed. The minor alleles of rs10482614 and rs41423247 (BclI) were observable at frequencies of 29.3 and 36.2 %, respectively. As previously reported, rs10482614 was in LD with BclI (Cramer’s association coefficient, V = 0.324, p value = 3.032e−08, r 2 = 0.16 and d’ = 0.77; Fig. 2). Although the presence of the minor allele of rs10482614 (G/A) removes a CpG dinucleotide, there was no significant link between the presence of the rs10482614 minor allele and methylation of the 1H promoter (p = 0.316).
The individual response to an external stressor is dependent on a panoply of factors. Here, we report the impact of GR promoter DNA methylation and sequence variants on the physiological response to stress. Increased GR 1F and 1H methylation levels were significantly associated with decreased baseline blood pressure. GR haplotype 2 (minor allele of BclI) carriers had a higher cortisol response to the seCPT. In addition, GR haplotype 2 carriers had higher heart rate and higher blood pressure independent of experimental group. Haplotype 3 carriers had a stronger heart rate decrease post stress. A major novel finding was that the GR BclI minor allele was associated with higher GR promoter 1H methylation.
The physiological response to stress involves the sympathetic nervous system (SNS) and the HPA axis. The cold pressor test (CPT), introduced by Hines and Brown, reliably increases blood pressure, a thermoregulatory reflex as well as a global activation of the sympathetic nervous system under standardised conditions . Physiological responses, including vasoconstriction, increased skin conductance , and elevated blood pressure [41, 42] are induced. The addition of a social evaluative component in the seCPT adds a substantial HPA axis activation. However, rapid elevations in blood pressure trigger baroreflex mechanisms counteracting the heart rate increase. Consequently, blood pressure is considered the appropriate measure of cardiovascular reactivity in the seCPT , although we observed differences in both blood pressure and heart rate. The well-defined timing of the seCPT allowed us to successfully examine the baseline, immediate post-stress period and the return to baseline. At baseline, there were clear genetic and epigenetic effects on blood pressure. The seCPT induced a strong SNS and HPA axis response, and both systems were affected principally by genomic variants. The return to baseline was predominantly influenced by the genomic sequence.
Genomic variants had a significant effect on cardiovascular parameters. GR haplotype 2 (minor allele of BclI) carriers had higher baseline, peak and recovery period heart rate, and haplotype 3 carriers (minor allele of TthIIII, NR3C1-I and 9beta) had a stronger heart rate decrease post stress, both independent of the experimental group. Both of these haplotypes have previously been explored in detail, corresponding to haplotypes 4 and 2 from Cao-Lei et al.  and Kumsta et al. . Haplotype 2 appears to play a central role in determining the cardiovascular stress response. However, BclI is an intronic polymorphism, 646 bp downstream of the common exon 2 that has generally been found to associate with increased GC sensitivity [11, 44], although the mechanisms are unknown. When the expanded haplotype 2  is considered, around half of the carriers should also carry the minor alleles of the functional rs3806855 and rs3806854 in promoter 1B and rs10482614 in promoter 1H. In vitro, all three minor alleles reduced promoter 1H and 1B activity between 50 and 80 % . Methylation of the entire 1B or 1H promoter had a similar effect, reducing promoter activity by up to 90 %. Logically, carrying haplotype 2 or having high promoter 1H methylation would have similar consequences including lower GR levels and increased cardiovascular stress reactivity and activity. Although only methylation level was associated with differential cardiovascular responses to seCPT, whereas BclI/haplotype 2 influenced heart rate independent of experimental group, there will be overlap in the mechanisms underlying their actions. Nevertheless, SNPs in a high LD with those investigated in this study might be regulators of methylation and physiological traits, especially since genetic variation that leads to methylation and expression variation at the same locus is not a rare phenomenon . We hypothesise that the decreased promoter methylation observed in haplotype 2 carriers represents a counterbalance to the potential deleterious effects of the BclI genotype.
There is a well-established genetic component to variability in DNA methylation. Methylation quantitative trait loci (mQTL) are single genetic variants, often SNPs that correlate, or are associated with, DNA methylation levels. mQTLs operate over distances as large as 5 kb, occurring for approximately 2 % of the measured CpGs and 9.5 % of the expressed regions . In contrast to Bromer et al. , we observed the BclI minor allele to correlate with high sum promoter 1F methylation levels. In our linear mixed effect model of the stress response, there was a significant interaction between methylation, genotype and cardiovascular activity. Sum methylation levels for promoter 1H and 1F + 1H were higher in women than men, and methylation levels were not normally distributed in either sex. Sex-specific DNA methylation profiles not unexpected as genome-wide levels are known to be higher in males [46, 47]. However, locus specific increases are not limited to males but have also been reported for women [48–51]. Similarly, increased age has been linked to both a reduction in global methylation levels, and dramatic genome-wide redistributions of 5-mC . However, given the narrow age distribution of our participants, this was not observed. Although only methylation level was associated with differential cardiovascular responses to seCPT, whereas BclI influenced heart rate independent of experimental group, there might be a functional overlap between the two. Haplotype 2 and increased 1H methylation would both be expected to decrease promoter activity not only representing a specific GR mQTL, but also an expression methylation quantitative trait locus (emQTL) and even further a physiological expression methylation trait locus integrating the cardiovascular and stress responses with both genetic variants and methylation levels. Previous emQTL reports have all covered single CpG dinucleotides. There is currently contradicting data on the functional relevance of such limited methylation changes . For the GR, we have previously shown that complete methylation throughout each proximal GR promoter efficiently inactivates them . Similarly, methylation of a smaller (around 125 bp) fragment containing multiple CpGs also has functional effects, reducing promoter activity to ~25 % of the control, unmethylated sequence . However, there is currently no evidence that methylation of a single CpG has functional consequences on GR expression. The importance of promoter-wide changes in DNA methylation is supported by recent clinical data from subjects suffering from posttraumatic stress disorder (PTSD). Whilst Lebonté et al. identified two CpGs in GR promoter 1F that associated with PTSD, Yehuda et al. nicely demonstrated that changes occurred promoter-wide [18, 54]. This is mirrored in both the rodent maternal care paradigm and the healthy human brain. Screening chromosome 18 that contains the rat GR, differential DNA methylation was observed in clusters across broad genomic regions . At the individual CpG dinucleotide level strong distance-dependent correlations were found , further supporting our interpretation that DNA methylation changes occur in clusters and levels at individual CpGs are inter-dependent. These data lead us to suggest that our emQTL, unlike previous reports, is between haplotype 2 and a functionally relevant cluster of methylated CpGs in promoter 1H some 3 kbp upstream of the investigated region.
The generalizability and relevance of DNA methylation in peripheral blood samples to other tissues may appear questionable, as patterns are both locus and tissue specific. However, depending on the origin of the methylation patterns, it is probable that peripheral blood methylation levels are epigenetic proxies that mirror patterns in individual tissues of the cardiovascular system or the HPA axis. There are two plausible, non-exclusive mechanisms for this. Firstly, peripheral epigenetic variations may be the results of systemically acting circulating epigenetic modifiers such as cortisol . Secondly, they may originate from a commonly programmed developmental precursor tissue. DNA methylation is established de novo during embryogenesis, when it is particularly susceptible to environmental influences. Epigenetic changes across primary germ layers occurring in this period will result in levels common to several differentiated tissues . We have observed a strong correlation in methylation levels between ectoderm-derived tissues such as the anterior pituitary and the adrenal gland , as well as throughout the different neural tube derived tissues throughout the human brain  supporting the latter hypothesis. The corollary to this is that peripheral methylation levels may also be proxies for functional difference in GC sensitivity in other tissues from the same developmental origins.
The observation that haplotypes 2 and 3 have specific and different cardiovascular effects suggests that they act through different pathways. This concords with prior evidence that the renal pressure-natriuresis system and acute sympathetic activation mechanisms influence baseline cardiovascular traits and cardiovascular reactivity, respectively . However, the role of GC and the GR in these mechanisms is unclear. In GC induced hypertension, pharmacological stimulation and receptor blocking data exclude direct GC/GR interactions [59, 60], suggesting indirect mechanisms such as oxidative stress or nitric oxide deficiency [61, 62]. Nevertheless, in vitro GC have significant effects on the NO system, including reducing endothelial and inducible NOS levels, reducing l-arginine and co-factor availability as well as inhibition of transmembrane l-arginine transport [63, 64]. Our data confirms this link between GC/GR and the cardiovascular response, albeit potentially via an indirect mechanism. There was a very clear link between SBP, to a lesser extent DBP, and methylation of promoter 1B and 1H. This was confirmed by the observation that GR haplotype 3 carriers had lower blood pressure after seCPT and a higher heart rate decrease. This implies that carriers of GR haplotype 3 may be protected against hypertension to some extent, even if there does not appear to be direct GC/GR involvement. Inversely, GR haplotype 2 had a significantly higher baseline heart rate. Indeed, the constituent BclI has been linked to hypertension . Cortisol secretion (AUCg) was increased uniquely in carriers of the BclI containing haplotype 2. Our observations on the functional effects of haplotypes 2 and 3 can be generalised from our highly homogenous population to other populations, as these haplotypes have identical functional effects irrespective of ethnicity .
A weakness of our study is the limited number of SNPs that were investigated; however, the principal haplotypes previously established in the literature were readily identified. Similarly, the relatively small number of donors was counterbalanced by the high homogeneity of the young, infrequent-smoking, undergraduate student population reducing confounding socioeconomic factors.
This is one of the first studies linking epigenetic modifications of the GR promoter, receptor genotype and physiological measures of the stress response. In the baseline period prior to the water task, there were clear genetic and epigenetic effects on blood pressure, particularly the BclI containing haplotype 2 and promoter 1F and 1H methylation. This was independent of the experimental group. The water task induced a strong cardiovascular and HPA axis response in the seCPT group and both systems were affected principally by the functional genetic variants. Methylation predicted lower SBP and DBP evolution over time in response to the water task. The return to baseline was predominantly influenced by the genomic sequence. The BclI polymorphism was associated with promoter 1H methylation levels. Promoter 1F methylation levels did not associate with any of the observed genetic variants, and as such are potentially influenced by the environment. Overall, we have shown that the induction and resolution of the stress response is controlled by an exquisite mix of genetic and epigenetic factors.
Participants were recruited from the University of Trier (Germany) via e-mail and poster advertisements as previously reported . Briefly, 232 healthy non- and low-frequency smokers (<5 cigarettes per day) with a body mass index between 19 and 25 kg/m2 were recruited, and 218 (115 women and 103 men) completed the experimental protocol. Subjects with an increased objective or subjective sensitivity to cold and any indication of circulatory disturbances or cardiovascular problems were excluded. All participants completed the validated German version of The Beck Depression Inventory (BDI-II), and donors with scores above 19, consistent with moderate depression, were excluded [67–69]. As previously reported, caffeinated and alcoholic drinks, physical exercise and meals were not permitted in the 3 h immediately preceding the experimental visit . All experiments were performed between 1:30 and 6 pm. In accordance to the declaration of Helsinki, the research was approved by the ethical committee of the medical association of Rhineland-Palatinate, and all participants gave their written informed consent.
Socially evaluated cold pressor test
The socially evaluated cold pressor test (seCPT) was performed as previously reported [43, 70]. Briefly, participants assigned to the seCPT group were asked to completely immerse their hand in ice-cold (2–3 °C) water. Participants assigned to the control group were asked to completely immerse their hand in isothermic (35–37 °C) water. Participants in the seCPT group were under the social surveillance of an experimenter; their perceptions of social evaluation, uncertainty and lack of control were enhanced by warning them that the procedure may be painful, not communicating the duration of immersion during the test and informing them that their performance would be recorded for subsequent facial expression analysis. Participants assigned to the control group were not under social surveillance, and no video camera was present. All participants were asked to remove their hand from the water after 3 min. The sampling schedule is outlined in Additional file 1: Table S2. Immediately before and after cessation of the experiment, participants were asked to make a subjective rating of arousal, stress, anxiety, tension and activity on visual analogue scales ranging from 0 (“not at all”) to 100 (“very much”) in 10-point increments. Saliva samples were collected using absorbent cotton rolls (Salivette, Sarstedt, Nuembrecht, Germany). Samples were stored at −20 °C until analysis. Salivary cortisol was measured in duplicate using a time-resolved fluorescence immunoassay . As previously reported, intra-assay and inter-assay coefficients of variance were 4.0–6.7 and 7.1–9.0 %, respectively . Heart rate (HR) and blood pressure (SBP, DBP) were measured throughout the experiment using the Dinamap System (Critikon; Tampa, FL, USA) with the cuff placed on the right upper arm.
DNA was extracted from EDTA anti-coagulated blood using the salting out protocol of Miller et al. . Genomic DNA concentration was measured on a NanoDrop 1000 spectrophotometer (NanoDrop Technologies, Rockland, DE, USA). DNA was stored at −20 °C prior to bisulfite modification and pyrosequencing or genotyping.
Bisulfite modification and pyrosequencing were performed in duplicate as previously reported [15, 20, 74]. Briefly, 400-ng genomic DNA was bisulphite converted using the EpiTect-Bisulfite Kit (Qiagen) according to the manufacturer’s protocol. Promoters 1F and 1H were subsequently amplified by PCR and quantitatively pyrosequenced as previously reported [15, 20]. Pyrosequencing was performed using a PyroMark ID system, and methylation levels of each CpG dinucleotide was analysed using the Pyro Q-CpG software (version 1.0.9, Biotage). Positive controls were generated by incubation of genomic DNA with SssI, and bisulphite conversion efficiency was calculated from the conversion rate of cytosine to thymidine when not immediately followed by a guanidine as previously described [15, 20].
Genotyping and haplotype construction
The GR polymorphisms TthIIII (rs10052957), NR3C1-I (rs10482605), the promoter 1H SNP (rs10482614), ER22/23EK (rs6189 and rs6190), BclI (rs41423247) and 9beta (rs6198) were genotyped using a single nucleotide primer extension reaction, for which specific primers for each SNP were used in the SNPStart Master Mix kit from Beckman Coulter and where fragments were analysed with the CEQ8000 Genetic Analysis System (Beckman Coulter, Inc., Germany). Detailed information about primer sequences, PCR conditions and purification methods are available in supplementary information (Additional file 1: Table S2). Sanger sequencing of promoter 1H was performed as previously described . All SNPs were tested for Hardy-Weinberg equilibrium. Linkage disequilibrium (LD) was assessed for all seven SNPs using Haploview 4.2 , and LD scores were expressed as D’. Individual haplotypes were reconstructed using PHASE, version 2.1 [76, 77] (http://stephenslab.uchicago.edu/software.html#phase), which uses an algorithm based on coalescence-based Bayesian haplotype inference for predicting haplotypes from genotype data, combining modelling strategy with computational strategies.
Data reduction and statistical analysis
GR genotyping data was reduced by dichotymizing the SNPs with low minor allele frequencies combining the hetero- and homozygous carriers of the minor allele in one group. Cardiovascular data (HR, SBP, DBP) was reduced by extracting the mean increase of dependent variables from baseline to the peak after the water task and mean decrease from the peak to the recovery period. These are referred to as e.g. “SBP peak”, “SBP increase” and “SBP decrease”. Baseline was considered as average value of three measurements before water task. Recovery period was calculated using the three measurements after the water task. Cortisol data was reduced to the area under the curve with respect to ground (AUCg) and increase (AUCi; ). All variables were tested for normality graphically using kernel density plots and normal Q-Q plots and numerically using the Shapiro and Kolmogorov-Smirnov tests and the values of kurtosis and skewedness from the corresponding functions of the R package “moments”. Principal component (PCA) and internal consistency analysis (Cronbach’s alpha) were performed on all questionnaire derived data. Spherical representations of a correlation matrix and variance inflation factors (VIF) were used to identify correlations and co-linearity between covariates and explanatory variables. Variables with VIF >5 are considered co-linear and excluded from all subsequent models and analyses. Confounding factors were evaluated in a bivariate analysis for association with methylation status and test group using linear regression and Pearson’s or Spearman’s correlations or the chi-squared test, respectively. Any variable showing a significant association (p < 0.05) was included in a linear mixed effect model as a covariate. Linear mixed effects model selection was based on maximum likelihood, and an autoregressive matrix was chosen for the residuals covariance structure. All the analyses were performed using R, version 3.0.1 (The R foundation for Statistical Computing) except genotype and haplotype analyses for which general linear models (GLMs) were computed using SPSS 20.0 to assess the between-subjects effect genotype as well as the interaction time × genotype × groups for the cortisol level. Differences were considered to be significant when p < 0.05 after suitable post hoc correction in all statistical analyses. The Bonferroni correction was used for all bivariate analyses, and Tukey’s HSD was used for the repeated measures ANOVA (from R package Tukey HSD).
area under the curve
diastolic blood pressure
systolic blood pressure
socially evaluated cold pressor test
single nucleotide polymorphism
This work was financially supported by grants from the Deutsche Forschungsgemeinschaft, Germany (GRK 1389/1), the Fonds National de la Recherche (C12/BM/3985792 “EpiPath”), the Luxembourg Institute of Health (LIH) and the Ministry of Higher Education and Research of Luxembourg N°1176135.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- de Kloet ER, Joels M, Holsboer F. Stress and the brain: from adaptation to disease. Nat Rev Neurosci. 2005;6(6):463–75. doi:10.1038/nrn1683.PubMedView ArticleGoogle Scholar
- De Geus EJ, Kupper N, Boomsma DI, Snieder H. Bivariate genetic modeling of cardiovascular stress reactivity: does stress uncover genetic variance? Psychosom Med. 2007;69(4):356–64. doi:10.1097/PSY.0b013e318049cc2d.PubMedView ArticleGoogle Scholar
- Witzmann SR, Turner JD, Meriaux SB, Meijer OC, Muller CP. Epigenetic regulation of the glucocorticoid receptor promoter 1(7) in adult rats. Epigenetics. 2012;7(11):1290–301. doi:10.4161/epi.22363.PubMedPubMed CentralView ArticleGoogle Scholar
- Weaver IC, Cervoni N, Champagne FA, D’Alessio AC, Sharma S, Seckl JR, et al. Epigenetic programming by maternal behavior. Nat Neurosci. 2004;7(8):847–54.PubMedView ArticleGoogle Scholar
- McGowan PO, Sasaki A, D’Alessio AC, Dymov S, Labonte B, Szyf M, et al. Epigenetic regulation of the glucocorticoid receptor in human brain associates with childhood abuse. Nat Neurosci. 2009;12(3):342–8. doi:10.1038/nn.2270.PubMedPubMed CentralView ArticleGoogle Scholar
- Turner JD, Muller CP. Structure of the glucocorticoid receptor (NR3C1) gene 5′ untranslated region: identification, and tissue distribution of multiple new human exon 1. J Mol Endocrinol. 2005;35(2):283–92.PubMedView ArticleGoogle Scholar
- Turner JD, Vernocchi S, Schmitz S, Muller CP. Role of the 5′-untranslated regions in post-transcriptional regulation of the human glucocorticoid receptor. Biochim Biophys Acta. 2014;1839(11):1051–61. doi:10.1016/j.bbagrm.2014.08.010.PubMedView ArticleGoogle Scholar
- Leenen FA, Vernocchi S, Hunewald OE, Schmitz S, Molitor AM, Muller CP, et al. Where does transcription start? 5′-RACE adapted to next-generation sequencing. Nucleic Acids Res. 2015. doi:10.1093/nar/gkv1328.PubMedGoogle Scholar
- Cao-Lei L, Leija SC, Kumsta R, Wust S, Meyer J, Turner JD, et al. Transcriptional control of the human glucocorticoid receptor: identification and analysis of alternative promoter regions. Hum Genet. 2011;129(5):533–43. doi:10.1007/s00439-011-0949-1.PubMedView ArticleGoogle Scholar
- Spijker AT, van Rossum EF. Glucocorticoid sensitivity in mood disorders. Neuroendocrinology. 2012;95(3):179–86. doi:10.1159/000329846.PubMedView ArticleGoogle Scholar
- van Rossum EF, Koper JW, van den Beld AW, Uitterlinden AG, Arp P, Ester W, et al. Identification of the BclI polymorphism in the glucocorticoid receptor gene: association with sensitivity to glucocorticoids in vivo and body mass index. Clin Endocrinol (Oxf). 2003;59(5):585–92.View ArticleGoogle Scholar
- van Rossum EF, Lamberts SW. Polymorphisms in the glucocorticoid receptor gene and their associations with metabolic parameters and body composition. Recent Prog Horm Res. 2004;59:333–57.PubMedView ArticleGoogle Scholar
- Derijk RH, Schaaf MJ, Turner G, Datson NA, Vreugdenhil E, Cidlowski J, et al. A human glucocorticoid receptor gene variant that increases the stability of the glucocorticoid receptor beta-isoform mRNA is associated with rheumatoid arthritis. J Rheumatol. 2001;28(11):2383–8.PubMedGoogle Scholar
- Breslin MB, Geng CD, Vedeckis WV. Multiple promoters exist in the human GR gene, one of which is activated by glucocorticoids. Mol Endocrinol. 2001;15(8):1381–95.PubMedView ArticleGoogle Scholar
- Cao-Lei L, Suwansirikul S, Jutavijittum P, Meriaux SB, Turner JD, Muller CP. Glucocorticoid receptor gene expression and promoter CpG modifications throughout the human brain. J Psychiatr Res. 2013;47(11):1597–607. doi:10.1016/j.jpsychires.2013.07.022.PubMedView ArticleGoogle Scholar
- Turner JD, Alt SR, Cao L, Vernocchi S, Trifonova S, Battello N, et al. Transcriptional control of the glucocorticoid receptor: CpG islands, epigenetics and more. Biochem Pharmacol. 2010;80(12):1860–8. doi:10.1016/j.bcp.2010.06.037.PubMedView ArticleGoogle Scholar
- Zaidi SK, Young DW, Montecino M, Lian JB, Stein JL, van Wijnen AJ, et al. Architectural epigenetics: mitotic retention of mammalian transcriptional regulatory information. Mol Cell Biol. 2010;30(20):4758–66. doi:10.1128/MCB.00646-10.PubMedPubMed CentralView ArticleGoogle Scholar
- Yehuda R, Flory JD, Bierer LM, Henn-Haase C, Lehrner A, Desarnaud F, et al. Lower methylation of glucocorticoid receptor gene promoter 1F in peripheral blood of veterans with posttraumatic stress disorder. Biol Psychiatry. 2015;77(4):356–64. doi:10.1016/j.biopsych.2014.02.006.PubMedView ArticleGoogle Scholar
- Oberlander TF, Weinberg J, Papsdorf M, Grunau R, Misri S, Devlin AM. Prenatal exposure to maternal depression, neonatal methylation of human glucocorticoid receptor gene (NR3C1) and infant cortisol stress responses. Epigenetics. 2008;3(2):97–106.PubMedView ArticleGoogle Scholar
- Alt SR, Turner JD, Klok MD, Meijer OC, Lakke EA, Derijk RH, et al. Differential expression of glucocorticoid receptor transcripts in major depressive disorder is not epigenetically programmed. Psychoneuroendocrinology. 2010;35(4):544–56. doi:10.1016/j.psyneuen.2009.09.001.PubMedView ArticleGoogle Scholar
- Klok MD, Alt SR, Irurzun Lafitte AJ, Turner JD, Lakke EA, Huitinga I, et al. Decreased expression of mineralocorticoid receptor mRNA and its splice variants in postmortem brain regions of patients with major depressive disorder. J Psychiatr Res. 2011;45(7):871–8. doi:10.1016/j.jpsychires.2010.12.002.PubMedView ArticleGoogle Scholar
- Edelman S, Shalev I, Uzefovsky F, Israel S, Knafo A, Kremer I, et al. Epigenetic and genetic factors predict women’s salivary cortisol following a threat to the social self. PLoS One. 2012;7(11):e48597. doi:10.1371/journal.pone.0048597.PubMedPubMed CentralView ArticleGoogle Scholar
- Tyrka AR, Price LH, Marsit C, Walters OC, Carpenter LL. Childhood adversity and epigenetic modulation of the leukocyte glucocorticoid receptor: preliminary findings in healthy adults. PLoS One. 2012;7(1):e30148. doi:10.1371/journal.pone.0030148.PubMedPubMed CentralView ArticleGoogle Scholar
- de Rooij SR, Costello PM, Veenendaal MV, Lillycrop KA, Gluckman PD, Hanson MA, et al. Associations between DNA methylation of a glucocorticoid receptor promoter and acute stress responses in a large healthy adult population are largely explained by lifestyle and educational differences. Psychoneuroendocrinology. 2012;37(6):782–8. doi:10.1016/j.psyneuen.2011.09.010.PubMedView ArticleGoogle Scholar
- Jones PA. Functions of DNA methylation: islands, start sites, gene bodies and beyond. Nat Rev Genet. 2012;13(7):484–92.PubMedView ArticleGoogle Scholar
- Mazzio EA, Soliman KF. Basic concepts of epigenetics: impact of environmental signals on gene expression. Epigenetics. 2012;7(2):119–30. doi:10.4161/epi.7.2.18764.PubMedPubMed CentralView ArticleGoogle Scholar
- Wan J, Oliver VF, Wang G, Zhu H, Zack DJ, Merbs SL, et al. Characterization of tissue-specific differential DNA methylation suggests distinct modes of positive and negative gene expression regulation. BMC Genomics. 2015;16:49. doi:10.1186/s12864-015-1271-4.PubMedPubMed CentralView ArticleGoogle Scholar
- Maunakea AK, Nagarajan RP, Bilenky M, Ballinger TJ, D’Souza C, Fouse SD, et al. Conserved role of intragenic DNA methylation in regulating alternative promoters. Nature. 2010;466(7303):253–7. doi:10.1038/nature09165.PubMedPubMed CentralView ArticleGoogle Scholar
- Laurent L, Wong E, Li G, Huynh T, Tsirigos A, Ong CT, et al. Dynamic changes in the human methylome during differentiation. Genome Res. 2010;20(3):320–31. doi:10.1101/gr.101907.109.PubMedPubMed CentralView ArticleGoogle Scholar
- Shukla S, Kavak E, Gregory M, Imashimizu M, Shutinoski B, Kashlev M, et al. CTCF-promoted RNA polymerase II pausing links DNA methylation to splicing. Nature. 2011;479(7371):74–9. doi:10.1038/nature10442.PubMedView ArticleGoogle Scholar
- Irizarry RA, Ladd-Acosta C, Wen B, Wu Z, Montano C, Onyango P, et al. The human colon cancer methylome shows similar hypo- and hypermethylation at conserved tissue-specific CpG island shores. Nat Genet. 2009;41(2):178–86. doi:10.1038/ng.298.PubMedPubMed CentralView ArticleGoogle Scholar
- Hellman A, Chess A. Gene body-specific methylation on the active X chromosome. Science. 2007;315(5815):1141–3. doi:10.1126/science.1136352.PubMedView ArticleGoogle Scholar
- Daskalakis NP, Yehuda R. Site-specific methylation changes in the glucocorticoid receptor exon 1F promoter in relation to life adversity: systematic review of contributing factors. Front Neurosci. 2014;8:369. doi:10.3389/fnins.2014.00369.PubMedPubMed CentralView ArticleGoogle Scholar
- Bell JT, Pai AA, Pickrell JK, Gaffney DJ, Pique-Regi R, Degner JF, et al. DNA methylation patterns associate with genetic and gene expression variation in HapMap cell lines. Genome Biol. 2011;12(1):R10. doi:10.1186/gb-2011-12-1-r10.PubMedPubMed CentralView ArticleGoogle Scholar
- Gertz J, Varley KE, Reddy TE, Bowling KM, Pauli F, Parker SL, et al. Analysis of DNA methylation in a three-generation family reveals widespread genetic influence on epigenetic regulation. PLoS Genet. 2011;7(8):e1002228. doi:10.1371/journal.pgen.1002228.PubMedPubMed CentralView ArticleGoogle Scholar
- Bromer C, Marsit CJ, Armstrong DA, Padbury JF, Lester B. Genetic and epigenetic variation of the glucocorticoid receptor (NR3C1) in placenta and infant neurobehavior. Dev Psychobiol. 2013;55(7):673–83. doi:10.1002/dev.21061.PubMedPubMed CentralGoogle Scholar
- Larra MF, Schulz A, Schilling TM, Ferreira de Sa DS, Best D, Kozik B, et al. Heart rate response to post-learning stress predicts memory consolidation. Neurobiol Learn Mem. 2014;109:74–81. doi:10.1016/j.nlm.2013.12.004.PubMedView ArticleGoogle Scholar
- Kumsta R, Moser D, Streit F, Koper JW, Meyer J, Wust S. Characterization of a glucocorticoid receptor gene (GR, NR3C1) promoter polymorphism reveals functionality and extends a haplotype with putative clinical relevance. Am J Med Genet B Neuropsychiatr Genet. 2009;150B(4):476–82. doi:10.1002/ajmg.b.30837.PubMedView ArticleGoogle Scholar
- Roy-Gagnon M-H, Weir MR, Sorkin JD, Ryan KA, Sack PA, Hines S, et al. Genetic influences on blood pressure response to the cold pressor test: results from the Heredity and Phenotype Intervention Heart Study. J Hypertens. 2008;26(4):729–36. doi:10.1097/HJH.0b013e3282f524b4.PubMedPubMed CentralView ArticleGoogle Scholar
- Buchanan TW, Tranel D, Adolphs R. Impaired memory retrieval correlates with individual differences in cortisol response but not autonomic response. Learn Mem. 2006;13(3):382–7. doi:10.1101/lm.206306.PubMedPubMed CentralView ArticleGoogle Scholar
- Peckerman A, Hurwitz BE, Saab PG, Llabre MM, McCabe PM, Schneiderman N. Stimulus dimensions of the cold pressor test and the associated patterns of cardiovascular response. Psychophysiology. 1994;31(3):282–90. doi:10.1111/j.1469-8986.1994.tb02217.x.PubMedView ArticleGoogle Scholar
- Lovallo W. The cold pressor test and autonomic function: a review and integration. Psychophysiology. 1975;12(3):268–82. doi:10.1111/j.1469-8986.1975.tb01289.x.PubMedView ArticleGoogle Scholar
- Schwabe L, Haddad L, Schachinger H. HPA axis activation by a socially evaluated cold-pressor test. Psychoneuroendocrinology. 2008;33(6):890–5. doi:10.1016/j.psyneuen.2008.03.001.PubMedView ArticleGoogle Scholar
- Manenschijn L, van den Akker EL, Lamberts SW, van Rossum EF. Clinical features associated with glucocorticoid receptor polymorphisms. An overview. Ann N Y Acad Sci. 2009;1179:179–98. doi:10.1111/j.1749-6632.2009.05013.x.PubMedView ArticleGoogle Scholar
- Wagner JR, Busche S, Ge B, Kwan T, Pastinen T, Blanchette M. The relationship between DNA methylation, genetic and expression inter-individual variation in untransformed human fibroblasts. Genome Biol. 2014;15(2):R37. doi:10.1186/gb-2014-15-2-r37.PubMedPubMed CentralView ArticleGoogle Scholar
- Fuke C, Shimabukuro M, Petronis A, Sugimoto J, Oda T, Miura K, et al. Age related changes in 5-methylcytosine content in human peripheral leukocytes and placentas: an HPLC-based study. Ann Hum Genet. 2004;68(Pt 3):196–204. doi:10.1046/j.1529-8817.2004.00081.x.PubMedView ArticleGoogle Scholar
- Shimabukuro M, Sasaki T, Imamura A, Tsujita T, Fuke C, Umekage T, et al. Global hypomethylation of peripheral leukocyte DNA in male patients with schizophrenia: a potential link between epigenetics and schizophrenia. J Psychiatr Res. 2007;41(12):1042–6. doi:10.1016/j.jpsychires.2006.08.006.PubMedView ArticleGoogle Scholar
- Sandovici I, Kassovska-Bratinova S, Loredo-Osti JC, Leppert M, Suarez A, Stewart R, et al. Interindividual variability and parent of origin DNA methylation differences at specific human Alu elements. Hum Mol Genet. 2005;14(15):2135–43. doi:10.1093/hmg/ddi218.PubMedView ArticleGoogle Scholar
- Sarter B, Long TI, Tsong WH, Koh WP, Yu MC, Laird PW. Sex differential in methylation patterns of selected genes in Singapore Chinese. Hum Genet. 2005;117(4):402–3. doi:10.1007/s00439-005-1317-9.PubMedView ArticleGoogle Scholar
- Eckhardt F, Lewin J, Cortese R, Rakyan VK, Attwood J, Burger M, et al. DNA methylation profiling of human chromosomes 6, 20 and 22. Nat Genet. 2006;38(12):1378–85. doi:10.1038/ng1909.PubMedPubMed CentralView ArticleGoogle Scholar
- El-Maarri O, Becker T, Junen J, Manzoor SS, Diaz-Lacava A, Schwaab R, et al. Gender specific differences in levels of DNA methylation at selected loci from human total blood: a tendency toward higher methylation levels in males. Hum Genet. 2007;122(5):505–14. doi:10.1007/s00439-007-0430-3.PubMedView ArticleGoogle Scholar
- Liu L, van Groen T, Kadish I, Li Y, Wang D, James SR, et al. Insufficient DNA methylation affects healthy aging and promotes age-related health problems. Clinical epigenetics. 2011;2(2):349–60. doi:10.1007/s13148-011-0042-6.PubMedPubMed CentralView ArticleGoogle Scholar
- Vinkers CH, Kalafateli AL, Rutten BP, Kas MJ, Kaminsky Z, Turner JD, et al. Traumatic stress and human DNA methylation: a critical review. Epigenomics. 2015;7(4):593–608. doi:10.2217/epi.15.11.PubMedView ArticleGoogle Scholar
- Labonte B, Azoulay N, Yerko V, Turecki G, Brunet A. Epigenetic modulation of glucocorticoid receptors in posttraumatic stress disorder. Transl Psychiatry. 2014;4:e368. doi:10.1038/tp.2014.3.PubMedPubMed CentralView ArticleGoogle Scholar
- Suderman M, McGowan PO, Sasaki A, Huang TC, Hallett MT, Meaney MJ, et al. Conserved epigenetic sensitivity to early life experience in the rat and human hippocampus. Proc Natl Acad Sci U S A. 2012;109 Suppl 2:17266–72. doi:10.1073/pnas.1121260109.PubMedPubMed CentralView ArticleGoogle Scholar
- Guintivano J, Kaminsky ZA. Role of epigenetic factors in the development of mental illness throughout life. Neurosci Res. 2014. doi:10.1016/j.neures.2014.08.003.PubMedGoogle Scholar
- Petronis A. Epigenetics and twins: three variations on the theme. Trends Genet. 2006;22(7):347–50. doi:10.1016/j.tig.2006.04.010.PubMedView ArticleGoogle Scholar
- Busjahn A, Faulhaber H-D, Viken RJ, Rose RJ, Luft FC. Genetic influences on blood pressure with the cold-pressor test: a twin study. J Hypertens. 1996;14(10):1195–9.PubMedView ArticleGoogle Scholar
- Whitworth JA, Mangos GJ, Kelly JJ. Cushing, cortisol, and cardiovascular disease. Hypertension. 2000;36(5):912–6.PubMedView ArticleGoogle Scholar
- Hattori T, Murase T, Iwase E, Takahashi K, Ohtake M, Tsuboi K, et al. Glucocorticoid-induced hypertension and cardiac injury: effects of mineralocorticoid and glucocorticoid receptor antagonism. Nagoya J Med Sci. 2013;75(1-2):81–92.PubMedPubMed CentralGoogle Scholar
- Plomin R, DeFries JC, Knopik VS, Neiderheiser J. Behavioral genetics. Palgrave Macmillan; 2013.
- Ong SL, Zhang Y, Whitworth JA. Reactive oxygen species and glucocorticoid-induced hypertension. Clin Exp Pharmacol Physiol. 2008;35(4):477–82. doi:10.1111/j.1440-1681.2008.04900.x.PubMedView ArticleGoogle Scholar
- Radomski MW, Palmer RM, Moncada S. Glucocorticoids inhibit the expression of an inducible, but not the constitutive, nitric oxide synthase in vascular endothelial cells. Proc Natl Acad Sci. 1990;87(24):10043–7.PubMedPubMed CentralView ArticleGoogle Scholar
- Simmons WW, Ungureanu-Longrois D, Smith GK, Smith TW, Kelly RA. Glucocorticoids regulate inducible nitric oxide synthase by inhibiting tetrahydrobiopterin synthesis and L-arginine transport. J Biol Chem. 1996;271(39):23928–37. doi:10.1074/jbc.271.39.23928.PubMedView ArticleGoogle Scholar
- Weaver JU, Hitman GA, Kopelman PG. An association between a BclI restriction fragment length polymorphism of the glucocorticoid receptor locus and hyperinsulinaemia in obese women. J Mol Endocrinol. 1992;9(3):295–300. doi:10.1677/jme.0.0090295.PubMedView ArticleGoogle Scholar
- Melcescu E, Griswold M, Xiang L, Belk S, Montgomery D, Bray M, et al. Prevalence and cardiometabolic associations of the glucocorticoid receptor gene polymorphisms N363S and BclI in obese and non-obese black and white Mississippians. Hormones (Athens). 2012;11(2):166–77.View ArticleGoogle Scholar
- Beck AT, Steer RA. Internal consistencies of the original and revised Beck Depression Inventory. J Clin Psychol. 1984;40(6):1365–7.PubMedView ArticleGoogle Scholar
- Hautzinger M, Bailer M, Worall H, Keller F. BDI-Beck-Depressions-Inventar. Bern: Huber; 1994.Google Scholar
- Beck AT, Steer RA, Brown GK. Manual for the Beck Depression Inventory-II. San Antonio: Psychological Corporation; 1996.Google Scholar
- Lass-Hennemann J, Kuehl LK, Schulz A, Oitzl MS, Schachinger H. Stress strengthens memory of first impressions of others’ positive personality traits. PLoS One. 2011;6(1):e16389. doi:10.1371/journal.pone.0016389.PubMedPubMed CentralView ArticleGoogle Scholar
- Dressendorfer RA, Kirschbaum C, Rohde W, Stahl F, Strasburger CJ. Synthesis of a cortisol-biotin conjugate and evaluation as a tracer in an immunoassay for salivary cortisol measurement. J Steroid Biochem Mol Biol. 1992;43(7):683–92.PubMedView ArticleGoogle Scholar
- Macedo JA, Hesse J, Turner JD, Meyer J, Hellhammer DH, Muller CP. Glucocorticoid sensitivity in fibromyalgia patients: decreased expression of corticosteroid receptors and glucocorticoid-induced leucine zipper. Psychoneuroendocrinology. 2008;33(6):799–809. doi:10.1016/j.psyneuen.2008.03.012.PubMedView ArticleGoogle Scholar
- Miller SA, Dykes DD, Polesky HF. A simple salting out procedure for extracting DNA from human nucleated cells. Nucleic Acids Res. 1988;16(3):1215.PubMedPubMed CentralView ArticleGoogle Scholar
- Turner JD, Pelascini LP, Macedo JA, Muller CP. Highly individual methylation patterns of alternative glucocorticoid receptor promoters suggest individualized epigenetic regulatory mechanisms. Nucleic Acids Res. 2008;36:7207–18.PubMedPubMed CentralView ArticleGoogle Scholar
- Barrett JC, Fry B, Maller J, Daly MJ. Haploview: analysis and visualization of LD and haplotype maps. Bioinformatics. 2005;21(2):263–5. doi:10.1093/bioinformatics/bth457.PubMedView ArticleGoogle Scholar
- Stephens M, Donnelly P. A comparison of Bayesian methods for haplotype reconstruction from population genotype data. Am J Hum Genet. 2003;73(5):1162–9. doi:10.1086/379378.PubMedPubMed CentralView ArticleGoogle Scholar
- Stephens M, Smith NJ, Donnelly P. A new statistical method for haplotype reconstruction from population data. Am J Hum Genet. 2001;68(4):978–89. doi:10.1086/319501.PubMedPubMed CentralView ArticleGoogle Scholar
- Pruessner JC, Kirschbaum C, Meinlschmid G, Hellhammer DH. Two formulas for computation of the area under the curve represent measures of total hormone concentration versus time-dependent change. Psychoneuroendocrinology. 2003;28(7):916–31.PubMedView ArticleGoogle Scholar