==== Front
BMC Public Health
BMC Public Health
BMC Public Health
1471-2458
BioMed Central London

37380964
16165
10.1186/s12889-023-16165-x
Research
L-shaped association of serum 25-hydroxyvitamin D with all-cause and cardiovascular mortality in older people with chronic kidney disease: results from the NHANES database prospective cohort study
http://orcid.org/0000-0002-4921-4355
Li Rugang rugangli1@163.com

1
Li Yang leeyoungmm@163.com

2
Fan Zhongcheng fanzc090606@163.com

3
Liu Zhaoqi qwe1148437587@163.com

4
Lin Juhua juhualin2023@163.com

4
He Min minhe20030701@126.com

14
1 grid.411679.c 0000 0004 0605 3373 Department of Nephrology, Affiliated Yuebei People’s Hospital of Shantou University Medical College, No. 133 South Huimin Road, Shaoguan, 512026 Guangdong China
2 Department of Nephrology, Haikou Municipal People’s Hospital and Central South University Xiangya Medical College Affiliated Hospital, Haikou, China
3 Department of Osteology, Haikou Municipal People’s Hospital and Central South University Xiangya Medical College Affiliated Hospital, Haikou, China
4 grid.410560.6 0000 0004 1760 3078 Guangdong Medical University, Zhanjiang, 524001 Guangdong China
28 6 2023
28 6 2023
2023
23 126018 4 2023
21 6 2023
© The Author(s) 2023
https://creativecommons.org/licenses/by/4.0/ 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.
Background

This study was conducted to assess the association of serum 25-hydroxyvitamin D [25(OH)D] concentrations with all-cause and cardiovascular disease (CVD) mortality in older people with chronic kidney disease (CKD) in the United States.

Methods

We identified 3230 CKD participants aged ≥ 60 years from the National Health and Nutrition Examination Survey (2001–2018). CKD was defined as an estimated glomerular filtration rate (eGFR) < 60 ml/min/1.73 m2. Mortality outcomes were determined by linkage to National Death Index (NDI) records through December 31, 2019. Restricted cubic spline based on Cox regression models were utilized to elucidate the nonlinear relationship between serum 25(OH)D concentrations and mortality in patients with CKD.

Results

During median 74 months of follow-up, 1615 all-cause death and 580 CVD death were recorded. We found an L-shaped association between serum 25(OH)D concentrations and all-cause and CVD mortality, reaching a plateau at 90 nmol/L. Accordingly, per one-unit increment in natural log-transformed 25(OH)D was associated with a 32% and 33% reduced risk of all-cause mortality (hazard ratio [HR] 0.68; 95%CI, 0.56 to 0.83) and CV mortality (HR 0.69; 95%CI, 0.49 to 0.97) in participants with serum 25(OH)D < 90 nmol/L, but no considerable difference was observed in participants with serum 25(OH)D ≥ 90 nmol/L. Compared with those in the deficiency group (< 50 nmol/L), insufficient (50 to < 75 nmol/L) and sufficient group (≥ 75 nmol/L) were significantly associated with lower all-cause mortality (HR,0.83; 95%CI, 0.71 to 0.97 and HR, 0.75; 95%CI, 0.64 to 0.89) and CV mortality (HR,0.87; 95%CI, 0.68 to 1.10 and HR, 0.77; 95%CI, 0.59 to < 1.0), respectively.

Conclusion

An L-shaped relationship between serum 25(OH)D levels with all-cause and CVD mortality was observed in elderly CKD patients in the United States. A 25(OH)D concentration of 90 nmol/L may be the target to reduce the risk of premature death.

Keywords

Chronic kidney disease
25-hydroxyvitamin D
Mortality
Cardiovascular diseases
Older people
Hainan provincial natural science foundation[819QN384, 820QN422] Li Yang http://dx.doi.org/10.13039/501100013138 Hainan Provincial Department of Health [19A200125, 20A200321] Li Yang issue-copyright-statement© BioMed Central Ltd., part of Springer Nature 2023
==== Body
pmcBackground

Chronic kidney disease (CKD) affects more than 10% of the general population worldwide, and this situation is more common among older people [1, 2]. CKD has become one of the major causes of death from non-communicable diseases worldwide [3]. It is critical to identify modifiable factors to prevent premature death from CKD in older adults.

