Regional Deformation of the Optic Nerve Head and Peripapillary Sclera During IOP Elevation

Exploring the Feasibility of Estimating Intraocular Pressure Using Vibrational Response of the Eye: A Methodological Approach

This study addresses the limitations of current tonometry techniques by exploring vibroacoustic properties for estimating intraocular pressure (IOP), a key diagnostic parameter for monitoring glaucoma-a significant risk factor for vision loss. Utilizing vivo porcine eyeballs, we investigated the relationship between IOP and the nonlinear vibration transfer function ratio (NVTFR). Through applying varying vibration levels and analyzing responses with transfer function analysis and univariate regression, we identified a strong negative correlation between NVTFR and IOP, evidenced by a Pearson correlation coefficient of -0.8111 and significant results from generalized linear model (GLM) regression (p-value < 0.001). These findings indicate the potential of NVTFR as a vital indicator of IOP changes. Our study highlights the feasibility of using vibroacoustic properties, specifically NVTFR, to measure IOP. While further refinement is necessary for in vivo application, this approach opens new possibilities for non-invasive and patient-friendly IOP monitoring, potentially enhancing ophthalmology diagnostic techniques and providing a foundation for future research and development in this critical area.

Biomechanical Correlations Between the Cornea and the Optic Nerve Head

Purpose

A thin cornea is a potent risk factor for glaucoma. The underlying mechanisms remain unexplained. It has been postulated that central corneal thickness (CCT) may be a surrogate for biomechanical parameters of the posterior eye. In this study, we aimed to explore correlations of biomechanical responses between the cornea and the optic nerve head (ONH) and the peripapillary sclera (PPS) to elevated intraocular pressure (IOP), the primary risk factor of glaucoma.

Methods

Inflation tests were performed in nine pairs of human donor globes. One eye of each pair was randomly assigned for cornea or posterior eye inflation. IOP was raised from 5 to 30 millimeters of mercury (mmHg) at 0.5 mmHg steps in the whole globe and the cornea or the ONH/PPS was imaged using a 50 MHz ultrasound probe. Correlation-based ultrasound speckle tracking was used to calculate tissue displacements and strains. Associations of radial, tangential, and shear strains at 30 mmHg between the cornea and the ONH or PPS were evaluated.

Results

Corneal shear strain was significantly correlated with ONH shear strain (R = 0.857, P = 0.003) and PPS shear strain (R = 0.724, P = 0.028). CCT was not correlated with any strains in the cornea, ONH, or PPS.

Conclusions

Our results suggested that an eye that experiences a larger shear strain in the cornea would likely experience a larger shear strain in its ONH and PPS at IOP elevations. The strong correlation between the cornea's and the ONH's shear response to IOP provides new insights and suggests a plausible explanation of the cornea's connection to glaucoma risk.

Fibrous finite element modeling of the optic nerve head region

The optic nerve head (ONH) region at the posterior pole of the eye is supported by a fibrous structure of collagen fiber bundles. Discerning how the fibrous structure determines the region biomechanics is crucial to understand normal physiology, and the roles of biomechanics on vision loss. The fiber bundles within the ONH structure exhibit complex three-dimensional (3D) organization and continuity across the various tissue components. Computational models of the ONH, however, usually represent collagen fibers in a homogenized fashion without accounting for their continuity across tissues, fibers interacting with each other and other fiber-specific effects in a fibrous structure. We present a fibrous finite element (FFE) model of the ONH that incorporates discrete collagen fiber bundles and their histology-based 3D organization to study ONH biomechanics as a fibrous structure. The FFE model was constructed using polarized light microscopy data of porcine ONH cryosections, representing individual fiber bundles in the sclera, dura and pia maters with beam elements and canal tissues as continuum structures. The FFE model mimics the histological in-plane orientation and width distributions of collagen bundles as well as their continuity across different tissues. Modeling the fiber bundles as linear materials, the FFE model predicts the nonlinear ONH response observed in an inflation experiment from the literature. The model also captures important microstructural mechanisms including fiber interactions and long-range strain transmission among bundles that have not been considered before. The FFE model presented here advances our understanding of the role of fibrous collagen structure in the ONH biomechanics. STATEMENT OF SIGNIFICANCE: The microstructure and mechanics of the optic nerve head (ONH) are central to ocular physiology. Histologically, the ONH region exhibits a complex continuous fibrous structure of collagen bundles. Understanding the role of the fibrous collagen structure on ONH biomechanics requires high-fidelity computational models previously unavailable. We present a computational model of the ONH that incorporates histology-based fibrous collagen structure derived from polarized light microscopy images. The model predictions agree with experiments in the literature, and provide insight into important microstructural mechanisms of fibrous tissue biomechanics, such as long-range strain transmission along fiber bundles. Our model can be used to study the microstructural basis of biomechanical damage and the effects of collagen remodeling in glaucoma.

Three-Dimensional Ultrasound Elastography Detects Age-Related Increase in Anterior Peripapillary Sclera and Optic Nerve Head Compression During IOP Elevation

Purpose

High-frequency ultrasound elastography offers a tool to resolve the complex and heterogeneous deformation through the full thickness of the optic nerve head (ONH) and peripapillary sclera (PPS). Using this tool, we quantified the three-dimensional deformation of the ONH and PPS in human donor eyes and evaluated age-associated changes.

Methods

The ONH and PPS in 15 human donor globes were imaged with a 50-MHz ultrasound probe while increasing intraocular pressure (IOP) from 15 to 30 mm Hg. Tissue displacements were obtained using correlation-based ultrasound speckle tracking. Three-dimensional spherical strains (radial, circumferential, meridional, and respective shear strains) were calculated for the ONH and PPS volumes segmented from three-dimensional ultrasound images. Age-related trends of different strains in each region of interest were explored.

Results

