Association of breast cancer risk, density, and stiffness: global tissue stiffness on breast MR elastography (MRE).

PURPOSE: Quantify in vivo biomechanical tissue properties in various breast densities and in average risk and high-risk women using Magnetic Resonance Imaging (MRI)/MRE and examine the association between breast biomechanical properties and cancer risk based on patient demographics and clinical data. METHODS: Patients with average risk or high-risk of breast cancer underwent 3.0 T breast MR imaging and elastography. Breast parenchymal enhancement (BPE), density (from most recent mammogram), stiffness, elasticity, and viscosity were recorded. Within each breast density group (non-dense versus dense), stiffness, elasticity, and viscosity were compared across risk groups (average versus high). Separately for stiffness, elasticity, and viscosity, a multivariable logistic regression model was used to evaluate whether the MRE parameter predicted risk status after controlling for clinical factors. RESULTS: 50 average risk and 86 high-risk patients were included. Risk groups were similar in age, density, and menopausal status. Among patients with dense breasts, mean stiffness, elasticity, and viscosity were significantly higher in high-risk patients (N = 55) compared to average risk patients (N = 34; all p < 0.001). Stiffness remained a significant predictor of risk status (OR = 4.26, 95% CI [1.96, 9.25]) even after controlling for breast density, BPE, age, and menopausal status. Similar results were seen for elasticity and viscosity. CONCLUSION: A structurally based, quantitative biomarker of tissue stiffness obtained from MRE is associated with differences in breast cancer risk in dense breasts. Tissue stiffness could provide a novel prognostic marker to help identify high-risk women with dense breasts who would benefit from increased surveillance and/or risk reduction measures.

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.

From cells to tissue: How cell scale heterogeneity impacts glioblastoma growth and treatment response.

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.

MetaMarker: a pipeline for de novo discovery of novel metagenomic biomarkers.

SUMMARY: We present MetaMarker, a pipeline for discovering metagenomic biomarkers from whole-metagenome sequencing samples. Different from existing methods, MetaMarker is based on a de novo approach that does not require mapping raw reads to a reference database. We applied MetaMarker on whole-metagenome sequencing of colorectal cancer (CRC) stool samples from France to discover CRC specific metagenomic biomarkers. We showed robustness of the discovered biomarkers by validating in independent samples from Hong Kong, Austria, Germany and Denmark. We further demonstrated these biomarkers could be used to build a machine learning classifier for CRC prediction. AVAILABILITY AND IMPLEMENTATION: MetaMarker is freely available at https://bitbucket.org/mkoohim/metamarker under GPLv3 license. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

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.

Speed Switch in Glioblastoma Growth Rate due to Enhanced Hypoxia-Induced Migration.

We analyze the wave speed of the Proliferation Invasion Hypoxia Necrosis Angiogenesis (PIHNA) model that was previously created and applied to simulate the growth and spread of glioblastoma (GBM), a particularly aggressive primary brain tumor. We extend the PIHNA model by allowing for different hypoxic and normoxic cell migration rates and study the impact of these differences on the wave-speed dynamics. Through this analysis, we find key variables that drive the outward growth of the simulated GBM. We find a minimum tumor wave-speed for the model; this depends on the migration and proliferation rates of the normoxic cells and is achieved under certain conditions on the migration rates of the normoxic and hypoxic cells. If the hypoxic cell migration rate is greater than the normoxic cell migration rate above a threshold, the wave speed increases above the predicted minimum. This increase in wave speed is explored through an eigenvalue and eigenvector analysis of the linearized PIHNA model, which yields an expression for this threshold. The PIHNA model suggests that an inherently faster-diffusing hypoxic cell population can drive the outward growth of a GBM as a whole, and that this effect is more prominent for faster-proliferating tumors that recover relatively slowly from a hypoxic phenotype. The findings presented here act as a first step in enabling patient-specific calibration of the PIHNA model.

