Langworthy, Jazlynn

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 It has been widely acknowledged that limited access to the human brain, along with the importance of its structures, presents significant challenges in obtaining brain tumor tissue. This poses difficulties for clinical decision-making, as healthcare professionals often have to rely solely on imaging to address complex problems. There is currently no definitive method to make these decisions other than relying on informed yet subjective human judgment. That said, imaging abnormalities are a result of a presence of particular tumor-related biology. Because of this, and the heterogeneity frequently observed in brain tumors, we have developed an image-localized biopsy protocol that has been employed for research purposes at Mayo Clinic Arizona for over five years. This multidisciplinary approach involves the collaborative efforts from neurosurgeons, neuro-oncologists, scrub-techs, nurses, radiologists, as well as multiple research labs, all working within a well-established pipeline. The implementation of this protocol has allowed us to obtain image-localized biopsies during surgical procedures, enabling the correlation of local imaging findings with corresponding tissue features through the utilization of machine learning models. The primary objective of this endeavor is to enhance our understanding of the underlying biological mechanisms of brain tumors based on patients' imaging data. This protocol has been successful for our research purposes, and with minor adjustments, it can be readily adopted by any hospital in order to collect more data about these incredibly heterogeneous diseases. While the use of neuronavigation during surgery has become commonplace in many medical institutions, localization of tissue samples is not. Utilizing pre-existing infrastructure in the operating room is a small burden with big potential positive outcomes. Sample locations can be useful to clinicians during tumor board discussions, to better assess disease heterogeneity. These minor refinements in the surgical workflow have the capacity to provide clinicians with additional clinically-relevant data without disrupting patient care.

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 AIMS High-grade glioma continues to have dismal survival with current standard-of-care treatment, owing in part to its intra- and inter-patient heterogeneity. Typical diagnostic biopsies are taken from the dense tumor core to determine the presence of abnormal cells and the status of a few key genes (e.g. IDH1, MGMT). However, the tumor core is typically resected, leaving behind possibly genetically, transcriptomically and/or phenotypically distinct invasive margins that repopulate the disease. As these remaining populations are the ones ultimately being treated, it is important to know their compositional differences from the tumor core. We aim to identify the phenotypic niches defined by the relative composition of key cellular populations and understand their variation amongst patients. METHOD We have established an image-localized research biopsy study, that samples from both the invasive margin and tumor core. From this protocol, we currently have 202 samples from 58 patients with available bulk RNA-Seq, collected between Mayo Clinic and Barrow Neurological Institute. Using a single-cell reference dataset from our collaborators at Columbia University, we used CIBERSORTx, a deconvolution method, to predict relative abundances of 7 normal, 6 glioma, and 5 immune cell states for each sample. RESULTS We find that these cell state abundances connect to patient survival and show regional differences. For example, proneural glioma states were higher in invasive regions, whereas proliferative and mesenchymal states were higher in the tumor core. CONCLUSIONS Our analysis demonstrates a need to characterize the residual tissue following glioma resection to better understand the recurrent disease.

Abstract Characterizing intra- and inter-tumoral heterogeneity of glioblastoma has historically relied on discrete classifications of malignant cell populations leaving immune and other cell populations, known to exist admixed with the malignant tumor cells, relatively neglected. Manifold learning algorithms can manage deconvolving multiple cell populations and are often used to track cell state transitions in single cell transcriptomics. We applied a manifold learning approach to TCGA microarray data (Nf525) and bulk transcriptomics of 134 image localized biopsies across 30 patients with primary and 9 with recurrent glioblastoma to further elucidate how to organize biopsies across a spectrum of possible tissue states. The algorithm revealed a low-dimensional manifold graph for which each biopsy lives across 3 polarizing tissue states - one that is associated with diffusely invaded brain, one that is enriched in mesenchymal genes, and one that is enriched in classical proliferative tumor signatures. We deconvolved the bulk transcriptomics of the image-localized biopsies to reveal the relative abundance of 18 malignant, immune, and other cell subpopulations in each biopsy. Overlaying the cellular decomposition onto the manifold graph visualizing the tissue state distributions revealed that transitions between states correlate with changes in cellular cohabitation composition. The tumor cellular cohabitation ecologies have the lowest diversity, as inferred by ecological measures such as Shannon entropy and evenness, at the distal poles of the graph when compared to the transitional arms. Further, we found that the relationship between imaging appearance of contrast enhancement on T1-weighted MRI and the biopsy cellular composition varies with sex and primary vs recurrent biopsy status. The limited spectrum of possible tissue states revealed by the manifold learning is suggestive of a limited continuum along which tumor and non-tumoral cell populations can cohabitate. Such a limited low-dimensional biological space may constrain the dynamics of tumor biology in a predictable manner.

