A Methodology for Individual-Specific Modeling of Rat Optic Nerve Head Biomechanics in Glaucoma

Biomechanical analysis of ocular diseases and its in vitro study methods

Ocular diseases are closely related to the physiological changes in the eye sphere and its contents. Using biomechanical methods to explore the relationship between the structure and function of ocular tissue is beneficial to reveal the pathological processes. Studying the pathogenesis of various ocular diseases will be helpful for the diagnosis and treatment of ocular diseases. We provide a critical review of recent biomechanical analysis of ocular diseases including glaucoma, high myopia, and diabetes. And try to summarize the research about the biomechanical changes in ocular tissues (e.g., optic nerve head, sclera, cornea, etc.) associated with those diseases. The methods of ocular biomechanics research in vitro in recent years are also reviewed, including the measurement of biomechanics by ophthalmic equipment, finite element modeling, and biomechanical analysis methods. And the preparation and application of microfluidic eye chips that emerged in recent years were summarized. It provides new inspiration and opportunity for the pathogenesis of eye diseases and personalized and precise treatment.

The role of elastin on the mechanical properties of the anterior leaflet in porcine tricuspid valves

Elastin is present in the extracellular matrix (ECM) of connective tissues, and its mechanical properties are well documented. In Marfan syndrome, however, the inability to properly code for the protein fibrillin-1 prematurely leads to the degradation and loss of elastin fiber integrity in the ECM. In this study, the role of elastin in the ECM of the anterior leaflet of the tricuspid valve was investigated by examining the biomechanical behavior of porcine leaflets before and after the application of the enzyme elastase. Five loading protocols were applied to the leaflet specimens in two groups (elastase-treated and control samples). The mechanical response following elastase application yielded a significantly stiffer material in both the radial and circumferential directions. At a physiological level of stress (85 kPa), the elastase group had an average strain of 26.21% and 6.32% in the radial and circumferential directions, respectively, at baseline prior to elastase application. Following elastase treatment, the average strain was 5.28% and 0.97% in the radial and circumferential directions, respectively. No statistically significant change was found in the control group following sham treatment with phosphate-buffered saline (PBS). Two-photon microscopy images confirmed that after the removal of elastin, the collagen fibers displayed a loss of undulation. With a significant reduction in radial compliance, the ability to withstand physiological loads may be compromised. As such, an extracellular matrix that is structurally deficient in elastin may hinder normal tricuspid valve function.

Vertical Position of the Central Retinal Vessel in the Optic Disc and Its Association With the Site of Visual Field Defects in Glaucoma

PURPOSE

The purpose of this study was to investigate the association between the vertical position of the central retinal vessel (CRV) within the optic nerve head (ONH) and the site of visual field defects (VFDs) in glaucoma.

DESIGN

Cross-sectional study.

METHODS

The vertical position of the CRV was identified in 134 glaucoma eyes and 61 normal eyes at the point at which CRV exited the lamina cribrosa (LC) onto the ONH surface, by using spectral-domain optical coherence tomography (exit position). The position was also identified at the entry point into the LC from the retrolaminar ONH region (entry position), which was little influenced by glaucomatous LC deformation, therefore close to the original position before the glaucoma development. Positions were compared among glaucoma eyes with different sites of VFDs, and between glaucoma and normal eyes.

RESULTS

In glaucoma eyes, the entry position of the CRV was in the superior ONH region in 63.0% of eyes with superior VFDs and in the inferior ONH region in 97.8% of eyes with inferior VFDs (P < .0001). The exit position exhibited a similar percentage. The vertical CRV positions were not significantly different between glaucoma and normal eyes, both at the entry and exit positions.

CONCLUSIONS

Eyes with CRVs in the superior ONH region were significantly more likely to form VFDs in the superior hemifields and vice versa. The vertical position of the CRV was little altered by the development of glaucoma. The original position of the CRV before the development of glaucoma may influence regional susceptibility to glaucomatous stress and may be useful in predicting initial sites of VFDs.

Individual-Specific Modeling of Rat Optic Nerve Head Biomechanics in Glaucoma

