Reverberant 3D optical coherence elastography maps the elasticity of individual corneal layers

Mouse brain elastography changes with sleep/wake cycles, aging, and Alzheimer's disease

Understanding the physiological processes in aging and how neurodegenerative disorders affect cognitive function is a high priority for advancing human health. One specific area of recently enabled research is the in vivo biomechanical state of the brain. This study utilized reverberant optical coherence elastography, a high-resolution elasticity imaging method, to investigate stiffness changes during the sleep/wake cycle, aging, and Alzheimer's disease in murine models. Four-dimensional scans of 44 wildtype mice, 13 mice with deletion of aquaporin-4 water channel, and 12 mice with Alzheimer-related pathology (APP/PS1) demonstrated that (1) cortical tissue became softer (on the order of a 10% decrease in shear wave speed) when young wildtype mice transitioned from wake to anesthetized, yet this effect was lost in aging and with mice overexpressing amyloid-β or lacking the water channel AQP4. (2) Cortical stiffness increased with age in all mice lines, but wildtype mice exhibited the most prominent changes as a function of aging. The study provides novel insight into the brain's biomechanics, the constraints of fluid flow, and how the state of brain activity affects basic properties of cortical tissues.

Spatial mapping of corneal biomechanical properties using wave-based optical coherence elastography

Quantifying the mechanical properties of the cornea can provide valuable insights into the occurrence and progression of keratoconus, as well as the effectiveness of corneal crosslinking surgery. This study presents a non-contact and non-invasive wave-based optical coherence elastography system that utilizes air-pulse stimulation to create a two-dimensional map of corneal elasticity. Homogeneous and dual concentration phantoms were measured with the sampling of 25 × 25 points over a 6.6 × 6.6 mm2 area, to verify the measurement capability for elastic mapping and the spatial resolution (0.91 mm). The velocity of elastic waves distribution of porcine corneas before and after corneal crosslinking surgery were further mapped, showing a significant change in biomechanics in crosslinked region. This system features non-invasiveness and high resolution, holding great potential for application in ophthalmic clinics.

Dolutegravir induces FOLR1 expression during brain organoid development

During the first month of pregnancy, the brain and spinal cord are formed through a process called neurulation. However, this process can be altered by low serum levels of folic acid, environmental factors, or genetic predispositions. In 2018, a surveillance study in Botswana, a country with a high incidence of human immunodeficiency virus (HIV) and lacking mandatory food folate fortification programs, found that newborns whose mothers were taking dolutegravir (DTG) during the first trimester of pregnancy had an increased risk of neural tube defects (NTDs). As a result, the World Health Organization and the U.S. Food and Drug Administration have issued guidelines emphasizing the potential risks associated with the use of DTG-based antiretroviral therapies during pregnancy. To elucidate the potential mechanisms underlying the DTG-induced NTDs, we sought to assess the potential neurotoxicity of DTG in stem cell-derived brain organoids. The gene expression of brain organoids developed in the presence of DTG was analyzed by RNA sequencing, Optical Coherence Tomography (OCT), Optical Coherence Elastography (OCE), and Brillouin microscopy. The sequencing data shows that DTG induces the expression of the folate receptor (FOLR1) and modifies the expression of genes required for neurogenesis. The Brillouin frequency shift observed at the surface of DTG-exposed brain organoids indicates an increase in superficial tissue stiffness. In contrast, reverberant OCE measurements indicate decreased organoid volumes and internal stiffness.

Noncontact elasticity measurement of hydrogels in a culture dish using reverberant optical coherence elastography

Estimating the elasticity of hydrogel phantoms in a cell culture plane is important for understanding the cell behavior in response to various types of mechanical stimuli. Hence, a noncontact tool for measuring the elastic properties of hydrogel phantoms in such three-dimensional cell cultures is required. A well-known method to determine the mechanical properties of hydrogels is the transient wave method. However, due to the multiple reflections of waves from the boundaries, a bigger cell culture plane or multiple directional filters may be required. In this study, we utilized reverberant shear wave elastography, which is based on the autocorrelation principle, to evaluate the shear wave speed in hydrogel samples within a culture dish. Numerical simulations were performed first to confirm the validity of the reverberant elastography method. Subsequently, we used this method to measure the wave speeds in hydrogel phantoms with different concentrations. Shear rheology tests were also performed, and their results were found to be in good agreement with the measured shear wave speeds. The proposed method could be useful for measuring the elasticity of tissues in tissue engineering applications in an inexpensive and noncontact manner.

Full-field noise-correlation elastography for in-plane mechanical anisotropy imaging

Elastography contrast imaging has great potential for the detection and characterization of abnormalities in soft biological tissues to help physicians in diagnosis. Transient shear-waves elastography has notably shown promising results for a range of clinical applications. In biological soft tissues such as muscle, high mechanical anisotropy implies different stiffness estimations depending on the direction of the measurement. In this study, we propose the evolution of a noise-correlation elastography approach for in-plane anisotropy mapping. This method is shown to retrieve anisotropy from simulation images before being validated on agarose anisotropic tissue-mimicking phantoms, and the first results on in-vivo biological fibrous tissues are presented.

Dual-convolutional neural network-enhanced strain estimation method for optical coherence elastography

Strain estimation is vital in phase-sensitive optical coherence elastography (PhS-OCE). In this Letter, we introduce a novel, to the best of our knowledge, method to improve strain estimation by using a dual-convolutional neural network (Dual-CNN). This approach requires two sets of PhS-OCE systems: a high-resolution system for high-quality training data and a cost-effective standard-resolution system for practical measurements. During training, high-resolution strain results acquired from the former system and the pre-existing strain estimation CNN serve as label data, while the narrowed light source-acquired standard-resolution phase results act as input data. By training a new network with this data, high-quality strain results can be estimated from standard-resolution PhS-OCE phase results. Comparison experiments show that the proposed Dual-CNN can preserve the strain quality even when the light source bandwidth is reduced by over 80%.

In Vivo Evaluation of Corneal Biomechanics Following Cross-Linking Surgeries Using Optical Coherence Elastography in a Rabbit Model of Keratoconus

Purpose

Validation of the feasibility of novel acoustic radiation force optical coherence elastography (ARF-OCE) for the evaluation of biomechanical enhancement of the in vivo model of keratoconus by clinical cross-linking (CXL) surgery.

Methods

Twelve in vivo rabbit corneas were randomly divided into two groups. Both groups were treated with collagenase type II, and a keratoconus model was obtained. Then, the two groups were treated with CXL procedures with different irradiation energy of 15 J and 30 J (CXL-15 J and CXL-30 J, respectively). An ARF-OCE probe with an ultrasmall ultrasound transducer was used to detect the biomechanical properties of cornea. An antisymmetric Lamb wave model was combined with the frequency dispersion relationship to achieve depth-resolved elastography.

