Skip to main content

Association between multiple vitamins and bone mineral density: a cross-sectional and population-based study in the NHANES from 2005 to 2006

Abstract

Background

Bone mineral density (BMD) alterations in response to multivitamin exposure were rarely studied. Our study assessed the association of coexposure to six types of vitamins (i.e., vitamins B12, B9, C, D, A and E) with BMD measurements in adults in the US.

Methods

Data were collected from participants aged ≥ 20 years (n = 2757) in the U.S. National Health and Nutrition Examination Surveys (NHANES) from 2005 to 2006. Multiple linear regression, restricted cubic splines, principal component analysis (PCA) and weighted quantile sum (WQS) regression were performed for statistical analysis.

Results

The circulating levels of vitamins B12 and C were positively associated with BMDs, and an inverted L-shaped exposure relationship was observed between serum vitamin C and BMDs. PCA identified two principal components: one for ‘water-soluble vitamins’, including vitamins B12, B9 and C, and one for ‘fat-soluble vitamins’, including vitamins A, D and E. The former was positively associated with total femur (β = 0.009, 95%CI: 0.004, 0.015) and femoral neck (β = 0.007, 95%CI: 0.002, 0.013) BMDs, and the latter was negatively associated with BMDs with non-statistical significance. The WQS index constructed for the six vitamins was significantly related to total femur (β = 0.010, 95%CI: 0.001, 0.018) and femoral neck (β = 0.008, 95%CI: 0.001, 0.015) BMDs, and vitamins B12 and C weighted the most. The WQS index was inversely related to BMDs with non-statistical significance, and vitamins E and A weighted the most.

Conclusion

Our findings suggested a positive association between water-soluble vitamin coexposure and BMD, and the association was mainly driven by vitamins B12 and C. Negative association between fat-soluble vitamin coexposure and BMD was indicated, mainly driven by vitamins E and A. An inverted L-shaped exposure relationship was found between vitamin C and BMD.

Peer Review reports

Introduction

Osteoporosis is one of the most serious bone diseases affecting elderly individuals worldwide. The International Osteoporosis Foundation reported more than 8.9 million fractures due to osteoporosis every year globally [1]. Osteoporosis may cause unforeseen fractures and other distressful events, such as pain, disability and huge financial burden [2]. The mortality rate is up to 20%–40% in the first year after the occurrence of a hip fracture [3]. Many factors, including ageing, being female, disrupted estrogen metabolism and lack of physical activity, are associated with bone mineral density (BMD) reduction and contribute to osteoporosis pathogenesis [4, 5].

In recent years, the influence of diet nutrients and vitamins on bone health have been extensively investigated. A combination of vitamins D and K or supplementation of vitamin E alone inhibited bone loss in rats and mice [6, 7]. An in vitro study showed that high levels of vitamins C and E as antioxidants can promote the proliferation of bovine osteoblasts and support bone regeneration [8]. Several population-based studies revealed that elevated amount of vitamin C was associated with increased BMD and decreased risk of fracture [2, 9]. A study from the United States of America indicated that vitamins B12 and B9 were beneficial for osteoporosis prevention in the elderly women [10]. However, a recent randomized clinical trial performed on 2919 elderly individuals demonstrated no reduction in fracture risk after daily vitamin B supplementations [11]. One study in a mouse model showed that vitamin A promoted the proliferation of osteoclasts and thereby reduced BMD [12]. A cross sectional study including 6481 Korean adults found that vitamin A was associated with low BMD and increased fracture risk [13]. However, a 16-year randomized control trial observed no increase in fracture risk after long-term supplementation with high dose of vitamin A [14].

These results indicate that vitamins play critical roles in bone health maintenance. However, most studies have focused on a single or two vitamins [13, 15]. A drawback in these studies is the lack of understanding of the associations of multivitamin coexposure with bone health due to simultaneous exposure to multivitamins from diverse sources for individuals in reality. As far as we know, no study has explored the association between multiple vitamins and BMD, and our study is an attempt to fill the knowledge gap. Current studies have adopted traditional regression models to investigate vitamin and bone associations, and thus collinearity often occur, particularly regression coefficient sign reversal and variance inflation [10, 11]. Therefore, appropriate approaches, such as weighted quantile sum (WQS) regression modelling and principal component analysis (PCA), are needed in exploring the associations of multivitamin coexposure with BMD measurements and the role of each vitamin in vitamin–BMD associations. WQS regression can be used in assessing the dependency of an outcome on a weighted index indicating coexposure to multiple relevant vitamins. Using the idea of dimensionality reduction, the PCA can be used in developing a set of comprehensive indicators from multiple correlated vitamins and maximizing the coverage of original vitamin information. These methods have been gradually used in estimate the association of various multiple exposures, such as metals, and dietary carotenoid intake with bone health [16, 17].

In this study, we utilized the data of the US nationwide population from the National Health and Nutrition Examination Survey (NHANES) to study the individual and mixed associations of six vitamins (B12, B9, C, D, A and E) on BMD in US adults aged ≥ 20 years by using multiple linear regression, restricted cubic splines, PCA models and WQS regression.

Materials and methods

Study population