A Mechanistic Investigation into Ischemia-Driven Distal Recurrence of Glioblastoma.

Glioblastoma (GBM) is the most aggressive primary brain tumor with a short median survival. Tumor recurrence is a clinical expectation of this disease and usually occurs along the resection cavity wall. However, previous clinical observations have suggested that in cases of ischemia following surgery, tumors are more likely to recur distally. Through the use of a previously established mechanistic model of GBM, the Proliferation Invasion Hypoxia Necrosis Angiogenesis (PIHNA) model, we explore the phenotypic drivers of this observed behavior. We have extended the PIHNA model to include a new nutrient-based vascular efficiency term that encodes the ability of local vasculature to provide nutrients to the simulated tumor. The extended model suggests sensitivity to a hypoxic microenvironment and the inherent migration and proliferation rates of the tumor cells are key factors that drive distal recurrence.

Learning Equations from Biological Data with Limited Time Samples.

Equation learning methods present a promising tool to aid scientists in the modeling process for biological data. Previous equation learning studies have demonstrated that these methods can infer models from rich datasets; however, the performance of these methods in the presence of common challenges from biological data has not been thoroughly explored. We present an equation learning methodology comprised of data denoising, equation learning, model selection and post-processing steps that infers a dynamical systems model from noisy spatiotemporal data. The performance of this methodology is thoroughly investigated in the face of several common challenges presented by biological data, namely, sparse data sampling, large noise levels, and heterogeneity between datasets. We find that this methodology can accurately infer the correct underlying equation and predict unobserved system dynamics from a small number of time samples when the data are sampled over a time interval exhibiting both linear and nonlinear dynamics. Our findings suggest that equation learning methods can be used for model discovery and selection in many areas of biology when an informative dataset is used. We focus on glioblastoma multiforme modeling as a case study in this work to highlight how these results are informative for data-driven modeling-based tumor invasion predictions.

A Deep Convolutional Neural Network for Annotation of Magnetic Resonance Imaging Sequence Type.

The explosion of medical imaging data along with the advent of big data analytics has launched an exciting era for clinical research. One factor affecting the ability to aggregate large medical image collections for research is the lack of infrastructure for automated data annotation. Among all imaging modalities, annotation of magnetic resonance (MR) images is particularly challenging due to the non-standard labeling of MR image types. In this work, we aimed to train a deep neural network to annotate MR image sequence type for scans of brain tumor patients. We focused on the four most common MR sequence types within neuroimaging: T1-weighted (T1W), T1-weighted post-gadolinium contrast (T1Gd), T2-weighted (T2W), and T2-weighted fluid-attenuated inversion recovery (FLAIR). Our repository contains images acquired using a variety of pulse sequences, sequence parameters, field strengths, and scanner manufacturers. Image selection was agnostic to patient demographics, diagnosis, and the presence of tumor in the imaging field of view. We used a total of 14,400 two-dimensional images, each visualizing a different part of the brain. Data was split into train, validation, and test sets (9600, 2400, and 2400 images, respectively) and sets consisted of equal-sized groups of image types. Overall, the model reached an accuracy of 99% on the test set. Our results showed excellent performance of deep learning techniques in predicting sequence types for brain tumor MR images. We conclude deep learning models can serve as tools to support clinical research and facilitate efficient database management.

Automatic 3D Nonlinear Registration of Mass Spectrometry Imaging and Magnetic Resonance Imaging Data.

