Shear wave dispersion ultrasonic vibrometry for measuring prostate shear stiffness and viscosity: an in vitro pilot study

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.

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.

Quantifying the Impact of Cancer on the Viscoelastic Properties of the Prostate Gland using a Quasi-Linear Viscoelastic Model

Pathological disorders can alter the mechanical properties of biological tissues, and studying such changes can help to better understand the disease progression. The prostate gland is no exception, as previous studies have shown that cancer can affect its mechanical properties. However, most of these studies have focused on the elastic properties of the tissue and have overlooked the impact of cancer on its viscous response. To address this gap, we used a quasi-linear viscoelastic model to investigate the impact of cancer on both the elastic and viscous characteristics of the prostate gland. By comparing the viscoelastic properties of segments influenced by cancer and those unaffected by cancer in 49 fresh prostates, removed within two hours after prostatectomy surgery, we were able to determine the influence of cancer grade and tumor volume on the tissue. Our findings suggest that tumor volume significantly affects both the elastic modulus and viscosity of the prostate (p-value less than 2%). Specifically, we showed that cancer increases Young's modulus and shear relaxation modulus by 20%. These results have implications for using mechanical properties of the prostate as a potential biomarker for cancer. However, developing an in vivo apparatus to measure these properties remains a challenge that needs to be addressed in future research. STATEMENT OF SIGNIFICANCE: This study is the first to explore how cancer impacts the mechanical properties of prostate tissues using a quasi-linear viscoelastic model. We examined 49 fresh prostate samples collected immediately after surgery and correlated their properties with cancer presence identified in pathology reports. Our results demonstrate a 20% change in the viscoelastic properties of the prostate due to cancer. We initially validated our approach using tissue-mimicking phantoms and then applied it to differentiate between cancerous and normal prostate tissues. These findings offer potential for early cancer detection by assessing these properties. However, conducting these tests in vivo remains a challenge for future research.

Full waveform inversion for arterial viscoelasticity

Objective. Arterial viscosity is emerging as an important biomarker, in addition to the widely used arterial elasticity. This paper presents an approach to estimate arterial viscoelasticity using shear wave elastography (SWE).Approach. While dispersion characteristics are often used to estimate elasticity from SWE data, they are not sufficiently sensitive to viscosity. Driven by this, we develop a full waveform inversion (FWI) methodology, based on directly matching predicted and measured wall velocity in space and time, to simultaneously estimate both elasticity and viscosity. Specifically, we propose to minimize an objective function capturing the correlation between measured and predicted responses of the anterior wall of the artery.Results. The objective function is shown to be well-behaving (generally convex), leading us to effectively use gradient optimization to invert for both elasticity and viscosity. The resulting methodology is verified with synthetic data polluted with noise, leading to the conclusion that the proposed FWI is effective in estimating arterial viscoelasticity.Significance. Accurate estimation of arterial viscoelasticity, not just elasticity, provides a more precise characterization of arterial mechanical properties, potentially leading to a better indicator of arterial health.

High frequency ultrasound vibrational shear wave elastography for preclinical research

Preclinical evaluation of novel therapies using models of cancer is an important tool in cancer research, where imaging can provide non-invasive tools to characterise the internal structure and function of tumours. The short propagation paths when imaging tumours and organs in small animals allow the use of high frequencies for both ultrasound and shear waves, providing the opportunity for high-resolution shear wave elastography and hence its use for studying the heterogeneity of tissue elasticity, where heterogeneity may be a predictor of tissue response. Here we demonstrate vibrational shear wave elastography (VSWE) using a mechanical actuator to produce high frequency (up to 1000 Hz) shear waves in preclinical tumours, an alternative to the majority of preclinical ultrasound SWE studies where an acoustic radiation force impulse is required to create a relatively low-frequency broad-band shear-wave pulse. We implement VSWE with a high frequency (17.8 MHz) probe running a focused line-by-line ultrasound imaging sequence which as expected was found to offer improved detection of 1000 Hz shear waves over an ultrafast planar wave imaging sequence in a homogenous tissue-mimicking phantom. We test the VSWE in anex vivotumour xenograft, demonstrating the ability to detect shear waves up to 10 mm from the contactor position at 1000 Hz. By reducing the kernel size used for shear wave speed estimation to 1 mm we are able to produce shear wave speed images with spatial resolution of this order. Finally, we present VSWE data from xenograft tumoursin vivo, demonstrating the feasibility of the technique in mice under isoflurane sedation. Mean shear wave speeds in the tumours are in good agreements with those reported by previous authors. Characterising the frequency dependence of shear wave speed demonstrates the potential to quantify the viscoelastic properties of tumoursin vivo.

Experimental evidence of shear waves in fractional viscoelastic rheological models

Fractional viscoelastic rheological models, such as the Kelvin Voigt Fractional Derivative model, have been proposed in the literature for modelling shear wave propagation in soft tissue. In this article, our previously developed wave propagation model for transluminal propagation based on a Kelvin Voigt Fractional Derivative wave equation is experimentally validated. The transluminal procedure uses the transmission and detection of shear waves through the luminal wall. The model was compared against high-speed camera observations in translucent elastography phantoms with similar viscoelastic properties to prostate tissue. An ad hoc cross-correlation procedure was used to reconstruct the angular displacement from the high-speed camera observations. Rheometry and shear wave elastography were used for characterising the shear wave velocity dispersion curve for the phantoms. Fractional viscoelastic properties were derived after fitting the dispersion curve to its analytical expression. Propagation features and amplitude spectra from simulations and high-speed camera observations were compared. The obtained results indicate that the model replicates the experimental observations with acceptable accuracy. The model presented here provides a useful tool to model transluminal procedures based on wave propagation and its interaction with the mechanical properties of the tissue outside the lumen.

Broadband-excitation-based mechanical spectroscopy of highly viscous tissue-mimicking phantoms

