Assessing the biomechanical properties of the porcine crystalline lens as a function of intraocular pressure with optical coherence elastography

Functional effects of the spatial-varying lens mechanical properties in accommodation

Lens biomechanical properties are critical for our eyes to accommodate. While it is well understood that lens mechanical properties change with age, different experimental techniques have been used over the years, with varying results on how the lens modulus changes. In this study, we developed a spatial-varying elasticity model to characterize the overall elastic modulus of the lens and establish its effect on accommodation. First, to validate the model, ex vivo porcine lenses underwent compression testing using biopsy punches of different diameters to change the percentage of nucleus within samples. Importantly, we found that, indeed, changing nucleus/cortex spatial ratio produces dramatic (∼7-fold) increase in overall sample modulus. Comparing the model with human lens spatial ratios, we demonstrate how changing spatial mechanics are more influential than peak modulus changes on overall elastic modulus. Next, in vivo clinical measurements of the spatial-varying lens modulus were used to generate a simplified mechanical-optical model of accommodation. We defined an ellipsoid lens with patient-derived modulus and geometry measurements, and a statics simulation and ray tracing analysis were performed through the deformed and undeformed lens. The resulting accommodation estimates agree with general accommodation expectations.

Storage-induced mechanical changes of porcine lenses assessed with optical coherence elastography and inverse finite element modeling

Introduction

In an effort of gaining a better understanding of the lens mechanics, ex vivo lenses samples are often used. Yet, ex vivo tissue might undergo important postmortem changes depending on the unavoidable preservation method employed. The purpose of this study was to assess how various storage conditions and the removal of the lens capsule affect the mechanical properties of ex vivo porcine lens samples.

Methods

A total of 81 freshly enucleated porcine eyes were obtained and divided into six groups and preserved differently. In the first three groups, the lens within the intact eye was preserved for 24 h by: (i) freezing at -80°C (n = 12), (ii) freezing at -20°C (n = 12), and (iii) refrigeration at +8°C (n = 12). In the remaining groups, the lenses were immediately extracted and treated as follows: (iv) kept intact, no storage (n = 12), (v) decapsulated, no storage (n = 21), and (vi) immersed in Minimum Essential Medium (MEM) at +8°C (n = 12) for 24 h. Frozen lenses were thawed at room temperature. Each lens was compressed between two glass lamella and subjected, first to a period of relaxation during which the compression force was recorded and second to an oscillating micro-compression while the deformation was recorded with a total of 256 subsequent B-scans via optical coherence tomography. The corresponding axial strain was retrieved via phase-sensitive image processing and subsequently used as input for an inverse finite element analysis (iFEA) to retrieve the visco-hyperelastic material properties of the lenses.

Results

After freezing at temperatures of -80°C and -20°C, the cortical strains increased by 14% (p = 0.01) and 34% (p < 0.001), and the nuclear strains decreased by 17% (p = 0.014) and 36% (p < 0.001), compared to the lenses tested immediately after postmortem, respectively. According to iFEA, this resulted from an increased ratio of the nuclear: cortical E-modulus (4.06 and 7.06) in -80°C and -20°C frozen lenses compared to fresh lenses (3.3). Decapsulation had the largest effect on the material constant C10, showing an increase both in the nucleus and cortex. Preservation of the intact eye in the refrigerator induced the least mechanical alterations in the lens, compared to the intact fresh condition.

Discussion

Combining iFEA with optical coherence elastography allowed us to identify important changes in the lens mechanics induced after different preserving ex vivo methods.

Optical coherence elastography and its applications for the biomechanical characterization of tissues

The biomechanical characterization of the tissues provides significant evidence for determining the pathological status and assessing the disease treatment. Incorporating elastography with optical coherence tomography (OCT), optical coherence elastography (OCE) can map the spatial elasticity distribution of biological tissue with high resolution. After the excitation with the external or inherent force, the tissue response of the deformation or vibration is detected by OCT imaging. The elastogram is assessed by stress-strain analysis, vibration amplitude measurements, and quantification of elastic wave velocities. OCE has been used for elasticity measurements in ophthalmology, endoscopy, and oncology, improving the precision of diagnosis and treatment of disease. In this article, we review the OCE methods for biomechanical characterization and summarize current OCE applications in biomedicine. The limitations and future development of OCE are also discussed during its translation to the clinic.

A Miniature Dual-Fiber Probe for Quantitative Optical Coherence Elastography

OBJECTIVE

