Material properties of the posterior human sclera

Revealing regional variations in scleral shear modulus in a rabbit eye model using multi-directional ultrasound optical coherence elastography

The mechanical properties of the sclera play a critical role in supporting the ocular structure and maintaining its shape. However, non-invasive measurements to quantify scleral biomechanics remain challenging. Recently introduced multi-directional optical coherence elastography (OCE) combined with an air-coupled ultrasound transducer for excitation of elastic surface waves was used to estimate phase speed and shear modulus in ex vivo rabbit globes (n = 7). The scleral phase speed (12.1 ± 3.2 m/s) was directional-dependent and higher than for corneal tissue (5.9 ± 1.4 m/s). In the tested locations, the sclera proved to be more anisotropic than the cornea by a factor of 11 in the maximum of modified planar anisotropy coefficient. The scleral shear moduli, estimated using a modified Rayleigh-Lamb wave model, showed significantly higher values in the circumferential direction (65.4 ± 31.9 kPa) than in meridional (22.5 ± 7.2 kPa); and in the anterior zone (27.3 ± 9.3 kPa) than in the posterior zone (17.8 ± 7.4 kPa). The multi-directional scanning approach allowed both quantification and radial mapping of estimated parameters within a single measurement. The results indicate that multi-directional OCE provides a valuable non-invasive assessment of scleral tissue properties that may be useful in the development of improved ocular models, the evaluation of potential myopia treatment strategies, and disease characterization and monitoring.

Biomechanical changes of tree shrew posterior sclera during experimental myopia, after retrobulbar vehicle injections, and crosslinking using genipin

Myopia is a common ocular condition characterized by biomechanical weakening revealed by increasing creep rate, cyclic softening scleral thinning, change of collagen fibril crimping, and excessive elongation of the posterior sclera resulting in blurred vision. Animal studies support scleral crosslinking as a potential treatment for myopia control by strengthening the weakened sclera and slowing scleral expansion. While multiple studies investigated aspects of the biomechanical weakening and strengthening effects in myopia and after scleral crosslinking, a comprehensive analysis of the underlying mechanical changes including the effect of vehicle injections is still missing. The purpose of this study was to provide a comprehensive analysis of biomechanical changes by scleral inflation testing in experimental myopia, after retrobulbar vehicle injections and scleral crosslinking using genipin in tree shrews. Our results suggest that biomechanical weakening in myopia involves an increased creep rate and higher strain levels at which collagen fibers uncrimp. Both weakening effects were reduced after scleral crosslinking using genipin at doses that were effective in slowing myopia progression. Vehicle injections increased mechanical hysteresis and had a small but significant effect on slowing myopia progression. Also, our results support scleral crosslinking as a potential treatment modality that can prevent or counteract scleral weakening effects in myopia. Furthermore, vehicle solutions may cause independent biomechanical effects, which should be considered when developing and evaluating scleral crosslinking procedures.

A simple computational model for scleral stiffness assessments via air-puff deformation OCT

Introduction: The mechanical properties of the sclera are related to its structural function, and changes to these properties are believed to contribute to pathologies such as myopia. Air-puff deformation imaging is a tool that uses an imaging system coupled with an air-puff excitation source to induce and measure deformation in a tissue in vivo. Typically used for the study of the cornea's mechanical properties and IOP, this tool has been proposed as a method to evaluate scleral stiffness. Methods: In this work, we present a computational model of the rabbit eye to assess scleral deformation under air-puff. Parametric studies were conducted to evaluate the effects of material properties, intraocular pressure, and other parameters on the deformation response. Output from the model was also compared to experimental measurements of air-puff deformation in rabbit eyes under varying IOP. Results: Central deformation response was found to be most influenced by material properties of the sclera (at site of air-puff and posterior), thickness, and IOP, whereas deformation profile was most influenced by material properties. Experimental and simulated IOP dependence were found to be similar (RMSE = 0.13 mm). Discussion: Scleral APDI could be a useful tool for quick in vivo assessment of scleral stiffness.

A Noninvasive Clinical Method to Measure in Vivo Mechanical Strains of the Lamina Cribrosa by OCT

Objective

To measure mechanical strain of the lamina cribrosa (LC) after intraocular pressure (IOP) change produced 1 week after a change in glaucoma medication.

Design

Cohort study.

Participants

Adult glaucoma patients (23 eyes, 15 patients) prescribed a change in IOP-lowering medication.

Intervention

Noninvasive OCT imaging of the eye.

Main Outcome Measures

Deformation calculated by digital volume correlation of OCT scans of the LC before and after IOP lowering by medication.

Results

Among 23 eyes, 17 eyes of 12 persons had IOP lowering ≥ 3 mmHg (reduced IOP group) with tensile anterior-posterior Ezz strain = 1.0% ± 1.1% (P = 0.003) and compressive radial strain (Err) = -0.3% ± 0.5% (P = 0.012; random effects models accounting inclusion of both eyes in some persons). Maximum in-plane principal (tensile) strain and maximum shear strain in the reduced-IOP group were as follows: Emax = 1.7% ± 1.0% and Γmax = 1.4% ± 0.7%, respectively (both P < 0.0001 vs. zero). Reduced-IOP group strains Emax and Γmax were significantly larger with greater % IOP decrease (P < 0.0001 and P < 0.0001, respectively). The compliances of the Ezz, Emax, and Γmax strain responses, defined as strain normalized by the IOP decrease, were larger with more abnormal perimetric mean deviation or visual field index values (all P ≤ 0.02). Strains were unrelated to age (all P ≥ 0.088). In reduced-IOP eyes, mean LC anterior border posterior movement was only 2.05 μm posteriorly (P = 0.052) and not related to % IOP change (P = 0.94, random effects models). Only Err was significantly related to anterior lamina depth change, becoming more negative with greater posterior LC border change (P = 0.015).

Conclusions

Lamina cribrosa mechanical strains can be effectively measured by changes in eye drop medication using OCT and are related to degree of visual function loss in glaucoma.

Financial Disclosures

Proprietary or commercial disclosure may be found in the Footnotes and Disclosures at the end of this article.

Effect of twisting of intravitreal injections on ocular bio-mechanics: a novel insight to ocular surgery

Although intravitreal (IVT) injections provide several advantages in treating posterior segment eye diseases, several associated challenges remain. The current study uses the finite element method (FEM) to highlight the effect of IVT needle rotation along the insertion axis on the reaction forces and deformation inside the eye. A comparison of the reaction forces at the eye's key locations has been made with and without rotation. In addition, a sensitivity analysis of various parameters, such as the needle's angular speed, insertion location, angle, gauge, shape, and intraocular pressure (IOP), has been carried out to delineate the individual parameter's effect on reaction forces during rotation. Results demonstrate that twisting the needle significantly reduces the reaction forces at the penetration location and throughout the needle travel length, resulting in quicker penetration. Moreover, ocular biomechanics are influenced by needle insertion location, angle, shape, size, and IOP. The reaction forces incurred by the patient may be reduced by using a bevel needle of the higher gauge when inserted close to the normal of the local scleral surface toward the orra serrata within the Pars Plana region. Results obtained from the current study can deepen the understanding of the twisting needle's interaction with the ocular tissue.

Postmortem Digital Image Correlation and Finite Element Modeling Demonstrate Posterior Scleral Deformations during Optic Nerve Adduction Tethering

