The old and new face of craniofacial research: How animal models inform human craniofacial genetic and clinical data

Dyrk1a is required for craniofacial development in Xenopus laevis

Loss of function variations in the dual specificity tyrosine-phosphorylation-regulated kinase 1 A (DYRK1A) gene are associated with craniofacial malformations in humans. Here we characterized the effects of deficient DYRK1A in craniofacial development using a developmental model, Xenopus laevis. Dyrk1a mRNA and protein were expressed throughout the developing head and both were enriched in the branchial arches which contribute to the face and jaw. Consistently, reduced Dyrk1a function, using dyrk1a morpholinos and pharmacological inhibitors, resulted in orofacial malformations including hypotelorism, altered mouth shape, slanted eyes, and narrower face accompanied by smaller jaw cartilage and muscle. Inhibition of Dyrk1a function resulted in misexpression of key craniofacial regulators including transcription factors and members of the retinoic acid signaling pathway. Two such regulators, sox9 and pax3 are required for neural crest development and their decreased expression corresponds with smaller neural crest domains within the branchial arches. Finally, we determined that the smaller size of the faces, jaw elements and neural crest domains in embryos deficient in Dyrk1a could be explained by increased cell death and decreased proliferation. This study is the first to provide insight into why craniofacial birth defects might arise in humans with variants of DYRK1A.

A review of the reproductive biology of mormyroid fishes: An emerging model for biomedical research

Mormyroidea is a superfamily of weakly electric African fishes with great potential as a model in a variety of biomedical research areas including systems neuroscience, muscle cell and craniofacial development, ion channel biophysics, and flagellar/ciliary biology. However, they are currently difficult to breed in the laboratory setting, which is essential for any tractable model organism. As such, there is a need to better understand the reproductive biology of mormyroids to breed them more reliably in the laboratory to effectively use them as a biomedical research model. This review seeks to (1) briefly highlight the biomedically relevant phenotypes of mormyroids and (2) compile information about mormyroid reproduction including sex differences, breeding season, sexual maturity, gonads, gametes, and courtship/spawning behaviors. We also highlight areas of mormyroid reproductive biology that are currently unexplored and/or have the potential for further investigation that may provide insights into more successful mormyroid laboratory breeding methods.

Craniofacial chondrogenesis in organoids from human stem cell-derived neural crest cells

Knowledge of cell signaling pathways that drive human neural crest differentiation into craniofacial chondrocytes is incomplete, yet essential for using stem cells to regenerate craniomaxillofacial structures. To accelerate translational progress, we developed a differentiation protocol that generated self-organizing craniofacial cartilage organoids from human embryonic stem cell-derived neural crest stem cells. Histological staining of cartilage organoids revealed tissue architecture and staining typical of elastic cartilage. Protein and post-translational modification (PTM) mass spectrometry and snRNA-seq data showed that chondrocyte organoids expressed robust levels of cartilage extracellular matrix (ECM) components: many collagens, aggrecan, perlecan, proteoglycans, and elastic fibers. We identified two populations of chondroprogenitor cells, mesenchyme cells and nascent chondrocytes, and the growth factors involved in paracrine signaling between them. We show that ECM components secreted by chondrocytes not only create a structurally resilient matrix that defines cartilage, but also play a pivotal autocrine cell signaling role in determining chondrocyte fate.

Selection Signatures Reveal Candidate Genes for the Cornish Rex Breed-Specific Phenotype

Many coat color, behavioral and morphological traits are specific and fixed across cat breeds, with several variants influencing these traits being common among different breeds. In the domestic cat, rexoid mutations have been documented in several breeds. In the Cornish Rex, four bp deletion in the LPAR6 gene has been found to cause a frame shift and a premature stop codon. In addition to the rexoid coat, Cornish Rex cats also have a characteristic head, ear shape and body type. Analysis of the selection signatures in the Cornish Rex genome revealed several regions that are under selective pressure. One of these is located in CFA B4, in the region where the ALX1 gene is located. The ALX1 gene in Burmese cats disrupts the cranial morphogenesis and causes brachycephaly in the heterozygous state. In our study, we confirmed the presence of a deletion in LPAR6 in 20 Cornish Rex and in four F1 hybrids between Cornish Rex and domestic cat. However, we did not confirm the presence of the deletion in ALX1 in Cornish Rex cats. Genome-wide selection signature analysis was performed using ROH islands and integrated haplotype score (iHS) statistics based on publicly available SNP array data of 11 Cornish Rex cats. The selection signatures were detected on chromosomes A1, A3, C2, B1, B4 and D1.

Dynamic enhancer landscapes in human craniofacial development

The genetic basis of human facial variation and craniofacial birth defects remains poorly understood. Distant-acting transcriptional enhancers control the fine-tuned spatiotemporal expression of genes during critical stages of craniofacial development. However, a lack of accurate maps of the genomic locations and cell type-resolved activities of craniofacial enhancers prevents their systematic exploration in human genetics studies. Here, we combine histone modification, chromatin accessibility, and gene expression profiling of human craniofacial development with single-cell analyses of the developing mouse face to define the regulatory landscape of facial development at tissue- and single cell-resolution. We provide temporal activity profiles for 14,000 human developmental craniofacial enhancers. We find that 56% of human craniofacial enhancers share chromatin accessibility in the mouse and we provide cell population- and embryonic stage-resolved predictions of their in vivo activity. Taken together, our data provide an expansive resource for genetic and developmental studies of human craniofacial development.

Dyrk1a is required for craniofacial development in Xenopus laevis

Loss of function mutations in the dual specificity tyrosine-phosphorylation-regulated kinase 1A (DYRK1A) gene are associated with craniofacial malformations in humans. Here we characterized the effects of deficient DYRK1A in craniofacial development using a developmental model, Xenopus laevis . Dyrk1a mRNA and protein was expressed throughout the developing head and was enriched in the branchial arches which contribute to the face and jaw. Consistently, reduced Dyrk1a function, using dyrk1a morpholinos and pharmacological inhibitors, resulted in orofacial malformations including hypotelorism, altered mouth shape, slanted eyes, and narrower face accompanied by smaller jaw cartilage and muscle. Inhibition of Dyrk1a function resulted in misexpression of key craniofacial regulators including transcription factors and members of the retinoic acid signaling pathway. Two such regulators, sox9 and pax3 are required for neural crest development and their decreased expression corresponds with smaller neural crest domains within the branchial arches. Finally, we determined that the smaller size of the faces, jaw elements and neural crest domains in embryos deficient in Dyrk1a could be explained by increased cell death and decreased proliferation. This study is the first to provide insight into why craniofacial birth defects might arise in humans with DYRK1A mutations.

Reflecting the human lip in vitro: Cleft lip skin and mucosa keratinocytes keep their identities

OBJECTIVES

