Association Between Alterations of the Choriocapillaris Microcirculation and Visual Function and Cone Photoreceptors in Patients With Diabetes

Abstract

Purpose: The purpose of this study was to investigate the association between the choriocapillaris microcirculation and the visual function and cone photoreceptor structure in patients with diabetes.

Methods: Thirteen control subjects and 26 patients with type 2 diabetes were recruited. The patients with diabetes were divided into three groups based on the grade of diabetic retinopathy (DR). The retinal sensitivity (RS) was evaluated using microperimetry. Cone photoreceptor mosaics were imaged using an adaptive optics retinal camera, and the cone heterogeneity packing index (HPi) was calculated. Optical coherence tomography angiography (OCTA) images of the choriocapillaris were obtained using spectral-domain OCTA, and the area of flow deficit (FD) was evaluated. All parameters were evaluated in the foveal and parafoveal areas.

Results: The study included four patients with diabetes without retinopathy, 12 patients with nonproliferative diabetic retinopathy (NPDR), and 10 patients with proliferative diabetic retinopathy (PDR). The foveal and parafoveal FDs were correlated significantly (fovea, r = −0.58; P = 0.046 and r = −0.82; P = 0.003; parafovea, r = −0.59; P = 0.044 and r = −0.72; and P = 0.019, respectively) with the RS in patients with NPDR and PDR, but not in control and no diabetic retinopathy (NDR) groups. There were no differences in the foveal HPi among the groups.

Conclusions: Impaired choriocapillaris microcirculation is associated with impaired visual function but not cone photoreceptor integrity in eyes with DR.

Diabetic retinopathy (DR), clinically defined as a microangiopathy in the retina, is a leading cause of visual loss.¹ Microangiopathy in diabetic eyes affects two microcirculatory systems that supply oxygen and nutrients to the neural retinal tissue: the retinal circulation, which supplies the inner retinal layer and about 10% to 15% of the outer retina,² and the choroidal circulation, which mainly covers the outer retinal layer, including the photoreceptor layer through the choriocapillaris.³ In addition to the traditional idea that retinal circulatory disruption is the main culprit in visual functional loss, accumulating evidence suggests that the choroidal circulation also may be affected in diabetic eyes and contribute to the outer retinal damage. Histologic studies have shown choriocapillaris abnormalities, such as vessel dropout, ischemia, and aneurysms, in diabetic eyes.⁴ In addition, in eyes with DR, disruption of the choroidal circulation is accompanied by loss of the ellipsoid zone,⁵ which is thought to represent the mitochondria in the photoreceptors.⁶ Therefore, impaired choroidal microcirculation also may affect the photoreceptor structures/functions in diabetic eyes. Taken together, choroidal circulatory disruption likely causes visual functional loss.^7,8 However, the association between the choroidal circulatory disturbances and the photoreceptor structure has not been investigated quantitatively, especially in diabetic eyes. 

Thus far, no studies have used a combination of high-resolution imaging modalities and high-spatial resolution functional studies to investigate the association between the choroidal circulation and visual function, although more attention is starting to be paid to the choroidal circulation. Optical coherence tomography angiography (OCTA), for example, shows the blood flow and density in the retina and choriocapillaris quantitatively. Using OCTA, clinical studies have reported structural alterations of the choroidal vasculature in a quantitative manner in eyes with DR^9–13 using an area of flow deficit (FD), which is an index of capillary loss in the choriocapillaris layer,^14–16 or a similar index.⁵ In addition, assessment of the regularity of the cone photoreceptors became possible with the use of an adaptive optic (AO) retinal imaging system, which can observe the cone photoreceptor mosaic in near histologic resolution. The heterogeneity packing index (HPi) can serve as a sensitive biomarker for abnormalities in the regularity of the cone photoreceptors.^17,18 Several studies have reported changes in the cone mosaic in patients with diabetes.¹⁷ Finally, the local macular visual function can be assessed reliably using microperimetry combined with fundus imaging, as this enables precise measurement as the retinal sensitivity (RS).^7,19 A previous study reported that eyes with DR had reduced RS in the area of retinal ischemia.²⁰ 

In the current study, to gain insight into the potential association among the choroidal microcirculation and the visual function and photoreceptor structures in patients with diabetes, we conducted a prospective study and investigated the aforementioned parameters in various stages of diabetes. 

Methods

Patient Recruitment

