Caterina R. Giner

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.

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

ABSTRACT 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.

Evaluation of Alternative High-Throughput Sequencing Methodologies for the Monitoring of Marine Picoplanktonic Biodiversity Based on rRNA Gene Amplicons

Sequencing of rRNA gene polymerase chain reaction amplicons (rRNA tags) is the most common approach for investigating microbial diversity. The recent development of high-throughput sequencing (HTS) technologies has enabled the exploration of microbial biodiversity at an unprecedented scale, greatly expanding our knowledge on the microbiomes of marine ecosystems. These approaches provide accurate, fast and cost efficient observations of the marine communities, and thus, may be suitable tools in biodiversity monitoring programs. To reach this goal, consistent and comparable methodologies must be used over time and within sites. Here, we have performed a cross-platform study of the two most common HTS methodologies, i.e. 454-pyrosequencing and Illumina tags to evaluate their usefulness in biodiversity monitoring and assessment of environmental status. Plankton biodiversity has been compared through both methodologies by sequencing the 16S and 18S rRNA genes of a set of samples collected in the coast of Barcelona (NW Mediterranean). The results show that, despite differences observed in the rare OTUs retrieved, both platforms provide a comparable view of the marine picoplankton communities. On a taxonomic level, there was an accurate overlap in the detected phyla between the two methods and the overall estimates of alpha- and betadiversity were comparable. In addition, we explored the concept of ‘indicator species’ and found that certain taxa (i.e. members of the Gammaproteobacteria among others) as well as the ratio between some phylogenetic groups (i.e. the ratio of Alphaproteobacteria / Gammaproteobacteria, Alteromonas / SAR11 and Alteromonas + Oceanospirillales / SAR11) have potential for being useful indicators of environmental status. The data show that implementing new protocols and identifying indicators of environmental status based on rRNA amplicon sequencing is feasible, and that is worth exploring whether the identified indices are universally applicable.

Different processes shape prokaryotic and picoeukaryotic assemblages in the sunlit ocean microbiome

The smallest members of the sunlit-ocean microbiome (prokaryotes and picoeukaryotes) participate in a plethora of ecosystem functions with planetary-scale effects. Understanding the processes determining the spatial turnover of this assemblage can help us better comprehend the links between microbiome species composition and ecosystem function. Ecological theory predicts that selection, dispersal and drift are main drivers of species distributions, yet, the relative quantitative importance of these ecological processes in structuring the surface-ocean microbiome is barely known. Here we quantified the role of selection, dispersal and drift in structuring surface-ocean prokaryotic and picoeukaryotic assemblages by using community DNA-sequence data collected during the global Malaspina expedition. We found that dispersal limitation was the dominant process structuring picoeukaryotic communities, while a balanced combination of dispersal limitation, selection and drift shaped prokaryotic counterparts. Subsequently, we determined the agents exerting abiotic selection as well as the spatial patterns emerging from the action of different ecological processes. We found that selection exerted via temperature had a strong influence on the structure of prokaryotic communities, particularly on species co-occurrences, a pattern not observed among communities of picoeukaryotes. Other measured abiotic variables had limited selective effects on microbiome structure. Picoeukaryotes presented a higher differentiation between neighbouring communities and a higher distance-decay when compared to prokaryotes, agreeing with their higher dispersal limitation. Finally, drift seemed to have a limited role in structuring the sunlit-ocean microbiome. The different predominance of ecological processes acting on particular subsets of the ocean microbiome suggests uneven responses to environmental change. SIGNIFICANCE STATEMENT The global ocean contains one of the largest microbiomes on Earth and changes on its structure can impact the functioning of the biosphere. Yet, we are far from understanding the mechanisms that structure the global ocean microbiome, that is, the relative importance of environmental selection, dispersal and random events (drift). We evaluated the role of these processes at the global scale, based on data derived from a circumglobal expedition and found that these ecological processes act differently on prokaryotes and picoeukaryotes, two of the main components of the ocean microbiome. Our work represents a significant contribution to understand the assembly of marine microbial communities, providing also insights on the links between ecological mechanisms, microbiome structure and ecosystem function.

Quantifying long-term predictability in microbial plankton dynamics

Determining predictability in community turnover is a key ecological question. In the microbial world, seasonality has been reported for communities inhabiting temperate zones, but not much is known on seasonality for individual species. Specifically, we have a vague understanding on the amount of species displaying predictability during temporal community turnover as well as on their dynamics. Here we developed a ‘Recurrence Index’ to quantify predictability in microbial species. Applying our index to 18S rDNA metabarcoding data from one of the longest temporal observatories of marine plankton we determined that 13% of the picoeukaryotic and 19% of the nanoeukaryotic species, accounting for about 40% of the community abundance in both fractions, feature predictable dynamics when sampled monthly during 10 years. Thus, most of the species analysed had unpredictable temporal abundance patterns. Altogether, we show that species with both predictable and unpredictable temporal dynamics can occur within the same seasonal microbial community.