Skip to main content

Association between iron metabolism and non-alcoholic fatty liver disease: results from the National Health and Nutrition Examination Survey (NHANES 2017–2018) and a controlled animal study

Abstract

Background

Iron metabolism may be involved in the pathogenesis of the non-alcoholic fatty liver disease (NAFLD). The relationship between iron metabolism and NAFLD has not been clearly established. This study aimed to clarify the relationship between biomarkers of iron metabolism and NAFLD.

Methods

Based on the National Health and Nutrition Examination Survey (NHANES), restricted cubic spline models and multivariable logistic regression were used to examine the association between iron metabolism [serum iron (SI), serum ferritin (SF), transferrin saturation (TSAT), and soluble transferrin receptor (sTfR)] and the risk for NAFLD. In addition, stratified subgroup analysis was performed for the association between TSAT and NAFLD. Moreover, serum TSAT levels were determined in male mice with NAFLD. The expression of hepcidin and ferroportin, vital regulators of iron metabolism, were analyzed in the livers of mice by quantitative real-time PCR (qRT-PCR) and patients with NAFLD by microarray collected from the GEO data repository.

Results

Patients with NAFLD showed decreased SI, SF, and TSAT levels and increased STfR levels based on the NHANES. After adjusting for confounding factors, TSAT was significantly negatively correlated with NAFLD. Of note, the relationship between TSAT and NAFLD differed in the four subgroups of age, sex, race, and BMI (P for interaction < 0.05). Consistently, mice with NAFLD exhibited decreased serum TSAT levels. Decreased hepcidin and increased ferroportin gene expression were observed in the livers of patients and mice with NAFLD.

Conclusion

Serum TSAT levels and hepatic hepcidin expression were decreased in both patients and mice with NAFLD. Among multiple biomarkers of iron metabolism, lower TSAT levels were significantly associated with a higher risk of NAFLD in the U.S. general population. These findings might provide new ideas for the prediction, diagnosis, and mechanistic exploration of NAFLD.

Introduction

Affecting ~ 25% of the general population worldwide, non-alcoholic fatty liver disease (NAFLD) has become the most common liver disease and is predicted to become increasingly prevalent, particularly among children and younger adults [1, 2]. NAFLD is a progressive disease that can further develop into non-alcoholic steatohepatitis (NASH), liver fibrosis, cirrhosis, and hepatocellular carcinoma, and it is always accompanied by diabetes, cardiovascular disease, and chronic kidney disease, all of which will lead to an increased risk of mortality [3, 4]. The pathogenesis of NAFLD is not fully understood, which is associated with multiplex risk factors [5,6,7]: metabolic risk factors, gut microbiome composition, genetic factors, epigenetic factors, environmental risk factors, etc.

In recent years, the impact of iron metabolism on NAFLD has attracted renewed attention due to the proposal of “ferroptosis”, an iron-dependent form of cell death [8, 9]. Elena et al. [3] discovered that variants of genes related to iron metabolism are associated with high ferritin levels and increased hepatic iron in an Italian cohort of patients with NAFLD. Jordi et al. [10] showed that iron status influences liver fat accumulation in NAFLD through the gut microbiome. And a rat model of NAFLD exhibited systemic iron deficiency and hepatic iron overload [11]. The above studies support the notion that dysregulated iron homeostasis plays a role in the pathogenesis of NAFLD. And it was reported that a disturbance in iron metabolism affects one-third of patients with NAFLD [12]. However, the association between iron metabolism including serum levels of iron (SI), ferritin (SF), transferrin saturation (TSAT), and soluble transferrin receptor (sTfR) and NAFLD in general population has not been clearly established. Also, whether those biomarkers of iron metabolism could provide certain diagnostic value for NAFLD remains unknown.

Iron, a critical part of the hemoglobin of human beings [12], is essential for oxygen transportation, energy formation, and many cellular functions including DNA synthesis and repair [13, 14]. Iron cannot be naturally created by the human body, but must be ingested through diet or supplements [15]. For men and postmenopausal women, the recommended dietary allowance (RDA) of iron is 8 mg/day, whereas for premenopausal women, it is 18 mg/day according to the Institute of Medicine (US) (Institute of Medicine, 2001). Dietary iron is absorbed mainly in the duodenum and upper jejunum but is not enough for daily needs. The recycled iron from senescent or damaged erythrocytes in the spleen, liver, and bone marrow also contributes to iron storage [16]. Hepcidin is the master regulator for maintaining iron homeostasis, inducing degradation and internalization of the iron exporter ferroportin to inhibit the release of iron from recycling macrophages and absorption from dietary sources, thereby indirectly reducing iron entry into the bloodstream [17, 18]. The liver is the major site of iron storage [19] and plays a central part in iron homeostasis by producing hepcidin [18]. Therefore, liver diseases such as NAFLD may be more susceptible to iron status.

The movement of iron between different cells and body tissues is mainly conducted by the transport protein transferrin, which plays a central role in iron metabolism. Each transferrin protein has two iron-binding domains that can reversibly bind two atoms of iron in a soluble nontoxic form [20, 21]. Serum iron (SI) consists of mostly transferrin-bound iron and a negligible amount of free non-transferrin-bound iron, which is toxic and was reported to promote oxidative stress [22, 23]. Ferritin is a major iron storage protein [24], and serum ferritin (SF) is the small quantity of ferritin circulating in the blood, which is regarded as a marker of iron stores in healthy individuals and those with early iron deficiency (ID). Transferrin saturation (TSAT), which indicates how many transferrin iron-binding sites are occupied, is considered an important biochemical marker of body iron status [25]. sTfR is a soluble form of transferrin receptor identified in serum [26] that reflects the demand for iron in cells and increases rapidly with the depletion of stored iron in the early stage of ID. Four biomarkers, SI, SF, TSAT, and sTfR, are indicators of iron metabolism and are commonly measured in clinics [25, 26].

Although several epidemiological studies have explored associations between biomarkers of iron metabolism and NAFLD, the conclusions were inconclusive and inconsistent. Jung et al. found that serum ferritin levels were positively associated with liver steatosis and fibrosis in the Korean general population [27]. Yang et al. [28] discovered that higher serum iron levels decreased the risk of NAFLD. And to the best of our knowledge, most studies have only evaluated the relationship between one index of iron metabolism and NAFLD. In this study, we comprehensively investigated the association between biomarkers of iron metabolism, including SI, SF, TSAT, and sTfR, and the prevalence of NAFLD in the general population of the U.S., utilizing data from the National Health and Nutrition Examination Survey (NHANES) collected from 2017 to 2018. And we have found that lower TSAT levels were significantly associated with a higher risk of NAFLD, which might provide additional information on biomarkers for the diagnosis of NAFLD. In addition, a controlled animal study was conducted to verify the results of the NHANES.

Materials and methods

Population-based human study (NHANES)

Study population

