Abstract
Glioblastomas (GBMs) are known to be complex tumors comprising multiple subpopulations of genetically-distinct cancer cells; it is thought that this genetic variation is a major factor in the lack of observed survival benefit of treatment regimes that target one of these subpopulations. The field of radiogenomics seeks to study correlations between MRI patterns and genetic features of GBM tumors. Spatial radiogenomic maps produced using machine-learning (ML) methods that are trained against information from image-localized patient biopsies identify regions where particular cancer sub-populations are predicted to occur within a GBM, thus non-invasively characterizing the regional genetic variability of these tumors. These tumor subpopulations may also interact with one another, in ways which may be of a competitive or cooperative nature to varying degrees. It is important to ascertain the nature of these interactions, as they may have implications for treatment response to targeted therapies, and characterization of the spatio-temporal dynamics of these co-evolving sub-populations will shed light on why some therapies fail. Here we combine mathematical modeling techniques and spatially-resolved radiogenomic maps to study the nature of these interactions between molecularly-distinct GBM subpopulations. We model the interactions between cell populations using a partial differential equation based formalism. The model is parameterized using radiogenomic ML maps from which we infer the nature of interactions between subpopulations. Furthermore, using maps as inputs, the model turns static maps into dynamic information, thus providing insight into how these subpopulations composing the tumor change over time and the effect this has on observed treatment response for individual patients.