Co-expression of Foxa.a, Foxd and Fgf9/16/20 defines a transient mesendoderm regulatory state in ascidian embryos

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.

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.

A comprehensive study of arthropod and onychophoran Fox gene expression patterns

Fox genes represent an evolutionary old class of transcription factor encoding genes that evolved in the last common ancestor of fungi and animals. They represent key-components of multiple gene regulatory networks (GRNs) that are essential for embryonic development. Most of our knowledge about the function of Fox genes comes from vertebrate research, and for arthropods the only comprehensive gene expression analysis is that of the fly Drosophila melanogaster. For other arthropods, only selected Fox genes have been investigated. In this study, we provide the first comprehensive gene expression analysis of arthropod Fox genes including representative species of all main groups of arthropods, Pancrustacea, Myriapoda and Chelicerata. We also provide the first comprehensive analysis of Fox gene expression in an onychophoran species. Our data show that many of the Fox genes likely retained their function during panarthropod evolution highlighting their importance in development. Comparison with published data from other groups of animals shows that this high degree of evolutionary conservation often dates back beyond the last common ancestor of Panarthropoda.

New developments in the biology of fibroblast growth factors

The fibroblast growth factor (FGF) family is composed of 18 secreted signaling proteins consisting of canonical FGFs and endocrine FGFs that activate four receptor tyrosine kinases (FGFRs 1-4) and four intracellular proteins (intracellular FGFs or iFGFs) that primarily function to regulate the activity of voltage-gated sodium channels and other molecules. The canonical FGFs, endocrine FGFs, and iFGFs have been reviewed extensively by us and others. In this review, we briefly summarize past reviews and then focus on new developments in the FGF field since our last review in 2015. Some of the highlights in the past 6 years include the use of optogenetic tools, viral vectors, and inducible transgenes to experimentally modulate FGF signaling, the clinical use of small molecule FGFR inhibitors, an expanded understanding of endocrine FGF signaling, functions for FGF signaling in stem cell pluripotency and differentiation, roles for FGF signaling in tissue homeostasis and regeneration, a continuing elaboration of mechanisms of FGF signaling in development, and an expanding appreciation of roles for FGF signaling in neuropsychiatric diseases. This article is categorized under: Cardiovascular Diseases > Molecular and Cellular Physiology Neurological Diseases > Molecular and Cellular Physiology Congenital Diseases > Stem Cells and Development Cancer > Stem Cells and Development.

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

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

Reprint of: Conditional specification of endomesoderm

Early in animal development many cells are conditionally specified based on observations that those cells can be directed toward alternate fates. The endomesoderm is so named because early specification produces cells that often have been observed to simultaneously express both early endoderm and mesoderm transcription factors. Experiments with these cells demonstrate that their progeny can directed entirely toward endoderm or mesoderm, whereas normally they establish both germ layers. This review examines the mechanisms that initiate the conditional endomesoderm state, its metastability, and the mechanisms that resolve that state into definitive endoderm and mesoderm.

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.

Conditional specification of endomesoderm

Early in animal development many cells are conditionally specified based on observations that those cells can be directed toward alternate fates. The endomesoderm is so named because early specification produces cells that often have been observed to simultaneously express both early endoderm and mesoderm transcription factors. Experiments with these cells demonstrate that their progeny can directed entirely toward endoderm or mesoderm, whereas normally they establish both germ layers. This review examines the mechanisms that initiate the conditional endomesoderm state, its metastability, and the mechanisms that resolve that state into definitive endoderm and mesoderm.

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.

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.

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

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

Using linkage logic theory to control dynamics of a gene regulatory network of a chordate embryo

