Journal of Ophthalmology
Volume 2018 (2018), Article ID 1089043, 8 pages
https://doi.org/10.1155/2018/1089043

Research Article

Association between Visual Acuity and Retinal Layer Metrics in Diabetics with and without Macular Edema

LakshmiPriya Rangaraju,¹ Xuejuan Jiang,² J. Jason McAnany,¹ Michael R. Tan,¹ Justin Wanek,¹ Norman P. Blair,¹ Jennifer I. Lim,¹ and Mahnaz Shahidi²

¹Department of Ophthalmology and Visual Sciences, University of Illinois at Chicago, Chicago, IL, USA
²Department of Ophthalmology, University of Southern California, Los Angeles, CA, USA

Correspondence should be addressed to Mahnaz Shahidi; mshahidi@usc.edu

Received 13 April 2018; Revised 26 July 2018; Accepted 7 August 2018; Published 3 October 2018

Academic Editor: Alejandro Cerviño

Copyright © 2018 LakshmiPriya Rangaraju et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

Purpose. Diabetes is known to cause alterations in retinal microvasculature and tissue that progressively lead to visual impairment. Optical coherence tomography (OCT) is useful for assessment of total retinal thickening due to diabetic macular edema (DME). In the current study, we determined associations between visual acuity (VA) and retinal layer thickness, reflectance, and interface disruption derived from enface OCT images in subjects with and without DME. Materials and Methods. Best corrected VA was measured and high-density OCT volume scans were acquired in 149 diabetic subjects. A previously established image segmentation method identified retinal layer interfaces and locations of visually indiscernible (disrupted) interfaces. Enface thickness maps and reflectance images of the nerve fiber layer (NFL), combined ganglion cell and inner plexiform layer (GCLIPL), inner nuclear layer (INL), outer plexiform layer (OPL), outer nuclear layer (ONL), photoreceptor outer segment layer (OSL), and retinal pigment epithelium (RPE) were generated in the central macular subfield. The associations among VA and retinal layer metrics were determined by multivariate linear regressions after adjusting for covariates (age, sex, race, HbA1c, diabetes type, and duration) and correcting for multiple comparisons. Results. In DME subjects, increased GCLIPL and OPL thickness and decreased OSL thickness were associated with reduced VA. Furthermore, increased NFL reflectance and decreased OSL reflectance were associated with reduced VA. Additionally, increased areas of INL and ONL interface disruptions were associated with reduced VA. In subjects without DME, increased INL thickness was associated with reduced VA, whereas in subjects without DME but with previous antivascular endothelium growth factor treatment, thickening of OPL was associated with reduced VA. Conclusions. Alterations in retinal layer thickness and reflectance metrics derived from enface OCT images were associated with reduced VA with and without presence of DME, suggestive of their potential for monitoring development, progression, and treatment of DME.

1. Introduction

Diabetes is known to cause alterations in retinal microvasculature and tissue that progressively lead to visual impairment. Indeed, diabetic retinopathy (DR) is the leading cause of vision loss in working-age adults [1]. One consequence of DR is the development of diabetic macular edema (DME) due to the accumulation of fluid within the central retinal tissue, which is a major contributor towards vision loss [2]. Reduction in visual acuity (VA) with progression of DR based on fundus photography in individuals with DME has been established [3]. Optical coherence tomography (OCT) is currently the standard of clinical care for detecting abnormalities in retinal structure and quantifying the extent of retinal thickening due to DME [4]. Furthermore, high-resolution OCT can also quantify subtle retinal thickening not discernible on clinical examination in individuals with mild DME [5]. Additionally, previous studies have shown methods for 3D OCT imaging [6] and repeatable retinal layer thickness measurements in healthy and multiple sclerosis patients using commercially available OCT instruments (Cirrus HD-OCT and Spectralis SD-OCT) [7].

Although the total retinal thickness in the central subfield has been shown to be correlated with VA, the association was weaker than that of individual retinal layers [7, 8]. Specifically, thinning of the nerve fiber layer (NFL) in individuals with minimal or no DR and thickening of the inner nuclear layer (INL) and outer plexiform layer (OPL) in individuals with DME have been reported [9–15]. In addition to these changes, reduced VA has been correlated with thinning of the ganglion cell layer + inner plexiform layer (GCLIPL) and the photoreceptor outer segment layer (OSL) in subjects with and without DME [16, 17].

