Abstract:Background: Glioblastoma (GBM) is adults' most common and fatally malignant brain tumor. The heterogeneity is the leading cause of treatment failure. However, the relationship between cellular heterogeneity, tumor microenvironment, and GBM progression is still elusive. Methods: Integrated analysis of single-cell RNA sequencing (scRNA-seq) and spatial transcriptome sequencing (stRNA-seq) of GBM were conducted to analyze the spatial tumor microenvironment. We investigated the subpopulation heterogeneity of malignant cells through gene set enrichment analyses, cell communications analyses, and pseudotime analyses. Significantly changed genes of the pseudotime analysis were screened to create a tumor progress-related gene risk score (TPRGRS) using Cox regression algorithms in the bulkRNA-sequencing(bulkRNA-seq) dataset. We combined the TPRGRS and clinical characteristics to predict the prognosis of patients with GBM. Furthermore, functional analysis was applied to uncover the underlying mechanisms of the TPRGRS. Results: GBM cells were accurately charted to their spatial locations and uncovered their spatial colocalization. The malignant cells were divided into five clusters with transcriptional and functional heterogeneity, including unclassified malignant cells and astrocyte-like, mesenchymal-like, oligodendrocytes-progenitor-like, and neural-progenitor-like malignant cells. Cell-cell communications analysis in scRNA-seq and stRNA-seq identified ligand-receptor pairs of the CXCL, EGF, FGF, and MIF signaling pathways as bridges implying that tumor microenvironment may cause malignant cells' transcriptomic adaptability and disease progression. Pseudotime analysis showed the differentiation trajectory of GBM cells from proneural to mesenchymal transition and identified genes or pathways that affect cell differentiation. TPRGRS could successfully divide patients with GBM in three datasets into high- and low-risk groups, which was proved to be a prognostic factor independent of routine clinicopathological characteristics. Functional analysis revealed the TPRGRS associated with growth factor binding, cytokine activity, signaling receptor activator activity functions, and oncogenic pathways. Further analysis revealed the association of the TPRGRS with gene mutations and immunity in GBM. Finally, the external datasets and qRT-PCR verified high expressions of the TPRGRS mRNAs in GBM cells. Conclusion: Our study provides novel insights into heterogeneity in GBM based on scRNA-seq and stRNA-seq data. Moreover, our study proposed a malignant cell transition-based TPRGRS through integrated analysis of bulkRNA-seq and scRNA-seq data, combined with the routine clinicopathological evaluation of tumors, which may provide more personalized drug regimens for GBM patients.

Single‐Cell RNA Sequencing Reveals Tumor Heterogeneity, Microenvironment, and Drug‐resistance Mechanisms of Recurrent Glioblastoma

Comprehensive Analysis of Multi-Omics Data of Recurrent Gliomas Identifies a Recurrence-Related Signature As a Novel Prognostic Marker.

Single-cell atlas reveals the immunosuppressive microenvironment and Treg cells landscapes in recurrent Glioblastoma

Relapse pathway of glioblastoma revealed by single-cell molecular analysis

The heterogeneous cellular states of glioblastoma stem cells revealed by single cell analysis

Integration analysis of single-cell and spatial transcriptomics reveal the cellular heterogeneity landscape in glioblastoma and establish a polygenic risk model

Tumor Heterogeneity and Molecular Characteristics of Glioblastoma Revealed by Single-Cell RNA-Seq Data Analysis

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

Single-cell RNA-seq highlights intratumoral heterogeneity in primary glioblastoma

Comprehensive portrait of recurrent glioblastoma multiforme in molecular and clinical characteristics

Single-cell Sequencing Data Analysis Reveals Heterogeneity of Glioblastoma Multiforme

Single-Cell Atlas Reveals Complexity of the Immunosuppressive Microenvironment of Initial and Recurrent Glioblastoma

Intratumor heterogeneity, microenvironment, and mechanisms of drug resistance in glioma recurrence and evolution

Single‐cell RNA‐seq reveals the invasive trajectory and molecular cascades underlying glioblastoma progression

Abstract 5482: Spatial and multi-omic interrogation of the primary and recurrent glioblastoma microenvironment

Decoding Key Cell Sub-Populations and Molecular Alterations in Glioblastoma at Recurrence by Single-Cell Analysis

Single-Cell Transcriptomics Revealed Subtype-Specific Tumor Immune Microenvironments in Human Glioblastomas

Whole Exome Sequencing Reveals the Genetic Heterogeneity and Evolutionary History of Primary Gliomas and Matched Recurrences

Pervasive structural heterogeneity rewires glioblastoma chromosomes to sustain patient-specific transcriptional programs

Intratumoral Cell Heterogeneity in Patient-Derived Glioblastoma Cell Lines Revealed by Single-Cell RNA-Sequencing

Molecular characteristics of single patient-derived glioma stem-like cells from primary and recurrent glioblastoma