Skip to main content

The association between fatty liver index and onset of diabetes: secondary analysis of a population-based cohort study

Abstract

Background

According to research, the fatty liver index (FLI) is associated with diabetes. However, few studies have been conducted to investigate the relationship between FLI and diabetes risk from various perspectives. This study comprehensively investigated the relationship between FLI and incident diabetes in a large Japanese population.

Methods

This retrospective cohort study included 14,280 participants from Murakami Memorial Hospital in Japan from 2004 to 2015. The independent and dependent variables are FLI and risk of type 2 diabetes mellitus (T2DM), respectively. To examine the link between FLI and incident T2DM, Cox proportional-hazards regression was employed. In addition, we performed a number of sensitivity studies to guarantee the validity of the results. Moreover, we conducted subgroup analyses.

Results

After adjusting covariates, the results showed that FLI was positively associated with the risk of T2DM (HR = 1.019, 95%CI: 1.012, 1.025). Additionally, the sensitivity analysis showed how reliable the outcomes were. And a stronger association between FLI and incident T2DM was observed in the regular exercisers (HR = 1.036, 95%CI: 1.019–1.053, P < 0.0001) and the population without ethanol consumption (HR = 1.028, 95%CI: 1.017–1.039, P < 0.0001). Besides, receiver operating characteristic (ROC) curve analysis showed that FLI was better than waist circumference, triglycerides, body mass index, and gamma-glutamyl transferase in predicting incident T2DM.

Conclusion

FLI is positively associated with incident T2DM.

Peer Review reports

Background

Diabetes is one of the twenty-first century's fastest-rising diseases. According to the International Diabetes Federation, 537 million people were diagnosed with diabetes mellitus worldwide in 2021, with the figure anticipated to rise to 783 million by 2045 [1]. As the most common chronic disease globally, diabetes and its complications significantly impact individuals and their families, national economies, and health systems. Diabetes was expected to cost the world USD 966 billion in direct health spending in 2021 [2]. Sadly, the mortality rate related to type 2 diabetes mellitus (T2DM) complications was high, and it was estimated that more than 6.7 million adults aged 20–79 would die from diabetes in 2021 worldwide [1]. Therefore, it is critical to better understand the diabetes risk factors that can be employed for diabetes prevention and screening.

Numerous studies have revealed that fatty liver is a risk factor for developing diabetes [3,4,5,6]. Nonalcoholic fatty liver disease (NAFLD) is driven by ectopic fat accumulation in the liver. Diagnosis of NAFLD is defined as the presence of hepatic steatosis, ballooning, and lobular inflammation with or without fibrosis [7]. Previous studies have shown that hepatic steatosis could instigate liver injury, inflammation, and insulin resistance (IR) [8, 9]. IR is a determining factor in the pathophysiology of T2DM [10]. The gold standard to diagnose hepatic steatosis is liver biopsy, but it has the limitations such as sampling variability, invasive nature, and high cost. Therefore, numerous non-invasive methods are used for diagnosis, including serum markers or imaging modalities such as ultrasound [11]. Fatty liver index (FLI) is a surrogate marker of hepatic steatosis, which is estimated from a formula that includes waist circumference (WC), body mass index (BMI), gamma-glutamyltransferase (GGT), and triglycerides (TG) [12]. This index could replace ultrasonography to diagnose the extent of hepatic steatosis and is a valid predictor of NAFLD [13]. Several studies have previously observed that FLI was a predictor of diabetes development [5, 13,14,15,16]. A meta-analysis [16] identified a positive association between FLI and the risk of diabetes mellitus incidence. However, because the majority of these studies did not use subgroup analysis, it was impossible to understand differences in the relationship between FLI and diabetes in different populations. Furthermore, the value of baseline FLI in predicting future diabetes occurrence, the optimal FLI threshold, and the sensitivity and specificity of predicting diabetes all require further investigation. In summary, more research was needed on the relationships between FLI and incident diabetes, such as the differences in the relationship between the two in different populations, the specific predictive value of FLI for diabetes, and so on. Therefore, this study aimed to comprehensively assess the relationship between FLI and incident diabetes in a large cohort of Japanese adults.

Methods

Data source and participants

