Delimiting operational taxonomic units for assessing ciliate environmental diversity using small-subunit rRNA gene sequences

Reproductive Tract Microbial Transitions from Late Gestation to Early Postpartum Using 16S rRNA Metagenetic Profiling in First-Pregnancy Heifers

Studies in recent years indicate that reproductive tract microbial communities are crucial for shaping mammals' health and reproductive outcomes. Following parturition, uterine bacterial contamination often occurs due to the open cervix, which may lead to postpartum uterine inflammatory diseases, especially in primiparous individuals. However, investigations into spatio-temporal microbial transitions in the reproductive tract of primigravid females remain limited. Our objective was to describe and compare the microbial community compositions in the vagina at late gestation and in the vagina and uterus at early postpartum in first-pregnancy heifers. Three swab samples were collected from 33 first-pregnancy Holstein Friesian heifers: one vaginal sample at gestation day 258 ± 4, and vaginal and uterine samples at postpartum day 7 ± 2. Each sample underwent 16S rRNA V4 region metagenetic analysis via Illumina MiSeq, with bioinformatics following Mothur MiSeq SOP. The reproductive tract bacterial communities were assigned to 1255 genus-level OTUs across 30 phyla. Dominant phyla, accounting for approximately 90% of the communities, included Proteobacteria, Firmicutes, Actinobacteria, Bacteroidetes, and Fusobacteria. However, the results revealed distinct shifts in microbial composition between the prepartum vagina (Vag-pre), postpartum vagina (Vag-post), and postpartum uterus (Utr-post). The Vag-pre and Utr-post microbial profiles were the most distinct. The Utr-post group had lower relative abundances of Proteobacteria but higher abundances of Bacteroidetes, Fusobacteria, and Tenericutes compared to Vag-pre, while Vag-post displayed intermediate values for these phyla, suggesting a transitional profile. Additionally, the Utr-post group exhibited lower bacterial richness and diversity compared to both Vag-pre and Vag-post. The unsupervised probabilistic Dirichlet Multinomial Mixtures model identified two distinct community types: most Vag-pre samples clustered into one type and Utr-post samples into another, while Vag-post samples were distributed evenly between the two. LEfSe analysis revealed distinct microbial profiles at the genus level. Overall, specific microbial markers were associated with anatomical and temporal transitions, revealing a dynamic microbial landscape during the first pregnancy and parturition. These differences highlight the complexity of these ecosystems and open new avenues for research in reproductive biology and microbial ecology.

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.

Overview of the Diversity, Phylogeny and Biogeography of Strombidiid Oligotrich Ciliates (Protista, Ciliophora), With a Brief Revision and a Key to the Known Genera

Strombidiids are common free-living ciliates that have colonized coastal and open oceanic waters across the world. In recent years, numerous new taxa and gene sequences of strombidiids have been reported, revealing a large diversity of both their morphologic and genetic features. Here, we compare the taxonomic characters of all genera in the family Strombidiidae, provide a key to their identification, and investigate their molecular phylogeny. In addition, we analyze their regional distribution based on faunal data accumulated in China and attempt to infer their global distribution based on SSU rRNA gene sequence data. The current work revises the systematics of strombidiids based on morphologic, phylogenetic, and biogeographic evidence and provides a genus-level review of marine strombidiids.

Improving eDNA-based protist diversity assessments using networks of amplicon sequence variants

Effective and precise grouping of highly similar sequences remains a major bottleneck in the evaluation of high-throughput sequencing datasets. Amplicon sequence variants (ASVs) offer a promising alternative that may supersede the widely used operational taxonomic units (OTUs) in environmental sequencing studies. We compared the performance of a recently developed pipeline based on the algorithm DADA2 for obtaining ASVs against a pipeline based on the algorithm SWARM for obtaining OTUs. Illumina-sequencing of 29 individual ciliate species resulted in up to 11 ASVs per species, while SWARM produced up to 19 OTUs per species. To improve the congruency between species diversity and molecular diversity, we applied sequence similarity networks (SSNs) for second-level sequence grouping into network sequence clusters (NSCs). At 100% sequence similarity in SWARM-SSNs, NSC numbers decreased from 7.9-fold overestimation without abundance filter, to 4.5-fold overestimation when an abundance filter was applied. For the DADA2-SSN approach, NSC numbers decreased from 3.5-fold to 3-fold overestimation. Rand index cluster analyses predicted best binning results between 97% and 94% sequence similarity for both DADA2-SSNs and SWARM-SSNs. Depending on the ecological questions addressed in an environmental sequencing study with protists we recommend ASVs as replacement for OTUs, best in combination with SSNs.

Are We Overestimating Protistan Diversity in Nature?

Documenting the immense diversity of single-celled, eukaryotic organisms (protists) has been a formidable challenge for ecologists. These species were originally defined by morphological criteria, but shortcomings of the morphospecies concept, and a bewildering array of sizes and cellular attributes, has made constructing a taxonomy that is useful for ecologists nearly impossible. Consequently, physiological and genetic information has been integrated to address these shortcomings, and to develop the framework of a unifying taxonomy. DNA sequence information, in particular, has revolutionized studies of protistan diversity. However, the exponential increase in sequence-based protistan species richness published from field surveys in recent years raises the question of whether we have moved beyond characterizing species-level diversity and begun to reveal intraspecies diversity. The answer to that question appears to be 'yes', at least for some protistan lineages. The need to document such microdiversity may be justified, but it is important for protistologists to recognize and acknowledge that possibility, and its consequences.

The ecology and diversity of microbial eukaryotes in geothermal springs

