Your new experience awaits. Try the new design now and help us make it even better

ORIGINAL RESEARCH article

Front. Endocrinol., 12 June 2025

Sec. Clinical Diabetes

Volume 16 - 2025 | https://doi.org/10.3389/fendo.2025.1557166

This article is part of the Research TopicPrevention and Treatment Advancements in Diabetic RetinopathyView all 10 articles

Universal nomogram for predicting referable diabetic retinopathy: a validated model for community and ophthalmic outpatient populations using easily accessible indicators

Niu Dongling&#x;Niu Dongling1†Kang Ziwei,&#x;Kang Ziwei1,2†Sun Juanling,Sun Juanling1,2Zhang LiZhang Li1Wang ChangWang Chang1Lei TingLei Ting1Liu Hongli*Liu Hongli1*Zhang Yanchun,*Zhang Yanchun1,2*
  • 1Xi’an People’s Hospital (Xi’an Fourth Hospital), Affiliated People’s Hospital of Northwest University, Xi’an, Shaanxi, China
  • 2Department of Fundus Surgery (Division IV), Shaanxi Eye Hospital, Xi’an, Shaanxi, China

Purpose: This study aimed to develop and validate a universal nomogram for predicting referable diabetic retinopathy (RDR) in type 2 diabetes mellitus (T2DM) patients, using easily accessible clinical indicators for both community and ophthalmic outpatient populations.

Methods: A cross-sectional study was conducted with 1,830 T2DM patients from 14 communities in Xi’an, Shaanxi, China. Participants completed questionnaires, underwent physical exams, and ophthalmic assessments. Univariate analysis and least absolute shrinkage and selection operator (LASSO) regression identified key predictors for RDR. A nomogram was developed using multivariable logistic regression. Model performance was evaluated through area under the curve (AUC), accuracy, precision, recall, F1 score, Youden index, calibration curves, and decision curve analysis (DCA). The dataset was split into training (80%) and test (20%) sets, with external validation using 123 T2DM outpatients from Shaanxi Eye Hospital.

Results: Seven key predictors were identified: serum creatinine, urea nitrogen, urine glucose, HbA1c, urinary microalbumin, diabetes duration, and systolic blood pressure. The nomogram exhibited moderate predictive accuracy, with AUCs of 0.730 (95% CI: 0.691–0.759), 0.767 (95% CI: 0.704–0.831), and 0.723 (95% CI: 0.610–0.835) for the training, test, and external validation sets, respectively. DCA showed that using the model is beneficial for threshold probabilities between 8% and 72%, supporting its broad clinical utility.

Conclusion: This nomogram, based on readily available clinical indicators, provides a reliable and scalable tool for predicting RDR risk in both community and ophthalmic settings. It offers a practical solution for early detection and personalized management of RDR, with broad applicability and clinical potential.

1 Introduction

Diabetes mellitus (DM) has reached epidemic proportions, currently affecting over 460 million people worldwide and projected to exceed 700 million by 2045 (1, 2). Diabetic retinopathy (DR), one of the most serious microvascular complications of DM, remains the leading cause of irreversible vision impairment in working-age adults (3, 4). Clinical strategies for DR management emphasize three pillars: (1) systemic risk factor control, (2) population-based screening, and (3) targeted treatment for vision-threatening DR (vtDR). However, despite technological advances in teleophthalmology and artificial intelligence–assisted image analysis, current screening programs are constrained by their dependence on dilated fundus examination or retinal photography, both of which require trained personnel and specialized equipment (5). This resource dependency is particularly limiting in primary care settings, where the majority of DR cases are first encountered (6). Diagnostic challenge is most critical for referable diabetic retinopathy (RDR)-typically defined as more than mild nonproliferative DR (NPDR) and/or any stage with macular edema (ME) as defined by the the International Clinical Diabetic Retinopathy (ICDR), namely, any retinal thickening, exudate, or microaneurysm within 1 disc diameter of the fovea (7). In many regions, especially those with limited resources, a mismatch between screening infrastructure and patient needs leads to delayed diagnosis and suboptimal care. Studies have shown that up to 50% of RDR cases remain undetected until advanced stages (8, 9), and that early identification combined with standardized treatment can reduce the risk of blindness by up to 98% (10). Community-based DR screening has consistently demonstrated better cost-effectiveness and accessibility compared to managing late-stage complications, yet its scalability remains challenged by logistical and financial constraints.

In this context, predictive models have emerged as promising adjuncts to improve DR risk stratification and pre-screening efficiency. A variety of risk engines, ranging from simple clinical scoring systems to machine learning–based algorithms, have been proposed for DR prediction (11, 12). Among them, nomograms have gained prominence for their ability to synthesize multiple predictors into intuitive, individualized risk assessments, facilitating evidence-based clinical decision-making (13). Nonetheless, existing DR prediction models face several notable limitations. First, most models rely heavily on cumulative glycemic exposure (e.g., diabetes duration and HbA1c), which, according to large epidemiological datasets, accounts for only ~11% of the variance in DR risk (14). Second, many models incorporate specialized diagnostic parameters, such as optical coherence tomography (OCT), diabetic peripheral neuropathy assessments, or ankle-brachial index, limiting their applicability in primary care or low-resource settings (1517). Third, generalizability is often undermined by methodological constraints: most models are developed and validated in single-center, hospital-based cohorts without external validation, which may introduce selection bias and reduce applicability to broader populations (1521).

To address these gaps, we developed a simplified yet robust RDR risk prediction nomogram based solely on routinely collected clinical data—including blood pressure, glycemic status, renal function markers, and urine glucose—from a large-scale community-based cohort in Xi’an, China. The model was then externally validated in an independent tertiary outpatient population, enabling assessment of its generalizability across healthcare levels and demographic subgroups. By focusing specifically on RDR rather than any DR, the model directly supports clinical referral decisions, which is critical for preventing vision-threatening disease progression. This approach emphasizes clinical usability, implementation feasibility, and translational relevance, especially in primary care environments where access to ophthalmologic resources is limited. Ultimately, our nomogram offers a scalable and cost-effective tool for early risk identification, optimized resource allocation, and enhanced integration of DR screening into routine diabetes management.

2 Methods

2.1 Data source and collection

This cross-sectional study was conducted from 2021 to 2023 across 14 community health service centers in Shaanxi Province, China, and at the ophthalmology outpatient department of Xi’an People’s Hospital (Xi’an Fourth Hospital). The inclusion criteria were: adults (age ≥ 18 years) with type 2 diabetes mellitus (T2DM), diagnosed according to the World Health Organization (22). All participants must have been permanent residents of the selected communities and listed in the community chronic disease registries. They should have undergone diabetic retinopathy screening during the study period. Finally, they must have had access to key clinical and laboratory data for analysis. The exclusion criteria were as follows: patients who were younger than 18 years, those diagnosed with type 1 or gestational diabetes, individuals with insufficient data regarding diabetes duration or other key variables, pregnant or lactating women, and patients with malignancies or severe systemic illnesses that could impact ocular status or data completeness (e.g., end-stage renal disease on dialysis). Additionally, patients with ocular co-morbidities such as severe cataracts, retinal vein occlusion, or advanced age-related macular degeneration were excluded, as these conditions could confound the assessment of diabetic retinopathy (DR). Patients who met the criteria were recruited consecutively, and comprehensive evaluations were conducted using a standardized questionnaire and examination protocol. The study adhered to the Declaration of Helsinki and was approved by the Ethics Committee of Xi’an People’s Hospital (Xi’an Fourth Hospital). All participants provided written informed consent prior to enrollment.

