Abstract:Craniofacial structures change dynamically in morphology during development through the coordinated regulation of various cellular molecules. However, it remains unclear how these complex mechanisms are regulated in a spatiotemporal manner. Here we applied natural cubic splines to model gene and microRNA (miRNA) expression from embryonic day (E) 10.5 to E14.5 in the proximal and distal regions of the maxillary processes to identify spatiotemporal patterns of gene and miRNA expression, followed by constructing corresponding regulatory networks. Three major groups of differentially expressed genes (DEGs) were identified, including 3,927 temporal, 314 spatial, and 494 spatiotemporal DEGs. Unsupervised clustering further resolved these spatiotemporal DEGs into 8 clusters with distinct expression patterns. Interestingly, we found 2 clusters of differentially expressed miRNAs: 1 had 80 miRNAs monotonically decreasing and the other had 97 increasing across developmental stages. To evaluate the phenotypic relevance of these DEGs during craniofacial development, we integrated data from the CleftGeneDB database and constructed the regulatory networks of genes related to orofacial clefts. Our analysis revealed 2 hub miRNAs, mmu-miR-325-3p and mmu-miR-384-5p, that repressed cleft-related genes Adamts3, Runx2, Fgfr2, Acvr1, and Edn2, while their expression increased over time. On the contrary, 2 hub miRNAs, mmu-miR-218-5p and mmu-miR-338-5p, repressed cleft-related genes Pbx2, Ermp1, Snai1, Tbx2, and Bmi1, while their expression decreased over time. Our experiments indicated that these miRNA mimics significantly inhibited cell proliferation in mouse embryonic palatal mesenchymal (MEPM) cells and O9-1 cells through the regulation of genes associated with cleft palate and validated the role of our regulatory networks in orofacial clefts. To facilitate interactive exploration of these data, we developed a user-friendly web tool to visualize the gene and miRNA expression patterns across developmental stages, as well as the regulatory networks (https://fyan.shinyapps.io/facebase_shiny/). Taken together, our results provide a valuable resource that serves as a reference map for future research in craniofacial development.

Expression analysis of genes including Zfhx4 in mice and zebrafish reveals a temporospatial conserved molecular basis underlying craniofacial development

Genome-wide Analysis of Facial Skeletal Regionalization in Zebrafish.

Zebrafish models of orofacial clefts

Transforming growth factor beta signaling and craniofacial development: modeling human diseases in zebrafish

Spatiotemporal MicroRNA-Gene Expression Network Related to Orofacial Clefts

Fox Proteins Are Modular Competency Factors For Facial Cartilage And Tooth Specification

Mouse Models of Orofacial Clefts: SHH and TGF-β Pathways

Identification of potential molecular mechanism related to craniofacial dysmorphism caused by FOXI3 deficiency

Single-cell Transcriptomic Signatures and Gene Regulatory Networks Modulated by Wls in Mammalian Midline Facial Formation and Clefts.

Duplicated zebrafish co-orthologs of parathyroid hormone-related peptide (PTHrP, Pthlh) play different roles in craniofacial skeletogenesis

Crucial and Overlapping Roles of Six1 and Six2 in Craniofacial Development.

Time-series expression profiles of mRNAs and lncRNAs during mammalian palatogenesis.

Examining craniofacial variation among crispant and mutant zebrafish models of human skeletal diseases

Using an aquatic model, Xenopus laevis, to uncover the role of chromodomain 1 in craniofacial disorders

Single-cell multiomics decodes regulatory programs for mouse secondary palate development

Integrative mRNA and miRNA Expression Profiles from Developing Zebrafish Head Highlight Brain-Preference Genes and Regulatory Networks

Functional and Genetic Analyses Unveil the Implication of CDC27 in Hemifacial Microsomia

A Shh-Foxf-Fgf18-Shh Molecular Circuit Regulating Palate Development

Fgfr2c~(C342Y/+) Mouse Model of Crouzon Syndrome in Craniosynostosis

Early embryogenesis in CHDFIDD mouse model reveals facial clefts and altered cranial neurogenesis

Genetic association and functional validation of ZFP36L2 in non-syndromic orofacial cleft subtypes