Results

Compared with the phase velocity of the Lamb wave in healthy corneas (approximately 3.96 ± 0.27 m/s), the phase velocity of the Lamb wave was lower in the keratoconus region (P < 0.05), with an average value of 3.12 ± 0.12 m/s. Moreover, the corneal stiffness increased after CXL treatment (P < 0.05), and the average phase velocity of the Lamb wave was 4.3 ± 0.19 m/s and 4.54 ± 0.13 m/s after CXL-15 J and CXL-30 J treatment.

Conclusions

The Young's moduli of the keratoconus regions were significantly lower than the healthy corneas. Moreover, the Young's modulus of the keratoconus regions was significantly higher after CXL-30 J treatment than after CXL-15 J treatment. We demonstrated that the ARF-OCE technique has great potential in screening keratoconus and guiding clinical CXL treatment.

Translational Relevance

This work accelerates the clinical translation of OCE systems using ultrasmall ultrasound transducers and is used to guide CXL procedures.

Photoacoustic micro-viscoelastography for mapping mechanocellular properties

Cellular biomechanical properties provide essential insights into biological functions regarding health and disease. Current measurements of the biomechanical properties of cells require physical contact with cells or pre-loading on the cells. Here, we have developed photoacoustic micro-viscoelastography (PAMVE), which utilizes the phase characteristics of photoacoustic (PA) response, for mapping mechanocellular properties in a load-free manner. PAMVE realizes the local viscoelasticity measurement on the macrophages and red blood cells with micrometer scale. Furthermore, PAMVE can successfully identify the adipose cell and skeletal muscle cell due to the difference in their composition-related biomechanical properties. PAMVE represents an irreplaceable option for interrogating characteristic mechanocellular properties, opening the possibility of studying cellular mechanobiology and pathophysiology.

Nonlinear fourth-order elastic characterization of the cornea using torsional wave elastography

Measuring the mechanical nonlinear properties of the cornea remains challenging due to the lack of consensus in the methodology and in the models that effectively predict its behaviour. This study proposed developing a procedure to reconstruct nonlinear fourth-order elastic properties of the cornea based on a mathematical model derived from the theory of Hamilton et al. and using the torsional wave elastography (TWE) technique. In order to validate its diagnostic capability of simulated pathological conditions, two different groups were studied, non-treated cornea samples (n=7), and ammonium hydroxide ([Formula: see text]) treated samples (n=7). All the samples were measured in-plane by a torsional wave device by increasing IOP from 5 to 25 mmHg with 5 mmHg steps. The results show a nonlinear variation of the shear wave speed with the IOP, with higher values for higher IOPs. Moreover, the shear wave speed values of the control group were higher than those of the treated group. The study also revealed significant differences between the control and treated groups for the Lamé parameter [Formula: see text] (25.9-6.52 kPa), third-order elastic constant A (215.09-44.85 kPa), and fourth-order elastic constant D (523.5-129.63 kPa), with p-values of 0.010, 0.024, and 0.032, respectively. These findings demonstrate that the proposed procedure can distinguish between healthy and damaged corneas, making it a promising technique for detecting diseases associated with IOP alteration, such as corneal burns, glaucoma, or ocular hypertension.

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.

Possible depth-resolved reconstruction of shear moduli in the cornea following collagen crosslinking (CXL) with optical coherence tomography and elastography

Corneal collagen crosslinking (CXL) is commonly used to prevent or treat keratoconus. Although changes in corneal stiffness induced by CXL surgery can be monitored with non-contact dynamic optical coherence elastography (OCE) by tracking mechanical wave propagation, depth dependent changes are still unclear if the cornea is not crosslinked through the whole depth. Here, phase-decorrelation measurements on optical coherence tomography (OCT) structural images are combined with acoustic micro-tapping (AµT) OCE to explore possible reconstruction of depth-dependent stiffness within crosslinked corneas in an ex vivo human cornea sample. Experimental OCT images are analyzed to define the penetration depth of CXL into the cornea. In a representative ex vivo human cornea sample, crosslinking depth varied from ∼100 µm in the periphery to ∼150 µm in the cornea center and exhibited a sharp in-depth transition between crosslinked and untreated areas. This information was used in an analytical two-layer guided wave propagation model to quantify the stiffness of the treated layer. We also discuss how the elastic moduli of partially CXL-treated cornea layers reflect the effective engineering stiffness of the entire cornea to properly quantify corneal deformation.

Multifocal acoustic radiation force-based reverberant optical coherence elastography for evaluation of ocular globe biomechanical properties

Significance

Quantifying the biomechanical properties of the whole eye globe can provide a comprehensive understanding of the interactions among interconnected ocular components during dynamic physiological processes. By doing so, clinicians and researchers can gain valuable insights into the mechanisms underlying ocular diseases, such as glaucoma, and design interventions tailored to each patient's unique needs.

Aim

The aim of this study was to evaluate the feasibility and effectiveness of a multifocal acoustic radiation force (ARF) based reverberant optical coherence elastography (RevOCE) technique for quantifying shear wave speeds in different ocular components simultaneously.

Approach

We implemented a multifocal ARF technique to generate reverberant shear wave fields, which were then detected using phase-sensitive optical coherence tomography. A 3D-printed acoustic lens array was employed to manipulate a collimated ARF beam generated by an ultrasound transducer, producing multiple focused ARF beams on mouse eye globes ex vivo. RevOCE measurements were conducted using an excitation pulse train consisting of 10 cycles at 3 kHz, followed by data processing to produce a volumetric map of the shear wave speed.

Results

The results show that the system can successfully generate reverberant shear wave fields in the eye globe, allowing for simultaneous estimation of shear wave speeds in various ocular components, including cornea, iris, lens, sclera, and retina. A comparative analysis revealed notable differences in wave speeds between different parts of the eye, for example, between the apical region of the cornea and the pupillary zone of the iris (p = 0.003 ). Moreover, the study also revealed regional variations in the biomechanical properties of ocular components as evidenced by greater wave speeds near the apex of the cornea compared to its periphery.

Conclusions

The study demonstrated the effectiveness of RevOCE based on a non-invasive multifocal ARF for assessing the biomechanical properties of the whole eyeball. The findings indicate the potential to provide a comprehensive understanding of the mechanical behavior of the whole eye, which could lead to improved diagnosis and treatment of ocular diseases.

Possible depth-resolved reconstruction of shear moduli in the cornea following collagen crosslinking (CXL) with optical coherence tomography and elastography

