Ciliates and the Rare Biosphere: A Review

Using digital PCR to predict ciliate abundance from ribosomal RNA gene copy numbers

Ciliates play a key role in most ecosystems. Their abundance in natural samples is crucial for answering many ecological questions. Traditional methods of quantifying individual species, which rely on microscopy, are often labour-intensive, time-consuming and can be highly biassed. As a result, we investigated the potential of digital polymerase chain reaction (dPCR) for quantifying ciliates. A significant challenge in this process is the high variation in the copy number of the taxonomic marker gene (ribosomal RNA [rRNA]). We first quantified the rRNA gene copy numbers (GCN) of the model ciliate, Paramecium tetraurelia, during different stages of the cell cycle and growth phases. The per-cell rRNA GCN varied between approximately 11,000 and 130,000, averaging around 50,000 copies per cell. Despite these variations in per-cell rRNA GCN, we found a highly significant correlation between GCN and cell numbers. This is likely due to the coexistence of different cellular stages in an uncontrolled (environmental) ciliate population. Thanks to the high sensitivity of dPCR, we were able to detect the target gene in a sample that contained only a single cell. The dPCR approach presented here is a valuable addition to the molecular toolbox in protistan ecology. It may guide future studies in quantifying and monitoring the abundance of targeted (even rare) ciliates in natural samples.

Insights into freshwater ciliate diversity through high throughput DNA metabarcoding

The freshwater bodies of India are highly biodiverse but still understudied, especially concerning ciliates. Ciliates constitute a significant portion of eukaryotic diversity and play crucial roles in microbial loops, nutrient recycling, and ecosystem maintenance. The present study aimed to elucidate ciliate diversity in three freshwater sites in the Delhi region of India: Okhla Bird Sanctuary (OBS), Sanjay Lake (SL), and Raj Ghat pond (RJ). This study represents the first investigation into the taxonomic diversity and richness of freshwater ciliates in India using a high-throughput DNA metabarcoding approach. For the analysis, total environmental DNA was extracted from the three freshwater samples, followed by sequencing of the 18S V4 barcode region and subsequent phylogenetic analyses. Operational taxonomic units (OTU) analyses revealed maximum species diversity in OBS (106), followed by SL (104) and RJ (99) sites. Ciliates from the classes Oligohymenophorea, Prostomatea, and Spirotrichea were dominant in the three sites. The study discusses the ability of the metabarcoding approach to uncover unknown and rare species. The study highlights the need for refined reference databases and cautious interpretation of the high-throughput sequencing-generated data while emphasizing the complementary nature of molecular and morphological approaches in studying ciliate diversity.

New contributions to the Cyrtophoria ciliates (Protista, Ciliophora): Establishment of new taxa and phylogenetic analyses using two ribosomal genes

Periphytic ciliates play a vital role in the material cycle and energy flow of microbial food web, however, their taxonomy and biodiversity are inadequately studied given their high species richness. Two new and one little known species, viz. Derouxella lembodes gen. et sp. nov., Cyrtophoron multivacuolatum sp. nov., and Cyrtophoron apsheronica Aliev, 1991, collected from coastal waters of China, were investigated using modern methods. Derouxella gen. nov. can be recognized by having dorsoventrally flattened body, a podite, one fragmented preoral kinety, two parallel circumoral kineties, and somatic kineties progressively shortened from right to left. Morphological classification and phylogenetic analyses based on nuclear small subunit ribosomal RNA (nSSU rRNA) and mitochondrial small subunit ribosomal RNA (mtSSU rRNA) gene sequence data inferred that Derouxella gen. nov. occupies an intermediate position between Hartmannulidae and Dysteriidae. Cyrtophoron multivacuolatum sp. nov. is characterized by large body size, the numbers of somatic kineties and nematodesmal rods, and having numerous contractile vacuoles. The genus Cyrtophoron and the poorly known species C. apsheronica were redefined. Even with the addition of newly obtained nSSU rRNA and mtSSU rRNA gene sequences of Cyrtophoron, the family Chlamydodontidae was still recovered as a monophyletic group, the monophyly of Cyrtophoron was supported too.

The Diversity Patterns of Rare to Abundant Microbial Eukaryotes Across a Broad Range of Salinities in a Solar Saltern

