Conventional frequency ultrasonic biomarkers of cancer treatment response in vivo

Apriori prediction of chemotherapy response in locally advanced breast cancer patients using CT imaging and deep learning: transformer versus transfer learning

Objective

Neoadjuvant chemotherapy (NAC) is a key element of treatment for locally advanced breast cancer (LABC). Predicting the response to NAC for patients with Locally Advanced Breast Cancer (LABC) before treatment initiation could be beneficial to optimize therapy, ensuring the administration of effective treatments. The objective of the work here was to develop a predictive model to predict tumor response to NAC for LABC using deep learning networks and computed tomography (CT).

Materials and methods

Several deep learning approaches were investigated including ViT transformer and VGG16, VGG19, ResNet-50, Res-Net-101, Res-Net-152, InceptionV3 and Xception transfer learning networks. These deep learning networks were applied on CT images to assess the response to NAC. Performance was evaluated based on balanced_accuracy, accuracy, sensitivity and specificity classification metrics. A ViT transformer was applied to utilize the attention mechanism in order to increase the weight of important part image which leads to better discrimination between classes.

Results

Amongst the 117 LABC patients studied, 82 (70%) had clinical-pathological response and 35 (30%) had no response to NAC. The ViT transformer obtained the best performance range (accuracy = 71 ± 3% to accuracy = 77 ± 4%, specificity = 86 ± 6% to specificity = 76 ± 3%, sensitivity = 56 ± 4% to sensitivity = 52 ± 4%, and balanced_accuracy=69 ± 3% to balanced_accuracy=69 ± 3%) depending on the split ratio of train-data and test-data. Xception network obtained the second best results (accuracy = 72 ± 4% to accuracy = 65 ± 4, specificity = 81 ± 6% to specificity = 73 ± 3%, sensitivity = 55 ± 4% to sensitivity = 52 ± 5%, and balanced_accuracy = 66 ± 5% to balanced_accuracy = 60 ± 4%). The worst results were obtained using VGG-16 transfer learning network.

Conclusion

Deep learning networks in conjunction with CT imaging are able to predict the tumor response to NAC for patients with LABC prior to start. A ViT transformer could obtain the best performance, which demonstrated the importance of attention mechanism.

Enhanced full-inversion-based ultrasound elastography for evaluating tumor response to neoadjuvant chemotherapy in patients with locally advanced breast cancer

PURPOSE

An enhanced ultrasound elastography technique is proposed for early assessment of locally advanced breast cancer (LABC) response to neoadjuvant chemotherapy (NAC).

METHODS

The proposed elastography technique inputs ultrasound radiofrequency data obtained through tissue quasi-static stimulation and adapts a strain refinement algorithm formulated based on fundamental principles of continuum mechanics, coupled with an iterative inverse finite element method to reconstruct the breast Young's modulus (E) images. The technique was explored for therapy response assessment using data acquired from 25 LABC patients before and at weeks 1, 2, and 4 after the NAC initiation (100 scans). The E ratio of tumor to the surrounding tissue was calculated at different scans and compared to the baseline for each patient. Patients' response to NAC was determined many months later using standard clinical and histopathological criteria.

RESULTS

Reconstructed E ratio changes obtained as early as one week after the NAC onset demonstrate very good separation between the two cohorts of responders and non-responders to NAC. Statistically significant differences were observed in the E ratio changes between the two patient cohorts at weeks 1 to 4 after treatment (p-value < 0.001; statistical power greater than 97%). A significant difference in axial strain ratio changes was observed only at week 4 (p-value = 0.01; statistical power = 76%). No significant difference was observed in tumor size changes at weeks 1, 2 or 4.

CONCLUSION

The proposed elastography technique demonstrates a high potential for chemotherapy response monitoring in LABC patients and superior performance compared to strain imaging.

Quantitative Ultrasound for Evaluation of Tumour Response to Ultrasound-Microbubbles and Hyperthermia

Objectives: Prior study has demonstrated the implementation of quantitative ultrasound (QUS) for determining the therapy response in breast tumour patients. Several QUS parameters quantified from the tumour region showed a significant correlation with the patient's clinical and pathological response. In this study, we aim to identify if there exists such a link between QUS parameters and changes in tumour morphology due to combined ultrasound-stimulated microbubbles (USMB) and hyperthermia (HT) using the breast xenograft model (MDA-MB-231). Method: Tumours grown in the hind leg of severe combined immuno-deficient mice were treated with permutations of USMB and HT. Ultrasound radiofrequency data were collected using a 25 MHz array transducer, from breast tumour-bearing mice prior and post-24-hour treatment. Result: Our result demonstrated an increase in the QUS parameters the mid-band fit and spectral 0-MHz intercept with an increase in HT duration combined with USMB which was found to be reflective of tissue structural changes and cell death detected using haematoxylin and eosin and terminal deoxynucleotidyl transferase dUTP nick end labelling stain. A significant decrease in QUS spectral parameters was observed at an HT duration of 60 minutes, which is possibly due to loss of nuclei by the majority of cells as confirmed using histology analysis. Morphological alterations within the tumour might have contributed to the decrease in backscatter parameters. Conclusion: The work here uses the QUS technique to assess the efficacy of cancer therapy and demonstrates that the changes in ultrasound backscatters mirrored changes in tissue morphology.

Implementation of Non-Invasive Quantitative Ultrasound in Clinical Cancer Imaging

Quantitative ultrasound (QUS) is a non-invasive novel technique that allows treatment response monitoring. Studies have shown that QUS backscatter variables strongly correlate with changes observed microscopically. Increases in cell death result in significant alterations in ultrasound backscatter parameters. In particular, the parameters related to scatterer size and scatterer concentration tend to increase in relation to cell death. The use of QUS in monitoring tumor response has been discussed in several preclinical and clinical studies. Most of the preclinical studies have utilized QUS for evaluating cell death response by differentiating between viable cells and dead cells. In addition, clinical studies have incorporated QUS mostly for tissue characterization, including classifying benign versus malignant breast lesions, as well as responder versus non-responder patients. In this review, we highlight some of the important findings of previous preclinical and clinical studies and expand the applicability and therapeutic benefits of QUS in clinical settings. We summarized some recent clinical research advances in ultrasound-based radiomics analysis for monitoring and predicting treatment response and characterizing benign and malignant breast lesions. We also discuss current challenges, limitations, and future prospects of QUS-radiomics.

Low-frequency ultrasound irradiation increases paclitaxel-induced sarcoma cells apoptosis and facilitates the transmembrane delivery of drugs

Sarcoma is a malignant tumor derived from interstitial tissues and requires comprehensive treatment including chemotherapy. Paclitaxel (PTX) is an active agent against sarcoma, but its effect is not sufficiently acceptable and needs to be improved. Low-frequency ultrasound (LFU) has been documented to improve the efficacy of drugs by inducing reversible changes in membrane permeability; however, the effects of the combined use of LFU and PTX for sarcoma tumors remain unclear and warrant further investigation. We investigated the effects of 30 kHz LFU treatment combined with PTX on sarcoma cells A-204 and HT-1080 by analyzing in vitro apoptosis and cell growth inhibition rates, and determined their antitumor effects by examining tumor weights with or without LFU in the S180 sarcoma xenograft model. Drug concentrations in the subcutaneous tumors were measured using high performance liquid chromatography (HPLC). LFU combined with PTX significantly induced cell apoptosis, and blocked the cell cycle of sarcoma cells in G2/M phase, and furthermore, inhibited the activation of JAK2/STAT3 signaling pathway. Meanwhile, LFU combined with PTX inhibited the expression of PD-L1 in vitro, suggesting the potential of enhanced antitumor immunity by this treatment. LFU combined with PTX significantly inhibited the growth of S180 tumors transplanted subcutaneously in Institute of Cancer Research (ICR) mice, and its enhanced effect may be associated with increased local concentrations of PTX in tumor tissues in vivo, with no significant adverse subsequences on body weight observed. We conclude that the combination of LFU and PTX has synergistic antitumor effects and is a candidate for subcutaneous treatment of sarcoma by further increasing the intracellular concentration of PTX.

Deep learning of quantitative ultrasound multi-parametric images at pre-treatment to predict breast cancer response to chemotherapy

