Genomic and phenotypic characterization of 26 novel marine bacterial strains with relevant biogeochemical roles and widespread presence across the global ocean.

Prokaryotes dominate global oceans and shape biogeochemical cycles, yet most taxa remain uncultured and uncharacterized as of today. Here we present the characterization of 26 novel marine bacterial strains from a large isolate collection obtained from Blanes Bay (NW Mediterranean) microcosm experiments made in the four seasons. Morphological, cultural, biochemical, physiological, nutritional, genomic, and phylogenomic analyses were used to characterize and phylogenetically place the novel isolates. The strains represent 23 novel bacterial species and six novel genera: three novel species pertaining to class Alphaproteobacteria (families Rhodobacteraceae and Sphingomonadaceae), six novel species and three new genera from class Gammaproteobacteria (families Algiphilaceae, Salinispheraceae, and Alteromonadaceae), 13 novel species and three novel genera from class Bacteroidia (family Flavobacteriaceae), and one new species from class Rhodothermia (family Rubricoccaceae). The bacteria described here have potentially relevant roles in the cycles of carbon (e.g., carbon fixation or energy production via proteorhodopsin), nitrogen (e.g., denitrification or use of urea), sulfur (oxidation of sulfur compounds), phosphorus (acquisition and use of different forms of phosphorus and remodeling of membrane phospholipids), and hydrogen (oxidation of hydrogen to obtain energy). We mapped the genomes of the presented strains to the Tara Oceans metagenomes to reveal that these strains were globally distributed, with those of the family Flavobacteriaceae being the most widespread and abundant, while Rhodothermia being the rarest and most localized. While molecular-only approaches are also important, our study stresses the importance of culturing as a powerful tool to further understand the functioning of marine bacterial communities.

Patterns of prokaryotic activity along the marine planktonic matter continuum.

Prokaryotic abundance and activity are commonly assessed by dividing them into two size-fractions: free-living and attached to particles. Nevertheless, as organic matter, essential for the growth of heterotrophic prokaryotes, is present in the environment in a continuum of sizes, from purely dissolved to large particles, defining the activity of the prokaryotic community would be more accurate by considering all the distinct size fractions. To achieve this, we measured prokaryotic abundance (PA), heterotrophic prokaryotic activity (as leucine incorporation) and extracellular enzyme activities at a coastal site in the NW Mediterranean Sea. We conducted measurements in both bulk seawater and size fractionated samples sequentially passing through 5 different filter types: 0.2-0.8-3-5-10 µm pore size. Our results indicate that the fraction <0.8 µm contained the highest percentage of cells (91.62 ±â€¯1.08 %) and leucine incorporation rates (72.21 ±â€¯3.51 %). Most of the extracellular enzyme activity appeared in the dissolved fraction (<0.2 µm; 19.82-79.41 %), yet the specific activity of the enzymes (per cell activity) was 100-1000 times higher in the particulate (>0.8 µm) than in the free-living (0.2-0.8 µm) fraction. The size fraction with highest specific activities for leucine incorporation and most of the enzyme activities (ß-glucosidase, esterase, Leu-aminopeptidase and alkaline phosphatase) was the 5-10 µm fraction. In contrast, the higher specific chitobiase activity in the >10 µm fraction, suggests that the prokaryotic community colonizing large particles might be more specialized in the hydrolysis of organic matter of zooplanktonic origin than the community colonizing smaller particles.

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.

High amino acid osmotrophic incorporation by marine eukaryotic phytoplankton revealed by click chemistry.

The osmotrophic uptake of dissolved organic compounds in the ocean is considered to be dominated by heterotrophic prokaryotes, whereas the role of planktonic eukaryotes is still unclear. We explored the capacity of natural eukaryotic plankton communities to incorporate the synthetic amino acid L-homopropargylglycine (HPG, analogue of methionine) using biorthogonal noncanonical amino acid tagging (BONCAT), and we compared it with prokaryotic HPG use throughout a 9-day survey in the NW Mediterranean. BONCAT allows to fluorescently identify translationally active cells, but it has never been applied to natural eukaryotic communities. We found a large diversity of photosynthetic and heterotrophic eukaryotes incorporating HPG into proteins, with dinoflagellates and diatoms showing the highest percentages of BONCAT-labelled cells (49 ± 25% and 52 ± 15%, respectively). Among them, pennate diatoms exhibited higher HPG incorporation in the afternoon than in the morning, whereas small (≤5 µm) photosynthetic eukaryotes and heterotrophic nanoeukaryotes showed the opposite pattern. Centric diatoms (e.g. Chaetoceros, Thalassiosira, and Lauderia spp.) dominated the eukaryotic HPG incorporation due to their high abundances and large sizes, accounting for up to 86% of the eukaryotic BONCAT signal and strongly correlating with bulk 3H-leucine uptake rates. When including prokaryotes, eukaryotes were estimated to account for 19-31% of the bulk BONCAT signal. Our results evidence a large complexity in the osmotrophic uptake of HPG, which varies over time within and across eukaryotic groups and highlights the potential of BONCAT to quantify osmotrophy and protein synthesis in complex eukaryotic communities.

