Skip to main content

ORIGINAL RESEARCH article

Front. Neurol., 22 August 2022
Sec. Stroke

Ischemic and haemorrhagic stroke risk estimation using a machine-learning-based retinal image analysis

Updated
\nYimin Qu,,&#x;Yimin Qu1,2,3Yuanyuan Zhuo&#x;Yuanyuan Zhuo4Jack Lee,&#x;Jack Lee1,2Xingxian HuangXingxian Huang4Zhuoxin YangZhuoxin Yang4Haibo YuHaibo Yu4Jinwen ZhangJinwen Zhang4Weiqu YuanWeiqu Yuan4Jiaman WuJiaman Wu5David OwensDavid Owens6Benny Zee,
Benny Zee1,2*
  • 1Centre for Clinical Research and Biostatistics, Jockey Club School of Public Health and Primary Care, The Chinese University of Hong Kong, Hong Kong, China
  • 2Centre for Clinical Trials and Biostatistics Lab, CUHK Shenzhen Research Institute, Shenzhen, China
  • 3School of Population Medicine and Public Health, Chinese Academy of Medical Science (CAMS) and Peking Union Medical College (PUMC), Beijing, China
  • 4Shenzhen Traditional Chinese Medicine Hospital, Shenzhen, China
  • 5Affiliated Shenzhen Maternity and Child Healthcare Hospital, Southern Medical University, Shenzhen, China
  • 6Swansea University, Wales, United Kingdom

Background: Stroke is the second leading cause of death worldwide, causing a considerable disease burden. Ischemic stroke is more frequent, but haemorrhagic stroke is responsible for more deaths. The clinical management and treatment are different, and it is advantageous to classify their risk as early as possible for disease prevention. Furthermore, retinal characteristics have been associated with stroke and can be used for stroke risk estimation. This study investigated machine learning approaches to retinal images for risk estimation and classification of ischemic and haemorrhagic stroke.

Study design: A case-control study was conducted in the Shenzhen Traditional Chinese Medicine Hospital. According to the computerized tomography scan (CT) or magnetic resonance imaging (MRI) results, stroke patients were classified as either ischemic or hemorrhage stroke. In addition, a control group was formed using non-stroke patients from the hospital and healthy individuals from the community. Baseline demographic and medical information was collected from participants' hospital medical records. Retinal images of both eyes of each participant were taken within 2 weeks of admission. Classification models using a machine-learning approach were developed. A 10-fold cross-validation method was used to validate the results.

Results: 711 patients were included, with 145 ischemic stroke patients, 86 haemorrhagic stroke patients, and 480 controls. Based on 10-fold cross-validation, the ischemic stroke risk estimation has a sensitivity and a specificity of 91.0% and 94.8%, respectively. The area under the ROC curve for ischemic stroke is 0.929 (95% CI 0.900 to 0.958). The haemorrhagic stroke risk estimation has a sensitivity and a specificity of 93.0% and 97.1%, respectively. The area under the ROC curve is 0.951 (95% CI 0.918 to 0.983).

Conclusion: A fast and fully automatic method can be used for stroke subtype risk assessment and classification based on fundus photographs alone.

Background

Stroke is one of the most important causes of morbidity and mortality worldwide. It is the second leading cause of death, accounting for 6.3 million deaths in 2015 worldwide (1). Despite the decreasing trend in China for stroke prevalence since the 1990s, the absolute number of deaths and the loss of disability-adjusted life-years keep increasing (2). Stroke has become the leading cause of mortality (3), with an age-standardized mortality rate of 114.8/100 000 person-years in 2013 (4). The overall stroke burden is exceptionally high in rural areas where medical resources are limited (4).

As therapeutic options are limited, especially in rural areas, feasible and effective screening strategies are needed to identify high-risk stroke patients. Traditional methods to assess stroke risk include ultrasound, computed tomography angiography (CTA), and magnetic resonance angiography (MRA). Ultrasound can evaluate vascular stenosis and assess blood flow velocity in the carotid artery. Yet, some research has reported that carotid stenosis is not a good enough tool for stroke screening since most stroke patients do not have moderate or high stenosis that could have been detected before an incidence of stroke (5, 6). CTA and MRA can detect extensive cerebrovascular abnormalities (7). These techniques are valid with high accuracy, yet the relatively high cost and invasive quality made them impossible to be used as screening tools. Recently, there are digital solutions to assess stroke risk for the purpose of prevention (8). However, these tools were derived from the Framingham Stroke Risk Score prediction algorithm and were enhanced to include additional lifestyle risk factors shown to be important for stroke and CVD occurrence. The additional factors may be a helpful indication of risk or a response to the outcome. The advantage is that the algorithm is easy to use but the accuracy remains a question. More study is needed to find better factors to raise the predictive accuracy. Therefore, there is an urgent need for additional techniques to detect the subtle changes, ideally at an early stage before the incidence occurs, so that prevention can be considered to avoid the damage.