25-hydroxyvitamin D [25(OH)D] is the main storage form of vitamin D and is involved in phosphate and calcium metabolism. Vitamin D deficiency is associated with musculoskeletal health, cardiovascular diseases, hypertension, autoimmune diseases, obesity, diabetes and depression [4–7]. A single-center retrospective study and post-hoc subgroup analysis of the LURIC study both found that low 25(OH)D levels were associated with increased all-cause mortality [8, 9]. Additionally, a meta-analysis of 10 studies found that higher 25(OH)D levels were associated with better survival in patients with CKD [10]. The study of NHANES III (Third National Health and Nutrition Examination Survey 1988–1994) showed non-dialysis CKD patients with low serum 25(OH)D levels had an increased risk of death compared to the high level group [11]. Lifestyles, the spectrum of disease, and treatments have dramatic changes during the last two decades. It is necessary to re-assess the relationship between 25(OH)D and risk of death. Additionally, the effectiveness of vitamin D supplementation on CVD and all-cause mortality is controversial. This could be attributed to uncertainty about the optimal serum 25(OH)D concentration necessary for intervention. Given that most published studies are small, relatively old data, and inadequate adjustments for important factors in CKD, further studies are needed to investigate the optimal vitamin D concentration for patients with CKD is crucial.

We conducted a large prospective cohort study, based on the National Health and Nutrition Examination Survey (NHANES) to assess the association between serum 25(OH)D levels with all-cause and CVD mortality, and to evaluate the optimal target level of vitamin D supplementation in elderly patients with CKD.

Methods

Study design and population

NHANES is a complex, multistage sampling design nationally representative study to assess the health and nutritional status of the noninstitutionalized civilian population in the U.S. From 1999, it became a continuous program, with every two years representing a cycle [12]. NHANES was performed by the Centers for Disease Control and Prevention (CDC), and approved by the institutional review board of the National Center of Health Statistics. All participants provided written informed consent.

We used nine cycles of the NHANES from 2001 to 2018 (vitamin D data were not available in the NHANES 1999–2000). There are 17,222 individuals aged ≥ 60 years who are eligible for follow up. 3646 participants with CKD diagnosis [estimated glomerular filtration rate (eGFR) < 60 ml/min/1.73 m2] and not dialysis in the previous 12 months were included in this analysis. Those with missing information on age, sex, race, education, BMI, biochemical test, comorbidities, co-medications (n = 354), or abnormal values (n = 62) were further excluded. Therefore, a total of 3230 participants entered in our present analysis.

Measurement of Serum 25(OH)D

In NHANES, the concentration of serum 25(OH)D was determined by DiaSorin radioimmunoassay kit (Stillwater, MN) in 2001–2006, and subsequently by standardized liquid chromatography-tandem mass spectrometry (LC–MS/MS) since the 2007–2008 cycle. In accordance with CDC recommendations, we used LC–MS/MS equivalent data for all analyses. The details of the conversion are specified elsewhere [13].

Ascertainment of mortality

Mortality from any cause and CVD was ascertained by linkage to the National Death Index through 31 December 2019. The ICD-10 was used to determine disease-specific death. CVD mortality was defined as the primary cause of death being any disease of the circulatory system (ICD-10 codes I00-I09, I11, I13, I20-I51, or I60–I69).

Assessment of covariates

Information on age, sex, race, education level, smoking status, recreational activity, comorbidities, and medication use was collected from household interviews using standardized questionnaires. Body weight, height were obtained when people participated in physical examinations at a mobile examination center. BMI was calculated as weight in kilograms divided by height in meters squared. Race was classified as non-Hispanic White, non-Hispanic Black, Mexican American, or other. Education level was categorized as less than high school, high school or equivalent, or college or above. Smoking status was classified as never smoker, former smoker, or current smoker. Physical activity was defined based on the three common domains of participation corresponding to the interviewer-administered Physical Activity Questionnaire (PAQ), including leisure and recreational activities, occupation-related physical activity, and transportation-related physical activity. This study focused only on leisure and recreational activities and defined as “yes” as long as participation in the activities [14].

In addition, creatinine, triglycerides, total cholesterol, UACR, and parathyroid hormone (PTH) were measured at baseline when the participants provided their blood samples. eGFR was calculated by the CKD Epidemiology Collaboration equation [15].

Diabetes was defined as self-reported doctor diagnosis of diabetes, use of insulin or oral hypoglycemic medication, fasting glucose ≥ 7.0 mmol/L, random glucose ≥ 11.1 mmol/L, or glycated hemoglobin A1c (HbA1c) ≥ 6.5%. Hypertension status was obtained from self-report, or systolic blood pressure ≥ 140 mmHg or diastolic blood pressure ≥ 90 mm Hg. CVD contains self report coronary heart disease, congestive heart failure, heart attack, stroke, and angina. CCI was calculated based on the literature published in 1994 [16].

Statistical analysis

