Radiomic Analysis of Pharmacokinetic Heterogeneity Within Tumor Based on the Unsupervised Decomposition of Dynamic Contrast-Enhanced MRI for Predicting Histological Characteristics of Breast Cancer

Preliminary study on DCE-MRI radiomics analysis for differentiation of HER2-low and HER2-zero breast cancer

Purpose

This study aims to evaluate the utility of radiomic features from dynamic contrast-enhanced magnetic resonance imaging (DCE-MRI) in distinguishing HER2-low from HER2-zero breast cancer.

Patients and methods

We retrospectively analyzed 118 MRI cases, including 78 HER2-low and 40 HER2-zero patients confirmed by immunohistochemistry or fluorescence in situ hybridization. From each DCE-MRI case, 960 radiomic features were extracted. These features were screened and reduced using intraclass correlation coefficient, Mann-Whitney U test, and least absolute shrinkage to establish rad-scores. Logistic regression (LR) assessed the model's effectiveness in distinguishing HER2-low from HER2-zero. A clinicopathological MRI characteristic model was constructed using univariate and multivariate analysis, and a nomogram was developed combining rad-scores with significant MRI characteristics. Model performance was evaluated using the receiver operating characteristic (ROC) curve, and clinical benefit was assessed with decision curve analysis.

Results

The radiomics model, clinical model, and nomogram successfully distinguished between HER2-low and HER2-zero. The radiomics model showed excellent performance, with area under the curve (AUC) values of 0.875 in the training set and 0.845 in the test set, outperforming the clinical model (AUC = 0.691 and 0.672, respectively). HER2 status correlated with increased rad-score and Time Intensity Curve (TIC). The nomogram outperformed both models, with AUC, sensitivity, and specificity values of 0.892, 79.6%, and 82.8% in the training set, and 0.886, 83.3%, and 90.9% in the test set.

Conclusions

The DCE-MRI-based nomogram shows promising potential in differentiating HER2-low from HER2-zero status in breast cancer patients.

MRI-Based Radiomics Methods for Predicting Ki-67 Expression in Breast Cancer: A Systematic Review and Meta-analysis

RATIONALE AND OBJECTIVES

The purpose of this systematic review and meta-analysis was to assess the quality and diagnostic accuracy of MRI-based radiomics for predicting Ki-67 expression in breast cancer.

MATERIALS AND METHODS

A systematic literature search was performed to find relevant studies published in different databases, including PubMed, Web of Science, and Embase up until March 10, 2023. All papers were independently evaluated for eligibility by two reviewers. Studies that matched research questions and provided sufficient data for quantitative synthesis were included in the systematic review and meta-analysis, respectively. The quality of the articles was assessed using Quality Assessment of Diagnostic Accuracy Studies 2 (QUADAS-2) and Radiomics Quality Score (RQS) tools. The predictive value of MRI-based radiomics for Ki-67 antigen in patients with breast cancer was assessed using pooled sensitivity (SEN), specificity, and area under the curve (AUC). Meta-regression was performed to explore the cause of heterogeneity. Different covariates were used for subgroup analysis.

RESULTS

31 studies were included in the systematic review; among them, 21 reported sufficient data for meta-analysis. 20 training cohorts and five validation cohorts were pooled separately. The pooled sensitivity, specificity, and AUC of MRI-based radiomics for predicting Ki-67 expression in training cohorts were 0.80 [95% CI, 0.73-0.86], 0.82 [95% CI, 0.78-0.86], and 0.88 [95%CI, 0.85-0.91], respectively. The corresponding values for validation cohorts were 0.81 [95% CI, 0.72-0.87], 0.73 [95% CI, 0.62-0.82], and 0.84 [95%CI, 0.80-0.87], respectively. Based on QUADAS-2, some risks of bias were detected for reference standard and flow and timing domains. However, the quality of the included article was acceptable. The mean RQS score of the included articles was close to 6, corresponding to 16.6% of the maximum possible score. Significant heterogeneity was observed in pooled sensitivity and specificity of training cohorts (I2 > 75%). We found that using deep learning radiomic methods, magnetic field strength (3 T vs. 1.5 T), scanner manufacturer, region of interest structure (2D vs. 3D), route of tissue sampling, Ki-67 cut-off, logistic regression for model construction, and LASSO for feature reduction as well as PyRadiomics software for feature extraction had a great impact on heterogeneity according to our joint model analysis. Diagnostic performance in studies that used deep learning-based radiomics and multiple MRI sequences (e.g., DWI+DCE) was slightly higher. In addition, radiomic features derived from DWI sequences performed better than contrast-enhanced sequences in terms of specificity and sensitivity. No publication bias was found based on Deeks' funnel plot. Sensitivity analysis showed that eliminating every study one by one does not impact overall results.