In addition to retinal thickness, alterations in retinal layer reflectance and interface have also been reported in DR. Specifically, in subjects with DME, reductions in photoreceptor outer segment length and disruptions of photoreceptor inner/outer segment junctions have been related to reduced VA [18, 19]. Additionally, inner retinal layer interface disruptions (visually indiscernible layer interfaces) and discontinuities in the inner segment/outer segment junction and in the external limiting membrane were associated with reduced VA [20–23]. However, most previous studies have examined retinal layer thickness, reflectance, and interface disruption from single OCT B-scans which limits localization of the spatial extent of retinal abnormalities. To better understand the spatial characteristics of retinal pathology, a method of retinal layer segmentation that generates three-dimensional outer retinal topography and reflectivity maps has been developed [24]. Although this approach provides more information than that obtained from single B-scans, the algorithms were based on images obtained using a spectral domain OCT prototype instrument. Using a commercially available OCT instrument, we have previously reported and validated methods for generating enface thickness maps and reflectance images of individual retinal layers from a high-density raster of OCT B-scan images and demonstrated alterations at different stages of DR [25–27]. In the current study, we identified individual retinal layers with thickness, reflectance, and interface disruption associated with VA in groups of DR subjects with and without DME.

2. Methods

2.1. Subjects

The research study was approved by an Institutional Review Board at the University of Illinois at Chicago and followed the tenets of the Declaration of Helsinki. Prior to enrollment, the study was explained to the subjects and informed consent was obtained. A total of 149 diabetic subjects participated in the study. All subjects underwent clinical examination by retinal specialists. Exclusion criteria were refractive error greater than 6 diopters of myopia, clinical diagnosis of glaucoma, age-related macular degeneration, retinal vascular occlusions or other conditions that could alter the anatomic integrity of the retina, history of intraocular surgery, cataract surgery performed less than 4 months prior to imaging, lens nuclear sclerosis score greater than 2+, or posterior subcapsular cataract concurrent with VA less than 20/20. One eye per subject was selected based on the exclusion criteria. If both eyes qualified, the eye with better image quality was included. Based on clinical examination by retina specialists, subjects’ eyes were classified as no DR (NDR; N = 51), nonproliferative DR (NPDR; N = 59), or proliferative DR (PDR; N = 39). The subjects were categorized into two subgroups, DME and no-DME, based on central subfield thickness (CST) being greater than 320 μm (males) and 304 μm (females) [28]. Twenty-eight subjects (NDR = 2; NPDR = 16; PDR = 10) had DME and 121 did not have DME at the time of imaging. Twenty-one of the 28 DME and 25 of the 121 no-DME subjects had previously received anti-VEGF therapy. Best-corrected VA was measured at a 4-meter distance using a retro-illuminated Early Treatment Diabetic Retinopathy Study (ETDRS) chart by an ophthalmic technician who was trained in the ETDRS protocol.

2.2. Image Acquisition

Spectral domain OCT (SDOCT) imaging of a retinal area of 20° × 15° centered on the fovea was performed using a commercially available instrument (Spectralis; Heidelberg Engineering, Heidelberg, Germany). A high-density SDOCT raster volume scan was generated from 73 raster horizontal B-scans (9 averaged frames) with a vertical spacing of 62 μm. Each B-scan consisted of 1024 A-scans and had a depth resolution of 3.9 μm.

2.3. Image Analysis

SDOCT B-scans were analyzed using our previously described automated image segmentation software based on graph theory and dynamic programming [25]. In brief, a graph was generated for all SDOCT B-scans with edge weights designated according to vertical gradients in the images. A horizontal path through the graph that minimized the total sum of the weights was found using Dijkstra’s algorithm and defined a line separating two retinal cell layers. Figure 1 displays eight retinal interfaces that were detected: (1) vitreous and NFL, (2) NFL and GCLIPL, (3) GCLIPL and INL, (4) INL and OPL, (5) OPL and outer nuclear layer (ONL), (6) ONL and OSL, (7) OSL and retinal pigment epithelium (RPE), and (8) RPE and choroid. Following our previously reported method, an operator was able to review the automated interface detection results by scrolling through all SDOCT B-scans. If necessary, errors were corrected by manually selecting the segmentation line that required adjustment and drawing a revised line corresponding to the visualized cell-layer interface. The program then regenerated an automated line by restricting the graph search area to a small vertical image region around the manually drawn line and recalculating the minimum graph cut solution, as previously described [25]. Locations of retinal layer interfaces that were not visually discernible (disruptions) due to gross abnormalities in retinal cell layer architecture were manually selected by the operator. The automated image segmentation method was previously validated [25] by demonstrating a high correlation with data provided by commercial instruments and also showing a decrease in retinal thickness with increased age, consistent with previous reports. Additionally, the average error rate obtained by the automated segmentation method was shown to be 7% in NPDR subjects [25].

Figure 1: Example of an OCT B-scan in a DR subject without DME, representing the eight segmented retinal layer interfaces.

