An early requirement for maternal FoxH1 during zebrafish gastrulation

Conservation of cis-Regulatory Syntax Underlying Deuterostome Gastrulation

Throughout embryonic development, the shaping of the functional and morphological characteristics of embryos is orchestrated by an intricate interaction between transcription factors and cis-regulatory elements. In this study, we conducted a comprehensive analysis of deuterostome cis-regulatory landscapes during gastrulation, focusing on four paradigmatic species: the echinoderm Strongylocentrotus purpuratus, the cephalochordate Branchiostoma lanceolatum, the urochordate Ciona intestinalis, and the vertebrate Danio rerio. Our approach involved comparative computational analysis of ATAC-seq datasets to explore the genome-wide blueprint of conserved transcription factor binding motifs underlying gastrulation. We identified a core set of conserved DNA binding motifs associated with 62 known transcription factors, indicating the remarkable conservation of the gastrulation regulatory landscape across deuterostomes. Our findings offer valuable insights into the evolutionary molecular dynamics of embryonic development, shedding light on conserved regulatory subprograms and providing a comprehensive perspective on the conservation and divergence of gene regulation underlying the gastrulation process.

Maternal vgll4a regulates zebrafish epiboly through Yap1 activity

Gastrulation in zebrafish embryos commences with the morphogenetic rearrangement of blastodermal cells, which undergo a coordinated spreading from the animal pole to wrap around the egg at the vegetal pole. This rearrangement, known as epiboly, relies on the orchestrated activity of maternal transcripts present in the egg, compensating for the gradual activation of the zygotic genome. Epiboly involves the mechano-transducer activity of yap1 but what are the regulators of yap1 activity and whether these are maternally or zygotically derived remain elusive. Our study reveals the crucial role of maternal vgll4a, a proposed Yap1 competitor, during zebrafish epiboly. In embryos lacking maternal/zygotic vgll4a (MZvgll4a), the progression of epiboly and blastopore closure is delayed. This delay is associated with the ruffled appearance of the sliding epithelial cells, decreased expression of yap1-downstream targets and transient impairment of the actomyosin ring at the syncytial layer. Our study also shows that, rather than competing with yap1, vgll4a modulates the levels of the E-cadherin/β-catenin adhesion complex at the blastomeres' plasma membrane and hence their actin cortex distribution. Taking these results together, we propose that maternal vgll4a acts at epiboly initiation upstream of yap1 and the E-cadherin/β-catenin adhesion complex, contributing to a proper balance between tissue tension/cohesion and contractility, thereby promoting a timely epiboly progression.

MCPIP1 functions as a safeguard of early embryonic development

Monocyte chemoattractant protein-induced protein 1 (MCPIP1), also called Regnase-1, is an RNase that has been described as a key negative modulator of inflammation. MCPIP1 also controls numerous tumor-related processes, such as proliferation, apoptosis and differentiation. In this study, we utilized a zebrafish model to investigate the role of Mcpip1 during embryogenic development. Our results demonstrated that during embryogenesis, the expression of the zc3h12a gene encoding Mcpip1 undergoes dynamic changes. Its transcript levels gradually increase from the 2-cell stage to the spherical stage and then decrease rapidly. We further found that ectopic overexpression of wild-type Mcpip1 but not the catalytically inactive mutant form resulted in an embryonic lethal phenotype in zebrafish embryos (24 hpf). At the molecular level, transcriptomic profiling revealed extensive changes in the expression of genes encoding proteins important in the endoplasmic reticulum stress response and in protein folding as well as involved in the formation of primary germ layer, mesendoderm and endoderm development, heart morphogenesis and cell migration. Altogether, our results demonstrate that the expression of zc3h12a must be tightly controlled during the first cell divisions of zebrafish embryos and that a rapid decrease in its mRNA expression is an important factor promoting proper embryo development.

FoxH1 Represses the Promoter Activity of cyp19a1a in the Ricefield Eel (Monopterus albus)

Forkhead box H1 (FoxH1) is a sexually dimorphic gene in Oreochromis niloticus, Oplegnathus fasciatus, and Acanthopagrus latus, indicating that it is essential for gonadal development. In the present study, the molecular characteristics and potential function of FoxH1 and the activation of the cyp19a1a promoter in vitro were evaluated in Monopterus albus. The levels of foxh1 in the ovaries were three times higher than those in the testes and were regulated by gonadotropins (Follicle-Stimulating Hormone and Human Chorionic Gonadotropin). FoxH1 colocalized with Cyp19a1a in the oocytes and granulosa cells of middle and late vitellogenic follicles. In addition, three FoxH1 binding sites were identified in the proximal promoter of cyp19a1a, namely, FH1 (-871/-860), FH2 (-535/-524), and FH3 (-218/-207). FoxH1 overexpression significantly attenuated the activity of the cyp19a1a promoter in CHO cells, and FH1/2 mutation increased promoter activity. Taken together, these results suggest that FoxH1 may act as an important regulator in the ovarian development of M. albus by repressing cyp19a1a promoter activity, which provides a foundation for the study of FoxH1 function in bony fish reproductive processes.

Intergenerational plasticity to cycling high temperature and hypoxia affects offspring stress responsiveness and tolerance in zebrafish

Predicted climate change-induced increases in heat waves and hypoxic events will have profound effects on fishes, yet the capacity of parents to alter offspring phenotype via non-genetic inheritance and buffer against these combined stressors is not clear. This study tested how prolonged adult zebrafish exposure to combined diel cycles of thermal stress and hypoxia affect offspring early survival and development, parental investment of cortisol and heat shock proteins (HSPs), larval offspring stress responses, and both parental and offspring heat and hypoxia tolerance. Parental exposure to the combined stressor did not affect fecundity, but increased mortality, produced smaller embryos and delayed hatching. The combined treatment also reduced maternal deposition of cortisol and increased embryo hsf1, hsp70a, HSP70, hsp90aa and HSP90 levels. In larvae, basal cortisol levels did not differ between treatments, but acute exposure to combined heat stress and hypoxia increased cortisol levels in control larvae with no effect on larvae from exposed parents. In contrast, whereas larval basal hsf1, hsp70a and hsp90aa levels differed between parental treatments, the combined acute stressor elicited similar transcriptional responses across treatments. Moreover, the combined acute stressor only induced a marked increase in HSP47 levels in the larvae derived from exposed parents. Finally, combined hypoxia and elevated temperatures increased both thermal and hypoxia tolerance in adults and conferred an increase in offspring thermal but not hypoxia tolerance. These results demonstrate that intergenerational acclimation to combined thermal stress and hypoxia elicit complex carryover effects on stress responsiveness and offspring tolerance with potential consequences for resilience.

