Introduction

Classic funduscopic β-parapapillary atrophy (β-PPA) is a region of visible large choroidal vessels and sclera due to the lack of retinal pigment epithelium (RPE). β-PPA can be histologically subclassified based on the presence or absence of Bruch’s membrane (BM): β-PPA with BM (PPA+BM; β zone) and β-PPA without BM (PPA−BM; γ zone). It has been known that PPA+BM might be due to age-related atrophy of the RPE and associated with the development of glaucoma1,2,3,4. Whereas, PPA−BM has been considered as a result of mechanical stretching in the peripapillary sclera and border tissue during axial elongation; hence a significant association between axial length and PPA−BM width has been consistently reported1,2,3,5. However, we previously reported that PPA+BM was also found in 20–30-year-old myopia subjects and suggested that PPA+BM in young myopia patients might be related to axial elongation, not age-related atrophy6. In accordance with our report, Lee et al.7 recently found an enlargement of PPA+BM as myopia progressed and proposed that PPA+BM found in children might be an entirely different entity from PPA+BM in elderly individuals, sharing an etiology with PPA−BM.

Several factors might contribute to β-PPA development during axial elongation, but the etiology of PPA+BM and PPA−BM enlargement remains elusive so far. It is unclear why some myopia patients exhibit only PPA−BM, whereas others with a similar degree of myopia exhibit both PPA+BM and PPA−BM. In this study, we hypothesized that morphologic features of the optic nerve head (ONH) and the posterior pole might give clues to the mechanisms of differential developments of β-PPA in myopic eyes. Therefore, the purpose of this work was to investigate the relationship between microstructure of β-PPA and the morphologic features of ONH and posterior pole in highly myopic eyes and to determine the mechanism of differential developments of β-PPA in myopic eyes.

Methods

Subjects

Subjects were recruited from the Young Myopia Study of Chonnam National University Hospital, an ongoing cross-sectional study that commenced in January 2018. The present 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 (CNUH-2018-147). All patients provided written informed consents before enrollment in the study.

The Young Myopia Study enrolled consecutive participants who visited the general eye clinic for medical check-up and met all of the inclusion criteria and none of the exclusion criteria. The details on the study protocol is also described in our recent report8. All subjects underwent comprehensive ophthalmologic examination consisting of the measurement of best-corrected visual acuity (BCVA), intraocular pressure (IOP) by Goldmann applanation tonometry, and refractive error by automated refraction. Anterior chamber angle assessment using gonioscopy was performed on all eyes. ONH and RNFL examination using color stereoscopic disc photography and red-free RNFL fundus photography, and Swedish Interactive Threshold Algorithm standard 30-2 perimetry with a Humphrey Field Analyzer (Carl Zeiss Meditec Inc., Dublin, CA, USA) were performed. Axial length, central corneal thickness, and corneal curvature were measured using optical low-coherence reflectometry (Lenstar; Haag-Streit AG, Koeniz, Switzerland). For all subjects, a detailed medical history was recorded.

The inclusion criteria were: (1) healthy subjects 20 to 35 years old, (2) highly myopic eyes with axial length greater than 26 mm, (3) astigmatism within ± 2 diopters (D), (4) BCVA of 20/25 or better, (5) IOP ≤ 21 mmHg, (6) normal anterior chamber angles, (7) nonglaucomatous ONHs on stereoscopic photographs (8) absence of any RNFL abnormalities on red-free fundus photographs, and (9) normal visual field (VF) results in both eyes. Because myopic refractive error can be affected by lenticular changes, and aging may increase the incidence of glaucoma, we excluded subjects older than 35 years. Normal VF presentation was defined as a glaucoma hemifield test result within normal limits, as well as mean and pattern standard deviation values associated with probabilities of normality higher than 5%. Some eyes had an enlarged blind spot associated with a large area of PPA, and such eyes were also included in this study.

Subjects were excluded if they had any of the following: (1) a family history of glaucoma in a first-degree relative, (2) history of intraocular or refractive surgery, (3) pathologic myopia (patch chorioretinal atrophy, lacquer crack lesions, intrachoroidal cavitations, or choroidal neovascularization), (4) other evidence of retinal pathology, (5) opaque media, or (6) poor-quality OCT images because of irregular tear film or poor cooperation.

