Are We Overestimating Protistan Diversity in Nature?

Harmonized coexistence of intragenomic variations in diatom Skeletonema strains

Metabarcoding analysis has been demonstrated to be an effective technology for monitoring diversity and dynamics of phytoplankton including Skeletonema species. Although molecular diversity uncovered in metabarcoding projects has generally been interpreted as sum of interspecies diversity and intraspecies diversity, accumulating evidence suggests that it also harbors unprecedentedly high levels of intra-genomic variations (IGVs). As up to thousands of amplicon sequence variants (ASVs) identified in a typical metabarcoding project can be annotated to be Skeletonema species, we hypothesize that substantial portions of these ASVs are contributed by IGVs. Here, the nature of IGVs in Skeletonema species was quantitatively analyzed by carrying out single-strain metabarcoding analysis of 18S rDNA V4 in 49 strains belonging to seven Skeletonema species. Results showed that each Skeletonema strain harbored a high level of IGVs as expected. While many Skeletonema strains each contained one dominant ASV and a substantial number of ASVs displaying much lower relative abundance, other Skeletonema strains each contained multiple ASVs with comparable or nearly equally abundances. Thus the co-existence of multiple dominant ASVs in a single cell indicated a tug-of-war of these variants in evolution, which may eventually result in harmonized coexistence of multiple dominant ASVs. A total of nine dominant ASVs and 652 non-dominant ASVs were found in 49 strains of seven Skeletonema species, indicating rich interspecies and intraspecies variations, and complex evolution of IGVs in genus of Skeletonema. The results confirmed that the extensive degree of IGVs was the main contributor to the high molecular diversity revealed by metabarcoding analysis. This study highlights the importance of quantitative characterization of IGVs in Skeletonema species for accurate interpretation of species diversity in metabarcoding analysis.

Intragenomic variations of the harmful algal bloom species Phaeocystis globosa through single-strain metabarcoding analysis

Recent studies demonstrated that the existence of high-levels intragenomic variations (IGVs) may lead to overinterpretation of species diversity and genetic diversity in metabarcoding analysis. In this study, IGVs of the V4 region of 18S rRNA gene (18S rDNA V4) in Phaeocystis globosa were ascertained through metabarcoding analysis results of 54 P. globosa strains. Each P. globosa strain harbored over 200 ASVs, in which one ASV was dominant with higher relative abundance than others. According to different dominant ASVs, 54 strains could be divided into two groups, suggesting the high genetic diversity of P. globosa. ASVs were shared in different strains, suggesting that IGVs were genuine existed, rather than sequencing errors. Metabarcoding analysis of field samples identified large numbers of IGVs. High levels of IGVs may be due to incomplete homogenization or existence of nonfunctional pseudogenes. An accurate understanding of IGVs is necessary for the correct interpretation of metabarcoding analysis results.

The Pathology of the Brain Eating Amoeba Naegleria fowleri

The genus Naegleria is a taxonomic subfamily consisting of 47 free-living amoebae. The genus can be found in warm aqueous or soil habitats worldwide. The species Naegleria fowleri is probably the best-known species of this genus. As a facultative parasite, the protist is not dependent on hosts to complete its life cycle. However, it can infect humans by entering the nose during water contact, such as swimming, and travel along the olfactory nerve to the brain. There it causes a purulent meningitis (primary amoebic meningoencephalitis or PAME). Symptoms are severe and death usually occurs within the first week. PAME is a frightening infectious disease for which there is neither a proven cure nor a vaccine. In order to contain the disease and give patients any chance to survival, action must be taken quickly. A rapid diagnosis is therefore crucial. PAME is diagnosed by the detection of amoebae in the liquor and later in the cerebrospinal fluid. For this purpose, CSF samples are cultured and stained and finally examined microscopically. Molecular techniques such as PCR or ELISA support the microscopic analysis and secure the diagnosis.

Comparative assessment of the intragenomic variations of dinoflagellate Tripos species through single-cell sequencing

Tripos is a large dinoflagellate genus widely distributed in the world's oceans. Morphology-based species identification is inconclusive due to high morphological intraspecific variability. Metabarcoding analysis has been demonstrated to be effective for species identification and tracking their spatiotemporal dynamics. However, accumulating evidence suggests high levels of intragenomic variations (IGVs) are common in many algae, leading to concerns about overinterpretation of molecular diversity in metabarcoding studies. In this project, we evaluated and compared IGVs in Tripos species by conducting the first high-throughput sequencing (HTS) of 18S rDNA V4 of Tripos single cells. High numbers of haplotypes (19-172) were identified in each of the 30 Tripos cells. Each cell contained one dominant haplotype with high relative abundance and many haplotypes with lower abundances. Thus, the presence of multiple minor haplotypes substantially overestimate the molecular diversity identified in metabarcoding analysis, which encompass not only interspecific and intraspecific diversities, but high levels of IGVs.

Species delimitation 4.0: integrative taxonomy meets artificial intelligence

Although species are central units for biological research, recent findings in genomics are raising awareness that what we call species can be ill-founded entities due to solely morphology-based, regional species descriptions. This particularly applies to groups characterized by intricate evolutionary processes such as hybridization, polyploidy, or asexuality. Here, challenges of current integrative taxonomy (genetics/genomics + morphology + ecology, etc.) become apparent: different favored species concepts, lack of universal characters/markers, missing appropriate analytical tools for intricate evolutionary processes, and highly subjective ranking and fusion of datasets. Now, integrative taxonomy combined with artificial intelligence under a unified species concept can enable automated feature learning and data integration, and thus reduce subjectivity in species delimitation. This approach will likely accelerate revising and unraveling eukaryotic biodiversity.

The global genetic diversity of planktonic foraminifera reveals the structure of cryptic speciation in plankton

The nature and extent of diversity in the plankton has fascinated scientists for over a century. Initially, the discovery of many new species in the remarkably uniform and unstructured pelagic environment appeared to challenge the concept of ecological niches. Later, it became obvious that only a fraction of plankton diversity had been formally described, because plankton assemblages are dominated by understudied eukaryotic lineages with small size that lack clearly distinguishable morphological features. The high diversity of the plankton has been confirmed by comprehensive metabarcoding surveys, but interpretation of the underlying molecular taxonomies is hindered by insufficient integration of genetic diversity with morphological taxonomy and ecological observations. Here we use planktonic foraminifera as a study model and reveal the full extent of their genetic diversity and investigate geographical and ecological patterns in their distribution. To this end, we assembled a global data set of ~7600 ribosomal DNA sequences obtained from morphologically characterised individual foraminifera, established a robust molecular taxonomic framework for the observed diversity, and used it to query a global metabarcoding data set covering ~1700 samples with ~2.48 billion reads. This allowed us to extract and assign 1 million reads, enabling characterisation of the structure of the genetic diversity of the group across ~1100 oceanic stations worldwide. Our sampling revealed the existence of, at most, 94 distinct molecular operational taxonomic units (MOTUs) at a level of divergence indicative of biological species. The genetic diversity only doubles the number of formally described species identified by morphological features. Furthermore, we observed that the allocation of genetic diversity to morphospecies is uneven. Only 16 morphospecies disguise evolutionarily significant genetic diversity, and the proportion of morphospecies that show genetic diversity increases poleward. Finally, we observe that MOTUs have a narrower geographic distribution than morphospecies and that in some cases the MOTUs belonging to the same morphospecies (cryptic species) have different environmental preferences. Overall, our analysis reveals that even in the light of global genetic sampling, planktonic foraminifera diversity is modest and finite. However, the extent and structure of the cryptic diversity reveals that genetic diversification is decoupled from morphological diversification, hinting at different mechanisms acting at different levels of divergence.

