Quantifying biomarkers of axonal degeneration in early glaucoma to find the disc at risk

Jul 28, 2023

Scientific Reports volume 12, Article number: 9366 (2022) Cite this article

1110 Accesses

1 Citations

1 Altmetric

Metrics details

To evaluate regional axonal-related parameters as a function of disease stage in primary open angle glaucoma (POAG) and visual field (VF) sensitivity. Spectral domain optical coherence tomography was used to acquire 20° scans of POAG (n = 117) or healthy control (n = 52) human optic nerve heads (ONHs). Region specific and mean nerve fibre layer (NFL) thicknesses, border NFL and peripapillary NFL, minimum rim width (MRW)/ area (MRA) and prelamina thickness; and volume were compared across POAG disease stages and with visual field sensitivity. Differences identified between early glaucoma (EG), preperimetric glaucoma (PG) and control (C) ONHs included thinner PG prelamina regions than in controls (p < 0.05). Mean border NFL was thinner in EG (p < 0.001) and PG (p = 0.049) compared to control eyes; and EG mean, and inferior and ST, border NFL was thinner than in PG (p < 0.01). Mean, superior and inferior PG peripapillary NFL were thinner than in controls (p < 0.05), and EG ST peripapillary NFL was thinner than in PG (p = 0.023). MRW differences included: PG SN and inferior less than in controls (p < 0.05); thinner EG mean regional, inferior, nasal, and ST MRW versus PG MRW (p < 0.05). Regional border NFL, peripapillary NFL, MRW, MRA, prelamina thickness (except centre, p = 0.127) and prelamina volume (p < 0.05) were significantly associated with VF mean deviation (MD). Novel axon-derived indices hold potential as biomarkers to detect early glaucoma and identify ONHs at risk.

Glaucoma remains the leading cause of irreversible blindness worldwide1,2 and it is estimated that over 70 million people are affected by the disease3,4,5. Intraocular pressure (IOP) is considered the major risk factor for primary open angle glaucoma (POAG)6,7,8, and current POAG treatment methods are aimed at lowering IOP via medical or surgical techniques9,10,11.

Irreversible vision loss in glaucoma is caused by the degeneration and/or loss of retinal ganglion cells (RGCs) and their axons12,13. RGC neuronal cell bodies are located within the ganglion cell layer, and their axons form the innermost retinal layer, the retinal nerve fibre layer (RNFL)14,15. Signs of glaucomatous optic neuropathy include characteristic changes in the appearance of the optic disc including enlargement of the optic cup and loss of the neuroretinal rim16,17,18,19 and RNFL defects due to loss of RGC axons20,21,22.

It is estimated that between 25 to 35% of RGCs and their axons are lost before irreversible visual field (VF) defects are detected using standard visual field tests23,24. Examination of the ONH and RNFL currently form the basis of glaucoma detection14,15,21. However, inconsistencies exist between clinicians in the subjective evaluation of the ONH to identify glaucomatous structural damage or disease progression25,26,27,28. It is likely that subtle glaucomatous disc changes are difficult to detect, identifying a need for novel optic disc assessment protocols in glaucoma. Since the introduction of optical coherence tomography (OCT)29,30 to the clinic setting, this technology has been used to quantify various RGC axon and disc related parameters to assess in vivo changes in glaucoma31,32,33,34.

RNFL measures have been reported to better detect progressive glaucoma disease, with inferior RNFL a better discriminator between healthy eyes and eyes with glaucomatous visual field loss than macula retinal thickness measures35,36. Furthermore, Lisboa et al.37 reported that RNFL measurements quantified by SD-OCT performed better than ONH and macula measures for detecting preperimetric glaucomatous damage38. Additionally, NFL thickness has been reported to provide a stronger correlation with visual field loss than optic disc cupping and neuroretinal rim area39. Jonas et al.40 reported sectorial neuroretinal rim loss, with pronounced inferior-temporal rim loss in eyes with modest glaucomatous damage (defined according to visual field Mean Deviation [MD]) and decreased superior-temporal rim width in eyes with moderate glaucoma. These studies support the use of regional RNFL and neuroretinal rim thickness measures as important clinical factors in glaucoma diagnosis and follow-up.

However, the use of RNFL thickness in glaucoma assessment does not provide sufficient information to stage glaucoma-related changes within the ONH41. Characterisation of changes to ONH structure such as prelamina cupping and lamina cribrosa alterations has been reviewed and proposed to enable more precise detection of early axonal loss in glaucoma42,43. Furthermore, peripapillary RNFL thickness measures in glaucoma disease do not take into consideration artefactual changes in RNFL that can arise owing to the angular relationship between RGC axons and the Bruch’s membrane opening (BMO) plane. To address this, Povazay et al.44 proposed neuroretinal rim measurements such as OCT-based minimum distance mapping to provide an optical correlate of true RNFL thickness around the ONH region.

Later, minimum rim width (MRW) and area (MRA) were reported to enhance glaucoma detection45,46, and correlate with RNFL and visual field loss47. However, Amini et al.48 reported that peripapillary NFL in glaucoma provided a stronger correlation with visual field sensitivity than MRW measures.

The optimal axonal parameter representative of RGC axon integrity to characterise glaucoma disease stage or be used in early detection of glaucomatous optic neuropathy remains undecided. Using high resolution OCT, this study aims to use novel indices of RGC axon integrity, including border NFL thickness (height above Bruch’s membrane termination) and prelamina volume; as well as regional analyses of other known indices; peripapillary NFL, MRW, MRA, and prelamina thickness, to characterise glaucomatous optic neuropathy as a function of POAG disease stage in attempt to determine an optimal index, or group index, for early detection of optic discs at risk of progressive damage.

This cross-sectional study included 59 POAG participants (mean age ± SD: 70.19 ± 8.43 years) and 26 control participants (mean age ± SD: 67.27 ± 5.06 years). Participants were recruited from the University Hospital of Wales, Cardiff and/or Cardiff University School of Optometry and Vision Sciences (CUS Optom). This study was approved by Wales Research Ethics Committee and CUS Optom ethics audit committee; investigations were carried out in accordance with the tenets and declaration of Helsinki. Participant demographics are presented in Table 1.

All participants underwent assessment of both eyes including best-corrected visual acuity, Goldmann applanation tonometry, slit-lamp biomicroscopy, and dilated stereoscopic examination of the ONH and fundus. Axial eye length was determined using an IOL-Master (Carl Zeiss, Meditec Inc, Germany). Central corneal thickness (CCT) and anterior chamber depth were measured with a Pentacam (Oculus Optikgeräte, Germany). Visual field status was determined using Humphrey Visual Field Analyser 24–2 SITA standard (Carl Zeiss Meditec Inc, Germany).

