Next Article in Journal
Circ_RNF13 Regulates the Stemness and Chemosensitivity of Colorectal Cancer by Transcriptional Regulation of DDX27 Mediated by TRIM24 Stabilization
Next Article in Special Issue
Dysfunctional Lipid Metabolism—The Basis for How Genetic Abnormalities Express the Phenotype of Aggressive Prostate Cancer
Previous Article in Journal
Approach to Contemporary Risk Assessment, Prevention and Management of Thrombotic Complications in Multiple Myeloma
Previous Article in Special Issue
The Role of Radiomics and AI Technologies in the Segmentation, Detection, and Management of Hepatocellular Carcinoma
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Review

Implementation of Non-Invasive Quantitative Ultrasound in Clinical Cancer Imaging

by
Deepa Sharma
1,2,3,*,†,
Laurentius Oscar Osapoetra
1,2,3,† and
Gregory J. Czarnota
1,2,3,*
1
Imaging Research and Physical Sciences, Sunnybrook Health Sciences Centre, Toronto, ON M4N 3M5, Canada
2
Department of Radiation Oncology, Sunnybrook Health Sciences Centre, Toronto, ON M4N 3M5, Canada
3
Departments of Medical Biophysics and Radiation Oncology, University of Toronto, Toronto, ON M5S, Canada
*
Authors to whom correspondence should be addressed.
These authors contributed equally to this work.
Cancers 2022, 14(24), 6217; https://doi.org/10.3390/cancers14246217
Submission received: 23 November 2022 / Accepted: 2 December 2022 / Published: 16 December 2022
(This article belongs to the Special Issue Feature Reviews in Cancer Therapy)

Abstract

:

Simple Summary

At present, quantitative ultrasound (QUS) is increasingly utilized in cancer imaging. Compared to other imaging modalities that take several days to weeks to assess treatment effectiveness, QUS can provide a rapid treatment evaluation. Its implementation has been documented in characterizing benign versus malignant breast lesions, assessing lymph nodes, predicting and monitoring tumor response, etc. Both preclinical and clinical studies have confirmed that changes in QUS parameters are directly correlated with tissue microstructural alterations. QUS parameters and textural analyses have widely been used to predict and monitor neoadjuvant chemotherapy (NAC) response in locally advanced breast cancer (LABC) patients. Thus, QUS methods have emerged as one of the most useful imaging techniques for the management of several tumor types.

Abstract

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.

1. Introduction

Over the years, multiple technologies have paved the way in the treatment of cancers. Diagnosing cancer is a multi-step, complex process. Diagnostic imaging procedures have made earlier detection of cancer possible. Some of the commonly used imaging techniques include computed tomography (CT), magnetic resonance imaging (MRI), positron emission tomography (PET), ultrasound, and X-ray [1,2]. Quantitative ultrasound (QUS) is one of such kinds that has gained popularity by offering many advantages over other modalities. It is a non-invasive device that is cost-effective, portable, and radiation-free [3]. No incisions or cuts are required during the imaging process, which makes it a safe and easily accessible device. The use of QUS has widely been explored both as a diagnostic and therapeutic imaging modality.
Several preclinical and clinical studies have implemented QUS for cancer prediction and monitoring [4,5,6,7,8,9]. Different parameters can be extracted from linear regression analysis using normalized power spectra of radio frequency (RF) data from ultrasound [10,11,12,13,14,15]. Linear-fit parameters include the mid-band fit (MBF), 0-MHz intercept (SI), and spectral slope (SS). The first two parameters are associated with the amount of backscattering, while the last one is associated with the scatter’s size [10,11,12]. In order to obtain spectral parameters that reflect more about acoustic scatterers, such as their concentration and size, theoretical acoustic scattering models can be fitted to the measured back-scattering function of the sample. This results in average scatterer diameter (ASD)/effective scatterer diameter (ESD)/effective scatterer size (ESS) and average acoustic concentration (AAC)/effective acoustic concentration (EAC)/effective scatterer concentration (ESC) parameters. These parameters (ASD and AAC) can be estimated using the spherical Gaussian scattering model (SGM) and the fluid-filled-sphere model (FFSM) form factor models to the ultrasonic backscatter coefficient (BSC) [13,14,15]. The ASD and AAC parameters correlate with the changes in scatterer sizes and concentrations, respectively. QUS spectral parametric imaging creates maps of quantitative parameters obtained through spectral analysis of ultrasound RF data. These maps reveal tissue intrinsic scattering properties. Tissue microstructures and their scattering properties are distinct between invasive and non-invasive tumors [16], and from responsive to treatment compared to those of non-responsive ones [17]. The QUS spectral parametric images can provide a surrogate delineation of tumor microstructures providing diagnostic [18,19] and prognostic potentials [20,21,22,23,24,25]. QUS spectral parametric maps are created using a sliding window technique. The window extracts a block of RF data for spectral analysis. The kernel is moved to all points in the region of interest (ROI), given some kernel overlap factors. For clinical applications, typically a 2 mm × 2 mm kernel size is used. This size amounts to approximately ten acoustic wavelengths. A Hanning gating function is applied in the axial direction to smooth out the RF segments. The Fast Fourier Transform (FFT) method is utilized to extract the spectral contents of each RF segment in the block. An average power spectrum associated with the kernel is then obtained by averaging RF spectra in the kernel.
A reference phantom technique is utilized to obtain tissue-dependent scattering components from the RF signal [26,27]. The spatially homogeneous reference phantom has well-characterized frequency-dependent back-scattering and attenuation functions. The phantom is scanned using the same ultrasound system and acquisition settings for extracting reference RF data. Corresponding to the RF block from the sample, the reference RF block is acquired from the same location in the reference phantom. The spectral normalization procedure divides the sample spectrum by the reference spectrum. Along with the corresponding attenuation functions and reference backscatter coefficient (BSC), this allows for the estimation of the tissue-dependent scattering function. Parametrization of the frequency-dependent scattering function results in the linear-fit and acoustic scattering parameters useful for diagnostics and prognostics [18,19,20,21,22,23,24,25].
These analyses are repeated for each point in the ROI to come up with parametric images of QUS spectral parameters. From the parametric maps, radiomics features can be extracted that include first-order statistical, morphological, and textural features. These radiomic features are potential imaging biomarkers for diagnosis and prognosis purposes [18,19,20,21,22,23,24,25]. Radiomics texture features rely on the hypothesis that tissue heterogeneity can be quantified through their surrogate spectral parametric maps (Figure 1).

1.1. Ultrasound Imaging of Cell Death in Tumor Response

Studies have confirmed a strong correlation between changes in cell nuclear structure and ultrasound parameters using both low- and high-frequency ultrasound. A study by Banihashemi et al. showed an increase in MBF and SS parameters in a time-dependent manner following photodynamic therapy (PDT) in an in vivo melanoma model. They reported 45% apoptotic cell death at 24 h, that corresponded to maximum SS increase, showing a value of 0.435 ± 0.07 dB/MHz. Similarly, an increase in MBF was observed between 12 to 24 h. At 48 h, around 50% of the cells were seen to have their nuclei disappeared as a result of late stage of apoptosis. The loss of cell nuclei at 48 h resulted in the decrement of ultrasound backscatter parameters [28]. Thus, a direct correlation between cell death and ultrasound parameters was established in this study. A similar outcomes were reported by Sadeghi-Naini et al. in a clinical study conducted with locally advanced breast cancer (LABC) patients receiving chemotherapy. They reported an increase in MBF and SI at week 4 of treatment. Responding patients demonstrated 9.1 ± 1.2 dBr increases in MBF as compared to the non-responder with the MBF value of 1.9 ± 1.1 dBr. Similarly, the SI variable increased in responding patients to 8.9 ± 1.9 dBr compared to 1.6 ± 0.9 dBr observed in non-responding patients. Additionally, histopathology data from the non-responder patient demonstrated a large compact mass in the mastectomy specimen; however, no such mass was observed in the specimen of responder patients. Thus, their result showed that patients responding to treatment showed significant changes in QUS parameters as compared to non-responders [5]. Moreover, parameters such as AAC and ASD have also been found to strongly correspond to cell death. A clinical study reported a substantial increase in AAC and ASD at weeks 1, 4, and 8 in LABC patients that responded to chemotherapy. Maximum increase was observed at week 8. On the contrary, no such increase was seen in non-responding patients [13]. Similar results were also reported with breast cancer xenograft following an exposure to chemotherapy. Higher AAC corresponded to the highest cell death, indicating 60% at 24 h. A strong correlation between AAC and cell death was demonstrated with (R2 SGM = 0.40) [14]. Thus, many endeavors have been made that confirm the utility of QUS and its parameters for monitoring treatment response, as well as for differentiating between clinically responding and non-responding patients. In the following sections, we discuss in detail some of the important findings that used QUS in clinical settings for cancer diagnosis, tissue characterization, treatment prediction and monitoring. The summary of QUS application and its utilization in preclinical and clinical studies is presented in Table 1.

