Acircularity and circularity indexes of the foveal avascular zone in high myopia

This study explored the association between foveal avascular zone (FAZ) parameters and high myopia using optical coherence tomography angiography. We divided 106 eyes of 106 patients into quartiles based on the axial length. The upper quartile was then defined as the high myopia group (n = 27), while the lower quartile was the non-high myopia group (n = 26). The areas and minor axis lengths of superficial and deep FAZ, the perimeters and major axis lengths of deep FAZ were significantly larger in eyes with high myopia than in eyes with non-high myopia (P < 0.05). Inversely, the subfoveal choroidal thickness was significantly thinner in eyes with high myopia than in those with non-high myopia. Linear regression analyses showed that no significant correlation was observed between FAZ areas and acircularity and circularity indexes of FAZ in non-high myopia group. Conversely, FAZ areas strongly correlated with acircularity and circularity indexes of FAZ in high myopia group. We found that an increase in the FAZ area in highly myopic eyes was accompanied by a significant variation in FAZ acircularity and circularity indexes. Further research should address whether these findings are associated with future disease development in highly myopic eyes.

Associations of the FAZ parameters of OCTA with FAZ area. Linear regression analysis adjusted for age and sex was conducted to examine whether the FAZ parameters of OCTA were correlated with FAZ area in the non-high myopia and high myopia groups. As shown in Table 3, the areas of the superficial and deep FAZ were not significantly correlated with the circularity index of the FAZ in the non-high myopia group. Conversely, the areas of the superficial and deep FAZ were significantly negatively correlated with circularity index (P = 0.001 in both layers) of the FAZ in the high myopia group. Likewise, the areas of the superficial and deep FAZ were not significantly correlated with the acircularity index of the FAZ in the non-high myopia group (Fig. 1A, B), but the areas of the superficial and deep FAZ were significantly positively correlated with acircularity index of the FAZ in the high myopia group (P < 0.001 in both layers) (Fig. 1C, D). In both groups, the areas of the superficial and deep FAZ were significantly positively correlated with the perimeters of the FAZ (P < 0.001 in both layers). Only the areas of the superficial FAZ in the non-high myopia group were significantly correlated with the axial ratio of the FAZ (P = 0.011). In addition, the areas of the superficial FAZ in the non-high myopia group were weakly correlated with the major axis lengths of the FAZ (P = 0.024), while other major and minor axis lengths of the FAZ were strongly associated with the FAZ areas (P < 0.001 in both groups). www.nature.com/scientificreports/

