Long Term
Long Term
https://doi.org/10.1210/clinem/dgae081
Advance access publication 19 March 2024
Clinical Research Article
Abstract
Context: Long-term glucocorticoid levels in scalp hair (HairGCs), including cortisol and the inactive form cortisone, represent the cumulative
systemic exposure to glucocorticoids over months. HairGCs have repeatedly shown associations with cardiometabolic and immune
parameters, but longitudinal data are lacking.
Design: We investigated 6341 hair samples of participants from the Lifelines cohort study for cortisol and cortisone levels and associated these
to incident cardiovascular diseases (CVD) during 5 to 7 years of follow-up. We computed the odds ratio (OR) of HairGC levels for incident CVD via
logistic regression, adjusting for classical cardiovascular risk factors, and performed a sensitivity analysis in subcohorts of participants < 60 years
and ≥ 60 years of age. We also associated HairGC levels to immune parameters (total leukocytes and subtypes).
Results: Hair cortisone levels (available in n = 4701) were independently associated with incident CVD (P < .001), particularly in younger
individuals (multivariate-adjusted OR 4.21, 95% CI 1.91-9.07 per point increase in 10-log cortisone concentration [pg/mg], P < .001). All
immune parameters except eosinophils were associated with hair cortisone (all multivariate-adjusted P < .05).
Conclusion: In this large, prospective cohort study, we found that long-term cortisone levels, measured in scalp hair, represent a relevant and
significant predictor for future CVD in younger individuals. These results highlight glucocorticoid action as possible treatment target for CVD
prevention, where hair glucocorticoid measurements could help identify individuals that may benefit from such treatments.
Key Words: hair glucocorticoids, cardiovascular diseases
Abbreviations: BMI, body mass index; CVD, cardiovascular disease; HairGC, glucocorticoid levels in scalp hair; HPA, hypothalamic-pituitary-adrenal;
IL-, interleukin; IQR, interquartile range; OR, odds ratio; WC, waist circumference.
The hypothalamic-pituitary-adrenal (HPA)-axis is crucial for that is present in nearly all cells of the human body (1). Cortisol
homeostasis and survival, as its end products, glucocorticoid can be converted to the inactive hormone cortisone, and vice ver
hormones, have pleiotropic effects throughout the body and sa, by 11-beta-hydroxysteroid dehydrogenase, which is one of
are important regulators of metabolism and immune function. the mechanisms that regulates the net glucocorticoid effect at tis
Upon physiological or psychological stress, the HPA-axis is acti sue level (2).
vated starting with release of corticotropin-releasing hormone Pathological disturbances in the amounts of circulating
(CRH) and adrenocorticotropic hormone (ACTH) from the glucocorticoids cause a variety of serious clinical symptoms.
hypothalamus and anterior pituitary, respectively. This subse In case of severe glucocorticoid excess, known as the
quently leads to increasing levels of circulating glucocorticoid Cushing syndrome, patients develop features resembling the
hormones. metabolic syndrome, including abdominal obesity, insulin re
In humans, the main active glucocorticoid hormone is cortisol. sistance, and hypertension (3). It has been hypothesized that
Cortisol binds to the glucocorticoid receptor, a nuclear receptor also in the general population increased glucocorticoid action
Received: 18 August 2023. Editorial Decision: 6 February 2024. Corrected and Typeset: 19 March 2024
© The Author(s) 2024. Published by Oxford University Press on behalf of the Endocrine Society.
This is an Open Access article distributed under the terms of the Creative Commons Attribution License (https://creativecommons.org/licenses/by/4.0/), which
permits unrestricted reuse, distribution, and reproduction in any medium, provided the original work is properly cited.
The Journal of Clinical Endocrinology & Metabolism, 2024, Vol. 109, No. 10 2521
contributes to metabolic derangements, perhaps to a lesser Manenschijn et al (11)). With alpha set at 0.05 and 1-β set at
extent (4). However, historically it has been difficult to estab 0.90, an estimated 1056 participants were required for each
lish a clear relation between the levels of circulating gluco quartile. Accounting for an expected response rate of 85%
corticoid and anthropometric measurements in the general (known from the Lifelines second visit) and for anticipated tech
population, and conflicting results were found when gluco nical difficulties in the HairGC measurements, we estimated
corticoid measurements as in blood, saliva, or urine were re that we required 6000 to 6500 hair samples for the current
lated to, for example, body mass index (BMI) and waist study.
circumference (WC) (5). However, all these measurements To enable future associations with genetic data, we selected
concern short-term or point measurements and may not the 6341 participants from whom hair samples from 2014 as
represent long-term exposure to glucocorticoids. An interest well as genome-wide association study (GWAS) data were avail
ing perspective may come from hair glucocorticoid (HairGC) able. The follow-up visit was scheduled to take place between
analyses, including the measurements of both cortisol and its October 2019 and April 2020, but eventually this was extended,
inactive form cortisone in proximal scalp hair. This is a reli as site visits were postponed due to COVID-19 measures. For
n 4758 1571
a
Hair cortisone levels, pg/mg, median [IQR] 5.19 [4.02, 7.09] 6.72 [5.15, 9.59] <.001 25.9
Hair cortisol levels (%)a < 2.5 pg/mg 1619 (76.4) 461 (70.9) .005 56.2
> 2.5 pg/mg 499 (23.6) 189 (29.1)
Age, years, mean (SD) 53.21 (10.06) 54.29 (10.88) <.001 0.2
BMI (kg/m2), mean (SD) 26.31 (4.62) 26.46 (3.64) .239 0.3
Waist circumference (cm), mean (SD) 87.94 (12.13) 96.10 (10.80) <.001 0.3
Corticosteroid use in past 3 months (%) no 3505 (80.6) 1197 (85.4) <.001 9.3
yes 846 (19.4) 204 (14.6)
Parametric data are presented as mean(SD), nonparametric data are presented as median (IQR; interquartile range).
Abbreviations: BMI, body mass index; HDL, high-density lipoprotein; LDL, low-density lipoprotein.
a
Flowchart of hair cortisol and cortisone results is given in Fig. 1A and 1B.
Other Relevant Covariates analyses (> 2.5 pg/mg or < 2.5 pg/mg). Descriptive statistics
By using questionnaires, we obtained self-reported data concern of the baseline measurements are reported as median and
ing education level (classified as low/medium/high according to interquartile range (IQR) for nonparametric data and mean
the standardized cutoffs applicable to the Dutch education sys ± SD for parametric data. We investigated via linear regres
tem, Standaard onderwijsindeling); alcohol use based on ques sion whether anthropometric data at baseline were associated
tions regarding the number of drinking days per week and the with hair cortisone and cortisol levels.
number of drinks per occasion (subsequently classified as First, we used univariate logistic regression to estimate the
“Non-drinkers,” “0-1 drinks per day,” “1-2 drinks per day,” influence of HairGC levels and classical cardiovascular risk
or “More than 2 drinks per day”); and past or current smoking factors on the outcome incident CVD (present/absent).
(yes/no). Participants were also asked whether they used any Next, we used stepwise logistic regression to estimate the in
corticosteroid-containing drugs in the last 3 months before tak fluence of HairGC levels on incident CVD. First, age was
ing the hair sample (yes/no) and whether they had been diag added to the model, including the interaction term
nosed with diabetes on any of the previous visits (yes/no). age*HairGC level (Model A). To investigate sex-differences,
we performed exploratory analysis in males and females sep
arately, but refrained from further sex-specific analyses due to
Statistical Analysis a lack of statistical power, particularly for males. Also, insuf
Analyses were performed using R Studio (Version 4.2.2). ficient data were available on menopausal status of females.
HairGC results below
√ the lower limit of quantification were The final analyses were corrected for age, including the
replaced by 2.5/ 2 (25), and because of non-normality they interaction term age*HairGC, sex, corticosteroid use, WC,
were 10-log transformed for all analyses. As the majority of current smoking, systolic blood pressure, high-density lipo
cortisol data were below the lower limit of quantification, protein (HDL)- and total cholesterol, having a past history
we used cortisol as a dichotomous variable in longitudinal of diabetes mellitus, and having a past history of CVD
The Journal of Clinical Endocrinology & Metabolism, 2024, Vol. 109, No. 10 2523
Reason exclusion n
Less than 7.5mg hair: 889
Other reasons*: 1819
Ion-rao: 220
Reason exclusion n
Signal to noise <10, co-eluon or no peak: 637
Reason exclusion n
Less than 7.5mg hair: 889
Other reasons*: 654
Ion-rao: 61
Reason exclusion n
Signal to noise <10, co-eluon or no peak: 46
Corsol <2.5pg/mg Quantave data:
n=171 n=4530
Figure 1. Flowchart of hair cortisol (A) and hair cortisone analysis (B). *Indicates problem with sample preparation, sample matrix, and/or
chromatography.
Figure 2. Heatmap showing the predicted probability of cardiovascular diseases during follow-up, according to the interaction between age at baseline
and hair cortisone levels (pg/mg, 10-log transformed) at baseline. Red color indicates high probability of new cardiovascular diseases (CVD), purple
indicates lower probability.
2524 The Journal of Clinical Endocrinology & Metabolism, 2024, Vol. 109, No. 10
Table 2. Results of the multivariate logistic regression for the and cortisone levels (adjusting for age, sex, BMI, and recent
outcome self-reported incident cardiovascular diseases during the 5 corticosteroid use), and whether they were independent pre
to 7 years of follow-up, for (a) hair cortisone (as continuous
predictor, including the interaction term hair cortisone*age) and (b)
dictors for incident CVD by adding them to Model B.
hair cortisol (dichotomized as above/below the lower limit of
quantification [2.5pg/mg])
Results
Characteristic ORa 95% CIa P value Hair Glucocorticoids and Anthropometrics at
Hair cortisone levels (total n = 2878, number of events = 214)
Baseline
Hair cortisone level (pg/mg), 83.3 2.36, 2855 .015
Baseline characteristics of the full cohort from whom hair
10-log samples were obtained (n = 6341), are shown in Table 1. In
Age, per 10 years 2.64 1.55, 4.51 <.001
brief, the population was predominantly female (75%), the
mean age was 53 years, and the mean BMI was 26 kg/m2.
Hair cortisone level (pg/mg), 0.52 0.28, 0.96 .038
10-log * Age, per 10 years Hair cortisol levels were below the lower limit of quantifica
Sex-Specific Analysis of quantification, having high hair cortisol (> 2.5 pg/mg)
In a separate univariate analysis within female participants only was in trend associated with future CVD (P = .072 for
(n = 2441, 165 CVD cases) hair cortisone levels were still signifi hair cortisol, P < .001 for age, and P = .075 for the inter
cantly associated with CVD (P < .001, OR = 3.05 [95% CI, action term cortisol*age; Table 2B). In the subgroup of par
01.56-5.79]). In males, the association lost significance ticipants with known cortisol data, cortisone was still
(n = 704, 78 CVD cases, P = .2, OR = 1.81 [95% CI, significant in univariate analysis but not in multivariate
0.65-4.91]). Due to lack of power, we refrained from further analyses according to Model A and B (n = 1847, P = .02
sex-specific analyses, but we used “sex” as a covariate in all for univariate analysis; n = 1847 and P = .11 for Model A;
multivariate models. and n = 1695 and P = .2 in Model B). In the subgroup of
participants without past CVD, hair cortisone levels were
Fully Adjusted Models still significantly associated with incident CVD (P < .019
for hair cortisone, P = .057 for hair cortisone*age when
In the fully adjusted model (using Model B, including classical
used in Model B).
cardiovascular risk factors), hair cortisone was still associated
with incident CVD (Table 2A, graphical representation of pre
dicted probability for incident CVD by hair cortisone levels Sensitivity Analysis in Age Categories
shown in Fig. 3A). Within the subgroup of individuals < 60 years, hair cortisone
Cortisol was not significantly associated with CVD in was strongly and significantly associated with incident
Model B when used on a continuous scale (P = .193 and CVD after adjustment for relevant covariates via Model B
P = .209 for the interaction term cortisol*age). However, (OR = 4.21 [95% CI, 1.91-9.07], P = .005, Table 3). In elder
when dichotomized in those below and above the lower limit ly individuals (≥ 60 years), hair cortisone levels were not
2526 The Journal of Clinical Endocrinology & Metabolism, 2024, Vol. 109, No. 10
Table 3. Results of the sensitivity analyses in subgroups based on (r = 0.025, P = .19). After adjustment for age, sex, BMI, and
age (categorized as < 60 years and ≥ 60 years) corticosteroid use, only lymphocytes were significantly associ
ated with hair cortisol (P = .04, Table 5).
Predicting incident cardiovascular diseases by baseline hair cortisone,
split by age category
Immune Parameters at Follow-Up
Characteristic OR 95% CI P value When total leukocytes, neutrophils, and monocytes were
< 60 years (total n = 2229, number of events = 125) added to Model B with cortisone as predictor, only monocytes
were a significant independent risk factor for incident CVD
Hair cortisone level (pg/mg), 10-log 4.21 1.91, 9.07 <.001
(OR = 3.37 [95% CI, 1.19-9.36], P = .021), leukocytes and
Sex, male 1.14 .71, 1.79 .6
neutrophils were not significant (Supplemental data 4A &
Waist circumference 1.01 .99, 1.02 .5 4B (24)). In the same model with dichotomized cortisol,
Any corticosteroid use in last 1.13 .68, 1.81 .6 none of the immune parameters were significant predictors
3 months, yes of incident CVD.
Total leukocytes (10 e-9/L) Neutrophil granulocytes Monocytes (10 e-9/L) Lymphocytes (10 e-9/L) Eosinophil granulocytes
(10 e-9/L) (10 e-9/L)
Hair cortisone levels (pg/mg), 10-log 0.48 (0.24-0.71) <.001 0.44 (0.28-0.60) <.001 0.04 (0.02-0.06) <.001 0.09 (0.00-0.17) .046 0.00 (−0.02-0.02) .781
Age, per 10 years −0.15 (−0.20-−0.10) <.001 −0.14 (−0.18-−0.11) <.001 0.01 (0.00-0.01) <.001 −0.01 (−0.03-0.01) .363 0.00 (−0.00-0.01 .483
Sex, male −0.01 (−0.14-0.11) .821 −0.05 (−0.14-0.04) .260 0.05 (0.04-0.07) <.001 −0.07 (−0.12-−0.02) .003 0.03 (0.02-0.03) <.001
Recent corticosteroid use, yes 0.27 (0.14-0.40) <.001 0.19 (0.09-0.28) <.001 0.02 (0.01-0.03) <.001 0.04 (−0.01-0.09) .127 0.03 (0.02-0.04) <.001
BMI (kg/m2) 0.07 (0.06-0.08) <.001 0.05 (0.04-0.06) <.001 0.00 (0.00-0.00) <.001 0.02 (0.01-0.02) <.001 0.00 (−0.00-0.00) .064
Observations 4243 4210 4210 4210 4210
Results of a multivariate linear regression model where hair cortisone levels (pg/mg), age, sex, recent corticosteroid use and body mass index (BMI) are included. Linear regression estimates and 95% confidence intervals are provided.
P values were significant at P < .05 (presented in bold).
The Journal of Clinical Endocrinology & Metabolism, 2024, Vol. 109, No. 10
Table 5. Cross-sectional linear regression models investigating the associations between hair cortisol levels and immune parameters, adjusting for age, sex, recent corticosteroid use and body mass
index (BMI, kg/m2)
Predictors Total Leukocytes (10 e-9/L) Neutrophil Granulocytes Monocyte (10 e-9/L) Lymphocytes (10 e-9/L) Eosinophil Granulocytes
(10 e-9/L) (10 e-9/L)
Hair cortisol levels (pg/mg), 10-log 0.36 (−0.04-0.76) .075 0.22 (−0.05-0.50) .112 0.01 (−0.02-0.04) .471 0.15 (0.01-0.30) .036 0.02 (−0.01-0.05) .158
Age, per 10 years −0.18 (−0.24-−0.11) <.001 −0.15 (−0.19-−0.10) <.001 0.01 (0.00-0.01) .001 −0.02 (−0.05-0.00) .058 0.00 (−0.00-0.01) .849
Sex, male 0.06 (−0.10-0.22) .499 0.02 (−0.09-0.13) .659 0.06 (0.05-0.07) <.001 −0.08 (−0.14-−0.03) .004 0.03 (0.01-0.04) <.001
Recent corticosteroid use, yes 0.29 (0.12-0.47) .001 0.21 (0.09-0.33) .001 0.02 (0.01-0.03) .005 0.03 (−0.03-0.10) .333 0.03 (0.02-0.04) <.001
2
BMI (kg/m ) 0.07 (0.06-0.09) <.001 0.05 (0.04-0.06) <.001 0.00 (0.00-0.00) <.001 0.02 (0.01-0.02) <.001 0.00 (0.00-0.00) .030
Observations 2514 2499 2499 2499 2499
Results of a multivariate linear regression model where hair cortisol levels (pg/mg), age, sex, recent corticosteroid use and body mass index (BMI) are included. Linear regression estimates and 95% confidence intervals are provided.
P values were significant at P < .05 (presented in bold).
2527
traditional risk factors. Of note, it is well established that found), or following psychological or physical stressors (eg,
mean daily glucocorticoid levels and hair glucocorticoid levels chronic pain, sleep deficit, endocrine disruptors) (4), and/or
generally increase with aging (10, 29). In this regard, it can be altered cortisol metabolism eg, in the liver or visceral adipose
hypothesized that higher glucocorticoid levels in younger indi tissue (35).
viduals may actually reflect underlying accelerated biological In Cushing disease, treating the hypercortisolism greatly re
aging, where these individuals resemble older individuals duces cardiovascular risk (36). However, until now, in people
with regard to CVD risk. This may be mediated by an acceler without Cushing disease there were only modest effects of drugs
ation of the age-related increase in activation of the immune targeting glucocorticoid action (37-39). It can be speculated that
system (“inflammaging”), which is highly relevant in the de HairGC measurements can be used to identify individuals at
velopment of atherosclerosis (30). Indeed, in the current risk who may benefit from pharmacological or non-
work we found significant cross-sectional associations be pharmacological glucocorticoid-lowering interventions or pre
tween various cellular immune parameters and glucocorticoid ventive strategies, ultimately decreasing cardiovascular risk.
levels. In the multivariate-adjusted longitudinal analyses, In summary, we demonstrate that long-term glucocorticoid
4. van der Valk ES, Savas M, van Rossum EFC. Stress and obesity: are 23. Odijk R, van den Berg S, Dorst K, de Rijke Y. An optimized LC-MS/
there more susceptible individuals? Curr Obes Rep. 2018;7(2): MS method for cortisol and cortisone in human hair. presented at:
193-203. Psychoneuroendocrinology 2019; Session 107(S) (2019) 1–81.
5. Rodriguez AC I, Epel ES, White ML, Standen EC, Seckl JR, Tomiyama 24. van der Valk ES, Mohseni M, Iyer AM, et al. Data from:
AJ. Hypothalamic-pituitary-adrenal axis dysregulation and cortisol Supplementary data to article “Long-term glucocorticoid exposure
activity in obesity: a systematic review. Psychoneuroendocrinology. and incident cardiovascular diseases—the Lifelines cohort”. Figshare.
Dec 2015;62:301-318. 2023. doi:10.6084/m9.figshare.24595611
6. Noppe G, Rijke YB, Dorst K, Akker ELT. LC-MS/MS-based meth 25. Harel O, Perkins N, Schisterman EF. The use of multiple imput
od for long-term steroid profiling in human scalp hair. Clin ation for data subject to limits of detection. Sri Lankan J Appl
Endocrinol (Oxf). 2015;83(2):162-166. Stat. 2014;5(4):227-246.
7. Short SJ, Stalder T, Marceau K, et al. Correspondence between hair 26. Crawford AA, Soderberg S, Kirschbaum C, et al. Morning plasma
cortisol concentrations and 30-day integrated daily salivary and cortisol as a cardiovascular risk factor: findings from prospective
weekly urinary cortisol measures. Psychoneuroendocrinology. cohort and Mendelian randomization studies. Eur J Endocrinol.
2016;71:12-18. 2019;181(4):429-438.
8. Mirzaian M, van Zundert SKM, Schilleman WF, et al. 27. van Bussel EF, Hoevenaar-Blom MP, Poortvliet RKE, et al.