Correlation between local stress and strain and lamina cribrosa connective tissue volume fraction in normal monkey eyes

Lamina Cribrosa Insertions Into the Sclera Are Sparser, Narrower, and More Slanted in the Anterior Lamina

Purpose

The lamina cribrosa (LC) depends on the sclera for support. The support must be provided through the LC insertions. Although a continuous insertion over the whole LC periphery is often assumed, LC insertions are actually discrete locations where LC collagenous beams meet the sclera. We hypothesized that LC insertions vary in number, size, and shape by quadrant and depth.

Methods

Coronal cryosections through the full LCs from six healthy monkey eyes were imaged using instant polarized light microscopy. The images were registered into a stack, on which we manually marked LC insertion outlines, nothing their position in-depth and quadrant (inferior, superior, nasal, or temporal). From the marks, we determined the insertion number, width, angle to the canal wall (90 degrees = perpendicular), and insertion ratio (fraction of LC periphery represented by insertions). Using linear mixed effect models, we determined if the insertion characteristics were associated with depth or quadrant.

Results

Insertions in the anterior LC were sparser, narrower, and more slanted than those in deeper LC (P values < 0.001). There were more insertions spanning a larger ratio of the canal wall in the middle LC than in the anterior and posterior (P values < 0.001). In the nasal quadrant, the insertion angles were significantly smaller (P < 0.001).

Conclusions

LC insertions vary substantially and significantly over the canal. The sparser, narrower, and more slanted insertions of the anterior-most LC may not provide the robust support afforded by insertions of the middle and posterior LC. These variations may contribute to the progressive deepening of the LC and regional susceptibility to glaucoma.

Scheimpflug-Based Corneal Biomechanical Analysis As A Predictor of Glaucoma in Eyes With High Myopia

Purpose

To address if corneal biomechanical behavior has a predictive value for the presence of glaucomatous optical neuropathy in eyes with high myopia.

Patients and Methods

This observational cross-sectional study included 209 eyes from 108 consecutive patients, divided into four groups: high myopia and primary open-angle glaucoma (POAG) - HMG, n = 53; high myopia without POAG - HMNG, n = 53; non-myopic with POAG - POAG, n = 50; non-myopic and non-POAG- NMNG, n = 53. Biomechanical assessment was made through a Scheimpflug-camera-based technology. Receiver operating characteristic curves were made for the discrimination between groups. Multivariable logistic regression models were performed to address the predictive value of corneal biomechanics for the presence of glaucoma.

Results

Areas Under the Receiver Operating Characteristic (AUROCs) above 0.6 were found in 6 parameters applied to discriminate between HMG and HMNG and six parameters to discriminate between POAG and NMNG. The biomechanical models with the highest power of prediction for the presence of glaucoma included 5 parameters with an AUROC of 0.947 for eyes with high myopia and 6 parameters with an AUROC of 0.857 for non-myopic eyes. In the final model, including all eyes, and adjusted for the presence of high myopia, the highest power of prediction for the presence of glaucoma was achieved including eight biomechanical parameters, with an AUROC of 0.917.

Conclusion

Corneal biomechanics demonstrated differences in eyes with glaucoma and mainly in myopic eyes. A biomechanical model based on multivariable logistic regression analysis and adjusted for high myopia was built, with an overall probability of 91.7% for the correct prediction of glaucomatous damage.

Computational study of the mechanical behavior of the astrocyte network and axonal compartments in the mouse optic nerve head

Glaucoma is a blinding disease characterized by the degeneration of the retinal ganglion cell (RGC) axons at the optic nerve head (ONH). A major risk factor for glaucoma is the intraocular pressure (IOP). However, it is currently impossible to measure the IOP-induced mechanical response of the axons of the ONH. The objective of this study was to develop a computational modeling method to estimate the IOP-induced strains and stresses in the axonal compartments in the mouse astrocytic lamina (AL) of the ONH, and to investigate the effect of the structural features on the mechanical behavior. We developed experimentally informed finite element (FE) models of six mouse ALs to investigate the effect of structure on the strain responses of the astrocyte network and axonal compartments to pressure elevation. The specimen-specific geometries of the FE models were reconstructed from confocal fluorescent images of cryosections of the mouse AL acquired in a previous study that measured the structural features of the astrocytic processes and axonal compartments. The displacement fields obtained from digital volume correlation in prior inflation tests of the mouse AL were used to determine the displacement boundary conditions of the FE models. We then applied Gaussian process regression to analyze the effects of the structural features on the strain outcomes simulated for the axonal compartments. The axonal compartments experienced, on average, 6 times higher maximum principal strain but 1800 times lower maximum principal stress compared to those experienced by the astrocyte processes. The strains experienced by the axonal compartments were most sensitive to variations in the area of the axonal compartments. Larger axonal compartments that were more vertically aligned, closer to the AL center, and with lower local actin area fraction had higher strains. Understanding the factors affecting the deformation in the axonal compartments will provide insights into mechanisms of glaucomatous axonal damage.

A direct fiber approach to model sclera collagen architecture and biomechanics

Sclera collagen fiber microstructure and mechanical behavior are central to eye physiology and pathology. They are also complex, and are therefore often studied using modeling. Most models of sclera, however, have been built within a conventional continuum framework. In this framework, collagen fibers are incorporated as statistical distributions of fiber characteristics such as the orientation of a family of fibers. The conventional continuum approach, while proven successful for describing the macroscale behavior of the sclera, does not account for the sclera fibers are long, interwoven and interact with one another. Hence, by not considering these potentially crucial characteristics, the conventional approach has only a limited ability to capture and describe sclera structure and mechanics at smaller, fiber-level, scales. Recent advances in the tools for characterizing sclera microarchitecture and mechanics bring to the forefront the need to develop more advanced modeling techniques that can incorporate and take advantage of the newly available highly detailed information. Our goal was to create a new computational modeling approach that can represent the sclera fibrous microstructure more accurately than with the conventional continuum approach, while still capturing its macroscale behavior. In this manuscript we introduce the new modeling approach, that we call direct fiber modeling, in which the collagen architecture is built explicitly by long, continuous, interwoven fibers. The fibers are embedded in a continuum matrix representing the non-fibrous tissue components. We demonstrate the approach by doing direct fiber modeling of a rectangular patch of posterior sclera. The model integrated fiber orientations obtained by polarized light microscopy from coronal and sagittal cryosections of pig and sheep. The fibers were modeled using a Mooney-Rivlin model, and the matrix using a Neo-Hookean model. The fiber parameters were determined by inversely matching experimental equi-biaxial tensile data from the literature. After reconstruction, the direct fiber model orientations agreed well with the microscopy data both in the coronal plane (adjusted R2 = 0.8234) and in the sagittal plane (adjusted R2 = 0.8495) of the sclera. With the estimated fiber properties (C10 = 5746.9 MPa; C01 = -5002.6 MPa, matrix shear modulus 200 kPa), the model's stress-strain curves simultaneously fit the experimental data in radial and circumferential directions (adjusted R2's 0.9971 and 0.9508, respectively). The estimated fiber elastic modulus at 2.16% strain was 5.45 GPa, in reasonable agreement with the literature. During stretch, the model exhibited stresses and strains at sub-fiber level, with interactions among individual fibers which are not accounted for by the conventional continuum methods. Our results demonstrate that direct fiber models can simultaneously describe the macroscale mechanics and microarchitecture of the sclera, and therefore that the approach can provide unique insight into tissue behavior questions inaccessible with continuum approaches.

Histology of the optic nerve head with special reference to the layer-specific distribution of composite fibers at and near the lamina cribrosa: An immunohistochemical study using specimens from elderly donated cadavers

BACKGROUND

This study aimed to demonstrate the composite fibers of the lamina cribrosa (LC) and their layer-specific distributions. The elastic fiber-rich septa, showing a cribriform arrangement in the optic nerve, may continue into the LC.

METHODS

Orbital content, including the long course of the optic nerve, was obtained from 25 elderly cadavers. Sagittal and cross-sections were prepared from each specimen. In addition to elastica Masson staining, immunohistochemistry was performed for elastin, glial fibrillary acidic protein (GFAP), S100 protein (S100), and CD68 in microglia.

RESULTS