Standard rheometers assess mechanical properties of viscoelastic samples up to 100 Hz, which often hinders the assessment of the local-scale dynamics. We demonstrate that high-frequency analysis can be achieved by inducing broadband waves and monitoring their media-dependent propagation using optical coherence tomography. Here, we present a new broadband wave analysis based on two-dimensional Fourier transformation. We validated this method by comparing the mechanical parameters to monochromatic excitation and a standard oscillatory test data. Our method allows for high-frequency mechanical spectroscopy, which could be used to investigate the local-scale dynamics of different biological tissues and the influence of diseases on their microstructure.

Wave Propagation in a Fractional Viscoelastic Tissue Model: Application to Transluminal Procedures

In this article, a wave propagation model is presented as the first step in the development of a new type of transluminal procedure for performing elastography. Elastography is a medical imaging modality for mapping the elastic properties of soft tissue. The wave propagation model is based on a Kelvin Voigt Fractional Derivative (KVFD) viscoelastic wave equation, and is numerically solved using a Finite Difference Time Domain (FDTD) method. Fractional rheological models, such as the KVFD, are particularly well suited to model the viscoelastic response of soft tissue in elastography. The transluminal procedure is based on the transmission and detection of shear waves through the luminal wall. Shear waves travelling through the tissue are perturbed after encountering areas of altered elasticity. These perturbations carry information of medical interest that can be extracted by solving the inverse problem. Scattering from prostate tumours is used as an example application to test the model. In silico results demonstrate that shear waves are satisfactorily transmitted through the luminal wall and that echoes, coming from reflected energy at the edges of an area of altered elasticity, which are feasibly detectable by using the transluminal approach. The model here presented provides a useful tool to establish the feasibility of transluminal procedures based on wave propagation and its interaction with the mechanical properties of the tissue outside the lumen.

Why Are Viscosity and Nonlinearity Bound to Make an Impact in Clinical Elastographic Diagnosis?

The adoption of multiscale approaches by the biomechanical community has caused a major improvement in quality in the mechanical characterization of soft tissues. The recent developments in elastography techniques are enabling in vivo and non-invasive quantification of tissues' mechanical properties. Elastic changes in a tissue are associated with a broad spectrum of pathologies, which stems from the tissue microstructure, histology and biochemistry. This knowledge is combined with research evidence to provide a powerful diagnostic range of highly prevalent pathologies, from birth and labor disorders (prematurity, induction failures, etc.), to solid tumors (e.g., prostate, cervix, breast, melanoma) and liver fibrosis, just to name a few. This review aims to elucidate the potential of viscous and nonlinear elastic parameters as conceivable diagnostic mechanical biomarkers. First, by providing an insight into the classic role of soft tissue microstructure in linear elasticity; secondly, by understanding how viscosity and nonlinearity could enhance the current diagnosis in elastography; and finally, by compounding preliminary investigations of those elastography parameters within different technologies. In conclusion, evidence of the diagnostic capability of elastic parameters beyond linear stiffness is gaining momentum as a result of the technological and imaging developments in the field of biomechanics.

Shear Waves Reveal Viscoelastic Changes in Skeletal Muscles After Hemispheric Stroke

We investigated alterations in material properties such as elasticity and viscoelasticity of stroke-affected muscles using ultrasound induced shear waves and mechanical models. We used acoustic radiation force to generate shear waves along fascicles of biceps muscles and measured their propagation velocity. The shear wave data were collected in muscles of 13 hemiplegic stroke survivors under passive conditions at 90°, 120°, and 150° elbow flexion angles. In a viscoelastic medium, as opposed to a purely elastic medium, the shear wave propagation velocity depends on the frequency content of the induced wave. Therefore, in addition to the shear wave group velocity (GpV), we also estimated a frequency-dependent phase velocity (PhV). We found significantly higher GpVs and PhVs in stroke-affected muscles ( ). The velocity data were used to estimate shear elasticity and viscosity using an elastic and viscoelastic material models. A pure elastic model showed increased shear elasticity in stroke-affected muscles ( ). The Voigt model estimates of viscoelastic properties were also significantly different between the stroke-impaired and non-impaired muscles. We observed significantly larger model-estimated viscosity values on the stroke-affected side at elbow flexion angles of 120° and 150°. Furthermore, the creep behavior (tissue strain resulting from the application of sudden constant stress) of the model was also different between muscles of the paretic and non-paretic side. We speculate that these changes are associated with the structural disruption of muscles after stroke and may potentially affect force generation from muscle fibers as well as transmission of force to tendons.

Kelvin-Voigt Parameters Reconstruction of Cervical Tissue-Mimicking Phantoms Using Torsional Wave Elastography

The reconstruction of viscous properties of soft tissues, and more specifically, of cervical tissue is a challenging problem. In this paper, a new method is proposed to reconstruct the viscoelastic parameters of cervical tissue-mimicking phantoms by a Torsional Wave Elastography (TWE) technique. The reconstruction method, based on a Probabilistic Inverse Problem (PIP) approach, is presented and experimentally validated against Shear Wave Elastography (SWE). The anatomy of the cervical tissue has been mimicked by means of a two-layer gelatine phantom that simulates the epithelial and connective layers. Five ad hoc oil-in-gelatine phantoms were fabricated at different proportion to test the new reconstruction technique. The PIP approach was used for reconstructing the Kelvin-Voigt (KV) viscoelastic parameters by comparing the measurements obtained from the TWE technique with the synthetic signals from a Finite Difference Time Domain (FDTD) KV wave propagation model. Additionally, SWE tests were realized in order to characterize the viscoelastic properties of each batch of gelatine. Finally, validation was carried out by comparing the KV parameters inferred from the PIP with those reconstructed from the shear wave dispersion curve obtained from the SWE measurements. In order to test the degree of agreement between both techniques, a Student's T-test and a Pearson's correlation study were performed. The results indicate that the proposed method is able to reconstruct the KV viscoelastic properties of the cervical tissue, for both the epithelial and connective layers, as well as the thickness of the first layer with acceptable accuracy.

