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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

Initiation of the zygotic genetic program in the ascidian embryo.

In embryos of ascidians, which are invertebrate chordates, the expression of a small number of genes from the zygotic genome begins between the 8- and 16-cell stages. They are considered the first zygotic genes to be expressed under the direct control of maternal factors. The initial transcriptional quiescence before the 8-cell stage is essential for establishing differential gene expression patterns between the animal and vegetal hemispheres, because these hemispheres are first segregated into distinct blastomeres at the 8-cell stage. While the mechanisms of this transcriptional silencing have not been understood, the mechanism by which maternal factors cooperatively establish the first differential gene expression patterns has been well understood. Here, we review the initial processes that occur until the 16-cell stage in ascidian embryos.

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.

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.

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.

ANISEED 2017: extending the integrated ascidian database to the exploration and evolutionary comparison of genome-scale datasets.

ANISEED (www.aniseed.cnrs.fr) is the main model organism database for tunicates, the sister-group of vertebrates. This release gives access to annotated genomes, gene expression patterns, and anatomical descriptions for nine ascidian species. It provides increased integration with external molecular and taxonomy databases, better support for epigenomics datasets, in particular RNA-seq, ChIP-seq and SELEX-seq, and features novel interactive interfaces for existing and novel datatypes. In particular, the cross-species navigation and comparison is enhanced through a novel taxonomy section describing each represented species and through the implementation of interactive phylogenetic gene trees for 60% of tunicate genes. The gene expression section displays the results of RNA-seq experiments for the three major model species of solitary ascidians. Gene expression is controlled by the binding of transcription factors to cis-regulatory sequences. A high-resolution description of the DNA-binding specificity for 131 Ciona robusta (formerly C. intestinalis type A) transcription factors by SELEX-seq is provided and used to map candidate binding sites across the Ciona robusta and Phallusia mammillata genomes. Finally, use of a WashU Epigenome browser enhances genome navigation, while a Genomicus server was set up to explore microsynteny relationships within tunicates and with vertebrates, Amphioxus, echinoderms and hemichordates.

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.

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.

Antagonism between ß-catenin and Gata.a sequentially segregates the germ layers of ascidian embryos.

Many animal embryos use nuclear ß-catenin (nß-catenin) during the segregation of endomesoderm (or endoderm) from ectoderm. This mechanism is thus likely to be evolutionarily ancient. In the ascidian embryo, nß-catenin reiteratively drives binary fate decisions between ectoderm and endomesoderm at the 16-cell stage, and then between endoderm and margin (mesoderm and caudal neural) at the 32-cell stage. At the 16-cell stage, nß-catenin activates endomesoderm genes in the vegetal hemisphere. At the same time, nß-catenin suppresses the DNA-binding activity of a maternal transcription factor, Gata.a, through a physical interaction, and Gata.a thereby activates its target genes only in the ectodermal lineage. In the present study, we found that this antagonism between nß-catenin and Gata.a also operates during the binary fate switch at the 32-cell stage. Namely, in marginal cells where nß-catenin is absent, Gata.a directly activates its target, Zic-r.b (ZicL), to specify the marginal cell lineages. Thus, the antagonistic action between nß-catenin and Gata.a is involved in two consecutive stages of germ layer segregation in ascidian embryos.