The LC beam usually had fewer elastic fibers than the septa, but an elastic fiber-rich zone was observed along the scleral flange. GFAP-positive fibers were rich in the prelaminar area, whereas S100-positive fibers were rich in all layers of the LC. Double-positive (GFAP+/S100+) fibers were present in the prelaminar area. In contrast, S100-single positive fibers were evident in the LC and retrolaminar areas and were likely to insert into a sclera-choroid border area. The density of macrophages and microglia was not different between the septa and LC. Individual variations were observed in the distribution and density of the nerve-associated fibrous tissues.

CONCLUSION

The LC beam was quite different from the septa in the composite fibers and architecture. Transverse fibers, dominant in the LC beam, corresponded to fibrous processes of astrocytes and other nerve-associated fibrous tissues. Many of these nerve elements suggest low mechanical properties of the LC.

Eye-specific 3D modeling of factors influencing oxygen concentration in the lamina cribrosa

Our goal was to identify the factors with the strongest influence on the minimum lamina cribrosa (LC) oxygen concentration as potentially indicative of conditions increasing hypoxia risk. Because direct measurement of LC hemodynamics and oxygenation is not yet possible, we developed 3D eye-specific LC vasculature models. The vasculature of a normal monkey eye was perfusion-labeled post-mortem. Serial cryosections through the optic nerve head were imaged using fluorescence and polarized light microscopy to visualize the vasculature and collagen, respectively. The vasculature within a 450 μm-thick region containing the LC - identified from the collagen, was segmented, skeletonized, and meshed for simulations. Using Monte Carlo sampling, 200 vascular network models were generated with varying vessel diameter, neural tissue oxygen consumption rate, inflow hematocrit, and blood pressures (arteriole, venule, anterior boundary, and posterior boundary). Factors were varied over ranges of baseline ±20% with uniform probability. For each model we first obtained the blood flow, and from this the neural tissue oxygen concentration. ANOVA was used to identify the factors with the strongest influence on the minimum (10th percentile) oxygen concentration in the LC. The three most influential factors were, in ranked order, vessel diameter, neural tissue oxygen consumption rate, and arteriole pressure. There was a strong interaction between vessel diameter and arteriole pressure whereby the impact of one factor was larger when the other factor was small. Our results show that, for the eye analyzed, conditions that reduce vessel diameter, such as vessel compression due to elevated intraocular pressure or gaze-induced tissue deformation, may particularly contribute to decreased LC oxygen concentration. More eyes must be analyzed before generalizing.

Microstructural Deformations Within the Depth of the Lamina Cribrosa in Response to Acute In Vivo Intraocular Pressure Modulation

Purpose

The lamina cribrosa (LC) is a leading target for initial glaucomatous damage. We investigated the in vivo microstructural deformation within the LC volume in response to acute IOP modulation while maintaining fixed intracranial pressure (ICP).

Methods

In vivo optic nerve head (ONH) spectral-domain optical coherence tomography (OCT) scans (Leica, Chicago, IL, USA) were obtained from eight eyes of healthy adult rhesus macaques (7 animals; ages = 7.9-14.4 years) in different IOP settings and fixed ICP (8-12 mm Hg). IOP and ICP were controlled by cannulation of the anterior chamber and the lateral ventricle of the brain, respectively, connected to a gravity-controlled reservoir. ONH images were acquired at baseline IOP, 30 mm Hg (H1-IOP), and 40 to 50 mm Hg (H2-IOP). Scans were registered in 3D, and LC microstructure measurements were obtained from shared regions and depths.

Results

Only half of the eyes exhibited LC beam-to-pore ratio (BPR) and microstructure deformations. The maximal BPR change location within the LC volume varied between eyes. BPR deformer eyes had a significantly higher baseline connective tissue volume fraction (CTVF) and lower pore aspect ratio (P = 0.03 and P = 0.04, respectively) compared to BPR non-deformer. In all eyes, the magnitude of BPR changes in the anterior surface was significantly different (either larger or smaller) from the maximal change within the LC (H1-IOP: P = 0.02 and H2-IOP: P = 0.004).

Conclusions

The LC deforms unevenly throughout its depth in response to IOP modulation at fixed ICP. Therefore, analysis of merely the anterior LC surface microstructure will not fully capture the microstructure deformations within the LC. BPR deformer eyes have higher CTVF than BPR non-deformer eyes.

Scleral Cross-Linking in Form-Deprivation Myopic Guinea Pig Eyes Leads to Glaucomatous Changes

Purpose

To investigate the potential glaucomatous changes caused by scleral cross-linking (CXL) in a guinea pig form-deprivation (FD) myopia model.

Methods

Eighty 4-week-old tricolor guinea pigs were divided into four groups: FD only, genipin CXL only, FD plus CXL, and control. Refractive error, axial length (AL), intraocular pressure (IOP), and structural and vasculature optic disc changes in optical coherence tomography (OCT) and OCT angiography (OCTA) were measured at baseline and day 21. CXL efficacy was evaluated by scleral rigidity Young's modulus values. Histological and molecular changes in the anterior chamber angle, retina, and sclera were assessed.

Results

Baseline parameters were similar among groups (P > 0.05). The FD plus CXL group at day 21 had the least increase of AL (0.14 ± 0.08 mm) and highest IOP elevation (31.5 ± 3.6 mmHg) compared with the FD-only group (AL: 0.68 ± 0.17 mm; IOP: 22.2 ± 2.6 mmHg) and the control group (AL: 0.24 ± 0.09 mm; IOP: 17.4 ± 1.8 mmHg) (all P < 0.001). OCT and OCTA parameters of the optic disc in the FD plus CXL group at day 21 showed glaucomatous changes and decreased blood flow signals. Sclera rigidity increased in the CXL and FD plus CXL groups. Advanced glycation end products deposited extensively in the retina, choroid, and sclera of FD plus CXL eyes.

Conclusions

CXL causes increased IOP and subsequent optic disc, anterior segment, and scleral changes while inhibiting myopic progression and axial elongation in FD guinea pig eyes. Therefore, applying CXL to control myopia raises safety concerns.

Dual Orexin Receptor Antagonist Attenuates Increases in IOP, ICP, and Translaminar Pressure Difference After Stimulation of the Hypothalamus in Rats

PURPOSE

Intraocular pressure (IOP) remains the only modifiable risk factor for glaucoma progression. Our previous discovery that stimulation of nuclei within the hypothalamus can modulate IOP, intracranial pressure (ICP), and translaminar pressure difference (TLPD) fluctuations led us to investigate this pathway further. Our purpose was to determine the role of orexin neurons, primarily located in the dorsomedial hypothalamus (DMH) and perifornical (PeF) regions of the hypothalamus, in modulating these pressures.

METHODS

Sprague Dawley rats were pretreated systemically with a dual orexin receptor antagonist (DORA-12) at 30 mg/Kg (n = 8), 10 mg/Kg (n = 8), or vehicle control (n = 8). The IOP, ICP, heart rate (HR), and mean arterial pressure (MAP) were recorded prior to and following excitation of the DMH/PeF using microinjection of the gamma-aminobutyric acid (GABA)A receptor antagonist bicuculline methiodide (BMI).

RESULTS

Administration of the DORA at 30 mg/Kg significantly attenuated peak IOP by 5.2 ± 3.6 mm Hg (P = 0.007). During the peak response period (8-40 minutes), the area under the curve (AUC) for the 30 mg/Kg DORA cohort was significantly lower than the control cohort during the same period (P = 0.04). IOP responses for peak AUC versus DORA dose, from 0 to 30 mg/Kg, were linear (R2 = 0.18, P = 0.04). The ICP responses during the peak response period (4-16 minutes) versus DORA dose were also linear (R2 = 0.24, P = 0.014). Pretreatment with DORA significantly decreased AUC for the TLPD following stimulation of the DMH/PeF (10 mg/kg, P = 0.045 and 30 mg/kg, P = 0.015).

CONCLUSIONS

DORAs have the potential to attenuate asynchronous changes in IOP and in ICP and to lessen the extent of TLPDs that may result from central nervous system (CNS) activation.

Molecular regulation of neuroinflammation in glaucoma: Current knowledge and the ongoing search for new treatment targets