Solar salterns are excellent artificial systems for examining species diversity and succession along salinity gradients. Here, the eukaryotic community in surface water of a Korean solar saltern (30 to 380 practical salinity units) was investigated from April 2019 to October 2020 using Illumina sequencing targeting the V4 and V9 regions of 18S rDNA. A total of 926 operational taxonomic units (OTUs) and 1,999 OTUs were obtained with the V4 and V9 regions, respectively. Notably, most of the OTUs were microbial eukaryotes, and the high-abundance groups (> 5% relative abundance (RA), Alveolata, Stramenopila, Archaeplastida, and Opisthokonta) usually accounted for > 90% of the total cumulative read counts and > 80% of all OTUs. Moreover, the high-abundance Alveolata (larger forms) and Stramenopila (smaller forms) groups displayed a significant inverse relationship, probably due to predator-prey interactions. Most of the low-abundance (0.1-5% RA) and rare (< 0.1% RA) groups remained small portion during the field surveys. Taxonomic novelty (at < 90% sequence identity) was high in the Amoebozoa, Cryptista, Haptista, Rhizaria, and Stramenopila groups (69.8% of all novel OTUs), suggesting the presence of a large number of hidden species in hypersaline environments. Remarkably, the high-abundance groups had little overlap with the other groups, implying the weakness of rare-to-prevalent community dynamics. The low-abundance Discoba group alone temporarily became the high-abundance group, suggesting that it is an opportunistic group. Overall, the composition and diversity of the eukaryotic community in hypersaline environments may be persistently stabilized, despite diverse disturbance events.

Ecology of planktonic ciliates in a changing world: Concepts, methods, and challenges

Plankton ecologists ultimately focus on forecasting, both applied and environmental outcomes. We review how appreciating planktonic ciliates has become central to these predictions. We explore the 350-year-old canon on planktonic ciliates and examine its steady progression, which has been punctuated by conceptual insights and technological breakthroughs. By reflecting on this process, we offer suggestions as to where future leaps are needed, with an emphasis on predicting outcomes of global warming. We conclude that in terms of climate change research: (i) climatic hotspots (e.g. polar oceans) require attention; (ii) simply adding ciliate measurements to zooplankton/phytoplankton-based sampling programs is inappropriate; (iii) elucidating the rare biosphere's functional ecology requires culture-independent genetic methods; (iv) evaluating genetic adaptation (microevolution) and population composition shifts is required; (v) contrasting marine and freshwaters needs attention; (vi) mixotrophy needs attention; (vii) laboratory and field studies must couple automated measurements and molecular assessment of functional gene expression; (viii) ciliate trophic diversity requires appreciation; and (ix) marrying gene expression and function, coupled with climate change scenarios is needed. In short, continued academic efforts and financial support are essential to achieve the above; these will lead to understanding how ciliates will respond to climate change, providing tools for forecasting.

A Case Study of the Morphological and Molecular Variation within a Ciliate Genus: Taxonomic Descriptions of Three Dysteria Species (Ciliophora, Cyrtophoria), with the Establishment of a New Species

Three Dysteria species, D. crassipes Claparède & Lachmann, 1859; D. brasiliensis Faria et al., 1922; and D. paracrassipes n. sp., were collected from subtropical coastal waters of the East China Sea, near Ningbo, China. The three species were studied based on their living morphology, infraciliature, and molecular data. The new species D. paracrassipes n. sp. is very similar to D. crassipes in most morphological features except the preoral kinety, which is double-rowed in the new species (vs. single-rowed in D. crassipes). The difference in the small ribosomal subunit sequences (SSU rDNA) between these two species is 56 bases, supporting the establishment of the new species. The Ningbo population of D. crassipes is highly similar in morphology to other known populations. Nevertheless, the SSU rDNA sequences of these populations are very different, indicating high genetic diversity and potentially cryptic species. Dysteria brasiliensis is cosmopolitan with many described populations worldwide and four deposited SSU rDNA sequences. The present work supplies morphological and molecular information from five subtropical populations of D. brasiliensis that bear identical molecular sequences but show significant morphological differences. The findings of this study provide an opportunity to improve understanding of the morphological and genetic diversity of ciliates.

High variability in SSU rDNA gene copy number among planktonic foraminifera revealed by single-cell qPCR