Discussion
In this study, we utilized OCTA to analyze the macular microvasculature and further explore the morphological changes of FAZ. We found that several FAZ parameters were larger in highly myopic eyes than in non-highly myopic eyes, and that the differences were even more pronounced in the deep FAZ layer. Moreover, the acircularity and circularity indexes of the FAZ were significantly correlated with FAZ enlargement in highly myopic eyes, whereas no significant relationships were found between these parameters in the non-highly myopic group. To our knowledge, this was the first report to determine FAZ deformation in high myopia. Several reports have investigated the FAZ in myopic eyes 5,[14][15][16]18,19 . Sung et al. reported lower peripapillary vessel density and larger superficial and deep FAZ areas in highly myopic eyes, as elucidated using OCTA 5   www.nature.com/scientificreports/ density in the superficial and deep layers were correlated with axial elongation 15 . All these studies reported that FAZ enlargement was accompanied by a reduction of retinal vessel density in the macular region, indicating that the size of the FAZ indirectly reflects variations in retinal perfusion. Although we did not measure vessel density in the present study, our results are in accordance with those of the previous reports. The mechanism of FAZ enlargement in myopic eyes remains unclear, although some hypotheses have been offered. First, retinal vascular trunk dragging during myopia progression may decrease retinal blood flow and further result in variation in the FAZ size. Two previous studies reported that the FAZ areas were significantly associated with optic nerve head tilt, which occurs during axial elongation 5,14 . Second, in eyes with axial elongation, overall macular thinning may cause reduced oxygen consumption, resulting in decreased retinal blood flow and increased FAZ area 20,21 . In the current study, the FAZ areas in highly myopic eyes were significantly larger than those in non-highly myopic eyes, and the difference was more significant in the deep FAZ layers, suggesting that the deep capillary plexus is more susceptible to myopia-related changes, possible because the arterial blood supply differs between the superficial and deep capillary plexuses. The oxygen and nutrition demand of the superficial retina is met by the central retinal artery, while that of the deep retina is met by the choroidal vascular system 22 . In our study, the subfoveal choroidal thickness was significantly thinner in eyes with high myopia than in those without high myopia. This result is in line with the study by Al-Sheikh et al. 23 who reported that choroidal thinning is more prominent in highly myopic eyes. Thus, the FAZ enlargement might be more pronounced in the deep retinal layer.
Numerous studies have researched FAZ deformation, but none have explored the morphological changes in the FAZ in high myopia. We speculated that either or both of the following mechanisms influence FAZ morphology: First, the irregular reduction in macular blood perfusion caused by the retinal vascular trunk dragging, with consequent FAZ deformation and enlargement during axial elongation. Second, as the retinal vascular trunk was dragged nasally in highly myopic eyes, the FAZ in the macula may have been distorted towards the ONH, leading to FAZ deformation. Recent prospective observational study on progressive myopia demonstrated that the central vascular trunk in the lamina cribrosa was dragged nasally during axial elongation 24 . Therefore, depending on geometrical change in the FAZ, it may be possible to identify myopia severity and macular perfusion status. Future longitudinal study should verify this.
So far, to our knowledge, there has been no study regarding the FAZ circularity index of high myopia and non-high myopia. However, a number of papers have been reported that FAZ circularity was decreased after suffering a period of pathophysiological process in several disease status [25][26][27][28][29] . They showed the FAZ shape relatively circular initially, and evolved from circular to acircular during disease progression. Figure 2A, B. are the scan www.nature.com/scientificreports/ images from a non-highly myopic eye. Because the eye was less affected by aforementioned two mechanisms, the shape was approach to circular. In contrast, Fig. 2C, D. are the representative images of highly myopic eye, and we can see the increase in FAZ area and decrease in circularity index. We infer that FAZ enlargement in non-high myopia is attribute to individual variation and not the result of authentic FAZ enlargement, hence the significant change of circularity index was not found in this group. As previous study showed that the circularity index is proportional to the area and inversely proportional to the square of the perimeter. In DR, the FAZ has a higher amount of irregular outlines, thus increasing the perimeter for a given surface value compared with a healthy control 26 . Likewise, the high myopization may not only cause an increase of the FAZ area, but also higher amount of irregular outlines and longer perimeter, ultimately lead to circularity index decrease. Of note, we found that the FAZ areas were not significantly associated with the acircularity or circularity index in non-highly myopic eyes. However, they were positively correlated with acircularity index and negatively correlated with circularity index in highly myopic eyes, suggesting that the FAZ area becomes larger and more acircular as myopia progresses. Variations in the FAZ acircularity/circularity indexes and FAZ size may have multiple clinical significances. Choi et al. reported that decreased FAZ circularity index may represent a disruption of the parafoveal capillary network in patients with glaucoma 25 . Kwon et al. showed that a decreased FAZ circularity index was significantly associated with central visual field defect, and that the FAZ size was significantly associated with the severity of central visual field defect in glaucomatous eyes 26 . Kim et al. also showed that decreased FAZ circularity is a good indicator of vascular dropout, and that it is associated with disease progression in vascular maculopathy 27 . Another study found a significant difference in acircularity index between controls and non-proliferative DR. This result could inform objective staging of the disease 13 , before the disease has progressed significantly, the retinal microvascular system in myopic eyes has already suffered structural damage, and discernible microvascular changes are visible in the FAZ. Given that FAZ acircularity is associated with visual field defects and macular vascular dropout, our results might explain the reason for the high prevalence of central scotoma in myopic glaucoma 30 .
Unlike other indicators such as area, perimeter, and length, the acircularity index can quantify FAZ geometry without the need for axial length measurements to correct for retinal magnification 13 . It follows that the www.nature.com/scientificreports/ acircularity index (including circularity index) could serve as indicators to predict pathological myopic changes in advance. Long-term monitoring of FAZ evolution using the acircularity and circularity indexes may be necessary to substantiate the clinical value of these metrics, which are not influenced by ocular magnification correction.
The limitations of the current study were as follows: The study was retrospective, and to evaluate the specific influence of axial myopia, we set further criteria, which resulted in a small sample size. A prospective big data study is needed to improve data accuracy. Furthermore, we did not measure the retinal thickness in detail. One previous study demonstrated that early foveal microcirculatory alterations in diabetic eyes were related to macular ganglion cell/inner plexiform layer thickness, regardless of the presence of DR 27 . Intraretinal thickness may therefore be related to FAZ deformation in myopic eyes. Future studies should investigate this relationship further to deepen the current understanding of myopia progression.
In conclusion, we observed FAZ enlargement in highly myopic eyes, especially in the deep retinal layer. This enlargement accompanied changes in the acircularity and circularity indexes of the FAZ in highly myopic eyes. Further research should address whether these findings are associated with future disease development in highly myopic eyes.

