Combinatorial roles for BMPs and Endothelin 1 in patterning the dorsal-ventral axis of the craniofacial skeleton

Integration of the sensory and skeletal systems: A classical perspective on neuromast-bone interactions

The question of how morphological features of the vertebrate face evolve over deep geologic time can be inferred, in part, from the fossil record. However, paleontological evidence does not inform the developmental processes governing facial form and patterning. Insight to this fundamental feature of evolution, therefore, must come from empirical studies in extant model systems. An essential developmental feature of vertebrate facial skeleton evolution is the integration between sensation and bone. In the classic literature, a longstanding observation is that the positions of certain facial bone primary ossification centers (i.e., bone initiation sites) co-localize to the precise positions of sensory organs called canal neuromasts. These organs can be found in stereotypical positions along the head and flank of extant and extinct specimens. Some early researchers proposed that a developmental relationship couples these sensory organs to underlying bones (or bony structures). The evidence to support this relationship, however, was sparse and many remained skeptical - or outright disputed - the existence of such an integrated relationship. Several researchers instead argued that co-localization of neuromasts and bone were merely coincidental. Here, we review the classic literature, summarizing contrasting claims of several authors, stretching back to 1882. This revealed that speculation around an integrative sensory-skeletal relationship fell into three general categories: the presence of a developmental relationship between the structures, the absence of a relationship, or the existence of an indirect relationship between the structures. In summarizing the classic literature, we present further evidence in support of a developmental relationship illuminated by contemporary findings. We propose sensory-skeletal integration is a key developmental mechanism that could be modeled in paleontological research to further inform processes of cranial evolution over deep geological time.

Nkx2.7 is a conserved regulator of craniofacial development

Craniofacial malformations arise from developmental defects in the head, face, and neck with phenotypes such as 22q11.2 deletion syndrome illustrating a developmental link between cardiovascular and craniofacial morphogenesis. NKX2-5 is a key cardiac transcription factor associated with congenital heart disease and mouse models of Nkx2-5 deficiency highlight roles in cardiac development. In zebrafish, nkx2.5 and nkx2.7 are paralogues in the NK4 family expressed in cardiomyocytes and pharyngeal arches. Despite shared cellular origins of cardiac and craniofacial tissues, the function of NK4 factors in head and neck patterning has not been elucidated. Molecular evolutionary analysis of NK4 genes shows that nkx2.5 and nkx2.7 are ohnologs resulting from whole genome duplication events. Nkx2.7 serves as a previously unappreciated regulator of branchiomeric muscle and cartilage formation for which nkx2.5 cannot fully compensate. Mechanistically, our results highlight that Nkx2.7 patterns the cranial neural crest and functions upstream of Endothelin1 to inhibit Notch signals. Together, our studies shed light on an evolutionarily conserved Nkx transcription factor with unique functions in vertebrate craniofacial development, advancing our understanding of congenital head and neck deformities.

The Gq/11 family of Gα subunits is necessary and sufficient for lower jaw development

Vertebrate jaw development is coordinated by highly conserved ligand-receptor systems such as the peptide ligand Endothelin 1 (Edn1) and Endothelin receptor type A (Ednra), which are required for patterning of lower jaw structures. The Edn1/Ednra signaling pathway establishes the identity of lower jaw progenitor cells by regulating expression of numerous patterning genes, but the intracellular signaling mechanisms linking receptor activation to gene regulation remain poorly understood. As a first step towards elucidating this mechanism, we examined the function of the Gq/11 family of Gα subunits in zebrafish using pharmacological inhibition and genetic ablation of Gq/11 activity, and transgenic induction of a constitutively active Gq protein in edn1-/- embryos. Genetic loss of Gq/11 activity fully recapitulated the edn1-/- phenotype, with genes encoding G11 being most essential. Furthermore, inducing Gq activity in edn1-/- embryos not only restored Edn1/Ednra-dependent jaw structures and gene expression signatures but also caused homeosis of the upper jaw structure into a lower jaw-like structure. These results indicate that Gq/11 is necessary and sufficient to mediate the lower jaw patterning mechanism for Ednra in zebrafish.

Loss of Atoh8 Affects Neurocranial and Axial Skeleton Development in Zebrafish

BACKGROUND

The basic helix-loop-helix (bHLH) transcription factor atonal homologue 8 (Atoh8) has been implicated in various developmental and physiological processes by means of transient knockdown and conditional knockout approaches in zebrafish, chick and mouse. Despite its demonstrated involvement in multiple tissues, the role of Atoh8 remains elusive in zebrafish. A recent permanent knockout study in zebrafish investigated the role of Atoh8 on the background of previous morpholino studies which demonstrated various developmental defects but could not find any of the morpholino-based effects in the mutant. In mice, a knockout study demonstrated involvement of the transcription factor in skeletal development, showing that disruption of the atoh8 gene results in reduction of skeletal size. We investigated a mutant fish line generated using clustered regularly interspaced short palindromic repeats (CRISPR)/CRISPR-associated protein 9 (Cas9) (CRISPR/Cas9)-technology for possible phenotypic effects on zebrafish skeletogenesis.

METHODS

Here, we present a CRISPR/Cas9-generated atoh8 permanent zebrafish mutant and investigate the phenotypic effects of the knockout on the developing zebrafish craniofacial and axial skeleton. We investigated the expression pattern of the gene in wildtype and conducted detailed morphometric analysis for a variety of bone and cartilage elements of the developing skeleton at 12 days post fertilisation (dpf) in zebrafish siblings from a heterozygous mating using detailed morphometric measurements and statistical analysis of the results.

RESULTS

Homozygous mutants are viable into late adulthood and show no overt morphological phenotype. Despite the prominent appearance of atoh8 signal in various embryonic and larval craniofacial and axial skeletal structures, detailed morphometric analysis revealed only subtle phenotypic effects of the mutation on skeletal development in zebrafish. We found the formation of the orbital cartilages of the developing neurocranium and the progress of chordacentra mineralisation to be negatively affected by loss of the transcription factor.

CONCLUSIONS

Despite the very subtle phenotypic effect of our mutation, we were able to show involvement of atoh8 in the skeletal development of zebrafish. We attribute the mild phenotype to a compensatory mechanism induced by nonsense-mediated degradation of messenger ribonucleic acid (mRNA) as suggested in the recent literature. The effect of atoh8-disruption on zebrafish skeletal development suggests that the loss of atoh8 cannot be compensated for at interfaces where more than one embryonic cell lineage contributes to bone and cartilage formation.

