General patterns of diversity in major marine microeukaryote lineages

A new sampling device for metabarcoding surveillance of port communities and detection of non-indigenous species

Metabarcoding techniques are revolutionizing studies of marine biodiversity. They can be used for monitoring non-indigenous species (NIS) in ports and harbors. However, they are often biased by inconsistent sampling methods and incomplete reference databases. Logistic constraints in ports prompt the development of simple, easy-to-deploy samplers. We tested a new device called polyamide mesh for ports organismal monitoring (POMPOM) with a high surface-to-volume ratio. POMPOMS were deployed inside a fishing and recreational port in the Mediterranean alongside conventional settlement plates. We also compiled a curated database with cytochrome oxidase (COI) sequences of Mediterranean NIS. COI metabarcoding of the communities settled in the POMPOMs captured a similar biodiversity than settlement plates, with shared molecular operational units (MOTUs) representing ca. 99% of reads. 38 NIS were detected in the port accounting for ca. 26% of reads. POMPOMs were easy to deploy and handle and provide an efficient method for NIS surveillance.

Inputs of seabird guano alter microbial growth, community composition and the phytoplankton-bacterial interactions in a coastal system

Seabird guano enters coastal waters providing bioavailable substrates for microbial plankton, but their role in marine ecosystem functioning remains poorly understood. Two concentrations of the water soluble fraction (WSF) of gull guano were added to different natural microbial communities collected in surface waters from the Ría de Vigo (NW Spain) in spring, summer, and winter. Samples were incubated with or without antibiotics (to block bacterial activity) to test whether gull guano stimulated phytoplankton and bacterial growth, caused changes in taxonomic composition, and altered phytoplankton-bacteria interactions. Alteromonadales, Sphingobacteriales, Verrucomicrobia and diatoms were generally stimulated by guano. Chlorophyll a (Chl a) concentration and bacterial abundance significantly increased after additions independently of the initial ambient nutrient concentrations. Our study demonstrates, for the first time, that the addition of guano altered the phytoplankton-bacteria interaction index from neutral (i.e. phytoplankton growth was not affected by bacterial activity) to positive (i.e. phytoplankton growth was stimulated by bacterial activity) in the low-nutrient environment occurring in spring. In contrast, when environmental nutrient concentrations were high, the interaction index changed from positive to neutral after guano additions, suggesting the presence of some secondary metabolite in the guano that is needed for phytoplankton growth, which would otherwise be supplied by bacteria.

Parasitic infections by Group II Syndiniales target selected dinoflagellate host populations within diverse protist assemblages in a model coastal pond

Protists are integral to marine food webs and biogeochemical cycles; however, there is a paucity of data describing specific ecological niches for some of the most abundant taxa in marker gene libraries. Syndiniales are one such group, often representing the majority of sequence reads recovered from picoplankton samples across the global ocean. However, the prevalence and impacts of syndinian parasitism in marine environments remains unclear. We began to address these critical knowledge gaps by generating a high-resolution time series (March-October 2018) in a productive coastal pond. Seasonal shifts in protist populations, including parasitic Syndiniales, were documented during periods of higher primary productivity and increased summer temperature-driven stratification. Elevated concentrations of infected hosts and free-living parasite spores occurred at nearly monthly intervals in July, August, and September. We suggest intensifying stratification during this period correlated with the increased prevalence of dinoflagellates that were parasitized by Group II Syndiniales. Infections in some protist populations were comparable to previously reported large single-taxon dinoflagellate blooms. Infection dynamics in Salt Pond demonstrated the propagation of syndinian parasites through mixed protist assemblages and highlighted patterns of host/parasite interactions that better reflect many other marine environments where single taxon blooms are uncommon.

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.

The dinoflagellate Alexandrium catenella producing only carbamate toxins may account for the seafood poisonings in Qinhuangdao, China

