Reduction of Glucocorticoid Receptor Function in Chronic Fatigue Syndrome
Glucocorticoid receptor (GR) function may have aetiopathogenic significance in chronic fatigue syndrome (CFS), via its essential role in mediating inflammatory responses as well as in hypothalamic-pituitary-adrenal axis regulation. GR function can be estimated ex vivo by measuring dexamethasone (dex) modulation of cytokine response to lipopolysaccharide (LPS), and in vivo using the impact of dex on cortisol levels. This study aimed to compare the GR function between CFS (), primary Sjögren’s syndrome (a disease group control) (), and sedentary healthy controls (HCs) (), and to investigate its relationship with clinical measures. In the GR ex vivo response assay, whole blood was diluted and incubated with LPS (to stimulate cytokine production), with or without 10 or 100 nanomolar concentrations of dex. Cytometric bead array (CBA) and flow cytometry enabled quantification of cytokine levels (TNFα, interleukin- (IL-) 6, and IL-10) in the supernatants. In the in vivo response assay, five plasma samples were taken for determination of total cortisol concentration using ELISA at half-hourly intervals on two consecutive mornings separated by ingestion of 0.5 mg of dex at 11 pm. The association of the data from the in vivo and ex vivo analyses with reported childhood adversity was also examined. CFS patients had reduced LPS-induced IL-6 and TNFα production compared to both control groups and reduced suppression of TNFα by the higher dose of dex compared to HCs. Cortisol levels, before or after dex, did not differ between CFS and HCs. Cortisol levels were more variable in CFS than HCs. In the combined group (CFS plus HC), cortisol concentrations positively and ex vivo GR function (determined by dex-mediated suppression of IL-10) negatively correlated with childhood adversity score. The results do not support the hypothesis that GR dysregulation is aetiopathogenic in CFS and suggest that current and future endocrine cross-sectional studies in CFS may be vulnerable to the confounding influence of childhood trauma which is likely increased by comorbid depression.
Chronic fatigue syndrome (CFS) has a prevalence of 2% in the UK [1, 2]. It is defined by profound, persistent, medically unexplained fatigue lasting at least 6 months, which is not caused by ongoing exertion, not significantly eased by rest, and is severe enough to cause considerable loss of function [3–5]. Alongside this are symptoms of inflammation, pain, cognitive deficits, and psychiatric and bowel problems . Often, biological tests and physical examinations are unremarkable. CFS affects all ages and the peak age of onset is 20–40. Full recovery is rare [2, 6] and comorbidity with depression is common.
Many putative causes of CFS have been investigated but the absence of an agreed pathogenesis impacts the development of effective diagnostics and treatments. It is likely that multiple factors contribute which involve a number of interacting biological, environmental, and psychosocial factors [1, 2, 5, 7–11].
The recognised temporal relationship between stressors and the onset and course of CFS suggests an aetiopathogenic role for systems controlling the stress response including the sympathetic nervous system and the hypothalamic-pituitary-adrenal (HPA) axis [4, 8, 12–17]. The functionality of the glucocorticoid receptor (GR), determined by its sensitivity, affinity, and density and by its interaction with transcription factors , is arguably the defining factor in HPA axis regulation  and responsible, in a large part, for basal concentrations of cortisol throughout the day. In HPA axis downregulation in response to the GR agonist, dexamethasone (dex) is the most commonly used in vivo method of determining GR function. An ex vivo technique dependent on the inhibitory effect of GR activation on cytokine release is also utilized [20, 21]. Cross-sectional studies in CFS tend to show basal hypocortisolaemia [22, 23], attenuated diurnal variation [4, 24], an attenuated response to activation by CRH or ACTH [16, 23, 25–28], an enhanced suppression by dex [10, 17, 28–30], and an enhanced dex-induced suppression of IL-6, TNFα, IL-10, and IL-4 synthesis [21, 31, 32] and of peripheral blood mononuclear cell proliferation . Genetic studies in CFS have shown the salience of functional single nuclear polymorphisms in NR3C1 [31, 32], which codes for GR, and have also shown hypomethylation of the 1F promotor region of this gene [33–35]. The endocrine findings however are not consistent, and elevated  or normal [18, 27, 25, 37] basal cortisol levels have also been reported, as has a normal cortisol response to wakening  and to psychosocial stress [21, 38]. There is also some evidence of an association of HPA axis dysregulation with symptom severity and poorer prognosis [15, 39] but only speculation about the mechanism through which HPA axis abnormalities may result in the symptoms of CFS; glucose supply , hypotension with associated reduced cerebral perfusion , and CRH-induced appetite and sleep disturbance have been considered .
The HPA axis interacts with many other systems, notably the immune system . Glucocorticoids (GCs) modulate immune responses by altering gene expression, transcription, translation, and protein secretion , either directly, by decreasing transcription of the genes which code for cytokines, or indirectly, by inhibiting proinflammatory transcription factors . GCs inhibit, with different sensitivities, cytokines such as IL-6, IL-1, and TNF (TNF the most, IL-6 the least) [42, 43]. Immune activation, with an increase in proinflammatory cytokine concentrations, including IL-6, TNFα, and IL-1 [2, 44, 45], has been rather inconsistently, for example , demonstrated in CFS, may be secondary to insufficient glucocorticoid signaling , and may result in pain, fatigue, cognitive deficits, and other symptoms which are characteristic of CFS [12, 16, 38, 46].
In order to examine the nature, extent, and impact of HPA axis dysregulation in CFS, we sought to compare GR function using both in vivo and ex vivo assessment. We therefore examined the HPA axis and immune system function in a sample of patients with CFS and in healthy comparators and in participants with the systemic autoimmune condition, primary Sjögren’s syndrome (pSS), who acted as disease group comparators.
Three groups were recruited. The study was carried out in accordance with the Declaration of Helsinki. The study design was approved by the Newcastle and North Tyneside Ethics Committee. All participants provided written informed consent. Participants were aged 22–68 years old. Exclusion criteria consisted of age < 18 years, a current or past axis I psychiatric diagnosis confirmed using the Structured Clinical Interview for DSM-IV [47, 48], and use, in the 72 hours prior to enrolment, of antihypertensives, antidepressants, or analgesics. Samples were collected as part of an MRC-funded cohort study (MRC MR/J002712/1). 48 participants with CFS (13 males (mean age 52.2) and 35 females (mean age = 44.9)) were recruited via the local CFS clinical service, all fulfilled the Fukuda diagnostic criteria, and had a mean FIS of 88 and CTQ of 32. Twenty healthy comparators (HC; 7 males (mean age = 43.1) and 13 females (mean age 44.9)) were recruited from a HC database, word of mouth, social media, and advertisement in the hospital (mean FIS = 4, CTQ = 29). HCs were age and sex matched to the patients, and attempts were made to match on activity levels using the Mean International Physical Activity Questionnaire (IPAQ) although the mean CFS IPAQ rating was “low” and the HC rating was “medium.” Primary Sjögren’s syndrome (pSS) patients () fulfilled the American European Consensus Group classification  and were recruited from the United Kingdom pSS Registry .
2.2. Symptom Assessment Tools
The CFS participants completed the Fatigue Impact Scale (FIS)  and the Childhood Trauma Questionnaire-Short Form (CTQ-SF) . The FIS quantifies individual perception of the impact that fatigue has on daily functioning (Fisk et al.) . There are 40 items, each scored on a 5-point Likert scale providing a continuous scale of 0–160. It comprises three subscales looking at the impact that fatigue has on physical (10 items: motivation, effort, stamina, and coordination), psychosocial (20 items: isolation, emotions, coping, and workload), and cognitive (10 items: concentration, memory, and thinking) functioning. A higher score indicates greater fatigue. The CTQ is a 28 item self-report scale which measures the frequency and severity of childhood adversity. It consists of five factors: emotional abuse, physical abuse, sexual abuse, emotional neglect, and physical neglect, and possesses good psychometric properties . Items are scored on a Likert scale, with responses ranging from 1 (“never true”) to 5 (“very often true”).
2.3. In Vivo Assessment of HPA Axis Function by Measurement of Cortisol Levels in Response to Low-Dose dex in CFS Compared to Healthy Controls
In the 48 CFS patients and the first 10 healthy controls (HCs), plasma samples were taken in lithium-heparin vacutainers at 30-minute intervals between 10 am and noon on two consecutive days (day 1 and day 2). At 11 pm on day 1, participants took oral dex (0.5 mg). Practical consideration meant that only the baseline blood (10 am day 1 sample) were taken for the other 10 HCs. Within one hour of collection, the blood was spun at 1600 g for 10 minutes at room temperature. Aliquots of plasma were extracted and stored at −80°C until analysis. Plasma cortisol concentrations were quantified using 15 lot-matched cortisol ELISA kits, supplied by Abcam and used according to the manufacturer’s protocol. The lower limit of cortisol detection was 2.44 ng/ml.
2.4. Ex Vivo Measurement of GR Function: The Glucocorticoid Receptor Response Assay
The GR response assay  utilised lithium-heparin-treated blood taken at 10 am on day 1. It was set up in sterile 48-well plates within 3 hours of collection. Blood was diluted 1/10 with room temperature RPMI 1640 containing penicillin-streptomycin and L-glutamine, mixed thoroughly by inversion, and then incubated for 24 hours at 37°C in a humidified atmosphere of 95% air and 5% CO2. There were four conditions, “null” (medium alone), “LPS” in which 50 μl of a 200 μg/ml LPS solution was added to the 400 μl of diluted blood (to stimulate cytokine production from cells), “dex10” in which 50 μl of a 100 nanomolar (nM) dex solution was also added, and “dex100” in which a tenfold stronger dex solution was added. After incubation, the assay plates were spun in a 4°C centrifuge for 10 minutes at 1000 RPM. 300 μl supernatant samples were then harvested for each condition, transferred to 0.6 ml microcentrifuge tubes, and stored at −20°C. Cytokine concentrations (TNFα, IL-6, and IL-10) were determined using Cytometric Bead Array (CBA) and flow cytometry according to the manufacturer’s instruction (BD Bioscience). Percentage suppression on dex 10 nM (% dex10) and on dex 100 nM (% dex100) was calculated using the following equation:
2.5. Baseline Cytokine Levels
Gel-based specimen tubes were used to collect a further serum sample at 10 am on day 1 for baseline measurement of a range of inflammatory markers. These were spun within 3 hours at 1600 RPM for 10 minutes at room temperature. 2 × 1 ml of serum was extracted and stored at −80°C until analysis which utilized the method described above, that is, CBA and flow cytometry. Two different dilutions of serum samples were required for the measurement of different cytokines; this was therefore conducted in 2 batches for CFS (), HC (), and pSS () participants.
2.6. Data and Statistical Analysis
Statistical tests were carried out using SPSS version 23 and the “R” statistical package. Graphs were produced using GraphPad Prism version 5.01 and MATLAB. values are two-tailed with significance set at .
Repeated measures ANOVA was conducted with time (5) as within and group (2) as between factors to examine cortisol concentrations on each of the two days. Area under the curve (AUC) was also calculated using trapezoid integration, both for day 1 and day 2. Specifically, AUC with respect to ground (AUCg), considered to be a measure of overall cortisol output including baseline activity, and AUC increase (AUCi), a putative measure of the sensitivity of GR to modulation, were calculated . The difference between AUCg on the 2 days was also calculated (delta AUCg). Shapiro-Wilk and QQ plots (data not shown) revealed that even after Box-Cox transformation, neither AUC nor cytokine data met the assumptions required for ANOVA; thus, nonparametric comparisons were used. Spearman correlations were conducted to examine the relationship between childhood adversity and endocrine parameters. All data, shown or not shown, is available for scrutiny upon request.
3.1. In Vivo Assessment of HPA Axis Function by Measurement of Cortisol Levels in Response to Low-Dose dex in CFS Compared to Healthy Controls
There was a significant effect of time on day one (, , ) but not on day two (, , ). No effect of group on either day (day one, , , ; day 2, , , ) (see Figure 1). Cortisol AUCs (g or i, day one, day two, or delta) did not differ ().
3.2. Ex Vivo Assessment of HPA Axis Function: The Glucocorticoid Receptor (GR) Response Assay
GR response assay blood was not taken for two CFS patients. One HC was removed due to an abnormally high null value and one CFS patient removed due to not stimulating sufficiently on LPS. Analysis was therefore conducted on samples taken from CFS (), HC (), and pSS () participants.
In the null sample, pSS participants had higher IL-10 levels than the CFS or HC participants but there were no differences between groups for IL-6 or TNFα. LPS induced a robust cytokine response, and after LPS, group differences were evident such that, for the positive cytokines IL-6 and TNFα, pSS participants had higher levels than CFS participants, who had higher levels than HCs. For IL-10, the difference was between pSS (higher) and HCs. Median cytokine levels were, on the whole, lower in the dex 10 nM samples and invariably in the dex 100 nM samples than the LPS alone samples. Percentage suppression was greater for samples incubated with dex 100 nM than those incubated with dex 10 nM. Percentage suppression with dex 10 nM or dex 100 nM was not different between the groups, except for a greater suppression in the dex 100 nM condition determined using TNFα (see Table 1).
Individual AUC and cytokine values were visualized using frequency density plots. Visual inspection suggested greater variability in patients compared to HCs with some patients showing less suppression (see Figure 2). In order to investigate this, we performed post hoc -tests of standard deviations derived from Bayesian hierarchical models of outcome measures under Gaussian (normal) priors. The advantage of this approach lies in its treatment of parameters as sampling variables rather than population attributes which allows us to formally compare modelled estimates of their values. Significance was assessed by comparing the 95% credible intervals of the posterior distributions. Vague priors were used for the mean and standard deviation, and analysis was carried out using the BEST package in the R statistical environment . Results are displayed in Table 1. There was no difference in standard deviation of day 1 AUCg distributions. There was a significant difference in standard deviations for day 1 AUCi (), day 2 AUCg (), and day 2 AUCi (). There was also a significant difference in standard deviations of TNF distributions (). The difference in IL-6 was marginal () though there was no significant difference in IL-10 distributions ().
3.3. The Relationship between Reported Adversity and HPA Axis Function in CFS and Controls
Correlation coefficients, in CFS participants, reveal a negative relationship between the CTQ score (and the emotional subscores) and cortisol AUC but no significant relationship with % suppression (after incubation with 10 nm dex). In HCs, there was no significant relationship between CTQ scores and cortisol AUC, but there was a positive relationship with percentage IL-6 suppression (see Table 2). In the combined sample, of CFS and HC (see Table 3), the significant correlations were between the CTQ total score and AUCg (positive) and IL-10 (negative).
3.4. Baseline Cytokine Levels
CFS patients showed reduced production of IP-10 and IL-12/23p40 compared to HC and of IP-10, MIP1α, IL-6, and IL-1β compared to pSS participants. pSS participants showed increased production of MIP1α and IL-6 compared to HV. A violin plot was designed post hoc using “R” statistical software to visualize variance between populations and determine whether subpopulations were present and is displayed in Figure 3.
We did not demonstrate a difference in cortisol levels between participants with CFS and healthy volunteers. This differs from the majority of cross-sectional HPA axis studies in this population (see “Introduction” and  for review). This difference may be related to the population; the sample here, for instance, was rigorously screened for comorbid depression, and Papadopoulos et al.  have previously demonstrated that dex-induced cortisol suppression differed only in CFS patients with comorbid depression or it may be a type II error consequent on the small sample size combined with the marked variation in cortisol levels in CFS as highlighted by the frequency density graph and the significantly greater cortisol variability in patients with CFS. The aetiopathogenic relevance of this variability is unknown but it suggests a lack of precision in cortisol regulation . The heterogeneity in cortisol concentrations may suggest clinical heterogeneity within the diagnostic grouping of CFS and emphasizes the impact of disparate and competing factors on GR function including current and previous stressors, the common use of antidepressants [20, 54] (even in those who have never met criteria for major depressive disorder) [20, 57], and the impact of a primary dysregulation of proinflammatory cytokines .
The baseline cytokine data emphasized the status of pSS as an inflammatory disorder. The ex vivo data revealed a reduced capacity for a proinflammatory cytokine response to LPS in CFS compared with HCs (and an increased responsivity compared to the pSS participants). It further revealed that (independent of group) incubation with dex, in a dose-dependent manner, as expected, suppressed cytokine release. The percent suppression of LPS-induced TNFα release by 100 nM solution of dex was less in CFS patients than HCs. This may be suggestive of reduced GR function in CFS but any such interpretation must be made with caution as the impact of 10 nM dex did not significantly differ between CFS and HCs, neither was a significant effect seen when IL-6 or IL-10 was used as the output variable. That TNFα was most sensitive to suppression by dex accords with the existing literature , is congruent with the theory that GCs may preferentially inhibit Th1 over Th2 cells , and suggests that TNFα may be the most appropriate cytokine for GR response assay studies in CFS. The variability in percentage suppression of cytokine levels by dex is also greater in CFS than healthy or pSS controls.
In CFS participants, there was a relationship between the score on the childhood trauma questionnaire and cortisol AUCs such that higher reported levels of early adversity correlated negatively with cortisol. Interestingly, a different pattern was seen in healthy volunteers in whom reported childhood adversity associated positively with dex-induced IL-6 suppression in the absence of an effect of cortisol concentrations. When the groups were combined to maximise power, a negative relationship between reported adversity and cortisol levels and a positive relationship with GR function (here shown using IL-10 not IL-6) are revealed.
The variability in cortisol has implications for the interpretation of existing and future endocrine cohort studies in CFS because of the associated risk of type I and type II errors; our data, for instance, would suggest that the proportion of participants in a sample who experienced childhood adversity will be expected to determine the likelihood that basal hypocortisolaemia will be shown. In addition to the CTQ total score, we report here also the emotional neglect subscale, having previously argued that the pervasive nature of emotional neglect ensures that it enacts the greater sustained impact on behavioural and endocrine function [60, 61].
There can be few who argue with the notion that, in the general population, early adversity, acting for instance through methylation or other epigenetic mechanisms, impacts GR function and so GR mediated negative feedback on the HPA axis and thus cortisol synthesis [62, 63] and, further, that this has relevance for understanding the pathophysiology of mood disorders . It is of interest, here, to consider the implication that this has for our understanding of the pathophysiology of CFS and for the interpretation of endocrine studies in this population. We have previously postulated that childhood adversity is not a risk factor for CFS per se, but it can appear to be because of the impact of comorbid or misdiagnosed depression . It has been further conjectured that comorbid depression may commonly confound CFS studies  and, just as this may lead to erroneous finding of increased rates of childhood adversity in CFS, similarly, it may explain the methylation pattern  including in the NR3C1-1F promoter region , the increased GR function (shown using the DST, the dex/CRH test , or ex vivo measures), and the basal hypocortisolaemia which have been (inconsistently ) shown in previous CFS studies. In this current study, CTQ scores were not greater in the CFS participants than HCs, and it is interesting to note that the basal cortisol or GR function as determined by post-dex cortisol or dex-induced suppression of cytokine synthesis was not convincingly different.
Despite our rigorous exclusion of those who met the diagnostic criteria for depression and the lack of difference in childhood adversity reported by CFS patients compared with HCs, there was a signal that HPA axis regulation was different in CFS; the variability of pre- and post-dex cortisol levels and of dex-induced cytokine suppression was increased in CFS, the proinflammatory cytokine response to LPS was attenuated, and the TNFα suppression by the larger dex dose was greater, and, whilst we do not want to make too much of this, the graph suggested (but the stats did not back up) the possibility that post-dex cortisol was lower in CFS than HCs. Further research is needed to understand the cause and significance of this data; this will need large, well-characterised groups and will need consideration to be given to the interacting networks of biological, psychological, and social factors.
The data used to support the findings of this study are available from the corresponding author upon request.
Conflicts of Interest
The authors declare that they have no conflicts of interest.
This study was funded by the Medical Research Council (MR/J002712/1), by ME Research UK, and by Action for ME. Thanks are due to the participants and to Heather Slater for the help with the standard deviation graphs.
A. J. Parker, S. Wessely, and A. J. Cleare, “The neuroendocrinology of chronic fatigue syndrome and fibromyalgia,” Psychological Medicine, vol. 31, no. 8, pp. 1331–1345, 2001.View at: Google Scholar
N. G. Klimas, F. R. Salvato, R. Morgan, and M. A. Fletcher, “Immunologic abnormalities in chronic fatigue syndrome,” Journal of Clinical Microbiology, vol. 28, no. 6, pp. 1403–1410, 1990.View at: Google Scholar
T. Theorell, V. Blomkvist, G. Lindh, and B. Evengard, “Critical life events, infections, and symptoms during the year preceding chronic fatigue syndrome (CFS): an examination of CFS patients and subjects with a nonspecific life crisis,” Psychosomatic Medicine, vol. 61, no. 3, pp. 304–310, 1999.View at: Publisher Site | Google Scholar
J. Visser, E. Lentjes, I. Haspels et al., “Increased sensitivity to glucocorticoids in peripheral blood mononuclear cells of chronic fatigue syndrome patients, without evidence for altered density or affinity of glucocorticoid receptors,” Journal of Investigative Medicine, vol. 49, no. 2, pp. 195–204, 2001.View at: Publisher Site | Google Scholar
U. M. Nater, E. Maloney, R. S. Boneva et al., “Attenuated morning salivary cortisol concentrations in a population-based study of persons with chronic fatigue syndrome and well controls,” The Journal of Clinical Endocrinology & Metabolism, vol. 93, no. 3, pp. 703–709, 2008.View at: Publisher Site | Google Scholar
J. Visser, W. Graffelman, B. Blauw et al., “LPS-induced IL-10 production in whole blood cultures from chronic fatigue syndrome patients is increased but supersensitive to inhibition by dexamethasone,” Journal of Neuroimmunology, vol. 119, no. 2, pp. 343–349, 2001.View at: Publisher Site | Google Scholar
J. Visser, B. Blauw, B. Hinloopen et al., “CD4 T lymphocytes from patients with chronic fatigue syndrome have decreased interferon-gamma production and increased sensitivity to dexamethasone,” The Journal of Infectious Diseases, vol. 177, no. 2, pp. 451–454, 1998.View at: Publisher Site | Google Scholar
E. Vangeel, F. Van Den Eede, T. Hompes et al., “Chronic fatigue syndrome and DNA hypomethylation of the glucocorticoid receptor gene promoter 1F region: associations with HPA axis hypofunction and childhood trauma,” Psychosomatic Medicine, vol. 77, no. 8, pp. 853–862, 2015.View at: Publisher Site | Google Scholar
B. N. Goertzel, C. Pennachin, L. de Souza Coelho, B. Gurbaxani, E. M. Maloney, and J. F. Jones, “Combinations of single nucleotide polymorphisms in neuroendocrine effector and receptor genes predict chronic fatigue syndrome,” Pharmacogenomics, vol. 7, no. 3, pp. 475–483, 2006.View at: Publisher Site | Google Scholar
A. Papadopoulos, M. Ebrecht, A. D. L. Roberts, L. Poon, N. Rohleder, and A. J. Cleare, “Glucocorticoid receptor mediated negative feedback in chronic fatigue syndrome using the low dose (0.5 mg) dexamethasone suppression test,” Journal of Affective Disorders, vol. 112, no. 1-3, pp. 289–294, 2009.View at: Publisher Site | Google Scholar
J. C. Buckingham, H. D. Loxley, H. C. Christian, and J. G. Philip, “Activation of the HPA axis by immune insults: roles and interactions of cytokines, eicosanoids, and glucocorticoids,” Pharmacology, Biochemistry, and Behavior, vol. 54, no. 1, pp. 285–298, 1996.View at: Publisher Site | Google Scholar
M. B. First, R. L. Spitzer, M. Gibbon, and J. B. W. Williams, Structured Clinical Interview for DSM-IV Axis I Disorders, Research Version. Biometrics Research, New York, 1997.
M. B. First, R. L. Spitzer, M. Gibbon, and W. JBW, “Structured Clinical Interview for DSM-IV Axis I disorders-Patient Edition,” in SCID-I/P, Version 2.0, Biometrics Research Department, New York, NY, USA, 1995.View at: Google Scholar
J. C. Pruessner, C. Kirschbaum, G. Meinlschmid, and D. H. Hellhammer, “Two formulas for computation of the area under the curve represent measures of total hormone concentration versus time-dependent change,” Psychoneuroendocrinology, vol. 28, no. 7, pp. 916–931, 2003.View at: Publisher Site | Google Scholar
L. A. Carvalho, B. A. Garner, T. Dew, H. Fazakerley, and C. M. Pariante, “Antidepressants, but not antipsychotics, modulate GR function in human whole blood: an insight into molecular mechanisms,” European Neuropsychopharmacology, vol. 20, no. 6, pp. 379–387, 2010.View at: Publisher Site | Google Scholar
S. Watson, B. M. Owen, P. Gallagher, A. J. Hearn, A. H. Young, and I. N. Ferrier, “Family history, early adversity and the hypothalamic-pituitary-adrenal (HPA) axis: mediation of the vulnerability to mood disorders,” Neuropsychiatric Disease and Treatment, vol. 3, no. 5, pp. 647–653, 2007.View at: Google Scholar