Eligibility was determined by two glaucoma specialists (S.W.P and M.S.S). The evaluators were blinded to all other patient and ocular data, and an eye was excluded from the study if a consensus was not reached. In cases where both a subject’s eyes met the inclusion criteria, one eye was randomly selected for the study.

Spectral-domain optical coherence tomography imaging

All participants underwent OCT imaging using SD-OCT (Heidelberg Spectralis SD-OCT; Spectralis software version 6.9.4; Heidelberg Engineering GmbH, Heidelberg, Germany). One experienced operator performed all the OCT scans. Magnification error was corrected using the formula provided by the manufacturer based on the results of keratometry and focus setting during image acquisition. OCT images with insufficient quality (typically truncated B-scans and quality score <30) were excluded.

Bruch’s membrane opening (BMO) area was measured using the Glaucoma Module Premium Edition (GMPE) software (Fig. 1A). The scan protocol composed of 24 equally spaced radial B-scans, each with 768 A-scans covering a 15° region centered on the optic disc. Twenty-five B-scans were averaged automatically for each scan location. The software automatically detected the 48 BMO points from the 24 radial scans. All the B-scans were manually checked, and BMO points were corrected when necessary. When BMO points were indiscernible, the points were fitted with a spline to derive a closed curve and a smooth contour line based on BMO points of adjacent B-scans8. If BMO points were indiscernible on four or more consecutive B-scan images, the eyes were excluded from the study. First, BMO points were checked by one experienced evaluator (M.S.S) and then reassessed and confirmed by a senior glaucoma specialist (S.W.P). In case of any discrepancy, the eyes were excluded from the analysis. BMO area was computed automatically by the built-in software. In this scan mode, the foveal pit and two BMO points in each of the two radial B-scans (that were perpendicular to each other) were automatically segmented to estimate the center of BMO and determine the FoBMO axis (the line connecting the center of the fovea and BMO center), which served as a reference for the scans. The FoBMO angle was defined as the angle between the FoBMO axis and the horizontal axis of the acquired image frame.

Figure 1
figure 1

Methods used for the ONH and CT analysis by using SD-OCT. (A) Bruch’s membrane opening (red dots) is identified on confocal scanning laser ophthalmoscopy and B-scan images. FoBMO axis (the line connecting the center of the fovea and BMO center) and FoBMO angle are also illustrated on the fundus image. Scan from EDI OCT scan image passing through the ONH (center, midsuperior, and midinferior) is also presented. The yellow arrow head indicates the RPE termination, red arrow head indicates the termination of BM, and blue arrow head indicates the optic disc margin. The distance from the edge of BM (red arrow head) to the margin of RPE termination (yellow arrow head) and to the optic disc margin (blue arrow head) are defined as PPA+BM width and PPA−BM width, respectively. Anterior lamina depth is defined as the distance between the sclerochoroid junction reference line (white solid line) and the anterior border of the LC (yellow dashed line). LC thickness is defined as the distance between the anterior and posterior borders (yellow dashed lines) of the LC. (B) Circular peripapillary scan of 4.1 mm is obtained. After the manual segmentation of BM (red solid line) and the sclerochoroidal border (green solid line), peripapillary CT was calculated. The global and the 6 Garway-Heath regional values (NS, N, NI, TI, T, and TS) was used. (C) EDI horizontal scan passing through the center of the fovea was used to measure subfoveal CT. ONH = optic nerve head; CT = choroidal thickness; BM = Bruch’s membrane; BMO = Bruch’s membrane opening; EDI = enhanced depth imaging; RPE = retinal pigment epithelium; PPA+BM = β-parapapillary atrophy with BM; PPA−BM = β-parapapillary atrophy without BM; LC = lamina cribrosa.

Optic nerve head measurements