Transcriptomic and metabolomic analyses revealed epiboly delayed mechanisms of 2,5-dichloro-1, 4-benuinone on zebrafish embryos

2,5-Dichloro-1,4-benzenediol (2,5-DCBQ) is a putative disinfection by-product that belongs to the halogenated benzoquinone class. However, its developmental toxicity and related mechanism remained unclarified. In our study, we used zebrafish embryos as the model and exposed them to graded concentrations of 2,5-DCBQ (100, 200, 300, 400 μg/L). We found that the rate of epiboly abnormalities increased significantly in a concentration-dependent manner. The results of whole-mount in situ hybridization (WISH) indicated that the expression patterns and levels of chordin (dorsoventral marker), foxa2 (endodermal marker), eve1 (ventral mesodermal marker), and foxb1a (ectodermal marker) were altered, suggesting that 2,5-DCBQ might affect the germ layer development of zebrafish embryos. Integrated transcriptomic and metabolomic analyses were adopted to explore the molecular mechanisms of embryonic developmental delays. The results showed that 2,5-DCBQ exposure induced 1163 differentially expressed genes (DEGs) and 37 differential metabolites (DEMs). Bioinformatic analysis enriched the most affected molecular pathways (Wnt signaling pathway, cell adhesion molecules, actin cytoskeleton regulation) and metabolic pathways (purine metabolism, aminoacyl-tRNA biosynthesis, arginine and proline metabolism) in zebrafish embryos. To summarize, our findings broadened the molecular mechanisms of 2,5-DCBQ embryotoxicity through multi-omics and bioinformatic analyses.

Origin, form and function of extraembryonic structures in teleost fishes

Teleost eggs have evolved a highly derived early developmental pattern within vertebrates as a result of the meroblastic cleavage pattern, giving rise to a polar stratified architecture containing a large acellular yolk and a small cellular blastoderm on top. Besides the acellular yolk, the teleost-specific yolk syncytial layer (YSL) and the superficial epithelial enveloping layer are recognized as extraembryonic structures that play critical roles throughout embryonic development. They provide enriched microenvironments in which molecular feedback loops, cellular interactions and mechanical signals emerge to sculpt, among other things, embryonic patterning along the dorsoventral and left-right axes, mesendodermal specification and the execution of morphogenetic movements in the early embryo and during organogenesis. An emerging concept points to a critical role of extraembryonic structures in reinforcing early genetic and morphogenetic programmes in reciprocal coordination with the embryonic blastoderm, providing the necessary boundary conditions for development to proceed. In addition, the role of the enveloping cell layer in providing mechanical, osmotic and immunological protection during early stages of development, and the autonomous nutritional support provided by the yolk and YSL, have probably been key aspects that have enabled the massive radiation of teleosts to colonize every ecological niche on the Earth. This article is part of the theme issue 'Extraembryonic tissues: exploring concepts, definitions and functions across the animal kingdom'.

Rab5ab-Mediated Yolk Cell Membrane Endocytosis Is Essential for Zebrafish Epiboly and Mechanical Equilibrium During Gastrulation

Morphogenesis in early embryos demands the coordinated distribution of cells and tissues to their final destination in a spatio-temporal controlled way. Spatial and scalar differences in adhesion and contractility are essential for these morphogenetic movements, while the role that membrane remodeling may play remains less clear. To evaluate how membrane turnover modulates tissue arrangements we studied the role of endocytosis in zebrafish epiboly. Experimental analyses and modeling have shown that the expansion of the blastoderm relies on an asymmetry of mechanical tension in the yolk cell generated as a result of actomyosin-dependent contraction and membrane removal. Here we show that the GTPase Rab5ab is essential for the endocytosis and the removal of the external yolk cell syncytial layer (E-YSL) membrane. Interfering in its expression exclusively in the yolk resulted in the reduction of yolk cell actomyosin contractility, the disruption of cortical and internal flows, a disequilibrium in force balance and epiboly impairment. We conclude that regulated membrane remodeling is crucial for directing cell and tissue mechanics, preserving embryo geometry and coordinating morphogenetic movements during epiboly.

Homozygous mutation of foxh1 arrests oogenesis causing infertility in female Nile tilapia†

Foxh1, a member of fox gene family, was first characterized as a transcriptional partner in the formation of the Smad protein complex. Recent studies have shown foxh1 is highly expressed in the cytoplasm of oocytes in both tilapia and mouse. However, its function in oogenesis remains unexplored. In the present study, foxh1-/- tilapia was created by CRISPR/Cas9. At 180 dah (days after hatching), the foxh1-/- XX fish showed oogenesis arrest and a significantly lower GSI. The transition of oocytes from phase II to phase III and follicle cells from one to two layers was blocked, resulting in infertility of the mutant. Transcriptomic analysis revealed that expression of genes involved in estrogen synthesis and oocyte growth were altered in the foxh1-/- ovaries. Loss of foxh1 resulted in significantly decreased Cyp19a1a and increased Cyp11b2 expression, consistent with significantly lower concentrations of serum estradiol-17β (E2) and higher concentrations of 11-ketotestosterone (11-KT). Moreover, administration of E2 rescued the phenotypes of foxh1-/- XX fish, as indicated by the appearance of phase III and IV oocytes and absence of Cyp11b2 expression. Taken together, these results suggest that foxh1 functions in the oocytes to regulate oogenesis by promoting cyp19a1a expression, and therefore estrogen production. Disruption of foxh1 may block the estrogen synthesis and oocyte growth, leading to the arrest of oogenesis and thus infertility in tilapia.

Assessing the effects of an acute exposure to worst-case concentration of Cry proteins on zebrafish using the embryotoxicity test and proteomics analysis

