Species delimitation for the molecular taxonomy and ecology of the widely distributed microbial eukaryote genus Euplotes (Alveolata, Ciliophora)

Diverse patterns of correspondence between protist metabarcodes and protist metagenome-assembled genomes

Two common approaches to study the composition of environmental protist communities are metabarcoding and metagenomics. Raw metabarcoding data are usually processed into Operational Taxonomic Units (OTUs) or amplicon sequence variants (ASVs) through clustering or denoising approaches, respectively. Analogous approaches are used to assemble metagenomic reads into metagenome-assembled genomes (MAGs). Understanding the correspondence between the data produced by these two approaches can help to integrate information between the datasets and to explain how metabarcoding OTUs and MAGs are related with the underlying biological entities they are hypothesised to represent. MAGs do not contain the commonly used barcoding loci, therefore sequence homology approaches cannot be used to match OTUs and MAGs. We made an attempt to match V9 metabarcoding OTUs from the 18S rRNA gene (V9 OTUs) and MAGs from the Tara Oceans expedition based on the correspondence of their relative abundances across the same set of samples. We evaluated several metrics for detecting correspondence between features in these two datasets and developed controls to filter artefacts of data structure and processing. After selecting the best-performing metrics, ranking the V9 OTU/MAG matches by their proportionality/correlation coefficients and applying a set of selection criteria, we identified candidate matches between V9 OTUs and MAGs. In some cases, V9 OTUs and MAGs could be matched with a one-to-one correspondence, implying that they likely represent the same underlying biological entity. More generally, matches we observed could be classified into 4 scenarios: one V9 OTU matches many MAGs; many V9 OTUs match many MAGs; many V9 OTUs match one MAG; one V9 OTU matches one MAG. Notably, we found some instances in which different OTU-MAG matches from the same taxonomic group were not classified in the same scenario, with all four scenarios possible even within the same taxonomic group, illustrating that factors beyond taxonomic lineage influence the relationship between OTUs and MAGs. Overall, each scenario produces a different interpretation of V9 OTUs, MAGs and how they compare in terms of the genomic and ecological diversity they represent.

Nontriplet feature of genetic code in Euplotes ciliates is a result of neutral evolution

The triplet nature of the genetic code is considered a universal feature of known organisms. However, frequent stop codons at internal mRNA positions in Euplotes ciliates ultimately specify ribosomal frameshifting by one or two nucleotides depending on the context, thus posing a nontriplet feature of the genetic code of these organisms. Here, we sequenced transcriptomes of eight Euplotes species and assessed evolutionary patterns arising at frameshift sites. We show that frameshift sites are currently accumulating more rapidly by genetic drift than they are removed by weak selection. The time needed to reach the mutational equilibrium is several times longer than the age of Euplotes and is expected to occur after a several-fold increase in the frequency of frameshift sites. This suggests that Euplotes are at an early stage of the spread of frameshifting in expression of their genome. In addition, we find the net fitness burden of frameshift sites to be noncritical for the survival of Euplotes. Our results suggest that fundamental genome-wide changes such as a violation of the triplet character of genetic code can be introduced and maintained solely by neutral evolution.

Ciliate Morpho-Taxonomy and Practical Considerations before Deploying Metabarcoding to Ciliate Community Diversity Surveys in Urban Receiving Waters