1.2. Clinical Applications of Classification Models Developed from QUS Spectral Parametric Images Using Machine Learning Approaches

Tumors Characterization

Apart from monitoring tumor response and treatment efficacy, the use of QUS in classifying benign versus malignant tumors has also been widely studied [54,61,62]. The first study by Sadeghi-Naini et al. utilized QUS spectral and textural analysis techniques to characterize breast lesions from 78 patients. Intra-lesion heterogeneity within tissue micro-structures was quantified using textural features and average-based mean-value parameters from each parametric map [54]. Later on, Osapoetra et al. expanded this study to a larger cohort by analyzing texture and texture-derivate within the peri-tumoral breast tissue [61]. Both of the studies only utilized the Gray Level Co-Occurrence Matrix (GLCM)-based texture method to characterize breast lesions. Furthermore, another study by Osapoetra et al. incorporated different texture methods, including GLCM, the Gray Level Run Length Matrix (GLRLM), and the Gray Level Size Zone Matrix (GLSZM). They applied these methods to QUS spectral parametric images to characterize breast lesions from the tumor core and tumor margin [62].
A study by Sadeghi-Naini et al. explored different texture features to discriminate between benign versus malignant breast lesions. Based on the size, density, and distribution of acoustic scatterers, the QUS texture parameters can be utilized to quantify intra-lesional heterogeneity, providing characterization and estimation of tissue microstructure. In their study, several texture features, including contrast (CON), correlation (COR), homogeneity (HOM), and energy (ENE), were generated using parametric maps of MBF, SS, SI, SAS, EAC, and ESD. These texture features were then used as biomarkers to classify breast lesions. The result revealed significant differences in benign versus malignant lesions with (p < 0.05) in most of the textural features, including MBF (CON), (COR), (HOM), SS (CON), (COR), (HOM), and SI (CON), (COR), (HOM), SAS (ENE), ESD (CON), (COR), and EAC (HOM), (ENE) [54].
QUS spectral parametric imaging and texture analyses have been applied for the characterization of breast lesions [61,62]. In this task, Osapoetra et al. developed classification models using radiomics features extracted from the parametric images, along with traditional machine learning classifiers, to distinguish malignant lesions from benign breast lesions [61]. In that study, the authors hypothesized that texture-derivate features exhibit significant discriminating power for developing a multi-feature classification model. Their cohort consisted of 204 patients with breast lesions (99 benign and 105 malignant). Figure 2 and Figure 3 show representative B-mode, ASD, AAC, MBF, SS, and SI parametric maps from the benign and malignant groups, respectively. Figure 4 depicts scatter and box plots of representative discriminating features. Texture-derivate features were demonstrated as ones with the most separation between the two groups. They reported the best classification performance of 90% sensitivity, 92% specificity, 91% accuracy, and 0.93 area under the curve (AUC) using features from the tumor core and margin [61].
In another study, different ways to quantify texture images were studied for the task of characterizing breast lesions. Classification models were developed using the GLCM, GLRLM, and GLSZM radiomics features of QUS spectral parametric images. They reported a classification accuracy from 73–91%, depending on tumor and margin attributes and classification algorithms, on a cohort of 193 patients with breast lesions (92 benign and 101 malignant cases) [62].
QUS spectroscopy have also been applied for the characterization of different types of cancer. Rohrbach et al. demonstrated an application of QUS in the non-invasive characterization of prostate cancers using a novel transrectal high-frequency ultrasound system (ExactVu micro-ultrasound Exact Imaging, Markham, ON, Canada). They developed a multi-variate model from ASD, AAC, envelope statistics, and prostate-specific antigen (PSA) features, achieving the best AUC of 0.81 using a linear discriminant classifier [55]. Furthermore, Goundan et al. reported a thyroid nodules (TNs) characterization study [63]. The cohort consisted of 225 TNs (24 malignant and 201 benign) from 208 patients. They developed a classification model using QUS spectral features that include ASD, AAC, MBF, and SI parameters, along with an envelope statistics feature in the Nakagami shape parameter. Their model achieved the best classification performance of 0.857 ± 0.033 receiver operating characteristics (ROC) AUC, which was comparable with the value achieved using the American College of Radiology Thyroid Imaging, Reporting and Data System (ACR TI-RADS) and American Thyroid Association (ATA) risk-stratification systems. Furthermore, they also demonstrated a reduction in the number of unnecessary fine-needle biopsies (FNBs) when the combination of QUS parameters were used separately or with the ACR TI-RADS [63].

1.3. LABC QUS Treatment Response Prediction

QUS spectral parametric imaging and texture analyses have been applied for the clinical application of response characterization to cancer treatments. The underlying hypothesis is that QUS spectral parameters are sensitive enough to detect microscopical changes in tumors. As the tumors respond to cancer therapeutics, these changes can be used as imaging biomarkers for developing classification models. Such models are potentially useful in clinical practice, as they can guide early interventions during the full course of cancer treatment. For example, in the neoadjuvant chemotherapy (NAC) management of LABC, early prediction of non-responder patients will allow clinicians to decide if other therapeutic agents might fare better for the objective of reducing the tumor size prior to the main surgical procedure to remove the LABC tumors.
LABC treatment response prediction models from QUS spectral parametric imaging and texture methods have been reported in a number of recent studies [20,21,22,23,25]. These studies investigated several aspects of the classification model for treatment response prediction that utilize radiomics features from higher-order GLCM-based textural features [21,24]; assessment of potential issues regarding the variability of the radiomics features as a result of different ultrasound systems in the clinic [20]; and model evaluation from multi-institutional data [18,21], utilization of a novel approach that extract radiomics features from intra-tumoral regions obtained through an unsupervised segmentation technique [22], and development of a deep learning-based model [23].
First, Dasgupta et al. studied the classification performance of models developed utilizing higher-order textural features (texture-derivative) [58]. The cohort consisted of 100 patients (83 responders and 17 non-responders) with LABC undergoing NAC. They extracted a pool of radiomics features, including the first-order statistical mean-value, GLCM-based texture, and the proposed GLCM-based texture-derivate features. The best set of features was determined through a sequential feature selection approach to develop classification models to identify non-responders from responders. They reported the best classification performance of 87% sensitivity, 81% specificity, 82% accuracy, and 0.86 AUC [58].
Second, Sannachi et al. addressed the potential issue of variations in the measured radiomics features as they were acquired across different ultrasound systems prevalent in clinics, along with their impacts on model classification performance. A GE-LOGIQ E9 (General Electric Healthcare, Milwaukee, WI, USA) and Ultrasonix-RP ((ULX) Ultrasonix Medical Corp., Richmond, BC, Canada) ultrasound system were utilized. They analyzed radiomics features from a cohort of 24 patients with LABC [20]. A previously developed classification model for treatment response monitoring was also applied to assess the performance of the smaller cohort [54]. They concluded that observed variations in data due to system-specific features were small, and the results of the prediction models were comparable across both ultrasound systems utilized Table 2. They found tissue heterogeneity to be the dominant factor causing variations in the measured radiomics features [20] (Figure 5 and Figure 6).
DiCenzo et al. studied the efficacy of the classification model for therapy response prediction using LABC data from multiple institutions. Their cohort consisted of 82 patients with LABC undergoing NAC (48 responders and 34 non-responders). The patients included were from multiple institutions. They developed classification models using a pool of pre-treatment radiomics features that include first-order statistical mean-value and GLCM-based texture features. They reported the best performance metric of 87% accuracy in predicting non-responders from responders [21]. Subsequently, Osapoetra et al. also assessed the classification model for therapy response prediction, utilizing more pre-treatment radiomics features that include higher-order GLCM-based texture on a similar multi-institutional dataset. The cohort consisted of 74 patients (42 responders and 32 non-responders) with LABC. Higher-order texture features have shown their discriminating power for breast lesions characterization datasets. They reported the best classification performance of 88% sensitivity, 78% specificity, 84% accuracy, and 0.86 AUC [18].
In the typical framework for building a classification model for treatment response prediction using radiomics texture features, the average value or average texture features were estimated from the whole QUS spectral parametric image. Taleghamar et al. proposed a different approach to utilizing spatial information from intra-tumoral parametric maps. Their approach first segmented a discrete number of regions of the tumor from the parametric maps, prior to obtaining average features. They utilized an unsupervised learning approach for segmentation using a hidden Markov random field (HMRF) expectation maximization (EM) algorithm. Their approach allows for the utilization of radiomics features only from a certain region of the intra-tumoral parametric maps. Their cohort consisted of 181 patients (138 responders and 43 non-responders) with LABC undergoing NAC. The radiomics features included first-order statistical mean-value and signal to noise ratio (SNR) from each QUS spectral parametric image at pre-treatment. They reported a classification model using the Adaptive Boosting (AdaBoost) technique, attaining 85.4% accuracy and 0.89 AUC on an independent test set [22].
Deep learning models learn how to represent data with multiple levels of abstraction using multiple processing layers [64]. They learn the features directly from data, in contrast to traditional machine learning involving feature engineering. These methods have set the state-of-the-art performance for tasks in computer vision, natural language processing, and many other domains [39]. This trend motivated a recent effort by Taleghamar et al. in utilizing the deep learning framework to boost the performance of the classification model for treatment response prediction [23]. In contrast to previous studies on deep learning for breast lesions characterization that use B-mode images, the inputs to the network are QUS spectral parametric images [22]. In that study, Taleghamar et al. implemented convolutional neural network (CNN) models for treatment response prediction [23]. Their approach utilized a modified residual network version 101 (ResNet) [65] and a modified residual attention network version 56 (RAN) [66] as the convolutional base. A densely connected classifier was attached on top of the base layer. The networks were trained from scratch using augmented QUS spectral parametric images from their cohort. The cohort consisted of 181 LABC patients (138 responders and 43 non-responders) undergoing NAC. They reported the best performance of 88% accuracy and 0.86 AUC on an independent test set [23].