The Heidelberg Spectralis OCT enhanced depth imaging (EDI) mode was used for the other ONH measurements (Fig. 1A). The details on the ONH measurements are also described in our previous reports6,8. The ONH was scanned by centering a 15° × 10° rectangular scan on the ONH. Each OCT volume consisted of 49 serial horizontal B-scans (4.5 mm lines; 50 images averaged) spaced at approximately 63 μm intervals. Infrared (IR) fundus images were acquired simultaneously by using a confocal scanning laser ophthalmoscope. Three sections that passed through the ONH in the center, midsuperior, and midinferior regions were selected. Mid-superior and mid-inferior locations were defined as the middle of the vertical line connecting the center to the superior and inferior margin of the optic disc. The closest horizontal B-scan frames to mid-superior and mid-inferior points were chosen for image analysis. All the study parameters were measured in each of these frames by two independent, blinded examiners (M.S.S and H.H).

Temporal β-zone PPA margin and disc margin were defined using IR fundus images. Temporal BMO points were identified in the horizontal EDI scan image. PPA width was measured using the synchronous viewing of the IR fundus image and the selected location on the OCT scan. β-PPA width was defined as the distance between the beginning of the retinal pigment epithelium (RPE) (i.e., temporal β-PPA margin) and temporal disc margin on each horizontal B-scan image. Based on the location of BM termination, the β-PPA was further divided into PPA+BM and PPA−BM. PPA+BM width was defined as the distance from the beginning of the RPE to BMO, and PPA−BM width was defined as the distance from the temporal disc margin to the beginning of BM. In case of the termination of BM is not clearly visible on the horizontal EDI B-scan, we reviewed the radial B-scans from GMPE mode to identify the result of BMO delineation, and referred to them for determination of BM termination points on the EDI B-scans. The measurement was performed using a previously described method6,7,8,9,10,11.

The sclerochoroidal junction reference plane was used to overcome the effect of choroidal thickness on lamina cribrosa (LC) depth measurements12,13. The sclerochoroidal junction reference line was defined as the line connecting two points of the anterior scleral surface located at 1750 μm from the center of BMO in each B-scan image. The vertical distance between the reference line and the anterior LC surface was measured at the center of the ONH and defined as anterior lamina depth. The anterior and posterior margins of the highly reflective region at the ONH vertical center in each B-scan image were used as the borders of the LC, and the perpendicular distance between these two borders was defined as LC thickness. In cases where the central retinal vessel trunk prevented visualization, measurements were performed on the temporal side.

The measurement was performed using a built-in caliper tool of the intrinsic OCT viewer, and average data of three horizontal B-scan images (center, midsuperior, and midinferior) were calculated and used in this study. The final analysis used the mean of the values obtained by the two examiners.

Detection of lamina cribrosa defect

The EDI OCT images of the ONH were reviewed carefully for detecting focal LC defects. A focal LC defect was defined as a loss of high reflectivity from the anterior-to-posterior border of the full-thickness LC on B-scan images. To avoid false-positives, the defects were required to have had a maximal diameter greater than 100 μm and a depth greater than 30 μm, and to have been present in two adjacent B-scans14,15. Shadows were differentiated from LC defects on the basis of their characteristic signal void behind the vessels and tissues. LC defect margins were independent of the location of the vessels and neural tissues. Images were reviewed by 2 experienced glaucoma specialists (M.S.S and S.W.P), and disagreements were addressed via discussion between the 2 evaluators to achieve consensus. An eye was excluded from study analyses if a consensus could not be reached.

Measurement of macular bruch’s membrane length

Using the scan running through the fovea and the center of the optic disc (from the posterior pole asymmetry analysis protocol of the Heidelberg Spectralis OCT), the macular BM lengths were measured16. Macular BM length was defined as the distance between the fovea to the end of BM, in the direction of the ONH17. The averaged data from two independent examiners (M.S.S and H.H) were used in this study.

Choroid thickness measurement

