Ascidian embryonic cells with properties of neural-crest cells and neuromesodermal progenitors of vertebrates.

Neural-crest cells and neuromesodermal progenitors (NMPs) are multipotent cells that are important for development of vertebrate embryos. In embryos of ascidians, which are the closest invertebrate relatives of vertebrates, several cells located at the border between the neural plate and the epidermal region have neural-crest-like properties; hence, the last common ancestor of ascidians and vertebrates may have had ancestral cells similar to neural-crest cells. However, these ascidian neural-crest-like cells do not produce cells that are commonly of mesodermal origin. Here we showed that a cell population located in the lateral region of the neural plate has properties resembling those of vertebrate neural-crest cells and NMPs. Among them, cells with Tbx6-related expression contribute to muscle near the tip of the tail region and cells with Sox1/2/3 expression give rise to the nerve cord. These observations and cross-species transcriptome comparisons indicate that these cells have properties similar to those of NMPs. Meanwhile, transcription factor genes Dlx.b, Zic-r.b and Snai, which are reminiscent of a gene circuit in vertebrate neural-crest cells, are involved in activation of Tbx6-related.b. Thus, the last common ancestor of ascidians and vertebrates may have had cells with properties of neural-crest cells and NMPs and such ancestral cells may have produced cells commonly of ectodermal and mesodermal origins.

A digital twin reproducing gene regulatory network dynamics of early Ciona embryos indicates robust buffers in the network.

How gene regulatory networks (GRNs) encode gene expression dynamics and how GRNs evolve are not well understood, although these problems have been studied extensively. We created a digital twin that accurately reproduces expression dynamics of 13 genes that initiate expression in 32-cell ascidian embryos. We first showed that gene expression patterns can be manipulated according to predictions by this digital model. Next, to simulate GRN rewiring, we changed regulatory functions that represented their regulatory mechanisms in the digital twin, and found that in 55 of 100 cases, removal of a single regulator from a conjunctive clause of Boolean functions did not theoretically alter qualitative expression patterns of these genes. In other words, we found that more than half the regulators gave theoretically redundant temporal or spatial information to target genes. We experimentally substantiated that the expression pattern of Nodal was maintained without one of these factors, Zfpm, by changing the upstream regulatory sequence of Nodal. Such robust buffers of regulatory mechanisms may provide a basis of enabling developmental system drift, or rewiring of GRNs without changing expression patterns of downstream genes, during evolution.

BMP signaling is required to form the anterior neural plate border in ascidian embryos.

Cranial neurogenic placodes have been considered vertebrate innovations. However, anterior neural plate border (ANB) cells of ascidian embryos share many properties with vertebrate neurogenic placodes; therefore, it is now believed that the last common ancestor of vertebrates and ascidians had embryonic structures similar to neurogenic placodes of vertebrate embryos. Because BMP signaling is important for specifying the placode region in vertebrate embryos, we examined whether BMP signaling is also involved in gene expression in the ANB region of ascidian embryos. Our data indicated that Admp, a divergent BMP family member, is mainly responsible for BMP signaling in the ANB region, and that two BMP-antagonists, Noggin and Chordin, restrict the domain, in which BMP signaling is activated, to the ANB region, and prevent it from expanding to the neural plate. BMP signaling is required for expression of Foxg and Six1/2 at the late gastrula stage, and also for expression of Zf220, which encodes a zinc finger transcription factor in late neurula embryos. Because Zf220 negatively regulates Foxg, when we downregulated Zf220 by inhibiting BMP signaling, Foxg was upregulated, resulting in one large palp instead of three palps (adhesive organs derived from ANB cells). Functions of BMP signaling in specification of the ANB region give further support to the hypothesis that ascidian ANB cells share an evolutionary origin with vertebrate cranial placodes.

Regulators specifying cell fate activate cell cycle regulator genes to determine cell numbers in ascidian larval tissues.

In animal development, most cell types stop dividing before terminal differentiation; thus, cell cycle control is tightly linked to cell differentiation programmes. In ascidian embryos, cell lineages do not vary among individuals, and rounds of the cell cycle are determined according to cell lineages. Notochord and muscle cells stop dividing after eight or nine rounds of cell division depending on their lineages. In the present study, we showed that a Cdk inhibitor, Cdkn1.b, is responsible for stopping cell cycle progression in these lineages. Cdkn1.b is also necessary for epidermal cells to stop dividing. In contrast, mesenchymal and endodermal cells continue to divide even after hatching, and Myc is responsible for maintaining cell cycle progression in these tissues. Expression of Cdkn1.b in notochord and muscle is controlled by transcription factors that specify the developmental fate of notochord and muscle. Likewise, expression of Myc in mesenchyme and endoderm is under control of transcription factors that specify the developmental fate of mesenchyme and endoderm. Thus, cell fate specification and cell cycle control are linked by these transcription factors.

