Your new experience awaits. Try the new design now and help us make it even better

ORIGINAL RESEARCH article

Front. Nutr., 04 November 2025

Sec. Nutritional Epidemiology

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

This article is part of the Research TopicClinical Research on Nutrition and Children's Health: Focusing on Genetics, Metabolism, and the Gut MicrobiomeView all 6 articles

Association of vitamin D levels with metabolic dysfunction-associated fatty liver disease in children aged 12–18 years


Xuejie Gao,&#x;Xuejie Gao1,2Yuyun Chen&#x;Yuyun Chen1Xinrui Wang,Xinrui Wang1,2Yuehang ChenYuehang Chen1Xiaoyan ChenXiaoyan Chen1Haibo Li
Haibo Li2*Hong Ye,
Hong Ye1,2*
  • 1Department of Digestive & Nutrition, Fujian Children's Hospital (Fujian Branch of Shanghai Children's Medical Center), College of Clinical Medicine for Obstetrics & Gynecology and Pediatrics, Fujian Medical University, Fuzhou, Fujian, China
  • 2Fujian Maternity and Child Health Hospital, College of Clinical Medicine for Obstetrics & Gynecology and Pediatrics, Fujian Medical University, Fuzhou, Fujian, China

Objective: This study examines the association between serum vitamin D levels and the prevalence of metabolic dysfunction-associated fatty liver disease (MAFLD) in adolescents, along with potential modifying factors.

Methods: Data from 950 adolescents aged 12–18 years in the National Health and Nutrition Examination Survey (NHANES) 2017–2018 were analyzed. MAFLD was defined using hepatic steatosis and metabolic dysfunction criteria. Serum 25(OH)D levels were measured, and weighted logistic regression and restricted cubic spline models were applied to assess their association with MAFLD risk. Stratified analyses were also conducted.

Results: Lower serum 25(OH)D levels were significantly associated with higher MAFLD risk (p < 0.001), showing a nonlinear dose-response relationship. Adolescents with 25(OH)D ≥ 75 nmol/L had a 57% lower risk of MAFLD compared to those with levels < 50 nmol/L. Stratified analysis indicated that the protective effect of vitamin D was more evident in individuals with higher retinol levels, though retinol alone was not significantly associated with MAFLD.

Conclusion: Vitamin D deficiency is significantly associated with MAFLD in adolescents, with a nonlinear dose-response relationship modulated by retinol status. These findings underscore the potential role of vitamin D in MAFLD prevention and provide a basis for further prospective or intervention studies.

Introduction

Metabolic dysfunction-associated fatty liver disease (MAFLD), formerly named non-alcoholic fatty liver disease (NAFLD), is a significant chronic liver disease worldwide. A recent meta-analysis revealed that the global prevalence of MAFLD is 30%, reflecting a 50% increase from 1990 to 2019 (1). This increase presents a notable public health burden (2). Notably, with the socioeconomic transformation and lifestyle changes, the prevalence of MAFLD in adolescents has shown a significant trend of younger age (3, 4). The primary pathological feature of MAFLD is fat deposition in hepatocytes (≥5%), driven by metabolic disorders (5). This pathological process is closely linked to abnormal liver function and significantly raises the risk of developing type 2 diabetes and cardiovascular diseases in adulthood through mechanisms such as insulin resistance and chronic inflammation (6). Therefore, MAFLD poses a long-term threat to public health.

As a fat-soluble molecule with steroid hormone characteristics (7), the multi-dimensional metabolic regulatory functions of vitamin D are gradually being uncovered, revealing its extensive physiological roles. Previous studies have established the central role of vitamin D in bone mineral metabolism, maintaining bone health by promoting intestinal calcium absorption and regulating bone remodeling (8, 9). Recent research has shown that vitamin D is broadly involved in physiological processes such as immune regulation, hormone secretion, and cell proliferation and differentiation (10). Vitamin D receptors are widely distributed in the liver, adipose tissue, and immune cells, where they play a crucial role in maintaining metabolic homeostasis by regulating glucose and lipid metabolism, suppressing the release of proinflammatory factors, and improving insulin sensitivity (11). In the liver, vitamin D is metabolized by 25-hydroxylase (CYP2R1) into 25-hydroxyvitamin D [25(OH)D] (12). Vitamin D status is typically classified according to serum 25(OH)D concentrations. Generally, vitamin D deficiency is defined as serum 25(OH)D levels below 20 ng/mL (50 nmol/L), while insufficiency refers to serum concentrations between 20 and 30 ng/mL (50–75 nmol/L). By contrast, vitamin D sufficiency is generally considered to be achieved when serum 25(OH)D levels exceed 30 ng/mL (75 nmol/L) (13).

Several studies have suggested that vitamin D deficiency may play a critical role in the development and progression of various metabolic diseases, including insulin resistance, obesity, hypertension, metabolic syndrome, and cardiovascular disease (1417). In addition, the association between vitamin D levels and NAFLD has attracted increasing attention (18). A cross-sectional study revealed that low serum 25(OH)D levels in the general population are associated with an increased risk of NAFLD (19). Furthermore, a study from Korea indicated that vitamin D deficiency may be considered an independent risk factor for suspected NAFLD in adolescents (20). Besides, obesity constitutes one of the primary risk factors for MAFLD. It significantly increases the risk of developing MAFLD by promoting hepatic fat accumulation and contributing to metabolic dysfunction (21). Concurrently, obesity is also associated with vitamin D deficiency (22), suggesting that vitamin D may mediate the effect of obesity on MAFLD pathogenesis. However, research on the association between MAFLD and vitamin D in adolescents aged 12–18 years remains relatively limited.

