Single-cell vs. bulk activity properties of coastal bacterioplankton over an annual cycle in a temperate ecosystem

Spatial and seasonal variability of picoplankton abundance and growth rates in the southern Bay of Biscay

Autotrophic and heterotrophic picoplankton play fundamental roles in marine food webs and biogeochemical cycles. However, their growth responses have seldom been jointly assessed, including many temperate regions such as the Bay of Biscay. There, previous studies have shown their relevance in carbon fluxes. We describe here the spatio-temporal variability of the abundances and growth rates of the picoplanktonic groups routinely distinguished by flow cytometry (Synechococcus and Prochlorococcus cyanobacteria, two groups of differently sized picoeukaryotes and two groups of heterotrophic bacteria distinguished by their relative nucleic acid content) in the central Cantabrian Sea (S Bay of Biscay). To that end, from February to December 2021 we collected surface water on 5 occasions from 6 stations distributed along the S Bay of Biscay (6-3°W) and incubated it after removing protistan grazers in order to determine their dynamics along the seasonal cycle as well as the inshore-offshore and the west-east gradients. Seasonal variations in initial and maximum abundances generally matched previous knowledge of the region but growth rates were more variable, with Prochlorococcus and high nucleic acid (HNA) bacteria showing the maximum values (up to 2 d-1) while negative growth was observed in one third of Synechococcus incubations. Temporal differences generally overrode differences along the inshore-offshore gradient in trophic status while in situ and maximum abundances of most of the groups generally decreased towards the east following the increase in stratification and lower nutrient availability. Responses to stratification suggest Prochlorococcus and low nucleic acid (LNA) cells may prevail among autotrophic and heterotrophic bacteria, respectively, in a warmer ocean.

Weekly variations of viruses and heterotrophic nanoflagellates and their potential impact on bacterioplankton in shallow waters of the central Red Sea

Bacterioplankton play a pivotal role in marine ecosystems. However, their temporal dynamics and underlying control mechanisms are poorly understood in tropical regions such as the Red Sea. Here, we assessed the impact of bottom-up (resource availability) and top-down (viruses and heterotrophic nanoflagellates) controls on bacterioplankton abundances by weekly sampling a coastal central Red Sea site in 2017. We monitored microbial abundances by flow cytometry together with a set of environmental variables including temperature, salinity, dissolved organic and inorganic nutrients and chlorophyll a. We distinguished five groups of heterotrophic bacteria depending on their physiological properties relative nucleic acid content, membrane integrity and cell-specific respiratory activity, two groups of Synechococcus cyanobacteria and three groups of viruses. Viruses controlled heterotrophic bacteria for most of the year, as supported by a negative correlation between their respective abundances and a positive one between bacterial mortality rates and mean viral abundances. On the contrary, heterotrophic nanoflagellates abundance covaried with that of heterotrophic bacteria. Heterotrophic nanoflagellates showed preference for larger bacteria from both the high and low nucleic acid content groups. Our results demonstrate that top-down control is fundamental in keeping heterotrophic bacterioplankton abundances low (< 5 × 10 ⁵ cells mL⁻¹) in Red Sea coastal waters.

Changes in Activity and Community Composition Shape Bacterial Responses to Size-Fractionated Marine DOM

To study the response of bacteria to different size-fractions of naturally occurring dissolved organic matter (DOM), a natural prokaryotic community from North Atlantic mesopelagic waters (1000 m depth) was isolated and grown in (i) 0.1-μm filtered seawater (CONTROL), (ii) the low-molecular-weight (<1 kDa) DOM fraction (L-DOM), and (iii) the recombination of high- (>1 kDa) and low-molecular-weight DOM fractions (H + L-DOM), to test the potential effect of ultrafiltration on breaking the DOM size continuum. Prokaryotic abundance and leucine incorporation were consistently higher in the H + L-DOM niche than in the L-DOM and CONTROL treatments, suggesting a different interaction with each DOM fraction and the disruption of the structural DOM continuum by ultrafiltration, respectively. Rhodobacterales (Alphaproteobacteria) and Flavobacteriales (Bacteroidetes) were particularly enriched in L-DOM and closely related to the colored DOM (CDOM) fraction, indicating the tight link between these groups and changes in DOM aromaticity. Conversely, some other taxa that were rare or undetectable in the original bacterial community were enriched in the H + L-DOM treatment (e.g., Alteromonadales belonging to Gammaproteobacteria), highlighting the role of the rare biosphere as a seed bank of diversity against ecosystem disturbance. The relationship between the fluorescence of protein-like CDOM and community composition of populations in the H + L-DOM treatment suggested their preference for labile DOM. Conversely, the communities growing on the L-DOM niche were coupled to humic-like CDOM, which may indicate their ability to degrade more reworked DOM and/or the generation of refractory substrates (as by-products of the respiration processes). Most importantly, L- and/or H + L-DOM treatments stimulated the growth of unique bacterial amplicon sequence variants (ASVs), suggesting the potential of environmental selection (i.e., changes in DOM composition and availability), particularly in the light of climate change scenarios. Taken together, our results suggest that different size-fractions of DOM induced niche-specialization and differentiation of mesopelagic bacterial communities.