Postmortem human eyes were subjected to optic nerve (ON) traction in adduction and elevated intraocular pressure (IOP) to investigate scleral surface deformations. We incrementally adducted 11 eyes (age 74.1 ± 9.3 years, standard deviation) from 26° to 32° under normal IOP, during imaging of the posterior globe, for analysis by three-dimensional digital image correlation (3D-DIC). In the same eyes, we performed uniaxial tensile testing in multiple regions of the sclera, ON, and ON sheath. Based on individual measurements, we analyzed eye-specific finite element models (FEMs) simulating adduction and IOP loading. Analysis of 3D-DIC showed that the nasal sclera up to 1 mm from the sheath border was significantly compressed during adduction. IOP elevation from 15 to 30 mmHg induced strains less than did adduction. Tensile testing demonstrated ON sheath stiffening above 3.4% strain, which was incorporated in FEMs of adduction tethering that was quantitatively consistent with changes in scleral deformation from 3D-DIC. Simulated IOP elevation to 30 mmHg did not induce scleral surface strains outside the ON sheath. ON tethering in incremental adduction from 26° to 32° compressed the nasal and stretched the temporal sclera adjacent to the ON sheath, more so than IOP elevation. The effect of ON tethering is influenced by strain stiffening of the ON sheath.

Two-dimensional elastic distribution imaging of the sclera using acoustic radiation force optical coherence elastography

The scleral elasticity is closely related with many ocular diseases, but the relevant research is still insufficient. Here, we utilized optical coherence elastography to carefully study biomechanical properties of the sclera at different positions and under different intraocular pressures. Meanwhile, elastic wave velocity and Young's modulus of each position were obtained using a phase velocity algorithm. Accordingly, the two-dimensional elasticity distribution image was achieved by mapping the Young's modulus values to the corresponding structure based on the relationship between the position and its Young's modulus. Therefore, elastic information in regions-of-interest can be read and compared directly from the scleral structure, indicating that our method may be a very useful tool to evaluate the elasticity of sclera and provide intuitive and reliable proof for diagnosis and research.

AI-based clinical assessment of optic nerve head robustness superseding biomechanical testing

BACKGROUND/AIMS

To use artificial intelligence (AI) to: (1) exploit biomechanical knowledge of the optic nerve head (ONH) from a relatively large population; (2) assess ONH robustness (ie, sensitivity of the ONH to changes in intraocular pressure (IOP)) from a single optical coherence tomography (OCT) volume scan of the ONH without the need for biomechanical testing and (3) identify what critical three-dimensional (3D) structural features dictate ONH robustness.

METHODS

316 subjects had their ONHs imaged with OCT before and after acute IOP elevation through ophthalmo-dynamometry. IOP-induced lamina cribrosa (LC) deformations were then mapped in 3D and used to classify ONHs. Those with an average effective LC strain superior to 4% were considered fragile, while those with a strain inferior to 4% robust. Learning from these data, we compared three AI algorithms to predict ONH robustness strictly from a baseline (undeformed) OCT volume: (1) a random forest classifier; (2) an autoencoder and (3) a dynamic graph convolutional neural network (DGCNN). The latter algorithm also allowed us to identify what critical 3D structural features make a given ONH robust.

RESULTS

All three methods were able to predict ONH robustness from a single OCT volume scan alone and without the need to perform biomechanical testing. The DGCNN (area under the curve (AUC): 0.76±0.08) outperformed the autoencoder (AUC: 0.72±0.09) and the random forest classifier (AUC: 0.69±0.05). Interestingly, to assess ONH robustness, the DGCNN mainly used information from the scleral canal and the LC insertion sites.

CONCLUSIONS

We propose an AI-driven approach that can assess the robustness of a given ONH solely from a single OCT volume scan of the ONH, and without the need to perform biomechanical testing. Longitudinal studies should establish whether ONH robustness could help us identify fast visual field loss progressors.

PRECIS

Using geometric deep learning, we can assess optic nerve head robustness (ie, sensitivity to a change in IOP) from a standard OCT scan that might help to identify fast visual field loss progressors.

The Effect of Scleral Buckle Surgery on Tonographic Outflow Facility, Positional Intraocular Pressure, and Ocular Biomechanics

PURPOSE

To investigate the in vivo effect of scleral buckle surgery on ocular biomechanics and aqueous humor dynamics.

DESIGN

Prospective observational cross-sectional study.

PARTICIPANTS

Nine patients with unilateral 360 degree encircling scleral buckles without vitrectomy for rhegmatogenous retinal detachments, between 3 and 39 months postoperative.

METHODS

All measurements were performed in both eyes of all participants. Intraocular pressure (IOP) was measured in the seated and supine positions using pneumatonometry. Outflow facility was measured using 2-minute weighted pneumatonography. Ocular rigidity coefficient was determined from the Friedenwald equations based on the difference in IOP with and without a weighted tonometer tip. The percentage change in IOP upon transitioning from seated to supine was calculated. Measurements for buckled and nonbuckled eyes were compared using paired Student t test of means.

MAIN OUTCOME MEASURES

Sitting and supine IOP and percentage difference between the 2 positions; outflow facility; ocular rigidity coefficient.

RESULTS

Seated IOP was similar between buckled and nonbuckled eyes (16.1 ± 2.5 vs. 16.7 ± 2.7 mmHg; P = 0.5) whereas supine IOP was lower in buckled eyes compared with nonbuckled eyes (18.7 ± 2.6 vs. 21.3 ± 2.5 mmHg; P = 0.008). The percentage increase in IOP upon change in body position from seated to supine was greater in nonbuckled eyes (17.4 ± 9.4% vs. 27.6 ± 9.5%; P = 0.005). Ocular rigidity coefficient was lower in buckled (9.9 × 10-3 ± 1.4 × 10-3 μL-1) vs. nonbuckled eyes (14.4 × 10-3 ± 3.1 × 10-3 μL-1; P = 0.006). Outflow facility was not significantly different in buckled and nonbuckled eyes.

CONCLUSIONS

Scleral buckling decreases ocular rigidity but does not affect outflow facility. This change in ocular biomechanics likely results in the attenuated IOP change from seated to supine position. Decreased ocular rigidity may also reduce IOP fluctuations and potentially reduce the risk for glaucoma progression.

FINANCIAL DISCLOSURE(S)

Proprietary or commercial disclosure may be found in the Footnotes and Disclosures at the end of this article.

Displacement of the Lamina Cribrosa With Acute Intraocular Pressure Increase in Brain-Dead Organ Donors

Purpose

To examine deformations of the optic nerve head (ONH) deep tissues in response to acute elevation of intraocular pressure (IOP).

Methods

Research-consented brain-dead organ donors underwent imaging by spectral domain optical coherence tomography (OCT). OCT imaging was repeated while the eye was sequentially maintained at manometric pressures of 10, 30, and 50 mm Hg. Radial scans of the ONH were automatically segmented by deep learning and quantified in three dimensions by a custom algorithm. Change in lamina cribrosa (LC) depth and choroidal thickness was correlated with IOP and age by linear mixed-effect models. LC depth was computed against commonly utilized reference planes.

Results

Twenty-six eyes from 20 brain-dead organ donors (age range, 22-62 years; median age, 43 years) were imaged and quantified. LC depth measured against a reference plane based on Bruch's membrane (BM), BM opening, and an anterior sclera canal opening plane showed both a reduction and an increase in LC depth with IOP elevation. LC depth universally increased in depth when measured against a sclera reference plane. Choroidal (-0.5222 µm/mm Hg, P < 0.001) and retinal nerve fiber layer thickness (-0.0717 µm/mm Hg, P < 0.001) significantly thinned with increasing IOP. The magnitude of LC depth change with IOP was significantly smaller with increasing age (P < 0.03 for all reference planes).