Feasibility study of using the dispersion of surface acoustic wave impulse for viscoelasticity characterization in tissue mimicking phantoms

The aim of this study was to investigate the feasibility of utilizing the phase velocity dispersion of impulse surface acoustic wave (SAW) for viscoelasticity characterization of soft materials. The focused ultrasound transducer and the phase-sensitive optical coherence tomography were applied as the impulse SAW inducer and tracker, respectively. Three types of liquid-paraffin-based cream-in-agar phantoms were tested. Phase velocity dispersion curve was extracted using a Fourier transform-based phase velocity analysis algorithm. Viscoelastic parameters were obtained by fitting the dispersion curve of SAW into Rayleigh wave dispersion equation. The estimated viscoelasticity was compared with that from spherical indenter, ramp-hold relaxation testing for validation. Both results show an increasing trend in the elasticity and decreasing trend in the viscosity with the concentration of liquid-paraffin-based cream increasing in the samples. The proposed method has the capability of evaluating the viscoelastic properties of homogeneous soft tissue. By combining viscoelastic parameters estimated from the proposed method, the dispersive SAW-impulse-based viscosity-compensated elastography could be further developed.

Novel application of three-dimensional shear wave elastography in the detection of clinically significant prostate cancer

The present study evaluated three-dimensional shear wave elastography (3D SWE) in the detection of clinically significant prostate cancer. Clinically significant prostate cancer was defined by a minimum of one biopsy core with a Gleason score of 3+4 or 6 with a maximum cancer core length >4 mm. Patients with serum prostate-specific antigen levels of 4.0-20.0 ng/ml who were suspected of having prostate cancer from multi-parametric magnetic resonance imaging (mpMRI) were prospectively recruited. The 3D SWE was performed pre-biopsy, after which patients underwent MRI-transrectal ultrasound image-guided targeted biopsies for cancer-suspicious lesions and 12-core systematic biopsies. The pathological biopsy results were compared with the mpMRI and 3D SWE images. A total of 12 patients who were suspected of having significant cancer on mpMRI were included. The median pre-biopsy PSA value was 5.65 ng/ml. Of the 12 patients, 10 patients were diagnosed as having prostate cancer. In the targeted biopsy lesions, there was a significant difference in Young's modulus between the cancer-detected area (median 64.1 kPa, n=20) and undetected area (median 30.8 kPa, n=8; P<0.0001). On evaluation of receiver operating characteristics, a cut-off value of the Young's modulus of 41.0 kPa was used for the detection of clinically significant cancer, with which the sensitivity, specificity, positive predictive value (PPV) and negative predictive value (NPV) of cancer detection were 58, 97, 86 and 87%, respectively. When combining this cut-off tissue elasticity value with Prostate Imaging Reporting and Data System (PI-RADS) scores, the sensitivity, specificity, positive predictive value and negative predictive value of cancer detection were improved to 70, 98, 91 and 92%, respectively. In the cancer-detected lesions, a significant correlation was identified between the tissue elasticity value of the lesions and Gleason score (r=0.898, P<0.0001). In conclusion, PI-RADS combined with measurement of Young's modulus by 3D SWE may improve the diagnosis of clinically significant prostate cancer.

Carotid Plaque Stiffness Measured with Supersonic Shear Imaging and Its Correlation with Serum Homocysteine Level in Ischemic Stroke Patients

Objective

To ascertain the feasibility of using shear wave velocity (SWV) in assessing the stiffness of carotid plaque by supersonic shear imaging (SSI) and explore preliminary clinical value for such evaluation.

Materials and Methods

Supersonic shear imaging was performed in 142 patients with ischemic stroke, including 76 males and 66 females with mean age of 66 years (range, 45-80 years). The maximum, minimum, and mean values of SWV were measured for 129 carotid plaques. SWVs were compared between echolucent and echogenic plaques. Correlations between SWVs and serum homocysteine levels were investigated. Based on neurological symptom, the surrogate marker of vulnerable plaque (VP), binary logistic regression was performed and area under curve (AUC) of homocysteine only and homocysteine combing SWV_mean was calculated respectively.

Results

Echogenic plaques (n = 51) had higher SWVs than echolucent ones (n = 78) (SWV_min 3.91 [3.24-4.17] m/s vs. 1.51 [1.04-1.94] m/s; SWV_mean, 4.29 [3.98-4.57] m/s vs. 2.09 [1.69-2.41] m/s; SWV_max, 4.67 [4.33-4.86] m/s vs. 2.62 [2.32-3.31] m/s all p values < 0.01). Pearson correlation analysis showed that stiffness of plaques was negatively correlated with homocysteine level. R values for SWV_min, SWV_mean, and SWV_max were -0.205, -0.213, and -0.199, respectively. Binary logistic regression analysis showed that sex (p = 0.008), low-density lipoprotein (p = 0.015), triglycerides (p = 0.011), SWV_mean (p = 0.004), and hyper-homocysteinemia (p = 0.010) were significantly associated with symptomatic ischemic stroke. Receiver operating characteristic curves revealed that SWV_mean combing serum homocysteine level (AUC = 0.67) presented better diagnostic value than serum homocysteine only (AUC = 0.60) for symptomatic ischemic stroke.

Conclusion

Supersonic shear imaging could be used to quantitatively evaluate stiffness of both echolucent and echogenic carotid plaques. More importantly, SWVs of plaques were not only correlated to serum homocysteine level, but also associated with symptomatic ischemic stroke, suggesting that SSI might be useful for understanding more about VP.

An ultrasound transient elastography system with coded excitation

BACKGROUND

Ultrasound transient elastography technology has found its place in elastography because it is safe and easy to operate. However, it's application in deep tissue is limited. The aim of this study is to design an ultrasound transient elastography system with coded excitation to obtain greater detection depth.