Linkage logic theory provides a mathematical criterion to control network dynamics by manipulating activities of a subset of network nodes, which are collectively called a feedback vertex set (FVS). Because many biological functions emerge from dynamics of biological networks, this theory provides a promising tool for controlling biological functions. By manipulating the activity of FVS molecules identified in a gene regulatory network (GRN) for fate specification of seven tissues in ascidian embryos, we previously succeeded in reproducing six of the seven cell types. Simultaneously, we discovered that the experimentally reconstituted GRN lacked information sufficient to reproduce muscle cells. Here, we utilized linkage logic theory as a tool to find missing edges in the GRN. Then, we identified a FVS from an updated version of the GRN and confirmed that manipulating the activity of this FVS was sufficient to induce all seven cell types, even in a multi-cellular environment. Thus, linkage logic theory provides tools to find missing edges in experimentally reconstituted networks, to determine whether reconstituted networks contain sufficient information to fulfil expected functions, and to reprogram cell fate.

A serpin is required for ectomesoderm, a hallmark of spiralian development

Among animals, diploblasts contain two germ layers, endoderm and ectoderm, while triploblasts have a distinct third germ layer called the mesoderm. Spiralians are a group of triploblast animals that have highly conserved development: they share the distinctive spiralian cleavage pattern as well as a unique source of mesoderm, the ectomesoderm. This population of mesoderm is distinct from endomesoderm and is considered a hallmark of spiralian development, but the regulatory network that drives its development is unknown. Here we identified ectomesoderm-specific genes in the mollusc Tritia (aka Ilyanassa) obsoleta through differential gene expression analyses comparing control and ectomesoderm-ablated embryos, followed by in situ hybridization of identified transcripts. We identified a Tritia serpin gene (ToSerpin1) that appears to be specifically expressed in the ectomesoderm of the posterior and head. Ablation of the 3a and 3b cells, which make most of the ectomesoderm, abolishes ToSerpin1 expression, consistent with its expression in these cells. Morpholino knockdown of ToSerpin1 causes ectomesoderm defects, most prominently in the muscle system of the larval head. This is the first gene identified that is specifically implicated in spiralian ectomesoderm development.

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.

Tunicate gastrulation

Tunicates are a diverse group of invertebrate marine chordates that includes the larvaceans, thaliaceans, and ascidians. Because of their unique evolutionary position as the sister group of the vertebrates, tunicates are invaluable as a comparative model and hold the promise of revealing both conserved and derived features of chordate gastrulation. Descriptive studies in a broad range of tunicates have revealed several important unifying traits that make them unique among the chordates, including invariant cell lineages through gastrula stages and an overall morphological simplicity. Gastrulation has only been studied in detail in ascidians such as Ciona and Phallusia, where it involves a simple cup-shaped gastrula driven primarily by endoderm invagination. This appears to differ significantly from vertebrate models, such as Xenopus, in which mesoderm convergent extension and epidermal epiboly are major contributors to involution. These differences may reflect the cellular simplicity of the ascidian embryo.

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

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

The Ciona myogenic regulatory factor functions as a typical MRF but possesses a novel N-terminus that is essential for activity

Electroporation-based assays were used to test whether the myogenic regulatory factor (MRF) of Ciona intestinalis (CiMRF) interferes with endogenous developmental programs, and to evaluate the importance of its unusual N-terminus for muscle development. We found that CiMRF suppresses both notochord and endoderm development when it is expressed in these tissues by a mechanism that may involve activation of muscle-specific microRNAs. Because these results add to a large body of evidence demonstrating the exceptionally high degree of functional conservation among MRFs, we were surprised to discover that non-ascidian MRFs were not myogenic in Ciona unless they formed part of a chimeric protein containing the CiMRF N-terminus. Equally surprising, we found that despite their widely differing primary sequences, the N-termini of MRFs of other ascidian species could form chimeric MRFs that were also myogenic in Ciona. This domain did not rescue the activity of a Brachyury protein whose transcriptional activation domain had been deleted, and so does not appear to constitute such a domain. Our results indicate that ascidians have previously unrecognized and potentially novel requirements for MRF-directed myogenesis. Moreover, they provide the first example of a domain that is essential to the core function of an important family of gene regulatory proteins, one that, to date, has been found in only a single branch of the family.

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

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

Transcriptional regulation of the Ciona Gsx gene in the neural plate