Demographic factors (age, sex), lifestyle factors (smoking and alcohol use, exercise habits), socioeconomic status (education level, occupation, income), and medical history (duration of diabetes, medications, comorbid conditions) were recorded. All measurements were performed by trained healthcare staff. Physical examinations included blood pressure measurements using an electronic sphygmomanometer; anthropometric measurements (height, weight, waist and hip circumferences) for body mass index (BMI) and waist-to-hip ratio calculation. These candidate predictors were selected based on known DR risk factors in the literature and their availability in routine practice (18).

Subsequent to a minimum fasting period of 8 hours, fasting blood and urine samples were collected in the morning. fasting blood and urine samples were obtained in the morning for laboratory analyses. Blood tests included fasting plasma glucose, glycated hemoglobin (HbA1c), and a full biochemical panel: total cholesterol, triglycerides, high-density lipoprotein cholesterol, low-density lipoprotein cholesterol, blood urea nitrogen (BUN), serum creatinine (Scr), uric acid, and estimated glomerular filtration rate (eGFR). A complete blood count was also done (white blood cell count and differential, red blood cell count, hemoglobin, hematocrit, etc.), from which the neutrophil-to-lymphocyte ratio (NLR) was derived. Urine tests included a standard urinalysis (which detects qualitative urine glucose and protein) and measurement of urinary microalbumin (mALB) or albumin-to-creatinine ratio if indicated. All laboratory assays were performed in the hospital’s clinical laboratory center using standardized methods and quality control. These candidate variables were selected based on known risk factors for DR reported in the literature and their availability in routine clinical practice. By focusing on commonly measured indicators, we aimed to maximize the model’s practicality in general healthcare settings.

Visual acuity was measured and intraocular pressure (IOP) was checked using a non-contact tonometer (Icare ic100, Icare Finland Oy, Finland). Slit-lamp biomicroscopy (KJ5S1, Suzhou Kangjie Medical Inc, China) was performed to assess the anterior segment. For fundus evaluation, we obtained digital fundus color photographs for each eye using the Horus DEC 200 device (Miis, Hsinchu, Taiwan), employing a 45-degree double field centered on the fovea and optic disc, with or without pharmacologic pupil dilation. Two experienced retinal specialists, masked to clinical data, independently diagnosed and graded DR by the fundus images in accordance with the International Clinical Diabetic Retinopathy (ICDR) severity scale. In cases of grading discrepancy, a third senior retinal specialist (≥20 years’ experience) adjudicated the findings. Patients were classified as having RDR or non-RDR based on the worst eye. Data were entered into EDC electronic databases, with range and logic checks to minimize entry errors.

2.2 Study design

This investigation analyzed 56 indicators from T2DM patients receiving care at 14 community health service centers in Xi’an, Shaanxi, China. Four variables exhibiting missing rates above 30% (postprandial blood glucose (PBG), hypertension, heart disease, and hyperlipidemia) were omitted from the analysis. For data with missing rates below 30%, imputation was performed using the random forest technique. The dataset from the community cohort (N = 1,830) was randomly split into a training set (80%, N = 1,464) for model development and an internal test set (20%, N = 366) for validation of performance. External validation was conducted using an additional 123 samples from Xi’an People’s Hospital (Xi’an Fourth Hospital), a tertiary hospital. This external test assesses the model’s generalizability to a different patient population and care setting. Figure 1 provides a detailed illustration of the research methodology.

Figure 1
www.frontiersin.org

Figure 1. Diagrammatic representation of the study protocol.

2.3 Statistical analysis

Statistical analyses were executed using R software (version 4.2.3, Vienna, Austria) and Python (version 3.7). Continuous variables were represented as mean ± standard deviation or median with interquartile range (IQR), while categorical variables were expressed as frequencies and percentages. Group comparisons for continuous variables were conducted using the t-test or Mann–Whitney U test as appropriate, whereas categorical variables were analyzed using the chi-square test. Prior to multivariable modeling, we systematically evaluated potential multicollinearity among predictors using both correlation matrices (visualized via heatmaps). Variables demonstrating high multicollinearity (Pearson correlation coefficient r > 0.9) were excluded from further analysis to avoid compromising model accuracy.

The univariate comparisons between the RDR group and non-RDR group within the training set, were first performed using appropriate tests (Student’s t-test or Mann-Whitney U for continuous variables, and chi-square for categorical variables) to identify factors significantly associated with RDR. All variables with P < 0.05 in univariate analysis were considered candidates for multivariable modeling. To further select predictors and avoid overfitting, the least absolute shrinkage and selection operator (LASSO) method were applied on the training data, using 10-fold cross-validation to determine the optimal regularization parameter (λ) (23). LASSO is a penalized regression technique that shrinks less important feature coefficients toward zero, effectively performing variable selection. We evaluated the path of coefficients and the cross-validated error to choose a parsimonious set of predictors (at the λ that minimizes the cross-validation error within one standard error of the minimum). By including all significant covariates in the multivariable model, we controlled for potential confounders; each predictor’s effect was thus estimated while adjusting for the others. Using the predictors retained by LASSO, we fitted a multivariable logistic regression model to the training set to estimate the probability of RDR. The regression coefficients were then used to construct a nomogram, a graphical prediction tool that assigns a point score to each predictor value and yields an individualized risk. The formula for the logistic model (logit function) was presented.

Model performance was evaluated using multiple metrics, including area under the receiver operating characteristic curve (AUROC), accuracy, precision, recall, F1 score, balanced accuracy, Youden index, calibration curves, and decision curve analysis (DCA). Discrimination ability was assessed via AUROC, with values exceeding 0.7 considered indicative of good performance (19). Calibration and clinical utility were evaluated using calibration curves and decision curve analysis, respectively (20, 24, 25). Model calibration, the agreement between predicted probabilities and observed outcomes, was assessed by calibration curves, where predictions were grouped into deciles and plotted against actual outcome proportions; a calibration curve close to the 45° line indicates good calibration. DCA examines the net benefit of using the model across a range of threshold probabilities at which a clinician would intervene. Net benefit is calculated by weighing true positives against false positives, relative to strategies of referring all patients vs. referring none. From the DCA, we identified the range of risk thresholds where the nomogram provides greater net benefit than “treat all” or “treat none” approaches. This range informs where the model is useful in practice. All statistical tests were two-tailed, with statistical significance set at P < 0.05.

3 Results

3.1 Baseline characteristics

A total of 1,830 patients with type 2 diabetes were included in the internal cohort (training + test), of whom 362 (19.78%) diagnosed as RDR. An 80:20 ratio was employed to randomly allocate samples into training and test sets. The training set encompassed 1,464 patients, including 294 (20.08%) RDR cases, while the test set comprised 366 individuals, of which 68 (18.58%) were RDR cases. The external validation cohort (N = 123) consisted of T2DM outpatients, with 23 RDR cases (18.70%) (Supplementary Table S1).