Cell models have shown great promise as tools for research, potentially providing intriguing alternatives to animal models. However, the original tissue characteristics must be maintained in culture, a fact that is often assumed, but seldom assessed. We aimed to follow the retention of the original tissue identities of cleft lip-derived skin and mucosa keratinocytes in vitro.

METHODS

Cleft lip-derived keratinocytes were isolated from discarded tissue along the cleft margins during cheiloplasty. Cell identities were assessed by immunohistochemistry and quantitative real-time PCR for tissue-specific markers and compared with native lip tissue. Moreover, keratinocytes were regularly analyzed for the retention of the original tissue characteristics by the aforementioned methods as well as by differentiation assays.

RESULTS

The various anatomical zones of the human lip could be distinguished using a panel of differentiation and functional-based markers. Using these markers, retention of the original tissue identities could be followed and confirmed in the corresponding primary keratinocytes in culture.

CONCLUSIONS

Our findings promote patient-derived cells retaining their original identities as astonishing and clinically relevant in vitro tools. Such cells allow a better molecular understanding of various lip-associated pathologies as well as their modeling in vitro, including but not restricted to orofacial clefts.

Validation of a zebrafish developmental defects assay as a qualified alternative test for its regulatory use following the ICH S5(R3) guideline

Zebrafish is a popular toxicology model and provides an ethically acceptable small-scale analysis system with the complexity of a complete organism. Our goal is to further validate this model for its regulatory use for reproductive and developmental defects by testing the compounds indicated in the "Guideline on detection of reproductive and developmental toxicity for human pharmaceuticals" (ICH S5(R3) guideline.) To determine the embryotoxic and developmental risk of the 32 reference compounds listed in the ICH S5(R3) guideline, the presence of morphological alterations in zebrafish embryos was analyzed at two different stages to calculateLC50 and EC50 values for each stage. Teratogenic Indexes were established as the ratio between LC50 and EC50 critical for the proper compound classification as teratogenic when it is ≥ 2. A total of three biological replicates have been conducted to study the reproducibility of the assay. The chemicals' concentration in the medium and internally in the zebrafish embryos was evaluated. In this study, the 3 negative compounds were properly categorized while 23 compounds out of the 29 reference ones (sensitivity of 79.31%) were classified as teratogenic in zebrafish. The 6 that had false-negative results were classified 4 as inconclusive, 1 as not toxic, and 1 compound resulted toxic for zebrafish embryos under testing conditions. After the bioavailability experiments, some of the obtained inconclusive results were refined. The developmental defects assay in zebrafish gives an accuracy of 89.66%, sensitivity of 88.46%, specificity and repeatability of 100% compared to mammals; therefore, this is a well-integrated strategy using New Alternative Methods, to minimize the use of animals in developmental toxicity studies.

Craniofacial developmental biology in the single-cell era

The evolution of a unique craniofacial complex in vertebrates made possible new ways of breathing, eating, communicating and sensing the environment. The head and face develop through interactions of all three germ layers, the endoderm, ectoderm and mesoderm, as well as the so-called fourth germ layer, the cranial neural crest. Over a century of experimental embryology and genetics have revealed an incredible diversity of cell types derived from each germ layer, signaling pathways and genes that coordinate craniofacial development, and how changes to these underlie human disease and vertebrate evolution. Yet for many diseases and congenital anomalies, we have an incomplete picture of the causative genomic changes, in particular how alterations to the non-coding genome might affect craniofacial gene expression. Emerging genomics and single-cell technologies provide an opportunity to obtain a more holistic view of the genes and gene regulatory elements orchestrating craniofacial development across vertebrates. These single-cell studies generate novel hypotheses that can be experimentally validated in vivo. In this Review, we highlight recent advances in single-cell studies of diverse craniofacial structures, as well as potential pitfalls and the need for extensive in vivo validation. We discuss how these studies inform the developmental sources and regulation of head structures, bringing new insights into the etiology of structural birth anomalies that affect the vertebrate head.

Cell Type- and Tissue-specific Enhancers in Craniofacial Development

The genetic basis of craniofacial birth defects and general variation in human facial shape remains poorly understood. Distant-acting transcriptional enhancers are a major category of non-coding genome function and have been shown to control the fine-tuned spatiotemporal expression of genes during critical stages of craniofacial development1-3. However, a lack of accurate maps of the genomic location and cell type-specific in vivo activities of all craniofacial enhancers prevents their systematic exploration in human genetics studies. Here, we combined histone modification and chromatin accessibility profiling from different stages of human craniofacial development with single-cell analyses of the developing mouse face to create a comprehensive catalogue of the regulatory landscape of facial development at tissue- and single cell-resolution. In total, we identified approximately 14,000 enhancers across seven developmental stages from weeks 4 through 8 of human embryonic face development. We used transgenic mouse reporter assays to determine the in vivo activity patterns of human face enhancers predicted from these data. Across 16 in vivo validated human enhancers, we observed a rich diversity of craniofacial subregions in which these enhancers are active in vivo. To annotate the cell type specificities of human-mouse conserved enhancers, we performed single-cell RNA-seq and single-nucleus ATAC-seq of mouse craniofacial tissues from embryonic days e11.5 to e15.5. By integrating these data across species, we find that the majority (56%) of human craniofacial enhancers are functionally conserved in mice, providing cell type- and embryonic stage-resolved predictions of their in vivo activity profiles. Using retrospective analysis of known craniofacial enhancers in combination with single cell-resolved transgenic reporter assays, we demonstrate the utility of these data for predicting the in vivo cell type specificity of enhancers. Taken together, our data provide an expansive resource for genetic and developmental studies of human craniofacial development.

Identification of de novo variants in nonsyndromic cleft lip with/without cleft palate patients with low polygenic risk scores

BACKGROUND

Nonsyndromic cleft lip with/without cleft palate (nsCL/P) is a congenital malformation of multifactorial etiology. Research has identified >40 genome-wide significant risk loci, which explain less than 40% of nsCL/P heritability. Studies show that some of the hidden heritability is explained by rare penetrant variants.

METHODS

To identify new candidate genes, we searched for highly penetrant de novo variants (DNVs) in 50 nsCL/P patient/parent-trios with a low polygenic risk for the phenotype (discovery). We prioritized DNV-carrying candidate genes from the discovery for resequencing in independent cohorts of 1010 nsCL/P patients of diverse ethnicities and 1574 population-matched controls (replication). Segregation analyses and rare variant association in the replication cohort, in combination with additional data (genome-wide association data, expression, protein-protein-interactions), were used for final prioritization.

CONCLUSION

In the discovery step, 60 DNVs were identified in 60 genes, including a variant in the established nsCL/P risk gene CDH1. Re-sequencing of 32 prioritized genes led to the identification of 373 rare, likely pathogenic variants. Finally, MDN1 and PAXIP1 were prioritized as top candidates. Our findings demonstrate that DNV detection, including polygenic risk score analysis, is a powerful tool for identifying nsCL/P candidate genes, which can also be applied to other multifactorial congenital malformations.