Measurements of soil protist richness and community composition are influenced by primer pair, annealing temperature, and bioinformatics choices

Protists are a diverse and understudied group of microbial eukaryotic organisms especially in terrestrial environments. Advances in molecular methods are increasing our understanding of the distribution and functions of these creatures; however, there is a vast array of choices researchers make including barcoding genes, primer pairs, PCR settings, and bioinformatic options that can impact the outcome of protist community surveys. Here, we tested four commonly used primer pairs targeting the V4 and V9 regions of the 18S rRNA gene using different PCR annealing temperatures and processed the sequences with different bioinformatic parameters in 10 diverse soils to evaluate how primer pair, amplification parameters, and bioinformatic choices influence the composition and richness of protist and non-protist taxa using Illumina sequencing. Our results showed that annealing temperature influenced sequencing depth and protist taxon richness for most primer pairs, and that merging forward and reverse sequencing reads for the V4 primer pairs dramatically reduced the number of sequences and taxon richness of protists. The data sets of primers that targeted the same 18S rRNA gene region (e.g., V4 or V9) had similar protist community compositions; however, data sets from primers targeting the V4 18S rRNA gene region detected a greater number of protist taxa compared to those prepared with primers targeting the V9 18S rRNA region. There was limited overlap of protist taxa between data sets targeting the two different gene regions (80/549 taxa). Together, we show that laboratory and bioinformatic choices can substantially affect the results and conclusions about protist diversity and community composition using metabarcoding.IMPORTANCEEcosystem functioning is driven by the activity and interactions of the microbial community, in both aquatic and terrestrial environments. Protists are a group of highly diverse, mostly unicellular microbes whose identity and roles in terrestrial ecosystem ecology have been largely ignored until recently. This study highlights the importance of choices researchers make, such as primer pair, on the results and conclusions about protist diversity and community composition in soils. In order to better understand the roles protist taxa play in terrestrial ecosystems, biases in methodological and analytical choices should be understood and acknowledged.

Time-series metagenomics reveals changing protistan ecology of a temperate dimictic lake

BACKGROUND

Protists, single-celled eukaryotic organisms, are critical to food web ecology, contributing to primary productivity and connecting small bacteria and archaea to higher trophic levels. Lake Mendota is a large, eutrophic natural lake that is a Long-Term Ecological Research site and among the world's best-studied freshwater systems. Metagenomic samples have been collected and shotgun sequenced from Lake Mendota for the last 20 years. Here, we analyze this comprehensive time series to infer changes to the structure and function of the protistan community and to hypothesize about their interactions with bacteria.

RESULTS

Based on small subunit rRNA genes extracted from the metagenomes and metagenome-assembled genomes of microeukaryotes, we identify shifts in the eukaryotic phytoplankton community over time, which we predict to be a consequence of reduced zooplankton grazing pressures after the invasion of a invasive predator (the spiny water flea) to the lake. The metagenomic data also reveal the presence of the spiny water flea and the zebra mussel, a second invasive species to Lake Mendota, prior to their visual identification during routine monitoring. Furthermore, we use species co-occurrence and co-abundance analysis to connect the protistan community with bacterial taxa. Correlation analysis suggests that protists and bacteria may interact or respond similarly to environmental conditions. Cryptophytes declined in the second decade of the timeseries, while many alveolate groups (e.g., ciliates and dinoflagellates) and diatoms increased in abundance, changes that have implications for food web efficiency in Lake Mendota.

CONCLUSIONS

We demonstrate that metagenomic sequence-based community analysis can complement existing efforts to monitor protists in Lake Mendota based on microscopy-based count surveys. We observed patterns of seasonal abundance in microeukaryotes in Lake Mendota that corroborated expectations from other systems, including high abundance of cryptophytes in winter and diatoms in fall and spring, but with much higher resolution than previous surveys. Our study identified long-term changes in the abundance of eukaryotic microbes and provided context for the known establishment of an invasive species that catalyzes a trophic cascade involving protists. Our findings are important for decoding potential long-term consequences of human interventions, including invasive species introduction. Video Abstract.

Integrating depth-dependent protist dynamics and microbial interactions in spring succession of a freshwater reservoir

BACKGROUND

Protists are essential contributors to eukaryotic diversity and exert profound influence on carbon fluxes and energy transfer in freshwaters. Despite their significance, there is a notable gap in research on protistan dynamics, particularly in the deeper strata of temperate lakes. This study aimed to address this gap by integrating protists into the well-described spring dynamics of Římov reservoir, Czech Republic. Over a 2-month period covering transition from mixing to established stratification, we collected water samples from three reservoir depths (0.5, 10 and 30 m) with a frequency of up to three times per week. Microbial eukaryotic and prokaryotic communities were analysed using SSU rRNA gene amplicon sequencing and dominant protistan groups were enumerated by Catalysed Reporter Deposition-Fluorescence in situ Hybridization (CARD-FISH). Additionally, we collected samples for water chemistry, phyto- and zooplankton composition analyses.

RESULTS

Following the rapid changes in environmental and biotic parameters during spring, protistan and bacterial communities displayed swift transitions from a homogeneous community to distinct strata-specific communities. A prevalence of auto- and mixotrophic protists dominated by cryptophytes was associated with spring algal bloom-specialized bacteria in the epilimnion. In contrast, the meta- and hypolimnion showcased a development of a protist community dominated by putative parasitic Perkinsozoa, detritus or particle-associated ciliates, cercozoans, telonemids and excavate protists (Kinetoplastida), co-occurring with bacteria associated with lake snow.

CONCLUSIONS

Our high-resolution sampling matching the typical doubling time of microbes along with the combined microscopic and molecular approach and inclusion of all main components of the microbial food web allowed us to unveil depth-specific populations' successions and interactions in a deep lentic ecosystem.

Genomic insights into the cellular specialization of predation in raptorial protists

BACKGROUND

Predation is a fundamental mechanism for organisms to acquire energy, and various species have evolved diverse tools to enhance their hunting abilities. Among protozoan predators, raptorial Haptorian ciliates are particularly fascinating as they possess offensive extrusomes known as toxicysts, which are rapidly discharged upon prey contact. However, our understanding of the genetic processes and specific toxins involved in toxicyst formation and discharge is still limited.

RESULTS

In this study, we investigated the predation strategies and subcellular structures of seven Haptoria ciliate species and obtained their genome sequences using single-cell sequencing technology. Comparative genomic analysis revealed distinct gene duplications related to membrane transport proteins and hydrolytic enzymes in Haptoria, which play a crucial role in the production and discharge of toxicysts. Transcriptomic analysis further confirmed the abundant expression of genes related to membrane transporters and cellular toxins in Haptoria compared to Trichostomatia. Notably, polyketide synthases (PKS) and L-amino acid oxidases (LAAO) were identified as potentially toxin genes that underwent extensive duplication events in Haptoria.

CONCLUSIONS

Our results shed light on the evolutionary and genomic adaptations of Haptorian ciliates for their predation strategies in evolution and provide insights into their toxic mechanisms.

