Glioblastomas are malignant tumors of the central nervous system hallmarked by subclonal diversity and dynamic adaptation amid developmental hierarchies. The source of dynamic reorganization within the spatial context of these tumors remains elusive. Here, we characterized glioblastomas by spatially resolved transcriptomics, metabolomics, and proteomics. By deciphering regionally shared transcriptional programs across patients, we infer that glioblastoma is organized by spatial segregation of lineage states and adapts to inflammatory and/or metabolic stimuli, reminiscent of the reactive transformation in mature astrocytes. Integration of metabolic imaging and imaging mass cytometry uncovered locoregional tumor-host interdependence, resulting in spatially exclusive adaptive transcriptional programs. Inferring copy-number alterations emphasizes a spatially cohesive organization of subclones associated with reactive transcriptional programs, confirming that environmental stress gives rise to selection pressure. A model of glioblastoma stem cells implanted into human and rodent neocortical tissue mimicking various environments confirmed that transcriptional states originate from dynamic adaptation to various environments.

Microglia are tissue-resident macrophages of the CNS that orchestrate local immune responses and contribute to several neurological and psychiatric diseases. Little is known about human microglia and how they orchestrate their highly plastic, context-specific adaptive responses during pathology. Here we combined two high-dimensional technologies, single-cell RNA-sequencing and time-of-flight mass cytometry, to identify microglia states in the human brain during homeostasis and disease. This approach enabled us to identify and characterize a previously unappreciated spectrum of transcriptional states in human microglia. These transcriptional states are determined by their spatial distribution, and they further change with aging and brain tumor pathology. This description of multiple microglia phenotypes in the human CNS may open promising new avenues for subset-specific therapeutic interventions.

Reactive astrocytes evolve after brain injury, inflammatory and degenerative diseases, whereby they undergo transcriptomic re-programming. In malignant brain tumors, their function and crosstalk to other components of the environment is poorly understood. Here we report a distinct transcriptional phenotype of reactive astrocytes from glioblastoma linked to JAK/STAT pathway activation. Subsequently, we investigate the origin of astrocytic transformation by a microglia loss-of-function model in a human organotypic slice model with injected tumor cells. RNA-seq based gene expression analysis of astrocytes reveals a distinct astrocytic phenotype caused by the coexistence of microglia and astrocytes in the tumor environment, which leads to a large release of anti-inflammatory cytokines such as TGFβ, IL10 and G-CSF. Inhibition of the JAK/STAT pathway shifts the balance of pro- and anti-inflammatory cytokines towards a pro-inflammatory environment. The complex interaction of astrocytes and microglia cells promotes an immunosuppressive environment, suggesting that tumor-associated astrocytes contribute to anti-inflammatory responses.

Despite recent advances in cancer immunotherapy, certain tumor types, such as Glioblastomas, are highly resistant due to their tumor microenvironment disabling the anti-tumor immune response. Here we show, by applying an in-silico multidimensional model integrating spatially resolved and single-cell gene expression data of 45,615 immune cells from 12 tumor samples, that a subset of Interleukin-10-releasing HMOX1⁺ myeloid cells, spatially localizing to mesenchymal-like tumor regions, drive T-cell exhaustion and thus contribute to the immunosuppressive tumor microenvironment. These findings are validated using a human ex-vivo neocortical glioblastoma model inoculated with patient derived peripheral T-cells to simulate the immune compartment. This model recapitulates the dysfunctional transformation of tumor infiltrating T-cells. Inhibition of the JAK/STAT pathway rescues T-cell functionality both in our model and in-vivo, providing further evidence of IL-10 release being an important driving force of tumor immune escape. Our results thus show that integrative modelling of single cell and spatial transcriptomics data is a valuable tool to interrogate the tumor immune microenvironment and might contribute to the development of successful immunotherapies.

Cancer-associated fibroblasts (CAFs) were presumed absent in glioblastoma given the lack of brain fibroblasts. Serial trypsinization of glioblastoma specimens yielded cells with CAF morphology and single-cell transcriptomic profiles based on their lack of copy number variations (CNVs) and elevated individual cell CAF probability scores derived from the expression of 9 CAF markers and absence of 5 markers from non-CAF stromal cells sharing features with CAFs. Cells without CNVs and with high CAF probability scores were identified in single-cell RNA-Seq of 12 patient glioblastomas. Pseudotime reconstruction revealed that immature CAFs evolved into subtypes, with mature CAFs expressing actin alpha 2, smooth muscle (ACTA2). Spatial transcriptomics from 16 patient glioblastomas confirmed CAF proximity to mesenchymal glioblastoma stem cells (GSCs), endothelial cells, and M2 macrophages. CAFs were chemotactically attracted to GSCs, and CAFs enriched GSCs. We created a resource of inferred crosstalk by mapping expression of receptors to their cognate ligands, identifying PDGF and TGF-β as mediators of GSC effects on CAFs and osteopontin and HGF as mediators of CAF-induced GSC enrichment. CAFs induced M2 macrophage polarization by producing the extra domain A (EDA) fibronectin variant that binds macrophage TLR4. Supplementing GSC-derived xenografts with CAFs enhanced in vivo tumor growth. These findings are among the first to identify glioblastoma CAFs and their GSC interactions, making them an intriguing target.