In this study, a novel deep learning-based methodology was investigated to predict breast cancer response to neo-adjuvant chemotherapy (NAC) using the quantitative ultrasound (QUS) multi-parametric imaging at pre-treatment. QUS multi-parametric images of breast tumors were generated using the data acquired from 181 patients diagnosed with locally advanced breast cancer and planned for NAC followed by surgery. The ground truth response to NAC was identified for each patient after the surgery using the standard clinical and pathological criteria. Two deep convolutional neural network (DCNN) architectures including the residual network and residual attention network (RAN) were explored for extracting optimal feature maps from the parametric images, with a fully connected network for response prediction. In different experiments, the features maps were derived from the tumor core only, as well as the core and its margin. Evaluation results on an independent test set demonstrate that the developed model with the RAN architecture to extract feature maps from the expanded parametric images of the tumor core and margin had the best performance in response prediction with an accuracy of 88% and an area under the receiver operating characteristic curve of 0.86. Ten-year survival analyses indicate statistically significant differences between the survival of the responders and non-responders identified based on the model prediction at pre-treatment and the standard criteria at post-treatment. The results of this study demonstrate the promising capability of DCNNs with attention mechanisms in predicting breast cancer response to NAC prior to the start of treatment using QUS multi-parametric images.

Ultrasound-Mediated Cancer Therapeutics Delivery using Micelles and Liposomes: A Review

BACKGROUND

Existing cancer treatment methods have many undesirable side effects that greatly reduce the quality of life of cancer patients.

OBJECTIVE

This review will focus on the use of ultrasound-responsive liposomes and polymeric micelles in cancer therapy.

METHODS

This review presents a survey of the literature regarding ultrasound-triggered micelles and liposomes using articles recently published in various journals, as well as some new patents in this field.

RESULTS

Nanoparticles have proven promising as cancer theranostic tools. Nanoparticles are selective in nature, have reduced toxicity, and controllable drug release patterns making them ideal carriers for anticancer drugs. Numerous nanocarriers have been designed to combat malignancies, including liposomes, micelles, dendrimers, solid nanoparticles, quantum dots, gold nanoparticles, and, more recently, metal-organic frameworks. The temporal and spatial release of therapeutic agents from these nanostructures can be controlled using internal and external triggers, including pH, enzymes, redox, temperature, magnetic and electromagnetic waves, and ultrasound. Ultrasound is an attractive modality because it is non-invasive, can be focused on the diseased site, and has a synergistic effect with anticancer drugs.

CONCLUSION

The functionalization of micellar and liposomal surfaces with targeting moieties and the use of ultrasound as a triggering mechanism can help improve the selectivity and enable the spatiotemporal control of drug release from nanocarriers.

Characterizing intra-tumor regions on quantitative ultrasound parametric images to predict breast cancer response to chemotherapy at pre-treatment

The efficacy of quantitative ultrasound (QUS) multi-parametric imaging in conjunction with unsupervised classification algorithms was investigated for the first time in characterizing intra-tumor regions to predict breast tumor response to chemotherapy before the start of treatment. QUS multi-parametric images of breast tumors were generated using the ultrasound radiofrequency data acquired from 181 patients diagnosed with locally advanced breast cancer and planned for neo-adjuvant chemotherapy followed by surgery. A hidden Markov random field (HMRF) expectation maximization (EM) algorithm was applied to identify distinct intra-tumor regions on QUS multi-parametric images. Several features were extracted from the segmented intra-tumor regions and tumor margin on different parametric images. A multi-step feature selection procedure was applied to construct a QUS biomarker consisting of four features for response prediction. Evaluation results on an independent test set indicated that the developed biomarker coupled with a decision tree model with adaptive boosting (AdaBoost) as the classifier could predict the treatment response of patient at pre-treatment with an accuracy of 85.4% and an area under the receiver operating characteristic (ROC) curve (AUC) of 0.89. In comparison, the biomarkers consisted of the features derived from the entire tumor core (without consideration of the intra-tumor regions), and the entire tumor core and the tumor margin could predict the treatment response of patients with an accuracy of 74.5% and 76.4%, and an AUC of 0.79 and 0.76, respectively. Standard clinical features could predict the therapy response with an accuracy of 69.1% and an AUC of 0.6. Long-term survival analyses indicated that the patients predicted by the developed model as responders had a significantly better survival compared to the non-responders. Similar findings were observed for the two response cohorts identified at post-treatment based on standard clinical and pathological criteria. The results obtained in this study demonstrated the potential of QUS multi-parametric imaging integrated with unsupervised learning methods in identifying distinct intra-tumor regions in breast cancer to characterize its responsiveness to chemotherapy prior to the start of treatment.

Prediction of chemotherapy response in breast cancer patients at pre-treatment using second derivative texture of CT images and machine learning

•

Textural and second derivative textural features of CT images can be used in conjunction with machine learning models to predict breast cancer response to chemotherapy prior to the start of treatment.

•

The proposed predictive model separates the patients at pre-treatment into two cohorts (responders/non-responders) with significantly different survival.

•

The proposed methodology is a step forward towards the precision oncology paradigm for breast cancer patients.

Although neoadjuvant chemotherapy (NAC) is a crucial component of treatment for locally advanced breast cancer (LABC), only about 70% of patients respond to it. Effective adjustment of NAC for individual patients can significantly improve survival rates of those resistant to standard regimens. Thus, the early prediction of NAC outcome is of great importance in facilitating a personalized paradigm for breast cancer therapeutics. In this study, quantitative computed tomography (qCT) parametric imaging in conjunction with machine learning techniques were investigated to predict LABC tumor response to NAC. Textural and second derivative textural (SDT) features of CT images of 72 patients diagnosed with LABC were analysed before the initiation of NAC to quantify intra-tumor heterogeneity. These quantitative features were processed through a correlation-based feature reduction followed by a sequential feature selection with a bootstrap 0.632+ area under the receiver operating characteristic (ROC) curve (AUC0.632+) criterion. The best feature subset consisted of a combination of one textural and three SDT features. Using these features, an AdaBoost decision tree could predict the patient response with a cross-validated AUC0.632+ accuracy, sensitivity and specificity of 0.88, 85%, 88% and 75%, respectively. This study demonstrates, for the first time, that a combination of textural and SDT features of CT images can be used to predict breast cancer response NAC prior to the start of treatment which can potentially facilitate early therapy adjustments.

A Brain Tumor Segmentation Framework Based on Outlier Detection Using One-Class Support Vector Machine

Accurate segmentation of brain tumors is a challenging task and also a crucial step in diagnosis and treatment planning for cancer patients. Magnetic resonance imaging (MRI) is the standard imaging modality for detection, characterization, treatment planning and outcome evaluation of brain tumors. MRI scans are usually acquired at multiple sessions before and after the treatment. An automatic segmentation framework is highly desirable to segment brain tumors in MR images as it streamlines the image-guided radiation therapy workflow considerably. Automatic segmentation of brain tumors also facilitates an incremental development of data-driven systems for therapy outcome prediction based on radiomics analysis. In this study, an outlier-detection-based segmentation framework is proposed to delineate brain tumors in magnetic resonance (MR) images automatically. The proposed method considers the tumor and edema pixels in an MR image as outliers compared to the pixels associated with the healthy tissue. The framework generates two outlier masks using independent one-class support vector machines that operate on post-contrast T1-weighted (T1w) and T2-weighted-fluid-attenuation-inversion-recovery (T2-FLAIR) images. The outlier masks are subsequently refined and fused using a number of morphological and logical operators to estimate a tumor mask for each image slice. The framework was constructed and evaluated using the MRI data acquired from 35 and 5 patients with brain metastasis, respectively. The obtained results demonstrated an average Dice similarity coefficient and Hausdorff distance of 0.84 ± 0.06 and 1.85 ± 0.48 mm, respectively, between the manual (ground truth) and automatic tumor contours, on the independent test set.

Optimization of Ultrasound Backscatter Spectroscopy to Assess Neurotoxic Effects of Anesthesia in the Newborn Non-human Primate Brain

Studies in animal models have revealed that long exposures to anesthetics can induce apoptosis in the newborn and young developing brain. These effects have not been confirmed in humans because of the lack of a non-invasive, practical in vivo imaging tool with the ability to detect these changes. Following the successful use of ultrasound backscatter spectroscopy (UBS) to monitor in vivo cell death in breast tumors, we aimed to use UBS to assess the neurotoxicity of the anesthetic sevoflurane (SEVO) in a non-human primate (NHP) model. Sixteen 2- to 7-day-old rhesus macaques were exposed for 5 h to SEVO. Ultrasound scanning was done with a phased array transducer on a clinical ultrasound scanner operated at 10 MHz. Data consisting of 10-15 frames of radiofrequency (RF) echo signals from coronal views of the thalamus were obtained 0.5 and 6.0 h after initiating exposure. The UBS parameter "effective scatterer size" (ESS) was estimated by fitting a scattering form factor (FF) model to the FF measured from RF echo signals. The approach involved analyzing the frequency dependence of the measured FF to characterize scattering sources and selecting the FF model based on a χ² goodness-of-fit criterion. To assess data quality, a rigorous acceptance criterion based on the analysis of prevalence of diffuse scattering (an assumption in the estimation of ESS) was established. ESS changes after exposure to SEVO were compared with changes in a control group of five primates for which ultrasound data were acquired at 0 and 10 min (no apoptosis expected). Over the entire data set, the average measured FF at 0.5 and 6.0 h monotonically decreased with frequency, justifying fitting a single FF over the analysis bandwidth. χ² values of a (inhomogeneous continuum) Gaussian FF model were one-fifth those of the discrete fluid sphere model, suggesting that a continuum scatterer model better represents ultrasound scattering in the young rhesus brain. After application of the data quality criterion, only 5 of 16 subjects from the apoptotic group and 5 of 5 subjects from the control group fulfilled the acceptance criteria. All subjects in the apoptotic group that passed the acceptance criterion exhibited a significant ESS reduction at 6.0 h. These changes (-6.4%, 95% Interquartile Range: -14.3% to -3.3%) were larger than those in the control group (-0.8%, 95% Interquartile Range: -2.0% to 1.5%]). Data with a low prevalence of diffuse scattering corresponded to possibly biased results. Thus, ESS has the potential to detect changes in brain microstructure related to anesthesia-induced apoptosis.

