Role of texture analysis in breast MRI as a cancer biomarker: A review

Texture analysis using non-contrast MRI to predict estrogen receptor status in breast cancer lesions: Is it feasible?

PURPOSE

We aimed to evaluate the feasibility of non-contrast-enhanced T1 sequence in texture analysis of breast cancer lesions to predict their estrogen receptor status.

METHODS

The study included 85 pathologically proven breast cancer lesions in 53 patients. Immunohistochemical studies were performed to determine the estrogen receptor status (ER). Lesions were divided into two groups: ER + ve status and ER-ve status. Texture analysis using the second-order analysis features [The Co-occurrence matrix (11 features)] was applied on both T1 and dynamic contrast-enhanced (DCE) MRI images for each lesion. Texture features gained from both T1 and DCE images were analyzed to obtain cut-off values using ROC curves to sort lesions according to their estrogen receptor status.

RESULTS

Angular second momentum and some of the entropy-based features showed statistically significant cut-off values in differentiation between the two groups [P-values for pre- and post-contrast images for AngSecMom (0.001, 0.008), sum entropy (0.003,0.005), and entropy (0.033,0.019), respectively]. On comparing the AUCs between pre- and post-contrast images, we found that differences were statistically insignificant. Sum of squares, sum variance, and sum average showed statistically significant cut-off points only on pre-contrast images [P-values for sum of squares (0.018), sum variance (0.024), and sum average (0.039)].

CONCLUSIONS

Texture analysis features showed promising results in predicting estrogen receptor status of breast cancer lesions on non-contrast T1 images.

Muscle texture features on preoperative MRI for diagnosis and assessment of severity of congenital muscular torticollis

OBJECTIVES

To develop an objective method based on texture analysis on MRI for diagnosis of congenital muscular torticollis (CMT).

MATERIAL AND METHODS

The T1- and T2-weighted imaging, Q-dixon, and T1-mapping MRI data of 38 children with CMT were retrospectively analyzed. The region of interest (ROI) was manually drawn at the level of the largest cross-sectional area of the SCM on the affected side. MaZda software was used to obtain the texture features of the T2WI sequences of the ROI in healthy and affected SCM. A radiomics diagnostic model based on muscle texture features was constructed using logistic regression analysis. Fatty infiltration grade was calculated by hematoxylin and eosin staining, and fibrosis ratio by Masson staining. Correlation between the MRI parameters and pathological indicators was analyzed.

RESULTS

There was positive correlation between fatty infiltration grade and mean value, standard deviation, and maximum value of the Q-dixon sequence of the affected SCM (correlation coefficients, 0.65, 0.59, and 0.58, respectively, P < 0.05).Three muscle texture features-S(2,2)SumAverg, S(3,3)SumVarnc, and T2WI extreme difference-were selected to construct the diagnostic model. The model showed significant diagnostic value for CMT (P < 0.05). The area under the curve of the multivariate conditional logistic regression model was 0.828 (95% confidence interval 0.735-0.922); the sensitivity was 0.684 and the specificity 0.868.

CONCLUSION

The radiomics diagnostic model constructed using T2WI muscle texture features and MRI signal values appears to have good diagnostic efficiency. Q-dixon sequence can reflect the fatty infiltration grade of CMT.

A multi-omics method for breast cancer diagnosis based on metabolites in exhaled breath, ultrasound imaging, and basic clinical information

Background and aims

Through a nested cohort study, we evaluated the diagnostic performance of breath-omics in differentiating between benign and malignant breast lesions, and assessed the diagnostic performance of a multi-omics approach that combines breath-omics, ultrasound radiomics, and clinic-omics in distinguishing between benign and malignant breast lesions.

Materials and methods

We recruited 1,723 consecutive patients who underwent an automated breast volume scanner (ABVS) examination. Breath samples were collected and analyzed by high-pressure photon ionization time-of-flight mass spectrometry (HPPI-TOF-MS) to obtain breath-omics features. 238 of 1,723 enrolled participants have received pathological confirmation of breast nodules finally. The breast lesions of the 238 participants were contoured manually based on ABVS images for ultrasound radiomics feature calculation. Then, single- and multi-omics models were constructed and evaluated for breast nodules diagnosis via five-fold cross-validation.

Results

The area under the curve (AUC) of the breath-omics model was 0.855. In comparison, the multi-omics model demonstrated superior diagnostic performance for breast cancer, with sensitivity, specificity, and AUC of 84.1 %, 89.9 %, and 0.946, respectively. The multi-omics performance was comparable to that of the Breast Imaging Reporting and Data System (BI-RADS) classification via senior ultrasound physician evaluation.

Conclusion

The multi-omics approach combining metabolites in exhaled breath, ultrasound imaging, and basic clinical information exhibits superior diagnostic performance and promises to be a non-invasive and reliable tool for breast cancer diagnosis.

Development and validation of radiomics nomogram for metastatic status of epithelial ovarian cancer

To develop and validate an enhanced CT-based radiomics nomogram for evaluating preoperative metastasis risk of epithelial ovarian cancer (EOC). One hundred and nine patients with histologically confirmed EOC were retrospectively enrolled. The volume of interest (VOI) was delineated in preoperative enhanced CT images, and 851 radiomics features were extracted. The radiomics features were selected by the least absolute shrinkage and selection operator (LASSO), and the rad-score was calculated using the formula of the radiomics label. A clinical model, radiomics model, and combined model were constructed using the logistic regression classification algorithm. Receiver operating characteristic (ROC) curve analysis and decision curve analysis (DCA) were used to evaluate the diagnostic performance of the models. Seventy-five patients (68.8%) were histologically confirmed to have metastasis. Eleven optimal radiomics features were retained by the LASSO algorithm to develop the radiomic model. The combined model for evaluating metastasis of EOC achieved area under the curve (AUC) values of 0.929 (95% CI 0.8593-0.9996) in the training cohort and 0.909 (95% CI 0.7921-1.0000) in the test cohort. To facilitate clinical use, a radiomic nomogram was built by combining the clinical characteristics with rad-score. The DCA indicated that the nomogram had the most significant net benefit when the threshold probability exceeded 15%, surpassing the benefits of both the treat-all and treat-none strategies. Compared with clinical model and radiomics model, the radiomics nomogram has the best diagnostic performance in evaluating EOC metastasis. The nomogram is a useful and convenient tool for clinical doctors to develop personalized treatment plans for EOC patients.

Radiogenomics Map-Based Molecular and Imaging Phenotypical Characterization in Localised Prostate Cancer Using Pre-Biopsy Biparametric MR Imaging

To create a radiogenomics map and evaluate the correlation between molecular and imaging phenotypes in localized prostate cancer (PCa), using radical prostatectomy histopathology as a reference standard. Radiomic features were extracted from T2-weighted (T2WI) and Apparent Diffusion Coefficient (ADC) images of clinically localized PCa patients (n = 15) across different Gleason score-based risk categories. DNA extraction was performed on formalin-fixed, paraffin-embedded (FFPE) samples. Gene expression analysis of androgen receptor expression, apoptosis, and hypoxia was conducted using the Chromosome Analysis Suite (ChAS) application and OSCHIP files. The relationship between gene expression alterations and textural features was assessed using Pearson's correlation analysis. Receiver operating characteristic (ROC) analysis was utilized to evaluate the predictive accuracy of the model. A significant correlation was observed between radiomic texture features and copy number variation (CNV) of genes associated with apoptosis, hypoxia, and androgen receptor (p-value ≤ 0.05). The identified radiomic features, including Sum Entropy ADC, Inverse Difference ADC, Sum Variance T2WI, Entropy T2WI, Difference Variance T2WI, and Angular Secondary Moment T2WI, exhibited potential for predicting cancer grade and biological processes such as apoptosis and hypoxia. Incorporating radiomics and genomics into a prediction model significantly improved the prediction of prostate cancer grade (clinically significant prostate cancer), yielding an AUC of 0.95. Radiomic texture features significantly correlate with genotypes for apoptosis, hypoxia, and androgen receptor expression in localised prostate cancer. Integration of these into the prediction model improved prediction accuracy of clinically significant prostate cancer.

A multi-instance tumor subtype classification method for small PET datasets using RA-DL attention module guided deep feature extraction with radiomics features

BACKGROUND

Positron emission tomography (PET) is extensively employed for diagnosing and staging various tumors, including liver cancer, lung cancer, and lymphoma. Accurate subtype classification of tumors plays a crucial role in formulating effective treatment plans for patients. Notably, lymphoma comprises subtypes like diffuse large B-cell lymphoma and Hodgkin's lymphoma, while lung cancer encompasses adenocarcinoma, small cell carcinoma, and squamous cell carcinoma. Similarly, liver cancer consists of subtypes such as cholangiocarcinoma and hepatocellular carcinoma. Consequently, the subtype classification of tumors based on PET images holds immense clinical significance. However, in clinical practice, the number of cases available for each subtype is often limited and imbalanced. Therefore, the primary challenge lies in achieving precise subtype classification using a small dataset.

METHOD

This paper presents a novel approach for tumor subtype classification in small datasets using RA-DL (Radiomics-DeepLearning) attention. To address the limited sample size, Support Vector Machines (SVM) is employed as the classifier for tumor subtypes instead of deep learning methods. Emphasizing the importance of texture information in tumor subtype recognition, radiomics features are extracted from the tumor regions during the feature extraction stage. These features are compressed using an autoencoder to reduce redundancy. In addition to radiomics features, deep features are also extracted from the tumors to leverage the feature extraction capabilities of deep learning. In contrast to existing methods, our proposed approach utilizes the RA-DL-Attention mechanism to guide the deep network in extracting complementary deep features that enhance the expressive capacity of the final features while minimizing redundancy. To address the challenges of limited and imbalanced data, our method avoids using classification labels during deep feature extraction and instead incorporates 2D Region of Interest (ROI) segmentation and image reconstruction as auxiliary tasks. Subsequently, all lesion features of a single patient are aggregated into a feature vector using a multi-instance aggregation layer.

RESULT