An outbreak of paralytic shellfish poisoning, recorded in April 2016 in Qinhuangdao China, was suspected to be caused by a toxic species in genus Alexandrium. Shortly after the poisoning outbreak, shellfish and net-concentrated phytoplankton samples were collected from the Bohai Sea, and analysed using high performance liquid chromatography coupled with fluorescence detection. Paralytic shellfish toxins (PSTs) were detected in both phytoplankton and shellfish samples, with similar toxin profiles dominated by carbamate toxins. High throughput sequencing data for phytoplankton samples collected previously in the coastal waters of Qinhuangdao were then analysed, and 8 operational taxonomic units (OTUs) were assigned to Alexandrium affine, A. andersonii/A. ostenfeldii, A. catenella, A. fraterculus, A. hiranoi/A. pseudogonyaulax, A. margalefii, A. pacificum and A. pohangense, among which A. catenella, A. pacificum and A. ostenfeldii could be potential producers of PSTs. During a cruise in 2019, three isolates of Alexandrium were established by cyst germination, and identified as A. catenella based on the sequences of the 28S ribosomal RNA gene (28S rDNA) D1-D2 region. Interestingly, all the three strains had the same toxin profile consisting of gonyautoxins 1, 3, 4 (GTX1, 3, 4) and neosaxitoxin (NEO). The toxin profile is similar to those of phytoplankton samples collected previously in the coastal waters of Qinhuangdao, but remarkably different from the general toxin profile of A. catenella dominated by N-sulfocarbamoyl toxins C1-2 in the Bohai Sea and the Yellow Sea. The results suggest that A. catenella is most likely to be the causative species of the poisoning outbreak in Qinhuangdao. As far as we know, this is the first report of A. catenella in the Bohai Sea producing PSTs dominated by high potent gonyautoxins GTX1-4. Occurrence of the highly toxic A. catenella will increase the risk of paralytic shellfish poisoning, which necessitates in-depth mechanism studies and increasing monitoring efforts.

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.

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.

Marked changes in diversity and relative activity of picoeukaryotes with depth in the world ocean

Microbial eukaryotes are key components of the ocean plankton. Yet, our understanding of their community composition and activity in different water layers of the ocean is limited, particularly for picoeukaryotes (0.2–3 µm cell size). Here, we examined the picoeukaryotic communities inhabiting different vertical zones of the tropical and subtropical global ocean: surface, deep chlorophyll maximum, mesopelagic (including the deep scattering layer and oxygen minimum zones), and bathypelagic. Communities were analysed by high-tthroughput sequencing of the 18S rRNA gene (V4 region) as represented by DNA (community structure) and RNA (metabolism), followed by delineation of Operational Taxonomic Units (OTUs) at 99% similarity. We found a stratification of the picoeukaryotic communities along the water column, with assemblages corresponding to the sunlit and dark ocean. Specific taxonomic groups either increased (e.g., Chrysophyceae or Bicosoecida) or decreased (e.g., Dinoflagellata or MAST-3) in abundance with depth. We used the rRNA:rDNA ratio of each OTU as a proxy of metabolic activity. The highest relative activity was found in the mesopelagic layer for most taxonomic groups, and the lowest in the bathypelagic. Altogether, we characterize the change in community structure and metabolic activity of picoeukaryotes with depth in the global ocean, suggesting a hotspot of activity in the mesopelagic.

Copper Affects Composition and Functioning of Microbial Communities in Marine Biofilms at Environmentally Relevant Concentrations

Copper (Cu) pollution in coastal areas is a worldwide threat for aquatic communities. This study aims to demonstrate the usefulness of the DNA metabarcoding analysis in order to describe the ecotoxicological effect of Cu at environmental concentrations on marine periphyton. Additionally, the study investigates if Cu-induced changes in community structure co-occurs with changes in community functioning (i.e., photosynthesis and community tolerance to Cu). Periphyton was exposed for 18 days to five Cu concentrations, between 0.01 and 10 μM, in a semi-static test. Diversity and community structure of prokaryotic and eukaryotic organisms were assessed by 16S and 18S amplicon sequencing, respectively. Community function was studied as impacts on algal biomass and photosynthetic activity. Additionally, we studied Pollution-Induced Community Tolerance (PICT) using photosynthesis as the endpoint. Sequencing results detected an average of 9,504 and 1,242 OTUs for 16S and 18S, respectively, reflecting the high biodiversity of marine periphytic biofilms. Eukaryotes represent the most Cu-sensitive kingdom, where effects were seen already at concentrations as low as 0.01 μM. The structure of the prokaryotic part of the community was impacted at slightly higher concentrations (0.06 μM), which is still in the range of the Cu concentrations observed in the area (0.08 μM). The current environmental quality standard for Cu of 0.07 μM therefore does not seem to be sufficiently protective for periphyton. Cu exposure resulted in a more Cu-tolerant community, which was accompanied by a reduced total algal biomass, increased relative abundance of diatoms and a reduction of photosynthetic activity. Cu exposure changed the network of associations between taxa in the communities. A total of 23 taxa, including taxa within Proteobacteria, Bacteroidetes, Stramenopiles, and Hacrobia, were identified as being particularly sensitive to Cu. DNA metabarcoding is presented as a sensitive tool for community-level ecotoxicological studies that allows to observe impacts simultaneously on a multitude of pro- and eukaryotic taxa, and therefore to identify particularly sensitive, non-cultivable taxa.

