Back to Journals » Diabetes, Metabolic Syndrome and Obesity » Volume 15

Development and External Validation of a Nomogram and a Risk Table for Prediction of Type 2 Diabetic Kidney Disease Progression Based on a Retrospective Cohort Study in China

Authors Gao YM , Feng ST, Yang Y, Li ZL, Wen Y, Wang B, Lv LL, Xing GL, Liu BC

Received 2 December 2021

Accepted for publication 2 March 2022

Published 14 March 2022 Volume 2022:15 Pages 799—811

DOI https://doi.org/10.2147/DMSO.S352154

Checked for plagiarism Yes

Review by Single anonymous peer review

Peer reviewer comments 2

Editor who approved publication: Prof. Dr. Juei-Tang Cheng



Yue-Ming Gao,1 Song-Tao Feng,1 Yang Yang,2 Zuo-Lin Li,1 Yi Wen,1 Bin Wang,1 Lin-Li Lv,1 Guo-Lan Xing,2 Bi-Cheng Liu1

1Institute of Nephrology, Zhongda Hospital, Southeast University School of Medicine, Nanjing, 210009, People’s Republic of China; 2Institute of Nephrology, The First Affiliated Hospital of Zhengzhou University, Zhengzhou, 450000, People’s Republic of China

Correspondence: Bi-Cheng Liu, Institute of Nephrology, Zhongda Hospital, Southeast University School of Medicine, 87 Dingjiaqiao Road, Nanjing, Jiangsu Province, 210009, People’s Republic of China, Tel +86-25-83262422, Email [email protected]

Purpose: Diabetic kidney disease (DKD) is the leading cause of end-stage renal disease worldwide. Risk assessment provides information about patient prognosis, contributing to the risk stratification of patients and the rational allocation of medical resources. We aimed to develop a model for individualized prediction of renal function decline in patients with type 2 DKD (T2DKD).
Patients and Methods: In a retrospective observational study, we followed 307 T2DKD patients and evaluated the determinants of 1) risk of doubling in serum creatinine (Scr), 2) risk of eGFR< 15 mL/min/1.73m2 using potential risk factors at baseline. A prediction model represented by a nomogram and a risk table was developed using Cox regression and externally validated in another cohort with 206 T2DKD patients. The discrimination and calibration of the prediction model were evaluated by the concordance index (C-index) and calibration curve, respectively.
Results: Four predictors were selected to establish the final model: Scr, urinary albumin/creatinine ratio, plasma albumin, and insulin treatment. The nomogram achieved satisfactory prediction performance, with a C-index of 0.791 [95% confidence interval (CI) 0.762– 0.820] in the derivation cohort and 0.793 (95% CI 0.746– 0.840) in the external validation cohort. Then, all predictors were scored according to their weightings. A risk table with the highest score of 11.5 was developed. The C-index of the risk table was 0.764 (95% CI: 0.731– 0.797), which was similar to the external validation cohort (0.763; 95% CI: 0.714– 0.812). Additionally, the patients were divided into two groups based on the risk table, and significant differences in the probability of outcome events were observed between the high-risk (score > 2) and low-risk (score ≤ 2) groups in the derivation and external validation cohorts (P < 0.001).
Conclusion: The nomogram and the risk table using readily available clinical parameters could be new tools for bedside prediction of renal function decline in T2DKD patients.

Keywords: diabetic kidney disease, type 2 diabetes, estimated glomerular filtration rate, progression, prediction model

Graphical Abstract:

Introduction

Diabetes mellitus (DM) is a global health challenge and has reached alarming levels.1 Diabetic kidney disease (DKD) is a major microvascular complication of DM, occurring in approximately half of all patients with type 2 DM (T2DM) and one-third with type 1 DM (T1DM).2 DKD is defined as an estimated glomerular filtration rate (eGFR)<60 mL/min/1.73m2 and/or urinary albumin/creatinine ratio (UACR)≥30 mg/g for more than three months caused by diabetes,3,4 which accounts for approximately 1/3 of the whole disease burden of chronic kidney disease (CKD) worldwide.5 In the past few decades, China has witnessed a dramatic increase in diabetes prevalence from an estimation of 0.67% in 1980 to 10.9% in 2013.6,7 Since 2011, the percentage of DKD in China has exceeded CKD related to glomerulonephritis, and the gap is growing.8 Despite advances over the past two decades in postponing the progression of DKD, it has become evident that DKD constitutes the leading cause of end-stage renal disease (ESRD) in many populations,9 which also contributes significantly to the excess all-cause and cardiovascular mortality in T2DM patients.10

