The present study used the RNA sequencing (RNA-seq) dataset to identify prognostic snoRNAs and construct a prognostic signature of The Cancer Genome Atla (TCGA) lower grade glioma (LGG) cohort, and comprehensive analysis of this signature.
RNA-seq dataset of 488 patients from TCGA LGG cohort were included in this study. Comprehensive analysis including function enrichment, gene set enrichment analysis (GSEA), immune infiltration, cancer immune microenvironment, and connectivity map (CMap) were used to evaluate the snoRNAs prognostic signature.
We identified 21 LGG prognostic snoRNAs and constructed a novel eleven-snoRNA prognostic signature for LGG patients. Survival analysis suggests that this signature is an independent prognostic risk factor for LGG, and the prognosis of LGG patients with a high-risk phenotype is poor (adjusted P = 0.003, adjusted hazard ratio = 2.076, 95% confidence interval = 1.290–3.340). GSEA and functional enrichment analysis suggest that this signature may be involved in the following biological processes and signaling pathways: such as cell cycle, Wnt, mitogen-activated protein kinase, janus kinase/signal transducer and activator of tran-ions, T cell receptor, nuclear factor-kappa B signaling pathway. CMap analysis screened out ten targeted therapy drugs for this signature: 15-delta prostaglandin J2, MG-262, vorinostat, 5155877, puromycin, anisomycin, withaferin A, ciclopirox, chloropyrazine and megestrol. We also found that high- and low-risk score phenotypes of LGG patients have significant differences in immune infiltration and cancer immune microenvironment.
The present study identified a novel eleven-snoRNA prognostic signature of LGG and performed a integrative analysis of its molecular mechanisms and relationship with tumor immunity.