CONCLUSION

This meta-analysis showed that MRI-based radiomics has a good diagnostic accuracy in differentiating breast cancer patients with high Ki-67 expression from low-expressing groups. However, the sensitivity and specificity of these methods still do not surpass 90%, restricting them from being used as a supplement to current pathological assessments (e.g., biopsy or surgery) to predict Ki-67 expression accurately.

Deconvolution-Based Pharmacokinetic Analysis to Improve the Prediction of Pathological Information of Breast Cancer

Pharmacokinetic (PK) parameters, revealing changes in the tumor microenvironment, are related to the pathological information of breast cancer. Tracer kinetic models (e.g., Tofts-Kety model) with a nonlinear least square solver are commonly used to estimate PK parameters. However, the method is sensitive to noise in images. To relieve the effects of noise, a deconvolution (DEC) method, which was validated on synthetic concentration-time series, was proposed to accurately calculate PK parameters from breast dynamic contrast-enhanced magnetic resonance imaging. A time-to-peak-based tumor partitioning method was used to divide the whole tumor into three tumor subregions with different kinetic patterns. Radiomic features were calculated from the tumor subregion and whole tumor-based PK parameter maps. The optimal features determined by the fivefold cross-validation method were used to build random forest classifiers to predict molecular subtypes, Ki-67, and tumor grade. The diagnostic performance evaluated by the area under the receiver operating characteristic curve (AUC) was compared between the subregion and whole tumor-based PK parameters. The results showed that the DEC method obtained more accurate PK parameters than the Tofts method. Moreover, the results showed that the subregion-based Ktrans (best AUCs = 0.8319, 0.7032, 0.7132, 0.7490, 0.8074, and 0.6950) achieved a better diagnostic performance than the whole tumor-based Ktrans (AUCs = 0.8222, 0.6970, 0.6511, 0.7109, 0.7620, and 0.5894) for molecular subtypes, Ki-67, and tumor grade. These findings indicate that DEC-based Ktrans in the subregion has the potential to accurately predict molecular subtypes, Ki-67, and tumor grade.

Efficient estimation of pharmacokinetic parameters from breast dynamic contrast-enhanced MRI based on a convolutional neural network for predicting molecular subtypes

Objective. Tracer kinetic models allow for estimating pharmacokinetic (PK) parameters, which are related to pathological characteristics, from breast dynamic contrast-enhanced magnetic resonance imaging. However, existing tracer kinetic models subject to inaccuracy are time-consuming for PK parameters estimation. This study aimed to accurately and efficiently estimate PK parameters for predicting molecular subtypes based on convolutional neural network (CNN).Approach. A CNN integrating global and local features (GL-CNN) was trained using synthetic data where known PK parameters map was used as the ground truth, and subsequently used to directly estimate PK parameters (volume transfer constantKtransand flux rate constantKep) map. The accuracy assessed by the peak signal-to-noise ratio (PSNR), structural similarity (SSIM), and concordance correlation coefficient (CCC) was compared between the GL-CNN and Tofts-based PK parameters in synthetic data. Radiomic features were calculated from the PK parameters map in 208 breast tumors. A random forest classifier was constructed to predict molecular subtypes using a discovery cohort (n= 144). The diagnostic performance evaluated on a validation cohort (n= 64) using the area under the receiver operating characteristic curve (AUC) was compared between the GL-CNN and Tofts-based PK parameters.Main results. The average PSNR (48.8884), SSIM (0.9995), and CCC (0.9995) between the GL-CNN-basedKtransmap and ground truth were significantly higher than those between the Tofts-basedKtransmap and ground truth. The GL-CNN-basedKtransobtained significantly better diagnostic performance (AUCs = 0.7658 and 0.8528) than the Tofts-basedKtransfor luminal B and HER2 tumors. The GL-CNN method accelerated the computation by speed approximately 79 times compared to the Tofts method for the whole breast of all patients.Significance. Our results indicate that the GL-CNN method can be used to accurately and efficiently estimate PK parameters for predicting molecular subtypes.

Artificial intelligence in breast imaging: potentials and challenges