Corneal collagen crosslinking (CXL) is commonly used to prevent or treat keratoconus. Although changes in corneal stiffness induced by CXL surgery can be monitored with non-contact dynamic optical coherence elastography (OCE) by tracking mechanical wave propagation, depth dependent changes are still unclear if the cornea is not crosslinked through the whole depth. Here, phase-decorrelation measurements on optical coherence tomography (OCT) structural images are combined with acoustic micro-tapping (A$\mu$T) OCE to explore possible reconstruction of depth-dependent stiffness within crosslinked corneas in an ex vivo human cornea sample. Experimental OCT images are analyzed to define the penetration depth of CXL into the cornea. In a representative ex vivo human cornea sample, crosslinking depth varied from $\sim 100\mu m$ in the periphery to $\sim 150\mu m$ in the cornea center and exhibited a sharp in-depth transition between crosslinked and untreated areas. This information was used in an analytical two-layer guided wave propagation model to quantify the stiffness of the treated layer. We also discuss how the elastic moduli of partially CXL-treated cornea layers reflect the effective engineering stiffness of the entire cornea to properly quantify corneal deformation.

Reverberant optical coherence elastography using multifocal acoustic radiation force

In this study, we introduce a multifocal acoustic radiation force source that combines an ultrasound transducer and a 3D-printed acoustic lens for application in reverberant optical coherence elastography (Rev-OCE). An array of plano-concave acoustic lenses, each with an 11.8 mm aperture diameter, were used to spatially distribute the acoustic energy generated by a 1 MHz planar ultrasound transducer, producing multiple focal spots on a target plane. These focal spots generate reverberant shear wave fields detected by the optical coherence tomography (OCT) system. The effectiveness of the multifocal Rev-OCE system in probing mechanical properties with high resolution is demonstrated in layered gelatin phantoms.

Debiased ambient vibrations optical coherence elastography to profile cell, organoid and tissue mechanical properties

The role of the mechanical environment in defining tissue function, development and growth has been shown to be fundamental. Assessment of the changes in stiffness of tissue matrices at multiple scales has relied mostly on invasive and often specialist equipment such as AFM or mechanical testing devices poorly suited to the cell culture workflow.In this paper, we have developed a unbiased passive optical coherence elastography method, exploiting ambient vibrations in the sample that enables real-time noninvasive quantitative profiling of cells and tissues. We demonstrate a robust method that decouples optical scattering and mechanical properties by actively compensating for scattering associated noise bias and reducing variance. The efficiency for the method to retrieve ground truth is validated in silico and in vitro, and exemplified for key applications such as time course mechanical profiling of bone and cartilage spheroids, tissue engineering cancer models, tissue repair models and single cell. Our method is readily implementable with any commercial optical coherence tomography system without any hardware modifications, and thus offers a breakthrough in on-line tissue mechanical assessment of spatial mechanical properties for organoids, soft tissues and tissue engineering.

Novel acoustic radiation force optical coherence elastography based on ultrasmall ultrasound transducer for biomechanics evaluation of in vivo cornea

We developed a novel acoustic radiation force optical coherence elastography (ARF-OCE) based on an ultrasmall ultrasound transducer for quantitative biomechanics evaluations of in vivo cornea. A custom single-sided meta-ultrasonic transducer with an outer diameter of 1.8 mm, focal spot diameter of 1.6 mm, central frequency of 930 kHz, and focal length of 0.8 mm was applied to excite the sample. The sample arm of the ARF-OCE system employed a three-dimensional printed holder that allowed for ultrasound excitation and ARF-OCE detection. The phase-resolved algorithm was combined with a Lamb wave model to depth-resolved evaluate corneal biomechanics after keratoconus and cross-linking treatments (CXL). The results showed that, compare to the healthy cornea, the Lamb wave velocity was significantly reduced in the keratoconus, increased in the cornea after CXL, and increased with cross-linked irradiation energy in the cornea. These results indicated the good clinical translation potential of the proposed novel ARF-OCE.

Fluid compartments influence elastography of the aging mouse brain

Objective. Elastography of the brain has the potential to reveal subtle but clinically important changes in the structure and composition as a function of age, disease, and injury.Approach. In order to quantify the specific effects of aging on mouse brain elastography, and to determine the key factors influencing observed changes, we applied optical coherence tomography reverberant shear wave elastography at 2000 Hz to a group of wild-type healthy mice ranging from young to old age.Main results. We found a strong trend towards increasing stiffness with age, with an approximately 30% increase in shear wave speed from 2 months to 30 months within this sampled group. Furthermore, this appears to be strongly correlated with decreasing measures of whole brain fluid content, so older brains have less water and are stiffer. Rheological models are applied, and the strong effect is captured by specific assignment of changes to the glymphatic compartment of the brain fluid structures along with a correlated change in the parenchymal stiffness.Significance. Short-term and longer-term changes in elastography measures may provide a sensitive biomarker of progressive and fine-scale changes in the glymphatic fluid channels and parenchymal components of the brain.

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.

Reverberant magnetic resonance elastographic imaging using a single mechanical driver

Reverberant elastography provides fast and robust estimates of shear modulus; however, its reliance on multiple mechanical drivers hampers clinical utility. In this work, we hypothesize that for constrained organs such as the brain, reverberant elastography can produce accurate magnetic resonance elastograms with a single mechanical driver. To corroborate this hypothesis, we performed studies on healthy volunteers (n= 3); and a constrained calibrated brain phantom containing spherical inclusions with diameters ranging from 4-18 mm. In both studies (i.e. phantom and clinical), imaging was performed at frequencies of 50 and 70 Hz. We used the accuracy and contrast-to-noise ratio performance metrics to evaluate reverberant elastograms relative to those computed using the established subzone inversion method. Errors incurred in reverberant elastograms varied from 1.3% to 16.6% when imaging at 50 Hz and 3.1% and 16.8% when imaging at 70 Hz. In contrast, errors incurred in subzone elastograms ranged from 1.9% to 13% at 50 Hz and 3.6% to 14.9% at 70 Hz. The contrast-to-noise ratio of reverberant elastograms ranged from 63.1 to 73 dB compared to 65 to 66.2 dB for subzone elastograms. The average global brain shear modulus estimated from reverberant and subzone elastograms was 2.36 ± 0.07 kPa and 2.38 ± 0.11 kPa, respectively, when imaging at 50 Hz and 2.70 ± 0.20 kPa and 2.89 ± 0.60 kPa respectively, when imaging at 70 Hz. The results of this investigation demonstrate that reverberant elastography can produce accurate, high-quality elastograms of the brain with a single mechanical driver.

Recent advances in optical elastography and emerging opportunities in the basic sciences and translational medicine [Invited]

Optical elastography offers a rich body of imaging capabilities that can serve as a bridge between organ-level medical elastography and single-molecule biophysics. We review the methodologies and recent developments in optical coherence elastography, Brillouin microscopy, optical microrheology, and photoacoustic elastography. With an outlook toward maximizing the basic science and translational clinical impact of optical elastography technologies, we discuss potential ways that these techniques can integrate not only with each other, but also with supporting technologies and capabilities in other biomedical fields. By embracing cross-modality and cross-disciplinary interactions with these parallel fields, optical elastography can greatly increase its potential to drive new discoveries in the biomedical sciences as well as the development of novel biomechanics-based clinical diagnostics and therapeutics.

