Relationship of Serum Vitamin D Concentrations and Allostatic Load as a Measure of Cumulative Biological Risk among the US Population: A Cross-Sectional Study

Abstract

Introduction

The allostatic load (AL) index is a multi-systemic measure of physiologic dysregulation known to be associated with chronic exposure to stress and adverse health outcomes. We examined the relationship between AL and serum 25-hydroxyvitamin D (25(OH)D) concentration in non-institutionalized US adults.

Methods

Data from the Third National Health and Nutrition Examination Survey (NHANES III, 1988–94) were used to calculate two versions of AL including 9 biomarkers and another two with 14 biomarkers (systolic and diastolic blood pressure, pulse rate, serum cholesterol, serum HDL-cholesterol, glycated hemoglobin, sex-specific waist-to-hip ratio, serum albumin, and serum C-reactive protein for AL1, and, additionally body mass index, serum triglyceride, serum creatinine, and serum herpes I & II antibodies for AL2), each set defined by predefined cut-offs or by quartiles. Serum vitamin D concentration was ranked into quartiles. Logistic regression, Poisson regression and linear regression were used to examine the association of serum 25(OH)D concentrations on AL, after adjusting for biological, physiological, socioeconomic, lifestyle, and health variables.

Results

Odds Ratios (OR) for high AL of the lowest 25(OH)D serum quartile were between 1.45 (95% CI: 1.28, 1.67) and 1.79 (95% CI: 1.39, 2.32) for the fully adjusted model, depending on AL version. Inverse relationships between vitamin D serum concentrations were observed for all AL versions and every adjustment. This relationship was consistent after stratification by sex, age or ethnic background. Sensitivity to low 25(OH)D concentrations was highest among the youngest group (20–39 years) with an OR of 2.11 (95% CI: 1.63, 2.73) for the lowest vitamin D quartile Q1.

Conclusions

Vitamin D had a consistent and statistically significant inverse association with all tested models of high AL, which remained consistent after adjusting for biological, socioeconomic, lifestyle and health variables. Our study adds evidence linking low 25(OH)D concentrations with poorer health, further-reaching than bone health.

About the Authors:

Regina Frei

Affiliation: Division of Chronic Disease Epidemiology, Epidemiology, Biostatistics and Prevention Institute (EBPI), University of Zurich, Zurich, Switzerland

Sarah R. Haile

Affiliation: Division of Chronic Disease Epidemiology, Epidemiology, Biostatistics and Prevention Institute (EBPI), University of Zurich, Zurich, Switzerland

Margot Mutsch

* E-mail: [email protected]

¶‡ MM and SR are joint last authors on this work.

Affiliation: Division of Chronic Disease Epidemiology, Epidemiology, Biostatistics and Prevention Institute (EBPI), University of Zurich, Zurich, Switzerland

Sabine Rohrmann

¶‡ MM and SR are joint last authors on this work.

Affiliation: Division of Chronic Disease Epidemiology, Epidemiology, Biostatistics and Prevention Institute (EBPI), University of Zurich, Zurich, Switzerland

Introduction

Chronic stress, such as repeated exposure to stressors–or exposure over an extended time period–is known to be associated with negative effects on human health [1]. Thereby, chronic stress has been linked with health conditions such as cardiovascular disease (CVD), metabolic, immunologic, psychological and cognitive disorders [2–5]. The universal use of the term “stress” in popular language makes it an imprecise word to describe how the body copes with changing environmental, psychosocial and physical circumstances. In order to find a more comprehensive and process-oriented term, “allostasis” was established, meaning literally “maintaining stability through change”. The allostasis model suggests that the goal of regulation is resilience to adapt to environmental demands [6]. To match physiological parameters of chronic stress with allostasis, the allostatic load (AL) index was introduced by McEwen and Stellar [1]. AL is a cumulative measure of biomarkers representing physiological stress in the organism, such as neuroendocrine, cardiovascular, immune and metabolic biomarkers [7]. High AL has been linked with several adverse health outcomes including all-cause mortality [7–9], cardiovascular diseases [7, 10], and diminished mental health [11].

The AL index was originally measured using 10 biomarkers: systolic and diastolic blood pressure, total cholesterol (TC), high-density-lipoprotein (HDL), glycated hemoglobin (HbA1c), waist-to-hip ratio (WHR), dehydroepiandrosterone sulfate (DHEA-S), urinary epinephrine, norepinephrine, and cortisol [7]. The biomarkers were dichotomized with the highest quartile coding for high risk (lowest quartile for DHEA-S and HDL) and summed for the AL index. More recent studies used variable numbers of biomarkers and different approaches to dichotomization [12, 13]. Currently, there is no consensus on either the choice of or method for dichotomizing the biomarkers which comprise the AL index.

Various studies have linked AL with social contexts, including socioeconomic status [14, 15], ethnicity [8, 16, 17], neighborhood [17, 18], gender [19], and social relationships [20]. Rosenberg et al. found an association between AL and serum carotenoid concentration among middle-aged US adults [21]. However, there is limited knowledge on the impact of lifestyle factors, especially nutrition, on AL.

Discussion of vitamin D has expanded considerably since studies have found its role far beyond bone health. Recent studies have focused on the impact of 25(OH)D on the immune response, respiratory health, diabetes risk, cancer, CVDs or obesity, as presented in several systematic reviews and meta-analyses [22–24].

According to Troesch et al., who analyzed the vitamin intake in several industrialized countries (including the US, Germany, UK, and the Netherlands), vitamin D is one of the critical vitamins for which intake is often below the recommended level [25]. Globally, 6.7% of the overall population is reported to be vitamin D deficient (concentration below 25 nmol/L), 37% have estimates below 50 nmol/L [26].

