这是indexloc提供的服务,不要输入任何密码
Skip to main content

Deep-learning prediction of cardiovascular outcomes from routine retinal images in individuals with type 2 diabetes

Abstract

Background

Prior studies have demonstrated an association between retinal vascular features and cardiovascular disease (CVD), however most studies have only evaluated a few simple parameters at a time. Our aim was to determine whether a deep-learning artificial intelligence (AI) model could be used to predict CVD outcomes from routinely obtained diabetic retinal screening photographs and to compare its performance to a traditional clinical CVD risk score.

Methods

We included 6127 individuals with type 2 diabetes without myocardial infarction or stroke prior to study entry. The cohort was divided into training (70%), validation (10%) and testing (20%) cohorts. Clinical 10-year CVD risk was calculated using the pooled cohort equation (PCE) risk score. A polygenic risk score (PRS) for coronary heart disease was also obtained. Retinal images were analysed using an EfficientNet-B2 network to predict 10-year CVD risk. The primary outcome was time to first major adverse CV event (MACE) including CV death, myocardial infarction or stroke.

Results

1241 individuals were included in the test cohort (mean PCE 10-year CVD risk 35%). There was a strong correlation between retinal predicted CVD risk and the PCE risk score (r = 0.66) but not the polygenic risk score (r = 0.05). There were 288 MACE events. Higher retina-predicted risk was significantly associated with increased 10-year risk of MACE (HR 1.05 per 1% increase; 95% CI 1.04–1.06, p < 0.001) and remained so after adjustment for the PCE and polygenic risk score (HR 1.03; 95% CI 1.02–1.04, p < 0.001). The retinal risk score had similar performance to the PCE (both AUC 0.697) and when combined with the PCE and polygenic risk score had significantly improved performance compared to the PCE alone (AUC 0.728). An increase in retinal-predicted risk within 3 years was associated with subsequent increased MACE likelihood.

Conclusions

A deep-learning AI model could accurately predict MACE from routine retinal screening photographs with a comparable performance to traditional clinical risk assessment in a diabetic cohort. Combining the AI-derived retinal risk prediction with a coronary heart disease polygenic risk score improved risk prediction. AI retinal assessment might allow a one-stop CVD risk assessment at routine retinal screening.

Graphical abstract

Introduction

Cardiovascular disease (CVD) is the leading cause of morbidity and mortality among individuals with type 2 diabetes mellitus [1]. Diabetes is associated with accelerated incidence and progression of CVD [2] and CVD prevention plays a key role in diabetes management [3]. Recent trials have demonstrated that sodium-glucose co-transporter 2 inhibitors (SGLT2i) [4] and glucagon-like peptide 1 receptor agonists (GLP1-RA) [5] reduce the incidence of major cardiovascular events in individuals with type 2 diabetes at increased CVD risk, hence opportunities for improved CVD risk prediction in these patients might help to select patients who would derive particular benefit from these treatments.

In many healthcare systems, including in the United Kingdom, individuals with diabetes undergo routine retinal screening to identify the presence of diabetic retinopathy or maculopathy [6]. Such screening photographs allow direct visualization of the retinal microvasculature that could be used to inform on CVD risk. Several studies have identified the presence of retinopathy as an independent CVD risk predictor [7,8,9,10]. Subsequent studies have found associations between quantitative retinal vascular parameters such as vessel diameter, tortuosity and fractal dimension and increased CVD risk [11,12,13,14,15,16]. Measurement of these parameters can however be time-consuming and is not routinely performed in clinical practice.

Retinal photographs are ideal for analysis by artificial intelligence (AI) methodology [17]. Use of AI can not only reduce analysis time, but also provides the opportunity to take full advantage of the data contained within the image as opposed to just one or two parameters. Recent work has identified the ability of AI to predict cardiovascular risk factors and outcomes from retinal images with high diagnostic accuracy [18,19,20,21]. There have been few cohort studies reporting the feasibility of using AI to determine an individual’s CVD risk from a retinal photograph. The retinal microvasculature, directly observable with comparatively inexpensive devices, might not only provide an assessment of environmental risk factors associated with CVD, but it may also be reflective of an individual’s genetic susceptibility to these factors. The concept of using a retinal photo, perhaps combined with a marker of genetic CVD risk, as a one-stop assessment of global CVD risk is certainly appealing [11]. Little work exists on whether retinal images from multiple visits could be used to provide a dynamic assessment of CVD risk.

The aim of this study was to use a deep-learning AI model to predict the incidence of major adverse cardiovascular events (MACE) from routinely obtained diabetes retinal screening (DRS) photographs in individuals without prior myocardial infarction or stroke, determine the model’s prognostic performance when combined with genetic risk scores and ascertain whether changes in retinal AI-predicted CVD risk could be used to predict outcome.

Methods

Study cohort

