Swanson, Kristin R

Abstract With mean survival from glioblastoma multiforme (GBM) barely edging over 15 months, it become imperative to develop novel tools informing on both progression and treatment efficacy. As the aggressive growth of GBM outstrips available resources, regions of hypoxia develop. Hypoxia can be clinically assessed through [18F]-flouromisonidazole (FMISO) PET. However, PET imaging suffers from low resolution and reconstruction artifact. We present the comparative and predictive results of a multi-stage model utilizing the mathematical biology of cancer and the physics of PET imaging. Our mathematical model of GBM characterizes cancer cells into the interacting phenotypes of normoxic, hypoxic, and necrotic cells supported by vasculature. The pharmacokinetic activity of FMISO is induced in the simulated tissue then acquired and reconstructed by an analytic simulator of the PET process. This multi-stage process creates a patient-specific virtual PET images with characteristic of the clinical PET scan. Extending previous work, we model the full three-dimensional dynamics of tumor progression, pharmacokinetic activity, and FMISO PET for the tumor kinetics of six GBM patients. Virtual and clinical FMISO PET images are compared. Overall hypoxic burden and spatial distribution show strong correspondence between virtual and clinical FMISO PET images for all six patients. Simulated FMISO PET dynamics, with tumor kinetics derived from routine clinical MRI scans, provide a unique and evolving tool giving insight into the biological connections between magnetic resonance and molecular imaging. Citation Format: Joshua Jacobs, Andrea Hawkins-Darrud, Robert Harrison, Sandra Johnston, Paul Kinahan, Kristin Swanson. Correlating magnetic resonance and molecular imaging using three dimensional untreated virtual control. [abstract]. In: Proceedings of the AACR Special Conference on Engineering and Physical Sciences in Oncology; 2016 Jun 25-28; Boston, MA. Philadelphia (PA): AACR; Cancer Res 2017;77(2 Suppl):Abstract nr A10.