Optical coherence elastography (OCE) allows for high resolution analysis of elastic tissue properties. However, due to the limited penetration of light into tissue, miniature probes are required to reach structures inside the body, e.g., vessel walls. Shear wave elastography relates shear wave velocities to quantitative estimates of elasticity. Generally, this is achieved by measuring the runtime of waves between two or multiple points. For miniature probes, optical fibers have been integrated and the runtime between the point of excitation and a single measurement point has been considered. This approach requires precise temporal synchronization and spatial calibration between excitation and imaging.

METHODS

We present a miniaturized dual-fiber OCE probe of 1 mm diameter allowing for robust shear wave elastography. Shear wave velocity is estimated between two optics and hence independent of wave propagation between excitation and imaging. We quantify the wave propagation by evaluating either a single or two measurement points. Particularly, we compare both approaches to ultrasound elastography.

RESULTS

Our experimental results demonstrate that quantification of local tissue elasticities is feasible. For homogeneous soft tissue phantoms, we obtain mean deviations of 0.15 ms-1 and 0.02 ms-1 for single-fiber and dual-fiber OCE, respectively. In inhomogeneous phantoms, we measure mean deviations of up to 0.54 ms-1 and 0.03 ms-1 for single-fiber and dual-fiber OCE, respectively.

CONCLUSION

We present a dual-fiber OCE approach that is much more robust in inhomogeneous tissues. Moreover, we demonstrate the feasibility of elasticity quantification in ex-vivo coronary arteries.

SIGNIFICANCE

This study introduces an approach for robust elasticity quantification from within the tissue.

Viscoelastic properties of porcine lenses using optical coherence elastography and inverse finite element analysis

The mechanical properties of the crystalline lens are crucial in determining the changes in lens shape that occur during the accommodation process and are also a major factor in the development of the two most prevalent age-related diseases of the lens, presbyopia and cataracts. However, a comprehensive understanding of these properties is currently lacking. Previous methods for characterizing the mechanical properties of the lens have been limited by the amount of data that could be collected during each test and the lack of complex material modeling. These limitations were mainly caused by the lack of imaging techniques that can provide data for the entire crystalline lens and the need for more complex models to describe the non-linear behavior of the lens. To address these issues, we characterized the mechanical properties of 13 porcine lenses during an ex vivo micro-controlled-displacement compression experiment using optical coherence elastography (OCE) and inverse finite element analysis (iFEA). OCE allowed us to quantify the internal strain distribution of the lens and differentiate between the different parts of the lens, while iFEA enabled us to implement an advanced material model to characterize the viscoelasticity of the lens nucleus and the relative stiffness gradient in the lens. Our findings revealed a pronounced and rapid viscoelastic behavior in the lens nucleus (g₁ = 0.39 ± 0.13, τ₁ = 5.01 ± 2.31 s) and identified the lens nucleus as the stiffest region, with a stiffness 4.42 ± 1.20 times greater than the anterior cortex and 3.47 ± 0.82 times greater than the posterior cortex. However, due to the complex nature of lens properties, it may be necessary to employ multiple tests simultaneously for a more comprehensive understanding of the crystalline lens.

The lens capsule significantly affects the viscoelastic properties of the lens as quantified by optical coherence elastography

The crystalline lens is a transparent, biconvex structure that has its curvature and refractive power modulated to focus light onto the retina. This intrinsic morphological adjustment of the lens to fulfill changing visual demands is achieved by the coordinated interaction between the lens and its suspension system, which includes the lens capsule. Thus, characterizing the influence of the lens capsule on the whole lens’s biomechanical properties is important for understanding the physiological process of accommodation and early diagnosis and treatment of lenticular diseases. In this study, we assessed the viscoelastic properties of the lens using phase-sensitive optical coherence elastography (PhS-OCE) coupled with acoustic radiation force (ARF) excitation. The elastic wave propagation induced by ARF excitation, which was focused on the surface of the lens, was tracked with phase-sensitive optical coherence tomography. Experiments were conducted on eight freshly excised porcine lenses before and after the capsular bag was dissected away. Results showed that the group velocity of the surface elastic wave, V, in the lens with the capsule intact (V=2.55±0.23 m/s) was significantly higher (p &lt; 0.001) than after the capsule was removed (V=1.19±0.25 m/s). Similarly, the viscoelastic assessment using a model that utilizes the dispersion of a surface wave showed that both Young’s modulus, E, and shear viscosity coefficient, η, of the encapsulated lens (E=8.14±1.10 kPa,η=0.89±0.093 Pa∙s) were significantly higher than that of the decapsulated lens (E=3.10±0.43 kPa,η=0.28±0.021 Pa∙s). These findings, together with the geometrical change upon removal of the capsule, indicate that the capsule plays a critical role in determining the viscoelastic properties of the crystalline lens.