Decades of research into the Bacteria and Archaea living in geothermal spring ecosystems have yielded great insight into the diversity of life and organismal adaptations to extreme environmental conditions. Surprisingly, while microbial eukaryotes (protists) are also ubiquitous in many environments, their diversity across geothermal springs has mostly been ignored. We used high-throughput sequencing to illuminate the diversity and structure of microbial eukaryotic communities found in 160 geothermal springs with broad ranges in temperature and pH across the Taupō Volcanic Zone in New Zealand. Protistan communities were moderately predictable in composition and varied most strongly across gradients in pH and temperature. Moreover, this variation mirrored patterns observed for bacterial and archaeal communities across the same spring samples, highlighting that there are similar ecological constraints across the tree of life. While extreme pH values were associated with declining protist diversity, high temperature springs harbored substantial amounts of protist diversity. Although protists are often overlooked in geothermal springs and other extreme environments, our results indicate that such environments can host distinct and diverse protistan communities.

The Hidden Diversity of Flagellated Protists in Soil

Protists are among the most diverse and abundant eukaryotes in soil. However, gaps between described and sequenced protist morphospecies still present a pending problem when surveying environmental samples for known species using molecular methods. The number of sequences in the molecular PR² database (∼130,000) is limited compared to the species richness expected (>1 million protist species) - limiting the recovery rate. This is important, since high throughput sequencing (HTS) methods are used to find associative patterns between functional traits, taxa and environmental parameters. We performed HTS to survey soil flagellates in 150 grasslands of central Europe, and tested the recovery rate of ten previously isolated and cultivated cercomonad species, among locally found diversity. We recovered sequences for reference soil flagellate species, but also a great number of their phylogenetically evaluated genetic variants, among rare and dominant taxa with presumably own biogeography. This was recorded among dominant (cercozoans, Sandona), rare (apusozoans) and a large hidden diversity of predominantly aquatic protists in soil (choanoflagellates, bicosoecids) often forming novel clades associated with uncultured environmental sequences. Evaluating the reads, instead of the OTUs that individual reads are usually clustered into, we discovered that much of this hidden diversity may be lost due to clustering.

Discrepancies Between Molecular and Morphological Databases of Soil Ciliates Studied for Temperate Grasslands of Central Europe

By measuring the change in soil protist communities, the effect of human land use on grasslands can be monitored to promote sustainable ecosystem functioning. Protists form the active link in the rhizosphere between the plant roots and higher trophic organisms; however, only few morphological species and their ecological values have yet been described in this context. To investigate the communicability between morphological and molecular databases used in the molecular barcoding of protists and in the biomonitoring of grassland soil, the present high-throughput sequencing (HTS) study (N=150) covered the area of central Europe (mesoscale) known to be well studied for ciliated protists. HTS delivered 2,404 unique reads identifying taxa in all major ciliophoran classes but exact reference matches were few. The study identified clear discrepancies between databases for well-studied taxa, where molecular databases contained multiple gene variants for single morphospecies of dominant taxa. Gene variants presented own biogeography - the eukaryotic microdiversity along gradients (e.g., land-use intensity, soil water). It is possible that many of the so called novel phylogenetic lineages and hidden diversity pointed out in environmental surveys could be evidence for the severe lack of molecular data for already known and morphologically described species, present in morphological databases.

Patterns and Drivers of Vertical Distribution of the Ciliate Community from the Surface to the Abyssopelagic Zone in the Western Pacific Ocean

The deep sea is one of the largest but least understood ecosystems on earth. Knowledge about the diversity and distribution patterns as well as drivers of microbial eukaryote (including ciliates) along the water column, particularly below the photic zone, is scarce. In this study, we investigated the diversity of pelagic ciliates, the main group of marine microeukaryotes, their vertical distribution from the surface to the abyssopelagic zone, as well as their horizontal distribution over a distance of 1,300 km in the Western Pacific Ocean, using high-throughput DNA and cDNA (complementary DNA) sequencing. No distance-decay relationship could be detected along the horizontal scale; instead, a distinct vertical distribution within the ciliate communities was revealed. The alpha diversity of the ciliate communities in the deep chlorophyll maximum (DCM) and the 200 m layer turned out to be significantly higher compared with the other water layers. The ciliate communities in the 200 m water layer appeared to be more similar to those in deeper layers from 1,000 m to about 5,000 m than to the surface and DCM ciliate communities. Dominant species in the bathypelagic and abyssopelagic zone, particularly some parasites, were also detected in the 200 m layer, but were almost absent in the surface layer. The 200 m layer, therefore, seems to be an important “species bank” for deep ocean layers. Statistical analyses further revealed significant effects of temperature and chlorophyll a on the partitioning of ciliate diversity, indicating that environmental factors are a stronger force in shaping marine pelagic ciliate communities than the geographic distance.

Yeasts dominate soil fungal communities in three lowland Neotropical rainforests

Forest soils typically harbour a vast diversity of fungi, but are usually dominated by filamentous (hyphae-forming) taxa. Compared to temperate and boreal forests, though, we have limited knowledge about the fungal diversity in tropical rainforest soils. Here we show, by environmental metabarcoding of soil samples collected in three Neotropical rainforests, that Yeasts dominate the fungal communities in terms of the number of sequencing reads and OTUs. These unicellular forms are commonly found in aquatic environments, and their hyperdiversity may be the result of frequent inundation combined with numerous aquatic microenvironments in these rainforests. Other fungi that are frequent in aquatic environments, such as the abundant Chytridiomycotina, were also detected. While there was low similarity in OTU composition within and between the three rainforests, the fungal communities in Central America were more similar to each other than the communities in South America, reflecting a general biogeographic pattern also seen in animals, plants and protists.

Ciliate diversity and distribution patterns in the sediments of a seamount and adjacent abyssal plains in the tropical Western Pacific Ocean

Background