Methods
Subjects. This retrospective, cross-sectional study was conducted according to the tenets of the Declaration of Helsinki and was approved by the Institutional Review Board of Chonnam National University Hospital. The participants were informed of the study objectives and written informed consent was obtained from all participants.
Healthy myopic volunteers were recruited consecutively from July 1, 2016 to September 31, 2016. A detailed medical history was recorded for each subject. All subjects underwent complete ophthalmic examination, including best-corrected visual acuity (BCVA), intraocular pressure (IOP) using Goldmann applanation tonometry, slit-lamp examination, refractive error assessment, fundus examination using color stereoscopic disc photography and red-free RNFL fundus photography, and Swedish Interactive Threshold Algorithm standard 30-2 perimetry using a Humphrey Field Analyzer (Carl Zeiss Meditec Inc., Dublin, California, USA). All IOP measurements were recorded between 4:00 PM and 7:00 PM. Axial length and central corneal thickness were measured using optical low-coherence reflectometry (Lenstar; Haag-Streit AG, Koeniz, Switzerland).
All volunteers were required to meet the following inclusion criteria: (1) age between 18 and 40 years; (2) astigmatism within ± 2 diopters; (3) BCVA of 20/25 or better; (4) IOP ≤ 21 mmHg; (5) normal optic nerve head with non-glaucomatous changes on stereoscopic photographs, including an intact neuroretinal rim without thinning, pit, or localized pallor; (6) absence of RNFL abnormalities on red-free fundus photographs; and (7) normal visual field results in both eyes. The exclusion criteria included (1) history of glaucoma or intraocular surgery; (2) IOP > 21 mmHg; (3) myopic degeneration, including peripheral retinal tears, intrachoroidal cavitations, or choroidal neovascularization; (4) history of systemic disease that could affect the eyes, such as diabetes mellitus, (5) other evidence of retinal pathology or opaque media. For cases in which both eyes of a subject met the inclusion criteria, one eye was chosen randomly.