After vitamin D is photosynthesized in the skin or ingested via food, it is hydroxylated in the liver to 25(OH)D (calcifediol) and then transported to the kidneys, where it is hydroxylated to 1,25(OH)2D (calcitriol), the active metabolite [27]. 1,25(OH)2D is a potent steroid hormone that regulates the immune system via the vitamin D receptor (VDR) which is present in many immune cell types such as monocytes, macrophages and dendritic cells [28]. In general, 1,25(OH)2D activates the innate immune system and inhibits adaptive immunity, and as a result it may serve a critical role in coordinating inflammatory and anti-inflammatory processes [29, 30]. Due to the interrelation of chronic stress, immune response, and inflammation, it has been suggested that sufficient concentrations of vitamin D may protect against chronic stress, assessed with a high AL index as a marker. In this cross-sectional study, we therefore aimed to investigate the relationship between AL and serum 25(OH)D concentrations. Based on our findings, we discuss the implications for the reduction of stress-related adverse health outcomes.

Methods

Allostatic Load

Allostatic load was calculated as a summary measure of multisystem dysregulation including cardiovascular, metabolic and inflammation-related biomarkers. We defined two AL scores: allostatic load 1 (AL1) calculated from 9 biomarkers and allostatic load 2 (AL2) calculated from 14 biomarkers. Because different methods for low/high risk calculation of biomarkers were used in previous studies, we assessed our biomarkers in two ways, once using predefined cut points and once analyzing the AL counts in quartiles (Table 1). This leads to four AL versions for this study: AL1 and AL2 calculated with the cut point method and AL1 and AL2 calculated with the quartile method.

[Figure omitted. See PDF.]

Table 1. Cut points for biomarkers and prevalence of high-risk biomarkers with cut point and quartile method.

https://doi.org/10.1371/journal.pone.0139217.t001

In AL1, we included systolic and diastolic blood pressure, pulse rate, serum cholesterol, serum HDL, glycated hemoglobin, WHR, serum albumin, and serum C-reactive protein (CRP) as it was done previously [14, 15]. For AL2, we additionally included body mass index (BMI) [31], serum triglyceride (TG) [32], serum creatinine [33], and serum herpes simplex virus I & II antibodies (H1A & H2A). Fig 1 provides an overview. Each of the biomarkers was dichotomized into “high risk = 1” and “low risk = 0”. The dichotomization was done with cut points, as shown in Table 1, or with quartiles, respectively. For the quartiles, the highest quartile coded for high risk for all biomarkers except albumin and HDL-cholesterol, where the lowest quartile represented high risk. Serum herpes I & II antibodies were binary variables with positive indicating as “high risk”. AL score-ranges were 0–9 for AL1 and 0–14 for AL2. For all four versions of AL, summary scores of 0 to 2 were defined as “low AL” and ≥3 as “high AL”, as done in previous studies [8, 16, 20].

[Figure omitted. See PDF.]

Fig 1. Procedure to establish the four allostatic load versions.

AL, allostatic load index, BMI, body mass index, HDL-cholesterol, high-density lipoprotein-cholesterol.

https://doi.org/10.1371/journal.pone.0139217.g001

Study Population

The data for the current study originates from the Third National Health and Nutrition Examination Survey (NHANES III, 1988–1994). NHANES is a cross-sectional study providing information on the distribution of various health and nutritional indicators and potential risk factors. The survey is representative of the non-institutionalized US population, with oversampling of African-Americans, Mexican-Americans and people aged 60 years and over. Additional information about NHANES III can be found elsewhere [34, 35].

The original sample size for this study was 20,050. Participants younger than 20 years (n = 1,225) or without valid serum vitamin D data (n = 3,076) were excluded. Furthermore, participants with missing data for any of the biological markers used for the calculation of the AL (n = 4,844 for AL1 and n = 15,057 for AL2) were excluded. After missing data exclusions, data from 14,213 subjects, (7,510 females) were available for calculating AL1 and 4,620 subjects (2,334 females) available for AL2. Some participants have more than one missing value(for example one is under 20 years old and has a missing value for 25(OH)D shows up in both "missing numbers"), and therefore subtraction leads to other results.

Serum Vitamin D 25(OH)D

Concentrations of serum vitamin D 25(OH)D were measured with a radioimmunoassay kit (DiaSorin, Inc., Stillwater, MN; normal range 22.5–94 nmol/L). Blood samples were collected in northern latitudes during summer and in southern latitudes during winter.

In the absence of a generally accepted definition of 25(OH)D deficiency, most investigators assessed serum concentrations higher than 50 nmol/L as sufficient [36]. Serum concentrations under 25 nmol/L have been associated with disorders of bone metabolism [26, 36]. Therefore, we used 25 nmol/L as a cut point for severe vitamin D deficiency and 50 nmol/L for moderate vitamin D deficiency. Furthermore, vitamin D serum concentrations were categorized into quartiles based on the distribution in the cohort with quartile 1 (Q1) as the lowest vitamin D concentration. We defined vitamin D supplementation as consumption of any vitamin D containing supplement for at least one month on at least 20 days per month.

Covariates

