Long-Term Remodeling Response in the Lamina Cribrosa Years after IOP Lowering by Suturelysis after Trabeculectomy
=================================================================================================================

* Cameron A. Czerpak
* Harry A. Quigley
* Thao D. Nguyen

## ABSTRACT

**Objective** To measure the active remodeling of the lamina cribrosa (LC) years after IOP lowering by suturelysis.

**Design** Cohort study.

**Participants** Glaucoma patients were imaged 20 minutes after laser suturelysis following trabeculectomy surgery and at their follow-up appointment 1-4 years later (16 image pairs; 15 persons).

**Intervention** Non-invasive optical coherence tomography (OCT) imaging of the eye.

**Main Outcomes** Deformation calculated by correlating OCT scans of the LC immediately after IOP lowering by suturelysis and those acquired years later (defined as *follow-up* strain).

**Results** Mean LC follow-up strain in the anteroposterior direction (*Ezz*) was 14.0 ± 21.3% (mean, standard deviation, p=0.03), while the LC anterior border moved 60.9 ± 54.6 μm into the eye (p=0.0006) on long-term, maintained IOP lowering. *Ezz* at follow-up was 14 times larger than the direct *Ezz* response to IOP lowering by suturelysis. There was a significant association between larger LC anterior movement and greater *Ezz* (p=0.004) at follow-up. Thinner retinal nerve fiber layer (RNFL) at suturelysis was associated with greater follow-up *Ezz* (p=0.04). Worsening visual field indexes during follow-up were associated with greater LC widening (positive remodeling *Eθθ*, p=0.02). Eyes with a greater counterclockwise twist (positive *Eθz*) at suturelysis had greater reversal clockwise twist at follow-up (negative remodeling *Eθz*, p=0.007).

**Conclusion** Follow-up strains and LC border position changes measured years after IOP lowering are far larger than the immediate strain response and LC border movement response to IOP lowering and indicate dramatic remodeling of the LC anatomical structure caused by IOP lowering and glaucoma progression. The remodeling includes a substantial increase in LC thickness and movement into the eye. Eyes with greater direct strain response to IOP-lowering strains, greater glaucoma damage at suturelysis, and greater worsening of visual field at follow-up experienced greater remodeling