Extensive intragenomic variations of the 18S rDNA V4 region in the toxigenic diatom species Pseudo-nitzschia multistriata revealed through high-throughput sequencing

Metabarcoding analysis is an effective technique for monitoring the domoic acid-producing Pseudo-nitzschia species in marine environments, uncovering high-levels of molecular diversity. However, such efforts may result in the overinterpretation of Pseudo-nitzschia species diversity, as molecular diversity not only encompasses interspecies and intraspecies diversities but also exhibits extensive intragenomic variations (IGVs). In this study, we analyzed the V4 region of the 18S rDNA of 30 strains of Pseudo-nitzschia multistriata collected from the coasts of China. The results showed that each P. multistriata strain harbored about a hundred of unique 18S rDNA V4 sequence varieties, of which each represented by a unique amplicon sequence variant (ASV). This study demonstrated the extensive degree of IGVs in P. multistriata strains, suggesting that IGVs may also present in other Pseudo-nitzschia species and other phytoplankton species. Understanding the scope and levels of IGVs is crucial for accurately interpreting the results of metabarcoding analysis.

The high molecular diversity in Noctiluca scintillans is dominated by intra-genomic variations revealed by single cell high-throughput sequencing of 18S rDNA V4

The application of high-throughput sequencing (HTS) technologies has revolutionized research on phytoplankton biodiversity by generating an unprecedented amount of molecular data in marine ecosystem surveys. However, high-level of molecular diversity uncovered in HTS-based metabarcoding analyses may lead to overinterpretation of phytoplankton diversity due to excessive intra-genomic variations (IGVs). The aims in this study are to explore the nature of phytoplankton molecular diversity and to test the hypothesis. We carried out single-cell metabarcoding analysis of 18S rDNA V4 sequences obtained in single Noctiluca scintillans cells isolated from various sites in coastal waters of China. Results showed that each single N. scintillans cell harbored a high level of IGVs with about 100 amplicon sequence variants (ASVs). The large numbers of non-dominant ASVs identified in N. scintillans cells, which might correspond to the larger numbers of ASVs annotated as N. scintillans and showed similar temporal dynamics in metabarcoding analyses, could inflate the inter-species diversity or intra-species genetic diversity. In addition, there were large numbers of additional ASVs that were not annotated as N. scintillans. These non-N. scintillans ASVs might represent diverse preys for N. scintillans, consistent with previous reports that N. scintillans may act as chance predator of a broad-spectrum preys. This single-cell study has unambiguously demonstrated that the existence of high levels of IGVs in N. scintillans and most likely many other phytoplankton species, demonstrating that the majority of the molecular diversity revealed in metabarcoding analysis, which were generally interpreted as the sum of inter-species diversity and intra-species diversity, actually included high levels of IGVs and should be interpreted with caution.

Predation has a significant impact on the complexity and stability of microbial food webs in subalpine lakes

IMPORTANCE

As an important part of microbial food webs, protists transfer organic carbon and nutrients to higher trophic levels in aquatic ecosystems. Protist predation often influences the abundance and composition of bacterial communities. However, we still do not understand whether and how predation affects the complexity and stability of microbial food webs. This study assessed the seasonal dynamic characteristics and driving factors of microbial food webs in terms of complexity and stability. Our findings have implications for future surveys to reveal the effects of climate and environmental changes.

Parasites of firebugs in Austria with focus on the "micro"-diversity of the cosmopolitan trypanosomatid Leptomonas pyrrhocoris

In this work, we investigated parasites of the firebug Pyrrhocoris apterus in Austria and demonstrated that in addition to the extensively studied Leptomonas pyrrhocoris, it can also be infected by Blastocrithidia sp. and by a mermithid, which for the first time has been characterized using molecular methods. This diversity can be explained by the gregarious lifestyle, as well as the coprophagous and cannibalistic behavior of the insect hosts that makes them susceptible to various parasites. In addition, we showed no tight association of the L. pyrrhocoris haplotypes and geographical locations (at least, considering the relatively small scale of locations in Austria) implying that the natural populations of L. pyrrhocoris are mixed due to the mobility of their firebug hosts.

Contemporary issues, current best practice and ways forward in soil protist ecology

Soil protists are increasingly studied due to a release from previous methodological constraints and the acknowledgement of their immense diversity and functional importance in ecosystems. However, these studies often lack sufficient depth in knowledge, which is visible in the form of falsely used terms and false- or over-interpreted data with conclusions that cannot be drawn from the data obtained. As we welcome that also non-experts include protists in their still mostly bacterial and/or fungal-focused studies, our aim here is to help avoid some common errors. We provide suggestions for current terms to use when working on soil protists, like protist instead of protozoa, predator instead of grazer, microorganisms rather than microflora and other terms to be used to describe the prey spectrum of protists. We then highlight some dos and don'ts in soil protist ecology including challenges related to interpreting 18S rRNA gene amplicon sequencing data. We caution against the use of standard bioinformatic settings optimized for bacteria and the uncritical reliance on incomplete and partly erroneous reference databases. We also show why causal inferences cannot be drawn from sequence-based correlation analyses or any sampling/monitoring, study in the field without thorough experimental confirmation and sound understanding of the biology of taxa. Together, we envision this work to help non-experts to more easily include protists in their soil ecology analyses and obtain more reliable interpretations from their protist data and other biodiversity data that, in the end, will contribute to a better understanding of soil ecology.

Intragenomic diversity of the V9 hypervariable domain in eukaryotes has little effect on metabarcoding

Metabarcoding revolutionized our understanding of diversity and ecology of microorganisms in different habitats. However, it is also associated with several inherent biases, one of which is associated with intragenomic diversity of a molecular barcode. Here, we compare intragenomic variability of the V9 region of the 18S rRNA gene in 19 eukaryotic phyla abundant in marine plankton. The level of intragenomic variability is comparable across all the phyla, and in most genomes and transcriptomes one V9 sequence and one OTU is predominant. However, most of the variability observed at the barcode level is probably caused by sequencing errors and is mitigated by using a denoising tool, DADA2. The SWARM algorithm commonly used in metabarcoding studies is not optimal for collapsing genuine and erroneous sequences into a single OTU, leading to an overestimation of diversity in metabarcoding data. For an unknown reason, SWARM inflates diversity of eupelagonemids more than that of other eukaryotes.

Assembly processes and functional diversity of marine protists and their rare biosphere

BACKGROUND

The mechanisms shaping the rare microbial biosphere and its role in ecosystems remain unclear. We developed an approach to study ecological patterns in the rare biosphere and use it on a vast collection of marine microbiomes, sampled in coastal ecosystems at a regional scale. We study the assembly processes, and the ecological strategies constituting the rare protistan biosphere. Using the phylogeny and morpho-trophic traits of these protists, we also explore their functional potential.

RESULTS

Taxonomic community composition remained stable along rank abundance curves. Conditionally rare taxa, driven by selection processes, and transiently rare taxa, with stochastic distributions, were evidenced along the rank abundance curves of all size-fractions. Specific taxa within the divisions Sagenista, Picozoa, Telonemia, and Choanoflagellida were rare across time and space. The distribution of traits along rank abundance curves outlined a high functional redundancy between rare and abundant protists. Nevertheless, trophic traits illustrated an interplay between the trophic groups of different size-fractions.

CONCLUSIONS