Abstract Purpose: Investigate the utility of patient-specific spatial predictions of tumor cell density from a bio-mathematical model. Introduction: Glioblastomas (GBMs) are the most malignant of all primary brain tumors. While it is known there is always a non-detectable portion of the tumor, current techniques of monitoring GBM progression, imaging and initial histological assessment, are not able to reliably estimate the tumor invasion past the enhancing region on T2-Weighted (T2W) imaging. Over the last two decades, a large effort has been made to create a simple patient-specific mathematical model of gliomas. The resulting model, referred to as the Proliferation-Invasion (PI) model, is based on two key parameters, the net growth rate, ρ, and the dispersal coefficient, D. In this model, the ratio of D/ρ is related to degree of invasion and the product D*ρ, is related to the speed of growth. The intuitive understanding provided by this model has been able to provide patient-specific understanding of disease kinetics enabling prediction of outcomes following surgical resection, radiation and the development of a prognostic response metric. Previous literature utilizing this model has been based on the assumption that what is seen on the pretreatment T1-Weighted contrast-enhanced (T1Gd) and T2W, images correspond to an 80% and 16% tumor cell density threshold respectively. This assumption allows for an estimate of D/ρ from a single time point of imaging. While these values were based on extensive experience, for ethical and technical reasons, they have never been rigorously investigated histologically. Recent technological advances have made it possible for surgeons to use an MRI to guide the acquisition of tissue making it possible to know with a good degree of accuracy where on the MR image the histological specimen comes from. Methods: Model Calibration : To estimate D/ρ for each patient, we assume abnormalities on the T1Gd and T2W images correspond to an 80% and 16% tumor cell density threshold respectively. We then utilize a Bayesian calibration approach based on adaptive grid refinement while holding the velocity constant to find the most likely value of D/ρ to match the observed radial measurements. Three-Dimensional Density Maps : Given a gray/white segmentation and an estimate for D/ρ, we can build a tumor cell density prediction in the patient's anatomy using the Eikonal equations and the modified Fast Marching Method (FMM) algorithm presented by Konukoglu et al. Patient Cohort : Eighteen patients were recruited with clinically suspected GBM undergoing preoperative stereotactic MRI for surgical resection with IRB approval Barrow Neurological Institute and Mayo Clinic in Arizona. Surgical Biopsy : Pre-operative conventional MRI, including T1Gd and T2W, was utilized to guide stereotactic biopsies. An average of 5–6 tissue specimens were acquired from each tumor by using stereotactic surgical localization, following the smallest possible diameter craniotomies to minimize brain shift. Histological Analysis : 4 μm tissue sections were stained with hematoxylin and eosin (H&E) for neuropathology review. H&E slides were reviewed by a neuropathologist. The percent tumor nuclei in the field was estimated. Results and Conclusion: Our analysis showed that the intuitive ordering of diffuse to nodular tumor profiles from the PI model is consistent with the histologic data. Further, utilizing the histologic data, limitations to a universal threshold cutoff of the T2/FLAIR region are specified and a better threshold cutoff value is suggested. We hope that this paper will lead to more studies of this kind resulting in more accurate ways to predict the spatial distribution of GBM tumor cells from medical imaging. Such advances could revolutionize surgical procedures and radiation therapy and may enable additional insights into tumor kinetic differences meaningful for treatment. Citation Format: Andrea Hawkins-Daarud, Lauren DeGirolamo, Joshua Jacobs, Kamala Clark-Swanson, Jennifer M. Eschbacher, Kris A. Smith, Peter Nakaji, Leslie C. Baxter, John P. Karis, Teresa Wu, J. Ross Mitchell, Jing Li, Leland Hu, Kristin R. Swanson. Histologic evidence for a bio-mathematical model of glioblastoma invasion. [abstract]. In: Proceedings of the AACR Special Conference on Engineering and Physical Sciences in Oncology; 2016 Jun 25-28; Boston, MA. Philadelphia (PA): AACR; Cancer Res 2017;77(2 Suppl):Abstract nr A08.

Abstract BACKGROUND Extent of resection (EOR) correlates with glioblastoma outcomes. Resectability and EOR depend on anatomical, clinical, and surgeon factors. Resectability likely influences outcome in and of itself, but an accurate measurement of resectability remains elusive. An understanding of resectability and the factors that influence it may provide a means to control a confounder in clinical trials and provide reference for decision making. OBJECTIVE To provide proof of concept of the use of the collective wisdom of experienced brain tumor surgeons in assessing glioblastoma resectability. METHODS We surveyed 13 academic tumor neurosurgeons nationwide to assess the resectability of newly diagnosed glioblastoma. Participants reviewed 20 cases, including digital imaging and communications in medicine-formatted pre- and postoperative magnetic resonance images and clinical vignettes. The selected cases involved a variety of anatomical locations and a range of EOR. Participants were asked about surgical goal, eg, gross total resection, subtotal resection (STR), or biopsy, and rationale for their decision. We calculated a “resectability index” for each lesion by pooling responses from all 13 surgeons. RESULTS Neurosurgeons’ individual surgical goals varied significantly ( P = .015), but the resectability index calculated from the surgeons’ pooled responses was strongly correlated with the percentage of contrast-enhancing residual tumor ( R = 0.817, P < .001). The collective STR goal predicted intraoperative decision of intentional STR documented on operative notes ( P < .01) and nonresectable residual ( P < .01), but not resectable residual. CONCLUSION In this pilot study, we demonstrate the feasibility of measuring the resectability of glioblastoma through crowdsourcing. This tool could be used to quantify resectability, a potential confounder in neuro-oncology clinical trials.

AbstractBackgroundGlioblastoma (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.MethodsWe 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).ResultsWe 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).ConclusionMRI and texture analysis can help characterize regional genetic heterogeneity, which offers potential diagnostic value under the paradigm of individualized oncology.