We initially examined 52 candidate variables for association with RDR. Comparison of baseline variables between RDR and non-RDR groups in the combined internal cohort revealed significant differences in multiple demographic and clinical factors. Patients with RDR tended to have lower educational attainment, a higher prevalence of glycosuria, and more frequent use of hypoglycemic drugs or insulin (HDI) compared to non-RDR patients. Poor fasting blood glucose control was more common in the RDR group, as were diabetic nephropathy (DN) and peripheral neuropathy (DPN) diagnoses. A history of hypoglycemic coma (HGC) was also more prevalent among RDR cases (P < 0.05). In terms of continuous variables, the RDR group had significantly higher median systolic blood pressure (SBP) and a higher inflammatory cell count (neutrophils and white blood cell count), reflected in a mildly elevated neutrophil-lymphocyte ratio (NLR). Markers of renal function were worse in RDR patients: median serum creatinine (SCr) and blood urea nitrogen (BUN) were higher, while estimated glomerular filtration rate (eGFR) was lower in RDR vs non-RDR. Likewise, urine microalbumin (mALB) was markedly elevated in the RDR group (median 52.9 vs 27.2 mg/L) and proteinuria (PRO) was more frequent (37.3% vs 23.6%; P < 0.001) in RDR patients. Glycemic indices differed as expected: RDR patients had higher HbA1c (median 8.0% vs 7.0%) and a greater waist-to-hip ratio (WHR) than non-RDR (both P < 0.05). Notably, RDR patients had longer diabetes duration (median 10.0 vs 6.5 years) and were diagnosed with diabetes at a younger age on average (first diagnosis age ~56.0 vs 60.4 years; P < 0.001). These results underscore that the RDR group had more adverse profiles in glycemic control, blood pressure, renal function, and diabetes chronicity (Table 1). Violin plots illustrate continuous variables between groups (Figure 2). A correlation analysis found that hematocrit and hemoglobin were highly collinear (r > 0.9) (Figure 3), so hemoglobin was excluded from further modeling to avoid multicollinearity.

Table 1
www.frontiersin.org

Table 1. Demographic and clinical characteristics of the training and test sets.

Figure 2
www.frontiersin.org

Figure 2. Comparative distribution analysis of continuous variables between training and test datasets.

Figure 3
www.frontiersin.org

Figure 3. Correlation matrix of variables within the training cohort.

3.2 Risk factors screening

The training set patients were classified into RDR and non-RDR groups. Univariate analysis in the training set confirmed that numerous factors were significantly associated with RDR status (P < 0.05), including higher SBP, BUN, NEU, NLR, SCr, mALB, HbA1c, and presence of glycosuria, as well as longer diabetes duration and histories of DN, DPN, poor fasting glycemic control, hypoglycemic coma, and proteinuria (Table 2). These variables were then subjected to LASSO regression analysis (Figure 4). Ten-fold cross-validation optimized the regularization parameter (λ). Figure 4 demonstrates that increasing λ intensifies variable shrinkage penalties, reducing the number of selected variables. Figure 4A illustrates coefficient changes for each variable as λ increases, along with predictor numbers. Figure 4B displays λ logarithmic values (horizontal axis), error values (vertical axis), and selected predictor numbers for each λ value. Dashed lines indicate optimal λ values (λ.min and λ.1se). Two Gaussian LASSO regression models were developed based on different λ criteria. Model I, based on λ.min (λ = 0.002), incorporated ten non-zero coefficient predictors: eGFR, age, SCr, BUN, GLU, HbA1c, mALB, duration, SBP, and NLR (Supplementary Table S2). Model II, using λ.1se (λ = 0.03), included seven non-zero coefficient predictors: SCr, BUN, GLU, HbA1c, mALB, duration, and SBP (Supplementary Table S3). Each of these predictors showed a positive association with RDR risk in multivariate analysis in Model II.

Table 2
www.frontiersin.org

Table 2. Univariate analysis of variables associated with RDR.

Figure 4
www.frontiersin.org

Figure 4. Feature selection utilizing LASSO regression. (A) Optimal lambda parameter identifying four variables with non-zero coefficients. (B) Visualization of partial likelihood deviance against log(lambda), with vertical dashed lines indicating the 1-SE criterion.

3.3 Construction and evaluation of nomogram

Using the coefficients, A nomogram was constructed to facilitate the visualization of the diagnostic model and predict the probability of RDR in an accessible manner. The performance of nomograms were evaluated from the following aspects.

Discrimination: Model I was developed based on the aforementioned 10 variables (Figure 5A). The AUCs were 0.730 (95% CI: 0.697–0.764), 0.762 (95% CI: 0.698–0.825), and 0.678 (95% CI: 0.564–0.793) for the training, test, and external validation sets, respectively (Figures 6A–C). Model II was constructed utilizing the seven variables mentioned above (Figure 5B). The AUCs were 0.730 (95% CI: 0.691–0.759) for the training set, 0.767 (95% CI: 0.704–0.831) for the test set, and 0.723 (95% CI: 0.610–0.835) for the external validation set (Figures 6A–C), indicating consistent moderate accuracy. The results indicate that both models achieve moderate discriminative power. For context, Model II’s AUCs were comparable to Model I in the internal data (~0.730 in training and ~0.760 in test) but notably higher on external validation (0.723 vs 0.678 for Model I), suggesting Model II generalized better to new patients.

Figure 5
www.frontiersin.org

Figure 5. Risk nomogram Model I (A) and Model II (B) for identifying RDR in T2DM patients.

Figure 6
www.frontiersin.org

Figure 6. Receiver operating characteristic curve analysis of RDR risk prediction Model I and Model II in the training set (A, D), test set (B, E) and external validation set (C, F).

Classification metrics: The detailed performance metrics for the two models across the three datasets are summarized in Table 3. Each model was evaluated in terms of accuracy, precision, recall, F1 score, balanced accuracy, Youden index, and cutoff, as presented in Table 4. Using the optimal probability cutoff derived from the training set (Youden’s index optimized at 0.186 for Model II), we evaluated the sensitivity, specificity, and predictive values. In the training set, Model II correctly identified 66.7% of RDR cases (sensitivity 66.7%) while correctly excluding 82.4% of non-RDR cases (specificity 82.4%). The overall accuracy in the training data was 81.6%, and the balanced accuracy (average of sensitivity and specificity) was 74.6%. Owing to the low prevalence of RDR, the positive predictive value (PPV) was modest (17.0%), whereas the negative predictive value (NPV) was very high (97.8%), indicating that a patient predicted as low-risk by the model is very likely to truly be non-RDR. The Youden index for the training model was 0.382, confirming a good trade-off between sensitivity and specificity at the chosen cutoff. When applied to the internal test set (N = 366, 18.6% RDR prevalence), Model II maintained similar performance. The AUC of 0.767 was essentially identical to training, and at the same probability threshold (0.186), sensitivity was ~56% and specificity ~84%. The PPV was ~20.6% and NPV ~96.3% in test set, consistent with expectations given the slightly lower RDR prevalence in the test sample. In the external validation set (N = 123, 18.7% RDR), the nomogram’s discrimination remained good (AUC 0.723). Applying the fixed 0.186 cutoff to this cohort yielded a high sensitivity of 92.8% with specificity of 26.3%. This indicates that the model identified nearly all actual RDR cases in the external sample (only ~7% false negatives), at the expense of more false positives among non-RDR patients. The PPV in external validation was 27.1%, while the NPV was 92.6%. The lower specificity and PPV in the external set reflect that the nomogram, calibrated on the internal cohort, tended to over-predict RDR for the new population. Nonetheless, the high NPV and sensitivity in external validation underscore that Model II is effective as a screening tool -it rarely misses true RDR cases, even if it flags some false positives. Given that missing an RDR case has greater clinical consequence (potentially untreated retinopathy) than over-referral, this operating profile is acceptable for a referral triage model. Overall, the model’s discrimination and classification metrics indicate robust performance, with Model II showing slightly improved generalizability relative to the more complex Model I.

