Regulation of facial morphogenesis by endothelin signaling: insights from mice and fish

Clinical and Genetic Correlation in Neurocristopathies: Bridging a Precision Medicine Gap

Neurocristopathies (NCPs) encompass a spectrum of disorders arising from issues during the formation and migration of neural crest cells (NCCs). NCCs undergo epithelial-mesenchymal transition (EMT) and upon key developmental gene deregulation, fetuses and neonates are prone to exhibit diverse manifestations depending on the affected area. These conditions are generally rare and often have a genetic basis, with many following Mendelian inheritance patterns, thus making them perfect candidates for precision medicine. Examples include cranial NCPs, like Goldenhar syndrome and Axenfeld-Rieger syndrome; cardiac-vagal NCPs, such as DiGeorge syndrome; truncal NCPs, like congenital central hypoventilation syndrome and Waardenburg syndrome; and enteric NCPs, such as Hirschsprung disease. Additionally, NCCs' migratory and differentiating nature makes their derivatives prone to tumors, with various cancer types categorized based on their NCC origin. Representative examples include schwannomas and pheochromocytomas. This review summarizes current knowledge of diseases arising from defects in NCCs' specification and highlights the potential of precision medicine to remedy a clinical phenotype by targeting the genotype, particularly important given that those affected are primarily infants and young children.

Identification and characterization of Hand2 upstream genomic enhancers active in developing stomach and limbs

BACKGROUND

The bHLH transcription factor HAND2 plays important roles in the development of the embryonic heart, face, limbs, and sympathetic and enteric nervous systems. To define how and when HAND2 regulates these developmental systems, requires understanding the transcriptional regulation of Hand2.

RESULTS

Remarkably, Hand2 is flanked by an extensive upstream gene desert containing a potentially diverse enhancer landscape. Here, we screened the regulatory interval 200 kb proximal to Hand2 for putative enhancers using evolutionary conservation and histone marks in Hand2-expressing tissues. H3K27ac signatures across embryonic tissues pointed to only two putative enhancer regions showing deep sequence conservation. Assessment of the transcriptional enhancer potential of these elements using transgenic reporter lines uncovered distinct in vivo enhancer activities in embryonic stomach and limb mesenchyme, respectively. Activity of the identified stomach enhancer was restricted to the developing antrum and showed expression within the smooth muscle and enteric neurons. Surprisingly, the activity pattern of the limb enhancer did not overlap Hand2 mRNA but consistently yielded a defined subectodermal anterior expression pattern within multiple transgenic lines.

CONCLUSIONS

Together, these results start to uncover the diverse regulatory potential inherent to the Hand2 upstream regulatory interval.

Alcian Blue Staining for Chondrocranium Development in Zebrafish

Craniofacial abnormalities are one of the most frequent birth malformations in humans, affecting around one in every thousand live births. The zebrafish (Danio rerio), a model organism that has seen increased usage in toxicological research in recent years, is ideal for assessing the effects of various chemicals on bone and cartilage structures. Chondrogenesis developed in zebrafish embryos by embryonic day 2, and supporting cartilage components are apparent at hatching (72 h post-fertilization). Individual cartilage may be observed using Alcian Blue staining as early as 2 days post-fertilization (dpf). The preferential binding of Alcian Blue causes the staining of zebrafish cartilage to acidic glycoproteins in an acidic solution (pH 2.2). In 72-120 hpf embryos, the cranial skeleton is easily visible after cartilage staining using Alcian Blue. Various cranial lengths and structures can be determined by measuring specific distances and angles to optimize the quantitative analysis of cranial malformations in zebrafish after exposure to various toxic agents. This chapter explains the Alcian Blue staining procedure to identify craniofacial cartilaginous structures in zebrafish embryos.

TFAP2 paralogs regulate midfacial development in part through a conserved ALX genetic pathway

Cranial neural crest development is governed by positional gene regulatory networks (GRNs). Fine-tuning of the GRN components underlies facial shape variation, yet how those networks in the midface are connected and activated remain poorly understood. Here, we show that concerted inactivation of Tfap2a and Tfap2b in the murine neural crest, even during the late migratory phase, results in a midfacial cleft and skeletal abnormalities. Bulk and single-cell RNA-seq profiling reveal that loss of both TFAP2 family members dysregulates numerous midface GRN components involved in midface morphogenesis, patterning and differentiation. Notably, Alx1, Alx3 and Alx4 (ALX) transcript levels are reduced, whereas ChIP-seq analyses suggest TFAP2 family members directly and positively regulate ALX gene expression. Tfap2a, Tfap2b and ALX co-expression in midfacial neural crest cells of both mouse and zebrafish implies conservation of this regulatory axis across vertebrates. Consistent with this notion, tfap2a zebrafish mutants present with abnormal alx3 expression patterns, Tfap2a binds ALX loci and tfap2a-alx3 genetic interactions are observed. Together, these data demonstrate TFAP2 paralogs regulate vertebrate midfacial development in part by activating expression of ALX transcription factor genes.

Endothelin signaling in development

Since the discovery of endothelin 1 (EDN1) in 1988, the role of endothelin ligands and their receptors in the regulation of blood pressure in normal and disease states has been extensively studied. However, endothelin signaling also plays crucial roles in the development of neural crest cell-derived tissues. Mechanisms of endothelin action during neural crest cell maturation have been deciphered using a variety of in vivo and in vitro approaches, with these studies elucidating the basis of human syndromes involving developmental differences resulting from altered endothelin signaling. In this Review, we describe the endothelin pathway and its functions during the development of neural crest-derived tissues. We also summarize how dysregulated endothelin signaling causes developmental differences and how this knowledge may lead to potential treatments for individuals with gene variants in the endothelin pathway.

Foxi3GFP and Foxi3CreER mice allow identification and lineage labeling of pharyngeal arch ectoderm and endoderm, and tooth and hair placodes

BACKGROUND

FOXI3 is a forkhead family transcription factor that is expressed in the progenitors of craniofacial placodes, epidermal placodes, and the ectoderm and endoderm of the pharyngeal arch region. Loss of Foxi3 in mice and pathogenic Foxi3 variants in dogs and humans cause a variety of craniofacial defects including absence of the inner ear, severe truncations of the jaw, loss or reduction in external and middle ear structures, and defects in teeth and hair.

RESULTS

To allow for the identification, isolation, and lineage tracing of Foxi3-expressing cells in developing mice, we targeted the Foxi3 locus to create Foxi3GFP and Foxi3CreER mice. We show that Foxi3GFP mice faithfully recapitulate the expression pattern of Foxi3 mRNA at all ages examined, and Foxi3CreER mice can trace the derivatives of pharyngeal arch ectoderm and endoderm, the pharyngeal pouches and clefts that separate each arch, and the derivatives of hair and tooth placodes.

CONCLUSIONS

Foxi3GFP and Foxi3CreER mice are new tools that will be of use in identifying and manipulating pharyngeal arch ectoderm and endoderm and hair and tooth placodes.

Identification of an adverse outcome pathway (AOP) for chemical-induced craniofacial anomalies using the transgenic zebrafish model