Our results suggest that rare and abundant protists are evolutionary closely related, most notably due to the high microdiversity found in the rare biosphere. We evidenced a succession of assembly processes and strategies of rarity along rank abundance curves that we hypothesize to be common to most microbiomes at the regional scale. Despite high functional redundancy in the rare protistan biosphere, permanently rare protists were evidenced, and they could play critical functions as bacterivores and decomposers from within the rare biosphere. Finally, changes in the composition of the rare protistan biosphere could be influenced by the trophic regime of aquatic ecosystems. Our work contributes to understanding the role of rare protists in microbiomes.

Daily dynamics of contrasting spring algal blooms in Santa Monica Bay (central Southern California Bight)

Protistan algae (phytoplankton) dominate coastal upwelling ecosystems where they form massive blooms that support the world's most important fisheries and constitute an important sink for atmospheric CO₂ . Bloom initiation is well understood, but the biotic and abiotic forces that shape short-term dynamics in community composition are still poorly characterized. Here, high-frequency (daily) changes in relative abundance dynamics of the metabolically active protistan community were followed via expressed 18S V4 rRNA genes (RNA) throughout two algal blooms during the spring of 2018 and 2019 in Santa Monica Bay (central Southern California Bight). A diatom bloom formed after wind-driven, nutrient upwelling events in both years, but different taxa dominated each year. Whereas diatoms bloomed following elevated nutrients and declined after depletion each year, a massive dinoflagellate bloom manifested under relatively low inorganic nitrogen conditions following diatom bloom senescence in 2019 but not 2018. Network analysis revealed associations between diatoms and cercozoan putative parasitic taxa and syndinean parasites during 2019 that may have influenced the demise of the diatoms, and the transition to a dinoflagellate-dominated bloom.

Does filter pore size introduce bias in DNA sequence-based plankton community studies?

The cell size of microbial eukaryotic plankton normally ranges from 0.2 to 200 μm. During the past decade, high-throughput sequencing of DNA has been revolutionizing their study on an unprecedented scale. Nonetheless, it is currently unclear whether we can accurately, effectively, and quantitatively depict the microbial eukaryotic plankton community using size-fractionated filtration combined with environmental DNA (eDNA) molecular methods. Here we assessed the microbial eukaryotic plankton communities with two filtering strategies from two subtropical reservoirs, that is one-step filtration (0.2–200 μm) and size-fractionated filtration (0.2–3 and 3–200 μm). The difference of 18S rRNA gene copy abundance between the two filtering treatments was less than 50% of the 0.2–200 μm microbial eukaryotic community for 95% of the total samples. Although the microbial eukaryotic plankton communities within the 0.2–200 μm and the 0.2–3 and 3–200 μm size fractions had approximately identical 18S rRNA gene copies, there were significant differences in their community composition. Furthermore, our results demonstrate that the systemic bias introduced by size-fractionation filtration has more influence on unique OTUs than shared OTUs, and the significant differences in abundance between the two eukaryotic plankton communities largely occurred in low-abundance OTUs in specific seasons. This work provides new insights into the use of size-fractionation in molecular studies of microbial eukaryotes populating the plankton.

Estuarine microbial networks and relationships vary between environmentally distinct communities

Microbial interactions have profound impacts on biodiversity, biogeochemistry, and ecosystem functioning, and yet, they remain poorly understood in the ocean and with respect to changing environmental conditions. We applied hierarchical clustering of an annual 16S and 18S amplicon dataset in the Skidaway River Estuary, which revealed two similar clusters for prokaryotes (Bacteria and Archaea) and protists: Cluster 1 (March-May and November-February) and Cluster 2 (June-October). We constructed co-occurrence networks from each cluster to explore how microbial networks and relationships vary between environmentally distinct periods in the estuary. Cluster 1 communities were exposed to significantly lower temperature, sunlight, NO3, and SiO4; only NH4 was higher at this time. Several network properties (e.g., edge number, degree, and centrality) were elevated for networks constructed with Cluster 1 vs. 2 samples. There was also evidence that microbial nodes in Cluster 1 were more connected (e.g., higher edge density and lower path length) compared to Cluster 2, though opposite trends were observed when networks considered Prokaryote-Protist edges only. The number of Prokaryote-Prokaryote and Prokaryote-Protist edges increased by >100% in the Cluster 1 network, mainly involving Flavobacteriales, Rhodobacterales, Peridiniales, and Cryptomonadales associated with each other and other microbial groups (e.g., SAR11, Bacillariophyta, and Strombidiida). Several Protist-Protist associations, including Bacillariophyta correlated with Syndiniales (Dino-Groups I and II) and an Unassigned Dinophyceae group, were more prevalent in Cluster 2. Based on the type and sign of associations that increased in Cluster 1, our findings indicate that mutualistic, competitive, or predatory relationships may have been more representative among microbes when conditions were less favorable in the estuary; however, such relationships require further exploration and validation in the field and lab. Coastal networks may also be driven by shifts in the abundance of certain taxonomic or functional groups. Sustained monitoring of microbial communities over environmental gradients, both spatial and temporal, is critical to predict microbial dynamics and biogeochemistry in future marine ecosystems.

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.

Hierarchical phylogenetic community assembly of soil protists in a temperate agricultural field

Protists are abundant, diverse and perform essential functions in soils. Protistan community structure and its change across time or space are traditionally studied at the species-level but the relative importance of the processes shaping these patterns depends on the taxon phylogenetic resolution. Using 18S rDNA amplicon data of the Cercozoa, a group of dominant soil protists, from an agricultural field in western Germany, we observed a turnover of relatively closely related taxa (from sequence variants to genus-level clades) across soil depth; while across soil habitats (rhizosphere, bulk soil, drilosphere) we observed turnover of relatively distantly related taxa, confirming Paracercomonadidae as a rhizosphere-associated clade. We extended our approach to show that closely related Cercozoa encounter divergent AM fungi across soil depth and that distantly related Cercozoa encounter closely related AM fungi across soil compartments. This study suggests that soil Cercozoa community assembly at the field-scale is driven by niche-based processes shaped by evolutionary legacy of adaptation to conditions primarily related to soil compartment, followed by soil layer, giving a deeper understanding on the selection pressures that shaped their evolution.

Genome architecture used to supplement species delineation in two cryptic marine ciliates

The purpose of this study is to determine which taxonomic methods can elucidate clear and quantifiable differences between two cryptic ciliate species, and to test the utility of genome architecture as a new diagnostic character in the discrimination of otherwise indistinguishable taxa. Two cryptic tintinnid ciliates, Schmidingerella arcuata and Schmidingerella meunieri, are compared via traditional taxonomic characters including lorica morphometrics, ribosomal RNA (rRNA) gene barcodes and ecophysiological traits. In addition, single-cell 'omics analyses (single-cell transcriptomics and genomics) are used to elucidate and compare patterns of micronuclear genome architecture between the congeners. The results include a highly similar lorica that is larger in S. meunieri, a 0%-0.5% difference in rRNA gene barcodes, two different and nine indistinguishable growth responses among 11 prey treatments, and distinct patterns of micronuclear genomic architecture for genes detected in both ciliates. Together, these results indicate that while minor differences exist between S. arcuata and S. meunieri in common indices of taxonomic identification (i.e., lorica morphology, DNA barcode sequences and ecophysiology), differences exist in their genomic architecture, which suggests potential genetic incompatibility. Different patterns of micronuclear architecture in genes shared by both isolates also enable the design of species-specific primers, which are used in this study as unique "architectural barcodes" to demonstrate the co-occurrence of both ciliates in samples collected from a NW Atlantic estuary. These results support the utility of genomic architecture as a tool in species delineation, especially in ciliates that are cryptic or otherwise difficult to differentiate using traditional methods of identification.

