hand2 and Dlx genes specify dorsal, intermediate and ventral domains within zebrafish pharyngeal arches

Pax9 drives development of the upper jaw but not teeth in zebrafish

Loss of dentition has occurred repeatedly throughout vertebrate evolution. Cyprinid fish, including zebrafish, form teeth only deep within the pharynx, not on their oral jaws. However, zebrafish still robustly express transcription factors associated with mammalian tooth development in the neural crest-derived mesenchyme surrounding the mouth. We investigated whether this expression is vestigial or whether these factors contribute to the formation of non-tooth mesenchymal structures in the oral region, using Pax9 as a test case. Zebrafish homozygous for two different pax9 mutant alleles develop the normal complement of pharyngeal teeth but fail to form the premaxilla bone, most of the maxilla, and nasal and maxillary barbels. Lack of most of the upper jaw complex does not preclude effective feeding in the laboratory environment. We observe a significant deficit of sp7:EGFP + osteoblasts and adjacent alx4a:DsRed+ condensing mesenchyme around the maxilla, and no accumulation of either in the premaxillary domain. Loss of pax9 may prevent osteoprogenitors from maintaining the state of condensation required for full osteogenic differentiation. We conclude that Pax9 is not unequivocally required for all vertebrate tooth development but instead may be involved in the development of a variety of organs forming through mesenchymal condensation around the mouth.

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.

Assessing candidate DLX-regulated genes in the first pharyngeal arch of chick embryos

BACKGROUND

Insights into the development and evolution of asymmetrical jaws will require an understanding of the gene regulatory networks that underpin the differential morphogenesis of the maxillary and mandibular domains of the first pharyngeal arch in a variety of gnathostomes. While a robust relationship has been demonstrated between jaw patterning and the Endothelin-Dlx gene axis, much less is known of the next level of genes in the jaw patterning hierarchy.

RESULTS

Several genes, whose expression depends on Dlx5 and/or Dlx6, have been identified in mice. Here, we examined the expression patterns of the chick orthologues of some of those genes, namely GSC, PITX1, HAND2, and GBX2, and tested their dependence on endothelin signaling to assess whether there is a conserved regulatory relationship between those genes in the chick embryo. To further validate these genes as direct DLX targets, we identified conserved non-coding sequences containing candidate DLX binding motifs and demonstrated DLX-responsiveness in vitro.

CONCLUSIONS

The evidence presented in this study combines to support the hypothesis that these four genes are direct targets of DLX transcription factors in the lower jaw-forming tissue.

Combining genetic and single-cell expression data reveals cell types and novel candidate genes for orofacial clefting

Non-syndromic cleft lip with/without cleft palate (nsCL/P) is one of the most common birth defects and has a multifactorial etiology. To date, over 45 loci harboring common risk variants have been identified. However, the effector genes at these loci, and the cell types that are affected by risk alleles, remain largely unknown. To address this, we combined genetic data from an nsCL/P genome-wide association study (GWAS) with single-cell RNA sequencing data obtained from the heads of unaffected human embryos. Using the recently developed single-cell disease relevance score (scDRS) approach, we identified two major cell types involved in nsCL/P development, namely the epithelium and the HAND2+ pharyngeal arches (PA). Combining scDRS with co-expression networks and differential gene expression analysis, we prioritized nsCL/P candidate genes, some of which were additionally supported by GWAS data (e.g., CTNND1, PRTG, RPL35A, RAB11FIP1, KRT19). Our results suggest that specific epithelial and PA sub-cell types are involved in nsCL/P development, and harbor a substantial fraction of the genetic risk for nsCL/P.

Structure and function of the larval teleost fish gill

The fish gill is a multifunctional organ that is important in multiple physiological processes such as gas transfer, ionoregulation, and chemoreception. This characteristic organ of fishes has received much attention, yet an often-overlooked point is that larval fishes in most cases do not have a fully developed gill, and thus larval gills do not function identically as adult gills. In addition, large changes associated with gas exchange and ionoregulation happen in gills during the larval phase, leading to the oxygen and ionoregulatory hypotheses examining the environmental constraint that resulted in the evolution of gills. This review thus focuses exclusively on the larval fish gill of teleosts, summarizing the development of teleost larval fish gills and its function in gas transfer, ionoregulation, and chemoreception, and comparing and contrasting it to adult gills where applicable, while providing some insight into the oxygen vs ionoregulatory hypotheses debate.

Fast-twitch myofibrils grow in proportion to Mylpf dosage in the zebrafish embryo

Muscle cells become stronger by expanding myofibrils, the chains of sarcomeres that produce contraction. Here we investigate how Mylpf (Myosin Light Chain Phosphorylatable Fast) abundance impacts myofibril assembly in fast-twitch muscle. The two zebrafish Mylpf genes (mylpfa and mylpfb) are exclusively expressed in fast-twitch muscle. We show that these cells initially produce six times more mylpfa mRNA and protein than mylpfb. The combined Mylpf protein dosage is necessary for and proportionate to fast-twitch myofibril growth in the embryo. Fast-twitch myofibrils are severely reduced in the mylpfa -/- mutant, leading to loss of high-speed movement; however, by persistent slow movement this mutant swims as far through time as its wild-type sibling. Although the mylpfb -/- mutant has normal myofibrils, myofibril formation fails entirely in the mylpfa -/- ;mylpfb -/- double mutant, indicating that the two genes are collectively essential to myofibril formation. Fast-twitch myofibril width is restored in the mylpfa -/- mutant by transgenic expression of mylpfa-GFP, mylpfb-GFP, and by human MYLPF-GFP to a degree corresponding linearly with GFP brightness. This correlate is inverted by expression of MYLPF alleles that cause Distal Arthrogryposis, which reduce myofibril size in proportion to protein abundance. These effects indicate that Mylpf dosage controls myofibril growth, impacting embryonic development and lifelong health.

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.