Disentangling biodiversity and community assembly effects on ecosystem function has always been an important topic in ecological research. The development and application of a DNA metabarcoding method has fundamentally changed the way we describe prokaryotic communities and estimate biodiversity. Compared to prokaryotes (bacteria and archaea), the eukaryotic microbes (unicellular eukaryotes) also fulfill extremely important ecological functions in different ecosystems regarding their intermediate trophic positions. For instance, ciliated microbes (accounting for a substantial portion of the diversity of unicellular eukaryotes) perform pivotal roles in microbial loops and are essential components in different ecosystems, especially in water purification processes. Therefore, the community composition of ciliated species has been widely utilized as a proxy for water quality and biological assessment in urban river ecosystems and WWTPs (wastewater treatment plants). Unfortunately, investigating the dynamic changes and compositions in ciliate communities relies heavily on existing morpho-taxonomical descriptions, which is limited by traditional microscopic approaches. To deal with this dilemma, we discuss the DNA-based taxonomy of ciliates, the relative merits and challenges of deploying its application using DNA metabarcoding for surveys of ciliate community diversity in urban waterbodies, and provide suggestions for minimizing relevant sources of biases in its implementation. We expect that DNA metabarcoding could untangle relationships between community assembly and environmental changes affecting ciliate communities. These analyses and discussions could offer a replicable method in support of the application of evaluating communities of ciliated protozoa as indicators of urban freshwater ecosystems.

Primary Structure and Coding Genes of Two Pheromones from the Antarctic Psychrophilic Ciliate, Euplotes focardii

In ciliates, diffusible cell type-specific pheromones regulate cell growth and mating phenomena acting competitively in both autocrine and heterologous fashion. In Euplotes species, these signaling molecules are represented by species-specific families of structurally homologous small, disulfide-rich proteins, each specified by one of a series of multiple alleles that are inherited without relationships of dominance at the mat-genetic locus of the germinal micronuclear genome, and expressed as individual gene-sized molecules in the somatic macronuclear genome. Here we report the 85-amino acid sequences and the full-length macronuclear nucleotide coding sequences of two pheromones, designated Ef-1 and Ef-2, isolated from the supernatant of a wild-type strain of a psychrophilic species of Euplotes, E. focardii, endemic to Antarctic coastal waters. An overall comparison of the determined E. focardii pheromone and pheromone-gene structures with their homologs from congeneric species provides an initial picture of how an evolutionary increase in the complexity of these structures accompanies Euplotes speciation.

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.

Diversity Patterns of Protists Are Highly Affected by Methods Disentangling Biological Variants: A Case Study in Oligotrich (s.l.) Ciliates

Protists are a dominant group in marine microplankton communities and play important roles in energy flux and nutrient cycling in marine ecosystems. Environmental sequences produced by high-throughput sequencing (HTS) methods are increasingly used for inferring the diversity and distribution patterns of protists. However, studies testing whether methods disentangling biological variants affect the diversity and distribution patterns of protists using field samples are insufficient. Oligotrich (s.l.) ciliates are one group of the abundant and dominant planktonic protists in coastal waters and open oceans. Using oligotrich (s.l.) ciliates in field samples as an example, the present study indicates that DADA2 performs better than SWARM, UNOISE, UPARSE, and UCLUST for inferring diversity patterns of oligotrich (s.l.) ciliates in the Pearl River Estuary and surrounding regions. UPARSE and UNOISE might underestimate species richness. SWARM might not be suitable for the resolution of alpha diversity owing to its rigorous clustering and sensitivity to sequence variations. UCLUST with 99% clustering threshold overestimates species richness, and the beta diversity pattern inferred by DADA2 is more reasonable than that of the other methods. Additionally, salinity is shown to be one of the key factors responsible for variations in the community distribution of ciliates, but infrequent marine-freshwater transitions occurred during evolutionary terms of this group.

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.

Chromosome organization and gene expansion in the highly fragmented genome of the ciliate Strombidium stylifer