METHODS

The ultrasound transient elastography system requires tissue vibration to be strictly synchronous with ultrasound detection. Therefore, an ultrasound transient elastography system with coded excitation was designed. A central component of this transient elastography system was an arbitrary waveform generator with multi-channel signals output function. This arbitrary waveform generator was used to produce the tissue vibration signal, the ultrasound detection signal and the synchronous triggering signal of the radio frequency data acquisition system. The arbitrary waveform generator can produce different forms of vibration waveform to induce different shear wave propagation in the tissue. Moreover, it can achieve either traditional pulse-echo detection or a phase-modulated or a frequency-modulated coded excitation. A 7-chip Barker code and traditional pulse-echo detection were programmed on the designed ultrasound transient elastography system to detect the shear wave in the phantom excited by the mechanical vibrator. Then an elasticity QA phantom and sixteen in vitro rat livers were used for performance evaluation of the two detection pulses.

RESULTS

The elasticity QA phantom's results show that our system is effective, and the rat liver results show the detection depth can be increased more than 1 cm. In addition, the SNR (signal-to-noise ratio) is increased by 15 dB using the 7-chip Barker coded excitation.

CONCLUSIONS

Applying 7-chip Barker coded excitation technique to the ultrasound transient elastography can increase the detection depth and SNR. Using coded excitation technology to assess the human liver, especially in obese patients, may be a good choice.

Mechanics of ultrasound elastography

Ultrasound elastography enables in vivo measurement of the mechanical properties of living soft tissues in a non-destructive and non-invasive manner and has attracted considerable interest for clinical use in recent years. Continuum mechanics plays an essential role in understanding and improving ultrasound-based elastography methods and is the main focus of this review. In particular, the mechanics theories involved in both static and dynamic elastography methods are surveyed. They may help understand the challenges in and opportunities for the practical applications of various ultrasound elastography methods to characterize the linear elastic, viscoelastic, anisotropic elastic and hyperelastic properties of both bulk and thin-walled soft materials, especially the in vivo characterization of biological soft tissues.

Inability of shear-wave elastography to distinguish malignant from benign prostate tissue - a comparison of biopsy, whole-mount sectioning and shear-wave elastography

Aim

This study was designed to assess the possible usefulness of shear-wave elastography in differentiating between benign and malignant tissue in prostate neoplasia.

Patients and methods

A total of 120 prostate tissue samples were obtained from 10 patients treated by radical prostatectomy and investigated pre-operatively by ultrasound elastography followed by directed biopsy. After resection, whole-mount sectioning and histological examination was performed. The predictions based on shear-wave elastography were compared with biopsy and histological results.

Results

The comparison between the results of shear-wave elastography and those of biopsy was performed by receiver operating characteristic analysis, which suggested an optimum cut-off tissue elasticity value of 50 kPa, in agreement with earlier studies aimed at distinguishing between benign and malignant tissue. However, the diagnostic selectivity (and thus the diagnostic power) was poor (area under the curve 0.527, which hardly differs from the value of 0.500 that would correspond to a complete lack of predictive power); furthermore, application of this cut-off value to the samples led to a sensitivity of only 74% and a specificity of only 43%. An analogous comparison between the results of shear-wave elastography and those of whole-mount histology, which itself is more reliable than biopsy, gave an even poorer diagnostic selectivity (sensitivity of 62%, specificity of 35%). Meaningful association with Gleason score was not found for D’Amico risk groups (p = 0.35).

Conclusions

The (negative) findings of this investigation add to the dissonance among results of studies investigating the possible value of shear-wave elastography as a diagnostic tool to identify malignant neoplasia. There is a clear need for further research to elucidate the diversity of study results and to identify the usefulness, if any, of the method in question.

One-dimensional nonlinear elastodynamic models and their local conservation laws with applications to biological membranes

The framework of incompressible nonlinear hyperelasticity and viscoelasticity is applied to the derivation of one-dimensional models of nonlinear wave propagation in fiber-reinforced elastic solids. Equivalence transformations are used to simplify the resulting wave equations and to reduce the number of parameters. Local conservation laws and global conserved quantities of the models are systematically computed and discussed, along with other related mathematical properties. Sample numerical solutions are presented. The models considered in the paper are appropriate for the mathematical description of certain aspects of the behavior of biological membranes and similar structures.

Shear-wave elasticity imaging of a liver fibrosis mouse model using high-frequency ultrasound

The objective of this study was to develop a high-frequency imaging platform for evaluating liver fibrosis in mice based on shear-wave elasticity imaging (SWEI). Although SWEI has been used to diagnose hepatic fibrosis clinically, it is performed at relatively low frequencies (<20 MHz). For preclinical ultrasound imaging in small animals, a high-frequency (>30 MHz) single-element transducer with mechanical scanning is often used. In this study we developed a new SWEI system based on a 40-MHz single-element transducer for imaging and a separate 20-MHz excitation transducer for producing the radiation force and the associated shear waves. Liver fibrosis was induced in ten C57BL/6 (B6) mice using carbon tetrachloride; the other ten mice served as the control group. Synchronizing the excitation beam (i.e., the beam from the excitation transducer) and the detection beam sequence (i.e., the beam from the imaging transducer) allows this mechanical-scanning setup to analyze the shear-wave dispersion relation. The liver viscoelastic properties were determined in vivo by measuring the shear-wave dispersion curve followed by fitting to the Voigt model. The mice were then killed and the fibrosis stage was evaluated (from F0 to F4) based on the METAVIR score. The measured mean values of liver elasticity and viscosity, respectively, ranged from 1.06 to 1.89 kPa and from 1.29 to 1.75 Pa∙s for normal F0 and fibrosis stages of F3 and F4. The Spearman coefficients for the correlations between the measured elasticity and viscosity at various fibrosis stages as assessed by the METAVIR score were 0.73 (p < 0.001) and 0.634 (p = 0.0013), respectively. We also found that the collagen content in the liver was linearly correlated with the measured elasticity (r(2) = 0.54, p < 0.001) and less strongly with the viscosity (r2 = 0.26, p = 0.022). Finally, the diagnosis performance of high-frequency SWEI was evaluated using multivariate receiver operating characteristic curve (ROC) analysis. The areas under the multivariate ROC curve for diagnosing fibrosis stages of F ≥ 3, F = 4, F0 vs. F3, F0 vs. F4, and F3 vs. F4 were 0.9, 0.98, 0.83, 1.0, and 0.96, respectively. Compared with traditional ROC analysis, an improved diagnosis performance was found for diagnosing fibrosis stages of F ≥ 3 and F0 vs. F3. These results demonstrate that the developed high-frequency SWEI platform can yield quantitative viscoelastic properties for diagnosing various fibrosis stages in mice. It is a promising tool for studying the progression of liver fibrosis in preclinical animal models both noninvasively and quantitatively.

