# Prediction of Cardiovascular Parameters With Supervised Machine Learning From Singapore “I” Vessel Assessment and OCT-Angiography: A Pilot Study

Louis Arnould, Charles Guenancia, Abderrahmane Bourredjem, Christine Binquet, Pierre-Henry Gabrielle, Pétra Eid, Florian Baudin, Ryo Kawasaki, Yves Cottin, Catherine Creuzot-Garcher, et al.

<figure><img src="image_1.png" /></figure>

## To cite this version:

Louis Arnould, Charles Guenancia, Abderrahmane Bourredjem, Christine Binquet, Pierre-Henry Gabrielle, et al. Prediction of Cardiovascular Parameters With Supervised Machine Learning From Singapore “I” Vessel Assessment and OCT-Angiography: A Pilot Study. Translational vision science & technology, 2021, 10 (13), pp.20. 10.1167/tvst.10.13.20. hal-03435632

HAL Id: hal-03435632

https://hal.science/hal-03435632

Submitted on 3 Dec 2021

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_2.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
# Prediction of Cardiovascular Parameters With Supervised Machine Learning From Singapore “I” Vessel Assessment and OCT-Angiography: A Pilot Study

Louis Arnould<fnref n="1" />⁻<fnref n="3" />, Charles Guenancia<fnref n="4" />,<fnref n="5" />, Abderrahmane Bourredjem<fnref n="2" />, Christine Binquet<fnref n="2" />, Pierre-Henry Gabrielle<fnref n="1" />,<fnref n="3" />, Pétra Eid<fnref n="1" />, Florian Baudin<fnref n="1" />, Ryo Kawasaki<fnref n="6" />, Yves Cottin<fnref n="4" />,<fnref n="5" />, Catherine Creuzot-Garcher<fnref n="1" />,<fnref n="3" />, and Sabir Jacquir<fnref n="7" />

**Correspondence:** Louis Arnould, Ophthalmology Department, University Hospital, 14 rue Paul Gaffarel, 21079 Dijon CEDEX, France. e-mail: louis.arnould@chu-dijon.fr

**Received:** March 1, 2021

**Accepted:** August 11, 2021

**Published:** November 12, 2021

**Keywords:** retina; supervised machine learning (ML), cardiovascular risk score, optical coherence tomography angiography (OCT-A)

**Citation:** 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(13):20, https://doi.org/10.1167/tvst.10.13.20

**Purpose:** Assessment of cardiovascular risk is the keystone of prevention in cardiovascular disease. The objective of this pilot study was to estimate the cardiovascular risk score (American Hospital Association [AHA] risk score, Syntax risk, and SCORE risk score) with machine learning (ML) model based on retinal vascular quantitative parameters.

**Methods:** We proposed supervised ML algorithm to predict cardiovascular parameters in patients with cardiovascular diseases treated in Dijon University Hospital using quantitative retinal vascular characteristics measured with fundus photography and optical coherence tomography – angiography (OCT-A) scans (alone and combined). To describe retinal microvascular network, we used the Singapore “I” Vessel Assessment (SIVA), which extracts vessel parameters from fundus photography and quantitative OCT-A retinal metrics of superficial retinal capillary plexus.

**Results:** The retinal and cardiovascular data of 144 patients were included. This paper presented a high prediction rate of the cardiovascular risk score. By means of the Naïve Bayes algorithm and SIVA + OCT-A data, the AHA risk score was predicted with 81.25% accuracy, the SCORE risk with 75.64% accuracy, and the Syntax score with 96.53% of accuracy.

**Conclusions:** Performance of these algorithms demonstrated in this preliminary study that ML algorithms applied to quantitative retinal vascular parameters with SIVA software and OCT-A were able to predict cardiovascular scores with a robust rate. Quantitative retinal vascular biomarkers with the ML strategy might provide valuable data to implement predictive model for cardiovascular parameters.

**Translational Relevance:** Small data set of quantitative retinal vascular parameters with fundus and with OCT-A can be used with ML learning to predict cardiovascular parameters.

<footnote n="1">¹ Ophthalmology Department, University Hospital, Dijon, France</footnote>

<footnote n="2">² INSERM, CIC1432, Clinical Epidemiology Unit, Dijon, France; Dijon University Hospital, Clinical Investigation Center, Clinical Epidemiology/Clinical Trials Unit, Dijon, France</footnote>

<footnote n="3">³ Centre des Sciences du Gout et de l'Alimentation, AgroSup Dijon, CNRS, INRAE, Université Bourgogne Franche-Comté, Dijon, France</footnote>

<footnote n="4">⁴ Cardiology Department, University Hospital, Dijon, France</footnote>

<footnote n="5">⁵ PEC 2, University Hospital, Dijon, France</footnote>

<footnote n="6">⁶ Department of Vision Informatics, Osaka University Graduate School of Medicine, Suita, Japan</footnote>

<footnote n="7">⁷ Université Paris-Saclay, CNRS, Institut des Neurosciences Paris-Saclay, Gif-sur-Yvette, France</footnote>
Introduction

Retinal vascular imaging is constantly improving. The retinal microvascular network can be thoroughly described with fundus photograph imaging analysis software, such as the Singapore “I” Vessel Assessment (SIVA).¹,² In addition, the quantitative description of retinal microvascularization was recently enhanced by the development of a new noninvasive technique: optical coherence tomography angiography (OCT-A).³ It has been shown that quantitative retinal vascular characteristics in both fundus photographs with SIVA software and metrics of retinal vascular density with OCT-A are associated with patients’ systemic vascular alteration,⁴ cardiovascular risk profile⁵,⁶ and cardio-vascular complications.⁷ In addition, deep learning algorithms have been described for diabetic retinopa-thy detection, retinopathy of prematurity screening, grading of age-related macular degeneration, and glaucoma detection.⁸⁻¹¹ Furthermore, Poplin et al. recently indicated that retinal imaging analysis of fundus photography using deep learning could be used to predict a wide range of cardiovascular risk factors.¹² However, deep learning is limited in some healthcare applications, particularly in a context of sparse data and real world clinical data.¹³ Thus, machine learn-ing (ML) methods could be trained more easily and provide better overall performance when compared to deep learning with a small data set. We hypothesized that classic ML algorithms could predict cardiovascu-lar risk scores and systemic parameters from quantita-tive retinal vascular data obtained with SIVA software and OCT-A. Better cardiovascular risk stratification is consequently of growing interest given that cardio-vascular disease (CVD) remains one of the leading causes of death worldwide.¹⁴ The purpose of this pilot study was to develop a prediction model with a supervised ML approach to cardiovascular parameters using retinal vascular characteristics measured with SIVA software and OCT-A (alone and combined). Our objective was to estimate the cardiovascular risk score (American Hospital Association [AHA] risk score, Syntax score, and SCORE risk) with this ML model.

