Citations

Abstract Glioblastoma (GBM) is a lethal, treatment-resistant brain cancer characterized by diffuse brain invasion, complex neuron-glioma signaling, and cellular, molecular, and spatial heterogeneity. Surgery is the current mainstay of treatment, aimed at maximizing tumor cytoreduction, decompressing the brain to improve neurological function, and providing material for analyses to inform prognostication and adjuvant treatment selection. To accurately model GBM surgical resection and tissue sampling in preclinical animal models, we developed the Maximal Safe Intracranial Surgery (MSIS) system. MSIS integrated a miniaturized neurosurgical resection device and enabled a safe, tunable, and reproducible approach for intracranial tumor resection, multi-regional biopsy, and longitudinal tissue sampling in murine models, with stereotactic, sterile, and high-viability tissue collection. MSIS was applied across five distinct orthotopic brain tumor models to reproducibly generate histologically and radiologically representative models of post-surgical GBM. Similar to humans, increasing the extent of resection correlated with improved animal survival, and recurrent growth patterns varied significantly by model type. Transcriptional profiles of matched primary vs. post-surgical recurrent tumors revealed enrichment of margin-associated transcriptional states consistent with neuron-glioma signaling. Overall, the MSIS system offers opportunities to model the standard-of-care, dissect the spatial and temporal dynamics of GBM, and pair therapeutic discovery with contextual tumor pathobiology.

Pharmacological plasminogen reduction enhanced liver regeneration experimentally and clinically. siRNA-induced plasminogen deficiency promoted hepatocyte proliferation after partial hepatectomy in mice, contrasting genetic deficiency models. In the HeLiX trial, tranexamic acid reduced post-hepatectomy liver failure odds, suggesting a novel therapeutic strategy.

Abstract Background Glioblastoma (GBM) is driven, in part, by glioma stem cells (GSCs), which sustain therapeutic resistance and recurrence. OLIG2 is a key regulator of GSC programs but the signaling pathways that cooperate with Oligodendrocyte transcription factor 2 (OLIG2) to maintain GSC states remain incompletely defined. Although Signal transducer and activator of transcription (STAT) signaling has been implicated in GSC biology, work to date has focused primarily on STAT3, and the contribution of STAT5 remains poorly defined. Here, we investigated whether STAT5B functionally cooperates with OLIG2 in sustaining GSC phenotypes. Methods Transcriptomic analyses, proximity ligation assays and biochemical studies were utilized to define OLIG2-STAT5B interactions in patient-derived GSCs. Functional consequences of OLIG2, STAT5B, BIRC5, and NUSAP1 suppression were assessed using small interfering RNA/short hairpin RNA-knockdown, stemness, and migration assays in EGFRvIII-expressing GSCs. In vivo tumor growth studies were performed using orthotopic xenograft models. Results OLIG2 and STAT5B expression were positively correlated in GBM datasets and patient-derived models. OLIG2 physically associated with STAT5B in EGFRvIII-expressing GSCs, and loss of OLIG2 reduced STAT5 activation, supporting a functional link between these factors. OLIG2 or STAT5B suppression induced convergent transcriptional changes, including downregulation of genes associated with GSC maintenance and proliferation, notably stemness regulators BIRC5 and NUSAP1. Silencing BIRC5 or NUSAP1 impaired GSC self-renewal and migration. In orthotopic models, STAT5B knockdown prolonged survival and phenocopied OLIG2 suppression and stemness loss. Conclusions These findings identify STAT5B as a functional partner of OLIG2 in GBM and define anOLIG2-STAT5B transcriptional axis that sustains GSC maintenance and malignant behavior. Targeting this axis may represent a therapeutic strategy to disrupt treatment-refractory glioma cell states.