Glaucoma is the second leading cause of blindness worldwide and is characterized by the death of retinal ganglion cells (RGCs), the cells that send vision information to the brain. Their axons exit the eye at the optic nerve head (ONH), the main site of damage in glaucoma. The importance of biomechanics in glaucoma is indicated by the fact that elevated intraocular pressure (IOP) is a causative risk factor for the disease. However, exactly how biomechanical insult leads to RGC death is not understood. Although rat models are widely used to study glaucoma, their ONH biomechanics have not been characterized in depth. Therefore, we aimed to do so through finite element (FE) modeling. Utilizing our previously described method, we constructed and analyzed ONH models with individual-specific geometry in which the sclera was modeled as a matrix reinforced with collagen fibers. We developed eight sets of scleral material parameters based on results from our previous inverse FE study and used them to simulate the effects of elevated IOP in eight model variants of each of seven rat ONHs. Within the optic nerve, highest strains were seen inferiorly, a pattern that was consistent across model geometries and model variants. In addition, changing the collagen fiber direction to be circumferential within the peripapillary sclera resulted in more pronounced decreases in strain than changing scleral stiffness. The results from this study can be used to interpret data from rat glaucoma studies to learn more about how biomechanics affects RGC pathogenesis in glaucoma.

Numerical Investigation on the Role of Mechanical Factors Contributing to Globe Flattening in States of Elevated Intracranial Pressure

Flattening of the posterior eye globe in the magnetic resonance (MR) images is a sign associated with elevated intracranial pressure (ICP), often seen in people with idiopathic intracranial hypertension (IIH). The exact underlying mechanisms of globe flattening (GF) are not fully known but mechanical factors are believed to play a role. In the present study, we investigated the effects of material properties and pressure loads on GF. For this purpose, we used a generic finite element model to investigate the deformation of the posterior eyeball. The degree of GF in numerical models and the significance of different mechanical factors on GF were characterized using an automated angle-slope technique and a statistical measure. From the numerical models, we found that ICP had the most important role in GF. We also showed that the angle-slope graphs pertaining to MR images from five people with high ICP can be represented numerically by manipulating the parameters of the finite element model. This numerical study suggests that GF observed in IIH patients can be accounted for by the forces caused by elevation of ICP from its normal level, while material properties of ocular tissues, such as sclera (SC), peripapillary sclera (PSC), and optic nerve (ON), would impact its severity.

Pressure-Induced Changes in Astrocyte GFAP, Actin, and Nuclear Morphology in Mouse Optic Nerve

Purpose

To conduct quantitative analysis of astrocytic glial fibrillary acidic protein (GFAP), actin and nuclei distribution in mouse optic nerve (ON) and investigate changes in the measured features after 3 days of ocular hypertension (OHT).

Method

Serial cross-sections of 3-day microbead-induced OHT and control ONs were fluorescently labelled and imaged using confocal microscope. Eighteen structural features were measured from the acquired images, including GFAP coverage, actin area fraction, process thickness, and aspect ratio of cell nucleus. The measured features were analyzed for variations with axial locations along ON and radial zones transverse to ON, as well as for the correlations with degree of intraocular pressure (IOP) change.

Results

The most significant changes in structural features after 3-day OHT occurred in the unmyelinated ON region (R1), and the changes were greater with greater IOP elevation. Although the GFAP, actin, axonal, and ON areas all increased in 3-day OHT ONs in R1 (P ≤ 0.004 for all), the area fraction of GFAP actually decreased (P = 0.02), the actin area fraction was stable and individual axon compartments were unchanged in size. Within R1, the number of nuclear clusters increased (P < 0.001), but the mean size of nuclear clusters was smaller (P = 0.02) and the clusters became rounder (P < 0.001). In all cross-sections of control ONs, astrocytic processes were thickest in the rim zone compared with the central and peripheral zones (P ≤ 0.002 for both), whereas the overall process width in R1 decreased after 3 days of OHT (P < 0.001).

Conclusions

The changes in structure elucidated IOP-generated alterations that underlie astrocyte mechanotranslational responses relevant to glaucoma.

Decellularizing the Porcine Optic Nerve Head: Toward a Model to Study the Mechanobiology of Glaucoma

Purpose

Studying the extracellular matrix (ECM) remodeling of the lamina cribrosa in vivo can be extremely challenging and costly. There exist very few options for studying optic nerve head (ONH) mechanobiology in vitro that are able to reproduce the complex anatomic and biomechanical environment of the ONH. Herein, we have developed a decellularization procedure that will enable more anatomically relevant and cost-efficient future studies of ECM remodeling of the ONH.

