Abstract
OBJECTIVE: To investigate the immune cell infiltration status in glioblastoma multiforme (GBM) and construct a novel prognostic risk model that can predict patients' prognosis. METHODS: The Cancer Genome Atlas (TCGA) database was used to obtain RNA-sequence information and relevant clinical data. We performed Pearson correlation, univariate Cox regression to screen m6A-related prognostic lncRNA. GMB patients' samples were separated into different clusters through the ConsensusClusterPlus package. The risk score model was established through LASSO regression analysis. Besides, KEGG pathway enrichment analysis was implemented. CIBERSORT algorithm was used to analyze the difference of 22 types of immune cell infiltration in different cluster of GBM patient. Cox regression analyses were used to verify the independence of the model and correlation analysis was performed to demonstrate the link between our model and clinical characteristics of GBM patients. Experiments were used to validate the differential expression of the model lncRNA in patients with different prognosis. RESULTS: 17 lncRNA related to prognosis were screened from 1021 m6A-related lncRNAs. Further, four m6A-related lncRNAs that were significantly correlated with GBM prognosis were selected to establish our prognostic risk model, which had excellent accuracy and can independently predict the prognosis of GBM patients. The infiltration fractions of T regulatory cells, T cells CD4 memory activated and neutrophils were positively associated with risk score, which suggested a significant relationship between the model and tumor immune microenvironment. CONCLUSION: The m6A-related RNA risk model offered potential for identifying biomarkers of therapy and predicting prognosis of GBM patients.