Search

Background and objectiveGlioblastoma (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.MethodsWe 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.ResultsWSO-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.ConclusionsThis 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.

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.

Glioblastomas are aggressive primary brain tumors known for their inter- and intratumor heterogeneity. This disease is uniformly fatal, with intratumor heterogeneity the major reason for treatment failure and recurrence. Just like the nature vs nurture debate, heterogeneity can arise from intrinsic or environmental influences. Whilst it is impossible to clinically separate observed behavior of cells from their environmental context, using a mathematical framework combined with multiscale data gives us insight into the relative roles of variation from different sources. To better understand the implications of intratumor heterogeneity on therapeutic outcomes, we created a hybrid agent-based mathematical model that captures both the overall tumor kinetics and the individual cellular behavior. We track single cells as agents, cell density on a coarser scale, and growth factor diffusion and dynamics on a finer scale over time and space. Our model parameters were fit utilizing serial MRI imaging and cell tracking data from ex vivo tissue slices acquired from a growth-factor driven glioblastoma murine model. When fitting our model to serial imaging only, there was a spectrum of equally-good parameter fits corresponding to a wide range of phenotypic behaviors. When fitting our model using imaging and cell scale data, we determined that environmental heterogeneity alone is insufficient to match the single cell data, and intrinsic heterogeneity is required to fully capture the migration behavior. The wide spectrum of in silico tumors also had a wide variety of responses to an application of an anti-proliferative treatment. Recurrent tumors were generally less proliferative than pre-treatment tumors as measured via the model simulations and validated from human GBM patient histology. Further, we found that all tumors continued to grow with an anti-migratory treatment alone, but the anti-proliferative/anti-migratory combination generally showed improvement over an anti-proliferative treatment alone. Together our results emphasize the need to better understand the underlying phenotypes and tumor heterogeneity present in a tumor when designing therapeutic regimens.

BackgroundTemozolomide (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 O6-methylguanine-DNA methyltransferase (MGMT) promoter methylation.Methods and findingsUsing 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.ConclusionsOur 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.

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.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

Malignant gliomas are incurable, primary brain neoplasms noted for their potential to extensively invade brain parenchyma. Current methods of clinical imaging do not elucidate the full extent of brain invasion, making it difficult to predict which, if any, patients are likely to benefit from gross total resection. Our goal was to apply a mathematical modeling approach to estimate the overall tumor invasiveness on a patient-by-patient basis and determine whether gross total resection would improve survival in patients with relatively less invasive gliomas.In 243 patients presenting with contrast-enhancing gliomas, estimates of the relative invasiveness of each patient's tumor, in terms of the ratio of net proliferation rate of the glioma cells to their net dispersal rate, were derived by applying a patient-specific mathematical model to routine pretreatment MR imaging. The effect of varying degrees of extent of resection on overall survival was assessed for cohorts of patients grouped by tumor invasiveness.We demonstrate that patients with more diffuse tumors showed no survival benefit (P = 0.532) from gross total resection over subtotal/biopsy, while those with nodular (less diffuse) tumors showed a significant benefit (P = 0.00142) with a striking median survival benefit of over eight months compared to sub-totally resected tumors in the same cohort (an 80% improvement in survival time for GTR only seen for nodular tumors).These results suggest that our patient-specific, model-based estimates of tumor invasiveness have clinical utility in surgical decision making. Quantification of relative invasiveness assessed from routinely obtained pre-operative imaging provides a practical predictor of the benefit of gross total resection.

Accurate clinical assessment of a patient's response to treatment for glioblastoma multiforme (GBM), the most malignant type of primary brain tumor, is undermined by the wide patient-to-patient variability in GBM dynamics and responsiveness to therapy. Using computational models that account for the unique geometry and kinetics of individual patients' tumors, we developed a method for assessing treatment response that discriminates progression-free and overall survival following therapy for GBM. Applying these models as untreated virtual controls, we generate a patient-specific "Days Gained" response metric that estimates the number of days a therapy delayed imageable tumor progression. We assessed treatment response in terms of Days Gained scores for 33 patients at the time of their first MRI scan following first-line radiation therapy. Based on Kaplan-Meier analyses, patients with Days Gained scores of 100 or more had improved progression-free survival, and patients with scores of 117 or more had improved overall survival. Our results demonstrate that the Days Gained response metric calculated at the routinely acquired first post-radiation treatment time point provides prognostic information regarding progression and survival outcomes. Applied prospectively, our model-based approach has the potential to improve GBM treatment by accounting for patient-to-patient heterogeneity in GBM dynamics and responses to therapy.

To demonstrate a method of generating patient-specific, biologically-guided radiotherapy dose plans and compare them to the standard-of-care protocol.We integrated a patient-specific biomathematical model of glioma proliferation, invasion and radiotherapy with a multiobjective evolutionary algorithm for intensity-modulated radiation therapy optimization to construct individualized, biologically-guided plans for 11 glioblastoma patients. Patient-individualized, spherically-symmetric simulations of the standard-of-care and optimized plans were compared in terms of several biological metrics.The integrated model generated spatially non-uniform doses that, when compared to the standard-of-care protocol, resulted in a 67% to 93% decrease in equivalent uniform dose to normal tissue, while the therapeutic ratio, the ratio of tumor equivalent uniform dose to that of normal tissue, increased between 50% to 265%. Applying a novel metric of treatment response (Days Gained) to the patient-individualized simulation results predicted that the optimized plans would have a significant impact on delaying tumor progression, with increases from 21% to 105% for 9 of 11 patients.Patient-individualized simulations using the combination of a biomathematical model with an optimization algorithm for radiation therapy generated biologically-guided doses that decreased normal tissue EUD and increased therapeutic ratio with the potential to improve survival outcomes for treatment of glioblastoma.

Many drugs investigated for the treatment of glioblastoma (GBM) have had disappointing clinical trial results. Efficacy of these agents is dependent on adequate delivery to sensitive tumor cell populations, which is limited by the blood-brain barrier (BBB). Additionally, tumor heterogeneity can lead to subpopulations of cells with different sensitivities to anti-cancer drugs, further impacting therapeutic efficacy. Thus, it may be important to evaluate the extent to which BBB limitations and heterogeneous sensitivity each contribute to a drug's failure. To address this challenge, we developed a minimal mathematical model to characterize these elements of overall drug response, informed by time-series bioluminescence imaging data from a treated patient-derived xenograft (PDX) experimental model. By fitting this mathematical model to a preliminary dataset in a series of nonlinear regression steps, we estimated parameter values for individual PDX subjects that correspond to the dynamics seen in experimental data. Using these estimates as a guide for parameter ranges, we ran model simulations and performed a parameter sensitivity analysis using Latin hypercube sampling and partial rank correlation coefficients. Results from this analysis combined with simulations suggest that BBB permeability may play a slightly greater role in therapeutic efficacy than relative drug sensitivity. Additionally, we discuss recommendations for future experiments based on insights gained from this model. Further research in this area will be vital for improving the development of effective new therapies for glioblastoma patients.

Abstract 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.

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.

Abstract 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.

Abstract 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.

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.

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.

Abstract 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

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.

Despite major drug discovery efforts, the therapeutic options for glioblastoma (GBM) remain inadequate. Here they analyze patient-derived xenograft model of GBM to quantitatively map distribution and cellular response to the EGFR inhibitor erlotinib, and report heterogeneous erlotinib delivery to intracranial tumors to be inadequate to inhibit EGFR signaling.

Glioblastoma, the most aggressive form of primary brain tumor is predominantly assessed with gadolinium-enhanced T1-weighted (T1Gd) and T2-weighted magnetic resonance imaging (MRI). Pixel intensity enhancement on the T1Gd image is understood to correspond to the gadolinium contrast agent leaking from the tumor-induced neovasculature, while hyperintensity on the T2/FLAIR images corresponds with edema and infiltrated tumor cells. None of these modalities directly show tumor cells; rather, they capture abnormalities in the microenvironment caused by the presence of tumor cells. Thus, assessing disease response after treatments impacting the microenvironment remains challenging through the obscuring lens of MR imaging. Anti-angiogenic therapies have been used in the treatment of gliomas with spurious results ranging from no apparent response to significant imaging improvement with the potential for extremely diffuse patterns of tumor recurrence on imaging and autopsy. Anti-angiogenic treatment normalizes the vasculature, effectively decreasing vessel permeability and thus reducing tumor-induced edema, drastically altering T2-weighted MRI. We extend a previously developed mathematical model of glioma growth to explicitly incorporate edema formation allowing us to directly characterize and potentially predict the effects of anti-angiogenics on imageable tumor growth. A comparison of simulated glioma growth and imaging enhancement with and without bevacizumab supports the current understanding that anti-angiogenic treatment can serve as a surrogate for steroids and the clinically-driven hypothesis that anti-angiogenic treatment may not have any significant effect on the growth dynamics of the overall tumor-cell populations. However, the simulations do illustrate a potentially large impact on the level of edematous extracellular fluid, and thus on what would be imageable on T2/FLAIR MR for tumors with lower proliferation rates.

Purpose: Despite the intra- and intertumoral heterogeneity seen in glioblastoma multiforme (GBM), there is little definitive data on the underlying cause of the differences in patient survivals. Serial imaging assessment of tumor growth allows quantification of tumor growth kinetics (TGK) measured in terms of changes in the velocity of radial expansion seen on imaging. Because a systematic study of this entire TGK phenotype-growth before treatment and during each treatment to recurrence -has never been coordinately studied in GBMs, we sought to identify whether patients cluster into discrete groups on the basis of their TGK., Patients and Methods: From our multi-institutional database, we identified 48 patients who underwent maximally safe resection followed by radiotherapy with imaging follow-up through the time of recurrence. The patients were then clustered into two groups through a k-means algorithm taking as input only the TGK before and during treatment., Results: There was a significant survival difference between the clusters ( P = .003). Paradoxically, patients among the long-lived cluster had significantly larger tumors at diagnosis ( P = .027) and faster growth before treatment ( P = .003) but demonstrated a better response to adjuvant chemotherapy ( P = .048). A predictive model was built to identify which cluster patients would likely fall into on the basis of information that would be available to clinicians immediately after radiotherapy (accuracy, 90.3%)., Conclusion: Dichotomizing the heterogeneity of GBMs into two populations-one faster growing yet more responsive with increased survival and one slower growing yet less responsive with shorter survival-suggests that many patients who receive standard-of-care treatments may get better benefit from select alternative treatments.