Microbes are dominant drivers of biogeochemical processes, yet drawing a global picture of functional diversity, microbial community structure, and their ecological determinants remains a grand challenge. We analyzed 7.2 terabases of metagenomic data from 243 Tara Oceans samples from 68 locations in epipelagic and mesopelagic waters across the globe to generate an ocean microbial reference gene catalog with >40 million nonredundant, mostly novel sequences from viruses, prokaryotes, and picoeukaryotes. Using 139 prokaryote-enriched samples, containing >35,000 species, we show vertical stratification with epipelagic community composition mostly driven by temperature rather than other environmental factors or geography. We identify ocean microbial core functionality and reveal that >73% of its abundance is shared with the human gut microbiome despite the physicochemical differences between these two ecosystems.

The Mediterranean Sea is a marine biodiversity hot spot. Here we combined an extensive literature analysis with expert opinions to update publicly available estimates of major taxa in this marine ecosystem and to revise and update several species lists. We also assessed overall spatial and temporal patterns of species diversity and identified major changes and threats. Our results listed approximately 17,000 marine species occurring in the Mediterranean Sea. However, our estimates of marine diversity are still incomplete as yet—undescribed species will be added in the future. Diversity for microbes is substantially underestimated, and the deep-sea areas and portions of the southern and eastern region are still poorly known. In addition, the invasion of alien species is a crucial factor that will continue to change the biodiversity of the Mediterranean, mainly in its eastern basin that can spread rapidly northwards and westwards due to the warming of the Mediterranean Sea. Spatial patterns showed a general decrease in biodiversity from northwestern to southeastern regions following a gradient of production, with some exceptions and caution due to gaps in our knowledge of the biota along the southern and eastern rims. Biodiversity was also generally higher in coastal areas and continental shelves, and decreases with depth. Temporal trends indicated that overexploitation and habitat loss have been the main human drivers of historical changes in biodiversity. At present, habitat loss and degradation, followed by fishing impacts, pollution, climate change, eutrophication, and the establishment of alien species are the most important threats and affect the greatest number of taxonomic groups. All these impacts are expected to grow in importance in the future, especially climate change and habitat degradation. The spatial identification of hot spots highlighted the ecological importance of most of the western Mediterranean shelves (and in particular, the Strait of Gibraltar and the adjacent Alboran Sea), western African coast, the Adriatic, and the Aegean Sea, which show high concentrations of endangered, threatened, or vulnerable species. The Levantine Basin, severely impacted by the invasion of species, is endangered as well.This abstract has been translated to other languages (File S1).

Viruses influence ecosystems by modulating microbial population size, diversity, metabolic outputs, and gene flow. Here, we use quantitative double-stranded DNA (dsDNA) viral-fraction metagenomes (viromes) and whole viral community morphological data sets from 43 Tara Oceans expedition samples to assess viral community patterns and structure in the upper ocean. Protein cluster cataloging defined pelagic upper-ocean viral community pan and core gene sets and suggested that this sequence space is well-sampled. Analyses of viral protein clusters, populations, and morphology revealed biogeographic patterns whereby viral communities were passively transported on oceanic currents and locally structured by environmental conditions that affect host community structure. Together, these investigations establish a global ocean dsDNA viromic data set with analyses supporting the seed-bank hypothesis to explain how oceanic viral communities maintain high local diversity.

No abstract

We analyzed heterotrophic, pelagic bacterial production and specific growth rate data from 57 studies conducted in fresh, marine and estuarine/coastal waters. Strong positive relationships were identified between 1) bacterial production and bacterial abundance and 2) bacterial production and algal biomass. The relationship between bacterial production and bacterial abundance was improved by also considering water temperature. The analysis of covariance model revealed consistent differences between fresh, marine and estuarine/coastal waters, with production consistently high in estuarine/coastal environments. The log-linear regression coefficient of abundance was not significantly different from 1.00, and this linear relationship permitted the use of specific growth rate (SGR in day(-1)) as a dependent variable. A strong relationship was identified between specific growth rate and temperature. This relationship differed slightly across the three habitats. A substantial portion of the residual variation from this relationship was accounted for by algal biomass, including the difference between marine and estuarine/coastal habitats. A small but significant difference between the fresh- and saltwater habitats remained. No significant difference between the chlorophyll effect in different habitats was identified. The model of SGR against temperature and chlorophyll was much weaker for freshwater than for marine environments. For a small subset of the data set, mean cell volume accounted for some of the residual variation in SGR. Pronounced seasonality, fluctuations in nutrient quality, and variation of the grazing environment may contribute to the unexplained variation in specific growth.