Benthic ciliates and the environmental factors shaping their distribution are far from being completely understood. Likewise, deep-sea systems are amongst the least understood ecosystems on Earth. In this study, using high-throughput DNA sequencing, we investigated the diversity and community composition of benthic ciliates in different sediment layers of a seamount and an adjacent abyssal plain in the tropical Western Pacific Ocean with water depths ranging between 813 m and 4566 m. Statistical analyses were used to assess shifts in ciliate communities across vertical sediment gradients and water depth.

Results

Nine out of 12 ciliate classes were detected in the different sediment samples, with Litostomatea accounting for the most diverse group, followed by Plagiopylea and Oligohymenophorea. The novelty of ciliate genetic diversity was extremely high, with a mean similarity of 93.25% to previously described sequences. On a sediment depth gradient, ciliate community structure was more similar within the upper sediment layers (0-1 and 9-10 cm) compared to the lower sediment layers (19-20 and 29-30 cm) at each site. Some unknown ciliate taxa which were absent from the surface sediments were found in deeper sediments layers. On a water depth gradient, the proportion of unique OTUs was between 42.2% and 54.3%, and that of OTUs shared by all sites around 14%. However, alpha diversity of the different ciliate communities was relatively stable in the surface layers along the water depth gradient, and about 78% of the ciliate OTUs retrieved from the surface layer of the shallowest site were shared with the surface layers of sites deeper than 3800 m. Correlation analyses did not reveal any significant effects of measured environmental factors on ciliate community composition and structure.

Conclusions

We revealed an obvious variation in ciliate community along a sediment depth gradient in the seamount and the adjacent abyssal plain and showed that water depth is a less important factor shaping ciliate distribution in deep-sea sediments unlike observed for benthic ciliates in shallow seafloors. Additionally, an extremely high genetic novelty of ciliate diversity was found in these habitats, which points to a hot spot for the discovery of new ciliate species.

Electronic supplementary material

The online version of this article (10.1186/s12866-017-1103-6) contains supplementary material, which is available to authorized users.

Strengths and Biases of High-Throughput Sequencing Data in the Characterization of Freshwater Ciliate Microbiomes

Molecular surveys of eukaryotic microbial communities employing high-throughput sequencing (HTS) techniques are rapidly supplanting traditional morphological approaches due to their larger data output and reduced bench work time. Here, we directly compare morphological and Illumina data obtained from the same samples, in an effort to characterize ciliate faunas from sediments in freshwater environments. We show how in silico processing affects the final outcome of our HTS analysis, providing evidence that quality filtering protocols strongly impact the number of predicted taxa, but not downstream conclusions such as biogeography patterns. We determine the abundance distribution of ciliates, showing that a small fraction of abundant taxa dominates read counts. At the same time, we advance reasons to believe that biases affecting HTS abundances may be significant enough to blur part of the underlying biological picture. We confirmed that the HTS approach detects many more taxa than morphological inspections, and highlight how the difference varies among taxonomic groups. Finally, we hypothesize that the two datasets actually correspond to different conceptions of "diversity," and consequently that neither is entirely superior to the other when investigating environmental protists.

The Protistan Microbiome of Grassland Soil: Diversity in the Mesoscale

