ORIGINAL RESEARCH article

Front. Nutr., 08 October 2025

Sec. Clinical Nutrition

Volume 12 - 2025 | https://doi.org/10.3389/fnut.2025.1597527

Associations of vitamin D with lipid metabolism, inflammation, and mortality vary by glycemic status and gender: a nationwide prospective study

  • 1. Clinical Nutrition Division, The Second Affiliated Hospital of Soochow University, Suzhou, China

  • 2. School of Public Health, Suzhou Medical College of Soochow University, Suzhou, China

Abstract

Background:

We aimed to investigate the associations of 25-hydroxyvitamin D (25OHD) with lipid, inflammatory, and injury biomarkers, and their potential mediating roles in vitamin D-related mortality.

Methods:

We analyzed data from 4,144 participants in the 2015–2018 National Health and Nutrition Examination Survey (NHANES). Glycemic status was classified using fasting plasma glucose and hemoglobin A1c levels. Key biomarkers were assessed, including visceral adiposity index (VAI) and systemic inflammation response index (SIRI). Multivariate linear regression and restricted cubic spline models examined associations.

Results:

In females, 25OHD was positively associated with HDL-cholesterol across all glycemic subgroups, and with triglycerides and total cholesterol in the normal glycemic subgroup (all p < 0.05). In males, inverse associations of 25OHD with triglycerides and VAI were most pronounced in the prediabetes group (p < 0.05). 25OHD was inversely associated with SIRI in normoglycemic individuals and prediabetic males. Mediation analyses revealed that SIRI partially mediated the inverse association between 25OHD and mortality only in prediabetic males, while lipid and injury biomarkers showed no significant mediation effects in any subgroup.

Conclusion:

Gender and glycemic status influence the associations between vitamin D and biomarkers, with inflammation potentially mediating the relationship between vitamin D and mortality in prediabetic males.

1 Introduction

As a classic fat-soluble secosteroid hormone, vitamin D plays a crucial role in calcium homeostasis and bone metabolism. Nevertheless, growing evidence over the past decades indicates that vitamin D’s influence extends far beyond skeletal health, encompassing a wide array of physiological processes, including glucose-lipid metabolism, inflammation response, and immune function (1, 2). Observational studies and meta-analyses have linked vitamin D deficiency to an increased risk of various non-communicable chronic diseases, including cardiovascular disease (3), type 2 diabetes (4), certain cancers (5), and even all-cause mortality (6). Consequently, maintaining optimal vitamin D levels has been proposed as a potential strategy for improving overall health and preventing chronic diseases (7).

Despite growing recognition of its multiple effects, the precise mechanisms by which vitamin D influences metabolic health remain incompletely understood. In particular, the extent to which vitamin D modulates distinct biological processes such as lipid metabolism, systemic inflammation, and cellular injury, as well as whether these factors mediate its effects on long-term health and mortality, remains an area of active investigation. A substantial body of evidence indicates that the effects of vitamin D on metabolic health are significantly modulated by both gender and glycemic status (8, 9). Therefore, a central, a priori hypothesis of our study was that these factors would act as critical effect modifiers, making a stratified analysis a necessary, hypothesis-driven component of our research design. A nuanced understanding of these specific interactions is particularly important, as the progression from normoglycemia to diabetes is tightly linked to disturbances in adiposity and systemic inflammation—two core components of overall metabolic health.

To better characterize these complex relationships, researchers are increasingly moving beyond single-marker analysis to gain a more holistic view of biological processes. This trend is evident across diverse fields, employing strategies that range from reference-free transcriptomic analyses in oncology (10) to the use of composite indices in metabolic research, such as the visceral adiposity index (VAI). The VAI, which integrates waist circumference (WC), body mass index (BMI), triglycerides (TG), and high-density lipoprotein cholesterol (HDL-C), provides a robust measure of metabolic dysfunction (11). Similarly, the systemic inflammation response index (SIRI) and systemic immune-inflammation index (SII), derived from complete blood counts, offer valuable insights into an individual’s systemic inflammatory burden (12). However, the extent to which vitamin D is associated with these indices—and whether these associations vary by gender and glycemic status—remains largely unexplored.

Therefore, in this study, we aimed to examine the associations between vitamin D status and lipid, inflammation, and cellular injury biomarkers, including VAI, SII, and SIRI. Additionally, we sought to determine whether these factors mediate the relationship between vitamin D and mortality risk, with a particular focus on gender- and glycemia-specific differences.

2 Materials and methods

2.1 Study design and participant selection

This study utilized baseline data from the National Health and Nutrition Examination Survey (NHANES) 2015–2018 cycles, with mortality outcomes ascertained prospectively through linkage to National Death Index (NDI) mortality files. This design allows for a temporal relationship to be established between baseline exposures and subsequent mortality. NHANES is a nationally representative program of studies designed to assess the health and nutritional status of adults and children in the United States, employing a complex, multistage probability sampling design (13). Data collection encompasses detailed questionnaires, standardized physical examinations, and extensive laboratory analyses. The flow chart of this study is presented in Figure 1. Initially, 6,227 participants aged over 12 years who provided blood samples were considered. Participants were excluded according to the following exclusion criteria: (1) absence of valid fasting glucose measurements (number = 773); (2) pregnancy status (number = 40); (3) missing serum 25-hydroxyvitamin D (25OHD) measurements (number = 25); (4) incomplete cardiometabolic or inflammatory biomarker data (number = 276); and (5) missing critical demographic or lifestyle information (age, race/ethnicity, anthropometric measurements, smoking status, alcohol consumption, or sedentary behavior; number = 969). Although the initial inclusion criterion allowed participants aged ≥ 12 years with available laboratory data, all individuals under 18 years of age were ultimately excluded due to missing data on key variables required for this analysis—such as fasting glucose, lipid and inflammatory markers, lifestyle factors, and mortality follow-up. As a result, the final analytic sample consisted of 4,144 adults aged 18 years and older.

Figure 1

2.2 Anthropometric and laboratory assessments

Body weight and height were measured by trained health technicians using standardized protocols. BMI was calculated as weight in kilograms divided by height in meters squared. WC was measured at the level just above the iliac crest using a flexible tape, with participants standing and breathing normally. Blood samples were collected following a standardized protocol and analyzed at laboratories certified by the Centers for Disease Control and Prevention. All assays were performed using validated enzymatic or hematological methods under rigorous quality control, as detailed in the NHANES Laboratory Procedures Manuals. Serum triglycerides and HDL-C were measured enzymatically using the Roche Modular P chemistry analyzer. The same platform was also used to measure injury markers including creatine phosphokinase (CPK), lactate dehydrogenase (LDH), alanine aminotransferase (ALT), aspartate aminotransferase (AST). Complete blood counts, including neutrophils, lymphocytes, monocytes, and platelets, were assessed using the Beckman Coulter DxH 800 hematology analyzer. All assays followed rigorous quality control and quality assurance procedures as outlined in the NHANES Laboratory Procedures Manuals. Full details of these examination and laboratory procedures are publicly available at https://www.cdc.gov/nchs/nhanes/ (Accessed on 29 August 2025).

2.3 Outcome assessment

Mortality status was determined by linkage of NHANES participants to the NDI through probabilistic matching using a combination of personal identifiers (e.g., Social Security number, name, date of birth). The publicly available Linked Mortality File (LMF) used in this study, which provides data through December 31, 2019, can be accessed at http://www.cdc.gov/nchs/linked-data/mortality-files/index.html (Accessed on 29 August 2025). From this file, we identified all-cause mortality as the primary outcome. The follow-up time was calculated from the date of the NHANES examination to the date of death or the end of the follow-up period (December 31, 2019), whichever came first.

2.4 Definition of diabetes status

Diabetes status was determined according to the American Diabetes Association criteria (14): diabetes was defined as a fasting plasma glucose (FPG) level ≥ 7.0 mmol/L or hemoglobin A1C (HbA1c) ≥ 6.5%; prediabetes was defined as 5.6 mmol/L ≤ FPG < 7.0 mmol/L or 5.7% ≤ HbA1c < 6.5%; and normoglycemic status was defined as FPG < 5.6 mmol/L or HbA1c < 5.7%.

2.5 Calculation of Indices

The following composite indices were calculated based on previously established formulas. SII, SIRI, and VAI are calculated as follows: SII = (platelet count × neutrophil count)/lymphocyte count (15); SIRI = (neutrophil count × monocyte count)/lymphocyte count (16); VAI (male) = (WC/[39.68 + 1.88 × BMI]) × (TG/1.03) × (1.31/HDL-C); VAI (female) = (WC/[36.58 + 1.89 × BMI]) × (TG/0.81) × (1.52/HDL-C) (17).