Patients with diabetes and healthy volunteers were recruited from the Department of Ophthalmology at Asahikawa Medical University Hospital between October 2018 and April 2020. The Institutional Review Board of Asahikawa Medical University approved this cross-sectional study, which adhered to the tenets of the Declaration of Helsinki. All patients provided written informed consent to participate in this research. 

Study Population and Examinations

This study included consecutive patients with type 2 diabetes mellitus that fulfilled the diagnosis based on the American Diabetes Association criteria.²¹ After the pupil was dilated with a drop containing tropicamide and phenylephrine hydrochloride (Santen Pharmaceutical Co., Osaka, Japan), each patient underwent a standard ophthalmic examination, including color fundus imaging using Optos 200Tx (Nikon, Tokyo, Japan), infrared fundus imaging at 30 × 30 degrees by Spectralis spectral-domain optical coherence tomography (OCT; Heidelberg Engineering GmbH, Heidelberg, Germany), spectral-domain OCT (RTVue-XR Avanti, Optovue Inc., Fremont, CA, USA), microperimetry, AO retinal imaging, and OCTA. For healthy volunteers, all examinations were performed except for AO retinal imaging. 

We excluded eyes with other retinal diseases, a history of vitrectomy or laser treatment of the macula, panretinal photocoagulation (PRP) treatment within 12 weeks, antivascular endothelial growth factor (VEGF) therapy within 6 months of the study, high myopia (>6 diopters), astigmatism (>3 diopters), long axial length (>26.0 mm), significant media opacity, movement or shadow artifacts in OCTA images, or diabetic macular edema within a 3 × 3-mm macular area, which was evaluated using spectral-domain OCT. 

An experienced ophthalmologist (A.S.) graded the DR images in a masked fashion based on the Early Treatment of Diabetic Retinopathy Study Severity Scale²² and patients were divided into three groups: no diabetic retinopathy (NDR), nonproliferative diabetic retinopathy (NPDR), and proliferative diabetic retinopathy (PDR). When both eyes of a patient had the same stage of DR, the right eye was enrolled, and if the eyes had different stages of DR, the eye with the more advanced stage was included. 

Microperimetry

The retinal function was assessed based on the RS values obtained by microperimetry (MAIA, Centervue, Padova, Italy). Twelve points were measured at both 1 degree (foveal area) and 3.5 degrees (parafoveal area) from the central fixation target (Fig. 1). The microperimetry examination was performed in a darkened room using the following parameters: Goldman III stimuli size, a background luminance corresponding to 4asb, and the 4-2 threshold strategy. 

Figure 1.

 

Evaluation pattern in microperimetry. The retinal sensitivity (RS) is measured at 24 points. The foveal and parafoveal RS values are defined as the average of the 12 points at 1 and 3.5 degrees that corresponded to 0.286 mm* and 1.001 mm† from the central fixation target, respectively.

AO Retinal Image Processing and Analysis

Cone mosaic images were acquired using the rtx1-e AO retinal flood-illumination camera (Imagine Eyes, Orsay, France). This AO camera uses infrared illumination (wavelength, 850 nm) and has a resolution of 3.5 µm, and the field of view is 4 × 4 degrees. Nine overlapping images were acquired to create a montage covering an area of 10 × 10 degrees from the central retina. The MosaicJ plugin (ImageJ, version 1.52, Image Processing and Analysis in Java; National Institutes of Health, Bethesda, MD, USA) was used to manually create a wide-field montage of the cone mosaic images²³ (Fig. 2A). At each microperimetry measurement point (Figs. 2B–D), we assessed the cone packing arrangement using Voronoi diagrams (Fig. 2F) with AOdetect software (AOdetect Mosaic XL, version 2.0 b18, Imagine Eyes) provided by the manufacturer. The algorithm successively uses Gaussian background subtraction, adaptive denoising, multiscale digital filtering, and regional maxima calculation to compute estimates of the signal level and position of any small and bright round spots in the image. After this initial segmentation, the spots that fulfilled both conditions were selected: their signal was significantly higher than the photon noise, and they were not too close (i.e. <3.2 µm) to spots with a higher signal level. These conditions were intended to prevent detection of rods near cones. Once the segmentation was complete, the software analyzed the spatial distribution of the detected cells and computed several metrics, including inter-cell spacing, cell density, and the number of nearest neighbors using classical Delaunay triangulation and Voronoi tessellation algorithms. First, we chose a sampling window (Fig. 2D) of 0.35 × 0.35 degrees equal to each microperimetry spot size. Second, the AOdetect software recognized and counts the cone cells automatically (Fig. 2F). When the software missed cones, two masked graders (T.M. and T.O.) inspected and manually corrected the images after automated detection of the cone cells (Fig. 2E). According to a previous report, we analyzed only cones agreed upon by both graders and we excluded cones that were not agreed upon by both graders.²⁴ After the corrections, a Voronoi diagram was displayed, and cone density was calculated automatically. From each Voronoi diagram, we calculated the HPi,¹⁷ which represents the fractional increase in four-sided and eight-sided nonhexagonal Voronoi tiles with regard to six-sided cells; a lower HPi indicates a greater deviation from the normal cone packing arrangement. All AO images were analyzed after correction for the axial length of the eye.²⁵ 

