Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

An efficient approach to estimate the risk of coronary artery disease for people living with HIV using machine-learning-based retinal image analysis

  • Grace Lui,

    Roles Conceptualization, Data curation, Formal analysis, Funding acquisition, Investigation, Methodology, Project administration, Resources, Supervision, Writing – original draft, Writing – review & editing

    Affiliations Department of Medicine and Therapeutics, Prince of Wales Hospital, The Chinese University of Hong Kong, Shatin, Hong Kong SAR, Stanley Ho Centre for Emerging Infectious Diseases, The Chinese University of Hong Kong, Shatin, Hong Kong SAR

  • Ho Sang Leung,

    Roles Data curation, Investigation, Writing – review & editing

    Affiliation Department of Imaging and Interventional Radiology, Prince of Wales Hospital, Shatin, Hong Kong SAR

  • Jack Lee,

    Roles Formal analysis, Methodology, Software, Validation, Writing – review & editing

    Affiliation Centre for Clinical Research and Biostatistics, The Jockey Club School of Public Health and Primary Care, The Chinese University of Hong Kong, Shatin, Hong Kong SAR

  • Chun Kwok Wong,

    Roles Investigation, Writing – review & editing

    Affiliation Department of Chemical Pathology, The Chinese University of Hong Kong, Shatin, Hong Kong SAR

  • Xinxin Li,

    Roles Formal analysis

    Affiliation Centre for Clinical Research and Biostatistics, The Jockey Club School of Public Health and Primary Care, The Chinese University of Hong Kong, Shatin, Hong Kong SAR

  • Mary Ho,

    Roles Investigation, Writing – review & editing

    Affiliation Department of Ophthalmology, Prince of Wales Hospital, Shatin, Hong Kong SAR

  • Vivian Wong,

    Roles Data curation, Project administration

    Affiliation Department of Medicine and Therapeutics, Prince of Wales Hospital, The Chinese University of Hong Kong, Shatin, Hong Kong SAR

  • Timothy Li,

    Roles Data curation, Project administration

    Affiliation Department of Medicine and Therapeutics, Prince of Wales Hospital, The Chinese University of Hong Kong, Shatin, Hong Kong SAR

  • Tracy Ho,

    Roles Data curation

    Affiliation Department of Medicine and Therapeutics, Prince of Wales Hospital, The Chinese University of Hong Kong, Shatin, Hong Kong SAR

  • Yin Yan Chan,

    Roles Data curation

    Affiliation Department of Medicine and Therapeutics, Prince of Wales Hospital, The Chinese University of Hong Kong, Shatin, Hong Kong SAR

  • Shui Shan Lee,

    Roles Investigation, Writing – review & editing

    Affiliation Stanley Ho Centre for Emerging Infectious Diseases, The Chinese University of Hong Kong, Shatin, Hong Kong SAR

  • Alex PW Lee,

    Roles Investigation, Writing – review & editing

    Affiliations Department of Medicine and Therapeutics, Prince of Wales Hospital, The Chinese University of Hong Kong, Shatin, Hong Kong SAR, Laboratory of Cardiac Imaging and 3D Printing, Li Ka Shing Institute of Health Science, The Chinese University of Hong Kong, Shatin, Hong Kong SAR

  • Ka Tak Wong,

    Roles Investigation, Writing – review & editing

    Affiliation Department of Imaging and Interventional Radiology, Prince of Wales Hospital, Shatin, Hong Kong SAR

  • Benny Zee

    Roles Conceptualization, Formal analysis, Funding acquisition, Methodology, Supervision, Writing – original draft, Writing – review & editing

    bzee@cuhk.edu.hk

    Affiliation Centre for Clinical Research and Biostatistics, The Jockey Club School of Public Health and Primary Care, The Chinese University of Hong Kong, Shatin, Hong Kong SAR

Abstract

Background

People living with HIV (PLWH) have increased risks of non-communicable diseases, especially cardiovascular diseases. Current HIV clinical management guidelines recommend regular cardiovascular risk screening, but the risk equation models are not specific for PLWH. Better tools are needed to assess cardiovascular risk among PLWH accurately.

Methods

We performed a prospective study to determine the performance of automatic retinal image analysis in assessing coronary artery disease (CAD) in PLWH. We enrolled PLWH with ≥1 cardiovascular risk factor. All participants had computerized tomography (CT) coronary angiogram and digital fundus photographs. The primary outcome was coronary atherosclerosis; secondary outcomes included obstructive CAD. In addition, we compared the performances of three models (traditional cardiovascular risk factors alone; retinal characteristics alone; and both traditional and retinal characteristics) by comparing the area under the curve (AUC) of receiver operating characteristic curves.

Results

Among the 115 participants included in the analyses, with a mean age of 54 years, 89% were male, 95% had undetectable HIV RNA, 45% had hypertension, 40% had diabetes, 45% had dyslipidemia, and 55% had obesity, 71 (61.7%) had coronary atherosclerosis, and 23 (20.0%) had obstructive CAD. The machine-learning models, including retinal characteristics with and without traditional cardiovascular risk factors, had AUC of 0.987 and 0.979, respectively and had significantly better performance than the model including traditional cardiovascular risk factors alone (AUC 0.746) in assessing coronary artery disease atherosclerosis. The sensitivity and specificity for risk of coronary atherosclerosis in the combined model were 93.0% and 93.2%, respectively. For the assessment of obstructive CAD, models using retinal characteristics alone (AUC 0.986) or in combination with traditional risk factors (AUC 0.991) performed significantly better than traditional risk factors alone (AUC 0.777). The sensitivity and specificity for risk of obstructive CAD in the combined model were 95.7% and 97.8%, respectively.