Methods

Porcine posterior poles were decellularized using a detergent and enzyme-based decellularization protocol. DNA quantification and histology were used to investigate the effectiveness of the protocol. We subsequently investigated the ability of a polyethylene glycol (PEG)-based hydrogel to restore the ONH's ability to hold pressure following decellularization. Anterior-posterior displacement of the decellularized and PEG treated ONH in a pressure bioreactor was used to evaluate the biomechanical response of the ONH.

Results

DNA quantification and histology confirmed decellularization using Triton X-100 at low concentration for 48 hours successfully reduced the cellular content of the tissue by 94.9% compared with native tissue while preserving the ECM microstructure and basal lamina of the matrix. Infiltrating the decellularized tissues with PEG 6000 and PEG 10,000 hydrogel restored their ability to hold pressure, producing displacements similar to those measured for the non-decellularized control samples.

Conclusions

Our decellularized ONH model is capable of producing scaffolds that are cell-free and maintain the native ECM microstructure.

Translational Relevance

This model represents a platform to study the mechanobiology in the ONH and potentially for glaucoma drug testing.

Biomechanical properties of the rat sclera obtained with inverse finite element modeling

It is widely accepted that biomechanics plays an important role in glaucoma pathophysiology, but the mechanisms involved are largely unknown. Rats are a common animal model of glaucoma, and finite element models are being developed to provide much-needed insight into the biomechanical environment of the posterior rat eye. However, material properties of rat ocular tissues, including the sclera, are currently unknown. Since the sclera plays a major role in posterior ocular biomechanics, our goal was to use inverse finite element modeling to extract rat scleral material properties. We first used digital image correlation to measure scleral surface displacement during whole-globe inflation testing. We modeled the sclera as a nonlinear material with embedded collagen fibers and then fit modeled displacements to experimental data using a differential evolution algorithm. Subject-specific models were constructed in which 3 parameters described the stiffness of the ground substance and collagen fibers in the posterior eye, and 16 parameters defined the primary orientation and alignment of fibers within eight scleral sub-regions. We successfully extracted scleral material properties for eight rat eyes. Model displacements recreated general patterns of the experimental displacements but did not always match local patterns. The fiber directions and fiber concentration parameters were highly variable, but on average, fibers were aligned circumferentially and were more aligned in the peripapillary sclera than in the peripheral sclera. The material properties determined here will be used to inform future finite element models of the rat posterior eye with the goal of elucidating the role of biomechanics in glaucoma pathophysiology.

Factors affecting optic nerve head biomechanics in a rat model of glaucoma

Glaucoma is the leading cause of irreversible blindness and is characterized by the death of retinal ganglion cells, which carry vision information from the retina to the brain. Although it is well accepted that biomechanics is an important part of the glaucomatous disease process, the mechanisms by which biomechanical insult, usually due to elevated intraocular pressure (IOP), leads to retinal ganglion cell death are not understood. Rat models of glaucoma afford an opportunity for learning more about these mechanisms, but the biomechanics of the rat optic nerve head (ONH), a primary region of damage in glaucoma, are only just beginning to be characterized. In a previous study, we built finite-element models with individual-specific rat ONH geometries. Here, we developed a parametrized model of the rat ONH and used it to perform a sensitivity study to determine the influence that six geometric parameters and 13 tissue material properties have on rat optic nerve biomechanical strains due to IOP elevation. Strain magnitudes and patterns in the parametrized model generally matched those from individual-specific models, suggesting that the parametrized model sufficiently approximated rat ONH anatomy. Similar to previous studies in human eyes, we found that scleral properties were highly influential: the six parameters with highest influence on optic nerve strains were optic nerve stiffness, IOP, scleral thickness, the degree of alignment of scleral collagen fibres, scleral ground substance stiffness and the scleral collagen fibre uncrimping coefficient. We conclude that a parametrized modelling strategy is an efficient approach that allows insight into rat ONH biomechanics. Further, scleral properties are important influences on rat ONH biomechanics, and additional efforts should be made to better characterize rat scleral collagen fibre organization.

A neuroglia-based interpretation of glaucomatous neuroretinal rim thinning in the optic nerve head