Covariates included biological variables (age, sex, race/ethnicity), socioeconomic variables (education, census region, urbanization, marital status, poverty-income ratio), lifestyle factors (alcohol consumption, smoking status, physical activity, diet) and health variables (self-reported general health, long-term medication). Age in years was included as a continuous factor. Race/ethnicity categories were non-Hispanic White, non-Hispanic Black, Mexican-American and other. Education was classified based on US Census criteria in “less than high-school”, “high-school graduate” and “college graduate”. Living environment was defined both geographically (Northeast, Midwest, South and West) and by type (i.e., as urban (central counties of metro areas with a population of 1 million or more) or rural). Marital status was coded as “married or living as married”, “widowed, divorced or separated” and “never married”. Poverty-income ratio is the ratio of the midpoint of the observed family income as the numerator and the poverty threshold, the age of the family reference person, and the calendar year in which the family was interviewed as the denominator. According to the NHANES III recommendation, poverty-income ratios below 1 were coded as “below poverty” and 1 and above as “at or above poverty” [35]. Alcohol consumption was based on participants’ statements about frequency of beer, wine or hard liquor intake, and was grouped into sex-adjusted categories of “non-drinker”, “moderate-drinker” or “excessive-drinker”, as described elsewhere [37]. The healthy-eating index is a measure of overall dietary quality and was treated as a continuous variable, as described elsewhere [38]. Smoking behavior was determined by questionnaire and grouped into “non-smoker”, “current-smoker” or “ex-smoker”. The frequency of moderate and vigorous physical activity was grouped into age-adjusted categories of “no activities”, “1–2 activities per week”, “3–4 activities per week” and “5+ activities per week”. Self-reported general health was determined by questionnaire and included “excellent”, “very good and good” and “fair and poor”. Long-term medications, including systemic drugs such as antihypertensives, antihyperlipidemics, blood glucose regulators, mental health medications, antineoplastics and others, were counted (“medication”) if they had been taken for more than 90 days, otherwise it was considered “no medication”.

Statistical Analyses

The current study results were weighted for the complex multistage probability sampling design and population estimation, with weighting factors provided by NHANES III, in order to account for the complex survey design, survey non-response, and post-stratification. Logistic regression was used to calculate the relationship of AL (high versus low) with individual serum vitamin D. AL could also be considered a count of the corresponding biomarkers, which was analyzed using Poisson regression. The results from zero-inflated Poisson regression, allowing more frequent 0 counts than otherwise expected, were similar to those of the standard Poisson regression, and are not shown. A further analysis considered AL as a continuous outcome examined using linear regression (S1 Table). Vitamin D was considered both as quartiles and as a continuous predictor. A basic model without adjustments was performed; further models adjusted in a nested fashion for biological factors, socioeconomic factors, lifestyle factors, and the final model additionally adjusted for health factors. All models were assessed for AL1 and AL2, both for the cut point definition and for the quartile definition. The association of vitamin D supplementation and AL as well as serum vitamin D concentration was examined using the χ2 test. Correlations between vitamin D and the individual biomarkers were estimated using the square root of the R² statistic from the corresponding linear regression models. Analyses were performed using IBM SPSS Statistics Version 22. Linear and (zero inflated) Poisson regression analyses were performed using R (version 3.2.1).

Results

Vitamin D serum concentrations ranged from 8.7 to 243.6 nmol/L with a mean concentration of 73.7 nmol/L. The mean among male participants was 77.6 nmol/L, females had a significantly lower mean concentration of 70.2 nmol/L (p<0.001). In our sample, 2.0% of all participants had a severe deficiency in 25(OH)D concentration and 22.5% a moderate deficiency.

AL scores ranged from 0–9 for AL1 with mean values of 1.8 (standard deviation SD: 1.44) for the cut point method, 1.9 (SD: 1.68) for the quartile method. Due to the additional biomarkers, AL scores ranged from 0–14 for AL2, with mean values of 3.0 (SD: 2.06) for the AL2 cut point method and 3.3 (SD: 2.30) for the quartile method, respectively. The prevalence of a high AL was 28.6% (32.5%) for AL1 cut point (quartile) and expectedly higher at 52.8% (57.2%) for AL2 cut point (quartile).

The high-risk prevalence of individual biomarkers according to the categorization method differed for some biomarkers between the methods used (Table 1). About 25% of participants (deviations due to weighting) were in the high-risk group of glycated hemoglobin, diastolic blood pressure, and serum creatinine, but only 5.4%, 6.4%, and 5.5%, respectively, based on the cut point method. On the other hand, 62.7% were in the high-risk group for WHR when using sex-specific cut points compared with approximately 25% using the quartile method. The highest correlations between vitamin D and individual biomarkers were found for glycated hemoglobin, albumin, body mass index and antibodies against herpes simplex viruses (Table 1).

A statistically significant inverse relationship was observed between vitamin D supplementation and 25(OH)D deficiency (p<0.001, chi-square test). However, there was no statistically significant association between vitamin D supplementation and AL1 cut point (p = 0.121) or AL1 quartile (p = 0.050).

Table 2 shows the baseline characteristics by AL (AL1 cut point) and vitamin D quartiles. The distribution towards the highest quartile for both sexes was caused by the NHANES III weight, which shifted more than 25% of participants into the highest vitamin D quartile. Ethnic background was strongly associated with vitamin D quartiles, as only 9.3% non-Hispanic White participants were in the lowest quartile, but 49.4% of non-Hispanic Blacks and 22.3% of Mexican-Americans. Higher education, urban living environment and better economic status were associated with lower AL and higher serum 25(OH)D concentrations. Non-smokers, moderate-drinkers and physically active participants had the lowest percentage of high AL, whereas ex-smokers, moderate drinkers and physically active participants had the highest percentage in the highest vitamin D quartile. Participants considering themselves having excellent general health and those who took no longterm medication had lower AL and higher 25(OH)D serum concentrations.

[Figure omitted. See PDF.]

Table 2. Baseline characteristics by AL1 cut point and 25(OH)D concentration quartiles.

https://doi.org/10.1371/journal.pone.0139217.t002

Vitamin D concentrations were significantly associated with AL (Table 3). For all four AL versions tested and for all adjustments, the odds ratios for high AL were statistically significantly higher for participants in the lowest vitamin D quartile compared to the highest quartile. In the fully adjusted model, the ORs of having high AL for those in the lowest 25(OH)D quartile were 1.45 (for AL1 quartiles: 95% CI: 1.28,1.67) and 1.79 (for AL2 cutpoints: 95% CI: 1.39, 2.32) depending on AL definition.

[Figure omitted. See PDF.]

Table 3. Regression results for AL and 25(OH)D concentration.

https://doi.org/10.1371/journal.pone.0139217.t003