Predicting kidney disease progression in DKD patients remains a significant clinical challenge. Although eGFR is a well-established prognostic biomarker for predicting ESRD and cardiovascular events in DKD patients, it is insufficient to predict disease progression in the early stages. Albuminuria also strongly predicts the progression of DKD, but it lacks specificity and sensitivity for predicting ESRD and eGFR decline.11,12 Moreover, albuminuria was reported to follow a remission/regression trajectory rather than a linearly progressive process, which might reduce its prognostic value.13 Furthermore, recent studies have produced an explosion of novel assay-based risk biomarkers for DKD, but their clinical use remains limited because of instability and high cost.14 Clearly, it is necessary to strengthen the risk assessment of type 2 DKD (T2DKD) progression, which will contribute to making clinical decisions and reasonably allocating medical resources.

In general, previous studies have proposed a vast number of risk factors for progressive DKD in T2DM patients.15–22 Given that a patient may have multiple risk factors. Thus, the establishment of multifactorial prediction models is of great significance to simplify T2DKD prognostication. Several prediction models were established to predict DKD progression in T2DM patients. Elley et al15 developed five Cox proportional hazards models to predict fatal or nonfatal ESRD events in patients with T2DM from the New Zealand Diabetes Cohort Study. These models included sex, ethnicity, age, diabetes duration, albuminuria, Scr, systolic blood pressure (SBP), glycosylated hemoglobin (HbA1c), smoking status, and previous cardiovascular disease status, with the highest C-index of 0.89. Jardine et al16 published an equation predicting the risk of major kidney-related events, including seven predictors: eGFR, UACR, SBP, HbA1c, diabetic retinopathy, sex, and educational attainment, with a C-index of 0.847. Wan et al18 developed new gender-specific models to provide a more accurate 5-year ESRD risk prediction for Chinese patients with T2DM than other existing models, which included 11–12 predictors. However, these models were aimed at patients with T2DM instead of DKD. Moreover, these models were presented in the form of complex equations, which may be difficult to be applied in routine clinical practice.

This retrospective observational study aimed to develop and externally validate a model to predict renal function decline in Chinese patients with T2DKD. The prediction model was represented by a nomogram and a risk table, which could be easily applied in clinical practice by using readily available clinical parameters.

Patients and Methods

Study Design

This observational, retrospective cohort study was conducted in accordance with the Transparent Reporting of a Multivariable Prediction Model for Individual Prognosis or Diagnosis (TRIPOD) statement.23

Study Population

In this retrospective cohort study, patients with T2DKD who were hospitalized in the Institute of Nephrology, Zhongda Hospital, Southeast University School of Medicine (baseline from June 2013 to June 2017) were screened. The inclusion criteria included: (1) eGFR<60 mL/min/1.73m2 and/or UACR≥30 mg/g for more than three months caused by T2DM; (2) the follow-up time was more than six months. The exclusion criteria included: (1) patients with eGFR<15 mL/min/1.73m2 or with renal replacement treatment; (2) patients with nondiabetic renal diseases, such as IgA nephropathy or membranous nephropathy; (3) patients with other systemic diseases involving the kidney, such as multiple myeloma, allergic purpura, vasculitis, or systemic lupus erythematosus; (4) baseline data were incomplete. Finally, 307 T2DKD patients were included as the derivation cohort. The flowchart for the selection of patients in the derivation cohort was shown in Supplement Figure 1. Another 206 hospitalized T2DKD patients in the First Hospital Affiliated with Zhengzhou University from January 2018 to January 2020 were enrolled as the external validation cohort. For each patient, the start of observation (origin) was the date of the first hospitalization. This study was approved by the Ethics Committee of Zhongda Hospital, Southeast University (2019ZDSYLL057-P01) and was in accordance with the 1964 Helsinki Declaration and its later amendments or comparable ethical standards. Informed consent was obtained from all participants included in the study.

Study Outcomes

The primary endpoint was defined as the doubling of Scr or eGFR<15 mL/min/1.73m2 up to January 2021 (end date), whatever came first. Doubling of Scr was defined as an increase of 100% and to at least 177 µmol/L, as suggested previously.24 For patients without outcomes, the end date was defined as the time of the latest hospitalization or the latest clinic visit. The follow-up time for each patient was calculated as the number of days between the origin and end date. The Chronic Kidney Disease Epidemiology Collaboration (CKD-EPI) equation was chosen to calculate eGFR based on different Scr values, age, and sex.25,26

Predictors

