MorphoSeq: Full Single-Cell Transcriptome Dynamics Up to Gastrulation in a Chordate.

Single-cell RNA sequencing (scRNA-seq) provides a leap forward in resolving cellular diversity and developmental trajectories but fails to comprehensively delineate the spatial organization and precise cellular makeup of individual embryos. Here, we reconstruct from scRNA-seq and light sheet imaging data a canonical digital embryo that captures the genome-wide gene expression trajectory of every single cell at every cell division in the 18 lineages up to gastrulation in the ascidian Phallusia mammillata. By using high-coverage scRNA-seq, we devise a computational framework that stratifies single cells of individual embryos into cell types without prior knowledge. Unbiased transcriptome data analysis mapped each cell's physical position and lineage history, yielding the complete history of gene expression at the genome-wide level for every single cell in a developing embryo. A comparison of individual embryos reveals both extensive reproducibility between symmetric embryo sides and a large inter-embryonic variability due to small differences in embryogenesis timing.

Extreme genome scrambling in marine planktonic Oikopleura dioica cryptic species.

Genome structural variations within species are rare. How selective constraints preserve gene order and chromosome structure is a central question in evolutionary biology that remains unsolved. Our sequencing of several genomes of the appendicularian tunicate Oikopleura dioica around the globe reveals extreme genome scrambling caused by thousands of chromosomal rearrangements, although showing no obvious morphological differences between these animals. The breakpoint accumulation rate is an order of magnitude higher than in ascidian tunicates, nematodes, Drosophila, or mammals. Chromosome arms and sex-specific regions appear to be the primary unit of macrosynteny conservation. At the microsyntenic level, scrambling did not preserve operon structures, suggesting an absence of selective pressure to maintain them. The uncoupling of the genome scrambling with morphological conservation in O. dioica suggests the presence of previously unnoticed cryptic species and provides a new biological system that challenges our previous vision of speciation in which similar animals always share similar genome structures.

A change in cis-regulatory logic underlying obligate versus facultative muscle multinucleation in chordates.

Vertebrates and tunicates are sister groups that share a common fusogenic factor, Myomaker (Mymk), that drives myoblast fusion and muscle multinucleation. Yet they are divergent in when and where they express Mymk. In vertebrates, all developing skeletal muscles express Mymk and are obligately multinucleated. In tunicates, Mymk is expressed only in post-metamorphic multinucleated muscles, but is absent from mononucleated larval muscles. In this study, we demonstrate that cis-regulatory sequence differences in the promoter region of Mymk underlie the different spatiotemporal patterns of its transcriptional activation in tunicates and vertebrates. Although in vertebrates myogenic regulatory factors (MRFs) such as MyoD1 alone are required and sufficient for Mymk transcription in all skeletal muscles, we show that transcription of Mymk in post-metamorphic muscles of the tunicate Ciona requires the combinatorial activity of MRF, MyoD and Early B-cell Factor (Ebf). This macroevolutionary difference appears to be encoded in cis, likely due to the presence of a putative Ebf-binding site adjacent to predicted MRF binding sites in the Ciona Mymk promoter. We further discuss how Mymk and myoblast fusion might have been regulated in the last common ancestor of tunicates and vertebrates, for which we propose two models.

Specification of distinct cell types in a sensory-adhesive organ important for metamorphosis in tunicate larvae.

The papillae of tunicate larvae contribute sensory, adhesive, and metamorphosis-regulating functions that are crucial for the biphasic lifestyle of these marine, non-vertebrate chordates. We have identified additional molecular markers for at least 5 distinct cell types in the papillae of the model tunicate Ciona, allowing us to further study the development of these organs. Using tissue-specific CRISPR/Cas9-mediated mutagenesis and other molecular perturbations, we reveal the roles of key transcription factors and signaling pathways that are important for patterning the papilla territory into a highly organized array of different cell types and shapes. We further test the contributions of different transcription factors and cell types to the production of the adhesive glue that allows for larval attachment during settlement, and to the processes of tail retraction and body rotation during metamorphosis. With this study, we continue working towards connecting gene regulation to cellular functions that control the developmental transition between the motile larva and sessile adult of Ciona.

