A Radiomic Approach for Evaluating Intra-Subgroup Heterogeneity in SHH and Group 4 Pediatric Medulloblastoma: A Preliminary Multi-Institutional Study.

Medulloblastoma (MB) is the most frequent malignant brain tumor in children with extensive heterogeneity that results in varied clinical outcomes. Recently, MB was categorized into four molecular subgroups, WNT, SHH, Group 3, and Group 4. While SHH and Group 4 are known for their intermediate prognosis, studies have reported wide disparities in patient outcomes within these subgroups. This study aims to create a radiomic prognostic signature, medulloblastoma radiomics risk (mRRisk), to identify the risk levels within the SHH and Group 4 subgroups, individually, for reliable risk stratification. Our hypothesis is that this signature can comprehensively capture tumor characteristics that enable the accurate identification of the risk level. In total, 70 MB studies (48 Group 4, and 22 SHH) were retrospectively curated from three institutions. For each subgroup, 232 hand-crafted features that capture the entropy, surface changes, and contour characteristics of the tumor were extracted. Features were concatenated and fed into regression models for risk stratification. Contrasted with Chang stratification that did not yield any significant differences within subgroups, significant differences were observed between two risk groups in Group 4 (p = 0.04, Concordance Index (CI) = 0.82) on the cystic core and non-enhancing tumor, and SHH (p = 0.03, CI = 0.74) on the enhancing tumor. Our results indicate that radiomics may serve as a prognostic tool for refining MB risk stratification, towards improved patient care.

Leveraging radiomics and machine learning to differentiate radiation necrosis from recurrence in patients with brain metastases.

OBJECTIVE: Radiation necrosis (RN) can be difficult to radiographically discern from tumor progression after stereotactic radiosurgery (SRS). The objective of this study was to investigate the utility of radiomics and machine learning (ML) to differentiate RN from recurrence in patients with brain metastases treated with SRS. METHODS: Patients with brain metastases treated with SRS who developed either RN or tumor reccurence were retrospectively identified. Image preprocessing and radiomic feature extraction were performed using ANTsPy and PyRadiomics, yielding 105 features from MRI T1-weighted post-contrast (T1c), T2, and fluid-attenuated inversion recovery (FLAIR) images. Univariate analysis assessed significance of individual features. Multivariable analysis employed various classifiers on features identified as most discriminative through feature selection. ML models were evaluated through cross-validation, selecting the best model based on area under the receiver operating characteristic (ROC) curve (AUC). Specificity, sensitivity, and F1 score were computed. RESULTS: Sixty-six lesions from 55 patients were identified. On univariate analysis, 27 features from the T1c sequence were statistically significant, while no features were significant from the T2 or FLAIR sequences. For clinical variables, only immunotherapy use after SRS was significant. Multivariable analysis of features from the T1c sequence yielded an AUC of 76.2% (standard deviation [SD] ± 12.7%), with specificity and sensitivity of 75.5% (± 13.4%) and 62.3% (± 19.6%) in differentiating radionecrosis from recurrence. CONCLUSIONS: Radiomics with ML may assist the diagnostic ability of distinguishing RN from tumor recurrence after SRS. Further work is needed to validate this in a larger multi-institutional cohort and prospectively evaluate it's utility in patient care.

Radiomic Analysis to Predict Histopathologically Confirmed Pseudoprogression in Glioblastoma Patients.