Cry1C, Cry1F and Cry1Ab are insecticidal proteins from Bacillus thuringiensis (Bt) which are expressed in transgenic crops. Given the entry of these proteins into aquatic environments, it is relevant to evaluate their impacts on aquatic organisms. In this work, we sought to evaluate the effects of Cry1C, Cry1F and Cry1Ab on zebrafish embryos and larvae of a predicted worst-case scenario concentration of these proteins (set to 1.1 mg/L). For that, we coupled a traditional toxicity approach (the zebrafish embryotoxicity test and dosage of enzymatic biomarkers) to gel free proteomics analysis. At the concentration tested, these proteins did not cause adverse effects in the zebrafish early life stages, either by verifying phenotypic endpoints of toxicity or alterations in representative enzymatic biomarkers (catalase, glutathione-S-tranferase and lactate-dehydrogenase). At the molecular level, the Cry proteins tested lead to very small changes in the proteome of zebrafish larvae. In a global way, these proteins upregulated the expression of vitellogenins. Besides that, Cry1C e Cry1F deregulated heterogeneous nuclear ribonucleoproteins (Hnrnpa0l and Hnrnpaba, respectively), implicated in mRNA processing and gene regulation. Overall, these data lead to the conclusion that Cry1C, Cry1F and Cry1Ab proteins, even at a very high concentration, have limited effects in the early stages of zebrafish life.

Analysis of zebrafish periderm enhancers facilitates identification of a regulatory variant near human KRT8/18

Genome-wide association studies for non-syndromic orofacial clefting (OFC) have identified single nucleotide polymorphisms (SNPs) at loci where the presumed risk-relevant gene is expressed in oral periderm. The functional subsets of such SNPs are difficult to predict because the sequence underpinnings of periderm enhancers are unknown. We applied ATAC-seq to models of human palate periderm, including zebrafish periderm, mouse embryonic palate epithelia, and a human oral epithelium cell line, and to complementary mesenchymal cell types. We identified sets of enhancers specific to the epithelial cells and trained gapped-kmer support-vector-machine classifiers on these sets. We used the classifiers to predict the effects of 14 OFC-associated SNPs at 12q13 near KRT18. All the classifiers picked the same SNP as having the strongest effect, but the significance was highest with the classifier trained on zebrafish periderm. Reporter and deletion analyses support this SNP as lying within a periderm enhancer regulating KRT18/KRT8 expression.

Setting up for gastrulation in zebrafish

Soon after fertilization the zebrafish embryo generates the pool of cells that will give rise to the germline and the three somatic germ layers of the embryo (ectoderm, mesoderm and endoderm). As the basic body plan of the vertebrate embryo emerges, evolutionarily conserved developmental signaling pathways, including Bmp, Nodal, Wnt, and Fgf, direct the nearly totipotent cells of the early embryo to adopt gene expression profiles and patterns of cell behavior specific to their eventual fates. Several decades of molecular genetics research in zebrafish has yielded significant insight into the maternal and zygotic contributions and mechanisms that pattern this vertebrate embryo. This new understanding is the product of advances in genetic manipulations and imaging technologies that have allowed the field to probe the cellular, molecular and biophysical aspects underlying early patterning. The current state of the field indicates that patterning is governed by the integration of key signaling pathways and physical interactions between cells, rather than a patterning system in which distinct pathways are deployed to specify a particular cell fate. This chapter focuses on recent advances in our understanding of the genetic and molecular control of the events that impart cell identity and initiate the patterning of tissues that are prerequisites for or concurrent with movements of gastrulation.

Mechanisms of zebrafish epiboly: A current view

Epiboly is a conserved gastrulation movement describing the thinning and spreading of a sheet or multi-layer of cells. The zebrafish embryo has emerged as a vital model system to address the cellular and molecular mechanisms that drive epiboly. In the zebrafish embryo, the blastoderm, consisting of a simple squamous epithelium (the enveloping layer) and an underlying mass of deep cells, as well as a yolk nuclear syncytium (the yolk syncytial layer) undergo epiboly to internalize the yolk cell during gastrulation. The major events during zebrafish epiboly are: expansion of the enveloping layer and the internal yolk syncytial layer, reduction and removal of the yolk membrane ahead of the advancing blastoderm margin and deep cell rearrangements between the enveloping layer and yolk syncytial layer to thin the blastoderm. Here, work addressing the cellular and molecular mechanisms as well as the sources of the mechanical forces that underlie these events is reviewed. The contribution of recent findings to the current model of epiboly as well as open questions and future prospects are also discussed.

Tris(1,3-dichloro-2-propyl) Phosphate Exposure During the Early-Blastula Stage Alters the Normal Trajectory of Zebrafish Embryogenesis

Tris(1,3-dichloro-2-propyl) phosphate (TDCIPP) is an organophosphate flame retardant used around the world. Within zebrafish, we previously showed that initiation of TDCIPP exposure during cleavage (0.75 h post-fertilization, hpf) results in epiboly disruption at 6 hpf, leading to dorsalized embryos by 24 hpf, a phenotype that mimics the effects of dorsomorphin (DMP), a bone morphogenetic protein (BMP) antagonist that dorsalizes embryos in the absence of epiboly defects. The objective of this study was to (1) investigate the role of BMP signaling in TDCIPP-induced toxicity during early embryogenesis, (2) identify other pathways and processes targeted by TDCIPP, and (3) characterize the downstream impacts of early developmental defects. Using zebrafish as a model, we first identified a sensitive window for TDCIPP-induced effects following exposure initiation at 0.75 hpf. We then investigated the effects of TDCIPP on the transcriptome during the first 24 h of development using mRNA sequencing and amplicon sequencing. Finally, we relied on whole-mount immunohistochemistry, dye-based labeling, and morphological assessments to study abnormalities later in embryonic development. Overall, our data suggest that the initiation of TDCIPP exposure during early blastula alters the normal trajectory of early embryogenesis by inducing gastrulation defects and aberrant germ-layer formation, leading to abnormal tissue and organ development within the embryo.

Contractility, differential tension and membrane removal lead zebrafish epiboly biomechanics

Precise tissue remodeling during development is essential for shaping embryos and optimal organ function. Epiboly is an early gastrulation event by which the blastoderm expands around the yolk to engulf it. Three different layers are involved in this process, an epithelial layer (the enveloping layer, EVL), the embryo proper, constituted by the deep cells (DCs), and the yolk cell. Although teleost epiboly has been studied for many years, a clear understanding of its mechanics was still missing. Here we present new information on the cellular, molecular and mechanical elements involved in epiboly that, together with some other recent data and upon comparison with previous biomechanical models, lets conclude that the expansion of the epithelia is passive and driven by active cortical contraction and membrane removal in the adjacent layer, the External Yolk Syncytial Layer (E-YSL). The isotropic actomyosin contraction of the E-YSL cortex generates an anisotropic stress pattern and a directional net movement consequence of the differences in the deformation response of the 2 opposites adjacent domains (EVL and the Yolk Cytoplasmic Layer - YCL). Contractility is accompanied by the local formation of membrane folds and its removal by Rab5ab dependent macropinocytosis. The increase in area of the epithelia during the expansion is achieved by cell-shape changes (flattening) responding to spherical geometrical cues. The counterbalance between the geometry of the embryo and forces dissipation among different elements is therefore essential for epiboly global coordination.