Conclusion

In this cohort of Asian PLWH at risk of cardiovascular diseases, retinal characteristics, either alone or combined with traditional risk factors, had superior performance in assessing coronary atherosclerosis and obstructive CAD.

Summary

People living with HIV in an Asian cohort with risk factors for cardiovascular disease had a high prevalence of coronary artery disease (CAD). A machine-learning-based retinal image analysis could increase the accuracy in assessing the risk of coronary atherosclerosis and obstructive CAD.

Background

Effective and durable anti-retroviral therapy allows us to witness tremendous improvement in life expectancy in people living with HIV (PLWH). However, substantial morbidity and mortality in PLWH are due to non-communicable diseases, such as cardiovascular diseases.

PLWH had a two-fold increased risk of cardiovascular diseases [1]. The global burden of cardiovascular diseases attributable to HIV has tripled since the 1990s, with sub-Saharan Africa and Asia-Pacific regions being the most affected areas [1]. Despite significant improvements in the management of HIV and its comorbidities over the last two decades. Recent population-based studies consistently showed that PLWH had higher cardiovascular and cerebrovascular disease risks [2, 3]. PLWH also carries a higher prevalence of risk factors for cardiovascular diseases, such as hypertension and diabetes [2]. Studies performed in Asia also showed that PLWH had high risks of coronary artery disease (CAD), partly attributable to the high prevalence of cardiovascular disease risk factors, suboptimal screening, and suboptimal management of these risk factors, especially in general low- and middle-income settings [4].

Current HIV clinical management guidelines recommend regular cardiovascular risk screening in PLWH; however, the best risk prediction model for PLWH is uncertain [5, 6]. Although various cardiovascular disease risk prediction functions are currently available. These algorithms have primarily been developed in non-HIV-infected populations. Therefore, they could not accurately predict risk in PLWH, possibly due to differences in pathogenesis underlying cardiovascular disease [7]. Moreover, many of these functions, including HIV-specific prediction models, have not been adequately validated in Asian populations of PLWH [4].

Retinal vascular characteristics have been recognized to encompass features associated with systemic diseases, such as diabetes and hypertension. Recent research has demonstrated a broader breadth of the application of retinal biomarkers for diagnostic, monitoring, and prognostic purposes in a wide range of chronic diseases [8]. Recently, retinal image characteristics have been shown to be closely linked with multiple cardiovascular risk factors and major cardiovascular events [9]. In particular, several retinal vascular characteristics, including arteriolar and venular calibre, curvature tortuosity, and branching complexity, were shown to have associations with CAD [10, 11].

Traditionally, manual interpretation of retinal images was heavily operator-dependent and time-consuming and was subjected to measurement error. Recently, the availability of automated models has shown superior accuracy and has transformed the practicality of adopting the assessment of retinal characteristics into routine clinical use [8, 12, 13]. However, most studies often involved a limited number of vascular characteristics [14]. In contrast, contemporary computing methods allow efficient and accurate measurements of a broad spectrum of retinal microvasculature characteristics. The retinal characteristics include vascular calibre, tortuosity, density, and branching complexity [15]. Retinal image analysis also has the advantages of being a simple, non-invasive procedure, requiring minimal operator training.

The application of retinal image analysis in assessing the risk of cardiovascular diseases has not yet been evaluated in PLWH. Also, HIV infection per se and associated opportunistic diseases cause various retinal abnormalities [16]. Therefore, studies showing a correlation between retinal vascular characteristics and cardiovascular diseases performed in non-HIV-infected populations might not be generalizable to PLWH. This study aimed to determine the prevalence of CAD among high-risk PLWH in a predominantly Asian population and determine the performance of machine-learning-based retinal image analysis in assessing the risk of CAD in PLWH compared to traditional risk prediction tools.

Methods

We performed a prospective study on CAD in an Asian population of PLWH with atherosclerotic risk factors at the Prince of Wales Hospital Infectious Diseases clinic in Hong Kong from February 2019 to February 2021. The primary outcome is coronary atherosclerosis; secondary outcomes include the presence of any coronary artery calcium (CAC), significant CAC, and obstructive CAD.

We enrolled PLWH aged 30 years or above with either chest pain or the presence of one or more risk factors for cardiovascular disease. Cardiovascular disease risk factors included hypertension, diabetes mellitus, dyslipidemia (defined by total cholesterol ≥6.2 mmol/L, HDL cholesterol ≤0.9 mmol/L, triglyceride ≥2.3 mmol/L, or use of the lipid-lowering drug) [17], current smoker, obesity (defined by body mass index ≥27 kg/m2) [18], and family history of CAD (defined by a first-degree relative with myocardial infarction before age 50 years) [17]. Patients with previously diagnosed CAD, creatinine clearance < 60mL/min, allergy to intravenous contrast, and pregnancy were excluded. All participants provided written informed consent. The study was approved by the Joint Chinese University of Hong Kong-New Territories East Cluster Clinical Research Ethics Committee.