In Vivo Evaluation of the Effects of SMILE with Different Amounts of Stromal Ablation on Corneal Biomechanics by Optical Coherence Elastography

This work aims to depth-resolved quantitatively analyze the effect of different stromal ablation amounts on the corneal biomechanical properties during small incision lenticule extraction (SMILE) using optical coherence elastography (OCE). A 4.5-MHz ultrasonic transducer was used to excite elastic waves in the corneal tissue. The OCE system combined with the antisymmetric Lamb wave model was employed to achieve a high-resolution, high-sensitivity, and depth-resolved quantitative detection of the corneal Young's modulus. Eighteen rabbits were randomly divided into three groups; each group had six rabbits. The first and second groups underwent -3D and -6D SMILE surgeries, and the third group was the control group, respectively. Young's modulus of the corneal cap and residual stromal bed (RSB) were both increased after SMILE, which shared the stress under intraocular pressure (IOP). Furthermore, the Young's modulus of both the corneal cap and RSB after 3D SMILE group were significantly lower than that in the -6D group, which indicated that the increases in the post-operative corneal Young's modulus were positively correlated with the amount of stromal ablation. The OCE system for quantitative spatial characterization of corneal biomechanical properties can provide useful information on the extent of safe ablation for SMILE procedures.

Theory of sleep/wake cycles affecting brain elastography

As elastography of the brain finds increasing clinical applications, fundamental questions remain about baseline viscoelastic properties of the brainin vivo. Furthermore, the underlying mechanisms of how and why elastographic measures can change over time are still not well understood. To study these issues, reverberant shear wave elastography using an optical coherence tomography scanner is implemented on a mouse model, both under awake conditions and in a sleep state where there are known changes in the glymphatic fluid flow system in the brain. We find that shear wave speed, a measure of stiffness, changes by approximately 12% between the two states, sleep versus awake, in the entire cortical brain imaging volume. Our microchannel flow model of biphasic (fluid plus solid) tissue provides a plausible rheological model based on the fractal branching vascular and perivascular system, plus a second parallel system representing the finer scale glymphatic fluid microchannels. By adjusting the glymphatic system fluid volume proportional to the known sleep/wake changes, we are able to approximately predict the measured shear wave speeds and their change with the state of the glymphatic system. The advantages of this model are that its main parameters are derived from anatomical measures and are linked to other major derivations of branching fluid structures including Murray's Law. The implications for clinical studies are that elastography of the brain is strongly influenced by the regulation or dysregulation of the vascular, perivascular, and glymphatic systems.

In vivo stiffness measurement of epidermis, dermis, and hypodermis using broadband Rayleigh-wave optical coherence elastography

Traveling-wave optical coherence elastography (OCE) is a promising technique to measure the stiffness of biological tissues. While OCE has been applied to relatively homogeneous samples, tissues with significantly varying elasticity through depth pose a challenge, requiring depth-resolved measurement with sufficient resolution and accuracy. Here, we develop a broadband Rayleigh-wave OCE technique capable of measuring the elastic moduli of the 3 major skin layers (epidermis, dermis, and hypodermis) reliably by analyzing the dispersion of leaky Rayleigh surface waves over a wide frequency range of 0.1-10 kHz. We show that a previously unexplored, high frequency range of 4-10 kHz is critical to resolve the thin epidermis, while a low frequency range of 0.2-1 kHz is adequate to probe the dermis and deeper hypodermis. We develop a dual bilayer-based inverse model to determine the elastic moduli in all 3 layers and verify its high accuracy with finite element analysis and skin-mimicking phantoms. Finally, the technique is applied to measure the forearm skin of healthy volunteers. The Young's modulus of the epidermis (including the stratum corneum) is measured to be ∼ 4 MPa at 4-10 kHz, whereas Young's moduli of the dermis and hypodermis are about 40 and 15 kPa, respectively, at 0.2-1 kHz. Besides dermatologic applications, this method may be useful for the mechanical analysis of various other layered tissues with sub-mm depth resolution. STATEMENT OF SIGNIFICANCE: To our knowledge, this is the first study that resolves the stiffness of the thin epidermis from the dermis and hypodermis, made possible by using high-frequency (4 - 10 kHz) elastic waves and optical coherence elastography. Beyond the skin, this technique may be useful for mechanical characterizations of various layered biomaterials and tissues.

Comprehensive experimental assessments of rheological models' performance in elastography of soft tissues

Elastography researchers have utilized several rheological models to characterize soft tissue viscoelasticity over the past thirty years. Due to the frequency-dependent behavior of viscoelastic parameters as well as the different techniques and frequencies employed in various studies of soft tissues, rheological models have value in standardizing disparate techniques via explicit mathematical representations. However, the important question remains: which of the several available models should be considered for widespread adoption within a theoretical framework? We address this by evaluating the performance of three well established rheological models to characterize ex vivo bovine liver tissues: the Kelvin-Voigt (KV) model as a 2-parameter model, and the standard linear solid (SLS) and Kelvin-Voigt fractional derivative (KVFD) models as 3-parameter models. The assessments were based on the analysis of time domain behavior (using stress relaxation tests) and frequency domain behavior (by measuring shear wave speed (SWS) dispersion). SWS was measured over a wide range of frequency from 1 Hz to 1 kHz using three different tests: (i) harmonic shear tests using a rheometer, (ii) reverberant shear wave (RSW) ultrasound elastography scans, and (iii) RSW optical coherence elastography scans, with each test targeting a distinct frequency range. Our results demonstrated that the KVFD model produces the only mutually consistent rendering of time and frequency domain data for liver. Furthermore, it reduces to a 2-parameter model for liver (correspondingly to a 2-parameter "spring-pot" or power-law model for SWS dispersion) and provides the most accurate predictions of the material viscoelastic behavior in time (>98% accuracy) and frequency (>96% accuracy) domains. STATEMENT OF SIGNIFICANCE: Rheological models are applied in quantifying tissues viscoelastic properties. This study is unique in presenting comprehensive assessments of rheological models.

Feasibility of a Deep Learning approach to estimate Shear Wave Speed using the framework of Reverberant Shear Wave Elastography: A numerical simulation study