Conclusions

LC depth changes with IOP reduce with age and are significantly affected by the reference plane of choice, which highlights a need for standardizing LC metrics to properly follow progressive remodeling of the loadbearing tissues of the ONH by OCT imaging and for the definition of a reference database.

2D or not 2D? Mapping the in-depth inclination of the collagen fibers of the corneoscleral shell

The collagen fibers of the corneoscleral shell play a central role in the eye mechanical behavior. Although it is well-known that these fibers form a complex three-dimensional interwoven structure, biomechanical and microstructural studies often assume that the fibers are aligned in-plane with the tissues. This is convenient as it removes the out-of-plane components and allows focusing on the 2D maps of in-plane fiber organization that are often quite complex. The simplification, however, risks missing potentially important aspects of the tissue architecture and mechanics. In the cornea, for instance, fibers with high in-depth inclination have been shown to be mechanically important. Outside the cornea, the in-depth fiber orientations have not been characterized, preventing a deeper understanding of their potential roles. Our goal was to characterize in-depth collagen fiber organization over the whole corneoscleral shell. Seven sheep whole-globe axial sections from eyes fixed at an IOP of 50 mmHg were imaged using polarized light microscopy to measure collagen fiber orientations and density. In-depth fiber orientation distributions and anisotropy (degree of fiber alignment) accounting for fiber density were quantified over the whole sclera and in 15 regions: central cornea, peripheral cornea, limbus, anterior equator, equator, posterior equator, posterior sclera and peripapillary sclera on both nasal and temporal sides. Orientation distributions were fitted using a combination of a uniform distribution and a sum of π-periodic von Mises distributions, each with three parameters: primary orientation μ, fiber concentration factor k, and weighting factor a. To study the features of fibers that are not in-plane, i.e., fiber inclination, we quantified the percentage of inclined fibers and the range of inclination angles (half width at half maximum of inclination angle distribution). Our measurements showed that the fibers were not uniformly in-plane but exhibited instead a wide range of in-depth orientations, with fibers significantly more aligned in-plane in the anterior parts of the globe. We found that fitting the orientation distributions required between one and three π-periodic von Mises distributions with different primary orientations and fiber concentration factors. Regions of the posterior globe, particularly on the temporal side, had a larger percentage of inclined fibers and a larger range of inclination angles than anterior and equatorial regions. Variations of orientation distributions and anisotropies may imply varying out-of-plane tissue mechanical properties around the eye globe. Out-of-plane fibers could indicate fiber interweaving, not necessarily long, inclined fibers. Effects of small-scale fiber undulations, or crimp, were minimized by using tissues from eyes at high IOPs. These fiber features also play a role in tissue stiffness and stability and are therefore also important experimental information.

A noninvasive clinical method to measure in vivo mechanical strains of the lamina cribrosa by optical coherence tomography

Objective

To measure mechanical strain of the lamina cribrosa (LC) after intraocular pressure (IOP) change produced one week after a change in glaucoma medication.

Design

Cohort study.

Participants

Adult glaucoma patients (23 eyes, 15 patients) prescribed a change in IOP-lowering medication.

Intervention

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

Main Outcomes

Deformation calculated by digital volume correlation of OCT scans of the LC before and after IOP lowering by medication.

Results

Among 23 eyes, 17 eyes of 12 persons had IOP lowering ≥ 3 mmHg (reduced IOP group) with tensile anterior-posterior E zz strain = 1.0% ± 1.1% (p = 0.003) and compressive radial strain ( E rr ) = -0.3% ± 0.5% (p=0.012; random effects models accounting inclusion of both eyes in some persons). Maximum in-plane principal (tensile) strain and maximum shear strain in the reduced IOP group were: E max = 1.7% ± 1.0% and Γ max = 1.4% ± 0.7%, respectively (both p<0.0001 versus zero). Reduced IOP group strains E max and Γ max were significantly larger with greater %IOP decrease (<0.0001, <0.0001). The compliance of the E zz , E max , and Γ max strain response, defined as strain normalized by the IOP decrease, were larger with more abnormal perimetric mean deviation or visual field index values (all p≥0.02). Strains were unrelated to age (all p≥0.088). In reduced IOP eyes, mean LC anterior border posterior movement was only 2.05 μm posteriorly (p=0.052) and not related to % IOP change (p=0.94, random effects models). Only E rr was significantly related to ALD change, becoming more negative with greater posterior LC border change (p=0.015).

Conclusion

LC mechanical strains can be effectively measured by changes in eye drop medication using OCT and are related to degree of visual function loss in glaucoma.

Trial Registration

ClinicalTrials.gov Identifier: NCT03267849.

How Myopia and Glaucoma Influence the Biomechanical Susceptibility of the Optic Nerve Head

PURPOSE

The purpose of this study was to assess optic nerve head (ONH) deformations following acute intraocular pressure (IOP) elevations and horizontal eye movements in control eyes, highly myopic (HM) eyes, HM eyes with glaucoma (HMG), and eyes with pathologic myopia (PM) alone or PM with staphyloma (PM + S).

METHODS

We studied 282 eyes, comprising of 99 controls (between +2.75 and -2.75 diopters), 51 HM (< -5 diopters), 35 HMG, 21 PM, and 75 PM + S eyes. For each eye, we imaged the ONH using spectral-domain optical coherence tomography (OCT) under the following conditions: (1) primary gaze, (2) 20 degrees adduction, (3) 20 degrees abduction, and (4) primary gaze with acute IOP elevation (to ∼35 mm Hg) achieved through ophthalmodynamometry. We then computed IOP- and gaze-induced ONH displacements and effective strains. Effective strains were compared across groups.

RESULTS

Under IOP elevation, we found that HM eyes exhibited significantly lower strains (3.9 ± 2.4%) than PM eyes (6.9 ± 5.0%, P < 0.001), HMG eyes (4.7 ± 1.8%, P = 0.04), and PM + S eyes (7.0 ± 5.2%, P < 0.001). Under adduction, we found that HM eyes exhibited significantly lower strains (4.8% ± 2.7%) than PM + S eyes (6.0 ± 3.1%, P = 0.02). We also found that eyes with higher axial length were associated with higher strains.

CONCLUSIONS

Our study revealed that eyes with HMG experienced significantly greater strains under IOP compared to eyes with HM. Furthermore, eyes with PM + S had the highest strains on the ONH of all groups.

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.

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

Purpose

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

Methods

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

Results

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

Conclusions

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

High-frequency ultrasound point-of-care device to quantify myopia-induced microstructural changes in the anterior sclera

PURPOSE

To develop a point-of-care (POC) device using high-frequency ultrasound (US) for evaluating microstructural changes in the anterior sclera associated with myopia.

METHODS

The proposed POC device must satisfy four primary requirements for effective clinical use: the measurement component is handheld; the software must be simple and provide real-time feedback; patient safety and health data security requirements set forth by relevant governing bodies must be satisfied and the measurement data must have sufficient signal-to-noise ratio (SNR) and repeatability. Radiofrequency (RF) echo data acquired by the POC device will be processed using our quantitative US methods to characterise tissue microstructure and biomechanical properties.

RESULTS

All stated requirements have been met in the developed POC device. The high-frequency transducer is housed in a custom, 3D-printed, pen-like holder that allows for easy measurements of the anterior sclera. Custom software provides a simple interface for data acquisition, real-time data display and secure data storage. Exposimetry measurements of the US pressure field indicate device compliance with United States Food and Drug Administration limits for ophthalmic US. In vivo measurements on a volunteer suggest the RF data SNR and acquisition consistency are suitable for quantitative analysis.

