In multicellular organisms, cell type diversity and fate depend on specific sets of transcript isoforms generated by post-transcriptional RNA processing. Here, we used Schmidtea mediterranea, a flatworm with extraordinary regenerative abilities and a large pool of adult stem cells, as an in vivo model to study the role of Uridyl-rich small nuclear RNAs (UsnRNAs), which participate in multiple RNA processing reactions including splicing, in stem cell regulation. We characterized the planarian UsnRNA repertoire, identified stem cell-enriched variants and obtained strong evidence for an increased rate of UsnRNA 3'-processing in stem cells compared to their differentiated counterparts. Consistently, components of the Integrator complex showed stem cell-enriched expression and their depletion by RNAi disrupted UsnRNA processing resulting in global changes of splicing patterns and reduced processing of histone mRNAs. Interestingly, loss of Integrator complex function disrupted both stem cell maintenance and regeneration of tissues. Our data show that the function of the Integrator complex in UsnRNA 3'-processing is conserved in planarians and essential for maintaining their stem cell pool. We propose that cell type-specific modulation of UsnRNA composition and maturation contributes to in vivo cell fate choices, such as stem cell self-renewal in planarians.

The posttranslational modification of proteins by the ubiquitination pathway is an important regulatory mechanism in eukaryotes. To date, however, studies on the evolutionary history of the proteins involved in this pathway have been restricted to E1 and E2 enzymes, whereas E3 studies have been focused mainly in metazoans and plants. To have a wider perspective, here we perform a genomic survey of the HECT family of E3 ubiquitin-protein ligases, an important part of this posttranslational pathway, in genomes from representatives of all major eukaryotic lineages. We classify eukaryotic HECTs and reconstruct, by phylogenetic analysis, the putative repertoire of these proteins in the last eukaryotic common ancestor (LECA). Furthermore, we analyze the diversity and complexity of protein domain architectures of HECTs along the different extant eukaryotic lineages. Our data show that LECA had six different HECTs and that protein expansion and N-terminal domain diversification shaped HECT evolution. Our data reveal that the genomes of animals and unicellular holozoans considerably increased the molecular and functional diversity of their HECT system compared with other eukaryotes. Other eukaryotes, such as the Apusozoa Thecanomas trahens or the Heterokonta Phytophthora infestans, independently expanded their HECT repertoire. In contrast, plant, excavate, rhodophyte, chlorophyte, and fungal genomes have a more limited enzymatic repertoire. Our genomic survey and phylogenetic analysis clarifies the origin and evolution of different HECT families among eukaryotes and provides a useful phylogenetic framework for future evolutionary studies of this regulatory pathway.

Foxp3+ regulatory T cells (Tregs) are potent suppressor cells, essential for the maintenance of immune homeostasis. Most Tregs develop in the thymus and are then released into the immune periphery. However, some Tregs populate the thymus and constitute a major subset of yet poorly understood cells. Here we describe a subset of thymus recirculating IL18R+ Tregs with molecular characteristics highly reminiscent of tissue-resident effector Tregs. Moreover, we show that IL18R+ Tregs are endowed with higher capacity to populate the thymus than their IL18R- or IL18R-/- counterparts, highlighting the key role of IL18R in this process. Finally, we demonstrate that IL18 signaling is critical for the induction of the key thymus-homing chemokine receptor - CCR6 on Tregs. Collectively, this study provides a detailed characterization of the mature Treg subsets in the mouse thymus and identifies a key role of IL18 signaling in controlling the CCR6-CCL20-dependent migration of Tregs into the thymus.

The discovery of giant viruses infecting eukaryotes from diverse ecosystems has revolutionized our understanding of the evolution of viruses and their impact on protist biology, yet knowledge on their replication strategies and transcriptome regulation remains limited. Here, we profile single-cell transcriptomes of the globally distributed microalga Emiliania huxleyi and its specific giant virus during infection. We detected profound heterogeneity in viral transcript levels among individual cells. Clustering single cells based on viral expression profiles enabled reconstruction of the viral transcriptional trajectory. Reordering cells along this path unfolded highly resolved viral genetic programs composed of genes with distinct promoter elements that orchestrate sequential expression. Exploring host transcriptome dynamics across the viral infection states revealed rapid and selective shutdown of protein-encoding nuclear transcripts, while the plastid and mitochondrial transcriptomes persisted into later stages. Single-cell RNA-seq opens a new avenue to unravel the life cycle of giant viruses and their unique hijacking strategies.

