Multiple-response regression analysis links magnetic resonance imaging features to de-regulated protein expression and pathway activity in lower grade glioma

[Radioproteomics for Discriminating the Activity and Inactivity of Immune Checkpoint Molecules in Breast Cancer]

PURPOSE

Radioproteomics studies investigating the relationship between lesion phenotype and proteins have been progressed. The purpose of this study was to develop a radioproteomics method for discriminating between active and inactive immune checkpoint molecules based on lesion phenotype.

METHODS

From the public database TCGA-BRCA, mRNA and fat suppression contrast-enhanced T1-weighted images of 49 patients with breast cancer were selected for the experiment. Using mRNA, we defined cases with active (10 cases) and inactive (39 cases) immune checkpoint molecules. To discriminate these cases using lesion phenotype, 275 radiomics features were measured from the tumor area. After selecting 3 radiomics features by using Lasso, logistic regression was employed to discriminate between active and inactive cases of immune checkpoint molecules.

RESULTS

Evaluation of ROC analysis showed that the AUC was 0.81.

CONCLUSION

Patients whose immune cell function is being braked by immune checkpoint molecules are likely to respond to immune checkpoint inhibitors when their activity is inhibited. Therefore, our results may be applied to predict the effects of immune checkpoint inhibitors in breast cancer treatment.

Multi-Parametric MRI-Based Radiomics Models for Predicting Molecular Subtype and Androgen Receptor Expression in Breast Cancer

Objective

To investigate whether radiomics features extracted from multi-parametric MRI combining machine learning approach can predict molecular subtype and androgen receptor (AR) expression of breast cancer in a non-invasive way.

Materials and Methods

Patients diagnosed with clinical T2–4 stage breast cancer from March 2016 to July 2020 were retrospectively enrolled. The molecular subtypes and AR expression in pre-treatment biopsy specimens were assessed. A total of 4,198 radiomics features were extracted from the pre-biopsy multi-parametric MRI (including dynamic contrast-enhancement T1-weighted images, fat-suppressed T2-weighted images, and apparent diffusion coefficient map) of each patient. We applied several feature selection strategies including the least absolute shrinkage and selection operator (LASSO), and recursive feature elimination (RFE), the maximum relevance minimum redundancy (mRMR), Boruta and Pearson correlation analysis, to select the most optimal features. We then built 120 diagnostic models using distinct classification algorithms and feature sets divided by MRI sequences and selection strategies to predict molecular subtype and AR expression of breast cancer in the testing dataset of leave-one-out cross-validation (LOOCV). The performances of binary classification models were assessed via the area under the receiver operating characteristic curve (AUC), accuracy, sensitivity, specificity, positive predictive value (PPV), and negative predictive value (NPV). And the performances of multiclass classification models were assessed via AUC, overall accuracy, precision, recall rate, and F1-score.

Results

A total of 162 patients (mean age, 46.91 ± 10.08 years) were enrolled in this study; 30 were low-AR expression and 132 were high-AR expression. HR+/HER2− cancers were diagnosed in 56 cases (34.6%), HER2+ cancers in 81 cases (50.0%), and TNBC in 25 patients (15.4%). There was no significant difference in clinicopathologic characteristics between low-AR and high-AR groups (P > 0.05), except the menopausal status, ER, PR, HER2, and Ki-67 index (P = 0.043, <0.001, <0.001, 0.015, and 0.006, respectively). No significant difference in clinicopathologic characteristics was observed among three molecular subtypes except the AR status and Ki-67 (P = <0.001 and 0.012, respectively). The Multilayer Perceptron (MLP) showed the best performance in discriminating AR expression, with an AUC of 0.907 and an accuracy of 85.8% in the testing dataset. The highest performances were obtained for discriminating TNBC vs. non-TNBC (AUC: 0.965, accuracy: 92.6%), HER2+ vs. HER2− (AUC: 0.840, accuracy: 79.0%), and HR+/HER2− vs. others (AUC: 0.860, accuracy: 82.1%) using MLP as well. The micro-AUC of MLP multiclass classification model was 0.896, and the overall accuracy was 0.735.

Conclusions

Multi-parametric MRI-based radiomics combining with machine learning approaches provide a promising method to predict the molecular subtype and AR expression of breast cancer non-invasively.

Exploratory Analysis of Qualitative MR Imaging Features for the Differentiation of Glioblastoma and Brain Metastases

OBJECTIVES

To identify qualitative VASARI (Visually AcceSIble Rembrandt Images) Magnetic Resonance (MR) Imaging features for differentiation of glioblastoma (GBM) and brain metastasis (BM) of different primary tumors.

MATERIALS AND METHODS

T1-weighted pre- and post-contrast, T2-weighted, and T2-weighted, fluid attenuated inversion recovery (FLAIR) MR images of a total of 239 lesions from 109 patients with either GBM or BM (breast cancer, non-small cell (NSCLC) adenocarcinoma, NSCLC squamous cell carcinoma, small-cell lung cancer (SCLC)) were included. A set of adapted, qualitative VASARI MR features describing tumor appearance and location was scored (binary; 1 = presence of feature, 0 = absence of feature). Exploratory data analysis was performed on binary scores using a combination of descriptive statistics (proportions with 95% binomial confidence intervals), unsupervised methods and supervised methods including multivariate feature ranking using either repeated fitting or recursive feature elimination with Support Vector Machines (SVMs).