Abstract Glioblastoma (GBM) is the most aggressive primary brain tumor with a median survival of 14 months. GBMs are challenging to treat due to their heterogeneous nature. It has also been seen that these tumors have sex differences in their cellular subtypes as well as imaging. Radiomics has the potential to provide a non-invasive, spatial understanding of genetic and epigenetic diversity in these complex tumors and to aid in treatment planning. We have an ongoing study to obtain image-localized biopsies from GBM patients, allowing us to complete radiomic analysis and make connections between immunohistochemistry (IHC) and magnetic resonance imaging (MRI) features. We sought to determine if the patterns on imaging were correlated with underlying tumor biology. We focused on immunohistochemistry (IHC) markers of key features of tumor biology including SOX2 for stem-like tumor cells, CD68 for immune response and Ki67 for proliferation kinetics. Our study included 38 patients with a total of 99 biopsies (bxs): 27 males with 77 bxs and 11 females with 22 bxs. Biopsies were sectioned and stained for the SOX2, CD68, and KI67 markers. We computed 18 first-order radiomic features at each biopsy location for patients’ multimodal MRIs: T1W, T1Gd, T2W, FLAIR, apparent diffusion coefficient, diffusion weighted imaging (DWI) and susceptibility weighted imaging. We then performed correlation analysis between each radiomic feature and marker abundance for each IHC stain. Overall, we found sex-distinct patterns connecting imaging with these IHC markers. For example, amongst female patients, DWI held more prominent correlations with SOX2 than in males. Whereas there were more correlations between CD68 IHC abundance and T1Gd imaging features in males compared to females. Taken together, the overall patterns connecting locoregional imaging features to these IHC markers showed sex-distinct patterns suggesting the potential for sex to be an important biological variable when interpreting the biology underlying imaging changes.

Abstract Background Magnetic Resonance Imaging (MRI) is the mainstay for neurosurgical oncology but not for informing us about glioma biology. An obstacle to developing MR-based glioma biomarkers is the absence of rigorous correlation between MRI features and glioma biology as assessed in multi-regional biopsies, within and across patients. Methods We directly addressed this obstacle by collating a unique cohort of 202 MRI-localized biopsies from 58 patients. We define a low-dimensional transcriptional pseudotime continuum along which heterogeneous high-grade glioma (HGG) samples organize both within and across patients. Results We observe three polarized transcriptional tissue states: infiltrated brain, immune/inflammatory, and proliferative associated with patterns of cohabitation of cellular subpopulations. The states and deconvolved populations show correlation with enhancement status on T1Gd MRI. Moreover, discrete MRI habitats, regions sharing common imaging features, defined as combinations of high or low signal intensity across multiparametric MRI revealed 14 MRI habitats. We order the MRI habitats according to the average pseudotime on the transcriptional continuum. We find that MRI habitats with low pseudotime (associated with early tumor development and diffusely invaded brain tissue) localized at the periphery of the tumor whilst high pseudotime either proliferative or immune/inflammatory states were towards the core of the lesion. We find that composition of MRI habitats are impacted by MGMT status and patient sex. Conclusion This suggests that ongoing aggregation of MRI-localized biopsies may augment our projection of biology onto MRI habitats to support tracking of the evolution of cellular ecologies within and across each patient’s tumor over time. Key Points Trajectory inference of HGG samples reveals a continuum connecting three polarized tissue states reflecting transitions in cellular population ecologies and gene set enrichment Combining multiparametric MRI features into regional MRI habitats reveal clinically-distinct patterns of cell co-habitation and tissue state MGMT methylation status shifts the cellular patterns of cohabitation within MRI habitats Importance of the Study Regional heterogeneity is a defining characteristic of high-grade glioma (HGG). Mapping the transcriptomic biology of multiregionally sampled MRI-localized biopsies from high-grade gliomas onto a low-dimensional space reveals a continuum between three primary tissue states: infiltrated brain, immune/inflammatory, and proliferative. Each state is associated with distinct cell cohabitation ecologies. Transitions between the tissue states coordinate with changes in tumor composition, revealing potential direct connections between tumor biology shifts and their manifestation on clinical imaging.