Table 3
www.frontiersin.org

Table 3. Detailed performance metrics for the three datasets in Model I and Model II.

Table 4
www.frontiersin.org

Table 4. Different indicators for evaluating the effectiveness of RDR risk prediction models.

Stepwise forward regression analysis identified ten independent predictors for Model I and seven independent predictors for Model II. Supplementary Tables S2 and S3 present the coefficients of the covariates in both models. In the multivariate analysis, two logistic regression models (Model I and Model II) were constructed for RDR prediction using different predictor combinations. An equation was derived to estimate the probability (P) of RDR occurrence based on the coefficients of the significant predictors as follows:

Model I:

  logit(P)=0.872+(0.002)×eGFR+0×mALB+0.050×Duration+0.012×SBP+0.064×NLR+0.021×Age+0.001×SCr+0.173×BUN+0.222×HbA1c+0.317×GLU

Model II:

 logit(P)=0.672+0×mALB+0.049×Duration+0.012×SBP+0.003×SCr+0.167×BUN+0.222×HbA1c+0.338×GLU

where logit(P) is the linear predictive value (log-odds of RDR). The P (probability) of RDR for a given patient is then obtained as:

P=elogit(P)1+elogit(P)

e is Euler’s number, approximately 2.71828.

Calibration: Calibration curves were employed to evaluate the accuracy of the prediction models, demonstrating that the predicted probabilities closely corresponded to the actual probabilities in Models I and II (Figures 7A, B).

Figure 7
www.frontiersin.org

Figure 7. The calibration curve (A, C) and decision curves analysis (B, D) of Model I and Model II in the training set.

Decision curve analysis: Decision curve analysis (DCA) indicated that the threshold probability for identifying referable diabetic retinopathy (RDR) ranged from 8% to 70% for Model I and 8% to 72% for Model II (Figures 7C, D). Within these respective ranges, the application of both nomogram models provides greater clinical net benefit compared with the alternative strategies of universally referring all patients (“treat-all”) or none (“treat-none”). In clinical practice, the threshold probability reflects the predicted risk at which a clinician would typically consider intervening (i.e., referring the patient to an ophthalmologist for detailed examination). The wide threshold probability range for Model II (8%-72%) indicates its applicability across diverse clinical scenarios.

For instance, if a clinician chooses a threshold of 20% predicted risk to trigger referral (a moderately conservative choice), the decision curve indicates the nomogram will improve patient outcomes by identifying high-risk individuals for timely specialist eye evaluation while avoiding unnecessary referrals in low-risk patients. In practical terms, applying Model II means that patients with a predicted RDR risk above the clinician’s chosen cutoff would be referred for prompt retinal examination (potentially preventing vision loss through early treatment), whereas those below the threshold could be monitored routinely in primary care. Thus, the nomogram provides clinical net benefit over a broad range of decision thresholds, supporting its value as a risk stratification tool for guiding referral decisions. Ultimately, the robust performance and fewer input variables required for Model II support its selection as the finalized nomogram. Its demonstrated utility across a wide range of clinical decision-making scenarios emphasizes its practical value for personalized diabetic retinopathy management, facilitating precise and efficient care delivery.

Risk stratification across DR severity: To further assess the model’s behavior, we analyzed the distribution of nomogram scores and predicted risk across different DR severity levels. Among patients who had RDR, higher risk scores were correlated with more advanced DR stages. In the internal cohort’s RDR cases (moderate NPDR, severe NPDR, or PDR), the nomogram score was significantly higher in severe NPDR cases compared to moderate NPDR (median points 68.97 vs 62.82, P = 0.006). Correspondingly, the median predicted probability of RDR in moderate NPDR was 23.5% (IQR 13.8–37.2%), rising to 30.5% (IQR 17.3–51.3%) in severe NPDR cases. This trend was even more evident in the external validation set’s RDR patients. In the external cohort, the nomogram’s median predicted risk was only 11.8% for those with moderate NPDR, compared to 22.0% in severe NPDR and 48.2% in PDR (overall P < 0.001 across stages) (Supplementary Tables S4-S8). Thus, the model assigned substantially higher risk probabilities to individuals with sight-threatening DR (severe NPDR or PDR) than to those with milder RDR. These findings indicate that Model II not only differentiates patients with any referable retinopathy from those without, but it also stratifies risk within the RDR spectrum, aligning with disease severity. In clinical terms, a higher nomogram score signals not just the presence of RDR but likely a more advanced stage of retinopathy, which could help prioritize patients by urgency. The consistency of this pattern in both internal and external sets supports the model’s validity and its potential utility in identifying patients at highest risk of advanced DR progression.

3.4 Clinical utility of the nomogram

A dynamic nomogram was constructed using the coefficients from Model II (Figure 8A) to evaluate and illustrate the risk of developing RDR in T2DM patients (Figure 8B). Numerical values were assigned to each risk factor. To determine the points for each factor, a perpendicular line was drawn from its corresponding value to the “Points” axis. The points from all factors were subsequently summed to calculate the Total Score, and a vertical line was drawn to the “Total Score” row to assess the probability of RDR occurrence. An exemplar patient from the external dataset is highlighted in red. The varying dimensions of the rectangles demonstrate the differences in the relative proportion of patients across each subgroup (Figure 8B).

Figure 8
www.frontiersin.org

Figure 8. (A) Referable diabetic retinopathy (RDR) risk nomogram. Red indicates that the variable is a risk factor for RDR, and blue indicates that the variable is a protective factor for RDR. The darker the color, the greater the influence of the variable on RDR. (B) Dynamic nomogram. A dynamic nomogram based on Model II was created to predict the risk of developing RDR in T2DM patients. According to the patient’s SCr (100 μmol/L), HbA1c (8.2%), SBP (138 mmHg), GLU(1, positive), BUN (10.7 mmol/L), mALB (2339.27 ug/mg), duration (24 years), the total points was 129, and the predicted probability of RDR was 0.936. Therefore, the patient had an 93.6% chance of developing RDR.

4 Discussion