We collected demographic and clinical information, including HIV-related clinical data and comorbidities. We measured body weight, height, waist circumference, and blood pressure. After at least 8 hours of fasting, we collected blood samples for glucose, HbA1c, insulin, cholesterol, triglyceride, creatinine, C reactive protein, fibrinogen, D-dimer, sCD14, sCD163, and adiponectin.

We calculated cardiovascular risk using four cardiovascular risk prediction functions for each participant. They included Framingham risk score for 10-year coronary heart disease risk [19], QRISK3 [20], Pooled Cohort ASCVD risk equations [21], and 10-year cardiovascular disease risk using the Data-collection on Adverse Effects of Anti-HIV Drugs (D:A:D) study (DAD) cohort risk prediction [22]. A prediction of <10% was considered a low risk of cardiovascular disease.

All participants underwent a coronary CT angiogram, which included a non-contrast CT for calcium scoring and subsequently contrast-enhanced CT with cardiac gating for coronary angiography. We obtained a CAC score from the non-contrast CT and quantified it using the Agatston method [23]. The presence of any CAC was defined as Agatston score > 0, and significant CAC was defined as Agatston score ≥ 100. We assessed the coronary plaque burden, including the presence, the site, composition (calcified, non-calcified, or mixed), and degree of stenosis of the plaques following the Society of Cardiovascular Computed Tomography guidelines [24]. In brief, stenosis was categorized as normal or minimal (0–25%), mild (26% - 50%), moderate (51% - 75%), severe non-subtotal occlusion (76% - 90%), and severe subtotal occlusion (91–99%). Coronary atherosclerosis was defined as the presence of any plaques in one or more coronary artery segments, while moderate to severe stenosis of the lumen was considered as obstructive CAD.

A trained research nurse acquired digital fundus photographs using a Canon CR2-AF non-mydriatic retinal camera from both eyes of each participant. We then assessed retinal characteristics, for example, retinal vessel measurements, arteriole-venous nicking, arteriole occlusion, haemorrhages, exudates, tortuosity, bifurcation coefficients, asymmetry of branches, and bifurcation angles. The definitions of the retinal parameters were previously presented in detail [11, 25].

Statistical methods

We performed these measurements using R (University of Auckland, Auckland) and Matlab (MathWorks, Massachusetts, USA) computer software. The methods included fractal analysis, high-order spectra analysis, and statistical texture analysis [25].

We presented descriptive statistics for the baseline characteristics. We compared demographic, clinical, and retinal characteristics between those with and without primary and secondary outcomes using an independent two-sample Student t-test and Mann Whitney U test for continuous variables and a chi-square test for categorical variables. Using stepwise logistic regression analyses, we determined the associations of traditional cardiovascular risk factors and retinal characteristics with the primary and secondary outcomes. Three different models were evaluated: Model 1 included clinical characteristics regarded as traditional cardiovascular risk factors in PLWH; model 2 included retinal characteristics only; model 3 included both the traditional cardiovascular risk factors and retinal characteristics. All covariates with a p-value less than 0.1 were kept in the final model. We compared the performances of different models by comparing the area under the curve (AUC) of receiver operating characteristic (ROC) curves using the Delong method [26] The sensitivity and specificity of the models will also be calculated.

For the classification analysis, we used machine learning and deep learning techniques. Using Matlab, we first applied transfer deep network ResNet50 convolutional neural network with retinal images as input, and the outputs were features generated at the layer of ‘’fc1000_softmax’’, based on pixels associated with the specific outcome status [27]. We also extracted the texture/spectrum/fractal-based features that are associated with the specific outcome by using the automatic retinal image analysis (ARIA) algorithm written in Matlab [28]. We then used the Glmnet approach to select significant features based on the penalised maximum likelihood by using R and Matlab [29, 30]. These refined features are highly associated with the specific outcome. Finally, we translated the features extracted from the aforementioned machine learning approaches to commonly used retinal characteristics measured from the images using ImageJ. This part of the analysis helped enhance our understanding of retinal characteristics that contribute to the classification and identification of the specific outcome and was performed with SPSS. For the validation, we applied a 10-fold cross-validation method by using a support vector machine (SVM) algorithm for testing datasets that were not used in the training of the model [30, 31]. This was performed by partitioning the dataset and using a subset to train the algorithm, and the remaining subset of data for testing. Each time we ran the cross-validation analysis, we used 10% of the data for testing that were not used at all in the training data. The advantage of this method is that the data used for testing in each run were excluded from the specific training models for the purpose of validation to reduce the problem of overfitting and overestimation of the sensitivity and specificity. Because cross-validation does not use all of the data to build a model, it is a commonly used method to prevent overfitting during training. Fig 1 shows the flowchart of the described methodology.

thumbnail
Fig 1. Flowchart of the method for the development of the classification models.

https://doi.org/10.1371/journal.pone.0281701.g001

Results

This study enrolled 120 participants during the recruitment period. Five participants were excluded from the analyses. Among them, 3 participants were excluded due to the unavailability of good quality retinal images, 1 participant was deceased, and 1 participant was lost to follow-up after recruitment without completing all study procedures. Among the 115 participants included in the analyses, the mean (±standard deviation) age was 54±10 years, 89% were male, the median (interquartile range/IQR) duration of HIV diagnosis was 12 (7–17) years, 95% had undetectable serum HIV RNA, and the median (IQR) CD4 count was 632 (451–840) cells/mm3. Among this cohort, 45% had hypertension, 40% had diabetes, 45% had dyslipidemia, and 55% had obesity. In addition, chest pain and dyspnea were present in 17% and 18% of participants. Their detailed demographic and clinical characteristics are shown in Tables 1 and 2.

