Multiyear analysis uncovers coordinated seasonality in stocks and composition of the planktonic food web in the Baltic Sea proper.

The planktonic realm from bacteria to zooplankton provides the baseline for pelagic aquatic food webs. However, multiple trophic levels are seldomly included in time series studies, hampering a holistic understanding of the influence of seasonal dynamics and species interactions on food web structure and biogeochemical cycles. Here, we investigated plankton community composition, focusing on bacterio-, phyto- and large mesozooplankton, and how biotic and abiotic factors correlate at the Linnaeus Microbial Observatory (LMO) station in the Baltic Sea from 2011 to 2018. Plankton communities structures showed pronounced dynamic shifts with recurring patterns. Summarizing the parts of the planktonic microbial food web studied here to total carbon, a picture emerges with phytoplankton consistently contributing > 39% while bacterio- and large mesozooplankton contributed ~ 30% and ~ 7%, respectively, during summer. Cyanophyceae, Actinobacteria, Bacteroidetes, and Proteobacteria were important groups among the prokaryotes. Importantly, Dinophyceae, and not Bacillariophyceae, dominated the autotrophic spring bloom whereas Litostomatea (ciliates) and Appendicularia contributed significantly to the consumer entities together with the more traditionally observed mesozooplankton, Copepoda and Cladocera. Our findings of seasonality in both plankton composition and carbon stocks emphasize the importance of time series analyses of food web structure for characterizing the regulation of biogeochemical cycles and appropriately constraining ecosystem models.

Naturally induced biphasic phytoplankton spring bloom reveals rapid and distinct substrate and bacterial community dynamics.

Phytoplankton spring blooms are typical features in coastal seas and provide heterotrophic bacteria with a rich blend of dissolved substrates. However, they are difficult to study in coastal seas in-situ. Here, we induced a phytoplankton spring bloom and followed its fate for 37 days in four 600 L-mesocosms. To specifically investigate the significance of phytoplankton-born dissolved organic carbon (DOC) we used artificial seawater with low DOC background and inoculated it with a 100 µm-prefiltered plankton community from the North Sea. A biphasic bloom developed, dominated by diatoms and Phaeocystis globosa respectively. In between, bacterial numbers peaked, followed by a peak in virus-like particles, implying that virus infection caused the collapse. Concentrations of dissolved free amino acids exhibited rapid changes, in particular during the diatom bloom and until the peak in bacterial abundance. Dissolved combined amino acids and neutral monosaccharides accumulated continuously, accounting for 22% of DOC as a mean and reaching levels as high as 44%. Bacterial communities were largely dominated by Bacteroidetes, especially the NS3a marine group (family Flavobacteriaceae), but Rhodobacteraceae and Gammaproteobacteria were also prominent members. Our study shows rapid organic matter and community composition dynamics that are hard to trace in natural coastal ecosystems.

Warmer temperatures favor slower-growing bacteria in natural marine communities.

Earth's life-sustaining oceans harbor diverse bacterial communities that display varying composition across time and space. While particular patterns of variation have been linked to a range of factors, unifying rules are lacking, preventing the prediction of future changes. Here, analyzing the distribution of fast- and slow-growing bacteria in ocean datasets spanning seasons, latitude, and depth, we show that higher seawater temperatures universally favor slower-growing taxa, in agreement with theoretical predictions of how temperature-dependent growth rates differentially modulate the impact of mortality on species abundances. Changes in bacterial community structure promoted by temperature are independent of variations in nutrients along spatial and temporal gradients. Our results help explain why slow growers dominate at the ocean surface, during summer, and near the tropics and provide a framework to understand how bacterial communities will change in a warmer world.

Quantification of metabolic niche occupancy dynamics in a Baltic Sea bacterial community.

