Background: Kawasaki disease (KD) is a pediatric inflammatory disorder causes coronary artery complications. The disease overlapping manifestations with a set of symptomatically like diseases such as bacterial and viral infections, juvenile idiopathic arthritis, Henoch-Schönlein purpura, infection of unknown etiology, group-A streptococcal and adenoviral infections, and incomplete KD could lead to misdiagnosis of the disease. Methods: In the present study, we applied weighted gene co-expression network analysis (WGCNA) to identify network modules of co-expressed genes in GSE73464 and also, limma package was used to identify the differentially expressed genes (DEGs) in KD expression arrays composed of GSE73464, GSE18606, GSE109351, and GSE68004. By merging the results of WGCNA and limma, we detected hub genes. Then, analyzed the peripheral blood mononuclear cells (PBMCs) of 16 patients and 8 control subjects using Real-Time Polymerase Chain Reaction (RT-PCR) to evaluate the previous results. Results: We assessed the diagnostic potency of the screened genes by plotting the area under curve (AUC). We finally identified 2 genes CASP5(Caspase 5) and CR1(Complement C3b/C4b Receptor 1) which were shown to potentially discriminate KD from other similar diseases and also from healthy people. Conclusions: The results of RT-PCR and AUC confirmed the diagnostic potentials of two suggested biomarkers for KD.