The dominant form of IOP-induced deformation in the ONH and PPS was radial compression. High-magnitude localized out-of-plane shear strains were also observed in both regions. Most strains were concentrated in the anterior one-half of the ONH and PPS. The magnitude of radial and volumetric strains increased with age in the anterior ONH and anterior PPS, indicating greater radial compression and volume loss during IOP elevation in older age.

Conclusions

The age-associated increase of radial compression, the predominant form of IOP-induced deformation in anterior ONH and PPS, may underlie age-associated glaucoma risk. High-frequency ultrasound elastography offers a useful tool to quantify all types of deformation comprehensively in all regions of ONH and PPS, which may improve our understanding of the biomechanical factors contributing to glaucoma risk.

An Inverse Method to Determine Mechanical Parameters of Porcine Vitreous Bodies Based on the Indentation Test

The vitreous body keeps the lens and retina in place and protects these tissues from physical insults. Existing studies have reported that the mechanical properties of vitreous body varied after liquefaction, suggesting mechanical properties could be effective parameters to identify vitreous liquefaction process. Thus, in this work, we aimed to propose a method to determine the mechanical properties of vitreous bodies. Fresh porcine eyes were divided into three groups, including the untreated group, the 24 h liquefaction group and the 48 h liquefaction group, which was injected collagenase and then kept for 24 h or 48 h. The indentation tests were carried out on the vitreous body in its natural location while the posterior segment of the eye was fixed in the container. A finite element model of a specimen undertaking indentation was constructed to simulate the indentation test with surface tension of vitreous body considered. Using the inverse method, the mechanical parameters of the vitreous body and the surface tension coefficient were determined. For the same parameter, values were highest in the untreated group, followed by the 24 h liquefaction group and the lowest in the 48 h liquefaction group. For C₁₀ in the neo-Hookean model, the significant differences were found between the untreated group and liquefaction groups. This work quantified vitreous body mechanical properties successfully using inverse method, which provides a new method for identifying vitreous liquefactions related studies.

Simultaneous Assessment of the Whole Eye Biomechanics Using Ultrasonic Elastography

OBJECTIVE

Current elastography techniques in the field of ophthalmology usually target one specific tissue, such as the cornea or the sclera. However, the eye is an inter-related organ, and some ocular diseases can alter the biomechanical properties of multiple anatomical structures. Hence, there is a need to develop an imaging tool that can non-invasively, quantitatively, and accurately characterize dynamic changes among these biomechanical properties.

METHODS

A high resolution ultrasound elastography system was developed to achieve this goal. The efficacy and accuracy of the system was first validated on tissue-mimicking phantoms while mechanical testing measurements served as the gold standard. Next, an in vivo elevated intraocular pressure (IOP) model was established in rabbits to further test our system. In particular, elastography measurements were obtained at 5 IOP levels, ranging from 10 mmHg to 30 mmHg in 5 mmHg increments. Spatial-temporal maps of the multiple ocular tissues (cornea, lens, iris, optic nerve head, and peripapillary sclera) were obtained.

RESULTS

The spatial-temporal maps were acquired simultaneously for the ocular tissues at the 5 different IOP levels. The statistical analysis of the elastic wave speed was presented for ocular tissues. Finally, the mapping for the elastic wave speed of each ocular component was acquired at each IOP level.

CONCLUSION

Our elastography system can concurrently assess the biomechanical properties of multiple ocular structures and detect changes in biomechanical properties associated with changes in IOP.

SIGNIFICANCE

This system provides a novel tool to measure and quantify the biomechanical properties of the whole eye.

The impact of horizontal eye movements versus intraocular pressure on optic nerve head biomechanics: A tridimensional finite element analysis study

It has been proposed that eye movements could be related to glaucoma development. This research aimed to compare the impact of intraocular pressure (IOP) versus horizontal duction on optic nerve head (ONH) strains. Thus, a tridimensional finite element model of the eye including the three tunics of the eye, all of the meninges, and the subarachnoid space (SAS) was developed using a series of medical tests and anatomical data. The ONH was divided into 22 subregions, and the model was subjected to 21 different eye pressures, as well as 24 different degrees of adduction and abduction ranging from 0.5° to 12°. Mean deformations were documented along anatomical axes and in principal directions. Additionally, the impact of tissue stiffness was assessed. The results show no statistically significant differences between the lamina cribrosa (LC) strains due to eye rotation and IOP variation. However, when assessing LC regions some experienced a reduction in principal strains following a 12° duction, while after the IOP reached 12 mmHg, all LC subzones showed an increase in strains. From an anatomical perspective, the effect on the ONH following 12° duction was opposite to that observed after a rise in IOP. Moreover, high strain dispersion inside the ONH subregions was obtained with lateral eye movements, which was not observed with increased IOP and variation. Finally, SAS and orbital fat stiffness strongly influenced ONH strains during eye movements, while SAS stiffness was also influential under ocular hypertension. Even if horizontal eye movements cause large ONH deformations, their biomechanical effect would be markedly distinct from that induced by IOP. It could be predicted that, at least in physiological conditions, their potential to cause axonal injury would not be so relevant. Thus, a causative role in glaucoma does not appear likely. By contrast, an important role of SAS would be expectable.

Regional differences of the sclera in the ocular hypertensive rat model induced by circumlimbal suture

PURPOSE

To describe the regional differences of the sclera in ocular hypertension (OHT) models with the inappropriate extension of the ocular axis.

METHODS

To discover the regional differences of the sclera at the early stage, OHT models were established using circumlimbal suture (CS) or sclerosant injection (SI). Axial length (AL) was measured by ultrasound and magnetic resonance imaging. The glaucoma-associated distinction was determined by intraocular pressure (IOP) and retrograde tracing of retinal ganglion cells (RGCs). The central thickness of the ganglion cell complex (GCC) was measured by optical coherence tomography. RGCs and collagen fibrils were detected using a transmission electron microscope, furthermore, anti-alpha smooth muscle actin (αSMA) was determined in the early stage after the operation.

RESULTS