Progress in molecular methods has enabled the monitoring of bacterial populations in time. Nevertheless, understanding community dynamics and its links with ecosystem functioning remains challenging due to the tremendous diversity of microorganisms. Conceptual frameworks that make sense of time series of taxonomically rich bacterial communities, regarding their potential ecological function, are needed. A key concept for organizing ecological functions is the niche, the set of strategies that enable a population to persist and define its impacts on the surroundings. Here we present a framework based on manifold learning to organize genomic information into potentially occupied bacterial metabolic niches over time. Manifold learning tries to uncover low-dimensional data structures in high-dimensional data sets that can be used to describe the data in reduced dimensions. We apply the method to re-construct the dynamics of putatively occupied metabolic niches using a long-term bacterial time series from the Baltic Sea, the Linnaeus Microbial Observatory (LMO). The results reveal a relatively low-dimensional space of occupied metabolic niches comprising groups of taxa with similar functional capabilities. Time patterns of occupied niches were strongly driven by seasonality. Some metabolic niches were dominated by one bacterial taxon, whereas others were occupied by multiple taxa, depending on the season. These results illustrate the power of manifold learning approaches to advance our understanding of the links between community composition and functioning in microbial systems. IMPORTANCE The increase in data availability of bacterial communities highlights the need for conceptual frameworks to advance our understanding of these complex and diverse communities alongside the production of such data. To understand the dynamics of these tremendously diverse communities, we need tools to identify overarching strategies and describe their role and function in the ecosystem in a comprehensive way. Here, we show that a manifold learning approach can coarse grain bacterial communities in terms of their metabolic strategies and that we can thereby quantitatively organize genomic information in terms of potentially occupied niches over time. This approach, therefore, advances our understanding of how fluctuations in bacterial abundances and species composition can relate to ecosystem functions and it can facilitate the analysis, monitoring, and future predictions of the development of microbial communities.

Spatial and biological oceanographic insights into the massive fish-killing bloom of the haptophyte Chrysochromulina leadbeateri in northern Norway.

A bloom of the fish-killing haptophyte Chrysochromulina leadbeateri in northern Norway during May and June 2019 was the most harmful algal event ever recorded in the region, causing massive mortalities of farmed salmon. Accordingly, oceanographic and biodiversity aspects of the bloom were studied in unprecedented detail, based on metabarcoding and physico-chemical and biotic factors related with the dynamics and distribution of the bloom. Light- and electron-microscopical observations of nanoplankton samples from diverse locations confirmed that C. leadbeateri was dominant in the bloom and the primary cause of associated fish mortalities. Cell counts by light microscopy and flow cytometry were obtained throughout the regional bloom within and adjacent to five fjord systems. Metabarcoding sequences of the V4 region of the 18S rRNA gene from field material collected during the bloom and a cultured isolate from offshore of Tromsøy island confirmed the species identification. Sequences from three genetic markers (18S, 28S rRNA gene and ITS region) verified the close if not identical genetic similarity to C. leadbeateri from a previous massive fish-killing bloom in 1991 in northern Norway. The distribution and cell abundance of C. leadbeateri and related Chrysochromulina species in the recent incident were tracked by integrating observations from metabarcoding sequences of the V4 region of the 18S rRNA gene. Metabarcoding revealed at least 14 distinct Chrysochromulina variants, including putative cryptic species. C. leadbeateri was by far the most abundant of these species, but with high intraspecific genetic variability. Highest cell abundance of up to 2.7 × 107 cells L - 1 of C. leadbeateri was found in Balsfjorden; the high cell densities were associated with stratification near the pycnocline (at ca. 12 m depth) within the fjord. The cell abundance of C. leadbeateri showed positive correlations with temperature, negative correlation with salinity, and a slightly positive correlation with ambient phosphate and nitrate concentrations. The spatio-temporal succession of the C. leadbeateri bloom suggests independent initiation from existing pre-bloom populations in local zones, perhaps sustained and supplemented over time by northeastward advection of the bloom from the fjords.

Seasonal Dynamics in Carbon Cycling of Marine Bacterioplankton Are Lifestyle Dependent.