Predictors were prespecified according to the clinical experience and current literature on established and potential predictors in the progression of DKD.14 Clinical characteristics, including age, sex, body mass index (BMI), SBP, diastolic blood pressure (DBP), duration of DM, smoking, family history of DM, medical history of diabetic retinopathy, hypertension, coronary heart disease (CHD) and stroke, and insulin and renin-angiotensin-aldosterone system inhibitor (RAASi) treatment were collected. Laboratory parameters considered were Scr, eGFR, hemoglobin, plasma albumin (ALB), fasting plasma glucose, total cholesterol, low-density lipoprotein cholesterol (LDL-C), high-density lipoprotein cholesterol (HDL-C), triglyceride, HbA1c, blood urea nitrogen (BUN), serum uric acid, serum cystatin C, 24-h urine protein, and UACR.

Data Collection

Baseline demographics, complications or comorbidities, treatment, laboratory parameters, and outcomes were collected in detail from electronic medical records of these patients. Data were extracted by two trained physicians (Yue-Ming Gao and Song-Tao Feng) and checked by the third experienced physician (Bin Wang). Patients were also contacted via telephone if there were missing data.

Statistical Analysis

Analyses were performed using R software, version 3.5.3 (R Project for Statistical Computing). We used the means ± standard deviation (SD) or medians [interquartile range (IQR)] to express continuous variables, whereas categorical variables were presented as frequencies (%). Continuous variables were compared using the independent sample t-test or the Mann–Whitney U-test, while categorical variables were compared using the chi-square or Fisher’s exact tests. P-values<0.05 (two-sided) were considered significant. Missing data were imputed using random forests-based missing data imputation (R package missForest).27

Cox proportional hazards regression was used to establish the prediction model. The backward step method was used to select the independent variables of the model. A restricted cubic spline (RCS) curve and a minimum p-value method were used to determine the optimal cut-off value of the continuous variables, which were then transformed into categorical variables. A nomogram and a risk table were used to illustrate the model.

The performance of the prediction model was assessed by discrimination and calibration.28 Discrimination was determined by the area under the receiver operating characteristic (ROC) curve (AUC), which ranged from 0.5 (no discrimination) to 1 (perfect discrimination).29 Calibration was determined by a visual calibration plot comparing the predicted and actual probability of T2DKD progression. In addition, patients were categorized into two risk groups based on the risk table, including low- and high-risk. Furthermore, a Kaplan–Meier curve was generated for each risk group to evaluate the predictive power of the risk table.

Results

Study Population and Outcomes

The derivation cohort consisted of 307 eligible patients with a median follow-up of 29 months, and there were 198 outcomes (64.50%). The baseline clinical data of the derivation cohort are shown in Table 1. Compared to patients without outcome, patients with outcomes had higher levels of SBP, DBP, Scr, BUN, total cholesterol, LDL-C, serum uric acid, serum cystatin C, 24-h urine protein and UACR, lower levels of eGFR, hemoglobin and plasma ALB, a higher proportion of insulin treatment, and lower proportion of RAASi treatment. Another 206 eligible patients with a median follow-up of 18 months were chosen as the external validation cohort, and there were 97 outcomes (47%).

Table 1 Baseline Clinical Characteristics of the Derivation Cohort

Selected Factors for Model

The hazard ratio (HR) and p-value (<0.1) of univariate regression analysis were used to select the potential risk factors (Table 2). Patients with higher levels of SBP (HR, 1.012; 95% CI, 1.005–1.019; p=0.001), total cholesterol (HR, 1.170; 95% CI, 1.078–1.269; p<0.001), serum uric acid (HR, 1.002; 95% CI, 1.001–1.003; p=0.001), Scr (HR, 1.010; 95% CI, 1.008–1.013; p<0.001) and UACR (HR, 1.029; 95% CI, 1.020–1.037; p<0.001), and a history of insulin use (HR, 1.490; 95% CI, 1.080–2.055; p=0.015) had a higher risk of renal function decline. Additionally, patients with higher levels of hemoglobin (HR, 0.976; 95% CI, 0.970–0.983; p<0.001) and eGFR (HR, 0.975; 95% CI, 0.968–0.982; p<0.001), and a history of RAASi use (HR, 0.576; 95% CI, 0.432–0.768; p<0.001) had a lower risk of progression. Other predictors, including DBP, LDL-C, BUN, 24-h urine protein, and serum cystatin C, are also statistically significant in univariate Cox regression analysis. Independent variables were selected by a backward stepwise multivariate Cox regression method with a p-value threshold of 0.05 (Figure 1A). Four variables, including the insulin treatment (HR, 1.50; 95% CI, 1.08–2.08; p=0.016), plasma ALB (HR, 0.40; 95% CI, 0.31–0.50; p<0.001), Scr (HR, 1.11; 95% CI, 1.08–1.13; p<0.001), and UACR (HR, 1.25; 95% CI, 1.13–1.39; p<0.001), were incorporated in the final model.