CONCLUSIONS

A fully functioning POC device using high-frequency US has been created for evaluating the microstructure of the anterior sclera. Planned studies using the POC device to scan the eyes of myopia patients will help clarify how the anterior sclera microstructure may be affected by myopia. If effective, this portable, inexpensive and user-friendly system could be an important part of routine eye examinations.

Corneal Hysteresis for the Diagnosis of Glaucoma and Assessment of Progression Risk: A Report by the American Academy of Ophthalmology

PURPOSE

To review the current published literature on the utility of corneal hysteresis (CH) to assist the clinician in the diagnosis of glaucoma or in the assessment of risk for disease progression in existing glaucoma patients.

METHODS

Searches of the peer-reviewed literature in the PubMed database were performed through July 2022. The abstracts of 423 identified articles were examined to exclude reviews and non-English articles. After inclusion and exclusion criteria were applied, 19 articles were selected, and the panel methodologist rated them for level of evidence. Eight articles were rated level I, and 5 articles were rated level II. The 6 articles rated level III were excluded.

RESULTS

Corneal hysteresis is lower in patients with primary open-angle glaucoma, primary angle-closure glaucoma, pseudoexfoliative glaucoma, and pseudoexfoliation syndrome compared with normal subjects. Interpretation of low CH in patients with high intraocular pressure (IOP) or on topical hypotensive medications is complicated by the influence of these parameters on CH measurements. However, CH is also lower in treatment-naïve, normal-tension glaucoma patients compared with normal subjects who have a similar IOP. In addition, lower CH is associated with an increased risk of progression of glaucoma based on visual fields or structural markers in open-angle glaucoma patients, including those with apparently well-controlled IOP.

CONCLUSIONS

Corneal hysteresis is lower in glaucoma patients compared with normal subjects, and lower CH is associated with an increased risk of disease progression. However, a causal relationship remains to be demonstrated. Nevertheless, measurement of CH complements current structural and functional assessments in determining disease risk in glaucoma suspects and patients.

FINANCIAL DISCLOSURE(S)

Proprietary or commercial disclosure may be found after the references.

Typical localised element-specific finite element anterior eye model

Purpose

The study presents an averaged anterior eye geometry model combined with a localised material model that is straightforward, appropriate and amenable for implementation in finite element (FE) modelling.

Methods

Both right and left eye profile data of 118 subjects (63 females and 55 males) aged 22-67 years (38.5 ± 7.6) were used to build an averaged geometry model. Parametric representation of the averaged geometry model was achieved through two polynomials dividing the eye into three smoothly connected volumes. This study utilised the collagen microstructure x-ray data of 6 ex-vivo healthy human eyes, 3 right eyes and 3 left eyes in pairs from 3 donors, 1 male and 2 females aged between 60 and 80 years, to build a localised element-specific material model for the eye.

Results

Fitting the cornea and the posterior sclera sections to a 5th-order Zernike polynomial resulted in 21 coefficients. The averaged anterior eye geometry model recorded a limbus tangent angle of 37° at a radius of 6.6 mm from the corneal apex. In terms of material models, the difference between the stresses generated in the inflation simulation up to 15 mmHg in the ring-segmented material model and localised element-specific material model were significantly different (p < 0.001) with the ring-segmented material model recording average Von-Mises stress 0.0168 ± 0.0046 MPa and the localised element-specific material model recording average Von-Mises stress 0.0144 ± 0.0025 MPa.

Conclusions

The study illustrates an averaged geometry model of the anterior human eye that is easy to generate through two parametric equations. This model is combined with a localised material model that can be used either parametrically through a Zernike fitted polynomial or non-parametrically as a function of the azimuth angle and the elevation angle of the eye globe. Both averaged geometry and localised material models were built in a way that makes them easy to implement in FE analysis without additional computation cost compared to the limbal discontinuity so-called idealised eye geometry model or ring-segmented material model.

The Effects of Negative Periocular Pressure on Biomechanics of the Optic Nerve Head and Cornea: A Computational Modeling Study

Purpose

The purpose of this study was to evaluate the effects of negative periocular pressure (NPP), and concomitant intraocular pressure (IOP) lowering, on the biomechanics of the optic nerve head (ONH) and cornea.

Methods

We developed a validated finite element (FE) model of the eye to compute tissue biomechanical strains induced in response to NPP delivered using the Multi-Pressure Dial (MPD) system. The model was informed by clinical measurements of IOP lowering and was based on published tissue properties. We also conducted sensitivity analyses by changing pressure loads and tissue properties.

Results

Application of -7.9 mmHg NPP decreased strain magnitudes in the ONH by c. 50% whereas increasing corneal strain magnitudes by c. 25%. Comparatively, a similar increase in corneal strain was predicted to occur due to an increase in IOP of 4 mmHg. Sensitivity studies indicated that NPP lowers strain in the ONH by reducing IOP and that these effects persisted over a range of tissue stiffnesses and spatial distributions of NPP.

Conclusions

NPP is predicted to considerably decrease ONH strain magnitudes. It also increases corneal strain but to an extent expected to be clinically insignificant. Thus, using NPP to lower IOP and hence decrease ONH mechanical strain is likely biomechanically beneficial for patients with glaucoma.

Translational Relevance

This study provides the first description of how NPP affects ONH biomechanics and explains the underlying mechanism of ONH strain reduction. It complements current empirical knowledge about the MPD system and guides future studies of NPP as a treatment for glaucoma.

Differing Associations between Optic Nerve Head Strains and Visual Field Loss in Patients with Normal- and High-Tension Glaucoma

PURPOSE

To study the associations between optic nerve head (ONH) strains under intraocular pressure (IOP) elevation with retinal sensitivity in patients with glaucoma.

DESIGN

Clinic-based cross-sectional study.

PARTICIPANTS

Two hundred twenty-nine patients with primary open-angle glaucoma (subdivided into 115 patients with high-tension glaucoma [HTG] and 114 patients with normal-tension glaucoma [NTG]).

METHODS

For 1 eye of each patient, we imaged the ONH using spectral-domain OCT under the following conditions: (1) primary gaze and (2) primary gaze with acute IOP elevation (to approximately 35 mmHg) achieved through ophthalmodynamometry. A 3-dimensional strain-mapping algorithm was applied to quantify IOP-induced ONH tissue strain (i.e., deformation) in each ONH. Strains in the prelaminar tissue (PLT), the retina, the choroid, the sclera, and the lamina cribrosa (LC) were associated (using linear regression) with measures of retinal sensitivity from the 24-2 Humphrey visual field test (Carl Zeiss Meditec). This was performed globally, then locally according to a previously published regionalization scheme.

MAIN OUTCOME MEASURES

Associations between ONH strains and values of retinal sensitivity from visual field testing.

RESULTS

For patients with HTG, we found (1) significant negative linear associations between ONH strains and retinal sensitivity (P < 0.001; on average, a 1% increase in ONH strains corresponded to a decrease in retinal sensitivity of 1.1 decibels [dB]), (2) that high-strain regions colocalized with anatomically mapped regions of high visual field loss, and (3) that the strongest negative associations were observed in the superior region and in the PLT. In contrast, for patients with NTG, no significant associations between strains and retinal sensitivity were observed except in the superotemporal region of the LC.

CONCLUSIONS