Bone morphogenetic protein signaling pathway- Ethanol interactions disrupt palate formation independent of gata3

Fetal Alcohol Spectrum Disorders (FASD) describes a wide array of neurological defects and craniofacial malformations, associated with ethanol teratogenicity. While there is growing evidence for a genetic component to FASD, little is known of the genes underlying these ethanol-induced defects. Along with timing and dosage, genetic predispositions may help explain the variability within FASD. From a screen for gene-ethanol interactions, we found that mutants for Bmp signaling components are ethanol-sensitive leading to defects in the zebrafish palate. Loss of Bmp signaling results in reductions in gata3 expression in the maxillary domain of the neural crest in the 1st pharyngeal arch, leading to palate defects while upregulation of human GATA3 rescues these defects. Here, we show that ethanol-treated Bmp mutants exhibit misshaped and/or broken trabeculae. Surprisingly, up regulation of GATA3 does not rescue ethanol-induced palate defects and gata3 expression was not altered in ethanol-treated Bmp mutants or dorsomorphin-treated larvae. Timing of ethanol sensitivity shows that Bmp mutants are ethanol sensitive from 10 to 18 hours post-fertilization (hpf), prior to Bmp's regulation of gata3 in palate formation. This is consistent with our previous work with dorsomorphin-dependent knock down of Bmp signaling from 10 to 18 hpf disrupting endoderm formation and subsequent jaw development. Overall, this suggests that ethanol disrupts Bmp-dependent palate development independent of and earlier than the role of gata3 in palate formation by disrupting epithelial development. Ultimately, these data demonstrate that zebrafish is a useful model to identify and characterize gene-ethanol interactions and this work will directly inform our understanding of FASD.

Bone Morphogenetic Protein signaling pathway - ethanol interactions disrupt palate formation independent of gata3

Fetal Alcohol Spectrum Disorders (FASD) describes a wide array of neurological defects and craniofacial malformations, associated with ethanol teratogenicity. While there is growing evidence for a genetic component to FASD, little is known of the genes underlying these ethanol-induced defects. Along with timing and dosage, genetic predispositions may help explain the variability within FASD. From a screen for gene-ethanol interactions, we found that mutants for Bmp signaling components are ethanol-sensitive leading to defects in the zebrafish palate. Loss of Bmp signaling results in reductions in gata3 expression in the maxillary domain of the neural crest in the 1st pharyngeal arch, leading to palate defects while upregulation of human GATA3 rescues these defects. Here, we show that ethanol-treated Bmp mutants exhibit misshaped and/or broken trabeculae. Surprisingly, up regulation of GATA3 does not rescue ethanol-induced palate defects and gata3 expression was not altered in ethanol-treated Bmp mutants or dorsomorphin-treated larvae. Timing of ethanol sensitivity shows that Bmp mutants are ethanol sensitive from 10-18 hours post-fertilization (hpf), prior to Bmp's regulation of gata3 in palate formation. This is consistent with our previous work with dorsomorphin-dependent knock down of Bmp signaling from 10-18 hpf disrupting endoderm formation and subsequent jaw development. Overall, this suggests that ethanol disrupts Bmp-dependent palate development independent of and earlier than the role of gata3 in palate formation by disrupting epithelial development. Ultimately, these data demonstrate that zebrafish is a useful model to identify and characterize gene-ethanol interactions and this work will directly inform our understanding of FASD.

Highlights

Bmp pathway mutants are ethanol sensitive resulting in palate defects. Ethanol disrupts Bmp-dependent palate development independent of gata3 . Timing of ethanol sensitivity suggests ethanol disrupts Bmp-dependent epithelial morphogenesis.

The Gq/11 family of Gα subunits is necessary and sufficient for lower jaw development

Vertebrate jaw development is coordinated by highly conserved ligand-receptor systems such as the peptide ligand Endothelin 1 (Edn1) and Endothelin receptor type A (Ednra), which are required for patterning of lower jaw structures. The Edn1/Ednra signaling pathway establishes the identity of lower jaw progenitor cells by regulating expression of numerous patterning genes, but the intracellular signaling mechanisms linking receptor activation to gene regulation remain poorly understood. As a first step towards elucidating this mechanism, we examined the function of the Gq/11 family of Gα subunits in zebrafish using pharmacological inhibition and genetic ablation of Gq/11 activity and transgenic induction of a constitutively active Gq protein in edn1 -/- embryos. Genetic loss of Gq/11 activity fully recapitulated the edn1 -/- phenotype, with genes encoding G11 being most essential. Furthermore, inducing Gq activity in edn1 -/- embryos not only restored Edn1/Ednra-dependent jaw structures and gene expression signatures but also caused homeosis of the upper jaw structure into a lower jaw-like structure. These results indicate that Gq/11 is necessary and sufficient to mediate the lower jaw patterning mechanism for Ednra in zebrafish.

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

Humans and other jawed vertebrates rely heavily on their craniofacial skeleton for eating, breathing, and communicating. As such, it is vital that the elements of the craniofacial skeleton develop properly during embryogenesis to ensure a high quality of life and evolutionary fitness. Indeed, craniofacial abnormalities, including cleft palate and craniosynostosis, represent some of the most common congenital abnormalities in newborns. Like many other organ systems, the development of the craniofacial skeleton is complex, relying on specification and migration of the neural crest, patterning of the pharyngeal arches, and morphogenesis of each skeletal element into its final form. These processes must be carefully coordinated and integrated. One way this is achieved is through the spatial and temporal deployment of cell signaling pathways. Recent studies conducted using the zebrafish model underscore the importance of the Transforming Growth Factor Beta (TGF-β) and Bone Morphogenetic Protein (BMP) pathways in craniofacial development. Although both pathways contain similar components, each pathway results in unique outcomes on a cellular level. In this review, we will cover studies conducted using zebrafish that show the necessity of these pathways in each stage of craniofacial development, starting with the induction of the neural crest, and ending with the morphogenesis of craniofacial elements. We will also cover human skeletal and craniofacial diseases and malformations caused by mutations in the components of these pathways (e.g., cleft palate, craniosynostosis, etc.) and the potential utility of zebrafish in studying the etiology of these diseases. We will also briefly cover the utility of the zebrafish model in joint development and biology and discuss the role of TGF-β/BMP signaling in these processes and the diseases that result from aberrancies in these pathways, including osteoarthritis and multiple synostoses syndrome. Overall, this review will demonstrate the critical roles of TGF-β/BMP signaling in craniofacial development and show the utility of the zebrafish model in development and disease.

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.

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.

