Abstract

Introduction. DN is a common complication of diabetes. However, diabetes combined with renal injury may involve DN or NDKD, with different treatment schemes. The purpose of our study was to determine the independent risk factors of DN and establish a risk score model to help differentiate DN and NDKD, providing a reference for clinical treatment. Methods. A total of 678 T2D patients who had undergone renal biopsy in four affiliated hospitals of Peking University were consecutively enrolled. Patients were assigned to the DN group and NDKD group according to histopathological results. Seventy percent of patients from PKUFH were randomly assigned to the training group, and the remaining 30% were assigned to the internal validation group. Patients from the other three centers were assigned to the external validation group. We used univariate and multivariate logistic regression analyses to identify independent risk factors of DN in the training group and conducted multivariate logistic regression analysis with these independent risk factors in the training group to find regression coefficients “” to establish a risk score model. Finally, we conducted internal and external validation of the model with ROC curves. Results. Diabetic retinopathy, diabetes years, , 24 h , and no hematuria were independent risk factors (), and each factor scored 2, 1, 1, 1, and 1. We assigned the patients to a low-risk group (0-1 points), a medium-risk group (2-3 points), and a high-risk group (4-6 points), representing unlikely DN, possibly DN, and a high probability of DN, respectively. The AUCs were 0.860, 0.924, and 0.855 for the training, internal validation, and external validation groups, respectively. Conclusion. The risk score model could help differentiate DN and NDKD in a noninvasive manner, reduce the number of renal biopsies, and provide a reference for clinical treatment.

1. Introduction

Globally, approximately 536.6 million (10.5%) people had diabetes in 2021, and the number is expected to rise to 783.2 million (12.2%) by 2045 [1]. Diabetic nephropathy (DN) develops in approximately 40% of patients with type 2 diabetes (T2D) and is the leading cause of end-stage kidney disease (ESKD) [2]. The rate of ESKD remains unchanged despite a decline in other diabetes complications [3]. A better understanding of the characteristics of DN is thus warranted to potentially prevent the progression from DN to ESKD.

In clinical situations, patients with T2D accompanied by chronic kidney disease (CKD) are commonly seen. Apart from DN, nondiabetic kidney disease (NDKD) accounts for a large proportion of cases. A review identified 40 studies (5304 data) worldwide from 1977 to 2019 that examined global renal biopsy in type 2 diabetes patients, and the overall prevalence rates of DN, NDKD, and DN plus NDKD were reported to be 41.3%, 40.6%, and 18.1%, respectively [4]. NDKD (including mixed form) manifests a wide spectrum of pathological lesions across the world, including focal segmental glomerulosclerosis (FSGS), membranous nephropathy (MN), and immunoglobulin A nephropathy (IgAN) [5]. The management of DN has evolved from a simple glycemic management perspective to a multifactorial intensive treatment perspective that is different from various NDKD management approaches [6]. Therefore, it is important to distinguish DN from all CKDs. At present, differentiation between DN and NDKD mainly depends on the gold standard—ultrasound-guided renal biopsy. However, it is an invasive examination, and many patients are unable to undergo it because of abnormal coagulation function, kidney atrophy, isolated kidney, or severe cardiopulmonary insufficiency. Thus, a noninvasive differentiation method is needed.

In our study, by retrospectively analyzing the clinical data and histopathological results of T2D patients, we compared the differences between the clinical characteristics of DN and NDKD, analyzed the clinical characteristics and risk factors of DN, and then established a DN risk score model for T2D patients to help clinicians better differentiate DN and NDKD, providing a reference for clinical treatment.

2. Materials and Methods

2.1. Patients

We consecutively recorded 1,101 patients with T2D who underwent percutaneous renal biopsy from January 1, 2018, to December 31, 2019, at Peking University First Hospital (PKUFH) (289 patients) and from January 1, 2012, to December 31, 2020, at Peking University People’s Hospital (PKUPH) (90 patients), Peking University Third Hospital (PKUTH) (259 patients), and Peking University International Hospital (PKUIH) (40 patients). Patients with type 1 diabetes, steroid diabetes and other special types of diabetes, acute kidney injury, autoimmune diseases such as systemic lupus erythematosus, Sjogren’s syndrome, amyloidosis, hepatitis B-associated glomerulonephritis, renal tumors, drug-related renal injury, urinary tract infections, and unclear diagnoses or missing data were excluded (Figure 1). Finally, 678 patients were enrolled in our study, including 295 DN patients and 383 NDKD patients based on histopathological results. Each center’s ethics committee approved the study after ethical review.