Analyzing and modeling rheological behavior of liver fibrosis in rats using shear viscoelastic moduli

The process of liver fibrosis changes the rheological properties of liver tissue. This study characterizes and compares liver fibrosis stages from F0 to F4 in rats in terms of shear viscoelastic moduli. Here two viscoelastic models, the Zener model and Voigt model, were applied to experimental data of rheometer tests and then values of elasticity and viscosity were estimated for each fibrosis stage. The results demonstrate that moderate fibrosis (≤F2) has a good correlation with liver viscoelasticity. The mean Zener elasticity E1 increases from (0.452±0.094) kPa (F0) to (1.311±0.717) kPa (F2), while the mean Voigt elasticity E increases from (0.618±0.089) kPa (F0) to (1.701±0.844) kPa (F2). The mean Zener viscosity increases from (3.499±0.186) Pa·s (F0) to (4.947±1.811) Pa·s (F2) and the mean Voigt viscosity increases from (3.379±0.316) Pa·s (F0) to (4.625±1.296) Pa·s (F2). Compared with viscosity, the elasticity shows smaller variations at stages F1 and F2 no matter what viscoelastic model is used. Therefore, the estimated elasticity is more effective than viscosity for differentiating the fibrosis stages from F0 to F2.

New aspects in shear-wave elastography of prostate cancer

AIM

This study was designed to evaluate the performance of shear-wave elastography as a diagnostic tool for prostate cancer in a larger cohort of patients than previously reported.

PATIENTS AND METHODS

Seventy-three patients with suspected prostate carcinoma were investigated by ultrasound elastography followed by directed biopsy. The elastographic and histological results for all biopsies were compared.

RESULTS

After exclusion of invalid and non-assessable results, 794 samples were obtained for which both a histological assessment and an elastometric result (tissue stiffness in kPa) were available: according to the histology 589 were benign and 205 were malignant. Tissue elasticity was found to be weakly correlated with patient's age, PSA level and gland volume. ROC analysis showed that, for the set of results acquired, elastometry did not fulfil literature claims that it could identify malignant neoplasia with high sensitivity and specificity. However, it did show promise in distinguishing between Gleason scores ≤6 and >6 when malignancy had already been identified. Unexpected observations were the finding of a smaller proportion of tumours in the lateral regions of the prostate than generally expected, and also the observation that the elasticity of benign prostate tissue is region-sensitive, the tissue being stiffest in the basal region and more elastic at the apex.

CONCLUSIONS

Shear-wave elastography was found to be a poor predictor of malignancy, but for malignant lesions an elasticity cut-off of 80 kPa allowed a fairly reliable distinction between lesions with Gleason ≤6 and those with Gleason >6. We demonstrate an increase in elasticity of benign prostate tissue from the basal to the apical region.

Multiparametric 3D in vivo ultrasound vibroelastography imaging of prostate cancer: Preliminary results

PURPOSE

Ultrasound-based solutions for diagnosis and prognosis of prostate cancer are highly desirable. The authors have devised a method for detecting prostate cancer using a vibroelastography (VE) system developed in our group and a tissue classification approach based on texture analysis of VE images.

METHODS

The VE method applies wide-band mechanical vibrations to the tissue. Here, the authors report on the use of this system for cancer detection and show that the texture of VE images characterized by the first and the second order statistics of the pixel intensities form a promising set of features for tissue typing to detect prostate cancer. The system was used to image patients prior to radical surgery. The removed specimens were sectioned and studied by an experienced histopathologist. The authors registered the whole-mount histology sections to the ultrasound images using an automatic registration algorithm. This enabled the quantitative evaluation of the performance of the authors' imaging method in cancer detection in an unbiased manner. The authors used support vector machine (SVM) classification to measure the cancer detection performance of the VE method. Regions of tissue of size 5 × 5 mm, labeled as cancer and noncancer based on automatic registration to histology slides, were classified using SVM.

RESULTS

The authors report an area under ROC of 0.81 ± 0.10 in cancer detection on 1066 tissue regions from 203 images. All cancer tumors in all zones were included in this analysis and were classified versus the noncancer tissue in the peripheral zone. This outcome was obtained in leave-one-patient-out validation.

CONCLUSIONS

The developed 3D prostate vibroelastography system and the proposed multiparametric approach based on statistical texture parameters from the VE images result in a promising cancer detection method.

The role of viscosity estimation for oil-in-gelatin phantom in shear wave based ultrasound elastography