Metabarcoding has become the workhorse of community ecology. Sequencing a taxonomically informative DNA fragment from environmental samples gives fast access to community composition across taxonomic groups, but it relies on the assumption that the number of sequences for each taxon correlates with its abundance in the sampled community. However, gene copy number varies among and within taxa, and the extent of this variability must therefore be considered when interpreting community composition data derived from environmental sequencing. Here we measured with single-cell qPCR the SSU rDNA gene copy number of 139 specimens of five species of planktonic foraminifera. We found that the average gene copy number varied between of ~4000 to ~50,000 gene copies between species, and individuals of the same species can carry between ~300 to more than 350,000 gene copies. This variability cannot be explained by differences in cell size and considering all plausible sources of bias, we conclude that this variability likely reflects dynamic genomic processes acting during the life cycle. We used the observed variability to model its impact on metabarcoding and found that the application of a correcting factor at species level may correct the derived relative abundances, provided sufficiently large populations have been sampled.

Taxonomy and Phylogeny of Three Species of Dysteria (Ciliophora, Phyllopharyngea) Including the Description of Dysteria ozakii nom. nov

In the last two decades, cyrtophorian ciliates have been revealed to demonstrate a high species diversity. But this group remains difficult to study, mainly because of their low abundance and relatively few taxonomically informative morphological characters. As a contribution to the taxonomy of cyrtophorians, here we investigate three Dysteria species based on their live morphology, ciliary pattern, and molecular phylogeny. Dysteria ozakii nom. nov. can be recognized by its elongate body shape and four right kineties, including three frontoventral kineties. A neotype has to be fixed for the species as no type materials were deposited. The other two species, D. brasiliensis Faria et al., 1922 and D. compressa (Gourret & Roeser, 1886) Kahl, 1931, are redescribed and supplementary information for each is supplied. Phylogenetic analyses based on small-subunit (SSU) rRNA gene sequences support the validity of the species. In addition, four species of bacterial epibionts were observed on the surface of the three Dysteria spp. The identities of these bacterial species are discussed based on the newly obtained 16S rRNA gene sequences.

Ciliate SSU-rDNA reference alignments and trees for phylogenetic placements of metabarcoding data

Although ciliates are one of the most dominant microbial eukaryotic groups in many environments, there is a lack of updated global ciliate alignments and reference trees that can be used for phylogenetic placement methods to analyze environmental metabarcoding data. Here we fill this gap by providing reference alignments and trees for those ciliates taxa with available SSU-rDNA sequences derived from identified species. Each alignment contains 478 ciliate and six outgroup taxa, and they were made using different masking strategies for alignment positions (unmasked, masked and masked except the hypervariable V4 region). We constrained the monophyly of the major ciliate groups based on the recently updated classification of protists and based on phylogenomic data. Taxa of uncertain phylogenetic position were kept unconstrained, except for Mesodinium species that we constrained to form a clade with the Litostomatea. These ciliate reference alignments and trees can be used to perform taxonomic assignments of metabarcoding data, discover novel ciliate clades, estimate species richness, and overlay measured ecological parameters onto the phylogenetic placements.

Protists in the Plant Microbiome: An Untapped Field of Research

Protists are mostly unicellular eukaryotes. Some protists are beneficial for plants, while others live as endosymbionts and can cause severe plant diseases. More detailed studies on plant-protist interactions exist only for plant pathogens and parasites. A number of protists live as inconspicuous endophytes and cause no visible disease symptoms, while others appear closely associated with the rhizosphere or phyllosphere of plants, but we still have only a vague understanding on their identities and functions. Here, we provide a protocol on how to assess the plant-associated protist community via Illumina-sequencing of ribosomal marker-amplicons and describe how to assign taxonomic affiliation to the obtained sequences.

Perspectives from Ten Years of Protist Studies by High-Throughput Metabarcoding

During the last decade, high-throughput metabarcoding became routine for analyzing protistan diversity and distributions in nature. Amid a multitude of exciting findings, scientists have also identified and addressed technical and biological limitations, although problems still exist for inference of meaningful taxonomic and ecological knowledge based on short DNA sequences. Given the extensive use of this approach, it is critical to settle our understanding on its strengths and weaknesses and to synthesize up-to-date methodological and conceptual trends. This article summarizes key scientific and technical findings, and identifies current and future directions in protist research that uses metabarcoding.

Rain-Fed Granite Rock Basins Accumulate a High Diversity of Dormant Microbial Eukaryotes