RESULTS

GBMs were found to involve all lobes of the cerebrum with a fronto-occipital gradient, often affected the corpus callosum (32.4%, 95% CI 19.1-49.2), and showed a strong preference for the right hemisphere (79.4%, 95% CI 63.2-89.7). BMs occurred most frequently in the frontal lobe (35.1%, 95% CI 28.9-41.9) and cerebellum (28.3%, 95% CI 22.6-34.8). The appearance of GBMs was characterized by preference for well-defined non-enhancing tumor margin (100%, 89.8-100), ependymal extension (52.9%, 36.7-68.5) and substantially less enhancing foci than BMs (44.1%, 28.9-60.6 vs. 75.1%, 68.8-80.5). Unsupervised and supervised analyses showed that GBMs are distinctively different from BMs and that this difference is driven by definition of non-enhancing tumor margin, ependymal extension and features describing laterality. Differentiation of histological subtypes of BMs was driven by the presence of well-defined enhancing and non-enhancing tumor margins and localization in the vision center. SVM models with optimal hyperparameters led to weighted F1-score of 0.865 for differentiation of GBMs from BMs and weighted F1-score of 0.326 for differentiation of BM subtypes.

CONCLUSION

VASARI MR imaging features related to definition of non-enhancing margin, ependymal extension, and tumor localization may serve as potential imaging biomarkers to differentiate GBMs from BMs.

TCGA-TCIA Impact on Radiogenomics Cancer Research: A Systematic Review

In the last decade, the development of radiogenomics research has produced a significant amount of papers describing relations between imaging features and several molecular 'omic signatures arising from next-generation sequencing technology and their potential role in the integrated diagnostic field. The most vulnerable point of many of these studies lies in the poor number of involved patients. In this scenario, a leading role is played by The Cancer Genome Atlas (TCGA) and The Cancer Imaging Archive (TCIA), which make available, respectively, molecular 'omic data and linked imaging data. In this review, we systematically collected and analyzed radiogenomic studies based on TCGA-TCIA data. We organized literature per tumor type and molecular 'omic data in order to discuss salient imaging genomic associations and limitations of each study. Finally, we outlined the potential clinical impact of radiogenomics to improve the accuracy of diagnosis and the prediction of patient outcomes in oncology.

Radiogenomics: bridging imaging and genomics

From diagnostics to prognosis to response prediction, new applications for radiomics are rapidly being developed. One of the fastest evolving branches involves linking imaging phenotypes to the tumor genetic profile, a field commonly referred to as "radiogenomics." In this review, a general outline of radiogenomic literature concerning prominent mutations across different tumor sites will be provided. The field of radiogenomics originates from image processing techniques developed decades ago; however, many technical and clinical challenges still need to be addressed. Nevertheless, increasingly accurate and robust radiogenomic models are being presented and the future appears to be bright.

Use of Genome-Scale Integrated Analysis to Identify Key Genes and Potential Molecular Mechanisms in Recurrence of Lower-Grade Brain Glioma

BACKGROUND The aim of this study was to identify gene signals for lower-grade glioma (LGG) and to assess their potential as recurrence biomarkers. MATERIAL AND METHODS An LGG-related mRNA sequencing dataset was downloaded from The Cancer Genome Atlas (TCGA) Informix. Multiple bioinformatics analysis methods were used to identify key genes and potential molecular mechanisms in recurrence of LGG. RESULTS A total of 326 differentially-expressed genes (DEGs), were identified from 511 primary LGG tumor and 18 recurrent samples. Gene ontology (GO) analysis revealed that the DEGs were implicated in cell differentiation, neuron differentiation, negative regulation of neuron differentiation, and cell proliferation in the forebrain. The Kyoto Encyclopedia of Genes and Genomes (KEGG) database suggests that DEGs are associated with proteoglycans in cancer, the Wnt signaling pathway, ECM-receptor interaction, the PI3K-Akt signaling pathway, transcriptional deregulation in cancer, and the Hippo signaling pathway. The hub DEGs in the protein-protein interaction network are apolipoprotein A2 (APOA2), collagen type III alpha 1 chain (COL3A1), collagen type I alpha 1 chain (COL1A1), tyrosinase (TYR), collagen type I alpha 2 chain (COL1A2), neurotensin (NTS), collagen type V alpha 1 chain (COL5A1), poly(A) polymerase beta (PAPOLB), insulin-like growth factor 2 mRNA-binding protein 1 (IGF2BP1), and anomalous homeobox (ANHX). GSEA revealed that the following biological processes may associated with LGG recurrence: cell cycle, DNA replication and repair, regulation of apoptosis, neuronal differentiation, and Wnt signaling pathway. CONCLUSIONS Our study demonstrated that hub DEGs may assist in the molecular understanding of LGG recurrence. These findings still need further molecular studies to identify the assignment of DEGs in LGG.