2.2. Data Collection

We collected patients’ demographics, medical history, laboratory results, and renal histopathological results through the medical record system. The data on age, sex, body mass index (BMI), diabetes duration, diabetic complications, other chronic disease histories—such as primary hypertension and dyslipidemia, hemoglobin concentration, urine occult blood, serum creatinine, estimated glomerular filtration rate (eGFR), triglycerides, total cholesterol, urinary microalbumin (MA), albumin creatinine ratio (ACR), 24-hour urinary total protein quantification (24 h UTP), light microscopy, direct immunofluorescence, electron microscopy, etc. are shown in Table 1.

2.3. Renal Histopathology

The biopsy specimens were independently reviewed and scored by two experienced renal pathologists. Light microscopic sections were stained with hematoxylin and eosin, periodic acid-Schiff, Masson’s trichrome, and periodic acid methenamine silver stain. The diagnostic standard of DN was the Tervaert grading standard [7].

2.4. Subgroup

According to the histopathological results, patients were assigned to a DN group and an NDKD group. The DN group contained 178 isolated DN (iDN) patients and 117 DN plus NDKD patients (Figure 1). NDKD included MN and IgA nephropathy. In total, there were 383 patients.

2.5. Statistical Analyses

Descriptive statistics are reported as percentages for categorical data, averages () for continuous, normally distributed data, and medians (interquartile range, IQR) for continuous, nonnormally distributed data. The respective between-group comparisons were performed using tests, tests, and the Mann–Whitney U tests. Independent associations between clinical variables and the risk of diabetic nephropathy were studied with logistic regression analysis. All data were analyzed using IBM SPSS Statistics for Windows version 20.0.

Seventy percent of patients from PKUFH were randomly assigned to the training group (202 patients), and the remaining 30% were assigned to the internal validation group (87 patients). Patients from the other three centers were assigned to the external validation group (389 patients).

We used univariate and multivariate logistic regression analyses to identify independent risk factors for DN in the training group (Table 2).

Multivariate logistic regression analysis was then conducted with the above independent risk factors in the training group to find the regression coefficients “” to establish a risk score model (Table 3). The values of different factors were determined according to the coefficients, and the total risk score was judged by the cumulative values of each factor. The risk score model was internally and externally validated using the validation groups with ROC curve and area under the ROC curve (AUC). The goodness-of-fit (GOF) of the model was tested by the Hosmer–Lemeshow method. All tests were conducted bilaterally, and indicated a statistically significant difference.

3. Results

3.1. General Data

There were 678 patients included in our final study, including 467 males and 211 females, with an average age of (range 15-81) years. The median duration of diabetes was 6 (range 0~33) years. The average diabetes onset age was (range 14~78) years. A total of 542 patients (80.1%) had hypertension, 396 patients (58.4%) had hyperlipidemia, 111 patients (16.4%) had coronary heart disease, and 99 patients (14.6%) had previous cerebral vascular disease. A total of 311 patients (45.9%) had a smoking history, 237 patients (35.0%) had a drinking history, and 226 patients (33.3%) had a family history of T2D. There were 167 patients (57.7%) whose eGFR was <60 ml/min/1.73 m2 and 81 patients (28.0%) whose eGFR was <30 ml/min/1.73 m2. A total of 198 patients (29.2%) had negative urine occult blood. A total of 639 patients (94.2%) had proteinuria in routine urine tests. A total of 106 patients (106/663, 16.0%) had a 24-hour UTP of less than 1 g, 160 patients (160/663, 24.1%) had a 24-hour UTP of ≥1 g and <3 g, and 397 patients (397/663, 59.9%) had a 24-hour UTP of ≥3 g.

There was no significant difference in each variable between the training group and the internal validation group, except for 24 h UTP () (data not shown). Differences between the DN and NDKD groups in the training group were similar to those in the total sample population, except for uric acid and drinking history. Data from the DN and NDKD groups of the training group, internal validation group, and external validation group are shown in Table 1.

3.2. Risk Factors

We explored the optimal cutoff values for all continuous variables and transformed them into binary variables. The specific cutoff values are shown in Table 2.

In the training group, we conducted univariate logistic regression and found significant differences in diabetic retinopathy, diabetes duration, eGFR, urea nitrogen, sex, diabetes onset age, 24 h UTP, hematuria, uric acid, smoking history, and hypertension history between the DN and NDKD groups.

Then, through multivariate logistic regression analysis, we found that diabetic retinopathy, diabetes years, , 24 h , and no hematuria were independent risk factors of DN () (Table 2).