Abstract Glioblastoma (GBM) is the most common and aggressive primary brain tumor in adults, with a dismal 5-year survival rate of ~6.5%. Vasogenic cerebral edema is a severe complication that occurs in the majority of GBM patients, leading to a significant increase in intracranial pressure, neurological deficits, and increased mortality rates. GBM-associated edema is almost exclusively managed by the corticosteroid dexamethasone, the use of which is associated with immunosuppression and interference with radiation therapy. While there is considerable research on edema formation, little is known about how edema resolution occurs in GBM, particularly in response to dexamethasone or other anti-edema therapies. This gap in understanding highlights the need for new approaches to treat edema effectively. The water channel aquaporin-4 (AQP4), a key player in the brain’s water homeostasis, is upregulated in GBM compared to the normal brain. To investigate the role of AQP4 in GBM-associated edema, we used the RCAS/tv-a system to generate de novo GBM with different driver mutations in Aqp4 knockout (KO) and wild-type (WT) mice. Our results indicated that Aqp4 KO tumor-bearing mice had shortened survival and increased edema content. Additionally, Aqp4 loss decreased astrocyte coverage specifically in the tumor core but not the tumor periphery. Our single cell RNA-seq data suggested that Aqp4 loss led to an astrocytic shift towards a less reactive, progenitor-like state, while also decreasing the myeloid and CD8+ T cell population in the tumor microenvironment. Interestingly, tumor-bearing mice treated with anti-VEGFA antibody, which has anti-edema properties, showed increased expression of Aqp4 and reduced vessel leakage. This suggests that AQP4 plays a critical role in mediating the effects of anti-VEGFA therapy. Overall, this study highlights the important role of AQP4 in the resolution of edema, while also revealing its additional effects on astrocyte activation and the immune microenvironment in GBM.

Abstract High-grade gliomas (HGG) have significant spatial heterogeneity, complicating therapeutic approaches. While multiparametric MRI (mpMRI) enables noninvasive assessment of this heterogeneity, most predictive models lack biological interpretability when used to classify or predict outcomes. Here, we developed an unsupervised framework that identifies novel mpMRI feature patterns through dimensional reduction, linking them to clinically relevant molecular signatures. We analyzed 18 advanced mpMRI features from 73 stereotactic biopsies across 20 HGG patients. Imaging data were normalized, followed by dimensional reduction using an iterative UMAP across varying parameters. Density-based clustering was used to generate stable image-derived clusters through consensus co-clustering. RNA-seq data from matched biopsies were analyzed using DESeq2 and GSEA to identify differentially expressed genes and enriched pathways. We then constructed a network linking specific imaging features to these transcriptional signatures. Validation was performed in an independent cohort, allowing us to extend clusters across full brain mpMRI. Unsupervised clustering of mpMRI features produced nine stable image-based clusters. These clusters were not fully explained by contrast enhancement, metabolic class, or tumor purity, shown by intraclass correlation (ICC) analyses, suggesting they capture novel biological variation. Transcriptomic analysis showed each cluster was associated with distinct gene expression patterns and enriched biological pathways, including neuronal differentiation, immune chemotaxis, and chromatin remodeling. Key imaging features (sRCBV, MTT, K2)—combined in cluster-specific patterns—enabled construction of an interpretable imaging–transcriptomic network. This was validated in an independent cohort. We then applied to full brain mpMRI, mapping spatially distinct regions to biological trends and patient outcomes. Addressing a research gap in imaging HGG, we present an interpretable network that connects specific patterns of mpMRI to clinically meaningful molecular signatures. This approach advances our understanding of tumor heterogeneity and supports future translational work in precision patient care.

Abstract BACKGROUND Glioblastoma (GBM) is a molecularly and phenotypically heterogeneous tumor, with regional invasion patterns driving treatment resistance and universal recurrence. In infiltrative regions, GBM cells adopt at least two predominant invasive phenotypes: Neuronal-enriched (NEU), with synaptic and neurodevelopment programs, and Glycolytic/Plurimetabolic (GPM), with heightened immune and metabolic activity. These subtypes differ in therapeutic vulnerability and microenvironmental interactions; features not captured by conventional contrast-enhanced (T1+C) MRI due to minimal blood-brain barrier (BBB) disruption. Our group sought to develop predictive models to distinguish variability of NEU/GPM subtypes within infiltrative GBM, utilizing graph convolutional networks (GCN), advanced multiparametric MRI (mpMRI), and spatially-matched image-localized biopsies with transcriptomic profiling. METHODS We analyzed 18 mpMRI features —including conventional, diffusion-weighted, and Dynamic Susceptibility Contrast (DSC) perfusion MRI— of 101 image-localized biopsies from MRI non-enhancing infiltrative margins of 64 GBM patient tumors. Features were hierarchically clustered per modality; each cluster was represented as a node in a Neo4j graph. Edges connected nodes sharing overlapping samples. Nodes were annotated with NEU or GPM subtype proportion. GCNs were trained to classify subtype-enriched nodes using graph topology and node features. RESULTS GCNs achieved high predictive performance on held-out test data (AUC = 0.95 for NEU, 0.88 for GPM; 70:30 train-test split) with similar 5-fold cross-validation performance and generalized well to an independent cohort (AUC > 0.85, n = 45). Importantly, performance was robust without relying on T1+C imaging, enabling prediction beyond BBB-disrupted regions. GCNs outperformed traditional classifiers (logistic regression, random forest). CONCLUSION Our results demonstrate that GCN-modeling of mpMRI features enables non-invasive prediction of invasive GBM subtypes at the tumor margin. Unlike convolutional neural networks (CNN), which rely on spatial continuity, GCNs leverage relational structure between heterogeneous imaging features, making them ideal for modeling sparse, non-grid data. This may support non-invasive regional vulnerability detection, surgical targeting, and improved response assessment.