Considering vitamin D concentration as a continuous variable showed similar results. In all cases, increasing vitamin D by 10 nmol/L was slightly but statistically significantly associated with a decrease in allostatic load (Table 3), with ORs ranging from 0.91 (for AL2 cutpoint: 95% CI: 0.88, 0.93) to 0.95 (for AL1 cutpoint: 95% CI: 0.94, 0.97; for AL1 quartile, 95% CI: 0.93, 0.98).

Examining allostatic load as the count of the corresponding biomarkers, subjects with the lowest quartile vitamin D had up to 1.23 times the AL score as those in the highest quartile (rate ratio for AL2 cutpoint, 95% CI: 1.15–1.31) (Table 3). Further considering vitamin D as a continuous covariate, vitamin D continues to show a slight but statistically significant negative association with allostatic load, with a rate ratio for AL1 cutpoint of 0.98 (95% CI: 0.98–0.99). Results from the linear regression analysis of allostatic load are found in the supplementary material (S1 Table).

Stratifyingthe analysis by sex showed no differences between male and female participants (Fig 2). No cut-off age that altered the relationship of vitamin D with AL was observed when logistic regression was performed by age categories, although sensitivity to low 25(OH)D concentrations was highest among the youngest group (20–39 years) with an OR of 2.11 (95% CI: 1.63–2.73) for the lowest vitamin D quartile Q1. Higher serum 25(OH)D was associated with lower AL among all ethnicities, but only non-Hispanic Whites showed significant ORs for all vitamin D quartiles compared to the highest quartile, with an OR of 1.46 for Q1 (95% CI: 1.22–1.74). In all other ethnic groups, the comparison between top versus bottom quartile was statistically insignificant.

[Figure omitted. See PDF.]

Fig 2. Logistic regression analysis of the allostatic load index and 25(OH)D concentrations stratified by biological covariates.

Results are stratified by sex, age, and ethnic background and are presented as ORs for high AL with 95% confidence intervals. Q, 25(OH)D concentrations in quartiles with Q1 lowest and Q4 highest quartile (reference).

https://doi.org/10.1371/journal.pone.0139217.g002

Among participants describing their general health as “excellent”, we observed no statistically significant association between low 25(OH)D concentration and low AL, with an OR of 0.66 (95% CI: 0.42; 1.04) for Q1. However, participants considering their general health as “very good and good” or “fair and poor”showed significant associations between low 25(OH)D concentrations and high AL.

When logistic regression was performed with long-term medication, this phenomenon did not recur as participants with no medication had an OR of 1.51 (95%CI: 1.26; 1.80) for Q1 (Fig 3). Adding single medication groups (such as antihypertensives, antihyperlipidemics, blood glucose regulators, mental health medications, and antineoplastics) instead of the combination of all of them did not change the results significantly. Also, adding vitamin D supplementation to the fully adjusted logistic regression model did not change ORs significantly (results not shown).

[Figure omitted. See PDF.]

Fig 3. Logistic regression analysis of the allostatic load index and 25(OH)D concentrations stratified by health covariates.

Results are stratified by self-reported health status, and long-term medication and are presented as ORs for high AL with 95% confidence intervals. Q, 25(OH)D concentrations in quartiles with Q1 lowest and Q4 highest quartile (reference).

https://doi.org/10.1371/journal.pone.0139217.g003

Discussion

In this cross-sectional study, we observed a statistically significant and consistent inverse relationship between serum 25(OH)D concentration and AL as a measure of cumulative biological risk. Various ways of assessing AL provided consistent results. Those results remained similar when the analyses were adjusted for biological, socioeconomic status, lifestyle and general health status. This is the first report examining the association of serum vitamin D concentrations and AL.

The observed relationship could be explained by the various functions that vitamin D has on immune homeostasis, and there is increasing evidence that higher serum 25(OH)D concentrations can alleviate inflammatory diseases. Our findings are consistent with the results of other studies showing an inverse association of serum 25(OH)D concentrations with health events or biomarkers we had included in our AL index, including obesity [39], hypertension [40, 41], insulin resistance [42–44], CRP concentrations [45], and dyslipidemia [46]. Weyland et al. reviewed the evidence to support a causal association between vitamin D status and cardiovascular disease risk according to Hill’s criteria and found that all relevant criteria for causal association were fulfilled [47]. However, an umbrella review of hundreds of studies showed that results differ between observational and intervention studies such, that the causal effect of vitamin D on mortality remains unclear [22]. The same conclusion needs to be drawn for e.g. cardiovascular diseases, cancer, or diabetes mellitus.

The association of vitamin D and AL remained consistent after stratification by several biological and health status variables. That females had a significantly lower mean 25(OH)D serum concentration than males did not affect the association between vitamin D and AL.

Vulnerability to high AL has been linked to minority ethnic status by many authors [16, 48], which we reproduced in our study (Table 2). We observed significantly higher 25(OH)D concentrations in non-Hispanic Whites (79.1 nmol/L) compared to non-Hispanic Blacks (48.2 nmol/L) and Mexican-Americans (62.6 nmol/L). This is mostly due to different effectiveness of vitamin D synthesis depending on skin pigmentation, but other factors, such as distinct outdoor activities or dietary patterns may contribute to these differences as well. However, we only observed statistically significant associations for non-Hispanic Whites (Fig 2).

Within our sample, the prevalence of a high AL was 13.2% for participants aged 20–39 years, 33.8% for 40–59 year aged and as high as 54.9% for participants older than 60 years. It is well known that AL increases (non-linearly) with age, which is plausible because AL measures the cumulative biological risk that increases with age [49]. Sensitivity to low 25(OH)D concentrations was highest among the youngest group (20–39 years) with an OR of 2.11 (95% CI: 1.63; 2.73) for the lowest vitamin D quartile Q1, but higher 25(OH)D concentrations were associated with lower OR for high AL within all age categories. By including all adult age groups in the analysis we were able to identify an otherwise neglected at-risk group, those aged 20–39 years. Further research is needed to understand the relevance of this finding for this age group.