Although free-living (FL) and particle-attached (PA) bacteria are recognized as ecologically distinct compartments of marine microbial food-webs, few, if any, studies have determined their dynamics in abundance, function (production, respiration and substrate utilization) and taxonomy over a yearly cycle. In the Baltic Sea, abundance and production of PA bacteria (defined as the size-fraction >3.0 µm) peaked over 3 months in summer (6 months for FL bacteria), largely coinciding with blooms of Chitinophagales (Bacteroidetes). Pronounced changes in the growth efficiency (range 0.05-0.27) of FL bacteria (defined as the size-fraction <3.0 µm) indicated the magnitude of seasonal variability of ecological settings bacteria experience. Accordingly, 16S rRNA gene analyses of bacterial community composition uncovered distinct correlations between taxa, environmental variables and metabolisms, including Firmicutes associated with elevated hydrolytic enzyme activity in winter and Verrucomicrobia with utilization of algal-derived substrates during summer. Further, our results suggested a substrate-controlled succession in the PA fraction, from Bacteroidetes using polymers to Actinobacteria and Betaproteobacteria using monomers across the spring to autumn phytoplankton bloom transition. Collectively, our findings emphasize pronounced seasonal changes in both the composition of the bacterial community in the PA and FL size-fractions and their contribution to organic matter utilization and carbon cycling. This is important for interpreting microbial ecosystem function-responses to natural and human-induced environmental changes.

Auxiliary Metabolic Gene Functions in Pelagic and Benthic Viruses of the Baltic Sea.

Marine microbial communities are facing various ecosystem fluctuations (e.g., temperature, organic matter concentration, salinity, or redox regimes) and thus have to be highly adaptive. This might be supported by the acquisition of auxiliary metabolic genes (AMGs) originating from virus infections. Marine bacteriophages frequently contain AMGs, which allow them to augment their host's metabolism or enhance virus fitness. These genes encode proteins for the same metabolic functions as their highly similar host homologs. In the present study, we analyzed the diversity, distribution, and composition of marine viruses, focusing on AMGs to identify their putative ecologic role. We analyzed viruses and assemblies of 212 publicly available metagenomes obtained from sediment and water samples across the Baltic Sea. In general, the virus composition in both compartments differed compositionally. While the predominant viral lifestyle was found to be lytic, lysogeny was more prevalent in sediments than in the pelagic samples. The highest proportion of AMGs was identified in the genomes of Myoviridae. Overall, the most abundantly occurring AMGs are encoded for functions that protect viruses from degradation by their hosts, such as methylases. Additionally, some detected AMGs are known to be involved in photosynthesis, 7-cyano-7-deazaguanine synthesis, and cobalamin biosynthesis among other functions. Several AMGs that were identified in this study were previously detected in a large-scale analysis including metagenomes from various origins, i.e., different marine sites, wastewater, and the human gut. This supports the theory of globally conserved core AMGs that are spread over virus genomes, regardless of host or environment.

Seasonal and Spatial Variations in Synechococcus Abundance and Diversity Throughout the Gullmar Fjord, Swedish Skagerrak.

The picophytoplankton Synechococcus is a globally abundant autotroph that contributes significantly to primary production in the oceans and coastal areas. These cyanobacteria constitute a diverse genus of organisms that have developed independent niche spaces throughout aquatic environments. Here, we use the 16S V3-V4 rRNA gene region and flow cytometry to explore the diversity of Synechococcus within the picophytoplankton community in the Gullmar Fjord, on the west coast of Sweden. We conducted a station-based 1-year time series and two transect studies of the fjord. Our analysis revealed that within the large number of Synechococcus amplicon sequence variants (ASVs; 239 in total), prevalent ASVs phylogenetically clustered with clade representatives in both marine subcluster 5.1 and 5.2. The near-surface composition of ASVs shifted from spring to summer, when a 5.1 subcluster dominated community developed along with elevated Synechococcus abundances up to 9.3 × 104 cells ml-1. This seasonal dominance by subcluster 5.1 was observed over the length of the fjord (25 km), where shifts in community composition were associated with increasing depth. Unexpectedly, the community shift was not associated with changes in salinity. Synechococcus abundance dynamics also differed from that of the photosynthetic picoeukaryote community. These results highlight how seasonal variations in environmental conditions influence the dynamics of Synechococcus clades in a high latitude threshold fjord.

Connectivity of Fennoscandian Shield terrestrial deep biosphere microbiomes with surface communities.

