The process of speciation, through which two reproductively isolated species emerge from an ancestral population, is often gradual, including a period during which gene flow still occurs, and hence speciation remains reversible. However, if isolation persists, the progressive accumulation of genetic incompatibilities ultimately locks the speciation process, reaching the point where hybrid offspring are inviable or sterile. In plants and animals, the establishment of species barriers is quite rapid: beyond a net genetic divergence of 0.02 substitutions per synonymous site, interspecific introgressions become extremely rare. Here we report that in ciliates, species barriers can be broken despite huge genetic divergence. We show that the genome of Paramecium sonneborni contains hundreds of DNA segments that were acquired by recurrent horizontal transfers from multiple very distant lineages (synonymous divergence ∼1 substitution/site). These segments approach the size of full-length chromosomes, which indicates that they were introgressed by inter-specific mating. These data show that in the ancestry of P. sonneborni , hybridization events with highly divergent species occurred repeatedly and led to viable and fertile offspring. We propose that this unexpected promiscuous sexual behavior is made possible by the fact that paramecia separate somatic and germline functions in two types of nuclei: the elimination of non-self DNA fragments during the development of the somatic genome allows individuals of a species to remain phenotypically the same, despite recurrent genetic flow from distant species into their germline.

Paramecia belong to the ciliate phylum of unicellular eukaryotes characterized by nuclear dimorphism. A diploid germline micronucleus (MIC) transmits genetic information across sexual generations. A polyploid transcriptionally active somatic macronucleus (MAC) develops at each sexual generation from a copy of the MIC through programmed DNA elimination (PDE) of > 30

Sexual reproduction in eukaryotes typically involves the merging of gametes, cells, or nuclei belonging to different mating classes, such as small and large gametes, mating types, or sexes. In this review, we combine theoretical results, recent genomic data, and documented changes in mating systems among ciliate, fungal, plant, algal, and animal lineages to identify general patterns in the evolution of mating classes and their underlying genetic architectures. We emphasize how interactions among distinct coexisting "layers" of mating classes within a species (e.g., anisogamy plus self-incompatibility systems) have shaped their evolution. We raise key unresolved questions, including why no species appears to maintain more than three coexisting layers of mating classes.

Whole-genome duplications (WGDs) have shaped the gene repertoire of many eukaryotic lineages. The redundancy created by WGDs typically results in a phase of massive gene loss. However, some WGD-derived paralogs are maintained over long evolutionary periods, and the relative contributions of different selective pressures to their maintenance are still debated. Previous studies have revealed a history of three successive WGDs in the lineage of the ciliate Paramecium tetraurelia and two of its sister species from the Paramecium aurelia complex. Here, we report the genome sequence and analysis of 10 additional P. aurelia species and 1 additional out group, revealing aspects of post-WGD evolution in 13 species sharing a common ancestral WGD. Contrary to the morphological radiation of vertebrates that putatively followed two WGD events, members of the cryptic P. aurelia complex have remained morphologically indistinguishable after hundreds of millions of years. Biases in gene retention compatible with dosage constraints appear to play a major role opposing post-WGD gene loss across all 13 species. In addition, post-WGD gene loss has been slower in Paramecium than in other species having experienced genome duplication, suggesting that the selective pressures against post-WGD gene loss are especially strong in Paramecium. A near complete lack of recent single-gene duplications in Paramecium provides additional evidence for strong selective pressures against gene dosage changes. This exceptional data set of 13 species sharing an ancestral WGD and 2 closely related out group species will be a useful resource for future studies on Paramecium as a major model organism in the evolutionary cell biology.

The Paramecium aurelia complex, a group of morphologically similar but sexually incompatible sibling species, is a unique example of the evolutionary plasticity of mating-type systems. Each species has two mating types, O (Odd) and E (Even). Although O and E types are homologous in all species, three different modes of determination and inheritance have been described: genetic determination by Mendelian alleles, stochastic developmental determination, and maternally inherited developmental determination. Previous work in three species of the latter kind has revealed the key roles of the E-specific transmembrane protein mtA and its highly specific transcription factor mtB: type O clones are produced by maternally inherited genome rearrangements that inactivate either mtA or mtB during development. Here we show, through transcriptome analyses in five additional species representing the three determination systems, that mtA expression specifies type E in all cases. We further show that the Mendelian system depends on functional and nonfunctional mtA alleles, and identify novel developmental rearrangements in mtA and mtB which now explain all cases of maternally inherited mating-type determination. Epistasis between these genes likely evolved from less specific interactions between paralogs in the P. aurelia common ancestor, after a whole-genome duplication, but the mtB gene was subsequently lost in three P. aurelia species which appear to have returned to an ancestral regulation mechanism. These results suggest a model accounting for evolutionary transitions between determination systems, and highlight the diversity of molecular solutions explored among sibling species to maintain an essential mating-type polymorphism in cell populations.