We obtained 360° circular RNFL measurements centered on BMO center. Among the three circular scans along the peripapillary circles (diameters of 3.5, 4.1, and 4.7 mm), the 4.1 mm diameter scans were analyzed to minimize the interference of a large PPA on the OCT scan path. For the measurement of the peripapillary choroidal thickness (CT), the upper and lower segmentation lines of the circular scan were manually delineated. The lines were adjusted to align with the inner scleral wall and posterior border of the RPE to define the outer and inner boundaries of the choroid, respectively (Fig. 1B). Peripapillary CT was automatically computed using the RNFL thickness sector algorithms (global, nasal-superior [NS], nasal [N], nasal-inferior [NI], temporal-inferior [TI], temporal [T], and temporal-superior [TS]). Subfoveal CT was measured from the horizontal EDI scan running through the fovea. The vertical distance from the hyper-scattering outer border of the RPE to the inner border of the sclera at the fovea was defined as subfoveal CT (Fig. 1C)8. Images in which the RPE and chorioscleral interface could not be clearly identified were excluded from the analysis. The averaged data from two independent examiners (M.S.S and H.H) were used in this study.

Statistical analysis

SPSS version 21.0 (SPSS, Chicago, IL, USA) was used for all statistical analyses. Agreement of PPA+BM width, PPA−BM width, LC thickness, anterior laminar depth, macular BM length, and peripapillary and subfoveal CT between two observers was assessed using the Bland-Altman method, which plots the means against differences18. The limits of agreement were defined as the mean differences of two measurements ± 1.96 standard deviations (SD) of the difference. The normality of distribution was verified using the Shapiro-Wilk normality test. Baseline characteristics were reported in counts and proportions or mean ± SD values as appropriate. Groups were compared using the chi-square test, Student’s t-test, or Mann-Whitney U test as appropriate. Linear regressions were used to search for associations between the microstructure of β-PPA and various ocular parameters. Standardized and unstandardized regression coefficients, with a 95% confidence level, were presented. When performing multivariate regression analysis, factors showing multicollinearity were not included as explanatory variables. All P values from linear regression analyses were adjusted to control the false discovery rate using the Benjamini-Hochberg procedure19. A P value of less than 0.05 was considered statistically significant.

Results

Among the eyes enrolled in the Young Myopia Study, 101 eyes of 101 subjects met the inclusion criteria with axial lengths greater than 26 mm. Of these, seven eyes were excluded because of indiscernible BMO on four or more consecutive B-scan images or discrepancy in BMO determination between the glaucoma specialists (M.S.S and S.W.P). Five eyes were additionally excluded because of inadequate visualization of the anterior LC surface, posterior LC surface, or both. None of the subjects was excluded because of the inability to determine RPE and chorioscleral interface. Finally, the remaining 89 eyes of the 89 subjects with high myopia were evaluated. Among the 89 eyes, automated detection of BMO points was found to be inaccurate in 57 eyes (64.04%), and manual corrections were made. Interobserver agreement, determined using Bland-Altman plots, in the measurements of PPA+BM width, PPA−BM width, LC thickness, macular BM length, and peripapillary and subfoveal CT for all the subjects showed no systematic differences in measurements (data not shown). Demographic and ocular characteristics are described in Table 1. The mean axial length of included eyes was 27.03 ± 0.88 mm (range, 26.04–29.76 mm).

Table 1 Demographic characteristics of the highly myopic eyes.

Table 2 shows the associations between the microstructure of β-PPA and ocular parameters after adjusting for age, sex, and axial length in highly myopic eyes. PPA−BM was positively correlated with BMO area (P < 0.001), and negatively correlated with subfoveal, global peripapillary, temporal peripapillary, and temporal-inferior peripapillary CT (P = 0.005, P = 0.042, P = 0.004, and P = 0.021, respectively). Whereas, BMO area was not associated with the width of PPA+BM. The large PPA+BM was significantly related to a thinner LC (P < 0.001), deeper anterior lamina surface (P = 0.002), larger FoBMO angle (P = 0.031), longer macular BM length (P < 0.001), and thinner subfoveal and temporal peripapillary CT (P = 0.009 and P = 0.005) in highly myopic eyes.

Table 2 Association between the microstructure of parapapillary atrophy and ocular parameters in highly myopic eyes.