Neuroinflammation relying on the inflammatory responses of glial cells has emerged as an impactful component of the multifactorial etiology of neurodegeneration in glaucoma. It has become increasingly evident that despite early adaptive and reparative features of glial responses, prolonged reactivity of the resident glia, along with the peripheral immune cells, create widespread toxicity to retinal ganglion cell (RGC) axons, somas, and synapses. As much as the synchronized responses of astrocytes and microglia to glaucoma-related stress or neuron injury, their bi-directional interactions are critical to build and amplify neuroinflammation and to dictate the neurodegenerative outcome. Although distinct molecular programs regulate somatic and axonal degeneration in glaucoma, inhibition of neurodegenerative inflammation can provide a broadly beneficial treatment strategy to rescue RGC integrity and function. Since inflammatory toxicity and mitochondrial dysfunction are converging etiological paths that can boost each other and feed into a vicious cycle, anti-inflammatory treatments may also offer a multi-target potential. This review presents an overview of the current knowledge on neuroinflammation in glaucoma with particular emphasis on the cell-intrinsic and cell-extrinsic factors involved in the reciprocal regulation of glial responses, the interdependence between inflammatory and mitochondrial routes of neurodegeneration, and the research aspects inspiring for prospective immunomodulatory treatments. With the advent of powerful technologies, ongoing research on molecular and functional characteristics of glial responses is expected to accumulate more comprehensive and complementary information and to rapidly move the field forward to safe and effective modulation of the glial pro-inflammatory activities, while restoring or augmenting the glial immune-regulatory and neurosupport functions.

Finite element modeling of the complex anisotropic mechanical behavior of the human sclera and pia mater

BACKGROUND AND OBJECTIVE

Accurate finite element (FE) simulation of the optic nerve head (ONH) depends on accurate mechanical properties of the load-bearing tissues. The peripapillary sclera in the ONH exhibits a depth-dependent, anisotropic, heterogeneous collagen fiber distribution. This study proposes a novel cable-in-solid modeling approach that mimics heterogeneous anisotropic collagen fiber distribution, validates the approach against published experimental biaxial tensile tests of scleral patches, and demonstrates its effectiveness in a complex model of the posterior human eye and ONH.

METHODS

A computational pipeline was developed that defines control points in the sclera and pia mater, distributes the depth-dependent circumferential, radial, and isotropic cable elements in the sclera and pia in a pattern that mimics collagen fiber orientation, and couples the cable elements and solid matrix using a mesh-free penalty-based cable-in-solid algorithm. A parameter study was performed on a model of a human scleral patch subjected to biaxial deformation, and computational results were matched to published experimental data. The new approach was incorporated into a previously published eye-specific model to test the method; results were then interpreted in relation to the collagen fibers' (cable elements) role in the resultant ONH deformations, stresses, and strains.

RESULTS

Results show that the cable-in-solid approach can mimic the full range of scleral mechanical behavior measured experimentally. Disregarding the collagen fibers/cable elements in the posterior eye model resulted in ∼20-60% greater tensile and shear stresses and strains, and ∼30% larger posterior deformations in the lamina cribrosa and peripapillary sclera.

CONCLUSIONS

The cable-in-solid approach can easily be implemented into commercial FE packages to simulate the heterogeneous and anisotropic mechanical properties of collagenous biological tissues.

Ocular biomechanics due to ground blast reinforcement

BACKGROUND AND OBJECTIVE

Bomb blast injuries exerts a shearing force on the air-tissue interfaces, causing devastating ocular injury from the blast wave. Improvised explosive devices (IEDs) are usually placed at different heights from the ground to induce more severe injury through ground blast reinforcement (GBR). However, there is still a lack of knowledge of the role of GBR and IED height from the ground on ocular biomechanics, and how they can affect the intraocular pressure (IOP) in the eye. This study aimed to estimate the IOP due to frontal IED explosion at different heights from the ground using a fluid-structure interaction model with and without GBR effects.

METHODS

A 2 kg IED was placed within 5 m of the victim at 5 different heights from the ground, including 0, 0.42, 0.85, 1.27, and 1.70 m. Two different blast formulations were used to simulate the IED explosion: (a) spherical airburst, with no amplification of the initial shock wave due to interaction with the ground-surface, and (b) hemispherical surface-burst, where the initial blast wave is immediately reflected and reinforced by the ground (GBR).

RESULTS

Results revealed that the blast wave due to GBR reaches to the skull prior to the IED blast itself. The GBR also reached to the skull ∼ 0.6 ms earlier when the IED was on the ground compared to the height of 1.70 m. The highest and lowest IOPs of ∼ 17,000 and ∼ 15,000 mmHg were observed at the IED heights of 1.70 and 0 m from the ground considering GBR. However, when the role of the GBR is ignored, IOP of ∼ 9,000 mmHg was observed regardless of the IED height from the ground. The deformation in the apex of the cornea was higher when considering the GBR (∼ 0.75 cm) versus no GBR (∼ 0.65 cm). Considering GBR led to higher stresses and strains in the sclera.

CONCLUSIONS

When the role of GBR was ignored, the results showed similar patterns and magnitudes of stresses and deformations in the skull and eye regardless of the height of the IED from the ground, which was not the case when GBR was considered. The findings of this study suggest the critical role of GBR in ocular blast simulations.

Modeling the biomechanics of the lamina cribrosa microstructure in the human eye

Glaucoma is among the leading causes of blindness worldwide that is characterized by irreversible damage to the retinal ganglion cell axons in the lamina cribrosa (LC) region of the optic nerve head (ONH), most often associated with elevated intraocular pressure (IOP). The LC is a porous, connective tissue structure that provides mechanical support to the axons as they exit the eye and the biomechanics of the LC microstructure likely play a crucial role in protecting the axons passing through it. There is a limited knowledge of the IOP-driven biomechanics of the LC microstructure, primarily due to its small size and the difficulty with imaging the LC both in vitro and in vivo. We present finite element (FE) models of three human eye posterior poles that include the LC microstructure and interspersed neural tissues (NT) composed of retinal axons that are constructed directly from segmented, binary images of the LC. These models were used to estimate the stresses and strains in the LC and NT for an acute IOP elevation from 0 to 45 mmHg and compared with identical models except that the LC was represented as a homogenized continuum material with either homogeneous isotropic neo-Hookean properties or heterogeneous properties derived from local connective tissue volume fraction (CTVF) and predominant LC beam orientation. Stresses and strains in the LC and NT microstructure were investigated, and results were compared against those from the models wherein the LC was represented as a homogenized continuum. The regionalized volumetric average stresses and strains showed that the microstructural model yielded similar patterns to our prior approach using an LC continuum representation with mapped LC CTVF/anisotropy, but the microstructural modeling approach allows analysis of the stresses and strains in the LC and NT separately. As expected, the LC beams carried most of the IOP load in the microstructural models but exhibited less strain, while the encapsulated NT exhibited lower stresses and much higher strains. Results also revealed that the continuum models underestimate the maximum strains in the LC beams and NT by a factor of 2-3. Microstructural modeling should provide greater insight into the biomechanical factors driving damage to the axons (NT) and LC connective tissue remodeling that occur in glaucoma. The methods presented are ideal for modeling any structure with a complex microstructure composed of different materials, such as trabecular bone, lung, and tissue engineering scaffolds such as decellularized LC. Matlab code for mesh generation from a segmented image stack of the microstructure is included as Supplemental Material. STATEMENT OF SIGNIFICANCE: Glaucoma is among the leading causes of blindness worldwide that is characterized by axon damage in the lamina cribrosa (LC) region of the eye. We present a new approach for finite element modeling the entire eye-specific 3D LC microstructure and the interspersed neural tissues, incorporated into an eye-specific posterior eye model that provides appropriate boundary and loading conditions. Results are presented for three human donor eyes, showing that prior modeling approaches underestimate the stresses and strains in the laminar microstructure. We constructed models from image stacks of the segmented microstructure (Matlab code included) using an approach that is ideal for modeling any structure with a complex microstructure composed of different materials, such as trabecular bone, lung, and tissue engineering scaffolds.

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.

So-Called Lamina Cribrosa Defects May Mitigate IOP-Induced Neural Tissue Insult

Purpose

The prevailing theory about the function of lamina cribrosa (LC) connective tissues is that they provide structural support to adjacent neural tissues. Missing connective tissues would compromise this support and therefore are regarded as “LC defects”, despite scarce actual evidence of their role. We examined how so-called LC defects alter IOP-related mechanical insult to the LC neural tissues.

Methods

We built numerical models incorporating LC microstructure from polarized light microscopy images. To simulate LC defects of varying sizes, individual beams were progressively removed. We then compared intraocular pressure (IOP)-induced neural tissue deformations between models with and without defects. To better understand the consequences of defect development, we also compared neural tissue deformations between models with partial and complete loss of a beam.

Results

The maximum stretch of neural tissues decreased non-monotonically with defect size. Maximum stretch in the model with the largest defect decreased by 40% in comparison to the model with no defects. Partial loss of a beam increased the maximum stretch of neural tissues in its adjacent pores by 162%, compared with 63% in the model with complete loss of a beam.

