ORIGINAL RESEARCH article
Sec. Aquatic Microbiology
Viral-Mediated Microbe Mortality Modulated by Ocean Acidification and Eutrophication: Consequences for the Carbon Fluxes Through the Microbial Food Web
- 1Biological Oceanography Laboratory, Austral Center for Scientific Research (CONICET), Ushuaia, Argentina
- 2Department of Marine Biology and Oceanography, Institut de Ciències del Mar (CSIC), Barcelona, Spain
- 3Center for Marine Ecosystems Research, School of Sciences, Edith Cowan University, Joondalup, WA, Australia
- 4Institució Catalana de Recerca i Estudis Avançats (ICREA), Barcelona, Spain
- 5Centre for Ecology and Evolution in Microbial Model Systems, Linnaeus University, Kalmar, Sweden
Anthropogenic carbon emissions are causing changes in seawater carbonate chemistry including a decline in the pH of the oceans. While its aftermath for calcifying microbes has been widely studied, the effect of ocean acidification (OA) on marine viruses and their microbial hosts is controversial, and even more in combination with another anthropogenic stressor, i.e., human-induced nutrient loads. In this study, two mesocosm acidification experiments with Mediterranean waters from different seasons revealed distinct effects of OA on viruses and viral-mediated prokaryotic mortality depending on the trophic state and the successional stage of the plankton community. In the winter bloom situation, low fluorescence viruses, the most abundant virus-like particle (VLP) subpopulation comprising mostly bacteriophages, were negatively affected by lowered pH with nutrient addition, while the bacterial host abundance was stimulated. High fluorescence viruses, containing cyanophages, were stimulated by OA regardless of the nutrient conditions, while cyanobacteria of the genus Synechococcus were negatively affected by OA. Moreover, the abundance of very high fluorescence viruses infecting small haptophytes tended to be lower under acidification while their putative hosts' abundance was enhanced, suggesting a direct and negative effect of OA on viral–host interactions. In the oligotrophic summer situation, we found a stimulating effect of OA on total viral abundance and the viral populations, suggesting a cascading effect of the elevated pCO2 stimulating autotrophic and heterotrophic production. In winter, viral lysis accounted for 30 ± 16% of the loss of bacterial standing stock per day (VMMBSS) under increased pCO2 compared to 53 ± 35% in the control treatments, without effects of nutrient additions while in summer, OA had no significant effects on VMMBSS (35 ± 20% and 38 ± 5% per day in the OA and control treatments, respectively). We found that phage production and resulting organic carbon release rates significantly reduced under OA in the nutrient replete winter situation, but it was also observed that high nutrient loads lowered the negative effect of OA on viral lysis, suggesting an antagonistic interplay between these two major global ocean stressors in the Anthropocene. In summer, however, viral-mediated carbon release rates were lower and not affected by lowered pH. Eutrophication consistently stimulated viral production regardless of the season or initial conditions. Given the relevant role of viruses for marine carbon cycling and the biological carbon pump, these two anthropogenic stressors may modulate carbon fluxes through their effect on viruses at the base of the pelagic food web in a future global change scenario.
Current anthropogenic carbon release rates are unprecedented over the last 66 million years (Zeebe et al., 2016). Since preindustrial time until the present, the concentration of carbon dioxide (CO2) in the atmosphere has increased from ~277 ppm (e.g., Joos and Spahni, 2008) to 410 ppm in 2019 (Dlugokencky and Tans, 2020) due to the burning of fossil fuels, cement manufacturing, and land use changes. The ocean has absorbed about 30% of these CO2 emissions since 1750, leading to changes in seawater carbonate chemistry and a decline in pH (Sabine et al., 2004). Preindustrial pH values in the surface ocean of about 8.2 are expected to decrease to 7.8 at the end of the twenty-first century in a high-CO2-emission, “business-as-usual” scenario (Bopp et al., 2013), a process referred to as ocean acidification (OA, Caldeira and Wickett, 2003). During the past century, anthropogenic inputs of nitrogen (N) and phosphorus (P) to coastal ecosystems via river discharge have augmented mostly due to the use of fertilizers in agriculture (Galloway et al., 2004). Increased nutrient and organic matter supply to coastal environments has led to coastal eutrophication, which constitutes a major threat to ocean health (Rabalais et al., 2009). Thus, in addition to OA, eutrophication is another major global change pressure that directly affects marine ecosystems in the Anthropocene (Diaz and Rosenberg, 2008; Doney et al., 2009). Synergistic (Cai et al., 2011) or antagonistic interplays between both processes have been suggested (Borges and Gypens, 2010; Malone and Newton, 2020).
The microbial food web, composed of nano- and microalgae, autotrophic picoplankton, heterotrophic protists, heterotrophic bacteria and archaea, and virioplankton, constitutes the base of marine food webs and is a key component in the transfer of carbon and the regeneration of micro- and macronutrients (Azam, 1998). It has been estimated that ~50% of the global primary production is generated in the oceans by phytoplankton photosynthesis (Field et al., 1998). Heterotrophic prokaryotes contribute greatly to the degradation and transformation of organic matter through their extensive metabolic activity and act as an intermediate link between the dissolved organic carbon (DOC) pool and higher trophic levels through the microbial loop (Azam et al., 1983). On the other hand, marine viruses play a pivotal role in the carbon cycle of marine planktonic food webs (Suttle, 2005). Through killing the hosts and releasing the content of the cytoplasm to the environment, viruses transfer particulate organic matter (POM) to the dissolved organic matter (DOM) pool (Wilhelm and Suttle, 1999; Sheik et al., 2014). Mathematical models have shown that this so-called viral shunt has a negative effect on zooplankton production through its link to the microbial loop, especially under oligotrophic regimes (Murray and Eldridge, 1994). However, by fueling the DOM pool, viral activity also stimulates bacterial production (Malits and Weinbauer, 2009) with synergetic effects on the microbial food web (Berdjeb et al., 2011), boosts the recycling of nutrients (Gobler et al., 1997; Poorvin et al., 2004), and stimulates autotrophic production (Weinbauer et al., 2011) and bacterial respiration (Fuhrman, 1999; Middelboe and Lyck, 2002; Bonilla-Findji et al., 2008). Overall, viral lysis shifts the food web toward a more regenerative pathway, reducing the carbon flux to higher trophic levels (Brussaard et al., 2008b). Along with the lytic viral reproduction where the phage injects its genetic material into the host in order to redirect its metabolism toward the production of new phages and cell lysis, the genome of the phage may also remain in the host in a dormant stage and replicate along with the host (lysogenic cycle), until the lytic cycle is induced by stress or other unfavorable conditions (Weinbauer, 2004). The lytic cycle is also induced when bacteria grow actively (Wilson and Mann, 1997), and lysogeny has been shown to be important at high host densities (Knowles et al., 2016). Both life styles (lysis, lysogeny) are involved in maintaining genetic diversity among the host community (Weinbauer and Rassoulzadegan, 2004; Suttle, 2007; Malits and Weinbauer, 2009) having direct implications on the bulk community enzymatic machinery and, consequently, on global biogeochemical cycles.
Ocean microbes are currently exposed to large depth, regional, seasonal, and even daily pH variations driven by biological and physical processes, e.g., CO2 fixation during a phytoplankton bloom increases the pH in the water (Joint et al., 2011). For example, at the ocean time series station ALOHA, in the oligotrophic Central Pacific, the significant long-term decreasing trend of surface water pH by 0.04 pH units from 1988 to 2007 is superimposed to a seasonal surface pH variability of up to 0.06 pH units (Dore et al., 2009). In addition, microbes are exposed to pH values as low as those predicted in the surface oceans for the end of the twenty-first century in the thermocline, where sinking organic matter decomposition by aerobic respiration results in a reduction in pH (Joint et al., 2011). Moreover, in coastal systems, where increased nutrient loading and eutrophication stimulate microbial activity, the observed reductions in pH greatly exceed the values expected from anthropogenic CO2 uptake alone (Provoost et al., 2010). These observations suggest that, even if OA does not lead to dramatic changes in the biogeochemical cycles driven by microbes due to their flexibility to adjust to pH changes (Joint et al., 2011), coastal eutrophication can amplify the effects of pH change and thus affect biochemical processes.
Experimental studies have documented OA-induced changes in the phytoplankton community with shifts toward small picoplankton (Meakin and Wyman, 2011; Brussaard et al., 2013; Spilling et al., 2016; Crawfurd et al., 2017; Schulz et al., 2017), with major implications for the C fluxes in the marine planktonic food webs (Worden et al., 2015). However, mesocosm experiments in the Baltic Sea revealed that, in a picoplankton dominated food web, reduced plankton community respiration rates under high partial pressure of carbon dioxide (pCO2) did not translate into increased carbon export (Spilling et al., 2016).
Mesocosm experiments in the Arctic (Roy et al., 2013; Sperling et al., 2013; Zhang et al., 2013) and North Sea (Newbold et al., 2012; Oliver et al., 2014) did not find evidence of a pCO2 effect on the bulk prokaryotic community composition. However, the combined effects of acidification and nutrient additions (Baltar et al., 2015) or warming (Lindh et al., 2013) did select for specific bacterial phylotypes evidencing the synergistic effects of human-induced perturbations on marine systems. Bacterial bulk production, enzyme activity, and growth rate were indeed enhanced under future atmospheric pCO2 levels (Grossart et al., 2006; Sala et al., 2016).
The effect of changes in pCO2/pH on lytic viral production and viral abundance has been less studied, and the few reports available so far provide contradictory results. Classically studied Escherichia coli bacteriophages such as T2 and T7 are indeed sensitive to changes in pH (Danovaro et al., 2011). Elevated pCO2 did not alter total viral abundance in mesocosm experiments in the North Sea (Rochelle-Newall et al., 2004) and the Arctic Ocean (Brussaard et al., 2013) but increased viral abundances in a Baltic Sea mesocosm study (Tsiola et al., 2017). In a mesocosm experiment, large (i.e., high fluorescence) virus abundance was higher in the control compared to elevated pCO2, and two specific large double-stranded DNA (dsDNA) viruses infecting the haptophytes Emiliania huxleyi and Crysochromulina ericina decreased in abundance with increasing pCO2 levels, suggesting changes in viral diversity (Larsen et al., 2008). Other studies in virus–host systems have shown a negative effect of OA on primary production of Synecchococcus and on the infection by its associated virus (Traving et al., 2013) but no effects on viral lysis of Micromonas pulsilla (Maat et al., 2014) with a potential positive feedback for carbon and nutrient cycling. No discernable effects of enhanced pCO2 on lytic viral production could be detected in other acidification mesocosm experiments in the Baltic Sea (Crawfurd et al., 2017), Mediterranean (Tsiola et al., 2017), and Arctic Ocean (Vaqué et al., 2019), while lysogeny increased with OA in combination with warming in that Arctic Ocean study (Vaqué et al., 2019). Nonetheless, the effect of decreased pH on lysogeny remains uncertain.
Thus, no conclusive data on the effect of ocean acidification on virioplankton and the consequences for microbial food web carbon fluxes exist, even less in combination with eutrophication. In this context, we conducted two mesocosm acidification experiments with microbial communities of the Mediterranean, in which we also added excess nutrients, in two different seasons, winter and summer, to investigate the combined effect of OA and eutrophication on the microbial food web. These experiments already demonstrated a shift to medium-sized phytoplankton under acidified but oligotrophic conditions (Sala et al., 2016), eutrophication effects on dissolved organic matter quality and composition regardless of changes in pH (Aparicio et al., 2016), synergistic effects of nutrient loading and OA on prokaryotic community structure (Baltar et al., 2015), and OA effects on prokaryoplankton gene expression (Bunse et al., 2016). In the present study, we followed viral abundances along with those of their hosts and assessed phage-mediated bacterial mortality and lysogeny in these mesocosm experiments in order to evaluate the potential implications of anthropogenic perturbations of marine systems on viruses and viral-mediated processes in the microbial food web.
Materials and Methods
Two experiments (referred to as WINTER and SUMMER hereafter) were performed with water collected from the Blanes Bay Microbial Observatory (BBMO), NW Mediterranean (Gasol et al., 2016) on February 17, 2010 and on July 6, 2011, respectively. For both experiments, surface water samples were mixed and quickly transferred to eight 200-L polyethylene mesocosms at the Institut de Ciències del Mar aquaria facilities in Barcelona (corresponding to day 0 of the experiments). Four experimental conditions were randomly assigned to duplicated containers: KB (control), KA (lowered pH), NB (nutrient amended), and NA (nutrient amended and lowered pH). Experiments were conducted in a temperature-controlled chamber, set at approximately in situ temperature (Table 1), and a combination of cool-white and grolux lamps was used to illuminate the mesocosms. Measured light intensity inside the containers was 121.3 ± 3.5 μmol m−2 s−1 in the WINTER and 140.8 ± 13.5 μmol m−2 s−1 in the SUMMER experiments. The light/dark cycle was set at 12:12 h. Further details can be found in Sala et al. (2016), Baltar et al. (2015), Aparicio et al. (2016), and Bunse et al. (2016).
Table 1. In situ and experimental conditions in the mesocosms after adding nutrients in experiments WINTER and SUMMER.
Nitrogen and phosphorus were added at Redfield ratio, while silicate was added in excess to ensure diatom growth (the final molar ratio for P/N/Si was ~1:16:30) before the lights were turned on, on February 18 and July 7, for experiments WINTER and SUMMER, respectively (corresponding to day 1). In KA and NA, pH was artificially lowered by bubbling, in a controlled way, small amounts of CO2 (99.9% purity) directly into the mesocosms. This was done every morning in order to maintain the levels of pH in the acidified tanks at around 0.25–0.30 pH units lower than the controls (KB and NB), which were bubbled with equivalent amounts of pure air to ensure similar turbulent conditions. The lowering in pH of the acidified treatments vs. the controls is equivalent to the values projected for the end of the twenty-first century following relatively pessimistic scenarios (Bopp et al., 2013). Seawater pH in the mesocosms was continuously monitored using glass electrodes (LL Ecotrode plus—Metrohm, calibrated every day with a Tris buffer, following standard procedures, Dickson et al., 2007) and recorded by a D130 data logger (Consort, Belgium). In addition, prior to each controlled addition of CO2, we performed precise measurements of the mesocosm's seawater pH using spectrophotometry (Clayton and Byrne, 1993) and alkalinity through a fast, single-point potentiometric titration (Pérez et al., 2000). Initial conditions are summarized in Table 1. To avoid sedimentation, we gently agitated the water during 5 min twice a day. Samples for chlorophyll, microbial abundances, and prokaryotic activity were taken daily (for details see Sala et al., 2016), while experiments for virus-mediated mortality of heterotrophic prokaryotes and the fraction of lytic and temperate viruses were performed at days 1, 5, and 8 for the WINTER and at days 0, 4, and 8 for the SUMMER experiment.
Inorganic Nutrients and Chlorophyll a
Samples for inorganic nutrients were kept frozen at −20°C until analysis, which were performed using a CFA Bran+Luebbe autoanalyzer following the methods described by Hansen and Koroleff (2007).
For chlorophyll a (chl a) analysis according to Yentsch and Menzel (1963), 50 ml of seawater was filtered through Whatman GF/F filters. Pigments were extracted in 90% acetone at 4°C for 24 h and determined by measuring their fluorescence using a Turner Designs fluorometer.
For picophytoplankton, fresh, unstained samples were analyzed in a Becton Dickinson FACSCalibur flow cytometer at high speed (about 100 μl min−1) following Marie et al. (2001). Phototrophic populations (Prochlorococcus, Synechococcus, small and large picoeukaryotes) were discriminated and enumerated according to their light scatter and specific autofluorescence properties.
Since the experiment took place with surface waters that do not harbor almost any Archaea (Alonso-Sáez et al., 2007), we are using the term bacteria to indicate both Bacteria and Archaea. Samples for heterotrophic bacterial abundance (BA) were fixed with paraformaldehyde (1%) and glutaraldehyde (0.05%), kept at room temperature for about 10 min and then flash frozen in liquid nitrogen. Within a few days, the samples were thawed, stained with SYBR Green I (Molecular Probes Inc.) for 10 min, and analyzed in a Becton Dickinson FACSCalibur flow cytometer as described previously (Marie et al., 1997; Gasol and Morán, 2015). Bacteria were determined in plots of 90° light scatter (SSC) vs. green DNA fluorescence (FL1). Differences in FL1 allowed to separate bacteria with low nucleic acid content (LNA) from those with high nucleic acid content (HNA) (Gasol et al., 1999).
For the abundance of virus-like particles (VLPs), samples (1 ml) were fixed with 0.2-μm-filtered glutaraldehyde (0.5% final concentration), incubated at 4°C for 15–30 min, and subsequently frozen in liquid nitrogen and stored at −80°C. Upon thawing, viruses were stained with SYBR Green I (Molecular Probes Inc.) for 10 min in the dark at 80°C and quantified after dilution with TE buffer [10 mM Tris, 1 mM ethylenediaminetetraacetic acid (EDTA), pH = 8] using flow cytometry and an optimized protocol (Marie et al., 1999; Brussaard, 2004). Four VLP populations were distinguished based on their signature in the cytometric plots of side scatter (SSC) vs. green fluorescence (FL1): low, medium, high, and very high green fluorescence VLPs (Supplementary Figure 1). The latter group was identified by its cytometric signature comparable to viruses infecting Pyramimonas orientalis and Phaeocystis pouchetii (Brussaard, 2004). The low and medium green fluorescence subpopulations represent mostly the numerically dominant bacteriophages (Brussaard, 2009). The sample flowrate was accurately calibrated following the protocol of Marie et al. (2001) and used to calculate the abundances of viruses.
Bulk heterotrophic bacterial production (BP) was estimated from 3H-leucine incorporation (Kirchman et al., 1985). For each sample, quadruplicate aliquots (1.2 ml) and two trichloroacetic acid (TCA)-killed controls were incubated with 40 nM 3H-leucine for about 1.5 h in the dark at in situ temperature. Subsequently, leucine incorporation was stopped by adding 120 μl of cold TC 50% to each replicate. Samples were stored at −20°C until processed following published methods (Smith and Azam, 1992). Leucine incorporated into bacterial biomass was converted to bacterial carbon production using the theoretical factor of 1.55 kg C mol−1 Leu that assumes no isotope dilution (Simon and Azam, 1989).
The numbers of actively respiring bacteria were determined using the fluorogenic tetrazolium dye 5-cyano-2,3-ditolyl tetrazolium chloride (CTC) labeling of highly active cells (Sherr et al., 1999; Sieracki et al., 1999). CTC was added to 1-ml subsamples at a final 5 mM concentration from a daily prepared 10× batch and incubated at in situ temperature for 3 h. After incubation, the samples were analyzed with a FACSCalibur flow cytometer as described in detail elsewhere (Gasol and Arístegui, 2007).
Lytic Viral Production, Fraction of Infected Cells and Lysogens, Viral-Mediated Mortality, and C Release Rates
Lytic viral production (VPL), the fraction of infected cells (FIC), induced viral production from lysogens (VPI), and the fraction of lysogenic cells (FLC) were estimated for selected samples using the virus reduction approach (VRA, Weinbauer et al., 2010). The underlying principle of the VRA is to reduce viral abundances in the water, thereby essentially preventing new viral infections. Thus, it is assumed that the viruses produced during the incubation originate from already infected cells. In order to eliminate viruses, bacteria in ~400 ml from one replicate of each treatment were concentrated using a tangential flow system with a peristaltic pump (Watson-Marlow 323) equipped with a 0.2-μm cartridge (VIVAFLOW 200). To obtain virus-free seawater, the 0.2-μm pore-size ultrafiltrate was passed through a 30-kDa cartridge (VIVAFLOW 200). The bacterial concentrates were brought up to the original volume with virus-free seawater and incubated in triplicates in 50-ml Falcon tubes (BD Biosciences) at in situ temperature and in the dark (Table 1) for 24 h. At time 0 of these VRA experiments, three additional tubes were amended with mitomycin C (MC, Sigma) at a final concentration of 1 μg ml−1 in order to induce the lytic cycle of lysogens; untreated samples served as controls (Paul and Weinbauer, 2010). Subsamples (1 ml) for VA and BA from each incubation were taken every 2–4 h during the first 6–12 h and after 24 h, fixed with glutaraldehyde (0.5% final concentration), incubated at 4°C for 15–30 min, subsequently frozen in liquid nitrogen, and stored at −80°C until counted by flow cytometry as described above. VPL was calculated as the increase in viral abundance over short time intervals (~4 h). An increase in viral abundance in the MC treatments represents VPL + VPI (Paul and Weinbauer, 2010; Weinbauer et al., 2010). VPL and VPI were corrected for the changes in the initial BA in the viral production assays with respect to in situ BA. The burst size (BS) was estimated following the approach of Wells and Deming (2006), i.e., dividing the number of viruses produced during the first hour of incubation by the concomitant decline of bacterial abundance. Dividing the number of produced phages by BS yields the number of lysed cells and gives an estimate of the fraction of infected cells (FIC) when divided by BA at the start of the experiments. Dividing the number of induced phages by BS and the BA at t0 gives an estimate of FLC. To obtain the rate of cell lysis, viral production corrected for in situ bacterial abundance was divided by the average estimated BS for WINTER and SUMMER, respectively. Lysis rates were used to calculate virus-mediated mortality of bacteria per day as a percentage of the bacterial standing stock (VMMBSS day−1). Carbon release rates through viral lysis were calculated based on VP by converting the number of lysed bacteria into carbon using a factor of 20 fg C cell−1 (Lee and Fuhrman, 1987).
All statistical analyses were performed with JMP 7.0 (SAS). The Shapiro–Wilk W-test was used to check for normal distribution of data. Analysis of covariance (ANCOVA) with time as a covariate was used to discern the effects of time from those of the lowered pH and nutrient amendment for time-averaged parameters. Kruskal–Wallis tests for non-normal distributions were used to evaluate the acidification effect for K and N treatments, separately. Spearman rank correlation for non-parametric data was performed to determine the relationships between the various measured parameters.
Initial temperature, light conditions, pH, and the concentrations of inorganic nutrients for the WINTER and SUMMER experiments are summarized in Table 1. For a detailed description of the experimental conditions and microplankton and bacterial dynamics in the experiments, see Baltar et al. (2015), Sala et al. (2016), and Aparicio et al. (2016). In WINTER, in situ nutrient concentrations (phosphate and silicate) were four to five times higher than in SUMMER. Initial experimental nitrogen concentrations in the N treatments were higher in the WINTER than in the SUMMER experiment, as they were added by multiplying the monthly average concentration measured in the BBMO during the last 10 years by a factor of 8 (Table 1). Due to the biological activity in the mesocosms, the pH in the experimental chambers of the WINTER experiment tended to increase with time and varied by up to 0.12 U (7.8 ± 0.05 in the KA, 7.9 ± 0.07 in the KB, 7.9 ± 0.12 in the NA, and 8.1 ± 0.12 in the NB treatments, Figure 1A). On the contrary, during SUMMER, the pH tended to decrease and was less variable (7.6 ± 0.08 in the KA, 7.9 ± 0.06 in the KB, 7.6 ± 0.07 in the NA, 7.9 ± 0.05 in the NB treatments, Figure 1B).
Figure 1. Temporal dynamics of (A,B) pH, (C,D) total chlorophyll a, (E,F) abundance of Synechococcus sp., (G,H) small picoeukaryotes, and (I,J) large picoeukaryotes during the experiments WINTER (left) and SUMMER (right). Bars indicate the maximum and minimum values of experimental duplicates; when not visible, they are contained within the symbol. Note the differences in the y-axis scales of chlorophyll a between seasons. KB, control; KA, lowered pH; NB, nutrient amended; NA, nutrient amended and lowered pH.
Chlorophyll a Concentration and Autotrophic Picoplankton
Average chl a concentration was six times higher in the WINTER (6.8 ± 8.9 μg L−1) than in the SUMMER experiment (1.1 ± 0.9 μg L−1, n = 159, P < 0.0001) and increased significantly in the nutrient-amended treatments with time and with respect to the controls (Table 2). Chl a peaked on days 7 and 6 in the WINTER and the SUMMER experiment, respectively, being almost 10 times higher in WINTER than in SUMMER (Figures 1C,D). In WINTER, chl a concentration increased significantly with acidification only in the K treatments (Table 3), while in SUMMER, chl a concentration was significantly stimulated in the low pH treatments, regardless the nutrient conditions (Figure 1D and Table 2).
Table 2. Results of an analysis of covariance (ANCOVA) with time as a covariate showing the significance of the effects of acidification and nutrient amendments during the time course (days 2–9) of the experiments WINTER and SUMMER.
Table 3. Spearman's rank correlation coefficients (ρ) and corresponding significance (P) of measured pH in the experimental mesocosm with other parameters shown for the experiments WINTER (n = 12 for VPL and VMM, n = 24 for CTC+ and n = 36 for the other parameters) and SUMMER (n = 14–15 for VPL and VMM, n = 26 for LNA and HNA, and n = 36 for the other parameters).
Autotrophic picoplankton growth was stimulated by nutrient addition anticipating the chlorophyll peak at least in WINTER (Figure 1 and Table 2). Cynanobacteria of the genus Synechococcus were negatively affected by acidification in WINTER and not affected in SUMMER (Figures 1E,F and Table 2). Prochlorococcus sp. was about one order of magnitude less abundant than Synechococcus sp. and not affected by OA (Supplementary Figure 2 and Table 2). The average abundance of small picoeukaryotes in WINTER (1.1 ± 0.1 × 104 ml−1) doubled that observed in SUMMER (0.6 ± 0.1 × 104 ml−1, Kruskal–Wallis, n = 143, P < 0.01) and increased significantly with nutrient additions in both experiments (Figures 1G,H and Table 2). The abundance of large picoeukaryotes was not significantly different between experiments (4.7 ± 4.4 × 103 and 3.4 ± 2.5 × 103 ml−1 in WINTER and SUMMER, respectively), increased with nutrient addition, and was stimulated by the lowered pH only in SUMMER (Figures 1I,J and Tables 1, 2).
Heterotrophic Bacterial Abundance and Activity
Time-integrated BA was significantly higher in the WINTER (2.3 ± 0.1 × 106 ml−1) than in the SUMMER experiment (0.7 ± 0.1 × 106 ml−1, Kruskal–Wallis, n = 144, P < 0.0001), and BA dynamics differed between seasons: in WINTER, BA increased 7.6 times from day 0 until day 4 without differences between treatments (Figure 2A), while in SUMMER, BA decreased initially and increased from day 3 in all treatments but more significantly in the nutrient-amended ones. An acidification effect was found only for the treatments without nutrient addition (Figure 2B and Table 2).
Figure 2. Temporal dynamics of (A,B) bacterial abundance, (C,D) bacterial heterotrophic production, and (E,F) the percentage of CTC-labeled cells, i.e., actively respiring bacteria during the experiments WINTER (left) and SUMMER (right). Bars indicate the maximum and minimum values of experimental duplicates; when not visible, they are contained within the symbol. Note the differences in the y-axis scales of bacterial abundance and % CTC+ cells between seasons. For treatment codes, see Figure 1.
Heterotrophic bacterial carbon production rates increased in WINTER from 0.26 μg C L−1 day−1 at day 0 to 24.44–29.09 μg C L−1 day−1 until day 3, while they decreased in SUMMER from initial 21.28 ± 2.91 μg C L−1 day−1 and 36.01 ± 5.18 μg C L−1 day−1 in the nutrient unamended and amended treatments, respectively, until days 2–3 (Figures 2C,D). BP was stimulated by nutrient addition but was not affected by lowered pH (Table 2) and was significantly related to the abundance of HNA bacteria except in the SUMMER treatments, where it was related to that of LNA bacteria (Supplementary Table 1). The percentage of CTC+ cells was significantly stimulated by acidification only in the K treatments of WINTER (ANCOVA, n = 40, P < 0.05, Figures 2E,F).
Abundance of the Viral Populations
Total VLP abundance increased after an initial stable phase from day 7 and 5 in experiments WINTER and SUMMER, respectively, and the same behavior was exhibited by the low subpopulation (Figures 3A,B); this increase was significantly higher in the N treatments than in the K treatments (Table 2). Total VLP abundance and the abundance of the MEDIUM subpopulation integrated over the experimental time were significantly lower in the WINTER than in the SUMMER experiment (Figures 3C,D, Kruskal–Wallis, n = 144, P < 0.0001).
Figure 3. Temporal dynamics of the abundance of the different virus-like particle (VLP) subpopulations: (A,B) low, (C,D) medium, (E,F) high, and (G) very high fluorescence VLP during experiments WINTER (left) and SUMMER (right). For treatment codes, see Figure 1.
In experiment WINTER, only the high subpopulation was stimulated by acidification (Figures 3E,F, Table 2). A cytometric group of very high fluorescence VLP (very high, Supplementary Figure 1) emerged in the second half of the WINTER experiment (Figure 3G). The abundance of the very high subpopulation was significantly stimulated by nutrient additions but significantly lower in the acidified treatments toward the end of the experiment (ANCOVA, n = 24, P < 0.05). Total VLP and the subpopulation of low fluorescence VLP, i.e., mostly bacteriophages, were negatively affected by lowered pH with nutrient addition (Table 3).
In the nutrient-deplete SUMMER experiment, in turn, total VLP abundance and the abundances of all subpopulations were correlated negatively to the measured pH in the mesocosms (Table 3), i.e., increased with acidification. For the growth phase of VLP (days 6–9), the lowered pH stimulated the low fluorescence subpopulation, i.e., the bacteriophages, regardless the nutrient addition (Figure 3B, ANCOVA, n = 15, P < 0.0001). In the nutrient treatments, the medium fluorescence subpopulation abundances, also comprising bacteriophages, were triggered in the acidified containers for the time interval 5–7 days (Figure 3D, ANCOVA, P = 0.0011).
Viral Production and Lysogeny
Initial lytic viral production rates (VPL) did not differ significantly between seasons and averaged 0.68 ± 0.32 × 107 VLP ml−1 day−1 but evolved differently between experiments (Figures 4A,B). In WINTER, VPL increased in all treatments but were significantly lower in the KA treatments (2.73 ± 0.65 × 107 VLP ml−1 day−1) with respect to KB (4.41 ± 0.56 × 107 VLP ml−1 day−1, pooled data for days 5 and 8, Kruskal–Wallis, n = 12, P < 0.01). VPL was significantly stimulated by nutrient addition (Kruskal–Wallis, n = 23, P < 0.05) with a trend of reduced VPL due to lowered pH (4.41 ± 1.64 × 107 and 5.63 ± 1.13 × 107 VLP ml−1 day−1 in NA and NB treatments, respectively, Kruskal–Wallis, n = 11, P = 0.05). Lysogenic production could only be induced at day 0 (0.69 ± 0.01 × 107 VLP ml−1 day−1) and in the NA treatments at day 5 (3.46 ± 2.96 × 107 VLP ml−1 day−1, Table 4).
Figure 4. (A,B) Lytic viral production (VPL), (C,D) viral-mediated mortality expressed as a percentage of the bacterial standing stock (VMMBSS), (E,F) carbon release rates (μg C L−1 day−1) through viral lysis calculated from VP from all treatments in the experiments WINTER (left) and SUMMER (right). Note the difference in the y-axis scales between seasons, in particular for VP. Bars indicate standard deviations. Stars indicate significant differences between treatments (P < 0.05). For treatment codes, see Figure 1.
Table 4. The fraction of lysogenic viral production (VPI) over total viral production, the fraction of infected cells (FIC), and the fraction of lysogenic cells (FLC) in the different treatments from the WINTER and SUMMER experiments.
In the SUMMER experiment, initial VPL rates decreased in the K treatments on days 4 and 8 (Figure 4B, Kruskal–Wallis, n = 15, P < 0.05) and were not affected by lowered pH. In the NB treatments, lysogenic production rates accounted for 100% of virus production, i.e., virus production could only be detected in the treatments where lysogens were induced; VPI were 0.30 ± 0.11 × 107 and 0.48 ± 0.13 × 107 VLP ml−1 day−1 at days 4 and 8, respectively, being lower or equal compared to the VPL rates in the acidified, i.e., NA, treatments. Moreover, VPI could only be induced at day 0 without nutrient addition (0.33 ± 0.11 × 107 VLP ml−1 day−1) and at day 4 in the KA treatments (0.19 ± 0.14 × 107 VLP ml−1 day−1, Table 4).
Calculated BS averaged 57 ± 50 in WINTER (15–226) and 16 ± 10 (4–43) in SUMMER. The fraction of infected bacteria (FIC) was low in both experiments and all treatments and accounted for 13.9 ± 18.9% in WINTER and 4.8 ± 3.2% in SUMMER. Only in WINTER, FIC decreased with nutrient additions (Kruskal–Wallis, n = 23, P < 0.05) and was not affected by acidification. The FLC also remained low throughout the experiments and treatments without discernable treatment effects (Table 4).
Viral-Mediated Bacterial Mortality and Carbon Fluxes
The viral-mediated loss of bacterial standing stock (VMMBSS day−1) did not differ between seasons (39 ± 28% in WINTER and 41 ± 22% in SUMMER) but was initially lower in WINTER (21 ± 5%) than in SUMMER (54 ± 30%, Figures 4C,D). VMMBSS was lowered by acidification only in WINTER (Figures 4C,D, Kruskal–Wallis, n = 23, P < 0.05). Carbon released from the particulate to the dissolved pool due to viral activity was four times higher in WINTER (15.0 ± 5.2 μg C L−1 day−1) than in SUMMER (3.9 ± 1.3 μg C L−1 day−1, Kruskal–Wallis, n = 40, P < 0.0001, Figures 4E,F). Only in WINTER viral-mediated carbon flux was significantly and negatively affected by acidification (Figure 4E, Kruskal–Wallis, n = 23, P < 0.05) and also stimulated by nutrient addition (Figure 4, Kruskal–Wallis, n = 23, P < 0.01).
After pooling all data from both seasons and all treatments, viral-mediated carbon flux increased with pH (Spearman ρ = 0.56, P = 0.0002, n = 40) and was positively correlated to chl a concentration and BA (Spearman ρ = 0.67 and 0.75, respectively, P < 0.0001, n = 40). However, when separating by seasons and nutrient treatments, the trend of increasing viral production and viral-mediated mortality with increasing pH was significant only in the K treatments of the WINTER experiment (Table 3). In addition, lytic viral production correlated significantly and positively with bacterial abundance and production only in WINTER without nutrient amendment (Supplementary Tables 1, 2).
Ex situ mesocosm experiments were performed to study the effect of acidification and eutrophication on virioplankton and viral activity under controlled laboratory conditions. Although controlled ex situ experimental conditions may be biased by avoiding processes that could have intervened in the natural environment, they allow to identify potential mechanisms in the response of the viral community to the studied anthropogenic stressors.
The addition of nutrients stimulated viral production regardless of the season or initial conditions. While in the nutrient-replete situation of the WINTER experiment, the abundance of viral populations and bacteriophage production tended to be negatively affected by OA; in the SUMMER experiment, we found a stimulating trend of lowered pH on viral abundance and production. Nutrient availability seems to play a major role in modulating the effect of increasing pCO2 on viruses, but the composition and physiological state of the initial plankton communities also modified OA effects on viral-mediated mortality and on the resulting carbon fluxes.
Effect of OA and Eutrophication on Phytoplankton Viruses
In the WINTER experiment, we found a population of very high fluorescence VLP, as identified by their signature in the cytometric plot of high side scatter and fluorescence, comparable to the cytometric signature of viruses infecting small haptophytes and chlorophytes (Brussaard, 2004). The fact that this population was emerging in the decline phase of the bloom of the “large picoeukaryotes” cytometric group, mostly constituted of small mixotrophic haptophytes (2–4 μm, Massana, pers. comm.), which are the most abundant nanoflagellates at the sampling site year-round (Unrein et al., 2014), suggests that these were actually their potential viruses. Very high fluorescence VLP could not be detected in the SUMMER experiment, where large picoeukaryotes were less abundant. The abundance of these VLP tended to be lower at lowered pH although their putative hosts' abundance was enhanced. Consistently, Larsen et al. (2008) found a decreased production of nanophytoplankton-specific viruses (e.g., infecting Emiliania huxleyi) coupled to increased host production under elevated pCO2, suggesting that lower viral production was not a secondary effect of reduced host production but a direct and negative effect on viral–host interactions. Maat et al. (2014) found viral lysis of M. pulsilla reduced under enhanced pCO2 and P-limited conditions. In addition, during an acidification experiment in the Arctic, very small sized picophytoplankton was stimulated by elevated levels of pCO2, while the viral group associated to the decline in the picophytoplankton bloom and apparently responsible for that decline was not found to be affected by pCO2 (Brussaard et al., 2013).
In WINTER, cyanobacteria of the genus Synechococcus were negatively affected by OA, while the subpopulation of high fluorescence VLP was stimulated regardless of the nutrient conditions (Table 2 and Figure 3E). The high subpopulation is generally constituted by cyanophages and viruses of picoeukaryotes (Brussaard, 2004; Brussaard et al., 2008a; Martinez et al., 2014). The positive correlation of this subpopulation to Synechococcus (Spearman ρ = 0.57, n = 72, P < 0.0001) and to Prochlorococcus (Spearman ρ = 0.41, n = 72, P = 0.0004) lends further support that this group was mainly comprised of cyanophages. Thus, the negative response of Synechococcus in WINTER could be the consequence of enhanced viral lysis in this season. Indeed, freshwater cyanophages have been shown to tolerate a broad range of pH (Suttle, 2000). Moreover, in the SUMMER experiment, where the pH changed over 0.6 U, these putative cyanophages were stimulated by acidification in both nutrient treatments (Table 3). There, however, Synechococcus was stimulated in the OA treatments, at least from day 5 in the nutrient control treatments, where pH reached the lowest values due to microbial activities (Figure 1). Several studies suggest a stimulating effect of increased pCO2 on cyanobacteria production due to energy savings related to the reduced use of carbon concentration mechanisms (Fu et al., 2007; Raven et al., 2012). In the case of the oligotrophic conditions of the SUMMER experiment, the enhanced production of putative cyanophages might appear as a secondary effect of enhanced host production under elevated levels of pCO2. In fact, studies with the model system cyanophage S-PM2 and its host Synechococcus sp. WH7803 revealed a negative effect of OA on the growth rate of Synecchococcus, which compromised biogenesis and the replication cycle of S-PM2 (Traving et al., 2013).
Effect of OA and Eutrophication on Bacteriophage Production
In WINTER, bacterial production increased during the first day, when bacteriophage production was significantly reduced under OA, and this was most apparent in the nutrient control treatments. Since bacterial host abundance was stimulated or not affected by pH in the nutrient amended and control treatments, respectively (Table 3), phage replication rates or phage–host interactions seem to be directly affected by acidification. In fact, metatranscriptome analysis from WINTER found that phage encoding genes most similar to Podoviridae phages infecting SAR11-clade bacteria were negatively affected by acidification, and this was more conspicuous in the non-amended treatments (Bunse et al., 2016).
If bacterial species resistant to viral infection are promoted by ocean acidification, this could affect phage production and total bacterial mortality. High throughput sequencing of 16S ribosomal RNA (rRNA) gen amplicons from the present study showed the synergistic effects of OA and nutrient amendment selected for some specific OTUs belonging to SAR86 (Baltar et al., 2015). However, SAR86 OTUs had only a minor contribution to total bacterial community and might not have been responsible for the reduction in phage-mediated bacterial mortality in the acidification treatments of WINTER (Figure 4C). Moderate changes in pH were shown to cause compositional shifts in the bacterial community from the North Sea in all seasons (Krause et al., 2012). Conversely, most other OA mesocosm experiments showed a lack of clear pCO2 effects on global bacterial community structure (Newbold et al., 2012; Roy et al., 2013; Sperling et al., 2013; Zhang et al., 2013; Oliver et al., 2014), and thus, changes in the bulk bacterial mortality are unlikely caused by the selection of lysis-resistant bacterial species through acidification.
Instead of looking at slight phylogenetic changes in the bacterial community as a response to OA, the metabolic state of the bacterial hosts, which is crucial for successful viral infection, proliferation, and bulk viral production, should be considered. Metatranscriptome analyses from the WINTER experiment showed that bacteria adjusted their gene expression patterns to acidification by enhancing the expression of genes encoding proton pumps, such as respiration complexes, proteorhodopsin, and membrane transporters in the nutrient control treatments (Bunse et al., 2016). Proton transport across membranes demands much energy at the expense of bacterial growth and growth efficiency. Since viral infection and BS rely on the bacterial growth rate and growth efficiency (Middelboe, 2000), viral production rates should decrease with lower bacterial growth efficiency. Actually, VP rates were significantly and positively correlated to bacterial heterotrophic production and abundance in the nutrient deplete controls (Supplementary Figure 1). There, the lowered pH stimulated the percentage of CTC-labeled cells but did not affect bacterial heterotrophic production rates (Figure 2), suggesting lower bacterial growth efficiency under OA as demonstrated from previous experiments (Motegi et al., 2013; Celussi et al., 2017) with implications for bacteriophage production.
Concerning lysogeny as a potential viral strategy to cope with unfavorable conditions (Weinbauer, 2004), we found a very low incidence of the lysogenic life style with no discernible effect of OA (Table 4) in accordance with a previous study from the Mediterranean (Tsiola et al., 2017). Excluding day 0, only 4 out of 12 measurements of viral production showed lysogenic production induced by mitomycin C (MC). Microcosm experiments to discern the effect of the viral life style on bacterial community composition have shown consistently lower viral and bacterial numbers in the MC treatments (Chen et al., 2019), suggesting harmful effects of the inducing reagent itself on bacterial growth, inhibiting certain bacterial clades (Hewson and Fuhrman, 2007). However, in the NB treatments of SUMMER, temperate viruses accounted for the totality of viral production. It has been suggested that lysogeny dominates in oligotrophic environments with low host abundances and contact rates (Weinbauer et al., 2003a) and that prophage induction occurs in rapidly growing cells (Wilson and Mann, 1997). Moreover, protistan grazing also increases bacterial growth rates by recycling nutrients (Simek et al., 2001; Weinbauer et al., 2003b), and this might induce lysogens into the lytic cycle. Unfortunately, grazing rates were not measured in the present study, but the bacterial dynamics and total phytoplankton biomass in SUMMER were very similar to a previous experiment from the same sampling site, where grazing accounted for up to 100% of the bacterial standing stock (Baltar et al., 2016). Interestingly, and similar to the present study, VLP reached considerably high numbers at the end of that experiment (Sandaa et al., 2009).
Modulation of OA Effects on the Microbial Community by Seasonality and Nutrient Availability
The response of the microbial community and their associated viruses to lowered pH differed considerably between seasons. Higher pCO2 did not stimulate auto- and heterotrophic production of the microbial community in a nutrient replete situation in WINTER, while it did stimulate the microbial community from an oligotrophic situation in SUMMER (Figures 1, 2; Sala et al., 2016). Besides the nutrient availability, the seasonal differences were most conspicuous in the different successional stages of the plankton community with a diatom-dominated phytoplankton bloom in winter and pico- and nanoflagellates prevailing in the summer plankton community (Sala et al., 2016).
In our experiments, total VLP abundance and the abundance of VLP subpopulations were intimately linked to microbial production dynamics, in agreement with Larsen et al. (2001, 2004). They were also stimulated consistently by acidification in both nutrient treatments in SUMMER, as a consequence of the cascading effect of OA on phytoplankton and bacterial production in the oligotrophic situation (Table 3). This finding contrasted with the nutrient-replete WINTER situation (Figure 4), where phage-mediated processes were significantly reduced by OA. These differences in the pH sensitivity of viral production and host production is in accordance with results from the Pelagic Ecosystem CO2 Enrichment studies (PeECE I–III). There, bacterial production was either not affected (Allgaier et al., 2008) or stimulated by lowered pH (Grossart et al., 2006). However, these experiments were performed under similar conditions, same site and approximately the same dates, but under different trophic conditions. This was translated into a lower total phytoplankton biomass. Riebesell et al. (2008) argued that the trophic state of the system, i.e., limitation by inorganic and organic nutrients, should influence, at least indirectly, the microbial community.
Indeed, the effects of lowered pH on viral abundance or phage-mediated bacterial mortality were more pronounced after the chlorophyll peak. Significant effects of acidification on the microbial community have been found to occur after the bloom phase (Arnosti et al., 2011; Sperling et al., 2013). During an in situ mesocosm experiment, elevated levels of pCO2 triggered higher diversity of the particle-associated bacterial community only after the breakdown of the phytoplankton bloom (Sperling et al., 2013) concomitant with increased dissolved and particulate primary production (Engel et al., 2013) and bacterial production (Piontek et al., 2013), but with decreased bacterial abundance (Brussaard et al., 2013). The latter was suggested to be due to enhanced viral lysis rates in a high pCO2 postbloom situation.
Implications for the Carbon Fluxes in the Microbial Food Web
In the experiments performed with a nano–picoplankton-dominated community, in the SUMMER oligotrophic situation, phage-mediated bacterial mortality was, in general, low without a clear effect of OA. Nonetheless, phage abundances were high. A 3-year study from the same sampling site revealed that a specific class of organic aggregates, the transparent exopolymeric particles (TEP), recurrently accumulate in the stratified, oligotrophic conditions of early summer (Ortega-Retuerta et al., 2018). Suspended organic aggregates are hotspots for bacterial growth (Azam and Long, 2001; Azam and Malfatti, 2007), and high viral abundance found on TEP (Mari et al., 2007) suggests high viral infection of TEP-attached bacteria. It must be pointed out that the virus production assays (VRAs) used in the present study exclude aggregate-attached bacteria due to prefiltration of the water and, consequently, phage production associated to TEPs. This might explain the discrepancy between the low phage production and high viral numbers observed. However, a simulation model of a Phaeocystis globosa bloom demonstrated that TEPs scavenge viral particles and reduce their infectivity (Ruardij et al., 2005). Moreover, infected host cells might be grazed selectively, lowering the contact rate between viruses and uninfected host cells (Ruardij et al., 2005). In fact, seasonal studies from the same sampling site show peaks of nanoflagellate abundance and relatively important grazing rates in summer (Boras et al., 2009; Unrein et al., 2014) suggesting a less efficient viral shunt in that season.
Instead, phage-mediated loss of bacterial biomass in WINTER was apparently the dominant mortality agent for bacteria there. Since bacterial abundance was higher in WINTER than in SUMMER, this may have increased the encounter probability with the host. In addition, it has been suggested that viruses are the prevailing mortality agent in eutrophic waters over protistan grazing (Steward et al., 1996), consistent with observations in nutrient-rich bloom/postbloom situations where the organic carbon in the microbial food web is mainly processed by the viral shunt (Malits et al., 2014). Viral lysis, by destroying the infected cell and releasing the content of the cytoplasm to the environment, mainly proteins and glucose (Boras et al., 2015), feeds the DOC pool and reduces the transfer of photosynthetically fixed carbon to higher trophic levels and to export; hence, retaining it in the basal levels of the planktonic food web (Sheik et al., 2014), increasing bacterial respiration (Middelboe and Lyck, 2002; Bonilla-Findji et al., 2008), and weakening the biological carbon pump, i.e., the flux of organic carbon to the deep ocean.
We found phage production and concomitant organic carbon release rates significantly reduced in the acidified treatments of the nutrient replete WINTER experiment, especially in the non-amended treatments. When comparing the organic carbon processed by the viral shunt with the total autotrophic biomass converted to organic carbon with a C/Chl a ratio estimated for bloom conditions in the Northwest Mediterranean (Gutiérrez-Rodríguez et al., 2010), viral lysis transferred 41 ± 5% of the phytoplankton standing stock to the DOC pool, compared to 18 ± 8% under acidification. Although carbon release rates calculated from VP might be overestimated due to the plausible stimulation of bacterial production in the VRA experiments (Weinbauer et al., 2009), it is evident that OA reduces the DOC supply for bacterial activity mediated by viral lysis (Malits and Weinbauer, 2009).
While no significant differences in the net accumulation of DOC and the protein-like DOM compounds could be detected between the acidified and non-acidified conditions in the present study (Aparicio et al., 2016), enhanced DOC production under elevated pCO2 concentrations has often been observed in published studies (Kim et al., 2011; Yoshimura et al., 2013; Paul et al., 2015). Mesocosm experiments in the Arctic Ocean have shown that DOC exudation by phytoplankton increased under acidification but without concomitant net DOC accumulation (Engel et al., 2013), suggesting that the increased bacterial biomass production (Piontek et al., 2013) was responsible for the DOC removal. Since bacterial heterotrophic production was not stimulated under OA in WINTER, reduced DOC release mediated by viral lysis under acidification could have been compensated by increased DOC production under lowered pH conditions (Kim et al., 2011; Engel et al., 2013; Yoshimura et al., 2013; Paul et al., 2015). This would have led to non-significant differences in net DOC accumulation between treatments, as we observed (Aparicio et al., 2016).
Seawater culture experiments with bacterioplankton communities from contrasting sites such as the coastal upwelling influenced Santa Barbara Channel, and the oligotrophic systems of the Sargasso Sea and the South Pacific Subtropical Gyre evidenced that elevated pCO2 led to greater DOC removal by bacteria growing on a range of added DOC types, from glucose to complex phytoplankton products (James et al., 2017). The authors inferred higher bacterioplankton respiration rates under OA from changes in DOC and bacterial biomass and suggested reduced vertical DOC export. In the present study, we indeed found that actively respiring bacterial cells labeled by CTC increased with decreasing pH, at least in the nutrient amended treatments (Table 3), in accordance with previous studies in the Mediterranean (Celussi et al., 2017). Together with the lower calcification rates of planktonic organisms and a consequently reduced particle flux to the deep ocean in an acidified environment (Engel et al., 2005), this results in a positive feedback to atmospheric pCO2. Reduced viral lysis under OA, however, should counterbalance these negative effects of lowered pH on the biological carbon pump, at least in bloom and postbloom situations where viral lysis plays a relevant role in carbon cycling through the microbial food web.
Our experimental results provide insight into how two major global change pressures, ocean acidification and eutrophication, modulate viral production and viral-mediated processes in the plankton community at different seasons. While nutrient loading consistently enhanced phage production regardless of the season, the effect of OA on viruses differed between seasons. In the nutrient-replete bloom/postbloom situation in WINTER, viral production and phage-mediated bacterial mortality were reduced under increased pCO2. In the oligotrophic SUMMER situation, where viral-mediated carbon fluxes were less important compared to WINTER, OA stimulated viruses through the cascading effect of elevated pCO2 on the autotrophic and heterotrophic production. Given the relevant role of viruses in marine carbon cycling and the biological carbon pump, these two anthropogenic stressors may modulate carbon fluxes through their effect on viruses at the base of the pelagic food web in future global change scenarios.
This is the first study that shows phage-mediated bacterial mortality leading to altered carbon fluxes in response to elevated pCO2 and nutrient additions simulating the combined effect of acidification and eutrophication. This work contributes to the understanding of how these two anthropogenic stressors modulate carbon fluxes at the base of the pelagic food web in future global change scenarios.
Data Availability Statement
The original contributions presented in the study are included in the article/Supplementary Materials, further inquiries can be directed to the corresponding authors.
AM, JB, VB, EC, CP, JP, JG, CM, MS, and DV conceived, planned, and carried out the mesocosm experiments. JB and AM performed the mortality experiments. AM analyzed viruses by flow cytometry and took the lead in writing the manuscript. All authors provided critical feedback and helped to shape the manuscript.
The financial support was provided by the Spanish Ministry of Economy and Competitiveness in the framework of the projects STORM (CTM2009–09352), DOREMI (CTM2012-34294), MANIFEST (CTM2012-32017), and PROTOS (CTM2009-08783). We are also grateful for the funding from the Generalitat de Catalunya through the Grup de Diversitat Microbiana en Ecosistemes Acuàtics (2014SGR/1591), Grup de Biogeoquímica Marina i Canvi Global (2017SGR1011) and the Grup d'Estructura i Funció de Xarxes Tròfiques Microbianes Planctòniques (2014SGR/1179) and a grant from the Spanish Ministry of Education (SB2010-0079) to AM. With the institutional support of the Severo Ochoa Center of Excellence accreditation (CEX2019-000928-S).
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
The authors thank Anderson Cabral and Clara Cardelús for research assistance, the Experimental Aquarium Zone (ZAE) of the ICM-CSIC for providing the installations to deploy the containers for the experiments, and the Analytical Chemistry laboratory staff of the ICM-CSIC for performing the nutrient analyses.
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2021.635821/full#supplementary-material
Allgaier, M., Riebesell, U., Vogt, M., Thyrhaug, R., and Grossart, H. P. (2008). Coupling of heterotrophic bacteria to phytoplankton bloom development at different pCO2 levels: a mesocosm study. Biogeosciences 5, 1007–1022. doi: 10.5194/bg-5-1007-2008
Alonso-Sáez, L., Balagué, V., Sà, E. L., Sánchez, O., González, J. M., Pinhassi, J., et al. (2007). Seasonality in bacterial diversity in north-west Mediterranean coastal waters: assessment through clone libraries, fingerprinting and FISH. FEMS Microbiol. Ecol. 60, 98–112. doi: 10.1111/j.1574-6941.2006.00276.x
Aparicio, F. L., Nieto-Cid, M., Borrull, E., Calvo, E., Pelejero, C., Sala, M. M., et al. (2016). Eutrophication and acidification: do they induce changes in the dissolved organic matter dynamics in the coastal Mediterranean Sea? Sci. Total Environ. 563–564, 179–189. doi: 10.1016/j.scitotenv.2016.04.108
Arnosti, C., Grossart, H. P., Mhling, M., Joint, I., and Passow, U. (2011). Dynamics of extracellular enzyme activities in seawater under changed atmospheric pCO2: a mesocosm investigation. Aquat. Microb. Ecol. 64, 285–298. doi: 10.3354/ame01522
Azam, F., Fenchel, T., Field, J. G., Gray, J. S., Meyer-Reil, L. A., and Thingstad, F. (1983). The ecological role of water-column microbes in the sea. Mar. Ecol. Prog. Ser. 10, 257–263. doi: 10.3354/meps010257
Baltar, F., Palovaara, J., Unrein, F., Catala, P., Hornák, K., Šimek, K., et al. (2016). Marine bacterial community structure resilience to changes in protist predation under phytoplankton bloom conditions. ISME J. 10, 568–581. doi: 10.1038/ismej.2015.135
Baltar, F., Palovaara, J., Vila-Costa, M., Salazar, G., Calvo, E., Pelejero, C., et al. (2015). Response of rare, common and abundant bacterioplankton to anthropogenic perturbations in a Mediterranean coastal site. FEMS Microbiol. Ecol. 91:fiv058. doi: 10.1093/femsec/fiv058
Berdjeb, L., Pollet, T., Domaizon, I., and Jacquet, S. (2011). Effect of grazers and viruses on bacterial community structure and production in two contrasting trophic lakes. BMC Microbiol. 11:88. doi: 10.1186/1471-2180-11-88
Bonilla-Findji, O., Malits, A., Lefèvre, D., Rochelle-Newall, E., Lemée, R., Weinbauer, M. G., et al. (2008). Viral effects on bacterial respiration, production and growth efficiency: consistent trends in the Southern Ocean and the Mediterranean Sea. Deep. Res. Part II Top. Stud. Oceanogr. 55, 559–912. doi: 10.1016/j.dsr2.2007.12.004
Bopp, L., Resplandy, L., Orr, J. C., Doney, S. C., Dunne, J. P., Gehlen, M., et al. (2013). Multiple stressors of ocean ecosystems in the 21st century: projections with CMIP5 models. Biogeosciences 10, 6225–6245. doi: 10.5194/bg-10-6225-2013
Boras, J. A., Sala, M. M., Vázquez-Domínguez, E., Weinbauer, M. G., and Vaqué, D. (2009). Annual changes of bacterial mortality due to viruses and protists in an oligotrophic coastal environment (NW Mediterranean). Env. Microbiol. 11, 1181–1193. doi: 10.1111/j.1462-2920.2008.01849.x
Boras, J. A., Vaqué, D., Maynou, F., Sà, E. L., Weinbauer, M. G., and Sala, M. M. (2015). Factors shaping bacterial phylogenetic and functional diversity in coastal waters of the NW Mediterranean Sea. Estuar. Coast. Shelf Sci. 154, 102–110. doi: 10.1016/j.ecss.2014.12.039
Borges, A. V., and Gypens, N. (2010). Carbonate chemistry in the coastal zone responds more strongly to eutrophication than ocean acidification. Limnol. Oceanogr. 55, 346–353. doi: 10.4319/lo.2010.55.1.0346
Brussaard, C. P. D. (2009). “Enumeration of bacteriophages using flow cytometry,” in Bacteriophages: Methods and Protocols, Vol. 1: Isolation, Characterization, and Interactions, eds M. R. J. Clokie and A. M. Kropinski (Totowa, NJ: Humana Press), 97–111.
Brussaard, C. P. D., Noordeloos, A. A. M., Witte, H., Collenteur, M. C. J., Schulz, K., Ludwig, A., et al. (2013). Arctic microbial community dynamics influenced by elevated CO2 levels. Biogeosciences 10, 719–731. doi: 10.5194/bg-10-719-2013
Brussaard, C. P. D., Timmermans, K. R., Uitz, J., and Veldhuis, M. J. W. (2008a). Virioplankton dynamics and virally induced phytoplankton lysis versus microzooplankton grazing southeast of the Kerguelen (Southern Ocean). Deep Sea Res. Part II Top. Stud. Oceanogr. 55, 752–765. doi: 10.1016/j.dsr2.2007.12.034
Brussaard, C. P. D., Wilhelm, S. W., Thingstad, F., Weinbauer, M. G., Bratbak, G., Heldal, M., et al. (2008b). Global-scale processes with a nanoscale drive: the role of marine viruses. ISME J. 2, 575–578. doi: 10.1038/ismej.2008.31
Bunse, C., Lundin, D., Karlsson, C. M. G., Akram, N., Vila-Costa, M., Palovaara, J., et al. (2016). Response of marine bacterioplankton pH homeostasis gene expression to elevated CO2. Nat. Clim. Chang. 6, 483–487. doi: 10.1038/nclimate2914
Cai, W.-J., Hu, X., Huang, W.-J., Murrell, M. C., Lehrter, J. C., Lohrenz, S. E., et al. (2011). Acidification of subsurface coastal waters enhanced by eutrophication. Nat. Geosci. 4:766. doi: 10.1038/ngeo1297
Celussi, M., Malfatti, F., Annalisa, F., Gazeau, F., Giannakourou, A., Pitta, P., et al. (2017). Ocean acidification effect on prokaryotic metabolism tested in two diverse trophic regimes in the Mediterranean Sea. Estuar. Coast. Shelf Sci. 186, 125–138. doi: 10.1016/j.ecss.2015.08.015
Chen, X., Ma, R., Yang, Y., Jiao, N., and Zhang, R. (2019). Viral regulation on bacterial community impacted by lysis-lysogeny switch: a microcosm experiment in eutrophic coastal waters. Front. Microbiol. 10:1763. doi: 10.3389/fmicb.2019.01763
Clayton, T. D., and Byrne, R. H. (1993). Spectrophotometric seawater pH measurements: total hydrogen ion concentration scale calibration of m-cresol purple and at-sea results. Deep Sea Res. Part I Oceanogr. Res. 40, 2115–2129. doi: 10.1016/0967-0637(93)90048-8
Crawfurd, K. J., Alvarez-Fernandez, S., Mojica, K. D. A., Riebesell, U., and Brussaard, C. P. D. (2017). Alterations in microbial community composition with increasing f CO2: a mesocosm study in the eastern Baltic Sea. Biogeosciences 14, 3831–3849. doi: 10.5194/bg-14-3831-2017
Danovaro, R., Corinaldesi, C., Dell'anno, A., Fuhrman, J. A., Middelburg, J. J., Noble, R. T., et al. (2011). Marine viruses and global climate change. FEMS Microbiol. Rev. 35, 993–1034. doi: 10.1111/j.1574-6976.2010.00258.x
Dlugokencky, E., and Tans, P. (2020). Trends in atmospheric carbon dioxide. National Oceanic & Atmospheric Administration, Earth System Research Laboratory (NOAA/ESRL). Available online at: http://www.esrl.noaa.gov/gmd/ccgg/trends/global.html
Dore, J. E., Lukas, R., Sadler, D. W., Church, M. J., and Karl, D. M. (2009). Physical and biogeochemical modulation of ocean acidification in the central North Pacific. Proc. Natl. Acad. Sci. U.S.A. 106, 12235–12240. doi: 10.1073/pnas.0906044106
Engel, A., Borchard, C., Piontek, J., Schulz, K. G., Riebesell, U., and Bellerby, R. (2013). CO2 increases 14C primary production in an Arctic plankton community. Biogeosciences 10, 1291–1308. doi: 10.5194/bg-10-1291-2013
Engel, A., Zondervan, I., Aerts, K., Beaufort, L., Benthien, A., Chou, L., et al. (2005). Testing the direct effect of CO2 concentration on a bloom of the coccolithophorid Emiliania huxleyi in mesocosm experiments. Limnol. Oceanogr. 50, 493–507. doi: 10.4319/lo.2005.50.2.0493
Field, C. B., Behrenfeld, M. J., Randerson, J. T., and Falkowski, P. (1998). Primary production of the biosphere: integrating terrestrial and oceanic components. Science 281, 237–240. doi: 10.1126/science.281.5374.237
Fu, F.-X., Warner, M. E., Zhang, Y., Feng, Y., and Hutchins, D. A. (2007). Effects of increased temperature and CO2 on photosynthesis, growth and elemental ratios in marine in marine Synechococcus and Prochlorococcus (cyanobacteria). J. Phycol. 43, 485–496. doi: 10.1111/j.1529-8817.2007.00355.x
Galloway, J. N., Dentener, F. J., Capone, D. G., Boyer, E. W., Howarth, R. W., Seitzinger, S. P., et al. (2004). Nitrogen cycles: past, present, and future. Biogeochemistry 70, 153–226. doi: 10.1007/s10533-004-0370-0
Gasol, J. M., Cardelús, C., Morán, X. A. G., Balagué, V., Forn, I., Marrasé, C., et al. (2016). Seasonal patterns in phytoplankton primary production and photosynthetic parameters in a coastal time-series station of the NW Mediterranean Sea. Scientia Marina, 80S1, 63–77. doi: 10.3989/scimar.04480.06E
Gasol, J. M., and Morán, X. A. G. (2015). “Flow cytometric determination of microbial abundances and its use to obtain indices of community structure and relative activity,” in Hydrocarbon and Lipid Microbiology Protocols, eds T. J. McGenity, K. N. Timmis, and B. Nogales (Berlin; Heidelberg: Springer Protocols Handbooks). doi: 10.1007/8623_2015_139
Gasol, J. M., Zweifel, U. L., Peters, F., Fuhrman, J. A., and Hågström, A. (1999). Significance of size and nucleic acid content heterogeneity as measured by flow cytometry in natural planktonic bacteria. Appl. Env. Microbiol. 65, 4475–4483. doi: 10.1128/AEM.65.10.4475-4483.1999
Gobler, C. J., Hutchins, D. A., Fisher, N. S., Cosper, E. M., and Sanudo-Wilhelmy, S. A. (1997). Release and bioavailability of C, N, P, Se., and Fe following viral lysis of a marine chrysophyte. Limnol. Ocean. 42, 1492–1504. doi: 10.4319/lo.19126.96.36.1992
Grossart, H. P., Allgaier, M., Passow, U., and Riebesell, U. (2006). Testing the effect of CO2 concentration on the dynamics of marine heterotrophic bacterioplankton. Limnol. Ocean. 51, 1–11. doi: 10.4319/lo.2006.51.1.0001
Gutiérrez-Rodríguez, A., Latasa, M., Estrada, M., Vidal, M., and Marrasé, C. (2010). Carbon fluxes through major phytoplankton groups during the spring bloom and post-bloom in the Northwestern Mediterranean Sea. Deep Sea Res. Part I Oceanogr. Res. Pap. 57, 486–500. doi: 10.1016/j.dsr.2009.12.013
Hansen, H. P., and Koroleff, F. (2007). “Determination of nutrients,” in Methods of Seawater Analysis, eds K. Grasshoff, K. Kremling, and M. Ehrhardt (John Wiley & Sons), 159–228. doi: 10.1002/9783527613984.ch10
James, A. K., Passow, U., Brzezinski, M. A., Parsons, R. J., Trapani, J. N., and Carlson, C. A. (2017). Elevated pCO2 enhances bacterioplankton removal of organic carbon. PLoS ONE 12:e0173145. doi: 10.1371/journal.pone.0173145
Joos, F., and Spahni, R. (2008). Rates of change in natural and anthropogenic radiative forcing over the past 20,000 years. Proc. Natl. Acad. Sci. U.S.A. 105, 1425 LP−1430. doi: 10.1073/pnas.0707386105
Kim, J.-M., Lee, K., Shin, K., Yang, E. J., Engel, A., Karl, D. M., et al. (2011). Shifts in biogenic carbon flow from particulate to dissolved forms under high carbon dioxide and warm ocean conditions. Geophys. Res. Lett. 38:L08612. doi: 10.1029/2011GL047346
Kirchman, D., K'Nees, E., and Hodson, R. (1985). Leucine incorporation and its potential as a measure of protein synthesis by bacteria in natural aquatic systems. Appl. Env. Microbiol. 49, 599–607. doi: 10.1128/AEM.49.3.599-607.1985
Knowles, B., Silveira, C. B., Bailey, B. A., Barott, K., Cantu, V. A., Cobián-Güemes, A. G., et al. (2016). Lytic to temperate switching of viral communities. Nature 531, 466–470. doi: 10.1038/nature17193
Krause, E., Wichels, A., Giménez, L., Lunau, M., Schilhabel, M. B., and Gerdts, G. (2012). Small changes in pH have direct effects on marine bacterial community composition: a microcosm approach. PLoS ONE 7:e47035. doi: 10.1371/journal.pone.0047035
Larsen, A., Castberg, T., Sandaa, R. A., Brussaard, C. P. D., Egge, J. K., Heldal, M., et al. (2001). Population dynamics and diversity of phytoplankton, bacteria and viruses in a seawater enclosure. Mar. Ecol. Prog. Ser. 221, 47–57. doi: 10.3354/meps221047
Larsen, A., Gro, A., Flaten, F., Sandaa, R. A., Castberg, T., Thyrhaug, R., et al. (2004). Spring phytoplankton bloom dynamics in Norwegian coastal waters: microbial community succession and diversity. Limnol. Oceanogr. 49, 180–190. doi: 10.4319/lo.2004.49.1.0180
Larsen, J. B., Larsen, A., Thyrhaug, R., Bratbak, G., and Sandaa, R. A. (2008). Response of marine viral populations to a nutrient induced phytoplankton bloom at different pCO2 levels. Biogeosciences 5, 523–533. doi: 10.5194/bg-5-523-2008
Lee, S., and Fuhrman, J. E. (1987). Relationship between biovolume and biomass of naturally derived marine bacterioplankton. Appl. Environ. Microbiol. 53, 1298–1303. doi: 10.1128/AEM.53.6.1298-1303.1987
Lindh, M. V., Riemann, L., Baltar, F., Romero-Oliva, C., Salomon, P. S., Granéli, E., et al. (2013). Consequences of increased temperature and acidification on bacterioplankton community composition during a mesocosm spring bloom in the Baltic Sea. Environ. Microbiol. Rep. 5, 252–262. doi: 10.1111/1758-2229.12009
Maat, D. S., Crawfurd, K. J., Timmermans, K. R., and Brussaard, C. P. D. (2014). Elevated CO2 and phosphate limitation favor Micromonas pusilla through stimulated growth and reduced viral impact. Appl. Environ. Microbiol. 80, 3119 LP−3127. doi: 10.1128/AEM.03639-13
Malits, A., Christaki, U., Obernosterer, I., and Weinbauer, M. G. G. (2014). Enhanced viral production and virus-mediated mortality of bacterioplankton in a natural iron-fertilized bloom event above the Kerguelen Plateau. Biogeosciences 11, 6841–6853. doi: 10.5194/bg-11-6841-2014
Mari, X., Kerros, M. E., and Weinbauer, M. G. (2007). Virus attachment to transparent exopolymeric particles along trophic gradients in the southwestern lagoon of New Caledonia. Appl. Env. Microbiol. 73, 5245–5252. doi: 10.1128/AEM.00762-07
Marie, D., Brussaard, C. P. D., Thyrhaug, R., Bratbak, G., and Vaulot, D. (1999). Enumeration of marine viruses in culture and natural samples by flow cytometry. Appl. Env. Microbiol. 65, 45–52. doi: 10.1128/AEM.65.1.45-52.1999
Marie, D., Partensky, F., Jacquet, S., and Vaulot, D. (1997). Enumeration and cell cycle analysis of natural populations of marine picoplankton by flow cytometry using the nucleic acid stain SYBR green I. Appl. Env. Microbiol. 63, 186–193. doi: 10.1128/AEM.63.1.186-193.1997
Marie, D., Partensky, F., Vaulot, D., and Brussaard, C. (2001). Enumeration of phytoplankton, bacteria, and viruses in marine samples. Curr. Protoc. Cytom. Chapter 11:Unit 11.11. doi: 10.1002/0471142956.cy1111s10
Motegi, C., Tanaka, T., Piontek, J., Brussaard, C. P. D., Gattuso, J. P., and Weinbauer, M. G. (2013). Effect of CO2 enrichment on bacterial metabolism in an Arctic fjord. Biogeosciences 10, 3285–3296. doi: 10.5194/bg-10-3285-2013
Murray, A. G., and Eldridge, P. M. (1994). Marine viral ecology: incorporation of bacteriophage into the microbial planktonic food web paradigm. J. Plankton Res. 16, 627–641. doi: 10.1093/plankt/16.6.627
Newbold, L. K., Oliver, A. E., Booth, T., Tiwari, B., Desantis, T., Maguire, M., et al. (2012). The response of marine picoplankton to ocean acidification. Env. Microbiol. 14, 2293–2307. doi: 10.1111/j.1462-2920.2012.02762.x
Oliver, A. E., Newbold, L. K., Whiteley, A. S., and van der Gast, C. J. (2014). Marine bacterial communities are resistant to elevated carbon dioxide levels. Environ. Microbiol. Rep. 6, 574–582. doi: 10.1111/1758-2229.12159
Ortega-Retuerta, E., Marrasé, C., Muñoz-Fernández, A., Sala, M. M., Simó, R., and Gasol, J. M. (2018). Seasonal dynamics of transparent exopolymer particles (TEP) and their drivers in the coastal NW Mediterranean Sea. Sci. Total Environ. 631–632, 180–190. doi: 10.1016/j.scitotenv.2018.02.341
Paul, A. J., Bach, L. T., Schulz, K. G., Boxhammer, T., Czerny, J., Achterberg, E. P., et al. (2015). Effect of elevated CO2 on organic matter pools and fluxes in a summer Baltic Sea plankton community. Biogeosciences. 12, 6181–6203. doi: 10.5194/bg-12-6181-2015
Paul, J. H., and Weinbauer, M. G. (2010). “Detection of lysogeny in marine environments,” in Manual of Aquatic Viral Ecology, eds C. Suttle, S. W. Wilhelm, and M. G. Weinbauer (Waco: ASLO), 1–8. doi: 10.4319/mave.2010.978-0-9845591-0-7.30
Piontek, J., Borchard, C., Sperling, M., Schulz, K. G., Riebesell, U., and Engel, A. (2013). Response of bacterioplankton activity in an Arctic fjord system to elevated pCO2: results from a mesocosm perturbation study. Biogeosciences 10, 297–314. doi: 10.5194/bg-10-297-2013
Poorvin, L., Rinta-Kanto, J. M., Hutchins, D. A., and Wilhelm, S. W. (2004). Viral release of iron and its bioavailability to marine plankton. Limnol. Oceanogr. 49, 1734–1741. doi: 10.4319/lo.2004.49.5.1734
Provoost, P., van Heuven, S., Soetaert, K., Laane, R. W. P. M., and Middelburg, J. J. (2010). Seasonal and long-term changes in pH in the Dutch coastal zone. Biogeosciences 7, 3869–3878. doi: 10.5194/bg-7-3869-2010
Raven, J. A., Giordano, M., Beardall, J., and Maberly, S. C. (2012). Algal evolution in relation to atmospheric CO2: carboxylases, carbon-concentrating mechanisms and carbon oxidation cycles. Philos. Trans. R. Soc. Lond. B. Biol. Sci. 367, 493–507. doi: 10.1098/rstb.2011.0212
Riebesell, U., Bellerby, R. G. J., Grossart, H. P., and Thingstad, F. (2008). Mesocosm CO2 perturbation studies: from organism to community level. Biogeosciences 5, 1157–1164. doi: 10.5194/bg-5-1157-2008
Rochelle-Newall, E., Delille, B., Frankignoulle, M., Gattuso, J. P., Jacquet, S., Riebesell, U., et al. (2004). Chromophoric dissolved organic matter in experimental mesocosms maintained under different pCO2 levels. Mar. Ecol. Prog. Ser. 272, 25–31. doi: 10.3354/meps272025
Roy, A. S., Gibbons, S. M., Schunck, H., Owens, S., Caporaso, J. G., Sperling, M., et al. (2013). Ocean acidification shows negligible impacts on high-latitude bacterial community structure in coastal pelagic mesocosms. Biogeosciences 10, 555–566. doi: 10.5194/bg-10-555-2013
Ruardij, P., Veldhuis, M. J. W., and Brussaard, C. P. D. (2005). Modeling the bloom dynamics of the polymorphic phytoplankter Phaeocystis globosa: impact of grazers and viruses. Harmful Algae 4, 941–963. doi: 10.1016/j.hal.2004.12.011
Sala, M. M., Aparicio, F. L., Balagué, V., Boras, J. A., Borrull, E., Cardelús, C., et al. (2016). Contrasting effects of ocean acidification on the microbial food web under different trophic conditions. ICES J. Mar. Sci. 73, 670–679. doi: 10.1093/icesjms/fsv130
Sandaa, R.-A. A., Gomez-Consarnau, L., Pinhassi, J., Riemann, L., Malits, A., Weinbauer, M. G. G., et al. (2009). Viral control of bacterial biodiversity - evidence from a nutrient-enriched marine mesocosm experiment. Env. Microbiol. 11, 2585–2597. doi: 10.1111/j.1462-2920.2009.01983.x
Schulz, K. G., Bach, L. T., Bellerby, R. G. J., Bermúdez, R., Büdenbender, J., Boxhammer, T., et al. (2017). Phytoplankton blooms at increasing levels of atmospheric carbon dioxide: experimental evidence for negative effects on prymnesiophytes and positive on small picoeukaryotes. Front. Mar. Sci. 4:64. doi: 10.3389/fmars.2017.00064
Sheik, A. R., Brussaard, C. P. D., Lavik, G., Lam, P., Musat, N., Krupke, A., et al. (2014). Responses of the coastal bacterial community to viral infection of the algae Phaeocystis globosa. ISME J. 8, 212–225. doi: 10.1038/ismej.2013.135
Sherr, B. F., del Giorgio, P. A., and Sherr, E. B. (1999). Estimating abundance and single-cell characteristics of actively respiring bacteria via the redox dye CTC. Aquat. Microb. Ecol. 18, 117–131. doi: 10.3354/ame018117
Sieracki, M. E., Cucci, T. L., and Nicinski, J. (1999). Flow cytometric analysis of 5-cyano-2,3-ditolyl tetrazolium chloride activity of marine bacterioplankton in dilution cultures. Appl. Env. Microbiol. 65, 2409–2417. doi: 10.1128/AEM.65.6.2409-2417.1999
Simek, K., Pernthaler, J., Weinbauer, M. G., Hornak, K., Dolan, J. R., Nedoma, J., et al. (2001). Changes in bacterial community composition and dynamics and viral mortality rates associated with enhanced flagellate grazing in a mesoeutrophic reservoir. Appl. Env. Microbiol. 67, 2723–2733. doi: 10.1128/AEM.67.6.2723-2733.2001
Sperling, M., Piontek, J., Gerdts, G., Wichels, A., Schunck, H., Roy, A. S., et al. (2013). Effect of elevated CO2 on the dynamics of particle-attached and free-living bacterioplankton communities in an Arctic fjord. Biogeosciences 10, 181–191. doi: 10.5194/bg-10-181-2013
Spilling, K., Paul, A. J., Virkkala, N., Hastings, T., Lischka, S., Stuhr, A., et al. (2016). Ocean acidification decreases plankton respiration: evidence from a mesocosm experiment. Biogeosciences 13, 4707–4719. doi: 10.5194/bg-13-4707-2016
Suttle, C. A. (2000). “Ecological, evolutionary, and geochemical consequences of viral infection of cyanobacteria and eukaryotic algae,” in Viral Ecology, ed C. J. Hurst (Cambridge: Academic Press), 247–296. doi: 10.1016/B978-012362675-2/50007-0
Traving, S. J., Clokie, M. R., and Middelboe, M. (2013). Increased acidification has a profound effect on the interactions between the cyanobacterium Synechococcus sp. WH7803 and its viruses. FEMS Microbiol. Ecol. 87, 133–141. doi: 10.1111/1574-6941.12199
Tsiola, A., Pitta, P., Giannakourou, A., Bourdin, G., Marro, S., Maugendre, L., et al. (2017). Ocean acidification and viral replication cycles: frequency of lytically infected and lysogenic cells during a mesocosm experiment in the NW Mediterranean Sea. Estuar. Coast. Shelf Sci. 186, 139–151. doi: 10.1016/j.ecss.2016.05.003
Vaqué, D., Lara, E., Arrieta, J. M., Holding, J., Sà, E. L., Hendriks, I. E., et al. (2019). Warming and CO2 enhance arctic heterotrophic microbial activity. Front. Microbiol. 10:494. doi: 10.3389/fmicb.2019.00494
Weinbauer, M. G., Arrieta, J. M., Griebler, C., and Herndl, G. J. (2009). Enhanced viral production and infection of bacterioplankton during an iron induced phytoplankton bloom in the Southern Ocean. Limnol. Oceanogr. 54, 774–784. doi: 10.4319/lo.2009.54.3.0774
Weinbauer, M. G., Bonilla-Findji, O., Chan, A. M., Dolan, J. R., Short, S. M., Simek, K., et al. (2011). Synechococcus growth in the ocean may depend on the lysis of heterotrophic bacteria. J. Plankton Res. 33, 1465–1476. doi: 10.1093/plankt/fbr041
Weinbauer, M. G., Brettar, I., and Höfle, M. G. (2003a). Lysogeny and virus-induced mortality of bacterioplankton in surface, deep, and anoxic marine waters. Limnol. Oceanogr. 48, 1457–1465. doi: 10.4319/lo.2003.48.4.1457
Weinbauer, M. G., Christaki, U., Nedoma, J., and Simek, K. (2003b). Comparing the effects of resource enrichment and grazing on viral production in a meso-eutrophic reservoir. Aquat. Microb. Ecol. 31, 137–144. doi: 10.3354/ame031137
Weinbauer, M. G., Rowe, J. M., and Wilhelm, S. W. (2010). “Determining rates of virus production in aquatic systems by the virus reduction approach,” in Manual of Aquatic Viral Ecology, eds C. Suttle, S. W. Wilhelm, and M. G. Weinbauer (Waco: ASLO), 1–8. doi: 10.4319/mave.2010.978-0-9845591-0-7.1
Wells, L. E., and Deming, J. W. (2006). Significance of bacterivory and viral lysis in bottom waters of Franklin Bay, Canadian Arctic, during winter. Aquat. Microb. Ecol. 43, 209–221. doi: 10.3354/ame043209
Worden, A. Z., Follows, M. J., Giovannoni, S. J., Wilken, S., Zimmerman, A. E., and Keeling, P. J. (2015). Rethinking the marine carbon cycle: factoring in the multifarious lifestyles of microbes. Science 347:1257594. doi: 10.1126/science.1257594
Yoshimura, T., Suzuki, K., Kiyosawa, H., Ono, T., Hattori, H., Kuma, K., et al. (2013). Impacts of elevated CO2 on particulate and dissolved organic matter production: microcosm experiments using iron-deficient plankton communities in open subarctic waters. J. Oceanogr. 69, 601–618. doi: 10.1007/s10872-013-0196-2
Keywords: ocean acidification, eutrophication, microbial food web, viral shunt, carbon fluxes
Citation: Malits A, Boras JA, Balagué V, Calvo E, Gasol JM, Marrasé C, Pelejero C, Pinhassi J, Sala MM and Vaqué D (2021) Viral-Mediated Microbe Mortality Modulated by Ocean Acidification and Eutrophication: Consequences for the Carbon Fluxes Through the Microbial Food Web. Front. Microbiol. 12:635821. doi: 10.3389/fmicb.2021.635821
Received: 30 November 2020; Accepted: 22 February 2021;
Published: 14 April 2021.
Edited by:Connie Lovejoy, Laval University, Canada
Reviewed by:Anne-Claire Baudoux, Centre National de la Recherche Scientifique (CNRS), France
Stéphan Jacquet, Institut National de Recherche pour l'agriculture, l'alimentation et l'environnement (INRAE), France
Copyright © 2021 Malits, Boras, Balagué, Calvo, Gasol, Marrasé, Pelejero, Pinhassi, Sala and Vaqué. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.