Craniofacial anomalies are one of the most frequent birth defects worldwide and are often caused by genetic and environmental factors such as pharmaceuticals and chemical agents. Although identifying adverse outcome pathways (AOPs) is a central issue for evaluating the teratogenicity, the AOP causing craniofacial anomalies has not been identified. Recently, zebrafish has gained interest as an emerging model for predicting teratogenicity because of high throughput, cost-effectiveness and availability of various tools for examining teratogenic mechanisms. Here, we established zebrafish sox10-EGFP reporter lines to visualize cranial neural crest cells (CNCCs) and have identified the AOPs for craniofacial anomalies. When we exposed the transgenic embryos to teratogens that were reported to cause craniofacial anomalies in mammals, CNCC migration and subsequent morphogenesis of the first pharyngeal arch were impaired at 24 hours post-fertilization. We also found that cell proliferation and apoptosis of the migratory CNCCs were disturbed, which would be key events of the AOP. From these results, we propose that our sox10-EGFP reporter lines serve as a valuable model for detecting craniofacial skeletal abnormalities, from early to late developmental stages. Given that the developmental process of CNCCs around this stage is highly conserved between zebrafish and mammals, our findings can be extrapolated to mammalian craniofacial development and thus help in predicting craniofacial anomalies in human.

Transcriptomic analysis reveals the role of SIX1 in mouse cranial neural crest patterning and bone development

BACKGROUND

Genetic variants of the transcription factor SIX1 and its co-factor EYA1 underlie 50% of Branchio-oto-renal syndrome (BOR) cases. BOR is characterized by craniofacial defects, including malformed middle ear ossicles leading to conductive hearing loss. In this work, we expand our knowledge of the Six1 gene regulatory network by using a Six1-null mouse line to assess gene expression profiles of E10.5 mandibular arches, which give rise to the neural crest (NC)-derived middle ear ossicles and lower jaw, via bulk RNA sequencing.

RESULTS

Our transcriptomic analysis led to the identification of 808 differentially expressed genes that are related to translation, NC cell differentiation, osteogenesis, and chondrogenesis including components of the WNT signaling pathway. As WNT signaling is a known contributor to bone development, we demonstrated that SIX1 is required for expression of the WNT antagonist Frzb in the mandibular arch, and determined that SIX1 expression results in repression of WNT signaling.

CONCLUSION

Our results clarify the mechanisms by which SIX1 regulates the development of NC-derived craniofacial elements that are altered in SIX1-associated disorders. In addition, this work identifies novel genes that could be causative to this birth defect and establishes a link between SIX1 and WNT signaling during patterning of NC cells.

Linking Vertebrate Gene Duplications to the New Head Hypothesis

Vertebrates have diverse morphologies and various anatomical novelties that set them apart from their closest invertebrate relatives. A conspicuous head housing a large brain, paired sense organs, and protected by a skeleton of cartilage and bone is unique to vertebrates and is a defining feature of this taxon. Gans and Northcutt (1980s) proposed that the evolution of this "new head" was dependent on two key developmental innovations: neural crest cells (NCCs) and ectodermal placodes. NCCs are migratory embryonic cells that form bone, cartilage, and neurons in the new head. Based on genome size, Ohno (1970s) proposed a separate hypothesis, stating that vertebrate genome content was quadrupled via two rounds (2R) of whole genome duplications (WGDs), and the surplus of genetic material potentiated vertebrate morphological diversification. While both hypotheses offer explanations for vertebrate success, it is unclear if, and how, the "new head" and "2R" hypotheses are linked. Here, we consider both hypotheses and evaluate the experimental evidence connecting the two. Overall, evidence suggests that while the origin of the NC GRN predates the vertebrate WGDs, these genomic events may have potentiated the evolution of distinct genetic subnetworks in different neural crest subpopulations. We describe the general composition of the NC GRN and posit that its increased developmental modularity facilitated the independent evolution of NC derivatives and the diversification of the vertebrate head skeleton. Lastly, we discuss experimental strategies needed to test whether gene duplications drove the diversification of neural crest derivatives and the "new head".

Auriculocondylar syndrome: Pathogenesis, clinical manifestations and surgical therapies

Auriculocondylar syndrome (ARCND) is a genetic and rare craniofacial condition caused by abnormal development of the first and second pharyngeal arches during the embryonic stage and is characterized by peculiar auricular malformations (question mark ears), mandibular condyle hypoplasia, micrognathia and other less-frequent features. GNAI3, PLCB4 and EDN1 have been identified as pathogenic genes in this syndrome so far, all of which are implicated in the EDN1-EDNRA signal pathway. Therefore, ARCND is genetically classified as ARCND1, ARCND2 and ARCND3 based on the mutations in GNAI3, PLCB4 and EDN1, respectively. ARCND is inherited in an autosomal dominant or recessive mode with significant intra- and interfamilial phenotypic variation and incomplete penetrance, rendering its diagnosis difficult and therapies individualized. To raise clinicians' awareness of the rare syndrome, we focused on the currently known pathogenesis, pathogenic genes, clinical manifestations and surgical therapies in this review.

Stripes and loss of color in ball pythons (Python regius) are associated with variants affecting endothelin signaling

Color patterns in nonavian reptiles are beautifully diverse, but little is known about the genetics and development of these patterns. Here, we investigated color patterning in pet ball pythons (Python regius), which have been bred to show color phenotypes that differ dramatically from the wildtype form. We report that several color phenotypes in pet animals are associated with putative loss-of-function variants in the gene encoding endothelin receptor EDNRB1: (1) frameshift variants in EDNRB1 are associated with conversion of the normal mottled color pattern to skin that is almost fully white, (2) missense variants affecting conserved sites of the EDNRB1 protein are associated with dorsal, longitudinal stripes, and (3) substitutions at EDNRB1 splice donors are associated with subtle changes in patterning compared to wildtype. We propose that these phenotypes are caused by loss of specialized color cells (chromatophores), with loss ranging from severe (fully white) to moderate (dorsal striping) to mild (subtle changes in patterning). Our study is the first to describe variants affecting endothelin signaling in a nonavian reptile and suggests that reductions in endothelin signaling in ball pythons can produce a variety of color phenotypes, depending on the degree of color cell loss.

TFAP2 paralogs regulate midfacial development in part through a conserved ALX genetic pathway

Cranial neural crest development is governed by positional gene regulatory networks (GRNs). Fine-tuning of the GRN components underly facial shape variation, yet how those in the midface are connected and activated remain poorly understood. Here, we show that concerted inactivation of Tfap2a and Tfap2b in the murine neural crest even during the late migratory phase results in a midfacial cleft and skeletal abnormalities. Bulk and single-cell RNA-seq profiling reveal that loss of both Tfap2 members dysregulated numerous midface GRN components involved in midface fusion, patterning, and differentiation. Notably, Alx1/3/4 (Alx) transcript levels are reduced, while ChIP-seq analyses suggest TFAP2 directly and positively regulates Alx gene expression. TFAP2 and ALX co-expression in midfacial neural crest cells of both mouse and zebrafish further implies conservation of this regulatory axis across vertebrates. Consistent with this notion, tfap2a mutant zebrafish present abnormal alx3 expression patterns, and the two genes display a genetic interaction in this species. Together, these data demonstrate a critical role for TFAP2 in regulating vertebrate midfacial development in part through ALX transcription factor gene expression.

Neural crest mechanosensors: Seeing old proteins in a new light