Validation experiments were conducted on three PET datasets, specifically the liver cancer dataset, lung cancer dataset, and lymphoma dataset. In the context of lung cancer, our proposed method achieved impressive performance with Area Under Curve (AUC) values of 0.82, 0.84, and 0.83 for the three-classification task. For the binary classification task of lymphoma, our method demonstrated notable results with AUC values of 0.95 and 0.75. Moreover, in the binary classification task of liver tumor, our method exhibited promising performance with AUC values of 0.84 and 0.86.

CONCLUSION

The experimental results clearly indicate that our proposed method outperforms alternative approaches significantly. Through the extraction of complementary radiomics features and deep features, our method achieves a substantial improvement in tumor subtype classification performance using small PET datasets.

Diagnosis of benign and malignant nodules with a radiomics model integrating features from nodules and mammary regions on DCE-MRI

Objectives

To establish a radiomics model for distinguishing between the benign and malignant mammary gland nodules via combining the features from nodule and mammary regions on DCE-MRI.

Methods

In this retrospective study, a total of 103 cases with mammary gland nodules (malignant/benign = 80/23) underwent DCE-MRI, and was confirmed by biopsy pathology. Features were extracted from both nodule region and mammary region on DCE-MRI. Three SVM classifiers were built for diagnosis of benign and malignant nodules as follows: the model with the features only from nodule region (N model), with the features only from mammary region (M model) and the model combining the features from nodule region and mammary region (NM model). The performance of models was evaluated with the area under the curve of receiver operating characteristic (AUC).

Results

One radiomic features is selected from nodule region and 3 radiomic features is selected from mammary region. Compared with N or M model, NM model exhibited the best performance with an AUC of 0.756.

Conclusions

Compared with the model only using the features from nodule or mammary region, the radiomics-based model combining the features from nodule and mammary region outperformed in the diagnosis of benign and malignant nodules.

Advancing Diagnostics and Patient Care: The Role of Biomarkers in Radiology

The integration of biomarkers into medical practice has revolutionized the field of radiology, allowing for enhanced diagnostic accuracy, personalized treatment strategies, and improved patient care outcomes. This review offers radiologists a comprehensive understanding of the diverse applications of biomarkers in medicine. By elucidating the fundamental concepts, challenges, and recent advancements in biomarker utilization, it will serve as a bridge between the disciplines of radiology and epidemiology. Through an exploration of various biomarker types, such as imaging biomarkers, molecular biomarkers, and genetic markers, I outline their roles in disease detection, prognosis prediction, and therapeutic monitoring. I also discuss the significance of robust study designs, blinding, power and sample size calculations, performance metrics, and statistical methodologies in biomarker research. By fostering collaboration between radiologists, statisticians, and epidemiologists, I hope to accelerate the translation of biomarker discoveries into clinical practice, ultimately leading to improved patient care.

Detecting Double Expression Status in Primary Central Nervous System Lymphoma Using Multiparametric MRI Based Machine Learning

BACKGROUND

Double expression lymphoma (DEL) is a subtype of primary central nervous system lymphoma (PCNSL) that often has a poor prognosis. Currently, there are limited noninvasive ways to detect protein expression.

PURPOSE

To detect DEL in PCNSL using multiparametric MRI-based machine learning.

STUDY TYPE

Retrospective.

POPULATION

Forty PCNSL patients were enrolled in the study among whom 17 were DEL (9 males and 8 females, 61.29 ± 14.14 years) and 23 were non-DEL (14 males and 9 females, 55.57 ± 14.16 years) with 59 lesions (28 DEL and 31 non-DEL).

FIELD STRENGTH/SEQUENCE

ADC map derived from DWI (b = 0/1000 s/mm2 ), fast spin echo T2WI, T2FLAIR, and contrast-enhanced T1 weighted imaging (T1CE) were collected at 3.0 T.

ASSESSMENT

Two raters manually segmented lesions by ITK-SNAP on ADC, T2WI, T2FLAIR and T1CE. A total of 2234 radiomics features from the tumor segmentation area were extracted. The t-test was conducted to filter the features, and elastic net regression algorithm combined with recursive feature elimination was used to calculate the essential features. Finally, 12 groups with combinations of different sequences were fitted to 6 classifiers, and the optimal models were selected.

STATISTICAL TESTS

Continuous variables were assessed by the t-test, while categorical variables were assessed by the non-parametric test. Interclass correlation coefficient tested variables' consistency. Sensitivity, specificity, accuracy F1-score, and area under the curve (AUC) were used to evaluate model performance.

RESULTS

DEL status could be identified to varying degrees with 72 models based on radiomics, and model performance could be improved by combining different sequences and classifiers. Both SVMlinear and logistic regression (LR) combined with four sequence group had similar largest AUCmean (0.92 ± 0.09 vs. 0.92 ± 0.05), and SVMlinear was considered as the optimal model in this study since the F1-score of SVMlinear (0.88) was higher than that of LR (0.83).

DATA CONCLUSION

Multiparametric MRI-based machine learning is promising in DEL detection.

EVIDENCE LEVEL

4 TECHNICAL EFFICACY STAGE: 2.

Assessment and risk prediction of frailty using texture-based muscle ultrasound image analysis and machine learning techniques

The purpose of this study was to evaluate texture-based muscle ultrasound image analysis for the assessment and risk prediction of frailty phenotype. This retrospective study of prospectively acquired data included 101 participants who underwent ultrasound scanning of the anterior thigh. Participants were subdivided according to frailty phenotype and were followed up for two years. Primary and secondary outcome measures were death and comorbidity, respectively. Forty-three texture features were computed from the rectus femoris and the vastus intermedius muscles using statistical methods. Model performance was evaluated by computing the area under the receiver operating characteristic curve (AUC) while outcome prediction was evaluated using regression analysis. Models developed achieved a moderate to good AUC (0.67 ≤ AUC ≤ 0.79) for categorizing frailty. The stepwise multiple logistic regression analysis demonstrated that they correctly classified 70-87% of the cases. The models were associated with increased comorbidity (0.01 ≤ p ≤ 0.18) and were predictive of death for pre-frail and frail participants (0.001 ≤ p ≤ 0.016). In conclusion, texture analysis can be useful to identify frailty and assess risk prediction (i.e. mortality) using texture features extracted from muscle ultrasound images in combination with a machine learning approach.

Radiomics Based on DCE-MRI for Predicting Response to Neoadjuvant Therapy in Breast Cancer

RATIONALE AND OBJECTIVES

To compare the value of radiomics and diameter% based on pre- and early-treatment dynamic enhanced MR (DCE-MRI) of the breast in predicting response to neoadjuvant therapy (NAT) in breast cancer and to construct a tool for early noninvasive prediction of NAT outcomes.

MATERIALS AND METHODS

Retrospective analysis of clinical and imaging data of 142 patients with primary invasive breast cancer who underwent DCE-MRI before and after two cycles of NAT at our institution. Enroled patients were randomly assigned in a 7:3 ratio to the training group and the test group. Patients were divided into pathological complete response (pCR) and non-pathological complete response groups based on surgical pathology findings after NAT. The maximum diameter relative regression values (Diameter%) before and after treatment were calculated and the conventional imaging Diameter% model was constructed. Based on pre- and early-NAT DCE-MRI, the optimal features of pre-NAT, early-NAT, and delta radiomics were screened using redundancy analysis, least absolute shrinkage, and selection operator methods to construct the corresponding radiomics model and calculate the Radscores. Indicators that were statistically significant in the univariate analysis of clinical data were further screened by stepwise regression and combined with Radscores to construct the fusion model. All models were evaluated and compared.

RESULTS

In the test set, the area under the curve (AUC) of the delta radiomics model (0.87) was higher than that of the pre-NAT, early-NAT radiomics models (0.57, 0.78) and the Diameter% model (0.83). The fusion model had the best efficacy in predicting pCR after NAT, with AUCs of 0.91 in the training and test sets. And its nomogram plot showed that Radscore of early-NAT radiomics had the greatest weight. In the test set, the fusion model and Delta radiomics model improved the efficacy of predicting pCR by 35.56% and 14.19%, respectively, compared to the Diameter% model (P = 0 and .039). Clinical decision curves showed the highest overall clinical benefit for the fusion model.

CONCLUSION

Radiomics, especially delta and early-NAT radiomics, may be potential biomarkers for early noninvasive prediction of NAT outcomes. And a fusion model constructed from meaningful clinicopathological indicators combined with radiomics can effectively predict NAT response.

Development of a prediction model for gross residual in high-grade serous ovarian cancer by combining preoperative assessments of abdominal and pelvic metastases and multiparametric MRI

RATIONALE AND OBJECTIVES

To preoperatively predict residual tumor (RT) in patients with high-grade serous ovarian carcinoma (HGSOC) via a radiomic-clinical nomogram.

METHODS

A total of 128 patients with advanced HGSOC were enrolled (training cohort: n=106; validation cohort: n=22). Serum cancer antigen-125 (CA125), serum human epididymis protein 4 (HE-4) level, and neutrophil-to-lymphocyte ratio (NLR) were obtained from the medical records. Metastases in abdomen and pelvis (MAP) of HGSOC patients was evaluated and scored based on preoperative abdominal and pelvic enhanced CT, MRI and/or PET-CT. A volume of interest (VOI) of each tumor was manually contoured along the boundary slice-by-slice. Radiomic features were extracted from the T2-weighted imaging (T2WI), diffusion-weighted imaging (DWI) and apparent diffusion coefficient (ADC) images. Univariate and multivariate analyses were used to determine the independent predictors of RT status. Least absolute shrinkage and selection operator (LASSO) logistic regression was performed to select optimal features and construct radiomic models. A radiomic-clinical nomogram incorporating radiomic signature and clinical parameters was developed and evaluated in training and validation cohorts.

RESULTS

MAP score (p = 0.002), HE-4 level (p = 0.001) and NLR (p = 0.008) were independent predictors of RT status. The final radiomic-clinical nomogram showed satisfactory prediction performance in training (AUC = 0.936), cross validation (AUC = 0.906) and separate validation cohorts (AUC = 0.900), and fitted well in calibration curves (p > 0.05). Decision curve further confirmed the clinical application value of the nomogram.

CONCLUSION

The proposed MRI-based radiomic-clinical nomogram achieved excellent preoperative prediction of the RT status in HGSOC.

Nanotechnology revolutionises breast cancer treatment: harnessing lipid-based nanocarriers to combat cancer cells