Chromosomes are well-organized carriers of genetic information in eukaryotes and are usually quite long, carrying hundreds and thousands of genes. Intriguingly, a clade of single-celled ciliates, Spirotrichea, feature nanochromosomes-also called "gene-sized chromosomes". These chromosomes predominantly carry only one gene, flanked by short telomere sequences. However, the organization and copy number variation of the chromosomes in these highly fragmented genomes remain unexplored in many groups of Spirotrichea, including the marine Strombidium. Using deep genome sequencing, we assembled the macronuclear genome of Strombidium stylifer into more than 18,000 nanochromosomes (~2.4 Kb long on average). Our results show that S. stylifer occupies an intermediate position during the evolutionary history of Strombidium lineage and experienced significant expansions in several gene families related to guanyl ribonucleotide binding. Based on the nucleotide distribution bias analysis and conserved motifs search in non-genic regions, we found that the subtelomeric regions have a conserved adenine-thymine (AT)-rich sequence motif. We also found that the copy number of nanochromosomes lacks precise regulation. This work sheds light on the unique features of chromosome structure in eukaryotes with highly fragmented genomes and reveals that a rather specialized evolutionary strategy at the genomic level has resulted in great diversity within the ciliated lineages.

Molecular phylogeny of a new gonostomatid ciliate revealing a discrepancy between interphasic and cell divisional patterns (Ciliophora, Hypotricha)

The morphology, ontogenesis, and molecular phylogeny of a new ciliate, Gonostomum koreanum n. sp., discovered in a terrestrial moss sample from South Korea, were investigated. Morphologically, it is characterized by a gonostomatid oral apparatus, two macronuclear nodules, six frontoventral rows, the two rightmost of which (frontoventral rows V and VI) extend posteriorly to near pretransverse and transverse cirri, and three dorsal kineties each with a single caudal cirrus posteriorly. The new species is easily confused with members of the genus Metagonostomum because of the long frontoterminal cirral row but differs mainly in the de novo (vs. intrakinetal) origin of anlage VI, a character found only in Gonostomum and Paragonostomum. To solve the discrepancy between the interphasic and ontogenetic patterns, we additionally performed morphological and multigene analyses on three gonostomatid species, namely Gonostomum koreanum n. sp. and its morphologically (M. gonostomoida) and ontogenetically (G. kuehnelti) most similar species. The multigene analyses show that the new species is closely related to G. kuehnelti and the core gonostomatids consists of five groups based on the origin of the frontoventral rows.

Estimating hybridization as a consequence of climatic fluctuations: A case study of two geometridae species

In some cases, the phylogenetic analysis based on the mitochondrial DNA (mtDNA) and the nuclear DNA (ncDNA) are discordant. There are three major causes of the discordance within insects, including hybridization, incomplete lineage sorting (ILS) and infection by Wolbachia. In this study, we used a combination of multilocus and coalescent analyses to explore these processes occurred during the evolutionary history of Limbatochlamys rosthorni Rothschild, 1894 and Limbatochlamys pararosthorni Han and Xue, 2005. The ncDNA phylogenetic tree supported two reciprocally monophyletic species, whereas the mtDNA results failed to reveal such a structure and revealed an extensive level of admixture between two species. Because of very low Wolbachia infection rates (<20%), we firstly excluded this reason for the discordance. The fixed nucleotide differences and large genetic distances (1.5-2.5%) at the ncDNA genes suggested that the lineage sorting process between these two species is nearly complete and two species have experienced a prolonged period of independent evolution. Thus, we secondly excluded ILS. Sharing haplotypes, mtDNA gene flow occurring and the transitional samples with morphological features supported hybridization. The distribution contraction during glaciations and postglacial distribution expansion might have facilitated hybridization. Taken together, our study indicates that the current genetic structure of L. rosthorni and L. pararosthorni is the results of contraction and fragmentation into separated refugia during glaciations, followed by postglacial expansion and admixture.

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.

Morphology, ultrastructure, genomics, and phylogeny of Euplotes vanleeuwenhoeki sp. nov. and its ultra-reduced endosymbiont "Candidatus Pinguicoccus supinus" sp. nov

