Identification of biomarkers and ceRNA network in glioblastoma through bioinformatic analysis and evaluation of potential prognostic values

BACKGROUND: Glioblastoma (GBM) is one of the most common and malignant primary brain tumors in adults, with high mortality rates and limited treatment. Based on bioinformatic analyses, this study aimed to identify biomarkers and relevant molecular pathways that may serve as potential targets for the treatment of GBM. METHODS: Expression profiles were downloaded from the National Center for Biotechnology Information (NCBI) Gene Expression Omnibus (GEO) database; nine GBM samples and three normal samples were extracted from the GSE104267 dataset. Differentially-expressed messenger RNA (mRNA) and long non-coding RNA (lncRNA) were screened from the preprocessed dataset. The clusterProfiler package in R was used to perform a biological process (BP) analysis of gene ontology (GO), and a Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis was performed separately in upregulated and downregulated groups. A competing endogenous RNA (ceRNA) network was constructed using Cytoscape. Based on data downloaded from The Cancer Genome Atlas (TCGA), Kaplan-Meier (K-M) survival curves were established. Real-time quantitative reverse transcription polymerase chain reaction (qRT-PCR) was performed to evaluate IL10RB antisense RNA 1 (IL10RB-AS1) expression in GBM tissue compared with that in normal brain tissue. RESULTS: A total of 253 differentially-expressed genes (DEGs) were obtained. Based on BP and KEGG enrichment annotation analyses, 11 lncRNA-related pathways were identified through function prediction analysis. A competing endogenous RNA (ceRNA) subnetwork, including 21 nodes and 29 regulatory pairs, was then constructed. Based on the clinical data of GBM in TCGA, one survival-related DEG, IL10RB-AS1, was identified using the log-rank statistical test. K-M survival curves of IL10RB-AS1 and expression levels of IL10RB-AS1 in both GBM and normal brain tissue were obtained. CONCLUSIONS: Through the combination of bioinformatic analyses, one survival-related differentially-expressed lncRNA, IL10RB-AS1, was identified. This, along with several related signaling pathways and ceRNA systems that were elucidated in GBM have potential prognostic value and might offer new possibilities for the treatment of GBM.