Patients and Methods

Study Design and Patients

This pilot study was an ancillary study of a previ-
ous pilot prospective cross-sectional study conducted
in Dijon University Hospital's Cardiology Inten-
sive Care Unit. The methodology of the EYE-MI

study and the patients’ baseline characteristics have
been detailed elsewhere.⁵ Briefly, from May 2016 to
May 2017, patients presenting with acute coronary
syndrome (ACS) were included. They were taken to
the ophthalmology department within the first 2 days
of their hospitalization for an examination of the
retinal microvasculature using OCT-A and fundus
photographs. The exclusion criteria were: retinal
disease (vascular occlusion, diabetic retinopathy, and
macular degeneration), patients under 18 years of age,
those under guardianship, patients with hemodynamic
instability, and patients without both retinal examina-
tions (SIVA and OCT-A). We also excluded severely
myopic eyes (axial length greater than 26 mm) because
this could affect retinal microvascular density.¹⁵ The
study was approved by the Dijon University Hospital
ethics committee and was registered as 2017-A02095-
48. It complied with the tenets of the Declaration of
Helsinki and a written informed consent was obtained
from the patients. We followed the Strengthening the
Reporting of Observational Studies in Epidemiology
(STROBE) statement according to the EQUATOR
Guidelines.¹⁶

Retinal Microvascular Image Acquisition and
Quantitative Analysis

After inclusion, the patients underwent an OCT-
A examination (CIRRUS HD-OCT, Model 5000;
Carl Zeiss Meditec AG) and 45 degree color retinal
photographs, centered on the optic disc, were obtained
with a fundus camera (TRC NW6S, Topcon, Tokyo,
Japan) for both eyes. This eye examination was
performed under mydriasis obtained with eye drops
containing tropicamide 0.5% (Thea, Clermont-
Ferrand, France). Axial length was measured using an
optical biometer (IOL Master; Carl Zeiss Meditec AG,
Jena, Germany). Quantitative OCT-A metrics were
obtained with the angiography software (Angioplex,
version 10; Carl Zeiss Meditec AG; Fig. 1A). We inves-
tigated the retinal vascular features in the superficial
capillary plexus (SCP) and the following measure-
ments were taken in the 3 × 3 mm angiograms: the
area of the foveal avascular zone (FAZ; mm²), perfu-
sion density (area, unitless), and vessel density (length,
mm⁻¹). These densities were both measured in the
inner and full Early Treatment Diabetic Retinopathy
Study (EDTRS) circle.

In addition, fundus photographs were anonymously
sent to the reading center in Yamagata University,
Japan (author R.K.), and a single trained grader
extracted retinal vessel characteristics with the SIVA
software (Fig. 1B). The computerized analysis of the
<figure><img src="image_3.png" /><figcaption>**Figure 1.** Retinal microvascular image. (A) Retinal superficial capillary plexus with optical coherence tomography angiography examinations. (B) Forty-five-degree color retinal photographs, centered on the optic disc.</figcaption></figure>

retinal vascular network computerized analysis was based on the analysis of vessels from the center of the optic disc and then to three successive zones corresponding to 0.5 (zone A), 1 (zone B), and 2 (zone C) disc diameter. The six largest arterioles and veins were analyzed.

Thus, only one eye (with both OCT-A and SIVA examinations) was retained for analysis for each participant and its selection followed the criteria described below: (1) fundus photograph and OCT-A of the right eye for participants born in even-numbered years and the left eye for those born in odd-numbered years; (2) in single-eye patients, the functional eye was selected; and (3) when SIVA or OCT-A were uninterpretable on one eye, the other one was retained for analysis. Only OCT-A images with a signal strength >7/10 were retained. For this ancillary study we kept 144 eyes with both OCT-A and SIVA examinations.

## Cardiovascular Data Collection

As described previously, we extracted cardiovascular data from medical records and observation sheets used by the Observatoire des Infarctus de Côte d'Or (RICO).⁵ The following data were collected: age, sex, high blood pressure, diabetes mellitus history, body mass index (BMI), hypercholesterolemia, and current smoking. From the above data, cardiovascular risk scores defined by the AHA (AHA risk score) for a high-risk population were calculated. The AHA risk score includes age, sex, the ethnic origin, the history of arterial hypertension and diabetes, active smoking, systolic, and diastolic arterial pressure and levels of total cholesterol and high-density lipoprotein (HDL)

cholesterol levels.¹⁷ The anatomic Syntax score, a risk stratification score for coronary lesions (length, bifurcation, diffuse disease, calcifications, thrombus, and total occlusion) was determined for all of the patients who underwent coronarography.¹⁸ We finally calculated the SCORE risk.¹⁹

## Pre-Processing of OCT-A and SIVA Data

The dispersion of the value of the different parameter is significant from one patient to another and from one parameter to another. The scale of values also differs. To provide a common referential for all data, these were normalized between 0 and 1. The SIVA and OCT-A settings used in this study are listed in Figures 1 and 2. Figure 2A and Figure 3A illustrate this variability and dispersion of values. Figure 2B and Figure 3B show the normalized data.

## Machine Learning Approach

The method used can be summarized by the following conceptual framework:

1. To build a prediction model of the cardiovascular risk score and parameters based on supervised machine learning approach. The process applied a set of known input (OCT-A and/or SIVA data) with response data (cardiovascular risk score and parameters) and created a model to produce reasonable predictions of these cardiovascular parameters. In this study, supervised learning was based on classification (K-nearest
<figure><img src="image_4.png" /><figcaption>**Figure 2.** Distribution of optical coherence tomography angiography parameters values 1: Foveal avascular zone; 2: Vessel inner; 3: Vessel full; 4: Perfusion inner; 5: Perfusion full. (A) Original data, (B) normalized data between 0 and 1.</figcaption></figure>

<figure><img src="image_5.png" /><figcaption>**Figure 3.** Distribution of Singapore "I" Vessel Assessment parameters values: 1: Biggest six arterioles in Zone B 2: Biggest six veins in Zone B 3: Arteriole - Venular Ratio of Zone B 4: Biggest six arterioles in Zone C 5: Biggest six veins in Zone C 6: Arteriole - Venular Ratio of Zone C 7: Fractal dimension total zone C 8: Fractal dimension arterioles zone C 9: Fractal dimension venules zone C 10: Simple tortuosity arteriole 11: Curvature tortuosity arteriole 12: Simple tortuosity venule 13: Curvature tortuosity venule 14: Simple tortuosity vessels 15: Curvature tortuosity vessels. (A) Original data, (B) normalized data between 0 and 1.</figcaption></figure>

neighbors, discriminant analysis, and Naïve Bayes) and regression (decision trees) techniques were used in order to develop these predictive models.

2. To use the predictive model obtained in the previous step to estimate the cardiovascular parameters of patients depending on OCT-A and SIVA retinal vascular features.

Further details on the methodology were provided in the Supplementary Material. The cardiovascular parameters were divided into two groups. The first group consisted of the cardiovascular risk score: AHA Risk score, Syntax score, and SCORE risk. Regarding the AHA Risk score in this study, it has not been used as a primary prevention risk assessment as it was calculated in a high-risk profile population. We chose

this score to attest the systemic vascular profile of our population because no therapeutic decision (aspirin, statins, and antihypertensive drugs) was made with regard to this calculation. This group accounted for the primary prediction goal. The second one included age, sex, high blood pressure history, diabetes mellitus history, hypercholesterolemia, current smoking status, and BMI.

In this study, each cardiovascular parameter was categorized using two or three labels (Score risk ≥ 5%, AHA risk ≥ 20%, Score syntax ≥ 33, age > 60 years, and BMI > 25 and > 30). To avoid bias in the choice of the learning group patients, patients were selected randomly. In our algorithm, the minimum number of patients required for the learning group was not defined. The learning group increased in increments of 10 patients. For each increment, the predictive model
was applied to the database in order to investigate the correlation between the prediction rate and the size of the learning group. All programs were written using Matlab software. These programs can be shared on request to the authors. For each data set, the prediction rates of the four ML algorithms were compared. After that, the prediction rates of the same algorithm were compared according to the three data sets (OCT versus SIVA versus OCT + SIVA). Paired tests with post hoc correction using Tukey pair wise multiple comparisons were used.²⁰ This analysis was done using the mixed procedure of SAS statistical software (version 9.4; SAS Institute, Inc., Cary, NC, USA).

## Results

In this study, the data of 144 patients were taken into account. Baseline characteristics of the study participants are presented in Table 1. The mean age

was 61.9 (±12.6) years old and 20.1% were female patients. There was no difference between participants and nonparticipants from the EYE-MI study.⁵ The prediction rates of cardiovascular parameters based on the K-nearest neighbor (KNN) approach and the Naïve Bayes approach are presented in Supplementary Figures S1 and S2. These figures show the value of the prediction rate depending on the number of learning patients.

Figures 4, 5, and 6 displayed the comparison of the prediction results from the four ML algorithms for the cardiovascular risk score.

As shown in Figures 4, 5, and 6, the prediction rate was dependent on the number of learning patients used to build the prediction model. The prediction rates of the four ML techniques are reported in Table 2 from OCT-A data, Table 3 from SIVA data, and Table 4 from OCT-A + SIVA data. There were significant differences between prediction rates obtained from the discriminant analysis, the KNN method, the Naïve Bayes approach, and the decision tree classification

<table><caption>**Table 1.** Baseline Characteristics Between Participants and Nonparticipants</caption><thead><tr><th></th><th>Global EYE-MI Population<br>N = 237</th><th>Participants<br>N = 144</th><th>Nonparticipants<br>N = 93</th><th>P Value</th></tr></thead><tbody><tr><td><strong>Cardiovascular characteristics</strong></td><td></td><td></td><td></td><td></td></tr><tr><td>Age, y</td><td>62.0 (±13.0)</td><td>61.9 (±12.6)</td><td>61.4 (±12.8)</td><td>0.643</td></tr><tr><td>Gender, female</td><td>51.0 (21.5)</td><td>29.0 (20.1)</td><td>22.0 (23.7)</td><td>0.630</td></tr><tr><td>Previous high blood pressure</td><td>121.0 (51.0)</td><td>74.0 (51.4)</td><td>47.0 (50.5)</td><td>1.0</td></tr><tr><td>Previous diabetes</td><td>53.0 (22.4)</td><td>35.0 (24.3)</td><td>18.0 (19.4)</td><td>0.463</td></tr><tr><td>Active smoking</td><td>67.0 (28.3)</td><td>40.0 (27.8)</td><td>27.0 (29.0)</td><td>0.578</td></tr><tr><td>Body mass index, m²/kg</td><td>26.7 (±5.7)</td><td>27.1 (±4.2)</td><td>27.3 (±4.7)</td><td>0.677</td></tr><tr><td>Hypercholesterolemia</td><td>96.0 (40.5)</td><td>61.0 (42.4)</td><td>35.0 (37.6)</td><td>0.556</td></tr><tr><td>Family history of CHD</td><td>80.0 (33.8)</td><td>52.0 (36.1)</td><td>28.0 (30.1)</td><td>0.416</td></tr><tr><td>Systolic pressure at admission, mm Hg</td><td>144.1 (±29.8)</td><td>143.0 (±29.9)</td><td>144.2 (±29.5)</td><td>0.864</td></tr><tr><td>Diastolic pressure at admission, mm Hg</td><td>84.2 (±18.7)</td><td>83.6 (±19.8)</td><td>84.4 (±17.2)</td><td>0.753</td></tr><tr><td>LVEF at admission, %</td><td>54.0 (±10.7)</td><td>53.9 (±11.1)</td><td>54.2 (±10.2)</td><td>0.835</td></tr><tr><td><strong>Vascular history</strong></td><td></td><td></td><td></td><td></td></tr><tr><td>Ischemic coronary heart disease</td><td>51.0 (21.5)</td><td>35.0 (24.3)</td><td>16.0 (17.2)</td><td>0.256</td></tr><tr><td>Carotid atheroma</td><td>10.0 (4.2)</td><td>7.0 (4.9)</td><td>3.0 (3.2)</td><td>0.779</td></tr><tr><td>Peripheral artery disease</td><td>12.0 (5.1)</td><td>8.0 (5.6)</td><td>4.0 (4.3)</td><td>0.899</td></tr><tr><td>Chronic kidney failure</td><td>7.0 (3.0)</td><td>6.0 (4.2)</td><td>1.0 (1.1)</td><td>0.327</td></tr><tr><td>Ischemic stroke</td><td>9.0 (3.8)</td><td>4.0 (2.8)</td><td>5.0 (5.4)</td><td>0.507</td></tr><tr><td><strong>Acute coronary syndrome</strong></td><td></td><td></td><td></td><td></td></tr><tr><td>STEMI</td><td>94.0 (39.7)</td><td>54.0 (37.5)</td><td>40.0 (43.0)</td><td>0.626</td></tr><tr><td>NSTEMI</td><td>113.0 (47.6)</td><td>70.0 (48.6)</td><td>43.0 (46.2)</td><td>0.645</td></tr><tr><td>Unstable angina</td><td>30.0 (12.7)</td><td>20.0 (13.9)</td><td>10.0 (10.8)</td><td>0.756</td></tr><tr><td><strong>Cardiovascular risk scores</strong></td><td></td><td></td><td></td><td></td></tr><tr><td>AHA Risk score</td><td>19.8 (±14.5)</td><td>18.7 (±14.7)</td><td>21.5 (±14.1)</td><td>0.145</td></tr><tr><td>Syntax score</td><td>11.6 (±9.5)</td><td>11.5 (±9.7)</td><td>11.9 (±9.2)</td><td>0.759</td></tr><tr><td>SCORE risk</td><td>3.5 (±2.8)</td><td>3.2 (2.3)</td><td>2.7 (±1.7)</td><td>0.364</td></tr></tbody></table>