DLX gene expression in the developing chick pharyngeal arches and relationship to endothelin signaling and avian jaw patterning

BACKGROUND

A hinged jaw that articulates with the skull base is a striking feature of the vertebrate head and has been greatly modified between, and within, vertebrate classes. Genes belonging to the DLX homeobox family are conserved mediators of local signaling pathways that distinguish the dorsal and ventral aspects of the first pharyngeal arch. Specifically, a subset of DLX genes are expressed in the cranial neural crest-derived mandibular ectomesenchyme in response to ventral endothelin signaling, an important step that confers the first arch with maxillary and mandibular identities. Downstream targets of DLX genes then execute the morphogenetic processes that lead to functional jaws. Identifying lineage-specific variations in DLX gene expression and the regulatory networks downstream of DLX action is necessary to understand how different kinds of jaws evolved.

RESULTS

Here, we describe and compare the expression of all six DLX genes in the chick pharyngeal arches, focusing on the period of active patterning in the first arch. Disruption of endothelin signaling results in the down-regulation of ventral-specific DLX genes and confirms their functional role in avian jaw patterning.

CONCLUSIONS

This expression resource will be important for comparative embryology and for identifying synexpression groups of DLX-regulated genes in the chick.

Contribution of Signal Transducer and Activator of Transcription 3 (STAT3) to Bone Development and Repair

Signal transducer and activator of transcription 3 (STAT3) is a transcription factor activated canonically by numerous cytokines and other factors, with significant roles in immunity, immune diseases, and cancer. It has also been implicated in several human skeletal disorders, with loss-of-function (LOF) mutations associated with aberrant skeletal development. To gain further insights, two zebrafish STAT3 lines were investigated: a complete LOF knockout (KO) mutant and a partial LOF mutant with the transactivation domain truncated (ΔTAD). Consistent with other studies, the KO mutants were smaller, with reduced length in early embryos exacerbated by a decreased growth rate from 5 days postfertilization (dpf). They displayed skeletal deformities that approached 80% incidence by 30 dpf, with a significant reduction in early bone but not cartilage formation. Further analysis additionally identified considerable abrogation of caudal fin regeneration, concomitant with a paucity of infiltrating macrophages and neutrophils, which may be responsible for this. Most of these phenotypes were also observed in the ΔTAD mutants, indicating that loss of canonical STAT3 signaling was the likely cause. However, the impacts on early bone formation and regeneration were muted in the ΔTAD mutant, suggesting the potential involvement of noncanonical functions in these processes.

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.

Craniofacial developmental biology in the single-cell era

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

T, Mishoe Hernandez L, DeLaurier A. Exposure to valproic acid (VPA) reproduces hdac1 loss of function phenotypes in zebrafish

Histone deacetylases are enzymes that remove acetyl groups from histone tails and are understood to act as repressors of transcriptional activity. Hdac1 has been previously shown to function in eye, pectoral fin, heart, liver, and pharyngeal skeletal development. We show that high doses of Valproic Acid (VPA) reproduce the hdac1 phenotype. We identify tbx5 genes as potential targets of Hdac1 in eye, pectoral fin, and heart development. Using timed exposures, we show that skeletal structures in the pharyngeal arches are impacted by VPA between 24-36 hours post-fertilization, indicating a role for Hdac1 during post-migration patterning, differentiation, or proliferation of cranial neural crest cells.

Loss of dlx5a/dlx6a Locus Alters Non-Canonical Wnt Signaling and Meckel's Cartilage Morphology

The dlx genes encode transcription factors that establish a proximal-distal polarity within neural crest cells to bestow a regional identity during craniofacial development. The expression regions of dlx paralogs are overlapping yet distinct within the zebrafish pharyngeal arches and may also be involved in progressive morphologic changes and organization of chondrocytes of the face. However, how each dlx paralog of dlx1a, dlx2a, dlx5a and dlx6a affects craniofacial development is still largely unknown. We report here that the average lengths of the Meckel's, palatoquadrate and ceratohyal cartilages in different dlx mutants were altered. Mutants for dlx5a-/- and dlx5i6-/-, where the entire dlx5a/dlx6a locus was deleted, have the shortest lengths for all three structures at 5 days post fertilization (dpf). This phenotype was also observed in 14 dpf larvae. Loss of dlx5i6 also resulted in increased proliferation of neural crest cells and expression of chondrogenic markers. Additionally, altered expression and function of non-canonical Wnt signaling were observed in these mutants suggesting a novel interaction between dlx5i6 locus and non-canonical Wnt pathway regulating ventral cartilage morphogenesis.

Cartilage diversification and modularity drove the evolution of the ancestral vertebrate head skeleton

The vertebrate head skeleton has evolved a myriad of forms since their divergence from invertebrate chordates. The connection between novel gene expression and cell types is therefore of importance in this process. The transformation of the jawed vertebrate (gnathostome) head skeleton from oral cirri to jointed jaw elements required a diversity of cartilages as well as changes in the patterning of these tissues. Although lampreys are a sister clade to gnathostomes, they display skeletal diversity with distinct gene expression and histologies, a useful model for addressing joint evolution. Specifically, the lamprey tissue known as mucocartilage has noted similarities with the jointed elements of the mandibular arch in jawed vertebrates. We thus asked whether the cells in lamprey mucocartilage and gnathostome joint tissue could be considered homologous. To do this, we characterized new genes that are involved in gnathostome joint formation and characterized the histochemical properties of lamprey skeletal types. We find that most of these genes are minimally found in mucocartilage and are likely later innovations, but we do identify new activity for gdf5/6/7b in both hyaline and mucocartilage, supporting its role as a chondrogenic regulator. Contrary to previous works, our histological assays do not find any perichondrial fibroblasts surrounding mucocartilage, suggesting that mucocartilage is non-skeletogenic tissue that is partially chondrified. Interestingly, we also identify new histochemical features of the lamprey otic capsule that diverge from normal hyaline. Paired with our new insights into lamprey mucocartilage, we propose a broader framework for skeletal evolution in which an ancestral soxD/E and gdf5/6/7 network directs mesenchyme along a spectrum of cartilage-like features.

