Viscoelastic imaging of breast tumor microenvironment with ultrasound

Amplitude-Modulation Frequency Optimization for Enhancing Harmonic Motion Imaging Performance of Breast Tumors in the Clinic

OBJECTIVE

Elastography images tissue mechanical responses and infers the underlying properties to aid diagnosis and treatment response monitoring. The estimation of absolute or relative tumor properties may vary with dimensions even when the mechanical properties remain constant. Harmonic motion imaging (HMI) uses amplitude-modulated (AM) focused ultrasound to interrogate the targeted tissue's viscoelastic properties. In this study, effects of AM frequencies on HMI were investigated in terms of inclusion relative stiffness and size estimation.

METHODS

AM frequencies from 200 to 600 Hz in steps of 100 Hz were considered using a 5.3-kPa phantom with cylindrical inclusions (Young's modulus: 22, 31, 44, 56 kPa, and diameter: 4.8, 8.1, 13.6, 19.8 mm) to optimize the performance of HMI in characterizing tumors with the same mechanical properties and of different dimensions.

RESULTS

Consistent displacement ratios (DRs) (17.5% variation) of the inclusion to background were obtained with 200-Hz AM for breast-tumor-mimicking inclusions albeit a suboptimal inclusion size estimation obtained. 400-Hz was otherwise used for small and low-contrast inclusions (4.8 mm, 22 or 31 kPa). A linear relationship (R2 = 0.9043) was found between the inverse DR at these frequencies and the Young's modulus ratio. 400 Hz obtained the most accurate inclusion size estimation with an overall estimation error on the lateral dimension of 0.5 mm. In vivo imaging of breast cancer patients (n = 5) was performed at 200 or 400 Hz.

CONCLUSION

The results presented herein indicate that the HMI AM frequency could be optimized adaptively in cases of different applications, i.e., at 200 or 400 Hz, depending on whether aimed for consistent DR measurement for tumor response assessment or tumor margin delineation for surgical planning. HMI may thus be capable of predicting the pathologic endpoint of tumors in response to neoadjuvant chemotherapy (NACT) as early as 3 weeks into treatment.

Mathematical Models for Ultrasound Elastography: Recent Advances to Improve Accuracy and Clinical Utility

Changes in biomechanical properties such as elasticity modulus, viscosity, and poroelastic features are linked to the health status of biological tissues. Ultrasound elastography is a non-invasive imaging tool that quantitatively maps these biomechanical characteristics for diagnostic and treatment monitoring purposes. Mathematical models are essential in ultrasound elastography as they convert the raw data obtained from tissue displacement caused by ultrasound waves into the images observed by clinicians. This article reviews the available mathematical frameworks of continuum mechanics for extracting the biomechanical characteristics of biological tissues in ultrasound elastography. Continuum-mechanics-based approaches such as classical viscoelasticity, elasticity, and poroelasticity models, as well as nonlocal continuum-based models, are described. The accuracy of ultrasound elastography can be increased with the recent advancements in continuum modelling techniques including hyperelasticity, biphasic theory, nonlocal viscoelasticity, inversion-based elasticity, and incorporating scale effects. However, the time taken to convert the data into clinical images increases with more complex models, and this is a major challenge for expanding the clinical utility of ultrasound elastography. As we strive to provide the most accurate imaging for patients, further research is needed to refine mathematical models for incorporation into the clinical workflow.

Shear wave elastography can stratify rectal cancer response to short-course radiation therapy

Rectal cancer is a deadly disease typically treated using neoadjuvant chemoradiotherapy followed by total mesorectal excision surgery. To reduce the occurrence of mesorectal excision surgery for patients whose tumors regress from the neoadjuvant therapy alone, conventional imaging, such as computed tomography (CT) or magnetic resonance imaging (MRI), is used to assess tumor response to neoadjuvant therapy before surgery. In this work, we hypothesize that shear wave elastography offers valuable insights into tumor response to short-course radiation therapy (SCRT)-information that could help distinguish radiation-responsive from radiation-non-responsive tumors and shed light on changes in the tumor microenvironment that may affect radiation response. To test this hypothesis, we performed elastographic imaging on murine rectal tumors (n = 32) on days 6, 10, 12, 16, 18, 20, 23, and 25 post-tumor cell injection. The study revealed that radiation-responsive and non-radiation-responsive tumors had different mechanical properties. Specifically, radiation-non-responsive tumors showed significantly higher shear wave speed SWS (p < 0.01) than radiation-responsive tumors 11 days after SCRT. Furthermore, there was a significant difference in shear wave attenuation (SWA) (p < 0.01) in radiation-non-responsive tumors 16 days after SCRT compared to SWA measured just one day after SCRT. These results demonstrate the potential of shear wave elastography to provide valuable insights into tumor response to SCRT and aid in exploring the underlying biology that drives tumors' responses to radiation.

The emerging promise of tumour mechanobiology in cancer treatment

Tumour cell biomechanics has lately came to the fore as a disparate feature that fosters cancer development and progression. Tumour mechanosensing entails a mechanical interplay amongst tumour cells, extracellular matrix (ECM) and cells of the tumour microenvironment (TME). Sensory receptors (mechanoceptors) detect changes of extracellular mechanical inputs such as various types of mechanical forces/stress and trigger oncogenic signalling pathways advocating for cancer initiation, growth, survival, angiogenesis, invasion, metastasis, and immune evasion. Moreover, alterations in ECM stiffness and potentiation of mechanostimulated transcriptional regulatory molecules (transcription factors/cofactors) have been shown to strongly correlate with resistance to anticancer drugs. On this basis, new mechanosensitive proteins emerge as potential therapeutic targets and/or biomarkers in cancer. Accordingly, tumour mechanobiology arises as a promising field that can potentially provide novel combinatorial regimens to reverse drug resistance, as well as offer unprecedented targeting approaches that may help to more effectively treat a large proportion of solid tumours and their complications. Here, we highlight recent findings regarding various aspects of tumour mechanobiology in the clinical setting and discuss evidence-based perspectives of developing diagnostic/prognostic tools and therapeutic approaches that exploit tumour-TME physical associations.

Extracellular matrix as a driver for intratumoral heterogeneity

The architecture of an organ is built through interactions between its native cells and its connective tissue consisting of stromal cells and the extracellular matrix (ECM). Upon transformation through tumorigenesis, such interactions are disrupted and replaced by a new set of intercommunications between malignantly transformed parenchyma, an altered stromal cell population, and a remodeled ECM. In this perspective, we propose that the intratumoral heterogeneity of cancer cell phenotypes is an emergent property of such reciprocal intercommunications, both biochemical and mechanical-physical, which engender and amplify the diversity of cell behavioral traits. An attempt to assimilate such findings within a framework of phenotypic plasticity furthers our understanding of cancer progression.

Predictive value of combining clinicopathological, multimodal ultrasonic characteristics in axillary lymph nodal metastasis burden of patients with cT1-2N0 breast cancer

OBJECTIVE

Since the ACOSOG Z0011 trial, the clinical examination of axillary lymph node-negative early breast cancer patients (cT1-2N0) can be used to predict the burden of axillary lymph nodes (ALNs) by axillary ultrasound (AUS). To improving diagnosis of axillary lymph node metastasis (ALNM), we assessed the value of combining clinicopathological, conventional ultrasound, SWE features in the cT1-2N0 breast cancer patients.

METHODS

Retrospective analysis of 285 patients with cT1-2N0 breast cancer who underwent preoperative ultrasound examination of the lesion and axillary, with shear wave elastography (SWE) of the lesions. According to the postoperative pathological results, they were divided into ≤2 metastatic ALNs group (low nodal burden, LNB) and >  2 metastatic ALNs group (high nodal burden, HNB). Binary logistic regression analysis was used to screen independent risk factors and establish prediction models. The best cut-off value of continuous variables is determined by the receiver operating characteristic curve, and the performance of the prediction model is evaluated.

RESULTS

Presence of lymphovascular invasion (OR = 7.966, P = 0.010), tumor size (OR = 2.485, P = 0.019), Emean of intratumor (OR = 0.939, P = 0.002) and cortical thickness of lymph node (OR = 9.277, P <  0.001) were independent risk predictors for HNB of cT1-2N0 Group. The predictive model of combined method had better performance in predicting HNB of cT1-2N0 compared with models based on SWE and conventional ultrasound alone (area under the curve: 0.824 vs 0.658, P <  0.001; 0.824 vs 0.789, P = 0.035).

CONCLUSIONS

The predictive models of combined method obtained from significant clinicopathological and ultrasonographic features can potentially improve the diagnosis and individual treatment of ALNM in patients with cT1-2N0 breast cancer.

Imaging the Local Nonlinear Viscoelastic Properties of Soft Tissues: Initial Validation and Expected Benefits

Imaging tissue mechanical properties has shown promise in noninvasive assessment of numerous pathologies. Researchers have successfully measured many linear tissue mechanical properties in laboratory and clinical settings. Currently, multiple complex mechanical effects such as frequency-dependence, anisotropy, and nonlinearity are being investigated separately. However, a concurrent assessment of these complex effects may enable more complete characterization of tissue biomechanics and offer improved diagnostic sensitivity. In this work, we report for the first time a method to map the frequency-dependent nonlinear parameters of soft tissues on a local scale. We recently developed a nonlinear elastography model that combines strain measurements from arbitrary tissue compression with radiation-force-based broadband shear wave speed (WS) measurements. Here, we extended this model to incorporate local measurements of frequency-dependent shear modulus. This combined approach provides a local frequency-dependent nonlinear parameter that can be obtained with arbitrary, clinically realizable tissue compression. Initial assessments using simulations and phantoms validate the accuracy of this approach. We also observed improved contrast in nonlinearity parameter at higher frequencies. Results from ex-vivo liver experiments show 32, 25, 34, and 38 dB higher contrast in elastograms than traditional linear elasticity, elastic nonlinearity, viscosity, and strain imaging methods, respectively. A lesion, artificially created by injection of glutaraldehyde into a liver specimen, showed a 59% increase in the frequency-dependent nonlinear parameter and a 17% increase in contrast ratio.

Shear wave elastography and pulsed doppler for breast lesions: Similar diagnostic performance and positively correlated stiffness and blood flow resistance

PURPOSE

To compare the diagnostic performance of shear wave elastography (SWE) and pulsed Doppler ultrasound in breast lesions, and to explore whether the quantitative SWE parameters correlated with pulsed Doppler ultrasound parameters.

MATERIALS AND METHODS

Seventy-nine patients with 79 breast lesions who had undergone conventional ultrasound, pulsed Doppler ultrasound and SWE examination were included. All of them underwent core needle biopsy or surgery within one week. Parameters including Emax (the maximum elastic modulus), Emean (mean elastic modulus), Emin (minimum elastic modulus), Esd (elastic modulus standard deviation), and RI (resistive index), PI (pulsatility index), PSV (peak systolic velocity) and EDV (end diastolic velocity) were obtained for statistical analysis.

RESULTS

Almost all SWE parameters were significantly different between benign and malignant breast lesions (P＜0.05). There was no significant difference between Esd and PI (P＞0.05), which had the best AUC among SWE and vascular parameters respectively (0.877 vs. 0.871). Emax showed a moderate correlation with PI (P = 0.000, r = 0.552) and RI (P = 0.000, r = 0.544), and Esd moderately correlated with PI (P = 0.000, r = 0.567) and RI (P = 0.000, r = 0.546). For the benign group, no parameters showed any significant correlation (P＞0.05), while for the malignant group, Emax and Esd also significantly correlated with PI or RI.

CONCLUSIONS

SWE and pulsed Doppler ultrasound had similar diagnostic efficacy for breast lesions. SWE and pulsed Doppler parameters were significantly correlated in breast lesions, especially in malignant ones, indicating the potential association between elastographic and vascular characteristics of breast tumors.

Shear wave imaging and classification using extended Kalman filter and decision tree algorithm

Shear wave ultrasound elastography is a quantitative imaging approach in soft tissues based on viscosity-elastic properties. Complex shear modulus (CSM) estimation is an effective solution to analyze tissues' physical properties for elasticity and viscosity based on the wavenumber and attenuation coefficient. CSM offers a way to detect and classify some types of soft tissues. However, CSM-based elastography inherits some obstacles, such as estimation precision and calculation complexity. This work proposes an approach for two-dimensional CSM estimation and soft tissue classification using the Extended Kalman Filter (EKF) and Decision Tree (DT) algorithm, named the EKF-DT approach. CSM estimation is obtained by applying EKF to exploit shear wave propagation at each spatial point. Afterward, the classification of tissues is done by a direct and efficient decision tree algorithm categorizing three types of normal, cirrhosis, and fibrosis liver tissues. Numerical simulation scenarios have been employed to illustrate the recovered quality and practicality of the proposed method's liver tissue classification. With the EKF, the estimated wave number and attenuation coefficient are close to the ideal values, especially the estimated wave number. The states of three liver tissue types were automatically classified by applying the DT coupled with two proposed thresholds of elasticity and viscosity: (2.310 kPa, 1.885 Pa.s) and (3.620 kPa 3.146 Pa.s), respectively. The proposed method shows the feasibility of CSM estimation based on the wavenumber and attenuation coefficient by applying the EKF. Moreover, the DT can automate the classification of liver tissue conditions by proposing two thresholds. The proposed EKF-DT method can be developed by 3D image reconstruction and empirical data before applying it in medical practice.

Imaging methods to evaluate tumor microenvironment factors affecting nanoparticle drug delivery and antitumor response

Standard small molecule and nanoparticulate chemotherapies are used for cancer treatment; however, their effectiveness remains highly variable. One reason for this variable response is hypothesized to be due to nonspecific drug distribution and heterogeneity of the tumor microenvironment, which affect tumor delivery of the agents. Nanoparticle drugs have many theoretical advantages, but due to variability in tumor microenvironment (TME) factors, the overall drug delivery to tumors and associated antitumor response are low. The nanotechnology field would greatly benefit from a thorough analysis of the TME factors that create these physiological barriers to tumor delivery and treatment in preclinical models and in patients. Thus, there is a need to develop methods that can be used to reveal the content of the TME, determine how these TME factors affect drug delivery, and modulate TME factors to increase the tumor delivery and efficacy of nanoparticles. In this review, we will discuss TME factors involved in drug delivery, and how biomedical imaging tools can be used to evaluate tumor barriers and predict drug delivery to tumors and antitumor response.

Utility of Ultrasound Strain Elastography to Differentiate Benign from Malignant Lesions of the Breast

Background

The purpose of this study was to determine the utility and diagnostic performance of strain elastography (SE) in differentiating benign from malignant lesions of the breast.

Methods

In this prospective study, 50 palpable breast masses in 50 patients were examined by mammography, B-mode ultrasound (US) and SE. Lesions were categorized using Breast Imaging Reporting and Data System (BIRADS) scoring based on mammographic and sonographic features. Elasticity scores were assessed on a five-point scale based on the distribution of strain, and the lesion size on SE imaging and B-mode (elasticity imaging/B mode [EI/B] ratio) was compared. Findings were correlated with the BIRADS assessment and diagnostic performance of sonoelastography was evaluated taking histopathology as reference standard.

Results

Histopathology revealed 29 (58%) malignant and 21 (42%) benign lesions. Infiltrative ductal carcinoma and fibroadenoma were the most common malignant and benign lesions, respectively. The sensitivity, specificity, positive predictive value, negative predictive value, and accuracy of SE was 100%, 76.1%, 85.2%, 100%, and 90%, respectively. Higher elasticity score was significantly associated with malignant histopathology (P < 0.00001). The mean EI/B ratio for malignant lesions was 1.36 ± 0.24 while that of benign lesions was 1.03 ± 0.30 (P = 0.000).

Conclusion

Real-time SE of the breast, with its superior sensitivity and specificity, could provide improved characterization of benign and malignant breast masses compared with mammography and conventional US. Due to greater diagnostic accuracy, SE can be an effective adjunctive tool to B-mode US in predicting malignancy of breast, as well as in reducing the need for biopsies in benign breast lesions.

Evaluation of internal and shell stiffness in the differential diagnosis of breast non-mass lesions by shear wave elastography

BACKGROUND

The diagnostic specificity of conventional ultrasound for breast non-mass lesions (NMLs) is low at approximately 21%-43%. Shear wave elastography (SWE) can distinguish benign from malignant lesions by evaluating the internal and peripheral stiffness. SWE has good reproducibility and high diagnostic efficacy. However, there are very few independent studies on the diagnostic value of SWE in breast NMLs.

AIM

To determine the value of SWE in the differential diagnosis of breast NMLs.

METHODS

This study enrolled a total of 118 patients with breast NMLs who underwent SWE examinations in the Beijing Shijitan Hospital Affiliated to Capital Medical University and The Second Hospital of Shandong University from January 2019 to January 2020. The internal elastic parameters of the lesions were recorded, including maximum (Emax), mean (Emean) and minimum elastic values and the standard deviation. The following peripheral parameters were noted: Presence of a “stiff rim” sign; Emax, and Emean elasticity values within 1 mm, 1.5 mm, 2 mm, 2.5 mm and 3 mm from the edge of NMLs. The receiver operating characteristic curve of each parameter was drawn, and the areas under the curve were calculated.

RESULTS

Emax, Emean and elastic values, and the standard deviation of the internal elastic values in malignant NMLs were significantly higher than those in benign NMLs (P < 0.05). The percentage with the “stiff rim” sign in malignant NMLs was significantly higher than that in the benign group (P < 0.05), and Emax and Emean at the shell of 1 mm, 1.5 mm, 2 mm, 2.5 mm and 3 mm in the malignant group were all higher than those in the benign group (P < 0.05). Of the surrounding elasticity values, Emax of the shell at 2.5 mm in malignant NMLs had maximum areas under the curve of 0.900, and the corresponding sensitivity and specificity were 94.57% and 85.86%, respectively.

CONCLUSION

The “stiff rim” sign and multiple quantitative elastic values within and around the lesion had good diagnostic performance in the differential diagnosis of breast NMLs. Emax in peripheral tissue had better diagnostic efficiency than other parameters.

Quantitative characterization of viscoelastic fracture induced by time-dependent intratumoral pressure in a 3D model tumor

In the tumor environment, interstitial pressure drives interstitial flow drainage from the tumor core to the lymphatic vessels. Recent studies have highlighted the key role of interstitial pressure in tumor development and cell migration. High intratumoral pressures, up to 60 mm Hg , have been reported in cancer patients. In a previous study, we showed that such pressure levels induce fracture in an experimental tumor model consisting of a microfluidic system holding a cellular aggregate. Here, we investigate and quantify the characteristics of tumor model fracture under a range of flow conditions. Our findings suggest a strong dependence of viscoelastic fracture behavior on the loading rate exerted by flow. The aggregate exhibits fragile fracture at high loading rates and ductile fracture at lower rates. The loading rate also modifies the permeability of the cellular aggregate, as well as the persistence time of the load required to induce fracture. The quantification parameters we propose here, evaluated for an in vitro model tumor without the extracellular matrix, could be applied to characterize tumor fracture under more realistic interstitial flow conditions.

Data-Driven Elasticity Imaging Using Cartesian Neural Network Constitutive Models and the Autoprogressive Method

Quasi-static elasticity imaging techniques rely on model-based mathematical inverse methods to estimate mechanical parameters from force-displacement measurements. These techniques introduce simplifying assumptions that preclude exploration of unknown mechanical properties with potential diagnostic value. We previously reported a data-driven approach to elasticity imaging using artificial neural networks (NNs) that circumvents limitations associated with model-based inverse methods. NN constitutive models can learn stress-strain behavior from force-displacement measurements using the autoprogressive (AutoP) method without prior assumptions of the underlying constitutive model. However, information about internal structure was required. We invented Cartesian NN constitutive models (CaNNCMs) that learn the spatial variations of material properties. We are presenting the first implementation of CaNNCMs trained with AutoP to develop data-driven models of 2-D linear-elastic materials. Both simulated and experimental force-displacement data were used as input to AutoP to show that CaNNCMs are able to model both continuous and discrete material property distributions with no prior information of internal object structure. Furthermore, we demonstrate that CaNNCMs are robust to measurement noise and can reconstruct reasonably accurate Young's modulus images from a sparse sampling of measurement data. CaNNCMs are an important step toward clinical use of data-driven elasticity imaging using AutoP.

A Feedback Loop between Hypoxia and Matrix Stress Relaxation Increases Oxygen-Axis Migration and Metastasis in Sarcoma

Upregulation of collagen matrix crosslinking directly increases its ability to relieve stress under the constant strain imposed by solid tumor, a matrix property termed stress relaxation. However, it is unknown how rapid stress relaxation in response to increased strain impacts disease progression in a hypoxic environment. Previously, it has been demonstrated that hypoxia-induced expression of the crosslinker procollagen-lysine, 2-oxoglutarate 5-dioxygenase 2 (PLOD2), in sarcomas has resulted in increased lung metastasis. Here, we show that short stress relaxation times led to increased cell migration along a hypoxic gradient in 3D collagen matrices, and rapid stress relaxation upregulated PLOD2 expression via TGFβ-SMAD2 signaling, forming a feedback loop between hypoxia and the matrix. Inhibition of this pathway led to a decrease in migration along the hypoxic gradients. In vivo, sarcoma primed in a hypoxic matrix with short stress relaxation time enhanced collagen fiber size and tumor density and increased lung metastasis. High expression of PLOD2 correlated with decreased overall survival in patients with sarcoma. Using a patient-derived sarcoma cell line, we developed a predictive platform for future personalized studies and therapeutics. Overall, these data show that the interplay between hypoxia and matrix stress relaxation amplifies PLOD2, which in turn accelerates sarcoma cell motility and metastasis. SIGNIFICANCE: These findings demonstrate that mechanical (stress relaxation) and chemical (hypoxia) properties of the tumor microenvironment jointly accelerate sarcoma motility and metastasis via increased expression of collagen matrix crosslinker PLOD2.

Reliability of Sonoelastography in Ductal Carcinoma

The objective of the study was to evaluate the reliability of sonoelastography in ductal carcinoma in patients in primary and secondary health care settings. Google Scholar, PubMed, Medline, Medscape, Wikipedia and NCBI were searched in January 2018 for all original research and review articles to identify relevant studies. Two reviewers selected the articles independently for based on the title and abstract. The selection criteria were sonoelastography accuracy for diagnosing ductal carcinoma as index text, B-mode sonography, and micropure imaging; surgery and histologic findings were used as reference text; and benign and malignant breast abnormalities and ductal carcinoma were used as target conditions. Two reviewers extracted the data on selected study characteristics, and the results were used to construct the tables and figures. Fifteen studies on ductal carcinoma were found. The overall sensitivity of sonoelastography in diagnosing ductal carcinoma was 85.7%, and the specificity, 79.8%. On the basis of the literature review, it was concluded that sonoelastography has high sensitivity and specificity in diagnosing ductal carcinoma.

Mechanisms and impact of altered tumour mechanics

The physical characteristics of tumours are intricately linked to the tumour phenotype and difficulties during treatment. Many factors contribute to the increased stiffness of tumours; from increased matrix deposition, matrix remodelling by forces from cancer cells and stromal fibroblasts, matrix crosslinking, increased cellularity, and the build-up of both solid and interstitial pressure. Increased stiffness then feeds back to increase tumour invasiveness and reduce therapy efficacy. Increased understanding of this interplay is offering new therapeutic avenues.

Complex mechanics of the heterogeneous extracellular matrix in cancer

The extracellular matrix (ECM) performs many critical functions, one of which is to provide structural and mechanical integrity, and many of the constituent proteins have clear mechanical roles. The composition and structural characteristics of the ECM are widely variable among different tissues, suiting diverse functional needs. In diseased tissues, particularly solid tumors, the ECM is complex and influences disease progression. Cancer and stromal cells can significantly influence the matrix composition and structure and thus the mechanical properties of the tumor microenvironment (TME). In this review, we describe the interactions that give rise to the structural heterogeneity of the ECM and present the techniques that are widely employed to measure ECM properties and remodeling dynamics. Furthermore, we review the tools for measuring the distinct nature of cell-ECM interactions within the TME.

Biomimetic tumor microenvironments based on collagen matrices

This review provides an overview of the current approaches to engineer defined 3D matrices for the investigation of tumor cell behaviorin vitro, with a focus on collagen-based fibrillar systems.

Elastographic Assessment of Xenograft Pancreatic Tumors

High tissue pressures prevent chemotherapeutics from reaching the parenchyma of pancreatic ductal adenocarcinoma, which makes it difficult to treat this aggressive disease. Researchers currently use invasive probes to monitor the effectiveness of pressure-reducing therapies, but this practice introduces additional complications. Here, we hypothesize that Young's modulus is a good surrogate for tissue pressure because collagen density and hyaluoronic acid, the key features of the tumor microenvironment responsible for high tissue pressures, also affect modulus elastograms. To corroborate this hypothesis, we used model-based quasi-static elastography to assess how the Young's modulus of naturally occurring AsPc-1 pancreatic tumors varies with collagen density and hyaluoronic acid concentration. We observed that Young's moduli of orthotopically grown xenograft tumors were 6 kPa (p < 0.05) higher than that of their subcutaneously grown counterparts. We also observed a strong correlation between Young's modulus and regions within the tumors with high collagen (R² ≈ 0.8) and hyaluoronic acid (R² ≈ 0.6) densities. These preliminary results indicate that hyaluronic acid and collagen density, features of the pancreatic ductal adenocarcinoma tumor microenvironment responsible for high tissue pressure, influence Young's modulus.

Viscoelasticity Mapping by Identification of Local Shear Wave Dynamics

Estimation of soft tissue elasticity is of interest in several clinical applications. For instance, tumors and fibrotic lesions are notoriously stiff compared with benign tissue. A fully quantitative measure of lesion stiffness can be obtained by shear wave (SW) elastography. This method uses an acoustic radiation force to produce laterally propagating SWs that can be tracked to obtain the velocity, which in turn is related to Young's modulus. However, not only elasticity, but also viscosity plays an important role in the propagation process of SWs. In fact, viscosity itself is a parameter of diagnostic value for the detection and characterization of malignant lesions. In this paper, we describe a new method that enables imaging viscosity from SW elastography by local model-based system identification. By testing the method on simulated data sets and performing in vitro experiments, we show that the ability of the proposed technique to generate parametric maps of the viscoelastic material properties from SW measurements, opening up new possibilities for noninvasive tissue characterization.

Model-based vascular elastography improves the detection of flow-induced carotid artery remodeling in mice

Increased arterial thickness measured with ultrasound correlates with future cardiovascular events, but conventional ultrasound imaging techniques cannot distinguish between intima, media, or atherosclerotic plaque in the carotid artery. In this work, we evaluated how well vascular elastography can detect intimal changes in a mouse model of carotid remodeling. We ligated the left external and internal branches of the carotid artery of male FVB mice and performed sham operations for 2 weeks. High-resolution ultrasound imaging accurately detected lower blood velocities and low blood volume flow in the carotid arteries after ligation in FVB mice. However, ultrasound could not detect differences in the carotid wall even at 2 weeks post-surgery. The Young’s modulus was measured based on displacements of the carotid artery wall, and Young’s modulus was 2-fold greater in shams at 1 week post ligation, and 3-fold greater 2 weeks after ligation. Finally, the higher Young’s modulus was most associated with higher intimal thickness but not medial or adventitial thickness as measured by histology. In conclusion, we developed a robust ultrasound-based elastography method for early detection of intimal changes in small animals.

Automated Compression Device for Viscoelasticity Imaging

Noninvasive measurement of tissue viscoelastic properties is gaining more attention for screening and diagnostic purposes. Recently, measuring dynamic response of tissue under a constant force has been studied for estimation of tissue viscoelastic properties in terms of retardation times. The essential part of such a test is an instrument that is capable of creating a controlled axial force and is suitable for clinical applications. Such a device should be lightweight, portable, and easy to use for patient studies to capture tissue dynamics under external stress. In this paper, we present the design of an automated compression device for studying the creep response of materials with tissue-like behaviors. The device can be used to apply a ramp-and-hold force excitation for a predetermined duration of time and it houses an ultrasound probe for monitoring the creep response of the underlying tissue. To validate the performance of the device, several creep tests were performed on tissue-mimicking phantoms, and the results were compared against those from a commercial mechanical testing instrument. Using a second-order Kelvin-Voigt model and surface measurement of the forces and displacements, retardation times T1 and T2 were estimated from each test. These tests showed strong agreement between our automated compression device and the commercial mechanical testing system, with an average relative error of 2.9% and 12.4%, for T1 and T2, respectively. Also, we present the application of compression device to measure local retardation times for four different phantoms with different size and stiffness.

Virtual Breast Quasi-static Elastography (VBQE)

Viscoelasticity Imaging (VEI) has been proposed to measure relaxation time constants for characterization of in vivo breast lesions. In this technique, an external compression force on the tissue being imaged is maintained for a fixed period of time to induce strain creep. A sequence of ultrasound echo signals is then utilized to generate time-resolved strain measurements. Relaxation time constants can be obtained by fitting local time-resolved strain measurements to a viscoelastic tissue model (e.g., a modified Kevin-Voigt model). In this study, our primary objective is to quantitatively evaluate the contrast transfer efficiency (CTE) of VEI, which contains useful information regarding image interpretations. Using an open-source simulator for virtual breast quasi-static elastography (VBQE), we conducted a case study of contrast transfer efficiency of VEI. In multiple three-dimensional (3D) numerical breast phantoms containing various degrees of heterogeneity, finite element (FE) simulations in conjunction with quasi-linear viscoelastic constitutive tissue models were performed to mimic data acquisition of VEI under freehand scanning. Our results suggested that there were losses in CTE, and the losses could be as high as -18 dB. FE results also qualitatively corroborated clinical observations, for example, artifacts around tissue interfaces.

Optimization of a Pixel-to-Pixel Curve-Fitting Method for Poroelastography Imaging

Ultrasound poroelastography is an imaging modality used to characterize the temporal behavior of soft tissue that can be modeled as a solid permeated by interconnected pores filled with liquid (poroelastic medium). It could be useful in the stage classification of lymphedema. Generally, time-constant models are applied to strain images, and precision of the fitting process, computational cost and versatility in response to changes in tissues properties are crucial aspects of clinical applications. In the work described here, we performed creep experiments on poroelastic phantoms and used rheologic models to visualize the changes in viscoelastic response associated with fluid mobility. We used the Levenberg-Marquardt algorithm as a fitting tool and performed parametric studies to improve its performance. On the basis of these studies, we proposed an optimization schema for the pixel-to-pixel curve-fitting process. We determined that the bimodal Kelvin-Voigt model describes efficiently the temporal evolution of the strain images in heterogeneous phantoms.

Design and Testing of a Single-Element Ultrasound Viscoelastography System for Point-of-Care Edema Quantification

Management of fluid overload in patients with end-stage renal disease represents a unique challenge to clinical practice because of the lack of accurate and objective measurement methods. Currently, peripheral edema is subjectively assessed by palpation of the patient's extremities, ostensibly a qualitative indication of tissue viscoelastic properties. New robust quantitative estimates of tissue fluid content would allow clinicians to better guide treatment, minimizing reactive treatment decision making. Ultrasound viscoelastography (UVE) can be used to estimate strain in viscoelastic tissue, deriving material properties that can help guide treatment. We are developing and testing a simple, low-cost UVE system using a single-element imaging transducer that is simpler and less computationally demanding than array-based systems. This benchtop validation study tested the feasibility of using the UVE system by measuring the mechanical properties of a tissue-mimicking material under large strains. We generated depth-dependent creep curves and viscoelastic parameter maps of time constants and elastic moduli for the Kelvin model of viscoelasticity. During testing, the UVE system performed well, with mean UVE-measured strain matching standard mechanical testing with maximum absolute errors ≤4%. Motion tracking revealed high correlation and signal-to-noise ratios, indicating that the system is reliable.

High-frequency ultrasound analysis of post-mitotic arrest cell death

Non-invasive monitoring of cancer cell death would permit rapid feedback on treatment response. One technique showing such promise is quantitative ultrasound. High-frequency ultrasound spectral radiofrequency analysis was used to study cell death in breast cancer cell samples. Quantitative ultrasound parameters, including attenuation, spectral slope, spectral 0-MHz-intercept, midband fit, and fitted parameters displayed significant changes with paclitaxel-induced cell death, corresponding to observations of morphological changes seen in histology and electron microscopy. In particular, a decrease in spectral slope from 0.24±0.07 dB/MHz to 0.04±0.09 dB/MHz occurred over 24 hours of treatment time and was identified as an ultrasound parameter capable of differentiating post-mitotic arrest cell death from classical apoptosis. The formation of condensed chromatin aggregates of 1 micron or greater in size increased the number of intracellular scatterers, consistent with a hypothesis that nuclear material is a primary source of ultrasound scattering in dying cells. It was demonstrated that the midband fit quantitatively correlated to cell death index, with a Pearson R-squared value of 0.99 at p<0.01. These results suggest that high-frequency ultrasound can not only qualitatively assess the degree of cancer cell death, but may be used to quantify the efficacy of chemotherapeutic treatments.

Implications of Lymphatic Transport to Lymph Nodes in Immunity and Immunotherapy

Adaptive immune response consists of many highly regulated, multistep cascades that protect against infection while preserving the health of autologous tissue. The proper initiation, maintenance, and resolution of such responses require the precise coordination of molecular and cellular signaling over multiple time and length scales orchestrated by lymphatic transport. In order to investigate these functions and manipulate them for therapy, a comprehensive understanding of how lymphatics influence immune physiology is needed. This review presents the current mechanistic understanding of the role of the lymphatic vasculature in regulating biomolecule and cellular transport from the interstitium, peripheral tissue immune surveillance, the lymph node stroma and microvasculature, and circulating lymphocyte homing to lymph nodes. This review also discusses the ramifications of lymphatic transport in immunity as well as tolerance and concludes with examples of how lymphatic-mediated targeting of lymph nodes has been exploited for immunotherapy applications.

Acoustic Radiation Force Impulse Technology in the Differential Diagnosis of Solid Breast Masses with Different Sizes: Which Features Are Most Efficient?

Purpose. To evaluate diagnostic performance of acoustic radiation force impulse (ARFI) technology for solid breast masses with different sizes and determine which features are most efficient. Materials and Methods. 271 solid breast masses in 242 women were examined with ARFI, and their shear wave velocities (SWVs), Virtual Touch tissue imaging (VTI) patterns, and area ratios (ARs) were measured and compared with their histopathological outcomes. Receiver operating characteristic curves (ROC) were calculated to assess diagnostic performance of ARFI for small masses (6–14 mm) and big masses (15–30 mm). Results. SWV of mass was shown to be positively associated with mass size (P < 0.001). For small masses, area under ROC (Az) of AR was larger than that of SWV (P < 0.001) and VTI pattern (P < 0.001); no significant difference was found between Az of SWV and that of VTI pattern (P = 0.906). For big masses, Az of VTI pattern was less than that of SWV (P = 0.008) and AR (P = 0.002); no significant difference was identified between Az of SWV and that of AR (P = 0.584). Conclusions. For big masses, SWV and AR are both efficient measures; nevertheless, for small masses, AR seems to be the best feature.

Lymph node biophysical remodeling is associated with melanoma lymphatic drainage

Tissue remodeling is a characteristic of many solid tumor malignancies including melanoma. By virtue of tumor lymphatic transport, remodeling pathways active within the local tumor microenvironment have the potential to be operational within lymph nodes (LNs) draining the tumor interstitium. Here, we show that lymphatic drainage from murine B16 melanomas in syngeneic, immune-competent C57Bl/6 mice is associated with LN enlargement as well as nonuniform increases in bulk tissue elasticity and viscoelasticity, as measured by the response of whole LNs to compression. These remodeling responses, which quickly manifest in tumor-draining lymph nodes (TDLNs) after tumor inoculation and before apparent metastasis, were accompanied by changes in matrix composition, including up to 3-fold increases in the abundance of soluble collagen and hyaluronic acid. Intranodal pressures were also significantly increased in TDLNs (+1 cmH2O) relative to both non-tumor-draining LNs (-1 cmH2O) and LNs from naive animals (-1 to 2 cmH2O). These data suggest that the reorganization of matrix structure, composition, and fluid microenvironment within LNs associated with tumor lymphatic drainage parallels remodeling seen in primary malignancies and has the potential to regulate the adhesion, proliferation, and signaling function of LN-resident cells involved in directing melanoma disease progression.

Bio-optics based sensation imaging for breast tumor detection using tissue characterization

The tissue inclusion parameter estimation method is proposed to measure the stiffness as well as geometric parameters. The estimation is performed based on the tactile data obtained at the surface of the tissue using an optical tactile sensation imaging system (TSIS). A forward algorithm is designed to comprehensively predict the tactile data based on the mechanical properties of tissue inclusion using finite element modeling (FEM). This forward information is used to develop an inversion algorithm that will be used to extract the size, depth, and Young's modulus of a tissue inclusion from the tactile data. We utilize the artificial neural network (ANN) for the inversion algorithm. The proposed estimation method was validated by a realistic tissue phantom with stiff inclusions. The experimental results showed that the proposed estimation method can measure the size, depth, and Young's modulus of a tissue inclusion with 0.58%, 3.82%, and 2.51% relative errors, respectively. The obtained results prove that the proposed method has potential to become a useful screening and diagnostic method for breast cancer.

Parametric imaging of viscoelasticity using optical coherence elastography

We demonstrate imaging of soft tissue viscoelasticity using optical coherence elastography. Viscoelastic creep deformation is induced in tissue using step-like compressive loading and the resulting time-varying deformation is measured using phase-sensitive optical coherence tomography. From a series of co-located B-scans, we estimate the local strain rate as a function of time, and parameterize it using a four-parameter Kelvin-Voigt model of viscoelastic creep. The estimated viscoelastic strain and time constant are used to visualize viscoelastic creep in 2D, dual-parameter viscoelastograms. We demonstrate our technique on six silicone tissue-simulating phantoms spanning a range of viscoelastic parameters. As an example in soft tissue, we report viscoelastic contrast between muscle and connective tissue in fresh, ex vivo rat gastrocnemius muscle and mouse abdominal transection. Imaging viscoelastic creep deformation has the potential to provide complementary contrast to existing imaging modalities, and may provide greater insight into disease pathology.

The Added Value of Color Parameters in Analyzing Elastographic Images of Ultrasound Detected Breast Focal Lesions

AIMS

The purpose of the study was to determine if the color quantitative analysis obtained on elastographic images of breast lesions could improve the benign-malignant differentiation, and also to identify some of the circumstances which would benefit most from such an analysis.

PATIENTS AND METHODS

The study design was a longitudinal prospective one, all data being acquired between May 2007 and September 2008. The US device used: Hitachi 8500 EUB machine with elastography option. For suspicious breast lesions histopathology was obtained by means of percutaneous biopsy or post-surgery. Studied color parameters (numeric values): average color (red, green, blue), color dispersion, average intensity, average hue, hue dispersion. Calculus modality: Image Processing Version 1.3, a program developed in collaboration with the Technical University of Cluj Napoca.

RESULTS

Seventy-one (71) women were selected for the study. A hundred and six circumscribed breast lesions were detected by means of ultrasound in the studied group. Five color parameters were independently associated with the histological diagnosis (AvgBlue, AvgGreen and AvgRed; DispRed and DispIntensity) with AvgBlue parameter making the most important contribution (p<0.0001); the greater the values of AvgBlue (more than 92), the higher the chances of malignancy and the greater the values of AvgGreen (more than 88), the higher the chances for a benign lesion.

CONCLUSION

High numeric values for Avg Blue (more than 92) would increase the probability of malignancy and thus recommend a more aggressive diagnostic management (biopsy), while high numeric values for AvgGreen (more than 88) would reassure the examiner to proceed conservatively with short interval or routine follow-ups.

The tactile sensation imaging system for embedded lesion characterization

Elasticity is an important indicator of tissue health, with increased stiffness pointing to an increased risk of cancer. We investigated a tissue inclusion characterization method for the application of early breast tumor identification. A tactile sensation imaging system (TSIS) is developed to capture images of the embedded lesions using total internal reflection principle. From tactile images, we developed a novel method to estimate that size, depth, and elasticity of the embedded lesion using 3-D finite-element-model-based forward algorithm, and neural-network-based inversion algorithm are employed. The proposed characterization method was validated by the realistic tissue phantom with inclusions to emulate the tumors. The experimental results showed that, the proposed characterization method estimated the size, depth, and Young's modulus of a tissue inclusion with 6.98%, 7.17%, and 5.07% relative errors, respectively. A pilot clinical study was also performed to characterize the lesion of human breast cancer patients using TSIS.

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.

Ultrasound elastography: advantages, limitations and artefacts of the different techniques from a study on a phantom

Ultrasound elastography is a technique currently under development. Its use in clinical practice is complicated because of the wide range of techniques used by the different manufacturers and the parameters proposed to characterise tissues. A comparative analysis on five ultrasound diagnostic systems has been performed on a calibrated elasticity phantom and demonstrated that: (1) all systems tested are reliable for simple qualitative analysis: is a nodule present and is it harder or softer than neighbouring tissues? (2) the deformation or hardness ratios between two regions are usually, however, not proportional to the theoretical ratios and only a binary analysis greater than 1 (harder) and less than 1 (softer) is reliable and could be used as a negative predictive value (NPV) for malignant lesions, as has been suggested by some authors; (3) finally, quantitative analysis using shear wave techniques performed variably, reliable measurements being obtained with only one of the systems. Measurements produced by these different systems must not be compared in clinical practice to monitor a patient and the threshold values proposed in the literature must only be used in an analysis carried out with the same system and same probe.

Real-time 1-D/2-D transient elastography on a standard ultrasound scanner using mechanically induced vibration

Transient elastography has been well established in the literature as a means of assessing the elasticity of soft tissue. In this technique, tissue elasticity is estimated from the study of the propagation of the transient shear waves induced by an external or internal source of vibration. Previous studies have focused mainly on custom single-element transducers and ultrafast scanners which are not available in a typical clinical setup. In this work, we report the design and implementation of a transient elastography system on a standard ultrasound scanner that enables quantitative assessment of tissue elasticity in real-time. Two new custom imaging modes are introduced that enable the system to image the axial component of the transient shear wave, in response to an externally induced vibration, in both 1-D and 2-D. Elasticity reconstruction algorithms that estimate the tissue elasticity from these transient waves are also presented. Simulation results are provided to show the advantages and limitations of the proposed system. The performance of the system is also validated experimentally using a commercial elasticity phantom.

Breast elastography: A literature review

Breast elastography is a new sonographic imaging technique which provides information on breast lesions in addition to conventional ultrasonography (US) and mammography. Elastography provides a noninvasive evaluation of the stiffness of a lesion. Today, two technical solutions are available for clinical use: strain elastography and shear wave elastography. Initial evaluations of these techniques in clinical trials suggest that they may substantially improve the possibility of differentiating benign from malignant breast lesions thereby limiting recourse to biopsy and considerably reducing the number of benign breast biopsy diagnoses. This article reviews the basics of this technique, how to perform the examination, image interpretation and the results of major clinical studies. Although elastography is easy to perform, training and technical knowledge are required in order to obtain images permitting a correct interpretation. This paper will highlight the technique and point out common pitfalls.

Breast elastography is a new sonographic imaging technique which provides information on breast lesions in addition to conventional ultrasonography (US) and mammography. Elastography provides a noninvasive evaluation of the stiffness of a lesion. Today, two technical solutions are available for clinical use: strain elastography and shear wave elastography. Initial evaluations of these techniques in clinical trials suggest that they may substantially improve the possibility of differentiating benign from malignant breast lesions thereby limiting recourse to biopsy and considerably reducing the number of benign breast biopsy diagnoses. This article reviews the basics of this technique, how to perform the examination, image interpretation and the results of major clinical studies. Although elastography is easy to perform, training and technical knowledge are required in order to obtain images permitting a correct interpretation. This paper will highlight the technique and point out common pitfalls.

Edge-preserving speckle texture removal by interference-based speckle filtering followed by anisotropic diffusion

Ultrasonography has an inherent noise pattern, called speckle, which is known to hamper object recognition for both humans and computers. Speckle noise is produced by the mutual interference of a set of scattered wavefronts. Depending on the phase of the wavefronts, the interference may be constructive or destructive, which results in brighter or darker pixels, respectively. We propose a filter that minimizes noise fluctuation while simultaneously preserving local gray level information. It is based on steps to attenuate the destructive and constructive interference present in ultrasound images. This filter, called interference-based speckle filter followed by anisotropic diffusion (ISFAD), was developed to remove speckle texture from B-mode ultrasound images, while preserving the edges and the gray level of the region. The ISFAD performance was compared with 10 other filters. The evaluation was based on their application to images simulated by Field II (developed by Jensen et al.) and the proposed filter presented the greatest structural similarity, 0.95. Functional improvement of the segmentation task was also measured, comparing rates of true positive, false positive and accuracy. Using three different segmentation techniques, ISFAD also presented the best accuracy rate (greater than 90% for structures with well-defined borders).

Strain-compounding technique with ultrasound Nakagami imaging for distinguishing between benign and malignant breast tumors

PURPOSE

The scatterer properties of breast tissues are related to the presence of collagen structures, while the elasticity properties of breast tissues depend on their structural organization; these two characteristics are functionally complementary in ultrasound-based tissue characterizations. This study investigated the use of a strain-compounding technique with Nakagami imaging to provide information associated with the scatterer and elasticity characteristics of tissues when attempting to identify benign and malignant breast tumors.

METHODS

The efficacy of the proposed method was tested by collecting raw data of ultrasound backscattered signals from 50 clinical cases (25 benign tumors and 25 malignant tumors, as verified by histology biopsies). The different strain conditions were created by applying manual compression. For each region in which breast tumors were suspected, estimates of the full width at half maximum (FWHM) from the Gaussian fitting curve for the Nakagami-parameter histogram in the strain-compounding Nakagami images were divided by those of the corresponding reference Nakagami images (uncompressed images); this parameter was denoted as the FWHM ratio. Receiver operating characteristic (ROC) curve analysis was adopted to assess the diagnostic performance.

RESULTS

The results demonstrated that the difference in scatterer distributions between before and after compounding was greater for benign tumors than for malignant tumors. The FWHM ratio estimates for benign and malignant tumors were 0.76 ± 0.14 and 0.96 ± 0.06 (mean ± standard deviation), respectively (p < 0.01). The mean area under the ROC curve using the FWHM ratio estimates was 0.92, with a 95% confidence interval of 0.83-1.00.

CONCLUSIONS

These findings indicate that the strain-compounding Nakagami imaging method based on the acquisition of multiple frames under different strain states could provide objective information that would improve the ability to classify benign and malignant breast tumors.

Breast imaging: how we manage diagnostic technology at a multidisciplinary breast center

This paper discusses the most important aspects and problems related to the management of breast cancer imaging, at a center specialized in breast pathology. We review the established and emerging diagnostic techniques, their indications, and peculiarities: digital mammography, CAD systems, and the recent digital breast tomosynthesis, ultrasound and complementary elastography, molecular imaging techniques, magnetic resonance imaging, advanced sequences (diffusion), and positron emission mammography (PEM). The adequate integration and rational management of these techniques is essential, but this is not always easy, in order to achieve a successful diagnosis.

Review of tissue simulating phantoms with controllable optical, mechanical and structural properties for use in optical coherence tomography

We review the development of phantoms for optical coherence tomography (OCT) designed to replicate the optical, mechanical and structural properties of a range of tissues. Such phantoms are a key requirement for the continued development of OCT techniques and applications. We focus on phantoms based on silicone, fibrin and poly(vinyl alcohol) cryogels (PVA-C), as we believe these materials hold the most promise for durable and accurate replication of tissue properties.

Imaging feedback of histotripsy treatments using ultrasound shear wave elastography

Histotripsy is a cavitation-based ultrasound therapy that mechanically fractionates soft solid tissues into fluid-like homogenates. This paper investigates the feasibility of imaging the tissue elasticity change during the histotripsy process as a tool to provide feedback for the treatments. The treatments were performed on agar tissue phantoms and ex vivo kidneys using 3-cycle ultrasound pulses delivered by a 750-kHz therapeutic array at peak negative/positive pressure of 17/108 MPa and a repetition rate of 50 Hz. Lesions with different degrees of damage were created with increasing numbers of therapy pulses from 0 to 2000 pulses per treatment location. The elasticity of the lesions was measured with ultrasound shear wave elastography, in which a quasi-planar shear wave was induced by acoustic radiation force generated by the therapeutic array, and tracked with ultrasound imaging at 3000 frames per second. Based on the shear wave velocity calculated from the sequentially captured frames, the Young's modulus was reconstructed. Results showed that the lesions were more easily identified on the shear wave velocity images than on B-mode images. As the number of therapy pulses increased from 0 to 2000 pulses/location, the Young's modulus decreased exponentially from 22.1 ± 2.7 to 2.1 ± 1.1 kPa in the tissue phantoms (R2 = 0.99, N = 9 each), and from 33.0 ± 7.1 to 4.0 ± 2.5 kPa in the ex vivo kidneys (R2 = 0.99, N = 8 each). Correspondingly, the tissues transformed from completely intact to completely fractionated as examined via histology. A good correlation existed between the lesions' Young's modulus and the degree of tissue fractionation as examined with the percentage of remaining structurally intact cell nuclei (R2 = 0.91, N = 8 each). These results indicate that lesions produced by histotripsy can be detected with high sensitivity using shear wave elastography. Because the decrease in the tissue elasticity corresponded well with the morphological and histological change, this study provides a basis for predicting the local treatment outcomes from tissue elasticity change.

Harmonic Motion Imaging (HMI) for Tumor Imaging and Treatment Monitoring

Palpation is an established screening procedure for the detection of several superficial cancers including breast, thyroid, prostate, and liver tumors through both self and clinical examinations. This is because solid masses typically have distinct stiffnesses compared to the surrounding normal tissue. In this paper, the application of Harmonic Motion Imaging (HMI) for tumor detection based on its stiffness as well as its relevance in thermal treatment is reviewed. HMI uses a focused ultrasound (FUS) beam to generate an oscillatory acoustic radiation force for an internal, non-contact palpation to internally estimate relative tissue hardness. HMI studies have dealt with the measurement of the tissue dynamic motion in response to an oscillatory acoustic force at the same frequency, and have been shown feasible in simulations, phantoms, ex vivo human and bovine tissues as well as animals in vivo. Using an FUS beam, HMI can also be used in an ideal integration setting with thermal ablation using high-intensity focused ultrasound (HIFU), which also leads to an alteration in the tumor stiffness. In this paper, a short review of HMI is provided that encompasses the findings in all the aforementioned areas. The findings presented herein demonstrate that the HMI displacement can accurately depict the underlying tissue stiffness, and the HMI image of the relative stiffness could accurately detect and characterize the tumor or thermal lesion based on its distinct properties. HMI may thus constitute a non-ionizing, cost-efficient and reliable complementary method for noninvasive tumor detection, localization, diagnosis and treatment monitoring.

Detection of Breast Lesions using an Automated Breast Volume Scanner System

OBJECTIVE: This study investigated the clinical utility of an automated breast volume scanner (ABVS) system for the detection of breast lesions. METHODS: The breasts of 81 patients referred for ultrasonographic examination were scanned using the ABVS system and handheld ultrasonography independently by two experienced examiners. The ABVS was used to perform scans of the breast in three directions (anteroposterior, lateral and medial), with the addition of further inferior and superior scans if necessary. The scanning data were then stored and automatically reconstructed. For hand-held ultrasonography the whole breast was scanned radially from the outside to the centre of the nipple. RESULTS: The numbers of lesions reported by the two examiners were 89 and 99, respectively, using the ABVS (not statistically significant), compared with 60 and 85, respectively, using handheld ultrasonography (statistically significant). CONCLUSIONS: The ABVS system is an operator-independent method for automated breast scanning. It detected more breast lesions and provided additional information for the diagnosis of intraductal and malignant lesions compared with hand-held ultrasonography.

Principles of ultrasound elastography

Many manufacturers of ultrasound systems offer elastography capabilities. Though the research literature describes techniques involving external actuators, commercial systems have preferred to adopt techniques that may be used with a single transducer. Techniques may be divided into those which measure strain and those which measure shear wave velocity. Strain elastography involves deformation of the tissue followed by imaging of the degree of compression or extension of the tissue. Strain elastography does not estimate tissue stiffness; however, the strain ratio may be used as a surrogate index of stiffness. Shear wave elastography provides true quantitative information on elastic modulus. This involves induction of shear waves, estimation of shear wave velocity c s and conversion to elastic modulus E using the equation E = 3 ρ cs2 where ρ is the density. The description of tissues as being purely elastic is simplistic. In practice they may also exhibit time-dependent viscous behaviour. Recent literature describe methods that have been developed for estimation of the viscoelastic behaviour from the change in strain with time or by estimation of the shear wave dispersion, a technique known as ‘shear wave spectroscopy’. These methods may become commercially available in the medium term, offering a new quantity (tissue viscosity) for diagnostic use.

Real-time quasi-static ultrasound elastography

Ultrasound elastography is a technique used for clinical imaging of tissue stiffness with a conventional ultrasound machine. It was first proposed two decades ago, but active research continues in this area to the present day. Numerous clinical applications have been investigated, mostly related to cancer imaging, and though these have yet to prove conclusive, the technique has seen increasing commercial and clinical interest. This paper presents a review of the most widely adopted, non-quantitative, techniques focusing on technical innovations rather than clinical applications. The review is not intended to be exhaustive, concentrating instead on placing the various techniques in context according to the authors' perspective of the field.