Mechanical forces exerted on neural crest cells control their collective migration and differentiation. This perspective discusses our current understanding of neural crest mechanotransduction during cell migration and differentiation. Additionally, we describe proteins that have mechanosensitive functions in other systems, such as mechanosensitive G-protein-coupled receptors, mechanosensitive ion channels, cell-cell adhesion, and cell-matrix-interacting proteins, and highlight that these same proteins have in the past been studied in neural crest development from a purely signaling point of view. We propose that future studies elucidate the mechanosensitive functions these receptors may play in neural crest development and integrate this with their known molecular role.

Candidate positive targets of LHX6 and LHX8 transcription factors in the developing upper jaw

Craniofacial development is controlled by a large number of genes, which interact with one another to form a complex gene regulatory network (GRN). Key components of GRN are signaling molecules and transcription factors. Therefore, identifying targets of core transcription factors is an important part of the overall efforts toward building a comprehensive and accurate model of GRN. LHX6 and LHX8 are transcription factors expressed in the oral mesenchyme of the first pharyngeal arch (PA1), and they are crucial regulators of palate and tooth development. Previously, we performed genome-wide transcriptional profiling and chromatin immunoprecipitation to identify target genes of LHX6 and LHX8 in PA1, and described a set of genes repressed by LHX. However, there has not been any discussion of the genes positively regulated by LHX6 and LHX8. In this paper, we revisited the above datasets to identify candidate positive targets of LHX in PA1. Focusing on those with known connections to craniofacial development, we performed RNA in situ hybridization to confirm the changes in expression in Lhx6;Lhx8 mutant. We also confirmed the binding of LHX6 to several putative enhancers near the candidate target genes. Together, we have uncovered novel connections between Lhx and other important regulators of craniofacial development, including Eya1, Barx1, Rspo2, Rspo3, and Wnt11.

Further delineation of Auriculocondylar syndrome based on 14 novel cases and reassessment of 25 published cases

Auriculocondylar syndrome (ACS) is a rare craniofacial disorder characterized by mandibular hypoplasia and an auricular defect at the junction between the lobe and helix, known as a "Question Mark Ear" (QME). Several additional features, originating from the first and second branchial arches and other tissues, have also been reported. ACS is genetically heterogeneous with autosomal dominant and recessive modes of inheritance. The mutations identified to date are presumed to dysregulate the endothelin 1 signalling pathway. Here we describe 14 novel cases and reassess 25 published cases of ACS through a questionnaire for systematic data collection. All patients harbour mutation(s) in PLCB4, GNAI3 or EDN1. This series of patients contributes to the characterization of additional features occasionally associated with ACS such as respiratory, costal, neurodevelopmental and genital anomalies, and provides management and monitoring recommendations. This article is protected by copyright. All rights reserved.

The crucial role of model systems in understanding the complexity of cell signaling in human neurocristopathies

Animal models are useful to study the molecular, cellular, and morphogenetic mechanisms underlying normal and pathological development. Cell-based study models have emerged as an alternative approach to study many aspects of human embryonic development and disease. The neural crest (NC) is a transient, multipotent, and migratory embryonic cell population that generates a diverse group of cell types that arises during vertebrate development. The abnormal formation or development of the NC results in neurocristopathies (NCPs), which are characterized by a broad spectrum of functional and morphological alterations. The impaired molecular mechanisms that give rise to these multiphenotypic diseases are not entirely clear yet. This fact, added to the high incidence of these disorders in the newborn population, has led to the development of systematic approaches for their understanding. In this article, we have systematically reviewed the ways in which experimentation with different animal and cell model systems has improved our knowledge of NCPs, and how these advances might contribute to the development of better diagnostic and therapeutic tools for the treatment of these pathologies. This article is categorized under: Congenital Diseases > Genetics/Genomics/Epigenetics Congenital Diseases > Stem Cells and Development Congenital Diseases > Molecular and Cellular Physiology Neurological Diseases > Genetics/Genomics/Epigenetics.

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.

GATA3 is essential for separating patterning domains during facial morphogenesis

Neural crest cells (NCCs) within the mandibular and maxillary prominences of the first pharyngeal arch are initially competent to respond to signals from either region. However, mechanisms that are only partially understood establish developmental tissue boundaries to ensure spatially correct patterning. In the 'hinge and caps' model of facial development, signals from both ventral prominences (the caps) pattern the adjacent tissues whereas the intervening region, referred to as the maxillomandibular junction (the hinge), maintains separation of the mandibular and maxillary domains. One cap signal is GATA3, a member of the GATA family of zinc-finger transcription factors with a distinct expression pattern in the ventral-most part of the mandibular and maxillary portions of the first arch. Here, we show that disruption of Gata3 in mouse embryos leads to craniofacial microsomia and syngnathia (bony fusion of the upper and lower jaws) that results from changes in BMP4 and FGF8 gene regulatory networks within NCCs near the maxillomandibular junction. GATA3 is thus a crucial component in establishing the network of factors that functionally separate the upper and lower jaws during development.

Genome sequences of Tropheus moorii and Petrochromis trewavasae, two eco-morphologically divergent cichlid fishes endemic to Lake Tanganyika

With more than 1000 species, East African cichlid fishes represent the fastest and most species-rich vertebrate radiation known, providing an ideal model to tackle molecular mechanisms underlying recurrent adaptive diversification. We add high-quality genome reconstructions for two phylogenetic key species of a lineage that diverged about ~ 3-9 million years ago (mya), representing the earliest split of the so-called modern haplochromines that seeded additional radiations such as those in Lake Malawi and Victoria. Along with the annotated genomes we analysed discriminating genomic features of the study species, each representing an extreme trophic morphology, one being an algae browser and the other an algae grazer. The genomes of Tropheus moorii (TM) and Petrochromis trewavasae (PT) comprise 911 and 918 Mbp with 40,300 and 39,600 predicted genes, respectively. Our DNA sequence data are based on 5 and 6 individuals of TM and PT, and the transcriptomic sequences of one individual per species and sex, respectively. Concerning variation, on average we observed 1 variant per 220 bp (interspecific), and 1 variant per 2540 bp (PT vs PT)/1561 bp (TM vs TM) (intraspecific). GO enrichment analysis of gene regions affected by variants revealed several candidates which may influence phenotype modifications related to facial and jaw morphology, such as genes belonging to the Hedgehog pathway (SHH, SMO, WNT9A) and the BMP and GLI families.

The power of zebrafish models for understanding the co-occurrence of craniofacial and limb disorders

Craniofacial and limb defects are two of the most common congenital anomalies in the general population. Interestingly, these defects are not mutually exclusive. Many patients with craniofacial phenotypes, such as orofacial clefting and craniosynostosis, also present with limb defects, including polydactyly, syndactyly, brachydactyly, or ectrodactyly. The gene regulatory networks governing craniofacial and limb development initially seem distinct from one another, and yet these birth defects frequently occur together. Both developmental processes are highly conserved among vertebrates, and zebrafish have emerged as an advantageous model due to their high fecundity, relative ease of genetic manipulation, and transparency during development. Here we summarize studies that have used zebrafish models to study human syndromes that present with both craniofacial and limb phenotypes. We discuss the highly conserved processes of craniofacial and limb/fin development and describe recent zebrafish studies that have explored the function of genes associated with human syndromes with phenotypes in both structures. We attempt to identify commonalities between the two to help explain why craniofacial and limb anomalies often occur together.

Evolution of the endothelin pathway drove neural crest cell diversification