PRDM1 DNA-binding zinc finger domain is required for normal limb development and is disrupted in split hand/foot malformation

Split hand/foot malformation (SHFM) is a rare limb abnormality with clefting of the fingers and/or toes. For many individuals, the genetic etiology is unknown. Through whole-exome and targeted sequencing, we detected three novel variants in a gene encoding a transcription factor, PRDM1, that arose de novo in families with SHFM or segregated with the phenotype. PRDM1 is required for limb development; however, its role is not well understood and it is unclear how the PRDM1 variants affect protein function. Using transient and stable overexpression rescue experiments in zebrafish, we show that the variants disrupt the proline/serine-rich and DNA-binding zinc finger domains, resulting in a dominant-negative effect. Through gene expression assays, RNA sequencing, and CUT&RUN in isolated pectoral fin cells, we demonstrate that Prdm1a directly binds to and regulates genes required for fin induction, outgrowth and anterior/posterior patterning, such as fgfr1a, dlx5a, dlx6a and smo. Taken together, these results improve our understanding of the role of PRDM1 in the limb gene regulatory network and identified novel PRDM1 variants that link to SHFM in humans.

Genetically engineered zebrafish as models of skeletal development and regeneration

Zebrafish (Danio rerio) are aquatic vertebrates with significant homology to their terrestrial counterparts. While zebrafish have a centuries-long track record in developmental and regenerative biology, their utility has grown exponentially with the onset of modern genetics. This is exemplified in studies focused on skeletal development and repair. Herein, the numerous contributions of zebrafish to our understanding of the basic science of cartilage, bone, tendon/ligament, and other skeletal tissues are described, with a particular focus on applications to development and regeneration. We summarize the genetic strengths that have made the zebrafish a powerful model to understand skeletal biology. We also highlight the large body of existing tools and techniques available to understand skeletal development and repair in the zebrafish and introduce emerging methods that will aid in novel discoveries in skeletal biology. Finally, we review the unique contributions of zebrafish to our understanding of regeneration and highlight diverse routes of repair in different contexts of injury. We conclude that zebrafish will continue to fill a niche of increasing breadth and depth in the study of basic cellular mechanisms of skeletal biology.

Dynamic regulation of gene expression and morphogenesis in the zebrafish embryo test after exposure to all-trans retinoic acid

The zebrafish embryotoxicity test (ZET) is widely used in developmental toxicology. The analysis of gene expression regulation in ZET after chemical exposure provides mechanistic information about the effects of chemicals on morphogenesis in the test. The gene expression response magnitude has been shown to change with exposure duration. The objective of this work is to study the effect of the exposure duration on the magnitude of gene expression changes in the all-trans retinoic acid (ATRA) signaling pathway in the ZET. Retinoic acid regulation is a key driver of morphogenesis and is therefore employed here as an indicator for the regulation of developmental genes. A teratogenic concentration of 7.5 nM of ATRA was given at 3 hrs post fertilization (hpf) for a range of exposure durations until 120 hrs of development. The expression of a selection of genes related to ATRA signaling and downstream developmental genes was determined. The highest magnitudes of gene expression regulation were observed after 2-24 hrs exposure with an optimal response after 4 hrs. Longer exposures showed a decrease in the gene expression response, although continued exposure to 120 hpf caused malformations and lethality. This study shows that assessment of gene expression regulation at early time points after the onset of exposure in the ZET may be optimal for the prediction of developmental toxicity. We believe these results could help optimize sensitivity in future studies with ZET.

A novel cis-regulatory element drives early expression of Nkx3.2 in the gnathostome primary jaw joint

The acquisition of movable jaws was a major event during vertebrate evolution. The role of NK3 homeobox 2 (Nkx3.2) transcription factor in patterning the primary jaw joint of gnathostomes (jawed vertebrates) is well known, however knowledge about its regulatory mechanism is lacking. In this study, we report a proximal enhancer element of Nkx3.2 that is deeply conserved in most gnathostomes but undetectable in the jawless hagfish and lamprey. This enhancer is active in the developing jaw joint region of the zebrafish Danio rerio, and was thus designated as jaw joint regulatory sequence 1 (JRS1). We further show that JRS1 enhancer sequences from a range of gnathostome species, including a chondrichthyan and mammals, have the same activity in the jaw joint as the native zebrafish enhancer, indicating a high degree of functional conservation despite the divergence of cartilaginous and bony fish lineages or the transition of the primary jaw joint into the middle ear of mammals. Finally, we show that deletion of JRS1 from the zebrafish genome using CRISPR/Cas9 results in a significant reduction of early gene expression of nkx3.2 and leads to a transient jaw joint deformation and partial fusion. Emergence of this Nkx3.2 enhancer in early gnathostomes may have contributed to the origin and shaping of the articulating surfaces of vertebrate jaws.

Coordinated patterning of zebrafish caudal fin symmetry by a central and two peripheral organizers

BACKGROUND

