Purpose: This study aims to introduce and clinically validate a new algorithm that can determine the biomechanical properties of the human cornea in vivo. Methods: A parametric study was conducted involving representative finite element models of human ocular globes with wide ranges of geometries and material biomechanical behavior. The models were subjected to different levels of intraocular pressure (IOP) and the action of external air puff produced by a non-contact tonometer. Predictions of dynamic corneal response under air pressure were analyzed to develop an algorithm that can predict the cornea's material behavior. The algorithm was assessed using clinical data obtained from 480 healthy participants where its predictions of material behavior were tested against variations in central corneal thickness (CCT), IOP and age, and compared against those obtained in earlier studies on ex-vivo human ocular tissue. Results: The algorithm produced a material stiffness parameter (Stress-Strain Index or SSI) that showed no significant correlation with both CCT (p > 0.05) and IOP (p > 0.05), but was significantly correlated with age (p < 0.01). The stiffness estimates and their variation with age were also significantly correlated (p < 0.01) with stiffness estimates obtained earlier in studies on ex-vivo human tissue. Conclusions: The study introduced and validated a new method for estimating the in vivo biomechanical behavior of healthy corneal tissue. The method can aid optimization of procedures that interfere mechanically with the cornea such as refractive surgeries and introduction of corneal implants.

It is thought that corneal surface topography may be stabilized by the angular orientation of out-of plane lamellae that insert into the anterior limiting membrane. In this study, micro-focus X-ray scattering data were used to obtain quantitative information about lamellar inclination (with respect to the corneal surface) and the X-ray scatter intensity throughout the depth of the cornea from the centre to the temporal limbus. The average collagen inclination remained predominantly parallel to the tissue surface at all depths. However, in the central cornea, the spread of inclination angles was greatest in the anterior-most stroma (reflecting the increased lamellar interweaving in this region), and decreased with tissue depth; in the peripheral cornea inclination angles showed less variation throughout the tissue thickness. Inclination angles in the deeper stroma were generally higher in the peripheral cornea, suggesting the presence of more interweaving in the posterior stroma away from the central cornea. An increase in collagen X-ray scatter was identified in a region extending from the sclera anteriorly until about 2 mm from the corneal centre. This could arise from the presence of larger diameter fibrils, probably of scleral origin, which are known to exist in this region. Incorporation of this quantitative information into finite-element models will further improve the accuracy with which they can predict the biomechanical response of the cornea to pathology and refractive procedures.

Corneal biomechanics play a fundamental role in the genesis and progression of corneal pathologies, such as keratoconus; in corneal remodeling after corneal surgery; and in affecting the measurement accuracy of glaucoma biomarkers, such as the intraocular pressure (IOP). Air-puff induced corneal deformation imaging reveals information highlighting normal and pathological corneal response to a non-contact mechanical excitation. However, current commercial systems are limited to monitoring corneal deformation only on one corneal meridian. Here, we present a novel custom-developed swept-source optical coherence tomography (SSOCT) system, coupled with a collinear air-puff excitation, capable of acquiring dynamic corneal deformation on multiple meridians. Backed by numerical simulations of corneal deformations, we propose two different scan patterns, aided by low coil impedance galvanometric scan mirrors that permit an appropriate compromise between temporal and spatial sampling of the corneal deformation profiles. We customized the air-puff module to provide an unobstructed SSOCT field of view and different peak pressures, air-puff durations, and distances to the eye. We acquired multi-meridian corneal deformation profiles (a) in healthy human eyes in vivo, (b) in porcine eyes ex vivo under varying controlled IOP, and (c) in a keratoconus-mimicking porcine eye ex vivo. We detected deformation asymmetries, as predicted by numerical simulations, otherwise missed on a single meridian that will substantially aid in corneal biomechanics diagnostics and pathology screening.

Purpose. The organization of scleral collagen helps to determine the eye's biomechanical response to intraocular pressure (IOP), and may therefore be important in glaucoma. This study provides a quantitative assessment of changes in scleral collagen fibril organization in bead-induced murine experimental glaucoma. Methods. Wide-angle X-ray scattering was used to study the effect of bead-induced glaucoma on posterior scleral collagen organization in one eye of 12 CD1 mice, with untreated fellow eyes serving as controls. Three collagen parameters were measured: the local preferred fibril directions, the degree of collagen anisotropy, and the total fibrillar collagen content. Results. The mouse sclera featured a largely circumferential orientation of fibrillar collagen with respect to the optic nerve head canal. Localized alteration to fibril orientations was evident in the inferior peripapillary sclera of bead-treated eyes. Collagen anisotropy was significantly (P < 0.05) reduced in bead-treated eyes in the superior peripapillary (Treated: 43 ± 8%; Control: 49 ± 6%) and midposterior (Treated: 39 ± 4%; Control: 43 ± 4%) sclera, and in the peripapillary region overall (Treated: 43 ± 6%; Control: 47 ± 3%). No significant differences in total collagen content were found between groups. Conclusions. Spatial changes in collagen fibril anisotropy occur in the posterior sclera of mice with bead-induced chronic IOP elevation and axonal damage. These results support the idea that dynamic changes in scleral form and structure play a role in the development of experimental glaucoma in mice, and potentially in human glaucoma.

Purpose: To introduce a new method to map the mechanical stiffness of healthy and keratoconic corneas.Methods: Numerical modeling based on the finite element method was used to carry out inverse analysis of simulated healthy and keratoconic corneas to determine the regional variation of mechanical stiffness across the corneal surface based on established trends in collagen fibril distribution. The Stress–Strain Index (SSI), developed and validated in an earlier study and presented as a parameter that can estimate the overall stress–strain behavior of corneal tissue, was adopted in this research as a measure of corneal stiffness. The regional variation of SSI across the corneal surface was estimated using inverse analysis while referring to the common features of collagen fibrils’ distribution obtained from earlier x-ray scattering studies. Additionally, for keratoconic corneas, a method relating keratoconic cone features and cornea’s refractive power to the reduction in collagen fibril density inside the cone was implemented in the development of SSI maps. In addition to the simulated cases, the study also included two keratoconus cases, for which SSI maps were developed.Results: SSI values varied slightly across corneal surface in the simulated healthy eyes. In contrast, both simulated and clinical keratoconic corneas demonstrated substantial reductions in SSI values inside the cone. These SSI reductions depended on the extent of the disease and increased with more considerable simulated losses in fibril density in the cone area. SSI values and their regional variation showed little change with changes in IOP, corneal thickness, and curvature.Conclusion: SSI maps provide an estimation of the regional variation of biomechanical stiffness across the corneal surface. The maps could be particularly useful in keratoconic corneas, demonstrating the dependence of corneal biomechanical behavior on the tissue’s microstructure and offering a tool to fundamentally understand the mechanics of keratoconus progression in individual patients.

A quantitative map of collagen fibril orientation across the human eyeball coat, including both the cornea and the sclera, has been obtained using a combination of synchrotron wide-angle X-ray scattering (WAXS) and three-dimensional point mapping. A macromolecular crystallography beamline, in a custom-modified fibre diffraction setup, was used to record the 1.6 nm intermolecular equatorial reflection from fibrillar collagen at 0.5 mm spatial resolution across a flat-mounted human eyeball coat. Fibril orientation, derived as an average measure of the tissue thickness, was quantified by extraction of the azimuthal distribution of WAXS scatter intensity. Vector plots of preferential fibre orientation were remapped onto an idealized eyeball surface using a custom-built numerical algorithm, to obtain a three-dimensional representation of the collagen fibril architecture.