2.6 Statistical analysis

Baseline participant characteristics were summarized using descriptive statistics. Data are presented as numbers (percentages) for categorical variables or medians (interquartile ranges) for continuous variables. Categorical variables were compared using the chi-square (χ2) test, while continuous variables were compared using the Kruskal–Wallis H test.

2.6.1 Multivariate linear regression

Associations between serum 25OHD levels and lipid, inflammatory, and cellular injury biomarkers were evaluated using multivariate linear regression models. The general form of the model is:

where Y denotes the dependent biomarker, X25OHD is the serum 25OHD level, Xi (i = 2, …, p) are adjustment covariates, βᵢ are the corresponding regression coefficients, p is the total number of predictors in the model, and ϵ is the random error term. Model 1 was adjusted for age, race/ethnicity, and BMI, and model 2 additionally adjusted for alcohol consumption, smoking status, and sedentary activity.

2.6.2 Restricted cubic splines (RCS)

To assess the potential nonlinear association between serum 25OHD and mortality, we fitted RCS within a Cox proportional hazards model:

where s(X) denotes the spline function of serum 25OHD, Zj is the j-th covariate, γj its coefficient, and q is the number of covariates. This corresponds to the proportional hazards form h(t∣X, Z) = h0 (t) exp.{s(X) + ∑γjZj}, where h0 (t) is the baseline hazard. We used three knots placed at the 10th, 50th, and 90th percentiles of the 25OHD distribution. Nonlinearity was tested using likelihood ratio tests, and the proportional hazards assumption was verified using Schoenfeld residuals.

2.6.3 Mediation analyses

To investigate whether biomarkers mediated the effect of 25OHD (X) on mortality (Y), we applied a counterfactual-based mediation framework. The mediator (M) model was:

and the outcome model for time-to-event mortality was:

where T is the survival time, Wⱼ denotes the j-th covariate (indexed from 1 to r, the total number of covariates); a quantifies the effect of X on M; b reflects the effect of M on the log-survival time given X; and c’ represents the direct effect of X on the log-survival time given M. The term σεY represents the error, where σ is the scale parameter of the survival distribution. Based on these models, which were fitted using a Weibull distribution in the survreg function, the average causal mediation effect (ACME) and average direct effect (ADE) were estimated using the mediation package in R with quasi-Bayesian simulations (1,000 iterations). The analysis allowed for an exposure-mediator interaction.

All statistical analyses were performed using R software (version 4.3.3). A two-sided p-value < 0.05 was considered statistically significant.

3 Results

3.1 Baseline characteristics of participants

As shown in Table 1, the study population comprised 1,391 normoglycemic, 2020 prediabetic, and 733 diabetic individuals. The median age increased significantly from 35.0 years in the normoglycemic group to 63.0 years in the diabetic group (p < 0.001). The proportion of females was significantly higher in the normoglycemic group (59.7%) compared with the diabetic group (46.2%; p < 0.001). Racial/ethnic composition also differed significantly, with a higher prevalence of non-Hispanic Whites among normoglycemic and prediabetic participants. Additionally, smoking prevalence was greater among diabetic participants (47.3%) compared with normoglycemic individuals (36.6%; p < 0.001), and BMI increased progressively across groups. Triglyceride was highest, while HDL-C was lowest in the diabetic group. Other laboratory indices—including SIRI, LDH, and ALT—were significantly elevated, whereas serum 25(OH)D levels were highest in the diabetic group.

Table 1

VariableNormalPrediabetesDiabetesp-value
Number1,3912,020733
Age, years35.0 (25.0–51.0)54.0 (39.0–66.0)63.0 (54.0–71.0)<0.001
Gender
 Female830 (59.7)948 (47.0)339 (46.2)<0.001
 Male561 (40.3)1,072 (53.1)394 (52.8)
Race/Ethnicity
 Mexican American188 (13.5)345 (17.1)137 (18.7)0.006
 Other Hispanic150 (10.8)228 (11.3)97 (13.2)
 Non-Hispanic White516 (37.1)705 (34.9)224 (30.6)
 Non-Hispanic Black291 (20.9)418 (20.7)164 (23.4)
 Other (including multiracial)246 (17.7)324 (16.0)111 (15.1)
Alcohol consumption
 Yes1,096 (78.8)1,599 (79.2)560 (76.4)0.286
 No295 (21.2)421 (20.8)173 (23.6)
Smoking status
 Yes509 (36.6)885 (43.8)347 (47.3)<0.001
 No882 (63.4)1,135 (56.2)386 (52.7)
Body mass index, kg/m226.0 (22.6–30.6)28.9 (25.2–33.7)31.2 (27.2–36.3)<0.001
Sedentary activity, minute300 (180–480)300 (180–480)300 (180–480)0.176
Triglycerides, mmol/L0.8 (0.6–1.2)1.1 (0.7–1.6)1.4 (1.0–1.9)<0.001
Total Cholesterol, mmol/L4.6 (4.0–5.3)4.9 (4.2–5.6)4.5 (3.9–5.3)<0.001
HDL-C, mmol/L1.4 (1.2–1.8)1.3 (1.1–1.6)1.2 (1.0–1.4)<0.001
LDL-C, mmol/L2.7 (2.2–3.3)2.9 (2.4–3.5)2.6 (1.9–3.2)<0.001
VAI0.9 (0.6–1.5)1.3 (0.8–2.2)1.9 (1.2–3.0)<0.001
RDW, %13.3 (12.9–13.9)13.6 (13.1–14.3)13.8 (13.2–14.6)<0.001
hs-CRP, mg/L1.34 (0.6–3.2)2.1 (0.9–4.8)2.9 (1.3–6.2)<0.001
SIRI0.87 (0.58–1.58)1.00 (0.65–1.47)1.19 (0.77–1.76)<0.001
SII406.30 (291.00–569.25)427.83 (303.61–609.71)459.00 (328.07–639.95)<0.001
CPK, IU/L107.0 (74.0–170.0)115.0 (75.0–181.0)101.0 (67.5–168.0)<0.001
LDH, IU/L133.0 (117.0–155.0)142.0 (123.0–164.0)140.0 (122.0–163.0)<0.001
ALT, U/L17.0 (13.0–24.0)20.0 (15.0–28.0)22.0 (16.0–31.0)<0.001
AST, U/L21.0 (17.0–25.0)21.0 (18.0–26.0)21.0 (17.0–27.0)<0.001
25OHD, nmol/L61.8 (44.8–78.9)63.3 (46.3–83.1)66.4 (45.8–87.3)<0.001

The baseline characteristics of NHANES 2015–2018 participants grouped by diabetes status.

Values are presented as count (percentage) or median (25–75 interquartile range). SII is calculated as the product of neutrophil count and platelet count divided by lymphocyte count. SIRI is calculated as the product of neutrophil count and monocyte count divided by lymphocyte count. 25OHD, 25-hydroxyvitamin D; ALT, alanine aminotransferase; AST, aspartate aminotransferase; CPK, creatine phosphokinase; HDL-C, high-density lipoprotein-cholesterol; hs-CRP, high-sensitivity C-reactive protein; LDH, lactate dehydrogenase; LDL-C, low-density lipoprotein-cholesterol; RDW, red cell distribution width; SII, systemic immune-inflammation index; SIRI, systemic inflammation response index; VAI, visceral adiposity index.

3.2 Associations of 25OHD with lipid, inflammation, and injury biomarkers

The associations between vitamin D and various lipid parameters were modified by both diabetes status and gender (Table 2). In normoglycemic females, 25OHD was positively associated with triglycerides, total cholesterol, and HDL-C (all p < 0.05). In prediabetes, 25OHD was positively associated with HDL-C in both genders, while inverse associations with triglycerides and VAI were evident only in males. Among diabetic participants, a consistent positive association between 25OHD and HDL-C was observed in females, with no significant associations in other lipid measures for either gender.

Table 2