The unicellular ancestor of animals had a complex repertoire of genes linked to multicellular processes. This suggests that changes in the regulatory genome, rather than in gene innovation, were key to the origin of animals. Here, we carry out multiple functional genomic assays in Capsaspora owczarzaki, the unicellular relative of animals with the largest known gene repertoire for transcriptional regulation. We show that changing chromatin states, differential lincRNA expression, and dynamic cis-regulatory sites are associated with life cycle transitions in Capsaspora. Moreover, we demonstrate conservation of animal developmental transcription-factor networks and extensive network interconnection in this premetazoan organism. In contrast, however, Capsaspora lacks animal promoter types, and its regulatory sites are small, proximal, and lack signatures of animal enhancers. Overall, our results indicate that the emergence of animal multicellularity was linked to a major shift in genome cis-regulatory complexity, most notably the appearance of distal enhancer regulation.

Nonaggregative multicellularity requires strict control of cell number. The Hippo signaling pathway coordinates cell proliferation and apoptosis and is a central regulator of organ size in animals. Recent studies have shown the presence of key members of the Hippo pathway in nonbilaterian animals, but failed to identify this pathway outside Metazoa. Through comparative analyses of recently sequenced holozoan genomes, we show that Hippo pathway components, such as the kinases Hippo and Warts, the coactivator Yorkie, and the transcription factor Scalloped, were already present in the unicellular ancestors of animals. Remarkably, functional analysis of Hippo components of the amoeboid holozoan Capsaspora owczarzaki, performed in Drosophila melanogaster, demonstrate that the growth-regulatory activity of the Hippo pathway is conserved in this unicellular lineage. Our findings show that the Hippo pathway evolved well before the origin of Metazoa and highlight the importance of Hippo signaling as a key developmental mechanism predating the origin of Metazoa.

Post-transcriptional regulatory mechanisms are of fundamental importance to form robust genetic networks, but their roles in stem cell pluripotency remain poorly understood. Here, we use freshwater planarians as a model system to investigate this and uncover a role for CCR4-NOT mediated deadenylation of mRNAs in stem cell differentiation. Planarian adult stem cells, the so-called neoblasts, drive the almost unlimited regenerative capabilities of planarians and allow their ongoing homeostatic tissue turnover. While many genes have been demonstrated to be required for these processes, currently almost no mechanistic insight is available into their regulation. We show that knockdown of planarian Not1, the CCR4-NOT deadenylating complex scaffolding subunit, abrogates regeneration and normal homeostasis. This abrogation is primarily due to severe impairment of their differentiation potential. We describe a stem cell specific increase in the mRNA levels of key neoblast genes after Smed-not1 knock down, consistent with a role of the CCR4-NOT complex in degradation of neoblast mRNAs upon the onset of differentiation. We also observe a stem cell specific increase in the frequency of longer poly(A) tails in these same mRNAs, showing that stem cells after Smed-not1 knock down fail to differentiate as they accumulate populations of transcripts with longer poly(A) tails. As other transcripts are unaffected our data hint at a targeted regulation of these key stem cell mRNAs by post-transcriptional regulators such as RNA-binding proteins or microRNAs. Together, our results show that the CCR4-NOT complex is crucial for stem cell differentiation and controls stem cell-specific degradation of mRNAs, thus providing clear mechanistic insight into this aspect of neoblast biology.

The evolution of metazoans from their unicellular ancestors was one of the most important events in the history of life. However, the cellular and genetic changes that ultimately led to the evolution of multicellularity are not known. In this study, we describe an aggregative multicellular stage in the protist Capsaspora owczarzaki, a close unicellular relative of metazoans. Remarkably, transition to the aggregative stage is associated with significant upregulation of orthologs of genes known to establish multicellularity and tissue architecture in metazoans. We further observe transitions in regulated alternative splicing during the C. owczarzaki life cycle, including the deployment of an exon network associated with signaling, a feature of splicing regulation so far only observed in metazoans. Our results reveal the existence of a highly regulated aggregative stage in C. owczarzaki and further suggest that features of aggregative behavior in an ancestral protist may had been co-opted to develop some multicellular properties currently seen in metazoans. DOI: http://dx.doi.org/10.7554/eLife.01287.001.

