FGF8/17/18 functions together with FGF9/16/20 during formation of the notochord in Ciona embryos

Early transcriptional similarities between two distinct neural lineages during ascidian embryogenesis

In chordates, the central nervous system arises from precursors that have distinct developmental and transcriptional trajectories. Anterior nervous systems are ontogenically associated with ectodermal lineages while posterior nervous systems are associated with mesoderm. Taking advantage of the well-documented cell lineage of ascidian embryos, we asked to what extent the transcriptional states of the different neural lineages become similar during the course of progressive lineage restriction. We performed single-cell RNA sequencing (scRNA-seq) analyses on hand-dissected neural precursor cells of the two distinct lineages, together with those of their sister cell lineages, with a high temporal resolution covering five successive cell cycles from the 16-cell to neural plate stages. A transcription factor binding site enrichment analysis of neural specific genes at the neural plate stage revealed limited evidence for shared transcriptional control between the two neural lineages, consistent with their different ontogenies. Nevertheless, PCA analysis and hierarchical clustering showed that, by neural plate stages, the two neural lineages cluster together. Consistent with this, we identified a set of genes enriched in both neural lineages at the neural plate stage, including miR-124, Celf3.a, Zic.r-b, and Ets1/2. Altogether, the current study has revealed genome-wide transcriptional dynamics of neural progenitor cells of two distinct developmental origins. Our scRNA-seq dataset is unique and provides a valuable resource for future analyses, enabling a precise temporal resolution of cell types not previously described from dissociated embryos.

Temporospatial hierarchy and allele-specific expression of zygotic genome activation revealed by distant interspecific urochordate hybrids

Zygotic genome activation (ZGA) is a universal process in early embryogenesis of metazoan, when the quiescent zygotic nucleus initiates global transcription. However, the mechanisms related to massive genome activation and allele-specific expression (ASE) remain not well understood. Here, we develop hybrids from two deeply diverged (120 Mya) ascidian species to symmetrically document the dynamics of ZGA. We identify two coordinated ZGA waves represent early developmental and housekeeping gene reactivation, respectively. Single-cell RNA sequencing reveals that the major expression wave exhibits spatial heterogeneity and significantly correlates with cell fate. Moreover, allele-specific expression occurs in a species- rather than parent-related manner, demonstrating the divergence of cis-regulatory elements between the two species. These findings provide insights into ZGA in chordates.

Diverse logics and grammar encode notochord enhancers

The notochord is a defining feature of all chordates. The transcription factors Zic and ETS regulate enhancer activity within the notochord. We conduct high-throughput screens of genomic elements within developing Ciona embryos to understand how Zic and ETS sites encode notochord activity. Our screen discovers an enhancer located near Lama, a gene critical for notochord development. Reversing the orientation of an ETS site within this enhancer abolishes expression, indicating that enhancer grammar is critical for notochord activity. Similarly organized clusters of Zic and ETS sites occur within mouse and human Lama1 introns. Within a Brachyury (Bra) enhancer, FoxA and Bra, in combination with Zic and ETS binding sites, are necessary and sufficient for notochord expression. This binding site logic also occurs within other Ciona and vertebrate Bra enhancers. Collectively, this study uncovers the importance of grammar within notochord enhancers and discovers signatures of enhancer logic and grammar conserved across chordates.

Quantitative Dissection of the Proximal Ciona brachyury Enhancer

A major goal in biology is to understand the rules by which cis-regulatory sequences control spatially and temporally precise expression patterns. Here we present a systematic dissection of the proximal enhancer for the notochord-specific transcription factor brachyury in the ascidian chordate Ciona. The study uses a quantitative image-based reporter assay that incorporates a dual-reporter strategy to control for variable electroporation efficiency. We identified and mutated multiple predicted transcription factor binding sites of interest based on statistical matches to the JASPAR binding motif database. Most sites (Zic, Ets, FoxA, RBPJ) were selected based on prior knowledge of cell fate specification in both the primary and secondary notochord. We also mutated predicted Brachyury sites to investigate potential autoregulation as well as Fos/Jun (AP1) sites that had very strong matches to JASPAR. Our goal was to quantitatively define the relative importance of these different sites, to explore the importance of predicted high-affinity versus low-affinity motifs, and to attempt to design mutant enhancers that were specifically expressed in only the primary or secondary notochord lineages. We found that the mutation of all predicted high-affinity sites for Zic, FoxA or Ets led to quantifiably distinct effects. The FoxA construct caused a severe loss of reporter expression whereas the Ets construct had little effect. A strong Ets phenotype was only seen when much lower-scoring binding sites were also mutated. This supports the enhancer suboptimization hypothesis proposed by Farley and Levine but suggests that it may only apply to some but not all transcription factor families. We quantified reporter expression separately in the two notochord lineages with the expectation that Ets mutations and RBPJ mutations would have distinct effects given that primary notochord is induced by Ets-mediated FGF signaling whereas secondary notochord is induced by RBPJ/Su(H)-mediated Notch/Delta signaling. We found, however, that ETS mutations affected primary and secondary notochord expression relatively equally and that RBPJ mutations were only moderately more severe in their effect on secondary versus primary notochord. Our results point to the promise of quantitative reporter assays for understanding cis-regulatory logic but also highlight the challenge of arbitrary statistical thresholds for predicting potentially important sites.

Single-cell analysis of cell fate bifurcation in the chordate Ciona

BACKGROUND

Inductive signaling interactions between different cell types are a major mechanism for the further diversification of embryonic cell fates. Most blastomeres in the model chordate Ciona robusta become restricted to a single predominant fate between the 64-cell and mid-gastrula stages. The deeply stereotyped and well-characterized Ciona embryonic cell lineages allow the transcriptomic analysis of newly established cell types very early in their divergence from sibling cell states without the pseudotime inference needed in the analysis of less synchronized cell populations. This is the first ascidian study to use droplet scRNAseq with large numbers of analyzed cells as early as the 64-cell stage when major lineages such as primary notochord first become fate restricted.

RESULTS AND CONCLUSIONS

