BioMed Research International
Volume 2021 (2021), Article ID 5527460, 12 pages
https://doi.org/10.1155/2021/5527460

Research Article

A Prediction Model Based on Noninvasive Indicators to Predict the 8-Year Incidence of Type 2 Diabetes in Patients with Nonalcoholic Fatty Liver Disease: A Population-Based Retrospective Cohort Study

Xintian Cai, Qing Zhu, Yuanyuan Cao, Shasha Liu, Mengru Wang, Ting Wu, Jing Hong, Ayguzal Ahmat, Xiayire Aierken, and Nanfang Li

Hypertension Center of People’s Hospital of Xinjiang Uygur Autonomous Region, Xinjiang Hypertension Institute, National Health Committee Key Laboratory of Hypertension Clinical Research, Urumqi, China

Correspondence should be addressed to Nanfang Li; lnanfang2016@sina.com

Received 30 January 2021; Accepted 8 May 2021; Published 15 May 2021

Academic Editor: Ali Abbara

Copyright © 2021 Xintian Cai et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

Background. The prevention of type 2 diabetes (T2D) and its associated complications has become a major priority of global public health. In addition, there is growing evidence that nonalcoholic fatty liver disease (NAFLD) is associated with an increased risk of diabetes. Therefore, the purpose of this study was to develop and validate a nomogram based on independent predictors to better assess the 8-year risk of T2D in Japanese patients with NAFLD. Methods. This is a historical cohort study from a collection of databases that included 2741 Japanese participants with NAFLD without T2D at baseline. All participants were randomized to a training cohort () and a validation cohort (). The data of the training cohort were analyzed using the least absolute shrinkage and selection operator method to screen the suitable and effective risk factors for Japanese patients with NAFLD. A cox regression analysis was applied to build a nomogram incorporating the selected features. The C-index, receiver operating characteristic curve (ROC), calibration plot, decision curve analysis, and Kaplan-Meier analysis were used to validate the discrimination, calibration, and clinical usefulness of the model. The results were reevaluated by internal validation in the validation cohort. Results. We developed a simple nomogram that predicts the risk of T2D for Japanese patients with NAFLD by using the parameters of smoking status, waist circumference, hemoglobin A1c, and fasting blood glucose. For the prediction model, the C-index of training cohort and validation cohort was 0.839 (95% confidence interval (CI), 0.804-0.874) and 0.822 (95% CI, 0.777-0.868), respectively. The pooled area under the ROC of 8-year T2D risk in the training cohort and validation cohort was 0.811 and 0.805, respectively. The calibration curve indicated a good agreement between the probability predicted by the nomogram and the actual probability. The decision curve analysis demonstrated that the nomogram was clinically useful. Conclusions. We developed and validated a nomogram for the 8-year risk of incident T2D among Japanese patients with NAFLD. Our nomogram can effectively predict the 8-year incidence of T2D in Japanese patients with NAFLD and helps to identify people at high risk of T2D early, thus contributing to effective prevention programs for T2D.

1. Introduction

The global prevalence of adult diabetes has increased rapidly in recent decades and has become a major public health problem [1]. Type 2 diabetes (T2D), the most common form of diabetes, is a chronic disease characterized by elevated blood glucose levels due to insufficient insulin production and insulin resistance [2]. In addition, there is growing evidence that nonalcoholic fatty liver disease (NAFLD) is associated with an increased risk of diabetes, independent of traditional diabetes risk factors [3, 4]. As a debilitating chronic epidemic, a core component of T2D prevention strategies is to identify individuals at high risk of T2D [5]. Studies have shown that early identification of people at high risk for T2D and timely lifestyle changes or pharmacological interventions can delay the process of β-cell failure and the development of T2D [6, 7]. Therefore, it is important to study the risk factors for the development of T2D in Japanese patients with NAFLD and to find an easy, reliable, and accurate screening tool for the identification of T2D high-risk groups in NAFLD. This will contribute to the effective implementation of T2D prevention programs in Japanese adults with NAFLD.

Risk prediction models have great potential in the decision-making process for the management of subhealthy populations and patients [8, 9]. Risk prediction models can help guide screening and interventions to predict the onset of diseases. At present, a variety of risk prediction models have been constructed to identify individuals at high risk of T2D, such as Leicester Risk Assessment [10], the Cambridge Risk Score [11], QDiabetes® Calculator [12], and the FINDRISC [13], but they all have a number of limitations. First, most do not take into account lifestyle changes such as physical activity, smoking, and alcohol consumption behaviors. Others are based on invasive and cost-effective data, or small-scale and inappropriate cohort selection. Others are based on short-term follow-up or lack of transparent reporting on the steps that produced the pattern. Most importantly, these diabetes prediction models are based on the general population and limited research focused on individuals at low risk.

The purpose of this paper is to develop a T2D risk prediction model for Japanese patients with NAFLD based on data from the NAGALA cohort study to better screen and assess the 8-year risk of developing T2D in high-risk nondiabetic patients.

2. Materials and Methods

2.1. Data Source

