Heartbeat-Induced Corneal Axial Displacement and Strain Measured by High Frequency Ultrasound Elastography in Human Volunteers

Keratectasia severity staging and progression assessment based on the biomechanical E-staging

Until recently, corneal topography has been the gold standard in detecting keratectasia and monitoring its progression. The recently introduced ABCD tomographic keratoconus staging system focuses on anterior ("A") and posterior ("B") radius of curvature, thinnest corneal thickness ("C"), best-corrected visual acuity with spectacles ("D") and is supplemented with the introduction of the biomechanical E-staging (BEST, "E"). The need for biomechanical staging arose from the fact of altered biomechanical characteristics of keratectasia in comparison to healthy corneas. Ectatic corneas usually exhibit a biomechanical weakening and greater deformation than healthy corneas when exposed to a biomechanical stressor such as a standardized air puff indentation as provided by the Corvis ST® (CST, Oculus, Wetzlar, Germany). The BEST is based on the linear term of the Corvis Biomechanical Index (CBI) and provides a biomechanical keratoconus severity staging and progression assessment within the CST software. This review traces the development of the BEST as an addition to the tomographic ABCD staging system and highlights its strengths and limitations when applied in daily practice for the detection, monitoring and progression assessment in keratectasia.

Biomechanical Correlations Between the Cornea and the Optic Nerve Head

Purpose

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

Methods

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

Results

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

Conclusions

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

Structural control of corneal transparency, refractive power and dynamics

The cornea needs to be transparent to visible light and precisely curved to provide the correct refractive power. Both properties are governed by its structure. Corneal transparency arises from constructive interference of visible light due to the relatively ordered arrangement of collagen fibrils in the corneal stroma. The arrangement is controlled by the negatively charged proteoglycans surrounding the fibrils. Small changes in fibril organisation can be tolerated but larger changes cause light scattering. Corneal keratocytes do not scatter light because their refractive index matches that of the surrounding matrix. When activated, however, they become fibroblasts that have a lower refractive index. Modelling shows that this change in refractive index significantly increases light scatter. At the microscopic level, the corneal stroma has a lamellar structure, the parallel collagen fibrils within each lamella making a large angle with those of adjacent lamellae. X-ray scattering has shown that the lamellae have preferred orientations in the human cornea: inferior-superior and nasal-temporal in the central cornea and circumferential at the limbus. The directions at the centre of the cornea may help withstand the pull of the extraocular muscles whereas the pseudo-circular arrangement at the limbus supports the change in curvature between the cornea and sclera. Elastic fibres are also present; in the limbus they contain fibrillin microfibrils surrounding an elastin core, whereas at the centre of the cornea, they exist as thin bundles of fibrillin-rich microfibrils. We present a model based on the structure described above that may explain how the cornea withstands repeated pressure changes due to the ocular pulse.

Measuring intraocular pressure with OCT: the first approach

The variability of corneal OCT speckle statistics is indirectly related to changes in corneal microstructure, which may be induced by intraocular pressure (IOP). A new approach is considered, which attempts to estimate IOP based on corneal speckle statistics in OCT images. An area (A) under trajectories of contrast ratio with respect to stromal depth was calculated. The proposed method was evaluated on OCT images from the ex-vivo study on porcine eyeballs and in-vivo study on human corneas. A statistically significant multivariate linear regression model was obtained from the ex-vivo study: IOP = 0.70 · A - 6.11, in which IOP was precisely controlled in the anterior chamber. The ex-vivo study showed good correlation between A and IOP (R = 0.628, at the least) whereas the in-vivo study showed poor correlation between A and clinical air-puff tonometry based estimates of IOP (R = 0.351, at the most), indicating substantial differences between the two studies. The results of the ex-vivo study show the potential for OCT speckle statistics to be utilized for measuring IOP using static corneal imaging that does not require corneal deformation. Nevertheless, further work is needed to validate this approach in living human corneas.

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.

Simultaneous Assessment of the Whole Eye Biomechanics Using Ultrasonic Elastography

OBJECTIVE

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

METHODS

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

RESULTS

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

CONCLUSION

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

SIGNIFICANCE

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

Corneal Biomechanics Losses Caused by Refractive Surgery

Recent advances, specifically in the understanding of the biomechanical properties of the cornea and its response to diseases and surgical interventions, have significantly improved the safety and surgical outcomes of corneal refractive surgery, whose popularity and demand continue to grow worldwide. However, iatrogenic keratectasia resulting from the deterioration in corneal biomechanics caused by surgical interventions, although rare, remains a global concern. On one hand, in vivo biomechanical evaluation, enabled by clinical imaging systems such as the ORA and the Corvis ST, has significantly improved the risk profiling of patients for iatrogenic keratectasia. That is despite the fact the biomechanical metrics provided by these systems are considered indicators of the cornea's overall stiffness rather than its intrinsic material properties. On the other hand, new surgical modalities including SMILE were introduced to offer superior biomechanical performance to LASIK, but this superiority could not be proven clinically, creating more myths than answers. The literature also includes sound evidence that tPRK provided the highest preservation of corneal biomechanics when compared to both LASIK and SMILE. The aim of this review is twofold; to discuss the importance of corneal biomechanical evaluation prior to refractive surgery, and to assess the current understanding of cornea's biomechanical deterioration caused by mainstream corneal refractive surgeries. The review has led to an observation that new imaging techniques, parameters and evaluation systems may be needed to reflect the true advantages of specific refractive techniques and when these advantages are significant enough to offer better protection against post-surgery complications.

