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

Convolutional Neural Networks Enable Direct Strain Estimation in Quasistatic Optical Coherence Elastography

Assessing the biomechanical properties of tissues can provide important information for disease diagnosis and therapeutic monitoring. Optical coherence elastography (OCE) is an emerging technology for measuring the biomechanical properties of tissues. Clinical translation of this technology is underway, and steps are being implemented to streamline data collection and processing. OCE data can be noisy, data processing can require significant manual tuning, and a single acquisition may contain gigabytes of data. In this work, we introduce a convolutional neural network-based method to translate raw OCE phase data to strain for quasistatic OCE that is ~40X faster than the conventional least squares approach by bypassing many intermediate data processing steps. The results suggest that a machine learning approach may be a valuable tool for fast, efficient, and accurate extraction of biomechanical information from raw OCE data.

Toward quantitative X-ray elastography of coronary arteries using flexural pulse waves

Dynamic elastography uses an imaging system to visualize the propagation of elastic waves, the speed of which is directly related to the elasticity felt by palpation. Very few studies have focused on X-ray elastography because of the technical challenges it poses: a planar image of an integration volume at a very slow sampling rate. We demonstrate that tracking a slow elastic wave guided along a one-dimensional structure could provide a possible solution. The recently discovered flexural pulse wave, which is naturally generated by heartbeats and propagates along arteries, is the perfect candidate for X-ray elastography. As it reflects the cardiovascular health of patients, arterial elasticity is a biomarker of high clinical interest. We first validate the method by measuring the elasticity in artery phantoms using X-ray. We then move on to data obtained in vivo on coronary arteries during a routine angiography examination. During coronary angiography, a catheter is used to inject an X-ray contrast dye into the patient's aorta. X-rays are then taken as the dye spreads through the coronary arteries. It shows the movement of the coronary arteries for a few seconds and provides an opportunity to follow the natural flexural pulse waves. The obtained Young's moduli for two patients are E = 38 ± 30 kPa and E = 38 ± 28 kPa, respectively. These preliminary results are expected to pave the way for X-ray elastography.

Automated chick gender determination using optical coherence tomography and deep learning

Chick gender classification is crucial for optimizing poultry production, yet traditional methods such as vent sexing and ultrasound remain limited by human expertise, labor intensity, and insufficient resolution. This study introduces a novel approach that integrates Optical Coherence Tomography (OCT) and deep learning to enable high-resolution, non-invasive chick sexing. Unlike conventional imaging techniques, OCT provides micrometer-scale visualization of cloacal structures, allowing precise differentiation between male and female chicks based on internal anatomical markers. We developed a custom convolutional neural network (CNN) optimized for OCT data, incorporating asymmetric image resizing and enhanced feature extraction to improve classification accuracy. Our model achieved 79 % accuracy, outperforming conventional architectures such as Inception (63 %) and VGG-16 (74 %), highlighting the importance of a tailored, domain-specific model. This is the first study to integrate OCT with deep learning for automated chick sexing, demonstrating a scalable, real-time alternative to expert-dependent vent sexing. With further advancements in imaging and machine learning, our approach has the potential to transform chick sexing in commercial hatcheries, reducing reliance on skilled labor while enhancing classification efficiency and precision.

Rapid Stiffness Mapping in Soft Biologic Tissues With Micrometer Resolution Using Optical Multifrequency Time-Harmonic Elastography

Rapid mapping of the mechanical properties of soft biological tissues from light microscopy to macroscopic imaging can transform fundamental biophysical research by providing clinical biomarkers to complement in vivo elastography. This work introduces superfast optical multifrequency time-harmonic elastography (OMTHE) to remotely encode surface and subsurface shear wave fields for generating maps of tissue stiffness with unprecedented detail resolution. OMTHE rigorously exploits the space-time propagation characteristics of multifrequency time-harmonic waves to address current limitations of biomechanical imaging and elastography. Key solutions are presented for stimulation, wave decoding, and stiffness reconstruction of shear waves at multiple harmonic frequencies, all tuned to provide consistent stiffness values across resolutions from microns to millimeters. OMTHE's versatility is demonstrated by simulations, phantoms, Bacillus subtilis biofilms, zebrafish embryos and adult zebrafish, reflecting the diversity of biological systems from a mechanics perspective. By zooming in on stiffness details from coarse to finer scales, OMTHE has the potential to advance mechanobiology and offers a way to perform biomechanics-based tissue histology that consistently matches in vivo time-harmonic elastography in patients.

Integrated Difference Autocorrelation: A Novel Approach to Estimate Shear Wave Speed in the Presence of Compression Waves

OBJECTIVE