We found significant negative associations between IOP-induced ONH strains and retinal sensitivity in a relatively large glaucoma cohort. Specifically, patients with HTG who experienced higher ONH strains were more likely to exhibit lower retinal sensitivities. Interestingly, this trend in general was less pronounced in patients with NTG, which could suggest a distinct pathophysiologic relationship between the two glaucoma subtypes.

Numerical simulation of viscoelastic fluid–structure interaction benchmarks and their application to the human eye

We present a numerical solution method for time-dependent viscoelastic fluid–structure interaction employing the arbitrary Lagrangian Eulerian framework. The derived monolithic variational formulation is discretized in time using the shifted Crank–Nicolson scheme and in space using the finite element method. For the linearisation we employ Newton’s method with exact Jacobians. The viscoelastic fluid is modelled either using the Oldroyd-B or a Burgers-type model. The elastic structures are non-linear hyperelastic materials. We validate the implementation on benchmark problems and numerically analyse the convergence for global mesh refinement and adaptive mesh refinement using the dual-weighted residual method. Furthermore we numerically analyse the influence of the viscoelasticity of the fluid on typical goal functionals like the drag, the lift and the displacement. The derived numerical solution method is applied to ophthalmology where we analyse the interaction of the viscoelastic vitreous with its surrounding elastic structures.

Bioengineering Human Pluripotent Stem Cell-Derived Retinal Organoids and Optic Vesicle-Containing Brain Organoids for Ocular Diseases

Retinal organoids are three-dimensional (3D) structures derived from human pluripotent stem cells (hPSCs) that mimic the retina's spatial and temporal differentiation, making them useful as in vitro retinal development models. Retinal organoids can be assembled with brain organoids, the 3D self-assembled aggregates derived from hPSCs containing different cell types and cytoarchitectures that resemble the human embryonic brain. Recent studies have shown the development of optic cups in brain organoids. The cellular components of a developing optic vesicle-containing organoids include primitive corneal epithelial and lens-like cells, retinal pigment epithelia, retinal progenitor cells, axon-like projections, and electrically active neuronal networks. The importance of retinal organoids in ocular diseases such as age-related macular degeneration, Stargardt disease, retinitis pigmentosa, and diabetic retinopathy are described in this review. This review highlights current developments in retinal organoid techniques, and their applications in ocular conditions such as disease modeling, gene therapy, drug screening and development. In addition, recent advancements in utilizing extracellular vesicles secreted by retinal organoids for ocular disease treatments are summarized.

Ocular Pulse Amplitude Correlates With Ocular Rigidity at Native IOP Despite the Variability in Intraocular Pulse Volume With Each Heartbeat

Purpose

The purpose of this study was to assess ocular coat mechanical behavior using controlled ocular microvolumetric injections (MVI) of 15 µL of balanced salt solution (BSS) infused over 1 second into the anterior chamber (AC) via a syringe pump.

Methods

Intraocular pressure (IOP) was continuously recorded at 200 Hz with a validated implantable IOP telemetry system in 7 eyes of 7 male rhesus macaques (nonhuman primates [NHPs]) during 5 MVIs in a series at native (3 trials), 15 and 20 mm Hg baseline IOPs, repeated in 2 to 5 sessions at least 2 weeks apart. Ocular rigidity coefficients (K) and ocular pulse volume (PV) were calculated for each trial. Data were averaged across sessions within eyes; PV was analyzed with a three-level nested ANOVA, and parameter relationships were analyzed with Pearson Correlation and linear regression.

Results

After MVI at native baseline IOP of 10.4 ± 1.6 mm Hg, IOP increased by 9.1 ± 2.8 mm Hg (∆IOP) at a 9.6 ± 2.7 mm Hg/s slope, ocular pulse amplitude (OPA) was 0.70 ± 0.13 mm Hg on average; the average K was 0.042 ± 0.010 µL-1 and average PV was 1.16 ± 0.43 µL. PV varied significantly between trials, days, and eyes (P ≤ 0.05). OPA was significantly correlated with K at native IOP: Pearson coefficients ranged from 0.71 to 0.83 (P ≤ 0.05) and R2 ranged from 0.50 to 0.69 (P ≤ 0.05) during the first trial.

Conclusions

The MVI-driven ∆IOP and slope can be used to assess ocular coat mechanical behavior and measure ocular rigidity.

Translational Relevance

Importantly, OPA at native IOP is correlated with ocular rigidity despite the significant variability in PV between heartbeats.

A minimal model of elastic instabilities in biological filament bundles

We present a model of a system of elastic fibres which exhibits complex, coupled, nonlinear deformations via a connecting elastic spring network. This model can capture physically observed deformations such as global buckling, pinching and internal collapse. We explore the transitions between these deformation modes numerically, using an energy minimization approach, highlighting how supported environments, or stiff outer sheath structures, favour internal structural collapse over global deformation. We then derive a novel analytic buckling criterion for the internal collapse of the system, a mode of structural collapse pertinent in many biological filament bundles such as the optic nerve bundle and microtubule bundles involved in cell abscission.

Materials for Implantable Surface Electrode Arrays: Current Status and Future Directions

Surface electrode arrays are mainly fabricated from rigid or elastic materials, and precisely manipulated ductile metal films, which offer limited stretchability. However, the living tissues to which they are applied are nonlinear viscoelastic materials, which can undergo significant mechanical deformation in dynamic biological environments. Further, the same arrays and compositions are often repurposed for vastly different tissues rather than optimizing the materials and mechanical properties of the implant for the target application. By first characterizing the desired biological environment, and then designing a technology for a particular organ, surface electrode arrays may be more conformable, and offer better interfaces to tissues while causing less damage. Here, the various materials used in each component of a surface electrode array are first reviewed, and then electrically active implants in three specific biological systems, the nervous system, the muscular system, and skin, are described. Finally, the fabrication of next-generation surface arrays that overcome current limitations is discussed.

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.

Method for the biomechanical analysis of aqueous veins and perilimbal sclera by three-dimensional photoacoustic imaging and strain field calculation

A method motivated by the eye's aqueous veins is described for the imaging and strain calculation within soft biological tissues. A challenge to the investigation of the biomechanics of the aqueous vein-perilimbal sclera tissue complex is resolution of tissue deformations as a function of intraocular pressure and the subsequent calculation of strain (a normalized measure of deformation). The method involves perfusion of the eye with a contrast agent during conduction of non-invasive, optical resolution photoacoustic microscopy. This imaging technique permits three-dimensional displacement measurements of tracked points on the inner walls of the veins which are used in a finite element model to determine the corresponding strains. The methods are validated against two standard strain measurement methods. Representative porcine globe perfusion experiments are presented that demonstrate the power of the method to determine complex strain fields in the veins dependent on intraocular pressure as well as vein anatomy. In these cases, veins are observed to move radially outward during increases in intraocular pressure and to possess significant spatial strain variation, possibly influenced by their branching patterns. To the authors' knowledge, these are the only such quantitative, data driven, calculations of the aqueous vein strains available in the open literature.

Estimation of scleral mechanical properties from air-puff optical coherence tomography

We introduce a method to estimate the biomechanical properties of the porcine sclera in intact eye globes ex vivo, using optical coherence tomography that is coupled with an air-puff excitation source, and inverse optimization techniques based on finite element modeling. Air-puff induced tissue deformation was determined at seven different locations on the ocular globe, and the maximum apex deformation, the deformation velocity, and the arc-length during deformation were quantified. In the sclera, the experimental maximum deformation amplitude and the corresponding arc length were dependent on the location of air-puff excitation. The normalized temporal deformation profile of the sclera was distinct from that in the cornea, but similar in all tested scleral locations, suggesting that this profile is independent of variations in scleral thickness. Inverse optimization techniques showed that the estimated scleral elastic modulus ranged from 1.84 ± 0.30 MPa (equatorial inferior) to 6.04 ± 2.11 MPa (equatorial temporal). The use of scleral air-puff imaging holds promise for non-invasively investigating the structural changes in the sclera associated with myopia and glaucoma, and for monitoring potential modulation of scleral stiffness in disease or treatment.

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.