We have downloaded the raw data uploaded by Okamura et al. from the “DATADRYAD” database (http://www.datadryad.org) for free. Since Okamura et al. [14] have granted the data Dryad website ownership of the raw data, we were able to use them for secondary data analysis according to different scientific hypotheses (Dryad data package: 10.5061/dryad.8q0p192).

2.2. Data Description

The raw data variables in the database file included baseline information, incident T2D, and follow-up time. The following variables were extracted: incident T2D, follow-up duration, gender, age, waist circumference (WC), BMI (body mass index), fatty liver, alcohol consumption, smoking status, exercise habits, aspartate transaminase (AST), alanine aminotransferase (ALT), high-density lipoprotein cholesterol (HDL-C), gamma-glutamyl transferase (GGT), triglycerides (TG), total cholesterol (TC), fasting plasma glucose (FPG), hemoglobin A1c (HbA1c), systolic blood pressure (SBP), and diastolic blood pressure (DBP).

2.3. Study Design and Participants

In this study, Okamura et al. [14] used the NAGALA (NAfld in the Gifu Area, Longitudinal Analysis) database to investigate the effect of obesity phenotype on the risk of developing T2D. Since most of the participants require repeated examinations, the researchers conducted a follow-up study of incident T2D diagnosed by blood tests and fatty liver diagnosed by abdominal ultrasound. This study was a secondary analysis of the open data from the NAGALA study. The inclusion criteria for the NAGALA study were described in detail in the original article [15]. Briefly, a total of 15744 participants were selected according to the following exclusion criteria: (1) a lack of important data; (2) the participants had hepatitis B or C virus or fatty liver disease; (3) alcohol intake exceeding 60 g per day for men or 40 g per day for women; (4) the participants took medication at baseline; (5) fasting plasma  mmol/L.

2.4. Data Collection and Measurements

In the initial study, a standardized self-administered questionnaire was used to survey all participants about their medical history and lifestyle factors, including physical activity, alcohol consumption, and smoking habits. Researchers assessed alcohol consumption by asking participants about the type and amount of alcohol consumed per week in the previous month and then estimated the average weekly alcohol consumption. Participants were divided into the following four groups: no or light drinking, <40 g/week; light drinking, 40-140 g/week; moderate drinking, 140-280 g/week; or heavy drinking, >280 g/week. The researchers also divided the participants into three groups based on their smoking status: never smokers, past smokers, or current smokers. Nonsmokers were defined as participants who never smoked, past smokers were defined as participants who smoked in the past but quit before the baseline examination, and current smokers were defined as participants who smoked at the time of the baseline examination. Researchers also investigated participants’ recreational and physical activities. The researchers defined regular exercisers as participants who participated in any type of exercise at least once a week.

2.5. Definitions

T2D was defined as  mmol/mol,  mg/dL, and/or self-reported diabetes during follow-up. NAFLD was defined as having fatty liver demonstrated by abdominal ultrasound.

2.6. Ethical Approval

As this study is the second analysis of existing anonymous data, informed consent of participants is not required. Published paper details the ethical permission [14].

2.7. Statistical Analyses

The study is consistent with the transparent report of the multivariate predictive model of individual prognosis or diagnosis (TRIPOD): the TRIPOD statement [16].

Statistical analyses were performed using R software (version 3.6.3; https://www.R-project.org). First, 2741 Japanese patients with NAFLD were randomly divided into a training cohort of 2058 and a validation cohort of 683 for external validation using the R-caret package, consistent with a theoretical ratio of 3 : 1. Data were expressed as deviation (normal distribution) or median (quartile) (skewed distribution) for continuous variables, and categorical variables were evaluated by calculating frequencies or percentages. Two-sample -tests were used to analyze differences between the training and validation cohorts for normally distributed continuous variables, the Wilcoxon rank-sum test for nonnormally distributed continuous variables, and the chi-square test for categorical variables. And then, the data of the training cohort were analyzed using the least absolute shrinkage and selection operator (LASSO) method to screen the suitable and effective risk factors for Japanese patients with NAFLD. LASSO regression is a method to simplify high-dimensional data. Features with nonzero coefficients can be selected in the LASSO regression model. Next, indicators selected in the LASSO regression model were included in the univariate and multivariate cox regression analysis of risk factors related to T2D, and the hazard ratio (HR) and 95% confidence interval (CI) were calculated. The results of the univariate and multivariate cox regression analyzes were visualized using forest plots. Finally, results of the multivariate cox regression analysis were used to construct a nomogram prediction model. In addition, a variety of validation methods were used to estimate the accuracy of the risk prediction model by using data from the training and validation cohorts, respectively. C-index and receiver operating characteristic (ROC) curve were used to quantify the discrimination performance of nomogram. We plotted and calculated calibration curves using the rms software package, which was used to evaluate the calibration of the T2D risk nomogram and accompanied by a Hosmer-Lemeshow test. Decision curve analysis was performed to determine the clinical application of the T2D risk prediction model: the proportion of true positive results minus the proportion of false positive results, and then, the relative risks of false positive and false negative results were weighted to obtain the net benefits of decision-making. Bootstraps for 1000 resample were performed on the ROC curve, C-index, calibration curve, and decision curve analysis to reduce overfitting deviation. Survival analysis was also performed using the Kaplan-Meier analysis between low-risk and high-risk groups according to the cut-off value of 50%, and the log-rank test was performed to compare survival variance in different groups. All statistical tests were two-sided, and values of <0.05 were considered significant.

3. Results

3.1. Baseline Characteristics of the Study Cohort

A total of 2741 participants were included in this study, of which 2058 were in the training cohort and 683 were in the validation cohort. A flow diagram of studying design is depicted in Figure 1. The overall incidence of T2D was 8.14% (223/2518). In the training and validation cohorts, the incidence of T2D was 157 (7.63%) and 66 (9.66%), respectively. The median follow-up time for the training cohort was 1865 days (quartile: 779-3445), and the median follow-up time for the validation cohort was 2073 days (quartile: 1054-3474). In addition, there were no significant differences observed between the two cohorts. Baseline characteristics of training and validation cohorts are summarized in Table 1.

Figure 1: Flow diagram of study design.

Table 1: Demographic and clinical characteristics of the study population in the training and validation cohorts.

Variables
Total
Training cohort
Validation cohort
value

No. of participants
2741
2058
683

Age (years)

0.054

BMI (kg/m²)

0.898

WC (cm)

0.706

ALT (IU/L)
27.00 (20.00-39.00)
27.00 (20.00-39.00)
26.00 (20.00-38.00)
0.843

AST (IU/L)

21.00 (17.00-26.00)
20.00 (16.00-25.50)
0.805

GGT (IU/L)
23.00 (17.00-35.00)
23.00 (17.00-35.00)
23.00 (16.00-34.00)
0.139

HDL-C (mmol/L)

0.564

TC (mmol/L)

0.559

TG (mmol/L)
1.25 (0.88-1.82)
1.25 (0.88-1.82)
1.25 (0.85-1.82)
0.526

HbA1c (mmol/mol)

0.648

FPG (mg/dL)

0.464

SBP (mmHg)

0.152

DBP (mmHg)

0.169

Gender ( (%))

0.635

    Female
486 (17.73%)
369 (17.93%)
117 (17.13%)

    Male
2255 (82.27%)
1689 (82.07%)
566 (82.87%)

Habit of exercise ( (%))

0.313

    No
2340 (85.37%)
1765 (85.76%)
575 (84.19%)

    Yes
401 (14.63%)
293 (14.24%)
108 (15.81%)

Alcohol consumption ( (%))

<0.001

    None
2088 (76.18%)
1574 (76.48%)
514 (75.26%)

    Light
286 (10.43%)
232 (11.27%)
54 (7.91%)

    Moderate
250 (9.12%)
179 (8.70%)
71 (10.40%)

    Heavy
117 (4.27%)
73 (3.55%)
44 (6.44%)

Smoking status ( (%))

0.022

    Never
1226 (44.73%)
944 (45.87%)
282 (41.29%)

    Past
726 (26.49%)
549 (26.68%)
177 (25.92%)

    Current
789 (28.79%)
565 (27.45%)
224 (32.80%)

Follow-up duration (days)
1902.00 (826.00-3460.00)
1865.00 (778.50-3444.75)
2073.00 (1053.50-3474.00)
0.195

Incident T2D ( (%))

0.092

    No
2518 (91.86%)
1901 (92.37%)
617 (90.34%)

    Yes
223 (8.14%)
157 (7.63%)
66 (9.66%)

Notes: data are presented as (%), , or median (IQR). Abbreviations: WC: waist circumference; BMI: body mass index; AST: aspartate aminotransferase; ALT: alanine aminotransferase; HDL-C: high-density lipoprotein cholesterol; GGT: gamma-glutamyl transferase; TG: triglyceride; TC: total cholesterol; FPG: fasting plasma glucose; HbA1c: hemoglobin A1c; DBP: diastolic blood pressure; SBP: systolic blood pressure; T2D: type 2 diabetes.

3.2. Characteristics of Selection by LASSO Regression Analysis

Seventeen potential risk factors were selected from demography and clinical characteristics and analyzed by LASSO regression (Figures 2(a) and 2(b)). Nonzero characteristic variables were selected based on the statistical approach of the LASSO regression model. Therefore, the number of potential variables was reduced from seventeen to four, including smoking status, WC, FPG, and HbA1c. Table 2 shows the specific coefficients corresponding to the variables of lambda.1se.

(a)

(b)

Figure 2: Demographic and clinical feature selection using the LASSO regression model. (a) Tenfold cross-validated error (first vertical line equals the minimum error, whereas the second vertical line shows the cross-validated error within 1 standard error of the minimum). (b) LASSO coefficient profiles of all the clinical features. A coefficient profile plot was produced against the log(lambda) sequence. Each of the different colored curves in the figure represents the trajectory of each independent variable coefficient. The vertical coordinate is the value of the coefficient, the lower horizontal coordinate is log(lambda), and the upper horizontal coordinate is the number of nonzero coefficients in the model. LASSO: least absolute shrinkage and selection operator; SE: standard error.

Table 2: Coefficients and lambda.1se value of the LASSO regression based on the training cohort.

Factors
Coefficients
Lambda.1se

Smoking status
0.065
0.018

WC (cm)
0.022

FPG (mg/dL)
0.052

HbA1c (mmol/mol)
0.199

Abbreviations: LASSO: least absolute shrinkage and selection operator; WC: waist circumference; FPG: fasting plasma glucose; HbA1c: hemoglobin A1c.

3.3. Univariate and Multivariate cox Regression Analysis in the Training Cohort

Univariate and multivariate cox regression analyses (Figure 3(a) of univariate analysis and Figure 3(b) of multivariate analysis) were performed on 2058 Japanese patients with NAFLD in the training cohort. The results showed that smoking status, WC, FPG, and HbA1c were considered to be independent predictors of T2D ().

Figure 3: Forest plot of univariable and multivariable cox regression analysis for risk factors associated with T2D. HR: hazard ratio; CI: confidence interval.

3.4. Development of the Individualized Prediction Model

We have combined the above four independent predictors into a predictive model and displayed it in the form of a nomogram. As shown in Figure 4, the nomogram is a quantitative and convenient tool. To obtain a personalized 8-year risk of T2D in Japanese patients with NAFLD, a vertical line was drawn from the values on the point scale to assess these points, which were then summed to obtain values for each variable. The sum includes the total score and matches the risk on the bottom axis.

Figure 4: Nomogram for predicting the 8-year risk of T2D in adults with NAFLD. To use the nomogram, an individual patient’s value is located on each variable axis, and a line is drawn upward to determine the number of points received for each variable value. The scores for all variables are then added to obtain the total score, and a vertical line is drawn from the total-points row to estimate the 8-year risk of T2D in adults with NAFLD at the lower line of the nomogram.

3.5. Performance of the Nomogram

The ROC curve and C-index were used to evaluate the discriminatory ability of the prediction model. For the prediction model, the pooled area under the ROC of the nomogram was 0.811 with a sensitivity and specificity of 77.49% and 72.36%, respectively, in the training cohort (Figure 5(a)). It was 0.805 in the validation cohort, with sensitivities and specificities of 69.88% and 75.59%, respectively (Figure 5(b)), which indicates a moderately good performance. The C-index of training cohort and validation cohort was 0.839 (95% CI, 0.804-0.874) and 0.822 (95% CI, 0.777-0.868), respectively. As shown in Table 3, the nomogram showed a good prediction model. The prediction model was calibrated using the calibration curve and the Hosmer-Lemeshow test. From the calibration curves, the prediction model showed a good fit in both the training and validation cohorts (Figures 6(a) and 6(b)). As shown by the Hosmer-Lemeshow test, there was good agreement between the predicted and actual probabilities in both the training and validation cohorts. The decision curve analysis of the training cohort (Figure 7(a)) and the validation cohort (Figure 7(b)) indicated that the application of the prediction model in Japanese patients with NAFLD to predict the risk of T2D incidence is more effective than the intervention-for-all-patients scheme. Each patient was divided into a high-risk or low-risk group according to the cut-off value of 50% predicted by nomogram. Kaplan-Meier survival analysis yielded a significant difference in T2D-free survival probability between the training cohort (Figure 8(a)) and the validation cohort (Figure 8(b)). This stratification could effectively discriminate the T2D prevalence outcomes of the two risk groups in the training and validation cohort.

(a)

(b)

Figure 5: ROC curve of the nomogram in the training and validation cohort. (a) ROC curve of the nomogram in the training cohort. (b) ROC curve of the nomogram in the validation cohort. ROC: receiver operating characteristic; AUC: area under the curve (bootstrap resampling ).

Table 3: C-index in the nomogram based on training cohort and validation cohort.

C-index (95% CI)
Dxy
aDxy
Variance
value
value

Training cohort
0.839 (0.804, 0.874)
0.678
0.678
0.036
19.055
<0.001
2058

Validation cohort
0.822 (0.777, 0.868)
0.644
0.644
0.047
13.827
<0.001
683

(a)

(b)

Figure 6: Calibration curves for the training and validation cohort models. (a) Calibration curve of the nomogram in the training cohort. (b) Calibration curve of the nomogram in the validation cohort. The red curve is a calibration curve corresponding to the actual situation. The blue curve represents the 95% CI range of the calibration curve (bootstrap resampling ).

(a)

(b)

Figure 7: Decision curve analysis of the nomogram in the training (a) and validation cohorts (b). The -axis stands the net benefit. The -axis indicates the threshold probability. The red line represents the nomogram. The black line displays the net benefit of the strategy of treating no patients. The gray line displays the net benefit of the strategy of treating all patients (bootstrap resampling ).

(a)

(b)

Figure 8: Kaplan-Meier curves of risk group stratification based on the predictor from nomogram prediction. Kaplan-Meier curves show the T2D-free survival probability of high risk (blue) and low risk (red) groups between the training cohort (a) and the validation cohort (b). Risk tables are also added.

4. Discussion

With rapid economic development, human lifestyles have changed dramatically worldwide. The prevalence and incidence of T2D are rapidly increasing worldwide [17]. T2D is associated with an increased risk of cardiovascular disease and premature death. It is the main cause of end-stage renal disease, blindness, and nontraumatic amputations resulting from microvascular complications, thus imposing a significant economic burden on society [18, 19]. The risk of T2D is strongly associated with lifestyle, nutritional status, and environmental factors [20]. Several large intervention studies have shown that lifestyle changes or pharmacological interventions targeting people at risk of T2D can effectively prevent or delay the onset of T2D and reduce the risk of death from T2D and its complications [21, 22]. The key to successful intervention is the early identification of people at high risk of T2D [18]. In recent years, numerous T2D risk prediction models have been developed and tested, such as the Australian Type 2 Diabetes Risk Assessment Tool [23], the Cambridge Risk Score [11], the Finnish Risk Score [13], the Framingham Diabetes Risk Score [24], and among others [25–27]. However, there is a lack of a risk prediction model based on cohort study data for the Japanese population, especially for patients with NAFLD. Most of the available T2D risk prediction models are only applicable to the target population, and direct application of these models constructed mainly from populations of European origin may underestimate the risk of developing T2D in the Japanese population. Based on data from the NAGALA cohort study, we aimed to develop a T2D risk prediction model for Japanese patients with NAFLD to identify individuals with high risk of T2D.

In this retrospective cohort study, we developed and validated a nomogram model using cost-effective and easily available parameters to predict the 8-year risk of T2D in Japanese patients with NAFLD and to help clinicians identify high-risk populations for T2D. In the training and validation cohorts, our nomogram has excellent prediction performance and also has excellent consistency on the calibration curve. The decision curve analysis illustrates the clinical application value of nomogram. To the best of our knowledge, this study is the first nomogram to use continuous values instead of segmented values to estimate the risk of type 2 diabetes in Japanese patients with NAFLD. In addition, the nomogram will be of great practical value due to its easily available parameters.

Our prediction model includes four parameters: smoking status, WC, FPG, and HbA1c. These variables identified as risk factors for T2D were consistent with previous studies [28, 29]. In our nomogram, current or past smokers have a higher risk of developing T2D than never-smokers. Numerous epidemiological studies have confirmed that smoking is associated not only with the occurrence of T2D but also with the increased risk of T2D hospitalization and mortality, and the risk increases in a dose-dependent manner with the increase of daily smoking [30–32]. According to the 2014 US secretary of health report, compared with nonsmokers, smoking increases the risk of T2D of active smokers by 30-40%, which indicates that smoking cessation should be emphasized as a basic public health strategy to combat the global diabetes epidemic [33]. The World Health Organization also recognizes that smoking is a preventable risk factor for T2D and agrees to avoid smoking/quit smoking as part of its lifestyle recommendations [30], although there is a lack of a complete understanding of the potential pathways of tobacco abuse, especially the mechanism of pancreatic beta cells. However, data from numerous clinical studies suggest that smoking and nicotine have effects on body composition, pancreatic beta cell function, and peripheral insulin sensitivity [20, 29, 34, 35].

Obesity has become a major global epidemic affecting more than 300 million people, and it is the most important risk factor leading to the onset of T2D [36, 37]. BMI has been used as a surrogate marker of obesity and as one of the predictor variables in most diabetes risk models [38–41]. However, BMI does not reflect central obesity. Compared with BMI, WC has a better predictive value for incident T2D, which is consistent with the results of this study [42–44]. WC is a simple anthropometric parameter of abdominal obesity; it is an indicator of central obesity, whereas BMI is an indicator of general obesity [45, 46]. Central obesity is a recognized risk factor for the metabolic syndrome and is strongly associated with the secretion of adipocytokines and inflammatory cytokines, all of which are strongly associated with an increased risk of developing T2D [47–49]. Elevated WC always leads to an accumulation of abdominal fat and a subsequent increase in free fatty acid levels [50]. Excess of circulating free fatty acids leads to insulin resistance by inhibiting insulin signalling and directly accelerating the rate of hepatic gluconeogenesis, coupled with desensitization of the hepatic regulatory loop involving fatty acids by hypothalamic sensing [51–53]. Therefore, waist circumference was used as the basis for the construction of the predictive model in this study.

The FPG level can reflect the secretion level and function of basal insulin [54]. Many epidemiological studies have shown that baseline FPG levels are highly predictive of T2D, and the elevated baseline FPG levels are closely related to the increased risk of developing T2D [55–58]. HbA1c is a comprehensive measurement method of circulating blood glucose level, which reflects the average blood glucose level over the previous 2-3 months and is used as the gold standard for long-term follow-up of blood glucose control [59, 60]. Compared with the oral glucose tolerance test (OGTT) and 2hPG, the measurement of HbA1c is faster and more convenient and can be measured at any time, regardless of the length of fasting or the composition of the previous meal [61, 62]. Epidemiological studies show that elevated HbA1c in nondiabetic adults is associated with T2D incidence rate, incidence rate of cardiovascular disease, and mortality [63–65]. A recent recommendation by the American Diabetes Association Committee on the Diagnosis and Classification of Diabetes Mellitus advocates the use of HbA1c as a practical and effective test method in the diagnosis of prediabetes to identify high-risk groups, and it may be cost-effective to carry out intensive lifestyle interventions to prevent T2D. The decision of the American Diabetes Association Committee is mainly based on the established association between HbA1c and microvascular diseases [59].

Although our nomogram performed well, several limitations warrant mention. First, we relied on FPG and HbA1c, rather than OGTT, to define incident T2D. However, OGTT is not feasible to carry out this test in a large sample. HbA1c does not need fasting, reflecting the long-term blood glycemic status. In addition, the International Expert Committee also recommended the use of HbA1c to diagnose diabetes [66]. Second, our validation cohort was derived from the same population as the training cohort, which may indicate that the findings were overly optimistic. Future research could verify the performance of this prediction model with other databases. Third, this large-scale cohort study was conducted in Japan. Therefore, whether the results of this study can be generalized to other ethnic groups and some specific groups, such as pregnant women and children, requires further validation by external cohorts. Fourth, the current assessment results might not be satisfactory in practice, and some novel biochemical markers or indicators, especially of genetic factors, could improve the performance of the prediction model in the future. Finally, this report is a second analysis based on the existing database. Although many confounding factors have been adjusted, these potential predictors were not included in our prediction model because data on socioeconomic status, lifestyle (except for smoking, drinking, and exercise habits), disease history (such as cardiovascular disease and chronic kidney disease), family health history (such as diabetes), and specific timing of diabetes diagnosis (especially for self-reported diabetes) were not collected in the database.

5. Conclusion

In summary, we developed and validated a nomogram for the 8-year risk of incident in T2D among Japanese patients with NAFLD, including smoking status, WC, FPG, and HbA1c. Our prediction model can effectively predict the 8-year incidence of T2D in Japanese patients with NAFLD and helps to identify people at high risk of T2D early, thus contributing to effective prevention programs for T2D.

Data Availability

All datasets generated and/or analyzed during the present study are included in this published article and available in Dryad (http://www.datadryad.org/).

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Authors’ Contributions

Xintian Cai and Nanfang Li conceived and designed the study, are involved in data management, wrote the first draft, and revised the final manuscript. Yuanyuan Cao, Qing Zhu, Jing Hong, and Shasha Liu screened data from the database. Ayguzal Ahmat, Xiayire Aierken, Ting Wu, and Mengru Wang contributed to the collection and analysis of the data and the preparation of figures and tables. All authors read and approved the final manuscript.

Acknowledgments

The study was supported by the National Key Research and Development Program (grant no. 2018YFC1311503).

References

D. J. Magliano, J. W. Sacre, J. L. Harding, E. W. Gregg, P. Z. Zimmet, and J. E. Shaw, “Young-onset type 2 diabetes mellitus -- implications for morbidity and mortality,” Nature Reviews. Endocrinology, vol. 16, no. 6, pp. 321–331, 2020.

G. G. Nair, E. S. Tzanakakis, and M. Hebrok, “Emerging routes to the generation of functional β-cells for diabetes mellitus cell therapy,” Nature Reviews. Endocrinology, vol. 16, no. 9, pp. 506–518, 2020.

H. Tilg, A. R. Moschen, and M. Roden, “NAFLD and diabetes mellitus,” Nature Reviews. Gastroenterology & Hepatology, vol. 14, no. 1, pp. 32–42, 2017.

J. M. Hazlehurst, C. Woods, T. Marjot, J. F. Cobbold, and J. W. Tomlinson, “Non-alcoholic fatty liver disease and diabetes,” Metabolism, Clinical and Experimental, vol. 65, no. 8, pp. 1096–1108, 2016.

W. D. Strain and P. M. Paldánius, “Diabetes, cardiovascular disease and the microcirculation,” Cardiovascular Diabetology, vol. 17, no. 1, p. 57, 2018.

R. Zhao, A. L. Hui, X. Yao, C. Zhang, and F. Shen, “Diabetes prevention, early intervention, and nondrug therapy,” Journal of Diabetes Research, vol. 2017, Article ID 5184250, 2 pages, 2017.

B. Brannick and S. Dagogo-Jack, “Prediabetes and cardiovascular disease: pathophysiology and interventions for prevention and risk reduction,” Endocrinology and Metabolism Clinics of North America, vol. 47, no. 1, pp. 33–50, 2018.

F. Lucaroni, D. Cicciarella Modica, M. Macino et al., “Can risk be predicted? An umbrella systematic review of current risk prediction models for cardiovascular diseases, diabetes and hypertension,” BMJ Open, vol. 9, no. 12, article e030234, 2019.

W. Jiang, J. Wang, X. Shen et al., “Establishment and validation of a risk prediction model for early diabetic kidney disease based on a systematic review and meta-analysis of 20 cohorts,” Diabetes Care, vol. 43, no. 4, pp. 925–933, 2020.

L. J. Gray, N. A. Taub, K. Khunti et al., “The Leicester Risk Assessment score for detecting undiagnosed type 2 diabetes and impaired glucose regulation for use in a multiethnic UK setting,” Diabetic Medicine, vol. 27, no. 8, pp. 887–895, 2010.

S. J. Griffin, P. S. Little, C. N. Hales, A. L. Kinmonth, and N. J. Wareham, “Diabetes risk score: towards earlier detection of type 2 diabetes in general practice,” Diabetes/Metabolism Research and Reviews, vol. 16, no. 3, pp. 164–171, 2000.

J. Hippisley-Cox, C. Coupland, J. Robson, A. Sheikh, and P. Brindle, “Predicting risk of type 2 diabetes in England and Wales: prospective derivation and validation of QDScore,” BMJ, vol. 338, no. mar17 2, p. b880, 2009.

J. Lindström and J. Tuomilehto, “The diabetes risk score: a practical tool to predict type 2 diabetes risk,” Diabetes Care, vol. 26, no. 3, pp. 725–731, 2003.

T. Okamura, Y. Hashimoto, M. Hamaguchi, A. Obora, T. Kojima, and M. Fukui, “Ectopic fat obesity presents the greatest risk for incident type 2 diabetes: a population-based longitudinal study,” International Journal of Obesity, vol. 43, no. 1, pp. 139–148, 2019.

M. Hamaguchi, T. Kojima, N. Takeda et al., “The metabolic syndrome as a predictor of nonalcoholic fatty liver disease,” Annals of Internal Medicine, vol. 143, no. 10, pp. 722–728, 2005.

G. S. Collins, J. B. Reitsma, D. G. Altman, and K. G. Moons, “Transparent reporting of a multivariable prediction model for individual prognosis or diagnosis (TRIPOD): the TRIPOD statement,” BMJ, vol. 350, p. g7594, 2015.

N. Lascar, J. Brown, H. Pattison, A. H. Barnett, C. J. Bailey, and S. Bellary, “Type 2 diabetes in adolescents and young adults,” The lancet. Diabetes & endocrinology, vol. 6, no. 1, pp. 69–80, 2018.

Y. Zheng, S. H. Ley, and F. B. Hu, “Global aetiology and epidemiology of type 2 diabetes mellitus and its complications,” Nature Reviews. Endocrinology, vol. 14, no. 2, pp. 88–98, 2018.

A. L. Gloyn and D. J. Drucker, “Precision medicine in the management of type 2 diabetes,” The lancet. Diabetes & endocrinology, vol. 6, no. 11, pp. 891–900, 2018.

H. Kolb and S. Martin, “Environmental/lifestyle factors in the pathogenesis and prevention of type 2 diabetes,” BMC Medicine, vol. 15, no. 1, p. 131, 2017.

M. Captieux, R. Prigge, S. Wild, and B. Guthrie, “Defining remission of type 2 diabetes in research studies: a systematic scoping review,” PLoS Medicine, vol. 17, no. 10, article e1003396, 2020.

N. Wang, F. Zhu, L. Chen, and K. Chen, “Proteomics, metabolomics and metagenomics for type 2 diabetes and its complications,” Life Sciences, vol. 212, pp. 194–202, 2018.

L. Chen, D. J. Magliano, B. Balkau et al., “AUSDRISK: an Australian type 2 diabetes risk assessment tool based on demographic, lifestyle and simple anthropometric measures,” The Medical Journal of Australia, vol. 192, no. 4, pp. 197–202, 2010.

S. G. Wannamethee, A. G. Shaper, L. Lennon, and R. W. Morris, “Metabolic syndrome vs Framingham risk score for prediction of coronary heart disease, stroke, and type 2 diabetes mellitus,” Archives of Internal Medicine, vol. 165, no. 22, pp. 2644–2650, 2005.

X. Cai, Q. Zhu, T. Wu et al., “Development and validation of a novel model for predicting the 5-year risk of type 2 diabetes in patients with hypertension: a retrospective cohort study,” BioMed Research International, vol. 2020, Article ID 9108216, 12 pages, 2020.

S. Reitmeier, S. Kiessling, T. Clavel et al., “Arrhythmic gut microbiome signatures predict risk of type 2 diabetes,” Cell Host & Microbe, vol. 28, no. 2, pp. 258–272.e6, 2020.

L. Wilkinson, N. Yi, T. Mehta, S. Judd, and W. T. Garvey, “Development and validation of a model for predicting incident type 2 diabetes using quantitative clinical data and a Bayesian logistic model: a nationwide cohort and modeling study,” PLoS Medicine, vol. 17, no. 8, article e1003232, 2020.

A. Weisman, G. S. Fazli, A. Johns, and G. L. Booth, “Evolving trends in the epidemiology, risk factors, and prevention of type 2 diabetes: a review,” The Canadian Journal of Cardiology, vol. 34, no. 5, pp. 552–564, 2018.

A. Rawshani, A. Rawshani, and S. Gudbjörnsdottir, “Smoking and other risk factors in type 2 diabetes,” The New England Journal of Medicine, vol. 379, no. 26, p. 2575, 2018.

J. Maddatu, E. Anderson-Baucum, and C. Evans-Molina, “Smoking and the risk of type 2 diabetes,” Translational research : the journal of laboratory and clinical medicine, vol. 184, pp. 101–107, 2017.

A. Pan, Y. Wang, M. Talaei, F. B. Hu, and T. Wu, “Relation of active, passive, and quitting smoking with incident type 2 diabetes: a systematic review and meta-analysis,” The lancet. Diabetes & endocrinology, vol. 3, no. 12, pp. 958–967, 2015.

S. Yuan and S. C. Larsson, “A causal relationship between cigarette smoking and type 2 diabetes mellitus: a Mendelian randomization study,” Scientific Reports, vol. 9, no. 1, article 19342, 2019.

G. W. Warren, A. J. Alberg, A. S. Kraft, and K. M. Cummings, “The 2014 surgeon General’s report: “the health consequences of smoking--50 years of progress”: a paradigm shift in cancer care,” Cancer, vol. 120, no. 13, pp. 1914–1916, 2014.

A. Pan, Y. Wang, M. Talaei, and F. B. Hu, “Relation of smoking with total mortality and cardiovascular events among patients with diabetes mellitus: a meta-analysis and systematic review,” Circulation, vol. 132, no. 19, pp. 1795–1804, 2015.

M. Śliwińska-Mossoń and H. Milnerowicz, “The impact of smoking on the development of diabetes and its complications,” Diabetes & Vascular Disease Research, vol. 14, no. 4, pp. 265–276, 2017.

S. N. Bhupathiraju and F. B. Hu, “Epidemiology of obesity and diabetes and their cardiovascular complications,” Circulation Research, vol. 118, no. 11, pp. 1723–1735, 2016.

A. Boles, R. Kandimalla, and P. H. Reddy, “Dynamics of diabetes and obesity: epidemiological perspective,” Biochimica et Biophysica Acta, Molecular Basis of Disease, vol. 1863, no. 5, pp. 1026–1036, 2017.

D. R. Leitner, G. Frühbeck, V. Yumuk et al., “Obesity and type 2 diabetes: two diseases with a need for combined treatment strategies - EASO can lead the way,” Obesity Facts, vol. 10, no. 5, pp. 483–492, 2017.

E. Sutkowska, K. Sutkowski, M. Sokołowski, E. Franek, and D. S. Sr, “Distribution of the highest plantar pressure regions in patients with diabetes and its association with peripheral neuropathy, gender, age, and BMI: one centre study,” Journal of Diabetes Research, vol. 2019, Article ID 7395769, 2019.

M. K. Svensson, M. Tyrberg, L. Nyström et al., “The risk for diabetic nephropathy is low in young adults in a 17-year follow-up from the Diabetes Incidence Study in Sweden (DISS). Older age and higher BMI at diabetes onset can be important risk factors,” Diabetes/Metabolism Research and Reviews, vol. 31, no. 2, pp. 138–146, 2015.

C. J. Chiu, S. L. Li, C. H. Wu, and Y. F. Du, “BMI trajectories as a harbinger of pre-diabetes or underdiagnosed diabetes: an 18-year retrospective cohort study in Taiwan,” Journal of General Internal Medicine, vol. 31, no. 10, pp. 1156–1163, 2016.

J. Wei, X. Liu, H. Xue, Y. Wang, and Z. Shi, “Comparisons of visceral adiposity index, body shape index, body mass index and waist circumference and their associations with diabetes mellitus in adults,” Nutrients, vol. 11, no. 7, p. 1580, 2019.

D. C. Seo, S. Choe, and M. R. Torabi, “Is waist circumference ≥ 102/88 cm better than body mass index ≥ 30 to predict hypertension and diabetes development regardless of gender, age group, and race/ethnicity? Meta-analysis,” Preventive Medicine, vol. 97, pp. 100–108, 2017.

M. Venkatrao, R. Nagarathna, S. S. Patil, A. Singh, S. K. Rajesh, and H. Nagendra, “A composite of BMI and waist circumference may be a better obesity metric in Indians with high risk for type 2 diabetes: an analysis of NMB-2017, a nationwide cross-sectional study,” Diabetes Research and Clinical Practice, vol. 161, p. 108037, 2020.

Q. Zeng, Y. He, S. Dong et al., “Optimal cut-off values of BMI, waist circumference and waist:height ratio for defining obesity in Chinese adults,” The British Journal of Nutrition, vol. 112, no. 10, pp. 1735–1744, 2014.

E. Oliveros, V. K. Somers, O. Sochor, K. Goel, and F. Lopez-Jimenez, “The concept of normal weight obesity,” Progress in Cardiovascular Diseases, vol. 56, no. 4, pp. 426–433, 2014.

U. Smith, “Abdominal obesity: a marker of ectopic fat accumulation,” The Journal of Clinical Investigation, vol. 125, no. 5, pp. 1790–1792, 2015.

F. Salehinia, H. Abdi, F. Hadaegh et al., “Abdominal obesity phenotypes and incident diabetes over 12 years of follow-up: the Tehran lipid and glucose study,” Diabetes Research and Clinical Practice, vol. 144, pp. 17–24, 2018.

F. Bragg, K. Tang, Y. Guo et al., “Associations of general and central adiposity with incident diabetes in Chinese men and women,” Diabetes Care, vol. 41, no. 3, pp. 494–502, 2018.

A. Tchernof and J. P. Després, “Pathophysiology of human visceral obesity: an update,” Physiological Reviews, vol. 93, no. 1, pp. 359–404, 2013.

T. S. Han, Y. Y. Al-Gindan, L. Govan, C. R. Hankey, and M. Lean, “Associations of BMI, waist circumference, body fat, and skeletal muscle with type 2 diabetes in adults,” Acta Diabetologica, vol. 56, no. 8, pp. 947–954, 2019.

G. Vazquez, S. Duval, D. R. Jacobs Jr., and K. Silventoinen, “Comparison of body mass index, waist circumference, and waist/hip ratio in predicting incident diabetes: a meta-analysis,” Epidemiologic Reviews, vol. 29, no. 1, pp. 115–128, 2007.

Y. Ren, Y. Liu, X. Sun et al., “Cohort study to determine the waist circumference cutoffs for predicting type 2 diabetes mellitus in rural China,” Diabetes/Metabolism Research and Reviews, vol. 34, no. 6, p. e3007, 2018.

J. Lu, J. He, M. Li et al., “Predictive value of fasting glucose, postload glucose, and hemoglobin A1c on risk of diabetes and complications in Chinese adults,” Diabetes Care, vol. 42, no. 8, pp. 1539–1548, 2019.

M. Bozorgmanesh, F. Hadaegh, S. Ghaffari, H. Harati, and F. Azizi, “A simple risk score effectively predicted type 2 diabetes in Iranian adult population: population-based cohort study,” European Journal of Public Health, vol. 21, no. 5, pp. 554–559, 2011.

A. G. Tabák, M. Jokela, T. N. Akbaraly, E. J. Brunner, M. Kivimäki, and D. R. Witte, “Trajectories of glycaemia, insulin sensitivity, and insulin secretion before diagnosis of type 2 diabetes: an analysis from the Whitehall II study,” Lancet, vol. 373, no. 9682, pp. 2215–2221, 2009.

A. Y. Strandberg, T. Pienimäki, K. H. Pitkälä, R. S. Tilvis, V. V. Salomaa, and T. E. Strandberg, “Comparison of normal fasting and one-hour glucose levels as predictors of future diabetes during a 34-year follow-up,” Annals of Medicine, vol. 45, no. 4, pp. 336–340, 2013.

G. Danaei, M. M. Finucane, Y. Lu et al., “National, regional, and global trends in fasting plasma glucose and diabetes prevalence since 1980: systematic analysis of health examination surveys and epidemiological studies with 370 country-years and 27 million participants,” Lancet, vol. 378, no. 9785, pp. 31–40, 2011.

I. S. Okosun, M. Davis-Smith, J. Paul Seale, and J. Ngulefac, “Applicability of a combination of hemoglobin A(1c) and fasting plasma glucose in population-based prediabetes screening,” Journal of Diabetes, vol. 4, no. 4, pp. 407–416, 2012.

W. Jia, “Standardising HbA1c-based diabetes diagnosis: opportunities and challenges,” Expert Review of Molecular Diagnostics, vol. 16, no. 3, pp. 343–355, 2016.

K. Chatzianagnostou, L. Vigna, S. di Piazza et al., “Low concordance between HbA1c and OGTT to diagnose prediabetes and diabetes in overweight or obesity,” Clinical Endocrinology, vol. 91, no. 3, pp. 411–416, 2019.

M. Janghorbani and M. Amini, “Incidence of type 2 diabetes by HbA1c and OGTT: the Isfahan diabetes prevention study,” Acta Diabetologica, vol. 49, no. 1, pp. S73–S79, 2012.

A. P. Carson, C. S. Fox, D. K. McGuire et al., “Low hemoglobin A1c and risk of all-cause mortality among US adults without diabetes,” Circulation. Cardiovascular Quality and Outcomes, vol. 3, no. 6, pp. 661–667, 2010.

J. Lu, W. Wang, M. Li et al., “Associations of hemoglobin A1c with cardiovascular disease and mortality in Chinese adults with diabetes,” Journal of the American College of Cardiology, vol. 72, no. 24, pp. 3224-3225, 2018.

G. Mancini, D. J. Maron, P. M. Hartigan et al., “Lifestyle, glycosylated hemoglobin A1c, and survival among patients with stable ischemic heart disease and diabetes,” Journal of the American College of Cardiology, vol. 73, no. 16, pp. 2049–2058, 2019.

International Expert Committee, “International expert committee report on the role of the A1C assay in the diagnosis of diabetes,” Diabetes Care, vol. 32, no. 7, pp. 1327–1334, 2009.