Species of Tetrahymena Identical by Small Subunit rRNA Gene Sequences are Discriminated by Mitochondrial Cytochrome c Oxidase I Gene Sequences

Integrative studies on three new freshwater Amphileptus species (Ciliophora, Pleurostomatida) discovered in northern China

The morphology and molecular phylogeny of freshwater pleurostomatid ciliates are insufficiently explored. In the present study, we investigated three new Amphileptus species discovered in Lake Weishan and its vicinity, northern China, using standard alpha-taxonomic methods. Amphileptus paracarchesii sp. nov. is characterized by a lateral fossa (groove) in the posterior body portion, four macronuclear nodules, contractile vacuoles distributed along the dorsal margin, and 4-6 left and 44-50 right somatic kineties. Amphileptus pilosus sp. nov. differs from congeners by having 4-14 macronuclear nodules, numerous contractile vacuoles scattered throughout the cytoplasm, and 22-31 left and 35-42 right somatic kineties. Amphileptus orientalis sp. nov. is characterized by two ellipsoidal macronuclear nodules, three ventral contractile vacuoles, and about four left and 31-35 right somatic kineties. Phylogenetic analyses of nuclear small subunit ribosomal DNA (SSU rDNA) sequences indicate that the family Amphileptidae might be monophyletic while the genus Amphileptus is paraphyletic, as Pseudoamphileptus macrostoma robustly groups with Amphileptus sp. Although deep phylogenetic relationships of amphileptids are poorly resolved, multiple well-delimited species groups are recognizable within the genus Amphileptus.

Supplementary Information

The online version contains supplementary material available at 10.1007/s42995-022-00143-0.

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.

DNA barcoding and coalescent-based delimitation of endosymbiotic clevelandellid ciliates (Ciliophora: Clevelandellida): a shift to molecular taxonomy in the inventory of ciliate diversity in panesthiine cockroaches

Phylogenetically distinct lineages may be hidden behind identical or highly similar morphologies. The phenomenon of morphological crypticity has been recently detected in symbiotic ciliates of the family Clevelandellidae, as multivariate and Fourier shape analyses failed to distinguish genetically distinct taxa. To address the question of species boundaries, the phylogenetic information contained in the rDNA cistron of clevelandellid ciliates, which had been isolated from the digestive tract of blaberid cockroaches, was studied using a multifaceted statistical approach. Multigene phylogenies revealed that the genus Clevelandella is paraphyletic containing members of the genus Paraclevelandia. To resolve the paraphyly of Clevelandella, two new genera, Anteclevelandella gen. nov. and Rhynchoclevelandella gen. nov., are proposed based on morphological synapomorphies and shared molecular characters. Multigene analyses and Bayesian species delimitation supported the existence of 13 distinct species within the family Clevelandellidae, eight of which represent new taxa. Moreover, two new Nyctotherus species were recognized within the clade that is sister to the Clevelandellidae. According to the present distance and network analyses, the first two domains of the 28S rRNA gene showed much higher power for species discrimination than the 18S rRNA gene and ITS region. Therefore, the former molecular marker was proposed to be a suitable group-specific barcode for the family Clevelandellidae.

A discovery of two new Tetrahymena species parasitizing slugs and mussels: morphology and multi-gene phylogeny of T. foissneri sp. n. and T. unionis sp. n

The presence of parasitic ciliates of the hymenostome genus Tetrahymena was examined in 150 mollusks belonging to six bivalve and 13 gastropod species in Slovakia, Central Europe. Tetrahymenids were detected only in two species, viz., in the invasive Lusitanian slug (Arion vulgaris) and in the native swollen river mussel (Unio tumidus). Although only 10.52% of the examined mollusk taxa were positive, their Tetrahymena infections were very intensive accounting for several hundreds of ciliates per host. Phylogenetic analyses of the 16S and 18S rRNA genes as well as of the barcoding region of the gene encoding for cytochrome c oxidase subunit I revealed that both isolates represent new taxa, T. foissneri sp. n. and T. unionis sp. n. The former species belongs to the 'borealis' clade and its nearest relative is T. limacis, a well-known parasite of slugs and snails. Besides molecular data, T. foissneri can be distinguished from T. limacis also morphologically by the body shape of the parasitic-phase form, dimensions of micronuclei, and the silverline system. On the other hand, T. unionis was classified within the 'paravorax' clade along with T. pennsylvaniensis, T. glochidiophila, and T. nigricans. Although these four species are genetically distinct, T. unionis could be morphologically separated only from T. nigricans by body shape and size. The present study suggests that both aquatic and terrestrial mollusks represent interesting hosts for the discovery of novel Tetrahymena lineages.