We identify 59 distinct cell states, including new subregions of the b-line neural lineage and the early induction of the tail tip epidermis. We find that 34 of these cell states are directly or indirectly dependent on MAPK-mediated signaling critical to early Ciona patterning. Most of the MAPK-dependent bifurcations are canalized with the signal-induced cell fate lost upon MAPK inhibition, but the posterior endoderm is unique in being transformed into a novel state expressing some but not all markers of both endoderm and muscle. Divergent gene expression between newly bifurcated sibling cell types is dominated by upregulation in the induced cell type. The Ets family transcription factor Elk1/3/4 is uniquely upregulated in nearly all the putatively direct inductions. Elk1/3/4 upregulation together with Ets transcription factor binding site enrichment analysis enables inferences about which bifurcations are directly versus indirectly controlled by MAPK signaling. We examine notochord induction in detail and find that the transition between a Zic/Ets-mediated regulatory state and a Brachyury/FoxA-mediated regulatory state is unexpectedly late. This supports a "broad-hourglass" model of cell fate specification in which many early tissue-specific genes are induced in parallel to key tissue-specific transcriptional regulators via the same set of transcriptional inputs.

FGF signaling acts on different levels of mesoderm development within Spiralia

FGF signaling is involved in mesoderm induction in members of deuterostomes (e.g. tunicates, hemichordates), but not in flies and nematodes, in which it has a role in mesoderm patterning and migration. However, we need comparable studies in other protostome taxa in order to decipher whether this mesoderm-inducing function of FGF extends beyond the lineage of deuterostomes. Here, we investigated the role of FGF signaling in mesoderm development in three species of lophophorates, a clade within the protostome group Spiralia. Our gene expression analyses show that the mesodermal molecular patterning is conserved between brachiopods and phoronids, but the spatial and temporal recruitment of transcription factors differs significantly. Moreover, the use of the inhibitor SU5402 demonstrates that FGF signaling is involved in different steps of mesoderm development, as well as in morphogenetic movements of gastrulation and axial elongation. Our findings suggest that the mesoderm-inducing role of FGF extends beyond the group of deuterostomes.

Using linkage logic theory to control dynamics of a gene regulatory network of a chordate embryo

Linkage logic theory provides a mathematical criterion to control network dynamics by manipulating activities of a subset of network nodes, which are collectively called a feedback vertex set (FVS). Because many biological functions emerge from dynamics of biological networks, this theory provides a promising tool for controlling biological functions. By manipulating the activity of FVS molecules identified in a gene regulatory network (GRN) for fate specification of seven tissues in ascidian embryos, we previously succeeded in reproducing six of the seven cell types. Simultaneously, we discovered that the experimentally reconstituted GRN lacked information sufficient to reproduce muscle cells. Here, we utilized linkage logic theory as a tool to find missing edges in the GRN. Then, we identified a FVS from an updated version of the GRN and confirmed that manipulating the activity of this FVS was sufficient to induce all seven cell types, even in a multi-cellular environment. Thus, linkage logic theory provides tools to find missing edges in experimentally reconstituted networks, to determine whether reconstituted networks contain sufficient information to fulfil expected functions, and to reprogram cell fate.

Brachyury controls Ciona notochord fate as part of a feed-forward network

The notochord is a defining feature of the chordates. The transcription factor Brachyury (Bra) is a key regulator of notochord fate but here we show that it is not a unitary master regulator in the model chordate Ciona Ectopic Bra expression only partially reprograms other cell types to a notochord-like transcriptional profile and a subset of notochord-enriched genes is unaffected by CRISPR Bra disruption. We identify Foxa.a and Mnx as potential co-regulators, and find that combinatorial cocktails are more effective at reprogramming other cell types than Bra alone. We reassess the network relationships between Bra, Foxa.a and other components of the notochord gene regulatory network, and find that Foxa.a expression in the notochord is regulated by vegetal FGF signaling. It is a direct activator of Bra expression and has a binding motif that is significantly enriched in the regulatory regions of notochord-enriched genes. These and other results indicate that Bra and Foxa.a act together in a regulatory network dominated by positive feed-forward interactions, with neither being a classically defined master regulator.

Molecular insights into deuterostome evolution from hemichordate developmental biology

Hemichordates, along with echinoderms and chordates, belong to the lineage of bilaterians called the deuterostomes. Their phylogenetic position as an outgroup to chordates provides an opportunity to investigate the evolutionary origins of the chordate body plan and reconstruct ancestral deuterostome characters. The body plans of the hemichordates and chordates are organizationally divergent making anatomical comparisons very challenging. The developmental underpinnings of animal body plans are often more conservative than the body plans they regulate, and offer a novel data set for making comparisons between morphologically divergent body architectures. Here I review the hemichordate developmental data generated over the past 20 years that further test hypotheses of proposed morphological affinities between the two taxa, but also compare the conserved anteroposterior, dorsoventral axial patterning programs and germ layer specification programs. These data provide an opportunity to determine which developmental programs are ancestral deuterostome or bilaterian innovations, and which ones occurred in stem chordates or vertebrates representing developmental novelties of the chordate body plan.

Ciona Brachyury proximal and distal enhancers have different FGF dose-response relationships

Many genes are regulated by two or more enhancers that drive similar expression patterns. Evolutionary theory suggests that these seemingly redundant enhancers must have functionally important differences. In the simple ascidian chordate Ciona, the transcription factor Brachyury is induced exclusively in the presumptive notochord downstream of lineage specific regulators and FGF-responsive Ets family transcription factors. Here we exploit the ability to finely titrate FGF signaling activity via the MAPK pathway using the MEK inhibitor U0126 to quantify the dependence of transcription driven by different Brachyury reporter constructs on this direct upstream regulator. We find that the more powerful promoter-adjacent proximal enhancer and a weaker distal enhancer have fundamentally different dose-response relationships to MAPK inhibition. The Distal enhancer is more sensitive to MAPK inhibition but shows a less cooperative response, whereas the Proximal enhancer is less sensitive and more cooperative. A longer construct containing both enhancers has a complex dose-response curve that supports the idea that the proximal and distal enhancers are moderately super-additive. We show that the overall expression loss from intermediate doses of U0126 is not only a function of the fraction of cells expressing these reporters, but also involves graded decreases in expression at the single-cell level. Expression of the endogenous gene shows a comparable dose-response relationship to the full length reporter, and we find that different notochord founder cells are differentially sensitive to MAPK inhibition. Together, these results indicate that although the two Brachyury enhancers have qualitatively similar expression patterns, they respond to FGF in quantitatively different ways and act together to drive high levels of Brachyury expression with a characteristic input/output relationship. This indicates that they are fundamentally not equivalent genetic elements.

