Guillem Salazar

Structure and function of the global ocean microbiome

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.

Determinants of community structure in the global plankton interactome

Species interaction networks are shaped by abiotic and biotic factors. Here, as part of the Tara Oceans project, we studied the photic zone interactome using environmental factors and organismal abundance profiles and found that environmental factors are incomplete predictors of community structure. We found associations across plankton functional types and phylogenetic groups to be nonrandomly distributed on the network and driven by both local and global patterns. We identified interactions among grazers, primary producers, viruses, and (mainly parasitic) symbionts and validated network-generated hypotheses using microscopy to confirm symbiotic relationships. We have thus provided a resource to support further research on ocean food webs and integrating biological components into ocean models.

Metagenomic 16S rDNA Illumina tags are a powerful alternative to amplicon sequencing to explore diversity and structure of microbial communities

Sequencing of 16S rDNA polymerase chain reaction (PCR) amplicons is the most common approach for investigating environmental prokaryotic diversity, despite the known biases introduced during PCR. Here we show that 16S rDNA fragments derived from Illumina-sequenced environmental metagenomes (mi tags) are a powerful alternative to 16S rDNA amplicons for investigating the taxonomic diversity and structure of prokaryotic communities. As part of the Tara Oceans global expedition, marine plankton was sampled in three locations, resulting in 29 subsamples for which metagenomes were produced by shotgun Illumina sequencing (ca. 700 Gb). For comparative analyses, a subset of samples was also selected for Roche-454 sequencing using both shotgun (m454 tags; 13 metagenomes, ca. 2.4 Gb) and 16S rDNA amplicon (454 tags; ca. 0.075 Gb) approaches. Our results indicate that by overcoming PCR biases related to amplification and primer mismatch, mi tags may provide more realistic estimates of community richness and evenness than amplicon 454 tags. In addition, mi tags can capture expected beta diversity patterns. Using mi tags is now economically feasible given the dramatic reduction in high-throughput sequencing costs, having the advantage of retrieving simultaneously both taxonomic (Bacteria, Archaea and Eukarya) and functional information from the same microbial community.

Ubiquitous Healthy Diatoms in the Deep Sea Confirm Deep Carbon Injection by the Biological Pump

The role of the ocean as a sink for CO2 is partially dependent on the downward transport of phytoplankton cells packaged within fast-sinking particles. However, whether such fast-sinking mechanisms deliver fresh organic carbon down to the deep bathypelagic sea and whether this mechanism is prevalent across the ocean requires confirmation. Here we report the ubiquitous presence of healthy photosynthetic cells, dominated by diatoms, down to 4,000 m in the deep dark ocean. Decay experiments with surface phytoplankton suggested that the large proportion (18%) of healthy photosynthetic cells observed, on average, in the dark ocean, requires transport times from a few days to a few weeks, corresponding to sinking rates (124–732 m d−1) comparable to those of fast-sinking aggregates and faecal pellets. These results confirm the expectation that fast-sinking mechanisms inject fresh organic carbon into the deep sea and that this is a prevalent process operating across the global oligotrophic ocean.

Global diversity and biogeography of deep-sea pelagic prokaryotes

The deep-sea is the largest biome of the biosphere, and contains more than half of the whole ocean’s microbes. Uncovering their general patterns of diversity and community structure at a global scale remains a great challenge, as only fragmentary information of deep-sea microbial diversity exists based on regional-scale studies. Here we report the first globally comprehensive survey of the prokaryotic communities inhabiting the bathypelagic ocean using high-throughput sequencing of the 16S rRNA gene. This work identifies the dominant prokaryotes in the pelagic deep ocean and reveals that 50% of the operational taxonomic units (OTUs) belong to previously unknown prokaryotic taxa, most of which are rare and appear in just a few samples. We show that whereas the local richness of communities is comparable to that observed in previous regional studies, the global pool of prokaryotic taxa detected is modest (~3600 OTUs), as a high proportion of OTUs are shared among samples. The water masses appear to act as clear drivers of the geographical distribution of both particle-attached and free-living prokaryotes. In addition, we show that the deep-oceanic basins in which the bathypelagic realm is divided contain different particle-attached (but not free-living) microbial communities. The combination of the aging of the water masses and a lack of complete dispersal are identified as the main drivers for this biogeographical pattern. All together, we identify the potential of the deep ocean as a reservoir of still unknown biological diversity with a higher degree of spatial complexity than hitherto considered.

Marked seasonality of aerobic anoxygenic phototrophic bacteria in the coastal NW Mediterranean Sea as revealed by cell abundance, pigment concentration and pyrosequencing of pufM gene.

The abundance and diversity of aerobic anoxygenic phototrophs (AAPs) were studied for a year cycle at the Blanes Bay Microbial Observatory (NW Mediterranean) and their potential links to an array of environmental variables were explored. Cell numbers were low in winter and peaked in summer, showing a marked seasonality that positively correlated with day length and light at the surface. Bacteriochlorophyll a concentration, their light-harvesting pigment, was only detected between April and October, and pigment cell quota showed large variations during this period. Pyrosequencing analysis of the pufM gene revealed that the most abundant operational taxonomic units (OTUs) were affiliated to phylogroup K (Gammaproteobacteria) and uncultured phylogroup C, although they were outnumbered by alphaproteobacterial OTUs in spring. Overall, richness was higher in winter than in summer, showing an opposite trend to abundance and day length. Clustering of samples by multivariate analyses showed a clear seasonality that suggests a succession of different AAP subpopulations over time. Temperature, chlorophyll a and day length were the environmental drivers that best explained the distribution of AAP assemblages. These results indicate that AAP bacteria are highly dynamic and undergo seasonal variations in diversity and abundance mostly dictated by environmental conditions as exemplified by light availability.

