Ultrasound Shear Wave Viscoelastography: Model-Independent Quantification of the Complex Shear Modulus

Ultrasound shear wave viscoelastography to characterize liver nodules

Purpose. To investigate the diagnostic performance of ultrasound (US)-based shear wave speed (SWS), shear wave attenuation (SWA), and combination of them as shear wave viscoelastography (SWVE) methods in patients undergoing US to characterize focal liver nodules.Materials and methods. In this prospective cross-sectional study, 70 patients with 72 nodules were enrolled. Investigational US and clinical magnetic resonance imaging (MRI) examinations were performed in all participants. The composite reference standard included MRI or histopathology to differentiate benign and malignant nodules. A linear discriminant analysis (LDA) was used to assess the combination of SWVE methods. Analyzes included Mann-WhitneyUtest, receiver operating characteristic analysis, and computation of sensitivity and specificity at the point that maximized the Youden index.Results. Mean SWS was significantly higher in malignant than benign nodules (2.49 ± 0.76 m s-1vs. 1.72 ± 0.70,p< 0.001), whereas SWA was lower (0.56 ± 0.30 vs. 1.10 ± 0.43 Np/m/Hz,p< 0.001). To differentiate between malignant and benign nodules, SWS with a threshold of 2.43 m s-1achieved a sensitivity of 0.54 (95% confidence interval (CI): 0.38-0.69) and a specificity of 0.88 (CI: 0.74-0.95). SWA with a threshold of 0.81 Np/m/Hz yielded a sensitivity of 0.81 (CI: 0.66-0.90) and a specificity of 0.74 (CI: 0.58-0.86). Combining these SWVE methods using a LDA resulted in a sensitivity of 0.81 (CI: 0.66-0.91) and a specificity of 0.86 (CI: 0.71-0.94).Conclusion. Malignant nodules had higher SWS and lower SWA than benign ones. The combination of SWS and SWA in a LDA classification algorithm increased the diagnostic performance.

Estimation of shear viscoelasticity via time-domain elastic full waveform inversion in ultrasound shear wave elastography

The elastography inversion process typically assumes local homogeneity or ignores viscosity, which can lead to inaccuracies. This study presents a methodology for estimating viscous and elastic properties of a heterogenous media based on time-domain elastic full waveform inversion using the limited data available in conventional ultrasound shear wave elastography (SWE). The proposed method is based on the elastodynamic equation and minimizes a squared residual misfit function. To improve the optimization process, a total viscoelastic search space is considered and a quasi-Newton optimization is employed, where gradients are approximated using the adjoint state method. Additionally, Tikhonov regularization is incorporated to handle noisy and sparse displacement data. The proposed method was evaluated across various scenarios simulating real-world experimental conditions, accounting for noise levels and temporal sampling sparsity of displacement field. Using noisy displacement data from a single slice of three-dimensional volume, imitating the ultrasound SWE, resulted in 2.91% error for elasticity estimation and 23.97% error for viscosity estimation across the field of view. This framework enables the estimation of viscous and elastic properties of heterogeneous media and has shown promising results. Additionally, the optimization was improved, and noisy and sparse displacement data were effectively addressed.

Shear wave elastography primer for the abdominal radiologist

PURPOSE

Shear wave elastography (SWE) provides a means for adding information about the mechanical properties of tissues to a diagnostic ultrasound examination. It is important to understand the physics and methods by which the measurements are made to aid interpretation of the results as they relate to disease processes.

METHODS

The components of how ultrasound is used to generate shear waves and make measurements of the induced motion are reviewed. The physics of shear wave propagation are briefly described for elastic and viscoelastic tissues. Additionally, shear wave propagation in homogeneous and inhomogeneous cases is addressed.

RESULTS

SWE technology has been implemented by many clinical vendors with different capabilities. Various quality metrics are used to define valid measurements based on aspects of the shear wave signals or wave velocity estimates.

CONCLUSION

There are many uses for SWE in abdominal imaging, but it is important to understand how the measurements are performed to gauge their utility for diagnosis of different conditions. Continued efforts to make the technology robust in complex clinical situations are ongoing, but many applications actively benefit from added information about tissue mechanical properties for a more holistic view of the patient for diagnosis or assessment of prognosis and treatment management.

Between-Visit Reproducibility of Shear Wave Viscoelastography in Volunteers and Patients With Metabolic Dysfunction-Associated Steatotic Liver Disease

OBJECTIVE

To assess the reproducibility of six ultrasound (US)-determined shear wave (SW) viscoelastography parameters for assessment of mechanical properties of the liver in volunteers and patients with biopsy-proven metabolic dysfunction-associated steatotic liver disease (MASLD) or metabolic dysfunction-associated steatohepatitis (MASH).

METHODS

This prospective, cross-sectional, institutional review board-approved study included 10 volunteers and 20 patients with MASLD or MASH who underwent liver US elastography twice, at least 2 weeks apart. SW speed (SWS), Young's modulus (E), shear modulus (G), SW attenuation (SWA), SW dispersion (SWD), and viscosity were computed from radiofrequency data recorded on a research US scanner. Linear mixed models were used to consider the sonographer on duty as a confounder. The reproducibility of measurements was assessed by intraclass correlation coefficient (ICC), coefficient of variation (CV), reproducibility coefficient (RDC), and Bland-Altman analyses.

RESULTS

The sonographer performing the exam had no impact on viscoelastic parameters (P > .05). ICCs of SWS, E, G, SWA, SWD, and viscosity were, respectively, 0.89 (95% confidence intervals [CI]: 0.79-0.95), 0.81 (95% CI: 0.79-0.95), 0.90 (95% CI: 0.80-0.95), 0.96 (95% CI: 0.93-0.98), 0.78 (95% CI: 0.60-0.89), and 0.90 (95% CI: 0.80-0.95); CVs were 11.9, 23.3, 24.2, 10.1, 29.0, and 32.2%; RDCs were 33.0, 64.5, 66.9, 27.7, 80.3, and 89.2%, and Bland-Altman mean biases and 95% limits of agreement were -0.05 (-0.45, 0.35) m/s, -0.61 (-5.33, 4.10) kPa, -0.25 (-2.06, 1.56) kPa, -0.01 (-0.27, 0.26) Np/m/Hz, -0.09 (-7.09, 6.91) m/s/kHz, and -0.33 (-2.60, 1.94) Pa/s, between the two visits.

CONCLUSION

US-determined viscoelastography parameters can be measured with high reproducibility and consistency between two visits 2 weeks apart on the same ultrasound machine.

A Comprehensive Review on the Viscoelastic Parameters Used for Engineering Materials, Including Soft Materials, and the Relationships between Different Damping Parameters