Sample weights, clustering, and stratification were incorporated into all analyses because of the complex sampling design of the NHANES. Participants were followed up to death or the date of 31 December 2019, whichever comes first. According to the Endocrine Society Clinical Practice Guidelines, vitamin D status was categorized into three groups: deficiency (< 50 nmol/L), insufficient (50 to < 75 nmol/L), and sufficient (≥ 75 nmol/L) [17]. Serum 25(OH)D concentrations were also analyzed as a continuous variable after natural log transformation. We estimated the relationship between serum 25(OH)D concentrations and all-cause mortality as well as CVD mortality using Cox proportion hazard model, with or without adjustment for age, sex, age, race, education, smoking status, BMI, and recreational activity, UACR, eGFR, total cholesterol, triglycerides, comorbidities (hypertension, diabetes, CVD, CCI), and co-medications. Additionally, a restricted cubic spline based on Cox regression, with 3 knots, was performed to test for linearity and explore the shape of the dose–response relation of serum 25(OH)D concentrations and mortality.

Stratified analyses were also conducted by sex, ethnicity (White or non-White), BMI (< 30.0, or ≥ 30.0 kg/m2), eGFR(< 45, ≥ 45 ml/min/1.73m2), and with or without CVD, hypertension and diabetes.

Two sensitivity analyses were performed to test the robustness of our findings. First, to reduce the potential reverse causal bias, we included participants with more than 2 years of follow-up. Second, considering the interrelationship of PTH and vitamin D status, PTH was further adjusted (only available in the NHANES 2003–2006, n = 737).

All analyses were carried out with R version 4.2.0, and a two-tailed P < 0.05 was considered statistically significant.

Results

Among the 3230 elderly CKD participants (mean age, 74.5 years; 52.1% female), 732 (23.9%) had deficient vitamin D (< 50 nmol/L) and 1061 (32.8%) had insufficient vitamin D (< 75 nmol/L). The baseline characteristics of the study population according to serum 25(OH)D status are shown in Table 1. Participants who had higher 25(OH)D levels were more likely to be female, non-Hispanic White; had higher education levels and recreational activity; were less likely to be current smoker and obese; and lower comorbidities.Table 1 Baseline characteristics of participants with CKD stratified by serum 25(OH)D concentrations

Variablesa	Serum 25(OH)D concentrations (nmol/L)b	P value	
 < 50 (n = 772)	50 to < 75 (n = 1061)	 ≥ 75 (n = 1397)	
Age, yr	73.9 (0.3)	74.5 (0.3)	74.6 (0.2)	0.21	
Female, n (%)	451 (67.6)	494 (53.6)	739 (59.8)	 < 0.001	
Race, n (%)				 < 0.001	
 White	357 (70.1)	721 (84.3)	1019 (87.6)		
 Black	233 (17.5)	146 ( 6.8)	185 ( 5.6)		
 Mexican American	99 (4.1)	84 (2.7)	65 (1.7)		
 Other	83 (8.3)	110 (6.2)	128 (5.2)		
Education, n (%)				 < 0.001	
 Low than high school	332 (34.5)	371 (27.2)	360 (19.4)		
 High school	182 (27.0)	274 (28.2)	374 (27.7)		
 College	258 (38.6)	416 (44.7)	663 (52.9)		
Smoke, n (%)				 < 0.001	
 Never	373 (48.1)	506 (49.5)	692 (52.2)		
 Former	297 (37.7)	466 (43.0)	617 (42.7)		
 Current	102 (14.2)	89 ( 7.5)	88 ( 5.2)		
BMI, kg/m2				 < 0.001	
  < 25	176 (22.3)	250 (22.6)	376 (26.2)		
 25–29	240 (29.8)	400 (37.0)	558 (38.5)		
  ≥ 30	356 (48.0)	411 (40.4)	463 (35.3)		
Recreational activity, n (%)	109 (13.4)	248 (25.6)	455 (36.3)	 < 0.001	
eGFR, ml/min/1.73m2	46.7 (0.5)	48.49 (0.4)	48.0 (0.3)	0.01	
UACR, mg/g	263.5 (46.1)	115.9 (18.0)	75.0 ( 9.6)	 < 0.001	
TCHO, mmol/L	4.9 (0.1)	5.0 (0.0)	4.9 (0.1)	0.12	
TG, mmol/L	1.8 (0.1)	1.9 (0.1)	1.8 (0.0)	0.31	
CCI, count	2.5 (0.1)	2.1 (0.1)	2.4 (0.1)	0.003	
Hypertension, n (%)	655 (83.3)	843 (78.2)	1137 (80.3)	0.1	
Diabetes, n (%)	348 (44.2)	388 (32.8)	491 (31.3)	 < 0.001	
CVD, n (%)	321 (42.9)	421 (37.5)	516 (34.7)	0.01	
Anti-hypertension, n (%)	579 (76.3)	751 (68.2)	1014 (68.8)	0.01	
Lipid-lowering drugs, n (%)	347 (48.2)	503 (45.9)	764 (55.0)	 < 0.001	
Hypoglycemic, n (%)	239 (29.8)	271 (23.6)	335 (21.6)	0.002	
Abbreviation: BMI Body mass index, UACR Urinary albumin creatinine ratio, eGFR Estimated glomerular filtration rate, TCHO Total cholesterol, TG Triglycerides, CVD Cardiovascular disease, CCI Charlson Comorbidity Index