Conclusions

Missing LC connective tissues can mitigate IOP-induced neural tissue insult, suggesting that the role of the LC connective tissues is more complex than simply fortifying against IOP. The numerical models further predict that partial loss of a beam is biomechanically considerably worse than complete loss of a beam, perhaps explaining why defects have been reported clinically but partial beams have not.

Corneal Biomechanical Assessment with Ultra-High-Speed Scheimpflug Imaging During Non-Contact Tonometry: A Prospective Review

BACKGROUND

In recent years, increasing interest has arisen in the application of data from corneal biomechanics in many areas of ophthalmology, particularly to assist in the detection of early corneal ectasia or ectasia susceptibility, to predict corneal response to surgical or therapeutic interventions and in glaucoma management. Technology has evolved and, recently, the Scheimpflug principle was associated with a non-contact air-puff tonometer, allowing a thorough analysis of corneal biomechanics and a biomechanically corrected intraocular pressure assessment, opening up new perspectives both in ophthalmology and in other medical areas. Data from corneal biomechanics assessment are being integrated in artificial intelligence models in order to increase its value in clinical practice.

OBJECTIVE

To review the state of the art in the field of corneal biomechanics assessment with special emphasis to the technology based on ultra-high-speed Scheimpflug imaging during non-contact tonometry.

SUMMARY

A meticulous literature review was performed until the present day. We used 136 published manuscripts as our references. Both information from healthy individuals and descriptions of possible associations with systemic diseases are described. Additionally, it exposed information regarding several fields of ocular pathology, from cornea and ocular surface through areas of refractive surgery and glaucoma until vascular and structural diseases of the chorioretinal unit.

Biomechanical properties of retina and choroid: a comprehensive review of techniques and translational relevance

Studying the biomechanical properties of biological tissue is crucial to improve our understanding of disease pathogenesis. The biomechanical characteristics of the cornea, sclera and the optic nerve head have been well addressed with an extensive literature and an in-depth understanding of their significance whilst, in comparison, knowledge of the retina and choroid is relatively limited. Knowledge of these tissues is important not only to clarify the underlying pathogenesis of a wide variety of retinal and vitreoretinal diseases, including age-related macular degeneration, hereditary retinal dystrophies and vitreoretinal interface diseases but also to optimise the surgical handling of retinal tissues and, potentially, the design and properties of implantable retinal prostheses and subretinal therapies. Our aim with this article is to comprehensively review existing knowledge of the biomechanical properties of retina, internal limiting membrane (ILM) and the Bruch’s membrane–choroidal complex (BMCC), highlighting the potential implications for clinical and surgical practice. Prior to this we review the testing methodologies that have been used both in vitro, and those starting to be used in vivo to aid understanding of their results and significance.

Analysis of the effects of finite element type within a 3D biomechanical model of a human optic nerve head and posterior pole

BACKGROUND AND OBJECTIVE

Biomechanical stresses and strains can be simulated in the optic nerve head (ONH) using the finite element (FE) method, and various element types have been used. This study aims to investigate the effects of element type on the resulting ONH stresses and strains.

METHODS

A single eye-specific model was constructed using 3D delineations of anatomic surfaces in a high-resolution, fluorescent, 3D reconstruction of a human posterior eye, then meshed using our simple meshing algorithm at various densities using 4- and 10-noded tetrahedral elements, as well as 8- and 20-noded hexahedral elements. A mesh-free approach was used to assign heterogeneous, anisotropic, hyperelastic material properties to the lamina cribrosa, sclera and pia. The models were subjected to elevated IOP of 45 mmHg after pre-stressing from 0 to 10 mmHg, and solved in the open-source FE package Calculix; results were then interpreted in relation to computational time and simulation accuracy, using the quadratic hexahedral model as the reference standard.

RESULTS

The 10-noded tetrahedral and 20R-noded hexahedral elements exhibited similar scleral canal and laminar deformations, as well as laminar and scleral stress and strain distributions; the quadratic tetrahedral models ran significantly faster than the quadratic hexahedral models. The linear tetrahedral and hexahedral elements were stiffer compared to the quadratic element types, yielding much lower stresses and strains in the lamina cribrosa.

CONCLUSIONS

Prior studies have shown that 20-noded hexahedral elements yield the most accurate results in complex models. Results show that 10-noded tetrahedral elements yield very similar results to 20-noded hexahedral elements and so they can be used interchangeably, with significantly lower computational time. Linear element types did not yield acceptable results.

THE INFERENCE OF THE CHANGES OF AXONAL TRANSPORT OF OPTIC NERVE BY DEFORMATIONS OF LAMINA CRIBROSA

Understanding the relationship between the changes in the axonal transport of the optic nerve (ON) and lamina cribrosa (LC) deformation will be helpful to estimate the degree of axonal transport block by measuring the LC deformation in vivo. First, the changes in the axonal transport of the ON were studied using an acute high intraocular pressure (IOP) model, which was established by perfusing saline water into the anterior chamber of cats. The IOP of cat was unilaterally elevated to and maintained at 30, 40, and 50[Formula: see text]mmHg. The axonal transport of the ON was examined by confocal laser scanning microscope. Then the deformations and stress distributions of the LC and ON were calculated using a three-dimensional finite element model of the LC microstructure including ON. The results showed axonal transport changes of ON increased with elevation of the IOPs. While Young’s modulus of the LC and ON were assumed as 0.1[Formula: see text]MPa and 0.03[Formula: see text]MPa, the numerical simulation results showed that LC had displacements of 0.02, 0.03, and 0.04[Formula: see text]mm backward at the IOPs of 30, 40, and 50[Formula: see text]mmHg, respectively. The calculated compressive strain applied to the ON were 0.0425, 0.0567, and 0.0709 under 30, 40, and 50[Formula: see text]mmHg IOP, respectively. The results of strain and stress analysis of LC and ON showed that the deformation of LC would compress the ON. The axonal transport abnormalities could be inferred by measuring the LC deformation in vivo.

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.

Morphometric, Hemodynamic, and Biomechanical Factors Influencing Blood Flow and Oxygen Concentration in the Human Lamina Cribrosa

Purpose

We developed a combined biomechanical and hemodynamic model of the human eye to estimate blood flow and oxygen concentration within the lamina cribrosa (LC) and rank the factors that influence LC oxygen concentration.

Methods

We generated 5000 finite-element eye models with detailed microcapillary networks of the LC and computed the oxygen concentration of the lamina retinal ganglion cell axons. For each model, we varied the intraocular pressure (IOP) from 10 mm Hg to 55 mm Hg in 5-mm Hg increments, the cerebrospinal fluid pressure (13 ± 2 mm Hg), cup depth (0.2 ± 0.1 mm), scleral stiffness (±20% of the mean values), LC stiffness (0.41 ± 0.2 MPa), LC radius (1.2 ± 0.12 mm), average LC pore size (5400 ± 2400 µm²), the microcapillary arrangement (radial, isotropic, or circumferential), and perfusion pressure (50 ± 9 mm Hg). Blood flow was assumed to originate from the LC periphery and drain via the central retinal vein. Finally, we performed linear regressions to rank the influence of each factor on the LC tissue oxygen concentration.

Results

LC radius and perfusion pressure were the most important factors in influencing the oxygen concentration of the LC. IOP was another important parameter, and eyes with higher IOP had higher compressive strain and slightly lower oxygen concentration. In general, superior–inferior regions of the LC had significantly lower oxygen concentration than the nasal–temporal regions, resulting in an hourglass pattern of oxygen deficiency.

Conclusions

To the best of our knowledge, this study is the first to implement a comprehensive hemodynamical model of the eye that accounts for the biomechanical forces and morphological parameters of the LC. The results provide further insight into the possible relationship of biomechanical and vascular pathways leading to ischemia-induced optic neuropathy.

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.

A broad perspective on the molecular regulation of retinal ganglion cell degeneration in glaucoma