Data from the NHANES (https://www.cdc.gov/nchs/nhanes/index.htm) was used, which is a cross-sectional survey administered by the Health Statistics Center of American Centers for Disease Control and Prevention (CDC) to noninstitutionalized US residents. A stratified and multistage sampling design was adopted for the selection of representative samples. Baseline data on health and nutritional status were obtained by in-person interviews, mobile physical examination and laboratory testing. The NHANES was approved by the Ethics Review Committee of American National Health Statistics Center. Each participant provided informed consent.

To date, the NHANES has been conducted for 11 rounds every 2 years since 1999. Given that the serum levels of the six vitamins, namely, B12, B9, C, D, A, and E, and the BMD measurement of the femur site were fully measured only in the NHANES 2005–2006, data from this cycle were used in this study. A total of 10,348 participants were included in the 2005–2006 cycle, and 5369 participants aged 20 years or below were excluded. In addition, we excluded 1734 participants who lacked information on BMD measurements and serum levels of the studied vitamins. Furthermore, we excluded 398 participants with missing covariates and 90 participants with osteoporosis treatments. Finally, a total of 2757 participants aged ≥ 20 years were included for the statistical analysis. The flow chart of participant exclusion is shown in Fig. 1.

Fig. 1
figure 1

Flow chart for the selection of US population in the study

Vitamin measurements

Serum samples were collected in parallel with in-person interviews. And serum specimens were then processed, stored at –30 ℃ and shipped to the Division of Laboratory Sciences and CDC for analysis. The serum levels of the six vitamins were detected. Vitamin B12 and red blood cell folate (i.e., vitamin B9) were measured by using a ‘Quantaphase II Folate/Vitamin B12’ radioassay kit (Bio-Rad Laboratories) (https://wwwn.cdc.gov/Nchs/Nhanes/2005-2006/FOLATE_D.htm; https://wwwn.cdc.gov/Nchs/Nhanes/2005–2006/B12_D.htm). Vitamins A (retinol) and E (α-tocopherol) were measured by high-performance liquid chromatography (HPLC) with photodiode array detection (https://wwwn.cdc.gov/Nchs/Nhanes/2005-2006/VITAEC_D.htm). Vitamin C (ascorbic acid) was measured by isocratic HPLC with electrochemical detection at 650 mV1 (https://wwwn.cdc.gov/Nchs/Nhanes/2005-2006/VIC_D.htm). Serum 25-OH-D was measured by using an equilibrium radioimmunoassay procedure (https://wwwn.cdc.gov/Nchs/Nhanes/2005-2006/VID_D.htm). The process of vitamin measurements was described in detail in the NHANES laboratory manuals (https://wwwn.cdc.gov/nchs/nhanes/continuousnhanes/manuals.aspx?BeginYear=2005).

The limit of detection (LOD) for vitamins B12, B9, A, D and E were 20 pg/mL, 0.05 ng/mL, 1.03 μg/dL, 3.75 nmol/L, and 40.67 μg/dL, respectively. The lower limit of detection (LLOD) for vitamin C was 0.03 mg/dL in two rounds of NHANES. For the lower LOD or LLOD values, an imputed value of LOD/LLOD divided by the square root of 2 was adopted.

BMD measurements

All participants received BMD measurements at the skeleton site of the femur via dual-energy X-ray absorptiometry (DXA) in the mobile examination center. Bone scans and scan analyses were performed using a Hologic QDR-4500A fan-beam densitometer (Hologic, Inc., Bedford, Massachusetts), and BMD measurement results were analyzed with Hologic Discovery v12.4 software (https://wwwn.cdc.gov/Nchs/Data/Nhanes/Dxa/DXXFEM_D.htm). Given that the total femur and femoral neck BMD measurements are relatively indicative of hip fracture and osteoporosis diagnosis in clinical setting, we used the total femur and femoral neck BMDs as outcome variables [18]. In the NHANES, BMD measurements were performed on the left femur and femoral neck region, unless answers to safety exclusion questions warrant the use of the right femoral neck.

Covariates

Covariates in adjustment included age, sex, race, family poverty-to-income ratio (PIR), education, body mass index (BMI), calcium concentration, physical activity, drinking status, smoking status, regular milk consumption, osteoporosis family history, daily prednisone/cortisone intake, prevalent diseases of diabetes, and hypertension. They were selected as covariates according to the previous literature [19,20,21,22,23] and biological considerations. Sociodemographic factors (age, sex, race, family PIR, education, and osteoporosis family history), lifestyle factors (physical activity, drinking status, smoking status, regular milk consumption, and prednisone/cortisone intake daily), and disease status (prevalent diseases of diabetes, and hypertension) were obtained by in-person interviews. BMI data were obtained by physical examination, and calcium concentration was obtained by laboratory analysis.

Age was categorized into three groups (20–39, 40–59 and ≥ 60 years), and modelled as a continuous variable. Sex was categorized as male or female. Race was classified into five groups (Mexican American, Hispanic, non-Hispanic white, non-Hispanic black, and other races, including non-Hispanic multiracial). Education was classified into three groups (lower than high school, high school, and above high school). The family PIR was classified by poverty threshold (values 0 ~ 0.99 indicating below poverty and those ≥ 1.00 indicating as or above poverty) [19]. The BMI was calculated as weight (kg) divided by height squared (m2) and classified as no obese(< 30 kg/m2) and obese (≥ 30 kg/m2). Total calcium was measured by ion activity in solution with Beckman Synchron LX20 (https://wwwn.cdc.gov/Nchs/Nhanes/2005-2006/BIOPRO_D.htm), and was treated as a continuous variable. Physical activity was classified into four groups (sedentary = ‘no regular physical activity’; insufficient = ‘ < 500 metabolic equivalent (MET)-minutes per week’; moderate = ‘500–1000 metabolic equivalent (MET)-minutes per week’; and high = ‘ > 1000 metabolic equivalent (MET)-minutes per week’) according to the 2008 Physical Activity Guidelines for Americans [20]. Drinking status was categorized into three groups (current = ‘ > 12 alcohol drink in a year”; ever = ‘ < 12 alcohol drinks in a year but > 12 alcohol drinks in life’; and never = ‘ < 12 alcohol drinks in a year and in life’) [21]. Smoking status was categorized into three groups (current = ‘ > 100 cigarettes in life and smoked every day or some days during the survey time’; ever = ‘ > 100 cigarettes in life and did not smoke during the survey time at all’; and never = ‘ < 100 cigarettes in life’) [22] Regular milk consumption was assessed by the question ‘Have you use regular milk 5 times per week?’ (current, never, and ever). Family history of osteoporosis was assessed by the question ‘Have your parents ever told by a health professional that they had osteoporosis?’ (yes, no). Prednisone/cortisone intake daily was assessed by the question ‘Have you ever taken prednisone or cortisone daily?’ (yes, no) [23]. Prevalent diabetes (yes and no), and hypertension (yes and no) were assessed by the questions ‘Have doctors told you have diabetes?’ and ‘Have doctors ever told you had high blood pressure?’, respectively [23].

Statistical analysis

The guidelines for the NHANES data analysis were retrieved from to (https://wwwn.cdc.gov/nchs/nhanes/tutorials/default.aspx). The summary statistics and distributions of BMD levels were computed for all the participants. Continuous variables were described by using mean and standard deviation (SD), and categorical variables were described by using frequencies and percentages. The concentrations of the six vitamins were described by using median and interquartile range (IQR) according to the variables of the participants’ characteristics. Differences in participants’ characteristics by sex were analyzed by using chi-square test (categorical variables) or t-test (continuous variables). Mutual correlations among the six vitamins were assessed by using Pearson correlation coefficients.

We used multiple linear regression to assess the dependency of total femur and femoral neck BMDs on each vitamin after adjusting all 15 covariates. Four groups by quartile (Q1, Q2, Q3, and Q4) were established for each vitamin, and the first group for each vitamin was used as the reference. Sampling weight was used in analyzing the data from complex survey. We also used restricted cubic splines to analyze the exposure relationship between log-transformed vitamins and BMDs.

WQS regression

WQS regression is an analytic method integrated with linear or logistic regression. It was used in estimating associations between multiple vitamin coexposure and BMD measurements [24]. The concentrations of the six vitamins were natural log-transformed because of their right-skewed distributions. WQS regression modelling consisted of two steps: i) the significance tests for the regression coefficient between the WQS index and outcome, and ii) the calculation of corresponding weights of important components in bootstrap samples. The weighted index was assessed with the following equation.

$$g\left(\mu \right)={\beta }_{0}+{\beta }_{1}\left(\sum_{i=1}^{c}{w}_{i}{q}_{i}\right)+{z}^{\mathrm{^{\prime}}}\varphi$$

Where \({\beta }_{0}\) and \({\beta }_{1}\) are the intercept and regression coefficient, respectively, \({z}^{^{\prime}}\) and \(\varphi\) are the vector of covariates and its coefficient, respectively, and \(c\) is the number of vitamins considered in the analysis (i.e., six in this study). The term \(\sum_{i=1}^{c}{w}_{i}{q}_{i}\) represented the sum of the entire weighted index, and the weights were constrained to sum to 1. The variable \(g\left(\mu \right)\) is a linear link function when outcomes were continuous (BMD). We split the data into two datasets (40% as training set and 60% as validation set), estimated the average of weights of each vitamin in the training dataset, and obtained \({\beta }_{1}\) in the validation dataset through 1000 times of bootstrapping.

We constrained the model was constrained to analyze the effect on outcomes in a single direction at one time. Given that we had no adequate knowledge of the directions of the multiple vitamins-BMD associations, we conducted the model with \({\beta }_{1}\) in positive and negative directions, and determined the two potentials of association directions. One of the directions was the assumption that the components of WQS index were all positively associated with BMD, whereas the other was the assumption that the components of WQS index were all negatively associated with BMD.

PCA

On the basis of a dimensionality reduction principle, PCA recombines several related original variables into a group of unrelated composite variables and extracts optimal composite variables to reflect original variable information as much as possible. PCA was used in deriving vitamin patterns on the basis of six log-transformed serum vitamins. A varimax rotation was applied to increase variance, obtain simple structure in the factor loading matrix and to save factor scores for each participant. To identify the number of components to be retained, we adopted the most widely used criterion (eigenvalues > 1.0). Two principal components (PCs) were retained for multiple linear regression analysis.

SAS 9.4 and R (4.0.4) were used for all statistical analyses. WQS regression was conducted with R packages ‘gWQS’ (version 3.0.3). All significance levels were set to less than 0.05.

Results

General characteristics of participants

The general characteristics of participants are presented in Table 1. A total of 2757 individuals (including 1230 women) were included in this study and had a mean age of 47.97 years. Men had a higher proportion of high physical activity and higher calcium level and were more likely to be alcohol drinkers, smokers and regular milk consumers compared with women (P value < 0.001). The mean BMD at the total femur and femoral neck was higher in men than in women (P value < 0.001). And Table S1 shows the overall distribution of six vitamins, and Table S2 shows the baseline levels of serum vitamins according to the baseline variables of the subjects.

Table 1 Characteristics of participants in the NHANES 2005–2006, stratified by sex

The mutual correlations among each two of the six vitamins were all statistically significant, except the correlation between vitamins A and B12, with a correlation coefficient (r) ranging from 0.03 to 0.41 (P value < 0.001, Fig. S1).

Association between each vitamin and BMD measurements

We assessed the association between the circulating level of each vitamin and BMDs through multivariable linear regression analysis (Table 2). High vitamin B12 quartiles were positively associated with total femur (β Q4 VS Q1 = 0.021, 95% CI: 0.006, 0.036; P value = 0.005) and femoral neck (β Q4 VS Q1 = 0.021, 95% CI: 0.007, 0.035; P value = 0.003) BMDs after adjusting all the mentioned covariates. High vitamin C quartiles were positively associated with total femur (β Q4 VS Q1 = 0.024, 95% CI: 0.009, 0.040; P value = 0.002) and femoral neck (β Q4 VS Q1 = 0.020, 95% CI: 0.005, 0.034; P value = 0.009) BMDs. The associations of vitamins B9, A, D and E levels with BMD measurements did not achieve significance.

Table 2 Associations between the six vitamins and BMDs among participants in NHANES 2005–2006

Figures 2 and 3 present the exposure–response relationship between serum vitamin and BMDs. An inverted L-shaped exposure–response relationship was observed for vitamin C and total femur, and femoral neck BMDs (P for nonlinearity < 0.05).

Fig. 2
figure 2

Adjusted exposure–response relationship between log-transformed vitamins and total femur BMD

Fig. 3
figure 3

Adjusted exposure–response relationship between log-transformed vitamins and femoral neck BMD

Multiple vitamin coexposure and BMD association assessed by PCA

Table S3 shows the factor loading matrix for the two retained PCs. We treated the two PCs as continuous variables and fitted the multiple linear regression to assess their associations with BMDs. The results are shown in Table 3. The factor 1, which had heavy loadings on ‘fat-soluble vitamins’ as vitamins E, A, and D, was associated with decreased total femur (β = 0.002, 95% CI: − 0.004, 0.008; P value = 0.480) and femoral neck (β =  − 0.004, 95% CI: − 0.010, 0.002; P value = 0.161) BMDs, although the effect values did not reach statistical significance. The factor 2, which had heavy loadings on the ‘water-soluble vitamins’ as vitamins B12, B9, and C, was associated with increased total femur (β = 0.009, 95% CI: 0.004, 0.015; P value = 0.001) and femoral neck (β = 0.007, 95% CI: 0.002, 0.013; P value = 0.005) BMDs.

Table 3 Association between the two major principal component and BMDs among participants in NHANES 2005–2006

Multiple vitamin coexposure and BMD association assessed by WQS regression

Table 4 summarizes the results of the WQS regression analysis of covariate-adjusted associations between multivitamin coexposure and BMD measurements. Specifically, the positive associations of statistical significance were found for the WQS indexes composed by the circulating levels of the six vitamins and total femur (βWQS = 0.010, 95% CI: 0.001, 0.018; P value = 0.021) and femoral neck (βWQS = 0.008, 95% CI: 0.001, 0.015; P value = 0.022) BMDs.

Table 4 Associations between the weighted quantile sum regression indexes composed for multiple vitamins and BMDs by WQS among participants in NHANES 2005–2006

The estimated weights for vitamins in composing BMD-associated WQS indexes are illustrated in Table S4 and Fig. 4. Vitamins B12 and C were predominant members for the positive association with total femur BMD (weights of 0.51 and 0.35, respectively) and femoral neck (weights of 0.65 and 0.32, respectively). Vitamins B12, C, and B9 were predominant members (weights of 0.47, 0.35, and 0.14, respectively) for the positive association with LS BMD, although of no significance was found. Vitamins E and A weighted the most in the multivitamin coexposure and total femur (weights of 0.65 and 0.20, respectively) and femoral neck (weights of 0.56 and 0.29, respectively) BMDs association, although of no significance was observed.

Fig. 4
figure 4

The WQS regression estimated weights of each of the six vitamins associated with total femur BMD (A) and femoral neck BMD (B) in the positive direction and total femur BMD (C) and femoral neck BMD (D) in the negative direction

Discussion

This study mainly aimed to determine the associations of coexposure with multiple vitamins (i.e., B12, B9, C, D, A, and E) and DXA-derived BMD measurements in US adults included in NHANES 2005–2006. Multiple linear regression analysis showed the positive associations of circulating vitamins B12 and C with BMDs, and an inverted L-shaped exposure–response relationship was observed between vitamin C and BMDs. WQS regression and PCA revealed a significant positive association between coexposure to water-soluble vitamins and BMDs, and the association was mainly driven by vitamins B12 and C. A negative association was found between coexposure to fat-soluble vitamins and BMDs, and was mainly driven by vitamins E and A.

Water-soluble vitamins and BMD

Water-soluble vitamins consist of vitamin B and C family members. The positive associations between vitamin B members (B12 and B9) and BMD were consistent with the results of previous studies [10]. WQS analysis showed that vitamin B12 had the highest weight in the positive associations between multivitamins and BMDs, indicating that vitamin B12 may have a beneficial effect on bone health. A randomized trial reported that supplementation with a high-dose mixture of B vitamins contributed to elevated BMDs [25]. Vitamin B was a cofactor in homocysteine metabolism, and vitamin B12 and B9 deficiency caused elevated plasma homocysteine, which may interfere with collagen cross-linking, and damage the bone matrix and result in decreased BMD [26]. Previous trial in adults showed that treatment with B vitamins can decrease plasma homocysteine levels [27]. A previous study found that no association between vitamin B after adjusting for vitamins C and A, inconsistent with our findings [28]. Assumptions showed that they may indirectly regulate bone metabolism by acting on other vitamins. Several studies reported that vitamin B supplementation had no effect that lowered homocysteine level or changed bone markers [11, 29]. However, these studies have small sample size and restricted subjects.

Also known as ascorbic acid, vitamin C is a water-soluble vitamin. In the present study, we found that vitamin C contributes greatly and positively to the association of multivitamin coexposure with BMDs. To date, many studies have demonstrated that vitamin C contributes to elevated BMDs [9, 15]. Vitamin C can stimulate types 1 and 2 collagen synthesis, limit the lifetimes of osteoclasts, and prevent the loss of osteoblast markers [30]. Moreover, vitamins C and D can synergistically enhance osteocalcin (OC) promoter activity and carboxylated OC synthesis and thereby improve bone mass [28]. However, a recent study found a strong correlation between vitamin C and OC after adjusting serum vitamin D, indicating that vitamin C plays a more direct and important role in OC [31]. Vitamin C, which has strong biological activity, promotes the absorption of calcium in the body, and alleviates the toxicity of external toxic substances, is a nutrient that benefits bone metabolism [2]. Clinically, vitamin C deficiency causes scurvy and a range of skeletal symptoms, such as osteoporosis [32]. We found an inverted L-shaped exposure–response relationship between vitamin C and BMD. This association may be related to oxalate accumulation due to a high level of vitamin C. A case report showed that a high level of vitamin C induced oxalate accumulation, and resulted in bone pain [33]. Differences among these studies may be related to the dual roles of vitamin C in osteoclasts, and thus more studies must be performed.

Fat-soluble vitamins and BMD

All the three statistical methods suggested significant negative association between vitamin E and BMD, consistent with the previous study [34]. WQS analysis shows that vitamin E had the highest weight in the negative associations between multivitamins and BMDs, indicating that vitamin E may have a damaging effect on bone health. In our study, a-tocopherol levels in sera were detected as the biomarkers of internal exposure to vitamin E. In mammals, a-tocopherol is the dominant isoform of vitamin E and is the major form of vitamin E supplements in the U.S. [35]. Serum a-tocopherol improved osteoclast activity, induced osteoclast fusion and differentiation and exerted an unfavorable effect on the bone [6]. However, a study with inconsistent conclusions showed that a-tocopherol was positively associated with BMD [36]. Therefore, studies are warranted to determine the association of vitamin E and its isomers with bone health.

As the main fat-soluble vitamins, vitamins A and D play interactive functions with each other. The binding of vitamin D receptor (VDR) and nuclear retinoid X receptor (RXR) forms the VDR/RXR complex, which further combines with vitamin D to exert vitamin D’s effect [28]. The RXR is also the receptor of vitamin A, and elevated levels of vitamin A may compete with vitamin D to bind to the VDR/RXR complex, thereby inhibiting vitamin D’s effect [13, 37]. On the one hand, WQS analysis and PCA revealed suggestive negative association of vitamin A with BMD, which is consistent with some epidemiologic studies [37, 38]. The possible mechanisms may be that its metabolic intermediate retinoic acid can increase osteoclast proliferation and enhance bone resorption, which weakens the bone [12]. On the other hand, our study found no significant positive association of serum vitamin D with BMD. To the best of our knowledge, vitamin D plays an important role in BMD by regulating calcium and phosphorus status and promoting bone mineralization [39]. Vitamin D deficiency can lead to osteomalacia [40], which would lead to decreased BMD. Meanwhile, a 1-year randomized control trial have reported that a dietary pattern with vitamin D reduced bone loss in the elderly with osteoporosis [41]. The difference in our results may be due to the inhibitory effect of vitamin A on vitamin D. Given that they can bind to osteocalcin-related receptor RXR [28], a competitive receptor relationship antagonizes some actions with each other. One study provided evidence that excessive VA may inhibit the favorable effect of vitamin D on BMD and potentiate the vitamin D deficiency–insufficiency effect [42]. Vitamins A and D are mutually restricted. Patients with severe vitamin D deficiency should avoid excessive intake of vitamin A when taking vitamin D as supplementary therapy. Interestingly, a 4-year cohort study found an inverted U-shaped association between vitamin A and BMD, indicating that low level of vitamin A reduces BMD [43]. Therefore, vitamin A may play dual roles on bone metabolism and exerts a positive effect on BMD at low levels and an inhibition effect on vitamin D at high levels. Thus, more studies are needed to determine its optimal level.

Strengths and limitations

Our study has several strengths. The WQS regression and PCA methods were used in determining the associations of multivitamin coexposure and BMD measurements, accounting for the collinearity of multivitamins. The WQS approach was used in evaluating potentially important contributors to the association between multivitamin co-exposure and BMD. Serum vitamin concentrations were used to represent vitamin exposure levels, which was more accurate and indicative of internal exposure level than food-frequency questionnaire survey.

Our study has some potential limitations. We used the data from a cross sectional survey, and thus result on the causal inference of multivitamin exposure on BMD was limited. Large prospective cohort studies and experimental studies are needed to accumulate evidence. We included six types of vitamins in this study and omitted other vitamins (e.g., vitamin K) due to the absence of eligible data. US participants were studied, and the generalization of our findings to other ethnic populations is limited.

Conclusion

The associations between exposure to multiple vitamins and BMDs were explored by using generalized linear regression, restricted cubic splines, WQS regression, and PCA methods with the NHANES 2005–2006 data of US adults aged ≥ 20 years. Vitamin B12 and C levels were associated positively, and an inverted L-shaped exposure–response relationship was observed between vitamin C and BMD. The water-soluble vitamins were significantly associated with increased BMDs, and association was mainly driven by B12 and C. The fat-soluble vitamins were associated with decreased BMDs, and the association was mainly driven by vitamin E and A. Our study provided evidence of the associations between multiple vitamin coexposure and BMDs.

Availability of data and materials

Publicly available datasets were analyzed in this study. This data can be found here: https://www.cdc.gov/nchs/nhanes/index.htm

References

  1. Johnell O, Kanis JA. An estimate of the worldwide prevalence and disability associated with osteoporotic fractures. Osteoporos Int. 2006;17(12):1726–33. https://doi.org/10.1007/s00198-006-0172-4.

    Article  CAS  PubMed  Google Scholar 

  2. Malmir H, Shab-Bidar S, Djafarian K. Vitamin C intake in relation to bone mineral density and risk of hip fracture and osteoporosis: a systematic review and meta-analysis of observational studies. Br J Nutr. 2018;119(8):847–58. https://doi.org/10.1017/S0007114518000430.

    Article  CAS  PubMed  Google Scholar 

  3. Leibson CL, Tosteson ANA, Gabriel SE, Ransom JE, Melton LJ. Mortality, disability, and nursing home use for persons with and without hip fracture: a population-based study. J Am Geriatr Soc. 2002;50(10):1644–50. https://doi.org/10.1046/j.1532-5415.2002.50455.x.

    Article  PubMed  Google Scholar 

  4. Levin VA, Jiang X, Kagan R. Estrogen therapy for osteoporosis in the modern era. Osteoporos Int. 2018;29(5):1049–55. https://doi.org/10.1007/s00198-018-4414-z.

    Article  CAS  PubMed  Google Scholar 

  5. Chen LR, Hou PH, Chen KH. Nutritional support and physical modalities for people with osteoporosis: current opinion. Nutrients. 2019;11(12):15. https://doi.org/10.3390/nu11122848.

    Article  Google Scholar 

  6. Fujita K, Iwasaki M, Ochi H, Fukuda T, Ma CS, Miyamoto T, Takitani K, et al. Vitamin E decreases bone mass by stimulating osteoclast fusion. Nat Med. 2012;18(4):589–94. https://doi.org/10.1038/nm.2659.

    Article  CAS  PubMed  Google Scholar 

  7. Iwamoto J, Yeh JK, Takeda T. Comparative effects of vitamin K and vitamin D supplementation on prevention of osteopenia in calcium-deficient young rats. J Musculoskelet Neuronal Interact. 2003;3(4):423–423. https://doi.org/10.1016/s8756-3282(03)00249-7.

    Article  Google Scholar 

  8. Urban K, Hohling HJ, Luttenberg B, Szuwart T, Plate U, Biomineralisation Res Unit. An in vitro study of osteoblast vitality influenced by the vitamins C and E. Head Face Med. 2012; 8: 10. https://doi.org/10.1186/1746-160X-8-25

  9. Brzezinska O, Lukasik Z, Makowska J, Walczak K. Role of Vitamin C in osteoporosis development and treatment-a literature review. Nutrients. 2020;12(8):22. https://doi.org/10.3390/nu12082394.

    Article  CAS  Google Scholar 

  10. Bailey RL, Looker AC, Lu ZH, Fan RZ, Eicher-Miller HA, Fakhouri TH, Gahche JJ, et al. B-vitamin status and bone mineral density and risk of lumbar osteoporosis in older females in the United States. Am J Clin Nutr. 2015;102(3):687–94. https://doi.org/10.3945/ajcn.115.108787.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  11. van Wijngaarden JP, Swart KM, Enneman AW, Dhonukshe-Rutten RA, van Dijk SC, Ham AC, Brouwer-Brolsma EM, et al. Effect of daily vitamin B-12 and folic acid supplementation on fracture incidence in elderly individuals with an elevated plasma homocysteine concentration: B-PROOF, a randomized controlled trial. Am J Clin Nutr. 2014;100(1):1578–86. https://doi.org/10.3945/ajcn.114.090043.

    Article  CAS  PubMed  Google Scholar 

  12. Hu LJ, Lind T, Sundqvist A, Jacobson A, Melhus K. Retinoic acid increases proliferation of human osteoclast progenitors and inhibits RANKL-stimulated osteoclast differentiation by suppressing RANK. PLoS ONE. 2010;5(10):10. https://doi.org/10.1371/journal.pone.0013305.

    Article  CAS  Google Scholar 

  13. Joo NS, Yang SW, Song BC, Yeum KJ. Vitamin A intake, serum Vitamin D and bone mineral density: analysis of the korea national health and nutrition examination survey (KNHANES, 2008–2011). Nutrients. 2015;7(3):1716–27. https://doi.org/10.3390/nu7031716.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  14. Ambrosini GL, Bremner AP, Reid A, Mackerras D, Alfonso H, Olsen NJ, Musk AW, et al. No dose-dependent increase in fracture risk after long-term exposure to high doses of retinol or beta-carotene. Osteoporos Int. 2013;24(4):1285–93. https://doi.org/10.1007/s00198-012-2131-6.

    Article  CAS  PubMed  Google Scholar 

  15. Deyhim F, Strong K, Deyhim N, Vandyousefi S, Stamatikos A, Faraji B. Vitamin C reverses bone loss in an osteopenic rat model of osteoporosis. Int J Vitam Nutr Res. 2018;88(1–2):58–64. https://doi.org/10.1024/0300-9831/a000486.

    Article  CAS  PubMed  Google Scholar 

  16. Wei Mh, Cui Y, Zhou HL, Song WJ, Wang Q. Associations of multiple metals with bone mineral density: a population-based study in US adults. Chemosphere. 2021;282:131150. https://doi.org/10.1016/j.chemosphere.2021.131150.

  17. Sugiura M, Nakamura M, Ogawa K, Ikoma Y, Ando F, Shimokata H, Yano M. Dietary patterns of antioxidant vitamin and carotenoid intake associated with bone mineral density: findings from post-menopausal Japanese female subjects. Osteoporos Int. 2011;22(1):143–52. https://doi.org/10.1007/s00198-010-1239-9.

    Article  CAS  PubMed  Google Scholar 

  18. Marshall D, Johnell O, Wedel H. Meta-analysis of how well measures of bone mineral density predict occurrence of osteoporotic fractures. BMJ. 1996;312(7041):1254–9. https://doi.org/10.1136/bmj.312.7041.1254.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  19. Cui Y, Zhou HL, Wei MH, Song WJ, Di DS, Zhang RY, et al. Multiple vitamin co-exposure and mortality risk: a prospective study. Clin Nutr. 2022;41(2):337–47. https://doi.org/10.1016/j.clnu.2021.12.010.

    Article  CAS  PubMed  Google Scholar 

  20. Tucker LA. Physical activity and telomere length in U.S. men and women: An NHANES investigation. Prev Med. 2017;100:145–51. https://doi.org/10.1016/j.ypmed.2017.04.027.

    Article  PubMed  Google Scholar 

  21. Zhou HL, Wei MH, Cui Y, Di DS, Song WJ, Zhang RY, et al. Association between water intake and mortality risk-evidence from a national prospective study. Front Nutr. 2022;9:822119. https://doi.org/10.3389/fnut.2022.822119.

  22. Zhang Y, Dong T, Hu W, Wang X, Xu B, Lin Z, et al. Association between exposure to a mixture of phenols, pesticides, and phthalates and obesity: comparison of three statistical models. Environ Int. 2019;123:325–36. https://doi.org/10.1016/j.envint.2018.11.076.

    Article  CAS  PubMed  Google Scholar 

  23. Fan Y, Ni S, Zhang H. Associations of copper intake with bone mineral density and osteoporosis in adults: data from the national health and nutrition examination survey. Biol Trace Elem Res. 2022;200(5):2062–8. https://pubmed.ncbi.nlm.nih.gov/34283365/.

  24. Carrico C, Gennings C, Wheeler DC, Factor-Litvak P. Characterization of weighted quantile sum regression for highly correlated data in a risk analysis setting. J Agric Biol Environ Stat. 2015;20(1):100–20. https://doi.org/10.1007/s13253-014-0180-3.

    Article  PubMed  Google Scholar 

  25. Herrmann M, Umanskaya N, Traber L, Schmidt- Gayk H, Menke W, Lanzer G, Lenhart M, et al. The effect of B-vitamins on biochemical bone turnover markers and bone mineral density in osteoporotic patients: a 1-year double blind placebo controlled trial. Clin Chem Lab Med. 2007;45(12):1785–92. https://doi.org/10.1515/CCLM.2007.352.

    Article  CAS  PubMed  Google Scholar 

  26. Thaler R, Zwerina J, Rumpler M, Spitzer S, Gamsjaeger S, Paschalis EP, Paschalis EP, et al. Homocysteine induces serum amyloid A3 in osteoblasts via unlocking RGD-motifs in collagen. FASEB J. 2013;27(2):446–63. https://doi.org/10.1096/fj.12-208058.

    Article  CAS  PubMed  Google Scholar 

  27. Albert CM, Cook NR, Gaziano JM, Elaine Z, MacFadyen J, Danielson E, Buring JE, et al. Effect of folic acid and B vitamins on risk of cardiovascular events and total mortality among women at high risk for cardiovascular disease: a randomized trial. JAMA. 2008;299(17):2027–36. https://doi.org/10.1001/jama.299.17.2027.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  28. Johnson NA, Chen BH, Sung SY, Liao CH, Hsiao WC, Chung WK, Hsieh CL. A novel targeting modality for renal cell carcinoma: human osteocalcin promoter-mediated gene therapy synergistically induced by vitamin C and vitamin D3. J Gene Med. 2010;12:892–903. https://doi.org/10.1002/jgm.1516.

    Article  CAS  PubMed  Google Scholar 

  29. Salari P, Abdollahi M, Heshmat R, Meybodi HA, Razi F. Effect of folic acid on bone metabolism: a randomized double blind clinical trial in postmenopausal osteoporotic women. Daru J Pharm Sci. 2014;22:7. https://doi.org/10.1186/s40199-014-0062-9.

    Article  CAS  Google Scholar 

  30. Hie M, Tsukamoto I. Vitamin C-deficiency stimulates osteoclastogenesis with an increase in RANK expression. J Nutr Biochem. 2011;22(2):164–71. https://doi.org/10.1016/j.jnutbio.2010.01.002.

    Article  CAS  PubMed  Google Scholar 

  31. Torbergsen AC, Watne LO, Wyller TB, Frihagen F, Strømsøe K, Bøhmer T, Mowe M. Micronutrients and the risk of hip fracture: case-control study. Clin Nutr. 2017;36(2):438–43. https://doi.org/10.1016/j.clnu.2015.12.014.

    Article  CAS  PubMed  Google Scholar 

  32. Snoddy AME, Buckley HR, Elliott GE, Standen VG, Arriaza BT, Halcrow SE. Macroscopic features of scurvy in human skeletal remains: a literature synthesis and diagnostic guide. Am J Phys Anthropol. 2018;167(4):876–95. https://doi.org/10.1002/ajpa.23699.

    Article  PubMed  PubMed Central  Google Scholar 

  33. Kennedy SS, Perilloux A, Pereira RC, Handelman G, Wesseling-Perry K, Salusky IB. Vitamin C overload may contribute to systemic oxalosis in children receiving dialysis. Pediatr Nephrol. 2021;36(2):435–41. https://doi.org/10.1007/s00467-020-04702-4.

    Article  PubMed  Google Scholar 

  34. Zhang J, Hu X, Zhang J. Associations between serum vitamin E concentration and bone mineral density in the US elderly population. Osteoporos Int. 2017;28(4):1245–53. https://doi.org/10.1007/s00198-016-3855-5.

    Article  CAS  PubMed  Google Scholar 

  35. Qing J, Stephan C, Shigenaga MK, Ames BN. γ- tocopherol, the major form of vitamin E in the US diet, deserves more attention. Am J Clin Nutr. 2001;74(6):714–22. https://doi.org/10.1093/ajcn/74.6.714.

    Article  Google Scholar 

  36. Yang TC, Duthie GG, Aucott LS, Macdonald HM. Vitamin E homologues alpha- and gamma-tocopherol are not associated with bone turnover markers or bone mineral density in peri- menopausal and post-menopausal women. Osteoporos Int. 2016;27(7):2281–90. https://doi.org/10.1007/s00198-015-3470-x.

    Article  CAS  PubMed  Google Scholar 

  37. Annu K, Cline C, Yasuda K, Ganguly S, Pesch A, Cooper B, Janke L, et al. Role of Vitamins A and D in BCR-ABL Arf-/- acute lymphoblastic leukemia. Sci Rep. 2020;10(1):2359. https://doi.org/10.1038/s41598-020-59101-4.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  38. Feskanich D, Singh V, Willett WC, Colditz GA. Vitamin A intake and hip fractures among postmenopausal women. JAMA. 2002;287:47–54. https://doi.org/10.1001/jama.287.1.47.

    Article  CAS  PubMed  Google Scholar 

  39. Michaelsson K, Lithell H, Vessby B, Melhus H. Serum retinol levels and the risk of fracture. N Engl J Med. 2003;348:287–94. https://doi.org/10.1056/NEJMoa021171.

    Article  CAS  PubMed  Google Scholar 

  40. Metzger M, Houillier P, Gauci C, Haymann JP, Flamant M, Thervet E, Boffa JJ, et al. Relation between circulating levels of 25(OH) Vitamin D and parathyroid hormone in chronic kidney disease: quest for a threshold. J Clin Endocrinol Metab. 2013;98(7):2922–8. https://doi.org/10.1210/jc.2013-1294.

    Article  CAS  PubMed  Google Scholar 

  41. Jennings A, Cashman KD, Gillings R, Cassidy A, Tang J, Fraser W, et al. A Mediterranean-like dietary pattern with vitamin D3 (10 µg/d) supplements reduced the rate of bone loss in older Europeans with osteoporosis at baseline: results of a 1-y randomized controlled trial. Am J Clin Nutr. 2018;108(3):633–40. https://pubmed.ncbi.nlm.nih.gov/30007343/.

  42. Holick MF, Binkley NC, Bischoff-Ferrari HA, Gordon CM, Hanley DA, Heaney RP, Murad MH, et al. Guidelines for preventing and treating Vitamin D deficiency and insufficiency revisited. J Clin Endocrinol Metab. 2012;97(4):1153–8. https://doi.org/10.1210/jc.2011-2601.

    Article  CAS  PubMed  Google Scholar 

  43. Metz AL, Walser MM, Olson WG. The interaction of dietary vitamin A and vitamin D related to skeletal development in the turkey poult. J Nutr. 1985;115(7):929–35. https://doi.org/10.1093/jn/115.7.929.

    Article  CAS  PubMed  Google Scholar 

Download references

Acknowledgements

Not applicable

Funding

This work was supported by the National Natural Science Foundation of China (Grant no. 82273711) and the Funding for Scientific Research Projects from Wuhan Municipal Health Commission (Grant no. WY22B06).

Author information

Authors and Affiliations

Authors

Contributions

Ruyi Zhang: Investigation, Methodology, Writing- Original draft preparation, Writing – Review & Editing. Qin Huang: Investigation, Methodology, Writing- Original draft preparation, Writing – Review & Editing. Guanhua Su: Investigation, Methodology, Writing- Original draft preparation, Writing – Review & Editing. Muhong Wei: Investigation, Methodology, Writing – Review & Editing. Yuan Cui: Investigation, Writing – Review & Editing. Haolong Zhou: Investigation, Writing – Review & Editing. Wenjing Song: Investigation, Writing – Review & Editing. Dongsheng Di: Investigation, Writing – Review & Editing. Junan Liu: Conceptualization, Methodology, Writing—Review & Editing. Qi Wang: Conceptualization, Methodology, Writing—Review & Editing. The author(s) read and approved the final manuscript.

Corresponding author

Correspondence to Qi Wang.

Ethics declarations

Ethics approval and consent to participate

The study was conducted according to the guidelines of the Declaration of Helsinki, and approved by the Ethics Review Committee of American National Health Statistics Center. Informed consent was obtained from all subjects involved in this study.

Consent for publication

Not applicable.

Competing interests

The authors declare no competing interests.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1:

Fig.S1. Pearson correlations among blood concentrations of the six studied vitamins (B12, B9,C, D, A and E) in this research. All P values <0.001 except for thecorrelation between vitamins A and B12. Table S1. Distribution of concentrationof the six studied vitamins among participants in NHANES 2005-2006. Table S2. Distribution of concentrationof the six studied vitamins accordingto the characteristics of participants in the NHANES 2005–2006.Table S3. Orthogonal rotated factor-loading matrix for serumvitamins among participants in NHANESs 2003-2006.Table S4. Weight of vitamins in WQS regression for the associationswith BMDs among participants in NHANES 2005-2006.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Zhang, R., Huang, Q., Su, G. et al. Association between multiple vitamins and bone mineral density: a cross-sectional and population-based study in the NHANES from 2005 to 2006. BMC Musculoskelet Disord 24, 113 (2023). https://doi.org/10.1186/s12891-023-06202-6

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12891-023-06202-6

Keywords