Optical coherence elastography for tissue characterization: a review

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.

Brillouin microscopy monitors rapid responses in subcellular compartments

Measurements and imaging of the mechanical response of biological cells are critical for understanding the mechanisms of many diseases, and for fundamental studies of energy, signal and force transduction. The recent emergence of Brillouin microscopy as a powerful non-contact, label-free way to non-invasively and non-destructively assess local viscoelastic properties provides an opportunity to expand the scope of biomechanical research to the sub-cellular level. Brillouin spectroscopy has recently been validated through static measurements of cell viscoelastic properties, however, fast (sub-second) measurements of sub-cellular cytomechanical changes have yet to be reported. In this report, we utilize a custom multimodal spectroscopy system to monitor for the very first time the rapid viscoelastic response of cells and subcellular structures to a short-duration electrical impulse. The cytomechanical response of three subcellular structures - cytoplasm, nucleoplasm, and nucleoli - were monitored, showing distinct mechanical changes despite an identical stimulus. Through this pioneering transformative study, we demonstrate the capability of Brillouin spectroscopy to measure rapid, real-time biomechanical changes within distinct subcellular compartments. Our results support the promising future of Brillouin spectroscopy within the broad scope of cellular biomechanics.

Disruption of Fuz in mouse embryos generates hypoplastic hindbrain development and reduced cranial nerve ganglia

BACKGROUND

The brain and spinal cord formation is initiated in the earliest stages of mammalian pregnancy in a highly organized process known as neurulation. Environmental or genetic interferences can impair neurulation, resulting in clinically significant birth defects known collectively as neural tube defects. The Fuz gene encodes a subunit of the CPLANE complex, a macromolecular planar polarity effector required for ciliogenesis. Ablation of Fuz in mouse embryos results in exencephaly and spina bifida, including dysmorphic craniofacial structures due to defective cilia formation and impaired Sonic Hedgehog signaling.

RESULTS

We demonstrate that knocking Fuz out during embryonic mouse development results in a hypoplastic hindbrain phenotype, displaying abnormal rhombomeres with reduced length and width. This phenotype is associated with persistent reduction of ventral neuroepithelial stiffness in a notochord adjacent area at the level of the rhombomere 5. The formation of cranial and paravertebral ganglia is also impaired in these embryos.

CONCLUSIONS

This study reveals that hypoplastic hindbrain development, identified by abnormal rhombomere morphology and persistent loss of ventral neuroepithelial stiffness, precedes exencephaly in Fuz ablated murine mutants, indicating that the gene Fuz has a critical function sustaining normal neural tube development and neuronal differentiation.

In vivo dynamic volumetric imaging of mouse testis and epididymis with optical coherence tomography†

The implementation of live imaging in reproductive research is crucial for studying the physiological dynamics. Sperm transport is a highly dynamic process regulated by tubular contractions and luminal flows within the male reproductive tract. However, due to the lack of imaging techniques to capture these dynamics in vivo, there is little information on the physiological and biomechanical regulation of sperm transport through the male reproductive tract. Here, we present a functional in vivo imaging approach using optical coherence tomography, enabling live, label-free, depth-resolved, three-dimensional, high-resolution visualization of the mouse testis and epididymis. With this approach, we spatiotemporally captured tubular contractility in mouse testis and epididymis, as well as microstructures of these reproductive organs. Our findings demonstrated that the contraction frequency varies significantly depending on the epididymal regions, suggesting the spatial regulation of epididymal contractility. Furthermore, we implemented quantitative measurements of the contraction wave and luminal transport through the epididymal duct, revealing the physiological dynamics within the male reproductive tract. The results show that the contraction wave propagates along the epididymal duct and the wave propagation velocity was estimated in vivo. In conclusion, this is the first study to develop in vivo dynamic volumetric imaging of the male reproductive tract, which allows for quantitative analysis of the dynamics associated with sperm transport. This study sets a platform for various studies investigating normal and abnormal male reproductive physiology as well as the pharmacological and environmental effects on reproductive functions in mouse models, ultimately contributing to a comprehensive understanding of male reproductive disorders.

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%.

Multifunctional adaptive optics optical coherence tomography allows cellular scale reflectometry, polarimetry, and angiography in the living human eye

Clinicians are unable to detect glaucoma until substantial loss or dysfunction of retinal ganglion cells occurs. To this end, novel measures are needed. We have developed an optical imaging solution based on adaptive optics optical coherence tomography (AO-OCT) to discern key clinical features of glaucoma and other neurodegenerative diseases at the cellular scale in the living eye. Here, we test the feasibility of measuring AO-OCT-based reflectance, retardance, optic axis orientation, and angiogram at specifically targeted locations in the living human retina and optic nerve head. Multifunctional imaging, combined with focus stacking and global image registration algorithms, allows us to visualize cellular details of retinal nerve fiber bundles, ganglion cell layer somas, glial septa, superior vascular complex capillaries, and connective tissues. These are key histologic features of neurodegenerative diseases, including glaucoma, that are now measurable in vivo with excellent repeatability and reproducibility. Incorporating this noninvasive cellular-scale imaging with objective measurements will significantly enhance existing clinical assessments, which is pivotal in facilitating the early detection of eye disease and understanding the mechanisms of neurodegeneration.

Photothermal optical coherence tomography for 3D live cell detection and mapping

Imaging cells in their 3D environment with molecular specificity is important to cell biology study. Widely used microscopy techniques, such as confocal microscopy, have limited imaging depth when probing cells in optically scattering media. Optical coherence tomography (OCT) can provide millimeter-level depth for imaging of highly scattering media but lacks the contrast to distinguish cells from extracellular matrix or to distinguish between different types of cells. Photothermal OCT (PT-OCT) is a promising technique to obtain molecular contrast at the imaging scale of OCT. Here, we report PT-OCT imaging of live, nanoparticle-labeled cells in 3D. In particular, we demonstrate detection and mapping of single cell in 3D without causing call death, and show the feasibility of 3D cell mapping through optical scattering media. This work presents live cell detection and mapping at an imaging scale that complements the major microscopy techniques, which is potentially useful to study cells in their 3D native or culture environment.

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.

Non-destructive Characterization of Skeletal Muscle Extracellular Matrix Morphology by Combining Optical Coherence Tomography (OCT) Imaging with Tissue Clearing