thumbnail
Table 1. Demographic and clinical characteristics in participants with and without coronary atherosclerosis1.

https://doi.org/10.1371/journal.pone.0281701.t001

thumbnail
Table 2. Demographic and clinical characteristics in participants with and without obstructive coronary artery disease (CAD)#.

https://doi.org/10.1371/journal.pone.0281701.t002

Seventy-one participants (61.7%) had coronary atherosclerosis; also, these 71 participants were found to have a presence of any CAC. Thirty-five participants (30.4%) had significant CAC, and 23 (20.0%) had obstructive CAD. Coronary atherosclerosis was associated with male gender, older age, dyslipidemia, hypertension, lower CD4:CD8 ratio, lower HDL cholesterol, and higher triglyceride (Table 1). Obstructive CAD was associated with older age, dyspnea, lower total cholesterol, lower LDL cholesterol, and higher sCD163 (Table 2). Significant CAC was associated with older age, dyslipidemia, and lower LDL cholesterol (S1 Table).

The retinal characteristics associated with coronary atherosclerosis, obstructive CAD, and significant CAC are shown in S2 and S3 Tables. Several retinal vascular characteristics were associated with coronary atherosclerosis. They include narrower arterioles, wider venules, and a lower degree of arteriolar branching,

The associations between traditional cardiovascular risk factors and retinal characteristics with coronary atherosclerosis are shown in Table 3. We further improve the classification using a machine-learning approach. The model, including traditional cardiovascular risk factors and retinal characteristics, had 93.0% sensitivity and 93.2% specificity (area under the ROC curve (AUC) was 0.987, 95% CI 0.973–1.00). The performance is similar to the model with retinal characteristics alone, with 91.5% sensitivity and 88.6% specificity (AUC 0.979, 95% CI 0.960–0.998). However, both models with retinal characteristics were better than the model with traditional cardiovascular risk factors alone, with 81.7% sensitivity and 54.5% specificity (AUC 0.746, 95% CI 0.652–0.841) in assessing the risk of coronary atherosclerosis (Fig 2A).

thumbnail
Fig 2.

a. AUC comparison between different models for the prediction of coronary atherosclerosis. (Model 1: Inclusion of traditional cardiovascular risk factors only; Model 2: Inclusion of retinal information only; Model 3: Inclusion of traditional cardiovascular risk factors and retinal information). b. AUC comparison between different models for the prediction of obstructive coronary artery disease. (Model 1: Inclusion of traditional cardiovascular risk factors only; Model 2: Inclusion of retinal information only; Model 3: Inclusion of traditional cardiovascular risk factors and retinal information).

https://doi.org/10.1371/journal.pone.0281701.g002

thumbnail
Table 3. Stepwise logistic regression analysis showing factors associated with coronary atherosclerosis in different models.

https://doi.org/10.1371/journal.pone.0281701.t003

The models showing the associations between traditional cardiovascular risk factors plus retinal characteristics for obstructive CAD are shown in Table 4. For assessing the risk of obstructive CAD, the model including retinal variables combined with traditional risk factors had a sensitivity and specificity of 95.7% and 97.8%, respectively (AUC 0.991, 95% CI 0.978–1.000). The model with retinal characteristics alone had sensitivity and specificity of 87.0% and 96.7%, respectively (AUC 0.986, 95% CI 0.970–1.000). Both models, including retinal characteristics, performed significantly better than the model with traditional risk factors alone and had sensitivity and specificity of 34.3% and 93.5%, respectively (AUC 0.777, 95% CI 0.672–0.883) (Fig 2B).

thumbnail
Table 4. Stepwise logistic regression analysis showing factors associated with obstructive coronary artery disease in different models.

https://doi.org/10.1371/journal.pone.0281701.t004

Likewise, models including retinal characteristics alone or those combined with traditional cardiovascular risk factors performed significantly better than those with traditional risk factors alone in assessing the risk of significant CAC (S4 Table and S1 Fig). An example of retinal vascular characteristics and their related measurements is shown in Fig 3.

thumbnail
Fig 3.

a & b show the retinal images of two participants with coronary artery disease. White arrow indicates arteriole occlusion. Black arrow indicates a venous vessel with significant vascular tortuosity. The dashed white arrow indicates a venous bifurcation; θ1 and θ2 specify the bifurcation coefficients of the two branching venules (MBCV).

https://doi.org/10.1371/journal.pone.0281701.g003

All conventional cardiovascular disease risk prediction functions had limited performance in assessing the risk of coronary atherosclerosis and obstructive CAD. Among those participants with coronary atherosclerosis, 45.1%, 57.7%, 52.1% and 42.3% were categorized as having low risk using Framingham risk score, QRISK3, Pooled Cohort ASCVD risk equations, and 10-year DAD cohort risk prediction, respectively. The corresponding figures for obstructive CAD among these groups were 30.4%, 52.2%, 34.8% and 34.8%, respectively. Moreover, ROC curve analyses showed that all of the prediction functions had AUC <0.7, with QRISK3 having the highest AUC for assessing both coronary atherosclerosis (AUC 0.667, 95% CI 0.566–0.767) and obstructive CAD (AUC 0.663, 95% CI 0.545–0.780) (S5 Table). For fairness in comparison, we used only the logistic regression model as a method for the comparison without using the machine-learning classification method. The models adopting the retinal characteristics alone (AUC 0.902, 95% CI 0.810–0.995) or in combination with the QRISK3 score (AUC 0.912, 95% CI 0.829–0.994) had a significantly better performance than QRISK3 score alone in assessing the risk of obstructive CAD (S6 and S7 Tables and S2 and S3 Figs).