A priori prediction of tumour response to neoadjuvant chemotherapy in breast cancer patients using quantitative CT and machine learning

Response to Neoadjuvant chemotherapy (NAC) has demonstrated a high correlation to survival in locally advanced breast cancer (LABC) patients. An early prediction of responsiveness to NAC could facilitate treatment adjustments on an individual patient basis that would be expected to improve treatment outcomes and patient survival. This study investigated, for the first time, the efficacy of quantitative computed tomography (qCT) parametric imaging to characterize intra-tumour heterogeneity and its application in predicting tumour response to NAC in LABC patients. Textural analyses were performed on CT images acquired from 72 patients before the start of chemotherapy to determine quantitative features of intra-tumour heterogeneity. The best feature subset for response prediction was selected through a sequential feature selection with bootstrap 0.632 + area under the receiver operating characteristic (ROC) curve (\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\mathrm{A}\mathrm{U}\mathrm{C}}_{0.632+}$$\end{document}AUC0.632+) as a performance criterion. Several classifiers were evaluated for response prediction using the selected feature subset. Amongst the applied classifiers an Adaboost decision tree provided the best results with cross-validated \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\mathrm{A}\mathrm{U}\mathrm{C}}_{0.632+}$$\end{document}AUC0.632+, accuracy, sensitivity and specificity of 0.89, 84%, 80% and 88%, respectively. The promising results obtained in this study demonstrate the potential of the proposed biomarkers to be used as predictors of LABC tumour response to NAC prior to the start of treatment.

Mean Scatterer Spacing Estimation Using Cepstrum-Based Continuous Wavelet Transform

The goal of this study was to develop an ultrasound (US) scatterer spacing estimation method using an enhanced cepstral analysis based on continuous wavelet transforms (CWTs). Simulations of backscattering media containing periodic and quasi-periodic scatterers were carried out to test the developed algorithm. Experimental data from HT-29 pellets and in vivo PC3 tumors were then used to estimate the mean scatterer spacing. For simulated media containing quasi-periodic scatterers at 1-mm and 100- [Formula: see text] spacing with 5% positional variation, the developed algorithm yielded a spacing estimation error of ~1% for 25- and 55-MHz US pulses. The mean scatterer spacing of HT-29 cell pellets (31.97 [Formula: see text]) was within 3% of the spacing obtained from histology and agreed with the predicted spacing from simulations based on the same pellets for both frequencies. The agreement extended to in vivo PC3 tumors estimation of the spacing with a variance of 1.68% between the spacing derived from the tumor histology and the application of the CWT to the experimental results. The developed technique outperformed the traditional cepstral methods as it can detect nonprominent peaks from quasi-random scatterer configurations. This work can be potentially used to detect morphological tissue changes during normal development or disease treatment.

Quantitative Ultrasound Monitoring of Breast Tumour Response to Neoadjuvant Chemotherapy: Comparison of Results Among Clinical Scanners

Quantitative ultrasound (QUS) techniques have been demonstrated to detect cell death in vitro and in vivo. Recently, multi-feature classification models have been incorporated into QUS texture-feature analysis methods to increase further the sensitivity and specificity of detecting treatment response in locally advanced breast cancer patients. To effectively incorporate these analytic methods into clinical applications, QUS and texture-feature estimations should be independent of data acquisition systems. The study here investigated the consistencies of QUS and texture-feature estimation techniques relative to several factors. These included the ultrasound system properties, the effects of tissue heterogeneity and the effects of these factors on the monitoring of response to neoadjuvant chemotherapy. Specifically, tumour-response-detection performance based on QUS and texture parameters using two clinical ultrasound systems was compared. Observed variations in data between the systems were small and the results exhibited good agreement in tumour response predictions obtained from both ultrasound systems. The results obtained in this study suggest that tissue heterogeneity was a dominant feature in the parameters measured with the two different ultrasound systems; whereas differences in ultrasound system beam properties only exhibited a minor impact on texture features. The McNemar statistical test performed on tumour response prediction results from the two systems did not reveal significant differences. Overall, the results in this study demonstrate the potential to achieve reliable and consistent QUS and texture-based analyses across different ultrasound imaging platforms.

An Automatic Framework for Segmentation of Brain Tumours at Follow-up Scans after Radiation Therapy

Brain metastasis is the most common intracranial malignancy with a poor overall survival (OS) after treatment. The standard stereotactic radiation therapy (SRT) planning procedure for brain metastasis requires delineating the tumour volume on magnetic resonance (MR) images. MR images are also acquired at multiple follow-up scans after SRT to monitor the treatment outcome through measuring changes in the physical dimensions of the tumour. Such measurements require manual segmentation of the tumour volume on multiple slices of several follow-up images which is tedious and impedes the SRT evaluation work flow considerably. In this study, an automatic framework was proposed to segment the tumour volume on longitudinal MR images acquired at standard follow-up scans after SRT. The multi-step segmentation framework was based on region growing and morphological snakes models that applied the standard SRT planning tumour contour as a basis to approximate the tumour shape and location at each follow-up scan for an accurate automatic segmentation of tumour volume. The framework was evaluated using the MR imaging data acquired from five patients prior to and at three follow-up scans after SRT. The preliminary results indicated that the Dice similarity coefficient between the ground truth tumour masks and their automatically segmented counterparts ranged between 0.84 and 0.90, while the average Dice coefficient for all the follow-up scans was 0.88. The results obtained implied a good potential of the proposed framework for being incorporated into the SRT treatment planning and evaluation systems as well as outcome prediction models.

Quantitative MRI Biomarkers of Stereotactic Radiotherapy Outcome in Brain Metastasis

About 20-40% of cancer patients develop brain metastases, causing significant morbidity and mortality. Stereotactic radiation treatment is an established option that delivers high dose radiation to the target while sparing the surrounding normal tissue. However, up to 20% of metastatic brain tumours progress despite stereotactic treatment, and it can take months before it is evident on follow-up imaging. An early predictor of radiation therapy outcome in terms of tumour local failure (LF) is crucial, and can facilitate treatment adjustments or allow for early salvage treatment. In this study, an MR-based radiomics framework was proposed to derive and investigate quantitative MRI (qMRI) biomarkers for the outcome of LF in brain metastasis patients treated with hypo-fractionated stereotactic radiation therapy (SRT). The qMRI biomarkers were constructed through a multi-step feature extraction/reduction/selection framework using the conventional MR imaging data acquired from 100 patients (133 lesions), and were applied in conjunction with machine learning techniques for outcome prediction and risk assessment. The results indicated that the majority of the features in the optimal qMRI biomarkers characterize the heterogeneity in the surrounding regions of tumour including edema and tumour/lesion margins. The optimal qMRI biomarker consisted of five features that predict the outcome of LF with an area under the curve (AUC) of 0.79, and a cross-validated sensitivity and specificity of 81% and 79%, respectively. The Kaplan-Meier analyses showed a statistically significant difference in local control (p-value < 0.0001) and overall survival (p = 0.01). Findings from this study are a step towards using qMRI for early prediction of local failure in brain metastasis patients treated with SRT. This may facilitate early adjustments in treatment, such as surgical resection or salvage radiation, that can potentially improve treatment outcomes. Investigations on larger cohorts of patients are, however, required for further validation of the technique.

A priori prediction of breast tumour response to chemotherapy using quantitative ultrasound imaging and artificial neural networks

We demonstrate the clinical utility of combining quantitative ultrasound (QUS) imaging of the breast with an artificial neural network (ANN) classifier to predict the response of breast cancer patients to neoadjuvant chemotherapy (NAC) administration prior to the start of treatment.