Although the physical properties of a structure, such as stiffness, can be determined using some statical tests, the identification of damping parameters requires a dynamic test. In general, both theoretical prediction and experimental identification of damping are quite difficult. There are many different techniques available for damping identification, and each method gives a different damping parameter. The dynamic indentation method, rheometry, atomic force microscopy, and resonant vibration tests are commonly used to identify the damping of materials, including soft materials. While the viscous damping ratio, loss factor, complex modulus, and viscosity are quite common to describe the damping of materials, there are also other parameters, such as the specific damping capacity, loss angle, half-power bandwidth, and logarithmic decrement, to describe the damping of various materials. Often, one of these parameters is measured, and the measured parameter needs to be converted into another damping parameter for comparison purposes. In this review, the theoretical derivations of different parameters for the description and quantification of damping and their relationships are presented. The expressions for both high damping and low damping are included and evaluated. This study is considered as the first comprehensive review article presenting the theoretical derivations of a large number of damping parameters and the relationships among many damping parameters, with a quantitative evaluation of accurate and approximate formulas. This paper could be a primary resource for damping research and teaching.

Characterizing dispersion in bovine liver using ARFI-based shear wave rheometry

Background:Dispersion presents both a challenge and a diagnostic opportunity in shear wave elastography (SWE).Shear Wave Rheometry(SWR) is an inversion technique for processing SWE data acquired using an acoustic radiation force impulse (ARFI) excitation. The main advantage of SWR is that it can characterize the shear properties of homogeneous soft media over a wide frequency range. Assumptions associated with SWR include tissue homogeneity, tissue isotropy, and axisymmetry of the ARFI excitation).Objective:Evaluate the validity of the SWR assumptions in ex vivo bovine liver.Approach:SWR was used to measure the shear properties of bovine liver tissue as function of frequency over a large frequency range. Assumptions associated with SWR (tissue homogeneity, tissue isotropy, and axisymmetry of the ARFI excitation) were evaluated through measurements performed at multiple locations and probe orientations. Measurements focused on quantities that would reveal violations of the assumptions.Main results:Measurements of shear properties were obtained over the 25-250 Hz range, and showed a 4-fold increase in shear storage modulus (from 1 to 4 kPa) and over a 10-fold increase in the loss modulus (from 0.2 to 3 kPa) over that decade-wide frequency range. Measurements under different conditions were highly repeatable, and model error was low in all cases.Significance and Conclusion:SWR depends on modeling the ARFI-induced shear wave as a full vector viscoelastic shear wave resulting from an axisymmetric source; it is agnostic to any specific rheological model. Despite this generality, the model makes three main simplifying assumptions. These results show that the modeling assumptions used in SWR are valid in bovine liver over a wide frequency band.

Ultrasound Viscoelastography by Acoustic Radiation Force: A State-of-the-Art Review

Ultrasound elastography (USE) is a promising tool for tissue characterization as several diseases result in alterations of tissue structure and composition, which manifest as changes in tissue mechanical properties. By imaging the tissue response to an applied mechanical excitation, USE mimics the manual palpation performed by clinicians to sense the tissue elasticity for diagnostic purposes. Next to elasticity, viscosity has recently been investigated as an additional, relevant, diagnostic biomarker. Moreover, since biological tissues are inherently viscoelastic, accounting for viscosity in the tissue characterization process enhances the accuracy of the elasticity estimation. Recently, methods exploiting different acquisition and processing techniques have been proposed to perform ultrasound viscoelastography. After introducing the physics describing viscoelasticity, a comprehensive overview of the currently available USE acquisition techniques is provided, followed by a structured review of the existing viscoelasticity estimators classified according to the employed processing technique. These estimators are further reviewed from a clinical usage perspective, and current outstanding challenges are discussed.

Probing Tissue Viscoelasticity With STL Ultrasound Shearwave Spectroscopy Using Cole-Cole Diagrams

OBJECTIVE

Viscoelasticity is mapped by dispersion in shearwave elastography. Incomplete spectral information of shearwaves is therefore used to estimate mechanical stiffness. We propose capturing the "full-waveform-information" of the shear wave spectra to better resolve complex shear modulus μ* (ω). Approach is validated on phantom models, animal tissues, and feasibility demonstrated on human post-delivery placenta.

METHODS

We captured robust estimates of μ* in ex-vivo livers subjected to water bath ablation, glutaraldehyde exposure and in the placenta.

RESULTS

Complex modulus at 200 Hz is more reflective of tissue stiffness than cross-correlation estimate. Bias increased in phantoms with higher gelatin (G) (0.65: 6% G) and oil (O) (0.58: 6% G and 40% O) concentration, compared to elastic phantoms with low stiffness (0.33: 3% G). Actual tissues also reported higher bias in cross-correlation estimate (rabbit liver: 0.61, porcine liver: 2.20, and human placenta: 0.63). Stiffness is sensitive to ablation temperature, where the overall modulus changed from 3.02 KPa at 16 °C to 2.75 KPa at 56 °C in water bath. With exposure to Glutaraldehyde, the overall modulus increased from 2.37 to 9.03 KPa. Reconstruction errors in the loss modulus decreased by 68% with the power law compared to a Maxwell model in porcine livers with Cole-Cole inverse fitting.

CONCLUSION

Omitting Shear wave attenuation leads to bias. Reconstruction of rheological response with a model is sensitive to its architecture and also the framework.

SIGNIFICANCE

We use "full spectral information" in ultrasound shear wave elastography to better map μ*(ω) changes in viscoelastic tissues.

Systematic quantification of differences in shear wave elastography estimates between linear-elastic and viscoelastic material assumptionsa)

Quantitative, accurate, and standardized metrics are important for reliable shear wave elastography (SWE)-based biomarkers. For over two decades, the linear-elastic material assumption has been employed in SWE modes. In recent years, viscoelasticity estimation methods have been adopted in a few clinical systems. The current study aims to systematically quantify differences in SWE estimates obtained using linear-elastic and viscoelastic material assumptions. An acousto-mechanical simulation framework of acoustic radiation force impulse-based SWE was created to elucidate the effect of material viscosity and shear modulus on SWE estimates. Shear modulus estimates exhibited errors up to 72% when a numerical viscoelastic phantom was assessed as linearly elastic. Shear modulus estimates of polyvinyl alcohol phantoms between rheometry and SWE following the Kelvin-Voigt viscoelastic model assumptions were not significantly different. However, the percentage difference in shear modulus estimates between rheometry and SWE using the linear-elastic assumption was 50.1%-62.1%. In ex vivo liver, the percentage difference in shear modulus estimates between linear-elastic and viscoelastic methods was 76.1%. These findings provide a direct and systematic quantification of the potential error introduced when viscoelastic tissues are imaged with SWE following the linear-elastic assumption. This work emphasizes the need to utilize viscoelasticity estimation methods for developing robust quantitative imaging biomarkers.

