Supersonic Shear Wave Elastography of Response to Anti-cancer Therapy in a Xenograft Tumor Model

Evaluating ECM stiffness and liver cancer radiation response via shear-wave elasticity in 3D culture models

BACKGROUND

The stiffness of the tumor microenvironment (TME) directly influences cellular behaviors. Radiotherapy (RT) is a common treatment for solid tumors, but the TME can impact its efficacy. In the case of liver cancer, clinical observations have shown that tumors within a cirrhotic, stiffer background respond less to RT, suggesting that the extracellular matrix (ECM) stiffness plays a critical role in the development of radioresistance.

METHODS

This study explored the effects of ECM stiffness and the inhibition of lysyl oxidase (LOX) isoenzymes on the radiation response of liver cancer in a millimeter-sized three-dimensional (3D) culture. We constructed a cube-shaped ECM-based millimeter-sized hydrogel containing Huh7 human liver cancer cells. By modulating the collagen concentration, we produced two groups of samples with different ECM stiffnesses to mimic the clinical scenarios of normal and cirrhotic livers. We used a single-transducer system for shear-wave-based elasticity measurement, to derive Young's modulus of the 3D cell culture to investigate how the ECM stiffness affects radiosensitivity. This is the first demonstration of a workflow for assessing radiation-induced response in a millimeter-sized 3D culture.

RESULTS

Increased ECM stiffness was associated with a decreased radiation response. Moreover, sonoporation-assisted LOX inhibition with BAPN (β-aminopropionitrile monofumarate) significantly decreased the initial ECM stiffness and increased RT-induced cell death. Inhibition of LOX was particularly effective in reducing ECM stiffness in stiffer matrices. Combining LOX inhibition with RT markedly increased radiation-induced DNA damage in cirrhotic liver cancer cells, enhancing their response to radiation. Furthermore, LOX inhibition can be combined with sonoporation to overcome stiffness-related radioresistance, potentially leading to better treatment outcomes for patients with liver cancer.

CONCLUSIONS

The findings underscore the significant influence of ECM stiffness on liver cancer's response to radiation. Sonoporation-aided LOX inhibition emerges as a promising strategy to mitigate stiffness-related resistance, offering potential improvements in liver cancer treatment outcomes.

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.

Preclinical Three-Dimensional Vibrational Shear Wave Elastography for Mapping of Tumour Biomechanical Properties In Vivo

Preclinical investigation of the biomechanical properties of tissues and their treatment-induced changes are essential to support drug-discovery, clinical translation of biomarkers of treatment response, and studies of mechanobiology. Here we describe the first use of preclinical 3D elastography to map the shear wave speed (cs), which is related to tissue stiffness, in vivo and demonstrate the ability of our novel 3D vibrational shear wave elastography (3D-VSWE) system to detect tumour response to a therapeutic challenge. We investigate the use of one or two vibrational sources at vibrational frequencies of 700, 1000 and 1200 Hz. The within-subject coefficients of variation of our system were found to be excellent for 700 and 1000 Hz and 5.4 and 6.2%, respectively. The relative change in cs measured with our 3D-VSWE upon treatment with an anti-vascular therapy ZD6126 in two tumour xenografts reflected changes in tumour necrosis. U-87 MG drug vs vehicle: Δcs = −24.7 ± 2.5 % vs 7.5 ± 7.1%, (p = 0.002) and MDA-MB-231 drug vs vehicle: Δcs = −12.3 ± 2.7 % vs 4.5 ± 4.7%, (p = 0.02). Our system enables rapid (<5 min were required for a scan length of 15 mm and three vibrational frequencies) 3D mapping of quantitative tumour viscoelastic properties in vivo, allowing exploration of regional heterogeneity within tumours and speedy recovery of animals from anaesthesia so that longitudinal studies (e.g., during tumour growth or following treatment) may be conducted frequently.

Multiparametric Ultrasound of Cervical Lymph Node Metastases in Head and Neck Cancer for Planning Non-Surgical Therapy