The deep biosphere is an energy constrained ecosystem yet fosters diverse microbial communities that are key in biogeochemical cycling. Whether microbial communities in deep biosphere groundwaters are shaped by infiltration of allochthonous surface microorganisms or the evolution of autochthonous species remains unresolved. In this study, 16S rRNA gene amplicon analyses showed that few groups of surface microbes infiltrated deep biosphere groundwaters at the Äspö Hard Rock Laboratory, Sweden, but that such populations constituted up to 49% of the microbial abundance. The dominant persisting phyla included Patescibacteria, Proteobacteria, and Epsilonbacteraeota. Despite the hydrological connection of the Baltic Sea with the studied groundwaters, infiltrating microbes predominantly originated from deep soil groundwater. Most deep biosphere groundwater populations lacked surface representatives, suggesting that they have evolved from ancient autochthonous populations. We propose that deep biosphere groundwater communities in the Fennoscandian Shield consist of selected infiltrated and indigenous populations adapted to the prevailing conditions.

Taxon-Specific Shifts in Bacterial and Archaeal Transcription of Dissolved Organic Matter Cycling Genes in a Stratified Fjord.

A considerable fraction of organic matter derived from photosynthesis in the euphotic zone settles into the ocean's interior and, as it progresses, is degraded by diverse microbial consortia that utilize a suite of extracellular enzymes and membrane transporters. Still, the molecular details that regulate carbon cycling across depths remain little explored. As stratification in fjords has made them attractive models to explore patterns in biological oceanography, we here analyzed bacterial and archaeal transcription in samples from five depth layers in the Gullmar Fjord, Sweden. Transcriptional variation over depth correlated with gradients in chlorophyll a and nutrient concentrations. Differences in transcription between sampling dates (summer and early autumn) were strongly correlated with ammonium concentrations, which potentially was linked with a stronger influence of (micro-)zooplankton grazing in summer. Transcriptional investment in carbohydrate-active enzymes (CAZymes) decreased with depth and shifted toward peptidases, partly a result of elevated CAZyme transcription by Flavobacteriales, Cellvibrionales, and Synechococcales at 2 to 25 m and a dominance of peptidase transcription by Alteromonadales and Rhodobacterales from 50 m down. In particular, CAZymes for chitin, laminarin, and glycogen were important. High levels of transcription of ammonium transporter genes by Thaumarchaeota at depth (up to 18% of total transcription), along with the genes for ammonia oxidation and CO2 fixation, indicated that chemolithoautotrophy contributed to the carbon flux in the fjord. The taxon-specific expression of functional genes for processing of the marine pool of dissolved organic matter and inorganic nutrients across depths emphasizes the importance of different microbial foraging mechanisms over spatiotemporal scales for shaping biogeochemical cycles. IMPORTANCE It is generally recognized that stratification in the ocean strongly influences both the community composition and the distribution of ecological functions of microbial communities, which in turn are expected to shape the biogeochemical cycling of essential elements over depth. Here, we used metatranscriptomics analysis to infer molecular detail on the distribution of gene systems central to the utilization of organic matter in a stratified marine system. We thereby uncovered that pronounced shifts in the transcription of genes encoding CAZymes, peptidases, and membrane transporters occurred over depth among key prokaryotic orders. This implies that sequential utilization and transformation of organic matter through the water column is a key feature that ultimately influences the efficiency of the biological carbon pump.

Composition and Seasonality of Membrane Transporters in Marine Picoplankton.

In this study, we examined transporter genes in metagenomic and metatranscriptomic data from a time-series survey in the temperate marine environment of the Baltic Sea. We analyzed the abundance and taxonomic distribution of transporters in the 3µm-0.2µm size fraction comprising prokaryotes and some picoeukaryotes. The presence of specific transporter traits was shown to be guiding the succession of these microorganisms. A limited number of taxa were associated with the dominant transporter proteins that were identified for the nine key substrate categories for microbial growth. Throughout the year, the microbial taxa at the level of order showed highly similar patterns in terms of transporter traits. The distribution of transporters stayed the same, irrespective of the abundance of each taxon. This would suggest that the distribution pattern of transporters depends on the bacterial groups being dominant at a given time of the year. Also, we find notable numbers of secretion proteins that may allow marine bacteria to infect and kill prey organisms thus releasing nutrients. Finally, we demonstrate that transporter proteins may provide clues to the relative importance of biogeochemical processes, and we suggest that virtual transporter functionalities may become important components in future population dynamics models.

Sweet spheres: succession and CAZyme expression of marine bacterial communities colonizing a mix of alginate and pectin particles.