Cryptic host-driven speciation of mobilid ciliates epibiotic on freshwater planarians

Mobilids are among the most taxonomically diverse but morphologically uniform groups of epibiotic ciliates. They attach to their hosts by means of an adhesive disc as harmless commensals such as Urceolaria, or as parasites causing significant economic loss such as some Trichodina species. We investigated the diversity, species boundaries, and phylogenetic relationships of mobilids associated with freshwater planarians, using 114 new sequences of two mitochondrial (16S rRNA gene and cytochrome c oxidase gene) and five nuclear (18S rRNA gene, ITS1-5.8S-ITS2 region, D1/D2 domains of 28S rRNA gene) markers. Although the morphological disparity of the isolated trichodinids and urceolariids was low, Bayesian coalescent analyses revealed the existence of five distinct evolutionary lineages/species given the seven molecular markers. The occurrence of mobilids perfectly correlated with their planarian hosts: Trichodina steinii and two Urceolaria mitra-like taxa were associated exclusively with the planarian Dugesia gonocephala, Trichodina polycelis sp. n. with the planarian Polycelis felina, and Trichodina schmidtea sp. n. with the planarian Schmidtea polychroa. Host organisms thus very likely constitute sharply isolated niches that might permit speciation of their epibiotic ciliates, even though no distinct morphological features appear to be recognizable among ciliates originating from different hosts.

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.

Re-discovery and novel contributions to morphology and multigene phylogeny of Myxophyllum steenstrupi (Ciliophora: Pleuronematida), an obligate symbiont of terrestrial pulmonates

Myxophyllum steenstrupi is a symbiotic ciliate living in the body slime and mantle cavity of terrestrial pulmonates (Gastropoda: Pulmonata). In the present study, M. steenstrupi was re-discovered after almost 30 years and characterized using an integrative morpho-molecular approach for the first time. Myxophyllum is distinguished by a broadly ovate, about 140 × 115 μm-sized body, a nuclear apparatus typically composed of seven macronuclear nodules and a single micronucleus, a central contractile vacuole, a shallow oral cavity situated in the posterior body region and dense somatic ciliature with extensive thigmotactic field. According to the present phylogenetic analyses of two mitochondrial and five nuclear markers, M. steenstrupi is classified in the predominantly free-living order Pleuronematida (Oligohymenophorea: Scuticociliatia). This order also encompasses other taxa isolated from molluscs and traditionally classified along with Myxophyllum in the order Thigmotrichida. The proper classifications of Myxophyllum was hampered by the dramatic remodelling of its oral apparatus (reduction of the paroral membrane and adoral organelles, formation of vestibular kineties), a transformation that was likely induced by its firm association with terrestrial gastropods. The present study also documents that various ciliate lineages independently became commensals or parasites of various aquatic and terrestrial molluscs.

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.

Multi-gene phylogeny of Tetrahymena refreshed with three new histophagous species invading freshwater planarians

Planarians represent an insufficiently explored group of aquatic invertebrates that might serve as hosts of histophagous ciliates belonging to the hymenostome genus Tetrahymena. During our extensive research on freshwater planarians, parasitic tetrahymenas were detected in two of the eight planarian species investigated, namely, in Dugesia gonocephala and Girardia tigrina. Using the 16S and 18S rRNA genes as well as the barcoding cytochrome oxidase subunit I, one ciliate species was identified as T. scolopax and three species were recognized as new forms: T. acanthophora, T. dugesiae, and T. nigricans. Thus, 25% of the examined planarian taxa are positive for Tetrahymena species and three of them represent new taxa, indicating a large undescribed ciliate diversity in freshwater planarians. According to phylogenetic analyses, histophagous tetrahymenas show a low phylogenetic host specificity. Although T. acanthophora, T. dugesiae, and T. scolopax clustered together within the "borealis" clade, the former species has been detected exclusively in G. tigrina, while the two latter species only in D. gonocephala. Tetrahymena nigricans, which has been isolated only from G. tigrina, was classified within the "paravorax" clade along with T. glochidiophila which feeds on glochidia. The present phylogenetic reconstruction of ancestral life strategies suggested that the last common ancestor of the family Tetrahymenidae was free-living, unlike the progenitor of the subclass Hymenostomatia which was very likely parasitic. Consequently, there were at least seven independent shifts back to parasitism/histophagy within Tetrahymena: one each in the "paravorax" and "australis" clades and at least five transfers back to parasitism in the "borealis" clade.