Table 2 Predictors Identified by Univariate Cox Regression Analysis

Figure 1 Forest plots of the predictors selected by multivariate Cox analysis in the derivation cohort. (A) Results of the multivariate Cox regression analysis of the derivation cohort; (B) the multivariate regression results after continuous variables transformed into categorical variables in the derivation cohort. aRepresented the level of ALB was between 25 and 35 g/L; brepresented the level ALB level was less than 25 g/L. Red squares and horizontal bars represent the overall estimates and 95% CIs.

Abbreviations: HR, hazard ratio; CI, confidence interval; Scr, serum creatinine; ALB, albumin; UACR, urinary albumin/creatinine ratio.

A Nomogram of the Final Prediction Model

A nomogram predicting the risk of renal function decline in T2DKD patients is shown in Figure 2. Four predictors, including Scr, insulin treatment, plasma ALB, and UACR, were incorporated in the nomogram. The total score could be calculated by adding up the assigned number of points for each predictor in the nomogram. A higher total score indicated a worse prognosis for each patient. For example, a patient with a normal Scr level (73 µmol/L), lower plasma ALB level (32 g/L), higher UACR level (256 mg/g), and insulin treatment would have a total score of 77.2, indicating a predicted one-year or two-year event-free survival rates of 89% and 78.6%, respectively.

Figure 2 A nomogram of the final prediction model.

Abbreviations: Scr, serum creatinine; ALB, albumin; UACR, urinary albumin/creatinine ratio.

A Risk Table of the Final Prediction Model

Predictors selected in the final model included three continuous variables (Scr, ALB, and UACR) and one categorical variable (insulin treatment). First, we used a minimum p-value approach to evaluate the optimal cut-off value of event-free survival to categorize continuous variables. Additionally, the linearity assumption in continuous variables was examined using RCS curves (Figure 3). Combined with the results of the minimum p-value method and RCS curves, and fully considering the convenience of clinical application, we categorized the above three continuous variables as follows: Scr (µmol/L; <180 or ≥180), plasma ALB (g/L; ≥35, [25–35), or < 25), UACR (mg/g; < 300 or ≥300). The multivariate regression results after continuous variables were transformed into categorical variables are shown in Figure 1B. Next, the score of each variable was assigned according to the HR (Supplement Table 1). Finally, each category was assigned a score, and a risk table was developed as follows: Scr (µmol/L; <180 scores 0 and ≥180 scores 4), plasma ALB (g/L; ≥35 scores 0, 25–35 scores 2, and <25 scores 4), UACR (mg/g; <300 scores 0 and ≥300 scores 2), insulin treatment (no scores 0 and yes scores 1.5), and the highest score was 11.5 (Table 3).

Table 3 A Simple Risk Table of the Final Prediction Model

Figure 3 Restricted cubic spline (RCS) curves of the continuous variables involved in the final model.

Abbreviations: Scr, serum creatinine; ALB, albumin; UACR, urinary albumin/creatinine ratio.

Performance of the Nomogram and the Risk Table

Based on the ROC analysis, the nomogram showed robust discrimination, with an AUC of 0.791 (95% CI 0.762–0.820) in the derivation cohort and 0.793 (95% CI 0.746–0.840) in the validation cohort (Figure 4A). Calibration curves of the nomogram showed that the 24-month risk of the study outcomes predicted by the nomogram was well consistent with the actual probabilities both in the derivation (Figure 5A) and the validation cohort (Figure 5B). The risk table also showed robust discrimination, with an AUC of 0.764 (95% CI 0.731–0.797) in the derivation cohort and 0.763 (95% CI 0.714–0.812) in the validation cohort (Figure 4B). Calibration curves of the risk table showed that the 24-month risk of the outcomes predicted by the risk table was well consistent with the actual probabilities both in the derivation (Figure 5C) and the validation cohort (Figure 5D).

Figure 4 Receiver operating characteristic (ROC) curve for the nomogram and the risk table. The solid blue line represented the derivation cohort, and the solid red line represented the validation cohort. (A) ROC curve for the nomogram. The AUC and its 95% CI were 0.791 (0.762–0.820) in the derivation cohort and 0.793 (0.746–0.840) in the validation cohort; (B) ROC curve for the risk table. The AUC and its 95% CI were 0.764 (0.731–0.797) in the derivation cohort and 0.763 (0.714–0.812) in the validation cohort.

Abbreviations: ROC, receiver operating characteristic; AUC, area under the ROC curve; CI, confidence interval.

