Single-Cell and Machine Learning-Based Pyroptosis-Related Gene Signature Predicts Prognosis and Immunotherapy Response in Glioblastoma

Liren Fang^1*Desheng Wang¹Fanlei Meng¹Yinzhi Wang¹Lu Feng^2*Hong Li¹

• ¹Department of Neurology, Second Hospital of Tianjin Medical University, Tianjin, China
• ²Taizhou Central Hospital, Taizhou, China

Background Glioblastoma (GBM) is the most aggressive primary malignancy of the central nervous system, characterized by profound heterogeneity and an immunosuppressive microenvironment, leading to dismal prognosis. Pyroptosis, an inflammatory form of programmed cell death, has been increasingly linked to tumor immunity and progression; however, its molecular roles and clinical implications in GBM remain insufficiently understood. Methods We integrated bulk transcriptome profiles from TCGA-GBM, CGGA, and GEO datasets with single-cell RNA sequencing data from GSE141383 and GSE223063. A comprehensive GBM single-cell atlas was constructed using Seurat and Harmony, and malignant epithelial cells were inferred via inferCNV. Pyroptosis activity was quantified by five complementary algorithms, while Monocle2 and Slingshot were employed for pseudotime trajectory reconstruction, and SCENIC was applied for transcription factor network analysis. Candidate prognostic genes identified from malignant epithelial subsets were further used to develop a Pyroptosis-Related Gene Signature (PRGS) through a systematic evaluation of ten machine learning algorithms and their combinations, with subsequent validation across multiple cohorts. Functional enrichment (GSVA, GSEA), tumor microenvironment estimation (ESTIMATE, ssGSEA), drug sensitivity prediction (GDSC2), and in vitro experiments were performed to characterize the biological and therapeutic relevance of PRGS, with MAP1B selected for experimental validation. Results Single-cell analyses revealed heterogeneous pyroptosis activity across GBM cell populations. Distinct ligand–receptor communications were observed between high-and low-pyroptosis groups, among which the SPP1-centered signaling axis showed pronounced remodeling, suggesting a pivotal role in tumor–immune crosstalk. Pseudotime and regulatory network analyses of malignant epithelial cells further delineated differentiation trajectories and transcriptional regulators. The PRGS, established by StepCox[both]+Ridge modeling, demonstrated robust prognostic stratification and predictive power across independent datasets. High PRGS scores were consistently associated with poorer survival outcomes, higher TIDE scores, and reduced IPS values, indicating enhanced immune evasion and attenuated immunotherapy benefit. Enrichment analyses highlighted that high PRGS tumors were linked to metabolic reprogramming and DNA repair pathways, whereas low PRGS tumors exhibited signatures of immune activation. Drug sensitivity analyses revealed distinct therapeutic vulnerabilities between subgroups. Functional assays confirmed that MAP1B promotes proliferation, migration, and invasion in GBM cells, reinforcing its oncogenic role.