Glaucoma is a complex neurodegenerative disease involving RGC axons, somas, and synapses at dendrites and axon terminals. Recent research advancements in the field have revealed a bigger picture of glaucomatous neurodegeneration that encompasses multiple stressors, multiple injury sites, multiple cell types, and multiple signaling pathways for asynchronous degeneration of RGCs during a chronic disease period. Optic nerve head is commonly viewed as the critical site of injury in glaucoma, where early injurious insults initiate distal and proximal signaling for axonal and somatic degeneration. Despite compartmentalized processes for degeneration of RGC axons and somas, there are intricate interactions between the two compartments and mechanistic overlaps between the molecular pathways that mediate degeneration in axonal and somatic compartments. This review summarizes the recent progress in the molecular understanding of RGC degeneration in glaucoma and highlights various etiological paths with biomechanical, metabolic, oxidative, and inflammatory components. Through this growing body of knowledge, the glaucoma community moves closer toward causative treatment of this blinding disease.

A Mesh-Free Approach to Incorporate Complex Anisotropic and Heterogeneous Material Properties into Eye-Specific Finite Element Models

Commercial finite element modeling packages do not have the tools necessary to effectively incorporate the complex anisotropic and heterogeneous material properties typical of the biological tissues of the eye. We propose a mesh-free approach to incorporate realistic material properties into finite element models of individual human eyes. The method is based on the idea that material parameters can be estimated or measured at so called control points, which are arbitrary and independent of the finite element mesh. The mesh-free approach approximates the heterogeneous material parameters at the Gauss points of each finite element while the boundary value problem is solved using the standard finite element method. The proposed method was applied to an eye-specific model a human posterior pole and optic nerve head. We demonstrate that the method can be used to effectively incorporate experimental measurements of the lamina cribrosa micro-structure into the eye-specific model. It was convenient to define characteristic material orientations at the anterior and posterior scleral surface based on the eye-specific geometry of each sclera. The mesh-free approach was effective in approximating these characteristic material directions with smooth transitions across the sclera. For the first time, the method enabled the incorporation of the complex collagen architecture of the peripapillary sclera into an eye-specific model including the recently discovered meridional fibers at the anterior surface and the depth dependent width of circumferential fibers around the scleral canal. The model results suggest that disregarding the meridional fiber region may lead to an underestimation of local strain concentrations in the retina. The proposed approach should simplify future studies that aim to investigate collagen remodeling in the sclera and optic nerve head or in other biological tissues with similar challenges.

Association between peripapillary scleral deformation and choroidal microvascular circulation in glaucoma

Peripapillary vessel density, which is reduced in eyes with glaucoma, has been proposed as a diagnostic tool for the desease and peripapillary choroidal microvasculature dropout(MvD) is considered one of pathophysiological manifestation of glaucomatous damage. However, little is known about the underlying pathogenic mechanism of dropout. According to recent studies, MvD is associated with structural changes in ONH structures. Therefore, we investigated the association between peripapillary scleral deformation and MvD. Data from 62 open-angle glaucoma (OAG) eyes with MvD and 36 eyes without MvD were analyzed in this study. And eyes with MvD were classified into two groups based on location: a juxtapapillary group and a non-juxtapapillary group for further analysis. More eyes with MvD had focal scleral deformation than did those without MvD (64.5% versus 2.8%; P < 0.001). Peripapillary choroidal thickness and focal scleral deformation were significantly associated with MvD. And juxtapapillary group was more associated with focal scleral deformation and coincidental RNFL defects than non-juxtapapillary groups. Peripapillary choroidal MvD was associated with the presence of scleral deformation, especially with juxtapapillary MvD, which was related to corresponding RNFL defects.

A Protective Eye Shield Reduces Limbal Strain and Its Variability During Simulated Sleep in Adults With Glaucoma

PURPOSE

To determine the effect of wearing a protective eye shield (mask) on limbal strain magnitude and variability in glaucoma eyes when sleeping with 1 side of the face down (FD) against a pillow.

METHODS

A prospective, randomized, interventional trial was conducted at the Wilmer Eye Institute with 36 glaucoma patients. A contact lens sensor measured limbal strain (output in equivalent millivolts) during intervals of up to 60 minutes in lateral decubitus, FD, and supine positions. Eighteen subjects wore a mask during 1 of 2 FD intervals, with randomized assignment of the interval. Data from additional trials with no mask were included in some analyses. In addition, some facial-feature dimensions from 3D scanned images of 23 subjects were compared with limbal strain data.

RESULTS

Wearing a mask trends toward a reduced mean change in contact lens sensor output (limbal strain) on moving to a FD positions [+34.1 mVeq, P=0.01 reduced by -22.3 mVeq, P=0.09 (n=36)]. Mask wearing reduced variability in strain while FD [-22.8 mVeq, P=0.04 (n=18)]. In eyes with past progressive visual field loss, the effect of the mask reduced mean strain change when moving to FD [-44.8 mVeq, P=0.02 (n=31)]. Longer corneal apex to nose-tip and to temple lengths were associated with reduced variability while FD [P=0.02 and 0.04, respectively (n=23)]. Treating both lengths as confounding factors increased statistical significance, particularly for analysis of the no-mask change in strain data moving to and from the FD position [P=0.004 to 0.002 and P=0.03 to 0.01 (n=23)].

CONCLUSION AND RELEVANCE

Wearing a mask reduced limbal strain and variation in limbal strain during simulated FD sleep, particularly in eyes with past field worsening, as did some facial features.

Primary Acute Angle-Closure Glaucoma: Three-Dimensional Reconstruction Imaging of Optic Nerve Heard Structure in Based on Optical Coherence Tomography (OCT)

Background

In glaucoma, the cup to plate ratio enlargement is a recognized pathological phenomenon. At present, the research on optic nerve in China and abroad mainly focuses on 2-dimensional research, and the measurement of 3-dimensional volume data is less well studied. Therefore, the recognition of 3-dimensional morphological changes is conducive to timely clinical intervention to prevent or slow down progressive vision loss.

Material/Methods

In this paper, optical coherence tomography (OCT) volume imaging technology was used to analyze and compare the morphological changes of primary acute angle-closure glaucoma in three-dimensional morphology, reconstruct the volume data of three-dimensional optic nerve head (ONH), and make morphological measurements.

Results

The rim width of the glaucoma group was significantly lower than that of the control group, and the average volume and intraocular pressure of the optic cup were significantly increased (P<0.05), while the rim width and intraocular pressure of the other group were not significantly changed (P>0.05).

Conclusions

We used three-dimensional reconstruction to identify OCT images between glaucoma patients and the control group with significant differences.

Quantification of collagen fiber structure using second harmonic generation imaging and two-dimensional discrete Fourier transform analysis: Application to the human optic nerve head

Second harmonic generation (SHG) microscopy is widely used to image collagen fiber microarchitecture due to its high spatial resolution, optical sectioning capabilities and relatively nondestructive sample preparation. Quantification of SHG images requires sensitive methods to capture fiber alignment. This article presents a two-dimensional discrete Fourier transform (DFT)-based method for collagen fiber structure analysis from SHG images. The method includes integrated periodicity plus smooth image decomposition for correction of DFT edge discontinuity artefact, avoiding the loss of peripheral image data encountered with more commonly used windowing methods. Outputted parameters are as follows: the collagen fiber orientation distribution, aligned collagen content and the degree of collagen fiber dispersion along the principal orientation. We demonstrate its application to determine collagen microstructure in the human optic nerve head, showing its capability to accurately capture characteristic structural features including radial fiber alignment in the innermost layers of the bounding sclera and a circumferential collagen ring in the mid-stromal tissue. Higher spatial resolution rendering of individual lamina cribrosa beams within the nerve head is also demonstrated. Validation of the method is provided in the form of correlative results from wide-angle X-ray scattering and application of the presented method to other fibrous tissues.

Age-Dependent Variation of Lamina Cribrosa Displacement During the Standardized Valsalva Maneuver

Based on biomechanical theory, lamina cribrosa (LC) displacement, the key component of progressive glaucomatous change, is presumed to be dependent on intraocular pressure (IOP) as well as tissue stiffness of LC. In the performance of the Valsalva maneuver, both IOP and cerebrospinal fluid pressure can increase. The present study investigated the age-dependent variation of LC displacement during the standardized Valsalva maneuver in healthy subjects. Sixty-three (63) eyes (age range: 20-76 years) were prospectively underwent IOP measurement and Cirrus HD-OCT optic disc scans before and during the standardized Valsalva maneuver. During the standardized Valsalva maneuver, the IOP significantly increased from 13.2 ± 2.9 mmHg to 18.6 ± 5.2 mmHg (P < 0.001). The maximal LC depth significantly decreased in the younger age groups (age: 20 s to 40 s) but not in the older age groups (age: over 50). The BMO distance did not change significantly. Younger age (P = 0.009), a smaller increase of IOP during the Valsalva maneuver (P = 0.002), and greater baseline maximal LC depth (P = 0.013) were associated with more anterior displacement of the LC during the standardized Valsalva maneuver. Taken together, age as well as translaminar pressure dynamics seems to play a crucial role in LC biomechanics.

