Newly designed foraminifera primers identify habitat-specific lineages through metabarcoding analyses.

Foraminifera include diverse shell-building lineages found in a wide array of aquatic habitats from the deep-sea to intertidal zones to brackish and freshwater ecosystems. Recent estimates of morphological and molecular foraminifera diversity have increased the knowledge of foraminiferal diversity, which is critical as these lineages are used as bioindicators of past and present environmental perturbation. However, a comparative analysis of foraminiferal biodiversity between their major habitats (freshwater, brackish, intertidal, and marine) is underexplored, particularly using molecular tools. Here, we present a metabarcoding survey of foraminiferal diversity across different ecosystems using newly designed foraminifera-specific primers that target the hypervariable regions of the foraminifera SSU-rRNA gene (~250-300 bp long). We tested these primer sets on four foraminifera species and then across several environments: the intertidal zone, coastal ecosystems, and freshwater vernal pools. We retrieved 655 operational taxonomic units (OTUs); the majority of which are undetermined taxa that have no closely matching sequences in the reference database. Furthermore, we identified 163 OTUs with distinct habitat preferences. Most of the observed OTUs belonged to lineages of single-chambered foraminifera, including poorly explored freshwater foraminifera which encompass a clade of Reticulomyxa-like forms. Our pilot study provides the community with an additional set of newly designed and taxon-specific primers to elucidate foraminiferal diversity across different habitats.

PhyloToL: A Taxon/Gene-Rich Phylogenomic Pipeline to Explore Genome Evolution of Diverse Eukaryotes.

Estimating multiple sequence alignments (MSAs) and inferring phylogenies are essential for many aspects of comparative biology. Yet, many bioinformatics tools for such analyses have focused on specific clades, with greatest attention paid to plants, animals, and fungi. The rapid increase in high-throughput sequencing (HTS) data from diverse lineages now provides opportunities to estimate evolutionary relationships and gene family evolution across the eukaryotic tree of life. At the same time, these types of data are known to be error-prone (e.g., substitutions, contamination). To address these opportunities and challenges, we have refined a phylogenomic pipeline, now named PhyloToL, to allow easy incorporation of data from HTS studies, to automate production of both MSAs and gene trees, and to identify and remove contaminants. PhyloToL is designed for phylogenomic analyses of diverse lineages across the tree of life (i.e., at scales of >100 My). We demonstrate the power of PhyloToL by assessing stop codon usage in Ciliophora, identifying contamination in a taxon- and gene-rich database and exploring the evolutionary history of chromosomes in the kinetoplastid parasite Trypanosoma brucei, the causative agent of African sleeping sickness. Benchmarking PhyloToL's homology assessment against that of OrthoMCL and a published paper on superfamilies of bacterial and eukaryotic organellar outer membrane pore-forming proteins demonstrates the power of our approach for determining gene family membership and inferring gene trees. PhyloToL is highly flexible and allows users to easily explore HTS data, test hypotheses about phylogeny and gene family evolution and combine outputs with third-party tools (e.g., PhyloChromoMap, iGTP).

Ongoing speciation within the diatom Fragilariopsis kerguelensis in the Southern Ocean?

While the identification of microbial eukaryotes using metabarcoding tools is now widespread, additional data are needed to confirm molecular observations, to mark the difference between species and population variants, and to better understand the biogeography of microbial eukaryotes. In this issue of Molecular Ecology, Postel et al not only use three molecular approaches to identify subgroups of Fragilariopsis kerguelensis but also use morphology and physiology to better understand the relationship between the three genotypes. They revealed that (a) the three genotypes of the diatom F. kerguelensis show negligible gene flux; (b) two of the genotypes are geographically isolated with different physiology but still able to crossbreed; and (c) the remaining genotype is omnipresent but reproductively isolated.

High Diversity of Testate Amoebae (Amoebozoa, Arcellinida) Detected by HTS Analyses in a New England Fen using Newly Designed Taxon-specific Primers.