Retinal vessels are the only visible vessels accessible by simple fundus photography (9). They have the same embryo origin and histological structure as cerebral vessels (1013). Retinal microvascular damages can reflect damage to cerebral microvasculature and neurons (14). It provides us with a convenient way to assess cardiovascular conditions. Previous studies have demonstrated that retinal characteristics contain valuable information for stroke risk assessment and conventional clinical variables (1520). In addition, retinal microvasculature may provide adequate information to explain the underlying pathophysiological changes of various stroke subtypes (21).

In addition to finding indicators to establish a model for stroke risk estimation, identifying stroke subtypes is also vital for guiding clinical treatment and management. Ischemic stroke is due to a lack of blood flow and accounts for about 80% of strokes. Haemorrhagic stroke is due to bleeding and accounts for about 20% of strokes (22). Stroke subtyping can have different purposes. First, classifying patients is needed for therapeutic decision-making in clinical practice. An ischemic stroke may be treatable with a medication that can break down the clot, such as aspirin. While a haemorrhagic stroke may benefit from surgery (23). Haemorrhagic stroke has a much higher death rate than ischemic stroke (24). The strategies for preventing haemorrhagic and ischemic stroke are similar but not the same due to different disease pathology (25). Ischemic stroke prevention requires a comprehensive approach to the variety of stroke risk factors a patient may encounter. Similarly, prevention for haemorrhagic stroke will have to target efforts against the vascular risk factors significant in the haemorrhage's etiology. For preventing ischemic stroke, platelet antiaggregant and anticoagulant medications are usually required. In contrast, some degree of avoidance of these same medications is an issue in preventing haemorrhagic stroke (25).

This study aimed to establish risk estimation models for ischemic and haemorrhagic stroke patients and contribute to the early classification of the two-stroke subtypes with retinal characteristics.

Study subjects

The stroke cases for this study were obtained from the Shenzhen Traditional Chinese Medicine (SZTCM) Hospital. Cases were defined as ischemic stroke patients and haemorrhagic stroke patients. Control subjects included patients with hypertension, dyslipidemia, or diabetes at the same hospital. In addition, we have increased the number of control subjects with healthy subjects from the community. The participants' inclusion criteria included good health status and ability to sit on a chair for retinal image taking; having clear basic information and disease diagnosis; having appropriate blood pressure measures, blood glucose, and blood lipids. Subjects with any of the following conditions were excluded from the study: cataracts or other eye diseases that affected retinal image taking, requiring close clinical monitoring, being too weak to comply with the research, and being unable or unwilling to comply with disease examination. For this study, we have obtained clinical research ethics approval from the Shenzhen Traditional Chinese Medicine Hospital Ethics Committee (Ref. No.: K2019-005-01) and the Joint Chinese University of Hong Kong—New Territories East Cluster Clinical Research Ethics Committee (Ref. No.: 2020.093). Baseline information, including age, sex, smoking, drinking, and medical history, was collected by trained doctors upon admission. Trained nurses measured patients' height, weight, and blood pressure on the first day of admission. Body mass index (BMI) was calculated as the weight (kg) divided by the square of the height (m). The blood samples were collected in the morning. A fully automatic analyser measured blood lipid levels of total cholesterol (TC), triglycerides (TG), high-density lipoprotein (HDL) cholesterol, low-density lipoprotein (LDL) cholesterol, fasting blood glucose level, glycated hemoglobin test (HbA1c), coagulation index levels of activated partial thromboplastin time (APTT), prothrombin time (PT), fibrinogen (FIB), thrombin time (TT.); and blood routine examination results including white blood cell (WBC), hemoglobin, blood platelet (PLT), and plateletcrit (PCT). Hypertension was defined as a systolic blood pressure of ≥ 140 mmHg or diastolic blood pressure of ≥ 90 mmHg or hypertension history. Diabetes mellitus was defined as fasting plasma glucose ≥ 6.1 mmol/L or HbA1c ≥ 5.8% or having a diabetes history. Dyslipidemia was defined as total cholesterol (TC) ≥ 240 mg/L, triglyceride (TG) ≥ 200 mg/dL, HDL <40 mg/dL, low-density lipoprotein cholesterol (LDL-c) ≥ 160 mg/dL, or as patients with a history of dyslipidemia. All patients underwent detailed radiographic evaluations, including a cranial magnetic resonance imaging (MRI) scan and a duplex color Doppler ultrasound or contrast-enhanced cranial magnetic resonance imaging angiography (MRA) (26). In addition, retinal photography was taken within 2 weeks of hospital admission.

Methods

Quantitative variables were expressed as the mean ± standard deviation, and categorical variables were expressed as counts with percentages. For univariate analysis, independent t-tests were conducted to compare continuous data between groups, and the chi-square tests were conducted for categorical data analysis. A fully automatic retinal image analysis for stroke subtypes was developed using R and Matlab computer software to estimate retinal microvascular characteristics and incorporate machine-learning techniques to estimate risks of ischemic and haemorrhagic strokes. The detailed methods of the automatic retinal imaging analysis method have been reported previously for studies related to cerebral magnetic resonance imaging (2729).

The odds ratios (OR) and 95% confidence interval (CI) were reported for variables in the model. To ensure the consistency of the models and to avoid overfitting, we have conducted a 10-fold cross-validation analysis. The sensitivity, specificity, and area under the receiver operating characteristic curve (AUC of ROC) were reported for each model. The Delong method was used to compare the difference between AUCs (30). P < 0.05 was considered as statistical significance.