In this study, we developed and validated a parsimonious nomogram (Model II) for predicting referable diabetic retinopathy in patients with T2DM. The model integrates seven readily available clinical and laboratory variables and demonstrated good discriminative performance. Importantly, we externally validated the model on an independent cohort, confirming its generalizability. The nomogram provides individual risk estimates for moderate-or-worse DR in a form that is easily interpretable by clinicians. Using a risk score tool, frontline providers can stratify diabetic patients by RDR risk and guide timely ophthalmology referrals. This approach addresses the pressing need for scalable, cost-effective DR screening strategies, as traditional fundus examination programs often fail to cover all at-risk patients. To our knowledge, this is one of the first RDR prediction models in a Chinese population that combines community health data with hospital data and undergoes external validation, supporting its broader applicability.

4.1 Innovation and comparison to prior models

Our study offers a significant advancement in diabetic retinopathy (DR) risk prediction by developing a clinically applicable nomogram that achieves robust discriminative performance (AUC ≈ 0.73) using only routinely collected clinical parameters, including HbA1c, systolic blood pressure, renal function markers, and urine glucose. Unlike prior models that often rely on specialized, costly, or technologically demanding indicators, our model is optimized for scalability and accessibility in real-world settings, particularly in primary care. A comparative analysis with existing models (Table 5) highlights this distinction. For instance, Ke et al. (19) utilized 2-hour C-peptide levels and nerve conduction assessments, achieving an AUC of 0.75 but at the cost of requiring non-routine testing (19). Wang et al. (29) reported higher predictive accuracy (AUC > 0.80) using insulin resistance indices and vitamin D levels, yet these biomarkers are not standard in most clinical workflows (20). Conversely, community-based models by Pan et al. (21) and Mo et al. (28) employed readily available variables (AUC ~0.70) but lacked external validation, limiting their generalizability (15, 21).

Table 5
www.frontiersin.org

Table 5. Comparison of predictive models for DR in T2DM patients.

In contrast, our nomogram is built upon a large, community-derived dataset and validated externally in an independent tertiary hospital cohort, ensuring consistent performance across diverse healthcare settings and patient populations. This dual-source design, encompassing both general T2DM populations with early or absent retinopathy and ophthalmic outpatients with more advanced disease, captures a broad spectrum of DR severity (Supplementary Tables S4–S8), enhancing model robustness and applicability.

A key strength of our model lies in its explicit focus on referable DR (RDR), a clinically actionable endpoint that directly informs ophthalmologic referral decisions. While many previous models predict general DR risk without regard to severity, our model offers immediate utility in clinical triage by distinguishing those at risk for vision-threatening stages. Furthermore, by employing LASSO regression for variable selection, we identified a parsimonious set of seven predictors that preserved model performance while enhancing interpretability and implementation feasibility.

The integration of community chronic disease registry data with hospital-based clinical records further underscores the model’s translational relevance and real-world alignment. This methodological approach bridges care levels, yielding a universal risk tool that can support both proactive identification of high-risk individuals in primary care and risk stratification of referred patients in ophthalmology settings.

In summary, Model II strikes a clinically meaningful balance between predictive accuracy, interpretability, and generalizability. Its exclusive reliance on universally obtainable variables, rigorous external validation, and focus on an actionable outcome distinguish it from existing models. This work contributes a practical, scalable framework for risk-based DR screening—particularly suited for integration into population health strategies and electronic medical record systems in resource-limited settings.

4.2 Pathophysiological and clinical implications

The composition of the final model carries implications for both understanding DR and managing it. Glycemic control and hypertension are well-established modifiable risk factors for DR. Elevated HbA1c and prolonged diabetes duration are associated with increased DR risk, reflecting chronic hyperglycemia’s detrimental effects on retinal microvasculature (23). Our findings echo landmark trials (DCCT, UKPDS) which showed that tighter glycemic control markedly reduces DR incidence and progression. For instance, every 1% reduction in HbA1c can lower the risk of DR complications by roughly 30%-35%% (27, 30). Notably, our model includes glycosuria as an independent predictor, suggesting that acute hyperglycemic episodes, indicated by urine glucose positivity, may contribute to retinal damage beyond chronic glycemic burden. This aligns with evidence that glycemic variability can induce oxidative stress and inflammatory responses in retinal tissues (31). Thus, the presence of glycosuria in a patient could identify those experiencing unrecognized hyperglycemic excursions, further elevating DR risk despite similar HbA1c levels. These findings align with the current clinical guidelines advocating tight blood glucose regulation as a cornerstone of DR risk reduction and overall diabetes management (32).

Renal function markers emerged as significant predictors in our nomogram, underscoring the interconnected nature of diabetic microvascular complications (33). In our study, elevated SCr, BUN, and mALB, were associated with higher RDR risk. These findings are consistent with studies that microalbuminuria is a marker for the presence and severity of DR (29, 34, 35). These findings support the concept that nephropathy and retinopathy share common pathogenic mechanisms, such as endothelial dysfunction and chronic inflammation (34, 36). Clinically, this underscores the value of monitoring renal status in diabetic patients not just to manage nephropathy, but also as a window into DR risk (37).

Hypertension, particularly elevated SBP, contributes to retinal capillary damage and ischemia, is a known exacerbating factor for DR. Our nomogram indicate increase in SBP was associated with higher RDR risk (38). Blood pressure control, particularly maintaining SBP in recommended ranges, has been shown to slow DR progression and is advocated by current guidelines. The epidemiologic data indicate that a 10 mmHg SBP reduction yields about a 13% decrease in DR risk. There is also evidence that fluctuations in blood pressure may impact DR progression. Effective blood pressure management is crucial in mitigating DR progression.

The prominence of glycemic control and blood pressure confirms that aggressive management of these modifiable factors remains central to preventing DR. Our findings reinforce existing guidelines that advocate tight control of HbA1c and blood pressure to reduce microvascular complications. Additionally, the inclusion of renal indicators (SCr, BUN, mALB) in the risk score highlights that clinicians should view abnormal kidney findings as a red flag for potential retinopathy. In practice, a diabetic patient with rising creatinine or significant albuminuria should prompt heightened vigilance for DR and perhaps expedited retinal screening. This interdependence of diabetic complications suggests a multidisciplinary care approach: close collaboration between endocrinologists, nephrologists, and ophthalmologists can facilitate early intervention across all complication domains.

4.3 Clinical application