Abstract Glioma stem cells (GSCs) are key drivers of glioblastoma (GBM) progression, therapy resistance, and recurrence, making them critical therapeutic targets. While the lineage-defining transcription factor OLIG2 is critical for GSC maintenance, the molecular co-regulators and downstream pathways it engages remain poorly defined. Here, we identify a previously unrecognized, functionally significant interaction between OLIG2 and STAT5B, uncovering a novel regulatory axis in GSCs. Using proximity ligation assay (PLA), we demonstrate that OLIG2 interacts with STAT5 in human and mouse GSC models, particularly in EGFRvIII-mutant cells. We find that OLIG2 preferentially partners with STAT5B and not STAT5A. Knockdown studies in GSCs reveal that OLIG2 and STAT5B co-regulate critical oncogenes BIRC5 and NUSAP1, both of which are required for sphere formation and migration capacity of EGFRvIII positive GSCs. In orthotopic xenograft models, knockdown of either STAT5B or OLIG2 significantly prolongs survival reduces expression of stemness and invasion markers, SOX2, VIMENTIN, and NUSAP1. In contrast, STAT5A knockdown fails to recapitulate these phenotypes, underscoring the role of OLIG2-STAT5B signaling in GSCs. These findings reveal a novel OLIG2-STAT5B transcriptional circuit that plays a role in GSC stemness maintenance, invasion, and tumor growth, particularly in EGFRvIII-driven gliomas. Targeting this axis may offer a new avenue for GBM treatment.

Abstract IDH-wildtype glioblastoma (GBM) is the most malignant primary adult brain tumor, with a five-year survival rate below 5%. Despite maximal treatment, including surgical resection, radiotherapy, and temozolomide (TMZ), recurrence is universal, driven by profound intratumoral heterogeneity, diffuse invasion, and therapy resistance. A hallmark of IDH-wildtype GBM is frequent alteration of epidermal growth factor receptor (EGFR), including amplification (50%), the EGFRvIII deletion variant (25%), and extracellular missense mutations such as G598V and A289V (10–15%). These missense mutations increase ligand affinity and receptor activity, yet their downstream signaling consequences and role in shaping the tumor phenotypes remain poorly defined. Recent pathway-based subclassifications identified four subtypes in IDH-wildtype GBM. Notably, neuronal (NEU) and glycolytic/plurimetabolic (GPM) subtypes predominate the invasive tumor margins. NEU programs are enriched in synaptic and neurodifferentiation signals, while GPM programs are associated with metabolic adaptation, angiogenesis, and immune evasion. Here, we show EGFR-G598V and EGFR-A289V mutations differentially regulate these invasive subtypes identities. Transcriptional analysis reveals EGFR G598V upregulates metabolic and inflammatory pathways consistent with the GPM phenotype, while EGFR A289V activates genes involved in neuronal differentiation and synaptic function, aligning with the NEU subtype. Phosphoproteomics and kinase-substrate analysis further uncovered distinct kinases to be functionally activated/inhibited in EGFR A289V-expressing cells and G598V cells. Mechanistic studies revealed EGFR A289V induces sustained EGFR phosphorylation at Y992, activating PLCγ2 via the SRC-family kinase LCK. This cascade drives calcium release, CAMK1 activation, and STAT5 phosphorylation, defining a novel EGFR A289V–LCK/PLCγ2/CAMK1/STAT5 axis. In addition, Y1068 phosphorylation is also constitutively active in EGFR A289V, resulting in increased BRAF activation. In contrast, EGFR G598V activates MAPK6 and CAMK2 inducing AKT2 activation, reinforcing GPM-associated signaling. These findings uncover mutation-specific signaling networks that drive subtype fate in invasive GBM and highlight unique pathway(s) as tractable targets for subtype-directed therapy in EGFR-mutant GBM.