Compared with the control group, the eyes in OHT models showed an increased IOP (P < 0.001 in the CS group, P = 0.001 in the SI group), growing AL (P = 0.026 in the CS group, P = 0.043 in the SI group), reduction of central RGCs (P < 0.001 in the CS group, P = 0.017 in the SI group), thinning central GCC (P < 0.001 in the CS group), and a distinctive expression of αSMA in the central sclera in the early 4-week stage after the operation (P = 0.002 in the CS group). Compared with the SI group, the eye in the CS group showed a significantly increased AL (7.1 ± 0.4 mm, P = 0.031), reduction of central RGCs (2121.1 ± 87.2 cells/mm², P = 0.001), thinning central GCC (71.4 ± 0.8 pixels, P = 0.015), and a distinctive expression of αSMA (P = 0.005). Additionally, ultrastructural changes in RGCs, scleral collagen fibers, and collagen crimp were observed in the different regions. Increased collagen volume fraction in the posterior segment of the eyeball wall (30.2 ± 3.1%, P = 0.022) was observed by MASSON staining in the CS group.

CONCLUSION

Regional differences of the sclera in the ocular hypertensive rat model induced by CS may provide a reference for further treatment of scleral-related eye disorders.

Real-time imaging of optic nerve head collagen microstructure and biomechanics using instant polarized light microscopy

Current tools lack the temporal or spatial resolution necessary to image many important aspects of the architecture and dynamics of the optic nerve head (ONH). We evaluated the potential of instant polarized light microscopy (IPOL) to overcome these limitations by leveraging the ability to capture collagen fiber orientation and density in a single image. Coronal sections through the ONH of fresh normal sheep eyes were imaged using IPOL while they were stretched using custom uniaxial or biaxial micro-stretch devices. IPOL allows identifying ONH collagen architectural details, such as fiber interweaving and crimp, and has high temporal resolution, limited only by the frame rate of the camera. Local collagen fiber orientations and deformations were quantified using color analysis and image tracking techniques. We quantified stretch-induced collagen uncrimping of lamina cribrosa (LC) and peripapillary sclera (PPS), and changes in LC pore size (area) and shape (convexity and aspect ratio). The simultaneous high spatial and temporal resolutions of IPOL revealed complex ONH biomechanics: i) stretch-induced local deformation of the PPS was nonlinear and nonaffine. ii) under load the crimped collagen fibers in the PPS and LC straightened, without torsion and with only small rotations. iii) stretch-induced LC pore deformation was anisotropic and heterogeneous among pores. Overall, with stretch the pores were became larger, more convex, and more circular. We have demonstrated that IPOL reveals details of collagen morphology and mechanics under dynamic loading previously out of reach. IPOL can detect stretch-induced collagen uncrimping and other elements of the tissue nonlinear mechanical behavior. IPOL showed changes in pore morphology and collagen architecture that will help improve understanding of how LC tissue responds to load.

Interplay between intraocular and intracranial pressure effects on the optic nerve head in vivo

Intracranial pressure (ICP) has been proposed to play an important role in the sensitivity to intraocular pressure (IOP) and susceptibility to glaucoma. However, the in vivo effects of simultaneous, controlled, acute variations in ICP and IOP have not been directly measured. We quantified the deformations of the anterior lamina cribrosa (ALC) and scleral canal at Bruch's membrane opening (BMO) under acute elevation of IOP and/or ICP. Four eyes of three adult monkeys were imaged in vivo with OCT under four pressure conditions: IOP and ICP either at baseline or elevated. The BMO and ALC were reconstructed from manual delineations. From these, we determined canal area at the BMO (BMO area), BMO aspect ratio and planarity, and ALC median depth relative to the BMO plane. To better account for the pressure effects on the imaging, we also measured ALC visibility as a percent of the BMO area. Further, ALC depths were analyzed only in regions where the ALC was visible in all pressure conditions. Bootstrap sampling was used to obtain mean estimates and confidence intervals, which were then used to test for significant effects of IOP and ICP, independently and in interaction. Response to pressure manipulation was highly individualized between eyes, with significant changes detected in a majority of the parameters. Significant interactions between ICP and IOP occurred in all measures, except ALC visibility. On average, ICP elevation expanded BMO area by 0.17 mm² at baseline IOP, and contracted BMO area by 0.02 mm² at high IOP. ICP elevation decreased ALC depth by 10 μm at baseline IOP, but increased depth by 7 μm at high IOP. ALC visibility decreased as ICP increased, both at baseline (-10%) and high IOP (-17%). IOP elevation expanded BMO area by 0.04 mm² at baseline ICP, and contracted BMO area by 0.09 mm² at high ICP. On average, IOP elevation caused the ALC to displace 3.3 μm anteriorly at baseline ICP, and 22 μm posteriorly at high ICP. ALC visibility improved as IOP increased, both at baseline (5%) and high ICP (8%). In summary, changing IOP or ICP significantly deformed both the scleral canal and the lamina of the monkey ONH, regardless of the other pressure level. There were significant interactions between the effects of IOP and those of ICP on LC depth, BMO area, aspect ratio and planarity. On most eyes, elevating both pressures by the same amount did not cancel out the effects. Altogether our results show that ICP affects sensitivity to IOP, and thus that it can potentially also affect susceptibility to glaucoma.

Strain by virtual extensometers and video-imaging optical coherence tomography as a repeatable metric for IOP-Induced optic nerve head deformations

PURPOSE

To determine if in vivo strain response of the Optic Nerve Head (ONH) to IOP elevation visualized using Optical Coherence Tomography (OCT) video imaging and quantified using novel virtual extensometers was able to be provided repeatable measurements of tissue specific deformations.

METHODS

