Single‐cell RNA sequencing reveals tumor heterogeneity, microenvironment, and drug‐resistance mechanisms of recurrent glioblastoma

Abstract Glioblastomas are highly heterogeneous brain tumors. Despite the availability of standard treatment for glioblastoma multiforme (GBM), i.e., Stupp protocol, which involves surgical resection followed by radiotherapy and chemotherapy, glioblastoma remains refractory to treatment and recurrence is inevitable. Moreover, the biology of recurrent glioblastoma remains unclear. Increasing evidence has shown that intratumoral heterogeneity and the tumor microenvironment contribute to therapeutic resistance. However, the interaction between intracellular heterogeneity and drug resistance in recurrent GBMs remains controversial. The aim of this study was to map the transcriptome landscape of cancer cells and the tumor heterogeneity and tumor microenvironment in recurrent and drug‐resistant GBMs at a single‐cell resolution and further explore the mechanism of drug resistance of GBMs. We analyzed six tumor tissue samples from three patients with primary GBM and three patients with recurrent GBM in which recurrence and drug resistance developed after treatment with the standard Stupp protocol using single‐cell RNA sequencing. Using unbiased clustering, nine major cell clusters were identified. Upregulation of the expression of stemness‐related and cell‐cycle‐related genes was observed in recurrent GBM cells. Compared with the initial GBM tissues, recurrent GBM tissues showed a decreased proportion of microglia, consistent with previous reports. Finally, vascular endothelial growth factor A expression and the blood–brain barrier permeability were high, and the O 6 ‐methylguanine DNA methyltransferase‐related signaling pathway was activated in recurrent GBM. Our results delineate the single‐cell map of recurrent glioblastoma, tumor heterogeneity, tumor microenvironment, and drug‐resistance mechanisms, providing new insights into treatment strategies for recurrent glioblastomas.