Feasibility of Phase Velocity Imaging Using Multi Frequency Oscillation-Shear Wave Elastography

OBJECTIVE

To assess viscoelasticity, a pathologically relevant biomarker, shear wave elastography (SWE) generally uses phase velocity (PV) dispersion relationship generated via pulsed acoustic radiation force (ARF) excitation pulse. In this study, a multi-frequency oscillation (MFO)- excitation pulse with higher weight to higher frequencies is proposed to generate PV images via the generation of motion with energy concentrated at the target frequencies in contrast to the broadband frequency motion generated in pulsed SWE (PSWE).

METHODS

The feasibility of MFO-SWE to generate PV images at 100 to 1000 Hz in steps of 100 Hz was investigated by imaging 6 and 70 kPa inclusions with 6.5 and 10.4 mm diameter and ex vivo bovine liver with and without the presence of an aberration layer and chicken muscle ex vivo, and 4T1 mouse breast tumor, in vivo with comparisons to PSWE.

RESULTS

MFO-SWE-derived CNR was statistically higher than PSWE for 6 kPa (both with and without aberration) and 70 kPa (with aberration) inclusions and derived SNR of the liver was statistically higher than PSWE at higher frequency (600-1000 Hz). Quantitatively, at 600-1000 Hz, MFO-SWE improved CNR of inclusions (without and with) aberration on an average by (8.2 and 156)% and of the tumor by 122%, respectively, and improved SNR of the liver (without and with) aberration by (20.2 and 51.5)% and of chicken muscle by 72%, respectively compared to the PSWE.

CONCLUSIONS AND SIGNIFICANCE

These results indicate the advantages of MFO-SWE to improve PV estimation at higher frequencies which could improve viscoelasticity quantification and feature delineation.

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.

Investigation of the relationship between tensile viscoelasticity and unloaded ultrasound shear wave measurements in ex vivo tendon

Mechanical properties of biological tissues are of key importance for proper function and in situ methods for mechanical characterization are sought after in the context of both medical diagnosis as well as understanding of pathophysiological processes. Shear wave elastography (SWE) and accompanying physical modelling methods provide valid estimates of stiffness in quasi-linear viscoelastic, isotropic tissue but suffer from limitations in assessing non-linear viscoelastic or anisotropic material, such as tendon. Indeed, mathematical modelling predicts the longitudinal shear wave velocity to be unaffected by the tensile but rather the shear viscoelasticity. Here, we employ a heuristic experimental testing approach to the problem to assess the most important potential confounders, namely tendon mass density and diameter, and to investigate associations between tendon tensile viscoelasticity with shear wave descriptors. Small oscillatory testing of animal flexor tendons at two baseline stress levels over a large frequency range comprehensively characterized tensile viscoelastic behavior. A broad set of shear wave descriptors was retrieved on the unloaded tendon based on high frame-rate plane wave ultrasound after applying an acoustic deformation impulse. Tensile modulus and strain energy dissipation increased logarithmically and linearly, respectively, with the frequency of the applied strain. Shear wave descriptors were mostly unaffected by tendon diameter but were highly sensitive to tendon mass density. Shear wave group and phase velocity showed no association with tensile elasticity or strain rate-stiffening but did show an association with tensile strain energy dissipation. The longitudinal shear wave velocity may not characterize tensile elasticity but rather tensile viscous properties of transversely isotropic collagenous tissues.

Shear-Wave Elastography and Viscosity PLUS for the Assessment of Peripheric Muscles in Healthy Subjects: A Pre- and Post-Contraction Study

Viscosity is a novel parameter, recently introduced in the use of elastographic techniques, correlating to shear-wave dispersion. The purpose of this study was to provide normal reference viscosity values for the peripheral muscles in healthy volunteers. This prospective study included 38 subjects who underwent US examinations between November 2021 and January 2022. Measurements were taken on the calf and the deltoid muscles in both pre- and post-contraction states. The age range was 21–29 years, with a median of 26 years. The SWE and ViPLUS values in the deltoid muscles were significantly higher than in the soleus muscles in both pre- and post-contraction sets (p = 0.002). There were statistically significant differences between the pre- and post-contraction values for both the SWE and ViPLUS values in the subgroup analysis. The ICC estimates and the 95% confidence intervals were based on a mean rating (k = 2), an absolute agreement, and a two-way random-effects model, demonstrating excellent agreement between the measurements taken by the two examiners.

Improved two-point frequency shift power method for measurement of shear wave attenuation

Quantitative assessment of mechanical properties of biological soft tissues is frequently evaluated using a noninvasive modality, called ultrasound shear wave elastography (SWE). SWE typically exerts an acoustic radiation force (ARF) to produce shear waves propagating in the lateral direction for which velocities and attenuations are measured. The tissue viscoelasticity is commonly studied by investigating the shear wave phase velocity curves. Viscoelastic tissue properties can also be characterized through utilizing various shear wave attenuation techniques. In this study, we propose an improved method for measuring the shear wave attenuation, called two-point frequency shift power (2P-FSP), which is an improved version of the two-point frequency shift (2P-FS) method. The technique is fully data driven and does not use a rheological model for mathematical modeling. The 2P-FSP method utilizes the power spectra frequency shift of shear waves measured at two spatial positions, which provides robustness to noise. The conceptual basis for the 2P-FSP is provided and tested with numerical and experimental data. We investigated how the location of the first signal and the distance interval between the two locations influence the shear wave attenuation measurement in the 2P-FSP technique. We utilized the 2P-FSP method on numerical phantom data generated using a finite-difference-based method in tissue-mimicking viscoelastic media. Moreover, we tested the 2P-FSP method with data from custom-made tissue-mimicking viscoelastic phantom experiments, and ex vivo porcine liver. We compared results from the proposed technique with results from 2P-FS and analytical values in the case of simulations. The results showed that the 2P-FSP method provides improved results over the 2P-FS technique for lower signal-to-noise ratio (SNR) and locations farther from the push location considered, and can be used to measure attenuation of viscoelastic soft tissues.

Deep learning in ultrasound elastography imaging: A review