Micro-CT analysis reveals the changes in bone mineral density in zebrafish craniofacial skeleton with age

Bone has multiple functions in animals, such as supporting the body for mobility. The zebrafish skeleton is composed of craniofacial and axial skeletons. It shares a physiological curvature and consists of a similar number of vertebrae as humans. Bone degeneration and malformations have been widely studied in zebrafish as human disease models. High-resolution imaging and different bone properties such as density and volume can be obtained using micro-computed tomography (micro-CT). This study aimed to understand the possible changes in the structure and bone mineral density (BMD) of the vertebrae and craniofacial skeleton with age (4, 12 and 24 months post fertilisation [mpf]) in zebrafish. Our data showed that the BMD in the vertebrae and specific craniofacial skeleton (mandibular arch, ceratohyal and ethmoid plate) of 12 and 24 mpf fish were higher than that of the 4 mpf fish. In addition, we found the age-dependent increase in BMD was not ubiquitously observed in facial bones, and such differences were not correlated with bone type. In summary, such additional information on the craniofacial skeleton could help in understanding bone development throughout the lifespan of zebrafish.

Treatment with an inhibitor of matrix metalloproteinase 9 or cathepsin K lengthens embryonic lower jaw bone

OBJECTIVES

Skeletal malocclusions are common, and severe malocclusions are treated by invasive surgeries. Recently, jaw bone length has been shown to be developmentally controlled by osteoclasts. Our objective was to determine the effect of inhibiting osteoclast-secreted proteolytic enzymes on lower jaw bone length of avian embryos by pharmacologically inhibiting matrix metalloproteinase-9 (MMP9) or cathepsin K (CTSK).

METHODS

Quail (Coturnix coturnix japonica) embryos were given a single dose of an inhibitor of MMP9 (iMMP9), an inhibitor CTSK (iCTSK), or vehicle at a developmental stage when bone deposition is beginning to occur. At a developmental stage when the viscerocranium is largely calcified, the heads were scanned via micro-computed tomography and reproducible landmarks were placed on 3D-reconstructed skulls; the landmark coordinates were used to quantify facial bone dimensions.

RESULTS

Approximately half of the quail given either iMMP9 or iCTSK demonstrated an overt lower jaw phenotype, characterized by longer lower jaw bones and a greater lower to upper jaw ratio than control embryos. Additionally, iMMP9-treated embryos exhibited a significant change in midface length and iCTSK-treated embryos had significant change in nasal bone length.

CONCLUSION

MMP9 and CTSK play a role in osteoclast-mediated determination of lower jaw bone length. Pharmacological inhibition of MMP9 or CTSK may be a promising therapeutic alternative to surgery for treating skeletal jaw malocclusions, but more preclinical research is needed prior to clinical translation.

Genome-wide analysis of copy-number variation in humans with cleft lip and/or cleft palate identifies COBLL1, RIC1, and ARHGEF38 as clefting genes

Cleft lip with or without cleft palate (CL/P) is a common birth defect with a complex, heterogeneous etiology. It is well established that common and rare sequence variants contribute to the formation of CL/P, but the contribution of copy-number variants (CNVs) to cleft formation remains relatively understudied. To fill this knowledge gap, we conducted a large-scale comparative analysis of genome-wide CNV profiles of 869 individuals from the Philippines and 233 individuals of European ancestry with CL/P with three primary goals: first, to evaluate whether differences in CNV number, amount of genomic content, or amount of coding genomic content existed within clefting subtypes; second, to assess whether CNVs in our cohort overlapped with known Mendelian clefting loci; and third, to identify unestablished Mendelian clefting genes. Significant differences in CNVs across cleft types or in individuals with non-syndromic versus syndromic clefts were not observed; however, several CNVs in our cohort overlapped with known syndromic and non-syndromic Mendelian clefting loci. Moreover, employing a filtering strategy relying on population genetics data that rare variants are on the whole more deleterious than common variants, we identify several CNV-associated gene losses likely driving non-syndromic clefting phenotypes. By prioritizing genes deleted at a rare frequency across multiple individuals with clefts yet enriched in our cohort of individuals with clefts compared to control subjects, we identify COBLL1, RIC1, and ARHGEF38 as clefting genes. CRISPR-Cas9 mutagenesis of these genes in Xenopus laevis and Danio rerio yielded craniofacial dysmorphologies, including clefts analogous to those seen in human clefting disorders.

Acetaminophen Disrupts the Development of Pharyngeal Arch-Derived Cartilage and Muscle in Zebrafish

Acetaminophen is a common analgesic, but its potential effects on early embryonic development are not well understood. Previous studies using zebrafish (Danio rerio) have described the effects of acetaminophen on liver development and physiology, and a few have described gross physiological and morphological defects. Using a high but non-embryonic lethal dose of acetaminophen, we probed for defects in zebrafish craniofacial cartilage development. Strikingly, acetaminophen treatment caused severe craniofacial cartilage defects, primarily affecting both the presence and morphology of pharyngeal arch-derived cartilages of the viscerocranium. Delaying acetaminophen treatment restored developing cartilages in an order correlated with their corresponding pharyngeal arches, suggesting that acetaminophen may target pharyngeal arch development. Craniofacial cartilages are derived from cranial neural crest cells; however, many neural crest cells were still seen along their expected migration paths, and most remaining cartilage precursors expressed the neural crest markers sox9a and sox10, then eventually col2a1 (type II collagen). Therefore, the defects are not primarily due to an early breakdown of neural crest or cartilage differentiation. Instead, apoptosis is increased around the developing pharyngeal arches prior to chondrogenesis, further suggesting that acetaminophen may target pharyngeal arch development. Many craniofacial muscles, which develop in close proximity to the affected cartilages, were also absent in treated larvae. Taken together, these results suggest that high amounts of acetaminophen can disrupt multiple aspects of craniofacial development in zebrafish.

16p12.1 Deletion Orthologs are Expressed in Motile Neural Crest Cells and are Important for Regulating Craniofacial Development in Xenopus laevis

Copy number variants (CNVs) associated with neurodevelopmental disorders are characterized by extensive phenotypic heterogeneity. In particular, one CNV was identified in a subset of children clinically diagnosed with intellectual disabilities (ID) that results in a hemizygous deletion of multiple genes at chromosome 16p12.1. In addition to ID, individuals with this deletion display a variety of symptoms including microcephaly, seizures, cardiac defects, and growth retardation. Moreover, patients also manifest severe craniofacial abnormalities, such as micrognathia, cartilage malformation of the ears and nose, and facial asymmetries; however, the function of the genes within the 16p12.1 region have not been studied in the context of vertebrate craniofacial development. The craniofacial tissues affected in patients with this deletion all derive from the same embryonic precursor, the cranial neural crest, leading to the hypothesis that one or more of the 16p12.1 genes may be involved in regulating neural crest cell (NCC)-related processes. To examine this, we characterized the developmental role of the 16p12.1-affected gene orthologs, polr3e, mosmo, uqcrc2, and cdr2, during craniofacial morphogenesis in the vertebrate model system, Xenopus laevis. While the currently-known cellular functions of these genes are diverse, we find that they share similar expression patterns along the neural tube, pharyngeal arches, and later craniofacial structures. As these genes show co-expression in the pharyngeal arches where NCCs reside, we sought to elucidate the effect of individual gene depletion on craniofacial development and NCC migration. We find that reduction of several 16p12.1 genes significantly disrupts craniofacial and cartilage formation, pharyngeal arch migration, as well as NCC specification and motility. Thus, we have determined that some of these genes play an essential role during vertebrate craniofacial patterning by regulating specific processes during NCC development, which may be an underlying mechanism contributing to the craniofacial defects associated with the 16p12.1 deletion.