Ciliates are unicellular eukaryotes with both a germline genome and a somatic genome in the same cytoplasm. The somatic macronucleus (MAC), responsible for gene expression, is not sexually transmitted but develops from a copy of the germline micronucleus (MIC) at each sexual generation. In the MIC genome of Paramecium tetraurelia, genes are interrupted by tens of thousands of unique intervening sequences called internal eliminated sequences (IESs), which have to be precisely excised during the development of the new MAC to restore functional genes. To understand the evolutionary origin of this peculiar genomic architecture, we sequenced the MIC genomes of 9 Paramecium species (from approximately 100 Mb in Paramecium aurelia species to >1.5 Gb in Paramecium caudatum). We detected several waves of IES gains, both in ancestral and in more recent lineages. While the vast majority of IESs are single copy in present-day genomes, we identified several families of mobile IESs, including nonautonomous elements acquired via horizontal transfer, which generated tens to thousands of new copies. These observations provide the first direct evidence that transposable elements can account for the massive proliferation of IESs in Paramecium. The comparison of IESs of different evolutionary ages indicates that, over time, IESs shorten and diverge rapidly in sequence while they acquire features that allow them to be more efficiently excised. We nevertheless identified rare cases of IESs that are under strong purifying selection across the aurelia clade. The cases examined contain or overlap cellular genes that are inactivated by excision during development, suggesting conserved regulatory mechanisms. Similar to the evolution of introns in eukaryotes, the evolution of Paramecium IESs highlights the major role played by selfish genetic elements in shaping the complexity of genome architecture and gene expression.

AbstractThe Paramecium aurelia complex, a group of morphologically similar but sexually incompatible sibling species, is a unique example of the evolutionary plasticity of mating-type systems. Each species has two mating types, O (Odd) and E (Even). Although O and E types are homologous in all species, three different modes of determination and inheritance have been described: genetic determination by Mendelian alleles, stochastic developmental determination, and maternally inherited developmental determination. Previous work in three species of the latter kind has revealed the key roles of the E-specific transmembrane protein mtA and its highly specific transcription factor mtB: type O clones are produced by maternally inherited genome rearrangements that inactivate either mtA or mtB during development. Here we show, through transcriptome analyses in five additional species representing the three determination systems, that mtA expression specifies type E in all cases. We further show that the Mendelian system depends on functional and nonfunctional mtA alleles, and identify novel developmental rearrangements in mtA and mtB which now explain all cases of maternally inherited mating-type determination. Epistasis between these genes likely evolved from less specific interactions between paralogs in the P. aurelia common ancestor, after a whole-genome duplication, but the mtB gene was subsequently lost in three P. aurelia species which appear to have returned to an ancestral regulation mechanism. These results suggest a model accounting for evolutionary transitions between determination systems, and highlight the diversity of molecular solutions explored among sibling species to maintain an essential mating-type polymorphism in cell populations.

Whole-Genome Duplications (WGDs) have shaped the gene repertoire of many eukaryotic lineages. The redundancy created by WGDs typically results in a phase of massive gene loss. However, some WGD-derived paralogs are maintained over long evolutionary periods and the relative contributions of different selective pressures to their maintenance is still debated. Previous studies have revealed a history of three successive WGDs in the lineage of the ciliate Paramecium tetraurelia and two of its sister species from the P. aurelia complex. Here, we report the genome sequence and analysis of 10 additional P. aurelia species and one additional outgroup, allowing us to track post-WGD evolution in 13 species that share a common ancestral WGD. We found similar biases in gene retention compatible with dosage constraints playing a major role opposing post-WGD gene loss across all 13 species. Interestingly we found that post-WGD gene loss was slower in Paramecium than in other species having experienced genome duplication, suggesting that the selective pressures against post-WGD gene loss are especially strong in Paramecium . We also report a lack of recent segmental duplications in Paramecium , which we interpret as additional evidence for strong selective pressures against individual genes dosage changes. Finally, we hope that this exceptional dataset of 13 species sharing an ancestral WGD and two closely related outgroup species will be a useful resource for future studies and will help establish Paramecium as a major model organism in the study of post-WGD evolution.