Deciphering phylogenetic relationships and delimiting species boundaries using a Bayesian coalescent approach in protists: A case study of the ciliate genus Spirostomum (Ciliophora, Heterotrichea)

The ciliate genus Spirostomum comprises eight morphospecies, inhabiting diverse aquatic environments worldwide, where they can be used as water quality indicators. Although Spirostomum species are relatively easily identified using morphological methods, the previous nuclear rDNA-based phylogenies indicated several conflicts in morphospecies delineation. Moreover, the single locus phylogenies and previous analytical approaches could not unambiguously resolve phylogenetic relationships among Spirostomum morphospecies. Here, we attempt to investigate species boundaries and evolutionary history of Spirostomum taxa, using 166 new sequences from multiple populations employing one mitochondrial locus (CO1 gene) and two nuclear loci (rRNA operon and alpha-tubulin gene). In accordance with previous studies, relationships among the eight Spirostomum morphospecies were poorly supported statistically in individual gene trees. To overcome this problem, we utilised for the first time in ciliates the Bayesian coalescent approach, which accounts for ancestral polymorphisms, incomplete lineage sorting, and recombination. This strategy enabled us to robustly resolve deep relationships between Spirostomum species and to support the hypothesis that taxa with compact macronucleus and taxa with moniliform macronucleus each form a distinct lineage. Bayesian coalescent-based delimitation analyses strongly statistically supported the traditional morphospecies concept but also indicated that there are two S. minus-like cryptic species and S. teres is non-monophyletic. Spirostomum teres was very likely defined by a set of ancestral features of lineages that also gave rise to S. yagiui and S. dharwarensis. However, molecular data from type populations of the morphospecies S. minus and S. teres are required to unambiguously resolve the taxonomic problems.

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.

Further analyses on the phylogeny of the subclass Scuticociliatia (Protozoa, Ciliophora) based on both nuclear and mitochondrial data

So far, the phylogenetic studies on ciliated protists have mainly based on single locus, the nuclear ribosomal DNA (rDNA). In order to avoid the limitations of single gene/genome trees and to add more data to systematic analyses, information from mitochondrial DNA sequence has been increasingly used in different lineages of ciliates. The systematic relationships in the subclass Scuticociliatia are extremely confused and largely unresolved based on nuclear genes. In the present study, we have characterized 72 new sequences, including 40 mitochondrial cytochrome oxidase c subunit I (COI) sequences, 29 mitochondrial small subunit ribosomal DNA (mtSSU-rDNA) sequences and three nuclear small subunit ribosomal DNA (nSSU-rDNA) sequences from 47 isolates of 44 morphospecies. Phylogenetic analyses based on single gene as well as concatenated data were performed and revealed: (1) compared to mtSSU-rDNA, COI gene reveals more consistent relationships with those of nSSU-rDNA; (2) the secondary structures of mtSSU-rRNA V4 region are predicted and compared in scuticociliates, which can contribute to discrimination of closely related species; (3) neither nuclear nor mitochondrial data support the monophyly of the order Loxocephalida, which may represent some divergent and intermediate lineages between the subclass Scuticociliatia and Hymenostomatia; (4) the assignments of thigmotrichids to the order Pleuronematida and the confused taxon Sulcigera comosa to the genus Histiobalantium are confirmed by mitochondrial genes; (5) both nuclear and mitochondrial data reveal that the species in the family Peniculistomatidae always group in the genus Pleuronema, suggesting that peniculistomatids are more likely evolved from Pleuronema-like ancestors; (6) mitochondrial genes support the monophyly of the order Philasterida, but the relationships among families of the order Philasterida remain controversial due to the discrepancies between their morphological and molecular data.

A study on the non-monophyletic genera Australothrix and Holostichides based on multigene and morphological analyses with a reexamination of type materials (Protozoa: Ciliophora)