Neural crest cells (NCCs) are migratory, multipotent embryonic cells that are unique to vertebrates and form an array of clade-defining adult features. The evolution of NCCs has been linked to various genomic events, including the evolution of new gene-regulatory networks^1,2, the de novo evolution of genes³ and the proliferation of paralogous genes during genome-wide duplication events⁴. However, conclusive functional evidence linking new and/or duplicated genes to NCC evolution is lacking. Endothelin ligands (Edns) and endothelin receptors (Ednrs) are unique to vertebrates^3,5,6, and regulate multiple aspects of NCC development in jawed vertebrates^7-10. Here, to test whether the evolution of Edn signalling was a driver of NCC evolution, we used CRISPR-Cas9 mutagenesis¹¹ to disrupt edn, ednr and dlx genes in the sea lamprey, Petromyzon marinus. Lampreys are jawless fishes that last shared a common ancestor with modern jawed vertebrates around 500 million years ago¹². Thus, comparisons between lampreys and gnathostomes can identify deeply conserved and evolutionarily flexible features of vertebrate development. Using the frog Xenopus laevis to expand gnathostome phylogenetic representation and facilitate side-by-side analyses, we identify ancient and lineage-specific roles for Edn signalling. These findings suggest that Edn signalling was activated in NCCs before duplication of the vertebrate genome. Then, after one or more genome-wide duplications in the vertebrate stem, paralogous Edn pathways functionally diverged, resulting in NCC subpopulations with different Edn signalling requirements. We posit that this new developmental modularity facilitated the independent evolution of NCC derivatives in stem vertebrates. Consistent with this, differences in Edn pathway targets are associated with differences in the oropharyngeal skeleton and autonomic nervous system of lampreys and modern gnathostomes. In summary, our work provides functional genetic evidence linking the origin and duplication of new vertebrate genes with the stepwise evolution of a defining vertebrate novelty.

Neural crest cells require Meis2 for patterning the mandibular arch via the Sonic hedgehog pathway

Cranial neural crest cells (cNCCs) originate in the anterior neural tube and populate pharyngeal arches in which they contribute to formation of bone and cartilage. This cell population also provides molecular signals for the development of tissues of non-neural crest origin, such as the tongue muscles, teeth enamel or gland epithelium. Here we show that the transcription factor Meis2 is expressed in the oral region of the first pharyngeal arch (PA1) and later in the tongue primordium. Conditional inactivation of Meis2 in cNCCs resulted in loss of Sonic hedgehog signalling in the oropharyngeal epithelium and impaired patterning of PA1 along the lateral-medial and oral-aboral axis. Failure of molecular specification of PA1, illustrated by altered expression of Hand1/2, Dlx5, Barx1, Gsc and other markers, led to hypoplastic tongue and ectopic ossification of the mandible. Meis2-mutant mice thus display craniofacial defects that are reminiscent of several human syndromes and patients with mutations in the Meis2 gene.

Loss-of-function of Endothelin receptor type A results in Oro-Oto-Cardiac syndrome

Craniofacial morphogenesis is regulated in part by signaling from the Endothelin receptor type A (EDNRA). Pathogenic variants in EDNRA signaling pathway components EDNRA, GNAI3, PCLB4, and EDN1 cause Mandibulofacial Dysostosis with Alopecia (MFDA), Auriculocondylar syndrome (ARCND) 1, 2, and 3, respectively. However, cardiovascular development is normal in MFDA and ARCND individuals, unlike Ednra knockout mice. One explanation may be that partial EDNRA signaling remains in MFDA and ARCND, as mice with reduced, but not absent, EDNRA signaling also lack a cardiovascular phenotype. Here we report an individual with craniofacial and cardiovascular malformations mimicking the Ednra ^-/- mouse phenotype, including a distinctive micrognathia with microstomia and a hypoplastic aortic arch. Exome sequencing found a novel homozygous missense variant in EDNRA (c.1142A>C; p.Q381P). Bioluminescence resonance energy transfer assays revealed that this amino acid substitution in helix 8 of EDNRA prevents recruitment of G proteins to the receptor, abrogating subsequent receptor activation by its ligand, Endothelin-1. This homozygous variant is thus the first reported loss-of-function EDNRA allele, resulting in a syndrome we have named Oro-Oto-Cardiac Syndrome. Further, our results illustrate that EDNRA signaling is required for both normal human craniofacial and cardiovascular development, and that limited EDNRA signaling is likely retained in ARCND and MFDA individuals. This work illustrates a straightforward approach to identifying the functional consequence of novel genetic variants in signaling molecules associated with malformation syndromes.

Selective breeding modifies mef2ca mutant incomplete penetrance by tuning the opposing Notch pathway

Deleterious genetic mutations allow developmental biologists to understand how genes control development. However, not all loss of function genetic mutants develop phenotypic changes. Many deleterious mutations only produce a phenotype in a subset of mutant individuals, a phenomenon known as incomplete penetrance. Incomplete penetrance can confound analyses of gene function and our understanding of this widespread phenomenon remains inadequate. To better understand what controls penetrance, we capitalized on the zebrafish mef2ca mutant which produces craniofacial phenotypes with variable penetrance. Starting with a characterized mef2ca loss of function mutant allele, we used classical selective breeding methods to generate zebrafish strains in which mutant-associated phenotypes consistently appear with low or high penetrance. Strikingly, our selective breeding for low penetrance converted the mef2ca mutant allele behavior from homozygous lethal to homozygous viable. Meanwhile, selective breeding for high penetrance converted the mef2ca mutant allele from fully recessive to partially dominant. Comparing the selectively-bred low- and high-penetrance strains revealed that the strains initially respond similarly to the mutation, but then gene expression differences between strains emerge during development. Thus, altered temporal genetic circuitry can manifest through selective pressure to modify mutant penetrance. Specifically, we demonstrate differences in Notch signaling between strains, and further show that experimental manipulation of the Notch pathway phenocopies penetrance changes occurring through selective breeding. This study provides evidence that penetrance is inherited as a liability-threshold trait. Our finding that vertebrate animals can overcome a deleterious mutation by tuning genetic circuitry complements other reported mechanisms of overcoming deleterious mutations such as transcriptional adaptation of compensatory genes, alternative mRNA splicing, and maternal deposition of wild-type transcripts, which are not observed in our system. The selective breeding approach and the resultant genetic circuitry change we uncovered advances and expands our current understanding of genetic and developmental resilience.

Some deleterious gene mutations only affect a subset of genetically mutant animals. This widespread phenomenon, known as mutant incomplete penetrance, complicates discovery of causative gene mutations in both model organisms and human disease. This study utilized the zebrafish mef2ca transcription factor mutant that produces craniofacial skeleton defects with incomplete penetrance. Selectively breeding zebrafish families for low- or high-penetrance mutants for many generations created different zebrafish strains with consistently low or high penetrance. Comparing these strains allowed us to gain insight into the mechanisms that control penetrance. Specifically, genes under the control of mef2ca are initially similarly expressed between the two strains, but differences between strains emerge during development. We found that genetic manipulation of these downstream genes mimics the effects of our selective breeding. Thus, selective breeding for penetrance can change the genetic circuitry downstream of the mutated gene. We propose that small differences in gene circuitry between individuals is one mechanism underlying susceptibility or resilience to genetic mutations.

Nucleoporin 62-Like Protein is Required for the Development of Pharyngeal Arches through Regulation of Wnt/β-Catenin Signaling and Apoptotic Homeostasis in Zebrafish

