Histone acetylation gene‐based biomarkers as novel markers of the immune microenvironment in glioblastoma
Authors: Yan Xiuyou , Feng Lu , Xu Zhengbao , Chen Wei , Yan Hao , Wu Panxing , Ding Chao , Zhu Xuanhao , Lu Yansong 
Published: April 2023

Background

Glioblastoma (GBM) is a primary malignant tumour with high intracranial morbidity, high malignancy and poor prognosis. Abnormal changes in histone acetylation are closely related to the occurrence and development of cancer. However, there is still a lack of systematic research on histone acetylation in GBM.

Methods

Whole‐transcriptome sequencing data and clinical data of GBM patients were obtained through the TCGA database. Single‐cell RNA‐sequencing (scRNA‐seq) data from GBM patients were obtained from GSE146711 in the Gene Expression Omnibus database. Cell descending fractionation was first performed for scRNA‐seq on GBM. The CellChat and PROGENy scores explore the impact of the histone acetylation pathway in GBM on intercellular chat and tumour pathways. The AddModuleScore function evaluates the enrichment score of histone acetylation in cells and divides them into high‐histone acetylation and low‐histone acetylation groups. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses were performed on the differential genes between different histone acetylation states, and the biological processes and pathways that may be affected by histone acetylation were evaluated. Based on this, a prognostic model was constructed using least absolute shrinkage and selection operator (LASSO) analysis, and survival analysis was performed to evaluate its prognostic performance. Finally, we also analysed the main effects of the constructed histone acetylation‐related model on GBM immune infiltration by multiple methods, and analysed the main mutation data of its different subgroups.

Results

GBM samples mainly include seven large cell populations: oligodendrocyte precursor cells (OPCs), myeloid, neoplastic, oligodendrocytes, astrocytes, vascular and neurons. Cellchat and ProgenY scores revealed that in GBM tumours, histone acetylation interacts closely with multiple immune cells and tumour pathways. GO and KEGG analyses revealed the main impact proteins and pathway correlates of histone acetylation. Five histone acetylation genes were screened using LASSO analysis and a prognostic model was constructed. The results revealed that prognostic models were significant in the prognostic stratification of patients in both the training and validation groups of GBM patients. Immune infiltration analysis revealed that the mechanism of histone acetylation in GBM may be related to the immune infiltration of multiple effector immune cells.

Conclusions

Our histone acetylation‐based biomarkers are closely associated with immune microenvironmental infiltration and functional mutations in multiple tumour pathways in GBM. This suggests that histone acetylation may reveal microscopic alterations in the tumour microenvironment, and may provide potential evidence and a research basis for the development of novel therapeutic targets for GBM. On this basis, a novel perspective on the spatial biology and immunological understanding of GBM is provided.