Competition between histone and transcription factor binding regulates the onset of transcription in zebrafish embryos

Upon fertilization, the genome of animal embryos remains transcriptionally inactive until the maternal-to-zygotic transition. At this time, the embryo takes control of its development and transcription begins. How the onset of zygotic transcription is regulated remains unclear. Here, we show that a dynamic competition for DNA binding between nucleosome-forming histones and transcription factors regulates zebrafish genome activation. Taking a quantitative approach, we found that the concentration of non-DNA-bound core histones sets the time for the onset of transcription. The reduction in nuclear histone concentration that coincides with genome activation does not affect nucleosome density on DNA, but allows transcription factors to compete successfully for DNA binding. In agreement with this, transcription factor binding is sensitive to histone levels and the concentration of transcription factors also affects the time of transcription. Our results demonstrate that the relative levels of histones and transcription factors regulate the onset of transcription in the embryo.

DOI:http://dx.doi.org/10.7554/eLife.23326.001

The DNA in a fertilized egg contains all the information required to form an animal’s body. In order for the animal to develop properly, particular genes encoded in the DNA are only active at specific times. The DNA is wrapped around proteins called histones, which allows the DNA to be tightly packed inside the cell. However, histones can block other proteins called transcription factors from binding to the DNA to activate the genes. Young embryos initially develop with all of their genes switched off, relying on the nutrients and other molecules provided by their mother. After some time, the embryo starts to switch on its own genes to take control of its own development, but it was not clear how this happens.

Joseph et al. investigated how genes are activated in zebrafish embryos, which are often used as models to study how animals develop. The experiments show that competition between histones and transcription factors for binding to DNA controls when genes are switched on. In young fish embryos, there are so many histones present that transcription factors have no opportunity to bind to DNA. Over time, however, the numbers of histones decrease, allowing transcription factors to bind to DNA and switch on genes.

Histones and transcription factors regulate the activity of genes throughout the life of the animal. Therefore, competition between these two types of protein may also control gene activity in other situations. A better understanding of how gene activity is controlled could allow researchers to more easily grow different types of cell in the laboratory or to reprogram specific cells in the body. As such, these new findings may aid the development of therapies to regenerate organs or tissues that have been damaged by injury or disease.

DOI:http://dx.doi.org/10.7554/eLife.23326.002

Expression and distribution of forkhead activin signal transducer 2 (FAST2) during follicle development in mouse ovaries and pre-implantation embryos

Xenopus forkhead activin signal transducer 1 (xFAST 1) was first characterized in Xenopus as the transcriptional partner for Smad proteins. FAST2, which is the xFAST 1 homologues in mouse, is expressed during early developmental stages of the organism. However, the function of FAST2 in mouse ovaries and pre-implantation embryos is unclear. Therefore, we investigated its expression during these processes. In postnatal mice, FAST2 was expressed in oocytes and thecal cells from postnatal day (PD) 1 to PD 21. In gonadotropin-induced immature mice, FAST2 was expressed in oocytes, thecal cells and newly formed corpora lutea (CLs), but was expressed at a lower level in degenerated CLs. Similar results were observed upon western blot analyses. In meloxicam-treated immature mice, ovulation was inhibited and FAST2 was expressed in thecal cells, luteinized granulosa cells and entrapped oocytes. Immunofluorescence results showed that FAST2 was expressed in the cytoplasm and nucleus but not the nucleolus from the zygote to 8-cell embryo stage, after which it was localized to the cytoplasm of the morulae and inner cell mass of the blastocysts. Taken together, these observations suggest that FAST2 is expressed in a cell-specific manner during ovarian follicle development, ovulation, luteinization and early embryonic development, and that FAST2 might play important roles in these physiological processes.

FoxH1 mediates a Grg4 and Smad2 dependent transcriptional switch in Nodal signaling during Xenopus mesoderm development

In the vertebrate blastula and gastrula the Nodal pathway is essential for formation of the primary germ layers and the organizer. Nodal autoregulatory feedback potentiates signaling activity, but mechanisms limiting embryonic Nodal ligand transcription are poorly understood. Here we describe a transcriptional switch mechanism mediated by FoxH1, the principle effector of Nodal autoregulation. FoxH1 contains a conserved engrailed homology (EH1) motif that mediates direct binding of groucho-related gene 4 (Grg4), a Groucho family corepressor. Nodal-dependent gene expression is suppressed by FoxH1, but enhanced by a FoxH1 EH1 mutant, indicating that the EH1 motif is necessary for repression. Grg4 blocks Nodal-induced mesodermal gene expression and Nodal autoregulation, suggesting that Grg4 limits Nodal pathway activity. Conversely, blocking Grg4 function in the ectoderm results in ectopic expression of Nodal target genes. FoxH1 and Grg4 occupy the Xnr1 enhancer, and Grg4 occupancy is dependent on the FoxH1 EH1 motif. Grg4 occupancy at the Xnr1 enhancer significantly decreases with Nodal activation or Smad2 overexpression, while FoxH1 occupancy is unaffected. These results suggest that Nodal-activated Smad2 physically displaces Grg4 from FoxH1, an essential feature of the transcriptional switch mechanism. In support of this model, when FoxH1 is unable to bind Smad2, Grg4 occupancy is maintained at the Xnr1 enhancer, even in the presence of Nodal signaling. Our findings reveal that FoxH1 mediates both activation and repression of Nodal gene expression. We propose that this transcriptional switch is essential to delimit Nodal pathway activity in vertebrate germ layer formation.

PIAS-like protein Zimp7 is required for the restriction of the zebrafish organizer and mesoderm development