One of the most common cancers that occur in females is breast cancer. Despite the significant leaps and bounds that have been made in treatment of breast cancer, the disease remains one of the leading causes of death among women and a major public health challenge. The therapeutic efficacy of chemotherapeutics is hindered by chemoresistance and toxicity. Nano-based lipid drug delivery systems offer controlled drug release, nanometric size and site-specific targeting. Breast cancer treatment includes surgery, chemotherapy and radiotherapy. Despite this, no single method of treatment for the condition is currently effective due to cancer stem cell metastasis and chemo-resistance. Therefore, the employment of nanocarrier systems is necessary in order to target breast cancer stem cells. This article addresses breast cancer treatment options, including modern treatment procedures such as chemotherapy, etc. and some innovative therapeutic options highlighting the role of lipidic nanocarriers loaded with chemotherapeutic drugs such as nanoemulsion, solid-lipid nanoparticles, nanostructured lipid carriers and liposomes, and their investigations have demonstrated that they can limit cancer cell growth, reduce the risk of recurrence, as well as minimise post-chemotherapy metastasis. This article also explores FDA-approved lipid-based nanocarriers, commercially available formulations, and ligand-based formulations that are being considered for further research.

Can magnetic resonance imaging texture analysis change the breast imaging reporting and data system category of breast lesions?

PURPOSE

This study aimed to reveal magnetic resonance imaging (MRI) texture analysis (TA)'s contribution to categorizing breast lesions according to the Breast Imaging-Reporting and Data System (BI-RADS) lexicon.

METHOD

Two hundred and seventeen women with BI-RADS category 3, 4, and 5 lesions on breast MRI were included in the study. For TA, the region of interest was drawn manually to encompass the entire lesion on the fat-suppressed T2W and the first post-contrast T1W images. To identify the independent predictors of breast cancer, multivariate logistic regression analyses were performed using texture parameters. Estimated benign and malignant groups were formed according to the TA regression model.

RESULTS

Texture parameters extracted from T2WI, including median, gray-level co-occurrence matrix (GLCM) contrast, GLCM correlation, GLCM joint entropy, GLCM sum entropy, and GLCM sum of squares, and parameters extracted from T1WI, including maximum, GLCM contrast, GLCM joint entropy, GLCM sum entropy, were independent predictors of breast cancer. In the estimated new groups according to the TA regression model, 19 (91%) of the benign 4a lesions were downgraded to BI-RADS category 3.

CONCLUSIONS

The addition of quantitative parameters obtained by MRI TA to BI-RADS criteria significantly increased the accuracy rate in differentiating benign and malignant breast lesions. When categorizing BI-RADS 4a lesions, the use of MRI TA in addition to conventional imaging findings may reduce unnecessary biopsy rates.

Improvement of semantic segmentation through transfer learning of multi-class regions with convolutional neural networks on supine and prone breast MRI images

Semantic segmentation of breast and surrounding tissues in supine and prone breast magnetic resonance imaging (MRI) is required for various kinds of computer-assisted diagnoses for surgical applications. Variability of breast shape in supine and prone poses along with various MRI artifacts makes it difficult to determine robust breast and surrounding tissue segmentation. Therefore, we evaluated semantic segmentation with transfer learning of convolutional neural networks to create robust breast segmentation in supine breast MRI without considering supine or prone positions. Total 29 patients with T1-weighted contrast-enhanced images were collected at Asan Medical Center and two types of breast MRI were performed in the prone position and the supine position. The four classes, including lungs and heart, muscles and bones, parenchyma with cancer, and skin and fat, were manually drawn by an expert. Semantic segmentation on breast MRI scans with supine, prone, transferred from prone to supine, and pooled supine and prone MRI were trained and compared using 2D U-Net, 3D U-Net, 2D nnU-Net and 3D nnU-Net. The best performance was 2D models with transfer learning. Our results showed excellent performance and could be used for clinical purposes such as breast registration and computer-aided diagnosis.

The use of radiomics and machine learning for the differentiation of chondrosarcoma from enchondroma

PURPOSE

To construct and compare machine learning models for differentiating chondrosarcoma from enchondroma using radiomic features from T1 and fat suppressed Proton density (PD) magnetic resonance imaging (MRI).

METHODS

Eighty-eight patients (57 with enchondroma, 31 with chondrosarcoma) were retrospectively included. Histogram matching and N4ITK MRI bias correction filters were applied. An experienced musculoskeletal radiologist and a senior resident in radiology performed manual segmentation. Voxel sizes were resampled. Laplacian of Gaussian filter and wavelet-based features were used. One thousand eight hundred eighty-eight features were obtained for each patient, with 944 from T1 and 944 from PD images. Sixty-four unstable features were removed. Seven machine learning models were used for classification.

RESULTS

Classification with all features showed neural network was the best model for both readers' datasets with area under the curve (AUC), classification accuracy (CA), and F1 score of 0.979, 0.984; 0.920, 0.932; and 0.889, 0.903, respectively. Four features, including one common to both readers, were selected using fast correlation based filter. The best performing models with selected features were gradient boosting for Fatih Erdem's dataset and neural network for Gülen Demirpolat's dataset with AUC, CA, and F1 score of 0.990, 0.979; 0.943, 0.955; 0.921, 0.933, respectively. Neural Network was the second-best model for FE's dataset based on AUC (0.984).

CONCLUSION

Using pathology as a gold standard, this study defined and compared seven well-performing models to distinguish enchondromas from chondrosarcomas and provided radiomic feature stability and reproducibility among the readers.

Radiogenomics Reveals Correlation between Quantitative Texture Radiomic Features of Biparametric MRI and Hypoxia-Related Gene Expression in Men with Localised Prostate Cancer

OBJECTIVES

To perform multiscale correlation analysis between quantitative texture feature phenotypes of pre-biopsy biparametric MRI (bpMRI) and targeted sequence-based RNA expression for hypoxia-related genes.

MATERIALS AND METHODS

Images from pre-biopsy 3T bpMRI scans in clinically localised PCa patients of various risk categories (n = 15) were used to extract textural features. The genomic landscape of hypoxia-related gene expression was obtained using post-radical prostatectomy tissue for targeted RNA expression profiling using the TempO-sequence method. The nonparametric Games Howell test was used to correlate the differential expression of the important hypoxia-related genes with 28 radiomic texture features. Then, cBioportal was accessed, and a gene-specific query was executed to extract the Oncoprint genomic output graph of the selected hypoxia-related genes from The Cancer Genome Atlas (TCGA). Based on each selected gene profile, correlation analysis using Pearson's coefficients and survival analysis using Kaplan-Meier estimators were performed.

RESULTS

The quantitative bpMR imaging textural features, including the histogram and grey level co-occurrence matrix (GLCM), correlated with three hypoxia-related genes (ANGPTL4, VEGFA, and P4HA1) based on RNA sequencing using the TempO-Seq method. Further radiogenomic analysis, including data accessed from the cBioportal genomic database, confirmed that overexpressed hypoxia-related genes significantly correlated with a poor survival outcomes, with a median survival ratio of 81.11:133.00 months in those with and without alterations in genes, respectively.

CONCLUSION

This study found that there is a correlation between the radiomic texture features extracted from bpMRI in localised prostate cancer and the hypoxia-related genes that are differentially expressed. The analysis of expression data based on cBioportal revealed that these hypoxia-related genes, which were the focus of the study, are linked to an unfavourable survival outcomes in prostate cancer patients.

Delta computed tomography radiomics features-based nomogram predicts long-term efficacy after neoadjuvant chemotherapy in advanced gastric cancer

BACKGROUND AND OBJECTIVE

No effective preoperative tool is available for predicting the prognosis of advanced gastric cancer (AGC) treated by neoadjuvant chemotherapy (NAC). We aimed to explore the association between change values ("delta") in the radiomic signatures of computed tomography (CT) (delCT-RS) before and after NAC for AGC and overall survival(OS).

METHODS AND DESIGN

A total of 132 AGC patients with AGC were studied as a training cohort in our center, and 45 patients from another center were used as an external validation set. A radiomic signatures-clinical-nomogram(RS-CN) was established using delCT-RS and preoperative clinical variables. The prediction performance of RS-CN was evaluated using the area under the receiver operating characteristic (ROC)curve (AUC values), time-dependent ROC, decision curve analysis(DCA) and C-index.

RESULTS

Multivariable Cox regression analyses showed that delCT-RS, cT-stage, cN-stage, Lauren-type and the value of variation of carcinoma embryonic antigen (CEA) between NAC were independent risk factors for 3-year OS of AGC. In the training cohort, RS-CN had a good prediction performance for OS (C-Index 0.73) and AUC values were significantly better than those of delCT-RS, ypTNM-stage and tumor regression grade(TRG) (0.827 vs 0.704 vs 0.749 vs 0.571, p < 0.001). DCA and time-dependent ROC of RS-CN were better than those of ypTNM stage, TRG grade and delCT-RS. The prediction performance of the validation set was equivalent to that of the training set. The cut-off (177.2) of RS-CN score was obtained from X-Tile software, a score of > 177.2 was defined as high-risk group(HRG), and scores of ≤ 177.2 were defined as the low-risk group(LRG). The 3-year OS and disease free survival(DFS) of patients in the LRG were significantly better than those in the HRG. Adjuvant chemotherapy(AC) can only significantly improve the 3-year OS and DFS of the LRG. (p < 0.05).

CONCLUSIONS

Our nomogram based on delCT-RS has good prediction of prognosis before surgery and helps identify patients that are most likely to benefit from AC. It works well in precise and individualised NAC in AGC.

Four-Dimensional Machine Learning Radiomics for the Pretreatment Assessment of Breast Cancer Pathologic Complete Response to Neoadjuvant Chemotherapy in Dynamic Contrast-Enhanced MRI

BACKGROUND

Breast cancer response to neoadjuvant chemotherapy (NAC) is typically evaluated through the assessment of tumor size reduction after a few cycles of NAC. In case of treatment ineffectiveness, this results in the patient suffering potentially severe secondary effects without achieving any actual benefit.

PURPOSE

To identify patients achieving pathologic complete response (pCR) after NAC by spatio-temporal radiomic analysis of dynamic contrast-enhanced (DCE) MRI images acquired before treatment.

STUDY TYPE

Single-center, retrospective.