The active free-living bathypelagic microbiome is largely dominated by rare surface taxa.

A persistent microbial seed bank is postulated to sustain the marine biosphere, and recent findings show that prokaryotic taxa present in the ocean's surface dominate prokaryotic communities throughout the water column. Yet, environmental conditions exert a tight control on the activity of prokaryotes, and drastic changes in these conditions are known to occur from the surface to deep waters. The simultaneous characterization of the total (DNA) and active (i.e. with potential for protein synthesis, RNA) free-living communities in 13 stations distributed across the tropical and subtropical global ocean allowed us to assess their change in structure and diversity along the water column. We observed that active communities were surprisingly more similar along the vertical gradient than total communities. Looking at the vertical connectivity of the active vs. the total communities, we found that taxa detected in the surface sometimes accounted for more than 75% of the active microbiome of bathypelagic waters (50% on average). These active taxa were generally rare in the surface, representing a small fraction of all the surface taxa. Our findings show that the drastic vertical change in environmental conditions leads to the inactivation and disappearance of a large proportion of surface taxa, but some surface-rare taxa remain active (or with potential for protein synthesis) and dominate the bathypelagic active microbiome.

Growth and mortality of aerobic anoxygenic phototrophs in the North Pacific Subtropical Gyre.

Aerobic anoxygenic phototrophic (AAP) bacteria harvest light energy using bacteriochlorophyll-containing reaction centers to supplement their mostly heterotrophic metabolism. While their abundance and growth have been intensively studied in coastal environments, much less is known about their activity in oligotrophic open ocean regions. Therefore, we combined in situ sampling in the North Pacific Subtropical Gyre, north of O'ahu island, Hawaii, with two manipulation experiments. Infra-red epifluorescence microscopy documented that AAP bacteria represented approximately 2% of total bacteria in the euphotic zone with the maximum abundance in the upper 50 m. They conducted active photosynthetic electron transport with maximum rates up to 50 electrons per reaction center per second. The in situ decline of bacteriochlorophyll concentration over the daylight period, an estimate of loss rates due to predation, indicated that the AAP bacteria in the upper 50 m of the water column turned over at rates of 0.75-0.90 d-1. This corresponded well with the specific growth rate determined in dilution experiments where AAP bacteria grew at a rate 1.05 ± 0.09 d-1. An amendment of inorganic nitrogen to obtain N:P = 32 resulted in a more than 10 times increase in AAP abundance over 6 days. The presented data document that AAP bacteria are an active part of the bacterioplankton community in the oligotrophic North Pacific Subtropical Gyre and that their growth was mostly controlled by nitrogen availability and grazing pressure.IMPORTANCEMarine bacteria represent a complex assembly of species with different physiology, metabolism, and substrate preferences. We focus on a specific functional group of marine bacteria called aerobic anoxygenic phototrophs. These photoheterotrophic organisms require organic carbon substrates for growth, but they can also supplement their metabolic needs with light energy captured by bacteriochlorophyll. These bacteria have been intensively studied in coastal regions, but rather less is known about their distribution, growth, and mortality in the oligotrophic open ocean. Therefore, we conducted a suite of measurements in the North Pacific Subtropical Gyre to determine the distribution of these organisms in the water column and their growth and mortality rates. A nutrient amendment experiment showed that aerobic anoxygenic phototrophs were limited by inorganic nitrogen. Despite this, they grew more rapidly than average heterotrophic bacteria, but their growth was balanced by intense grazing pressure.

Marine picoplankton metagenomes and MAGs from eleven vertical profiles obtained by the Malaspina Expedition.

The Ocean microbiome has a crucial role in Earth's biogeochemical cycles. During the last decade, global cruises such as Tara Oceans and the Malaspina Expedition have expanded our understanding of the diversity and genetic repertoire of marine microbes. Nevertheless, there are still knowledge gaps regarding their diversity patterns throughout depth gradients ranging from the surface to the deep ocean. Here we present a dataset of 76 microbial metagenomes (MProfile) of the picoplankton size fraction (0.2-3.0 µm) collected in 11 vertical profiles covering contrasting ocean regions sampled during the Malaspina Expedition circumnavigation (7 depths, from surface to 4,000 m deep). The MProfile dataset produced 1.66 Tbp of raw DNA sequences from which we derived: 17.4 million genes clustered at 95% sequence similarity (M-GeneDB-VP), 2,672 metagenome-assembled genomes (MAGs) of Archaea and Bacteria (Malaspina-VP-MAGs), and over 100,000 viral genomic sequences. This dataset will be a valuable resource for exploring the functional and taxonomic connectivity between the photic and bathypelagic tropical and sub-tropical ocean, while increasing our general knowledge of the Ocean microbiome.

Disentangling microbial networks across pelagic zones in the tropical and subtropical global ocean.

Microbial interactions are vital in maintaining ocean ecosystem function, yet their dynamic nature and complexity remain largely unexplored. Here, we use association networks to investigate possible ecological interactions in the marine microbiome among archaea, bacteria, and picoeukaryotes throughout different depths and geographical regions of the tropical and subtropical global ocean. Our findings reveal that potential microbial interactions change with depth and geographical scale, exhibiting highly heterogeneous distributions. A few potential interactions were global, meaning they occurred across regions at the same depth, while 11-36% were regional within specific depths. The bathypelagic zone had the lowest proportion of global associations, and regional associations increased with depth. Moreover, we observed that most surface water associations do not persist in deeper ocean layers despite microbial vertical dispersal. Our work contributes to a deeper understanding of the tropical and subtropical global ocean interactome, which is essential for addressing the challenges posed by global change.

Global biogeography of the smallest plankton across ocean depths.

Tiny ocean plankton (picoplankton) are fundamental for the functioning of the biosphere, but the ecological mechanisms shaping their biogeography were partially understood. Comprehending whether these microorganisms are structured by niche versus neutral processes is relevant in the context of global change. We investigate the ecological processes (selection, dispersal, and drift) structuring global-ocean picoplanktonic communities inhabiting the epipelagic (0 to 200 meters), mesopelagic (200 to 1000 meters), and bathypelagic (1000 to 4000 meters) zones. We found that selection decreased, while dispersal limitation increased with depth, possibly due to differences in habitat heterogeneity and dispersal barriers such as water masses and bottom topography. Picoplankton ß-diversity positively correlated with environmental heterogeneity and water mass variability, but this relationship tended to be weaker for eukaryotes than for prokaryotes. Community patterns were more pronounced in the Mediterranean Sea, probably because of its cross-basin environmental heterogeneity and deep-water isolation. We conclude that different combinations of ecological mechanisms shape the biogeography of the ocean microbiome across depths.

Top abundant deep ocean heterotrophic bacteria can be retrieved by cultivation.

Traditional culture techniques usually retrieve a small fraction of the marine microbial diversity, which mainly belong to the so-called rare biosphere. However, this paradigm has not been fully tested at a broad scale, especially in the deep ocean. Here, we examined the fraction of heterotrophic bacterial communities in photic and deep ocean layers that could be recovered by culture-dependent techniques at a large scale. We compared 16S rRNA gene sequences from a collection of 2003 cultured heterotrophic marine bacteria with global 16S rRNA metabarcoding datasets (16S TAGs) covering surface, mesopelagic and bathypelagic ocean samples that included 16 of the 23 samples used for isolation. These global datasets represent 60 322 unique 16S amplicon sequence variants (ASVs). Our results reveal a significantly higher proportion of isolates identical to ASVs in deeper ocean layers reaching up to 28% of the 16S TAGs of the bathypelagic microbial communities, which included the isolation of 3 of the top 10 most abundant 16S ASVs in the global bathypelagic ocean, related to the genera Sulfitobacter, Halomonas and Erythrobacter. These isolates contributed differently to the prokaryotic communities across different plankton size fractions, recruiting between 38% in the free-living fraction (0.2-0.8 µm) and up to 45% in the largest particles (20-200 µm) in the bathypelagic ocean. Our findings support the hypothesis that sinking particles in the bathypelagic act as resource-rich habitats, suitable for the growth of heterotrophic bacteria with a copiotroph lifestyle that can be cultured, and that these cultivable bacteria can also thrive as free-living bacteria.

Expanding success in the isolation of abundant marine bacteria after reduction in grazing and viral pressure and increase in nutrient availability.