The ONHs of 5 eyes from 5 non-human primates (NHPs) were imaged by Spectralis OCT. A vertical and a horizontal B-scan of the ONH were continuously recorded for 60 s at 6 Hz (video imaging mode) during IOP elevation from 10 to 30 mmHg. Imaging was repeated over three imaging sessions. The 2D normal strain was computed by template-matching digital image correlation using virtual extensometers. ANOVA F-test (F) was used to compare inter-eye, inter-session, and inter-tissue variability for the prelaminar, Bruch's membrane opening (BMO), lamina cribrosa (LC) and choroidal regions (against variance the error term). F-test of the ratio between inter-eye to inter-session variability was used to test for strain repeatability across imaging sessions (FIS).

RESULTS

Variability of strain across imaging session (F = 0.7263, p = 0.4855) and scan orientation was not significant (F = 1.053, p = 0.3066). Inter session variability of strain was significantly lower than inter-eye variability (FIS = 22.63, p = 0.0428) and inter-tissue variability (FIS = 99.33 p = 0.00998). After IOP elevation, strain was highest in the choroid (-18.11%, p < 0.001), followed by prelaminar tissue (-11.0%, p < 0.001), LC (-3.79%, p < 0.001), and relative change in BMO diameter (-0.57%, p = 0.704).

CONCLUSIONS

Virtual extensometers applied to video-OCT were sensitive to the eye-specific and tissue-specific mechanical response of the ONH to IOP and were repeatable across imaging sessions.

Biomechanics of the optic nerve head and sclera in canine glaucoma: A brief review

Glaucoma is a leading cause of irreversible blindness, a progressive optic neuropathy with retinal ganglion cell (RGC) death beginning in the optic nerve head (ONH). A primary risk factor for developing glaucoma is elevated intraocular pressure (IOP). Reducing IOP is the only treatment proven to be effective at delaying disease progression. Nevertheless, even when patients have their IOP reduced, the majority of them continue to lose vision. There are, in both humans and dogs, significant interindividual variabilities in susceptibilities to IOP-induced optic nerve damage. Vision loss progresses much more slowly in Beagles with open-angle glaucoma (OAG) caused by ADAMTS10 mutation. This can be attributed to the mutation-related altered ocular biomechanical properties. The principal site of optic nerve (ON) damage in glaucoma is the ONH. It is suggested that the biomechanical properties of the ONH and the surrounding peripapillary sclera (PPS) contribute to glaucoma development and progression. As far as the beneficial biomechanical properties of the ONH and PPS for a decreased susceptibility and slow progression of glaucoma, data are inconsistent and conflicting. Recent biomechanical studies on beagles with ADAMTS10 mutation demonstrated that the mutant dogs have mechanically weak posterior sclera. This weakness was associated with a reduced collagen density and a lower proportion of insoluble collagen. These changes, observed before glaucoma development, were considered intrinsic characteristics caused by the mutation rather than a secondary effect of IOP elevation. Further studies of ADAMTS10-OAG may elucidate the effects of altered biomechanical properties of ONH and PPS in determining the glaucoma progression.

Non-invasive Clinical Measurement of Ocular Rigidity and Comparison to Biomechanical and Morphological Parameters in Glaucomatous and Healthy Subjects

Purpose: To assess ocular rigidity using dynamic optical coherence tomography (OCT) videos in glaucomatous and healthy subjects, and to evaluate how ocular rigidity correlates with biomechanical and morphological characteristics of the human eye.

Methods: Ocular rigidity was calculated using Friedenwald's empirical equation which estimates the change in intraocular pressure (IOP) produced by volumetric changes of the eye due to choroidal pulsations with each heartbeat. High-speed OCT video was utilized to non-invasively measure changes in choroidal volume through time-series analysis. A control-case study design was based on 23 healthy controls and 6 glaucoma cases. Multiple diagnostic modalities were performed during the same visit including Spectralis OCT for nerve head video, Pascal Dynamic Contour Tonometry for IOP and ocular pulse amplitude (OPA) measurement, Corvis ST for measuring dynamic biomechanical response, and Pentacam for morphological characterization.

Results: Combining glaucoma and healthy cohorts (n = 29), there were negative correlations between ocular rigidity and axial length (Pearson R = −0.53, p = 0.003), and between ocular rigidity and anterior chamber volume (R = −0.64, p = 0.0002). There was a stronger positive correlation of ocular rigidity and scleral stiffness (i.e., stiffness parameter at the highest concavity [SP-HC]) (R = 0.62, p = 0.0005) compared to ocular rigidity and corneal stiffness (i.e., stiffness parameter at the first applanation [SP-A1]) (R = 0.41, p = 0.033). In addition, there was a positive correlation between ocular rigidity and the static pressure-volume ratio (P/V ratio) (R = 0.72, p < 0.0001).

Conclusions: Ocular rigidity was non-invasively assessed using OCT video and OPA in a clinic setting. The significant correlation of ocular rigidity with biomechanical parameters, SP-HC and P/V ratio, demonstrated the validity of the ocular rigidity measurement. Ocular rigidity is driven to a greater extent by scleral stiffness than corneal stiffness. These in vivo methods offer an important approach to investigate the role of ocular biomechanics in glaucoma.

Effect of different preconditioning protocols on the viscoelastic inflation response of the posterior sclera