Based on the current research advancements in the potential role of vitamin D in metabolic diseases, this study proposes the scientific hypothesis that adolescents with lower 25(OH)D levels may be more susceptible to MAFLD. Given the new nomenclature and diagnostic criteria for MAFLD, the population with MAFLD differs from that with NAFLD in many respects (23). Therefore, it is necessary to further clarify the characteristics and influencing factors of MAFLD. This study utilizes data from the 2017–2018 National Health and Nutrition Examination Survey (NHANES) to investigate the association between serum 25(OH)D status and MAFLD, and to further explore potential effect modifiers.

Material and methods

Study population

The National Health and Nutrition Examination Survey (NHANES) represents a continuous research initiative designed to evaluate health and dietary conditions across a diverse, non-institutionalized U.S. population through its sophisticated, multistage probability sampling framework (24). For this investigation, we utilized data from the 2017–2018 NHANES cycle, specifically chosen for its inclusion of liver ultrasound transient electrography assessments. Our analysis focused on adolescents aged 12–18 years. Participants lacking essential data points were excluded, including those without liver ultrasound transient electrography results and missing critical metabolic dysfunction indicators such as body mass index (BMI), waist circumference (WC), high-density lipoprotein (HDL) -cholesterol levels, blood pressure (BP) readings (both systolic and diastolic), hemoglobin A1C (HbA1C), and fasting plasma glucose (FPG) measurements (Figure 1). The National Center for Health Statistics (NCHS) oversees the administration and maintenance of NHANES. The NCHS Institutional Review Board/Ethic Review Board granted approval for the NHANES protocols (Continuation of Protocol #2011-17, effective until 26 October 2017; Protocol #2018-01, effective from 26 October 2017). All participants provided written informed consent.

Figure 1
Flowchart illustrating participant selection from NHANES 2017-2018. Begins with 9,254 participants, excludes 8,203 aged over eighteen or under twelve, leaving 1,051 aged twelve to eighteen. Excludes 95 missing liver ultrasound data, resulting in 956. Further excludes 6 with missing biochemical data, yielding a final sample of 950.

Figure 1. Flowchart of participant inclusion and exclusion.

Data collection

Study subjects provided self-reported information regarding their age, gender, racial-ethnic background (classified as non-Hispanic White, non-Hispanic Black people, Hispanic, or other racial groups), family income (<$100,000, $100,000 and Over), season of blood collection (November 1 through April 30, or May 1 through October 31), secondhand smoke exposure, days physically active (≥5 Days/Week, at least 60 min/day), and prior medical conditions. Anthropometric measurements such as stature (in centimeters), body mass (in kilograms), and abdominal circumference (in centimeters) were obtained during the Mobile Examination Center (MEC) assessment; BMI was derived by dividing the weight in kilograms by the square of height in meters.

Biological sample analysis

Venous blood samples were collected from participants at the NHANES MEC facility. The concentrations of 25(OH)D2 and 25(OH)D3 were determined through high-performance liquid chromatography coupled with tandem mass spectrometry (HPLC-MS/MS). Quantitative estimation was performed by comparing the response ratio (analytic peak area/internal standard peak area) of unknown samples with that of known analytic concentrations in calibration solutions (25). The total serum vitamin D [25(OH)D] level was calculated as the combined concentrations of 25(OH)D2 and 25(OH)D3. FPG, lipid profiles (HDL-C, triglycerides [TG], total cholesterol [TC], low-density lipoprotein Cholesterol [LDL-C]), liver and renal function indicators (alanine aminotransferase [ALT], gamma-glutamyl transferase [GGT], creatinine) were measured using standardized laboratory methods. Retinol levels in serum were analyzed using an adapted high-performance liquid chromatography method with photodiode array detection, complemented by spectrophotometric techniques for quantitative assessment.

Defining MAFLD

Given the lack of explicit international consensus, hepatic steatosis was identified using a median controlled attenuation parameter (CAP) score of ≥248 dB/m, as this threshold demonstrated reliable diagnostic accuracy in a recent meta-analysis of individual patient data (26). Following the latest expert consensus (27), MAFLD in adolescents was defined by the coexistence of hepatic steatosis (assessed via CAP) along with one or more of the following conditions: overweight or obesity [BMI > 1 SD above the World Health Organization growth reference median (28) or as a waist circumference > 90th percentile for age and sex (29)], prediabetes and diabetes [(1) a previous diagnosis of diabetes; (2) a HbA1c level ≥ 5.7% (48 mmol/mol); or (3) a FPG ≥ 100 mg/dL (30)], and the presence of two or more metabolic disturbances. These metabolic disturbances encompassed increase BP (a systolic BP > 130 mmHg and/or a diastolic BP > 85 mmHg), triglyceride concentrations ≥ 150 mg/dL, HDL-cholesterol levels below 40 mg/dL, and a triglycerides-to-HDL cholesterol ratio exceeding 2.25 (with adult MAFLD criteria extended to adolescents aged 16 years and above).

Statistical analysis

Given a complex, multi-stage sample design used in the NHANES, we applied appropriate sample weights in all analyses to account for clustering, stratification, non-response, and oversampling population. Continuous variables were expressed presented as weighted mean (standard error), categorical variables were presented as unweighted frequencies (%). Chi-squared test with Rao & Scott's second-order correction; Wilcoxon rank-sum test for complex survey samples.

Weighted restricted cubic spline regression models with three knots were used to further examine the nonlinear relationships of 25(OH)D and retinol with MAFLD. Non-linearity tested by including a quadratic term in the regression models. Survey-weight adjusted multivariable logistic regressions were performed to determine the independent associations between the 25(OH)D or retinol (both continuous variables and dichotomous variable) and MAFLD. Confounders were selected based on clinical significance and prior epidemiological evidence. Pre-specified subgroup analyses were conducted. Interaction across subgroups was tested using the likelihood ratio test. Moreover, interaction spline was produced by using the estimated odds ratio (OR) for MAFLD of a logistic model with a 25(OH)D—retinol level interaction term. Missing covariates were addressed using imputation by chained equations (31). All covariates had missing rates below 5%, except for physical activity (missing 14.5%). Sensitivity analyses were conducted through complete-case exclusion. Furthermore, an additional sensitivity analysis was performed to assess the robustness of findings to unmeasured confounding using the E-value methodology of VanderWeele and Ding (32).

A 2-tailed p ≤ 0.05 was considered statistically significant, and the analyses were performed using the “survey”, “rms”, and “interaction RCS” package of R software (Version 4.2.2, http://www.R-project.org, The R Foundation) and Free Statistics analysis platform (Version 2.1.1, Beijing, China, http://www.clinicalscientists.cn/freestatistics).

Results

Characteristics of study participants

Table 1 outlines the characteristics of the study participants (N = 950), categorized by 25(OH)D status. The overall mean age was 15.03 years, with a slight variation across 25(OH)D categories. The majority of participants were non-Hispanic Black people (50.60%), followed by non-Hispanic White (12.92%). Ethnic distribution significantly differed by 25(OH)D status, with Hispanic-Mexican Americans and non-Hispanic Black peoples more prevalent in the lower 25(OH)D categories. Overweight and abdominal obesity were more common in participants with lower 25(OH)D levels, while metabolic abnormalities and steatosis showed a trend toward higher prevalence in these groups. Notably, the prevalence of MAFLD was significantly higher in individuals with lower 25(OH)D levels, highlighting a potential association between vitamin D deficiency and liver health.

Table 1
www.frontiersin.org

Table 1. Characteristics of study participants, divided by 25(OH)D status.

Vitamin D levels and MAFLD

The nonlinear curve fitting analysis established a significant association between 25(OH)D levels and MAFLD risk (Figure 2). The overall association was highly significant (p < 0.001), with a notable p-value for non-linearity at 0.002, suggesting a complex relationship rather than a simple linear trend. To further investigate this association, we conducted logistic regression analysis examining 25(OH)D levels and MAFLD prevalence in U.S. adolescents (Table 2). The analysis revealed a significant inverse dose-response relationship: each 5 nmol/L increase in 25(OH)D was associated with a 10% reduction in MAFLD risk after adjusting for potential confounders (adjusted OR = 0.90, 95% CI: 0.85~0.95, p < 0.001). When examining specific vitamin D thresholds, adolescents with 25(OH)D levels <50 nmol/L (reference group) exhibited the highest MAFLD risk, while those with optimal levels (≥ 75 nmol/L) showed a remarkable 57% risk reduction (OR = 0.43, 95% CI: 0.22~0.82, p = 0.011). The dose-response nature of this protective association was further confirmed by significant trend analysis results (p for trend <0.05).

Figure 2
A graph showing the relationship between 25(OH)D levels in nmol/L (x-axis) and odds ratio with 95% confidence interval (y-axis). The curve indicates a non-linear decline, crossing the reference line at an odds ratio of 1 around a value of 54.7 nmol/L. Shaded areas represent confidence intervals. Vertical dashed lines highlight values at 50 and 75 nmol/L. Statistical values are provided: non-linearity p-value is 0.002, the overall p-value is less than 0.001.

Figure 2. The dose-response associations between 25(OH)D levels and MAFLD. Analyses were adjusted for age, sex, ethnicity, education, season of blood collection, income, secondhand smoke, physically active, low-density lipoprotein cholesterol, alanine aminotransferase, gamma-glutamyl transferase, and creatinine.

Table 2
www.frontiersin.org

Table 2. The associations of 25(OH)D level with MAFLD in US adolescents.

Subgroup analyses and interactions

A significant retinol-dependent association between 25(OH)D levels and MAFLD was observed in interacting fitted curves, the 25(OH)D had a significantly stronger association with MAFLD at higher retinol levels (p for interaction <0.05, Figure 3). Specifically, in stratified analysis, children have a higher retinol nutritional status exhibited an inverse relationship between high 25(OH)D levels and MAFLD (OR = 0.22), contrasting with no significant effect for low retinol levels groups (OR = 0.65) (Table 3). Besides, age group, sex, ethnicity, overweight, abdominal obesity, season of blood collection, and physically active did not influence the relationship (p for interaction > 0.05).

Figure 3
Histogram displaying the relationship between retinol levels in micrograms per deciliter and the odds ratio (OR) for MAFLD per one standard deviation increase in 25(OH)D. Blue and red bars represent retinol distribution; a yellow line with a shaded area indicates the OR trend and confidence interval. Retinol levels range from 20 to 80, with OR values from 0.0 to 3.0.

Figure 3. Odds ratio (OR) of vitamin D for MAFLD as a function of retinol levels. OR for MAFLD per 1-SD increase in the 25(OH)D (P for interaction <0.05). Estimates derived from adjusted logistic models with restricted cubic splines using knots at tertiles (33rd, 67th).

Table 3
www.frontiersin.org

Table 3. Associations of 25(OH)D with MAFLD, divided by age, sex and retinol status.

Retinol levels and MAFLD

In light of the combined health effects of vitamin A and vitamin D nutritional status in children, this study further analyzed the relationship between retinol levels and MAFLD. However, the results did not reach statistical significance. In the restricted cubic spline analysis, the overall p = 0.608 (Supplementary Figure 1). In the logistic regression analysis (Supplementary Table 1), the adjusted OR for retinol per 5 μg/dL increase was 1.03 (95% CI: 0.93-1.14, p = 0.562). Quartile analysis revealed similar trends, with Q1 as the reference group; Q2 (OR = 0.90, 95% CI: 0.40–2.03, p = 0.746), Q3 (OR = 1.25, 95% CI: 0.52–3.03, p = 0.517), and Q4 (OR = 1.09, 95% CI: 0.49–2.43, p = 0.790) showed no significant differences. The trend analysis supported these findings (p for trend > 0.05).

Discussion

This study, based on NHANES 2017–2018 data, comprehensively evaluated the association between serum 25(OH)D levels and MAFLD risk in adolescents. The results revealed an inverse association between serum 25(OH)D levels and MAFLD risk, demonstrating a non-linear dose-response relationship that was modified by factors such as retinol levels. The results of this study will provide new scientific evidence for the potential role of vitamin D in the onset and progression of MAFLD.

Our study found significant differences in serum 25(OH)D levels across ethnic and socioeconomic groups. Higher concentrations were most common in Non-Hispanic Black people individuals, those with lower family incomes, and people who were not overweight or abdominally obese. Notably, we observed a significantly higher prevalence of MAFLD in populations with lower serum 25(OH)D levels, with MAFLD prevalence progressively decreasing as 25(OH)D levels increased. Similar findings have been reported in previous studies; a prospective study demonstrated that lower 25(OH)D concentrations were significantly associated with NAFLD in adolescents, and this association was independent of obesity (33). Additionally, an observational study conducted in Brazil found a high prevalence of vitamin D deficiency and insufficiency among obese adolescents, which was closely associated with metabolic changes (34).Vitamin D may exert a protective effect against hepatic fat deposition by improving insulin sensitivity (35), inhibiting inflammatory responses, and regulating lipid metabolism (36). Therefore, maintaining adequate vitamin D status may be an important metabolic regulatory factor for reducing the risk of MAFLD in adolescents. Additionally, race and family income have been shown to be closely related to the risk of MAFLD (37, 38). Thus, these factors were adjusted for as covariates in our study. The results demonstrated a non-linear dose-response relationship between serum 25(OH)D levels and MAFLD risk in adolescents, which is consistent with previous research. A prospective cohort study in China also found a non-linear negative association between serum 25(OH)D concentrations and all-cause mortality in patients with MAFLD or NAFLD (39). Our findings further demonstrate that vitamin D deficiency exerts a more pronounced adverse effect on MAFLD risk, whereas its protective efficacy appears to stabilize beyond a specific 25(OH)D concentration threshold. This suggests a potential “threshold effect” of vitamin D in regulating metabolic homeostasis. Adequate levels are essential for protection, but higher concentrations may not provide extra benefits. These findings highlight the need for targeted vitamin D interventions. It is more effective to focus on deficient groups rather than supplementing those who already have adequate levels.

This study explored the moderating effects of age levels on the relationship between vitamin D and MAFLD. Age-stratified analysis showed that the protective effect of vitamin D was more significant in adolescents aged ≥ 15 years, while no statistically significant association was observed in the group <15 years old. This age-dependent differential effect may be attributed to pubertal development changes and hormonal fluctuations (40). The adolescent period represents a critical window of metabolic transformation, characterized by substantial increases in growth hormone and sex hormone levels, coupled with progressive insulin resistance (41). These physiological changes potentially enhance susceptibility to MAFLD, thereby increasing hepatic fat accumulation risk (42). Indeed, puberty has been identified as a high-risk period for the development and progression of obesity and metabolic dysregulation (43). Vitamin D acts as a metabolic modulator and may enhance protection by regulating insulin sensitivity, inflammation, and lipid metabolism (44, 45).

Interestingly, we further found differential associations between vitamin D and MAFLD across retinol concentration subgroups. Specifically, among individuals with higher retinol levels, the protective effect of vitamin D against MAFLD was significantly enhanced, while no similar effect was observed in those with lower retinol levels. This observation implies that retinol may enhance vitamin D metabolic regulatory capacity through synergistic interactions, consequently mitigating MAFLD risk. As the main active form of vitamin A, retinol is an essential micronutrient for mammals, primarily stored in the liver, kidneys, and pulmonary tissues (46). Its physiological significance extends to the regulation of cellular differentiation, immune modulation, and lipid homeostasis (47). Studies have confirmed that retinol plays an important role in regulating liver lipid metabolism and inflammatory response (48, 49). A Dutch study (49) demonstrated that retinol ameliorates lipid dysmetabolism and inflammation through modulation of hepatic signaling pathways, providing mechanistic insights into its potential metabolic protective properties. Although retinol can enhance the protective effect of vitamin D, this study did not find a significant independent association between retinol and MAFLD risk. This suggests that retinol may not directly influence MAFLD pathogenesis but rather exerts its metabolic protective effects indirectly through vitamin D pathway modulation. This interpretation aligns with Zeng et al.'s finding (50), which indicates that retinol enhances vitamin D's effects mainly by boosting vitamin D receptor activity, rather than through direct involvement in MAFLD pathways. In summary, retinol may play a synergistic role in vitamin D-mediated metabolic homeostasis regulation by enhancing the activity of the vitamin D signaling pathway, but its preventive effect on MAFLD alone may be limited. Therefore, future nutritional intervention strategies should pay more attention to the interaction between fat-soluble vitamins, especially the synergistic supplementation of vitamin A and vitamin D, in order to more effectively reduce the risk of MAFLD in adolescents.

Based on the nationally representative data from NHANES, this study systematically evaluated the relationship between vitamin D levels and the risk of MAFLD in adolescents and further explored the potential interactions of age and retinol. The results highlighted the potential protective role of vitamin D in the risk of MAFLD in individuals with high retinol levels, where the metabolic protective effects were more significant. However, this study still has some limitations. First, the cross-sectional nature of this study fundamentally limits causal inference, as reverse causality cannot be ruled out. Consequently, caution is warranted when interpreting these findings. Therefore, future prospective cohort studies or intervention studies should be conducted to further clarify the causal role of vitamin D in the occurrence and development of MAFLD and the related biological mechanisms. Second, this study did not investigate factors such as dietary vitamin D intake, sunlight exposure, socioeconomic indicators, and lifestyle factors. Although multiple confounding variables were adjusted for in the statistical models, certain potential confounders remained incompletely controlled. However, we performed sensitivity analyses using the E-value methodology to assess the robustness of our findings to unmeasured confounding. These analyses indicated that our results are less likely to be reversed by unadjusted confounding factors. Third, the ethnic distribution of the sample limits the generalizability of our findings. Specifically, non-Hispanic Black peoples represented 50.6% of participants and Hispanic Mexican Americans 23.1%, while Asian and other ethnic groups were underrepresented. As vitamin D metabolism varies with skin pigmentation and cultural practices, these results should be generalized to other populations with caution. Finally, it should be noted that the observed associations of 25(OH)D with MAFLD, divided by retinol status, may be subject to false positives attributable to multiple testing. In light of this, we strongly recommend that future research prioritize the use of prospective validation cohorts in conjunction with mechanistic studies. This approach will enable a more robust exploration of the association between 25(OH)D and MAFLD.

In conclusion, our investigation revealed a significant inverse correlation between 25(OH)D concentrations and MAFLD risk in adolescent populations, characterized by a nonlinear dose-response relationship. Notably, a trend of enhanced protective effects of vitamin D was observed in the subgroup with higher retinol levels. This exploratory finding suggests that retinol may potentiate the biological activity of vitamin D through synergistic mechanisms, providing a preliminary hypothesis for future research on fat-soluble vitamin interactions. However, these observations require further validation through prospective studies combined with mechanistic experiments. Our results indicate that maintaining adequate vitamin D levels may have positive implications for MAFLD prevention in adolescents. While focusing on specific nutrient supplementation, future research could further explore the potential value of synergistic effects among multiple micronutrients. Such comprehensive nutritional intervention approaches may offer more comprehensive strategic references for the prevention and management of metabolic diseases in adolescents.

Data availability statement

The original contributions presented in the study are included in the article/Supplementary material, further inquiries can be directed to the corresponding authors.

Ethics statement

The studies involving humans were approved by The National Center for Health Statistics Institutional Review Board/Ethic Review Board. The studies were conducted in accordance with the local legislation and institutional requirements. Written informed consent for participation in this study was provided by the participants' legal guardians/next of kin.

Author contributions

XG: Conceptualization, Data curation, Validation, Methodology, Writing – original draft, Software. YuyC: Conceptualization, Validation, Methodology, Writing – original draft, Software. XW: Software, Writing – review & editing, Validation, Data curation. YueC: Validation, Supervision, Software, Writing – review & editing, Investigation. XC: Writing – original draft, Supervision, Visualization, Data curation, Investigation. HL: Supervision, Methodology, Conceptualization, Writing – review & editing, Project administration. HY: Conceptualization, Supervision, Writing – review & editing, Methodology, Funding acquisition, Project administration.

Funding

The author(s) declare that financial support was received for the research and/or publication of this article. This work was supported by the province-level special subsidy funds for the finance department in Fujian Province (No: Fujian Finance Index (2023) 834). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

Acknowledgments

We are grateful to thank all the participants for their valuable contributions.

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 author(s) 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.

Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fnut.2025.1615851/full#supplementary-material

References

1. Younossi ZM, Golabi P, Paik JM, Henry A, Van Dongen C, Henry L. The global epidemiology of nonalcoholic fatty liver disease (Nafld) and nonalcoholic steatohepatitis (Nash): a systematic review. Hepatology. (2023) 77:1335–47. doi: 10.1097/HEP.0000000000000004

PubMed Abstract | Crossref Full Text | Google Scholar

2. Eslam M, Newsome PN, Sarin SK, Anstee QM, Targher G, Romero-Gomez M, et al. A new definition for metabolic dysfunction-associated fatty liver disease: an international expert consensus statement. J Hepatol. (2020) 73:202–9. doi: 10.1016/j.jhep.2020.07.045

PubMed Abstract | Crossref Full Text | Google Scholar

3. Engin A. The definition and prevalence of obesity and metabolic syndrome. Adv Exp Med Biol. (2017) 960:1–17. doi: 10.1007/978-3-319-48382-5_1

PubMed Abstract | Crossref Full Text | Google Scholar

4. Ge X, Zheng L, Wang M, Du Y, Jiang J. Prevalence trends in non-alcoholic fatty liver disease at the global, regional and national levels, 1990-2017: a population-based observational study. BMJ Open. (2020) 10:e036663. doi: 10.1136/bmjopen-2019-036663

PubMed Abstract | Crossref Full Text | Google Scholar

5. Wang MF, Wan B, Wu YL, Huang JF, Zhu YY Li YB. Clinic-pathological features of metabolic associated fatty liver disease with hepatitis B virus infection. World J Gastroenterol. (2021) 27:336–44. doi: 10.3748/wjg.v27.i4.336

PubMed Abstract | Crossref Full Text | Google Scholar

6. Yang K, Song M. New Insights into the pathogenesis of metabolic-associated fatty liver disease (Mafld): gut-liver-heart crosstalk. Nutrients. (2023) 15:3970. doi: 10.3390/nu15183970

PubMed Abstract | Crossref Full Text | Google Scholar

7. Eliades M, Spyrou E. Vitamin D: a new player in non-alcoholic fatty liver disease? World J Gastroenterol. (2015) 21:1718–27. doi: 10.3748/wjg.v21.i6.1718

PubMed Abstract | Crossref Full Text | Google Scholar

8. Turner AG, Anderson PH, Morris HA. Vitamin D and bone health. Scand J Clin Lab Invest Suppl. (2012) 243:65–72.

Google Scholar

9. Ebeling PR. Vitamin D and bone health: epidemiologic studies. Bonekey Rep. (2014) 3:511. doi: 10.1038/bonekey.2014.6

PubMed Abstract | Crossref Full Text | Google Scholar

10. Kwok RM, Torres DM, Harrison SA. Vitamin D and nonalcoholic fatty liver disease (Nafld): is it more than just an association? Hepatology. (2013) 58:1166–74. doi: 10.1002/hep.26390

PubMed Abstract | Crossref Full Text | Google Scholar

11. Szymczak-Pajor I, Drzewoski J, Sliwinska A. The molecular mechanisms by which vitamin D prevents insulin resistance and associated disorders. Int J Mol Sci. (2020) 21:6644. doi: 10.3390/ijms21186644

PubMed Abstract | Crossref Full Text | Google Scholar

12. Jones G, Prosser DE, Kaufmann M. The activating enzymes of vitamin D metabolism (25-and 1α-Hydroxylases). In: Vitamin D. London: Elsevier (2018). p. 57–79.

Google Scholar

13. Holick MF. Vitamin D deficiency. N Engl J Med. (2007) 357:266–81. doi: 10.1056/NEJMra070553

PubMed Abstract | Crossref Full Text | Google Scholar

14. Trimarco V, Manzi MV, Mancusi C, Strisciuglio T, Fucile I, Fiordelisi A, et al. Insulin resistance and vitamin D deficiency: a link beyond the appearances. Front Cardiovasc Med. (2022) 9:859793. doi: 10.3389/fcvm.2022.859793

PubMed Abstract | Crossref Full Text | Google Scholar

15. Mozos I, Marginean O. Links between vitamin D deficiency and cardiovascular diseases. Biomed Res Int. (2015) 2015:109275. doi: 10.1155/2015/109275

PubMed Abstract | Crossref Full Text | Google Scholar

16. Al-Sumaih I, Johnston B, Donnelly M, O'Neill C. The relationship between obesity, diabetes, hypertension and vitamin D deficiency among saudi arabians aged 15 and over: results from the saudi health interview survey. BMC Endocr Disord. (2020) 20:81. doi: 10.1186/s12902-020-00562-z

PubMed Abstract | Crossref Full Text | Google Scholar

17. Gradillas-García A, Álvarez J, Rubio JA, de Abajo FJ. Relationship between vitamin D deficiency and metabolic syndrome in adult population of the community of madrid. Endocrinología y Nutrición. (2015) 62:180–7. doi: 10.1016/j.endoen.2015.04.004

PubMed Abstract | Crossref Full Text | Google Scholar

18. Chung GE, Kim D, Kwak MS, Yang JI, Yim JY, Lim SH, et al. The serum vitamin D level is inversely correlated with nonalcoholic fatty liver disease. Clin Mol Hepatol. (2016) 22:146–51. doi: 10.3350/cmh.2016.22.1.146

PubMed Abstract | Crossref Full Text | Google Scholar

19. Ciardullo S, Muraca E, Cannistraci R, Perra S, Lattuada G, Perseghin G. Low 25 (Oh) vitamin D levels are associated with increased prevalence of nonalcoholic fatty liver disease and significant liver fibrosis. Diabetes Metab Res Rev. (2023) 39:e3628. doi: 10.1002/dmrr.3628

PubMed Abstract | Crossref Full Text | Google Scholar

20. Cho YH, Kim JW, Shim JO, Yang HR, Chang JY, Moon JS, et al. Association between vitamin D deficiency and suspected nonalcoholic fatty liver disease in an adolescent population. Pediatr Gastroenterol Hepatol Nutr. (2019) 22:233–41. doi: 10.5223/pghn.2019.22.3.233

PubMed Abstract | Crossref Full Text | Google Scholar

21. Gutiérrez-Cuevas J, Santos A, Armendariz-Borunda J. Pathophysiological molecular mechanisms of obesity: a link between Mafld and Nash with cardiovascular diseases. Int J Mol Sci. (2021) 22:11629. doi: 10.3390/ijms222111629

PubMed Abstract | Crossref Full Text | Google Scholar

22. Li H, Xiao P, Cheng H, Zhao X, Yan Y, Liu J, et al. Central body fat deposits are associated with poor vitamin D status in Chinese children and adolescents. Nutrition. (2022) 99–100:111651. doi: 10.1016/j.nut.2022.111651

PubMed Abstract | Crossref Full Text | Google Scholar

23. Gofton C, Upendran Y, Zheng MH, George J. Mafld: how is it different from Nafld? Clin Mol Hepatol. (2023) 29:S17–S31. doi: 10.3350/cmh.2022.0367

PubMed Abstract | Crossref Full Text | Google Scholar

24. National Center for Health Statistics. National Health and Nutrition Examination Survey. (2023). Available online at: http://www.cdc.gov/nchs/nhanes/ (Accessed September 1, 2023).

Google Scholar

25. CDC. Laboratory Procedure Manual. (2023). Available online at: https://wwwn.cdc.gov/nchs/data/nhanes/public/2017/labmethods/VID-J-MET-508.pdf (Accessed September 1, 2023).

Google Scholar

26. Karlas T, Petroff D, Sasso M, Fan JG Mi YQ, de Lédinghen V, et al. Individual patient data meta-analysis of controlled attenuation parameter (Cap) technology for assessing steatosis. J Hepatol. (2017) 66:1022–30. doi: 10.1016/j.jhep.2016.12.022

PubMed Abstract | Crossref Full Text | Google Scholar

27. Eslam M, Alkhouri N, Vajro P, Baumann U, Weiss R, Socha P, et al. Defining paediatric metabolic (dysfunction)-associated fatty liver disease: an international expert consensus statement. Lancet Gastroenterol Hepatol. (2021) 6:864–73. doi: 10.1016/S2468-1253(21)00183-7

PubMed Abstract | Crossref Full Text | Google Scholar

28. World Health Organization Bmi-for-Age Charts (5-19years) (2022). Available online at: https://www.who.int/tools/growth-reference-data-for-5to19-years/indicators/bmi-for-age (Accessed February 21, 2022).

Google Scholar

29. Xi B, Zong X, Kelishadi R, Litwin M, Hong YM, Poh BK, et al. International waist circumference percentile cutoffs for central obesity in children and adolescents aged 6 to 18 years. J Clin Endocrinol Metab. (2020) 105:e1569–83. doi: 10.1210/clinem/dgaa578

PubMed Abstract | Crossref Full Text | Google Scholar

30. American Diabetes Association Professional Practice Committee. 2. classification and diagnosis of diabetes: standards of medical care in diabetes-2022. Diabetes Care (2022) 45:S17–s38. doi: 10.2337/dc22-S002

PubMed Abstract | Crossref Full Text | Google Scholar

31. van Buuren S, Groothuis-Oudshoorn K. Mice: multivariate imputation by chained equations in R. J Statist Softw. (2011) 45:1–67. doi: 10.18637/jss.v045.i03

Crossref Full Text | Google Scholar

32. VanderWeele TJ, Ding P. Sensitivity analysis in observational research: introducing the E-value. Ann Intern Med. (2017) 167:268–74. doi: 10.7326/M16-2607

PubMed Abstract | Crossref Full Text | Google Scholar

33. Black LJ, Jacoby P, She Ping-Delfos WC, Mori TA, Beilin LJ, Olynyk JK, et al. Low serum 25-hydroxyvitamin D concentrations associate with non-alcoholic fatty liver disease in adolescents independent of adiposity. J Gastroenterol Hepatol. (2014) 29:1215–22. doi: 10.1111/jgh.12541

PubMed Abstract | Crossref Full Text | Google Scholar

34. Teixeira JS, Bull Ferreira Campos A, Cordeiro A, Pereira SE, Saboya CJ, Ramalho A. Vitamin D nutritional status and its relationship with metabolic changes in adolescents and adults with severe obesity. Nutr Hosp. (2018) 35:847–53. doi: 10.20960/nh.1657

PubMed Abstract | Crossref Full Text | Google Scholar

35. Barchetta I, Cimini FA, Cavallo MG. Vitamin D and metabolic dysfunction-associated fatty liver disease (Mafld): an update. Nutrients. (2020) 12:3302. doi: 10.3390/nu12113302

PubMed Abstract | Crossref Full Text | Google Scholar

36. Bozic M, Guzman C, Benet M, Sanchez-Campos S, Garcia-Monzon C, Gari E, et al. Hepatocyte vitamin D receptor regulates lipid metabolism and mediates experimental diet-induced steatosis. J Hepatol. (2016) 65:748–57. doi: 10.1016/j.jhep.2016.05.031

PubMed Abstract | Crossref Full Text | Google Scholar

37. Shaheen M, Schrode KM, Pan D, Kermah D, Puri V, Zarrinpar A, et al. Sex-specific differences in the association between race/ethnicity and Nafld among us population. Front Med. (2021) 8:795421. doi: 10.3389/fmed.2021.795421

PubMed Abstract | Crossref Full Text | Google Scholar

38. Lee EY, Nguyen VH, Cheung R, Nguyen MH. Trends of chronic liver diseases by income level and socioeconomic factors in the united states: a population-based study. Aliment Pharmacol Ther. (2024) 60:1374–87. doi: 10.1111/apt.18242

PubMed Abstract | Crossref Full Text | Google Scholar

39. Zhang JJ Yu HC, Li Y, Zhang YB, Geng TT, Lu Q, et al. Association between Serum 25-hydroxy vitamin D concentrations and mortality among individuals with metabolic dysfunction-associated fatty liver disease: a prospective cohort study. Am J Clin Nutr. (2022) 116:1409–17. doi: 10.1093/ajcn/nqac260

PubMed Abstract | Crossref Full Text | Google Scholar

40. Pires LV, Gonzalez-Gil EM, Anguita-Ruiz A, Bueno G, Gil-Campos M, Vazquez-Cobela R, et al. The vitamin D decrease in children with obesity is associated with the development of insulin resistance during puberty: the pubmep study. Nutrients. (2021) 13:4488. doi: 10.3390/nu13124488

PubMed Abstract | Crossref Full Text | Google Scholar

41. Scapaticci S, D'Adamo E, Mohn A, Chiarelli F, Giannini C. Non-alcoholic fatty liver disease in obese youth with insulin resistance and type 2 diabetes. Front Endocrinol. (2021) 12:639548. doi: 10.3389/fendo.2021.639548

PubMed Abstract | Crossref Full Text | Google Scholar

42. Sakurai Y, Kubota N, Yamauchi T, Kadowaki T. Role of insulin resistance in Mafld. Int J Mol Sci. (2021) 22:4156. doi: 10.3390/ijms22084156

PubMed Abstract | Crossref Full Text | Google Scholar

43. Cui X, Sun X, Li Q, Chen Z. Changes in blood glucose and lipid metabolism levels in children with central precocious puberty and its correlation with obesity. Front Pediatr. (2024) 12:1488522. doi: 10.3389/fped.2024.1488522

PubMed Abstract | Crossref Full Text | Google Scholar

44. Argano C, Mirarchi L, Amodeo S, Orlando V, Torres A, Corrao S. 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

PubMed Abstract | Crossref Full Text | Google Scholar

45. Namakin K, Hosseini M, Zardast M, Mohammadifard M. Vitamin D effect on ultrasonography and laboratory indices and biochemical indicators in the blood: an interventional study on 12 to 18-year-old children with fatty liver. Pediatr Gastroenterol Hepatol Nutr. (2021) 24:187–96. doi: 10.5223/pghn.2021.24.2.187

PubMed Abstract | Crossref Full Text | Google Scholar

46. Murat G, Sule A. Understanding the role of vitamin a and its precursors in the immune system. Nutr Clinique et Métabolisme. (2022) 36:89–98. doi: 10.1016/j.nupar.2021.10.002

Crossref Full Text | Google Scholar

47. Chen W, Chen G. The roles of vitamin a in the regulation of carbohydrate, lipid, and protein metabolism. J Clin Med. (2014) 3:453–79. doi: 10.3390/jcm3020453

PubMed Abstract | Crossref Full Text | Google Scholar

48. Petiz LL, Kunzler A, Bortolin RC, Gasparotto J, Matté C, Moreira JCF, et al. Role of vitamin A oral supplementation on oxidative stress and inflammatory response in the liver of trained rats. Appl Physiol Nutr Metab. (2017) 42:1192–200. doi: 10.1139/apnm-2017-0193

PubMed Abstract | Crossref Full Text | Google Scholar

49. Saeed A, Dullaart RPF, Schreuder T, Blokzijl H, Faber KN. Disturbed vitamin A metabolism in non-alcoholic fatty liver disease (Nafld). Nutrients. (2017) 10:29. doi: 10.3390/nu10010029

PubMed Abstract | Crossref Full Text | Google Scholar

50. Zeng Y, Luo M, Pan L, Chen Y, Guo S, Luo D, et al. Vitamin D signaling maintains intestinal innate immunity and gut microbiota: potential intervention for metabolic syndrome and Nafld. Am J Physiol Gastrointest Liver Physiol. (2020) 318:G542–G53. doi: 10.1152/ajpgi.00286.2019

PubMed Abstract | Crossref Full Text | Google Scholar

Keywords: MAFLD, vitamin D, adolescents, NHANES, retinol status

Citation: Gao X, Chen Y, Wang X, Chen Y, Chen X, Li H and Ye H (2025) Association of vitamin D levels with metabolic dysfunction-associated fatty liver disease in children aged 12–18 years. Front. Nutr. 12:1615851. doi: 10.3389/fnut.2025.1615851

Received: 22 April 2025; Accepted: 14 October 2025;
Published: 04 November 2025.

Edited by:

Getu Gizaw, Jimma University, Ethiopia

Reviewed by:

Ruchika Maurya, Medical University of South Carolina, United States
Renying Xu, Shanghai Jiao Tong University, China

Copyright © 2025 Gao, Chen, Wang, Chen, Chen, Li and Ye. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Haibo Li, aGFpYm9saTg5QDE2My5jb20=; Hong Ye, eWhmanNmeTk5QDE2My5jb20=

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.