Screening genes crucial for pediatric pilocytic astrocytoma using weighted gene coexpression network analysis combined with methylation data analysis

Abstract

To identify novel genes associated with pediatric pilocytic astrocytoma (PA) for better understanding the molecular mechanism underlying the pediatric PA pathogenesis. Gene expression profile data of GSE50161 and GSE44971 and the methylation data of GSE44684 were downloaded from Gene Expression Omnibus. The differentially expressed genes (DEGs) between PA and normal control samples were screened using the limma package in R, and then used to construct weighted gene coexpression network (WGCN) using the WGCN analysis (WGCNA) package in R. Significant modules of DEGs were selected using the clustering analysis. Function enrichment analysis of the DEGs in significant modules were performed using the WGCNA package and clusterprofiler package in R. Correlation between methylation sites of DEGs and PA was analyzed using the CpGassoc package in R. Totally, 3479 DEGs were screened in PA samples. Thereinto, 3424 DEGs were used to construct the WGCN. Several significant modules of DEGs were selected based on the WGCN, in which the turquoise module was positively related to PA, whereas blue module was negatively related to PA. DEGs (for example, DOCK2 (dedicator of cytokinesis 2), DOCK8 and FCGR2A (Fc fragment of IgG, low affinity IIa)) in blue module were mainly involved in Fc gamma R-mediated phagocytosis pathway and natural killer cell-mediated cytotoxicity pathway. Methylations of 14 DEGs among the top 30 genes in blue module were related to PA. Our data suggest that DOCK2, DOCK8 and FCGR2A may represent potential therapeutic targets in PA that merits further investigation.