1.4. Head and Neck QUS Treatment Response Prediction

Another application of QUS radiomics was for building classification models to predict treatment response in 59 patients (22 early responders (ER), 29 late-responders (LR), and 8 progressive diseases (PD)) with head and neck squamous cell carcinoma (HNSCC) undergoing radical radiotherapy (RT) with or without concurrent chemotherapy. The classification tasks identified non-responder (NR) from responders (R). Furthermore, they also separated NR into LR and those with PD [19]. Table 3 tabulates the classification metrics in the classification of ER from partial or NR. Table 4 summarizes the classification performance in the classification of LR versus persistent or PD. Figure 7 depicts a representative decision boundary, along with points from the two groups in a three-dimensional space of the selected best three features [19].
Tran et al. [57] reported classification models for predicting treatment response to radical RT in patients with head and neck malignancies. Their cohort consisted of 36 patients, with 14 complete responders (CR) and 22 partial responders (PR). Classification models were built using the best combination of attributes from a pool of radiomics features. The radiomics features consisted of first-order statistical mean-value and GLCM-based texture features of QUS spectral parametric images obtained at 24 h, 1 week, and 4 weeks post-treatment. They reported the best classification accuracies of 80, 86, and 85% at 24 h, week 1, and week 4, respectively [57].

2. Recurrence Prediction

2.1. Recurrence Prediction of Head and Neck Cancer Using Radiomics of QUS Spectral Parametric Imaging

Radiomics of QUS spectral parametric imaging have been applied for a clinical application in the identification of recurrence in 51 patients (17 recurrences vs. 34 non-recurrence) with head-neck squamous cell carcinoma (HNSCC) treated with RT with or without concurrent chemotherapy [24,59]. Dasgupta et al. implemented traditional machine learning in building a classification model for recurrence prediction. A set of pre-treatment radiomics features that included first-order statistics mean-value and GLCM-based textures were extracted from the parametric images. They reported the highest classification performance of 76% sensitivity, 71% specificity, 75% accuracy, and 0.74 AUC with the KNN (k-nearest neighbors) classifier [24].
Fatima et al. presented a radiomics study that develop a classification model for predicting recurrence in a cohort of 51 patients (17 recurrences vs. 34 non-recurrence) with HNSCC. In that study, radiomics imaging features from QUS spectral parametric images included mean-value first-order statistics and GLCM-based textures. Radiomics features from pre-treatment were subtracted from those collected at week 1 and week 4 after starting RT. They developed classification models using delta-features and reported classification performance of 80% and 82% accuracy in predicting recurrence using changes in the radiomics features at week 1 and week 4 of therapy, respectively [59].

2.2. Recurrence Prediction of Locally-Advanced Breast Cancer Using Radiomics of QUS Spectral Parametric Imaging

Dasgupta et al. reported a classification model to predict recurrence using radiomics features of QUS spectral parametric images obtained at pre-treatment [25]. Their cohort consisted of 83 patients (28 recurrences versus 55 non-recurrence) with LABC undergoing NAC. The pool of radiomics features consisted of first-order statistical mean-value, GLCM-based texture, and GLCM-based higher-order textures. They reported the best classification performance of 71% sensitivity, 87% specificity, 82% accuracy, and 0.76 AUC using an support vector machine (SVM) classifier. Using the same data, they also evaluated different targets in predicting recurrence. They reported a five-year recurrence-free survival of 83% and 54% (p = 0.003) and a five-year overall survival of 85% and 74% (p = 0.083) using the SVM model [25].
Bhardwaj et al. developed a classification model for predicting disease relapses in a cohort of 83 patients with LABC (28 recurrences and 55 non-recurrence) treated with NAC [60]. The pool of radiomics features included first-order statistical mean-value, GLCM-based textures, and GLCM-based higher-order textures from QUS spectral parametric images. They reported the best model from baseline and changes in week 4 features, achieving 87% sensitivity, 75% specificity, 81% accuracy, and 0.83 AUC [60].

3. Conclusions

In summary, QUS and its spectroscopic variables can be used to monitor treatment response, as well as to characterize responders versus non-responders and benign versus malignant tumor lesions and so on. Along with this, QUS is one of the fastest standard methods to predict treatment outcome. While most of the traditional methods incorporate tumor size measurement as a gold standard for response monitoring, that takes up to several days to weeks or even months. In contrast, imaging techniques such as QUS can provide details of tissue morphological alterations within a few hours. Thus, understanding the tumor heterogeneity and treatment response using QUS in advance might allow treatment switch if necessary at the earliest possible time.
Application of radiomics to QUS spectral parametric imaging has found a plethora of utilization in tumors characterization, treatment response characterization, monitoring, and prediction of disease relapses. Despite these promising avenues, caution needs to be exercised in regards to its implementation. These studies rely on the concept of learning from labeled data and the evaluation of generalized performance on new unseen data. Ongoing efforts are being made to refine the model building and evaluation strategies in order to eliminate potential pitfalls related to leakage of information from training to validation or test sets.
Radiomics features of QUS spectral parametric images are potential attributes for developing robust classification models for computer-aided diagnosis (CAD) systems with immense potential for its future incorporation clinically to evaluate cancer treatment responses.

Limitations and Future Directions