Light supports cell-integrity and growth rates of taxonomically diverse coastal photoheterotrophs

Despite the widespread distribution of proteorhodopsin (PR)-containing bacteria in the oceans, the use of light-derived energy to promote bacterial growth has only been shown in a few bacterial isolates, and there is a paucity of data describing the metabolic effects of light on environmental photoheterotrophic taxa. Here, we assessed the effects of light on the taxonomic composition, cell integrity and growth responses of microbial communities in monthly incubations between spring and autumn under different environmental conditions. The photoheterotrophs expressing PR in situ were dominated by Pelagibacterales and SAR116 in July and November, while members of Euryarchaeota, Gammaproteobacteria and Bacteroidetes dominated the PR expression in spring. Cell-membrane integrity decreased under dark conditions throughout most of the assessment, with maximal effects in summer, under low-nutrient conditions. A positive effect of light on growth was observed in one incubation (out of nine), coinciding with a declining phytoplankton bloom. Light-enhanced growth was found in Gammaproteobacteria (Alteromonadales) and Bacteroidetes (Polaribacter and Tenacibaculum). Unexpectedly, some Pelagibacterales also exhibited higher growth rates under light conditions. We propose that the energy harvested by PRs helps to maintain cell viability in dominant coastal photoheterotrophic oligotrophs while promoting the growth of some widespread taxa benefiting from the decline of phytoplankton blooms.

Responses of physiological groups of tropical heterotrophic bacteria to temperature and dissolved organic matter additions: food matters more than warming

Compared to higher latitudes, tropical heterotrophic bacteria may be less responsive to warming because of strong bottom‐up control. In order to separate both drivers, we determined the growth responses of bacterial physiological groups to temperature after adding dissolved organic matter (DOM) from mangroves, seagrasses and glucose to natural seawater from the Great Barrier Reef. Low (LNA) and high (HNA) nucleic acid content, membrane‐intact (Live) and membrane‐damaged (Dead) plus actively respiring (CTC+) cells were monitored for 4 days. Specific growth rates of the whole community were significantly higher (1.9 day^‐1) in the mangrove treatment relative to the rest (0.2–0.4 day^‐1) at in situ temperature and their temperature dependence, estimated as activation energy, was also consistently higher. Strong bottom‐up control was suggested in the other treatments. Cell size depended more on DOM than temperature. Mangrove DOM resulted in significantly higher contributions of Live, HNA and CTC+ cells to total abundance, while the seagrass leachate reduced Live cells below 50%. Warming significantly decreased Live and CTC+ cells contributions in most treatments. Our results suggest that only in the presence of highly labile compounds, such as mangroves DOM, can we anticipate increases in heterotrophic bacteria biomass in response to warming in tropical regions.

Regulation of Low and High Nucleic Acid Fluorescent Heterotrophic Prokaryote Subpopulations and Links to Viral-Induced Mortality Within Natural Prokaryote-Virus Communities