In a study of ciliate diversity, we collected two new species of the genera Australothrix and Holostichides. Based on nuclear ribosomal genes, our study shows that these two genera are genetically non-monophyletic. To clarify the issues of the non-monophyly, we reexamined type material of H. heterotypicus, H. terrae, and Birojimia soyaensis. Based on multigene and morphometric analyses, Australothrix and Holostichides are clearly non-monophyletic even in mitochondrial CO1 gene trees. The multigene analyses show a clade composed of A. lineae sp. nov., H. terrae, and B. soyaensis, suggesting that a cytopharynx with argyrophilic structures might be their synapomorphy. A list of species with this type of cytopharynx is provided for further studies to either accept or reject this hypothesis. Based on the reexamination, we discriminate the morphologically nearly identical (cryptic/sibling) species H. obliquocirratus sp. nov. from H. heterotypicus and they show significant genetic dissimilarities in the multigene trees. They showed only few morphological (non-quantitative) differences and thus distinguishing them morphologically needs careful investigation.

Hidden genomic evolution in a morphospecies-The landscape of rapidly evolving genes in Tetrahymena

A morphospecies is defined as a taxonomic species based wholly on morphology, but often morphospecies consist of clusters of cryptic species that can be identified genetically or molecularly. The nature of the evolutionary novelty that accompanies speciation in a morphospecies is an intriguing question. Morphospecies are particularly common among ciliates, a group of unicellular eukaryotes that separates 2 kinds of nuclei—the silenced germline nucleus (micronucleus [MIC]) and the actively expressed somatic nucleus (macronucleus [MAC])—within a common cytoplasm. Because of their very similar morphologies, members of the Tetrahymena genus are considered a morphospecies. We explored the hidden genomic evolution within this genus by performing a comprehensive comparative analysis of the somatic genomes of 10 species and the germline genomes of 2 species of Tetrahymena. These species show high genetic divergence; phylogenomic analysis suggests that the genus originated about 300 million years ago (Mya). Seven universal protein domains are preferentially included among the species-specific (i.e., the youngest) Tetrahymena genes. In particular, leucine-rich repeat (LRR) genes make the largest contribution to the high level of genome divergence of the 10 species. LRR genes can be sorted into 3 different age groups. Parallel evolutionary trajectories have independently occurred among LRR genes in the different Tetrahymena species. Thousands of young LRR genes contain tandem arrays of exactly 90-bp exons. The introns separating these exons show a unique, extreme phase 2 bias, suggesting a clonal origin and successive expansions of 90-bp–exon LRR genes. Identifying LRR gene age groups allowed us to document a Tetrahymena intron length cycle. The youngest 90-bp exon LRR genes in T. thermophila are concentrated in pericentromeric and subtelomeric regions of the 5 micronuclear chromosomes, suggesting that these regions act as genome innovation centers. Copies of a Tetrahymena Long interspersed element (LINE)-like retrotransposon are very frequently found physically adjacent to 90-bp exon/intron repeat units of the youngest LRR genes. We propose that Tetrahymena species have used a massive exon-shuffling mechanism, involving unequal crossing over possibly in concert with retrotransposition, to create the unique 90-bp exon array LRR genes.

Genomic comparison of ten morphologically very similar species of ciliate from the genus Tetrahymena reveals how parallel microevolutionary processes have shaped their genomes and created unique genes through retrotransposition.

Techniques and tools for species identification in ciliates: a review

Ciliates are highly divergent unicellular eukaryotic organisms with nuclear dualism and a highly specialized ciliary pattern. They inhabit all biotopes and play crucial roles in regulating microbial food webs as they prey on bacteria, protists and even on microscopic animals. Nevertheless, subtle morphological differences and tiny sizes hinder proper species identification for many ciliates. In the present review, an attempt has been made to elaborate the various approaches used by modern day ciliate taxonomists for species identification. The different approaches involved in taxonomic characterization of ciliates such as classical (using live-cell observations, staining techniques, etc.), molecular (involving various marker genes) and statistical (delimitation of cryptic species) methods have been reviewed. Ecological and behavioural aspects in species identification have also been discussed. In present-day taxonomy, it is important to use a ‘total evidence’ approach in identifying ciliates, relying on both classical and molecular information whenever possible. This integrative approach will help in the mergence of classical methods with modern-day tools for comprehensive species description in future.

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.

Assessing SSU rDNA Barcodes in Foraminifera: A Case Study using Bolivina quadrata