Diabetes statusGenderVariableCrude modelAdjusted model 1Adjusted model 2
βPβ1P1β2P2
NormalFemaleTriglyceride, mmol/L0.002<0.0010.0010.0460.0010.030
Total Cholesterol, mmol/L0.007<0.0010.0030.0270.0030.030
HDL-C, mmol/L0.004<0.0010.002<0.0010.002<0.001
LDL-C, mmol/L0.0020.0150.0000.9550.0000.907
VAI0.0010.3550.0010.3310.0010.286
MaleTriglyceride, mmol/L0.0000.6870.0000.8640.0000.950
Total Cholesterol, mmol/L0.0010.650−0.0020.234−0.0020.206
HDL-C, mmol/L0.0000.8330.0000.6250.0000.573
LDL-C, mmol/L0.0010.633−0.0020.299−0.0020.263
VAI<0.0010.8440.0010.4360.0020.372
PrediabetesFemaleTriglyceride, mmol/L0.0010.1790.0000.5090.0000.634
Total Cholesterol, mmol/L0.004<0.0010.0000.9360.0000.882
HDL-C, mmol/L0.003<0.0010.0010.0040.0010.006
LDL-C, mmol/L0.0000.626−0.0010.339−0.0010.389
VAI−0.0010.306<0.0010.785<0.0010.811
MaleTriglyceride, mmol/L−0.0020.016−0.0030.003−0.0030.004
Total Cholesterol, mmol/L−0.0020.138−0.0010.653−0.0010.640
HDL-C, mmol/L0.002<0.0010.002<0.0010.002<0.001
LDL-C, mmol/L−0.0030.002−0.0010.230−0.0010.232
VAI−0.005<0.001−0.0040.003−0.0040.004
DiabetesFemaleTriglyceride, mmol/L0.0000.7700.0000.9840.0000.893
Total Cholesterol, mmol/L0.0010.7840.0020.3730.0020.375
HDL-C, mmol/L0.003<0.0010.0020.0020.0020.002
LDL-C, mmol/L−0.0020.230<0.0010.8740.0000.858
VAI−0.0040.097−0.0010.703−0.0010.798
MaleTriglyceride, mmol/L−0.0010.447−0.0010.620−0.0010.604
Total Cholesterol, mmol/L−0.0050.013−0.0010.620−0.0010.604
HDL-C, mmol/L0.0010.1710.0000.6810.0000.664
LDL-C, mmol/L−0.0050.003−0.0010.707−0.0010.687
VAI−0.0030.3920.0010.8450.0010.872

The association of serum 25OHD levels with lipid biomarkers in NHANES 2015–2018 participants.

The regression coefficients β and p-values are derived from linear regression analysis. Model 1, adjusted for age, race/ethnicity, and body mass index; model 2, adjusted for age, race/ethnicity, body mass index, alcohol consumption, smoking status, and minutes sedentary activity. 25OHD, 25-hydroxyvitamin D; HDL-C, high-density lipoprotein-cholesterol; LDL-C, low-density lipoprotein-cholesterol; VAI, visceral adiposity index. For models with VAI as the outcome, BMI was not included as an adjustment variable to avoid collinearity.

As presented in Table 3, in normoglycemic individuals, 25OHD was inversely associated with SIRI in females, and with both SII and SIRI in males after adjustment (all p < 0.05), while no significant association was found with high-sensitivity C-reactive protein (hs-CRP) or red cell distribution width (RDW) in certain subgroups. Similarly, in males, the serum 25OHD level was significantly and inversely associated with SII and SIRI after adjustment. In prediabetic males, the inverse association with SIRI remained significant (p = 0.032), whereas such associations were not observed in females. In the diabetic group, initial positive associations in males were attenuated after adjustment.

Table 3

Diabetes statusGenderVariableCrude modelAdjusted model 1Adjusted model 2
βPβ1P1β2P2
NormalFemaleRDW, %−0.0050.002−0.0020.211−0.0020.199
hs-CRP, mg/L−0.0130.043−0.0050.519−0.0040.546
SIRI0.0000.499−0.0020.006−0.0020.008
SII−0.2820.294−0.3350.281−0.3270.293
MaleRDW, %0.0000.875−0.0010.386−0.0010.417
hs-CRP, mg/L−0.0080.415−0.0120.278−0.0120.287
SIRI0.0010.340−0.0040.006−0.0050.005
SII−0.2980.464−1.0480.019−1.1420.011
PrediabetesFemaleRDW, %−0.007<0.001−0.0020.139−0.0020.137
hs-CRP, mg/L−0.035<0.001−0.0160.112−0.0150.136
SIRI0.0020.0210.0000.8790.0000.636
SII0.3510.3200.1580.6950.2440.543
MaleRDW, %0.0010.4320.0000.7190.0000.830
hs-CRP, mg/L0.0010.876−0.0040.574−0.0030.747
SIRI0.0020.105−0.0030.016−0.0030.032
SII0.6370.050−0.2480.488−0.2060.564
DiabetesFemaleRDW, %−0.0030.182−0.0010.579−0.0020.479
hs-CRP, mg/L−0.0420.075−0.0190.434−0.0220.361
SIRI−0.0010.446−0.0030.086−0.0030.074
SII−0.5280.357−0.7600.213−0.7860.199
MaleRDW, %−0.0010.667−0.0030.233−0.0030.242
hs-CRP, mg/L−0.0110.5260.0030.8910.0030.887
SIRI0.0050.0030.0000.8380.0000.830
SII1.0570.0380.2530.6520.2440.662

The association of serum 25OHD levels with inflammation biomarkers in NHANES 2015–2018 participants.

The regression coefficients β and p-values are derived from linear regression analysis. Model 1, adjusted for age, race/ethnicity, and body mass index; model 2, adjusted for age, race/ethnicity, body mass index, alcohol consumption, smoking status, and minutes sedentary activity. hs-CRP, high-sensitivity C-reactive protein; RDW, red cell distribution width; SII, systemic immune-inflammation index; SIRI, systemic inflammation response index.

The association between 25OHD and cellular injury is demonstrated in Table 4. Normoglycemic females exhibited a positive association between 25OHD and AST (p = 0.005). Normoglycemic males showed an inverse association between 25OHD and ALT in the adjusted model 1 (p = 0.042), however, this association lost significance upon further adjustment in the adjusted model 2 (p = 0.052). In prediabetic males and diabetic individuals, 25OHD was inversely associated with CPK in crude analyses, but these associations were not maintained after adjustment.

Table 4

Diabetes statusGenderVariableCrude modelAdjusted model 1Adjusted model 2
βPβ1P1β2P2
NormalFemaleCPK, IU/L−0.1290.3960.1880.2860.1770.317
LDH, IU/L0.0990.0080.0360.3830.0240.560
ALT, U/L0.0350.0520.0360.0930.0360.093
AST, U/L0.042<0.0010.0330.0070.0340.005
MaleCPK, IU/L−0.3050.4140.8670.0330.7880.054
LDH, IU/L0.0190.7370.0140.819−0.0350.564
ALT, U/L−0.0430.101−0.0580.042−0.0550.052
AST, U/L0.0080.640−0.0100.611−0.0100.589
PrediabetesFemaleCPK, IU/L−0.1360.0700.0910.2680.0760.352
LDH, IU/L0.132<0.0010.0810.0560.0560.070
ALT, U/L0.0000.9920.0230.1990.0250.174
AST, U/L0.0080.7730.0050.8660.0090.763
MaleCPK, IU/L−1.8840.003−0.6280.381−0.6470.368
LDH, IU/L0.0580.1380.0190.6620.0080.850
ALT, U/L−0.083<0.001−0.0310.120−0.0300.132
AST, U/L−0.0340.018−0.0210.189−0.0200.220
DiabetesFemaleCPK, IU/L−0.3210.029−0.2180.161−0.2220.154
LDH, IU/L0.0050.933−0.0470.421−0.0570.322
ALT, U/L0.0030.9170.0140.5910.0160.520
AST, U/L−0.0030.903−0.0100.665−0.0090.709
MaleCPK, IU/L−0.6820.0180.1120.7110.1070.725
LDH, IU/L0.0750.2070.0840.2040.0860.192
ALT, U/L−0.0510.1060.0230.5010.0220.526
AST, U/L−0.0030.9000.0160.5100.0160.522

The association of serum 25OHD levels with cellular injury biomarkers in NHANES 2015–2018 participants.

The regression coefficients β and p-values are derived from linear regression analysis. Model 1, adjusted for age, race/ethnicity, and body mass index; model 2, adjusted for age, race/ethnicity, body mass index, alcohol consumption, smoking status, and minutes sedentary activity. ALT, alanine aminotransferase; AST, aspartate aminotransferase; CPK, creatine phosphokinase; LDH, lactate dehydrogenase.

3.3 Vitamin D and total mortality