Purpose: Pseudoprogression mimicking recurrent glioblastoma remains a diagnostic challenge that may adversely confound or delay appropriate treatment or clinical trial enrollment. We sought to build a radiomic classifier to predict pseudoprogression in patients with primary isocitrate dehydrogenase wild type glioblastoma. Methods and Materials: We retrospectively examined a training cohort of 74 patients with isocitrate dehydrogenase wild type glioblastomas with brain magnetic resonance imaging including dynamic contrast enhanced T1 perfusion before resection of an enhancing lesion indeterminate for recurrent tumor or pseudoprogression. A recursive feature elimination random forest classifier was built using nested cross-validation without and with O6-methylguanine-DNA methyltransferase status to predict pseudoprogression. Results: A classifier constructed with cross-validation on the training cohort achieved an area under the receiver operating curve of 81% for predicting pseudoprogression. This was further improved to 89% with the addition of O6-methylguanine-DNA methyltransferase status into the classifier. Conclusions: Our results suggest that radiomic analysis of contrast T1-weighted images and magnetic resonance imaging perfusion images can assist the prompt diagnosis of pseudoprogression. Validation on external and independent data sets is necessary to verify these advanced analyses, which can be performed on routinely acquired clinical images and may help inform clinical treatment decisions.

Volumetric analysis of IDH-mutant lower-grade glioma: a natural history study of tumor growth rates before and after treatment.

BACKGROUND: Lower-grade gliomas (LGGs) with isocitrate dehydrogenase 1 and/or 2 (IDH1/2) mutations have long survival times, making evaluation of treatment efficacy difficult. We investigated the volumetric growth rate of IDH mutant gliomas before and after treatment with established glioma therapies to determine whether a significant change in growth rate could be documented and perhaps be used in the future to evaluate treatment response to investigational agents in LGG trials. METHODS: In this multicenter retrospective study, 230 adult patients with IDH1/2 mutated LGGs (World Health Organization grade II or III) undergoing surgery, radiation, or chemotherapy for progressive non-enhancing tumor were identified. Subjects were required to have 3 MRI scans containing T2/fluid attenuated inversion recovery imaging spanning a minimum of 6 months prior to treatment. A mixed-effect model was used to estimate tumor growth prior to treatment. A subset of 95 patients who received chemotherapy, radiotherapy, or chemoradiotherapy and had 2 posttreatment imaging time points available were evaluated for change in pre- and posttreatment volumetric growth rates using a piecewise mixed model. RESULTS: The pretreatment volumetric growth rate across all 230 patients was 27.37%/180 days (95% CI: [23.36%, 31.51%]). In the 95 patients with both pre- and posttreatment scans available, there was a significant difference in volumetric growth rates before (26.63%/180 days, 95% CI: [19.31%, 34.40%]) and after treatment (-15.24% /180 days, 95% CI: [-21.37%, -8.62%]) (P < 0.0001). The growth rates for patient subgroup with 1p/19q codeletion (N = 118) was significantly slower than the rate of the 1p/19q non-codeleted group (N = 68) (22.84% vs 35.49%, P = 0.0108). CONCLUSION: In this study, we evaluated the growth rates of IDH mutant gliomas before and after standard therapy. Further study is needed to establish whether a change in growth rate is associated with patient survival and its use as a surrogate endpoint in clinical trials for IDH mutant LGGs.

MRI radiomic features are associated with survival in melanoma brain metastases treated with immune checkpoint inhibitors.

BACKGROUND: Melanoma brain metastases historically portend a dismal prognosis, but recent advances in immune checkpoint inhibitors (ICIs) have been associated with durable responses in some patients. There are no validated imaging biomarkers associated with outcomes in patients with melanoma brain metastases receiving ICIs. We hypothesized that radiomic analysis of magnetic resonance images (MRIs) could identify higher-order features associated with survival. METHODS: Between 2010 and 2019, we retrospectively reviewed patients with melanoma brain metastases who received ICI. After volumes of interest were drawn, several texture and edge descriptors, including first-order, Haralick, Gabor, Sobel, and Laplacian of Gaussian (LoG) features were extracted. Progression was determined using Response Assessment in Neuro-Oncology Brain Metastases. Univariate Cox regression was performed for each radiomic feature with adjustment for multiple comparisons followed by Lasso regression and multivariate analysis. RESULTS: Eighty-eight patients with 196 total brain metastases were identified. Median age was 63.5 years (range, 19-91 y). Ninety percent of patients had Eastern Cooperative Oncology Group performance status of 0 or 1 and 35% had elevated lactate dehydrogenase. Sixty-three patients (72%) received ipilimumab, 11 patients (13%) received programmed cell death protein 1 blockade, and 14 patients (16%) received nivolumab plus ipilimumab. Multiple features were associated with increased overall survival (OS), and LoG edge features best explained the variation in outcome (hazard ratio: 0.68, P = 0.001). In multivariate analysis, a similar trend with LoG was seen, but no longer significant with OS. Findings were confirmed in an independent cohort. CONCLUSION: Higher-order MRI radiomic features in patients with melanoma brain metastases receiving ICI were associated with a trend toward improved OS.