aAll estimates accounted for complex survey designs. Continuous variables were expressed as mean (standard error). Categorical variables were expressed as number (percent)

bAll intervals are left closed and right open

During a median follow-up of 74 months (interquartile range, 38–119 months), 1615 all-cause deaths were documented, including 580 CV deaths. The association between serum 25(OH)D concentrations (Fig. 1) with the risk of all-cause and CV mortality followed a L-shape (P for nonlinearity 0.002 and 0.01), with the smallest Akaike information criterion (AIC) of three knots. Accordingly, in the threshold effect analysis, per one-unit increment in natural log-transformed 25(OH)D was associated with a 32% and 33% reduced risk of all-cause mortality (hazard ratio [HR] 0.68; 95%CI, 0.56 to 0.83) and CV mortality (HR 0.69; 95%CI, 0.49 to 0.97) in participants with serum 25(OH)D < 90 nmol/L, but no considerable difference was observed in participants with serum 25(OH)D ≥ 90 nmol/L (Table 2).Fig. 1 Relation of serum 25(OH)D concentrations with all-cause and CVD mortality. Hazard ratio (HR) was adjusted for age, sex, race, education, smoking status, BMI, recreational activity, UACR, eGFR, total cholesterol, triglycerides, comorbidities (hypertension, diabetes, CVD, CCI), and co-medications. Abbreviation: CVD, Cardiovascular disease

Table 2 Threshold effect analyses of serum 25(OH)D with all-cause and CVD mortality using two-piecewise regression models

Serum 25(OH)D, nmol/L	Crude HR (95%CI)a	Adjusted HR (95%CI)ab	
 < 90 (n = 2381)	
 All-cause mortality	0.62 (0.54,0.72)	0.68 (0.56,0.83)	
 CVD mortality	0.67 (0.50,0.88)	0.69 (0.49,0.97)	
 ≥ 90 (n = 849)	
 All-cause mortality	0.83 (0.31,2.23)	1.43 (0.55,3.72)	
 CVD mortality	1.18 (0.30,4.62)	2.66 (0.68,10.44)	
Abbreviation: CVD Cardiovascular disease

aPer one-unit increment in natural log-transformed 25(OH)D

bAdjusted for age, sex, race, education, smoking status, BMI, recreational activity, UACR, eGFR, total cholesterol, triglycerides, comorbidities (hypertension, diabetes, CVD, CCI), and co-medications

In addition, we also included the serum 25(OH)D concentrations as a categorical variable divided into three levels, deficiency (< 50 nmol/L), insufficient (50 to < 75 nmol/L) and sufficient (≥ 75 nmol/L), and the sufficient category including the inflection point (90 nmol/L). Compared with those in the deficiency group, the insufficient and sufficient groups were significantly associated with lower all-cause mortality (HR,0.83; 95%CI, 0.71 to 0.97 and HR, 0.75; 95%CI, 0.64 to 0.89) and CV mortality (HR,0.87; 95%CI, 0.68 to 1.10 and HR, 0.77; 95%CI, 0.59 to < 1.0), respectively (Table 3).Table 3 The association between serum 25(OH)D concentrations with all-cause and CVD mortality

	Serum 25(OH)D concentrations (nmol/L)	
 < 50 (n = 772)	50 to < 75 (n = 1061)	 ≥ 75 (n = 1397)	
All-cause mortality	
 Number of deaths	485	593	537	
 Crude HR (95%CI)	reference	0.79 (0.69,0.90)	0.68 (0.59,0.79)	
 Adjusted HR (95%CI)a	reference	0.83 (0.71,0.97)	0.75 (0.64,0.89)	
CVD mortality	
 Number of deaths	164	228	188	
 Crude HR (95%CI)	reference	0.86 (0.68,1.09)	0.69 (0.54,0.88)	
 Adjusted HR (95%CI)a	reference	0.87 (0.68,1.10)	0.77 (0.59, < 1.0)	
Abbreviation: CVD Cardiovascular disease

aAdjusted for age, sex, race, education, smoking status, BMI, recreational activity, UACR, eGFR, total cholesterol, triglycerides, comorbidities (hypertension, diabetes, CVD, CCI), and co-medications