Rain fed granite rock basins are ancient geological landforms of worldwide distribution and structural simplicity. They support habitats that can switch quickly from terrestrial to aquatic along the year. Diversity of animals and plants, and the connexion between communities in different basins have been widely explored in these habitats, but hardly any research has been carried out on microorganisms. The aim of this study is to provide the first insights on the diversity of eukaryotic microbial communities from these environments. Due to the ephemeral nature of these aquatic environments, we predict that the granitic basins should host a high proportion of dormant microeukaryotes. Based on an environmental DNA diversity survey, we reveal diverse communities with representatives of all major eukaryotic taxonomic supergroups, mainly composed of a diverse pool of low abundance OTUs. Basin communities were very distinctive, with alpha and beta diversity patterns non-related to basin size or spatial distance respectively. Dissimilarity between basins was mainly characterised by turnover of OTUs. The strong microbial eukaryotic heterogeneity observed among the basins may be explained by a complex combination of deterministic factors (diverging environment in the basins), spatial constraints, and randomness including founder effects. Most interestingly, communities contain organisms that cannot coexist at the same time because of incompatible metabolic requirements, thus suggesting the existence of a pool of dormant organisms whose activity varies along with the changing environment. These organisms accumulate in the pools, which turns granitic rock into high biodiversity microbial islands whose conservation and study deserve further attention.

Bacterial and ciliate biofilm community structure at different spatial levels of a salt lake meta-community

Meta-communities are assembled along an ecological scale that determines local and regional diversity. Spatial patterns have been detected in planktonic bacterial communities at distances <20 m, but little is known about the occurrence of similar variation for other microbial groups and changes in microbial meta-community assembly at different levels of a meta-community. To examine this variation, the biofilm of eight saline ponds were used to investigate processes shaping diversity within ponds (β) and between ponds (δ). Bacterial and ciliate communities were assessed using ARISA and T-RFLP respectively, while diversity partitioning methods were used to examine the importance of taxonomic turnover and variation partitioning was used to distinguish spatial from environmental determinants. The results show that turnover is important for determining β- and δ-diversity of biofilms. Spatial factors are important drivers of bacterial β-diversity but were unimportant for ciliate β-diversity. Environmental variation was a strong determinant of bacterial and ciliate δ-diversity, suggesting sorting processes are important for assembling pond communities. Determinants of diversity in bacteria are not universal for ciliates, suggesting higher functional redundancy of bacteria or the greater niche breadth of ciliates may be important in discriminating assembly processes between the two organisms.

Techniques and tools for species identification in ciliates: a review

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

Protist Biodiversity and Biogeography in Lakes From Four Brazilian River-Floodplain Systems

The biodiversity and biogeography of protists inhabiting many ecosystems have been intensely studied using different sequencing approaches, but tropical ecosystems are relatively under-studied. Here, we sampled planktonic waters from 32 lakes associated with four different river-floodplains systems in Brazil, and sequenced the DNA using a metabarcoding approach with general eukaryotic primers. The lakes were dominated by the largely free-living Discoba (mostly the Euglenida), Ciliophora, and Ochrophyta. There was low community similarity between lakes even within the same river-floodplain. The protists inhabiting these floodplain systems comprise part of the large and relatively undiscovered diversity in the tropics.

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.

Common and Rare Taxa of Planktonic Ciliates: Influence of Flood Events and Biogeographic Patterns in Neotropical Floodplains

After much discussion about the cosmopolitan nature of microbes, the great issue nowadays is to identify at which spatial extent microorganisms may display biogeographic patterns and if temporal variation is important in altering those patterns. Here, planktonic ciliates were sampled from shallow lakes of four Neotropical floodplains, distributed over a spatial extent of ca. 3000 km, during high and low water periods, along with several abiotic and biotic variables potentially affecting the ciliate community. We found that common ciliate species were more associated with environmental gradients and rare species were more related to spatial variables; however, this pattern seemed to change depending on the temporal and spatial scales considered. Environmental gradients were more important in the high waters for both common and rare species. In low waters, common species continued to be mainly driven by environmental conditions, but rare species were more associated with the spatial component, suggesting dispersal limitation likely due to differences in dispersal ability and ecological tolerance of species. We also found that common and rare species were related to different environmental variables, suggesting different ecological niches. At the largest spatial extents, rare species showed clear biogeographic patterns.