Cardiopharyngeal deconstruction and ancestral tunicate sessility.

A central question in chordate evolution is the origin of sessility in adult ascidians, and whether the appendicularian complete free-living style represents a primitive or derived condition among tunicates1. According to the 'a new heart for a new head' hypothesis, the evolution of the cardiopharyngeal gene regulatory network appears as a pivotal aspect to understand the evolution of the lifestyles of chordates2-4. Here we show that appendicularians experienced massive ancestral losses of cardiopharyngeal genes and subfunctions, leading to the 'deconstruction' of two ancestral modules of the tunicate cardiopharyngeal gene regulatory network. In ascidians, these modules are related to early and late multipotency, which is involved in lineage cell-fate determination towards the first and second heart fields and siphon muscles. Our work shows that the deconstruction of the cardiopharyngeal gene regulatory network involved the regressive loss of the siphon muscle, supporting an evolutionary scenario in which ancestral tunicates had a sessile ascidian-like adult lifestyle. In agreement with this scenario, our findings also suggest that this deconstruction contributed to the acceleration of cardiogenesis and the redesign of the heart into an open-wide laminar structure in appendicularians as evolutionary adaptations during their transition to a complete pelagic free-living style upon the innovation of the food-filtering house5.

Regulation of anterior neurectoderm specification and differentiation by BMP signaling in ascidians.

The most anterior structure of the ascidian larva is made of three palps with sensory and adhesive functions essential for metamorphosis. They derive from the anterior neural border and their formation is regulated by FGF and Wnt. Given that they also share gene expression profiles with vertebrate anterior neural tissue and cranial placodes, their study should shed light on the emergence of the unique vertebrate telencephalon. We show that BMP signaling regulates two phases of palp formation in Ciona intestinalis. During gastrulation, the anterior neural border is specified in a domain of inactive BMP signaling, and activating BMP prevented its formation. During neurulation, BMP defines ventral palp identity and indirectly specifies the inter-papilla territory separating the ventral and dorsal palps. Finally, we show that BMP has similar functions in the ascidian Phallusia mammillata, for which we identified novel palp markers. Collectively, we provide a better molecular description of palp formation in ascidians that will be instrumental for comparative studies.

Genome-wide analysis of early vascular tunic repair and regeneration for Botrylloides digenesis reveals striking similarities to human wound healing.

The early stages of regeneration after injury are similar to those of wound healing. The ascidian Botrylloides diegensis can regenerate an entire adult from a small fragment of vascular tunic following the removal of all zooids in an injury-induced regeneration model. We investigated the molecular and cellular changes following injury to determine the differences between the healing process and the initiation of whole-body regeneration (WBR). We conducted transcriptome analysis at specific time points during regeneration and wound healing to identify differentially expressed genes (DEGs) and the unique biological processes associated with each state. Our findings revealed 296 DEGs at 10 h post-injury (hpi), with 71 highly expressed in healed tissue and 225 expressed during the WBR process. These DEGs were predicted to play roles in tissue reorganization, integrin signaling, extracellular matrix organization, and the innate immune system. Pathway analysis of the upregulated genes in the healed tunic indicated functional enrichment related to tissue repair, as has been observed in other species. Additionally, we examined the cell types in the tunic and ampullae in both tissue states using histology and in situ hybridization for six genes identified by transcriptome analysis. We observed strong mRNA expression in cells within the WBR tunic, and in small RNA-positive granules near the tunic edge. We hypothesized that many of these genes function in the compaction of the ampullae tunic, which is a pivotal process for WBR and dormancy in B. diegensis, and in an immune response. These findings establish surprising similarities between ascidian regeneration and human wound healing, emphasizing the potential for future investigations into human regenerative and repair mechanisms. This study provides valuable insights into the gene sets specifically activated during regeneration compared to wound healing, shedding light on the divergent activities of these processes.

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.

Pax3/7 gene function in Oikopleura dioica supports a neuroepithelial-like origin for its house-making Fol territory.

