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

First molecular evidence of hybridization in endosymbiotic ciliates (Protista, Ciliophora)

Hybridization is an important evolutionary process that can fuel diversification via formation of hybrid species or can lead to fusion of previously separated lineages by forming highly diverse species complexes. We provide here the first molecular evidence of hybridization in wild populations of ciliates, a highly diverse group of free-living and symbiotic eukaryotic microbes. The impact of hybridization was studied on the model of Plagiotoma, an obligate endosymbiont of the digestive tube of earthworms, using split decomposition analyses and species networks, 2D modeling of the nuclear rRNA molecules and compensatory base change analyses as well as multidimensional morphometrics. Gene flow slowed down and eventually hampered the diversification of Lumbricus-dwelling plagiotomids, which collapsed into a single highly variable biological entity, the P. lumbrici complex. Disruption of the species boundaries was suggested also by the continuum of morphological variability in the phenotypic space. On the other hand, hybridization conspicuously increased diversity in the nuclear rDNA cistron and somewhat weakened the host structural specificity of the P. lumbrici complex, whose members colonize a variety of phylogenetically closely related anecic and epigeic earthworms. By contrast, another recorded species, P. aporrectodeae sp. n., showed no signs of introgression, no variability in the rDNA cistron, and very high host specificity. These contrasting eco-evolutionary patterns indicate that hybridization might decrease the alpha-diversity by dissolving species boundaries, weaken the structural host specificity by broadening ecological amplitudes, and increase the nuclear rDNA variability by overcoming concerted evolution within the P. lumbrici species complex.

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.

Paramecium bursaria—A Complex of Five Cryptic Species: Mitochondrial DNA COI Haplotype Variation and Biogeographic Distribution

Ciliates are a diverse protistan group and many consist of cryptic species complexes whose members may be restricted to particular biogeographic locations. Mitochondrial genes, characterized by a high resolution for closely related species, were applied to identify new species and to distinguish closely related morphospecies. In the current study, we analyzed 132 sequences of COI mtDNA fragments obtained from P. bursaria species collected worldwide. The results allowed, for the first time, to generate a network of COI haplotypes and demonstrate the relationships between P. bursaria strains, as well as to confirm the existence of five reproductively isolated haplogroups. The P. bursaria haplogroups identified in the present study correspond to previously reported syngens (R1, R2, R3, R4, and R5), thus we decided to propose the following binominal names for each of them: P. primabursaria, P. bibursaria, P. tribursaria, P. tetrabursaria, and P. pentabursaria, respectively. The phylogeographic distribution of P. bursaria species showed that P. primabursaria and P. bibursaria were strictly Eurasian, except for two South Australian P. bibursaria strains. P. tribursaria was found mainly in Eastern Asia, in two stands in Europe and in North America. In turn, P. tetrabursaria was restricted to the USA territory, whereas P. pentabursaria was found in two European localities.

Genetic Diversity and Phylogeny of the Genus Euplotes (Protozoa, Ciliophora) Revealed by the Mitochondrial CO1 and Nuclear Ribosomal Genes

Nuclear ribosomal and mitochondrial genes have been utilized individually or in combination to identify known species and discriminate closely related species. However, compared with metazoans, genetic diversity within the ciliate order Euplotida is poorly known. The aim of this study is to investigate how much nucleotide sequence divergence occurs within Euplotes. A total of 14 new gene sequences, comprising four SSU rDNA and 10 CO1 (including three species for the first time) were obtained. Phylogenetic analyses were carried out based on sequences of two DNA fragments from the same 27 isolates. We found that CO1 revealed a larger interspecific divergence than the SSU rRNA gene, thus demonstrating a higher resolution for separating congeners. Genetic distances differ significantly at the species level. Euplotes balteatus was revealed to have a large intraspecific variation at two loci, while E. vannus showed different levels of haplotype variability, which appeared as a polyphyletic cluster on the CO1 tree. These high genetic divergences suggest the presence of more cryptic species. By contrast, the CO1 gene showed low variability within E. raikovi, appearing as monophyletic clusters, which indicates that this species could be identified based on this gene. Conclusively, CO1 is a suitable marker for the study of genetic diversity within Euplotes, and increased taxon sampling gives an opportunity to screen relationships among members of this genus. Additionally, current data present no clear biogeographical pattern for Euplotes.