There are potential limitations associated with the usage of QUS. Since ultrasound is a handheld device, its reproducibility remains a question as compared to other techniques, such as mammography or MRI, that are automatic with high accuracy. Another important point is that most of the studies conducted using QUS are retrospectively performed at a single center, which suggests that the model might not be fully stable and there might also be a lack of robustness regarding the system. Although normalization and attenuation correction procedures compensate for instrument-specific components in the RF signal, there is a need to perform standardized measurements across institutions, clinical trials, and ultrasound system vendors. Such studies will corroborate the reproducibility and conformity of QUS backscatter measurements across possible sources of variations. Anderson et al. reported an interlaboratory study that independently measured ultrasound BSC from several well-characterized phantom materials [43]. In that study, the phantom materials were constructed with known scatterers’ distribution, allowing a direct comparison of the estimated BSCs with those obtained from acoustic scattering theory. The phantoms were grouped into low-attenuating and tissue-like attenuating materials. Ultrasonic backscattering characterizations of these phantoms were performed over frequency analysis ranging from 1 to 12 MHz. Each institution utilized several ultrasonic transducers with different characteristics, including nominal center frequencies, analysis bandwidths, aperture numbers, and focal lengths. The authors observed an excellent agreement between the BSC estimates from both laboratories. Their study demonstrates the ability to accurately estimate parameters derived from the BSC, including linear fit parameters, along with the scattering parameters. Thus, more inter-institutional studies should be conducted, which not only will provide large sample sizes, but also provide a variety of datasets that can improve the accuracy of the models.

References

Author Contributions