We have previously observed the predominant expression of nucleoporin 62-like (Nup62l) mRNA in the pharyngeal region of zebrafish, which raises the question whether Nup62l has important implications in governing the morphogenesis of pharyngeal arches (PA) in zebrafish. Herein, we explored the functions of Nup62l in PA development. The disruption of Nup62l with a CRISPR/Cas9-dependent gene knockout approach led to defective PA, which was characterized by a thinned and shortened pharyngeal region and a significant loss of pharyngeal cartilages. During pharyngeal cartilage formation, prechondrogenic condensation and chondrogenic differentiation were disrupted in homozygous nup62l-mutants, while the specification and migration of cranial neural crest cells (CNCCs) were unaffected. Mechanistically, the impaired PA region of nup62l-mutants underwent extensive apoptosis, which was mainly dependent on activation of p53-dependent apoptotic pathway. Moreover, aberrant activation of a series of apoptotic pathways in nup62l-mutants is closely associated with the inactivation of Wnt/β-catenin signaling. Thus, these findings suggest that the regulation of Wnt/β-catenin activity by Nup62l is crucial for PA formation in zebrafish.

The role of GPCRs in bone diseases and dysfunctions

The superfamily of G protein-coupled receptors (GPCRs) contains immense structural and functional diversity and mediates a myriad of biological processes upon activation by various extracellular signals. Critical roles of GPCRs have been established in bone development, remodeling, and disease. Multiple human GPCR mutations impair bone development or metabolism, resulting in osteopathologies. Here we summarize the disease phenotypes and dysfunctions caused by GPCR gene mutations in humans as well as by deletion in animals. To date, 92 receptors (5 glutamate family, 67 rhodopsin family, 5 adhesion, 4 frizzled/taste2 family, 5 secretin family, and 6 other 7TM receptors) have been associated with bone diseases and dysfunctions (36 in humans and 72 in animals). By analyzing data from these 92 GPCRs, we found that mutation or deletion of different individual GPCRs could induce similar bone diseases or dysfunctions, and the same individual GPCR mutation or deletion could induce different bone diseases or dysfunctions in different populations or animal models. Data from human diseases or dysfunctions identified 19 genes whose mutation was associated with human BMD: 9 genes each for human height and osteoporosis; 4 genes each for human osteoarthritis (OA) and fracture risk; and 2 genes each for adolescent idiopathic scoliosis (AIS), periodontitis, osteosarcoma growth, and tooth development. Reports from gene knockout animals found 40 GPCRs whose deficiency reduced bone mass, while deficiency of 22 GPCRs increased bone mass and BMD; deficiency of 8 GPCRs reduced body length, while 5 mice had reduced femur size upon GPCR deletion. Furthermore, deficiency in 6 GPCRs induced osteoporosis; 4 induced osteoarthritis; 3 delayed fracture healing; 3 reduced arthritis severity; and reduced bone strength, increased bone strength, and increased cortical thickness were each observed in 2 GPCR-deficiency models. The ever-expanding number of GPCR mutation-associated diseases warrants accelerated molecular analysis, population studies, and investigation of phenotype correlation with SNPs to elucidate GPCR function in human diseases.

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.

Endothelin receptor Aa regulates proliferation and differentiation of Erb-dependent pigment progenitors in zebrafish

Skin pigment patterns are important, being under strong selection for multiple roles including camouflage and UV protection. Pigment cells underlying these patterns form from adult pigment stem cells (APSCs). In zebrafish, APSCs derive from embryonic neural crest cells, but sit dormant until activated to produce pigment cells during metamorphosis. The APSCs are set-aside in an ErbB signaling dependent manner, but the mechanism maintaining quiescence until metamorphosis remains unknown. Mutants for a pigment pattern gene, parade, exhibit ectopic pigment cells localised to the ventral trunk, but also supernumerary cells restricted to the Ventral Stripe. Contrary to expectations, these melanocytes and iridophores are discrete cells, but closely apposed. We show that parade encodes Endothelin receptor Aa, expressed in the blood vessels, most prominently in the medial blood vessels, consistent with the ventral trunk phenotype. We provide evidence that neuronal fates are not affected in parade mutants, arguing against transdifferentiation of sympathetic neurons to pigment cells. We show that inhibition of BMP signaling prevents specification of sympathetic neurons, indicating conservation of this molecular mechanism with chick and mouse. However, inhibition of sympathetic neuron differentiation does not enhance the parade phenotype. Instead, we pinpoint ventral trunk-restricted proliferation of neural crest cells as an early feature of the parade phenotype. Importantly, using a chemical genetic screen for rescue of the ectopic pigment cell phenotype of parade mutants (whilst leaving the embryonic pattern untouched), we identify ErbB inhibitors as a key hit. The time-window of sensitivity to these inhibitors mirrors precisely the window defined previously as crucial for the setting aside of APSCs in the embryo, strongly implicating adult pigment stem cells as the source of the ectopic pigment cells. We propose that a novel population of APSCs exists in association with medial blood vessels, and that their quiescence is dependent upon Endothelin-dependent factors expressed by the blood vessels.

BMP signaling is required for nkx2.3-positive pharyngeal pouch progenitor specification in zebrafish

Pharyngeal pouches, a series of outpocketings that bud from the foregut endoderm, are essential to the formation of craniofacial skeleton as well as several important structures like parathyroid and thymus. However, whether pharyngeal pouch progenitors exist in the developing gut tube remains unknown. Here, taking advantage of cell lineage tracing and transgenic ablation technologies, we identified a population of nkx2.3⁺ pouch progenitors in zebrafish embryos and demonstrated an essential requirement of ectodermal BMP2b for their specification. At early somite stages, nkx2.3⁺ cells located at lateral region of pharyngeal endoderm give rise to the pouch epithelium except a subpopulation expressing pdgfαa rather than nkx2.3. A small-scale screen of chemical inhibitors reveals that BMP signaling is necessary to specify these progenitors. Loss-of-function analyses show that BMP2b, expressed in the pharyngeal ectoderm, actives Smad effectors in endodermal cells to induce nkx2.3⁺ progenitors. Collectively, our study provides in vivo evidence for the existence of pouch progenitors and highlights the importance of BMP2b signaling in progenitor specification.

Pharyngeal pouches are essential to the formation of craniofacial skeleton as well as several important structures like parathyroid and thymus, but whether their progenitors exist in the developing gut tube remains unknown. Our study provide in vivo evidence that, in the early somite stages, nkx2.3⁺ cells are present in the lateral pharyngeal endoderm and give rise to the pouch epithelium. We further reveal that ectodermal BMP2b is essential for the activation of Smad effectors in endodermal cells, thereby facilitating pouch progenitor specification. Collectively, our discoveries shed new light on the cellular and molecular mechanisms of pharyngeal pouch development.

Transcriptomics reveals complex kinetics of dorsal-ventral patterning gene expression in the mandibular arch