Abstract Introduction: High-grade gliomas (HGG) are aggressive brain tumors with heterogeneous cell populations, including invasive regions driving disease progression and recurrence. Our previous work integrated molecular profiling with spatially matched advanced multiparametric MRI in a cohort of HGG patients with multi-regional sampling, uncovering distinct biological signatures linked to individual advanced MRI phenotypes. Yet, the relationships between MRI phenotypes—and how this potentially informs underlying biology—remains unexplored. Here, we leverage an unsupervised approach to dimensional reduction for a broad array of MRI features and compare regional transcriptomic signatures in a cohort of image-localized HGG samples. Methods: From our existing data, we selected 20 HGG patients with 73 multi-regional tumor biopsies, each including 18 advanced MRI parameter calculations. With an iterative approach, data was projected using uniform manifold approximation (UMAP) at multiple resolutions and clustered, producing high-confidence stable clusters. We used random forest modeling to assess cluster associations with previously identified signatures. Cluster-specific transcriptomic profiles were detected by DESeq2 and subsequent gene-set enrichment analysis. We further evaluated Pearson correlations and statistical fit with MRI image parameters by performing pathway analysis on genes with significant correlations by image feature, across clusters. Results: 8 clusters emerged from our analyses. Most advanced MRI features exhibited statistically significant differences across clusters, forming unique patterns of feature behavior. Notably, while these image-based clusters did not correlate well with metabolic profiles, contrast enhancement, or genomics, they did demonstrate distinct and HGG-relevant transcriptomic profiles, with shared and divergent signatures. Across clusters, we also identified several correlations between specific MRI parameters and gene pathway trends. These highlighted combined patterns of multiple MRI features, offering a spectrum of regional phenotypes imperceptible with single feature analysis. We generalize these results to an image-based decision matrix to detect biologic signatures for validation. Conclusion: Our integrative analysis of MRI mapping stratifies imaged-based patterns to non-invasively assess molecular signatures within HGG. This approach will help elucidate MRI biomarkers that may predict glioma behavior and support future individualized treatment approaches. Citation Format: Matthew W. Kenaston, Matthew Flick, Kasthuri Kannan, Michele Ceccarelli, Nhan Tran, Leland S. Hu. Integrative analysis of regional MRI phenotypes in high-grade glioma: Dimensional reduction of multi-omics imaging and molecular profiles [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2025; Part 1 (Regular Abstracts); 2025 Apr 25-30; Chicago, IL. Philadelphia (PA): AACR; Cancer Res 2025;85(8_Suppl_1):Abstract nr 5776.

Abstract BACKGROUND Standard of care surgery for high grade glioma (HGG) employs Magnetic Resonance Imaging to identify the contrast enhancing (CE) core for resection, leaving a non-enhancing (NE) invasive rim of tumor which contributes to recurrence. Multi-parametric MRI (mpMRI) has recently been used to map regional biological heterogeneity of neuronal and immune signatures in HGG, and continues to offer a promising strategy for non-invasive surveillance of tumor progression. Here, we present a graph network analysis relating mpMRI to biological signatures in the NE tumor rim and show the utility of communities of imaging features for discriminating tumor phenotypes. METHODS 53 mpMRI features across 101 spatially-localized RNA-sequenced biopsy samples from 64 patients with IDH-wt HGG were represented as nodes within a Neo4j graph network. Modularity optimization was used to identify communities of densely connected nodes. Communities of nodes representing high proportions of NE samples (one-tailed test) were kept for further analysis. Gene set enrichment analysis (GSEA) on samples within selected imaging communities was analyzed for biological process enrichment using clusterProfiler. Single sample GSEA (ssGSEA) using ssGSEA2 was used to annotate neuronal and immune cell states across samples within selected communities. mpMRI feature selection within each community identified the most relevant imaging features for isolating samples classified into neuronal and immune cell states. RESULTS 11 imaging communities gathered a high proportion of NE samples. Three NE communities showed significant enrichment of neuronal signatures, while one showed significant enrichment of immune pathways. One neuronal-predominant community showed unique upregulation of 14 pathways related to neuronal migration/astrocyte projection. CONCLUSION Network imaging communities can identify clusters of NE glioma samples showing contrasting neuronal and immune signatures, including distinct upregulation of migration and projection pathways. Mapping these imaging feature communities to specific biological processes implicated in glioma progression may aid treatment planning and prognostication for glioma patients.