Functional comparison of human ACVR1 and zebrafish Acvr1l FOP-associated variants in embryonic zebrafish

BACKGROUND

Fibrodysplasia ossificans progressiva (FOP), a rare disease characterized by progressive heterotopic ossification of muscle and connective tissues, is caused by autosomal dominant activating mutations in the type I receptor, ACVR1/ALK2. The classic human FOP variant, ACVR1R206H , shows increased bone morphogenetic protein (BMP) signaling and activation by activins.

RESULTS

Here, we performed in vivo functional characterization of human ACVR1R206H and orthologous zebrafish Acvr1lR203H using early embryonic zebrafish dorsoventral patterning as a phenotypic readout for receptor activity. Our results showed that human ACVR1R206H and zebrafish Acvr1lR203H exhibit functional differences in early embryonic zebrafish, and that human ACVR1R206H retained its signaling activity in the absence of a ligand-binding domain (LBD). We also showed, for the first time, that zebrafish Acvr2ba/Acvr2bb receptors are required for human ACVR1R206H signaling in early embryonic zebrafish.

CONCLUSIONS

Together, these data provide new insight into ACVR1R206H signaling pathways that may facilitate the design of new and effective therapies for FOP patients.

Myoneurin regulates BMP signaling by competing with Ppm1a for Smad binding

A delicate balance of BMP activity is critical for tissue formation and organogenesis. However, the mechanical molecular details in ensuring the proper duration and intensity of BMP signaling have yet to be fully elucidated. Here, we identified a zebrafish mutant with a disrupted gene encoding for the BTB/POZ and zinc finger protein myoneurin (Mynn). mynn^−/− mutants exhibited severe loss of pharyngeal cartilage elements, owing to poor proliferation, blocked differentiation, and low viability of cranial neural crest cells. Depletion of mynn in both zebrafish embryos and mammalian cells led to a reduction of the BMP signal activity. Mechanistically, Mynn interacts with Smad proteins in the nucleus, thereby disrupting the association between Smad protein and the phosphatase Ppm1a. Ultimately, this interaction prevents Smad dephosphorylation. More broadly, our findings may provide a new strategy to balance BMP signal activity via competitive binding of Mynn and Ppm1a to Smad proteins during pharyngeal cartilage formation.

•

mynn gene is essential for pharyngeal cartilage development

•

mynn is required for the proliferation, differentiation, and survival of the CNCCs

•

Mynn has an evolutionarily conserved function in supporting BMP signal

•

Mynn maintains BMP signal activity by competing with Ppm1a for Smad binding

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.

Metal transporter Slc30a1 controls pharyngeal neural crest differentiation via the zinc-Snai2-Jag1 cascade

The pharyngeal arch (PA) is a neural crest (NC)-derived organ that is transiently developed during embryogenesis and is required for the subsequent development of various tissues. However, the role of zinc during PA differentiation from NC progenitor cells is unknown. Here, we found that the metal transporters Slc30a1a and Slc30a1b mediate zinc homeostasis during PA differentiation. Slc30a1-deficient zebrafish develop zinc accumulation in NC cells, with increased expression of stemness markers and PA dorsal genes, and SMART-seq analyses revealed that the genes snai2 and jag1b may serve as downstream targets. Furthermore, functional studies showed that knocking down either snai2 or jag1b rescues PA development in Slc30a1-deficient zebrafish. Notably, we identified the double zinc-finger domain in the transcription factor Snai2 as a zinc-responsive element that regulates jag1b expression. Our findings indicate that the Slc30a1/zinc-snai2-jag1b axis is an essential regulatory network controlling PA differentiation, shedding new light on the function of zinc homeostasis in maintaining NC cell stemness and multipotency in vertebrates.

New locus underlying auriculocondylar syndrome (ARCND): 430 kb duplication involving TWIST1 regulatory elements

Background

Auriculocondylar syndrome (ARCND) is a rare genetic disease that affects structures derived from the first and second pharyngeal arches, mainly resulting in micrognathia and auricular malformations. To date, pathogenic variants have been identified in three genes involved in the EDN1-DLX5/6 pathway (PLCB4, GNAI3 and EDN1) and some cases remain unsolved. Here we studied a large unsolved four-generation family.

Methods

We performed linkage analysis, resequencing and Capture-C to investigate the causative variant of this family. To test the pathogenicity of the CNV found, we modelled the disease in patient craniofacial progenitor cells, including induced pluripotent cell (iPSC)-derived neural crest and mesenchymal cells.

Results

This study highlights a fourth locus causative of ARCND, represented by a tandem duplication of 430 kb in a candidate region on chromosome 7 defined by linkage analysis. This duplication segregates with the disease in the family (LOD score=2.88) and includes HDAC9, which is located over 200 kb telomeric to the top candidate gene TWIST1. Notably, Capture-C analysis revealed multiple cis interactions between the TWIST1 promoter and possible regulatory elements within the duplicated region. Modelling of the disease revealed an increased expression of HDAC9 and its neighbouring gene, TWIST1, in neural crest cells. We also identified decreased migration of iPSC-derived neural crest cells together with dysregulation of osteogenic differentiation in iPSC-affected mesenchymal stem cells.

Conclusion

Our findings support the hypothesis that the 430 kb duplication is causative of the ARCND phenotype in this family and that deregulation of TWIST1 expression during craniofacial development can contribute to the phenotype.

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.

The Mandibular and Hyoid Arches-From Molecular Patterning to Shaping Bone and Cartilage