Mouse models in palate development and orofacial cleft research: Understanding the crucial role and regulation of epithelial integrity in facial and palate morphogenesis

Cleft lip and cleft palate are common birth defects resulting from genetic and/or environmental perturbations of facial development in utero. Facial morphogenesis commences during early embryogenesis, with cranial neural crest cells interacting with the surface ectoderm to form initially partly separate facial primordia consisting of the medial and lateral nasal prominences, and paired maxillary and mandibular processes. As these facial primordia grow around the primitive oral cavity and merge toward the ventral midline, the surface ectoderm undergoes a critical differentiation step to form an outer layer of flattened and tightly connected periderm cells with a non-stick apical surface that prevents epithelial adhesion. Formation of the upper lip and palate requires spatiotemporally regulated inter-epithelial adhesions and subsequent dissolution of the intervening epithelial seam between the maxillary and medial/lateral nasal processes and between the palatal shelves. Proper regulation of epithelial integrity plays a paramount role during human facial development, as mutations in genes encoding epithelial adhesion molecules and their regulators have been associated with syndromic and non-syndromic orofacial clefts. In this chapter, we summarize mouse genetic studies that have been instrumental in unraveling the mechanisms regulating epithelial integrity and periderm differentiation during facial and palate development. Since proper epithelial integrity also plays crucial roles in wound healing and cancer, understanding the mechanisms regulating epithelial integrity during facial development have direct implications for improvement in clinical care of craniofacial patients.

E-liquids and vanillin flavoring disrupts retinoic acid signaling and causes craniofacial defects in Xenopus embryos

Environmental teratogens such as smoking are known risk factors for developmental disorders such as cleft palate. While smoking rates have declined, a new type of smoking, called vaping is on the rise. Vaping is the use of e-cigarettes to vaporize and inhale an e-liquid containing nicotine and food-like flavors. There is the potential that, like smoking, vaping could also pose a danger to the developing human. Rather than waiting for epidemiological and mammalian studies, we have turned to an aquatic developmental model, Xenopus laevis, to more quickly assess whether e-liquids contain teratogens that could lead to craniofacial malformations. Xenopus, like zebrafish, has the benefit of being a well-established developmental model and has also been effective in predicting whether a chemical could be a teratogen. We have determined that embryonic exposure to dessert flavored e-liquids can cause craniofacial abnormalities, including an orofacial cleft in Xenopus. To better understand the underlying mechanisms contributing to these defects, transcriptomic analysis of the facial tissues of embryos exposed to a representative dessert flavored e-liquid vapor extract was performed. Analysis of differentially expressed genes in these embryos revealed several genes associated with retinoic acid metabolism or the signaling pathway. Consistently, retinoic acid receptor inhibition phenocopied the craniofacial defects as those embryos exposed to the vapor extract of the e-liquid. Such malformations also correlated with a group of common differentially expressed genes, two of which are associated with midface birth defects in humans. Further, e-liquid exposure sensitized embryos to forming craniofacial malformations when they already had depressed retinoic acid signaling. Moreover, 13-cis-retinoic acid treatment could significantly reduce the e-liquid induced malformation in the midface. Such results suggest the possibility of an interaction between retinoic acid signaling and e-liquid exposure. One of the most popular and concentrated flavoring chemicals in dessert flavored e-liquids is vanillin. Xenopus embryos exposed to this chemical closely resembled embryos exposed to dessert-like e-liquids and a retinoic acid receptor antagonist. In summary, we determined that e-liquid chemicals, in particular vanillin, can cause craniofacial defects potentially by dysregulating retinoic acid signaling. This work warrants the evaluation of vanillin and other such flavoring additives in e-liquids on mammalian development.

Zebrafish models of fetal alcohol spectrum disorders

Fetal alcohol spectrum disorder (FASD) describes a wide range of structural deficits and cognitive impairments. FASD impacts up to 5% of children born in the United States each year, making ethanol one of the most common teratogens. Due to limitations and ethical concerns, studies in humans are limited in their ability to study FASD. Animal models have proven critical in identifying and characterizing the mechanisms underlying FASD. In this review, we will focus on the attributes of zebrafish that make it a strong model in which to study ethanol-induced developmental defects. Zebrafish have several attributes that make it an ideal model in which to study FASD. Zebrafish produced large numbers of externally fertilized, translucent embryos. With a high degree of genetic amenability, zebrafish are at the forefront of identifying and characterizing the gene-ethanol interactions that underlie FASD. Work from multiple labs has shown that embryonic ethanol exposures result in defects in craniofacial, cardiac, ocular, and neural development. In addition to structural defects, ethanol-induced cognitive and behavioral impairments have been studied in zebrafish. Building upon these studies, work has identified ethanol-sensitive loci that underlie the developmental defects. However, analyses show there is still much to be learned of these gene-ethanol interactions. The zebrafish is ideally suited to expand our understanding of gene-ethanol interactions and their impact on FASD. Because of the conservation of gene function between zebrafish and humans, these studies will directly translate to studies of candidate genes in human populations and allow for better diagnosis and treatment of FASD.

Factors Affecting Outcomes of Craniofacial Clinical Trials

PURPOSE

This study aims to provide a cross-sectional view of craniofacial trials and to identify factors associated with completion, publication, and trial longevity.

MATERIALS AND METHODS

This is a cross-sectional study of cleft and craniofacial clinical trials registered with ClinicalTrials.gov between September 1999 and April 2020. Predictor variables included funding source, trial design, trial location, number of recruitment sites, and investigator specialty.Study outcomes were completion status, publication status, and trial duration. Univariate comparisons and multivariate regression models were calculated for each outcome.

RESULTS

The final sample included 179 clinical trials pertaining to craniofacial care. Nearly all trials were single-center (86.5%), and roughly half of trials were interventional (57.0%) or conducted in the United States (40.5%). No single specialty predominated, although plastic surgery (13.4%) was the most common investigator specialty. The completion rate was 82.7%, the publication rate was 40.8%, and the mean trial duration was 39.1 months. Interventional design (odds ratio = 0.30, P = 0.02) and United States location (odds ratio = 0.15, P < 0.01) were each independently associated with lower odds of trial completion. Trial longevity was independently associated with the National Institute of Health-funding (P < 0.01) and multicenter design (P < 0.01).