Impact of image preprocessing on the scanner dependence of multi-parametric MRI radiomic features and covariate shift in multi-institutional glioblastoma datasets.

Recent advances in radiomics have enhanced the value of medical imaging in various aspects of clinical practice, but a crucial component that remains to be investigated further is the robustness of quantitative features to imaging variations and across multiple institutions. In the case of MRI, signal intensity values vary according to the acquisition parameters used, yet no consensus exists on which preprocessing techniques are favorable in reducing scanner-dependent variability of image-based features. Hence, the purpose of this study was to assess the impact of common image preprocessing methods on the scanner dependence of MRI radiomic features in multi-institutional glioblastoma multiforme (GBM) datasets. Two independent GBM cohorts were analyzed: 50 cases from the TCGA-GBM dataset and 111 cases acquired in our institution, and each case consisted of 3 MRI sequences viz. FLAIR, T1-weighted, and T1-weighted post-contrast. Five image preprocessing techniques were examined: 8-bit global rescaling, 8-bit local rescaling, bias field correction, histogram standardization, and isotropic resampling. A total of 420 features divided into eight categories representing texture, shape, edge, and intensity histogram were extracted. Two distinct imaging parameters were considered: scanner manufacturer and scanner magnetic field strength. Wilcoxon tests identified features robust to the considered acquisition parameters under the selected image preprocessing techniques. A machine learning-based strategy was implemented to measure the covariate shift between the analyzed datasets using features computed using the aforementioned preprocessing methods. Finally, radiomic scores (rad-scores) were constructed by identifying features relevant to patients' overall survival after eliminating those impacted by scanner variability. These were then evaluated for their prognostic significance through Kaplan-Meier and Cox hazards regression analyses. Our results demonstrate that overall, histogram standardization contributes the most in reducing radiomic feature variability as it is the technique to reduce the covariate shift for three feature categories and successfully discriminate patients into groups of different survival risks.

Reliability of tumor segmentation in glioblastoma: Impact on the robustness of MRI-radiomic features.

PURPOSE: The use of radiomic features as biomarkers of treatment response and outcome or as correlates to genomic variations requires that the computed features are robust and reproducible. Segmentation, a crucial step in radiomic analysis, is a major source of variability in the computed radiomic features. Therefore, we studied the impact of tumor segmentation variability on the robustness of MRI radiomic features. METHOD: Fluid-attenuated inversion recovery (FLAIR) and contrast-enhanced T1-weighted (T1WICE ) MRI of 90 patients diagnosed with glioblastoma were segmented using a semiautomatic algorithm and an interactive segmentation with two different raters. We analyzed the robustness of 108 radiomic features from five categories (intensity histogram, gray-level co-occurrence matrix, gray-level size-zone matrix (GLSZM), edge maps, and shape) using intra-class correlation coefficient (ICC) and Bland and Altman analysis. RESULTS: Our results show that both segmentation methods are reliable with ICC ≥ 0.96 and standard deviation (SD) of mean differences between the two raters (SDdiffs ) ≤ 30%. Features computed from the histogram and co-occurrence matrices were found to be the most robust (ICC ≥ 0.8 and SDdiffs  ≤ 30% for most features in these groups). Features from GLSZM were shown to have mixed robustness. Edge, shape, and GLSZM features were the most impacted by the choice of segmentation method with the interactive method resulting in more robust features than the semiautomatic method. Finally, features computed from T1WICE and FLAIR images were found to have similar robustness when computed with the interactive segmentation method. CONCLUSION: Semiautomatic and interactive segmentation methods using two raters are both reliable. The interactive method produced more robust features than the semiautomatic method. We also found that the robustness of radiomic features varied by categories. Therefore, this study could help motivate segmentation methods and feature selection in MRI radiomic studies.