The results are displayed as *n* (%) for categorical variables and as mean and standard deviation M (±SD) for continuous variables.

CHD, Cardiovascular and heart disease; LVEF, Left ventricular ejection fraction; NSTEMI, Non ST-Elevation Myocardial Infarction; STEMI, ST-Elevation Myocardial Infarction. *P* value for comparison between participants and non-participants.
<figure><img src="image_6.png" /><figcaption>**Figure 4.** Prediction model for cardiovascular risk score with optical coherence tomography angiography data according to four machine learning algorithms (K-nearest neighbors [KNN], discriminant analysis, Naïve Bayes, and decision tree).</figcaption></figure>

<figure><img src="image_7.png" /><figcaption>**Figure 5.** Prediction model for cardiovascular risk score with Singapore “I” Vessel Assessment data according to four machine learning algorithms (K-nearest neighbors [KNN], discriminant analysis, Naïve Bayes, and decision tree).</figcaption></figure>

<figure><img src="image_8.png" /><figcaption>**Figure 6.** Prediction model for cardiovascular risk score with combined optical coherence tomography angiography + Singapore “I” Vessel Assessment data according to four machine learning algorithms (K-nearest neighbors [KNN], discriminant analysis, Naïve Bayes, and decision tree).</figcaption></figure>

with the Type 3 Tests of Fixed Effects. When using multiple test adjustment, two algorithms (Naïve Bayes and KNN) were associated with higher prediction rates compared to the discriminant analysis and the decision tree classification (Supplementary Table).

The ranges of accuracy were for KNN and OCT-A + SIVA data (0.25–0.97), OCT-A data (0.31–0.98), and SIVA data (0.24–0.98). For the Naïve Bayes approach the range of accuracy were OCT-A + SIVA data (0.24–0.98), OCT-A data (0.29–0.98), and SIVA data
<table><caption>**Table 2.** Prediction Rate of the Four Machine Learning Techniques Using the Optical Coherence Tomography Angiography Data (n = 144)</caption><thead><tr><th></th><th>Discriminant Analysis</th><th>KNN</th><th>Naïve Bayes</th><th>Decision Tree</th><th>P Value</th></tr></thead><tbody><tr><td>Syntax score</td><td>95.19 ± 0.74</td><td>96.13 ± 1.08</td><td>96.23 ± 1.88</td><td>95.19 ± 0.19</td><td>0.01</td></tr><tr><td>AHA risk score</td><td>74.40 ± 1.33</td><td>76.09 ± 3.08</td><td>80.31 ± 4.14</td><td>74.65 ± 2.34</td><td>&lt;0.01</td></tr><tr><td>SCORE risk</td><td>73.02 ± 1.36</td><td>76.16 ± 6.33</td><td>76.19 ± 5.30</td><td>73.36 ± 3.24</td><td>0.02</td></tr><tr><td>Age</td><td>63.05 ± 1.62</td><td>70.34 ± 8.55</td><td>73.07 ± 8.26</td><td>68.60 ± 6.42</td><td>&lt;0.01</td></tr><tr><td>Sex, male</td><td>77.63 ± 5.47</td><td>77.43 ± 7.49</td><td>77.63 ± 5.47</td><td>77.28 ± 8.27</td><td>&lt;0.01</td></tr><tr><td>High blood pressure history</td><td>53.62 ± 3.68</td><td>67.56 ± 8.30</td><td>69.79 ± 10.19</td><td>62.45 ± 7.85</td><td>&lt;0.01</td></tr><tr><td>Diabetes mellitus history</td><td>75.45 ± 2.30</td><td>78.08 ± 5.83</td><td>75.45 ± 2.30</td><td>73.91 ± 5.17</td><td>&lt;0.01</td></tr><tr><td>Hypercholesterolemia</td><td>56.15 ± 3.82</td><td>70.44 ± 11.16</td><td>66.47 ± 9.06</td><td>59.72 ± 6.41</td><td>&lt;0.01</td></tr><tr><td>Current smoking</td><td>46.43 ± 4.55</td><td>56.00 ± 12.70</td><td>59.52 ± 13.06</td><td>48.86 ± 10.08</td><td>&lt;0.01</td></tr><tr><td>Body mass index</td><td>40.53 ± 9.32</td><td>58.23 ± 11.95</td><td>58.98 ± 12.97</td><td>51.14 ± 10.86</td><td>&lt;0.01</td></tr></tbody></table>

The prediction rates (%) are displayed as mean and standard deviation (M ± SD).

