Revealing the potential pathogenesis of glioma by utilizing a glioma associated protein-protein interaction network.

This study aims to explore the potential mechanism of glioma through bioinformatic approaches. The gene expression profile (GSE4290) of glioma tumor and non-tumor samples was downloaded from Gene Expression Omnibus database. A total of 180 samples were available, including 23 non-tumor and 157 tumor samples. Then the raw data were preprocessed using robust multiarray analysis, and 8,890 differentially expressed genes (DEGs) were identified by using t-test (false discovery rate < 0.0005). Furthermore, 16 known glioma related genes were abstracted from Genetic Association Database. After mapping 8,890 DEGs and 16 known glioma related genes to Human Protein Reference Database, a glioma associated protein-protein interaction network (GAPN) was constructed. In addition, 51 sub-networks in GAPN were screened out through Molecular Complex Detection (score ≥ 1), and sub-network 1 was found to have the closest interaction (score = 3). What' more, for the top 10 sub-networks, Gene Ontology (GO) enrichment analysis (p value < 0.05) was performed, and DEGs involved in sub-network 1 and 2, such as BRMS1L and CCNA1, were predicted to regulate cell growth, cell cycle, and DNA replication via interacting with known glioma related genes. Finally, the overlaps of DEGs and human essential, housekeeping, tissue-specific genes were calculated (p value = 1.0, 1.0, and 0.00014, respectively) and visualized by Venn Diagram package in R. About 61% of human tissue-specific genes were DEGs as well. This research shed new light on the pathogenesis of glioma based on DEGs and GAPN, and our findings might provide potential targets for clinical glioma treatment.

Gene co-expression network and function modules in three types of glioma.

The aim of the present study was to identify the disease­associated genes and their functions involved in the development of three types of glioma (astrocytoma, glioblastoma and oligodendroglioma) with DNA microarray technology, and to analyze their differences and correlations. First, the gene expression profile GSE4290 was downloaded from the Gene Expression Omnibus database, then the probe­level data were pre­processed and the differentially expressed genes (DEGs) were identified with limma package in R language. Gene functions of the selected DEGs were further analyzed with the Database for Annotation, Visualization and Integrated Discovery. After the co­expression network of DEGs was constructed by Cytoscape, the functional modules were mined and enrichment analysis was performed, and then the similarities and differences between any two types of glioma were compared. A total of 1151 genes between normal and astrocytoma tissues, 684 genes between normal and malignant glioma tissues, and 551 genes between normal and oligodendroglioma tissues were filtered as DEGs, respectively. By constructing co­expression networks of DEGs, a total of 77232, 455 and 987 interactions were involved in the differentially co­expressed networks of astrocytoma, oligodendroglioma and glioblastoma, respectively. The functions of DEGs were consistent with the modules in astrocytoma, glioblastoma and oligodendroglioma, which were mainly enriched in neuron signal transmission, immune responses and synthesis of organic acids, respectively. Model functions of astrocytoma and glioblastoma were similar (mainly related with immune response), while the model functions of oligodendroglioma differed markedly from that of the other two types. The identification of the associations among these three types of glioma has potential clinical utility for improving the diagnosis of different types of glioma in the future. In addition, these results have marked significance in studying the underlying mechanisms, distinguishing between normal and cancer tissues, and examining novel therapeutic strategies for patients with glioma.