Microbial communities inhabiting a multipond solar saltern were analysed and compared using SSU rRNA polymerase chain reaction (PCR)-based fingerprintings carried out in parallel by four laboratories. A salinity gradient from seawater (3.7%) to NaCl precipitation (37%) was studied for Bacteria, Archaea and Eukarya, and laboratories applied their own techniques and protocols on the same set of samples. Members of all three domains were retrieved from all salt concentrations. Three fingerprinting techniques were used: denaturing gradient gel electrophoresis (DGGE), ribosomal internal spacer analysis (RISA), and terminal-restriction fragments length polymorphism (T-RFLP). In addition, each laboratory used its own biomass collection method and DNA extraction protocols. Prokaryotes were addressed using DGGE and RISA with different 'domain-specific' primers sets. Eukaryotes were analysed by one laboratory using DGGE and T-RFLP, but targeting the same 18S rDNA site. Fingerprints were compared through cluster analysis and non-metric multidimensional scaling plots. This exercise allowed fast comparison of microbial assemblages and determined to what extent the picture provided by each laboratory was similar to those of others. Formation of two main, salinity-based groups of samples in prokaryotes (4-15% and 22-37% salinity) was consistent for all the laboratories. When other clusters appeared, this was a result of the particular technique and the protocol used in each case, but more affected by the primers set used. Eukaryotic microorganisms changed more from pond to pond; 4-5% and 8-37% salinity were but the two main groups detected. Archaea showed the lowest number of bands whereas Eukarya showed the highest number of operational taxonomic units (OTUs) in the initial ponds. Artefacts appeared in the DGGE from ponds with extremely low microbial richness. On the other hand, different 16S rDNA fragments with the same restriction or internal transcribed spacer (ITS) length were the main limitations for T-RFLP and RISA analyses, respectively, in ponds with the highest OTUs richness. However, although the particular taxonomic composition could vary among protocols, the general structure of the microbial assemblages was maintained.

SummaryThe ocean is home to myriad small planktonic organisms that underpin the functioning of marine ecosystems. However, their spatial patterns of diversity and the underlying drivers remain poorly known, precluding projections of their responses to global changes. Here we investigate the latitudinal gradients and global predictors of plankton diversity across archaea, bacteria, eukaryotes, and major virus clades using both molecular and imaging data from Tara Oceans. We show a decline of diversity for most planktonic groups toward the poles, mainly driven by decreasing ocean temperatures. Projections into the future suggest that severe warming of the surface ocean by the end of the 21st century could lead to tropicalization of the diversity of most planktonic groups in temperate and polar regions. These changes may have multiple consequences for marine ecosystem functioning and services and are expected to be particularly significant in key areas for carbon sequestration, fisheries, and marine conservation.Video Abstract

Microbes drive ecosystems under constraints imposed by viruses. However, a lack of virus genome information hinders our ability to answer fundamental, biological questions concerning microbial communities. Here we apply single-virus genomics (SVGs) to assess whether portions of marine viral communities are missed by current techniques. The majority of the here-identified 44 viral single-amplified genomes (vSAGs) are more abundant in global ocean virome data sets than published metagenome-assembled viral genomes or isolates. This indicates that vSAGs likely best represent the dsDNA viral populations dominating the oceans. Species-specific recruitment patterns and virome simulation data suggest that vSAGs are highly microdiverse and that microdiversity hinders the metagenomic assembly, which could explain why their genomes have not been identified before. Altogether, SVGs enable the discovery of some of the likely most abundant and ecologically relevant marine viral species, such as vSAG 37-F6, which were overlooked by other methodologies.