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.

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.

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.

How much temporal recurrence is present in microbial assemblages is still an unanswered ecological question. Even though marked seasonal changes have been reported for whole microbial communities, less is known on the dynamics and seasonality of individual taxa. Here, we aim at understanding microbial recurrence at three different levels: community, taxonomic group and operational taxonomic units (OTUs). For that, we focused on a model microbial eukaryotic community populating a long‐term marine microbial observatory using 18S rRNA gene data from two organismal size fractions: the picoplankton (0.2–3 µm) and the nanoplankton (3–20 µm). We have developed an index to quantify recurrence in particular taxa. We found that community structure oscillated systematically between two main configurations corresponding to winter and summer over the 10 years studied. A few taxonomic groups such as Mamiellophyceae or MALV‐III presented clear recurrence (i.e., seasonality), whereas 13%–19% of the OTUs in both size fractions, accounting for ~40% of the relative abundance, featured recurrent dynamics. Altogether, our work links long‐term whole community dynamics with that of individual OTUs and taxonomic groups, indicating that recurrent and non‐recurrent changes characterize the dynamics of microbial assemblages.

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.

Background The ocean microbiota modulates global biogeochemical cycles and changes in its configuration may have largescale 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 the previous, 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.

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. IMPORTANCEDirect 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. Protists are key components of marine ecosystems, being major players in the global respiration and production budgets (1, 2) and playing central roles in marine food webs (3). Despite their importance and ubiquity, it was only during the past decade that environmental studies, based on molecular (i.e., culture-independent) techniques, revealed an unsuspected protist diversity in a large variety of marine ecosystems (4-13). These studies we...

Surveying microbial diversity and function is accomplished by combining complementary molecular tools. Among them, metagenomics is a PCR free approach that contains all genetic information from microbial assemblages and is today performed at a relatively large scale and reasonable cost, mostly based on very short reads. Here, we investigated the potential of metagenomics to provide taxonomic reports of marine microbial eukaryotes. We prepared a curated database with reference sequences of the V4 region of 18S rDNA clustered at 97% similarity and used this database to extract and classify metagenomic reads. More than half of them were unambiguously affiliated to a unique reference whilst the rest could be assigned to a given taxonomic group. The overall diversity reported by metagenomics was similar to that obtained by amplicon sequencing of the V4 and V9 regions of the 18S rRNA gene, although either one or both of these amplicon surveys performed poorly for groups like Excavata, Amoebozoa, Fungi and Haptophyta. We then studied the diversity of picoeukaryotes and nanoeukaryotes using 91 metagenomes from surface down to bathypelagic layers in different oceans, unveiling a clear taxonomic separation between size fractions and depth layers. Finally, we retrieved long rDNA sequences from assembled metagenomes that improved phylogenetic reconstructions of particular groups. Overall, this study shows metagenomics as an excellent resource for taxonomic exploration of marine microbial eukaryotes.