Contrasting prevalence of selection and drift in the community structuring of bacteria and microbial eukaryotes

Whether or not communities of microbial eukaryotes are structured in the same way as bacteria is a general and poorly explored question in ecology. Here, we investigated this question in a set of planktonic lake microbiotas in Eastern Antarctica that represent a natural community ecology experiment. Most of the analysed lakes emerged from the sea during the last 6000 years, giving rise to waterbodies that originally contained marine microbiotas and that subsequently evolved into habitats ranging from freshwater to hypersaline. We show that habitat diversification has promoted selection driven by the salinity gradient in bacterial communities (explaining ∼ 72% of taxa turnover), while microeukaryotic counterparts were predominantly structured by ecological drift (∼72% of the turnover). Nevertheless, we also detected a number of microeukaryotes with specific responses to salinity, indicating that albeit minor, selection has had a role in the structuring of specific members of their communities. In sum, we conclude that microeukaryotes and bacteria inhabiting the same communities can be structured predominantly by different processes. This should be considered in future studies aiming to understand the mechanisms that shape microbial assemblages.

Consistent patterns of high alpha and low beta diversity in tropical parasitic and free-living protists

Tropical animals and plants are known to have high alpha diversity within forests, but low beta diversity between forests. By contrast, it is unknown whether microbes inhabiting the same ecosystems exhibit similar biogeographic patterns. To evaluate the biogeographies of tropical protists, we used metabarcoding data of species sampled in the soils of three lowland Neotropical rainforests. Taxa-area and distance-decay relationships for three of the dominant protist taxa and their subtaxa were estimated at both the OTU and phylogenetic levels, with presence-absence and abundance-based measures. These estimates were compared to null models. High local alpha and low regional beta diversity patterns were consistently found for both the parasitic Apicomplexa and the largely free-living Cercozoa and Ciliophora. Similar to animals and plants, the protists showed spatial structures between forests at the OTU and phylogenetic levels, and only at the phylogenetic level within forests. These results suggest that the biogeographies of macro- and micro-organismal eukaryotes in lowland Neotropical rainforests are partially structured by the same general processes. However, and unlike the animals and plants, the protist OTUs did not exhibit spatial structures within forests, which hinders our ability to estimate the local and regional diversity of protists in tropical forests.

Surface ocean metabarcoding confirms limited diversity in planktonic foraminifera but reveals unknown hyper-abundant lineages

Since the advent of DNA metabarcoding surveys, the planktonic realm is considered a treasure trove of diversity, inhabited by a small number of abundant taxa, and a hugely diverse and taxonomically uncharacterized consortium of rare species. Here we assess if the apparent underestimation of plankton diversity applies universally. We target planktonic foraminifera, a group of protists whose known morphological diversity is limited, taxonomically resolved and linked to ribosomal DNA barcodes. We generated a pyrosequencing dataset of ~100,000 partial 18S rRNA foraminiferal sequences from 32 size fractioned photic-zone plankton samples collected at 8 stations in the Indian and Atlantic Oceans during the Tara Oceans expedition (2009–2012). We identified 69 genetic types belonging to 41 morphotaxa in our metabarcoding dataset. The diversity saturated at local and regional scale as well as in the three size fractions and the two depths sampled indicating that the diversity of foraminifera is modest and finite. The large majority of the newly discovered lineages occur in the small size fraction, neglected by classical taxonomy. These unknown lineages dominate the bulk [>0.8 µm] size fraction, implying that a considerable part of the planktonic foraminifera community biomass has its origin in unknown lineages.

Large-scale ocean connectivity and planktonic body size