<table><caption>**Table 3.** Prediction Rate of the Four Machine Learning Techniques Using the Singapore "I" Vessel Assessment Data (n = 144)</caption><thead><tr><th></th><th>Discriminant Analysis</th><th>KNN</th><th>Naïve Bayes</th><th>Decision Tree</th><th>P Value</th></tr></thead><tbody><tr><td>Syntax score</td><td>94.94 ± 0.74</td><td>95.83 ± 1.19</td><td>96.28 ± 1.21</td><td>95.34 ± 0.42</td><td>&lt;0.01</td></tr><tr><td>AHA risk score</td><td>72.82 ± 4.68</td><td>79.22 ± 5.55</td><td>80.61 ± 4.24</td><td>74.50 ± 6.62</td><td>&lt;0.01</td></tr><tr><td>SCORE risk</td><td>65.33 ± 3.25</td><td>70.54 ± 8.56</td><td>74.36 ± 6.17</td><td>66.42 ± 8.87</td><td>&lt;0.01</td></tr><tr><td>Age</td><td>54.56 ± 4.65</td><td>66.27 ± 10.06</td><td>54.56 ± 4.65</td><td>63.69 ± 7.32</td><td>&lt;0.01</td></tr><tr><td>Sex, male</td><td>79.66 ± 0.74</td><td>79.22 ± 12.86</td><td>84.47 ± 6.14</td><td>78.77 ± 7.89</td><td>0.09</td></tr><tr><td>High blood pressure history</td><td>59.03 ± 4.73</td><td>64.83 ± 10.74</td><td>67.71 ± 8.51</td><td>66.52 ± 7.08</td><td>&lt;0.01</td></tr><tr><td>Diabetes mellitus history</td><td>74.21 ± 3.62</td><td>80.36 ± 5.31</td><td>82.29 ± 6.60</td><td>72.87 ± 6.27</td><td>&lt;0.01</td></tr><tr><td>Hypercholesterolemia</td><td>60.22 ± 3.90</td><td>69.84 ± 9.27</td><td>69.94 ± 9.55</td><td>63.89 ± 8.85</td><td>&lt;0.01</td></tr><tr><td>Current smoking</td><td>44.94 ± 5.60</td><td>54.02 ± 15.92</td><td>59.62 ± 16.82</td><td>50.79 ± 11.88</td><td>&lt;0.01</td></tr><tr><td>Body mass index</td><td>47.47 ± 5.25</td><td>59.87 ± 10.84</td><td>60.86 ± 14.48</td><td>53.08 ± 9.48</td><td>&lt;0.01</td></tr></tbody></table>

The prediction rates (%) are displayed as mean and standard deviation (M ± SD).

(0.39–0.97). For both strategies, sensitivity and specificity ranged from 0 to 1. The KNN and the Naïve Bayes approaches more accurately predicted the three cardiovascular risk scores compared to

the discriminant analysis and the decision tree approaches. Finally, we compared the prediction rate with SIVA, OCT-A, and OCT-A + SIVA data for these two algorithms (Naïve Bayes

<table><caption>**Table 4.** Prediction Rate of the Four Machine Learning Techniques Using the Both Optical Coherence Tomography Angiography and the Singapore "I" Vessel Assessment Data (n = 144)</caption><thead><tr><th></th><th>Discriminant Analysis</th><th>KNN</th><th>Naïve Bayes</th><th>Decision Tree</th><th>P Value</th></tr></thead><tbody><tr><td>Syntax score</td><td>95.14 ± 0.82</td><td>95.83 ± 1.25</td><td>96.53 ± 1.25</td><td>95.19 ± 0.51</td><td>&lt;0.01</td></tr><tr><td>AHA risk score</td><td>73.61 ± 3.13</td><td>74.95 ± 3.76</td><td>81.25 ± 3.84</td><td>73.46 ± 6.68</td><td>&lt;0.01</td></tr><tr><td>SCORE risk</td><td>72.22 ± 1.77</td><td>75.25 ± 5.61</td><td>75.64 ± 5.96</td><td>73.12 ± 4.20</td><td>0.03</td></tr><tr><td>Age</td><td>61.36 ± 6.49</td><td>70.63 ± 6.95</td><td>70.14 ± 10.55</td><td>69.20 ± 7.46</td><td>&lt;0.01</td></tr><tr><td>Sex, male</td><td>77.83 ± 6.87</td><td>77.73 ± 11.59</td><td>86.46 ± 5.08</td><td>78.37 ± 8.83</td><td>&lt;0.01</td></tr><tr><td>High blood pressure history</td><td>60.37 ± 4.91</td><td>67.46 ± 10.94</td><td>66.96 ± 9.63</td><td>67.81 ± 7.97</td><td>&lt;0.01</td></tr><tr><td>Diabetes mellitus history</td><td>76.24 ± 3.36</td><td>79.86 ± 4.58</td><td>83.48 ± 5.50</td><td>75.10 ± 5.71</td><td>&lt;0.01</td></tr><tr><td>Hypercholesterolemia</td><td>61.61 ± 5.09</td><td>71.03 ± 12.16</td><td>70.09 ± 9.82</td><td>66.12 ± 8.68</td><td>&lt;0.01</td></tr><tr><td>Current smoking</td><td>48.16 ± 6.58</td><td>58.28 ± 13.67</td><td>59.97 ± 15.54</td><td>50.94 ± 12.27</td><td>&lt;0.01</td></tr><tr><td>Body mass index</td><td>47.82 ± 4.65</td><td>59.57 ± 12.28</td><td>60.37 ± 14.67</td><td>51.89 ± 10.12</td><td>&lt;0.01</td></tr></tbody></table>

The prediction rates (%) are displayed as mean and standard deviation (M ± SD).
<table><caption>**Table 5.** Comparison of Prediction of Cardiovascular Parameters With Machine Learning With OCT-A, SIVA, and OCT-A + SIVA</caption><thead><tr><th></th><th>Reference</th><th>Compared Strategy</th><th>Estimate</th><th>Standard Error</th><th>Adjusted P Value</th></tr></thead><tbody><tr><td><strong>Naïve Bayes</strong></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Age</td><td></td><td></td><td></td><td></td><td></td></tr><tr><td></td><td>OCT-A</td><td>OCT-A + SIVA</td><td>2.93</td><td>1.48</td><td>0.14</td></tr><tr><td></td><td>SIVA</td><td>OCT-A + SIVA</td><td>-2.08</td><td>1.48</td><td>0.35</td></tr><tr><td>High blood pressure history</td><td></td><td></td><td></td><td></td><td></td></tr><tr><td></td><td>OCT-A</td><td>OCT-A + SIVA</td><td>2.83</td><td>0.69</td><td>&lt;0.01</td></tr><tr><td></td><td>SIVA</td><td>OCT-A + SIVA</td><td>0.74</td><td>0.69</td><td>0.53</td></tr><tr><td>Hypercholesterolemia</td><td></td><td></td><td></td><td></td><td></td></tr><tr><td></td><td>OCT-A</td><td>OCT-A + SIVA</td><td>-3.62</td><td>1.21</td><td>0.02</td></tr><tr><td></td><td>SIVA</td><td>OCT-A + SIVA</td><td>-0.15</td><td>1.21</td><td>0.99</td></tr><tr><td><strong>KNN</strong></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>AHA risk score</td><td></td><td></td><td></td><td></td><td></td></tr><tr><td></td><td>OCT-A</td><td>OCT-A + SIVA</td><td>1.14</td><td>1.19</td><td>0.61</td></tr><tr><td></td><td>SIVA</td><td>OCT-A + SIVA</td><td>4.27</td><td>1.19</td><td>&lt;0.01</td></tr><tr><td>SCORE risk</td><td></td><td></td><td></td><td></td><td></td></tr><tr><td></td><td>OCT-A</td><td>OCT-A + SIVA</td><td>0.94</td><td>1.40</td><td>0.78</td></tr><tr><td></td><td>SIVA</td><td>OCT-A + SIVA</td><td>-4.71</td><td>1.40</td><td>&lt;0.01</td></tr><tr><td>Age</td><td></td><td></td><td></td><td></td><td></td></tr><tr><td></td><td>OCT-A</td><td>OCT-A + SIVA</td><td>-0.2976</td><td>1.5351</td><td>0.98</td></tr><tr><td></td><td>SIVA</td><td>OCT-A + SIVA</td><td>-4.3651</td><td>1.5351</td><td>0.02</td></tr></tbody></table>

OCT-A, optical coherence tomography angiography; SIVA, the Singapore "I" Vessel Assessment.

and KNN) and each cardiovascular parameters (Table 5).