Caudal fin symmetry characterizes teleosts and likely contributes to their evolutionary success. However, the coordinated development and patterning of skeletal elements establishing external symmetry remains incompletely understood. We explore the spatiotemporal emergence of caudal skeletal elements in zebrafish to consider evolutionary and developmental origins of caudal fin symmetry.

RESULTS

Transgenic reporters and skeletal staining reveal that the hypural diastema-defining gap between hypurals 2 and 3 forms early and separates progenitors of two plates of connective tissue. Two sets of central principal rays (CPRs) synchronously, sequentially, and symmetrically emerge around the diastema. The two dorsal- and ventral-most rays (peripheral principal rays, PPRs) arise independently and earlier than adjacent CPRs. Muscle and tendon markers reveal that different muscles attach to CPR and PPR sets.

CONCLUSIONS

We propose that caudal fin symmetry originates from a central organizer that establishes the hypural diastema and bidirectionally patterns surrounding tissue into two plates of connective tissue and two mirrored sets of CPRs. Further, two peripheral organizers unidirectionally specify PPRs, forming a symmetric "composite" fin derived from three fields. Distinct CPR and PPR ontogenies may represent developmental modules conferring ray identities, muscle connections, and biomechanical properties. Our model contextualizes mechanistic studies of teleost fin morphological variation.

Photo-enhanced toxicity of crude oil on early developmental stages of Atlantic cod (Gadus morhua)

Photo-enhanced toxicity of crude oil is produced by exposure to ultraviolet (UV) radiation. Atlantic cod (Gadus morhua) embryos were exposed to crude oil with and without UV radiation (290-400 nm) from 3 days post fertilization (dpf) until 6 dpf. Embryos from the co-exposure experiment were continually exposed to UV radiation until hatching at 11 dpf. Differences in body burden levels and cyp1a expression in cod embryos were observed between the exposure regimes. High doses of crude oil produced increased mortality in cod co-exposed embryos, as well as craniofacial malformations and heart deformities in larvae from both experiments. A higher number of differentially expressed genes (DEGs) and pathways were revealed in the co-exposure experiment, indicating a photo-enhanced effect of crude oil toxicity. Our results provide mechanistic insights into crude oil and photo-enhanced crude oil toxicity, suggesting that UV radiation increases the toxicity of crude oil in early life stages of Atlantic cod.

Regeneration of Jaw Joint Cartilage in Adult Zebrafish

The poor intrinsic repair capacity of mammalian joint cartilage likely contributes to the high incidence of arthritis worldwide. Adult zebrafish can regenerate many structures that show limited or no healing capacity in mammals, including the jawbone. To test whether zebrafish can also regenerate damaged joints, we developed a surgical injury model in which the zebrafish jaw joint is destabilized via transection of the major jaw joint ligament, the interopercular-mandibular (IOM). Unilateral transection of the IOM ligament in 1-year-old fish resulted in an initial reduction of jaw joint cartilage by 14 days, with full regeneration of joint cartilage by 28 days. Joint cartilage regeneration involves the re-entry of articular chondrocytes into the cell cycle and the upregulated expression of sox10, a marker of developing chondrocytes in the embryo that becomes restricted to a subset of joint chondrocytes in adults. Genetic ablation of these sox10-expressing chondrocytes shows that they are essential for joint cartilage regeneration. To uncover the potential source of new chondrocytes during joint regeneration, we performed single-cell RNA sequencing of the uninjured adult jaw joint and identified multiple skeletal, connective tissue, and fibroblast subtypes. In particular, we uncovered a joint-specific periosteal population expressing coch and grem1a, with the jaw joint chondrocytes marked by grem1a expression during regeneration. Our findings demonstrate the capacity of zebrafish to regenerate adult joint cartilage and identify candidate cell types that can be tested for their roles in regenerative response.

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.

A comprehensive series of Irx cluster mutants reveals diverse roles in facial cartilage development

Proper function of the vertebrate skeleton requires the development of distinct articulating embryonic cartilages. Irx transcription factors are arranged in co-regulated clusters that are expressed in the developing skeletons of the face and appendages. IrxB cluster genes are required for the separation of toes in mice and formation of the hyoid joint in zebrafish, yet whether Irx genes have broader roles in skeletal development remains unclear. Here, we perform a comprehensive loss-of-function analysis of all 11 Irx genes in zebrafish. We uncover conserved requirements for IrxB genes in formation of the fish and mouse scapula. In the face, we find a requirement for IrxAb genes and irx7 in formation of anterior neural crest precursors of the jaw, and for IrxBa genes in formation of endodermal pouches and gill cartilages. We also observe extensive joint loss and cartilage fusions in animals with combinatorial losses of Irx clusters, with in vivo imaging revealing that at least some of these fusions arise through inappropriate chondrogenesis. Our analysis reveals diverse roles for Irx genes in the formation and later segmentation of the facial skeleton.

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.

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.

Transgene-mediated skeletal phenotypic variation in zebrafish