Co-axial acoustic-based optical coherence vibrometry probe for the quantification of resonance frequency modes in ocular tissue

We present a co-axial acoustic-based optical coherence vibrometry probe (CoA-OCV) for vibro-acoustic resonance quantification in biological tissues. Sample vibrations were stimulated via a loudspeaker, and pre-compensation was used to calibrate the acoustic spectrum. Sample vibrations were measured via phase-sensitive swept-source optical coherence tomography (OCT). Resonance frequencies of corneal phantoms were measured at varying intraocular pressures (IOP), and dependencies on Young´s Modulus (E), phantom thickness and IOP were observed. Cycling IOP revealed hysteresis. For E = 0.3 MPa, resonance frequencies increased with IOP at a rate of 3.9, 3.7 and 3.5 Hz/mmHg for varied thicknesses and 1.7, 2.5 and 2.8 Hz/mmHg for E = 0.16 MPa. Resonance frequencies increased with thickness at a rate of 0.25 Hz/µm for E = 0.3 MPa, and 0.40 Hz/µm for E = 0.16 MPa. E showed the most predominant impact in the shift of the resonance frequencies. Full width at half maximum (FWHM) of the resonance modes increased with increasing thickness and decreased with increasing E. Only thickness and E contributed to the variance of FWHM. In rabbit corneas, resonance frequencies of 360-460 Hz were observed. The results of the current study demonstrate the feasibility of CoA-OCV for use in future OCT-V studies.

High-Frequency Ultrasound Elastography to Assess the Nonlinear Elastic Properties of the Cornea and Ciliary Body

Mechanical properties of the anterior anatomical structures of the eye, such as the cornea and ciliary body, play a key role in the ocular function and homeostasis. However, measuring the biomechanical properties of the anterior ocular structures, especially deeper structures, such as the ciliary body, remains a challenge due to the lack of high-resolution imaging tools. Herein, we implement a mechanical shaker-based high-frequency ultrasound elastography technique that can track the induced elastic wave propagation to assess the linear and nonlinear elastic properties of anterior ocular structures. The findings of this study advance our understanding of the role of anterior ocular structures in the pathogenesis of different ocular disorders, such as glaucoma.

Elasticity measurements of ocular anterior and posterior segments using optical coherence elastography

The changes of biomechanical properties, especially the elasticity of the ocular tissues, are closely related to some ophthalmic diseases. Currently, the ophthalmic optical coherence elastography (OCE) systems are dedicated either to the anterior segment or to the retina. The elasticity measurements of the whole eye remain challenging. Here we demonstrated an acoustic radiation force optical coherence elastography (ARF-OCE) method to quantify the elasticity of the cornea and the retina. The experiment results show that the Young's moduli of the cornea and the retina were 16.66 ± 6.51 kPa and 207.96 ± 4.75 kPa, respectively. Our method can measure the elasticity of the anterior segment and the posterior segment, and provides a powerful tool to enhance ophthalmology research.

Wave-based optical coherence elastography: The 10-year perspective

After 10 years of progress and innovation, optical coherence elastography (OCE) based on the propagation of mechanical waves has become one of the major and the most studied OCE branches, producing a fundamental impact in the quantitative and nondestructive biomechanical characterization of tissues. Preceding previous progress made in ultrasound and magnetic resonance elastography; wave-based OCE has pushed to the limit the advance of three major pillars: (1) implementation of novel wave excitation methods in tissues, (2) understanding new types of mechanical waves in complex boundary conditions by proposing advance analytical and numerical models, and (3) the development of novel estimators capable of retrieving quantitative 2D/3D biomechanical information of tissues. This remarkable progress promoted a major advance in answering basic science questions and the improvement of medical disease diagnosis and treatment monitoring in several types of tissues leading, ultimately, to the first attempts of clinical trials and translational research aiming to have wave-based OCE working in clinical environments. This paper summarizes the fundamental up-to-date principles and categories of wave-based OCE, revises the timeline and the state-of-the-art techniques and applications lying in those categories, and concludes with a discussion on the current challenges and future directions, including clinical translation research.