## Discussion

In this cross-sectional pilot study, we applied a prediction model of cardiovascular parameters from ML algorithms using retinal vascular characteristics measured with SIVA software and OCT-A (alone and combined). Overall, we observed that this approach may be effective to procure a moderate to robust prediction rate for a specific cardiovascular data set.

To the best of our knowledge, this is the first study to explore the potential interest of using quantitative parameters of the retinal microvascular network measured by means of SIVA software and OCT-A to predict the cardiovascular risk factors burden of high cardiovascular risk profile patients with a non-deep learning model. In this study, we focused on patients with a history of cardiovascular disease (inclusion criteria: ACS). OCT-A and retinal vascular analysis based on fundus photographs with software, such as SIVA, provide us with a large amount of quantitative data. Supervised classifiers help us to

properly take advantage of all these quantitative data. Although deep learning based on convolutional neural networks in ophthalmology is commonly used to detect glaucoma, aged-related macular degeneration and diabetic retinopathy,²¹⁻²⁴ the sparsity of OCT-A and SIVA data added to our relatively small sample size make it difficult to apply this type of deep learning algorithm in this study.²⁵,²⁶ In this regard, our approach of cardiovascular prediction with retinal vascular biomarker is novel and helpful.

Publications concerning the prediction of the cardiovascular risk profile with ML algorithm are steadily increasing in number.²⁷,²⁸ Alaa et al. recently investigated the cardiovascular disease risk prediction with ML methods.²⁹ They conducted a large prospective cohort study and analyzed data on 423,604 participants without cardiovascular disease at baseline in UK Biobank. They used more than 470 variables, notably on health and medical history, dietary and nutritional information, and sociodemographics. They found that their ML algorithm significantly improved the accuracy of cardiovascular disease risk prediction compared to gold standard scoring systems based on conventional risk factors (Framingham score).³⁰

Comparison with Poplin et al.'s results (prediction for age only area under the concentration curve [AUC])
= 0.66 [0.61 to 0.71] and for the risk SCORE AUC = 0.72 [0.67 to 0.76]) is impossible because our strategy of using the quantitative dataset and the supervised ML model is significantly different.¹² The UK Biobank provides a tremendous amount of information in a very large population, making it possible to predict cardiovascular parameters and more recently hemoglobin concentration.³¹ Moreover, Gerrits et al. recently published a deep learning algorithm to predict cardiometabolic risk factors based on 12,000 retinal images from 3000 participants.³² Recently, Cheung et al. suggested creating a “retinal vessel CVD risk score” based on the artificial intelligence system.³³ This deep learning project will have to overcome differences in terms of imaging device, image quality, and protocol standardization. Our approach is more pragmatic and could be reserved for smaller centers with a limited data set. Furthermore, in our algorithms, we were able to include different kinds of retinal image parameters: SIVA and OCT-A combined and not only retinal vessel caliber.

In this study, prediction rates were moderate to high for all the cardiovascular risk factors, which was particularly true for diabetic status, the Syntax risk score, AHA risk score, and SCORE risk score. Even if patients with diabetes did not present with diabetic retinopathy, the prediction rates were high with the four algorithms for SIVA, OCT-A alone, and SIVA and OCT-A combined. These findings demonstrate that vascular changes are significant before incidence of diabetic retinopathy.³⁴⁻³⁶ The present study is also valuable because of the good prediction rate not only for cardiovascular risk factors but also for well-established cardiovascular scores with Naïve Bayes and SIVA + OCT-A data: the AHA risk score 81.25%, SCORE risk 75.64%, and 96.53% for the Syntax score.

We compared the four algorithms with multiple test adjustments. Overall, Naïve Bayes and KNN gave a higher prediction rate compared to the discriminant analysis and the decision tree classification for every cardiovascular parameter. Comparison of SIVA, OCT-A alone, or combined did not demonstrate clear superiority of one technique over the others. Very different retinal vascular characteristics could be extracted from OCT-A and SIVA. With OCT-A, one could extract the fovea and capillary network features, whereas with SIVA from fundus photographs represent the analysis of venules and arterioles from the center of the optic disc and then to three successive peripheral zones. We aimed to combine both in order to give a thorough description of the retinal vascular network as a whole. Interestingly, SIVA and OCT-A combined did not give better predictive rate compared to SIVA or OCT-A alone. We could hypothesize that predictive vascular