Tracing the evolutionary history of blood cells to the unicellular ancestor of animals.

Blood cells are thought to have emerged as phagocytes in the common ancestor of animals followed by the appearance of novel blood cell lineages such as thrombocytes, erythrocytes, and lymphocytes, during evolution. However, this speculation is not based on genetic evidence and it is still possible to argue that phagocytes in different species have different origins. It also remains to be clarified how the initial blood cells evolved; whether ancient animals have solely developed de novo programs for phagocytes or they have inherited a key program from ancestral unicellular organisms. Here, we traced the evolutionary history of blood cells, and cross-species comparison of gene expression profiles revealed that phagocytes in various animal species and Capsaspora (C.) owczarzaki, a unicellular organism, are transcriptionally similar to each other. We also found that both phagocytes and C. owczarzaki share a common phagocytic program, and that CEBPα is the sole transcription factor highly expressed in both phagocytes and C. owczarzaki. We further showed that the function of CEBPα to drive phagocyte program in nonphagocytic blood cells has been conserved in tunicate, sponge, and C. owczarzaki. We finally showed that, in murine hematopoiesis, repression of CEBPα to maintain nonphagocytic lineages is commonly achieved by polycomb complexes. These findings indicate that the initial blood cells emerged inheriting a unicellular organism program driven by CEBPα and that the program has also been seamlessly inherited in phagocytes of various animal species throughout evolution.

A Manually Curated Gene Model Set for an Ascidian, Ciona robusta (Ciona intestinalis Type A).

Gene/transcript model sets predicted from decoded genome sequences are an important resource for a wide range of biological studies. Accuracy of gene models is therefore critical for deducing accurate conclusions. Computationally predicted models are sometimes inconsistent with experimental data from cDNA clones and RNA-sequencing. In an ascidian, Ciona robusta (Ciona intestinalis type A), a manually curated gene/transcript model set, which was constructed using an assembly in which 68% of decoded sequences were associated with chromosomes, had been used during the last decade. Recently a new genome assembly was published, in which over 95% of decoded sequences are associated with chromosomes. In the present study, we provide a high-quality version of the gene/transcript model set for the latest assembly. Because the Ciona genome has been used in a variety of studies such as developmental biological studies, evolutionary studies, and physiological studies, the current gene/transcript model set provides a fundamental biological resource.

Protein kinases and protein phosphatases encoded in the Ciona robusta genome.

Protein kinases (PKs) and protein phosphatases (PPs) regulate the phosphorylation of proteins that are involved in a variety of biological processes. To study such biological processes systematically, it is important to know the whole repertoire of PKs and PPs encoded in a genome. In the present study, we surveyed the genome of an ascidian (Ciona robusta or Ciona intestinalis type A) to comprehensively identify the genes that encoded PKs and PPs. Because ascidians belong to the sister group of vertebrates, a comparison of the whole repertoire of PKs and PPs of ascidians with those of vertebrates may help to delineate the complements of these proteins that were present in the last common ancestor of these two groups of animals. Our results show that the repertory of PPs was much more expanded in vertebrates than the repertory of PKs. We also showed that approximately 75% of PKs and PPs were expressed during development from eggs to larvae. Thus, the present study provides catalogs for PKs and PPs encoded in the ascidian genome. These catalogs will be useful for systematic studies of biological processes that involve phosphorylation and for evolutionary studies of the origin of vertebrates.

Tbx15/18/22 shares a binding site with Tbx6-r.b to maintain expression of a muscle structural gene in ascidian late embryos.