In multivariate linear regression model, to explore the independent factors associated with the microstructure of β-PPA, the PPA+BM was negatively associated with LC thickness (P = 0.003) and temporal peripapillary CT (P = 0.034), and positively associated with anterior lamina depth (P < 0.001), FoBMO angle (P = 0.034), and macular BM length (P = 0.008). The overall R2 value derived from the multivariate linear regression model for PPA+BM width was 0.472 (Table 3). As for PPA−BM, the multivariate model showed that only BMO area (P = 0.001) was independently associated with the PPA−BM width in highly myopic eyes. The R2 value for the model was 0.247 (Table 4).

Table 3 Associations between enlargement of PPA+BM width and ocular parameters in highly myopic eyes.
Table 4 Associations between enlargement of PPA−BM width and ocular parameters in highly myopic eyes.

We further analyzed the eyes with β-PPA (81 eyes) after dividing them into two groups according to the presence of PPA+BM: eyes with PPA+BM group (n = 62) and eyes without PPA+BM group (with PPA−BM only; n = 19). Among 62 eyes with PPA+BM, 3 eyes had only PPA+BM, and 59 eyes had both PPA+BM and PPA−BM. The mean width of PPA+BM in 3 eyes having only PPA+BM was 199.33 ± 9.71 μm. The mean width of PPA+BM and PPA−BM in 59 eyes with both PPA+BM and PPA−BM in their β-PPA area was 346.90 ± 160.26 μm and 415.51 ± 178.91 μm, respectively. In comparison with eyes with PPA−BM only, those with PPA+BM exhibited a significantly thinner LC (P = 0.009), deeper anterior lamina surface (P = 0.029), larger FoBMO angle (P = 0.007), longer macular BM length (P = 0.010), and thinner subfoveal CT (P = 0.037) (Table 5). Figures 2 and 3 present representative cases of highly myopic eyes with a different microstructure of β-PPA.

Table 5 Comparison of characteristics between highly myopic eyes with PPA+BM and without PPA+BM (with PPA−BM only).
Figure 2
figure 2

Images from myopic eyes with similar axial length and different microstructure of β-PPA. (A) A myopic eye of 29-year-old man with axial length of 26.35 mm. The width of PPA+BM and PPA−BM were 462 μm and 367 μm, respectively. The BMO area was 2.31 mm2. The EDI OCT image demonstrating a thinner LC and deeper anterior lamina surface compared with those from B. (B) A myopic eye 25-year-old man with axial length of 26.27 mm. The width of PPA+BM and PPA−BM were 157 μm and 485 μm, respectively. The BMO area was 3.05 mm2. The EDI OCT image demonstrating thicker LC thickness and shallower anterior lamina depth compared with those from A.

Figure 3
figure 3

Representative cases showing relationship between the macular BM length and microstructure of β-PPA. (A) A myopic eye of 27-year-old-woman with PPA−BM only. The macular BM length was 3466 μm. (B) A myopic eyes of 26-year-old woman with both PPA+BM and PPA−BM. The macular BM length was 4435 μm. Although the eye of case A has longer axial length, the eye of case B has longer macular BM length. PPA+BM = β-parapapillary atrophy with BM; PPA−BM = β-parapapillary atrophy without BM.

Discussion

In this study, we evaluated 89 highly myopic eyes with axial length greater than 26 mm, and β-PPA was found in 81 of 89 eyes (91.01%). The microstructure of β-PPA was different among the 81 eyes; 19 eyes had only PPA−BM, whereas 62 eyes had PPA+BM in their β-PPA area. Importantly, we observed that anatomical parameters associated with the width of PPA+BM and PPA−BM were markedly different, after adjusting for significant confounders, age, sex, and axial length2,3. The large PPA+BM was associated with a thin LC, deep anterior lamina surface, long macular BM length, close FoBMO axis relative to the horizontal axis, and temporal peripapillary CT thinning. Meanwhile, the width of PPA−BM was positively associated with BMO area. To our knowledge, this is the first work thoroughly investigating the ONH characteristics which have a different β-PPA microstructure among highly myopic subjects. Earlier researches cannot be compared with the current study since earlier works included a wide range of age and refractive error groups2,3,20. Our conclusions are clinically important given the importance of β-PPA in myopic eyes, and furthermore, this research might provide clues to a association between PPA+BM and glaucoma susceptibility in myopic eyes.