Intra-genomic rRNA gene variability of Nassellaria and Spumellaria (Rhizaria, Radiolaria) assessed by Sanger, MinION and Illumina sequencing

Ribosomal RNA (rRNA) genes are known to be valuable markers for the barcoding of eukaryotic life and its phylogenetic classification at various taxonomic levels. The large-scale exploration of environmental microbial diversity through metabarcoding approaches have been focused mainly on the V4 and V9 regions of the 18S rRNA gene. The accurate interpretation of such environmental surveys is hampered by technical (e.g., PCR and sequencing errors) and biological biases (e.g., intra-genomic variability). Here we explored the intra-genomic diversity of Nassellaria and Spumellaria specimens (Radiolaria) by comparing Sanger sequencing with Illumina and Oxford Nanopore Technologies (MinION). Our analysis determined that intra-genomic variability of Nassellaria and Spumellaria is generally low, yet some Spumellaria specimens showed two different copies of the V4 with <97% similarity. From the different sequencing methods, Illumina showed the highest number of contaminations (i.e., environmental DNA, cross-contamination, tag-jumping), revealed by its high sequencing depth; and MinION showed the highest sequencing rate error (~14%). Yet the long reads produced by MinION (~2900 bp) allowed accurate phylogenetic reconstruction studies. These results highlight the requirement for a careful interpretation of Illumina based metabarcoding studies, in particular regarding low abundant amplicons, and open future perspectives towards full-length rDNA environmental metabarcoding surveys.

Diversity of Microbial Eukaryotes Along the West Antarctic Peninsula in Austral Spring

During a cruise from October to November 2019, along the West Antarctic Peninsula, between 64.32 and 68.37°S, we assessed the diversity and composition of the active microbial eukaryotic community within three size fractions: micro- (> 20 μm), nano- (20-5 μm), and pico-size fractions (5-0.2 μm). The communities and the environmental parameters displayed latitudinal gradients, and we observed a strong similarity in the microbial eukaryotic communities as well as the environmental parameters between the sub-surface and the deep chlorophyll maximum (DCM) depths. Chlorophyll concentrations were low, and the mixed layer was shallow for most of the 17 stations sampled. The richness of the microplankton was higher in Marguerite Bay (our southernmost stations), compared to more northern stations, while the diversity for the nano- and pico-plankton was relatively stable across latitude. The microplankton communities were dominated by autotrophs, mostly diatoms, while mixotrophs (phototrophs-consuming bacteria and kleptoplastidic ciliates, mostly alveolates, and cryptophytes) were the most abundant and active members of the nano- and picoplankton communities. While phototrophy was the dominant trophic mode, heterotrophy (mixotrophy, phagotrophy, and parasitism) tended to increase southward. The samples from Marguerite Bay showed a distinct community with a high diversity of nanoplankton predators, including spirotrich ciliates, and dinoflagellates, while cryptophytes were observed elsewhere. Some lineages were significantly related-either positively or negatively-to ice coverage (e.g., positive for Pelagophyceae, negative for Spirotrichea) and temperature (e.g., positive for Cryptophyceae, negative for Spirotrichea). This suggests that climate changes will have a strong impact on the microbial eukaryotic community.

Diversity and ecology of Radiolaria in modern oceans

Among the many inhabitants of planktonic communities, several lineages have biomineralized intricate skeletons. These have existed for millions of years and include the Radiolaria, a group of marine protists, many of which bear delicate mineral skeletons of different natures. Radiolaria are well known for their paleontological signatures, but little is known about the ecology of modern assemblages. They are found from polar to tropical regions, in the sunlit layers of the ocean down to the deep and cold bathypelagic. They are closely involved in the biogeochemical cycles of silica, carbon and strontium sulfate, carrying important amounts of such elements to the deep ocean. However, relatively little is known on the actual extent of genetic diversity or biogeographic patterns. The rapid emergence and acceptance of molecular approaches have nevertheless led to major advances in our understanding of diversity within and evolutionary relationships between major radiolarian groups. Here, we review the state of knowledge relating to the classification, diversity and ecology of extant radiolarian orders, highlighting the substantial gaps in our understanding of the extent of their contribution to marine biodiversity and their role in marine food webs.

Fine-scale genetic diversity and putative ecotypes of oxymonad protists coinhabiting the hindgut of Reticulitermes speratus

The hindgut of lower termites is generally coinhabited by multiple morphologically identifiable protist species. However, it is unclear how many protist species truly coexist in this miniaturized environment, and moreover, it is difficult to define the fundamental unit of protist diversity. Species delineation of termite gut protists has therefore been guided without a theory-based concept of species. Here, we focused on the hindgut of the termite Reticulitermes speratus, where 10 or 11 morphologically distinct oxymonad cell types, that is, morphospecies, coexist. We elucidated the phylogenetic structure of all co-occurring oxymonads and addressed whether their diversity can be explained by the "ecotype" hypothesis. Oxymonad-specific 18S rRNA gene amplicon sequencing analyses of whole-gut samples, combined with single-cell 18S rRNA sequencing of the oxymonad morphospecies, identified 210 one-nucleotide-level variants. The phylogenetic analysis of these variants revealed the presence of microdiverse clusters typically within 1% sequence divergence. Each known oxymonad morphospecies comprised one to several monophyletic or paraphyletic microdiverse clusters. Using these sequence data sets, we conducted computational simulation to predict the rates of ecotype formation and periodic selection, and to demarcate putative ecotypes. Our simulations suggested that the oxymonad genetic divergence is constrained primarily by strong selection, in spite of limited population size and possible bottlenecks during intergenerational transmission. A total of 33 oxymonad ecotypes were predicted, and most of the putative ecotypes were consistently detected among different colonies and host individuals. These findings provide a possible theoretical basis for species diversity and underlying mechanisms of coexistence of termite gut protists.

Diversity and ecology of protists revealed by metabarcoding

Protists are the dominant eukaryotes in the biosphere where they play key functional roles. While protists have been studied for over a century, it is the high-throughput sequencing of molecular markers from environmental samples - the approach of metabarcoding - that has revealed just how diverse, and abundant, these small organisms are. Metabarcoding is now routine to survey environmental diversity, so data have rapidly accumulated from a multitude of environments and at different sampling scales. This mass of data has provided unprecedented opportunities to study the taxonomic and functional diversity of protists, and how this diversity is organised in space and time. Here, we use metabarcoding as a common thread to discuss the state of knowledge in protist diversity research, from technical considerations of the approach to important insights gained on diversity patterns and the processes that might have structured this diversity. In addition to these insights, we conclude that metabarcoding is on the verge of an exciting added dimension thanks to the maturation of high-throughput long-read sequencing, so that a robust eco-evolutionary framework of protist diversity is within reach.

Phagotrophic protists (protozoa) in Antarctic terrestrial ecosystems: diversity, distribution, ecology, and best research practices

