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.

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.

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.

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.

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.

ANISEED 2019: 4D exploration of genetic data for an extended range of tunicates.

ANISEED (https://www.aniseed.cnrs.fr) is the main model organism database for the worldwide community of scientists working on tunicates, the vertebrate sister-group. Information provided for each species includes functionally-annotated gene and transcript models with orthology relationships within tunicates, and with echinoderms, cephalochordates and vertebrates. Beyond genes the system describes other genetic elements, including repeated elements and cis-regulatory modules. Gene expression profiles for several thousand genes are formalized in both wild-type and experimentally-manipulated conditions, using formal anatomical ontologies. These data can be explored through three complementary types of browsers, each offering a different view-point. A developmental browser summarizes the information in a gene- or territory-centric manner. Advanced genomic browsers integrate the genetic features surrounding genes or gene sets within a species. A Genomicus synteny browser explores the conservation of local gene order across deuterostome. This new release covers an extended taxonomic range of 14 species, including for the first time a non-ascidian species, the appendicularian Oikopleura dioica. Functional annotations, provided for each species, were enhanced through a combination of manual curation of gene models and the development of an improved orthology detection pipeline. Finally, gene expression profiles and anatomical territories can be explored in 4D online through the newly developed Morphonet morphogenetic browser.

Iterative and Complex Asymmetric Divisions Control Cell Volume Differences in Ciona Notochord Tapering.

The notochord of the invertebrate chordate Ciona forms a tapered rod at tailbud stages consisting of only 40 cylindrical cells in a single-file column. This tapered shape involves differences in notochord cell volume along the anterior-posterior axis. Here, we quantify sibling cell volume asymmetry throughout the developing notochord and find that there are distinctive patterns of unequal cleavage in all 4 bilateral pairs of A-line primary notochord founder cells and also in the B-line-derived secondary notochord founder cells. A quantitative model confirms that the observed patterns of unequal cleavage are sufficient to explain all the anterior-posterior variation in notochord cell volume. Many examples are known of cells that divide asymmetrically to give daughter cells of different size and fate. Here, by contrast, a series of subtle but iterative and finely patterned asymmetric divisions controls the shape of an entire organ. Quantitative 3D analysis of cell shape and spindle positioning allows us to infer multiple cellular mechanisms driving these unequal cleavages, including polarized displacements of the mitotic spindle, contributions from the shape of the mother cell, and late changes occurring between anaphase and abscission that potentially involve differential cortical contractility. We infer differential use of these mechanisms between different notochord blastomeres and also between different rounds of cell division. These results demonstrate a new role for asymmetric division in directly shaping a developing organ and point toward complex underlying mechanisms.

Multiple inputs into a posterior-specific regulatory network in the Ciona notochord.

The gene regulatory networks underlying Ciona notochord fate specification and differentiation have been extensively investigated, but the regulatory basis for regionalized expression within the notochord is not understood. Here we identify three notochord-expressed genes, C11.331, C12.115 and C8.891, with strongly enriched expression in the secondary notochord cells at the posterior tip of the tail. C11.331 and C12.115 share a distinctive expression pattern that is highly enriched in the secondary notochord lineage but also graded within that lineage with the strongest expression at the posterior tip. Both genes show similar responses to pharmacological perturbations of Wnt and FGF signaling, consistent with an important role for Wnt and FGF ligands expressed at the tail tip. Reporter analysis indicates that the C11.331 cis-regulatory regions are extensively distributed, with multiple non-overlapping regions conferring posterior notochord-enriched expression. Fine-scale analysis of a minimal cis-regulatory module identifies discrete positive and negative elements including a strong silencer. Truncation of the silencer region leads to increased expression in the primary notochord, indicating that C11.331 expression is influenced by putative regulators of primary versus secondary notochord fate. The minimal CRM contains predicted ETS, GATA, LMX and Myb sites, all of which lead to reduced expression in secondary notochord when mutated. These results show that the posterior-enriched notochord expression of C11.331 depends on multiple inputs, including Wnt and FGF signals from the tip of the tail, multiple notochord-specific regulators, and yet-to-be identified regulators of regional identity within the notochord.

The Ciona Notochord Gene Regulatory Network.

Complex gene regulatory networks are at the heart of cell fate specification and differentiation. The simple chordate Ciona has remarkable advantages for the dissection of these regulatory networks, including a stereotypically chordate but extremely small and simple embryo, a streamlined and compact genome, and highly efficient transgenesis by electroporation. Here we use the Ciona notochord as an example of how these characteristics can be exploited to understand both the early network controlling cell fate as well as the tissue-specific network controlling notochord differentiation and morphogenesis.

A temperature-adjusted developmental timer for precise embryonic staging.

Developmental biology research depends on careful staging of developing embryos, but the rate of development is extremely sensitive to the temperature at which embryos are raised. It is not always practical to grow embryos at a precisely controlled temperature, so here we describe a simple, inexpensive device based on an Arduino-compatible microprocessor and temperature sensor that provides a metric of developmental time that compensates for changes in temperature. The underlying assumption is that the rate of development will be linear with respect to temperature over an organism's thermal tolerance range. The device measures the ambient temperature and integrates effective degree-minutes over time. For convenience, this is displayed to the user as a temperature-adjusted standard developmental time. In initial testing we have found the device to be extremely helpful for fixing Ciona embryos during precise developmental windows.

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

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

Functional and evolutionary insights from the Ciona notochord transcriptome.

The notochord of the ascidian Ciona consists of only 40 cells, and is a longstanding model for studying organogenesis in a small, simple embryo. Here, we perform RNAseq on flow-sorted notochord cells from multiple stages to define a comprehensive Ciona notochord transcriptome. We identify 1364 genes with enriched expression and extensively validate the results by in situ hybridization. These genes are highly enriched for Gene Ontology terms related to the extracellular matrix, cell adhesion and cytoskeleton. Orthologs of 112 of the Ciona notochord genes have known notochord expression in vertebrates, more than twice as many as predicted by chance alone. This set of putative effector genes with notochord expression conserved from tunicates to vertebrates will be invaluable for testing hypotheses about notochord evolution. The full set of Ciona notochord genes provides a foundation for systems-level studies of notochord gene regulation and morphogenesis. We find only modest overlap between this set of notochord-enriched transcripts and the genes upregulated by ectopic expression of the key notochord transcription factor Brachyury, indicating that Brachyury is not a notochord master regulator gene as strictly defined.

Dynamics of cell polarity in tissue morphogenesis: a comparative view from Drosophila and Ciona.

Tissues in developing embryos exhibit complex and dynamic rearrangements that shape forming organs, limbs, and body axes. Directed migration, mediolateral intercalation, lumen formation, and other rearrangements influence the topology and topography of developing tissues. These collective cell behaviors are distinct phenomena but all involve the fine-grained control of cell polarity. Here we review recent findings in the dynamics of polarized cell behavior in both the Drosophila ovarian border cells and the Ciona notochord. These studies reveal the remarkable reorganization of cell polarity during organ formation and underscore conserved mechanisms of developmental cell polarity including the Par/atypical protein kinase C (aPKC) and planar cell polarity pathways. These two very different model systems demonstrate important commonalities but also key differences in how cell polarity is controlled in tissue morphogenesis. Together, these systems raise important, broader questions on how the developmental control of cell polarity contributes to morphogenesis of diverse tissues across the metazoa.

Covert Prepatterning of a Cell Division Wave.

A directional wave of mitosis, progressing posterior to anterior across the epidermis, is important for neural tube closure in the invertebrate chordate Ciona intestinalis. In this issue of Developmental Cell, Ogura and Sasakura (2016) show that the patterning of this wave unexpectedly has complex origins in the previous cell cycle.

Reciprocal and dynamic polarization of planar cell polarity core components and myosin.

The Ciona notochord displays planar cell polarity (PCP), with anterior localization of Prickle (Pk) and Strabismus (Stbm). We report that a myosin is polarized anteriorly in these cells and strongly colocalizes with Stbm. Disruption of the actin/myosin machinery with cytochalasin or blebbistatin disrupts polarization of Pk and Stbm, but not of myosin complexes, suggesting a PCP-independent aspect of myosin localization. Wash out of cytochalasin restored Pk polarization, but not if done in the presence of blebbistatin, suggesting an active role for myosin in core PCP protein localization. On the other hand, in the pk mutant line, aimless, myosin polarization is disrupted in approximately one third of the cells, indicating a reciprocal action of core PCP signaling on myosin localization. Our results indicate a complex relationship between the actomyosin cytoskeleton and core PCP components in which myosin is not simply a downstream target of PCP signaling, but also required for PCP protein localization.

Stochasticity and stereotypy in the Ciona notochord.

Fate mapping with single cell resolution has typically been confined to embryos with completely stereotyped development. The lineages giving rise to the 40 cells of the Ciona notochord are invariant, but the intercalation of those cells into a single-file column is not. Here we use genetic labeling methods to fate map the Ciona notochord with both high resolution and large sample sizes. We find that the ordering of notochord cells into a single column is not random, but instead shows a distinctive signature characteristic of mediolaterally-biased intercalation. We find that patterns of cell intercalation in the notochord are somewhat stochastic but far more stereotyped than previously believed. Cell behaviors vary by lineage, with the secondary notochord lineage being much more constrained than the primary lineage. Within the primary lineage, patterns of intercalation reflect the geometry of the intercalating tissue. We identify the latest point at which notochord morphogenesis is largely stereotyped, which is shortly before the onset of mediolateral intercalation and immediately after the final cell divisions in the primary lineage. These divisions are consistently oriented along the AP axis. Our results indicate that the interplay between stereotyped and stochastic cell behaviors in morphogenesis can only be assessed by fate mapping experiments that have both cellular resolution and large sample sizes.

Quantitative and in toto imaging in ascidians: working toward an image-centric systems biology of chordate morphogenesis.

Developmental biology relies heavily on microscopy to image the finely controlled cell behaviors that drive embryonic development. Most embryos are large enough that a field of view with the resolution and magnification needed to resolve single cells will not span more than a small region of the embryo. Ascidian embryos, however, are sufficiently small that they can be imaged in toto with fine subcellular detail using conventional microscopes and objectives. Unlike other model organisms with particularly small embryos, ascidians have a chordate embryonic body plan that includes a notochord, hollow dorsal neural tube, heart primordium and numerous other anatomical details conserved with the vertebrates. Here we compare the size and anatomy of ascidian embryos with those of more traditional model organisms, and relate these features to the capabilities of both conventional and exotic imaging methods. We review the emergence of Ciona and related ascidian species as model organisms for a new era of image-based developmental systems biology. We conclude by discussing some important challenges in ascidian imaging and image analysis that remain to be solved.

Exploiting the extraordinary genetic polymorphism of ciona for developmental genetics with whole genome sequencing.

Studies in tunicates such as Ciona have revealed new insights into the evolutionary origins of chordate development. Ciona populations are characterized by high levels of natural genetic variation, between 1 and 5%. This variation has provided abundant material for forward genetic studies. In the current study, we make use of deep sequencing and homozygosity mapping to map spontaneous mutations in outbred populations. With this method we have mapped two spontaneous developmental mutants. In Ciona intestinalis we mapped a short-tail mutation with strong phenotypic similarity to a previously identified mutant in the related species Ciona savignyi. Our bioinformatic approach mapped the mutation to a narrow interval containing a single mutated gene, α-laminin3,4,5, which is the gene previously implicated in C. savignyi. In addition, we mapped a novel genetic mutation disrupting neural tube closure in C. savignyi to a T-type Ca(2+) channel gene. The high efficiency and unprecedented mapping resolution of our study is a powerful advantage for developmental genetics in Ciona, and may find application in other outbred species.

Anterior-posterior regionalized gene expression in the Ciona notochord.

BACKGROUND: In the simple ascidian chordate Ciona, the signaling pathways and gene regulatory networks giving rise to initial notochord induction are largely understood and the mechanisms of notochord morphogenesis are being systematically elucidated. The notochord has generally been thought of as a non-compartmentalized or regionalized organ that is not finely patterned at the level of gene expression. Quantitative imaging methods have recently shown, however, that notochord cell size, shape, and behavior vary consistently along the anterior-posterior (AP) axis. RESULTS: Here we screen candidate genes by whole mount in situ hybridization for potential AP asymmetry. We identify 4 genes that show non-uniform expression in the notochord. Ezrin/radixin/moesin (ERM) is expressed more strongly in the secondary notochord lineage than the primary. CTGF is expressed stochastically in a subset of notochord cells. A novel calmodulin-like gene (BCamL) is expressed more strongly at both the anterior and posterior tips of the notochord. A TGF-ß ortholog is expressed in a gradient from posterior to anterior. The asymmetries in ERM, BCamL, and TGF-ß expression are evident even before the notochord cells have intercalated into a single-file column. CONCLUSIONS: We conclude that the Ciona notochord is not a homogeneous tissue but instead shows distinct patterns of regionalized gene expression.

3D-printed microwell arrays for Ciona microinjection and timelapse imaging.

Ascidians such as Ciona are close chordate relatives of the vertebrates with small, simple embryonic body plans and small, simple genomes. The tractable size of the embryo offers considerable advantages for in toto imaging and quantitative analysis of morphogenesis. For functional studies, Ciona eggs are considerably more challenging to microinject than the much larger eggs of other model organisms such as zebrafish and Xenopus. One of the key difficulties is in restraining the eggs so that the microinjection needle can be easily introduced and withdrawn. Here we develop and test a device to cast wells in agarose that are each sized to hold a single egg. This injection mold is fabricated by micro-resolution stereolithography with a grid of egg-sized posts that cast corresponding wells in agarose. This 3D printing technology allows the rapid and inexpensive testing of iteratively refined prototypes. In addition to their utility in microinjection, these grids of embryo-sized wells are also valuable for timelapse imaging of multiple embryos.