Using a 6 MHz ultrasound system, radiofrequency (RF) ultrasound data were acquired from 100 patients with biopsy-confirmed locally advanced breast cancer prior to the start of NAC. Quantitative ultrasound mean parameter intensity and texture features were computed from the tumour core and margin, and were compared to the clinical/pathological response and 5-year recurrence-free survival (RFS) of patients. A multi-parametric QUS model in conjunction with an ANN classifier predicted patient response with 96 ± 6% accuracy, and a 0.96 ± 0.08 area under the receiver operating characteristic curve (AUC), compared to 65 ± 10 % accuracy and 0.67 ± 0.14 AUC achieved using a K-Nearest Neighbour (KNN) algorithm. A separate ANN model predicted patient RFS with 85 ± 7% accuracy, and a 0.89 ± 0.11 AUC, whereas the KNN methodology achieved a 58 ± 6 % accuracy and a 0.64 ± 0.09 AUC.

The application of ANN for classifying patient response based on tumour QUS features performs well in terms of predicting response to chemotherapy. The findings here provide a framework for developing personalized a priori chemotherapy selection for patients that are candidates for NAC, potentially resulting in improved patient treatment outcomes and prognosis.

Quantitative ultrasound and apoptotic death in the neonatal primate brain

Apoptosis is triggered in the developing mammalian brain by sedative, anesthetic or antiepileptic drugs during late gestation and early life. Whether human children are vulnerable to this toxicity mechanism remains unknown, as there are no imaging techniques to capture it. Apoptosis is characterized by distinct structural features, which affect the way damaged tissue scatters ultrasound compared to healthy tissue. We evaluated whether apoptosis, triggered by the anesthetic sevoflurane in the brains of neonatal rhesus macaques, can be detected using quantitative ultrasound (QUS). Neonatal (n = 15) rhesus macaques underwent 5 h of sevoflurane anesthesia. QUS images were obtained through the sagittal suture at 0.5 and 6 h. Brains were collected at 8 h and examined immunohistochemically to analyze apoptotic neuronal and oligodendroglial death. Significant apoptosis was detected in white and gray matter throughout the brain, including the thalamus. We measured a change in the effective scatterer size (ESS), a QUS biomarker derived from ultrasound echo signals obtained with clinical scanners, after sevoflurane-anesthesia in the thalamus. Although initial inclusion of all measurements did not reveal a significant correlation, when outliers were excluded, the change in the ESS between the pre- and post-anesthesia measurements correlated strongly and proportionally with the severity of apoptotic death. We report for the first time in vivo changes in QUS parameters, which may reflect severity of apoptosis in the brains of infant nonhuman primates. These findings suggest that QUS may enable in vivo studies of apoptosis in the brains of human infants following exposure to anesthetics, antiepileptics and other brain injury mechanisms.

Computer-Aided Theragnosis Based on Tumor Volumetric Information in Breast Cancer

OBJECTIVE

A computer-assisted technology has recently been proposed for the assessment of therapeutic responses to neoadjuvant chemotherapy in patients with locally advanced breast cancer (LABC). The system, however, extracted features from individual scans in a tumor irrespective of its relation to the other scans of the same patient, ignoring the volumetric information. This study addresses this problem by introducing a novel engineered texton-based method in order to account for volumetric information in the design of textural descriptors to represent tumor scans.

METHODS

A noninvasive computer-aided-theragnosis (CAT) system was developed by employing multiparametric QUS spectral and backscatter coefficient maps. The proceeding was composed of two subdictionaries: one built on the "pretreatment" and another on "week " scans, where was 1, 4, or 8. The learned dictionary of each patient was subsequently used to compute the model (histogram of textons) for each scan of the patient. Advanced machine learning techniques including a kernel-based dissimilarity measure to estimate the distances between "pretreatment" and "mid-treatment" scans as an indication of treatment effectiveness, learning from imbalanced data, and supervised learning were subsequently employed on the texton-based features.

RESULTS

The performance of the CAT system was tested using statistical tests of significance and leave-one-subject-out (LOSO) classification on 56 LABC patients. The proposed texton-based CAT system indicated significant differences in changes between the responding and nonresponding patient populations and achieved high accuracy, sensitivity, and specificity in discriminating between the two patient groups early after the start of treatment, i.e., on weeks 1 and 4 of several months of treatment. Specifically, the CAT system achieved the area under curve of 0.81, 0.83, and 0.85 on weeks 1, 4, and 8, respectively.

CONCLUSION

The proposed texton-based CAT system accounted for the volumetric information in "pretreatment" and "mid-treatment" scans of each patient. It was demonstrated that this attribute of the CAT system could boost its performance compared to the cases that the features were extracted from solely individual scans.

Chemical exchange saturation transfer MRI to assess cell death in breast cancer xenografts at 7T

Purpose

Detecting cell death and predicting tumor response early in a course of chemotherapy could help optimize treatment regimens and improve clinical outcomes. Chemical exchange saturation transfer (CEST) MRI was investigated in vivo to study properties that may be able to detect cancer death.

Results

Using a magnetization transfer ratio (MTR) cutoff of 0.12 at 1.8 ppm was able to differentiate between viable tumor and cell death regions. Comparison of MTR values at this frequency showed significant differences (p < 0.0001) between viable tumor and cell death regions, matching patterns seen on histology. Using this cutoff, the mean increase in cell death index (± standard error of the mean) after chemotherapy was 4 ± 4%, 10% ± 7%, 10 ± 8%, and 4 ± 9% at 4, 8, 12, and 24 h, respectively.

Conclusions

CEST MRI can detect cell death in MDA-231 xenografts but further work is needed to characterize the clinical applications of this finding. Maximum response to chemotherapy occurred at 8-12 h after chemotherapy injection in this in vivo tumor model.

Materials and Methods

Breast cancer xenografts (MDA-MB-231) were scanned using 7 T MRI before and after chemotherapy. As a measure of CEST effect at 0.5 µT saturation amplitude, MTR values at frequency offsets of 1.8 and -3.3 ppm were evaluated. CEST signals after chemotherapy treatment were compared to cell-death histopathology of tumors.

Ultrasonic histogram assessment of early response to concurrent chemo-radiotherapy in patients with locally advanced cervical cancer: a feasibility study

PURPOSE

To monitor early response for locally advanced cervical cancers undergoing concurrent chemo-radiotherapy (CCRT) by ultrasonic histogram.

METHODS

B-mode ultrasound examinations were performed at 4 time points in thirty-four patients during CCRT. Six ultrasonic histogram parameters were used to assess the echogenicity, homogeneity and heterogeneity of tumors.

RESULTS

I_peak increased rapidly since the first week after therapy initiation, whereas W_low, W_high and A_high changed significantly at the second week. The average ultrasonic histogram progressively moved toward the right and converted into more symmetrical shape.

CONCLUSION

Ultrasonic histogram could be served as a potential marker to monitor early response during CCRT.

Tumour Vascular Shutdown and Cell Death Following Ultrasound-Microbubble Enhanced Radiation Therapy

High-dose radiotherapy effects are regulated by acute tumour endothelial cell death followed by rapid tumour cell death instead of canonical DNA break damage. Pre-treatment with ultrasound-stimulated microbubbles (USMB) has enabled higher-dose radiation effects with conventional radiation doses. This study aimed to confirm acute and longitudinal relationships between vascular shutdown and tumour cell death following radiation and USMB in a wild type murine fibrosarcoma model using in vivo imaging.

Methods: Tumour xenografts were treated with single radiation doses of 2 or 8 Gy alone, or in combination with low-/high-concentration USMB. Vascular changes and tumour cell death were evaluated at 3, 24 and 72 h following therapy, using high-frequency 3D power Doppler and quantitative ultrasound spectroscopy (QUS) methods, respectively. Staining using in situ end labelling (ISEL) and cluster of differentiation 31 (CD31) of tumour sections were used to assess cell death and vascular distributions, respectively, as gold standard histological methods.

Results: Results indicated a decrease in the power Doppler signal of up to 50%, and an increase of more than 5 dBr in cell-death linked QUS parameters at 24 h for tumours treated with combined USMB and radiotherapy. Power Doppler and quantitative ultrasound results were significantly correlated with CD31 and ISEL staining results (p < 0.05), respectively. Moreover, a relationship was found between ultrasound power Doppler and QUS results, as well as between micro-vascular densities (CD31) and the percentage of cell death (ISEL) (R² 0.5-0.9).

Conclusions: This study demonstrated, for the first time, the link between acute vascular shutdown and acute tumour cell death using in vivo longitudinal imaging, contributing to the development of theoretical models that incorporate vascular effects in radiation therapy. Overall, this study paves the way for theranostic use of ultrasound in radiation oncology as a diagnostic modality to characterize vascular and tumour response effects simultaneously, as well as a therapeutic modality to complement radiation therapy.

Ultrasonic RF time series for early assessment of the tumor response to chemotherapy