**Trial Registration** [ClinicalTrials.gov](https://ClinicalTrials.gov) Identifier: [NCT03267849](http://medrxiv.org/lookup/external-ref?link_type=CLINTRIALGOV&access_num=NCT03267849&atom=%2Fmedrxiv%2Fearly%2F2023%2F06%2F03%2F2023.05.30.23290730.atom)

Key words
*   remodeling
*   strain
*   glaucoma
*   optic nerve head
*   digital volume correlation

## INTRODUCTION

The lamina cribrosa (LC) is the primary site of retinal ganglion cell axon damage in glaucoma within the optic nerve head (ONH).1–4 Increases in intraocular pressure (IOP) from glaucoma increase biomechanical stress at the ONH, producing deformations that can be quantified by measuring strains, which estimate expansions, contractions, and twisting motions. Strains in the ONH can lead to detrimental effects on astrocytes, axons, and nutritional blood flow.5–9 The ONH response to IOP may provide biomarkers for the incidence and progressive worsening of glaucoma.10

The structures and collagen network in the ONH have been shown to remodel in human glaucoma eyes11–13 and experimental monkey glaucoma.14–17 The LC alterations include a thinner, more curved, and more excavated shape.18 The pores within the LC network are filled in by astrocytes that migrate into pores.19 IOP affects the cells within the ONH, and its neurons, glia and fibroblasts have been shown to have mechanosensitive properties.20–24

Glaucoma treatment to reduce IOP has been demonstrated to be beneficial, but the cellular events that occur in the ONH after IOP lowering have not been studied in detail. We and others have conducted *in vivo* research to study the events following IOP lowering using optical coherence tomography (OCT). Sequential OCT imaging can measure ONH mechanical behavior by comparing images at higher and lower IOP. Many previous studies measured only the depth change of the LC anterior border25 after IOP change, demonstrating that LC depth decreases weeks to months after IOP lowering.26–28 Shortly after an average of 12 mmHg of IOP lowering, we found that the LC moved an average of 1.33 ± 6.26 µm posteriorly (out of the eye) when measured 20 minutes after suturelysis following trabeculectomy surgery. The time sequence of LC border movement after trabeculectomy, studied by Kadziauskiene et al.29, was greatest 3-10 days after trabeculectomy (40 μm anterior), followed by a slower, further anterior LC border movement during the next 6 months, as also found by Krzyżanowska-Berkowska et al.30, 31 and others.32 This suggests that the LC border experiences at least 2 phases of change after IOP lowering: an initial period of passive biomechanics, potentially due to viscoelastic response of its extracellular matrix, followed by a longer, slower period of tissue remodeling involving cellular responses.

Thus, the movement of the LC border is a limited surrogate for the dynamic events affecting LC mechanics. Deformations within the LC tissues may provide more accurate assessment of stretching, twisting, and compression that retinal ganglion cell axons experience. We and others have quantified the mechanical deformations from short-term IOP change within the optic nerve by measuring strain using digital volume correlation (DVC).33–42 In this time frame, the LC expanded in thickness by 1% and contracted in radius by 0.2%. By contrast, Girard et al.37 measured LC strains of 8.6% at 50 days after IOP lowering. No study to date measured mechanical ONH strains years after IOP lowering. In this study, we utilized DVC analysis to compare ONH images taken immediately after IOP lowering by suturelysis with images of the same eyes 1-4 years later, producing what are defined here as remodeling strains of the internal LC tissue.

## METHODS

### Experimental Participants

This study was approved by the institutional review board of the Johns Hopkins School of Medicine and can be found at [ClinicalTrials.gov](https://ClinicalTrials.gov) as [NCT03267849](http://medrxiv.org/lookup/external-ref?link_type=CLINTRIALGOV&access_num=NCT03267849&atom=%2Fmedrxiv%2Fearly%2F2023%2F06%2F03%2F2023.05.30.23290730.atom). Written informed consent was obtained from every patient prior to OCT imaging and IOP measurements and the study adheres to the tenets of the Declaration of Helsinki.

Sixteen eyes from 15 patients were reimaged 1-4 years (2.04 ± 0.98 years) after IOP lowering using spectral-domain OCT (SD-OCT, Spectralis, Heidelberg Engineering, Heidelberg, Germany) and the same procedure as used in the prior published study.42 All patients had undergone trabeculectomy and a postoperative suturelysis procedure and participated in our previous study.42 Patient age at the time of suturelysis ranged from 54 to 86 years old (mean 67.6 ± 8.3 years). One patient underwent suturelysis in both the left and right eyes and both were included here. Mean deviation (MD) and visual field index (VFI) were captured at routine visits using Humphrey Field Analyzer (HFA) 24-2 Swedish Interactive Threshold Algorithm (SITA Standard) field tests (Zeiss Meditec, Dublin, CA) within 3 months of suturelysis imaging and the time of follow-up imaging. Retinal nerve fiber layer (RNFL) thickness was measured using Cirrus OCT (Zeiss) at the time of suturelysis. Fewer than half of the patients completed RNFL measurements at the time of the follow-up imaging and therefore change in RNFL over time was not included in the analysis. MD at suturelysis ranged from severe to mild damage (−27.43 to - 1.32 dB; mean (standard deviation) = −10.24 ± 6.82 dB). The distribution of damage was: 4 mild damage eyes (MD > −6), 6 moderate damage eyes (−12 dB < MD < −6 dB), and 6 severe damage eyes (MD < −12 dB). Mean visual field MD change over time from surgery to later imaging was −1.11 ± 3.98 dB.

### OCT Imaging

IOP was measured before each imaging session using the iCare tonometer (iCare Finland Oy). IOP after suturelysis ranged from 3 to 23 mmHg (12.88 ± 6.39 mmHg) and IOP at the 1-4 year follow up ranged from 3 to 19 mmHg (8.75 ± 4.58 mmHg), for a modest average decrease in IOP of 4.13 ± 7.89 mmHg (p=0.05). We also recorded applanation IOP measurements from up to 3 visits before trabeculectomy, with mean = 19.85 ± 7.07 mmHg, range: 7.33 – 37.0 mmHg. In addition, we recorded all IOP measurements from suturelysis to follow-up imaging (minimum, maximum, mean, range, standard deviation, and variance). Corneal curvature was measured using the IOL Master instrument (Zeiss Meditec) and input to Spectralis for each imaging session. Each patient was imaged using the ONH-RC scan pattern consisting of 24 radial SD-OCT scans with enhanced depth imaging on (Figure 1). The scans were centered on the ONH and were 768 x 495 pixels with 7.5° between adjacent scans in the radial (R), anterior-posterior (Z), and theta (Θ) directions, respectively. The resolution in R varied due to keratometry corrections (5.33-6.51 μm/pixel) and Z resolution remained constant at 3.87 μm/pixel. The resolution in 0 was approximately 108 μm/pixel at the BMO, with resolution improving towards the center of the LC. Three imaging sets were taken at the same session to estimate baseline and correlation errors using a previously published method.34, 42 Before processing with DVC, the natural speckle pattern of the OCT images was enhanced using contrast-limited adaptive histogram equalization (CLAHE) in FIJI (ImageJ) (Figure 1). Noise was reduced in areas of low signal using a gamma correction of 1.75 following CLAHE. Images were then imported into MATLAB 2019a (Mathworks) to reconstruct the images into a 3- dimensional matrix of 8-bit intensity values.

![Figure 1:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2023/06/03/2023.05.30.23290730/F1.medium.gif)

[Figure 1:](http://medrxiv.org/content/early/2023/06/03/2023.05.30.23290730/F1)

Figure 1: 
A) OCT image of the ONH after suturelysis, B) 1 year later.

### Image Segmentation

All 24 radial scans were manually segmented by a trained observer (C.A.C.) to mark the anterior border of the LC and Bruch’s membrane opening (BMO) on both sides of the ONH. The anterior LC markings were fit using piecewise linear interpolation, and a second parallel curvilinear border was drawn 250 μm posterior to represent the posterior LC. The LC region was defined as the area between the two LC borders and within BMO delineated by vertical lines.42 250 μm was chosen as the posterior border because few eyes have a visible change in LC signal to clearly mark the posterior border of the LC. The chosen thickness was based on histological measurements of the human LC.43

### Digital Volume Correlation

We have published in detail the DVC method used to measure strain and strain error in the LC.34, 42 DVC was applied to correlate the OCT volumes acquired at the follow-up appointment years after suturelysis and the OCT volumes acquired 20 minutes after suturelysis to calculate the follow-up deformation. In contrast, the direct deformation response to IOP lowering was calculated by correlating the OCT volumes taken 20 minutes after suturelysis and the OCT volume acquired before suturelysis. The deformation at follow-up includes anatomical changes caused by the active remodeling response to glaucoma and the initial IOP lowering and by the direct deformation response to a mild decrease in IOP during the follow-up period.

DVC is a method that tracks the natural speckle pattern of the OCT images in three dimensions (R, Z, Θ) from a before image to an after image to measure mechanical deformations (displacements). The 3-D volume of the first image is broken down into subvolumes, and the algorithm searches the second volume for the best match to the before volume. At every location the algorithm tries to match the volumes, and it measures how well the volumes correlated (correlation-coefficient). The algorithm then records what displacements in R, Z, and Θ were needed to move the subvolume from the before to the after location with the highest correlation-coefficient. The DVC algorithm used in this study is an extension of Bar-Kochba et al.44 to analyze radial scans and calculate displacements *Ur* (radial), *Uz* (axial), and *Uθ* (circumferential) in the cylindrical coordinate system. The algorithm subvolume begins large at 256 x 256 x 32 pixels (R, Z, Θ) and refines to 2 x 2 x 1 pixels to better measure large differences in displacements. A correlation coefficient indicates the quality of the subvolume matching, with low values primarily occurring due to image noise, low signal, or blood vessel shadows. A correlation coefficient of 0.055 was used to remove displacements from dark regions with no contrast.34, 42 The average percent of correlation in the LC was 49.2 ± 16% across the 16 image sets.

The output of DVC is a field of displacements. We then took the gradient of the displacement field to measure the relative deformations (expansions, contractions, and shears) within the field (strain). LC depth change following suturelysis was measured from the displacement of the anterior LC border *Uz* with respect to the 2 BMO points. Displacements were averaged across the entire LC border, rather than choosing the central most point in the LC to measure displacement. Strain and LC depth change outcomes will be referred to as *follow-up*, meaning they measure the remodeling strains or displacements from after suturelysis to the *follow-up* imaging session years later. Strain and LC-depth change measured in our previous study42 with these eyes will be referred to as *after-lysis*, meaning they measure the strain from before suturelysis to 20 minutes *after suturelysis*.

The three normal (*Ezz*, *Err*, *Eθθ*) and three shear strains (*Ezθ*, *Erz*, *Erθ*) of the Green-Lagrange strain tensor were measured using DVC in 3-dimensions: anterior-posterior strain (*Ezz*), radial strain (*Err*), circumferential strain (*Eθθ*), and shear strains cylindrical twist strain (*Ezθ*), in plane shear strain (*Erz*), and shear strain that transforms a circular cross section into a twisted ellipse (*Erθ*). Positive values for normal strains represent elongation, while negative values represent contraction. Shear strains are deformations in two directions. A shear strain applied to a rectangle transforms the shape into a diamond. The maximum principal strain (*Emax*) and maximum shear strain (*Γmax*) were calculated in the RZ plane (within the radial scan) to reduce the impact of the lower resolution of image data between the radial slices.

The potential error in strains were estimate in two ways: baseline error and correlation error. Baseline error is calculated by comparing 2 image sets taken immediately after the other and applying DVC. Ideally, baseline error would be 0; however, image quality and patient movement or misalignment produce some disparity. Correlation error was measured by manually applying a known displacement to an OCT image, then strains are calculated between the original and distorted images, and DVC applied. The difference between strain calculated from the known distortion and that after applying the displacement is the correlation error, again ideally 0. However, low image quality and the presence of blood vessels overlying the LC contribute to correlation error.

For the present images, the average overall baseline and correlation errors (*Uz*) were −0.21 ± 0.73 μm and 0.62 ± 0.12 μm, respectively (Supplement Tables S1, S2). The average baseline and correlation strain errors for strain *Ezz* were 0.134 ± 0.167% and 0.531 ± 0.122%, and for *Err* were 0.021 ± 0.106% and −0.118 ± 0.079%, respectively. The strains and error equations have previously been published.34, 42

One patient had sustained low IOP (mean = 4 mmHg) with stable function despite some chorioretinal wrinkling, which required use of radial scans from two separate follow-up imaging sets to register with the after-lysis scans. Measurement of LC border movement was excluded for this patient as the measurement relies on the BMO position.

### Analysis

We used t-tests to determine the significance of the mean LC follow-up strains compared to no change and to baseline and correlation errors. Relationships between the follow-up strain, IOP, LC border movement, glaucoma damage (MD, VFI, and RNFL thickness), and after-lysis strain were analyzed using linear regression.

## RESULTS

### LC strains and depth change

At follow-up, the LC border moved 60.9 +/- 54.6 μm into the eye (anterior) (p=0.0007) and LC follow-up *Ezz* strain was 14.0 +/- 21% (tensile, or increased in thickness, p=0.02; Figure 1, Table 1). Follow-up *Emax* and *ρmax* were predominately influenced by the large follow-up *Ezz* measurement and were 17.7 +/- 20.3% (p=0.003) and 10.7 +/- 9.74% (p=0.0005), respectively. Only follow-up *Ezz* outcomes were greater than baseline error (p=0.02) and correlation error (p=0.02) (Supplement Table S3). The other follow-up components of the strain tensor were 1 to 2 orders of magnitude smaller than follow-up *Ezz* (Table 1). Follow-up radial strain *Err* was nearly 0 (p=0.98) and the LC contracted circumferentially but not significantly (follow-up *Eθθ* = −0.278 +/- 2.44%; p=0.66). Comparing the follow-up strains to the after-lysis strains published previously for this subset of eyes42 showed that the follow-up *Ezz* was 14 times larger than after-lysis *Ezz* (0=0.02) and the follow-up LC border anterior movement was 124 times larger than 20 minutes after-lysis (p=0.0007). Follow-up *Emax* and *Γmax* were also 8 (p=0.003) and 6 (p=0.0005) times larger than after-lysis *Emax* and *Γmax*.

![Figure 2:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2023/06/03/2023.05.30.23290730/F2.medium.gif)

[Figure 2:](http://medrxiv.org/content/early/2023/06/03/2023.05.30.23290730/F2)

Figure 2: 
A representative nasal-temporal single B-scan from one patient illustrates a strongly positive follow-up *Ezz*, indicating a 14% increase in LC thickness one year after IOP-lowering.

View this table:
[Table 1:](http://medrxiv.org/content/early/2023/06/03/2023.05.30.23290730/T1)

Table 1: 
Comparison of mean strains measured 20 minutes after suturelysis and follow-up strains from 1- 4 years later. ALD = anterior LC depth change (minus indicates anterior movement) (n=16 eyes).

### LC depth change with strain

We compared follow-up strains and follow-up LC border movement, which showed that the LC expands as the LC border moves into the eye (anterior). Greater movement of the LC border into the eye at follow-up was associated with more tensile follow-up *Ezz* (R2=0.48, p=0.004), greater *Emax* (R2=0.51, p=0.003), and greater *Γmax* (R2=0.54, p=0.002) (Figure 3) (Supplement Table S4). Follow-up LC border movement was not associated with follow-up *Err*, *Eθθ*, or shear strains (all p>0.11).

![Figure 3:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2023/06/03/2023.05.30.23290730/F3.medium.gif)

[Figure 3:](http://medrxiv.org/content/early/2023/06/03/2023.05.30.23290730/F3)

Figure 3: 
Linear regression analysis showed that greater anterior movement of the LC border was associated with greater follow-up *Ezz* (LC thickening)

### LC after-lysis and IOP change

A smaller IOP decrease immediately after suturelysis was associated with greater follow-up movement of the LC border in the anterior direction (R2=0.27, p=0.046) (Supplement Table S5). We next compared the mean IOP for visits prior to trabeculectomy to the immediately after-lysis IOP and related this difference to LC outcomes. The IOP decrease from before trabeculectomy to after-lysis was not significantly associated with follow-up *Eθθ* (R2=0.23, p=0.06), *Err* (R2=0.20, p=0.08), or *Ezz* (R2=0.0095, p=0.72) (Supplement Table S6). The IOP change from suturelysis to long-term follow-up was not associated with follow-up strain or LC border movement (p≥0.1) (Supplement Table S7). The time from suturelysis to follow-up imaging (1-4 years) was also not associated with the follow-up strain response (p>0.06) (Supplement Table S8). Regarding the IOP history at all appointments following suturelysis, neither maximum, minimum, range, mean, standard deviation nor variance of IOP were associated with follow-up strain or LC border movement (all p>0.06) (Supplement Tables S9- S14).

### LC after-lysis and glaucomatous neuropathy

We investigated the effects of glaucoma damage on the follow-up strain using the RNFL thickness, MD, and VFI at suturelysis, and the MD and VFI change from after-lysis to long-term follow-up. Eyes with more glaucoma damage (thinner RNFL) experienced greater follow-up *Ezz* (more expansion; R2=0.24, p=0.05; Figure 4), while thinner RNFL was nearly significantly associated with more anterior movement of the LC border (R2=0.24, p=0.06) (Supplement Table S15). MD and VFI at suturelysis were not associated with follow-up strain or LC border movement (p>0.3) (Supplement Tables S16-S17). Patients with greater MD decrease from suturelysis to follow-up experienced more tensile follow-up *Eθθ* (R2=0.31, p=0.02) and positive *Erθ* (R2=0.28, p=0.03) (Supplement Table S18). Patients with greater decrease in VFI from suturelysis to follow up experienced more tensile *Eθθ* (R2=0.33, p=0.02) (Figure 5) (Supplement Table S19).

![Figure 4:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2023/06/03/2023.05.30.23290730/F4.medium.gif)

[Figure 4:](http://medrxiv.org/content/early/2023/06/03/2023.05.30.23290730/F4)

Figure 4: 
Eyes with more glaucoma damage (thinner RNFL) at suturelysis experienced greater follow-up *Ezz* (LC thickening)

![Figure 5:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2023/06/03/2023.05.30.23290730/F5.medium.gif)

[Figure 5:](http://medrxiv.org/content/early/2023/06/03/2023.05.30.23290730/F5)

Figure 5: 
Patients with worse progression (decrease in VFI from suturelysis to follow-up) experienced more tensile follow-up *Eθθ* (widening).

### Comparison of after-lysis and follow-up strains

Some after-lysis strains were associated with follow-up strains. Eyes with a greater counterclockwise twist (positive *Eθz*) after-lysis had greater clockwise twist (negative *Eθz*) at follow-up (R2=0.42, p=0.007) (Figure 6) (Supplement Table S20). Eyes with a greater after-lysis shear strain (*Erθ*) had a borderline association with greater MD decrease during the follow-up period (R2=0.16, p=0.12) (Supplement Table S21).

![Figure 6:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2023/06/03/2023.05.30.23290730/F6.medium.gif)

[Figure 6:](http://medrxiv.org/content/early/2023/06/03/2023.05.30.23290730/F6)

Figure 6: 
Eyes with a greater twist *Eθz* had greater negative follow-up *Eθz*.

## DISCUSSION

While many publications documented the dynamic nature of LC border movement after IOP lowering, only one study, to our knowledge, has measured LC strains at some time after IOP lowering.37 We measured strains both immediately after lowering and years later, with the hypothesis that detailed measurement of LC behavior throughout its structure would be more informative than documentation of its border alone, as has been frequently reported.26, 45 Our data show that follow-up LC strains were up to 14 times larger than after-lysis strains, with the most obvious increase in strain direction being anterior-posterior (*Ezz*), indicating axial expansion. While follow-up radial, circumferential, and shear strains were larger than their after-lysis counterparts, they remained small compared to follow-up *Ezz* strain.

We found that more movement of the LC border in the anterior direction was associated with larger *Emax* and *Γmax*, similar to *Ezz*. It is certainly consistent with intuitive logic that the LC would move anteriorly if it became thicker. However, as we previously reported34, 42, LC border movement is not associated with all components of strain. Furthermore, it is more likely that data on internal LC mechanical change would provide information on the expansions, contractions, and twisting motions transmitted to the axons that are thought to underlie glaucoma injury.

Strains measured 20 minutes after-lysis42 were associated with follow-up strains. Counterclockwise twist after-lysis was predictive of clockwise twisting at follow-up. Eyes with more glaucoma damage exhibited a more compliant strain response after-lysis and eyes with greater glaucoma damage at lysis experienced more axial expansion at follow-up. Together, these results suggest that the strain response and mechanical properties of the LC may provide biomarkers for the medium term remodeling process that could contribute to glaucoma damage. We hypothesize that strains measured over a short period of time (20 minutes after suturelysis or 1 week after IOP lowering by drops medication) could predict which patients will have greater glaucoma progressive worsening. Greater after-lysis shear strain had a borderline significant association with greater MD decrease at follow-up (p=0.12). It may be that a larger study will find that LC biomechanics offer a biomarker for the early diagnosis and progression of glaucoma.

A large follow-up *Ezz* may indicate the LC is returning toward a state prior to damage, given that the LC is thinner and more curved with glaucoma damage.18 However, the eyes with the largest *Ezz* and LC border movement had the thinnest RNFL thickness at the time of suturelysis (p=0.05, p=0.06 respectively). We previously reported that the short-term strain response (after-lysis) was more compliant in eyes with greater glaucoma damage, as measured by thinner RNF, more negative MD, and lower VFI.42 The result may be due to the thinner structure of a glaucoma LC, changes in the material properties of the LC with glaucoma damage, or a combination of both factors. While there has been some evidence that the sclera, overall, is less compliant with experimental glaucoma and in human glaucoma eyes,46–53 we determined that in vitro mechanical responses are more compliant in eyes with greater pre-mortem glaucoma damage.54 Furthermore, recent histological studies found that eyes with greater axon loss had thinner LC connective tissue beams.19 Taken together, these facts indicate that LC compliance is either greater in eyes destined to develop worse damage, or, that the damage process produces greater compliance. Prospective, longitudinal study of the predictive value of LC strains and their change over time are needed to determine which of these possibilities is more likely.

In our previous study, we defined LC compliance as strain per mm Hg IOP lowering, determining that greater compliance was associated with worse glaucoma damage.42 Since the IOP at follow-up was slightly, but not greatly lower than that immediately after suturelysis, we could not generate such a compliance measure for the eyes studied here. We did find that the IOP decrease at suturelysis was associated with greater follow-up LC border movement. However, there is a wide range of IOP change at suturelysis, which depends on a number of short-term factors in the immediate post-operative period. Generally, one might expect a larger mechanical stress (larger IOP change) to cause more movement in the short-term, but how this would influence more remote behavior is problematic.20–24

The main limitations of this study are its relatively small sample size, variable follow-up imaging times, and the variability and slight decrease in IOP at follow-up imaging. Initial patients in this study were imaged in 2018 and the repeat imaging was only undertaken after validation of the method.34, 42 Linear regression of the strain measures compared to time between did not show suturelysis and follow-up imaging trended toward a significant associations (p>0.06), and a larger study with multiple imaging time points is indicated. We plan to address both limitations by quantifying the changes in strain over time in eyes following eye drop treatment at consistent time intervals. In a subsequent report, we will demonstrate that valid strains can be measured one week after either starting or temporarily stopping eye drop treatment.

## Supporting information

Supplemental Tables [[supplements/290730_file02.pdf]](pending:yes)

## Data Availability

All data produced in the present study are available upon reasonable request to the authors

## Footnotes

*   **Conflict of Interest:** The authors have no conflicting relationships related to the devices used in this research. Heidelberg Engineering donated the imaging device, but were not involved with data analysis or interpretation, nor with drafting of the manuscript.

*   **Financial Support:** This research was funded in part by NIH NEI R01 02120 (Dr. Quigley), Brightfocus G2021012S, P30 01765 (Wilmer Institute Core Grant), Research to Prevent Blindness, the A. Edward Maumenee Professorship (Dr. Quigley), and BrightFocus.

*   **Supplemental Information:** This article contains additional online-only material. The following should appear online-only: Supplement Tables S1-S21

## Abbreviations

LC
:   lamina cribrosa
ONH
:   optic nerve head
PPS
:   peripapillary sclera
SD-OCT
:   spectral domain optical coherence tomography
RNFL
:   retinal nerve fiber layer
MD
:   mean deviation
VFI
:   visual field index
IOP
:   intraocular pressure.

*   Received May 30, 2023.
*   Revision received May 30, 2023.
*   Accepted June 3, 2023.


*   © 2023, Posted by Cold Spring Harbor Laboratory

The copyright holder for this pre-print is the author. All rights reserved. The material may not be redistributed, re-used or adapted without the author's permission.

## REFERENCES

1.  1.Anderson DR, Hendrickson A. Effect of intraocular pressure on rapid axoplasmic transport in monkey optic nerve. Invest Ophthalmol Vis Sci. 1974;13(10).
    
    

2.  2.Vrabec F. Glaucomatous cupping of the human optic disk - A neuro-histologic study. Albr von Graefes Arch für Klin und Exp Ophthalmol. 1976;198(3). doi:10.1007/BF00410715
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1007/BF00410715&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=1083166&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F06%2F03%2F2023.05.30.23290730.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=A1976BK42800004&link_type=ISI) 

3.  3.Quigley HA, Addicks EM, Green WR, Maumenee AE. Optic nerve damage in human glaucoma. II. The site of injury and susceptibility to damage. Arch Ophthalmol (Chicago, Ill 1960). 1981;99(4):635-649.
    
    

4.  4.Quigley HA, Addicks EM. Regional Differences in the Structure of the Lamina Cribrosa and Their Relation to Glaucomatous Optic Nerve Damage. Arch Ophthalmol. 1981;99(1). doi:10.1001/archopht.1981.03930010139020
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1001/archopht.1981.03930010139020&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=7458737&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F06%2F03%2F2023.05.30.23290730.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=A1981KZ84900019&link_type=ISI) 

5.  5.Minckler DS, Bunt AH, Johanson GW. Orthograde and retrograde axoplasmic transport during acute ocular hypertension in the monkey. Invest Ophthalmol Vis Sci. 1977;16(5).
    
    

6.  6.Gaasterland D, Tanishima T, Kuwabara T. Axoplasmic flow during chronic experimental glaucoma. I. Light and electron microscopic studies of the monkey optic nervehead during development of glaucomatous cupping. Invest Ophthalmol Vis Sci. 1978;17(9).
    
    

7.  7.Burgoyne CF, Crawford Downs J, Bellezza AJ, Francis Suh JK, Hart RT. The optic nerve head as a biomechanical structure: A new paradigm for understanding the role of IOP-related stress and strain in the pathophysiology of glaucomatous optic nerve head damage. Prog Retin Eye Res. 2005;24(1). doi:10.1016/j.preteyeres.2004.06.001
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.preteyeres.2004.06.001&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=15555526&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F06%2F03%2F2023.05.30.23290730.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000225852400002&link_type=ISI) 

8.  8.Sigal IA, Ethier CR. Biomechanics of the optic nerve head. Exp Eye Res. 2009;88:799–807.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.exer.2009.02.003&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=19217902&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F06%2F03%2F2023.05.30.23290730.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000265778000025&link_type=ISI) 

9.  9.Quigley HA. The contribution of the sclera and lamina cribrosa to the pathogenesis of glaucoma: Diagnostic and treatment implications. In: Progress in Brain Research. Vol 220.; 2015. doi:10.1016/bs.pbr.2015.04.003
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/bs.pbr.2015.04.003&link_type=DOI) 

10. 10.Quigley HA, Cone FE. Development of diagnostic and treatment strategies for glaucoma through understanding and modification of scleral and lamina cribrosa connective tissue. Cell Tissue Res. 2013;353(2). doi:10.1007/s00441-013-1603-0
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1007/s00441-013-1603-0&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=23535950&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F06%2F03%2F2023.05.30.23290730.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000321951900004&link_type=ISI) 

11. 11.Quigley HA, Dorman-Pease ME, Brown AE. Quantitative study of collagen and elastin of the optic nerve head and sclera in human and experimental monkey glaucoma. Curr Eye Res. 1991;10(9). doi:10.3109/02713689109013884
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.3109/02713689109013884&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=1790718&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F06%2F03%2F2023.05.30.23290730.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=A1991GG04300010&link_type=ISI) 