Multimodal integration between mass spectrometry imaging (MSI) and radiology-established modalities such as magnetic resonance imaging (MRI) would allow the investigations of key questions in complex biological systems such as the central nervous system. Such integration would provide complementary multiscale data to bridge the gap between molecular and anatomical phenotypes, potentially revealing new insights into molecular mechanisms underlying anatomical pathologies presented on MRI. Automatic coregistration between 3D MSI/MRI is a computationally challenging process due to dimensional complexity, MSI data sparsity, lack of direct spatial-correspondences, and nonlinear tissue deformation. Here, we present a new computational approach based on stochastic neighbor embedding to nonlinearly align 3D MSI to MRI data, identify and reconstruct biologically relevant molecular patterns in 3D, and fuse the MSI datacube to the MRI space. We demonstrate our method using multimodal high-spectral resolution matrix-assisted laser desorption ionization (MALDI) 9.4 T MSI and 7 T in vivo MRI data, acquired from a patient-derived, xenograft mouse brain model of glioblastoma following administration of the EGFR inhibitor drug of Erlotinib. Results show the distribution of some identified molecular ions of the EGFR inhibitor erlotinib, a phosphatidylcholine lipid, and cholesterol, which were reconstructed in 3D and mapped to the MRI space. The registration quality was evaluated on two normal mouse brains using the Dice coefficient for the regions of brainstem, hippocampus, and cortex. The method is generic and can therefore be applied to hyperspectral images from different mass spectrometers and integrated with other established in vivo imaging modalities such as computed tomography (CT) and positron emission tomography (PET).

The 2019 mathematical oncology roadmap.

Whether the nom de guerre is Mathematical Oncology, Computational or Systems Biology, Theoretical Biology, Evolutionary Oncology, Bioinformatics, or simply Basic Science, there is no denying that mathematics continues to play an increasingly prominent role in cancer research. Mathematical Oncology-defined here simply as the use of mathematics in cancer research-complements and overlaps with a number of other fields that rely on mathematics as a core methodology. As a result, Mathematical Oncology has a broad scope, ranging from theoretical studies to clinical trials designed with mathematical models. This Roadmap differentiates Mathematical Oncology from related fields and demonstrates specific areas of focus within this unique field of research. The dominant theme of this Roadmap is the personalization of medicine through mathematics, modelling, and simulation. This is achieved through the use of patient-specific clinical data to: develop individualized screening strategies to detect cancer earlier; make predictions of response to therapy; design adaptive, patient-specific treatment plans to overcome therapy resistance; and establish domain-specific standards to share model predictions and to make models and simulations reproducible. The cover art for this Roadmap was chosen as an apt metaphor for the beautiful, strange, and evolving relationship between mathematics and cancer.

Lesion Dynamics Under Varying Paracrine PDGF Signaling in Brain Tissue.

Paracrine PDGF signaling is involved in many processes in the body, both normal and pathological, including embryonic development, angiogenesis, and wound healing as well as liver fibrosis, atherosclerosis, and cancers. We explored this seemingly dual (normal and pathological) role of PDGF mathematically by modeling the release of PDGF in brain tissue and then varying the dynamics of this release. Resulting simulations show that by varying the dynamics of a PDGF source, our model predicts three possible outcomes for PDGF-driven cellular recruitment and lesion growth: (1) localized, short duration of growth, (2) localized, chronic growth, and (3) widespread chronic growth. Further, our model predicts that the type of response is much more sensitive to the duration of PDGF exposure than the maximum level of that exposure. This suggests that extended duration of paracrine PDGF signal during otherwise normal processes could potentially lead to lesions having a phenotype consistent with pathologic conditions.

Simulating PDGF-Driven Glioma Growth and Invasion in an Anatomically Accurate Brain Domain.