It is known that changes in the mechanical properties of tissues are associated with the onset and progression of certain diseases. Ultrasound elastography is a technique to characterize tissue stiffness using ultrasound imaging either by measuring tissue strain using quasi-static elastography or natural organ pulsation elastography, or by tracing a propagated shear wave induced by a source or a natural vibration using dynamic elastography. In recent years, deep learning has begun to emerge in ultrasound elastography research. In this review, several common deep learning frameworks in the computer vision community, such as multilayer perceptron, convolutional neural network, and recurrent neural network are described. Then, recent advances in ultrasound elastography using such deep learning techniques are revisited in terms of algorithm development and clinical diagnosis. Finally, the current challenges and future developments of deep learning in ultrasound elastography are prospected. This article is protected by copyright. All rights reserved.

The Revisited Frequency-Shift Method for Shear Wave Attenuation Computation and Imaging

Ultrasound (US) shear wave (SW) elastography has been widely studied and implemented on clinical systems to assess the elasticity of living organs. Imaging of SW attenuation reflecting viscous properties of tissues has received less attention. A revisited frequency shift (R-FS) method is proposed to improve the robustness of SW attenuation imaging. Performances are compared with the FS method that we originally proposed and with the two-point frequency shift (2P-FS) and attenuation measuring US SW elastography (AMUSE) methods. In the proposed R-FS method, the shape parameter of the gamma distribution fitting SW spectra is assumed to vary with distance, in contrast to FS. Second, an adaptive random sample consensus (A-RANSAC) line fitting method is used to prevent outlier attenuation values in the presence of noise. Validation was made on ten simulated phantoms with two viscosities (0.5 and 2 Pa [Formula: see text]) and different noise levels (15 to -5 dB), two experimental homogeneous gel phantoms, and six in vivo liver acquisitions on awake ducks (including three normal and three fatty duck livers). According to the conducted experiments, R-FS revealed mean reductions in coefficients of variation (CV) of 62.6% on simulations, 62.5% with phantoms, and 62.3% in vivo compared with FS. Corresponding reductions compared with 2P-FS were 45.4%, 77.1%, and 62.0%, respectively. Reductions in normalized root-mean-square errors for simulations were 63.9% and 48.7% with respect to FS and 2P-FS, respectively.

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.

Nonlinear Characterization of Tissue Viscoelasticity With Acoustoelastic Attenuation of Shear Waves

Shear-wave elastography (SWE) measures shear-wave speed (SWS), which is related to the underlying shear modulus of soft tissue. SWS in soft tissue changes depending on the amount of external strain that soft tissue is subjected to due to the acoustoelastic (AE) phenomenon. In the literature, variations of SWS as a function of applied uniaxial strain were used for nonlinear characterization, assuming soft tissues to be elastic, although soft tissues are indeed viscoelastic in nature. Hence, nonlinear characterization using SWS alone is insufficient. In this work, we use SWS together with shear-wave attenuation (SWA) during incremental quasi-static compressions in order to derive biomechanical characterization based on the AE theory in terms of well-defined storage and loss moduli. As part of this study, we also quantify the effect of applied strain on measurements of SWS and SWA since such confounding effects need to be taken into account when using SWS and/or SWA, e.g., for staging a disease state, while such effects can also serve as an additional imaging biomarker. Our results from tissue-mimicking phantoms with varying oil percentages and ex vivo porcine liver experiments demonstrate the feasibility of our proposed methods. In both experiments, SWA was observed to decrease with applied strain. For 10% compression in ex vivo livers, shear-wave attenuation decreased, on average, by 28% (93 Np/m), while SWS increased, on average, by 20% (0.26 m/s).

Shear Wave Elastography: A Review on the Confounding Factors and Their Potential Mitigation in Detecting Chronic Kidney Disease

Early detection of chronic kidney disease is important to prevent progression of irreversible kidney damage, reducing the need for renal transplantation. Shear wave elastography is ideal as a quantitative imaging modality to detect chronic kidney disease because of its non-invasive nature, low cost and portability, making it highly accessible. However, the complexity of the kidney architecture and its tissue properties give rise to various confounding factors that affect the reliability of shear wave elastography in detecting chronic kidney disease, thus limiting its application to clinical trials. The objective of this review is to highlight the confounding factors presented by the complex properties of the kidney, in addition to outlining potential mitigation strategies, along with the prospect of increasing the versatility and reliability of shear wave elastography in detecting chronic kidney disease.

Viscoelasticity Imaging of Biological Tissues and Single Cells Using Shear Wave Propagation

Changes in biomechanical properties of biological soft tissues are often associated with physiological dysfunctions. Since biological soft tissues are hydrated, viscoelasticity is likely suitable to represent its solid-like behavior using elasticity and fluid-like behavior using viscosity. Shear wave elastography is a non-invasive imaging technology invented for clinical applications that has shown promise to characterize various tissue viscoelasticity. It is based on measuring and analyzing velocities and attenuations of propagated shear waves. In this review, principles and technical developments of shear wave elastography for viscoelasticity characterization from organ to cellular levels are presented, and different imaging modalities used to track shear wave propagation are described. At a macroscopic scale, techniques for inducing shear waves using an external mechanical vibration, an acoustic radiation pressure or a Lorentz force are reviewed along with imaging approaches proposed to track shear wave propagation, namely ultrasound, magnetic resonance, optical, and photoacoustic means. Then, approaches for theoretical modeling and tracking of shear waves are detailed. Following it, some examples of applications to characterize the viscoelasticity of various organs are given. At a microscopic scale, a novel cellular shear wave elastography method using an external vibration and optical microscopy is illustrated. Finally, current limitations and future directions in shear wave elastography are presented.

Phase Velocity Estimation With Expanded Bandwidth in Viscoelastic Phantoms and Tissues

Ultrasound shear wave elastography (SWE) is a technique used to measure mechanical properties to evaluate healthy and pathological soft tissues. SWE typically employs an acoustic radiation force (ARF) to generate laterally propagating shear waves that are tracked in the spatiotemporal domains, and algorithms are used to estimate the wave velocity. The tissue viscoelasticity is often examined through analyzing the shear wave phase velocity dispersion curves, which is the variation of phase velocity with frequency or wavelength. A number of available methods to estimate dispersion exist, which can differ in resolution and variance. Moreover, most of these techniques reconstruct dispersion curves for a limited frequency band. In this work, we propose a novel method used for dispersion curve calculation. Our unique approach uses a generalized Stockwell transformation combined with a slant frequency-wavenumber analysis (GST-SFK). We tested the GST-SFK method on numerical phantom data generated using a finite-difference-based method in tissue-mimicking viscoelastic media. In addition, we evaluated the method on numerical shear wave motion data with different amounts of white Gaussian noise added. Additionally, we performed tests on data from custom-made tissue-mimicking viscoelastic phantom experiments, ex vivo porcine liver measurements, and in vivo liver tissue experiments. We compared results from our method with two other techniques used for estimating shear wave phase velocity: the two-dimensional Fourier transform (2D-FT) and the eigenvector (EV) method. Tests carried out revealed that the GST-SFK method provides dispersion curve estimates with lower errors over a wider frequency band in comparison to the 2D-FT and EV methods. In addition, the GST-SFK provides expanded bandwidth by a factor of two or more to be used for phase velocity estimation, which is meaningful for a tissue dispersion analysis in vivo.

