Regulatory blueprint for a chordate embryo

Transcriptomes of a fast-developing chordate uncover drastic differences in transcription factors and localized maternal RNA composition compared with those of ascidians

The larvacean Oikopleura dioica is a fast-developing chordate because of its small number of cells (∼4500 in juveniles) and rapid development to complete morphogenesis by 10 h after fertilization. Strikingly, most of its blastomeres are restricted to give rise to a single cell-type by the 32-cell stage of embryogenesis, unlike cell fate determination at the 110-cell stage in ascidians. In this study, RNA-sequencing (RNA-seq) revealed non-canonical properties of O. dioica: (1) an initial zygotic gene expression of 950 genes at the 16- to 32-cell stage; (2) 25 transcription factors (TFs) are expressed in the 32-cell stage (fewer than half of the TFs underlying gene regulatory networks in ascidian embryogenesis were lost or not expressed); (3) five maternal mRNAs localized in the vegetal-posterior blastomeres in animal and vegetal hemispheres; and (4) three maternal mRNAs localized in the small vegetal pole region of unfertilized eggs. These observations indicate that this fast-developing chordate lacks the first phase of development in ascidians: fertilization-driven ooplasmic movements that drive postplasmic RNAs toward the vegetal pole. These data have been deposited in ANISEED (https://www.aniseed.fr/) as transcriptome resources.

A master regulatory loop that activates genes in a temporally coordinated manner in muscle cells of ascidian embryos

Ascidian larval muscle cells present a classic example of autonomous development. A regulatory mechanism for these cells has been extensively investigated, and the regulatory gene circuit has been documented from maternal factors to a muscle-specific gene. In the present study, we comprehensively identified genes expressed specifically in ascidian muscle cells, and found that all of them are under control of a positive regulatory loop of Tbx6-r.b and Mrf, the core circuit identified previously. We also found that several transcription factors under control of the Tbx6-r.b/Mrf regulatory loop exhibited various temporal expression profiles, which are probably important for creating functional muscle cells. These results, together with results of previous studies, provide an exhaustive view of the regulatory system enabling autonomous development of ascidian larval muscle cells. It shows that the Tbx6-r.b/Mrf regulatory loop, but not a single gene, serves a 'master' regulatory function. This master regulatory loop not only controls spatial gene expression patterns, but also governs temporal expression patterns in ascidian muscle cells.

A simple inland culture system provides insights into ascidian post-embryonic developmental physiology

Maintenance and breeding of experimental organisms are fundamental to life sciences, but both initial and running costs, and hands-on zootechnical demands can be challenging for many laboratories. Here, we first aimed to further develop a simple protocol for reliable inland culture of tunicate model species of the Ciona genus. We cultured both Ciona robusta and Ciona intestinalis in controlled experimental conditions, with a focus on dietary variables, and quantified growth and maturation parameters. From statistical analysis of these standardized datasets, we gained insights into the post-embryonic developmental physiology of Ciona and inferred an improved diet and culturing conditions for sexual maturation. We showed that body length is a critical determinant of both somatic and sexual maturation, which suggests the existence of systemic control mechanisms of resource allocation towards somatic growth or maturation and supports applying size selection as a predictor of reproductive fitness in our inland culture to keep the healthiest animals at low density in the system. In the end, we successfully established a new protocol, including size selection, to promote both sperm and egg production. Our protocol using small tanks will empower researchers to initiate inland Ciona cultures with low costs and reduced space constraints.

Fine-Tuned Expression of Evolutionarily Conserved Signaling Molecules in the Ciona Notochord

The notochord is an axial structure required for the development of all chordate embryos, from sea squirts to humans. Over the course of more than half a billion years of chordate evolution, in addition to its structural function, the notochord has acquired increasingly relevant patterning roles for its surrounding tissues. This process has involved the co-option of signaling pathways and the acquisition of novel molecular mechanisms responsible for the precise timing and modalities of their deployment. To reconstruct this evolutionary route, we surveyed the expression of signaling molecules in the notochord of the tunicate Ciona, an experimentally amenable and informative chordate. We found that several genes encoding for candidate components of diverse signaling pathways are expressed during notochord development, and in some instances, display distinctive regionalized and/or lineage-specific patterns. We identified and deconstructed notochord enhancers associated with TGF-β and Ctgf, two evolutionarily conserved signaling genes that are expressed dishomogeneously in the Ciona notochord, and shed light on the cis-regulatory origins of their peculiar expression patterns.

Single-cell Transcriptomic Studies Unveil Potential Nodes of the Notochord Gene Regulatory Network

Transcription factors (TFs) are DNA-binding proteins able to modulate the timing, location, and levels of gene expression by binding to regulatory DNA regions. Therefore, the repertoire of TFs present in the genome of a multicellular organism and the expression of variable constellations of TFs in different cellular cohorts determine the distinctive characteristics of developing tissues and organs. The information on tissue-specific assortments of TFs, their cross-regulatory interactions, and the genes/regulatory regions targeted by each TF is summarized in gene regulatory networks (GRNs), which provide genetic blueprints for the specification, development, and differentiation of multicellular structures. In this study, we review recent transcriptomic studies focused on the complement of TFs expressed in the notochord, a distinctive feature of all chordates. We analyzed notochord-specific datasets available from organisms representative of the three chordate subphyla, and highlighted lineage-specific variations in the suite of TFs expressed in their notochord. We framed the resulting findings within a provisional evolutionary scenario, which allows the formulation of hypotheses on the genetic/genomic changes that sculpted the structure and function of the notochord on an evolutionary scale.

Neural crest lineage in the protovertebrate model Ciona

Neural crest cells are multipotent progenitors that produce defining features of vertebrates such as the 'new head'1. Here we use the tunicate, Ciona, to explore the evolutionary origins of neural crest since this invertebrate chordate is among the closest living relatives of vertebrates2-4. Previous studies identified two potential neural crest cell types in Ciona, sensory pigment cells and bipolar tail neurons5,6. Recent findings suggest that  bipolar tail neurons are homologous to cranial sensory ganglia rather than derivatives of neural crest7,8. Here we show that the pigment cell lineage also produces neural progenitor cells that form regions of the juvenile nervous system following metamorphosis. Neural progenitors are also a major derivative of neural crest in vertebrates, suggesting that the last common ancestor of tunicates and vertebrates contained a multipotent progenitor population at the neural plate border. It would therefore appear that a key property of neural crest, multipotentiality, preceded the emergence of vertebrates.

Early transcriptional similarities between two distinct neural lineages during ascidian embryogenesis

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

Specification and survival of post-metamorphic branchiomeric neurons in a non-vertebrate chordate

Tunicates are the sister group to the vertebrates, yet most species have a life cycle split between swimming larva and sedentary adult phases. During metamorphosis, larval neurons are replaced by adult-specific ones. The regulatory mechanisms underlying this replacement remain largely unknown. Using tissue-specific CRISPR/Cas9-mediated mutagenesis in the tunicate Ciona, we show that orthologs of conserved hindbrain and branchiomeric neuron regulatory factors Pax2/5/8 and Phox2 are required to specify the 'neck', a cellular compartment set aside in the larva to give rise to cranial motor neuron-like neurons post-metamorphosis. Using bulk and single-cell RNA-sequencing analyses, we characterize the transcriptome of the neck downstream of Pax2/5/8. We present evidence that neck-derived adult ciliomotor neurons begin to differentiate in the larva and persist through metamorphosis, contrary to the assumption that the adult nervous system is formed after settlement and the death of larval neurons during metamorphosis. Finally, we show that FGF signaling during the larval phase alters the patterning of the neck and its derivatives. Suppression of FGF converts neck cells into larval neurons that fail to survive metamorphosis, whereas prolonged FGF signaling promotes an adult neural stem cell-like fate.

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.

Development and circuitry of the tunicate larval Motor Ganglion, a putative hindbrain/spinal cord homolog

The Motor Ganglion (MG) is a small collection of neurons that control the swimming movements of the tunicate tadpole larva. Situated at the base of the tail, molecular and functional comparisons suggest that may be a homolog of the spinal cord and/or hindbrain ("rhombospinal" region) of vertebrates. Here we review the most current knowledge of the development, connectivity, functions, and unique identities of the neurons that comprise the MG, drawn mostly from studies in Ciona spp. The simple cell lineages, minimal cellular composition, and comprehensively mapped "connectome" of the Ciona MG all make this an excellent model for studying the development and physiology of motor control in aquatic larvae.

Cellular remodeling and JAK inhibition promote zygotic gene expression in the Ciona germline

Transcription control is a major determinant of cell fate decisions in somatic tissues. By contrast, early germline fate specification in numerous vertebrate and invertebrate species relies extensively on RNA-level regulation, exerted on asymmetrically inherited maternal supplies, with little-to-no zygotic transcription. However delayed, a maternal-to-zygotic transition is nevertheless poised to complete the deployment of pre-gametic programs in the germline. Here, we focus on early germline specification in the tunicate Ciona to study zygotic genome activation. We first demonstrate that a peculiar cellular remodeling event excludes localized postplasmic Pem-1 mRNA, which encodes the general inhibitor of transcription. Subsequently, zygotic transcription begins in Pem-1-negative primordial germ cells (PGCs), as revealed by histochemical detection of elongating RNA Polymerase II, and nascent Mef2 transcripts. In addition, we uncover a provisional antagonism between JAK and MEK/BMPRI/GSK3 signaling, which controls the onset of zygotic gene expression, following cellular remodeling of PGCs. We propose a 2-step model for the onset of zygotic transcription in the Ciona germline and discuss the significance of germ plasm dislocation and remodeling in the context of developmental fate specification.

Quantitative proteome dynamics across embryogenesis in a model chordate

The evolution of gene expression programs underlying the development of vertebrates remains poorly characterized. Here, we present a comprehensive proteome atlas of the model chordate Ciona, covering eight developmental stages and ∼7,000 translated genes, accompanied by a multi-omics analysis of co-evolution with the vertebrate Xenopus. Quantitative proteome comparisons argue against the widely held hourglass model, based solely on transcriptomic profiles, whereby peak conservation is observed during mid-developmental stages. Our analysis reveals maximal divergence at these stages, particularly gastrulation and neurulation. Together, our work provides a valuable resource for evaluating conservation and divergence of multi-omics profiles underlying the diversification of vertebrates.

Cis-regulatory interfaces reveal the molecular mechanisms underlying the notochord gene regulatory network of Ciona

Tissue-specific gene expression is fundamental in development and evolution, and is mediated by transcription factors and by the cis-regulatory regions (enhancers) that they control. Transcription factors and their respective tissue-specific enhancers are essential components of gene regulatory networks responsible for the development of tissues and organs. Although numerous transcription factors have been characterized from different organisms, the knowledge of the enhancers responsible for their tissue-specific expression remains fragmentary. Here we use Ciona to study the enhancers associated with ten transcription factors expressed in the notochord, an evolutionary hallmark of the chordate phylum. Our results illustrate how two evolutionarily conserved transcription factors, Brachyury and Foxa2, coordinate the deployment of other notochord transcription factors. The results of these detailed cis-regulatory analyses delineate a high-resolution view of the essential notochord gene regulatory network of Ciona, and provide a reference for studies of transcription factors, enhancers, and their roles in development, disease, and evolution.

Molecular control of cellulosic fin morphogenesis in ascidians

BACKGROUND

The tunicates form a group of filter-feeding marine animals closely related to vertebrates. They share with them a number of features such as a notochord and a dorsal neural tube in the tadpole larvae of ascidians, one of the three groups that make tunicates. However, a number of typical chordate characters have been lost in different branches of tunicates, a diverse and fast-evolving phylum. Consequently, the tunic, a sort of exoskeleton made of extracellular material including cellulose secreted by the epidermis, is the unifying character defining the tunicate phylum. In the larva of ascidians, the tunic differentiates in the tail into a median fin (with dorsal and ventral extended blades) and a caudal fin.

RESULTS

Here we have performed experiments in the ascidian Phallusia mammillata to address the molecular control of tunic 3D morphogenesis. We have demonstrated that the tail epidermis medio-lateral patterning essential for peripheral nervous system specification also controls tunic elongation into fins. More specifically, when tail epidermis midline identity was abolished by BMP signaling inhibition, or CRISPR/Cas9 inactivation of the transcription factor coding genes Msx or Klf1/2/4/17, median fin did not form. We postulated that this genetic program should regulate effectors of tunic secretion. We thus analyzed the expression and regulation in different ascidian species of two genes acquired by horizontal gene transfer (HGT) from bacteria, CesA coding for a cellulose synthase and Gh6 coding for a cellulase. We have uncovered an unexpected dynamic history of these genes in tunicates and high levels of variability in gene expression and regulation among ascidians. Although, in Phallusia, Gh6 has a regionalized expression in the epidermis compatible with an involvement in fin elongation, our functional studies indicate a minor function during caudal fin formation only.

CONCLUSIONS

Our study constitutes an important step in the study of the integration of HGT-acquired genes into developmental networks and a cellulose-based morphogenesis of extracellular material in animals.

Ovarian Follicle Development in Ascidians

Ovarian follicle development is an essential process for continuation of sexually reproductive animals, and is controlled by a wide variety of regulatory factors such as neuropeptides and peptide hormones in the endocrine, neuroendocrine, and nervous systems. Moreover, while some molecular mechanisms underlying follicle development are conserved, others vary among species. Consequently, follicle development processes are closely related to the evolution and diversity of species. Ciona intestinalis type A (Ciona rubusta) is a cosmopolitan species of ascidians, which are the closest relative of vertebrates. However, unlike vertebrates, ascidians are not endowed with the hypothalamus-pituitary-gonadal axis involving pituitary gonadotropins and sexual steroids. Combined with the phylogenetic position of ascidians as the closest relative of vertebrates, such morphological and endocrine features suggest that ascidians possess both common and species-specific regulatory mechanisms in follicle development. To date, several neuropeptides have been shown to participate in the growth of vitellogenic follicles, oocyte maturation of postvitellogenic follicles, and ovulation of fully mature follicles in a developmental stage-specific fashion. Furthermore, recent studies have shed light on the evolutionary processes of follicle development throughout chordates. In this review, we provide an overview of the neuropeptidergic molecular mechanism in the premature follicle growth, oocyte maturation, and ovulation in Ciona, and comparative views of the follicle development processes of mammals and teleosts.

Cereblon influences the timing of muscle differentiation in Ciona tadpoles

Thalidomide has a dark history as a teratogen, but in recent years, its derivates have been shown to function as potent chemotherapeutic agents. These drugs bind cereblon (CRBN), the substrate receptor of an E3 ubiquitin ligase complex, and modify its degradation targets. Despite these insights, remarkably little is known about the normal function of cereblon in development. Here, we employ Ciona, a simple invertebrate chordate, to identify endogenous Crbn targets. In Ciona, Crbn is specifically expressed in developing muscles during tail elongation before they acquire contractile activity. Crbn expression is activated by Mrf, the ortholog of MYOD1, a transcription factor important for muscle differentiation. CRISPR/Cas9-mediated mutations of Crbn lead to precocious onset of muscle contractions. By contrast, overexpression of Crbn delays contractions and is associated with decreased expression of contractile protein genes such as troponin. This reduction is possibly due to reduced Mrf protein levels without altering Mrf mRNA levels. Our findings suggest that Mrf and Crbn form a negative feedback loop to control the precision of muscle differentiation during tail elongation.

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.

Specification and survival of post-metamorphic branchiomeric neurons in the hindbrain of a non-vertebrate chordate

Tunicates are the sister group to the vertebrates, yet most species have a life cycle split between swimming larva and sedentary adult phases. During metamorphosis, larval neurons are largely replaced by adult-specific ones. Yet the regulatory mechanisms underlying this neural replacement remain largely unknown. Using tissue-specific CRISPR/Cas9-mediated mutagenesis in the tunicate Ciona, we show that orthologs of conserved hindbrain and branchiomeric neuron regulatory factors Pax2/5/8 and Phox2 are required to specify the "Neck", a compartment of cells set aside in the larva to give rise to cranial motor neuron-like neurons in the adult. Using bulk and single-cell RNAseq analyses, we also characterize the transcriptome of the Neck downstream of Pax2/5/8. Surprisingly, we find that Neck-derived adult ciliomotor neurons begin to differentiate in the larva, contrary to the long-held assumption that the adult nervous system is formed only after settlement and the death of larval neurons during metamorphosis. Finally, we show that manipulating FGF signaling during the larval phase alters the patterning of the Neck and its derivatives. Suppression of FGF converts Neck cells into larval neurons that fail to survive metamorphosis, while prolonged FGF signaling promotes an adult neural stem cell-like fate instead.

Mathematical modeling for developmental processes

We review several mathematical models and concepts in developmental biology that have been established over the last decade. (1) Feedback vertex set: Ascidian embryos contain cells of seven types, and cell fate is controlled by ~100 interacting genes. The "feedback vertex set" of the directed graph of the gene regulatory network consists of a small number of genes. By experimentally manipulating them, we can differentiate cells into any cell type. (2) Tissue deformation: Describing morphological changes in tissues and relating them to gene expression and other cellular processes is key in understanding morphogenesis. Expansion and anisotropy of the tissue are described by a "deformation tensor" at each location. A study on chick limb bud formation revealed that both the volume growth rate and anisotropy in deformation differed significantly between locations and stages. (3) Mechanobiology: Forces operating on each cell may alter cell shape and gene expression, which may subsequently exert forces on their surroundings. Measurements of force, tissue shape, and gene expression help us understand autonomous tissue deformation. (4) Adaptive design of development: An optimal growth schedule in fluctuating environments explains the growth response to starvation in Drosophila larvae. Adaptive placement of morphogen sources makes development robust to noises.

Controlling complex dynamical systems based on the structure of the networks

Progress of molecular biology resulted in the accumulation of information on biomolecular interactions, which are complex enough to be termed as networks. Dynamical behavior generated by complex network systems is considered to be the origin of the biological functions. One of the largest missions in modern life science is to obtain logical understanding for the dynamics of complex systems based on experimentally identified networks. However, a network does not provide sufficient information to specify dynamics explicitly, i.e. it lacks information of mathematical formulae of functions or parameter values. One has to develop mathematical models under assumptions of functions and parameter values to know the detail of dynamics of network systems. In this review, on the other hand, we introduce our own mathematical theory to understand the behavior of biological systems from the information of regulatory networks alone. Using the theory, important aspects of dynamical properties can be extracted from networks. Namely, key factors for observing/controlling the whole dynamical system are determined from network structure alone. We also show an application of the theory to a real biological system, a gene regulatory network for cell-fate specification in ascidian. We demonstrate that the system was completely controllable by experimental manipulations of the key factors identified by the theory from the information of network alone. This review article is an extended version of the Japanese article, Controlling Cell-Fate Specification System Based on a Mathematical Theory of Network Dynamics, published in SEIBUTSU BUTSURI Vol. 60, p. 349-351 (2020).

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.

Zic-r.b controls cell numbers in Ciona embryos by activating CDKN1B

The control of cell numbers and the establishment of cell types are two processes that are essential in early embryonic development. We have a reasonable understanding of how these processes occur individually, but we have considerably less sophisticated understanding of how these processes are linked. Tunicates have fixed cell lineages with predictable cell cycles, making them well suited to investigate these processes. In the ascidian Ciona, we show that the transcription factor Zic-r.b, known to be involved in establishing several cell types in early development also activates the expression of the cell cycle inhibitor CDKN1B. Zic-r.b is a major missing component of the cell division clock establishing specific cell numbers. We also show that a larvacean homolog of Zic-r.b is expressed one cell cycle earlier than its Ciona counterpart. The early expression in larvaceans may explain why they have half as many notochord cells as ascidians and may illustrate a general mechanism to evolve changes in morphology.

Diverse logics and grammar encode notochord enhancers

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

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.

Admp regulates tail bending by controlling ventral epidermal cell polarity via phosphorylated myosin localization in Ciona

Ventral tail bending, which is transient but pronounced, is found in many chordate embryos and constitutes an interesting model of how tissue interactions control embryo shape. Here, we identify one key upstream regulator of ventral tail bending in embryos of the ascidian Ciona. We show that during the early tailbud stages, ventral epidermal cells exhibit a boat-shaped morphology (boat cell) with a narrow apical surface where phosphorylated myosin light chain (pMLC) accumulates. We further show that interfering with the function of the BMP ligand Admp led to pMLC localizing to the basal instead of the apical side of ventral epidermal cells and a reduced number of boat cells. Finally, we show that cutting ventral epidermal midline cells at their apex using an ultraviolet laser relaxed ventral tail bending. Based on these results, we propose a previously unreported function for Admp in localizing pMLC to the apical side of ventral epidermal cells, which causes the tail to bend ventrally by resisting antero-posterior notochord extension at the ventral side of the tail.

DMRT Transcription Factors in the Control of Nervous System Sexual Differentiation

Sexual phenotypic differences in the nervous system are one of the most prevalent features across the animal kingdom. The molecular mechanisms responsible for sexual dimorphism throughout metazoan nervous systems are extremely diverse, ranging from intrinsic cell autonomous mechanisms to gonad-dependent endocrine control of sexual traits, or even extrinsic environmental cues. In recent years, the DMRT ancient family of transcription factors has emerged as being central in the development of sex-specific differentiation in all animals in which they have been studied. In this review, we provide an overview of the function of Dmrt genes in nervous system sexual regulation from an evolutionary perspective.

On the evolutionary origins and regionalization of the neural crest

The neural crest is a vertebrate-specific embryonic stem cell population that gives rise to a vast array of cell types throughout the animal body plan. These cells are first born at the edges of the central nervous system, from which they migrate extensively and differentiate into multiple cellular derivatives. Given the unique set of structures these cells comprise, the origin of the neural crest is thought to have important implications for the evolution and diversification of the vertebrate clade. In jawed vertebrates, neural crest cells exist as distinct subpopulations along the anterior-posterior axis. These subpopulations differ in terms of their respective differentiation potential and cellular derivatives. Thus, the modern neural crest is characterized as multipotent, migratory, and regionally segregated throughout the embryo. Here, we retrace the evolutionary origins of the neural crest, from the appearance of conserved regulatory circuitry in basal chordates to the emergence of neural crest subpopulations in higher vertebrates. Finally, we discuss a stepwise trajectory by which these cells may have arisen and diversified throughout vertebrate evolution.

Highly distinct genetic programs for peripheral nervous system formation in chordates

Background

Vertebrates develop their peripheral nervous system (PNS) from transient unique embryonic structures, the neural crest, and the ectodermal placodes that are located at the border of the forming central nervous system. By contrast, in the invertebrate chordates, amphioxus and ascidians, a large part of the PNS originates at the opposite of the embryo, in the ventral ectoderm. In both groups, a biphasic mechanism regulates ventral PNS formation: high BMP levels specify a neurogenic territory within which glutamatergic epidermal sensory neuron formation is controlled by the Notch pathway. Given these similarities and the phylogenetic relationships within chordates, it is likely that ventral PNS is an ancestral feature in chordates and that it has been lost in vertebrates.

Results

In order to get insights into the molecular control of ventral PNS formation and to test the hypothesis of their homology and potential contribution to the emergence of vertebrate PNS, we undertook a close comparison of ventral PNS formation in the ascidian Phallusia mammillata and the amphioxus Branchiostoma lanceolatum. Using timed RNA-seq series, we identified novel markers of the ventral PNS during different phases of its development in both species. By extensively determining the expression of paralogous and orthologous genes, we observed that only a minority of genes have a shared expression in the ventral PNS. However, a large fraction of ventral PNS orthologous genes are expressed in the dorsally forming PNS of vertebrates.

Conclusions

Our work has significantly increased the molecular characterization of ventral PNS formation in invertebrate chordates. The low observed conservation of gene expression in the ventral PNS suggests that the amphioxus and ascidian ventral PNS are either not homologous, or alternatively extensive drift has occurred in their regulatory mechanisms following a long period (600 My) of separate evolution and accelerated evolution in the ascidian lineage. The homology to genes expressed in the dorsally forming PNS of vertebrates suggests that ancestral sensory neurons gene networks have been redeployed in vertebrates.

ERK signaling dissolves ERF repression condensates in living embryos

Phase separation underlies the organization of the nucleus, including the biogenesis of nucleoli and the packaging of heterochromatin. Here we explore the regulation of transcription factor condensates involved in gene repression by ERK signaling in gastrulating embryos of a simple proto-vertebrate (Ciona). ERK signaling induces nuclear export of the transcriptional repressor Ets-2 repressive factor (ERF), which has been linked to various human developmental disorders. Using high-resolution imaging, we show that ERF is localized within discrete nuclear condensates that dissolve upon ERK activation. Interestingly, we observe dynamic pulses of assembly and dissociation during interphase, providing visualization of a nuclear phase separation process regulated by cell signaling. We discuss the implications of these observations for producing sharp on/off switches in gene activity and suppressing noise in cell-cell signaling events.

Quantitative Dissection of the Proximal Ciona brachyury Enhancer

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

Xbp1 and Brachyury establish an evolutionarily conserved subcircuit of the notochord gene regulatory network

Gene regulatory networks coordinate the formation of organs and structures that compose the evolving body plans of different organisms. We are using a simple chordate model, the Ciona embryo, to investigate the essential gene regulatory network that orchestrates morphogenesis of the notochord, a structure necessary for the proper development of all chordate embryos. Although numerous transcription factors expressed in the notochord have been identified in different chordates, several of them remain to be positioned within a regulatory framework. Here, we focus on Xbp1, a transcription factor expressed during notochord formation in Ciona and other chordates. Through the identification of Xbp1-downstream notochord genes in Ciona, we found evidence of the early co-option of genes involved in the unfolded protein response to the notochord developmental program. We report the regulatory interplay between Xbp1 and Brachyury, and by extending these results to Xenopus, we show that Brachyury and Xbp1 form a cross-regulatory subcircuit of the notochord gene regulatory network that has been consolidated during chordate evolution.

Ebf Activates Expression of a Cholinergic Locus in a Multipolar Motor Ganglion Interneuron Subtype in Ciona

Conserved transcription factors termed “terminal selectors” regulate neuronal sub-type specification and differentiation through combinatorial transcriptional regulation of terminal differentiation genes. The unique combinations of terminal differentiation gene products in turn contribute to the functional identities of each neuron. One well-characterized terminal selector is COE (Collier/Olf/Ebf), which has been shown to activate cholinergic gene batteries in C. elegans motor neurons. However, its functions in other metazoans, particularly chordates, is less clear. Here we show that the sole COE ortholog in the non-vertebrate chordate Ciona robusta, Ebf, controls the expression of the cholinergic locus VAChT/ChAT in a single dorsal interneuron of the larval Motor Ganglion, which is presumed to be homologous to the vertebrate spinal cord. We propose that, while the function of Ebf as a regulator of cholinergic neuron identity conserved across bilaterians, its exact role may have diverged in different cholinergic neuron subtypes (e.g., interneurons vs. motor neurons) in chordate-specific motor circuits.

Omics Studies for the Identification of Ascidian Peptides, Cognate Receptors, and Their Relevant Roles in Ovarian Follicular Development

Omics studies contribute to the elucidation of genomes and profiles of gene expression. In the ascidian Ciona intestinalis Type A (Ciona robusta), mass spectrometry (MS)-based peptidomic studies have detected numerous Ciona-specific (nonhomologous) neuropeptides as well as Ciona homologs of typical vertebrate neuropeptides and hypothalamic peptide hormones. Candidates for cognate G protein-coupled receptors (GPCRs) for these peptides have been found in the Ciona transcriptome by two ways. First, Ciona homologous GPCRs of vertebrate counterparts have been detected by sequence homology searches of cognate transcriptomes. Second, the transcriptome-derived GPCR candidates have been used for machine learning-based systematic prediction of interactions not only between Ciona homologous peptides and GPCRs but also between novel Ciona peptides and GPCRs. These data have ultimately led to experimental evidence for various Ciona peptide-GPCR interactions. Comparative transcriptomics between the wildtype and Ciona vasopressin (CiVP) gene-edited Ciona provide clues to the biological functions of CiVP in ovarian follicular development and whole body growth. Furthermore, the transcriptomes of follicles treated with peptides, such as Ciona tachykinin and cionin (a Ciona cholecystokinin homolog), have revealed key regulatory genes for Ciona follicle growth, maturation, and ovulation, eventually leading to the verification of essential and novel molecular mechanisms underlying these biological events. These findings indicate that omics studies, combined with artificial intelligence and single-cell technologies, pave the way for investigating in greater details the nervous, neuroendocrine, and endocrine systems of ascidians and the molecular and functional evolution and diversity of peptidergic regulatory networks throughout chordates.

Developmental biology of the larvacean Oikopleura dioica: Genome resources, functional screening, and imaging

The larvacean Oikopleura dioica is a cosmopolitan planktonic chordate and is closely related to vertebrates. It is characterized by a tadpole-shaped morphology with notochord flanked by muscle in the tail and brain on the dorsal side, a short life cycle of five days, a compact genome of approximately 56 Mb, a simple and transparent body with a small number of cells (~4000 in functional juveniles), invariant embryonic cell lineages, and fast development that ensures complete morphogenesis and organ formation 10 h after fertilization. With these features, this marine chordate is a promising and advantageous animal model in which genetic manipulation is feasible. In this review, we introduce relevant resources and modern techniques that have been developed: (1) Genome and transcriptomes. Oikopleura dioica has the smallest genome among non-parasitic metazoans. Its genome databases have been generated using three geographically distant O. dioica populations, and several intra-species sequence differences are becoming evident; (2) Functional genetic knockdown techniques. Comprehensive screening of genes is feasible using ovarian microinjection and double-strand DNA-induced gene knockdown; and (3) Live imaging of embryos and larvae. Application of these techniques has uncovered novel aspects of development, including meiotic cell arrest, left-right patterning, epidermal cell patterning, and mouth formation involving the connection of ectoderm and endoderm sheets. Oikopleura dioca has become very useful for developmental and evolutionary studies in chordates.

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.

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

BACKGROUND

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

RESULTS AND CONCLUSIONS

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

A Preliminary Single-Cell RNA-Seq Analysis of Embryonic Cells That Express Brachyury in the Amphioxus, Branchiostoma japonicum

Among chordate taxa, the cephalochordates diverged earlier than urochordates and vertebrates; thus, they retain unique, primitive developmental features. In particular, the amphioxus notochord has muscle-like properties, a feature not seen in urochordates or vertebrates. Amphioxus contains two Brachyury genes, Bra1 and Bra2. Bra2 is reportedly expressed in the blastopore, notochord, somites, and tail bud, in contrast to a low level of Bra1 expression only in notochord. To distinguish the expression profiles of the two Brachyury genes at the single-cell level, we carried out single-cell RNA-seq (scRNA-seq) analysis using the amphioxus, Branchiostoma japonicum. This scRNA-seq analysis classified B. japonicum embryonic cells into 15 clusters at developmental stages from midgastrula to early swimming larva. Brachyury was expressed in cells of clusters 4, 5, 8, and 9. We first confirmed that cluster 8 comprises cells that form somites since this cluster specifically expresses four myogenic factor genes. Cluster 9 contains a larger number of cells with high levels of Bra2 expression and a smaller number of cells with Bra1 expression. Simultaneous expression in cluster 9 of tool-kit genes, including FoxA, Goosecoid, and hedgehog, showed that this cluster comprises cells that form the notochord. Expression of Bra2, but not Bra1, in cells of clusters 4 and 5 at the gastrula stage together with expression of Wnt1 and Caudal indicates that clusters 4 and 5 comprise cells of the blastopore, which contiguously form the tail bud. In addition, Hox1, Hox3, and Hox4 were highly expressed in Bra2-expressing clusters 4, 5, 8, and 9 in a temporally coordinated manner, suggesting roles of anterior Hox genes in specification of mesodermal organs, including somites, notochord, and tail bud. This scRNA-seq analysis therefore highlights differences between the two Brachyury genes in relation to embryonic regions in which they are expressed and their levels of expression. Bra2 is the ancestral Brachyury in amphioxus, since expression in the blastopore is shared with other deuterostomes. On the other hand, Bra1 is a duplicate copy and likely evolved a supplementary function in notochord and somite formation in the Branchiostoma lineage.

Evolution of Developmental Programs for the Midline Structures in Chordates: Insights From Gene Regulation in the Floor Plate and Hypochord Homologues of Ciona Embryos

In vertebrate embryos, dorsal midline tissues, including the notochord, the prechordal plate, and the floor plate, play important roles in patterning of the central nervous system, somites, and endodermal tissues by producing extracellular signaling molecules, such as Sonic hedgehog (Shh). In Ciona, hedgehog.b, one of the two hedgehog genes, is expressed in the floor plate of the embryonic neural tube, while none of the hedgehog genes are expressed in the notochord. We have identified a cis-regulatory region of hedgehog.b that was sufficient to drive a reporter gene expression in the floor plate. The hedgehog.b cis-regulatory region also drove ectopic expression of the reporter gene in the endodermal strand, suggesting that the floor plate and the endodermal strand share a part of their gene regulatory programs. The endodermal strand occupies the same topographic position of the embryo as does the vertebrate hypochord, which consists of a row of single cells lined up immediately ventral to the notochord. The hypochord shares expression of several genes with the floor plate, including Shh and FoxA, and play a role in dorsal aorta development. Whole-embryo single-cell transcriptome analysis identified a number of genes specifically expressed in both the floor plate and the endodermal strand in Ciona tailbud embryos. A Ciona FoxA ortholog FoxA.a is shown to be a candidate transcriptional activator for the midline gene battery. The present findings suggest an ancient evolutionary origin of a common developmental program for the midline structures in Olfactores.

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.

Heterologous functional expression of ascidian Nav1 channels and close relationship with the evolutionary ancestor of vertebrate Nav channels

Voltage-gated sodium channels (Nav1s) are responsible for the initiation and propagation of action potentials in neurons, muscle, and endocrine cells. Many clinically used drugs such as local anesthetics and antiarrhythmics inhibit Nav1s, and a variety of inherited human disorders are caused by mutations in Nav1 genes. Nav1s consist of the main α subunit and several auxiliary β subunits. Detailed information on the structure–function relationships of Nav1 subunits has been obtained through heterologous expression experiments and analyses of protein structures. The basic properties of Nav1s, including their gating and ion permeation, were classically described in the squid giant axon and other invertebrates. However, heterologous functional expression of Nav1s from marine invertebrates has been unsuccessful. Ascidians belong to the Urochordata, a sister group of vertebrates, and the larval central nervous system of ascidians shows a similar plan to that of vertebrates. Here, we report the biophysical properties of ascidian Ciona Nav1 (CiNav1a) heterologously expressed in Xenopus oocytes. CiNav1a exhibited tetrodotoxin-insensitive sodium currents with rapid gating kinetics of activation and inactivation. Furthermore, consistent with the fact that the Ciona genome lacks orthologous genes to vertebrate β subunits, the human β1 subunit did not influence the gating properties when coexpressed with CiNav1a. Interestingly, CiNav1a contains an ankyrin-binding motif in the II–III linker, which can be targeted to the axon initial segment of mammalian cortical neurons. Our findings provide a platform to gain insight into the evolutionary and biophysical properties of Nav1s, which are important for the development of targeted therapeutics.

Neuromesodermal Lineage Contribution to CNS Development in Invertebrate and Vertebrate Chordates

Ascidians are invertebrate chordates and the closest living relative to vertebrates. In ascidian embryos a large part of the central nervous system arises from cells associated with mesoderm rather than ectoderm lineages. This seems at odds with the traditional view of vertebrate nervous system development which was thought to be induced from ectoderm cells, initially with anterior character and later transformed by posteriorizing signals, to generate the entire anterior-posterior axis of the central nervous system. Recent advances in vertebrate developmental biology, however, show that much of the posterior central nervous system, or spinal cord, in fact arises from cells that share a common origin with mesoderm. This indicates a conserved role for bi-potential neuromesoderm precursors in chordate CNS formation. However, the boundary between neural tissue arising from these distinct neural lineages does not appear to be fixed, which leads to the notion that anterior-posterior patterning and neural fate formation can evolve independently.

Expression and Functional Analyses of Ectodermal Transcription Factors FoxJ-r, SoxF, and SP8/9 in Early Embryos of the Ascidian Halocynthia roretzi

The spatiotemporal expression of zygotic genes is regulated by transcription factors, which mediate cell fate decision and morphogenesis. Investigation of the expression patterns and their transcriptional regulatory relationships is crucial to understand embryonic development. Staged RNA-seq of the ascidian Halocynthia roretzi has previously shown that nine genes encoding transcription factors are transiently expressed at the blastula stage, which is the stage at which cell fates are specified and differentiation starts. Six of these transcription factors have already been found to play important roles during early development. However, the functions of the other transcription factors (FoxJ-r, SoxF, and SP8/9) remain unknown. The study of the spatial and temporal expression patterns showed that all three genes were expressed in the animal hemisphere as early as the 16-cell stage. This is likely due to transcription factor genes that are expressed in the vegetal hemisphere, which have been extensively and comprehensively analyzed in previous studies of ascidians. Functional analyses using FoxJ-r morphants showed that they resulted in the disruption of laterality and the absence of epidermal mono-cilia, suggesting FoxJ-r functions in cilia formation and, consequently, in the generation of left-right asymmetry, as observed in vertebrates. SoxF knockdown resulted in incomplete epiboly by the ectoderm during gastrulation, while SP8/9 knockdown showed no phenotype until the tailbud stage in the present study, although it was expressed during blastula stages. Our results indicate that transcription factor genes expressed at the cleavage stages play roles in diverse functions, and are not limited to cell fate specification.

Properties of repression condensates in living Ciona embryos

Recent studies suggest that transcriptional activators and components of the pre-initiation complex (PIC) form higher order associations-clusters or condensates-at active loci. Considerably less is known about the distribution of repressor proteins responsible for gene silencing. Here, we develop an expression assay in living Ciona embryos that captures the liquid behavior of individual nucleoli undergoing dynamic fusion events. The assay is used to visualize puncta of Hes repressors, along with the Groucho/TLE corepressor. We observe that Hes.a/Gro puncta have the properties of viscous liquid droplets that undergo limited fusion events due to association with DNA. Hes.a mutants that are unable to bind DNA display hallmarks of liquid-liquid phase separation, including dynamic fusions of individual condensates to produce large droplets. We propose that the DNA template serves as a scaffold for the formation of Hes condensates, but limits the spread of transcriptional repressors to unwanted regions of the genome.

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.

Evolution of new cell types at the lateral neural border

During the course of evolution, animals have become increasingly complex by the addition of novel cell types and regulatory mechanisms. A prime example is represented by the lateral neural border, known as the neural plate border in vertebrates, a region of the developing ectoderm where presumptive neural and non-neural tissue meet. This region has been intensively studied as the source of two important embryonic cell types unique to vertebrates-the neural crest and the ectodermal placodes-which contribute to diverse differentiated cell types including the peripheral nervous system, pigment cells, bone, and cartilage. How did these multipotent progenitors originate in animal evolution? What triggered the elaboration of the border during the course of chordate evolution? How is the lateral neural border patterned in various bilaterians and what is its fate? Here, we review and compare the development and fate of the lateral neural border in vertebrates and invertebrates and we speculate about its evolutionary origin. Taken together, the data suggest that the lateral neural border existed in bilaterian ancestors prior to the origin of vertebrates and became a developmental source of exquisite evolutionary change that frequently enabled the acquisition of new cell types.

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

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

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

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

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.