12. 12.Quigley H, Pease ME, Thibault D. Change in the appearance of elastin in the lamina cribrosa of glaucomatous optic nerve heads. Graefes Arch Clin Exp Ophthalmol. 1994;232(5):257–261. doi:10.1007/BF00194473
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1007/BF00194473&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=8045433&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F06%2F03%2F2023.05.30.23290730.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=A1994NH43100001&link_type=ISI) 

13. 13.Jones HJ, Girard MJ, White N, et al. Quantitative analysis of three-dimensional fibrillar collagen microstructure within the normal, aged and glaucomatous human optic nerve head. J R Soc Interface. 2015;12(106). doi:10.1098/rsif.2015.0066
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1098/rsif.2015.0066&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=25808336&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F06%2F03%2F2023.05.30.23290730.atom) 

14. 14.Yang H, Williams G, Crawford Downs J, et al. Posterior (outward) migration of the lamina cribrosa and early cupping in monkey experimental glaucoma. Invest Ophthalmol Vis Sci. 2011;52(10). doi:10.1167/iovs.11-7448
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NDoiaW92cyI7czo1OiJyZXNpZCI7czoxMDoiNTIvMTAvNzEwOSI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDIzLzA2LzAzLzIwMjMuMDUuMzAuMjMyOTA3MzAuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 