Some morphologically distinguishable hypotrich ciliates share identical 18S rRNA gene sequences – taxonomic insights from a case study on Oxytricha species (Protista, Ciliophora)

Modern taxonomic studies of ciliated protozoa require both morphological and molecular data. One dilemma is how to distinguish morphologically similar species with few nucleotide differences in the widely used marker, the 18S rRNA gene. In the present study, two Oxytricha species were morphologically documented using light and electron microscopy. The mitochondrial cytochrome c oxidase subunit I (COI) gene and a fragment of the rRNA gene covering the 18S-ITS1-5.8S-ITS2-28S rRNA gene regions were sequenced. Phylogenetic analyses of all available Oxytricha granulifera-related populations were performed to reveal the internal relationships of this group. We described a new species, Oxytricha atypica sp. nov., distinguished from its congeners by having seven postoral ventral cirri resulting from the additional fragmentation of anlage V during ontogenesis. Although their 18S rRNA genes differ by only one nucleotide, divergence of the COI gene is as high as 11.8% between O. atypica and the closely related species, O. granulifera. All but one of the COI nucleotide substitutions were synonymous. We documented the highly conserved nature of the 18S rRNA gene in the morphospecies of Oxytricha. Based on these findings, we speculate that O. granulifera contains cryptic species or morphospecies needing further characterization, and new insights for the taxonomy of hypotrich ciliates are also discussed.

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.

Ciliates as model organisms for the ecotoxicological risk assessment of heavy metals: A meta-analysis

Ciliates are key components of aquatic ecosystems, significantly contributing to the decomposition of organic matter and energy transfer to higher trophic levels. They are considered good biological indicators of chemical pollution and relatively sensitive to heavy metal contamination. In this study, we performed a meta-analysis of the available toxicity data of heavy metals and ciliates to assess: (1) the sensitivity of freshwater ciliates to different heavy metals, (2) the relative sensitivity of ciliates in comparison to the standard test species used in ecotoxicological risk assessment, and (3) the difference in sensitivity across ciliate taxa. Our study shows that the tolerance of ciliates to heavy metals varies notably, which is partly influenced by differences in methodological conditions across studies. Ciliates are, in general, sensitive to Mercury > Cadmium > Copper > Zinc > Lead > Chromium. Also, this study shows that most ciliates are more tolerant to heavy metal pollution than the standard test species used in ecotoxicological risk assessments, i.e., Raphidocelis subcapitata, Daphnia magna, and Onchornyncus mykiss. Threshold concentrations derived from toxicity data for these species is expected to confer sufficient protection for the vast majority of ciliate species. Our data analysis also shows that the most commonly tested ciliate species, Paramecium caudatum and Tetrahymena thermophila, are not necessarily the most sensitive ones to heavy metal pollution. Finally, this study stresses the importance of developing standard toxicity test protocols for ciliates, which could lead to a better comprehension of the toxicological impact of heavy metals and other contaminants to ciliate species.

Paramecium Diversity and a New Member of the Paramecium aurelia Species Complex Described from Mexico

Paramecium (Ciliophora) is an ideal model organism to study the biogeography of protists. However, many regions of the world, such as Central America, are still neglected in understanding Paramecium diversity. We combined morphological and molecular approaches to identify paramecia isolated from more than 130 samples collected from different waterbodies in several states of Mexico. We found representatives of six Paramecium morphospecies, including the rare species Paramecium jenningsi, and Paramecium putrinum, which is the first report of this species in tropical regions. We also retrieved five species of the Paramecium aurelia complex, and describe one new member of the complex, Paramecium quindecaurelia n. sp., which appears to be a sister species of Paramecium biaurelia. We discuss criteria currently applied for differentiating between sibling species in Paramecium. Additionally, we detected diverse bacterial symbionts in some of the collected ciliates.

