Seasonality of biogeochemically relevant microbial genes in a coastal ocean microbiome.

Microbes drive the biogeochemical cycles of marine ecosystems through their vast metabolic diversity. While we have a fairly good understanding of the spatial distribution of these metabolic processes in various ecosystems, less is known about their seasonal dynamics. We investigated the annual patterns of 21 biogeochemical relevant functions in an oligotrophic coastal ocean site by analysing the presence of key genes, analysing high-rank gene taxonomy and the dynamics of nucleotide variants. Most genes presented seasonality: photoheterotrophic processes were enriched during spring, phosphorous-related genes were dominant during summer, coinciding with potential phosphate limitation, and assimilatory nitrate reductases appeared mostly during summer and autumn, correlating negatively with nitrate availability. Additionally, we identified the main taxa driving each function at each season and described the role of underrecognized taxa such as Litoricolaceae in carbon fixation (rbcL), urea degradation (ureC), and CO oxidation (coxL). Finally, the seasonality of single variants of some families presented a decoupling between the taxonomic abundance patterns and the functional gene patterns, implying functional specialization of the different genera. Our study unveils the seasonality of key biogeochemical functions and the main taxonomic groups that harbour these relevant functions in a coastal ocean ecosystem.

Global diversity and distribution of aerobic anoxygenic phototrophs in the tropical and subtropical oceans.

The aerobic anoxygenic phototrophic (AAP) bacteria are common in most marine environments but their global diversity and biogeography remain poorly characterized. Here, we analyzed AAP communities across 113 globally-distributed surface ocean stations sampled during the Malaspina Expedition in the tropical and subtropical ocean. By means of amplicon sequencing of the pufM gene, a genetic marker for this functional group, we show that AAP communities along the surface ocean were mainly composed of members of the Halieaceae (Gammaproteobacteria), which were adapted to a large range of environmental conditions, and of different clades of the Alphaproteobacteria, which seemed to dominate under particular circumstances, such as in the oligotrophic gyres. AAP taxa were spatially structured within each of the studied oceans, with communities from adjacent stations sharing more taxonomic similarities. AAP communities were composed of a large pool of rare members and several habitat specialists. When compared to the surface ocean prokaryotic and picoeukaryotic communities, it appears that AAP communities display an idiosyncratic global biogeographical pattern, dominated by selection processes and less influenced by dispersal limitation. Our study contributes to the understanding of how AAP communities are distributed in the horizontal dimension and the mechanisms underlying their distribution across the global surface ocean.

Bacterial responses to background organic pollutants in the northeast subarctic Pacific Ocean.

Thousands of man-made synthetic chemicals are released to oceans and compose the anthropogenic dissolved organic carbon (ADOC). Little is known about the effects of this chronic pollution on marine microbiome activities. In this study, we measured the pollution level at three sites in the Northeast Subarctic Pacific Ocean (NESAP) and investigated how mixtures of three model families of ADOC at different environmentally relevant concentrations affected naturally occurring marine bacterioplankton communities' structure and metabolic functioning. The offshore northernmost site (North) had the lowest concentrations of hydrocarbons, as well as organophosphate ester plasticizers, contrasting with the two other continental shelf sites, the southern coastal site (South) being the most contaminated. At North, ADOC stimulated bacterial growth and promoted an increase in the contribution of some Gammaproteobacteria groups (e.g. Alteromonadales) to the 16 rRNA pool. These groups are described as fast responders after oil spills. In contrast, minor changes in South microbiome activities were observed. Gene expression profiles at Central showed the coexistence of ADOC degradation and stress-response strategies to cope with ADOC toxicities. These results show that marine microbial communities at three distinct domains in NESAP are influenced by background concentrations of ADOC, expanding previous assessments for polar and temperate waters.

Differential recruitment of opportunistic taxa leads to contrasting abilities in carbon processing by bathypelagic and surface microbial communities.