Spectral-domain optical coherence tomography. All subjects underwent macular imaging with
Avanti spectral-domain OCT (RTVue-XR Avanti; Optovue, Fremont, California, USA). The fovea was defined as the area within the central 1-mm ring and the parafoveal area within the central 3-mm ring of the Early Treatment Diabetic Retinopathy Study grid. The distance between the internal limiting membrane (ILM) and the middle of the retinal pigment epithelium (RPE) as the full retinal layer thickness. Subfoveal choroidal thickness was evaluated by the enhanced depth imaging (EDI) technique of the Avanti spectral-domain OCT device. From the 12-mm-length EDI scan running through the fovea, vertical distance from the hyper-scattering outer border of the RPE to the inner border of the sclera at the center of the fovea was measured and it was defined as subfoveal choroidal thickness.
Spectral-domain optical coherence tomography angiography. The spectral-domain OCT instrument used was the RTVue Avanti with prototype AngioVue® OCTA software (Optovue, Inc., Fremont, CA). This instrument operates at ~ 840 nm wavelength and conducts 70,000 A-scans per second to acquire OCTA volumes consisting of two repeated B-scans from 304 sequential, uniformly spaced locations. Each B-scan consists of 304 A-scans, giving a total of 2 × 304 × 304 A-scans per acquisition, with a total acquisition time of approximately 3 s and an axial optical resolution of ~ 5 μm. Macular 3 × 3-mm region scans centered on the fovea were captured from each subject. OCTA perfusion maps were generated using the split-spectrum amplitude decorrelation angiography (SSADA) algorithm, as previously described 31 . The Optovue software built into the Avanti OCT automatically segmented the inner retinal blood vessels into the superficial and deep capillary plexus layers. The boundaries of the superficial capillary plexus extended from 3 mm below the ILM to 15 mm below the inner plexiform layer. The deep capillary plexus extended from 15 to 70 mm below the IPL.

FAZ measurements.
To investigate the morphological changes in the FAZ, the following measurements of FAZ superficial and deep layers were made, as described in detail elsewhere: FAZ area, FAZ perimeter, acircularity index, circularity index, axial ratio, major axis length, and minor axis length 13,32 . In the superficial layer, the area and perimeter length of relatively regular FAZ could be measured automatically with ImageJ software 1.52a (National Institutes of Health, USA), then we manually measured the major axis length and minor axis length after outlining FAZ manually using ImageJ software (Fig. 3). In the deep layer, we firstly outlined the irregular FAZ manually, then measured the area and perimeter length automatically and the major axis length and minor axis length manually using ImageJ software (Fig. 4) www.nature.com/scientificreports/ would also measure it according to the method of measuring deep layer. Each OCTA image contains 605 × 605 pixel TIF was used for analysis with ImageJ software. Littman and the modified Bennett formulas were used to correct ocular magnification induced by axial variation on the area of FAZ 33 . A perfectly circular FAZ has an acircularity index equal to 1.0, with deviations from a circular shape leading to an increase in this metric. Acircularity is defined by the following equation: Circularity is a shape descriptor that can mathematically indicate the degree of similarity of the FAZ to a perfect circle. A value of 1.0 designates a perfect circle; as the circularity value decreases, the shape is increasingly less circular. Circularity is defined by the following equation: The following parameters were determined from a best-fit ellipse of the FAZ: the length of the major and minor axes, and the axial ratio, which is calculated using the following equation: To decrease the measurement error, the average data from two independent examiners (H.P. and H.Z.) were used statistically in this study.
Statistical analysis. IBM SPSS Statistics for Windows, Version 25.0 (IBM Corp., Armonk, NY, USA) was used for statistical analysis. The inter-observer agreement of the FAZ parameters was evaluated by calculating the ICC. The Mann-Whitney U test was used to compare the FAZ parameters between non-high myopia and Acircularity = perimeter of FAZ/perimeter of the circle with equal area Circularity = 4π × area / perimeter 2 Axial ratio = major axis length of FAZ/minor axis length of FAZ www.nature.com/scientificreports/ high myopia groups. Linear regression analysis was used to investigate whether the areas of the superficial and deep FAZ were correlated with other FAZ parameters in each group. P-values < 0.05 were considered statistically significant.

Data availability
The datasets generated during the current study are available from the corresponding author on reasonable request. www.nature.com/scientificreports/ Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http:// creat iveco mmons. org/ licen ses/ by/4. 0/.