Figure 2.

 

Analysis of cone mosaic image using AO. (A, B) A cone montage image and a microperimetry image are superimposed on an infrared fundus image, respectively. (C) The images in A and B are superimposed. (D) The evaluation points on a cone mosaic image. The yellow square is the sampling window for the analysis in E. The inner green circle indicates the foveal area, and the surrounding blue circle, which has an outer diameter of 2.5 mm and an inner diameter of 1 mm, indicates the parafoveal area in OCTA analysis. (E) The AOdetect software is used to analyze the cone packing arrangement. A sampling window of 0.35 × 0.35 degrees corresponds to the microperimetry spot size at the point of the microperimetry measurement. The AOdetect software recognizes and counts the cone cells automatically. Crosses indicate automatically counted cones; yellow dots, additional cones; green dashed circles; excluded cells. (F) A Voronoi diagram of the cone arrangement in E. The number of the closest neighbors (nn) to the cone cells (%) is determined and then the HPi is calculated using the following formula: 6-sided cell - (4-sided cell + 8-sided cell)/100.

OCTA Image Processing and Analysis

We acquired 3 × 3-mm images centered on the fovea using the RTVue-XR Avanti system (Optovue Inc., Fremont, CA, USA) with split-spectrum amplitude-decorrelation angiography software.²⁶ The foveal avascular zone (FAZ) parameters (i.e. FAZ area and FAZ perimeter), were calculated automatically by the built-in software (AngioVue version 2018.0.0.18). To analyze the choriocapillaris microcirculation, two sets of high-quality OCTA images were obtained from each eye (image quality, >7/10). After image acquisition, we excluded eyes with movement or shadow artifacts on the OCTA images. The retinal pigment epithelium (RPE) layer was segmented automatically by the built-in segmentation tool in the review software. The choriocapillaris slab then was segmented with a thickness of 10 µm under 31 µm from the RPE layer.¹⁴ Thereafter, the following processing was performed using the ImageJ software. Two choriocapillaris images generated from two different OCTA cube scan sets were stacked to create a two-frame image and were registered before multiple image averaging.²⁷ First, the two-frame stack of superficial slabs and choriocapillaris slabs was registered, according to previous reports.^28,29 After registration, the two-frame stacks of the choriocapillaris and superficial slabs were compounded into one image by projecting the average intensity (Fig. 3A), and a central square area of 284 × 284 pixels was cropped. The averaged choriocapillaris slabs then were binarized for quantitative image analysis of the FD using the Phansalkar method (radius, 15 pixels).¹⁴ To exclude shadows from confounding the analysis, a thresholded image of the large superficial retinal vessels was merged with the binarized choriocapillaris image (Fig. 3A) according to previous reports.^14,30 After binarization, a 2.5-mm-diameter circle centered on each microperimetry measurement point was cropped. A 1-mm-diameter circle image was used as the foveal area, and a ring-shaped area excluding the 1-mm central circle was used for the parafoveal area, similar to a previous study.¹⁵ The terms “foveal” and “parafoveal” were defined, respectively, by the area of the center and the area of the inner circle in the Early Treatment Diabetic Retinopathy Study (ETDRS) grid. The parafoveal area of the ETDRS grid generally is donut-shaped with an outer diameter of 3 mm and an inner diameter of 1 mm. However, the outer diameter of 2.5 mm was used in the current study to prevent the boundary edge of the OCTA images from being included in the selected area of analysis. The ratio of black pixels in the area, except for the superficial blood vessels, was calculated as the FD (%) (Fig. 4). All OCTA parameters were analyzed after correction for the axial length of the eye.³¹ 