Filopodia are fine actin-based cellular projections used for both environmental sensing and cell motility, and they are essential organelles for metazoan cells. In this study, we reconstruct the origin of metazoan filopodia and microvilli. We first report on the evolutionary assembly of the filopodial molecular toolkit and show that homologs of many metazoan filopodial components, including fascin and myosin X, were already present in the unicellular or colonial progenitors of metazoans. Furthermore, we find that the actin crosslinking protein fascin localizes to filopodia-like structures and microvilli in the choanoflagellate Salpingoeca rosetta. In addition, homologs of filopodial genes in the holozoan Capsaspora owczarzaki are upregulated in filopodia-bearing cells relative to those that lack them. Therefore, our findings suggest that proteins essential for metazoan filopodia and microvilli are functionally conserved in unicellular and colonial holozoans and that the last common ancestor of metazoans bore a complex and specific filopodial machinery.

Myosins are key components of the eukaryotic cytoskeleton, providing motility for a broad diversity of cargoes. Therefore, understanding the origin and evolutionary history of myosin classes is crucial to address the evolution of eukaryote cell biology. Here, we revise the classification of myosins using an updated taxon sampling that includes newly or recently sequenced genomes and transcriptomes from key taxa. We performed a survey of eukaryotic genomes and phylogenetic analyses of the myosin gene family, reconstructing the myosin toolkit at different key nodes in the eukaryotic tree of life. We also identified the phylogenetic distribution of myosin diversity in terms of number of genes, associated protein domains and number of classes in each taxa. Our analyses show that new classes (i.e., paralogs) and domain architectures were continuously generated throughout eukaryote evolution, with a significant expansion of myosin abundance and domain architectural diversity at the stem of Holozoa, predating the origin of animal multicellularity. Indeed, single-celled holozoans have the most complex myosin complement among eukaryotes, with paralogs of most myosins previously considered animal specific. We recover a dynamic evolutionary history, with several lineage-specific expansions (e.g., the myosin III-like gene family diversification in choanoflagellates), convergence in protein domain architectures (e.g., fungal and animal chitin synthase myosins), and important secondary losses. Overall, our evolutionary scheme demonstrates that the ancestral eukaryote likely had a complex myosin repertoire that included six genes with different protein domain architectures. Finally, we provide an integrative and robust classification, useful for future genomic and functional studies on this crucial eukaryotic gene family.

Recently, a functional IP3R ortholog (CO.IP3R-A) capable of IP3-induced Ca2+ release has been discovered in Capsaspora owczarzaki, a close unicellular relative to Metazoa. In contrast to mammalian IP3Rs, CO.IP3R-A is not modulated by Ca2+, ATP or PKA. Protein-sequence analysis revealed that CO.IP3R-A contained a putative binding site for anti-apoptotic Bcl-2, although Bcl-2 was not detected in Capsaspora owczarzaki and only appeared in Metazoa. Here, we examined whether human Bcl-2 could form a complex with CO.IP3R-A channels and modulate their Ca2+-flux properties using ectopic expression approaches in a HEK293 cell model in which all three IP3R isoforms were knocked out. We demonstrate that human Bcl-2 via its BH4 domain could functionally interact with CO.IP3R-A, thereby suppressing Ca2+ flux through CO.IP3R-A channels. The BH4 domain of Bcl-2 was sufficient for interaction with CO.IP3R-A channels. Moreover, mutating the Lys17 of Bcl-2's BH4 domain, the residue critical for Bcl-2-dependent modulation of mammalian IP3Rs, abrogated Bcl-2's ability to bind and inhibit CO.IP3R-A channels. Hence, this raises the possibility that a unicellular ancestor of animals already had an IP3R that harbored a Bcl-2-binding site. Bcl-2 proteins may have evolved as controllers of IP3R function by exploiting this pre-existing site, thereby counteracting Ca2+-dependent apoptosis.

Stony corals are colonial cnidarians that sustain the most biodiverse marine ecosystems on Earth: coral reefs. Despite their ecological importance, little is known about the cell types and molecular pathways that underpin the biology of reef-building corals. Using single-cell RNA sequencing, we define over 40 cell types across the life cycle of Stylophora pistillata. We discover specialized immune cells, and we uncover the developmental gene expression dynamics of calcium-carbonate skeleton formation. By simultaneously measuring the transcriptomes of coral cells and the algae within them, we characterize the metabolic programs involved in symbiosis in both partners. We also trace the evolution of these coral cell specializations by phylogenetic integration of multiple cnidarian cell type atlases. Overall, this study reveals the molecular and cellular basis of stony coral biology.