Preoperative MRI-radiomics features improve prediction of survival in glioblastoma patients over MGMT methylation status alone.

BACKGROUND: Glioblastoma (GBM) is the most common malignant central nervous system tumor, and MGMT promoter hypermethylation in this tumor has been shown to be associated with better prognosis. We evaluated the capacity of radiomics features to add complementary information to MGMT status, to improve the ability to predict prognosis. METHODS: 159 patients with untreated GBM were included in this study and divided into training and independent test sets. 286 radiomics features were extracted from the magnetic resonance images acquired prior to any treatments. A least absolute shrinkage selection operator (LASSO) selection followed by Kaplan-Meier analysis was used to determine the prognostic value of radiomics features to predict overall survival (OS). The combination of MGMT status with radiomics was also investigated and all results were validated on the independent test set. RESULTS: LASSO analysis identified 8 out of the 286 radiomic features to be relevant which were then used for determining association to OS. One feature (edge descriptor) remained significant on the external validation cohort after multiple testing (p=0.04) and the combination with MGMT identified a group of patients with the best prognosis with a survival probability of 0.61 after 43 months (p=0.0005). CONCLUSION: Our results suggest that combining radiomics with MGMT is more accurate in stratifying patients into groups of different survival risks when compared to with using these predictors in isolation. We identified two subgroups within patients who have methylated MGMT: one with a similar survival to unmethylated MGMT patients and the other with a significantly longer OS.

Technical Note: Extension of CERR for computational radiomics: A comprehensive MATLAB platform for reproducible radiomics research.

PURPOSE: Radiomics is a growing field of image quantitation, but it lacks stable and high-quality software systems. We extended the capabilities of the Computational Environment for Radiological Research (CERR) to create a comprehensive, open-source, MATLAB-based software platform with an emphasis on reproducibility, speed, and clinical integration of radiomics research. METHOD: The radiomics tools in CERR were designed specifically to quantitate medical images in combination with CERR's core functionalities of radiological data import, transformation, management, image segmentation, and visualization. CERR allows for batch calculation and visualization of radiomics features, and provides a user-friendly data structure for radiomics metadata. All radiomics computations are vectorized for speed. Additionally, a test suite is provided for reconstruction and comparison with radiomics features computed using other software platforms such as the Insight Toolkit (ITK) and PyRadiomics. CERR was evaluated according to the standards defined by the Image Biomarker Standardization Initiative. CERR's radiomics feature calculation was integrated with the clinically used MIM software using its MATLAB® application programming interface. RESULTS: The CERR provides a comprehensive computational platform for radiomics analysis. Matrix formulations for the compute-intensive Haralick texture resulted in speeds that are superior to the implementation in ITK 4.12. For an image discretized into 32 bins, CERR achieved a speedup of 3.5 times over ITK. The CERR test suite enabled the successful identification of programming errors as well as genuine differences in radiomics definitions and calculations across the software packages tested. CONCLUSION: The CERR's radiomics capabilities are comprehensive, open-source, and fast, making it an attractive platform for developing and exploring radiomics signatures across institutions. The ability to both choose from a wide variety of radiomics implementations and to integrate with a clinical workflow makes CERR useful for retrospective as well as prospective research analyses.