Topographic analysis of macular choriocapillaris flow deficits in diabetic retinopathy using swept–source optical coherence tomography angiography
International Journal of Retina and Vitreous volume 6, Article number: 6 (2020)
The purpose of this study was to investigate the association between diabetic retinopathy (DR) severity and macular choriocapillaris (CC) flow deficit percentage (FD %) in different macular regions using swept-source optical coherence tomography angiography (SS-OCTA).
Diabetic patients with SS-OCTA images were graded by severity and retrospectively assessed. CC FD % was calculated in four different regions of the OCTA image: inner, middle, outer, and full-field region. The generalized estimating equations (GEE) approach for clustered eye data was used to determine effect size and significance of age and disease severity on FD % for each region.
160 eyes from 90 total diabetic patients met inclusion criteria. Out of 90 patients, 33 had no DR, 17 had mild nonproliferative DR (NPDR), 8 had moderate NPDR, 10 had severe NPDR and 22 had proliferative DR. Age and DR severity had a significant positive association with FD % for each region studied with a greater effect in the two centermost regions. The increase in flow deficit percentage per year of age by region was: inner 0.12 (p < 0.001), middle 0.09 (p < 0.001), outer 0.05 (p < 0.001, full-field 0.06 (p < 0.001). The increase in flow deficit percentage per increase in diabetic retinopathy severity stage by region was: inner 0.65 (p < 0.0087), middle 0.56 (p < 0.0012), outer 0.33 (p < 0.045), full-field 0.36 (p < 0.018).
Topographic analysis of the CC FD % in diabetic eyes suggests that CC flow impairment corresponds to DR severity, with all studied regions of the CC significantly affected. There was greater regional impairment due to age and disease severity in the inner and middle regions.
Diabetic retinopathy (DR) is the most common microvascular complication of diabetes, and its diagnosis and treatment are critical to ophthalmic practice [1, 2]. In addition to retinopathy, diabetic choroidopathy has also been recognized as part of the vascular inflammatory effect of diabetes [3, 4]. Prior studies show that diabetic choroidopathy can be assessed clinically with optical coherence tomography (OCT), demonstrating decreased choroidal thickness in diabetic patients as well as decreased capillary perfusion density in the choriocapillaris (CC) [5,6,7,8]. The ability to analyze the choroid, and specifically the choriocapillaris using OCT angiography (OCTA), has led to new descriptions of the pathophysiologic effects of diabetes [9, 10].
The CC is a single layer of capillaries that supplies the retinal pigment epithelium and photoreceptors . It is a dense network that is unique among capillary beds in its increased caliber and permeability . The recent increased interest in the choriocapillaris results directly from advances in OCTA, a technology that uses OCT motion signal to image depth-resolved blood flow. Previous studies of the CC were limited to post-mortem histology as other methods of imaging the vasculature, such as fluorescein or ICG angiography, cannot easily visualize the small, fenestrated vessels of the CC . In addition, OCT imaging of the CC has historically been challenging because of the signal scattering effect of the RPE . New technologies in OCT such as swept-source laser sources have increased the quality of images as swept-source does not suffer from the sensitivity roll-off of spectral-domain technology and penetrates into deeper layers using infrared wavelengths. This has allowed better imaging of the CC, approaching the resolution of histologic methods . With swept-source OCTA (SS-OCTA) the CC structure and functional perfusion may now be studied in vivo .
OCTA imaging of the CC has already allowed for a better understanding of CC changes associated with natural as well as pathologic processes. For instance, the percent of CC area without flow (called flow deficit percentage) in the choriocapillaris has been shown to be significantly affected by age [16, 17]. Specifically, CC flow voids increase with age and the association rises with proximity to the central fovea, with more deficits in the central region [18, 19]. Qualitatively, it has also been shown that CC flow impairment is present in early stages of diabetes sometime even before the appearance of DR on clinical exam, and that long-standing diabetes is correlated with CC drop-out on histology [13, 20]. While there is evidence that CC loss exists in eyes with DR, it is not known whether this degeneration follows a pattern or affects all areas equally. The aim of this study is to examine the relationship between the severity of DR and CC loss. Specifically, it asks how DR severity affects different regions of the CC, given that confounders such as age also affect these regions.
This study included diabetic subjects who were imaged with SS-OCTA in the retina clinic at the New England Eye Center, Boston, MA between January 2017 and July 2018. The study was performed in accordance with the Health Insurance Portability and Accountability Act and adhered to the principles of the Declaration of Helsinki. The subjects signed informed consents and the research was approved by the Institutional Review Board for Tufts Medical Center.
Exclusion criteria comprised patients with a concurrent diagnosis of AMD or other retinal or choroidal vascular disorders, laser in the field of the image, or poor image quality/segmentation. Image quality exclusion was qualitative, rather than based on a signal strength cutoff. Images were graded by an expert grader (author A. Yasin Alibhai MD) into the following groups using fundus imaging and fluorescein angiography when available: diabetes with no diabetic retinopathy (DMnoDR), mild non-proliferative diabetic retinopathy (NPDR), moderate NPDR, severe NPDR, and proliferative diabetic retinopathy (PDR) according to Early Treatment of Diabetic Retinopathy Study (ETDRS) methodology . History of focal laser, pan-retinal photocoagulation (PRP), or anti-vascular endothelial growth factor (anti-VEGF) injection was recorded. Hemoglobin A1c (HbA1c) from within 1 year of imaging was collected.
Images were acquired using the PLEX Elite 9000 SS-OCTA (Carl Zeiss Meditec Inc., Dublin, CA, USA). The PLEX Elite 9000 has an A-scan rate of 100,000 scans per second with a central wavelength of 1060 nm, axial resolution of ~ 6 μm, lateral resolution of ~ 14 μm, and scans at a depth of 3 mm. The standard 6 mm x 6 mm protocol was used, centered on the fovea. This scan protocol consists of 500 B-scans with 500 A-scans per B-scan. The automated segmentation algorithm on the PLEX was used to segment the choriocapillaris slab from 29 µm to 49 µm posterior to the automatically-fitted retinal pigment epithelium line, as is standard on the device [15, 22]. The resulting angiography en face image was exported for analysis. Some scans were further excluded due to poor imaging (for example, motion artifact, poor segmentation of the choriocapillaris by the automated algorithm, or laser scar in the field of the image).
Image processing and analysis
All image processing was performed using a custom Matlab program (MathWorks, Natick, MA, USA). Flow deficit images were computed by binarizing the en face OCTA images using a fixed global threshold (threshold value of 0.2, with OCTA data values taking on values in the [0, 1] range); this value was subjectively chosen, independent of this study, on the basis of the qualitative appearance of flow deficits in OCTA data of healthy eyes. Figure 1 depicts representative cases for each DR stage of the binarized images. To reduce the influence of noise, all flow deficits having an area of less than 500 µm2 were removed from the analysis. Flow deficit percentage for a given region is calculated as the area of non-flow in the region divided by the total area of the region multiplied by 100. Based on the ETDRS grid, flow deficit percentage was calculated for the Inner (0.5 mm radius circle), Middle (0.5 mm–1.5 mm annulus), Outer (1.5 mm–3 mm annulus) and Full-Field regions for each 6 mm × 6 mm binarized image (Fig. 2). These methods described are similar to those used in the literature [18, 23].
Statistical analysis was done using the open source statistical program R (R Foundation for Statistical Computing, Vienna, Austria) . The specific packages and scripts implemented were ggplot2 for graphing, geepack for GEE analysis, and MESS for the quasi-Akaike Information Criterion (QIC) function, which allows for model confirmation [25,26,27]. Mean FD % in each region by age group and DR severity were calculated. Mean HbA1c was calculated for each severity group. This study included paired eye data from subjects, which introduced dependency in the outcome measured that could be resolved by analyzing the data with a generalized estimating equations (GEE) approach. The GEE analysis is an extension of the generalized linear model that produces a model of the effects, and their significance levels, on a dependent variable (in this case FD %); specifically it is able to take into account correlated data such as paired-eye data. The GEE analysis determined the coefficient and significance of the contribution to the model for the different variables of interest, age and severity . An exchangeable correlation structure was chosen using both a priori assumptions and confirmed using QIC testing [27,28,29]. Each patient served as a cluster with correlated outcome measures (left eye, right eye). Disease stage was coded as on a Likert scale and was treated as a continuous variable. A p-value of .05 was used as a significance cut-off for all applicable measures.
Out of 186 eyes, 26 eyes were excluded due to poor image quality or coexisting unrelated retinal disease. 160 eyes from 90 diabetic patients with SS-OCTA imaging met inclusion criteria. The mean age was 58.14 years. 51 were male and 39 were female. HbA1c values were available for 71 of the 90 patients. Study population characteristics are outlined in Table 1. Of the total 160 eyes, 19 eyes had a history of PRP, 21 eyes had a history of focal laser with scarring outside of the area of scanning, and 37 eyes had a history of anti-VEGF injection.
Using simple averages, uncorrected for other variables, age and DR Severity showed increasing FD % with proximity to the fovea (Figs. 3 and 4). Using the GEE approach, which takes account for paired eye correlation and other included variables, age and DR severity were independently shown to have a significant positive association with FD % with greater effect size in the central two regions (Table 2). The coefficient can be thought of as the amount FD % would increase given a one unit increase in the covariate (age or DR Severity). Other GEE models tested included additional covariates such as gender, history of VEGF injection, history of PRP, history of focal laser, presence of edema, and eye laterality. None of these additional covariates showed significance, in agreement with a previous study . Thus only results for the model including age and DR severity are reported. The α correlation values are non-zero, showing correlation between eyes of the same subject, which justifies using the GEE approach.
Improvements in OCTA technology have only recently allowed for detailed, in vivo study of the CC. This study affirms previous findings of the effects of age and DR on the CC, and offers new insights into the effects of DR severity on flow deficit percentage in the CC. An important tool used in this study was the generalized estimating equation which allows for the use of fellow eyes, increasing the power of the study as well as the ability to identify the significance and effect size of multiple covariates .
This study agrees with previous findings in both OCTA and histology that suggest the CC degenerates with age [16,17,18,19]. Our results suggest that aging is significantly associated with an increase in flow deficit percentage of the CC, and this finding increases with proximity to the central fovea. The effect of age on FD % in the innermost regions studied is roughly twice that in the outer region and full-field measures. Given that multiple studies have found this effect, it becomes clear that future studies must control for these factors to correctly interpret findings.
While diabetes is a systemic disease, traditionally the ophthalmologic focus has been on its retinal effects. Diabetes as a choroidopathy is relatively under-studied partly because, until recently, detecting changes in the choroidal vasculature was difficult (limited to angiography and histology). OCTA and enhanced depth imaging (EDI) on OCT have extended the ophthalmologist’s view into the choroid . Using OCTA, Choi et al. and Forte et al. reported regions of CC impairment in subjects with diabetic retinopathy [20, 32]. These authors demonstrated that most subjects, even with non-clinical DR, have CC alterations. Dodo et al. found that non-perfusion of the CC is correlated with diabetic retinopathy and visual acuity. The group concluded that choroidal nonperfusion is likely to part of the pathogenesis of decreased vision and is associated with a disrupted photoreceptor layer .
This study found that DR severity had a significant positive association with FD %. That is, as diabetic retinopathy worsened there was a corresponding increasing deficit in the CC. As with age, this effect also increased with proximity to the fovea. Thus, for a progression of one stage in severity the flow deficit percentage increased (worsened) for the two inner regions more than the outer and full field regions. As an example, the inner region coefficient size was 1.8 times that of the effect for the full field image (0.65 versus 0.36, respectively).
This study’s results are qualified by several limitations. Drawbacks of this study include its retrospective nature, lack of age-matched controls, and the relative low numbers of patients in the moderate and severe NPDR groups. Additionally, this study’s results are valid for the effect of DR severity but do not provide a granular pairwise analysis between stages. This means that a progression from no DR to mild NPDR is treated the same as a progression from severe NPDR to PDR. The GEE analysis is an extension of the generalized linear model (GLM), but the effect of disease stage may be such that a non-linear model would better capture the relationship between DR severity and FD %. However, the strengths of the GEE allowed for accounting for correlation between paired eyes, as previously discussed. Moreover, a linear model still allows for conclusions to be made about the direction and magnitude of effects. Future studies should address these concerns as well as look prospectively at the effect of treatment (both pharmaceutical and lifestyle modifications) on FD %. Disease duration should also be accounted for as this has been reported to be an even stronger predictor of disease progression and would likely be a significant covariate of FD % as well .
Swept-source optical coherence tomography angiography imaging of diabetic eyes reveals increased flow deficits in the choriocapillaris with increasing proximity to the central macula that correlate with more severe disease. This study also reaffirms previous findings regarding the topographic effects of age on flow deficit in the choriocapillaris.
Availability of data and materials
The datasets used and/or analyzed during the current study are available from the corresponding author on reasonable request.
Anti-vascular endothelial growth factor
Diabetes with no diabetic retinopathy
- FD %:
Flow deficit percentage
Early Treatment of Diabetic Retinopathy Study
Enhanced depth imaging
Non-proliferative diabetic retinopathy
Optical coherence tomography
Optical coherence tomography angiography
Swept-source optical coherence tomography angiography
Generalized estimating equations
Antonetti DA, Klein R, Gardner TW. Diabetic retinopathy. N Engl J Med. 2012;366:1227–39.
Duh EJ, Sun JK, Stitt AW. Diabetic retinopathy: current understanding, mechanisms, and treatment strategies. JCI Insight. 2017. https://doi.org/10.1172/jci.insight.93751.
Hidayat AA, Fine BS. Diabetic choroidopathy. Light and electron microscopic observations of seven cases. Ophthalmology. 1985;92:512–22.
Lutty GA. Diabetic choroidopathy. Vis Res. 2017;139:161–7.
Querques G, Lattanzio R, Querques L, Del Turco C, Forte R, Pierro L, et al. Enhanced depth imaging optical coherence tomography in type 2 diabetes. Invest Ophthalmol Vis Sci. 2012;53:6017–24.
Adhi M, Brewer E, Waheed NK, Duker JS. Analysis of morphological features and vascular layers of choroid in diabetic retinopathy using spectral-domain optical coherence tomography. JAMA Ophthalmol. 2013;131:1267–74.
Gildea D. The diagnostic value of optical coherence tomography angiography in diabetic retinopathy: a systematic review. Int Ophthalmol. 2018. https://doi.org/10.1007/s10792-018-1034-8.
Nagaoka T, Kitaya N, Sugawara R, Yokota H, Mori F, Hikichi T, et al. Alteration of choroidal circulation in the foveal region in patients with type 2 diabetes. Br J Ophthalmol. 2004;88:1060–3.
Ferrara D, Waheed NK, Duker JS. Investigating the choriocapillaris and choroidal vasculature with new optical coherence tomography technologies. Prog Retin Eye Res. 2016;52:130–55.
Wang JC, Laíns I, Providência J, Armstrong GW, Santos AR, Gil P, et al. Diabetic choroidopathy: choroidal vascular density and volume in diabetic retinopathy with swept–source optical coherence tomography. Am J Ophthalmol. 2017;184:75–83.
Tan B, Chua J, Barathi VA, Baskaran M, Chan A, Lin E, et al. Quantitative analysis of choriocapillaris in non-human primates using swept-source optical coherence tomography angiography (SS-OCTA). Biomed Opt Expr. 2019;10:356–71.
Agarwal A. Gass’ atlas of macular diseases e-book. Edinburgh: Elsevier health sciences; 2011.
Fryczkowski AW, Sato SE, Hodes BL. Changes in the diabetic choroidal vasculature: scanning electron microscopy findings. Ann Ophthalmol. 1988;20:299–305.
Choi W, Mohler KJ, Potsaid B, Lu CD, Liu JJ, Jayaraman V, et al. Choriocapillaris and choroidal microvasculature imaging with ultrahigh speed OCT angiography. PLoS ONE. 2013;8:e81499.
Lauermann JL, Eter N, Alten F. Optical coherence tomography angiography offers new insights into choriocapillaris perfusion. Ophthalmologica. 2018;239:74–84.
Spaide RF. Choriocapillaris flow features follow a power law distribution: implications for characterization and mechanisms of disease progression. Am J Ophthalmol. 2016;170:58–67.
Ramrattan RS, van der Schaft TL, Mooy CM, de Bruijn WC, Mulder PG, de Jong PT. Morphometric analysis of Bruch’s membrane, the choriocapillaris, and the choroid in aging. Invest Ophthalmol Vis Sci. 1994;35:2857–64.
Zheng F, Zhang Q, Shi Y, Russell JF, Motulsky EH, Banta JT, et al. Age–dependent changes in the macular choriocapillaris of normal eyes imaged with swept–source oct angiography. Am J Ophthalmol. 2019. https://doi.org/10.1016/j.ajo.2018.12.025.
Nassisi M, Baghdasaryan E, Tepelus T, Asanad S, Borrelli E, Sadda SR. Topographic distribution of choriocapillaris flow deficits in healthy eyes. PLoS ONE. 2018;13:e0207638.
Forte R, Haulani H, Jürgens I. Quantitative and qualitative analysis of the three capillary plexuses and choriocapillaris in patients with type 1 and type 2 diabetes mellitus without clinical signs of diabetic retinopathy: a prospective pilot study. Retina. 2018. https://doi.org/10.1097/IAE.0000000000002376.
Early Treatment Diabetic Retinopathy Study Research Group. Grading diabetic retinopathy from stereoscopic color fundus photographs—an extension of the modified airlie house classification. ETDRS report number 10. Ophthalmology. 1991;98:786–806.
Sugano Y, Sekiryu T, Furuta M, Tomita R, Shintake H, Maehara H, et al. Morphometrical evaluation of the choriocapillaris imaged by swept-source optical coherence tomography angiography. Clin Ophthalmol. 2018;12:2267–76.
Moult EM, et al. Spatial distribution of choriocapillaris impairment in eyes with choroidal neovascularization secondary to age–elated macular degeneration: a quantitative oct angiography study. Retina. 2019. https://doi.org/10.1097/iae.0000000000002556.
R Core Team. R: A Language and Environment for Statistical Computing [Internet]. Vienna, Austria: R Foundation for Statistical Computing; 2013. http://www.R-project.org/.
Højsgaard S, Halekoh U, Yan J. The R package geepack for generalized estimating equations. J Stat Softw. 2005;15:1–11.
Wickham H. ggplot2: elegant graphics for data analysis. 1st ed. New York: Springer; 2009.
Ekstrm C. MESS package | R Documentation [Internet]. [cited 2019 Jun 5]. https://www.rdocumentation.org/packages/MESS/versions/0.5.5.
Ying G-S, Maguire MG, Glynn R, Rosner B. Tutorial on biostatistics: linear regression analysis of continuous correlated eye data. Ophthalmic Epidemiol. 2017;24:130–40.
Pan W. Akaike’s information criterion in generalized estimating equations. Biometrics. 2001;57:120–5.
Liang K-Y, Zeger SL. Longitudinal data analysis using generalized linear models. Biometrika Narnia. 1986;73:13–22.
Hua R, Liu L, Wang X, Chen L. Imaging evidence of diabetic choroidopathy in vivo: angiographic pathoanatomy and choroidal-enhanced depth imaging. PLoS ONE. 2013;8:e83494.
Choi W, Waheed NK, Moult EM, Adhi M, Lee B, De Carlo T, et al. Ultrahigh speed swept source optical coherence tomography angiography of retinal and choriocapillaris alterations in diabetic patients with and without retinopathy. Retina. 2017;37:11–21.
Dodo Y, Suzuma K, Ishihara K, Yoshitake S, Fujimoto M, Yoshitake T, et al. Clinical relevance of reduced decorrelation signals in the diabetic inner choroid on optical coherence tomography angiography. Sci Rep. 2017;7:5227.
Klein R, Knudtson MD, Lee KE, Gangnon R, Klein BEK. The Wisconsin epidemiologic study of diabetic retinopathy: XXII the twenty-five-year progression of retinopathy in persons with type 1 diabetes. Ophthalmology. 2008;115:1859–68.
A preliminary version of these results was presented at the International Retinal Imaging Symposium held at Doheny Eye Institute UCLA, Los Angeles CA. March 16, 2019.
This work was supported by the Macula Vision Research Foundation (West Conshohocken, PA, USA), the Massachusetts Lions Clubs (Belmont, MA, USA), the National Institutes of Health (Grant Number 5-R01-EY011289-31), the Air Force Office of Scientific Research (Grant Number FA9550-15-1-0473), the Champalimaud Vision Award (Lisbon, Portugal), the Beckman-Argyros Award in Vision Research (Irvine, CA, USA), and the Yale School of Medicine Medical Student Fellowship (New Haven, CT, USA).
Ethics approval and consent to participate
This study was approved by the Institutional Review Board for Tufts Medical Center and all patients provided informed consent.
Consent for publication
AI: Topcon Medical Systems, Inc. (Tokyo, Japan) (S), Nidek Medical Products, Inc. (Gamagori, Aichi, Japan) (S) CRB: Carl Zeiss Meditec, Inc. (Dublin, CA, USA) (S), Genentech (South San Francisco, CA, USA), Novartis (Basel, Switzerland) (S), Allergan (Dublin, Ireland) (C) JGF: Topcon Medical Systems Inc. (Oakland, NJ, USA) (C, F), Optovue, Inc. (Fremont, CA, USA) (C, R, S), Carl Zeiss Meditec, Inc. (Dublin, CA, USA) (R) JSD: Carl Zeiss Meditec, Inc. (Dublin, CA, USA) (C, F), Optovue, Inc. (Fremont, CA, USA) (C, F) NKW: Macula Vision Research Foundation (West Conshohocken, PA, USA) (F), Topcon Medical Systems, Inc. (Oakland, NJ, USA) (F), Nidek Medical Products, Inc. (Fremont, CA, USA) (S, F), Optovue, Inc. (Fremont, CA, USA) (C), Carl Zeiss Meditec, Inc. (Dublin, CA, USA) (F) ER: Regeneron Pharmaceuticals, Inc. (Tarrytown, NY, USA) (F).
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Gendelman, I., Alibhai, A.Y., Moult, E.M. et al. Topographic analysis of macular choriocapillaris flow deficits in diabetic retinopathy using swept–source optical coherence tomography angiography. Int J Retin Vitr 6, 6 (2020). https://doi.org/10.1186/s40942-020-00209-0