Histologic validation of optical coherence tomography-based three-dimensional morphometric measurements of the human optic nerve head: Methodology and preliminary results

PURPOSE

To compare the three-dimensional (3D) morphology of the deep load-bearing structures of the human optic nerve head (ONH) as revealed in vivo by spectral domain optical coherence tomography (SDOCT) with ex vivo quantitative 3D histology.

METHODS

SDOCT imaging of the ONH was performed in six eyes from three brain-dead organ donors on life-support equipment awaiting organ procurement (in vivo conditions). Following organ procurement (ex vivo conditions), the eyes were enucleated and underwent a pars plana vitrectomy followed by pressurization to physiologic IOP and immersion fixation. Ex vivo ONH morphology was obtained from high-fidelity episcopic fluorescent 3D reconstruction. Morphologic parameters of the observed ONH canal geometry and peripapillary choroid, as well as the shape, visibility and depth of the lamina cribrosa were compared between ex vivo and in vivo measurements using custom software to align, scale, and manually delineate the different regions of the ONH.

RESULTS

There was significant correspondence between in vivo and ex vivo measurements of the depth and shape of the lamina cribrosa, along with the size and shape of Bruch's membrane opening (BMO) and anterior scleral canal opening (ASCO). Weaker correspondence was observed for choroidal thickness; as expected, a thinner choroid was seen ex vivo due to loss of blood volume upon enucleation (-79.9%, p < 0.001). In addition, the lamina was shallower (-32.3%, p = 0.0019) and BMO was smaller ex vivo (-3.38%, p = 0.026), suggesting post mortem shrinkage of the fixed tissue. On average, while highly variable, only 31% of the anterior laminar surface was visible in vivo with SDOCT (p < 0.001).

CONCLUSIONS

Morphologic parameters by SDOCT imaging of the deep ONH showed promising correspondence to histology metrics. Small but significant shrinkage artifact, along with large effects of exsanguination of the choroid, was seen in the ex vivo reconstructions of fixed tissues that may impact the quantification of ex vivo histoarchitecture, and this should be considered when developing models and biomarkers based on ex vivo imaging of fixed tissue. Lack of visibly of most of the lamina surface in SDOCT images is an important limitation to metrics and biomarkers based on in vivo images of the ONH deep tissues.

Automated MRI-based quantification of posterior ocular globe flattening and recovery after long-duration spaceflight

BACKGROUND/OBJECTIVES

Spaceflight associated neuro-ocular syndrome (SANS), a health risk related to long-duration spaceflight, is hypothesized to result from a headward fluid shift that occurs with the loss of hydrostatic pressure gradients in weightlessness. Shifts in the vascular and cerebrospinal fluid compartments alter the mechanical forces at the posterior eye and lead to flattening of the posterior ocular globe. The goal of the present study was to develop a method to quantify globe flattening observed by magnetic resonance imaging after spaceflight.

SUBJECTS/METHODS

Volumetric displacement of the posterior globe was quantified in 10 astronauts at 5 time points after spaceflight missions of ~6 months.

RESULTS

Mean globe volumetric displacement was 9.88 mm³ (95% CI 4.56-15.19 mm³, p < 0.001) on the first day of assessment after the mission (R[return]+ 1 day); 9.00 mm³ (95% CI 3.73-14.27 mm³, p = 0.001) at R + 30 days; 6.53 mm³ (95% CI 1.24-11.83 mm³, p < 0.05) at R + 90 days; 4.45 mm³ (95% CI -0.96 to 9.86 mm³, p = 0.12) at R + 180 days; and 7.21 mm³ (95% CI 1.82-12.60 mm³, p < 0.01) at R + 360 days.

CONCLUSIONS

There was a consistent inward displacement of the globe at the optic nerve, which had only partially resolved 1 year after landing. More pronounced globe flattening has been observed in previous studies of astronauts; however, those observations lacked quantitative measures and were subjective in nature. The novel automated method described here allows for detailed quantification of structural changes in the posterior globe that may lead to an improved understanding of SANS.

Dynamic Optical Coherence Elastography of the Anterior Eye: Understanding the Biomechanics of the Limbus

Purpose

Currently, the biomechanical properties of the corneo-scleral limbus when the eye-globe deforms are largely unknown. The purpose of this study is to evaluate changes in elasticity of the cornea, sclera, and limbus when subjected to different intraocular pressures (IOP) using wave-based optical coherence elastography (OCE). Special attention was given to the elasticity changes of the limbal region with respect to the elasticity variations in the neighboring corneal and scleral regions.

Methods

Continuous harmonic elastic waves (800 Hz) were mechanically induced in the sclera near the corneo-sclera limbus of in situ porcine eye-globes (n = 8). Wave propagation was imaged using a phase-sensitive optical coherence tomography system (PhS-OCT). The eyes were subjected to five different IOP-levels (10, 15, 20, 30, and 40 mm Hg), and spatially distributed propagation velocities were calculated along corneal, limbal, and scleral regions. Finite element analysis (FEA) of the same regions under the same excitation conditions were conducted for further validation of results.

Results

FEA demonstrated that the stiffness of the heterogeneous cornea-limbus-sclera transition can be characterized by phase velocity measurements of the elastic waves produced at 800 Hz in the anterior eye. Experimental results revealed that the wave speed in the limbus (cL = 6.5 m/s) is between the cornea (cc = 2.9 m/s) and sclera (cs = 10.0 m/s) at a physiological IOP level (15 mm Hg) and rapidly increases as the IOP level is increased, even surpassing the wave speed in the sclera. Finally, the change in elastic wave speed in the limbus (ΔcL∼18.5 m/s) was greater than in the cornea (Δcc ∼12.6 m/s) and sclera (Δcs∼8.1 m/s) for the same change in IOP.

Conclusions

We demonstrated that wave-based OCE can be utilized to assess limbus biomechanical properties. Moreover, experimental evidence showed that the corneo-scleral limbus is highly nonlinear compared to the cornea and sclera when the eye-globe is deformed by an increase of IOP. This may suggest that the limbus has enough structural flexibility to stabilize anterior eye shape during IOP changes.

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.

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.

Investigation of the Optic Nerve Head Morphology Influence to the Optic Nerve Head Biomechanics - Patient Specific Model

Glaucoma is associated with damage and death of optic nerve fibers within the Lamina Cribrosa (LC) region of the Optic Nerve Head. The pathogenesis of the disease is unclear, and the anterior LC surface morphology of different individuals can be one of the possible contributor of glaucoma development and progression. The current study evaluates the relationship between the LC surface curvature and distribution of shear stresses on the LC surface. The patient-specific reconstructed ocular model was developed and analyzed in a finite element analysis software. In addition, the effect of elastic modulus of both sclera and LC on the shear stress was examined. Results showed that there is a correlation between the shear stress distribution and the curvature of the anterior LC surface. This finding highlights the potential significance of the LC morphology on the distribution of LC shear stress and require further investigation.

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.

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

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

In Vivo Measurements of Prelamina and Lamina Cribrosa Biomechanical Properties in Humans

Purpose

To develop and use a custom virtual fields method (VFM) to assess the biomechanical properties of human prelamina and lamina cribrosa (LC) in vivo.