The assembly of the neuronal and other major cell type programs occurred early in animal evolution. We can reconstruct this process by studying non-bilaterians like placozoans. These small disc-shaped animals not only have nine morphologically described cell types and no neurons but also show coordinated behaviors triggered by peptide-secreting cells. We investigated possible neuronal affinities of these peptidergic cells using phylogenetics, chromatin profiling, and comparative single-cell genomics in four placozoans. We found conserved cell type expression programs across placozoans, including populations of transdifferentiating and cycling cells, suggestive of active cell type homeostasis. We also uncovered fourteen peptidergic cell types expressing neuronal-associated components like the pre-synaptic scaffold that derive from progenitor cells with neurogenesis signatures. In contrast, earlier-branching animals like sponges and ctenophores lacked this conserved expression. Our findings indicate that key neuronal developmental and effector gene modules evolved before the advent of cnidarian/bilaterian neurons in the context of paracrine cell signaling.

Planarian adult stem cells (pASCs) or neoblasts represent an ideal system to study the evolution of stem cells and pluripotency as they underpin an unrivaled capacity for regeneration. We wish to understand the control of differentiation and pluripotency in pASCs and to understand how conserved, convergent or divergent these mechanisms are across the Bilateria. Here we show the planarian methyl-CpG Binding Domain 2/3 (mbd2/3) gene is required for pASC differentiation during regeneration and tissue homeostasis. The genome does not have detectable levels of 5-methylcytosine (5(m)C) and we find no role for a potential DNA methylase. We conclude that MBD proteins may have had an ancient role in broadly controlling animal stem cell pluripotency, but that DNA methylation is not involved in planarian stem cell differentiation.

The origin of the eukaryotic cell is one of the most important transitions in the history of life. However, the emergence and early evolution of eukaryotes remains poorly understood. Recent data have shown that the last eukaryotic common ancestor (LECA) was much more complex than previously thought. The LECA already had the genetic machinery encoding the endomembrane apparatus, spliceosome, nuclear pore, and myosin and kinesin cytoskeletal motors. It is unclear, however, when the functional regulation of these cellular components evolved. Here, we address this question by analyzing the origin and evolution of the ubiquitin (Ub) signaling system, one of the most important regulatory layers in eukaryotes. We delineated the evolution of the whole Ub, Small-Ub-related MOdifier (SUMO), and Ub-fold modifier 1 (Ufm1) signaling networks by analyzing representatives from all major eukaryotic, bacterial, and archaeal lineages. We found that the Ub toolkit had a pre-eukaryotic origin and is present in three extant archaeal groups. The pre-eukaryotic Ub toolkit greatly expanded during eukaryogenesis, through massive gene innovation and diversification of protein domain architectures. This resulted in a LECA with essentially all of the Ub-related genes, including the SUMO and Ufm1 Ub-like systems. Ub and SUMO signaling further expanded during eukaryotic evolution, especially labeling and delabeling enzymes responsible for substrate selection. Additionally, we analyzed protein domain architecture evolution and found that multicellular lineages have the most complex Ub systems in terms of domain architectures. Together, we demonstrate that the Ub system predates the origin of eukaryotes and that a burst of innovation during eukaryogenesis led to a LECA with complex posttranslational regulation.

The emergence and diversification of cell types is a leading factor in animal evolution. So far, systematic characterization of the gene regulatory programs associated with cell type specificity was limited to few cell types and few species. Here, we perform whole-organism single-cell transcriptomics to map adult and larval cell types in the cnidarian Nematostella vectensis, a non-bilaterian animal with complex tissue-level body-plan organization. We uncover eight broad cell classes in Nematostella, including neurons, cnidocytes, and digestive cells. Each class comprises different subtypes defined by the expression of multiple specific markers. In particular, we characterize a surprisingly diverse repertoire of neurons, which comparative analysis suggests are the result of lineage-specific diversification. By integrating transcription factor expression, chromatin profiling, and sequence motif analysis, we identify the regulatory codes that underlie Nematostella cell-specific expression. Our study reveals cnidarian cell type complexity and provides insights into the evolution of animal cell-specific genomic regulation.