Multiparametric in vivo ultrasound shear wave viscoelastography on farm-raised fatty duck livers: human radiology imaging applied to food sciences

Nine mulard ducks that were being raised for foie gras (steatosis) production went through in vivo shear wave (SW) elastography imaging of their liver during the force-feeding period to investigate changes in liver tissue characteristics. A total of 4 imaging sessions at an interval of 3 to 4 d were conducted at the farm on each animal. Three ducks were sacrificed at the second, third, and fourth imaging sessions for histopathology analysis of all animals at these time points. Six SW elastography parameters were evaluated: SW speed, SW attenuation, SW dispersion, Young's modulus, viscosity, and shear modulus. Shear waves of different frequencies propagate with different phase velocities. Thus, SW speed and other dependent parameters such as Young's modulus, viscosity, and shear modulus were computed at 2 frequencies: 75 and 202 Hz. Each parameter depicted a statistically significant trend along the force-feeding process (P-values between 0.001 and 0.0001). The fat fraction of the liver increased over the 12-day period of feeding. All parameters increased monotonically over time at 75 Hz, whereas modal relations were seen at 202 Hz. Shear wave dispersion measured between 75 and 202 Hz depicted a plateau from day 5. Based on this validation, proposed imaging methods are aimed to be used in the future on naturally fed ducks and geese.

Fat and fibrosis as confounding cofactors in viscoelastic measurements of the liver

Elastography provides significant information on staging of fibrosis in patients with liver disease and may be of some value in assessing steatosis. However, there remain questions as to the role of steatosis and fibrosis as cofactors influencing the viscoelastic measurements of liver tissues, particularly shear wave speed (SWS) and shear wave attenuation (SWA). In this study, by employing the theory of composite elastic media as well as two independent experimental measurements on oil-in-gelatin phantoms and also finite element simulations, it is consistently shown that fat and fibrosis jointly influence the SWS and SWA measurements. At a constant level of fat, fibrosis stages can influence the SWA by factors of 2-4. Moreover, the rate of increase in SWA with increasing fat is strongly influenced by the stages of fibrosis; softer background cases (low fibrosis stages) have higher rate of SWA increase with fat than those with stiffer moduli (higher fibrosis stages). Meanwhile, SWS results are influenced by the presence of fat, however the degree of variability is more subtle. The results indicate the importance of jointly considering fat and fibrosis as contributors to SWS and SWA measurements in complex liver tissues and in the design and interpretation of clinical trials.

Resonant acoustic rheometry for non-contact characterization of viscoelastic biomaterials

Resonant Acoustic Rheometry (RAR) is a new, non-contact technique to characterize the mechanical properties of soft and viscoelastic biomaterials, such as hydrogels, that are used to mimic the extracellular matrix in tissue engineering. RAR uses a focused ultrasound pulse to generate a microscale perturbation at the sample surface and tracks the ensuing surface wave using pulse-echo ultrasound. The frequency spectrum of the resonant surface waves is analyzed to extract viscoelastic material properties. In this study, RAR was used to characterize fibrin, gelatin, and agarose hydrogels. Single time point measurements of gelled samples with static mechanical properties showed that RAR provided consistent quantitative data and measured intrinsic material characteristics independent of ultrasound parameters. RAR was also used to longitudinally track dynamic changes in viscoelastic properties over the course of fibrin gelation, revealing distinct phase and material property transitions. Application of RAR was verified using finite element modeling and the results were validated against rotational shear rheometry. Importantly, RAR circumvents some limitations of conventional rheology methods and can be performed in a high-throughput manner using conventional labware. Overall, these studies demonstrate that RAR can be a valuable tool to noninvasively quantify the viscoelastic mechanical properties of soft hydrogel biomaterials.

The Physical Microenvironment of Tumors: Characterization and Clinical Impact

The tumor microenvironment plays a critical role in tumorigenesis and metastasis. As tightly controlled extracellular matrix homeostasis is lost during tumor progression, a dysregulated extracellular matrix can significantly alter cellular phenotype and drive malignancy. Altered physical properties of the tumor microenvironment alter cancer cell behavior, limit delivery and efficacy of therapies, and correlate with tumorigenesis and patient prognosis. The physical features of the extracellular matrix during tumor progression have been characterized; however, a wide range of methods have been used between studies and cancer types resulting in a large range of reported values. Here, we discuss the significant mechanical and structural properties of the tumor microenvironment, summarizing their reported values and clinical impact across cancer type and grade. We attempt to integrate the values in the literature to identify sources of reported differences and commonalities to better understand how aberrant extracellular matrix dynamics contribute to cancer progression. An intimate understanding of altered matrix properties during malignant transformation will be crucial in effectively detecting, monitoring, and treating cancer.

Viscoelastic Response Ultrasound Derived Relative Elasticity and Relative Viscosity Reflect True Elasticity and Viscosity: In Silico and Experimental Demonstration

Viscoelastic response (VisR) ultrasound characterizes the viscoelastic properties of tissue by fitting acoustic radiation force (ARF)-induced displacements in the region of ARF excitation to a 1-D mass-spring-damper (MSD) model. Elasticity and viscosity are calculated separately but relative to the applied ARF amplitude. We refer to these parameters as "relative elasticity (RE)" and "relative viscosity (RV)." We herein test the hypothesis that RE and RV linearly correlate to true elasticity and viscosity in tissue. VisR imaging was simulated in 144 homogeneous viscoelastic materials with varying elasticities and viscosities. Derived RE linearly correlated with material elasticity and varied by an average of 2.52% when the material viscosity changed from 0.1 to 1.3 Pa · s. Derived RV linearly correlated with material viscosity but varied by an average of 102.5% when material elasticity changed from 3.33 to 20 kPa. The effect of elasticity on RV measurement was compensated using the slope of the linear relationship between RV and natural frequency ( ωtextn ). After compensation, RV [Formula: see text] (elasticity compensated RV) linearly correlated with material viscosity and varied by less than 1.00% on average when the modeled shear elastic modulus changed from 3.3 to 20 kPa. In addition to elasticity compensation, variation in ARF amplitude over depth was compensated, yielding REDC and [Formula: see text]. REDC and [Formula: see text] successfully contrasted elastic and viscous inclusions, respectively, in three simulated phantoms. Experimentally, in the homogeneous oil-in-gelatin phantoms and excised livers, REDC linearly correlated with shear wave dispersion ultrasound vibrometry (SDUV) derived shear elastic modulus, and [Formula: see text] linearly correlated with SDUV-derived shear viscosity. In excised livers containing viscoelastic oil-in-gelatin inclusions, the inclusions were successfully contrasted from the liver background by both REDC and [Formula: see text]. These results suggest that RE and RV are relevant for qualitatively assessing the elastic and viscous properties of tissue.

