Despite aggressive multimodal treatment, glioblastoma (GBM), a grade IV primary brain tumor, still portends a poor prognosis with a median overall survival of 12–16 months. The complexity of GBM treatment mainly lies in the inter- and intra-tumoral heterogeneity, which largely contributes to the treatment-refractory and recurrent nature of GBM. By paving the road towards the development of personalized medicine for GBM patients, the cancer genome atlas classification scheme of GBM into distinct transcriptional subtypes has been considered an invaluable approach to overcoming this heterogeneity. Among the identified transcriptional subtypes, the mesenchymal subtype has been found associated with more aggressive, invasive, angiogenic, hypoxic, necrotic, inflammatory, and multitherapy-resistant features than other transcriptional subtypes. Accordingly, mesenchymal GBM patients were found to exhibit worse prognosis than other subtypes when patients with high transcriptional heterogeneity were excluded. Furthermore, identification of the master mesenchymal regulators and their downstream signaling pathways has not only increased our understanding of the complex regulatory transcriptional networks of mesenchymal GBM, but also has generated a list of potent inhibitors for clinical trials. Importantly, the mesenchymal transition of GBM has been found to be tightly associated with treatment-induced phenotypic changes in recurrence. Together, these findings indicate that elucidating the governing and plastic transcriptomic natures of mesenchymal GBM is critical in order to develop novel and selective therapeutic strategies that can improve both patient care and clinical outcomes. Thus, the focus of our review will be on the recent advances in the understanding of the transcriptome of mesenchymal GBM and discuss microenvironmental, metabolic, and treatment-related factors as critical components through which the mesenchymal signature may be acquired. We also take into consideration the transcriptomic plasticity of GBM to discuss the future perspectives in employing selective therapeutic strategies against mesenchymal GBM.

In isocitrate dehydrogenase wildtype glioblastoma (GBM), cellular heterogeneity across and within tumors may drive therapeutic resistance. Here we analyzed 121 primary and recurrent GBM samples from 59 patients using single-nucleus RNA sequencing and bulk tumor DNA sequencing to characterize GBM transcriptional heterogeneity. First, GBMs can be classified by their broad cellular composition, encompassing malignant and nonmalignant cell types. Second, in each cell type we describe the diversity of cellular states and their pathway activation, particularly an expanded set of malignant cell states, including glial progenitor cell-like, neuronal-like and cilia-like. Third, the remaining variation between GBMs highlights three baseline gene expression programs. These three layers of heterogeneity are interrelated and partially associated with specific genetic aberrations, thereby defining three stereotypic GBM ecosystems. This work provides an unparalleled view of the multilayered transcriptional architecture of GBM. How this architecture evolves during disease progression is addressed in the companion manuscript by Spitzer et al.

Several immunotherapy clinical trials in recurrent glioblastoma have reported long-term survival benefits in 10–20% of patients. Here we perform genomic analysis of tumor tissue from recurrent WHO grade IV glioblastoma patients acquired prior to immunotherapy intervention. We report that very low tumor mutation burden is associated with longer survival after recombinant polio virotherapy or after immune checkpoint blockade in recurrent glioblastoma patients. A relationship between tumor mutation burden and survival is not observed in cohorts of immunotherapy naïve newly diagnosed or recurrent glioblastoma patients. Transcriptomic analyses reveal an inverse relationship between tumor mutation burden and enrichment of inflammatory gene signatures in cohorts of recurrent, but not newly diagnosed glioblastoma tumors, implying that a relationship between tumor mutation burden and tumor-intrinsic inflammation evolves upon recurrence.

The evolution of isocitrate dehydrogenase (IDH)-wildtype glioblastoma (GBM) after standard-of-care therapy remains poorly understood. Here we analyzed matched primary and recurrent GBMs from 59 patients using single-nucleus RNA sequencing and bulk DNA sequencing, assessing the longitudinal evolution of the GBM ecosystem across layers of cellular and molecular heterogeneity. The most consistent change was a lower malignant cell fraction at recurrence and a reciprocal increase in glial and neuronal cell types in the tumor microenvironment (TME). The predominant malignant cell state differed between most matched pairs, but no states were exclusive or highly enriched in either time point, nor was there a consistent longitudinal trajectory across the cohort. Nevertheless, specific trajectories were enriched in subsets of patients. Changes in malignant state abundances mirrored changes in TME composition and baseline profiles, reflecting the co-evolution of the GBM ecosystem. Our study provides a blueprint of GBM’s diverse longitudinal trajectories and highlights the treatment and TME modifiers that shape them.