The Zmiz2 (Zimp7) protein and its homolog Zmiz1 (Zimp10) were initially identified in humans as androgen receptor co-activators. Sequence analysis revealed the presence of an SP-RING/Miz domain, which is highly conserved in members of the PIAS family and confers SUMO-conjugating activity. Zimp7 has been shown to interact with components of the Wnt/β-Catenin signaling pathway and with Brg1 and BAF57, components of the ATP-dependent mammalian SWI/SNF-like BAF chromatin-remodeling complexes. In this work, we analyze the role of zygotic Zimp7 in zebrafish development. We describe evidence indicating that Zimp7 is required for mesoderm development and dorsoventral patterning. Morpholino-mediated reduction of zygotic Zimp7 produced axial mesodermal defects that were preceded by up-regulation of organizer genes such as bozozok, goosecoid and floating head at the onset of gastrulation and by down-regulation of the ventral markers vox, vent and eve1 indicating loss of the ventrolateral mesoderm. Consistently, embryos overexpressing zimp7 RNA exhibited midline defects such as loss of forebrain and cyclopia accompanied by transcriptional changes directly opposite of those found in the morphants. In addition, the patterning of ventralized embryos produced by the overexpression of vox and vent was restored by a reduction of Zimp7 activity. Altogether, our findings indicate that Zimp7 is involved in transcriptional regulation of factors that are essential for patterning in the dorsoventral axis.

Zebrafish Rnf111 is encoded by multiple transcripts and is required for epiboly progression and prechordal plate development

Arkadia (also known as RING finger 111) encodes a nuclear E3 ubiquitin ligase that targets intracellular effectors and modulators of TGFβ/Nodal-related signaling for polyubiquitination and proteasome-dependent degradation. In the mouse, loss of Arkadia results in early embryonic lethality, with defects attributed to compromised Nodal signaling. Here, we report the isolation of zebrafish arkadia/rnf111, which is represented by 5 transcript variants. arkadia/rnf111 is broadly expressed during the blastula and gastrula stages, with eventual enrichment in the anterior mesendoderm, including the prechordal plate. Morpholino knockdown experiments reveal an unexpected role for Arkadia/Rnf111 in both early blastula organization and epiboly progression. Using a splice junction morpholino, we present additional evidence that arkadia/rnf111 transcript variants containing a 3' alternative exon are specifically required for epiboly progression in the late gastrula. This result suggests that arkadia/rnf111 transcript variants encode functionally relevant protein isoforms that provide additional intracellular flexibility and regulation to the Nodal signaling pathway.

Genome-wide view of TGFβ/Foxh1 regulation of the early mesendoderm program

Nodal/TGFβ signaling regulates diverse biological responses. By combining RNA-seq on Foxh1 and Nodal signaling loss-of-function embryos with ChIP-seq of Foxh1 and Smad2/3, we report a comprehensive genome-wide interaction between Foxh1 and Smad2/3 in mediating Nodal signaling during vertebrate mesendoderm development. This study significantly increases the total number of Nodal target genes regulated by Foxh1 and Smad2/3, and reinforces the notion that Foxh1-Smad2/3-mediated Nodal signaling directly coordinates the expression of a cohort of genes involved in the control of gene transcription, signaling pathway modulation and tissue morphogenesis during gastrulation. We also show that Foxh1 may function independently of Nodal signaling, in addition to its role as a transcription factor mediating Nodal signaling via Smad2/3. Finally, we propose an evolutionarily conserved interaction between Foxh1 and PouV, a mechanism observed in Pou5f1-mediated regulation of pluripotency in human embryonic stem and epiblast cells.

Global identification of Smad2 and Eomesodermin targets in zebrafish identifies a conserved transcriptional network in mesendoderm and a novel role for Eomesodermin in repression of ectodermal gene expression

Background

Nodal signalling is an absolute requirement for normal mesoderm and endoderm formation in vertebrate embryos, yet the transcriptional networks acting directly downstream of Nodal and the extent to which they are conserved is largely unexplored, particularly in vivo. Eomesodermin also plays a role in patterning mesoderm and endoderm in vertebrates, but its mechanisms of action and how it interacts with the Nodal signalling pathway are still unclear.

Results

Using a combination of expression analysis and chromatin immunoprecipitation with deep sequencing (ChIP-seq) we identify direct targets of Smad2, the effector of Nodal signalling in blastula stage zebrafish embryos, including many novel target genes. Through comparison of these data with published ChIP-seq data in human, mouse and Xenopus we show that the transcriptional network driven by Smad2 in mesoderm and endoderm is conserved in these vertebrate species. We also show that Smad2 and zebrafish Eomesodermin a (Eomesa) bind common genomic regions proximal to genes involved in mesoderm and endoderm formation, suggesting Eomesa forms a general component of the Smad2 signalling complex in zebrafish. Combinatorial perturbation of Eomesa and Smad2-interacting factor Foxh1 results in loss of both mesoderm and endoderm markers, confirming the role of Eomesa in endoderm formation and its functional interaction with Foxh1 for correct Nodal signalling. Finally, we uncover a novel role for Eomesa in repressing ectodermal genes in the early blastula.

Conclusions

Our data demonstrate that evolutionarily conserved developmental functions of Nodal signalling occur through maintenance of the transcriptional network directed by Smad2. This network is modulated by Eomesa in zebrafish which acts to promote mesoderm and endoderm formation in combination with Nodal signalling, whilst Eomesa also opposes ectoderm gene expression. Eomesa, therefore, regulates the formation of all three germ layers in the early zebrafish embryo.

Electronic supplementary material

The online version of this article (doi:10.1186/s12915-014-0081-5) contains supplementary material, which is available to authorized users.

Microinjection manipulations in the elucidation of Xenopus brain development

Microinjection has a long and distinguished history in Xenopus and has been used to introduce a surprisingly diverse array of agents into embryos by both intra- and intercellular means. In addition to nuclei, investigators have variously injected peptides, antibodies, biologically active chemicals, lineage markers, mRNA, DNA, morpholinos, and enzymes. While enumerating many of the different microinjection approaches that can be taken, we will focus upon the mechanical operations and options available to introduce mRNA, DNA, and morpholinos intracellularly into early stage embryos for the study of neurogenesis.

Dynamin-dependent maintenance of epithelial integrity is essential for zebrafish epiboly

Epiboly, the thinning and spreading of one tissue over another, is a widely employed morphogenetic movement that is essential for the development of many organisms. In the zebrafish embryo, epiboly describes the coordinated vegetal movement of the deep cells, enveloping layer (EVL) and yolk syncytial layer (YSL) to engulf the yolk cell. Recently, we showed that the large GTPase Dynamin plays a fundamental role in epiboly in the early zebrafish embryo. Because Dynamin plays a well-described role in vesicle scission during endocytosis, we predicted that Dynamin might regulate epiboly through participating in bulk removal of the yolk cell membrane ahead of the advancing margin, a proposed part of the epiboly motor. Unexpectedly, we found that Dynamin function was dispensable in the yolk cell and instead, it was required to maintain the epithelial integrity of the EVL during epiboly. Here, we present a model describing the maintenance of EVL integrity, which is required for the proper generation and transmission of tension during epiboly. Furthermore, we discuss the role of Dynamin-mediated regulation of ezrin-radixin-moesin (ERM) family proteins in the maintenance of epithelial integrity.