When and where genes are expressed is controlled by regulatory DNA regions known as enhancers. Genes often have multiple enhancers that control expression in different cell types or embryonic regions, but there are also genes that have multiple enhancers that control near-identical expression patterns. These ‘shadow’ enhancers are common features of many animal genomes, but it is unclear to what extent they are truly identical in function. Here we studied a pair of shadow enhancers for the notochord-specific gene Brachyury in the simple model chordate Ciona that are both directly activated by the same signaling pathway. We titrated the activity of this pathway with graded doses of a pharmacological inhibitor and measured the effects in quantitative enhancer assays. We found that the two enhancers had significant differences in sensitivity and cooperativity to the same shared regulator and are thus not identical in function when assessed quantitatively. We also identified subtle differences in sensitivity to this upstream signal between different notochord precursor cells.

Contact area-dependent cell communication and the morphological invariance of ascidian embryogenesis

Marine invertebrate ascidians display embryonic reproducibility: Their early embryonic cell lineages are considered invariant and are conserved between distantly related species, despite rapid genomic divergence. Here, we address the drivers of this reproducibility. We used light-sheet imaging and automated cell segmentation and tracking procedures to systematically quantify the behavior of individual cells every 2 minutes during Phallusia mammillata embryogenesis. Interindividual reproducibility was observed down to the area of individual cell contacts. We found tight links between the reproducibility of embryonic geometries and asymmetric cell divisions, controlled by differential sister cell inductions. We combined modeling and experimental manipulations to show that the area of contact between signaling and responding cells is a key determinant of cell communication. Our work establishes the geometric control of embryonic inductions as an alternative to classical morphogen gradients and suggests that the range of cell signaling sets the scale at which embryonic reproducibility is observed.

The notochord gene regulatory network in chordate evolution: Conservation and divergence from Ciona to vertebrates

The notochord is a structure required for support and patterning of all chordate embryos, from sea squirts to humans. An increasing amount of information on notochord development and on the molecular strategies that ensure its proper morphogenesis has been gleaned through studies in the sea squirt Ciona. This invertebrate chordate offers a fortunate combination of experimental advantages, ranging from translucent, fast-developing embryos to a compact genome and impressive biomolecular resources. These assets have enabled the rapid identification of numerous notochord genes and cis-regulatory regions, and provide a rather unique opportunity to reconstruct the gene regulatory network that controls the formation of this developmental and evolutionary chordate landmark. This chapter summarizes the morphogenetic milestones that punctuate notochord formation in Ciona, their molecular effectors, and the current knowledge of the gene regulatory network that ensures the accurate spatial and temporal orchestration of these processes.

A Nodal/Eph signalling relay drives the transition from apical constriction to apico-basal shortening in ascidian endoderm invagination

Gastrulation is the first major morphogenetic event during animal embryogenesis. Ascidian gastrulation starts with the invagination of 10 endodermal precursor cells between the 64- and late 112-cell stages. This process occurs in the absence of endodermal cell division and in two steps, driven by myosin-dependent contractions of the acto-myosin network. First, endoderm precursors constrict their apex. Second, they shorten apico-basally, while retaining small apical surfaces, thereby causing invagination. The mechanisms that prevent endoderm cell division, trigger the transition between step 1 and step 2, and drive apico-basal shortening have remained elusive. Here, we demonstrate a conserved role for Nodal and Eph signalling during invagination in two distantly related ascidian species, Phallusia mammillata and Ciona intestinalis. Specifically, we show that the transition to step 2 is triggered by Nodal relayed by Eph signalling. Additionally, our results indicate that Eph signalling lengthens the endodermal cell cycle, independently of Nodal. Finally, we find that both Nodal and Eph signals are dispensable for endoderm fate specification. These results illustrate commonalities as well as differences in the action of Nodal during ascidian and vertebrate gastrulation.

Tunicate gastrulation

Tunicates are a diverse group of invertebrate marine chordates that includes the larvaceans, thaliaceans, and ascidians. Because of their unique evolutionary position as the sister group of the vertebrates, tunicates are invaluable as a comparative model and hold the promise of revealing both conserved and derived features of chordate gastrulation. Descriptive studies in a broad range of tunicates have revealed several important unifying traits that make them unique among the chordates, including invariant cell lineages through gastrula stages and an overall morphological simplicity. Gastrulation has only been studied in detail in ascidians such as Ciona and Phallusia, where it involves a simple cup-shaped gastrula driven primarily by endoderm invagination. This appears to differ significantly from vertebrate models, such as Xenopus, in which mesoderm convergent extension and epidermal epiboly are major contributors to involution. These differences may reflect the cellular simplicity of the ascidian embryo.

Transcriptional regulation of the Ciona Gsx gene in the neural plate

The ascidian neural plate consists of a defined number of identifiable cells organized in a grid of rows and columns, representing a useful model to investigate the molecular mechanisms controlling neural patterning in chordates. Distinct anterior brain lineages are specified via unique combinatorial inputs of signalling pathways with Nodal and Delta-Notch signals patterning along the medial-lateral axis and FGF/MEK/ERK signals patterning along the anterior-posterior axis of the neural plate. The Ciona Gsx gene is specifically expressed in the a9.33 cells in the row III/column 2 position of anterior brain lineages, characterised by a combinatorial input of Nodal-OFF, Notch-ON and FGF-ON. Here, we identify the minimal cis-regulatory element (CRE) of 376 bp, which can recapitulate the early activation of Gsx. We show that this minimal CRE responds in the same way as the endogenous Gsx gene to manipulation of FGF- and Notch-signalling pathways and to overexpression of Snail, a mediator of Nodal signals, and Six3/6, which is required to demarcate the anterior boundary of Gsx expression at the late neurula stage. We reveal that sequences proximal to the transcription start site include a temporal regulatory element required for the precise transcriptional onset of gene expression. We conclude that sufficient spatial and temporal information for Gsx expression is integrated in 376 bp of non-coding cis-regulatory sequences.