Figure 5 The calibration curves for the nomogram and the risk table in the derivation and validation cohort. (A) The calibration curves for the nomogram in the derivation cohort; (B) the calibration curves for the nomogram in the validation cohort; (C) the calibration curves for the risk table in the derivation cohort; (D) the calibration curves for the risk table in the validation cohort. The calibration plot showed the agreement between the predicted probability (x-axis) and the actual probability (y-axis) of the 24-month risk of the study outcomes. A perfect prediction would correspond to the 45° grey dotted line. Spike histograms on the top of each picture reflected the number of T2DKD patients with a predicted probability corresponding to the x-axis value.

Risk Stratification of Risk Table

Based on the risk table, we used the minimum p-value method to separate the total scores and fully considered the clinical application. Patients were finally categorized into two risk groups, including low- (score≤2) and high-risk (score>2). In the low-risk group, the probability of patients who had outcome events was 30.6% in the derivation cohort and 27.7% in the validation cohort, which increased to 69.4% and 72.3% in the high-risk group, respectively (Figure 6C). The prevalence of the event-free survival rate was calculated, and a Kaplan–Meier curve was generated in each risk group both in the derivation (Figure 6A) and the validation cohort (Figure 6B). Compared with the low-risk group, the high-risk group had significantly lower event-free survival rates. Compared with the low-risk group, the HRs of the probability of outcome events in the high-risk group were 3.87 (95% CI: 2.63–5.68) in the derivation cohort and 5.90 (95% CI: 2.73–12.73) in the validation cohort, respectively (p<0.001).

Figure 6 Risk stratification and Kaplan–Meier curve of each risk group in the derivation and validation cohort. (A) Kaplan–Meier curve in the derivation cohort; (B) Kaplan–Meier curve in the validation cohort; (C) risk stratification based on the risk table in the derivation and the validation cohort.

Discussion

In recent years, risk prediction for CKD progression has raised significant attention, and several prediction instruments were developed to guide risk stratification using either traditional regression or novel statistical methods.30–35 However, these prediction models focused on the general CKD population, not on patients with DKD. Due to the much higher risk of progressing to ESRD in DKD patients than other types of CKD,2 the performance of these models is needed to be further elucidated, and prediction models aimed at evaluating the risk of DKD progression are needed.

In this study, we demonstrated that the higher UACR level, the lower plasma ALB level, the higher Scr level and insulin treatment significantly increased the risk of renal function decline in patients with T2DKD. Using these predictors, we developed and externally validated a nomogram and a risk table containing variables obtained routinely in clinical practice, which could serve as practicable bedside tools to identify T2DKD patients at high risk of renal function decline. Moreover, based on the risk table, we could easily divide T2DKD patients into different risk categories, which would be helpful to guide clinical decision making.

Albuminuria is a traditional marker of DKD progression. In predicting the risk of progression to hard renal endpoints, studies based on large cohorts have included albuminuria as an essential predictor in their prediction models. In a retrospective analysis based on the results of the large multicenter New Zealand Diabetes Cohort Study, Elley et al15 developed a series of renal risk models to predict the 5-year risk of ESRD in T2DM patients (baseline median eGFR of 75 mL/min/1.73m2), and they found albuminuria was a strong predictor. In addition, based on the results of the ADVANCE study, Jardine et al16 established a risk score for major kidney-related events in T2DM patients (mean eGFR was 74.6 mL/min/1.73m2), and UACR was one of the most important predictors in the final model. In the Chinese population, Wan et al18 developed a gender-specific 5-year ESRD risk prediction model among 149,333 adults with T2DM, and UACR was identified as an important predictor. Through the Cox regression analysis, we also confirmed the importance of UACR in predicting hard renal endpoints. It is noteworthy that the previous studies were focused on T2DM patients with a preserved renal function instead of DKD. However, we included patients with T2DKD (median eGFR of 46.72 mL/min/1.73m2) in our derivation cohort, which could better characterize renal function decline in T2DKD patients.

Previous studies also investigated risk factors for eGFR decline in T2DM patients. In a prospective observational cohort study, Zoppini et al17 followed 1682 Caucasian patients with T2DM and preserved kidney function and identified that albuminuria was the strongest predictor of annual eGFR decline. Additionally, Lorenzo et al36 also reported that albuminuria was the strongest predictor of faster annual eGFR decline in 153 Caucasian T2DM patients with moderate to advanced CKD. According to the results of a World Health Organization multinational study, the incidence of proteinuria was higher in Chinese diabetic patients when compared with other populations.37,38 Therefore, our study provides further evidence that albuminuria is an essential risk factor of renal function decline in Chinese patients with T2DKD. The results of our study emphasized the importance of UACR in predicting renal function decline in T2DKD patients, which indicates that regular UACR assessment should be recommended as part of the routine test for identifying high-risk patients.