The ascidian larval tail contains muscle cells for swimming. Most of these muscle cells differentiate autonomously. The genetic program behind this autonomy has been studied extensively and the genetic cascade from maternal factors to initiation of expression of a muscle structural gene, Myl.c, has been uncovered; Myl.c expression is directed initially by transcription factor Tbx6-r.b at the 64-cell stage and then by the combined actions of Tbx6-r.b and Mrf from the gastrula to early tailbud stages. In the present study, we showed that transcription of Myl.c continued in late tailbud embryos and larvae, although a fusion protein of Tbx6-r.b and GFP was hardly detectable in late tailbud embryos. A knockdown experiment, reporter assay, and in vitro binding assay indicated that an essential cis-regulatory element of Myl.c that bound Tbx6-r.b in early embryos bound Tbx15/18/22 in late embryos to maintain expression of Myl.c. We also found that Tbx15/18/22 was controlled by Mrf, which constitutes a regulatory loop with Tbx6-r.b. Therefore, our data indicated that Tbx15/18/22 was activated initially under control of this regulatory loop as in the case of Myl.c, and then Tbx15/18/22 maintained the expression of Myl.c after Tbx6-r.b had disappeared. RNA-sequencing of Tbx15/18/22 morphant embryos revealed that many muscle structural genes were regulated similarly by Tbx15/18/22. Thus, the present study revealed the mechanisms of maintenance of transcription of muscle structural genes in late embryos in which Tbx15/18/22 takes the place of Tbx6-r.b.

Cell geometry, signal dampening, and a bimodal transcriptional response underlie the spatial precision of an ERK-mediated embryonic induction.

Precise control of lineage segregation is critical for the development of multicellular organisms, but our quantitative understanding of how variable signaling inputs are integrated to activate lineage-specific gene programs remains limited. Here, we show how precisely two out of eight ectoderm cells adopt neural fates in response to ephrin and FGF signals during ascidian neural induction. In each ectoderm cell, FGF signals activate ERK to a level that mirrors its cell contact surface with FGF-expressing mesendoderm cells. This gradual interpretation of FGF inputs is followed by a bimodal transcriptional response of the immediate early gene, Otx, resulting in its activation specifically in the neural precursors. At low levels of ERK, Otx is repressed by an ETS family transcriptional repressor, ERF2. Ephrin signals are critical for dampening ERK activation levels across ectoderm cells so that only neural precursors exhibit above-threshold levels, evade ERF repression, and "switch on" Otx transcription.

The gene regulatory system for specifying germ layers in early embryos of the simple chordate.

In animal embryos, gene regulatory networks control the dynamics of gene expression in cells and coordinate such dynamics among cells. In ascidian embryos, gene expression dynamics have been dissected at the single-cell resolution. Here, we revealed mathematical functions that represent the regulatory logics of all regulatory genes expressed at the 32-cell stage when the germ layers are largely specified. These functions collectively explain the entire mechanism by which gene expression dynamics are controlled coordinately in early embryos. We found that regulatory functions for genes expressed in each of the specific lineages contain a common core regulatory mechanism. Last, we showed that the expression of the regulatory genes became reproducible by calculation and controllable by experimental manipulations. Thus, these regulatory functions represent an architectural design for the germ layer specification of this chordate and provide a platform for simulations and experiments to understand the operating principles of gene regulatory networks.

Two distinct motifs for Zic-r.a drive specific gene expression in two cell lineages.

Zic-r.a, a maternal transcription factor, specifies posterior fate in ascidian embryos. However, its direct target, Tbx6-r.b, does not contain typical Zic-r.a-binding sites in its regulatory region. Using an in vitro selection assay, we found that Zic-r.a binds to sites dissimilar to the canonical motif, by which it activates Tbx6-r.b in a sub-lineage of muscle cells. These sites with non-canonical motifs have weak affinity for Zic-r.a; therefore, it activates Tbx6-r.b only in cells expressing Zic-r.a abundantly. Meanwhile, we found that Zic-r.a expressed zygotically in late embryos activates neural genes through canonical sites. Because different zinc-finger domains of Zic-r.a are important for driving reporters with canonical and non-canonical sites, it is likely that the non-canonical motif is not a divergent version of the canonical motif. In other words, our data indicate that the non-canonical motif represents a motif distinct from the canonical motif. Thus, Zic-r.a recognizes two distinct motifs to activate two sets of genes at two timepoints in development. This article has an associated 'The people behind the papers' interview.

Chromosomal Inversion Polymorphisms in Two Sympatric Ascidian Lineages.