When considering relationships between genotype and phenotype we frequently ignore the fact that the genome of a typical animal, notably including that of a fish and a human, harbours a huge amount of foreign DNA. Such DNA, in the form of transposable elements, can affect genome function in a major way, and transgene biology needs to be included in our understanding of the genome. Here we examine an unexpected phenotypic effect of the chromosomally integrated transgene fli1a-F-hsp70l:Gal4VP16 that serves as a model for transgene function generally. We examine larval fras1 mutant zebrafish (Danio rerio). Gal4VP16 is a potent transcriptional activator that is already well known for toxicity and mediating unusual transcriptional effects. In the presence of the transgene, phenotypes in the neural crest-derived craniofacial skeleton, notably fusions and shape changes associated with loss of function fras1 mutations, are made more severe, as we quantify by scoring phenotypic penetrance, the fraction of mutants expressing the trait. A very interesting feature is that the enhancements are highly specific for fras1 mutant phenotypes, occurring in the apparent absence of more widespread changes. Except for the features due to the fras1 mutation, the transgene-bearing larvae appear generally healthy and to be developing normally. The transgene behaves as a genetic partial dominant: a single copy is sufficient for the enhancements, yet, for some traits, two copies may exert a stronger effect. We made new strains bearing independent insertions of the fli1a-F-hsp70l:Gal4VP16 transgene in new locations in the genome, and observed increased severities of the same phenotypes as observed for the original insertion. This finding suggests that sequences within the transgene, for example Gal4VP16, are responsible for the enhancements, rather than the effect on neighbouring host sequences (such as an insertional mutation). The specificity and biological action underlying the traits are subjects of considerable interest for further investigation, as we discuss. Our findings show that work with transgenes needs to be undertaken with caution and attention to detail.

The alx3 gene shapes the zebrafish neurocranium by regulating frontonasal neural crest cell differentiation timing

During craniofacial development, different populations of cartilage- and bone-forming cells develop in precise locations in the head. Most of these cells are derived from pluripotent cranial neural crest cells and differentiate with distinct developmental timing and cellular morphologies. The mechanisms that divide neural crest cells into discrete populations are not fully understood. Here, we use single-cell RNA sequencing to transcriptomically define different populations of cranial neural crest cells. We discovered that the gene family encoding the Alx transcription factors is enriched in the frontonasal population of neural crest cells. Genetic mutant analyses indicate that alx3 functions to regulate the distinct differentiation timing and cellular morphologies among frontonasal neural crest cell subpopulations. This study furthers our understanding of how genes controlling developmental timing shape craniofacial skeletal elements.

Zebrafish model for spondylo-megaepiphyseal-metaphyseal dysplasia reveals post-embryonic roles of Nkx3.2 in the skeleton

The regulated expansion of chondrocytes within growth plates and joints ensures proper skeletal development through adulthood. Mutations in the transcription factor NKX3.2 underlie spondylo-megaepiphyseal-metaphyseal dysplasia (SMMD), which is characterized by skeletal defects including scoliosis, large epiphyses, wide growth plates and supernumerary distal limb joints. Whereas nkx3.2 knockdown zebrafish and mouse Nkx3.2 mutants display embryonic lethal jaw joint fusions and skeletal reductions, respectively, they lack the skeletal overgrowth seen in SMMD patients. Here, we report adult viable nkx3.2 mutant zebrafish displaying cartilage overgrowth in place of a missing jaw joint, as well as severe dysmorphologies of the facial skeleton, skullcap and spine. In contrast, cartilage overgrowth and scoliosis are absent in rare viable nkx3.2 knockdown animals that lack jaw joints, supporting post-embryonic roles for Nkx3.2. Single-cell RNA-sequencing and in vivo validation reveal increased proliferation and upregulation of stress-induced pathways, including prostaglandin synthases, in mutant chondrocytes. By generating a zebrafish model for the skeletal overgrowth defects of SMMD, we reveal post-embryonic roles for Nkx3.2 in dampening proliferation and buffering the stress response in joint-associated chondrocytes.

The landscape of regulatory genes in brain-wide neuronal phenotypes of a vertebrate brain

Multidimensional landscapes of regulatory genes in neuronal phenotypes at whole-brain levels in the vertebrate remain elusive. We generated single-cell transcriptomes of ~67,000 region- and neurotransmitter/neuromodulator-identifiable cells from larval zebrafish brains. Hierarchical clustering based on effector gene profiles ('terminal features') distinguished major brain cell types. Sister clusters at hierarchical termini displayed similar terminal features. It was further verified by a population-level statistical method. Intriguingly, glutamatergic/GABAergic sister clusters mostly expressed distinct transcription factor (TF) profiles ('convergent pattern'), whereas neuromodulator-type sister clusters predominantly expressed the same TF profiles ('matched pattern'). Interestingly, glutamatergic/GABAergic clusters with similar TF profiles could also display different terminal features ('divergent pattern'). It led us to identify a library of RNA-binding proteins that differentially marked divergent pair clusters, suggesting the post-transcriptional regulation of neuron diversification. Thus, our findings reveal multidimensional landscapes of transcriptional and post-transcriptional regulators in whole-brain neuronal phenotypes in the zebrafish brain.

Genetic Factors in Nonsyndromic Orofacial Clefts

Orofacial clefts (OFCs) are the most common congenital birth defects in humans and immediately recognized at birth. The etiology remains complex and poorly understood and seems to result from multiple genetic and environmental factors along with gene-environment interactions. It can be classified into syndromic (30%) and nonsyndromic (70%) clefts. Nonsyndromic OFCs include clefts without any additional physical or cognitive deficits. Recently, various genetic approaches, such as genome-wide association studies (GWAS), candidate gene association studies, and linkage analysis, have identified multiple genes involved in the etiology of OFCs. This article provides an insight into the multiple genes involved in the etiology of OFCs. Identification of specific genetic causes of clefts helps in a better understanding of the molecular pathogenesis of OFC. In the near future, it helps to provide a more accurate diagnosis, genetic counseling, personalized medicine for better clinical care, and prevention of OFCs.

Spatiotemporal cellular movement and fate decisions during first pharyngeal arch morphogenesis