The mandibular and hyoid arches collectively make up the facial skeleton, also known as the viscerocranium. Although all three germ layers come together to assemble the pharyngeal arches, the majority of tissue within viscerocranial skeletal components differentiates from the neural crest. Since nearly one third of all birth defects in humans affect the craniofacial region, it is important to understand how signalling pathways and transcription factors govern the embryogenesis and skeletogenesis of the viscerocranium. This review focuses on mouse and zebrafish models of craniofacial development. We highlight gene regulatory networks directing the patterning and osteochondrogenesis of the mandibular and hyoid arches that are actually conserved among all gnathostomes. The first part of this review describes the anatomy and development of mandibular and hyoid arches in both species. The second part analyses cell signalling and transcription factors that ensure the specificity of individual structures along the anatomical axes. The third part discusses the genes and molecules that control the formation of bone and cartilage within mandibular and hyoid arches and how dysregulation of molecular signalling influences the development of skeletal components of the viscerocranium. In conclusion, we notice that mandibular malformations in humans and mice often co-occur with hyoid malformations and pinpoint the similar molecular machinery controlling the development of mandibular and hyoid arches.

Variation in phenotypes from a Bmp-Gata3 genetic pathway is modulated by Shh signaling

We sought to understand how perturbation of signaling pathways and their targets generates variable phenotypes. In humans, GATA3 associates with highly variable defects, such as HDR syndrome, microsomia and choanal atresia. We previously characterized a zebrafish point mutation in gata3 with highly variable craniofacial defects to the posterior palate. This variability could be due to residual Gata3 function, however, we observe the same phenotypic variability in gata3 null mutants. Using hsp:GATA3-GFP transgenics, we demonstrate that Gata3 function is required between 24 and 30 hpf. At this time maxillary neural crest cells fated to generate the palate express gata3. Transplantation experiments show that neural crest cells require Gata3 function for palatal development. Via a candidate approach, we determined if Bmp signaling was upstream of gata3 and if this pathway explained the mutant's phenotypic variation. Using BRE:d2EGFP transgenics, we demonstrate that maxillary neural crest cells are Bmp responsive by 24 hpf. We find that gata3 expression in maxillary neural crest requires Bmp signaling and that blocking Bmp signaling, in hsp:DN-Bmpr1a-GFP embryos, can phenocopy gata3 mutants. Palatal defects are rescued in hsp:DN-Bmpr1a-GFP;hsp:GATA3-GFP double transgenic embryos, collectively demonstrating that gata3 is downstream of Bmp signaling. However, Bmp attenuation does not alter phenotypic variability in gata3 loss-of-function embryos, implicating a different pathway. Due to phenotypes observed in hypomorphic shha mutants, the Sonic Hedgehog (Shh) pathway was a promising candidate for this pathway. Small molecule activators and inhibitors of the Shh pathway lessen and exacerbate, respectively, the phenotypic severity of gata3 mutants. Importantly, inhibition of Shh can cause gata3 haploinsufficiency, as observed in humans. We find that gata3 mutants in a less expressive genetic background have a compensatory upregulation of Shh signaling. These results demonstrate that the level of Shh signaling can modulate the phenotypes observed in gata3 mutants.

Conserved and unique transcriptional features of pharyngeal arches in the skate (Leucoraja erinacea) and evolution of the jaw

The origin of the jaw is a long-standing problem in vertebrate evolutionary biology. Classical hypotheses of serial homology propose that the upper and lower jaw evolved through modifications of dorsal and ventral gill arch skeletal elements, respectively. If the jaw and gill arches are derived members of a primitive branchial series, we predict that they would share common developmental patterning mechanisms. Using candidate and RNAseq/differential gene expression analyses, we find broad conservation of dorsoventral (DV) patterning mechanisms within the developing mandibular, hyoid, and gill arches of a cartilaginous fish, the skate (Leucoraja erinacea). Shared features include expression of genes encoding members of the ventralizing BMP and endothelin signaling pathways and their effectors, the joint markers nkx3.2 and gdf5 and prochondrogenic transcription factor barx1, and the dorsal territory marker pou3f3. Additionally, we find that mesenchymal expression of eya1/six1 is an ancestral feature of the mandibular arch of jawed vertebrates, whereas differences in notch signaling distinguish the mandibular and gill arches in skate. Comparative transcriptomic analyses of mandibular and gill arch tissues reveal additional genes differentially expressed along the DV axis of the pharyngeal arches, including scamp5 as a novel marker of the dorsal mandibular arch, as well as distinct transcriptional features of mandibular and gill arch muscle progenitors and developing gill buds. Taken together, our findings reveal conserved patterning mechanisms in the pharyngeal arches of jawed vertebrates, consistent with serial homology of their skeletal derivatives, as well as unique transcriptional features that may underpin distinct jaw and gill arch morphologies.

Endothelial struts enable the generation of large lumenized blood vessels de novo

De novo blood vessel formation occurs through coalescence of endothelial cells (ECs) into a cord-like structure, followed by lumenization either through cell-^1-3 or cord-hollowing^4-7. Vessels generated in this manner are restricted in diameter to one or two ECs, and these models fail to explain how vasculogenesis can form large-diameter vessels. Here, we describe a model for large vessel formation that does not require a cord-like structure or a hollowing step. In this model, ECs coalesce into a network of struts in the future lumen of the vessel, a process dependent upon bone morphogenetic protein signalling. The vessel wall forms around this network and consists initially of only a few patches of ECs. To withstand external forces and to maintain the shape of the vessel, strut formation traps erythrocytes into compartments to form a rigid structure. Struts gradually prune and ECs from struts migrate into and become part of the vessel wall. Experimental severing of struts resulted in vessel collapse, disturbed blood flow and remodelling defects, demonstrating that struts enable the patency of large vessels during their formation.

Pharmacological Manipulation of Early Zebrafish Skeletal Development Shows an Important Role for Smad9 in Control of Skeletal Progenitor Populations

Osteoporosis and other conditions associated with low bone density or quality are highly prevalent, are increasing as the population ages and with increased glucocorticoid use to treat conditions with elevated inflammation. There is an unmet need for therapeutics which can target skeletal precursors to induce osteoblast differentiation and osteogenesis. Genes associated with high bone mass represent interesting targets for manipulation, as they could offer ways to increase bone density. A damaging mutation in SMAD9 has recently been associated with high bone mass. Here we show that Smad9 labels groups of osteochondral precursor cells, which are not labelled by the other Regulatory Smads: Smad1 or Smad5. We show that Smad9⁺ cells are proliferative, and that the Smad9⁺ pocket expands following osteoblast ablation which induced osteoblast regeneration. We further show that treatment with retinoic acid, prednisolone, and dorsomorphin all alter Smad9 expression, consistent with the effects of these drugs on the skeletal system. Taken together these results demonstrate that Smad9⁺ cells represent an undifferentiated osteochondral precursor population, which can be manipulated by commonly used skeletal drugs. We conclude that Smad9 represents a target for future osteoanabolic therapies.

