Factors based on optical coherence tomography correlated with vision impairment in diabetic patients

The aim of this study is to evaluate the relationship between retinal structures and visual acuity in diabetic patients using optical coherence tomography (OCT) and OCT angiography (OCTA). This study was a retrospective observational study conducted at a single medical center in Japan. Evaluation of retinal images was analyzed using spectral domain OCT. Twelve factors including central retinal thickness, length of disorganization of retinal inner layer (DRIL), number of inner hyperreflective foci, number of outer hyperreflective foci, height of intraretinal fluid, height of subretinal fluid, length of external limiting membrane disruption, length of external ellipsoid zone (EZ) disruption, vessel density of superficial capillary plexus (SCP), foveal avascular zone (FAZ) area, and FAZ circularity were analyzed based on OCT/OCTA findings. Multivariate analysis was used to investigate the OCT-based factors that could be correlated with poor visual acuity in treatment-naïve diabetic eyes. A total of 183 eyes of 123 diabetic patients with type 2 diabetes (mean age 61.9 ± 12.3 years, 66 men and 57 women) and 62 eyes of 55 control subjects (mean age 64.4 ± 12.5 years, 15 men and 40 women) was enrolled in this study. Multiple regression analysis showed that OCT-based factors correlated with visual acuity were length of DRIL (β = 0.24, P < 0.01), length of EZ disruption (β = 0.35, P < 0.001), and FAZ circularity (β =  − 0.14, P < 0.05). The other factors showed no significant correlation. In conclusion, the length of DRIL, length of EZ disruption, and FAZ circularity measured by OCT were identified as related factors for visual impairment in treatment-naïve diabetic eyes.