Thin Lamina Cribrosa Beams Have Different Collagen Microstructure Than Thick Beams

Purpose

To compare the collagen microstructural crimp characteristics between thin and thick lamina cribrosa (LC) beams.

Methods

Seven eyes from four sheep were fixed at 5 mm Hg IOP in 10% formalin. For each eye, one to three coronal cryosections through the LC were imaged with polarized light microscopy and analyzed to visualize the LC and determine collagen fiber microstructure. For every beam, we measured its width and three characteristics of the crimp of its collagen fibers: waviness, tortuosity, and amplitude. Linear mixed effects models were used to test whether crimp characteristics were associated with the LC beam width.

Results

For each eye and over all the eyes, LC beam width was positively associated with crimp waviness and tortuosity, and negatively associated with crimp amplitude (P's < 0.0001). Thin beams, average width 13.11 μm, had average (SD) waviness, tortuosity, and amplitude of 0.27 (0.17) radians, 1.017 (0.028) and 1.88 (1.41) μm, respectively. For thick beams, average width 26.10 μm, these characteristics were 0.33 (0.18) radians, 1.025 (0.037) and 1.58 (1.36) μm, respectively.

Conclusions

Our results suggest heterogeneity in LC beam mechanical properties. Thin beams were less wavy than their thicker counterparts, suggesting that thin beams may stiffen at lower IOP than thick beams. This difference may allow thin beams to support similar amounts of IOP-induced force as thicker beams, thus providing a similar level of structural support to the axons at physiologic IOP, despite the differences in width. Measurements of beam-level mechanical properties are needed to confirm these predictions.

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

Purpose

To measure the deformation of the porcine optic nerve head (ONH) and peripapillary sclera (PPS) in response to intraocular pressure (IOP) elevation.

Methods

High-frequency ultrasound was used to image the ONH and PPS of 12 porcine eyes during ex vivo inflation testing from 5 to 30 mm Hg. A speckle tracking algorithm was used to compute tissue displacements in the anterior-posterior direction and expansion of the scleral canal. Through-thickness, in-plane, and shear strains were calculated within the ONH. Regional displacements and strains were analyzed and compared.

Results

The ONH and PPS showed overall posterior displacement in response to IOP elevation. Posterior displacement of the ONH was larger than and strongly correlated with the posterior displacement of the PPS throughout inflation testing. Scleral canal expansion was much smaller and leveled off quicker than ONH posterior displacement as IOP increased. Through-thickness compression was concentrated in the anterior ONH, which also experienced larger in-plane and shear strains than the posterior ONH. Within the anterior ONH, all three strains were significantly higher in the periphery compared with the center, with the shear strain exhibiting the greatest difference between the two regions.

Conclusions

High-resolution ultrasound speckle tracking revealed the full-thickness mechanical response of the posterior eye to IOP elevation. A mismatch in posterior displacement was found between the ONH and PPS, and regional analyses showed a concentration of strains within the periphery of the anterior porcine ONH. These deformation patterns may help in understanding IOP-associated optic nerve damage and glaucoma susceptibility.

Cerebrospinal Fluid Pressure: Revisiting Factors Influencing Optic Nerve Head Biomechanics

Purpose

To model the sensitivity of the optic nerve head (ONH) biomechanical environment to acute variations in IOP, cerebrospinal fluid pressure (CSFP), and central retinal artery blood pressure (BP).

Methods

We extended a previously published numerical model of the ONH to include 24 factors representing tissue anatomy and mechanical properties, all three pressures, and constraints on the optic nerve (CON). A total of 8340 models were studied to predict factor influences on 98 responses in a two-step process: a fractional factorial screening analysis to identify the 16 most influential factors, followed by a response surface methodology to predict factor effects in detail.

Results

The six most influential factors were, in order: IOP, CON, moduli of the sclera, lamina cribrosa (LC) and dura, and CSFP. IOP and CSFP affected different aspects of ONH biomechanics. The strongest influence of CSFP, more than twice that of IOP, was on the rotation of the peripapillary sclera. CSFP had similar influence on LC stretch and compression to moduli of sclera and LC. On some ONHs, CSFP caused large retrolamina deformations and subarachnoid expansion. CON had a strong influence on LC displacement. BP overall influence was 633 times smaller than that of IOP.

Conclusions

Models predict that IOP and CSFP are the top and sixth most influential factors on ONH biomechanics. Different IOP and CSFP effects suggest that translaminar pressure difference may not be a good parameter to predict biomechanics-related glaucomatous neuropathy. CON may drastically affect the responses relating to gross ONH geometry and should be determined experimentally.

Long-Term Shape, Curvature, and Depth Changes of the Lamina Cribrosa after Trabeculectomy

PURPOSE

To evaluate changes in lamina cribrosa (LC) shape, curvature, and depth after trabeculectomy.

DESIGN

Prospective, observational case series.

PARTICIPANTS

A total of 112 patients (118 eyes) with open- or closed-angle glaucoma undergoing trabeculectomy.

METHODS

The optic nerve head was imaged using enhanced depth imaging spectral-domain OCT before trabeculectomy and at 6 follow-up visits throughout the first postoperative year. The anterior LC surface and Bruch's membrane opening were marked in the serial horizontal B scans for the analysis of LC parameters using Morphology 1.0 software. Postoperative morphologic LC changes were assessed.

MAIN OUTCOME MEASURES

The postoperative LC global shape index (GSI), nasal-temporal (N-T) and superior-inferior (S-I) curvatures, and mean and sectoral LC depth (LCD).

RESULTS

The mean LC GSI increased only during the early postoperative period (P = 0.02), resulting in a change toward the saddle-rut shape. There was a flattening of the LC curvature in N-T (P < 0.001) and S-I (P = 0.003) meridians 12 months after trabeculectomy. A shallowing of the mean and sectoral LCD from baseline was significant throughout the entire follow-up period (P < 0.001) and progressed up to postoperative month 6. Twenty-eight patients showed a deepening of the LC from baseline in at least 1 visit. Eyes with shallower LCD compared with baseline responded to intraocular pressure (IOP) reduction with greater movement anteriorly than eyes with deeper LCD (P = 0.002). Greater IOP reduction (P = 0.007), less retinal nerve fiber layer thinning over the year (P = 0.003), and more superiorly-inferiorly curved baseline LC (P = 0.001) were associated with an increase in GSI. Younger age and IOP reduction were related to LC shallowing (P < 0.001, P = 0.002) and N-T flattening (P < 0.001).

CONCLUSIONS

In most eyes, trabeculectomy resulted in long-term flattening and shallowing of the LC. However, in some eyes, LC deepened from baseline. Change in LC global shape appeared to be temporal. Reduction in IOP plays an important role in the early phase of LC change; however, in the later phase, LC remodeling may play a crucial role in view of stable IOP.

In vivo assessment of changes in corneal hysteresis and lamina cribrosa position during acute intraocular pressure elevation in eyes with markedly asymmetrical glaucoma

Purpose

To investigate the biomechanical response of the cornea, lamina cribrosa (LC), and prelaminar tissue (PT) to an acute intraocular pressure (IOP) increase in patients with markedly asymmetrical glaucoma and in healthy controls.

Patients and methods

A total of 24 eyes of 12 patients with markedly asymmetrical primary open-angle glaucoma (POAG) and 12 eyes of 12 healthy patients were examined with spectral-domain optical coherence tomography (SD-OCT) and ocular response analyzer (ORA) at baseline and during acute IOP elevation by means of an ophthalmodynamometer. The displacement of the LC and PT and the change in corneal hysteresis (CH) and corneal resistance factor (CRF) were evaluated.

Results

Following a mean IOP increase of 12.3±2.4 mmHg, eyes with severe glaucoma demonstrated an overall mean anterior displacement of the LC (−6.58±26.09 µm) as opposed to the posterior laminar displacement in eyes with mild glaucoma (29.08±19.28 µm) and in healthy eyes (30.3±10.9; p≤0.001 and p=0.001, respectively). The PT displaced posteriorly during IOP elevation in all eyes. The CH decreased in eyes with severe glaucoma during IOP elevation (from 9.30±3.65 to 6.92±3.04 mmHg; p=0.012), whereas the CRF increased markedly in eyes with mild glaucoma (from 8.61±2.30 to 12.38±3.64; p=0.002) and in eyes with severe glaucoma (from 9.02±1.48 to 15.20±2.06; p=0.002). The increase in CRF correlated with the anterior displacement of the LC in eyes with severe glaucoma.

Conclusion

Eyes with severe glaucoma exhibited a mean overall anterior displacement of the anterior laminar surface, while eyes with mild glaucoma and healthy eyes showed a posterior displacement of the LC during IOP elevation. The CH decreased significantly from baseline only in eyes with severe glaucoma, but the CRF increased significantly in all glaucomatous eyes. The CRF increase correlated with the anterior displacement of the LC in eyes with severe glaucoma.

Imaging of the lamina cribrosa and its role in glaucoma: a review

The lamina cribrosa of the optic nerve head serves two contrasting roles; it must be porous to allow retinal ganglion cell axons to pass through, and yet at the same time, it must also provide adequate structural support to withstand the stresses and strains across it. Improvements in imaging such as optical coherence tomography image capture and image processing have allowed detailed in vivo studies of lamina cribrosa macro- and micro-architectural characteristics. This has aided our understanding of the optic nerve head as a complex biomechanical structure. In this review, we first aim to frame the biomechanical considerations of lamina cribrosa in a clinical context; in doing so, we also explore the concept of the translaminar pressure difference. Second, we aim to highlight the technological advances in imaging the lamina cribrosa and its accompanying clinical implications, and future directions in this quickly progressing field.

Spaceflight-Induced Intracranial Hypertension and Visual Impairment: Pathophysiology and Countermeasures

Visual impairment intracranial pressure (VIIP) syndrome is considered an unexplained major risk for future long-duration spaceflight. NASA recently redefined this syndrome as Spaceflight-Associated Neuro-ocular Syndrome (SANS). Evidence thus reviewed supports that chronic, mildly elevated intracranial pressure (ICP) in space (as opposed to more variable ICP with posture and activity on Earth) is largely accounted for by loss of hydrostatic pressures and altered hemodynamics in the intracranial circulation and the cerebrospinal fluid system. In space, an elevated pressure gradient across the lamina cribrosa, caused by a chronic but mildly elevated ICP, likely elicits adaptations of multiple structures and fluid systems in the eye which manifest themselves as the VIIP syndrome. A chronic mismatch between ICP and intraocular pressure (IOP) in space may acclimate the optic nerve head, lamina cribrosa, and optic nerve subarachnoid space to a condition that is maladaptive to Earth, all contributing to the pathogenesis of space VIIP syndrome. Relevant findings help to evaluate whether artificial gravity is an appropriate countermeasure to prevent this seemingly adverse effect of long-duration spaceflight.

Lamina Cribrosa Pore Shape and Size as Predictors of Neural Tissue Mechanical Insult

Purpose

The purpose of this study was to determine how the architecture of the lamina cribrosa (LC) microstructure, including the shape and size of the lamina pores, influences the IOP-induced deformation of the neural tissues within the LC pores using computational modeling.

Methods

We built seven specimen-specific finite element models of LC microstructure with distinct nonlinear anisotropic properties for LC beams and neural tissues based on histological sections from three sheep eyes. Changes in shape (aspect ratio and convexity) and size (area and perimeter length) due to IOP-induced hoop stress were calculated for 128 LC pores. Multivariate linear regression was used to determine if pore shape and size were correlated with the strain in the pores. We also compared the microstructure models to a homogenized model built following previous approaches.

Results

The LC microstructure resulted in focal tensile, compressive, and shear strains in the neural tissues of the LC that were not predicted by homogenized models. IOP-induced hoop stress caused pores to become larger and more convex; however, pore aspect ratio did not change consistently. Peak tensile strains within the pores were well predicted by a linear regression model considering the initial convexity (negative correlation, P < 0.001), aspect ratio (positive correlation, P < 0.01), and area (negative correlation, P < 0.01). Significant correlations were also found when considering the deformed shape and size of the LC pores.

Conclusions

The deformation of the LC neural tissues was largely dependent on the collagenous LC beams. Simple measures of LC pore shape and area provided good estimates of neural tissue biomechanical insult.

Formalin Fixation and Cryosectioning Cause Only Minimal Changes in Shape or Size of Ocular Tissues

Advances in imaging have made it increasingly common to study soft tissues without first embedding them in plastic or paraffin and without using labels or stains. The process, however, usually still involves fixation and cryosectioning, which could deform the tissues. Our goal was to quantify the morphological changes of ocular tissues caused by formalin fixation and cryosectioning. From each of 6 porcine eyes, 4 regions were obtained: cornea, equatorial and posterior sclera, and posterior pole containing the optic nerve head. Samples were imaged using visible light microscopy fresh, 1-minute and 24-hours post-fixation, and post-cryosectioning. Effects were assessed by 14 parameters representing sample size and shape. Overall, formalin fixation and sectioning caused only minimal changes to the ocular tissues, with average percentage parameter differences of 0.1%, 1%, and 1.2% between fresh and post-fixing by 1 minute, 24 hours, and post-cryosectioning, respectively. Parameter changes were not directional, and were only weakly dependent on the duration of fixation and the region of the eye. These results demonstrate that formalin fixation and cryosectioning are good choices for studying ocular tissue morphology and structure, as they do not cause the large tissue shrinkage or distortions typically associated with other, more complicated, techniques.

Effects of collagen microstructure and material properties on the deformation of the neural tissues of the lamina cribrosa

It is widely considered that intraocular pressure (IOP)-induced deformation within the neural tissue pores of the lamina cribrosa (LC) contributes to neurodegeneration and glaucoma. Our goal was to study how the LC microstructure and mechanical properties determine the mechanical insult to the neural tissues within the pores of the LC. Polarized light microscopy was used to measure the collagen density and orientation in histology sections of three sheep optic nerve heads (ONH) at both mesoscale (4.4μm) and microscale (0.73μm) resolutions. Mesoscale fiber-aware FE models were first used to calculate ONH deformations at an IOP of 30mmHg. The results were then used as boundary conditions for microscale models of LC regions. Models predicted large insult to the LC neural tissues, with 95th percentile 1st principal strains ranging from 7 to 12%. Pores near the scleral boundary suffered significantly higher stretch compared to pores in more central regions (10.0±1.4% vs. 7.2±0.4%; p=0.014; mean±SD). Variations in material properties altered the minimum, median, and maximum levels of neural tissue insult but largely did not alter the patterns of pore-to-pore variation, suggesting these patterns are determined by the underlying structure and geometry of the LC beams and pores. To the best of our knowledge, this is the first computational model that reproduces the highly heterogeneous neural tissue strain fields observed experimentally.

STATEMENT OF SIGNIFICANCE

The loss of visual function associated with glaucoma has been attributed to sustained mechanical insult to the neural tissues of the lamina cribrosa due to elevated intraocular pressure. Our study is the first computational model built from specimen-specific tissue microstructure to consider the mechanics of the neural tissues of the lamina separately from the connective tissue. We found that the deformation of the neural tissue was much larger than that predicted by any recent microstructure-aware models of the lamina. These results are consistent with recent experimental data and the highest deformations were found in the region of the lamina where glaucomatous damage first occurs. This study provides new insight into the complex biomechanical environment within the lamina.

The connective tissue phenotype of glaucomatous cupping in the monkey eye - Clinical and research implications

In a series of previous publications we have proposed a framework for conceptualizing the optic nerve head (ONH) as a biomechanical structure. That framework proposes important roles for intraocular pressure (IOP), IOP-related stress and strain, cerebrospinal fluid pressure (CSFp), systemic and ocular determinants of blood flow, inflammation, auto-immunity, genetics, and other non-IOP related risk factors in the physiology of ONH aging and the pathophysiology of glaucomatous damage to the ONH. The present report summarizes 20 years of technique development and study results pertinent to the characterization of ONH connective tissue deformation and remodeling in the unilateral monkey experimental glaucoma (EG) model. In it we propose that the defining pathophysiology of a glaucomatous optic neuropathy involves deformation, remodeling, and mechanical failure of the ONH connective tissues. We view this as an active process, driven by astrocyte, microglial, fibroblast and oligodendrocyte mechanobiology. These cells, and the connective tissue phenomena they propagate, have primary and secondary effects on retinal ganglion cell (RGC) axon, laminar beam and retrolaminar capillary homeostasis that may initially be "protective" but eventually lead to RGC axonal injury, repair and/or cell death. The primary goal of this report is to summarize our 3D histomorphometric and optical coherence tomography (OCT)-based evidence for the early onset and progression of ONH connective tissue deformation and remodeling in monkey EG. A second goal is to explain the importance of including ONH connective tissue processes in characterizing the phenotype of a glaucomatous optic neuropathy in all species. A third goal is to summarize our current efforts to move from ONH morphology to the cell biology of connective tissue remodeling and axonal insult early in the disease. A final goal is to facilitate the translation of our findings and ideas into neuroprotective interventions that target these ONH phenomena for therapeutic effect.

Roles of exosomes in the normal and diseased eye

Exosomes are nanometer-sized vesicles that are released by cells in a controlled fashion and mediate a plethora of extra- and intercellular activities. Some key functions of exosomes include cell-cell communication, immune modulation, extracellular matrix turnover, stem cell division/differentiation, neovascularization and cellular waste removal. While much is known about their role in cancer, exosome function in the many specialized tissues of the eye is just beginning to undergo rigorous study. Here we review current knowledge of exosome function in the visual system in the context of larger bodies of data from other fields, in both health and disease. Additionally, we discuss recent advances in the exosome field including use of exosomes as a therapeutic vehicle, exosomes as a source of biomarkers for disease, plus current standards for isolation and validation of exosome populations. Finally, we use this foundational information about exosomes in the eye as a platform to identify areas of opportunity for future research studies.

Deformation of the Lamina Cribrosa and Optic Nerve Due to Changes in Cerebrospinal Fluid Pressure

Purpose

Cerebrospinal fluid pressure (CSFp) changes are involved or implicated in various ocular conditions including glaucoma, idiopathic intracranial hypertension, and visual impairment and intracranial pressure syndrome. However, little is known about the effects of CSFp on lamina cribrosa and retrolaminar neural tissue (RLNT) biomechanics, potentially important in these conditions. Our goal was to use an experimental approach to visualize and quantify the deformation of these tissues as CSFp increased.

Methods

The posterior eye and RLNT of porcine eyes (n = 3) were imaged using synchrotron radiation phase-contrast micro-computed tomography (PC μCT) at an intraocular pressure of 15 mm Hg and CSFps of 4, 10, 20, and 30 mm Hg. Scans of each tissue region were acquired at each CSFp step and analyzed using digital volume correlation to determine 3-dimensional tissue deformations.

Results

Elevating CSFp increased the strain in the lamina cribrosa and RLNT of all three specimens, with the largest strains occurring in the RLNT. Relative to the baseline CSFp of 4 mm Hg, at 30 mm Hg, the lamina cribrosa experienced a mean first and third principal strain of 4.4% and −3.5%, respectively. The corresponding values for the RLNT were 9.5% and −9.1%.

Conclusions

CSFp has a significant impact on the strain distributions within the lamina cribrosa and, more prominently, within the RLNT. Elevations in CSFp were positively correlated with increasing deformations in each region and may play a role in ocular pathologies linked to changes in CSFp.

Lamina Cribrosa Microarchitecture in Monkey Early Experimental Glaucoma: Global Change

Purpose

The purpose of this study was to characterize experimental glaucoma (EG) versus control eye differences in lamina cribrosa (LC), beam diameter (BD), pore diameter (PD), connective tissue volume fraction (CTVF), connective tissue volume (CTV), and LC volume (LV) in monkey early EG.

Methods

Optic nerve heads (ONHs) of 14 unilateral EG and 6 bilateral normal (BN) monkeys underwent three-dimensional reconstruction and LC beam segmentation. Each beam and pore voxel was assigned a diameter based on the largest sphere that contained it before transformation to a common cylinder with inner, middle, and outer layers. Full-thickness and layer averages for BD, PD, CTVF, CTV, and LV were calculated for each ONH. Beam diameter and PD distributions for each ONH were fit to a gamma distribution and summarized by scale and shape parameters. Experimental glaucoma and depth effects were assessed for each parameter by linear mixed-effects (LME) modeling. Animal-specific EG versus control eye differences that exceeded the maximum intereye difference among the six BN animals were considered significant.

Results

Overall EG eye mean PD was 12.8% larger (28.2 ± 5.6 vs. 25.0 ± 3.3 μm), CTV was 26.5% larger (100.06 ± 47.98 vs. 79.12 ± 28.35 × 10⁶ μm³), and LV was 40% larger (229.29 ± 98.19 vs. 163.63 ± 39.87 × 10⁶ μm³) than control eyes (P ≤ 0.05, LME). Experimental glaucoma effects were significantly different by layer for PD (P = 0.0097) and CTVF (P < 0.0001). Pore diameter expanded consistently across all PDs. Experimental glaucoma eye-specific parameter change was variable in magnitude and direction.

Conclusions

Pore diameter, CTV, and LV increase in monkey early EG; however, EG eye-specific change is variable and includes both increases and decreases in BD and CTVF.

Displacement of the Lamina Cribrosa in Response to Acute Intraocular Pressure Elevation in Normal Individuals of African and European Descent

Purpose

To assess if the in vivo mechanical displacement of the anterior laminar cribrosa surface (ALCS) as a response of an acute elevation in intraocular pressure (IOP) differs in individuals of European (ED) and African descent (AD).

Methods

Spectral-domain optical coherence tomography (SDOCT) scans were obtained from 24 eyes of 12 individuals of AD and 18 eyes of 9 individuals of ED at their normal baseline IOP and after 60 seconds IOP elevation using ophthalmodynamometry. Change in depth (displacement) of the LC and to the prelaminar tissue (PLT) were computed in association with the change (delta) in IOP (Δ IOP), race, age, corneal thickness, corneal rigidity (ocular response analyzer [ORA]), and axial.

Results

In the ED group for small IOP elevations (Δ IOP < 12 mm Hg), the ALCS initially displaced posteriorly but for larger increase of IOP an anterior displacement of the lamina followed. Inversely, in the AD group the ALCS did not show a significant posterior displacement for small Δ IOP, while for larger IOP increases the ALCS significantly displaced posteriorly. Posterior displacement of the lamina cribrosa (LC) was also significantly correlated with longer axial length, higher corneal thickness, and ORA parameters. Prelaminar tissue posteriorly displaced for any magnitude of Δ IOP, in both groups.

Conclusions

The African descent group demonstrated a greater acute posterior bowing of the LC after adjustment for age, axial length, Bruch's membrane opening (BMO) area, and ORA parameters. Greater PLT posterior displacement was also seen in the AD group with increasing IOP, which was tightly correlated with the displacement of the LC.

SIMULATING THE EFFECTS OF ELEVATED INTRAOCULAR PRESSURE ON OCULAR STRUCTURES USING A GLOBAL FINITE ELEMENT MODEL OF THE HUMAN EYE

Elevated intraocular pressure (IOP) may be the primary risk factor to the development of glaucoma. Finite element (FE) modeling is commonly considered as an effective method to quantitatively analyze pathogenesis of glaucoma. Recent researches focus on establishing partial human eye models. A refined global human eye model was developed using ANSYS software to investigate the correlation between IOP elevation and biomechanical responses. First, the pressure transferring process according to IOP elevation in the whole eye was analyzed to simulate the effects of IOP elevation on glaucoma. Then, the biomechanical responses of the anterior eye segment under various pressure differences between the anterior and posterior chambers (AC and PC) were analyzed to simulate posterior nonadhesion of iris and posterior synechia. This global eye model not only simulated the responses of elevated IOP on ocular structures, but also revealed the process of pressure transferring among each tissue from the anterior eye segment to the optic nerve head (ONH) region. The local mechanical characteristics of the ocular structures obtained from the global model agreed with previous findings. This global model may shed light on the studies of multifactorial glaucoma.

Anisotropic Finite Element Modeling Based on a Harmonic Field for Patient-Specific Sclera

Purpose. This study examined the influence of anisotropic material for human sclera. Method. First, the individual geometry of patient-specific sclera was reproduced from a laser scan. Then, high quality finite element modeling of individual sclera was performed using a convenient automatic hexahedral mesh generator based on harmonic field and integrated with anisotropic material assignment function. Finally, comparison experiments were designed to investigate the effects of anisotropy on finite element modeling of sclera biomechanics. Results. The experimental results show that the presented approach can generate high quality anisotropic hexahedral mesh for patient-specific sclera. Conclusion. The anisotropy shows significant differences for stresses and strain distribution and careful consideration should be given to its use in biomechanical FE studies.