Histology is an essential step to visualize and analyze the microstructure of any biological tissue; however, histological processing is often irreversible, and histological samples are unable to be imaged or tested further. In this work, a novel non-destructive protocol is proposed for morphological analysis of skeletal muscles, combining Optical Coherence Tomography (OCT) imaging with Tissue Clearing. Imaging combining OCT and Propylene Glycol (PG) as a tissue-clearing agent, was performed on rat tail and extensor digitorum longus (EDL) muscle. The results show that the extracellular matrix morphology of skeletal muscles, including muscular fibers and the whole microstructure architecture were clearly identified. PG improved OCT imaging as measured by image quality metric Contrast Per Pixel CPP (increases by 3.9%), Naturalness Image Quality Evaluator NIQE (decreases by 23%), and Volume of Interest VOI size (higher for CPP and lower for NIQE values). The tendon microstructure was observed with less precision, as collagen fibers could not be clearly detected. The reversibility of the optical effects of the PG on the immersed tissue (in a Phosphate-Buffered Saline solution) was studied comparing native and rehydrated OCT image acquisition from a single EDL sample. Optical properties and microstructure visibility (CPP and NIQE) have been recovered to 99% of the native sample values. Moreover, clearing process caused shrinkage of the tissue recovered to 86% of the original width. Future work will aim to employ the proposed experimental protocol to identify the local mechanical properties of biological tissues.

Disruption of Fuz in mouse embryos generates hypoplastic hindbrain development and reduced cranial nerve ganglia

The formation of the brain and spinal cord is initiated in the earliest stages of mammalian pregnancy in a highly organized process known as neurulation. Convergent and extension movements transforms a flat sheet of ectodermal cells into a narrow and elongated line of neuroepithelia, while a major source of Sonic Hedgehog signaling from the notochord induces the overlying neuroepithelial cells to form two apposed neural folds. Afterward, neural tube closure occurs by synchronized coordination of the surface ectoderm and adjacent neuroepithelial walls at specific axial regions known as neuropores. Environmental or genetic interferences can impair neurulation resulting in neural tube defects. The Fuz gene encodes a subunit of the CPLANE complex, which is a macromolecular planar polarity effector required for ciliogenesis. Ablation of Fuz in mouse embryos results in exencephaly and spina bifida, including dysmorphic craniofacial structures due to defective cilia formation and impaired Sonic Hedgehog signaling. In this work, we demonstrate that knocking Fuz out during embryonic mouse development results in a hypoplastic hindbrain phenotype, displaying abnormal rhombomeres with reduced length and width. This phenotype is associated with persistent loss of ventral neuroepithelial stiffness, in a notochord adjacent area at the level of the rhombomere 5, preceding the development of exencephaly in Fuz ablated mutants. The formation of cranial and paravertebral ganglia is also impaired in these embryos, indicating that Fuz has a critical function sustaining normal neural tube development and neuronal differentiation.

SIGNIFICANCE STATEMENT

Neural tube defects (NTDs) are a common cause of disability in children, representing the second most common congenital structural malformation in humans following only congenital cardiovascular malformations. NTDs affect approximately 1 to 2 pregnancies per 1000 births every year worldwide, when the mechanical forces folding the neural plate fails to close at specific neuropores located anteriorly (cranial) or posteriorly (caudal) along the neural tube, in a process known as neurulation, which happens throughout the third and fourth weeks of human pregnancy.

The impact of amplitude modulation frequency in harmonic motion imaging on inclusion characterization

OBJECTIVE

Ultrasound elasticity imaging techniques aim to provide a non-invasive characterization of tissue mechanical properties to detect pathological changes and monitor disease progression. Harmonic motion imaging (HMI) is an ultrasound-based elasticity imaging technique that utilizes an oscillatory acoustic radiation force to induce localized displacements and estimate relative tissue stiffness. Previous studies have applied a low amplitude modulation (AM) frequency of 25 or 50 Hz in HMI to assess the mechanical properties of different tissue types. In this study, we investigate the dependence of AM frequency in HMI and whether the frequency can be adjusted based on the size and mechanical properties of the underlying medium for enhanced image contrast and inclusion detection.

METHODS

A tissue-mimicking phantom with embedded inclusions at different sizes and stiffnesses was imaged within a range of AM frequencies from 25 to 250 Hz at 25-Hz step size.

DISCUSSION

The AM frequency at which the maximum contrast and CNR are achieved depends on the size and stiffness of the inclusions. A general trend shows that contrast and CNR peak at higher frequencies for smaller inclusions. In addition, for some inclusions with the same size but different stiffnesses, the optimized AM frequency increases with the stiffness of the inclusion. Nevertheless, there is a shift between the frequencies at which the contrast peaks and those with maximum CNR. Finally, in agreement with the phantom findings, imaging an ex-vivo human specimen with a 2.7-cm breast tumor at a range of AM frequencies showed that the highest contrast and CNR are achieved at the AM frequency of 50 Hz.

CONCLUSION

These findings indicate that the AM frequency can be optimized in different applications of HMI, especially in the clinic, for improved detection and characterization of tumors with different geometries and mechanical properties.

Corneal Surface Wave Propagation Associated with Intraocular Pressures: OCT Elastography Assessment in a Simplified Eye Model

Assessing corneal biomechanics in vivo has long been a challenge in the field of ophthalmology. Despite recent advances in optical coherence tomography (OCT)-based elastography (OCE) methods, controversy remains regarding the effect of intraocular pressure (IOP) on mechanical wave propagation speed in the cornea. This could be attributed to the complexity of corneal biomechanics and the difficulties associated with conducting in vivo corneal shear-wave OCE measurements. We constructed a simplified artificial eye model with a silicone cornea and controllable IOPs and performed surface wave OCE measurements in radial directions (54-324°) of the silicone cornea at different IOP levels (10-40 mmHg). The results demonstrated increases in wave propagation speeds (mean ± STD) from 6.55 ± 0.09 m/s (10 mmHg) to 9.82 ± 0.19 m/s (40 mmHg), leading to an estimate of Young's modulus, which increased from 145.23 ± 4.43 kPa to 326.44 ± 13.30 kPa. Our implementation of an artificial eye model highlighted that the impact of IOP on Young's modulus (ΔE = 165.59 kPa, IOP: 10-40 mmHg) was more significant than the effect of stretching of the silicone cornea (ΔE = 15.79 kPa, relative elongation: 0.98-6.49%). Our study sheds light on the potential advantages of using an artificial eye model to represent the response of the human cornea during OCE measurement and provides valuable insights into the impact of IOP on wave-based OCE measurement for future in vivo corneal biomechanics studies.