Genomic data for less than one quarter of ∼1.8 million named species on earth exist in public databases like GenBank. Little information exists on the estimated one million small sized (1-100μm) heterotrophic nanoflagellates and ciliates and their taxa-area relationship. We analyzed environmental DNA from 150 geo-referenced grassland plots representing topographical and land-use ranges typical for Central Europe. High through-put barcoding allowed the identification of operational taxonomic units (OTUs) at species level, with high pairwise identity to reference sequences (≥99.7%), but also the identification of sequences at the genus (≥97%) and class (≥80%) taxonomic level. Species richness analyses revealed, on average, 100 genus level OTUs (332 unique individual read (UIR) and 56 class level OTUs per gram of soil sample in the mesoscale (1-1000km). Database shortfalls were highlighted by increased uncertain taxonomic lineages at lower resolution (≥80% sequence identity). No single barcode occurred ubiquitously across all sites. Taxa-area relationships indicated that OTUs spread over the entire mesoscale were more similar than in the local scale and increased land-use (fertilization, mowing and grazing) promoted taxa-area separation. Only a small fraction of sequences strictly matched reference library sequences, suggesting a large protistan "dark matter" in soil which warrants further research.

Distribution of Ciliates in Intertidal Sediments across Geographic Distances: A Molecular View

The estimation of the diversity and geographic distribution of protists in particular ciliates has long been an ongoing debate. We estimated the distribution of ciliates in intertidal sediments with geographic distance ranging from 10cm to 1,000km, using high-throughput DNA sequencing. Higher community similarity, in particular that of rare ciliate taxa, occurred between replicates than between sites. Statistical analyses showed that the community dissimilarities of both abundant and moderately abundant ciliates had correlations with geographic distance at the scales of 50m, 200km and 1000km, indicating a distance-decay relationship. No significant correlation was observed for rare taxa. The analysis of Bray-Curtis similarity showed a higher pairwise community similarity in abundant taxa than in moderately abundant taxa, and that in rare taxa was the lowest. Abundant taxa usually can disperse in a wider range than rare taxa, though their distribution is restricted by distance to some extent. Rare taxa are potentially more sensitive to changing environments, but no distance-decay relationship could be observed. The data indicate contrasting patterns of geographic distribution of the abundant and rare ciliate taxa and a weak distance-decay relationship for relatively abundant taxa at a scale over 50m.

Phylogenetically Structured Differences in rRNA Gene Sequence Variation among Species of Arbuscular Mycorrhizal Fungi and Their Implications for Sequence Clustering

Arbuscular mycorrhizal (AM) fungi form mutualisms with plant roots that increase plant growth and shape plant communities. Each AM fungal cell contains a large amount of genetic diversity, but it is unclear if this diversity varies across evolutionary lineages. We found that sequence variation in the nuclear large-subunit (LSU) rRNA gene from 29 isolates representing 21 AM fungal species generally assorted into genus- and species-level clades, with the exception of species of the genera Claroideoglomus and Entrophospora However, there were significant differences in the levels of sequence variation across the phylogeny and between genera, indicating that it is an evolutionarily constrained trait in AM fungi. These consistent patterns of sequence variation across both phylogenetic and taxonomic groups pose challenges to interpreting operational taxonomic units (OTUs) as approximations of species-level groups of AM fungi. We demonstrate that the OTUs produced by five sequence clustering methods using 97% or equivalent sequence similarity thresholds failed to match the expected species of AM fungi, although OTUs from AbundantOTU, CD-HIT-OTU, and CROP corresponded better to species than did OTUs from mothur or UPARSE. This lack of OTU-to-species correspondence resulted both from sequences of one species being split into multiple OTUs and from sequences of multiple species being lumped into the same OTU. The OTU richness therefore will not reliably correspond to the AM fungal species richness in environmental samples. Conservatively, this error can overestimate species richness by 4-fold or underestimate richness by one-half, and the direction of this error will depend on the genera represented in the sample.

IMPORTANCE

Arbuscular mycorrhizal (AM) fungi form important mutualisms with the roots of most plant species. Individual AM fungi are genetically diverse, but it is unclear whether the level of this diversity differs among evolutionary lineages. We found that the amount of sequence variation in an rRNA gene that is commonly used to identify AM fungal species varied significantly between evolutionary groups that correspond to different genera, with the exception of two genera that are genetically indistinguishable from each other. When we clustered groups of similar sequences into operational taxonomic units (OTUs) using five different clustering methods, these patterns of sequence variation caused the number of OTUs to either over- or underestimate the actual number of AM fungal species, depending on the genus. Our results indicate that OTU-based inferences about AM fungal species composition from environmental sequences can be improved if they take these taxonomically structured patterns of sequence variation into account.

Comparison of three clustering approaches for detecting novel environmental microbial diversity

Discovery of novel diversity in high-throughput sequencing studies is an important aspect in environmental microbial ecology. To evaluate the effects that amplicon clustering methods have on the discovery of novel diversity, we clustered an environmental marine high-throughput sequencing dataset of protist amplicons together with reference sequences from the taxonomically curated Protist Ribosomal Reference (PR²) database using three de novo approaches: sequence similarity networks, USEARCH, and Swarm. The potentially novel diversity uncovered by each clustering approach differed drastically in the number of operational taxonomic units (OTUs) and in the number of environmental amplicons in these novel diversity OTUs. Global pairwise alignment comparisons revealed that numerous amplicons classified as potentially novel by USEARCH and Swarm were more than 97% similar to references of PR². Using shortest path analyses on sequence similarity network OTUs and Swarm OTUs we found additional novel diversity within OTUs that would have gone unnoticed without further exploiting their underlying network topologies. These results demonstrate that graph theory provides powerful tools for microbial ecology and the analysis of environmental high-throughput sequencing datasets. Furthermore, sequence similarity networks were most accurate in delineating novel diversity from previously discovered diversity.

Morphological Description and Molecular Phylogeny of Two Species of Levicoleps (Ciliophora, Prostomatida), L. taehwae nov. spec. and L. biwae jejuensis nov. subspec., collected in Korea

Two colepid ciliates, Levicoleps taehwae nov. spec. and L. biwae jejuensis nov. subspec., were collected from the brackish water of the Taehwa River and a small freshwater pond in Jeju Island, South Korea, respectively. Their living morphology, infraciliature, and small subunit (SSU) rRNA gene sequences were determined using standard methods. Barrel-shaped L. taehwae nov. spec. is a small ciliate with an average size of 45 × 25 μm in vivo, about 15 ciliary rows each composed of 12 monokinetids and two perioral dikinetids, and two 20 μm-long caudal cilia. The sequence length and GC content of the SSU rRNA gene are 1,669 bp, 44.5%. This novel species is similar in body size to Coleps hirtus, and has six armor tiers and hirtus-type tier plates, and the same number of ciliary rows as C. hirtus; however, it can be distinguished from the latter by the absence of armor spines and its sequence similarity of SSU rRNA gene is about 92.8% which indicates that it is a distinct form. Levicoleps biwae jejuensis nov. subspec., is a medium colepid ciliate which has a barrel-shaped body, about 22 somatic kineties and 16 transverse ciliary rows, three mini adoral organelles, and four 15 μm-long caudal cilia. The sequence length and GC content of the SSU rRNA gene are 1,666 bp and 44.4%.

Swarm v2: highly-scalable and high-resolution amplicon clustering

Previously we presented Swarm v1, a novel and open source amplicon clustering program that produced fine-scale molecular operational taxonomic units (OTUs), free of arbitrary global clustering thresholds and input-order dependency. Swarm v1 worked with an initial phase that used iterative single-linkage with a local clustering threshold (d), followed by a phase that used the internal abundance structures of clusters to break chained OTUs. Here we present Swarm v2, which has two important novel features: (1) a new algorithm for d = 1 that allows the computation time of the program to scale linearly with increasing amounts of data; and (2) the new fastidious option that reduces under-grouping by grafting low abundant OTUs (e.g., singletons and doubletons) onto larger ones. Swarm v2 also directly integrates the clustering and breaking phases, dereplicates sequencing reads with d = 0, outputs OTU representatives in fasta format, and plots individual OTUs as two-dimensional networks.

New Primers Targeting Full-Length Ciliate 18S rRNA Genes and Evaluation of Dietary Effect on Rumen Ciliate Diversity in Dairy Cows

Analysis of the full-length 18S rRNA gene sequences of rumen ciliates is more reliable for taxonomical classification and diversity assessment than the analysis of partial hypervariable regions only. The objective of this study was to develop new oligonucleotide primers targeting the full-length 18S rRNA genes of rumen ciliates, and to evaluate the effect of different sources of dietary fiber (corn stover or a mixture of alfalfa hay and corn silage) and protein (mixed rapeseed, cottonseed, and/or soybean meals) on rumen ciliate diversity in dairy cows. Primers were designed based on a total of 137 previously reported ciliate 18S rRNA gene sequences. The 3'-terminal sequences of the newly designed primers, P.1747r_2, P.324f, and P.1651r, demonstrated >99% base coverage. Primer pair D (P.324f and P.1747r_2) was selected for the cloning and sequencing of ciliate 18S rRNA genes because it produced a 1423-bp amplicon, and did not amply the sequences of other eukaryotic species, such as yeast. The optimal species-level cutoff value for distinguishing between the operational taxonomic units of different ciliate species was 0.015. The phylogenetic analysis of full-length ciliate 18S rRNA gene sequences showed that distinct ciliate profiles were induced by the different sources of dietary fiber and protein. Dasytricha and Entodinium were the predominant genera in the ruminal fluid of dairy cattle, and Dasytricha was significantly more abundant in cows fed with corn stover than in cows fed with alfalfa hay and corn silage.

Mining environmental high-throughput sequence data sets to identify divergent amplicon clusters for phylogenetic reconstruction and morphotype visualization

Environmental high-throughput sequencing (envHTS) is a very powerful tool, which in protistan ecology is predominantly used for the exploration of diversity and its geographic and local patterns. We here used a pyrosequenced V4-SSU rDNA data set from a solar saltern pond as test case to exploit such massive protistan amplicon data sets beyond this descriptive purpose. Therefore, we combined a Swarm-based blastn network including 11 579 ciliate V4 amplicons to identify divergent amplicon clusters with targeted polymerase chain reaction (PCR) primer design for full-length small subunit of the ribosomal DNA retrieval and probe design for fluorescence in situ hybridization (FISH). This powerful strategy allows to benefit from envHTS data sets to (i) reveal the phylogenetic position of the taxon behind divergent amplicons; (ii) improve phylogenetic resolution and evolutionary history of specific taxon groups; (iii) solidly assess an amplicons (species') degree of similarity to its closest described relative; (iv) visualize the morphotype behind a divergent amplicons cluster; (v) rapidly FISH screen many environmental samples for geographic/habitat distribution and abundances of the respective organism and (vi) to monitor the success of enrichment strategies in live samples for cultivation and isolation of the respective organisms.

Bioinformatic Amplicon Read Processing Strategies Strongly Affect Eukaryotic Diversity and the Taxonomic Composition of Communities

Amplicon read sequencing has revolutionized the field of microbial diversity studies. The technique has been developed for bacterial assemblages and has undergone rigorous testing with mock communities. However, due to the great complexity of eukaryotes and the numbers of different rDNA copies, analyzing eukaryotic diversity is more demanding than analyzing bacterial or mock communities, so studies are needed that test the methods of analyses on taxonomically diverse natural communities. In this study, we used 20 samples collected from the Baltic Sea ice, slush and under-ice water to investigate three program packages (UPARSE, mothur and QIIME) and 18 different bioinformatic strategies implemented in them. Our aim was to assess the impact of the initial steps of bioinformatic strategies on the results when analyzing natural eukaryotic communities. We found significant differences among the strategies in resulting read length, number of OTUs and estimates of diversity as well as clear differences in the taxonomic composition of communities. The differences arose mainly because of the variable number of chimeric reads that passed the pre-processing steps. Singleton removal and denoising substantially lowered the number of errors. Our study showed that the initial steps of the bioinformatic amplicon read processing strategies require careful consideration before applying them to eukaryotic communities.

High diversity of protistan plankton communities in remote high mountain lakes in the European Alps and the Himalayan mountains

We analyzed the genetic diversity (V4 region of the 18S rRNA) of planktonic microbial eukaryotes in four high mountain lakes including two remote biogeographic regions (the Himalayan mountains and the European Alps) and distinct habitat types (clear and glacier-fed turbid lakes). The recorded high genetic diversity in these lakes was far beyond of what is described from high mountain lake plankton. In total, we detected representatives from 66 families with the main taxon groups being Alveolata (55.0% OTUs 97%, operational taxonomic units), Stramenopiles (34.0% OTUs 97%), Cryptophyta (4.0% OTUs 97%), Chloroplastida (3.6% OTUs 97%) and Fungi (1.7% OTUs 97%). Centrohelida, Choanomonada, Rhizaria, Katablepharidae and Telonema were represented by <1% OTUs 97%. Himalayan lakes harbored a higher plankton diversity compared to the Alpine lakes (Shannon index). Community structures were significantly different between lake types and biogeographic regions (Fisher exact test, P < 0.01). Network analysis revealed that more families of the Chloroplastida (10 vs 5) and the Stramenopiles (14 vs 8) were found in the Himalayan lakes than in the Alpine lakes and none of the fungal families was shared between them. Biogeographic aspects as well as ecological factors such as water turbidity may structure the microbial eukaryote plankton communities in such remote lakes.

Swarm: robust and fast clustering method for amplicon-based studies

Popular de novo amplicon clustering methods suffer from two fundamental flaws: arbitrary global clustering thresholds, and input-order dependency induced by centroid selection. Swarm was developed to address these issues by first clustering nearly identical amplicons iteratively using a local threshold, and then by using clusters' internal structure and amplicon abundances to refine its results. This fast, scalable, and input-order independent approach reduces the influence of clustering parameters and produces robust operational taxonomic units.

Intracellular diversity of the V4 and V9 regions of the 18S rRNA in marine protists (radiolarians) assessed by high-throughput sequencing

Metabarcoding is a powerful tool for exploring microbial diversity in the environment, but its accurate interpretation is impeded by diverse technical (e.g. PCR and sequencing errors) and biological biases (e.g. intra-individual polymorphism) that remain poorly understood. To help interpret environmental metabarcoding datasets, we investigated the intracellular diversity of the V4 and V9 regions of the 18S rRNA gene from Acantharia and Nassellaria (radiolarians) using 454 pyrosequencing. Individual cells of radiolarians were isolated, and PCRs were performed with generalist primers to amplify the V4 and V9 regions. Different denoising procedures were employed to filter the pyrosequenced raw amplicons (Acacia, AmpliconNoise, Linkage method). For each of the six isolated cells, an average of 541 V4 and 562 V9 amplicons assigned to radiolarians were obtained, from which one numerically dominant sequence and several minor variants were found. At the 97% identity, a diversity metrics commonly used in environmental surveys, up to 5 distinct OTUs were detected in a single cell. However, most amplicons grouped within a single OTU whereas other OTUs contained very few amplicons. Different analytical methods provided evidence that most minor variants forming different OTUs correspond to PCR and sequencing artifacts. Duplicate PCR and sequencing from the same DNA extract of a single cell had only 9 to 16% of unique amplicons in common, and alignment visualization of V4 and V9 amplicons showed that most minor variants contained substitutions in highly-conserved regions. We conclude that intracellular variability of the 18S rRNA in radiolarians is very limited despite its multi-copy nature and the existence of multiple nuclei in these protists. Our study recommends some technical guidelines to conservatively discard artificial amplicons from metabarcoding datasets, and thus properly assess the diversity and richness of protists in the environment.

SSU rDNA sequence diversity and seasonally differentiated distribution of nanoplanktonic ciliates in neritic Bohai and Yellow Seas as revealed by T-RFLP

Nanociliates have been frequently found to be important players in the marine microbial loop, however, little is known about their diversity and distribution in coastal ecosystems. We investigated the molecular diversity and distribution patterns of nanoplanktonic oligotrich and choreotrich (OC) ciliates in surface water of three neritic basins of northern China, the South Yellow Sea (SYS), North Yellow Sea (NYS), and Bohai Sea (BS) in June and November 2011. SSU rRNA gene clone libraries generated from three summertime samples (sites B38, B4 and H8) were analyzed and revealed a large novel ribotype diversity, of which many were low-abundant phylotypes belonging to the subclass Oligotrichia, but divergent from described morphospecies. Based on the data of terminal-restriction fragment length polymorphism (T-RFLP) analysis of all 35 samples, we found that the T-RF richness was generally higher in the SYS than in the BS, and negatively correlated with the molar ratio of P to Si. Overall, multidimensional scaling and permutational multivariate analysis of variance of the community turnover demonstrated a distinct seasonal pattern but no basin-to-basin differentiation across all samples. Nevertheless, significant community differences among basins were recognized in the winter dataset. Mantel tests showed that the environmental factors, P:Si ratio, water temperature and concentration of dissolved oxygen (DO), determined the community across all samples. However, both biogeographic distance and environment shaped the community in winter, with DO being the most important physicochemical factor. Our results indicate that the stoichiometric ratio of P:Si is a key factor, through which the phytoplankton community may be shaped, resulting in a cascade effect on the diversity and community composition of OC nanociliates in the N-rich, Si-limited coastal surface waters, and that the Yellow Sea Warm Current drives the nanociliate community, and possibly the microbial food webs, in the coastal ecosystem in winter.

Investigating microbial eukaryotic diversity from a global census: insights from a comparison of pyrotag and full-length sequences of 18S rRNA genes

Next-generation DNA sequencing (NGS) approaches are rapidly surpassing Sanger sequencing for characterizing the diversity of natural microbial communities. Despite this rapid transition, few comparisons exist between Sanger sequences and the generally much shorter reads of NGS. Operational taxonomic units (OTUs) derived from full-length (Sanger sequencing) and pyrotag (454 sequencing of the V9 hypervariable region) sequences of 18S rRNA genes from 10 global samples were analyzed in order to compare the resulting protistan community structures and species richness. Pyrotag OTUs called at 98% sequence similarity yielded numbers of OTUs that were similar overall to those for full-length sequences when the latter were called at 97% similarity. Singleton OTUs strongly influenced estimates of species richness but not the higher-level taxonomic composition of the community. The pyrotag and full-length sequence data sets had slightly different taxonomic compositions of rhizarians, stramenopiles, cryptophytes, and haptophytes, but the two data sets had similarly high compositions of alveolates. Pyrotag-based OTUs were often derived from sequences that mapped to multiple full-length OTUs at 100% similarity. Thus, pyrotags sequenced from a single hypervariable region might not be appropriate for establishing protistan species-level OTUs. However, nonmetric multidimensional scaling plots constructed with the two data sets yielded similar clusters, indicating that beta diversity analysis results were similar for the Sanger and NGS sequences. Short pyrotag sequences can provide holistic assessments of protistan communities, although care must be taken in interpreting the results. The longer reads (>500 bp) that are now becoming available through NGS should provide powerful tools for assessing the diversity of microbial eukaryotic assemblages.

Molecular diversity reveals previously undetected air-dispersed protist colonists in a Mediterranean area

The molecular diversity of air-dispersed protists was examined through the 18S rRNA gene clone library construction in air samples and samples from experimental water containers passively collecting air-dispersed microorganisms, from July 2007 till October 2008 in three different sites of Northern Greece. The majority of the samplings took place in an urban industrialized coastal city (Thessaloniki). In all the samples, a total of 29 unique phylotypes were detected belonging to 10 known major taxonomic groups. The most abundant phylotypes were affiliated to known taxa of Ciliophora and Chlorophyceae, commonly found in various habitats. Additionally, various previously unnoticed and under-studied taxa, such as Bicosoecida, Oomycetes and Labyrinthulomycetes, were detected. These taxa are potentially important in ecological processes, through dispersal and colonization of various habitats. Multivariate statistical analysis associated the most abundant phylotypes with rainfall, suggesting that rain is a favorable means for reposition of air-dispersed protists. This is the first study investigating the molecular diversity of air-dispersed protists, including algae and heterotrophic protists.

Placing environmental next-generation sequencing amplicons from microbial eukaryotes into a phylogenetic context

Nucleotide positions in the hypervariable V4 and V9 regions of the small subunit (SSU)-rDNA locus are normally difficult to align and are usually removed before standard phylogenetic analyses. Yet, with next-generation sequencing data, amplicons of these regions are all that are available to answer ecological and evolutionary questions that rely on phylogenetic inferences. With ciliates, we asked how inclusion of the V4 or V9 regions, regardless of alignment quality, affects tree topologies using distinct phylogenetic methods (including PairDist that is introduced here). Results show that the best approach is to place V4 amplicons into an alignment of full-length Sanger SSU-rDNA sequences and to infer the phylogenetic tree with RAxML. A sliding window algorithm as implemented in RAxML shows, though, that not all nucleotide positions in the V4 region are better than V9 at inferring the ciliate tree. With this approach and an ancestral-state reconstruction, we use V4 amplicons from European nearshore sampling sites to infer that rather than being primarily terrestrial and freshwater, colpodean ciliates may have repeatedly transitioned from terrestrial/freshwater to marine environments.

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.

General patterns of diversity in major marine microeukaryote lineages

Microeukaryotes have vital roles for the functioning of marine ecosystems, but still some general characteristics of their current diversity and phylogeny remain unclear. Here we investigated both aspects in major oceanic microeukaryote lineages using 18S rDNA (V4–V5 hypervariable regions) sequences from public databases that derive from various marine environmental surveys. A very carefully and manually curated dataset of 8291 Sanger sequences was generated and subsequently split into 65 taxonomic groups (roughly to Class level based on KeyDNATools) prior to downstream analyses. First, we calculated genetic distances and clustered sequences into Operational Taxonomic Units (OTUs) using different distance cut-off levels. We found that most taxonomic groups had a maximum pairwise genetic distance of 0.25. Second, we used phylogenetic trees to study general evolutionary patterns. These trees confirmed our taxonomic classification and served to run Lineage Through Time (LTT) plots. LTT results indicated different cladogenesis dynamics across groups, with some displaying an early diversification and others a more recent one. Overall, our study provides an improved description of the microeukaryote diversity in the oceans in terms of genetic differentiation within groups as well as in the general phylogenetic structure. These results will be important to interpret the large amount of sequence data that is currently generated by High Throughput Sequencing technologies.

Phylogenetic investigations on ten genera of tintinnid ciliates (Ciliophora: Spirotrichea: Tintinnida), based on small subunit ribosomal RNA gene sequences

Tintinnida is a diverse taxon that accommodates over 1,500 morphospecies, which is an important component of marine planktonic food webs. However, evolutionary relationships of tintinnids are poorly known because molecular data of most groups within this order are lacking. In our study, the small subunit (SSU) rRNA genes representing 10 genera, 5 families of Tintinnida were sequenced, including the first SSU rRNA gene sequences for Coxliella, Dadayiella, Epiplocyloides, and Protorhabdonella, and phylogenetic trees were constructed to assess their intergeneric relationships. Phylogenies inferred from different methods showed that (1) Three newly sequenced Eutintinnus species fell into Eutintinnus clade forming a sister group to the clade containing Amphorides, Steenstrupiella, Amphorellopsis, and Salpingella; (2) Surprisingly, the genetic distances between Amphorides amphora and Amphorellopsis acuta population 1 was even smaller than that between the two populations of Amphorellopsis acuta, casting doubt on the validity of Amphorides and Amphorellopsis as presently defined; (3) The SSU rRNA sequences of Dadayiella ganymedes and Parundella aculeata were almost identical. Therefore, Parundella ganymedes novel combination is proposed; (4) Coxliella, which is currently assigned within Metacylididae, branched instead with some Tintinnopsis species. Furthermore, the validation of Coxliella, which was considered to be a "questionable" genus, was confirmed based on evidences from morphology, ecology, and molecular data; (5) Protorhabdonella and Rhabdonella showed rather low intergeneric distance and grouped together with strong support suggesting that Rhabdonellidae is a well-defined taxon; and (6) Epiplocyloides branched with species in Cyttarocylididae indicating their close relationship.

Morphological, ultrastructural, and molecular characterization of Euplotidium rosati n. sp. (Ciliophora, Euplotida) from Guam

We combined morphological (i.e. live, stained, scanning, and transmission electron microscopy) with morphometric and molecular analysis to describe a ciliate species collected from shallow reefs in Guam, grown, and maintained in our laboratory. The species was recognized as a member of Euplotidium, and compared with established species of the genus: Euplotidium itoi Ito 1958; Euplotidium psammophilus (Vacelet 1961) Borror 1972; Euplotidium arenarium Magagnini and Nobili 1964; Euplotidium helgae Hartwig 1980; Euplotidium prosaltans Tuffrau 1985, and Euplotidium smalli Lei, Choi and Xu, 2002. To obtain more elements to compare the species, new morphometric data and additional SSU rRNA gene sequences of E. itoi and of E. arenarium are reported. On the basis of this comparison, we established the new species Euplotidium rosati that has a cirral pattern composed of 12 frontoventral and six transverse cirri, and lacks the left marginal cirrus. Euplotidium rosati harbors on its dorsal surface epixenosomes, the peculiar extrusive symbionts described in other Euplotidium species. The whole body of our observations together with the analysis of the data available in the literature leads us to propose a redefinition of the genus. The results may also be useful to clarify the tangled relationship between Euplotidium and Gastrocirrhus.

Diversity patterns of uncultured Haptophytes unravelled by pyrosequencing in Naples Bay

Haptophytes are a key phylum of marine protists, including ~300 described morphospecies and 80 morphogenera. We used 454 pyrosequencing on large subunit ribosomal DNA (LSU rDNA) fragments to assess the diversity from size-fractioned plankton samples collected in the Bay of Naples. One group-specific primer set targeting the LSU rDNA D1/D2 region was designed to amplify Haptophyte sequences from nucleic acid extracts (total DNA or RNA) of two size fractions (0.8-3 or 3-20 μm) and two sampling depths [subsurface, at 1 m, or deep chlorophyll maximum (DCM) at 23 m]. 454 reads were identified using a database covering the entire Haptophyta diversity currently sequenced. Our data set revealed several hundreds of Haptophyte clusters. However, most of these clusters could not be linked to taxonomically known sequences: considering OTUs(97%) (clusters build at a sequence identity level of 97%) on our global data set, less than 1% of the reads clustered with sequences from cultures, and less than 12% clustered with reference sequences obtained previously from cloning and Sanger sequencing of environmental samples. Thus, we highlighted a large uncharacterized environmental genetic diversity, which clearly shows that currently cultivated species poorly reflect the actual diversity present in the natural environment. Haptophyte community appeared to be significantly structured according to the depth. The highest diversity and evenness were obtained in samples from the DCM, and samples from the large size fraction (3-20 μm) taken at the DCM shared a lower proportion of common OTUs(97%) with the other samples. Reads from the species Chrysoculter romboideus were notably found at the DCM, while they could be detected at the subsurface. The highest proportion of totally unknown OTUs(97%) was collected at the DCM in the smallest size fraction (0.8-3 μm). Overall, this study emphasized several technical and theoretical barriers inherent to the exploration of the large and largely unknown diversity of unicellular eukaryotes.

Act together-implications of symbioses in aquatic ciliates

Mutual interactions in the form of symbioses can increase the fitness of organisms and provide them with the capacity to occupy new ecological niches. The formation of obligate symbioses allows for rapid evolution of new life forms including multitrophic consortia. Microbes are important components of many known endosymbioses and their short generation times and strong potential for genetic exchange may be important drivers of speciation. Hosts provide endo- and ectosymbionts with stable, nutrient-rich environments, and protection from grazers. This is of particular importance in aquatic ecosystems, which are often highly variable, harsh, and nutrient-deficient habitats. It is therefore not surprising that symbioses are widespread in both marine and freshwater environments. Symbioses in aquatic ciliates are good model systems for exploring symbiont-host interactions. Many ciliate species are globally distributed and have been intensively studied in the context of plastid evolution. Their relatively large cell size offers an ideal habitat for numerous microorganisms with different functional traits including commensalism and parasitism. Phagocytosis facilitates the formation of symbiotic relationships, particularly since some ingested microorganisms can escape the digestion. For example, photoautotrophic algae and methanogens represent endosymbionts that greatly extend the biogeochemical functions of their hosts. Consequently, symbiotic relationships between protists and prokaryotes are widespread and often result in new ecological functions of the symbiotic communities. This enables ciliates to thrive under a wide range of environmental conditions including ultraoligotrophic or anoxic habitats. We summarize the current understanding of this exciting research topic to identify the many areas in which knowledge is lacking and to stimulate future research by providing an overview on new methodologies and by formulating a number of emerging questions in this field.

Plankton microorganisms coinciding with two consecutive mass fish kills in a newly reconstructed lake

Lake Karla, Greece, was dried up in 1962 and its refilling started in 2009. We examined the Cyanobacteria and unicellular eukaryotes found during two fish kill incidents, in March and April 2010, in order to detect possible causative agents. Both microscopic and molecular (16S/18S rRNA gene diversity) identification were applied. Potentially toxic Cyanobacteria included representatives of the Planktothrix and Anabaena groups. Known toxic eukaryotes or parasites related to fish kill events were Prymnesium parvum and Pfiesteria cf. piscicida, the latter being reported in an inland lake for the second time. Other potentially harmful microorganisms, for fish and other aquatic life, included representatives of Fungi, Mesomycetozoa, Alveolata, and Heterokontophyta (stramenopiles). In addition, Euglenophyta, Chlorophyta, and diatoms were represented by species indicative of hypertrophic conditions. The pioneers of L. Karla's plankton during the first months of its water refilling process included species that could cause the two observed fish kill events.

Marine protistan diversity

Protists have fascinated microbiologists since their discovery nearly 350 years ago. These single-celled, eukaryotic species span an incredible range of sizes, forms, and functions and, despite their generally diminutive size, constitute much of the genetic diversity within the domain Eukarya. Protists in marine ecosystems play fundamental ecological roles as primary producers, consumers, decomposers, and trophic links in aquatic food webs. Much of our knowledge regarding the diversity and ecological activities of these species has been obtained during the past half century, and only within the past few decades have hypotheses depicting the evolutionary relationships among the major clades of protists attained some degree of consensus. This recent progress is attributable to the development of genetic approaches, which have revealed an unexpectedly large diversity of protists, including cryptic species and previously undescribed clades of protists. New genetic tools now exist for identifying protistan species of interest and for reexamining long-standing debates regarding the biogeography of protists. Studies of protistan diversity provide insight regarding how species richness and community composition contribute to ecosystem function. These activities support the development of predictive models that describe how microbial communities will respond to natural or anthropogenically mediated changes in environmental conditions.