Discussion

In this study, we evaluated machine-learning-based retinal imaging analysis to assess the risk of CAD in PLWH. In this cohort of Asian PLWH at risk of developing cardiovascular diseases, retinal characteristics, either alone or combined with traditional risk factors, enhanced the performance in assessing the risk of coronary atherosclerosis and obstructive CAD. In contrast, the performance of traditional cardiovascular risk prediction functions has a lot of room to improve.

This study showed that 62% of at-risk PLWH had coronary atherosclerosis, and 20% had obstructive CAD. The most extensive cohort study involving PLWH in Asia had recently demonstrated that traditional cardiovascular risk factors, including older age, hypertension, dyslipidemia, and high body mass index, were the major contributing factors for the development of cardiovascular disease diseases among PLWH in the region [32]. Our study further confirmed the high prevalence of CAD and its contributing risk factors among PLWH, highlighting the importance of prevention and treatment of these risk factors among PLWH. Moreover, an accurate tool to predict cardiovascular risk in these at-risk individuals is highly recommended.

We demonstrated that all of the evaluated cardiovascular risk prediction functions had limited accuracy in assessing the risks of coronary atherosclerosis and obstructive CAD in Asian PLWH. Evidence from current literature suggests that currently available cardiovascular risk prediction functions have suboptimal performance among PLWH. For example, Framingham risk scores tended to underestimate the prevalence of CAD and other atherosclerotic cardiovascular diseases among PLWH in the United States across all risk groups [33] while overestimating cardiovascular risk in European [34] and Asian [35] populations of PLWH. The available HIV-specific DAD cohort risk-prediction model was developed in a primarily European cohort [6]. It predicted a lower proportion having a high risk of cardiovascular disease among Asian PLWH than other prediction functions [3537]. However, to the best of our knowledge, it had not yet been validated in any Asian populations of PLWH. Better tools to accurately assess the risk of CAD and other cardiovascular diseases among Asian PLWH are preferred.

Currently available studies performed in general HIV-uninfected populations have demonstrated the association of retinal vascular characteristics with multiple cardiovascular risk factors, including body mass index, smoking, hypertension, diabetes, and dyslipidemia [12, 13, 38, 39]. To supplement, these studies included multi-ethnic populations, including China [12], Britain [13], and the Middle East [38, 39]. In particular, smaller arteriolar widths, larger venular widths, and increased arteriolar and venular tortuosity were associated with cardiovascular risk factors [13]. Furthermore, using deep learning on retinal images, it was able to predict the presence of CAC and stratify cardiovascular disease risk in multi-ethnic populations [40]. In addition, the retinal vascular density and vessel branching complexity were able to predict higher mortality incidents in populations with high rates of hypertension and diabetes [15].

Regarding the prediction models of CAD, longitudinal studies showed that arteriolar narrowing, venular widening and fractal dimension accurately predicted incident CAD and CAD mortality in population-based studies [14, 41]. Patients hospitalized with acute coronary syndrome had lower retinal inner vessel length density and perfusion density than controls with lower cardiovascular risk [42]. In another study from China, patients with stable CAD had lower vessel density in superficial capillary plexus and deep capillary plexus; vessel density was also associated with CAD severity [43]. To further understand the pathogenesis of retinal vessel disease, recent genome-wide association analyses identified loci associated with retinal microvascular architecture, including genes associated with inflammatory, chemokine and angiogenesis pathways [15]. Such associations suggested that the retina may provide a window for identifying pathogenic processes underlying cardiovascular diseases in both HIV-uninfected populations and PLWH.

Our study further showed that a machine-learning-based retinal image analysis could improve the accuracy of assessing CAD risk among PLWH. In addition, it enhanced the performance of several widely adopted cardiovascular risk prediction functions in PLWH. This observation further supports that the retinal vessels provide an in-vivo examination of vascular sequelae secondary to CAD risk factors, including hypertension and diabetes. In studies performed in HIV-uninfected populations, retinal image analysis, when coupled with risk prediction functions, had better performance than risk prediction models alone in cardiovascular risk stratification [40] and prediction of cardiovascular mortality [44]. Retinal image analysis can potentially enhance the current risk prediction functions in cardiovascular disease risk stratification among PLWH.

The performance of retinal image analysis was best in assessing the risk of obstructive CAD in our cohort. While retinal characteristics can act as systemic biomarkers [8], the results from our study supported that retinal characteristics would be most helpful in identifying at-risk PLWH with obstructive CAD. This group of patients may benefit from more stringent cardiovascular risk factor control and coronary interventions.