Two-dimensional (2D) dynamic vibration optical coherence elastography (DV-OCE) for evaluating mechanical properties: a potential application in tissue engineering

Mechanical properties in tissues are an important indicator because they are associated with disease states. One of the well-known excitation sources in optical coherence elastography (OCE) to determine mechanical properties is acoustic radiation force (ARF); however, a complicated focusing alignment cannot be avoided. Another excitation source is a piezoelectric (PZT) stack to obtain strain images via compression, which can affect the intrinsic mechanical properties of tissues in tissue engineering. In this study, we report a new technique called two-dimensional (2D) dynamic vibration OCE (DV-OCE) to evaluate 2D wave velocities without tedious focusing alignment procedures and is a non-contact method with respect to the samples. The three-dimensional (3D) Fourier transform was utilized to transfer the traveling waves (x, y, t) into 3D k-space (k_x, k_y, f). A spatial 2D wavenumber filter and multi-angle directional filter were employed to decompose the waves with omni-directional components into four individual traveling directions. The 2D local wave velocity algorithm was used to calculate a 2D wave velocity map. Six materials, two homogeneous phantoms with 10 mm thickness, two homogeneous phantoms with 2 mm thickness, one heterogeneous phantom with 2 mm diameter inclusion and an ex vivo porcine kidney, were examined in this study. In addition, the ARF-OCE was used to evaluate wave velocities for comparison. Numerical simulations were performed to validate the proposed 2D dynamic vibration OCE technique. We demonstrate that the experimental results were in a good agreement with the results from ARF-OCE (transient OCE) and numerical simulations. Our proposed 2D dynamic vibration OCE could potentially pave the way for mechanical evaluation in tissue engineering and for laboratory translation with easy-to-setup and contactless advantages.

Characterization of natural frequencies from nanoscale tissue oscillations using dynamic optical coherence elastography

We demonstrate the use of OCT-based elastography for soft-tissue characterization using natural frequency oscillations. Sub-micrometer to sub-nanometer oscillations were induced in tissue phantoms and human cornea in vivo by perpendicular air-pulse stimulation and observed by common-path OCT imaging (sensitivity: 0.24 nm). Natural frequency and damping ratio were acquired in temporal and frequency domains using a single degree of freedom method. The dominant natural frequency was constant for different stimulation pressures (4-32 Pa) and measured distances (0.3-5.3 mm), and decreased as the sample thickness increased. The dominant natural frequencies of 0.75-2% agar phantoms were 127-774 Hz (mean coefficient of variation [CV]: 0.9%), and correlated with the square root of Young's moduli (16.5-117.8 kPa, mean CV: 5.8%). These preliminary studies show repeatable in vivo corneal natural frequency measurements (259 Hz, CV: 1.9%). This novel OCE approach can distinguish tissues and materials with different mechanical properties using the small-amplitude tissue oscillation features, and is suitable for characterizing delicate tissues in vivo such as the eye.

Multimodal quantitative optical elastography of the crystalline lens with optical coherence elastography and Brillouin microscopy

Assessing the biomechanical properties of the crystalline lens can provide crucial information for diagnosing disease and guiding precision therapeutic interventions. Existing noninvasive methods have been limited to global measurements. Here, we demonstrate the quantitative assessment of the elasticity of crystalline lens with a multimodal optical elastography technique, which combines dynamic wave-based optical coherence elastography (OCE) and Brillouin microscopy to overcome the drawbacks of individual modalities. OCE can provide direct measurements of tissue elasticity rapidly and quantitatively, but it is a challenge to image transparent samples such as the lens because this technique relies on backscattered light. On the other hand, Brillouin microscopy can map the longitudinal modulus with micro-scale resolution in transparent samples. However, the relationship between Brillouin-deduced modulus and Young's modulus is not straightforward and sample dependent. By combining these two techniques, we can calibrate Brillouin measurements with OCE, based on the same sample, allowing us to completely map the Young's modulus of the crystalline lens. The combined system was first validated with tissue-mimicking gelatin phantoms of varying elasticities (N = 9). The OCE data was used to calibrate the Brillouin shift measurements and subsequently map the Young's modulus of the phantoms. After validation, OCE and Brillouin measurements were performed on ex-vivo porcine lenses (N = 6), and the Young's modulus of the lenses was spatially mapped. The results show a strong correlation between Young's moduli measured by OCE and longitudinal moduli measured by Brillouin microscopy. The correlation coefficient R was 0.98 for the phantoms and 0.94 for the lenses, respectively. The mean Young's modulus of the anterior and posterior lens was 1.98 ± 0.74 kPa and 2.93 ± 1.13 kPa, respectively, and the Young's modulus of the lens nucleus was 11.90 ± 2.94 kPa. The results presented in this manuscript open a new way for truly quantitative biomechanical mapping of optically transparent (or low scattering) tissues in 3D.