Phagotrophic protists (formerly protozoa) are a highly diverse, polyphyletic grouping of generally unicellular, heterotrophic eukaryotes that are key regulators of the soil microbiome. The biodiversity and ecology of soil phagotrophic protists are still largely uncharacterized, especially in the Antarctic, which possesses some of the harshest terrestrial environments known and potentially many physiologically unique and scientifically interesting species. Antarctic soil systems are also highly limited in terms of moisture, temperature, and carbon, and the resulting reduced biological complexity can facilitate fine-tuned investigation of the drivers and functioning of microbial communities. To facilitate and encourage future research into protist biodiversity and ecology, especially in context of the broader functioning of Antarctic terrestrial communities, I review the biodiversity, distribution, and ecology of Antarctic soil phagotrophic protists. Biodiversity appears to be highly structured by region and taxonomic group, with the Antarctic Peninsula having the highest taxonomic diversity and ciliates (Ciliophora) being the most diverse taxonomic group. However, richness estimates are likely skewed by disproportionate sampling (over half of the studies are from the peninsula), habitat type bias (predominately moss-associated soils), investigator bias (toward ciliates and the testate amoeba morphogroup), and methodological approach (toward cultivation and morphological identification). To remedy these biases, a standardized methodology using both morphological and molecular identification and increased emphasis on microflagellate and naked amoeba morphogroups is needed. Additionally, future research should transition away from biodiversity survey studies to dedicated ecological studies that emphasize the function, ecophysiology, endemicity, dispersal, and impact of abiotic drivers beyond moisture and temperature.

Seasonal and spatial variation of protist communities from reef water and open ocean water in patchy coral reef areas of a semi-enclosed bay

Protists are an important component of the marine ecosystem and play an essential role in material cycle and energy flow, but the distribution of protists in coral reefs have not been fully studied. In this study, high-throughput amplicon sequencing technology was used to study the biodiversity and community structure of protists from coral reefs and open sea areas, with the typical semi-enclosed bay Daya Bay as the research field. There were significant seasonal differences in the dominant phyla of protists, biodiversity index values and βeta diversity (P < 0.05) but no significant differences in the different sampling areas (P > 0.05). The topological parameters of the co-occurrence network showed the protist co-occurrence network in the open sea had more complex interactions and stronger stability than in the coral reef areas because of the hydrodynamics, waves, and relatively poor nutrients. Redundancy analysis and the Mantel test showed that the structure of the protist community was affected by seawater temperature, pH, salinity, and dissolved oxygen. This study analysed the temporal and spatial differences in protists in the coral reef and open sea areas of Daya Bay to provide important information for the study of protist biodiversity and community structure in semi-enclosed bays.

Exploring the Use of Environmental DNA (eDNA) to Detect Animal Taxa in the Mesopelagic Zone

Animal biodiversity in the ocean’s vast mesopelagic zone is relatively poorly studied due to technological and logistical challenges. Environmental DNA (eDNA) analyses show great promise for efficiently characterizing biodiversity and could provide new insight into the presence of mesopelagic species, including those that are missed by traditional net sampling. Here, we explore the utility of eDNA for identifying animal taxa. We describe the results from an August 2018 cruise in Slope Water off the northeast United States. Samples for eDNA analysis were collected using Niskin bottles during five CTD casts. Sampling depths along each cast were selected based on the presence of biomass as indicated by the shipboard Simrad EK60 echosounder. Metabarcoding of the 18S V9 gene region was used to assess taxonomic diversity. eDNA metabarcoding results were compared with those from net-collected (MOCNESS) plankton samples. We found that the MOCNESS sampling recovered more animal taxa, but the number of taxa detected per liter of water sampled was significantly higher in the eDNA samples. eDNA was especially useful for detecting delicate gelatinous animals which are undersampled by nets. We also detected eDNA changes in community composition with depth, but not with sample collection time (day vs. night). We provide recommendations for applying eDNA-based methods in the mesopelagic including the need for studies enabling interpretation of eDNA signals and improvement of barcode reference databases.

CARD-FISH in the Sequencing Era: Opening a New Universe of Protistan Ecology

Phagotrophic protists are key players in aquatic food webs. Although sequencing-based studies have revealed their enormous diversity, ecological information on in situ abundance, feeding modes, grazing preferences, and growth rates of specific lineages can be reliably obtained only using microscopy-based molecular methods, such as Catalyzed Reporter Deposition-Fluorescence in situ Hybridization (CARD-FISH). CARD-FISH is commonly applied to study prokaryotes, but less so to microbial eukaryotes. Application of this technique revealed that Paraphysomonas or Spumella-like chrysophytes, considered to be among the most prominent members of protistan communities in pelagic environments, are omnipresent but actually less abundant than expected, in contrast to little known groups such as heterotrophic cryptophyte lineages (e.g., CRY1), cercozoans, katablepharids, or the MAST lineages. Combination of CARD-FISH with tracer techniques and application of double CARD-FISH allow visualization of food vacuole contents of specific flagellate groups, thus considerably challenging our current, simplistic view that they are predominantly bacterivores. Experimental manipulations with natural communities revealed that larger flagellates are actually omnivores ingesting both prokaryotes and other protists. These new findings justify our proposition of an updated model of microbial food webs in pelagic environments, reflecting more authentically the complex trophic interactions and specific roles of flagellated protists, with inclusion of at least two additional trophic levels in the nanoplankton size fraction. Moreover, we provide a detailed CARD-FISH protocol for protists, exemplified on mixo- and heterotrophic nanoplanktonic flagellates, together with tips on probe design, a troubleshooting guide addressing most frequent obstacles, and an exhaustive list of published probes targeting protists.

A freshwater radiation of diplonemids

Diplonemids are considered marine protists and have been reported among the most abundant and diverse eukaryotes in the world oceans. Recently we detected the presence of freshwater diplonemids in Japanese deep freshwater lakes. However, their distribution and abundances in freshwater ecosystems remain unknown. We assessed abundance and diversity of diplonemids from several geographically distant deep freshwater lakes of the world by amplicon-sequencing, shotgun metagenomics and catalysed reporter deposition-fluorescent in situ hybridization (CARD-FISH). We found diplonemids in all the studied lakes, albeit with low abundances and diversity. We assembled long 18S rRNA sequences from freshwater diplonemids and showed that they form a new lineage distinct from the diverse marine clades. Freshwater diplonemids are a sister-group to a marine clade, which are mainly isolates from coastal and bay areas, suggesting a recent habitat transition from marine to freshwater habitats. Images of CARD-FISH targeted freshwater diplonemids suggest they feed on bacteria. Our analyses of 18S rRNA sequences retrieved from single-cell genomes of marine diplonemids show they encode multiple rRNA copies that may be very divergent from each other, suggesting that marine diplonemid abundance and diversity both have been overestimated. These results have wider implications on assessing eukaryotic abundances in natural habitats by using amplicon-sequencing alone.

Modeling the ecological status response of rivers to multiple stressors using machine learning: A comparison of environmental DNA metabarcoding and morphological data