Gliomas are the most common of all primary brain tumors. They are characterized by their diffuse infiltration of the brain tissue and are uniformly fatal, with glioblastoma being the most aggressive form of the disease. In recent years, the over-expression of platelet-derived growth factor (PDGF) has been shown to produce tumors in experimental rodent models that closely resemble this human disease, specifically the proneural subtype of glioblastoma. We have previously modeled this system, focusing on the key attribute of these experimental tumors-the "recruitment" of oligodendroglial progenitor cells (OPCs) to participate in tumor formation by PDGF-expressing retrovirally transduced cells-in one dimension, with spherical symmetry. However, it has been observed that these recruitable progenitor cells are not uniformly distributed throughout the brain and that tumor cells migrate at different rates depending on the material properties in different regions of the brain. Here we model the differential diffusion of PDGF-expressing and recruited cell populations via a system of partial differential equations with spatially variable diffusion coefficients and solve the equations in two spatial dimensions on a mouse brain atlas using a flux-differencing numerical approach. Simulations of our in silico model demonstrate qualitative agreement with the observed tumor distribution in the experimental animal system. Additionally, we show that while there are higher concentrations of OPCs in white matter, the level of recruitment of these plays little role in the appearance of "white matter disease," where the tumor shows a preponderance for white matter. Instead, simulations show that this is largely driven by the ratio of the diffusion rate in white matter as compared to gray. However, this ratio has less effect on the speed of tumor growth than does the degree of OPC recruitment in the tumor. It was observed that tumor simulations with greater degrees of recruitment grow faster and develop more nodular tumors than if there is no recruitment at all, similar to our prior results from implementing our model in one dimension. Combined, these results show that recruitment remains an important consideration in understanding and slowing glioma growth.

Mathematical model of perineural tumor spread: a pilot study.

BACKGROUND: Perineural spread (PNS) of pelvic cancer along the lumbosacral plexus is an emerging explanation for neoplastic lumbosacral plexopathy (nLSP) and an underestimated source of patient morbidity and mortality. Despite the increased incidence of PNS, these patients are often times a clinical conundrum-to diagnose and to treat. Building on previous results in modeling glioblastoma multiforme (GBM), we present a mathematical model for predicting the course and extent of the PNS of recurrent tumors. METHODS: We created three-dimensional models of perineurally spreading tumor along the lumbosacral plexus from consecutive magnetic resonance imaging scans of two patients (one each with prostate cancer and cervical cancer). We adapted and applied a previously reported mathematical model of GBM to progression of tumor growth along the nerves on an anatomical model obtained from a healthy subject. RESULTS: We were able to successfully model and visualize perineurally spreading pelvic cancer in two patients; average growth rates were 60.7 mm/year for subject 1 and 129 mm/year for subject 2. The model correlated well with extent of PNS on MRI scans at given time points. CONCLUSIONS: This is the first attempt to model perineural tumor spread and we believe that it provides a glimpse into the future of disease progression monitoring. Every tumor and every patient are different, and the possibility to report treatment response using a unified scale-as "days gained"-will be a necessity in the era of individualized medicine. We hope our work will serve as a springboard for future connections between mathematics and medicine.

Patient-specific mathematical neuro-oncology: using a simple proliferation and invasion tumor model to inform clinical practice.

Glioblastoma multiforme (GBM) is the most common malignant primary brain tumor associated with a poor median survival of 15-18 months, yet there is wide heterogeneity across and within patients. This heterogeneity has been the source of significant clinical challenges facing patients with GBM and has hampered the drive toward more precision or personalized medicine approaches to treating these challenging tumors. Over the last two decades, the field of Mathematical Neuro-oncology has grown out of desire to use (often patient-specific) mathematical modeling to better treat GBMs. Here, we will focus on a series of clinically relevant results using patient-specific mathematical modeling. The core model at the center of these results incorporates two hallmark features of GBM, proliferation [Formula: see text] and invasion (D), as key parameters. Based on routinely obtained magnetic resonance images, each patient's tumor can be characterized using these two parameters. The Proliferation-Invasion (PI) model uses [Formula: see text] and D to create patient-specific growth predictions. The PI model, its predictions, and parameters have been used in a number of ways to derive biological insight. Beyond predicting growth, the PI model has been utilized to identify patients who benefit from different surgery strategies, to prognosticate response to radiation therapy, to develop a treatment response metric, and to connect clinical imaging features and genetic information. Demonstration of the PI model's clinical relevance supports the growing role for it and other mathematical models in routine clinical practice.

Biologically-informed deep neural networks provide quantitative assessment of intratumoral heterogeneity in post-treatment glioblastoma.