D.S., L.O.O. and G.J.C.; writing—original draft preparation, reviewing and editing. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Fass, L. Imaging and cancer: A review. Mol. Oncol. 2008, 2, 115–152. [Google Scholar] [CrossRef] [PubMed]
  2. Frangioni, J.V. New Technologies for Human Cancer Imaging. J. Clin. Oncol. 2008, 26, 4012–4021. [Google Scholar] [CrossRef] [PubMed]
  3. Oelze, M.L.; Mamou, J. Review of Quantitative Ultrasound: Envelope Statistics and Backscatter Coefficient Imaging and Contributions to Diagnostic Ultrasound. IEEE Trans. Ultrason. Ferroelectr. Freq. Control. 2016, 63, 336–351. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Sadeghi-Naini, A.; Papanicolau, N.; Falou, O.; Tadayyon, H.; Lee, J.; Zubovits, J.; Sadeghian, A.; Karshafian, R.; Al-Mahrouki, A.; Giles, A.; et al. Low-frequency quantitative ultrasound imaging of cell death in vivo. Med. Phys. 2013, 40, 082901. [Google Scholar] [CrossRef] [Green Version]
  5. Sadeghi-Naini, A.; Papanicolau, N.; Falou, O.; Zubovits, J.; Dent, R.; Verma, S.; Trudeau, M.; Boileau, J.F.; Spayne, J.; Iradji, S.; et al. Quantitative Ultrasound Evaluation of Tumor Cell Death Response in Locally Advanced Breast Cancer Patients Receiving Chemotherapy. Clin. Cancer Res. 2013, 19, 2163–2174. [Google Scholar] [CrossRef] [Green Version]
  6. Kim, H.C.; Al-Mahrouki, A.; Gorjizadeh, A.; Sadeghi-Naini, A.; Karshafian, R.; Czarnota, G.J. Quantitative Ultrasound Characterization of Tumor Cell Death: Ultrasound-Stimulated Microbubbles for Radiation Enhancement. PLoS ONE 2014, 9, e102343. [Google Scholar] [CrossRef]
  7. Sadeghi-Naini, A.; Zhou, S.; Gangeh, M.J.; Jahedmotlagh, Z.; Falou, O.; Ranieri, S.; Azrif, M.; Giles, A.; Czarnota, G.J. Quantitative evaluation of cell death response in vitro and in vivo using conventional-frequency ultrasound. Oncoscience 2015, 2, 716–726. [Google Scholar] [CrossRef] [Green Version]
  8. Tran, W.T.; Sannachi, L.; Papanicolau, N.; Tadayyon, H.; Al Mahrouki, A.; El Kaffas, A.; Gorjizadeh, A.; Lee, J.; Czarnota, G.J. Quantitative ultrasound imaging of therapy response in bladder cancer in vivo. Oncoscience 2016, 3, 122–133. [Google Scholar] [CrossRef] [Green Version]
  9. Sharma, D.; Osapoetra, L.O.; Faltyn, M.; Do, N.N.A.; Giles, A.; Stanisz, M.; Sannachi, L.; Czarnota, G.J. Quantitative ultrasound characterization of therapy response in prostate cancer in vivo. Am. J. Transl. Res. 2021, 13, 4437–4449. [Google Scholar]
  10. Sadeghi-Naini, A.; Falou, O.; Tadayyon, H.; Al-Mahrouki, A.; Tran, W.; Papanicolau, N.; Kolios, M.C.; Czarnota, G.J. Conventional Frequency Ultrasonic Biomarkers of Cancer Treatment Response In Vivo. Transl. Oncol. 2013, 6, 234–243. [Google Scholar] [CrossRef] [Green Version]
  11. Lizzi, F.L.; Ostromogilsky, M.; Feleppa, E.J.; Rorke, M.C.; Yaremko, M.M. Relationship of Ultrasonic Spectral Parameters to Features of Tissue Microstructure. IEEE Trans. Ultrason. Ferroelectr. Freq. Control. 1987, 34, 319–329. [Google Scholar] [CrossRef] [PubMed]
  12. Lizzi, F.L.; Astor, M.; Feleppa, E.J.; Shao, M.; Kalisz, A. Statistical framework for ultrasonic spectral parameter imaging. Ultrasound Med. Biol. 1997, 23, 1371–1382. [Google Scholar] [CrossRef] [PubMed]
  13. Sannachi, L.; Tadayyon, H.; Sadeghi-Naini, A.; Tran, W.; Gandhi, S.; Wright, F.; Oelze, M.; Czarnota, G. Non-invasive evaluation of breast cancer response to chemotherapy using quantitative ultrasonic backscatter parameters. Med. Image Anal. 2015, 20, 224–236. [Google Scholar] [CrossRef] [PubMed]
  14. Tadayyon, H.; Sannachi, L.; Sadeghi-Naini, A.; Al-Mahrouki, A.; Tran, W.T.; Kolios, M.C.; Czarnota, G.J. Quantification of Ultrasonic Scattering Properties of In Vivo Tumor Cell Death in Mouse Models of Breast Cancer. Transl. Oncol. 2015, 8, 463–473. [Google Scholar] [CrossRef] [Green Version]
  15. Ghoshal, G.; Luchies, A.C.; Blue, J.P.; Oelze, M.L. Temperature dependent ultrasonic characterization of biological media. J. Acoust. Soc. Am. 2011, 130, 2203–2211. [Google Scholar] [CrossRef] [Green Version]
  16. Frank, G.A.; Danilova, N.V.; Andreeva, Y.Y.; Nefedova, N.A. WHO classification of tumors of the breast, 2012. Arkhiv Patol. 2013, 75, 53–63. [Google Scholar]
  17. Ogston, K.N.; Miller, I.D.; Payne, S.; Hutcheon, A.W.; Sarkar, T.K.; Smith, I.; Schofield, A.; Heys, S.D. A new histological grading system to assess response of breast cancers to primary chemotherapy: Prognostic significance and survival. Breast 2003, 12, 320–327. [Google Scholar] [CrossRef]
  18. Osapoetra, L.O.; Sannachi, L.; Quiaoit, K.; Dasgupta, A.; DiCenzo, D.; Fatima, K.; Wright, F.; Dinniwell, R.; Trudeau, M.; Gandhi, S.; et al. A priori prediction of response in multicentre locally advanced breast cancer (LABC) patients using quantitative ultrasound and derivative texture methods. Oncotarget 2021, 12, 81–94. [Google Scholar] [CrossRef]
  19. Osapoetra, L.O.; Dasgupta, A.; DiCenzo, D.; Fatima, K.; Quiaoit, K.; Saifuddin, M.; Karam, I.; Poon, I.; Husain, Z.; Tran, W.T.; et al. Assessment of clinical radiosensitivity in patients with head-neck squamous cell carcinoma from pre-treatment quantitative ultrasound radiomics. Sci. Rep. 2021, 11, 6117. [Google Scholar] [CrossRef]
  20. Sannachi, L.; Gangeh, M.; Naini, A.-S.; Bhargava, P.; Jain, A.; Tran, W.T.; Czarnota, G.J. Quantitative Ultrasound Monitoring of Breast Tumour Response to Neoadjuvant Chemotherapy: Comparison of Results Among Clinical Scanners. Ultrasound Med. Biol. 2020, 46, 1142–1157. [Google Scholar] [CrossRef]
  21. DiCenzo, D.; Quiaoit, K.; Fatima, K.; Bhardwaj, D.; Sannachi, L.; Gangeh, M.; Sadeghi-Naini, A.; Dasgupta, A.; Kolios, M.C.; Trudeau, M.; et al. Quantitative ultrasound radiomics in predicting response to neoadjuvant chemotherapy in patients with locally advanced breast cancer: Results from multi-institutional study. Cancer Med. 2020, 9, 5798–5806. [Google Scholar] [CrossRef] [PubMed]
  22. Taleghamar, H.; Moghadas-Dastjerdi, H.; Czarnota, G.J.; Sadeghi-Naini, A. Characterizing intra-tumor regions on quantitative ultrasound parametric images to predict breast cancer response to chemotherapy at pre-treatment. Sci. Rep. 2021, 11, 14865. [Google Scholar] [CrossRef] [PubMed]
  23. Taleghamar, H.; Jalalifar, S.A.; Czarnota, G.J.; Sadeghi-Naini, A. Deep learning of quantitative ultrasound multi-parametric images at pre-treatment to predict breast cancer response to chemotherapy. Sci. Rep. 2022, 12, 2244. [Google Scholar] [CrossRef] [PubMed]
  24. Dasgupta, A.; Fatima, K.; DiCenzo, D.; Bhardwaj, D.; Quiaoit, K.; Saifuddin, M.; Karam, I.; Poon, I.; Husain, Z.; Tran, W.T.; et al. Quantitative ultrasound radiomics in predicting recurrence for patients with node-positive head-neck squamous cell carcinoma treated with radical radiotherapy. Cancer Med. 2021, 10, 2579–2589. [Google Scholar] [CrossRef] [PubMed]
  25. Dasgupta, A.; Bhardwaj, D.; DiCenzo, D.; Fatima, K.; Osapoetra, L.O.; Quiaoit, K.; Saifuddin, M.; Brade, S.; Trudeau, M.; Gandhi, S.; et al. Radiomics in predicting recurrence for patients with locally advanced breast cancer using quantitative ultrasound. Oncotarget 2021, 12, 2437–2448. [Google Scholar] [CrossRef]
  26. Yao, L.X.; Zagzebski, J.A.; Madsen, E.L. Backscatter Coefficient Measurements Using a Reference Phantom to Extract Depth-Dependent Instrumentation Factors. Ultrason. Imaging 1990, 12, 58–70. [Google Scholar] [CrossRef]
  27. Labyed, Y.; Bigelow, T.A.; McFarlin, B.L. Estimate of the attenuation coefficient using a clinical array transducer for the detection of cervical ripening in human pregnancy. Ultrasonics 2011, 51, 34–39. [Google Scholar] [CrossRef] [Green Version]
  28. Banihashemi, B.; Vlad, R.; Debeljevic, B.; Giles, A.; Kolios, M.C.; Czarnota, G.J.; Sharkey, R.M.; Karacay, H.; Litwin, S.; Rossi, E.A.; et al. Ultrasound Imaging of Apoptosis in Tumor Response: Novel Preclinical Monitoring of Photodynamic Therapy Effects. Cancer Res. 2008, 68, 8590–8596. [Google Scholar] [CrossRef] [Green Version]
  29. Czarnota, G.J.; Kolios, M.C.; Vaziri, H.; Benchimol, S.; Ottensmeyer, F.P.; Sherar, M.D.; Hunt, J.W. Ultrasonic biomicroscopy of viable, dead and apoptotic cells. Ultrasound Med. Biol. 1997, 23, 961–965. [Google Scholar] [CrossRef]
  30. Czarnota, G.J.; Kolios, M.C.; Abraham, J.; Portnoy, M.; Ottensmeyer, F.P.; Hunt, J.W.; Sherar, M.D. Ultrasound imaging of apoptosis: High-resolution non-invasive monitoring of programmed cell death in vitro, in situ and in vivo. Br. J. Cancer 1999, 81, 520–527. [Google Scholar] [CrossRef]
  31. Feleppa, E.J.; Ennis, R.D.; Schiff, P.B.; Wuu, C.-S.; Kalisz, A.; Ketterling, J.; Urban, S.; Liu, T.; Fair, W.R.; Porter, C.R.; et al. Spectrum-analysis and neural networks for imaging to detect and treat prostate cancer. Ultrason. Imaging 2001, 23, 135–146. [Google Scholar] [CrossRef] [PubMed]
  32. Feleppa, E.J.; Ennis, R.D.; Schiff, P.B.; Wuu, C.-S.; Kalisz, A.; Ketterling, J.; Urban, S.; Liu, T.; Fair, W.R.; Porter, C.R.; et al. Ultrasonic spectrum-analysis and neural-network classification as a basis for ultrasonic imaging to target brachytherapy of prostate cancer. Brachytherapy 2002, 1, 48–53. [Google Scholar] [CrossRef] [PubMed]
  33. Kolios, M.C.; Czarnota, G.J.; Lee, M.; Hunt, J.W.; Sherar, M.D. Ultrasonic spectral parameter characterization of apoptosis. Ultrasound Med. Biol. 2002, 28, 589–597. [Google Scholar] [CrossRef] [PubMed]
  34. Oelze, M.L.; Zachary, J.F.; O’Brien, W.D. Characterization of tissue microstructure using ultrasonic backscatter: Theory and technique for optimization using a Gaussian form factor. J. Acoust. Soc. Am. 2002, 112, 1202–1211. [Google Scholar] [CrossRef] [Green Version]
  35. Oelze, M.L.; Zachary, J.F.; O’Brien, W.D. Parametric imaging of rat mammary tumors in vivo for the purposes of tissue characterization. J. Ultrasound Med. 2002, 21, 1201–1210. [Google Scholar] [CrossRef]
  36. Oelze, M.L.; O’Brien, W.D.; Blue, J.P.; Zachary, J.F. Differentiation and Characterization of Rat Mammary Fibroadenomas and 4T1 Mouse Carcinomas Using Quantitative Ultrasound Imaging. IEEE Trans. Med. Imaging 2004, 23, 764–771. [Google Scholar] [CrossRef]
  37. Feleppa, E.J.; Porter, C.R.; Ketterling, J.; Lee, P.; Dasgupta, S.; Urban, S.; Kalisz, A. Recent Developments in Tissue-Type Imaging (TTI) for Planning and Monitoring Treatment of Prostate Cancer. Ultrason. Imaging 2004, 26, 163–172. [Google Scholar] [CrossRef]
  38. Vlad, R.M.; Czarnota, G.J.; Giles, A.; Sherar, M.D.; Hunt, J.W.; Kolios, M.C. High-frequency ultrasound for monitoring changes in liver tissue during preservation. Phys. Med. Biol. 2005, 50, 197–213. [Google Scholar] [CrossRef] [Green Version]
  39. Tunis, A.S.; Czarnota, G.J.; Giles, A.; Sherar, M.D.; Hunt, J.W.; Kolios, M.C. Monitoring structural changes in cells with high-frequency ultrasound signal statistics. Ultrasound Med. Biol. 2005, 31, 1041–1049. [Google Scholar] [CrossRef]
  40. Oelze, M.L.; Zachary, J.F. Examination of cancer in mouse models using high-frequency quantitative ultrasound. Ultrasound Med. Biol. 2006, 32, 1639–1648. [Google Scholar] [CrossRef]
  41. Vlad, R.M.; Alajez, N.M.; Giles, A.; Kolios, M.C.; Czarnota, G.J. Quantitative Ultrasound Characterization of Cancer Radiotherapy Effects In Vitro. IInt. J. Radiat. Oncol. Biol. Phys. 2008, 72, 1236–1243. [Google Scholar] [CrossRef] [PubMed]
  42. Vlad, R.M.; Brand, S.; Giles, A.; Kolios, M.C.; Czarnota, G.J. Quantitative Ultrasound Characterization of Responses to Radiotherapy in Cancer Mouse Models. Clin. Cancer Res. 2009, 15, 2067–2075. [Google Scholar] [CrossRef] [PubMed]
  43. Anderson, J.J.; Herd, M.-T.; King, M.R.; Haak, A.; Hafez, Z.T.; Song, J.; Oelze, M.L.; Madsen, E.L.; Zagzebski, J.A.; O’Brien, J.W.D.; et al. Interlaboratory Comparison of Backscatter Coefficient Estimates for Tissue-Mimicking Phantoms. Ultrason. Imaging 2010, 32, 48–64. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  44. Mamou, J.; Coron, A.; Hata, M.; Machi, J.; Yanagihara, E.; Laugier, P.; Feleppa, E.J. Three-Dimensional High-Frequency Characterization of Cancerous Lymph Nodes. Ultrasound Med. Biol. 2010, 36, 361–375. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  45. Mamou, J.; Coron, A.; Oelze, M.L.; Saegusa-Beecroft, E.; Hata, M.; Lee, P.; Machi, J.; Yanagihara, E.; Laugier, P.; Feleppa, E.J. Three-Dimensional High-Frequency Backscatter and Envelope Quantification of Cancerous Human Lymph Nodes. Ultrasound Med. Biol. 2011, 37, 345–357. [Google Scholar] [CrossRef] [Green Version]
  46. Lee, J.; Karshafian, R.; Papanicolau, N.; Giles, A.; Kolios, M.C.; Czarnota, G.J. Quantitative Ultrasound for the Monitoring of Novel Microbubble and Ultrasound Radiosensitization. Ultrasound Med. Biol. 2012, 38, 1212–1221. [Google Scholar] [CrossRef]
  47. Nam, K.; Zagzebski, J.A.; Hall, T.J. Quantitative Assessment of In Vivo Breast Masses Using Ultrasound Attenuation and Backscatter. Ultrason. Imaging 2013, 35, 146–161. [Google Scholar] [CrossRef] [Green Version]
  48. Saegusa-Beecroft, E.; Machi, J.; Mamou, J.; Hata, M.; Coron, A.; Yanagihara, E.T.; Yamaguchi, T.; Oelze, M.L.; Laugier, P.; Feleppa, E.J. Three-dimensional quantitative ultrasound for detecting lymph node metastases. J. Surg. Res. 2013, 183, 258–269. [Google Scholar] [CrossRef] [Green Version]
  49. Lavarello, R.J.; Ridgway, W.R.; Sarwate, S.S.; Oelze, M.L. Characterization of Thyroid Cancer in Mouse Models Using High-Frequency Quantitative Ultrasound Techniques. Ultrasound Med. Biol. 2013, 39, 2333–2341. [Google Scholar] [CrossRef] [Green Version]
  50. Tadayyon, H.; Sadeghi-Naini, A.; Wirtzfeld, L.; Wright, F.C.; Czarnota, G. Quantitative ultrasound characterization of locally advanced breast cancer by estimation of its scatterer properties. Med. Phys. 2014, 41, 012903. [Google Scholar] [CrossRef]
  51. Ghoshal, G.; Kemmerer, J.P.; Karunakaran, C.; Abuhabsah, R.; Miller, R.J.; Sarwate, S.; Oelze, M.L. Quantitative ultrasound imaging for monitoring in situ high-intensity focused ultrasound exposure. Ultrason. Imaging 2014, 36, 239–255. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  52. Ghoshal, G.; Kemmerer, J.P.; Karunakaran, C.; Miller, R.J.; Oelze, M.L. Quantitative Ultrasound for Monitoring High-Intensity Focused Ultrasound Treatment In Vivo. IEEE Trans. Ultrason. Ferroelectr. Freq. Control. 2016, 63, 1234–1242. [Google Scholar] [CrossRef] [PubMed]
  53. Pasternak, M.M.; Wirtzfeld, L.A.; Kolios, M.C.; Czarnota, G.J. High-Frequency Ultrasound Analysis of Post-Mitotic Arrest Cell Death. Oncoscience 2016, 3, 109–121. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  54. Sadeghi-Naini, A.; Suraweera, H.; Tran, W.T.; Hadizad, F.; Bruni, G.; Rastegar, R.F.; Curpen, B.; Czarnota, G.J. Breast-Lesion Characterization using Textural Features of Quantitative Ultrasound Parametric Maps. Sci. Rep. 2017, 7, 13638. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  55. Rohrbach, D.; Wodlinger, B.; Wen, J.; Mamou, J.; Feleppa, E. High-Frequency Quantitative Ultrasound for Imaging Prostate Cancer Using a Novel Micro-Ultrasound Scanner. Ultrasound Med. Biol. 2018, 44, 1341–1354. [Google Scholar] [CrossRef]
  56. Quiaoit, K.; DiCenzo, D.; Fatima, K.; Bhardwaj, D.; Sannachi, L.; Gangeh, M.; Sadeghi-Naini, A.; Dasgupta, A.; Kolios, M.C.; Trudeau, M.; et al. Quantitative ultrasound radiomics for therapy response monitoring in patients with locally advanced breast cancer: Multi-institutional study results. PLoS ONE 2020, 15, e0236182. [Google Scholar] [CrossRef]
  57. Tran, W.T.; Suraweera, H.; Quiaoit, K.; DiCenzo, D.; Fatima, K.; Jang, D.; Bhardwaj, D.; Kolios, C.; Karam, I.; Poon, I.; et al. Quantitative ultrasound delta-radiomics during radiotherapy for monitoring treatment responses in head and neck malignancies. Futur. Sci. OA 2020, 6, 1–15. [Google Scholar] [CrossRef]
  58. Dasgupta, A.; Brade, S.; Snnachi, L.; Quiaoit, K.; Fatima, K.; DiCenzo, D.; Osapoetra, L.O.; Saifuddin, M.; Trudeau, M.; Gandhi, S.; et al. Quantitative ultrasound radiomics using texture derivatives in prediction of treatment response to neo-adjuvant chemotherapy for locally advanced breast cancer. Oncotarget 2020, 11, 3782–3792. [Google Scholar] [CrossRef]
  59. Fatima, K.; Dasgupta, A.; DiCenzo, D.; Kolios, C.; Quiaoit, K.; Saifuddin, M.; Sandhu, M.; Bhardwaj, D.; Karam, I.; Poon, I.; et al. Ultrasound delta-radiomics during radiotherapy to predict recurrence in patients with head and neck squamous cell carcinoma. Clin. Transl. Radiat. Oncol. 2021, 28, 62–70. [Google Scholar] [CrossRef]
  60. Bhardwaj, D.; Dasgupta, A.; DiCenzo, D.; Brade, S.; Fatima, K.; Quiaoit, K.; Trudeau, M.; Gandhi, S.; Eisen, A.; Wright, F.; et al. Early Changes in Quantitative Ultrasound Imaging Parameters during Neoadjuvant Chemotherapy to Predict Recurrence in Patients with Locally Advanced Breast Cancer. Cancers 2022, 14, 1247. [Google Scholar] [CrossRef]
  61. Osapoetra, L.O.; Sannachi, L.; DiCenzo, D.; Quiaoit, K.; Fatima, K.; Czarnota, G.J. Breast lesion characterization using Quantitative Ultrasound (QUS) and derivative texture methods. Transl. Oncol. 2020, 13, 100827. [Google Scholar] [CrossRef] [PubMed]
  62. Osapoetra, L.O.; Chan, W.; Tran, W.; Kolios, M.C.; Czarnota, G.J. Comparison of methods for texture analysis of QUS parametric images in the characterization of breast lesions. PLoS ONE 2020, 15, e0244965. [Google Scholar] [CrossRef] [PubMed]
  63. Goundan, P.N.; Mamou, J.; Rohrbach, D.; Smith, J.; Patel, H.; Wallace, K.D.; Feleppa, E.J.; Lee, S.L. A Preliminary Study of Quantitative Ultrasound for Cancer-Risk Assessment of Thyroid Nodules. Front. Endocrinol. 2021, 12, 627698. [Google Scholar] [CrossRef] [PubMed]
  64. Najafabadi, M.M.; Villanustre, F.; Khoshgoftaar, T.M.; Seliya, N.; Wald, R.; Muharemagic, E. Deep learning applications and challenges in big data analytics. J. Big Data 2015, 2, 1. [Google Scholar] [CrossRef] [Green Version]
  65. He, K.; Zhang, X.; Ren, S.; Sun, J. Deep Residual Learning for Image Recognition. In Proceedings of the 2016 IEEE Conference on Computer Vision and Pattern Recognition (CVPR), Las Vegas, NV, USA, 27–30 June 2016; pp. 770–778. [Google Scholar]
  66. Wang, F.; Jiang, M.; Qian, C.; Yang, S.; Li, C.; Zhang, H.; Wang, X.; Tang, X. Residual Attention Network for Image Classification. In Proceedings of the 2017 IEEE Conference on Computer Vision and Pattern Recognition (CVPR), Honolulu, HI, USA, 21–26 July 2017; pp. 6450–6458. [Google Scholar] [CrossRef]