Multiple inputs into a posterior-specific regulatory network in the Ciona notochord

The gene regulatory networks underlying Ciona notochord fate specification and differentiation have been extensively investigated, but the regulatory basis for regionalized expression within the notochord is not understood. Here we identify three notochord-expressed genes, C11.331, C12.115 and C8.891, with strongly enriched expression in the secondary notochord cells at the posterior tip of the tail. C11.331 and C12.115 share a distinctive expression pattern that is highly enriched in the secondary notochord lineage but also graded within that lineage with the strongest expression at the posterior tip. Both genes show similar responses to pharmacological perturbations of Wnt and FGF signaling, consistent with an important role for Wnt and FGF ligands expressed at the tail tip. Reporter analysis indicates that the C11.331 cis-regulatory regions are extensively distributed, with multiple non-overlapping regions conferring posterior notochord-enriched expression. Fine-scale analysis of a minimal cis-regulatory module identifies discrete positive and negative elements including a strong silencer. Truncation of the silencer region leads to increased expression in the primary notochord, indicating that C11.331 expression is influenced by putative regulators of primary versus secondary notochord fate. The minimal CRM contains predicted ETS, GATA, LMX and Myb sites, all of which lead to reduced expression in secondary notochord when mutated. These results show that the posterior-enriched notochord expression of C11.331 depends on multiple inputs, including Wnt and FGF signals from the tip of the tail, multiple notochord-specific regulators, and yet-to-be identified regulators of regional identity within the notochord.

Reiterative use of FGF signaling in mesoderm development during embryogenesis and metamorphosis in the hemichordate Ptychodera flava

BACKGROUND

Mesoderm is generally considered to be a germ layer that is unique to Bilateria, and it develops into diverse tissues, including muscle, and in the case of vertebrates, the skeleton and notochord. Studies on various deuterostome animals have demonstrated that fibroblast growth factor (FGF) signaling is required for the formation of many mesodermal structures, such as vertebrate somites, from which muscles are differentiated, and muscles in sea urchin embryos, suggesting an ancient role of FGF signaling in muscle development. However, the formation of trunk muscles in invertebrate chordates is FGF-independent, leading to ambiguity about this ancient role in deuterostomes. To further understand the role of FGF signaling during deuterostome evolution, we investigated the development of mesodermal structures during embryogenesis and metamorphosis in Ptychodera flava, an indirect-developing hemichordate that has larval morphology similar to echinoderms and adult body features that are similar to chordates.

RESULTS

Here we show that genes encoding FGF ligands, FGF receptors and transcription factors that are known to be involved in mesoderm formation and myogenesis are expressed dynamically during embryogenesis and metamorphosis. FGF signaling at the early gastrula stage is required for the specification of the mesodermal cell fate in P. flava. The mesoderm cells are then differentiated stepwise into the hydroporic canal, the pharyngeal muscle and the muscle string; formation of the last two muscular structures are controlled by FGF signaling. Moreover, augmentation of FGF signaling during metamorphosis accelerated the process, facilitating the transformation from cilia-driven swimming larvae into muscle-driven worm-like juveniles.

CONCLUSIONS

Our data show that FGF signaling is required for mesoderm induction and myogenesis in the P. flava embryo, and it is reiteratively used for the morphological transition during metamorphosis. The dependence of muscle development on FGF signaling in both planktonic larvae and sand-burrowing worms supports its ancestral role in deuterostomes.

Controlling Cell Fate Specification System by Key Genes Determined from Network Structure

Network structures describing regulation between biomolecules have been determined in many biological systems. Dynamics of molecular activities based on such networks are considered to be the origin of many biological functions. Recently, it has been proved mathematically that key nodes for controlling dynamics in networks are identified from network structure alone. Here, we applied this theory to a gene regulatory network for the cell fate specification of seven tissues in the ascidian embryo and found that this network, which consisted of 92 factors, had five key molecules. By controlling the activities of these key molecules, the specific gene expression of six of seven tissues observed in the embryo was successfully reproduced. Since this method is applicable to all nonlinear dynamic systems, we propose this method as a tool for controlling gene regulatory networks and reprogramming cell fates.

An FGF-driven feed-forward circuit patterns the cardiopharyngeal mesoderm in space and time

In embryos, multipotent progenitors divide to produce distinct progeny and express their full potential. In vertebrates, multipotent cardiopharyngeal progenitors produce second-heart-field-derived cardiomyocytes, and branchiomeric skeletal head muscles. However, the mechanisms underlying these early fate choices remain largely elusive. The tunicate Ciona emerged as an attractive model to study early cardiopharyngeal development at high resolution: through two asymmetric and oriented divisions, defined cardiopharyngeal progenitors produce distinct first and second heart precursors, and pharyngeal muscle (aka atrial siphon muscle, ASM) precursors. Here, we demonstrate that differential FGF-MAPK signaling distinguishes between heart and ASM precursors. We characterize a feed-forward circuit that promotes the successive activations of essential ASM determinants, Hand-related, Tbx1/10 and Ebf. Finally, we show that coupling FGF-MAPK restriction and cardiopharyngeal network deployment with cell divisions defines the timing of gene expression and permits the emergence of diverse cell types from multipotent progenitors.

Regulatory Principles Governing Tissue Specificity of Developmental Enhancers

Transcriptional enhancers are short segments of genomic DNA (50 bp to 1 kb in length) that can work over long distances (≥1 Mb) to regulate gene expression in specific cells and tissues. Genomic assays have identified on the order of 400,000 to one million putative enhancers in the human genome (e.g., ENCODE Consortium). This suggests that a typical gene is regulated by tens of enhancers, ensuring stringent regulation of gene expression in response to a variety of intrinsic and external signals. Despite the discovery of the first transcriptional enhancer more than 30 years ago, we know surprisingly little about how enhancers regulate gene expression. In particular, the relationship between primary DNA sequence and enhancer specificity remains obscure. Here we summarize recent high-throughput studies in whole embryos aimed at the systematic identification of the sequence and organizational constraints underlying enhancer function and specificity.