The ascidian neural plate consists of a defined number of identifiable cells organized in a grid of rows and columns, representing a useful model to investigate the molecular mechanisms controlling neural patterning in chordates. Distinct anterior brain lineages are specified via unique combinatorial inputs of signalling pathways with Nodal and Delta-Notch signals patterning along the medial-lateral axis and FGF/MEK/ERK signals patterning along the anterior-posterior axis of the neural plate. The Ciona Gsx gene is specifically expressed in the a9.33 cells in the row III/column 2 position of anterior brain lineages, characterised by a combinatorial input of Nodal-OFF, Notch-ON and FGF-ON. Here, we identify the minimal cis-regulatory element (CRE) of 376 bp, which can recapitulate the early activation of Gsx. We show that this minimal CRE responds in the same way as the endogenous Gsx gene to manipulation of FGF- and Notch-signalling pathways and to overexpression of Snail, a mediator of Nodal signals, and Six3/6, which is required to demarcate the anterior boundary of Gsx expression at the late neurula stage. We reveal that sequences proximal to the transcription start site include a temporal regulatory element required for the precise transcriptional onset of gene expression. We conclude that sufficient spatial and temporal information for Gsx expression is integrated in 376 bp of non-coding cis-regulatory sequences.

Expression and prognosis analyses of forkhead box A (FOXA) family in human lung cancer

Despite advances in early diagnosis and treatment, cancer still remains the major reason of mortality worldwide. The forkhead box A (FOXA) family is reported to participate in diverse human diseases. However, little is known about their expression and prognostic values in human lung cancer. Herein, we conducted a detailed cancer vs. normal analysis. The mRNA expression levels of FOXA family in numerous kind of cancers, including lung cancer, were analyzed using the Oncomine and GEPIA database. We observed that the mRNA expression levels of FOXA1, and FOXA3 were all increased while FOXA2 were decreased in most cancers compared with normal tissues, especially in lung cancer. Moreover, the expression levels of FOXA1, and FOXA3 are also highly expressed, while FOXA2 were decreased in almost all cancer cell lines, particularly in lung cancer cell lines, analyzing by Cancer Cell Line Encyclopedia (CCLE) and EMBL-EBI databases. Furthermore, the LinkedOmics database was used to evaluate the prognostic values, indicating that higher expression of FOXA1, FOXA3 indicated a poor overall survival (OS), while increased FOXA2 revealed a better OS in lung cancer. To conclusion, FOXA family showed significant expression differences between cancer and normal tissues, especially lung cancer, and FOXA1, FOXA3 could be promising prognostic biomarkers for lung cancer.

Ascidian Zic Genes

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

Controlling Cell Fate Specification System by Key Genes Determined from Network Structure

Network structures describing regulation between biomolecules have been determined in many biological systems. Dynamics of molecular activities based on such networks are considered to be the origin of many biological functions. Recently, it has been proved mathematically that key nodes for controlling dynamics in networks are identified from network structure alone. Here, we applied this theory to a gene regulatory network for the cell fate specification of seven tissues in the ascidian embryo and found that this network, which consisted of 92 factors, had five key molecules. By controlling the activities of these key molecules, the specific gene expression of six of seven tissues observed in the embryo was successfully reproduced. Since this method is applicable to all nonlinear dynamic systems, we propose this method as a tool for controlling gene regulatory networks and reprogramming cell fates.

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

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

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

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

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

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

The forkhead transcription factor UNC-130/FOXD integrates both BMP and Notch signaling to regulate dorsoventral patterning of the C. elegans postembryonic mesoderm