Enface thickness maps and reflectance images were generated for each of 7 retinal layers (NFL, GCLIPL, INL, OPL, ONL, OSL, and RPE) based on segmentation of the 8 retinal interfaces in the SDOCT B-scans. Regions containing retinal layer interface disruption were not assigned thickness or reflectance values. Retinal layer metrics of thickness, reflectance, and areas of interface disruption were evaluated in the ETDRS central subfield (1 mm diameter) [29]. Mean thickness metrics were calculated for each layer (NFLT, GCLIPLT, INLT, OPLT, ONLT, and OSLT). Reflectance ratio metrics were also calculated in the central subfield for each layer (NFLR, GCLIPLR, INLR, OPLR, ONLR, and OSLR) as the mean intensity of each layer divided by the mean intensity of the RPE (layerIntensity/RPEIntensity).

Percent areas of layer interface disruption relative to the total central ETDRS subfield area (NFLd, INLd, ONLd, and RPEd) were calculated. NFLd was calculated based on the interface disruption of both vitreous/NFL and NFL/GCLIPL interfaces relative to the total central subfield area. Similarly, INLd was calculated based on disruptions of both GCLIPL/INL and INL/OPL interfaces; ONLd was calculated based on disruptions of both OPL/ONL and ONL/OSL interfaces; RPEd was calculated based on disruptions of both OSL/RPE and RPE/choroid interfaces.

2.4. Statistical Analysis

Associations of retinal layer thickness, reflectance, and percent area of interface disruptions were determined by multivariate general linear models, adjusting for age, sex, race, diabetes type, diabetes duration, and HbA1C level. Diagnostics of model assumptions were performed on all models and residuals for VA followed a normal distribution. Potential influential data points were also identified; sensitivity analyses excluding those data points were performed and the regression results were not changed materially. With a sample size of 28, the statistical power for multivariate linear regression (with 7 variables) was 80% to detect a partial correlation of 0.56 or greater, and with a sample size of 121, the power was 80% to detect a partial correlation of 0.26 or greater. For the associations between VA and retinal layer interface disruptions, Kruskal–Wallis tests were also performed and results were similar to those obtained by analysis of variance. All statistical tests were conducted using SAS 9.4 (SAS Institute Inc., Cary, NC). All P-values were from two-sided tests. For the associations between VA and layer-specific measurements, significance was accepted at to account for multiple comparisons using Bonferroni correction.

3. Results

Demographic details and ocular characteristics of subjects are presented in Table 1. Mean CST in all subjects was 283 ± 56 µm (N = 149). As expected, CST was greater in subjects with DME (371 ± 70 µm; N = 28), compared with subjects without DME (263 ± 23 µm; N = 121) (P < 0.001). Log MAR VA of all subjects was 0.04 ± 0.13. Subjects with DME had worse VA (0.14 ± 0.16 log MAR), compared with subjects with no DME (0.02 ± 0.11 log MAR) (P < 0.001). Examples of enface retinal layer thickness maps and reflectance images (not normalized to RPE reflectance) in a DR subject with DME are displayed in Figure 2. On average, ONL and NFL have the largest and smallest thickness, respectively. Due to the foveal depression, GCLIPL has minimal thickness in the center of the fovea. The ONL and OSL have the lowest and highest reflectance, respectively. Regions of retinal layer interface disruptions are represented by black and yellow on the thickness and reflectance layer maps, respectively.

Table 1: Study population and ocular characteristics.

Total sample size = 149

Population characteristics

Age, mean (SD), years
56.0 (11.8)

Female sex, number (%)
87 (58.4%)

Race/ethnicity, number (%)

    White
25 (16.8%)

    African American
77 (51.7%)

    Hispanic or Latino
42 (28.2%)

    Asian
5 (3.4%)

DM type, number (%)

    1
12 (8.1%)

    2
137 (91.9%)

DM duration, mean (SD), years
16.3 (10.1)

Glycated hemoglobin level, mean (SD), %
7.9 (1.7)

Ocular characteristics

Study eye, number (%)

    Right
96 (64%)

    Left
53 (36%)

Spherical equivalent refractive error, mean (SD), D
−0.66 (1.74)

Visual acuity, mean (SD), log MAR units
0.04 (0.13)

Diabetic retinopathy stage, number (%)

    No DR
51 (34%)

    NPDR
59 (40%)

    PDR
39 (27%)

Presence of DME, number (%)
28 (18.8%)

Central subfield thickness, mean (SD), μm
283 (56)

DM = diabetes mellitus; DR = diabetic retinopathy; PDR = proliferative diabetic retinopathy; NPDR = nonproliferative diabetic retinopathy; DME = diabetic macular edema; SD = standard deviation.

