Regulatory blueprint for a chordate embryo

Transcription Factors of the bHLH Family Delineate Vertebrate Landmarks in the Nervous System of a Simple Chordate

Tunicates are marine invertebrates whose tadpole-like larvae feature a highly simplified version of the chordate body plan. Similar to their distant vertebrate relatives, tunicate larvae develop a regionalized central nervous system and form distinct neural structures, which include a rostral sensory vesicle, a motor ganglion, and a caudal nerve cord. The sensory vesicle contains a photoreceptive complex and a statocyst, and based on the comparable expression patterns of evolutionarily conserved marker genes, it is believed to include proto-hypothalamic and proto-retinal territories. The evolutionarily conserved molecular fingerprints of these landmarks of the vertebrate brain consist of genes encoding for different transcription factors, and of the gene batteries that they control, and include several members of the bHLH family. Here we review the complement of bHLH genes present in the streamlined genome of the tunicate Ciona robusta and their current classification, and summarize recent studies on proneural bHLH transcription factors and their expression territories. We discuss the possible roles of bHLH genes in establishing the molecular compartmentalization of the enticing nervous system of this unassuming chordate.

The Cis-Regulatory Code for Kelch-like 21/30 Specific Expression in Ciona robusta Sensory Organs

The tunicate Ciona robusta is an emerging model system to study the evolution of the nervous system. Due to their small embryos and compact genomes, tunicates, like Ciona robusta, have great potential to comprehend genetic circuitry underlying cell specific gene repertoire, among different neuronal cells. Their simple larvae possess a sensory vesicle comprising two pigmented sensory organs, the ocellus and the otolith. We focused here on Klhl21/30, a gene belonging to Kelch family, that, in Ciona robusta, starts to be expressed in pigmented cell precursors, becoming specifically maintained in the otolith precursor during embryogenesis. Evolutionary analyses demonstrated the conservation of Klhl21/30 in all the chordates. Cis-regulatory analyses and CRISPR/Cas9 mutagenesis of potential upstream factors, revealed that Klhl21/30 expression is controlled by the combined action of three transcription factors, Mitf, Dmrt, and Msx, which are downstream of FGF signaling. The central role of Mitf is consistent with its function as a fundamental regulator of vertebrate pigment cell development. Moreover, our results unraveled a new function for Dmrt and Msx as transcriptional co-activators in the context of the Ciona otolith.

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.

Fgf-driven Tbx protein activities directly induce myf5 and myod to initiate zebrafish myogenesis

Skeletal muscle derives from dorsal mesoderm formed during vertebrate gastrulation. Fibroblast growth factor (Fgf) signalling cooperates with Tbx transcription factors to promote dorsal mesoderm formation, but their role in myogenesis has been unclear. Using zebrafish, we show that dorsally derived Fgf signals act through Tbx16 and Tbxta to induce slow and fast trunk muscle precursors at distinct dorsoventral positions. Tbx16 binds to and directly activates the myf5 and myod genes, which are required for commitment to myogenesis. Tbx16 activity depends on Fgf signalling from the organiser. In contrast, Tbxta is not required for myf5 expression, but binds a specific site upstream of myod that is not bound by Tbx16 and drives (dependent on Fgf signals) myod expression in adaxial slow precursors, thereby initiating trunk myogenesis. After gastrulation, when similar muscle cell populations in the post-anal tail are generated from tailbud, declining Fgf signalling is less effective at initiating adaxial myogenesis, which is instead initiated by Hedgehog signalling from the notochord. Our findings suggest a hypothesis for ancestral vertebrate trunk myogenic patterning and how it was co-opted during tail evolution to generate similar muscle by new mechanisms.

This article has an associated ‘The people behind the papers’ interview.

Highlighted Article: Tbx16 and Tbxta activate myf5 and myod directly during the earliest myogenesis in zebrafish, and Fgf signalling acts through Tbx16 to drive myogenesis in trunk but not tail.

MorphoSeq: Full Single-Cell Transcriptome Dynamics Up to Gastrulation in a Chordate

Single-cell RNA sequencing (scRNA-seq) provides a leap forward in resolving cellular diversity and developmental trajectories but fails to comprehensively delineate the spatial organization and precise cellular makeup of individual embryos. Here, we reconstruct from scRNA-seq and light sheet imaging data a canonical digital embryo that captures the genome-wide gene expression trajectory of every single cell at every cell division in the 18 lineages up to gastrulation in the ascidian Phallusia mammillata. By using high-coverage scRNA-seq, we devise a computational framework that stratifies single cells of individual embryos into cell types without prior knowledge. Unbiased transcriptome data analysis mapped each cell’s physical position and lineage history, yielding the complete history of gene expression at the genome-wide level for every single cell in a developing embryo. A comparison of individual embryos reveals both extensive reproducibility between symmetric embryo sides and a large inter-embryonic variability due to small differences in embryogenesis timing.

•

Integration of scRNA-seq and imaging data yield a canonical digital embryo

•

Cell type classification without prior knowledge

•

De novo reconstruction of the lineage history and spatial organization of the embryo

•

Timing differences contribute to inter-embryo variability in gene expression

Identification of an LPS-Induced Chemo-Attractive Peptide from Ciona robusta