In the subgroup analysis, the association between 25(OH)D and all-cause mortality was consistent across subgroups stratified by sex, race, eGFR, presence of CVD, hypertension and diabetes. Although, the interactions for BMI were lower than 0.05, due to decreased tendency in the similar directionality of the associations in those with BMI more than 30 kg/m2, suggesting that there are more factors related to death in obese people besides vitamin D (Fig. 2).Fig. 2 Association between serum 25(OH)D concentrations and all-cause mortality in various subgroups. Hazard ratio (HR) was adjusted for age, sex, race, education, smoking status, BMI, recreational activity, UACR, eGFR, total cholesterol, triglycerides, comorbidities (hypertension, diabetes, CVD, CCI), and co-medications. Abbreviation: eGFR, estimated glomerular filtration rate; CVD, Cardiovascular disease

In the sensitivity analyses, compared with the deficiency level, the results of the insufficient group (HR,0.76; 95%CI, 0.59 to 0.99) and sufficient group (HR,0.65; 95%CI, 0.47 to 0.92) did not significantly change when we included participants with more than 2 years of follow-up. Similarly, when further adjusted for PTH, the results also did not materially change (HR, 0.74; 95%CI, 0.57 to 0.95) and (HR,0.68; 95%CI, 0.49 to 0.94).

Discussion

This large prospective cohort study of US older people with CKD showed a L-shaped association between serum 25(OH)D concentrations with all-cause and CVD mortality, with the plateau at 90 nmol/L. This association was independent of traditional risk factors, and was consistently observed across subgroups stratified by eGFR (< 45, ≥ 45 ml/min/1.73m2), gender, BMI (< 30, ≥ 30 kg/m2), race (White, non-White) and with or without CVD, diabetes, and hypertension.

The relationship between 25(OH)D concentration and mortality remains a topic of interest to clinicians and the public in patients with CKD. Our study showed a high prevalence of insufficient 25(OH)D (32.8%) in patients with CKD, indicating widespread vitamin D deficiency in those people. Consistent with previous research [8, 10, 11, 18], our multivariate analysis showed that higher serum 25(OH)D concentrations were independently associated with reduced all-cause and cardiovascular mortality.

The definition of sufficiency and toxicity thresholds for vitamin D remain unclear [19]. A retrospective cohort study showed an independent association between 25(OH)D levels of less than 37.5 nmol/L and all-cause mortality in non-dialysis CKD patients [8]. In a previous NHANES study [11], subgroups with 25(OH)D levels less than 37.5 nmol/L had a higher risk of all-cause mortality compared to subgroups with 25(OH)D levels higher than 75 nmol/L with or without multiple adjustment. We prepared serum 25(OH)D < 50 nmol/L as the reference group, and a significant reduction of all-cause and CVD mortality in the 50–75 nmol/L group and ≥ 75 nmol/L group was observed. Moreover, we observed an L-shaped relationship between serum 25(OH)D concentration with all-cause and CVD mortality, similar to a large general population survey in the United Kingdom [20]. Their study showed an L-shaped relationship between genetically predicted 25-(OH)D and all-cause mortality, with a sharp increase risk of death at concentrations below 50 nmol/L. In another cohort study (UK Biobank, EPIC-CVD and two Copenhagen population studies) utilizing population-based Mendelian randomization analysis also showed a non-linear dose–response relationship of vitamin D with coronary heart disease, stroke and all-cause mortality [21]. They observed an inverse correlation between genetically predicted 25(OH)D concentrations and all-cause mortality up to 40 nmol/L. This association is usually nonlinear. Mortality decreases with increasing concentrations of 25(OH)D until a certain inflection point, beyond which further mortality declines are not observed. An inflection point at 90 nmol/L was observed in this older CKD population, beyond which no further reduction of death was shown. The finding suggest that vitamin D requirements is differs in various disease, and older patients with CKD may need relatively higher levels of vitamin D. Mendelian randomization studies require consideration of instrumental variables. Because genes are randomized and polymorphic, this may lead to large differences in study results across populations. However, our research is a nationally representative survey with a complex, multi-stage sampling design. The variables are simple, standardized and comprehensive, therefore the findings of the study may be more robust and representative.

There is few evidence shown beneficial effects of vitamin D supplement on clinical outcomes. In the general population, vitamin D supplementation has not significantly reduced mortality from all-causes or CVD [22, 23]. One possibility is that people with vitamin D deficiency or severe deficiency were excluded from the study. Furthermore, the effects of vitamin D supplement on all-cause and CVD mortality in CKD individuals are controversial [19, 24]. This may be because there is no further reduction in mortality when serum vitamin D concentrations reach a certain threshold. Our analysis of data from a nationally representative population-based survey showed an L-shaped relationship between serum 25(OH)D concentrations and all-cause and CVD mortality, which could partially explain these controversies. Future randomized controlled trials need to fully account for this L-shaped relationship.