Figure 3.

 

Analysis of FD in the choriocapillaris (CC) using OCTA. (A) Two superficial and choriocapillaris slabs are registered and averaged, and a central square area of 284 × 284 pixels is cropped. The averaged choriocapillaris slab is binarized using the Phansalkar method and merged with a masked image of the superficial vessels. (B) The averaged superficial and binarized CC slabs are superimposed on an infrared (IR) fundus image. (C) Two merged images in B are superimposed. A microperimetry (MP) central measurement point is marked on the superimposed image (yellow dot). Images other than the binarized choriocapillaris slab and the central point (yellow dot) then are deleted.

Figure 4.

 

Foveal and parafoveal areas of FDs in a binarized choriocapillaris slab. (A) Circles 2.5 mm in diameter (blue circle) and 1 mm in diameter (green circle) are selected on a binarized choriocapillaris image and the central microperimetry point (yellow dot) is marked. (B) A 2.5-mm circle is cropped from the entire image and a 1-mm circle then is cropped from the 2.5-mm circle. The residual ring-shaped image is the parafoveal area. (C) The cropped 1-mm circle is the foveal area.

Image Overlay

Using the retinal vasculature as a guide, the AO montage image, a microperimetry image, and OCTA images then were overlaid manually over the infrared fundus images to facilitate correlations among the RS, HPi, and FD, respectively (Figs. 2, 3B, 3C). Twelve cone mosaic and microperimetry measurement points were arranged within each 1-mm-diameter central area and the 2.5-mm donut-shaped area. The averaged values of the HPi and RS within these areas were defined similarly to the OCTA analysis areas of FD, that is, as foveal and parafoveal (Fig. 2D). These processes were performed by GNU Image Manipulation Program (version 2.10.10, GIMP Development Team). 

Statistical Analysis

To obtain the foveal and parafoveal RS and HPi values, we calculated the average of 12 points within the foveal and parafoveal areas. The categorical values were compared using Fisher's exact test. One-way analysis of variance was performed to compare the FD, RS, and HPi among the eyes with NDR, NPDR, and PDR at the foveal and parafoveal areas using the Holm post hoc test. The relationships between the FD and HPi and between the FD and RS were determined with Pearson correlation coefficients at the foveal and parafoveal areas. Statistical analyses were performed using EZR³² (Jichi Medical University Saitama Medical Center), a graphic user interface for R (The R Foundation for Statistical Computing, version 3.6.0), and a modified version of the R commander (version 2.5-3). All measurements are expressed as the means ± standard errors. All P values were based on two-sided tests; P < 0.05 was considered statistically significant. 

Results

Thirteen healthy participants and 26 patients with diabetes were recruited for this study, 4 patients with NDR, 12 with NPDR, and 10 with PDR. Patients with NPDR were significantly (P = 0.033) older than the PDR group, and a significantly (P = 0.034) increased number of patients had a history of PRP to treat advanced DR. Other parameters, including the FAZ area and FAZ perimeter in the inner retina, did not differ significantly. The cone density in the foveal area was significantly (P = 0.044) decreased in patients with PDR compared to patients with NDR (P = 0.043; Table). 

Table.

 

Patient Characteristics

 

At both the foveal and parafoveal areas, the FD was significantly (P = 0.003 and P = 0.006, respectively) higher in eyes with PDR than in the control eyes (P = 0.001 and P = 0.005, respectively), and the FD increased numerically with advanced DR, although the difference did not reach significance. The parafoveal FD was significantly (P < 0.001) higher than the foveal FD only in patients with NPDR (Fig. 5A). Overall, the RS decreased with advanced disease stage. The foveal and parafoveal RS values in NPDR and PDR were significantly lower compared with the control and NDR values. The foveal RS was significantly (P = 0.033) lower than the parafoveal RS in patients with PDR. The HPi, especially at the parafoveal area, was decreased significantly in those with DR (i.e. NPDR and PDR) compared to those with NDR (P = 0.012 and P = 0.037, respectively), but there was no difference between NPDR and PDR. There also were no differences in the foveal HPi values among all groups. The foveal HPi was significantly (P = 0.049) lower than the parafoveal HPi in patients with NDR (Fig. 6). 