CONCLUSIONS

Craniofacial trials are multidisciplinary and have a high rate of completion. Although most existing trials were conducted at only a single-center, multicenter efforts significantly increased trial longevity without compromising completion and publication rates. Given the diverse array of conditions and lines of inquiry that compose craniofacial care, it is reassuring that collaboration did not negatively affect trial outcomes.

A Living Cell Repository of the Cranio-/Orofacial Region to Advance Research and Promote Personalized Medicine

The prevalence of congenital anomalies in newborns is estimated to be as high as 6%, many of which involving the cranio-/orofacial region. Such malformations, including several syndromes, are usually identified prenatally, at birth, or rarely later in life. The lack of clinically relevant human cell models of these often very rare conditions, the societal pressure to avoid the use of animal models and the fact that the biological mechanisms between rodents and human are not necessarily identical, makes studying cranio-/orofacial anomalies challenging. To overcome these limitations, we are developing a living cell repository of healthy and diseased cells derived from the cranio-/orofacial region. Ultimately, we aim to make patient-derived cells, which retain the molecular and genetic characteristics of the original anomaly or disease in vitro, available for the scientific community. We report our efforts in establishing a human living cell bank derived from the cranio-/orofacial region of otherwise discarded tissue samples, detail our strategy, processes and quality checks. Such specific cell models have a great potential for discovery and translational research and might lead to a better understanding and management of craniofacial anomalies for the benefit of all affected individuals.

An Irf6-Esrp1/2 regulatory axis controls midface morphogenesis in vertebrates

Irf6 and Esrp1 are important for palate development across vertebrates. In zebrafish, we found that irf6 regulates the expression of esrp1 We detailed overlapping Irf6 and Esrp1/2 expression in mouse orofacial epithelium. In zebrafish, irf6 and esrp1/2 share expression in periderm, frontonasal ectoderm and oral epithelium. Genetic disruption of irf6 and esrp1/2 in zebrafish resulted in cleft of the anterior neurocranium. The esrp1/2 mutant also developed cleft of the mouth opening. Lineage tracing of cranial neural crest cells revealed that the cleft resulted not from migration defect, but from impaired chondrogenesis. Analysis of aberrant cells within the cleft revealed expression of sox10, col1a1 and irf6, and these cells were adjacent to krt4+ and krt5+ cells. Breeding of mouse Irf6; Esrp1; Esrp2 compound mutants suggested genetic interaction, as the triple homozygote and the Irf6; Esrp1 double homozygote were not observed. Further, Irf6 heterozygosity reduced Esrp1/2 cleft severity. These studies highlight the complementary analysis of Irf6 and Esrp1/2 in mouse and zebrafish, and identify a unique aberrant cell population in zebrafish expressing sox10, col1a1 and irf6 Future work characterizing this cell population will yield additional insight into cleft pathogenesis.

Zebrafish Models of Craniofacial Malformations: Interactions of Environmental Factors

The zebrafish is an appealing model organism for investigating the genetic (G) and environmental (E) factors, as well as their interactions (GxE), which contribute to craniofacial malformations. Here, we review zebrafish studies on environmental factors involved in the etiology of craniofacial malformations in humans including maternal smoking, alcohol consumption, nutrition and drug use. As an example, we focus on the (cleft) palate, for which the zebrafish ethmoid plate is a good model. This review highlights the importance of investigating ExE interactions and discusses the variable effects of exposure to environmental factors on craniofacial development depending on dosage, exposure time and developmental stage. Zebrafish also promise to be a good tool to study novel craniofacial teratogens and toxin mixtures. Lastly, we discuss the handful of studies on gene–alcohol interactions using mutant sensitivity screens and reverse genetic techniques. We expect that studies addressing complex interactions (ExE and GxE) in craniofacial malformations will increase in the coming years. These are likely to uncover currently unknown mechanisms with implications for the prevention of craniofacial malformations. The zebrafish appears to be an excellent complementary model with high translational value to study these complex interactions.

Cellular and molecular mechanisms in the development of a cleft lip and/or cleft palate; insights from zebrafish (Danio rerio)

Human cleft lip and/or palate (CLP) are immediately recognizable congenital abnormalities of the face. Lip and palate develop from facial primordia through the coordinated activities of ectodermal epithelium and neural crest cells (NCCs) derived from ectomesenchyme tissue. Subtle changes in the regulatory mechanisms of NCC or ectodermal epithelial cells can result in CLP. Genetic and environmental contributions or a combination of both play a significant role in the progression of CLP. Model organisms provide us with a wealth of information in understanding the pathophysiology and genetic etiology of this complex disease. Small teleost, zebrafish (Danio rerio) is one of the popular model in craniofacial developmental biology. The short generation time and large number of optically transparent, easily manipulated embryos increase the value of zebrafish to identify novel candidate genes and gene regulatory networks underlying craniofacial development. In addition, it is widely used to identify the mechanisms of environmental teratogens and in therapeutic drug screening. Here, we discuss the value of zebrafish as a model to understand epithelial and NCC induced ectomesenchymal cell activities during early palate morphogenesis and robustness of the zebrafish in modern research on identifying the genetic and environmental etiological factors of CLP.

Knockout of the gene encoding the extracellular matrix protein SNED1 results in early neonatal lethality and craniofacial malformations

BACKGROUND

The extracellular matrix (ECM) is a fundamental component of multicellular organisms that orchestrates developmental processes and controls cell and tissue organization. We previously identified the novel ECM protein SNED1 as a promoter of breast cancer metastasis and showed that its level of expression negatively correlated with breast cancer patient survival. Here, we sought to identify the roles of SNED1 during murine development.

RESULTS

We generated two novel Sned1 knockout mouse strains and showed that Sned1 is essential since homozygous ablation of the gene led to early neonatal lethality. Phenotypic analysis of the surviving knockout mice revealed a role for SNED1 in the development of craniofacial and skeletal structures since Sned1 knockout resulted in growth defects, nasal cavity occlusion, and craniofacial malformations. Sned1 is widely expressed in embryos, notably by cell populations undergoing epithelial-to-mesenchymal transition, such as the neural crest cells. We further show that mice with a neural-crest-cell-specific deletion of Sned1 survive, but display facial anomalies partly phenocopying the global knockout mice.

CONCLUSIONS

Our results demonstrate requisite roles for SNED1 during development and neonatal survival. Importantly, the deletion of 2q37.3 in humans, a region that includes the SNED1 locus, has been associated with facial dysmorphism and short stature.

Modelling the developmental spliceosomal craniofacial disorder Burn-McKeown syndrome using induced pluripotent stem cells