Understanding the ecological status response of rivers to multiple stressors is a precondition for river restoration and management. However, this requires the collection of appropriate data, including environmental variables and the status of aquatic organisms, and analysis via a suitable model that captures the nonlinear relationships between ecological status and various stressors. The morphological approach has been the standard data collection method employed for establishing the status of aquatic organisms. However, this approach is very laborious and restricted to a specific set of organisms. Recently, an environmental DNA (eDNA) metabarcoding data approach has been developed that is far more efficient than the morphological approach and potentially applicable to an unlimited set of organisms. However, it remains unclear how well eDNA metabarcoding data reflects the impacts of environmental stressors on aquatic ecosystems compared with morphological data, which is essential for clarifying the potential applications of eDNA metabarcoding data in the ecological monitoring and management of rivers. The present work addresses this issue by modeling organism diversity based on three indices with respect to multiple environmental variables in both the catchment and reach scales. This is done by corresponding support vector machine (SVM) models constructed from eDNA metabarcoding and morphological data on 24 sampling locations in the Taizi River basin, China. According to the mean absolute percent error (MAPE) between the measured diversity index values and the index values predicted by the SVM models, the SVM models constructed from eDNA metabarcoding data (MAPE = 3.87) provide more accurate predictions than the SVM models constructed from morphological data (MAPE = 28.36), revealing that the eDNA metabarcoding data better reflects environmental conditions. In addition, the sensitivity of SVM model predictions of the ecological indices for both catchment-scale and reach-scale stressors is evaluated, and the stressors having the greatest impact on the ecological status of rivers are identified. The results demonstrate that the ecological status of rivers is more sensitive to environmental stressors at the reach scale than to stressors at the catchment scale. Therefore, our study is helpful in exploring the potential applications of eDNA metabarcoding data and SVM modeling in the ecological monitoring and management of rivers.

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.

Comparing Morphological and Molecular Estimates of Species Diversity in the Freshwater Genus Synura (Stramenopiles): A Model for Understanding Diversity of Eukaryotic Microorganisms

We performed a comparison of molecular and morphological diversity in a freshwater colonial genus Synura (Chrysophyceae, Stramenopiles), using the island of Newfoundland (Canada) as a case study. We examined the morphological species diversity in collections from 79 localities, and compared these findings to diversity based on molecular characters for 150 strains isolated from the same sites. Of 27 species or species-level lineages identified, only one third was recorded by both molecular and morphological techniques, showing both approaches are complementary in estimating species diversity within this genus. Eight taxa, each representing young evolutionary lineages, were recovered only by sequencing of isolated colonies, whereas ten species were recovered only microscopically. Our complex investigation, involving both morphological and molecular examinations, indicates that our knowledge of Synura diversity is still poor, limited only to a few well-studied areas. We revealed considerable cryptic diversity within the core S. petersenii and S. leptorrhabda lineages. We further resolved the phylogenetic position of two previously described taxa, S. kristiansenii and S. petersenii f. praefracta, propose species-level status for S. petersenii f. praefracta, and describe three new species, S. vinlandica, S. fluviatilis, and S. cornuta. Our findings add to the growing body of literature detailing distribution patterns observed in the genus, ranging from cosmopolitan species, to highly restricted taxa, to species such as S. hibernica found along coastal regions on multiple continents. Finally, our study illustrates the usefulness of combining detailed morphological information with gene sequence data to examine species diversity within chrysophyte algae.

Unraveling the Diversity of Eukaryotic Microplankton in a Large and Deep Perialpine Lake Using a High Throughput Sequencing Approach

The structure of microbial communities, microalgae, heterotrophic protozoa and fungi contributes to characterize food webs and productivity and, from an anthropogenic point of view, the qualitative characteristics of water bodies. Traditionally, in freshwater environments many investigations have been directed to the study of pelagic microalgae (“phytoplankton”) and periphyton (i.e., photosynthetic and mixotrophic protists) through the use of light microscopy (LM). While the number of studies on bacterioplankton communities have shown a substantial increase after the advent of high-throughput sequencing (HTS) approaches, the study of the composition, structure, and spatio-temporal patterns of microbial eukaryotes in freshwater environments was much less widespread. Moreover, the understanding of the correspondence between the relative phytoplankton abundances estimated by HTS and LM is still incomplete. Taking into account these limitations, this study examined the biodiversity and seasonality of the community of eukaryotic microplankton in the epilimnetic layer of a large and deep perialpine lake (Lake Garda) using HTS. The analyses were carried out at monthly frequency during 2014 and 2015. The results highlighted the existence of a rich and well diversified community and the presence of numerous phytoplankton taxa that were never identified by LM in previous investigations. Furthermore, the relative abundances of phytoplankton estimated by HTS and LM showed a significant relationship at different taxonomic ranks. In the 2 years of investigation, the temporal development of the whole micro-eukaryotic community showed a clear non-random and comparable distribution pattern, with the main taxonomic groups coherently distributed in the individual seasons. In perspective, the results obtained in this study highlight the importance of HTS approaches in assessing biodiversity and the relative importance of the main protist groups along environmental gradients, including those caused by anthropogenic impacts (e.g., eutrophication and climate change).

Rapid succession drives spring community dynamics of small protists at Helgoland Roads, North Sea

The dynamics of diatoms and dinoflagellates have been monitored for many decades at the Helgoland Roads Long-Term Ecological Research site and are relatively well understood. In contrast, small-sized eukaryotic microbes and their community changes are still much more elusive, mainly due to their small size and uniform morphology, which makes them difficult to identify microscopically. By using next-generation sequencing, we wanted to shed light on the Helgoland planktonic community dynamics, including nano- and picoplankton, during a spring bloom. We took samples from March to May 2016 and sequenced the V4 region of the 18S rDNA. Our results showed that mixotrophic and heterotrophic taxa were more abundant than autotrophic diatoms. Dinoflagellates dominated the sequence assemblage, and several small-sized eukaryotic microbes like Haptophyta, Choanoflagellata, Marine Stramenopiles and Syndiniales were identified. A diverse background community including taxa from all size classes was present during the whole sampling period. Five phases with several communities were distinguished. The fastest changes in community composition took place in phase 3, while the communities from phases 1 to 5 were more similar to each other despite contrasting environmental conditions. Synergy effects of next-generation sequencing and traditional methods may be exploited in future long-term observations.

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.

A Primer for Microbiome Time-Series Analysis

Time-series can provide critical insights into the structure and function of microbial communities. The analysis of temporal data warrants statistical considerations, distinct from comparative microbiome studies, to address ecological questions. This primer identifies unique challenges and approaches for analyzing microbiome time-series. In doing so, we focus on (1) identifying compositionally similar samples, (2) inferring putative interactions among populations, and (3) detecting periodic signals. We connect theory, code and data via a series of hands-on modules with a motivating biological question centered on marine microbial ecology. The topics of the modules include characterizing shifts in community structure and activity, identifying expression levels with a diel periodic signal, and identifying putative interactions within a complex community. Modules are presented as self-contained, open-access, interactive tutorials in R and Matlab. Throughout, we highlight statistical considerations for dealing with autocorrelated and compositional data, with an eye to improving the robustness of inferences from microbiome time-series. In doing so, we hope that this primer helps to broaden the use of time-series analytic methods within the microbial ecology research community.

Disentangling the mechanisms shaping the surface ocean microbiota

BACKGROUND

The ocean microbiota modulates global biogeochemical cycles and changes in its configuration may have large-scale consequences. Yet, the underlying ecological mechanisms structuring it are unclear. Here, we investigate how fundamental ecological mechanisms (selection, dispersal and ecological drift) shape the smallest members of the tropical and subtropical surface-ocean microbiota: prokaryotes and minute eukaryotes (picoeukaryotes). Furthermore, we investigate the agents exerting abiotic selection on this assemblage as well as the spatial patterns emerging from the action of ecological mechanisms. To explore this, we analysed the composition of surface-ocean prokaryotic and picoeukaryotic communities using DNA-sequence data (16S- and 18S-rRNA genes) collected during the circumglobal expeditions Malaspina-2010 and TARA-Oceans.

