Objective: Parkinson's disease (PD) is the most common neurodegenerative movement disorder and immune-mediated mechanism is considered to be crucial to pathogenesis. Here, we investigated the role of humoral immune regulatory molecules in the pathogenesis of PD.Methods: Firstly, we performed a series of bioinformatic analyses utilizing the expression profile of the peripheral blood mononuclear cell (PBMC) obtained from the GEO database (GSE100054, GSE49126, and GSE22491) to identify differentially expressed genes related to humoral immune regulatory mechanisms between PD and healthy controls. Subsequently, we verified the results using quantitative polymerase chain reaction (Q-PCR) and enzymelinked immunosorbent assay (ELISA) in clinical blood specimen. Lastly, receiver operating characteristic (ROC) curve analysis was performed to determine the diagnostic effects of verified molecules.Results: We obtained 13 genes that were mainly associated with immunerelated biological processes in PD using bioinformatic analysis. Then, we selected PPBP, PROS1, and LCN2 for further exploration. Fascinatingly, our experimental results don't always coincide with the expression profile. PROS1 and LCN2 plasma levels were significantly higher in PD patients compared to controls (p < 0.01 and p < 0.0001). However, the PPBP plasma level and expression in the PBMC of PD patients was significantly decreased compared to controls (p < 0.01 and p < 0.01). We found that PPBP, PROS1, and LCN2 had an area under the curve (AUC) of 0.663 (95%CI: 0.551-0.776), 0.674 (95%CI: 0.569-0.780), and 0.885 (95%CI: 0.814-0.955). Furthermore, in the biological process analysis of gene ontology (GO), the three molecules were all involved in humoral immune response (GO:0006959).