The Small Subunit Ribosomal RNA gene (SSU rDNA) is a widely used tool to reconstruct phylogenetic relationships among foraminiferal species. Recently, the highly variable regions of this gene have been proposed as DNA barcodes to identify foraminiferal species. However, the resolution of these barcodes has not been well established, yet. In this study, we evaluate four SSU rDNA hypervariable regions (37/f, 41/f, 43/e, and 45/e) as DNA barcodes to distinguish among species of the genus Bolivina, with particular emphasis on Bolivina quadrata for which ten new sequences (KY468817-KY468826) were obtained during this study. Our analyses show that a single SSU rDNA hypervariable sequence is insufficient to resolve all Bolivina species and that some regions (37/f and 41/f) are more useful than others (43/e and 45/e) to distinguish among closely related species. In addition, polymorphism analyses reveal a high degree of variability. In the context of barcoding studies, these results emphasize the need to assess the range of intraspecific variability of DNA barcodes prior to their application to identify foraminiferal species in environmental samples; our results also highlight the possibility that a longer SSU rDNA region might be required to distinguish among species belonging to the same taxonomic group (i.e. genus).

The Green Tetrahymena utriculariae n. sp. (Ciliophora, Oligohymenophorea) with Its Endosymbiotic Algae (Micractinium sp.), Living in Traps of a Carnivorous Aquatic Plant

The genus Tetrahymena (Ciliophora, Oligohymenophorea) probably represents the best studied ciliate genus. At present, more than forty species have been described. All are colorless, i.e. they do not harbor symbiotic algae, and as aerobes they need at least microaerobic habitats. Here, we present the morphological and molecular description of the first green representative, Tetrahymena utriculariae n. sp., living in symbiosis with endosymbiotic algae identified as Micractinium sp. (Chlorophyta). The full life cycle of the ciliate species is documented, including trophonts and theronts, conjugating cells, resting cysts and dividers. This species has been discovered in an exotic habitat, namely in traps of the carnivorous aquatic plant Utricularia reflexa (originating from Okavango Delta, Botswana). Green ciliates live as commensals of the plant in this anoxic habitat. Ciliates are bacterivorous, however, symbiosis with algae is needed to satisfy cell metabolism but also to gain oxygen from symbionts. When ciliates are cultivated outside their natural habitat under aerobic conditions and fed with saturating bacterial food, they gradually become aposymbiotic. Based on phylogenetic analyses of 18S rRNA and mitochondrial cox1 genes T. utriculariae forms a sister group to Tetrahymena thermophila.

Uncovering Cryptic Diversity in Two Amoebozoan Species Using Complete Mitochondrial Genome Sequences

The Amoebozoa are a major eukaryotic lineage that encompasses a wide range of amoeboid organisms. The group is understudied from a systematic perspective: molecular tools have only been applied in the last 15 yr. Hence, there is an undersampling of both genes and taxa in the group especially compared to plants, animals, and fungi. Here, we present the complete mitochondrial genomes of two ubiquitous and abundant morpho-species (Acanthamoeba castellanii and Vermamoeba vermiformis). Both have mitochondrial genomes of close relatives previously available, enabling insights into recent divergences at a genomic scale, while simultaneously offering comparisons with divergence estimates obtained from traditionally used single genes, SSU rDNA and cox1. The newly sequenced mt genomes are significantly divergent from their previously sequenced conspecifics (A. castellannii 16.4% divergence at nucleotide level and 10.4% amino acid; V. vermiformis 21.6% and 13.1%, respectively), while divergence at the small subunit ribosomal DNA is below 1% within both species. Morphological analyses determined that these lineages are indistinguishable from their previously sequenced counterparts. Phylogenetic reconstructions using 26 mt genes also indicate a level of divergence that is comparable to divergence among species, while reconstructions using the small subunit ribosomal DNA (SSU rDNA) do not. In addition, we demonstrate that between closely related taxa, there are high levels of synteny, which can be explored for primer design to obtain larger fragments than the traditional barcoding genes. We conclude that, although most systematic work has relied on SSU, this gene alone can severely underestimate diversity. Thus, we suggest that the mt genome emerges as an alternative for unraveling the lower level phylogenetic relationships of Amoebozoa.

A morphogenetic survey on ciliate plankton from a mountain lake pinpoints the necessity of lineage-specific barcode markers in microbial ecology