There are several possible mechanisms explaining the association between lower 25(OH)D levels and an increased risk of mortality. Numerous diseases that cause death have been linked to low vitamin D levels, such as dyslipidemia, hypertension, and diabetes [25–27]. Moreover, low vitamin D levels are associated with coronary artery calcification, thickening of the carotid intima-media and other cardiovascular disease [28–30]. These conditions may therefore constitute a risk of mortality, while high levels of vitamin D may reduce signal between renin, angiotensin and aldosterone, improve endothelial dysfunction and modulate immune function [31–33]. By inhibiting macrophage migration and cholesterol uptake, suppressing foam cell formation, and reversing atherogenic cholesterol metabolism, it attenuates atherosclerosis [34, 35]. There is evidence that secondary hyperparathyroidism associated with low vitamin D levels increases all-cause and CVD mortality in multiple ways. In addition, it activates the RAS system, worsens hypertension, and exacerbates anemia by exacerbating abnormalities in calcium and phosphorus metabolism and thus vascular calcification [36–38]. Nevertheless, further mechanistic studies are needed to better understand how vitamin D prevents all-cause and CVD mortality in patients with CKD.

Strengths of our study include the prospective study design, relatively large sample size, and a nationally representative sample of older Americans with CKD, which facilitates replication. As well, this study had sufficient statistical power due to the substantial number of deaths during a long period of time. Moreover, we adjusted for multiple potential confounders based on the understanding of elderly population, that improved the validity and robustness of our conclusions. There are also some limitations to be considered. First, because it was an observational study, residual confounding factors are possible. Second, serum 25(OH)D concentrations were measured only once, which may underestimate the true association of interest [39]. Third, the current study did not assess vitamin D metabolite ratios (e.g., the ratio of serum 25[OH]D to 1,25[OH]2D could be used as a surrogate for 24-hydroxylase) and genetic variation in genes related to vitamin metabolism [19]. Third, this study defined CKD as eGFR < 60 ml/min/1.73 m2. It should be cautious when generalizing the finding to patients with albuminuria and hematuria alone. Finally, the role of psychosocial stress or genetic susceptibility cannot be excluded in the current study.

Conclusions

In a nationally representative sample of older people with CKD, we found a L-shaped association between serum 25(OH)D concentrations with all-cause and CVD mortality, with the plateau at 90 nmol/L. The threshold of 90 nmol/L may be a target for interventions to reduce the risk of death.

Abbreviations

CKD Chronic kidney disease

CVD Cardiovascular disease

25(OH)D 25-Hydroxyvitamin D

NHANES National Health and Nutrition Examination Survey

NDI National Death Index

NHANES III Third National Health and Nutrition Examination Survey

CDC Centers for Disease Control and Prevention

95% CI 95% Confidence interval

eGFR Estimated glomerular filtration rate

HR Hazard ratio

PIR Family poverty income ratio

BMI Body mass index

UACR Urinary albumin creatinine ratio

eGFR Estimated glomerular filtration rate

TCHO Total cholesterol

TG Triglycerides

CVD Cardiovascular disease

CCI Charlson Comorbidity Index

Acknowledgements

The authors thank the participants of the NHANES databases. Thanks to Zhang Jing (Shanghai Tongren Hospital) for developing the nhanesR package, which makes it easier to explore the NHANES database.

Authors’ contributions

R.G.L., Y.L. and M.H. conceived and designed the study; Y.L., Z.Q.L, and J.H.L. analyzed the data; R.G.L. and Y.L. wrote the paper. R.G.L. and M.H. edited the manuscript; Z.C.F. revised and polished a draft of manuscript. All authors approved the submitted and final versions.

Funding

This work was supported by Hainan provincial natural science foundation [819QN384, 820QN422]; Hainan Provincial Department of Health [19A200125, 20A200321].

Availability of data and materials

The datasets generated and analyzed in the current study are available at NHANES website: https://www.cdc.gov/nchs/nhanes/index.htm.

Declarations

Ethics approval and consent to participate