We evaluated patients from the Genetics of Diabetes Audit and Research in Tayside Scotland (GoDARTS) study. Full cohort details have been previously described [22]. GoDARTS is a case-control cohort study conducted in Tayside, Scotland and includes 10,149 individuals with type 2 diabetes and 8,157 individuals without type 2 diabetes at the time of recruitment. A comprehensive set of baseline data were obtained at recruitment including clinical and lifestyle parameters, and consent was given for electronic health record linkage for past and future clinical events and healthcare contacts such as prescribing, laboratory tests, eye screening, hospitalisations and deaths. A blood sample was obtained at the time of recruitment for genotyping, with genome-wide association performed using various genotyping arrays as previously described. The GoDARTS study and electronic health record (EHR) linkage has been approved by the East of Scotland Research Ethics Committee. The EHR is fully anonymized and provided to researchers through robust information governance protocols administered by the Health Informatics Centre (HIC) Trusted Research Environment (TRE), including overarching research ethics approval for studies conducted within the TRE.

We studied patients with type 2 diabetes with no history of myocardial infarction (MI) or stroke (ICD codes I21-I23 and I60-I63) prior to the date of the first available retinal photograph (which was taken as the study entry date for this analysis). This date was also used to determine baseline demographic and clinical characteristics in order to calculate clinical 10-year cardiovascular risk using the Pooled Cohort Equation [23]. Clinical variables nearest the study entry date were used. Where there was more than one measure of a variable (e.g. cholesterol, blood pressure) we used the median over the 3 years prior to study entry. A genome wide polygenic risk score (PRS) for coronary heart disease (CHD) was constructed from the genome-wide genotyping data available in the GoDARTS bioresource using published data as described previously by Khera et al. [11, 24]. This PRS was taken from a GWAS involving over 180,000 individuals and included over 6.6 million genetic variants. The PRS was z-standardised.

Retinal imaging