Polysaccharide particles are important substrates and microhabitats for marine bacteria. However, substrate-specific bacterial dynamics in mixtures of particle types with different polysaccharide composition, as likely occurring in natural habitats, are undescribed. Here, we studied the composition, functional diversity and gene expression of marine bacterial communities colonizing a mix of alginate and pectin particles. Amplicon, metagenome and metatranscriptome sequencing revealed that communities on alginate and pectin particles significantly differed from their free-living counterparts. Unexpectedly, microbial dynamics on alginate and pectin particles were similar, with predominance of amplicon sequence variants (ASVs) from Tenacibaculum, Colwellia, Psychrobium and Psychromonas. Corresponding metagenome-assembled genomes (MAGs) expressed diverse alginate lyases, several colocalized in polysaccharide utilization loci. Only a single, low-abundant MAG showed elevated transcript abundances of pectin-degrading enzymes. One specific Glaciecola ASV dominated the free-living fraction, possibly persisting on particle-derived oligomers through different glycoside hydrolases. Elevated ammonium uptake and metabolism signified nitrogen as an important factor for degrading carbon-rich particles, whereas elevated methylcitrate and glyoxylate cycles suggested nutrient limitation in surrounding waters. The bacterial preference for alginate, whereas pectin primarily served as colonization scaffold, illuminates substrate-driven dynamics within mixed polysaccharide pools. These insights expand our understanding of bacterial niche specialization and the biological carbon pump in macroalgae-rich habitats.

Ecosystem-wide metagenomic binning enables prediction of ecological niches from genomes.

The genome encodes the metabolic and functional capabilities of an organism and should be a major determinant of its ecological niche. Yet, it is unknown if the niche can be predicted directly from the genome. Here, we conduct metagenomic binning on 123 water samples spanning major environmental gradients of the Baltic Sea. The resulting 1961 metagenome-assembled genomes represent 352 species-level clusters that correspond to 1/3 of the metagenome sequences of the prokaryotic size-fraction. By using machine-learning, the placement of a genome cluster along various niche gradients (salinity level, depth, size-fraction) could be predicted based solely on its functional genes. The same approach predicted the genomes' placement in a virtual niche-space that captures the highest variation in distribution patterns. The predictions generally outperformed those inferred from phylogenetic information. Our study demonstrates a strong link between genome and ecological niche and provides a conceptual framework for predictive ecology based on genomic data.

High Frequency Multi-Year Variability in Baltic Sea Microbial Plankton Stocks and Activities.

Marine bacterioplankton are essential in global nutrient cycling and organic matter turnover. Time-series analyses, often at monthly sampling frequencies, have established the paramount role of abiotic and biotic variables in structuring bacterioplankton communities and productivities. However, fine-scale seasonal microbial activities, and underlying biological principles, are not fully understood. We report results from four consecutive years of high-frequency time-series sampling in the Baltic Proper. Pronounced temporal dynamics in most investigated microbial variables were observed, including bacterial heterotrophic production, plankton biomass, extracellular enzyme activities, substrate uptake rate constants of glucose, pyruvate, acetate, amino acids, and leucine, as well as nutrient limitation bioassays. Spring blooms consisting of diatoms and dinoflagellates were followed by elevated bacterial heterotrophic production and abundances. During summer, bacterial productivity estimates increased even further, coinciding with an initial cyanobacterial bloom in early July. However, bacterial abundances only increased following a second cyanobacterial bloom, peaking in August. Uptake rate constants for the different measured carbon compounds varied seasonally and inter-annually and were highly correlated to bacterial productivity estimates, temperature, and cyanobacterial abundances. Further, we detected nutrient limitation in response to environmental conditions in a multitude of microbial variables, such as elevated productivities in nutrient bioassays, changes in enzymatic activities, or substrate preferences. Variations among biotic variables often occurred on time scales of days to a few weeks, yet often spanning several sampling occasions. Such dynamics might not have been captured by sampling at monthly intervals, as compared to more predictable transitions in abiotic variables such as temperature or nutrient concentrations. Our study indicates that high resolution analyses of microbial biomass and productivity parameters can help out in the development of biogeochemical and food web models disentangling the microbial black box.

Marine Bacterioplankton Seasonal Succession Dynamics.

