A computational approach based on weighted gene co-expression network analysis for biomarkers analysis of Parkinson’s disease and construction of diagnostic model

Background: Parkinson’s disease (PD) is a common age-related chronic neurodegenerative disease. There is currently no affordable, effective, and less invasive test for PD diagnosis. Metabolite profiling in blood and blood-based gene transcripts is thought to be an ideal method for diagnosing PD.

Aim: In this study, the objective is to identify the potential diagnostic biomarkers of PD by analyzing microarray gene expression data of samples from PD patients.

Methods: A computational approach, namely, Weighted Gene Co-expression Network Analysis (WGCNA) was used to construct co-expression gene networks and identify the key modules that were highly correlated with PD from the GSE99039 dataset. The Least Absolute Shrinkage and Selection Operator (LASSO) regression analysis was performed to identify the hub genes in the key modules with strong association with PD. The selected hub genes were then used to construct a diagnostic model based on logistic regression analysis, and the Receiver Operating Characteristic (ROC) curves were used to evaluate the efficacy of the model using the GSE99039 dataset. Finally, Reverse Transcription-Polymerase Chain Reaction (RT-PCR) was used to validate the hub genes.

Results: WGCNA identified two key modules associated with inflammation and immune response. Seven hub genes, LILRB1, LSP1, SIPA1, SLC15A3, MBOAT7, RNF24, and TLE3 were identified from the two modules and used to construct diagnostic models. ROC analysis showed that the diagnostic model had a good diagnostic performance for PD in the training and testing datasets. Results of the RT-PCR experiments showed that there were significant differences in the mRNA expression of LILRB1, LSP1, and MBOAT7 among the seven hub genes.