Abstract

The relationship between standardized serum 25-hydroxyvitamin D (25(OH)D) concentration and incident anemia in the United States (U.S.) is unclear. The purpose of our study was to examine the association between serum 25(OH)D and anemia risk. We performed a cross-sectional analysis of the U.S. population participating in the National Health and Nutrition Examination Survey (NHANES) between 2001 and 2018. A generalized linear model and restricted cubic spline (RCS) plot curve were constructed to assess the relationship between serum 25(OH)D concentration and anemia incidence. Additionally, the association between serum 25(OH)D concentration and red blood cell (RBC) count and hemoglobin (HB) levels was investigated using generalized additive models with smooth functions. Subgroup analysis also was performed. A total of 29933 individuals were included in our research. After adjusting for known confounding variables, compared with the lowest quartile, the odds ratios (ORs) with 95% confidence intervals (CIs) for association of serum 25(OH)D with anemia across the second, third, and fourth quartiles were 0.735 (0.651, 0.829), 0.527 (0.461, 0.602), and 0.696 (0.611, 0.792), respectively. Serum 25(OH)D concentration was associated with anemia risk in a U-shaped pattern, as shown by an RCS plot ( for nonlinearity <0.001). In addition, RBC count and Hb levels initially increased and then decreased as serum 25(OH)D levels increased. Serum 25(OH)D concentration and risk of anemia are associated with a U-shaped curve in the U.S. general population. Serum 25(OH)D concentration in the range 59.7–70.3 nmol/l was associated with anemia incidence <1. Therefore, the risk of anemia can be reduced by close monitoring and appropriate vitamin D supplementation.

1. Introduction

Worldwide, anemia affects one out of every three people and is estimated to represent a greater global disease burden than major depression and chronic respiratory illnesses combined [1, 2]. Symptoms of anemia include fatigue, weakness, cognitive decline, and even death, and the condition is more common in women [3]. Anemia pathogenesis involves an imbalance between red blood cell loss and production, which may be due to ineffective or inadequate erythropoiesis and/or excessive red blood cell loss [1]. The most common conditions causing anemia are iron deficiency, thalassemia and hemoglobinopathies, folate deficiency, and parasitic diseases [4].

As a fat-soluble vitamin, vitamin D is naturally found in very few foods; the canonical role of vitamin D is to regulate the balance of calcium and phosphorus metabolism and minerals in bones [5]. Vitamin D deficiency has been linked to anemia [69], since vitamin D is found in the bone marrow and can promote red blood cell production [1012]. The main ring form of vitamin D, 25(OH)D is created by hydroxylation of vitamin D, which can be taken in from exogenous sources or synthesized in the liver [13]. In recent years, the relationship between 25(OH)D and anemia has attracted increasing attention. Many previous studies have shown a correlation between 25(OH)D and anemia in chronic kidney disease as well as sickle cell anemia [1418]; however, few recent studies have explored the link between serum 25(OH)D levels and anemia risk, and previous studies had small sample sizes or concomitant disease interference.

Due to the detrimental effects of anemia, recognizing risk factors for anemia development and devising measures to promptly avoid or control its negative consequences would be highly advantageous. According to epidemiological research, the relationship between serum 25(OH)D and anemia risk in the United States (U.S.) population remains unclear. Therefore, in this study, we analyzed data from the Nutrition and Health Examination Survey (NHANES) 2001–2018 to investigate the link between serum 25(OH)D and anemia incidence in the U.S. general population.

2. Materials and Methods

2.1. Study Population