High-frequency ultrasound detects biomechanical weakening in keratoconus with lower stiffness at higher grade

In vivo biomechanical characterization of the cornea remains a challenge. We have developed a high-frequency ultrasound elastography method, the ocular pulse elastography (OPE), to measure corneal axial displacement (CAD) induced by the ocular pulse. Here we compared CAD and a stiffness index derived from CAD between keratoconus patients and normal controls. We also explored the trend of these parameters with keratoconus grade. Twenty normal subjects and twenty keratoconus patients were recruited in this study. Corneal topography, tomography, intraocular pressure (IOP) and ocular pulse amplitude (OPA) were obtained in each measured eye. The cornea’s heartbeat-induced cyclic axial displacements were measured by high-frequency (50 MHz) ultrasound. A corneal stiffness index (CSI) was derived from CAD normalized against OPA. CAD and CSI were compared between normal and keratoconus groups, and across keratoconus grades. Keratoconus corneas had significantly greater CAD and lower CSI than normal controls (p’s<0.01). Both parameters correlated strongly with grade, in which CAD increased significantly (p = 0.002) and CSI decreased significantly (p = 0.011) with grade. These results suggested a biomechanical weakening in keratoconus which worsens at higher disease severity. This study also demonstrated the ability of high-frequency ultrasound elastography to provide a safe, quick, and accurate evaluation of the cornea’s biomechanical condition in vivo. The OPE-measured biomechanical metrics, when integrated with existing diagnostic criteria, may aid the decision-making in the early and definitive diagnosis and staging of keratoconus.

Multimodal Heartbeat and Compression Optical Coherence Elastography for Mapping Corneal Biomechanics

The biomechanical properties of the cornea have a profound influence on the health, structural integrity, and function of the eye. Understanding these properties may be critical for diagnosis and identifying disease pathogenesis. This work demonstrates how two different elastography techniques can be combined for a multimodal approach to measuring corneal biomechanical properties. Heartbeat optical coherence elastography (Hb-OCE) and compression OCE were performed simultaneously to measure the stiffness of the cornea in an in vivo rabbit model. Measurements were further performed after collagen crosslinking to demonstrate how the combined technique can be used to measure changes in corneal stiffness and map mechanical contrast. The results of this work further suggest that measurements from Hb-OCE and compression OCE are comparable, meaning that Hb-OCE and compression OCE may be used interchangeably despite distinct differences in both techniques.

Temporal phase evolution OCT for measurement of tissue deformation in the human retina in-vivo

We demonstrate the use of temporal phase evolution (TPE-) OCT methods to evaluate retinal tissue deformation in-vivo over time periods of several seconds. A custom built spectral domain (SD)-OCT system with an integrated retinal tracker, ensuring stable imaging with sub-speckle precision, was used for imaging. TPE-OCT measures and images phase differences between an initial reference B-scan and each of the subsequent B-scans of the evaluated temporal sequence. In order to demonstrate the precision and repeatability of the measurements, retinal nerve fiber (RNF) tissue deformations induced by retinal vessels pulsating with the heartbeat were analyzed in several healthy subjects. We show TPE maps (M-scans of phase evolution as a function of position along B-scan trace vs. time) of wrapped phase data and corresponding deformation maps in selected regions of the RNF layer (RNFL) over the course of several cardiac cycles. A reproducible phase pattern is seen at each heartbeat cycle for all imaged volunteers. RNF tissue deformations near arteries and veins up to ∼ 1.6 µm were obtained with an average precision for a single pixel of about 30 nm. Differences of motion induced by arteries and veins are also investigated.

Regional variation of corneal stromal deformation measured by high-frequency ultrasound elastography

The cornea's mechanical response to intraocular pressure elevations may alter in ectatic diseases such as keratoconus. Regional variations of mechanical deformation in normal and keratoconus eyes during intraocular pressure elevation have not been well-characterized. We applied a high-frequency ultrasound elastography technique to characterize the regional deformation of normal and keratoconus human corneas through the full thickness of corneal stroma. A cross-section centered at the corneal apex in 11 normal and 2 keratoconus human donor eyes was imaged with high-frequency ultrasound during whole globe inflation from 5 to 30 mmHg. An ultrasound speckle tracking algorithm was used to compute local tissue displacements. Radial, tangential, and shear strains were mapped across the imaged cross-section. Strains in the central (1 mm surrounding apex) and paracentral (1 to 4 mm from apex) regions were analyzed in both normal and keratoconus eyes. Additional regional analysis was performed in the eye with severe keratoconus presenting significant thinning and scarring. Our results showed that in normal corneas, the central region had significantly smaller tangential stretch than the paracentral region, and that within the central region, the magnitudes of radial and shear strains were significantly larger than that of tangential strain. The eye with mild keratoconus had similar shear strain but substantially larger radial strains than normal corneas, while the eye with severe keratoconus had similar overall strains as in normal eyes but marked regional heterogeneity and large strains in the cone region. These findings suggested regional variation of mechanical responses to intraocular pressure elevation in both normal and keratoconus corneas, and keratoconus appeared to be associated with mechanical weakening in the cone region, especially in resisting radial compression. Comprehensive characterization of radial, tangential, and shear strains through corneal stroma may provide new insights to understand the biomechanical alterations in keratoconus.