Analyses of high-throughput environmental sequencing data have become the 'gold-standard' to address fundamental questions of microbial diversity, ecology and biogeography. Findings that emerged from sequencing are, e.g. the discovery of the extensive 'rare microbial biosphere' and its potential function as a seed-bank. Even though applied since several years, results from high-throughput environmental sequencing have hardly been validated. We assessed how well pyrosequenced amplicons [the hypervariable eukaryotic V4 region of the small subunit ribosomal RNA (SSU rRNA) gene] reflected morphotype ciliate plankton. Moreover, we assessed if amplicon sequencing had the potential to detect the annual ciliate plankton stock. In both cases, we identified significant quantitative and qualitative differences. Our study makes evident that taxon abundance distributions inferred from amplicon data are highly biased and do not mirror actual morphotype abundances at all. Potential reasons included cell losses after fixation, cryptic morphotypes, resting stages, insufficient sequence data availability of morphologically described species and the unsatisfying resolution of the V4 SSU rRNA fragment for accurate taxonomic assignments. The latter two underline the necessity of barcoding initiatives for eukaryotic microbes to better and fully exploit environmental amplicon data sets, which then will also allow studying the potential of seed-bank taxa as a buffer for environmental changes.

A new tetrahymena (ciliophora, oligohymenophorea) from groundwater of cape town, South Africa

The identification of species within the genus Tetrahymena is known to be difficult due to their essentially identical morphology, the occurrence of cryptic and sibling species and the phenotypic plasticity associated with the polymorphic life cycle of some species. We have combined morphology and molecular biology to describe Tetrahymena aquasubterranea n. sp. from groundwater of Cape Town, Republic of South Africa. The phylogenetic analysis compares the cox1 gene sequence of T. aquasubterranea with the cox1 gene sequences of other Tetrahymena species and uses the interior-branch test to improve the resolution of the evolutionary relationships. This showed a considerable genetic divergence of T. aquasubterranea to its next relative, T. farlyi, of 9.2% (the average cox1 divergence among bona fide species of Tetrahymena is ~ 10%). Moreover, the analysis also suggested a sister relationship between T. aquasubterranea and a big clade comprising T. farleyi, T. tropicalis, T. furgasoni and T. mobilis. The morphological data available for these species show that they share with T. aquasubterranea a pyriformis-like life style and at least two of them, T. farleyi and T. mobilis, a similar type II silverline pattern consisting of primary and secondary meridians. Tetrahymena aquasubterranea exhibits a biphasic life cycle with trophonts and theronts, is amicronucleate, and feeds on bacteria.

Development of a DNA barcoding system for the Ixodida (Acari: Ixodida)

To control the spread of tick-borne diseases, there is an urgent need to develop a reliable technique that can distinguish different species of ticks. DNA barcoding has been proved to be a powerful tool to identify species of arthropods, but this technique has not yet been developed for identifying ticks. Here, we screened and analyzed 1082 sequences of ticks from BOLD system and GenBank, consisting of 647 16S, 325 COI, and 110 18S. These sequences are reported in previous studies and considered to be correctly identified at the species level. Through the analyses of genetic divergences and neighbor-joining (NJ) phylogenetic relationships between the species of ticks, our results show that COI and 16S are reliable in discriminating species of ticks and the 18S could discriminate ticks at the genera level. New universal primers for 16S, 18S, and COI of ticks were designed and a DNA barcoding system for the Ixodida was developed. To assess the performance of this system, 57 specimens of ticks were collected within China. Our results show that DNA barcoding system could correctly identify the species of specimens in adult and subadult stages. This system would assist non-taxonomists to conveniently identify the species of Ixodida based on DNA sequences rather than morphological traits. However, there are still serious deficiencies in the information of 16S and COI of some species of ticks, and additional research is needed to resolve this problem.

Restricted distribution and limited gene flow in the model ciliate Tetrahymena thermophila

The biogeography of microbial eukaryotes has long been debated, but few phylogeographic data have been available to assess whether protists tend to have ubiquitous or endemic distributions. We addressed this issue in the ciliate Tetrahymena thermophila, a highly successful model system in cell and molecular biology. We found that this species has a distribution that is restricted to the Eastern United States, with high diversity in the northeast and low diversity across the rest of its distribution. We find high levels of population subdivision, low rates of migration and significant isolation by distance, supporting the moderate endemicity model of protist biogeography. This restricted gene flow may be a result of small population size, which would reduce the probability of migration events, or the inability to establish after migration. This work lays the foundation for T. thermophila to become a valuable model system for studying population biology.

Identification of Paramecium bursaria syngens through molecular markers--comparative analysis of three loci in the nuclear and mitochondrial DNA