Background: We aimed to evaluate multiparametric ultrasound, to achieve a better understanding of the baseline characteristics of suspected cervical lymph node metastases in head and neck cancer before induction chemotherapy or chemoradiation. Methods: From February 2020 to April 2021, our complete ultrasound examination protocol was carried out on clinically evident malignant lymph nodes of histologically proven HNSCC in the pre-therapeutic setting. Results: A total of 13 patients were eligible for analysis. Using elastography, irregular clear hardening in areas in the center of the lymph node could be detected in all cases. Elastographic Q-analysis showed a significantly softer cortex compared to the center and surrounding tissue. The time–intensity curve analysis showed high values for the area under the curve and a short time-to-peak (fast wash-in) in all cases compared to the surrounding tissue. A parametric evaluation of contrast enhanced the ultrasound in the early arterial phase and showed an irregular enhancement from the margin in almost all investigated lymph nodes. These results show that the implementation of comprehensive, multiparametric ultrasound is suitable for classifying suspected lymph node metastasis more precisely than conventional ultrasound alone in the pre-therapeutic setting of HNSCC. Thus, these parameters may be used for improvements in the re-staging after chemoradiation or neoadjuvant therapy monitoring, respectively.

Non-Invasive Ultrasonic Description of Tumor Evolution

Simple Summary

During tumor evolution, heterogeneous structural and functional changes occur in the tumor microenvironment. These complex changes have pro- or anti-tumorigenesis effects and have an impact on therapy efficiency. Therefore, the tumor microenvironment needs to be non-invasively characterized over time. The aim of this preclinical work is to compare the sensitivity of modifications occurring during tumor evolution of volume, immunohistochemistry and non-invasive quantitative ultrasound parameters (Shear Wave Elastography and dynamic Contrast-Enhanced Ultrasound) and to study the link between them. The complementary evaluation over time of multiple morphological and functional parameters during tumor growth underlines the need to integrate histological, morphological, functional, and, ultimately, genomic information into models that can consider the temporal and spatial variability of features to better understand tumor evolution.

Abstract

Purpose: There is a clinical need to better non-invasively characterize the tumor microenvironment in order to reveal evidence of early tumor response to therapy and to better understand therapeutic response. The goals of this work are first to compare the sensitivity to modifications occurring during tumor growth for measurements of tumor volume, immunohistochemistry parameters, and emerging ultrasound parameters (Shear Wave Elastography (SWE) and dynamic Contrast-Enhanced Ultrasound (CEUS)), and secondly, to study the link between the different parameters. Methods: Five different groups of 9 to 10 BALB/c female mice with subcutaneous CT26 tumors were imaged using B-mode morphological imaging, SWE, and CEUS at different dates. Whole-slice immunohistological data stained for the nuclei, T lymphocytes, apoptosis, and vascular endothelium from these tumors were analyzed. Results: Tumor volume and three CEUS parameters (Time to Peak, Wash-In Rate, and Wash-Out Rate) significantly changed over time. The immunohistological parameters, CEUS parameters, and SWE parameters showed intracorrelation. Four immunohistological parameters (the number of T lymphocytes per mm² and its standard deviation, the percentage area of apoptosis, and the colocalization of apoptosis and vascular endothelium) were correlated with the CEUS parameters (Time to Peak, Wash-In Rate, Wash-Out Rate, and Mean Transit Time). The SWE parameters were not correlated with the CEUS parameters nor with the immunohistological parameters. Conclusions: US imaging can provide additional information on tumoral changes. This could help to better explore the effect of therapies on tumor evolution, by studying the evolution of the parameters over time and by studying their correlations.

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

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

Early assessment of shear wave elastography parameters foresees the response to neoadjuvant chemotherapy in patients with invasive breast cancer

BACKGROUND

Early prediction of tumor response to neoadjuvant chemotherapy (NACT) is crucial for optimal treatment and improved outcome in breast cancer patients. The purpose of this study is to investigate the role of shear wave elastography (SWE) for early assessment of response to NACT in patients with invasive breast cancer.

METHODS

