Zebrafish furin mutants reveal intricacies in regulating Endothelin1 signaling in craniofacial patterning

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.

Tending to the Facial Surfaces of a Mathematical Biology Head-Scratcher: Why Does the Head of the Sea Turtle Natator depressus Resemble a Convex Zygomorphic Dodecahedron?

Two convex polyhedra that markedly resemble the head of the flatback sea turtle hatchling are identified. The first example is a zygomorphic tetragonal dodecahedron, while the other, an even better matching structure, is a related tetradecahedron, herein speculated to arise from this particular dodecahedron via known mechanisms gleaned from studies of the behavior of foams. A segmented, biomorphic, convex polyhedral model to address cephalic topology is thus presented stemming from solid geometry, anatomical observations, and a recently computed densest local packing arrangement of fifteen slightly oblate spheroids in which fourteen oblate spheroids surround a central such spheroid. This particular array of oblate spheroids shares salient structural features with the aforementioned dodecahedron. Successful testing of the model has been achieved by converting this array of fifteen oblate spheroids constructed with putty to the cephaloid dodecahedron in a process involving ventral elongation induced by stretching in the anterior direction along the anteroposterior axis (convergent extension). During convergent extension, the two left most anterolateral oblate spheroids that are in direct contact with the ventral spheroid of the array merge into a single lateral facet of the incipient dodecahedron, while the corresponding two right such oblate spheroids do the same. Thus, the fourteen outer oblate spheroids of the array give rise to the twelve facets of the finalized dodecahedron, while the central oblate spheroid remnant assumes an interior dodecahedral position. The hypothetical dodecahedron to tetradecahedron transformation entails the collapse of a tetravalent vertex (which is known to occur in foams as part of a T1 transition) followed by bilateral facet splitting. Remarkably, a model stipulating that convexity is to be retained in connection with this sequence of steps necessitates that the starting dodecahedral template undergoes modification to become a tetradecahedron in possession of precisely the highly ordered feature found at the top of the head of numerous specimens of the flatback sea turtle hatchling, namely, a fused medial pentagon-heptagon pair in the form of a pentagonal frontal scute and heptagonal frontoparietal scute. Such a possible new instance of geometric biomorphy, taken together with the correct anticipation of the cephalic pentagon-heptagon pair, might serve to instill further confidence in renewed efforts to shed light on morphogenesis with foam embryo models.

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.

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.

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.

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.

Translational control of furina by an RNA regulon is important for left-right patterning, heart morphogenesis and cardiac valve function

Heart development is a complex process that requires asymmetric positioning of the heart, cardiac growth and valve morphogenesis. The mechanisms controlling heart morphogenesis and valve formation are not fully understood. The pro-convertase FurinA functions in heart development across vertebrates. How FurinA activity is regulated during heart development is unknown. Through computational analysis of the zebrafish transcriptome, we identified an RNA motif in a variant FurinA transcript harbouring a long 3' untranslated region (3'UTR). The alternative 3'UTR furina isoform is expressed prior to organ positioning. Somatic deletions in the furina 3'UTR lead to embryonic left-right patterning defects. Reporter localisation and RNA-binding assays show that the furina 3'UTR forms complexes with the conserved RNA-binding translational repressor, Ybx1. Conditional ybx1 mutant embryos show premature and increased Furin reporter expression, abnormal cardiac morphogenesis and looping defects. Mutant ybx1 hearts have an expanded atrioventricular canal, abnormal sino-atrial valves and retrograde blood flow from the ventricle to the atrium. This is similar to observations in humans with heart valve regurgitation. Thus, the furina 3'UTR element/Ybx1 regulon is important for translational repression of FurinA and regulation of heart development.

A novel insight on SARS-CoV-2 S-derived fragments in the control of the host immunity