Figure 2: Enface thickness maps and reflectance images in a DR subject with DME. Top (left to right): ETDRS central subfield thickness maps of NFL, GCLIPL, INL, OPL, ONL, OSL, and RPE retinal layers. Bottom (left to right): reflectance images of NFL, GCLIPL, INL, OPL, ONL, OSL, and RPE retinal layers. Regions of retinal layer interface disruptions are represented by black and yellow on retinal layer thickness and reflectance maps, respectively. In these regions of disrupted interfaces, thickness and reflectance values were not assigned.

3.1. Retinal Layer Thickness

Mean retinal layer thickness and regression coefficients obtained in DME subjects after adjusting for age, sex, race, diabetes type, diabetes duration, and HbA1c are listed in Table 2. Increased GCLIPLT and OPLT (coefficients ≥ 0.03 log MAR/10 micron) and decreased OSLT (coefficient = −0.117 log MAR/10 micron) were associated with reduced VA (P < 0.008). That is, each 10 micron increase in GCLIPLT and OPLT resulted in an approximate 1 or 2 letter loss of VA, whereas each 10 micron decrease in OSLT resulted in approximately 5 letters of VA loss. Mean retinal layer thickness and regression coefficients obtained in no-DME subjects after adjusting for age, sex, race, diabetes type, diabetes duration, and HbA1c are listed in Table 3. Increased INLT was associated with reduced VA (coefficient = 0.048 log MAR/10 micron). Mean retinal layer thickness and regression coefficients obtained in no-DME subjects with a history of anti-VEGF treatment after adjusting for age, sex, race, diabetes type, diabetes duration, and HbA1c are listed in Table 4. Increased OPLT was associated with reduced VA (coefficient = 0.061 log MAR/10 micron).

Table 2: Multivariate associations of retinal layer thickness and reflectance ratio with logMAR visual acuity based on data from DME subjects (N = 28).

Metrics
Thickness (microns), mean ± SD
Regression coefficient (95% CI) [1]

NFLT
21 ± 6
0.074 (−0.020, 0.167)

GCLIPLT
43 ± 26
0.030 (0.013, 0.048)

INLT
55 ± 57
0.014 (0.003, 0.024)

OPLT
24 ± 8
0.101 (0.043, 0.160)

ONLT
159 ± 57
−0.007 (−0.018, 0.003)

OSLT
42 ± 9
−0.117 (−0.165, −0.069)

Metrics
Reflectance ratio, mean ± SD
Regression coefficient (95% CI) [1]

NFLR
0.63 ± 0.12
0.077 (0.040, 0.114)

GCLIPLR
0.65 ± 0.07
0.085 (−0.003, 0.173)

INLR
0.49 ± 0.08
−0.028 (−0.123, 0.068)

OPLR
0.53 ± 0.08
0.057 (−0.035, 0.150)

ONLR
0.35 ± 0.04
0.012 (−0.130, 0.154)

OSLR
0.76 ± 0.11
−0.069 (−0.117, −0.021)

. ¹Regression coefficients represent logMAR change with 10 micron increase in thickness or 0.1 increase in reflectance ratio, after adjusting for age, sex, race, diabetes type, diabetes duration, and HbA1c. NFL = nerve fiber layer; GCLIPL = ganglion cell layer + inner plexiform layer; INL = inner nuclear layer; OPL = outer plexiform layer; ONL = outer nuclear layer; OSL = photoreceptor outer segment layer; RPE = retinal pigment epithelium. Suffix T = thickness; suffix R = reflectance ratio.

Table 3: Multivariate associations of retinal layer thickness and reflectance ratio with logMAR visual acuity based on data from subjects without DME (N = 121).

Metrics
Thickness (microns), mean ± SD
Regression coefficient (95% CI) [1]

NFLT
18 ± l3
0.004 (−0.060, 0.068)

GCLIPLT
28 ± 12
−0.003 (−0.021, 0.016)

INLT
21 ± 8
0.048 (0.023, 0.073)

OPLT
21 ± 8
0.029 (0.001, 0.056)

ONLT
114 ± 22
−0.002 (−0.012, 0.007)

OSLT
44 ± 7
−0.008 (−0.036, 0.019)

Metrics
Reflectance ratio, mean ± SD
Regression coefficient (95% CI) [1]

NFLR
0.54 ± 0.07
0.029 (0, 0.057)

GCLIPLR
0.61 ± 0.07
0.021 (−0.010, 0.052)

INLR
0.50 ± 0.06
0.010 (−0.025, 0.044)

OPLR
0.52 ± 0.06
0.009 (−0.026, 0.045)

ONLR
0.37 ± 0.05
0.043 (0.001, 0.085)

OSLR
0.81 ± 0.07
−0.007 (−0.037, 0.023)