Figure 2 presents the non-linear associations between serum 25OHD levels and total mortality across different diabetes statuses and genders. Among normoglycemic, prediabetic, and female diabetic participants, mortality risk declined with increasing serum 25OHD levels up to a certain threshold, after which the association plateaued. Notably, a statistically significant inverse association was observed only in male prediabetic participants. Although an inverted U-shaped relationship was observed in male diabetic participants, this association did not reach statistical significance.

Figure 2

3.4 Mediation analyses of lipids, inflammation, and injury

The results of mediation analyses are summarized in Tables 5–7. In normoglycemic and diabetic participants, lipid biomarkers did not significantly mediate the relationship between 25OHD and mortality in either gender. Among prediabetics, significant ADEs for total cholesterol, HDL-C, and low-density lipoprotein-cholesterol (LDL-C) were observed in females. For inflammatory markers, only prediabetic males exhibited a significant ACME for SIRI (p = 0.042), with no significant mediation observed in other groups. Similarly, mediation analyses for cellular injury markers revealed no significant mediation effects in normoglycemic or diabetic participants, with only prediabetic females showing significant ADEs in LDH, ALT, and AST.

Table 5

Diabetes statusGenderLipid biomarkersADE (95%CI)p-valueACME (95%CI)p-valueTotal effect (95%CI)p-value
NormalFemaleTriglyceride, mmol/L−1.18e-04 (−5.09e-04–0.00)0.510−4.84e-06 (−2.17e-05–0.00)0.490−1.23e-04 (−5.29e-04–0.00)0.500
Total Cholesterol, mmol/L−1.18e-04 (−5.11e-04–0.00)0.520−4.65e-06 (−4.41e-05–0.00)0.790−1.23e-04 (−5.14e-04–0.00)0.500
HDL-C, mmol/L−1.06e-04 (−5.13e-04–0.00)0.640−1.66e-05 (−7.12e-05–0.00)0.380−1.23e-04 (−5.20e-04–0.00)0.590
LDL-C, mmol/L−1.23e-04 (−5.11e-04–0.00)0.490−5.70e-08 (−1.42e-05–0.00)0.990−1.23e-04 (−5.15e-04–0.00)0.490
VAI−1.71e-04 (−5.59e-04–0.00)0.410−1.87e-06 (−9.79e-06–0.00)0.570−1.73e-04 (−5.62e-04–0.00)0.400
MaleTriglyceride, mmol/L−1.03e-04 (−8.73e-04–0.00)0.700−5.12e-07 (−3.09e-05–0.00)0.930−1.03e-04 (−8.83e-04 –0.00)0.700
Total Cholesterol, mmol/L−1.15e-04 (−8.83e-04–0.00)0.7801.15e-05 (−1.73e-05–0.00)0.450−1.03e-04 (−8.88e-04–0.00)0.800
HDL-C, mmol/L−1.02e-04 (−8.84e-04–0.00)0.780−9.38e-07 (−6.89e-05–0.00)0.960−1.03e-04 (−8.75e-04–0.00)0.800
LDL-C, mmol/L−1.18e-04 (−9.67e-04–0.00)0.6901.50e-05 (−1.54e-05–0.00)0.410−1.03e-04 (−9.24e-04–0.00)0.720
VAI−8.96e-05 (−8.56e-04–0.00)0.740−2.28e-06 (−3.17e-05–0.00)0.850−9.19e-05 (−8.56e-04–0.00)0.750
PrediabetesFemaleTriglyceride, mmol/L−4.60e-04 (−9.04e-04–0.00)0.0568.97e-06 (−2.65e-05–0.00)0.602−4.51e-04 (−9.10e-04–0.00)0.060
Total Cholesterol, mmol/L−4.48e-04 (−9.33e-04–0.00)0.040−2.89e-06 (−4.49e-05–0.00)0.880−4.51e-04 (−9.25e-04–0.00)0.044
HDL-C, mmol/L−4.68e-04 (−9.15e-04–0.00)0.0241.77e-05 (−1.73e-05–0.00)0.370−4.51e-04 (−8.99e-04–0.00)0.034
LDL-C, mmol/L−4.67e-04 (−9.38e-04–0.00)0.0381.62e-05 (−2.04e-05–0.00)0.418−4.51e-04 (−9.19e-04–0.00)0.040
VAI−3.06e-04 (−6.78e-04–0.00)0.0801.94e-05 (−4.28e-06–0.00)0.120−2.87e-04 (−6.49e-04–0.00)0.100
MaleTriglyceride, mmol/L−3.31e-04 (−9.76e-04–0.00)0.3201.90e-05 (−3.30e-05–0.00)0.290−3.12e-04 (−9.54e-04–0.00)0.330
Total Cholesterol, mmol/L−3.21e-04 (−9.90e-04–0.00)0.3409.29e-06 (−2.73e-05–0.00)0.660−3.12e-04 (−9.69e-04–0.00)0.350
HDL-C, mmol/L−3.09e-04 (−9.59e-04–0.00)0.370−2.53e-06 (−9.26e-05–0.00)0.920−3.12e-04 (−9.48e-04–0.00)0.350
LDL-C, mmol/L−3.36e-04 (−1.02e-03–0.00)0.3202.40e-05 (−1.41e-05–0.00)0.220−3.12e-04 (−1.01e-03–0.00)0.350
VAI−2.90e-04 (−9.86e-04–0.00)0.3606.42e-06 (−4.48e-05–0.00)0.720−2.83e-04 (−9.64e-04–0.00)0.370
DiabetesFemaleTriglyceride, mmol/L−4.65e-04 (−1.13e-03–0.00)0.230−5.95e-06 (−1.13e-04–0.00)0.920−4.70e-04 (−1.19e-03–0.00)0.240
Total Cholesterol, mmol/L−4.74e-04 (−1.24e-03–0.00)0.2203.17e-06 (−4.76e-05–0.00)0.840−4.70e-04 (−1.24e-03–0.00)0.230
HDL-C, mmol/L−6.71e-04 (−1.45e-03–0.00)0.0802.00e-04 (−4.33e-05–0.00)0.200−4.70e-04 (−1.27e-03–0.00)0.210
LDL-C, mmol/L−4.72e-04 (−1.23e-03–0.00)0.2501.42e-06 (−4.40e-05–0.00)0.970−4.70e-04 (−1.21e-03–0.00)0.250
VAI−5.49e-04 (−1.40e-03–0.00)0.1602.94e-05 (−3.28e-05–0.00)0.350−5.19e-04 (−1.39e-03–0.00)0.180
MaleTriglyceride, mmol/L−3.33e-04 (−1.13e-03–0.00)0.4704.55e-06 (−4.94e-05–0.00)0.890−3.29e-04 (−1.13e-03–0.00)0.470
Total Cholesterol, mmol/L−3.28e-04 (−1.18e-03–0.00)0.410−9.97e-07 (−5.45e-05–0.00)0.960−3.29e-04 (−1.19e-03–0.00)0.400
HDL-C, mmol/L−3.36e-04 (−1.12e-03–0.00)0.4206.83e-06 (−4.72e-05–0.00)0.800−3.29e-04 (−1.10e-03–0.00)0.430
LDL-C, mmol/L−3.29e-04 (−1.15e-03–0.00)0.440−2.81e-07 (−5.58e-05–0.00)1.000−3.29e-04 (−1.17e-03–0.00)0.440
VAI−4.09e-04 (−1.22e-03–0.00)0.360−1.10e-06 (−4.38e-05–0.00)0.940−4.10e-04 (−1.23e-03–0.00)0.350

The mediation effect of lipid biomarkers in the association between 25OHD and mortality in NHANES 2015–2018 participants.

25OHD, 25-hydroxyvitamin D; ACME, average causal mediated effect; ADE, average direct effect; CI, confidence interval; HDL-C, high-density lipoprotein-cholesterol; LDL-C, low-density lipoprotein-cholesterol; VAI, visceral adiposity index.

Table 6

