Abstract:Objective: IDH-wildtype glioblastoma (GBM) is the most aggressive brain tumor in adults. Notable for substantial inter- and intra-tumoral heterogeneity, GBM is highly resistant to standard-of-care treatment. Emerging spatial omics technologies, such as GeoMx Digital Spatial Profiling (DSP), offer novel approaches to interrogate the molecular features of regionally defined cellular ecosystems in cancer. Here, we leverage DSP to study whole-transcriptome gene expression signatures in a spatially defined manner across five IDH-wildtype GBMs. Methods: Formalin-fixed paraffin-embedded tumors were assessed histologically and three spatially distinct areas from each tumor were used to construct a custom "spatial" tissue microarray. From each set of three areas, we identified nine spatially distinct regions of interest (ROIs) for bulk transcriptomic analysis according to two histopathological features: distribution (cellular tumor core or tumor periphery) and vascular proximity (near or far). We evaluated normalized target expression data statistically using linear mixed modeling (LMM) with Benjamini-Hochberg correction. Results: 44 profiled ROIs and 8531 transcriptomic targets were analyzed following quality control trimming. 3D principal component analysis and unsupervised hierarchical clustering analysis showed clustering according to tumor specimen as well as distribution (core versus periphery), thereby substantiating the presence of inter- and intra-tumoral transcriptional heterogeneity. LMM-based differential gene expression analysis identified distinct molecular signatures when comparing core and periphery. Specifically, we observed elevated log2-transformed expression of migration and proliferation markers including LAMC1 (p=0.00297) and FABP7 (p=0.00434) in core ROIs and neurodevelopmental markers including NTRK2 (p=0.00734) and DAAM2 (p=0.00263) in periphery ROIs. Additionally, reactome analysis demonstrated considerable enrichment for the aberrant programmed cell death pathway (R-HSA-9645723, coverage=71.84%) in core ROIs and the MHC class II antigen presentation pathway (R-HSA-2132295, coverage=69.11%) in periphery ROIs. Interestingly, differential gene expression analysis did not reveal any biologically significant targets when subsetting for vascular proximity within cellular tumor core ROIs. However, reactome analysis of this regional subset demonstrated upregulation of various cellular proliferation pathways (G2/M transition, APC/C activity, DNA synthesis, etc.) within vessel-near ROIs when compared with vessel-far ROIs. Conclusions: We preliminarily characterize the molecular heterogeneity of GBM using a spatially guided approach. Ongoing analyses are focused on dataset validation, mixed cell deconvolution, further regional genotyping, and potential biomarkers. Citation Format: Adam Luo, Kelly Casella, Shivang Sharma, Ezra Baraban, Calixto-Hope Lucas, Eugene Shenderov. Utilizing whole-transcriptome digital spatial profiling for glioblastoma clinical biomarker discovery [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2024; Part 1 (Regular s); 2024 Apr 5-10; San Diego, CA. Philadelphia (PA): AACR; Cancer Res 2024;84(6_Suppl) nr 3766.

Leveraging multi-omics data to infer regulators of mRNA 3' end processing in glioblastoma

Alternative Polyadenylation in Glioblastoma Multiforme and Changes in Predicted RNA Binding Protein Profiles.

The Non-N-6-Methyladenosine Epitranscriptome Patterns and Characteristics of Tumor Microenvironment Infiltration and Mesenchymal Transition in Glioblastoma

RNA Modification Patterns Based on Major RNA Modifications Define Tumor Microenvironment Characteristics in Glioblastoma

The RNA-binding protein SERBP1 functions as a novel oncogenic factor in glioblastoma by bridging cancer metabolism and epigenetic regulation

RNA N6-Methyladenosine Regulator-Mediated Methylation Modifications Pattern and Immune Infiltration Features in Glioblastoma

Identification and classification of glioma subtypes based on RNA-binding proteins

Abstract 5446: The role of snoRNAs in GBM cells, GSCs, neuronal precursor cells (NPCs), astrocytes, and normal brain tissues

Abstract B038: Accurate quantitative profiling of rna modifications and their associations with glioblastoma

Single-cell Multi-Omics Sequencing Uncovers Region-Specific Plasticity of Glioblastoma for Complementary Therapeutic Targeting

Molecular Insights and Prognosis Associated With RBM8A in Glioblastoma

A novel model-based identification and validation of pivotal RNA-binding protein in the evolution of glioma

Identification of Hub Genes and Regulatory Factors of Glioblastoma Multiforme Subgroups by RNA-seq Data Analysis

Transcriptome-wide analysis of circRNA and RBP profiles and their molecular and clinical relevance for GBM

Identification of prognostic biomarkers in glioblastoma using a long non-coding RNA-mediated, competitive endogenous RNA network

Epco-19. Multiomic Profiling Of Epigenetic Heterogeneity In Glioblastoma At Single-Cell Resolution Using Multiomic Single-Cell Cut&Tag

RBBP6 maintains glioblastoma stem cells through CPSF3-dependent alternative polyadenylation

Spatial multi-omics defines a shared infiltrative signature in glioblastoma at the resection margin

Abstract 3766: Utilizing whole-transcriptome digital spatial profiling for glioblastoma clinical biomarker discovery

Exploring the role of m6A methylation regulators in glioblastoma multiforme and their impact on the tumor immune microenvironment

Systematic integration of m6A regulators and autophagy-related genes in combination with long non-coding RNAs predicts survival in glioblastoma multiforme