Global patterns of planktonic diversity are mainly determined by the dispersal of propagules with ocean currents. However, the role that abundance and body size play in determining spatial patterns of diversity remains unclear. Here we analyse spatial community structure - β-diversity - for several planktonic and nektonic organisms from prokaryotes to small mesopelagic fishes collected during the Malaspina 2010 Expedition. β-diversity was compared to surface ocean transit times derived from a global circulation model, revealing a significant negative relationship that is stronger than environmental differences. Estimated dispersal scales for different groups show a negative correlation with body size, where less abundant large-bodied communities have significantly shorter dispersal scales and larger species spatial turnover rates than more abundant small-bodied plankton. Our results confirm that the dispersal scale of planktonic and micro-nektonic organisms is determined by local abundance, which scales with body size, ultimately setting global spatial patterns of diversity.

Body size is hypothesised to be a major determinant of β-diversity in passively-dispersing marine organisms. Here, Villarino et al. show that plankton body size determines rates of dispersal along marine currents, with shorter dispersal and higher species spatial turnover in larger organisms.

The diversity and biogeography of abundant and rare intertidal marine microeukaryotes explained by environment and dispersal limitation

Benthic microeukaryotes are key ecosystem drivers in marine sandy beaches, an important and dynamic environment; however, little is known about their diversity and biogeography on a large spatial scale. Here, we investigated the community composition and geographical distributions of benthic microeukaryotes using high-throughput sequencing of the 18S rRNA gene and quantified the contributions of environmental factors and spatial separation on the distribution patterns of both rare and abundant taxa. We collected 36 intertidal samples at 12 sandy beaches from four regions that spanned distances from 0.001 to 12,000 km. We found 12,890 operational taxonomic units (OTUs; 97% sequence identity level) including members of all eukaryotic super-groups and several phyla of uncertain position. Arthropoda and Diatomeae dominated the sequence reads in abundance, but Ciliophora and Discoba were the most diverse groups across all samples. About one-third of the OTUs could not be definitively classified at a similarity level of 80%, supporting the view that a large number of rare and minute marine species may have escaped previous characterization. We found generally similar geographical patterns for abundant and rare microeukaryotic sub-communities, and both showed a significant distance-decay similarity trend. Variation partitioning showed that both rare and abundant sub-communities exhibited a slightly stronger response to environmental factors than spatial (distance) factors. However, the abundant sub-community was strongly correlated with variations in spatial, environmental and sediment grain size factors (66% of variance explained), but the rare assemblage was not (16%). This suggests that different or more complex mechanisms generate and maintain diversity in the rare biosphere in this habitat.

Vertical distribution of major photosynthetic picoeukaryotic groups in stratified marine waters

Photosynthetic picoeukaryotes (PPEs) are fundamental contributors to oceanic primary production and form diverse communities dominated by prymnesiophytes, chlorophytes, pelagophytes and chrysophytes. Here, we studied the vertical distribution of these major groups in two offshore regions of the northern Iberian Peninsula during summer stratification. We performed a fine-scale vertical sampling (every ∼2 m) across the DCM and used fluorescence in situ hybridization (FISH) to determine the PPE composition and to explore the possible segregation of target groups in the light, nutrient and temperature gradients. Chlorophytes, pelagophytes and prymnesiophytes, in this order of abundance, accounted for the total PPEs recorded by flow cytometry in the Avilés canyon, and for more than half in the Galicia Bank, whereas chrysophytes were undetected. Among the three detected groups, often the prymnesiophytes were dominant in biomass. In general, all groups were present throughout the water column with abundance peaks around the DCM, but their distributions differed: pelagophytes were located deeper than the other two groups, chlorophytes presented two peaks and prymnesiophytes exhibited surface abundances comparable to those at the DCM. This study offers first indications that the vertical distribution of different PPE groups is heterogeneous within the DCM.

Chronic Polyaromatic Hydrocarbon (PAH) Contamination Is a Marginal Driver for Community Diversity and Prokaryotic Predicted Functioning in Coastal Sediments

Benthic microorganisms are key players in the recycling of organic matter and recalcitrant compounds such as polyaromatic hydrocarbons (PAHs) in coastal sediments. Despite their ecological importance, the response of microbial communities to chronic PAH pollution, one of the major threats to coastal ecosystems, has received very little attention. In one of the largest surveys performed so far on coastal sediments, the diversity and composition of microbial communities inhabiting both chronically contaminated and non-contaminated coastal sediments were investigated using high-throughput sequencing on the 18S and 16S rRNA genes. Prokaryotic alpha-diversity showed significant association with salinity, temperature, and organic carbon content. The effect of particle size distribution was strong on eukaryotic diversity. Similarly to alpha-diversity, beta-diversity patterns were strongly influenced by the environmental filter, while PAHs had no influence on the prokaryotic community structure and a weak impact on the eukaryotic community structure at the continental scale. However, at the regional scale, PAHs became the main driver shaping the structure of bacterial and eukaryotic communities. These patterns were not found for PICRUSt predicted prokaryotic functions, thus indicating some degree of functional redundancy. Eukaryotes presented a greater potential for their use as PAH contamination biomarkers, owing to their stronger response at both regional and continental scales.