A conspicuous cell-shape phenotype known as "screwy" was reported to result from mutations at two or three uncharacterized loci in the ciliate Paramecium tetraurelia. Here, we describe a new screwy mutation, Spinning Top, which appeared spontaneously in the cross of an unrelated mutant with reference strain 51. The macronuclear (MAC) genome of the Spinning Top mutant is shown to lack a ~28.5-kb segment containing 18 genes at the end of one chromosome, which appears to result from a collinear deletion in the micronuclear (MIC) genome. We tested several candidate genes from the deleted locus by dsRNA-induced silencing in wild-type cells, and identified a single gene responsible for the phenotype. This gene, named Spade, encodes a 566-aa glutamine-rich protein with a C2HC zinc finger. Its silencing leads to a fast phenotype switch during vegetative growth, but cells recover a wild-type phenotype only 5-6 divisions after silencing is stopped. We analyzed 5 independently-obtained mutant alleles of the Sc1 locus, and concluded that all of them also lack the Spade gene and a number of neighboring genes in the MAC and MIC genomes. Mapping of the MAC deletion breakpoints revealed two different positions among the 5 alleles, both of which differ from the Spinning Top breakpoint. These results suggest that this MIC chromosome region is intrinsically unstable in strain 51.

ParameciumDB (https://paramecium.i2bc.paris-saclay.fr) is a community model organism database for the genome and genetics of the ciliate Paramecium. ParameciumDB development relies on the GMOD (www.gmod.org) toolkit. The ParameciumDB web site has been publicly available since 2006 when the P. tetraurelia somatic genome sequence was released, revealing that a series of whole genome duplications punctuated the evolutionary history of the species. The genome is linked to available genetic data and stocks. ParameciumDB has undergone major changes in its content and website since the last update published in 2011. Genomes from multiple Paramecium species, especially from the P. aurelia complex, are now included in ParameciumDB. A new modern web interface accompanies this transition to a database for the whole Paramecium genus. Gene pages have been enriched with orthology relationships, among the Paramecium species and with a panel of model organisms across the eukaryotic tree. This update also presents expert curation of Paramecium mitochondrial genomes.

In the ciliate Paramecium tetraurelia, functional genes are reconstituted during development of the somatic macronucleus through the precise excision of ∼45 000 single-copy Internal Eliminated Sequences (IESs), thought to be the degenerate remnants of ancient transposon insertions. Like introns, IESs are marked only by a weak consensus at their ends. How such a diverse set of sequences is faithfully recognized and precisely excised remains unclear: specialized small RNAs have been implicated, but in their absence up to ∼60% of IESs are still correctly excised. To get further insight, we designed a mutagenesis screen based on the hypersensitivity of a specific excision event in the mtA gene, which determines mating types. Unlike most IES-containing genes, the active form of mtA is the unexcised one, allowing the recovery of hypomorphic alleles of essential IES recognition/excision factors. Such is the case of one mutation recovered in the Piwi gene PTIWI09, a key player in small RNA-mediated IES recognition. Another mutation identified a novel protein with a C2H2 zinc finger, mtGa, which is required for excision of a small subset of IESs characterized by enrichment in a 5-bp motif. The unexpected implication of a sequence-specific factor establishes a new paradigm for IES recognition and/or excision.

Background: Most eukaryotic genes are subject to alternative splicing (AS), which may contribute to the production of protein variants or to the regulation of gene expression via nonsense-mediated messenger RNA (mRNA) decay (NMD). However, a fraction of splice variants might correspond to spurious transcripts and the question of the relative proportion of splicing errors to functional splice variants remains highly debated. Results: We propose a test to quantify the fraction of AS events corresponding to errors. This test is based on the fact that the fitness cost of splicing errors increases with the number of introns in a gene and with expression level. We analyzed the transcriptome of the intron-rich eukaryote Paramecium tetraurelia. We show that in both normal and in NMD-deficient cells, AS rates strongly decrease with increasing expression level and with increasing number of introns. This relationship is observed for AS events that are detectable by NMD as well as for those that are not, which invalidates the hypothesis of a link with the regulation of gene expression. Our results show that in genes with a median expression level, 92-98% of observed splice variants correspond to errors. We observed the same patterns in human transcriptomes and we further show that AS rates correlate with the fitness cost of splicing errors. Conclusions: These observations indicate that genes under weaker selective pressure accumulate more maladaptive substitutions and are more prone to splicing errors. Thus, to a large extent, patterns of gene expression variants simply reflect the balance between selection, mutation, and drift.

While sex is an ancient and highly conserved eukaryotic invention, self-incompatibility systems such as mating types or sexes appear to be derived limitations that show considerable evolutionary plasticity. Within a single class of ciliates, Paramecium and Tetrahymena species have long been known to present a wide variety of mating type numbers and modes of inheritance, but only recently have the genes involved been identified. Although similar transmembrane proteins mediate self/nonself recognition in both ciliates, the mechanisms of mating type determination differ widely, ranging from Mendelian systems to developmental nuclear differentiation, either stochastic or maternally inherited. The non-Mendelian systems rely on programmed editing of the germline genome that occurs during differentiation of the somatic nucleus, and they have co-opted different DNA recombination mechanisms—some previously unknown. Here we review the recent molecular advances and some remaining unsolved questions and discuss the possible implications of these diverse mechanisms for inbreeding/outbreeding balance regulation.

BACKGROUND:DNA elimination is developmentally programmed in a wide variety of eukaryotes, including unicellular ciliates, and leads to the generation of distinct germline and somatic genomes. The ciliate Paramecium tetraurelia harbors two types of nuclei with different functions and genome structures. The transcriptionally inactive micronucleus contains the complete germline genome, while the somatic macronucleus contains a reduced genome streamlined for gene expression. During development of the somatic macronucleus, the germline genome undergoes massive and reproducible DNA elimination events. Availability of both the somatic and germline genomes is essential to examine the genome changes that occur during programmed DNA elimination and ultimately decipher the mechanisms underlying the specific removal of germline-limited sequences.RESULTS:We developed a novel experimental approach that uses flow cell imaging and flow cytometry to sort subpopulations of nuclei to high purity. We sorted vegetative micronuclei and macronuclei during development of P. tetraurelia. We validated the method by flow cell imaging and by high throughput DNA sequencing. Our work establishes the proof of principle that developing somatic macronuclei can be sorted from a complex biological sample to high purity based on their size, shape and DNA content. This method enabled us to sequence, for the first time, the germline DNA from pure micronuclei and to identify novel transposable elements. Sequencing the germline DNA confirms that the Pgm domesticated transposase is required for the excision of all ~45,000 Internal Eliminated Sequences. Comparison of the germline DNA and unrearranged DNA obtained from PGM-silenced cells reveals that the latter does not provide a faithful representation of the germline genome.CONCLUSIONS:We developed a flow cytometry-based method to purify P. tetraurelia nuclei to high purity and provided quality control with flow cell imaging and high throughput DNA sequencing. We identified 61 germline transposable elements including the first Paramecium retrotransposons. This approach paves the way to sequence the germline genomes of P. aurelia sibling species for future comparative genomic studies.

In various organisms, an efficient RNAi response can be triggered by feeding cells with bacteria producing double-stranded RNA (dsRNA) against an endogenous gene. However, the detailed mechanisms and natural functions of this pathway are not well understood in most cases. Here, we studied siRNA biogenesis from exogenous RNA and its genetic overlap with endogenous RNAi in the ciliate Paramecium tetraurelia by high-throughput sequencing. Using wild-type and mutant strains deficient for dsRNA feeding we found that high levels of primary siRNAs of both strands are processed from the ingested dsRNA trigger by the Dicer Dcr1, the RNA-dependent RNA polymerases Rdr1 and Rdr2 and other factors. We further show that this induces the synthesis of secondary siRNAs spreading along the entire endogenous mRNA, demonstrating the occurrence of both 3'-to-5' and 5'-to-3' transitivity for the first time in the SAR clade of eukaryotes (Stramenopiles, Alveolates, Rhizaria). Secondary siRNAs depend on Rdr2 and show a strong antisense bias; they are produced at much lower levels than primary siRNAs and hardly contribute to RNAi efficiency. We further provide evidence that the Paramecium RNAi machinery also processes single-stranded RNAs from its bacterial food, broadening the possible natural functions of exogenously induced RNAi in this organism.

Understanding the links between genetic, epigenetic and non-genetic factors throughout the lifespan and across generations and their role in disease susceptibility and disease progression offer entirely new avenues and solutions to major problems in our society. To overcome the numerous challenges, we have come up with nine major conclusions to set the vision for future policies and research agendas at the European level.

In most eukaryotes, small RNA-mediated gene silencing pathways form complex interacting networks. In the ciliate Paramecium tetraurelia, at least two RNA interference (RNAi) mechanisms coexist, involving distinct but overlapping sets of protein factors and producing different types of short interfering RNAs (siRNAs). One is specifically triggered by high-copy transgenes, and the other by feeding cells with double-stranded RNA (dsRNA)-producing bacteria. In this study, we designed a forward genetic screen for mutants deficient in dsRNA-induced silencing, and a powerful method to identify the relevant mutations by whole-genome sequencing. We present a set of 47 mutant alleles for five genes, revealing two previously unknown RNAi factors: a novel Paramecium-specific protein (Pds1) and a Cid1-like nucleotidyl transferase. Analyses of allelic diversity distinguish non-essential and essential genes and suggest that the screen is saturated for non-essential, single-copy genes. We show that non-essential genes are specifically involved in dsRNA-induced RNAi while essential ones are also involved in transgene-induced RNAi. One of the latter, the RNA-dependent RNA polymerase RDR2, is further shown to be required for all known types of siRNAs, as well as for sexual reproduction. These results open the way for the dissection of the genetic complexity, interconnection, mechanisms and natural functions of RNAi pathways in P. tetraurelia.