Participant inclusion criteria included a mean spherical correction within ± 6.00 dioptres, and less than 3.00 cylindrical correction. Control participants were required to have best-corrected visual acuity of 0.2 logMAR or better, intraocular pressure (IOP) ≤ 21 mmHg, normal visual fields, a normal appearance to the optic disc, and no ocular or systemic conditions affecting the ONH. Eyes with significant cataract or ocular conditions other than glaucoma that could affect the ONH structure or visual fields were excluded. Diagnosis of POAG by a Consultant Ophthalmologist was based upon typical glaucomatous optic disc features including diffuse or focal thinning/notching of the NRR16,49, IOP > 21 mmHg prior to topical treatment, and characteristic visual field defects50,51. POAG participants that displayed characteristic ON changes without visual field defects were classified as preperimetric glaucoma (PG). Participants were divided into three POAG groups according to visual field mean deviation: PG (no visual field defect); early glaucoma (EG; visual field mean deviation better than − 6 dB); and moderate-advanced glaucoma (MAG; visual field mean deviation worse than − 6 dB)52,53. Participants were excluded from the study if they displayed unreliable visual field testing, defined as fixation losses over 20% or false positive/negative errors over 20%54,55.

Spectral domain optical coherence tomography (SD-OCT) scans (20° scan angle) centred on the ONH, comprising 512 × 512 × 1024 pixels, were acquired in both eyes of each participant using long wavelength SD-OCT. The latter was a custom-built research OCT device, with a superluminescent diode light source (1 micron ASE module, NP-Photonics, Tuscan, US) with a centre wavelength of 1040 nm (bandwidth 70 nm), fibre-optically connected to a 20:80 beam-splitter serving the sample and reference arm respectively56. The system used a grating-based spectrometer with a Goodrich SUI-LDH-1.7 camera and operated at 47,000 a-scans per second with an axial resolution of ~ 7 µm and transverse resolution of ~ 15 µm57,58,59. Prior to OCT image acquisition, the power of the OCT system was recorded to ensure it was less than 2.5mW at the cornea, below the maximum power limit for a 10-s duration exposure60,61.

Acquired spectral data in FD1 file format were converted to 16-bit TIFF image format using custom-written software, OCT1_FD1_v2.2 (MATLAB 2014b, Math-works, US). OCT tomograms were then aligned within the 3D image datasets using stack registration using the Fiji ImageJ (version 1.52a, National Institutes of Health, USA); Schindelin et al.62 plugin StackReg63. OCT image noise was reduced using a 3D median filter. Image brightness and contrast were adjusted to redistribute pixel intensities for optimal visualisation of the ONH. Image pixel to distance calibration was performed, accounting for the effect of axial eye length and tissue refractive index64,65, as described by Terry et al.59.

Diagrams of optic nerve heads, indicating the position of axonal parameter measurements can be seen in Fig. 1a–d. The latter can also be observed in representative OCT tomograms of optic nerve heads within the groups: control, preperimetric glaucoma, early glaucoma and moderate-advanced glaucoma, with corresponding visual field plots in Fig. 1e.

(a) Schematic optic nerve head (ONH) diagram depicting (a) measurement of border nerve fibre layer (bNFL); (b) peripapillary nerve fibre layer (pNFL) measured 1.7 mm from ONH centre; (c) minimum rim width (MRW); (d) prelamina thickness calculated as the difference between anterior lamina cribrosa (LC) and prelamina surface depths relative to Bruch’s membrane opening (BMO). (e) Representative OCT tomograms (inferior-nasal to superior-temporal axis) of optic nerve heads in Control, Preperimetric, Early Glaucoma and Moderate-advanced Glaucoma groups, with corresponding visual field plots. Axon-related parameters: border NFL (red arrows), position of peripapillary NFL measurement (yellow arrow heads), Bruch’s membrane opening (orange line), MRW (blue arrows), prelamina thickness (green arrows). Scale bar represents 500 µm. Visual field plots are normal for control and preperimetric groups; note that the black region on right hand-side represents the physiological blind spot. Increasing visual field loss is indicated, from early glaucoma to moderate-advanced groups, by the increased areas of grey and black.

Each OCT dataset was sliced radially at 45° intervals around the ONH centre to produce four OCT tomograms through the 3D ONH image datasets, with orientations: superior – inferior, nasal – temporal, and superior nasal – inferior temporal, and superior temporal – inferior nasal, enabling measurement of 8 regions. Within each OCT tomogram, Bruch’s membrane terminations were identified and used to create a reference plane across Bruch’s membrane opening (BMO) for subsequent measurements.

A novel measure, NFL thickness at the ONH border (border NFL) was measured as the vertical distance from each Bruch’s membrane termination to the NFL surface (Fig. 1a). Peripapillary NFL was measured 1.7 mm from the BMO centre in each tomogram (Fig. 1b). Minimum rim width (MRW) was defined as the minimum distance from each Bruch’s membrane termination to the inner limiting membrane as described by Povazay et al.44, as shown in Fig. 1c. Minimum rim area (MRA) was calculated, using MRW values, as described by Gardiner et al.47.

Regional (superior, inferior, nasal, temporal, superior temporal, inferior temporal, superior temporal, inferior nasal and centre) prelamina thickness was measured within each ONH OCT dataset. The central region was measured at the midpoint of the BMO reference plane within each OCT tomogram, with other regions being measured at the mid-points between BMO centre and BMO terminations. Prelamina thickness was calculated by subtracting prelamina surface depth from anterior lamina cribrosa surface depth, both measured axially with respect to BMO reference plane (Fig. 1d).

Prelamina volume of each 3D ONH OCT image dataset was quantified using Amira software (version 6.0, Thermo Fisher Scientific, UK). First a reference BMO plane area was created by generate a 2D surface from landmarks placed around the ONH border. Then, using a similar a protocol with the addition of the anterior prelamina, 3D volumetric surface was generated to quantify the optic cup volume. Next, this was repeated using the anterior lamina cribrosa surface, and the optic cup volume was subtracted from this volume to calculate prelamina volume.

To overcome inter-eye correlation with the fellow eye, a common approach in ophthalmological research is either to average the data from fellow eyes, or to collect data from one eye only66. However, Murdoch et al.67,68 have suggested that it is ‘a waste’ to discard fellow eye data, and that appropriate statistical techniques, such as general linear mixed-effects models, can be applied to avoid data autocorrelation or having to discard data from the fellow eye68,69. Incorporating data from both eyes of a participant has the advantages of enhanced statistical power, more interpretable regression coefficients, greater precision of estimation and less sensitivity to missing data66,68,70.

