3D Characterization of corneal deformation using ultrasound speckle tracking

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.

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.

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.

IOP-induced regional displacements in the optic nerve head and correlation with peripapillary sclera thickness

Mechanical insult induced by intraocular pressure (IOP) is likely a driving force in the disease process of glaucoma. This study aimed to evaluate regional displacements in human optic nerve head (ONH) and peripapillary tissue (PPT) in response to acute IOP elevations, and their correlations with morphological characteristics of the posterior eye. Cross-sectional (2D) images of the ONH and PPT in 14 globes of 14 human donors were acquired with high-frequency ultrasound during whole globe inflation from 5 to 30 mm Hg. High-frequency ultrasound has a spatial resolution of tens of micrometers and is capable of imaging through the ONH and PPT thickness. Tissue displacements were calculated using a correlation-based speckle tracking algorithm for a dense matrix of kernels covering the 2D imaging plane. The ONH was manually segmented in the ultrasound B-mode images acquired at 5 mmHg based on echogenicity. The lamina cribrosa (LC) boundaries were visible in eight of the fourteen eyes and the LC region was segmented using a semi-automated superpixel-based method. The ONH had larger radial displacement than the PPT in all tested eyes and the difference increased with increasing IOP. A significant negative correlation was found between ONH-PPT displacement difference and PPT thickness (p < 0.05), while no significant correlations were found between ONH-PPT displacement difference and other morphological parameters including PPT radius of curvature, scleral canal size, LC thickness and anterior LC surface depth. Within the ONH, the radial displacement decreased in the region anterior to and across LC but not in the region posterior to LC. Finite element models using simplified geometry and material properties confirmed the role of LC in reducing the overall ONH radial displacements, but did not predict the displacement gradient change observed experimentally. These results suggested that a thinner PPT may be associated with a larger relative posterior motion of the ONH with respect to the surrounding PPT and the LC may play a major role in preventing excessive posterior displacement of ONH during acute IOP elevations.

Mechanical Stability of Cryopreserved Split-Thickness Tectonic Corneal Grafts

PURPOSE

To evaluate the influence of cryopreservation on the pressure-strain relationship of microkeratome dissected anterior stromal grafts (ASGs).

METHODS

Partial thickness ASGs were created from 7 pairs of human corneas and randomized to immediate grafting or grafting after 3 months of cryopreservation at -80°C into a whole globe ex vivo corneal perforation model. High frequency ultrasound speckle tracking was used to calculate the cross-sectional axial and lateral strains in each graft at increasing intraocular pressure (IOP) from 5 to 30 mm Hg. The mean axial and lateral strains were compared between the paired groups.

RESULTS

The mean axial and lateral strains were not significantly different between the cryopreserved and noncryopreserved ASGs. The mean lateral strains at 30 mm Hg in the noncryopreserved and cryopreserved grafts were 2.4% ± 2.1% and 1.4% ± 0.7% (P = 0.294), respectively. The mean axial strains at 30 mm Hg in the noncryopreserved and cryopreserved grafts were -7.8% ± 3.3% and -5.5% ± 3.0% (P = 0.198), respectively. A linear pressure-strain relationship was found for all grafts at physiologic IOP.

CONCLUSIONS

ASGs cryopreserved at -80°C maintain their IOP-strain relationship compared with noncryopreserved ASGs at physiologic pressures, supporting the potential use of cryopreserved human corneal stroma for patch grafting procedures.

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.

Ocular Pulse Elastography: Imaging Corneal Biomechanical Responses to Simulated Ocular Pulse Using Ultrasound

Purpose

In vivo evaluation of corneal biomechanics holds the potential for improving diagnosis and management of ocular diseases. We aimed to develop an ocular pulse elastography (OPE) technique to quantify corneal strains generated by naturally occurring pulsations of the intraocular pressure (IOP) using high-frequency ultrasound.

Methods

Simulated ocular pulses were induced in whole porcine and human donor globes to investigate the effects of physiologic variations in baseline IOP, ocular pulse amplitude, and frequency on corneal strains. Ocular pulse-induced strains were measured in additional globes before and after UVA-riboflavin-induced corneal crosslinking. The central cornea in each eye was imaged with a 50-MHz ultrasound imaging system and correlation-based speckle tracking of radiofrequency data was used to calculate tissue displacements and strains.