Methods

Clinical data of 20 healthy, 20 ocular hypertensive (OHT), 20 primary open-angle glaucoma, and 16 primary angle-closure glaucoma eyes were analyzed. For each eye, the intraocular pressure (IOP) and optical coherence tomography (OCT) images of the optic nerve head (ONH) were acquired at the normal state and after acute IOP elevation. The IOP-induced deformation of the ONH was obtained from the OCT volumes using a three-dimensional tracking algorithm and fed into the VFM to extract the biomechanical properties of the prelamina and the LC in vivo. Statistical measurements and P values from the Mann-Whitney-Wilcoxon tests were reported.

Results

The average shear moduli of the prelamina and the LC were 64.2 ± 36.1 kPa and 73.1 ± 46.9 kPa, respectively. The shear moduli of the prelamina of healthy subjects were significantly lower than those of the OHT subjects. Comparisons between healthy and glaucoma subjects could not be made robustly due to a small sample size.

Conclusions

We have developed a methodology to assess the biomechanical properties of human ONH tissues in vivo and provide preliminary comparisons in healthy and OHT subjects. Our proposed methodology may be of interest for glaucoma management.

Integration of polarized spatial frequency domain imaging (pSFDI) with a biaxial mechanical testing system for quantification of load-dependent collagen architecture in soft collagenous tissues

Collagen fiber networks provide the structural strength of tissues, such as tendons, skin and arteries. Quantifying the fiber architecture in response to mechanical loads is essential towards a better understanding of the tissue-level mechanical behaviors, especially in assessing disease-driven functional changes. To enable novel investigations into these load-dependent fiber structures, a polarized spatial frequency domain imaging (pSFDI) device was developed and, for the first time, integrated with a biaxial mechanical testing system. The integrated instrument is capable of a wide-field quantification of the fiber orientation and the degree of optical anisotropy (DOA), representing the local degree of fiber alignment. The opto-mechanical instrument''s performance was assessed through uniaxial loading on tendon tissues with known collagen fiber microstructures. Our results revealed that the bulk fiber orientation angle of the tendon tissue changed minimally with loading (median ± 0.5*IQR of 52.7° ± 3.3° and 51.9° ± 3.3° under 0 and 3% longitudinal strains, respectively), whereas on a micro-scale, the fibers became better aligned with the direction of loading: the DOA (mean ± SD) increased from 0.149 ± 0.032 to 0.198 ± 0.056 under 0 and 3% longitudinal strains, respectively, p < 0.001. The integrated instrument was further applied to study two representative mitral valve anterior leaflet (MVAL) tissues subjected to various biaxial loads. The fiber orientations within these representative MVAL tissue specimens demonstrated noticeable heterogeneity, with the local fiber orientations dependent upon the sample, the spatial and transmural locations, and the applied loading. Our results also showed that fibers were generally better aligned under equibiaxial (DOA = 0.089 ± 0.036) and circumferentially-dominant loading (DOA = 0.086 ± 0.037) than under the radially-dominant loading (DOA = 0.077 ± 0.034), indicating circumferential predisposition. These novel findings exemplify a deeper understanding of the load-dependent collagen fiber microstructures obtained through the use of the integrated opto-mechanical instrument. STATEMENT OF SIGNIFICANCE: In this study, a novel quantitative opto-mechanical system was developed by combining a polarized Spatial Frequency Domain Imaging (pSFDI) device with a biaxial mechanical tester. The integrated system was used to quantify the load-dependent collagen fiber microstructures in representative tendon and mitral valve anterior leaflet (MVAL) tissues. Our results revealed that MVAL's fiber architectures exhibited load-dependent spatial and transmural heterogeneities, suggesting further microstructural complexity than previously reported in heart valve tissues. These novel findings were possible through the system's ability to, for the first time, capture the load-dependent collagen architecture in the mitral valve anterior leaflet tissue over a wide field of view (e.g., 10 × 10 mm for the MVAL tissue specimens). Such capabilities afford unique future opportunities to improve patient outcomes through concurrent mechanical and microstructural assessments of healthy and diseased tissues in conditions such as heart valve regurgitation and calcification.

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.

Full-field displacement measurement of corneoscleral shells by combining multi-camera speckle interferometry with 3D shape reconstruction

Changes in the biomechanical properties of the connective tissue of the eye occur with age and underlie the development of several ocular diseases, such as glaucoma, myopia, and keratoconus. The biomechanical dynamics of ocular connective tissue are measured by ex vivo inflation testing, in which intraocular pressure (IOP) is varied and optical methods are used to produce maps of corneal and scleral displacement. Current optical methods are limited by acquisition rate, occlusions, poor spatial resolution, and insufficient 3D mapping. We developed an interferometric optical method integrates four-camera electronic speckle pattern interferometry (ESPI) and a novel three-dimensional (3D) shape reconstruction process to measure shape and full-field mechanical deformations of corneal and scleral shells during ex vivo inflation testing. Each camera provides accurate measurements of the laser beam phase related to deformations of the specimen surface; a multi-view stereovision method generates the shape of the specimen and a functional form that links every pixel of a given camera to 3D points on the specimen's visible surface. In this way, dynamic deformations of the specimen are localized, with quantification of the time-dependent 3D displacements of the specimen at nanometric accuracy. The ESPI-3D system is suitable for analyzing scleral deformation and morphological changes caused by time-varying IOP.

Measurement of Ocular Compliance Using iPerfusion

The pressure-volume relationship of the eye is determined by the biomechanical properties of the corneoscleral shell and is classically characterised by Friedenwald's coefficient of ocular rigidity or, alternatively, by the ocular compliance (OC), defined as dV/dP. OC is important in any situation where the volume (V) or pressure (P) of the eye is perturbed, as occurs during several physiological and pathological processes. However, accurately measuring OC is challenging, particularly in rodents. We measured OC in 24 untreated enucleated eyes from 12 C57BL/6 mice using the iPerfusion system to apply controlled pressure steps, whilst measuring the time-varying flow rate into the eye. Pressure and flow data were analysed by a “Discrete Volume” (integrating the flow trace) and “Step Response” method (fitting an analytical solution to the pressure trace). OC evaluated at 13 mmHg was similar between the two methods (Step Response, 41 [37, 46] vs. Discrete Volume, 42 [37, 48] nl/mmHg; mean [95% CI]), although the Step Response Method yielded tighter confidence bounds on individual eyes. OC was tightly correlated between contralateral eyes (R² = 0.75, p = 0.0003). Following treatment with the cross-linking agent genipin, OC decreased by 40 [33, 47]% (p = 0.0001; N = 6, Step Response Method). Measuring OC provides a powerful tool to assess corneoscleral biomechanics in mice and other species.

Inflation experiments and inverse finite element modelling of posterior human sclera

The complexity of inverse finite element modelling methods used in ocular biomechanics research has significantly increased in recent years in order to produce material parameters that capture microscale tissue behaviour. This study presents a more accessible method for researchers to optimise sclera material parameters for use in finite element studies where macroscale sclera displacements are required. Five human donor sclerae aged between 36 and 72 years were subjected to cycles of internal pressure up to 61 mmHg using a custom-built inflation rig. Displacements were measured using a laser beam and two cameras through a digital image correlation algorithm. Specimen-specific finite element models incorporating regional thickness variation and sclera surface topography were divided into six circumferential regions. An inverse finite element procedure was used to optimise Ogden material parameters for each region. The maximum root mean squared (RMS) error between the numerical and experimental displacements within individual specimens was 17.5 µm. The optimised material parameters indicate a gradual reduction in material stiffness (as measured by the tangent modulus) from the equator to the posterior region at low-stress levels up to 0.005 MPa. The variation in stiffness between adjacent regions became gradually less apparent and statistically insignificant at higher stresses. The study demonstrated how inflation testing combined with inverse modelling could be used to effectively characterise regional material properties capable of reproducing global sclera displacements. The material properties were found to vary between specimens, and it is expected that age could be a contributing factor behind this variation.

