Charting the single-cell and spatial landscape of IDH-wild-type glioblastoma with GBmap

Background. Glioblastoma (GB), particularly IDH-wild type, is the most aggressive brain malignancy with a dismal prognosis. Despite advances in molecular profiling, the complexity of its tumor microenvironment and spatial organization remains poorly understood. This study aimed to create a comprehensive single-cell and spatial atlas of GB to unravel its cellular heterogeneity, spatial architecture, and clinical relevance. Methods. We integrated single-cell RNA sequencing data from 26 datasets, encompassing over 1.1 million cells from 240 patients, to construct GBmap, a harmonized single-cell atlas. High-resolution spatial transcriptomics was employed to map the spatial organization of GB tissues. We developed the Tumor Structure Score (TSS) to quantify tumor organization and correlated it with patient outcomes. Results. We showcase the applications of GBmap for reference mapping, transfer learning, and biological discoveries. GBmap revealed extensive cellular heterogeneity, identifying rare populations such as tumor-associated neutrophils and homeostatic microglia. Spatial analysis uncovered 7 distinct tumor niches, with hypoxia-dependent niches strongly associated with poor prognosis. The TSS demonstrated that highly organized tumors, characterized by well-defined vasculature and hypoxic niches, correlated with worse survival outcomes. Conclusions. This study provides a comprehensive resource for understanding glioblastoma heterogeneity and spatial organization. GBmap and the TSS provide an integrative view of tumor architecture in GB, highlighting hypoxia-driven niches that may represent avenues for further investigation. Our resource can facilitate exploratory analyses and hypothesis generation to better understand disease progression.