Zebrafish Dynamin is required for maintenance of enveloping layer integrity and the progression of epiboly

Epiboly, the first morphogenetic cell movement that occurs in the zebrafish embryo, is the process by which the blastoderm thins and spreads to engulf the yolk cell. This process requires the concerted actions of the deep cells, the enveloping layer (EVL) and the extra-embryonic yolk syncytial layer (YSL). The EVL is mechanically coupled to the YSL which acts as an epiboly motor, generating the force necessary to draw the blastoderm towards the vegetal pole though actomyosin flow and contraction of the actomyosin ring. However, it has been proposed that the endocytic removal of yolk cell membrane just ahead of the advancing blastoderm may also play a role. To assess the contribution of yolk cell endocytosis in driving epiboly movements, we used a combination of drug- and dominant-negative-based approaches to inhibit Dynamin, a large GTPase with a well-characterized role in vesicle scission. We show that Dynamin-dependent endocytosis in the yolk cell is dispensable for epiboly of the blastoderm. However, global inhibition of Dynamin function revealed that Dynamin plays a fundamental role within the blastoderm during epiboly, where it maintains epithelial integrity and the transmission of tension across the EVL. The epithelial defects were associated with disrupted tight junctions and a striking reduction of cortically localized phosphorylated ezrin/radixin/moesin (P-ERM), key regulators of epithelial integrity in other systems. Furthermore, we show that Dynamin maintains EVL and promotes epiboly progression by antagonizing Rho A activity.

Fluorescent activated cell sorting (FACS) combined with gene expression microarrays for transcription enrichment profiling of zebrafish lateral line cells

Transgenic lines carrying fluorescent reporter genes like GFP have been of great value in the elucidation of developmental features and physiological processes in various animal models, including zebrafish. The lateral line (LL), which is a fish specific superficial sensory organ, is an emerging organ model for studying complex cellular processes in the context of the whole living animal. Cell migration, mechanosensory cell development/differentiation and regeneration are some examples. This sensory system is made of superficial and sparse small sensory patches called neuromasts, with less than 50 cells in any given patch. The paucity of cells is a real problem in any effort to characterize those cells at the transcriptional level. We describe here a method which we applied to efficiently separate subpopulation of cells of the LL, using two distinct stable transgenic zebrafish lines, Tg(cldnb:gfp) and Tg(tnks1bp1:EGFP). In both cases, the GFP positive (GFP+) cells were separated from the remainder of the animal by using a Fluorescent Activated Cell Sorter (FACS). The transcripts of the GFP+ cells were subsequently analyzed on gene expression microarrays. The combination of FACS and microarrays is an efficient method to establish a transcriptional signature for discrete cell populations which would otherwise be masked in whole animal preparation.

Incomplete splicing, cell division defects, and hematopoietic blockage in dhx8 mutant zebrafish

BACKGROUND

Vertebrate hematopoiesis is a complex developmental process that is controlled by genes in diverse pathways. To identify novel genes involved in early hematopoiesis, we conducted an ENU (N-ethyl-N-nitrosourea) mutagenesis screen in zebrafish. The mummy (mmy) line was investigated because of its multiple hematopoietic defects.

RESULTS

Homozygous mmy embryos lacked circulating blood cell types and were dead by 30 hr post-fertilization (hpf). The mmy mutants did not express myeloid markers and had significantly decreased expression of progenitor and erythroid markers in primitive hematopoiesis. Through positional cloning, we identified a truncation mutation in dhx8 in the mmy fish. dhx8 is the zebrafish ortholog of the yeast splicing factor prp22, which is a DEAH-box RNA helicase. mmy mutants had splicing defects in many genes, including several hematopoietic genes. mmy embryos also showed cell division defects as characterized by disorganized mitotic spindles and formation of multiple spindle poles in mitotic cells. These cell division defects were confirmed by DHX8 knockdown in HeLa cells.

CONCLUSIONS

Together, our results confirm that dhx8 is involved in mRNA splicing and suggest that it is also important for cell division during mitosis. This is the first vertebrate model for dhx8, whose function is essential for primitive hematopoiesis in developing embryos.

Nodal-dependent mesendoderm specification requires the combinatorial activities of FoxH1 and Eomesodermin

Vertebrate mesendoderm specification requires the Nodal signaling pathway and its transcriptional effector FoxH1. However, loss of FoxH1 in several species does not reliably cause the full range of loss-of-Nodal phenotypes, indicating that Nodal signals through additional transcription factors during early development. We investigated the FoxH1-dependent and -independent roles of Nodal signaling during mesendoderm patterning using a novel recessive zebrafish FoxH1 mutation called midway, which produces a C-terminally truncated FoxH1 protein lacking the Smad-interaction domain but retaining DNA–binding capability. Using a combination of gel shift assays, Nodal overexpression experiments, and genetic epistasis analyses, we demonstrate that midway more accurately represents a complete loss of FoxH1-dependent Nodal signaling than the existing zebrafish FoxH1 mutant schmalspur. Maternal-zygotic midway mutants lack notochords, in agreement with FoxH1 loss in other organisms, but retain near wild-type expression of markers of endoderm and various nonaxial mesoderm fates, including paraxial and intermediate mesoderm and blood precursors. We found that the activity of the T-box transcription factor Eomesodermin accounts for specification of these tissues in midway embryos. Inhibition of Eomesodermin in midway mutants severely reduces the specification of these tissues and effectively phenocopies the defects seen upon complete loss of Nodal signaling. Our results indicate that the specific combinations of transcription factors available for signal transduction play critical and separable roles in determining Nodal pathway output during mesendoderm patterning. Our findings also offer novel insights into the co-evolution of the Nodal signaling pathway, the notochord specification program, and the chordate branch of the deuterostome family of animals.