A bipolar role of the transcription factor ERG for cnidarian germ layer formation and apical domain patterning

Germ layer formation and axial patterning are biological processes that are tightly linked during embryonic development of most metazoans. In addition to canonical WNT, it has been proposed that ERK-MAPK signaling is involved in specifying oral as well as aboral territories in cnidarians. However, the effector and the molecular mechanism underlying latter phenomenon is unknown. By screening for potential effectors of ERK-MAPK signaling in both domains, we identified a member of the ETS family of transcription factors, Nverg that is bi-polarily expressed prior to gastrulation. We further describe the crucial role of NvERG for gastrulation, endomesoderm as well as apical domain formation. The molecular characterization of the obtained NvERG knock-down phenotype using previously described as well as novel potential downstream targets, provides evidence that a single transcription factor, NvERG, simultaneously controls expression of two different sets of downstream targets, leading to two different embryonic gene regulatory networks (GRNs) in opposite poles of the developing embryo. We also highlight the molecular interaction of cWNT and MEK/ERK/ERG signaling that provides novel insight into the embryonic axial organization of Nematostella, and show a cWNT repressive role of MEK/ERK/ERG signaling in segregating the endomesoderm in two sub-domains, while a common input of both pathways is required for proper apical domain formation. Taking together, we build the first blueprint for a global cnidarian embryonic GRN that is the foundation for additional gene specific studies addressing the evolution of embryonic and larval development.

FGF and canonical Wnt signaling cooperate to induce paraxial mesoderm from tailbud neuromesodermal progenitors through regulation of a two-step epithelial to mesenchymal transition

Mesoderm induction begins during gastrulation. Recent evidence from several vertebrate species indicates that mesoderm induction continues after gastrulation in neuromesodermal progenitors (NMPs) within the posteriormost embryonic structure, the tailbud. It is unclear to what extent the molecular mechanisms of mesoderm induction are conserved between gastrula and post-gastrula stages of development. Fibroblast growth factor (FGF) signaling is required for mesoderm induction during gastrulation through positive transcriptional regulation of the T-box transcription factor brachyury We find in zebrafish that FGF is continuously required for paraxial mesoderm (PM) induction in post-gastrula NMPs. FGF signaling represses the NMP markers brachyury (ntla) and sox2 through regulation of tbx16 and msgn1, thereby committing cells to a PM fate. FGF-mediated PM induction in NMPs functions in tight coordination with canonical Wnt signaling during the epithelial to mesenchymal transition (EMT) from NMP to mesodermal progenitor. Wnt signaling initiates EMT, whereas FGF signaling terminates this event. Our results indicate that germ layer induction in the zebrafish tailbud is not a simple continuation of gastrulation events.

Differential temporal control of Foxa.a and Zic-r.b specifies brain versus notochord fate in the ascidian embryo

In embryos of an invertebrate chordate, Ciona intestinalis, two transcription factors, Foxa.a and Zic-r.b, are required for specification of the brain and the notochord, which are derived from distinct cell lineages. In the brain lineage, Foxa.a and Zic-r.b are expressed with no temporal overlap. In the notochord lineage, Foxa.a and Zic-r.b are expressed simultaneously. In the present study, we found that the temporally non-overlapping expression of Foxa.a and Zic-r.b in the brain lineage was regulated by three repressors: Prdm1-r.a (formerly called BZ1), Prdm1-r.b (BZ2) and Hes.a. In morphant embryos of these three repressor genes, Foxa.a expression was not terminated at the normal time, and Zic-r.b was precociously expressed. Consequently, Foxa.a and Zic-r.b were expressed simultaneously, which led to ectopic activation of Brachyury and its downstream pathways for notochord differentiation. Thus, temporal controls by transcriptional repressors are essential for specification of the two distinct fates of brain and notochord by Foxa.a and Zic-r.b Such a mechanism might enable the repeated use of a limited repertoire of transcription factors in developmental gene regulatory networks.

Co-expression of Foxa.a, Foxd and Fgf9/16/20 defines a transient mesendoderm regulatory state in ascidian embryos

In many bilaterian embryos, nuclear β-catenin (nβ-catenin) promotes mesendoderm over ectoderm lineages. Although this is likely to represent an evolutionary ancient developmental process, the regulatory architecture of nβ-catenin-induced mesendoderm remains elusive in the majority of animals. Here, we show that, in ascidian embryos, three nβ-catenin transcriptional targets, Foxa.a, Foxd and Fgf9/16/20, are each required for the correct initiation of both the mesoderm and endoderm gene regulatory networks. Conversely, these three factors are sufficient, in combination, to produce a mesendoderm ground state that can be further programmed into mesoderm or endoderm lineages. Importantly, we show that the combinatorial activity of these three factors is sufficient to reprogramme developing ectoderm cells to mesendoderm. We conclude that in ascidian embryos, the transient mesendoderm regulatory state is defined by co-expression of Foxa.a, Foxd and Fgf9/16/20.

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

Syntax compensates for poor binding sites to encode tissue specificity of developmental enhancers

Transcriptional enhancers are short segments of DNA that switch genes on and off in response to a variety of intrinsic and extrinsic signals. Despite the discovery of the first enhancer more than 30 y ago, the relationship between primary DNA sequence and enhancer activity remains obscure. In particular, the importance of "syntax" (the order, orientation, and spacing of binding sites) is unclear. A high-throughput screen identified synthetic notochord enhancers that are activated by the combination of ZicL and ETS transcription factors in Ciona embryos. Manipulation of these enhancers elucidated a "regulatory code" of sequence and syntax features for notochord-specific expression. This code enabled in silico discovery of bona fide notochord enhancers, including those containing low-affinity binding sites that would be excluded by standard motif identification methods. One of the newly identified enhancers maps upstream of the known enhancer that regulates Brachyury (Ci-Bra), a key determinant of notochord specification. This newly identified Ci-Bra shadow enhancer contains binding sites with very low affinity, but optimal syntax, and therefore mediates surprisingly strong expression in the notochord. Weak binding sites are compensated by optimal syntax, whereas enhancers containing high-affinity binding affinities possess suboptimal syntax. We suggest this balance has obscured the importance of regulatory syntax, as noncanonical binding motifs are typically disregarded by enhancer detection methods. As a result, enhancers with low binding affinities but optimal syntax may be a vastly underappreciated feature of the regulatory genome.