Figure 1. A summary of QUS spectroscopy workflow. QUS spectroscopy extracts scattering components from the RF signal that uniquely characterize tissue samples. The spectral normalization procedure removes instrument-dependent components from the RF signal, while the attenuation correction procedure compensates loss in the acoustic signal as it propagates through intervening tissue layers. Parametrization of the attenuation-corrected NPS or the BSC allows for tumors characterization, providing diagnostic and prognostic values. These parameters can then be utilized to develop a classification model for diagnostic and prognostic purposes. BSC: backscatter coefficient, NPS: normalized power spectrum, QUS: quantitative ultrasound.
Figure 1. A summary of QUS spectroscopy workflow. QUS spectroscopy extracts scattering components from the RF signal that uniquely characterize tissue samples. The spectral normalization procedure removes instrument-dependent components from the RF signal, while the attenuation correction procedure compensates loss in the acoustic signal as it propagates through intervening tissue layers. Parametrization of the attenuation-corrected NPS or the BSC allows for tumors characterization, providing diagnostic and prognostic values. These parameters can then be utilized to develop a classification model for diagnostic and prognostic purposes. BSC: backscatter coefficient, NPS: normalized power spectrum, QUS: quantitative ultrasound.
Cancers 14 06217 g001
Figure 2. Representative QUS spectral parametric images from benign and malignant group. (A) benign breast lesions (left columns) (B) malignant breast lesions (right columns). The color bars present the range for ASD of 160 μm, AAC of 70 dB/cm3, MBF of 44 dB, SS of 10 dB/MHz, and SI of 70 dB. Scale bar: 1 cm. ASD: average scatterer diameter, AAC: average acoustic concentration, MBF: mid-band fit, SS: spectral slope, SI: 0-MHz intercept. Reprinted with permission from: (Figure 1) [61].
Figure 2. Representative QUS spectral parametric images from benign and malignant group. (A) benign breast lesions (left columns) (B) malignant breast lesions (right columns). The color bars present the range for ASD of 160 μm, AAC of 70 dB/cm3, MBF of 44 dB, SS of 10 dB/MHz, and SI of 70 dB. Scale bar: 1 cm. ASD: average scatterer diameter, AAC: average acoustic concentration, MBF: mid-band fit, SS: spectral slope, SI: 0-MHz intercept. Reprinted with permission from: (Figure 1) [61].
Cancers 14 06217 g002
Figure 3. Representative texture maps representing local quantification of image texture. (A) benign breast lesions (left columns) (B) malignant breast lesions (right columns). Scale bar: 1 cm. MBF: mid-band fit. Reprinted with permission from: (Figure 2) [61].
Figure 3. Representative texture maps representing local quantification of image texture. (A) benign breast lesions (left columns) (B) malignant breast lesions (right columns). Scale bar: 1 cm. MBF: mid-band fit. Reprinted with permission from: (Figure 2) [61].
Cancers 14 06217 g003
Figure 4. Box and scatter plots of representative radiomics features that demonstrate statistically significant differences. Statistical significant, are represented as (* p < 0.05), (** p < 0.01), and (*** p < 0.001). B: benign breast lesions, M: malignant breast lesions, ASD: average scatterer diameter, AAC: average acoustic concentration, MBF: mid-band fit, SS: spectral slope, SI: 0-MHz intercept, CON: contrast, COR: correlation, ENE: energy, HOM: homogeneity. Reprinted with permission from: (Figure 3) [61].
Figure 4. Box and scatter plots of representative radiomics features that demonstrate statistically significant differences. Statistical significant, are represented as (* p < 0.05), (** p < 0.01), and (*** p < 0.001). B: benign breast lesions, M: malignant breast lesions, ASD: average scatterer diameter, AAC: average acoustic concentration, MBF: mid-band fit, SS: spectral slope, SI: 0-MHz intercept, CON: contrast, COR: correlation, ENE: energy, HOM: homogeneity. Reprinted with permission from: (Figure 3) [61].
Cancers 14 06217 g004
Figure 5. The root means square deviation of differences in QUS parameters: MBF and SS, MBF-texture parameters: MBF-CON, MBF-COR, MBF-ENE, and MBF-HOM due to measurement uncertainty, variations in ultrasound systems, and tissue heterogeneity. These results suggest that the inherent tissue heterogeneity was the most dominant contributor to variations in the estimated radiomics features of QUS spectral parametric images. Overall, measurement uncertainty and variations in clinical ultrasound systems contribute less than the tissue heterogeneity component. MBF: mid-band fit, SS: spectral slope, CON: contrast, COR: correlation, ENE: energy, HOM: homogeneity. Reprinted with permission from: (Figure 3) [20].
Figure 5. The root means square deviation of differences in QUS parameters: MBF and SS, MBF-texture parameters: MBF-CON, MBF-COR, MBF-ENE, and MBF-HOM due to measurement uncertainty, variations in ultrasound systems, and tissue heterogeneity. These results suggest that the inherent tissue heterogeneity was the most dominant contributor to variations in the estimated radiomics features of QUS spectral parametric images. Overall, measurement uncertainty and variations in clinical ultrasound systems contribute less than the tissue heterogeneity component. MBF: mid-band fit, SS: spectral slope, CON: contrast, COR: correlation, ENE: energy, HOM: homogeneity. Reprinted with permission from: (Figure 3) [20].
Cancers 14 06217 g005
Figure 6. Bland-Altman plots comparing radiomics features obtained using an Ultrasonix-RP (Ultrasonix Medical Corp., Richmond, BC, Canada) and a GEGE-LOGIQ E9 (General Electric Healthcare, Milwaukee, WI, USA) ultrasound system. The limits of agreement were indicated by the mean difference and ± standard deviation. MBF: mid-band fit, SS: spectral slope, CON: contrast, COR: correlation, ENE: energy, HOM: homogeneity, ULX: ultrasonix RP L14-5/60, GE: GE-LOGIQ 9 L-D. Reprinted with permission from: (Figure 6) [20].
Figure 6. Bland-Altman plots comparing radiomics features obtained using an Ultrasonix-RP (Ultrasonix Medical Corp., Richmond, BC, Canada) and a GEGE-LOGIQ E9 (General Electric Healthcare, Milwaukee, WI, USA) ultrasound system. The limits of agreement were indicated by the mean difference and ± standard deviation. MBF: mid-band fit, SS: spectral slope, CON: contrast, COR: correlation, ENE: energy, HOM: homogeneity, ULX: ultrasonix RP L14-5/60, GE: GE-LOGIQ 9 L-D. Reprinted with permission from: (Figure 6) [20].
Cancers 14 06217 g006
Figure 7. A representative nonlinear decision boundary obtained from an SVM model using the RBF kernel to separate complete responders from those with partial/non-responders. The model consisted of the best combination of three radiomics features obtained through a forward sequential feature selection with cross-validation. SVM: support vector machines, RBF: radial basis functions, MBF: mid-band fit, CON: contrast, ENE: energy, HOM: homogeneity. Reprinted with permission from: (Figure 3) [19].
Figure 7. A representative nonlinear decision boundary obtained from an SVM model using the RBF kernel to separate complete responders from those with partial/non-responders. The model consisted of the best combination of three radiomics features obtained through a forward sequential feature selection with cross-validation. SVM: support vector machines, RBF: radial basis functions, MBF: mid-band fit, CON: contrast, ENE: energy, HOM: homogeneity. Reprinted with permission from: (Figure 3) [19].
Cancers 14 06217 g007
Table 1. Summary of QUS application and utilization in preclinical and clinical studies.
Table 1. Summary of QUS application and utilization in preclinical and clinical studies.
QUS FrequencyImplementationReferences
50 MHzCell death characterization (in vitro)[29]
40 MHzMonitoring treatment response (in vitro, in situ and in vivo)[30]
7 MHzTissue characterization (clinical)[31]
5.75 MHzTissue characterization (clinical)[32]
30 MHz and 34 MHzCell death characterization (in vitro)[33]
4–12 MHzTissue characterization (in vivo)[34]
8 MHzTissue characterization (in vivo)[35]
8.5 MHz and 20 MHzTissue characterization (in vivo)[36]
7.5 MHzTissue characterization (clinical)[37]
40 MHzTissue characterization (in vivo)[38]
20 MHzExamining cell structural changes (in vitro)[39]
20 MHzTissue characterization (in vivo)[40]
26 MHzMonitoring treatment response (in vivo)[28]
20 MHzMonitoring treatment response (in vitro)[41]
20 MHzMonitoring treatment response (in vivo)[42]
1–12 MHzTissue characterization (phantoms) (clinical)[43]
25.6 MHzLymph nodes characterization (clinical)[44]
25.6 MHzLymph nodes characterization (clinical)[45]
25 MHzMonitoring treatment response (in vivo)[46]
10 MHz and 15 MHzTissue characterization (phantoms) (clinical)[47]
7 MHzMonitoring treatment response (clinical)[5]
25.6 MHzDetecting lymph node metastases (clinical)[48]
∼7 MHz and 20 MHzTreatment response monitoring (in vivo)[4]
40 MHzTissue characterization (in vivo)[49]
6 MHzTissue characterization (clinical)[50]
25 MHzMonitoring treatment response (in vivo)[6]
6 MHzMonitoring treatment response (in situ)[51]
7 MHzMonitoring treatment response (clinical)[13]
~7 MHzMonitoring treatment response (in vitro, in vivo)[7]
7 MHz and 20 MHzMonitoring treatment response (in vivo)[14]
6 MHzMonitoring treatment response (in vivo)[52]
25 MHzCell death characterization (in vitro)[53]
25 MHzMonitoring treatment response (in vivo)[8]
~6 MHzTissue characterization (clinical)[54]
29 MHzTissue characterization (clinical)[55]
6.5 MHz and 6.9 MHzPredicting treatment response (clinical)[21]
6.3 MHz and 7 MHzMonitoring treatment response (clinical)[56]
8 MHzMonitoring treatment response (clinical)[57]
7 MHzPredicting treatment response (clinical)[58]
6.5 MHzPredicting tumor recurrence (clinical)[59]
6.5 MHzPredicting treatment response (clinical)[19]
7 MHzPredicting tumor recurrence (clinical)[25]
7 MHzPredicting tumor recurrence (clinical)[60]
Table 2. Summary of classification performance in the classification of non-responder from responder at weeks 1, 4, and 8 after the initiation of chemotherapy using data from two ultrasound systems. The asterisk (*) represents a statistically significant difference (p < 0.05) carried between treatment response predictions from two ultrasound systems. The classification model for treatment response monitoring was based on a previously reported study. US: ultrasound, ULX: ULX L14- 5/60, GE: GE-LOGIQ E9. Adapted with permission from [20].
Table 2. Summary of classification performance in the classification of non-responder from responder at weeks 1, 4, and 8 after the initiation of chemotherapy using data from two ultrasound systems. The asterisk (*) represents a statistically significant difference (p < 0.05) carried between treatment response predictions from two ultrasound systems. The classification model for treatment response monitoring was based on a previously reported study. US: ultrasound, ULX: ULX L14- 5/60, GE: GE-LOGIQ E9. Adapted with permission from [20].
Scan TimeUS SystemSensitivity [%]Specificity [%]Accuracy [%]McNemar p *
Week 1ULX60.050.058.80.752
GE60.050.058.8
Week 4ULX78.966.777.30.545
GE64.366.769.1
Week 8ULX71.4100.073.30.683
GE71.4100.073.3
Table 3. Response at 3 months: Complete responder (early responder) (n = 22) vs. partial/non-responder (n = 37). Representative performance of classification models for identifying complete responders from partial or non-responder. The analysis indicated the SVM-RBF model performed the best in this particular task. FLD: Fischer’s linear discriminant analysis, KNN: k-nearest neighbour, SVM-RBF: support vector machine-radial based function, MBF: mid-band fit, SI: 0-MHz intercept, SS: spectral slope, ASD: average scatterer diameter, AAC: average acoustic concentration, COR: correlation, CON: contrast, HOM: homogeneity, ENE: energy. Adapted with permission from [19].
Table 3. Response at 3 months: Complete responder (early responder) (n = 22) vs. partial/non-responder (n = 37). Representative performance of classification models for identifying complete responders from partial or non-responder. The analysis indicated the SVM-RBF model performed the best in this particular task. FLD: Fischer’s linear discriminant analysis, KNN: k-nearest neighbour, SVM-RBF: support vector machine-radial based function, MBF: mid-band fit, SI: 0-MHz intercept, SS: spectral slope, ASD: average scatterer diameter, AAC: average acoustic concentration, COR: correlation, CON: contrast, HOM: homogeneity, ENE: energy. Adapted with permission from [19].
ClassifierSensitivitySpecificityAccuracyAUCSelected Features
FLD7381780.75MBF-HOM-CON, MBF, SI-CON-ENE
KNN7384800.80SS-COR-COR, MBF-ENE-HOM
SVM-RBF8695920.91MBF-HOM-CON, MBF-ENE-CON, ASD-HOM-ENE
Table 4. Final Response. Late responder (n = 29) vs. persistent/progressive disease (n = 8). FLD: Fischer’s linear discriminant analysis, KNN: k-nearest neighbour, SVM-RBF: support vector machine-radial based function, MBF: mid-band fit, SI: 0-MHz intercept, SS: spectral slope, ASD: average scatterer diameter, AAC: average acoustic concentration, CON: contrast, COR: correlation, ENE: energy, HOM: homogeneity. Adapted with permission from [19].
Table 4. Final Response. Late responder (n = 29) vs. persistent/progressive disease (n = 8). FLD: Fischer’s linear discriminant analysis, KNN: k-nearest neighbour, SVM-RBF: support vector machine-radial based function, MBF: mid-band fit, SI: 0-MHz intercept, SS: spectral slope, ASD: average scatterer diameter, AAC: average acoustic concentration, CON: contrast, COR: correlation, ENE: energy, HOM: homogeneity. Adapted with permission from [19].
ClassifierSensitivitySpecificityAccuracyAUCSelected Features
FLD86100890.92AAC-ENE-HOM, AAC-HOM-CON
KNN9388920.90AAC-HOM, ASD-ENE-HOM
SVM-RBF9788950.97SS, SS-HOM-CON
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Sharma, D.; Osapoetra, L.O.; Czarnota, G.J. Implementation of Non-Invasive Quantitative Ultrasound in Clinical Cancer Imaging. Cancers 2022, 14, 6217. https://doi.org/10.3390/cancers14246217

AMA Style

Sharma D, Osapoetra LO, Czarnota GJ. Implementation of Non-Invasive Quantitative Ultrasound in Clinical Cancer Imaging. Cancers. 2022; 14(24):6217. https://doi.org/10.3390/cancers14246217

Chicago/Turabian Style

Sharma, Deepa, Laurentius Oscar Osapoetra, and Gregory J. Czarnota. 2022. "Implementation of Non-Invasive Quantitative Ultrasound in Clinical Cancer Imaging" Cancers 14, no. 24: 6217. https://doi.org/10.3390/cancers14246217

APA Style

Sharma, D., Osapoetra, L. O., & Czarnota, G. J. (2022). Implementation of Non-Invasive Quantitative Ultrasound in Clinical Cancer Imaging. Cancers, 14(24), 6217. https://doi.org/10.3390/cancers14246217

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop