A glycosylation-related gene signature predicts prognosis, immune microenvironment infiltration, and drug sensitivity in glioma

Abstract

Glioma represents the most common primary cancer of the central nervous system in adults. Glycosylation is a prevalent post-translational modification that occurs in eukaryotic cells, leading to a wide array of modifications on proteins. We obtained the clinical information, bulk RNA-seq data, and single-cell RNA sequencing (scRNA-seq) from The Cancer Genome Atlas (TCGA), Chinese Glioma Genome Atlas (CGGA), Gene Expression Omnibus (GEO), and Repository of Molecular Brain Neoplasia Data (Rembrandt) databases. RNA sequencing data for normal brain tissues were accessed from the Genotype-Tissue Expression (GTEx) database. Then, the glycosylation genes that were differentially expressed were identified and further subjected to variable selection using a least absolute shrinkage and selection operator (LASSO)-regularized Cox model. We further conducted enrichment analysis, qPCR, nomogram, and single-cell transcriptome to detect the glycosylation signature. Drug sensitivity analysis was also conducted. A five-gene glycosylation signature (CHPF2, PYGL, GALNT13, EXT2, and COLGALT2) classified patients into low- or high-risk groups. Survival analysis, qPCR, ROC curves, and stratified analysis revealed worse outcomes in the high-risk group. Furthermore, GSEA and immune infiltration analysis indicated that the glycosylation signature has the potential to predict the immune response in glioma. In addition, four drugs (crizotinib, lapatinib, nilotinib, and topotecan) showed different responses between the two risk groups. Glioma cells had been classified into seven lines based on single-cell expression profiles. The five-gene glycosylation signature can accurately predict the prognosis of glioma and may offer additional guidance for immunotherapy.