Marked seasonality and high spatial variation in estuarine ciliates are driven by exchanges between the 'abundant' and 'intermediate' biospheres

We examined the spatial and temporal variability of ciliate community in a subtropical estuary by rRNA and rDNA-based high throughput sequencing of 97 samples collected along the entire salinity gradient at two-month intervals in 2014. Community divided statistically into three groups: freshwater (salinity < 0.5‰), oligohaline and mesohaline (0.5‰ < salinity < 18‰), and polyhaline and euhaline (18‰ < salinity < 40‰). Across all three groups, salinity explained most of the community variability. Within each group, seasonal shifts in community formed cool (spring and winter) and warm (summer and autumn) subgroups, indicating that spatial variability overrode seasonal changes in determining community composition. Cool and warm groups showed opposite associations with temperature and prey proxies, suggesting distinct seasonal niche separation. The community reassembly of cool and warm groups was essentially due to transitions between intermediate (with relative abundance of 0.01–1%) and abundant (with relative abundance > 1%) OTUs. Further analyses demonstrated that the intermediate group not only encompassed comparable OTU richness to that of the total community and maintained high metabolic activity but also had the highest proportion in transition, either to abundance or rarity, thus offering a first view on how it varies across space and time and revealing the essential role it played in maintaining stability and functionality within the community.

Transition boundaries for protistan species turnover in hypersaline waters of different biogeographic regions

The identification of environmental barriers which govern species distribution is a fundamental concern in ecology. Even though salt was previously identified as a major transition boundary for micro- and macroorganisms alike, the salinities causing species turnover in protistan communities are unknown. We investigated 4.5 million high-quality protistan metabarcodes (V4 region of the SSU rDNA) obtained from 24 shallow salt ponds (salinities 4%-44%) from South America and Europe. Statistical analyses of protistan community profiles identified four salinity classes, which strongly selected for different protistan communities: 4-9%, 14-24%, 27-36% and 38-44%. The proportion of organisms unknown to science is highest in the 14-24% salinity class, showing that environments within this salinity range are an unappreciated reservoir of as yet undiscovered organisms. Distinct higher-rank taxon groups dominated in the four salinity classes in terms of diversity. As increasing salinities require different cellular responses to cope with salt, our results suggest that different evolutionary lineages of protists have evolved distinct haloadaptation strategies. Salinity appears to be a stronger selection factor for the structuring of protistan communities than geography. Yet, we find a higher degree of endemism in shallow salt ponds compared with less isolated ecosystems such as the open ocean. Thus, rules for biogeographic structuring of protistan communities are not universal, but depend on the ecosystem under consideration.

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.

Benthic protists: the under-charted majority

Marine protist diversity inventories have largely focused on planktonic environments, while benthic protists have received relatively little attention. We therefore hypothesize that current diversity surveys have only skimmed the surface of protist diversity in marine sediments, which may harbor greater diversity than planktonic environments. We tested this by analyzing sequences of the hypervariable V4 18S rRNA from benthic and planktonic protist communities sampled in European coastal regions. Despite a similar number of OTUs in both realms, richness estimations indicated that we recovered at least 70% of the diversity in planktonic protist communities, but only 33% in benthic communities. There was also little overlap of OTUs between planktonic and benthic communities, as well as between separate benthic communities. We argue that these patterns reflect the heterogeneity and diversity of benthic habitats. A comparison of all OTUs against the Protist Ribosomal Reference database showed that a higher proportion of benthic than planktonic protist diversity is missing from public databases; similar results were obtained by comparing all OTUs against environmental references from NCBI's Short Read Archive. We suggest that the benthic realm may therefore be the world's largest reservoir of marine protist diversity, with most taxa at present undescribed.

Declines in both redundant and trace species characterize the latitudinal diversity gradient in tintinnid ciliates