For the classification models, we used machine-learning and deep learning techniques. Using Matlab, we first applied a transfer net “ResNet50” convolutional neural network with retinal images as input. The outputs were features generated based on pixels associated with stroke subtype status. We also extracted the texture/fractal/spectrum-related features (such as high order spectra and fractal dimensions) associated with stroke subtypes using the automatic retinal image analysis (ARIA) algorithm written in Matlab (31). We then used the glmnet approach to select the most important subset of features based on the penalized maximum likelihood using R and Matlab. These refined features are highly associated with stroke subtypes. Finally, we translated the features extracted from the above machine-learning approaches to commonly used retinal characteristics measured from the images using ImageJ. This part of the analysis, performed with SPSS, helped enhance our understanding of retinal characteristics that contribute to the classification and identification of specific stroke subtypes.

Retinal parameters estimation

Canon non–mydriatic retinal camera (Canon-CR2) was used to capture the retinal color image using a 45° field of view centered on the fovea. The retinal characteristics measurement tools previously developed as part of the ARIA algorithm were used to estimate the parameters of retinal vessels. The retinal characteristics measurement previously developed as part of the Automatic Retinal Image Analysis (ARIA) algorithm was used to estimate the parameters of retinal vessels (31, 32). The following is a brief description:

Retinal vessel measurements: Our estimates were based on the formula developed by Knudtson et al. (33) to describe the retinal vessel measurements into the central retinal artery equivalent (CRAE) and central retinal vein equivalent (CRVE). Arteriole-venous nipping and arteriole occlusion: The sign of arteriole-venous nipping was marked as the narrowing of the venule at the crossing point of arteriole. The arteriole occlusions (Aocclusion) were presented as thread-like arterioles when the blood inside the arterioles was stopped by emboli. Hemorrhages and exudates: Status of hemorrhages and exudates were recorded as either present or absent. Hemorrhages and exudates were key determinants for the severity of diabetic retinopathy as they were found to be associated with stroke in other studies. Tortuosity: Tortuosity was assessed by visual grading of one fovea-centered and one disc-centered fundus image from each image. The grading levels for retinal arterial tortuosity were either predominantly straight arteries or mild to severe tortuosity with at least one inflection of at least one major artery. Bifurcation coefficients (BC): Bifurcation coefficient (BC) is the ratio of the sum of the cross-sectional areas of the daughter vessels of a bifurcation to that of the parent stem. The means of the bifurcation coefficient of arterioles (BCA) and venules (BCV) were used. Asymmetry of branches and bifurcation angles: Asymmetry index (AI) is the ratio of diameters of two daughter branches. The AI was calculated as AI=D1/D2, where D1 and D2 were smaller and larger branches, respectively. The mean of the three sets of AI of arterioles (Aasymmetry) and venules (Vasymmetry) was used. The angle between two daughter branches of the same branches studied in the BC was measured. The centerline of two branches was drawn, and the angle was calculated to represent the branching angle. The mean of the bifurcation angles of arterioles (Aangle), and mean of bifurcation angles of venules (Vangle) from the three sets of vessels in one retinal image were used for the analysis.

Clinical risk factors estimation

In addition to the retinal microvascular characteristics, we used machine-learning techniques to estimate important clinical risk factors and distinguish stroke subtypes. Previous studies reported several clinical characteristics differences between haemorrhagic stroke and ischemic stroke (34). For example, Zhang et al. (35) reported that ischemic stroke patients are significantly older (p < 0.001), have a higher proportion of family history of stroke (p = 0.01), obesity (p < 0.001), diabetes (p = 0.004), TIA (p = 0.017), atrial fibrillation (p = 0.002), lower level of HDL (p = 0.001), and carotid atheroma (p = 0.002). At the same time, haemorrhagic stroke patients have a higher proportion of males (p = 0.023), alcohol drinking (p = 0.003), hypertension (p = 0.003), and increased WBC (p < 0.001). Our study would use retinal images to estimate the clinical risk factors and compare the control, haemorrhagic, and ischemic stroke groups to provide further insight into the retinal image analysis for stroke subtypes classification.

Results

Seven hundred eleven patients were enrolled, including 145 ischemic stroke patients, 86 haemorrhagic patients, and 480 controls. Among the 480 controls, 123 came from the Shenzhen TCM Hospital and 357 from healthy volunteers in the community. Descriptive statistics for stroke subtype (ischemic stroke / haemorrhagic stroke) and control related to baseline information and cardiovascular risk factors are shown in Table 1. For the comparison between ischemic stroke and control groups, age, systolic and diastolic blood pressure, and hypertension were significantly higher, but the proportion of males was significantly smaller. The same pattern occurred for the haemorrhagic stroke compared to the control group, except that significantly more males were in the haemorrhagic stroke group.

TABLE 1
www.frontiersin.org

Table 1. Comparison of clinical characteristics between control, ischemic and hemorrhagic stroke.