Self-reported general health has been shown to be a reliable predictor for mortality, even when controlling for other health-related variables [50, 51] and a recent study showed an association with AL [52]. Among the different health dimensions, vitality, in particular,was important for individuals to assess their self-reported general health [53]. Interestingly, for participants who considered their health status as “excellent,” no statistically significant association between higher 25(OH)D concentrations and AL was observed (Fig 3).

When long-term medication was used as an indicator of health (instead of self-reported general health), this effect did not appear for participants that took no medication; instead we observed a significant relationship between high vitamin D concentrations and AL.

We tested four different AL versions because there are no accepted or validated standards for this practice (number and type of biomarkers to be included, method of dichotomization). The AL1 cut point is equivalent to the method used in several current studies [8, 14, 17, 21]. Quartiles instead of cut points were used for the dichotomization of the biomarkers [48], to account for an internal reference instead of external ones. To investigate the impact of additional biomarkers, we added BMI, TG, H1A, H2A and serum creatinine to both models. Thereby, we were consistent with other studies and were able to: allow for comparisons (BMI, TG), target supplemental markers of the immune system (H1A, H2A) as well as the cardio-renal system and control for a potential relationship between renal function and serum 25(OH)D concentrations [54]. We decided not to include weights for each factor as there are no standards for that.

ORs for high AL between the four tested versions were rather similar (Table 3). However, the AL2 cut point differed in terms of gradual reduction of OR from Q4 to Q1, as it showed very similar values for Q1 and Q2 and a drop-off to Q3, which were explained by TG and creatinine. Despite this, a significant difference between the cut point and the quartile risk categorization was observed when looking at individual biomarkers; in the cut point model, a higher proportion of about 25% of the participants was allocated to the high-risk category. For the specific selection of biomarkers studied here, these differences cancelled each other out when summed to the AL index. However, for another selection of biomarkers this may not necessarily be the case, and other methods of risk classification could lead to different outcomes.

In the current study neither the number of biomarkers nor the method of dichotomization changed the results in a significant way. Due to the exclusion of participants with missing values for one of the biomarkers, the sample size was reduced from 14,213 participants for AL1 to 4,620 for AL2, which made AL2 the weaker tool for the statistical analyses of subgroups. Additional modeling as a further validity measure, such as examining the vitamin D concentration as a continuous variable and analyzing the allostatic load as a count of the corresponding biomarkers, provided the same statistically significant inverse relationship between vitamin D and AL.

The lack of agreement among prior studies on the inclusion of different biomarkers makes it difficult to compare results of different studies. Two recent reviews on measuring AL identified a total of 39 and 51different biomarkers, respectively, that were used to calculate the AL index [12, 13]. Further studies are needed to examine which combination of biomarkers represents the most valid and robust aggregate predictor of adverse health outcomes. We provide some evidence to include 25(OH)D concentrations in AL scores, but further assessments are needed.

Furthermore, AL, self-reported general health and vitamin D concentrations represent relevant markers for a health assessment as they are all associated with all-cause mortality. Within this study, their relationships were shown to be statistically significant. It remains to be seen whether preventive measures to influence these parameters might have an effect on one another as well as on health events.

Up to now, evidence that vitamin D supplementation is indicated to treat or prevent several chronic diseases is inconsistent and potential mechanisms are only partially elucidated. In spite of the widespread fortification of various foods with vitamin D, supplementation and supply from sunlight exposure, the prevalence of deficiency increased between 1988–1994 and 2001–2002 and remained on the same concentration until 2005–2006 in US adults [55]. The optimal serum concentration of 25(OH)D is still under debate and no general definitions for deficiencies are available. It might be that for prevention of immune-related diseases, higher vitamin D concentrations are needed than for bone and mineral homeostasis [47, 56]. Furthermore, it is still an open question whether a low 25(OH)D serum concentration is the cause or effect of disease. In the near future, results of four new long-term and large-scale supplementation studies are expected on the effect of vitamin D on cancer, cardiovascular diseases and diabetes, which will hopefully fill some of the current knowledge gaps [57].

This study provides no knowledge about the causal relationship of serum 25(OH)D concentrations and AL due to its cross-sectional design. 25(OH)D is considered indicative of an individual’s vitamin D status, whereas 1,25(OH)₂D, the active metabolite, is homeostatically controlled in the blood [58]. Thus, only 25(OH)D, but not 1,25(OH)2D concentrations, have been taken into account and the latter could have a cross-regulating effect [30].

Our study has several strengths. The study population is nationally representative for non-institutionalized US adults, which, therefore, enhances the generalizability of the results. The inclusion of a wide variety of biological, demographic, socioeconomic status, lifestyle, and health status-associated covariates reduces the risk of confounding. To our knowledge, AL has not previously been assessed with such a broad spectrum of covariates. The study contributes to a limited body of research on the relationship of cumulative stress exposure, AL as a health risk indicator, and vitamin D status.

Conclusions

Allostatic load is a powerful tool to assess the harmful effects of chronic stress to the human body on a multisystem level. There is an urgent need of agreement on the method and the biomarkers to calculate the AL index. Vitamin D has a statistically significant inverse association with all tested models of high AL, which was persistent after adjusting for biological, socioeconomic, lifestyle and health variables. Our study adds to the increasing evidence linking low vitamin D serum concentrations to adverse health outcomes, beyond its role in bone health. Prospective studies need to address the question whether low 25(OH)D concentrations are the cause or the outcome of high AL. If vitamin D supplementation can be shown to lower AL, vitamin D could become a relevant preventive measure to reduce the harmful effects of chronic stress in humans.

Supporting Information

[Figure omitted. See PDF.]

S1 Table. Linear regression results for allostatic load (AL) and 25-hydroxy-vitamin D (vitamin D).