The practical utility of this RDR risk model is highlighted by its ease of use and the decision curve analysis. All predictors in Model II are regularly collected in routine diabetes care, making the nomogram readily implementable without special equipment. A clinician or health worker can simply input a patient’s values into the nomogram (or a calculator or web tool based on it) to obtain an individualized risk estimate. This has significant implications for diabetes management and screening: it enables a shift toward risk-based referral for retinal exams. For example, using our model, a general practitioners or endocrinologists could identify a subset of diabetics who have RDR, say, >20% predicted probability of referable DR-those patients would be prioritized for prompt dilated fundus examination by an ophthalmologist. Conversely, patients scoring below a low threshold (perhaps <10% risk) might be safe to defer immediate specialist referral and continue periodic monitoring, especially in resource-constrained settings. By applying such a strategy, healthcare systems could optimize resource allocation, focusing ophthalmologic services on those most likely to benefit and reducing unnecessary referrals of low-risk individuals. The decision curve analysis confirmed that using the nomogram for referral decisions would confer net benefit across a broad range of threshold choices (from very sensitive to more specific criteria). This flexibility means the tool can adapt to different clinical policies: whether one prefers to “catch” as many cases as possible (lower threshold) or avoid false positives (higher threshold), the model adds value over no model at all. An illustrative scenario is a moderate threshold of 20–30% risk: at that level, the DCA suggests substantial clinical benefit in terms of earlier DR detection and prevention of vision loss. Thus, the nomogram can be tailored to the context – for instance, a community screening program might use a low cutoff to maximize sensitivity, whereas a specialist clinic triaging referrals might use a higher cutoff to ensure specificity. In all cases, the individualized risk estimate fosters a more nuanced approach than one-size-fits-all screening. Moreover, because the model is simple, it could be integrated into electronic health record systems to automatically flag high-risk patients during routine visits (e.g., via an embedded risk score or alert). This kind of integration would make risk stratification seamless and could prompt clinicians in real time to arrange retinal examinations for those flagged at high risk. From a patient education standpoint, the nomogram could be a useful counseling tool. For example, showing a patient how their risk score would drop with a lower HbA1c or blood pressure (nomogram scenarios) might motivate improved adherence to therapy. In essence, our model encapsulates the multifactorial nature of DR risk into a tangible score, which can help communicate risk to patients and drive home the importance of holistic risk factor control. Overall, the real-world application of Model II is as a decision support tool that bridges the gap between primary diabetes care and specialty ophthalmology, ensuring that each patient’s need for retinopathy screening is assessed on the basis of objective risk.

4.4 Limitations and future directions

Despite the strengths of our study, several limitations warrant consideration. First, our nomogram was developed and validated using cross-sectional datasets derived from community-based screenings in Northwest China. This design limits the model’s ability to predict future incidence of RDR, as it captures only a single time-point diagnosis rather than disease progression over time. Furthermore, voluntary participation in these screenings may introduce selection bias, as participants may differ from the broader diabetic population in health behavior, socioeconomic status, or access to care. As such, the generalizability of our findings to other regions, ethnicities, or healthcare systems may be limited. Future prospective, multicenter studies across geographically and demographically diverse populations are essential to improve external validity and ensure broader applicability.

Second, the model primarily incorporates static clinical variables, whereas diabetes is a dynamic condition influenced by time-varying factors such as glycemic variability, blood pressure fluctuations, and renal function changes. Incorporating longitudinal data—such as visit-to-visit HbA1c trajectories, systolic blood pressure variability, or progressive microalbuminuria—may enhance predictive granularity and support the development of temporally adaptive risk models. Such refinements would be particularly valuable for identifying patients at imminent risk of DR progression, enabling more timely interventions.

Third, the absence of certain potentially important predictors—such as dietary habits, physical activity, family history, and genetic markers—may reduce the comprehensiveness of our risk stratification approach. Integrating these variables in future studies could improve the personalization and precision of DR risk prediction, particularly in subgroups with atypical risk profiles.

Fourth, while our model was constructed using logistic regression for its interpretability and ease of implementation, advanced machine learning (ML) methods, including gradient boosting, support vector machines, and deep neural networks, have shown superior predictive performance in some DR studies. Nevertheless, ML models often lack transparency, are resource-intensive, and may pose barriers to routine clinical deployment. Head-to-head comparisons of conventional versus ML-based models—focusing not only on accuracy but also on interpretability, clinical utility, and integration feasibility—will be critical for guiding model selection in real-world practice.

Finally, translating prediction models into clinical benefit requires seamless integration with healthcare systems. Embedding the nomogram within electronic health record (EHR) platforms would enable real-time risk assessment by automatically extracting and updating relevant patient data (e.g., HbA1c, blood pressure, renal indices). Such integration could generate dynamic alerts or risk scores within clinical workflows, facilitating timely referral of high-risk patients and supporting targeted screening strategies. In resource-constrained environments, this approach could significantly optimize allocation of ophthalmologic services and reduce preventable vision loss.

In conclusion, future efforts should focus on prospective validation, incorporation of time-dependent and novel predictors, comparative evaluation of modeling strategies, and EHR-based implementation. These advances will further improve the nomogram’s predictive precision, clinical relevance, and scalability, thereby enhancing its role in personalized, risk-based diabetic retinopathy screening and management.

5 Conclusion

In conclusion, we have developed a clinically useful and interpretable nomogram for referable DR risk that relies on seven routine variables reflecting glycemic control, blood pressure, and kidney health. Model II offers moderate but consistent accuracy in identifying patients with sight-threatening retinopathy, as evidenced by internal and external validations. The model’s strength lies in its simplicity and strong grounding in pathophysiology, it leverages the interrelated risk factors that drive microvascular damage in diabetes. For clinicians, this tool provides a quick risk assessment that can support decision-making on patient referrals for ophthalmologic evaluation. For patients and health systems, its use could mean earlier detection of vision-threatening DR and more efficient use of specialist resources. The broad threshold range yielding net benefit indicates that the nomogram is robust across various clinical risk tolerances. Future refinements, including prospective studies and integration of dynamic risk factors, could further enhance its performance. Nonetheless, even in its current form, this risk score represents an innovative step toward personalized DR screening, helping bridge the gap between primary diabetes care and preventive ophthalmology. By identifying high-risk individuals before irreversible eye damage occurs, such a model can contribute to reducing the burden of diabetic blindness and improving outcomes through timely intervention.

Data availability statement

The original contributions presented in the study are included in the article/Supplementary Material. Further inquiries can be directed to the corresponding authors.

Ethics statement

The studies involving humans were approved by Xi’an People’s Hospital (Xi’an Fourth Hospital). The studies were conducted in accordance with the local legislation and institutional requirements. The participants provided their written informed consent to participate in this study.

Author contributions

ND: Funding acquisition, Writing – review & editing, Conceptualization, Methodology, Software, Writing – original draft. KZ: Software, Writing – review & editing, Data curation, Investigation, Visualization. SJ: Data curation, Investigation, Software, Validation, Writing – review & editing. ZL: Data curation, Investigation, Resources, Writing – review & editing. WC: Data curation, Investigation, Writing – review & editing. LT: Formal analysis, Methodology, Software, Validation, Writing – review & editing. LH: Formal analysis, Methodology, Software, Data curation, Funding acquisition, Investigation, Resources, Supervision, Writing – review & editing. ZY: Funding acquisition, Data curation, Formal analysis, Investigation, Resources, Supervision, Project administration, Validation, Writing – review & editing.

Funding

The author(s) declare that financial support was received for the research and/or publication of this article. This study was funded by the Shaanxi Provincial Key R&D Program (2021SF-162), Xi’an Science and Technology Plan Major Research Project (201805104YX12SF38(3)), International Diabetes Exchange and Practice Special Foundation (Z-2017-26-2302), Scientific Research Incubation Foundation of Xi’an People’s Hospital (LH-18, LH-11), Youth Program of Shaanxi Natural Science Basic Research Project (2023-JC-QN-0885), Scientific Research Program of Xi’an Health Commission (2024qn14), and General Research of Science and Technology of Xi’an (24YXYJ093).