3.3. Risk Score Model

Multivariate logistic regression analysis was reconducted with the above 5 independent risk factors in the training group. According to the regression coefficient “,” 1 point was assigned for every ~1.1 increase (Table 3). According to , diabetic retinopathy counted as 2 points, diabetes years counted as 1 point, counted as 1 point, 24 h counted as 1 point and no hematuria also counted as 1 point (Table 3). The total risk score was between 0 and 6. The patients were classified into low-risk (0-1 point), medium-risk (2-3 points), and high-risk (4-6 points) groups, representing unlikely DN, possibly DN, and high probability of DN, respectively. In the training group, the DN rates were 15.1% (11/73) in the low-risk group, 54.2% (38/70) in the medium-risk group, and 91.5% (54/59) in the high-risk group (Table 4).

3.4. Model Validation
3.4.1. Internal Validation

To verify the diagnostic efficiency of the model, the sensitivity and specificity of different cutoff points were calculated. In the internal validation group, when the cutoff point was 2, the sensitivity of the diagnosis of DN reached 0.957, the specificity was 0.725, and the Jordan index was 0.682. When the cutoff point was 4, the specificity of the diagnosis of DN reached 1. By drawing the ROC curve, (, 95% CI 0.871-0.976) (Figure 2(b)). The GOF was good according to the Hosmer–Lemeshow test (, ). In the internal validation group, the DN rates were 6.5% (2/31) in the low-risk group, 64.5% (20/31) in the medium-risk group, and 100% (25/25) in the high-risk group (Table 4).

3.4.2. External Validation

To further verify the diagnostic efficacy of the risk score model, we enrolled 389 patients from three other hospitals from 2012 to 2020 and applied the risk score model for external validation. In the external validation group, when the cutoff point was 2, the sensitivity of the diagnosis of DN reached 0.906, and the specificity was 0.596. When the cutoff point was 4, the specificity of the diagnosis of DN reached 0.936. By drawing the ROC curve, the (, 95% CI 0.814-0.895) (Figure 2(c)). The Hosmer–Lemeshow test was further used to test the GOF of the model, which showed that the model had a good fitting effect (, ). In the external validation group, the DN rates were 8.5% (13/153) in the low-risk group, 41.6% (57/137) in the medium-risk group, and 81.9% (68/83) in the high-risk group (Table 4).

4. Discussion

According to a previous systematic review, not all T2D patients undergoing renal puncture had DN, 40.6% had NDKD, and 18.1% had DN plus NDKD [4]. Treatment plans for kidney diseases differ with different causes. In DN, the focus is more on the management of blood glucose and strengthening as well as the comprehensive management of related metabolic diseases, and the daily nutritional treatment of DN patients is also different from that of patients with other CKDs [6, 8]. In addition, DN progresses to ESKD faster than other kidney diseases [9]. Therefore, it is important to distinguish DN from all CKDs to guide clinical treatment.

In our study, hepatitis B-associated glomerulonephritis, drug-related renal injury, and other kidney diseases that were clinically easy to diagnose were excluded. However, there are still many NDKDs that are difficult to identify without renal biopsy, such as MN, IgA nephropathy, and ischemic nephropathy. In our study, NDKD accounted for 56.5% of the total population, and DN plus NDKD accounted for 17.3%, similar to previous reports.

As previous studies have shown [10, 11], the well-known predictors for the presence of DN (diabetic retinopathy, diabetes duration, and no hematuria) were also evident in our study. We found that the cutoff value of diabetes duration was 5 years, which is conducive to the earlier identification of DN compared with the “10 years” [12] or “15 years” [13] proposed in previous studies. This study found that was also an indicator, which may be related to the fact that the kidney function of DN patients deteriorates more rapidly than that of patients with other CKDs, thus reflecting a more rapid development and accounting for more cases of ESKD [12, 14, 15]. Nephrotic syndrome is a relatively common clinical manifestation of DN [16, 17], mainly manifested as massive proteinuria, which is consistent with the 24 h found in our study, but massive proteinuria alone cannot indicate DN and may also signify membranous nephropathy, etc. [16], so the combination of the above five indicators is conducive to more accurate identification of DN. It has a higher accuracy than diabetic retinopathy or diabetes duration alone [18, 19].