Intratumoral heterogeneity poses a significant challenge to the diagnosis and treatment of glioblastoma (GBM). This heterogeneity is further exacerbated during GBM recurrence, as treatment-induced reactive changes produce additional intratumoral heterogeneity that is ambiguous to differentiate on clinical imaging. There is an urgent need to develop non-invasive approaches to map the heterogeneous landscape of histopathological alterations throughout the entire lesion for each patient. We propose to predictively fuse Magnetic Resonance Imaging (MRI) with the underlying intratumoral heterogeneity in recurrent GBM using machine learning (ML) by leveraging image-localized biopsies with their associated locoregional MRI features. To this end, we develop BioNet, a biologically-informed neural network model, to predict regional distributions of three tissue-specific gene modules: proliferating tumor, reactive/inflammatory cells, and infiltrated brain tissue. BioNet offers valuable insights into the integration of multiple implicit and qualitative biological domain knowledge, which are challenging to describe in mathematical formulations. BioNet performs significantly better than a range of existing methods on cross-validation and blind test datasets. Voxel-level prediction maps of the gene modules by BioNet help reveal intratumoral heterogeneity, which can improve surgical targeting of confirmatory biopsies and evaluation of neuro-oncological treatment effectiveness. The non-invasive nature of the approach can potentially facilitate regular monitoring of the gene modules over time, and making timely therapeutic adjustment. These results also highlight the emerging role of ML in precision medicine.

Androgen loss weakens anti-tumor immunity and accelerates brain tumor growth.

Many cancers, including glioblastoma (GBM), have a male-biased sex difference in incidence and outcome. The underlying reasons for this sex bias are unclear but likely involve differences in tumor cell state and immune response. This effect is further amplified by sex hormones, including androgens, which have been shown to inhibit anti-tumor T cell immunity. Here, we show that androgens drive anti-tumor immunity in brain tumors, in contrast to its effect in other tumor types. Upon castration, tumor growth was accelerated with attenuated T cell function in GBM and brain tumor models, but the opposite was observed when tumors were located outside the brain. Activity of the hypothalamus-pituitary-adrenal gland (HPA) axis was increased in castrated mice, particularly in those with brain tumors. Blockade of glucocorticoid receptors reversed the accelerated tumor growth in castrated mice, indicating that the effect of castration was mediated by elevated glucocorticoid signaling. Furthermore, this mechanism was not GBM specific, but brain specific, as hyperactivation of the HPA axis was observed with intracranial implantation of non-GBM tumors in the brain. Together, our findings establish that brain tumors drive distinct endocrine-mediated mechanisms in the androgen-deprived setting and highlight the importance of organ-specific effects on anti-tumor immunity.

Biologically informed deep neural networks provide quantitative assessment of intratumoral heterogeneity in post treatment glioblastoma.

Intratumoral heterogeneity poses a significant challenge to the diagnosis and treatment of recurrent glioblastoma. This study addresses the need for non-invasive approaches to map heterogeneous landscape of histopathological alterations throughout the entire lesion for each patient. We developed BioNet, a biologically-informed neural network, to predict regional distributions of two primary tissue-specific gene modules: proliferating tumor (Pro) and reactive/inflammatory cells (Inf). BioNet significantly outperforms existing methods (p < 2e-26). In cross-validation, BioNet achieved AUCs of 0.80 (Pro) and 0.81 (Inf), with accuracies of 80% and 75%, respectively. In blind tests, BioNet achieved AUCs of 0.80 (Pro) and 0.76 (Inf), with accuracies of 81% and 74%. Competing methods had AUCs lower or around 0.6 and accuracies lower or around 70%. BioNet's voxel-level prediction maps reveal intratumoral heterogeneity, potentially improving biopsy targeting and treatment evaluation. This non-invasive approach facilitates regular monitoring and timely therapeutic adjustments, highlighting the role of ML in precision medicine.

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.