Shear wave based ultrasound elastography utilizes mechanical excitation or acoustic radiation force to induce shear waves in deep tissue. The tissue response is monitored to obtain elasticity information about the tissue. During the past two decades, tissue elasticity has been extensively studied and has been used in clinical disease diagnosis. However, biological soft tissues are viscoelastic in nature. Therefore, they should be simultaneously characterized in terms of elasticity and viscosity. In this study, two shear wave-based elasticity imaging methods, shear wave dispersion ultrasound vibrometry (SDUV) and acoustic radiation force impulsive (ARFI) imaging, were compared. The discrepancy between the measurements obtained by the two methods was analyzed, and the role of viscosity was investigated. To this end, four types of gelatin phantoms containing 0%, 20%, 30% and 40% castor oil were fabricated to mimic different viscosities of soft tissue. For the SDUV method, the shear elasticity μ1 was 3.90 ± 0.27 kPa, 4.49 ± 0.16 kPa, 2.41 ± 0.33 kPa and 1.31 ± 0.09 kPa; and the shear viscosity μ2 was 1.82 ± 0.31 Pa•s, 2.41 ± 0.35 Pa•s, 2.65 ± 0.13 Pa•s and 2.89 ± 0.14 Pa•s for 0%, 20%, 30% and 40% oil, respectively in both cases. For the ARFI measurements, the shear elasticity μ was 7.30 ± 0.20 kPa, 8.20 ± 0.31 kPa, 7.42 ± 0.21 kPa and 5.90 ± 0.36 kPa for 0%, 20%, 30% and 40% oil, respectively. The SDUV results demonstrated that the elasticity first increased from 0% to 20% oil and then decreased for the 30% and 40% oil. The viscosity decreased consistently as the concentration of castor oil increased from 0% to 40%. The elasticity measured by ARFI showed the same trend as that of the SDUV but exceeded the results measured by SDUV. To clearly validate the impact of viscosity on the elasticity estimation, an independent measurement of the elasticity and viscosity by dynamic mechanical analysis (DMA) was conducted on these four types of gelatin phantoms and then compared with SDUV and ARFI results. The shear elasticities obtained by DMA (3.44 ± 0.31 kPa, 4.29 ± 0.13 kPa, 2.05 ± 0.29 kPa and 1.06 ± 0.18 kPa for 0%, 20%, 30% and 40% oil, respectively) were lower than those by SDUV, whereas the shear viscosities obtained by DMA (2.52 ± 0.32 Pa·s, 3.18 ± 0.12 Pa·s, 3.98 ± 0.19 Pa·s and 4.90 ± 0.20 Pa·s for 0%, 20%, 30% and 40% oil, respectively) were greater than those obtained by SDUV. However, the DMA results showed that the trend in the elasticity and viscosity data was the same as that obtained from the SDUV and ARFI. The SDUV results demonstrated that adding castor oil changed the viscoelastic properties of the phantoms and resulted in increased dispersion of the shear waves. Viscosity can provide important and independent information about the inner state of the phantoms, in addition to the elasticity. Because the ARFI method ignores the dispersion of the shear waves, namely viscosity, it may bias the estimation of the true elasticity. This study sheds further light on the significance of the viscosity measurements in shear wave based elasticity imaging methods.

Quantitative analysis of liver fibrosis in rats with shearwave dispersion ultrasound vibrometry: comparison with dynamic mechanical analysis

Ultrasonic elastography, a non-invasive technique for assessing the elasticity properties of tissues, has shown promising results for disease diagnosis. However, biological soft tissues are viscoelastic in nature. Shearwave dispersion ultrasound vibrometry (SDUV) can simultaneously measure the elasticity and viscosity of tissue using shear wave propagation speeds at different frequencies. In this paper, the viscoelasticity of rat livers was measured quantitatively by SDUV for normal (stage F0) and fibrotic livers (stage F2). Meanwhile, an independent validation study was presented in which SDUV results were compared with those derived from dynamic mechanical analysis (DMA), which is the only mechanical test that simultaneously assesses the viscoelastic properties of tissue. Shear wave speeds were measured at frequencies of 100, 200, 300 and 400 Hz with SDUV and the storage moduli and loss moduli were measured at the frequency range of 1-40 Hz with DMA. The Voigt viscoelastic model was used in the two methods. The mean elasticity and viscosity obtained by SDUV ranged from 0.84±0.13 kPa (F0) to 1.85±0.30 kPa (F2) and from 1.12±0.11 Pa s (F0) to 1.70±0.31 Pa s (F2), respectively. The mean elasticity and viscosity derived from DMA ranged from 0.62±0.09 kPa (F0) to 1.70±0.84 kPa (F2) and from 3.38±0.32 Pa s (F0) to 4.63±1.30 Pa s (F2), respectively. Both SDUV and DMA demonstrated that the elasticity of rat livers increased from stage F0 to F2, a finding which was consistent with previous literature. However, the elasticity measurements obtained by SDUV had smaller differences than those obtained by DMA, whereas the viscosities obtained by the two methods were obviously different. We suggest that the difference could be related to factors such as tissue microstructure, the frequency range, sample size and the rheological model employed. For future work we propose some improvements in the comparative tests between SDUV and DMA, such as enlarging the harmonic frequency range of the shear wave to highlight the role of viscosity, finding an appropriate rheological model to improve the accuracy of tissue viscoelasticity estimations.

Theoretical prediction of ultrasound elastography for detection of early osteoarthritis

Ultrasound elastography could be used as a new noninvasive technique for detecting early osteoarthritis. As the first critical step, this study theoretically predicted the excitation power and the measurement errors in detecting cartilage detect. A finite element model was used to simulate wave propagation of elastography in the cartilage. The wave was produced by a force F, and the wave speed C was calculated. The normal cartilage model was used to define the relationship between the wave speed and elastic modulus. Various stiffness values were simulated. F = 10 N with a duration of 0.5 ms was required for having measurable deformation (10 μm) at the distal site. The deformation had a significant rise when the wave crossed the defect. The relationship between the wave speed and elastic parameters was found as C = 1.57 × (E)/(2 × ρ(1+μ)))^1/2, where E was the elastic modulus, μ was Poisson's ratio, and ρ was the density. For the simulated defect with an elastic modulus of 7 MPa which was slightly stiffer than the normal cartilage, the measurement error was 0.1 MPa. The results suggested that, given the simulated conditions, this new technique could be used to detect the defect in early osteoarthritis.