In all statistical analyses performed in this study, data of each parameter (specific region or mean), acquired from both eyes (where possible) of each participant were used. To account for the inter-eye correlation of data between eyes of the same participant, general linear mixed-effects statistical models were constructed including a repeated measures component. The ‘lme4’ package was used to fit linear mixed-effects regression models within RStudio, version 1.2.1335 (http://cran.r-project.org/package=lme4). For each parameter measured e.g., superior border NFL, a linear mixed-effects regression model was developed. To optimise each model, a stepwise deletion of fixed effects (namely participant age, stage of glaucoma, axial eye length, anterior chamber depth, refractive error, central corneal thickness and intraocular pressure) was performed to determine the association between each fixed effect and the ONH parameter. Only significantly associated variables (p < 0.05) were included in regression models to account for effects on the axonal-related parameter. Inter-group differences (i.e., between the different stages of glaucoma) for each parameter were determined using Tukey post-hoc pairwise comparisons using ‘emmeans’ (http://cran.r- project.org/package = emmeans). Statistical differences were determined at p < 0.05.

Data normality was determined using histograms and the Shapiro–Wilk test with significance assumed at p < 0.05. Associations between each parameter and visual field mean deviation were determined using Pearson’s correlation coefficient.

Approved by South-East Wales Research Ethics Committee Panel C; reference number: 10/WSE03/24 and Cardiff University School of Optometry and Vision Sciences ethical committee; project number: 1299.

Informed written consent granted prior to participation.

None.

Regional border NFL inter-glaucoma stage differences and associations with visual field mean deviation are shown in Fig. 2a and b, respectively. Mean border NFL in the PG group was significantly less than in control eyes (p = 0.049). Compared to PG ONHs, average (p < 0.001), inferior (p = 0.007), and superior temporal (p < 0.001) border NFL were significantly less in EG ONHs. Border NFL was significantly less in MAG than in EG ONHs in all regions (p < 0.05), except temporal (p = 0.493). In all regions analysed, border NFL was significantly negatively correlated with VF MD (p < 0.001, Fig. 2b).

Regional border NFL as a function of glaucoma disease stage (a); black point represents mean and error bars indicate 95% confidence intervals, and (b) visual field status (VF MD). Blue line represents regression line and grey shading indicates 95% confidence intervals. Red text indicates significant Pearson’s correlation (r) values at p < 0.05. *p < 0.05, **p < 0.01, ***p < 0.001.

Regional peripapillary NFL inter-glaucoma stage differences and associations with VF MD are shown in Fig. 3a and b, respectively. Mean (p = 0.002), superior (p < 0.001), and inferior (p = 0.033) peripapillary NFL were significantly thinner in PG ONHs than in control eyes, and EG superior temporal peripapillary NFL was thinner than in PG ONHs (p = 0.023). MAG peripapillary NFL was thinner than EG in all regions (p < 0.05), except for superior (p = 0.674), superior nasal (p = 0.443), and inferior nasal (p = 0.135) peripapillary NFL. Peripapillary NFL was significantly negatively associated with VF MD in all regions (p < 0.001, Fig. 3b).

Regional pNFL as a function of glaucoma disease stage (a); black point represents mean and error bars indicate 95% confidence intervals, and (b) visual field status (VF MD). Blue line represents regression line and grey shading indicates 95% confidence intervals. Red text indicates significant Pearson’s correlation (r) values at p < 0.05. *p < 0.05, **p < 0.01, ***p < 0.001.

Regional MRW inter-glaucoma stage differences and associations with VF MD are shown in Fig. 4a and b, respectively. MRW was significantly thinner in PG than controls in the inferior (p = 0.027) and superior nasal (p = 0.037) regions. PG MRW was greater than in EG group in the inferior (p = 0.042), nasal (p = 0.014) and superior temporal (p = 0.036) regions, including the regional mean (p = 0.014, Fig. 4a). MRW was less in MAG, compared to EG in all regions (p < 0.05), except nasal (p = 0.208) and temporal (p = 0.547). MRW was significantly negatively associated with VF MD in all regions (p < 0.001, Fig. 4b).

Regional MRW as a function of glaucoma disease stage (a); black point represents mean and error bars indicate 95% confidence intervals, and (b) visual field status (VF MD). Blue line represents regression line and grey shading indicates 95% confidence intervals. Red text indicates significant Pearson’s correlation (r) values at p < 0.05. *p < 0.05, **p < 0.01, ***p < 0.001.

MRA did not significantly differ between PG and control eyes in any region (p > 0.300). Mean MRA (p = 0.011), and MRAs in nasal (p = 0.020), inferior temporal (p = 0.036), superior temporal (p = 0.026), and inferior nasal (p = 0.031) regions were significantly less in EG compared to PG. Mean MRA (p = 0.012), and superior (p < 0.001), and inferior (p = 0.012) MRA, were significantly less in MAG, compared to EG ONHs (Fig. 5a). MRA regions and mean MRA were significantly negatively correlated with VF MD (p < 0.01, Fig. 5b).

Regional MRA as a function of glaucoma disease stage (a); black point represents mean and error bars indicate 95% confidence intervals, and (b) visual field status (VF MD). Blue line represents regression line and grey shading indicates 95% confidence intervals. Red text indicates significant Pearson’s correlation (r) values at p < 0.05. *p < 0.05, **p < 0.01, ***p < 0.001.

The prelamina was thinner in PG than control eyes in all ONH regions (p < 0.05). No differences were identified between PG and EG (p > 0.20), or between MAG and EG prelamina thicknesses in any region (p > 0.07, Fig. 6a). A significant negative correlation between VF MD and prelamina thickness in all ONH regions (except centre (p = 0.127, Fig. 6b) was found (p < 0.05).

Regional prelamina thickness as a function of glaucoma disease stage (a); black point represents mean and error bars indicate 95% confidence intervals, and (b) visual field status (VFMD). Blue line represents regression line and grey shading indicates 95% confidence intervals. Red text indicates significant Pearson’s correlation (r) values at p < 0.05. *p < 0.05, **p < 0.01, ***p < 0.001.

Prelamina volume (Fig. 7a) did not significantly differ between controls and PG ONHs (p = 0.585), or between PG and EG ONHs (p = 0.681). However, prelamina volume was significantly less in MAG ONHs than in EG (p = 0.007), PG (p = 0.001) and control eyes (p < 0.001). A significant negative correlation was found between prelamina volume and VF MD (p < 0.001, Fig. 7b).

Prelamina volume as a function of glaucoma disease stage (a); black point represents mean and error bars indicate 95% confidence intervals, and (b) visual field status (VF MD). Blue line represents regression line and grey shading indicates 95% confidence intervals. Red text indicates significant Pearson’s correlation (r) values at p < 0.05. *p < 0.05, **p < 0.01, ***p < 0.001.

An association between increasing age and a decrease in mean peripapillary NFL, border NFL, MRW and MRA was determined (p < 0.01). Age was negatively associated with border NFL, MRW, and MRA in all regions (p < 0.05), except nasal (p = 0.245) and inferior nasal (p = 0.606) border NFL, nasal (p = 0.061) and inferior nasal (p = 0.412) MRW, and temporal (p = 0.108) and inferior nasal (p = 0.279) MRA. Additionally, age was found to have a negative association with prelamina thickness in the inferior and superior nasal regions (p < 0.05). Increasing axial length was negatively associated with mean, superior, and inferior MRW (p < 0.05), and border NFL in all regions (p < 0.05), except nasal (p = 0.083), temporal (p = 0.088), and inferior temporal (p = 0.191). CCT showed a positive association with the superior, superior temporal and inferior nasal prelaminar thickness (p < 0.05).

All association factors were accounted for in the relevant statistical models used to analyse the respective axonal parameters, described above, with respect to glaucoma disease stage.

Since vision loss in glaucoma cannot be recovered, early detection and diagnosis is essential. OCT evaluation of RNFL thickness has been increasingly adopted as an objective clinical measure for detection and assessment of glaucoma34,71. This study aimed to further investigate in vivo indices of RGC axons that alter in the earliest stages of glaucoma disease or may hold potential to characterise different stages of disease. This is the first study to quantify in vivo measurements of border NFL thickness and prelamina volume, and evaluate regional measures of prelamina thickness as a function of glaucoma disease stage.

The most important finding is that significant differences in border NFL, peripapillary NFL, MRW and prelamina thickness could be seen between control and PG ONHs, demonstrating that alterations in RGC axonal parameters can be quantified and detected in vivo, prior to permanent vision loss. Our finding consistent with previous studies reporting structural damage to the RNFL before clinically detectable loss of vision24,49,72. Indeed, Kerrigan-Baumrind et al.23 reported a loss of 25% to 35% of RGCs and their axons was associated with abnormalities detected by automated visual field testing.

The prelamina thickness was significantly less in PG ONHs compared with controls in all regions analysed. Average border NFL and peripapillary NFL, superior and inferior peripapillary NFL, as well as inferior and superior nasal MRW significantly decreased in PG ONHs, compared to controls. However, MRA did not significantly differ between PG ONHs and controls in any region analysed. Since border NFL, peripapillary NFL, and MRW were also significantly less in EG ONHs, compared to control ONHs, their potential use in identification of early-stage glaucoma is indicated.

Additionally, various specific regions of border NFL, peripapillary NFL, MRW, and MRA significantly differed between PG, EG, and MAG ONHs, suggesting that other axon-related indices may be useful in determining disease stage. Here, superior temporal peripapillary NFL in PG ONHs was thicker than in EG ONHs, and peripapillary NFL thickness differed between EG and MAG in all regions, except the superior, superior nasal, and inferior nasal regions, consistent with previous studies that reported peripapillary NFL to have good diagnostic sensitivity for glaucoma detection73,74,75, and also indicate disease progression76,77. In vivo measurements of peripapillary NFL have been shown to be preferable for glaucoma detection over macula retinal thickness parameters37,78,79,80. Furthermore, Sung et al.81 reported that peripapillary NFL thickness (automated measure by Cirrus SD-OCT) outperformed rim area, cup volume, and vertical cup-disc ratio for glaucoma discrimination, particularly in early glaucoma, whilst in advanced glaucoma, rim area and peripapillary NFL were comparable. Similarly, in moderate-advanced glaucoma participants (average VF MD: − 10.4 ± 8.5 dB), Mwanza et al.33 reported no difference in glaucoma diagnostic ability between ONH parameters and peripapillary NFL thickness measures acquired using the Cirrus SD-OCT device.

Gardiner et al.46 suggested that MRW and MRA may be more sensitive parameters for early glaucoma detection, with peripapillary NFL perhaps preferable for monitoring structural change. However, in this study, MRA did not significantly differ between PG and control eyes in any region analysed. We suggest that parameters such as border NFL, peripapillary NFL, and MRW may be better indicators of early glaucoma onset, whilst border NFL, peripapillary NFL, MRW, and MRA also provide insight into stage of disease.

Our study is the first to quantify specific regions of prelamina thickness and prelamina volume in vivo in human POAG. A significant reduction in regional prelamina thickness and volume was associated with decreasing VF sensitivity. Prelamina thickness was significantly less in PG and EG ONHs compared to controls, although did not significantly differ between PG and EG, or between EG and MAG in any ONH region. This suggests that specific regions of the prelamina thins prior to clinically detectable vision loss. Although no differences in prelamina volume was observed between early glaucoma stages or controls, prelamina volume significantly correlated with VF MD indicating a potential role in clinical monitoring of disease progression. Our findings are consistent with prelamina compression in early stages of glaucoma, followed by neural loss with advancement of glaucomatous optic neuropathy.

Multivariate analyses showed age and axial length to have significant associations with axon-related parameters, and that these factors should be accounted for when applied in a clinical setting for glaucoma assessment. The finding is not surprising in view of the known age-related decline in the density of RGC axons82 at a rate of 0.5% per year83 and reduced sensitivity across the visual field occurs as part of the normal ageing process84. The association of axial length with thinner average, superior, and inferior MRW, and most regions of border NFL could be attributed to larger eyes having larger ONHs85,86. MRA adjusts for ocular size and did not show any association with axial length. Additionally, since we report a positive association between CCT and prelaminar thickness in some regions, CCT should also be factored into analyses of OCT-derived measures of prelamina thickness. While prelamina measurements show promise, these measurements can be confounded by the ONH vasculature which results in significant shadowing and OCT signal attenuation.

In summary, here we report on novel quantitative measures of the optic nerve head as a function of different degrees of glaucomatous optic neuropathy as determined by classification of disease stage. Of the nine regions of prelamina thickness measured, all were able to distinguish between ONHs in the preperimetric group or early glaucoma group compared to healthy control optic nerve heads, indicating that prelamina thickness was an important parameter for detection of glaucoma before vision loss was detected. Additionally, the eight specific regions of border NFL were also able to differentiate between those ONHs with early glaucoma and healthy controls, whilst the mean value for border NFL was able to distinguish between successive stages of glaucoma, including between preperimetric glaucomatous optic nerve heads and controls, providing another novel measure for early detection prior to vision loss, and also for monitoring of disease stage. The latter was consistent with the border NFL having a significant correlation with visual field sensitivity, which was stronger than that of prelamina thickness. Prelamina volume, the third novel measure used in this study; was not able to distinguish between control and early-stage glaucoma.

Glaucoma diagnosis remains a challenge due to the slowly progressive nature of the disease and the sensitivity for testing for preperimetric glaucoma (reviewed by87. Our data reported here on measurements of specific regions of border NFL and prelamina thickness (previously not reported on) holds potential for a clinical role in prediction of optic nerves at risk of developing glaucomatous optic neuropathy and characterisation of disease stage. Additionally, their structure–function association with visual field sensitivity indicates their use in monitoring of disease stage. The concept of the novel border NFL thickness holds potential for clinical utility as it is derived as a measure from the RPE/Bruch’s membrane opening, considered to be a stable reference plane. The latter is recommended as a robust OCT-based disc margin88, thereby making parameters easier to segment (International Nomenclature for OCT panel89) easier to segment than traditional circumpapillary RNFL measures, as well as reliable even in tilted discs90.

All authors have full control of all primary data, and they agree to allow Scientific Reports to review their data upon request to the corresponding author.

Tham, Y. C. et al. Global prevalence of glaucoma and projections of glaucoma burden through 2040 a systematic review and meta-analysis. Ophthalmology 121(11), 2081–2090. https://doi.org/10.1016/j.ophtha.2014.05.013 (2014).

Article PubMed Google Scholar

Flaxman, S. R. et al. Global causes of blindness and distance vision impairment 1990–2020: a systematic review and meta-analysis. Lancet Glob. Health 5(12), E1221–E1234. https://doi.org/10.1016/s2214-109x(17)30393-5 (2017).

Article PubMed Google Scholar

Quigley, H. A. & Broman, A. T. The number of people with glaucoma worldwide in 2010 and 2020. Br. J. Ophthalmol. 90(3), 262–267. https://doi.org/10.1136/bjo.2005.081224 (2006).

Article CAS PubMed PubMed Central Google Scholar

Weinreb, R. N. et al. Primary open-angle glaucoma. Nat. Rev. Dis. Primers. 2, 1–19. https://doi.org/10.1038/nrdp.2016.67 (2016).

Article Google Scholar

Bourne, R. R. A. et al. Number of people blind or visually impaired by glaucoma worldwide and in world regions 1990–2010: a meta-analysis. Plos One https://doi.org/10.1371/journal.pone.0162229 (2016).

Article PubMed PubMed Central Google Scholar

Davis, B. M. et al. Glaucoma: the retina and beyond. Acta Neuropathol. 132(6), 807–826. https://doi.org/10.1007/s00401-016-1609-2 (2016).

Article PubMed PubMed Central Google Scholar

Gordon, M. O. et al. The Ocular Hypertension Treatment Study - Baseline factors that predict the onset of primary open-angle glaucoma. Arch. Ophthalmol. 120(6), 714–720 (2002).

Article PubMed Google Scholar

Leske, M. C. et al. Risk factors for incident open-angle glaucoma. Ophthalmology 115(1), 85–93. https://doi.org/10.1016/j.ophtha.2007.03.017 (2008).

Article PubMed Google Scholar

Garway-Heath, D. F. et al. Latanoprost for open-angle glaucoma (UKGTS): a randomised, multicentre, placebo-controlled trial. Lancet 385(9975), 1295–1304. https://doi.org/10.1016/s0140-6736(14)62111-5 (2015).

Article CAS PubMed Google Scholar

Crowston, J. G. & Weinreb, R. N. Glaucoma medication and aqueous humor dynamics. Curr. Opin. Ophthalmol. 16(2), 94–100 (2005).

Article PubMed Google Scholar

Heijl, A. et al. Reduction of intraocular pressure and glaucoma progression - results from the early manifest glaucoma trial. Arch. Ophthalmol. 120(10), 1268–1279 (2002).

Article PubMed Google Scholar

Morgan, J. E. Retina ganglion cell degeneration in glaucoma: an opportunity missed? A review. Clin. Exp. Ophthalmol. 40(4), 364–368. https://doi.org/10.1111/j.1442-9071.2012.02789.x (2012).

Article PubMed Google Scholar

Nickells, R. W. The cell and molecular biology of glaucoma: mechanisms of retinal ganglion cell death. Invest. Ophthalmol. Vis. Sci. 53(5), 2476–2481. https://doi.org/10.1167/iovs.12-9483h (2012).

Article CAS PubMed PubMed Central Google Scholar

Weinreb, R. N., Aung, T. & Medeiros, F. A. The pathophysiology and treatment of glaucoma: a review. JAMA 311(18), 1901–1911. https://doi.org/10.1001/jama.2014.3192 (2014).

Article CAS PubMed PubMed Central Google Scholar

Jonas, J. B. et al. Glaucoma. Lancet 390(10108), 2183–2193. https://doi.org/10.1016/s0140-6736(17)31469-1 (2017).

Article PubMed Google Scholar

Jonas, J. B., Gusek, G. C. & Naumann, G. O. H. Optic disc morphometry in chronic primary open-angle glaucoma. 1. Morphometric intrapapillary characteristics. Graefes Arch. Clin. Exp. Ophthalmol. 226(6), 522–530. https://doi.org/10.1007/bf02169199 (1988).

Article CAS PubMed Google Scholar

Schuster, A. K. et al. The diagnosis and treatment of glaucoma. Deutsches Arzteblatt International 117(13), 225–234. https://doi.org/10.3238/arztebl.2020.0225 (2020).

Article PubMed Google Scholar

Quigley, H. A. Glaucoma. Lancet 377(9774), 1367–1377. https://doi.org/10.1016/s0140-6736(10)61423-7 (2011).

Article PubMed Google Scholar

Budde, W.M. and J.B. Jonas, Morphology of the optic disc in glaucoma. I. Primary open-angle glaucomas. Klinische Monatsblatter Fur Augenheilkunde, 1999. 215(4): p. 211–220. https://doi.org/10.1055/s-2008-1034702

Bowd, C. et al. The retinal nerve fiber layer thickness in ocular hypertensive, normal, and glaucomatous eyes with optical coherence tomography. Arch. Ophthalmol. 118(1), 22–26 (2000).

Article CAS PubMed Google Scholar

Yu, M. et al. Risk of visual field progression in glaucoma patients with progressive retinal nerve fiber layer thinning a 5-year prospective study. Ophthalmology 123(6), 1201–1210. https://doi.org/10.1016/j.ophtha.2016.02.017 (2016).

Article PubMed Google Scholar

Miki, A. et al. Rates of retinal nerve fiber layer thinning in glaucoma suspect eyes. Ophthalmology 121(7), 1350–1358. https://doi.org/10.1016/j.ophtha.2014.01.017 (2014).

Article PubMed Google Scholar

Kerrigan-Baumrind, L. A. et al. Number of ganglion cells in glaucoma eyes compared with threshold visual field tests in the same persons. Invest. Ophthalmol. Vis. Sci. 41(3), 741–748 (2000).

CAS PubMed Google Scholar

Harwerth, R. S. et al. Linking structure and function in glaucoma. Prog. Retin. Eye Res. 29(4), 249–271. https://doi.org/10.1016/j.preteyeres.2010.02.001 (2010).

Article CAS PubMed PubMed Central Google Scholar

Hong, S. W. et al. Glaucoma specialist optic disc margin, rim margin, and rim width discordance in glaucoma and glaucoma suspect eyes. Am. J. Ophthalmol. 192, 65–76. https://doi.org/10.1016/j.ajo.2018.04.022 (2018).

Article PubMed PubMed Central Google Scholar

Varma, R., Steinmann, W. C. & Scott, I. U. Expert agreement in evaluating the optic disk for glaucoma. Ophthalmology 99(2), 215–221 (1992).

Article CAS PubMed Google Scholar

Jampel, H. D. et al. Agreement among glaucoma specialists in assessing progressive disc changes from photographs in open-angle glaucoma patients. Am. J. Ophthalmol. 147(1), 39–44. https://doi.org/10.1016/j.ajo.2008.07.023 (2009).

Article PubMed Google Scholar

Breusegem, C. et al. Agreement and accuracy of non-expert ophthalmologists in assessing glaucomatous changes in serial stereo optic disc photographs. Ophthalmology 118(4), 742–746. https://doi.org/10.1016/j.ophtha.2010.08.019 (2011).

Article PubMed Google Scholar

Drexler, W. & Fujimoto, J. G. State-of-the-art retinal optical coherence tomography. Prog. Retin. Eye Res. 27(1), 45–88. https://doi.org/10.1016/j.preteyeres.2007.07.005 (2008).

Article PubMed Google Scholar

Schmitt, J. M. Optical coherence tomography (OCT): A review. IEEE J. Sel. Top. Quantum Electron. 5(4), 1205–1215. https://doi.org/10.1109/2944.796348 (1999).

Article ADS CAS Google Scholar

Sharma, P. et al. Diagnostic tools for glaucoma detection and management. Surv. Ophthalmol. 53, S17–S32. https://doi.org/10.1016/j.survophthal.2008.08.003 (2008).

Article PubMed PubMed Central Google Scholar

Loewen, N. A. et al. Combining measurements from three anatomical areas for glaucoma diagnosis using Fourier-domain optical coherence tomography. Br. J. Ophthalmol. 99(9), 1224–1229. https://doi.org/10.1136/bjophthalmol-2014-305907 (2015).

Article PubMed Google Scholar

Mwanza, J. C. et al. Ability of Cirrus HD-OCT optic nerve head parameters to discriminate normal from glaucomatous eyes. Ophthalmology 118(2), 241-U34. https://doi.org/10.1016/j.ophtha.2010.06.036 (2011).

Article PubMed Google Scholar

Fortune, B. Optical coherence tomography evaluation of the optic nerve head neuro-retinal rim in glaucoma. Clin. Exp. Optom. 102(3), 286–290. https://doi.org/10.1111/cxo.12833 (2019).

Article PubMed Google Scholar

Medeiros, F. A. et al. Evaluation of retinal nerve fiber layer, optic nerve head, and macular thickness measurements for glaucoma detection using optical coherence tomography. Am. J. Ophthalmol. 139(1), 44–55. https://doi.org/10.1016/j.ajo.2004.08.069 (2005).

Article PubMed Google Scholar

Medeiros, F. A. et al. Detection of glaucoma progression with stratus OCT retinal nerve fiber layer, optic nerve head, and macular thickness measurements. Invest. Ophthalmol. Vis. Sci. 50(12), 5741–5748. https://doi.org/10.1167/iovs.09-3715 (2009).

Article PubMed Google Scholar

Lisboa, R. et al. Comparison of different spectral domain OCT scanning protocols for diagnosing preperimetric glaucoma. Invest. Ophthalmol. Vis. Sci. 54(5), 3417–3425. https://doi.org/10.1167/iovs.13-11676 (2013).

Article PubMed PubMed Central Google Scholar

Lisboa, R. et al. Diagnosing preperimetric glaucoma with spectral domain optical coherence tomography. Ophthalmology 119(11), 2261–2269. https://doi.org/10.1016/j.ophtha.2012.06.009 (2012).

Article PubMed Google Scholar

Schuman, J. S. et al. Quantification of nerve-fiber layer thickness in normal and glaucomatous eyes using optical coherence tomography - a pilot-study. Arch. Ophthalmol. 113(5), 586–596. https://doi.org/10.1001/archopht.1995.01100050054031 (1995).

Article CAS PubMed Google Scholar

Jonas, J. B., Fernandez, M. C. & Sturmer, J. Pattern of glaucomatous neuroretinal rim loss. Ophthalmology 100(1), 63–68 (1993).

Article CAS PubMed Google Scholar

Medeiros, F. A. et al. A combined index of structure and function for staging glaucomatous damage. Arch. Ophthalmol. 130(9), 1107–1116. https://doi.org/10.1001/archophthalmol.2012.827 (2012).

Article PubMed PubMed Central Google Scholar

Burgoyne, C. F. & Downs, J. C. Premise and prediction - How optic nerve head biomechanics underlies the susceptibility and clinical behavior of the aged optic nerve head. J. Glaucoma 17(4), 318–328. https://doi.org/10.1097/IJG.0b013e31815a343b (2008).

Article PubMed PubMed Central Google Scholar

Burgoyne, C. F. A biomechanical paradigm for axonal insult within the optic nerve head in aging and glaucoma. Exp. Eye Res. 93(2), 120–132. https://doi.org/10.1016/j.exer.2010.09.005 (2011).

Article CAS PubMed Google Scholar

Povazay, B. et al. Minimum distance mapping using three-dimensional optical coherence tomography for glaucoma diagnosis. J. Biomed. Opt. https://doi.org/10.1117/1.2773736 (2007).

Article PubMed Google Scholar

Chauhan, B. C. et al. Enhanced detection of open-angle glaucoma with an anatomically accurate optical coherence tomography-derived neuroretinal rim parameter. Ophthalmology 120(3), 535–543. https://doi.org/10.1016/j.ophtha.2012.09.055 (2013).

Article MathSciNet PubMed Google Scholar

Gardiner, S. K. et al. Structural measurements for monitoring change in glaucoma: comparing retinal nerve fiber layer thickness with minimum rim width and area. Invest. Ophthalmol. Vis. Sci. 56(11), 6886–6891. https://doi.org/10.1167/iovs.15-16701 (2015).

Article PubMed PubMed Central Google Scholar

Gardiner, S. K. et al. A method to estimate the amount of neuroretinal rim tissue in glaucoma: comparison with current methods for measuring rim area. Am. J. Ophthalmol. 157(3), 540–549. https://doi.org/10.1016/j.ajo.2013.11.007 (2014).

Article PubMed Google Scholar

Amini, N. et al. Structure-function relationships in perimetric glaucoma: comparison of minimum-rim width and retinal nerve fiber layer parameters. Invest. Ophthalmol. Vis. Sci. 58(11), 4623–4631. https://doi.org/10.1167/iovs.17-21936 (2017).

Article PubMed PubMed Central Google Scholar

Quigley, H. A., Addicks, E. M. & Green, W. R. Optic nerve damage in human glaucoma. III. Quantitative correlation of nerve fiber loss and visual field defect in glaucoma, ischemic neuropathy, papilledema, and toxic neuropathy. Arch. Ophthalmol. 100(1), 135–46 (1982).

Article CAS PubMed Google Scholar

Kitazawa, Y. & Yamamoto, T. Glaucomatous visual field defects: Their characteristics and how to detect them. Clin. Neurosci. 4(5), 279–283 (1997).

CAS PubMed Google Scholar

Keltner, J. L. et al. The association between glaucomatous visual fields and optic features in the ocular treatment study nerve head hypertension. Ophthalmology 113(9), 1603–1612. https://doi.org/10.1016/j.ophtha.2006.05.061 (2006).

Article PubMed Google Scholar

Budenz, D. L. et al. Comparison of glaucomatous visual field defects using standard full threshold and Swedish interactive threshold algorithms. Arch. Ophthalmol. 120(9), 1136–1141 (2002).

Article PubMed Google Scholar

Mills, R. P. et al. Categorizing the stage of glaucoma from pre-diagnosis to end-stage disease. Am. J. Ophthalmol. 141(1), 24–30. https://doi.org/10.1016/j.ajo.2005.07.044 (2006).

Article PubMed Google Scholar

Vingrys, A. J. & Demirel, S. False-response monitoring during automated perimetry. Optom. Vis. Sci. 75(7), 513–517. https://doi.org/10.1097/00006324-199807000-00020 (1998).

Article CAS PubMed Google Scholar

Newkirk, M. R. et al. Assessment of false positives with the humphrey field analyzer II perimeter with the SITA algorithm. Invest. Ophthalmol. Vis. Sci. 47(10), 4632–4637. https://doi.org/10.1167/iovs.05-1598 (2006).

Article PubMed Google Scholar

Povazay, B. et al. Three-dimensional optical coherence tomography at 1050 nm versus 800 nm in retinal pathologies: enhanced performance and choroidal penetration in cataract patients. J. Biomed. Opt. 12(4), 1–7. https://doi.org/10.1117/1.2773728 (2007).

Article Google Scholar

Esmaeelpour, M. et al. Three-dimensional 1060-nm OCT: choroidal thickness maps in normal subjects and improved posterior segment visualization in cataract patients. Invest. Ophthalmol. Vis. Sci. 51(10), 5260–5266. https://doi.org/10.1167/iovs.10-5196 (2010).

Article PubMed Google Scholar

Wood, A. et al. Retinal and choroidal thickness in early age-related macular degeneration. Am. J. Ophthalmol. 152(6), 1030–1038. https://doi.org/10.1016/j.ajo.2011.05.021 (2011).

Article PubMed Google Scholar

Terry, L. et al. Automated retinal layer segmentation using spectral domain optical coherence tomography: evaluation of inter-session repeatability and agreement between devices. PLoS ONE 11(9), 1–15. https://doi.org/10.1371/journal.pone.0162001 (2016).

Article CAS Google Scholar

ANSI, Safe use of lasers & safe use of optical fiber communications. American National Standards Institute - Z136 Committee. Orlando Laser Institute America, 2000: p. 168–173.

ICNIRP, Revision of guidelines on limits of exposure to laser radiation of wavelengths between 400nm and 1.4um. Health Physics, 2000. 79(9): p. 431–440.

Schindelin, J. et al. Fiji: an open-source platform for biological-image analysis. Nat. Methods 9(7), 676–682. https://doi.org/10.1038/nmeth.2019 (2012).

Article CAS PubMed Google Scholar

Thevenaz, P., Ruttimann, U. E. & Unser, M. A pyramid approach to subpixel registration based on intensity. IEEE Trans. Image Process. 7(1), 27–41. https://doi.org/10.1109/83.650848 (1998).

Article ADS CAS PubMed Google Scholar

Littmann, H. Determination of the real size of an object on the fundus of the living eye. Klin. Monatsbl. Augenheilkd. 180(4), 286–289. https://doi.org/10.1055/s-2008-1055068 (1982).

Article CAS PubMed Google Scholar

Bennett, A. G., Rudnicka, A. R. & Edgar, D. F. Improvements on littmann method of determining the size of retinal features by fundus photography. Graefes Arch. Clin. Exp. Ophthalmol. 232(6), 361–367. https://doi.org/10.1007/bf00175988 (1994).

Article CAS PubMed Google Scholar

Yu, F. & Afifi, A. A. Descriptive statistics in ophthalmic research. Am. J. Ophthalmol. 147(3), 389–391. https://doi.org/10.1016/j.ajo.2008.12.018 (2009).

Article PubMed Google Scholar

Murdoch, I. E., Morris, S. S. & Cousens, S. N. People and eyes: statistical approaches in ophthalmology. Br. J. Ophthalmol. 82(8), 971–973. https://doi.org/10.1136/bjo.82.8.971 (1998).

Article CAS PubMed PubMed Central Google Scholar

Burton, P., Gurrin, L. & Sly, P. Extending the simple linear regression model to account for correlated responses: an introduction to generalized estimating equations and multi-level mixed modelling. Stat. Med. 17(11), 1261–1291. https://doi.org/10.1002/(sici)1097-0258(19980615)17:11%3c1261::aid-sim846%3e3.0.co;2-z (1998).

3.0.co;2-z" data-track-action="article reference" href="https://doi.org/10.1002%2F%28sici%291097-0258%2819980615%2917%3A11%3C1261%3A%3Aaid-sim846%3E3.0.co%3B2-z" aria-label="Article reference 68" data-doi="10.1002/(sici)1097-0258(19980615)17:113.0.co;2-z">Article CAS PubMed Google Scholar

Fan, Q., Teo, Y. Y. & Saw, S. M. Application of advanced statistics in ophthalmology. Invest. Ophthalmol. Vis. Sci. 52(9), 6059–6065. https://doi.org/10.1167/iovs.10-7108 (2011).

Article PubMed Google Scholar

Glynn, R. J. & Rosner, B. Accounting for the correlation between fellow eyes in regression analysis. Arch. Ophthalmol. 110(3), 381–387. https://doi.org/10.1001/archopht.1992.01080150079033 (1992).

Article CAS PubMed Google Scholar

Na, J. H. et al. Detection of glaucomatous progression by spectral-domain optical coherence tomography. Ophthalmology 120(7), 1388–1395. https://doi.org/10.1016/j.ophtha.2012.12.014 (2013).

Article PubMed Google Scholar

Harwerth, R. S. & Quigley, H. A. Visual field defects and retinal ganglion cell losses in patients with glaucoma. Arch. Ophthalmol. 124(6), 853–859. https://doi.org/10.1001/archopht.124.6.853 (2006).

Article PubMed PubMed Central Google Scholar

Leung, C. K. S. et al. Retinal nerve fiber layer imaging with spectral-domain optical coherence tomography analysis of the retinal nerve fiber layer map for glaucoma detection. Ophthalmology 117(9), 1684–1691. https://doi.org/10.1016/j.ophtha.2010.01.026 (2010).

Article PubMed Google Scholar

Li, S. F. et al. Evaluation of optic nerve head and retinal nerve fiber layer in early and advance glaucoma using frequency-domain optical coherence tomography. Graefes Arch. Clin. Exp. Ophthalmol. 248(3), 429–434. https://doi.org/10.1007/s00417-009-1241-0 (2010).

Article PubMed Google Scholar

Oddone, F. et al. Influence of disc size on optic nerve head versus retinal nerve fiber layer assessment for diagnosing glaucoma. Ophthalmology 118(7), 1340–1347. https://doi.org/10.1016/j.ophtha.2010.12.017 (2011).

Article PubMed Google Scholar

Leung, C. K. S. et al. Evaluation of retinal nerve fiber layer progression in glaucoma: a study on optical coherence tomography guided progression analysis. Invest. Ophthalmol. Vis. Sci. 51(1), 217–222. https://doi.org/10.1167/iovs.09-3468 (2010).

Article MathSciNet PubMed Google Scholar

Bussel, I. I., Wollstein, G. & Schuman, J. S. OCT for glaucoma diagnosis, screening and detection of glaucoma progression. Br. J. Ophthalmol. 98, 15–19. https://doi.org/10.1136/bjophthalmol-2013-304326 (2014).

Article Google Scholar

Sung, K. R. et al. Macular assessment using optical coherence tomography for glaucoma diagnosis. Br. J. Ophthalmol. 96(12), 1452–1455. https://doi.org/10.1136/bjophthalmol-2012-301845 (2012).

Article PubMed Google Scholar

Fang, Y. A. et al. Diagnostic capability of Fourier-Domain optical coherence tomography in early primary open angle glaucoma. Chin. Med. J. 123(15), 2045–2050. https://doi.org/10.3760/cma.j.issn.0366-6999.2010.15.017 (2010).

Article PubMed Google Scholar

Oddone, F. et al. Macular versus retinal nerve fiber layer parameters for diagnosing manifest glaucoma a systematic review of diagnostic accuracy studies. Ophthalmology 123(5), 939–949. https://doi.org/10.1016/j.ophtha.2015.12.041 (2016).

Article PubMed Google Scholar

Sung, K. R., Na, J. H. & Lee, Y. Glaucoma diagnostic capabilities of optic nerve head parameters as determined by Cirrus HD optical coherence tomography. J. Glaucoma 21(7), 498–504. https://doi.org/10.1097/IJG.0b013e318220dbb7 (2012).

Article PubMed Google Scholar

Leung, C. K. S. et al. Retinal nerve fiber layer imaging with spectral-domain optical coherence tomography: a prospective analysis of age-related loss. Ophthalmology 119(4), 731–737. https://doi.org/10.1016/j.ophtha.2011.10.010 (2012).

Article PubMed Google Scholar

Harwerth, R. S., Wheat, J. L. & Rangaswamy, N. V. Age-related losses of retinal ganglion cells and axons. Invest. Ophthalmol. Vis. Sci. 49(10), 4437–4443. https://doi.org/10.1167/iovs.08-1753 (2008).

Article PubMed Google Scholar

Johnson, C. A., Adams, A. J. & Lewis, R. A. Evidence for a neural basis of age-related visual field loss in normal observers. Invest. Ophthalmol. Vis. Sci. 30(9), 2056–2064 (1989).

CAS PubMed Google Scholar

Oliveira, C. et al. Axial length and optic disc size in normal eyes. Br. J. Ophthalmol. 91(1), 37–39. https://doi.org/10.1136/bjo.2006.102061 (2007).

Article MathSciNet CAS PubMed Google Scholar

Jonas, J. B. & Xu, L. Histological changes of high axial myopia. Eye 28(2), 113–117. https://doi.org/10.1038/eye.2013.223 (2014).

Article CAS PubMed Google Scholar

Mwanza, J. C. & Budenz, D. L. New developments in optical coherence tomography imaging for glaucoma. Curr. Opin. Ophthalmol. 29(2), 121–129. https://doi.org/10.1097/icu.0000000000000452 (2018).

Article PubMed Google Scholar

Hwang, Y. H., Kim, M. K. & Ahn, S. I. Consistency of Bruch membrane opening detection as determined by optical coherence tomography. J. Glaucoma 25(11), 873–878. https://doi.org/10.1097/ijg.0000000000000448 (2016).

Article PubMed Google Scholar

Staurenghi, G. et al. Proposed lexicon for anatomic landmarks in normal posterior segment spectral-domain optical coherence tomography. Ophthalmology 121(8), 1572–1578. https://doi.org/10.1016/j.ophtha.2014.02.023 (2014).

Article PubMed Google Scholar

Lee, E. J. et al. Glaucoma diagnostic ability of the new circumpapillary retinal nerve fiber layer thickness analysis based on bruch’s membrane opening. Invest. Ophthalmol. Vis. Sci. 57(10), 4194–4204. https://doi.org/10.1167/iovs.16-19578 (2016).

Article CAS PubMed Google Scholar

Download references

Funded by Welsh Government through Health and Care Research Wales Health Studentship Award: HS-16–39 (stipend, fees, and research support for RLB PhD study). The views expressed are those of the authors and not necessarily those of Health and Care Research Wales or Welsh Government. Wellcome Trust Institutional Strategic Support Fund: 204824/Z/16/Z (salary for RLB). College of Optometrists (PhD studentship: stipend, fees, and research support for BEF PhD study). IGA/UKEGS award (salary for KEM, research support for development of methodology). Medical Research Council G0800547 (development of OCT).

N. White is deceased.

School of Optometry and Vision Sciences, Cardiff University, Cardiff, UK

R. L. Bartlett, B. E. Frost, K. E. Mortlock, J. R. Fergusson, N. White, J. E. Morgan, R. V. North & J. Albon

Cardiff Institute for Tissue Engineering and Repair, Cardiff University, Cardiff, UK

R. L. Bartlett, K. E. Mortlock, J. R. Fergusson, N. White & J. Albon

Vivat Scientia Bioimaging Laboratories, Cardiff University, Cardiff, UK

R. L. Bartlett, K. E. Mortlock, J. R. Fergusson, N. White, R. V. North & J. Albon

School of Biosciences, Cardiff University, Cardiff, UK

B. E. Frost

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

J.A., R.N., K.M., N.W., and J.M. contributed to the study conception. All authors contributed to the study design and development of techniques. Material preparation, data collection and analysis were performed by R.B., B. F., and J. F. The first draft of the manuscript was written by R. B. and all authors commented on previous versions of the manuscript. All authors read and approved the final manuscript.

Correspondence to J. Albon.

The authors declare no competing interests.

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

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://creativecommons.org/licenses/by/4.0/.

Reprints and Permissions

Bartlett, R.L., Frost, B.E., Mortlock, K.E. et al. Quantifying biomarkers of axonal degeneration in early glaucoma to find the disc at risk. Sci Rep 12, 9366 (2022). https://doi.org/10.1038/s41598-022-12036-4

Download citation

Received: 02 September 2021

Accepted: 31 March 2022

Published: 07 June 2022

DOI: https://doi.org/10.1038/s41598-022-12036-4

Anyone you share the following link with will be able to read this content:

Sorry, a shareable link is not currently available for this article.

Provided by the Springer Nature SharedIt content-sharing initiative

Molecular Biology Reports (2022)

By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.