Stage of DR NDR (n = 92) NPDR (n = 67) PDR (n = 24) P value
Age (   www.nature.com/scientificreports/ Significance of OCT-based factors affecting visual acuity. DRIL. DRIL was characterized by unidentifiable boundaries of the ganglion cell layer (GCL)-inner plexiform layer (IPL) complex, inner nuclear layer (INL), and outer plexiform layer (OPL), correlated with worse visual acuity in DME, and the changes in DRIL predicted subsequent visual outcome 24,25 . Ocular disorders caused by retinal vascular dysfunction are likely to complicate disruption of the blood-retinal barrier (BRB), followed by extracellular fluid accumulation in the intraretinal or subretinal space 26 . Although the mechanism of DRIL formation is unknown, it is hypothesized that DM-induced microvascular damages inside the retina may represent structural deformation identified as DRIL based on OCT images 27 . In addition, the PDR may form a fibrovascular membrane that creates traction on the macula. Therefore, multiple mechanisms associated with vascular abnormality and mechanical stress may be involved in DRIL formation. This can lead to deformation or breakage of the synaptic junction between photoreceptors and retinal ganglion cells. Furthermore, although a histological observation of DRIL has not been demonstrated, DRIL reflects damage of various cells such as Müller cells, bipolar cells, horizontal cells, and amacrine cells, and is believed to affect the visual outcome of diabetic eyes 25 . Our results, showing the correlation with the length of DRIL and visual acuity, support the findings of previous reports, and the widespread presence of DRIL may contribute to visual dysfunction in treatment-naïve diabetic eyes.
EZ disruption. It is indisputable that the retinal photoreceptor layer can be accurately evaluated by examining the integrity of the EZ using OCT, which corresponds to the inner segment of photoreceptors, indicating that EZ integrity has been shown to be strongly correlated with visual acuity 28 . Although the mechanism of photoreceptor degeneration is not yet well defined, several processes may lead to photoreceptor damage. Destructed BRB associated with DM may induce extravasation of blood components and inflammatory cell infiltration, causing morphological abnormalities in the macula 29 . Uji et al. reported that HF is considered a lipoprotein, a precursor of hard exudates, which was incorporated by photoreceptors or macrophages, leading to poor visual acuity 30 . Murakami et al. suggested that incomplete EZ often exists under the cyst space, and that the cyst space of the outer plexiform layer may contribute to photoreceptor damage in DME 31 . In our study, EZ disruption was found in 10.9% of all DM patients, 10.4% in NPDR and 54.2% in PDR, which occurred more frequently in the advanced stages of the disease. Our study also confirmed that the patients showing high frequency of EZ disruption had significantly poorer visual acuity.
FAZ circularity. DR is primarily a disease of the retinal vasculature, and OCTA provides valuable information about vascular conditions not available with structural OCT. FAZ is a specialized area of the retina containing the highest density of cone photoreceptors and high oxygen consumption 32 . As the FAZ of the healthy eye usually exhibits a circular or elliptical shape, deviations from this shape are common in vascular pathology 33,34 . Balaratnasingam and colleagues used OCTA to reveal that FAZ area correlates with visual acuity in patients with DR, and proposed FAZ as a biomarker of visual function 27 . FAZ circularity is defined as the ratio of the FAZ boundary to the perimeter of a circle of equal area 33 . Krawitz and colleagues found an association between the disrupted circularity and DR severity, but reported no correlation with visual acuity 35 . On the other hand, Tang et al. demonstrated that FAZ area was not a sensitive marker that correlates with central vision, while decreased FAZ circularity is associated with impaired visual function 36 . In our study, we showed a correlation between visual acuity and FAZ circularity, in which eyes with poor visual acuity had a more irregular shape of FAZ. Alipour et al. explained that studies using the digital curvelet transform revealed that the FAZ of DR deviated from the mildly wavy boundaries found in healthy controls, which became apparent in severely diseased eyes 37 . Defects in the FAZ margin are mainly due to capillary loss and vascular remodeling. The FAZ circularity index allows quantification of disruption of the terminal capillaries in the fovea, which may be a more relevant measure of visual acuity. In addition, macular edema can result in decreased vascular elasticity due to mechanical stretching that results in blockage of blood vessels 38 , and can also change FAZ contours 39 . Therefore, assessing both the size and shape of FAZ is probably important for detecting pathological alteration of the macula in DR.

Changes in outer and inner retina.
In summary, our study showed that the following 3 OCT-based factors could play important roles in impaired visual function in diabetic eyes: length of DRIL, EZ disruption, and FAZ circularity. In addition, these factors had a significant correlation with each other in treatment-naïve diabetic eyes. However, there is a lack of information on how these factors correlate with each other. This study further discusses three categories: (1) DRIL and FAZ circularity, (2) DRIL and EZ disruption, and (3) EZ disruption and FAZ circularity.  41 . The metabolic and oxidative stresses of diabetes lead to increased sensitivity of ganglion cells, leading to neuronal loss 42 . After that, DRIL occurs with the severity of retinopathy, which may be negatively correlated with FAZ circularity in SPC. In fact, the frequency of DRIL and FAZ circularity found in this study was significantly high as DR severity increased. Taken together, the relationship between DRIL and FAZ circularity could be associated with both neurovascular damages during DR progression, although segmentation errors in OCTA might exhibit.
DRIL and EZ disruption. The sub-analysis of this study showed correlation with DRIL and EZ disruption, suggesting that the extent of DRIL may also be involved in the retinal outer layer in DR. It remains unknown how different cell types die during DR progression. DR is associated with loss of pericytes and endothelial cells in the retinal vasculature, leading to loss of the capillary and BRB 43 . The deep capillary plexus (DCP) including the INL and the OPL are probably important for nutrition to photoreceptor synapses, as the retinal circulation is responsible for approximately 15% of the oxygen supply to photoreceptor inner segments in the darkness 44 . Therefore, photoreceptor cells in the presence of DRIL may be susceptible to ischemic insults in DCP. The presence of DRIL indicates an anatomical disruption of visual transduction pathways in the medial retina and may interfere with the transmission of visual information from photoreceptors to ganglion cells. The anatomical changes observed in this OCT may cause the association with DRIL and EZ disruption.
EZ disruption and FAZ circularity. It has also been shown that diabetes is associated with increased loss of retinal neurons 45 . Retinal glial cells including astrocytes, Müller cells and microglia are involved in structural support and maintenance of homeostasis in the retina 46 . Under hyperglycemic stress, microglia are activated, various cytokines are secreted, and astrocytes and Müller cells are involved with the amplification of the inflammatory response by the production of proinflammatory cytokines 47 . Previous studies have reported that FAZ measurements of both SCP and DCP correlate with EZ and IZ disruption, and photoreceptor morphology and integrity are susceptible to macular ischemia 48 . Long-term or more severe macular ischemia not only contributes to the onset of macular edema, but may also cause retinal atrophy due to photoreceptor damage. It remains open to debate whether diabetic retinal neuropathy is the cause or effect of microangiopathy under chronic hyperglycemia. Barber proposed two hypotheses to explain the relationship between neurodegenerative and microangiopathy changes 49 . First, loss of BRB integrity manifests itself as increased vascular permeability, leading to edema and neuronal loss due to uncontrolled extracellular fluid composition in the retina. Alternatively, DM directly affects metabolism in the neural retina, increasing apoptosis of cells in the neural retina causing a breakdown of the BRB. Despite the close relationship between neurodegeneration and microangiopathy, the progression of DR is ultimately irreversible when untreated. To protect the eyesight of diabetic patients, it is necessary to find ways to prevent the gradual loss of neurons in the retina. Further studies are required to prove the validated association between OCT-based retinal structural abnormalities and visual function in diabetic eyes.

Limitations
We are aware of some limitations of this study. First, this report includes a retrospective design and a limited number of subjects. A larger sample size in prospective studies is important to reconfirm the results. Second, it was impossible to investigate the VD of the DCP, as the OCTA parameters automatically calculated in AngioPlex relate only to SCP. However, this technique produces projection artifacts in most cases, so SCP analysis is more accurate than DCP one 50 . Third, if the capillary flow rate detected by the OCTA algorithm is below a threshold, the terminal capillaries may not be fully visible. Therefore, measurements of FAZ area and VD may be affected in selected cases. Fourth, demographics such as age, duration of diabetes, HbA1c levels, and hypertension were not used in the regression analysis. This needs to ensure a larger sample size and will reveal a correlation with vision loss. Fifth, the biomorphometry performed in this study utilized data from a single visit. Therefore, it is impossible to use this result to judge the risk of disease progression or to predict response to treatment for DR. Further studies are needed on the natural history of DR and the correlation between post-treatment OCT-based factors and visual acuity. Sixth, we could not examine the degree of pathological alteration and visual impairment in the outer and inner retinas confirmed by OCT images. We need to secure more samples in the future to investigate how parafoveal signaling is related to the degree of visual impairment. Seventh, it was impossible to examine OCT-based variables that correlate with visual acuity in each severity grade of DR, since the subjects were limited to treatment-naïve eyes. Further case collection is required to perform multivariate analysis in each group.

Conclusions
The length of DRIL, length of EZ disruption, and FAZ circularity measured by OCT were identified as candidate factors for visual impairment in treatment-naïve diabetic eyes.

Materials and methods
Patients. This study was a retrospective observational study of patients who visited Teine Keijinkai Hospital from May 2019 to July 2020. The research plan of this study in compliance with the principles of the Helsinki Declaration was approved by Medical Ethics Committee in Teine Keijinkai Hospital (Registration number: 2-020041-00). Written informed consent from each patient was waived because this is a retrospective observational study. Type 2 DM patients and non-DM controls in the Japanese population were included in this study. All type 2 DM patients have been confirmed by a diabetic specialist in Teine Keijinkai Hospital General Internal Medicine and Infectious Diseases Internal Medicine Department or a nearby physician. All subjects underwent a comprehensive ophthalmic examination including best-corrected visual acuity converted to logarithm of the minimum angle of resolution (logMAR), Goldmann applanation tonometry, biomicroscopy of the anterior segment using a slit lamp examination, ophthalmoscopy of the posterior segment with dilated pupils, and axial length measured using optical biometry (IOL Master Zeiss; Jena, Germany) and SD-OCT (Cirrus HD-OCT, Carl Zeiss Meditec, Dublin, CA). Exclusion criteria were eyes with spherical power greater than -5 diopters, cylindrical astigmatism with cylinder power greater than 3 diopters, ocular axial length being greater than 26 mm, cataracts or vitreous hemorrhages, other fundus diseases, and a history of other fundus diseases. Moreover, eyes with a history of previous topical treatments including retinal laser photocoagulation, local vascular endothelial growth factor (VEGF), and/or triamcinolone acetonide injection, or ocular surgeries including vitrectomy or cataract surgery were excluded. Funduscopic examination including the periphery of the retina was performed with a pharmacologically dilated pupil, and the severity of retinopathy was classified into three major grades based on the International Clinical Diabetic Retinopathy and Diabetic Macula Edema Disease Severity Scales 4 : no DR; mild/moderate/severe NPDR; and PDR. Systemic characteristics of diabetic patients were recorded: duration of diabetes, HbA1c, HT, SBP and DBP, eGFR, DL, and TC.
OCT and OCTA image acquisition. The AngioPlex OCTA software featured an OCT microangiography algorithm that incorporated differences in both phase and intensity information, and a FastTrac retinal tracking system to reduce motion artifacts 51 . OCT imaging used three Scan acquisition options: Macular Cube Scan 200 × 200, single horizontal HD 1 Line 100 × (100 × averaged), and Angiography Scan 6 × 6. They were obtained by an experienced technician and only scans with a signal strength of more than 7/10 (out of 10) were evaluated. To avoid fixation errors, all scans were focused on the fovea and the segmentation errors were corrected manually. Each measurement was independently evaluated by two masked investigators (H.E. and H.T.). When there was disagreement, the evaluation was discussed to reach a consensus.
OCT-based factors. The presence of several morphological features was evaluated from B-scan and C-scan images ( Figs. 1 and 2); CRT, DRIL, HF, IRF, SRF, ELM, and EZ. Quantitative assessment of CRT was obtained from the central subfield of the Early Treatment Diabetic Retinopathy Study (ETDRS) grid with automated measurements (Fig. 1). Measurement of each lesion was measured manually using a caliper tool mounted on the instrument. DRIL was defined as a horizontal range where the boundaries of the GCL-IPL complex, INL, and OPL could not be identified ( Fig. 2A,D). HF was counted manually in each scan (Fig. 2B,E). To exclude hard exudates and microaneurysms from the analysis, we defined them as having the following morphological features: (1) reflectivity similar to that of nerve fiber layer; (2) absence of back-shadowing; and (3) < 30 μm diam- www.nature.com/scientificreports/ eter. HF was divided into two groups based on their locations in the retinal layers: the inner retinal layer was defined as the layer between the inner limiting membrane (ILM) and the lower boundary of the outer plexiform layer, and the outer retinal layer was defined as the layer between the outer nuclear layer and the retinal pigment epithelium (RPE)/Bruch's membrane boundaries. IRF was defined as the vertical extent of the hypo-reflection space observed in the outer nuclear layer, and when multiple IRF were present, the IRF closest to the center was selected (Fig. 2B,E). SRF was defined as the vertical distance in the largest hypo-reflective space between the sensory retina and the RPE inner edge (Fig. 2B,E). The extent of ELM and EZ destruction was defined as the horizontal extent of the obscured line (Fig. 2C,F). The extent of DRIL lesions, the number of HF, the presence of IRF and SRF, the extent of ELM and EZ disruption were evaluated in a 1500-μm wide region of the B-scan image centered on the fovea (Fig. 2).
OCTA parameters. OCTA images of the SCP were generated using an automated software algorithm ( Fig. 3A,D). The upper boundary of the SCP slab is the ILM and the lower boundary is an approximation of the IPL inferred by the following Eq. (1): where Z, T, and OPL indicate the boundary location, thickness, and outer plexiform layer, respectively. Several quantitative parameters were evaluated from OCTA images automatically generated from Angiography scan 6 × 6 mm slabs according to previous reports: VD, FAZ area, FAZ circularity 52 . VD was defined as the total area of perfused vasculature per unit area in the measurement area, determined by summing the number of pixels with flow signal and dividing the total number by the number of pixels in the study area. Quantitative assessment of VD was obtained from the central (1 mm circle) of the ETDRS grid with automated measurements (Fig. 3B,E). FAZ area was defined as the area included within the FAZ boundary (Fig. 3C,F)

Data availability
The dataset generated and analyzed in this study are available from the corresponding author upon reasonable request.