Larvacean tunicates feature a spectacular innovation not seen in other animals - the trunk oikoplastic epithelium (OE). This epithelium produces a house, a large and complex extracellular structure used for filtering and concentrating food particles. Previously we identified several homeobox transcription factor genes expressed during early OE patterning. Among these are two Pax3/7 copies that we named pax37A and pax37B. The vertebrate homologs, PAX3 and PAX7 are involved in developmental processes related to neural crest and muscles. In the ascidian tunicate Ciona intestinalis, Pax3/7 plays a role in the development of cells deriving from the neural plate border, including trunk epidermal sensory neurons and tail nerve cord neurons, as well as in the neural tube closure. Here we have investigated the roles of Oikopleura dioica pax37A and pax37B in the development of the OE, by using CRISPR-Cas9 mutant lines and analyzing scRNA-seq data from wild-type animals. We found that pax37B but not pax37A is essential for the differentiation of cell fields that produce the food concentrating filter of the house: the anterior Fol, giant Fol and Nasse cells. Trajectory analysis supported a neuroepithelial-like or a preplacodal ectoderm transcriptional signature in these cells. We propose that the highly specialized secretory epithelial cells of the Fol region either maintained or evolved neuroepithelial features. This is supported by a fragmented gene regulatory network involved in their development that also operates in ascidian epidermal neurons.

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.

IFI16 phase separation via multi-phosphorylation drives innate immune signaling.

The interferon inducible protein 16 (IFI16) is a prominent sensor of nuclear pathogenic DNA, initiating innate immune signaling and suppressing viral transcription. However, little is known about mechanisms that initiate IFI16 antiviral functions or its regulation within the host DNA-filled nucleus. Here, we provide in vitro and in vivo evidence to establish that IFI16 undergoes liquid-liquid phase separation (LLPS) nucleated by DNA. IFI16 binding to viral DNA initiates LLPS and induction of cytokines during herpes simplex virus type 1 (HSV-1) infection. Multiple phosphorylation sites within an intrinsically disordered region (IDR) function combinatorially to activate IFI16 LLPS, facilitating filamentation. Regulated by CDK2 and GSK3ß, IDR phosphorylation provides a toggle between active and inactive IFI16 and the decoupling of IFI16-mediated cytokine expression from repression of viral transcription. These findings show how IFI16 switch-like phase transitions are achieved with temporal resolution for immune signaling and, more broadly, the multi-layered regulation of nuclear DNA sensors.

Two distinct evolutionary conserved neural degeneration pathways characterized in a colonial chordate.

Colonial tunicates are marine organisms that possess multiple brains simultaneously during their colonial phase. While the cyclical processes of neurogenesis and neurodegeneration characterizing their life cycle have been documented previously, the cellular and molecular changes associated with such processes and their relationship with variation in brain morphology and individual (zooid) behavior throughout adult life remains unknown. Here, we introduce Botryllus schlosseri as an invertebrate model for neurogenesis, neural degeneration, and evolutionary neuroscience. Our analysis reveals that during the weekly colony budding (i.e., asexual reproduction), prior to programmed cell death and removal by phagocytes, decreases in the number of neurons in the adult brain are associated with reduced behavioral response and significant change in the expression of 73 mammalian homologous genes associated with neurodegenerative disease. Similarly, when comparing young colonies (1 to 2 y of age) to those reared in a laboratory for ∼20 y, we found that older colonies contained significantly fewer neurons and exhibited reduced behavioral response alongside changes in the expression of 148 such genes (35 of which were differentially expressed across both timescales). The existence of two distinct yet apparently related neurodegenerative pathways represents a novel platform to study the gene products governing the relationship between aging, neural regeneration and degeneration, and loss of nervous system function. Indeed, as a member of an evolutionary clade considered to be a sister group of vertebrates, this organism may be a fundamental resource in understanding how evolution has shaped these processes across phylogeny and obtaining mechanistic insight.

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.

Two distinct motifs for Zic-r.a drive specific gene expression in two cell lineages.

