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.

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.

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.

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.

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.

Whole-organ cell shape analysis reveals the developmental basis of ascidian notochord taper.

Here we use in toto imaging together with computational segmentation and analysis methods to quantify the shape of every cell at multiple stages in the development of a simple organ: the notochord of the ascidian Ciona savignyi. We find that cell shape in the intercalated notochord depends strongly on anterior-posterior (AP) position, with cells in the middle of the notochord consistently wider than cells at the anterior or posterior. This morphological feature of having a tapered notochord is present in many chordates. We find that ascidian notochord taper involves three main mechanisms: Planar Cell Polarity (PCP) pathway-independent sibling cell volume asymmetries that precede notochord cell intercalation; the developmental timing of intercalation, which proceeds from the anterior and posterior towards the middle; and the differential rates of notochord cell narrowing after intercalation. A quantitative model shows how the morphology of an entire developing organ can be controlled by this small set of cellular mechanisms.

Ciona genetics.

Ascidians, such as Ciona, are invertebrate chordates with simple embryonic body plans and small, relatively non-redundant genomes. Ciona genetics is in its infancy compared to many other model systems, but it provides a powerful method for studying this important vertebrate outgroup. Here we give basic methods for genetic analysis of Ciona, including protocols for controlled crosses both by natural spawning and by the surgical isolation of gametes; the identification and propagation of mutant lines; and strategies for positional cloning.

The ascidian mouth opening is derived from the anterior neuropore: reassessing the mouth/neural tube relationship in chordate evolution.

The relative positions of the brain and mouth are of central importance for models of chordate evolution. The dorsal hollow neural tube and the mouth have often been thought of as developmentally distinct structures that may have followed independent evolutionary paths. In most chordates however, including vertebrates and ascidians, the mouth primordia have been shown to fate to the anterior neural boundary. In ascidians such as Ciona there is a particularly intimate relationship between brain and mouth development, with a thin canal connecting the neural tube lumen to the mouth primordium at larval stages. This so-called neurohypophyseal canal was previously thought to be a secondary connection that formed relatively late, after the independent formation of the mouth primordium and the neural tube. Here we show that the Ciona neurohypophyseal canal is present from the end of neurulation and represents the anteriormost neural tube, and that the future mouth opening is actually derived from the anterior neuropore. The mouth thus forms at the anterior midline transition between neural tube and surface ectoderm. In the vertebrate Xenopus, we find that although the mouth primordium is not topologically continuous with the neural tube lumen, it nonetheless forms at this same transition point. This close association between the mouth primordium and the anterior neural tube in both ascidians and amphibians suggests that the evolution of these two structures may be more closely linked than previously appreciated.

Chongmague reveals an essential role for laminin-mediated boundary formation in chordate convergence and extension movements.

Although cell intercalation driven by non-canonical Wnt/planar cell polarity (PCP) pathway-dependent mediolateral cell polarity is important for notochord morphogenesis, it is likely that multiple mechanisms shape the notochord as it converges and extends. Here we show that the recessive short-tailed Ciona savignyi mutation chongmague (chm) has a novel defect in the formation of a morphological boundary around the developing notochord. chm notochord cells initiate intercalation normally, but then fail to maintain their polarized cell morphology and migrate inappropriately to become dispersed in the larval tail. This is unlike aimless (aim), a mutation in the PCP pathway component Prickle, which has a severe defect in early mediolateral intercalation but forms a robust notochord boundary. Positional cloning identifies chm as a mutation in the C. savignyi ortholog of the vertebrate alpha 3/4/5 family of laminins. Cs-lamalpha3/4/5 is highly expressed in the developing notochord, and Cs-lamalpha3/4/5 protein is specifically localized to the outer border of the notochord. Notochord convergence and extension, reduced but not absent in both chm and aim, are essentially abolished in the aim/aim; chm/chm double mutant, indicating that laminin-mediated boundary formation and PCP-dependent mediolateral intercalation are each able to drive a remarkable degree of tail morphogenesis in the absence of the other. These mechanisms therefore initially act in parallel, but we also find that PCP signaling has an important later role in maintaining the perinotochordal/intranotochordal polarity of Cs-lamalpha3/4/5 localization.

A second canon. Functions and mechanisms of beta-catenin-independent Wnt signaling.

More is becoming known about so-called noncanonical Wnt pathways that signal independently of beta-catenin. Here we review recent developments in both the functions and mechanisms of noncanonical Wnt signaling. We also discuss some unresolved and vexing questions. How many noncanonical Wnt pathways are there? How extensive are the parallels between Drosophila planar polarization and vertebrate convergence and extension? Last, we will outline some challenges and difficulties we foresee for this exciting but still very young field.

Zebrafish prickle, a modulator of noncanonical Wnt/Fz signaling, regulates gastrulation movements.

In addition to the canonical Wnt/beta-catenin signaling pathway, at least two noncanonical Wnt/Fz pathways have been described: the planar cell polarity (PCP) pathway in Drosophila [1] and the Wnt/calcium pathway in vertebrate embryos [2]. Recent work suggests that a vertebrate pathway homologous to the PCP pathway acts to regulate the convergent extension movements of gastrulation [3-7]. To further test this hypothesis, we have identified two zebrafish homologs of the Drosophila PCP gene prickle (pk) [8], both of which show discrete and dynamic expression patterns during gastrulation. Both gain and loss of pk1 function cause defects in convergent extension. Pk1 localizes to both the cytoplasm and the cell membrane, and its normal localization is partially dependent on its C-terminal prenylation motif. At the cell membrane, Pk1 is frequently localized asymmetrically around the cell and can colocalize with the signaling molecule Dishevelled (Dsh). In overexpression assays, Pk1 is able to activate AP-1-mediated transcription and inhibit activation of Wnt/beta-catenin signaling. Like noncanonical Wnts [9-10], overexpression of Pk1 increases the frequency of calcium transients in zebrafish blastulae. Our results support the idea that a vertebrate PCP pathway regulates gastrulation movements and suggest that there is overlap between the PCP and Wnt/calcium pathways.