. ¹Regression coefficients represent logMAR change with 10 micron increase in thickness or 0.1 increase in reflectance ratio, after adjusting for age, sex, race, diabetes type, diabetes duration, and HbA1c. NFL = nerve fiber layer; GCLIPL = ganglion cell layer + inner plexiform layer; INL = inner nuclear layer; OPL = outer plexiform layer; ONL = outer nuclear layer; OSL = photoreceptor outer segment layer; RPE = retinal pigment epithelium. Suffix T = thickness; suffix R = reflectance ratio.

Table 4: Multivariate associations of retinal layer thickness and reflectance ratio with logMAR visual acuity based on data from subjects without DME but with a history of anti-VEGF treatment (N = 25).

Metrics
Thickness (microns), mean ± SD
Regression coefficient (95% CI) [1]

NFLT
18 ± 3
0.132 (−0.054, 0.318)

GCLIPLT
29 ± 11
0.012 (−0.034, 0.057)

INLT
28 ± 10
0.046 (0.007, 0.086)

OPLT
21 ± 12
0.061 (0.028, 0.093)

ONLT
117 ± 39
−0.009 (−0.02, 0.002)

OSLT
41 ± 8
−0.021 (−0.078, 0.036)

Metrics
Reflectance ratio, mean ± SD
Regression coefficient (95% CI) [1]

NFLR
0.55 ± 0.09
0.031 (−0.037, 0.098)

GCLIPLR
0.61 ± 0.08
0.017 (−0.056, 0.09)

INLR
0.49 ± 0.07
−0.002 (−0.079, 0.076)

OPLR
0.51 ± 0.07
−0.014 (−0.09, 0.061)

ONLR
0.38 ± 0.06
0.01 (−0.068, 0.088)

OSLR
0.79 ± 0.07
0.003 (−0.079, 0.085)

. ¹Regression coefficients represent logMAR change with 10 micron increase in thickness or 0.1 increase in reflectance ratio, after adjusting for age, sex, race, diabetes type, diabetes duration, and HbA1c. NFL = nerve fiber layer; GCLIPL = ganglion cell layer + inner plexiform layer; INL = inner nuclear layer; OPL = outer plexiform layer; ONL = outer nuclear layer; OSL = photoreceptor outer segment layer; RPE = retinal pigment epithelium. Suffix T = thickness; suffix R = reflectance ratio.

3.2. Retinal Layer Reflectance

Mean retinal layer reflectance and regression coefficients obtained in DME subjects after adjusting for age, sex, race, diabetes type, diabetes duration, and HbA1c are listed in Table 2. Increased NFLR (coefficient = 0.077 log MAR/0.1 reflectance ratio) and decreased OSLR (coefficient = −0.069 log MAR/0.1 reflectance ratio) were associated with reduced VA (P < 0.008). Mean retinal layer reflectance and regression coefficients obtained in no-DME subjects after adjusting for age, sex, race, diabetes type, diabetes duration, and HbA1c are listed in Table 3. Retinal layer reflectance changes were not associated with reduced VA. Similarly, in no-DME subjects with a history of anti-VEGF treatment, there was no significant association between retinal layer reflectance and VA (Table 4).

3.3. Retinal Interface Disruption

Both increased INLd and ONLd were associated with worse VA after adjusting for age, sex, race, diabetes type, diabetes duration, and HbA1c (). The mean VA values for subjects with 0%, <10%, and ≥10% INLd were 0.02 ± 0.12 (N = 124), 0.09 ± 0.13 (N = 13), and 0.18 ± 0.14 (N = 12) log MAR, respectively (). The mean VA values in subjects with 0%, <10%, and ≥10% ONLd were 0.02 ± 0.12 (N = 125), 0.11 ± 0.15 (N = 11), and 0.16 ± 0.14 (N = 13) log MAR, respectively (). When the analyses were conducted separately for DME and no-DME subjects, there was no significant association between VA and INLd or ONLd ().

4. Discussion

In the current study, we determined associations between VA and retinal layer thickness and reflectance in DR subjects with and without DME. In DME subjects, GCLIPL and OPL thickening and OSL thinning, as well as hyper-reflectance of the NFL and hypo-reflectance of the OSL, were associated with reduced VA. Moreover, VA was reduced in regions of disrupted retinal interfaces bordering the INL and ONL. Interestingly, in subjects without DME, thickening of the INL was associated with reduced VA, whereas in the subset of these subjects who had previous anti-VEGF treatment, thickening of OPL was associated with reduced VA.