Acoustic radiation force optical coherence elastography: A preliminary study on biomechanical properties of trabecular meshwork

Evaluating biomechanical properties of trabecular meshwork (TM) is of great significance for understanding the mechanism of aqueous humor circulation and its relationship to some eye diseases such as glaucoma; however, there is almost no relevant study due to the lack of clinical measurement tool. In this paper, an acoustic radiation force optical coherence elastography (ARF-OCE) system is developed with the advantages of noninvasive detection, high resolution, high sensitivity, and high-speed imaging, by which elastic modulus of the porcine and human TMs is accurately quantified. As the first OCE imaging of TM, our study demonstrates that ARF-OCE may be an effective approach to advance the research of diseases related to aqueous humor circulation.

In vivo corneal elastography: A topical review of challenges and opportunities

Clinical measurement of corneal biomechanics can aid in the early diagnosis, progression tracking, and treatment evaluation of ocular diseases. Over the past two decades, interdisciplinary collaborations between investigators in optical engineering, analytical biomechanical modeling, and clinical research has expanded our knowledge of corneal biomechanics. These advances have led to innovations in testing methods (ex vivo, and recently, in vivo) across multiple spatial and strain scales. However, in vivo measurement of corneal biomechanics remains a long-standing challenge and is currently an active area of research. Here, we review the existing and emerging approaches for in vivo corneal biomechanics evaluation, which include corneal applanation methods, such as ocular response analyzer (ORA) and corneal visualization Scheimpflug technology (Corvis ST), Brillouin microscopy, and elastography methods, and the emerging field of optical coherence elastography (OCE). We describe the fundamental concepts, analytical methods, and current clinical status for each of these methods. Finally, we discuss open questions for the current state of in vivo biomechanics assessment techniques and requirements for wider use that will further broaden our understanding of corneal biomechanics for the detection and management of ocular diseases, and improve the safety and efficacy of future clinical practice.

Reconstructing the Spatial Distribution of the Relative Shear Modulus in Quasi-static Ultrasound Elastography: Plane Stress Analysis

Quasi-static ultrasound elastography (QSUE) is an imaging technique that mainly provides axial strain maps of tissues when the latter are subjected to compression. In this article, a method for reconstructing the relative shear modulus distribution within a linear elastic and isotropic medium, in QSUE, is introduced. More specifically, the plane stress inverse problem is considered. The proposed method is based on the variational formulation of the equilibrium equations and on the choice of adapted discretization spaces, and only requires displacement fields in the analyzed media to be determined. Results from plane stress and 3-D numerical simulations, as well as from phantom experiments, showed that the method is able to reconstruct the different regions within a medium, with shear modulus contrasts that unambiguously reveal whether inclusions are stiffer or softer than the surrounding material. More specifically, for the plane stress simulations, inclusion-to-background modulus ratios were found to be very accurately estimated, with an error lower than 3%. For the 3-D simulations, for which the plane stress conditions are no longer satisfied, these ratios were, as expected, less accurate, with an error that remained lower than 10% for two of the three cases analyzed but was around 34% for the last case. Concerning the phantom experiments, a comparison with a shear wave elastography technique from a clinical ultrasound scanner was also made. Overall, the inclusion-to-background shear modulus ratios obtained with our approach were found to be closer to those given by the phantom manufacturer than the ratios provided by the clinical system.

Motion-tracking Brillouin microscopy for in-vivo corneal biomechanics mapping

Corneal biomechanics play a critical role in maintaining corneal shape and thereby directly influence visual acuity. However, direct corneal biomechanical measurement in-vivo with sufficient accuracy and a high spatial resolution remains an open need. Here, we developed a three-dimensional (3D) motion-tracking Brillouin microscope for in-vivo corneal biomechanics mapping. The axial tracking utilized optical coherence tomography, which provided a tracking accuracy better than 3 µm. Meanwhile, 10 µm lateral tracking was achieved by tracking pupils with digital image processing. The 3D tracking enabled reconstruction of depth-dependent Brillouin distribution with a high spatial resolution. This superior technical performance enabled the capture of high-quality mechanical mapping in vivo even while the subject was breathing normally. Importantly, we improved Brillouin spectral measurements to achieve relative accuracy better than 0.07% verified by rubidium absorption frequencies, with 0.12% stability over 2000 seconds. These specifications finally yield the Brillouin measurement sensitivity that is required to detect ophthalmology-relevant corneal biomechanical properties.

In situ assessment of lens elasticity with noncontact optical coherence elastography

Lens biomechanics has great potential for application in clinical diagnostics and treatment monitoring of presbyopia and cataracts. However, current approaches to lens elastography do not meet the desired safety or sensitivity for clinical application. In this regard, we propose a noncontact optical coherence elastography (OCE) method to facilitate quantitative in situ imaging of lens elasticity. Elastic waves induced by air-pulse stimulation on the limbus propagate to the lens and are then imaged using custom-built swept-source optical coherence tomography to obtain the elastic wave velocity and Young's modulus. The proposed OCE method was first validated by comparing the results of in situ and in vitro measurements of porcine lenses. The results demonstrate that the Young's modulus measured in situ was highly consistent with that measured in vitro and had an intraclass correlation coefficient of 0.988. We further investigated the elastic changes induced by cold storage and microwave heating. During 36-hour cold storage, the mean Young's modulus gradually increased (from 5.62 ± 1.24 kPa to 11.40 ± 2.68 kPa, P < 0.0001, n = 9) along with the formation of nuclear opacities. 15-second microwave heating caused a greater increase in the mean Young's modulus (from 6.86 ± 1.21 kPa to 25.96 ± 8.64 kPa, P < 0.0025, n = 6) without apparent cataract formation. Accordingly, this study reports the first air-pulse OCE measurements of in situ lenses, which quantified the loss of lens elasticity during simulated cataract development with good repeatability and sensitivity, thus enhancing the potential for adoption of lens biomechanics in the clinic.

