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

Ciliate diversity in rodrigo de freitas lagoon (Rio de Janeiro, Brazil) from an integrative standpoint

The Rodrigo de Freitas Lagoon is a highly eutrophic lacustrine system and has one of the longest histories of exploration and anthropic alteration in Brazil. Despite its relevance, limited studies explored the diversity of micro-eukaryotes in the lagoon. Ciliates (Alveolata, Ciliophora) are overlooked in environmental microbiology, especially in tropical and subtropical ecosystems, resulting in limited knowledge about their diversity and functional relevance in South American habitats, particularly in coastal lagoons. To fill this gap, here we investigated the diversity of ciliates in a brackish coastal lagoon in an urban area of Rio de Janeiro, Brazil, applying and comparing the performance of morphological and metabarcoding approaches. The metabarcoding analysis, based on high-throughput sequencing of the hipervariable region V4 of the 18S rRNA genes detected 37 molecular operational taxonomic units (MOTUs) assigned to Ciliophora, representing only about a half (56.9%) of the diversity detected by microscopy, which counted 65 ciliate morphotypes. The most representative classes in both approaches were Spirotrichea and Oligohymenophorea. The metabarcoding analysis revealed that 35.3% of the ciliate MOTUs had less than 97% similarity to available sequences in the NCBI database, indicating that more than one-third of these MOTUs potentially represents still not represented or undescribed ciliate species in current databases. Our findings indicate that metabarcoding techniques can significantly enhance the comprehension of ciliate diversity in tropical environments, but the scarcity of reference sequences of brackish ciliates in molecular databases represents a challenge to the taxonomic assignment of the MOTUs. This study provides new insights into the diversity of ciliates in a threatened coastal lagoon, revealing a vast array of still unknown and rare ciliate taxonomic units in tropical environments.

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.

Information Scale Correction for Varying Length Amplicons Improves Eukaryotic Microbiome Data Integration

The integration and reanalysis of big data provide valuable insights into microbiome studies. However, the significant difference in information scale between amplicon data poses a key challenge in data analysis. Therefore, reducing batch effects is crucial to enhance data integration for large-scale molecular ecology data. To achieve this, the information scale correction (ISC) step, involving cutting different length amplicons into the same sub-region, is essential. In this study, we used the Hidden Markov model (HMM) method to extract 11 different 18S rRNA gene v4 region amplicon datasets with 578 samples in total. The length of the amplicons ranged from 344 bp to 720 bp, depending on the primer position. By comparing the information scale correction of amplicons with varying lengths, we explored the extent to which the comparability between samples decreases with increasing amplicon length. Our method was shown to be more sensitive than V-Xtractor, the most popular tool for performing ISC. We found that near-scale amplicons exhibited no significant change after ISC, while larger-scale amplicons exhibited significant changes. After the ISC treatment, the similarity among the data sets improved, especially for long amplicons. Therefore, we recommend adding ISC processing when integrating big data, which is crucial for unlocking the full potential of microbial community studies and advancing our knowledge of 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.

Ciliate Diversity From Aquatic Environments in the Brazilian Atlantic Forest as Revealed by High-Throughput DNA Sequencing

Rainforest aquatic ecosystems include complex habitats with scarce information on their unicellular eukaryote diversity and community structure. We have investigated the diversity of ciliates in freshwater and brackish environments along the Brazilian Atlantic Forest, based on the hypervariable V4 region of the 18S-rDNA obtained by high-throughput DNA sequencing. Our analyses detected 409 ciliate taxonomic units (OTUs), mostly attributed to the classes Oligohymenophorea and Spirotrichea. A total of 11 classes, 12 subclasses, 112 genera, and 144 species were reported. We found the following: (a) the ciliate communities are more diverse in freshwater- than in Atlantic Forest-associated brackish environments; (b) the ciliate communities are composed by a small amount of highly abundant OTUs, but a high number of low-abundant or rare OTUs; (c) nearly one-third of the ciliate OTUs share less than 97% sequence identity to reference sequences and (d) phylogenetic inference supports the hypothesis that the V4 region of the Ciliophora 18S-rDNA is a suitable marker for accurate evolutionary inferences at class level. Our results showed that a considerable fraction of the HTS-detected diversity of ciliates from Brazilian Atlantic Forest is not represented in the currently available molecular databases.

Bathymetric gradient shapes the community composition rather than the species richness of deep-sea benthic ciliates

The bathymetric gradient is one of the most important factors that regulate the distribution of life. However, community variations of benthic ciliates along bathymetric gradients in the deep sea remain rather unexplored. In this study, we hypothesize that in the deep sea, the bathymetric gradient shapes the benthic ciliate community composition rather than the species richness. Here, we evaluated the distribution patterns and drivers of benthic ciliate communities of an abyssal plain, a seamount, and a trench with water depths ranging from 800 m down to 6600 m by high throughput eDNA sequencing and statistical analyses. We observed no significant correlation between ciliate operated taxonomic unit (OTU) richness and water depth. A meta-analysis, which combined our previously published data from the neritic habitats, supports the notion that water depth exceeding 800 m has little effect on the richness of benthic ciliate species. In contrast, the composition of deep-sea ciliate communities was significantly distinct in different habitats along the bathymetric gradients. A SourceTracker analysis revealed extremely low connectivity among ciliate communities along the bathymetric gradients. More than 95% of the community dissimilarity in the deep-sea floor was attributed to species replacement, which might be caused by environmental sorting or historical constraints. Furthermore, the observed community variations could be ascribed more to water depth than to geographic distance. The findings imply that the strong force of environmental sorting along the bathymetric gradients and the low connectivity among the ciliate communities might lead to an isolated evolution. This could shape the community composition rather than the species richness, which is mainly determined by the limited nutrient availability and the extreme environmental conditions in the deep sea.

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.

EukRef-Ciliophora: a manually curated, phylogeny-based database of small subunit rRNA gene sequences of ciliates

High-throughput sequencing (HTS) surveys, among the most common approaches currently used in environmental microbiology, require reliable reference databases to be correctly interpreted. The EukRef Initiative (eukref.org) is a community effort to manually screen available small subunit (SSU) rRNA gene sequences and produce a public, high-quality and informative framework of phylogeny-based taxonomic annotations. In the context of EukRef, we present a database for the monophyletic phylum Ciliophora, one of the most complex, diverse and ubiquitous protist groups. We retrieved more than 11 500 sequences of ciliates present in GenBank (28% from identified isolates and 72% from environmental surveys). Our approach included the inference of phylogenetic trees for every ciliate lineage and produced the largest SSU rRNA tree of the phylum Ciliophora to date. We flagged approximately 750 chimeric or low-quality sequences, improved the classification of 70% of GenBank entries and enriched environmental and literature metadata by 30%. The performance of EukRef-Ciliophora is superior to the current SILVA database in classifying HTS reads from a global marine survey. Comprehensive outputs are publicly available to make the new tool a useful guide for non-specialists and a quick reference for experts.

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.