POPULATION

A total of 251 DCE-MRI pretreatment images of breast cancer patients.

FIELD STRENGTH/SEQUENCE

1.5 T/3 T, T1-weighted DCE-MRI.

ASSESSMENT

Tumor and peritumoral regions were segmented, and 348 radiomic features that quantify texture temporal variation, enhancement kinetics heterogeneity, and morphology were extracted. Based on subsets of features identified through forward selection, machine learning (ML) logistic regression models were trained separately with all images and stratifying on cancer molecular subtype and validated with leave-one-out cross-validation.

STATISTICAL TESTS

Feature significance was assessed using the Mann-Whitney U-test. Significance of the area under the receiver operating characteristics (ROC) curve (AUC) of the ML models was assessed using the associated 95% confidence interval (CI). Significance threshold was set to 0.05, adjusted with Bonferroni correction.

RESULTS

Nine features related to texture temporal variation and enhancement kinetics heterogeneity were significant in the discrimination of cases achieving pCR vs. non-pCR. The ML models achieved significant AUC of 0.707 (all cancers, n = 251, 59 pCR), 0.824 (luminal A, n = 107, 14 pCR), 0.823 (luminal B, n = 47, 15 pCR), 0.844 (HER2 enriched, n = 25, 11 pCR), 0.803 (triple negative, n = 72, 19 pCR).

DATA CONCLUSIONS

Differences in imaging phenotypes were found between complete and noncomplete responders. Furthermore, ML models trained per cancer subtype achieved high performance in classifying pCR vs. non-pCR cases. They may, therefore, have potential to help stratify patients according to the level of response predicted before treatment, pending further validation with larger prospective cohorts.

EVIDENCE LEVEL

4 TECHNICAL EFFICACY: Stage 4.

Clinical-Radiomics Nomogram from T1W, T1CE, and T2FS MRI for Improving Diagnosis of Soft-Tissue Sarcoma

PURPOSE

To compare values of multiparametric magnetic resonance imaging (MRI) sequences and propose clinical-radiomics nomogram for diagnosis of soft-tissue sarcoma (STS).

PROCEDURES

This study enrolled 148 patients from Dec. 2017 to Feb. 2021. All patients underwent T1-weighted (T1W), contrast-enhanced T1-weighted (T1CE), and T2-weighted fat-suppressed (T2FS) MRI scans. A total of 1967 radiomic features were extracted from the segmented regions of interest (ROIs) in each MRI sequence. Highly diagnostic radiomic features were selected with Mann-Whitney U test, elastic net, and Akaike's information criterion (AIC) based on MRI images. Logistical regression was used to build Rad scores. Clinical factors were analyzed using the chi-square test or Mann-Whitney U test. The performance of the Rad scores was judged using the area under the receiver operating characteristic area under the curve (ROC AUC), sensitivity, specificity, and accuracy. The nomogram was developed by integrating the Rad score and the most important clinical factor.

RESULTS

By combining the three MRI sequences, the Rad-Com was developed consisting of twelve features selected by with Mann-Whitney U test, elastic net, and AIC: four from T1W, three from TICE, and five from T2FS MRI. The margin (P < 0.05) demonstrated a statistically significant difference between patients with benign and malignant soft-tissue tumors (STT). The nomogram was constructed by integrating the Rad-Com and margin, which yielded favorable diagnostic AUCs of 0.919 (sensitivity (Sen) = 0.784, specificity (Spe) = 0.936) and 0.913 (Sen = 0.923, Spe = 0.792) in the training and validation cohort.

CONCLUSION

The proposed nomogram may have potential as a noninvasive marker for STS diagnosis.

MRI-based delta-radiomic features for prediction of local control in liver lesions treated with stereotactic body radiation therapy

Real-time magnetic resonance image guided stereotactic ablative radiotherapy (MRgSBRT) is used to treat abdominal tumors. Longitudinal data is generated from daily setup images. Our study aimed to identify delta radiomic texture features extracted from these images to predict for local control in patients with liver tumors treated with MRgSBRT. Retrospective analysis of an IRB-approved database identified patients treated with MRgSBRT for primary liver and secondary metastasis histologies. Daily low field strength (0.35 T) images were retrieved, and the gross tumor volume was identified on each image. Next, images' gray levels were equalized, and 39 s-order texture features were extracted. Delta-radiomics were calculated as the difference between feature values on the initial scan and after delivered biological effective doses (BED, α/β = 10) of 20 Gy and 40 Gy. Then, features were ranked by the Gini Index during training of a random forest model. Finally, the area under the receiver operating characteristic curve (AUC) was estimated using a bootstrapped logistic regression with the top two features. We identified 22 patients for analysis. The median dose delivered was 50 Gy in 5 fractions. The top two features identified after delivery of BED 20 Gy were gray level co-occurrence matrix features energy and gray level size zone matrix based large zone emphasis. The model generated an AUC = 0.9011 (0.752-1.0) during bootstrapped logistic regression. The same two features were selected after delivery of a BED 40 Gy, with an AUC = 0.716 (0.600-0.786). Delta-radiomic features after a single fraction of SBRT predicted local control in this exploratory cohort. If confirmed in larger studies, these features may identify patients with radioresistant disease and provide an opportunity for physicians to alter management much sooner than standard restaging after 3 months. Expansion of the patient database is warranted for further analysis of delta-radiomic features.

Fractal and textural imaging identify new subgroups of patients with colorectal cancer based on biophysical properties of the cancer cells

Colorectal cancer (CRC) can been sub-divided, based on the generation of tertiary lymphoid structures (TLS), into CRC with a Crohn's like lymphoid reaction (CLR) representing de novo formation of TLSs or CRC lacking TLSs that show Diffuse Inflammatory infiltration (DII). The association between TLS, early treatment initiation and longer survival highlights the need for deeper patient stratification that could lead to more targeted therapies. We hypothesized that such stratification might be achieved by using digital image analyses. Here we retrospectively analyzed 35 CRC patient samples classified as CLR or DII by digital analysis, focusing on the parameters Fractal dimension, Lacunarity and the textural features Angular second momentum, Correlation, Inverse difference momentum and Entropy. Significant differences in the grades of these parameters between the two patient groups provided preliminary data that additional biophysical information can divide CRC into at least 3 subgroups which encompass CLR and DII. Additional studies are needed to test if this sub-classification aids in the selection of targeted therapy for patients with CRC.

Combining Clinicopathology, IVIM-DWI and Texture Parameters for a Nomogram to Predict Treatment Response to Neoadjuvant Chemoradiotherapy in Locally Advanced Rectal Cancer Patients

Objectives

This study aimed to create a nomogram for the risk prediction of neoadjuvant chemoradiotherapy (nCRT) resistance in locally advanced rectal cancer (LARC).

Methods

Clinical data in this retrospective study were collected from a total of 135 LARC patients admitted to our hospital from June 2016 to December 2020. After screening by inclusion and exclusion criteria, 62 patients were included in the study. Texture analysis (TA) was performed on T2WI and DWI images. Patients were divided into response group (CR+PR) and no-response group (SD+PD) according to efficacy assessment. Multivariate analysis was performed on clinicopathology, IVIM-DWI and texture parameters for screening of independent predictors. A nomogram was created and model fit and clinical net benefit were assessed.

Results

Multivariate analysis of clinicopathology parameters showed that the differentiation and T stage were independent predictors (OR values were 14.516 and 11.589, resp.; P<0.05). Multivariate analysis of IVIM-DWI and texture parameters showed that f value and Rads-score were independent predictors (OR values were 0.855, 2.790, resp.; P<0.05). In this study, clinicopathology together with IVIM-DWI and texture parameters showed the best predictive efficacy (AUC=0.979). The nomogram showed good predictive performance and stability in identifying high-risk LARC patients who are resistant to nCRT (C-index=0.979). Decision curve analyses showed that the nomogram had the best clinical net benefit. Ten-fold cross-validation results showed that the average AUC value was 0.967, and the average C-index was 0.966.

Conclusions

The nomogram combining the differentiation, T stage, f value and Rads-score can effectively estimate the risk of nCRT resistance in patients with LARC.

Association between breast cancer's prognostic factors and 3D textural features of non-contrast-enhanced T1 weighted breast MRI

OBJECTIVES

The aim of this exploratory study was to evaluate whether three-dimensional texture analysis (3D-TA) features of non-contrast-enhanced T₁ weighted MRI associate with traditional prognostic factors and disease-free survival (DFS) of breast cancer.

METHODS

3D-T₁ weighted images from 78 patients with 81 malignant histopathologically verified breast lesions were retrospectively analysed using standard-size volumes of interest. Grey-level co-occurrence matrix (GLCM)-based features were selected for statistical analysis. In statistics the Mann-Whitney U and the Kruskal-Wallis tests, the Cox proportional hazards model and the Kaplan-Meier method were used.

RESULTS

Tumours with higher histological grade were significantly associated with higher contrast (1 voxel: p = 0.033, 2 voxels: p = 0.036). All the entropy parameters showed significant correlation with tumour grade (p = 0.015-0.050) but there were no statistically significant associations between other TA parameters and tumour grade. The Nottingham Prognostic Index (NPI) was correlated with contrast and sum entropy parameters. A higher sum variance TA parameter was a significant predictor of shorter DFS.

CONCLUSION

Texture parameters, assessed by 3D-TA from non-enhanced T₁ weighted images, indicate tumour heterogeneity but have limited independent prognostic value. However, they are associated with tumour grade, NPI, and DFS. These parameters could be used as an adjunct to contrast-enhanced TA parameters.

ADVANCES IN KNOWLEDGE

3D-TA of non-contrast enhanced T₁ weighted breast MRI associates with tumour grade, NPI, and DFS. The use of non-contrast 3D-TA parameters in adjunct with contrast-enhanced 3D-TA parameters warrants further research.

Imaging Tool for Predicting Renal Clear Cell Carcinoma Fuhrman Grade: Comparing R.E.N.A.L. Nephrometry Score and CT Texture Analysis

Background

Clear cell renal cell carcinoma (ccRCC) is the most common renal malignant tumor. Preoperative imaging boasts advantages in diagnosing and choosing treatment methods for ccRCC.

Purpose