This is the first attempt to resolve the phylogenetic relationship between different syngens of Paramecium bursaria and to investigate at a molecular level the intraspecific differentiation of strains originating from very distant geographical locations. Herein we introduce a new collection of five P. bursaria syngens maintained at St Petersburg State University, as the international collection of syngens was lost in the 1960s. To analyze the degree of speciation within Paramecium bursaria, we examined 26 strains belonging to five different syngens from distant and geographically isolated localities using rDNA (ITS1-5.8S-ITS2-5'LSU) fragments, mitochondrial cytochrome c oxidase subunit I (COI), and H4 gene fragments. It was shown that P. bursaria strains of the same syngens cluster together in all three inferred molecular phylogenies. The genetic diversity among the studied P. bursaria strains based on rDNA sequences was rather low. The COI divergence of Paramecium bursaria was also definitely lower than that observed in the Paramecium aurelia complex. The nucleotide sequences of the H4 gene analyzed in the present study indicate the extent of genetic differences between the syngens of Paramecium bursaria. Our study demonstrates the diagnostic value of molecular markers, which are important tools in the identification of Paramecium bursaria syngens.

A fuzzy-set-theory-based approach to analyse species membership in DNA barcoding

Reliable assignment of an unknown query sequence to its correct species remains a methodological problem for the growing field of DNA barcoding. While great advances have been achieved recently, species identification from barcodes can still be unreliable if the relevant biodiversity has been insufficiently sampled. We here propose a new notion of species membership for DNA barcoding-fuzzy membership, based on fuzzy set theory-and illustrate its successful application to four real data sets (bats, fishes, butterflies and flies) with more than 5000 random simulations. Two of the data sets comprise especially dense species/population-level samples. In comparison with current DNA barcoding methods, the newly proposed minimum distance (MD) plus fuzzy set approach, and another computationally simple method, 'best close match', outperform two computationally sophisticated Bayesian and BootstrapNJ methods. The new method proposed here has great power in reducing false-positive species identification compared with other methods when conspecifics of the query are absent from the reference database.

Small subunit ribosomal RNA and mitochondrial cytochrome c oxidase subunit 1 gene sequences of 21 strains of the parasitic scuticociliate Miamiensis avidus (Ciliophora, Scuticociliatia)

The scuticociliate Miamiensis avidus is a histophagous parasite that causes high mortality in cultured marine fishes. Small subunit ribosomal RNA (SSU rRNA) and mitochondrial cytochrome c oxidase subunit 1 (cox1) genes were analyzed for 21 strains of M. avidus isolated from diseased olive flounder (Paralichthys olivaceus), ridged-eye flounder (Pleuronichthys cornutus), and spotted knifejaw (Oplegnathus fasciatus) in Korea and Japan (collected in 2003-2007). Analysis of SSU rRNA gene sequences (1,759 bp) indicates they are very conserved with less than 0.17% (3 nucleotides) differences suggesting that SSU rRNA are useful to identify M. avidus; however, the cox1 gene (900 bp) has higher variations with intraspecific divergences up to 5.67% (51 nucleotides). A distance tree of cox1 gene sequences based on a neighbor-joining analysis can separate 21 strains into five cox1 types (two heterogeneous clusters and three individual branches). The cox1-type matches with serotype of strains but do not reflect geographical origins, host species, or pathogenicity.

Evidence for cryptic speciation in Carchesium polypinum Linnaeus, 1758 (Ciliophora: Peritrichia) inferred from mitochondrial, nuclear, and morphological markers

Protist diversity is currently a much debated issue in eukaryotic microbiology. Recent evidence suggests that morphological and genetic diversity might be decoupled in some groups of protists, including ciliates, and that these organisms might be much more diverse than their morphology implies. We sought to assess the genetic and morphological diversity of Carchesium polypinum, a widely distributed peritrich ciliate. The mitochondrial marker cytochrome c oxidase subunit I and the nuclear small subunit ribosomal RNA were used to examine genetic diversity. For the morphological assessment, live microscopy and Protargol staining were used. The mitochondrial marker revealed six robust, deeply diverging, and strongly supported clades, while the nuclear gene was congruent for three of these clades. There were no major differences among individuals from the different clades in any of the morphological features examined. Thus, the underlying genetic diversity in C. polypinum is greater than what its morphology suggests, indicating that morphology and genetics are not congruent in this organism. Furthermore, because the clades identified by the mitochondrial marker are so genetically diverse and are confirmed by a conserved nuclear marker in at least three cases, we propose that C. polypinum be designated as a "cryptic species complex." Our results provide another example where species diversity can be underestimated in microbial eukaryotes when using only morphological criteria to estimate species richness.

Expanding character sampling for ciliate phylogenetic inference using mitochondrial SSU-rDNA as a molecular marker

Molecular systematics of ciliates, particularly at deep nodes, has largely focused on increasing taxon sampling using the nuclear small subunit rDNA (nSSU-rDNA) locus. These previous analyses have generally been congruent with morphologically-based classifications, although there is extensive non-monophyly at many levels. However, caution is needed in interpreting these results as nSSU-rDNA is just a single molecular marker. Here the mitochondrial small subunit rDNA (mtSSU-rDNA) is evaluated for deep ciliate nodes using the Colpodea as an example. Overall, well-supported nodes in the mtSSU-rDNA and concatenated topologies are well supported in the nSSU-rDNA topology; e.g., the non-monophyly of the Cyrtolophosidida. The two moderately- to well-supported incongruences between the loci are the placement of the Sorogenida andColpoda aspera.Our analyses of mtSSU-rDNA support the conclusion, originally derived from nSSU-rDNA, that the morphological characters used in taxonomic circumscriptions of the Colpodea represent a mixture of ancestral and derived states. This demonstration of the efficacy of the mtSSU-rDNA will enable phylogenetic reconstructions of deep nodes in the ciliate tree of life to move from a single-locus to a multi-locus approach.

An SSU rDNA barcoding approach to the diversity of marine interstitial cercozoans, including descriptions of four novel genera and nine novel species

Environmental DNA surveys have revealed a great deal of hidden diversity within the Cercozoa. An investigation into the biodiversity of heterotrophic flagellates in marine benthic habitats of British Columbia, Canada, demonstrated the presence of several undescribed taxa with morphological features that resemble the cercozoan genera Cryothecomonas and Protaspis. Nine novel species of marine interstitial cercozoans are described that are distributed into five genera, four of which are new. Phylogenetic analyses of small subunit rDNA sequences derived from two uncultured isolates of Protaspis obliqua and nine novel cercozoan species (within four novel genera) provided organismal anchors that helped establish the cellular identities of several different environmental sequence clades. These data, however, also showed that the rarity of distinctive morphological features in cryomonads, and other groups of cercozoans, makes the identification and systematics of the group very difficult. Therefore, a DNA barcoding approach was applied as a diagnostic tool for species delimitation that used a 618 bp region at the 5' end of the SSU rDNA sequence. Nucleotide sequence analysis of this region showed high intergeneric sequence divergences of about 7% and very low intraspecific sequence divergences of 0-0.5%; phylogenetic analyses inferred from this barcoding region showed very similar tree topologies to those inferred from the full-length sequence of the gene. Overall, this study indicated that the 618 bp barcoding region of SSU rDNA sequences is a useful molecular signature for understanding the biodiversity and interrelationships of marine benthic cercozoans.

Relative diversity and community structure of ciliates in stream biofilms according to molecular and microscopy methods

Ciliates are an important component of aquatic ecosystems, acting as predators of bacteria and protozoa and providing nutrition for organisms at higher trophic levels. Understanding of the diversity and ecological role of ciliates in stream biofilms is limited, however. Ciliate diversity in biofilm samples from four streams subject to different impacts by human activity was assessed using microscopy and terminal restriction fragment length polymorphism (T-RFLP) analysis of 18S rRNA sequences. Analysis of 3' and 5' terminal fragments yielded very similar estimates of ciliate diversity. The diversity detected using microscopy was consistently lower than that suggested by T-RFLP analysis, indicating the existence of genetic diversity not apparent by morphological examination. Microscopy and T-RFLP analyses revealed similar relative trends in diversity between different streams, with the lowest level of biofilm-associated ciliate diversity found in samples from the least-impacted stream and the highest diversity in samples from moderately to highly impacted streams. Multivariate analysis provided evidence of significantly different ciliate communities in biofilm samples from different streams and seasons, particularly between a highly degraded urban stream and less impacted streams. Microscopy and T-RFLP data both suggested the existence of widely distributed, resilient biofilm-associated ciliates as well as ciliate taxa restricted to sites with particular environmental conditions, with cosmopolitan taxa being more abundant than those with restricted distributions. Differences between ciliate assemblages were associated with water quality characteristics typical of urban stream degradation and may be related to factors such as nutrient availability and macroinvertebrate communities. Microscopic and molecular techniques were considered to be useful complementary approaches for investigation of biofilm ciliate communities.