Ultrasound radio-frequency (RF) time series have been shown to carry tissue typing information. To evaluate the potential of RF time series for early prediction of tumor response to chemotherapy, 50MCF-7 breast cancer-bearing nude mice were randomized to receive cisplatin and paclitaxel (treatment group; n = 26) or sterile saline (control group; n = 24). Sequential ultrasound imaging was performed on days 0, 3, 6, and 8 of treatment to simultaneously collect B-mode images and RF data. Six RF time series features, slope, intercept, S1, S2, S3, and S4, were extracted during RF data analysis and contrasted with microstructural tumor changes on histopathology. Chemotherapy administration reduced tumor growth relative to control on days 6 and 8. Compared with day 0, intercept, S1, and S2 were increased while slope was decreased on days 3, 6, and 8 in the treatment group. Compared with the control group, intercept, S1, S2, S3, and S4 were increased, and slope was decreased, on days 3, 6, and 8 in the treatment group. Tumor cell density decreased significantly in the latter on day 3. We conclude that ultrasonic RF time series analysis provides a simple way to noninvasively assess the early tumor response to chemotherapy.

Breast-Lesion Characterization using Textural Features of Quantitative Ultrasound Parametric Maps

This study evaluated, for the first time, the efficacy of quantitative ultrasound (QUS) spectral parametric maps in conjunction with texture-analysis techniques to differentiate non-invasively benign versus malignant breast lesions. Ultrasound B-mode images and radiofrequency data were acquired from 78 patients with suspicious breast lesions. QUS spectral-analysis techniques were performed on radiofrequency data to generate parametric maps of mid-band fit, spectral slope, spectral intercept, spacing among scatterers, average scatterer diameter, and average acoustic concentration. Texture-analysis techniques were applied to determine imaging biomarkers consisting of mean, contrast, correlation, energy and homogeneity features of parametric maps. These biomarkers were utilized to classify benign versus malignant lesions with leave-one-patient-out cross-validation. Results were compared to histopathology findings from biopsy specimens and radiology reports on MR images to evaluate the accuracy of technique. Among the biomarkers investigated, one mean-value parameter and 14 textural features demonstrated statistically significant differences (p < 0.05) between the two lesion types. A hybrid biomarker developed using a stepwise feature selection method could classify the legions with a sensitivity of 96%, a specificity of 84%, and an AUC of 0.97. Findings from this study pave the way towards adapting novel QUS-based frameworks for breast cancer screening and rapid diagnosis in clinic.

Chemotherapy-Response Monitoring of Breast Cancer Patients Using Quantitative Ultrasound-Based Intra-Tumour Heterogeneities

Anti-cancer therapies including chemotherapy aim to induce tumour cell death. Cell death introduces alterations in cell morphology and tissue micro-structures that cause measurable changes in tissue echogenicity. This study investigated the effectiveness of quantitative ultrasound (QUS) parametric imaging to characterize intra-tumour heterogeneity and monitor the pathological response of breast cancer to chemotherapy in a large cohort of patients (n = 100). Results demonstrated that QUS imaging can non-invasively monitor pathological response and outcome of breast cancer patients to chemotherapy early following treatment initiation. Specifically, QUS biomarkers quantifying spatial heterogeneities in size, concentration and spacing of acoustic scatterers could predict treatment responses of patients with cross-validated accuracies of 82 ± 0.7%, 86 ± 0.7% and 85 ± 0.9% and areas under the receiver operating characteristic (ROC) curve of 0.75 ± 0.1, 0.80 ± 0.1 and 0.89 ± 0.1 at 1, 4 and 8 weeks after the start of treatment, respectively. The patients classified as responders and non-responders using QUS biomarkers demonstrated significantly different survivals, in good agreement with clinical and pathological endpoints. The results form a basis for using early predictive information on survival-linked patient response to facilitate adapting standard anti-cancer treatments on an individual patient basis.

High frequency ultrasound imaging and simulations of sea urchin oocytes

High frequency ultrasound backscatter signals from sea urchin oocytes were measured using a 40 MHz transducer and compared to numerical simulations. The Faran scattering model was used to calculate the ultrasound scattered from single oocytes in suspension. The urchin oocytes are non-nucleated with uniform size and biomechanical properties; the backscatter from each cell is similar and easy to simulate, unlike typical nucleated mammalian cells. The time domain signal measured from single oocytes in suspension showed two distinct peaks, and the power spectrum was periodic with minima spaced approximately 10 MHz apart. Good agreement to the Faran scattering model was observed. Measurements from tightly packed oocyte cell pellets showed similar periodic features in the power spectra, which was a result of the uniform size and consistent biomechanical properties of the cells. Numerical simulations that calculated the ultrasound scattered from individual oocytes within a three dimensional volume showed good agreement to the measured signals and B-scan images. A cepstral analysis of the signal was used to calculate the size of the cells, which was 78.7 μm (measured) and 81.4 μm (simulated). This work supports the single scattering approximation, where ultrasound is discretely scattered from single cells within a bulk homogeneous sample, and that multiple scattering has a negligible effect. This technique can be applied towards understanding the complex scattering behaviour from heterogeneous tissues.

Low-frequency ultrasound radiosensitization and therapy response monitoring of tumors: an in vivo study

A new framework has been introduced in this paper for tumor radiosensitization and therapy response monitoring using low-frequency ultrasound. Human fibrosarcoma xenografts grown in severe combined immunodeficiency (SCID) mice (n = 108) were treated using ultrasound-stimulated microbubbles at various concentration and exposed to different doses of radiation. Low-frequency ultrasound radiofrequency (RF) data were acquired from tumors prior to and at different times after treatment. Quantitative ultrasound (QUS) techniques were applied to generate spectral parametric maps of tumors. Textural analysis were performed to quantify spatial heterogeneities within QUS parametric maps. A hybrid model was developed using multiple regression analysis to predict extent of histological tumor cell death non-invasively based on QUS spectral and textural biomarkers. Results of immunohistochemistry on excised tumor sections demonstrated increases in cell death with higher concentration of microbubbles and radiation dose. Quantitative ultrasound results indicated changes that paralleled increases in histological cell death. Specifically, the hybrid QUS biomarker demonstrated a good correlation with extent of tumor cell death observed from immunohistochemistry. A linear discriminant analysis applied in conjunction with the receiver operating characteristic (ROC) curve analysis indicated that the hybrid QUS biomarker can classify tumor cell death fractions with an area under the curve of 91.2. The results obtained in this research suggest that low-frequency ultrasound can concurrently be used to enhance radiation therapy and evaluate tumor response to treatment.

A priori Prediction of Neoadjuvant Chemotherapy Response and Survival in Breast Cancer Patients using Quantitative Ultrasound

Quantitative ultrasound (QUS) can probe tissue structure and analyze tumour characteristics. Using a 6-MHz ultrasound system, radiofrequency data were acquired from 56 locally advanced breast cancer patients prior to their neoadjuvant chemotherapy (NAC) and QUS texture features were computed from regions of interest in tumour cores and their margins as potential predictive and prognostic indicators. Breast tumour molecular features were also collected and used for analysis. A multiparametric QUS model was constructed, which demonstrated a response prediction accuracy of 88% and ability to predict patient 5-year survival rates (p = 0.01). QUS features demonstrated superior performance in comparison to molecular markers and the combination of QUS and molecular markers did not improve response prediction. This study demonstrates, for the first time, that non-invasive QUS features in the core and margin of breast tumours can indicate breast cancer response to neoadjuvant chemotherapy (NAC) and predict five-year recurrence-free survival.

Effect of chromatin structure on quantitative ultrasound parameters

High-frequency ultrasound (~20 MHz) techniques were investigated using in vitro and ex vivo models to determine whether alterations in chromatin structure are responsible for ultrasound backscatter changes in biological samples. Acute myeloid leukemia (AML) cells and their isolated nuclei were exposed to various chromatin altering treatments. These included 10 different ionic environments, DNA cleaving and unfolding agents, as well as DNA condensing agents. Raw radiofrequency (RF) data was used to generate quantitative ultrasound parameters from spectral and form factor analyses. Chromatin structure was evaluated using electron microscopy. Results indicated that trends in quantitative ultrasound parameters mirrored trends in biophysical chromatin structure parameters. In general, higher ordered states of chromatin compaction resulted in increases to ultrasound paramaters of midband fit, spectral intercept, and estimated scatterer concentration, while samples with decondensed forms of chromatin followed an opposite trend. Experiments with isolated nuclei demonstrated that chromatin changes alone were sufficient to account for these observations. Experiments with ex vivo samples indicated similar effects of chromatin structure changes. The results obtained in this research provide a mechanistic explanation for ultrasound investigations studying scattering from cells and tissues undergoing biological processes affecting chromatin.

Computer aided prognosis for cell death categorization and prediction in vivo using quantitative ultrasound and machine learning techniques

PURPOSE