Ultrasound vibrometry using orthogonal- frequency-based vibration pulses

New vibration pulses are developed for shear wave generation in a tissue region with preferred spectral distributions for ultrasound vibrometry applications. The primary objective of this work is to increase the frequency range of detectable harmonics of the shear wave. The secondary objective is to reduce the required peak intensity of transmitted pulses that induce the vibrations and shear waves. Unlike the periodic binary vibration pulses, the new vibration pulses have multiple pulses in one fundamental period of the vibration. The pulses are generated from an orthogonal-frequency wave composed of several sinusoidal signals, the amplitudes of which increase with frequency to compensate for higher loss at higher frequency in tissues. The new method has been evaluated by studying the shear wave propagation in in vitro chicken and swine liver. The experimental results show that the new vibration pulses significantly increase tissue vibration with a reduced peak ultrasound intensity, compared with the binary vibration pulses.

Multi-push (MP) acoustic radiation force (ARF) ultrasound for assessing tissue viscoelasticity, in vivo

Acoustic radiation force (ARF) ultrasound is a method of elastographic imaging in which micron-scale tissue displacements, induced and tracked by ultrasound, reflect clinically relevant tissue mechanical properties. Our laboratory has recently shown that tissue viscoelasticity is assessed using the novel Multi-Push (MP) ARF method. MP ARF applies the Voigt model for viscoelastic materials and compares the displacements achieved by successive ARF excitations to qualitatively or quantitatively represent the relaxation time for constant stress, which is a direct descriptor of the viscoelastic response of the tissue. We have demonstrated MP ARF in custom viscoelastic tissue mimicking materials and implemented the method in vivo in canine muscle and human renal allografts, with strong spatial correlation between MP ARF findings and histochemical features and previously reported mechanical changes with renal disease. These data support that noninvasive MP ARF is capable of clinically relevant assessment of tissue viscoelastic properties.

Acoustic waves in medical imaging and diagnostics

Up until about two decades ago acoustic imaging and ultrasound imaging were synonymous. The term ultrasonography, or its abbreviated version sonography, meant an imaging modality based on the use of ultrasonic compressional bulk waves. Beginning in the 1990s, there started to emerge numerous acoustic imaging modalities based on the use of a different mode of acoustic wave: shear waves. Imaging with these waves was shown to provide very useful and very different information about the biological tissue being examined. We discuss the physical basis for the differences between these two basic modes of acoustic waves used in medical imaging and analyze the advantages associated with shear acoustic imaging. A comprehensive analysis of the range of acoustic wavelengths, velocities and frequencies that have been used in different imaging applications is presented. We discuss the potential for future shear wave imaging applications.

Acoustic radiation force elasticity imaging in diagnostic ultrasound

The development of ultrasound-based elasticity imaging methods has been the focus of intense research activity since the mid-1990s. In characterizing the mechanical properties of soft tissues, these techniques image an entirely new subset of tissue properties that cannot be derived with conventional ultrasound techniques. Clinically, tissue elasticity is known to be associated with pathological condition and with the ability to image these features in vivo; elasticity imaging methods may prove to be invaluable tools for the diagnosis and/or monitoring of disease. This review focuses on ultrasound-based elasticity imaging methods that generate an acoustic radiation force to induce tissue displacements. These methods can be performed noninvasively during routine exams to provide either qualitative or quantitative metrics of tissue elasticity. A brief overview of soft tissue mechanics relevant to elasticity imaging is provided, including a derivation of acoustic radiation force, and an overview of the various acoustic radiation force elasticity imaging methods.

Loss tangent and complex modulus estimated by acoustic radiation force creep and shear wave dispersion

Elasticity imaging methods have been used to study tissue mechanical properties and have demonstrated that tissue elasticity changes with disease state. In current shear wave elasticity imaging methods typically only shear wave speed is measured and rheological models, e.g. Kelvin-Voigt, Maxwell and Standard Linear Solid, are used to solve for tissue mechanical properties such as the shear viscoelastic complex modulus. This paper presents a method to quantify viscoelastic material properties in a model-independent way by estimating the complex shear elastic modulus over a wide frequency range using time-dependent creep response induced by acoustic radiation force. This radiation force induced creep method uses a conversion formula that is the analytic solution of a constitutive equation. The proposed method in combination with shearwave dispersion ultrasound vibrometry is used to measure the complex modulus so that knowledge of the applied radiation force magnitude is not necessary. The conversion formula is shown to be sensitive to sampling frequency and the first reliable measure in time according to numerical simulations using the Kelvin-Voigt model creep strain and compliance. Representative model-free shear complex moduli from homogeneous tissue mimicking phantoms and one excised swine kidney were obtained. This work proposes a novel model-free ultrasound-based elasticity method that does not require a rheological model with associated fitting requirements.

A Review of Shearwave Dispersion Ultrasound Vibrometry (SDUV) and its Applications

Measurement of tissue elasticity has emerged as an important advance in medical imaging and tissue characterization. However, soft tissue is inherently a viscoelastic material. One way to characterize the viscoelastic material properties of a material is to measure shear wave propagation velocities within the material at different frequencies and use the dispersion of the velocities, or variation with frequency, to solve for the material properties. Shearwave Dispersion Ultrasound Vibrometry (SDUV) is an ultrasound-based technique that uses this feature to characterize the viscoelastic nature of soft tissue. This method has been used to measure the shear elasticity and viscosity in various types of soft tissues including skeletal muscle, cardiac muscle, liver, kidney, prostate, and arterial vessels. This versatile technique provides measurements of viscoelastic material properties with high spatial and temporal resolution, which can be used for assessing these properties in normal and pathologic tissues. The goals of this paper are to 1) give an overview of viscoelasticity and shear wave velocity dispersion, 2) provide a history of the development of the SDUV method, and 3) survey applications for SDUV that have been previously reported.