Environmental Sequencing Provides Reasonable Estimates of the Relative Abundance of Specific Picoeukaryotes

High-throughput sequencing (HTS) is revolutionizing environmental surveys of microbial diversity in the three domains of life by providing detailed information on which taxa are present in microbial assemblages. However, it is still unclear how the relative abundance of specific taxa gathered by HTS correlates with cell abundances. Here, we quantified the relative cell abundance of 6 picoeukaryotic taxa in 13 planktonic samples from 6 European coastal sites using epifluorescence microscopy on tyramide signal amplification-fluorescence in situ hybridization preparations. These relative abundance values were then compared with HTS data obtained in three separate molecular surveys: 454 sequencing of the V4 region of the 18S ribosomal DNA (rDNA) using DNA and RNA extracts (DNA-V4 and cDNA-V4) and Illumina sequencing of the V9 region (cDNA-V9). The microscopic and molecular signals were generally correlated, indicating that a relative increase in specific 18S rDNA was the result of a large proportion of cells in the given taxa. Despite these positive correlations, the slopes often deviated from 1, precluding a direct translation of sequences to cells. Our data highlighted clear differences depending on the nucleic acid template or the 18S rDNA region targeted. Thus, the molecular signal obtained using cDNA templates was always closer to relative cell abundances, while the V4 and V9 regions gave better results depending on the taxa. Our data support the quantitative use of HTS data but warn about considering it as a direct proxy of cell abundances.

IMPORTANCE

Direct studies on marine picoeukaryotes by epifluorescence microscopy are problematic due to the lack of morphological features and due to the limited number and poor resolution of specific phylogenetic probes used in fluorescence in situ hybridization (FISH) routines. As a consequence, there is an increasing use of molecular methods, including high-throughput sequencing (HTS), to study marine microbial diversity. HTS can provide a detailed picture of the taxa present in a community and can reveal diversity not evident using other methods, but it is still unclear what the meaning of the sequence abundance in a given taxon is. Our aim is to investigate the correspondence between the relative HTS signal and relative cell abundances in selected picoeukaryotic taxa. Environmental sequencing provides reasonable estimates of the relative abundance of specific taxa. Better results are obtained when using RNA extracts as the templates, while the region of 18S ribosomal DNA had different influences depending on the taxa assayed.

Comparison of three clustering approaches for detecting novel environmental microbial diversity

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

Large variability of bathypelagic microbial eukaryotic communities across the world's oceans

In this work, we study the diversity of bathypelagic microbial eukaryotes (0.8-20 μm) in the global ocean. Seawater samples from 3000 to 4000 m depth from 27 stations in the Atlantic, Pacific and Indian Oceans were analyzed by pyrosequencing the V4 region of the 18S ribosomal DNA. The relative abundance of the most abundant operational taxonomic units agreed with the results of a parallel metagenomic analysis, suggesting limited PCR biases in the tag approach. Although rarefaction curves for single stations were seldom saturated, the global analysis of all sequences together suggested an adequate recovery of bathypelagic diversity. Community composition presented a large variability among samples, which was poorly explained by linear geographic distance. In fact, the similarity between communities was better explained by water mass composition (26% of the variability) and the ratio in cell abundance between prokaryotes and microbial eukaryotes (21%). Deep diversity appeared dominated by four taxonomic groups (Collodaria, Chrysophytes, Basidiomycota and MALV-II) appearing in different proportions in each sample. Novel diversity amounted to 1% of the pyrotags and was lower than expected. Our study represents an essential step in the investigation of bathypelagic microbial eukaryotes, indicating dominating taxonomic groups and suggesting idiosyncratic assemblages in distinct oceanic regions.

Comparative study of the validity of three regions of the 18S-rRNA gene for massively parallel sequencing-based monitoring of the planktonic eukaryote community

