Identification of Core Biomarkers Associated with Outcome in Glioma: Evidence from Bioinformatics Analysis

Glioma is the most common neoplasm of the central nervous system (CNS); the progression and outcomes of which are affected by a complicated network of genes and pathways. We chose a gene expression profile of GSE66354 from GEO database to search core biomarkers during the occurrence and development of glioma. A total of 149 samples, involving 136 glioma and 13 normal brain tissues, were enrolled in this article. 1980 differentially expressed genes (DEGs) including 697 upregulated genes and 1283 downregulated genes between glioma patients and healthy individuals were selected using GeoDiver and GEO2R tool. Then, gene ontology (GO) analysis as well as Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis were carried out using the Database for Annotation, Visualization and Integrated Discovery (DAVID). Moreover, Cytoscape with Search Tool for the Retrieval of Interacting Genes (STRING) and Molecular Complex Detection (MCODE) plug-in was employed to imagine protein-protein interaction (PPI) of these DEGs. The upregulated genes were enriched in cell cycle, ECM-receptor interaction, and p53 signaling pathway, while the downregulated genes were enriched in retrograde endocannabinoid signaling, glutamatergic synapse, morphine addiction, GABAergic synapse, and calcium signaling pathway. Subsequently, 4 typical modules were discovered by the PPI network utilizing MCODE software. Besides, 15 hub genes were chosen according to the degree of connectivity, including TP53, CDK1, CCNB1, and CCNB2, the Kaplan-Meier analysis of which was further identified. In conclusion, this bioinformatics analysis indicated that DEGs and core genes, such as TP53, might influence the development of glioma, especially in tumor proliferation, which were expected to be promising biomarkers for diagnosis and treatment of glioma.