Cranial neural crest (CNC) cells contribute to different cell types during embryonic development. It is unknown whether postmigratory CNC cells undergo dynamic cellular movement and how the process of cell fate decision occurs within the first pharyngeal arch (FPA). Our investigations demonstrate notable heterogeneity within the CNC cells, refine the patterning domains, and identify progenitor cells within the FPA. These progenitor cells undergo fate bifurcation that separates them into common progenitors and mesenchymal cells, which are characterized by Cdk1 and Spry2/Notch2 expression, respectively. The common progenitors undergo further bifurcations to restrict them into osteogenic/odontogenic and chondrogenic/fibroblast lineages. Disruption of a patterning domain leads to specific mandible and tooth defects, validating the binary cell fate restriction process. Different from the compartment model of mandibular morphogenesis, our data redefine heterogeneous cellular domains within the FPA, reveal dynamic cellular movement in time, and describe a sequential series of binary cell fate decision-making process.

Evolution of vertebrate gill covers via shifts in an ancient Pou3f3 enhancer

Whereas the gill chambers of jawless vertebrates open directly into the environment, jawed vertebrates evolved skeletal appendages that drive oxygenated water unidirectionally over the gills. A major anatomical difference between the two jawed vertebrate lineages is the presence of a single large gill cover in bony fishes versus separate covers for each gill chamber in cartilaginous fishes. Here, we find that these divergent patterns correlate with the pharyngeal arch expression of Pou3f3 orthologs. We identify a deeply conserved Pou3f3 arch enhancer present in humans through sharks but undetectable in jawless fish. Minor differences between the bony and cartilaginous fish enhancers account for their restricted versus pan-arch expression patterns. In zebrafish, mutation of Pou3f3 or the conserved enhancer disrupts gill cover formation, whereas ectopic pan-arch Pou3f3b expression generates ectopic skeletal elements resembling the multimeric covers of cartilaginous fishes. Emergence of this Pou3f3 arch enhancer >430 Mya and subsequent modifications may thus have contributed to the acquisition and diversification of gill covers and respiratory strategies during gnathostome evolution.

Development of the anterior lateral line system through local tissue-tissue interactions in the zebrafish head

BACKGROUND

The distribution of sensory organs is important for detecting environmental signals efficiently. The mechanosensory receptors of the lateral line system, neuromasts, are stereotypically distributed over the head and body surface of fish, although how neuromasts arise in these predetermined positions during development remains unclear.

RESULTS

We investigated the development of the anterior lateral line (ALL) system in zebrafish head. The ALL neuromasts formed in the predetermined positions through proliferation and differentiation of (a) nonmigratory lateral line primordia, (b) migratory primordia, (c) interneuromast cells connecting preexisting neuromasts, and (d) budding primordia. We demonstrated that R-spondin2 (Rspo2), an activator of Wnt/β-catenin signaling, is required for the development of a particular set of neuromasts associated with hyomandibular cartilage. Further genetic analyses suggested that Rspo2, which emanates from the hyoid mesenchyme, acts on the adjacent neuromast progenitor cells to stimulate their proliferation through activating Wnt/β-catenin signaling.

CONCLUSION

This study has revealed novel mechanisms for neuromast positioning through local tissue-tissue interactions, providing insights into the development and evolution of the vertebrate head.

Mutant dlx3b disturbs normal tooth mineralization and bone formation in zebrafish

Background

Tricho-dento-osseous (TDO) syndrome is an autosomal dominant disorder characterized by anomalies in hair, teeth and bone (OMIM190320). Various mutations of Distal-Less 3 (DLX3) gene are found to be responsible for human TDO. The aim of this study was to investigate effects of DLX3 on tooth and bone development using a zebrafish model.

Methods

The dlx3b mutant zebrafish lines were established using the gene targeting tool transcription activator-like effector nuclease (TALEN). Micro-computed tomography was used to render the three-dimensional skeletal structures of mutant fishes. The pharyngeal bone along with connected teeth was isolated and stained by Alizarine Red S, then observed under stereomicroscope. Scanning electron microscopy (SEM) and energy dispersive spectrometer (EDS) were used to examine the tooth surface morphology and mineral composition. Quantitative real-time PCR was used to analyze gene expression.

Results

A moderate curvature of the spine toward the dorsal side was found at the early larval stages, appearing in 86 out of 100 larvae in dlx3b^-/- group as compared to 3 out of 99 in the dlx3b^+/+ group. At the adult stage, three of the thirty dlx3b^-/- homozygotes exhibited prominent abnormal curvature in the spine. SEM revealed morphological surface changes in pharyngeal teeth enameloid, accompanied by a decrease in the mineral content detected by EDS. Furthermore, specific secretory calcium-binding phosphoprotein (SCPP) genes, including odam, scpp9, spp1, scpp1, and scpp5 were significantly downregulated in dlx3b mutants.

Conclusion

The findings of this study suggest that dlx3b is critical for enamel mineralization and bone formation in zebrafish. Moreover, the discovery of the downregulation of SCPP genes in dlx3b mutants sheds new light on the molecular mechanisms underlying TDO syndrome.

Molecular mechanisms underlying nuchal hump formation in dolphin cichlid, Cyrtocara moorii

East African cichlid fishes represent a model to tackle adaptive changes and their connection to rapid speciation and ecological distinction. In comparison to bony craniofacial tissues, adaptive morphogenesis of soft tissues has been rarely addressed, particularly at the molecular level. The nuchal hump in cichlids fishes is one such soft-tissue and exaggerated trait that is hypothesized to play an innovative role in the adaptive radiation of cichlids fishes. It has also evolved in parallel across lakes in East Africa and Central America. Using gene expression profiling, we identified and validated a set of genes involved in nuchal hump formation in the Lake Malawi dolphin cichlid, Cyrtocara moorii. In particular, we found genes differentially expressed in the nuchal hump, which are involved in controlling cell proliferation (btg3, fosl1a and pdgfrb), cell growth (dlk1), craniofacial morphogenesis (dlx5a, mycn and tcf12), as well as regulators of growth-related signals (dpt, pappa and socs2). This is the first study to identify the set of genes associated with nuchal hump formation in cichlids. Given that the hump is a trait that evolved repeatedly in several African and American cichlid lineages, it would be interesting to see if the molecular pathways and genes triggering hump formation follow a common genetic track or if the trait evolved in parallel, with distinct mechanisms, in other cichlid adaptive radiations and even in other teleost fishes.