Reverberant Shear Wave Elastography (RSWE) is an ultrasound elastography technique that offers great advantages, however, current estimators generate underestimations and time-consuming issues. As well, the involvement of Deep Learning into the medical imaging field with new tools to assess complex problems, makes it a great candidate to serve as a new approach for a RSWE estimator. This work addresses the application of a Deep Neural Network (DNN) for the estimation of Shear Wave Speed (SWS) maps from particle velocity using numerically simulated data. The architecture of the proposed network is based on a U-Net, which works with a custom loss function specifically adopted for the reconstruction task. Four DNNs were trained using four different databases: clean, noisy, acquired at variable frequency, and noisy and acquired at variable frequency data. After the training of the DNNs, the predicted SWS maps were evaluated based on different metrics related to segmentation, regression and similarity of images. The model for clean data showed better results with a Mean Absolute Error (MAE) of 0.011, Mean Square Error(MSE) of 0.001, modified Intersection over Union (mIoU) of 98.4%, Peak Signal to Noise Ratio (PSNR) of 32.925 and a Structural Similarity Index Measure (SSIM) of 0.99, for 250 (size of Testing Sets); while the other models delivered SSIM in the range of 0.87 to 0.96. It was concluded that noisy and clean data could be effectively handled by the model, while the other ones still need enhancement. Clinical Relevance- This work is focused on the application of a Deep Learning approach to accurately asses the Shear Wave Speed in numerical simulations of Reverberant Shear Wave Elastography approach. This novel estimator could be useful for future clinical experiments specially with real time applications to determine the status of living tissue such as detection of malignant or benign tumors located in breast cervix prostate or skin and in the diagnosis of other pathologies such us liver fibrosis.

Mechanical validation of viscoelastic parameters for different interface pressures using the Kelvin-Voigt fractional derivative model

The knowledge of the biomechanical properties of tissues is useful for different applications such as disease diagnosis and treatment monitoring. Reverberant Shear Wave Elastography (RSWE) is an approach that has reduced the restrictions on wave generation to characterize the shear wave velocity over a range of frequencies. This approach is based on the generation of a reverberant field that is generated by the reflections of waves from inhomogeneities and tissue boundaries that exist in the tissue. The Kelvin-Voigt Fractional Derivative model is commonly used to characterize elasticity and viscosity of soft tissue when using shear wave ultrasound elatography. These viscoelastic characteristics can be then validated using mechanical measurements (MM) such as stress relaxation. During RSWE acquisition, the effect of interface pressure, induced by pushing the probe on the skin through the gel pad, on the viscous and elastic characteristics of tissue can be investigated. However, the effect of interface pressure on the validity of the extracted viscous and elastic characteristics was not investigated before. Therefore, the purpose of this study was to compare the estimation of the viscoelastic parameters at different thickness of gel pad against the viscoelastic characteristics obtained from MM. The experiments were conducted in a tissue-mimicking phantom. The results confirm that the relaxed elastic constant (μ₀) can be depreciated. In addition, a higher congruence was found in the viscous parameter (η_α) estimated at 6 and 7 mm. On the other hand, a difference in the order of fractional derivative (α) was found.

Introduction to optical coherence elastography: tutorial

Optical coherence elastography (OCE) has seen rapid growth since its introduction in 1998. The past few decades have seen tremendous advancements in the development of OCE technology and a wide range of applications, including the first clinical applications. This tutorial introduces the basics of solid mechanics, which form the foundation of all elastography methods. We then describe how OCE measurements of tissue motion can be used to quantify tissue biomechanical parameters. We also detail various types of excitation methods, imaging systems, acquisition schemes, and data processing algorithms and how various parameters associated with each step of OCE imaging can affect the final quantitation of biomechanical properties. Finally, we discuss the future of OCE, its potential, and the next steps required for OCE to become an established medical imaging technology.

Assessing corneal cross-linking with reverberant 3D optical coherence elastography

SIGNIFICANCE

Corneal cross-linking (CXL) is a well-known procedure for treating certain eye disorders such as keratoconus. However, characterization of the biomechanical changes in the cornea as a result of this procedure is still under active research. Specifically, there is a clinical need for high-resolution characterization of individual corneal layers.

AIM

A high-resolution elastography method in conjunction with a custom optical coherence tomography system is used to track these biomechanical changes in individual corneal layers. Pre- and post-treatment analysis for both low-dose and high-dose CXL experiments are performed.

APPROACH

A recently developed elastography technique that utilizes the theory of reverberant shear wave fields, with optical coherence tomography as the modality, is applied to pig corneas ex vivo to evaluate elasticity changes associated with corneal CXL. Sets of low-dose and high-dose CXL treatments are evaluated before and after treatments with three pairs of pig corneas per experiment.

RESULTS

The reverberant three-dimensional (3D) optical coherence elastography (OCE) technique can identify increases in elasticity associated with both low-dose and high-dose CXL treatments. There is a notable graphical difference between low-dose and high-dose treatments. In addition, the technique is able to identify which layers of the cornea are potentially affected by the CXL procedure and provides insight into the nonlinearity of the elasticity changes.

CONCLUSIONS

The reverberant 3D OCE technique can identify depth-resolved changes in elasticity of the cornea associated with CXL procedures. This method could be translated to assess and monitor CXL efficacy in various clinical settings.

Optical Coherence Tomography as an Alternative Tool for Evaluating the Effects of Glyphosate on Hybrid Catfish (Clarias gariepinus ×fig Clarias macrocephalus)

•

The use of OCT to evaluate the effects of glyphosate on hybrid catfish was studied.

•

OCT can capture tissue damages in hybrid catfish caused by glyphosate.

•

OCT entails a short preparation time, simple analysis, and chemical free.

•

OCT has potential for real-time or field-based study in aquaculture.

Glyphosate contamination in fresh water is a major problem in agricultural countries. It affects many vital organs in freshwater organisms that are important in the food chain. Hence, the effects of glyphosate on living organism organs are of particular interest. However, several existing techniques for evaluating the effect of glyphosate on aquatic organisms require stained tissue. To study organ tissue with minimal processing, alternative technique is demanded. Here, we investigated the used of optical coherence tomography (OCT) as an alternative tool for ex vivo evaluation of the effect of glyphosate on organ tissues of aquatic organisms, i.e., hybrid catfish (Clarias gariepinus × Clarias macrocephalus). The targeted samples were organ tissues from the brain, gill, and liver of hybrid catfish after glyphosate exposion at concentration of 10 mg L⁻¹ for 24 h. The alteration was then verified by histology, and immunohistochemistry. The study found that all three techniques provide correlated results. We observed that OCT clearly showed damage to the brain and gill tissues of glyphosate-exposed hybrid catfish. However, the alteration in liver tissue was observable but not clear for this low concentration of exposure. The results from histology and immunohistochemistry confirmed the effect of glyphosate on brain, gill, and liver tissues of hybrid catfish. The results suggest that all three techniques could be used to examine the effects of glyphosate exposure in hybrid catfish. However, the choice of a suitable technique depends upon the purpose of the study.

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.

Viscoelastic properties' characterization of corneal stromal models using non-contact surface acoustic wave optical coherence elastography (SAW-OCE)

