BACKGROUND AND PURPOSE: Differentiating pseudoprogression, a radiation-induced treatment effect, from tumor progression on imaging is a substantial challenge in glioblastoma management. Unfortunately, guidelines set by the Response Assessment in Neuro-Oncology criteria are based solely on bidirectional diametric measurements of enhancement observed on T1WI and T2WI/FLAIR scans. We hypothesized that quantitative 3D shape features of the enhancing lesion on T1WI, and T2WI/FLAIR hyperintensities (together called the lesion habitat) can more comprehensively capture pathophysiologic differences across pseudoprogression and tumor recurrence, not appreciable on diametric measurements alone. MATERIALS AND METHODS: A total of 105 glioblastoma studies from 2 institutions were analyzed, consisting of a training (n = 59) and an independent test (n = 46) cohort. For every study, expert delineation of the lesion habitat (T1WI enhancing lesion and T2WI/FLAIR hyperintense perilesional region) was obtained, followed by extraction of 30 shape features capturing 14 “global” contour characteristics and 16 “local” curvature measures for every habitat region. Feature selection was used to identify most discriminative features on the training cohort, which were evaluated on the test cohort using a support vector machine classifier. RESULTS: The top 2 most discriminative features were identified as local features capturing total curvature of the enhancing lesion and curvedness of the T2WI/FLAIR hyperintense perilesional region. Using top features from the training cohort (training accuracy = 91.5%), we obtained an accuracy of 90.2% on the test set in distinguishing pseudoprogression from tumor progression. CONCLUSIONS: Our preliminary results suggest that 3D shape attributes from the lesion habitat can differentially express across pseudoprogression and tumor progression and could be used to distinguish these radiographically similar pathologies.

The concept of tumor field effect implies that cancer is a systemic disease with its impact way beyond the visible tumor confines. For instance, in Glioblastoma (GBM), an aggressive brain tumor, the increase in intracranial pressure due to tumor burden often leads to brain herniation and poor outcomes. Our work is based on the rationale that highly aggressive tumors tend to grow uncontrollably, leading to pronounced biomechanical tissue deformations in the normal parenchyma, which when combined with local morphological differences in the tumor confines on MRI scans, will comprehensively capture tumor field effect. Specifically, we present an integrated MRI-based descriptor, radiomic-Deformation and Textural Heterogeneity (r-DepTH). This descriptor comprises measurements of the subtle perturbations in tissue deformations throughout the surrounding normal parenchyma due to mass effect. This involves non-rigidly aligning the patients’ MRI scans to a healthy atlas via diffeomorphic registration. The resulting inverse mapping is used to obtain the deformation field magnitudes in the normal parenchyma. These measurements are then combined with a 3D texture descriptor, Co-occurrence of Local Anisotropic Gradient Orientations (COLLAGE), which captures the morphological heterogeneity and infiltration within the tumor confines, on MRI scans. In this work, we extensively evaluated r-DepTH for survival risk-stratification on a total of 207 GBM cases from 3 different cohorts (Cohort 1 ( ${n}{1} = {53}$ ), Cohort 2 ( ${n}{2} = {75}$ ), and Cohort 3 ( ${n}{3} = {79}$ )), where each of these three cohorts was used as a training set for our model separately, and the other two cohorts were used for testing, independently, for each training experiment. When employing Cohort 1 for training, r-DepTH yielded Concordance indices (C-indices) of 0.7 and 0.65, hazard ratios (HR) and Confidence Intervals (CI) of 10 (6 − 19) and 5 (3 − 8) on Cohorts 2 and 3, respectively. Similarly, training on Cohort 2 yielded C-indices of 0.6 and 0.7, HR and CI of 1 (0.7 − 2) and 3 (2 − 5) on Cohorts 1 and 3, respectively. Finally, training on Cohort 3 yielded C-indices of 0.75 and 0.63, HR and CI of 24 (10 − 57) and 12 (6 − 21) on Cohorts 1 and 2, respectively. Our results show that r-DepTH descriptor may serve as a comprehensive and a robust MRI-based prognostic marker of disease aggressiveness and survival in solid tumors.

To evaluate nnU-Net–based segmentation models for automated delineation of medulloblastoma tumors on multi-institutional MRI scans. This retrospective study included 78 pediatric patients (52 male, 26 female), with ages ranging from 2 to 18 years, with medulloblastomas, from three different sites (28 from hospital A, 18 from hospital B, and 32 from hospital C), who had data available from three clinical MRI protocols (gadolinium-enhanced T1-weighted, T2-weighted, and fluid-attenuated inversion recovery). The scans were retrospectively collected from the year 2000 until May 2019. Reference standard annotations of the tumor habitat, including enhancing tumor, edema, and cystic core plus nonenhancing tumor subcompartments, were performed by two experienced neuroradiologists. Preprocessing included registration to age-appropriate atlases, skull stripping, bias correction, and intensity matching. The two models were trained as follows: (a) the transfer learning nnU-Net model was pretrained on an adult glioma cohort (n = 484) and fine-tuned on medulloblastoma studies using Models Genesis and (b) the direct deep learning nnU-Net model was trained directly on the medulloblastoma datasets, across fivefold cross-validation. Model robustness was evaluated on the three datasets when using different combinations of training and test sets, with data from two sites at a time used for training and data from the third site used for testing. Analysis on the three test sites yielded Dice scores of 0.81, 0.86, and 0.86 and 0.80, 0.86, and 0.85 for tumor habitat; 0.68, 0.84, and 0.77 and 0.67, 0.83, and 0.76 for enhancing tumor; 0.56, 0.71, and 0.69 and 0.56, 0.71, and 0.70 for edema; and 0.32, 0.48, and 0.43 and 0.29, 0.44, and 0.41 for cystic core plus nonenhancing tumor for the transfer learning and direct nnU-Net models, respectively. The models were largely robust to site-specific variations. nnU-Net segmentation models hold promise for accurate, robust automated delineation of medulloblastoma tumor subcompartments, potentially leading to more effective radiation therapy planning in pediatric medulloblastoma. Keywords: Pediatrics, MR Imaging, Segmentation, Transfer Learning, Medulloblastoma, nnU-Net, MRI Supplemental material is available for this article. © RSNA, 2024 See also the commentary by Rudie and Correia de Verdier in this issue.