For the retinal characteristics, CRAE and CRVE, AVR, and bifurcation coefficients were significantly smaller in both the ischemic and haemorrhagic stroke groups. The other retinal characteristics such as bifurcation angles, asymmetry, tortuosity, nipping, hemorrhages, occlusion, and exudates have significantly larger values in the stroke-subtype groups than in control (Tables 2, 3). These results show many differences in retinal characteristics among the control group and the ischemic and haemorrhagic stroke groups.

TABLE 2
www.frontiersin.org

Table 2. Comparison of retinal characteristics between control and ischemic stroke.

TABLE 3
www.frontiersin.org

Table 3. Comparison of retinal characteristics between control and hemorrhagic stroke.

Table 4 shows the risk factors for control, ischemic stroke and haemorrhagic stroke. Comparing the three groups concerning clinical characteristics estimated from retinal images is to demonstrate that the retinal images contain information for the classification of stroke subtypes based on known significant clinical variables. In our study, we found that ischemic stroke patients who are older (p = 0.001) have more diabetes (p < 0.001) and carotid atherosclerosis (p < 0.001) than haemorrhagic stroke. In addition, both haemorrhagic and ischemic strokes have significantly more males, a higher proportion of patients with hypertension, atrial fibrillation (AF), lacunar infarct, and carotid atherosclerosis.

TABLE 4
www.frontiersin.org

Table 4. The distribution of estimated Clinical risk factors for the three classes of control, Ischemic stroke and Hemorrhagic Stroke.

For the classification analysis, Figure 1 shows the flow chart for the methods. We have analyzed the classification performance between using retinal characteristics alone vs. clinical characteristics using logistic regression. Delong's method was used to compare the AUCs of models. The results show that retinal characteristics performed significantly better than clinical characteristics alone (p < 0.001). The AUC for ischemic stroke based on clinical and retinal variables was 0.88 (95% CI of 0.84, 0.92) and 0.98 (95% CI of 0.97, 0.99), respectively (Figure 2). The AUC for haemorrhagic stroke based on clinical and retinal variables were 0.91(95% CI 0.87, 0.95) and 0.98 (95% CI of 0.97, 1.00), respectively (Figure 3).

FIGURE 1
www.frontiersin.org

Figure 1. Flowchart of the method for the development of the classification model.

FIGURE 2
www.frontiersin.org

Figure 2. ROC curves of the model using only clinical characterisitics and model using only retinal characteristics for ischemic stroke risk estimation.

FIGURE 3
www.frontiersin.org

Figure 3. ROC curves of the model using only clinical characteristics and model using only retinal characterisitcs for haemorrhagic stroke risk estimation.

For the ischemic stroke classification model, the 10-fold cross-validation gives sensitivity and specificity of 91.0% and 94.8%, respectively. The area under the ROC for ischemic stroke based on the 10-fold cross-validation analysis was 0.929 (95% CI of 0.900–0.958). The box plot for the probability of ischemic stroke is shown in Figure 4. For haemorrhagic stroke, the sensitivity and specificity were 93.0% and 97.1%, respectively. The area under the ROC for haemorrhagic stroke based on the 10-fold cross-validation was 0.951 (95% CI of 0.918–0.983). The box plot for the probability of haemorrhagic stroke is shown in Figure 5. Since there is an age difference between the control and the stroke groups, we carried out further investigation. Among the 480 controls in this study, 123 came from the same SZTCM hospital and 357 from healthy volunteers in the community. The average age of the 123 controls from SZTCM hospital was 52.13, similar to the haemorrhagic stroke patients. If we only use these 123 controls as the control group, the risk estimation models also perform well. The classification model for ischemic stroke vs. control had a sensitivity of 90.63%, a specificity of 91.56%, and an AUC of 0.98. The classification model for haemorrhagic stroke vs. control had a sensitivity of 92.97%, a specificity of 85.56%, and an AUC of 0.98. This result demonstrated the robustness of the models regardless of age. Both the sensitivity and the specificity of the two-stroke subtypes have high accuracy. Still, it is also essential to know if they can discriminate the two-stroke subtypes. We demonstrated these classification models for stroke subtypes have good discrimination power using the 30% validation portion of the data. In the scatter-plot, the classification models discriminate the probability of ischemic stroke, haemorrhagic stroke and control subjects in three apparent clusters, as shown in Figure 6. The figure can be divided into four quadrants if we draw a horizontal line for the ischemic stroke (y-axis) and a vertical line for haemorrhagic stroke (x-axis). The ischemic stroke patients were clustered in the top left-hand quadrant, with a high probability of ischemic stroke but a low probability of haemorrhagic stroke. The haemorrhagic stroke patients clustered in the lower right-hand quadrant, with a high probability of haemorrhagic stroke but a low probability of ischemic stroke. The control subjects clustered around the origin in the lower left-hand quadrant, with a low probability for both ischemic and haemorrhagic strokes.

FIGURE 4
www.frontiersin.org

Figure 4. Box plot for the probability of ischemic stroke.

FIGURE 5
www.frontiersin.org

Figure 5. Box plot for the probability of haemorrhagic stroke.

FIGURE 6
www.frontiersin.org

Figure 6. Scatter plot for both ischemic and haemorrhagic strokes.

When we evaluated the stroke subtypes together, the error rates for all three groups are shown in Table 5. For ischemic stroke patients, 9/145 (6.2%) were misclassified to control, and 3/145 (2.1%) were misclassified to haemorrhagic stroke. For haemorrhagic stroke, 4/86 (4.7%) were misclassified to control, and 11/86 (12.8%) were misclassified to ischemic stroke. For the control, the error rate was 29/480 (6%), with 28/480 (5.8%) misclassified as ischemic stroke and 1/480 (0.2%) misclassified as haemorrhagic stroke.

TABLE 5
www.frontiersin.org

Table 5. Ten-fold cross-validation of all three groups.

Discussion

The performance of the risk estimation models for ischemic and haemorrhagic strokes was excellent. We have a sensitivity and specificity of 91.0% and 94.8% for ischemic stroke classification and 93.0% and 97.1% for haemorrhagic stroke classification, respectively. These results showed that retinal characteristics were highly efficient in classifying stroke subtypes. In addition, we can now evaluate both risks of stroke subtypes longitudinally to study how they are related during their development.

The retinal characteristics are of significant interest as markers of stroke since they can be directly visualized via ophthalmoscopy (27). Previous studies have shown that retinal vascular changes vary according to stroke subtypes (17, 3638). Clinical risk factors provide general associations for stroke risk estimation, but they do not classify ischemic and haemorrhagic stroke separately. However, stroke prevention and clinical management rely on accurate risk estimation and classification. Many routine preventive treatments for ischemic stroke, including antiplatelet therapy, anticoagulants, and statins, have been noted to generate a higher risk for haemorrhagic stroke. For example, aspirin for platelet therapy would increase haemorrhagic stroke risk (3942). Study on the use of clopidogrel for antiaggregant yielded a similar result (43). Anticoagulation such as warfarin for stroke prevention also increases the risk of intracerebral hemorrhage (34, 4446). As a result, knowing the risk of stroke subtypes at an early stage is highly advantageous.

We further employed the machine-learning method to estimate clinical risk factors using retinal images and showed the differences between control, ischemic stroke, and haemorrhagic stroke agree with previous literature for stroke subtypes based on clinical risk factors. This analysis demonstrated that the retinal image contains information on the clinical variables that contributed to the model classification (35). Other investigations using retinal image information for clinical application are starting to appear. For example, a recent study shows that multifractals of the retinal vessel can be used to predict pial collateral status for patients with ischemic stroke (47).

Finally, we can use the estimated stroke subtypes risks as a target for designing a health and wellness plan from a prevention point of view. Prevention trials or lifestyle intervention studies are now feasible with the retinal image analysis approach.

Limitations

There are several limitations to this study. First, we did not have a separate data set for model validation in this research. Thus, we have conducted the 10-fold cross-validation, and the results showed that the performance of the models was stable on the training data set and the cross-validation data set. Second, the sample size was relatively small, which may affect the statistical power of the classification. Third, the haemorrhagic stroke cases in this study are mainly intracerebral hemorrhages. We have no subarachnoid hemorrhage sample in this study. Therefore, the classification may not apply to subarachnoid haemorrhagic stroke cases. Finally, since this is a case-control study, we cannot establish the temporal relationship if the retinal changes before the onset of the stroke.

Future direction

This study is a pioneer study with a potential future clinical application where we can apply the results from retinal imaging to the hospital Accident and Emergency (A&E) Department as a screening tool for stroke risk including subtypes classification. The management of ischemic and hemorrhagic strokes is very different. Retinal imaging is fast and convenient, it will provide crucial information for the A&E operation and help prioritize patients' specific needs for CT or MRI confirmation in a timely fashion.

Data availability statement

The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.

Ethics statement

The studies involving human participants were reviewed and approved by Joint CUHK-NTEC Clinical Research Ethics Committee. The patients/participants provided their written informed consent to participate in this study.

Author contributions

YQ and BZ contributed to the study design and writing of the initial draft of the manuscript. YQ, JZ, and JW contributed to the data collection. YZ, ZY, and HY contributed to data collection and their interpretation. DO and JW provided clinical advice and interpretation. JL and BZ carried out methodological development and retinal image analysis. JL and YQ provided statistical analysis and prepared figures and tables. All authors reviewed the manuscript, made significant contributions, and approved the submitted version.

Funding

This study was supported by the General Research Fund (GRF) of the Hong Kong Research Grant Council (No. 14139116); National Natural Science Foundation of China (No. 81803952); Science Technology and Innovation Commission of Shenzhen Municipality (KCXFZ20201221173208024).

Conflict of interest

Authors BZ and JL have a patent “Method and device for retinal image analysis” licensed to Health View Bioanalytic Limited and received royalties through the Chinese University of Hong Kong. BZ and JL are founders and shareholders of Health View Bioanalytic Limited, Bioanalytic Holdings Limited, and Bioanalytic International Holdings Limited.

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

Publisher's note

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

References

1. Townsend N, Bhatnagar P, Wilkins E, Wickramasinghe K, Rayner M. Cardiovascular Disease Statistics 2015. United Kingdom: Biritish Heart Foundation. (2015).

Google Scholar