In DME subjects, increased GCLIPLT and OPLT were associated with decreased VA. The relationship between increased OPLT and decreased VA is consistent with previous studies, which reported a correlation between increased edema and reduced VA [20]. Our findings are also consistent with other reports of the presence of cystoid spaces or increased OPL thickness, although the relation to VA was not reported in these studies [9, 11]. Decreased OSLT was associated with reduced VA, consistent with the findings of a previous study [18]. This change in OSLT suggests that photoreceptor degeneration, possibly secondary to macular edema, may contribute to vision loss. Further studies are needed to confirm our findings, given the limited number of DME subjects in this study, and to determine whether treatment based on thickening of specific retinal layers can potentially improve visual outcomes.

In subjects without DME, increased INLT was associated with reduced VA. This finding suggests that thickening of the INL can occur and perhaps precede the development of DME. Of note, the regression coefficient for the association between INLT and VA in no-DME subjects with a history of anti-VEGF treatment was similar to the value in subjects without DME, though it did not reach statistical significance after adjusting for multiple comparisons, possibly due to the smaller sample size. Interestingly, in no-DME subjects with a history of anti-VEGF treatment, thickening of OPL was associated with reduced VA, similar to the finding in DME subjects. Future longitudinal studies are needed to investigate thickness alterations in this retinal layer as an early indicator for recurrent DME or as a marker for the adequacy of DME treatment.

The current study demonstrated associations between changes in retinal layer reflectance and reduced VA. In the DME subjects, decreased OSLR and increased NFLR were associated with reduced VA. The association of OSLR and VA is consistent with a previous study that reported reduced reflectance of the photoreceptor layer when cystoid spaces were present in the OPL [30] and other studies that showed a correlation between the continuity of the photoreceptor inner and outer segment (ellipsoid) and VA [7, 19–21, 31, 32]. However, these previous studies evaluated the inner and outer segment junction visibility based on a single B-scan or a low-density raster of B-scans. In the current study, enface reflectance images of the OSL were obtained from high-density OCT B-scans, thereby providing a more accurate localization of the spatial extent of reduced OSL reflectance. To our knowledge, an association between increased NFLR and reduced VA has not been previously reported.

Several previous studies have reported methods for segmentation of different retinal layers in DR [33–39]. The accuracy of automated segmentation of retinal layers in healthy and DR subjects relies on the presence of distinct interfaces between layers. However, in DME subjects who have severe pathology, these interfaces cannot be clearly identified, even by expert human evaluation. The method presented in the current study allowed quantitative measurement of areas in which retinal layer interfaces were visually indiscernible due to pathologies and macular edema, thus providing a useful metric for assessing retinal integrity. The finding of an association between reduced VA and disrupted INL and ONL interfaces is consistent with previous studies that reported a correlation between combined inner retinal interface disruptions (NFL/GCLIPL, GCLIPL/INL/OPL, and OPL/ONL) and VA, although they did not evaluate disruptions in other layer interfaces (ONL/OSL and OSL/RPE) [22, 23].

5. Conclusion

Quantitative assessment of retinal layer integrity by enface OCT imaging may be clinically relevant for monitoring the progression of pathologies due to disease or their resolution following treatment. Concurrent assessment of thickness, reflectance, and interface disruption of individual retinal layers by enface OCT imaging provides a comprehensive approach for identifying anatomic outcomes that may be useful for monitoring development, progression, and treatment efficacy of DME.

Data Availability

The data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

This work was supported by the NIH grants DK010439 (MS), EY001792 (UIC Core), and EY026004 (JM), a Senior Scientific Investigator Award (MS), a Dolly Green Special Scholar Award (JM), and an unrestricted departmental grant from Research to Prevent Blindness.

References

S. Wild, A. Green, R. Sicree, and H. King, “Global prevalence of diabetes: estimates for the year 2000 and projections for 2030,” Diabetes Care, vol. 27, no. 5, pp. 1047–1053, 2004.

R. Lee, T. Y. Wong, and C. Sabanayagam, “Epidemiology of diabetic retinopathy, diabetic macular edema and related vision loss,” Eye and Vision, vol. 2, no. 1, p. 17, 2015.

B. Bengtsson, A. Heijl, and E. Agardh, “Visual fields correlate better than visual acuity to severity of diabetic retinopathy,” Diabetologia, vol. 48, no. 12, pp. 2494–2500, 2005.

L. M. Sakata, J. Deleon-Ortega, V. Sakata, and C. A. Girkin, “Optical coherence tomography of the retina and optic nerve-a review,” Clinical & Experimental Ophthalmology, vol. 37, no. 1, pp. 90–99, 2009.

J. C Brown, S. D. Solomon, S. B. Bressler, A. P. Schachat, C. DiBernardo, and N. M. Bressler, “Detection of diabetic foveal edema: contact lens biomicroscopy compared with optical coherence tomography,” Archives of Ophthalmology, vol. 122, no. 3, pp. 330–335, 2004.

