BackgroundAlthough increasing evidence has supported that Hirschsprung disease (HSCR) is the risk factor for children developing Crohn's disease (CD), the common mechanism of its co-occurrence remains unknown. The purpose of this study is to further explore the underlying mechanism and biomarkers for the co-occurrence of HSCR and CD. MethodsThe Gene Expression Omnibus (GEO) database was used to obtain gene expression profiles for CD (GSE95095) and HSCR (GSE98502). Following the identification of the shared differentially expressed genes (DEGs) of CD and HSCR, functional annotation, protein-protein interaction (PPI) network creation, and module assembly were performed to discover hub genes. RT-qPCR was performed to validate the expression of the hub genes in HSCR samples. The receiver operating characteristic (ROC) curve was utilized to assess the accuracy of the hub genes as biomarkers in predicting CD in both the training dataset and test dataset. ResultsA total of 103 common DEGs (50 downregulated genes and 53 upregulated genes) were chosen for further investigation. The importance of chemokines and cytokines in these two disorders is highlighted by functional analysis. MCODE plug identified three important modules, which functionally enriched the immune system process. Finally, nine hub genes were identified using cytoHubba, including IL1B, IL10, CXCL10, ICAM1, EGR1, FCGR3A, S100A12, S100A9, and FPR1. The nine hub genes were mainly enriched in immune- and inflammation-related pathways. External data profiles and RT-qPCR confirmed the expression of the nine hub genes in HSCR and CD. ROC analysis revealed that the nine hub genes had a strong diagnostic value. ConclusionOur study reveals the common pathogenesis of HSCR and CD. These hub genes and diagnostic models may provide novel insight for the diagnosis and treatment of HSCR complicated with CD.