The craniofacial developmental disorder Burn-McKeown Syndrome (BMKS) is caused by biallelic variants in the pre-messenger RNA splicing factor gene TXNL4A/DIB1. The majority of affected individuals with BMKS have a 34 base pair deletion in the promoter region of one allele of TXNL4A combined with a loss-of-function variant on the other allele, resulting in reduced TXNL4A expression. However, it is unclear how reduced expression of this ubiquitously expressed spliceosome protein results in craniofacial defects during development. Here we reprogrammed peripheral mononuclear blood cells from a BMKS patient and her unaffected mother into induced pluripotent stem cells (iPSCs) and differentiated the iPSCs into induced neural crest cells (iNCCs), the key cell type required for correct craniofacial development. BMKS patient-derived iPSCs proliferated more slowly than both mother- and unrelated control-derived iPSCs, and RNA-Seq analysis revealed significant differences in gene expression and alternative splicing. Patient iPSCs displayed defective differentiation into iNCCs compared to maternal and unrelated control iPSCs, in particular a delay in undergoing an epithelial-to-mesenchymal transition (EMT). RNA-Seq analysis of differentiated iNCCs revealed widespread gene expression changes and mis-splicing in genes relevant to craniofacial and embryonic development that highlight a dampened response to WNT signalling, the key pathway activated during iNCC differentiation. Furthermore, we identified the mis-splicing of TCF7L2 exon 4, a key gene in the WNT pathway, as a potential cause of the downregulated WNT response in patient cells. Additionally, mis-spliced genes shared common sequence properties such as length, branch point to 3’ splice site (BPS-3’SS) distance and splice site strengths, suggesting that splicing of particular subsets of genes is particularly sensitive to changes in TXNL4A expression. Together, these data provide the first insight into how reduced TXNL4A expression in BMKS patients might compromise splicing and NCC function, resulting in defective craniofacial development in the embryo.

Effects of Multi-Generational Soft Diet Consumption on Mouse Craniofacial Morphology

Variations in craniofacial morphology may arise as a result of adaptation to different environmental factors such as soft diet (SD), which lessens functional masticatory load. Prior studies have shown that changes in the masticatory muscle function associated with a switch to short-term SD led to changes in craniofacial morphology and alveolar bone architecture. However, the long-term effects of SD and the associated adaptive changes in craniofacial shape are unclear. Our novel study set out to profile prospective skull changes in mice fed with SDs over multiple generations using three-dimensional (3D) geometric morphometric analysis (GMA). Our results revealed that short-term SD consumption led to a significant decrease in craniofacial size, along with numerous shape changes. Long-term SD consumption over 15 continuous generations was not associated with changes in craniofacial size; however, shape analysis revealed mice with shortened crania and mandibles in the anteroposterior dimension, as well as relative widening in the transverse dimension compared to the average shape of all mice analyzed in our study. Moreover, changes in shape and size associated with different functional loads appeared to be independent - shape changes persisted after diets were switched for one generation, whereas size decreased after one generation and then returned to baseline size. Our study is the first to study the role of prolonged, multi-generational SD consumption in the determination of craniofacial size and shape.

What Do Animal Models Teach Us About Congenital Craniofacial Defects?

The formation of the head and face is a complex process which involves many different signaling cues regulating the migration, differentiation, and proliferation of the neural crest. This highly complex process is very error-prone, resulting in craniofacial defects in nearly 10,000 births in the United States annually. Due to the highly conserved mechanisms of craniofacial development, animal models are widely used to understand the pathogenesis of various human diseases and assist in the diagnosis and generation of preventative therapies and treatments. Here, we provide a brief background of craniofacial development and discuss several rare diseases affecting craniofacial bone development. We focus on rare congenital diseases of the cranial bone, facial jaw bones, and two classes of diseases, ciliopathies and RASopathies. Studying the animal models of these rare diseases sheds light not only on the etiology and pathology of each disease, but also provides meaningful insights towards the mechanisms which regulate normal development of the head and face.

Extracellular Matrix Composition and Remodeling: Current Perspectives on Secondary Palate Formation, Cleft Lip/Palate, and Palatal Reconstruction

Craniofacial development comprises a complex process in humans in which failures or disturbances frequently lead to congenital anomalies. Cleft lip with/without palate (CL/P) is a common congenital anomaly that occurs due to variations in craniofacial development genes, and may occur as part of a syndrome, or more commonly in isolated forms (non-syndromic). The etiology of CL/P is multifactorial with genes, environmental factors, and their potential interactions contributing to the condition. Rehabilitation of CL/P patients requires a multidisciplinary team to perform the multiple surgical, dental, and psychological interventions required throughout the patient's life. Despite progress, lip/palatal reconstruction is still a major treatment challenge. Genetic mutations and polymorphisms in several genes, including extracellular matrix (ECM) genes, soluble factors, and enzymes responsible for ECM remodeling (e.g., metalloproteinases), have been suggested to play a role in the etiology of CL/P; hence, these may be considered likely targets for the development of new preventive and/or therapeutic strategies. In this context, investigations are being conducted on new therapeutic approaches based on tissue bioengineering, associating stem cells with biomaterials, signaling molecules, and innovative technologies. In this review, we discuss the role of genes involved in ECM composition and remodeling during secondary palate formation and pathogenesis and genetic etiology of CL/P. We also discuss potential therapeutic approaches using bioactive molecules and principles of tissue bioengineering for state-of-the-art CL/P repair and palatal reconstruction.

Animal models of gene-alcohol interactions

Most birth defects arise from complex interactions between multiple genetic and environmental factors. However, our current understanding of how these interactions and their contributions affect birth defects remains incomplete. Human studies are limited in their ability to identify the fundamental causes of birth defects due to ethical and practical limitations. Animal models provide a great number of resources not available to human studies and they have been critical in advancing our understanding of birth defects and the complex interactions that underlie them. In this review, we discuss the use of animal models in the context of gene-environment interactions that underlie birth defects. We focus on alcohol which is the most common environmental factor associated with birth defects. Prenatal alcohol exposure leads to a wide range of cognitive impairments and structural deficits broadly termed fetal alcohol spectrum disorders (FASD). We discuss the broad impact of prenatal alcohol exposure on the developing embryo and elaborate on the current state of gene-alcohol interactions. Additionally, we discuss how animal models have informed our understanding of the genetics of FASD. Ultimately, these topics will provide insight into the use of animal models in understanding gene-environment interactions and their subsequent impact on birth defects.

Hunting for genes that shape human faces: Initial successes and challenges for the future

There is ample evidence from heritability studies, genetic syndromes and experimental animal models that facial morphology is strongly influenced by genes. In this brief review, we present an up-to-date overview of the efforts to identify genes associated with the size and shape of human facial features. We discuss recent methodological advances that have led to breakthroughs, but also the multitude of challenges facing the field. We offer perspective on possible applications of this line of research, particularly in the context of the precision genomics movement.

Neural crest development: insights from the zebrafish