Documentation of a new hypotrich species in the family Amphisiellidae, Lamtostyla gui n. sp. (Protista, Ciliophora) using a multidisciplinary approach

An integrated approach considering both morphologic and molecular data is now required to improve biodiversity estimations and provide more robust systematics interpretations in hypotrichs, a highly differentiated group of ciliates. In present study, we document a new hypotrich species, Lamtostyla gui n. sp., collected from Chongming wetland, Shanghai, China, based on investigations using living observation, protargol staining, scanning and transmission electron microscopy, and gene sequencing. The new species is mainly recognized by having a short amphisiellid median cirral row composed of four cirri, three frontoventral cirri, three dorsal kinetids, four to eight macronuclear nodules, and small colorless cortical granules distributed as rosettes around dorsal bristles. Transmission electron microscope observation finds the associated microtubules of cirri and pharyngeal discs of L. gui are distinct from those in other hypotrichs. Morphogenesis of this species indicates that parental adoral membranelles retained intact or partial renewed is a potential feature to separate Lamtostyla granulifera-group and Lamtostyla lamottei-group. Phylogenetic analysis based on small subunit ribosomal RNA (rRNA) gene shows that this molecular marker is not useful to resolve phylogenetic relationships of the genus Lamtostyla, as well as many other hypotrichous taxa. We additionally characterize the internal transcribed spacers (ITS) region and the almost complete large subunit rRNA, which will be essential for future studies aimed at solving phylogenetic problems of Lamtostyla, or even the family Amphisiellidae. As a final remark, the critical screening of GenBank using ITS genes of our organism allows us to recognize a large amount of hypotrichous sequences have been misclassified as fungi. This observation suggests that hypotrichs could be frequently found in fungi-rich environment and overlooked by fungal specialists.

Ciliate Environmental Diversity Can Be Underestimated by the V4 Region of SSU rDNA: Insights from Species Delimitation and Multilocus Phylogeny of Pseudokeronopsis (Protist, Ciliophora)

Metabarcoding and high-throughput sequencing methods have greatly improved our understanding of protist diversity. Although the V4 region of small subunit ribosomal DNA (SSU-V4 rDNA) is the most widely used marker in DNA metabarcoding of eukaryotic microorganisms, doubts have recently been raised about its suitability. Here, using the widely distributed ciliate genus Pseudokeronopsis as an example, we assessed the potential of SSU-V4 rDNA and four other nuclear and mitochondrial markers for species delimitation and phylogenetic reconstruction. Our studies revealed that SSU-V4 rDNA is too conservative to distinguish species, and a threshold of 97% and 99% sequence similarity detected only one and three OTUs, respectively, from seven species. On the basis of the comparative analysis of the present and previously published data, we proposed the multilocus marker including the nuclear 5.8S rDNA combining the internal transcribed spacer regions (ITS1-5.8S-ITS2) and the hypervariable D2 region of large subunit rDNA (LSU-D2) as an ideal barcode rather than the mitochondrial cytochrome c oxidase subunit 1 gene, and the ITS1-5.8S-ITS2 as a candidate metabarcoding marker for ciliates. Furthermore, the compensating base change and tree-based criteria of ITS2 and LSU-D2 were useful in complementing the DNA barcoding and metabarcoding methods by giving second structure and phylogenetic evidence.

Population Genetics of Paramecium Mitochondrial Genomes: Recombination, Mutation Spectrum, and Efficacy of Selection