This study is aimed at building models based on R.E.N.A.L. nephrometry score (RNS) and CT texture analysis (CTTA) to estimate the Fuhrman grade of ccRCC and comparing the advantages and disadvantages of the two models.

Materials and Methods

143 patients with pathologically confirmed ccRCC were enrolled. All patients were stratified into Fuhrman low-grade and high-grade groups with complete CT data and R.E.N.A.L. nephrometry scores. CTTA features were extracted from the ROI delineated at the largest tumor level, and RNS and CTTA features were included in the logistic regression model, respectively.

Results

RNS model constructed based on multivariate logistic regression analysis showed that 3 pts for R-scores, 2 pts for E-scores, and 3 pts for L-scores were significant indicators to predict high-grade ccRCC, the AUC of RNS model was 0.911, and the sensitivity and specificity were 71.11% and 83.67%, respectively. The CTTA-model confirmed energy, kurtosis, and entropy as independent predictive factors, and the AUC of CTTA model was 0.941, with an optimal sensitivity and specificity of 84.44% and 93.88%.

Conclusions

R.E.N.A.L. nephrometry score has a certain provocative effect on the Fuhrman pathological grading of ccRCC. As a potential emerging technology, CTTA is expected to replace R.E.N.A.L. nephrometry score in evaluating patients' Fuhrman classification, and this approach might become an available method for assisting clinicians in choosing appropriate operation.

MRI Radiomics for Assessment of Molecular Subtype, Pathological Complete Response, and Residual Cancer Burden in Breast Cancer Patients Treated With Neoadjuvant Chemotherapy

RATIONALE AND OBJECTIVES

There are limited data on pretreatment imaging features that can predict response to neoadjuvant chemotherapy (NAC). To extract volumetric pretreatment MRI radiomics features and assess corresponding associations with breast cancer molecular subtypes, pathological complete response (pCR), and residual cancer burden (RCB) in patients treated with NAC.

MATERIALS AND METHODS

In this IRB-approved study, clinical and pretreatment MRI data from patients with biopsy-proven breast cancer who received NAC between September 2009 and July 2016 were retrospectively analyzed. Tumors were manually identified and semi-automatically segmented on first postcontrast images. Morphological and three-dimensional textural features were computed, including unfiltered and filtered image data, with spatial scaling factors (SSF) of 2, 4, and 6 mm. Wilcoxon rank-sum tests and area under the receiver operating characteristic curve were used for statistical analysis.

RESULTS

Two hundred and fifty nine patients with unilateral breast cancer, including 73 (28.2%) HER2+, 112 (43.2%) luminal, and 74 (28.6%) triple negative breast cancers (TNBC), were included. There was a significant difference in the median volume (p = 0.008), median longest axial tumor diameter (p = 0.009), and median longest volumetric diameter (p = 0.01) among tumor subtypes. There was also a significant difference in minimum signal intensity and entropy among the tumor subtypes with SSF = 4 mm (p = 0.009 and p = 0.02 respectively) and SSF = 6 mm (p = 0.007 and p < 0.001 respectively). Additionally, sphericity (p = 0.04) in HER2+ tumors and entropy with SSF = 2, 4, 6 mm (p = 0.004, 0.02, 0.047 respectively) in luminal tumors were significantly associated with pCR. Multiple features demonstrated significant association (p < 0.05) with pCR in TNBC and with RCB in luminal tumors and TNBC, with standard deviation of intensity with SSF = 6 mm achieving the highest AUC (AUC = 0.734) for pCR in TNBC.

CONCLUSION

MRI radiomics features are associated with different molecular subtypes of breast cancer, pCR, and RCB. These features may be noninvasive imaging biomarkers to identify cancer subtype and predict response to NAC.

Breast Cancer Classification on Multiparametric MRI - Increased Performance of Boosting Ensemble Methods

Introduction: This study aims to assess the utility of Boosting ensemble classification methods for increasing the diagnostic performance of multiparametric Magnetic Resonance Imaging (mpMRI) radiomic models, in differentiating benign and malignant breast lesions. Methods: The dataset includes mpMR images of 140 female patients with mass-like breast lesions (70 benign and 70 malignant), consisting of Dynamic Contrast Enhanced (DCE) and T2-weighted sequences, and the Apparent Diffusion Coefficient (ADC) calculated from the Diffusion Weighted Imaging (DWI) sequence. Tumor masks were manually defined in all consecutive slices of the respective MRI volumes and 3D radiomic features were extracted with the Pyradiomics package. Feature dimensionality reduction was based on statistical tests and the Boruta wrapper. Hierarchical Clustering on Spearman's rank correlation coefficients between features and Random Forest classification for obtaining feature importance, were implemented for selecting the final feature subset. Adaptive Boosting (AdaBoost), Gradient Boosting (GB), Extreme Gradient Boosting (XGBoost) and Light Gradient Boosting Machine (LightGBM) classifiers, were trained and tested with bootstrap validation in differentiating breast lesions. A Support Vector Machine (SVM) classifier was also exploited for comparison. The Receiver Operator Characteristic (ROC) curves and DeLong's test were utilized to evaluate the classification performances. Results: The final feature subset consisted of 5 features derived from the lesion shape and the first order histogram of DCE and ADC images volumes. XGboost and LGBM achieved statistically significantly higher average classification performances [AUC = 0.95 and 0.94 respectively], followed by Adaboost [AUC = 0.90], GB [AUC = 0.89] and SVM [AUC = 0.88]. Conclusion: Overall, the integration of Ensemble Learning methods within mpMRI radiomic analysis can improve the performance of computer-assisted diagnosis of breast cancer lesions.

Vascularity and Dynamic Contrast-Enhanced Breast Magnetic Resonance Imaging

Angiogenesis is a key step in the initiation and progression of an invasive breast cancer. High microvessel density by morphological characterization predicts metastasis and poor survival in women with invasive breast cancers. However, morphologic characterization is subject to variability and only can evaluate a limited portion of an invasive breast cancer. Consequently, breast Magnetic Resonance Imaging (MRI) is currently being evaluated to assess vascularity. Recently, through the new field of radiomics, dynamic contrast enhanced (DCE)-MRI is being used to evaluate vascular density, vascular morphology, and detection of aggressive breast cancer biology. While DCE-MRI is a highly sensitive tool, there are specific features that limit computational evaluation of blood vessels. These include (1) DCE-MRI evaluates gadolinium contrast and does not directly evaluate biology, (2) the resolution of DCE-MRI is insufficient for imaging small blood vessels, and (3) DCE-MRI images are very difficult to co-register. Here we review computational approaches for detection and analysis of blood vessels in DCE-MRI images and present some of the strategies we have developed for co-registry of DCE-MRI images and early detection of vascularization.

Prediction of Clinically Significant Cancer Using Radiomics Features of Pre-Biopsy of Multiparametric MRI in Men Suspected of Prostate Cancer

Background: Texture features based on the spatial relationship of pixels, known as the gray-level co-occurrence matrix (GLCM), may play an important role in providing the accurate classification of suspected prostate cancer. The purpose of this study was to use quantitative imaging parameters of pre-biopsy multiparametric magnetic resonance imaging (mpMRI) for the prediction of clinically significant prostate cancer. Methods: This was a prospective study, recruiting 200 men suspected of having prostate cancer. Participants were imaged using a protocol-based 3T MRI in the pre-biopsy setting. Radiomics parameters were extracted from the T2WI and ADC texture features of the gray-level co-occurrence matrix were delineated from the region of interest. Radical prostatectomy histopathology was used as a reference standard. A Kruskal–Wallis test was applied first to identify the significant radiomic features between the three groups of Gleason scores (i.e., G1, G2 and G3). Subsequently, the Holm–Bonferroni method was applied to correct and control the probability of false rejections. We compared the probability of correctly predicting significant prostate cancer between the explanatory GLCM radiomic features, PIRADS and PSAD, using the area under the receiver operation characteristic curves. Results: We identified the significant difference in radiomic features between the three groups of Gleason scores. In total, 12 features out of 22 radiomics features correlated with the Gleason groups. Our model demonstrated excellent discriminative ability (C-statistic = 0.901, 95%CI 0.859–0.943). When comparing the probability of correctly predicting significant prostate cancer between explanatory GLCM radiomic features (Sum Variance T2WI, Sum Entropy T2WI, Difference Variance T2WI, Entropy ADC and Difference Variance ADC), PSAD and PIRADS via area under the ROC curve, radiomic features were 35.0% and 34.4% more successful than PIRADS and PSAD, respectively, in correctly predicting significant prostate cancer in our patients (p < 0.001). The Sum Entropy T2WI score had the greatest impact followed by the Sum Variance T2WI. Conclusion: Quantitative GLCM texture analyses of pre-biopsy MRI has the potential to be used as a non-invasive imaging technique to predict clinically significant cancer in men suspected of having prostate cancer.

The Potential Value of Texture Analysis Based on Dynamic Contrast-Enhanced MR Images in the Grading of Breast Phyllode Tumors

Purpose

To explore the value of texture analysis (TA) based on dynamic contrast-enhanced MR (DCE-MR) images in the differential diagnosis of benign phyllode tumors (BPTs) and borderline/malignant phyllode tumors (BMPTs).

Methods

A total of 47 patients with histologically proven phyllode tumors (PTs) from November 2012 to March 2020, including 26 benign BPTs and 21 BMPTs, were enrolled in this retrospective study. The whole-tumor texture features based on DCE-MR images were calculated, and conventional imaging findings were evaluated according to the Breast Imaging Reporting and Data System (BI-RADS). The differences in the texture features and imaging findings between BPTs and BMPTs were compared; the variates with statistical significance were entered into logistic regression analysis. The receiver operating characteristic (ROC) curve was used to assess the diagnostic performance of models from image-based analysis, TA, and the combination of these two approaches.

Results

Regarding texture features, three features of the histogram, two features of the gray-level co-occurrence matrix (GLCM), and three features of the run-length matrix (RLM) showed significant differences between the two groups (all p < 0.05). Regarding imaging findings, however, only cystic wall morphology showed significant differences between the two groups (p = 0.014). The areas under the ROC curve (AUCs) of image-based analysis, TA, and the combination of these two approaches were 0.687 (95% CI, 0.518–0.825, p = 0.014), 0.886 (95% CI, 0.760–0.960, p < 0.0001), and 0.894 (95% CI, 0.754–0.970, p < 0.0001), respectively.

