# One year structural and functional glaucoma progression after trabeculectomy

Jacqueline Chua, Aistė Kadziauskienė, Damon Wong, Rimvydas Ašoklis, Eugenijus Lesinskas, Nguyen Duc Quang, Rachel Chong, Bingyao Tan, Michaël Girard, Jean Martial Mari, et al.

<figure><img src="image_1.png" /></figure>

## To cite this version:

Jacqueline Chua, Aistė Kadziauskienė, Damon Wong, Rimvydas Ašoklis, Eugenijus Lesinskas, et al. 2020, 10 (1), 10.1038/s41598-020-59792-9. [DOI](https://doi.org/10.1038/s41598-020-59792-9) [HAL](https://www.hal.science/hal-03133394v1) [arXiv](https://arxiv.org/abs/2009.09377)

HAL Id: hal-03133394

https://hal.science/hal-03133394v1

Submitted on 15 Sep 2025

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_2.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
OPEN

# One year structural and functional glaucoma progression after trabeculectomy

Jacqueline Chua<sup>1,2</sup>, Aistė Kadziauskienė<sup>3,4</sup>, Damon Wong<sup>5,6</sup>, Rimvydas Ašoklis<sup>3,4</sup>, Eugenijus Lesinskas<sup>3,4</sup>, Nguyen Duc Quang<sup>1</sup>, Rachel Chong<sup>1,2</sup>, Bingyao Tan<sup>1,5</sup>, Michaël J. A. Girard<sup>1,7</sup>, Jean Martial Mari<sup>1,8</sup>, Jonathan G. Crowston<sup>1,2</sup>, Tin Aung<sup>1,2,9</sup> & Leopold Schmetterer<sup>1,2,6,10,11*</sup>

We evaluated the changes in visual field mean deviation (VF MD) and retinal nerve fibre layer (RNFL) thickness in glaucoma patients undergoing trabeculectomy. One hundred patients were examined with VF and spectral-domain optical coherence tomography (OCT) before trabeculectomy and 4 follow-up visits over one year. Linear mixed models were used to investigate factors associated with VF and RNFL. VF improved during the first 3 months of follow-up (2.55 ± 1.06 dB/year) and worsened at later visits (−1.14 ± 0.29 dB/year). RNFL thickness reduced by −4.21 ± 0.25 µm/year from 1st month of follow-up. Eyes with an absence of initial VF improvement (β = 0.64; 0.30–0.98), RNFL thinning (β = 0.15; 0.08–0.23), increasing intraocular pressure (IOP; β = −0.11; −0.18 to −0.03) and severe glaucoma (β = −10.82; −13.61 to −8.02) were associated with VF deterioration. Eyes with VF deterioration (β = 0.19; 0.08–0.29), increasing IOP (β = −0.09; −0.17 to −0.01), and moderate (β = −6.33; −12.17 to −0.49) or severe glaucoma (β = −19.58; −24.63 to −14.52) were associated with RNFL thinning. Changes in RNFL structure and function occur over a 1-year follow-up period after trabeculectomy. Early VF improvement is more likely to occur in patients with mild/moderate glaucoma, whereas those with severe glaucoma show greater decline over one year. Our findings indicate that progression is observable using OCT, even in late-stage glaucoma.

Clinically, reduction of intraocular pressure (IOP) remains the only proven treatment for glaucoma<sup>1–3</sup>. Landmark glaucoma randomized clinical trials have confirmed the value of lowering IOP in delaying visual field (VF) progression<sup>1–5</sup>. Interestingly, recent studies<sup>6–8</sup> have reported that the surgical lowering of IOP in glaucoma can result in the improvement of VF, which lends support to the preliminary data reported by Katz and Spaeth, who demonstrated reversal of structural change and field loss in newly diagnosed patients after commencement of treatment<sup>9,10</sup>.

Apart from functional changes, others have examined structural changes after trabeculectomy and their potential association with functional changes. The association between structural glaucomatous damage with VF decline after trabeculectomy is complex. The MoreFlow Medical Research Council 5-Fluorouracil (5-FU) study reported that out of 250 eyes, 20 eyes showed structural progression, 35 eyes showed VF progression, and 15 eyes showed both over a five years period<sup>11</sup>. In the Collaborative Initial Glaucoma Treatment Study (CIGTS), VF worsening was significantly associated with enlargement of the optic cup, but reversal of cupping was not associated

<footnote n="1"><sup>1</sup>Singapore Eye Research Institute, Singapore National Eye Centre, Singapore, Singapore. <sup>2</sup>Academic Clinical Program, Duke-NUS Medical School, Singapore, Singapore. <sup>3</sup>Clinic of Ears, Nose, Throat and Eye Diseases, Institute of Clinical Medicine, Faculty of Medicine, Vilnius University, Vilnius, Lithuania. <sup>4</sup>Department of Eye Diseases, Vilnius University Hospital Santaros Klinikos, Vilnius, Lithuania. <sup>5</sup>SERI-NTU Advanced Ocular Engineering (STANCE) Program, Singapore, Singapore. <sup>6</sup>Department of Ophthalmology, Lee Kong Chian School of Medicine, Nanyang Technological University, Singapore, Singapore. <sup>7</sup>Ophthalmic Engineering & Innovation Laboratory, Department of Biomedical Engineering, National University of Singapore, Singapore, Singapore. <sup>8</sup>GePaSud Laboratory, University of French Polynesia, Tahiti, French Polynesia. <sup>9</sup>Yong Loo Lin School of Medicine, National University of Singapore, Singapore, Singapore. <sup>10</sup>Department of Clinical Pharmacology, Medical University of Vienna, Vienna, Austria. <sup>11</sup>Center for Medical Physics and Biomedical Engineering, Medical University of Vienna, Vienna, Austria. *email: leopold.schmetterer@seri.com.sg</footnote>
<figure><img src="image_3.png" /><figcaption>**Figure 1.** Identification of eligible glaucoma patients who underwent trabeculectomy. Of the 124 patients enrolled, 100 were included in the analysis. Of which, majority were examined at 7 visits whereas 10 missed one follow-up visit. Visual field tests were performed for 5 visits (visit 1, 4–7) and OCT/IOP were performed at all 7 visits.</figcaption></figure>

with improvement of the VF¹². Other studies reported an early increase in retinal nerve fiber layer (RNFL) thickness after trabeculectomy, but did not investigate the relation with VF changes¹³⁻¹⁷. Also, these studies have used the older time-domain optical coherence tomography (OCT)¹⁴,¹⁶, stereoscopic optic disc photographs¹²,¹⁵ scanning laser polarimetry¹⁷, or Heidelberg Retina Tomography (HRT)¹¹,¹³, which required a manually drawn contour line and thus was more prone to operator induced variability than the newer spectral-domain OCT.

The purpose of the current study was to assess the impact of trabeculectomy and associated IOP lowering on rates of VF and RNFL thickness progression, and secondly to determine factors related to these changes following trabeculectomy.

## Results

Initially, 130 glaucomatous eyes of 124 patients undergoing trabeculectomy were enrolled (Fig. 1). Of these, 24 patients (25 eyes) were excluded because of pre-perimetric glaucoma (3 eyes), poor OCT quality (6 eyes), previous failed trabeculectomy (2 eyes), postoperative complications (1 eye), and inadequate number of VFs/OCTs measurements or duration of follow-up (13 eyes). Finally, we analyzed the data of 100 patients (105 eyes). Ten of these missed one follow-up visit.

The demographic and clinical baseline characteristics of the patients are presented in Table 1. The mean (±standard deviation) age of the patients' age was 67.28 ± 8.92 years and 47% female. Baseline VF MD value was -14.68 ± 8.43 dB (range, -7.18 to -20.55 dB), and RNFL thickness was 53.38 ± 13.62 µm (range, 43 to 60 µm). Majority of these patients had exfoliation syndrome and severe glaucoma.

The overall mean IOP before surgery was 27.36 ± 6.69 mmHg and 9.29 ± 3.98 mmHg after surgery, which represents a mean decrease of 66% (P < 0.001; Fig. 2A). VF MD improved during the first 3 months of follow-up (2.55 ± 1.06 dB/year; P = 0.016) and worsened at later follow-up visits (-1.14 ± 0.29 dB/year; P < 0.001; Fig. 2B and Table 2). RNFL thickness reduced from 1st month of follow-up onwards and the mean RNFL change rate was -4.21 ± 0.25 µm/year; P < 0.001 (Fig. 2C and Table 2).

## Multivariate Analysis

Factors associated with changes in VF MD during two periods were explored: (1) baseline to 3rd month (early, Tables 3) and (2) postoperative 3rd month to 12th month (late, Table 4). Eyes that tended to have RNFL thinning (β = 0.14, 95%CI, 0.06 to 0.23), lesser reduction in post-operatives IOP (β = -0.08, 95%CI, -0.12 to -0.04) and severe glaucoma (β = -11.88, 95%CI, -15.11 to -8.64; Table 3) were independently associated with earlier VF MD worsening. Eyes that tended to have an earlier VF MD worsening (β = 0.64, 95%CI, 0.30 to 0.98), RNFL thinning (β = 0.15, 95%CI, 0.08 to 0.23), lesser reduction in post-operatives IOP (β = -0.11, 95%CI, -0.18 to -0.03) and severe glaucoma (β = -10.82, 95%CI, -13.61 to -8.02; Table 4) were independently associated with later VF MD worsening. Figure 3 shows that those having an improvement in VF MD in early phase are less likely to have a worsening of their VF in later stages. Baseline IOP and lamina cribrosa features were not associated with VF MD change (P > 0.05; Tables 3 and 4).

We then examined the factors associated with changes in RNFL thickness during two periods: (1) baseline to 3rd month (early, Table 5) and (2) postoperative 3rd month to 12th month (late, Table 6). Eyes that tended to have a more pronounced worsening in VF (β = 0.26, 95%CI, 0.05 to 0.46), surprisingly, greater reduction in post-operatives IOP (β = 0.16, 95%CI, 0.11 to 0.22), pseudoexfoliative glaucoma than angle closure glaucoma (β = 13.88, 95%CI, 7.96 to 19.80), moderate (β = -9.31, 95%CI, -14.76 to -3.85) or severe glaucoma (β = -20.31, 95%CI, -25.75 to -14.87) than mild glaucoma and spherical cup shaped lamina cribrosa (β = 5.92, 95%CI, 0.53 to 11.31; Table 5) were independently associated with RNFL thinning in the early stage. Eyes that
<table><caption>**Table 1.** Patient Demographic and Clinical Characteristics. BCVA = best-corrected visual acuity; D = diopters; dB = decibels; IOP = intraocular pressure; MD = mean deviation; RNFL = retinal nerve fiber layer; VF = visual field. Values expressed as mean ± standard deviation, unless otherwise indicated.</caption><thead><tr><th>Baseline Characteristic</th><th>Glaucoma Eyes (n = 105 Eyes, 100 Participants)</th></tr></thead><tbody><tr><td>Age, years</td><td>67.28 ± 8.92</td></tr><tr><td>Gender (female:male)</td><td>47:53</td></tr><tr><td>BCVA, decimal scale</td><td>0.68 ± 0.25</td></tr><tr><td>Refractive error, D</td><td>-0.40 ± 1.85</td></tr><tr><td>Central corneal thickness, µm</td><td>519.45 ± 31.22</td></tr><tr><td>Axial eye length, mm</td><td>23.60 ± 0.93</td></tr><tr><td>IOP, mmHg</td><td>27.36 ± 6.69</td></tr><tr><td>VF MD, dB</td><td>-14.68 ± 8.43</td></tr><tr><td>Global RNFL thickness, µm</td><td>53.38 ± 13.62</td></tr><tr><td><strong>Lamina cribrosa</strong></td><td></td></tr><tr><td>Curvatures, mm<sup>-1</sup></td><td>0.59 ± 0.41</td></tr><tr><td>GSI</td><td>-0.62 ± 0.19</td></tr><tr><td>Depth, µm</td><td>451.84 ± 133.35</td></tr><tr><td><strong>Glaucoma subtype, n (%)</strong></td><td></td></tr><tr><td>Pseudoexfoliative</td><td>78 (74)</td></tr><tr><td>Primary open-angle</td><td>19 (18)</td></tr><tr><td>Primary angle closure</td><td>8 (8)</td></tr><tr><td><strong>Glaucoma severity, n (%)</strong></td><td></td></tr><tr><td>Mild</td><td>16 (15)</td></tr><tr><td>Moderate</td><td>17 (16)</td></tr><tr><td>Severe</td><td>72 (69)</td></tr></tbody></table>

tended to have a more pronounced worsening in VF (β = 0.19, 95%CI, 0.08 to 0.29), increasing post-operatives IOP (β = -0.09, 95%CI, -0.17 to -0.01), pseudoexfoliative glaucoma than angle closure glaucoma (β = 12.65, 95%CI, 6.42 to 18.88), moderate (β = -6.33, -12.17 to -0.49) or severe glaucoma (β = -19.58, 95%CI, -24.63 to -14.52) than mild glaucoma and decreasing lamina cribrosa depth (β = 0.014, 95%CI, 0.004 to 0.025; Table 6) were independently associated with RNFL thinning in the later stage.

## Discussion

The current study investigated the rates of progression of VF MD and RNFL thickness and factors related to changes in the VF MD and RNFL thickness in a clinical sample of glaucoma patients undergoing trabeculectomy. After surgical reduction in IOP, we identified an initial improvement in VF followed by a subsequent decline over 9 months. We also found that the later decline in VF MD correlated with the initial changes in VF MD, degree of IOP lowering and RNFL thinning. Finally, individuals with severe glaucoma had a worse decline in both VF and RNFL thickness. Understanding VF/RNFL thickness changes over time and their related factors may be relevant when managing patients with severe glaucoma who are at greatest risk of lifetime blindness.

A critical issue is whether changes in VF over the short term as observed in the present study are real or related to fluctuations. We have demonstrated an initial VF gain after surgical intervention that could also reflect a regression to the mean (RTM) rather than an actual VF change. There are reasons to believe that the initial VF improvement is real and not simply due to variability. Firstly, our participants underwent to two VF tests at baseline which would mitigate VF measurement inaccuracy. In addition, we selected the second VF report for analysis, which should hopefully reduce learning effects and reduce measurement variability. We further ran a separate analysis using the first VF report and significance of the measured effect size and directions remained similar with the model reported in Table 3. Secondly, those individuals were measured with worse baseline VF MD did not demonstrate a more distinct improvement in MD at three-month follow-up, which is suggestive of RTM. Instead, the VF improvement was seen markedly in eyes with mild or moderate glaucoma. Finally, we perceived significant correlations between VF improvement with IOP control and RNFL thickness, which would not be likely from random fluctuation associated with RTM phenomenon. The VF improvement over the 3 months after surgery was significantly related to RNFL thickness measurement, an objective and reproducible OCT biomarker, strongly indicates that the early VF improvement observed in this study is neither related to RTM nor with VF fluctuations.

Others have corroborated that VF improvement after glaucoma surgery as a real likelihood⁶⁻⁹. Spaeth was the first to report the effect of IOP lowering and the initial gain in VF following either trabeculectomy or argon laser trabeculoplasty⁹. Recent studies have reported short- and long-term VF improvement after glaucoma surgical intervention in individuals with mild to moderate glaucoma⁶,⁷.

We too observed an initial enhancement in VF in the first 3 months followed by a significant decline in the subsequent 9 months. This was seen even when IOP control remained good over the entire study period. However, it was important to state that higher IOP remained a risk factor for VF progression (Table 4). The most
<figure><img src="image_4.png" /><figcaption>**Figure 2.** Changes in (A) intraocular pressure, (2) visual field mean deviation and (3) retinal nerve fibre layer thickness after trabeculectomy. *Data shown at each time point are expressed as mean and P values were obtained using comparison with the previous follow-up using linear mixed model.</figcaption></figure>

<table><caption>**Table 2.** Rate of visual field mean deviation and retinal nerve fibre layer thickness change after trabeculectomy. dB = decibels; MD = mean deviation; RNFL = retinal nerve fiber layer; VF = visual field. Linear mixed model was used to adjust for age, scan quality of OCT, follow-up duration and intracorrelation between eyes for each patient.</caption><thead><tr><th>Characteristics</th><th>Rate of Change per Year (mean ± SE)</th><th>P value</th></tr></thead><tbody><tr><td>VF MD, dB</td><td></td><td></td></tr><tr><td>All visits</td><td>-0.43 ± 0.22</td><td><strong>0.049</strong></td></tr><tr><td>Baseline to 3 months postoperative</td><td>2.55 ± 1.06</td><td><strong>0.016</strong></td></tr><tr><td>Postoperative 3 months to 12 months</td><td>-1.14 ± 0.29</td><td><strong>&lt;0.001</strong></td></tr><tr><td>Global RNFL thickness, µm</td><td></td><td></td></tr><tr><td>Baseline to 12 months</td><td>-4.21 ± 0.25</td><td><strong>&lt;0.001</strong></td></tr></tbody></table>

important risk factor was, however, advanced glaucoma damage with some patients losing as much as 7 dB per year. The association of this VF loss with structural progression indicates that this is not due to cataract formation after surgery. Those patients having an initial improvement in VF were less likely to experience a VF decline
<table><caption>**Table 3.** Factors related to changes in visual field mean deviation (Early; Baseline to 3 months postoperative). IOP = intraocular pressure; RNFL = retinal nerve fiber layer; VF = visual field. Linear mixed model was used to adjust for patient cluster. *Multivariate model - adjusted for age, gender, quality of OCT, change in RNFL, change in IOP, baseline IOP, central corneal thickness, glaucoma severity and change in lamina cribrosa depth in the same model. †Those with severe glaucoma had -11.95 dB reduction in VF MD (95% CI, -14.83, -9.08; P < 0.001) in the univariate model and -9.27 dB reduction in VF MD (95% CI, -12.12, -6.42; P < 0.001) in the multivariate model compared with those with moderate glaucoma. ‡Changes in RNFL thickness (in µm) or IOP (in mmHg) or curvatures (mm - 1) or GSI or depth (µm) for 1 dB improvement in visual field mean deviation at early phase.</caption><thead><tr><th rowspan="3">Characteristics</th><th colspan="4">Per 1 dB change in visual field mean deviation</th></tr><tr><th colspan="2">Univariate</th><th colspan="2">Multivariate *</th></tr><tr><th>β (95%CI)</th><th>P value</th><th>β (95%CI)</th><th>P value</th></tr></thead><tbody><tr><td>Change in RNFL thickness, µm<sup>†</sup></td><td>0.29 (0.21, 0.37)</td><td><strong>&lt;0.001</strong></td><td>0.14 (0.06, 0.23)</td><td><strong>0.001</strong></td></tr><tr><td>Change in IOP, mmHg<sup>†</sup></td><td>-0.05 (-0.09, -0.02)</td><td><strong>0.001</strong></td><td>-0.08 (-0.12, -0.04)</td><td><strong>&lt;0.001</strong></td></tr><tr><td>Baseline IOP, mmHg</td><td>-0.17 (-0.41, 0.07)</td><td><strong>0.171</strong></td><td>-0.12 (-0.27, 0.02)</td><td>0.097</td></tr><tr><td>Central corneal thickness, µm</td><td>0.06 (0.01, 0.11)</td><td><strong>0.034</strong></td><td>0.03 (-0.01, 0.06)</td><td>0.105</td></tr><tr><td colspan="5"><strong>Glaucoma subtype</strong></td></tr><tr><td>Pseudoexfoliative</td><td colspan="4">Reference</td></tr><tr><td>Primary open-angle</td><td>-2.50 (-6.73, 1.73)</td><td>0.246</td><td>—</td><td>—</td></tr><tr><td>Primary angle-closure</td><td>0.17 (-5.96, 6.30)</td><td>0.957</td><td>—</td><td>—</td></tr><tr><td colspan="5"><strong>Glaucoma severity †</strong></td></tr><tr><td>Mild</td><td colspan="2">Reference</td><td colspan="2">Reference</td></tr><tr><td>Moderate</td><td>-3.61 (-7.32, 0.10)</td><td>0.056</td><td>-2.60 (-5.97, 0.76)</td><td>0.129</td></tr><tr><td>Severe</td><td>-15.56 (-18.51, -12.61)</td><td><strong>&lt;0.001</strong></td><td>-11.88 (-15.11, -8.64)</td><td><strong>&lt;0.001</strong></td></tr><tr><td colspan="5"><strong>Lamina cribrosa</strong></td></tr><tr><td>Changes in curvatures, mm<sup>-1</sup><sup>‡</sup></td><td>-1.08 (-2.35, 0.20)</td><td>0.099</td><td>—</td><td>—</td></tr><tr><td>Changes in GSI<sup>‡</sup></td><td>0.93 (-3.17, 5.03)</td><td>0.658</td><td>—</td><td>—</td></tr><tr><td>Changes in depth, µm<sup>‡</sup></td><td>-0.01 (-0.01, -0.01)</td><td><strong>0.039</strong></td><td>0.01 (-0.01, 0.01)</td><td>0.764</td></tr></tbody></table>

during the following 9 months. When comparing the VF changes over one year in the present study with those obtained in other trials¹,¹⁸⁻²², the more advanced glaucomatous damage in our population needs to be considered.

An important finding is that RNFL thickness significantly declined over the one-year study period. This was seen in spite of the mostly severe glaucoma cases and the very thin RNFL thickness at baseline. Cross-sectional studies have reported the “floor effect” of RNFL thickness at approximately 50 µm with Spectralis OCT, the level at which point no additional thinning can be distinguished²³⁻²⁵. Longitudinal studies have, however, shown that structural progression can still be seen in late stage disease and that the floor effects only occur at thinner RNFL values, which is in agreement with the present study²⁶,²⁷. Also, postmortem studies suggested that a few surviving ganglion cells remain detected in cases of advanced glaucoma²⁸. Hence, our data indicate that RNFL measurements can still be used to monitor VF progression in late stage glaucoma after trabeculectomy.

We²⁹ and others³⁰,³¹ have previously reported that the decline in RNFL thickness is dependent on the post-operative positioning of the lamina cribrosa. Also, we found that changes in VF were neither associated with the position or the shape of the lamina cribrosa. This makes lamina cribrosa parameters an unlikely candidate for sufficient structural progression analysis after trabeculectomy. This is in contrary to a study by Ha and co-workers, where they showed that the baseline mean anterior laminar cribrosa insertion depth was correlated with the rate of VF deterioration³². One plausible explanation may be related to the severity of glaucomatous damage. Ha *et al.* mostly had eyes with early glaucoma whereas ours were composed of eyes with late glaucoma. Hence, the behavior of the lamina cribrosa on VF worsening may be confounded by the glaucoma severity.

**Study strengths and limitations.** Strengths of the study include the prospective study design, large number of patients with severe glaucoma and spectral-domain OCT imaging. Compared with earlier studies, spectral-domain OCT is objective and reproducible and thus less prone to operator induced variability³³,³⁴. Limitations of the study include the relatively short observation period, small number of patients with mild or moderate disease and the lack of cataract staging after surgery. Approximately 48% of the eyes had mild cataract, 6% were pseudophakia, whereas the remaining 46% of the eyes did not have any diagnosis of lens disorders. Even though the severity of cataract was not graded, we would like to reiterate that we only included participants with good quality OCT scans. Those with poor quality OCT scans were excluded. Moreover, most of the patients had pseudoexfoliative glaucoma and it is unclear to which degree the results can be applicable to other study populations. Last, non-standard VF reliability criteria were used (e.g. 33% false positive errors instead of 15%). We selected a wider criteria of VF reliability in order to include more patients for analysis and because many patients
<table><caption>**Table 4.** Factors related to changes in visual field mean deviation (Late; Postoperative 3 months to 12 months). IOP = intraocular pressure; RNFL = retinal nerve fiber layer; VF = visual field. Linear mixed model was used to adjust for patient cluster. *Multivariate model - adjusted for age, gender, quality of OCT, change in RNFL, change in IOP, baseline IOP, central corneal thickness and glaucoma severity in the same model. † Those with severe glaucoma had -12.22 dB reduction in VF MD (95% CI, -15.18, -9.25; P < 0.001) in the univariate model and -8.39 dB reduction in VF MD (95% CI, -10.82, -5.95; P < 0.001) in the multivariate model compared with those with moderate glaucoma. ‡ Changes in early phase of VF MD (dB) or RNFL thickness (in µm) or IOP (in mmHg) or curvatures (mm<sup>-1</sup>) or GSI or depth (µm) for 1 dB improvement in visual field mean deviation in the late phase.</caption><thead><tr><th rowspan="3">Characteristics</th><th colspan="4">Per 1 dB change in visual field mean deviation</th></tr><tr><th colspan="2">Univariate</th><th colspan="2">Multivariate*</th></tr><tr><th>β (95%CI)</th><th>P value</th><th>β (95%CI)</th><th>P value</th></tr></thead><tbody><tr><td>Change in VF MD, dB<sup>‡</sup><br>(Postoperative 3 months minus baseline)</td><td>0.76 (0.17, 1.35)</td><td><strong>0.011</strong></td><td>0.64 (0.30, 0.98)</td><td><strong>&lt;0.001</strong></td></tr><tr><td>Change in RNFL thickness, µm<sup>‡</sup></td><td>0.37 (0.29, 0.44)</td><td><strong>&lt;0.001</strong></td><td>0.15 (0.08, 0.23)</td><td><strong>&lt;0.001</strong></td></tr><tr><td>Change in IOP, mmHg<sup>‡</sup></td><td>-0.12 (-0.19, -0.04)</td><td><strong>0.003</strong></td><td>-0.11 (-0.18, -0.03)</td><td><strong>0.009</strong></td></tr><tr><td>Baseline IOP, mmHg</td><td>-0.12 (-0.37, 0.13)</td><td>0.333</td><td>-0.04 (-0.18, 0.10)</td><td>0.615</td></tr><tr><td>Central corneal thickness, µm</td><td>0.05 (0.01, 0.11)</td><td><strong>0.041</strong></td><td>0.01 (-0.03, 0.03)</td><td>0.924</td></tr><tr><td><strong>Glaucoma subtype</strong></td><td></td><td></td><td></td><td></td></tr><tr><td>Pseudoexfoliative</td><td>Reference</td><td></td><td></td><td></td></tr><tr><td>Primary open-angle</td><td>-2.49 (-6.83, 1.85)</td><td>0.261</td><td>—</td><td>—</td></tr><tr><td>Primary angle-closure</td><td>0.01 (-6.28, 6.30)</td><td>0.997</td><td>—</td><td>—</td></tr><tr><td><strong>Glaucoma severity<sup>†</sup></strong></td><td></td><td></td><td></td><td></td></tr><tr><td>Mild</td><td>Reference</td><td></td><td>Reference</td><td></td></tr><tr><td>Moderate</td><td>-3.80 (-7.63, 0.04)</td><td>0.052</td><td>-2.43 (-5.37, 0.52)</td><td>0.106</td></tr><tr><td>Severe</td><td>-16.01 (-19.06, -12.97)</td><td><strong>&lt;0.001</strong></td><td>-10.82 (-13.61, -8.02)</td><td><strong>&lt;0.001</strong></td></tr><tr><td><strong>Lamina cribrosa</strong></td><td></td><td></td><td></td><td></td></tr><tr><td>Changes in curvatures, mm<sup>-1</sup><sup>†</sup></td><td>-0.23 (-1.72, 1.26)</td><td>0.762</td><td>—</td><td>—</td></tr><tr><td>Changes in GSI<sup>‡</sup></td><td>0.83 (-1.22, 2.87)</td><td>0.429</td><td>—</td><td>—</td></tr><tr><td>Changes in depth, µm<sup>‡</sup></td><td>-0.01 (-0.01, 0.01)</td><td>0.466</td><td>—</td><td>—</td></tr></tbody></table>

had late stage disease and therefore were unable to provide more reliable perimetry data. In spite of this limitation, we still observed VF changes over a period of time that were correlated to structural changes as assessed with OCT and therefore are likely to represent true changes rather than fluctuations.

In summary, we have shown that measurable changes occur in both VF MD and RNFL thickness after trabeculectomy. In the early phase, we observed an improvement in VF that was more likely to occur in patients with mild-moderate disease and less likely to occur in patients with RNFL thinning. From 3-12 months, a decline in VF was observed. Risk factors for the VF deterioration were increasing post-operatives IOP, late stage disease, RNFL thinning and absence of VF improvement in the first three months. Our data also indicate that even in late stage glaucoma, structural progression can be monitored with RNFL thickness measurement using spectral-domain OCT.

## Methods

This was a longitudinal study of consecutive patients, who were undergoing trabeculectomy at the Vilnius University Hospital Santaros Klinikos (2014-2017)²⁹. Approval for conducting the study was obtained from the Vilnius Regional Biomedical Research Ethics Committee and all study procedures adhered to recommendations of the Declaration of Helsinki. Written informed consents were obtained from participants. Inclusion criteria were defined as (1) clinical diagnosis of primary or secondary glaucoma, (2) trabeculectomy indicated because of progressing glaucoma or high risk of glaucoma progression due to high IOP (defined as IOP that is higher than patient's target IOP), (3) best corrected Snellen visual acuity of ≥ 0.1, (4) refractive error from -6.0 D to +6.0 D of sphere and ± 3.0 D of cylinder. Patients were excluded if they had any prior intraocular surgery, except phacoemulsification with intraocular lens implantation, other ophthalmological or neurological diseases affecting the VF, or poor image quality because of opaque ocular media. Glaucoma was defined based on the presence of glaucomatous optic neuropathy (neuroretinal rim thinning, notching or RNFL defects) with or without associated glaucomatous VF defect. A glaucomatous VF defect was defined as glaucoma hemifield test of standard automated perimetry outside normal limits and/or a cluster of at least three contiguous points on the pattern deviation plot with P < 5% and one with P < 1% probability of being normal or pattern standard deviation of <5%. The VF test was considered reliable if false positive and false negative errors were <33% and fixation losses <20%.

**Examination procedures.** Patients underwent baseline examination including, assessment of best-corrected visual acuity and autorefractometry (Topcon KR-1 Auto Kerato-Refractometer, Topcon Medical Systems, USA), slit-lamp biomicroscopy, dilated stereoscopic examination of fundus, and partial optical coherence interferometry (IOL Master, Carl Zeiss Meditec, Dublin, CA, USA) The mean of two IOP measurements
<figure><img src="image_5.png" /><figcaption>**Figure 3.** Scatterplot of the change in visual field mean deviation (VF MD) from baseline to 3rd month follow-up versus the change in VF MD from 1 year after surgery compared with before, stratified by glaucoma severity. Eyes having mild glaucoma were indicated with dark filled circle, those having moderate glaucoma in green hollow diamond and severe glaucoma in black hollow triangle. The VF before trabeculectomy and the VF postoperatively (3 months after surgery) were used for the calculation of VF worsening (r = 0.32; P = 0.005).</figcaption></figure>

<table><caption>**Table 5.** Factors related to changes in retinal nerve fibre layer thickness (Early; Baseline to 3 months postoperative). IOP = intraocular pressure; MD = mean deviation; RNFL = retinal nerve fiber layer; VF = visual field. Linear mixed model was used to adjust for patient cluster. *Multivariate model - adjusted for age, gender, quality of OCT, change in VF MD, change in IOP, baseline IOP, glaucoma severity and lamina cribrosa features (curved, GSI and depth) in the same model. ‡Changes in early phase of VF MD (dB) or IOP (in mmHg) or curvatures (mm<sup>-1</sup>) or GSI or depth (µm) for 1 µm increase in retinal nerve fibre layer thickness at early phase.</caption><thead><tr><th rowspan="3">Characteristics</th><th colspan="4">Per 1 µm change in retinal nerve fibre layer thickness</th></tr><tr><th colspan="2">Univariate</th><th colspan="2">Multivariate*</th></tr><tr><th>β (95%CI)</th><th>P value</th><th>β (95%CI)</th><th>P value</th></tr></thead><tbody><tr><td>Change in VF MD, dB<sup>†</sup><br>(Postoperative baseline to 3 months)</td><td>0.74 (0.55, 0.94)</td><td><strong>&lt;0.001</strong></td><td>0.26 (0.05, 0.46)</td><td><strong>0.015</strong></td></tr><tr><td>Change in IOP, mmHg<sup>†</sup></td><td>0.07 (0.04, 0.11)</td><td><strong>&lt;0.001</strong></td><td>0.16 (0.11, 0.22)</td><td><strong>&lt;0.001</strong></td></tr><tr><td>Baseline IOP, mmHg</td><td>0.09 (−0.28, 0.47)</td><td>0.621</td><td>0.18 (−0.05, 0.42)</td><td>0.125</td></tr><tr><td>Central corneal thickness, µm</td><td>0.06 (−0.02, 0.14)</td><td>0.122</td><td></td><td></td></tr><tr><td><strong>Glaucoma subtype</strong></td><td></td><td></td><td></td><td></td></tr><tr><td>Pseudoexfoliative</td><td>Reference</td><td></td><td>Reference</td><td></td></tr><tr><td>Primary open-angle</td><td>−1.46 (−7.84, 4.91)</td><td>0.653</td><td>3.26 (−0.92, 7.42)</td><td>0.126</td></tr><tr><td>Primary angle-closure</td><td>10.26 (1.00, 19.52)</td><td><strong>0.030</strong></td><td>13.88 (7.96, 19.80)</td><td><strong>&lt;0.001</strong></td></tr><tr><td><strong>Glaucoma severity</strong></td><td></td><td></td><td></td><td></td></tr><tr><td>Mild</td><td>Reference</td><td></td><td>Reference</td><td></td></tr><tr><td>Moderate</td><td>−7.48 (−13.86, −1.10)</td><td><strong>0.022</strong></td><td>−9.31 (−14.76, −3.85)</td><td><strong>0.001</strong></td></tr><tr><td>Severe</td><td>−22.89 (−27.95, −17.82)</td><td><strong>&lt;0.001</strong></td><td>−20.31 (−25.75, −14.87)</td><td><strong>&lt;0.001</strong></td></tr><tr><td><strong>Lamina cribrosa</strong></td><td></td><td></td><td></td><td></td></tr><tr><td>Changes in curvatures, mm<sup>−1</sup><sup>‡</sup></td><td>3.47 (2.20, 4.74)</td><td><strong>&lt;0.001</strong></td><td>−2.13 (−4.49, 0.22)</td><td>0.076</td></tr><tr><td>Changes in GSI<sup>‡</sup></td><td>5.52 (3.02, 8.01)</td><td><strong>&lt;0.001</strong></td><td>5.92 (0.53, 11.31)</td><td><strong>0.031</strong></td></tr><tr><td>Changes in depth, µm<sup>‡</sup></td><td>0.03 (0.02, 0.03)</td><td><strong>&lt;0.001</strong></td><td>0.004 (−0.008, 0.016)</td><td>0.514</td></tr></tbody></table>

using Goldmann applanation tonometry (GAT, Haag-Streit AG, Switzerland), spaced one minute apart was calculated. If two measurements differed by more than 2 mmHg, we took a third reading, and averaged the two closest values. One ophthalmologist examined the patients.

Automated perimetry using 30-2 Swedish Interactive Threshold Algorithm Standard strategy (Humphrey visual field analyzer, Carl Zeiss Meditec, Dublin, CA, USA) were performed at five visits (before the trabeculectomy and postoperatively 3 months, 6 months, 9 months and twelve months). Baseline VF was done twice to eliminate the known learning effect and the second VF report was used for current analysis if it was reliable. Glaucoma severity was staged based on the standard automated perimetry using the Hodapp-Parrish-Anderson criteria³⁵: mild glaucoma (mean deviation (MD)) > −6 dB, less than 25% of points are depressed <5% and less than 10 points are depressed <1% on a pattern deviation plot, all points in the central 5° with sensitivity ≥15 dB),
<table><caption>**Table 6.** Factors related to changes in retinal nerve fibre layer thickness (Late; Postoperative 3 months to 12 months). IOP = intraocular pressure; MD = mean deviation; RNFL = retinal nerve fibre layer; VF = visual field. Linear mixed model was used to adjust for patient cluster. *Multivariate model - adjusted for age, gender, quality of OCT, change in VF MD, change in IOP, baseline IOP, glaucoma severity and lamina cribrosa features (curved, GSI and depth) in the same model. ‡Changes in late phase of VF MD (dB) or IOP (in mmHg) or curvatures (mm<sup>-1</sup>) or GSI or depth (µm) for 1 µm increase in retinal nerve fibre layer thickness at late phase.</caption><thead><tr><th rowspan="3">Characteristics</th><th colspan="4">Per 1 µm change in retinal nerve fibre layer thickness</th></tr><tr><th colspan="2">Univariate</th><th colspan="2">Multivariate<sup>*</sup></th></tr><tr><th>Coefficients (95%CI)</th><th>P value</th><th>Coefficients (95%CI)</th><th>P value</th></tr></thead><tbody><tr><td>Change in VF MD, dB<sup>‡</sup> (Postoperative 3 months to 12 months)</td><td>0.30 (0.19, 0.40)</td><td>&lt;0.001</td><td>0.19 (0.08, 0.29)</td><td>0.001</td></tr><tr><td>Change in IOP, mmHg<sup>‡</sup></td><td>0.07 (0.04, 0.11)</td><td>&lt;0.001</td><td>-0.09 (-0.17, -0.01)</td><td>0.026</td></tr><tr><td>Baseline IOP, mmHg</td><td>0.09 (-0.28, 0.47)</td><td>0.621</td><td>0.08 (-0.17, 0.32)</td><td>0.538</td></tr><tr><td>Central corneal thickness, µm</td><td>0.06 (-0.01, 0.14)</td><td>0.108</td><td></td><td></td></tr><tr><td colspan="5"><strong>Glaucoma subtype</strong></td></tr><tr><td>Pseudoexfoliative</td><td>Reference</td><td></td><td>Reference</td><td></td></tr><tr><td>Primary open-angle</td><td>-1.46 (-7.84, 4.91)</td><td>0.653</td><td>2.57 (-1.96, 7.10)</td><td>0.267</td></tr><tr><td>Primary angle-closure</td><td>10.26 (1.00, 19.52)</td><td>0.030</td><td>12.65 (6.42, 18.88)</td><td>&lt;0.001</td></tr><tr><td colspan="5"><strong>Glaucoma severity</strong></td></tr><tr><td>Mild</td><td>Reference</td><td></td><td>Reference</td><td></td></tr><tr><td>Moderate</td><td>-7.48 (-13.86, -1.10)</td><td>0.022</td><td>-6.33 (-12.17, -0.49)</td><td>0.034</td></tr><tr><td>Severe</td><td>-22.89 (-27.95, -17.82)</td><td>&lt;0.001</td><td>-19.58 (-24.63, -14.52)</td><td>&lt;0.001</td></tr><tr><td colspan="5"><strong>Lamina cribrosa</strong></td></tr><tr><td>Changes in curvatures, mm<sup>-1</sup><sup>‡</sup></td><td>3.47 (2.20, 4.74)</td><td>&lt;0.001</td><td>-0.01 (-1.57, 1.54)</td><td>0.988</td></tr><tr><td>Changes in GSI<sup>‡</sup></td><td>5.52 (3.02, 8.01)</td><td>&lt;0.001</td><td>1.00 (-1.17, 3.16)</td><td>0.367</td></tr><tr><td>Changes in depth, µm<sup>‡</sup></td><td>0.03 (0.02, 0.03)</td><td>&lt;0.001</td><td>0.014 (0.004, 0.025)</td><td>0.007</td></tr></tbody></table>

moderate glaucoma (MD > −12dB, less than 50% of points are depressed <5% and less than 20 points are depressed <1% on the pattern deviation plot, only one hemifield having a point in the central 5° with sensitivity <15 dB, no points within 5° of fixation with sensitivity of 0 dB), severe glaucoma (MD < −12dB, more than 50% of points are depressed <5% and more than 20 points are depressed <1% on the pattern deviation plot, points within the central 5° with sensitivity <15 dB in both hemifields, at least one point with sensitivity of 0 dB within 5° of fixation).

All the patients underwent a limbal-based trabeculectomy, with or without adjunctive 5-fluorouracil, following the same surgical protocol by one of four surgeons. Subsequently, needling with 5-fluorouracil was performed if failure of the filtrating bleb occurred. Only patients with reduced postsurgical IOP continued the study.

**Spectral-domain optical coherence tomography.** We evaluated the RNFL thickness and lamina cribrosa features using the spectral-domain OCT (Heidelberg Spectralis, Heidelberg Engineering, Dossenheim, Germany) at seven visits (before the trabeculectomy and postoperatively 3–10 days, one month, three months, six months, nine months and twelve months). A 15 × 10° rectangle scan was centered on the optic nerve head. Each OCT volume consisted of 49 serial horizontal B scans (4.5 mm long lines, 40 images averaged) spaced at approximately 63 µm intervals. At least two OCT scans were taken and the one with the best quality was chosen. Images with a quality score ≤15 were excluded. The baseline OCT scan was set as a reference and all subsequent scans were done adherent to it. Potential magnification error was avoided by entering the corneal curvature and refraction of the eye before the OCT scanning. The RNFL thickness was measured automatically from the circumferential SD-OCT scan of 3.4 mm diameter centered at the ONH (single circle B scan of 12°, 100 images averaged).

To visualize the features of the lamina cribrosa, OCT images were enhanced using adaptive compensation (Reflectivity software, version 3.4, Ophthalmic Engineering & Innovation Laboratory, National University of Singapore, Singapore) and the measurements of lamina cribrosa were extracted using Morphology 1.0 software (Ophthalmic Engineering & Innovation Laboratory, National University of Singapore, Singapore)<sup>36–38</sup>. Curvature was expressed as the values in mm<sup>-1</sup>, negative values describing posteriorly curved lamina cribrosa and positive values indicating anteriorly curved lamina cribrosa<sup>39</sup>. Global shape index ranged between −1 and 1 and corresponds to a transition from spherical cup (GSI = −1) through a symmetric saddle-shaped LC (GSI = 0) to spherical cap (anteriorly curved LC; GSI = 1). The mean LCD was calculated as the mean depth of all points on the surface of the lamina cribrosa.

**Statistical analyses.** Baseline values were defined as those before trabeculectomy. The primary outcome variables were the postoperatively changes in VF MD and RNFL thickness. Rates of progression were calculated from longitudinal data using linear mixed models adjusting for baseline age, sex, follow-up duration, and scan quality of OCT at each visit and accounting for correlation between eyes. Associations between clinical factors (independent variables) and changes in VF MD and RNFL thickness (dependent variables) were assessed by
linear mixed models to account for longitudinal data. In addition to age, sex baseline IOP, and scan quality of OCT, variables with P < 0.1 from the univariate model were included in the multivariate model. Descriptive statistics were calculated as the mean and standard deviation (SD) and number (%). We used the Wilcoxon signed-rank sum test to compare the preoperative and postoperative measurements of IOP, VF MD and RNFL thickness. A P value of less than 0.05 was considered statistically significant. Analysis was performed using STATA 12.1 (StataCorp LP, College Station, Texas).

## Data availability

The datasets generated during and/or analyzed during the current study are not publicly available due to the terms of consent to which the participants agreed but are available from the corresponding author on reasonable request.

Received: 4 September 2019; Accepted: 4 February 2020;

Published online: 18 February 2020

## References

1. The Advanced Glaucoma Intervention Study (AGIS): 7. The relationship between control of intraocular pressure and visual field deterioration. The AGIS Investigators. *American journal of ophthalmology* **130**, 429–440 (2000).

2. Lichter, P. R. et al. Interim clinical outcomes in the Collaborative Initial Glaucoma Treatment Study comparing initial treatment randomized to medications or surgery. *Ophthalmology* **108**, 1943–1953 (2001).

3. The effectiveness of intraocular pressure reduction in the treatment of normal-tension glaucoma. Collaborative Normal-Tension Glaucoma Study Group. *American journal of ophthalmology* **126**, 498–505 (1998).

4. Heijl, A. et al. Reduction of intraocular pressure and glaucoma progression: results from the Early Manifest Glaucoma Trial. *Archives of ophthalmology* **120**, 1268–1279 (2002).

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

6. Wright, T. M., Goharian, I., Gardiner, S. K., Sehi, M. & Greenfield, D. S. Short-term enhancement of visual field sensitivity in glaucomatous eyes following surgical intraocular pressure reduction. *American journal of ophthalmology* **159**, 378–385 e371, https://doi.org/10.1016/j.ajo.2014.11.012 (2015).

7. Caprioli, J. et al. Trabeculectomy Can Improve Long-Term Visual Function in Glaucoma. *Ophthalmology* **123**, 117–128, https://doi.org/10.1016/j.ophtha.2015.09.027 (2016).

8. Musch, D. C. et al. Visual field improvement in the collaborative initial glaucoma treatment study. *American journal of ophthalmology* **158**, 96–104 e102, https://doi.org/10.1016/j.ajo.2014.04.003 (2014).

9. Spaeth, G. L. The effect of change in intraocular pressure on the natural history of glaucoma: lowering intraocular pressure in glaucoma can result in improvement of visual fields. *Transactions of the ophthalmological societies of the United Kingdom* **104**(Pt 3), 256–264 (1985).

10. Katz, L. J., Spaeth, G. L., Cantor, L. B., Poryzees, E. M. & Steinmann, W. C. Reversible optic disk cupping and visual field improvement in adults with glaucoma. *American journal of ophthalmology* **107**, 485–492 (1989).

11. Kotecha, A. et al. Optic disc and visual field changes after trabeculectomy. *Investigative ophthalmology & visual science* **50**, 4693–4699, https://doi.org/10.1167/iovs.08-3115 (2009).

12. Parrish, R. K. 2nd et al. Five-year follow-up optic disc findings of the Collaborative Initial Glaucoma Treatment Study. *American journal of ophthalmology* **147**, 717–724 e711, https://doi.org/10.1016/j.ajo.2008.10.007 (2009).

13. Figus, M. et al. Short-term changes in the optic nerve head and visual field after trabeculectomy. *Eye* **25**, 1057–1063, https://doi.org/10.1038/eye.2011.119 (2011).

14. Aydin, A., Wollstein, G., Price, L. L., Fujimoto, J. G. & Schuman, J. S. Optical coherence tomography assessment of retinal nerve fiber layer thickness changes after glaucoma surgery. *Ophthalmology* **110**, 1506–1511, https://doi.org/10.1016/S0161-6420(03)00493-7 (2003).

15. Sogano, S., Tomita, G. & Kitazawa, Y. Changes in retinal nerve fiber layer thickness after reduction of intraocular pressure in chronic open-angle glaucoma. *Ophthalmology* **100**, 1253–1258 (1993).

16. Raghu, N., Pandav, S. S., Kaushik, S., Ichhpujani, P. & Gupta, A. Effect of trabeculectomy on RNFL thickness and optic disc parameters using optical coherence tomography. *Eye* (Lond) **26**, 1131–1137, https://doi.org/10.1038/eye.2012.115 (2012).

17. Yamada, N., Tomita, G., Yamamoto, T. & Kitazawa, Y. Changes in the nerve fiber layer thickness following a reduction of intraocular pressure after trabeculectomy. *J Glaucoma* **9**, 371–375 (2000).

18. Membrey, W. L., Bunce, C., Poinosawmy, D. P., Fitzke, F. W. & Hitchings, R. A. Glaucoma surgery with or without adjunctive antiproliferatives in normal tension glaucoma: 2 Visual field progression. *The British journal of ophthalmology* **85**, 696–701 (2001).

19. Musch, D. C. et al. Intraocular pressure control and long-term visual field loss in the Collaborative Initial Glaucoma Treatment Study. *Ophthalmology* **118**, 1766–1773, https://doi.org/10.1016/j.ophtha.2011.01.047 (2011).

20. Folgar, F. A. et al. Effect of successful and partly successful filtering surgery on the velocity of glaucomatous visual field progression. *Journal of glaucoma* **21**, 615–618, https://doi.org/10.1097/IJG.0b013e31821db409 (2012).

21. Baril, C. et al. Rates of glaucomatous visual field change after trabeculectomy. *The British journal of ophthalmology* **101**, 874–878, https://doi.org/10.1136/bjophthalmol-2016-308948 (2017).

22. Naito, T. et al. Effect of trabeculectomy on visual field progression in Japanese progressive normal-tension glaucoma with intraocular pressure <15 mmHg. *PloS one* **12**, e0184096, https://doi.org/10.1371/journal.pone.0184096 (2017).

23. Mwanza, J. C. et al. Retinal nerve fibre layer thickness floor and corresponding functional loss in glaucoma. *The British journal of ophthalmology* **99**, 732–737, https://doi.org/10.1136/bjophthalmol-2014-305745 (2015).

24. Mwanza, J. C. et al. Residual and Dynamic Range of Retinal Nerve Fiber Layer Thickness in Glaucoma: Comparison of Three OCT Platforms. *Investigative ophthalmology & visual science* **56**, 6344–6351, https://doi.org/10.1167/iovs.15-17248 (2015).

25. Miraftabi, A. et al. Macular SD-OCT Outcome Measures: Comparison of Local Structure-Function Relationships and Dynamic Range. *Investigative ophthalmology & visual science* **57**, 4815–4823, https://doi.org/10.1167/iovs.16-19648 (2016).

26. Bowd, C., Zangwill, L. M., Weinreb, R. N., Medeiros, F. A. & Belghith, A. Estimating Optical Coherence Tomography Structural Measurement Floors to Improve Detection of Progression in Advanced Glaucoma. *American journal of ophthalmology* **175**, 37–44, https://doi.org/10.1016/j.ajo.2016.11.010 (2017).

27. Lavinsky, F. et al. Can Macula and Optic Nerve Head Parameters Detect Glaucoma Progression in Eyes with Advanced Circumpapillary Retinal Nerve Fiber Layer Damage? *Ophthalmology* **125**, 1907–1912, https://doi.org/10.1016/j.ophtha.2018.05.020 (2018).

28. Pavlidis, M., Stupp, T., Naskar, R., Cengiz, C. & Thanos, S. Retinal ganglion cells resistant to advanced glaucoma: a postmortem study of human retinas with the carbocyanine dye DiI. *Investigative ophthalmology & visual science* **44**, 5196–5205 (2003).
29. Kadziauskiene, A. et al. Long-Term Shape, Curvature, and Depth Changes of the Lamina Cribrosa after Trabeculectomy. *Ophthalmology* **125**, 1729-1740, https://doi.org/10.1016/j.ophtha.2018.05.011 (2018).

30. Krzyzanowska-Berkowska, P., Melinska, A., Helemejko, I. & Robert Iskander, D. Evaluating displacement of lamina cribrosa following glaucoma surgery. *Graefes Arch Clin Exp Ophthalmol* **256**, 791-800, https://doi.org/10.1007/s00417-018-3920-1 (2018).

31. Krzyzanowska-Berkowska, P., Czajor, K., Helemejko, I. & Iskander, D. R. Relationship between the rate of change in lamina cribrosa depth and the rate of retinal nerve fiber layer thinning following glaucoma surgery. *PLoS One* **13**, e0206040, https://doi.org/10.1371/journal.pone.0206040 (2018).

32. Ha, A. et al. Baseline Lamina Cribrosa Curvature and Subsequent Visual Field Progression Rate in Primary Open-Angle Glaucoma. *Ophthalmology* **125**, 1898-1906, https://doi.org/10.1016/j.ophtha.2018.05.017 (2018).

33. Leung, C. K. Optical Coherence Tomography Imaging for Glaucoma - Today and Tomorrow. *Asia Pac. J Ophthalmol (Phila)* **5**, 11-16, https://doi.org/10.1097/APO.0000000000000179 (2016).

34. Mwanza, J. C., Warren, J. L. & Budenz, D. L. Utility of combining spectral domain optical coherence tomography structural parameters for the diagnosis of early Glaucoma: a mini-review. *Eye and vision* **5**, 9, https://doi.org/10.1186/s40662-018-0101-6 (2018).

35. Susanna, R. Jr. & Vessani, R. M. Staging glaucoma patient: why and how? *Open Ophthalmol J* **3**, 59-64, https://doi.org/10.2174/1874364100903020059 (2009).

36. Girard, M. J., Strouthidis, N. G., Ethier, C. R. & Mari, J. M. Shadow removal and contrast enhancement in optical coherence tomography images of the human optic nerve head. *Investigative ophthalmology & visual science* **52**, 7738-7748, https://doi.org/10.1167/iovs.10-6925 (2011).

37. Girard, M. J. et al. Lamina cribrosa visibility using optical coherence tomography: comparison of devices and effects of image enhancement techniques. *Investigative ophthalmology & visual science* **56**, 865-874, https://doi.org/10.1167/iovs.14-14903 (2015).

38. Mari, J. M., Strouthidis, N. G., Park, S. C. & Girard, M. J. Enhancement of lamina cribrosa visibility in optical coherence tomography images using adaptive compensation. *Investigative ophthalmology & visual science* **54**, 2238-2247, https://doi.org/10.1167/iovs.12-11327 (2013).

39. Thakku, S. G. et al. A Global Shape Index to Characterize Anterior Lamina Cribrosa Morphology and Its Determinants in Healthy Indian Eyes. *Investigative ophthalmology & visual science* **56**, 3604-3614, https://doi.org/10.1167/iovs.15-16707 (2015).

## Acknowledgements

The study is funded by National Medical Research Council (NMRC CG/C010A/2017 and OFIRG/0048/2017 and TA/MOH-000249-00/2018) and the Duke-NUS Medical School (Duke-NUS-KP(Coll)/2018/0009A), Singapore.
The sponsor or funding organization had no role in the design or conduct of this research.

## Author contributions

J.C., A.K. and L.S. conceived and designed the study. A.K., R.A. and E.L. collected the data. J.C., D.W., N.D.Q., R.C., B.T., M.J.A.G., J.M.M., J.G.C., T.A. and L.S. analyzed and interpreted the data. J.C. and L.S. wrote the main manuscript text. All authors reviewed the manuscript.

## Competing interests

The authors declare no competing interests.

## Additional information

Correspondence and requests for materials should be addressed to L.S.

Reprints and permissions information is available at www.nature.com/reprints.

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

<figure><img src="image_6.png" /></figure>

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 license, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons license 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 license, visit http://creativecommons.org/licenses/by/4.0/.

© The Author(s) 2020