The nuclear 18S-rRNA gene has been used as a metabarcoding marker in massively parallel sequencing (MPS)-based environmental surveys for plankton biodiversity research. However, different hypervariable regions have been used in different studies, and their utility has been debated among researchers. In this study, detailed investigations into 18S-rRNA were carried out; we investigated the effective number of sequences deposited in international nucleotide sequence databases (INSDs), the amplification bias, and the amplicon sequence variability among the three variable regions, V1-3, V4-5 and V7-9, using in silico polymerase chain reaction (PCR) amplification based on INSDs. We also examined the primer universality and the taxonomic identification power, using MPS-based environmental surveys in the Sea of Okhotsk, to determine which region is more useful for MPS-based monitoring. The primer universality was not significantly different among the three regions, but the number of sequences deposited in INSDs was markedly larger for the V4-5 region than for the other two regions. The sequence variability was significantly different, with the highest variability in the V1-3 region, followed by the V7-9 region, and the lowest variability in the V4-5 region. The results of the MPS-based environmental surveys showed significantly higher identification power in the V1-3 and V7-9 regions than in the V4-5 region, but no significant difference was detected between the V1-3 and V7-9 regions. We therefore conclude that the V1-3 region will be the most suitable for future MPS-based monitoring of natural eukaryote communities, as the number of sequences deposited in INSDs increases.

Marine protist diversity in European coastal waters and sediments as revealed by high-throughput sequencing

Although protists are critical components of marine ecosystems, they are still poorly characterized. Here we analysed the taxonomic diversity of planktonic and benthic protist communities collected in six distant European coastal sites. Environmental deoxyribonucleic acid (DNA) and ribonucleic acid (RNA) from three size fractions (pico-, nano- and micro/mesoplankton), as well as from dissolved DNA and surface sediments were used as templates for tag pyrosequencing of the V4 region of the 18S ribosomal DNA. Beta-diversity analyses split the protist community structure into three main clusters: picoplankton-nanoplankton-dissolved DNA, micro/mesoplankton and sediments. Within each cluster, protist communities from the same site and time clustered together, while communities from the same site but different seasons were unrelated. Both DNA and RNA-based surveys provided similar relative abundances for most class-level taxonomic groups. Yet, particular groups were overrepresented in one of the two templates, such as marine alveolates (MALV)-I and MALV-II that were much more abundant in DNA surveys. Overall, the groups displaying the highest relative contribution were Dinophyceae, Diatomea, Ciliophora and Acantharia. Also, well represented were Mamiellophyceae, Cryptomonadales, marine alveolates and marine stramenopiles in the picoplankton, and Monadofilosa and basal Fungi in sediments. Our extensive and systematic sequencing of geographically separated sites provides the most comprehensive molecular description of coastal marine protist diversity to date.

Estimating Protistan Diversity Using High-Throughput Sequencing

Sequencing hypervariable regions from the 18S rRNA gene is commonly employed to characterize protistan biodiversity, yet there are concerns that short reads do not provide the same taxonomic resolution as full-length sequences. A total of 7,432 full-length sequences were used to perform an in silico analysis of how sequences of various lengths and target regions impact downstream ecological interpretations. Sequences that were longer than 400 nucleotides and included the V4 hypervariable region generated results similar to those derived from full-length 18S rRNA gene sequences. Present high-throughput sequencing capabilities are approaching protistan diversity estimation comparable to whole gene sequences.

Diversity and distribution of unicellular opisthokonts along the European coast analysed using high-throughput sequencing

The opisthokonts are one of the major super groups of eukaryotes. It comprises two major clades: (i) the Metazoa and their unicellular relatives and (ii) the Fungi and their unicellular relatives. There is, however, little knowledge of the role of opisthokont microbes in many natural environments, especially among non-metazoan and non-fungal opisthokonts. Here, we begin to address this gap by analysing high-throughput 18S rDNA and 18S rRNA sequencing data from different European coastal sites, sampled at different size fractions and depths. In particular, we analyse the diversity and abundance of choanoflagellates, filastereans, ichthyosporeans, nucleariids, corallochytreans and their related lineages. Our results show the great diversity of choanoflagellates in coastal waters as well as a relevant representation of the ichthyosporeans and the uncultured marine opisthokonts (MAOP). Furthermore, we describe a new lineage of marine fonticulids (MAFO) that appears to be abundant in sediments. Taken together, our work points to a greater potential ecological role for unicellular opisthokonts than previously appreciated in marine environments, both in water column and sediments, and also provides evidence of novel opisthokont phylogenetic lineages. This study highlights the importance of high-throughput sequencing approaches to unravel the diversity and distribution of both known and novel eukaryotic lineages.

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

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