Neuroretinal rim thinning (NRR) is a characteristic glaucomatous optic disc change. However, the precise mechanism of the rim thinning has not been completely elucidated. This review focuses on the structural role of the glioarchitecture in the formation of the glaucomatous NRR thinning. The NRR is a glia-framed structure, with honeycomb geometry and mechanically reinforced astrocyte processes along the transverse plane. When neural damage selectively involves the neuron and spares the glia, the gross structure of the tissue is preserved. The disorganization and loss of the glioarchitecture are the two hallmarks of optic nerve head (ONH) remodeling in glaucoma that leads to the thinning of NRR tissue upon axonal loss. This is in contrast to most non-glaucomatous optic neuropathies with optic disc pallor where hypertrophy of the glioarchitecture is associated with the seemingly absent optic disc cupping. Arteritic anterior ischemic optic neuropathy is an exception where pan-necrosis of ONH tissue leads to NRR thinning. Milder ischemia indicates selective neuronal loss that spares glia in non-arteritic anterior ischemic optic neuropathy. The biological reason is the heterogeneous glial response determined by the site, type, and severity of the injury. The neuroglial interpretation explains how the cellular changes underlie the clinical findings. Updated understandings on glial responses illustrate the mechanical, microenvironmental, and microglial modulation of activated astrocytes in glaucoma. Findings relevant to the possible mechanism of the astrocyte death in advanced glaucoma are also emerging. Ultimately, a better understanding of glaucomatous glial response may lead to glia-targeting neuroprotection in the future.

Compressive mechanical properties of rat and pig optic nerve head

Glaucoma is the leading cause of irreversible blindness worldwide. Elevated intraocular pressure (IOP), the primary risk factor for glaucoma, is thought to induce abnormally high strains in optic nerve head (ONH) tissues, which ultimately result in retinal ganglion cell damage and vision loss. The mechanisms by which excessive deformations result in vision loss remain incompletely understood. The ability of computational and in vitro models of the ONH to provide insight into these mechanisms, in many cases, depends on our ability to replicate the physiological environment, which in turn requires knowledge of tissue biomechanical properties. The majority of mechanical data published to date regarding the ONH has been obtained from tensile testing, yet compression has been shown to be the main mode of deformation in the ONH under elevated IOP. We have thus tested pig and rat ONH tissue using unconfined cyclic compression. The material constants C₁, obtained from fitting the stress vs. strain data with a neo-Hookean material model, were 428 [367, 488] Pa and 64 [53, 76] Pa (mean [95% Confidence Interval]) for pig and rat optic nerve head, respectively. Additionally, we investigated the effects of strain rate and tissue storage on C1 values. These data will inform future efforts to understand and replicate the in vivo biomechanical environment of the ONH.

Matching the LenStar optical biometer to A-Scan ultrasonography for use in small animal eyes with application to tree shrews

We describe an analysis strategy to obtain ultrasonography-matched axial dimensions of small animal eyes using the LenStar biometer. The LenStar optical low-coherence reflectometer is an attractive device for animal research due to its high precision, non-invasiveness, and the ability to measure the axial dimensions of cornea, anterior chamber, lens, vitreous chamber, and axial length. However, this optical biometer was designed for clinical applications in human eyes and its internal analysis provides inaccurate values when used on small eyes due to species-dependent differences in refractive indices and relative axial dimensions. The LenStar uses a near infrared light source to measure optical path lengths (OPLs) that are converted by the LenStar's EyeSuite software into geometrical lengths (GLs) based on the refractive indices and axial dimensions of the human eye. We present a strategy that extracts the OPLs, determines refractive indices specific for the small animal eye of interest and then calculates corrected GLs. The refractive indices are obtained by matching the LenStar values to ultrasonography values in the same eyes. As compared to ultrasounography, we found that the internal calculations of the LenStar underestimate the axial dimensions of all ocular compartments of the tree shrew eye: anterior segment depth by 6.17±4.50%, lens thickness by 1.37±3.06%, vitreous chamber depth by 29.23±2.35%, and axial length by 10.62±1.75%. Using tree shrew-specific refractive indices, the axial dimensions closely matched those measured by ultrasonography for each compartment. Our analysis strategy can be easily translated to other species by obtaining a similar paired data set using ultrasonography and LenStar, and applying our step by step procedures.