Testate (shell-building) amoebae, such as the Arcellinida (Amoebozoa), are useful bioindicators for climate change. Though past work has relied on morphological analyses to characterize Arcellinida diversity, genetic analyses revealed the presence of multiple cryptic species underlying morphospecies. Here, we design and deploy Arcellinida-specific primers for the SSU-rDNA gene to assess the community composition on the molecular level in a pilot study of two samplings from a New England fen: (1) 36-cm horizontal transects and vertical cores; and (2) 26-m horizontal transects fractioned into four size classes (2-10, 10-35, 35-100, and 100-300 µm). Analyses of these data show the following: (1) a considerable genetic diversity within Arcellinida, much of which comes from morphospecies lacking sequences on GenBank; (2) communities characterized by DNA (i.e. active + quiescent) are distinct from those characterized by RNA (i.e. active, indicator of biomass); (3) active communities on the surface tend to be more similar to one another than to core communities, despite considerable heterogeneity; and (4) analyses of communities fractioned by size find some lineages (OTUs) that are abundant in disjunct size categories, suggesting the possibility of life-history stages. Together, these data demonstrate the potential of these primers to elucidate the diversity of Arcellinida communities in diverse habitats.

Microbial Diversity in the Eukaryotic SAR Clade: Illuminating the Darkness Between Morphology and Molecular Data.

Despite their diversity and ecological importance, many areas of the SAR-Stramenopila, Alveolata, and Rhizaria-clade are poorly understood as the majority (90%) of SAR species lack molecular data and only 5% of species are from well-sampled families. Here, we review and summarize the state of knowledge about the three major clades of SAR, describing the diversity within each clade and identifying synapomorphies when possible. We also assess the "dark area" of SAR: the morphologically described species that are missing molecular data. The majority of molecular data for SAR lineages are characterized from marine samples and vertebrate hosts, highlighting the need for additional research effort in areas such as freshwater and terrestrial habitats and "non-vertebrate" hosts. We also describe the paucity of data on the biogeography of SAR species, and point to opportunities to illuminate diversity in this major eukaryotic clade. See also the video abstract here: https://youtu.be/_VUXqaX19Rw.

Unexpected biodiversity of ciliates in marine samples from below the photic zone.

Marine microbial eukaryotes play critical roles in planktonic food webs and have been described as most diverse in the photic zone where productivity is high. We used high-throughput sequencing (HTS) to analyse the spatial distribution of planktonic ciliate diversity from shallow waters (<30 m depth) to beyond the continental shelf (>800 m depth) along a 163 km transect off the coast of New England, USA. We focus on ciliates in the subclasses Oligotrichia and Choreotrichia (class Spirotrichea), as these taxa are major components of marine food webs. We did not observe the decrease of diversity below the photic zone expected based on productivity and previous analyses. Instead, we saw an increase of diversity with depth. We also observed that the ciliate communities assessed by HTS cluster by depth layer and degree of water column stratification, suggesting that community assembly is driven by environmental factors. Across our samples, abundant OTUs tend to match previously characterized morphospecies while rare OTUs are more often undescribed, consistent with the idea that species in the rare biosphere remain to be characterized by microscopy. Finally, samples taken below the photic zone also reveal the prevalence of two uncharacterized (i.e. lacking sequenced morphospecies) clades - clusters X1 and X2 - that are enriched within the nano-sized fraction (2-10 µm) and are defined by deletions within the region of the SSU-rDNA analysed here. Together, these data reinforce that we still have much to learn about microbial diversity in marine ecosystems, especially in deep-waters that may be a reservoir for rare species and uncharacterized taxa.

Pyrosequencing for assessing diversity of eukaryotic microbes: analysis of data on marine planktonic ciliates and comparison with traditional methods.

Assessing microbial diversity requires analysis of all three domains of life, including eukaryotic microbes. We examined the diversity of two ecologically important clades of microbial eukaryotes, ciliates in the subclasses Oligotrichia and Choreotrichia (class Spirotrichea), by comparing pyrosequencing to Sanger-sequenced clone libraries and microscopy. Using samples from a large temperate estuary (Long Island Sound, USA), we gained three major insights. First, richness estimates varied by up to one order of magnitude either using different criteria for pyrosequence processing or among pyrosequencing, cloning and microscopy, while taxon identification was almost always coherent. Error-correcting algorithms for pyrosequences ('denoising') reduced discrepancies in richness but also removed known morphospecies from the data. Second, although most of the pyrosequenced operational taxonomic units (OTUs) were distributed within known orders and families, we found evidence of a previously uncharacterized or unknown clade even in these ciliate lineages that have a rich history of morphological description. Third, pyrosequencing allowed the detection of OTUs that were either dominant or extremely rare in different samples. Our findings confirm the potential of pyrosequencing for quantifying microbial diversity, but also highlight the importance of careful evaluation of pyrosequence processing for using this method to address ecological questions.

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.

