Identification of a single-dose, low flip angle based CBV threshold for fractional tumor burden (FTB) mapping in recurrent glioblastoma

BACKGROUND AND PURPOSE

DSC-MRI can be used to generate fractional tumor burden (FTB) maps, via application of relative CBV thresholds, to spatially differentiate glioblastoma recurrence from post treatment radiation effects (PTRE). Image-localized histopathology was previously used to validate FTB maps derived from a reference DSC-MRI protocol using preload, a moderate flip angle (MFA, 60°) and post-processing leakage correction. Recently, a DSC-MRI protocol with a low flip angle (LFA, 30°) with no preload was shown to provide leakage-corrected RCBV equivalent to the reference protocol. This study aims to identify the RCBV thresholds for the LFA protocol that generate the most accurate FTB maps, concordant with those obtained from the reference MFA protocol.

MATERIALS AND METHODS

Fifty-two patients with grade IV GBM who had prior surgical resection and received chemotherapy and radiotherapy were included in the study. Two sets of DSC-MRI data were collected sequentially first using LFA protocol with no preload, which served as the preload for the subsequent MFA protocol. Standardized relative CBV maps (sRCBV) were obtained for each patient and co-registered with the anatomical post-contrast T1-weighted images. The reference MFA-based FTB maps were computed using previously published sRCBV thresholds (1.0 and 1.56). An ROC analysis was conducted to identify the optimal, voxelwise LFA sRCBV thresholds, and the sensitivity, specificity, and accuracy of the LFA-based FTB maps were computed with respect to the MFA-based reference.

RESULTS

The mean sRCBV values of tumors across patients exhibited strong agreement (CCC = 0.99) between the two protocols. Using the ROC analysis, the optimal lower LFA threshold that accurately distinguishes PTRE from tumor recurrence was found to be 1.0 (sensitivity: 87.77%; specificity: 90.22%), equivalent to the ground truth. To identify aggressive tumor regions, the ROC analysis identified an upper LFA threshold of 1.37 (sensitivity: 90.87%; specificity: 91.10%) for the reference MFA threshold of 1.56.

CONCLUSION

For LFA-based FTB maps, a sRCBV threshold of 1.0 and 1.37 can differentiate PTRE from recurrent tumor. FTB maps aids in surgical planning, guiding pathological diagnosis and treatment strategies in the recurrent setting. This study further confirms the reliability of single-dose LFA-based DSC-MRI.

ABBREVIATIONS

LFA = low flip angle; MFA = moderate flip angle; sRCBV = standardized relative cerebral blood volume; FTB = fractional tumor burden; PTRE = post treatment radiation effects; ROC = receiver operating characteristics; CCC = concordance correlation coefficient.

The ISMRM Open Science Initiative for Perfusion Imaging (OSIPI): Results from the OSIPI-Dynamic Contrast-Enhanced challenge

PURPOSE

K trans $$ {K}^{\mathrm{trans}} $$ has often been proposed as a quantitative imaging biomarker for diagnosis, prognosis, and treatment response assessment for various tumors. None of the many software tools forK trans $$ {K}^{\mathrm{trans}} $$ quantification are standardized. The ISMRM Open Science Initiative for Perfusion Imaging-Dynamic Contrast-Enhanced (OSIPI-DCE) challenge was designed to benchmark methods to better help the efforts to standardizeK trans $$ {K}^{\mathrm{trans}} $$ measurement.

METHODS

A framework was created to evaluateK trans $$ {K}^{\mathrm{trans}} $$ values produced by DCE-MRI analysis pipelines to enable benchmarking. The perfusion MRI community was invited to apply their pipelines forK trans $$ {K}^{\mathrm{trans}} $$ quantification in glioblastoma from clinical and synthetic patients. Submissions were required to include the entrants'K trans $$ {K}^{\mathrm{trans}} $$ values, the applied software, and a standard operating procedure. These were evaluated using the proposedOSIP I gold $$ \mathrm{OSIP}{\mathrm{I}}_{\mathrm{gold}} $$ score defined with accuracy, repeatability, and reproducibility components.

RESULTS

Across the 10 received submissions, theOSIP I gold $$ \mathrm{OSIP}{\mathrm{I}}_{\mathrm{gold}} $$ score ranged from 28% to 78% with a 59% median. The accuracy, repeatability, and reproducibility scores ranged from 0.54 to 0.92, 0.64 to 0.86, and 0.65 to 1.00, respectively (0-1 = lowest-highest). Manual arterial input function selection markedly affected the reproducibility and showed greater variability inK trans $$ {K}^{\mathrm{trans}} $$ analysis than automated methods. Furthermore, provision of a detailed standard operating procedure was critical for higher reproducibility.

CONCLUSIONS

This study reports results from the OSIPI-DCE challenge and highlights the high inter-software variability withinK trans $$ {K}^{\mathrm{trans}} $$ estimation, providing a framework for ongoing benchmarking against the scores presented. Through this challenge, the participating teams were ranked based on the performance of their software tools in the particular setting of this challenge. In a real-world clinical setting, many of these tools may perform differently with different benchmarking methodology.

Quantifying intra-tumoral genetic heterogeneity of glioblastoma toward precision medicine using MRI and a data-inclusive machine learning algorithm

BACKGROUND AND OBJECTIVE

Glioblastoma (GBM) is one of the most aggressive and lethal human cancers. Intra-tumoral genetic heterogeneity poses a significant challenge for treatment. Biopsy is invasive, which motivates the development of non-invasive, MRI-based machine learning (ML) models to quantify intra-tumoral genetic heterogeneity for each patient. This capability holds great promise for enabling better therapeutic selection to improve patient outcome.

METHODS

We proposed a novel Weakly Supervised Ordinal Support Vector Machine (WSO-SVM) to predict regional genetic alteration status within each GBM tumor using MRI. WSO-SVM was applied to a unique dataset of 318 image-localized biopsies with spatially matched multiparametric MRI from 74 GBM patients. The model was trained to predict the regional genetic alteration of three GBM driver genes (EGFR, PDGFRA and PTEN) based on features extracted from the corresponding region of five MRI contrast images. For comparison, a variety of existing ML algorithms were also applied. Classification accuracy of each gene were compared between the different algorithms. The SHapley Additive exPlanations (SHAP) method was further applied to compute contribution scores of different contrast images. Finally, the trained WSO-SVM was used to generate prediction maps within the tumoral area of each patient to help visualize the intra-tumoral genetic heterogeneity.

RESULTS

WSO-SVM achieved 0.80 accuracy, 0.79 sensitivity, and 0.81 specificity for classifying EGFR; 0.71 accuracy, 0.70 sensitivity, and 0.72 specificity for classifying PDGFRA; 0.80 accuracy, 0.78 sensitivity, and 0.83 specificity for classifying PTEN; these results significantly outperformed the existing ML algorithms. Using SHAP, we found that the relative contributions of the five contrast images differ between genes, which are consistent with findings in the literature. The prediction maps revealed extensive intra-tumoral region-to-region heterogeneity within each individual tumor in terms of the alteration status of the three genes.

CONCLUSIONS

This study demonstrated the feasibility of using MRI and WSO-SVM to enable non-invasive prediction of intra-tumoral regional genetic alteration for each GBM patient, which can inform future adaptive therapies for individualized oncology.

Advanced Imaging in the Diagnosis and Response Assessment of High-Grade Glioma: AJR Expert Panel Narrative Review

This AJR Expert Panel Narrative explores the current status of advanced MRI and PET techniques for the post-therapeutic response assessment of high-grade adult-type gliomas, focusing on ongoing clinical controversies in current practice. Discussed techniques that complement conventional MRI and aid the differentiation of recurrent tumor from post-treatment effects include DWI and diffusion tensor imaging; perfusion MRI techniques including dynamic susceptibility contrast (DSC), dynamic contrast-enhanced MRI, and arterial spin labeling; MR spectroscopy including assessment of 2-hydroxyglutarate (2HG) concentration; glucose- and amino acid (AA)-based PET; and amide proton transfer imaging. Updated criteria for Response Assessment in Neuro-Oncology are presented. Given the abundant supporting clinical evidence, the panel supports a recommendation that routine response assessment after HGG treatment should include perfusion MRI, particularly given the development of a consensus recommended DSC-MRI protocol. Although published studies support 2HG MRS and AA PET, these techniques' widespread adoption will likely require increased availability (for 2HG MRS) or increased insurance funding in the United States (for AA PET). The article concludes with a series of consensus opinions from the author panel, centered on the clinical integration of the advanced imaging techniques into posttreatment surveillance protocols.

Image-localized biopsy mapping of brain tumor heterogeneity: A single-center study protocol