A Subdomain Method for Mapping the Heterogeneous Mechanical Properties of the Human Posterior Sclera

Although strongly correlated with elevated intraocular pressure, primary open-angle glaucoma (POAG) occurs in normotensive eyes. Mechanical properties of the sclera around the optic nerve head (ONH) may play a role in this disparity. The purpose of this study is to present an automated inverse mechanics based approach to determine the distribution of heterogeneous mechanical properties of the human sclera as derived from its surface deformations arising from pressure inflation experiments. The scleral shell of a 78 year old European Descent male donor eye was utilized to demonstrate the method; the sclera was coated with a speckle pattern on the outer surface and was subjected to inflation pressures of 5, 15, 30, and 45 mmHg. The speckle pattern was imaged at each pressure, and a displacement field was calculated for each pressure step using a previously described sequential digital image correlation (S-DIC) technique. The fiber splay and fiber orientation of the sclera collagen were determined experimentally, and the thickness across the scleral globe was determined using micro CT images. The displacement field from the inflation test was used to calculate the strain and also used as an input for inverse mechanics to determine the heterogeneity of material properties. The scleral geometry was divided into subdomains using the first principal strain. The Holzapfel anisotropic material parameters of matrix and fiber stiffness were estimated within each individual subdomain using an inverse mechanics approach by minimizing the sum of the square of the residuals between the computational and experimental displacement fields. The mean and maximum error in displacement across all subdomains were 8.9 ± 3.0 μm and 13.2 μm, respectively. The full pressure-inflation forward mechanics experiment was done using subdomain-specific mechanical properties on the entire scleral surface. The proposed approach is effective in determining the distribution of heterogeneous mechanical properties of the human sclera in a user-independent manner. Our research group is currently utilizing this approach to better elucidate how scleral stiffness influences those at high risk for POAG.

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.

Estimating three-dimensional outflow and pressure gradients within the human eye

In this paper we set the previously reported pressure-dependent, ordinary differential equation outflow model by Smith and Gardiner for the human eye, into a new three-dimensional (3D) porous media outflow model of the eye, and calibrate model parameters using data reported in the literature. Assuming normal outflow through anterior pathways, we test the ability of 3D flow model to predict the pressure elevation with a silicone oil tamponade. Then assuming outflow across the retinal pigment epithelium is normal, we test the ability of the 3D model to predict the pressure elevation in Schwartz-Matsuo syndrome. For the first time we find the flow model can successfully model both conditions, which helps to build confidence in the validity and accuracy of the 3D pressure-dependent outflow model proposed here. We employ this flow model to estimate the translaminar pressure gradient within the optic nerve head of a normal eye in both the upright and supine postures, and during the day and at night. Based on a ratio of estimated and measured pressure gradients, we define a factor of safety against acute interruption of axonal transport at the laminar cribrosa. Using a completely independent method, based on the behaviour of dynein molecular motors, we compute the factor of safety against stalling the dynein molecule motors, and so compromising retrograde axonal transport. We show these two independent methods for estimating factors of safety agree reasonably well and appear to be consistent. Taken together, the new 3D pressure-dependent outflow model proves itself to capable of providing a useful modeling platform for analyzing eye behaviour in a variety of physiological and clinically useful contexts, including IOP elevation in Schwartz-Matsuo syndrome and with silicone oil tamponade, and potentially for risk assessment for optic glaucomatous neuropathy.

Analysis of X-ray scattering microstructure data for implementation in numerical simulations of ocular biomechanical behaviour

This study aimed to analyse microstructure data on the density and orientation of collagen fibrils in whole eye globes and to propose an effective method for the preparation of data for use in numerical simulations of the eye’s biomechanical performance. Wide-angle X-ray scattering was applied to seven healthy ex-vivo human eyes. Each eye was dissected into an anterior and a posterior cup, and radial incisions were used to flatten the tissue before microstructure characterisation. A method was developed to use the microstructure data obtained for the dissected tissue to build realistic 3D maps of fibril density and orientation covering the whole eye globe. At the central cornea, 61.5±2.3% of fibrils were aligned within 45° sectors surrounding the two orthogonal directions. In contrast, more than one-third of the total fibril content was concentrated along the circumferential direction at the limbus (37.0±2.4%) and around the optic nerve head (34.8±2.1%). The insertion locations of the four recti muscles exhibited a preference in the meridional direction near the equator (38.6±3.9%). There was also a significant difference in fibril density between the limbus and other regions (ratio = 1.91±0.45, p <0.01 at the central cornea and ratio = 0.80±0.21, p <0.01 at the posterior pole). Characterisation of collagen fibril density and orientation across the whole ocular surface has been possible but required the use of a technique that involved tissue dissection and hence caused tissue damage. The method presented in this paper aimed to minimise the effect of dissection on the quality of obtained data and was successful in identifying fibril distribution trends that were compatible with earlier studies, which concentrated on localised areas of the ocular globe.

Experimental myopia increases and scleral crosslinking using genipin inhibits cyclic softening in the tree shrew sclera

PURPOSE

Myopia progression is thought to involve biomechanical weakening of the sclera, which leads to irreversible deformations and axial elongation of the eye. Scleral crosslinking has been proposed as a potential treatment option for myopia control by strengthening the mechanically weakened sclera. The biomechanical mechanism by which the sclera weakens during myopia and strengthens after crosslinking is not fully understood. Here, we assess the effect of lens-induced myopia and exogenous crosslinking using genipin on the inelastic mechanical properties of the tree shrew sclera measured by cyclic tensile tests.

METHODS

Cyclic tensile tests were performed on 2-mm wide scleral strips at physiological loading conditions (50 cycles, 0-3.3 g, 30 s cycle^-1 ). Two scleral strips were obtained from each eye of juvenile tree shrews exposed to two different visual conditions: normal and 4 days of monocular -5 D lens wear to accelerate scleral remodelling and induce myopia. Scleral strips were mechanically tested at three alternative conditions: immediately after enucleation; after incubation in phosphate buffered saline (PBS) for 24 h at 37°C; and after incubation for 24 h in PBS supplemented with genipin at a low cytotoxicity concentration (0.25 mm). Cyclic softening was defined as the incremental strain increase from one cycle to the next.

RESULTS

-5D lens treatment significantly increased the cyclic softening response of the sclera when compared to contralateral control eyes (0.10% ± 0.029%, mean ± standard error, P = 0.037). Exogenous crosslinking of the lens treated sclera significantly decreased the cyclic softening response (-0.12% ± 0.014%, P = 2.2 × 10^-5 ). Contrary to all other groups, the genipin-cross-linked tissue did not exhibit cyclic softening significantly different from zero within the 50-cycle test.

CONCLUSIONS

Results indicated that cyclic tensile loading leads to an inelastic, cyclic softening of the juvenile tree shrew sclera. The softening rate increased during lens-induced myopia and was diminished after genipin crosslinking. This finding suggests that axial elongation in myopia may involve a biomechanical weakening mechanism that increased the cyclic softening response of the sclera, which was inhibited by scleral crosslinking using genipin.