Mesoderm patterning and morphogenesis in the polychaete Alitta virens (Spiralia, Annelida): Expression of mesodermal markers Twist, Mox, Evx and functional role for MAP kinase signaling

Mesoderm represents the evolutionary youngest germ layer and forms numerous novel tissues in bilaterian animals. Despite the established conservation of the gene regulatory networks that drive mesoderm differentiation (e.g. myogenesis), mechanisms of mesoderm specification are highly variable in distant model species. Thus, broader phylogenetic sampling is required to reveal common features of mesoderm formation across bilaterians. Here we focus on a representative of Spiralia, the marine annelid Alitta virens, whose mesoderm development is still poorly investigated on the molecular level. We characterize three novel early mesodermal markers for A. virens - Twist, Mox, and Evx - which are differentially expressed within the mesodermal lineages. The Twist mRNA is ubiquitously distributed in the fertilized egg and exhibits specific expression in endomesodermal- and ectomesodermal-founder cells at gastrulation. Twist is expressed around the blastopore and later in a segmental metameric pattern. We consider this expression to be ancestral, and in support of the enterocoelic hypothesis of mesoderm evolution. We also revealed an early pattern of the MAPK activation in A. virens that is different from the previously reported pattern in spiralians. Inhibition of the MAPK pathway by U0126 disrupts the metameric Twist and Mox expression, indicating an early requirement of the MAPK cascade for proper morphogenesis of endomesodermal tissues.

Suboptimization of developmental enhancers

Transcriptional enhancers direct precise on-off patterns of gene expression during development. To explore the basis for this precision, we conducted a high-throughput analysis of the Otx-a enhancer, which mediates expression in the neural plate of Ciona embryos in response to fibroblast growth factor (FGF) signaling and a localized GATA determinant. We provide evidence that enhancer specificity depends on submaximal recognition motifs having reduced binding affinities ("suboptimization"). Native GATA and ETS (FGF) binding sites contain imperfect matches to consensus motifs. Perfect matches mediate robust but ectopic patterns of gene expression. The native sites are not arranged at optimal intervals, and subtle changes in their spacing alter enhancer activity. Multiple tiers of enhancer suboptimization produce specific, but weak, patterns of expression, and we suggest that clusters of weak enhancers, including certain "superenhancers," circumvent this trade-off in specificity and activity.

Chondroitin 6-O-sulfotransferases are required for morphogenesis of the notochord in the ascidian embryo

BACKGROUND

Chondroitin sulfate (CS) is a sulfated polysaccharide chain that binds to various core proteins to form proteoglycans. The amount and position of sulfate groups in CS are variable among different tissues, and are determined by specific sulfotransferases. Although the ascidians are the closest relatives of vertebrates, the functions of their sulfotransferases have not been studied.

RESULTS

The genome of the ascidian Ciona intestinalis contains eight genes encoding proteins similar to chondroitin 6-O-sulfotransferases (C6STs), which appear to have independently diverged in the ascidian lineage during evolution. Among them, Ci-C6ST-like1 and Ci-C6ST-like7 were predominantly expressed in the developing notochord. In addition, they were weakly expressed in the neural tube. The disruption of either one of them affected the convergent extension movement of notochordal cells. Presumptive notochord cells coming from both sides of the embryo did not intercalate. The results suggest that both of them are necessary. In some cases, the anterior neural tube failed to close. Forced expression of Ci-C6ST-like1 or Ci-C6ST-like7 in the notochord restored the normal intercalation of notochordal cells, indicating that the effects of morpholino oligos are specific.

CONCLUSIONS

Ci-C6ST-like1 and Ci-C6ST-like7 are required for the morphogenesis of the notochord in the ascidian embryo.

Insulin-like genes in ascidians: findings in Ciona and hypotheses on the evolutionary origins of the pancreas

Insulin plays an extensively characterized role in the control of sugar metabolism, growth and homeostasis in a wide range of organisms. In vertebrate chordates, insulin is mainly produced by the beta cells of the endocrine pancreas, while in non-chordate animals insulin-producing cells are mainly found in the nervous system and/or scattered along the digestive tract. However, recent studies have indicated the notochord, the defining feature of the chordate phylum, as an additional site of expression of insulin-like peptides. Here we show that two of the three insulin-like genes identified in Ciona intestinalis, an invertebrate chordate with a dual life cycle, are first expressed in the developing notochord during embryogenesis and transition to distinct areas of the adult digestive tract after metamorphosis. In addition, we present data suggesting that the transcription factor Ciona Brachyury is involved in the control of notochord expression of at least one of these genes, Ciona insulin-like 2. Finally, we review the information currently available on insulin-producing cells in ascidians and on pancreas-related transcription factors that might control their expression.

Ephrin-mediated restriction of ERK1/2 activity delimits the number of pigment cells in the Ciona CNS

Recent evidence suggests that ascidian pigment cells are related to neural crest-derived melanocytes of vertebrates. Using live-imaging, we determine a revised cell lineage of the pigment cells in Ciona intestinalis embryos. The neural precursors undergo successive rounds of anterior-posterior (A-P) oriented cell divisions, starting at the blastula 64-cell stage. A previously unrecognized fourth A-P oriented cell division in the pigment cell lineage leads to the generation of the post-mitotic pigment cell precursors. We provide evidence that MEK/ERK signals are required for pigment cell specification until approximately 30min after the final cell division has taken place. Following each of the four A-P oriented cell divisions, ERK1/2 is differentially activated in the posterior sister cells, into which the pigment cell lineage segregates. Eph/ephrin signals are critical during the third A-P oriented cell division to spatially restrict ERK1/2 activation to the posterior daughter cell. Targeted inhibition of Eph/ephrin signals results in, at neurula stages, anterior expansion of both ERK1/2 activation and a pigment cell lineage marker and subsequently, at larval stages, supernumerary pigment cells. We discuss the implications of these findings with respect to the evolution of the vertebrate neural crest.