information is already high with one single device and that extra retinal quantitative data are not incremental. Multimodal imaging could be useful in future studies to predict continuous cardiovascular parameters and not with label (two or three in this study).

At the moment, machine and deep learning algorithms are more focused on fundus photographs.¹⁰,³⁷ We showed that applying ML Bayesian classifier on OCT-A quantitative data about retinal microvascularization (solo or combined) could also be of interest in predictive models.

We acknowledge several limitations to this study. First, one should remain very cautious regarding these findings given the small size of the input data set (144 patients, 5 retinal vascular parameters on OCT-A and 15 parameters with SIVA software). Second, the robust prediction rates could be related to the definition of two or three labels for each cardiovascular parameter. Third, we only considered one eye per patient. This selection was established on the basis of the quality of images creating a selection bias within each patient. Furthermore, if retinal quantitative parameters had been available for both eyes, we could have used one eye from the participant to train the model and used the other eye for validation. In future studies, we should consider each eye individually. Fourth, very high image acquisition quality is mandatory for this type of study. The predictive model performance could be impaired by the technical limitations of each imaging device (artifacts, segmentation abnormalities, and signal strength). Jammal et al. proposed a deep learning algorithm to detect segmentation errors on OCT scans for retinal nerve fiber layer measurement.³⁸ In future studies, this kind of approach could help us to improve the quality of the images' data set. Fifth, this was a cross-sectional study. Future cardiovascular events were unknown. As a consequence, we were not able to evaluate the incremental value of the ML predictive model based on retinal parameters compared to the usual risk score. Sixth, regarding OCT-A parameters, we only collected foveal vasculature structure (foveal avascular zone, vessel inner and full, and perfusion inner and full) by means of the Angioplex software (version 10; Carl Zeiss Meditec AG). Additional OCT-A retinal biomarkers, such as fractal dimension, could in the future improve our prediction rates. Seventh, we used the same data set for training the algorithms and for testing their prediction accuracy. In future studies, we could use an external dataset to improve our accuracy. Finally, our algorithm was applied to a population with a high cardiovascular risk profile. At the moment, it could not be replicated in a healthy population for primary prevention. Furthermore, a comparison group of healthy patients
could help in future longitudinal studies to improve the algorithms' performance.

In conclusion, these preliminary findings demonstrate that ML algorithms applied to quantitative retinal vascular parameters with SIVA software and OCT-A show a good predictive performance of cardiovascular risk factors and cardiovascular risk scores. Quantitative retinal vascular biomarkers combined with an artificial intelligence strategy might be valuable data to implement a predictive model for cardiovascular parameters.

## Acknowledgments

The authors are grateful to Linda Northrup who copyedited this manuscript and to Arnaud Attye for his careful review.

The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.

## Credit authorship contribution statement

Louis Arnould: Conceptualization, investigation, and writing of the original draft.

Charles Guenancia: Conceptualization, investigation, and writing of the original draft.

Abderrahmane Bourredjem: Formal analysis and methodology.

Christine Binquet: Project administration and supervision.

Pierre-Henry Gabrielle: Writing-review and editing.

Pétra Eid: Writing-review and editing.

Florian Baudin: Writing-review and editing.

Ryo Kawasaki: Investigation.

Yves Cottin: Project administration.

Catherine Creuzot-Garcher: Project administration and supervision.

Sabir Jacquir: Formal analysis, methodology, software, and validation.

Disclosure: L. Arnould, None; C. Guenancia, None; A. Bourredjem, None; C. Binquet, None; P.-H.

Gabrielle, None; P. Eid, None; F. Baudin, None; R. Kawasaki, None; Y. Cottin, None; C. Creuzot-Garcher, None; S. Jacquir, None

## References

1. Liew G, Wang JJ, Cheung N, et al. The retinal vasculature as a fractal: methodology, reliability, and relationship to blood pressure. *Ophthalmology*. 2008;115:1951–1956.

2. Li H, Hsu W, Lee ML, Wong TY. Automatic grading of retinal vessel caliber. *IEEE Trans Biomed Eng*. 2005;52:1352–1355.

3. Spaide RF, Klancnik JM, Cooney MJ. Retinal vascular layers imaged by fluorescein angiography and optical coherence tomography angiography. *JAMA Ophthalmol*. 2015;133:45–50.

4. Hannappe M-A, Arnould L, Méloux A, et al. Vascular density with optical coherence tomography angiography and systemic biomarkers in low and high cardiovascular risk patients. *Sci Rep*. 2020;10:16718.

5. Arnould L, Guenancia C, Azemar A, et al. The EYE-MI Pilot Study: A Prospective Acute Coronary Syndrome Cohort Evaluated With Retinal Optical Coherence Tomography Angiography. *Invest Ophthalmol Vis Sci*. 2018;59:4299–4306.

6. Arnould L, Binquet C, Guenancia C, et al. 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.

7. Alan G, Guenancia C, Arnould L, et al. Retinal Vascular Density as A Novel Biomarker of Acute Renal Injury after Acute Coronary Syndrome. *Sci Rep*. 2019;9:1–9.

8. Medeiros FA, Jammal AA, Thompson AC. From Machine to Machine: An OCT-Trained Deep Learning Algorithm for Objective Quantification of Glaucomatous Damage in Fundus Photographs. *Ophthalmology*. 2019;126:513–521.

9. Burlina PM, Joshi N, Pekala M, Pacheco KD, Freund DE, Bressler NM. Automated Grading of Age-Related Macular Degeneration From Color Fundus Images Using Deep Convolutional Neural Networks. *JAMA Ophthalmol*. 2017;135:1170–1176.

10. Gulshan V, Peng L, Coram M, et al. Development and Validation of a Deep Learning Algorithm for Detection of Diabetic Retinopathy in Retinal Fundus Photographs. *JAMA*. 2016;316:2402–2410.
11. Ting DSW, Wu W-C, Toth C. Deep learning for retinopathy of prematurity screening. *Br J Ophthalmol*. 2019;103:577-579.