RESULTS

We found that the two main components of the tropical and subtropical surface-ocean microbiota, prokaryotes and picoeukaryotes, appear to be structured by different ecological mechanisms. Picoeukaryotic communities were predominantly structured by dispersal-limitation, while prokaryotic counterparts appeared to be shaped by the combined action of dispersal-limitation, selection and drift. Temperature-driven selection appeared as a major factor, out of a few selected factors, influencing species co-occurrence networks in prokaryotes but not in picoeukaryotes, indicating that association patterns may contribute to understand ocean microbiota structure and response to selection. Other measured abiotic variables seemed to have limited selective effects on community structure in the tropical and subtropical ocean. Picoeukaryotes displayed a higher spatial differentiation between communities and a higher distance decay when compared to prokaryotes, consistent with a scenario of higher dispersal limitation in the former after considering environmental heterogeneity. Lastly, random dynamics or drift seemed to have a more important role in structuring prokaryotic communities than picoeukaryotic counterparts.

CONCLUSIONS

The differential action of ecological mechanisms seems to cause contrasting biogeography, in the tropical and subtropical ocean, among the smallest surface plankton, prokaryotes and picoeukaryotes. This suggests that the idiosyncrasy of the main constituents of the ocean microbiota should be considered in order to understand its current and future configuration, which is especially relevant in a context of global change, where the reaction of surface ocean plankton to temperature increase is still unclear. Video Abstract.

Tracing the Origin of Planktonic Protists in an Ancient Lake

Ancient lakes are among the most interesting models for evolution studies because their biodiversity is the result of a complex combination of migration and speciation. Here, we investigate the origin of single celled planktonic eukaryotes from the oldest lake in the world—Lake Baikal (Russia). By using 18S rDNA metabarcoding, we recovered 1414 Operational Taxonomic Units (OTUs) belonging to protists populating surface waters (1–50 m) and representing pico/nano-sized cells. The recovered communities resembled other lacustrine freshwater assemblages found elsewhere, especially the taxonomically unclassified protists. However, our results suggest that a fraction of Baikal protists could belong to glacial relicts and have close relationships with marine/brackish species. Moreover, our results suggest that rapid radiation may have occurred among some protist taxa, partially mirroring what was already shown for multicellular organisms in Lake Baikal. We found 16% of the OTUs belonging to potential species flocks in Stramenopiles, Alveolata, Opisthokonta, Archaeplastida, Rhizaria, and Hacrobia. Putative flocks predominated in Chrysophytes, which are highly diverse in Lake Baikal. Also, the 18S rDNA of a number of species (7% of the total) differed >10% from other known sequences. These taxa as well as those belonging to the flocks may be endemic to Lake Baikal. Overall, our study points to novel diversity of planktonic protists in Lake Baikal, some of which may have emerged in situ after evolutionary diversification.

Factors shaping community patterns of protists and bacteria on a European scale

Factors shaping community patterns of microorganisms are controversially discussed. Physical and chemical factors certainly limit the survival of individual taxa and maintenance of diversity. In recent years, a contribution of geographic distance and dispersal barriers to distribution patterns of protists and bacteria has been demonstrated. Organismic interactions such as competition, predation and mutualism further modify community structure and maintenance of distinct taxa. Here, we address the relative importance of these different factors in shaping protists and bacterial communities on a European scale using high-throughput sequencing data obtained from lentic freshwater ecosystems. We show that community patterns of protists are similar to those of bacteria. Our results indicate that cross-domain organismic factors are important variables with a higher influence on protists as compared with bacteria. Abiotic physical and chemical factors also contributed significantly to community patterns. The contribution of these latter factors was higher for bacteria, which may reflect a stronger biogeochemical coupling. The contribution of geographical distance was similar for both microbial groups.

Cryptic species in the parasitic Amoebophrya species complex revealed by a polyphasic approach

As critical primary producers and recyclers of organic matter, the diversity of marine protists has been extensively explored by high-throughput barcode sequencing. However, classification of short metabarcoding sequences into traditional taxonomic units is not trivial, especially for lineages mainly known by their genetic fingerprints. This is the case for the widespread Amoebophrya ceratii species complex, parasites of their dinoflagellate congeners. We used genetic and phenotypic characters, applied to 119 Amoebophrya individuals sampled from the same geographic area, to construct practical guidelines for species delineation that could be applied in DNA/RNA based diversity analyses. Based on the internal transcribed spacer (ITS) regions, ITS2 compensatory base changes (CBC) and genome k-mer comparisons, we unambiguously defined eight cryptic species among closely related ribotypes that differed by less than 97% sequence identity in their SSU rDNA. We then followed the genetic signatures of these parasitic species during a three-year survey of Alexandrium minutum blooms. We showed that these cryptic Amoebophrya species co-occurred and shared the same ecological niche. We also observed a maximal ecological fitness for parasites having narrow to intermediate host ranges, reflecting a high cost for infecting a broader host range. This study suggests that a complete taxonomic revision of these parasitic dinoflagellates is long overdue to understand their diversity and ecological role in the marine plankton.

The global-scale distributions of soil protists and their contributions to belowground systems

Protists are ubiquitous in soil, where they are key contributors to nutrient cycling and energy transfer. However, protists have received far less attention than other components of the soil microbiome. We used amplicon sequencing of soils from 180 locations across six continents to investigate the ecological preferences of protists and their functional contributions to belowground systems. We complemented these analyses with shotgun metagenomic sequencing of 46 soils to validate the identities of the more abundant protist lineages. We found that most soils are dominated by consumers, although parasites and phototrophs are particularly abundant in tropical and arid ecosystems, respectively. The best predictors of protist composition (primarily annual precipitation) are fundamentally distinct from those shaping bacterial and archaeal communities (namely, soil pH). Some protists and bacteria co-occur globally, highlighting the potential importance of these largely undescribed belowground interactions. Together, this study allowed us to identify the most abundant and ubiquitous protists living in soil, with our work providing a cross-ecosystem perspective on the factors structuring soil protist communities and their likely contributions to soil functioning.

Unassigned diversity of planktonic foraminifera from environmental sequencing revealed as known but neglected species

Most research on extant planktonic foraminifera has been directed towards larger species (>0.150 mm) which can be easily manipulated, counted and yield enough calcite for geochemical analyses. This has drawn attention towards the macroperforate clade and created an impression of their numerical and ecological dominance. Drawing such conclusions from the study of such “giants” is a dangerous path. There were times in the evolutionary history of planktonic foraminifera when all species were smaller than 0.1 mm and indeed numerous small taxa, mainly from the microperforate clade, have been formally described from the modern plankton. The significance of these small, obscure and neglected species is poorly characterized and their relationship to the newly discovered hyperabundant but uncharacterized lineages of planktonic foraminifera in metabarcoding datasets is unknown. To determine, who is hiding in the metabarcoding datasets, we carried out an extensive sequencing of 18S rDNA targeted at small and obscure species. The sequences of the newly characterized small and obscure taxa match many of the previously uncharacterized lineages found in metabarcoding data. This indicates that most of the modern diversity in planktonic foraminifera has been taxonomically captured, but the role of the small and neglected taxa has been severely underestimated.