The latitudinal diversity gradient is a well-known biogeographic pattern. However, rarely considered is how a cline in species richness may be reflected in the characteristics of species assemblages. Fewer species may equal fewer distinct ecological types, or declines in redundancy (species functionally similar to one another) or fewer trace species, those occurring in very low concentrations. We focused on tintinnid ciliates of the microzooplankton in which the ciliate cell is housed inside a species-specific lorica or shell. The size of lorica oral aperture, the lorica oral diameter (LOD), is correlated with a preferred prey size and maximum growth rate. Consequently, species of a distinct LOD are distinct in key ecologic characteristics, whereas those of a similar LOD are functionally similar or redundant species. We sampled from East Sea/Sea of Japan to the High Arctic Sea. We determined abundance distributions of biological species and also ecological types by grouping species in LOD size-classes, sets of ecologically similar species. In lower latitudes there are more trace species, more size-classes and the dominant species are accompanied by many apparently ecologically similar species, presumably able to replace the dominant species, at least with regard to the size of prey exploited. Such redundancy appears to decline markedly with latitude in assemblages of tintinnid ciliates. Furthermore, the relatively small species pools of the northern high latitude assemblages suggest a low capacity to adapt to changing conditions.

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.

Testing ecological theories with sequence similarity networks: marine ciliates exhibit similar geographic dispersal patterns as multicellular organisms

BACKGROUND

High-throughput sequencing technologies are lifting major limitations to molecular-based ecological studies of eukaryotic microbial diversity, but analyses of the resulting millions of short sequences remain a major bottleneck for these approaches. Here, we introduce the analytical and statistical framework of sequence similarity networks, increasingly used in evolutionary studies and graph theory, into the field of ecology to analyze novel pyrosequenced V4 small subunit rDNA (SSU-rDNA) sequence data sets in the context of previous studies, including SSU-rDNA Sanger sequence data from cultured ciliates and from previous environmental diversity inventories.

RESULTS

Our broadly applicable protocol quantified the progress in the description of genetic diversity of ciliates by environmental SSU-rDNA surveys, detected a fundamental historical bias in the tendency to recover already known groups in these surveys, and revealed substantial amounts of hidden microbial diversity. Moreover, network measures demonstrated that ciliates are not globally dispersed, but are structured by habitat and geographical location at intermediate geographical scale, as observed for bacteria, plants, and animals.

CONCLUSIONS

Currently available 'universal' primers used for local in-depth sequencing surveys provide little hope to exhaust the significantly higher ciliate (and most likely microbial) diversity than previously thought. Network analyses such as presented in this study offer a promising way to guide the design of novel primers and to further explore this vast and structured microbial diversity.

Protistan diversity in a permanently stratified meromictic lake (Lake Alatsee, SW Germany)

Protists play a crucial role for ecosystem function(ing) and oxygen is one of the strongest barriers against their local dispersal. However, protistan diversity in freshwater habitats with oxygen gradients received very little attention. We applied high-throughput sequencing of the V9 region (18S rRNA gene) to provide a hitherto unique spatiotemporal analysis of protistan diversity along the oxygen gradient of a freshwater meromictic lake (Lake Alatsee, SW Germany). In the mixolimnion, the communities experienced most seasonal structural changes, with Stramenopiles dominating in autumn and Dinoflagellata in summer. The suboxic interface supported the highest diversity, but only 23 OTUs95% (mainly Euglenozoa, after quality check and removal of operational taxonomic units (OTUs) with less than three sequences) were exclusively associated with this habitat. Eukaryotic communities in the anoxic monimolimnion showed the most stable seasonal pattern, with Chrysophyta and Bicosoecida being the dominant taxa. Our data pinpoint to the ecological role of the interface as a short-term 'meeting point' for protists, contributing to the coupling of the mixolimnion and the monimolimnion. Our analyses of divergent genetic diversity suggest a high degree of previously undescribed OTUs. Future research will have to reveal if this result actually points to a high number of undescribed species in such freshwater habitats.

Comparing High-throughput Platforms for Sequencing the V4 Region of SSU-rDNA in Environmental Microbial Eukaryotic Diversity Surveys

High-throughput sequencing platforms are continuing to increase resulting read lengths, which is allowing for a deeper and more accurate depiction of environmental microbial diversity. With the nascent Reagent Kit v3, Illumina MiSeq now has the ability to sequence the eukaryotic hyper-variable V4 region of the SSU-rDNA locus with paired-end reads. Using DNA collected from soils with analyses of strictly- and nearly identical amplicons, here we ask how the new Illumina MiSeq data compares with what we can obtain with Roche/454 GS FLX with regard to quantity and quality, presence and absence, and abundance perspectives. We show that there is an easy qualitative transition from the Roche/454 to the Illumina MiSeq platforms. The ease of this transition is more nuanced quantitatively for low-abundant amplicons, although estimates of abundances are known to also vary within platforms.