Diabetes statusGenderInflammation biomarkersADE (95%CI)p-valueACME (95%CI)p-valueTotal Effect (95%CI)p-value
NormalFemaleRDW, %−1.24e-04 (−4.99e-04–0.00)0.4808.94e-07 (−2.46e-05–0.00)0.920−1.23e-04 (−4.99e-04–0.00)0.490
hs-CRP, mg/L−1.26e-04 (−5.15e-04–0.00)0.5102.75e-06 (−5.92e-06–0.00)0.520−1.23e-04 (−5.15e-04–0.00)0.530
SIRI−1.02e-04 (−4.80e-04–0.00)0.640−2.09e-05 (−7.20e-05–0.00)0.250−1.23e-04 (−5.01e-04–0.00)0.580
SII−1.24e-04 (−5.03e-04–0.00)0.5201.48e-06 (−8.69e-06–0.00)0.820−1.23e-04 (−5.00e-04–0.00)0.520
MaleRDW, %−9.76e-05 (−9.12e-04–0.00)0.740−1.27e-05 (−7.15e-05–0.00)0.870−1.03e-04(−9.27e-04–0.00)0.730
hs-CRP, mg/L−9.05e-05 (−8.43e-04–0.00)0.810−7.21e-05 (−2.69e-04–0.00)0.380−1.03e-04(−8.65e-04–0.00)0.790
SIRI−1.43e-04 (−8.93e-04–0.00)0.6304.02e-05 (4.02e-05–0.00)0.280−1.03e-04 (−8.56e-04–0.00)0.690
SII−1.60e-04 (−1.00e-03–0.00)0.6105.70e-05 (−1.29e-05–0.00)0.150−1.03e-04 (−9.44e-04–0.00)0.690
PrediabetesFemaleRDW, %−4.32e-04 (−9.22e-04–0.00)0.082−1.85e-05 (−5.03e-05–0.00)0.158−4.51e-04 (−9.48e-04–0.00)0.074
hs-CRP, mg/L−4.50e-04 (−8.92e-04–0.00)0.036−4.38e-07 (−2.03e-05–0.00)0.920−4.51e-04 (−9.02e-04–0.00)0.038
SIRI−4.63e-04 (−9.07e-04–0.00)0.0521.26e-05 (−5.79e-05–0.00)0.720−4.51e-04 (−8.88e-04–0.00)0.054
SII−4.64e-04 (−9.11e-04–0.00)0.0401.36e-05 (−7.92e-05–0.00)0.802−4.51e-04 (−9.06e-04–0.00)0.042
MaleRDW, %−3.05e-04 (−9.64e-04–0.00)0.370−6.67e-06 (−1.20e-04–0.00)0.780−3.12e-04(−9.77e-04–0.00)0.350
hs-CRP, mg/L−3.02e-04 (−9.60e-04–0.00)0.380−9.75e-06 (−8.63e-05–0.00)0.710−3.12e-04 (−9.93e-04–0.00)0.380
SIRI−2.42e-04 (−8.53e-04–0.00)0.462−7.01e-05 (−2.21e-04–0.00)0.042−3.12e-04 (−9.63e-04–0.00)0.366
SII−3.01e-04 (−9.57e-04–0.00)0.340−1.06e-05 (−7.36e-05–0.00)0.580−3.12e-04 (−9.72e-04–0.00)0.320
DiabetesFemaleRDW, %−4.28e-04 (−1.17e-03–0.00)0.270−4.23e-05 (−1.73e-04–0.00)0.450−4.70e-04 (−1.22e-03–0.00)0.230
hs-CRP, mg/L−4.55e-04 (−1.23e-03–0.00)0.240−1.55e-05 (−6.48e-05–0.00)0.460−4.70e-04 (−1.23e-03–0.00)0.220
SIRI−4.29e-04 (−1.18e-03–0.00)0.280−4.13e-05 (−1.28e-04–0.00)0.290−4.70e-04 (−1.20e-03–0.00)0.220
SII−4.63e-04 (−1.19e-03–0.00)0.240−7.52e-06 (−6.30e-05–0.00)0.830−4.70e-04 (−1.19e-03–0.00)0.210
MaleRDW, %−2.57e-04 (−1.12e-03–0.00)0.500−7.21e-05 (−2.69e-04–0.00)0.340−3.29e-04 (−1.21e-03–0.00)0.380
hs-CRP, mg/L−3.44e-04 (−1.13e-03–0.00)0.4401.47e-05 (−2.05e-04–0.00)0.940−3.29e-04 (−1.11e-03–0.00)0.420
SIRI−3.08e-04 (−1.10e-03–0.00)0.430−2.12e-05 (−2.66e-04–0.00)0.790−3.29e-04 (−1.17e-03–0.00)0.400
SII−3.56e-04 (−1.17e-03–0.00)0.430−2.7e-05 (−1.17e-04–0.00)0.700−3.29e-04 (−1.17e-03–0.00)0.450

The mediation effect of inflammation biomarkers in the association between 25OHD and mortality in NHANES 2015–2018 participants.

25OHD, 25-hydroxyvitamin D; ACME, average causal mediated effect; ADE, average direct effect; CI, confidence interval; HDL-C, high-density lipoprotein-cholesterol; LDL-C, low-density lipoprotein-cholesterol; RDW, red cell distribution width; SII, systemic immune-inflammation index; SIRI, systemic inflammation response index.

Table 7

Diabetes statusGenderVariableADE of Vit D and 95% CIp-valueACME of variablesp-valueTotal effectp-value
NormalFemaleCPK, IU/L−1.23e-04 (−4.86e-04–0.00)0.5104.21e-07 (−1.18e-05–0.00)0.950−1.23e-04 (−4.84e-04–0.00)0.500
LDH, IU/L−1.23e-04 (−5.04e-04–0.00)0.4905.89e-07 (−6.44e-06–0.00)0.900−1.23e-04 (−5.01e-04–0.00)0.490
ALT, U/L−1.19e-04 (−5.19e-04–0.00)0.520−3.52e-06 (−4.62e-05–0.00)0.310−1.23e-04 (−5.42e-04–0.00)0.490
AST, U/L−1.24e-04 (−5.00e-04–0.00)0.5501.46e-06 (−3.44e-05–0.00)0.980−1.23e-04 (−4.85e-04–0.00)0.540
MaleCPK, IU/L−9.47e-05 (−8.47e-04–0.00)0.780−8.55e-06 (−6.00e-05–0.00)0.660−1.03e-04 (−8.44e-04–0.00)0.760
LDH, IU/L−9.82e-05 (−8.61e-04–0.00)0.700−5.00e-06 (−4.05e-05–0.00)0.810−1.03e-04 (−8.60e-04–0.00)0.700
ALT, U/L−9.12e-05 (−9.04e-04–0.00)0.790−1.20e-05 (−1.04e-04–0.00)0.700−1.03e-04 (−9.10e-04–0.00)0.760
AST, U/L−7.80e-05 (−8.62e-04–0.00)0.860−2.52e-05 (−1.75e-04–0.00)0.590−1.03e-04 (−8.78e-04–0.00)0.810
PrediabetesFemaleCPK, IU/L−4.58e-04 (−8.99e-04–0.00)0.0587.14e-06 (−9.05e-06–0.00)0.516−4.51e-04 (−8.92e-04–0.00)0.064
LDH, IU/L−4.65e-04 (−9.32e-04–0.00)0.0421.44e-05 (−8.05e-06–0.00)0.242−4.51e-04 (−9.16e-04–0.00)0.048
ALT, U/L−4.49e-04 (−8.96e-04–0.00)0.046−2.14e-06 (−1.75e-05–0.00)0.728−4.51e-04 (−9.03e-04–0.00)0.048
AST, U/L−4.51e-04 (−9.14e-04–0.00)0.0361.07e-07 (−3.48e-06–0.00)0.738−4.51e-04 (−9.10e-04–0.00)0.040
MaleCPK, IU/L−3.11e-04 (−9.79e-04–0.00)0.370−8.44e-07 (−5.56e-06–0.00)0.690−3.12e-04 (−9.79e-04–0.00)0.370
LDH, IU/L−3.16e-04 (−9.88e-04–0.00)0.3303.73e-06 (−9.78e-05–0.00)0.960−3.12e-04 (−9.84e-04–0.00)0.320
ALT, U/L−2.88e-04 (−9.62e-04–0.00)0.410−2.43e-05 (−8.09e-05–0.00)0.230−3.12e-04 (−9.92e-04–0.00)0.380
AST, U/L−2.83e-04 (−9.38e-04–0.00)0.390−2.86e-05 (−1.09e-04–0.00)0.300−3.12e-04 (−9.93e-04–0.00)0.370
DiabetesFemaleCPK, IU/L−4.69e-04 (−1.21e-03–0.00)0.190−1.91e-06 (−4.97e-05–0.00)0.920−4.70e-04 (−1.22e-03–0.00)0.190
LDH, IU/L−4.42e-04 (−1.16e-03–0.00)0.240−2.89e-05 (−1.51e-04–0.00)0.490−4.70e-04 (−1.20e-03–0.00)0.220
ALT, U/L−4.66e-04 (−1.17e-03–0.00)0.220−4.45e-06 (−3.86e-05–0.00)0.750−4.70e-04 (−1.18e-03–0.00)0.220
AST, U/L−4.67e-04 (−1.20e-03–0.00)0.240−3.43e-06 (−2.62e-05–0.00)0.990−4.70e-04 (−1.19e-03–0.00)0.230
MaleCPK, IU/L−3.14e-04 (−1.14e-03–0.00)0.520−1.51e-05 (−1.52e-04–0.00)0.840−3.29e-04 (−1.15e-03–0.00)0.500
LDH, IU/L−3.12e-04 (−1.14e-03–0.00)0.430−1.7e-05 (−1.04e-04–0.00)0.650−3.29e-04 (−1.12e-03–0.00)0.390
ALT, U/L−3.20e-04 (−1.23e-03–0.00)0.440−9.29e-06 (−7.32e-05–0.00)0.620−3.29e-04 (−1.22e-03–0.00)0.430
AST, U/L−3.20e-04 (−1.18e-03–0.00)0.440−9.30e-06 (−7.15e-05–0.00)0.710−3.29e-04 (−1.17e-03–0.00)0.420