At present, a one-size-fits-all approach is typically used for cancer therapy in patients. This is mainly because there is no current imaging-based clinical standard for the early assessment and monitoring of cancer treatment response. Here, the authors have developed, for the first time, a complete computer-aided-prognosis (CAP) system based on multiparametric quantitative ultrasound (QUS) spectroscopy methods in association with texture descriptors and advanced machine learning techniques. This system was used to noninvasively categorize and predict cell death levels in fibrosarcoma mouse tumors treated using ultrasound-stimulated microbubbles as novel endothelial-cell radiosensitizers.

METHODS

Sarcoma xenograft tumor-bearing mice were treated using ultrasound-stimulated microbubbles, alone or in combination with x-ray radiation therapy, as a new antivascular treatment. Therapy effects were assessed at 2-3, 24, and 72 h after treatment using a high-frequency ultrasound. Two-dimensional spectral parametric maps were generated using the power spectra of the raw radiofrequency echo signal. Subsequently, the distances between "pretreatment" and "post-treatment" scans were computed as an indication of treatment efficacy, using a kernel-based metric on textural features extracted from 2D parametric maps. A supervised learning paradigm was used to either categorize cell death levels as low, medium, or high using a classifier, or to "continuously" predict the levels of cell death using a regressor.

RESULTS

The developed CAP system performed at a high level for the classification of cell death levels. The area under curve of the receiver operating characteristic was 0.87 for the classification of cell death levels to both low/medium and medium/high levels. Moreover, the prediction of cell death levels using the proposed CAP system achieved a good correlation (r = 0.68,  p < 0.001) with histological cell death levels as the ground truth. A statistical test of significance between individual treatment groups with the corresponding control group demonstrated that the predicted levels indicated the same significant changes in cell death as those indicated by the ground-truth levels.

CONCLUSIONS

The technology developed in this study addresses a gap in the current standard of care by introducing a quality control step that generates potentially actionable metrics needed to enhance treatment decision-making. The study establishes a noninvasive framework for quantifying levels of cancer treatment response developed preclinically in tumors using QUS imaging in conjunction with machine learning techniques. The framework can potentially facilitate the detection of refractory responses in patients to a certain cancer treatment early on in the course of therapy to enable switching to more efficacious treatments.

High-frequency ultrasound detection of cell death: Spectral differentiation of different forms of cell death in vitro

High frequency quantitative ultrasound techniques were investigated to characterize different forms of cell death in vitro. Suspension-grown acute myeloid leukemia cells were treated to cause apoptosis, oncosis, mitotic arrest, and heat-induced death. Samples were scanned with 20 and 40 MHz ultrasound and assessed histologically in terms of cellular structure. Frequency-domain analysis of 20 MHz ultrasound data demonstrated midband fit changes of 6.0 ± 0.7 dBr, 6.2 ± 1.8 dBr, 4.0 ± 1.0 dBr and −4.6 ± 1.7 dBr after 48-hour cisplatinum-induced apoptosis, 48-hour oncotic decay, 36-hour colchicine-induced mitotic arrest, and heat treatment compared to control, respectively. Trends from 40 MHz ultrasound were similar. Spectral slope changes obtained from 40 MHz ultrasound data were reflective of alterations in cell and nucleus size. Chromatin pyknosis or lysis trends suggested that the density of nuclear material may be responsible for observed changes in ultrasound backscatter. Flow cytometry analysis confirmed the modes of cell death and supported midband fit trends in ultrasound data. Scatterer-size and concentration estimates obtained from a fluid-filled sphere form factor model further corresponded with spectral analysis and histology. Results indicate quantitative ultrasound spectral analysis may be used for probing anti-cancer response and distinguishing various modes of cell death in vitro.

Early detection of chemotherapy-refractory patients by monitoring textural alterations in diffuse optical spectroscopic images

PURPOSE

Changes in textural characteristics of diffuse optical spectroscopic (DOS) functional images, accompanied by alterations in their mean values, are demonstrated here for the first time as early surrogates of ultimate treatment response in locally advanced breast cancer (LABC) patients receiving neoadjuvant chemotherapy (NAC). NAC, as a standard component of treatment for LABC patient, induces measurable heterogeneous changes in tumor metabolism which were evaluated using DOS-based metabolic maps. This study characterizes such inhomogeneous nature of response development, by determining alterations in textural properties of DOS images apparent at early stages of therapy, followed later by gross changes in mean values of these functional metabolic maps.

METHODS

Twelve LABC patients undergoing NAC were scanned before and at four times after treatment initiation, and tomographic DOS images were reconstructed at each time. Ultimate responses of patients were determined clinically and pathologically, based on a reduction in tumor size and assessment of residual tumor cellularity. The mean-value parameters and textural features were extracted from volumetric DOS images for several functional and metabolic parameters prior to the treatment initiation. Changes in these DOS-based biomarkers were also monitored over the course of treatment. The measured biomarkers were applied to differentiate patient responses noninvasively and compared to clinical and pathologic responses.

RESULTS

Responding and nonresponding patients demonstrated different changes in DOS-based textural and mean-value parameters during chemotherapy. Whereas none of the biomarkers measured prior the start of therapy demonstrated a significant difference between the two patient populations, statistically significant differences were observed at week one after treatment initiation using the relative change in contrast/homogeneity of seven functional maps (0.001<p<0.049), and mean value of water content in tissue (p=0.010). The cross-validated sensitivity and specificity of these parameters at week one of therapy ranged between 80%-100% and 67%-100%, respectively. Higher levels of statistically significant differences were exhibited at week four after start of treatment, with cross-validated sensitivities and specificities ranging between 80% and 100% for three textural and three mean-value parameters. The combination of the textural and mean-value parameters in a "hybrid" profile could better separate the two patient populations early on during a course of treatment, with cross-validated sensitivities and specificities of up to 100% (p=0.001).

CONCLUSIONS

The results of this study suggest that alterations in textural characteristics of DOS images, in conjunction with changes in their mean values, can classify noninvasively the ultimate clinical and pathologic response of LABC patients to chemotherapy, as early as one week after start of their treatment. This provides a basis for using DOS imaging as a tool for therapy personalization.

Quantitative ultrasound assessment of breast tumor response to chemotherapy using a multi-parameter approach

PURPOSE

This study demonstrated the ability of quantitative ultrasound (QUS) parameters in providing an early prediction of tumor response to neoadjuvant chemotherapy (NAC) in patients with locally advanced breast cancer (LABC).

METHODS

Using a 6-MHz array transducer, ultrasound radiofrequency (RF) data were collected from 58 LABC patients prior to NAC treatment and at weeks 1, 4, and 8 of their treatment, and prior to surgery. QUS parameters including midband fit (MBF), spectral slope (SS), spectral intercept (SI), spacing among scatterers (SAS), attenuation coefficient estimate (ACE), average scatterer diameter (ASD), and average acoustic concentration (AAC) were determined from the tumor region of interest. Ultrasound data were compared with the ultimate clinical and pathological response of the patient's tumor to treatment and patient recurrence-free survival.

RESULTS

Multi-parameter discriminant analysis using the κ-nearest-neighbor classifier demonstrated that the best response classification could be achieved using the combination of MBF, SS, and SAS, with an accuracy of 60 ± 10% at week 1, 77 ± 8% at week 4 and 75 ± 6% at week 8. Furthermore, when the QUS measurements at each time (week) were combined with pre-treatment (week 0) QUS values, the classification accuracies improved (70 ± 9% at week 1, 80 ± 5% at week 4, and 81 ± 6% at week 8). Finally, the multi-parameter QUS model demonstrated a significant difference in survival rates of responding and non-responding patients at weeks 1 and 4 (p=0.035, and 0.027, respectively).

CONCLUSION

This study demonstrated for the first time, using new parameters tested on relatively large patient cohort and leave-one-out classifier evaluation, that a hybrid QUS biomarker including MBF, SS, and SAS could, with relatively high sensitivity and specificity, detect the response of LABC tumors to NAC as early as after 4 weeks of therapy. The findings of this study also suggested that incorporating pre-treatment QUS parameters of a tumor improved the classification results. This work demonstrated the potential of QUS and machine learning methods for the early assessment of breast tumor response to NAC and providing personalized medicine with regards to the treatment planning of refractory patients.

Assessment of tumor response to radiation and vascular targeting therapy in mice using quantitative ultrasound spectroscopy

PURPOSE

It is now recognized that the tumor vasculature is in part responsible for regulating tumor responses to radiation therapy. However, the extent to which radiation-based vascular damage contributes to tumor cell death remains unknown. In this work, quantitative ultrasound spectroscopy (QUS) methods were used to investigate the acute responses of tumors to radiation-based vascular treatments.

METHODS

