Circadian rhythm disruption is typical in Parkinson’s disease (PD) early stage, and it plays an important role in the prognosis of the treatment effect in the advanced stage of PD. There is growing evidence that circadian rhythm genes can influence development of PD. Therefore, this study explored specific regulatory mechanism of circadian genes (C-genes) in PD through bioinformatic approaches.
Differentially expressed genes (DEGs) between PD and control samples were identified from GSE22491 using differential expression analysis. The key model showing the highest correlation with PD was derived through WGCNA analysis. Then, DEGs, 1,288 C-genes and genes in key module were overlapped for yielding differentially expressed C-genes (DECGs), and they were analyzed for LASSO and SVM-RFE for yielding critical genes. Meanwhile, from GSE22491 and GSE100054, receiver operating characteristic (ROC) was implemented on critical genes to identify biomarkers, and Gene Set Enrichment Analysis (GSEA) was applied for the purpose of exploring pathways involved in biomarkers. Eventually, immune infiltrative analysis was applied for understanding effect of biomarkers on immune microenvironment, and therapeutic drugs which could affect biomarkers expressions were also predicted. Finally, we verified the expression of the genes by q-PCR.
Totally 634 DEGs were yielded between PD and control samples, and MEgreen module had the highest correlation with PD, thus it was defined as key model. Four critical genes (AK3, RTN3, CYP4F2, and LEPR) were identified after performing LASSO and SVM-RFE on 18 DECGs. Through ROC analysis, AK3, RTN3, and LEPR were identified as biomarkers due to their excellent ability to distinguish PD from control samples. Besides, biomarkers were associated with Parkinson’s disease and other functional pathways.
Through bioinformatic analysis, the circadian rhythm related biomarkers were identified (AK3, RTN3 and LEPR) in PD, contributing to studies related to PD treatment.