The 'DATADRYAD' database, which was accessible at www.Datadryad.org and offered free raw data downloads, was used to collect the data. In this study, we used the Dryad data package based on the Dryad Terms of Service [17] (Dryad data package: Takuro Okamura et al. (2018) Data from: Ectopic fat obesity presents the greatest risk for incident type 2 diabetes: a population-based longitudinal study. Dryad Digital Repository. https://doi.org/10.1038/s41366-018-0076-3). Variables contained in the database file were as follows: sex, age, BMI, diastolic blood pressure (DBP), WC, systolic blood pressure (SBP), glycated hemoglobin, total cholesterol (TC), fasting plasma glucose (FPG), high-density lipoprotein cholesterol (HDL-C), GGT, ethanol consumption, alanine aminotransferase (ALT), TG, aspartate aminotransferase (AST), smoking status, fatty liver, exercise, follow-up days and censor of T2DM during follow up. The original study's author waived all copyright and associated privileges for this data. As a result, we were able to use these data without violating the authors' rights for our secondary analysis. This study was performed according to the Declaration of Helsinki. The Murakami Memorial Hospital Ethics Committee approved this study, and all subjects provided informed permission [17].

The original data were gathered from the database of Murakami Memorial Hospital in Japan's NAGALA (NAfld in the Gifu Area, Longitudinal Analysis) [17]. In order to minimize selection bias, the participants were collected consecutively from Murakami Memorial Hospital in Japan. Their identity information was encoded into an untraceable code to ensure the participants’ privacy. Inclusion criteria were: participants who took the physical test between 2004 and 2015, and had completed at least two physical examinations. The exclusion criteria were as follows in the original research: (1) participants diagnosed with type 2 diabetes (n = 323) or with fasting plasma glucose was over 6.1 mmol/L at baseline (n = 808), (2) participants with known liver disease, such as hepatitis B or C virus (n = 416), (3) anyone who took any medication at baseline examination (n = 2321), (4) those who drank heavily (ethanol intake greater than 60 g/day for men and 40 g/day for women) (n = 739), (5) participants with a missed value of covariates, including abdominal ultrasonography, exercise, alcohol intake or height (n = 863) [17]. And then, FLI was calculated as mentioned by Bedogni et al. [12]: FLI = 100/(1 + e−z); z = 0.953 × ln[TG (mg/dL)] + 0.139 × BMI (kg/m2) + 0.718 × ln[GGT (U/L)] + 0.053 × WC (cm)-15.745. For this further study, we excluded men with ethanol consumption over 30 g/day and women with ethanol consumption over 20 g/day (n = 1184).

Study design and measurement of variables

The design of the study has been recorded elsewhere [17]. As mentioned in the NAGALA study, this study's clinical baseline information was collected through a standardized self-administered questionnaire, including physical activity, alcohol, smoking habits, and medical history. Ethanol consumption was evaluated by the mean ethanol intake of participants per week during the prior month. Smoking status was categorized into current smoker, ex-smoker, or non-smoker. To categorize participants into non- or regular exercisers, participants’ recreational and sports activities were investigated. Regular exercisers were defined as those who reported any exercise more than once a week [18]. BMI was computed as follows: BMI (kg/m2) = body weight (kg)/height2 (m2). After a night of fasting, venous blood was collected for hematological indicators testing, including TC, GGT, TG, ALT, HDL-C, AST, glycated hemoglobin, and FPG. Fatty liver was detected by the findings of abdominal ultrasonography conducted by competent technicians [19]. Finally, FLI was the target-independent variable obtained at baseline, calculated as Bedogni et al. described [12]. Besides, the dependent variable was incident T2DM acquired during follow-up. Because it was a retrospective study, it reduced the possibility of observation and selection bias.

Ascertainment of incident T2DM

T2DM was defined as fasting plasma glucose ≥ 7 mmol/L, glycated hemoglobin ≥ 6.5%, or self-reported during the follow-up period. The patients would be reviewed on the day of the diagnosis of T2DM or the day of the last visit, whichever comes first [20]. Participants who were unable to be reached for follow-up would still be included in the study [21].

Statistical analysis

First, we dealt with missing values of covariables. In the present study, only TG had missing values. The number of participants due to missing TG values is only 11 (0.07%), so we used the average of TG to impute the missing values [22].

Second, we stratified the participants by quartiles of FLI. Categorical variables were expressed as percentages of a particular group, while continuous variables with normal and skewed distributions were expressed as means with standard deviations or medians with interquartile ranges. Next, the Chi-square, Kruskal–Wallis test, and one-way ANOVA were used for the four groups' categorical, skewed continuous, and normal continuous variables. Person-years of follow-up were computed from the date of the baseline interview until the occurrence of T2DM or the date of the follow-up interview [23]. The incidence rates were given in terms of cumulative incidence and person-year incidence [24]. The Kaplan–Meier method was used to calculate survival estimates and time-to-event variables. A log-rank test was employed to assess the Kaplan–Meier survival probability among FLI groups [25].

Third, and foremost, the proportion hypothesis was examined. The Cox proportional hazards regression model was then used to estimate the risk ratios (HRs) and 95% confidence intervals (CIs) of incident T2DM after meeting the proportion criteria. The results of unadjusted, minimally adjusted ( SBP, gender, smoking status, DBP, age, ethanol consumption, and exercise), and fully adjusted ( DBP, ALT, smoking status, SBP, exercise, AST, ethanol consumption, age, glycated hemoglobin, TC, gender, HDL-C, fatty liver) analyses were simultaneously shown based on the STROBE statement [21]. When the covariances were included in the model and the hazard ratio changed by 10% or more, we made the necessary adjustments [26].

Additionally, we conducted a series of sensitivity analyses to ensure the robustness of the data analysis [27]. FLI was transformed into a categorical variable, and the P-value for the trend was determined. The goal of the test was to confirm the results of considering FLI as a continuous variable and to examine the probability of non-linearity. We excluded participants with fatty liver or any alcohol consumers in other sensitivity analyses. Besides, in order to guarantee the reliability of the results, we also employed generalized additive models (GAM) to include the continuity covariate as a curve into the equation [28]. Additionally, by computing E-values, we investigated the possibility of unmeasured confounding between FLI and incident T2DM [29].

Fourth, we stratified Cox proportional hazard models to investigate the robustness of the results across different subgroups (SBP, DBP, FPG, exercise, age, fatty liver, ethanol consumption, sex, smoking status). Firstly, we converted the continuous variable SBP (< 140 mmHg, ≥ 140 mmHg), DBP (< 90 mmHg, ≥ 90 mmHg), age (< 30, 30 to 40, 40 to 50, 50 to 60, ≥ 60) [30], ethanol consumption (= 0 g/week, > 0 g/week) to a categorical variable based on the clinical cut point while converted FPG based on binary. Secondly, in addition to the stratification factor itself, we adjusted each stratification for all factors (DBP, sex, ALT, SBP, AST, fatty liver, glycated hemoglobin, HDL-C, ethanol consumption, age, smoking status, TC, and exercise). Lastly, the likelihood ratio test was used to test for interaction in models with and without interaction terms [31, 32].

Finally, the ability of FLI, WC, TG, BMI, and GGT to predict the risk of T2DM was estimated using a receiver operating characteristic (ROC) curve. We also calculated net reclassification improvement (NRI) and integrated discrimination improvement (IDI) to compare the predictive values among different models.

The statistical software programs R (http://www.R-project.org, The R Foundation) and Empower-Stats (http://www.empowerstats.com, X&Y Solutions, Inc., Boston, MA) were used for all these studies. P values that were less than 0.05 (two-sided) were regarded as statistically significant.

Results

Initially, the study recruited 20,944 participants; afterward, 6,664 participants were excluded, and 14,280 persons (47.90% women and 52.10% men) remained for data analysis (Fig. 1). The average age was 43.53 ± 8.89 years old. During a mean follow-up period of 5.38 years, 324 (2.27%) individuals developed T2DM. The mean FLI was 12.61 ± 16.43, and the mean FPG was 5.15 ± 0.41 mmol/L.

Fig. 1
figure 1

Flowchart of study participants

Baseline characteristics of the study participants

We assigned participants into subgroups using FLI quartiles (< 2.19,2.19–5.70,5.70–15.97, ≥ 15.97) and showed baseline characteristics of the total population in Table 1. In the highest FLI group, people generally had higher BMI, age, DBP, SBP, ethanol consumption, WC, AST, GGT, ALT, TG, TC, glycated hemoglobin, FPG, and higher rates of current smokers and ex-smokers. Besides, the group (FLI ≥ 15.97) had a higher proportion of men and fatty liver.

Table 1 The Baseline Characteristics of participants presented by quartiles of FLI

Figure 2 showed the distribution of FLI levels. It had a skewed distribution and was in the range of 0.10 to 99.13. This result also showed that the majority of the FLI was less than 60. Participants were separated into two groups based on whether or not they developed T2DM in the future. Figure S1 depicts the FLI values in the two groups. The results showed that the distribution level of FLI was higher in the T2DM group, whereas it was comparatively low in the T2DM-free group. Male individuals had a greater incidence rate of T2DM than female subjects in age stratification by ten intervals, regardless of age group (Fig. 3). We also observed that the incidence of T2DM increased with age in female participants.

Fig. 2
figure 2

Distribution of FLI. It presented a skewed distribution while being in the range from 0.10 to 99.13

Fig. 3
figure 3

T2DM incidence of age stratification by 10 intervals

The incidence rate of incident T2DM

According to Table 2, 324 people had T2DM overall. All individuals' combined incidence rate was 375.23 per 100,000 person-years. Specifically, the incidence rates of the four FLI groups were 68.39, 122.69, 328.75, and 980.35 per 100,000 person-years, respectively. Participants with high FLI had a greater cumulative incidence (P < 0.001 for trend) than those with low FLI. The cumulative incidence of total incident T2DM and each FLI group was 2.269% (2.025%-2.513%), 0.420% (0.208%-0.632%), 0.728% (0.449%-1.007%), 1.989% (1.531%-2.447%), and 5.938% (5.163%-6.714%), respectively.

Table 2 Incidence rate of incident T2DM presented by quartiles of FLI

Univariate analysis

The results of the univariate analysis were shown in Table 3. By univariate Cox proportional hazard model, we found that exercise and ethanol consumption were not significantly associated with T2DM. We also found that HDL-C was negatively associated with incident T2DM. Univariate analysis illustrated that age, BMI, WC, ALT, SBP, GGT, glycated hemoglobin, TG, FPG, TC, AST, and DBP were positively associated with the risk of T2DM. Meanwhile, men, current and ex-smokers, and participants with fatty liver had a higher risk of developing T2DM.

Table 3 The results of univariate Cox proportional risk regression model for riks factors of T2DM

Figure 4 displayed the Kaplan–Meier survival curves for the probability of T2DM-free survival stratified by FLI groups. The probability of surviving without developing T2DM varied considerably between FLI groups (log-rank test, p < 0.0001). The probability of surviving without T2DM rapidly fell as FLI increased, revealing the top group at greatest risk for developing T2DM.

Fig. 4
figure 4

Kaplan–Meier event-free survival curve. Kaplan–Meier event-free survival curve. Kaplan–Meier analysis of incident T2DM-free survival based on FLI groups (log-rank, P < 0.0001)

The results of the relationship between FLI and incident T2DM

The relationship between FLI and T2DM events was examined using the Cox proportional hazards regression model. In Table 4, we simultaneously displayed the three adjusted models and the non-adjusted models. FLI was found to be positively linked with incident T2DM in the crude model (HR = 1.041, 95% CI:1.038–1.045, P < 0.0001). The minimum adjusted model (adjusted DBP, age, SBP, ethanol consumption, sex, smoking status, exercise) did not show any discernible changes in the outcome (HR = 1.040, 95%CI: 1.035–1.045). A positive association between FLI and incident T2DM was also observed after accounting for the entire model (adjusted ALT, exercise, TC, age, HDL-C, AST, DBP, sex, SBP, alcohol intake, smoking status, glycated hemoglobin, fatty liver) (HR = 1.019, 95%CI: 1.012–1.025, P < 0.0001). The findings indicated that a unit rise in FLI raised the risk of T2DM by 1.9%.

Table 4 Relationship between FLI and the incident T2DM in different Cox proportional risk models

Sensitivity analysis

The FLI was turned into a categorical variable for sensitivity analysis (Quartile). The findings of Model II revealed that the probability of acquiring T2DM rose with FLI quartiles 2, 3, and 4 compared to FLI quartile 1[HR 1.307 (0.685, 2.495), 1.799 (0.979, 3.308), and 2.262 (1.171, 4.368)]. In the entire model, the top quartile experienced a 1.26 times increase in T2DM risk compared to the bottom FLI quartile. And the T2DM risk trend in FLI quartiles was significant (P for trend = 0.00626) (Table 4).

In addition, we included the continuity covariate into the equation as a curve using a GAM. Model III results indicated the consistency of the findings with the fully adjusted model (HR = 1.020, 95%CI: 1.013–1.027, P < 0.00001), demonstrating the robustness of the findings (Table 4). The results of Model III also revealed that the probability of acquiring T2DM rose with FLI quartiles 2, 3, and 4 compared to FLI quartile 1[HR 1.261 (0.650, 2.446), 1.631 (0.854, 3.115), and 2.012 (1.003, 4.035)].

Furthermore, we calculated an E-value to evaluate the sensitivity to unmeasured confounding variables. The value of E was 1.16. The E-value was bigger than the relative risk of unmeasured confounders and incident T2DM, indicating that unmeasured or unknown confounders had minimal influence on the association between FLI and T2DM risk.

Furthermore, other sensitivity analyses excluded patients with fatty liver. After controlling for confounding factors, we found that the relationship between FLI and incident T2DM was not statistically significant in people without fatty liver (HR = 1.003, 95% CI:0.988 to 1.018). (Table 5). For sensitivity analysis, we additionally omitted any alcohol drinkers. After controlling for gender, SBP, ALT, DBP, TC, AST, HDL-C, smoking status, Glycated hemoglobin, exercise, age, and fatty liver, the results showed that FLI was still linked with incident T2DM (HR = 1.028, 95% CI:1.016 to 1.041, P < 0.00001) (Table 5). The findings of the sensitivity analysis revealed that the link between FLI and the risk of T2DM was very strong.

Table 5 Relationship between FLI and incident T2DM by Cox proportional risk models in different sensitivity analyses

The results of subgroup analyses

We used subgroup analysis to take into account other influencing factors, such as SBP, that could affect the findings on the connection between FLI and incident T2DM. We used FPG, exercise, age, ethanol consumption, fatty liver, sex, smoking status, SBP and DBP as the stratification variables to detect the trend of effect sizes in these variables (Table 6). Table 6 showed that exercise and ethanol consumption could modify the relationship between FLI and incident T2DM (P for interaction < 0.05). And a stronger association between FLI and incident T2DM was observed in the regular exercisers (HR = 1.036, 95%CI: 1.019–1.053, P < 0.0001) and the population without ethanol consumption (HR = 1.028, 95%CI: 1.017–1.039, P < 0.0001). In contrast, the weaker association between FLI and incident T2DM was probed in the people who was consuming alcohol (HR = 1.015, 95%CI:1.008, 1.022) and not exercising regularly (HR = 1.016, 95%CI:1.009, 1.023). These findings revealed a robust association between FLI and incident T2DM in the majority of subgroups.

Table 6 Effect size of FLI on T2DM in prespecified and exploratory subgroups

The results of the ROC curve analysis

In addition, we drew a ROC curve to measure the ability of FLI, WC, BMI, TG, and GGT to predict the risk of T2DM (Fig. 5). The areas under the curve (AUC) of each variable were as follows: GGT: 0.709 < TG: 0.733 < BMI: 0.749 < WC: 0.753 < FLI: 0.789 (Table 7). The highest Youden index of GGT, TG, BMI, WC, and FLI was 0.3365, 0.3711, 0.3792, 0.3941, 0.4463, and the corresponding optimal cut-off value was 17.5000, 0.8637, 23.5285, 83.1500, 13.5821, respectively. Besides, we calculated the net reclassification improvement (NRI) and integrated discrimination improvement (IDI) to compare the predictive values among the different models. Compared with WC, the IDI and NRI of the FLI were all increased(Table S1). Combining the AUC, Youden index, IDI, and NRI, we concluded that the predictive ability of FLI to incident T2DM was better than that of other variables (Table 7, Table S1).

Fig. 5
figure 5

The results of ROC curve analysis for measuring the ability of FLI, WC, BMI, TG, and GGT to predict the risk of T2DM. The areas under the curve (AUC) of each variable were as follows: GGT: 0.709 < TG: 0.733 < BMI: 0.749 < WC: 0.753 < FLI: 0.789

Table 7 Predictive value FLI, WC, TG, BMI, and GGT for predicting T2DM

Discussion

Our findings demonstrated that FLI was positively linked with the development of T2DM after controlling for other variables. Moreover, we conducted subgroup analysis to better understand the trend of FLI and the incidence of T2DM in different populations. Our study detected a stronger association in people with regular exercise and the population without ethanol consumption. On the contrary, a weaker relation between FLI and incident T2DM was observed in people who were consuming alcohol and not exercising regularly.

FLI is a simple and valuable clinical biomarker of hepatic steatosis because it only requires simple biochemical (GGT and TG) and anthropometric (WC and BMI) measurements. So FLI had been adopted by the European Association for the Study of the Liver (EASL), the European Association for the Study of Diabetes (EASD), and the European Association for the Study of Obesity (EASO) NAFLD guidelines as the preferred diagnostic tool for larger-scale screening studies [33]. It is accurate in detecting fatty liver. This index could also replace ultrasonography to diagnose the degree of liver steatosis [12]. Several preceding studies had proposed that FLI could predict the risk of diabetes [34,35,36,37,38,39,40], which agreed with our conclusions. Moreover, we constructed a ROC curve to estimate the ability of WC, BMI, TG, FLI, and GGT to predict the risk of T2DM, and we found that FLI’s AUC and Youden index was better than any other component of FLI.

Researchers have clarified the associations between FLI and the risk of diabetes [39, 40]. Several previous studies mentioned that high FLI increases the risk of developing diabetes mellitus [41,42,43]. For instance, FLI > 60 is independently related to the occurrence of T2DM after three years of follow-up in the Spanish PREDAPS study regardless of sex, age, or education level [44]. A recently published meta-analysis[16]with a sample of 70,918(27 studies) participants reported a direct relationship between the highest grade of FLI and an increased diabetes risk, which was not affected by continent, sex, and the quality of study after subgroup analysis. Besides, Wargny et al. [45] found that low FLI value was independently related to prediabetes reversion. FLI strongly predicted the risk of diabetes in a 5-year, prospective, observational study of 389 individuals with prediabetes, independently of traditional risk factors for diabetes mellitus such as sex, age, FPG, or glycated hemoglobin. We obtained a similar result from the Cox proportional hazard regression model, suggesting a significant and robust association between FLI and incident T2DM. It was worth noting that most participants in the present study had an FLI of less than 60. Firstly, the Japanese population might have different distribution of FLI compared to Caucasian populations due to different anthropometric characteristics. In addition, this study excluded participants with FPG > 6.1 mmol/L and heavy drinking (ethanol consumption of more than 30 g per day for men and 20 g per day for women). While this study also found that alcohol consumption and FPG increased with the increase in FLI. In addition, we found that they were also risk factors for diabetes in the present study. However, the PREDAPS study found FLI was independently associated with the incidence of T2DM when it was above 60 [44]. Obviously, the baseline alcohol consumption and FPG levels in the PREDAPS study were higher than those in the present study because their study population was pre-diabetes participants. Incident diabetes was positively associated with FLI, according to our research. The sensitivity analysis, however, revealed that this association is still present among those who do not use alcohol. Despite the majority of research participants having FLI levels below 60, the results served as a guide for treatment interventions to lower the risk of diabetes. This discovery would further expand the results of the PREDAPS study.

Surprisingly, by univariate Cox proportional hazard model, we found that exercise was not significantly associated with T2DM. However, as can be seen in Table 3, although exercise was not significantly associated with T2DM, there was a trend to decrease the risk of T2DM. Moreover, we did a subgroup analysis, and the interaction test was statistically significant. We found that in the population with regular exercise, the relationship between FLI and T2DM was enhanced. The reason might be that in patients with exercise, the level of other diabetes risk factors such as BMI and lipid profile decreased, so the effect of FLI on diabetes was relatively enhanced.

There are some strengths in our study. First, compared with most previous studies of the same kind, our sample size is relatively large. Second, since this study was observational, it was susceptible to potential confounding factors. We made rigorous statistical adjustments to reduce the interference of confounding factors to ensure the reliability of the results. Third, to ensure the validity of the findings, we performed a number of sensitivity analyses. These included categorizing the FLI, using a GAM to incorporate the continuity covariate as a curve into the equation, calculating E-values to investigate the possibility of unmeasured confounding, and re-examining the relationship between FLI and incident T2DM after excluding alcohol consumers and participants with fatty liver. Fourth, the effect modifier factor analysis allowed us to use the data better and draw conclusions steadily in the different subgroups of our study. Fifth, to assess the accuracy of the FLI, TG, WC, GGT, and BMI in predicting the risk of T2DM, we created a ROC curve.

There are some limitations in the current study because of the original study [17]. First, since all subjects are of Japanese descent, this study is not generalisable to all population. More research is required to ascertain the prevalence of diabetes, the relationship between FLI and diabetes, and more in people with various genetic backgrounds. The findings of this study might not be generalizable, given that subjects with a significant drinking habit, viral hepatitis, or any medication usage at baseline were omitted. Second, because this study is based on a secondary analysis of previously published data, it is not possible to make adjustments for variables like insulin concentration, low-density lipoprotein cholesterol, very low-density lipoprotein cholesterol, tumor necrosis factor, and interleukin-6 that were not part of the original dataset. To evaluate the possible impact of unmeasured confounders, we estimated the E-value and observed that unmeasured confounders were unlikely to explain the results. Third, censoring by death during follow-up was not included in the data from the initial study. Possible participant fatalities during follow-up were unavoidable for such a large sample of subjects. We can plan our research and gather information on potential future death censorship. So, using a competitive risk model, we may examine the connection between FLI and diabetes. Fourth, due to the lack of oral glucose tolerance tests, the incidence of diabetes may be underestimated. However, due to practical reasons and logistics, performing oral glucose tolerance tests on all participants is not feasible. Fifth, the study used ultrasonography that could not detect steatohepatitis (SH), which diagnosis could be obtained only by biopsy [12]. However, for obvious ethical reasons, an SH score will never be available in a representative sample of the general population. We can use other non-invasive scores to diagnose NAFLD for further studies, such as NAFLD Fibrosis score, Hepatic Steatosis index, FIB-4, APRI, lipid accumulation product. Sixth, the FLI and other parameters were only assessed at baseline in the current investigation, and FLI changes over time, information on first-line laboratory tests and repetition were not considered. In the future, we can think about structuring our studies or working with other researchers to gather as many data points as possible, such as details on the development of FLI during patient follow-up. Finally, rather than establishing a causal relationship between FLI and diabetes risk, this retrospective observational analysis presented an association inference.As a result, our findings should be interpreted with caution and confirmed by future research.

Conclusion

This study found a positive connection between FLI and incident T2DM in the Japanese population. As a result, aberrant FLI assists in identifying the Japanese population at high risk of T2DM, which would allow doctors to plan and implement suitable care methods ahead of time. And it would allow people to change their lifestyles in advance to lower the occurrence of T2DM by reminding them of the increased risk of diabetes during follow-up.

Availability of data and materials

Data can be downloaded from the ‘DATADRYAD’ database (www.Datadryad.org). All data generated or analysed during this study are included in the supplementary information files.

Abbreviations

ALT:

Alanine aminotransferase

AUC:

Area under the curve

AST:

Aspartate aminotransferase

DBP:

Diastolic blood pressure

BMI:

Body mass index

FPG:

Fasting plasma glucose

FLI:

Fatty liver index

GAM:

Generalized additive models

GGT:

Gammaglutamyltransferase

HDL-C:

High-density lipoprotein cholesterol

IDI:

Integrated discrimination improvement

IR:

Insulin resistance

NRI:

Net reclassification improvement

ROC:

Receiver operating characteristic

SBP:

Systolic blood pressure

SH:

Steatohepatitis

T2DM:

Type 2 diabetes mellitus

TC:

Total cholesterol

TG:

Triglyceride

WC:

Waist circumference

References

  1. Magliano DJ, Boyko EJDA. IDF DIABETES ATLAS. Brussels: International Diabetes Federation; 2021.

    Google Scholar 

  2. Sun H, Saeedi P, Karuranga S, Pinkepank M, Ogurtsova K, Duncan BB, et al. IDF Diabetes Atlas: Global, regional and country-level diabetes prevalence estimates for 2021 and projections for 2045. Diabetes Res Clin Pr. 2022;183:109119.

    Article  Google Scholar 

  3. Gastaldelli A, Kozakova M, Hojlund K, Flyvbjerg A, Favuzzi A, Mitrakou A, et al. Fatty liver is associated with insulin resistance, risk of coronary heart disease, and early atherosclerosis in a large European population. Hepatology. 2009;49:1537–44.

    Article  CAS  PubMed  Google Scholar 

  4. Okamura T, Hashimoto Y, Hamaguchi M, Obora A, Kojima T, Fukui M. Effect of alcohol consumption and the presence of fatty liver on the risk for incident type 2 diabetes: a population-based longitudinal study. BMJ Open Diabetes Res Care. 2020;8:e1629.

    Article  Google Scholar 

  5. Lee SM, Kwak SH, Koo JN, Oh IH, Kwon JE, Kim BJ, et al. Non-alcoholic fatty liver disease in the first trimester and subsequent development of gestational diabetes mellitus. Diabetologia. 2019;62:238–48.

    Article  CAS  PubMed  Google Scholar 

  6. Cho HJ, Hwang S, Park JI, Yang MJ, Hwang JC, Yoo BM, et al. Improvement of nonalcoholic fatty liver disease reduces the risk of type 2 diabetes mellitus. Gut Liver. 2019;13:440–9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  7. Pouwels S, Sakran N, Graham Y, Leal A, Pintar T, Yang W, et al. Non-alcoholic fatty liver disease (NAFLD): a review of pathophysiology, clinical management and effects of weight loss. Bmc Endocr Disord. 2022;22:63.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  8. Fujii H, Kawada N. The role of insulin resistance and diabetes in nonalcoholic fatty liver disease. Int J Mol Sci. 2020;21:3863.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  9. Tilg H, Adolph TE, Dudek M, Knolle P. Non-alcoholic fatty liver disease: the interplay between metabolism, microbes and immunity. Nat Metab. 2021;3:1596–607.

    Article  CAS  PubMed  Google Scholar 

  10. Sampath KA, Maiya AG, Shastry BA, Vaishali K, Ravishankar N, Hazari A, et al. Exercise and insulin resistance in type 2 diabetes mellitus: A systematic review and meta-analysis. Ann Phys Rehabil Med. 2019;62:98–103.

    Article  Google Scholar 

  11. Papatheodoridi M, Cholongitas E. Diagnosis of Non-alcoholic Fatty Liver Disease (NAFLD): Current concepts. Curr Pharm Des. 2018;24:4574–86.

    Article  CAS  PubMed  Google Scholar 

  12. Bedogni G, Bellentani S, Miglioli L, Masutti F, Passalacqua M, Castiglione A, et al. The fatty liver index: a simple and accurate predictor of hepatic steatosis in the general population. Bmc Gastroenterol. 2006;6:33.

    Article  PubMed  PubMed Central  Google Scholar 

  13. Koralegedara IS, Warnasekara JN, Rathnayake A, Dayaratne KG, Agampodi SB. Fatty Liver Index is a valid predictor of non-alcoholic fatty liver disease (NAFLD) in pregnancy. BMJ Open Gastroenterol. 2022;9:e000913.

    Article  PubMed  PubMed Central  Google Scholar 

  14. Giorda C, Forlani G, Manti R, Mazzella N, De Cosmo S, Rossi MC ,et al. Occurrence over time and regression of nonalcoholic fatty liver disease in type 2 diabetes. Diabetes Metab Res Rev. 2017;33. https://doi.org/10.1002/dmrr.2878.

  15. Busquets-Cortes C, Bennasar-Veny M, Lopez-Gonzalez AA, Fresneda S, Abbate M, Yanez AM. Utility of Fatty Liver Index to predict reversion to normoglycemia in people with prediabetes. PLoS ONE. 2021;16:e249221.

    Article  Google Scholar 

  16. Movahedian M, Rahmani J, HashemiNazari SS, Mohamadi S, Naik G, Hekmatdoost A. Fatty liver index and risk of diabetes incidence: A systematic review and dose-response meta-analysis of cohort studies. Prim Care Diabetes. 2020;14:577–83.

    Article  PubMed  Google Scholar 

  17. Okamura T, Hashimoto Y, Hamaguchi M, Obora A, Kojima T, Fukui M. Ectopic fat obesity presents the greatest risk for incident type 2 diabetes: a population-based longitudinal study. Int J Obesity. 2019;43:139–48.

    Article  Google Scholar 

  18. Ryu S, Chang Y, Kim DI, Kim WS, Suh BS. gamma-Glutamyltransferase as a predictor of chronic kidney disease in nonhypertensive and nondiabetic Korean men. Clin Chem. 2007;53:71–7.

    Article  CAS  PubMed  Google Scholar 

  19. Hamaguchi M, Kojima T, Itoh Y, Harano Y, Fujii K, Nakajima T, et al. The severity of ultrasonographic findings in nonalcoholic fatty liver disease reflects the metabolic syndrome and visceral fat accumulation. Am J Gastroenterol. 2007;102:2708–15.

    Article  PubMed  Google Scholar 

  20. Harreiter J. Roden M [Diabetes mellitus-Definition, classification, diagnosis, screening and prevention (Update 2019)]. Wien Klin Wochenschr. 2019;131:6–15.

    Article  PubMed  Google Scholar 

  21. Yokoyama M, Watanabe T, Otaki Y, Takahashi H, Arimoto T, Shishido T, et al. Association of the aspartate aminotransferase to alanine aminotransferase ratio with bnp level and cardiovascular mortality in the general population: the yamagata study 10-year follow-up. Dis Markers. 2016;2016:4857917.

    Article  PubMed  PubMed Central  Google Scholar 

  22. Erviti J, Alonso A, Oliva B, Gorricho J, Lopez A, Timoner J, et al. Oral bisphosphonates are associated with increased risk of subtrochanteric and diaphyseal fractures in elderly women: a nested case-control study. Bmj Open. 2013;3:e002091.

    Article  PubMed  PubMed Central  Google Scholar 

  23. Guo C, Qin P, Li Q, Zhang D, Tian G, Liu D, et al. Association between mean arterial pressure and risk of type 2 diabetes mellitus: The Rural Chinese Cohort Study. Prim Care Diabetes. 2020;14:448–54.

    Article  PubMed  Google Scholar 

  24. Qin H, Chen Z, Zhang Y, Wang L, Ouyang P, Cheng L, et al. Triglyceride to high-density lipoprotein cholesterol ratio is associated with incident diabetes in men: A retrospective study of Chinese individuals. J Diabetes Invest. 2020;11:192–8.

    Article  CAS  Google Scholar 

  25. Robson ME, Tung N, Conte P, Im SA, Senkus E, Xu B, et al. OlympiAD final overall survival and tolerability results: Olaparib versus chemotherapy treatment of physician’s choice in patients with a germline BRCA mutation and HER2-negative metastatic breast cancer. Ann Oncol. 2019;30:558–66.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  26. Vandenbroucke JP, von Elm E, Altman DG, Gotzsche PC, Mulrow CD, Pocock SJ, et al. Strengthening the Reporting of Observational Studies in Epidemiology (STROBE): explanation and elaboration. Int J Surg. 2014;12:1500–24.

    Article  PubMed  Google Scholar 

  27. Sun D, Li W, Zhang H, Li Y, Zhang Q. Inverted U-shaped relationship between body mass index and multivessel lesions in Chinese patients with myocardial infarction: a cross-sectional study. J Int Med Res. 2020;48:1220732372.

    Article  Google Scholar 

  28. Zhu F, Chen C, Zhang Y, Chen S, Huang X, Li J, et al. Elevated blood mercury level has a non-linear association with infertility in U.S. women: Data from the NHANES 2013-2016. Reprod Toxicol. 2020;91:53–8.

    Article  CAS  PubMed  Google Scholar 

  29. Haneuse S, VanderWeele TJ, Arterburn D. Using the E-Value to assess the potential effect of unmeasured confounding in observational studies. JAMA. 2019;321:602–3.

    Article  PubMed  Google Scholar 

  30. Zou Y, Sheng G, Yu M, Xie G. The association between triglycerides and ectopic fat obesity: An inverted U-shaped curve. PLoS ONE. 2020;15:e243068.

    Article  Google Scholar 

  31. Mullee A, Romaguera D, Pearson-Stuttard J, Viallon V, Stepien M, Freisling H, et al. Association between soft drink consumption and mortality in 10 European countries. Jama Intern Med. 2019;179:1479–90.

    Article  PubMed  Google Scholar 

  32. Keidel D, Anto JM, Basagana X, Bono R, Burte E, Carsin AE, et al. The Role of Socioeconomic Status in the Association of Lung Function and Air Pollution-A Pooled Analysis of Three Adult ESCAPE Cohorts. Int J Environ Res Public Health. 2019;16:1901.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  33. Easl EAFT, Easd EAFT, Easo EAFT. EASL-EASD-EASO Clinical Practice Guidelines for the management of non-alcoholic fatty liver disease. J Hepatol. 2016;64:1388–402.

    Article  Google Scholar 

  34. Nishi T, Babazono A, Maeda T, Imatoh T, Une H. Evaluation of the fatty liver index as a predictor for the development of diabetes among insurance beneficiaries with prediabetes. J Diabetes Investig. 2015;6:309–16.

    Article  PubMed  Google Scholar 

  35. Balkau B, Lange C, Vol S, Fumeron F, Bonnet F. Nine-year incident diabetes is predicted by fatty liver indices: the French D.E.S.I.R. study. Bmc Gastroenterol. 2010;10:56.

    Article  PubMed  PubMed Central  Google Scholar 

  36. Ruckert IM, Heier M, Rathmann W, Baumeister SE, Doring A, Meisinger C. Association between markers of fatty liver disease and impaired glucose regulation in men and women from the general population: the KORA-F4-study. PLoS ONE. 2011;6:e22932.

    Article  PubMed  PubMed Central  Google Scholar 

  37. Yadav D, Choi E, Ahn SV, Koh SB, Sung KC, Kim JY, et al. Fatty liver index as a simple predictor of incident diabetes from the KoGES-ARIRANG study. Medicine (Baltimore). 2016;95:e4447.

    Article  CAS  PubMed  Google Scholar 

  38. Jung CH, Lee WJ, Hwang JY, Yu JH, Shin MS, Lee MJ, et al. Assessment of the fatty liver index as an indicator of hepatic steatosis for predicting incident diabetes independently of insulin resistance in a Korean population. Diabet Med. 2013;30:428–35.

    Article  CAS  PubMed  Google Scholar 

  39. Olubamwo OO, Virtanen JK, Pihlajamaki J, Tuomainen TP. Association of fatty liver index with risk of incident type 2 diabetes by metabolic syndrome status in an Eastern Finland male cohort: a prospective study. BMJ Open. 2019;9:e26949.

    Article  Google Scholar 

  40. Hirata A, Sugiyama D, Kuwabara K, Hirata T, Tsutatani H, Funamoto M, et al. Fatty liver index predicts incident diabetes in a Japanese general population with and without impaired fasting glucose. Hepatol Res. 2018;48:708–16.

    Article  CAS  PubMed  Google Scholar 

  41. van den Berg EH, Flores-Guerrero JL, Gruppen EG, de Borst MH, Wolak-Dinsmore J, Connelly MA, et al. Non-Alcoholic Fatty Liver Disease and Risk of Incident Type 2 Diabetes: Role of Circulating Branched-Chain Amino Acids. Nutrients. 2019;11:705.

  42. Jager S, Jacobs S, Kroger J, Stefan N, Fritsche A, Weikert C, et al. Association between the fatty liver index and risk of type 2 diabetes in the EPIC-Potsdam study. PLoS ONE. 2015;10:e124749.

    Article  Google Scholar 

  43. Bennet L, Groop L, Franks PW. Country of birth modifies the association of fatty liver index with insulin action in Middle Eastern immigrants to Sweden. Diabetes Res Clin Pr. 2015;110:66–74.

    Article  CAS  Google Scholar 

  44. Franch-Nadal J, Caballeria L, Mata-Cases M, Mauricio D, Giraldez-Garcia C, Mancera J, et al. Fatty liver index is a predictor of incident diabetes in patients with prediabetes: The PREDAPS study. PLoS ONE. 2018;13:e198327.

    Article  Google Scholar 

  45. Wargny M, Smati S, Pichelin M, Bigot-Corbel E, Authier C, Dierry V, et al. Fatty liver index is a strong predictor of changes in glycemic status in people with prediabetes: The IT-DIAB study. PLoS ONE. 2019;14:e221524.

    Article  Google Scholar 

Download references

Acknowledgements

Not applicable

Funding

This study was supported in part by the Discipline Construction Ability Enhancement Project of the Shenzhen Health Commission (SZXJ2017031).

Author information

Authors and Affiliations

Authors

Contributions

Yanqiong ZHU and Haofei HU contributed to the study concept and design, researched and interpreted the data and drafted the manuscript. Yumei WU and Yinhua RAO researched data and reviewed the manuscript. Qixiang LI, Xuehui DUAN and Guopeng YAO oversaw the progress of the project, contributed to the discussion and reviewed the manuscript. As the work's guarantors, Hekun YIN, Ziyi LUO, and Haofei HU had complete access to all the study's data and are accountable for the data's integrity and the analysis' accuracy. The final manuscript was read and approved by all writers.

Corresponding authors

Correspondence to Hekun Yin or Ziyi Luo.

Ethics declarations

Ethics approval and consent to participate

This study was performed according to the Declaration of Helsinki. And in the previously published article [17], Okamura, et al. had clearly stated that: the study was approved by the ethics committee of Murakami Memorial Hospital. Informed consent was obtained from all participants.

Consent for publication

Not applicable.

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.

Data.

Additional file 2:

Table S1. Model performance of FLI compared with WC. Figure S1. Data visualization of FLI of all participants from the T2DM and non-T2DM groups.

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

Zhu, Y., Hu, H., Wu, Y. et al. The association between fatty liver index and onset of diabetes: secondary analysis of a population-based cohort study. BMC Public Health 23, 679 (2023). https://doi.org/10.1186/s12889-023-15442-z

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12889-023-15442-z

Keywords