High throughput miRNA sequencing and bioinformatics analysis identify the mesenchymal cell proliferation and apoptosis related miRNAs during fetal mice palate development

BACKGROUND: Palatogenesis requires a precise spatiotemporal regulation of gene expression. Recent studies indicate that microRNAs (miRNAs) are key factors in normal palatogenesis. The present study aimed to explain the regulatory mechanisms of miRNAs during palate development. METHODS: Pregnant ICR mice were choose at embryonic day 10.5 (E10.5). Hemotoxylin and eosin (H&E) staining was used to observe the morphological changes during the development of palatal process at embryonic day (E)13.5, E14.0, E14.5, E15.0 and E15.5. The fetal palatal tissues were collected at E13.5, E14.0, E14.5 and E15.0 to explore miRNA expression and function by high throughput sequencing and bioinformatic analysis. Mfuzz cluster analysis was used to look for miRNAs related to the fetal mice palate formation. The target genes of miRNAs were predicted by miRWalk. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis was performed base on target genes. The mesenchymal cell proliferation and apoptosis related miRNAs-genes networks were predicted and constructed using miRWalk and Cytoscape software. The expression of mesenchymal cell proliferation and apoptosis related miRNAs at the E13.5, E14.0, E14.5, and E15.0 was detected by a quantitative real-time PCR (RT-qPCR) assay. RESULTS: H&E staining found that the palatal process grows vertically along the sides of the tongue at E13.5, the position of the tongue begins to descend and the bilateral palatal processes rise above the tongue at E14.0, the palatal process grows horizontally at E14.5, there is palatal contact fusion at E15.0, and the palatal suture disappeared at E15.5. Nine clusters of miRNA expression changes were identified in the fetal mice palate formation progression, including two reducing trends, two rising trends and five disordered trends. Next, the heatmap showed the miRNA expression from Clusters 4, 6, 9, 12 in the E13.5, E14.0, E14.5 and E15.0 groups. GO functional and KEGG pathway enrichment analysis found target genes of miRNAs in clusters involved in regulation of mesenchymal phenotype and the mitogen-activated protein kinase (MAPK) signaling pathway. Next, mesenchymal phenotype related miRNA-genes networks were constructed. The heatmap showing that the mesenchymal phenotype related miRNA expression of Clusters 4, 6, 9 and 12 at E13.5, E14.0, E14.5 and E15.0. Furthermore, the mesenchymal cell proliferation and apoptosis related miRNA-gene networks were identified in Clusters 6 and 12, including mmu-miR-504-3p-Hnf1b, etc. The expression level of mesenchymal cell proliferation and apoptosis related miRNAs at the E13.5, E14.0, E14.5, and E15.0 was verified by a RT-qPCR assay. CONCLUSIONS: For the first time, we identified that clear dynamic miRNA expression during palate development. Furthermore, we demonstrated that mesenchymal cell proliferation and apoptosis related miRNAs, genes and the MAPK signaling pathway are important during fetal mice palate development.