Hu, Leland

Background/Objectives: White matter (WM) tract detection is critical in the presurgical planning of tumor resection. However, standard-of-care imaging techniques including T1-weighted, T2-weighted, and Diffusion Tensor Imaging (DTI) often fail to identify WM tracts within edematous regions. In T1/T2-weighted imaging, edema increases extracellular water and reduces tissue contrast, and in diffusion-weighted imaging, edema elevates isotropic diffusion, reducing sensitivity to anisotropic diffusion along WM tracts. Advanced biophysical diffusion modeling techniques such as Neurite Orientation Dispersion and Density Imaging (NODDI) and the Standard Model (SM) address this limitation by compartmentalizing the diffusion signal into free-water, intra-neurite, and extra-neurite contributions. Here, we test if biophysical multi-compartment models can robustly identify WM tracts and recover tractography streamlines within edematous regions. Methods: In this study, we use multi-shell diffusion-weighted MRI data obtained from patients with meningiomas—a pathology allowing for isolation of the effects of edema without the confounding effects of tumor cell invasion. We compared FA from standard and free-water-corrected DTI, the orientation dispersion index (ODI) from NODDI, and P2 (a scalar descriptor of fiber orientation coherence) from the SM fODF in edematous and unaffected contralateral WM regions. As a proof of concept, we visually evaluated the tractography performance across models. Results: Our results show that (1 − ODI) and P2 values in edema remained close to within-subject contralateral measurements, contrasting with substantial reductions in FA and FW-FA. (1 − ODI) showed a small but statistically significant increase in edema (~8%, p = 0.02), while P2 was unchanged. Conclusions: These results highlight the potential of biophysical diffusion models for preoperative mapping in edema.

Abstract M.D.-Ph.D. programs in the United States have traditionally followed a “2- n -2” curricular model, with the graduate phase occurring between the pre-clerkship and clerkship portions of medical training. While well established, this format can limit trainee autonomy in shaping their physician-scientist training trajectories. In response, some programs have introduced a “3- n -1” model, allowing students to complete clerkships before beginning Ph.D. training. Our institution implemented multiple flexible curricula in 2017. Understanding why trainees choose one pathway is important as programs consider implementing more adaptable training structures. To investigate these factors, we surveyed M.D.-Ph.D. students at our institution, which offers multiple flexible curricular alternatives, about contributions to their curricular decisions. Responses supported that trainees weigh considerations across medical education, scientific development, and integrated physician-scientist training domains. Although the traditional pathway was a popular option, most students pursued one of the flexible pathways. Our findings suggest that introducing flexibility in the timing of undergraduate medical education and graduate training can support diverse educational, logistical, and personal priorities while maintaining the rigor of physician-scientist training. Offering multiple pathways empowers trainees to design trajectories that best fit their needs. Continued longitudinal studies are needed to assess the long-term impacts of flexible curricula on physician-scientist career outcomes.

Lay Summary Accurately differentiating actual tumor progression from treatment-related changes, such as pseudoprogression or radiation-induced injury, remains one of the most critical diagnostic challenges in neuro-oncology for both primary and secondary brain tumors. Conventional MRI often lacks specificity because it primarily reflects blood-brain barrier disruption rather than viable tumor tissue. Consequently, advanced techniques like perfusion-weighted MRI and amino acid PET have gained significant attention. While perfusion-weighted MRI provides essential insights into blood flow and is widely available, accumulating evidence suggests that amino acid PET can be more reliable in identifying treatment-related changes. Since both modalities offer complementary information, there is an ongoing debate regarding their relative roles and the optimal way to integrate them into clinical practice. Combining these tools aims to provide more reliable diagnoses, ultimately preventing unnecessary interventions and ensuring that true tumor growth is detected as early as possible.

