Abstract
The present study aimed to identify differentially regulated genes between the peritumoral brain zone (PBZ) and tumor core (TC) of glioblastoma (GBM), to elucidate the underlying molecular mechanisms and provide a target for the treatment of tumors. The GSE13276 and GSE116520 datasets were downloaded from the Gene Expression Omnibus (GEO) database. Differentially expressed genes (DEGs) for the PBZ and TC were obtained using the GEO2R tool. The bioinformatics and evolutionary genomics online tool Venn was used to identify common DEGs between the two datasets. The Database for Annotation, Visualization, and Integrated Discovery online tool was used to analyze enriched pathways of the Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) databases. The Search Tool for the Retrieval of Interacting Genes/Proteins online tool was used to construct a protein-protein interaction (PPI) network of DEGs. Hub genes were identified using Cytohubba, a plug-in for Cytoscape. The Gene Expression Profiling Interactive Analysis (GEPIA) database was utilized to perform survival analysis. In total, 75 DEGs, including 12 upregulated and 63 downregulated genes, were identified. In the GO term analysis, these DEGs were mainly enriched in 'regulation of angiogenesis' and 'central nervous system development'. Furthermore, in the KEGG pathway analysis, the DEGs were mainly enriched in 'bladder cancer' and 'endocytosis'. When filtering the results of the PPI network analysis using Cytohubba, a total of 10 hub genes, including proteolipid protein 1, myelin associated oligodendrocyte basic protein, contactin 2, myelin oligodendrocyte glycoprotein, myelin basic protein, myelin associated glycoprotein, SRY-box transcription factor 10, C-X-C motif chemokine ligand 8 (CXCL8), vascular endothelial growth factor A (VEGFA) and plasmolipin, were identified. These hub genes were further subjected to GO term and KEGG pathway analysis, and were revealed to be enriched in 'central nervous system development', 'bladder cancer' and 'rheumatoid arthritis'. These hub genes were used to perform survival analysis using the GEPIA database, and it was determined that VEGFA and CXCL8 were significantly associated with a reduction in the overall survival of patients with GBM. In conclusion, the results suggest that the recurrence of GBM is associated with high gene expression levels VEGFA and CXCL8, and the development of the central nervous system.