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.

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.

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.

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.

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.

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.

FGF9/16/20 and Wnt-5alpha signals are involved in specification of secondary muscle fate in embryos of the ascidian, Halocynthia roretzi.

The tail muscle cells of the ascidian tadpole larva originate from two different lineages, the B- (primary) and A- and b- (secondary) line blastomeres of the eight-cell stage embryo. The primary muscle cells assume muscle fate cell-autonomously with the involvement of a localized muscle determinant, macho-1. On the other hand, fate determination of secondary muscle cells is a non-cell-autonomous process that depends on cellular interactions. In this paper, we investigated the mechanisms underlying fate specification of secondary muscle cells in Halocynthia roretzi. We found that FGF and Wnt5 signals were required. In contrast, the Nodal signal, which is required for specification of A-line muscle cells in another ascidian, Ciona intestinalis, was not necessary for the formation of any secondary muscle cells in Halocynthia embryo. Therefore, Halocynthia and Ciona show distinctly different mechanisms for generation of the secondary lineages, despite the fact that embryogenesis appears very similar between these species. We also found that the mechanisms involved in specification of A- and b-line muscle cells were distinct in that the required timing of the FGF signal for the A-line muscle cells preceded that for the b-line. Moreover, the inducer blastomeres for specification of these two lineages were different.

A bHLH transcription factor gene, Twist-like 1, is essential for the formation of mesodermal tissues of Ciona juveniles.

Ascidian larval mesenchyme cells, comprising about 900 cells, are derived from the A7.6, B8.5 and B7.7 blastomere pairs in the 110-cell embryo. Previous studies showed that the properties of mesenchyme cells are not uniform among the three lines in embryos of Ciona savignyi and Ciona intestinalis. After metamorphosis, the larval mesenchyme cells form the mesodermal tissues or organs of the adult body. In the present study, the developmental fates of A7.6-, B8.5- and B7.7-line mesenchyme cells were traced using DiI to determine the origins of juvenile mesodermal tissues of C. savignyi. It was demonstrated that each of the A7.6-, B8.5- and B7.7-line mesenchyme cells is distributed in different positions of the larval trunk, and then give rise to the different mesodermal tissues of juveniles. Twist-like 1 is a transcription factor gene essential for the specification of larval mesenchyme cells. Knockdown of this gene with specific morpholino antisense oligonucleotides affected not only the specification of larval mesenchyme cells, but also the formation of most of the mesodermal tissues of juveniles. The juvenile mesodermal tissues in the Twist-like 1-knockdown specimen were never compensated by the surrounding tissues. The present results therefore indicate that Twist-like 1 is required for the differentiation of most mesodermal precursors of adults.

Three distinct lineages of mesenchymal cells in Ciona intestinalis embryos demonstrated by specific gene expression.

The ascidian embryonic mesenchyme, comprising about 900 cells, forms mesodermal tissues or organs of the adult body after metamorphosis. The mesenchyme originates from the A7.6 [trunk lateral cells (TLCs)], B7.7, and B8.5 blastomeres of the 110-cell stage embryo. Previous studies showed that FGF9/16/20 is required for specification of the mesenchyme in Ciona embryos and that two different (A7.6 and B8.5/B7.7) but partially overlapping molecular mechanisms are associated with the expression of a basic helix-loop-helix (bHLH) transcription factor gene, Twist-like1, in the mesenchymal precursors, which triggers the differentiation process of mesenchyme cells. In the present study, we examined whether the three embryonic lineages express the same mesenchyme-specific structural genes under the control of a common mechanism or whether the three lineages are characterized by the expression of genes specific to each of the lineages. We characterized nine mesenchyme-specific genes in Ciona embryos and found that five were expressed in A7.6/B8.5/B7.7, two in B8.5/B7.7, and two in B7.7 only. FGF9/16/20 and Twist-like1 were required for the expression of all the mesenchyme-specific genes, except for three A7.6/B8.5/B7.7-specific genes in A7.6 progenitors. Overexpression of FGF9/16/20 or Twist-like1 upregulated the expression of A7.6/B8.5/B7.7- and B8.5/B7.7-specific genes, while it downregulated the expression of B7.7-specific genes. These results provide evidence that the differentiation of each of the three mesenchyme lineages of Ciona embryos is characterized by the expression of a specific set of genes, whose expression is controlled differentially.

A genomewide survey of developmentally relevant genes in Ciona intestinalis. II. Genes for homeobox transcription factors.