Conclusion

TA based on DCE-MR images has potential in differentiating BPTs and BMPTs.

Classification of MR-Detected Additional Lesions in Patients With Breast Cancer Using a Combination of Radiomics Analysis and Machine Learning

Objective

This study was conducted in order to investigate the feasibility of using radiomics analysis (RA) with machine learning algorithms based on breast magnetic resonance (MR) images for discriminating malignant from benign MR-detected additional lesions in patients with primary breast cancer.

Materials and Methods

One hundred seventy-four MR-detected additional lesions (benign, n = 86; malignancy, n = 88) from 158 patients with ipsilateral primary breast cancer from a tertiary medical center were included in this retrospective study. The entire data were randomly split to training (80%) and independent test sets (20%). In addition, 25 patients (benign, n = 21; malignancy, n = 15) from another tertiary medical center were included for the external test. Radiomics features that were extracted from three regions-of-interest (ROIs; intratumor, peritumor, combined) using fat-saturated T1-weighted images obtained by subtracting pre- from postcontrast images (SUB) and T2-weighted image (T2) were utilized to train the support vector machine for the binary classification. A decision tree method was utilized to build a classifier model using clinical imaging interpretation (CII) features assessed by radiologists. Area under the receiver operating characteristic curve (AUROC), accuracy, sensitivity, and specificity were used to compare the diagnostic performance.

Results

The RA models trained using radiomics features from the intratumor-ROI showed comparable performance to the CII model (accuracy, AUROC: 73.3%, 69.6% for the SUB RA model; 70.0%, 75.1% for the T2 RA model; 73.3%, 72.0% for the CII model). The diagnostic performance increased when the radiomics and CII features were combined to build a fusion model. The fusion model that combines the CII features and radiomics features from multiparametric MRI data demonstrated the highest performance with an accuracy of 86.7% and an AUROC of 91.1%. The external test showed a similar pattern where the fusion models demonstrated higher levels of performance compared with the RA- or CII-only models. The accuracy and AUROC of the SUB+T2 RA+CII model in the external test were 80.6% and 91.4%, respectively.

Conclusion

Our study demonstrated the feasibility of using RA with machine learning approach based on multiparametric MRI for quantitatively characterizing MR-detected additional lesions. The fusion model demonstrated an improved diagnostic performance over the models trained with either RA or CII alone.

Mammography-based Radiomics in Breast Cancer: A Scoping Review of Current Knowledge and Future Needs

RATIONALE AND OBJECTIVES

Breast cancer is a highly complex heterogeneous disease. Current validated prognostic factors (e.g., histological grade, lymph node involvement, receptor status, and proliferation index), as well as multigene tests (e.g., Oncotype DX and PAM50) are helpful to describe breast cancer characteristics and predict the chance of recurrence risk and survival. Nevertheless, they are invasive and cannot capture a complete heterogeneity of the entire breast tumor resulting in up to 30% of patients being either over- or under-treated for breast cancer. Furthermore, multigene testings are time consuming and expensive. Radiomics is emerging as a reliable, accurate, non-invasive, and cost-effective approach of using quantitative image features to classify breast cancer characteristics and predict patient outcomes. Several recent radiomics reviews have been conducted in breast cancer, however, specific mammography-based radiomics studies have not been well discussed. This scoping review aims to assess and summarize the current evidence on the potential usefulness of mammography-based (i.e., digital mammography, digital breast tomosynthesis, and contrast-enhanced mammography) radiomics in predicting factors that describe breast cancer characteristics, recurrence, and survival.

MATERIALS AND METHODS

PubMed database and eligible text reference were searched using relevant keywords to identify studies published between 2015 and December 19, 2020. Studies collected were screened and assessed based on the inclusion and exclusion criteria.

RESULTS

Eighteen eligible studies were included and organized into three main sections: radiomics predicting breast cancer characteristics, radiomics predicting breast cancer recurrence and survival, and radiomics integrating with clinical data. Majority of publications reported retrospective studies while three studies examined prospective cohorts. Encouraging results were reported, suggesting the potential clinical value of mammography-based radiomics. Further efforts are required to standardize radiomics approaches and catalogue reproducible and relevant mammographic radiomic features. The role of integrating radiomics with other information is discussed.

CONCLUSION

The potential role of mammography-based radiomics appears promising but more efforts are required to further evaluate its reliability as a routine clinical tool.

Radiomics in breast MRI: current progress toward clinical application in the era of artificial intelligence

Breast magnetic resonance imaging (MRI) is the most sensitive imaging modality for breast cancer diagnosis and is widely used clinically. Dynamic contrast-enhanced MRI is the basis for breast MRI, but ultrafast images, T2-weighted images, and diffusion-weighted images are also taken to improve the characteristics of the lesion. Such multiparametric MRI with numerous morphological and functional data poses new challenges to radiologists, and thus, new tools for reliable, reproducible, and high-volume quantitative assessments are warranted. In this context, radiomics, which is an emerging field of research involving the conversion of digital medical images into mineable data for clinical decision-making and outcome prediction, has been gaining ground in oncology. Recent development in artificial intelligence has promoted radiomics studies in various fields including breast cancer treatment and numerous studies have been conducted. However, radiomics has shown a translational gap in clinical practice, and many issues remain to be solved. In this review, we will outline the steps of radiomics workflow and investigate clinical application of radiomics focusing on breast MRI based on published literature, as well as current discussion about limitations and challenges in radiomics.

Texture Analysis and its Applications in Biomedical Imaging: A Survey

Texture analysis describes a variety of image analysis techniques that quantify the variation in intensity and pattern. This paper provides an overview of several texture analysis approaches addressing the rationale supporting them, their advantages, drawbacks, and applications. This surveys emphasis is in collecting and categorising over five decades of active research on texture analysis. Brief descriptions of different approaches are presented along with application examples. From a broad range of texture analysis applications, this surveys final focus is on biomedical image analysis. An up-to-date list of biological tissues and organs in which disorders produce texture changes that may be used to spot disease onset and progression is provided. Finally, the role of texture analysis methods as biomarkers of disease is summarised.

Features from MRI texture analysis associated with survival outcomes in triple-negative breast cancer patients

PURPOSE

The purpose of the study is to evaluate the associations between intratumoral or peritumoral textural features derived from pretreatment magnetic resonance imaging (MRI) and recurrence-free survival (RFS) in triple-negative breast cancer (TNBC) patients.

METHODS

Forty-three patients with TNBC who underwent preoperative MRI between February 2008 and March 2014 were included. We performed two-dimensional texture analysis on the intratumoral or peritumoral region of interest (ROI) on axial of T2-weighted image (T2WI), dynamic contrast-enhanced (DCE)-MRI and DCE-MRI subtraction images. We also analyzed histopathological data. Cox proportional hazards models were used to investigate associations with survival outcomes.

RESULTS

Twelve of the 43 patients (27.9%) had recurrence disease, at a median of 32.5 months follow-up (1.4-61.5 months). In univariate analysis, nine texture features in T2WI and DCE-MRI subtraction images were significantly associated with RFS. In multivariate analysis, intratumoral difference entropy in DCE-MRI subtraction images in the initial phase (hazard ratio 11.71; 95% confidence interval (CI) [1.41, 97.00]; p value 0.023) and, peritumoral difference variance in DCE-MRI subtraction images in the delayed phase (hazard ratio 9.60; 95% CI [1.98, 46.51]; p value 0.005), were both independently associated with RFS. Moreover, multivariate analysis revealed the presence of lymphovascular invasion as independently associated with RFS (hazard ratio 8.13; 95% CI [2.16, 30.30]; p value 0.002).

CONCLUSIONS

At pretreatment MRI, an intratumoral and peritumoral quantitative approach using texture analysis has the potential to serve as a prognostic marker in patients with TNBC.

Machine learning classification of texture features of MRI breast tumor and peri-tumor of combined pre- and early treatment predicts pathologic complete response

Purpose

This study used machine learning classification of texture features from MRI of breast tumor and peri-tumor at multiple treatment time points in conjunction with molecular subtypes to predict eventual pathological complete response (PCR) to neoadjuvant chemotherapy.

Materials and method

This study employed a subset of patients (N = 166) with PCR data from the I-SPY-1 TRIAL (2002–2006). This cohort consisted of patients with stage 2 or 3 breast cancer that underwent anthracycline–cyclophosphamide and taxane treatment. Magnetic resonance imaging (MRI) was acquired pre-neoadjuvant chemotherapy, early, and mid-treatment. Texture features were extracted from post-contrast-enhanced MRI, pre- and post-contrast subtraction images, and with morphological dilation to include peri-tumoral tissue. Molecular subtypes and Ki67 were also included in the prediction model. Performance of classification models used the receiver operating characteristics curve analysis including area under the curve (AUC). Statistical analysis was done using unpaired two-tailed t-tests.

Results

Molecular subtypes alone yielded moderate prediction performance of PCR (AUC = 0.82, p = 0.07). Pre-, early, and mid-treatment data alone yielded moderate performance (AUC = 0.88, 0.72, and 0.78, p = 0.03, 0.13, 0.44, respectively). The combined pre- and early treatment data markedly improved performance (AUC = 0.96, p = 0.0003). Addition of molecular subtypes improved performance slightly for individual time points but substantially for the combined pre- and early treatment (AUC = 0.98, p = 0.0003). The optimal morphological dilation was 3–5 pixels. Subtraction of post- and pre-contrast MRI further improved performance (AUC = 0.98, p = 0.00003). Finally, among the machine-learning algorithms evaluated, the RUSBoosted Tree machine-learning method yielded the highest performance.

Conclusion

AI-classification of texture features from MRI of breast tumor at multiple treatment time points accurately predicts eventual PCR. Longitudinal changes in texture features and peri-tumoral features further improve PCR prediction performance. Accurate assessment of treatment efficacy early on could minimize unnecessary toxic chemotherapy and enable mid-treatment modification for patients to achieve better clinical outcomes.

Value of Quantitative CTTA in Differentiating Malignant From Benign Bosniak III Renal Lesions on CT Images

OBJECTIVE

The aim of this study was to investigate whether computed tomography texture analysis can differentiate malignant from benign Bosniak III renal lesions on computed tomography (CT) images.