Our understanding of the neural crest, a key vertebrate innovation, is built upon studies of multiple model organisms. Early research on neural crest cells (NCCs) was dominated by analyses of accessible amphibian and avian embryos, with mouse genetics providing complementary insights in more recent years. The zebrafish model is a relative newcomer to the field, yet it offers unparalleled advantages for the study of NCCs. Specifically, zebrafish provide powerful genetic and transgenic tools, coupled with rapidly developing transparent embryos that are ideal for high-resolution real-time imaging of the dynamic process of neural crest development. While the broad principles of neural crest development are largely conserved across vertebrate species, there are critical differences in anatomy, morphogenesis, and genetics that must be considered before information from one model is extrapolated to another. Here, our goal is to provide the reader with a helpful primer specific to neural crest development in the zebrafish model. We focus largely on the earliest events-specification, delamination, and migration-discussing what is known about zebrafish NCC development and how it differs from NCC development in non-teleost species, as well as highlighting current gaps in knowledge.

Initiation and early growth of the skull vault in zebrafish

The zebrafish offers powerful advantages as a model system for examining the growth of the skull vault and the formation of cranial sutures. The zebrafish is well suited for large-scale genetic screens, available in large numbers, and continual advances in genetic engineering facilitate precise modeling of human genetic disorders. Most importantly, zebrafish are continuously accessible for imaging during critical periods of skull formation when both mouse and chick are physically inaccessible. To establish a foundation of information on the dynamics of skull formation, we performed a longitudinal study based on confocal microscopy of individual live transgenic zebrafish. Discrete events occur at stereotyped stages in overall growth, with little variation in timing among individuals. The frontal and parietal bones initiate as small clusters of cells closely associated with cartilage around the perimeter of the skull, prior to metamorphosis and the transition to juvenile fish. Over a period of ~30 days, the frontal and parietal bones grow towards the apex of the skull and meet to begin suture formation. To aid in visualization, we have generated interactive three-dimensional models based on the imaging data, with annotated cartilage and bone elements. We propose a framework to conceptualize development of bones of the skull vault in three phases: initiation in close association with cartilage; rapid planar growth towards the apex of the skull; and finally overlapping to form sutures. Our data provide an important framework for comparing the stages and timing of skull development across model organisms, and also a baseline for the examination of zebrafish mutants affecting skull development. To facilitate these comparative analyses, the raw imaging data and the models are available as an online atlas through the FaceBase consortium (facebase.org).

Vitamin D Receptor Signaling Regulates Craniofacial Cartilage Development in Zebrafish

Vitamin D plays essential roles in supporting the skeletal system. The active form of vitamin D functions through the vitamin D receptor (VDR). A hereditary vitamin-D-resistant rickets with facial dysmorphism has been reported, but the involvement of VDR signaling during early stages of craniofacial development remains to be elucidated. The present study investigated whether VDR signaling is implicated in zebrafish craniofacial cartilage development using a morpholino-based knockdown approach. Two paralogous VDR genes, vdra and vdrb, have been found in zebrafish embryos. Loss-of-vdra has no discernible effect on cartilage elements, whereas loss-of-vdrb causes reduction and malformation of craniofacial cartilages. Disrupting both vdra and vdrb leads to more severe defects or complete loss of cartilage. Notably, knockdown of vdrb results in elevated expression of follistatin a (fsta), a bone morphogenetic protein (BMP) antagonist, in the adjacent pharyngeal endoderm. Taken together, these findings strongly indicate that VDR signaling is required for early craniofacial cartilage development in zebrafish.

Craniofacial, orofacial and dental disorders: the role of the RAS/ERK pathway

Deviations from the precisely coordinated programme of human head development can lead to craniofacial and orofacial malformations often including a variety of dental abnormalities too. Although the aetiology is still unknown in many cases, during the last decades different intracellular signalling pathways have been genetically linked to specific disorders. Among these pathways, the RAS/extracellular signal-regulated kinase (ERK) signalling cascade is the focus of this review since it encompasses a large group of genes that when mutated cause some of the most common and severe developmental anomalies in humans. We present the components of the RAS/ERK pathway implicated in craniofacial and orodental disorders through a series of human and animal studies. We attempt to unravel the specific molecular targets downstream of ERK that act on particular cell types and regulate key steps in the associated developmental processes. Finally we point to ambiguities in our current knowledge that need to be clarified before RAS/ERK-targeting therapeutic approaches can be implemented.

Facial Genetics: A Brief Overview

Historically, craniofacial genetic research has understandably focused on identifying the causes of craniofacial anomalies and it has only been within the last 10 years, that there has been a drive to detail the biological basis of normal-range facial variation. This initiative has been facilitated by the availability of low-cost hi-resolution three-dimensional systems which have the ability to capture the facial details of thousands of individuals quickly and accurately. Simultaneous advances in genotyping technology have enabled the exploration of genetic influences on facial phenotypes, both in the present day and across human history. There are several important reasons for exploring the genetics of normal-range variation in facial morphology.     - Disentangling the environmental factors and relative parental biological contributions to heritable traits can help to answer the age-old question "why we look the way that we do?"     - Understanding the etiology of craniofacial anomalies; e.g., unaffected family members of individuals with non-syndromic cleft lip/palate (nsCL/P) have been shown to differ in terms of normal-range facial variation to the general population suggesting an etiological link between facial morphology and nsCL/P.     - Many factors such as ancestry, sex, eye/hair color as well as distinctive facial features (such as, shape of the chin, cheeks, eyes, forehead, lips, and nose) can be identified or estimated using an individual's genetic data, with potential applications in healthcare and forensics.     - Improved understanding of historical selection and adaptation relating to facial phenotypes, for example, skin pigmentation and geographical latitude.     - Highlighting what is known about shared facial traits, medical conditions and genes.

Translational Genetic Modelling of 3D Craniofacial Dysmorphology: Elaborating the Facial Phenotype of Neurodevelopmental Disorders Through the "Prism" of Schizophrenia

Purpose of Review

In the context of human developmental conditions, we review the conceptualisation of schizophrenia as a neurodevelopmental disorder, the status of craniofacial dysmorphology as a clinically accessible index of brain dysmorphogenesis, the ability of genetically modified mouse models of craniofacial dysmorphology to inform on the underlying dysmorphogenic process and how geometric morphometric techniques in mutant mice can extend quantitative analysis.

Recent Findings

Mutant mice with disruption of neuregulin-1, a gene associated meta-analytically with risk for schizophrenia, constitute proof-of-concept studies of murine facial dysmorphology in a manner analogous to clinical studies in schizophrenia. Geometric morphometric techniques informed on the topography of facial dysmorphology and identified asymmetry therein.

Summary

Targeted disruption in mice of genes involved in individual components of developmental processes and analysis of resultant facial dysmorphology using geometric morphometrics can inform on mechanisms of dysmorphogenesis at levels of incisiveness not possible in human subjects.