Our study has several limitations. First, we had a relatively small sample size and involved only PLWH with cardiovascular risk factors from Asia. The cross-sectional study design precluded the evaluation of prediction of incident cardiovascular events by retinal image analysis. Future studies should evaluate the performance of retinal image analysis in predicting CAD in a more diverse population of PLWH with different ethnicities and levels of cardiovascular risk. Previous studies have shown differences in vascular calibre and fractal dimension among Asian ethnic groups [45]. Finally, this study has not involved the use of optical coherence tomography angiography (OCTA), which can provide more detailed retinal and choroidal characteristics analyses, which have also been shown to be associated with CAD [46].

In conclusion, we have demonstrated that machine-learning-based retinal image analysis accurately assesses the risk of coronary atherosclerosis and obstructive CAD among PLWH with risk factors for cardiovascular diseases. This tool should be further validated in more diverse populations of PLWH for the adoption in clinical practice for CAD risk stratification.

Supporting information

S1 Table. Demographic and clinical characteristics in participants with and without significant coronary artery calcium.

https://doi.org/10.1371/journal.pone.0281701.s001

(TIF)

S2 Table. Retinal characteristics in participants with and without coronary atherosclerosis or obstructive CAD.

https://doi.org/10.1371/journal.pone.0281701.s002

(TIF)

S3 Table. Retinal characteristics in participants with and without significant coronary artery calcium.

https://doi.org/10.1371/journal.pone.0281701.s003

(TIF)

S4 Table. Stepwise logistic regression analysis shows factors associated with significant coronary artery calcium in different models.

https://doi.org/10.1371/journal.pone.0281701.s004

(TIFF)

S5 Table. Area under curve (AUC) of cardiovascular risk prediction functions for coronary atherosclerosis and obstructive coronary artery disease.

https://doi.org/10.1371/journal.pone.0281701.s005

(TIFF)

S6 Table. Stepwise logistic regression analysis showing risk prediction function and retinal characteristics associated with coronary atherosclerosis in different models.

https://doi.org/10.1371/journal.pone.0281701.s006

(TIFF)

S7 Table. Stepwise logistic regression analysis showing risk prediction function and retinal characteristics associated with obstructive coronary artery disease in different models.

https://doi.org/10.1371/journal.pone.0281701.s007

(TIFF)

S1 Fig. Machine-learning analysis for AUC comparison between different models in the prediction of significant CAC (Model 1: Inclusion of traditional cardiovascular risk factors only; Model 2: Inclusion of retinal information only; Model 3: Inclusion of traditional cardiovascular risk factors and retinal information).

https://doi.org/10.1371/journal.pone.0281701.s008

(TIFF)

S2 Fig. Machine-learning analysis for AUC comparison between different models including risk prediction function and retinal information in predicting coronary atherosclerosis.

(Model 1: Inclusion of cardiovascular risk prediction function with the highest AUC only; Model 2: Inclusion of retinal information only; Model 3: Inclusion of cardiovascular risk prediction score and retinal characteristics).

https://doi.org/10.1371/journal.pone.0281701.s009

(TIFF)

S3 Fig. Machine-learning analysis for AUC comparison between different models including risk prediction score and retinal information in predicting obstructive coronary artery disease.

(Model 1: Inclusion of cardiovascular risk prediction score with the highest AUC only; Model 2: Inclusion of retinal information only; Model 3: Inclusion of cardiovascular risk prediction score and retinal characteristics).

https://doi.org/10.1371/journal.pone.0281701.s010

(TIFF)