METHODS

This retrospective case-control study included 45 patients/lesions (22 benign and 23 malignant lesions) with Bosniak III renal lesions who underwent CT examination. Axial image slices in the unenhanced phase, corticomedullary phase, and nephrographic phase were selected and delineated manually. Computed tomography texture analysis was performed on each lesion during these 3 phases. Histogram-based, gray-level co-occurrence matrix, and gray-level run-length matrix features were extracted using open-source software and analyzed. In addition, receiver operating characteristic curve was constructed, and the area under the receiver operating characteristic curve (AUC) of each feature was constructed.

RESULTS

Of the 33 extracted features, 16 features showed significant differences (P < 0.05). Eight features were significantly different between the 2 groups after Holm-Bonferroni correction, including 3 histogram-based, 4 gray-level co-occurrence matrix, and 1 gray-level run-length matrix features (P < 0.01). The texture features resulted in the highest AUC of 0.769 ± 0.074. Renal cell carcinomas were labeled with a higher degree of lesion gray-level disorder and lower lesion homogeneity, and a model incorporating the 3 most discriminative features resulted in an AUC of 0.846 ± 0.058.

CONCLUSIONS

The results of this study showed that CT texture features were related to malignancy in Bosniak III renal lesions. Computed tomography texture analysis might help in differentiating malignant from benign Bosniak III renal lesions on CT images.

Texture Analysis Using Semiquantitative Kinetic Parameter Maps from DCE-MRI: Preoperative Prediction of HER2 Status in Breast Cancer

Objective

To evaluate whether texture features derived from semiquantitative kinetic parameter maps based on breast dynamic contrast-enhanced magnetic resonance imaging (DCE-MRI) can determine human epidermal growth factor receptor 2 (HER2) status of patients with breast cancer.

Materials and Methods

This study included 102 patients with histologically confirmed breast cancer, all of whom underwent preoperative breast DCE-MRI and were enrolled retrospectively. This cohort included 48 HER2-positive cases and 54 HER2-negative cases. Seven semiquantitative kinetic parameter maps were calculated on the lesion area. A total of 55 texture features were extracted from each kinetic parameter map. Patients were randomly divided into training (n = 72) and test (n = 30) sets. The least absolute shrinkage and selection operator (LASSO) was used to select features in the training set, and then, multivariate logistic regression analysis was conducted to establish the prediction models. The classification performance was evaluated by receiver operating characteristic (ROC) analysis.

Results

Among the seven prediction models, the model with features extracted from the early signal enhancement ratio (ESER) map yielded an area under the ROC curve (AUC) of 0.83 in the training set (sensitivity of 70.59%, specificity of 92.11%, and accuracy of 81.94%), and the highest AUC of 0.83 in the test set (sensitivity of 57.14%, specificity of 100.00%, and accuracy of 80.00%). The model with features extracted from the slope of signal intensity (SI_slope) map yielded the highest AUC of 0.92 in the training set (sensitivity of 82.35%, specificity of 97.37%, and accuracy of 90.28%), and an AUC of 0.79 in the test set (sensitivity of 92.86%, specificity of 68.75%, and accuracy of 80.00%).

Conclusions

Texture features derived from kinetic parameter maps, calculated based on breast DCE-MRI, have the potential to be used as imaging biomarkers to distinguish HER2-positive and HER2-negative breast cancer.

Preoperative prediction of axillary lymph node metastasis in patients with breast cancer based on radiomics of gray-scale ultrasonography

Background

To investigate the performance of a radiomics model based on gray-scale ultrasonography (US) for the preoperative non-invasive prediction of ipsilateral axillary lymph node (ALN) metastasis in patients with breast cancer (BC).

Methods

A total of 192 pathologically confirmed BC patients were included in this study. The training set was comprised of 132 patients from hospital 1 and the test set was comprised of 60 patients from hospital 2. All patients underwent US before percutaneous core biopsy and the results of ALN status reported by a radiologist with 12 years of experience were recorded. Radiomic features were extracted from the gray-scale US images. Max-relevance and min-redundancy (MRMR) and least absolute shrinkage and selection operator (LASSO) were used for data dimension reduction and feature selection. A radiomics model was constructed using LASSO and was validated using the leave group out cross-validation (LGOCV) method. The performance of the model was validated with receiver operating characteristic (ROC), calibration curve, and decision curve analysis.

Results

A total of 860 features were extracted from the gray-scale US images of each breast lesion, and 9 radiomic features were selected for model construction. The area under the curve (AUC), sensitivity, and specificity of the model for predicting ALN metastasis were 0.85, 78.9%, and 77.3% in the training set and 0.65, 68.0%, and 79.4% in the test set, respectively. The prediction performance of the model was significantly higher than that of the radiologist (AUC: 0.85 vs. 0.59, P<0.01) in the training set and was slightly higher than that of the radiologist (AUC: 0.65 vs. 0.63, P>0.05) in the test set.

Conclusions

The non-invasive radiomics model has the ability to predict ALN metastasis for patients with breast cancer and may outperform US-reported ALN status performed by the radiologist.

Radiomics Nomograms Based on Non-enhanced MRI and Clinical Risk Factors for the Differentiation of Chondrosarcoma from Enchondroma

BACKGROUND

Differentiating chondrosarcoma from enchondroma using conventional MRI remains challenging. An effective method for accurate preoperative diagnosis could affect the management and prognosis of patients.

PURPOSE

To validate and evaluate radiomics nomograms based on non-enhanced MRI and clinical risk factors for the differentiation of chondrosarcoma from enchondroma.

STUDY TYPE

Retrospective.

POPULATION

A total of 103 patients with pathologically confirmed chondrosarcoma (n = 53) and enchondroma (n = 50) were randomly divided into training (n = 68) and validation (n = 35) groups.

FIELD STRENGTH/SEQUENCE

Axial non-contrast-enhanced T1-weighted images (T1WI) and fat-suppressed T2-weighted images (T2WI-FS) were acquired at 3.0 T.

ASSESSMENT

Clinical risk factors (sex, age, and tumor location) and diagnosis assessment based on morphologic MRI by three radiologists were recorded. Three radiomics signatures were established based on the T1WI, T2WI-FS, and T1WI + T2WI-FS sequences. Three clinical radiomics nomograms were developed based on the clinical risk factors and three radiomics signatures.

STATISTICAL TESTS

The area under the receiver operating characteristic curve (AUC) was used to evaluate the performance of radiomics signatures and clinical radiomics nomograms.

RESULTS

Tumor location was an important clinical risk factor (P < 0.05). The radiomics signature based on T1WI and T1WI + T2WI-FS features performed better than that based on T2WI-FS in the validation group (AUC in the validation group: 0.961, 0.938, and 0.833, respectively; P < 0.05). In the validation group, the three clinical radiomics nomograms (T1WI, T2WI-FS, and T1WI + T2WI-FS) achieved AUCs of 0.938, 0.935, and 0.954, respectively. In all patients, the clinical radiomics nomogram based on T2WI-FS (AUC = 0.967) performed better than that based on T2WI-FS (AUC = 0.901, P < 0.05).

DATA CONCLUSION

The proposed clinical radiomics nomogram showed promising performance in differentiating chondrosarcoma from enchondroma.

LEVEL OF EVIDENCE

4 TECHNICAL EFFICACY: Stage 2.

Radiomics based on multisequence magnetic resonance imaging for the preoperative prediction of peritoneal metastasis in ovarian cancer

OBJECTIVES

To develop a radiomics signature based on multisequence magnetic resonance imaging (MRI) to preoperatively predict peritoneal metastasis (PM) in ovarian cancer (OC).

METHODS

Eighty-nine patients with OC were divided into a training cohort including patients (n = 54) with a single lesion and a validation cohort including patients (n = 35) with bilateral lesions. Radiomics features were extracted from the T2-weighted images (T2WIs), fat-suppressed T2WIs, multi-b-value diffusion-weighted images (DWIs), and corresponding parametric maps. A radiomics signature and nomogram incorporating the radiomics signature and clinical predictors were developed and validated on the training and validation cohorts, respectively.

RESULTS

The radiomics signature generated by 6 selected features showed a favorable discriminatory ability to predict PM in OC with an area under the curve (AUC) of 0.963 in the training cohort and an AUC of 0.928 in the validation cohort. The nomogram, comprising the radiomics signature, pelvic fluid, and CA-125 level, showed more favorable discrimination with an AUC of 0.969 in the training cohort and 0.944 in the validation cohort. Net reclassification index with values of 0.548 in the training cohort and 0.500 in the validation cohort.

CONCLUSION

Radiomics signature based on multisequence MRI serves as an effective quantitative approach to predict PM in OC patients. A nomogram of radiomics signature and clinical predictors could further improve the prediction ability of PM in patients with OC.

KEY POINTS

• Multisequence MRI-based radiomics showed a favorable discriminatory ability to predict PM in OC. • The nomogram incorporating the radiomics signature and clinical predictors was clinically useful to preoperatively predict PM in patients with OC.

Muscle Composition Analysis of Ultrasound Images: A Narrative Review of Texture Analysis

Skeletal muscle composition, often characterized by the degree of intramuscular adipose tissue, deteriorates with aging and disease and contributes to impairments in function and metabolism. Ultrasound can provide surrogate measures of muscle composition through measurement of echo intensity; however, there are several limitations associated with its analysis. More complex image processing features, broadly known as texture analysis, can also provide surrogates of muscle composition and may circumvent some of the limitations associated with muscle echo intensity. Here, texture features from the intensity histogram, gray-level co-occurrence matrix, run-length matrix, local binary pattern, blob analysis, texture anisotropy index and wavelet analysis are discussed. The purpose of this review was to provide a conceptual understanding of texture analysis as it pertains to muscle composition of ultrasound images.

Value of conventional magnetic resonance imaging texture analysis in the differential diagnosis of benign and borderline/malignant phyllodes tumors of the breast

BACKGROUND

The purpose of this study was to determine the potential value of magnetic resonance imaging (MRI) texture analysis (TA) in differentiating between benign and borderline/malignant phyllodes tumors of the breast.

METHODS