Two-Point Frequency Shift Method for Shear Wave Attenuation Measurement

Ultrasound shear wave elastography (SWE) is an increasingly used noninvasive modality for quantitative evaluation of tissue mechanical properties. SWE typically uses an acoustic radiation force to produce laterally propagating shear waves that are tracked in the spatial and temporal domains, in order to obtain the wave velocity. One of the ways to study the viscoelasticity is through studying the shear wave phase velocity dispersion curves. Shear wave attenuation can be also characterized in viscoelastic tissues with methods that use multiple lateral data samples. In this article, we present an alternative method for measuring the shear wave attenuation without using a rheological model two-point frequency shift (2P-FS). The technique uses information related to the amplitude spectra FS of shear waves measured at only two lateral locations. The theoretical basis for the 2P-FS is derived and validated. We examined how the first signal position and the distance between the two locations affect the shear wave attenuation estimation in the 2P-FS method. We tested this new method on digital phantom data created using the local interaction simulation approach (LISA) in viscoelastic media. Moreover, we tested data acquired from custom-made tissue-mimicking viscoelastic phantom experiments and ex vivo porcine liver measurements. We compared results from the 2P-FS method with the other two techniques used for assessing a shear wave attenuation: the FS-based method and the attenuation-measuring ultrasound shear wave elastography (AMUSE) technique. In addition, we evaluated the 2P-FS algorithm with different levels of added white Gaussian noise to the shear wave particle velocity using numerical phantoms. Tests conducted showed that the 2P-FS method gives robust results based on only two measurements and can be used to measure attenuation of viscoelastic soft tissues.

Spectral Quantification of Nonlinear Elasticity Using Acoustoelasticity and Shear-Wave Dispersion

Tissue biomechanical properties are known to be sensitive to pathological changes. Accordingly, various techniques have been developed to estimate tissue mechanical properties. Shear-wave elastography (SWE) measures shear-wave speed (SWS) in tissues, which can be related to shear modulus. Although viscosity or stress-strain nonlinearity may act as confounder of SWE, their explicit characterization may also provide additional information about tissue composition as a contrast modality. Viscosity can be related to frequency dispersion of SWS, which can be characterized using multi-frequency measurements, herein called spectral SWE (SSWE). Additionally, nonlinear shear modulus can be quantified and parameterized based on SWS changes with respect to applied stress, a phenomenon called acoustoelasticity (AE). In this work, we characterize the nonlinear parameters of tissue as a function of excitation frequency by utilizing both AE and SSWE together. For this, we apply incremental amounts of quasi-static stress on a medium, while imaging and quantifying SWS dispersion via SSWE. Results from phantom and ex vivo porcine liver experiments demonstrate the feasibility of measuring frequency-dependent nonlinear parameters using the proposed method. SWS propagation in porcine liver tissue was observed to change from 1.8 m/s at 100 Hz to 3.3 m/s at 700 Hz, while increasing by approximately 25% from a strain of 0% to 12% across these frequencies.

Reconstruction of Viscosity Maps in Ultrasound Shear Wave Elastography

Change in viscoelastic properties of biological tissues may often be symptomatic of a dysfunction that can be correlated to tissue pathology. Shear wave elastography is an imaging method mainly used to assess stiffness but with the potential to measure viscoelasticity of biological tissues. This can enable tissue characterization; and thus, can be used as a marker to improve diagnosis of pathological lesions. In this study, a frequency-shift method based framework is presented for the reconstruction of viscosity by analyzing the spectral properties of acoustic radiation force-induced shear waves. The aim of the study was to investigate the feasibility of viscosity reconstruction maps in homogeneous as well as heterogeneous samples. Experiments were performed in four in vitro phantoms, two ex vivo porcine liver samples, two ex vivo fatty duck liver samples, and one in vivo fatty goose liver. Successful viscosity maps were reconstructed in homogeneous and heterogeneous phantoms with embedded mechanical inclusions having different geometries. Quantitative values of viscosity obtained for two porcine liver tissues, two fatty duck liver samples, and one goose fatty liver were (mean ± SD) 0.61 ± 0.21, 0.52 ± 0.35; 1.28 ± 0.54, 1.36 ± 0.73, and 1.67 ± 0.70 Pa.s, respectively.

Shear wave propagation in viscoelastic media: validation of an approximate forward model

Many approaches to elastography incorporate shear waves; in some systems these are produced by acoustic radiation force (ARF) push pulses. Understanding the shape and decay of propagating shear waves in lossy tissues is key to obtaining accurate estimates of tissue properties, and so analytical models have been proposed. In this paper, we reconsider a previous analytical model with the goal of obtaining a computationally straightforward and efficient equation for the propagation of shear waves from a focal push pulse. Next, this model is compared with an experimental optical coherence tomography (OCT) system and with finite element models, in two viscoelastic materials that mimic tissue. We find that the three different cases-analytical model, finite element model, and experimental results-demonstrate reasonable agreement within the subtle differences present in their respective conditions. These results support the use of an efficient form of the Hankel transform for both lossless (elastic) and lossy (viscoelastic) media, and for both short (impulsive) and longer (extended) push pulses that can model a range of experimental conditions.

Analysis of Transient Shear Wave in Lossy Media

The propagation of shear waves from impulsive forces is an important topic in elastography. Observations of shear wave propagation can be obtained with numerous clinical imaging systems. Parameter estimations of the shear wave speed in tissues, and more generally the viscoelastic parameters of tissues, are based on some underlying models of shear wave propagation. The models typically include specific choices of the spatial and temporal shape of the impulsive force and the elastic or viscoelastic properties of the medium. In this work, we extend the analytical treatment of 2-D shear wave propagation in a biomaterial. The approach applies integral theorems relevant to the solution of the generalized Helmholtz equation, and does not depend on a specific rheological model of the tissue's viscoelastic properties. Estimators of attenuation and shear wave speed are derived from the analytical solutions, and these are applied to an elastic phantom, a viscoelastic phantom and in vivo liver using a clinical ultrasound scanner. In these samples, estimated shear wave group velocities ranged from 1.7 m/s in the liver to 2.5 m/s in the viscoelastic phantom, and these are lower-bounded by independent measurements of phase velocity.