Deep Tissue Characterization with Optical Coherence Elastography: A Comparison of Different Methods

The measurement of the biomechanical properties of the skin is of great interest since these properties play an important role in the development of several diseases such as skin cancer and systemic sclerosis. In this direction, several diagnostic tools have been developed to analyze the mechanical properties of the skin. Optical coherence elastography (OCE) is one of the emerging imaging techniques used for the characterization of the mechanical properties of the tissue quantitatively. In systemic sclerosis patients, the measurement of the mechanical properties of the deeper skin layers is desirable compared to the superficial layers. There are several variants of OCE that exist, but it is still not clear which method is more suitable for the measurement of the mechanical properties of the deeper tissue. In this work, we tested three common methods, the pulsed excitation method, the continuous wave excitation method, and the resonant frequency method, for the measurement of the mechanical properties of the deeper layers in the tissue. We found out that the pulsed wave excitation method provides the most reliable measurements in the shortest possible time compared to the other two methods.

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

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

Imaging of Single Transducer-Harmonic Motion Imaging-Derived Displacements at Several Oscillation Frequencies Simultaneously

Mapping of mechanical properties, dependent on the frequency of motion, is relevant in diagnosis, monitoring treatment response, or intra-operative surgical resection planning. While shear wave speeds at different frequencies have been described elsewhere, the effect of frequency on the "on-axis" acoustic radiation force (ARF)-induced displacement has not been previously investigated. Instead of generating single transducer-harmonic motion imaging (ST-HMI)-derived peak-to-peak displacement (P2PD) image at a particular frequency, a novel multi-frequency excitation pulse is proposed to generate P2PD images at 100-1000 Hz simultaneously. The performance of the proposed excitation pulse is compared with the ARFI by imaging 16 different inclusions (Young's moduli of 6, 9, 36, 70 kPa and diameters of 1.6, 2.5, 6.5, and 10.4 mm) embedded in an 18 kPa background. Depending on inclusion size and stiffness, the maximum CNR and contrast were achieved at different frequencies and were always higher than ARFI. The frequency, at which maximum CNR and contrast were achieved, increased with stiffness for fixed inclusion's size and decreased with size for fixed stiffness. In vivo feasibility is tested by imaging a 4T1 breast cancer mouse tumor on Day 6, 12, and 19 post-injection of tumor cells. Similar to phantoms, the CNR of ST-HMI images was higher than ARFI and increased with frequency for the tumor on Day 6. Besides, P2PD at 100-1000 Hz indicated that the tumor became stiffer with respect to the neighboring non-cancerous tissue over time. These results indicate the importance of using a multi-frequency excitation pulse to simultaneously generate displacement at multiple frequencies to better delineate inclusions or tumors.

Intravascular optical coherence elastography

Optical coherence elastography (OCE), a functional extension of optical coherence tomography (OCT), visualizes tissue strain to deduce the tissue's biomechanical properties. In this study, we demonstrate intravascular OCE using a 1.1 mm motorized catheter and a 1.6 MHz Fourier domain mode-locked OCT system. We induced an intraluminal pressure change by varying the infusion rate from the proximal end of the catheter. We analysed the pixel-matched phase change between two different frames to yield the radial strain. Imaging experiments were carried out in a phantom and in human coronary arteries in vitro. At an imaging speed of 3019 frames/s, we were able to capture the dynamic strain. Stiff inclusions in the phantom and calcification in atherosclerotic plaques are associated with low strain values and can be distinguished from the surrounding soft material, which exhibits elevated strain. For the first time, circumferential intravascular OCE images are provided side by side with conventional OCT images, simultaneously mapping both the tissue structure and stiffness.

Mouse embryo phenotyping with optical coherence tomography

With the explosion of gene editing tools in recent years, there has been a much greater demand for mouse embryo phenotyping, and traditional methods such as histology and histochemistry experienced a methodological renaissance as they became the principal tools for phenotyping. However, it is important to explore alternative phenotyping options to maximize time and resources and implement volumetric structural analysis for enhanced investigation of phenotypes. Cardiovascular phenotyping, in particular, is important to perform in vivo due to the dramatic structural and functional changes that occur in heart development over relatively short periods of time. Optical coherence tomography (OCT) is one of the most exciting advanced imaging techniques emerging within the field of developmental biology, and this review provides a summary of how it is currently being implemented in mouse embryo investigations and phenotyping. This review aims to provide an understanding of the approaches used in optical coherence tomography and how they can be applied in embryology and developmental biology, with the overall aim of bridging the gap between biology and technology.

Following the Beat: Imaging the Valveless Pumping Function in the Early Embryonic Heart

In vertebrates, the coordinated beat of the early heart tube drives cardiogenesis and supports embryonic growth. How the heart pumps at this valveless stage marks a fascinating problem that is of vital significance for understanding cardiac development and defects. The developing heart achieves its function at the same time as continuous and dramatic morphological changes, which in turn modify its pumping dynamics. The beauty of this muti-time-scale process also highlights its complexity that requires interdisciplinary approaches to study. High-resolution optical imaging, particularly fast, four-dimensional (4D) imaging, plays a critical role in revealing the process of pumping, instructing numerical modeling, and enabling biomechanical analyses. In this review, we aim to connect the investigation of valveless pumping mechanisms with the recent advancements in embryonic cardiodynamic imaging, facilitating interactions between these two areas of study, in hopes of encouraging and motivating innovative work to further understand the early heartbeat.

Longitudinal assessment of the effect of alkali burns on corneal biomechanical properties using optical coherence elastography

Eye injury due to alkali burn is a severe ocular trauma that can profoundly affect corneal structure and function, including its biomechanical properties. Here, we assess the changes in the mechanical behavior of mouse corneas in response to alkali-induced injury by conducting longitudinal measurements using optical coherence elastography (OCE). A non-contact air-coupled ultrasound transducer was used to induce elastic waves in control and alkali-injured mouse corneas in vivo, which were imaged with phase-sensitive optical coherence tomography. Corneal mechanical properties were estimated using a modified Rayleigh-Lamb wave model, and results show that Young's modulus of alkali-burned corneas were significantly greater than that of their healthy counterparts on days 7 (p = 0.029) and 14 (p = 0.026) after injury. These findings, together with the changes in the shear viscosity coefficient postburn, indicate that the mechanical properties of the alkali-burned cornea are significantly modulated during the wound healing process.

In Vivo Determination of the Human Corneal Elastic Modulus Using Vibrational Optical Coherence Tomography

Purpose

To determine the in vivo elastic modulus of the human cornea using vibrational optical coherence tomography (VOCT).

Methods

Vibrational analysis coupled with optical coherence tomography (OCT) was used to obtain the resonant frequency (RF) and elastic modulus of corneal structural components. VOCT corneal thickness values were measured using OCT images and correlated with corneal thickness determined with Pentacam (Oculus, Wetzlar, Germany). Moduli were obtained at two locations: central cornea (CC) and inferior cornea (IC). Measurements were obtained with and without anesthetic eye drops to assess their effect on the modulus measurements.

Results

VOCT thickness values correlated positively (R² = 0.97) and linearly (y = 1.039x–16.89) with those of Pentacam. Five RF peaks (1–5) were present, although their presence was variable across eyes. The RF for peaks 1 to 5 in the CC and IC ranged from 73.5 ± 4.9 to 239 ± 3 Hz and 72.1 ± 6.3 to 238 ± 4 Hz, respectively. CC and IC moduli for peaks 1 to 5 ranged from 1.023 ± 0.104 to 6.87 ± 0.33 MPa and 0.98 ± 0.15 to 6.52 ± 0.79 MPa, respectively. Topical anesthesia did not significantly alter the modulus (P > 0.05 for all), except for peak 2 in the CC (P < 0.05).

Conclusions

This pilot study demonstrates the utility of VOCT as an in vivo, noninvasive technology to measure the elastic modulus in human corneas. The structural origin of these moduli is hypothesized based on previous reports, and further analyses are necessary for confirmation.

Translational Relevance

This work presents VOCT as a novel approach to assess the in vivo elastic modulus of the cornea, an indicator of corneal structural integrity and health.

Measuring Human Corneal Stromal Biomechanical Properties Using Tensile Testing Combined With Optical Coherence Tomography

Purpose: To investigate the ex vivo elastic modulus of human corneal stroma using tensile testing with optical coherence tomography (OCT) imaging and its correlation with in vivo measurements using corneal visualization Scheimpflug technology.

Methods: Twenty-four corneal specimens extracted from stromal lenticules through small incision lenticule extraction were cut into strips for uniaxial tensile tests. In vivo corneal biomechanical responses were evaluated preoperatively using the corneal visualization Scheimpflug technology (CorVis ST). The correlation of the elastic modulus with clinical characteristics and dynamic corneal response parameters were analyzed using Spearman’s correlation analysis.

Results: The mean low strain tangent modulus (LSTM) of the human corneal stroma was 0.204 ± 0.189 (range 0.010–0.641) MPa, and high strain tangent modulus (HSTM) 5.114 ± 1.958 (range 2.755–9.976) MPa. Both LSTM (r = 0.447, p = 0.029) and HSTM (r = 0.557, p = 0.005) were positively correlated with the stress-strain index (SSI). LSTM was also positively correlated with the A1 deflection length (r = 0.427, p = 0.037) and A1 deflection area (r = 0.441, p = 0.031). HSTM was positively correlated with spherical equivalent (r = 0.425, p = 0.038).

Conclusions: The correlation of corneal elastic modulus with A1 deflection parameters and SSI may indicate a relationship between these parameters and tissue elasticity. The HSTM decreased with the degree of myopia. Combining tensile test with OCT may be a promising approach to assess corneal biomechanical properties.

Multi-class classification of breast tissue using optical coherence tomography and attenuation imaging combined via deep learning

We demonstrate a convolutional neural network (CNN) for multi-class breast tissue classification as adipose tissue, benign dense tissue, or malignant tissue, using multi-channel optical coherence tomography (OCT) and attenuation images, and a novel Matthews correlation coefficient (MCC)-based loss function that correlates more strongly with performance metrics than the commonly used cross-entropy loss. We hypothesized that using multi-channel images would increase tumor detection performance compared to using OCT alone. 5,804 images from 29 patients were used to fine-tune a pre-trained ResNet-18 network. Adding attenuation images to OCT images yields statistically significant improvements in several performance metrics, including benign dense tissue sensitivity (68.0% versus 59.6%), malignant tissue positive predictive value (PPV) (79.4% versus 75.5%), and total accuracy (85.4% versus 83.3%), indicating that the additional contrast from attenuation imaging is most beneficial for distinguishing between benign dense tissue and malignant tissue.

Analysis of strain estimation methods in phase-sensitive compression optical coherence elastography

In compression optical coherence elastography (OCE), deformation is quantified as the local strain at each pixel in the OCT field-of-view. A range of strain estimation methods have been demonstrated, yet it is unclear which method provides the best performance. Here, we analyze the two most prevalent strain estimation methods used in phase-sensitive compression OCE, i.e., weighted least squares (WLS) and the vector method. We introduce a framework to compare strain imaging metrics, incorporating strain sensitivity, strain signal-to-noise ratio (SNR), strain resolution, and strain accuracy. In addition, we propose a new phase unwrapping algorithm in OCE, fast phase unwrapping (FPU), and combine it with WLS, termed WLS_FPU. Using the framework, we compare this new strain estimation method with both a current implementation of WLS that incorporates weighted phase unwrapping (WPU), termed WLS_WPU, and the vector method. Our analysis reveals that the three methods provide similar strain sensitivity, strain SNR, and strain resolution, but that WLS_FPU extends the dynamic range of accurate, measurable local strain, e.g., measuring a strain of 2.5 mɛ with ∼4% error, that is ×11 and ×15 smaller than the error measured using WLS_WPU and the vector method, respectively. We also demonstrate, for the first time, the capability to detect sub-resolution contrast in compression OCE, i.e., changes in strain occurring within the strain axial resolution, and how this contrast varies between the different strain estimation methods. Lastly, we compare the performance of the three strain estimation methods on mouse skeletal muscle and human breast tissue and demonstrate that WLS_FPU avoids strain imaging artifacts resulting from phase unwrapping errors in WLS_WPU and provides improved contrast over the other two methods.

Single-cell all-optical coherence elastography with optical tweezers

The elastic properties of cells are important for many of their functions, however the development of label free noninvasive cellular elastography method is a challenging topic. We present a novel single-cell all-optical coherence elastography method that combines optical tweezers producing mechanical excitation on the cell membrane or organelle and phase-sensitive optical coherence microscopy measuring sample response and determining its mechanical properties. The method allows living cells imaging with a lateral resolution of 0.5 μm and an axial resolution up to 10 nm, making it possible to detect nanometer displacements of the cell organelles and to record the propagation of mechanical wave along the cell membrane in response to optical tweezers excitation. We also demonstrate applicability of the method on single living red blood cells, yeast and cancer cells. The all-optical nature of the method developed makes it a promising and easily applicable tool for studying cellular and subcellular mechanics in vivo.

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.

Mechanical modeling and characterization of human skin: A review

This paper reviews the advances made in recent years on modeling approaches and experimental techniques to characterize the mechanical properties of human skin. The skin is the largest organ of the human body that has a complex multi-layered structure with different mechanical behaviors. The mechanical properties of human skin play an important role in distinguishing between healthy and unhealthy skin. Furthermore, knowing these mechanical properties enables computer simulation, skin research, clinical studies, as well as diagnosis and treatment monitoring of skin diseases. This paper reviews the recent efforts on modeling skin using linear, nonlinear, viscoelastic, and anisotropic materials. The work also focuses on aging effects, microstructure analysis, and non-invasive methods for skin testing. A detailed explanation of the skin structure and numerical models, such as finite element models, are discussed in this work. This work also compares different experimental methods that measure the mechanical properties of human skin. The work reviews the experimental results in the literature and shows how the mechanical properties of human skin vary with the skin sites, the layers, and the structure of human skin. The paper also discusses how state-of-the-art technology can advance skin research.

Real-Time Strain and Elasticity Imaging in Phase-Sensitive Optical Coherence Elastography Using a Computationally Efficient Realization of the Vector Method

We present a real-time realization of OCT-based elastographic mapping local strains and distribution of the Young’s modulus in biological tissues, which is in high demand for biomedical usage. The described variant exploits the principle of Compression Optical Coherence Elastography (C-OCE) and uses processing of phase-sensitive OCT signals. The strain is estimated by finding local axial gradients of interframe phase variations. Instead of the popular least-squares method for finding these gradients, we use the vector approach, one of its advantages being increased computational efficiency. Here, we present a modified, especially fast variant of this approach. In contrast to conventional correlation-based methods and previously used phase-resolved methods, the described method does not use any search operations or local calculations over a sliding window. Rather, it obtains local strain maps (and then elasticity maps) using several transformations represented as matrix operations applied to entire complex-valued OCT scans. We first elucidate the difference of the proposed method from the previously used correlational and phase-resolved methods and then describe the proposed method realization in a medical OCT device, in which for real-time processing, a “typical” central processor (e.g., Intel Core i7-8850H) is sufficient. Representative examples of on-flight obtained elastographic images are given. These results open prospects for broad use of affordable OCT devices for high-resolution elastographic vitalization in numerous biomedical applications, including the use in clinic.

Characterization of the Biomechanical Properties of Skin Using Vibrational Optical Coherence Tomography: Do Changes in the Biomechanical Properties of Skin Stroma Reflect Structural Changes in the Extracellular Matrix of Cancerous Lesions?

Early detection of skin cancer is of critical importance since the five-year survival rate for early detected skin malignancies is 99% but drops to 27% for cancer that has spread to distant lymph nodes and other organs. Over 2.5 million benign skin biopsies (55% of the total) are performed each year in the US at an alarming cost of USD ~2.5 B. Therefore there is an unmet need for novel non-invasive diagnostic approaches to better differentiate between cancerous and non-cancerous lesions, especially in cases when there is a legitimate doubt that a biopsy may be required. The purpose of this study is to determine whether the differences in the extracellular matrices among normal skin, actinic keratosis (AK), basal cell carcinoma (BCC) and squamous cell carcinoma (SCC) can be assessed non-invasively using vibrational optical coherence tomography (VOCT). VOCT is a new diagnostic technology that uses infrared light and audible sound applied transversely to tissue to measure the resonant frequencies and elastic moduli of cells, dermal collagen, blood vessels and fibrous tissue in skin and lesion stroma without physically touching the skin. Our results indicate that the cellular, vascular and fibrotic resonant frequency peaks are altered in AK, BCC and SCC compared to those peaks observed in normal skin and can serve as physical biomarkers defining the differences between benign and cancerous skin lesions. The resonant frequency is increased from a value of 50 Hz in normal skin to a value of about 80 Hz in pre- and cancerous lesions. A new vascular peak is seen at 130 Hz in cancerous lesions that may reflect the formation of new tumor blood vessels. The peak at 260 Hz is similar to that seen in the skin of a subject with Scleroderma and skin wounds that have healed. The peak at 260 Hz appears to be associated with the deposition of large amounts of stiff fibrous collagen in the stroma surrounding cancerous lesions. Based on the results of this pilot study, VOCT can be used to non-invasively identify physical biomarkers that can help differentiate between benign and cancerous skin lesions. The appearance of new stiff cellular, fragile new vessels, and stiff fibrous material based on resonant frequency peaks and changes in the extracellular matrix can be used as a fingerprint of pre- and cancerous skin lesions.

Method for the biomechanical analysis of aqueous veins and perilimbal sclera by three-dimensional photoacoustic imaging and strain field calculation

A method motivated by the eye's aqueous veins is described for the imaging and strain calculation within soft biological tissues. A challenge to the investigation of the biomechanics of the aqueous vein-perilimbal sclera tissue complex is resolution of tissue deformations as a function of intraocular pressure and the subsequent calculation of strain (a normalized measure of deformation). The method involves perfusion of the eye with a contrast agent during conduction of non-invasive, optical resolution photoacoustic microscopy. This imaging technique permits three-dimensional displacement measurements of tracked points on the inner walls of the veins which are used in a finite element model to determine the corresponding strains. The methods are validated against two standard strain measurement methods. Representative porcine globe perfusion experiments are presented that demonstrate the power of the method to determine complex strain fields in the veins dependent on intraocular pressure as well as vein anatomy. In these cases, veins are observed to move radially outward during increases in intraocular pressure and to possess significant spatial strain variation, possibly influenced by their branching patterns. To the authors' knowledge, these are the only such quantitative, data driven, calculations of the aqueous vein strains available in the open literature.

Magnetic Resonance Elastography Reconstruction for Anisotropic Tissues

Elastography has become widely used clinically for characterising changes in soft tissue mechanics that are associated with altered tissue structure and composition. However, some soft tissues, such as muscle, are not isotropic as is assumed in clinical elastography implementations. This limits the ability of these methods to capture changes in anisotropic tissues associated with disease. The objective of this study was to develop and validate a novel elastography reconstruction technique suitable for estimating the linear viscoelastic mechanical properties of transversely isotropic soft tissues. We derived a divergence-free formulation of the governing equations for acoustic wave propagation through a linearly transversely isotropic viscoelastic material, and transformed this into a weak form. This was then implemented into a finite element framework, enabling the analysis of wave input data and tissue structural fibre orientations, in this case based on diffusion tensor imaging. To validate the material constants obtained with this method, numerous in silico phantom experiments were run which encompassed a range of variations in wave input directions, material properties, fibre structure and noise. The method was also tested on ex vivo muscle and in vivo human volunteer calf muscles, and compared with a previous curl-based inversion method. The new method robustly extracted the transversely isotropic shear moduli (G_⊥^', G_∥^', G^″) from the in silico phantom tests with minimal bias, including in the presence of experimentally realistic levels of noise in either fibre orientation or wave data. This new method performed better than the previous method in the presence of noise. Anisotropy estimates from the ex vivo muscle phantom agreed well with rheological tests. In vivo experiments on human calf muscles were able to detect increases in muscle shear moduli with passive muscle stretch. This new reconstruction method can be applied to quantify tissue mechanical properties of anisotropic soft tissues, such as muscle, in health and disease.

Polarisation optics for biomedical and clinical applications: a review

Many polarisation techniques have been harnessed for decades in biological and clinical research, each based upon measurement of the vectorial properties of light or the vectorial transformations imposed on light by objects. Various advanced vector measurement/sensing techniques, physical interpretation methods, and approaches to analyse biomedically relevant information have been developed and harnessed. In this review, we focus mainly on summarising methodologies and applications related to tissue polarimetry, with an emphasis on the adoption of the Stokes-Mueller formalism. Several recent breakthroughs, development trends, and potential multimodal uses in conjunction with other techniques are also presented. The primary goal of the review is to give the reader a general overview in the use of vectorial information that can be obtained by polarisation optics for applications in biomedical and clinical research.

Application of Magnetic Resonance Imaging in Liver Biomechanics: A Systematic Review

MRI-based biomechanical studies can provide a deep understanding of the mechanisms governing liver function, its mechanical performance but also liver diseases. In addition, comprehensive modeling of the liver can help improve liver disease treatment. Furthermore, such studies demonstrate the beginning of an engineering-level approach to how the liver disease affects material properties and liver function. Aimed at researchers in the field of MRI-based liver simulation, research articles pertinent to MRI-based liver modeling were identified, reviewed, and summarized systematically. Various MRI applications for liver biomechanics are highlighted, and the limitations of different viscoelastic models used in magnetic resonance elastography are addressed. The clinical application of the simulations and the diseases studied are also discussed. Based on the developed questionnaire, the papers' quality was assessed, and of the 46 reviewed papers, 32 papers were determined to be of high-quality. Due to the lack of the suitable material models for different liver diseases studied by magnetic resonance elastography, researchers may consider the effect of liver diseases on constitutive models. In the future, research groups may incorporate various aspects of machine learning (ML) into constitutive models and MRI data extraction to further refine the study methodology. Moreover, researchers should strive for further reproducibility and rigorous model validation and verification.

Lymph vessels visualization from optical coherence tomography data using depth-resolved attenuation coefficient calculation

Multimodal optical coherent tomography grows popularity with researchers and clinicians over the past decade. One of the modalities is lymphangiography, which allows visualization of the lymphatic vessel networks within optical coherence tomography (OCT) imaging volume. In the present study, it is shown that lymphatic vessel visualization obtained from the depth-resolved attenuation coefficient distributions, corrected for the noise, shows improved contrast and detail in comparison with previously proposed approaches. We also argue that the two most popular approaches for lymphatic vessel visualization, namely simple intensity thresholding and vesselness calculation based on local Hessian matrix eigenvalues, imply different definitions of the lymphatic vessel's appearance in the OCT volume and lead to the different networks.

Sensing Inhalation Injury-Associated Changes in Airway Wall Compliance by Anatomic Optical Coherence Elastography

Quantitative methods for assessing the severity of inhalation (burn) injury are needed to aid in treatment decisions. We hypothesize that it is possible to assess the severity of injuries on the basis of differences in the compliance of the airway wall. Here, we demonstrate the use of a custom-built, endoscopic, anatomic optical coherence elastography (aOCE) system to measure airway wall compliance. The method was first validated using airway phantoms, then performed on ex vivo porcine tracheas under varying degrees of inhalation (steam) injury. A negative correlation between aOCE-derived compliance and severity of steam injuries is found, and spatially-resolved compliance maps reveal regional heterogeneity in airway properties.

Measuring the elastic modulus of soft culture surfaces and three-dimensional hydrogels using atomic force microscopy

Growing interest in exploring mechanically mediated biological phenomena has resulted in cell culture substrates and 3D matrices with variable stiffnesses becoming standard tools in biology labs. However, correlating stiffness with biological outcomes and comparing results between research groups is hampered by variability in the methods used to determine Young's (elastic) modulus, E, and by the inaccessibility of relevant mechanical engineering protocols to most biology labs. Here, we describe a protocol for measuring E of soft 2D surfaces and 3D hydrogels using atomic force microscopy (AFM) force spectroscopy. We provide instructions for preparing hydrogels with and without encapsulated live cells, and provide a method for mounting samples within the AFM. We also provide details on how to calibrate the instrument, and give step-by-step instructions for collecting force-displacement curves in both manual and automatic modes (stiffness mapping). We then provide details on how to apply either the Hertz or the Oliver-Pharr model to calculate E, and give additional instructions to aid the user in plotting data distributions and carrying out statistical analyses. We also provide instructions for inferring differential matrix remodeling activity in hydrogels containing encapsulated single cells or organoids. Our protocol is suitable for probing a range of synthetic and naturally derived polymeric hydrogels such as polyethylene glycol, polyacrylamide, hyaluronic acid, collagen, or Matrigel. Although sample preparation timings will vary, a user with introductory training to AFM will be able to use this protocol to characterize the mechanical properties of two to six soft surfaces or 3D hydrogels in a single day.

Feasibility of Harmonic Motion Imaging Using a Single Transducer: In Vivo Imaging of Breast Cancer in a Mouse Model and Human Subjects

Harmonic motion imaging (HMI) interrogates the mechanical properties of tissues by simultaneously generating and tracking harmonic oscillation using focused ultrasound and imaging transducers, respectively. Instead of using two transducers, the objective of this work is to develop a single transducer HMI (ST-HMI) to both generate and track harmonic motion at "on-axis" to the force for facilitating data acquisition. In ST-HMI, the amplitude-modulated force was generated by modulating excitation pulse duration and tracking of motion was performed by transmitting tracking pulses interleaved between excitation pulses. The feasibility of ST-HMI was performed by imaging two elastic phantoms with three inclusions (N = 6) and comparing it with acoustic radiation force impulse (ARFI) imaging, in vivo longitudinal monitoring of 4T1, orthotropic breast cancer mice (N = 4), and patients (N = 3) with breast masses in vivo. Six inclusions with Young's moduli of 8, 10, 15, 20, 40, and 60 kPa were embedded in a 5 kPa background. The ST-HMI-derived peak-to-peak displacement (P2PD) successfully detected all inclusions with [Formula: see text] of the linear regression between the P2PD ratio of background to inclusion versus Young's moduli ratio of inclusion to background. The contrasts of 10 and 15 kPa inclusions were higher in ST-HMI than ARFI-derived images. In the mouse study, the median P2PD ratio of tumor to non-cancerous tissues was 3.0, 5.1, 6.1, and 7.7 at 1, 2, 3, and 4 weeks post-injection of the tumor cells, respectively. In the clinical study, ST-HMI detected breast masses including fibroadenoma, pseudo angiomatous stromal hyperplasia, and invasive ductal carcinoma with a P2PD ratio of 1.37, 1.61, and 1.78, respectively. These results indicate that ST-HMI can assess the mechanical properties of tissues via generation and tracking of harmonic motion "on-axis" to the ARF. This study is the first step towards translating ST-HMI in clinics.

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.

Optical coherence viscometry

We report a technique, named optical coherence viscometry (OCV), to measure the viscosity of Newtonian fluids in a noncontact manner. According to linear wave theory with small amplitudes, capillary waves are associated with fluid mechanical properties. To perform this measurement and avoid the overdamped effects of capillary waves in viscous fluids, transient acoustic radiation force was applied to generate capillary waves. Within a very limited field-of-view using optical coherence tomography, wave motion acquired in the time domain was analyzed using Fourier methods to study the wave velocity dispersion and attenuation relationships for capillary waves, which can reduce the fluid quantity drastically into tissue culture scale. We measure the viscosities of water, water-glycerol solutions with three concentrations, and biological plasma using the proposed OCV and compare the experimental results to theoretical calculations. OCV is sensitive to wave perturbations and can be a promising technique for measuring the viscosity of biological fluids and could be applied in future applications for measurements for lipid membranes in cell biology and tissue engineering investigation.

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.

Characterizing thrombus with multiple red blood cell compositions by optical coherence tomography attenuation coefficient

Embolectomy is one of the emergency procedures performed to remove emboli. Assessing the composition of human blood clots is an important diagnostic factor and could provide guidance for an appropriate treatment strategy for interventional physicians. Immunostaining has been used to identity compositions of clots as a gold-standard procedure, but it is time-consuming and cannot be performed in situ. Here, we proposed that the optical attenuation coefficient of optical coherence tomography (OCT) can be a reliable indicator as a new imaging modality to differentiate clot compositions. Fifteen human blood clots with multiple red blood cell (RBC) compositions from 21% to 95% were prepared using healthy human whole blood. A homogeneous gelatin phantom experiment and numerical simulation based on the Lambert-Beer's law were examined to verify the validity of the attenuation coefficient estimation. The results displayed that optical attenuation coefficients were strongly correlated with RBC compositions. We reported that attenuation coefficients could be a promising biomarker to guide the choice of an appropriate interventional device in a clinical setting and assist in characterizing blood clots.

Computational 4D-OCM for label-free imaging of collective cell invasion and force-mediated deformations in collagen

Traction force microscopy (TFM) is an important family of techniques used to measure and study the role of cellular traction forces (CTFs) associated with many biological processes. However, current standard TFM methods rely on imaging techniques that do not provide the experimental capabilities necessary to study CTFs within 3D collective and dynamic systems embedded within optically scattering media. Traction force optical coherence microscopy (TF-OCM) was developed to address these needs, but has only been demonstrated for the study of isolated cells embedded within optically clear media. Here, we present computational 4D-OCM methods that enable the study of dynamic invasion behavior of large tumor spheroids embedded in collagen. Our multi-day, time-lapse imaging data provided detailed visualizations of evolving spheroid morphology, collagen degradation, and collagen deformation, all using label-free scattering contrast. These capabilities, which provided insights into how stromal cells affect cancer progression, significantly expand access to critical data about biophysical interactions of cells with their environment, and lay the foundation for future efforts toward volumetric, time-lapse reconstructions of collective CTFs with TF-OCM.

Strain and elasticity imaging in compression optical coherence elastography: The two-decade perspective and recent advances

Quantitative mapping of deformation and elasticity in optical coherence tomography has attracted much attention of researchers during the last two decades. However, despite intense effort it took ~15 years to demonstrate optical coherence elastography (OCE) as a practically useful technique. Similarly to medical ultrasound, where elastography was first realized using the quasi-static compression principle and later shear-wave-based systems were developed, in OCE these two approaches also developed in parallel. However, although the compression OCE (C-OCE) was proposed historically earlier in the seminal paper by J. Schmitt in 1998, breakthroughs in quantitative mapping of genuine local strains and the Young's modulus in C-OCE have been reported only recently and have not yet obtained sufficient attention in reviews. In this overview, we focus on underlying principles of C-OCE; discuss various practical challenges in its realization and present examples of biomedical applications of C-OCE. The figure demonstrates OCE-visualization of complex transient strains in a corneal sample heated by an infrared laser beam.