Jackson, Pamela

Abstract Glioma stem cells (GSCs) are considered a major driver of glioblastoma (GBM) progression and are highly resistant to standard cytotoxic treatments. BMP4 has been shown to promote GSC differentiation, enhance radiosensitivity, slow tumor growth, and extend survival in animal models. Despite this promise, BMP4 has yet to achieve clinical impact, owing largely to heterogeneous and nonlinear responses across preclinical experimental systems. To elucidate how BMP4 could function as an effective targeted differentiation therapy in GBM, we develop a mathematical model that describes the growth of a GBM tumor via a hierarchy of GSCs, progenitor and terminally differentiated cells. We parameterize our model using new radiotherapy and proliferation assay data (with and without BMP4 exposure) from twelve patient-derived GSC lines. This integration of model and data allows us, for the first time, to quantitatively capture patient-specific heterogeneity in BMP4 sensitivity. We perform global sensitivity analysis on the model, identifying proliferation rate and GSC self-renewal sensitivity as key determinants of BMP4 efficacy. These parameters act as model-derived biomarkers that distinguish BMP4-responsive tumors from non-responsive ones. In an in silico analysis across a broad cohort of virtual patients, we find that continuous BMP4 delivery from surgical resection through radiotherapy consistently outperforms a single-dose strategy. Virtual clinical trials further show that, without stratification using these model-derived biomarkers, BMP4 yields little observable therapeutic benefit. In contrast, selecting patients with more proliferative, BMP4-responsive GSCs markedly increases the likelihood of observing a significant treatment effect. Citation Format: Nicholas Harbour, Lee Curtin, Loizos Michaelides, Matthew E. Hubbard, Pamela Jackson, Vinitha Rani, Rajappa Kenchappa, Virginea Farias, Anna Carrano, Markus Owen, Alfredo Quinones-Hinojosa, Kristin Swanson. Virtual clinical trials of BMP4 differentiation therapy [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2026; Part 1 (Regular Abstracts); 2026 Apr 17-22; San Diego, CA. Philadelphia (PA): AACR; Cancer Res 2026;86(7 Suppl):Abstract nr 6831.

Abstract High-grade glioma (HGG) continues to have a dismal prognosis owing in part to its heterogeneous and invasive nature both within and between patients. The most aggressive form of high-grade glioma, called glioblastoma, portends a median survival of around 15 months. Like many cancers, high-grade glioma rewires metabolism towards its own benefit, shifting away from the typical reliance of brain tissue on glucose towards other sources such as fatty acids and acetate. Using a cohort of image-localized high-grade glioma biopsies collected from diverse imageable regions with associated transcriptomics (58 patients, 202 biopsies), we are able to assess the spatial transcriptomic landscape of HGG. We have previously shown that the transcriptome of HGG can be characterized into a continuous transition between states, and we have shown that the cellular ecosystem follows a similar pattern, using a combination of trajectory inference and cellular deconvolution. These patterns are strongly associated with the non-enhancing and contrast-enhancing regions of the tumors, which are typically considered for treatment planning and progression assessment. Here, we applied differential expression of hallmark pathways to determine how metabolism factors into this paradigm. As expected, we find overall higher levels of metabolism-associated pathway expression in the contrast enhancing regions compared to the non-enhancing regions. Interestingly, we found metabolic expression differs between the ecologies associated with transcriptional subtypes of high-grade glioma. For example, fatty acid (p=0.003), xenobiotic (p<0.0001) and heme (p<0.0001) metabolism pathways were significantly upregulated for mesenchymal tumor samples, while significantly downregulated for proliferative tumor samples (p=0.02, p<0.0001, p<0.0001, respectively). As we continue to further understand the complex spatial and temporally heterogeneous landscape of these aggressive tumors, it is important to consider the metabolites on which they depend, to better inform prognosis and develop new treatments.

Abstract High-grade gliomas (HGGs) are aggressive, infiltrative brain tumors with poor survival outcomes and limited therapeutic options. Their spatial and molecular heterogeneity necessitates large, well-curated datasets integrating imaging, clinical, and genomic information to enable precision medicine and drive novel treatment development. Using the Medtronic StealthStation neuronavigation system, we collected over 1,000 spatially localized intraoperative biopsy samples from more than 250 patients undergoing resection for glioma at Mayo Clinic and partnering sites. Each biopsy site was recorded with 3D coordinates, labeled in real time, flash-frozen, and annotated with metadata including collection time, anatomical location, and surgical plan. Samples were processed and entered into a secure database along with the corresponding sampling coordinates, frozen pathology, whole-exome and RNA sequencing data, 5-ALA fluorescence status, patient demographics, and T1, T1Gd, T2, and T2-FLAIR MRI. This multimodal dataset supports radiogenomic research by enabling precise mapping of imaging features to molecular characteristics. A subset of this data, including 268 biopsies from 52 patients (25 male, 27 female) with reported intraoperative 5-ALA fluorescence, was used to develop a radiomics model predicting locoregional patterns of 5-ALA positivity based on preoperative MRI. Image features were extracted from the biopsy sample coordinates across multiple MRI sequences and paired with corresponding 5-ALA fluorescence status to train the model for regional fluorescence prediction. The final model achieved 83% accuracy on the training set and 86% on the validation set using a 70-15-15 data split; this included 208 samples for training and 42 for validation, stratified by patient to ensure no overlap of samples from the same individual across sets. These findings demonstrate the feasibility of predicting 5-ALA fluorescence using radiomic features from standard preoperative MRI, offering a promising noninvasive strategy for improving intraoperative decision-making and enabling more precise, personalized resections in patients with high-grade gliomas.

Abstract High grade gliomas (HGGs) are aggressive brain tumors that are difficult to remove completely, leading to poor survival rates. Contrast enhancement on T1Gd magnetic resonance imaging (MRI) has traditionally been used to guide resection, however it often underestimates true tumor burden, adding a layer of uncertainty for surgeons during pre-operative planning. Fluorescence guided surgery using 5-aminolevulinic acid (5-ALA) has emerged as an intraoperative solution and understanding it in conjunction with MRI features can address the limitations of routine T1Gd imaging during surgical planning. Here we investigate sex-specific patterns in HGGs using 5-ALA fluorescence and imaging features. We analyzed a patient dataset of 202 multiregional MRI-localized biopsies from 42 HGG patients (20 male and 22 female). Data included patient demographics, 5-ALA status (positive or negative), MRI annotations (CE = contrast enhancement, NE = non-enhancement, BAT = brain around tumor non-enhancement), and survival outcomes. Statistical analyses included Wilcoxon rank-sum tests, Kruskal-Wallis tests, and chi-squared tests to assess differences in imaging annotations and 5-ALA positive vs negative status among sex cohorts. Among all samples, CE and BAT were associated with higher 5-ALA positivity, while NE was associated with higher 5-ALA negativity (p = 0.02119). Female patients were less likely to have detectable 5-ALA in their multiregional biopsy samples than males (p < 0.0001). CE at the biopsy location was more likely to be indicative of 5-ALA positivity in males than females (p = 0.001927); whereas NE was associated with higher 5-ALA negativity in females than males (p = 0.03426). BAT was associated with higher 5-ALA positivity in males than females (p = 0.03895). Our study reveals sex-specific patterns in 5-ALA fluorescence and image features among HGG patients. These unique features may serve as potential pre-surgical markers for tumor regions that are more likely to benefit from 5-ALA-guided resection.

Abstract The intra- and inter-patient heterogeneity in high-grade glioma (HGG) continues to contribute to its poor prognosis. Clinical biopsies are often harvested from limited regions and typically are not image localized. Thus, they fail to capture the diversity within tumor regions, immune expression or normal cell abundances that play key roles in tumor development. It is important to gain an understanding of these subpopulation ecologies, their spatial resolution, and interactions between them that may then be exploited for future therapeutic benefit. Further, these may differ by patient characteristics such as sex, age at diagnosis and treatment status. Using an ongoing image-localized biopsy collection protocol, we have so far evaluated the bulk transcriptomics of 202 multi-regional biopsies from 58 patients to characterize HGG heterogeneity. These samples were processed through Monocle, a reverse graph embedding algorithm that groups samples into states and orders them along developmental trajectories. Deconvolution methods were previously used to predict relative abundances of 7 normal, 6 glioma, and 5 immune cell subpopulations for each sample, which we have now overlaid on the Monocle graph. Monocle classified HGG into 4 main states along a three-pronged trajectory. These states reveal distinct population ecologies with associated enriched gene pathways. We also note significant immune pathway enrichments that differ between state and patient-reported sex. Together, these algorithms reveal a simple transcriptomic trajectory that helps us understand the development and evolution of HGG. Characterizing the in vivo diversity within and between high grade gliomas is important for understanding prognosis, stratifying future treatments and ultimately improving patient outcome.

Abstract Magnetic resonance imaging (MRI) is key to clinically managing brain tumor patients, however connecting the biology to imaging remains challenging. We previously developed a two-compartment model of MRI signal intensity to quantitatively estimate relative edema abundance from T2-weighted MRIs. Using this model, we previously identified the fatty acid metabolism (FAM) and oxidative phosphorylation (OxPhos) pathways as sex-distinct, with both pathways amplified for high edema in males and low edema in females. The purpose of this project was to further delineate sex-distinct biology associated with MRI-estimated brain tumor edema abundance. We analyzed 179 multiregional samples (Female: 75; Male: 104) from 55 high grade glioma patients (Female: 21; Male: 34) for bulk RNA-Seq. Patients’ pre-surgical multiparametric MRIs were preprocessed and segmented for abnormal regions and normal tissue. Utilizing the segmentations and preprocessed images we estimated the relative fractions of extracellular and intracellular space based on the edema mathematical model. Samples were characterized by their edema scores, which were analyzed using differential expression for high and low edema, gene set enrichment analysis (GSEA) using MSigDB hallmarks, and leading edge interpretation. Based on transcriptomic leading edge analyses of high and low edema samples from sex-separated cohorts, we identified common and sex-distinct enrichment of the FAM and OxPhos pathways underlying edema patterns. Of the OxPhos pathway leading edge genes, 45 were common, 36 were unique to females, and 66 were unique to males. Of the FAM pathway leading edge genes, 8 were common, 39 were unique to females, and 29 were unique to males. Notably, expression of both IDH1 and IDH2 were increased for males in regions of high edema in the OxPhos pathway. IDH3a was decreased for females in regions of low edema in the OxPhos pathway. These data suggest that there may be sex-distinct metabolism underlying MRI measurable edema formation.

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

Abstract Magnetic resonance imaging (MRI) is a centerpiece of clinical management of brain tumor patients, yet the biology underlying even commonly interpreted imaging changes is unclear. The purpose of this project was to understand what biology is associated with brain tumor edema abundance and if that association was sex distinct. We have developed a two-compartment model of MRI signal intensity to quantitatively estimate relative edema abundance from T2-weighted MRIs. Multiparametric MRIs were preprocessed with bias field correction and intensity normalization. The images were co-registered and segmented for abnormal regions and normal tissue. The segmentations and preprocessed images were utilized in the edema mathematical model to estimate the fraction of extracellular space (Fecs) and fraction of intracellular space (Fics). Through an ongoing image-localized biopsy collection protocol, we analyzed 179 samples (Female: 75; Male: 104) from 55 patients (Female: 21; Male: 34) for bulk RNA-Seq. Samples were characterized by their edema Fecs scores and underwent differential expression, gene set enrichment analysis using MSigDB hallmarks and cellular deconvolution using CIBERSORTx. Through differential expression and gene set enrichment analyses of these spatial biopsies, we have identified sex-distinct gene expression and pathways corresponding to regions of high vs low edema. Within highly edematous regions, male samples were amplified for respiratory pathways, while female samples were suppressed. Connecting changes on MRIs with molecular markers from image-localized biopsies provides an opportunity to identify the biological drivers of brain tumor associated imageable edema. Specifically, elucidating the sex-distinct patterns connecting changes on MRI with cellular populations and molecular pathways could be used to better interpret imaging of patients and interpreting treatment response. Citation Format: Pamela R. Jackson, Lee Curtin, Sara Ranjbar, Kyle W. Singleton, Maciej M. Mrugala, Richard S. Zimmerman, Bernard R. Bendok, Peter Canoll, Kristin R. Swanson. Sex-distinct patterns of molecular pathways associated with brain tumor edema [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2024; Part 1 (Regular Abstracts); 2024 Apr 5-10; San Diego, CA. Philadelphia (PA): AACR; Cancer Res 2024;84(6_Suppl):Abstract nr 6168.

Abstract High-grade glioma continues to have dismal survival owing in part to its intra- and inter-patient heterogeneity. Standard clinical protocol collects tumor samples with the aim of providing or confirming a diagnosis and determining the status of a few key genes (e.g. IDH1, MGMT). However, this protocol is unable to capture the diversity within tumor regions, immune expression or normal cell abundances that play key roles in the development of the disease. To overcome this, during surgery we collect image-localized multi-regional biopsies to characterize this disease heterogeneity. Data collection is ongoing, and we currently have 202 samples from 58 patients with available bulk RNA-Seq. With a single-cell reference dataset from Columbia University, we used CIBERSORTx, a deconvolution method, to predict relative abundances of 7 normal, 6 glioma, and 5 immune cell subpopulations for each sample. We used Cox proportional hazard models with first-order statistics (mean, minimum, maximum) of abundances within patients to determine whether these are prognostic, then used TCGA RNA-Seq data to validate these findings. We found that one glioma proneural subtype was significantly beneficial for patient survival relative to other glioma subtypes across all statistics (and showed significance in TCGA). Proliferative and mesenchymal glioma subtypes also showed significance for one or two of the calculated statistics. Oligodendrocyte progenitor cell abundances were consistently significantly beneficial for patient survival in our cohort, while an increased abundance of abnormal (reactive) astrocytes is associated with poor prognosis. We also ran these analyses within invasive margin and core tumor regions to determine location-specific population abundances associated with patient survival. In conclusion, understanding the in vivo diversity of cellular subpopulations within high grade glioma is important for treatment stratification and patient benefit.

Abstract Identification of key phenotypic regions such as necrosis, contrast enhancement, and edema on magnetic resonance imaging (MRI) is important for understanding disease evolution and treatment response in patients with glioma. Manual delineation is time intensive and not feasible for a clinical workflow. Automating phenotypic region segmentation overcomes many issues with manual segmentation, however, current glioma segmentation datasets focus on pre-treatment, diagnostic scans, where treat-ment effects and surgical cavities are not present. Thus, existing automatic segmenta-tion models are not applicable to post-treatment imaging that is used for longitudinal evaluation of care. Here, we present a comparison of three-dimensional convolutional neural networks (nnU-Net architecture) trained on large temporally defined pre-treatment, post-treatment, and mixed cohorts. We used a total of 1563 imaging timepoints from 854 patients curated from 13 different institutions as well as diverse public data sets to understand the capabilities and limitations of automatic segmenta-tion on glioma images with different phenotypic and treatment appearance. We as-sessed the performance of models using Dice coefficients on test cases from each group comparing predictions with manual segmentations generated by trained techni-cians. We demonstrate that training a combined model can be as effective as models trained on just one temporal group. The results highlight the importance of a diverse training set, that includes images from the course of disease and with effects from treatment, in the creation of a model that can accurately segment glioma MRIs at multiple treatment time points.