Different factors affect the way dissolved organic matter (DOM) is processed in the ocean water column, including environmental conditions and the functional capabilities of the communities. Recent studies have shown that bathypelagic prokaryotes are metabolically flexible, but whether this versatility translates into a higher ability to process DOM has been barely explored. Here we performed a multifactorial transplant experiment to compare the growth, activity and changes in DOM quality in surface and bathypelagic waters inoculated with either surface or bathypelagic prokaryotic communities. The effect of nutrient additions to surface waters was also explored. Despite no differences in the cell abundance of surface and deep ocean prokaryotes were observed in any of the treatments, in surface waters with nutrients the heterotrophic production of surface prokaryotes rapidly decreased. Conversely, bathypelagic communities displayed a sustained production throughout the experiment. Incubations with surface prokaryotes always led to a significant accumulation of recalcitrant compounds, which did not occur with bathypelagic prokaryotes, suggesting they have a higher ability to process DOM. These contrasting abilities could be explained by the recruitment of a comparatively larger number of opportunistic taxa within the bathypelagic assemblages, which likely resulted in a broader community capability of substrate utilization.

Diverse patterns of correspondence between protist metabarcodes and protist metagenome-assembled genomes.

Two common approaches to study the composition of environmental protist communities are metabarcoding and metagenomics. Raw metabarcoding data are usually processed into Operational Taxonomic Units (OTUs) or amplicon sequence variants (ASVs) through clustering or denoising approaches, respectively. Analogous approaches are used to assemble metagenomic reads into metagenome-assembled genomes (MAGs). Understanding the correspondence between the data produced by these two approaches can help to integrate information between the datasets and to explain how metabarcoding OTUs and MAGs are related with the underlying biological entities they are hypothesised to represent. MAGs do not contain the commonly used barcoding loci, therefore sequence homology approaches cannot be used to match OTUs and MAGs. We made an attempt to match V9 metabarcoding OTUs from the 18S rRNA gene (V9 OTUs) and MAGs from the Tara Oceans expedition based on the correspondence of their relative abundances across the same set of samples. We evaluated several metrics for detecting correspondence between features in these two datasets and developed controls to filter artefacts of data structure and processing. After selecting the best-performing metrics, ranking the V9 OTU/MAG matches by their proportionality/correlation coefficients and applying a set of selection criteria, we identified candidate matches between V9 OTUs and MAGs. In some cases, V9 OTUs and MAGs could be matched with a one-to-one correspondence, implying that they likely represent the same underlying biological entity. More generally, matches we observed could be classified into 4 scenarios: one V9 OTU matches many MAGs; many V9 OTUs match many MAGs; many V9 OTUs match one MAG; one V9 OTU matches one MAG. Notably, we found some instances in which different OTU-MAG matches from the same taxonomic group were not classified in the same scenario, with all four scenarios possible even within the same taxonomic group, illustrating that factors beyond taxonomic lineage influence the relationship between OTUs and MAGs. Overall, each scenario produces a different interpretation of V9 OTUs, MAGs and how they compare in terms of the genomic and ecological diversity they represent.

Growth rates of marine prokaryotes are extremely diverse, even among closely related taxa.

Marine prokaryotes play crucial roles in ocean biogeochemical cycles, being their contribution strongly influenced by their growth rates. Hence, elucidating the variability and phylogenetic imprint of marine prokaryotes' growth rates are crucial for better determining the role of individual taxa in biogeochemical cycles. Here, we estimated prokaryotic growth rates at high phylogenetic resolution in manipulation experiments using water from the northwestern Mediterranean Sea. Experiments were run in the four seasons with different treatments that reduced growth limiting factors: predators, nutrient availability, viruses, and light. Single-amplicon sequence variants (ASVs)-based growth rates were calculated from changes in estimated absolute abundances using total prokaryotic abundance and the proportion of each individual ASV. The trends obtained for growth rates in the different experiments were consistent with other estimates based on total cell-counts, catalyzed reporter deposition fluorescence in situ hybridization subcommunity cell-counts or metagenomic-operational taxonomic units (OTUs). Our calculations unveil a broad range of growth rates (0.3-10 d-1) with significant variability even within closely related ASVs. Likewise, the impact of growth limiting factors changed over the year for individual ASVs. High numbers of responsive ASVs were shared between winter and spring seasons, as well as throughout the year in the treatments with reduced nutrient limitation and viral pressure. The most responsive ASVs were rare in the in situ communities, comprising a large pool of taxa with the potential to rapidly respond to environmental changes. Essentially, our results highlight the lack of phylogenetic coherence in the range of growth rates observed, and differential responses to the various limiting factors, even for closely related taxa.