Potential of zebrafish as a model to characterise MicroRNA profiles in mechanically mediated joint degeneration

Mechanically mediated joint degeneration and cartilage dyshomeostasis is implicated in highly prevalent diseases such as osteoarthritis. Increasingly, MicroRNAs are being associated with maintaining the normal state of cartilage, making them an exciting and potentially key contributor to joint health and disease onset. Here, we present a summary of current in vitro and in vivo models which can be used to study the role of mechanical load and MicroRNAs in joint degeneration, including: non-invasive murine models of PTOA, surgical models which involve ligament transection, and unloading models based around immobilisation of joints or removal of load from the joint through suspension. We also discuss how zebrafish could be used to advance this field, namely through the availability of transgenic lines relevant to cartilage homeostasis and the ability to accurately map strain through the cartilage, enabling the response of downstream MicroRNA targets to be followed dynamically at a cellular level in areas of high and low strain.

BMP2 signalling activation enhances bone metastases of non-small cell lung cancer

Distant metastases occur when non‐small cell lung cancer (NSCLC) is at late stages. Bone metastasis is one of the most frequent metastases of NSCLC and leads to poor prognosis. It has been reported that high expression of BMP2 in NSCLC correlates with poor survival, but whether BMP2 contributes to NSCLC bone metastasis remains largely unknown. The activation of BMP signalling is found in metastatic bone tumours of mice Lewis lung carcinoma and predicts poor survival in human NSCLC. BMP2 signalling activation can enhance bone metastasis of Lewis lung carcinoma. Moreover, BMP2 secreted by stroma fibroblasts can promote the migration and invasion of NSCLC cells. Besides, in combination with pre‐osteoblast and LLCs, BMP2 could enhance the differentiation of macrophages into osteoclasts to play roles in the osteolytic mechanism of NSCLC bone metastasis. Interestingly, NSCLC cells can also enrich BMP2 to pre‐osteoblasts to function in the osteoblastic mechanism. Our results firstly demonstrate the detailed mechanisms about what roles BMP2 signalling play in enhancing NSCLC bone metastases. These findings provide a new potential therapy choice for preventing bone metastases of NSCLC via the inhibition of BMP2 signalling.

A Rare Mutation in SMAD9 Associated With High Bone Mass Identifies the SMAD-Dependent BMP Signaling Pathway as a Potential Anabolic Target for Osteoporosis

Novel anabolic drug targets are needed to treat osteoporosis. Having established a large national cohort with unexplained high bone mass (HBM), we aimed to identify a novel monogenic cause of HBM and provide insight into a regulatory pathway potentially amenable to therapeutic intervention. We investigated a pedigree with unexplained HBM in whom previous sequencing had excluded known causes of monogenic HBM. Whole exome sequencing identified a rare (minor allele frequency 0.0023), highly evolutionarily conserved missense mutation in SMAD9 (c.65T>C, p.Leu22Pro) segregating with HBM in this autosomal dominant family. The same mutation was identified in another two unrelated individuals both with HBM. In silico protein modeling predicts the mutation severely disrupts the MH1 DNA-binding domain of SMAD9. Affected individuals have bone mineral density (BMD) Z-scores +3 to +5, mandible enlargement, a broad frame, torus palatinus/mandibularis, pes planus, increased shoe size, and a tendency to sink when swimming. Peripheral quantitative computed tomography (pQCT) measurement demonstrates increased trabecular volumetric BMD and increased cortical thickness conferring greater predicted bone strength; bone turnover markers are low/normal. Notably, fractures and nerve compression are not found. Both genome-wide and gene-based association testing involving estimated BMD measured at the heel in 362,924 white British subjects from the UK Biobank Study showed strong associations with SMAD9 (PGWAS = 6 × 10-16 ; PGENE = 8 × 10-17 ). Furthermore, we found Smad9 to be highly expressed in both murine cortical bone-derived osteocytes and skeletal elements of zebrafish larvae. Our findings support SMAD9 as a novel HBM gene and a potential novel osteoanabolic target for osteoporosis therapeutics. SMAD9 is thought to inhibit bone morphogenetic protein (BMP)-dependent target gene transcription to reduce osteoblast activity. Thus, we hypothesize SMAD9 c.65T>C is a loss-of-function mutation reducing BMP inhibition. Lowering SMAD9 as a potential novel anabolic mechanism for osteoporosis therapeutics warrants further investigation. © 2019 The Authors. Journal of Bone and Mineral Research published by American Society for Bone and Mineral Research.

Zebrafish: A Resourceful Vertebrate Model to Investigate Skeletal Disorders

Animal models are essential tools for addressing fundamental scientific questions about skeletal diseases and for the development of new therapeutic approaches. Traditionally, mice have been the most common model organism in biomedical research, but their use is hampered by several limitations including complex generation, demanding investigation of early developmental stages, regulatory restrictions on breeding, and high maintenance cost. The zebrafish has been used as an efficient alternative vertebrate model for the study of human skeletal diseases, thanks to its easy genetic manipulation, high fecundity, external fertilization, transparency of rapidly developing embryos, and low maintenance cost. Furthermore, zebrafish share similar skeletal cells and ossification types with mammals. In the last decades, the use of both forward and new reverse genetics techniques has resulted in the generation of many mutant lines carrying skeletal phenotypes associated with human diseases. In addition, transgenic lines expressing fluorescent proteins under bone cell- or pathway- specific promoters enable in vivo imaging of differentiation and signaling at the cellular level. Despite the small size of the zebrafish, many traditional techniques for skeletal phenotyping, such as x-ray and microCT imaging and histological approaches, can be applied using the appropriate equipment and custom protocols. The ability of adult zebrafish to remodel skeletal tissues can be exploited as a unique tool to investigate bone formation and repair. Finally, the permeability of embryos to chemicals dissolved in water, together with the availability of large numbers of small-sized animals makes zebrafish a perfect model for high-throughput bone anabolic drug screening. This review aims to discuss the techniques that make zebrafish a powerful model to investigate the molecular and physiological basis of skeletal disorders.

Novel Ethanol-Sensitive Mutants Identified in an F3 Forward Genetic Screen

