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.

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.

BackgroundGlioblastoma (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.Methods159 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.ResultsLASSO 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).ConclusionOur 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.

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.

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.

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.

We implemented and evaluated a multiple resolution residual network (MRRN) for multiple normal organs-at-risk (OAR) segmentation from computed tomography (CT) images for thoracic radiotherapy treatment (RT) planning. Our approach simultaneously combines feature streams computed at multiple image resolutions and feature levels through residual connections. The feature streams at each level are updated as the images are passed through various feature levels. We trained our approach using 206 thoracic CT scans of lung cancer patients with 35 scans held out for validation to segment the left and right lungs, heart, esophagus, and spinal cord. This approach was tested on 60 CT scans from the open-source AAPM Thoracic Auto-Segmentation Challenge dataset. Performance was measured using the Dice Similarity Coefficient (DSC). Our approach outperformed the best-performing method in the grand challenge for hard-to-segment structures like the esophagus and achieved comparable results for all other structures. Median DSC using our method was 0.97 (interquartile range [IQR]: 0.97-0.98) for the left and right lungs, 0.93 (IQR: 0.93-0.95) for the heart, 0.78 (IQR: 0.76-0.80) for the esophagus, and 0.88 (IQR: 0.86-0.89) for the spinal cord.

analysis to predict histopathologically confirmed pseudoprogression in glioblastoma patients, Advances in Radiation Oncology (2022), doi: