Dynamics of Gene Loss following Ancient Whole-Genome Duplication in the Cryptic Paramecium Complex.

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.

Paramecium (Oligohymenophorea, Ciliophora) diversity in Thailand sheds light on the genus biogeography and reveals new phylogenetic lineages.

Paramecium (Ciliophora, Oligohymenophorea) is a good model to study ciliate biogeography. Extensive sampling mainly in northern hemisphere has led to 16 valid morphological species description thus far. However, a majority of hard-to-reach regions, including South East Asia, are underinvestigated. Our study combined traditional morphological and molecular approaches to reveal the biodiversity of Paramecium in Thailand from more than 110 samples collected in 10 provinces. Representatives of seven morphological species were identified from our collection, including the rare species, such as P. gigas and P. jenningsi. Additionally, we detected five different sibling species of the P. aurelia complex, described a new cryptic species P. hiwatashii n. sp. phylogenetically related to P. caudatum, and discovered a potentially new genetic species of the P. bursaria species complex. We also documented a variety of bacterial cytoplasmic symbionts from at least nine monoclonal cultures of Paramecium.

'Candidatus Sarmatiella mevalonica' endosymbiont of the ciliate Paramecium provides insights on evolutionary plasticity among Rickettsiales.

Members of the bacterial order Rickettsiales are obligatorily associated with a wide range of eukaryotic hosts. Their evolutionary trajectories, in particular concerning the origin of shared or differential traits among distant sub-lineages, are still poorly understood. Here, we characterized a novel Rickettsiales bacterium associated with the ciliate Paramecium tredecaurelia and phylogenetically related to the Rickettsia genus. Its genome encodes significant lineage-specific features, chiefly the mevalonate pathway gene repertoire, involved in isoprenoid precursor biosynthesis. Not only this pathway has never been described in Rickettsiales, it also is very rare among bacteria, though typical in eukaryotes, thus likely representing a horizontally acquired trait. The presence of these genes could enable an efficient exploitation of host-derived intermediates for isoprenoid synthesis. Moreover, we hypothesize the reversed reactions could have replaced canonical pathways for producing acetyl-CoA, essential for phospholipid biosynthesis. Additionally, we detected phylogenetically unrelated mevalonate pathway genes in metagenome-derived Rickettsiales sequences, likely indicating evolutionary convergent effects of independent horizontal gene transfer events. Accordingly, convergence, involving both gene acquisitions and losses, is highlighted as a relevant evolutionary phenomenon in Rickettsiales, possibly favoured by plasticity and comparable lifestyles, representing a potentially hidden origin of other more nuanced similarities among sub-lineages.

Genetic diversity of Paramecium species on the basis of multiple loci analysis and ITS secondary structure models.

Among ciliates, Paramecium has become a privileged model for the study of "species problem" particularly in the case of the "Paramecium aurelia complex" that has been intensely investigated. Despite extensive studies, the taxonomy of Paramecium is still challenging. The major problem is an uneven sampling of Paramecium with relatively few representatives of each species. To investigate species from the less discovered region (Pakistan), 10 isolates of Paramecium species including a standing-alone FT8 strain previously isolated by some of us were subjected to molecular characterization. Fragments of 18S recombinant DNA (rDNA), ITS1-5.8S-ITS2-5'LSU rDNA, cytochrome c oxidase subunit II, and hsp70 genes were used as molecular markers for phylogenetic analysis of particular isolates. The nucleotide sequences of polymerase chain reaction products of all markers were compared with the available sequences of relevant markers of other Paramecium species from GenBank. Phylogenetic trees based on all molecular markers showed that all the nine strains had a very close relationship with Paramecium primaurelia except for the FT8 strain. FT8 consistently showed its unique position in comparison to all other species in the phylogenetic trees. Available sequences of internal transcribed spacer 1 (ITS1) and ITS2 and some other ciliate sequences from GenBank were used for the construction of secondary models. Two highly conserved helices supported by compensatory base changes among all ciliates of ITS2 secondary structures were found similar to other eukaryotes. Therefore, the most conserved 120 to 180 base pairs regions were identified for their comparative studies. We found that out of the three helices in ITS1 structure, helix B was more conserved in Paramecium species. Despite various substitutions in the primary sequence, it was observed that secondary structures of ITS1 and ITS2 could be helpful in interpreting the phylogenetic relationships both at species as well as at generic level.

The first case of microsporidiosis in Paramecium.

A new microsporidian species, Globosporidium paramecii gen. nov., sp. nov., from Paramecium primaurelia is described on the basis of morphology, fine structure, and SSU rRNA gene sequence. This is the first case of microsporidiosis in Paramecium reported so far. All observed stages of the life cycle are monokaryotic. The parasites develop in the cytoplasm, at least some part of the population in endoplasmic reticulum and its derivates. Meronts divide by binary fission. Sporogonial plasmodium divides by rosette-like budding. Early sporoblasts demonstrate a well-developed exospore forming blister-like structures. Spores with distinctive spherical shape are dimorphic in size (3.7 ± 0.2 and 1.9 ± 0.2 µm). Both types of spores are characterized by a thin endospore, a short isofilar polar tube making one incomplete coil, a bipartite polaroplast, and a large posterior vacuole. Experimental infection was successful for 5 of 10 tested strains of the Paramecium aurelia species complex. All susceptible strains belong to closely related P. primaurelia and P. pentaurelia species. Phylogenetic analysis placed the new species in the Clade 4 of Microsporidia and revealed its close relationship to Euplotespora binucleata (a microsporidium from the ciliate Euplotes woodruffi), to Helmichia lacustris and Mrazekia macrocyclopis, microsporidia from aquatic invertebrates.

Genome-defence small RNAs exapted for epigenetic mating-type inheritance.

In the ciliate Paramecium, transposable elements and their single-copy remnants are deleted during the development of somatic macronuclei from germline micronuclei, at each sexual generation. Deletions are targeted by scnRNAs, small RNAs produced from the germ line during meiosis that first scan the maternal macronuclear genome to identify missing sequences, and then allow the zygotic macronucleus to reproduce the same deletions. Here we show that this process accounts for the maternal inheritance of mating types in Paramecium tetraurelia, a long-standing problem in epigenetics. Mating type E depends on expression of the transmembrane protein mtA, and the default type O is determined during development by scnRNA-dependent excision of the mtA promoter. In the sibling species Paramecium septaurelia, mating type O is determined by coding-sequence deletions in a different gene, mtB, which is specifically required for mtA expression. These independently evolved mechanisms suggest frequent exaptation of the scnRNA pathway to regulate cellular genes and mediate transgenerational epigenetic inheritance of essential phenotypic polymorphisms.

High-Throughput Sequencing of the 16S rRNA Gene as a Survey to Analyze the Microbiomes of Free-Living Ciliates Paramecium.

Ciliates are the largest group of ubiquitous aquatic bacterivorous protists, and many species are easily cultivated. However, only few studies reported prokaryotic communities naturally associated with ciliate cells. Herein, we analyzed the microbiome composition of several strains of Paramecium (Ciliophora) originating from different locations and belonging to two morpho-species by high-throughput sequencing (HTS) of the 16S rRNA gene. Possible reasons of HTS results bias were addressed comparing DNA libraries obtained using different primers and different number of ciliate cells. Microbiomes associated with ciliates and their environments were always significantly different by prokaryotic taxonomic composition and bacterial richness. There were also pronounced differences between Paramecium strains. Interestingly, potentially pathogenic bacteria were revealed in Paramecium microbiomes.

Evaluation of Enrichment Protocols for Bacterial Endosymbionts of Ciliates by Real-Time PCR.

Large-scale studies on obligate bacterial endosymbionts may frequently require preliminary purification and enrichment protocols, which are often elaborate to set up and to evaluate, especially if the host organism is a protist. The purpose of this study was to develop a real-time PCR-based strategy and employ it for assessing two of such enrichment protocols for Holospora caryophila, hosted by the ciliate Paramecium. Four SSU rRNA gene-targeted real-time PCR assays were designed, which allowed to compare the amount of H. caryophila to other organisms, namely the host, its food bacterium (Raoultella planticola), and free-living bacteria present in the culture medium. By the use of the real-time PCR assays in combination, it was possible to conclude that the "cell fractionation" protocol was quite successful in the enrichment of the symbiont, while the "Percoll gradient" protocol will need further refinements to be fully repeatable. The proposed approach has the potential to facilitate and encourage future studies on the yet underexplored field of bacterial endosymbionts of ciliates and other protists. It can also find valuable applications for experimental questions other than those tested, such as fast and precise assessment of symbiont abundance in natural populations and comparison among multiple coexisting symbionts.

The size of DNA molecules and chromatin organization in the macronucleus of the ciliate Didinium nasutum (Ciliophora).

Pulsed-field gel electrophoresis (PFGE) was applied to analyze the molecular karyotype of the ciliate Didinium nasutum. The data obtained indicate that D. nasutum belongs to the ciliate species with subchromosomal macronuclear genome organization. No short "gene-sized" DNA molecules were detected. Macronuclear DNAs formed a continuous spectrum from 50 kbp to approximately 1,000 kbp in size with a peak plateau between 250 and 400 kbp. The macronuclear DNA molecules were packed into chromatin bodies of 80-265 nm in size. Comparison of the PFGE and electron microscopic data shows that most if not all chromatin bodies contain more than one DNA molecule.

Paramecium putrinum (Ciliophora, Protozoa): the first insight into the variation of two DNA fragments - molecular support for the existence of cryptic species.