Acknowledgments

We extend our gratitude to all project participants and the investigators for their contributions to data collection.

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Generative AI statement

The author(s) declare that no Generative AI was used in the creation of this manuscript.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fendo.2025.1557166/full#supplementary-material

References

1. Zheng Y, Ley SH, and Hu FB. Global aetiology and epidemiology of type 2 diabetes mellitus and its complications. Nat Rev Endocrinol (2018) 14:88–98. doi: 10.1038/nrendo.2017.151

PubMed Abstract | Crossref Full Text | Google Scholar

2. Sun H, Saeedi P, Karuranga S, Pinkepank M, Ogurtsova K, Duncan BB, et al. IDF Diabetes Atlas: Global, regional and country-level diabetes prevalence estimates for 2021 and projections for 2045. Diabetes Res Clin Pract (2022) 183:109119. doi: 10.1016/j.diabres.2021.109119

PubMed Abstract | Crossref Full Text | Google Scholar

3. Antonetti DA, Klein R, and Gardner TW. Diabetic retinopathy. N Engl J Med (2012) 366:1227–39. doi: 10.1056/NEJMra1005073

PubMed Abstract | Crossref Full Text | Google Scholar

4. Gudlavalleti VS, Shukla R, Batchu T, Malladi B, and Gilbert C. Public health system integration of avoidable blindness screening and management, India. Bull World Health Organ (2018) 96:705–15. doi: 10.2471/BLT.18.212167

PubMed Abstract | Crossref Full Text | Google Scholar

5. Vujosevic S, Aldington SJ, Silva P, Hernandez C, Scanlon P, Peto T, et al. Screening for diabetic retinopathy: new perspectives and challenges. Lancet Diabetes Endocrinol (2020) 8:337–47. doi: 10.1016/S2213-8587(19)30411-5

PubMed Abstract | Crossref Full Text | Google Scholar

6. Watson M, McCluskey PJ, Grigg JR, Kanagasingam Y, Daire J, and Estai M. Barriers and facilitators to diabetic retinopathy screening within Australian primary care. BMC Fam Pract (2021) 22:239. doi: 10.1186/s12875-021-01586-7

PubMed Abstract | Crossref Full Text | Google Scholar

7. Abramoff MD, Folk JC, Han DP, Walker JD, Williams DF, Russell SR, et al. Automated analysis of retinal images for detection of referable diabetic retinopathy. JAMA Ophthalmol (2013) 131:351–57. doi: 10.1001/jamaophthalmol.2013.1743

PubMed Abstract | Crossref Full Text | Google Scholar

8. Wong TY, Sun J, Kawasaki R, Ruamviboonsuk P, Gupta N, Lansingh VC, et al. Guidelines on diabetic eye care: The international council of ophthalmology recommendations for screening, follow-up, referral, and treatment based on resource settings. Ophthalmology (2018) 125:1608–22. doi: 10.1016/j.ophtha.2018.04.007

PubMed Abstract | Crossref Full Text | Google Scholar

9. Emamipour S, van der Heijden A, Nijpels G, Elders P, Beulens J, Postma MJ, et al. A personalised screening strategy for diabetic retinopathy: a cost-effectiveness perspective. Diabetologia (2020) 63:2452–61. doi: 10.1007/s00125-020-05239-9

PubMed Abstract | Crossref Full Text | Google Scholar

10. Lin KY, Hsih WH, Lin YB, Wen CY, and Chang TJ. Update in the epidemiology, risk factors, screening, and treatment of diabetic retinopathy. J Diabetes Investig (2021) 12:1322–25. doi: 10.1111/jdi.13480

PubMed Abstract | Crossref Full Text | Google Scholar

11. Gong D, Fang L, Cai Y, Chong I, Guo J, Yan Z, et al. Development and evaluation of a risk prediction model for diabetes mellitus type 2 patients with vision-threatening diabetic retinopathy. Front Endocrinol (Lausanne) (2023) 14:1244601. doi: 10.3389/fendo.2023.1244601

PubMed Abstract | Crossref Full Text | Google Scholar

12. Zhang C, Zhou L, Ma M, Yang Y, Zhang Y, and Zha X. Dynamic nomogram prediction model for diabetic retinopathy in patients with type 2 diabetes mellitus. BMC Ophthalmol (2023) 23:186. doi: 10.1186/s12886-023-02925-1

PubMed Abstract | Crossref Full Text | Google Scholar

13. Balachandran VP, Gonen M, Smith JJ, and DeMatteo RP. Nomograms in oncology: more than meets the eye. Lancet Oncol (2015) 16:e173–80. doi: 10.1016/S1470-2045(14)71116-7

PubMed Abstract | Crossref Full Text | Google Scholar

14. Hirsch IB and Brownlee M. Beyond hemoglobin A1c–need for additional markers of risk for diabetic microvascular complications. Jama (2010) 303:2291–92. doi: 10.1001/jama.2010.785

PubMed Abstract | Crossref Full Text | Google Scholar

15. Yang H, Xia M, Liu Z, Xing Y, Zhao W, Li Y, et al. Nomogram for prediction of diabetic retinopathy in patients with type 2 diabetes mellitus: A retrospective study. J Diabetes Complications (2022) 36:108313. doi: 10.1016/j.jdiacomp.2022.108313

PubMed Abstract | Crossref Full Text | Google Scholar

16. Yang J and Jiang S. Development and validation of a model that predicts the risk of diabetic retinopathy in type 2 diabetes mellitus patients. Acta Diabetol (2023) 60:43–51. doi: 10.1007/s00592-022-01973-1

PubMed Abstract | Crossref Full Text | Google Scholar

17. Li Y, Li C, Zhao S, Yin Y, Zhang X, and Wang K. Nomogram for prediction of diabetic retinopathy among type 2 diabetes population in Xinjiang, China. Diabetes Metab Syndr Obes (2022) 15:1077–89. doi: 10.2147/DMSO.S354611

PubMed Abstract | Crossref Full Text | Google Scholar

18. Chen X, Xie Q, Zhang X, Lv Q, Liu X, and Rao H. Nomogram prediction model for diabetic retinopathy development in type 2 diabetes mellitus patients: A retrospective cohort study. J Diabetes Res (2021) 2021:3825155. doi: 10.1155/2021/3825155

PubMed Abstract | Crossref Full Text | Google Scholar

19. Ke J, Li K, and Cao B. A nomogram for predicting vision-threatening diabetic retinopathy among mild diabetic retinopathy patients: A case-control and prospective study of type 2 diabetes. Diabetes Metab Syndr Obes (2023) 16:275–83. doi: 10.2147/DMSO.S394607

PubMed Abstract | Crossref Full Text | Google Scholar

20. Wang Q, Zeng N, Tang H, Yang X, Yao Q, Zhang L, et al. Diabetic retinopathy risk prediction in patients with type 2 diabetes mellitus using a nomogram model. Front Endocrinol (Lausanne) (2022) 13:993423. doi: 10.3389/fendo.2022.993423

PubMed Abstract | Crossref Full Text | Google Scholar