Background: Previously published work has demonstrated that the LPS injection of Ciona robusta leads to the overexpression of a truncated form of an immune-related mRNA (C8short) by means of Ciona robusta (CR) alternative polyadenylation (APA) (CR-APA). Methods: The 3D structure of the C8short-derived Ciona robusta chemo-attractive peptide (CrCP) was evaluated by homology modeling. The biological activity of the CrCP was studied in vitro using a primary human dermal cell line (HuDe). Real-Time PCR was used to investigate the expression levels of genes involved in cell motility. NF-κB signaling was studied by western blotting. Results: In silico modeling showed that CrCP displayed structural characteristics already reported for a short domain of the vertebrate CRK gene, suggesting its possible involvement in cell migration mechanisms. In vitro assays demonstrated that CrCP was capable of inducing the motility of HuDe cells in both wound healing and chemo-attractive experiments. qPCR demonstrated the capability of CrCP to modulate the expression of the matrix metalloproteinase-7 (MMP-7) and E-cadherin genes. Finally, western blot analysis demonstrated that treatment with CrCP induced activation of the NF-κB signaling pathway. Conclusion: Our results describe the characterization of the 3D structure and chemo-attractive activity of an LPS-induced CrCP peptide from Ciona robusta.

High-Throughput Protein Production Combined with High- Throughput SELEX Identifies an Extensive Atlas of Ciona robusta Transcription Factor DNA-Binding Specificities

Transcription factors (TFs) control gene transcription, binding to specific DNA motifs located in cis-regulatory elements across the genome. The identification of TF-binding motifs is thus an important aspect to understand the role of TFs in gene regulation. SELEX, Systematic Evolution of Ligands by EXponential enrichment, is an efficient in vitro method, which can be used to determine the DNA-binding specificity of TFs. Thanks to the development of high-throughput (HT) DNA cloning system and protein production technology, the classical SELEX assay has be extended to high-throughput scale (HT-SELEX).We report here the detailed protocol for the cloning, production, and purification of 420 Ciona robusta DNA BD. 263 Ciona robusta TF DNA-binding domain proteins were purified in milligram quantities and analyzed by HT-SELEX. The identification of 139 recognition sequences generates an atlas of protein-DNA-binding specificities that is crucial for the understanding of the gene regulatory network (GRN) of Ciona robusta. Overall, our analysis suggests that the Ciona robusta repertoire of sequence-specific transcription factors comprises less than 500 genes. The protocols for high-throughput protein production and HT-SELEX described in this article for the study of Ciona robusta TF DNA-binding specificity are generic and have been successfully applied to a wide range of TFs from other species, including human, mouse, and Drosophila.

A Nearly Complete Genome of Ciona intestinalis Type A (C. robusta) Reveals the Contribution of Inversion to Chromosomal Evolution in the Genus Ciona

Since its initial publication in 2002, the genome of Ciona intestinalis type A (Ciona robusta), the first genome sequence of an invertebrate chordate, has provided a valuable resource for a wide range of biological studies, including developmental biology, evolutionary biology, and neuroscience. The genome assembly was updated in 2008, and it included 68% of the sequence information in 14 pairs of chromosomes. However, a more contiguous genome is required for analyses of higher order genomic structure and of chromosomal evolution. Here, we provide a new genome assembly for an inbred line of this animal, constructed with short and long sequencing reads and Hi-C data. In this latest assembly, over 95% of the 123 Mb of sequence data was included in the chromosomes. Short sequencing reads predicted a genome size of 114-120 Mb; therefore, it is likely that the current assembly contains almost the entire genome, although this estimate of genome size was smaller than previous estimates. Remapping of the Hi-C data onto the new assembly revealed a large inversion in the genome of the inbred line. Moreover, a comparison of this genome assembly with that of Ciona savignyi, a different species in the same genus, revealed many chromosomal inversions between these two Ciona species, suggesting that such inversions have occurred frequently and have contributed to chromosomal evolution of Ciona species. Thus, the present assembly greatly improves an essential resource for genome-wide studies of ascidians.

The genetic program to specify ectodermal cells in ascidian embryos

The ascidian belongs to the sister group of vertebrates and shares many features with them. The gene regulatory network (GRN) controlling gene expression in ascidian embryonic development leading to the tadpole larva has revealed evolutionarily conserved gene circuits between ascidians and vertebrates. These conserved mechanisms are indeed useful to infer the original developmental programs of the ancestral chordates. Simultaneously, these studies have revealed which gene circuits are missing in the ascidian GRN; these gene circuits may have been acquired in the vertebrate lineage. In particular, the GRN responsible for gene expression in ectodermal cells of ascidian embryos has revealed the genetic programs that regulate the regionalization of the brain, formation of palps derived from placode-like cells, and differentiation of sensory neurons derived from neural crest-like cells. We here discuss how these studies have given insights into the evolution of these traits.

Finding cell-specific expression patterns in the early Ciona embryo with single-cell RNA-seq

Single-cell RNA-seq has been established as a reliable and accessible technique enabling new types of analyses, such as identifying cell types and studying spatial and temporal gene expression variation and change at single-cell resolution. Recently, single-cell RNA-seq has been applied to developing embryos, which offers great potential for finding and characterising genes controlling the course of development along with their expression patterns. In this study, we applied single-cell RNA-seq to the 16-cell stage of the Ciona embryo, a marine chordate and performed a computational search for cell-specific gene expression patterns. We recovered many known expression patterns from our single-cell RNA-seq data and despite extensive previous screens, we succeeded in finding new cell-specific patterns, which we validated by in situ and single-cell qPCR.