Paramecium putrinum (Claparede & Lachmann 1858) is one of the smallest (80-140 µm long) species of the genus Paramecium. Although it commonly occurs in freshwater reservoirs, no molecular studies of P. putrinum have been conducted to date. Herein we present an assessment of molecular variation in 27 strains collected from widely separated populations by using two selected DNA fragments (ITS1-5.8S-ITS2-5'LSU rDNA and COI mtDNA). Both the trees and haplotype networks reconstructed for both genome fragments show that the studied strains of P. putrinum form five main haplogroups. The mean distance between the studied strains is p-distance=0.007/0.068 (rDNA/COI) and exhibits similar variability as that between P. bursaria syngens. Based on these data, one could hypothesize that the clusters revealed in the present study may correspond to previously reported syngens and that there are at least five cryptic species within P. putrinum.

Host association and intracellularity evolved multiple times independently in the Rickettsiales.

The order Rickettsiales (Alphaproteobacteria) encompasses multiple diverse lineages of host-associated bacteria, including pathogens, reproductive manipulators, and mutualists. Here, in order to understand how intracellularity and host association originated in this order, and whether they are ancestral or convergently evolved characteristics, we built a large and phylogenetically-balanced dataset that includes de novo sequenced genomes and a selection of published genomic and metagenomic assemblies. We perform detailed functional reconstructions that clearly indicates "late" and parallel evolution of obligate host-association in different Rickettsiales lineages. According to the depicted scenario, multiple independent horizontal acquisitions of transporters led to the progressive loss of biosynthesis of nucleotides, amino acids and other metabolites, producing distinct conditions of host-dependence. Each clade experienced a different pattern of evolution of the ancestral arsenal of interaction apparatuses, including development of specialised effectors involved in the lineage-specific mechanisms of host cell adhesion and/or invasion.

Holospora-like bacteria "Candidatus Gortzia yakutica" and Preeria caryophila: Ultrastructure, promiscuity, and biogeography of the symbionts.

Two already known representatives of Holospora-like bacteria, "Candidatus Gortzia yakutica" from Paramecium putrinum and Preeria caryophila, originally retrieved from the Paramecium aurelia complex, were found in new hosts: Paramecium nephridiatum and Paramecium polycaryum, respectively. In the present study, these bacteria were investigated using morphological and molecular methods. For "Ca. G. yakutica", the first details of the electron microscopic structure in the main and new hosts were provided. Regarding Pr. caryophila, the ultrastructural description of this species was implemented by several features previously unknown, such as the so called "membrane cluster" dividing periplasm from cytoplasm and fine composition of infectious forms before and during its releasing from the infected macronucleus. The new combinations of these Holospora-like bacteria with ciliate hosts were discussed from biogeographical and ecological points of view. Host specificity of symbionts as a general paradigm was critically reviewed as well.

Variation in ribosomal and mitochondrial DNA sequences demonstrates the existence of intraspecific groups in Paramecium multimicronucleatum (Ciliophora, Oligohymenophorea).

This is the first phylogenetic study of the intraspecific variability within Paramecium multimicronucleatum with the application of two-loci analysis (ITS1-5.8S-ITS2-5'LSU rDNA and COI mtDNA) carried out on numerous strains originated from different continents. The species has been shown to have a complex structure of several sibling species within taxonomic species. Our analysis revealed the existence of 10 haplotypes for the rDNA fragment and 15 haplotypes for the COI fragment in the studied material. The mean distance for all of the studied P. multimicronucleatum sequence pairs was p=0.025/0.082 (rDNA/COI). Despite the greater variation of the COI fragment, the COI-derived tree topology is similar to the tree topology constructed on the basis of the rDNA fragment. P. multimicronucleatum strains are divided into three main clades. The tree based on COI fragment analysis presents a greater resolution of the studied P. multimicronucleatum strains. Our results indicate that the strains of P. multimicronucleatum that appear in different clades on the trees could belong to different syngens.

Species-Specific Duplication of Surface Antigen Genes in Paramecium.

Paramecium is a free-living ciliate that undergoes antigenic variation and still the functions of these variable surface antigen coats in this non-pathogenic ciliate remain elusive. Only a few surface antigen genes have been described, mainly in the two model species P. tetraurelia strain 51 and P. primaurelia strain 156. Given the lack of suitable sequence data to allow for phylogenetics and deeper sequence comparisons, we screened the genomes of six different Paramecium species for serotype genes and isolated 548 candidates. Our approach identified the subfamilies of the isogenes of individual serotypes that were mostly represented by intrachromosomal gene duplicates. These showed different duplication levels, and chromosome synteny suggested rather young duplication events after the emergence of the P. aurelia species complex, indicating a rapid evolution of surface antigen genes. We were able to identify the different subfamilies of the surface antigen genes with internal tandem repeats, which showed consensus motifs across species. The individual isogene families showed additional consensus motifs, indicating that the selection pressure holds individual amino acids constant in these repeats. This may be a hint of the receptor function of these antigens rather than a presentation of random epitopes, generating the variability of these surface molecules.

Cryptic Diversity in Paramecium multimicronucleatum Revealed with a Polyphasic Approach.

Paramecium (Ciliophora) systematics is well studied, and about twenty morphological species have been described. The morphological species may include several genetic species. However, molecular phylogenetic analyses revealed that the species diversity within Paramecium could be even higher and has raised a problem of cryptic species whose statuses remain uncertain. In the present study, we provide the morphological and molecular characterization of two novel Paramecium species. While Paramecium lynni n. sp., although morphologically similar to P. multimicronucleatum, is phylogenetically well separated from all other Paramecium species, Paramecium fokini n. sp. appears to be a cryptic sister species to P. multimicronucleatum. The latter two species can be distinguished only by molecular methods. The number and structure of micronuclei, traditionally utilized to discriminate species in Paramecium, vary not only between but also within each of the three studied species and, thus, cannot be considered a reliable feature for species identification. The geographic distribution of the P. multimicronucleatum and P. fokini n. sp. strains do not show defined patterns, still leaving space for a role of the geographic factor in initial speciation in Paramecium. Future findings of new Paramecium species can be predicted from the molecular data, while morphological characteristics appear to be unstable and overlapping at least in some species.

'Candidatus Gromoviella agglomerans', a novel intracellular Holosporaceae parasite of the ciliate Paramecium showing marked genome reduction.

Holosporales are an alphaproteobacterial lineage encompassing bacteria obligatorily associated with multiple diverse eukaryotes. For most representatives, little is known on the interactions with their hosts. In this study, we characterized a novel Holosporales symbiont of the ciliate Paramecium polycaryum. This bacterium inhabits the host cytoplasm, frequently forming quite large aggregates. Possibly due to such aggregates, host cells sometimes displayed lethal division defects. The symbiont was also able to experimentally stably infect another Paramecium polycaryum strain. The bacterium is phylogenetically related with symbionts of other ciliates and diplonemids, forming a putatively fast-evolving clade within the family Holosporaceae. Similarly to many close relatives, it presents a very small genome (<600 kbp), and, accordingly, a limited predicted metabolism, implying a heavy dependence on Paramecium, thanks also to some specialized membrane transporters. Characterized features, including the presence of specific secretion systems, are overall suggestive of a mild parasitic effect on the host. From an evolutionary perspective, a potential ancestral trend towards pronounced genome reduction and possibly linked to parasitism could be inferred, at least among fast-evolving Holosporaceae, with some lineage-specific traits. Interestingly, similar convergent features could be observed in other host-associated lineages, in particular Rickettsiales among Alphaproteobacteria.

Evolutionary Plasticity of Mating-Type Determination Mechanisms in Paramecium aurelia Sibling Species.

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.

Genetic diversity in the Paramecium aurelia species complex.

Current understanding of the population genetics of free-living unicellular eukaryotes is limited, and the amount of genetic variability in these organisms is still a matter of debate. We characterized-reproductively and genetically-worldwide samples of multiple Paramecium species belonging to a cryptic species complex, Paramecium aurelia, whose species have been shown to be reproductively isolated. We found that levels of genetic diversity both in the nucleus and in the mitochondrion are substantial within groups of reproductively compatible P. aurelia strains but drop considerably when strains are partitioned according to their phylogenetic groupings. Our study reveals the existence of discrepancies between the mating behavior of a number of P. aurelia strains and their multilocus genetic profile, a controversial finding that has major consequences for both the current methods of species assignment and the species problem in the P. aurelia complex.

Electrophoretic karyotype polymorphism of sibling species of the Paramecium aurelia complex.

Variability of karyotypes is one of the main mechanisms of speciation in organisms. Electrophoretic karyotypes of the macronucleus (MAC) obtained by pulsed-field gel electrophoresis were compared for 86 strains of all 15 sibling species of the Paramecium aurelia complex in order to determine if karyotype differences corresponded to biological species boundaries. Because the electrophoretic karyotype of the MAC reflects indirectly the frequency and distribution of fragmentation sites in the micronuclear (MIC) chromosomes, any change in MAC electrophoretic karyotype may be a marker of certain chromosomal mutations in the MIC. Thirteen main variants of electrophoretic MAC karyotypes were observed in this species complex. Ten of them appeared to correspond to biological species, while the three other variants characterized several species each. Intraspecific polymorphism was observed for several species: in some cases a certain variant of MAC karyotype was specific for all strains from the same part of the world. Distribution of the MAC karyotype variants along molecular phylogenetic trees of the P. aurelia complex shows that isolation of each species or group of species of this complex was accompanied by divergence in the molecular organization of the genome.

Species of the Paramecium aurelia complex in Russia: new stands and overall distribution.

New stands of Paramecium biaurelia, P. triaurelia, P. tetraurelia, P. pentaurelia, P. novaurelia, and P. dodecaurelia were recorded in Russia. Especially interesting is the record of P. novaurelia in Vladivostok, Russian Far East, as it is a very rare species outside of Europe. The distribution of species of the Paramecium aurelia complex in Eurasia with emphasis on findings in Russia is discussed.