The mandibular or first pharyngeal arch forms the upper and lower jaws in all gnathostomes. A gene regulatory network that defines ventral, intermediate, and dorsal domains along the dorsal-ventral (D-V) axis of the arch has emerged from studies in zebrafish and mice, but the temporal dynamics of this process remain unclear. To define cell fate trajectories in the arches we have performed quantitative gene expression analyses of D-V patterning genes in pharyngeal arch primordia in zebrafish and mice. Using NanoString technology to measure transcript numbers per cell directly we show that, in many cases, genes expressed in similar D-V domains and induced by similar signals vary dramatically in their temporal profiles. This suggests that cellular responses to D-V patterning signals are likely shaped by the baseline kinetics of target gene expression. Furthermore, similarities in the temporal dynamics of genes that occupy distinct pathways suggest novel shared modes of regulation. Incorporating these gene expression kinetics into our computational models for the mandibular arch improves the accuracy of patterning, and facilitates temporal comparisons between species. These data suggest that the magnitude and timing of target gene expression help diversify responses to patterning signals during craniofacial development.

Modeling craniofacial development reveals spatiotemporal constraints on robust patterning of the mandibular arch

How does pattern formation occur accurately when confronted with tissue growth and stochastic fluctuations (noise) in gene expression? Dorso-ventral (D-V) patterning of the mandibular arch specifies upper versus lower jaw skeletal elements through a combination of Bone morphogenetic protein (Bmp), Endothelin-1 (Edn1), and Notch signaling, and this system is highly robust. We combine NanoString experiments of early D-V gene expression with live imaging of arch development in zebrafish to construct a computational model of the D-V mandibular patterning network. The model recapitulates published genetic perturbations in arch development. Patterning is most sensitive to changes in Bmp signaling, and the temporal order of gene expression modulates the response of the patterning network to noise. Thus, our integrated systems biology approach reveals non-intuitive features of the complex signaling system crucial for craniofacial development, including novel insights into roles of gene expression timing and stochasticity in signaling and gene regulation.

News from the endothelin-3/EDNRB signaling pathway: Role during enteric nervous system development and involvement in neural crest-associated disorders

The endothelin system is a vertebrate-specific innovation with important roles in regulating the cardiovascular system and renal and pulmonary processes, as well as the development of the vertebrate-specific neural crest cell population and its derivatives. This system is comprised of three structurally similar 21-amino acid peptides that bind and activate two G-protein coupled receptors. In 1994, knockouts of the Edn3 and Ednrb genes revealed their crucial function during development of the enteric nervous system and melanocytes, two neural-crest derivatives. Since then, human and mouse genetics, combined with cellular and developmental studies, have helped to unravel the role of this signaling pathway during development and adulthood. In this review, we will summarize the known functions of the EDN3/EDNRB pathway during neural crest development, with a specific focus on recent scientific advances, and the enteric nervous system in normal and pathological conditions.

Neurocristopathies: New insights 150 years after the neural crest discovery

The neural crest (NC) is a transient, multipotent and migratory cell population that generates an astonishingly diverse array of cell types during vertebrate development. These cells, which originate from the ectoderm in a region lateral to the neural plate in the neural fold, give rise to neurons, glia, melanocytes, chondrocytes, smooth muscle cells, odontoblasts and neuroendocrine cells, among others. Neurocristopathies (NCP) are a class of pathologies occurring in vertebrates, especially in humans that result from the abnormal specification, migration, differentiation or death of neural crest cells during embryonic development. Various pigment, skin, thyroid and hearing disorders, craniofacial and heart abnormalities, malfunctions of the digestive tract and tumors can also be considered as neurocristopathies. In this review we revisit the current classification and propose a new way to classify NCP based on the embryonic origin of the affected tissues, on recent findings regarding the molecular mechanisms that drive NC formation, and on the increased complexity of current molecular embryology techniques.

Essential Role of Nr2f Nuclear Receptors in Patterning the Vertebrate Upper Jaw

The jaw is central to the extensive variety of feeding and predatory behaviors across vertebrates. The bones of the lower but not upper jaw form around an early-developing cartilage template. Whereas Endothelin1 patterns the lower jaw, the factors that specify upper-jaw morphology remain elusive. Here, we identify Nuclear Receptor 2f genes (Nr2fs) as enriched in and required for upper-jaw formation in zebrafish. Combinatorial loss of Nr2fs transforms maxillary components of the upper jaw into lower-jaw-like structures. Conversely, nr2f5 misexpression disrupts lower-jaw development. Genome-wide analyses reveal that Nr2fs repress mandibular gene expression and early chondrogenesis in maxillary precursors. Rescue of lower-jaw defects in endothelin1 mutants by reducing Nr2f dosage further demonstrates that Nr2f expression must be suppressed for normal lower-jaw development. We propose that Nr2fs shape the upper jaw by protecting maxillary progenitors from early chondrogenesis, thus preserving cells for later osteogenesis.

AP-2α and AP-2β cooperatively orchestrate homeobox gene expression during branchial arch patterning

The evolution of a hinged moveable jaw with variable morphology is considered a major factor behind the successful expansion of the vertebrates. DLX homeobox transcription factors are crucial for establishing the positional code that patterns the mandible, maxilla and intervening hinge domain, but how the genes encoding these proteins are regulated remains unclear. Herein, we demonstrate that the concerted action of the AP-2α and AP-2β transcription factors within the mouse neural crest is essential for jaw patterning. In the absence of these two proteins, the hinge domain is lost and there are alterations in the size and patterning of the jaws correlating with dysregulation of homeobox gene expression, with reduced levels of Emx, Msx and Dlx paralogs accompanied by an expansion of Six1 expression. Moreover, detailed analysis of morphological features and gene expression changes indicate significant overlap with various compound Dlx gene mutants. Together, these findings reveal that the AP-2 genes have a major function in mammalian neural crest development, influencing patterning of the craniofacial skeleton via the DLX code, an effect that has implications for vertebrate facial evolution, as well as for human craniofacial disorders.

Question mark ear deformity-revisited

We report on two unusual cases of Cosman (question mark) ear; both required modifications of the standard techniques for surgical correction. The first patient presented with a unilateral question mark ear and concurrent ear prominence and bulging of the cartilage of the anti-helix. Simultaneous correction was done using a combination of cartilage suturing/scoring (for the prominence and the cartilage bulge) as well as Al-Qattan’s “v-y skin flap-cartilage graft-z-plasty” technique (for the correction of the ear cleft deformity). The second patient had aurico-condylar syndrome with bilateral ear deformity and complete separation of ear lobes from the external ear. Staged transposition followed by Al-Qattan’s technique resulted in a satisfactory outcome.

Level of evidence: Level V, therapeutic study.

Negative regulation of endothelin signaling by SIX1 is required for proper maxillary development

Jaw morphogenesis is a complex event mediated by inductive signals that establish and maintain the distinct developmental domains required for formation of hinged jaws, the defining feature of gnathostomes. The mandibular portion of pharyngeal arch 1 is patterned dorsally by Jagged-Notch signaling and ventrally by endothelin receptor A (EDNRA) signaling. Loss of EDNRA signaling disrupts normal ventral gene expression, the result of which is homeotic transformation of the mandible into a maxilla-like structure. However, loss of Jagged-Notch signaling does not result in significant changes in maxillary development. Here we show in mouse that the transcription factor SIX1 regulates dorsal arch development not only by inducing dorsal Jag1 expression but also by inhibiting endothelin 1 (Edn1) expression in the pharyngeal endoderm of the dorsal arch, thus preventing dorsal EDNRA signaling. In the absence of SIX1, but not JAG1, aberrant EDNRA signaling in the dorsal domain results in partial duplication of the mandible. Together, our results illustrate that SIX1 is the central mediator of dorsal mandibular arch identity, thus ensuring separation of bone development between the upper and lower jaws.