Breast cancer, which is the most common type of malignant tumor among humans, is a leading cause of death in females. Standard treatment strategies, including neoadjuvant chemotherapy, surgery, postoperative chemotherapy, targeted therapy, endocrine therapy, and radiotherapy, are tailored for individual patients. Such personalized therapies have tremendously reduced the threat of breast cancer in females. Furthermore, early imaging screening plays an important role in reducing the treatment cycle and improving breast cancer prognosis. The recent innovative revolution in artificial intelligence (AI) has aided radiologists in the early and accurate diagnosis of breast cancer. In this review, we introduce the necessity of incorporating AI into breast imaging and the applications of AI in mammography, ultrasonography, magnetic resonance imaging, and positron emission tomography/computed tomography based on published articles since 1994. Moreover, the challenges of AI in breast imaging are discussed.

Radiomics characterization of tissues in an animal brain tumor model imaged using dynamic contrast enhanced (DCE) MRI

Here, we investigate radiomics-based characterization of tumor vascular and microenvironmental properties in an orthotopic rat brain tumor model measured using dynamic-contrast-enhanced (DCE) MRI. Thirty-two immune compromised-RNU rats implanted with human U-251N cancer cells were imaged using DCE-MRI (7Tesla, Dual-Gradient-Echo). The aim was to perform pharmacokinetic analysis using a nested model (NM) selection technique to classify brain regions according to vasculature properties considered as the source of truth. A two-dimensional convolutional-based radiomics analysis was performed on the raw-DCE-MRI of the rat brains to generate dynamic radiomics maps. The raw-DCE-MRI and respective radiomics maps were used to build 28 unsupervised Kohonen self-organizing-maps (K-SOMs). A Silhouette-Coefficient (SC), k-fold Nested-Cross-Validation (k-fold-NCV), and feature engineering analyses were performed on the K-SOMs' feature spaces to quantify the distinction power of radiomics features compared to raw-DCE-MRI for classification of different Nested Models. Results showed that eight radiomics features outperformed respective raw-DCE-MRI in prediction of the three nested models. The average percent difference in SCs between radiomics features and raw-DCE-MRI was: 29.875% ± 12.922%, p < 0.001. This work establishes an important first step toward spatiotemporal characterization of brain regions using radiomics signatures, which is fundamental toward staging of tumors and evaluation of tumor response to different treatments.

New Cluster Analysis Method for Quantitative Dynamic Contrast-Enhanced MRI Assessing Tumor Heterogeneity Induced by a Tumor-Microenvironmental Ameliorator (E7130) Treatment to a Breast Cancer Mouse Model

BACKGROUND

Dynamic contrast-enhanced magnetic resonance imaging (DCE-MRI) can provide insight into tumor perfusion. However, a method that can quantitatively measure the intra-tumor distribution of tumor voxel clusters with a distinct range of K^trans and v_e values remains insufficiently explored.

HYPOTHESIS

Two-dimensional cluster analysis may quantify the distribution of a tumor voxel subregion with a distinct range of K^trans and v_e values in human breast cancer xenografts.

STUDY TYPE

Prospective longitudinal study.

ANIMAL MODEL

Twenty-two female athymic nude mice with MCF-7 xenograft, treated with E7130, a tumor-microenvironmental ameliorator, or saline.

FIELD STRENGTH/SEQUENCE

9.4 Tesla, turbo rapid acquisition with relaxation enhancement, and spoiled gradient-echo sequences.

ASSESSMENT

We performed two-dimensional k-means clustering to identify tumor voxel clusters with a distinct range of K^trans and v_e values on Days 0, 2, and 5 after treatment, calculated the ratio of the number of tumor voxels in each cluster to the total number of tumor voxels, and measured the normalized distances defined as the ratio of the distance between each tumor voxel and the nearest tumor margin to a tumor radius.

STATISTICAL TESTS

Unpaired t-tests, Dunnett's multiple comparison tests, and Chi-squared test were used.

RESULTS

The largest and second largest clusters constituted 44.4% and 27.5% of all tumor voxels with cluster centroid values of K^trans at 0.040 min^-1 and 0.116 min^-1 , and v_e at 0.131 and 0.201, respectively. At baseline (Day 0), the average normalized distances for the largest and second largest clusters were 0.33 and 0.24, respectively. E7130-treated group showed the normalized distance of the initial largest cluster decreasing to 0.25, while that of the second largest cluster increasing to 0.31. Saline-treated group showed no change.

DATA CONCLUSION

A two-dimensional cluster analysis might quantify the spatial distribution of a tumor subregion with a distinct range of K^trans and v_e values.

LEVEL OF EVIDENCE

1 TECHNICAL EFFICACY: Stage 1.