Using Zebrafish to Test the Genetic Basis of Human Craniofacial Diseases

Genome-wide association studies (GWASs) opened an innovative and productive avenue to investigate the molecular basis of human craniofacial disease. However, GWASs identify candidate genes only; they do not prove that any particular one is the functional villain underlying disease or just an unlucky genomic bystander. Genetic manipulation of animal models is the best approach to reveal which genetic loci identified from human GWASs are functionally related to specific diseases. The purpose of this review is to discuss the potential of zebrafish to resolve which candidate genetic loci are mechanistic drivers of craniofacial diseases. Many anatomic, embryonic, and genetic features of craniofacial development are conserved among zebrafish and mammals, making zebrafish a good model of craniofacial diseases. Also, the ability to manipulate gene function in zebrafish was greatly expanded over the past 20 y, enabling systems such as Gateway Tol2 and CRISPR-Cas9 to test gain- and loss-of-function alleles identified from human GWASs in coding and noncoding regions of DNA. With the optimization of genetic editing methods, large numbers of candidate genes can be efficiently interrogated. Finding the functional villains that underlie diseases will permit new treatments and prevention strategies and will increase understanding of how gene pathways operate during normal development.

The Role of Noncoding Genetic Variation in Isolated Orofacial Clefts

In the past decade, medical genetic research has generated multiple discoveries, many of which were obtained via genome-wide association studies (GWASs). A major GWAS finding is that the majority of risk variants for complex traits map to noncoding regions. This has resulted in a paradigm shift in terms of the interpretation of human genomic sequence variation, with more attention now being paid to what was previously termed “junk DNA.” Translation of genetic findings into biologically meaningful results requires 1) large-scale and cell-specific efforts to annotate non-protein–coding regions and 2) the integration of comprehensive genomic data sets. However, this represents an enormous challenge, particularly in the case of human traits that arise during embryonic development, such as orofacial clefts (OFCs). OFC is a multifactorial trait and ranks among the most common of all human congenital malformations. These 2 attributes apply in particular to its isolated forms (nonsyndromic OFC [nsOFC]). Although genetic studies (including GWASs) have yielded novel insights into the genetic architecture of nsOFC, few data are available concerning causality and affected biological pathways. Reasons for this deficiency include the complex genetic architecture at risk loci and the limited availability of functional data sets from human tissues that represent relevant embryonic sites and time points. The present review summarizes current knowledge of the role of noncoding regions in nsOFC etiology. We describe the identification of genetic risk factors for nsOFC and several of the approaches used to identify causal variants at these loci. These strategies include the use of biological and genetic information from public databases, the assessment of the full spectrum of genetic variability within 1 locus, and comprehensive in vitro and in vivo experiments. This review also highlights the role of the emerging research field “functional genomics” and its increasing contribution to our biological understanding of nsOFC.

Systems biology of facial development: contributions of ectoderm and mesenchyme

The rapid increase in gene-centric biological knowledge coupled with analytic approaches for genomewide data integration provides an opportunity to develop systems-level understanding of facial development. Experimental analyses have demonstrated the importance of signaling between the surface ectoderm and the underlying mesenchyme are coordinating facial patterning. However, current transcriptome data from the developing vertebrate face is dominated by the mesenchymal component, and the contributions of the ectoderm are not easily identified. We have generated transcriptome datasets from critical periods of mouse face formation that enable gene expression to be analyzed with respect to time, prominence, and tissue layer. Notably, by separating the ectoderm and mesenchyme we considerably improved the sensitivity compared to data obtained from whole prominences, with more genes detected over a wider dynamic range. From these data we generated a detailed description of ectoderm-specific developmental programs, including pan-ectodermal programs, prominence- specific programs and their temporal dynamics. The genes and pathways represented in these programs provide mechanistic insights into several aspects of ectodermal development. We also used these data to identify co-expression modules specific to facial development. We then used 14 co-expression modules enriched for genes involved in orofacial clefts to make specific mechanistic predictions about genes involved in tongue specification, in nasal process patterning and in jaw development. Our multidimensional gene expression dataset is a unique resource for systems analysis of the developing face; our co-expression modules are a resource for predicting functions of poorly annotated genes, or for predicting roles for genes that have yet to be studied in the context of facial development; and our analytic approaches provide a paradigm for analysis of other complex developmental programs.

Mouth development

A mouth is present in all animals, and comprises an opening from the outside into the oral cavity and the beginnings of the digestive tract to allow eating. This review focuses on the earliest steps in mouth formation. In the first half, we conclude that the mouth arose once during evolution. In all animals, the mouth forms from ectoderm and endoderm. A direct association of oral ectoderm and digestive endoderm is present even in triploblastic animals, and in chordates, this region is known as the extreme anterior domain (EAD). Further support for a single origin of the mouth is a conserved set of genes that form a 'mouth gene program' including foxA and otx2. In the second half of this review, we discuss steps involved in vertebrate mouth formation, using the frog Xenopus as a model. The vertebrate mouth derives from oral ectoderm from the anterior neural ridge, pharyngeal endoderm and cranial neural crest (NC). Vertebrates form a mouth by breaking through the body covering in a precise sequence including specification of EAD ectoderm and endoderm as well as NC, formation of a 'pre-mouth array,' basement membrane dissolution, stomodeum formation, and buccopharyngeal membrane perforation. In Xenopus, the EAD is also a craniofacial organizer that guides NC, while reciprocally, the NC signals to the EAD to elicit its morphogenesis into a pre-mouth array. Human mouth anomalies are prevalent and are affected by genetic and environmental factors, with understanding guided in part by use of animal models. WIREs Dev Biol 2017, 6:e275. doi: 10.1002/wdev.275 For further resources related to this article, please visit the WIREs website.

Leveraging Online Resources to Prioritize Candidate Genes for Functional Analyses: Using the Fetal Testis as a Test Case

With each new microarray or RNA-seq experiment, massive quantities of transcriptomic information are generated with the purpose to produce a list of candidate genes for functional analyses. Yet an effective strategy remains elusive to prioritize the genes on these candidate lists. In this review, we outline a prioritizing strategy by taking a step back from the bench and leveraging the rich range of public databases. This in silico approach provides an economical, less biased, and more effective solution. We discuss the publicly available online resources that can be used to answer a range of questions about a gene. Is the gene of interest expressed in the system of interest (using expression databases)? Where else is this gene expressed (using added-value transcriptomic resources)? What pathways and processes is the gene involved in (using enriched gene pathway analysis and mouse knockout databases)? Is this gene correlated with human diseases (using human disease variant databases)? Using mouse fetal testis as an example, our strategies identified 298 genes annotated as expressed in the fetal testis. We cross-referenced these genes to existing microarray data and narrowed the list down to cell-type-specific candidates (35 for Sertoli cells, 11 for Leydig cells, and 25 for germ cells). Our strategies can be customized so that they allow researchers to effectively and confidently prioritize genes for functional analysis.