- 1Department of Endocrinology, Dalian Municipal Central Hospital, Dalian, China
- 2Infervision Institute of Research, Beijing, China
- 3Graduate School of Art and Science, Yale University, New Haven, CT, United States
Objective: To construct and validate prediction models for the risk of diabetic retinopathy (DR) in patients with type 2 diabetes mellitus.
Methods: Patients with type 2 diabetes mellitus hospitalized over the period between January 2010 and September 2018 were retrospectively collected. Eighteen baseline demographic and clinical characteristics were used as predictors to train five machine-learning models. The model that showed favorable predictive efficacy was evaluated at annual follow-ups. Multi-point data of the patients in the test set were utilized to further evaluate the model’s performance. We also assessed the relative prognostic importance of the selected risk factors for DR outcomes.
Results: Of 7943 collected patients, 1692 (21.30%) developed DR during follow-up. Among the five models, the XGBoost model achieved the highest predictive performance with an AUC, accuracy, sensitivity, and specificity of 0.803, 88.9%, 74.0%, and 81.1%, respectively. The XGBoost model’s AUCs in the different follow-up periods were 0.834 to 0.966. In addition to the classical risk factors of DR, serum uric acid (SUA), low-density lipoprotein cholesterol (LDL-C), total cholesterol (TC), estimated glomerular filtration rate (eGFR), and triglyceride (TG) were also identified to be important and strong predictors for the disease. Compared with the clinical diagnosis method of DR, the XGBoost model achieved an average of 2.895 years prior to the first diagnosis.
Conclusion: The proposed model achieved high performance in predicting the risk of DR among patients with type 2 diabetes mellitus at each time point. This study established the potential of the XGBoost model to facilitate clinicians in identifying high-risk patients and making type 2 diabetes management-related decisions.
Introduction
DR is the most common microvascular complication of diabetes mellitus. It has been demonstrated to be a leading cause of preventable blindness in the working-age population in most countries (1).
The American Academy of Ophthalmology (AAO), in 2019, stated that the prevalence of DR among diabetic patients worldwide is about 34.6%. 10.2% (28 million) diabetic patients suffer from vision-threatening DR (2). Effective management of DR requires a deep understanding of the predisposing factors, early diagnosis, and timely therapeutic intervention. Early identification of patients at risk of developing DR is the key to effective intervention, which is significant in reducing the progression of DR and thereby reducing the risk of blindness (3). Moreover, individual patients can be stratified according to different risk levels and get optimal treatment. Due to no typical symptoms in the early stage of the disease, however, most patients with DR may not seek medical evaluation until progression to the proliferative stage, resulting in irreversible visual damage (4). Therefore, methods for accurate prediction of the risk of DR are in urgent need.
At present, several DR risk prediction models based on cross-sectional studies have been developed (5–9). Deep learning algorithms were also applied (10). Although these models can predict the occurrence of DR at the index date, they cannot predict DR occurrence and development of the same patient at designated time points in the future. This will obviously restrict their clinical application. Similarly, based on the clinical characteristics related to the occurrence or development of DR, several models have been developed for optimization of the screening interval in DR screening (11, 12). However, due to the small number of cases in the studies, the proposed models have not been fully validated so far.
On the other hand, several studies investigated the pathogenesis and risk factors of DR to provide guidance for DR management. Epidemiological studies have shown that age, course of diabetes, hemoglobin A1c(HbA1c), fasting blood glucose (FBG), blood pressure, blood lipids, body mass index (BMI), smoking, proteinuria, and several others are all risk factors for DR (13, 14). Among them, duration of diabetes and hyperglycemia were demonstrated as strong risk factors for the occurrence and development of DR (15, 16). However, patients without DR were hardly unusual among those suffering from diabetes for a long time (17). The influences of other factors in DR occurrence also need to be proven. Further studies are required to elucidate the correlation and thus construct standard procedures for the management of this disease.
In this retrospective cohort study, we collected electronic health record data from hospitalized patients with type 2 diabetes and established the prediction model for future risks of DR based on a machine learning (ML) algorithm. To our knowledge, this is the first study to predict the occurrence of DR at each follow-up time point in up to 10 years. We also explore the risk factors that may affect the occurrence of DR and hope this work can provide a basis for further studies concerning the prevention and management of DR.
Materials And Methods
Study Subjects
The study protocol was approved by the Ethics Committee of Dalian Medical University Affiliated with the Central Hospital of Dalian. The ethics committee waived the requirement of written informed consent for all patients. The data did not contain any direct patient identifiers. The study adhered to the tenets of the Declaration of Helsinki.
Datasets
In this retrospective cohort study, we recruited inpatients admitted to the Department of endocrinology of Dalian Medical University Affiliated with the Central Hospital of Dalian from January 2010 to September 2018. Patients who met the following criteria were included in the current study: 1) Age≥18 years with the diagnosis of type 2 diabetes; 2) be hospitalized at least once over the follow-up period; 3) no presence of DR at the time of the first hospitalization. Patients with other types of diabetes (type 1 diabetes, special type diabetes, or gestational diabetes) or previous eye diseases (cataract, glaucoma, or other eye diseases) were excluded. With the criteria mentioned above, 7943 patients were selected for this study, including 1692 patients diagnosed with DR in the follow-up period (DR group) and 6521 control ones without DR (non-DR group). These patients were randomly divided into a training set (n=5559) and a test set (n=2384).
Diagnostic Criteria
Diagnostic criteria of type 2 DM are according to 1999 WHO diagnostic criteria (18). DR was examined by dilated fundus examination by ophthalmologists and endocrinologists together. Diagnostic criteria of DR case meet the Diabetic Retinopathy preferred practice pattern (PPP)-2019 guideline (19). The grading of DR in this study was based on the International Clinical Diabetic Retinopathy Severity Scales (20): Grace-1 for no apparent retinopathy; Grace-2 for mild non-proliferative diabetic retinopathy (NPDR), which includes microaneurysms only; Grace-3 for moderate NPDR, which includes more than just microaneurysms but less than severe nonproliferative DR; Grace-4 for severe NPDR (any of the following can be diagnosed as Grace-4: more than 20 intraretinal hemorrhages in each of 4 quadrants, definite venous beading in 2+ quadrants, Prominent intraretinal microvascular abnormalities in 1+ quadrant); and Grace-5 for proliferative diabetic retinopathy (PDR). The individual diagnoses were described according to the diagnosis and staging of the worse one between the two eyes.
Variable Selection
Patients’ information on demographics, medical history, medication, and laboratory data was extracted from the EHR at baseline and used as candidate predictor variables for developing DR (Tables S2, S3). These variables are commonly associated with the risk of DR, including 11 categorical variables and 20 continuous variables. The smoking status was defined as: 1) Current smokers: Patients who smoked more than 100 cigarettes and had not quit smoking at the index date; 2) Ex-smokers: Patients who smoked more than 100 cigarettes and stopped smoking 30 days before index date; 3) Never smokers: Patients who had never smoked. The drinking status was defined as: 1) Current drinkers: Patients who have been drinking or have not abstained from alcohol for more than one year; 2) Ex-drinkers: Patients who did not drink or had abstained from drinking for more than one year; 3) Never-drinks: Patients who had never drunk. Hypertension was defined as a systolic blood pressure of ≥140 mmHg and/or a diastolic blood pressure ≥90 mmHg, measured twice or more on different days, whether antihypertensive drugs were used or not (21). Blood samples were collected the next morning after hospital admission with at least 12-h fasting. HbA1c was measured using the glycosylated hemoglobin analyzer (TOSOH company of Japan, HLC-723G8). Four items of blood lipids (high-density lipoprotein cholesterol (HDL-C), TC, LDL-C, and TG), alanine transaminase (ALT), aspartate transaminase (AST), gamma-glutamyl transpeptidase(γ-GT), serum creatinine (Scr), SUA, and FBG were detected by automatic biochemical analyzer (Siemens, Germany, ADVIA2400 biochemical system). The eGFR was calculated according to the CKD-EPI formula (22). BMI was calculated as weight (kg) divided by the square of height in meters (m2); Blood pressure was measured on the right arm in sitting position three times consecutively at 5-min intervals, with the mean of the three measurements used for analysis.
Machine Learning Models Construction
Machine learning-based algorithms were selected according to the criteria below. First, the algorithms should show evaluation indicators based on mixed data type, including numerical variables and categorical variables. Second, the algorithms should have a wide range of applications and a history of successful usage in related fields. Based on the above criteria, we selected five machine learning algorithms, Random Forest (RF) (23), Extreme Gradient Boosting (XGBoost) (24), Logistic Regression (LR) (25), Support Vector Machine (SVM) (26), and K-Nearest Neighbor (K-NN) (27), to build the models with all the variables as predictors in this study. We used the GridSearchCV method to select the optimal key hyper-parameters of the five models as shown in Table S3, and other hyper-parameters were set as default. The models were then tested and internally validated by fivefold cross-validation.
Statistical Analysis
All statistical analyses were performed by using R statistical and computing software version 4.0.2 (http://www.r-project.org/). Chi-square test was used for categorical variables, and two-sample t-test for continuous variables to compare the variables between DR group and non-DR group. Two-tailed hypothesis testing was used. P<0.05 indicated a statistically significant result.
Performance of the proposed models was assessed mainly in terms of sensitivity, specificity, positive predictive value (PPV), negative predictive value (NPV), and accuracy. We assessed these metrics at the optimized prediction threshold identified by the Youden index. To compare performance across the models, we also calculated the area under the receiver operating characteristic curve (AUC). Significance testing and LASSO penalty were employed to identify clinical features that are associated with a high risk of DR. Nomograms were built based on the results of multivariate logistic regression analysis. The performance of the nomograms was measured using Harrell’s concordance index (C-index).
Results
Study Population
Baseline characteristics of 7943 patients with type 2 diabetes, including 1692 patients diagnosed as DR (DR group) and 6521 without DR signs (non-DR group) during the follow-up period, are presented in Table S4. The mean age of the patients was 63.54 years, and 45.7% were women. Mean follow-up time was 3.139 ± 2.243 years, and mean duration of diabetes mellitus was 1.449 ± 2.756 years. To further validate the performance of the model, 829 multipoint data from different follow-up time points were collected from 508 patients in the test set (supplementary test set).
Prediction Performances of the ML Models
The performances of the prediction models were assessed in terms of predefined evaluation metrics and ROC (Table 1 and Figure 1). The XGBoost model outperformed the other models with the highest AUC (0.913; 95% confidence interval (CI), 0.901-0.925), highest accuracy (79.9%; CI, 76.9%-83.5%), highest sensitivity (90.2%; CI, 84.8%-94.9%), highest specificity (77.1&; CI, 72.1%- 82.7%), highest PPV (51.6%; CI, 47.9%-57.5%) and highest NPV (96.7%; CI, 95.2%-98.1%).
Table 1 The performance metrics of the cross-validated machine learning algorithms on the test data.
Prediction Performance of the XGBoost Model in Different Follow-Up Periods
The XGBoost model was further assessed as a potential tool to predict the occurrence of DR (Table 2). According to follow-up intervals, the test set was divided into eight categories, respectively, each year of 1 to 8 years and more than 8 years. The highest AUC of 0.966 (95% CI, 0.955-0.977) and highest accuracy of 87.7% (95% CI, 85.5%-93.3%) were achieved in patients who were followed up within the timeframe of 1 to 2 years.
Information Contributions of the Features
The importance scores calculated by XGBoost are shown in Figure 2. Among all the features, the top five ones (HbA1c, duration, follow-up time, FBG, and age) contributed the most to identifying cases with DR risk in type 2 diabetes patients. In addition, the importance scores calculated by RF, which is the next-best performing model, are shown in Figure S1.
Subgroups Analysis Based on Follow-Up Period and DR Severity Level
The XGBoost model was further evaluated on the supplementary test set. Multipoint data of the same patient were used as new cases. The model showed an improved performance with an AUC, accuracy, sensitivity, specificity, PPV, and NPV of 0.922 (95%CI, 0.912-0.932), 81.5% (95% CI, 80.4%-84.2%), 92.6% (95% CI, 87.1%-94.2%), 76.6% (95% CI, 74.9%-82.4%), 63.6% (95% CI, 62.1%-68.8%), and 95.9% (95% CI, 93.5%-96.8%) respectively. Detailed attributes are presented in Table S5. Age, diuretic, fibrates, eGFR, HbA1c, duration, and follow-up time were significantly different between the cases with correct predictions and those with incorrect predictions. We then measured the time to diagnosis. The XGBoost model achieved an average of 2.895 ± 2.104 years prior to the first diagnosis by the clinical diagnosis method.
Subgroups analyses based on follow-up period and DR severity level were respectively conducted to assess the potential application of the model in clinical scenarios (Supplementary Figure S4). With the increase in follow-up time, the true positive prediction rate of the XGBoost model increases (Supplementary Figure S4A). The subgroup of patients followed up within the time frame of 7 to 8 years exhibited the highest positive prediction rate (100.0%). Data of all the others were higher than 87%. As for different DR severity subgroups, the model performed better for grade 3, reaching the highest positive prediction rate of 97.56% (Supplementary Figure S4B).
Nomograms for Predicting DR Risk of Patients With Type 2 Diabetes
Nomograms were developed to depict the association between the clinical variables and the occurrence of DR. After feature selection, we obtained 16 features including γGT, AST, SBP, eGFR, age, HbA1c, FBG, duration, follow-up time, insulin, diuretic, statins, fibrates, hypertension, smoking status, and drinking status. A classic LR model and an LR model integrated with the machine learning output were constructed (Supplementary Figure S2, S3). The C-index of the integrated LR model was 0.921, whereas that of the classic LR model was 0.836.
Discussion
In this study, we constructed five ML-based models for the risk of developing DR among patients with type 2 diabetes using EHR data, the XGboost model of which demonstrated the best performance. We further evaluated the predictive value of this model and found that it can accurately predict the occurrence of DR at each time point over 10 years. According to feature importance analysis, five risk factors, SUA, LDL-C, TC, eGFR, and TG, were recognized as important indicators of DR for the first time. Moreover, analysis of data from multiple time points revealed that predicting DR risk with the XGBoost model would obviously reduce the amount of time needed for the diagnosis of DR.
XGBoost algorithm is known for its good scalability and high running speed, which is thus commonly used and developed in recent years (28–30). Specifically, the XGBoost model provides the following advantages. First, the XGBoost algorithm reduces the order of magnitude of features while retaining all the features by introducing regularization terms into the objective function to avoid over-fitting. Secondly, the algorithm also used the random forest algorithm for reference, which can not only reduce the over-fitting problem but also reduce the computational complexity. We hypothesized that these reasons would be responsible for the XGBoost algorithm’s outstanding prediction ability for DR risk compared with the other four methods adopted in this study.
There were several studies gaining insights into DR risk prediction, some of which developed ML-based models. However, most of them were cross-sectional studies. In a study in Iran (31), 3734 diabetic patients were included to build a logistic regression model with an AUC of 0.704. The input variables included sex, age, diabetes duration, BMI, blood pressure (BP), HbA1c, FPS, cholesterol, and triglycerides. In the cross-sectional study by Ein hoh et al. (7) in South Korea, a LASSO prediction model was constructed for predicting the risk of DR, reaching an AUC of 0.82 and an accuracy of 75.2%. 490 diabetic patients were selected for the study, and several variables (age, sex, smoking history, drinking history, waist circumference, BMI, physical exercise status, medication history, blood pressure, and relevant laboratory results) were used as input. In another cross-sectional study in Taiwan Province of China (6), 536 patients with type 2 diabetes were selected and 10 predictors, including systolic blood pressure, diastolic blood pressure, BMI, age, gender, duration of diabetes, family history of diabetes, whether there was a blood glucose test, whether there was exercise, and whether there was insulin treatment, were included. Then, four prediction models (SVM, DT, artificial neural network, and LR) were constructed, respectively. The model based on SVM showed the best performance with an accuracy of 79.5% and an AUC of 0.839. Mo et al. (5). included 4170 diabetic patients in a cross-sectional study to predict the risk of DR. Seven input variables, including age, diabetes duration, postprandial blood glucose, HbA1c, urine creatinine, urine microalbumin, and systolic blood pressure were used to construct a multivariate logistic regression model. The model showed moderate predictive ability with an AUC of 0.715 in the validation set. The predictive models mentioned above were developed to identify patients with DR rather than predict the occurrence of DR in the future. Therefore, they can hardly be adapted to the scenarios that require pre-diagnosis screening and prompt intervention. To our knowledge, only one follow-up study was performed for predicting the occurrence of DR (32). The study was conducted with 5034 type 2 diabetes patients not affected by retinopathy at the time of the recruitment and a median follow-up time of 1.2 years. A Cox risk prediction model was constructed with seven predictors including diabetes duration, HbA1c, systolic blood pressure, diastolic blood pressure, proteinuria, creatinine clearance, and diabetes drug treatment, and showed a helpful predictive ability (C-index=0.746). However, this model cannot meet all the clinical requirements of long-term management of type 2 diabetes and associated complications considering that it can only predict the occurrence of DR in the next 1 to 4 years.
The XGBoost model we constructed performed well in different follow-up periods (AUC=0.834-0.966), which indicates that it could predict whether and when a patient with type 2 diabetes will develop DR in the next 1 to 8 years or even over 8 years. Results also show the proposed model is effective to demonstrate in patients who had DR of different severity levels during follow-up. For the individuals with moderate NPDR, the model exhibited the highest performance. As a transitional stage, moderate NPDR can progress and advance to vision-threatening retinopathies. Intervention methods such as intravitreous anti-VEGF treatment can be applied to secure moderate NPDR patients. As the patient progresses to severe NPDR and PDR, it is hard to get the same treatment effect. Moreover, the method requires no extra laboratory tests since the model input was drawn from demographic and clinical characteristics and routine test results. Therefore, the model will facilitate clinicians risk-stratifying patients with type 2 diabetes. For patients at risk of DR, clinicians can provide adequate follow-up management and effective therapeutic intervention. While for the low-risk population, visit frequency will be appropriately reduced, and thus personal and societal healthcare burdens will be reduced accordingly. We further developed nomograms to depict the association between the clinical variables and the probabilities of DR. The nomogram integrating the XGBoost-based machine learning output achieved a higher predictive performance, which provides an intuitive way to interpret the model and shows its potential to be a clinical decision support tool.
To identify the key features contributing to the pre-diagnosis of DR, the importance scores of the risk factors were calculated by XGBoost. According to the feature importance analysis, HbA1c, diabetes duration, age, and FBG were highly ranked, which corresponds to the current clinical perceptions. In addition, we found that SUA, LDL-C, TC, eGFR, and TG may have the potential to be strong predictors for DR. Moreover, we further calculated the importance scores of the risk factors by RF, and found that the results were similar to and complementary to those from XGBoost. In previous studies, Krizova et al. found that an elevated level of vitreous uric acid was closely related to the development of DR, and uric acid concentration was closely related to the degree of DR progression (33). Zhu et al. found that SUA in the vitreous activates the expression of retinal inflammatory factors through the Notch signaling pathway, which will induce oxidative stress and inflammatory response and thus promote the occurrence and development of DR (34). Obesity is a metabolic disease that is associated with insulin resistance and diabetes mellitus (35). In the Wisconsin Epidemiologic Study of Diabetic Retinopathy (WESDR), it was found that high TC levels could increase the risk of retinal hard exudation in patients with type 2 diabetes (36). High LDL-C and high TC/LDL-C levels were associated with retinal macular edema in SN-DREAMS studies (37). Besides, statin lipid-lowering drugs could significantly reduce the risk of DR (38). Since DR is a microvascular complication of diabetes, eGFR can be used as an important biochemical indicator reflecting DR. Chen et al. found that higher eGFR was positively correlated with the risk of DR (39). Therefore, the results presented in our study indicated that these risk factors could be used as early and effective predictors of DR for type 2 diabetes patients before recognizable symptoms appear.
There are still several limitations to this study. First, the prediction models were constructed with the risk factors extracted from EHR, except for those with excessive missing data including the ones that might be related to the occurrence of DR (e.g., hip circumference and neck circumference). Secondly, individualized management approaches for type 2 diabetes and associated complications were not fully described. The monitoring and tracking features were also not collected during the follow-up period. The above information will be recorded and used to optimize the XGBoost model in the future. In addition, this is a single-center retrospective study. We are planning to collect cases from several institutions and conduct a prospective study to validate the generalization ability of the proposed model.
In summary, we developed and evaluated an ML-based model for predicting the risk of DR. The results suggest that it is possible to pre-diagnose DR without fundus images. We also demonstrated the potential of applying XGBoost models to facilitate clinicians in accurately identifying the high-risk population for DR and formulating patient-specific management strategies, thereby reducing the occurrence and development of DR.
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 human participants were reviewed and approved by Ethics Committee of Dalian Central Hospital. The ethics committee waived the requirement of written informed consent for participation.
Author Contributions
XHL, ZG, and XYL conceived of and designed the study. YZ contributed to data collection and wrote the manuscript. MD contributed to the conceptualization of the project and edited the manuscript. SL contributed to discussion and revision of the manuscript. HY contributed to model building and wrote the corresponding manuscript. WC contributed to the modeling scheme design and technical support. MZ, PL, and QY contributed to part of manuscript proofreading. XHL, ZG, and XYL are the guarantor of this project and had full access to all of the data in the study and accuracy of the data analysis and reviewed the manuscript. All authors contributed to the article and approved the submitted version.
Funding
This work was funded by the National key R&D Program of China (2016YFC0901200), the Dalian Medical Science Research Project Grant Support (2011007).
Conflict of Interest
MD and WC are employed by Infervision.
The remaining 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.
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.
Acknowledgments
We would like to thank Yidu Cloud (Beijing) Technology Co, Ltd for their assistance in data searching, extraction and processing.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fendo.2022.876559/full#supplementary-material
References
1. Kempen JH, O'Colmain BJ, Leske MC, Haffner SM, Klein R, Moss SE. The Prevalence of Diabetic Retinopathy Among Adults in the United States. Arch Ophthalmol (2004) 122(4):552–63. doi: 10.1001/archopht.122.4.552
2. Flaxel CJ, Adelman RA, Bailey ST, Fawzi A, Lim JI, Vemulakonda GA. Diabetic Retinopathy Preferred Practice Pattern®. Ophthalmology (2020) 127(1):66–145. doi: 10.1016/j.ophtha.2019.09.025
3. Liew G, Michaelides M, Bunce C. A Comparison of the Causes of Blindness Certifications in England and Wales in Working Age Adults (16-64 Years), 1999-2000 With 2009-2010. BMJ Open (2014) 4(2):e004015. doi: 10.1136/bmjopen-2013-004015
4. Jampol LM, Glassman AR, Sun J. Evaluation and Care of Patients With Diabetic Retinopathy. N Engl J Med; (2020) 382(17):1629–37. doi: 10.1056/NEJMra1909637
5. Mo R, Shi R, Hu Y, 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
6. Tsao H-Y, Chan P-Y, Su EC-Y. Predicting Diabetic Retinopathy and Identifying Interpretable Biomedical Features Using Machine Learning Algorithms. BMC Bioinf (2018) 19(Suppl 9):283. doi: 10.1186/s12859-018-2277-0
7. Oh E, Yoo TK, Park E-C. Diabetic Retinopathy Risk Prediction for Fundus Examination Using Sparse Learning: A Cross-Sectional Study. BMC Med Inf decision making (2013) 13:106. doi: 10.1186/1472-6947-13-106
8. Ogunyemi O, Kermah D. Machine Learning Approaches for Detecting Diabetic Retinopathy From Clinical and Public Health Records. AMIA Annu Symp Proc (2015) 2015:983–90.
9. Ogunyemi OI, Gandhi M, Tayek C. Predictive Models for Diabetic Retinopathy From Non-Image Teleretinal Screening Data. AMIA Summits Trans Sci Proc (2019) 2019:472–7.
10. Bora A, Balasubramanian S, Babenko B, Virmani S, Venugopalan S, Mitani A. Predicting the Risk of Developing Diabetic Retinopathy Using Deep Learning. Lancet Digit Health (2021) 3(1):e10–9. doi: 10.1016/S2589-7500(20)30250-8
11. Mehlsen J, Erlandsen M, Poulsen PL, Bek T. Individualized Optimization of the Screening Interval for Diabetic Retinopathy: A New Model. Acta Ophthalmol (2012) 90(2):109–14. doi: 10.1111/j.1755-3768.2010.01882.x
12. Dagliati A, Marini S, Sacchi L, Cogni G, Teliti M, Tibollo V. Machine Learning Methods to Predict Diabetes Complications. J Diabetes Sci Technol (2018) 12(2):295–302. doi: 10.1177/1932296817706375
13. Ting DSW, Cheung GCM, Wong TY. Diabetic Retinopathy: Global Prevalence, Major Risk Factors, Screening Practices and Public Health Challenges: A Review. Clin Exp Ophthalmol (2016) 44(4):260–77. doi: 10.1111/ceo.12696
14. Varma R, Macias GL, Torres M, Klein R, Peña FY, Azen SP. Biologic Risk Factors Associated With Diabetic Retinopathy: The Los Angeles Latino Eye Study. Ophthalmology (2007) 114(7):1332–40. doi: 10.1016/j.ophtha.2006.10.023
15. Nathan DM, Genuth S, Lachin J, Cleary P, Crofford O, Davis M. The Effect of Intensive Treatment of Diabetes on the Development and Progression of Long-Term Complications in Insulin-Dependent Diabetes Mellitus. N Engl J Med (1993) 329(14):977–86. doi: 10.1056/NEJM1993093032914
16. UK Prospective Diabetes Study (UKPDS) Group. Intensive Blood-Glucose Control With Sulphonylureas or Insulin Compared With Conventional Treatment and Risk of Complications in Patients With Type 2 Diabetes (UKPDS 33). Lancet (1998) 352(9131):837–53. doi: 10.1016/S0140-6736(98)07019-6
17. Keenan HA, Costacou T, Sun JK, Doria A, Cavellerano J, Coney J. Clinical Factors Associated With Resistance to Microvascular Complications in Diabetic Patients of Extreme Disease Duration: The 50-Year Medalist Study. Diabetes Care (2007) 30(8):1995–7. doi: 10.2337/dc06-2222
18. Alberti KG, Zimmet PZ. Definition, Diagnosis and Classification of Diabetes Mellitus and its Complications. Part 1: Diagnosis and Classification of Diabetes Mellitus Provisional Report of a WHO Consultation. Diabetes Med (1998) 15(7):539–53. doi: 10.1002/(SICI)1096-9136(199807)15:7<539::AID-DIA668>3.0.CO;2-S
19. Niroomand M, Afsar J, Hosseinpanah F, Afrakhteh M, Farzaneh F, Serahati S. Comparison of the International Association of Diabetes in Pregnancy Study Group Criteria With the Old American Diabetes Association Criteria for Diagnosis of Gestational Diabetes Mellitus. Int J Endocrinol Metab (2019) 17(4):e88343. doi: 10.5812/ijem.88343
20. Wilkinson CP, Ferris FL 3rd, Klein RE, Lee PP, Agardh CD, Davis M. Proposed International Clinical Diabetic Retinopathy and Diabetic Macular Edema Disease Severity Scales. Ophthalmology (2003) 110(9):1677–82. doi: 10.1016/S0161-6420(03)00475-5
21. Williams B, Mancia G, Spiering W, Agabiti Rosei E, Azizi M, Burnier M. 2018 ESC/ESH Guidelines for the Management of Arterial Hypertension’. Eur Heart J (2018) 39(33):3021–104. doi: 10.1093/eurheartj/ehy339
22. Levey AS, Stevens LA, Schmid CH, Zhang Y, Castro AF III, Feldman HI. A New Equation to Estimate Glomerular Filtration Rate. Ann Internal Med (2009) 150(9):604–12. doi: 10.7326/0003-4819-150-9-200905050-00006
24. Chen T, Guestrin C. XGBoost: A Scalable Tree Boosting System. KDD '16: Proceedings of the 22nd ACM SIGKDD International Conference on Knowledge Discovery and Data Mining. (2016), 785–94. doi: 10.1145/2939672.2939785
26. Cortes C, Vapnik V. Support-Vector Networks. Mach Learn (1995) 20(3):273–97. doi: 10.1007/BF00994018
27. Cover T, Hart P. Nearest Neighbor Pattern Classification. IEEE Trans Inf Theory January (1967) 13(1):21–7. doi: 10.1109/TIT.1967.1053964
28. Hou N, Li M, He L, Xie B, Wang L, Zhang R. Predicting 30-Days Mortality for MIMIC-III Patients With Sepsis-3: A Machine Learning Approach Using XGboost. J Transl Med (2020) 18(1):462. doi: 10.1186/s12967-020-02620-5
29. Ogunleye A, Wang QG. XGBoost Model for Chronic Kidney Disease Diagnosis. IEEE/ACM Trans Comput Biol Bioinform. (2020) 17(6):2131–40. doi: 10.1109/TCBB.2019.2911071
30. Zhang Y, Feng T, Wang S, Dong R, Yang J, Su J. A Novel XGBoost Method to Identify Cancer Tissue-Of-Origin Based on Copy Number Variations. Front Genet (2020) 11:585029. doi: 10.3389/fgene.2020.585029
31. Hosseini SM, Maracy MR, Amini M, Baradaran HR. A Risk Score Development for Diabetic Retinopathy Screening in Isfahan-Iran. J Res Med Sci (2009) 14(2):105–10.
32. Semeraro F, Parrinello G, Cancarini A, Pasquini L, Zarra E, Cimino A. Predicting the Risk of Diabetic Retinopathy in Type 2 Diabetic Patients. J Diabetes its Complications (2011) 25(5):292–7. doi: 10.1016/j.jdiacomp.2010.12.002
33. Krizova L, Kalousova M, Kubena A, Benakova H, Zima T, Kovarik Z. Increased Uric Acid and Glucose Concentrations in Vitreous and Serum of Patients With Diabetic Macular Oedema. Ophthalmic Res (2011) 46(2):73–9. doi: 10.1159/000322994
34. Zhu D-D, Wang Y-Z, Zou C, She X-P, Zheng Z. The Role of Uric Acid in the Pathogenesis of Diabetic Retinopathy Based on Notch Pathway’. Biochem Biophys Res Commun (2018) 503(2):921–9. doi: 10.1016/j.bbrc.2018.06.097
35. Li C, Qu L, Farragher C, Vella A, Zhou B. MicroRNA Regulated Macrophage Activation in Obesity. J Trans Internal Med (2019) 7(2):46. doi: 10.2478/jtim-2019-0011
36. LeCaire TJ, Palta M, Klein R, Klein BE, Cruickshanks KJ. Assessing Progress in Retinopathy Outcomes in Type 1 Diabetes: Comparing Findings From the Wisconsin Diabetes Registry Study and the Wisconsin Epidemiologic Study of Diabetic Retinopathy. Diabetes Care (2013) 36(3):631–7. doi: 10.2337/dc12-0863
37. Raman R, Ganesan S, Pal SS, Gella LL, Kulothungan V, Sharma T. Incidence and Progression of Diabetic Retinopathy in Urban India: Sankara Nethralaya-Diabetic Retinopathy Epidemiology and Molecular Genetics Study (SN-DREAMS II), Report 1. Ophthalmic Epidemiol (2017) 24(5):294–302. doi: 10.1080/09286586.2017.129025
38. Kang EY-C, Chen T-H, Garg SJ, Sun C-C, Kang J-H, Wu W-C. Association of Statin Therapy With Prevention of Vision-Threatening Diabetic Retinopathy. JAMA Ophthalmol (2019) 137(4):363–71. doi: 10.1001/jamaophthalmol.2018.6399
Keywords: type 2 diabetes mellitus, diabetic retinopathy, machine learning, XGBoost model, cohort study
Citation: Zhao Y, Li X, Li S, Dong M, Yu H, Zhang M, Chen W, Li P, Yu Q, Liu X and Gao Z (2022) Using Machine Learning Techniques to Develop Risk Prediction Models for the Risk of Incident Diabetic Retinopathy Among Patients With Type 2 Diabetes Mellitus: A Cohort Study. Front. Endocrinol. 13:876559. doi: 10.3389/fendo.2022.876559
Received: 15 February 2022; Accepted: 12 April 2022;
Published: 17 May 2022.
Edited by:
Subrata Chakrabarti, Western University, CanadaReviewed by:
Matthew Cecchini, London Health Sciences Centre, CanadaChristina A. Castellani, Western University, Canada
Copyright © 2022 Zhao, Li, Li, Dong, Yu, Zhang, Chen, Li, Yu, Liu and Gao. 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: Xuhan Liu, TGl1eHVoYW4yMDIyMDIxNEAxNjMuY29t; Zhengnan Gao, Z2FvMjAwOEAxNjMuY29t
†These authors have contributed equally to this work and share first authorship