In a prospective study, 62 patients with biopsy-proven invasive breast cancer were enrolled. Three SWE studies were conducted on each patient: before, at mid-course, and after NACT but before surgery. A new parameter, mass characteristic frequency (fmass), along with SWE measurements and mass size was obtained from each SWE study visit. The clinical biomarkers were acquired from the pre-NACT core-needle biopsy. The efficacy of different models, generated with the leave-one-out cross-validation, in predicting response to NACT was shown by the area under the receiver operating characteristic curve and the corresponding sensitivity and specificity.

RESULTS

A significant difference was found for SWE parameters measured before, at mid-course, and after NACT between the responders and non-responders. The combination of Emean2 and mass size (s2) gave an AUC of 0.75 (0.95 CI 0.62-0.88). For the ER+ tumors, the combination of Emean_ratio1, s1, and Ki-67 index gave an improved AUC of 0.84 (0.95 CI 0.65-0.96). For responders, fmass was significantly higher during the third visit.

CONCLUSIONS

Our study findings highlight the value of SWE estimation in the mid-course of NACT for the early prediction of treatment response. For ER+ tumors, the addition of Ki-67improves the predictive power of SWE. Moreover, fmass is presented as a new marker in predicting the endpoint of NACT in responders.

Stiffness increases with myofibroblast content and collagen density in mesenchymal high grade serous ovarian cancer

Women diagnosed with high-grade serous ovarian cancers (HGSOC) are still likely to exhibit a bad prognosis, particularly when suffering from HGSOC of the Mesenchymal molecular subtype (50% cases). These tumors show a desmoplastic reaction with accumulation of extracellular matrix proteins and high content of cancer-associated fibroblasts. Using patient-derived xenograft mouse models of Mesenchymal and Non-Mesenchymal HGSOC, we show here that HGSOC exhibit distinct stiffness depending on their molecular subtype. Indeed, tumor stiffness strongly correlates with tumor growth in Mesenchymal HGSOC, while Non-Mesenchymal tumors remain soft. Moreover, we observe that tumor stiffening is associated with high stromal content, collagen network remodeling, and MAPK/MEK pathway activation. Furthermore, tumor stiffness accompanies a glycolytic metabolic switch in the epithelial compartment, as expected based on Warburg's effect, but also in stromal cells. This effect is restricted to the central part of stiff Mesenchymal tumors. Indeed, stiff Mesenchymal tumors remain softer at the periphery than at the core, with stromal cells secreting high levels of collagens and showing an OXPHOS metabolism. Thus, our study suggests that tumor stiffness could be at the crossroad of three major processes, i.e. matrix remodeling, MEK activation and stromal metabolic switch that might explain at least in part Mesenchymal HGSOC aggressiveness.

Evaluation of therapeutic effect of targeting nanobubbles conjugated with NET-1 siRNA by shear wave elastography: an in vivo study of hepatocellular carcinoma bearing mice model

This study aimed at investigating the tumor stiffness of hepatocellular carcinoma (HCC) bearing mice model in vivo to evaluate the therapeutic efficacy of targeting nanobubbles (TNBS) conjugated with NET-1 siRNA (NET-1 siRNA-TNBS). Also tested whether shear wave elastography (SWE) could demonstrate the pathological tumor changes and used to monitor therapeutic efficacy as a noninvasive method. The HCC bearing mice model was established by injecting human HCC cell line (HepG2). The mice were then divided into three groups randomly, and were treated with TNBS conjugated with NET-1 siRNA, TNBS conjugated with negative control gene, and saline as control. US-SWE was performed for three times. SWE values of all the tumors in three groups were increased with tumor growth. Emax was correlated with tumor size (p < .05). NET-1 gene (treatment group) significantly delayed the growth of tumor size compared to other two groups (p < .0001), showing a significantly increased Emax (p < .05). Immunohistochemical results showed that the NET-1 protein expression was significantly lower than the negative control and blank groups. In conclusion, TNBS conjugated with NET-1 siRNA inhibited tumor growth and prolonged the life of experimental animals. SWE provided a noninvasive and real time imaging method to detect the changes in tumor development.