Zic-r.a, a maternal transcription factor, specifies posterior fate in ascidian embryos. However, its direct target, Tbx6-r.b, does not contain typical Zic-r.a-binding sites in its regulatory region. Using an in vitro selection assay, we found that Zic-r.a binds to sites dissimilar to the canonical motif, by which it activates Tbx6-r.b in a sub-lineage of muscle cells. These sites with non-canonical motifs have weak affinity for Zic-r.a; therefore, it activates Tbx6-r.b only in cells expressing Zic-r.a abundantly. Meanwhile, we found that Zic-r.a expressed zygotically in late embryos activates neural genes through canonical sites. Because different zinc-finger domains of Zic-r.a are important for driving reporters with canonical and non-canonical sites, it is likely that the non-canonical motif is not a divergent version of the canonical motif. In other words, our data indicate that the non-canonical motif represents a motif distinct from the canonical motif. Thus, Zic-r.a recognizes two distinct motifs to activate two sets of genes at two timepoints in development. This article has an associated 'The people behind the papers' interview.

The past, the recent, and the ongoing evolutionary processes of the worldwide invasive ascidian Styela plicata.

Invasive species are one of the main threats to global biodiversity and, within marine ecosystems, tunicates feature some prominent examples. Styela plicata is an ascidian species inhabiting harbours in all temperate oceans and seas, thus being considered a thriving invasive species. However, this species' adaptive mechanisms, introduction history, and population structure have never been completely elucidated. Here, by genotyping 87 S. plicata individuals from 18 localities worldwide with 2b-RADseq, we confirm the global presence of four chromosome inversions, demonstrate population structuring on this species, detect local adaptation signals, and infer historical demographic events. We show that North Carolina individuals constitute an unrelated population, Atlanto-Mediterranean and Pacific localities form their own genetic clusters with substructuring, being the most evident the split between northern and southern Atlantic localities. The locality of South Carolina presents an intermediate genetic position between North Carolina and the other two groups pointing to a hybrid origin with recurrent gene flow. We generate and test demographic models, providing evidence of two independent introduction events to the Atlantic and Pacific, and an admixture that originated the population of South Carolina. Finally, we identify candidate loci for adaptation, with functions involved with cell processes, metabolism, development, and ion transport, among others. Overall, this study highlights the complex historical processes of S. plicata, which have led this species to its current distribution, population structure, and local adaptation footprint in oceans worldwide.

Model of neural induction in the ascidian embryo.

How cell specification can be controlled in a reproducible manner is a fundamental question in developmental biology. In ascidians, a group of invertebrate chordates, geometry plays a key role in achieving this control. Here, we use mathematical modeling to demonstrate that geometry dictates the neural-epidermal cell fate choice in the 32-cell stage ascidian embryo by a two-step process involving first the modulation of ERK signaling and second, the expression of the neural marker gene, Otx. The model describes signal transduction by the ERK pathway that is stimulated by FGF and attenuated by ephrin, and ERK-mediated control of Otx gene expression, which involves both an activator and a repressor of ETS-family transcription factors. Considering the measured area of cell surface contacts with FGF- or ephrin-expressing cells as inputs, the solutions of the model reproduce the experimental observations about ERK activation and Otx expression in the different cells under normal and perturbed conditions. Sensitivity analyses and computations of Hill coefficients allow us to quantify the robustness of the specification mechanism controlled by cell surface area and to identify the respective role played by each signaling input. Simulations also predict in which conditions the dual control of gene expression by an activator and a repressor that are both under the control of ERK can induce a robust ON/OFF control of neural fate induction.

Complex mammalian-like haematopoietic system found in a colonial chordate.

Haematopoiesis is an essential process that evolved in multicellular animals. At the heart of this process are haematopoietic stem cells (HSCs), which are multipotent and self-renewing, and generate the entire repertoire of blood and immune cells throughout an animal's life1. Although there have been comprehensive studies on self-renewal, differentiation, physiological regulation and niche occupation in vertebrate HSCs, relatively little is known about the evolutionary origin and niches of these cells. Here we describe the haematopoietic system of Botryllus schlosseri, a colonial tunicate that has a vasculature and circulating blood cells, and interesting stem-cell biology and immunity characteristics2-8. Self-recognition between genetically compatible B. schlosseri colonies leads to the formation of natural parabionts with shared circulation, whereas incompatible colonies reject each other3,4,7. Using flow cytometry, whole-transcriptome sequencing of defined cell populations and diverse functional assays, we identify HSCs, progenitors, immune effector cells and an HSC niche, and demonstrate that self-recognition inhibits allospecific cytotoxic reactions. Our results show that HSC and myeloid lineage immune cells emerged in a common ancestor of tunicates and vertebrates, and also suggest that haematopoietic bone marrow and the B. schlosseri endostyle niche evolved from a common origin.