Figure 5.

 

Each parameter of the foveal and parafoveal areas in each group of control and patients with DR. (A) The FD values. (B) The RS values. Statistical analyses between the areas and among the patients with DR were performed using the paired t-test and one-way analysis of variance with post hoc (Holm test), respectively. The mean value ± SD is indicated for each group. *P < 0.05, **P < 0.01. dB, decibels.

Figure 6.

 

The HPi values of the foveal and parafoveal areas in each group of patients with DR. Statistical analyses between the areas and among the patients with DR were performed using the paired t-test and one-way analysis of variance with post hoc (Holm test), respectively. The mean value ± SD is indicated for each group. *P < 0.05, **P < 0.01. dB, decibels.

No association was seen between the RS and FD in NDR; however, a significant correlation was found between the FD and RS in patients with NPDR and PDR both at the foveal and parafoveal areas (fovea, r = −0.58; P = 0.046 and r = −0.82; P = 0.003; parafovea, r = −0.59; P = 0.044 and r = −0.72; and P = 0.019; Figs. 7A–F). There was no significant correlation between the FD and RS in the controls (foveal, r = 0.05; P = 0.87, parafoveal, r = −0.22; P = 0.47). The OCTA parameters for retinal circulation, including the FAZ area and FAZ perimeter in the inner retina, did not differ. In contrast, the association between the HPi and FD was poor. In both the foveal and parafoveal areas, there was no significant correlation between the FD and HPi in any group (Figs. 8A–F). 

Figure 7.

 

The FD and RS values in patients with diabetes. There is no significant correlation (A, D) in NDR, but there are significant differences (B, E) in NPDR (fovea; r = −0.58, P = 0.046; parafovea; r = −0.59, P = 0.044) and (C, F) in PDR (fovea; r = −0.82, P = 0.003; parafovea; r = −0.72, P = 0.019). dB, decibels.

Figure 8.

 

The association between the FD and HPi in patients with diabetes. In the foveal area, there is no significant correlation (A–C) in all groups. In the parafoveal area, the FD is not correlated significantly with the HPi (D–F) across all groups.

Discussion

In the current study, the foveal and parafoveal FD were increased significantly in patients with PDR compared with the control participants. The current study corroborates the results of a previous study that reported that the capillary loss in the choriocapillaris tended to increase in advanced DR.¹¹ However, in contrast to a previous study,⁵ the parafoveal choriocapillaris microcirculation did not differ significantly across all diabetes mellitus groups, suggesting that the parafoveal choriocapillaris microcirculation may not depend on DR progression. 

Regarding the FD, there was no significant difference between NDR and NPDR in the foveal area. Analysis of the choriocapillaris FD may be affected by several factors. In addition to DR severity, DR treatment may affect the choriocapillaris. Anti-VEGF therapy can decrease the choriocapillaris capillary density.³³ In contrast, PRP reportedly increases the choriocapillaris microcirculation at the foveal area in patients with severe NPDR and PDR, as shown by laser Doppler flowmetry.³⁴ Because some current patients with NPDR had a medical history of PRP and injection of anti-VEGF therapy, the most plausible explanation for the absence of a difference in the FD between NDR and NPDR is that the FD in these patients was confounded by treatments. 

The current study identified significant differences in the RS among control and all DR groups in the foveal and parafoveal areas, and the RS tended to decrease in advanced DR. Consistent with our results, several studies in which the MP-1 microperimeter (Nidek, Aichi, Japan) was used reported that the foveal and parafoveal RS values were associated with progression of DR.^35,36 In addition, compared with the parafoveal RS, the foveal RS was reduced in patients with PDR in the current study. A previous study reported that the severity of foveal ischemia was related to outer retinal damage³⁷; notably eyes with PDR are prone to develop foveal ischemia.³⁸ 