Nkx2.5 regulates endothelin converting enzyme-1 during pharyngeal arch patterning

In gnathostomes, dorsoventral (D-V) patterning of neural crest cells (NCC) within the pharyngeal arches is crucial for the development of hinged jaws. One of the key signals that mediate this process is Endothelin-1 (EDN1). Loss of EDN1 binding to the Endothelin-A receptor (EDNRA) results in loss of EDNRA signaling and subsequent facial birth defects in humans, mice and zebrafish. A rate-limiting step in this crucial signaling pathway is the conversion of immature EDN1 into a mature active form by Endothelin converting enzyme-1 (ECE1). However, surprisingly little is known about how Ece1 transcription is induced or regulated. We show here that Nkx2.5 is required for proper craniofacial development in zebrafish and acts in part by upregulating ece1 expression. Disruption of nkx2.5 in zebrafish embryos results in defects in both ventral and dorsal pharyngeal arch-derived elements, with changes in ventral arch gene expression consistent with a disruption in Ednra signaling. ece1 mRNA rescues the nkx2.5 morphant phenotype, indicating that Nkx2.5 functions through modulating Ece1 expression or function. These studies illustrate a new function for Nkx2.5 in embryonic development and provide new avenues with which to pursue potential mechanisms underlying human facial disorders.

Novel adverse outcome pathways revealed by chemical genetics in a developing marine fish

Crude oil spills are a worldwide ocean conservation threat. Fish are particularly vulnerable to the oiling of spawning habitats, and crude oil causes severe abnormalities in embryos and larvae. However, the underlying mechanisms for these developmental defects are not well understood. Here, we explore the transcriptional basis for four discrete crude oil injury phenotypes in the early life stages of the commercially important Atlantic haddock (Melanogrammus aeglefinus). These include defects in (1) cardiac form and function, (2) craniofacial development, (3) ionoregulation and fluid balance, and (4) cholesterol synthesis and homeostasis. Our findings suggest a key role for intracellular calcium cycling and excitation-transcription coupling in the dysregulation of heart and jaw morphogenesis. Moreover, the disruption of ionoregulatory pathways sheds new light on buoyancy control in marine fish embryos. Overall, our chemical-genetic approach identifies initiating events for distinct adverse outcome pathways and novel roles for individual genes in fundamental developmental processes.

Signalling pathways in trophic skeletal development and morphogenesis: Insights from studies on teleost fish

During the development of the vertebrate feeding apparatus, a variety of complicated cellular and molecular processes participate in the formation and integration of individual skeletal elements. The molecular mechanisms regulating the formation of skeletal primordia and their development into specific morphological structures are tightly controlled by a set of interconnected signalling pathways. Some of these pathways, such as Bmp, Hedgehog, Notch and Wnt, are long known for their pivotal roles in craniofacial skeletogenesis. Studies addressing the functional details of their components and downstream targets, the mechanisms of their interactions with other signals as well as their potential roles in adaptive morphological divergence, are currently attracting considerable attention. An increasing number of signalling pathways that had previously been described in different biological contexts have been shown to be important in the regulation of jaw skeletal development and morphogenesis. In this review, I provide an overview of signalling pathways involved in trophic skeletogenesis emphasizing studies of the most species-rich group of vertebrates, the teleost fish, which through their evolutionary history have undergone repeated episodes of spectacular trophic diversification.

Embryonic expression of endothelins and their receptors in lamprey and frog reveals stem vertebrate origins of complex Endothelin signaling

Neural crest cells (NCCs) are highly patterned embryonic cells that migrate along stereotyped routes to give rise to a diverse array of adult tissues and cell types. Modern NCCs are thought to have evolved from migratory neural precursors with limited developmental potential and patterning. How this occurred is poorly understood. Endothelin signaling regulates several aspects of NCC development, including their migration, differentiation, and patterning. In jawed vertebrates, Endothelin signaling involves multiple functionally distinct ligands (Edns) and receptors (Ednrs) expressed in various NCC subpopulations. To test the potential role of endothelin signaling diversification in the evolution of modern, highly patterned NCC, we analyzed the expression of the complete set of endothelin ligands and receptors in the jawless vertebrate, the sea lamprey (Petromyzon marinus). To better understand ancestral features of gnathostome edn and ednr expression, we also analyzed all known Endothelin signaling components in the African clawed frog (Xenopus laevis). We found that the sea lamprey has a gnathsotome-like complement of edn and ednr duplicates, and these genes are expressed in patterns highly reminiscent of their gnathostome counterparts. Our results suggest that the duplication and specialization of vertebrate Endothelin signaling coincided with the appearance of highly patterned and multipotent NCCs in stem vertebrates.

Viable Ednra Y129F mice feature human mandibulofacial dysostosis with alopecia (MFDA) syndrome due to the homologue mutation

Animal models resembling human mutations are valuable tools to research the features of complex human craniofacial syndromes. This is the first report on a viable dominant mouse model carrying a non-synonymous sequence variation within the endothelin receptor type A gene (Ednra c.386A>T, p.Tyr129Phe) derived by an ENU mutagenesis program. The identical amino acid substitution was reported recently as disease causing in three individuals with the mandibulofacial dysostosis with alopecia (MFDA, OMIM 616367) syndrome. We performed standardized phenotyping of wild-type, heterozygous, and homozygous Ednra^Y129F mice within the German Mouse Clinic. Mutant mice mimic the craniofacial phenotypes of jaw dysplasia, micrognathia, dysplastic temporomandibular joints, auricular dysmorphism, and missing of the squamosal zygomatic process as described for MFDA-affected individuals. As observed in MFDA-affected individuals, mutant Ednra^Y129F mice exhibit hearing impairment in line with strong abnormalities of the ossicles and further, reduction of some lung volumetric parameters. In general, heterozygous and homozygous mice demonstrated inter-individual diversity of expression of the craniofacial phenotypes as observed in MFDA patients but without showing any cleft palates, eyelid defects, or alopecia. Mutant Ednra^Y129F mice represent a valuable viable model for complex human syndromes of the first and second pharyngeal arches and for further studies and analysis of impaired endothelin 1 (EDN1)–endothelin receptor type A (EDNRA) signaling. Above all, Ednra^Y129F mice model the recently published human MFDA syndrome and may be helpful for further disease understanding and development of therapeutic interventions.

Dysregulation of Wnt-Signaling and a Candidate Set of miRNAs Underlie the Effect of Metformin on Neural Crest Cell Development

Neural crest cells (NCC) are a population of epithelial cells that arise from the dorsal tube and undergo epithelial-mesenchymal transition (EMT) eventually generating tissues from peripheral nervous system, melanocytes, craniofacial cartilage, and bone. The antidiabetic drug metformin reportedly inhibits EMT in physiological conditions like cancer and fibrosis. We hypothesize that perturbation of EMT may also contribute to developmental disabilities associated with neural crest (NC) development. To understand the molecular network underlying metformin action during NC formation, we first differentiated murine embryonic stem (ES) cells into NCC and characterized them by demonstrating spatiotemporal regulation of key markers. Metformin treatment prompted a delay in delamination of NCC by inhibiting key markers like Sox-1, Sox-9, HNK-1, and p-75. We then revealed that metformin impedes Wnt axis, a major signaling pathway active during NC formation via DVL-3 inhibition and impairment in nuclear translocation of β-catenin. Concomitantly we identified and tested a candidate set of miRNAs that play a crucial role in NC cell fate determination. Further studies involving loss and gain of function confirmed that NCC specifiers like Sox-1 and Sox-9 are direct targets of miR-200 and miR-145, respectively and that they are essentially modulated by metformin. Our in vitro findings were strongly supported by in vivo studies in zebrafish. Given that metformin is a widely used drug, for the first time we demonstrate that it can induce a delayed onset of developmental EMT during NC formation by interfering with canonical Wnt signaling and mysregulation of miR-145 and miR-200.