Tumor xenografts (MDA-MB-231) were treated with single radiation doses of 2 or 8 Gy alone, or in combination with pharmacological agents that modulate vascular radiosensitivity. The midband fit, the slope, and the 0-MHz intercept QUS parameters were obtained from a linear-regression fit to the averaged power spectrum of frequency-dependent ultrasound backscatter and were used to quantify acute tumor responses following treatment administration. Power spectrums were extracted from raw volumetric radio-frequency ultrasound data obtained before and 24 h following treatment administration. These parameters have previously been correlated to tumor cell death. Staining using in situ end labeling, carbonic anhydrase 9 and cluster of differentiation 31 of tumor sections were used to assess cell death, oxygenation, and vasculature distributions, respectively.

RESULTS

Results indicate a significant midband fit QUS parameter increases of 3.2 ± 0.3 dBr and 5.4 ± 0.5 dBr for tumors treated with 2 and 8 Gy radiation combined with the antiangiogenic agent Sunitinib, respectively. In contrast, tumors treated with radiation alone demonstrated a significant midband fit increase of 4.4 ± 0.3 dBr at 8 Gy only. Preadministration of basic fibroblast growth factor, an endothelial radioprotector, acted to minimize tumor response following single large doses of radiation. Immunohistochemical analysis was in general agreement with QUS findings; an R(2) of 0.9 was observed when quantified cell death was correlated with changes in midband fit.

CONCLUSIONS

Results from QUS analysis presented in this study confirm that acute tumor response is linked to a vascular effect following high doses of radiation therapy. Overall, this is in agreement with previous reports suggesting that acute tumor radiation response is regulated by a vascular-driven response. Data also suggest that Sunitinib may enhance tumor radiosensitivity through a vascular remodeling process, and that QUS may be sensitive to changes in tissue properties associated with vascular remodeling. Finally, the work also demonstrates the ability of QUS methods to monitor response to radiation-based vascular strategies.

Quantitative ultrasound imaging of therapy response in bladder cancer in vivo

Background and Aims

Quantitative ultrasound (QUS) was investigated to monitor bladder cancer treatment response in vivo and to evaluate tumor cell death from combined treatments using ultrasound-stimulated microbubbles and radiation therapy.

Methods

Tumor-bearing mice (n=45), with bladder cancer xenografts (HT- 1376) were exposed to 9 treatment conditions consisting of variable concentrations of ultrasound-stimulated Definity microbubbles [nil, low (1%), high (3%)], combined with single fractionated doses of radiation (0 Gy, 2 Gy, 8 Gy). High frequency (25 MHz) ultrasound was used to collect the raw radiofrequency (RF) data of the backscatter signal from tumors prior to, and 24 hours after treatment in order to obtain QUS parameters. The calculated QUS spectral parameters included the mid-band fit (MBF), and 0-MHz intercept (SI) using a linear regression analysis of the normalized power spectrum.

Results and Conclusions

