Abstract
Glucocerebrosidase (GBA) mutations occur frequently in Parkinson's disease (PD) patients. This study aims to identify potential crucial genes and pathways associated with GBA mutations in patients with PD and to further analyze new molecular mechanisms related to the occurrence of gene mutations from the perspective of bioinformatics. Gene expression profiles of datasets GSE53424 and GSE99142 were acquired from the Gene Expression Ominibus database. Differentially expressed genes (DEGs) were detected, using the 'limma' package in R, comparing IDI-PD 1 (idiopathic PD patients) and GBA-PD 1 [PD patients with heterozygous GBA mutations (GBA N370S)] group samples. The functions of top modules were assessed using the DAVID, whereas gene ontology and Kyoto Encyclopedia of Genes and Genomes pathway enrichment analyses were performed. Protein-protein interaction networks were assembled with Cytoscape software and separated into subnetworks using the Molecular Complex Detection Algorithm. Data from GSE53424 and GSE99142 were also extracted to verify our findings. There were 283 DEGs identified in PD patients heterozygous for GBA mutations. Module analysis revealed that GBA mutations in PD patients were associated with significant pathways, including Calcium signaling pathway, Rap1 signaling pathway and Cytokine-cytokine receptor interaction. Hub genes of the two modules were corticotropin-releasing hormone (CRH) and Melatonin receptor 1B (MTNR1B). The expression of CRH was downregulated, whereas that of MTNR1B was upregulated in PD patients with GBA mutations. The expression of CRH and MTNR1B has diagnostic value for PD patients with heterozygous GBA mutations. Novel DEGs and pathways identified herein might provide new insights into the underlying molecular mechanisms of heterozygous GBA mutations in PD patients.