The mediation effect of injury biomarkers in the association between 25OHD and mortality in NHANES 2015–2018 participants.

25OHD, 25-hydroxyvitamin D; ACME, average causal mediated effect; ADE, average direct effect; ALT, alanine aminotransferase; AST, aspartate aminotransferase; CI, confidence interval; CPK, creatine phosphokinase; HDL-C, high-density lipoprotein-cholesterol; LDH, lactate dehydrogenase; LDL-C, low-density lipoprotein-cholesterol.

4 Discussion

This nationwide cohort study provides a novel and comprehensive exploration of the complex associations between vitamin D status, lipid metabolism, inflammation, cellular injury, and mortality, with a key focus on identifying gender-specific variations across different stages of glucose dysregulation. This nuanced approach allows for the identification of subgroups that may benefit most from targeted interventions in contrast to a one-size-fits-all approach.

The pronounced gender-specific associations we observed between vitamin D, lipid parameters, and inflammatory biomarkers are consistent with previous research, although specific findings across studies vary (18–20). For instance, Li et al. (20) reported that 25OHD levels were inversely correlated with LDL-C and total cholesterol (TC) in men aged over 50, but positively correlated with LDL-C and total cholesterol in women before age 50. In contrast, our study found that serum 25OHD was not associated with LDL-C in the overall population and was positively associated with total cholesterol only in normoglycemic females. These gender differences may be attributed to several factors. Firstly, sex hormones such as estrogen and testosterone play critical roles in lipid metabolism and immune regulation, which may influence how vitamin D exerts its biological effects. An animal study has shown that estrogen enhances vitamin D signaling by modulating hepatic estrogen receptor expression and vitamin D metabolic enzymes (21). In ovariectomized mice, concurrent deficiency of estrogen and vitamin D led to exacerbated hepatic steatosis, increased expression of lipogenesis-related genes [e.g., sterol regulatory element-binding protein 1c (SREBP1c), fatty acid synthase (FASn)], and elevated inflammatory cytokines [e.g., interleukin-6 (IL-6), tumor necrosis factor-α (TNF-α)], suggesting a synergistic role of estrogen and vitamin D in regulating lipid homeostasis and inflammation. Secondly, genetic factors may contribute to the observed disparities, as genes involved in vitamin D metabolism, vitamin D receptor (VDR) expression, and downstream signaling pathways might exhibit sex-specific polymorphisms or expression patterns. For instance, studies in rats have shown that the hepatic expression of Vdr (the gene encoding VDR) and Cyp2r1 (a key 25-hydroxylase involved in vitamin D activation) is significantly higher in females compared with males (22). Such genetically determined differences in VDR levels and the activity of vitamin D metabolizing enzymes can logically lead to distinct intracellular vitamin D signaling capacities and quantitatively different physiological responses to vitamin D between sexes. Thirdly, differences in body composition, particularly body fat distribution, may affect vitamin D bioavailability since vitamin D is fat-soluble and can be sequestered in adipose tissue. Research by Cominacini et al. (23) indicates that visceral adipose tissue (VAT) particularly can trap vitamin D, thereby diminishing its metabolic availability. Their study also observed that distinct body fat distribution patterns in males and females correlated differently with vitamin D levels, suggesting these sex-specific adiposity characteristics contribute to differential vitamin D storage and bioavailability between men and women.

The association of vitamin D status with mortality risk among populations has been evaluated for several years; however, conclusions remain uncertain due to contradictory findings (24). In our study, mortality risk slightly declined with increasing serum 25OHD levels up to a threshold among normoglycemic, prediabetic, and female diabetic participants—aligning with several studies reporting inverse associations between vitamin D levels and mortality to a certain degree (6, 25–28). However, some studies have not found a significant association between blood 25OHD levels and all-cause mortality (29, 30). These discrepancies extend to intervention studies on vitamin D supplementation, where some studies reported beneficial effects on mortality (31, 32), while others did not (33, 34). Variations in study design, population characteristics, vitamin D dosage and form, as well as follow-up duration, may all contribute to these inconsistent findings (24, 34).

Vitamin D appears to exert its protective effects primarily by modulating immune pathways, such as regulating pro-inflammatory cytokines like TNF-α, IL-6, and interleukin-1β (IL-1β) (35). This regulation is often mediated through the VDR, which can interfere with pro-inflammatory signaling cascades, such as the nuclear factor-κB (NF-κB) pathway, a central regulator of cytokine gene expression (36). The potential links between vitamin D and lipid parameters, inflammation, cellular damage, and mortality prompted us to investigate whether these factors mediate the association between vitamin D and mortality. To mitigate the limitations of single biomarkers, we incorporated composite indices including VAI, SII, and SIRI for a more comprehensive assessment. Regarding the significant mediation by SIRI in prediabetic males, this finding should be interpreted with caution. Although the ACME was statistically significant, the effect size was modest, suggesting that SIRI may play only a limited role in the observed association. Therefore, the result should be viewed as a potential direction for future research rather than conclusive evidence of a clinically meaningful pathway. Larger-scale epidemiological research and well-designed experimental studies are warranted to further elucidate the underlying mechanisms and validate these findings across diverse populations.

Contrary to expectations, our mediation analyses revealed that neither lipid-related nor injury-related biomarkers significantly mediated the association between vitamin D and mortality across subgroups. This absence of mediation is particularly noteworthy, suggesting that inflammation may represent a more dominant and statistically detectable pathway than lipid modulation in the vitamin D–mortality relationship. Nevertheless, the lack of significant findings does not rule out potential involvement of these pathways, as previous studies have reported inconsistent results. For example, Huang et al. (37) found no association between plasma 25OHD and lipid metabolism markers, and reported no significant causal relationship. In contrast, Lhilali et al. (38) observed that low 25OHD levels were associated with elevated lipid profiles and atherogenic indices. Similar inconsistencies were found in studies exploring injury-related biomarkers. Bello et al.’s (39) meta-analysis concluded that vitamin D supplementation had no significant effect on LDH or creatine kinase, whereas Reichrath et al. (40) reported that average serum 25OHD levels ≥ 10 ng/mL were associated with significantly lower LDH levels, with a 3.86 U/L decrease per 1 ng/mL increase in 25OHD. Taken together, these findings highlight the heterogeneous and context-dependent nature of vitamin D’s biological effects, underscoring the need for further mechanistic investigations in diverse populations and clinical settings.