The preoperative MRI data of 25 patients with benign phyllodes tumors (BPTs) and 19 patients with borderline/malignant phyllodes tumors (BMPTs) were retrospectively analyzed. A gray-level histogram and gray-level cooccurrence matrix (GLCM) were used for TA with fat-suppressed T2-weighted imaging (FS-T2WI), diffusion-weighted imaging (DWI), apparent diffusion coefficient (ADC) images, and 2- and 7-min postcontrast T1W images on dynamic contrast-enhanced MRI (DCE-T1WI_2min and DCE-T1WI_7min) between BPTs and BMPTs. Independent sample t-test and Mann-Whitney U test were performed for intergroup comparison. A regression model was established by using binary logistic regression analysis, and receiver operating characteristic (ROC) curve analysis was carried out to evaluate diagnostic efficiency.

RESULTS

For ADC images, the texture parameters angular second moment (ASM), correlation, contrast, entropy and the minimum gray values of ADC images (ADC_Minimum) showed significant differences between the BPT group and BMPT group (all p<0.05). The parameter entropy of FS-T2WI and the maximum gray values and kurtosis of the tumor solid region of DCE-T1WI_7min also showed significant differences between these two groups. Except for ADC_Minimum, angular second moment of FS-T2WI (FS-T2WI_ASM), and the maximum gray values of DCE-T1WI_7min (DCE-T1WI_7min-Maximum) of the tumor solid region, the AUC values of other positive texture parameters mentioned above were greater than 0.75. Binary logistic regression analysis demonstrated that the contrast of ADC images (ADC_Contrast) and entropy of FS-T2WI (FS-T2WI_Entropy) could be considered independent texture variables for the differential diagnosis of BPTs and BMPTs. Combined, the AUC of these parameters was 0.891 (95% CI: 0.793-0.988), with a sensitivity of 84.2% and a specificity of up to 89.0%.

CONCLUSION

Texture analysis could be helpful in improving the diagnostic efficacy of conventional MR images in differentiating BPTs and BMPTs.

Additive value of texture analysis based on breast MRI for distinguishing between benign and malignant non-mass enhancement in premenopausal women

BACKGROUND

Non-mass enhancement (NME) is a diagnostic dilemma and highly reliant on the experience of the radiologists. Texture analysis (TA) could serve as an objective method to quantify lesion characteristics. However, it remains unclear what role TA plays in a predictive model based on routine MRI characteristics. The purpose of this study was to explore the value of TA in distinguishing between benign and malignant NME in premenopausal women.

METHODS

Women in whom NME was histologically proven (n = 147) were enrolled (benign: 58; malignant: 89) was retrospective. Then, 102 and 45 patients were classified as the training and validation groups, respectively. Scanning sequences included Fat-suppressed T2-weighted and fat-suppressed contrast-enhanced T1-weighted which were acquired on a 1.5T MRI system. Clinical and routine MR characteristics (CRMC) were evaluated by two radiologists according to the Breast Imaging and Reporting and Data system (2013). Texture features were extracted from all post-contrast sequences in the training group. The combination model was built and then assessed in the validation group. Pearson's chi-square test and Mann-Whitney U test were used to compare categorical variables and continuous variables, respectively. Logistic regression analysis and receiver operating characteristic curve were employed to assess the diagnostic performance of CRMC, TA, and their combination model in NME diagnosis.

RESULTS

The combination model showed superior diagnostic performance in differentiating between benign and malignant NME compared to that of CRMC or TA alone (AUC, 0.887 vs 0.832 vs 0.74). Moreover, compared to CRMC, the model showed high specificity (72.5% vs 80%). The results obtained in the validation group confirmed the model was promising.

CONCLUSIONS

With the combined use of TA and CRMC could afford an improved diagnostic performance in differentiating between benign and malignant NME.

Functional 4-D clustering for characterizing intratumor heterogeneity in dynamic imaging: evaluation in FDG PET as a prognostic biomarker for breast cancer

Purpose

Probe-based dynamic (4-D) imaging modalities capture breast intratumor heterogeneity both spatially and kinetically. Characterizing heterogeneity through tumor sub-populations with distinct functional behavior may elucidate tumor biology to improve targeted therapy specificity and enable precision clinical decision making.

Methods

We propose an unsupervised clustering algorithm for 4-D imaging that integrates Markov-Random Field (MRF) image segmentation with time-series analysis to characterize kinetic intratumor heterogeneity. We applied this to dynamic FDG PET scans by identifying distinct time-activity curve (TAC) profiles with spatial proximity constraints. We first evaluated algorithm performance using simulated dynamic data. We then applied our algorithm to a dataset of 50 women with locally advanced breast cancer imaged by dynamic FDG PET prior to treatment and followed to monitor for disease recurrence. A functional tumor heterogeneity (FTH) signature was then extracted from functionally distinct sub-regions within each tumor. Cross-validated time-to-event analysis was performed to assess the prognostic value of FTH signatures compared to established histopathological and kinetic prognostic markers.

Results

Adding FTH signatures to a baseline model of known predictors of disease recurrence and established FDG PET uptake and kinetic markers improved the concordance statistic (C-statistic) from 0.59 to 0.74 (p = 0.005). Unsupervised hierarchical clustering of the FTH signatures identified two significant (p < 0.001) phenotypes of tumor heterogeneity corresponding to high and low FTH. Distributions of FDG flux, or Ki, were significantly different (p = 0.04) across the two phenotypes.

Conclusions

Our findings suggest that imaging markers of FTH add independent value beyond standard PET imaging metrics in predicting recurrence-free survival in breast cancer and thus merit further study.

Supplementary Information

The online version contains supplementary material available at 10.1007/s00259-021-05265-8.

Heterogeneity analysis of MRI T2 maps for measurement of early tumor response to radiotherapy

External beam radiotherapy (XRT) is a widely used cancer treatment, yet responses vary dramatically among patients. These differences are not accounted for in clinical practice, partly due to a lack of sensitive early response biomarkers. We hypothesize that quantitative magnetic resonance imaging (MRI) measures reflecting tumor heterogeneity can provide a sensitive and robust biomarker of early XRT response. MRI T2 mapping was performed every 72 hours following 10 Gy dose XRT in two models of pancreatic cancer propagated in the hind limb of mice. Interquartile range (IQR) of tumor T2 was presented as a potential biomarker of radiotherapy response compared with tumor growth kinetics, and biological validation was performed through quantitative histology analysis. Quantification of tumor T2 IQR showed sensitivity for detection of XRT-induced tumor changes 72 hours after treatment, outperforming T2-weighted and diffusion-weighted MRI, with very good robustness. Histological comparison revealed that T2 IQR provides a measure of spatial heterogeneity in tumor cell density, related to radiation-induced necrosis. Early IQR changes were found to correlate to subsequent tumor volume changes, indicating promise for treatment response prediction. Our preclinical findings indicate that spatial heterogeneity analysis of T2 MRI can provide a translatable method for early radiotherapy response assessment. We propose that the method may in future be applied for personalization of radiotherapy through adaptive treatment paradigms.

Differences in tumour heterogeneity based on dynamic contrast-enhanced MRI between tumour and peritumoural stroma for predicting Ki-67 status of invasive ductal carcinoma

AIM

To evaluate and compare the heterogeneity of intratumour and peritumour areas in the prediction of Ki-67 of invasive ductal carcinoma (IDC) and the predictive accuracy of different contrast frames based on dynamic contrast-enhanced magnetic resonance imaging (DCE-MRI).

MATERIALS AND METHODS

This study included 88 patients with histologically confirmed IDC with 57 patients with high Ki-67 status and 31 patients with low Ki-67 status. All patients underwent DCE-MRI before surgery. A grey-level co-occurrence matrix (GLCM) was performed on slice-matched images from six frames by drawing the region of the interest (ROI) on the inner and outer regions of the tumours. The correlations between texture characteristics and Ki-67 status of lesions were analysed, using the Mann-Whitney test and receiver operating characteristic curve analysis.

RESULTS

In the high-Ki-67 group, the entropy was significantly higher than that of the low-Ki-67 group (p<0.001). The entropy obtained, based on the tumour boundary as a band-like area inside and outside at the first post-contrast series, revealed the highest receiver operating characteristic (AUC = 0.765). In the multivariate analysis, a higher entropy value (>4.305; p<0.001) remained independently associated with a high-Ki-67 status after adjustment for menopausal status, tumour size, histologic grade, oestrogen receptor (ER) status, and progesterone receptor (PR) status. The other parameters did not show significant differences between the high- and low-Ki-67 groups.

CONCLUSION

Heterogeneity analysis based on DCE-MRI could discriminate between high- and low-Ki-67 status. Texture characteristics from the band-like region inside and outside the tumour boundary could predict the Ki-67 status and showed higher accuracy.

Machine Learning Based on Multi-Parametric MRI to Predict Risk of Breast Cancer

Purpose

Machine learning (ML) can extract high-throughput features of images to predict disease. This study aimed to develop nomogram of multi-parametric MRI (mpMRI) ML model to predict the risk of breast cancer.

Methods

The mpMRI included non-enhanced and enhanced T1-weighted imaging (T1WI), T2-weighted imaging (T2WI), apparent diffusion coefficient (ADC), K ^trans, K _ep, V _e, and V _p. Regions of interest were annotated in an enhanced T1WI map and mapped to other maps in every slice. 1,132 features and top-10 principal components were extracted from every parameter map. Single-parametric and multi-parametric ML models were constructed via 10 rounds of five-fold cross-validation. The model with the highest area under the curve (AUC) was considered as the optimal model and validated by calibration curve and decision curve. Nomogram was built with the optimal ML model and patients' characteristics.

Results

This study involved 144 malignant lesions and 66 benign lesions. The average age of patients with benign and malignant lesions was 42.5 years old and 50.8 years old, respectively, which were statistically different. The sixth and fourth principal components of K ^trans had more importance than others. The AUCs of K ^trans, K _ep, V _e and V _p, non-enhanced T1WI, enhanced T1WI, T2WI, and ADC models were 0.86, 0.81, 0.81, 0.83, 0.79, 0.81, 0.84, and 0.83 respectively. The model with an AUC of 0.90 was considered as the optimal model which was validated by calibration curve and decision curve. Nomogram for the prediction of breast cancer was built with the optimal ML models and patient age.

Conclusion

Nomogram could improve the ability of breast cancer prediction preoperatively.