https://doi.org/10.1371/journal.pone.0139217.s001

(DOCX)

Author Contributions

Conceived and designed the experiments: SR, MM. Analyzed the data: RF, MM, SR, SH. Wrote the paper: RF, SH, MM, SR.

Citation: Frei R, Haile SR, Mutsch M, Rohrmann S (2015) Relationship of Serum Vitamin D Concentrations and Allostatic Load as a Measure of Cumulative Biological Risk among the US Population: A Cross-Sectional Study. PLoS ONE 10(10): e0139217. https://doi.org/10.1371/journal.pone.0139217

References

1. McEwen BS, Stellar E. Stress and the Individual. Mechanisms Leading to Disease. Arch Intern Med. 1993;153:2093–101. doi: papers2://publication/uuid/60C31CC4-454C-4295-911E-3E42E0C4E338. pmid:8379800

2. Marin M-F, Lord C, Andrews J, Juster R-P, Sindi S, Arsenault-Lapierre Gev, et al. Chronic stress, cognitive functioning and mental health. Neurobiol Learn Mem. 2011;96(4):583–95. doi: papers2://publication/doi/10.1016/j.nlm.2011.02.016. pmid:21376129

3. Tamashiro KL, Sakai RR, Shively CA, Karatsoreos IN, Reagan LP. Chronic stress, metabolism, and metabolic syndrome. Stress. 2011;14:468–74. doi: papers2://publication/doi/10.3109/10253890.2011.606341. pmid:21848434

4. Low CA, Salomon K, Matthews KA. Chronic life stress, cardiovascular reactivity, and subclinical cardiovascular disease in adolescents. Psychosom Med. 2009;71:927–31. doi: papers2://publication/doi/10.1097/PSY.0b013e3181ba18ed. pmid:19737856

5. Bose M, Olivan B, Laferrere B. Stress and obesity: the role of the hypothalamic–pituitary–adrenal axis in metabolic disease. Curr Opin Endocrinol Diabetes Obes. 2009;16(5):340–6. doi: papers2://publication/doi/10.1097/MED.0b013e32832fa137.Stress. pmid:19584720

6. Sterling P. Principles of allostasis: optimal design, predictive regulation, pathophysiology and rational therapeutics. Allostasis, Homeostasis, and the Costs of Physiological Adaptation: MIT Press; 2003. p. 17–64.

7. Seeman TE, Singer BH, Rowe JW, Horwitz RI, McEwen BS. Price of adaptation—allostatic load and its health consequences. MacArthur studies of successful aging. Arch Intern Med. 1997;157:2259–68. doi: papers2://publication/doi/10.1001/archinte.1997.00440400111013. pmid:9343003

8. Borrell LN, Dallo FJ, Nguyen N. Racial/ethnic disparities in all-cause mortality in U.S. adults: the effect of allostatic load. Public Health Rep. 2010;125(6):810–6. doi: papers2://publication/uuid/328B79FA-0445-41AE-9502-D26BA366F408. pmid:21121226

9. Gruenewald TL, Seeman TE, Ryff CD, Karlamangla AS, Singer BH. Combinations of biomarkers predictive of later life mortality. Proc Natl Acad Sci USA. 2006;103:14158–63. doi: papers2://publication/doi/10.1073/pnas.0606215103. pmid:16983099

10. Sabbah W, Watt RG, Sheiham A, Tsakos G. Effects of allostatic load on the social gradient in ischaemic heart disease and periodontal disease: evidence from the Third National Health and Nutrition Examination Survey. J Epidemiol Community Health. 2008;62:415–20. doi: papers2://publication/doi/10.1136/jech.2007.064188. pmid:18413454

11. McEwen BS. Mood disorders and allostatic load. Biol Psychiatry. 2003;54(3):200–7. doi: papers2://publication/doi/10.1016/S0006-3223(03)00177-X. pmid:12893096

12. Juster R- P, McEwen BS, Lupien SJ. Allostatic load biomarkers of chronic stress and impact on health and cognition. Neurosci Biobehav Rev. 2010;35(1):2–16. doi: papers2://publication/doi/10.1016/j.neubiorev.2009.10.002. pmid:19822172

13. Mauss D, Li J, Schmidt B, Angerer P, Jarczok MN. Measuring allostatic load in the workforce—a systematic review. Ind Health. 2014. doi: papers2://publication/uuid/B557C76A-5AE8-42D3-8254-D7041D9D9026.

14. Seeman T, Merkin SS, Crimmins E, Koretz B, Charette S, Karlamangla A. Education, income and ethnic differences in cumulative biological risk profiles in a national sample of US adults: NHANES III (1988–1994). Soc Sci Med. 2008;66:72–87. doi: papers2://publication/doi/10.1016/j.socscimed.2007.08.027. pmid:17920177

15. Levine ME, Crimmins EM. A comparison of methods for assessing mortality risk. Am J Hum Biol. 2014;26(6):768–76. doi: pmid:25088793; PubMed Central PMCID: PMC4286244.

16. Peek MK, Cutchin MP, Salinas JJ, Sheffield KM, Eschbach K, Stowe RP, et al. Allostatic load among non-Hispanic Whites, non-Hispanic Blacks, and people of Mexican origin: effects of ethnicity, nativity, and acculturation. Am J Public Health. 2010;100:940–6. doi: papers2://publication/doi/10.2105/AJPH.2007.129312. pmid:19834005

17. Merkin SS, Basurto-D a vila R, Karlamangla A, Bird CE, Lurie N, Escarce J, et al. Neighborhoods and cumulative biological risk profiles by race/ethnicity in a national sample of U.S. adults: NHANES III. Ann Epidemiol. 2009;19(3):194–201. doi: papers2://publication/doi/10.1016/j.annepidem.2008.12.006. pmid:19217002