15. 15.He L, Yang H, Gardiner SK, et al. Longitudinal detection of optic nerve head changes by spectral domain optical coherence tomography in early experimental glaucoma. Invest Ophthalmol Vis Sci. 2013;55(1). doi:10.1167/iovs.13-13245
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NDoiaW92cyI7czo1OiJyZXNpZCI7czo4OiI1NS8xLzU3NCI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDIzLzA2LzAzLzIwMjMuMDUuMzAuMjMyOTA3MzAuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 

16. 16.Ivers KM, Sredar N, Patel NB, et al. In vivo changes in lamina cribrosa microarchitecture and optic nerve head structure in early experimental glaucoma. PLoS One. 2015;10(7). doi:10.1371/journal.pone.0134223
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.pone.0134223&link_type=DOI) 

17. 17.Strouthidis NG, Fortune B, Yang H, Sigal IA, Burgoyne CF. Longitudinal change detected by spectral domain optical coherence tomography in the optic nerve head and peripapillary retina in experimental glaucoma. Invest Ophthalmol Vis Sci. 2011;52(3). doi:10.1167/iovs.10-5599
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NDoiaW92cyI7czo1OiJyZXNpZCI7czo5OiI1Mi8zLzEyMDYiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyMy8wNi8wMy8yMDIzLjA1LjMwLjIzMjkwNzMwLmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 