Isolation of microorganisms is a useful approach to gathering knowledge about their genomic properties, physiology, and ecology, in addition to allowing the characterization of novel taxa. We performed an extensive isolation effort on samples from seawater manipulation experiments that were carried out during the four astronomical seasons in a coastal site of the northwest Mediterranean to evaluate the impact of grazing, viral mortality, resource competition reduction, and light presence/absence on bacterioplankton growth. Isolates were retrieved using two growth media, and their full 16S rRNA gene was sequenced to assess their identity and calculate their culturability across seasons and experimental conditions. A total of 1,643 isolates were obtained, mainly affiliated to the classes Gammaproteobacteria (44%), Alphaproteobacteria (26%), and Bacteroidia (17%). Isolates pertaining to class Gammaproteobacteria were the most abundant in all experiments, while Bacteroidia were preferentially enriched in the treatments with reduced grazing. Sixty-one isolates had a similarity below 97% to cultured taxa and are thus putatively novel. Comparison of isolate sequences with 16S rRNA gene amplicon sequences from the same samples showed that the percentage of reads corresponding to isolates was 21.4% within the whole data set, with dramatic increases in the summer virus-reduced (71%) and diluted (47%) treatments. In fact, we were able to isolate the top 10 abundant taxa in several experiments and from the whole data set. We also show that top-down and bottom-up controls differentially affect taxa in terms of culturability. Our results indicate that culturing marine bacteria using agar plates can be successful in certain ecological situations. IMPORTANCE Bottom-up and top-down controls greatly influence marine microbial community composition and dynamics, which in turn have effects on their culturability. We isolated a high amount of heterotrophic bacterial strains from experiments where seawater environmental conditions had been manipulated and found that decreasing grazing and viral pressure as well as rising nutrient availability are key factors increasing the success in culturing marine bacteria. Our data hint at factors influencing culturability and underpin bacterial cultures as a powerful way to discover new taxa.

Shifts in bacterioplankton community structure between dry and wet seasons in a tropical estuary strongly affected by riverine discharge.

Estuaries are among the most productive ecosystems in the world and are highly dynamic due to the interaction of freshwater and seawater, which results in strong spatial gradients in physico-chemical conditions. Bacterioplankton play a central role in these systems, driving the fluxes of carbon and energy, and being central for contaminant removal in human-impacted areas. Most studies on bacterioplankton dynamics have been carried out in temperate estuaries, and they show that salinity is a major factor driving bacterioplankton distribution. Tropical estuaries, although largely understudied, experience drastic variations in river discharge between the dry and the rainy seasons, influencing the spatial distribution of the salinity gradient and thus likely impacting bacterioplankton communities. Using Illumina sequencing of the 16S rRNA gene, here we studied bacterial communities from the Nicoya's Gulf (Costa Rica), a large tropical estuary characterized by high riverine discharge during the rainy season, to explore seasonal changes in the spatial distribution and connectivity of these communities along the Gulf. Our results show pronounced differences in bacterial diversity and community structure between seasons and zones within the estuary (the shallow upper Gulf, the middle zone and the lower zone, located in the marine end of the estuary). Bacterial communities from the different regions were more similar during the rainy season, suggesting a larger degree of microbial connectivity likely driven by the fast water circulation fueled by the riverine discharge. In the dry season, Enterobacteriales and Cyanobacteria dominated bacterial communities, whereas in the rainy season Alphaproteobacteria was the dominant group. These contrasting seasonal trends were consistent with the seasonal variations observed in bacterial assemblages during a year at a single station in the upper region of the Gulf. We conclude that the Gulf is highly dynamic in both the spatial and temporal scale and that bacterioplankton communities are strongly influenced by the riverine and tidal inputs during both seasons. This study sheds light on the sources of variability in the structure of bacterial communities in tropical estuarine systems, an understudied type of aquatic ecosystem, and sets the basis to design further comprehensive studies on their microbial diversity.

Water mass age structures the auxiliary metabolic gene content of free-living and particle-attached deep ocean viral communities.