This study benefits from a large, nationally representative sample, allowing for robust analysis and generalizability to the broader United States population. By further stratifying our analysis by gender and glycemic status, we were able to identify important group-specific associations and different mediation pathways that would likely be missed in an analysis of the whole population. Taken together, our findings highlight why a ‘one-size-fits-all’ approach to vitamin D may be insufficient. Instead, our results support the need for a more personalized approach, where vitamin D interventions could be tailored to an individual’s specific characteristics, such as their gender and glycemic state, to achieve better health outcomes. Nevertheless, several limitations in our study still need to be considered. First and foremost, while mortality outcomes were ascertained prospectively, baseline vitamin D levels and various biomarkers—including potential mediators—were measured cross-sectionally. Therefore, the temporal order required for causal inference is lacking. As such, the observed associations and estimated mediation effects should be interpreted with caution. Future longitudinal studies assessing changes in both vitamin D status and biomarker profiles over time are warranted to further elucidate and confirm these potentially mediated relationships; additionally, the reliance on single time-point 25OHD measurements—common in many observational studies including NHANES—may not capture long-term vitamin D status. This introduces potential exposure misclassification due to seasonal fluctuations, supplement use, and lifestyle variation, which may not only bias internal estimates but also complicate comparison across studies. These limitations should be considered when interpreting the observed associations; furthermore, despite adjusting for key confounders, the potential for residual confounding remains. Our models could not fully account for several unmeasured variables, such as detailed dietary habits, specific physical activity patterns, or the precise dosage and duration of vitamin D supplementation. These lifestyle factors could correlate with both vitamin D status and cardiometabolic health, thereby potentially confounding the observed associations. Therefore, while our findings are robust to the included covariates, we cannot entirely exclude the possibility that these unmeasured lifestyle factors may have contributed to the observed associations; and finally, as all participants were from the United States, caution is warranted when generalizing the findings to other populations.

In conclusion, our study demonstrates that the associations of vitamin D with lipid metabolism, inflammation, and cellular injury are strongly influenced by gender and glycemic status. Systemic inflammation partially mediates the relationship between vitamin D and mortality in prediabetic males. Although further research is needed to confirm these findings, our study suggests that vitamin D intervention strategies may need to be tailored based on an individual’s gender and glycemic status.

Statements

Data availability statement

Publicly available datasets were analyzed in this study. This data can be found here: https://wwwn.cdc.gov/nchs/nhanes/; www.cdc.gov/nchs/linked-data/mortality-files/index.html.

Ethics statement

The studies involving humans were approved by National Center for Health Statistics (NCHS) Ethics Review Board (ERB). The studies were conducted in accordance with the local legislation and institutional requirements. The participants provided their written informed consent to participate in this study.

Author contributions

MC: Project administration, Writing – review & editing, Funding acquisition, Conceptualization, Methodology. XJ: Visualization, Formal analysis, Software, Methodology, Writing – original draft. XX: Writing – original draft, Data curation, Funding acquisition, Validation, Formal analysis. SZ: Writing – original draft, Validation. YS: Writing – original draft, Data curation. YY: Writing – original draft, Software, Visualization. YG: Funding acquisition, Writing – original draft, Data curation, Validation. YH: Supervision, Conceptualization, Writing – review & editing, Project administration, Funding acquisition.

Funding

The author(s) declare that financial support was received for the research and/or publication of this article. This research was funded by the Suzhou Science and Education Strengthening Health Youth Project, grant numbers KJXW2021011 and KJXW2022014; the Suzhou Medical Innovation and Application Research Project, grant number SKY2023119; the Extracurricular Academic Research Fund for Undergraduates, grant number KY2024266B; the Priority Academic Program Development of Jiangsu Higher Education Institutions.

Acknowledgments

The authors appreciate the efforts of the NHANES staff, researchers, and participants whose contributions have made this valuable dataset accessible for public health research.

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Generative AI statement

The authors declare that no Gen AI was used in the creation of this manuscript.

Any alternative text (alt text) provided alongside figures in this article has been generated by Frontiers with the support of artificial intelligence and reasonable efforts have been made to ensure accuracy, including review by the authors wherever possible. If you identify any issues, please contact us.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

    Glossary

  • 25OHD

    25-hydroxyvitamin D

  • ALT

    alanine aminotransferase

  • AST

    aspartate aminotransferase

  • BMI

    body mass index

  • CPK

    creatine phosphokinase

  • FASn

    fatty acid synthase

  • FPG

    fasting plasma glucose

  • HbA1c

    hemoglobin A1c

  • HDL-C

    high-density lipoprotein cholesterol

  • hs-CRP

    high-sensitivity C-reactive protein

  • IL-1β

    interleukin-1β

  • IL-6

    interleukin-6

  • LDH

    lactate dehydrogenase

  • LDL-C

    low-density lipoprotein cholesterol

  • NF-κB

    nuclear factor-κB

  • NDI

    national death index

  • NHANES

    National Health and Nutrition Examination Survey

  • RDW

    red cell distribution width

  • SII

    systemic immune-inflammation index

  • SIRI

    systemic inflammation response index

  • SREBP1c

    sterol regulatory element-binding protein 1c

  • TC

    total cholesterol

  • TG

    triglycerides

  • TNF-α

    tumor necrosis factor-α

  • VAI

    visceral adiposity index

  • WC

    waist circumference