18. Bird CE, Seeman T, Escarce JeJ, Basurto-D a vila R, Finch BK, Dubowitz , et al. Neighbourhood socioeconomic status and biological 'wear and tear' in a nationally representative sample of US adults. J Epidemiol Community Health. 2010;64:860–5. doi: papers2://publication/doi/10.1136/jech.2008.084814. pmid:19759056

19. Juster R-P, Lupien S. A sex- and gender-based analysis of allostatic load and physical complaints. Gend Med. 2012;9:511–23. doi: papers2://publication/doi/10.1016/j.genm.2012.10.008. pmid:23158957

20. Seeman T, Glei D, Goldman N, Weinstein M, Singer B, Lin YH. Social relationships and allostatic load in Taiwanese elderly and near elderly. Soc Sci Med. 2004;59:2245–57. doi: papers2://publication/doi/10.1016/j.socscimed.2004.03.027. pmid:15450701

21. Rosenberg N, Park CG, Eldeirawi K. Relationship of serum carotenoid concentrations with allostatic load as a measure of chronic stress among middle-aged adults in the USA. Public Health Nutr. 2014;(24):1–9. doi: papers2://publication/doi/10.1017/S1368980014000056.

22. Theodoratou E, Tzoulaki I, Zgaga L, Ioannidis JPA. Vitamin D and multiple health outcomes: umbrella review of systematic reviews and meta-analyses of observational studies and randomised trials. BMJ. 2014;348:2035–54. doi: papers2://publication/doi/10.1136/bmj.g2035.

23. Autier P, Boniol M, Pizot C, Mullie P. Vitamin D status and ill health: a systematic review. Lancet Diabetes Endocrinol. 2013;2(1):76–89. doi: papers2://publication/doi/10.1016/S2213-8587(13)70165-7. pmid:24622671

24. Schöttker B, Ball D, Gellert C, Brenner H. Serum 25-hydroxyvitamin D levels and overall mortality. A systematic review and meta-analysis of prospective cohort studies. Ageing Res Rev. 2013;12(2):708–18. doi: papers2://publication/doi/10.1016/j.arr.2012.02.004. pmid:22343489

25. Troesch B, Hoeft B, McBurney M, Eggersdorfer M, Weber P. Dietary surveys indicate vitamin intakes below recommendations are common in representative Western countries. Br J Nutr. 2012;108(4):692–8. doi: papers2://publication/doi/10.1017/S0007114512001808. pmid:22691229

26. Hilger J, Friedel A, Herr R, Rausch T, Roos F, Wahl Da, et al. A systematic review of vitamin D status in populations worldwide. Br J Nutr. 2014;111(1):23–45. doi: papers2://publication/doi/10.1017/S0007114513001840. pmid:23930771

27. Chen TC, Chimeh F, Lu Z, Mathieu J, Person KS, Zhang A, et al. Factors that influence the cutaneous synthesis and dietary sources of vitamin D. Arch Biochem Biophys. 2007;460:213–7. doi: papers2://publication/doi/10.1016/j.abb.2006.12.017. pmid:17254541

28. White JH. Vitamin D signaling, infectious diseases, and regulation of innate immunity. Infect Immun. 2008;76(9):3837–43. doi: papers2://publication/doi/10.1128/IAI.00353-08. pmid:18505808

29. Hayes CE, Nashold FE, Spach KM, Pedersen LB. The immunological functions of the vitamin D endocrine system. Cell Mol Biol. 2003;49(2):277–300. doi: papers2://publication/uuid/5B9A96B0-FC59-4F87-8C73-01409988CE7D. pmid:12887108

30. Mangin M, Sinha R, Fincher K. Inflammation and vitamin D: the infection connection. Inflamm Res. 2014:803–19. doi: papers2://publication/doi/10.1007/s00011-014-0755-z. pmid:25048990

31. WHO Expert Committee. Physical status: the use and interpretation of anthropometry. WHO: 1995 pp 1–463 Report No: 854. 1995.

32. National Cholesterol Education Program Expert P. Report of the National Cholesterol Education Program Expert Panel on Detection, Evaluation and Treatment of High Blood Cholesterol in Adults. JAMA. 2001;285(19):2486–97. doi: papers2://publication/uuid/504FD4D2-02BB-49AA-AC49-B1FD15CC9954. pmid:11368702

33. Hwang A-C, Peng L-N, Wen Y-W, Tsai Y-W, Chang L-C, Chiou S-T, et al. Predicting all-cause and cause-specific mortality by static and dynamic measurements of allostatic load: a 10-year population-based cohort study in taiwan. J Am Med Dir Assoc. 2014;15(7):490–6. doi: papers2://publication/doi/10.1016/j.jamda.2014.02.001. pmid:24631353

34. Anon. Plan and operation of the Third National Health and Nutrition Examination Survey, 1988–94. Series 1: programs and collection procedures. Vital Health Stat 1: 1–407. US Department of Health and Human Service, 1994.

35. Anon. Analytic and Reporting Guidelines: The Third National Health and Nutrition Examination Survey, NHANES III (1988–94). Available: http://www.cdc.gov/nchs/nhanes/nh3rrm.htm. Accessed 1996.

36. van Schoor NM, Lips P. Worldwide vitamin D status. Best Pract Res Clin Endocrinol Metab. 2011;25:671–80. doi: papers2://publication/doi/10.1016/j.beem.2011.06.007. pmid:21872807

37. Ford ES, Zhao G, Tsai J, Li C. Low-risk lifestyle behaviors and all-cause mortality: findings from the National Health and Nutrition Examination Survey III Mortality Study. Am J Public Health. 2011;101(10):1922–9. doi: papers2://publication/doi/10.2105/AJPH.2011.300167. pmid:21852630

38. Anon. Third National Health and Nutrition Examination Survey (NHANES III), 1988–94. NHANES III Healthy Eating Index Data File. 2000.