Application of a forward model of axisymmetric shear wave propagation in viscoelastic media to shear wave elastography

A simple but general solution of Navier's equation for axisymmetric shear wave propagation in a homogeneous isotropic viscoelastic medium is presented. It is well-suited for use as a forward model for some acoustic radiation force impulse based shear wave elastography applications because it does not require precise knowledge of the strength of the source, nor its spatial or temporal distribution. Instead, it depends on two assumptions: (1) the source distribution is axisymmetric and confined to a small region near the axis of symmetry, and (2) the propagation medium is isotropic and homogeneous. The model accounts for the vector polarization of shear waves and exactly represents geometric spreading of the shear wavefield, whether spherical, cylindrical, or neither. It makes no assumption about the frequency dependence of material parameters, i.e., it is material-model independent. Validation using measured shear wavefields excited by acoustic radiation force in a homogeneous gelatin sample show that the model accounts for well over 90% of the measured wavefield "energy." An optimal fit of the model to simulated shear wavefields with noise in a homogeneous viscoelastic medium enables estimation of both the shear storage modulus and shear wave attenuation to within 1%.

Impact of Acoustic Radiation Force Excitation Geometry on Shear Wave Dispersion and Attenuation Estimates

Shear wave elasticity imaging (SWEI) characterizes the mechanical properties of human tissues to differentiate healthy from diseased tissue. Commercial scanners tend to reconstruct shear wave speeds for a region of interest using time-of-flight methods reporting a single shear wave speed (or elastic modulus) to the end user under the assumptions that tissue is elastic and shear wave speeds are not dependent on the frequency content of the shear waves. Human tissues, however, are known to be viscoelastic, resulting in dispersion and attenuation. Shear wave spectroscopy and spectral methods have been previously reported in the literature to quantify shear wave dispersion and attenuation, commonly making an assumption that the acoustic radiation force excitation acts as a cylindrical source with a known geometric shear wave amplitude decay. This work quantifies the bias in shear dispersion and attenuation estimates associated with making this cylindrical wave assumption when applied to shear wave sources with finite depth extents, as commonly occurs with realistic focal geometries, in elastic and viscoelastic media. Bias is quantified using analytically derived shear wave data and shear wave data generated using finite-element method models. Shear wave dispersion and attenuation bias (up to 15% for dispersion and 41% for attenuation) is greater for more tightly focused acoustic radiation force sources with smaller depths of field relative to their lateral extent (height-to-width ratios <16). Dispersion and attenuation errors associated with assuming a cylindrical geometric shear wave decay in SWEI can be appreciable and should be considered when analyzing the viscoelastic properties of tissues with acoustic radiation force source distributions with limited depths of field.

Acoustic radiation force induced resonance elastography of coagulating blood: theoretical viscoelasticity modeling and ex-vivo experimentation

Deep vein thrombosis is a common vascular disease that can lead to pulmonary embolism and death. The early diagnosis and clot age staging are important parameters for reliable therapy planning. This article presents an acoustic radiation force induced resonance elastography method for the viscoelastic characterization of clotting blood. The physical concept of this method relies on the mechanical resonance of the blood clot occurring at specific frequencies. Resonances are induced by focusing ultrasound beams inside the sample under investigation. Coupled to an analytical model of wave scattering, the ability of the proposed method to characterize the viscoelasticity of a mimicked venous thrombosis in the acute phase is demonstrated. Experiments with a gelatin-agar inclusion sample of known viscoelasticity are performed for validation and establishment of the proof of concept. In addition, an inversion method is applied in-vitro for the kinetic monitoring of the blood coagulation process of six human blood samples obtained from two volunteers. The computed elasticity and viscosity values of blood samples at the end of the 90 min kinetics were estimated at 411 ± 71 Pa and 0.25 ± 0.03 Pa.s for volunteer #1, and 387 ± 35 Pa and 0.23 ± 0.02 Pa.s for volunteer #2, respectively. The proposed method allowed reproducible time-varying thrombus viscoelastic measurements from samples having physiological dimensions.

Viscoelastic properties of human pancreatic tumors and in vitro constructs to mimic mechanical properties

Pancreatic ductal adenocarcinoma (PDAC) is almost universally fatal, in large part due to a protective fibrotic barrier generated by tumor-associated stromal (TAS) cells. This barrier is thought to promote cancer cell survival and confounds attempts to develop effective therapies. We present a 3D in vitro system that replicates the mechanical properties of the PDAC microenvironment, representing an invaluable tool for understanding the biology of the disease. Mesoscale indentation quantified viscoelastic metrics of resected malignant tumors, inflamed chronic pancreatitis regions, and histologically normal tissue. Both pancreatitis (2.15 ± 0.41 kPa, Mean ± SD) and tumors (5.46 ± 3.18 kPa) exhibit higher Steady-State Modulus (SSM) than normal tissue (1.06 ± 0.25 kPa; p < .005). The average viscosity of pancreatitis samples (63.2 ± 26.7 kPa·s) is significantly lower than that of both normal tissue (252 ± 134 kPa·s) and tumors (349 ± 222 kPa·s; p < .005). To mimic this remodeling behavior, PDAC and TAS cells were isolated from human PDAC tumors. Conditioned medium from PDAC cells was used to culture TAS-embedded collagen hydrogels. After 7 days, TAS-embedded gels in control medium reached SSM (1.45 ± 0.12 kPa) near normal pancreas, while gels maintained with conditioned medium achieved higher SSM (3.38 ± 0.146 kPa) consistent with tumors. Taken together, we have demonstrated an in vitro system that recapitulates in vivo stiffening of PDAC tumors. In addition, our quantification of viscoelastic properties suggests that elastography algorithms incorporating viscosity may be able to more accurately distinguish between pancreatic cancer and pancreatitis.

STATEMENT OF SIGNIFICANCE

Understanding tumor-stroma crosstalk in pancreatic ductal adenocarcinoma (PDAC) is challenged by a lack of stroma-mimicking model systems. To design appropriate models, pancreatic tissue must be characterized with a method capable of evaluating in vitro models as well. Our indentation-based characterization tool quantified the distinct viscoelastic signatures of inflamed resections from pancreatitis, tumors from PDAC, and otherwise normal tissue to inform development of mechanically appropriate engineered tissues and scaffolds. We also made progress toward a 3D in vitro system that recapitulates mechanical properties of tumors. Our in vitro model of stromal cells in collagen and complementary characterization system can be used to investigate mechanisms of cancer-stroma crosstalk in PDAC and to propose and test innovative therapies.

Quantitative characterization of viscoelastic behavior in tissue-mimicking phantoms and ex vivo animal tissues