PPA−BM, sometimes termed as the gamma zone, is thought to be the result from stretching of peripapillary border tissue during globe elongation. The strong association between axial length and PPA−BM width has been well documented1,2,3,5,20,21. However, this relationship seems to be non-linear in clinical practice, and Jonas et al.1 described the prevalence and size of PPA−BM increased steeply, starting at an axial length of 26.5 mm, in their histomorphometric investigations of human globes. Recently, Zhang et al.20,21 demonstrated a steep increase in BMO diameters beyond an axial length of 26.0 mm and reported that such BMO enlargement was significantly associated with the large PPA−BM in myopic eyes. Similarly, we found a significant positive association between the PPA−BM width and BMO area, after adjusting for significant confounders. Following the speculation by Zhang et al.20, we suggest that development and enlargement of the PPA−BM may be due to the disproportional enlargement between BMO and the peripapillary sclera opening. When the BMO enlarges during axial elongation, the peripapillary sclera opening might also enlarge, but due to the LC anchored at the peripapillary sclera flange, the peripapillary sclera opening enlargement cannot keep pace with the BMO enlargement. Regarding the reason for the location of PPA−BM (most PPA−BM is found at the temporal side of the ONH), possibly the concomitant nasal shifting of LC might cause the asymmetric stretching of nasal and temporal peripapillary border tissue, thus resulting in enlargement of the PPA−BM mostly at the temporal side of the ONH. Recently, Kim et al.7,22,23 introduced the concept of LC nasal shifting during myopia progression after the prospective longitudinal observation of the ONH in South Korean children. However, this concept is only a hypothesis and it has not been scientifically confirmed, thus future longitudinal investigations might be needed to confirm our assumption. Meanwhile, the BMO area was not associated with the width of PPA+BM, which implies that enlargement of PPA−BM and PPA+BM might be the result of different mechanisms during axial elongation.

It has been suggested that PPA+BM results from age-related degeneration of the RPE-BM complex. Curcio et al.4 previously demonstrated that the length of BM bared of RPE cells was significantly wider in the oldest eyes in their histomorphometric study. However, we have shown that PPA+BM can also be found in young myopic eyes. Our results are in line with recent reports6,7. The mechanism of RPE atrophy in the peripapillary region with aging is not fully determined. It is generally thought to be caused by insufficient blood supply due to a thinner choroid24,25. Regarding the PPA+BM found in young myopic eyes, Lee et al.7 postulated that RPE migration over the BM or the BM stretching, rather than the ischemic insults, might be associated with the PPA+BM enlargement during axial elongation. They speculated that mechanical stress on the peripapillary area during axial elongation might be involved in the enlargement of the PPA+BM in highly myopic eyes. Our result, which shows a significant positive relationship between the macular BM length and PPA+BM width, might support their speculation.

Jonas et al.17 showed that BM length in the macular lesion was not affected by axial elongation. They suggested that the axial elongation-related increase in the optic disc-fovea distance predominantly occurred through an enlargement of PPA−BM, whereas the macular BM does not expand. In line with their findings, there was no significant relationship between the axial length and macular BM length (r = −0.100, P = 0.351) in this study. Interestingly, we observed a statistically significant relationship between the macular BM length and the PPA+BM width among highly myopic eyes in multivariate analysis. Our findings might indicate that BM can be stretched during axial elongation, which may affect PPA+BM enlargement in highly myopic eyes. The close relationship between the FoBMO angle and PPA+BM width might be explained by macular BM stretching, since the FoBMO axis relative to the ONH could be affected by the degree and direction of macular BM stretching.

Interestingly, we observed that eyes with large PPA+BM represent characteristic structural changes in LC. The larger PPA+BM was significantly related to a thinner LC and deeper anterior lamina surface. Previous investigations revealed LC thinning in highly myopic eyes, and it was speculated that axial elongation-associated thinning of the LC might be due to enlargement of BMO26. However, we could not found a significant correlation between the LC thickness and BMO area (r = 0.089, P = 0.406) in our highly myopic subjects. Our finding shows that the ONH biomechanics are complex in myopic eyes, and we believe that the individual biomechanical properties of LC and adjacent peripapillary tissue might contribute to this complexity. Further investigation is needed to clarify the mechanism of LC changes found in eyes with a large PPA+BM.