39. Bouillon R, Carmeliet G, Lieben L, Watanabe M, Perino A, Auwerx J, et al. Vitamin D and energy homeostasis-of mice and men. Nat Rev Endocrinol. 2014;10(2):79–87. doi: papers2://publication/doi/10.1038/nrendo.2013.226. pmid:24247221

40. Tamez H, Thadhani RI. Vitamin D and hypertension: an update and review. Curr Opin Nephrol Hypertens. 2012;21:492–9. doi: papers2://publication/doi/10.1097/MNH.0b013e3283557bf0. pmid:22820371

41. Pilz S, Gaksch M, O'Hartaigh Bia, Tomaschitz A, M a rz W. The role of vitamin D deficiency in cardiovascular disease: Where do we stand in 2013? Arch Toxicol. 2013;87(12):2083–103. doi: papers2://publication/doi/10.1007/s00204-013-1152-z. pmid:24173581

42. Van Belle TL, Gysemans C, Mathieu C. Vitamin D and diabetes: The odd couple. Trends Endocrinol Metab. 2013;24(11):561–8. doi: papers2://publication/doi/10.1016/j.tem.2013.07.002. pmid:23972673

43. Dong J-Y, Zhang W-G, Chen JJ, Zhang Z-L, Han S-F, Qin L-Q. Vitamin D intake and risk of type 1 diabetes: a meta-analysis of observational studies. Nutrients. 2013;5:3551–62. doi: papers2://publication/doi/10.3390/nu5093551. pmid:24036529

44. Mitri J, Muraru MD, Pittas AG. Vitamin D and type 2 diabetes: a systematic review. Eur J Clin Nutr. 2011;65:1005–15. doi: papers2://publication/doi/10.1038/ejcn.2011.118. pmid:21731035

45. Chen N, Wan Z, Han S-F, Li B-Y, Zhang Z-L, Qin L-Q. Effect of vitamin D supplementation on the level of circulating high-sensitivity C-reactive protein: a meta-analysis of randomized controlled trials. Nutrients. 2014;6(6):2206–16. doi: papers2://publication/doi/10.3390/nu6062206. pmid:24918698

46. Zhang MC, Li HX, Liu HM, Lei H, Han L, Gao M, et al. Serum vitamin D is low and inversely associated with LDL cholesterol in the Kazak ethnic population: a cross-sectional study. Med Sci Monit. 2014;20:1274–83. doi: papers2://publication/doi/10.12659/MSM.890930. pmid:25064492

47. Weyland PG, Grant WB, Howie-Esquivel J. Does Sufficient Evidence Exist to Support a Causal Association between Vitamin D Status and Cardiovascular Disease Risk? An Assessment Using Hill's Criteria for Causality. Nutrients. 2014;6(9):3403–30. doi: papers2://publication/doi/10.3390/nu6093403. pmid:25184368

48. Geronimus AT, Hicken M, Keene D, Bound J. "Weathering" and age patterns of allostatic load scores among blacks and whites in the United States. Am J Public Health. 2006;96(5):826–33. doi: papers2://publication/doi/10.2105/AJPH.2004.060749. pmid:16380565

49. Crimmins EM, Johnston M, Hayward M, Seeman T. Age differences in allostatic load: an index of physiological dysregulation. Exp Gerontol. 2003;38(7):731–4. doi: papers2://publication/uuid/DEC4B95A-D075-48F8-945B-CC44FB07D39B. pmid:12855278

50. Zhao J, Yiengprugsawan V, Seubsman S- A, Kelly M, Bain C, Sleigh A. Self-reported health and subsequent mortality: an analysis of 767 deaths from a large Thai cohort study. BMC Public Health. 2014;14(1):860. doi: papers2://publication/doi/10.1186/1471-2458-14-860.

51. Idler EL, Benyamini Y. Self-rated health and mortality: a review of twenty-seven community studies. J Health Soc Behav. 1997;38:21–37. doi: papers2://publication/doi/10.2307/2955359. pmid:9097506

52. Vie TL, Hufthammer KO, Holmen TL, Meland E, Breidablik HJ. Is self-rated health a stable and predictive factor for allostatic load in early adulthood? Findings from the Nord Trøndelag Health Study (HUNT). Soc Sci Med. 2014;117:1–9. doi: papers2://publication/doi/10.1016/j.socscimed.2014.07.019. pmid:25016460

53. Au N, Johnston DW. Self-assessed health: What does it mean and what does it hide? Soc Sci Med. 2014;121C:21–8. doi: papers2://publication/doi/10.1016/j.socscimed.2014.10.007.

54. Tak YJ, Lee JG, Song SH, Kim YJ, Lee S, Jung DW, et al. The Relationship Between the Level of Serum 25-Hydroxyvitamin D and Renal Function in Patients Without Chronic Kidney Disease: A Cross-Sectional Study. J Ren Nutr. 2014. doi: papers2://publication/doi/10.1053/j.jrn.2014.07.007.

55. Looker AC, Dawson-Hughes B, Calvo MS, Gunter EW, Sahyoun NR. Serum 25-hydroxyvitamin D status of adolescents and adults in two seasonal subpopulations from NHANES III. Bone. 2002;30(5):771–7. doi: papers2://publication/uuid/DDE3C5CF-A5B9-4369-BE44-1198CFDC9971. pmid:11996918

56. Hewison M. An update on vitamin D and human immunity. Clin Endocrinol. 2012;76(3):315–25. doi: papers2://publication/doi/10.1111/j.1365-2265.2011.04261.x.

57. Bendik I, Friedel A, Roos FF, Weber P, Eggersdorfer M. Vitamin D: a critical and essential micronutrient for human health. Front Physiol. 2014;5:248. doi: papers2://publication/doi/10.3389/fphys.2014.00248. pmid:25071593

58. Zerwekh JE. Blood biomarkers of vitamin D status. Am J Clin Nutr. 2008;87(4):1087S–91S. pmid:18400739.