Brain cancers pose a novel set of difficulties due to the limited accessibility of human brain tumor tissue. For this reason, clinical decision-making relies heavily on MR imaging interpretation, yet the mapping between MRI features and underlying biology remains ambiguous. Standard (clinical) tissue sampling fails to capture the full heterogeneity of the disease. Biopsies are required to obtain a pathological diagnosis and are predominantly taken from the tumor core, which often has different traits to the surrounding invasive tumor that typically leads to recurrent disease. One approach to solving this issue is to characterize the spatial heterogeneity of molecular, genetic, and cellular features of glioma through the intraoperative collection of multiple image-localized biopsy samples paired with multi-parametric MRIs. We have adopted this approach and are currently actively enrolling patients for our 'Image-Based Mapping of Brain Tumors' study. Patients are eligible for this research study (IRB #16-002424) if they are 18 years or older and undergoing surgical intervention for a brain lesion. Once identified, candidate patients receive dynamic susceptibility contrast (DSC) perfusion MRI and diffusion tensor imaging (DTI), in addition to standard sequences (T1, T1Gd, T2, T2-FLAIR) at their presurgical scan. During surgery, sample anatomical locations are tracked using neuronavigation. The collected specimens from this research study are used to capture the intra-tumoral heterogeneity across brain tumors including quantification of genetic aberrations through whole-exome and RNA sequencing as well as other tissue analysis techniques. To date, these data (made available through a public portal) have been used to generate, test, and validate predictive regional maps of the spatial distribution of tumor cell density and/or treatment-related key genetic marker status to identify biopsy and/or treatment targets based on insight from the entire tumor makeup. This type of methodology, when delivered within clinically feasible time frames, has the potential to further inform medical decision-making by improving surgical intervention, radiation, and targeted drug therapy for patients with glioma.

Revealing the biology behind MRI signatures in high grade glioma

Magnetic resonance imaging (MRI) measurements are routinely collected during the treatment of high-grade gliomas (HGGs) to characterize tumor boundaries and guide surgical tumor resection. Using spatially matched MRI and transcriptomics we discovered HGG tumor biology captured by MRI measurements. We strategically overlaid the spatially matched omics characterizations onto a pre-existing transcriptional map of glioblastoma multiforme (GBM) to enhance the robustness of our analyses. We discovered that T1+C measurements, designed to capture vasculature and blood brain barrier (BBB) breakdown and subsequent contrast extravasation, also indirectly reveal immune cell infiltration. The disruption of the vasculature and BBB within the tumor creates a permissive infiltrative environment that enables the transmigration of anti-inflammatory macrophages into tumors. These relationships were validated through histology and enrichment of genes associated with immune cell transmigration and proliferation. Additionally, T2-weighted (T2W) and mean diffusivity (MD) measurements were associated with angiogenesis and validated using histology and enrichment of genes involved in neovascularization. Furthermore, we establish an unbiased approach for identifying additional linkages between MRI measurements and tumor biology in future studies, particularly with the integration of novel MRI techniques. Lastly, we illustrated how noninvasive MRI can be used to map HGG biology spatially across a tumor, and this provides a platform to develop diagnostics, prognostics, or treatment efficacy biomarkers to improve patient outcomes.

An image-based modeling framework for predicting spatiotemporal brain cancer biology within individual patients

Imaging is central to the clinical surveillance of brain tumors yet it provides limited insight into a tumor's underlying biology. Machine learning and other mathematical modeling approaches can leverage paired magnetic resonance images and image-localized tissue samples to predict almost any characteristic of a tumor. Image-based modeling takes advantage of the spatial resolution of routine clinical scans and can be applied to measure biological differences within a tumor, changes over time, as well as the variance between patients. This approach is non-invasive and circumvents the intrinsic challenges of inter- and intratumoral heterogeneity that have historically hindered the complete assessment of tumor biology and treatment responsiveness. It can also reveal tumor characteristics that may guide both surgical and medical decision-making in real-time. Here we describe a general framework for the acquisition of image-localized biopsies and the construction of spatiotemporal radiomics models, as well as case examples of how this approach may be used to address clinically relevant questions.

Integrated molecular and multiparametric MRI mapping of high-grade glioma identifies regional biologic signatures

Sampling restrictions have hindered the comprehensive study of invasive non-enhancing (NE) high-grade glioma (HGG) cell populations driving tumor progression. Here, we present an integrated multi-omic analysis of spatially matched molecular and multi-parametric magnetic resonance imaging (MRI) profiling across 313 multi-regional tumor biopsies, including 111 from the NE, across 68 HGG patients. Whole exome and RNA sequencing uncover unique genomic alterations to unresectable invasive NE tumor, including subclonal events, which inform genomic models predictive of geographic evolution. Infiltrative NE tumor is alternatively enriched with tumor cells exhibiting neuronal or glycolytic/plurimetabolic cellular states, two principal transcriptomic pathway-based glioma subtypes, which respectively demonstrate abundant private mutations or enrichment in immune cell signatures. These NE phenotypes are non-invasively identified through normalized K2 imaging signatures, which discern cell size heterogeneity on dynamic susceptibility contrast (DSC)-MRI. NE tumor populations predicted to display increased cellular proliferation by mean diffusivity (MD) MRI metrics are uniquely associated with EGFR amplification and CDKN2A homozygous deletion. The biophysical mapping of infiltrative HGG potentially enables the clinical recognition of tumor subpopulations with aggressive molecular signatures driving tumor progression, thereby informing precision medicine targeting.

A cross-sectional study to test equivalence of low- versus intermediate-flip angle dynamic susceptibility contrast MRI measures of relative cerebral blood volume in patients with high-grade gliomas at 1.5 Tesla field strength

Introduction

1.5 Tesla (1.5T) remain a significant field strength for brain imaging worldwide. Recent computer simulations and clinical studies at 3T MRI have suggested that dynamic susceptibility contrast (DSC) MRI using a 30° flip angle ("low-FA") with model-based leakage correction and no gadolinium-based contrast agent (GBCA) preload provides equivalent relative cerebral blood volume (rCBV) measurements to the reference-standard acquisition using a single-dose GBCA preload with a 60° flip angle ("intermediate-FA") and model-based leakage correction. However, it remains unclear whether this holds true at 1.5T. The purpose of this study was to test this at 1.5T in human high-grade glioma (HGG) patients.

Methods

This was a single-institution cross-sectional study of patients who had undergone 1.5T MRI for HGG. DSC-MRI consisted of gradient-echo echo-planar imaging (GRE-EPI) with a low-FA without preload (30°/P-); this then subsequently served as a preload for the standard intermediate-FA acquisition (60°/P+). Both normalized (nrCBV) and standardized relative cerebral blood volumes (srCBV) were calculated using model-based leakage correction (C+) with IBNeuro™ software. Whole-enhancing lesion mean and median nrCBV and srCBV from the low- and intermediate-FA methods were compared using the Pearson's, Spearman's and intraclass correlation coefficients (ICC).

Results

Twenty-three HGG patients composing a total of 31 scans were analyzed. The Pearson and Spearman correlations and ICCs between the 30°/P-/C+ and 60°/P+/C+ acquisitions demonstrated high correlations for both mean and median nrCBV and srCBV.

Conclusion

Our study provides preliminary evidence that for HGG patients at 1.5T MRI, a low FA, no preload DSC-MRI acquisition can be an appealing alternative to the reference standard higher FA acquisition that utilizes a preload.

Image-based models of T-cell distribution identify a clinically meaningful response to a dendritic cell vaccine in patients with glioblastoma

BACKGROUND

Glioblastoma is an extraordinarily heterogeneous tumor, yet the current treatment paradigm is a "one size fits all" approach. Hundreds of glioblastoma clinical trials have been deemed failures because they did not extend median survival, but these cohorts are comprised of patients with diverse tumors. Current methods of assessing treatment efficacy fail to fully account for this heterogeneity.

METHODS

Using an image-based modeling approach, we predicted T-cell abundance from serial MRIs of patients enrolled in the dendritic cell (DC) vaccine clinical trial. T-cell predictions were quantified in both the contrast-enhancing and non-enhancing regions of the imageable tumor, and changes over time were assessed.

RESULTS

A subset of patients in a DC vaccine clinical trial, who had previously gone undetected, were identified as treatment responsive and benefited from prolonged survival. A mere two months after initial vaccine administration, responsive patients had a decrease in model-predicted T-cells within the contrast-enhancing region, with a simultaneous increase in the T2/FLAIR region.

CONCLUSIONS

In a field that has yet to see breakthrough therapies, these results highlight the value of machine learning in enhancing clinical trial assessment, improving our ability to prospectively prognosticate patient outcomes, and advancing the pursuit towards individualized medicine.

Primary central nervous system lymphoma

Primary central nervous system lymphoma (PCNSL) is a diffuse large B cell lymphoma in which the brain, spinal cord, leptomeninges and/or eyes are exclusive sites of disease. Pathophysiology is incompletely understood, although a central role seems to comprise immunoglobulins binding to self-proteins expressed in the central nervous system (CNS) and alterations of genes involved in B cell receptor, Toll-like receptor and NF-κB signalling. Other factors such as T cells, macrophages or microglia, endothelial cells, chemokines, and interleukins, probably also have important roles. Clinical presentation varies depending on the involved regions of the CNS. Standard of care includes methotrexate-based polychemotherapy followed by age-tailored thiotepa-based conditioned autologous stem cell transplantation and, in patients unsuitable for such treatment, consolidation with whole-brain radiotherapy or single-drug maintenance. Personalized treatment, primary radiotherapy and only supportive care should be considered in unfit, frail patients. Despite available treatments, 15-25% of patients do not respond to chemotherapy and 25-50% relapse after initial response. Relapse rates are higher in older patients, although the prognosis of patients experiencing relapse is poor independent of age. Further research is needed to identify diagnostic biomarkers, treatments with higher efficacy and less neurotoxicity, strategies to improve the penetration of drugs into the CNS, and roles of other therapies such as immunotherapies and adoptive cell therapies.

Identification of single-dose, dual-echo based CBV threshold for fractional tumor burden mapping in recurrent glioblastoma

Background

Relative cerebral blood volume (rCBV) obtained from dynamic susceptibility contrast (DSC) MRI is widely used to distinguish high grade glioma recurrence from post treatment radiation effects (PTRE). Application of rCBV thresholds yield maps to distinguish between regional tumor burden and PTRE, a biomarker termed the fractional tumor burden (FTB). FTB is generally measured using conventional double-dose, single-echo DSC-MRI protocols; recently, a single-dose, dual-echo DSC-MRI protocol was clinically validated by direct comparison to the conventional double-dose, single-echo protocol. As the single-dose, dual-echo acquisition enables reduction in the contrast agent dose and provides greater pulse sequence parameter flexibility, there is a compelling need to establish dual-echo DSC-MRI based FTB mapping. In this study, we determine the optimum standardized rCBV threshold for the single-dose, dual-echo protocol to generate FTB maps that best match those derived from the reference standard, double-dose, single-echo protocol.

Methods

The study consisted of 23 high grade glioma patients undergoing perfusion scans to confirm suspected tumor recurrence. We sequentially acquired single dose, dual-echo and double dose, single-echo DSC-MRI data. For both protocols, we generated leakage-corrected standardized rCBV maps. Standardized rCBV (sRCBV) thresholds of 1.0 and 1.75 were used to compute single-echo FTB maps as the reference for delineating PTRE (sRCBV < 1.0), tumor with moderate angiogenesis (1.0 < sRCBV < 1.75), and tumor with high angiogenesis (sRCBV > 1.75) regions. To assess the sRCBV agreement between acquisition protocols, the concordance correlation coefficient (CCC) was computed between the mean tumor sRCBV values across the patients. A receiver operating characteristics (ROC) analysis was performed to determine the optimum dual-echo sRCBV threshold. The sensitivity, specificity, and accuracy were compared between the obtained optimized threshold (1.64) and the standard reference threshold (1.75) for the dual-echo sRCBV threshold.

Results

The mean tumor sRCBV values across the patients showed a strong correlation (CCC = 0.96) between the two protocols. The ROC analysis showed maximum accuracy at thresholds of 1.0 (delineate PTRE from tumor) and 1.64 (differentiate aggressive tumors). The reference threshold (1.75) and the obtained optimized threshold (1.64) yielded similar accuracy, with slight differences in sensitivity and specificity which were not statistically significant (1.75 threshold: Sensitivity = 81.94%; Specificity: 87.23%; Accuracy: 84.58% and 1.64 threshold: Sensitivity = 84.48%; Specificity: 84.97%; Accuracy: 84.73%).

Conclusions

The optimal sRCBV threshold for single-dose, dual-echo protocol was found to be 1.0 and 1.64 for distinguishing tumor recurrence from PTRE; however, minimal differences were observed when using the standard threshold (1.75) as the upper threshold, suggesting that the standard threshold could be used for both protocols. While the prior study validated the agreement of the mean sRCBV values between the protocols, this study confirmed that their voxel-wise agreement is suitable for reliable FTB mapping. Dual-echo DSC-MRI acquisitions enable robust single-dose sRCBV and FTB mapping, provide pulse sequence parameter flexibility and should improve reproducibility by mitigating variations in preload dose and incubation time.

Brain metastases: A Society for Neuro-Oncology (SNO) consensus review on current management and future directions

Brain metastases occur commonly in patients with advanced solid malignancies. Yet, less is known about brain metastases than cancer-related entities of similar incidence. Advances in oncologic care have heightened the importance of intracranial management. Here, in this consensus review supported by the Society for Neuro-Oncology (SNO), we review the landscape of brain metastases with particular attention to management approaches and ongoing efforts with potential to shape future paradigms of care. Each coauthor carried an area of expertise within the field of brain metastases and initially composed, edited, or reviewed their specific subsection of interest. After each subsection was accordingly written, multiple drafts of the manuscript were circulated to the entire list of authors for group discussion and feedback. The hope is that the these consensus guidelines will accelerate progress in the understanding and management of patients with brain metastases, and highlight key areas in need of further exploration that will lead to dedicated trials and other research investigations designed to advance the field.

Knowledge-infused Global-Local Data Fusion for Spatial Predictive Modeling in Precision Medicine

The automated capability of generating spatial prediction for a variable of interest is desirable in various science and engineering domains. Take Precision Medicine of cancer as an example, in which the goal is to match patients with treatments based on molecular markers identified in each patient's tumor. A substantial challenge, however, is that the molecular markers can vary significantly at different spatial locations of a tumor. If this spatial distribution could be predicted, the precision of cancer treatment could be greatly improved by adapting treatment to the spatial molecular heterogeneity. This is a challenging task because no technology is available to measure the molecular markers at each spatial location within a tumor. Biopsy samples provide direct measurement, but they are scarce/local. Imaging, such as MRI, is global, but it only provides proxy/indirect measurement. Also available are mechanistic models or domain knowledge, which are often approximate or incomplete. This paper proposes a novel machine learning framework to fuse the three sources of data/information to generate spatial prediction, namely the knowledge-infused global-local data fusion (KGL) model. A novel mathematical formulation is proposed and solved with theoretical study. We present a real-data application of predicting the spatial distribution of Tumor Cell Density (TCD)-an important molecular marker for brain cancer. A total of 82 biopsy samples were acquired from 18 patients with glioblastoma, together with 6 MRI contrast images from each patient and biological knowledge encoded by a PDE simulator-based mechanistic model called Proliferation-Invasion (PI). KGL achieved the highest prediction accuracy and minimum prediction uncertainty compared with a variety of competing methods. The result has important implications for providing individualized, spatially-optimized treatment for each patient.

Weakly Supervised Skull Stripping of Magnetic Resonance Imaging of Brain Tumor Patients

Automatic brain tumor segmentation is particularly challenging on magnetic resonance imaging (MRI) with marked pathologies, such as brain tumors, which usually cause large displacement, abnormal appearance, and deformation of brain tissue. Despite an abundance of previous literature on learning-based methodologies for MRI segmentation, few works have focused on tackling MRI skull stripping of brain tumor patient data. This gap in literature can be associated with the lack of publicly available data (due to concerns about patient identification) and the labor-intensive nature of generating ground truth labels for model training. In this retrospective study, we assessed the performance of Dense-Vnet in skull stripping brain tumor patient MRI trained on our large multi-institutional brain tumor patient dataset. Our data included pretreatment MRI of 668 patients from our in-house institutional review board-approved multi-institutional brain tumor repository. Because of the absence of ground truth, we used imperfect automatically generated training labels using SPM12 software. We trained the network using common MRI sequences in oncology: T1-weighted with gadolinium contrast, T2-weighted fluid-attenuated inversion recovery, or both. We measured model performance against 30 independent brain tumor test cases with available manual brain masks. All images were harmonized for voxel spacing and volumetric dimensions before model training. Model training was performed using the modularly structured deep learning platform NiftyNet that is tailored toward simplifying medical image analysis. Our proposed approach showed the success of a weakly supervised deep learning approach in MRI brain extraction even in the presence of pathology. Our best model achieved an average Dice score, sensitivity, and specificity of, respectively, 94.5, 96.4, and 98.5% on the multi-institutional independent brain tumor test set. To further contextualize our results within existing literature on healthy brain segmentation, we tested the model against healthy subjects from the benchmark LBPA40 dataset. For this dataset, the model achieved an average Dice score, sensitivity, and specificity of 96.2, 96.6, and 99.2%, which are, although comparable to other publications, slightly lower than the performance of models trained on healthy patients. We associate this drop in performance with the use of brain tumor data for model training and its influence on brain appearance.

Advanced MRI Protocols to Discriminate Glioma From Treatment Effects: State of the Art and Future Directions

In the follow-up treatment of high-grade gliomas (HGGs), differentiating true tumor progression from treatment-related effects, such as pseudoprogression and radiation necrosis, presents an ongoing clinical challenge. Conventional MRI with and without intravenous contrast serves as the clinical benchmark for the posttreatment surveillance imaging of HGG. However, many advanced imaging techniques have shown promise in helping better delineate the findings in indeterminate scenarios, as posttreatment effects can often mimic true tumor progression on conventional imaging. These challenges are further confounded by the histologic admixture that can commonly occur between tumor growth and treatment-related effects within the posttreatment bed. This review discusses the current practices in the surveillance imaging of HGG and the role of advanced imaging techniques, including perfusion MRI and metabolic MRI.

Shape matters: morphological metrics of glioblastoma imaging abnormalities as biomarkers of prognosis

Lacunarity, a quantitative morphological measure of how shapes fill space, and fractal dimension, a morphological measure of the complexity of pixel arrangement, have shown relationships with outcome across a variety of cancers. However, the application of these metrics to glioblastoma (GBM), a very aggressive primary brain tumor, has not been fully explored. In this project, we computed lacunarity and fractal dimension values for GBM-induced abnormalities on clinically standard magnetic resonance imaging (MRI). In our patient cohort (n = 402), we connect these morphological metrics calculated on pretreatment MRI with the survival of patients with GBM. We calculated lacunarity and fractal dimension on necrotic regions (n = 390), all abnormalities present on T1Gd MRI (n = 402), and abnormalities present on T2/FLAIR MRI (n = 257). We also explored the relationship between these metrics and age at diagnosis, as well as abnormality volume. We found statistically significant relationships to outcome for all three imaging regions that we tested, with the shape of T2/FLAIR abnormalities that are typically associated with edema showing the strongest relationship with overall survival. This link between morphological and survival metrics could be driven by underlying biological phenomena, tumor location or microenvironmental factors that should be further explored.

Evaluation of single bolus, dual-echo dynamic susceptibility contrast MRI protocols in brain tumor patients

Relative cerebral blood volume (rCBV) obtained from dynamic susceptibility contrast (DSC) MRI is adversely impacted by contrast agent leakage in brain tumors. Using simulations, we previously demonstrated that multi-echo DSC-MRI protocols provide improvements in contrast agent dosing, pulse sequence flexibility, and rCBV accuracy. The purpose of this study is to assess the in-vivo performance of dual-echo acquisitions in patients with brain tumors (n = 59). To verify pulse sequence flexibility, four single-dose dual-echo acquisitions were tested with variations in contrast agent dose, flip angle, and repetition time, and the resulting dual-echo rCBV was compared to standard single-echo rCBV obtained with preload (double-dose). Dual-echo rCBV was comparable to standard double-dose single-echo protocols (mean (standard deviation) tumor rCBV 2.17 (1.28) vs. 2.06 (1.20), respectively). High rCBV similarity was observed (CCC = 0.96), which was maintained across both flip angle (CCC = 0.98) and repetition time (CCC = 0.96) permutations, demonstrating that dual-echo acquisitions provide flexibility in acquisition parameters. Furthermore, a single dual-echo acquisition was shown to enable quantification of both perfusion and permeability metrics. In conclusion, single-dose dual-echo acquisitions provide similar rCBV to standard double-dose single-echo acquisitions, suggesting contrast agent dose can be reduced while providing significant pulse sequence flexibility and complementary tumor perfusion and permeability metrics.

Consensus recommendations for MRI and PET imaging of primary central nervous system lymphoma: guideline statement from the International Primary CNS Lymphoma Collaborative Group (IPCG)

Advanced molecular and pathophysiologic characterization of primary central nervous system lymphoma (PCNSL) has revealed insights into promising targeted therapeutic approaches. Medical imaging plays a fundamental role in PCNSL diagnosis, staging, and response assessment. Institutional imaging variation and inconsistent clinical trial reporting diminishes the reliability and reproducibility of clinical response assessment. In this context, we aimed to: (1) critically review the use of advanced positron emission tomography (PET) and magnetic resonance imaging (MRI) in the setting of PCNSL; (2) provide results from an international survey of clinical sites describing the current practices for routine and advanced imaging, and (3) provide biologically based recommendations from the International PCNSL Collaborative Group (IPCG) on adaptation of standardized imaging practices. The IPCG provides PET and MRI consensus recommendations built upon previous recommendations for standardized brain tumor imaging protocols (BTIP) in primary and metastatic disease. A biologically integrated approach is provided to addresses the unique challenges associated with the imaging assessment of PCNSL. Detailed imaging parameters facilitate the adoption of these recommendations by researchers and clinicians. To enhance clinical feasibility, we have developed both “ideal” and “minimum standard” protocols at 3T and 1.5T MR systems that will facilitate widespread adoption.

Determination of posterolateral oropharyngeal wall thickness and the potential implications for transoral surgical margins in tonsil cancer

BACKGROUND

Margins in transoral surgery for tonsil cancer can be limited by oropharyngeal wall thickness (OWT), but the normal range is not well established.

METHODS

In 240 noncancer subjects, OWT was measured bilaterally in the vicinity of the tonsils with MRI. Statistical analysis was performed to assess for interaction of age, sex, location, and obesity.

RESULTS

Mean(SD) OWT measured 3.4(0.6) mm posteriorly, 3.7(2.0) mm between the styloglossus and stylopharyngeus, and 5.3(0.8) mm laterally. OWT was greater in men, correlated with obesity, decreased posteriorly and laterally in the 60-80 versus 40-59 year age groups, and increased when styloglossus/stylopharyngeus were closer. OWT was <5 mm in 36.7%-97.9% of locations, with the largest percentage below this threshold located posteriorly.

CONCLUSIONS

OWT is frequently <5 mm, particularly in the posterior and intermuscular areas, suggesting that a smaller surgical margin may need to be accepted in transoral tonsil cancer surgery for anatomic reasons.

Uncertainty quantification in the radiogenomics modeling of EGFR amplification in glioblastoma

Radiogenomics uses machine-learning (ML) to directly connect the morphologic and physiological appearance of tumors on clinical imaging with underlying genomic features. Despite extensive growth in the area of radiogenomics across many cancers, and its potential role in advancing clinical decision making, no published studies have directly addressed uncertainty in these model predictions. We developed a radiogenomics ML model to quantify uncertainty using transductive Gaussian Processes (GP) and a unique dataset of 95 image-localized biopsies with spatially matched MRI from 25 untreated Glioblastoma (GBM) patients. The model generated predictions for regional EGFR amplification status (a common and important target in GBM) to resolve the intratumoral genetic heterogeneity across each individual tumor-a key factor for future personalized therapeutic paradigms. The model used probability distributions for each sample prediction to quantify uncertainty, and used transductive learning to reduce the overall uncertainty. We compared predictive accuracy and uncertainty of the transductive learning GP model against a standard GP model using leave-one-patient-out cross validation. Additionally, we used a separate dataset containing 24 image-localized biopsies from 7 high-grade glioma patients to validate the model. Predictive uncertainty informed the likelihood of achieving an accurate sample prediction. When stratifying predictions based on uncertainty, we observed substantially higher performance in the group cohort (75% accuracy, n = 95) and amongst sample predictions with the lowest uncertainty (83% accuracy, n = 72) compared to predictions with higher uncertainty (48% accuracy, n = 23), due largely to data interpolation (rather than extrapolation). On the separate validation set, our model achieved 78% accuracy amongst the sample predictions with lowest uncertainty. We present a novel approach to quantify radiogenomics uncertainty to enhance model performance and clinical interpretability. This should help integrate more reliable radiogenomics models for improved medical decision-making.

The Practical Application of Emerging Technologies Influencing the Diagnosis and Care of Patients With Primary Brain Tumors

Over the past decade, a variety of new and innovative technologies has led to important advances in the diagnosis and management of patients with primary malignant brain tumors. New approaches to surgical navigation and tumor localization, advanced imaging to define tumor biology and treatment response, and the widespread adoption of a molecularly defined integrated diagnostic paradigm that complements traditional histopathologic diagnosis continue to impact the day-to-day care of these patients. In the neuro-oncology clinic, discussions with patients about the role of tumor treating fields (TTFields) and the incorporation of next-generation sequencing (NGS) data into therapeutic decision-making are now a standard practice. This article summarizes newer applications of technology influencing the pathologic, neuroimaging, neurosurgical, and medical management of patients with malignant primary brain tumors.

Assessment of Prognostic Value of Cystic Features in Glioblastoma Relative to Sex and Treatment With Standard-of-Care

Glioblastoma (GBM) is the most aggressive primary brain tumor and can have cystic components, identifiable through magnetic resonance imaging (MRI). Previous studies suggest that cysts occur in 7–23% of GBMs and report mixed results regarding their prognostic impact. Using our retrospective cohort of 493 patients with first-diagnosis GBM, we carried out an exploratory analysis on this potential link between cystic GBM and survival. Using pretreatment MRIs, we manually identified 88 patients with GBM that had a significant cystic component at presentation and 405 patients that did not. Patients with cystic GBM had significantly longer overall survival and were significantly younger at presentation. Within patients who received the current standard of care (SOC) (N = 184, 40 cystic), we did not observe a survival benefit of cystic GBM. Unexpectedly, we did not observe a significant survival benefit between this SOC cystic cohort and patients with cystic GBM diagnosed before the standard was established (N = 40 with SOC, N = 19 without SOC); this significant SOC benefit was clearly observed in patients with noncystic GBM (N = 144 with SOC, N = 111 without SOC). When stratified by sex, the survival benefit of cystic GBM was only preserved in male patients (N = 303, 47 cystic). We report differences in the absolute and relative sizes of imaging abnormalities on MRI and the prognostic implication of cysts based on sex. We discuss hypotheses for these differences, including the possibility that the presence of a cyst could indicate a less aggressive tumor.

Deep neural network to locate and segment brain tumors outperformed the expert technicians who created the training data

Purpose: Deep learning (DL) algorithms have shown promising results for brain tumor segmentation in MRI. However, validation is required prior to routine clinical use. We report the first randomized and blinded comparison of DL and trained technician segmentations.

Approach: We compiled a multi-institutional database of 741 pretreatment MRI exams. Each contained a postcontrast T1-weighted exam, a T2-weighted fluid-attenuated inversion recovery exam, and at least one technician-derived tumor segmentation. The database included 729 unique patients (470 males and 259 females). Of these exams, 641 were used for training the DL system, and 100 were reserved for testing. We developed a platform to enable qualitative, blinded, controlled assessment of lesion segmentations made by technicians and the DL method. On this platform, 20 neuroradiologists performed 400 side-by-side comparisons of segmentations on 100 test cases. They scored each segmentation between 0 (poor) and 10 (perfect). Agreement between segmentations from technicians and the DL method was also evaluated quantitatively using the Dice coefficient, which produces values between 0 (no overlap) and 1 (perfect overlap).

Results: The neuroradiologists gave technician and DL segmentations mean scores of 6.97 and 7.31, respectively (p<0.00007). The DL method achieved a mean Dice coefficient of 0.87 on the test cases.

Conclusions: This was the first objective comparison of automated and human segmentation using a blinded controlled assessment study. Our DL system learned to outperform its “human teachers” and produced output that was better, on average, than its training data.

Consensus recommendations for a dynamic susceptibility contrast MRI protocol for use in high-grade gliomas

Despite the widespread clinical use of dynamic susceptibility contrast (DSC) MRI, DSC-MRI methodology has not been standardized, hindering its utilization for response assessment in multicenter trials. Recently, the DSC-MRI Standardization Subcommittee of the Jumpstarting Brain Tumor Drug Development Coalition issued an updated consensus DSC-MRI protocol compatible with the standardized brain tumor imaging protocol (BTIP) for high-grade gliomas that is increasingly used in the clinical setting and is the default MRI protocol for the National Clinical Trials Network. After reviewing the basis for controversy over DSC-MRI protocols, this paper provides evidence-based best practices for clinical DSC-MRI as determined by the Committee, including pulse sequence (gradient echo vs spin echo), BTIP-compliant contrast agent dosing (preload and bolus), flip angle (FA), echo time (TE), and post-processing leakage correction. In summary, full-dose preload, full-dose bolus dosing using intermediate (60°) FA and field strength-dependent TE (40-50 ms at 1.5 T, 20-35 ms at 3 T) provides overall best accuracy and precision for cerebral blood volume estimates. When single-dose contrast agent usage is desired, no-preload, full-dose bolus dosing using low FA (30°) and field strength-dependent TE provides excellent performance, with reduced contrast agent usage and elimination of potential systematic errors introduced by variations in preload dose and incubation time.

Identifying the spatial and temporal dynamics of molecularly-distinct glioblastoma sub-populations

Glioblastomas (GBMs) are the most aggressive primary brain tumours and have no known cure. Each individual tumour comprises multiple sub-populations of genetically-distinct cells that may respond differently to targeted therapies and may contribute to disappointing clinical trial results. Image-localized biopsy techniques allow multiple biopsies to be taken during surgery and provide information that identifies regions where particular sub-populations occur within an individual GBM, thus providing insight into their regional genetic variability. These sub-populations may also interact with one another in a competitive or cooperative manner; it is important to ascertain the nature of these interactions, as they may have implications for responses to targeted therapies. We combine genetic information from biopsies with a mechanistic model of interacting GBM sub-populations to characterise the nature of interactions between two commonly occurring GBM sub-populations, those with EGFR and PDGFRA genes amplified. We study population levels found across image-localized biopsy data from a cohort of 25 patients and compare this to model outputs under competitive, cooperative and neutral interaction assumptions. We explore other factors affecting the observed simulated sub-populations, such as selection advantages and phylogenetic ordering of mutations, which may also contribute to the levels of EGFR and PDGFRA amplified populations observed in biopsy data.

Evaluating the Use of rCBV as a Tumor Grade and Treatment Response Classifier Across NCI Quantitative Imaging Network Sites: Part II of the DSC-MRI Digital Reference Object (DRO) Challenge

We have previously characterized the reproducibility of brain tumor relative cerebral blood volume (rCBV) using a dynamic susceptibility contrast magnetic resonance imaging digital reference object across 12 sites using a range of imaging protocols and software platforms. As expected, reproducibility was highest when imaging protocols and software were consistent, but decreased when they were variable. Our goal in this study was to determine the impact of rCBV reproducibility for tumor grade and treatment response classification. We found that varying imaging protocols and software platforms produced a range of optimal thresholds for both tumor grading and treatment response, but the performance of these thresholds was similar. These findings further underscore the importance of standardizing acquisition and analysis protocols across sites and software benchmarking.

Sex-specific impact of patterns of imageable tumor growth on survival of primary glioblastoma patients

Background

Sex is recognized as a significant determinant of outcome among glioblastoma patients, but the relative prognostic importance of glioblastoma features has not been thoroughly explored for sex differences.

Methods

Combining multi-modal MR images, biomathematical models, and patient clinical information, this investigation assesses which pretreatment variables have a sex-specific impact on the survival of glioblastoma patients (299 males and 195 females).

Results

Among males, tumor (T1Gd) radius was a predictor of overall survival (HR = 1.027, p = 0.044). Among females, higher tumor cell net invasion rate was a significant detriment to overall survival (HR = 1.011, p < 0.001). Female extreme survivors had significantly smaller tumors (T1Gd) (p = 0.010 t-test), but tumor size was not correlated with female overall survival (p = 0.955 CPH). Both male and female extreme survivors had significantly lower tumor cell net proliferation rates than other patients (M p = 0.004, F p = 0.001, t-test).

Conclusion

Despite similar distributions of the MR imaging parameters between males and females, there was a sex-specific difference in how these parameters related to outcomes.

Image-based metric of invasiveness predicts response to adjuvant temozolomide for primary glioblastoma

Background

Temozolomide (TMZ) has been the standard-of-care chemotherapy for glioblastoma (GBM) patients for more than a decade. Despite this long time in use, significant questions remain regarding how best to optimize TMZ therapy for individual patients. Understanding the relationship between TMZ response and factors such as number of adjuvant TMZ cycles, patient age, patient sex, and image–based tumor features, might help predict which GBM patients would benefit most from TMZ, particularly for those whose tumors lack O⁶–methylguanine–DNA methyltransferase (MGMT) promoter methylation.

Methods and findings

Using a cohort of 90 newly–diagnosed GBM patients treated according to the standard of care, we examined the relationships between several patient and tumor characteristics and volumetric and survival outcomes during adjuvant chemotherapy. Volumetric changes in MR imaging abnormalities during adjuvant therapy were used to assess TMZ response. T1Gd volumetric response is associated with younger patient age, increased number of TMZ cycles, longer time to nadir volume, and decreased tumor invasiveness. Moreover, increased adjuvant TMZ cycles corresponded with improved volumetric response only among more nodular tumors, and this volumetric response was associated with improved survival outcomes. Finally, in a subcohort of patients with known MGMT methylation status, methylated tumors were more diffusely invasive than unmethylated tumors, suggesting the improved response in nodular tumors is not driven by a preponderance of MGMT methylated tumors.

Conclusions

Our finding that less diffusely invasive tumors are associated with greater volumetric response to TMZ suggests patients with these tumors may benefit from additional adjuvant TMZ cycles, even for those without MGMT methylation.

Performance of Standardized Relative CBV for Quantifying Regional Histologic Tumor Burden in Recurrent High-Grade Glioma: Comparison against Normalized Relative CBV Using Image-Localized Stereotactic Biopsies

BACKGROUND AND PURPOSE

Perfusion MR imaging measures of relative CBV can distinguish recurrent tumor from posttreatment radiation effects in high-grade gliomas. Currently, relative CBV measurement requires normalization based on user-defined reference tissues. A recently proposed method of relative CBV standardization eliminates the need for user input. This study compares the predictive performance of relative CBV standardization against relative CBV normalization for quantifying recurrent tumor burden in high-grade gliomas relative to posttreatment radiation effects.

MATERIALS AND METHODS

We recruited 38 previously treated patients with high-grade gliomas (World Health Organization grades III or IV) undergoing surgical re-resection for new contrast-enhancing lesions concerning for recurrent tumor versus posttreatment radiation effects. We recovered 112 image-localized biopsies and quantified the percentage of histologic tumor content versus posttreatment radiation effects for each sample. We measured spatially matched normalized and standardized relative CBV metrics (mean, median) and fractional tumor burden for each biopsy. We compared relative CBV performance to predict tumor content, including the Pearson correlation (r), against histologic tumor content (0%-100%) and the receiver operating characteristic area under the curve for predicting high-versus-low tumor content using binary histologic cutoffs (≥50%; ≥80% tumor).

RESULTS

Across relative CBV metrics, fractional tumor burden showed the highest correlations with tumor content (0%-100%) for normalized (r = 0.63, P < .001) and standardized (r = 0.66, P < .001) values. With binary cutoffs (ie, ≥50%; ≥80% tumor), predictive accuracies were similar for both standardized and normalized metrics and across relative CBV metrics. Median relative CBV achieved the highest area under the curve (normalized = 0.87, standardized = 0.86) for predicting ≥50% tumor, while fractional tumor burden achieved the highest area under the curve (normalized = 0.77, standardized = 0.80) for predicting ≥80% tumor.

CONCLUSIONS

Standardization of relative CBV achieves similar performance compared with normalized relative CBV and offers an important step toward workflow optimization and consensus methodology.

Imaging of intratumoral heterogeneity in high-grade glioma

High-grade glioma (HGG), and particularly Glioblastoma (GBM), can exhibit pronounced intratumoral heterogeneity that confounds clinical diagnosis and management. While conventional contrast-enhanced MRI lacks the capability to resolve this heterogeneity, advanced MRI techniques and PET imaging offer a spectrum of physiologic and biophysical image features to improve the specificity of imaging diagnoses. Published studies have shown how integrating these advanced techniques can help better define histologically distinct targets for surgical and radiation treatment planning, and help evaluate the regional heterogeneity of tumor recurrence and response assessment following standard adjuvant therapy. Application of texture analysis and machine learning (ML) algorithms has also enabled the emerging field of radiogenomics, which can spatially resolve the regional and genetically distinct subpopulations that coexist within a single GBM tumor. This review focuses on the latest advances in neuro-oncologic imaging and their clinical applications for the assessment of intratumoral heterogeneity.

Days gained response discriminates treatment response in patients with recurrent glioblastoma receiving bevacizumab-based therapies

Background

Accurate assessments of patient response to therapy are a critical component of personalized medicine. In glioblastoma (GBM), the most aggressive form of brain cancer, tumor growth dynamics are heterogenous across patients, complicating assessment of treatment response. This study aimed to analyze days gained (DG), a burgeoning model-based dynamic metric, for response assessment in patients with recurrent GBM who received bevacizumab-based therapies.

Methods

DG response scores were calculated using volumetric tumor segmentations for patients receiving bevacizumab with and without concurrent cytotoxic therapy (N = 62). Kaplan–Meier and Cox proportional hazards analyses were implemented to examine DG prognostic relationship to overall (OS) and progression-free survival (PFS) from the onset of treatment for recurrent GBM.

Results

In patients receiving concurrent bevacizumab and cytotoxic therapy, Kaplan–Meier analysis showed significant differences in OS and PFS at DG cutoffs consistent with previously identified values from newly diagnosed GBM using T1-weighted gadolinium-enhanced magnetic resonance imaging (T1Gd). DG scores for bevacizumab monotherapy patients only approached significance for PFS. Cox regression showed that increases of 25 DG on T1Gd imaging were significantly associated with a 12.5% reduction in OS hazard for concurrent therapy patients and a 4.4% reduction in PFS hazard for bevacizumab monotherapy patients.

Conclusion

DG has significant meaning in recurrent therapy as a metric of treatment response, even in the context of anti-angiogenic therapies. This provides further evidence supporting the use of DG as an adjunct response metric that quantitatively connects treatment response and clinical outcomes.

ENvironmental Dynamics Underlying Responsive Extreme Survivors (ENDURES) of Glioblastoma: A Multidisciplinary Team-based, Multifactorial Analytical Approach

Although glioblastoma (GBM) is a fatal primary brain cancer with short median survival of 15 months, a small number of patients survive >5 years after diagnosis; they are known as extreme survivors (ES). Because of their rarity, very little is known about what differentiates these outliers from other patients with GBM. For the purpose of identifying unknown drivers of extreme survivorship in GBM, the ENDURES consortium (ENvironmental Dynamics Underlying Responsive Extreme Survivors of GBM) was developed. This consortium is a multicenter collaborative network of investigators focused on the integration of multiple types of clinical data and the creation of patient-specific models of tumor growth informed by radiographic and histologic parameters. Leveraging our combined resources, the goals of the ENDURES consortium are 2-fold: (1) to build a curated, searchable, multilayered repository housing clinical and outcome data on a large cohort of ES patients with GBM; and (2) to leverage the ENDURES repository for new insights into tumor behavior and novel targets for prolonging survival for all patients with GBM. In this article, the authors review the available literature and discuss what is already known about ES. The authors then describe the creation of their consortium and some preliminary results.

Integration of machine learning and mechanistic models accurately predicts variation in cell density of glioblastoma using multiparametric MRI

Glioblastoma (GBM) is a heterogeneous and lethal brain cancer. These tumors are followed using magnetic resonance imaging (MRI), which is unable to precisely identify tumor cell invasion, impairing effective surgery and radiation planning. We present a novel hybrid model, based on multiparametric intensities, which combines machine learning (ML) with a mechanistic model of tumor growth to provide spatially resolved tumor cell density predictions. The ML component is an imaging data-driven graph-based semi-supervised learning model and we use the Proliferation-Invasion (PI) mechanistic tumor growth model. We thus refer to the hybrid model as the ML-PI model. The hybrid model was trained using 82 image-localized biopsies from 18 primary GBM patients with pre-operative MRI using a leave-one-patient-out cross validation framework. A Relief algorithm was developed to quantify relative contributions from the data sources. The ML-PI model statistically significantly outperformed (p < 0.001) both individual models, ML and PI, achieving a mean absolute predicted error (MAPE) of 0.106 ± 0.125 versus 0.199 ± 0.186 (ML) and 0.227 ± 0.215 (PI), respectively. Associated Pearson correlation coefficients for ML-PI, ML, and PI were 0.838, 0.518, and 0.437, respectively. The Relief algorithm showed the PI model had the greatest contribution to the result, emphasizing the importance of the hybrid model in achieving the high accuracy.

Moving Toward a Consensus DSC-MRI Protocol: Validation of a Low-Flip Angle Single-Dose Option as a Reference Standard for Brain Tumors

BACKGROUND AND PURPOSE

DSC-MR imaging using preload, intermediate (60°) flip angle and postprocessing leakage correction has gained traction as a standard methodology. Simulations suggest that DSC-MR imaging with flip angle = 30° and no preload yields relative CBV practically equivalent to the reference standard. This study tested this hypothesis in vivo.

MATERIALS AND METHODS

Eighty-four patients with brain lesions were enrolled in this 3-institution study. Forty-three patients satisfied the inclusion criteria. DSC-MR imaging (3T, single-dose gadobutrol, gradient recalled-echo-EPI, TE = 20-35 ms, TR = 1.2-1.63 seconds) was performed twice for each patient, with flip angle = 30°-35° and no preload (P-), which provided preload (P+) for the subsequent intermediate flip angle = 60°. Normalized relative CBV and standardized relative CBV maps were generated, including postprocessing with contrast agent leakage correction (C+) and without (C-) contrast agent leakage correction. Contrast-enhancing lesion volume, mean relative CBV, and contrast-to-noise ratio obtained with 30°/P-/C-, 30°/P-/C+, and 60°/P+/C- were compared with 60°/P+/C+ using the Lin concordance correlation coefficient and Bland-Altman analysis. Equivalence between the 30°/P-/C+ and 60°/P+/C+ protocols and the temporal SNR for the 30°/P- and 60°/P+ DSC-MR imaging data was also determined.

RESULTS

Compared with 60°/P+/C+, 30°/P-/C+ had closest mean standardized relative CBV (P = .61), highest Lin concordance correlation coefficient (0.96), and lowest Bland-Altman bias (μ = 1.89), compared with 30°/P-/C- (P = .02, Lin concordance correlation coefficient = 0.59, μ = 14.6) and 60°/P+/C- (P = .03, Lin concordance correlation coefficient = 0.88, μ = -10.1) with no statistical difference in contrast-to-noise ratios across protocols. The normalized relative CBV and standardized relative CBV were statistically equivalent at the 10% level using either the 30°/P-/C+ or 60°/P+/C+ protocols. Temporal SNR was not significantly different for 30°/P- and 60°/P+ (P = .06).

CONCLUSIONS

Tumor relative CBV derived from low-flip angle, no-preload DSC-MR imaging with leakage correction is an attractive single-dose alternative to the higher dose reference standard.

Evaluating Multisite rCBV Consistency from DSC-MRI Imaging Protocols and Postprocessing Software Across the NCI Quantitative Imaging Network Sites Using a Digital Reference Object (DRO)

Relative cerebral blood volume (rCBV) cannot be used as a response metric in clinical trials, in part, because of variations in biomarker consistency and associated interpretation across sites, stemming from differences in image acquisition and postprocessing methods (PMs). This study leveraged a dynamic susceptibility contrast magnetic resonance imaging digital reference object to characterize rCBV consistency across 12 sites participating in the Quantitative Imaging Network (QIN), specifically focusing on differences in site-specific imaging protocols (IPs; n = 17), and PMs (n = 19) and differences due to site-specific IPs and PMs (n = 25). Thus, high agreement across sites occurs when 1 managing center processes rCBV despite slight variations in the IP. This result is most likely supported by current initiatives to standardize IPs. However, marked intersite disagreement was observed when site-specific software was applied for rCBV measurements. This study's results have important implications for comparing rCBV values across sites and trials, where variability in PMs could confound the comparison of therapeutic effectiveness and/or any attempts to establish thresholds for categorical response to therapy. To overcome these challenges and ensure the successful use of rCBV as a clinical trial biomarker, we recommend the establishment of qualifying and validating site- and trial-specific criteria for scanners and acquisition methods (eg, using a validated phantom) and the software tools used for dynamic susceptibility contrast magnetic resonance imaging analysis (eg, using a digital reference object where the ground truth is known).

Accurate Patient-Specific Machine Learning Models of Glioblastoma Invasion Using Transfer Learning

BACKGROUND AND PURPOSE

MR imaging-based modeling of tumor cell density can substantially improve targeted treatment of glioblastoma. Unfortunately, interpatient variability limits the predictive ability of many modeling approaches. We present a transfer learning method that generates individualized patient models, grounded in the wealth of population data, while also detecting and adjusting for interpatient variabilities based on each patient's own histologic data.

MATERIALS AND METHODS

We recruited patients with primary glioblastoma undergoing image-guided biopsies and preoperative imaging, including contrast-enhanced MR imaging, dynamic susceptibility contrast MR imaging, and diffusion tensor imaging. We calculated relative cerebral blood volume from DSC-MR imaging and mean diffusivity and fractional anisotropy from DTI. Following image coregistration, we assessed tumor cell density for each biopsy and identified corresponding localized MR imaging measurements. We then explored a range of univariate and multivariate predictive models of tumor cell density based on MR imaging measurements in a generalized one-model-fits-all approach. We then implemented both univariate and multivariate individualized transfer learning predictive models, which harness the available population-level data but allow individual variability in their predictions. Finally, we compared Pearson correlation coefficients and mean absolute error between the individualized transfer learning and generalized one-model-fits-all models.

RESULTS

Tumor cell density significantly correlated with relative CBV (r = 0.33, P < .001), and T1-weighted postcontrast (r = 0.36, P < .001) on univariate analysis after correcting for multiple comparisons. With single-variable modeling (using relative CBV), transfer learning increased predictive performance (r = 0.53, mean absolute error = 15.19%) compared with one-model-fits-all (r = 0.27, mean absolute error = 17.79%). With multivariate modeling, transfer learning further improved performance (r = 0.88, mean absolute error = 5.66%) compared with one-model-fits-all (r = 0.39, mean absolute error = 16.55%).

CONCLUSIONS

Transfer learning significantly improves predictive modeling performance for quantifying tumor cell density in glioblastoma.

Is the blood-brain barrier really disrupted in all glioblastomas? A critical assessment of existing clinical data

The blood-brain barrier (BBB) excludes the vast majority of cancer therapeutics from normal brain. However, the importance of the BBB in limiting drug delivery and efficacy is controversial in high-grade brain tumors, such as glioblastoma (GBM). The accumulation of normally brain impenetrant radiographic contrast material in essentially all GBM has popularized a belief that the BBB is uniformly disrupted in all GBM patients so that consideration of drug distribution across the BBB is not relevant in designing therapies for GBM. However, contrary to this view, overwhelming clinical evidence demonstrates that there is also a clinically significant tumor burden with an intact BBB in all GBM, and there is little doubt that drugs with poor BBB permeability do not provide therapeutically effective drug exposures to this fraction of tumor cells. This review provides an overview of the clinical literature to support a central hypothesis: that all GBM patients have tumor regions with an intact BBB, and cure for GBM will only be possible if these regions of tumor are adequately treated.

Optimization of Acquisition and Analysis Methods for Clinical Dynamic Susceptibility Contrast MRI Using a Population-Based Digital Reference Object

BACKGROUND AND PURPOSE

The accuracy of DSC-MR imaging CBV maps in glioblastoma depends on acquisition and analysis protocols. Multisite protocol heterogeneity has challenged standardization initiatives due to the difficulties of in vivo validation. This study sought to compare the accuracy of routinely used protocols using a digital reference object.

MATERIALS AND METHODS

The digital reference object consisted of approximately 10,000 simulated voxels recapitulating typical signal heterogeneity encountered in vivo. The influence of acquisition and postprocessing methods on CBV reliability was evaluated across 6912 parameter combinations, including contrast agent dosing schemes, pulse sequence parameters, field strengths, and postprocessing methods. Accuracy and precision were assessed using the concordance correlation coefficient and coefficient of variation.

RESULTS

Across all parameter space, the optimal protocol included full-dose contrast agent preload and bolus, intermediate (60°) flip angle, 30-ms TE, and postprocessing with a leakage-correction algorithm (concordance correlation coefficient = 0.97, coefficient of variation = 6.6%). Protocols with no preload or fractional dose preload and bolus using these acquisition parameters were generally less robust. However, a protocol with no preload, full-dose bolus, and low (30°) flip angle performed very well (concordance correlation coefficient = 0.93, coefficient of variation = 8.7% at 1.5T and concordance correlation coefficient = 0.92, coefficient of variation = 8.2% at 3T).

CONCLUSIONS

Schemes with full-dose preload and bolus maximize CBV accuracy and reduce variability, which could enable smaller sample sizes and more reliable detection of CBV changes in clinical trials. When a lower total contrast agent dose is desired, use of a low flip angle, no preload, and full-dose bolus protocol may provide an attractive alternative.

Neuroimaging abnormalities in patients with Cowden syndrome: Retrospective single-center study

Background

We retrospectively reviewed the neuroimaging findings of patients with Cowden syndrome and determined their frequency in a single cohort.

Methods

Electronic medical records were queried from January 1999 to January 2017 to identify patients who fit the clinical criteria for diagnosis of Cowden syndrome with or without a documented PTEN mutation. Patients with brain MRI examinations were then identified.

Results

We retrospectively identified 44 patients with Cowden syndrome, 22 of whom had neuroimaging for review. Eleven (50%) had Lhermitte-Duclos disease, 4 (18.1%) had meningiomas, 13 (59.1%) had at least one developmental venous anomaly, 3 had cavernous malformations, 2 had evidence of dural arteriovenous fistula, 7 had increased white matter signal abnormalities relative to age (31.8%), 4 had prominent perivascular spaces, cerebellar tonsillar ectopia was present in 7 of 21 (33.3%), and 1 had cortical malformation.

Conclusions

It is important to recognize that in addition to Lhermitte-Duclos disease, other intracranial findings such as multiple venous anomalies, meningiomas, greater than expected white matter signal abnormality, prominent perivascular spaces, and cortical malformations may warrant a thorough evaluation for Cowden syndrome in the appropriate clinical setting. We further recommend that this broader spectrum of intracranial abnormalities be considered for addition to the Cowden syndrome diagnostic criteria at the time of next revision.

Characterizing the Influence of Preload Dosing on Percent Signal Recovery (PSR) and Cerebral Blood Volume (CBV) Measurements in a Patient Population With High-Grade Glioma Using Dynamic Susceptibility Contrast MRI

With DSC-MRI, contrast agent leakage effects in brain tumors can either be leveraged for percent signal recovery (PSR) measurements or be adequately resolved for accurate relative cerebral blood volume (rCBV) measurements. Leakage effects can be dimished by administration of a preload dose before imaging and/or specific postprocessing steps. This study compares the consistency of both PSR and rCBV measurements as a function of varying preload doses in a retrospective analysis of 14 subjects with high-grade gliomas. The scans consisted of 6 DSC-MRI scans during 6 sequential bolus injections (0.05 mmol/kg). Mean PSR was calculated for tumor and normal-appearing white matter regions of interest. DSC-MRI data were corrected for leakage effects before computing mean tumor rCBV. Statistical differences were seen across varying preloads for tumor PSR (P value = 4.57E-24). Tumor rCBV values did not exhibit statistically significant differences across preloads (P value = .14) and were found to be highly consistent for clinically relevant preloads (intraclass correlation coefficient = 0.93). For a 0.05 mmol/kg injection bolus and pulse sequence parameters used, the highest PSR contrast between normal-appearing white matter and tumor occurs when no preload is used. This suggests that studies using PSR as a biomarker should acquire DSC-MRI data without preload. The finding that leakage-corrected rCBV values do not depend on the presence or dose of preload contradicts that of previous studies with dissimilar acquisition protocols. This further confirms the sensitivity of rCBV to preload dosing schemes and pulse sequence parameters and highlights the importance of standardization efforts for achieving multisite rCBV consistency.

Dynamic Susceptibility Contrast-MRI Quantification Software Tool: Development and Evaluation

Relative cerebral blood volume (rCBV) is a magnetic resonance imaging biomarker that is used to differentiate progression from pseudoprogression in patients with glioblastoma multiforme, the most common primary brain tumor. However, calculated rCBV depends considerably on the software used. Automating all steps required for rCBV calculation is important, as user interaction can lead to increased variability and possible inaccuracies in clinical decision-making. Here, we present an automated tool for computing rCBV from dynamic susceptibility contrast-magnetic resonance imaging that includes leakage correction. The entrance and exit bolus time points are automatically calculated using wavelet-based detection. The proposed tool is compared with 3 Food and Drug Administration-approved software packages, 1 automatic and 2 requiring user interaction, on a data set of 43 patients. We also evaluate manual and automated white matter (WM) selection for normalization of the cerebral blood volume maps. Our system showed good agreement with 2 of the 3 software packages. The intraclass correlation coefficient for all comparisons between the same software operated by different people was >0.880, except for FuncTool when operated by user 1 versus user 2. Little variability in agreement between software tools was observed when using different WM selection techniques. Our algorithm for automatic rCBV calculation with leakage correction and automated WM selection agrees well with 2 out of the 3 FDA-approved software packages.

Radiogenomics to characterize regional genetic heterogeneity in glioblastoma

BACKGROUND

Glioblastoma (GBM) exhibits profound intratumoral genetic heterogeneity. Each tumor comprises multiple genetically distinct clonal populations with different therapeutic sensitivities. This has implications for targeted therapy and genetically informed paradigms. Contrast-enhanced (CE)-MRI and conventional sampling techniques have failed to resolve this heterogeneity, particularly for nonenhancing tumor populations. This study explores the feasibility of using multiparametric MRI and texture analysis to characterize regional genetic heterogeneity throughout MRI-enhancing and nonenhancing tumor segments.

METHODS

We collected multiple image-guided biopsies from primary GBM patients throughout regions of enhancement (ENH) and nonenhancing parenchyma (so called brain-around-tumor, [BAT]). For each biopsy, we analyzed DNA copy number variants for core GBM driver genes reported by The Cancer Genome Atlas. We co-registered biopsy locations with MRI and texture maps to correlate regional genetic status with spatially matched imaging measurements. We also built multivariate predictive decision-tree models for each GBM driver gene and validated accuracies using leave-one-out-cross-validation (LOOCV).

RESULTS

We collected 48 biopsies (13 tumors) and identified significant imaging correlations (univariate analysis) for 6 driver genes: EGFR, PDGFRA, PTEN, CDKN2A, RB1, and TP53. Predictive model accuracies (on LOOCV) varied by driver gene of interest. Highest accuracies were observed for PDGFRA (77.1%), EGFR (75%), CDKN2A (87.5%), and RB1 (87.5%), while lowest accuracy was observed in TP53 (37.5%). Models for 4 driver genes (EGFR, RB1, CDKN2A, and PTEN) showed higher accuracy in BAT samples (n = 16) compared with those from ENH segments (n = 32).

CONCLUSION

MRI and texture analysis can help characterize regional genetic heterogeneity, which offers potential diagnostic value under the paradigm of individualized oncology.

Genetic diversity in the germplasm of black pepper determined by EST-SSR markers

This study aimed to assess genetic diversity in the germplasm of black pepper from around the world using SSR markers from EST. In total, 13 markers were selected and successfully amplified the target loci across the black pepper germplasm. All the EST-SSR markers showed high levels of polymorphisms with an average polymorphism information content of 0.93. The genetic similarity coefficients among all accessions ranged from 0.724 to 1.000, with an average of 0.867. These results indicated that black pepper germplasms possess a complex genetic background and high genetic diversity. Based on a cluster analysis, 148 black pepper germplasms were grouped in two major clades: the Neotropics and the Asian tropics. Peperomia pellucida was grouped separately and distantly from all other accessions. These results generally agreed with the genetic and geographic distances. However, the Asian tropics clade did not cluster according to their geographic origins. In addition, compared with the American accessions, the Asian wild accessions and cultivated accessions grouped together, indicating a close genetic relationship. This verified the origin of black pepper. The newly developed EST-SSRs are highly valuable resources for the conservation of black pepper germplasm diversity and for black pepper breeding.

Multi-Parametric MRI and Texture Analysis to Visualize Spatial Histologic Heterogeneity and Tumor Extent in Glioblastoma

Background

Genetic profiling represents the future of neuro-oncology but suffers from inadequate biopsies in heterogeneous tumors like Glioblastoma (GBM). Contrast-enhanced MRI (CE-MRI) targets enhancing core (ENH) but yields adequate tumor in only ~60% of cases. Further, CE-MRI poorly localizes infiltrative tumor within surrounding non-enhancing parenchyma, or brain-around-tumor (BAT), despite the importance of characterizing this tumor segment, which universally recurs. In this study, we use multiple texture analysis and machine learning (ML) algorithms to analyze multi-parametric MRI, and produce new images indicating tumor-rich targets in GBM.

Methods

We recruited primary GBM patients undergoing image-guided biopsies and acquired pre-operative MRI: CE-MRI, Dynamic-Susceptibility-weighted-Contrast-enhanced-MRI, and Diffusion Tensor Imaging. Following image coregistration and region of interest placement at biopsy locations, we compared MRI metrics and regional texture with histologic diagnoses of high- vs low-tumor content (≥80% vs <80% tumor nuclei) for corresponding samples. In a training set, we used three texture analysis algorithms and three ML methods to identify MRI-texture features that optimized model accuracy to distinguish tumor content. We confirmed model accuracy in a separate validation set.

Results

We collected 82 biopsies from 18 GBMs throughout ENH and BAT. The MRI-based model achieved 85% cross-validated accuracy to diagnose high- vs low-tumor in the training set (60 biopsies, 11 patients). The model achieved 81.8% accuracy in the validation set (22 biopsies, 7 patients).

Conclusion

Multi-parametric MRI and texture analysis can help characterize and visualize GBM’s spatial histologic heterogeneity to identify regional tumor-rich biopsy targets.

Impact of Software Modeling on the Accuracy of Perfusion MRI in Glioma

BACKGROUND AND PURPOSE

Relative cerebral blood volume, as measured by T2*-weighted dynamic susceptibility-weighted contrast-enhanced MRI, represents the most robust and widely used perfusion MR imaging metric in neuro-oncology. Our aim was to determine whether differences in modeling implementation will impact the correction of leakage effects (from blood-brain barrier disruption) and the accuracy of relative CBV calculations as measured on T2*-weighted dynamic susceptibility-weighted contrast-enhanced MR imaging at 3T field strength.

MATERIALS AND METHODS

This study included 52 patients with glioma undergoing DSC MR imaging. Thirty-six patients underwent both non-preload dose- and preload dose-corrected DSC acquisitions, with 16 patients undergoing preload dose-corrected acquisitions only. For each acquisition, we generated 2 sets of relative CBV metrics by using 2 separate, widely published, FDA-approved commercial software packages: IB Neuro and nordicICE. We calculated 4 relative CBV metrics within tumor volumes: mean relative CBV, mode relative CBV, percentage of voxels with relative CBV > 1.75, and percentage of voxels with relative CBV > 1.0 (fractional tumor burden). We determined Pearson (r) and Spearman (ρ) correlations between non-preload dose- and preload dose-corrected metrics. In a subset of patients with recurrent glioblastoma (n = 25), we determined receiver operating characteristic area under the curve for fractional tumor burden accuracy to predict the tissue diagnosis of tumor recurrence versus posttreatment effect. We also determined correlations between rCBV and microvessel area from stereotactic biopsies (n = 29) in 12 patients.

RESULTS

With IB Neuro, relative CBV metrics correlated highly between non-preload dose- and preload dose-corrected conditions for fractional tumor burden (r = 0.96, ρ = 0.94), percentage > 1.75 (r = 0.93, ρ = 0.91), mean (r = 0.87, ρ = 0.86), and mode (r = 0.78, ρ = 0.76). These correlations dropped substantially with nordicICE. With fractional tumor burden, IB Neuro was more accurate than nordicICE in diagnosing tumor versus posttreatment effect (area under the curve = 0.85 versus 0.67) (P < .01). The highest relative CBV-microvessel area correlations required preload dose and IB Neuro (r = 0.64, ρ = 0.58, P = .001).

CONCLUSIONS

Different implementations of perfusion MR imaging software modeling can impact the accuracy of leakage correction, relative CBV calculation, and correlations with histologic benchmarks.

Biopsy validation of 18F-DOPA PET and biodistribution in gliomas for neurosurgical planning and radiotherapy target delineation: results of a prospective pilot study

BACKGROUND

Delineation of glioma extent for surgical or radiotherapy planning is routinely based on MRI. There is increasing awareness that contrast enhancement on T1-weighted images (T1-CE) may not reflect the entire extent of disease. The amino acid tracer (18)F-DOPA (3,4-dihydroxy-6-[18F] fluoro-l-phenylalanine) has a high tumor-to-background signal and high sensitivity for glioma imaging. This study compares (18)F-DOPA PET against conventional MRI for neurosurgical biopsy targeting, resection planning, and radiotherapy target volume delineation.

METHODS

Conventional MR and (18)F-DOPA PET/CT images were acquired in 10 patients with suspected malignant brain tumors. One to 3 biopsy locations per patient were chosen in regions of concordant and discordant (18)F-DOPA uptake and MR contrast enhancement. Histopathology was reviewed on 23 biopsies. (18)F-DOPA PET was quantified using standardized uptake values (SUV) and tumor-to-normal hemispheric tissue (T/N) ratios.

RESULTS

Pathologic review confirmed glioma in 22 of 23 biopsy specimens. Thirteen of 16 high-grade biopsy specimens were obtained from regions of elevated (18)F-DOPA uptake, while T1-CE was present in only 6 of those 16 samples. Optimal (18)F-DOPA PET thresholds corresponding to high-grade disease based on histopathology were calculated as T/N > 2.0. In every patient, (18)F-DOPA uptake regions with T/N > 2.0 extended beyond T1-CE up to a maximum of 3.5 cm. SUV was found to correlate with grade and cellularity.

CONCLUSIONS

(18)F-DOPA PET SUV(max) may more accurately identify regions of higher-grade/higher-density disease in patients with astrocytomas and will have utility in guiding stereotactic biopsy selection. Using SUV-based thresholds to define high-grade portions of disease may be valuable in delineating radiotherapy boost volumes.