Exploring miR-9 Involvement in Ciona intestinalis Neural Development Using Peptide Nucleic Acids

The microRNAs are small RNAs that regulate gene expression at the post-transcriptional level and can be involved in the onset of neurodegenerative diseases and cancer. They are emerging as possible targets for antisense-based therapy, even though the in vivo stability of miRNA analogues is still questioned. We tested the ability of peptide nucleic acids, a novel class of nucleic acid mimics, to downregulate miR-9 in vivo in an invertebrate model organism, the ascidian Ciona intestinalis, by microinjection of antisense molecules in the eggs. It is known that miR-9 is a well-conserved microRNA in bilaterians and we found that it is expressed in epidermal sensory neurons of the tail in the larva of C. intestinalis. Larvae developed from injected eggs showed a reduced differentiation of tail neurons, confirming the possibility to use peptide nucleic acid PNA to downregulate miRNA in a whole organism. By identifying putative targets of miR-9, we discuss the role of this miRNA in the development of the peripheral nervous system of ascidians.

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.

ANISEED 2019: 4D exploration of genetic data for an extended range of tunicates

ANISEED (https://www.aniseed.cnrs.fr) is the main model organism database for the worldwide community of scientists working on tunicates, the vertebrate sister-group. Information provided for each species includes functionally-annotated gene and transcript models with orthology relationships within tunicates, and with echinoderms, cephalochordates and vertebrates. Beyond genes the system describes other genetic elements, including repeated elements and cis-regulatory modules. Gene expression profiles for several thousand genes are formalized in both wild-type and experimentally-manipulated conditions, using formal anatomical ontologies. These data can be explored through three complementary types of browsers, each offering a different view-point. A developmental browser summarizes the information in a gene- or territory-centric manner. Advanced genomic browsers integrate the genetic features surrounding genes or gene sets within a species. A Genomicus synteny browser explores the conservation of local gene order across deuterostome. This new release covers an extended taxonomic range of 14 species, including for the first time a non-ascidian species, the appendicularian Oikopleura dioica. Functional annotations, provided for each species, were enhanced through a combination of manual curation of gene models and the development of an improved orthology detection pipeline. Finally, gene expression profiles and anatomical territories can be explored in 4D online through the newly developed Morphonet morphogenetic browser.

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.

Foxg specifies sensory neurons in the anterior neural plate border of the ascidian embryo

Foxg constitutes a regulatory loop with Fgf8 and plays an important role in the development of anterior placodes and the telencephalon in vertebrate embryos. Ascidians, which belong to Tunicata, the sister group of vertebrates, develop a primitive placode-like structure at the anterior boundary of the neural plate, but lack a clear counterpart of the telencephalon. In this animal, Foxg is expressed in larval palps, which are adhesive organs with sensory neurons. Here, we show that Foxg begins to be expressed in two separate rows of cells within the neural plate boundary region under the control of the MAPK pathway to pattern this region. However, Foxg is not expressed in the brain, and we find no evidence that knockdown of Foxg affects brain formation. Our data suggest that recruitment of Fgf to the downstream of Foxg might have been a critical evolutionary event for the telencephalon in the vertebrate lineage.

Vertebrate telencephalon formation requires Foxg-Fgf8 cross-regulation, but while ascidians express Foxg in the neural plate, they lack a telencephalon. Here the authors show that Foxg loss does not affect ascidian brain formation, indicating that telencephalon evolution required recruitment of Fgf downstream of Foxg.

Comprehensive single-cell transcriptome lineages of a proto-vertebrate

Ascidian embryos highlight the importance of cell lineages in animal development. As simple proto-vertebrates they also provide insights into the evolutionary origins of novel cell types, such as cranial placodes and neural crest. To build upon these efforts we have determined single cell transcriptomes for more than 90,000 cells spanning the entirety of Ciona intestinalis development, from the onset of gastrulation to swimming tadpoles. This represents an average of over 12-fold coverage for every cell at every stage of development, owing to the small cell numbers of ascidian embryos. Single cell transcriptome trajectories were used to construct “virtual” cell lineage maps and provisional gene networks for nearly 40 different neuronal subtypes comprising the larval nervous system. We summarize several applications of these datasets, including annotating the synaptome of swimming tadpoles and tracing the evolutionary origin of novel cell types such as the vertebrate telencephalon.

Regulation and evolution of muscle development in tunicates

For more than a century, studies on tunicate muscle formation have revealed many principles of cell fate specification, gene regulation, morphogenesis, and evolution. Here, we review the key studies that have probed the development of all the various muscle cell types in a wide variety of tunicate species. We seize this occasion to explore the implications and questions raised by these findings in the broader context of muscle evolution in chordates.

Evolution of Snail-mediated regulation of neural crest and placodes from an ancient role in bilaterian neurogenesis

A major challenge in vertebrate evolution is to identify the gene regulatory mechanisms that facilitated the origin of neural crest cells and placodes from ancestral precursors in invertebrates. Here, we show in lamprey, a primitively jawless vertebrate, that the transcription factor Snail is expressed simultaneously throughout the neural plate, neural plate border, and pre-placodal ectoderm in the early embryo and is then upregulated in the CNS throughout neurogenesis. Using CRISPR/Cas9-mediated genome editing, we demonstrate that Snail plays functional roles in all of these embryonic domains or their derivatives. We first show that Snail patterns the neural plate border by repressing lateral expansion of Pax3/7 and activating nMyc and ZicA. We also present evidence that Snail is essential for DlxB-mediated establishment of the pre-placodal ectoderm but is not required for SoxB1a expression during formation of the neural plate proper. At later stages, Snail regulates formation of neural crest-derived and placode-derived PNS neurons and controls CNS neural differentiation in part by promoting cell survival. Taken together with established functions of invertebrate Snail genes, we identify a pan-bilaterian mechanism that extends to jawless vertebrates for regulating neurogenesis that is dependent on Snail transcription factors. We propose that ancestral vertebrates deployed an evolutionarily conserved Snail expression domain in the CNS and PNS for neurogenesis and then acquired derived functions in neural crest and placode development by recruitment of regulatory genes downstream of neuroectodermal Snail activity. Our results suggest that Snail regulatory mechanisms in vertebrate novelties such as the neural crest and placodes may have emerged from neurogenic roles that originated early in bilaterian evolution.

Antero-posterior ectoderm patterning by canonical Wnt signaling during ascidian development

Wnt/β-catenin signaling is an ancient pathway in metazoans and controls various developmental processes, in particular the establishment and patterning of the embryonic primary axis. In vertebrates, a graded Wnt activity from posterior to anterior endows cells with positional information in the central nervous system. Recent studies in hemichordates support a conserved role for Wnt/β-catenin in ectoderm antero-posterior patterning at the base of the deuterostomes. Ascidians are marine invertebrates and the closest relatives of vertebrates. By combining gain- and loss-of-function approaches, we have determined the role of Wnt/β-catenin in patterning the three ectoderm derivatives of the ascidian Ciona intestinalis, central nervous system, peripheral nervous system and epidermis. Activating Wnt/β-catenin signaling from gastrulation led to a dramatic transformation of the ectoderm with a loss of anterior identities and a reciprocal anterior extension of posterior identities, consistent with studies in other metazoans. Surprisingly, inhibiting Wnt signaling did not produce a reciprocal anteriorization of the embryo with a loss of more posterior identities like in vertebrates and hemichordate. Epidermis patterning was overall unchanged. Only the identity of two discrete regions of the central nervous system, the anteriormost and the posteriormost regions, were under the control of Wnt. Finally, the caudal peripheral nervous system, while being initially Wnt dependent, formed normally. Our results show that the Ciona embryonic ectoderm responds to Wnt activation in a manner that is compatible with the proposed function for this pathway at the base of the deuterostomes. However, possibly because of its fast and divergent mode of development that includes extensive use of maternal determinants, the overall antero-posterior patterning of the Ciona ectoderm is Wnt independent, and Wnt/β-catenin signaling controls the formation of some sub-domains. Our results thus indicate that there has likely been a drift in the developmental systems controlling ectoderm patterning in the lineage leading to ascidians.

The Wnt/β-catenin pathway is a system of cell-cell communication. It has an ancient origin in animals and plays multiple roles during embryogenesis and adult life. In particular, it is involved in determining, in the vertebrate embryo, the identity of the different parts of the body and their relative positions along the antero-posterior axis. We have investigated in an ascidian (or sea squirt) species, a marine invertebrate that is closely related to vertebrates, whether this pathway had a similar role. Like in vertebrates, activating Wnt/β-catenin led to a posteriorization of the embryo with a loss of anterior structures. By contrast, unlike vertebrates, ascidian embryos formed rather normally following Wnt/β-catenin inactivation. Since hemichordates (or acorn worms), earlier divergent invertebrates, use Wnt/β-catenin in a manner comparable to vertebrates, it is in the ascidian lineage that changes have occurred. Consequently, ascidians build an antero-posterior axis, very similarly organized to that of vertebrates, but in a different way.

The regulatory pathway from genes directly activated by maternal factors to muscle structural genes in ascidian embryos

Striated muscle cells in the tail of ascidian tadpole larvae differentiate cell-autonomously. Although several key regulatory factors have been identified, the genetic regulatory pathway is not fully understood; comprehensive understanding of the regulatory pathway is essential for accurate modeling in order to deduce principles for gene regulatory network dynamics, and for comparative analysis on how ascidians have evolved the cell-autonomous gene regulatory mechanism. Here, we reveal regulatory interactions among three key regulatory factors, Zic-r.b, Tbx6-r.b and Mrf, and elucidate the mechanism by which these factors activate muscle structural genes. We reveal a cross-regulatory circuit among these regulatory factors, which maintains the expression of Tbx6-r.b and Mrf during gastrulation. Although these two factors combinatorially activate muscle structural genes in late-stage embryos, muscle structural genes are activated mainly by Tbx6-r.b before gastrulation. Time points when expression of muscle structural genes become first detectable are strongly correlated with the degree of Tbx6-r.b occupancy. Thus, the genetic pathway, starting with Tbx6-r.b and Zic-r.b, which are activated by maternal factors, and ending with expression of muscle structural genes, has been revealed.

Positioning a multifunctional basic helix-loop-helix transcription factor within the Ciona notochord gene regulatory network

In a multitude of organisms, transcription factors of the basic helix-loop-helix (bHLH) family control the expression of genes required for organ development and tissue differentiation. The functions of different bHLH transcription factors in the specification of nervous system and paraxial mesoderm have been widely investigated in various model systems. Conversely, the knowledge of the role of these regulators in the development of the axial mesoderm, the embryonic territory that gives rise to the notochord, and the identities of their target genes, remain still fragmentary. Here we investigated the transcriptional regulation and target genes of Bhlh-tun1, a bHLH transcription factor expressed in the developing Ciona notochord as well as in additional embryonic territories that contribute to the formation of both larval and adult structures. We describe its possible role in notochord formation, its relationship with the key notochord transcription factor Brachyury, and suggest molecular mechanisms through which Bhlh-tun1 controls the spatial and temporal expression of its effectors.

Evolution of embryonic cis-regulatory landscapes between divergent Phallusia and Ciona ascidians

Ascidian species of the Phallusia and Ciona genera are distantly related, their last common ancestor dating several hundred million years ago. Although their genome sequences have extensively diverged since this radiation, Phallusia and Ciona species share almost identical early morphogenesis and stereotyped cell lineages. Here, we explored the evolution of transcriptional control between P. mammillata and C. robusta. We combined genome-wide mapping of open chromatin regions in both species with a comparative analysis of the regulatory sequences of a test set of 10 pairs of orthologous early regulatory genes with conserved expression patterns. We find that ascidian chromatin accessibility landscapes obey similar rules as in other metazoa. Open-chromatin regions are short, highly conserved within each genus and cluster around regulatory genes. The dynamics of chromatin accessibility and closest-gene expression are strongly correlated during early embryogenesis. Open-chromatin regions are highly enriched in cis-regulatory elements: 73% of 49 open chromatin regions around our test genes behaved as either distal enhancers or proximal enhancer/promoters following electroporation in Phallusia eggs. Analysis of this datasets suggests a pervasive use in ascidians of "shadow" enhancers with partially overlapping activities. Cross-species electroporations point to a deep conservation of both the trans-regulatory logic between these distantly-related ascidians and the cis-regulatory activities of individual enhancers. Finally, we found that the relative order and approximate distance to the transcription start site of open chromatin regions can be conserved between Ciona and Phallusia species despite extensive sequence divergence, a property that can be used to identify orthologous enhancers, whose regulatory activity can partially diverge.

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.

Enhancer activities of amphioxus Brachyury genes in embryos of the ascidian, Ciona intestinalis

The notochord and somites are distinctive chordate structures. The T-box transcription factor gene, Brachyury, is expressed in notochord and plays a pivotal role in its formation. In the cephalochordate, Branchiostoma floridae, Brachyury is duplicated into BfBra1 and BfBra2, which are expressed in the somite-formation region as well. In a series of experiments to elucidate the regulatory machinery of chordate Brachyury expression, we carried out a lacZ reporter assay of BfBra in embryos of the urochordate, Ciona intestinalis. Vista analyses suggest the presence of conserved non-coding sequences, not only in the 5'-upstream, but also in the 3'-downstream and in introns of BfBra. We found that: (1) 5'-upstream sequences of both BfBra1 and BfBra2 promote lacZ expression in muscle cells, (2) 3'-downstream sequences have enhancer activity that promotes lacZ expression in notochord cells, and (3) introns of BfBra2 and BfBra1 exhibit lacZ expression preferentially in muscle and notochord cells. These results suggest shared cephalochordate Brachyury enhancer machinery that also works in urochordates. We discuss the results in relation to evolutionary modification of Brachyury expression in formation of chordate-specific organs characteristic of each lineage.

Evolution of the bilaterian mouth and anus

It is widely held that the bilaterian tubular gut with mouth and anus evolved from a simple gut with one major gastric opening. However, there is no consensus on how this happened. Did the single gastric opening evolve into a mouth, with the anus forming elsewhere in the body (protostomy), or did it evolve into an anus, with the mouth forming elsewhere (deuterostomy), or did it evolve into both mouth and anus (amphistomy)? These questions are addressed by the comparison of developmental fates of the blastopore, the opening of the embryonic gut, in diverse animals that live today. Here we review comparative data on the identity and fate of blastoporal tissue, investigate how the formation of the through-gut relates to the major body axes, and discuss to what extent evolutionary scenarios are consistent with these data. Available evidence indicates that stem bilaterians had a slit-like gastric opening that was partially closed in subsequent evolution, leaving open the anus and most likely also the mouth, which would favour amphistomy. We discuss remaining difficulties, and outline directions for future research.

Outstanding questions in developmental ERK signaling

The extracellular signal-regulated kinase (ERK) pathway leads to activation of the effector molecule ERK, which controls downstream responses by phosphorylating a variety of substrates, including transcription factors. Crucial insights into the regulation and function of this pathway came from studying embryos in which specific phenotypes arise from aberrant ERK activation. Despite decades of research, several important questions remain to be addressed for deeper understanding of this highly conserved signaling system and its function. Answering these questions will require quantifying the first steps of pathway activation, elucidating the mechanisms of transcriptional interpretation and measuring the quantitative limits of ERK signaling within which the system must operate to avoid developmental defects.

Ascidian Zic Genes

Ascidians are tunicates, which constitute the sister group of vertebrates. The ascidian genome contains two Zic genes, called Zic-r.a (also called Macho-1) and Zic-r.b (ZicL). The latter is a multi-copy gene, and the precise copy number has not yet been determined. Zic-r.a is maternally expressed, and soon after fertilization Zic-r.a mRNA is localized in the posterior pole of the zygote. Zic-r.a protein is translated there and is involved in specification of posterior fate; in particular it is important for specification of muscle fate. Zic-r.a is also expressed zygotically in neural cells of the tailbud stage. On the other hand, Zic-r.b is first expressed in marginal cells of the vegetal hemisphere of 32-cell embryos and then in neural cells that contribute to the central nervous system during gastrulation. Zic-r.b is required first for specification of mesodermal tissues and then for specification of the central nervous system. Their upstream and downstream genetic pathways have been studied extensively by functional assays, which include gene knockdown and chromatin immunoprecipitation assays. Thus, ascidian Zic genes play central roles in specification of mesodermal and neural fates.

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.

Distinct regulation of Snail in two muscle lineages of the ascidian embryo achieves temporal coordination of muscle development

The transcriptional repressor Snail is required for proper differentiation of the tail muscle of ascidian tadpole larvae. Two muscle lineages (B5.1 and B6.4) contribute to the anterior tail muscle cells, and are consecutively separated from a transcriptionally quiescent germ cell lineage at the 16- and 32-cell stages. Concomitantly, cells of these lineages begin to express Tbx6.b (Tbx6-r.b) at the 16- and 32-cell stages, respectively. Meanwhile, Snail expression begins in these two lineages simultaneously at the 32-cell stage. Here, we show that Snail expression is regulated differently between these two lineages. In the B5.1 lineage, Snail was activated through Tbx6.b, which is activated by maternal factors, including Zic-r.a. In the B6.4 lineage, the MAPK pathway was cell-autonomously activated by a constitutively active form of Raf, enabling Zic-r.a to activate Snail independently of Tbx6.b As a result, Snail begins to be expressed at the 32-cell stage simultaneously in these two lineages. Such shortcuts might be required for coordinating developmental programs in embryos in which cells become separated progressively from stem cells, including germline cells.

Depletion of Maternal Cyclin B3 Contributes to Zygotic Genome Activation in the Ciona Embryo

Most animal embryos display a delay in the activation of zygotic transcription during early embryogenesis [1]. This process is thought to help coordinate rapid increases in cell number during early development [2]. The timing of zygotic genome activation (ZGA) during the maternal-to-zygotic transition (MZT) remains uncertain despite extensive efforts. We explore ZGA in the simple protovertebrate, Ciona intestinalis. Single-cell RNA sequencing (RNA-seq) assays identified Cyclin B3 (Ccnb3) as a putative mediator of ZGA. Maternal Ccnb3 transcripts rapidly diminish in abundance during the onset of zygotic transcription at the 8-cell and 16-cell stages. Disruption of Ccnb3 activity results in precocious activation of zygotic transcription, while overexpression abolishes normal activation. These observations suggest that the depletion of maternal Cyclin B3 products is a critical component of the MZT and ZGA. We discuss evidence that this mechanism might play a conserved role in the MZT of other metazoans, including mice and humans.

Reporter Analyses Reveal Redundant Enhancers that Confer Robustness on Cis-Regulatory Mechanisms

Reporter analyses of Hox1 and Brachyury (Bra) genes have revealed examples of redundant enhancers that provide regulatory robustness. Retinoic acid (RA) activates through an RA-response element the transcription of Hox1 in the nerve cord of the ascidian Ciona intestinalis. We also found a weak RA-independent neural enhancer within the second intron of Hox1. The Hox1 gene in the larvacean Oikopleura dioica is also expressed in the nerve cord. The O. dioica genome, however, does not contain the RA receptor-encoding gene, and the expression of Hox1 has become independent of RA. We have found that the upstream sequence of the O. dioica Hox1 was able to activate reporter gene expression in the nerve cord of the C. intestinalis embryo, suggesting that an RA-independent regulatory system in the nerve cord might be common in larvaceans and ascidians. This RA-independent redundant regulatory system may have facilitated the Oikopleura ancestor losing RA signaling without an apparent impact on Hox1 expression domains. On the other hand, vertebrate Bra is expressed in the ventral mesoderm and notochord, whereas its ascidian ortholog is exclusively expressed in the notochord. Fibroblast growth factor (FGF) induces Bra in the ventral mesoderm in vertebrates, whereas it induces Bra in the notochord in ascidians. Disruption of the FGF signal does not completely silence Bra expression in ascidians, suggesting that FGF-dependent and independent enhancers might comprise a redundant regulatory system in ascidians. The existence of redundant enhancers, therefore, provides regulatory robustness that may facilitate the acquisition of new expression domains.

Dynamics of two key maternal factors that initiate zygotic regulatory programs in ascidian embryos

In animal embryos, transcription is repressed for a definite period of time after fertilization. In the embryo of the ascidian, Ciona intestinalis (type A; or Ciona robusta), transcription of regulatory genes is repressed before the 8- or 16-cell stages. This initial transcriptional quiescence is important to enable the establishment of initial differential gene expression patterns along the animal-vegetal axis by maternal factors, because the third cell division separates the animal and vegetal hemispheres into distinct blastomeres. Indeed, maternal transcription factors directly activate zygotic gene expression by the 16-cell stage; Tcf7/β-catenin activates genes in the vegetal hemisphere, and Gata.a activates genes in the animal hemisphere. In the present study, we revealed the dynamics of Gata.a and β-catenin, and expression profiles of their target genes precisely. β-catenin began to translocate into the nuclei at the 16-cell stage, and thus expression of β-catenin targets began at the 16-cell stage. Although Gata.a is abundantly present before the 8-cell stage, transcription of Gata.a targets was repressed at and before the 4-cell stage, and their expression began at the 8-cell stage. Transcription of the β-catenin targets may be repressed by the same mechanism in early embryos, because β-catenin targets were not expressed in 4-cell embryos treated with a GSK inhibitor, in which β-catenin translocated to the nuclei. Thus, these two maternal factors have different dynamics, which establish the pre-pattern for zygotic genetic programs in 16-cell embryos.

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.

Formation of adult organs through metamorphosis in ascidians

The representative characteristic of ascidians is their vertebrate-like, tadpole shape at the larval stage. Ascidians lose the tadpole shape through metamorphosis to become adults with a nonmotile, sessile body and a shape generally considered distinct from that of vertebrates. Solitary ascidians including Ciona species are extensively studied to understand the developmental mechanisms of ascidians, and to compare these mechanisms with their counterparts in vertebrates. In these ascidian species, the digestive and circulatory systems are not well developed in the larval trunk and the larvae do not take food. This is in contrast with the inner conditions of vertebrate tadpoles, which have functional organs comparable to those of adults. The adult organs and tissues of these ascidians become functional during metamorphosis that is completed quickly, suggesting that the ascidian larvae of solitary species are a transient stage of development. We here discuss how the cells and tissues in the ascidian larval body are converted into those of adults. The hearts of ascidians and vertebrates use closely related cellular and molecular mechanisms that suggest their shared origin. Hox genes of ascidians are essential for forming adult endodermal structures. To fully understand the development and evolution of chordates, a complete elucidation of the mechanisms underlying the adult tissue/organ formation of ascidians will be needed. WIREs Dev Biol 2018, 7:e304. doi: 10.1002/wdev.304 This article is categorized under: Comparative Development and Evolution > Body Plan Evolution Early Embryonic Development > Development to the Basic Body Plan.

Organ specific gene expression in the regenerating tail of Macrostomum lignano

Temporal and spatial characterization of gene expression is a prerequisite for the understanding of cell-, tissue-, and organ-differentiation. In a multifaceted approach to investigate gene expression in the tail plate of the free-living marine flatworm Macrostomum lignano, we performed a posterior-region-specific in situ hybridization screen, RNA sequencing (RNA-seq) of regenerating animals, and functional analyses of selected tail-specific genes. The in situ screen revealed transcripts expressed in the antrum, cement glands, adhesive organs, prostate glands, rhabdite glands, and other tissues. Next we used RNA-seq to characterize temporal expression in the regenerating tail plate revealing a time restricted onset of both adhesive organs and copulatory apparatus regeneration. In addition, we identified three novel previously unannotated genes solely expressed in the regenerating stylet. RNA interference showed that these genes are required for the formation of not only the stylet but the whole male copulatory apparatus. RNAi treated animals lacked the stylet, vesicula granulorum, seminal vesicle, false seminal vesicle, and prostate glands, while the other tissues of the tail plate, such as adhesive organs regenerated normally. In summary, our findings provide a large resource of expression data during homeostasis and regeneration of the morphologically complex tail regeneration and pave the way for a better understanding of organogenesis in M. lignano.

Conserved gene regulatory module specifies lateral neural borders across bilaterians

The lateral neural plate border (NPB), the neural part of the vertebrate neural border, is composed of central nervous system (CNS) progenitors and peripheral nervous system (PNS) progenitors. In invertebrates, PNS progenitors are also juxtaposed to the lateral boundary of the CNS. Whether there are conserved molecular mechanisms determining vertebrate and invertebrate lateral neural borders remains unclear. Using single-cell-resolution gene-expression profiling and genetic analysis, we present evidence that orthologs of the NPB specification module specify the invertebrate lateral neural border, which is composed of CNS and PNS progenitors. First, like in vertebrates, the conserved neuroectoderm lateral border specifier Msx/vab-15 specifies lateral neuroblasts in Caenorhabditis elegans Second, orthologs of the vertebrate NPB specification module (Msx/vab-15, Pax3/7/pax-3, and Zic/ref-2) are significantly enriched in worm lateral neuroblasts. In addition, like in other bilaterians, the expression domain of Msx/vab-15 is more lateral than those of Pax3/7/pax-3 and Zic/ref-2 in C. elegans Third, we show that Msx/vab-15 regulates the development of mechanosensory neurons derived from lateral neural progenitors in multiple invertebrate species, including C. elegans, Drosophila melanogaster, and Ciona intestinalis We also identify a novel lateral neural border specifier, ZNF703/tlp-1, which functions synergistically with Msx/vab-15 in both C. elegans and Xenopus laevis These data suggest a common origin of the molecular mechanism specifying lateral neural borders across bilaterians.

Expected Number of Fixed Points in Boolean Networks with Arbitrary Topology

Boolean network models describe genetic, neural, and social dynamics in complex networks, where the dynamics depend generally on network topology. Fixed points in a genetic regulatory network are typically considered to correspond to cell types in an organism. We prove that the expected number of fixed points in a Boolean network, with Boolean functions drawn from probability distributions that are not required to be uniform or identical, is one, and is independent of network topology if only a feedback arc set satisfies a stochastic neutrality condition. We also demonstrate that the expected number is increased by the predominance of positive feedback in a cycle.

Differential gene expression along the animal-vegetal axis in the ascidian embryo is maintained by a dual functional protein Foxd

In many animal embryos, a specific gene expression pattern is established along the animal-vegetal axis soon after zygotic transcription begins. In the embryo of the ascidian Ciona intestinalis, soon after the division that separates animal and vegetal hemispheres into distinct blastomeres, maternal Gata.a and β-catenin activate specific genes in the animal and vegetal blastomeres, respectively. On the basis of these initial distinct gene expression patterns, gene regulatory networks promote animal cells to become ectodermal tissues and vegetal cells to become endomesodermal tissues and a part of the nerve cord. In the vegetal hemisphere, β-catenin directly activates Foxd, an essential transcription factor gene for specifying endomesodermal fates. In the present study, we found that Foxd also represses the expression of genes that are activated specifically in the animal hemisphere, including Dmrt1, Prdm1-r.a (Bz1), Prdm1-r.b (Bz2), and Otx. A reporter assay showed that Dmrt1 expression was directly repressed by Foxd, and a chromatin immunoprecipitation assay showed that Foxd was bound to the upstream regions of Dmrt1, Prdm1-r.a, Prdm1-r.b, and Otx. Thus, Foxd has a dual function of activating specific gene expression in the vegetal hemisphere and of repressing the expression of genes that are normally expressed in the animal hemisphere. This dual function stabilizes the initial patterning along the animal-vegetal axis by β-catenin and Gata.a.

In embryogenesis of most animals, a specific gene expression pattern is established along the animal-vegetal axis first. In the embryo of the ascidian Ciona intestinalis, the activity of the maternal factor Gata.a is suppressed by β-catenin, which is active only in the vegetal hemisphere, and thereby these two factors activate specific genes in the animal and vegetal blastomeres, respectively. We found that a gene encoding a transcription factor, Foxd, which is a direct target of β-catenin, works as a promoter for endomesodermal fate and an inhibitor for ectodermal fate. In the ascidian embryo, the animal-vegetal axis initially established by the maternal factors is not stable enough for subsequent developmental processes, and needs to be maintained by Foxd. Thus, the animal hemisphere fate is suppressed first by the maternal factor β-catenin, and then by Foxd, which is activated by β-catenin. The primary embryonic axis is not stable initially, and stabilized by a transcription factor, which is expressed differentially along the axis.

How complexity increases in development: An analysis of the spatial-temporal dynamics of Gene expression in Ciona intestinalis

The increase in complexity in an embryo over developmental time is perhaps one of the most intuitive processes of animal development. It is also intuitive that the embryo becomes progressively compartmentalized over time and space. In spite of this intuitiveness, there are no systematic attempts to quantify how this occurs. Here, we present a quantitative analysis of the compartmentalization and spatial complexity of Ciona intestinalis over developmental time by analyzing thousands of gene expression spatial patterns from the ANISEED database. We measure compartmentalization in two ways: as the relative volume of expression of genes and as the disparity in gene expression between body parts. We also use a measure of the curvature of each gene expression pattern in 3D space. These measures show a similar increase over time, with the most dramatic change occurring from the 112-cell stage to the early tailbud stage. Combined, these measures point to a global pattern of increase in complexity in the Ciona embryo. Finally, we cluster the different regions of the embryo depending on their gene expression similarity, within and between stages. Results from this clustering analysis, which partially correspond to known fate maps, provide a global quantitative overview about differentiation and compartmentalization between body parts at each developmental stage.

Patterning of brain precursors in ascidian embryos

In terms of their embryonic origins, the anterior and posterior parts of the ascidian central nervous system (CNS) are associated with distinct germ layers. The anterior part of the sensory vesicle, or brain, originates from ectoderm lineages following a neuro-epidermal binary fate decision. In contrast, a large part of the remaining posterior CNS is generated following neuro-mesodermal binary fate decisions. Here, we address the mechanisms that pattern the anterior brain precursors along the medial-lateral axis (future ventral-dorsal) at neural plate stages. Our functional studies show that Nodal signals are required for induction of lateral genes, including Delta-like, Snail, Msxb and Trp Delta-like/Notch signalling induces intermediate (Gsx) over medial (Meis) gene expression in intermediate cells, whereas the combinatorial action of Snail and Msxb prevents the expression of Gsx in lateral cells. We conclude that despite the distinct embryonic lineage origins within the larval CNS, the mechanisms that pattern neural precursors are remarkably similar.

Tfap2 and Sox1/2/3 cooperatively specify ectodermal fates in ascidian embryos

Epidermis and neural tissues differentiate from the ectoderm in animal embryos. Although epidermal fate is thought to be induced in vertebrate embryos, embryological evidence has indicated that no intercellular interactions during early stages are required for epidermal fate in ascidian embryos. To test this hypothesis, we determined the gene regulatory circuits for epidermal and neural specification in the ascidian embryo. These circuits started with Tfap2-r.b and Sox1/2/3, which are expressed in the ectodermal lineage immediately after zygotic genome activation. Tfap2-r.b expression was diminished in the neural lineages upon activation of fibroblast growth factor signaling, which is known to induce neural fate, and sustained only in the epidermal lineage. Tfap2-r.b specified the epidermal fate cooperatively with Dlx.b, which was activated by Sox1/2/3 This Sox1/2/3-Dlx.b circuit was also required for specification of the anterior neural fate. In the posterior neural lineage, Sox1/2/3 activated Nodal, which is required for specification of the posterior neural fate. Our findings support the hypothesis that the epidermal fate is specified autonomously in ascidian embryos.