Despite all efforts to combat the pandemic of COVID-19, we are still living with high numbers of infected persons, an overburdened health care system, and the lack of an effective and definitive treatment. Understanding the pathophysiology of the disease is crucial for the development of new technologies and therapies for the best clinical management of patients. Since the manipulation of the whole virus requires a structure with an adequate level of biosafety, the development of alternative technologies, such as the synthesis of peptides from viral proteins, is a possible solution to circumvent this problem. In addition, the use and validation of animal models is of extreme importance to screen new drugs and to compress the organism's response to the disease. Peptides derived from recombinant S protein from SARS-CoV-2 were synthesized and validated by in silico, in vitro and in vivo methodologies. Macrophages and neutrophils were challenged with the peptides and the production of inflammatory mediators and activation profile were evaluated. These peptides were also inoculated into the swim bladder of transgenic zebrafish larvae at 6 days post fertilization (dpf) to mimic the inflammatory process triggered by the virus, which was evaluated by confocal microscopy. In addition, toxicity and oxidative stress assays were also developed. In silico and molecular dynamics assays revealed that the peptides bind to the ACE2 receptor stably and interact with receptors and adhesion molecules, such as MHC and TCR, from humans and zebrafish. Macrophages stimulated with one of the peptides showed increased production of NO, TNF-α and CXCL2. Inoculation of the peptides in zebrafish larvae triggered an inflammatory process marked by macrophage recruitment and increased mortality, as well as histopathological changes, similarly to what is observed in individuals with COVID-19. The use of peptides is a valuable alternative for the study of host immune response in the context of COVID-19. The use of zebrafish as an animal model also proved to be appropriate and effective in evaluating the inflammatory process, comparable to humans.

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.

Control of cranial ectomesenchyme fate by Nr2f nuclear receptors

Certain cranial neural crest cells are uniquely endowed with the ability to make skeletal cell types otherwise only derived from mesoderm. As these cells migrate into the pharyngeal arches, they downregulate neural crest specifier genes and upregulate so-called ectomesenchyme genes that are characteristic of skeletal progenitors. Although both external and intrinsic factors have been proposed as triggers of this transition, the details remain obscure. Here, we report the Nr2f nuclear receptors as intrinsic activators of the ectomesenchyme program: zebrafish nr2f5 single and nr2f2;nr2f5 double mutants show marked delays in upregulation of ectomesenchyme genes, such as dlx2a, prrx1a, prrx1b, sox9a, twist1a and fli1a, and in downregulation of sox10, which is normally restricted to early neural crest and non-ectomesenchyme lineages. Mutation of sox10 fully rescued skeletal development in nr2f5 single but not nr2f2;nr2f5 double mutants, but the initial ectomesenchyme delay persisted in both. Sox10 perdurance thus antagonizes the recovery but does not explain the impaired ectomesenchyme transition. Unraveling the mechanisms of Nr2f function will help solve the enduring puzzle of how cranial neural crest cells transition to the skeletal progenitor state.

A missense mutation in the proprotein convertase gene furinb causes hepatic cystogenesis during liver development in zebrafish

Hepatic cysts are fluid-filled lesions in the liver that are estimated to occur in 5% of the population. They may cause hepatomegaly and abdominal pain. Progression to secondary fibrosis, cirrhosis, or cholangiocarcinoma can lead to morbidity and mortality. Previous studies of patients and rodent models have associated hepatic cyst formation with increased proliferation and fluid secretion in cholangiocytes, which are partially due to impaired primary cilia. Congenital hepatic cysts are thought to originate from faulty bile duct development, but the underlying mechanisms are not fully understood. In a forward genetic screen, we identified a zebrafish mutant that developed hepatic cysts during larval stages. The cyst formation was not due to changes in biliary cell proliferation, bile secretion, or impairment of primary cilia. Instead, time-lapse live imaging data showed that the mutant biliary cells failed to form interconnecting bile ducts because of defects in motility and protrusive activity. Accordingly, immunostaining revealed a disorganized actin and microtubule cytoskeleton in the mutant biliary cells. By whole-genome sequencing, we determined that the cystic phenotype in the mutant was caused by a missense mutation in the furinb gene, which encodes a proprotein convertase. The mutation altered Furinb localization and caused endoplasmic reticulum (ER) stress. The cystic phenotype could be suppressed by treatment with the ER stress inhibitor 4-phenylbutyric acid and exacerbated by treatment with the ER stress inducer tunicamycin. The mutant liver also exhibited increased mammalian target of rapamycin (mTOR) signaling. Treatment with mTOR inhibitors halted cyst formation at least partially through reducing ER stress. Conclusion: Our study has established a vertebrate model for studying hepatic cystogenesis and illustrated the contribution of ER stress in the disease pathogenesis.

Adar-mediated A-to-I editing is required for embryonic patterning and innate immune response regulation in zebrafish