Taxonomy is the science of defining and naming groups of biological organisms based on shared characteristics and, more recently, on evolutionary relationships. With the birth of novel genomics/bioinformatics techniques and the increasing interest in microbiome studies, a further advance of taxonomic discipline appears not only possible but highly desirable. The present work proposes a new approach to modern taxonomy, consisting in the inclusion of novel descriptors in the organism characterization: (1) the presence of associated microorganisms (e.g.: symbionts, microbiome), (2) the mitochondrial genome of the host, (3) the symbiont genome. This approach aims to provide a deeper comprehension of the evolutionary/ecological dimensions of organisms since their very first description. Particularly interesting, are those complexes formed by the host plus associated microorganisms, that in the present study we refer to as "holobionts". We illustrate this approach through the description of the ciliate Euplotes vanleeuwenhoeki sp. nov. and its bacterial endosymbiont "Candidatus Pinguicoccus supinus" gen. nov., sp. nov. The endosymbiont possesses an extremely reduced genome (~ 163 kbp); intriguingly, this suggests a high integration between host and symbiont.

Two new species of Euplotes with cirrotype-9, Euplotes foissneri sp. nov. and Euplotes warreni sp. nov. (Ciliophora, Spirotrichea, Euplotida), from the coasts of Patagonia: implications from their distant, early and late branching in the Euplotes phylogenetic tree

Two new Euplotes species have been isolated from cold shallow sandy sediments of the extreme Southern Chilean coasts: Euplotes foissneri sp. nov., from a low-salinity site at Puerto Natales on the Pacific coast, and Euplotes warreni sp. nov., from a marine site at Punta Arenas on the Atlantic coast. Euplotes foissneri has a medium body size (53×36 µm in vivo), a dorsal surface marked by six prominent ridges, a double dargyrome, six dorsal and two ventrolateral kineties, a buccal field extending to about 3/4 of the body length, an adoral zone composed of 28-32 membranelles, and nine fronto-ventral, five transverse and two or three caudal cirri. The bulky, hook-, horseshoe- or 3-shaped macronucleus is associated with one sub-spherical micronucleus. The central body region hosts taxonomically unidentified endosymbiotic eubacteria. Euplotes warreni has a small body size (39×27 µm in vivo), a smooth dorsal surface marked by three deep grooves, a double dargyrome, four dorsal and two ventrolateral kineties, a buccal field extending to about 2/3 of the body length, an adoral zone composed of 23-25 adoral membranelles, and nine fronto-ventral, five transverse and three caudal cirri. The macronucleus is hook- or C-shaped and associated with one spherical micronucleus. Endosymbiotic bacteria belonging to the genus Francisella reside preferentially in the anterior cell region. Both species lack the fronto-ventral cirrus numbered 'V/2', whereby their cirrotype-9 conforms to the so-called 'pattern I', which is the basic distinctive trait of the genus Euplotopsis Borror and Hill, 1995. Phylogenetic analyses of small subunit rRNA gene sequences, however, classify E. warreni into its own early branching clade and E. foissneri into a late branching clade. This indicates a polyphyletic nature and taxonomic inconsistency of the genus Euplotopsis, which was erected to include Euplotes species with cirrotype-9 pattern I.

Environmental Parameters and Substrate Type Drive Microeukaryotic Community Structure During Short-Term Experimental Colonization in Subtropical Eutrophic Freshwaters

Microeukaryotes are key components of aquatic ecosystems and play crucial roles in aquatic food webs. However, influencing factors and potential assembly mechanisms for microeukaryotic community on biofilms are rarely studied. Here, those of microeukaryotic biofilms in subtropical eutrophic freshwaters were investigated for the first time based on 2,585 operational taxonomic units (OTUs) from 41 samples, across different environmental conditions and substrate types. Following conclusions were drawn: (1) Environmental parameters were more important than substrate types in structuring microeukaryotic community of biofilms in subtropical eutrophic freshwaters. (2) In the fluctuating river, there was a higher diversity of OTUs and less predictability of community composition than in the stable lake. Sessile species were more likely to be enriched on smooth surfaces of glass slides, while both free-swimming and attached organisms occurred within holes inside PFUs (polyurethane foam units). (3) Both species sorting and neutral process were mechanisms for assembly of microeukaryotic biofilms, but their importance varied depending on different habitats and substrates. (4) The effect of species sorting was slightly higher than the neutral process in river biofilms due to stronger environmental filtering. Species sorting was a stronger force structuring communities on glass slides than PFUs with more niche availability. Our study sheds light on assembly mechanisms for microeukaryotic community on different habitat and substrate types, showing that the resulting communities are determined by both sets of variables, in this case primarily habitat type. The balance of neutral process and species sorting differed between habitats, but the high alpha diversity of microeukaryotes in both led to similar sets of lifecycle traits being selected for in each case.