Genomics-informed models reveal extensive stretches of coastline under threat by an ecologically dominant invasive species.

Explaining why some species are widespread, while others are not, is fundamental to biogeography, ecology, and evolutionary biology. A unique way to study evolutionary and ecological mechanisms that either limit species' spread or facilitate range expansions is to conduct research on species that have restricted distributions. Nonindigenous species, particularly those that are highly invasive but have not yet spread beyond the introduced site, represent ideal systems to study range size changes. Here, we used species distribution modeling and genomic data to study the restricted range of a highly invasive Australian marine species, the ascidian Pyura praeputialis This species is an aggressive space occupier in its introduced range (Chile), where it has fundamentally altered the coastal community. We found high genomic diversity in Chile, indicating high adaptive potential. In addition, genomic data clearly showed that a single region from Australia was the only donor of genotypes to the introduced range. We identified over 3,500 km of suitable habitat adjacent to its current introduced range that has so far not been occupied, and importantly species distribution models were only accurate when genomic data were considered. Our results suggest that a slight change in currents, or a change in shipping routes, may lead to an expansion of the species' introduced range that will encompass a vast portion of the South American coast. Our study shows how the use of population genomics and species distribution modeling in combination can unravel mechanisms shaping range sizes and forecast future range shifts of invasive species.

Dynamics of Chromatin Opening across Larval Development in the Urochordate Ascidian Ciona savignyi.

Ascidian larvae undergo tail elongation and notochord lumenogenesis, making them an ideal model for investigating tissue morphogenesis in embryogenesis. The cellular and mechanical mechanisms of these processes have been studied; however, the underlying molecular regulatory mechanism remains to be elucidated. In this study, assays for transposase-accessible chromatin using sequencing (ATAC-seq) and RNA sequencing (RNA-seq) were applied to investigate potential regulators of the development of ascidian Ciona savignyi larvae. Our results revealed 351 and 138 differentially accessible region genes through comparisons of ATAC-seq data between stages 21 and 24 and between stages 24 and 25, respectively. A joint analysis of RNA-seq and ATAC-seq data revealed a correlation between chromatin accessibility and gene transcription. We further verified the tissue expression patterns of 12 different genes. Among them, Cs-matrix metalloproteinase 24 (MMP24) and Cs-krüppel-like factor 5 (KLF5) were highly expressed in notochord cells. Functional assay results demonstrated that both genes are necessary for notochord lumen formation and expansion. Finally, we performed motif enrichment analysis of the differentially accessible regions in different tailbud stages and summarized the potential roles of these motif-bearing transcription factors in larval development. Overall, our study found a correlation between gene expression and chromatin accessibility and provided a vital resource for understanding the mechanisms of the development of ascidian embryos.

Tunicate Eco-Evo-Devo laboratory in IMS-METU.

I completed my undergraduate education in Atatürk University, Education Faculty, Biology Department. Then pursued my graduate education at the Biology Department of Mersin University. Both my master's and PhD theses were on the biological and population genetics features of various fish species. My initial encounter with tunicates dates back to my Postdoc at Israel Oceanographic and Limnologic Research Institute (IOLR) in 2011, where I was working on a DNA barcoding project. During that time, the entire institute was actively engaged in research on tunicates, and discussions during lunchtime often revolved around this fascinating group of organisms. Prof. Rinkevich usually only spoke seriously about tunicate biology but 1 day he told me "You know Botryllus schlosseri is riding horse in Black Sea coasts of Turkiye." I was totally surprised and was trying to understand the meaning of this comment from a scientific perspective. He then showed me the picture of a B. schlosseri colony attached to a seahorse. Following several more Postdoc experiences, I began working as a Principal Investigator at Institute of Marine Sciences, Middle East Technical University (IMS-METU) in 2017. Since then, my team and I have been working on tunicate biodiversity, evolutionary biology, genomics, DNA barcoding, metabarcoding, metabolomics, whole-body regeneration (WBR) and aging related pathways.