The NHAENES database provides a large, sophisticated, stratified, ongoing analysis of nutrition and health data for the entire U.S. population [29]. NHANES data (2017–2018) were used to investigate the relationship between iron metabolism and NAFLD. Those participants who lacked serum iron (SI), serum ferritin (SF), transferrin saturation (TSAT), soluble transferrin (sTfR), and NAFLD data were excluded from this study (n = 918). There was a total of 5483 individuals considered. The National Center for Health Statistics Research Ethics Review Board approved all protocols and informed consent was obtained from all subjects. More information on the survey design, methodology, and data can be found on the NHANES website (https://www.cdc.gov/nchs/nhanes/).

NAFLD outcomes

NAFLD was defined using the U.S. fatty liver index (FLI), a well-validated diagnostic index [30], which was employed utilizing NHANES III data and calculated as an equation according to a previous study [31, 32] that included information on body mass index (BMI), gamma glutamyl transferase (GGT), triglycerides (TG), and waist circumference. All the information was collected concurrently with the status of iron metabolism. NAFLD was defined as an FLI score of ≥ 60. The FLI formula is expressed as follows:

$${\text{FLI}} = ({\text{e}}^{{0.953 * \ln ({\text{TG}}) + 0.139 * {\text{BMI}} + 0.718 * \ln ({\text{GGT}}) + 0.053 * {\text{waist circumference}}^{ - 15.745} }} )/(1 + {\text{e}}^{{0.953 * \ln ({\text{TG}}) + 0.139 * {\text{BMI}} + 0.718 * \ln ({\text{GGT}}) + 0.053 * {\text{waist circumference}}^{ - 15.745} }} ) * 100.$$

SI, SF, sTfR, and TSAT measurements

The serum iron concentration was measured using the DCX-800 system, which is a timed endpoint method [33]. Acetic acid releases iron from transferrin during the process, and hydroxylamine and thioglycolate reduce it to the ferrous state. FerroZine Iron Reagent quickly complexes with ferrous irons. At a fixed-time interval, the system tracks changes in absorbance at 560 nm. The concentration of iron in the sample is exactly proportional to this change in absorbance. Using the Roche/Hitachi 912 Clinical Analyzer, immunoturbidimetry was the technique of choice for measuring ferritin [34]. Antigen/antibody complexes were created when latex-bound ferritin antibodies interacted with the antigen in the sample. This was turbidimetrically measured after agglutination. The complexes produced were measured at 700 nm and were proportional to the ferritin content (primary wavelength). The method principle for the measurement of soluble transferrin receptor (sTfR) was a particle-enhanced immunoturbidimetric assay that used Roche kits for the Cobas® c501 clinical analyzer [35]. The antigen in the sample reacted with latex particles coated with anti-sTfR antibodies to generate an antigen/antibody combination. The precipitate was identified photometrically after agglutination. Total iron binding capacity (TIBC) was calculated indirectly using the unsaturated iron binding capacity (UIBC) method [36], and the transferrin saturation value was calculated as (iron/TIBC) × 100%. The biomarkers of iron metabolism (SI, SF, TSAT, and sTfR) were divided into quartiles for the assessment of the possible association between iron metabolism and incident NAFLD. The instrument measures SI, SF, TSAT, and sTfR in a range of 11–481 µg/dL, 2–1090 µg/L, 1.4–95% and 1–34.2 mg/L, respectively. Those lower than the detection limits value were missing values and were deleted in this study. The NHANES website (https://wwwn.cdc.gov/nchs/data/nhanes/2017-2018/manuals/2017_MEC_Laboratory_Procedures_Manual.pdf) describes a more thorough processing process for SI, SF, TSAT, and sTfR.

Covariates

This research included a number of covariates: demographic data (age, sex, race/ethnicity, family poverty income ratio (PIR), marital status, education level), dietary data (mean energy intake, protein intake, folic acid intake, vitamin B12 intake, vitamin C intake, iron intake), questionnaire data (hypertension, diabetes mellitus (DM), coronary heart disease (CHD), congestive heart failure (CHF), angina pectoris, heart attack, stroke, smoking, alcohol use, and physical activity (PA, which was collected from the Physical Activity questionnaire (PAQ) in NHANES and categorized into four groups according the intensity of PA: No, Moderate, Both (participants who had a combination of moderate-intensity and vigorous-intensity physical activity,), and Vigorous)), examination data (BMI and waist circumference), laboratory data (total cholesterol (TC), triglycerides (TG), high-density lipoprotein cholesterol (HDL cholesterol, HDL-C) and haemoglobin (HB), urinary albumin creatinine ratio (uACR), alanine aminotransferase (ALT), aspartate aminotransferase (AST), GGT, fasting glucose, fasting insulin, glycosylated haemoglobin (HbA1c), high-sensitivity C-reactive protein (hsCRP), estimated glomerular filtration rate (eGFR), uric acid (UA), blood urea nitrogen (BUN), and serum creatinine (Scr)). Detailed covariate information is available publicly from the NHANES database (http://www.cdc.gov/nchs/nhanes/).

Statistical analysis

The NHANES estimations were all based on sample weights [37]. All analyses were performed in version 3.6.4 of R (R Foundation for Statistical Computing, Vienna, Austria) and version 22.0 of SPSS (SPSS Inc., Chicago, IL, USA). Continuous variables are presented as the means ± standard deviations, and categorical variables are expressed as numbers (n) and percentages (%). To investigate the relationship between iron metabolism and NAFLD, multivariable logistic regression was used. First, Model 1 was adjusted for age and sex. Second, based on Model 1, race/ethnicity, level of education, marital status, family PIR, hypertension, DM, smoking status, and drinking status were further adjusted for (Model 2). Finally, Model 3 was updated as our main model and included Model 2 variables plus BMI, waist circumference, PA, mean energy intake, protein intake, folic acid intake, vitamin B12 intake, vitamin C intake, iron intake, the complication of CHD, CHF, angina pectoris, heart attack, stroke, and TC, TG, HDL-C, HB, HbA1c, hsCRP, ALT, AST, GGT, eGFR, uACR, UA, BUN, and Scr. Subgroup analyses were used to evaluate the relationship between iron status and NAFLD based on age, sex, race/ethnicity, hypertension, DM, and BMI. We used regularization technique (Least absolute shrinkage and selection operator (LASSO) regression) to solve the potential overfitting problem. In addition, multivariate stepwise regression analysis was performed to deal with multiple testing problem. P-value < 0.05 was considered statistically significant.

Laboratory-based controlled animal study

Five-week-old male C57BL/6J mice were purchased from Beijing Vital River Laboratory Animal Technology Co., Ltd., treated according to the guidelines of the China Pharmaceutical University Institutional Animal Care and Use Committee (approval number: 2019-03-001), and all mouse studies were reported according to the ARRIVE guidelines [38]. Mice were housed in a specific pathogen-free (SPF) facility and adapted to the new conditions for 2 days before the experiment. Mice were maintained in a temperature-controlled (22–23 °C) room with a 12:12-h light/dark cycle.

Animals and experimental protocols

To induce NAFLD, mice were fed a high-fat diet (D12492, 60% kcal from fat; research diet) and fructose (Y0002132, 2.31 g/100 ml; sigma) water (high-fat-fructose diet) for 30 weeks (n = 5). Mice fed normal chow served as controls (n = 5). When the mice were sacrificed, blood was collected from the inferior vena cava, and the entire middle lobe of the liver was fixed with formalin. The largest lobe in the middle of the liver was divided into 8 small pieces after cutting the edge, and one of those 8 pieces was directly added to TRIzol (9109, TaKaRa, Dalian, Liaoning, China) for RNA extraction.

Assays for biomarkers of iron metabolism

Transferrin saturation (TSAT) was defined as the ratio of serum iron (SI) and total iron-binding capacity (TIBC). Unsaturated iron-binding capacity (UIBC) is calculated by subtracting SI from TIBC [39]. SI and TIBC were measured using commercial kits (all from Nanjing Jiancheng Bioengineering Institute) according to standard procedures [11, 40]. Serum ferritin, sTfR, and transferrin were determined using ELISA kits (Jiangsu Meibiao Biotechnology Co., Ltd) according to the manufacturer’s instruction [41]. Serum ferrous irons (Fe2+) serum was measured using the phenanthroline colorimetric method [40] by colorimetric assay kit purchased from Elabscience following the manufacturer’s instructions.

Measurements of hepatic TG, TC and MDA

Three pieces of liver were selected for the quantification of triglyceride (TG), total cholesterol (TC), and malonic dialdehyde (MDA) respectively, which were determined by commercial kits (Nanjing Jiancheng Bioengineering Inst ferrous iron itute) according to the manufacturer’s protocols as previously described [5, 42]. Thiobarbituric acid (TBA) method was used for the determination of MDA, which was based on that MDA can react with TBA at high temperature and acidity to produce the red-brown product with the maximum absorption peak at 532 nm [42].

Histological analysis

The fixed liver lobes were embedded in paraffin, then sectioned at 5 µm and stained with hematoxylin–eosin (H&E) according to the standard protocol described [43,44,45]. Scores for steatosis, inflammation, and ballooning were assessed by a four-member research team according to the NAFLD Activity Score (NAS) system [46].

RNA extraction and quantitative real-time PCR (qRT‒PCR) analysis

Total RNA extraction and cDNA reversion were conducted according to methods described previously [5, 7, 47]. qRT-PCR was performed using SYBR Premix (Vazyme, Q331-02, China) according to the manufacturer’s instructions for the ABI StepOnePlus real-time PCR machine (Applied Biosystems). Custom primers were designed for mouse hepcidin (forward, 5′-TTGCGATACCAATGCAGAAGA-3′; reverse, 5′-GATGTGGCTCTAGGCTATGTTTTG-3′), ferroportin (forward, 5′-TTGCAGGAGTCATTGCTGCTA-3′; reverse: 5′-TGGAGTTCTGCACACCATT-3′), and 18S (forward, 5′-AGTCCCTGCCCTTTGTACACA-3′; reverse, 5′-CGATCCCAGGGCCTCACTA-3′). All real-time PCR primers were synthesized by Sangon Biotech (Shanghai, China). Gene expression was normalized to 18S and calculated as previously described [5, 7, 47].

Statistical analysis

The results are presented as the means ± SEMs and were statistically analysed by the unpaired Student's t test or Pearson product-moment correlation coefficients. Differences were considered significant when the P-value was < 0.05.

Microarray and RNA-seq data

Hepatic gene expression was compared among 20 patients with simple steatosis, 19 with non-alcoholic steatohepatitis (NASH), and 24 healthy controls in the GSE89632 dataset [48] (https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE89632). In the GSE126848 dataset, RNA sequencing was performed on liver biopsies obtained from healthy normal weight (n = 14) and healthy obese (n = 12) individuals, and simple steatosis (n = 15) and NASH (n = 16) patients [49] (https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE126848). In the GSE185051 dataset, liver biopsies from 51 paediatric NASH patients and 5 normal subjects were analysed by RNA sequencing [50] (https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE185051). 12 biopsy diagnosed NASH patients and 5 different subjects for healthy control groups were included in the GSE24807 dataset [51] (https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE24807).

To obtain these data, we used the package 'GEOquery' installed in R, especially the function getGEO(),  importing the dataset GSE number, and obtained the list object of the corresponding dataset. Then, the pData() and expr() functions were used to view the sample information and expression matrix, respectively. The chip platform used in the dataset was confirmed by viewing the web page and checking the list of objects. The gene annotation information of the platform was obtained by using the GEOquery software package and matched with the array information. The gene symbol was used to replace the probe name. Genes with |log2-fold changes (FCs)|> 0.5 and P-value < 0.05 were considered differentially expressed genes (DEGs). Then, the hepcidin and ferroportin genes were selected, and the expression information of these genes in each sample and the grouping information of each sample were extracted and stored in an Excel table for subsequent analysis. The differences between different groups were compared by the Mann–Whitney test and considered statistically significant at the P-value < 0.05.

Results

Population-based human study (NHANES)

Baseline characteristics

The population-weighted characteristics of this study are presented in Table 1. The final analysis consisted of 5483 individuals after the removal of people with missing iron metabolism index (SI, SF, TSAT, and sTfR) and NAFLD information. The average age of the entire population was 44.425 ± 0.622 years. The prevalence of NAFLD was 35.19% in this group. In addition, the SI, SF, TSAT, and sTfR averages among these 5483 individuals were 89.073 ± 1.107 µg/dL, 138.120 ± 3.320 µg/L, 27.489 ± 0.291%, and 3.211 ± 0.047 mg/L, respectively. Compared with participants without NAFLD, participants with NAFLD showed significantly decreased levels of SI (P < 0.001) and TSAT (P < 0.001) and a significantly increased level of sTfR (P < 0.001) (Table 1). However, SF levels did not differ significantly between the two groups. (P = 0.342), even though it was slightly decreased in NAFLD groups on average. Among 918 participants with missing data, 12.5% were Mexcan American, 9.6% were Other Hispanic, 25.2% were Non-Hispanic Black, 32.4% were Non-Hispanic White, and 20.4% were Other race. Those under the age of 60 accounted for 67.7%, and individuals 60 years of age or older accounted for 32.3%. In addition, those were male accounted for 46.1%, and individuals were female accounted for 53.9%. Finally, we also compared the characteristics of the populations between the those with missing values and those without in Additional file 5: Table S1.

Table 1 Characteristics of the study population

In healthy individuals, the normal ranges SI, SF, and TSAT in adults are 55–185 µg/dL, 15–250 µg/L and 14.4–46.8% respectively [52, 53]. As for the normal range of sTfR, different measurement methods provide different values due to the lack of uniform standards [26, 54]. The iron status index was calculated and divided into quartiles in this study. The SF ranged from < 42.50, 42.50 to 93.30, 93.31 to 179.00, and > 179.00 in quartile 1 (Q1), quartile 2 (Q2), quartile 3 (Q3), and quartile 4 (Q4), respectively; the SI ranged from < 63.00, 63.00 to 84.00, 84.01 to 108.00, and > 108.00 in quartile 1 (Q1), quartile 2 (Q2), quartile 3 (Q3), and quartile 4 (Q4), respectively; the TSAT ranged from < 19.00, 19.00 to 26.00, 26.01 to 34.00, and > 34.00 in quartile 1 (Q1), quartile 2 (Q2), quartile 3 (Q3), and quartile 4 (Q4), respectively; and the sTfR ranged from < 2.49, 2.49 to 2.99, 3.00 to 3.68, and > 3.68 in quartile 1 (Q1), quartile 2 (Q2), quartile 3 (Q3), and quartile 4 (Q4), respectively (Table 2).

Table 2 Adjusted ORs for associations between iron metabolism and the prevalence of NAFLD

Association between iron status and NAFLD

The potential relationship between each iron status index and NAFLD risk was examined using restricted cubic splines with a three-knot model. There was a linear and negative correlation between SI, SF, as well as TSAT and NAFLD (Fig. 1a–c), whereas there was a linear positive correlation between sTfR and NAFLD (Fig. 1d). In addition, the correlation heatmap of biomarkers of iron metabolism and NAFLD revealed that SI was negatively related with NAFLD (r =  − 0.13), SF was negatively related with NAFLD (r =  − 0.01), TSAT was negatively related with NAFLD (r =  − 0.13), and sTfR were positively related with NAFLD (r = 0.06) (Additional file 1: Fig. S1). Table 2 illustrates the results of the multivariate logistic regression analyses for the iron metabolism index (SI, SF, TSAT, and sTfR) and NAFLD. After adjustment for confounding factors, compared with people in the Q1 group, the odds ratios (ORs) with 95% confidence intervals (CIS) for NAFLD were 1.007 (0.844, 1.203), 0.994 (0.825, 1.197), and 0.842 (0.688, 1.030) for SI; 1.093 (0.905, 1.320), 1.026 (0.837, 1.258), and 0.973 (0.779, 1.214) for SF; 0.905 (0.758, 1.079), 0.987 (0.821, 1.186), and 0.706 (0.575, 0.868) for TSAT; and 1.031 (0.961, 1.235), 1.096 (0.916, 1.312), and 1.151 (0.955, 1.387) for sTfR in the Q2, Q3 and Q4 groups, respectively. In multivariate logistic regression model, TSAT levels were significantly inversely correlated with NAFLD. In addition, iron metabolism markers (SI, SF, TSAT, and sTfR) were used as continuous variables to explore their association with NAFLD (Additional file 6: Table S2). A total of 43 features were reduced to 13 potential predictors, including TSAT, on the basis of 5483 individuals, and were features with nonzero coefficients in the LASSO regression model (Additional file 2: Fig. S2a, b). We used the other 12 features (age, BMI, waist, Alt, Ast, HDL-C, HbA1c, protein intake, sex, education level, DM, and alcohol user) as covariates screened by LASSO regression analysis to explore the association between TSAT and NAFLD. The RCS plot also showed a liner and negative relationship between TSAT and NAFLD (Additional file 3: Fig. S3a). Additionally, by using multivariate stepwise regression analysis, we reduced 43 features to 13 potential predictors, including TSAT (Additional file 7: Table S3). Then, we used the other 12 covariates (age, BUN, smoker, waist, Alt, HbA1C, sex, TC, alcohol user, HDL-C, education level, and dietary Vitamin C intake) as covariates screened by multivariate stepwise regression analysis to explore the relationship between TATS and NAFLD. The RCS plot also showed a liner and negative relationship between TSAT and NAFLD (Additional file 3: Fig. S3b).

Fig. 1
figure 1

Restricted cubic spline plots of the association between the indicators of iron metabolism and the risk of NAFLD based on the NHANES database. The levels of serum iron (a), ferritin (b), TSAT (c), and sTfR (d) were log10 transformed. Analyses were adjusted for age, sex, education level, race/ethnicity, family poverty-income ratio, smoking, hypertension, diabetes mellitus, body mass index, waist circumference, physical activity, total energy intake, high-density lipoprotein-cholesterol, triglycerides, total cholesterol, uric acid, and hemoglobin. Solid and dashed lines represent the log-transformed odds ratios and the corresponding 95% confidence intervals

Subgroup analysis

Subgroup analyses stratified by age, sex, race, hypertension, DM, and BMI were further conducted (Table 3). As shown in Table 3, the association between TSAT and NAFLD was statistically significant in those subgroups of individuals aged < 60 years (P for trend = 0.009), female (P for trend = 0.019), Mexican American (P for trend = 0.018), participants without hypertension (P for trend = 0.003) or DM (P for trend = 0.013), and individuals with a BMI of ≥ 30 kg/m2 (P for trend = 0.033). Notably, the association differed in the four subgroups of age, sex, race, and BMI (P for interaction < 0.05).

Table 3 Subgroups analysis for the associations of TSAT with the prevalence of NAFLD

Laboratory-based controlled animal study

Since we found a significant inverse relationship between TSAT and NAFLD based on the NHANES survey, we constructed a high-fat-fructose diet NAFLD model in an animal study, which most closely resembled human NAFLD [55]. According to the statistical information in the review article [55], we selected male mice for the experiment, excluding the unstable factors that may be caused by estrogen or progesterone in female mice. Compared with mice fed a normal chow diet, mice fed a high-fat-fructose diet for 30 weeks showed significant weight gain, as well as significant increases in TG and TC levels in the liver (Additional file 4: Fig. S4a,b). Biopsy of liver tissue is the gold standard for the diagnosis of NAFLD [56]. As shown in Fig. 2a, liver sections from mice fed a high-fat-fructose diet exhibited more substantial liver fat accumulation, inflammatory cell infiltration, and hepatocyte distension, as quantified by NAS scores (Fig. 2b), compared with mice fed a normal chow diet. These results demonstrated the successful establishment of the NAFLD model in mice. We then determined the level of serum TSAT and found that it was significantly decreased in mice with NAFLD (Fig. 2c), which was consistent with the NHANES database. However, SI levels did not show the same downwards trend among mice with NAFLD as in the NHANES database. For the assessment of iron status in human, SI was not very informative because SI values vary greatly and even fluctuate hourly [11]. Both TIBC and UIBC were significantly elevated, as shown in Fig. 2d, which associated with a decrease in TSAT levels (as TSAT was defined as the ratio of SI in TIBC). Next, we evaluated serum ferritin and sTfR in the mouse model by ELISA and found that neither of them had significant changes (Fig. 2e–f). However, serum transferrin was significantly increased, which was consistent with the increase in TIBC (Fig. 2g). In addition, we measured serum ferrous irons (Fe2+) and found it was increased significantly in mice with NAFLD (Fig. 2h). Fe2+ is the key to producing the Fenton reaction [57], which might be related to excessive oxidative stress and liver injure proved by increased hepatic MDA levels in NAFLD (Additional file 4: Fig. S4c).

Fig. 2
figure 2

Decreased serum TSAT and hepatic hepcidin expression were shown in mice with NAFLD. aj Mice were fed a normal chow or high-fat-fructose diet for 30 weeks (n = 5 for each group). a Representative H&E staining of the liver sections. Scale bars: 50 µm. b NAS scores for quantifying H&E staining. c Levels of serum TSAT. d Levels of serum SI, TIBC and UIBC. e Levels of serum ferritin. f Levels of serum sTfR. g Levels of serum transferrin. h Levels of serum ferrous iron (Fe2+). i Relative mRNA levels of hepcidin and ferroportin in the livers of the mice in the indicated groups. j Correlation of liver hepcidin expression with serum TSAT levels

Hepcidin is a major regulator of iron metabolism, mainly by regulating the degradation of the iron exporter ferroportin [17, 18]. Under iron deficiency, the expression of hepcidin is reduced to ensure that more iron enters the blood and maintains iron balance [58]. We evaluated the mRNA expression of hepcidin in the liver and found that it was significantly decreased, while that of ferroportin was increased (Fig. 2i). The correlation between TSAT and hepcidin is not very clear [59, 60]; we examined this association and found that TSAT and hepcidin were positively correlated (Fig. 2j).

Microarray and RNA-seq data

Next, we investigated the transcriptional levels of hepcidin and ferroportin in the livers of NAFLD patients using Gene Expression Omnibus (GEO) datasets. Compared with healthy individuals, patients with NAFLD had decreased hepatic hepcidin (Fig. 3a) and increased hepatic ferroportin (Fig. 3b) gene expression [48,49,50,51].

Fig. 3
figure 3

Human hepcidin and ferroportin gene expression acquired from the Gene Expression Omnibus datasets. a Human hepcidin expression in the GSE89632 and GSE126848 databases. b Human ferroportin expression in the GSE89632, GSE185051 and GSE24807 databases

In summary, the current study comprehensively investigated the association between biomarkers of iron metabolism, including SI, SF, TSAT, and sTfR, and the prevalence of NAFLD in the general U.S. population. We found that decreased levels of SI, SF, and TSAT were associated with a high risk of NAFLD, while increased sTfR levels were associated with an elevated prevalence of NAFLD. We then used multivariate logistic regression analyses to assess the association between these four indicators and NAFLD and found that only the correlation between TSAT and NAFLD was statistically significant. The association was statistically significant among individuals younger than 60 years, woman, Mexican American, and those who had a BMI ≥ 30 kg/m2 (P for trend < 0.05). In addition, a controlled animal study was conducted to verify the results of the NHANES. Compared with mice fed normal chow, mice with NAFLD fed a high-fat-fructose diet showed decreased serum TSAT levels. Moreover, we analyzed the expression of hepcidin and ferroportin, vital regulators of iron metabolism, and observed decreased hepcidin and increased ferroportin gene expression in the livers of patients and mice with NAFLD. We clearly elucidated the relationship between iron metabolism indicators and NAFLD.

Discussion

The current study investigated the association between iron metabolism and the prevalence of NAFLD in a large-scale American population and a controlled animal study. Based on the NHANES database, decreased serum TSAT levels were found in patients with NAFLD and were significantly correlated with a higher risk of NAFLD. A subgroup analysis for the correlation between TSAT and NAFLD stratified by age, sex, hypertension, DM, and BMI was further performed and exhibited different associations in the four subgroups. Consistent with the results from the NHANES, mice with NAFLD also exhibited decreased serum TSAT levels, which was correlated with reduced liver hepcidin gene expression. Additionally, microarray data showed that hepcidin is decreased at the transcriptional level in NAFLD patients. These findings indicate the dysregulation of iron metabolism in NAFLD and might provide additional information on biomarkers for the diagnosis of NAFLD in the general population.

Iron, as an essential trace element in human beings, is abundant on earth. However, the bioavailability of iron is hampered by its ability to form highly insoluble oxides, with iron deficiency being the most common nutritional disorder, affecting billions of people worldwide [61]. In contrast to the high prevalence of iron deficiency, iron overload is less prevalent but is often observed among patients with NAFLD [62]. Excess iron is toxic, as it exacerbates oxidative stress and severely damages cells and tissues via the Fenton reaction [63]. Among the four indicators (SI, SF, TSAT and sTfR) of iron metabolism, TSAT was significantly decreased both in mice and patients with NAFLD. The mRNA level of the hepatic hepcidin, key regulator of iron metabolism, was also significantly decreased in mice and patients with NAFLD. The changes in TSAT and hepcidin indicated the dysregulation of iron metabolism in NAFLD. Among NAFLD patients of the NHANES, the average value of TSAT was 24.897 ± 0.378, which was still within the normal range. A decrease in TSAT levels indicated a reduction in the proportion of transferrin-bound iron. Transferrin is a class of glycoproteins that are easily modified by glycosylation [64]. Against the background of NAFLD, some transferrin may be abnormally glycosylated [65] and cannot bind and transport iron, resulting in an increase in unavailable iron, such as ferrous iron (Fe2+). In the redox cycle, Fe2+ can react with H2O2 in the Fenton reaction to produce a hydroxyl radical (•OH) that causes oxidative damage [57]. We examined the level of Fe2+ in mice in the NAFLD model and found a significant increase in the level of Fe2+ compared to control mice, which further aggravated oxidative stress, liver damage, and the development of NAFLD.

Growing evidence from recent studies has revealed an association between iron metabolism and NAFLD. Studies showed that serum ferritin levels were positively associated with NAFLD in an Italian cohort of patients [3] and the Korean general population [27], while the present study showed a nonsignificant negative association between serum ferritin and the risk of NAFLD in the general population of the United States. This difference might be a limitation in that the samples were only from the U.S. population collected by the NHANES 2017–2018. In addition, ferritin has been reported to be an acute-phase reactant protein [66, 67] whose expression increases during inflammation and may be limited in the assessment of iron stores. Although sTfR is less affected by inflammation, there is a lack of accuracy and little broad consensus regarding the use of sTfR in evaluating iron status [68]. TSAT, which is less affected by inflammation [69], is a more useful biomarker in assessing total body iron status [70] and has a certain value in providing prognosis and prediction information. It has been reported that low TSAT levels are associated with an increased risk of BMI [23] and increased mortality among patients with chronic heart failure (CHF) who have normal serum iron levels [71]. In addition, Campodonico et al. suggested that a TSAT level < 20% identifies heart failure patients with the poorest survival rate [72]. Kovesdy et al. found that decreased TSAT correlated with ID has become the most common nutritional disorder, frequently occurring in chronic inflammatory diseases and affecting a subset of patients with CHF, chronic kidney disease (CKD), and inflammatory bowel disease (IBD) [73]. NAFLD is accompanied by chronic inflammation when it progresses to the more severe form of NASH [7]. In our study, decreased TSAT levels were significantly associated with a high risk of NAFLD.

A significant association between TSAT and NAFLD was found among participants who aged < 60 years, female, Mexican American, without hypertension, without DM, and who had a BMI ≥ 30 kg/m2 (P for trend < 0.05) shown by subgroup analysis. It should be noted that statistical significance of the association in the Mexican American, possibly due to the overall younger age of participants in the Mexican population as compared with other races. Additionally, the association of TSAT and NAFLD differed by age, sex, self-reported race, and BMI (P for interaction < 0.05).

In an animal NAFLD model [11], rats fed a high-fructose diet showed decreased serum levels of TSAT, SI, and SF, showing systemic iron deficiency. The authors also observed iron dysregulation in the livers of these mice [11]. The liver is essential for iron sensing, transport, and regulation and orchestrates systemic iron balance by secreting hepcidin [18, 74]. In the circumstance of decreased iron status, the liver can sense this change and reduce the expression of hepcidin, thereby increasing iron bioavailability [17, 75]. It was reported that a reduced liver hepcidin mRNA level was associated with low TSAT levels [60]. Consistently, in our study, hepatic hepcidin mRNA levels were decreased in mice and patients with NAFLD and correlated with decreased serum TSAT levels in mice. In NAFLD, possibly due to a decrease in TSAT, the liver reduces hepcidin mRNA expression in response to this change, reducing ferroportin degradation, allowing more iron to be transported into the blood, and increasing iron bioavailability [17, 75]. Therefore, we concluded that decreased serum TSAT levels and reduced hepatic hepcidin expression indicate decreased systemic iron status and hepatic iron dysregulation in NAFLD, which might provide additional information on biomarkers for the diagnosis of NAFLD.

The current cross-sectional study enrolled a large number of adult participants, comprehensively analyzed the relationship between iron metabolism indicators and NAFLD and showed a significant negative correlation between serum TSAT and NAFLD, which may provide some information for the diagnosis of NAFLD and guidance for appropriate iron supplementation for patients with NAFLD. Nevertheless, there still exist some limitations that cannot be ignored. First, NAFLD status was defined using the U.S. fatty liver index rather than biopsy, which might contribute to inclusion bias among NAFLD patients. Second, a cross-sectional study design does not allow the establishment of causality. Thus, longitudinal studies are needed to determine whether low TSAT levels could be a good biomarker of NAFLD or if it is just a consequence of liver function alteration. Additionally, unmeasured confounders associated with NAFLD might have affected the results despite adjusting for covariates in the regression models. In the future, the impact of iron metabolism on NAFLD and its specific manifestations in the general population requires further prospective studies.

Conclusion

Decreased serum TSAT levels and hepatic hepcidin expression were observed in both patients and mice with NAFLD. Decreased serum TSAT levels are associated with a higher risk of NAFLD in the U.S. general population. The association was statistically significant among participants who were younger than 60 years, woman, and Mexican American, and those who had a BMI ≥ 30 kg/m2.

Availability of data and materials

The data used in this study are available upon reasonable request to the corresponding author, as well as in the NHANES (https://www.cdc.gov/nchs/nhanes/).

Abbreviations

NAFLD:

Non-alcoholic fatty liver disease

NASH:

Non-alcoholic steatohepatitis

NHANES:

National Health and Nutrition Examination Surveys

FLI:

Fatty Liver Index

BMI:

Body Mass Index

SI:

Serum iron

TSAT:

Transferrin saturation

sTfR:

Soluble transferrin receptor

SF:

Serum ferritin

Family PIR:

Family poverty income ratio

DM:

Diabetes mellitus

PA:

Physical activity

Hb:

Hemoglobin

hs -CRP:

High-sensitivity C-reactive protein

HbA1c:

Glycosylated hemoglobin

ALT:

Alanine aminotransferase

AST:

Aspartate aminotransferase

GGT:

Gamma-glutamyl transpeptidase

TC:

Total cholesterol

TG:

Triglycerides

HDL-C:

High density lipoprotein-cholesterol

BUN:

Blood urea nitrogen

UA:

Uric acid

Scr:

Serum creatinine

eGFR:

Estimated glomerular filtration rate

uACR:

Urinary albumin creatinine ratio

MDA:

Malonic dialdehyde

References:

  1. Forlano R, et al. Non-alcoholic fatty liver disease and vascular disease. Curr Vasc Pharmacol. 2021;19(3):269–79.

    Article  CAS  Google Scholar 

  2. Castera L. Diagnosis of non-alcoholic fatty liver disease/non-alcoholic steatohepatitis: non-invasive tests are enough. Liver Int. 2018;38(Suppl 1):67–70.

    Article  Google Scholar 

  3. Corradini E, et al. Ceruloplasmin gene variants are associated with hyperferritinemia and increased liver iron in patients with NAFLD. J Hepatol. 2021;75(3):506–13.

    Article  CAS  Google Scholar 

  4. Byrne CD, Targher G. NAFLD: a multisystem disease. J Hepatol. 2015;62(1 Suppl):S47-64.

    Article  Google Scholar 

  5. Li Y, et al. Amlodipine, an anti-hypertensive drug, alleviates non-alcoholic fatty liver disease by modulating gut microbiota. Br J Pharmacol. 2022;179(9):2054–77.

    Article  CAS  Google Scholar 

  6. Qian M, et al. Coordinated changes of gut microbiome and lipidome differentiates nonalcoholic steatohepatitis (NASH) from isolated steatosis. Liver Int. 2020;40(3):622–37.

    Article  CAS  Google Scholar 

  7. Wang L, et al. YIPF6 controls sorting of FGF21 into COPII vesicles and promotes obesity. Proc Natl Acad Sci USA. 2019;116(30):15184–93.

    Article  CAS  Google Scholar 

  8. Li J, et al. Ferroptosis: past, present and future. Cell Death Dis. 2020;11(2):88.

    Article  Google Scholar 

  9. Wu J, et al. Ferroptosis in liver disease: new insights into disease mechanisms. Cell Death Discov. 2021;7(1):276.

    Article  CAS  Google Scholar 

  10. Mayneris-Perxachs J, et al. Iron status influences non-alcoholic fatty liver disease in obesity through the gut microbiome. Microbiome. 2021;9(1):104.

    Article  CAS  Google Scholar 

  11. Wang C, et al. A high-fructose diet in rats induces systemic iron deficiency and hepatic iron overload by an inflammation mechanism. J Food Biochem. 2021;45(1):e13578.

    Article  CAS  Google Scholar 

  12. Galaris D, Barbouti A, Pantopoulos K. Iron homeostasis and oxidative stress: an intimate relationship. Biochim Biophys Acta Mol Cell Res. 2019;1866(12):118535.

    Article  CAS  Google Scholar 

  13. Puig S, et al. The elemental role of iron in DNA synthesis and repair. Metallomics. 2017;9(11):1483–500.

    Article  Google Scholar 

  14. Morales M, Xue X. Targeting iron metabolism in cancer therapy. Theranostics. 2021;11(17):8412–29.

    Article  CAS  Google Scholar 

  15. Coad J, Pedley K. Iron deficiency and iron deficiency anemia in women. Scand J Clin Lab Invest Suppl. 2014;244:82–9 (discussion 89).

    Article  Google Scholar 

  16. Moustarah F, Mohiuddin SS. Dietary iron. In: StatPearls. Treasure Island (FL): StatPearls Publishing Copyright©2022, StatPearls Publishing LLC; 2022.

  17. Ganz T, Nemeth E. Hepcidin and iron homeostasis. Biochim Biophys Acta. 2012;1823(9):1434–43.

    Article  CAS  Google Scholar 

  18. Wang CY, Babitt JL. Liver iron sensing and body iron homeostasis. Blood. 2019;133(1):18–29.

    Article  CAS  Google Scholar 

  19. Zhang Z, et al. Ferroportin1 in hepatocytes and macrophages is required for the efficient mobilization of body iron stores in mice. Hepatology. 2012;56(3):961–71.

    Article  CAS  Google Scholar 

  20. Mizutani K, Toyoda M, Mikami B. X-ray structures of transferrins and related proteins. Biochim Biophys Acta. 2012;1820(3):203–11.

    Article  CAS  Google Scholar 

  21. Gkouvatsos K, Papanikolaou G, Pantopoulos K. Regulation of iron transport and the role of transferrin. Biochim Biophys Acta. 2012;1820(3):188–202.

    Article  CAS  Google Scholar 

  22. Silva AMN, Rangel M. The (bio)chemistry of non-transferrin-bound iron. Molecules. 2022;27(6):1784.

    Article  CAS  Google Scholar 

  23. Pilar Vaquero M, et al. Diabesity negatively affects transferrin saturation and iron status. The DICARIVA study. Diabetes Res Clin Pract. 2021;172:108653.

    Article  CAS  Google Scholar 

  24. Knovich MA, et al. Ferritin for the clinician. Blood Rev. 2009;23(3):95–104.

    Article  CAS  Google Scholar 

  25. Elsayed ME, Sharif MU, Stack AG. Transferrin saturation: a body iron biomarker. Adv Clin Chem. 2016;75:71–97.

    Article  CAS  Google Scholar 

  26. Beguin Y. Soluble transferrin receptor for the evaluation of erythropoiesis and iron status. Clin Chim Acta. 2003;329(1–2):9–22.

    Article  CAS  Google Scholar 

  27. Jung JY, et al. Serum ferritin level is associated with liver steatosis and fibrosis in Korean general population. Hepatol Int. 2019;13(2):222–33.

    Article  Google Scholar 

  28. Yang HH, et al. Serum iron and risk of nonalcoholic fatty liver disease and advanced hepatic fibrosis in US adults. Sci Rep. 2021;11(1):10387.

    Article  CAS  Google Scholar 

  29. Ezzati TM, et al. Sample design: third National Health and Nutrition Examination Survey. Vital Health Stat. 1992;2(113):1–35.

    Google Scholar 

  30. Ruhl CE, Everhart JE. Fatty liver indices in the multiethnic United States National Health and Nutrition Examination Survey. Aliment Pharmacol Ther. 2015;41(1):65–76.

    Article  CAS  Google Scholar 

  31. Bedogni G, et al. The Fatty Liver Index: a simple and accurate predictor of hepatic steatosis in the general population. BMC Gastroenterol. 2006;6:33.

    Article  Google Scholar 

  32. Kobyliak N, et al. A multi-strain probiotic reduces the fatty liver index, cytokines and aminotransferase levels in NAFLD patients: evidence from a randomized clinical trial. J Gastrointestin Liver Dis. 2018;27(1):41–9.

    Article  Google Scholar 

  33. Chen YJ, et al. Serum iron and risk of diabetic retinopathy. Nutrients. 2020;12(8):2297.

    Article  CAS  Google Scholar 

  34. Hao H, et al. The association between METS-IR and serum ferritin level in United States female: a cross-sectional study based on NHANES. Front Med (Lausanne). 2022;9:925344.

    Article  Google Scholar 

  35. Zhu S, et al. Increased serum soluble transferrin receptor levels were associated with high prevalence of cardiovascular diseases: insights from the National Health and Nutrition Examination Survey 2017–2018. Front Cell Dev Biol. 2022;10:874846.

    Article  Google Scholar 

  36. Åsberg A, et al. Diagnosing empty iron stores in women: unbound iron binding capacity (UIBC) versus soluble transferrin receptor (sTFR). Scand J Clin Lab Invest. 2021;81(2):137–41.

    Article  Google Scholar 

  37. Heeringa SG, et al. Attributable fraction estimation from complex sample survey data. Ann Epidemiol. 2015;25(3):174–8.

    Article  Google Scholar 

  38. Guidelines for the treatment of animals in behavioural research and teaching. Anim Behav 2000; 59(1): 253–7.

  39. Faruqi A, Mukkamalla SKR. Iron binding capacity. In: StatPearls. Treasure Island (FL): StatPearls Publishing Copyright©2022, StatPearls Publishing LLC; 2022.

  40. Xie LH, et al. Mechanisms of an increased level of serum iron in gamma-irradiated mice. Radiat Environ Biophys. 2016;55(1):81–8.

    Article  CAS  Google Scholar 

  41. Zhang R, et al. Aggregation-induced emission probe for specific turn-on quantification of soluble transferrin receptor: an important disease marker for iron deficiency Anemia and kidney diseases. Anal Chem. 2018;90(2):1154–60.

    Article  CAS  Google Scholar 

  42. Mao GX, et al. Antiaging effect of pine pollen in human diploid fibroblasts and in a mouse model induced by D-galactose. Oxid Med Cell Longev. 2012;2012:750963.

    Article  Google Scholar 

  43. Wang L, et al. Intestinal REG3 lectins protect against alcoholic steatohepatitis by reducing mucosa-associated microbiota and preventing bacterial translocation. Cell Host Microbe. 2016;19(2):227–39.

    Article  CAS  Google Scholar 

  44. Mazagova M, et al. Commensal microbiota is hepatoprotective and prevents liver fibrosis in mice. FASEB J. 2015;29(3):1043–55.

    Article  CAS  Google Scholar 

  45. Wang L, et al. Nod2 deficiency protects mice from cholestatic liver disease by increasing renal excretion of bile acids. J Hepatol. 2014;60(6):1259–67.

    Article  CAS  Google Scholar 

  46. Kleiner DE, et al. Design and validation of a histological scoring system for nonalcoholic fatty liver disease. Hepatology. 2005;41(6):1313–21.

    Article  Google Scholar 

  47. Qian M, et al. Aryl hydrocarbon receptor deficiency in intestinal epithelial cells aggravates alcohol-related liver disease. Cell Mol Gastroenterol Hepatol. 2022;13(1):233–56.

    Article  CAS  Google Scholar 

  48. Arendt BM, et al. Altered hepatic gene expression in nonalcoholic fatty liver disease is associated with lower hepatic n-3 and n-6 polyunsaturated fatty acids. Hepatology. 2015;61(5):1565–78.

    Article  CAS  Google Scholar 

  49. Suppli MP, et al. Hepatic transcriptome signatures in patients with varying degrees of nonalcoholic fatty liver disease compared with healthy normal-weight individuals. Am J Physiol Gastrointest Liver Physiol. 2019;316(4):G462–72.

    Article  CAS  Google Scholar 

  50. Yao K, et al. Transcriptomic profiling of a multiethnic pediatric NAFLD cohort reveals genes and pathways associated with disease. Hepatol Commun. 2022;6(7):1598–610.

    Article  CAS  Google Scholar 

  51. Feng Q, et al. Increased apolipoprotein A5 expression in human and rat non-alcoholic fatty livers. Pathology. 2015;47(4):341–8.

    Article  CAS  Google Scholar 

  52. Card RT, Brown GM, Valberg LS. Serum iron and iron-binding capacity in normal subjects. Can Med Assoc J. 1964;90(10):618–22.

    CAS  Google Scholar 

  53. Garcia-Casal MN, et al. Serum or plasma ferritin concentration as an index of iron deficiency and overload. Cochrane Datab Syst Rev. 2021;5(5):Cd011817.

    Google Scholar 

  54. WHO [Internet]. Iron deficiency anaemia: assessment, prevention and control. In: A guide for programme managers. 2021.

  55. Im YR, et al. A systematic review of animal models of NAFLD finds high-fat, high-fructose diets most closely resemble human NAFLD. Hepatology. 2021;74(4):1884–901.

    Article  CAS  Google Scholar 

  56. Perakakis N, Stefanakis K, Mantzoros CS. The role of omics in the pathophysiology, diagnosis and treatment of non-alcoholic fatty liver disease. Metabolism. 2020;111:154320.

    Article  CAS  Google Scholar 

  57. Kajarabille N, Latunde-Dada GO. Programmed cell-death by ferroptosis: antioxidants as mitigators. Int J Mol Sci. 2019;20(19):4968.

    Article  CAS  Google Scholar 

  58. Sebastiani G, Wilkinson N, Pantopoulos K. Pharmacological targeting of the hepcidin/ferroportin axis. Front Pharmacol. 2016;7:160.

    Article  Google Scholar 

  59. Min HK, et al. Hepcidin, iron indices and bone mineral metabolism in non-dialysis chronic kidney disease. Nephrol Dial Transplant. 2020;35(1):147–54.

    CAS  Google Scholar 

  60. Guo P, et al. Hepcidin, an antimicrobial peptide is downregulated in ceruloplasmin-deficient mice. Peptides. 2009;30(2):262–6.

    Article  CAS  Google Scholar 

  61. González-Domínguez Á, et al. Iron metabolism in obesity and metabolic syndrome. Int J Mol Sci. 2020;21(15):5529.

    Article  Google Scholar 

  62. Datz C, Müller E, Aigner E. Iron overload and non-alcoholic fatty liver disease. Minerva Endocrinol. 2017;42(2):173–83.

    Article  Google Scholar 

  63. Corradini E, Buzzetti E, Pietrangelo A. Genetic iron overload disorders. Mol Aspects Med. 2020;75:100896.

    Article  CAS  Google Scholar 

  64. Ghanbari Z, et al. The effect of glycosylation on the transferrin structure: a molecular dynamic simulation analysis. J Theor Biol. 2016;404:73–81.

    Article  CAS  Google Scholar 

  65. Hülsmeier AJ, et al. Glycosylation site occupancy in health, congenital disorder of glycosylation and fatty liver disease. Sci Rep. 2016;6:33927.

    Article  Google Scholar 

  66. Modares Mousavi SR, et al. Correlation between serum ferritin level and histopathological disease severity in non-alcoholic fatty liver disease. Middle East J Dig Dis. 2018;10(2):90–5.

    Article  Google Scholar 

  67. Kell DB, Pretorius E. Serum ferritin is an important inflammatory disease marker, as it is mainly a leakage product from damaged cells. Metallomics. 2014;6(4):748–73.

    Article  CAS  Google Scholar 

  68. Wish JB. Assessing iron status: beyond serum ferritin and transferrin saturation. Clin J Am Soc Nephrol. 2006;1(Suppl 1):S4-8.

    Article  CAS  Google Scholar 

  69. Szőke D, Panteghini M. Diagnostic value of transferrin. Clin Chim Acta. 2012;413(15–16):1184–9.

    Article  Google Scholar 

  70. Cacoub P, Vandewalle C, Peoc’h K. Using transferrin saturation as a diagnostic criterion for iron deficiency: a systematic review. Crit Rev Clin Lab Sci. 2019;56(8):526–32.

    Article  CAS  Google Scholar 

  71. Moliner P, et al. Clinical correlates and prognostic impact of impaired iron storage versus impaired iron transport in an international cohort of 1821 patients with chronic heart failure. Int J Cardiol. 2017;243:360–6.

    Article  Google Scholar 

  72. Campodonico J, et al. Prognostic role of transferrin saturation in heart failure patients. Eur J Prev Cardiol. 2021;28(15):1639–46.

    Article  Google Scholar 

  73. Cappellini MD, et al. Iron deficiency across chronic inflammatory conditions: international expert opinion on definition, diagnosis, and management. Am J Hematol. 2017;92(10):1068–78.

    Article  CAS  Google Scholar 

  74. Graham RM, et al. Liver iron transport. World J Gastroenterol. 2007;13(35):4725–36.

    Article  CAS  Google Scholar 

  75. Li Y, et al. Regulation of iron homeostasis and related diseases. Mediat Inflamm. 2020;2020:6062094.

    Article  Google Scholar 

Download references

Acknowledgements

We thank all the volunteers who took part in the NHANES.

Funding

This study was supported by the ‘Double First-Class’ University Project (Grants CPU2018GF10 and CPU2018GY31 to L.W.).

Author information

Authors and Affiliations

Authors

Contributions

XZ was responsible for the conception and study design, animal experiments, and manuscript writing; RZ made great contributions to the statistical analysis, chart modifications, and full-text format adjustments. SX contributed to data acquisition and analysis; LW verified the integrity of the data and contributed to the study supervision and review. The version submitted for publication was authorized by all authors. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Lirui Wang.

Ethics declarations

Ethics approval and consent to participate

The NHANES and GEO are public databases. The patients involved in the database received ethical approval. Users can download relevant data for free for research and publication purposes.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher's Note

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

Supplementary Information

Additional file 1: Figure S1

The correlation heatmap of biomarkers of iron metabolism and NAFLD.

Additional file 2: Figure S2

LASSO regression model screening potential predictors of NAFLD. a LASSO regression model cross-validation plot. Draw a vertical line at the optimum with the minimum criterion and 1se of the minimum criterion. When λ = 0.020, we get 13 variables, including TATS, for further analysis. b Coefficient profile plot of predictors. Finally, 13 variables including TSAT, were selected at the optimal lambda. LASSO, least absolute shrinkage and selection operator.

Additional file 3: Figure S3

Restricted cubic spline plots of the association between the TATS and the risk of NAFLD. a Analyses were adjusted for age, BMI, waist, Alt, Ast, HDL-C, HbA1c, protein intake, sex, education level, DM, and alcohol user. b Analyses were adjusted for age, BUN, smoker, waist, Alt, HbA1C, sex, TC, alcohol user, HDL-C, education level, and dietary Vitamin C intake. Solid and dashed lines represent the log-transformed odds ratios and the corresponding 95% confidence intervals. The level of TSAT was log10 transformed.

Additional file 4: Figure S4

The establishment of the NAFLD model in mice. a–c Mice were fed a normal chow or high-fat-fructose diet for 30 weeks (n = 5 for each group). a Body weight of mice after 30 weeks of feeding. b Levels of hepatic TG (Left) and TC (Right). c Levels of hepatic MDA.

Additional file 5: Table S1

The characteristics of the populations between the those with missing values and those without.

Additional file 6: Table S2

Multivariable logistic regression of the associations between Iron metabolism and NAFLD.

Additional file 7: Table S3

Multivariable stepwise regression analysis for the association between TSAT and NAFLD.

Rights and permissions

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

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Zhang, X., Zuo, R., Xiao, S. et al. Association between iron metabolism and non-alcoholic fatty liver disease: results from the National Health and Nutrition Examination Survey (NHANES 2017–2018) and a controlled animal study. Nutr Metab (Lond) 19, 81 (2022). https://doi.org/10.1186/s12986-022-00715-y

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12986-022-00715-y

Keywords