Viscoelastic characterization of the tissue-engineered corneal stromal model is important for our understanding of the cell behaviors in the pathophysiologic altered corneal extracellular matrix (ECM). The effects of the interactions between stromal cells and different ECM characteristics on the viscoelastic properties during an 11-day culture period were explored. Collagen-based hydrogels seeded with keratocytes were used to replicate human corneal stroma. Keratocytes were seeded at 8 × 10³ cells per hydrogel and with collagen concentrations of 3, 5 and 7 mg/ml. Air-pulse-based surface acoustic wave optical coherence elastography (SAW-OCE) was employed to monitor the changes in the hydrogels' dimensions and viscoelasticity over the culture period. The results showed the elastic modulus increased by 111%, 56% and 6%, and viscosity increased by 357%, 210% and 25% in the 3, 5 and 7 mg/ml hydrogels, respectively. To explain the SAW-OCE results, scanning electron microscope was also performed. The results confirmed the increase in elastic modulus and viscosity of the hydrogels, respectively, arose from increased fiber density and force-dependent unbinding of bonds between collagen fibers. This study reveals the influence of cell-matrix interactions on the viscoelastic properties of corneal stromal models and can provide quantitative guidance for mechanobiological investigations which require collagen ECM with tuneable viscoelastic properties.

Optical coherence tomography

Optical coherence tomography (OCT) is a non-contact method for imaging the topological and internal microstructure of samples in three dimensions. OCT can be configured as a conventional microscope, as an ophthalmic scanner, or using endoscopes and small diameter catheters for accessing internal biological organs. In this Primer, we describe the principles underpinning the different instrument configurations that are tailored to distinct imaging applications and explain the origin of signal, based on light scattering and propagation. Although OCT has been used for imaging inanimate objects, we focus our discussion on biological and medical imaging. We examine the signal processing methods and algorithms that make OCT exquisitely sensitive to reflections as weak as just a few photons and that reveal functional information in addition to structure. Image processing, display and interpretation, which are all critical for effective biomedical imaging, are discussed in the context of specific applications. Finally, we consider image artifacts and limitations that commonly arise and reflect on future advances and opportunities.

Plantar Soft Tissue Characterization Using Reverberant Shear Wave Elastography: A Proof-of-Concept Study

Plantar soft tissue stiffness provides relevant information on biomechanical characteristics of the foot. Therefore, appropriate monitoring of foot elasticity could be useful for diagnosis, treatment or health care of people with complex pathologies such as a diabetic foot. In this work, the reliability of reverberant shear wave elastography (RSWE) applied to plantar soft tissue was investigated. Shear wave speed (SWS) measurements were estimated at the plantar soft tissue at the first metatarsal head, the third metatarsal head and the heel from both feet in five healthy volunteers. Experiments were repeated for a test-retest analysis with and without the use of gel pad using a mechanical excitation frequency range between 400 and 600 Hz. Statistical analysis was performed to evaluate the reliability of the SWS estimations. In addition, the results were compared against those obtained with a commercially available shear wave-based elastography technique, supersonic imaging (SSI). The results indicate a low coefficient of variation for test-retest experiments with gel pad (median: 5.59%) and without gel pad (median: 5.83%). Additionally, the values of the SWS measurements increase at higher frequencies (median values: 2.11 m/s at 400 Hz, 2.16 m/s at 450 Hz, 2.24 m/s at 500 Hz, 2.21 m/s at 550 Hz and 2.31 m/s at 600 Hz), consistent with previous reports at lower frequencies. The SWSs at the plantar soft tissue at the first metatarsal head, third metatarsal head and heel were found be significantly (p<0.05) different, with median values of 2.42, 2.16 and 2.03 m/s, respectively which indicates the ability of the method to differentiate between shear wave speeds at different anatomical locations. The results indicated better elastographic signal-to-noise ratios with RSWE compared to SSI because of the artifacts presented in the SWS generation. These preliminary results indicate that the RSWE approach can be used to estimate the plantar soft tissue elasticity, which may have great potential to better evaluate changes in biomechanical characteristics of the foot.

Practical settings for shear wave speed estimation using the framework of Reverberant Shear Wave Elastography: A numerical simulation study

Reverberant shear wave elastography (RSWE) has become a promising approach to quantifying soft tissues' viscoelastic properties by the propagating shear wave speed (SWS) estimation based on the particle velocity autocorrelation. In this work, three different practical settings were evaluated for the SWS estimation by numerical simulations of an isotropic, homogenous, and elastic medium: first, the 2D representation of the particle velocity, second, the spatial autocorrelation computation, and third, the selection of the curve fitting domain. We conclude that the 2D autocorrelation function using the Wiener-Khinchin theorem provides up to 127 times faster results than traditional autocorrelation methods. Additionally, we state that extracting the magnitude and phase from the Fourier transform of the temporal domain, applying the 2D-autocorrelation on a mobile square window sized at least two wavelengths, and fitting the monotonically decreasing part of the autocorrelation profile's central lobe results in more accurate (13.2% of bias) and precise (5.3% of CV) estimations than other practical settings.Clinical relevance- Affections in soft tissues' biomechanical properties are related to pathologies, such as tumor cancer, muscular degenerative diseases, or fibrosis. These changes are quantified by the SWS and its derived viscoelastic parameters. RSWE is a promising approach for their characterization. In this work, we evaluated alternative elections of practical settings within the methodology. Numerical simulations indicate they lead to faster and more reliable local SWS estimations than conventional settings.

Validations of the microchannel flow model for characterizing vascularized tissues

The microchannel flow model postulates that stress-strain behavior in soft tissues is influenced by the time constants of fluid-filled vessels related to Poiseuille’s law. A consequence of this framework is that changes in fluid viscosity and changes in vessel diameter (through vasoconstriction) have a measurable effect on tissue stiffness. These influences are examined through the theory of the microchannel flow model. Then, the effects of viscosity and vasoconstriction are demonstrated in gelatin phantoms and in perfused tissues, respectively. We find good agreement between theory and experiments using both a simple model made from gelatin and from living, perfused, placental tissue ex vivo.

Reverberant shear wave phase gradients for elastography

Reverberant shear wave fields are produced when multiple sources and multiple reflections establish a complex three-dimensional wave field within an organ. The expected values are assumed to be isotropic across all directions and the autocorrelation functions for velocity are expressed in terms of spherical Bessel functions. These results provide the basis for adroit implementations of elastography from imaging systems that can map out the internal velocity or displacement of tissues during reverberant field excitations. By examining the phase distribution of the reverberant field, additional estimators can be derived. In particular, we demonstrate that the reverberantphase gradientis shown to be proportional to the local value of wavenumber. This phase estimator is less sensitive to imperfections in the reverberant field distribution and requires a smaller support window, relative to earlier estimators based on autocorrelation. Applications are shown in simulations, phantoms, andin vivoliver.