Abstract Purpose: Despite advances in molecular diagnostics, glioblastoma (GBM) remains incurable, with recurrence driven by residual invasive tumor cells that escape surgical resection and therapy. Current genomic profiling to guide therapy for recurrent disease typically relies on a single-site clinical biopsy (SSCB) obtained from the contrasting-enhancing (CE) tumor core - the region targeted for maximal resection and most visible on Magnetic Resonance Imaging (MRI). However, this region poorly represents the non-enhancing (NE) invasive tumor margins, which harbor treatment-resistant and genetically distinct cell populations responsible for disease recurrence. GBM exhibits extensive spatial heterogeneity, raising concerns that a single tissue may underrepresent the true actionable landscape of the tumor. Experimental Procedures: We analyzed 302 biopsies from 66 GBM patients who underwent multiregional whole exome and RNA sequencing across CE and NE tumor regions. Of these, 17 patients also had standard-of-care SSCB molecular profiling utilizing clinically available genomic sequencing panels outside of the multiregional protocol. Variants that have been targeted in ongoing or successful human clinical trials were annotated as actionable. We compared actionability profiles between SSCB and multiregional sampling, to determine the gain of actionable variant detection across the two sampling strategies. We performed an in-silico biopsy simulation to model the coverage of these alterations across the tumor landscape with increasing sampling frequency in the multiregional protocol, versus the coverage from clinical sequencing panels on SSCBs. Results: Clinical sampling identified a median of one actionable mutation per patient (range 0-3). In contrast, multiregional profiling revealed additional potentially actionable variants in over 70% of patients, particularly within NE infiltrative regions that are not captured by conventional resection. Notably, nine potentially actionable copy-number variants, including PDGFRA amplification, were consistently absent from SSCB results and underrepresented in most clinical panels. In silico biopsy simulations demonstrated that actionable variant detection continues to rise even after sampling five additional regions, reflecting spatial segregation of therapeutic targets. Conclusions: These findings demonstrate that a single sample from the CE tumor core fails to capture the biologically diverse invasive tumor populations responsible for GBM recurrence. Multiregional sequencing uncovers clinically meaningful alterations that would otherwise remain undetected with conventional sampling, highlighting the risk of incomplete molecular characterization in current clinical biopsy workflows. Citation Format: Enmar Harfoush, Taylor Weiskittel, Mylan Blomquist, Dustin Grief, Aadel A. Chaudhuri, Leland S. Hu, Nhan L. Tran, Shannon Fortin-Ensign. Beyond a single core biopsy: Spatial genomic profiling improves therapeutic target detection in glioblastoma [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 706.

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 Glioblastoma (GBM) invasion is shaped by dynamic clonal evolution within the spatial and cellular architecture of the brain. To investigate this process, we integrated multiparametric MRI, whole-exome sequencing (WES), bulk RNA-seq and single-nucleus RNA-seq (snRNA-seq) from 69 IDH-wildtype GBM patients, using spatially resolved biopsies from contrast-enhancing (CE) core tumor and non-enhancing (NE) invasive regions. Phylogenetic analysis revealed distinct clonal architectures. CE-enriched clones carried truncal mutations in conanical drivers of cell cycle and RTK signaling (TP53, PTEN, EGFR, PDGFRA/B), and private alterations in DNA repair genes (CDH11, PAX7, RAD51B). In contrast, NE-enriched clones exhibited private mutations in MAPK and Hippo pathways components (STAG2, BRAF, MAP2K3, CTNNA1), suggesting region-specific evolutionary trajectories. Furthermore, snRNAseq cell-cell communication analysis uncovered enriched ligand-receptor signaling in NE regions between malignant cells and brain-resident populations, including PTN—PTPRZ1 (myeloid), NRG-ERBB4 and SLIT-ROBO (neuronal, astrocytic, and oligodendrocytic). Bulk RNAseq deconvolution confirmed increased representation of neurons and glial cells in biopsies harboring NE-specific clones. These findings define a unique population of GBM clones localized at the infiltrative margin, marked by private MAPK/Hippo mutations and enhanced intracellular communication with the surrounding brain microenvironment. The recurrent engagement of neural-glial signaling pathways suggests that these NE-associated clones co-opt normal cellular networks to support invasion, immune evasion, and persistence following surgical resection. This study establishes a mechanistic link between spatial clonal evolution and tumor-host communication at the GBM infiltrative region and provide a framework for targeting cell-cell signaling circuits that sustain invasive tumor ecosystems and escape therapeutic 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.