Study protocols for NHANES were approved by the NCHS Ethnics Review Board (ERB) (Protocol #98–12, #2006–06, #2011–17 and #2018–01, https://www.cdc.gov/nchs/nhanes/irba98.htm). All participants signed the informed consent before participating in the study.

All methods were carried out in accordance with relevant guidelines and regulations.

Consent for publication

Not applicable.

Competing interests

The authors declare no competing interests.

Publisher’s Note

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

Rugang Li, Yang Li, and Zhongcheng Fan contributed equally to this work.
==== Refs
References

1. Jager KJ Kovesdy C Langham R Rosenberg M Jha V Zoccali C A single number for advocacy and communication-worldwide more than 850 million individuals have kidney diseases Kidney Int 2019 96 1048 1050 10.1016/j.kint.2019.07.012 31582227
2. Kovesdy CP Epidemiology of chronic kidney disease: an update 2022 Kidney Int Suppl 2011 2022 12 7 11 10.1016/j.kisu.2021.11.003
3. Johansen KL Chertow GM Gilbertson DT Ishani A Israni A Ku E US renal data system 2022 annual data report: epidemiology of kidney disease in the United States Am J Kidney Dis 2023 81 A8 11 10.1053/j.ajkd.2022.12.001 36822739
4. Holick MF The vitamin D deficiency pandemic: approaches for diagnosis, treatment and prevention Rev Endocr Metab Disord 2017 18 153 165 10.1007/s11154-017-9424-1 28516265
5. 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 g2035 10.1136/bmj.g2035 24690624
6. Al Mheid I Quyyumi AA Vitamin D and cardiovascular disease: controversy unresolved J Am Coll Cardiol 2017 70 89 100 10.1016/j.jacc.2017.05.031 28662812
7. Sîrbe C Rednic S Grama A Pop TL An update on the effects of Vitamin D on the immune system and autoimmune diseases Int J Mol Sci 2022 23 9784 10.3390/ijms23179784 36077185
8. Navaneethan SD Schold JD Arrigain S Jolly SE Jain A Schreiber MJ Low 25-hydroxyvitamin D levels and mortality in non-dialysis-dependent CKD Am J Kidney Dis 2011 58 536 543 10.1053/j.ajkd.2011.04.028 21816525
9. Pilz S Tomaschitz A Friedl C Amrein K Drechsler C Ritz E Vitamin D status and mortality in chronic kidney disease Nephrol Dial Transplant 2011 26 3603 3609 10.1093/ndt/gfr076 21378153
10. Pilz S Iodice S Zittermann A Grant WB Gandini S Vitamin D status and mortality risk in CKD: a meta-analysis of prospective studies Am J Kidney Dis 2011 58 374 382 10.1053/j.ajkd.2011.03.020 21636193
11. Mehrotra R Kermah DA Salusky IB Wolf MS Thadhani RI Chiu Y-W Chronic kidney disease, hypovitaminosis D, and mortality in the United States Kidney Int 2009 76 977 983 10.1038/ki.2009.288 19657329
12. NHANES - National Health and Nutrition Examination Survey Homepage n.d. https://www.cdc.gov/nchs/nhanes/index.htm (Accessed 25 March 2023).
13. NHANES - Analytical Note for 25-Hydroxyvitamin D Data Analysis n.d. https://wwwn.cdc.gov/nchs/nhanes/vitamind/analyticalnote.aspx?b=2013&e=2014&d=VID_H. (Accessed 25 March 2023).
14. PAQ_E n.d. https://wwwn.cdc.gov/Nchs/Nhanes/2007-2008/PAQ_E.htm#PAQ650 (Accessed 30 March 2023).
15. Levey AS Stevens LA Schmid CH Zhang YL Castro AF Feldman HI A new equation to estimate glomerular filtration rate Ann Intern Med 2009 150 604 612 10.7326/0003-4819-150-9-200905050-00006 19414839
16. Charlson M Szatrowski TP Peterson J Gold J Validation of a combined comorbidity index J Clin Epidemiol 1994 47 1245 1251 10.1016/0895-4356(94)90129-5 7722560
17. Holick MF Binkley NC Bischoff-Ferrari HA Gordon CM Hanley DA Heaney RP Evaluation, treatment, and prevention of vitamin D deficiency: an Endocrine Society clinical practice guideline J Clin Endocrinol Metab 2011 96 1911 1930 10.1210/jc.2011-0385 21646368
18. Barreto DV Barreto FC Liabeuf S Temmar M Boitte F Choukroun G Vitamin D affects survival independently of vascular calcification in chronic kidney disease Clin J Am Soc Nephrol 2009 4 1128 1135 10.2215/CJN.00260109 19443628
19. Melamed ML Chonchol M Gutiérrez OM Kalantar-Zadeh K Kendrick J Norris K The Role of Vitamin D in CKD Stages 3 to 4: Report of a Scientific Workshop Sponsored by the National Kidney Foundation Am J Kidney Dis 2018 72 834 845 10.1053/j.ajkd.2018.06.031 30297082
20. Sutherland JP Zhou A Hyppönen E Vitamin D deficiency increases mortality risk in the UK Biobank : a nonlinear mendelian randomization study Ann Intern Med 2022 175 1552 1559 10.7326/M21-3324 36279545
21. Emerging Risk Factors Collaboration/EPIC-CVD/Vitamin D Studies Collaboration Estimating dose-response relationships for vitamin D with coronary heart disease, stroke, and all-cause mortality: observational and Mendelian randomisation analyses Lancet Diabetes Endocrinol 2021 9 837 46 10.1016/S2213-8587(21)00263-1 34717822
22 Zhang Y Fang F Tang J Jia L Feng Y Xu P Association between vitamin D supplementation and mortality: systematic review and meta-analysis BMJ 2019 366 l4673 10.1136/bmj.l4673 31405892
23. Scragg R Stewart AW Waayer D Lawes CMM Toop L Sluyter J Effect of monthly high-dose vitamin D supplementation on cardiovascular disease in the vitamin D assessment study: a randomized clinical trial JAMA Cardiol 2017 2 608 616 10.1001/jamacardio.2017.0175 28384800
24. Mann MC Hobbs AJ Hemmelgarn BR Roberts DJ Ahmed SB Rabi DM Effect of oral vitamin D analogs on mortality and cardiovascular outcomes among adults with chronic kidney disease: a meta-analysis Clin Kidney J 2015 8 41 48 10.1093/ckj/sfu122 25713709
25. Mousa H Elrayess MA Diboun I Jackson SK Zughaier SM Metabolomics profiling of vitamin D status in relation to dyslipidemia Metabolites 2022 12 771 10.3390/metabo12080771 36005643
26. Vimaleswaran KS Cavadino A Berry DJ LifeLines Cohort Study investigatorsJorde R Dieffenbach AK Association of vitamin D status with arterial blood pressure and hypertension risk: a mendelian randomisation study Lancet Diabetes Endocrinol 2014 2 719 29 10.1016/S2213-8587(14)70113-5 24974252
27. Song Y Wang L Pittas AG Del Gobbo LC Zhang C Manson JE Blood 25-hydroxy vitamin D levels and incident type 2 diabetes: a meta-analysis of prospective studies Diabetes Care 2013 36 1422 1428 10.2337/dc12-0962 23613602
28. de la Guía-Galipienso F Martínez-Ferran M Vallecillo N Lavie CJ Sanchis-Gomar F Pareja-Galeano H Vitamin D and cardiovascular health Clin Nutr 2021 40 2946 2957 10.1016/j.clnu.2020.12.025 33397599
29. Pilz S Verheyen N Gruebler MR Tomaschitz A Maerz W Vitamin D and cardiovascular disease prevention Nat Rev Cardiol 2016 13 404 417 10.1038/nrcardio.2016.73 27150190
30. Shen H Bielak LF Ferguson JF Streeten EA Yerges-Armstrong LM Liu J Association of the vitamin D metabolism gene CYP24A1 with coronary artery calcification Arterioscler Thromb Vasc Biol 2010 30 2648 2654 10.1161/ATVBAHA.110.211805 20847308
31. Li YC Kong J Wei M Chen Z-F Liu SQ Cao L-P 1,25-Dihydroxyvitamin D(3) is a negative endocrine regulator of the renin-angiotensin system J Clin Invest 2002 110 229 238 10.1172/JCI15219 12122115
32. Kim D-H Meza CA Clarke H Kim J-S Hickner RC Vitamin D and endothelial function Nutrients 2020 12 575 10.3390/nu12020575 32098418
33. Prietl B Treiber G Pieber TR Amrein K Vitamin D and immune function Nutrients 2013 5 2502 2521 10.3390/nu5072502 23857223
34. Oh J Weng S Felton SK Bhandare S Riek A Butler B 1,25(OH)2 vitamin d inhibits foam cell formation and suppresses macrophage cholesterol uptake in patients with type 2 diabetes mellitus Circulation 2009 120 687 698 10.1161/CIRCULATIONAHA.109.856070 19667238
35. Surdu AM Pînzariu O Ciobanu D-M Negru A-G Căinap S-S Lazea C Vitamin D and its role in the lipid metabolism and the development of atherosclerosis Biomedicines 2021 9 172 10.3390/biomedicines9020172 33572397
36. Lee JH O’Keefe JH Bell D Hensrud DD Holick MF Vitamin D deficiency an important, common, and easily treatable cardiovascular risk factor? J Am Coll Cardiol 2008 52 1949 1956 10.1016/j.jacc.2008.08.050 19055985
37. Kassi E Adamopoulos C Basdra EK Papavassiliou AG Role of vitamin D in atherosclerosis Circulation 2013 128 2517 2531 10.1161/CIRCULATIONAHA.113.002654 24297817
38. Khundmiri SJ Murray RD Lederer E PTH and Vitamin D Compr Physiol 2016 6 561 601 10.1002/cphy.c140071 27065162
39 Clarke R Shipley M Lewington S Youngman L Collins R Marmot M Underestimation of risk associations due to regression dilution in long-term follow-up of prospective studies Am J Epidemiol 1999 150 341 53 10.1093/oxfordjournals.aje.a010013 10453810