Adaptive incremental method for strain estimation in phase-sensitive optical coherence elastography

We proposed an adaptive incremental method for the cumulative strain estimation in phase-sensitive optical coherence elastography. The method firstly counts the amount of phase noise points by mapping a binary noise map. After the noise threshold value is preset, the interframe interval is adaptively adjusted in terms of the phase noise ratio. Finally, the efficient estimation of cumulative strain is implemented by reducing the cumulative number. Since the level of phase noise is related to the different strain rates in accordance with the speckle decorrelation, the proposed method can estimate the large strains with high computation efficiency as well as signal-to-noise ratio (SNR) enhancement in nonlinear change of sample deformations. Real experiments of visualizing polymerization shrinkage with nonlinear change of deformations were performed to prove the superiority of adaptive incremental method in estimating the large strains. The proposed method expands the practicability of the incremental method in more complex scenes.

Time-of-flight and noise-correlation-inspired algorithms for full-field shear-wave elastography using digital holography

SIGNIFICANCE

Quantitative stiffness information can be a powerful aid for tumor or fibrosis diagnosis. Currently, very promising elastography approaches developed for non-contact biomedical imaging are based on transient shear-waves imaging. Transient elastography offers quantitative stiffness information by tracking the propagation of a wave front. The most common method used to compute stiffness from the acquired propagation movie is based on shear-wave time-of-flight calculations.

AIM

We introduce an approach to transient shear-wave elastography with spatially coherent sources, able to yield full-field quantitative stiffness maps with reduced artifacts compared to typical artifacts observed in time-of-flight.

APPROACH

A noise-correlation algorithm developed for passive elastography is adapted to spatially coherent narrow or any band sources. This noise-correlation-inspired (NCi) method is employed in parallel with a classic time-of-flight approach. Testing is done on simulation images, experimental validation is conducted with a digital holography setup on controlled homogeneous samples, and full-field quantitative stiffness maps are presented for heterogeneous samples and ex-vivo biological tissues.

RESULTS

The NCi approach is first validated on simulations images. Stiffness images processed by the NCi approach on simulated inclusions display significantly less artifacts than with a time-of-flight reconstruction. The adaptability of the NCi algorithm to narrow or any band shear-wave sources was tested successfully. Experimental testing on homogeneous samples demonstrates similar values for both the time-of-flight and the NCi approach. Soft inclusions in agarose sample could be resolved using the NCi method and feasibility on ex-vivo biological tissues is presented.

CONCLUSIONS

The presented NCi approach was successful in computing quantitative full-field stiffness maps with narrow and broadband source signals on simulation and experimental images from a digital holography setup. Results in heterogeneous media show that the NCi approach could provide stiffness maps with less artifacts than with time-of-flight, demonstrating that a NCi algorithm is a promising approach for shear-wave transient elastography with spatially coherent sources.

Measurement of In Vivo Biomechanical Changes Attributable to Epithelial Removal in Keratoconus Using a Noncontact Tonometer

PURPOSE

To compare the biomechanical properties of the cornea after epithelial removal in eyes with keratoconus undergoing corneal cross-linking.

METHODS

Prospective interventional case series at a university hospital tertiary referral center. Corneal biomechanical properties were measured in patients with keratoconus undergoing corneal cross-linking, immediately before and after epithelial debridement by using a dynamic ultrahigh-speed Scheimpflug camera equipped with a noncontact tonometer.

RESULTS

The study comprised 45 eyes of 45 patients with a mean age of 19.6 ± 4.9 years (range 14-34). The cornea was found to be 23.7 ± 15.7 μm thinner after epithelial removal (P < 0.01). Corneal stiffness was reduced after epithelial removal as demonstrated by a significant decrease of parameters such as stiffness parameter A1 (12.31, P < 0.01), stiffness parameter-highest concavity (2.25, P < 0.01), A1 length (0.13 mm, P = 0.04), highest concavity radius of curvature (0.26 mm, P = 0.01), highest concavity time (0.22 ms, P = 0.04) and an increase in A1 velocity (-0.01 m/s, P = 0.01), A1 deformation amplitude (-0.03 mm, P ≤ 0.01), A1 deflection length (-0.32 mm, P < 0.01), A2 deformation amplitude (-0.03 mm, P = 0.01), and A2 deflection length (-1.00 mm, P < 0.01). There were no significant differences in biomechanical intraocular pressure (0.15 mm Hg, P = 0.78), deformation amplitude (0.03, P = 0.54), maximum inverse radius (-0.01 mm, P = 0.57), and whole eye movement length (-0.02 mm, P = 0.12).

CONCLUSIONS

Dynamic ultrahigh-speed Scheimpflug camera equipped with a noncontact tonometer offers an alternative method for in vivo measurements of the epithelial layer's contribution to corneal biomechanical properties. Our results suggest that corneal epithelium may play a more significant role in corneal biomechanical properties in patients with keratoconus than previously described.

Ultrasound Shear Wave Elastography and Transient Optical Coherence Elastography: Side-by-Side Comparison of Repeatability and Accuracy

Objective: We compare the repeatability and accuracy of ultrasound shear wave elastography (USE) and transient optical coherence elastography (OCE). Methods: Elastic wave speed in gelatin phantoms and chicken breast was measured with USE and OCE and compared with uniaxial mechanical compression testing. Intra- and Inter-repeatability were analyzed using Bland-Altman plots and intraclass correlation coefficients (ICC). Results: OCE and USE differed from uniaxial testing by a mean absolute percent error of 8.92% and 16.9%, respectively, across eight phantoms of varying stiffness. Upper and lower limits of agreement for intrasample repeatability for USE and OCE were ±0.075 m/s and −0.14 m/s and 0.13 m/s, respectively. OCE and USE both had ICCs of 0.9991. In chicken breast, ICC for USE was 0.9385 and for OCE was 0.9924. Conclusion: OCE and USE can detect small speed changes and give comparable measurements. These measurements correspond well with uniaxial testing.

Compressional Optical Coherence Elastography of the Cornea

Assessing the biomechanical properties of the cornea is crucial for detecting the onset and progression of eye diseases. In this work, we demonstrate the application of compression-based optical coherence elastography (OCE) to measure the biomechanical properties of the cornea under various conditions, including validation in an in situ rabbit model and a demonstration of feasibility for in vivo measurements. Our results show a stark increase in the stiffness of the corneas as IOP was increased. Moreover, UV-A/riboflavin corneal collagen crosslinking (CXL) also dramatically increased the stiffness of the corneas. The results were consistent across 4 different scenarios (whole CXL in situ, partial CXL in situ, whole CXL in vivo, and partial CXL in vivo), emphasizing the reliability of compression OCE to measure corneal biomechanical properties and its potential for clinical applications.