Although the foveal HPi did not differ significantly across all groups, the parafoveal HPi in patients with NPDR and PDR was significantly lower than that of NDR, suggesting that the parafoveal but not the foveal cone regularity may decrease as the DR developments. In patients with NDR, a previous study found that the thickness of the foveal photoreceptor layer decreased significantly compared with the control.³⁹ Regarding a regional difference, a previous report showed that the foveal cone regularity was higher than the parafoveal cone regularity in normal subjects.⁴⁰ In contrast, the current study identified a significant difference in the cone regularity between the foveal and parafoveal areas in patients with NDR; the foveal cone mosaic may be more irregular in patients with NDR compared to the parafoveal cone mosaic. Moreover, a previous study using AO-scanning laser ophthalmoscopy indicated that parafoveal cone regularity decreased in DR,⁴¹ which was consistent with our result. Therefore, early changes in the cone mosaic in diabetic eyes might be detected at the foveal area. However, regarding the cone measurements, a previous study reported that cone reflectivity fluctuates and there are cones with low reflectivity even in normal eyes.⁴² In the current study, the automated analysis software often missed those cones with low reflectivity. However, those missed cones were unclear and difficult to correct completely manually from the area of low reflectivity in addition to the automatically selected cones. Therefore, we corrected only “the clearly detected cells” by two graders. Those might affect the cone measurements and lead to underestimation of the cone density. Therefore, the HPi values might be affected by this limitation of the current study. Regarding the region at 3.5 degrees, there are many rods surrounding the cone,⁴³ and a cluster of neighboring rods can blur together and look like a single cone, because the rtx1 only has a resolution of up to 3.5 µm possibly causing overestimations of the numbers of cones in the periphery. The cone density and HPi values may have been affected by these limitations. 

The current study found a significant correlation between the choriocapillaris microcirculation and RS, irrespective of the area evaluated in the NPDR and PDR groups, but not in the control and NDR groups, suggesting that visual dysfunction is associated with choriocapillaris ischemia in eyes with DR. Experimentally, rat photoreceptors are specifically vulnerable to hypoxia.⁴⁴ A clinical study also indicated that loss of the choriocapillaris was related to the functional visual impairment in DR.¹¹ In conditions in which retinal homeostasis is impaired by development and progression of DR,^45,46 photoreceptor function may be highly vulnerable to changes in the choroidal circulation in eyes with NPDR and PDR. We speculated that choriocapillaris ischemia contributes to visual dysfunction through photoreceptor malfunction in eyes with DR. We suspect that the RS measured with the current protocol may not be sufficiently sensitive to evaluate minute changes in eyes with NDR, because there was almost no difference in the RS between the control and NDR groups. 

No correlation was found between the choriocapillaris microcirculation and cone regularity in both areas of NPDR and PDR despite an association between the choriocapillaris microcirculation and RS, suggesting that the effect of choriocapillaris flow may differ between the RS and cone mosaic. The structural damage to the cone photoreceptors in eyes with DR persists over time⁴⁷; cumulative damage resulting from a long-term abnormality of the choroidal microcirculation may result in irreversible changes in the cone mosaic in DR. However, the RS is dynamic and reversible (e.g. it improves after anti-VEGF treatment of macular edema).^48,49 Therefore, we believe it is rational to assume that the changes in the cone mosaic may not necessarily depend on the choriocapillaris microcirculation at the time of the examination. 

This study had several limitations. The study sample size was relatively small and the study design was cross-sectional. We were unable to determine causality. The HPi might be affected by the resolution of the rtx1 instrument. In addition, patients with NPDR were older than in the PDR group. The age difference may have affected our results, although a previous report did not show an association between reduced choriocapillaris flow and age in patients with DR.⁵ Many patients had undergone previous treatments either of anti-VEGF therapy or PRP. Larger-scale studies matching the effects of confounding factors affecting the chorioretinal blood flow, such as drugs⁵⁰ and systemic parameters,⁵¹ are needed to confirm our findings. Even in NDR eyes, the range of FD varied, and some eyes had more than 40% both in the foveal and parafoveal areas. It may be helpful to follow the eyes with NDR and NPDR in a future longitudinal study to observe if the FD in the choriocapillaris predicts the worsening or progression of DR. 

In summary, the current study identified an association between changes in the choriocapillaris microcirculation and the RS and cone photoreceptor mosaic in patients with diabetes. We found that changes in the choriocapillaris microcirculation are associated with visual function in eyes with DR. We believe these findings provide valuable insight into the pathophysiology of functional impairment and structural damage in patients with diabetes. 

Acknowledgments

The authors thank Nicolas Chateau, PhD (Imagine Eyes, Orsay, France) for providing the AOdetect software. 

Disclosure: T. Ro-Mase, None; S. Ishiko, None; T. Omae, None; A. Ishibazawa, None; A. Shimouchi, None; A. Yoshida, None