To reconstruct the evolutionary origin of multicellular animals from their unicellular ancestors, the genome sequences of diverse unicellular relatives are essential. However, only the genome of the choanoflagellate Monosiga brevicollis has been reported to date. Here we completely sequence the genome of the filasterean Capsaspora owczarzaki, the closest known unicellular relative of metazoans besides choanoflagellates. Analyses of this genome alter our understanding of the molecular complexity of metazoans' unicellular ancestors showing that they had a richer repertoire of proteins involved in cell adhesion and transcriptional regulation than previously inferred only with the choanoflagellate genome. Some of these proteins were secondarily lost in choanoflagellates. In contrast, most intercellular signalling systems controlling development evolved later concomitant with the emergence of the first metazoans. We propose that the acquisition of these metazoan-specific developmental systems and the co-option of pre-existing genes drove the evolutionary transition from unicellular protists to metazoans.

Eukaryotic microbes have three primary mechanisms for obtaining nutrients and energy: phagotrophy, photosynthesis and osmotrophy. Traits associated with the latter two functions arose independently multiple times in the eukaryotes. The Fungi successfully coupled osmotrophy with filamentous growth, and similar traits are also manifested in the Pseudofungi (oomycetes and hyphochytriomycetes). Both the Fungi and the Pseudofungi encompass a diversity of plant and animal parasites. Genome-sequencing efforts have focused on host-associated microbes (mutualistic symbionts or parasites), providing limited comparisons with free-living relatives. Here we report the first draft genome sequence of a hyphochytriomycete 'pseudofungus'; Hyphochytrium catenoides Using phylogenomic approaches, we identify genes of recent viral ancestry, with related viral derived genes also present on the genomes of oomycetes, suggesting a complex history of viral coevolution and integration across the Pseudofungi. H. catenoides has a complex life cycle involving diverse filamentous structures and a flagellated zoospore with a single anterior tinselate flagellum. We use genome comparisons, drug sensitivity analysis and high-throughput culture arrays to investigate the ancestry of oomycete/pseudofungal characteristics, demonstrating that many of the genetic features associated with parasitic traits evolved specifically within the oomycete radiation. Comparative genomics also identified differences in the repertoire of genes associated with filamentous growth between the Fungi and the Pseudofungi, including differences in vesicle trafficking systems, cell-wall synthesis pathways and motor protein repertoire, demonstrating that unique cellular systems underpinned the convergent evolution of filamentous osmotrophic growth in these two eukaryotic groups.

Sequence variants in cis-acting enhancers are important for polygenic disease, but their role in Mendelian disease is poorly understood. Redundancy between enhancers that regulate the same gene is thought to mitigate the pathogenic impact of enhancer mutations. Recent findings, however, have shown that loss-of-function mutations in a single enhancer near PTF1A cause pancreas agenesis and neonatal diabetes. Using mouse and human genetic models, we show that this enhancer activates an entire PTF1A enhancer cluster in early pancreatic multipotent progenitors. This leading role, therefore, precludes functional redundancy. We further demonstrate that transient expression of PTF1A in multipotent progenitors sets in motion an epigenetic cascade that is required for duct and endocrine differentiation. These findings shed insights into the genome regulatory mechanisms that drive pancreas differentiation. Furthermore, they reveal an enhancer that acts as a regulatory master key and is thus vulnerable to pathogenic loss-of-function mutations.

In contrast to transcriptional regulation, the function of alternative splicing (AS) in stem cells is poorly understood. In mammals, MBNL proteins negatively regulate an exon program specific of embryonic stem cells; however, little is known about the in vivo significance of this regulation. We studied AS in a powerful in vivo model for stem cell biology, the planarian Schmidtea mediterranea. We discover a conserved AS program comprising hundreds of alternative exons, microexons and introns that is differentially regulated in planarian stem cells, and comprehensively identify its regulators. We show that functional antagonism between CELF and MBNL factors directly controls stem cell-specific AS in planarians, placing the origin of this regulatory mechanism at the base of Bilaterians. Knockdown of CELF or MBNL factors lead to abnormal regenerative capacities by affecting self-renewal and differentiation sets of genes, respectively. These results highlight the importance of AS interactions in stem cell regulation across metazoans.