Bioinformatics analysis of microarray data to reveal the pathogenesis of diffuse intrinsic pontine glioma

BACKGROUND: Diffuse intrinsic pontine glioma (DIPG) is the main cause of pediatric brain tumor death. This study was designed to identify key genes associated with DIPG. METHODS: The gene expression profile GSE50021, which consisted of 35 pediatric DIPG samples and 10 normal brain samples, was downloaded from the Gene Expression Omnibus database. Differentially expressed genes (DEGs) were identified by limma package. Functional and pathway enrichment analyses were performed by the DAVID tool. Protein–protein interaction (PPI) network, and transcription factor (TF)–microRNA (miRNA)–target gene network were constructed using Cytoscape. Moreover, the expression levels of several genes were validated in human glioma cell line U251 and normal glia HEB cells through real-time polymerase chain reaction (PCR). RESULTS: A total of 378 DEGs were screened (74 up-regulated and 304 down-regulated genes). In the PPI network, GRM1, HTR2A, GRM7 and GRM2 had higher degrees. Besides, GRM1 and HTR2A were significantly enriched in the neuroactive ligand–receptor interaction pathway, and calcium signaling pathway. In addition, TFAP2C was a significant down-regulated functional gene and hsa-miR-26b-5p had a higher degree in the TF-miRNA-target gene network. PCR analysis revealed that GRM7 and HTR2A were significantly downregulated while TFAP2C was upregulated in U251 cells compared with that in HEB cells (p < 0.001). GRM2 was not detected in cells. CONCLUSIONS: GRM1 and HTR2A might function in DIPG through the neuroactive ligand–receptor interaction pathway and the calcium signaling pathway. Furthermore, the TFAP2C and hsa-miR-26b-5p might play important roles in the development and progression mechanisms of DIPG.