Preconditioning by repeated cyclic loads is routinely used in ex vivo mechanical testing of soft biological tissues. The goal of preconditioning is to achieve a steady and repeatable mechanical response and to measure material properties that are representative of the in vivo condition. Preconditioning protocols vary across studies, and their effect on the viscoelastic response of tested soft tissue is typically not reported or analyzed. We propose a methodology to systematically analyze the preconditioning process with application to inflation testing. We investigated the effect of preconditioning on the viscoelastic inflation response of tree shrew posterior sclera using two preconditioning protocols: (i) continuous cyclic loading-unloading without rest and (ii) cyclic loading-unloading with 15-min rest between cycles. Posterior scleral surface strain was measured using three-dimensional Digital Image Correlation (3D-DIC). We used five variables of characterizing features of the stress-strain loop curve to compare the two preconditioning protocols. Our results showed protocol-dependent differences in the tissue response during preconditioning and at the preconditioned state. Incorporating a resting time between preconditioning cycles significantly decreased the number of cycles (10.5 ± 2.9 cycles vs. 3.1 ± 0.5 cycles, p < 0.001) but increased the total time (15.8 ± 4.4 min vs. 51.2 ± 8.3 min, p < 0.001) needed to reach preconditioned state. At the preconditioned state, 2 of 5 characteristic variables differed significantly between protocols: hysteresis loop area (difference=0.023 kJ/m³, p = 0.0020) and elastic modulus at high IOPs (difference=24.0 MPa, p = 0.0238). Our results suggest that the analysis of the preconditioning process is an essential part of inflation experiments and a prerequisite to properly characterize the tissue viscoelastic response. Furthermore, material properties obtained at the preconditioned state can be impacted by the resting time used during preconditioning and may not be directly compared across studies if the resting time varies by 15 min between studies. STATEMENT OF SIGNIFICANCE: Although applying a preconditioning protocol by repeated cyclic loads is common practice in ex vivo mechanical characterization of soft tissues, the tissue response is typically not reported or analyzed, and the protocol's potential effect on the response remains unclear. This is partially caused by lack of a standardized methodology to precondition soft tissues. We present the first systematic analysis of two representative preconditioning protocols used during inflation testing in ocular biomechanics. Our results show protocol-dependent differences in the viscoelastic response during the preconditioning process and at the preconditioned state. Consequently, the analysis of the preconditioning response represents an essential part of mechanical testing and a prerequisite to properly characterize the tissue viscoelastic response. The effect of preconditioning on the preconditioned state response must be considered when comparing results across studies with different preconditioning protocols.

Ultrasonic elastography to assess biomechanical properties of the optic nerve head and peripapillary sclera of the eye

PURPOSE

To quantitatively investigate both optic nerve head (ONH) and peripapillary sclera (PPS) biomechanical properties of porcine eyes through an ultrasonic elastography imaging system in response to both increasing and decreasing intraocular pressure (IOP).

METHODS

The Young's modulus of the ONH and PPS were assessed using our high resolution ultrasonic imaging system which utilized a mechanical shaker to induce shear waves and an off-axis aligned 40 MHz needle transducer to track micron-level displacement along the direction of wave propagation. In this study, imaging on a total of 8 ex vivo porcine eyes preloaded with IOPs from 6 mmHg to 30 mmHg was performed. To have a better understanding of the effect of varying IOP on biomechanics, both increasing and decreasing IOPs were investigated.

RESULTS

The increase of the Young's modulus of ONH (92.4 ± 13.9 kPa at 6 mmHg to 224.7 ± 71.1 kPa at 30 mmHg) and PPS (176.8 ± 14.3 kPa at 6 mmHg to 573.5 ± 64.4 kPa at 30 mmHg) following IOP elevation could be observed in the reconstructed Young's modulus of the shear wave elasticity (SWE) imaging while the B-mode structural images remained almost unchanged. In addition, for the same IOP level, both ONH and PPS have a tendency to be stiffer with decreasing IOP as compared to increasing IOP.

CONCLUSIONS

Our results demonstrate the feasibility of using our ultrasonic elastography system to investigate the stiffness mapping of posterior eye with high resolution in both increasing and decreasing IOPs, making this technique potentially useful for glaucoma.

A pilot study for intraocular pressure measurements based on vibroacoustic parameters

The present study aimed to identify vibroacoustic properties associated with intraocular pressure (IOP) changes and to suggest a new way to measure the IOP based on these properties. Ten ex vivo porcine eyeballs were used in this study. Each eyeball was fixated in a central hole of a Styrofoam block, and vibration applied to the Styrofoam block was transmitted to the eyeball. An accelerometer directly attached to the eyeball measured the vibration response. Excitations and measurements were performed for 1 s, and the excitation magnitude was varied for the same signal in repeat measurements. A 30-gauge needle was inserted into the anterior chamber of the eyeball to inject a balanced salt solution, and the height of the bottle was adjusted to adjust the IOP. A tonometer was used under identical conditions to measure the IOP five times, and the mean value was determined for further analyses. The measurements showed that the parameters resonance frequency and change in the magnitude of the vibration response (CMVR) increased with rising IOP values. The CMVR was highly correlated with the IOP (p-value < 0.0001). A linear mixed effects model (LMM) was used as a statistical analysis method. We confirmed that vibroacoustic properties of the eyeball are correlated with IOP changes. It is expected that the CMVR will serve as a new parameter for IOP measurements. Thus, in the future, continuous IOP measurements would be easily performed using the CMVR.

IOP-induced regional displacements in the optic nerve head and correlation with peripapillary sclera thickness

Mechanical insult induced by intraocular pressure (IOP) is likely a driving force in the disease process of glaucoma. This study aimed to evaluate regional displacements in human optic nerve head (ONH) and peripapillary tissue (PPT) in response to acute IOP elevations, and their correlations with morphological characteristics of the posterior eye. Cross-sectional (2D) images of the ONH and PPT in 14 globes of 14 human donors were acquired with high-frequency ultrasound during whole globe inflation from 5 to 30 mm Hg. High-frequency ultrasound has a spatial resolution of tens of micrometers and is capable of imaging through the ONH and PPT thickness. Tissue displacements were calculated using a correlation-based speckle tracking algorithm for a dense matrix of kernels covering the 2D imaging plane. The ONH was manually segmented in the ultrasound B-mode images acquired at 5 mmHg based on echogenicity. The lamina cribrosa (LC) boundaries were visible in eight of the fourteen eyes and the LC region was segmented using a semi-automated superpixel-based method. The ONH had larger radial displacement than the PPT in all tested eyes and the difference increased with increasing IOP. A significant negative correlation was found between ONH-PPT displacement difference and PPT thickness (p < 0.05), while no significant correlations were found between ONH-PPT displacement difference and other morphological parameters including PPT radius of curvature, scleral canal size, LC thickness and anterior LC surface depth. Within the ONH, the radial displacement decreased in the region anterior to and across LC but not in the region posterior to LC. Finite element models using simplified geometry and material properties confirmed the role of LC in reducing the overall ONH radial displacements, but did not predict the displacement gradient change observed experimentally. These results suggested that a thinner PPT may be associated with a larger relative posterior motion of the ONH with respect to the surrounding PPT and the LC may play a major role in preventing excessive posterior displacement of ONH during acute IOP elevations.

