Preliminary identification of key miRNAs, signaling pathways, and genes associated with Hirschsprung's disease by analysis of tissue microRNA expression profiles

Background: Hirschsprung's disease (HSCR) is a congenital gut motility disorder of infants, and if left untreated, it is fatal to the affected infants. This study aimed to identify key microRNAs (miRNAs), signaling pathways and genes involved in the pathogenesis of HSCR. Methods: The miRNA microarray dataset GSE77296 was downloaded. Nine colon tissue samples were available: six from HSCR patients and three matched control samples. Differentially expressed miRNAs (DEMs) were identified after data preprocessing. Target genes of the selected upregulated and downregulated DEMs were predicted. In addition, functional enrichment analyses for the selected DEMs and target genes were conducted. Finally, interaction networks between the DEMs and target genes were constructed. Results: A total of 162 DEMs (73 upregulated and 89 downregulated) were obtained. A total of 2511 DEM-target gene pairs for the 40 selected DEMs were identified, including 1645 pairs for the upregulated DEMs and 866 pairs for the downregulated DEMs. The upregulated DEM miR-141-3p and down-regulated DEM miR-30a-3p were identified as key miRNAs by Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment and network analyses. Besides, KEGG pathway enrichment analysis revealed that pathways in cancer and the mitogen-activated protein kinase (MAPK) signaling pathway were key pathways. The key genes frizzled class receptor 3 (FZD3) and docking protein 6 (DOK6) were obtained through the DEM-target gene interaction networks. Conclusion: Two key miRNAs (miR-141-3p and miR-30a-3p), the MAPK signaling pathway and two key genes (FZD3 and DOK6) were implicated in the pathogenesis of HSCR.