The evolution of mitochondrial genomes and their population-genetic environment among unicellular eukaryotes are understudied. Ciliate mitochondrial genomes exhibit a unique combination of characteristics, including a linear organization and the presence of multiple genes with no known function or detectable homologs in other eukaryotes. Here we study the variation of ciliate mitochondrial genomes both within and across 13 highly diverged Paramecium species, including multiple species from the P. aurelia species complex, with four outgroup species: P. caudatum, P. multimicronucleatum, and two strains that may represent novel related species. We observe extraordinary conservation of gene order and protein-coding content in Paramecium mitochondria across species. In contrast, significant differences are observed in tRNA content and copy number, which is highly conserved in species belonging to the P. aurelia complex but variable among and even within the other Paramecium species. There is an increase in GC content from ∼20% to ∼40% on the branch leading to the P. aurelia complex. Patterns of polymorphism in population-genomic data and mutation-accumulation experiments suggest that the increase in GC content is primarily due to changes in the mutation spectra in the P. aurelia species. Finally, we find no evidence of recombination in Paramecium mitochondria and find that the mitochondrial genome appears to experience either similar or stronger efficacy of purifying selection than the nucleus.

Intra-population genetic diversity and its effects on outlining genetic diversity of ciliate populations: Using Paramecium multimicronucleatum as an example

Questions regarding ciliate distribution (endemism vs. cosmopolitanism) and degree of genetic diversity (high vs. low) remain unsettled, even when the same organism is under investigation. Presence of genes with high copy number and amplification of non-dominant haplotypes might account for the observed discordance in these studies. Herein, we used direct PCR and cloning sequencing to examine intra-population sequence diversity and its effect on assessments of phylogeography of Paramecium multimicronucleatum. Totally, 381 ITS1-5.8S rDNA-ITS2-28S rDNA and 304 mitochondrial cytochrome oxidase subunit I (COI) gene sequences were generated for 18 populations of P. multimicronucleatum. The following results were obtained: (1) Direct sequencing of PCR products captured the dominant ITS and LSU haplotypes, indicating that it is an appropriate strategy for constructing phylogeography of large-scale spatial populations. (2) Deep cloning was deemed more appropriate for the COI gene for population level studies, as direct sequencing could not easily capture the dominant haplotypes. (3) No endemic populations of P. multinucleatum were noted, indicating origin from a single founder population. (4) Nuclear genetic diversity within temporal populations was high, but only the dominant haplotypes seemed to be passed on to subsequent generations.

Evaluation of the molecular variability and characteristics of Paramecium polycaryum and Paramecium nephridiatum, within subgenus Cypriostomum (Ciliophora, Protista)

Although some Paramecium species are suitable research objects in many areas of life sciences, the biodiversity structure of other species is almost unknown. In the current survey, we present a molecular analysis of 60 Cypriostomum strains, which for the first time allows for the study of intra- and interspecific relationships within that subgenus, as well as the assessment of the biogeography patterns of its morphospecies. Analysis of COI mtDNA variation revealed three main clades (separated from each other by approximately 130 nucleotide substitutions), each one with internal sub-clusters (differing by 30 to 70 substitutions - a similar range found between P. aurelia cryptic species and P. bursaria syngens). The first clade is represented exclusively by P. polycaryum; the second one includes only four strains identified as P. calkinsi. The third cluster seems to be paraphyletic, as it includes P. nephridiatum, P. woodruffi, and Eucandidatus P. hungarianum. Some strains, previously identified as P. calkinsi, had COI sequences identical or very similar to P. nephridiatum ones. Morphological reinvestigation of several such strains revealed common morphological features with P. nephridiatum. The paper contains new information concerning speciation within particular species, i.e. existence of cryptic species within P. polycaryum (three) and in P. nephridiatum (six).

Utility of mitochondrial CO1 sequences for species discrimination of Spirotrichea ciliates (Protozoa, Ciliophora)