Retinal photography was captured during routine annual diabetic retinal screening. The Scottish Diabetic Retinal Screening programme uses standardized protocols (details available at https://www.ndrs.scot.nhs.uk/). In brief, all the retinal images have 45° field of view and are macula-centred. All photographs had been previously assessed by a trained ophthalmologist for the presence of retinopathy or maculopathy for clinical purposes at the time of photography. We used the earliest available retinal photograph to train and test the DL model, with images from both eyes used.

Image pre-processing and data split

The DRS images used had 13 different image sizes, three of which formed the vast majority (98.9%) of the data: 2336 × 3504 pixels (92.7%), 2304 × 3456 (4%) and 1696 × 2544 (2.2%). Smaller dimensions accounted for only 1% of the images. There were also significant variations in these images in terms of luminosity, color, focus and quality. To reduce the effects of these variations in training a DL model we pre-processed the images in three steps as described previously [25, 26]. First, we fit a circle on the original retinal image to remove the non-retinal (black background) region. Second, we resized the image to a standard size (260 × 260) and performed contrast limited adaptive histogram equalization (CLAHE) on R, G, and B color channels separately. Finally, pixel intensities were normalized to [0, 1]. A sample retinal image before and after pre-processing is shown in Fig. 1A and B.

Fig. 1
figure 1

Example of pre-processing result (top) and Average of Grad-CAM heatmaps for PCE risk score prediction superimposed on a sample left eye and right eye retina (bottom). A Original retinal image; B pre-processed image; C heatmap for left eye retinas; D heatmap for right eye retinas

The cohort was randomly split into three groups: 70% training, 10% validation and 20% testing. To avoid information leak, all retinal images belonging to same individual were allocated to one group only.

Deep learning architecture and training

A deep convolutional neural network (CNN) was used to predict risk of MACE from the retinal images. We adopted EfficientNet-B2 [27], part of the EfficientNet family of architectures, a well-tested model that has been reported to achieve excellent performance in image classification tasks [28]. The EfficientNet-B2 architecture was modified by replacing its fully connected layer with a global average pooling layer followed by a single output node with sigmoid activation. The convolutional layers were unchanged. We used 260 × 260 input images following recommendations for best performance suggested by the authors of the EfficientNet-B2 [27]. As transfer learning, we initialized the weights of all the convolutional layers with pre-trained weights on ImageNet [28], allowing us to avoid the expensive task of training the network from scratch [29]. The model contained a total of 7.7 M trainable parameters and all of these were trained.

The training split was used to train the model parameters, with the validation split used to avoid overfitting of the model on the trained dataset. Binary cross-entropy was used as the loss function. We employed the Nadam optimizer from Keras (https://keras.io/api/optimizers/Nadam/) with an initial learning rate of 0.001, reduced by a factor of 0.1 if the validation loss did not improve within 5 consecutive epochs (https://keras.io/api/callbacks/reduce_lr_on_plateau/). The minimum learning rate was set to 10− 5. The model was trained with a batch size 32 for 50 epochs, a number we empirically observed to be sufficient for convergence due to transfer learning. A random horizontal flip and random rotation were used for image augmentation of the training images. During validation, early stopping (https://keras.io/api/callbacks/early_stopping/) was implemented, halting training if there was no improvement in validation loss for 20 consecutive epochs. The weights set with best validation performance were saved for the final testing phase.

All work was carried out in the University of Dundee HIC Trusted Research Environment [30] on a NVIDIA TITAN Xp GPU. Python 3.6 was used for code development with libraries opencv [31], scikit-learn [32] for image processing, and Keras 2.2.2[33] with tensorflow 1.9.0 as back-end for building [34], training and testing DL models.

Heatmap visualization

We used gradient-based class activation mapping (grad-CAM) [35] to visualize the regions in the input images that contained critical information for MACE prediction. We applied grad-CAM to the last convolutional layer of the modified EfficientNet-B2.

Prediction of 10-year MACE

The ML model was used to calculate the predicted 10-year MACE (predMACE10) from the date of the retinal image in the test data. Since MACE events are linked to each individual (rather than specifically to the left or right eye), we included only individuals with retinal images from both eyes in the test dataset. Following previous work that demonstrated improved predictive accuracy by averaging results from multiple left and right eye retinal images of the same individual [36], we calculated an individual-level MACE prediction by averaging the predMACE10 scores from both the left and right eye images.

As well as prediction of MACE from baseline retinal images we also performed a longitudinal analysis to evaluate changes in retinal predicted MACE risk over time. We evaluated retinal images from all individuals in the test dataset who had a further retinal photograph available for analysis within 3 years of the baseline image. We calculated the difference between the predMACE10 score from the baseline retinal image to the subsequent image. To account for the second image being at different timepoints within the 3-year window we calculated the rate of change as follows:

$$ \begin{gathered} \quad Predicted~MACE_{{10}} ~difference \hfill \\ = ~\frac{{predMACE_{{10}} ~from~last~available~image - predMACE_{{10}} ~from~first~available~image}}{{Duration~between~last~and~first~image~capture~dates}} \hfill \\ \end{gathered} $$

Thus positive values suggested an increased MACE risk, while negative values indicated a reduced MACE risk.

Clinical outcomes

Clinical events from the date of study entry to December 2018 were obtained from healthcare electronic medical records from the GoDARTS cohort as previously described [11]. The primary outcome for this analysis was time to first MACE, defined as cardiovascular death, non-fatal myocardial infarction or non-fatal stroke using ICD codes. Causes of death were obtained from the General Register of Scotland, with any ICD-10 code from I00-I99 within the first two causes of death being classified as a cardiovascular death. Follow-up was limited to 10 years to correspond to the duration of risk prediction from the PCE risk score.

Statistical analysis

Continuous variables are reported here as mean ± standard deviation or median and interquartile range as appropriate and categorical variables as number and percentage.

Following recent common practice [18, 36,37,38], we used the area under receiver operating characteristic curve (AUC; ROC) for 10-year MACE as our primary evaluation to analyse the model’s performance on test data, defined in detail elsewhere [25]. To assess the statistical significance of the model performance, non-parametric bootstrap sampling was used on test data. We used 2000 random samples with replacement from the test data, where each sample size was the same as that of the test data. We computed evaluation metrics from each bootstrap sample and 95% confidence intervals (CI) was reported from the distribution obtained from all bootstrap samples.

Associations between predMACE10, PCE and PRS were assessed using Pearson’s correlation coefficient. To evaluate the association between the predMACE10 and clinical outcomes, we used Cox proportional hazard regression analysis and the Kaplan-Meier (KM) estimator used on tertiles of predMACE10 from the test data. To quantify the association further, we performed with multivariable adjustment for clinical and genetic cardiovascular risk calculated by the PCE and the coronary heart disease polygenic risk score.

To determine the comparative prognostic performance of predMACE10 we calculated the area under the curve for the association with MACE at 10 years. To determine the association between the predMACE10 change and clinical outcomes, the KM estimator was used on the top 20% and the bottom 80% of the computed predMACE10 difference from the test data, with Cox regression used to quantify the association with 10-year MACE.

Survival analysis was carried out using R (3.2.5) with ‘survival’ package and ‘survminer’ package (for the plots). All tests were two-sided and a p value < 0.05 considered significant.

Results

Baseline characteristics

In total, of the 10,149 individuals with type 2 diabetes in GoDARTS we included 6,127 individuals who had available retinal photographs (total 12,859 retinal photographs), divided into ~ 70% training (4281 individuals with 8990 images), ~ 10% validation (605 individuals with 1271 images) and ~ 20% testing (1241 individuals with 2598 images). Full baseline characteristics are described in Table 1. There were no significant differences in the data distribution between the training, validation and test cohorts. In the overall cohort, the mean age at baseline was 67 ± 11 years, with 55% of the cohort being male. As would be expected from a diabetes cohort, clinical cardiovascular risk was high, with a mean PCE 10-year cardiovascular disease risk score of 34 ± 20%.

Table 1 Baseline characteristics of whole cohort and data splits

Cardiovascular disease risk prediction using deep learning

For training, validation and testing the EfficientNet-B2 model, retinal images of the left and right eyes at baseline (earliest available image) were considered.

There was no significant difference in model performance between males and females. For prediction of 10-year MACE the model achieved an AUC of 0.691 (95% CI 0.666–0.715) on the complete test dataset; 0.688 (0.657–0.719) using left eye images only; 0.689 (0.660–0.718) using right eye images only and 0.697 (0.663–0.731) using the mean of both eyes. Given the similar performance using all four methods we used the individual level (i.e. the mean of both eyes for each individual) predMACE10 for further analyses.

The grad-CAM heatmaps indicated that the optic disc, macula and vasculature were the most important retinal features for predicting 10-year MACE. Figure 1C and D show average of grad-CAM heatmaps computed for all left eye and right eye retinal images for CVD risk prediction in the test dataset, superimposed on a sample retinal image to illustrate the anatomical correspondence for highlighted areas.

Association between retinal predicted MACE, clinical CVD risk and CHD PRS

Retinal predMACE10 appeared to correlate predominantly with clinical rather than genetic CVD risk. There was a strong correlation between retinal predMACE10 and actual PCE CVD risk (r = 0.66). There was no correlation between retinal predMACE10 and the coronary heart disease PRS (r = 0.05). Similarly, there was little correlation between the PRS and the actual PCE CVD risk (r=-0.02) (Fig. 2).

Fig. 2
figure 2

Correlations between individual-level retinal predicted 10-year MACE, clinical CVD risk and genetic CHD risk. A Predicted 10-year MACE vs. CHD PRS. B Predicted 10-year MACE vs. PCE risk score. Red dots represent females and blue dots represent males

Retinal predicted MACE and clinical outcome

In the test cohort, over a median follow-up duration of 9.1 years, there were 288 MACE events, including 244 CV deaths, 71 non-fatal MIs and 50 non-fatal strokes within 10 years from the date of retinal imaging. Higher retinal predMACE10 was significantly associated with increased 10-year risk of MACE (HR 1.05 per 1% increase; 95% CI 1.04–1.06, p < 0.001). Individuals in the highest tertile of retinal predMACE10 (predicted 10-year MACE 28–59%) had a significantly higher risk of MACE than those in the lower (2–15%) and middle (15–28%) tertiles (compared to lowest tertile - middle tertile HR 2.42; 95% CI 1.67–3.51, p < 0.001; highest tertile HR 5.63; 95% CI 3.99–7.94, p < 0.001) (Fig. 3).

Fig. 3
figure 3

Major adverse cardiovascular events stratified by retinal predicted CVD risk. Kaplan-Meier curves showing MACE-free survival at 10 years stratified by tertiles of retinal predicted CVD risk derived from individual-level prediction at baseline

Incremental prognostic value of retinal PredMACE10 in addition to PCE risk score and CHD PRS

After adjustment for the CHD polygenic risk score, higher retinal predMACE10 remained significantly associated with increased 10-year risk of MACE incidence (HR per 1% increase 1.05; 95% CI 1.04–1.06, p < 0.001). Similarly, higher retinal predMACE10 remained significantly associated with higher risk of 10-year MACE (HR per 1% increase 1.03; 95% CI 1.02–1.04, p < 0.001) after adjustment for the PCE risk score. In a model incorporating both the CHD polygenic risk score and the PCE risk score retinal predMACE10 remained independently associated with MACE at 10 years (HR per 1% increase 1.03; 95% CI 1.02–1.04, p < 0.001).

The retinal predMACE10 score demonstrated similar performance to the PCE risk score for prediction of MACE (AUC—retinal predMACE10 0.697 (95% CI 0.663–0.731); PCE 0.697 (95% CI 0.663–0.730) and better performance than the CHD polygenic risk score alone (AUC 0.575; 95% CI 0.538–0.613). There was a small but significant improvement in risk prediction when the PCE risk and CHD polygenic risk scores were added to the retinal predMACE10 score (0.728; 95% CI 0.695–0.761).

Longitudinal changes in retinal predicted MACE and subsequent clinical outcomes

1,072 individuals in the test cohort had an available second retinal photograph within 3 years of the baseline image. In this group, there were 226 MACE events after the second image. The median change in retinal predMACE10 was + 2.4% (IQR − 0.6–5.4%).

Individuals in the top quintile of predMACE10 difference (i.e. those with the largest increase in predMACE10 between the first and second retinal images, ranging from an increase of 3.2–37% per year) had a significantly higher risk of incident MACE (HR 1.54, 95% CI 1.14–2.08, p = 0.005) compared to the rest of the cohort (Fig. 4).

Fig. 4
figure 4

Major Adverse Cardiovascular Events Based on Progression of Retinal Predicted CVD Risk

Kaplan-Meier curves showing MACE-free survival in individuals who had a second retinal photograph, comparing individuals with most accelerated progression of retinal predicted CVD risk (top 20%) and those with slower progression or reduced retinal predicted CVD risk (bottom 80%).

Discussion

We have identified a number of key findings. First, we have shown that it is feasible to use a deep learning approach to predict CVD risk in individuals with type 2 diabetes without CVD at baseline from routine clinical retinal photographs. Second, we have shown that this retinal CVD risk prediction model (predMACE10) performs similarly to a traditional clinical risk score and is able to predict major adverse cardiovascular events independent of both traditional clinical and genomic CVD risk. Third, we have shown that the addition of retinal predMACE10 to genomic risk provides incremental prognostic value for prediction of MACE compared to a traditional clinical risk score. Finally, we have shown that increases in retinal predMACE10 are also associated with increased risk of future MACE, raising the possibility that this model could be used in routine annual diabetic eye screening.

There has been significant research interest in using deep-learning approaches to predict CVD risk and outcomes from retinal photography [39]. Previous studies have shown that AI can be used to predict CVD risk factors and outcomes such as mortality with high diagnostic accuracy [18,19,20]. This concept is particularly attractive: if CVD risk could be ascertained from routine clinical retinal images, this could avoid the need for multiple healthcare attendances. Our work confirms, in line with other studies, that a deep learning methodology is able to predict CVD risk, but, to our best knowledge, ours is the first study to compare retinal CVD risk prediction and the Pooled Cohort Equation clinical ASCVD risk score, demonstrating similar performance. Since retinal photographs can be relatively easily obtained but require expertise to interpret [40], an AI approach could allow adoption of rapid retinal CVD risk assessment in routine clinical care, including primary care (optometrists), saving time on regular blood pressure and cholesterol monitoring (for example) and reducing congestion in secondary care [41].

It has long been recognised that there is a significant genetic component to CVD risk, however beyond a simple binary assessment of the presence of a family member with established CVD, most clinical risk scores do not take this into account. Over the past decade advances in genotyping have allowed development of whole-genome polygenic risk scores that can predict coronary heart disease and CVD events [24, 42, 43]. As we have previously shown [11], individual markers of CVD risk can be combined with a polygenic risk score to provide incremental prediction of CVD risk. AI systems could make retinal assessment less time-consuming and more reproducible, and in this study the combination of retinal imaging and a polygenic risk score had incremental prognostic value over a standard clinical risk score.

The reproducibility of AI could provide further benefit if used as a serial assessment of CVD risk. We have, for the first time, shown that particularly high-risk patients can be identified by tracking longitudinal progression of retinal predMACE10. This might allow clinical implementation as it provides a measure of CVD risk that could be used as a treatment target. Furthermore, our algorithm was able to predict CVD risk with similar accuracy regardless of whether a single eye photograph or both were used, and whether the left or right eye was used. Again, this could aid clinical applicability as up to 20% patients will not have optimal imaging [44].

Our study does have limitations. It predominantly included individuals aged between 60 and 75; further work is required to investigate performance in older and younger individuals. As we used clinically available electronic health record data not all measurements were taken at the same time as retinal photographs were performed. Furthermore, the second retinal imaging assessment was not performed at the same time in each individual as these were performed as per clinical need—further studies should explore a standardized protocol for retinal image follow-up. Although we used electronic health record data, because of the unique patient identifier used in the Scottish healthcare system we are able to follow all individuals who stay within Scotland for clinical outcomes and vital status throughout life. Image size was reduced to 260 × 260, as is common in deep learning systems due to limitation of computational resources; larger images might improve performance and future work could assess whether would provide additional prognostic value. Our cohort was predominantly Caucasian; so further work is required to ensure that this approach can be used in other ethnicities. Further work is also required to externally validate and compare our ML model to other available algorithms.

In conclusion, we have shown that a deep-learning artificial intelligence methodology is able to accurately predict CVD risk from routine diabetic retinopathy screening photographs with a comparable performance to traditional clinical risk assessment. Longitudinal assessment of retinal risk also predicted adverse CVD events. Implementation of this AI approach alongside retinal screening programmes might be an efficient method of CVD risk assessment. Combining the AI-derived retinal risk with a coronary heart disease polygenic risk score improved risk prediction and underlines the future potential of a one-stop CVD risk assessment derived from retinal screening. Future work should be performed to help translate our exploratory findings into routine clinical practice.

Data availability

Patient electronic healthcare data was collected with approval from NHS Tayside Health Board by the Health Informatics Centre, University of Dundee (HIC). It was de-identified and linked with required fields specifically for this study, however, is still considered sensitive data and cannot be made publicly available. Requests for access can be made to HIC with the appropriate approvals.

Abbreviations

AI:

Artificial intelligence

CLAHE:

Contrast limited adaptive histogram equalization

CNN:

Convolutional neural network

CVD:

Cardiovascular disease

DL:

Deep learning

DRS:

Diabetes retinal screening

EHR:

Electronic health records

GLP1:

RA–glucagon–like peptide 1 receptor agonists

GoDARTS:

Genetics of diabetes audit and research in Tayside Scotland

Grad:

CAM–gradient–based class activation mapping

GWAS:

Genome–wide association study

HIC:

Health informatics centre

ICD:

International classification of diseases

MACE:

Major adverse cardiovascular events

MI:

Myocardial infarction

PCE:

Pooled cohort equations

predMACE10 :

Retinal 10–year MACE prediction

PRS:

Polygenic risk score

SGLT2i:

Sodium–glucose co–transporter 2 inhibitors

References

  1. Low Wang CC, Hess CN, Hiatt WR, Goldfine AB. Clinical update: cardiovascular disease in diabetes mellitus: atherosclerotic cardiovascular disease and heart failure in type 2 diabetes mellitus—mechanisms, management, and clinical considerations. Circulation. 2016;133:2459–502.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  2. Booth GL, Kapral MK, Fung K, Tu JV. Relation between age and cardiovascular disease in men and women with diabetes compared with non-diabetic people: a population-based retrospective cohort study. Lancet. 2006;368:29–36.

    Article  PubMed  Google Scholar 

  3. American Diabetes Association Professional Practice C. 10. Cardiovascular disease and risk management: standards of medical care in diabetes-2022. Diabetes Care. 2022;45:S144-S174.

  4. McGuire DK, Shih WJ, Cosentino F, Charbonnel B, Cherney DZI, Dagogo-Jack S, Pratley R, Greenberg M, Wang S, Huyck S, Gantz I, Terra SG, Masiukiewicz U, Cannon CP. Association of SGLT2 inhibitors with cardiovascular and kidney outcomes in patients with type 2 diabetes: a meta-analysis. JAMA Cardiol. 2021;6:148–58.

    Article  PubMed  Google Scholar 

  5. Sattar N, Lee MMY, Kristensen SL, Branch KRH, Del Prato S, Khurmi NS, Lam CSP, Lopes RD, McMurray JJV, Pratley RE, Rosenstock J, Gerstein HC. Cardiovascular, mortality, and kidney outcomes with GLP-1 receptor agonists in patients with type 2 diabetes: a systematic review and meta-analysis of randomised trials. Lancet Diabetes Endocrinol. 2021;9:653–62.

    Article  CAS  PubMed  Google Scholar 

  6. Amoaku WM, Ghanchi F, Bailey C, Banerjee S, Banerjee S, Downey L, Gale R, Hamilton R, Khunti K, Posner E, Quhill F, Robinson S, Setty R, Sim D, Varma D, Mehta H. Diabetic retinopathy and diabetic macular oedema pathways and management: UK consensus working group. Eye. 2020;34:1–51.

    Article  PubMed  PubMed Central  Google Scholar 

  7. Cheung N, Wang JJ, Klein R, Couper DJ, Sharrett AR, Wong TY. Diabetic retinopathy and the risk of coronary heart disease: the atherosclerosis risk in communities study. Diabetes Care. 2007;30:1742–6.

    Article  PubMed  Google Scholar 

  8. Cheung N, Wang JJ, Rogers SL, Brancati F, Klein R, Sharrett AR, Wong TY, Investigators AS. Diabetic retinopathy and risk of heart failure. J Am Coll Cardiol. 2008;51:1573–8.

    Article  PubMed  Google Scholar 

  9. Wong KH, Hu K, Peterson C, Sheibani N, Tsivgoulis G, Majersik JJ, de Havenon AH. Diabetic retinopathy and risk of stroke: a secondary analysis of the ACCORD eye study. Stroke. 2020;51:3733–6.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  10. Mordi IR, Tee A, Palmer CN, McCrimmon RJ, Doney ASF, Lang CC. Microvascular disease and heart failure with reduced and preserved ejection fraction in type 2 diabetes. ESC Heart Fail. 2020;7:1168–77.

    Article  PubMed  PubMed Central  Google Scholar 

  11. Mordi IR, Trucco E, Syed MG, MacGillivray T, Nar A, Huang Y, George G, Hogg S, Radha V, Prathiba V, Anjana RM, Mohan V, Palmer CNA, Pearson ER, Lang CC, Doney ASF. Prediction of major adverse cardiovascular events from retinal, clinical, and genomic data in individuals with type 2 diabetes: a population cohort study. Diabetes Care. 2022;45:710–6.

    Article  CAS  PubMed  Google Scholar 

  12. Liew G, Mitchell P, Rochtchina E, Wong TY, Hsu W, Lee ML, Wainwright A, Wang JJ. Fractal analysis of retinal microvasculature and coronary heart disease mortality. Eur Heart J. 2011;32:422–9.

    Article  PubMed  Google Scholar 

  13. Sandoval-Garcia E, McLachlan S, Price AH, MacGillivray TJ, Strachan MWJ, Wilson JF, Price JF. Retinal arteriolar tortuosity and fractal dimension are associated with long-term cardiovascular outcomes in people with type 2 diabetes. Diabetologia. 2021;64:2215–27.

    Article  PubMed  PubMed Central  Google Scholar 

  14. McGeechan K, Liew G, Macaskill P, Irwig L, Klein R, Klein BE, Wang JJ, Mitchell P, Vingerling JR, Dejong PT, Witteman JC, Breteler MM, Shaw J, Zimmet P, Wong TY. Meta-analysis: retinal vessel caliber and risk for coronary heart disease. Ann Intern Med. 2009;151:404–13.

    Article  PubMed  PubMed Central  Google Scholar 

  15. Zekavat SM, Raghu VK, Trinder M, Ye Y, Koyama S, Honigberg MC, Yu Z, Pampana A, Urbut S, Haidermota S, O’Regan DP, Zhao H, Ellinor PT, Segre AV, Elze T, Wiggs JL, Martone J, Adelman RA, Zebardast N, Del Priore L, Wang JC, Natarajan P. Deep learning of the retina enables phenome-and genome-wide analyses of the microvasculature. Circulation. 2022;145:134–50.

    Article  PubMed  Google Scholar 

  16. Arnould L, Binquet C, Guenancia C, Alassane S, Kawasaki R, Daien V, Tzourio C, Kawasaki Y, Bourredjem A, Bron A, Creuzot-Garcher C. Association between the retinal vascular network with Singapore I Vessel Assessment (SIVA) software, cardiovascular history and risk factors in the elderly: the Montrachet study, population-based study. PLoS ONE. 2018;13:e0194694.

    Article  PubMed  PubMed Central  Google Scholar 

  17. Wagner SK, Fu DJ, Faes L, Liu X, Huemer J, Khalid H, Ferraz D, Korot E, Kelly C, Balaskas K, Denniston AK, Keane PA. Insights into systemic disease through retinal imaging-based oculomics. Transl Vis Sci Technol. 2020;9:6.

    Article  PubMed  PubMed Central  Google Scholar 

  18. Poplin R, Varadarajan AV, Blumer K, Liu Y, McConnell MV, Corrado GS, Peng L, Webster DR. Prediction of cardiovascular risk factors from retinal fundus photographs via deep learning. Nat Biomed Eng. 2018;2:158–64.

    Article  PubMed  Google Scholar 

  19. Rim TH, Lee CJ, Tham YC, Cheung N, Yu M, Lee G, Kim Y, Ting DSW, Chong CCY, Choi YS, Yoo TK, Ryu IH, Baik SJ, Kim YA, Kim SK, Lee SH, Lee BK, Kang SM, Wong EYM, Kim HC, Kim SS, Park S, Cheng CY, Wong TY. Deep-learning-based cardiovascular risk stratification using coronary artery calcium scores predicted from retinal photographs. Lancet Digit Health. 2021;3:e306–16.

    Article  CAS  PubMed  Google Scholar 

  20. Chang J, Ko A, Park SM, Choi S, Kim K, Kim SM, Yun JM, Kang U, Shin IH, Shin JY, Ko T, Lee J, Oh BL, Park KH. Association of cardiovascular mortality and deep learning-funduscopic atherosclerosis score derived from retinal fundus images. Am J Ophthalmol. 2020;217:121–30.

    Article  PubMed  Google Scholar 

  21. Arnould L, Guenancia C, Bourredjem A, Binquet C, Gabrielle PH, Eid P, Baudin F, Kawasaki R, Cottin Y, Creuzot-Garcher C, Jacquir S. Prediction of cardiovascular parameters with supervised machine learning from Singapore I vessel assessment and OCT-angiography: a pilot study. Transl Vis Sci Technol. 2021;10:20.

    Article  PubMed  PubMed Central  Google Scholar 

  22. Hébert HL, Shepherd B, Milburn K, Veluchamy A, Meng W, Carr F, Donnelly LA, Tavendale R, Leese G, Colhoun HM, Dow E, Morris AD, Doney AS, Lang CC, Pearson ER, Smith BH, Palmer CNA. Cohort profile: genetics of diabetes audit and research in Tayside Scotland (GoDARTS). Int J Epidemiol. 2018;47:380–j381.

    Article  PubMed  Google Scholar 

  23. Arnett DK, Blumenthal RS, Albert MA, Buroker AB, Goldberger ZD, Hahn EJ, Himmelfarb CD, Khera A, Lloyd-Jones D, McEvoy JW, Michos ED, Miedema MD, Munoz D, Smith SC Jr., Virani SS, Williams KA, Sr., Yeboah J, Ziaeian B. 2019 ACC/AHA guideline on the primary prevention of cardiovascular disease: a report of the American college of cardiology/American heart association task force on clinical practice guidelines. Circulation. 2019;140:e596–646.

    PubMed  PubMed Central  Google Scholar 

  24. Khera AV, Chaffin M, Aragam KG, Haas ME, Roselli C, Choi SH, Natarajan P, Lander ES, Lubitz SA, Ellinor PT, Kathiresan S. Genome-wide polygenic scores for common diseases identify individuals with risk equivalent to monogenic mutations. Nat Genet. 2018;50:1219–24.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  25. Syed MG, Doney A, George G, Mordi I, Trucco E. Are cardiovascular risk scores from genome and retinal image complementary? A deep learning investigation in a diabetic cohort. 2021:109–18.

  26. Syed MG. Investigating the retina as a source of biomarkers for systemic conditions using artificial intelligence. 2023.

  27. Tan M, Le Q. EfficientNet: rethinking model scaling for convolutional neural networks. Paper presented at: proceedings of the 36th international conference on machine learning; 2019; Proceedings of machine learning research.

  28. Russakovsky O, Deng J, Su H, Krause J, Satheesh S, Ma S, Huang Z, Karpathy A, Khosla A, Bernstein M, Berg AC, Fei-Fei L. ImageNet large scale visual recognition challenge. Int J Comput Vis. 2015;115:211–52.

    Article  Google Scholar 

  29. Yu Y, Lin H, Meng J, Wei X, Guo H, Zhao Z. Deep transfer learning for modality classification of medical images. Information. 2017;8:91.

    Article  Google Scholar 

  30. HIC. Health Informatics Centre Services. 2021, June 2.

  31. Bradski G, Kaehler A, OpenCV. Dr Dobb’s J Softw Tools. 2000;3.

  32. Pedregosa F, Varoquaux G, Gramfort A, Michel V, Thirion B, Grisel O, Blondel M, Prettenhofer P, Weiss R, Dubourg V. Scikit-learn: machine learning in python. J Mach Learn Res. 2011;12:2825–30.

    Google Scholar 

  33. Chollet F. keras. 2015.

  34. Qubvel EN. 2021, June 21.

  35. Selvaraju RR, Cogswell M, Das A, Vedantam R, Parikh D, Batra D. Grad-cam: Visual explanations from deep networks via gradient-based localization. Proceedings of the IEEE international conference on computer vision. 2017:618–626.

  36. Gerrits N, Elen B, Van Craenendonck T, Triantafyllidou D, Petropoulos IN, Malik RA, De Boever P. Age and sex affect deep learning prediction of cardiometabolic risk factors from retinal images. Sci Rep. 2020;10:1–9.

    Article  Google Scholar 

  37. Rim TH, Lee G, Kim Y, Tham Y-C, Lee CJ, Baik SJ, Kim YA, Yu M, Deshmukh M, Lee BK. Prediction of systemic biomarkers from retinal photographs: development and validation of deep-learning algorithms. Lancet Digit Health. 2020;2:e526–36.

    Article  PubMed  Google Scholar 

  38. Kim YD, Noh KJ, Byun SJ, Lee S, Kim T, Sunwoo L, Lee KJ, Kang S-H, Park KH, Park SJ. Effects of hypertension, diabetes, and smoking on age and sex prediction from retinal fundus images. Sci Rep. 2020;10:1–14.

    Google Scholar 

  39. Wong DYL, Lam MC, Ran A, Cheung CY. Artificial intelligence in retinal imaging for cardiovascular disease prediction: current trends and future directions. Curr Opin Ophthalmol. 2022;33:440–6.

    Article  PubMed  Google Scholar 

  40. Fenner BJ, Wong RLM, Lam WC, Tan GSW, Cheung GCM. Advances in retinal imaging and applications in diabetic retinopathy screening: a review. Ophthalmol Ther. 2018;7:333–46.

    Article  PubMed  PubMed Central  Google Scholar 

  41. Mordi I, Trucco E. The eyes as a window to the heart: looking beyond the horizon. Br J Ophthalmol. 2022;106:1627–8.

    Article  PubMed  Google Scholar 

  42. Riveros-Mckay F, Weale ME, Moore R, Selzam S, Krapohl E, Sivley RM, Tarran WA, Sorensen P, Lachapelle AS, Griffiths JA, Saffari A, Deanfield J, Spencer CCA, Hippisley-Cox J, Hunter DJ, O’Sullivan JW, Ashley EA, Plagnol V, Donnelly P. Integrated polygenic tool substantially enhances coronary artery disease prediction. Circ Genom Precis Med. 2021;14:e003304.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  43. Inouye M, Abraham G, Nelson CP, Wood AM, Sweeting MJ, Dudbridge F, Lai FY, Kaptoge S, Brozynska M, Wang T, Ye S, Webb TR, Rutter MK, Tzoulaki I, Patel RS, Loos RJF, Keavney B, Hemingway H, Thompson J, Watkins H, Deloukas P, Di Angelantonio E, Butterworth AS, Danesh J, Samani NJ, Group UKBCCCW. Genomic risk prediction of coronary artery disease in 480,000 adults: implications for primary prevention. J Am Coll Cardiol. 2018;72:1883–93.

    Article  PubMed  PubMed Central  Google Scholar 

  44. Scanlon PH, Malhotra R, Thomas G, Foy C, Kirkpatrick JN, Lewis-Barned N, Harney B, Aldington SJ. The effectiveness of screening for diabetic retinopathy by digital imaging photography and technician ophthalmoscopy. Diabet Med. 2003;20:467–74.

    Article  CAS  PubMed  Google Scholar 

Download references

Acknowledgements

Not applicable.

Author information

Authors and Affiliations

Authors

Contributions

MS, ET, AD and IM conceived the idea for the study. MS and IM performed the analysis and wrote the initial manuscript. ET, AD, EP, RM, CL, MM and CP provided critical revision.

Corresponding author

Correspondence to Ify R. Mordi.

Ethics declarations

Ethics approval and consent to participate

Data analysis was conducted within the Health Informatics Centre Trusted Research Environment under its overarching ethics approval for approved researchers to conduct research within their secure environment (East of Scotland Research Ethics Committee reference: 18/ES/0126). As the study data are de-identified, consent from individual patients was not required.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Syed, M.G., Trucco, E., Mookiah, M.R.K. et al. Deep-learning prediction of cardiovascular outcomes from routine retinal images in individuals with type 2 diabetes. Cardiovasc Diabetol 24, 3 (2025). https://doi.org/10.1186/s12933-024-02564-w

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12933-024-02564-w

Keywords