Phenology and ecological role of aerobic anoxygenic phototrophs in freshwaters.

BACKGROUND: Aerobic anoxygenic phototrophic (AAP) bacteria are heterotrophic bacteria that supply their metabolism with light energy harvested by bacteriochlorophyll-a-containing reaction centers. Despite their substantial contribution to bacterial biomass, microbial food webs, and carbon cycle, their phenology in freshwater lakes remains unknown. Hence, we investigated seasonal variations of AAP abundance and community composition biweekly across 3 years in a temperate, meso-oligotrophic freshwater lake. RESULTS: AAP bacteria displayed a clear seasonal trend with a spring maximum following the bloom of phytoplankton and a secondary maximum in autumn. As the AAP bacteria represent a highly diverse assemblage of species, we followed their seasonal succession using the amplicon sequencing of the pufM marker gene. To enhance the accuracy of the taxonomic assignment, we developed new pufM primers that generate longer amplicons and compiled the currently largest database of pufM genes, comprising 3633 reference sequences spanning all phyla known to contain AAP species. With this novel resource, we demonstrated that the majority of the species appeared during specific phases of the seasonal cycle, with less than 2% of AAP species detected during the whole year. AAP community presented an indigenous freshwater nature characterized by high resilience and heterogenic adaptations to varying conditions of the freshwater environment. CONCLUSIONS: Our findings highlight the substantial contribution of AAP bacteria to the carbon flow and ecological dynamics of lakes and unveil a recurrent and dynamic seasonal succession of the AAP community. By integrating this information with the indicator of primary production (Chlorophyll-a) and existing ecological models, we show that AAP bacteria play a pivotal role in the recycling of dissolved organic matter released during spring phytoplankton bloom. We suggest a potential role of AAP bacteria within the context of the PEG model and their consideration in further ecological models.

Seasonal niche differentiation among closely related marine bacteria.

Bacteria display dynamic abundance fluctuations over time in marine environments, where they play key biogeochemical roles. Here, we characterized the seasonal dynamics of marine bacteria in a coastal oligotrophic time series station, tested how similar the temporal niche of closely related taxa is, and what are the environmental parameters modulating their seasonal abundance patterns. We further explored how conserved the niche is at higher taxonomic levels. The community presented recurrent patterns of seasonality for 297 out of 6825 amplicon sequence variants (ASVs), which constituted almost half of the total relative abundance (47%). For certain genera, niche similarity decreased as nucleotide divergence in the 16S rRNA gene increased, a pattern compatible with the selection of similar taxa through environmental filtering. Additionally, we observed evidence of seasonal differentiation within various genera as seen by the distinct seasonal patterns of closely related taxa. At broader taxonomic levels, coherent seasonal trends did not exist at the class level, while the order and family ranks depended on the patterns that existed at the genus level. This study identifies the coexistence of closely related taxa for some bacterial groups and seasonal differentiation for others in a coastal marine environment subjected to a strong seasonality.

Performance of the melting seawater-ice elution method on the metabarcoding characterization of benthic protist communities.