Collagen fiber interweaving is central to sclera stiffness

The mechanical properties of the microstructural components of sclera are central to eye physiology and pathology. Because these parameters are extremely difficult to measure directly, they are often estimated using inverse-modeling matching deformations of macroscopic samples measured experimentally. Although studies of sclera microstructure show collagen fiber interweaving, current models do not account for this interweaving or the resulting fiber-fiber interactions, which might affect parameter estimates. Our goal was to test the hypothesis that constitutive parameters estimated using inverse modeling differ if models account for fiber interweaving and interactions. We developed models with non-interweaving or interweaving fibers over a wide range of volume fractions (36-91%). For each model, we estimated fiber stiffness using inverse modeling matching biaxial experimental data of human sclera. We found that interweaving increased the estimated fiber stiffness. When the collagen volume fraction was 64% or less, the stiffness of interweaving fibers was about 1.25 times that of non-interweaving fibers. For higher volume fractions, the ratio increased substantially, reaching 1.88 for a collagen volume fraction of 91%. Simulating a model (interweaving/non-interweaving) using the fiber stiffness estimated from the other model produced substantially different behavior, far from that observed experimentally. These results show that estimating microstructural component mechanical properties is highly sensitive to the assumed interwoven/non-interwoven architecture. Moreover, the results suggest that interweaving plays an important role in determining the structural stiffness of sclera, and potentially of other soft tissues in which the collagen fibers interweave. STATEMENT OF SIGNIFICANCE: The collagen fibers of sclera are interwoven, but numerical models do not account for this interweaving or the resulting fiber-fiber interactions. To determine if interweaving matters, we examined the differences in the constitutive model parameters estimated using inverse modeling between models with interweaving and non-interweaving fibers. We found that the estimated stiffness of the interweaving fibers was up to 1.88 times that of non-interweaving fibers, and that the estimate increased with collagen volume fraction. Our results suggest that fiber interweaving is a fundamental characteristic of connective tissues, additional to anisotropy, density and orientation. Better characterization of interweaving, and of its mechanical effects is likely central to understanding microstructure and biomechanics of sclera and other soft tissues.

Effect of Changing Heart Rate on the Ocular Pulse and Dynamic Biomechanical Behavior of the Optic Nerve Head

Purpose

To study the effect of changing heart rate on the ocular pulse and the dynamic biomechanical behavior of the optic nerve head (ONH) using a comprehensive mathematical model.

Methods

In a finite element model of a healthy eye, a biphasic choroid consisted of a solid phase with connective tissues and a fluid phase with blood, and the lamina cribrosa (LC) was viscoelastic as characterized by a stress-relaxation test. We applied arterial pressures at 18 ocular entry sites (posterior ciliary arteries), and venous pressures at four exit sites (vortex veins). In the model, the heart rate was varied from 60 to 120 bpm (increment: 20 bpm). We assessed the ocular pulse amplitude (OPA), pulse volume, ONH deformations, and the dynamic modulus of the LC at different heart rates.

Results

With an increasing heart rate, the OPA decreased by 0.04 mm Hg for every 10 bpm increase in heart rate. The ocular pulse volume decreased linearly by 0.13 µL for every 10 bpm increase in heart rate. The storage modulus and the loss modulus of the LC increased by 0.014 and 0.04 MPa, respectively, for every 10 bpm increase in heart rate.

Conclusions

In our model, the OPA, pulse volume, and ONH deformations decreased with an increasing heart rate, whereas the LC became stiffer. The effects of blood pressure/heart rate changes on ONH stiffening may be of interest for glaucoma pathology.

Role of radially aligned scleral collagen fibers in optic nerve head biomechanics

Collagen fibers organized circumferentially around the canal in the peripapillary sclera are thought to provide biomechanical support to the sensitive tissues within the optic nerve head (ONH). Recent studies have demonstrated the existence of a family of fibers in the innermost sclera organized radially from the scleral canal. Our goal was to determine the role of these radial fibers in the sensitivity of scleral canal biomechanics to acute increases in intraocular pressure (IOP). Following the same general approach of previous parametric sensitivity studies, we created nonlinear generic finite element models of a posterior pole with various combinations of radial and circumferential fibers at an IOP of 0 mmHg. We then simulated the effects of normal and elevated IOP levels (15 and 30 mmHg). We monitored four IOP-induced geometric changes: peripapillary sclera stretch, scleral canal displacement, lamina cribrosa displacement, and scleral canal expansion. In addition, we examined the radial (maximum tension) and through-thickness (maximum compression) strains within the ONH tissues. Our models predicted that: 1) radial fibers reduced the posterior displacement of the lamina, especially at elevated IOP; 2) radial fibers reduced IOP-induced radial strain within the peripapillary sclera and retinal tissue; and 3) a combination of radial and circumferential fibers maintained strains within the ONH at a level similar to those conferred by circumferential fibers alone. In conclusion, radial fibers provide support for the posterior globe, additional to that provided by circumferential fibers. Most importantly, a combination of both fiber families can better protect ONH tissues from excessive IOP-induced deformation than either alone.

Three-Dimensional Inflation Response of Porcine Optic Nerve Head Using High-Frequency Ultrasound Elastography