Two parasitic ciliates (Protozoa: Ciliophora: Phyllopharyngea) isolated from respiratory-mucus of an unhealthy beluga whale: characterization, phylogeny and an assessment of morphological adaptations

Ciliates occur in the blowholes of marine mammals, but our understanding of their biology is poor. Consequently, we investigated an infestation of ciliates in an unhealthy, captive beluga whale that was exhibiting accelerated breathing, leukocytosis and expulsion of unusually large amounts of viscous sputum. This sputum contained ~104 ciliates mL-1 (when healthy, numbers were ten- to 100-fold lower). One known ciliate species, Planilamina ovata, is fully characterized, and a new species, Kyaroikeus paracetarius sp. nov., is here described. The new species is established based on its larger number of left kineties over its only congener. Sequences of small-subunit rDNA, large-subunit rDNA and ITS1-5.8S-ITS2 regions of these two taxa were used in phylogenetic analyses, inferring that Kyaroikeus and Planilamina have close affinity with the free-living family Dysteriidae, contradicting their morphology-based assignment to the family Kyaroikeidae. We suggest that Kyaroikeidae be relegated to subfamily status. Finally, by comparing parasitic species with free-living taxa, we suggest how these ciliates have adapted to their unique environment and how they may have initially invaded the host. We provide essential data and concepts for the continued evaluation of ciliate-parasites in whale blowholes.

Assessing the utility of Hsp90 gene for inferring evolutionary relationships within the ciliate subclass Hypotricha (Protista, Ciliophora)

BACKGROUND

Although phylogenomic analyses are increasingly used to reveal evolutionary relationships among ciliates, relatively few nuclear protein-coding gene markers have been tested for their suitability as candidates for inferring phylogenies within this group. In this study, we investigate the utility of the heat-shock protein 90 gene (Hsp90) as a marker for inferring phylogenetic relationships among hypotrich ciliates.

RESULTS

A total of 87 novel Hsp90 gene sequences of 10 hypotrich species were generated. Of these, 85 were distinct sequences. Phylogenetic analyses based on these data showed that: (1) the Hsp90 gene amino acid trees are comparable to the small subunit rDNA tree for recovering phylogenetic relationships at the rank of class, but lack sufficient phylogenetic signal for inferring evolutionary relationships at the genus level; (2) Hsp90 gene paralogs are recent and therefore unlikely to pose a significant problem for recovering hypotrich clades; (3) definitions of some hypotrich orders and families need to be revised as their monophylies are not supported by various gene markers; (4) The order Sporadotrichida is paraphyletic, but the monophyly of the "core" Urostylida is supported; (5) both the subfamily Oxytrichinae and the genus Urosoma seem to be non-monophyletic, but monophyly of Urosoma is not rejected by AU tests.

CONCLUSIONS

Our results for the first time demonstrate that the Hsp90 gene is comparable to SSU rDNA for recovering phylogenetic relationships at the rank of class, and its paralogs are unlikely to pose a significant problem for recovering hypotrich clades. This study shows the value of careful gene marker selection for phylogenomic analyses of ciliates.

Characterization of Simple Sequence Repeats (SSRs) in Ciliated Protists Inferred by Comparative Genomics