Adenosine deaminases (ADARs) catalyze the deamination of adenosine to inosine, also known as A-to-I editing, in RNA. Although A-to-I editing occurs widely across animals and is well studied, new biological roles are still being discovered. Here, we study the role of A-to-I editing in early zebrafish development. We demonstrate that Adar, the zebrafish orthologue of mammalian ADAR1, is essential for establishing the antero-posterior and dorso-ventral axes and patterning. Genome-wide editing discovery reveals pervasive editing in maternal and the earliest zygotic transcripts, the majority of which occurred in the 3'-UTR. Interestingly, transcripts implicated in gastrulation as well as dorso-ventral and antero-posterior patterning are found to contain multiple editing sites. Adar knockdown or overexpression affect gene expression by 12 hpf. Analysis of adar-/- zygotic mutants further reveals that the previously described role of Adar in mammals in regulating the innate immune response is conserved in zebrafish. Our study therefore establishes distinct maternal and zygotic functions of RNA editing by Adar in embryonic patterning along the zebrafish antero-posterior and dorso-ventral axes, and in the regulation of the innate immune response, respectively.

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.

Involvement of the zebrafish trrap gene in craniofacial development

Trrap (transformation/transcription domain-associated protein) is a component shared by several histone acetyltransferase (HAT) complexes and participates in transcriptional regulation and DNA repair; however, the developmental functions of Trrap in vertebrates are not fully understood. Recently, it has been reported that human patients with genetic mutations in the TRRAP gene show various symptoms, including facial dysmorphisms, microcephaly and global developmental delay. To investigate the physiological functions of Trrap, we established trrap gene-knockout zebrafish and examined loss-of-function phenotypes in the mutants. The trrap zebrafish mutants exhibited smaller eyes and heads than the wild-type zebrafish. The size of the ventral pharyngeal arches was reduced and the mineralization of teeth was impaired in the trrap mutants. Whole-mount in situ hybridization analysis revealed that dlx3 expression was narrowly restricted in the developing ventral pharyngeal arches, while dlx2b expression was diminished in the trrap mutants. These results suggest that trrap zebrafish mutants are useful model organisms for a human disorder associated with genetic mutations in the human TRRAP gene.

Protease-dependent defects in N-cadherin processing drive PMM2-CDG pathogenesis

The genetic bases for the congenital disorders of glycosylation (CDG) continue to expand, but how glycosylation defects cause patient phenotypes remains largely unknown. Here, we combined developmental phenotyping and biochemical studies in a potentially new zebrafish model (pmm2^sa10150) of PMM2-CDG to uncover a protease-mediated pathogenic mechanism relevant to craniofacial and motility phenotypes in mutant embryos. Mutant embryos had reduced phosphomannomutase activity and modest decreases in N-glycan occupancy as detected by matrix-assisted laser desorption ionization mass spectrometry imaging. Cellular analyses of cartilage defects in pmm2^sa10150 embryos revealed a block in chondrogenesis that was associated with defective proteolytic processing, but seemingly normal N-glycosylation, of the cell adhesion molecule N-cadherin. The activities of the proconvertases and matrix metalloproteinases responsible for N-cadherin maturation were significantly altered in pmm2^sa10150 mutant embryos. Importantly, pharmacologic and genetic manipulation of proconvertase activity restored matrix metalloproteinase activity, N-cadherin processing, and cartilage pathology in pmm2^sa10150 embryos. Collectively, these studies demonstrate in CDG that targeted alterations in protease activity create a pathogenic cascade that affects the maturation of cell adhesion proteins critical for tissue development.

Proprotein convertase furina is required for heart development in zebrafish

Cardiac looping and trabeculation are key processes during cardiac chamber maturation. However, the underlying mechanisms remain incompletely understood. Here, we report the isolation, cloning and characterization of the proprotein convertase furina from the cardiovascular mutant loft in zebrafish. loft is an ethylnitrosourea-induced mutant and has evident defects in the cardiac outflow tract, heart looping and trabeculation, the craniofacial region and pharyngeal arch arteries. Positional cloning revealed that furina mRNA was barely detectable in loft mutants, and loft failed to complement the TALEN-induced furina mutant pku338, confirming that furina is responsible for the loft mutant phenotypes. Mechanistic studies demonstrated that Notch reporter Tg(tp1:mCherry) signals were largely eliminated in mutant hearts, and overexpression of the Notch intracellular domain partially rescued the mutant phenotypes, probably due to the lack of Furina-mediated cleavage processing of Notch1b proteins, the only Notch receptor expressed in the heart. Together, our data suggest a potential post-translational modification of Notch1b proteins via the proprotein convertase Furina in the heart, and unveil the function of the Furina-Notch1b axis in cardiac looping and trabeculation in zebrafish, and possibly in other organisms.

Zebrafish as a Translational Model: An Experimental Alternative to Study the Mechanisms Involved in Anosmia and Possible Neurodegenerative Aspects of COVID-19?