Characterization of the biomechanical behavior of the optic nerve head (ONH) in response to intraocular pressure (IOP) elevation is important for understanding glaucoma susceptibility. In this study, we aimed to develop and validate a three-dimensional (3D) ultrasound elastographic technique to obtain mapping and visualization of the 3D distributive displacements and strains of the ONH and surrounding peripapillary tissue (PPT) during whole globe inflation from 15 to 30 mmHg. 3D scans of the posterior eye around the ONH were acquired through full tissue thickness with a high-frequency ultrasound system (50 MHz). A 3D cross-correlation-based speckle-tracking algorithm was used to compute tissue displacements at ∼30,000 kernels distributed within the region of interest (ROI), and the components of the strain tensors were calculated at each kernel by using least square estimation of the displacement gradients. The accuracy of displacement calculation was evaluated using simulated rigid-body translation on ultrasound radiofrequency (RF) data obtained from a porcine posterior eye. The accuracy of strain calculation was evaluated using finite element (FE) models. Three porcine eyes were tested showing that ONH deformation was heterogeneous with localized high strains. Substantial radial (i.e., through-thickness) compression was observed in the anterior ONH and out-of-plane (i.e., perpendicular to the surface of the shell) shear was shown to concentrate in the vicinity of ONH/PPT border. These preliminary results demonstrated the feasibility of this technique to achieve comprehensive 3D evaluation of the mechanical responses of the posterior eye, which may provide mechanistic insights into the regional susceptibility in glaucoma.

The inflation response of the human lamina cribrosa and sclera: Analysis of deformation and interaction

This study investigated the inflation response of the lamina cribrosa (LC) and adjacent peripapillary sclera (PPS) in post-mortem human eyes with no history of glaucoma. The posterior sclera of 13 human eyes from 7 donors was subjected to controlled pressurization between 5-45 mmHg. A laser-scanning microscope (LSM) was used to image the second harmonic generation (SHG) response of collagen and the two-photon fluorescent (TPF) response of elastin within the volume of the LC and PPS at each pressure. Image volumes were analyzed using digital volume correlation (DVC) to calculate the three-dimensional (3D) deformation field between pressures. The LC exhibited larger radial strain, Err, and maximum principal strain, Emax, (p < 0.0001) and greater posterior displacement (p=0.0007) compared to the PPS between 5-45 mmHg, but had similar average circumferential strain, Eθθ, and maximum shear strain, Γmax. The Emax and Γmax were highest near the LC-PPS interface and lowest in the nasal quadrant of both tissues. Larger LC area was associated with smaller Emax in the peripheral LC and larger Emax in the central LC (p ≤ 0.01). The Emax, Γmax, and Eθθ in the inner PPS increased with increasing strain in adjacent LC regions (p ≤ 0.001). Smaller strains in the PPS were associated with a larger difference in the posterior displacement between the PPS and central LC (p < 0.0001 for Emax and Err), indicating that a stiffer pressure-strain response of the PPS is associated with greater posterior bowing of the LC. STATEMENT OF SIGNIFICANCE: Glaucoma causes vision loss through progressive damage of the retinal ganglion axons at the lamina cribrosa (LC), a connective tissue structure that supports the axons as they pass through the eye wall. It is hypothesized that strains caused by intraocular pressure may initiate this damage and that these strains are modulated by the combined deformation of the LC and adjacent peripapillary sclera (PPS). In this study we present a method to measure the pressure-induced 3D displacement and strain field in the LC and PPS simultaneously. Regional strain variation in the LC and PPS was investigated and compared and strains were analyzed for associations with age, LC area, LC strain magnitude, and LC posterior motion relative to the PPS.

Continuous intraocular pressure monitoring in patients with obstructive sleep apnea syndrome using a contact lens sensor

Purpose

To analyse nocturnal intraocular pressure (IOP) fluctuations in patients with obstructive sleep apnea syndrome (OSAS) using a contact lens sensor (CLS) and to identify associations between the OSAS parameters determined by polysomnographic study (PSG) and IOP changes.

Method

Prospective, observational study. Twenty participants suspected of having OSAS were recruited. During PSG study, IOP was monitored using a CLS placed in the eye of the patient. The patients were classified according to the apnea-hypopnea index (AHI) in two categories, severe (>30) or mild/moderate (<30) OSAS. We evaluated several parameters determined by the IOP curves, including nocturnal elevations (acrophase) and plateau times in acrophase (PTs) defined by mathematical and visual methods.

Results

The IOP curves exhibited a nocturnal acrophase followed by PTs of varying extents at which the IOP remained higher than daytime measurement with small variations. We found significant differences in the length of the PTs in patients with severe OSAS compared to those with mild/moderate disease (P = 0.032/P = 0.028). We found a positive correlation between PTs and OSAS severity measured by the total number of apneic events (r = 0.681/0.751 P = 0.004/0.001) and AHI (r = 0.674/0.710, P = 0.004/0.002). Respiratory-related arousal and oxygen saturation also were associated significantly with the IOP PT length.

Conclusions

Periods of nocturnal IOP elevation lasted longer in severe OSAS patients than those with mild/moderate OSAS and correlate with the severity of the disease. The length of the nocturnal PT is also associated to respiratory parameters altered in patients with OSAS.

Quantitative confocal optical coherence elastography for evaluating biomechanics of optic nerve head using Lamb wave model

The mechanosensitivity of the optic nerve head (ONH) plays a pivotal role in the pathogenesis of glaucoma. Characterizing elasticity of the ONH over changing physiological pressure may provide a better understanding of how changes in intraocular pressure (IOP) lead to changes in the mechanical environment of the ONH. Optical coherence elastography (OCE) is an emerging technique that can detect tissue biomechanics noninvasively with both high temporal and spatial resolution compared with conventional ultrasonic elastography. We describe a confocal OCE system in measuring ONH elasticity in vitro, utilizing a pressure inflation setup in which IOP is controlled precisely. We further utilize the Lamb wave model to fit the phase dispersion curve during data postprocessing. We present a reconstruction of Young's modulus of the ONH by combining our OCE system with a Lamb wave model for the first time. This approach enables the quantification of Young's modulus of the ONH, which can be fit using a piecewise polynomial to the corresponding IOP.