BACKGROUND: Viruses play important roles in the ocean's biogeochemical cycles. Yet, deep ocean viruses are one of the most under-explored fractions of the global biosphere. Little is known about the environmental factors that control the composition and functioning of their communities or how they interact with their free-living or particle-attached microbial hosts. RESULTS: We analysed 58 viral communities associated with size-fractionated free-living (0.2-0.8 µm) and particle-attached (0.8-20 µm) cellular metagenomes from bathypelagic (2150-4018 m deep) microbiomes obtained during the Malaspina expedition. These metagenomes yielded 6631 viral sequences, 91% of which were novel, and 67 represented high-quality genomes. Taxonomic classification assigned 53% of the viral sequences to families of tailed viruses from the order Caudovirales. Computational host prediction associated 886 viral sequences to dominant members of the deep ocean microbiome, such as Alphaproteobacteria (284), Gammaproteobacteria (241), SAR324 (23), Marinisomatota (39), and Chloroflexota (61). Free-living and particle-attached viral communities had markedly distinct taxonomic composition, host prevalence, and auxiliary metabolic gene content, which led to the discovery of novel viral-encoded metabolic genes involved in the folate and nucleotide metabolisms. Water mass age emerged as an important factor driving viral community composition. We postulated this was due to changes in quality and concentration of dissolved organic matter acting on the host communities, leading to an increase of viral auxiliary metabolic genes associated with energy metabolism among older water masses. CONCLUSIONS: These results shed light on the mechanisms by which environmental gradients of deep ocean ecosystems structure the composition and functioning of free-living and particle-attached viral communities. Video Abstract.

A Metagenomic and Amplicon Sequencing Combined Approach Reveals the Best Primers to Study Marine Aerobic Anoxygenic Phototrophs.

Studies based on protein-coding genes are essential to describe the diversity within bacterial functional groups. In the case of aerobic anoxygenic phototrophic (AAP) bacteria, the pufM gene has been established as the genetic marker for this particular functional group, although available primers are known to have amplification biases. We review here the existing primers for pufM gene amplification, design new ones, and evaluate their phylogenetic coverage. We then use samples from contrasting marine environments to evaluate their performance. By comparing the taxonomic composition of communities retrieved with metagenomics and with different amplicon approaches, we show that the commonly used PCR primers are biased towards the Gammaproteobacteria phylum and some Alphaproteobacteria clades. The metagenomic approach, as well as the use of other combinations of the existing and newly designed primers, show that these groups are in fact less abundant than previously observed, and that a great proportion of pufM sequences are affiliated to uncultured representatives, particularly in the open ocean. Altogether, the framework developed here becomes a better alternative for future studies based on the pufM gene and, additionally, serves as a reference for primer evaluation of other functional genes.

Disentangling temporal associations in marine microbial networks.

BACKGROUND: Microbial interactions are fundamental for Earth's ecosystem functioning and biogeochemical cycling. Nevertheless, they are challenging to identify and remain barely known. Omics-based censuses are helpful in predicting microbial interactions through the statistical inference of single (static) association networks. Yet, microbial interactions are dynamic and we have limited knowledge of how they change over time. Here, we investigate the dynamics of microbial associations in a 10-year marine time series in the Mediterranean Sea using an approach inferring a time-resolved (temporal) network from a single static network. RESULTS: A single static network including microbial eukaryotes and bacteria was built using metabarcoding data derived from 120 monthly samples. For the decade, we aimed to identify persistent, seasonal, and temporary microbial associations by determining a temporal network that captures the interactome of each individual sample. We found that the temporal network appears to follow an annual cycle, collapsing, and reassembling when transiting between colder and warmer waters. We observed higher association repeatability in colder than in warmer months. Only 16 associations could be validated using observations reported in literature, underlining our knowledge gap in marine microbial ecological interactions. CONCLUSIONS: Our results indicate that marine microbial associations follow recurrent temporal dynamics in temperate zones, which need to be accounted for to better understand the functioning of the ocean microbiome. The constructed marine temporal network may serve as a resource for testing season-specific microbial interaction hypotheses. The applied approach can be transferred to microbiome studies in other ecosystems. Video Abstract.

Impact of particle flux on the vertical distribution and diversity of size-fractionated prokaryotic communities in two East Antarctic polynyas.

Antarctic polynyas are highly productive open water areas surrounded by ice where extensive phytoplankton blooms occur, but little is known about how these surface blooms influence carbon fluxes and prokaryotic communities from deeper waters. By sequencing the 16S rRNA gene, we explored the vertical connectivity of the prokaryotic assemblages associated with particles of three different sizes in two polynyas with different surface productivity, and we linked it to the magnitude of the particle export fluxes measured using thorium-234 (234Th) as particle tracer. Between the sunlit and the mesopelagic layers (700 m depth), we observed compositional changes in the prokaryotic communities associated with the three size-fractions, which were mostly dominated by Flavobacteriia, Alphaproteobacteria, and Gammaproteobacteria. Interestingly, the vertical differences between bacterial communities attached to the largest particles decreased with increasing 234Th export fluxes, indicating a more intense downward transport of surface prokaryotes in the most productive polynya. This was accompanied by a higher proportion of surface prokaryotic taxa detected in deep particle-attached microbial communities in the station with the highest 234Th export flux. Our results support recent studies evidencing links between surface productivity and deep prokaryotic communities and provide the first evidence of sinking particles acting as vectors of microbial diversity to depth in Antarctic polynyas, highlighting the direct influence of particle export in shaping the prokaryotic communities of mesopelagic waters.