Simple sequence repeats (SSRs) are prevalent in the genomes of all organisms. They are widely used as genetic markers, and are insertion/deletion mutation hotspots, which directly influence genome evolution. However, little is known about such important genomic components in ciliated protists, a large group of unicellular eukaryotes with extremely long evolutionary history and genome diversity. With recent publications of multiple ciliate genomes, we start to get a chance to explore perfect SSRs with motif size 1-100 bp and at least three motif repeats in nine species of two ciliate classes, Oligohymenophorea and Spirotrichea. We found that homopolymers are the most prevalent SSRs in these A/T-rich species, with AAA (lysine, charged amino acid; also seen as an SSR with one-adenine motif repeated three times) being the codons repeated at the highest frequencies in coding SSR regions, consistent with the widespread alveolin proteins rich in lysine repeats as found in Tetrahymena. Micronuclear SSRs are universally more abundant than the macronuclear ones of the same motif-size, except for the 8-bp-motif SSRs in extensively fragmented chromosomes. Both the abundance and A/T content of SSRs decrease as motif-size increases, while the abundance is positively correlated with the A/T content of the genome. Also, smaller genomes have lower proportions of coding SSRs out of all SSRs in Paramecium species. This genome-wide and cross-species analysis reveals the high diversity of SSRs and reflects the rapid evolution of these simple repetitive elements in ciliate genomes.

Comparative Studies on the Polymorphism and Copy Number Variation of mtSSU rDNA in Ciliates (Protista, Ciliophora): Implications for Phylogenetic, Environmental, and Ecological Research

While nuclear small subunit ribosomal DNA (nSSU rDNA) is the most commonly-used gene marker in studying phylogeny, ecology, abundance, and biodiversity of microbial eukaryotes, mitochondrial small subunit ribosomal DNA (mtSSU rDNA) provides an alternative. Recently, both copy number variation and sequence variation of nSSU rDNA have been demonstrated for diverse organisms, which can contribute to misinterpretation of microbiome data. Given this, we explore patterns for mtSSU rDNA among 13 selected ciliates (representing five classes), a major component of microbial eukaryotes, estimating copy number and sequence variation and comparing to that of nSSU rDNA. Our study reveals: (1) mtSSU rDNA copy number variation is substantially lower than that for nSSU rDNA; (2) mtSSU rDNA copy number ranges from 1.0 × 10⁴ to 8.1 × 10⁵; (3) a most common sequence of mtSSU rDNA is also found in each cell; (4) the sequence variation of mtSSU rDNA are mainly indels in poly A/T regions, and only half of species have sequence variation, which is fewer than that for nSSU rDNA; and (5) the polymorphisms between haplotypes of mtSSU rDNA would not influence the phylogenetic topology. Together, these data provide more insights into mtSSU rDNA as a powerful marker especially for microbial ecology studies.

Comparative Transcriptome Analyses during the Vegetative Cell Cycle in the Mono-Cellular Organism Pseudokeronopsis erythrina (Alveolata, Ciliophora)

Studies focusing on molecular mechanisms of cell cycles have been lagging in unicellular eukaryotes compared to other groups. Ciliates, a group of unicellular eukaryotes, have complex cell division cycles characterized by multiple events. During their vegetative cell cycle, ciliates undergo macronuclear amitosis, micronuclear mitosis, stomatogenesis and somatic cortex morphogenesis, and cytokinesis. Herein, we used the hypotrich ciliate Pseudokeronopsis erythrina, whose morphogenesis has been well studied, to examine molecular mechanisms of ciliate vegetative cell cycles. Single-cell transcriptomes of the growth (G) and cell division (D) stages were compared. The results showed that (i) More than 2051 significantly differentially expressed genes (DEGs) were detected, among which 1545 were up-regulated, while 256 were down-regulated at the D stage. Of these, 11 randomly picked DEGs were validated by reverse transcription quantitative polymerase chain reaction (RT-qPCR); (ii) Enriched DEGs during the D stage of the vegetative cell cycle of P. erythrina were involved in development, cortex modifications, and several organelle-related biological processes, showing correspondence of molecular evidence to morphogenetic changes for the first time; (iii) Several individual components of molecular mechanisms of ciliate vegetative division, the sexual cell cycle and cellular regeneration overlap; and (iv) The P. erythrina cell cycle and division have the same essential components as other eukaryotes, including cyclin-dependent kinases (CDKs), cyclins, and genes closely related to cell proliferation, indicating the conserved nature of this biological process. Further studies are needed focusing on detailed inventory and gene interactions that regulate specific ciliated cell-phase events.

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.

