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

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

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

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

Distinct and redundant roles for zebrafish her genes during mineralization and craniofacial patterning.

The Notch pathway is a cell-cell communication system which is critical for many developmental processes, including craniofacial development. Notch receptor activation induces expression of several well-known canonical targets including those encoded by the hes and her genes in mammals and zebrafish, respectively. The function of these genes, individually and in combination, during craniofacial development is not well understood. Here, we used zebrafish genetics to investigate her9 and her6 gene function during craniofacial development. We found that her9 is required for osteoblasts to efficiently mineralize bone, while cartilage is largely unaffected. Strikingly, gene expression studies in her9 mutants indicate that although progenitor cells differentiate into osteoblasts at the appropriate time and place, they fail to efficiently lay down mineralized matrix. This mineralization role of her9 is likely independent of Notch activation. In contrast, her9 also functions redundantly with her6 downstream of Jagged1b-induced Notch activation during dorsoventral craniofacial patterning. These studies disentangle distinct and redundant her gene functions during craniofacial development, including an unexpected, Notch independent, requirement during bone mineralization.

Variable paralog expression underlies phenotype variation.

Human faces are variable; we look different from one another. Craniofacial disorders further increase facial variation. To understand craniofacial variation and how it can be buffered, we analyzed the zebrafish mef2ca mutant. When this transcription factor encoding gene is mutated, zebrafish develop dramatically variable craniofacial phenotypes. Years of selective breeding for low and high penetrance of mutant phenotypes produced strains that are either resilient or sensitive to the mef2ca mutation. Here, we compared gene expression between these strains, which revealed that selective breeding enriched for high and low mef2ca paralog expression in the low- and high-penetrance strains, respectively. We found that mef2ca paralog expression is variable in unselected wild-type zebrafish, motivating the hypothesis that heritable variation in paralog expression underlies mutant phenotype severity and variation. In support, mutagenizing the mef2ca paralogs, mef2aa, mef2b, mef2cb, and mef2d demonstrated modular buffering by paralogs. Specifically, some paralogs buffer severity while others buffer variability. We present a novel, mechanistic model for phenotypic variation where variable, vestigial paralog expression buffers development. These studies are a major step forward in understanding the mechanisms of facial variation, including how some genetically resilient individuals can overcome a deleterious mutation.

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.

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.

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.

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.

Shifting Zebrafish Lethal Skeletal Mutant Penetrance by Progeny Testing.

Zebrafish mutant phenotypes are often incompletely penetrant, only manifesting in some mutants. Interesting phenotypes that inconsistently appear can be difficult to study, and can lead to confounding results. The protocol described here is a straightforward breeding paradigm to increase and decrease penetrance in lethal zebrafish skeletal mutants. Because lethal mutants cannot be selectively bred directly, the classic selective breeding strategy of progeny testing is employed. This method also includes protocols for Kompetitive Allele Specific PCR (KASP) genotyping zebrafish and staining larval zebrafish cartilage and bone. Applying the husbandry strategy described here can increase the penetrance of an interesting skeletal phenotype enabling more reproducible results in downstream applications. In addition, decreasing the mutant penetrance through this selective breeding strategy can reveal the developmental processes that most crucially require the function of the mutated gene. While the skeleton is specifically considered here, we propose that this methodology will be useful for all zebrafish mutant lines.

Ligament versus bone cell identity in the zebrafish hyoid skeleton is regulated by mef2ca.

Heightened phenotypic variation among mutant animals is a well-known, but poorly understood phenomenon. One hypothetical mechanism accounting for mutant phenotypic variation is progenitor cells variably choosing between two alternative fates during development. Zebrafish mef2cab1086 mutants develop tremendously variable ectopic bone in their hyoid craniofacial skeleton. Here, we report evidence that a key component of this phenotype is variable fate switching from ligament to bone. We discover that a 'track' of tissue prone to become bone cells is a previously undescribed ligament. Fate-switch variability is heritable, and comparing mutant strains selectively bred to high and low penetrance revealed differential mef2ca mutant transcript expression between high and low penetrance strains. Consistent with this, experimental manipulation of mef2ca mutant transcripts modifies the penetrance of the fate switch. Furthermore, we discovered a transposable element that resides immediately upstream of the mef2ca locus and is differentially DNA methylated in the two strains, correlating with differential mef2ca expression. We propose that variable transposon epigenetic silencing underlies the variable mef2ca mutant bone phenotype, and could be a widespread mechanism of phenotypic variability in animals.

Pharyngeal morphogenesis requires fras1-itga8-dependent epithelial-mesenchymal interaction.

Both Fras1 and Itga8 connect mesenchymal cells to epithelia by way of an extracellular 'Fraser protein complex' that functions in signaling and adhesion; these proteins are vital to the development of several vertebrate organs. We previously found that zebrafish fras1 mutants have craniofacial defects, specifically, shortened symplectic cartilages and cartilage fusions that spare joint elements. During a forward mutagenesis screen, we identified a new zebrafish mutation, b1161, that we show here disrupts itga8, as confirmed using CRISPR-generated itga8 alleles. fras1 and itga8 single mutants and double mutants have similar craniofacial phenotypes, a result expected if loss of either gene disrupts function of the Fraser protein complex. Unlike fras1 mutants or other Fraser-related mutants, itga8 mutants do not show blistered tail fins. Thus, the function of the Fraser complex differs in the craniofacial skeleton and the tail fin. Focusing on the face, we find that itga8 mutants consistently show defective outpocketing of a late-forming portion of the first pharyngeal pouch, and variably express skeletal defects, matching previously characterized fras1 mutant phenotypes. In itga8 and fras1 mutants, skeletal severity varies markedly between sides, indicating that both mutants have increased developmental instability. Whereas fras1 is expressed in epithelia, we show that itga8 is expressed complementarily in facial mesenchyme. Paired with the observed phenotypic similarity, this expression indicates that the genes function in epithelial-mesenchymal interactions. Similar interactions between Fras1 and Itga8 have previously been found in mouse kidney, where these genes both regulate Nephronectin (Npnt) protein abundance. We find that zebrafish facial tissues express both npnt and the Fraser gene fibrillin2b (fbn2b), but their transcript levels do not depend on fras1 or itga8 function. Using a revertible fras1 allele, we find that the critical window for fras1 function in the craniofacial skeleton is between 1.5 and 3 days post fertilization, which coincides with the onset of fras1-dependent and itga8-dependent morphogenesis. We propose a model wherein Fras1 and Itga8 interact during late pharyngeal pouch morphogenesis to sculpt pharyngeal arches through epithelial-mesenchymal interactions, thereby stabilizing the developing craniofacial skeleton.

Competition between Jagged-Notch and Endothelin1 Signaling Selectively Restricts Cartilage Formation in the Zebrafish Upper Face.

The intricate shaping of the facial skeleton is essential for function of the vertebrate jaw and middle ear. While much has been learned about the signaling pathways and transcription factors that control facial patterning, the downstream cellular mechanisms dictating skeletal shapes have remained unclear. Here we present genetic evidence in zebrafish that three major signaling pathways - Jagged-Notch, Endothelin1 (Edn1), and Bmp - regulate the pattern of facial cartilage and bone formation by controlling the timing of cartilage differentiation along the dorsoventral axis of the pharyngeal arches. A genomic analysis of purified facial skeletal precursors in mutant and overexpression embryos revealed a core set of differentiation genes that were commonly repressed by Jagged-Notch and induced by Edn1. Further analysis of the pre-cartilage condensation gene barx1, as well as in vivo imaging of cartilage differentiation, revealed that cartilage forms first in regions of high Edn1 and low Jagged-Notch activity. Consistent with a role of Jagged-Notch signaling in restricting cartilage differentiation, loss of Notch pathway components resulted in expanded barx1 expression in the dorsal arches, with mutation of barx1 rescuing some aspects of dorsal skeletal patterning in jag1b mutants. We also identified prrx1a and prrx1b as negative Edn1 and positive Bmp targets that function in parallel to Jagged-Notch signaling to restrict the formation of dorsal barx1+ pre-cartilage condensations. Simultaneous loss of jag1b and prrx1a/b better rescued lower facial defects of edn1 mutants than loss of either pathway alone, showing that combined overactivation of Jagged-Notch and Bmp/Prrx1 pathways contribute to the absence of cartilage differentiation in the edn1 mutant lower face. These findings support a model in which Notch-mediated restriction of cartilage differentiation, particularly in the second pharyngeal arch, helps to establish a distinct skeletal pattern in the upper face.

Epigenetic regulation of hematopoiesis by DNA methylation.

During embryonic development, cell type-specific transcription factors promote cell identities, while epigenetic modifications are thought to contribute to maintain these cell fates. Our understanding of how genetic and epigenetic modes of regulation work together to establish and maintain cellular identity is still limited, however. Here, we show that DNA methyltransferase 3bb.1 (dnmt3bb.1) is essential for maintenance of hematopoietic stem and progenitor cell (HSPC) fate as part of an early Notch-runx1-cmyb HSPC specification pathway in the zebrafish. Dnmt3bb.1 is expressed in HSPC downstream from Notch1 and runx1, and loss of Dnmt3bb.1 activity leads to reduced cmyb locus methylation, reduced cmyb expression, and gradual reduction in HSPCs. Ectopic overexpression of dnmt3bb.1 in non-hematopoietic cells is sufficient to methylate the cmyb locus, promote cmyb expression, and promote hematopoietic development. Our results reveal an epigenetic mechanism supporting the maintenance of hematopoietic cell fate via DNA methylation-mediated perdurance of a key transcription factor in HSPCs.

Patterns of variation and covariation in the shapes of mandibular bones of juvenile salmonids in the genus Oncorhynchus.

What is the nature of evolutionary divergence of the jaw skeleton within the genus Oncorhynchus? How can two associated bones evolve new shapes and still maintain functional integration? Here, we introduce and test a "concordance" hypothesis, in which an extraordinary matching of the evolutionary shape changes of the dentary and angular articular serves to preserve their fitting together. To test this hypothesis, we examined morphologies of the dentary and angular articular at parr (juvenile) stage, and at three levels of biological organization­between salmon and trout, between sister species within both salmon and trout, and among three types differing in life histories within one species, Oncorhynchus mykiss. The comparisons show bone shape divergences among the groups at each level; morphological divergence between salmon and trout is marked even at this relatively early life history stage. We observed substantial matching between the two mandibular bones in both pattern and amount of shape variation, and in shape covariation across species. These findings strongly support the concordance hypothesis, and reflect functional and/or developmental constraint on morphological evolution. We present evidence for developmental modularity within both bones. The locations of module boundaries were predicted from the patterns of evolutionary divergences, and for the dentary, at least, would appear to facilitate its functional association with the angular articular. The modularity results suggest that development has biased the course of evolution.

Role of mef2ca in developmental buffering of the zebrafish larval hyoid dermal skeleton.

Phenotypic robustness requires a process of developmental buffering that is largely not understood, but which can be disrupted by mutations. Here we show that in mef2ca(b1086) loss of function mutant embryos and early larvae, development of craniofacial hyoid bones, the opercle (Op) and branchiostegal ray (BR), becomes remarkably unstable; the large magnitude of the instability serves as a positive attribute to learn about features of this developmental buffering. The OpBR mutant phenotype variably includes bone expansion and fusion, Op duplication, and BR homeosis. Formation of a novel bone strut, or a bone bridge connecting the Op and BR together occurs frequently. We find no evidence that the phenotypic stability in the wild type is provided by redundancy between mef2ca and its co-ortholog mef2cb, or that it is related to the selector (homeotic) gene function of mef2ca. Changes in dorsal-ventral patterning of the hyoid arch also might not contribute to phenotypic instability in mutants. However, subsequent development of the bone lineage itself, including osteoblast differentiation and morphogenetic outgrowth, shows marked variation. Hence, steps along the developmental trajectory appear differentially sensitive to the loss of buffering, providing focus for the future study.

edn1 and hand2 Interact in early regulation of pharyngeal arch outgrowth during zebrafish development.

Endothelin-1 (Edn1) signaling provides a critical input to development of the embryonic pharygneal arches and their skeletal derivatives, particularly the articulating joints and the ventral skeleton including the lower jaw. Previous work in zebrafish has mostly focused on the role of Edn1 in dorsal-ventral (DV) patterning, but Edn1 signaling must also regulate tissue size, for with severe loss of the pathway the ventral skeleton is not only mispatterned, but is also prominently hypoplastic--reduced in size. Here we use mutational analyses to show that in the early pharyngeal arches, ventral-specific edn1-mediated proliferation of neural crest derived cells is required for DV expansion and outgrowth, and that this positive regulation is counterbalanced by a negative one exerted through a pivotal, ventrally expressed Edn1-target gene, hand2. We also describe a new morphogenetic cell movement in the ventral first arch, sweeping cells anterior in the arch to the region where the lower jaw forms. This movement is negatively regulated by hand2 in an apparently edn1-independent fashion. These findings point to complexity of regulation by edn1 and hand2 at the earliest stages of pharyngeal arch development, in which control of growth and morphogenesis can be genetically separated.

barx1 represses joints and promotes cartilage in the craniofacial skeleton.

The evolution of joints, which afford skeletal mobility, was instrumental in vertebrate success. Here, we explore the molecular genetics and cell biology that govern jaw joint development. Genetic manipulation experiments in zebrafish demonstrate that functional loss, or gain, of the homeobox-containing gene barx1 produces gain, or loss, of joints, respectively. Ectopic joints in barx1 mutant animals are present in every pharyngeal segment, and are associated with disrupted attachment of bone, muscles and teeth. We find that ectopic joints develop at the expense of cartilage. Time-lapse experiments suggest that barx1 controls the skeletal precursor cell choice between differentiating into cartilage versus joint cells. We discovered that barx1 functions in this choice, in part, by regulating the transcription factor hand2. We further show that hand2 feeds back to negatively regulate barx1 expression. We consider the possibility that changes in barx1 function in early vertebrates were among the key innovations fostering the evolution of skeletal joints.

FishFace: interactive atlas of zebrafish craniofacial development at cellular resolution.

BACKGROUND: The vertebrate craniofacial skeleton may exhibit anatomical complexity and diversity, but its genesis and evolution can be understood through careful dissection of developmental programs at cellular resolution. Resources are lacking that include introductory overviews of skeletal anatomy coupled with descriptions of craniofacial development at cellular resolution. In addition to providing analytical guidelines for other studies, such an atlas would suggest cellular mechanisms underlying development. DESCRIPTION: We present the Fish Face Atlas, an online, 3D-interactive atlas of craniofacial development in the zebrafish Danio rerio. Alizarin red-stained skulls scanned by fluorescent optical projection tomography and segmented into individual elements provide a resource for understanding the 3D structure of the zebrafish craniofacial skeleton. These data provide the user an anatomical entry point to confocal images of Alizarin red-stained zebrafish with transgenically-labelled pharyngeal arch ectomesenchyme, chondrocytes, and osteoblasts, which illustrate the appearance, morphogenesis, and growth of the mandibular and hyoid cartilages and bones, as viewed in live, anesthetized zebrafish during embryonic and larval development. Confocal image stacks at high magnification during the same stages provide cellular detail and suggest developmental and evolutionary hypotheses. CONCLUSION: The FishFace Atlas is a novel learning tool for understanding craniofacial skeletal development, and can serve as a reference for a variety of studies, including comparative and mutational analyses.

Notch signaling--constantly on the move.

The Notch pathway is a highly conserved and ubiquitous signaling system that functions in determining a diverse array of cell fates and regulates many cellular processes during embryonic development and throughout adulthood. Links to cancer, stroke and Alzheimer's disease underscore the need to define the molecular basis of Notch activation. Notch signaling is induced through direct cell-cell interactions that promote receptor activation following engagement with a membrane-bound Delta, Serrate, Lag-2 (DSL) ligand on adjacent cells. Cells take on distinct fates because Notch signaling is consistently activated in only one of the two interacting cells, highlighting the importance of establishing and maintaining signaling polarity. Studies in flies and worms have identified positive and negative transcriptional feedback mechanisms that amplify small differences in Notch and DSL ligand expression to bias which cells send or receive signals. However, endocytosis by signal-sending and signal-receiving cells also appears critical for directing and regulating Notch activation. In particular, endocytosis and membrane trafficking of DSL ligands, Notch and modulators can determine the competence of cells to send or receive signals that ensure reproducibility in generating cell types regulated by Notch signaling.