Viscoelasticity of soft tissue is often related to pathology, and therefore, has become an important diagnostic indicator in the clinical assessment of suspect tissue. Surgeons, particularly within head and neck subsites, typically use palpation techniques for intra-operative tumor detection. This detection method, however, is highly subjective and often fails to detect small or deep abnormalities. Vibroacoustography (VA) and similar methods have previously been used to distinguish tissue with high-contrast, but a firm understanding of the main contrast mechanism has yet to be verified. The contributions of tissue mechanical properties in VA images have been difficult to verify given the limited literature on viscoelastic properties of various normal and diseased tissue. This paper aims to investigate viscoelasticity theory and present a detailed description of viscoelastic experimental results obtained in tissue-mimicking phantoms (TMPs) and ex vivo tissues to verify the main contrast mechanism in VA and similar imaging modalities. A spherical-tip micro-indentation technique was employed with the Hertzian model to acquire absolute, quantitative, point measurements of the elastic modulus (E), long term shear modulus (η), and time constant (τ) in homogeneous TMPs and ex vivo tissue in rat liver and porcine liver and gallbladder. Viscoelastic differences observed between porcine liver and gallbladder tissue suggest that imaging modalities which utilize the mechanical properties of tissue as a primary contrast mechanism can potentially be used to quantitatively differentiate between proximate organs in a clinical setting. These results may facilitate more accurate tissue modeling and add information not currently available to the field of systems characterization and biomedical research.

Ultrafast imaging of cell elasticity with optical microelastography

In wave physics, and especially seismology, uncorrelated vibrations could be exploited using “noise correlation” tools to reconstruct images of a medium. By using a high-frequency vibration, a high-speed tracking device, and a reconstruction technique based on temporal correlations of travelling waves we conceptualized an optical microelastography technique to map elasticity of internal cellular structures. This technique, unlike other methods, can provide an elasticity image in less than a millisecond, thus opening the possibility of studying dynamic cellular processes and elucidating new mechanocellular properties. We call this proposed technique “cell quake elastography.”

Elasticity is a fundamental cellular property that is related to the anatomy, functionality, and pathological state of cells and tissues. However, current techniques based on cell deformation, atomic force microscopy, or Brillouin scattering are rather slow and do not always accurately represent cell elasticity. Here, we have developed an alternative technique by applying shear wave elastography to the micrometer scale. Elastic waves were mechanically induced in live mammalian oocytes using a vibrating micropipette. These audible frequency waves were observed optically at 200,000 frames per second and tracked with an optical flow algorithm. Whole-cell elasticity was then mapped using an elastography method inspired by the seismology field. Using this approach we show that the elasticity of mouse oocytes is decreased when the oocyte cytoskeleton is disrupted with cytochalasin B. The technique is fast (less than 1 ms for data acquisition), precise (spatial resolution of a few micrometers), able to map internal cell structures, and robust and thus represents a tractable option for interrogating biomechanical properties of diverse cell types.

Forward and inverse viscoelastic wave scattering by irregular inclusions for shear wave elastography

Inversion methods in shear wave elastography use simplifying assumptions to recover the mechanical properties of soft tissues. Consequently, these methods suffer from artifacts when applied to media containing strong stiffness contrasts, and do not provide a map of the viscosity. In this work, the shear wave field recorded inside and around an inclusion was used to estimate the viscoelastic properties of the inclusion and surrounding medium, based on an inverse problem approach assuming local homogeneity of both media. An efficient semi-analytical method was developed to model the scattering of an elastic wave by an irregular inclusion, based on a decomposition of the field by Bessel functions and on a decomposition of the boundaries as Fourier series. This model was validated against finite element modeling. Shear waves were experimentally induced by acoustic radiation force in soft tissue phantoms containing stiff and soft inclusions, and the displacement field was imaged at a high frame rate using plane wave imaging. A nonlinear least-squares algorithm compared the model to the experimental data and adjusted the geometrical and mechanical parameters. The estimated shear storage and loss moduli were in good agreement with reference measurements, as well as the estimated inclusion shape. This approach provides an accurate estimation of geometry and viscoelastic properties for a single inclusion in a homogeneous background in the context of radiation force elastography.

A Frequency-Shift Method to Measure Shear-Wave Attenuation in Soft Tissues

In vivo quantification of shear-wave attenuation in soft tissues may help to better understand human tissue rheology and lead to new diagnostic strategies. Attenuation is difficult to measure in acoustic radiation force elastography because the shear-wave amplitude decreases due to a combination of diffraction and viscous attenuation. Diffraction correction requires assuming a cylindrical wavefront and an isotropic propagation medium, which may not be the case in some applications. In this paper, the frequency-shift method, used in ultrasound imaging and seismology, was adapted for shear-wave attenuation measurement in elastography. This method is not sensitive to diffraction effects. For a linear frequency dependence of the attenuation, a closed-form relation was obtained between the decrease in the peak frequency of the gamma-distributed wave amplitude spectrum and the attenuation coefficient of the propagation medium. The proposed method was tested against a plane-wave reference method in homogeneous agar-gelatin phantoms with 0%, 10%, and 20% oil concentrations, and hence different attenuations of 0.117, 0.202, and 0.292 [Formula: see text]/Hz, respectively. Applicability to biological tissues was demonstrated with two ex vivo porcine liver samples (0.79 and 1.35 [Formula: see text]/Hz) and an in vivo human muscle, measured along (0.43 [Formula: see text]/Hz) and across (1.77 [Formula: see text]/Hz) the tissue fibers. In all cases, the data supported the assumptions of a gamma-distributed spectrum for the source and linear frequency attenuation for the tissue. This method provides tissue attenuation, which is relevant diagnostic information to model viscosity, in addition to shear-wave velocity used to assess elasticity. Data processing is simple and could be performed automatically in real time for clinical applications.

Current status of musculoskeletal application of shear wave elastography

Ultrasonography (US) is a very powerful diagnostic modality for the musculoskeletal system due to the ability to perform real-time dynamic high-resolution examinations with the Doppler technique. In addition to acquiring morphologic data, we can now obtain biomechanical information by quantifying the elasticity of the musculoskeletal structures with US elastography. The earlier diagnosis of degeneration and the ability to perform follow-up evaluations of healing and the effects of treatment are possible. US elastography enables a transition from US-based inspection to US-based palpation in order to diagnose the characteristics of tissue. Shear wave elastography is considered the most suitable type of US elastography for the musculoskeletal system. It is widely used for tendons, ligaments, and muscles. It is important to understand practice guidelines in order to enhance reproducibility. Incorporating viscoelasticity and overcoming inconsistencies among manufacturers are future tasks for improving the capabilities of US elastography.