References

  • 1.

    AoTKikutaJIshiiM. The effects of vitamin D on immune system and inflammatory diseases. Biomolecules. (2021) 11:1624. doi: 10.3390/biom11111624

  • 2.

    SurduAMPînzariuOCiobanuDMNegruAGCainapSSLazeaCet al. vitamin d and its role in the lipid metabolism and the development of atherosclerosis. Biomedicine. (2021) 9:172. doi: 10.3390/biomedicines9020172

  • 3.

    KamimuraDYimerWKShahAMMentzRJOshunbadeAHamidAet al. Vitamin D Levels in black americans and the association with left ventricular remodeling and incident heart failure with preserved ejectin fraction: the Jackson Heart Study. J Card Fail. (2023) 29:150–7. doi: 10.1016/j.cardfail.2022.07.049

  • 4.

    BerridgeMJ. Vitamin D deficiency and diabetes. Biochem J. (2017) 474:1321–32. doi: 10.1042/BCJ20170042

  • 5.

    ZhuKKnuimanMDivitiniMHungJLimEMCookeBRet al. Lower serum 25-hydroxyvitamin D is associated with colorectal and breast cancer, but not overall cancer risk: a 20-year cohort study. Nutr Res. (2019) 67:100–7. doi: 10.1016/j.nutres.2019.03.010

  • 6.

    ParkKYHanKYDHwangHSParkHKParkK. Serum 25-hydroxyvitamin D concentrations are inversely associated with all-cause mortality among Koreans: a nationwide cohort study. Nutr Res. (2023) 113:49–58. doi: 10.1016/j.nutres.2023.02.008

  • 7.

    HolickMF. revisiting vitamin d guidelines: a critical appraisal of the literature. Endocr Pract. (2024) 30:1227–41. doi: 10.1016/j.eprac.2024.10.011

  • 8.

    WierzbickaAOczkowiczM. Sex differences in vitamin D metabolism, serum levels and action. Br J Nutr. (2022) 128:2115–30. doi: 10.1017/S0007114522000149

  • 9.

    AlaidarousTAAlkahtaniNMAljuraibanGSAbulmeatyMMA. Impact of the glycemic control and duration of type 2 diabetes on vitamin D level and cardiovascular disease risk. J Diabetes Res. (2020) 2020:1–10. doi: 10.1155/2020/8431976

  • 10.

    EralpBSeferE. Reference-free inferring of transcriptomic events in cancer cells on single-cell data. BMC Cancer. (2024) 24:607. doi: 10.1186/s12885-024-12331-5

  • 11.

    BijariMJangjooSEmamiNRajiSMottaghiMMoallemRet al. The accuracy of visceral adiposity index for the screening of metabolic syndrome: a systematic review and meta-analysis. Int J Endocrinol. (2021) 2021:1–14. doi: 10.1155/2021/6684627

  • 12.

    LiJHeDYuJZChenSHWuQChengZXet al. dynamic status of sii and siri alters the risk of cardiovascular diseases: evidence from Kailuan Cohort Study. J Inflamm Res. (2022) 15:5945–57. doi: 10.2147/JIR.S378309

  • 13.

    SeligmanHKLeviRAdebiyiVOColeman-JensenAGuthrieJFFrongilloEA. Assessing and monitoring nutrition security to promote healthy dietary intake and outcomes in the United States. Annu Rev Nutr. (2023) 43:409–29. doi: 10.1146/annurev-nutr-062222-023359

  • 14.

    ElsayedNACommPPMccoyRGAleppoGBalapattabiKBeverlyEAet al. 2. Diagnosis and classification of diabetes: standards of care in diabetes-2025. Diabetes Care. (2025) 48:S27–49. doi: 10.2337/dc25-S002

  • 15.

    CandemirMKiziltunçENurkoçSSahinarslanA. Relationship between systemic immune-inflammation index (SII) and the severity of stable coronary artery disease. Angiology. (2021) 72:575–81. doi: 10.1177/0003319720987743

  • 16.

    LiXWCuiLXuHY. Association between systemic inflammation response index and chronic kidney disease: a population-based study. Front Endocrinol. (2024) 15:1329256. doi: 10.3389/fendo.2024.1329256

  • 17.

    Parra-SotoSBoonporJLynskeyNArayaCHoFPellJPet al. Association between visceral adiposity index and cancer risk in the UK Biobank cohort. Cancer. (2025) 131:e35576. doi: 10.1002/cncr.35576

  • 18.

    AlmesriNDasNSAliMEGumaaKGihaHA. Gender-dependent association of vitamin D deficiency with obesity and hypercholesterolemia (LDLC) in adults. Endoc Metab Immune Disord. (2020) 20:425–36. doi: 10.2174/1871530319666191009154528

  • 19.

    AbudawoodMTabassumHAnsarSAlmosaKSobkiSAliMNet al. Assessment of gender-related differences in vitamin D levels and cardiovascular risk factors in Saudi patients with type 2 diabetes mellitus. Saudi J Biol Sci. (2018) 25:31–6. doi: 10.1016/j.sjbs.2017.04.001

  • 20.

    LiXTLiuYNWangJYChenXReichetzederCElitokSet al. Vitamin D Is associated with lipid metabolism: a sex- and age-dependent analysis of a large outpatient cohort. Nutrients. (2024) 16:3936. doi: 10.3390/nu16223936

  • 21.

    BorgesCCBringhentiIMandarim-de-LacerdaCAAguilaMB. Vitamin D deficiency aggravates the liver metabolism and inflammation in ovariectomized mice. Biomed Pharmacother. (2018) 107:878–88. doi: 10.1016/j.biopha.2018.08.075

  • 22.

    OczkowiczMSzymczykBSwiatkiewiczMFurgal-DzierzukIKoseniukAWierzbickaAet al. Analysis of the effect of vitamin D supplementation and sex on and gene expression in Wistar rats’ tissues. J Steroid Biochem Mol Biol. (2021) 212:105918. doi: 10.1016/j.jsbmb.2021.105918

  • 23.

    CominaciniMFumaneriABalleriniLBraggioMValentiMTCarbonareLD. Unraveling the connection: visceral adipose tissue and vitamin D levels in obesity. Nutrients. (2023) 15:4259. doi: 10.3390/nu15194259

  • 24.

    HeathAKKimIYHodgeAMEnglishDRMullerDC. Vitamin D status and mortality: a systematic review of observational studies. Int J Environ Res Public Health. (2019) 16:383. doi: 10.3390/ijerph16030383

  • 25.

    LuoWXuDZhangJZhouYYangQLvQet al. Low serum 25-hydroxyvitamin D levels are associated with increased cardiovascular morbidity and mortality. Postgrad Med. (2023) 135:93–101. doi: 10.1080/00325481.2022.2161250

  • 26.

    WanZZGuoJYPanAChenCLiuLGLiuG. Association of serum 25-hydroxyvitamin D concentrations with all-cause and cause-specific mortality among individuals with diabetes. Diabetes Care. (2021) 44:350–7. doi: 10.2337/dc20-1485

  • 27.

    JinXRXiongSZJuSYZengYYanLJLYaoY. Serum 25-hydroxyvitamin D, albumin, and mortality among Chinese older adults: a population-based longitudinal study. J Clin Endocrinol Metab. (2020) 105:2762–70. doi: 10.1210/clinem/dgaa349

  • 28.

    YeHWLiYXLiuSMZhangXFLiangHZWangYet al. Association between serum 25-hydroxyvitamin D and vitamin D dietary supplementation and risk of all-cause and cardiovascular mortality among adults with hypertension. Nutr J. (2024) 23:33. doi: 10.1186/s12937-024-00914-8

  • 29.

    LinSWChenWFanJHDawseySMTaylorPRQiaoYLet al. Prospective study of serum 25-hydroxyvitamin D concentration and mortality in a Chinese population. Am J Epidemiol. (2012) 176:1043–50. doi: 10.1093/aje/kws285

  • 30.

    JevalikarGMithalASinghASharmaRFarooquiKJMahendruSet al. Lack of association of baseline 25-hydroxyvitamin D levels with disease severity and mortality in Indian patients hospitalized for COVID-19. Sci Rep. (2021) 11:6258. doi: 10.1038/s41598-021-85809-y

  • 31.

    ShaSNguyenTMNKuzniaSNiedermaierTZhuABrennerHet al. Real-world evidence for the effectiveness of vitamin D supplementation in reduction of total and cause-specific mortality. J Intern Med. (2022) 293:384–97. doi: 10.1111/joim.13578

  • 32.

    AcharyaPDaliaTRankaSSethiPOniOASafarovaMSet al. The effects of vitamin D supplementation and 25-hydroxyvitamin D levels on the risk of myocardial infarction and mortality. J Endocr Soc. (2021) 5:bvab124. doi: 10.1210/jendso/bvab124

  • 33.

    MansonJECookNRLeeIMChristenWBassukSSMoraSet al. Vitamin D supplements and prevention of cancer and cardiovascular disease. N Engl J Med. (2019) 380:33–44. doi: 10.1056/NEJMoa1809944

  • 34.

    ZhangYFangFTangJJJiaLFengYNXuPet al. Association between vitamin D supplementation and mortality: systematic review and meta-analysis. BMJ. (2019) 366:l4673. doi: 10.1136/bmj.l4673

  • 35.

    XiaohuaGDongdongLXiaotingNShuopingCFeixiaSHuajunYet al. Severe Vitamin D deficiency is associated with increased expression of inflammatory cytokines in painful diabetic peripheral neuropathy. Front Nutr. (2021) 8:612068. doi: 10.3389/fnut.2021.612068

  • 36.

    ArganoCMirarchiLAmodeoSOrlandoVTorresACorraoS. The role of vitamin D and its molecular bases in insulin resistance, diabetes, metabolic syndrome, and cardiovascular disease: state of the art. Int J Mol Sci. (2023) 24:15485. doi: 10.3390/ijms242015485

  • 37.

    HuangTAfzalSYuCQGuoYBianZYangLet al. Vitamin D and cause-specific vascular disease and mortality: a Mendelian randomisation study involving 99,012 Chinese and 106,911 European adults. BMC Med. (2019) 17:160. doi: 10.1186/s12916-019-1401-y

  • 38.

    LhilaliIZouineNGodderisLEl MidaouiAEl JaafariSFilali-ZegzoutiY. Relationship between Vitamin D insufficiency, lipid profile and atherogenic indices in healthy women aged 18–50 years. Eur J Invest Health Psychol Educ. (2024) 14:2337–57. doi: 10.3390/ejihpe14080155

  • 39.

    BelloHJCaballero-GarcíaAPérez-ValdecantosDRocheENoriegaDCCórdova-MartínezA. Effects of vitamin D in post-exercise muscle recovery. A systematic review and meta-analysis. Nutrients. (2021) 13:4013. doi: 10.3390/nu13114013

  • 40.

    ReichrathJBiersackFWagenpfeilSSchoepeJPfoehlerCSaternusRet al. Low vitamin D status predicts poor clinical outcome in advanced melanoma treated with immune checkpoint or BRAF/MEK inhibitors: a prospective non-interventional side-by-side analysis. Front Oncol. (2022) 12:839816. doi: 10.3389/fonc.2022.839816

Summary

Keywords

vitamin D, dyslipidemia, systemic inflammation, oxidative stress, cardiometabolic risk, gender differences, NHANES, glycemic control

Citation

Cai M, Jiang X, Xu X, Zhao S, Sun Y, Yang Y, Gu Y and Huang Y (2025) Associations of vitamin D with lipid metabolism, inflammation, and mortality vary by glycemic status and gender: a nationwide prospective study. Front. Nutr. 12:1597527. doi: 10.3389/fnut.2025.1597527

Received

21 March 2025

Accepted

17 September 2025

Published

08 October 2025

Volume

12 - 2025

Edited by

Kathryn Hart, University of Surrey, United Kingdom

Reviewed by

Alessandro De Oliveira, Universidade Federal de São João del-Rei, Brazil

Emre Sefer, Özyeğin University, Türkiye

Yuan Tian, Affiliated Hospital of Shandong University of Traditional Chinese Medicine, China

Muhammad Mazhar, Shenzhen University, China

Updates

Copyright

*Correspondence: Meizhi Cai, Yifan Huang,

†These authors have contributed equally to this work

Disclaimer

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.

Outline

Figures

Cite article

Copy to clipboard


Export citation file


Share article

Article metrics