p120RasGAP mediates ephrin/Eph-dependent attenuation of FGF/ERK signals during cell fate specification in ascidian embryos

ERK1/2 MAP kinase exhibits a highly dynamic activation pattern in developing embryos, which largely depends on fibroblast growth factor (FGF) signals. In ascidian embryos, FGF-dependent activation of ERK1/2 occurs differentially between sister cells during marginal zone and neural lineage patterning. Selective attenuation of FGF signals by localised ephrin/Eph signals accounts for this differential ERK activation, which controls the binary fate choice of each sibling cell pair. Here, we show that p120 Ras GTPase-activating protein (p120RasGAP) is a crucial mediator of these ephrin/Eph signals. First, inhibition of p120RasGAP has a similar effect to inhibition of ephrin/Eph function during marginal zone and neural patterning. Second, p120RasGAP acts epistatically to ephrin/Eph signals. Third, p120RasGAP physically associates with Eph3 in an ephrin-dependent manner. This study provides the first in vivo evidence that the functional association between Eph and RasGAP controls the spatial extent of FGF-activated ERK.

Ascidians as excellent models for studying cellular events in the chordate body plan

The larvae of non-vertebrate chordate ascidians consist of countable numbers of cells. With this feature, ascidians provide us with excellent models for studying cellular events in the construction of the chordate body. This review discusses the recent observations of morphogenetic movements and cell cycles and divisions along with tissue specifications during ascidian embryogenesis. Unequal cleavages take place at the posterior blastomeres during the early cleavage stages of ascidians, and the structure named the centrosome-attracting body restricts the position of the nuclei near the posterior pole to achieve the unequal cleavages. The most-posterior cells differentiate into the primordial germ cells. The gastrulation of ascidians starts as early as the 110-cell stage. During gastrulation, the endodermal cells show two-step changes in cell shape that are crucial for gastrulation. The ascidian notochord is composed of only 40 cells. The 40 cells align to form a single row by an event named the convergent extension, and then the notochord cells undergo vacuolation to transform the notochord into a single hollowed tube. The strictly restricted number of notochord cells is achieved by the regulated number of cell divisions coupled with the differentiation of the cells conducted by a key transcription factor, Brachyury. The dorsally located neural tube is a characteristic of chordates. During the closure of the ascidian neural tube, the epidermis surrounding the neural plate moves toward the midline to close the neural fold. This morphogenetic movement is allowed by an elongation of interphase in the epidermal cell cycles.

From notochord formation to hereditary chordoma: the many roles of Brachyury

Chordoma is a rare, but often malignant, bone cancer that preferentially affects the axial skeleton and the skull base. These tumors are both sporadic and hereditary and appear to occur more frequently after the fourth decade of life; however, modern technologies have increased the detection of pediatric chordomas. Chordomas originate from remnants of the notochord, the main embryonic axial structure that precedes the backbone, and share with notochord cells both histological features and the expression of characteristic genes. One such gene is Brachyury, which encodes for a sequence-specific transcription factor. Known for decades as a main regulator of notochord formation, Brachyury has recently gained interest as a biomarker and causative agent of chordoma, and therefore as a promising therapeutic target. Here, we review the main characteristics of chordoma, the molecular markers, and the clinical approaches currently available for the early detection and possible treatment of this cancer. In particular, we report on the current knowledge of the role of Brachyury and of its possible mechanisms of action in both notochord formation and chordoma etiogenesis.

FGF signaling induces mesoderm in the hemichordate Saccoglossus kowalevskii

FGFs act in vertebrate mesoderm induction and also play key roles in early mesoderm formation in ascidians and amphioxus. However, in sea urchins initial characterizations of FGF function do not support a role in early mesoderm induction, making the ancestral roles of FGF signaling and mechanisms of mesoderm specification in deuterostomes unclear. In order to better characterize the evolution of mesoderm formation, we have examined the role of FGF signaling during mesoderm development in Saccoglossus kowalevskii, an experimentally tractable representative of hemichordates. We report the expression of an FGF ligand, fgf8/17/18, in ectoderm overlying sites of mesoderm specification within the archenteron endomesoderm. Embryological experiments demonstrate that mesoderm induction in the archenteron requires contact with ectoderm, and loss-of-function experiments indicate that both FGF ligand and receptor are necessary for mesoderm specification. fgf8/17/18 gain-of-function experiments establish that FGF8/17/18 is sufficient to induce mesoderm in adjacent endomesoderm. These experiments suggest that FGF signaling is necessary from the earliest stages of mesoderm specification and is required for all mesoderm development. Furthermore, they suggest that the archenteron is competent to form mesoderm or endoderm, and that FGF signaling from the ectoderm defines the location and amount of mesoderm. When considered in a comparative context, these data support a phylogenetically broad requirement for FGF8/17/18 signaling in mesoderm specification and suggest that FGF signaling played an ancestral role in deuterostome mesoderm formation.

Extending the family table: Insights from beyond vertebrates into the regulation of embryonic development by FGFs

Since the discovery of fibroblast growth factors (FGFs) much focus has been placed on elucidating the roles for each vertebrate FGF ligand, receptor, and regulating molecules in the context of vertebrate development, human disorders and cancer. Studies in human, mouse, frog, chick, and zebrafish have made great contributions to our understanding of the role of FGFs in specific processes. However, in recent years, as more genomes are sequenced, information is becoming available from many non-vertebrate models and a more complete picture of the FGF superfamily as a whole is emerging. In some cases, less redundancy in these FGF signaling systems may allow for more mechanistic insights. Studies in sea anemones have highlighted how ancient FGF signaling is and helped provide insight into the evolution of the FGF gene family. Work in nematodes has shown that different splice forms can be used for functional specificity in invertebrate FGF signaling. Comparing FGFs between urochordates and vertebrates as well as between different insect species reveals important clues into the process of gene loss, duplication and subfunctionalization of FGFs throughout evolution. Finally, comparing all members of the FGF ligand superfamily reveals variability in many properties, which may point to a feature of FGFs as being highly adaptable with regards to protein structure and signaling mechanism. Further studies on FGF signaling outside of vertebrates is likely to continue to complement work in vertebrates by contributing additional insights to the FGF field and providing unexpected information that could be used for medical applications.