Bacterioplankton (bacteria and archaea) are indispensable regulators of global element cycles owing to their unique ability to decompose and remineralize dissolved organic matter. These microorganisms in surface waters worldwide exhibit pronounced seasonal succession patterns, governed by physicochemical factors (e.g., light, climate, and nutrient loading) that are determined by latitude and distance to shore. Moreover, we emphasize that the effects of large-scale factors are modulated regionally, and over shorter timespans (days to weeks), by biological interactions including molecule exchanges, viral lysis, and grazing. Thus the interplay and scaling between factors ultimately determine the success of particular bacterial populations. Spatiotemporal surveys of bacterioplankton community composition provide the necessary frame for interpreting how the distinct metabolisms encoded in the genomes of different bacteria regulate biogeochemical cycles.

Spatio-Temporal Interdependence of Bacteria and Phytoplankton during a Baltic Sea Spring Bloom.

In temperate systems, phytoplankton spring blooms deplete inorganic nutrients and are major sources of organic matter for the microbial loop. In response to phytoplankton exudates and environmental factors, heterotrophic microbial communities are highly dynamic and change their abundance and composition both on spatial and temporal scales. Yet, most of our understanding about these processes comes from laboratory model organism studies, mesocosm experiments or single temporal transects. Spatial-temporal studies examining interactions of phytoplankton blooms and bacterioplankton community composition and function, though being highly informative, are scarce. In this study, pelagic microbial community dynamics (bacteria and phytoplankton) and environmental variables were monitored during a spring bloom across the Baltic Proper (two cruises between North Germany to Gulf of Finland). To test to what extent bacterioplankton community composition relates to the spring bloom, we used next generation amplicon sequencing of the 16S rRNA gene, phytoplankton diversity analysis based on microscopy counts and population genotyping of the dominating diatom Skeletonema marinoi. Several phytoplankton bloom related and environmental variables were identified to influence bacterial community composition. Members of Bacteroidetes and Alphaproteobacteria dominated the bacterial community composition but the bacterial groups showed no apparent correlation with direct bloom related variables. The less abundant bacterial phyla Actinobacteria, Planctomycetes, and Verrucomicrobia, on the other hand, were strongly associated with phytoplankton biomass, diatom:dinoflagellate ratio, and colored dissolved organic matter (cDOM). Many bacterial operational taxonomic units (OTUs) showed high niche specificities. For example, particular Bacteroidetes OTUs were associated with two distinct genetic clusters of S. marinoi. Our study revealed the complexity of interactions of bacterial taxa with inter- and intraspecific genetic variation in phytoplankton. Overall, our findings imply that biotic and abiotic factors during spring bloom influence bacterial community dynamics in a hierarchical manner.

Stimulation of growth by proteorhodopsin phototrophy involves regulation of central metabolic pathways in marine planktonic bacteria.

Proteorhodopsin (PR) is present in half of surface ocean bacterioplankton, where its light-driven proton pumping provides energy to cells. Indeed, PR promotes growth or survival in different bacteria. However, the metabolic pathways mediating the light responses remain unknown. We analyzed growth of the PR-containing Dokdonia sp. MED134 (where light-stimulated growth had been found) in seawater with low concentrations of mixed [yeast extract and peptone (YEP)] or single (alanine, Ala) carbon compounds as models for rich and poor environments. We discovered changes in gene expression revealing a tightly regulated shift in central metabolic pathways between light and dark conditions. Bacteria showed relatively stronger light responses in Ala compared with YEP. Notably, carbon acquisition pathways shifted toward anaplerotic CO2 fixation in the light, contributing 31 ± 8% and 24 ± 6% of the carbon incorporated into biomass in Ala and YEP, respectively. Thus, MED134 was a facultative double mixotroph, i.e., photo- and chemotrophic for its energy source and using both bicarbonate and organic matter as carbon sources. Unexpectedly, relative expression of the glyoxylate shunt genes (isocitrate lyase and malate synthase) was >300-fold higher in the light--but only in Ala--contributing a more efficient use of carbon from organic compounds. We explored these findings in metagenomes and metatranscriptomes and observed similar prevalence of the glyoxylate shunt compared with PR genes and highest expression of the isocitrate lyase gene coinciding with highest solar irradiance. Thus, regulatory interactions between dissolved organic carbon quality and central metabolic pathways critically determine the fitness of surface ocean bacteria engaging in PR phototrophy.