In recent years, there have been intermittent studies on risk score models for DN diagnosis both domestically and internationally. Liu et al. [20] examined 200 patients with T2D who underwent renal biopsy and constructed a new diagnostic model as follows: . Validation tests determined that the accuracy of the new model was 90.9%. Zhou et al. [21] screened 110 diabetic patients with overt proteinuria but no severe renal failure for renal biopsy, and the diagnostic model was constructed as follows: . Zhang et al. [22] screened 1,030 patients with T2D who had undergone renal biopsy and selected ten variables to build the model using two machine learning methods. The models constructed by Liu and Zhou have a common problem, which is that the PDN formula is complex and not convenient for clinical application. It needs to be converted into a computer program to reflect its value. The model constructed by Zhang et al. contains ten variables and requires complex computer programs.

In our model, there were five variables categorized as 2 diabetes-related (diabetes duration and diabetic retinopathy) and 3 kidney-related (eGFR, 24 h UTP, and hematuria) items; these were all common indicators in clinical practice, even in most grassroots medical institutions. Our model score can be quickly calculated and ranked using simple manual methods. The model was verified internally and externally (spatially and temporally).

In our model, patients in the low-risk group had a 6.5% to 15.1% probability of DN, patients in the medium-risk group had a 41.6% to 64.5% probability of DN, and patients in the high-risk group had an 81.9% to 100% probability of DN. For high-risk group patients, it may be possible to reduce the number of renal biopsies.

All enrolled patients in our study were managed according to the latest guidelines during hospitalization. However, our study is a retrospective study with a wide time span. Since 2012, diabetes treatment drugs have also been updated. In recent years, with the application of kidney protective drugs such as sodium-glucose cotransporter-2 inhibitors (SGLT-2i) and mineralocorticoid receptor antagonists (MRAs) in clinical practice in China, the cutoff point values of DN duration in our model may change. For the same reason, although the main indications for renal biopsy in T2D patients involved in our study were rapid onset or progression of albuminuria or sudden onset of nephrotic syndrome, rapid eGFR decline, glomerular hematuria, active urine sediment, suspicion of other systemic diseases, and patients with type 1 diabetes, short diabetes duration and absence of retinopathy, etc., these indicators may vary slightly with the update of guidelines and individual complex situations that may arise in clinical practice [11, 23, 24]. These issues also exist in previous research models by Liu et al. [2022]. Meanwhile, since DN patients with microalbuminuria do not undergo routine renal biopsy, this study found that, compared to our clinical experience, the DN group exhibited higher 24-hour UTP and poorer renal function. These factors cannot be avoided but may lead to selection bias when enrolling patients.

In addition, Pafundi et al. [25] showed that urine albumin excretion is associated with the highest cardiorenal risk. Therefore, early identification and intervention for DN with MA are more important. Reviewing the renal histopathological characteristics of T2D patients with MA alone and establishing a model may make the model more meaningful for the early diagnosis of DN and facilitate clinical treatment.

5. Conclusions

We found that diabetic retinopathy, diabetes years, , 24 h , and no hematuria were independent risk factors for DN. Our study constructed a simple risk score model to identify DN patients from T2D patients with any renal injury. The risk score model was convenient to calculate and memorize and can help differentiate DN and NDKD in a noninvasive manner, reduce the number of renal biopsies, and provide a reference for clinical treatment.

Data Availability

Due to the limitations of the “Regulations on the Management of Human Genetic Resources of the People’s Republic of China,” this study is unable to provide open raw data.

Ethical Approval

The study protocol was reviewed and approved by the ethics committees of Peking University First Hospital (2020keyan052) and all other institutions, and the study was in accordance with the medical ethics standards of China.

Conflicts of Interest

All authors declare that they have no conflicts of interest relevant to this article.

Authors’ Contributions

Junqing Zhang contributed to the study’s conception and design. All authors contributed to the material preparation and data collection. Yuetong Zhao contributed to the data analysis. The first draft of the manuscript was written by Yuetong Zhao, and all authors commented on previous versions of the manuscript. All authors reviewed the manuscript and approved the final version for submission. Junqing Zhang is the guarantor of this work and, as such, had full access to all data in the study and takes responsibility for the integrity of the data and the accuracy of the data analysis.

Acknowledgments

The skilled technical assistance of Qinghao Liu (Peking University First Hospital) is gratefully acknowledged. The authors also acknowledge all physicians and nurses at each center participating in the enrollment of patients. This study was supported by grants from the National Natural Science Foundation of China (81970701); the National Key Research and Development Program (2016YFC1305401 and 2016YFC1305405); the Beijing Science and Technology Plan Project (D171100002817001); and Project 2020BD040 supported by PKU-Baidu Fund.