Multiple signaling pathways function combinatorially to form and pattern the primary tissue layers of almost all organisms, by interacting with each other and by utilizing different pathway components to perform specific roles. Here we investigated the combinatorial aspects of the Nodal signaling pathway, which is essential for proper induction of mesoderm and endoderm in vertebrates. We identified a new mutation in the zebrafish FoxH1 gene, which encodes a Nodal pathway transcription factor, a protein that responds to Nodal signals to carry out the pathway's cellular functions by regulating target gene expression. Using this mutation, we determined that FoxH1 acts in a combinatorial fashion with two other transcription factors, called Mixer and Eomesodermin, to carry out all roles of the Nodal pathway during early development. Through genetic manipulation, we were able to identify the discrete functions regulated by different combinations of these three transcription factors. Our results indicate that the availability of specific Nodal-responsive transcription factors dictates the functions of the Nodal pathway in specific areas of the developing embryo. Our work also provides evidence that the FoxH1 family of transcription factors evolved concomitantly, and perhaps causally, with the chordate branch of animals, to which all vertebrates including humans belong.

The zebrafish transcriptome during early development

Background

The transition from fertilized egg to embryo is accompanied by a multitude of changes in gene expression, and the transcriptional events that underlie these processes have not yet been fully characterized. In this study RNA-Seq is used to compare the transcription profiles of four early developmental stages in zebrafish (Danio rerio) on a global scale.

Results

An average of 79 M total reads were detected from the different stages. Out of the total number of reads 65% - 73% reads were successfully mapped and 36% - 44% out of those were uniquely mapped. The total number of detected unique gene transcripts was 11187, of which 10096 were present at 1-cell stage. The largest number of common transcripts was observed between 1-cell stage and 16-cell stage. An enrichment of gene transcripts with molecular functions of DNA binding, protein folding and processing as well as metal ion binding was observed with progression of development. The sequence data (accession number ERP000635) is available at the European Nucleotide Archive.

Conclusion

Clustering of expression profiles shows that a majority of the detected gene transcripts are present at steady levels, and thus a minority of the gene transcripts clusters as increasing or decreasing in expression over the four investigated developmental stages. The three earliest developmental stages were similar when comparing highly expressed genes, whereas the 50% epiboly stage differed from the other three stages in the identity of highly expressed genes, number of uniquely expressed genes and enrichment of GO molecular functions. Taken together, these observations indicate a major transition in gene regulation and transcriptional activity taking place between the 512-cell and 50% epiboly stages, in accordance with previous studies.

Requirement of Npc1 and availability of cholesterol for early embryonic cell movements in zebrafish

Niemann-Pick disease, type C (NP-C), often associated with Niemann-Pick disease, type C1 (NPC1) mutations, is a cholesterol-storage disorder characterized by cellular lipid accumulation, neurodegeneration, and reduced steroid production. To study NPC1 function in vivo, we cloned zebrafish npc1 and analyzed its gene expression and activity by reducing Npc1 protein with morpholino (MO)-oligonucleotides. Filipin staining in npc1-morphant cells was punctate, suggesting abnormal accumulation of cholesterol. Developmentally, reducing Npc1 did not disrupt early cell fate or survival; however, early morphogenetic movements were delayed, and the actin cytoskeleton network was abnormal. MO-induced defects were rescued with ectopic expression of mouse NPC1, demonstrating functional gene conservation, and by treatments with steroids pregnenolone or dexamethasone, suggesting that reduced steroidogenesis contributed to abnormal cell movements. Cell death was found in anterior tissues of npc1 morphants at later stages, consistent with findings in mammals. Collectively, these studies show that npc1 is required early for proper cell movement and cholesterol localization and later for cell survival.

Molecular dissection of the migrating posterior lateral line primordium during early development in zebrafish

Background

Development of the posterior lateral line (PLL) system in zebrafish involves cell migration, proliferation and differentiation of mechanosensory cells. The PLL forms when cranial placodal cells delaminate and become a coherent, migratory primordium that traverses the length of the fish to form this sensory system. As it migrates, the primordium deposits groups of cells called neuromasts, the specialized organs that contain the mechanosensory hair cells. Therefore the primordium provides both a model for studying collective directional cell migration and the differentiation of sensory cells from multipotent progenitor cells.

Results

Through the combined use of transgenic fish, Fluorescence Activated Cell Sorting and microarray analysis we identified a repertoire of key genes expressed in the migrating primordium and in differentiated neuromasts. We validated the specific expression in the primordium of a subset of the identified sequences by quantitative RT-PCR, and by in situ hybridization. We also show that interfering with the function of two genes, f11r and cd9b, defects in primordium migration are induced. Finally, pathway construction revealed functional relationships among the genes enriched in the migrating cell population.

Conclusions

Our results demonstrate that this is a robust approach to globally analyze tissue-specific expression and we predict that many of the genes identified in this study will show critical functions in developmental events involving collective cell migration and possibly in pathological situations such as tumor metastasis.

poky/chuk/ikk1 is required for differentiation of the zebrafish embryonic epidermis

An epidermis surrounds all vertebrates, forming a water barrier between the external environment and the internal space of the organism. In the zebrafish, the embryonic epidermis consists of an outer enveloping layer (EVL) and an inner basal layer that have distinct embryonic origins. Differentiation of the EVL requires the maternal effect gene poky/ikk1 in EVL cells prior to establishment of the basal layer. This requirement is transient and maternal Ikk1 is sufficient to allow establishment of the EVL and formation of normal skin in adults. Similar to the requirement for Ikk1 in mouse epidermis, EVL cells in poky mutants fail to exit the cell cycle or express specific markers of differentiation. In spite of the similarity in phenotype, the molecular requirement for Ikk1 is different between mouse and zebrafish. Unlike the mouse, EVL differentiation requires functioning Poky/Ikk1 kinase activity but does not require the HLH domain. Previous work suggested that the EVL was a transient embryonic structure, and that maturation of the epidermis required replacement of the EVL with cells from the basal layer. We show here that the EVL is not lost during embryogenesis but persists to larval stages. Our results show that while the requirement for poky/ikk1 is conserved, the differences in molecular activity indicate that diversification of an epithelial differentiation program has allowed at least two developmental modes of establishing a multilayered epidermis in vertebrates.

The tight junction component Claudin E is required for zebrafish epiboly

Zebrafish epiboly results in the thinning and spreading of the blastoderm to cover the yolk cell and close the blastopore. The extra-embryonic yolk syncytial layer (YSL) tows the blastoderm vegetally during epiboly by means of its tight junction attachments to the enveloping layer (EVL). Claudins are the major transmembrane protein components of tight junctions. Here, we focus on the function of Claudin E (Cldne), which is expressed specifically in the EVL. Morpholino knock-down of cldne produced a highly penetrant epiboly delay. Our analysis suggested that the EVL margin, which is attached to the YSL, was under reduced tension in morphant embryos. We propose that local variation in the strength of EVL-YSL attachment in morphant embryos resulted in slow and uneven advancement of the EVL and blastoderm. Our work is the first to demonstrate that Claudins are important for zebrafish epiboly.