Lamina cribrosa pore movement during acute intraocular pressure rise

Purpose

To assess changes in the position of lamina cribrosa pores (LCPs) induced by acute intraocular pressure (IOP) elevation.

Methods

A prospective observational study. Acute angle-closure suspects who underwent the 2-hour dark room prone provocative test (DRPPT) were included. At baseline and within 5 min after the DRPPT end, tonometry, fundus photography and optical coherence tomography were performed. Optic disc photos taken before and after the DRPPT were aligned and moving distance of each visible LCP was measured (LCPMD).

Results

38 eyes from 27 participants (age: 52.5±10.8 years) were included. The IOP rose from 16.7±3.2 mm Hg at baseline to 23.9±4.3 mm Hg at the DRPPT end. The mean lateral LCPMD was 28.1±14.6 µm (range: 5.0–77.2 µm), which increased with higher IOP rise (p=0.01) and deeper optic cup (p=0.02) in multivariate analysis. The intralamina range and SD of the LCPMD increased with younger age (p=0.01 and p=0.02, respectively) and with wider optic cup (p=0.01 and p=0.02, respectively). The LCP movements were headed to the superior direction in 12 (33%) eyes, inferior direction in 10 (28%) eyes, temporal direction in 9 (25%) eyes, and nasal direction in 5 (14%) eyes.

Conclusions

IOP rise is associated with LCP movements in the frontal plane, which are more pronounced with higher IOP rise and deeper optic cup. The intralamina variability in the IOP rise-associated LCPMD increased with younger age and wider optic cup. IOP variation-associated lateral LCP movements may be of interest to elucidate glaucomatous optic nerve damage.

Mechanical Deformation of Human Optic Nerve Head and Peripapillary Tissue in Response to Acute IOP Elevation

Purpose

To measure the deformation of the human optic nerve head (ONH) and peripapillary tissue (PPT) in response to acute intraocular pressure (IOP) elevation.

Methods

The ONH and PPT of 14 human donor globes were imaged with high-frequency ultrasonography during inflation testing from 5 to 30 mm Hg. A correlation-based speckle tracking algorithm was used to compute tissue displacements, and the through-thickness, in-plane, and shear strains were calculated by using least-squares strain estimation methods. The ONH and PPT were segmented along the anterior-posterior direction and the nasal-temporal direction. Regional displacements and strains were analyzed and compared.

Results

The ONH displaced more posteriorly than the PPT in response to an acute IOP increase. Scleral canal expansion was minimal but correlated with ONH posterior displacement at all IOP levels. Through-thickness compression was concentrated in the anterior of both the ONH and the PPT. Shear was concentrated in the vicinity of the canal with higher shear in the peripheral ONH than the central ONH and higher shear in the PPT near the scleral canal than that further away from the canal.

Conclusions

High-resolution ultrasound speckle tracking showed a displacement mismatch between the ONH and the PPT, larger compressive strains in the direction of IOP loading in the anterior ONH and PPT, and higher shear strains in the periphery of ONH in response to acute IOP elevation in the human eye. These findings delineate the deformation patterns within and around the ONH and may help understand IOP-associated optic nerve damage.

Scleral structure and biomechanics

As the eye's main load-bearing connective tissue, the sclera is centrally important to vision. In addition to cooperatively maintaining refractive status with the cornea, the sclera must also provide stable mechanical support to vulnerable internal ocular structures such as the retina and optic nerve head. Moreover, it must achieve this under complex, dynamic loading conditions imposed by eye movements and fluid pressures. Recent years have seen significant advances in our knowledge of scleral biomechanics, its modulation with ageing and disease, and their relationship to the hierarchical structure of the collagen-rich scleral extracellular matrix (ECM) and its resident cells. This review focuses on notable recent structural and biomechanical studies, setting their findings in the context of the wider scleral literature. It reviews recent progress in the development of scattering and bioimaging methods to resolve scleral ECM structure at multiple scales. In vivo and ex vivo experimental methods to characterise scleral biomechanics are explored, along with computational techniques that combine structural and biomechanical data to simulate ocular behaviour and extract tissue material properties. Studies into alterations of scleral structure and biomechanics in myopia and glaucoma are presented, and their results reconciled with associated findings on changes in the ageing eye. Finally, new developments in scleral surgery and emerging minimally invasive therapies are highlighted that could offer new hope in the fight against escalating scleral-related vision disorder worldwide.

Modeling the Mechanical Parameters of Glaucoma

Glaucoma is a major eye disease characterized by a progressive loss of retinal ganglion cells (RGCs). Biomechanical forces as a result of hydrostatic pressure and strain play a role in this disease. Decreasing intraocular pressure is the only available therapy so far, but is not always effective and does not prevent blindness in many cases. There is a need for drugs that protect RGCs from dying in glaucoma; to develop these, we need valid glaucoma and drug screening models. Since in vivo models are unsuitable for screening purposes, we focus on in vitro and ex vivo models in this review. Many groups have studied pressure and strain model systems to mimic glaucoma, to investigate the molecular and cellular events leading to mechanically induced RGC death. Therefore, the focus of this review is on the different mechanical model systems used to mimic the biomechanical forces in glaucoma. Most models use either cell or tissue strain, or fluid- or gas-controlled hydrostatic pressure application and apply it to the relevant cell types such as trabecular meshwork cells, optic nerve head astrocytes, and RGCs, but also to entire eyes. New model systems are warranted to study concepts and test experimental compounds for the development of new drugs to protect vision in glaucoma patients. Impact Statement The outcome of currently developed models to investigate mechanically induced retinal ganglion cell death by applying different mechanical strains varies widely. This suggests that a robust glaucoma model has not been developed yet. However, a comprehensive overview of current developments is not available. In this review, we have therefore assessed what has been done before and summarized the available knowledge in the field, which can be used to develop improved models for glaucoma research.