References

  1. 1. Shah ASV, Stelzle D, Lee KK, Beck EJ, Alam S, Clifford S, et al. Global Burden of Atherosclerotic Cardiovascular Disease in People Living with the Human Immunodeficiency Virus: A Systematic Review and Meta-Analysis. Circulation. 2018.
  2. 2. Gooden TE, Gardner M, Wang J, Jolly K, Lane DA, Benjamin LA, et al. Incidence of cardiometabolic diseases in people living with and without HIV in the UK: a population-based matched cohort study. J Infect Dis. 2021.
  3. 3. Alonso A, Barnes AE, Guest JL, Shah A, Shao IY, Marconi V. HIV Infection and Incidence of Cardiovascular Diseases: An Analysis of a Large Healthcare Database. Journal of the American Heart Association. 2019;8(14):e012241. pmid:31266386
  4. 4. Choi JY, Lui GCY, Liao CT, Yang CJ. Managing cardiovascular risk in people living with HIV in Asia—where are we now? HIV Med. 2021.
  5. 5. Ryom L, Cotter A, De Miguel R, Béguelin C, Podlekareva D, Arribas J, et al. 2019 update of the European AIDS Clinical Society Guidelines for treatment of people living with HIV version 10.0. HIV Medicine.n/a(n/a).
  6. 6. Feinstein MJ, Hsue PY, Benjamin LA, Bloomfield GS, Currier JS, Freiberg MS, et al. Characteristics, Prevention, and Management of Cardiovascular Disease in People Living With HIV: A Scientific Statement From the American Heart Association. Circulation. 2019;140(2):e98–e124. pmid:31154814
  7. 7. D’Agostino RB, Sr. Cardiovascular risk estimation in 2012: lessons learned and applicability to the HIV population. J Infect Dis. 2012;205 Suppl 3(Suppl 3):S362–7.
  8. 8. Ranchod TM. Systemic retinal biomarkers. Current opinion in ophthalmology. 2021;32(5):439–44. pmid:34231533
  9. 9. Poplin R, Varadarajan AV, Blumer K, Liu Y, McConnell MV, Corrado GS, et al. Prediction of cardiovascular risk factors from retinal fundus photographs via deep learning. Nature biomedical engineering. 2018;2(3):158–64. pmid:31015713
  10. 10. Wang SB, Mitchell P, Liew G, Wong TY, Phan K, Thiagalingam A, et al. A spectrum of retinal vasculature measures and coronary artery disease. Atherosclerosis. 2018;268:215–24. pmid:29050745
  11. 11. Guo VY, Chan JC, Chung H, Ozaki R, So W, Luk A, et al. Retinal Information is Independently Associated with Cardiovascular Disease in Patients with Type 2 diabetes. Sci Rep. 2016;6:19053. pmid:26754623
  12. 12. Zhang L, Yuan M, An Z, Zhao X, Wu H, Li H, et al. Prediction of hypertension, hyperglycemia and dyslipidemia from retinal fundus photographs via deep learning: A cross-sectional study of chronic diseases in central China. PLoS One. 2020;15(5):e0233166. pmid:32407418
  13. 13. Owen CG, Rudnicka AR, Welikala RA, Fraz MM, Barman SA, Luben R, et al. Retinal Vasculometry Associations with Cardiometabolic Risk Factors in the European Prospective Investigation of Cancer-Norfolk Study. Ophthalmology. 2019;126(1):96–106. pmid:30075201
  14. 14. Seidelmann SB, Claggett B, Bravo PE, Gupta A, Farhad H, Klein BE, et al. Retinal Vessel Calibers in Predicting Long-Term Cardiovascular Outcomes: The Atherosclerosis Risk in Communities Study. Circulation. 2016;134(18):1328–38. pmid:27682886
  15. 15. Zekavat SM, Raghu VK, Trinder M, Ye Y, Koyama S, Honigberg MC, et al. Deep Learning of the Retina Enables Phenome- and Genome-wide Analyses of the Microvasculature. Circulation. 2021. pmid:34743558
  16. 16. Stewart MW. Human immunodeficiency virus and its effects on the visual system. Infect Dis Rep. 2012;4(1):e25. pmid:24470932
  17. 17. Friis-Moller N, Weber R, Reiss P, Thiebaut R, Kirk O, d’Arminio Monforte A, et al. Cardiovascular disease risk factors in HIV patients—association with anti-retroviral therapy. Results from the DAD study. AIDS. 2003;17(8):1179–93. pmid:12819520
  18. 18. Consultation WHOE. Appropriate body-mass index for Asian populations and its implications for policy and intervention strategies. Lancet. 2004;363(9403):157–63. pmid:14726171
  19. 19. Wilson PW D’Agostino RB, Levy D, Belanger AM, Silbershatz H, Kannel WB. Prediction of coronary heart disease using risk factor categories. Circulation. 1998;97(18):1837–47. pmid:9603539
  20. 20. Hippisley-Cox J, Coupland C, Brindle P. Development and validation of QRISK3 risk prediction algorithms to estimate future risk of cardiovascular disease: prospective cohort study. BMJ. 2017;357:j2099. pmid:28536104
  21. 21. Goff DC jr , Lloyd-Jones DM, Bennett G, Coady S, D’Agostino RB Sr., Gibbons R, et al, American College of Cardiology/American Heart Association Task Force on Practice G. 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. J Am Coll Cardiol. 2014;63(25 Pt B):2935–59.
  22. 22. Friis-Moller N, Ryom L, Smith C, Weber R, Reiss P, Dabis F, et al. An updated prediction model of the global risk of cardiovascular disease in HIV-positive persons: The Data-collection on Adverse Effects of Anti-HIV Drugs (D:A:D) study. Eur J Prev Cardiol. 2016;23(2):214–23. pmid:25882821
  23. 23. Agatston AS, Janowitz WR, Hildner FJ, Zusmer NR, Viamonte M Jr., Detrano R. Quantification of coronary artery calcium using ultrafast computed tomography. J Am Coll Cardiol. 1990;15(4):827–32.
  24. 24. Leipsic J, Abbara S, Achenbach S, Cury R, Earls JP, Mancini GJ, et al. SCCT guidelines for the interpretation and reporting of coronary CT angiography: a report of the Society of Cardiovascular Computed Tomography Guidelines Committee. J Cardiovasc Comput Tomogr. 2014;8(5):342–58. pmid:25301040
  25. 25. Zhuo Y, Yu H, Yang Z, Zee B, Lee J, Kuang L. Prediction Factors of Recurrent Stroke among Chinese Adults Using Retinal Vasculature Characteristics. J Stroke Cerebrovasc Dis. 2017;26(4):679–85. pmid:28233623
  26. 26. DeLong ER, DeLong DM, DL C-P. Comparing the areas under two or more correlated receiver operating characteristic curves: a nonparametric approach. Biometrics 1988; 44:837–45. pmid:3203132
  27. 27. He, Kaiming, Xiangyu Zhang, Shaoqing Ren, and Jian Sun. "Deep residual learning for image recognition." In Proceedings of the IEEE conference on computer vision and pattern recognition, pp. 770–778. 2016. https://arxiv.org/pdf/1512.03385v1.pdf
  28. 28. Zee B, Lee J, Li Q, “Method and device for retinal image analysis”, Patent No. US8787638 B2, granted on 22 July 2014. https://patents.google.com/patent/US20120257164A1/en
  29. 29. Fan R.-E., Chen P.-H., and Lin C.-J. “Working set selection using second order information for training support vector machines.” Journal of Machine Learning Research, Vol. 6, 2005, pp. 1889–1918. https://dl.acm.org/doi/pdf/10.5555/1046920.1194907.
  30. 30. Hastie T., Tibshirani R., and Friedman J. The Elements of Statistical Learning, Second Edition. NY: Springer, 2008. https://web.stanford.edu/~hastie/Papers/ESLII.pdf
  31. 31. Scholkopf B., and Smola A. Learning with Kernels: Support Vector Machines, Regularization, Optimization and Beyond, Adaptive Computation and Machine Learning. Cambridge, MA: The MIT Press, 2002. https://dl.acm.org/doi/book/10.5555/559923
  32. 32. Bijker R, Jiamsakul A, Uy E, Kumarasamy N, Ditango R, Chaiwarith R, et al. Cardiovascular disease-related mortality and factors associated with cardiovascular events in the TREAT Asia HIV Observational Database (TAHOD). HIV Med. 2019;20(3):183–91. pmid:30620108
  33. 33. Triant VA, Perez J, Regan S, Massaro JM, Meigs JB, Grinspoon SK, et al. Cardiovascular Risk Prediction Functions Underestimate Risk in HIV Infection. Circulation. 2018;137(21):2203–14. pmid:29444987
  34. 34. Krikke M, Hoogeveen RC, Hoepelman AI, Visseren FL, Arends JE. Cardiovascular risk prediction in HIV-infected patients: comparing the Framingham, atherosclerotic cardiovascular disease risk score (ASCVD), Systematic Coronary Risk Evaluation for the Netherlands (SCORE-NL) and Data Collection on Adverse Events of Anti-HIV Drugs (D:A:D) risk prediction models. HIV Med. 2016;17(4):289–97. pmid:26268806
  35. 35. Edwards-Jackson N, Kerr S, Tieu H, Ananworanich J, Hammer S, Ruxrungtham K, et al. Cardiovascular risk assessment in persons with HIV infection in the developing world: comparing three risk equations in a cohort of HIV-infected Thais. ?2011;12(8):510–5.
  36. 36. Wu PY, Chen MY, Sheng WH, Hsieh SM, Chuang YC, Cheng A, et al. Estimated risk of cardiovascular disease among the HIV-positive patients aged 40 years or older in Taiwan. J Microbiol Immunol Infect. 2019;52(4):549–55. pmid:31036484
  37. 37. Marbaniang IP, Kadam D, Suman R, Gupte N, Salvi S, Patil S, et al. Cardiovascular risk in an HIV-infected population in India. Heart Asia. 2017;9(2):e010893. pmid:29467833
  38. 38. Van Craenendonck T, Gerrits N, Buelens B, Petropoulos IN, Shuaib A, Standaert A, et al. Retinal microvascular complexity comparing mono- and multifractal dimensions in relation to cardiometabolic risk factors in a Middle Eastern population. Acta ophthalmologica. 2021;99(3):e368–e77. pmid:32940010
  39. 39. Gerrits N, Elen B, Craenendonck TV, Triantafyllidou D, Petropoulos IN, Malik RA, et al. Age and sex affect deep learning prediction of cardiometabolic risk factors from retinal images. Sci Rep. 2020;10(1):9432. pmid:32523046
  40. 40. Rim TH, Lee CJ, Tham YC, Cheung N, Yu M, Lee G, et al. Deep-learning-based cardiovascular risk stratification using coronary artery calcium scores predicted from retinal photographs. The Lancet Digital health. 2021;3(5):e306–e16. pmid:33890578
  41. 41. Liew G, Mitchell P, Rochtchina E, Wong TY, Hsu W, Lee ML, et al. Fractal analysis of retinal microvasculature and coronary heart disease mortality. Eur Heart J. 2011;32(4):422–9. pmid:21138936
  42. 42. Hannappe MA, Arnould L, Méloux A, Mouhat B, Bichat F, Zeller M, et al. Vascular density with optical coherence tomography angiography and systemic biomarkers in low and high cardiovascular risk patients. Sci Rep. 2020;10(1):16718. pmid:33028913
  43. 43. Zhong P, Li Z, Lin Y, Peng Q, Huang M, Jiang L, et al. Retinal microvasculature impairments in patients with coronary artery disease: An optical coherence tomography angiography study. Acta ophthalmologica. 2021. pmid:33629471
  44. 44. Chang J, Ko A, Park SM, Choi S, Kim K, Kim SM, et al. Association of Cardiovascular Mortality and Deep Learning-Funduscopic Atherosclerosis Score derived from Retinal Fundus Images. American journal of ophthalmology. 2020;217:121–30. pmid:32222370
  45. 45. Li X, Wong WL, Cheung CY, Cheng CY, Ikram MK, Li J, et al. Racial differences in retinal vessel geometric characteristics: a multi-ethnic study in healthy Asians. Invest Ophthalmol Vis Sci. 2013;54(5):3650–6. pmid:23652487
  46. 46. Monteiro-Henriques I, Rocha-Sousa A, Barbosa-Breda J. Optical coherence tomography angiography changes in cardiovascular systemic diseases and risk factors: A Review. Acta ophthalmologica. 2021. pmid:33783129