AN OVERVIEW OF ELASTOGRAPHY - AN EMERGING BRANCH OF MEDICAL IMAGING

From times immemorial manual palpation served as a source of information on the state of soft tissues and allowed detection of various diseases accompanied by changes in tissue elasticity. During the last two decades, the ancient art of palpation gained new life due to numerous emerging elasticity imaging (EI) methods. Areas of applications of EI in medical diagnostics and treatment monitoring are steadily expanding. Elasticity imaging methods are emerging as commercial applications, a true testament to the progress and importance of the field.In this paper we present a brief history and theoretical basis of EI, describe various techniques of EI and, analyze their advantages and limitations, and overview main clinical applications. We present a classification of elasticity measurement and imaging techniques based on the methods used for generating a stress in the tissue (external mechanical force, internal ultrasound radiation force, or an internal endogenous force), and measurement of the tissue response. The measurement method can be performed using differing physical principles including magnetic resonance imaging (MRI), ultrasound imaging, X-ray imaging, optical and acoustic signals.Until recently, EI was largely a research method used by a few select institutions having the special equipment needed to perform the studies. Since 2005 however, increasing numbers of mainstream manufacturers have added EI to their ultrasound systems so that today the majority of manufacturers offer some sort of Elastography or tissue stiffness imaging on their clinical systems. Now it is safe to say that some sort of elasticity imaging may be performed on virtually all types of focal and diffuse disease. Most of the new applications are still in the early stages of research, but a few are becoming common applications in clinical practice.

Assessing the viscoelastic properties of thrombus using a solid-sphere-based instantaneous force approach

The viscoelastic properties of thrombus play a significant role when the clot closes a leak in a vessel of the blood circulation. The common method used to measure the viscoelastic properties of a clot employs a rheometer but this might be unsuitable due to the clot fiber network being broken up by excessive deformation. This study assessed the feasibility of using a novel acoustic method to assess the viscoelastic properties of blood clots. This method is based on monitoring the motion of a solid sphere in a blood clot induced by an applied instantaneous force. Experiments were performed in which a solid sphere was displaced by a 1 MHz single-element focused transducer, with a 20 MHz single-element focused transducer used to track this displacement. The spatiotemporal behavior of the sphere displacement was used to determine the viscoelastic properties of the clot. The experimental system was calibrated by measuring the viscoelastic modulus of gelatin using different types of solid spheres embedded in the phantoms and, then, the shear modulus and viscosity of porcine blood clots with hematocrits of 0% (plasma), 20% and 40% were assessed. The viscoelastic modulus of each clot sample was also measured directly by a rheometer for comparison. The results showed that the shear modulus increased from 173 ± 52 (mean ± SD) Pa for 40%-hematocrit blood clots to 619.5 ± 80.5 Pa for plasma blood clots, while the viscosity decreased from 0.32 ± 0.07 Pa∙s to 0.16 ± 0.06 Pa∙s, respectively, which indicated that the concentration of red blood cells and the amount of fibrinogen are the main determinants of the clot viscoelastic properties.

On the effects of reflected waves in transient shear wave elastography

In recent years, novel quantitative techniques have been developed to provide noninvasive and quantitative stiffness images based on shear wave propagation. Using radiation force and ultrafast ultrasound imaging, the supersonic shear imaging technique allows one to remotely generate and follow a transient plane shear wave propagating in vivo in real time. The tissue shear modulus, i.e., its stiffness, can then be estimated from the shear wave local velocity. However, because the local shear wave velocity is estimated using a time-of- flight approach, reflected shear waves can cause artifacts in the estimated shear velocity because the incident and reflected waves propagate in opposite directions. Such effects have been reported in the literature as a potential drawback of elastography techniques based on shear wave speed, particularly in the case of high stiffness contrasts, such as in atherosclerotic plaque or stiff lesions. In this letter, we present our implementation of a simple directional filter, previously used for magnetic resonance elastography, which separates the forward- and backward-propagating waves to solve this problem. Such a directional filter could be applied to many elastography techniques based on the local estimation of shear wave speed propagation, such as acoustic radiation force imaging (ARFI), shearwave dispersion ultrasound vibrometry (SDUV), needle-based elastography, harmonic motion imaging, or crawling waves when the local propagation direction is known and high-resolution spatial and temporal data are acquired.

Instantaneous axial force of a high-order Bessel vortex beam of acoustic waves incident upon a rigid movable sphere

The present investigation examines the instantaneous force resulting from the interaction of an acoustical high-order Bessel vortex beam (HOBVB) with a rigid sphere. The rigid sphere case is important in fluid dynamics applications because it perfectly simulates the interaction of instantaneous sound waves in a reduced gravity environment with a levitated spherical liquid soft drop in air. Here, a closed-form solution for the instantaneous force involving the total pressure field as well as the Bessel beam parameters is obtained for the case of progressive, stationary and quasi-stationary waves. Instantaneous force examples for progressive waves are computed for both a fixed and a movable rigid sphere. The results show how the instantaneous force per unit cross-sectional surface and unit pressure varies versus the dimensionless frequency ka (k is the wave number in the fluid medium and a is the sphere's radius), the half-cone angle β and the order m of the HOBVB. It is demonstrated here that the instantaneous force is determined only for (m,n) = (0,1) (where n is the partial-wave number), and vanishes for m>0 because of symmetry. In addition, the instantaneous force and normalized amplitude velocity results are computed and compared with those of a rigid immovable (fixed) sphere. It is shown that they differ significantly for ka values below 5. The proposed analysis may be of interest in the analysis of instantaneous forces on spherical particles for particle manipulation, filtering, trapping and drug delivery. The presented solutions may also serve as a method for comparison to other solutions obtained by strictly numerical or asymptotic approaches.