FGF signalling: diverse roles during early vertebrate embryogenesis

Fibroblast growth factor (FGF) signalling has been implicated during several phases of early embryogenesis, including the patterning of the embryonic axes, the induction and/or maintenance of several cell lineages and the coordination of morphogenetic movements. Here, we summarise our current understanding of the regulation and roles of FGF signalling during early vertebrate development.

Early zygotic expression of transcription factors and signal molecules in fully dissociated embryonic cells of Ciona intestinalis: A microarray analysis

Specification of early embryonic cells of animals is established by maternally provided factors and interactions of neighboring cells. The present study addressed a question of autonomous versus non-autonomous specification of embryonic cells by using the Ciona intestinalis embryo, in particular the genetic cascade of zygotic expression of transcription factor genes responsible for notochord specification. To examine this issue, we combined the classic experiment of continuous dissociation of embryonic cells with the modern technique of oligonucleotide-based microarrays. We measured early zygotic expression of 389 core transcription factors genes and 118 major signal molecule genes in embryonic cells that were fully dissociated from the first cleavage. Our results indicated that even if cells are free from contact with neighbors, the major transcription factor genes that have primary roles in embryonic cell specification commence their zygotic expression at the same time as in normal embryos. Dissociation of embryonic cells did not affect extracellular signal-regulated kinases (ERK) activity. Although normal embryos treated with U0126 failed to express Bra and Twist-like-1, dissociated embryonic cells treated with U0126 expressed the genes. These results are discussed in relation to the grade of autonomous versus non-autonomous genetic cascades that are responsible for the specification of early Ciona embryonic cells.

Gene-regulatory networks in the Ciona embryos

Ascidians belong to the subphylum Urochordata or Tunicata, which is the sister group of the vertebrates. The simple architecture of the ascidian larva represents the basic chordate body plan. Recent analyses have shown many instances of developmental mechanisms conserved during evolution, while these studies have also revealed a much larger number of instances of divergence. However, to precisely determine the degree of conservation and divergence, that is, how many ways are used to make tadpole-like larvae, we need a systems-level understanding of development. Because animal development is organized by the genome and the minimal functional unit of development is a cell, comprehensiveness and single-cell resolution are necessary for a systems-biological understanding of the development. In the ascidian Ciona intestinalis, gene-regulatory networks responsible for the embryonic development have been studied on a genome-wide scale and at single-cell resolution. The simplicity and compactness of the genome facilitates genome-wide studies. In the Ciona genome, only approximately 670 transcription factor genes are encoded, and their expression profiles during the embryonic development have been analyzed. Gene-knockdown analyses of the transcription factor genes expressed during the embryonic development have been performed. The simplicity of the embryo permits these analyses to be done at single-cell resolution. Actually, these simple embryos are now being modeled in the computer, which allows us to understand the gene-regulatory networks very precisely in three dimensions.

The evolution of cell types in animals: emerging principles from molecular studies

Cell types are fundamental units of multicellular life but their evolution is obscure. How did the first cell types emerge and become distinct in animal evolution? What were the sets of cell types that existed at important evolutionary nodes that represent eumetazoan or bilaterian ancestors? How did these ancient cell types diversify further during the evolution of organ systems in the descending evolutionary lines? The recent advent of cell type molecular fingerprinting has yielded initial insights into the evolutionary interrelationships of cell types between remote animal phyla and has allowed us to define some first principles of cell type diversification in animal evolution.

FGF signalling controls formation of the apical sensory organ in the cnidarian Nematostella vectensis

Fibroblast growth factor (FGF) signalling regulates essential developmental processes in vertebrates and invertebrates, but its role during early metazoan evolution remains obscure. Here, we analyse the function of FGF signalling in a non-bilaterian animal, the sea anemone Nematostella vectensis. We identified the complete set of FGF ligands and FGF receptors, of which two paralogous FGFs (NvFGFa1 and NvFGFa2) and one FGF receptor (NvFGFRa) are specifically coexpressed in the developing apical organ, a sensory structure located at the aboral pole of ciliated larvae from various phyla. Morpholino-mediated knockdown experiments reveal that NvFGFa1 and NvFGFRa are required for the formation of the apical organ, whereas NvFGFa2 counteracts NvFGFRa signalling to prevent precocious and ectopic apical organ development. Marker gene expression analysis shows that FGF signalling regulates local patterning in the aboral region. Furthermore, NvFGFa1 activates its own expression and that of the antagonistic NvFGFa2, thereby establishing positive- and negative-feedback loops. Finally, we show that loss of the apical organ upon NvFGFa1 knockdown blocks metamorphosis into polyps. We propose that the control of the development of sensory structures at the apical pole of ciliated larvae is an ancestral function of FGF signalling.

Ephrin signaling establishes asymmetric cell fates in an endomesoderm lineage of the Ciona embryo

Mesodermal tissues arise from diverse cell lineages and molecular strategies in the Ciona embryo. For example, the notochord and mesenchyme are induced by FGF/MAPK signaling, whereas the tail muscles are specified autonomously by the localized determinant, Macho-1. A unique mesoderm lineage, the trunk lateral cells, develop from a single pair of endomesoderm cells, the A6.3 blastomeres, which form part of the anterior endoderm, hematopoietic mesoderm and muscle derivatives. MAPK signaling is active in the endoderm descendants of A6.3, but is absent from the mesoderm lineage. Inhibition of MAPK signaling results in expanded expression of mesoderm marker genes and loss of endoderm markers, whereas ectopic MAPK activation produces the opposite phenotype: the transformation of mesoderm into endoderm. Evidence is presented that a specific Ephrin signaling molecule, Ci-ephrin-Ad, is required to establish asymmetric MAPK signaling in the endomesoderm. Reducing Ci-ephrin-Ad activity via morpholino injection results in ectopic MAPK signaling and conversion of the mesoderm lineage into endoderm. Conversely, misexpression of Ci-ephrin-Ad in the endoderm induces ectopic activation of mesodermal marker genes. These results extend recent observations regarding the role of Ephrin signaling in the establishment of asymmetric cell fates in the Ciona notochord and neural tube.