2. Naghavi M, Abajobir AA, Cristiana A, Abbas K, Abd–Allah F, Abera SF, et al. Global, regional, and national age-sex specific mortality for 264 causes of death, 1980–2016: a systematic analysis for the global burden of disease study 2016. Lancet. (2017) 390:1151–210. doi: 10.1016/S0140-6736(17)32152-9

PubMed Abstract | CrossRef Full Text | Google Scholar

3. Global regional and and national age–sex specific all-cause and cause-specific mortality for 240 causes of death 1990–2013: 1990–2013: a systematic analysis for the global burden of disease study 2013. Lancet. (2015) 385:117–71. doi: 10.1016/S0140-6736(14)61682-2

PubMed Abstract | CrossRef Full Text | Google Scholar

4. Wang W, Jiang B, Sun H, Ru X, Sun D, Wang L, et al. Prevalence, incidence, and mortality of stroke in china: results from a nationwide population-based survey of 480 687 adults. Circulation. (2017) 135:759–71. doi: 10.1161/CIRCULATIONAHA.116.025250

PubMed Abstract | CrossRef Full Text | Google Scholar

5. Smith-Bindman R, Bibbins-Domingo K, USPSTF. Recommendations for screening for carotid stenosis to prevent stroke—the need for more data. JAMA Netw Open. (2021) 4:e2036218–18. doi: 10.1001/jamanetworkopen.2020.36218

PubMed Abstract | CrossRef Full Text | Google Scholar

6. Naylor AR. Why is the management of asymptomatic carotid disease so controversial? Surgeon. (2015) 13:34–43. doi: 10.1016/j.surge.2014.08.004

PubMed Abstract | CrossRef Full Text | Google Scholar

7. Yew KS, Cheng E. Acute stroke diagnosis. Am Fam Physician. (2009) 80:33–40.

PubMed Abstract | Google Scholar

8. Feigin VL, Owolabi M, Hankey GJ, Pandian J, Martins SC. Digital health in primordial and primary stroke prevention: a systematic review. Stroke. (2022) 53:1008–19. doi: 10.1161/STROKEAHA.121.036400

PubMed Abstract | CrossRef Full Text | Google Scholar

9. Semerád L, Drahanský M. Retinal Vascular Characteristics. In: Uhl A, Busch C, Marcel S, et al., editor. Handbook of Vascular Biometrics. Cham: Springer International Publishing (2020). p. 309–54.

Google Scholar

10. Patton N, Aslam T, Macgillivray T, Pattie A, Deary IJ, Dhillon B, et al. Retinal vascular image analysis as a potential screening tool for cerebrovascular disease: a rationale based on homology between cerebral and retinal microvasculatures. J Anat. (2005) 206:319–48. doi: 10.1111/j.1469-7580.2005.00395.x

PubMed Abstract | CrossRef Full Text | Google Scholar

11. Risau W. Mechanisms of angiogenesis. Nature. (1997) 386:671–74. doi: 10.1038/386671a0

PubMed Abstract | CrossRef Full Text | Google Scholar

12. Hughes S, Yang H, Chan-Ling T. Vascularization of the human fetal retina: roles of vasculogenesis and angiogenesis. Invest Ophthal Visual Sci. (2000) 41:1217–28.

PubMed Abstract | Google Scholar

13. Delaey C, Van de Voorde J. Regulatory mechanisms in the retinal and choroidal circulation. Ophthalmic Res. (2000) 32:249–56. doi: 10.1159/000055622

PubMed Abstract | CrossRef Full Text | Google Scholar

14. Zhang Y, Shi C, Chen Y, Wang W, Huang S, Han Z, et al. Retinal structural and microvascular alterations in different acute ischemic stroke subtypes. J Ophthalmol. (2020) 2020:8850309–09. doi: 10.1155/2020/8850309

PubMed Abstract | CrossRef Full Text | Google Scholar

15. Wong TY, Klein R, Couper DJ, Cooper LS, Shahar E, Hubbard LD, et al. Retinal microvascular abnormalities and incident stroke: the atherosclerosis risk in communities study. Lancet. (2001) 358:1134–40. doi: 10.1016/S0140-6736(01)06253-5

PubMed Abstract | CrossRef Full Text | Google Scholar

16. Mitchell P, Wang JJ, Wong TY, Smith W, Klein R, Leeder SR, et al. Retinal microvascular signs and risk of stroke and stroke mortality. Neurology. (2005) 65:1005–09. doi: 10.1212/01.wnl.0000179177.15900.ca

PubMed Abstract | CrossRef Full Text | Google Scholar

17. Lindley RI, Wang JJ, Wong M-C, et al. Retinal microvasculature in acute lacunar stroke: a cross-sectional study. Lancet Neurol. (2009) 8:628–34. doi: 10.1016/S1474-4422(09)70131-0

PubMed Abstract | CrossRef Full Text | Google Scholar

18. Ikram MK, de Jong FJ, Bos MJ, Vingerling JR, Hofman A, Koudstaal PJ, et al. Retinal vessel diameters and risk of stroke. Rotterdam Study. (2006) 66:1339–43. doi: 10.1212/01.wnl.0000210533.24338.ea

PubMed Abstract | CrossRef Full Text | Google Scholar