Estimation of Vascular Permeability in Irregularly Shaped Cancers Using Ultrasound Poroelastography

OBJECTIVE

Vascular permeability (VP) is a mechanical parameter which plays an important role in cancer initiation, metastasis, and progression. To date, there are only a few non-invasive methods that can be used to image VP in solid tumors. Most of these methods require the use of contrast agents and are expensive, limiting widespread use.

METHODS

In this paper, we propose a new method to image VP in tumors, which is based on the use of ultrasound poroelastography. Estimation of VP by poroelastography requires knowledge of the Young's modulus (YM), Poisson's ratio (PR), and strain time constant (TC) in the tumors. In our method, we find the ellipse which best fits the tumor (regardless of its shape) using an eigen-system-based fitting technique and estimate the YM and PR using Eshelby's elliptic inclusion formulation. A Fourier method is used to estimate the axial strain TC, which does not require any initial guess and is highly robust to noise.

RESULTS

It is demonstrated that the proposed method can estimate VP in irregularly shaped tumors with an accuracy of above [Formula: see text] using ultrasound simulation data with signal-to-noise ratio of 20 dB or higher. In vivo feasibility of the proposed technique is demonstrated in an orthotopic mouse model of breast cancer.

CONCLUSION

The proposed imaging method can provide accurate and localized estimation of VP in cancers non-invasively and cost-effectively.

SIGNIFICANCE

Accurate and non-invasive assessment of VP can have a significant impact on diagnosis, prognosis, and treatment of cancers.

Design and Demonstration of a Configurable Imaging Platform for Combined Laser, Ultrasound, and Elasticity Imaging

This paper introduces a configurable combined laser, ultrasound, and elasticity (CLUE) imaging platform. The CLUE platform enables imaging sequences capable of simultaneously providing quantitative acoustic, optical, and mechanical contrast for comprehensive diagnosis and monitoring of complex diseases, such as cancer. The CLUE imaging platform was developed on a Verasonics ultrasound scanner integrated with a pulsed laser, and it was designed to be modular and scalable to allow researchers to create their own specific imaging sequences efficiently. The CLUE imaging platform and sequence were demonstrated in a tissue-mimicking phantom containing a stiff inclusion labeled with optically-activated nanodroplets and in an ex vivo mouse spleen. We have shown that CLUE imaging can simultaneously capture multi-functional imaging signals providing quantitative information on tissue.

VTIQ evaluates antitumor effects of NET-1 siRNA by UTMD in HCC xenograft models

The present study used a virtual touch tissue imaging and quantification (VTIQ) method to investigate the change in elasticity in xenograft tumor tissue models following silencing of the neuroepithelial-transforming protein 1 (NET-1) gene by ultrasound-targeted microbubble destruction (UTMD). A total of 24 xenograft models were established by subcutaneous injection of human hepatocellular carcinoma SMMC-7721 cells in BALB/c female nude mice. Then, NET-1 small interfering RNA (siRNA)-conjugated nanobubbles and a glypican-3 antibody were synthesized. The mean and maximum shear wave speed (SWS_mean and SWS_max) in the tumor tissue were measured prior to, during, and following therapy using VTIQ. The growth of the tumor size and survival time were recorded. The levels of NET-1 protein were evaluated by immunohistochemical staining. In addition, tumor, liver and kidney tissues of the nude mice were collected to confirm whether gene transfection treatment was toxic in vivo. In the UTMD delivery gene group, SWS_mean was correlated with the maximum diameter of the tumor (r=0.9806, P=0.0194). The immunohistochemical staining data indicated that the level of NET-1 protein in the treated groups was significantly decreased compared with those in the control groups. Additionally, no structural damage was observed in the nude mice liver and kidney tissues following treatment. Therefore, VTIQ measurement identified potential changes in the elastic properties of the tumors, which in turn may be associated with the stages of tumor development. The delivery method, UTMD, improves the antitumor effects of NET-1 siRNA and supports gene transfection as a promising therapeutic strategy.