Our study also demonstrated that risk factors including hypoalbuminemia and increase of Scr in addition to albuminuria predicted renal function decline in T2DKD patients. Although correlated with each other, every factor is irreplaceable and indispensable in our final model, and one could not be fully predicted by the others. Previously, serum albumin was reported as a predictor of ESRD in T2DM patients with nephropathy in the RENAAL study,39 and it was independent of and complementary to albuminuria in the final risk score. However, the risk score derived from the RENAAL study was lack of discrimination statistics and external validation. Our study evaluated the discrimination determined by the AUC and then validated our model in an external cohort, further confirming the irreplaceable role of serum albumin in predicting renal function decline in T2DKD patients.

Insulin treatment was also identified as an independent risk factor for T2DKD progression. Previous studies suggested that patients with worse glycemic control were more likely to initiate insulin therapy, and insulin users were inclined to suffer from longer duration of DM, greater incidence of cardiovascular disease, CKD, macroalbuminuria, as well as obesity and dyslipidemia.40,41 Our finding also indicated relatively poorer glycemic control among insulin users compared to patients without insulin treatment (median HbA1c 7.80% IQR [6.60, 9.30] % vs 6.90% IQR [6.13, 8.20] %, Z=−3.197, p=0.001). This finding was in consistent with a previous study conducted in 729 Japanese T2DM patients with preserved kidney function and normoalbuminuria, which demonstrated that the degree of chronic hyperglycemia plays a vital role in aggravating the rapid decline of annual eGFR during a 3-year follow-up period.42 Our results further indicated that insulin treatment might be a risk factor for renal function decline in T2DKD, even with advanced CKD.

Concerning the clinical utility of this model, we developed a nomogram to predict the one-year and two-year event-free survival rates. The nomogram behaved good performance with a C-index of 0.791 in the derivation cohort and 0.793 in the external validation cohort, showing good discriminatory performance, and the calibration was good as well. Then, we transformed the model into a simple risk table. For the continuous variables included in the model, the cut-off values were confirmed by a minimum p-value method and RCS curves, taking into full consideration the significance of the clinical practice. Although with some inevitable loss of precision, the performance remained satisfactory.

There are several limitations to this study. Firstly, although our model was externally validated in another center in China with satisfactory performance, the sample size is relatively small with a relatively short follow-up period. Future studies should include a more extensive study population with a longer follow-up time. Secondly, our study may suffer from bias in report and selection due to the retrospective study design. Research with a prospective study design in the future is warranted. Finally, DKD was mainly clinically diagnosed in our participants with diabetes. Future studies could optimize the inclusion criteria and stratify renal function to achieve a more accurate prediction.

Conclusion

Based on a retrospective cohort study, we developed a risk prediction model for T2DKD progression that integrates four easily accessible clinical risk factors (Scr, ALB, UACR, and insulin treatment) to form a practicable nomogram and risk table. Both the nomogram and risk table showed satisfactory prediction performance in the derivation and the external validation cohort. They could serve as quickly useful bedside tools for predicting the risk of renal function decline, thereby making a more rational clinical decision for patients with T2DKD.

Data Sharing Statement

The datasets used to support the findings of this study are available from the corresponding author upon request.

Ethics Approval and Consent to Participate

This study was approved by the Ethics Committee of Zhongda Hospital, Southeast University (2019ZDSYLL057-P01) and was in accordance with the 1964 Helsinki Declaration and its later amendments or comparable ethical standards. Informed consent was obtained from all participants included in the study.

Acknowledgments

The authors are grateful to the contributors and clinical subjects of this study.

Author Contributions

All authors made a significant contribution to the work reported, whether that is in the conception, study design, execution, acquisition of data, analysis and interpretation, or in all these areas; took part in drafting, revising or critically reviewing the article; gave final approval of the version to be published; have agreed on the journal to which the article has been submitted; and agree to be accountable for all aspects of the work.

Funding

This work was supported by the National Key Research and Development Program [grant number 2018YFC1314000].

Disclosure

The authors report no conflicts of interest in this work.

References

1. Sun H, Saeedi P, Karuranga S, et al. IDF diabetes atlas: global, regional and country-level diabetes prevalence estimates for 2021 and projections for 2045. Diabetes Res Clin Pract. 2022;183:109119.

2. Pugliese G, Penno G, Natali A, et al. Diabetic kidney disease: new clinical and therapeutic issues. Joint position statement of the Italian Diabetes Society and the Italian Society of Nephrology on “The natural history of diabetic kidney disease and treatment of hyperglycemia in patients with type 2 diabetes and impaired renal function”. Nutr Metab Cardiovasc Dis. 2019;29(11):1127–1150.