Background

Fetal alcohol spectrum disorders (FASD) collectively refer to all deleterious outcomes due to prenatal alcohol exposures. Alterations to the face are common phenotypes in FASD. While alcohol exposure is the underlying cause of FASD, many variables modify the outcomes of such exposures. Genetic risk is one such variable, yet we still have a limited understanding of the nature of the genetic loci mediating susceptibility to FASD.

Methods

We employed ENU‐based random mutagenesis in zebrafish to identify mutations that enhanced the teratogenicity of ethanol (EtOH). F3 embryos obtained from 126 inbred F2 families were exposed to 1% EtOH in the medium (approximately 41 mM tissue levels). Zebrafish stained with Alcian Blue and Alizarin Red were screened for qualitative alterations to the craniofacial skeleton between 4 and 7 days postfertilization (dpf).

Results

In all, we recovered 6 EtOH‐sensitive mutants, 5 from the genetic screen itself and one as a background mutation in one of our wild‐type lines. Each mutant has a unique EtOH‐induced phenotype relative to the other mutant lines. All but 1 mutation appears to be recessive in nature, and only 1 mutant, au29, has apparent craniofacial defects in the absence of EtOH. To validate the genetic screen, we genetically mapped au29 and found that it carries a mutation in a previously uncharacterized gene, si:dkey‐88l16.3.

Conclusions

The phenotypes of these EtOH‐sensitive mutants differ from those in previous characterizations of gene–EtOH interactions. Thus, each mutant is likely to provide novel insights into EtOH teratogenesis. Given that most of these mutants only have craniofacial defects in the presence of EtOH and our mapping of au29, it is also likely that many of the mutants will be previously uncharacterized. Collectively, our findings point to the importance of unbiased genetic screens in the identification, and eventual characterization, of risk alleles for FASD.

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.

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.

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.

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.

ClC-7 Regulates the Pattern and Early Development of Craniofacial Bone and Tooth

Human CLCN7 encodes voltage-gated chloride channel 7 (ClC-7); mutations of CLCN7 lead to osteopetrosis which is characterized by increased bone mass and impaired osteoclast function. In our previous clinical practice, we noticed that osteopetrosis patients with CLCN7 mutations had some special deformities in craniofacial morphology and tooth dysplasia. It is unclear whether these phenotypes are the typical features of CLCN7 involved osteopetrosis and whether ClC-7 could regulate the development of craniofacial bone and tooth in some signaling pathways. Methods: First, we collected 80 osteopetrosis cases from the literature and compared their craniofacial and dental phenotypes. Second, four osteopetrosis pedigrees with CLCN7 mutations were recruited from our clinic for gene testing and clinical analysis of their craniofacial and dental phenotypes. Third, we used a zebrafish model with clcn7 morpholino treatment to detect the effects of ClC-7 deficiency on the development of craniofacial and dental phenotypes. General observation, whole mount alcian blue and alizarin red staining, whole mount in situ hybridization, scanning electron microscope observation, lysoSensor staining, Q-PCR and western blotting were performed to observe the in vivo characteristics of craniofacial bone and tooth changes. Fourth, mouse marrow stromal cells were further primarily cultured to detect ClC-7 related mRNA and protein changes using siRNA, Q-PCR and western blotting. Results: Over 84% of osteopetrosis patients in the literature had some typical craniofacial and tooth phenotypes, including macrocephaly, frontal bossing, and changes in shape and proportions of facial skeleton, and these unique features are more severe and frequent in autosomal recessive osteopetrosis than in autosomal dominant osteopetrosis patients. Our four pedigrees with CLCN7 mutations confirmed the aforementioned clinical features. clcn7 knockdown in zebrafish reproduced the craniofacial cartilage defects and various dental malformations combined the decreased levels of col10a1, sp7, dlx2b, eve1, and cx43. Loss of clcn7 function resulted in lysosomal storage in the brain and jaw as well as downregulated cathepsin K (CTSK). The craniofacial phenotype severity also presented a dose-dependent relationship with the levels of ClC-7 and CTSK. ClC-7/CTSK further altered the balance of TGF-β/BMP signaling pathway, causing elevated TGF-β-like Smad2 signals and reduced BMP-like Smad1/5/8 signals in clcn7 morphants. SB431542 inhibitor of TGF-β pathway partially rescued the aforementioned craniofacial bone and tooth defects of clcn7 morphants. The ClC-7 involved CTSK/BMP and SMAD changes were also confirmed in mouse bone marrow stromal cells. Conclusion: These findings highlighted the vital role of clcn7 in zebrafish craniofacial bone and tooth development and mineralization, revealing novel insights for the causation of osteopetrosis with CLCN7 mutations. The mechanism chain of ClC-7/CTSK/ TGF-β/BMP/SMAD might explain the typical craniofacial bone and tooth changes in osteopetrosis as well as pycnodysostosis patients.

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.

Zebrafish as an Emerging Model for Osteoporosis: A Primary Testing Platform for Screening New Osteo-Active Compounds

Osteoporosis is metabolic bone disease caused by an altered balance between bone anabolism and catabolism. This dysregulated balance is responsible for fragile bones that fracture easily after minor falls. With an aging population, the incidence is rising and as yet pharmaceutical options to restore this imbalance is limited, especially stimulating osteoblast bone-building activity. Excitingly, output from large genetic studies on people with high bone mass (HBM) cases and genome wide association studies (GWAS) on the population, yielded new insights into pathways containing osteo-anabolic players that have potential for drug target development. However, a bottleneck in development of new treatments targeting these putative osteo-anabolic genes is the lack of animal models for rapid and affordable testing to generate functional data and that simultaneously can be used as a compound testing platform. Zebrafish, a small teleost fish, are increasingly used in functional genomics and drug screening assays which resulted in new treatments in the clinic for other diseases. In this review we outline the zebrafish as a powerful model for osteoporosis research to validate potential therapeutic candidates, describe the tools and assays that can be used to study bone homeostasis, and affordable (semi-)high-throughput compound testing.

An Adult Zebrafish Model of Fibrodysplasia Ossificans Progressiva

