Background: Preeclampsia is a form of hypertension in pregnancy, which induced by complicated factors. However, the pathogenesis of the disease is unclear. The present study was aimed to discover the critical biomarkers associated with the occurrence and development of preeclampsia. Methods:Gene data profile GSE75010 was downloaded from the Gene Expression Omnibus (GEO) database and used as discovery cohort to establish a WGCNA network determining significant modules which associated with clinical traits. Subsequently, functional enrichment analysis, pathway analysis and protein-protein interaction (PPI) network construction were performed on the core genes in significant modules to identify hub genes. Then, gene data profile GSE25906 was used as verified cohort to determine their diagnostic value of hub genes. The protein expression levels of these hub genes in preeclampsia and control placental tissues were verified using immunohistochemistry method. Finally, GSEA was performed to analyze their enrichment pathways. Results: Total 33 co-expression modules were identified after the establishment of WGCNA, of which 4 gene modules were identified as significant modules because they were related to multiple (>3) clinical traits. Total 75 core genes in significant modules were analyzed, and results showed that they were mainly enriched in adaptive immune response (Gene Ontology term) and platelet activation (Kyoto Encyclopedia of Genes and Genomes term). Finally, a total of 5 genes including TYROBP, PLEK, LCP2, HCK, ITGAM were identified as hub genes which scored high in PPI network and had high diagnostic value. Furthermore, the protein level of these 5 genes in placental tissues of preeclampsia was lower than that of the control group. Moreover, these 5 genes were all enriched in 17 pathways, including autoimmunity pathway. Conclusions:These 5 genes (TYROBP, PLEK, LCP2, HCK, ITGAM) may be closely related to the pathogenesis of preeclampsia, which may also help the diagnosis and therapy of preeclampsia.