It is well known that the choroid is markedly thinner in highly myopic eyes than in emmetropic eyes27,28. Fang et al.29 demonstrated that the choroidal change does not occur symmetrically and that there is a sequence for thinning. They suggested that this progressive and continuous choroidal thinning plays a crucial role in the development of various degrees of myopic maculopathy. Based on their observation, choroidal thinning starts from the temporal peripapillary region to the fovea and then enlarges toward the entire posterior pole. Considering that the choroidal thinning in high myopia is mechanical strain-dependent, such a non-uniform thinning of the choroid suggests that there might be an asymmetry in ocular growth patterns during myopia progression. Also, it indicates that temporal peripapillary and subfoveal choroid regions are more susceptible to mechanical strain than the other regions in myopic eyes. In line with their findings, in our study, both subfoveal and temporal peripapillary CT showed a tendency for thinning as the width of PPA+BM and PPA−BM became larger. However, the thinning was more remarkable in eyes with large PPA+BM. In multivariate analysis, temporal peripapillary CT was significantly associated with the width of PPA+BM, and eyes with PPA+BM exhibited significantly thinner subfoveal CT compared with eyes with PPA−BM only. Our findings suggest that eyes with large PPA+BM might have undergone more strain on the region between the fovea and the ONH during eyeball growth. Moreover, this hypothesis would be supported by macular BM lengthening in eyes with large PPA+BM.

Combining all these findings, we propose that the differential development of β-PPA in highly myopic eyes might result from different degrees of mechanical strain and their topographical variation during axial elongation. Large BMO area was associated with the large PPA−BM. Whereas large PPA+BM was associated with the longer macular BM, thinner subfoveal and peripapillary choroid, and some characteristic LC changes. If the mechanical strain is too high thus BMO widening is not enough to release this strain, the LC might be affected (thinning and deepening) and subsequently, BM lengthening and subfoveal choroid thinning might occur. These anatomical changes might be associated with the individual difference in biomechanical property of BM during axial elongation.

Limitations of our study include the modest sample size. However, by employing strict inclusion criteria regarding age and axial length, we could meticulously exclude the influence of age-related degeneration of the parapapillary area and investigate the exclusive influence of myopia. Another limitation is related to the difficulties of imaging BMO and deep ONH structures. The determination of BMO might be inaccurate in highly myopic eyes, and the anterior and posterior lamina surface could not be visualized in some healthy eyes30. However, in this study, we strictly excluded such patients having inadequate visibility of BMO or LC (11.88%). Third, the anterior laminar depth was measured at the center of the ONH. We did not use the BMO center as measurement location, because BMO center is often close to the disc edge in the tilted optic disc and, thus, considered to be not suitable for measuring anterior laminar depth. However, detection of ONH center by IR fundus image is sometimes difficult in myopic eyes. There is a possibility that the uncertain location of the ONH center might have caused inaccurate anterior laminar depth measurement. Fourth, the associations between the width of PPA+BM and LC changes could be affected by IOP, since LC is highly sensitive to IOP. However, the relationships between the width of PPA+BM and LC thickness and anterior laminar depth were remained statistically significant after adjusting for age, sex, axial length, and IOP (β = −0.402, P < 0.001 and β = 0.285, P = 0.008, respectively), thus the IOP is not likely to affect the main outcome of our study. Finally, this study is cross-sectional. Therefore, it is uncertain whether the previously described characteristics are the results of axial elongation. For example, Fang et al.29 reported a progressive choroidal thinning with aging in young, highly myopic adults (4.91 μ/year). The subfoveal and peripapillary choroidal thinning found in eyes with enlarged PPA+BM could be a result of the aging process, not a result of axial elongation induced mechanical strain. A future longitudinal study could help to clarify this study’s findings.

In conclusion, we have documented structural differences in the ONH and posterior pole in highly myopic eyes with a different β-PPA microstructure. Large PPA+BM was associated with the longer macular BM and structural changes in the LC. Future studies might be needed to elucidate the association between these structural features and development of glaucoma in myopic eyes.