18. 18.Quigley HA, Varma R, Tielsch JM, Katz J, Sommer A, Gilbert DL. The relationship between optic disc area and open-angle glaucoma: The Baltimore eye survey. J Glaucoma. 1999;8(6).
    
    

19. 19.Guan C, Pease ME, Quillen S, et al. Quantitative Microstructural Analysis of Cellular and Tissue Remodeling in Human Glaucoma Optic Nerve Head. Investig Ophthalmol Vis Sci. 2022;63(11). doi:10.1167/IOVS.63.11.18
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1167/IOVS.63.11.18&link_type=DOI) 

20. 20.Choi HJ, Sun D, Jakobs TC. Astrocytes in the optic nerve head express putative mechanosensitive channels. Mol Vis. 2015;21.
    
    

21. 21.Rogers R, Dharsee M, Ackloo S, Flanagan JG. Proteomics analyses of activated human optic nerve head lamina cribrosa cells following biomechanical strain. Invest Ophthalmol Vis Sci. 2012;53(7). doi:10.1167/iovs.11-8480
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NDoiaW92cyI7czo1OiJyZXNpZCI7czo5OiI1My83LzM4MDYiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyMy8wNi8wMy8yMDIzLjA1LjMwLjIzMjkwNzMwLmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 

22. 22.Quillen S, Schaub J, Quigley H, Pease M, Korneva A, Kimball E. Astrocyte responses to experimental glaucoma in mouse optic nerve head. PLoS One. 2020;15(8 August 2020). doi:10.1371/journal.pone.0238104
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.pone.0238104&link_type=DOI) 

