Identification and Validation of Prognostic Genes Related to Histone Lactylation Modification in Glioblastoma: An Integrated Analysis of Transcriptome and Single-cell RNA Sequencing.

Abstract

Background: The impact of histone lactylation modification (HLM) on glioblastoma (GBM) progression is not well understood. This study aimed to identify HLM-associated prognostic genes in GBM and explore their mechanisms of action. Methods: The presence and role of lactylation in glioma clinical tissue samples and its correlation with GBM progression were analysed through immunohistochemical staining and Western blotting. Sequencing data for GBM were obtained from publicly available databases. An initial correlation analysis was performed between global HLM levels and GBM. Differentially expressed HLM-related genes (HLMRGs) in GBM were identified by intersecting differentially expressed genes (DEGs) from the TCGA-GBM dataset, key module genes derived from weighted gene coexpression network analysis (WGCNA), and previously reported HLMRGs. Prognostic genes were subsequently identified through univariate Cox regression and least absolute shrinkage and selection operator (LASSO) regression analyses, which formed the basis for constructing a risk prediction model. Associations between HLMRGs and GBM were further evaluated via single-cell RNA sequencing (scRNA-seq) datasets. Complementary analyses, including functional enrichment, immune infiltration, somatic mutation, and nomogram-based survival prediction, were conducted alongside in vitro experiments. Additionally, drug sensitivity and Chinese medicine prediction analyses were performed to identify potential therapeutic agents for GBM. Results: We detected a significant increase in global lactylation levels in GBM, which correlated with patient survival. We identified 227 differentially expressed HLMRGs from the intersection of 3,343 differentially expressed genes and 948 key module genes, indicating strong prognostic potential. Notably, genes such as SNCAIP, TMEM100, NLRP11, HOXC11, and HOXD10 were highly expressed in GBM. Functional analysis suggested that HLMRGs are involved primarily in pathways related to cytokine‒cytokine receptor interactions, cell cycle regulation, and cellular interactions, including microglial differentiation states. Further connections were established between HLMRGs and infiltrating immune cells, particularly type 1 T helper (Th1) cells, as well as mutations in genes such as PTEN. The potential therapeutic agents identified included ATRA and Can Sha. Conclusion: The HLM-related gene risk prediction model shows substantial promise for improving patient management in GBM, providing crucial insights for clinical prognostic evaluations and immunotherapeutic approaches in GBM.