Assessment of candidate genes and genetic heterogeneity in human non syndromic orofacial clefts specifically non syndromic cleft lip with or without palate

Non syndromic orofacial clefts specifically non-syndromic cleft lip/palate are one of the most common craniofacial malformation among birth defects in human having multifactorial etiology with an incidence of 1:700/1000. On the basis of association with other congenital malformations or their presence as isolated anomaly, OFC can be classified as syndromic (30%) and nonsyndromic (70%) respectively. The major cause of disease demonstrates complex interplay between genetic and environmental factors. The pathogenic mechanism of underlying factors have been provided by different genetic studies on large-scale with significant recent advances in genotyping technologies usually based on linkage or genome wide association studies (GWAS). On the basis of recent studies, new tools to identify causative genes involved in NSCL/P reported approximately more than 30 genetic risk loci that are responsible for pathogenesis of facial deformation. Despite these findings, it is still uncertain that how much of variance in NSCL/P predisposing factors can be explain by identified risk loci, as they all together accounts for only 20%-25% of NSCL/P heritability. So there is need of further findings about the problem of rare low frequency coding variants and other missing responsive factors or genetic modifiers. This review will described those potential genes and loci reported in different studies whose involvement in pathogenesis of nonsyndromic OFC has wide scientific evidence.

Functional morphogenesis from embryos to adults: Late development shapes trophic niche in coral reef damselfishes

The damselfishes are one of the dominant coral reef fish lineages. Their ecological diversification has involved repeated transitions between pelagic feeding using fast bites and benthic feeding using forceful bites. A highly-integrative approach that combined gene expression assays, shape analyses, and high-speed video analyses was used to examine the development of trophic morphology in embryonic, larval, juvenile, and adult damselfishes. The anatomical characters that distinguish pelagic-feeding and benthic-feeding species do not appear until after larval development. Neither patterns of embryonic jaw morphogenesis, larval skull shapes nor larval bite mechanics significantly distinguished damselfishes from different adult trophic guilds. Analyses of skull shape and feeding performance identified two important transitions in the trophic development of a single species (the orange clownfish; Amphiprion percula): (a) a pronounced transformation in feeding mechanics during metamorphosis; and (b) more protracted cranial remodeling over the course of juvenile development. The results of this study indicate that changes in postlarval morphogenesis have played an important role in damselfish evolution. This is likely to be true for other fish lineages, particularly if they consist of marine species, the majority of which have planktonic larvae with different functional requirements for feeding in comparison to their adult forms.

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.

Zebrafish Models of Human Skeletal Disorders: Embryo and Adult Swimming Together

Danio rerio (zebrafish) is an elective model organism for the study of vertebrate development because of its high degree of homology with human genes and organs, including bone. Zebrafish embryos, because of the optical clarity, small size, and fast development, can be easily used in large-scale mutagenesis experiments to isolate mutants with developmental skeletal defects and in high-throughput screenings to find new chemical compounds for the ability to revert the pathological phenotype. On the other hand, the adult zebrafish represents another powerful resource for pathogenic and therapeutic studies about adult human bone diseases. In fish, some characteristics such as bone turnover, reparation, and remodeling of the adult bone tissue cannot be found at the embryonic stage. Several pathological models have been established in adult zebrafish such as bone injury models, osteoporosis, and genetic diseases such as osteogenesis imperfecta. Given the growing interest for metabolic diseases and their complications, adult zebrafish models of type 2 diabetes and obesity have been recently generated and analyzed for bone complications using scales as model system. Interestingly, an osteoporosis-like phenotype has been found to be associated with metabolic alterations suggesting that bone complications share the same mechanisms in humans and fish. Embryo and adult represent powerful resources in rapid development to study bone physiology and pathology from different points of view.

Muscle precursor cell movements in zebrafish are dynamic and require Six family genes

Muscle precursors need to be correctly positioned during embryonic development for proper body movement. In zebrafish, a subset of hypaxial muscle precursors from the anterior somites undergo long-range migration, moving away from the trunk in three streams to form muscles in distal locations such as the fin. We mapped long-distance muscle precursor migrations with unprecedented resolution using live imaging. We identified conserved genes necessary for normal precursor motility (six1a, six1b, six4a, six4b and met). These genes are required for movement away from somites and later to partition two muscles within the fin bud. During normal development, the middle muscle precursor stream initially populates the fin bud, then the remainder of this stream contributes to the posterior hypaxial muscle. When we block fin bud development by impairing retinoic acid synthesis or Fgfr function, the entire stream contributes to the posterior hypaxial muscle indicating that muscle precursors are not committed to the fin during migration. Our findings demonstrate a conserved muscle precursor motility pathway, identify dynamic cell movements that generate posterior hypaxial and fin muscles, and demonstrate flexibility in muscle precursor fates.

hdac4 mediates perichondral ossification and pharyngeal skeleton development in the zebrafish

Background

Histone deacetylases (HDACs) are epigenetic factors that function to repress gene transcription by removing acetyl groups from the N-terminal of histone lysines. Histone deacetylase 4 (HDAC4), a class IIa HDAC, has previously been shown to regulate the process of endochondral ossification in mice via repression of Myocyte enhancer factor 2c (MEF2C), a transcriptional activator of Runx2, which in turn promotes chondrocyte maturation and production of bone by osteoblasts.