23. 23.Xia J, Lim JC, Lu W, et al. Neurons respond directly to mechanical deformation with pannexin-mediated ATP release and autostimulation of P2X 7 receptors. J Physiol. 2012;590(10). doi:10.1113/jphysiol.2012.227983
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1113/jphysiol.2012.227983&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=22411013&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F06%2F03%2F2023.05.30.23290730.atom) 

24. 24.Kirwan RP, Fenerty CH, Crean J, Wordinger RJ, Clark AF, O’Brien CJ. Influence of cyclical mechanical strain on extracellular matrix gene expression in human lamina cribrosa cells in vitro. Mol Vis. 2005;11. doi:10.1016/s0021-9290(06)84552-5
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/s0021-9290(06)84552-5&link_type=DOI) 

25. 25.Rhodes LA, Huisingh C, Johnstone J, et al. Variation of laminar depth in normal eyes with age and race. Invest Ophthalmol Vis Sci. 2014;55(12). doi:10.1167/iovs.14-15251
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NDoiaW92cyI7czo1OiJyZXNpZCI7czoxMDoiNTUvMTIvODEyMyI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDIzLzA2LzAzLzIwMjMuMDUuMzAuMjMyOTA3MzAuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 

26. 26.Lee EJ, Kim TW, Weinreb RN. Reversal of lamina cribrosa displacement and thickness after trabeculectomy in glaucoma. Ophthalmology. 2012;119(7). doi:10.1016/j.ophtha.2012.01.034
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.ophtha.2012.01.034&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=22464141&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F06%2F03%2F2023.05.30.23290730.atom) 