In shear wave elastography (SWE), the aim is to measure the velocity of shear waves, however unwanted compression waves and bulk tissue motion pose challenges in evaluating tissue stiffness. Conventional approaches often struggle to discriminate between shear and compression waves, leading to inaccurate shear wave speed (SWS) estimation. In this study, we propose a novel approach known as the integrated difference autocorrelation (IDA) estimator to accurately estimate reverberant SWS in the presence of compression waves and noise.

METHODS

The IDA estimator, unlike conventional techniques, computes the subtraction of velocity between neighboring particles, effectively minimizing the impact of long wavelength compression waves and other wide-area movements such as those caused by respiration. We evaluated the effectiveness of IDA by: (1) using k-Wave simulations of a branching cylinder in a soft background, (2) using ultrasound elastography on a breast phantom, (3) using ultrasound elastography in the human liver-kidney region, and (4) using magnetic resonance elastography (MRE) on a brain phantom.

RESULTS

By applying IDA to unfiltered contaminated wave fields of simulation and elastography experiments, the estimated SWSs are in good agreement with the ground truth values (i.e., less than 2% error for the simulation, 9% error for ultrasound elastography of the breast phantom and 19% error for MRE).

CONCLUSION

Our results demonstrate that IDA accurately estimates SWS, revealing the existence of a lesion, even in the presence of strong compression waves.

SIGNIFICANCE

IDA exhibits consistency in SWS estimation across different modalities and excitation scenarios, highlighting its robustness and potential clinical utility.

Shear Wave Optical Coherence Elastography and Structural Analysis of the Postnatal Mouse Cornea Through Development

The cornea is a transparent lens at the forefront of the eye, serving as a structural barrier that protects the eye and provides the majority of the eye's refractive power. In this study, we utilized the postnatal mouse eye to characterize corneal structure and biomechanics. Between postnatal day (PN) 6 and PN24, we observed that elastic wave speeds are highest at PN6, gradually decrease through PN24, and then start to rise in adulthood (6 months). We found that corneal thickness is uncoupled from elastic wave speed and that the content and organization of the cornea primarily influence its mechanical properties. To the best of our knowledge, this work is the first detailed assessment of the postnatal mouse cornea's structure and biomechanics and warrants further investigation into the dynamic properties of the postnatal eye.

Asynchronous, semi-reverberant elastography

Optical coherence elastography measures elasticity-a property correlated with pathologies such as tumors due to fibrosis, atherosclerosis due to heterogeneous plaque composition, and ocular diseases such as keratoconus and glaucoma. Wave-based elastography, including reverberant elastography, leverages the properties of shear waves traveling through tissue primarily to infer shear modulus. These methods have already seen significant development over the past decade. However, existing implementations in OCT require robust synchronization of shear wave excitation with imaging, complicating widespread clinical adoption. We present a method for complete recovery of the harmonic shear wave field in an asynchronous, conventional frame-rate, raster-scanning OCT system by modeling raster-scanning as an amplitude modulation of the displacement field. This technique recovers the entire spatially and temporally coherent complex valued shear wave field from just two B-scans, while reducing the time scale for sensitivity to motion from minutes to tens of milliseconds. To the best of our knowledge, this work represents the first successful demonstration of reverberant elastography on a human subject in vivo with a conventional frame-rate, raster-scanning OCT system, greatly expanding opportunity for widespread translation.

Phase-Domain Photoacoustic Mechanical Imaging for Quantitative Elastography and Viscography

The role and importance of mechanical properties of cells and tissues in pathophysiological processes have widely been acknowledged. However, current elastography techniques most based on transverse elastic waves, diminish the translation of wave speed into elastic modulus due to its limited wave propagation direction. Here, we propose phase-domain photoacoustic mechanical imaging (PD-PAMI), leveraging the initial time and phase response characteristics of an omnidirectional photoacoustic elastic wave to quantitatively extract elastic and viscous moduli. Theoretical simulations and experiment on tissue-mimicking phantoms with different levels of viscoelastic properties were conducted to validate the approach with a precision in elasticity and viscosity estimation of 4.6% and 6.6%, respectively. The trans-scale viscoelasticity mappings over three length scales-covering cell, tissue section, and in vivo organ, were provided to demonstrate the scalability of the technique with different implementations of PD-PAMI. Experiments on animal models of breast tumour and atherosclerosis reveal that PD-PAMI technique enables effective monitoring of the viscoelastic parameters for examinations of the diseases involved with the variations in collagen or lipid composition and in inflammation level. PD-PAMI technique opens new perspectives of conventional PA imaging and provides new technical way for biomechanical imaging, prefiguring potential clinical applications in mechanopathology-involved disease diagnosis.

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.