Particle‐association lifestyle is a phylogenetically conserved trait in bathypelagic prokaryotes

The free‐living (FL) and particle‐attached (PA) marine microbial communities have repeatedly been proved to differ in their diversity and composition in the photic ocean and also recently in the bathypelagic ocean at a global scale. However, although high taxonomic ranks exhibit preferences for a PA or FL mode of life, it remains poorly understood whether two clear lifestyles do exist and how these are distributed across the prokaryotic phylogeny. We studied the FL (<0.8 μm) and PA (0.8–20 μm) prokaryotes at 30 stations distributed worldwide within the bathypelagic oceanic realm (2150–4000 m depth) using high‐throughput sequencing of the small subunit ribosomal RNA gene (16S rRNA). A high proportion of the bathypelagic prokaryotes were mostly found either attached to particles or freely in the surrounding water but rarely in both types of environments. In particular, this trait was deeply conserved through their phylogeny, suggesting that the deep‐ocean particles and the surrounding water constitute two highly distinct niches and that transitions from one to the other have been rare at an evolutionary timescale. As a consequence, PA and FL communities had clear alpha‐ and beta‐diversity differences that exceeded the global‐scale geographical variation. Our study organizes the bathypelagic prokaryotic diversity into a reasonable number of ecologically coherent taxa regarding their association with particles, a first step for understanding which are the microbes responsible for the processing of the dissolved and particulate pools of organic matter that have a very different biogeochemical role in the deep ocean.

Biogeography of the uncultured marine picoeukaryote MAST-4: temperature-driven distribution patterns

The MAST-4 (marine stramenopile group 4) is a widespread uncultured picoeukaryote that makes up an important fraction of marine heterotrophic flagellates. This group has low genetic divergence and is composed of a small number of putative species. We combined ARISA (automated ribosomal intergenic spacer analysis) and ITS (Internal Transcribed Spacer) clone libraries to study the biogeography of this marine protist, examining both spatial and temporal trends in MAST-4 assemblages and associated environmental factors. The most represented MAST-4 clades appeared adapted to different temperature ranges, and their distributions did not suggest clear geographical barriers for dispersal. Distant samples sharing the same temperature had very similar assemblages, especially in cold temperatures, where only one clade, E1, dominated. The most highly represented clades, A and E1, showed very little differentiation between populations from distant geographical regions. Within a single site, temporal variation also followed patterns governed by temperature. Our results contribute to the general discussion on microbial biogeography by showing strong environmental selection for some picoeukaryotes in the marine environment.

Cyanobacterial symbionts diverged in the late Cretaceous towards lineage-specific nitrogen fixation factories in single-celled phytoplankton

The unicellular cyanobacterium UCYN-A, one of the major contributors to nitrogen fixation in the open ocean, lives in symbiosis with single-celled phytoplankton. UCYN-A includes several closely related lineages whose partner fidelity, genome-wide expression and time of evolutionary divergence remain to be resolved. Here we detect and distinguish UCYN-A1 and UCYN-A2 lineages in symbiosis with two distinct prymnesiophyte partners in the South Atlantic Ocean. Both symbiotic systems are lineage specific and differ in the number of UCYN-A cells involved. Our analyses infer a streamlined genome expression towards nitrogen fixation in both UCYN-A lineages. Comparative genomics reveal a strong purifying selection in UCYN-A1 and UCYN-A2 with a diversification process ∼91 Myr ago, in the late Cretaceous, after the low-nutrient regime period occurred during the Jurassic. These findings suggest that UCYN-A diversified in a co-evolutionary process, wherein their prymnesiophyte partners acted as a barrier driving an allopatric speciation of extant UCYN-A lineages.

Weak Coherence in Abundance Patterns Between Bacterial Classes and Their Constituent OTUs Along a Regulated River

Deductions about the ecology of high taxonomic bacterial ranks (i.e., phylum, class, order) are often based on their abundance patterns, yet few studies have quantified how accurately variations in abundance of these bacterial groups represent the dynamics of the taxa within them. Using 454-pyrosequencing of the 16S rRNA gene, we investigated whether the changes in abundance of six dominant bacterial classes (Actinobacteria, Beta-/Alpha-/Gamma-proteobacteria, Flavobacteria, and Sphingobacteria) along a large dam-regulated river are reflected by those of their constituent Operational Taxonomic Units (OTUs; 97% similarity level). The environmental impact generated by the reservoirs promoted clear compositional shifts in all bacterial classes that resulted from changes in the abundance of individual OTUs rather than from the appearance of new taxa along the river. Abundance patterns at the class level represented the dynamics of only a small but variable proportion of their constituting OTUs, which were not necessarily the most abundant ones. Within most classes, we detected sub-groups of OTUs showing contrasting responses to reservoir-induced environmental changes. Overall, we show that the patterns observed at the class level fail to capture the dynamics of a significant fraction of their constituent members, calling for caution when the ecological attributes of high-ranks are to be interpreted.