27. 27.Park HYL, Shin HY, Jung KI, Park CK. Changes in the lamina and prelamina after intraocular pressure reduction in patients with primary open-angle glaucoma and acute primary angle-closure. Invest Ophthalmol Vis Sci. 2014;55(1). doi:10.1167/iovs.12-10329
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NDoiaW92cyI7czo1OiJyZXNpZCI7czo4OiI1NS8xLzIzMyI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDIzLzA2LzAzLzIwMjMuMDUuMzAuMjMyOTA3MzAuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 

28. 28.Yoshikawa M, Akagi T, Hangai M, et al. Alterations in the neural and connective tissue components of glaucomatous cupping after glaucoma surgery using swept-source optical coherence tomography. Invest Ophthalmol Vis Sci. 2014;55(1). doi:10.1167/iovs.13-11897
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NDoiaW92cyI7czo1OiJyZXNpZCI7czo4OiI1NS8xLzQ3NyI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDIzLzA2LzAzLzIwMjMuMDUuMzAuMjMyOTA3MzAuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 

29. 29.Kadziauskienė A, Jašinskienė E, Ašoklis R, et al. Long-Term Shape, Curvature, and Depth Changes of the Lamina Cribrosa after Trabeculectomy. Ophthalmology. 2018;125(11). doi:10.1016/j.ophtha.2018.05.011
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.ophtha.2018.05.011&link_type=DOI) 

30. 30.Krzyżanowska-Berkowska P, Melińska A, Helemejko I, Robert Iskander D. Evaluating displacement of lamina cribrosa following glaucoma surgery. Graefe’s Arch Clin Exp Ophthalmol. 2018;256(4). doi:10.1007/s00417-018-3920-1
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1007/s00417-018-3920-1&link_type=DOI) 

31. 31.Iskander R, Krzyżanowska-Berkowska P, Czajor K, Helemejko I. Relationship between the rate of change in lamina cribrosa depth and the rate of retinal nerve fiber layer thinning following glaucoma surgery. PLoS One. 2018;13(11). doi:10.1371/journal.pone.0206040
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.pone.0206040&link_type=DOI) 

32. 32.Park DY, Son S, Cha SC. Comparison of lamina cribrosa depth shallowing after trabeculectomy between primary open-angle glaucoma and exfoliation glaucoma. Sci Rep. 2022;12(1). doi:10.1038/s41598-022-19785-2
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41598-022-19785-2&link_type=DOI) 

33. 33.Girard MJA, Strouthidis NG, Desjardins A, Mari JM, Ethier CR. In vivo optic nerve head biomechanics: Performance testing of a threedimensional tracking algorithm. J R Soc Interface. 2013;10(87). doi:10.1098/rsif.2013.0459
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1098/rsif.2013.0459&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=23883953&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F06%2F03%2F2023.05.30.23290730.atom) 

34. 34.Midgett DE, Quigley HA, Nguyen TD. In vivo characterization of the deformation of the human optic nerve head using optical coherence tomography and digital volume correlation. Acta Biomater. 2019;96:385–399. doi:10.1016/j.actbio.2019.06.050
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.actbio.2019.06.050&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F06%2F03%2F2023.05.30.23290730.atom) 

35. 35.Beotra MR, Wang X, Tun TA, et al. In vivo three-dimensional lamina cribrosa strains in healthy, ocular hypertensive, and glaucoma eyes following acute intraocular pressure elevation. Invest Ophthalmol Vis Sci. 2018;59(1). doi:10.1167/iovs.17-21982
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1167/iovs.17-21982&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F06%2F03%2F2023.05.30.23290730.atom) 

36. 36.Zhang L, Beotra MR, Baskaran M, et al. In vivo measurements of prelamina and lamina cribrosa biomechanical properties in humans. Invest Ophthalmol Vis Sci. 2020;61(3). doi:10.1167/iovs.61.3.27
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1167/iovs.61.3.27&link_type=DOI) 

37. 37.Girard MJA, Beotra MR, Chin KS, et al. In Vivo 3-Dimensional Strain Mapping of the Optic Nerve Head Following Intraocular Pressure Lowering by Trabeculectomy. Ophthalmology. 2016;123(6). doi:10.1016/j.ophtha.2016.02.008
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.ophtha.2016.02.008&link_type=DOI) 

38. 38.Wang X, Beotra MR, Tun TA, et al. In vivo 3-dimensional strain mapping confirms large optic nerve head deformations following horizontal eye movements. Invest Ophthalmol Vis Sci. 2016;57(13). doi:10.1167/iovs.16-20560
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1167/iovs.16-20560&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F06%2F03%2F2023.05.30.23290730.atom) 

39. 39.Sigal IA, Wei J, Yang B, et al. Measuring in-vivo and in-situ ex-vivo the 3D deformation of the lamina cribrosa microstructure under elevated intraocular pressure. In:; 2018. doi:10.1117/12.2291243
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1117/12.2291243&link_type=DOI) 

40. 40.Tran H, Grimm J, Wang B, et al. Mapping in-vivo optic nerve head strains caused by intraocular and intracranial pressures. In: Optical Elastography and Tissue Biomechanics IV. Vol 10067.; 2017. doi:10.1117/12.2257360
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1117/12.2257360&link_type=DOI) 