Flow cytometric analysis of marine prokaryotes routinely reveals two distinct clusters of heterotrophic cells referred to as high nucleic acid fluorescent (HNA) and low nucleic acid fluorescent (LNA) populations. Evidence suggests that these may represent physiologically and ecologically distinct prokaryote populations. According to the "kill the winner" hypothesis, viral lysis reduces the efficiency of the microbial loop by decreasing the biomass and activity of the most abundant and active members of a population (i.e., competition specialist). Thus, viral-induced mortality may vary according to the physiology of HNA and LNA cells, with implications for the marine carbon cycle. Here, the abundance and production of heterotrophic prokaryotic populations were assessed in the North Atlantic during two phases of the annual plankton cycle and related to bottom-up (i.e., organic carbon variability) and top-down processes (i.e., viral abundance and lytic production). Our results demonstrate that the relative abundance of HNA and LNA heterotrophic cells and heterotrophic prokaryote production vary according to organic carbon variability in the water column, which can be strongly influenced by the physical eddy field (i.e., type of eddy: cyclonic, anticyclonic, or no eddy). In addition, the abundance and lytic production of virus subpopulations were correlated with  the cellular production and abundance of heterotrophic HNA and LNA prokaryote communities. Our data suggest group- and activity-specific linkages between hosts and viruses (i.e., HNA-V1 and LNA-V2). Specifically, V1 had a greater contribution to total viral production (i.e., 2.6-fold higher than V2 viruses), similar to their putative host. Finally, we explore potential implications of group- and activity-specific linkages between host and virus groups on the flux of carbon through the microbial food web.

Low Abundances but High Growth Rates of Coastal Heterotrophic Bacteria in the Red Sea

Characterized by some of the highest naturally occurring sea surface temperatures, the Red Sea remains unexplored regarding the dynamics of heterotrophic prokaryotes. Over 16 months, we used flow cytometry to characterize the abundance and growth of four physiological groups of heterotrophic bacteria: membrane-intact (Live), high and low nucleic acid content (HNA and LNA) and actively respiring (CTC+) cells in shallow coastal waters. Chlorophyll a, dissolved organic matter (DOC and DON) concentrations, and their fluorescent properties were also measured as proxies of bottom-up control. We performed short-term incubations (6 days) with the whole microbial community (Community treatment), and with the bacterial community only after removing predators by filtration (Filtered treatment). Initial bacterial abundances ranged from 1.46 to 4.80 × 105 cells mL-1. Total specific growth rates in the Filtered treatment ranged from 0.76 to 2.02 d-1. Live and HNA cells displayed similar seasonal patterns, with higher values during late summer and fall (2.13 and 2.33 d-1, respectively) and lower in late spring (1.02 and 1.01 d-1, respectively). LNA cells were outgrown by the other physiological groups (0.33-1.08 d-1) while CTC+ cells (0.28-1.85 d-1) showed weaker seasonality. The Filtered treatment yielded higher bacterial abundances than the Community treatment in all but 2 of the incubations, and carrying capacities peaked in November 2016 (1.04 × 106 cells mL-1), with minimum values (3.61 × 105 cells mL-1) observed in May 2017. The high temperatures experienced from May through October 2016 (33.4 ± 0.4°C) did not constrain the growth of heterotrophic bacteria. Indeed, bacterial growth efficiencies were positively correlated with environmental temperature, reflecting the presence of more labile compounds (high DON concentrations resulting in lower C:N ratios) in summer. The overall high specific growth rates and the consistently higher carrying capacities in the Filtered treatment suggest that strong top-down control by protistan grazers was the likely cause for the low heterotrophic bacteria abundances.

Temperature dependences of growth rates and carrying capacities of marine bacteria depart from metabolic theoretical predictions

Using the metabolic theory of ecology (MTE) framework, we evaluated over a whole annual cycle the monthly responses to temperature of the growth rates (μ) and carrying capacities (K) of heterotrophic bacterioplankton at a temperate coastal site. We used experimental incubations spanning 6ºC with bacterial physiological groups identified by flow cytometry according to membrane integrity (live), nucleic acid content (HNA and LNA) and respiratory activity (CTC+). The temperature dependence of μ at the exponential phase of growth was summarized by the activation energy (E), which was variable (-0.52 to 0.72 eV) but followed a seasonal pattern, only reaching the hypothesized value for aerobic heterotrophs of 0.65 eV during the spring bloom for the most active bacterial groups (live, HNA, CTC+). K (i.e. maximum experimental abundance) peaked at 4 × 10(6) cells mL(-1) and generally covaried with μ but, contrary to MTE predictions, it did not decrease consistently with temperature. In the case of live cells, the responses of μ and K to temperature were positively correlated and related to seasonal changes in substrate availability, indicating that the responses of bacteria to warming are far from homogeneous and poorly explained by MTE at our site.