Abiotic selection of microbial genome size in the global ocean.

Strong purifying selection is considered a major evolutionary force behind small microbial genomes in the resource-poor photic ocean. However, very little is currently known about how the size of prokaryotic genomes evolves in the global ocean and whether patterns reflect shifts in resource availability in the epipelagic and relatively stable deep-sea environmental conditions. Using 364 marine microbial metagenomes, we investigate how the average genome size of uncultured planktonic prokaryotes varies across the tropical and polar oceans to the hadal realm. We find that genome size is highest in the perennially cold polar ocean, reflecting elongation of coding genes and gene dosage effects due to duplications in the interior ocean microbiome. Moreover, the rate of change in genome size due to temperature is 16-fold higher than with depth up to 200 m. Our results demonstrate how environmental factors can influence marine microbial genome size selection and ecological strategies of the microbiome.

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.

A robust approach to estimate relative phytoplankton cell abundances from metagenomes.

Phytoplankton account for >45% of global primary production, and have an enormous impact on aquatic food webs and on the entire Earth System. Their members are found among prokaryotes (cyanobacteria) and multiple eukaryotic lineages containing chloroplasts. Genetic surveys of phytoplankton communities generally consist of PCR amplification of bacterial (16S), nuclear (18S) and/or chloroplastic (16S) rRNA marker genes from DNA extracted from environmental samples. However, our appreciation of phytoplankton abundance or biomass is limited by PCR-amplification biases, rRNA gene copy number variations across taxa, and the fact that rRNA genes do not provide insights into metabolic traits such as photosynthesis. Here, we targeted the photosynthetic gene psbO from metagenomes to circumvent these limitations: the method is PCR-free, and the gene is universally and exclusively present in photosynthetic prokaryotes and eukaryotes, mainly in one copy per genome. We applied and validated this new strategy with the size-fractionated marine samples collected by Tara Oceans, and showed improved correlations with flow cytometry and microscopy than when based on rRNA genes. Furthermore, we revealed unexpected features of the ecology of these ecosystems, such as the high abundance of picocyanobacterial aggregates and symbionts in the ocean, and the decrease in relative abundance of phototrophs towards the larger size classes of marine dinoflagellates. To facilitate the incorporation of psbO in molecular-based surveys, we compiled a curated database of >18,000 unique sequences. Overall, psbO appears to be a promising new gene marker for molecular-based evaluations of entire phytoplankton communities.

Functional responses of key marine bacteria to environmental change - toward genetic counselling for coastal waters.

Coastal ecosystems deteriorate globally due to human-induced stress factors, like nutrient loading and pollution. Bacteria are critical to marine ecosystems, e.g., by regulating nutrient cycles, synthesizing vitamins, or degrading pollutants, thereby providing essential ecosystem services ultimately affecting economic activities. Yet, until now bacteria are overlooked both as mediators and indicators of ecosystem health, mainly due to methodological limitations in assessing bacterial ecosystem functions. However, these limitations are largely overcome by the advances in molecular biology and bioinformatics methods for characterizing the genetics that underlie functional traits of key bacterial populations - "key" in providing important ecosystem services, being abundant, or by possessing high metabolic rates. It is therefore timely to analyze and define the functional responses of bacteria to human-induced effects on coastal ecosystem health. We posit that categorizing the responses of key marine bacterial populations to changes in environmental conditions through modern microbial oceanography methods will allow establishing the nascent field of genetic counselling for our coastal waters. This requires systematic field studies of linkages between functional traits of key bacterial populations and their ecosystem functions in coastal seas, complemented with systematic experimental analyses of the responses to different stressors. Research and training in environmental management along with dissemination of results and dialogue with societal actors are equally important to ensure the role of bacteria is understood as fundamentally important for coastal ecosystems. Using the responses of microorganisms as a tool to develop genetic counselling for coastal ecosystems can ultimately allow for integrating bacteria as indicators of environmental change.