Intraocular Pressure Monitoring Following Islet Transplantation to the Anterior Chamber of the Eye

Intraocular islet transplantation was investigated as a new procedure to treat diabetes. The development of this procedure requires close monitoring of the function of both eye and islet graft. We developed a soft, smart contact lens to monitor the intraocular pressure and applied this for noninvasive monitoring in association with the intraocular islet transplantation in diabetes. A strain sensor inside the lens can detect detailed changes in intraocular pressure by focusing the strain only in the desired, selective area of the contact lens. In addition, this smart contact lens can transmit the real-time value of the intraocular pressure wirelessly using an antenna. The wireless measurement of intraocular pressure that was obtained using this contact lens had a high correlation with the intraocular pressure measured by a rebound tonometer, thereby proving the good accuracy of the contact lens sensor. In the initial period, a slight elevation of intraocular pressure was observed, but the pressure returned to normal in the initial period after the transplantation. This type of monitoring will provide important information on potential changes in the intraocular pressure associated with the transplantation procedure, and it enables appropriate clinical safety steps to be taken, if needed.

Advances in Doppler optical coherence tomography and angiography

Since the first demonstration of Doppler optical coherence tomography (OCT) in 1997, several functional extensions of Doppler OCT have been developed, including velocimetry, angiogram, and optical coherence elastography. These functional techniques have been widely used in research and clinical applications, particularly in ophthalmology. Here, we review the principles, representative methods, and applications of different Doppler OCT techniques, followed by discussion on the innovations, limitations, and future directions of each of these techniques.

Simultaneously imaging and quantifying in vivo mechanical properties of crystalline lens and cornea using optical coherence elastography with acoustic radiation force excitation

The crystalline lens and cornea comprise the eye's optical system for focusing light in human vision. The changes in biomechanical properties of the lens and cornea are closely associated with common diseases, including presbyopia and cataract. Currently, most in vivo elasticity studies of the anterior eye focus on the measurement of the cornea, while lens measurement remains challenging. To better understand the anterior segment of the eye, we developed an optical coherence elastography system utilizing acoustic radiation force excitation to simultaneously assess the elasticities of the crystalline lens and the cornea in vivo. A swept light source was integrated into the system to provide an enhanced imaging range that covers both the lens and the cornea. Additionally, the oblique imaging approach combined with orthogonal excitation also improved the image quality. The system was tested through first ex vivo and then in vivo experiments using a rabbit model. The elasticities of corneal and lens tissue in an excised normal whole-globe and a cold cataract model were measured to reveal that cataractous lenses have a higher Young's modulus. Simultaneous in vivo elasticity measurements of the lens and cornea were performed in a rabbit model to demonstrate the correlations between elasticity and intraocular pressure and between elasticity and age. To the best of our knowledge, we demonstrated the first in vivo elasticity of imaging of both the lens and cornea using acoustic radiation force-optical coherence elastography, thereby providing a potential powerful clinical tool to advance ophthalmic research in disorders affecting the lens and the cornea.

Optical coherence elastography of cold cataract in porcine lens

Cataract is one of the most prevalent causes of blindness around the world. Understanding the mechanisms of cataract development and progression is important for clinical diagnosis and treatment. Cold cataract has proven to be a robust model for cataract formation that can be easily controlled in the laboratory. There is evidence that the biomechanical properties of the lens can be significantly changed by cataract. Therefore, early detection of cataract, as well as evaluation of therapies, could be guided by characterization of lenticular biomechanical properties. In this work, we utilized optical coherence elastography (OCE) to monitor the changes in biomechanical properties of ex vivo porcine lenses during formation of cold cataract. Elastic waves were induced in the porcine lenses by a focused micro air-pulse while the lenses were cooled, and the elastic wave velocity was translated to Young's modulus of the lens. The results show an increase in the stiffness of the lens due to formation of the cold cataract (from 11.3  ±  3.4 to 21.8  ±  7.8  kPa). These results show a relation between lens opacity and stiffness and demonstrate that OCE can assess lenticular biomechanical properties and may be useful for detecting and potentially characterizing cataracts.