Localized disease heterogeneity on imaging extracted via radiomics approaches have recently been associated with disease prognosis and treatment response. Traditionally, radiomics analyses leverage texture operators to derive voxel- or region-wise feature values towards quantifying subtle variations in image appearance within a region-of-interest (ROI). With the goal of mining additional voxel-wise texture patterns from radiomic expression maps, we introduce a new RADIomic Spatial TexturAl descripTor (RADISTAT). This was driven by the hypothesis that quantifying spatial organization of texture patterns within an ROI could allow for better capturing interactions between different tissue classes present in a given region; thus enabling more accurate characterization of disease or response phenotypes. RADISTAT involves: (a) robustly identifying sub-compartments of low, intermediate, and high radiomic expression (i.e. heterogeneity) in a feature map and (b) quantifying spatial organization of sub-compartments via graph interactions. RADISTAT was evaluated in two clinically challenging problems: (1) discriminating nodal/distant metastasis from metastasis-free rectal cancer patients on post-chemoradiation T2w MRI, and (2) distinguishing tumor progression from pseudo-progression in glioblastoma multiforme using post-chemoradiation T1w MRI. Across over 800 experiments, RADISTAT yielded a consistent discriminatory signature for tumor progression (GBM) and disease metastasis (RCa); where its sub-compartments were associated with pathologic tissue types (fibrosis or tumor, determined via fusion of MRI and pathology). In a multi-institutional setting for both clinical problems, RADISTAT resulted in higher classifier performance (11% improvement in AUC, on average) compared to radiomic descriptors. Furthermore, combining RADISTAT with radiomic descriptors resulted in significantly improved performance compared to using radiomic descriptors alone.

Recent advances in artificial intelligence have greatly impacted the field of medical imaging and vastly improved the development of computational algorithms for data analysis. In the field of pediatric neuro-oncology, radiomics, the process of obtaining high-dimensional data from radiographic images, has been recently utilized in applications including survival prognostication, molecular classification, and tumor type classification. Similarly, radiogenomics, or the integration of radiomic and genomic data, has allowed for building comprehensive computational models to better understand disease etiology. While there exist excellent review articles on radiomics and radiogenomic pipelines and their applications in adult solid tumors, in this review article, we specifically review these computational approaches in the context of pediatric medulloblastoma tumors. Based on our systematic literature research via PubMed and Google Scholar, we provide a detailed summary of a total of 15 articles that have utilized radiomic and radiogenomic analysis for survival prognostication, tumor segmentation, and molecular subgroup classification in the context of pediatric medulloblastoma. Lastly, we shed light on the current challenges with the existing approaches as well as future directions and opportunities with using these computational radiomic and radiogenomic approaches for pediatric medulloblastoma tumors.

Introduction: Medulloblastoma (MB) is a malignant, heterogenous brain tumor. Advances in molecular profiling have led to identifying four molecular subgroups of MB (WNT, SHH, Group 3, Group 4), each with distinct clinical behaviors. We hypothesize that (1) aggressive MB tumors, growing heterogeneously, induce pronounced local structural deformations in the surrounding parenchyma, and (b) these local deformations as captured on Gadolinium (Gd)-enhanced-T1w MRI are independently associated with molecular subgroups, as well as overall survival in MB patients.Methods: In this work, a total of 88 MB studies from 2 institutions were analyzed. Following tumor delineation, Gd-T_1w scan for every patient was registered to a normal age-specific T_1w-MRI template via deformable registration. Following patient-atlas registration, local structural deformations in the brain parenchyma were obtained for every patient by computing statistics from deformation magnitudes obtained from every 5mm annular region, 0 < d < 60 mm, where d is the distance from the tumor infiltrating edge.Results: Multi-class comparison via ANOVA yielded significant differences between deformation magnitudes obtained for Group 3, Group 4, and SHH molecular subgroups, observed up to 60-mm outside the tumor edge. Additionally, Kaplan-Meier survival analysis showed that the local deformation statistics, combined with the current clinical risk-stratification approaches (molecular subgroup information and Chang’s classification), could identify significant differences between high-risk and low-risk survival groups, achieving better performance results than using any of these approaches individually.Discussion: These preliminary findings suggest there exists significant association of our tumor-induced deformation descriptor with overall survival in MB, and that there could be an added value in using the proposed radiomic descriptor along with the current risk classification approaches, towards more reliable risk assessment in pediatric MB.