The Coronavirus disease-2019 (COVID-19) presents a variability of clinical symptoms, ranging from asymptomatic to severe respiratory and systemic conditions. In a cohort of patients, the Severe Acute Respiratory Syndrome Coronavirus (SARS-CoV-2), beyond the classical respiratory manifestations, induces anosmia. Evidence has suggested SARS-CoV-2-induced anosmia can be the result of neurodegeneration of the olfactory pathway. Neurologic symptoms associated with COVID-19 have been reported; however, the precise mechanism and possible long-lasting effects remain poorly investigated. Preclinical models are valuable tools for describing and testing new possible treatments for neurologic disorders. In this way, the zebrafish (Danio rerio) organism model represents an attractive tool in the field of neuroscience, showing economic and logistic advantages besides genetic and physiologic similarities with mammalian, including the brain structure and functions. Besides, its external embryonic development, high availability of eggs, and fast development allows easy genetic manipulation and fast replications. In the present review, we suggest that the zebrafish model can be advantageous to investigate the neurologic features of COVID-19.

Agaricus bisporus-Derived Glucosamine Hydrochloride Facilitates Skeletal Injury Repair through Bmp Signaling in Zebrafish Osteoporosis Model

Glucosamine hydrochloride (GAH), one of the most basic and important derivatives of chitin, is obtained by hydrolysis of chitin in concentrated hydrochloric acid. At present, little is known about how GAH functions in skeletal development. In this report, we demonstrate that GAH, extracted from the cell wall of Agaricus bisporus, acts in a dose-dependent manner to promote not only cartilage and bone development in larvae but also caudal fin regeneration in adult fish. Furthermore, GAH treatment causes a significant increase in expression of bone-related marker genes, indicating its important role in promoting skeletal development. We show that in both larval and adult osteoporosis models induced by high iron osteogenic defects are significantly ameliorated after treatment with GAH, which regulates expression of a series of bone-related genes. Finally, we demonstrate that GAH promotes skeletal development and injury repair through bone morphogenetic protein (Bmp) signaling, and it works at the downstream of the receptor level. Taken together, our findings not only provide a strong research foundation and strategy for the screening of natural osteoporosis drugs and product development using a zebrafish model but also establish the potential for the development of Agaricus bisporus-derived GAH as a new drug for osteoporosis treatment.

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.

Family-effects in the epigenomic response of red blood cells to a challenge test in the European sea bass (Dicentrarchus labrax, L.)

Abstract

Background

In fish, minimally invasive blood sampling is widely used to monitor physiological stress with blood plasma biomarkers. As fish blood cells are nucleated, they might be a source a potential new markers derived from ‘omics technologies. We modified the epiGBS (epiGenotyping By Sequencing) technique to explore changes in genome-wide cytosine methylation in the red blood cells (RBCs) of challenged European sea bass (Dicentrarchus labrax), a species widely studied in both natural and farmed environments.

Results

We retrieved 501,108,033 sequencing reads after trimming, with a mean mapping efficiency of 73.0% (unique best hits). Minor changes in RBC methylome appeared to manifest after the challenge test and a family-effect was detected. Only fifty-seven differentially methylated cytosines (DMCs) close to 51 distinct genes distributed on 17 of 24 linkage groups (LGs) were detected between RBCs of pre- and post-challenge individuals. Thirty-seven of these genes were previously reported as differentially expressed in the brain of zebrafish, most of them involved in stress coping differences. While further investigation remains necessary, few DMC-related genes associated to the Brain Derived Neurotrophic Factor, a protein that favors stress adaptation and fear memory, appear relevant to integrate a centrally produced stress response in RBCs.

Conclusion

Our modified epiGBS protocol was powerful to analyze patterns of cytosine methylation in RBCs of D. labrax and to evaluate the impact of a challenge using minimally invasive blood samples. This study is the first approximation to identify epigenetic biomarkers of exposure to stress in fish.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07420-9.

A familial PLCB4 mutation causing auriculocondylar syndrome 2 with variable severity