Next-generation environmental diversity surveys of foraminifera: preparing the future

Foraminifera are commonly defined as marine testate protists, and their diversity is mainly assessed on the basis of the morphology of their agglutinated or mineralized tests. Diversity surveys based on environmental DNA (eDNA) have dramatically changed this view by revealing an unexpected diversity of naked and organic-walled lineages as well as detecting foraminiferal lineages in soil and freshwater environments. Moreover, single-cell analyses have allowed discrimination among genetically distinctive types within almost every described morphospecies. In view of these studies, the foraminiferal diversity appeared to be largely underestimated, but its accurate estimation was impeded by the low speed and coverage of a cloning-based eDNA approach. With the advent of high-throughput sequencing (HTS) technologies, these limitations disappeared in favor of exhaustive descriptions of foraminiferal diversity in numerous samples. Yet, the biases and errors identified in early HTS studies raised some questions about the accuracy of HTS data and their biological interpretation. Among the most controversial issues affecting the reliability of HTS diversity estimates are (1) the impact of technical and biological biases, (2) the sensitivity and specificity of taxonomic sequence assignment, (3) the ability to distinguish rare species, and (4) the quantitative interpretation of HTS data. Here, we document the lessons learned from previous HTS surveys and present the current advances and applications focusing on foraminiferal eDNA. We discuss the problems associated with HTS approaches and predict the future trends and avenues that hold promises for surveying foraminiferal diversity accurately and efficiently.

Picoeukaryotic diversity and distribution in the subtropical-tropical South China Sea

Little is known regarding the diversity and distribution of picoeukaryotes in the northwestern Pacific Ocean, although these organisms are vital components of their environment. Here, we used a culture-independent approach to assess the 18S rDNA diversity of picoeukaryotes at six sampling sites along a transcontinental section of the South China Sea. The Alveolata group comprised 58.6% of the clones and was mainly represented by the novel marine alveolates (MALV)-I (18.8%) and MALV-II (30.6%), corresponding to 66.5% of all operational taxonomic units. Sequences affiliated with seven clades of the novel marine stramenopiles (MAST) were widely distributed in different clone libraries. We report an entirely new group representing the deepest evolutionary branch of the Hacrobia; this finding suggests the existence of novel picoeukaryotes at a high taxonomic level. Many phylotypes could not be taxonomically assigned, indicating the presence of numerous previously unknown groups. Horizontally, picoeukaryotic assemblages in the coastal water characterized with the rare occurrence of MALV-I were distinct from offshore communities. Vertically, MAST-4 were mainly retrieved in surface waters; however, the Radiolaria (Rhizaria) were mainly detected in clone libraries from depths of 60 m. Our findings further emphasize the immense diversity of picoeukaryotes, especially in the subtropical-tropical northwestern Pacific Ocean.

Patterns of rare and abundant marine microbial eukaryotes

BACKGROUND

Biological communities are normally composed of a few abundant and many rare species. This pattern is particularly prominent in microbial communities, in which most constituent taxa are usually extremely rare. Although abundant and rare subcommunities may present intrinsic characteristics that could be crucial for understanding community dynamics and ecosystem functioning, microbiologists normally do not differentiate between them. Here, we investigate abundant and rare subcommunities of marine microbial eukaryotes, a crucial group of organisms that remains among the least-explored biodiversity components of the biosphere. We surveyed surface waters of six separate coastal locations in Europe, independently considering the picoplankton, nanoplankton, and microplankton/mesoplankton organismal size fractions.

RESULTS

Deep Illumina sequencing of the 18S rRNA indicated that the abundant regional community was mostly structured by organismal size fraction, whereas the rare regional community was mainly structured by geographic origin. However, some abundant and rare taxa presented similar biogeography, pointing to spatiotemporal structure in the rare microeukaryote biosphere. Abundant and rare subcommunities presented regular proportions across samples, indicating similar species-abundance distributions despite taxonomic compositional variation. Several taxa were abundant in one location and rare in other locations, suggesting large oscillations in abundance. The substantial amount of metabolically active lineages found in the rare biosphere suggests that this subcommunity constitutes a diversity reservoir that can respond rapidly to environmental change.

CONCLUSIONS

We propose that marine planktonic microeukaryote assemblages incorporate dynamic and metabolically active abundant and rare subcommunities, with contrasting structuring patterns but fairly regular proportions, across space and time.

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

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