Fibrodysplasia ossificans progressiva (FOP) is a rare human skeletal disease caused by constitutively activating mutations in the gene ACVR1, which encodes a type I BMP/TGFβ family member receptor. FOP is characterized by progressive heterotopic ossification (HO) of fibrous tissues, including skeletal muscle, tendons, and ligaments, as well as malformation of the big toes, vertebral fusions, and osteochondromas. Surgical interventions in patients often result in enhanced HO, which can exacerbate rather than improve diagnostic outcomes. As a result of these difficulties, a variety of animal models are needed to study human FOP. Here we describe the methods for creating and characterizing zebrafish conditionally expressing Acvr1l^Q204D, the first adult zebrafish model for FOP.

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.

Endodermal pouch-expressed dmrt2b is important for pharyngeal cartilage formation

Pharyngeal pouches, a series of outpocketings derived from the foregut endoderm, are essential for craniofacial skeleton formation. However, the molecular mechanisms underlying endodermal pouch-regulated head cartilage development are not fully understood. In this study, we find that zebrafish dmrt2b, a gene encoding Doublesex- and Mab-3-related transcription factor, is specifically expressed in endodermal pouches and required for normal pharyngeal cartilage development. Loss of dmrt2b doesn't affect cranial neural crest (CNC) specification and migration, but leads to prechondrogenic condensation defects by reducing cxcl12b expression after CNC cell movement into the pharyngeal arches. Moreover, dmrt2b inactivation results in reduced proliferation and impaired differentiation of CNC cells. We also show that dmrt2b suppresses crossveinless 2 expression in endodermal pouches to maintain BMP/Smad signaling in the arches, thereby facilitating CNC cell proliferation and chondrogenic differentiation. This work provides insight into how transcription factors expressed in endodermal pouches regulate pharyngeal skeleton development through tissue-tissue interactions.

Maternal Bisphenol A exposure impaired endochondral ossification in craniofacial cartilage of rare minnow (Gobiocypris rarus) offspring

Bisphenol A (BPA), an endocrine disrupting compound, is present in the aquatic environment. BPA can mimic estrogen and cause adverse effects on development and reproduction in different organisms. As epigenetic modifications due to BPA exposure have been reported, the interest on the effects of this chemical has increased. To assess the potential effects of maternal BPA exposure on offspring bone development, adult Gobiocypris rarus (G. rarus) females were exposed to 15 μg L^-1 and 225 μg L^-1 BPA for 21 days. Eggs were collected after artificial spawning and fertilized with the fresh milt of non-exposed male fish. The offspring were raised in clean water and randomly selected for examination at different development stages. Our results showed that specific effects including poor quality of the embryos, increased malformation (bent spine and tail), and delayed craniofacial cartilage ossification of the larvae. Additionally, the transcripts of ossification related genes were significantly downregulated in offspring, and the lysyloxidase activity decreased. The present study demonstrated the maternal-mediated skeleton toxicity of BPA and its adverse effects on G. rarus.

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.

Evolution and development of the fish jaw skeleton

The evolution of the jaw represents a key innovation in driving the diversification of vertebrate body plans and behavior. The pharyngeal apparatus originated as gill bars separated by slits in chordate ancestors to vertebrates. Later, with the acquisition of neural crest, pharyngeal arches gave rise to branchial basket cartilages in jawless vertebrates (agnathans), and later bone and cartilage of the jaw, jaw support, and gills of jawed vertebrates (gnathostomes). Major events in the evolution of jaw structure from agnathans to gnathostomes include axial regionalization of pharyngeal elements and formation of a jaw joint. Hox genes specify the anterior-posterior identity of arches, and edn1, dlx, hand2, Jag1b-Notch2 signaling, and Nr2f factors specify dorsal-ventral identity. The formation of a jaw joint, an important step in the transition from an un-jointed pharynx in agnathans to a hinged jaw in gnathostomes involves interaction between nkx3.2, hand2, and barx1 factors. Major events in jaw patterning between fishes and reptiles include changes to elements of the second pharyngeal arch, including a loss of opercular and branchiostegal ray bones and transformation of the hyomandibula into the stapes. Further changes occurred between reptiles and mammals, including the transformation of the articular and quadrate elements of the jaw joint into the malleus and incus of the middle ear. Fossils of transitional jaw phenotypes can be analyzed from a developmental perspective, and there exists potential to use genetic manipulation techniques in extant taxa to test hypotheses about the evolution of jaw patterning in ancient vertebrates. This article is categorized under: Comparative Development and Evolution > Evolutionary Novelties Early Embryonic Development > Development to the Basic Body Plan Comparative Development and Evolution > Body Plan Evolution.

Injury of Adult Zebrafish Expressing Acvr1lQ204D Does Not Result in Heterotopic Ossification

Fibrodysplasia Ossificans Progressiva (FOP) is a rare, autosomal dominant genetic disorder in humans characterized by the gradual ossification of fibrous tissues, including skeletal muscle, tendons, and ligaments. In humans, mutations in the Type I BMP/TGFβ family member receptor gene, ACVR1, are associated with FOP. Zebrafish acvr1l, previously known as alk8, is the functional ortholog of human ACVR1. We previously created and characterized the first adult zebrafish model for FOP by generating animals harboring heat shock-inducible mCherry-tagged constitutively active Acvr1l (Q204D). Since injury is a known trigger for heterotopic ossification (HO) development in human FOP patients, in this study, we investigated several injury models in Acvr1l^Q204D-expressing zebrafish and the subsequent formation of HO. We performed studies of Activin A injection, cardiotoxin (CTX) injection, and caudal fin clip injury. We found that none of these methods resulted in HO formation at the site of injury. However, some of the cardiotoxin-injected and caudal fin-clipped animals did exhibit HO at distant sites, including the body cavity and along the spine. We describe these results in the context of new and exciting reports on FOP, and discuss future studies to better understand the etiology and progression of this disease.

Fox proteins are modular competency factors for facial cartilage and tooth specification

Facial form depends on the precise positioning of cartilage, bone, and tooth fields in the embryonic pharyngeal arches. How complex signaling information is integrated to specify these cell types remains a mystery. We find that modular expression of Forkhead domain transcription factors (Fox proteins) in the zebrafish face arises through integration of Hh, Fgf, Bmp, Edn1 and Jagged-Notch pathways. Whereas loss of C-class Fox proteins results in reduced upper facial cartilages, loss of F-class Fox proteins results in distal jaw truncations and absent midline cartilages and teeth. We show that Fox proteins are required for Sox9a to promote chondrogenic gene expression. Fox proteins are sufficient in neural crest-derived cells for cartilage development, and neural crest-specific misexpression of Fox proteins expands the cartilage domain but inhibits bone. These results support a modular role for Fox proteins in establishing the competency of progenitors to form cartilage and teeth in the face.

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.