19. Yin Wong T. Is retinal photography useful in the measurement of stroke risk? Lancet Neurol. (2004) 3:179–83. doi: 10.1016/S1474-4422(04)00682-9

PubMed Abstract | CrossRef Full Text | Google Scholar

20. Klein, R. Klein, BEK. Jensen, SC. Moss, SE. Meuer, SM. Retinal, emboli and stroke: the beaver dam eye study. Arch Ophthal. (1999) 117:1063–68. doi: 10.1001/archopht.117.8.1063

PubMed Abstract | CrossRef Full Text | Google Scholar

21. Lindley RI. Retinal microvascular signs: a key to understanding the underlying pathophysiology of different stroke subtypes? Inter J Stroke: Off J Inter Stroke Soc. (2008) 3:297–305. doi: 10.1111/j.1747-4949.2008.00215.x

PubMed Abstract | CrossRef Full Text | Google Scholar

22. Virani SS, Alonso A, Benjamin EJ, Bittencourt MS, Callaway CW, Carson AP, et al. Heart disease and stroke statistics-2020 update: a report from the american heart association. Circulation. (2020) 141:e139–596. doi: 10.1161/CIR.0000000000000757

PubMed Abstract | CrossRef Full Text | Google Scholar

23. Donnan GA, Fisher M, Macleod M, Davis SM, et al. Stroke. Lancet. (2008) 371:1612–23. doi: 10.1016/S0140-6736(08)60694-7

PubMed Abstract | CrossRef Full Text | Google Scholar

24. Andersen KK, Olsen TS, Dehlendorff C, Kammersgaard LP, et al. Hemorrhagic and ischemic strokes compared. Stroke. (2009) 40:2068–72. doi: 10.1161/STROKEAHA.108.540112

PubMed Abstract | CrossRef Full Text | Google Scholar

25. Fisher M, Vasilevko V, Cribbs DH. Mixed cerebrovascular disease and the future of stroke prevention. Transl Stroke Res. (2012) 3(Suppl. 1):39–51. doi: 10.1007/s12975-012-0185-6

PubMed Abstract | CrossRef Full Text | Google Scholar

26. Brott T, Adams HP, Olinger CP, Marler JR, Barsan WG, Biller J, et al. Measurements of acute cerebral infarction: a clinical examination scale. Stroke. (1989) 20:864–70. doi: 10.1161/01.STR.20.7.864

PubMed Abstract | CrossRef Full Text | Google Scholar

27. Moss HE. Retinal vascular changes are a marker for cerebral vascular diseases. Curr Neurol Neurosci Rep. (2015) 15:40–40. doi: 10.1007/s11910-015-0561-1

PubMed Abstract | CrossRef Full Text | Google Scholar

28. Lau A, Mok V, Lee J, Fan Y, Zeng J, Lam B., et al. “Retinal image analytics detects white matter hyperintensities in healthy adults”. Ann Clin Transl Neurol. (2019) 6:98–105. doi: 10.1002/acn3.688

PubMed Abstract | CrossRef Full Text | Google Scholar

29. Zee B, Wong Y, Lee J, Fan Y, Zeng J, Lam B, et al. “Machine-learning method for localisation of cerebral white matter hyperintensities in healthy adults based on retinal images”. Brain Commun. (2021) 3doi: 10.1093/braincomms/fcab124

PubMed Abstract | CrossRef Full Text | Google Scholar

30. DeLong ER, DeLong DM, Clarke-Pearson DL. Comparing the areas under two or more correlated receiver operating characteristic curves: a nonparametric approach. Biometrics. (1988) 44:837–45. doi: 10.2307/2531595

PubMed Abstract | CrossRef Full Text | Google Scholar

31. Zee B, Lee J, Li Q. Method and Device for Retinal Image Analysis. US Patent No US8787638 B2. Washington, DC: U.S. Patent and Trademark Office (2014).

Google Scholar

32. Zee B, Lee J, Li Q, Mok V, Li Q, Kong A, et al. “Stroke risk assessment for the community by automatic retinal image analysis using fundus photograph. Qual Prim Care. (2016) 24:114–24.

Google Scholar

33. Knudtson MD, Lee KE, Hubbard LD, Wong TY, Klein R, Klein BE. Revised formulas for summarizing retinal vessel diameters. Curr Eye Res. (2003) 27:143–9. doi: 10.1076/ceyr.27.3.143.16049

PubMed Abstract | CrossRef Full Text | Google Scholar

34. Morotti A, Goldstein JN. Anticoagulant-associated intracerebral hemorrhage. Brain Hemorrhages. (2020) 1:89–94. doi: 10.1016/j.hest.2020.01.001

CrossRef Full Text | Google Scholar

35. Zhang J, Wang Y, Wang GN, Sun H, Sun T, Shi JQ, et al. Clinical factors in patients with ischemic versus hemorrhagic stroke in East China. World J Emerg Med. (2011) 2:18–23. doi: 10.5847/wjem.j.1920-8642.2011.01.003

PubMed Abstract | CrossRef Full Text | Google Scholar