Incubation and grazing effects on spirotrich ciliate diversity inferred from molecular analyses of microcosm experiments.

We used an experimental approach of analyzing marine microcosms to evaluate the impact of both predation (top-down) and food resources (bottom-up) on spirotrich ciliate communities. To assess the diversity, we used two molecular methods-denaturing gradient gel electrophoresis (DGGE) and high-throughput sequencing (HTS). We carried out two types of experiments to measure top-down (adult copepods as predators) and bottom-up effects (phytoplankton as food resources) on the spirotrich ciliates. We observed both strong incubation effects (untreated controls departed from initial assessment of diversity) and high variability across replicates within treatments, particularly for the bottom-up experiments. This suggests a rapid community turn-over during incubation and differential susceptibility to the effects of experimental manipulation. Despite the variability, our analyses reveal some broad patterns such as (1) increasing adult copepod predator abundance had a greater impact on spirotrich ciliates than on other microbial eukaryotes; (2) there was no evidence for strong food selection by the dominant spirotrich ciliates.

Distribution of Abundant and Active Planktonic Ciliates in Coastal and Slope Waters Off New England.

Despite their important role of linking microbial and classic marine food webs, data on biogeographical patterns of microbial eukaryotic grazers are limited, and even fewer studies have used molecular tools to assess active (i.e., those expressing genes) community members. Marine ciliate diversity is believed to be greatest at the chlorophyll maximum, where there is an abundance of autotrophic prey, and is often assumed to decline with depth. Here, we assess the abundant (DNA) and active (RNA) marine ciliate communities throughout the water column at two stations off the New England coast (Northwest Atlantic)-a coastal station 43 km from shore (40 m depth) and a slope station 135 km off shore (1,000 m). We analyze ciliate communities using a DNA fingerprinting technique, Denaturing Gradient Gel Electrophoresis (DGGE), which captures patterns of abundant community members. We compare estimates of ciliate communities from SSU-rDNA (abundant) and SSU-rRNA (active) and find complex patterns throughout the water column, including many active lineages below the photic zone. Our analyses reveal (1) a number of widely-distributed taxa that are both abundant and active; (2) considerable heterogeneity in patterns of presence/absence of taxa in offshore samples taken 50 m apart throughout the water column; and (3) three distinct ciliate assemblages based on position from shore and depth. Analysis of active (RNA) taxa uncovers biodiversity hidden to traditional DNA-based approaches (e.g., clone library, rDNA amplicon studies).

Patchiness of Ciliate Communities Sampled at Varying Spatial Scales along the New England Shelf.

Although protists (microbial eukaryotes) provide an important link between bacteria and Metazoa in food webs, we do not yet have a clear understanding of the spatial scales on which protist diversity varies. Here, we use a combination of DNA fingerprinting (denaturant gradient gel electrophoresis or DGGE) and high-throughput sequencing (HTS) to assess the ciliate community in the class Spirotrichea at varying scales of 1-3 km sampled in three locations separated by at least 25 km-offshore, midshelf and inshore-along the New England shelf. Analyses of both abundant community (DGGE) and the total community (HTS) members reveal that: 1) ciliate communities are patchily distributed inshore (i.e. the middle station of a transect is distinct from its two neighboring stations), whereas communities are more homogeneous among samples within the midshelf and offshore stations; 2) a ciliate closely related to Pelagostrobilidium paraepacrum 'blooms' inshore and; 3) environmental factors may differentially impact the distributions of individual ciliates (i.e. OTUs) rather than the community as a whole as OTUs tend to show distinct biogeographies (e.g. some OTUs are restricted to the offshore locations, some to the surface, etc.). Together, these data show the complexity underlying the spatial distributions of marine protists, and suggest that biogeography may be a property of ciliate species rather than communities.