Methods & Materials

In this study, we generated two zebrafish lines with mutations in hdac4 using CRISPR/Cas9 and analyzed mutants for skeletal phenotypes and expression of genes known to be affected by Hdac4 expression.

Results

Lines have insertions causing a frameshift in a proximal exon of hdac4 and a premature stop codon. Mutations are predicted to result in aberrant protein sequence and a truncated protein, eliminating the Mef2c binding domain and Hdac domain. Zygotic mutants from two separate lines show a significant increase in ossification of pharyngeal ceratohyal cartilages at 7 days post fertilization (dpf) (p < 0.01, p < 0.001). At 4 dpf, mutant larvae have a significant increase of expression of runx2a and runx2b in the ceratohyal cartilage (p < 0.05 and p < 0.01, respectively). A subset of maternal-zygotic (mz) mutant and heterozygote larvae (40%) have dramatically increased ossification at 7 dpf compared to zygotic mutants, including formation of a premature anguloarticular bone and mineralization of the first and second ceratobranchial cartilages and symplectic cartilages, which normally does not occur until fish are approximately 10 or 12 dpf. Some maternal-zygotic mutants and heterozygotes show loss of pharyngeal first arch elements (25.9% and 10.2%, respectively) and neurocranium defects (30.8% and 15.2%, respectively). Analysis of RNA-seq mRNA transcript levels and in situ hybridizations from zygotic stages to 75–90% epiboly indicates that hdac4 is highly expressed in early embryos, but diminishes by late epiboly, becoming expressed again in larval stages.

Discussion

Loss of function of hdac4 in zebrafish is associated with increased expression of runx2a and runx2b targets indicating that a role for hdac4 in zebrafish is to repress activation of ossification of cartilage. These findings are consistent with observations of precocious cartilage ossification in Hdac4 mutant mice, demonstrating that the function of Hdac4 in skeletal development is conserved among vertebrates. Expression of hdac4 mRNA in embryos younger than 256–512 cells indicates that there is a maternal contribution of hdac4 to the early embryo. The increase in ossification and profound loss of first pharyngeal arch elements and anterior neurocranium in a subset of maternal-zygotic mutant and heterozygote larvae suggests that maternal hdac4 functions in cartilage ossification and development of cranial neural crest-derived structures.

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.

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.

Functional testing of a human PBX3 variant in zebrafish reveals a potential modifier role in congenital heart defects

Whole-genome and exome sequencing efforts are increasingly identifying candidate genetic variants associated with human disease. However, predicting and testing the pathogenicity of a genetic variant remains challenging. Genome editing allows for the rigorous functional testing of human genetic variants in animal models. Congenital heart defects (CHDs) are a prominent example of a human disorder with complex genetics. An inherited sequence variant in the human PBX3 gene (PBX3 p.A136V) has previously been shown to be enriched in a CHD patient cohort, indicating that the PBX3 p.A136V variant could be a modifier allele for CHDs. Pbx genes encode three-amino-acid loop extension (TALE)-class homeodomain-containing DNA-binding proteins with diverse roles in development and disease, and are required for heart development in mouse and zebrafish. Here, we used CRISPR-Cas9 genome editing to directly test whether this Pbx gene variant acts as a genetic modifier in zebrafish heart development. We used a single-stranded oligodeoxynucleotide to precisely introduce the human PBX3 p.A136V variant in the homologous zebrafish pbx4 gene (pbx4 p.A131V). We observed that zebrafish that are homozygous for pbx4 p.A131V are viable as adults. However, the pbx4 p.A131V variant enhances the embryonic cardiac morphogenesis phenotype caused by loss of the known cardiac specification factor, Hand2. Our study is the first example of using precision genome editing in zebrafish to demonstrate a function for a human disease-associated single nucleotide variant of unknown significance. Our work underscores the importance of testing the roles of inherited variants, not just de novo variants, as genetic modifiers of CHDs. Our study provides a novel approach toward advancing our understanding of the complex genetics of CHDs.

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.

Bapx1 upregulation is associated with ectopic mandibular cartilage development in amphibians

BACKGROUND

The emergence of novel structures during evolution is crucial for creating variation among organisms, but the underlying processes which lead to the emergence of evolutionary novelties are poorly understood. The gnathostome jaw joint is such a novelty, and the incorporation of bapx1 expression into the intermediate first pharyngeal arch may have played a major role in the evolution of this joint. Knockdown experiments revealed that loss of bapx1 function leads to the loss of the jaw joint, because Meckel's cartilage and the palatoquadrate fuse during development. We used Xenopus laevis and Ambystoma mexicanum to further investigate the function of bapx1 in amphibians. Bapx1 expression levels were upregulated through the use of Ly-294,002 and we investigated the potential consequences of the enhanced bapx1 expression in amphibians to test the hypothesized joint inducing function of bapx1.

RESULTS

We show that Ly-294,002 upregulates bapx1 expression in vivo. Additionally, ectopic mandibular arch derived cartilages develop after Ly-294,002 treatment. These ectopic cartilages are dorsoventrally oriented rods situated lateral to the palatoquadrate. The development of these additional cartilages did not change the muscular arrangement of mandibular arch-derived muscles.

CONCLUSIONS

Development of additional mandibular cartilages is not unusual in larval anurans. Therefore, changes in the bapx1 expression during evolution may have been the reason for the development of several additional cartilages in the larval anuran jaw. Furthermore, our observations imply a joint-promoting function of bapx1, which further substantiates its hypothetical role in the evolution of the gnathostome jaw joint.