3. KDOQI. KDOQI clinical practice guidelines and clinical practice recommendations for diabetes and chronic kidney disease. Am J Kidney Dis. 2007;49(2 Suppl 2):S12–S154.

4. Jiang G, Luk AOY, Tam CHT, et al. Progression of diabetic kidney disease and trajectory of kidney function decline in Chinese patients with type 2 diabetes. Kidney Int. 2019;95(1):178–187.

5. GBD Chronic Kidney Disease Collaboration. Global, regional, and national burden of chronic kidney disease, 1990–2017: a systematic analysis for the Global Burden of Disease Study 2017. Lancet. 2020;395(10225):709–733.

6. Wang L, Gao P, Zhang M, et al. Prevalence and ethnic pattern of diabetes and prediabetes in China in 2013. JAMA. 2017;317(24):2515–2523.

7. Wang L, Peng W, Zhao Z, et al. Prevalence and treatment of diabetes in China, 2013–2018. JAMA. 2021;326(24):2498–2506.

8. Zhang L, Long J, Jiang W, et al. Trends in chronic kidney disease in China. N Engl J Med. 2016;375(9):905–906.

9. Cheng HT, Xu X, Lim PS, Hung KY. Worldwide epidemiology of diabetes-related end-stage renal disease, 2000–2015. Diabetes Care. 2021;44(1):89–97.

10. González-Pérez A, Saez M, Vizcaya D, Lind M, Garcia Rodriguez L. Incidence and risk factors for mortality and end-stage renal disease in people with type 2 diabetes and diabetic kidney disease: a population-based cohort study in the UK. BMJ Open Diabetes Res Care. 2021;9(1):e002146.

11. Cea Soriano L, Johansson S, Stefansson B, Rodríguez LA. Cardiovascular events and all-cause mortality in a cohort of 57,946 patients with type 2 diabetes: associations with renal function and cardiovascular risk factors. Cardiovasc Diabetol. 2015;14:38.

12. Colhoun HM, Marcovecchio ML. Biomarkers of diabetic kidney disease. Diabetologia. 2018;61(5):996–1011.

13. Jiang W, Wang J, Shen X, 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. 2020;43(4):925–933.

14. Radcliffe NJ, Seah JM, Clarke M, MacIsaac RJ, Jerums G, Ekinci EI. Clinical predictive factors in diabetic kidney disease progression. J Diabetes Investig. 2017;8(1):6–18.

15. Elley CR, Robinson T, Moyes SA, et al. Derivation and validation of a renal risk score for people with type 2 diabetes. Diabetes Care. 2013;36(10):3113–3120.

16. Jardine MJ, Hata J, Woodward M, et al. Prediction of kidney-related outcomes in patients with type 2 diabetes. Am J Kidney Dis. 2012;60(5):770–778.

17. Zoppini G, Targher G, Chonchol M, et al. Predictors of estimated GFR decline in patients with type 2 diabetes and preserved kidney function. Clin J Am Soc Nephrol. 2012;7(3):401–408.

18. Wan EYF, Fong DYT, Fung CSC, et al. Prediction of new onset of end stage renal disease in Chinese patients with type 2 diabetes mellitus - a population-based retrospective cohort study. BMC Nephrol. 2017;18:257.

19. Fernández-Juárez G, Luño J, Barrio V, et al. 25 (OH) vitamin D levels and renal disease progression in patients with type 2 diabetic nephropathy and blockade of the renin-angiotensin system. Clin J Am Soc Nephrol. 2013;8(11):1870–1876.

20. Chang YH, Lei CC, Lin KC, Chang DM, Hsieh CH, Lee YJ. Serum uric acid level as an indicator for CKD regression and progression in patients with type 2 diabetes mellitus-a 4.6-year cohort study. Diabetes Metab Res Rev. 2016;32(6):557–564.

21. Cheng Y, Shang J, Liu D, Xiao J, Zhao Z. Development and validation of a predictive model for the progression of diabetic kidney disease to kidney failure. Ren Fail. 2020;42(1):550–559.

22. Chang YH, Chang DM, Lin KC, Hsieh CH, Lee YJ. High-density lipoprotein cholesterol and the risk of nephropathy in type 2 diabetic patients. Nutr Metab Cardiovasc Dis. 2013;23(8):751–757.

23. Moons KG, Altman DG, Reitsma JB, et al. Transparent Reporting of a multivariable prediction model for Individual Prognosis or Diagnosis (TRIPOD): explanation and elaboration. Ann Intern Med. 2015;162(1):W1–73.