The proper development of a multicellular organism requires precise spatial and temporal coordination of cell intrinsic and cell extrinsic regulatory mechanisms. Both Notch signaling and bone morphogenetic protein (BMP) signaling function to regulate the proper development of the C. elegans postembryonic mesoderm. We have identified the C. elegans FOXD transcription factor UNC-130 as a major target functioning downstream of both BMP signaling and Notch signaling to regulate dorsoventral patterning of the postembryonic mesoderm. We showed that unc-130 expression in the postembryonic M lineage is asymmetric: its absence of expression in the dorsal side of the M lineage requires the antagonism of BMP signaling by the zinc finger transcription factor SMA-9, while its expression in the ventral side of the M lineage is activated by LIN-12/Notch signaling. We further showed that the regulation of UNC-130 expression by BMP signaling and Notch signaling is specific to the M lineage, as the ventral expression of UNC-130 in the embryonically-derived bodywall muscles was not affected in either BMP pathway or Notch pathway mutants. Finally, we showed that the function of UNC-130 in the M lineage is independent of UNC-129, a gene previously shown to function downstream of and be repressed by UNC-130 for axon guidance. Our studies uncovered a new function of UNC-130/FOXD in the C. elegans postembryonic mesoderm, and identify UNC-130 as a critical factor that integrates two independent spatial cues for the proper patterning and fate specification of the C. elegans postembryonic mesoderm.

FOXA transcriptional factor modulates insect susceptibility to Bacillus thuringiensis Cry1Ac toxin by regulating the expression of toxin-receptor ABCC2 and ABCC3 genes

Cry toxins produced by Bacillus thuringiensis (Bt) are insecticidal proteins widely used in insect control. Recently, it was shown that ATP-binding cassette transporter proteins (ABC) such as ABCC2, ABCC3, ABCG1 and ABCA2 are implicated in the insecticidal action of Cry toxins as putative receptors. However, the transcriptional regulators involved in the expression of ABC transporter genes remain unknown. Sequence analysis of promoter regions of ABCC2 gene from Helicoverpa armigera and ABCC3 gene from Spodoptera litura Sl-HP cultured cells, revealed the potential participation of Forkhead box protein A (FOXA), a transcription factor that regulates the expression of genes through remodeling chromatin. To determine if FOXA was involved in regulating expression of ABCC2 and ABCC3 genes, the expression of FOXA, ABCC2 and ABCC3 was compared in Sl-HP cells that are sensitive to Cry1Ac toxin with those in S. frugiperda Sf9 cells that are not sensitive to the toxin. Expression levels of those genes were significantly higher in Sl-HP than in Sf9 cells. Transient expression of FOXA in Sf9 cells activated ABCC2 and ABCC3 transcription, which directly correlated with enhanced Cry1Ac-susceptibility in these cells. Silencing of FOXA gene expression by RNAi in H. armigera larvae resulted in a decreased expression of ABCC2 and ABCC3 without affecting expression of other Cry toxin receptor genes such as alkaline phosphatase, aminopeptidase or cadherin. Silencing of FOXA gene expression also resulted in a Cry1Ac-tolerant phenotype since lower mortality and higher pupation rate were observed in diet containing Cry1Ac protoxin in comparison with the control group. These results demonstrate that FOXA up-regulates expression of the Cry1Ac-toxin receptor ABCC2 and ABCC3 genes, and that lower FOXA expression correlates with tolerance to Cry toxin in cell lines and in lepidopteran larvae.

Cas9-mediated excision of Nematostella brachyury disrupts endoderm development, pharynx formation and oral-aboral patterning

The mesoderm is a key novelty in animal evolution, although we understand little of how the mesoderm arose. brachyury, the founding member of the T-box gene family, is a key gene in chordate mesoderm development. However, the brachyury gene was present in the common ancestor of fungi and animals long before mesoderm appeared. To explore ancestral roles of brachyury prior to the evolution of definitive mesoderm, we excised the gene using CRISPR/Cas9 in the diploblastic cnidarian Nematostella vectensis Nvbrachyury is normally expressed in precursors of the pharynx, which separates endoderm from ectoderm. In knockout embryos, the pharynx does not form, embryos fail to elongate, and endoderm organization, ectodermal cell polarity and patterning along the oral-aboral axis are disrupted. Expression of many genes both inside and outside the Nvbrachyury expression domain is affected, including downregulation of Wnt genes at the oral pole. Our results point to an ancient role for brachyury in morphogenesis, cell polarity and the patterning of both ectodermal and endodermal derivatives along the primary body axis.

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

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

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

Patterning of brain precursors in ascidian embryos

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

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

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