Genome-wide analysis of facial skeletal regionalization in zebrafish

Patterning of the facial skeleton involves the precise deployment of thousands of genes in distinct regions of the pharyngeal arches. Despite the significance for craniofacial development, how genetic programs drive this regionalization remains incompletely understood. Here we use combinatorial labeling of zebrafish cranial neural crest-derived cells (CNCCs) to define global gene expression along the dorsoventral axis of the developing arches. Intersection of region-specific transcriptomes with expression changes in response to signaling perturbations demonstrates complex roles for Endothelin 1 (Edn1) signaling in the intermediate joint-forming region, yet a surprisingly minor role in ventralmost regions. Analysis of co-variance across multiple sequencing experiments further reveals clusters of co-regulated genes, with in situ hybridization confirming the domain-specific expression of novel genes. We then created loss-of-function alleles for 12 genes and uncovered antagonistic functions of two new Edn1 targets, follistatin a (fsta) and emx2, in regulating cartilaginous joints in the hyoid arch. Our unbiased discovery and functional analysis of genes with regional expression in zebrafish arch CNCCs reveals complex regulation by Edn1 and points to novel candidates for craniofacial disorders.

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.

A Zebrafish Model of Human Fibrodysplasia Ossificans Progressiva

Fibrodysplasia ossificans progressiva (FOP) is a rare, autosomal dominant genetic disorder in humans characterized by explosive inflammatory response to injury leading to gradual ossification within fibrous tissues, including skeletal muscle, tendons, and ligaments. A variety of animal models are needed to study and understand the etiology of human FOP. To address this need, here we present characterizations of the first adult zebrafish model for FOP. In humans, activating mutations in the Type I BMP/TGFβ family member receptor, ACVR1, are associated with FOP. Zebrafish acvr1l, previously known as alk8, is the functional ortholog of human ACVR1, and has been studied extensively in the developing zebrafish embryo, where it plays a role in early dorsoventral patterning. Constitutively active and dominant negative mutations in zebrafish acvr1l cause early lethal defects. Therefore, to study roles for activating acvr1l mutations in adult zebrafish, we created transgenic animals expressing mCherry-tagged, heat-shock-inducible constitutively active Acvr1l, Acvr1l^Q204D, to investigate phenotypes in juvenile and adult zebrafish. Our studies showed that adult zebrafish expressing heat-shock-induced Acvr1l^Q204D develop a number of human FOP-like phenotypes, including heterotopic ossification lesions, spinal lordosis, vertebral fusions, and malformed pelvic fins. Together, these results suggest that transgenic zebrafish expressing heat-shock-inducible Acvr1l^Q204D can serve as a model for human FOP.

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.

Wdr68 Mediates Dorsal and Ventral Patterning Events for Craniofacial Development

Birth defects are among the leading causes of infant mortality and contribute substantially to illness and long-term disability. Defects in Bone Morphogenetic Protein (BMP) signaling are associated with cleft lip/palate. Many craniofacial syndromes are caused by defects in signaling pathways that pattern the cranial neural crest cells (CNCCs) along the dorsal-ventral axis. For example, auriculocondylar syndrome is caused by impaired Endothelin-1 (Edn1) signaling, and Alagille syndrome is caused by defects in Jagged-Notch signaling. The BMP, Edn1, and Jag1b pathways intersect because BMP signaling is required for ventral edn1 expression that, in turn, restricts jag1b to dorsal CNCC territory. In zebrafish, the scaffolding protein Wdr68 is required for edn1 expression and subsequent formation of the ventral Meckel’s cartilage as well as the dorsal Palatoquadrate. Here we report that wdr68 activity is required between the 17-somites and prim-5 stages, that edn1 functions downstream of wdr68, and that wdr68 activity restricts jag1b, hey1, and grem2 expression from ventral CNCC territory. Expression of dlx1a and dlx2a was also severely reduced in anterior dorsal and ventral 1^st arch CNCC territory in wdr68 mutants. We also found that the BMP agonist isoliquiritigenin (ISL) can partially rescue lower jaw formation and edn1 expression in wdr68 mutants. However, we found no significant defects in BMP reporter induction or pSmad1/5 accumulation in wdr68 mutant cells or zebrafish. The Transforming Growth Factor Beta (TGF-β) signaling pathway is also known to be important for craniofacial development and can interfere with BMP signaling. Here we further report that TGF-β interference with BMP signaling was greater in wdr68 mutant cells relative to control cells. To determine whether interference might also act in vivo, we treated wdr68 mutant zebrafish embryos with the TGF-β signaling inhibitor SB431542 and found partial rescue of edn1 expression and craniofacial development. While ISL treatment failed, SB431542 partially rescued dlx2a expression in wdr68 mutants. Together these findings reveal an indirect role for Wdr68 in the BMP-Edn1-Jag1b signaling hierarchy and dorso-anterior expression of dlx1a/2a.

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.

Noggin is required for first pharyngeal arch differentiation in the frog Xenopus tropicalis

The dorsal ventral axis of vertebrates requires high BMP activity for ventral development and inhibition of BMP activity for dorsal development. Presumptive dorsal regions of the embryo are protected from the ventralizing activity of BMPs by the secretion of BMP antagonists from the mesoderm. Noggin, one such antagonist, binds BMP ligands and prevents them from binding their receptors, however, a unique role for Noggin in amphibian development has remained unclear. Previously, we used zinc-finger nucleases to mutagenize the noggin locus in Xenopus tropicalis. Here, we report on the phenotype of noggin mutant frogs as a result of breeding null mutations to homozygosity. Early homozygous noggin mutant embryos are indistinguishable from wildtype siblings, with normal neural induction and neural tube closure. However, in late tadpole stages mutants present severe ventral craniofacial defects, notably a fusion of Meckel's cartilage to the palatoquadrate cartilage. Consistent with a noggin loss-of-function, mutants show expansions of BMP target gene expression and the mutant phenotype can be rescued with transient BMP inhibition. These results demonstrate that in amphibians, Noggin is dispensable for early embryonic patterning but is critical for cranial skeletogenesis.