41. 41.Fazio MA, Clark ME, Bruno L, Girkin CA. In vivo optic nerve head mechanical response to intraocular and cerebrospinal fluid pressure: imaging protocol and quantification method. Sci Rep. 2018;8(1). doi:10.1038/s41598-018-31052-x
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41598-018-31052-x&link_type=DOI) 

42. 42.Czerpak CA, Kashaf MS, Zimmerman BK, Quigley HA, Nguyen TD. The Strain Response to Intraocular Pressure Decrease in the Lamina Cribrosa of Patients with Glaucoma. Ophthalmol Glaucoma. 2022. doi:[https://doi.org/10.1016/j.ogla.2022.07.005](https://doi.org/10.1016/j.ogla.2022.07.005)
    
    

43. 43.Quigley HA, Hohman RM, Addicks EM, Massof RW, Richard Green W. Morphologic changes in the lamina cribrosa correlated with neural loss in open-angle glaucoma. Am J Ophthalmol. 1983;95(5). doi:10.1016/0002-9394(83)90389-6
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/0002-9394(83)90389-6&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=6846459&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F06%2F03%2F2023.05.30.23290730.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=A1983QR12600015&link_type=ISI) 

44. 44.Bar-Kochba E, Toyjanova J, Andrews E, Kim KS, Franck C. A Fast Iterative Digital Volume Correlation Algorithm for Large Deformations. Exp Mech. 2015;55(1). doi:10.1007/s11340-014-9874-2
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1007/s11340-014-9874-2&link_type=DOI) 

45. 45.Lanzagorta-Aresti A, Perez-Lopez M, Palacios-Pozo E, Davo-Cabrera J. Relationship between corneal hysteresis and lamina cribrosa displacement after medical reduction of intraocular pressure. Br J Ophthalmol. 2017;101(3). doi:10.1136/bjophthalmol-2015-307428
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MTI6ImJqb3BodGhhbG1vbCI7czo1OiJyZXNpZCI7czo5OiIxMDEvMy8yOTAiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyMy8wNi8wMy8yMDIzLjA1LjMwLjIzMjkwNzMwLmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 

46. 46.Hommer A, Fuchsjäger-Mayrl G, Resch H, Vass C, Garhofer G, Schmetterer L. Estimation of ocular rigidity based on measurement of pulse amplitude using pneumotonometry and fundus pulse using laser interferometry in glaucoma. Invest Ophthalmol Vis Sci. 2008;49(9). doi:10.1167/iovs.07-1342
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NDoiaW92cyI7czo1OiJyZXNpZCI7czo5OiI0OS85LzQwNDYiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyMy8wNi8wMy8yMDIzLjA1LjMwLjIzMjkwNzMwLmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 

47. 47.Wang J, Freeman EE, Descovich D, et al. Estimation of ocular rigidity in glaucoma using ocular pulse amplitude and pulsatile choroidal blood flow. Invest Ophthalmol Vis Sci. 2013;54(3). doi:10.1167/iovs.12-9841
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NDoiaW92cyI7czo1OiJyZXNpZCI7czo5OiI1NC8zLzE3MDYiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyMy8wNi8wMy8yMDIzLjA1LjMwLjIzMjkwNzMwLmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 

48. 48.Xiao C, Huang X, Li M, Wu J, Huang H. Evaluation of the biomechanical properties of the posterior eye using Shear wave Elastography in patients with increased intracranial pressure. Med Ultrason. 2021;23(3). doi:10.11152/mu-2940
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.11152/mu-2940&link_type=DOI) 

49. 49.Zeimer RC, Ogura Y. The Relation Between Glaucomatous Damage and Optic Nerve Head Mechanical Compliance. Arch Ophthalmol. 1989;107(8). doi:10.1001/archopht.1989.01070020298042
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1001/archopht.1989.01070020298042&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=2527025&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F06%2F03%2F2023.05.30.23290730.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=A1989AK11700034&link_type=ISI) 

50. 50.Coudrillier B, Tian J, Alexander S, Myers KM, Quigley HA, Nguyen TD. Biomechanics of the human posterior sclera: Age-and glaucoma-related changes measured using inflation testing. Invest Ophthalmol Vis Sci. 2012;53(4):1714–1728. doi:10.1167/iovs.11-8009
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NDoiaW92cyI7czo1OiJyZXNpZCI7czo5OiI1My80LzE3MTQiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyMy8wNi8wMy8yMDIzLjA1LjMwLjIzMjkwNzMwLmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 

51. 51.Sigal IA, Grimm JL, Jan NJ, Reid K, Minckler DS, Brown DJ. Eye-specific iop-induced displacements and deformations of human lamina cribrosa. Invest Ophthalmol Vis Sci. 2014;55(1). doi:10.1167/iovs.13-12724
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NDoiaW92cyI7czo1OiJyZXNpZCI7czo2OiI1NS8xLzEiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyMy8wNi8wMy8yMDIzLjA1LjMwLjIzMjkwNzMwLmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 

52. 52.Behkam R, Kollech HG, Jana A, et al. Racioethnic differences in the biomechanical response of the lamina cribrosa. Acta Biomater. 2019;88. doi:10.1016/j.actbio.2019.02.028
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.actbio.2019.02.028&link_type=DOI) 

53. 53.Midgett DE, Jefferys JL, Quigley HA, Nguyen TD. The inflation response of the human lamina cribrosa and sclera: Analysis of deformation and interaction. Acta Biomater. 2020;106. doi:10.1016/j.actbio.2020.01.049
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.actbio.2020.01.049&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F06%2F03%2F2023.05.30.23290730.atom) 

54. 54.Midgett D, Liu B, Ling YTT, Jefferys JL, Quigley HA, Nguyen TD. The effects of glaucoma on the pressure-induced strain response of the human lamina cribrosa. Invest Ophthalmol Vis Sci. 2020;61(4). doi:10.1167/iovs.61.4.41
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1167/iovs.61.4.41&link_type=DOI)