Spatio-temporal patterns of zooplankton in a main-stem dam affected tributary: a case study in the Xiangxi River of the Three Gorges Reservoir, China

As the ecologically important recipient channels for riverine ecosystems, tributaries provide unique microhabitats for microorganisms, among which zooplankton constitutes the most important heterotrophic organisms. In particular, the reduced water velocity caused by dams is more favorable for zooplankton development; therefore, dammed rivers are expected to support extremely diverse and abundant zooplankton communities and notably different spatiotemporal distribution patterns. So far, however, only very few molecular studies support these assumptions. Using high-throughput sequencing, a high number of 350 operational taxonomic units (OTUs; 97% cutoff) were retrieved from 30 samples collected in the Xiangxi River, the nearest large tributary upstream of the Three Gorges Dam. Zooplankton did not show significant spatial distribution in the channel. Instead, the community structures varied significantly over sampling dates, corroborating the seasonal patterns found in lakes and ponds in the subtropical zone. As expected, the community compositions were deterministically governed by environmental filtering processes (phylogenetic clustering), in which water velocity appeared to be much less important than other investigated environmental factors. Moreover, most of the detected phylotypes (OTUs) had a relatively high (>90%) sequence similarity to previously deposited sequences, suggesting a mediocre degree of genetic novelty within the zooplankton communities in the Xiangxi River.

Genome analyses of the new model protist Euplotes vannus focusing on genome rearrangement and resistance to environmental stressors