Auriculocondylar syndrome (ARCND, MIM #614669, #602483, and #615706); also known as ''question-mark ear syndrome'' or ''dysgnathia complex'', is a rare craniofacial malformation of first and second branchial arches with a prevalence of <1/1,000,000. It is characterized by a distinctive auricular malformation (question mark ear (QME)) and highly variable mandibular anomalies. Variants found in PLCB4, GNAI3, and in EDN1 genes are responsible for >90% of tested ARCND patients. Whole exome sequencing in a multigenerational Egyptian kindred with high intrafamilial variability revealed a known heterozygous missense variant in PLCB4 (NM_000933.3:c.1862G>A:p.(Arg621His)). This report increases the number of molecularly characterized ARCND patients to 29 and emphasizes the highly variable clinical presentation within families.

Zebrafish Pax1a and Pax1b are required for pharyngeal pouch morphogenesis and ceratobranchial cartilage development

Pharyngeal arches are derived from all three germ layers and molecular interactions among the tissue types are required for proper development of subsequent pharyngeal cartilages; however, the mechanisms underlying this process are not fully described. Here we report that in zebrafish, Pax1a and Pax1b have overlapping and essential functions in pharyngeal pouch morphogenesis and subsequent ceratobranchial cartilage development. Both pax1a and pax1b are co-expressed in pharyngeal pouches, and time-lapse imaging of a novel Tg(pax1b:eGFP) enhancer trap line further revealed the sequential segmental development of pharyngeal pouches. Zebrafish pax1a^-/-; pax1b^-/- double mutant embryos generated by CRISPR-Cas9 mutagenesis exhibit unsegmented pharyngeal pouches 2-5 with small outpocketings. Endodermal expression of fgf3, tbx1 and edn1 is also absent in pharyngeal pouches 2-5 at 36 h post fertilization (hpf). Loss of ceratobranchial cartilage 1-4 and reduced or absent expression of dlx2a and hand2 in the pharyngeal arches 3-6 are observed in CRISPR mutant and morphant embryos that are deficient in both zebrafish pax1a and pax1b at 96 or 36 hpf. These results suggest that zebrafish Pax1a and Pax1b both regulate pharyngeal pouch morphogenesis by modulating expression of fgf3 and tbx1. Furthermore, our data support a model wherein endodermal Pax1a and Pax1b act through Fgf3 and Tbx-Edn1 signaling to non-autonomously regulate the development of ceratobranchial cartilage via expression of dlx2a and hand2.

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.

Fibroblast Growth Factor Receptors Function Redundantly During Zebrafish Embryonic Development

Fibroblast growth factor (Fgf) signaling regulates many processes during development. In most cases, one tissue layer secretes an Fgf ligand that binds and activates an Fgf receptor (Fgfr) expressed by a neighboring tissue. Although studies have identified the roles of specific Fgf ligands during development, less is known about the requirements for the receptors. We have generated null mutations in each of the five fgfr genes in zebrafish. Considering the diverse requirements for Fgf signaling throughout development, and that null mutations in the mouse Fgfr1 and Fgfr2 genes are embryonic lethal, it was surprising that all zebrafish homozygous mutants are viable and fertile, with no discernable embryonic defect. Instead, we find that multiple receptors are involved in coordinating most Fgf-dependent developmental processes. For example, mutations in the ligand fgf8a cause loss of the midbrain-hindbrain boundary, whereas, in the fgfr mutants, this phenotype is seen only in embryos that are triple mutant for fgfr1a;fgfr1b;fgfr2, but not in any single or double mutant combinations. We show that this apparent fgfr redundancy is also seen during the development of several other tissues, including posterior mesoderm, pectoral fins, viscerocranium, and neurocranium. These data are an essential step toward defining the specific Fgfrs that function with particular Fgf ligands to regulate important developmental processes in zebrafish.

A genome-wide map of circular RNAs in adult zebrafish

Circular RNAs (circRNAs) are transcript isoforms generated by back-splicing of exons and circularisation of the transcript. Recent genome-wide maps created for circular RNAs in humans and other model organisms have motivated us to explore the repertoire of circular RNAs in zebrafish, a popular model organism. We generated RNA-seq data for five major zebrafish tissues - Blood, Brain, Heart, Gills and Muscle. The repertoire RNA sequence reads left over after reference mapping to linear transcripts were used to identify unique back-spliced exons utilizing a split-mapping algorithm. Our analysis revealed 3,428 novel circRNAs in zebrafish. Further in-depth analysis suggested that majority of the circRNAs were derived from previously well-annotated protein-coding and long noncoding RNA gene loci. In addition, many of the circular RNAs showed extensive tissue specificity. We independently validated a subset of circRNAs using polymerase chain reaction (PCR) and divergent set of primers. Expression analysis using quantitative real time PCR recapitulate selected tissue specificity in the candidates studied. This study provides a comprehensive genome-wide map of circular RNAs in zebrafish tissues.

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.

Abnormalities in gray matter microstructure in young adults with 22q11.2 deletion syndrome

BACKGROUND

22q11.2 Deletion Syndrome (22q11DS) is a genetic, neurodevelopmental disorder characterized by a chromosomal deletion and a distinct cognitive profile. Although abnormalities in the macrostructure of the cortex have been identified in individuals with 22q11DS, it is not known if there are additional microstructural changes in gray matter regions in this syndrome, and/or if such microstructural changes are associated with cognitive functioning.

METHODS

This study employed a novel diffusion MRI measure, the Heterogeneity of Fractional Anisotropy (HFA), to examine variability in the microstructural organization of the cortex in healthy young adults (N = 30) and those with 22q11DS (N = 56). Diffusion MRI, structural MRI, clinical and cognitive data were acquired.

RESULTS

Compared to controls, individuals with 22q11DS evinced increased HFA in cortical association (p = .003, d = 0.86) and paralimbic (p < .0001, d = 1.2) brain areas, whereas no significant differences were found between the two groups in primary cortical brain areas. Additionally, increased HFA of the right paralimbic area was associated with poorer performance on tests of response inhibition, i.e., the Stroop Test (rho = -0.37 p = .005) and the Gordon Diagnostic System Vigilance Commission (rho = -0.41 p = .002) in the 22q11DS group. No significant correlations were found between HFA and cognitive abilities in the healthy control group.

CONCLUSIONS

These findings suggest that cortical microstructural disorganization may be a neural correlate of response inhibition in individuals with 22q11DS. Given that the migration pattern of neural crest cells is disrupted at the time of early brain development in 22q11DS, we hypothesize that these neural alterations may be neurodevelopmental in origin, and reflect cortical dysfunction associated with cognitive deficits.

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.

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

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

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

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

Genome-wide analysis of facial skeletal regionalization in zebrafish

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

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

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

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

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

Building and maintaining joints by exquisite local control of cell fate

We owe the flexibility of our bodies to sophisticated articulations between bones. Establishment of these joints requires the integration of multiple tissue types: permanent cartilage that cushions the articulating bones, synovial membranes that enclose a lubricating fluid-filled cavity, and a fibrous capsule and ligaments that provide structural support. Positioning the prospective joint region involves establishment of an "interzone" region of joint progenitor cells within a nascent cartilage condensation, which is achieved through the interplay of activators and inhibitors of multiple developmental signaling pathways. Within the interzone, tight regulation of BMP and TGFβ signaling prevents the hypertrophic maturation of joint chondrocytes, in part through downstream transcriptional repressors and epigenetic modulators. Synovial cells then acquire further specializations through expression of genes that promote lubrication, as well as the formation of complex structures such as cavities and entheses. Whereas genetic investigations in mice and humans have uncovered a number of regulators of joint development and homeostasis, recent work in zebrafish offers a complementary reductionist approach toward understanding joint positioning and the regulation of chondrocyte fate at joints. The complexity of building and maintaining joints may help explain why there are still few treatments for osteoarthritis, one of the most common diseases in the human population. A major challenge will be to understand how developmental abnormalities in joint structure, as well as postnatal roles for developmental genes in joint homeostasis, contribute to birth defects and degenerative diseases of joints. WIREs Dev Biol 2017, 6:e245. doi: 10.1002/wdev.245 For further resources related to this article, please visit the WIREs website.

Wdr68 Mediates Dorsal and Ventral Patterning Events for Craniofacial Development

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

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.

The exquisite functions of the vertebrate face require the precise formation of its underlying bones. Remarkably, many of the genes required to shape the facial skeleton are the same from fish to man. In this study, we use the powerful zebrafish system to understand how the skeletal components of the face acquire different shapes during development. To do so, we analyze a series of mutants that disrupt patterning of the facial skeleton, and then assess how the genes affected in these mutants control cell fate in skeletal progenitor cells. From these genetic studies, we found that several pathways converge to control when and where progenitor cells commit to a cartilage fate, thus controlling the size and shape of cartilage templates for the later-arising bones. Our work thus reveals how regulating the timing of when progenitor cells make skeleton helps to shape the bones of the zebrafish face. As mutations in many of the genes studied are implicated in human craniofacial defects, differences in the timing of progenitor cell differentiation may also explain the wonderful diversity of human faces.

Establishment of proprotein convertase, furinA knocked-out lines in medaka, Oryzias latipes, and unique form of medaka furin-like prorprotein convertase (mflPC)

Furin is a member of the subtilisin-like proprotein convertase family. Medaka furin-like proprotein convertase (mflPC), a unique form of medaka FurinA (mFurinA) (GenBank accession no. AB092685.1) was cloned from the ovary cDNA library. Compared to human furin (GenBank accession no. NM_002569.3) and mFurinA in the structural motif of mflPC, only the catalytic domain and the N-terminal region of the P domain are highly conserved, but more C-terminal domains are truncated. Based on our research, there three forms of furin, mFurinA, mflPC and mFurinB that exist in medaka. These three genes are expressed in the developing embryos and ubiquitously in adult tissues. To investigate the function of mFurinA and mflPC, as a first step, mFurinA KO lines were established. The mFurinA KO larvae with abnormal phenotypes exhibit edema, abnormal body fluid accumulation in the pericardial and yolk sacs, enlarged hearts, clogged blood vessels, structurally weak eyes, and a very short life. The data suggests that abnormal processing of TGF-β may be one of the causes of these disorders. FurinA KO medaka is a good model for the study of human diseases such as Fraser Syndrome and Marfan syndrome. The creation of human genomic disorder models using recently advanced genome editing procedures informs us of the function of key molecules and their role in causing equivalent human disorders and will be useful as a tool to identify the mechanisms involved.

Zebrafish Craniofacial Development: A Window into Early Patterning

The formation of the face and skull involves a complex series of developmental events mediated by cells derived from the neural crest, endoderm, mesoderm, and ectoderm. Although vertebrates boast an enormous diversity of adult facial morphologies, the fundamental signaling pathways and cellular events that sculpt the nascent craniofacial skeleton in the embryo have proven to be highly conserved from fish to man. The zebrafish Danio rerio, a small freshwater cyprinid fish from eastern India, has served as a popular model of craniofacial development since the 1990s. Unique strengths of the zebrafish model include a simplified skeleton during larval stages, access to rapidly developing embryos for live imaging, and amenability to transgenesis and complex genetics. In this chapter, we describe the anatomy of the zebrafish craniofacial skeleton; its applications as models for the mammalian jaw, middle ear, palate, and cranial sutures; the superior imaging technology available in fish that has provided unprecedented insights into the dynamics of facial morphogenesis; the use of the zebrafish to decipher the genetic underpinnings of craniofacial biology; and finally a glimpse into the most promising future applications of zebrafish craniofacial research.

Sonic hedgehog from pharyngeal arch 1 epithelium is necessary for early mandibular arch cell survival and later cartilage condensation differentiation

BACKGROUND

Morphogenesis of vertebrate craniofacial skeletal elements is dependent on a key cell population, the cranial neural crest cells (NCC). Cranial NCC are formed dorsally in the cranial neural tube and migrate ventrally to form craniofacial skeletal elements as well as other tissues. Multiple extracellular signaling pathways regulate the migration, survival, proliferation, and differentiation of NCC.

RESULTS

In this study, we demonstrate that Shh expression in the oral ectoderm and pharyngeal endoderm is essential for mandibular development. We show that a loss of Shh in these domains results in increased mesenchymal cell death in pharyngeal arch 1 (PA1) after NCC migration. This increased cell death can be rescued in utero by pharmacological inhibition of p53. Furthermore, we show that epithelial SHH is necessary for the early differentiation of mandibular cartilage condensations and, therefore, the subsequent development of Meckel's cartilage, around which the dentary bone forms. Nonetheless, a rescue of the cell death phenotype does not rescue the defect in cartilage condensation formation.

CONCLUSIONS

Our results show that SHH produced by the PA1 epithelium is necessary for the survival of post-migratory NCC, and suggests a key role in the subsequent differentiation of chondrocytes to form Meckel's cartilage.

Cre recombinase-regulated Endothelin1 transgenic mouse lines: novel tools for analysis of embryonic and adult disorders

Endothelin-1 (EDN1) influences both craniofacial and cardiovascular development and a number of adult physiological conditions by binding to one or both of the known endothelin receptors, thus initiating multiple signaling cascades. Animal models containing both conventional and conditional loss of the Edn1 gene have been used to dissect EDN1 function in both embryos and adults. However, while transgenic Edn1 over-expression or targeted genomic insertion of Edn1 has been performed to understand how elevated levels of Edn1 result in or exacerbate disease states, an animal model in which Edn1 over-expression can be achieved in a spatiotemporal-specific manner has not been reported. Here we describe the creation of Edn1 conditional over-expression transgenic mouse lines in which the chicken β-actin promoter and an Edn1 cDNA are separated by a strong stop sequence flanked by loxP sites. In the presence of Cre, the stop cassette is removed, leading to Edn1 expression. Using the Wnt1-Cre strain, in which Cre expression is targeted to the Wnt1-expressing domain of the central nervous system (CNS) from which neural crest cells (NCCs) arise, we show that stable chicken β-actin-Edn1 (CBA-Edn1) transgenic lines with varying EDN1 protein levels develop defects in NCC-derived tissues of the face, though the severity differs between lines. We also show that Edn1 expression can be achieved in other embryonic tissues utilizing other Cre strains, with this expression also resulting in developmental defects. CBA-Edn1 transgenic mice will be useful in investigating diverse aspects of EDN1-mediated-development and disease, including understanding how NCCs achieve and maintain a positional and functional identity and how aberrant EDN1 levels can lead to multiple physiological changes and diseases.

Nodal signaling range is regulated by proprotein convertase-mediated maturation

Tissue patterning is established by extracellular growth factors or morphogens. Although different theoretical models explaining specific patterns have been proposed, our understanding of tissue pattern establishment in vivo remains limited. In many animal species, left-right patterning is governed by a reaction-diffusion system relying on the different diffusivity of an activator, Nodal, and an inhibitor, Lefty. In a genetic screen, we identified a zebrafish loss-of-function mutant for the proprotein convertase FurinA. Embryological and biochemical experiments demonstrate that cleavage of the Nodal-related Spaw proprotein into a mature form by FurinA is required for Spaw gradient formation and activation of Nodal signaling. We demonstrate that FurinA is required cell-autonomously for the long-range signaling activity of Spaw and no other Nodal-related factors. Combined in silico and in vivo approaches support a model in which FurinA controls the signaling range of Spaw by cleaving its proprotein into a mature, extracellular form, consequently regulating left-right patterning.

The proprotein convertase subtilisin/kexin furinA regulates zebrafish host response against Mycobacterium marinum

Tuberculosis is a chronic bacterial disease with a complex pathogenesis. An effective immunity against Mycobacterium tuberculosis requires both the innate and adaptive immune responses, including proper T helper (Th) type 1 cell function. FURIN is a proprotein convertase subtilisin/kexin (PCSK) enzyme, which is highly expressed in Th1 type cells. FURIN expression in T cells is essential for maintaining peripheral immune tolerance, but its role in the innate immunity and infections has remained elusive. Here, we utilized Mycobacterium marinum infection models in zebrafish (Danio rerio) to investigate how furin regulates host responses against mycobacteria. In steady-state furinAtd204e/+ fish reduced furinA mRNA levels associated with low granulocyte counts and elevated Th cell transcription factor expressions. Silencing furin genes reduced the survival of M. marinum-infected zebrafish embryos. A mycobacterial infection upregulated furinA in adult zebrafish, and infected furinAtd204e/+ mutants exhibited a proinflammatory phenotype characterized by elevated tumor necrosis factor a (tnfa), lymphotoxin alpha (lta) and interleukin 17a/f3 (il17a/f3) expression levels. The enhanced innate immune response in the furinAtd204e/+ mutants correlated with a significantly decreased bacterial burden in a chronic M. marinum infection model. Our data show that upregulated furinA expression can serve as a marker for mycobacterial disease, since it inhibits early host responses and consequently promotes bacterial growth in a chronic infection.

Mmp25β facilitates elongation of sensory neurons during zebrafish development

Matrix metalloproteinases (MMPs) are a large and complex family of zinc-dependent endoproteinases widely recognized for their roles in remodeling the extracellular matrix (ECM) during embryonic development, wound healing, and tissue homeostasis. Their misregulation is central to many pathologies, and they have therefore been the focus of biomedical research for decades. These proteases have also recently emerged as mediators of neural development and synaptic plasticity in vertebrates, however, understanding of the mechanistic basis of these roles and the molecular identities of the MMPs involved remains far from complete. We have identified a zebrafish orthologue of mmp25 (a.k.a. leukolysin; MT6-MMP), a membrane-type, furin-activated MMP associated with leukocytes and invasive carcinomas, but which we find is expressed by a subset of the sensory neurons during normal embryonic development. We detect high levels of Mmp25β expression in the trigeminal, craniofacial, and posterior lateral line ganglia in the hindbrain, and in Rohon-Beard cells in the dorsal neural tube during the first 48 h of embryonic development. Knockdown of Mmp25β expression with morpholino oligonucleotides results in larvae that are uncoordinated and insensitive to touch, and which exhibit defects in the development of sensory neural structures. Using in vivo zymography, we observe that Mmp25β morphant embryos show reduced Type IV collagen degradation in regions of the head traversed by elongating axons emanating from the trigeminal ganglion, suggesting that Mmp25β may play a pivotal role in mediating ECM remodeling in the vicinity of these elongating axons.

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.