12. Poplin R, Varadarajan AV, Blumer K, et al. Prediction of cardiovascular risk factors from retinal fundus photographs via deep learning. *Nature Biomedical Eng*. 2018;2:158.

13. Chen D, Liu S, Kingsbury P, et al. Deep learning and alternative learning strategies for retrospective real-world clinical data. *NPJ Digit Med*. 2019. Available at: https://www.ncbi.nlm.nih.gov/pmc/articles/PMC6550223/.

14. Cardiovascular diseases (CVDs). Available at https://www.who.int/news-room/fact-sheets/detail/cardiovascular-diseases-(cvds).

15. Li M, Yang Y, Jiang H, et al. Retinal Microvascular Network and Microcirculation Assessments in High Myopia. *Am J Ophthalmol*. 2017;174:56-67.

16. von Elm E, Altman DG, Egger M, et al. The Strengthening the Reporting of Observational Studies in Epidemiology (STROBE) statement: guidelines for reporting observational studies. *J Clin Epidemiol*. 2008;61:344-349.

17. Goff DC, Lloyd-Jones DM, Bennett G, et al. 2013 ACC/AHA guideline on the assessment of cardiovascular risk: a report of the American College of Cardiology/American Heart Association Task Force on Practice Guidelines. *Circulation*. 2014;129:S49-S73.

18. Kappetein AP, Dawkins KD, Mohr FW, et al. Current percutaneous coronary intervention and coronary artery bypass grafting practices for three-vessel and left main coronary artery disease. Insights from the SYNTAX run-in phase. *Eur J Cardiothorac Surg*. 2006;29:486-491.

19. Conroy RM, Pyörälä K, Fitzgerald AP, et al. Estimation of ten-year risk of fatal cardiovascular disease in Europe: the SCORE project. *Eur Heart J*. 2003;24:987-1003.

20. Tukey JW. The Problem of Multiple Comparisons. H.I. Braun. 1953. (London, England, UK: Chapman & Hall; vol. 8).

21. Bowd C, Hao J, Tavares IM, Medeiros FA, Zangwill LM, Lee T-W, et al. Bayesian machine learning classifiers for combining structural and functional measurements to classify healthy and glaucomatous eyes. *Invest Ophthalmol Vis Sci*. 2008;49:945-953.

22. Goldbaum MH, Sample PA, Chan K, et al. Comparing machine learning classifiers for diagnosing glaucoma from standard automated perimetry. *Invest Ophthalmol Vis Sci*. 2002;43:162-169.

23. Retraction: A Novel Image Recuperation Approach for Diagnosing and Ranking Retinopathy Disease Level Using Diabetic Fundus Image. *PLoS One*. 2017. Available at: https://www.ncbi.nlm.nih.gov/pmc/articles/PMC5515455/.

24. Fraccaro P, Nicolo M, Bonetto M, et al. Combining macula clinical signs and patient characteristics for age-related macular degeneration diagnosis: a machine learning approach. *BMC Ophthalmol*. 2015;15:10.

25. Null Guosheng Hu, Null Xiaojiang Peng, Null Yongxin Yang, Hospedales TM, Verbeek J. Frankenstein: Learning Deep Face Representations Using Small Data. *IEEE Trans Image Process*. 2018;27:293-303.

26. Chartrand G, Cheng PM, Vorontsov E, et al. Deep Learning: A Primer for Radiologists. *Radiographics*. 2017;37(7):2113-2131.

27. Ahmad T, Lund LH, Rao P, et al. Machine Learning Methods Improve Prognostication, Identify Clinically Distinct Phenotypes, and Detect Heterogeneity in Response to Therapy in a Large Cohort of Heart Failure Patients. *J Am Heart Assoc*. 2018;7:e008081.

28. Ambale-Venkatesh B, Yang X, Wu CO, et al. Cardiovascular Event Prediction by Machine Learning: The Multi-Ethnic Study of Atherosclerosis. *Circ Res*. 2017;121:1092-1101.

29. Alaa AM, Bolton T, Di Angelantonio E, Rudd JHF, van der Schaar M. Cardiovascular disease risk prediction using automated machine learning: A prospective study of 423,604 UK Biobank participants. *PLoS One*. 2019. Available at: https://www.ncbi.nlm.nih.gov/pmc/articles/PMC6519796/.

30. D'Agostino RB, Vasan RS, Pencina MJ, et al. General cardiovascular risk profile for use in primary care: the Framingham Heart Study. *Circulation*. 2008;117:743-753.

31. Mitani A, Huang A, Venugopalan S, et al. Detection of anaemia from retinal fundus images via deep learning. *Nature Biomedical Engineering*. 2020;4:18-27.

32. Gerrits N, Elen B, Craenendonck TV, et al. Age and sex affect deep learning prediction of cardiometabolic risk factors from retinal images. *Sci Rep*. 2020;10:9432.

33. Nature Research Bioengineering Community. A deep-learning system for the assessment of CVD risk via the measurement of retinal-vessel calibre. *Nature Research Bioengineering Community*. 2020. Available at: http://bioengineeringcommunity.nature.com/positions/a-deep-learning-system-for-the-assessment-of-
cvd-risk-via-the-measurement-of-retinal-vessel-calibre.

34. de Carlo TE, Chin AT, Bonini Filho MA, et al. Detection of microvascular changes in eyes of patients with diabetes but not clinical diabetic retinopathy using optical coherence tomography angiography. *Retina*. 2015;35:2364–2370.

35. Thompson IA, Durrani AK, Patel S. Optical coherence tomography angiography characteristics in diabetic patients without clinical diabetic retinopathy. *Eye*. 2019;33:648–652.

36. Leontidis G, Caliva F, Al-Diri B, Hunter A. Study of the retinal vascular changes between

the early stages of diabetes and first year of diabetic retinopathy. *Invest Ophthalmol Vis Sci*. 2015;56:2015.

37. Abràmoff MD, Folk JC, Han DP, et al. Automated analysis of retinal images for detection of referable diabetic retinopathy. *JAMA Ophthalmol*. 2013;131:351–357.

38. Jammal AA, Thompson AC, Ogata NG, et al. Detecting Retinal Nerve Fibre Layer Segmentation Errors on Spectral Domain-Optical Coherence Tomography with a Deep Learning Algorithm. *Sci Rep*. 2019;9:1–9.