Cumulative ligand activity of NODAL mutations and modifiers are linked to human heart defects and holoprosencephaly

The cyclopic and laterality phenotypes in model organisms linked to disturbances in the generation or propagation of Nodal-like signals are potential examples of similar impairments resulting in birth defects in humans. However, the types of gene mutation(s) and their pathogenetic combinations in humans are poorly understood. Here we describe a mutational analysis of the human NODAL gene in a large panel of patients with phenotypes compatible with diminished NODAL ligand function. Significant reductions in the biological activity of NODAL alleles are detected among patients with congenital heart defects (CHD), laterality anomalies (e.g. left-right mis-specification phenotypes), and only rarely holoprosencephaly (HPE). While many of these NODAL variants are typical for family-specific mutations, we also report the presence of alleles with significantly reduced activity among common population variants. We propose that some of these common variants act as modifiers and contribute to the ultimate phenotypic outcome in these patients; furthermore, we draw parallels with strain-specific modifiers in model organisms to bolster this interpretation.

Multiple embryo time-lapse imaging of zebrafish development

Understanding the dynamics of developmental and cellular processes requires documentation of their changes with appropriate temporal and spatial resolution. Furthermore, simultaneous recording from a population of embryos under identical conditions allows statistical estimates of precision and variability to be made. This chapter describes a protocol for time-lapse microscopy of multiple embryos in parallel developing under tightly controlled conditions. This method is currently best suited to follow tissue-scale morphogenetic movements with temporal resolution in the minute range, for hours or even days. Applications of the method include the comparison of the dynamics of a process of interest between groups of wild-type embryos and their mutant siblings or between embryos treated with different chemical compounds. Temperature control allows for the investigation of the temperature dependence of a process of interest.

Maternal Interferon Regulatory Factor 6 is required for the differentiation of primary superficial epithelia in Danio and Xenopus embryos

Early in the development of animal embryos, superficial cells of the blastula form a distinct lineage and adopt an epithelial morphology. In different animals, the fate of these primary superficial epithelial (PSE) cells varies, and it is unclear whether pathways governing segregation of blastomeres into the PSE lineage are conserved. Mutations in the gene encoding Interferon Regulatory Factor 6 (IRF6) are associated with syndromic and non-syndromic forms of cleft lip and palate, consistent with a role for Irf6 in development of oral epithelia, and mouse Irf6 targeted null mutant embryos display abnormal differentiation of oral epithelia and skin. In Danio rerio (zebrafish) and Xenopus laevis (African clawed frog) embryos, zygotic irf6 transcripts are present in many epithelial tissues including the presumptive PSE cells and maternal irf6 transcripts are present throughout all cells at the blastula stage. Injection of antisense oligonucleotides with ability to disrupt translation of irf6 transcripts caused little or no effect on development. By contrast, injection of RNA encoding a putative dominant negative Irf6 caused epiboly arrest, loss of gene expression characteristic of the EVL, and rupture of the embryo at late gastrula stage. The dominant negative Irf6 disrupted EVL gene expression in a cell autonomous fashion. These results suggest that Irf6 translated in the oocyte or unfertilized egg suffices for early development. Supporting the importance of maternal Irf6, we show that depletion of maternal irf6 transcripts in X. laevis embryos leads to gastrulation defects and rupture of the superficial epithelium. These experiments reveal a conserved role for maternally-encoded Irf6 in differentiation of a simple epithelium in X. laevis and D. rerio. This epithelium constitutes a novel model tissue in which to explore the Irf6 regulatory pathway.

Transcriptional profiling of endogenous germ layer precursor cells identifies dusp4 as an essential gene in zebrafish endoderm specification

A major goal for developmental biologists is to define the behaviors and molecular contents of differentiating cells. We have devised a strategy for isolating cells from diverse embryonic regions and stages in the zebrafish, using computer-guided laser photoconversion of injected Kaede protein and flow cytometry. This strategy enabled us to perform a genome-wide transcriptome comparison of germ layer precursor cells. Mesendoderm and ectoderm precursors cells isolated by this method differentiated appropriately in transplantation assays. Microarray analysis of these cells reidentified known genes at least as efficiently as previously reported strategies that relied on artificial mesendoderm activation or inhibition. We also identified a large set of uncharacterized mesendoderm-enriched genes as well as ectoderm-enriched genes. Loss-of-function studies revealed that one of these genes, the MAP kinase inhibitor dusp4, is essential for early development. Embryos injected with antisense morpholino oligonucleotides that targeted Dusp4 displayed necrosis of head tissues. Marker analysis during late gastrulation revealed a specific loss of sox17, but not of other endoderm markers, and analysis at later stages revealed a loss of foregut and pancreatic endoderm. This specific loss of sox17 establishes a new class of endoderm specification defect.

Reduced NODAL signaling strength via mutation of several pathway members including FOXH1 is linked to human heart defects and holoprosencephaly

Abnormalities of embryonic patterning are hypothesized to underlie many common congenital malformations in humans including congenital heart defects (CHDs), left-right disturbances (L-R) or laterality, and holoprosencephaly (HPE). Studies in model organisms suggest that Nodal-like factors provide instructions for key aspects of body axis and germ layer patterning; however, the complex genetics of pathogenic gene variant(s) in humans are poorly understood. Here we report our studies of FOXH1, CFC1, and SMAD2 and summarize our mutational analysis of three additional components in the human NODAL-signaling pathway: NODAL, GDF1, and TDGF1. We identify functionally abnormal gene products throughout the pathway that are clearly associated with CHD, laterality, and HPE. Abnormal gene products are most commonly detected in patients within a narrow spectrum of isolated conotruncal heart defects (minimum 5%-10% of subjects), and far less commonly in isolated laterality or HPE patients (approximately 1% for each). The difference in the mutation incidence between these groups is highly significant. We show that apparent gene dosage discrepancies between humans and model organisms can be reconciled by considering a broader combination of sequence variants. Our studies confirm that (1) the genetic vulnerabilities inferred from model organisms with defects in Nodal signaling are indeed analogous to humans; (2) the molecular analysis of an entire signaling pathway is more complete and robust than that of individual genes and presages future studies by whole-genome analysis; and (3) a functional genomics approach is essential to fully appreciate the complex genetic interactions necessary to produce these effects in humans.