Abstract
Glioblastoma is a heterogeneous tumor for which effective treatment options are limited and often insufficient. Few studies have examined the intratumoral transcriptional and proteomic heterogeneity of the glioblastoma microenvironment to characterize the spatial distribution of potential molecular and cellular therapeutic immunooncology targets. We applied an integrated multimodal approach comprised of NanoString GeoMx Digital Spatial Profiling, single-cell RNA-seq (scRNA-seq), and expert neuropathologic assessment to characterize archival formalin-fixed paraffin-embedded glioblastoma specimens. Clustering analysis and spatial cluster maps highlighted the intratumoral heterogeneity of each specimen. Mixed cell deconvolution analysis revealed that neoplastic and vascular cells were the prominent cell types throughout each specimen, with macrophages, oligodendrocyte precursors, neurons, astrocytes, and oligodendrocytes present in lower abundance and illustrated the regional distribution of the respective cellular enrichment scores. The spatial resolution of the actionable immunotherapeutic landscape showed that robust B7H3 gene and protein expression was broadly distributed throughout each specimen and identified STING and VISTA as potential targets. Lastly, we uncovered remarkable variability in VEGFA expression and discovered unanticipated associations between VEGFA, endothelial cell markers, hypoxia, and the expression of immunoregulatory genes, indicative of regionally distinct immunosuppressive microdomains. This work provides an early demonstration of the ability of an integrated panel-based spatial biology approach to characterize and quantify the intrinsic molecular heterogeneity of the glioblastoma microenvironment.