Patterns and processes in microbial biogeography: do molecules and morphologies give the same answers?

Our knowledge on microbial biogeography depends on the way we define and study diversity. In contrast to most microbes, some protist lineages have conspicuous structures that allow comparisons of diversity concepts and measures-those based on molecules and those based on morphology. We analyzed a group of shell-bearing planktonic ciliates, the tintinnids, in a coast-to-ocean gradient using high-throughput sequencing and microscopy. First, we compared molecular operational taxonomic units (OTUs) and morphospecies in terms of assemblage composition, distribution and relationships with the environment. OTUs revealed potentially novel and rare taxa, while morphospecies showed clearer correlations with environmental factors, and both approaches coincided in supporting a coastal versus oceanic pattern. Second, we explored which processes influence assembly across the environmental gradient examined. Assemblage fluctuations were associated with significant distance-decay and changes in morphospecies size and prey proxies, thus suggesting niche partitioning as a key structuring mechanism. Our conclusion is that molecules and morphologies generally agreed, but they provided complementary data, the first revealing hidden diversity, and the latter making better connections between distribution patterns and ecological processes. This highlights the importance of linking genotypes and phenotypes (using multidisciplinary analyses and/or reliable databases of barcoded species), to understand the diversity, biogeography and ecological roles of microbes.

Diversity of diversity: conceptual and methodological differences in biodiversity estimates of eukaryotic microbes as compared to bacteria.

Recent advances such as high-throughput sequencing (HTS) have changed conceptions about the magnitude of diversity on Earth. This is especially true for microbial lineages, which have seen the discovery of great numbers of rare forms in places such as the human gut as well as diverse environments (e.g., freshwater, marine, and soil). Given the differences in perceptions of diversity for bacterial and eukaryotic microbes, including divergent species concepts, HTS tools used to eliminate errors and population-level variation in bacteria may not be appropriate for microbial eukaryotes and may eliminate valid species from the data. We discuss here how the nature of biodiversity varies among microbial groups and the extent to which HTS tools designed for bacteria are useful for eukaryotes.

Microplanktonic community structure in a coastal system relative to a Phaeocystis bloom inferred from morphological and tag pyrosequencing methods.

BACKGROUND: Massive phytoplankton blooms, like the recurrent Phaeocystis proliferation observed every year in the Eastern English Channel (EEC), have a significant influence on the overall planktonic community structure and their food web dynamics. As well as being an important area for local fisheries, the EEC is an ideal ecosystem for work on microbial diversity. This is because, although its environmental context is relatively complex, it is reasonably well understood due to several years of monitoring and morphological observations of its planktonic organisms. The objective of our study was to better understand the under-explored microbial eukaryotic diversity relative to the Phaeocystis bloom. METHODOLOGY AND PRINCIPAL FINDINGS: The community structure of microplankton (diatoms, haptophytes, ciliates and dinoflagellates) was studied through morphological observations and tag pyrosequencing. During the annual Phaeocystis spring bloom, the phytoplankton biomass increased by 34-fold, while the microzooplankton biomass showed a 4-fold increase, representing on average about 4.6% of the biomass of their phytoplankton prey. Tag pyrosequencing unveiled an extensive diversity of Gymnodiniaceae, with G. spirale and G. fusiformis representing the most abundant reads. An extended diversity of Phaeocystales, with partial 18S rDNA genes sequence identity as low as 85% was found, with taxa corresponding to P. globosa, but also to unknown Phaeocystaceae. CONCLUSIONS: Morphological analyses and pyrosequencing were generally in accordance with capturing frequency shifts of abundant taxa. Tag pyrosequencing allowed highlighting the maintenance of microplankton diversity during the Phaeocystis bloom and the increase of the taxa presenting low number of reads (minor taxa) along with the dominant ones in response to biotic and/or abiotic changing conditions. Although molecular approaches have enhanced our perception on diversity, it has come to light that the challenge of modelling and predicting ecological change requires the use of different complementary approaches, to link taxonomic data with the functional roles of microbes in biogeochemical cycles.