C. Hitzenberger, P. Trost, P. W. Lo, and Q. Zhou, “Three-dimensional imaging of the human retina by high-speed optical coherence tomography,” Optics Express, vol. 11, no. 21, pp. 2753–2761, 2003.

T. Otani, Y. Yamaguchi, and S. Kishi, “Correlation between visual acuity and foveal microstructural changes in diabetic macular edema,” Retina, vol. 30, no. 5, pp. 774–780, 2010.

D. J. Browning, A. R. Glassman, L. P. Aiello et al., “Relationship between optical coherence tomography-measured central retinal thickness and visual acuity in diabetic macular edema,” Ophthalmology, vol. 114, no. 3, pp. 525–536, 2007.

F. Bandello, A. N. Tejerina, S. Vujosevic et al., “Retinal layer location of increased retinal thickness in eyes with subclinical and clinical macular edema in diabetes type 2,” Ophthalmic Research, vol. 54, no. 3, pp. 112–117, 2015.

J. Chhablani, A. Sharma, A. Goud et al., “Neurodegeneration in type 2 diabetes: evidence from spectral-domain optical coherence tomography,” Investigative Ophthalmology & Visual Science, vol. 56, no. 11, pp. 6333–6338, 2015.

T. Murakami and N. Yoshimura, “Structural changes in individual retinal layers in diabetic macular edema,” Journal of Diabetes Research, vol. 2013, Article ID 920713, 11 pages, 2013.

P. H. Peng, H. S. Lin, and S. Lin, “Nerve fibre layer thinning in patients with preclinical retinopathy,” Canadian Journal of Ophthalmology, vol. 44, no. 4, pp. 417–422, 2009.

H. W. Van Dijk, F. D. Verbraak, P. H. Kok et al., “Early neurodegeneration in the retina of type 2 diabetic patients,” Investigative Opthalmology & Visual Science, vol. 53, no. 6, pp. 2715–2719, 2012.

A. Verma, R. Raman, K. Vaitheeswaran et al., “Does neuronal damage precede vascular damage in subjects with type 2 diabetes mellitus and having no clinical diabetic retinopathy?” Ophthalmic Research, vol. 47, no. 4, pp. 202–207, 2012.

S. Vujosevic and E. Midena, “Retinal layers changes in human preclinical and early clinical diabetic retinopathy support early retinal neuronal and Muller cells alterations,” Journal of Diabetes Research, vol. 2013, Article ID 905058, 8 pages, 2013.

T. Alasil, P. A. Keane, J. F. Updike et al., “Relationship between optical coherence tomography retinal parameters and visual acuity in diabetic macular edema,” Ophthalmology, vol. 117, no. 12, pp. 2379–2386, 2010.

S. Bonnin, R. Tadayoni, A. Erginay, P. Massin, and B. Dupas, “Correlation between ganglion cell layer thinning and poor visual function after resolution of diabetic macular edema,” Investigative Ophthalmology & Visual Science, vol. 56, no. 2, pp. 978–982, 2015.

F. Forooghian, P. F. Stetson, S. A. Meyer et al., “Relationship between photoreceptor outer segment length and visual acuity in diabetic macular edema,” Retina, vol. 30, no. 1, pp. 63–70, 2010.

A. S. Maheshwary, S. F. Oster, R. M. Yuson, L. Cheng, F. Mojana, and W. R. Freeman, “The association between percent disruption of the photoreceptor inner segment-outer segment junction and visual acuity in diabetic macular edema,” American Journal of Ophthalmology, vol. 150, no. 1, pp. 63.e61–67.e61, 2010.

T. Murakami, K. Nishijima, A. Sakamoto, M. Ota, T. Horii, and N. Yoshimura, “Association of pathomorphology, photoreceptor status, and retinal thickness with visual acuity in diabetic retinopathy,” American Journal of Ophthalmology, vol. 151, no. 2, pp. 310–317, 2011.

H. J. Shin, S. H. Lee, H. Chung, and H. C. Kim, “Association between photoreceptor integrity and visual outcome in diabetic macular edema,” Graefe’s Archive for Clinical and Experimental Ophthalmology, vol. 250, no. 1, pp. 61–70, 2012.

J. K. Sun, M. M. Lin, J. Lammer et al., “Disorganization of the retinal inner layers as a predictor of visual acuity in eyes with center-involved diabetic macular edema,” JAMA Ophthalmology, vol. 132, no. 11, pp. 1309–1316, 2014.

J. K. Sun, S. H. Radwan, A. Z. Soliman et al., “Neural retinal disorganization as a robust marker of visual acuity in current and resolved diabetic macular edema,” Diabetes, vol. 64, no. 7, pp. 2560–2570, 2015.