Homeobox-containing genes play crucial roles in various developmental processes, including body-plan specification, pattern formation and cell-type specification. The present study searched the draft genome sequence and cDNA/EST database of the basal chordate Ciona intestinalis to identify 83 homeobox-containing genes in this animal. This number of homeobox genes in the Ciona genome is smaller than that in the Caenorhabditis elegans, Drosophila melanogaster, human and mouse genomes. Of the 83 genes, 76 have possible human orthologues and 7 may be unique to Ciona. The ascidian homeobox genes were classified into 11 classes, including Hox class, NK class, Paired class, POU class, LIM class, TALE class, SIX class, Prox class, Cut class, ZFH class and HNF1 class, according to the classification scheme devised for known homeobox genes. As to the Hox cluster, the Ciona genome contains single copies of each of the paralogous groups, suggesting that there is a single Hox cluster, if any, but genes orthologous to Hox7, 8, 9 and 11 were not found in the genome. In addition, loss of genes had occurred independently in the Ciona lineage and was noticed in Gbx of the EHGbox subclass, Sax, NK3, Vax and vent of the NK class, Cart, Og9, Anf and Mix of the Paired class, POU-I, III, V and VI of the POU class, Lhx6/7 of the LIM class, TGIF of the TALE class, Cux and SATB of the Cut class, and ZFH1 of the ZFH class, which might have reduced the number of Ciona homeobox genes. Interestingly, one of the newly identified Ciona intestinalis genes and its vertebrate counterparts constitute a novel subclass of HNF1 class homeobox genes. Furthermore, evidence for the gene structures and expression of 54 of the 83 homeobox genes was provided by analysis of ESTs, suggesting that cDNAs for these 54 genes are available. The present data thus reveal the repertoire of homeodomain-containing transcription factors in the Ciona genome, which will be useful for future research on the development and evolution of chordates.

A genomewide survey of developmentally relevant genes in Ciona intestinalis. VIII. Genes for PI3K signaling and cell cycle.

Cell growth and cell divisions are two fundamental biological processes for cells in multi-cellular organisms. The molecules involved in these biological processes are highly conserved within eukaryotes, including plants and unicellular organisms such as yeast. However, some regulatory molecules seem to be innovated during animal evolution. Therefore, to understand how the ubiquitous systems have evolved or have been conserved, we examined genes for the phosphoinositide 3-kinase (PI3K) pathway that is important for cell growth, and genes for cell cycle regulation in the genome of Ciona intestinalis. It was found that the Ciona intestinalis genome contains all the essential constituents of the PI3K pathway. In addition, the class IB PI3K catalytic and regulatory subunits, which had not previously been known in animals other than mammals, were found in the Ciona genome. Similarly, all essential cyclins and CDKs were found in the Ciona genome, while cyclin G and cyclin L were likely to be independently lost in the ascidian lineage, which may be dispensable for the cell cycle. Cyclin F, which was previously known only in vertebrates, was not found in the Ciona genome. Therefore, this gene was probably innovated during the evolution of vertebrates to be involved in vertebrate-specific cell cycle regulation. Since Ciona is regarded as one of the most primitive extant chordates, the present analysis gives us an insight into how these fundamental biological genes are evolved or are conserved during chordate evolution.

In situ screening for genes expressed preferentially in secondary mesenchyme cells of sea urchin embryos.

In sea urchin embryos, four types of non-skeletogenic mesodermal cells are derived from secondary mesenchyme cells (SMCs). Although determining the complete lineage of SMCs is currently a high-priority goal, specific markers for each type of SMC-derived cell in Hemicentrotus pulcherrimus are unavailable. To identify genes preferentially expressed in the various types of SMC-derived cells, we constructed a cDNA library of the archenteron isolated from late gastrulae. Both the 5' and 3' ends of 1,050 cDNAs randomly selected from 7,500 picked clones were sequenced. Based on the sequence at the 3' end, the cDNAs were grouped into 671 independent clusters. Among these, 605 clusters were analysed by whole-mount in situ hybridisation; 28% (170 clusters) exhibited differential expression patterns, while 24% were ubiquitously expressed and 48% did not show any staining. Of 170 clusters showing differential expression patterns, 33 clusters were differentially expressed in SMC-derived cells. From these clusters, several genes were obtained that were specifically or predominantly expressed in each type of SMCs, including coelomic pouch cells in which specific expression patterns have not been reported previously, and hence will be useful for lineage studies. Furthermore, in situ hybridisation revealed the existence of a new type or subpopulation of SMCs distributed sparsely in the blastocoel.

Specification and differentiation processes of secondary mesenchyme-derived cells in embryos of the sea urchin Hemicentrotus pulcherrimus.

Four types of mesoderm cells (pigment cells, blastocoelar cells, coelomic pouch cells and circumesophageal muscle cells) are derived from secondary mesenchyme cells (SMC) in sea urchin embryos. To gain information on the specification and differentiation processes of SMC-derived cells, we studied the exact number and division cycles of each type of cell in Hemicentrotus pulcherrimus. Numbers of blastocoelar cells, coelomic pouch cells and circumesophageal muscle fibers were 18.0 +/- 2.0 (36 h post-fertilization (h.p.f.)), 23.0 +/- 2.5 (36 h.p.f.) and 9.5 +/- 1.3 (60 h.p.f.), respectively, whereas the number of pigment cells ranged from 40 to 60. From the diameters of blastocoelar cells and coelomic pouch cells, the numbers of division cycles were elucidated; these two types of cells had undertaken 11 rounds of cell division by the prism stage, somewhat earlier than pigment cells. To determine the relationship among the four types of cells, we tried to alter the number of pigment cells with chemical treatment and found that CH3COONa increased pigment cells without affecting embryo morphology. Interestingly, the number of blastocoelar cells became smaller in CH3COONa-treated embryos. In contrast, blastocoelar cells were markedly increased with NiCl2 treatment, whereas the number of pigment cells was markedly decreased. The number of coelomic pouch cells and circumesophageal muscle fibers was not affected with these treatments, indicating that coelomic pouch and muscle cells are specified independently of, or at much later stages, than pigment and blastocoelar cells.