As a model organism for studies of cell and environmental biology, the free-living and cosmopolitan ciliate Euplotes vannus shows intriguing features like dual genome architecture (i.e., separate germline and somatic nuclei in each cell/organism), "gene-sized" chromosomes, stop codon reassignment, programmed ribosomal frameshifting (PRF) and strong resistance to environmental stressors. However, the molecular mechanisms that account for these remarkable traits remain largely unknown. Here we report a combined analysis of de novo assembled high-quality macronuclear (MAC; i.e., somatic) and partial micronuclear (MIC; i.e., germline) genome sequences for E. vannus, and transcriptome profiling data under varying conditions. The results demonstrate that: (a) the MAC genome contains more than 25,000 complete "gene-sized" nanochromosomes (~85 Mb haploid genome size) with the N50 ~2.7 kb; (b) although there is a high frequency of frameshifting at stop codons UAA and UAG, we did not observe impaired transcript abundance as a result of PRF in this species as has been reported for other euplotids; (c) the sequence motif 5'-TA-3' is conserved at nearly all internally-eliminated sequence (IES) boundaries in the MIC genome, and chromosome breakage sites (CBSs) are duplicated and retained in the MAC genome; (d) by profiling the weighted correlation network of genes in the MAC under different environmental stressors, including nutrient scarcity, extreme temperature, salinity and the presence of ammonia, we identified gene clusters that respond to these external physical or chemical stimulations, and (e) we observed a dramatic increase in HSP70 gene transcription under salinity and chemical stresses but surprisingly, not under temperature changes; we link this temperature-resistance to the evolved loss of temperature stress-sensitive elements in regulatory regions. Together with the genome resources generated in this study, which are available online at Euplotes vannus Genome Database (http://evan.ciliate.org), these data provide molecular evidence for understanding the unique biology of highly adaptable microorganisms.

Further analyses of variation of ribosome DNA copy number and polymorphism in ciliates provide insights relevant to studies of both molecular ecology and phylogeny

Sequence-based approaches, such as analyses of ribosome DNA (rDNA) clone libraries and high-throughput amplicon sequencing, have been used extensively to infer evolutionary relationships and elucidate the biodiversity in microbial communities. However, recent studies demonstrate both rDNA copy number variation and intra-individual (intra-genomic) sequence variation in many organisms, which challenges the application of the rDNA-based surveys. In ciliates, an ecologically important clade of microbial eukaryotes, rDNA copy number and sequence variation are rarely studied. In the present study, we estimate the intraindividual small subunit rDNA (SSU rDNA) copy number and sequence variation in a wide range of taxa covering nine classes and 18 orders of the phylum Ciliophora. Our studies reveal that: (i) intra-individual sequence variation of SSU rDNA is ubiquitous in all groups of ciliates detected and the polymorphic level varies among taxa; (ii) there is a most common version of SSU rDNA sequence in each cell that is highly predominant and may represent the germline micronuclear template; (iii) compared with the most common version, other variant sequences differ in only 1-3 nucleotides, likely generated during macronuclear (somatic) amplification; (iv) the intra-cell sequence variation is unlikely to impact phylogenetic analyses; (v) the rDNA copy number in ciliates is highly variable, ranging from 10³ to 10⁶, with the highest record in Stentor roeselii. Overall, these analyses indicate the need for careful consideration of SSU rDNA variation in analyses of the role of ciliates in ecosystems.

Time-course analysis of nuclear events during conjugation in the marine ciliate Euplotes vannus and comparison with other ciliates (Protozoa, Ciliophora)

Ciliates represent a morphologically and genetically distinct group of single-celled eukaryotes that segregate germline and somatic functions into two types of nuclei and exhibit complex cytogenetic events during the sexual process of conjugation, which is under the control of the so-called "mating type systems". Studying conjugation in ciliates may provide insight into our understanding of the origins and evolution of sex and fertilization. In the present work, we studied in detail the sexual process of conjugation using the model species Euplotes vannus, and compared these nuclear events with those occurring in other ciliates. Our results indicate that in E. vannus: 1) conjugation requires about 75 hours to complete: the longest step is the development of the new macronucleus (ca. 64h), followed by the nuclear division of meiosis I (5h); the mitotic divisions usually take only 2h; 2) there are three prezygotic divisions (mitosis and meiosis I and II), and two of the eight resulting nuclei become pronuclei; 3) after the exchange and fusion of the pronuclei, two postzygotic divisions occur; two of the four products differentiate into the new micronucleus and macronucleus, respectively, and the parental macronucleus degenerates completely; 4) comparison of the nuclear events during conjugation in different ciliates reveals that there are generally three prezygotic divisions while the number of postzygotic divisions is highly variable. These results can serve as reference to investigate the mating type system operating in this species and to analyze genes involved in the different steps of the sexual process.

New contribution to the species-rich genus Euplotes: Morphology, ontogeny and systematic position of two species (Ciliophora; Euplotia)

The morphology, ontogeny and phylogeny of two Euplotes species, E. estuarinus sp. nov. and a population of E. platystoma Dragesco and Dragesco-Kernéis, 1986, both collected from tropical brackish waters in south China, were investigated based on living morphology, ciliary pattern and molecular data. Euplotes estuarinus sp. nov. is small (about 60 × 40 μm in vivo), has a dargyrome of the double-eurystomus type, and the transverse cirri are arranged in two groups, with two left and three right ones. The original description of the poorly known species, E. platystoma, is brief, and the species was never investigated using live observation and molecular methods Hence, we provided a detailed redescription. Some stages of their morphogenesis were observed which proceed in the same pattern as in their congeners. The new species E. estuarinus sp. nov. clusters with E. curdsi, differing only by 1 bp in their SSU rRNA gene sequences, which is likely due to the recent speciation event and the limited resolution of the SSU rRNA gene at species level in this group as the two species are clearly morphologically distinct.