Community analysis of high- and low-nucleic acid-containing bacteria in NW Mediterranean coastal waters using 16S rDNA pyrosequencing

The ecological significance of the marine bacterial populations distinguishable by flow cytometry on the basis of the fluorescence (FL) of their nucleic acid (NA) content and proxies of cell size (such as side scatter, SSC) remains largely unknown. Some studies have suggested that cells with high NA (HNA) content and high SSC (HS) represent the active members of the community, while the low NA (LNA) cells are inactive members of the same phylogenetic groups. But group-specific activity measurements and phylogenetic assignment after cell sorting have suggested this is not be the case, particularly in open-ocean communities. To test the extent to which the different NA subgroups are similar, and consequently the extent to which they likely have similar ecological and biogeochemical roles in the environment, we analysed the phylogenetic composition of three populations after cell sorting [high NA-high SC (HNA-HS), high NA-low SC (HNA-LS), low NA (LNA)] by 454 pyrosequencing in two contrasting periods of the year in NW Mediterranean coastal waters (BBMO, Blanes Bay Microbial Observatory) where these three populations have recurrent seasonal patterns. Statistical analyses showed that summer and winter samples were significantly different and, importantly, the sorted populations within a sample were composed of different taxa. The majority of taxa were associated with one NA fraction only, and the degree of overlap (i.e. OTUs present simultaneously in 2 fractions) between HNA and LNA and between summer and winter communities was very small. Rhodobacterales, SAR116 and Bacteroidetes contributed primarily to the HNA fraction, whereas other groups such as SAR11 and SAR86 contributed largely to the LNA fractions. Gammaproteobacteria other than SAR86 showed less preference for one particular NA fraction. An increase in diversity was observed from the LNA to the HNA-HS fraction for both sample dates. Our results suggest that, in Blanes Bay, flow cytometric signatures of natural communities track their phylogenetic composition.

Relationships between coastal bacterioplankton growth rates and biomass production: comparison of leucine and thymidine uptake with single-cell physiological characteristics

Specific growth rates of heterotrophic bacterioplankton have been frequently estimated from in situ bacterial production (BP) to biomass (BB) ratios, using a series of assumptions that may result in serious discrepancies with values obtained from predator-free cultures. Here, we used both types of approaches together with a comprehensive assessment of single-cell physiological characteristics (membrane integrity, nucleic acid content, and active respiration) of coastal bacterioplankton during a complete annual cycle (February 2007-January 2008) in the southern Bay of Biscay off Xixón, Spain. Both leucine and thymidine incorporation rates were used in conjunction with empirical tracer to carbon or cells conversion factors (eCFs) to accurately derive BP. Leu and TdR incorporation rates covaried year-round, as did the corresponding eCFs at 0 and 50 m depth. eCFs peaked in autumn, with mean annual values close to the theoretical ones (3.4 kg C mol Leu(-1) and 2.0 × 10(18) cells mol TdR(-1)). Bacterial abundance (0.2-1.5 × 10(6) cells L(-1)) showed a bimodal distribution with maxima in May and October and minima in March. Live (membrane-intact) cells dominated year-round (79-97%), with high nucleic acid cells (42-88%) and actively respiring bacteria (CTC+, 1-16%) showing distinct surface maxima in April and July, respectively. BB (557-1,558 mg C m(-2)) and BP (7-139 mg C m(-2) day(-1)) presented two distinct peaks in spring and autumn, both of similar size due to a strong upwelling event observed in September. Specific growth rates (0.35-3.8 day(-1)) were one order of magnitude higher in predator-free incubations than bacterial turnover rates derived from integrated BP:BB ratios (0.01-0.16 and 0.01-0.09 day(-1), for Leu and TdR, respectively) and were not correlated, probably due to a significant contribution of low activity cells to total standing stocks. The Leu:TdR molar ratio averaged for the water column (6.6-25.5) decreased significantly with higher integrated BB, indicating that low standing stocks tend to present unbalanced growth. Discrepancies about the true magnitude of specific growth rates must be solved before fully appreciating the role of bacteria in the ocean carbon cycle.