Ciliates are a diverse species group of the Protozoa, and nuclear and mitochondrial genes have been utilized to discover new species and discriminate closely related species. The mitochondrial cytochrome c oxidase subunit 1 (CO1) gene has been used to discriminate metazoan species and has also been applied for some groups in the phylum Ciliophora. However, it is difficult to produce a universal primer as a standard barcode, because unlike metazoans, mitochondrial DNA sequences of ciliates are long and highly variable. Therefore, to design the new primer set, we sequenced the mitochondrial genomes of two pseudokeronopsids in the class Spirotrichea using next-generation sequencing technology (HiSeq™ 2000). Based on putative CO1 gene fragments of the pseudokeronopsids, we designed the new primer set and successfully sequenced the CO1 of 69 populations representing 47 species (five orders, 14 families, and 27 genera). We found that CO1 showed higher resolution for separating congeneric species than did nuclear SSU rRNA gene sequences, and we identified some putative cryptic species.

Delimiting Species Boundaries within a Paraphyletic Species Complex: Insights from Morphological, Genetic, and Molecular Data on Paramecium sonneborni (Paramecium aurelia species complex, Ciliophora, Protozoa)

The demarcation of boundaries between protist species is often problematic because of the absence of a uniform species definition, the abundance of cryptic diversity, and the occurrence of convergent morphology. The ciliates belonging to the Paramecium aurelia complex, consisting of 15 species, are a good model for such systematic and evolutionary studies. One member of the complex is P. sonneborni, previously known only from one stand in Texas (USA), but recently found in two new sampling sites in Cyprus (creeks running to Salt Lake and Oroklini Lake near Larnaca). The studied Paramecium sonneborni strains (from the USA and Cyprus) reveal low viability in the F1 and F2 generations of interstrain hybrids and may be an example of ongoing allopatric speciation. Despite its molecular distinctiveness, we postulate that P. sonneborni should remain in the P. aurelia complex, making it a paraphyletic taxon. Morphological studies have revealed that some features of the nuclear apparatus of P. sonneborni correspond to the P. aurelia spp. complex, while others are similar to P. jenningsi and P. schewiakoffi. The observed discordance indicates rapid splitting of the P. aurelia-P. jenningsi-P. schewiakoffi group, in which genetic, morphological, and molecular boundaries between species are not congruent.

Genetic differentiation of the mitochondrial cytochrome oxidase C subunit I gene in genus Paramecium (Protista, Ciliophora)

Background

The mitochondrial cytochrome c oxidase subunit I (COI) gene is being used increasingly for evaluating inter- and intra-specific genetic diversity of ciliated protists. However, very few studies focus on assessing genetic divergence of the COI gene within individuals and how its presence might affect species identification and population structure analyses.

Methodology/Principal findings

We evaluated the genetic variation of the COI gene in five Paramecium species for a total of 147 clones derived from 21 individuals and 7 populations. We identified a total of 90 haplotypes with several individuals carrying more than one haplotype. Parsimony network and phylogenetic tree analyses revealed that intra-individual diversity had no effect in species identification and only a minor effect on population structure.

Conclusions

Our results suggest that the COI gene is a suitable marker for resolving inter- and intra-specific relationships of Paramecium spp.

Multiple lines of evidence shed light on the occurrence of paramecium (ciliophora, oligohymenophorea) in bromeliad tank water

Phytotelmata are vegetal structures that hold water from the rain, and organic matter from the forest and the soil, resulting in small, compartmentalized bodies of water, which provide an essential environment for the establishment and development of many organisms. These microenvironments generally harbor endemic species, but many organisms that are found in lakes and rivers, are also present. Here, we report, for the first time, the occurrence of the ciliate genus Paramecium in the tank of the bromeliad species Aechmaea distichantha. The identification of the Paramecium species was performed based on live observations, protargol impregnation, scanning electronic microscopy, and sequencing of the 18s rRNA. The absence of Paramecium from bromeliad tank water was highlighted in several earlier investigations, and may be due to the fact that this species is unable to make cysts. The occurrence of Paramecium multimicronucleatum in our samples may be explained by the proximity between the bromeliads and the river, a potential source of the species. Further, we also believe that the counting methodology used in our study provides a more accurate analysis of the species diversity, since we investigated all samples within a maximum period of 6 h after sampling, allowing minimum loss of specimens.