36. Baker ML, Hand PJ, Liew G, Wong TY, Rochtchina E, Mitchell P, et al. Retinal microvascular signs may provide clues to the underlying vasculopathy in patients with deep intracerebral hemorrhage. Stroke. (2010) 41:618–23. doi: 10.1161/STROKEAHA.109.569764

PubMed Abstract | CrossRef Full Text | Google Scholar

37. Cheung N, Liew G, Lindley RI, Liu EY, Wang JJ, Hand P, et al. Retinal fractals and acute lacunar stroke. Ann Neurol. (2010) 68:107–11. doi: 10.1002/ana.22011

PubMed Abstract | CrossRef Full Text | Google Scholar

38. Doubal FN, MacGillivray TJ, Patton N, Dhillon B, Dennis MS, Wardlaw JM. Fractal analysis of retinal vessels suggests that a distinct vasculopathy causes lacunar stroke. Neurology. (2010) 74:1102–07. doi: 10.1212/WNL.0b013e3181d7d8b4

PubMed Abstract | CrossRef Full Text | Google Scholar

39. He J, Whelton PK, Vu B, Klag MJ. Aspirin and risk of hemorrhagic stroke: a meta-analysis of randomised controlled trials. JAMA. (1998) 280:1930–5. doi: 10.1001/jama.280.22.1930

PubMed Abstract | CrossRef Full Text | Google Scholar

40. Ge L, Niu G, Han X, Gao Y, Wu Q, Wu H, et al. Aspirin treatment increases the risk of cerebral microbleeds. Can J Neurol Sci. (2011) 38:863–8. doi: 10.1017/S0317167100012440

PubMed Abstract | CrossRef Full Text | Google Scholar

41. Vernooij MW, Haag MD, van der Lugt A, Hofman A, Krestin GP, Stricker BH, et al. Use of antithrombotic drugs and the presence of cerebral microbleeds: the Rotterdam Scan Study. Arch Neurol. (2009) 66:714–20. doi: 10.1001/archneurol.2009.42

PubMed Abstract | CrossRef Full Text | Google Scholar

42. Biffi A, Halpin A, Towfighi A, Gilson A, Busl K, Rost N, et al. Aspirin and recurrent intracerebral hemorrhage in cerebral amyloid angiopathy. Neurology. (2010) 75:693–98. doi: 10.1212/WNL.0b013e3181eee40f

PubMed Abstract | CrossRef Full Text | Google Scholar

43. CAPRIE Steering Committee. A randomised, blinded, trial of clopidogrel versus aspirin in patients at risk of ischaemic events (CAPRIE). Lancet. (1996) 348:1329–39. doi: 10.1016/S0140-6736(96)09457-3

PubMed Abstract | CrossRef Full Text | Google Scholar

44. Ueno H, Naka H, Ohshita T, Kondo K, Nomura E, Ohtsuki T, et al. Association between cerebral microbleeds on T2*-weighted MR images and recurrent hemorrhagic stroke in patients treated with warfarin following ischemic stroke. AJNR Am J Neuroradiol. (2008) 29:1483–6. doi: 10.3174/ajnr.A1120

PubMed Abstract | CrossRef Full Text | Google Scholar

45. Lee S-H, Ryu W-S, Roh J-K. Cerebral microbleeds are a risk factor for warfarin-related intracerebral hemorrhage. Neurology. (2009) 72:171–76. doi: 10.1212/01.wnl.0000339060.11702.dd

PubMed Abstract | CrossRef Full Text | Google Scholar

46. Lovelock CE, Cordonnier C, Naka H, Al-Shahi Salman R, Sudlow CL., Edinburgh Stroke Study Group, et al. Antithrombotic drug use, cerebral microbleeds, and intracerebral hemorrhage: a systematic review of published and unpublished studies. Stroke. (2010) 41:1222–28. doi: 10.1161/STROKEAHA.109.572594

PubMed Abstract | CrossRef Full Text | Google Scholar

47. Khan A, De Boever P, Gerrits N, Akhtar N, Saqqur M, Ponirakis G, et al. Retinal vessel multifractals predict pial collateral status in patients with acute ischemic stroke. PLoS ONE. (2022) 17:e0267837. doi: 10.1371/journal.pone.0267837

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: stroke subtypes classification, retinal image analysis, ischemic stroke, haemorrhagic stroke, machine-learning method

Citation: Qu Y, Zhuo Y, Lee J, Huang X, Yang Z, Yu H, Zhang J, Yuan W, Wu J, Owens D and Zee BC-Y (2022) Ischemic and haemorrhagic stroke risk estimation using a machine-learning-based retinal image analysis. Front. Neurol. 13:916966. doi: 10.3389/fneur.2022.916966

Received: 15 April 2022; Accepted: 01 August 2022;
Published: 22 August 2022.

Edited by:

Leonardo Roever MHS, Federal University of Uberlandia, Brazil

Reviewed by:

Rayaz A. Malik, Weill Cornell Medicine- Qatar, Qatar
JoonNyung Heo, Severance Hospital, South Korea

Copyright © 2022 Qu, Zhuo, Lee, Huang, Yang, Yu, Zhang, Yuan, Wu, Owens and Zee. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Benny Zee, bzee@cuhk.edu.hk

These authors have contributed equally to this work

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