Massive amplicon sequencing approaches to characterize the diversity of microbial eukaryotes in sediments are scarce and controls about the effects introduced by different methods to recover DNA are lacking. In this study, we compare the performance of the melting seawater-ice elution method on the characterization of benthic protist communities by 18S rRNA gene metabarcoding with results obtained by direct cell lysis and DNA purification from sediments. Even though the most abundant operational taxonomic units were recovered by both methods, eluted samples yielded higher richness than samples undergoing direct lysis. Both treatments allowed recovering the same taxonomic groups, although we observed significant differences in terms of relative abundance for some of them. Dinoflagellata and Ciliophora strongly dominated the community in eluted samples (> 80% reads). In directly lysed samples, they only represented 37%, while groups like Fungi and Ochrophytes were highly represented (> 20% reads respectively). Our results show that the elution process yields a higher protist richness estimation, most likely as a result of the higher sample volume used to recover organisms as compared to commonly used volumes for direct benthic DNA purification. Motile groups, like dinoflagellates and ciliates, are logically more enriched during the elution process.

Seasonal impact of grazing, viral mortality, resource availability and light on the group-specific growth rates of coastal Mediterranean bacterioplankton.

Estimation of prokaryotic growth rates is critical to understand the ecological role and contribution of different microbes to marine biogeochemical cycles. However, there is a general lack of knowledge on what factors control the growth rates of different prokaryotic groups and how these vary between sites and along seasons at a given site. We carried out several manipulation experiments during the four astronomical seasons in the coastal NW Mediterranean in order to evaluate the impact of grazing, viral mortality, resource competition and light on the growth and loss rates of prokaryotes. Gross and net growth rates of different bacterioplankton groups targeted by group-specific CARD-FISH probes and infrared microscopy (for aerobic anoxygenic phototrophs, AAP), were calculated from changes in cell abundances. Maximal group-specific growth rates were achieved when both predation pressure and nutrient limitation were experimentally minimized, while only a minimal effect of viral pressure on growth rates was observed; nevertheless, the response to predation removal was more remarkable in winter, when the bacterial community was not subjected to nutrient limitation. Although all groups showed increases in their growth rates when resource competition as well as grazers and viral pressure were reduced, Alteromonadaceae consistently presented the highest rates in all seasons. The response to light availability was generally weaker than that to the other factors, but it was variable between seasons. In summer and spring, the growth rates of AAP were stimulated by light whereas the growth of the SAR11 clade (likely containing proteorhodopsin) was enhanced by light in all seasons. Overall, our results set thresholds on bacterioplankton group-specific growth and mortality rates and contribute to estimate the seasonally changing contribution of various bacterioplankton groups to the function of microbial communities. Our results also indicate that the least abundant groups display the highest growth rates, contributing to the recycling of organic matter to a much greater extent than what their abundances alone would predict.

Long-term seasonal and interannual variability of marine aerobic anoxygenic photoheterotrophic bacteria.

We studied the long-term temporal dynamics of the aerobic anoxygenic phototrophic (AAP) bacteria, a relevant functional group in the coastal marine microbial food web, using high-throughput sequencing of the pufM gene coupled with multivariate, time series and co-occurrence analyses at the Blanes Bay Microbial Observatory (NW Mediterranean). Additionally, using metagenomics, we tested whether the used primers captured accurately the seasonality of the most relevant AAP groups. Phylogroup K (Gammaproteobacteria) was the greatest contributor to community structure over all seasons, with phylogroups E and G (Alphaproteobacteria) being prevalent in spring. Diversity indices showed a clear seasonal trend, with maximum values in winter, which was inverse to that of AAP abundance. Multivariate analyses revealed sample clustering by season, with a relevant proportion of the variance explained by day length, temperature, salinity, phototrophic nanoflagellate abundance, chlorophyll a, and silicate concentration. Time series analysis showed robust rhythmic patterns of co-occurrence, but distinct seasonal behaviors within the same phylogroup, and even within different amplicon sequence variants (ASVs) conforming the same operational taxonomic unit (OTU). Altogether, our results picture the AAP assemblage as highly seasonal and recurrent but containing ecotypes showing distinctive temporal niche partitioning, rather than being a cohesive functional group.