Stress-Strain Index Map: A New Way to Represent Corneal Material Stiffness

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.

An Intraocular Pressure Measurement Technique Based on Acoustic Radiation Force Using an Ultrasound Transducer: A Feasibility Study

High intraocular pressure (IOP) is one of the major risk factors for glaucoma, and thus accurate IOP measurements should be performed to diagnose and treat glaucoma early. In this study, a novel technique for measuring the IOP based on acoustic radiation force was proposed, and its potential was experimentally demonstrated. The proposed technique uses the acoustic radiation force to generate axial displacement on the ocular surface while simultaneously measuring the degree of deformation. In order to verify that the ocular displacement induced by the acoustic radiation force is related to the IOP, the experiment was conducted by fabricating a 5 MHz single element transducer and gelatin phantoms with different stiffness values. Our experimental results show that there is a close relationship between the ocular displacement by the acoustic radiation force and the IOP obtained by a commercial tonometer. Therefore, the proposed acoustic radiation force technique can be a promising candidate for measuring the IOP.

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.

In Vivo Human Corneal Shear-wave Optical Coherence Elastography

SIGNIFICANCE

A novel imaging technology, dynamic optical coherence elastography (OCE), was adapted for clinical noninvasive measurements of corneal biomechanics.

PURPOSE

Determining corneal biomechanical properties is a long-standing challenge. Elasticity imaging methods have recently been developed and applied for clinical evaluation of soft tissues in cancer detection, atherosclerotic plaque evaluation, surgical guidance, and more. Here, we describe the use of dynamic OCE to characterize mechanical wave propagation in the human cornea in vivo, thus providing a method for clinical determination of corneal biomechanical properties.

METHODS

High-resolution phase-sensitive optical coherence tomography imaging was combined with microliter air-pulse tissue stimulation to perform dynamic elasticity measurements in 18 eyes of nine participants. Low-pressure (0.1 mmHg), spatiotemporally discreet (150 μm, 800 μs) tissue stimulation produced submicron-scale tissue deformations that were measured at multiple positions over a 1-mm2 area. Surface wave velocity was measured and used to determine tissue stiffness. Elastic wave propagation velocity was measured and evaluated as a function of IOP and central corneal thickness.

RESULTS

Submicron corneal surface displacement amplitude (range, 0.005 to 0.5 μm) responses were measured with high sensitivity (0.24 nm). Corneal elastic wave velocity ranged from 2.4 to 4.2 m/s (mean, 3.5; 95% confidence interval, 3.2 to 3.8 m/s) and was correlated with central corneal thickness (r = 0.64, P < .001) and IOP (r = 0.52, P = .02).

CONCLUSIONS

Phase-sensitive optical coherence tomography imaging combined with microliter air-pulse mechanical tissue stimulation has sufficient detection sensitivity to observe submicron elastic wave propagation in corneal tissue. These measurements enable in vivo corneal stiffness determinations that will be further studied for use with disease detection and for monitoring clinical interventions.

4D deep learning for real-time volumetric optical coherence elastography

PURPOSE

Elasticity of soft tissue provides valuable information to physicians during treatment and diagnosis of diseases. A number of approaches have been proposed to estimate tissue stiffness from the shear wave velocity. Optical coherence elastography offers a particularly high spatial and temporal resolution. However, current approaches typically acquire data at different positions sequentially, making it slow and less practical for clinical application.

METHODS

We propose a new approach for elastography estimations using a fast imaging device to acquire small image volumes at rates of 831 Hz. The resulting sequence of phase image volumes is fed into a 4D convolutional neural network which handles both spatial and temporal data processing. We evaluate the approach on a set of image data acquired for gelatin phantoms of known elasticity.

RESULTS

Using the neural network, the gelatin concentration of unseen samples was predicted with a mean error of 0.65 ± 0.81 percentage points from 90 subsequent volumes of phase data only. We achieve a data acquisition and data processing time of under 12 ms and 22 ms, respectively.

CONCLUSIONS

We demonstrate direct volumetric optical coherence elastography from phase image data. The approach does not rely on particular stimulation or sampling sequences and allows the estimation of elastic tissue properties of up to 40 Hz.

Confocal air-coupled ultrasonic optical coherence elastography probe for quantitative biomechanics

We present an air-coupled ultrasonic radiation force probe co-focused with a phase-sensitive optical coherence tomography (OCT) system for quantitative wave-based elastography. A custom-made 1 MHz spherically focused piezoelectric transducer with a concentric 10 mm wide circular opening allowed for confocal micro-excitation of waves and phase-sensitive OCT imaging. Phantom studies demonstrated the capabilities of this probe to produce quasi-harmonic excitation up to 4 kHz for generation of elastic waves. Experimental results in ocular tissues showed highly detailed 2D and 3D elasticity mapping using this approach with great potential for clinical translation.

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.

In vivo noninvasive measurement of spatially resolved corneal elasticity in human eyes using Lamb wave optical coherence elastography

Current elastography techniques are limited in application to accurately assess spatially resolved corneal elasticity in vivo for human eyes. The air-puff optical coherence elastography (OCE) with an eye motion artifacts correction algorithm is developed to distinguish the in vivo cornea vibration from the eye motion and visualize the Lamb wave propagation clearly in healthy subjects. Based on the Lamb wave model, the phase velocity dispersion curve in the high-frequency is calculated to obtain spatially resolved corneal elasticity accurately with high repeatability. It is found that the corneal elasticity has regional variations and is correlated with intraocular pressure, which suggests that the method has the potential to provide noninvasive measurement of spatially resolved corneal elasticity in clinical practice.

Nearly-incompressible transverse isotropy (NITI) of cornea elasticity: model and experiments with acoustic micro-tapping OCE

The cornea provides the largest refractive power for the human visual system. Its stiffness, along with intraocular pressure (IOP), are linked to several pathologies, including keratoconus and glaucoma. Although mechanical tests can quantify corneal elasticity ex vivo, they cannot be used clinically. Dynamic optical coherence elastography (OCE), which launches and tracks shear waves to estimate stiffness, provides an attractive non-contact probe of corneal elasticity. To date, however, OCE studies report corneal moduli around tens of kPa, orders-of-magnitude less than those (few MPa) obtained by tensile/inflation testing. This large discrepancy impedes OCE's clinical adoption. Based on corneal microstructure, we introduce and fully characterize a nearly-incompressible transversely isotropic (NITI) model depicting corneal biomechanics. We show that the cornea must be described by at least two shear moduli, contrary to current single-modulus models, decoupling tensile and shear responses. We measure both as a function of IOP in ex vivo porcine cornea, obtaining values consistent with both tensile and shear tests. At pressures above 30 mmHg, the model begins to fail, consistent with non-linear changes in cornea at high IOP.