There were maximal increases in QUS parameters following treatments with high concentration microbubbles combined with 8 Gy radiation: (ΔMBF = +6.41 ± 1.40 (±SD) dBr and SI= + 7.01 ± 1.20 (±SD) dBr. Histological data revealed increased cell death, and a reduction in nuclear size with treatments, which was mirrored by changes in quantitative ultrasound parameters. QUS demonstrated markers to detect treatment effects in bladder tumors in vivo.

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.

Multiparametric monitoring of chemotherapy treatment response in locally advanced breast cancer using quantitative ultrasound and diffuse optical spectroscopy

PURPOSE

This study evaluated pathological response to neoadjuvant chemotherapy using quantitative ultrasound (QUS) and diffuse optical spectroscopy imaging (DOSI) biomarkers in locally advanced breast cancer (LABC).

MATERIALS AND METHODS

The institution's ethics review board approved this study. Subjects (n = 22) gave written informed consent prior to participating. US and DOSI data were acquired, relative to the start of neoadjuvant chemotherapy, at weeks 0, 1, 4, 8 and preoperatively. QUS parameters including the mid-band fit (MBF), 0-MHz intercept (SI), and the spectral slope (SS) were determined from tumor ultrasound data using spectral analysis. In the same patients, DOSI was used to measure parameters relating to tumor hemoglobin and composition. Discriminant analysis and receiver-operating characteristic (ROC) analysis was used to classify clinical and pathological response during treatment and to estimate the area under the curve (AUC). Additionally, multivariate analysis was carried out for pairwise QUS/DOSI parameter combinations using a logistic regression model.

RESULTS

Individual QUS and DOSI parameters, including the (SI), oxy-hemoglobin (HbO2), and total hemoglobin (HbT) were significant markers for response after one week of treatment (p < 0.01). Multivariate (pairwise) combinations increased the sensitivity, specificity and AUC at this time; the SI + HbO2 showed a sensitivity/specificity of 100%, and an AUC of 1.0.

CONCLUSIONS

QUS and DOSI demonstrated potential as coincident markers for treatment response and may potentially facilitate response-guided therapies. Multivariate QUS and DOSI parameters increased the sensitivity and specificity of classifying LABC patients as early as one week after treatment.

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

Our objective was to determine if supersonic shear wave elastography (SSWE) can detect changes in stiffness of a breast cancer model under therapy. A human invasive carcinoma was implanted in 22 mice. Eleven were treated with an anti-angiogenic therapy and 11 with glucose for 24 d. Tumor volume and stiffness were assessed during 2 wk before treatment and 0, 7, 12, 20 and 24 d after the start of therapy using SSWE. Pathology was assessed after 12 and 24 d of treatment. We found that response to therapy was associated with early softening of treated tumors only, resulting in a significant difference from non-treated tumors after 12 d of treatment (p = 0.03). On pathology, large areas of necrosis were observed at 12 d in treated tumors. Although treatment was still effective, treated tumors subsequently stiffened during a second phase of the treatment (days 12-24), with a small amount of necrosis observed on pathology on day 24. In conclusion, SSWE was able to measure changes in the stiffness of tumors in response to anti-cancer treatment. However, stiffness changes associated with good response to treatment may change over time, and increased stiffness may also reflect therapy efficacy.

Quantification of Ultrasonic Scattering Properties of In Vivo Tumor Cell Death in Mouse Models of Breast Cancer

INTRODUCTION: Quantitative ultrasound parameters based on form factor models were investigated as potential biomarkers of cell death in breast tumor (MDA-231) xenografts treated with chemotherapy. METHODS: Ultrasound backscatter radiofrequency data were acquired from MDA-231 breast cancer tumor–bearing mice (n = 20) before and after the administration of chemotherapy drugs at two ultrasound frequencies: 7 MHz and 20 MHz. Radiofrequency spectral analysis involved estimating the backscatter coefficient from regions of interest in the center of the tumor, to which form factor models were fitted, resulting in estimates of average scatterer diameter and average acoustic concentration (AAC). RESULTS: The ∆AAC parameter extracted from the spherical Gaussian model was found to be the most effective cell death biomarker (at the lower frequency range, r² = 0.40). At both frequencies, AAC in the treated tumors increased significantly (P = .026 and .035 at low and high frequencies, respectively) 24 hours after treatment compared with control tumors. Furthermore, stepwise multiple linear regression analysis of the low-frequency data revealed that a multiparameter quantitative ultrasound model was strongly correlated to cell death determined histologically posttreatment (r² = 0.74). CONCLUSION: The Gaussian form factor model–based scattering parameters can potentially be used to track the extent of cell death at clinically relevant frequencies (7 MHz). The 20-MHz results agreed with previous findings in which parameters related to the backscatter intensity (i.e., AAC) increased with cell death. The findings suggested that, in addition to the backscatter coefficient parameter ∆AAC, biological features including tumor heterogeneity and initial tumor volume were important factors in the prediction of cell death response.

Early prediction of therapy responses and outcomes in breast cancer patients using quantitative ultrasound spectral texture

Early alterations in textural characteristics of quantitative ultrasound spectral parametric maps, in conjunction with changes in their mean values, are demonstrated here, for the first time, to be capable of predicting ultimate clinical/pathologic responses of breast cancer patients to chemotherapy. Mechanisms of cell death, induced by chemotherapy within tumor, introduce morphological alterations in cancerous cells, resulting in measurable changes in tissue echogenicity. We have demonstrated that the development of such changes is reflected in early alterations in textural characteristics of quantitative ultrasound spectral parametric maps, followed by consequent changes in their mean values. The spectral/textural biomarkers derived on this basis have been demonstrated as non-invasive surrogates of breast cancer chemotherapy response. Particularly, spectral biomarkers sensitive to the size and concentration of acoustic scatterers could predict treatment response of patients with up to 80% of sensitivity and specificity (p=0.050), after one week within 3-4 months of chemotherapy. However, textural biomarkers characterizing heterogeneities in distribution of acoustic scatterers, could differentiate between treatment responding and non-responding patients with up to 100% sensitivity and 93% specificity (p=0.002). Such early prediction permits offering effective alternatives to standard treatment, or switching to a salvage therapy, for refractory patients.

Quantitative ultrasound spectroscopic imaging for characterization of disease extent in prostate cancer patients

Three-dimensional quantitative ultrasound spectroscopic imaging of prostate was investigated clinically for the noninvasive detection and extent characterization of disease in cancer patients and compared to whole-mount, whole-gland histopathology of radical prostatectomy specimens. Fifteen patients with prostate cancer underwent a volumetric transrectal ultrasound scan before radical prostatectomy. Conventional-frequency (~5MHz) ultrasound images and radiofrequency data were collected from patients. Normalized power spectra were used as the basis of quantitative ultrasound spectroscopy. Specifically, color-coded parametric maps of 0-MHz intercept, midband fit, and spectral slope were computed and used to characterize prostate tissue in ultrasound images. Areas of cancer were identified in whole-mount histopathology specimens, and disease extent was correlated to that estimated from quantitative ultrasound parametric images. Midband fit and 0-MHz intercept parameters were found to be best associated with the presence of disease as located on histopathology whole-mount sections. Obtained results indicated a correlation between disease extent estimated noninvasively based on midband fit parametric images and that identified histopathologically on prostatectomy specimens, with an r(2) value of 0.71 (P<.0001). The 0-MHz intercept parameter demonstrated a lower level of correlation with histopathology. Spectral slope parametric maps offered no discrimination of disease. Multiple regression analysis produced a hybrid disease characterization model (r(2)=0.764, P<.05), implying that the midband fit biomarker had the greatest correlation with the histopathologic extent of disease. This work demonstrates that quantitative ultrasound spectroscopic imaging can be used for detecting prostate cancer and characterizing disease extent noninvasively, with corresponding gross three-dimensional histopathologic correlation.

Review of ultrasound image guidance in external beam radiotherapy: I. Treatment planning and inter-fraction motion management

In modern radiotherapy, verification of the treatment to ensure the target receives the prescribed dose and normal tissues are optimally spared has become essential. Several forms of image guidance are available for this purpose. The most commonly used forms of image guidance are based on kilovolt or megavolt x-ray imaging. Image guidance can also be performed with non-harmful ultrasound (US) waves. This increasingly used technique has the potential to offer both anatomical and functional information.This review presents an overview of the historical and current use of two-dimensional and three-dimensional US imaging for treatment verification in radiotherapy. The US technology and the implementation in the radiotherapy workflow are described. The use of US guidance in the treatment planning process is discussed. The role of US technology in inter-fraction motion monitoring and management is explained, and clinical studies of applications in areas such as the pelvis, abdomen and breast are reviewed. A companion review paper (O'Shea et al 2015 Phys. Med. Biol. submitted) will extensively discuss the use of US imaging for intra-fraction motion quantification and novel applications of US technology to RT.

Noninvasive characterization of locally advanced breast cancer using textural analysis of quantitative ultrasound parametric images

PURPOSE: The identification of tumor pathologic characteristics is an important part of breast cancer diagnosis, prognosis, and treatment planning but currently requires biopsy as its standard. Here, we investigated a noninvasive quantitative ultrasound method for the characterization of breast tumors in terms of their histologic grade, which can be used with clinical diagnostic ultrasound data. METHODS: Tumors of 57 locally advanced breast cancer patients were analyzed as part of this study. Seven quantitative ultrasound parameters were determined from each tumor region from the radiofrequency data, including mid-band fit, spectral slope, 0-MHz intercept, scatterer spacing, attenuation coefficient estimate, average scatterer diameter, and average acoustic concentration. Parametric maps were generated corresponding to the region of interest, from which four textural features, including contrast, energy, homogeneity, and correlation, were determined as further tumor characterization parameters. Data were examined on the basis of tumor subtypes based on histologic grade (grade I versus grade II to III). RESULTS: Linear discriminant analysis of the means of the parametric maps resulted in classification accuracy of 79%. On the other hand, the linear combination of the texture features of the parametric maps resulted in classification accuracy of 82%. Finally, when both the means and textures of the parametric maps were combined, the best classification accuracy was obtained (86%). CONCLUSIONS: Textural characteristics of quantitative ultrasound spectral parametric maps provided discriminant information about different types of breast tumors. The use of texture features significantly improved the results of ultrasonic tumor characterization compared to conventional mean values. Thus, this study suggests that texture-based quantitative ultrasound analysis of in vivo breast tumors can provide complementary diagnostic information about tumor histologic characteristics.

Quantification of heterogeneity as a biomarker in tumor imaging: a systematic review

BACKGROUND

Many techniques are proposed for the quantification of tumor heterogeneity as an imaging biomarker for differentiation between tumor types, tumor grading, response monitoring and outcome prediction. However, in clinical practice these methods are barely used. This study evaluates the reported performance of the described methods and identifies barriers to their implementation in clinical practice.

METHODOLOGY

The Ovid, Embase, and Cochrane Central databases were searched up to 20 September 2013. Heterogeneity analysis methods were classified into four categories, i.e., non-spatial methods (NSM), spatial grey level methods (SGLM), fractal analysis (FA) methods, and filters and transforms (F&T). The performance of the different methods was compared.

PRINCIPAL FINDINGS

Of the 7351 potentially relevant publications, 209 were included. Of these studies, 58% reported the use of NSM, 49% SGLM, 10% FA, and 28% F&T. Differentiation between tumor types, tumor grading and/or outcome prediction was the goal in 87% of the studies. Overall, the reported area under the curve (AUC) ranged from 0.5 to 1 (median 0.87). No relation was found between the performance and the quantification methods used, or between the performance and the imaging modality. A negative correlation was found between the tumor-feature ratio and the AUC, which is presumably caused by overfitting in small datasets. Cross-validation was reported in 63% of the classification studies. Retrospective analyses were conducted in 57% of the studies without a clear description.

CONCLUSIONS

In a research setting, heterogeneity quantification methods can differentiate between tumor types, grade tumors, and predict outcome and monitor treatment effects. To translate these methods to clinical practice, more prospective studies are required that use external datasets for validation: these datasets should be made available to the community to facilitate the development of new and improved methods.

Quantitative ultrasound characterization of tumor cell death: ultrasound-stimulated microbubbles for radiation enhancement

The aim of this study was to assess the efficacy of quantitative ultrasound imaging in characterizing cancer cell death caused by enhanced radiation treatments. This investigation focused on developing this ultrasound modality as an imaging-based non-invasive method that can be used to monitor therapeutic ultrasound and radiation effects. High-frequency (25 MHz) ultrasound was used to image tumor responses caused by ultrasound-stimulated microbubbles in combination with radiation. Human prostate xenografts grown in severe combined immunodeficiency (SCID) mice were treated using 8, 80, or 1000 µL/kg of microbubbles stimulated with ultrasound at 250, 570, or 750 kPa, and exposed to 0, 2, or 8 Gy of radiation. Tumors were imaged prior to treatment and 24 hours after treatment. Spectral analysis of images acquired from treated tumors revealed overall increases in ultrasound backscatter intensity and the spectral intercept parameter. The increase in backscatter intensity compared to the control ranged from 1.9±1.6 dB for the clinical imaging dose of microbubbles (8 µL/kg, 250 kPa, 2 Gy) to 7.0±4.1 dB for the most extreme treatment condition (1000 µL/kg, 750 kPa, 8 Gy). In parallel, in situ end-labelling (ISEL) staining, ceramide, and cyclophilin A staining demonstrated increases in cell death due to DNA fragmentation, ceramide-mediated apoptosis, and release of cyclophilin A as a result of cell membrane permeabilization, respectively. Quantitative ultrasound results indicated changes that paralleled increases in cell death observed from histology analyses supporting its use for non-invasive monitoring of cancer treatment outcomes.

Categorizing extent of tumor cell death response to cancer therapy using quantitative ultrasound spectroscopy and maximum mean discrepancy

Quantitative ultrasound (QUS) spectroscopic techniques in conjunction with maximum mean discrepancy (MMD) have been proposed to detect, and to classify noninvasively the levels of cell death in response to cancer therapy administration in tumor models. Evaluation of xenograft tumor responses to cancer treatments were carried out using conventional-frequency ultrasound at different times after chemotherapy exposure. Ultrasound data were analyzed using spectroscopic techniques and multi-parametric QUS spectral maps were generated. MMD was applied as a distance criterion, measuring alterations in each tumor in response to chemotherapy, and the extent of cell death was classified into less/more than 20% and 40% categories. Statistically significant differences were observed between "pre-" and "post-treatment" groups at different times after chemotherapy exposure, suggesting a high capability of proposed framework for detecting tumor response noninvasively. Promising results were also obtained for categorizing the extent of cell death response in each tumor using the proposed framework, with gold standard histological quantification of cell death as ground truth. The best classification results were obtained using MMD when applied on histograms of QUS parametric maps. In this case, classification accuracies of 84.7% and 88.2% were achieved for categorizing extent of tumor cell death into less/more than 20% and 40%, respectively.