24. Rossing K, Christensen PK, Hovind P, Tarnow L, Rossing P, Parving HH. Progression of nephropathy in type 2 diabetic patients. Kidney Int. 2004;66(4):1596–1605.

25. Levey AS, Stevens LA, Schmid CH, et al. A new equation to estimate glomerular filtration rate. Ann Intern Med. 2009;150(9):604–612. Erratum in: Ann Intern Med. 2011;155(6):408.

26. Michels WM, Grootendorst DC, Verduijn M, Elliott EG, Dekker FW, Krediet RT. Performance of the Cockcroft-Gault, MDRD, and new CKD-EPI formulas in relation to GFR, age, and body size. Clin J Am Soc Nephrol. 2010;5(6):1003–1009.

27. Stekhoven DJ, Bühlmann P. MissForest–non-parametric missing value imputation for mixed-type data. Bioinformatics. 2012;28(1):112–118.

28. Iasonos A, Schrag D, Raj GV, Panageas KS. How to build and interpret a nomogram for cancer prognosis. J Clin Oncol. 2008;26(8):1364–1370.

29. Bandos AI, Rockette HE, Song T, Gur D. Area under the free-response ROC curve (FROC) and a related summary index. Biometrics. 2009;65(1):247–256.

30. Johnson ES, Thorp ML, Yang X, Charansonney OL, Smith DH. Predicting renal replacement therapy and mortality in CKD. Am J Kidney Dis. 2007;50(4):559–565.

31. Johnson ES, Thorp ML, Platt RW, Smith DH. Predicting the risk of dialysis and transplant among patients with CKD: a retrospective cohort study. Am J Kidney Dis. 2008;52(4):653–660.

32. Goto M, Wakai K, Kawamura T, Ando M, Endoh M, Tomino Y. A scoring system to predict renal outcome in IgA nephropathy: a nationwide 10-year prospective cohort study. Nephrol Dial Transplant. 2009;24(10):3068–3074.

33. Tangri N, Stevens LA, Griffith J, et al. A predictive model for progression of chronic kidney disease to kidney failure. JAMA. 2011;305(15):1553–1559.

34. Tangri N, Inker LA, Hiebert B, et al. A dynamic predictive model for progression of CKD. Am J Kidney Dis. 2017;69(4):514–520.

35. Segal Z, Kalifa D, Radinsky K, et al. Machine learning algorithm for early detection of end-stage renal disease. BMC Nephrol. 2020;21(1):518.

36. Lorenzo V, Saracho R, Zamora J, Rufino M, Torres A. Similar renal decline in diabetic and non-diabetic patients with comparable levels of albuminuria. Nephrol Dial Transplant. 2010;25(3):835–841.

37. Diabetes Drafting Group. Prevalence of small vessel and large vessel disease in diabetic patients from 14 centres. The World Health Organisation multinational study of vascular disease in diabetics. Diabetologia. 1985;28:615–640.

38. Chi ZS, Lee ET, Lu M, Keen H, Bennett PH. Vascular disease prevalence in diabetic patients in China: standardised comparison with the 14 centres in the WHO multinational study of vascular disease in diabetes. Diabetologia. 2001;44(Suppl 2):S82–86.

39. Keane WF, Zhang Z, Lyle PA, et al. Risk scores for predicting outcomes in patients with type 2 diabetes and nephropathy: the RENAAL study. Clin J Am Soc Nephrol. 2006;1(4):761–767.

40. Wong K, Glovaci D, Malik S, et al. Comparison of demographic factors and cardiovascular risk factor control among U.S. adults with type 2 diabetes by insulin treatment classification. J Diabetes Complications. 2012;26(3):169–174.

41. Bouchi R, Babazono T, Yoshida N, et al. Association of albuminuria and reduced estimated glomerular filtration rate with incident stroke and coronary artery disease in patients with type 2 diabetes. Hypertens Res. 2010;33(12):1298–1304.

42. Yokoyama H, Kanno S, Takahashi S, et al. Determinants of decline in glomerular filtration rate in nonproteinuric subjects with or without diabetes and hypertension. Clin J Am Soc Nephrol. 2009;4(9):1432–1440.

Creative Commons License © 2022 The Author(s). This work is published and licensed by Dove Medical Press Limited. The full terms of this license are available at https://www.dovepress.com/terms.php and incorporate the Creative Commons Attribution - Non Commercial (unported, v3.0) License. By accessing the work you hereby accept the Terms. Non-commercial uses of the work are permitted without any further permission from Dove Medical Press Limited, provided the work is properly attributed. For permission for commercial use of this work, please see paragraphs 4.2 and 5 of our Terms.