M. Wojtkowski, B. L. Sikorski, I. Gorczynska et al., “Comparison of reflectivity maps and outer retinal topography in retinal disease by 3-D Fourier domain optical coherence tomography,” Optics Express, vol. 17, no. 5, pp. 4189–4207, 2009.

J. Wanek, N. P. Blair, F. Y. Chau, Y. I. Lim, Y. I. Leiderman, and M. Shahidi, “Alterations in retinal layer thickness and reflectance at different stages of diabetic retinopathy by enface optical coherence tomography,” Investigative Ophthalmology & Visual Science, vol. 57, no. 9, pp. OCT341–OCT347, 2016.

F. Mohammad, J. Wanek, R. Zelkha, J. I. Lim, J. Chen, and M. Shahidi, “A method for en face OCT imaging of subretinal fluid in age-related macular degeneration,” Journal of Ophthalmology, vol. 2014, Article ID 720243, 6 pages, 2014.

J. Wanek, R. Zelkha, J. I. Lim, and M. Shahidi, “Feasibility of a method for en face imaging of photoreceptor cell integrity,” American Journal of Ophthalmology, vol. 152, no. 5, pp. 807.e1–814.e1, 2011.

K. V. Chalam, S. B. Bressler, A. R. Edwards et al., “Retinal thickness in people with diabetes and minimal or no diabetic retinopathy: Heidelberg Spectralis optical coherence tomography,” Investigative Ophthalmology & Visual Science, vol. 53, no. 3, pp. 8154–8161, 2012.

ETDRS Research Group, “Early Treatment Diabetic Retinopathy Study design and baseline patient characteristics. ETDRS report number 7,” Ophthalmology, vol. 98, no. 5, pp. 741–756, 1991.

T. Murakami, K. Nishijima, T. Akagi et al., “Optical coherence tomographic reflectivity of photoreceptors beneath cystoid spaces in diabetic macular edema,” Investigative Ophthalmology & Visual Science, vol. 53, no. 3, pp. 1506–1511, 2012.

A. Sakamoto, K. Nishijima, M. Kita, H. Oh, A. Tsujikawa, and N. Yoshimura, “Association between foveal photoreceptor status and visual acuity after resolution of diabetic macular edema by pars plana vitrectomy,” Graefe’s Archive for Clinical and Experimental Ophthalmology, vol. 247, no. 10, pp. 1325–1330, 2009.

Y. Shen, K. Liu, and X. Xu, “Correlation between visual function and photoreceptor integrity in diabetic macular edema: spectral-domain optical coherence tomography,” Current Eye Research, vol. 41, no. 3, pp. 391–399, 2016.

F. Mohammad, R. Ansari, J. Wanek, A. Francis, and M. Shahidi, “Feasibility of level-set analysis of enface OCT retinal images in diabetic retinopathy,” Biomedical Optics Express, vol. 6, no. 5, pp. 1904–1918, 2015.

Y. Huang, R. P. Danis, J. W. Pak et al., “Development of a semi-automatic segmentation method for retinal OCT images tested in patients with diabetic macular edema,” PLoS One, vol. 8, no. 12, Article ID e82922, 2013.

M. K. Garvin, M. D. Abramoff, X. Wu, S. R. Russell, T. L. Burns, and M. Sonka, “Automated 3-D intraretinal layer segmentation of macular spectral-domain optical coherence tomography images,” IEEE Transactions on Medical Imaging, vol. 28, no. 9, pp. 1436–1447, 2009.

K. A. Vermeer, J. Van Der Schoot, H. G. Lemij, and J. F. De Boer, “Automated segmentation by pixel classification of retinal layers in ophthalmic OCT images,” Biomedical Optics Express, vol. 2, no. 6, pp. 1743–1756, 2011.

S. J. Chiu, M. J. Allingham, P. S. Mettu, S. W. Cousins, J. A. Izatt, and S. Farsiu, “Kernel regression based segmentation of optical coherence tomography images with diabetic macular edema,” Biomedical Optics Express, vol. 6, no. 4, pp. 1172–1194, 2015.

S. J. Chiu, X. T. Li, P. Nicholas, C. A. Toth, J. A. Izatt, and S. Farsiu, “Automatic segmentation of seven retinal layers in SDOCT images congruent with expert manual segmentation,” Optics Express, vol. 18, no. 18, pp. 19413–19428, 2010.

S. J. Chiu, C. A. Toth, C. Bowes Rickman, J. A. Izatt, and S. Farsiu, “Automatic segmentation of closed-contour features in ophthalmic images using graph theory and dynamic programming,” Biomedical Optics Express, vol. 3, no. 5, pp. 1127–1140, 2012.