Chromosomal rearrangements can reduce fitness of heterozygotes and can thereby prevent gene flow. Therefore, such rearrangements can play a role in local adaptation and speciation. In particular, inversions are considered to be a major potential cause for chromosomal speciation. There are two closely related, partially sympatric lineages of ascidians in the genus Ciona, which we call type-A and type-B animals in the present study. Although these invertebrate chordates are largely isolated reproductively, hybrids can be found in wild populations, suggesting incomplete prezygotic barriers. Although the genome of type-A animals has been decoded and widely used, the genome for type-B animals has not been decoded at the chromosomal level. In the present study, we sequenced the genomes of two type-B individuals from different sides of the English Channel (in the zone of sympatry with type-A individuals) and compared them at the chromosomal level with the type-A genome. Although the overall structures were well conserved between type A and type B, chromosomal alignments revealed many inversions differentiating these two types of Ciona; it is probable that the frequent inversions have contributed to separation between these two lineages. In addition, comparisons of the genomes between the two type-B individuals revealed that type B had high rates of inversion polymorphisms and nucleotide polymorphisms, and thus type B might be in the process of differentiation into multiple new types or species. Our results suggest an important role of inversions in chromosomal speciation of these broadcasting spawners.

Cis-regulatory code for determining the action of Foxd as both an activator and a repressor in ascidian embryos.

In early embryos of Ciona, an invertebrate chordate, the animal-vegetal axis is established by the combinatorial actions of maternal factors. One target of these maternal factors, Foxd, is specifically expressed in the vegetal hemisphere and stabilizes the animal-vegetal axis by activating vegetal hemisphere-specific genes and repressing animal hemisphere-specific genes. This dual functionality is essential for the embryogenesis of early ascidian embryos; however, the mechanism by which Foxd can act as both a repressor and an activator is unknown. Here, we identify a Foxd binding site upstream of Lhx3/4, which is activated by Foxd, and compare it with a repressive Foxd binding site upstream of Dmrt.a. We found that activating sites bind Foxd with low affinity while repressive sites bind Foxd with high affinity. Reporter assays confirm that this qualitative difference between activating and repressive Foxd binding sites is sufficient to change Foxd functionality. We therefore conclude that the outcome of Foxd transcriptional regulation is encoded in cis-regulatory elements.

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.

A single-cell analysis of the molecular lineage of chordate embryogenesis.

Progressive unfolding of gene expression cascades underlies diverse embryonic lineage development. Here, we report a single-cell RNA sequencing analysis of the complete and invariant embryonic cell lineage of the tunicate Ciona savignyi from fertilization to the onset of gastrulation. We reconstructed a developmental landscape of 47 cell types over eight cell cycles in the wild-type embryo and identified eight fate transformations upon fibroblast growth factor (FGF) inhibition. For most FGF-dependent asymmetric cell divisions, the bipotent mother cell displays the gene signature of the default daughter fate. In convergent differentiation of the two notochord lineages, we identified additional gene pathways parallel to the master regulator T/Brachyury Last, we showed that the defined Ciona cell types can be matched to E6.5-E8.5 stage mouse cell types and display conserved expression of limited number of transcription factors. This study provides a high-resolution single-cell dataset to understand chordate early embryogenesis and cell lineage differentiation.

A gene regulatory network for cell fate specification in Ciona embryos.

Ascidian embryos are used as a model system in developmental biology due to their unique properties, including their invariant cell division patterns, being comprised of a small number of cells and tissues, the feasibility of their experimental manipulation, and their simple and compact genome. These properties have provided an opportunity for examining the gene regulatory network at the single cell resolution and at a genome-wide scale. This article summarizes when and where each regulatory gene is expressed in early ascidian embryos, and the extent to which the gene regulatory network explains each gene expression.

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.

Gata is ubiquitously required for the earliest zygotic gene transcription in the ascidian embryo.

In ascidian embryos, the earliest transcription from the zygotic genome begins between the 8-cell and 16-cell stages. Gata.a, a maternally expressed Gata transcription factor, activates target genes specifically in the animal hemisphere, whereas the complex of ß-catenin and Tcf7 antagonizes the activity of Gata.a and activates target genes specifically in the vegetal hemisphere. Here, we show that genes zygotically expressed at the 16-cell stage have significantly more Gata motifs in their upstream regions. These genes included not only genes with animal hemisphere-specific expression but also genes with vegetal hemisphere-specific expression. On the basis of this finding, we performed knockdown experiments for Gata.a and reporter assays, and found that Gata.a is required for the expression of not only genes with animal hemisphere-specific expression, but also genes with vegetal hemisphere-specific expression. Our data indicated that weak Gata.a activity that cannot induce animal hemisphere-specific expression can allow ß-catenin/Tcf7 targets to be expressed in the vegetal cells. Because genes zygotically expressed at the 32-cell stage also had significantly more Gata motifs in their upstream regions, Gata.a function may not be limited to the genes expressed specifically in the animal or vegetal hemispheres at the 16-cell stage, and Gata.a may play an important role in the earliest transcription of the zygotic genome.

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.

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.