Ultrasound elastography as an imaging biomarker for detection of early tumor response to chemotherapy in a murine breast cancer model: a feasibility study

OBJECTIVE

This study investigated the feasibility of using strain elastography (SE) and real time shear wave elastography (RT-SWE) to evaluate early tumor response to cytotoxic chemotherapy in a murine xenograft breast cancer tumor model.

METHODS

MCF-7 breast cancer-bearing nude mice were treated with either cisplatin 2 mg kg^-1 plus paclitaxel 10 mg kg^-1 (treatment group) or sterile saline (control group) once daily for 5 days. The tumor elasticity was measured by SE or RT-SWE before and after therapy. Tumor cell density was assessed by hematoxylin and eosin staining, and the ratio of collagen fibers in the tumor was evaluated by Van Gieson staining. The correlation between tumor elasticity, as determined by SE and SWE, as well as the pathological tumor responses were analyzed.

RESULTS

Chemotherapy significantly attenuated tumor growth compared to the control treatment (p < 0.05). Chemotherapy also significantly increased tumor stiffness (p < 0.05) and significantly decreased (p < 0.05) tumor cell density compared with the control. Moreover, chemotherapy significantly increased the ratio of collagen fibers (p < 0.05). Tumor stiffness was positively correlated with the ratio of collagen fibers but negatively correlated with tumor cell density.

CONCLUSION

The study suggests that ultrasound elastography by SE and SWE is a feasible tool for assessing early responses of breast cancer to chemotherapy in our murine xenograft model. Advances in knowledge: This study showed that the tumor elasticity determined by ultrasound elastography could be a feasible imaging biomarker for assessing very early therapeutic responses to chemotherapy.

Evaluation of Antivascular Combretastatin A4 P Efficacy Using Supersonic Shear Imaging Technique of Ectopic Colon Carcinoma CT26

A recent ultrasound imaging technique-shear wave elastography-showed its ability to image and quantify the mechanical properties of biological tissues, such as prostate or liver tissues. In the present study this technique was used to evaluate the relationship among tumor growth, stiffness and reduction of treatment with combretastatin (CA4 P) in allografted colon tumor CT26 in mice. During 12 d, CT26 tumor growth (n = 52) was imaged by ultrasound, and shear modulus was quantified, showing a good correlation between tumor volume and stiffness (r = 0.59). The treatment was initiated at d 12 and monitored every d during 4 d. Following the treatment, the tumor volume had decreased, while the elasticity of the tumor volume remained steady throughout the treatment. After segmentation using the shear modulus map, a detailed analysis showed a decrease in the stiffness after treatment. This reduction in the mechanical properties was shown to correlate with tissue reorganization, particularly, fibrosis and necrosis, assessed by histology.

Correlation of Ultrasound Shear Wave Elastography with Pathological Analysis in a Xenografic Tumour Model

The objective of this study was to evaluate the potential value of ultrasound (US) shear wave elastography (SWE) in assessing the relative change in elastic modulus in colorectal adenocarcinoma xenograft models in vivo and investigate any correlation with histological analysis. We sought to test whether non-invasive evaluation of tissue stiffness is indicative of pathological tumour changes and can be used to monitor therapeutic efficacy. US-SWE was performed in tumour xenografts in 15 NCr nude immunodeficient mice, which were treated with either the cytotoxic drug, Irinotecan, or saline as control. Ten tumours were imaged 48 hours post-treatment and five tumours were imaged for up to five times after treatment. All tumours were harvested for histological analysis and comparison with elasticity measurements. Elastic (Young's) modulus prior to treatment was correlated with tumour volume (r = 0.37, p = 0.008). Irinotecan administration caused significant delay in the tumour growth (p = 0.02) when compared to control, but no significant difference in elastic modulus was detected. Histological analysis revealed a significant correlation between tumour necrosis and elastic modulus (r = -0.73, p = 0.026). SWE measurement provided complimentary information to other imaging modalities and could indicate potential changes in the mechanical properties of tumours, which in turn could be related to the stages of tumour development.