Results

Ocular pulse-induced corneal strains followed the cyclic changes of IOP. Both baseline IOP and ocular pulse amplitude had a significant influence on strain magnitude. Variations in pulse frequency within the normal human heart rate range did not introduce detectable changes in corneal strains. A significant decrease of corneal strain, as quantified by the OPE technique, was observed after corneal crosslinking. The extent of corneal stiffening (i.e., strain reduction) seemed to correlate with the initial strain magnitude.

Conclusions

This ex vivo study demonstrated the feasibility of the OPE method to quantify corneal strains generated by IOP pulsation and detect changes associated with corneal crosslinking treatment.

Translational Relevance

Integrating in vivo measurement of IOP and ocular pulse amplitude, the OPE method may lead to a new clinical tool for safe and quick biomechanical evaluations of the cornea.

Corneal Hydration Control during Ex Vivo Experimentation Using Poloxamers

Purpose: The purpose of this study was to develop an effective treatment method using poloxamers to restore and maintain physiological hydration in postmortem porcine and human corneas during ex vivo experimentation, and to compare corneal inflation response with or without treatment.Materials and Methods: Corneal buttons obtained from whole globes (n = 30 porcine, n = 8 human) were treated with various concentrations of poloxamer 188 (P188, a synthetic macromolecule surfactant) for 24 hrs to identify the concentration that would return the cornea to near-physiological hydration (i.e. H = 3.2). Whole globes (n = 12 porcine, n = 16 human) were also used to monitor central corneal thickness (CCT) during deswelling treatment. Inflation testing from 5 to 30 mmHg was performed in the porcine globes and a subset of human globes to characterize the mechanical response of the cornea after treatment.Results: Physiological hydration was obtained after 24 hrs immersion in 3.25% P188 for porcine corneas and 4.25% P188 treatment for human corneas. CCT was stabilized and returned to physiological levels after 24 hrs of treatment in 3.25% P188 in porcine (891 ± 66 µm) and 4.25% P188 in human (574 ± 34 µm) whole globes. Corneal axial strains at 30 mmHg were significantly larger at physiological hydration than in swollen cornea in both porcine (-6.42%±1.50% vs. -3.64%±1.05%, p = .004) and human (-2.85%±0.09% in vs. -1.53%±0.27%, p = .031) eyes.Conclusions: Our results suggest that P188 treatment was effective in restoring and maintaining near physiological corneal hydration during ex vivo testing, and hydration appeared to significantly impact corneal inflation response in both porcine and human eyes.

Mechanical Deformation of Human Optic Nerve Head and Peripapillary Tissue in Response to Acute IOP Elevation

Purpose

To measure the deformation of the human optic nerve head (ONH) and peripapillary tissue (PPT) in response to acute intraocular pressure (IOP) elevation.

Methods

The ONH and PPT of 14 human donor globes were imaged with high-frequency ultrasonography during inflation testing from 5 to 30 mm Hg. A correlation-based speckle tracking algorithm was used to compute tissue displacements, and the through-thickness, in-plane, and shear strains were calculated by using least-squares strain estimation methods. The ONH and PPT were segmented along the anterior-posterior direction and the nasal-temporal direction. Regional displacements and strains were analyzed and compared.

Results

The ONH displaced more posteriorly than the PPT in response to an acute IOP increase. Scleral canal expansion was minimal but correlated with ONH posterior displacement at all IOP levels. Through-thickness compression was concentrated in the anterior of both the ONH and the PPT. Shear was concentrated in the vicinity of the canal with higher shear in the peripheral ONH than the central ONH and higher shear in the PPT near the scleral canal than that further away from the canal.

Conclusions

High-resolution ultrasound speckle tracking showed a displacement mismatch between the ONH and the PPT, larger compressive strains in the direction of IOP loading in the anterior ONH and PPT, and higher shear strains in the periphery of ONH in response to acute IOP elevation in the human eye. These findings delineate the deformation patterns within and around the ONH and may help understand IOP-associated optic nerve damage.