Developmental basis of phenotypic integration in two Lake Malawi cichlids

Background

Cichlid fishes from the Rift Lakes of East Africa have undergone the most spectacular adaptive radiations in vertebrate history. Eco-morphological adaptations in lakes Victoria, Malawi and Tanganyika have resulted in a vast array of skull shapes and sizes, yet primary axes of morphological variation are conserved in all three radiations, prominently including the size of the preorbital region of the skull. This conserved pattern suggests that development may constrain the trajectories of cichlid head morphological evolution.

Results

Here, we (1) present a comparative analysis of adult head morphology in two sand-dweller cichlids from Lake Malawi with preorbital size differences representative of the main axis of variation among the three lakes and (2) analyze the ontogeny of shape and size differences by focusing on known developmental modules throughout the head. We find that (1) developmental differences between the two species correlate with known developmental modules; (2) differences in embryonic cartilage development result in phenotypically integrated changes among all bones derived from a single cartilage, while differences in dermal bone development tend to influence isolated regions within a bone; and lastly (3) species-specific morphologies appear in the embryo as subtle differences, which become progressively amplified throughout ontogeny. We propose that this amplification takes place at skeletal growth zones, the locations and shapes of which are patterned during embryogenesis.

Conclusions

This study is the most anatomically comprehensive analysis of the developmental differences underlying cichlid skull evolution in the Rift Lakes of East Africa. The scale of our analysis reveals previously unnoticed correlations between developmental modules and patterns of phenotypic integration. We propose that the primary axes of morphological variation among East African cichlid adaptive radiations are constrained by the hierarchical modularity of the teleost head skeleton.

Postotic and preotic cranial neural crest cells differently contribute to thyroid development

Thyroid development and formation vary among species, but in most species the thyroid morphogenesis consists of five stages: specification, budding, descent, bilobation and folliculogenesis. The detailed mechanisms of these stages have not been fully clarified. During early development, the cranial neural crest (CNC) contributes to the thyroid gland. The removal of the postotic CNC (corresponding to rhombomeres 6, 7 and 8, also known as the cardiac neural crest) results in abnormalities of the cardiovascular system, thymus, parathyroid glands, and thyroid gland. To investigate the influence of the CNC on thyroid bilobation process, we divided the CNC into two regions, the postotic CNC and the preotic CNC (from the mesencephalon to rhombomere 5) regions and examined. We found that preotic CNC-ablated embryos had a unilateral thyroid lobe, and confirmed the presence of a single lobe or the absence of lobes in postotic CNC-ablated chick embryos. The thyroid anlage in each region-ablated embryos was of a normal size at the descent stage, but at a later stage, the thyroid in preotic CNC-ablated embryos was of a normal size, conflicting with a previous report in which the thyroid was reduced in size in the postotic CNC-ablated embryos. The postotic CNC cells differentiated into connective tissues of the thyroid in quail-to-chick chimeras. In contrast, the preotic CNC cells did not differentiate into connective tissues of the thyroid. We found that preotic CNC cells encompassed the thyroid anlage from the specification stage to the descent stage. Finally, we found that endothelin-1 and endothelin type A receptor-knockout mice and bosentan (endothelin receptor antagonist)-treated chick embryos showed bilobation anomalies that included single-lobe formation. Therefore, not only the postotic CNC, but also the preotic CNC plays an important role in thyroid morphogenesis.

Regulating Craniofacial Development at the 3' End: MicroRNAs and Their Function in Facial Morphogenesis

Defects in craniofacial development represent a majority of observed human birth defects, occurring at a rate as high as 1:800 live births. These defects often occur due to changes in neural crest cell (NCC) patterning and development and can affect non-NCC-derived structures due to interactions between NCCs and the surrounding cell types. Proper craniofacial development requires an intricate array of gene expression networks that are tightly controlled spatiotemporally by a number of regulatory mechanisms. One of these mechanisms involves the action of microRNAs (miRNAs), a class of noncoding RNAs that repress gene expression by binding to miRNA recognition sequences typically located in the 3' UTR of target mRNAs. Recent evidence illustrates that miRNAs are crucial for vertebrate facial morphogenesis, with changes in miRNA expression leading to facial birth defects, including some in complex human syndromes such as 22q11 (DiGeorge Syndrome). In this review, we highlight the current understanding of miRNA biogenesis, the roles of miRNAs in overall craniofacial development, the impact that loss of miRNAs has on normal development and the requirement for miRNAs in the development of specific craniofacial structures, including teeth. From these studies, it is clear that miRNAs are essential for normal facial development and morphogenesis, and a potential key in establishing new paradigms for repair and regeneration of facial defects.

Candidate Gene Analyses of Skeletal Variation in Malocclusion

This study evaluated associations between craniofacial candidate genes and skeletal variation in patients with malocclusion. Lateral cephalometric radiographs of 269 untreated adults with skeletal classes I, II, and III malocclusion were digitized with 14 landmarks. Two-dimensional coordinates were analyzed using Procrustes fit and principal component (PC) analysis to generate continuous malocclusion phenotypes. Skeletal class classifications (I, II, or III) were used as a categorical phenotype. Individuals were genotyped for 198 single-nucleotide polymorphisms (SNPs) in 71 craniofacial genes and loci. Phenotype-genotype associations were tested via multivariate linear regression for continuous phenotypes and multinomial logistic regression for skeletal malocclusion class. PC analysis resulted in 4 principal components (PCs) explaining 69% of the total skeletal facial variation. PC1 explained 32.7% of the variation and depicted vertical discrepancies ranging from skeletal deep to open bites. PC1 was associated with a SNP near PAX5 (P = 0.01). PC2 explained 21.7% and captured horizontal maxillomandibular discrepancies. PC2 was associated with SNPs upstream of SNAI3 (P = 0.0002) and MYO1H (P = 0.006). PC3 explained 8.2% and captured variation in ramus height, body length, and anterior cranial base orientation. PC3 was associated with TWIST1 (P = 0.000076). Finally, PC4 explained 6.6% and detected variation in condylar inclination as well as symphysis projection. PC4 was associated with PAX7 (P = 0.007). Furthermore, skeletal class II risk increased relative to class I with the minor alleles of SNPs in FGFR2 (odds ratio [OR] = 2.1, P = 0.004) and declined with SNPs in EDN1 (OR = 0.5, P = 0.007). Conversely, skeletal class III risk increased versus class I with SNPs in FGFR2 (OR 2.2, P = 0.005) and COL1A1 (OR = 2.1, P = 0.008) and declined with SNPs in TBX5 (OR = 0.5, P = 0.014). PAX5, SNAI3, MYO1H, TWIST1, and PAX7 are associated with craniofacial skeletal variation among patients with malocclusion, while FGFR2, EDN1, TBX5, and COL1A1 are associated with type of skeletal malocclusion.