NHANES is an American cross-sectional survey that collects data on the health and nutrition of the general population through stratified multistage random sampling (https://www.cdc.gov/nchs/nhanes/). Eight two-year cycles (2001–2018) of data from NHANES were used for this analysis. Participants with insufficient serum 25(OH)D or anemia data were excluded (n = 10632 and 146, respectively). Additionally, participants <18 years old and with missing demographic or biochemical data were excluded (n = 41606). Finally, this research included a total of 29933 individuals. Data on parathyroid hormone (PTH) levels, history of osteoporosis, and arthritis or rheumatism obtained by the affiliated hospital of Xuzhou Medical University between 2012 and 2022 were also included and analyzed in this study. NHANES was authorized by the National Center for Health Statistics study ethical review board, and each participant provided signed written informed permission [19]. All tests were taken at a mobile testing facility on-site.

2.2. Serum 25(OH)D Concentrations

Blood samples were collected during the examination, centrifuged, divided, and frozen at –70°C on-site, then shipped on dry ice to a central laboratory, where they were stored at –70°C until analysis. After acetonitrile-based extraction, a radioimmunoassay kit (DiaSorin, Stillwater, MN) was used to measure serum 25(OH)D levels at the National Center for Environmental Health (Atlanta, GA); however, due to concern about bias and imprecision of the DiaSorin radioimmunoassay (RIA), the Centers for Disease Control (CDC) developed regression equations to convert RIA values to liquid chromatography with tandem mass spectrometry (LC-MS/MS) equivalents for NHANES 2001–2006. In addition, the CDC laboratory analyzed serum 25(OH)D metabolites from 2007 to 2018, using LC-MS/MS, and calculated total serum 25(OH)D (nmol/L) as the sum of 25(OH)D3 and 25(OH)D2, excluding the C3-epi-25(OH)D3 metabolite. Detailed information and procedures are described on the NHANES website (https://wwwn.cdc.gov/nchs/nhanes/vitamind/analyticalnote.aspx).

2.3. Covariates

The following covariates were included in the study: age, sex, race/ethnicity, family poverty income ratio (PIR), education level, marital status, hypertension, diabetes mellitus (DM), smoking status, and alcohol consumption status; history of coronary heart disease (CHD), congestive heart failure (CHF), angina pectoris, heart attack, stroke, chronic kidney diseases (CKD), and osteoporosis; body mass index (BMI); waist circumference; folic acid, vitamin B12, vitamin C, iron; and serum iron, calcium, phosphorus, and PTH; and number of days with arthritis or rheumatism. In addition, red blood cell (RBC) numbers and hemoglobin (Hb) levels of participants were collected. During the home interview, the following data were self-reported by participants: age, sex, ethnicity, education level, marital status, smoking status, alcohol consumption status, and dietary intake; history of CHD, CHF, angina pectoris, heart attack, stroke, CKD, and osteoporosis; and number of days with arthritis or rheumatism. Data on serum iron, calcium, phosphorus, and PTH were obtained from laboratory tests. More information about the variables analyzed in this research can be found at https://www.cdc.gov/nchs/nhanes/.

2.4. Anemia Ascertainment

The World Health Organization defines anemia as hemoglobin level <13 g/dL in men and <12 g/dL in women [20]. A Beckman Coulter DxH 800 instrument in the NHANES mobile examination center produces a complete blood count from blood specimens, which provided blood cell distributions for all participants. Refer to the Laboratory Method Files section for a detailed description of the laboratory methods used.

2.5. Statistical Analysis

All analyses were performed using R version 3.6.4 (R Foundation for Statistical Computing, Vienna, Austria) and Stata version 13.0 (Stata Corporation, College Station, TX, USA). Samples with missing covariate data were excluded from this study. -value <0.05 was regarded as statistically significant. Serum 25(OH)D levels were divided into quartiles, and the lowest quartile served as the reference group (Q1). All estimates were calculated accounting for NHANES sample weights. Continuous variables are expressed as mean (standard deviation) and categorical variables are presented as numbers (%). Weighted linear regression models (continuous variables) and weighted chi-square tests (categorical variables) were used to assess the significance of differences. Multivariate logistic regression analysis was used to investigate the relationship between serum 25(OH)D levels and anemia risk. Model 1 was adjusted for age and sex. Model 2 was adjusted for model 1 variables plus race/ethnicity, education level, marital status, family PIR, smoking status, alcohol consumption status, history of hypertension, and DM. Finally, model 3 was adjusted for model 2 variables plus history of CHD, CHF, angina pectoris, heart attack, stroke, and CKD; BMI; waist circumference; folic acid, vitamin B12, vitamin C, iron; and serum iron, calcium, and phosphorus, as our core model.

3. Results

3.1. Baseline Characteristics

Table 1 shows the baseline characteristics of the research participants. The incidence of anemia in this group was 8.3%. The characteristics of the participants were subclassified based on serum 25(OH)D quartiles (Q1: 6.31–44.4 nmol/L; Q2: 44.5–60.6 nmol/L; Q3: 60.7–77.7 nmol/L; and Q4: 77.8–422 nmol/L). There were significant differences in age, sex, race, family PIR, education level, marital status, smoking status, alcohol consumption status, hypertension, DM, CHF, stroke, CKD, BMI, waist circumference, RBC, Hb, serum iron, calcium, phosphorus, mean energy intake, protein intake, folic acid intake, vitamin B12 intake, and iron intake among the Q1, Q2, Q3, and Q4 groups. Compared with Q1, Q3, and Q4 groups, the Q2 group had a lower proportion of participants with hypertension, heart attack, stroke, and CKD, as well as the highest RBC counts; however, compared with the Q1, Q2, and Q4 groups, a lower proportion of individuals in the Q3 group had DM, CHF, angina, and anemia, while they had the highest levels of Hb, mean energy intake, protein intake, folic acid intake, vitamin B12 intake, and iron intake. In addition, compared with Q1, Q2, and Q3, participants in Q4 were older and had the highest levels of family PIR, calcium, phosphorus, serum iron, and vitamin C intake, and had lower BMI and waist circumference. The basic characteristics of the 959 participants from the affiliated hospital of Xuzhou Medical University are shown in Supplementary Table 1, including the levels of PTH, history of osteoporosis, and number of days with arthritis or rheumatism.

3.2. Association between Serum 25(OH)D and Anemia

The findings of multivariate logistic regression analysis of the relationship between serum 25(OH)D and anemia are presented in Table 2. After adjusting for interfering factors, the odds ratios (ORs) with 95 percent confidence intervals (CIs) for association of serum 25(OH)D with anemia across quartiles two, three, and four were 0.735 (0.651, 0.829), 0.527 (0.461, 0.602), and 0.696 (0.611, 0.792), compared with Q1. As shown by the restricted cubic spline (RCS) plot, serum 25(OH)D was associated with anemia prevalence with a U-shaped curve ( for nonlinearity <0.001, Figure 1). As serum 25(OH)D concentrations increased, the risk of anemia initially decreased significantly; risk of anemia was lowest when serum 25(OH)D concentrations reached 65.0 nmol/L, then the curve showed an upward trend. In the fully adjusted model including PTH, history of osteoporosis, and number of days with arthritis or rheumatism, the ORs and 95 percent CIs for anemia incidence across the serum 25(OH)D quartiles were 0.771 (0.449, 1.326), 0.499 (0.267, 0.932), and 0.596 (0.323, 1.097), compared with the first quartile of serum 25(OH)D (Supplementary Table 2). Additionally, RCS plots constructed after adjusting for PTH, history of osteoporosis, and number of days with arthritis or rheumatism as covariates also displayed U-shaped curves (Supplementary Figure 1).

3.3. Associations of Serum 25(OH)D with RBC Count and Hb Level

Generalized additive models with smooth functions suggested nonlinear relationships between serum 25(OH)D and risk of anemia. As serum 25(OH)D levels increased, both RBC count and Hb level first increased and then decreased (Figures 2(a) and 2(b)). Additionally, generalized additive models with smooth functions generated following adjustment for PTH, history of osteoporosis, and number of days with arthritis or rheumatism as covariates also showed an increase, after declining (Supplementary Figures 2(a) and 2(b)).

3.4. Subgroup Analyses

Subgroup analyses, stratified by age, sex, race/ethnicity, hypertension, DM, and obesity were conducted to estimate the associations between serum 25(OH)D and incident anemia. The results demonstrated U-shaped relationships between serum 25(OH)D and risk of anemia in individuals of all ages, male or female, who were Mexican American, Other Hispanic, Non-Hispanic Black, or Non-Hispanic White, without hypertension, without DM, and never participated, with or without hypertension, with or without diabetes, and whether or not obese (Table 3). There were significant interactions of age, sex, race/ethnicity, hypertension, DM, and obesity with the associations with serum 25(OH)D in subgroup analyses. In addition, we performed subgroup analysis stratified by age, sex, hypertension, DM, and obesity, to examine the correlation between serum 25(OH)D and the risk of anemia, after adjusting for variables, including PTH, history of osteoporosis, and number of days with arthritis or rheumatism (Supplementary Table 3).

4. Discussion

In recent years, studies of the relationship between 25(OH)D and anemia have focused on sickle cell anemia, and have demonstrated that 25(OH)D supplementation is helpful in treatment of sickle cell anemia [16, 18, 2123]; however, the overall relationship between 25(OH)D and anemia, particularly its trend, has not been clarified. Our study revealed a U-shaped relationship between serum 25(OH)D levels and risk of anemia, and with increased 25(OH)D level, hemoglobin level, and RBC counts the association showed a downward trend, suggesting that caution should be applied in the use of 25(OH)D supplementation to treat anemia. Insufficient 25(OH)D levels can reduce local calcitriol production in the bone marrow and limit erythropoiesis. As calcitriol has a synergistic effect with endogenous erythropoietin, it can up-regulate the expression of erythropoietin receptors on erythrocyte progenitor cells [2426]; however, there is no clear explanation for the relationship between higher levels of 25(OH)D and increased risk of anemia, which appears contradictory to the established effect of 25(OH)D. We consider that this anomalous outcome may be attributable to concomitant diseases in participants, many of which led to end-stage metabolic disorders at physiological hormonal and cellular levels. It is likely that 25(OH)D inhibits other hematopoietic processes through a feedback system, thereby causing a decrease in red blood cell counts and hemoglobin levels, as well as an increase in the risk of anemia. According to our study, when serum 25(OH)D level is controlled between 59.7 and 70.3 nmol/l, the risk of anemia can be effectively reduced; this finding requires further confirmation in basic experiments.

Anemia is considered an independent cause of morbidity and mortality in the elderly, and its incidence varies with sex, race, and age [27, 28]. Anemia can reduce the ability of adults to work and affect the growth and development of children and adolescents, and is more common in women [29, 30]. Similarly, 25(OH)D may have different basal levels and sensitivities in different populations. For example, Atkinson et al. showed that black children have a lower 25(OH)D threshold than white children [8]. Our subgroup analysis also suggested significant interactions of anemia with age, sex, race/ethnicity, hypertension, diabetes, obesity, and serum 25(OH)D. For the Hispanics (other) ethnic category, the 25(OH)D threshold was significantly lower than that for other races.

Immune diseases are a cause of anemia, and are often accompanied by abnormal metabolism of vitamin D and PTH, leading to osteoporosis [31, 32]. Further, low levels of 25(OH)D may increase the risk of immune-related diseases because they exert an immunomodulatory effect through nuclear vitamin D receptor expressed by antigen-presenting cells and activated T and B cells [33, 34]. Therefore, we screened the database for data from patients with immune diseases and hyperparathyroidism; after adjusting for these covariates, the results were generally consistent with those above, confirming that arthritis and rheumatism are not additional risk factors for reduced plasma 25(OH)D concentrations [35]. Notably, 25(OH)D levels in patients with arthritis and rheumatism were higher than those in the whole group when anemia risk was in the trough, which may be due to the effects of the pathological processes underlying arthritis or rheumatism on anemia, or the regulation of 25(OH)D by PTH; however, neither significantly affected the overall trend relationship between anemia and 25(OH)D. Therefore, the use of vitamin D supplementation to treat anemia needs to be personalized for different populations, particularly in patients with comorbidities that may affect vitamin D metabolism. In summary, our study shows that there is no simple linear relationship between 25(OH)D and risk of anemia, and that age, sex, race, hypertension, diabetes, arthritis, rheumatism, and obesity all interact with 25(OH)D levels. Further basic and prospective experiments are needed to further explore the relationship between 25(OH)D and the mechanism underlying anemia development.

Our study has some limitations. First, we used the NHANES public database for this analysis, which covers the years 2001 to 2018. To verify our findings, participants from other nations should be recruited. Second, a retrospective study has the disadvantage of introducing bias to some relevant results. Third, a study of the mechanisms underlying serum 25(OH)D levels and anemia incidence is also necessary.

5. Conclusion

In conclusion, the relationship between serum 25(OH)D and risk of anemia presents as a U-shaped curve. An inflection point for serum 25(OH)D was observed in our study; the incidence of anemia was lowest when the serum 25(OH)D level was 65.0 nmol/l, and anemia risk was <1 when serum 25(OH)D levels were between 59.7 and 70.3 nmol/l. Therefore, with close monitoring and adequate vitamin D supplementation, the risk of anemia can be reduced.

Data Availability

The survey data are publicly available on the Internet for data users and researchers throughout the world (https://www.cdc.gov/nchs/nhanes/).

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Authors’ Contributions

Yang Zhang and Xiaotong Wang contributed to hypothesis development and manuscript preparation. Jia Wang contributed to study design. Weiwei Hu, Xiaolu Song, Ding Yuan, and Xianliang Yan undertook data analyses. Yang Zhang drafted and revised the manuscript. All authors approved the final draft of the manuscript for publication. Yang Zhang, Xiaotong Wang, and Jia Wang contributed equally to this work.

Acknowledgments

This work was not specifically funded. The authors thank the staff and the participants of the NHANES study for their valuable contributions, Xiue Du and Qiang Guo for their contributions of charts modifications and the full-text format adjustments.

Supplementary Materials

Supplementary Figure 1: Restricted cubic spline plot of the association between serum 25(OH)D and the incidence of anemia. Supplementary Figure 2: Associations of serum 25(OH)D with hemoglobin levels and red blood cell counts. (a) Association between serum 25(OH)D and hemoglobin level. (b) Association between serum 25(OH)D and red blood cell count. Supplementary Table 1: Study population data according to serum 25(OH)D quartiles. Supplementary Table 2: Adjusted ORs for associations between serum 25(OH)D and the risk of anemia. Supplementary Table 3: Subgroup analysis for associations between serum 25(OH)D and the risk of anemia. (Supplementary Materials)