21. Pan H, Sun J, Luo X, Ai H, Zeng J, Shi R, et al. A risk prediction model for type 2 diabetes mellitus complicated with retinopathy based on machine learning and its application in health management. Front Med (Lausanne) (2023) 10:1136653. doi: 10.3389/fmed.2023.1136653

PubMed Abstract | Crossref Full Text | Google Scholar

22. Gregg EW, Buckley J, Ali MK, Davies J, Flood D, Mehta R, et al. Improving health outcomes of people with diabetes: target setting for the WHO global diabetes compact. Lancet (2023) 401:1302–12. doi: 10.1016/S0140-6736(23)00001-6

PubMed Abstract | Crossref Full Text | Google Scholar

23. Larroumet A, Rigo M, Lecocq M, Delyfer MN, Korobelnik JF, Monlun M, et al. Previous dramatic reduction of HbA1c and retinopathy in type 2 diabetes. J Diabetes Complications (2020) 34:107604. doi: 10.1016/j.jdiacomp.2020.107604

PubMed Abstract | Crossref Full Text | Google Scholar

24. Austin PC, Harrell FJ, and van Klaveren D. Graphical calibration curves and the integrated calibration index (ICI) for survival models. Stat Med (2020) 39:2714–42. doi: 10.1002/sim.8570

PubMed Abstract | Crossref Full Text | Google Scholar

25. Fitzgerald M, Saville BR, and Lewis RJ. Decision curve analysis. Jama (2015) 313:409–10. doi: 10.1001/jama.2015.37

PubMed Abstract | Crossref Full Text | Google Scholar

26. Tao R, Yu X, Lu J, Wang Y, Lu W, Zhang Z, et al. A deep learning nomogram of continuous glucose monitoring data for the risk prediction of diabetic retinopathy in type 2 diabetes. Phys Eng Sci Med (2023) 46:813–25. doi: 10.1007/s13246-023-01254-3

PubMed Abstract | Crossref Full Text | Google Scholar

27. Liu Y, Zhou Z, Wang Z, Yang H, Zhang F, and Wang Q. Construction and clinical validation of a nomogram-based predictive model for diabetic retinopathy in type 2 diabetes. Am J Transl Res (2023) 15:6083–94.

Google Scholar

28. Mo R, Shi R, Hu Y, and Hu F. Nomogram-based prediction of the risk of diabetic retinopathy: A retrospective study. J Diabetes Res (2020) 2020:7261047. doi: 10.1155/2020/7261047

PubMed Abstract | Crossref Full Text | Google Scholar

29. Wang L, Jin L, Wang W, Gong X, Li Y, Li W, et al. Association of renal function with diabetic retinopathy and macular oedema among Chinese patients with type 2 diabetes mellitus. Eye (Lond) (2023) 37:1538–44. doi: 10.1038/s41433-022-02173-5

PubMed Abstract | Crossref Full Text | Google Scholar

30. Stratton IM, Adler AI, Neil HA, Matthews DR, Manley SE, Cull CA, et al. Association of glycaemia with macrovascular and microvascular complications of type 2 diabetes (UKPDS 35): prospective observational study. Bmj (2000) 321:405–12. doi: 10.1136/bmj.321.7258.405

PubMed Abstract | Crossref Full Text | Google Scholar

31. Zhai L, Lu J, Cao X, Zhang J, Yin Y, and Tian H. Association between the variability of glycated hemoglobin and retinopathy in patients with type 2 diabetes mellitus: A meta-analysis. Horm Metab Res (2023) 55:103–13. doi: 10.1055/a-1931-4400

PubMed Abstract | Crossref Full Text | Google Scholar

32. Diabetes. AA. 6. glycemic targets: Standards of medical care in diabetes-2021. Diabetes Care (2021) 44:S73–84. doi: 10.2337/dc21-S006

PubMed Abstract | Crossref Full Text | Google Scholar

33. Hwang H, Kim JY, Oh TK, Chae JB, and Kim DY. Relationship between clinical features of diabetic retinopathy and systemic factors in patients with newly diagnosed type II diabetes mellitus. J Korean Med Sci (2020) 35:e179. doi: 10.3346/jkms.2020.35.e179

PubMed Abstract | Crossref Full Text | Google Scholar

34. Romero-Aroca P, Baget-Bernaldiz M, Reyes-Torres J, Fernandez-Ballart J, Plana-Gil N, Mendez-Marin I, et al. Relationship between diabetic retinopathy, microalbuminuria and overt nephropathy, and twenty-year incidence follow-up of a sample of type 1 diabetic patients. J Diabetes Complications (2012) 26:506–12. doi: 10.1016/j.jdiacomp.2012.06.010

PubMed Abstract | Crossref Full Text | Google Scholar

35. Park HC, Lee YK, Cho A, Han CH, Noh JW, Shin YJ, et al. Diabetic retinopathy is a prognostic factor for progression of chronic kidney disease in the patients with type 2 diabetes mellitus. PloS One (2019) 14:e220506. doi: 10.1371/journal.pone.0220506

PubMed Abstract | Crossref Full Text | Google Scholar

36. Rao X, Wan M, Qiu C, and Jiang C. Role of cystatin C in renal damage and the optimum cut-off point of renal damage among patients with type 2 diabetes mellitus. Exp Ther Med (2014) 8:887–92. doi: 10.3892/etm.2014.1815

PubMed Abstract | Crossref Full Text | Google Scholar

37. Ng WY, Teo BW, Tai ES, Sethi S, Lamoureux E, Tien YW, et al. Cystatin C, chronic kidney disease and retinopathy in adults without diabetes. Eur J Prev Cardiol (2016) 23:1413–20. doi: 10.1177/2047487316637182

PubMed Abstract | Crossref Full Text | Google Scholar

38. Semeraro F, Parrinello G, Cancarini A, Pasquini L, Zarra E, Cimino A, et al. Predicting the risk of diabetic retinopathy in type 2 diabetic patients. J Diabetes Complications (2011) 25:292–97. doi: 10.1016/j.jdiacomp.2010.12.002

PubMed Abstract | Crossref Full Text | Google Scholar

Keywords: referable diabetic retinopathy, nomogram construction, community screening, risk stratification, easily accessible indicators

Citation: Dongling N, Ziwei K, Juanling S, Li Z, Chang W, Ting L, Hongli L and Yanchun Z (2025) Universal nomogram for predicting referable diabetic retinopathy: a validated model for community and ophthalmic outpatient populations using easily accessible indicators. Front. Endocrinol. 16:1557166. doi: 10.3389/fendo.2025.1557166

Received: 08 January 2025; Accepted: 09 May 2025;
Published: 12 June 2025.

Edited by:

Tomislav Bulum, Medical School University of Zagreb, Croatia

Reviewed by:

Thiago Gonçalves dos Santos Martins, Federal University of São Paulo, Brazil
Hossam El-Din Moustafa, Mansoura University, Egypt

Copyright © 2025 Dongling, Ziwei, Juanling, Li, Chang, Ting, Hongli and Yanchun. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Zhang Yanchun, emhhbmd5YW5jaHVuMTIzOUBtZWQubnd1LmVkdS5jbg==; Liu Hongli, bGl1aG9uZ2xpQHhqdHUuZWR1LmNu

These authors have contributed equally to this work and shared first authorship

Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.