- Original Article
- Open Access
Viral dynamics in two trophically different areas in the Central Adriatic Sea
Helgoland Marine Researchvolume 71, Article number: 22 (2017)
To understand the activity of marine viruses, experiments on viral production, viral decay and the percentage of lytic and lysogenic bacterial cells among the total number of bacterial cells were carried out seasonally at two stations in the Adriatic Sea with different trophic conditions. Additionally, we are providing an insight on the enrichment with dissolved and particulate organic matter by viral lysis in the studied area. Viral production was higher at the coastal station than at the open-sea station. Viral decay rates were also higher at the coastal sea station than at the open-sea station, and accounted for approximately 40% of viral production at both investigated stations. The percentage of lysogenic infection was lower than that of lytical infection, which indicates the prevalence of the lytic cycle at both stations. Viruses had a significant influence on bacterial mortality through high daily removal of the bacterial standing stock at the coastal and open-sea station. The viruses contributed to the restoration of dissolved organic carbon, nitrogen and phosphorus in the microbial loop by lysing the bacterial cells at the studied stations. All the above suggest that viruses are important in the microbial food web and an important factor in the control of bacterial populations within the study area.
Viruses are the most abundant and ubiquitous component of marine microbial plankton, and the major cause of bacterial mortality, since bacteria are thought to be their main hosts [1,2,3,4,5]. Viral abundance exceeds the abundance of their bacterial host by approximately 25 times, but this value may vary, generally from 3 to 100, due to different environmental conditions [3, 6]. In nutrient-rich and more productive environments, the difference is expected to be greater in favour of viral abundance . Viruses lyse 10–50% of bacterial standing stock per day on average . This value can be even higher, depending mostly on the host density and productivity, but also on environmental conditions, such as temperature, salinity and UV radiation [3, 4, 6, 7]. Additionally, virus-mediated bacterial mortality can even match grazer-mediated mortality in different marine environments [6, 8]. The lysis of host cells can have a great effect on nutrient and biogeochemical cycling, by releasing large amounts of dissolved and particulate organic matter into seawater. Furthermore, the ecological effects of viral lysis are manifested via changes in bacterial community structure on the principle of the “Killing The Winner” Model, which has been reviewed previously [4, 6, 9, 10].
Viral activity and, therefore, bacterial mortality, can be affected by local trophic conditions [11,12,13,14], since marine environments with higher trophic levels usually support higher abundance and productivity of bacteria, which may result in an increase in viral numbers too. Different trophic conditions might also have a substantial impact on the way in which viruses replicate in the marine environment. In the lytic cycle, after infection, viruses use the host metabolism to replicate, which usually results in the release of the new viruses as a result of bursting of the host cell. In the lysogenic cycle, the viral genome integrates with the host genome and doubles along with the host. The incorporated “prophage” is dormant until induction occurs (when the host is stressed), and the lytic cycle is initiated. Several authors have indicated that lysogeny is the most favourable way of viral replication in environments with low host abundance and activity, as opposed to environments with high host abundance and activity where a lytic replication cycle prevails [13, 15,16,17]. However, new findings suggest that lysogeny could also be favoured in environments with increased host density [18, 19].
Although the distribution of viruses and their relationship to other components of the microbial community in the Adriatic is well documented [11, 20,21,22,23,24,25], the activity of viruses, has been poorly studied so far. Previous research on viral dynamics has mostly been carried out in the northern Adriatic [13, 26, 27], whereas no data from other geographical areas of the Adriatic are available. The main goal of this study is to determine the potential influence of environmental trophic conditions on viral dynamics in the Central Adriatic, since evidently the nutrient-rich and more productive environments favour higher abundance of heterotrophic bacteria and, therefore, higher abundance and activity of viruses. To achieve this goal, experiments on viral production, viral decay and the percentage of lytical and lysogenic infection of bacterial cells among the total number of bacterial cells were carried out at two stations in trophically different areas of the Central Adriatic Sea. Additionally, from the viral production experiments, we estimated the role of viruses in the release of dissolved and particulate organic matter.
Materials and methods
The experiments were performed seasonally (from October 2010 to July 2011) at in situ temperatures onboard the R/V BIOS DVA at two stations with different trophic status, situated in the Central Adriatic. The coastal sea station is located in Kaštela Bay, which is under the influence of the karstic river Jadro and is mesotrophic, whereas the open-sea station is located near the island of Vis, and is characterized as oligotrophic  (Fig. 1). Samples for the experiments, described in detail below, were collected four times with 5-L Niskin bottles at three depths (at the surface, the thermocline layer and the bottom), and were later transferred to 2-L sterile containers for further analyses. To detect the layer of the thermocline, water temperature and salinity were recorded prior to sampling, using a SeaBird 25 CTD profiler. If the thermocline was not formed, then samples were taken from the depth of 10 m and are considered as the depth of the thermocline in the results.
The abundance of marine bacteria was determined by flow cytometry as described in Marie et al. . Samples (1 mL) were stained with SybrGreen I (stock solution diluted 100×), and analysed using a Beckman Coulter EPICS XL-MCL (high flow rate from 1.0 to 1.2 μL s−1).
Bacterial cell production
Bacterial cell production was determined using the 3H-thymidine incorporation technique . Conversion factors for bacterial production were calculated from the bacterial cell number and 3H-thymidine incorporation during bacterial growth in 1 μm pre-filtered seawater : CF = (N2 − N1)/3H, where N1 and N2 represent the numbers of bacteria at the beginning and the end of the experiment, respectively, and 3H is the integrated 3H-thymidine incorporation rate during the experiment.
The abundance of marine viruses was determined as described in Noble and Fuhrman , with slight modifications. Collected samples were processed immediately without any fixative to avoid potential underestimation of virus like particles counts using formaldehyde . Samples (2 mL) were filtered through 0.02-µm filters (Anodisc; diameter: 25 mm; Al2O3, Whatman) and stained with SYBR Green I (stock solution diluted 300×). Filters were incubated in the dark for 20 min and mounted on glass slides with a drop of 50% phosphate buffer (6.7 mM, pH 7.8) and 50% glycerol, containing 0.5% ascorbic acid. Slides were stored at − 20 °C until analysis (5–10 days). Viral counts were obtained by epifluorescence microscopy (Olympus BX 51, equipped with a blue excitation filter) under 1250× magnification (objective 100×, ocular 12.5×), and are expressed as virus-like particles (vlp) per mL.
To estimate viral production, decay, and the fraction of lysogenic cells, water samples were collected from three depths (as mentioned previously) and were transferred to 2-L sterile containers for subsequent different treatments.
The virus dilution technique described by Wilhelm et al.  was used to determine viral production rates and bacterial losses due to viruses. Triplicate samples (50 mL) were transferred to sterile glass bottles and mixed with 100 mL virus-free (filtered through 0.02 µm pore sized filters) seawater and were incubated in the dark at the in situ temperature. Parallel untreated samples were also incubated as a control. Subsamples of 2 mL were taken at the beginning of the experiment, and every 3 h for a 12-h period to determine viral abundance, as previously described. Viral production rates were estimated using “VIPCAL” software from each time period with a net increase in viral abundance . Turnover rates were determined by dividing viral production rates by viral abundance. The percentage of lytically infected bacterial cells was calculated from the ratio of lysed bacterial cells and the bacterial abundance as described by Weinbauer et al. .
The burst size (BS), or the number of viruses released per cell due to viral lysis, was estimated from time-course experiments of viral production, as described in Wells and Deming , Middelboe and Lyck , and Boras et al. . In brief, an increase in viral abundance during the incubation period was divided by a decrease in bacterial abundance for the same time period as that of the viral production experiment.
Viral decay rates were obtained by monitoring the changes in viral abundance in triplicate samples that were filtered through 0.2-µm polycarbonate filters to exclude bacteria and > 0.2 µm particles . Filtered water (200 mL) was incubated in the dark at in situ temperature. Subsamples (2 mL) were taken every 3 h for 12 h and samples were immediately processed for viral counts as described previously. The decay rates were calculated from the linear regression of log-transformed data where the slope of the line is the decay constant.
Fraction of lysogenic cells
The fraction of lysogenic cells was estimated using mytomicin C (SigmaChemical Co, No. M-0503), the most effective phage inducing agent [13, 41]. Collected seawater (150 mL) was placed in incubation bottles with mytomicin C (1 µg/mL final concentration). The samples were incubated for 12 h, in parallel with untreated water samples that served as control. Subsamples (2 mL) were collected at the beginning, after 6 h and subsequently, every 3 h until the end of the experiment. Viral and bacterial counts, as well as bacterial production were determined as described previously. The results were expressed as a percentage of lysogenic cells, as described by Paul and Weinbauer . Briefly, subtraction of the viral direct counts in the experiment and viral direct counts in the control was divided by the average BS and then by bacterial abundance at the beginning of the experiment.
Virus-mediated bacterial mortality
Virus-mediated bacterial mortality (VMM) was estimated as a percentage of bacterial standing stock loss per day by dividing the lysis rate of bacteria by the number of bacteria in the original sample. The lysis rate of bacteria was calculated from the ratio between viral production and burst size as described by Luef et al. .
Release of organic carbon, nitrogen and phosphorus by viral lysis
The amount of released carbon and nitrogen was determined by multiplying the lysis rate of bacterial cells and the mean amounts of carbon and nitrogen contained in a bacterial cell. For the released carbon content, we used a value of 20 fg C cell−1 and for nitrogen, 4 fg N cell−1 [35, 42,43,44,45]. The release rate of phosphorus was calculated from the ratio of carbon and phosphorus, as explained by Cole et al. . The amounts of carbon, nitrogen and phosphorus released, are expressed as µg carbon/nitrogen/phosphorus per litre per day.
The correlations between parameters were expressed as Pearson correlation coefficients. Analysis of variance (ANOVA) and t tests were used to assess the differences in microbiological parameters throughout the water column and between investigated stations.
Bacterial abundance and production
The distribution of bacteria and bacterial production at the studied area are shown in Table 1. Bacterial abundance at the coastal station ranged from 0.40 to 0.85 × 106 cells mL−1, with a mean of 0.58 ± 0.14 × 106 cells mL−1. At the open-sea station, bacterial abundance was two-fold lower and ranged from 0.25 to 0.34 × 106 cells mL−1 with a mean of 0.30 ± 0.03 × 106 cells mL−1. A significant difference in bacterial abundances between sampled stations (t test, p < 0.05) was observed. Furthermore, a significant difference between layers and seasons was observed for the coastal- and open-sea station (p < 0.05, ANOVA), with higher bacterial abundances in the surface layer during the warmer part of the year. Bacterial production was higher at the coastal rather than at the open sea station (0.19 ± 0.12 × 104 and 0.12 ± 0.03 × 104 cells mL−1 h−1, respectively). Average bacterial turnover at the coastal station was 6.60 ± 4.11, and 7.15 ± 2.76 day−1 at the open sea station.
Viral abundance and activity
Viral abundance ranged from 10.69 to 16.54 × 106 vlp mL−1 with a mean of 13.89 ± 1.97 × 106 vlp mL−1 at the coastal sea station, and was significantly higher (t test, p < 0.05) than at the open sea station where it ranged from 6.22 to 11.67 × 106 vlp mL−1 (mean value was 7.80 ± 1.85 × 106 vlp mL−1). Viral abundance was higher at the coastal station than at the open-sea station at all three investigated depths (p < 0.05, ANOVA) (Fig. 2).
The virus-to-bacteria ratio (VBR) was similar at both stations, with a mean value of about 25 (Table 1). The VBR was negatively correlated with bacterial abundance at the costal station (r = − 0.79, n = 12, p < 0.05), and positively correlated with viral abundance at the open sea station (r = 0.87, n = 9, p < 0.05). Viral production was higher at the coastal sea station where it ranged from 0.32 to 1.87 × 106 vlp mL−1 h−1, with a mean of 0.75 ± 0.44 × 106 vlp mL−1 h−1 (Table 2). At the open-sea station, viral production ranged from 0.15 to 0.95 × 106 vlp mL−1 h−1, with a mean of 0.49 ± 0.26 × 106 vlp mL−1 h−1. Viral turnover time varied from 5.01 to 56.01 h with a mean of 25.93 ± 14.55 h at the coastal sea station, and from 6.81 to 58.62 h with a mean of 20.62 ± 12.60 h at the open-sea station. Faster turnover time was determined during the warmer part of the year at the depths of the thermocline, when viral production was also the highest (Table 2). The BS ranged from 12 to 65 viruses per bacterial cell (mean = 41 ± 17 viruses per bacterial cell) at the coastal station, and from 18 to 101 viruses per bacterial cell (mean = 62 ± 29 viruses per bacterial cell) at the open sea station, and were statistically different (t test, p < 0.05). Viral decay values were higher at the coastal sea station (mean = of 0.25 ± 0.08 × 106 vlp mL−1 h−1) than at the open sea station (mean = 0.14 ± 0.03 × 106 vlp mL−1 h−1) (Table 2; Fig. 2).
Virus-mediated mortality of bacteria, fraction of lytical and lysogenic bacterial infection
The virus-mediated mortality of bacterial standing stock (VMM) varied between 30.68 and 145.50% day−1 (mean = 78.81 ± 31.54% day−1) at the coastal sea station, and between 50.68 and 78.81% day−1 (mean = 63.75 ± 8.80% day−1) at the open sea station (Table 2). At both stations, maximum values were found in July. The percentage of lytically infected bacterial cells was slightly higher at the coastal station (mean = 25.56 ± 12.75%) then at the open sea station (mean = 22.39 ± 5.47%) (Table 2; Fig. 3). The highest mean value (48.00%) was determined in February at the coastal station (depth 10 m), and in July (33.01%) at the open sea station (bottom layer). The percentage of lysogenic bacterial infection was low, with similar mean values at both studied stations (3.20 ± 3.02% at the coastal station and 3.14 ± 1.94% at the open-sea station) (Table 2; Fig. 3).
Release of organic carbon, nitrogen and phosphorus by viral lysis
The amounts of organic carbon, nitrogen and phosphorus released by viral lysis varied significantly between stations (t test, p < 0.05), and were greater at the coastal station (mean values = 9.26 ± 4.21 µg C L−1 day−1, 1.85 ± 0.84 µg N L−1 day−1, 0.37 ± 0.17 µg P L1 day−1) than at the open sea station (mean values = 3.80 ± 0.60 µg C L−1 day−1, 0.76 ± 0.12 µg N L−1 day−1, 0.15 ± 0.02 µg P L−1 day−1) (Table 3).
Since the investigated sites are located in the central Adriatic area, which is characterized by the different trophic conditions of the coastal and open sea waters [24, 25, 47,53,49], the differences in viral and bacterial abundance and activity are expected. Bacterial abundance was lower at the open than at the coastal sea station, which is consistent with the results previously found for the North [11, 13] and Mid Adriatic [25, 28]. Bacterial production followed a similar pattern to that of abundance, which is in accordance with previously published data for the same area , for the North Adriatic  and the Western Mediterranean Sea . The abundance of marine viruses was significantly lower (t test, p < 0.05) at the open-sea station than at the costal station, as previously described by Ordulj et al. . A similar pattern was observed in the Northern Adriatic [13, 21, 22, 50]. It is known that viral and bacterial abundances change according to the productivity of the aquatic environment . A higher abundance of marine bacteria is usually followed by a higher abundance of marine viruses, and thus, viral production .
The virus-to-bacteria ratio (VBR) was similar at both stations. VBR values are within the range previously determined for the central Adriatic . The similarities in VBR values between the two investigated stations could be explained by the higher production of viruses and their higher decay at the coastal station compared to the open sea station. Parrika et al.  noted that in the open and offshore waters, in less productive waters generally, VBR values are higher than in the usual eutrophic coastal waters, mostly because of vivory (grazing by the HNF), adsorption by particulate matter, degradation by heat-labile organic matter, UV radiation and temperature. Since the coastal station in this study is not eutrophic, but rather mesotrophic , similarities in VBR are expected. Lower than average VBR values were observed at the coastal-sea station during February at the bottom, and July in the surface layer, which might have been caused by low temperature and high UV radiation, respectively. Lower than average VBR values were also observed at the open-sea station during the warmer part of the year, which could be the result of higher UV radiation exposure of the surface layer. It is known that marine viruses are susceptible to UV radiation, especially in the surface layer [40, 52, 53]. Additionally, low temperature can affect virus-host interactions, which can result in less successful viral proliferation, and thus low abundances . Furthermore, since viral abundance changes over the course of days and weeks, it is possible that temporal changes could substantially affect the VBR .
The virus dilution technique  has become the preferred method for determining the activity of marine viruses . A dilution of viral density ensures the low occurrence of new phage infections and reduces possible viral losses due to new infection, predation and enzymatic degradation [34, 36]. Viral production showed the same pattern as viral abundance (Table 2; Fig. 2), and was higher at the coastal sea station. The highest values were observed at the depth of the thermocline during the warmer part of the year (June, July), at both investigated stations.
Viral production at the coastal sea station was within the range of data previously reported for the Gulf of Trieste , but was three-fold lower than the values estimated for the eutrophic area close to the river Po . These results are not surprising, since the estuarine area of the river Po has a higher trophic level than Kaštela Bay, due to higher nutrient input from the river Po . At the open sea station, viral production was similar to that found for the north-western Mediterranean , which was previously characterized as an oligotrophic area with low viral abundance .
At the studied stations, average bacterial turnover was 7 day−1, while average viral turnover was 1 day−1. In general, a bacterial community with lower density and lower bacterial turnover can sustain specific viral production, turnover and thus relatively high viral abundance because it can be more valuable to viruses, with a potential longer latent period, which results in more progeny per infected bacterial cell . This is in accordance with our findings, i.e. higher BS at the open sea station (41 ± 17 for coastal and 62 ± 29 on open sea station). On the other hand, lower viral decay at the open sea station could also be one of the reasons for the higher virus abundance and equal average VBR value determined at both stations. We also need to mention that although the dilution technique applied for the estimation of viral production currently is the most reliable method, it is based on sample processing that can potentially lead to inaccurate estimates . Additionally, while measuring viral production we did not consider that cyanobacterial phages could be present at both stations since the abundance of cyanobacteria there varies between 103 and 104 ml−1  and, therefore, cyanobacterial phages could affect the viral production estimates.
BS was slightly higher at the open-sea station, although BS usually decreases from eutrophic to oligotrophic areas [11, 13]. Since a VBR value of 25 was determined for the open-sea area, relatively higher BS values allowed for a greater number of free viral particles in the water column. With higher BS values, average turnover time was shorter at the open-sea station than at the coastal station (20.62 ± 12.60–25.93 ± 14.55 h, respectively). It should be noted that turnover time in the northern Adriatic  was even shorter, which corresponds with the higher trophic state of that area, compared to the area studied here. BS was positively correlated with viral production at the open sea station (r = 0.95, n = 9, p < 0.05). The maximum values for both investigated stations were determined at the depth of thermocline (Table 2). This is significant because the thermocline layer represents a barrier where nutrients are retained , making it more productive thus favouring higher viral production, and therefore higher burst size. Short turnover time and higher BS at the coastal and open-sea stations were followed by the highest values of viral production.
The results of viral decay were similar to the results determined for the mesotrophic and oligotrophic area of the north Adriatic basin . Because viruses decayed more slowly than they were produced, viral abundance in the water column was positively affected. At the coastal and open-sea stations, viral decay accounted for 43.15 ± 25.19 and 36.34 ± 21.91% of viral production, respectively, on average. The result for the open-sea station in this study is slightly different from the data published by Bongiorni et al. , where viral decay only partially balanced viral production (1.6%) at the oligotrophic open-sea station. The observed difference can be explained by two- to three-fold greater abundance of viruses in the water column, lower viral production, and higher production of bacteria at every investigated depth of the open-sea station in the middle Adriatic. Additionally, the vicinity of the mainland to the open-sea station can lead to a higher content of heat-labile and colloidal particles that can cause higher removal of viruses [40, 61]. By using the filtration based approach for the estimation of viral decay, we might have partially underestimated viral decay rates since bacteria, suspended matter and other organic molecules are removed during the process, and other physical factors such as UV radiation, were not considered either. However, we did observe different rates at the coastal and open sea stations, which were characterized by different trophic conditions ; therefore, our data can be considered, although they need to be confirmed by other approaches and methods of analysis.
VMM values were relatively high, which is in accordance with previously published data for the Adriatic and Mediterranean Sea [11, 16, 20, 26, 39]. The highest average value of VMM was observed in the thermocline layer for both investigated stations (91.20 ± 37.90% at the coastal and 67.81 ± 6.94% at the open-sea station). These results indicate that viruses are active components of microbial ecosystems and play a significant role in shaping bacterial communities. Several reports have stated that VMM can match or even exceed the mortality of bacteria caused by bacterivores, primarily heterotrophic nanoflagellates (HNF), which are the main predators of bacteria [4, 8, 9, 62]. According to our data, the removal rate of bacterial cells by viruses at the coastal-sea station ranged from 0.64 to 3.43 × 104 cells mL−1 h−1, with a mean value of 1.96 ± 0.92 × 104 cells mL−1 h−1, which is two-fold lower than the mean rate of bacterial removal by HNF, as previously calculated by Šolić and Krstulović [63, 64]. However, this research shows that viruses can remove a substantial part of the bacterial standing stock and that viruses, together with HNF grazing, generally play an important role in controlling the bacterial population. It is important to mention that in 1993, when grazing experiments were carried out by Šolić and Krstulović [63, 64], Kaštela Bay was characterised by a higher trophic level than when our research was conducted . This suggests that HNF grazing levels could be lower today.
The highest value of lysogenic infection was determined at the coastal station in February, in the bottom layer (11.78%), during low temperatures, although there was no statistically significant correlation between lysogeny and temperature. The highest value of lysogenic infection at the open-sea station was determined in July, at the surface (6.28%), which could have been caused by higher exposure to UV radiation during sampling. These results suggest that viruses predominantly caused lytic rather than lysogenic infection, which is confirmed by the high percentage of lytically infected bacteria at the coastal and open-sea station (Fig. 3), as mentioned before. The lytic cycle prevails over the lysogenic cycle in marine ecosystems probably because the production of viruses during the lytic cycle is dependant on the rate of multiplication of the viral genome in the host cell, whereas the lysogenic cycle depends largely on the growth rate of the host . Lysogeny is often the preferred strategy of viral replication in areas characterized by low bacterial abundance and production . Various authors have suggested that in order to maintain the lytic cycle in the marine ecosystems, viral density should be at least 104 viruses per mL [15, 48, 58, 66], and these conditions were found at the investigated stations. The percentage of lytically infected cells varies according to the trophic status of an area [4, 67]; therefore, higher values at a coastal sea station are expected.
The results for the amounts of organic carbon, nitrogen and phosphorus released by viral lysis fall within the range of those found in the northern Adriatic [13, 26], and are higher than the results reported for the oligotrophic Western Mediterranean [16, 39]. By lyzing a large proportion of bacterial cells per day, viruses have a substantial impact on nutrient cycling in the marine environment . While bacterial grazing by HNF and ciliates transfers particulate organic matter to higher trophic levels, viral lysis transforms bacterial cells into a pool of particulate and dissolved organic matter, which then becomes available to bacterial cells and other microorganisms. Viruses could be a significant source of released dissolved organic carbon, nitrogen and phosphorus, particularly at the coastal station where average release rates were higher (Table 3). Furthermore, these results indicate the importance of viruses in nutrient cycling, especially of phosphorus, in P-limited  and oligotrophic areas such as the Adriatic Sea [69, 70].
This study provides insights into the viral dynamics in the waters of the central Adriatic. The viral production and decay rates were higher in the coastal waters due to the higher trophic condition than in the open waters. The viruses represent an extremely active component of the microbial food network and notably influence the bacterial community, since they are a significant cause of bacterial mortality in the studied area. These results suggest that viruses predominantly cause lytic rather than lysogenic infections in the central Adriatic. Ultimately, this study shows that viral activity is an important source of organic nutrients, especially phosphorus, in the P-limited Adriatic Sea.
dissolved organic carbon
dissolved organic nitrogen
dissolved organic phosphorus
Thingstad TF, Heldal M, Bratbak G, Dundas I. Are viruses important partners in pelagic food webs? TREE. 1993;8(6):209–13.
Bratbak G, Thingstad F, Heldal M. Viruses and the microbial loop. Microb Ecol. 1994;28:209–21.
Fuhrman JA. Marine viruses and their biogeochemical and ecological effects. Nature. 1999;399:541–8.
Weinbauer MG. Ecology of prokaryotic viruses. FEMS Microbiol Rev. 2004;28:127–81.
Suttle CA. Viruses in the sea. Nature. 2005;437:356–61.
Jacquet S, Miki T, Noble RT, Peduzzi P, Wilhelm S. Viruses in aquatic ecosystems: important advancements of the last 20 years and prospects for the future in the field of microbial oceanography and limnology. Adv Ocean Limnol. 2010;1:71–101.
Sime-Ngando T. Environmental bacteriophages: viruses of microbes in aquatic ecosystems. Front Microbiol. 2014;5:355.
Fuhrman JA, Noble RT. Viruses and protists cause similar bacterial mortality in costal seawater. Limnol Oceanogr. 1995;40:1236–42.
Wommack KE, Colwell RR. Virioplankton: viruses in aquatic ecosystems. Microbiol Mol Biol Rev. 2000;64:2180–5.
Suttle CA. Marine viruses—major players in the global ecosystem. Nat Rev Microbiol. 2007;5:801–12.
Weinbauer MG, Fuks D, Puskarić S, Peduzzi P. Distribution of viruses and dissolved DNA along a coastal trophic gradient in the Northern Adriatic Sea. Appl Environ Microbiol. 1993;59:4074–82.
Noble RT, Fuhrman JA. Rapid virus production and removal as measured with fluorescently labeled viruses as tracers. Appl Environ Microbiol. 2000;66(9):3790–7.
Bongiorni L, Magagnini M, Armeni M, Noble RT, Danovaro R. Viral production, decay rates, and life strategies along a trophic gradient in the north Adriatic. Appl Environ Microbiol. 2005;71:6644–50.
Malits A, Christaki U, Obernosterer I, Weinbauer MG. Enhanced viral production and virus-mediated mortality of bacterioplankton in a natural iron-fertilized bloom event above the Kerguelen Plateau. Biogeosciences. 2014;11:6841–53.
Weinbauer MG, Suttle CA. Potential significance of lysogeny to bacteriophage production and bacterial mortality in coastal waters of the Gulf of Mexico. Aquat Microb Ecol. 1996;62:4374–80.
Weinbauer MG, Brettar I, Hofle MG. Lysogeny and virus-induced mortality of bacterioplankton in surface, deep, and anoxic marine waters. Limnol Oceanogr. 2003;48(4):1457–65.
Paul JH, Weinbauer MG. Detection of lysogeny in marine environments. MAVE. 2010;4:30–3.
Knowles B, Silveira CB, Bailey BA, Barott K, Cantu VA, Cobián-Güemes AG, Coutinho FH, Dinsdale EA, Felts B, Furby KA, George EE. Lytic to temperate switching of viral communities. Nature. 2016;531(7595):466–70.
Wigington CH, Sonderegger D, Brussard CPD, Buchan A, Finke JF, Fuhrman JA, Lennon JT, Middelboe M, Suttle CA, Stock C, Wilson WH, Wommack KE, Wilhelm SW, Weitz JS. Re-examination of the relationship between marine virus and microbial cell abundances. Nat Microbiol. 2016;1:15024.
Weinbauer MG, Fuks D, Puskaric S, Peduzzi P. Diel, seasonal and depth related variability of viruses and dissolved DNA in the Northern Adriatic Sea. Microb Ecol. 1995;30:25–41.
Weinbauer MG, Peduzzi P. Frequency, size and distribution of bacteriophages in different marine bacterial morphothypes. Mar Ecol Prog Ser. 1994;108:11–20.
Weinbauer MG, Peduzzi P. Significance of viruses versus heterotrophic nanoflagellates for controlling bacterial abundance in the northern Adriatic Sea. J Plankton Res. 1995;17:1851–6.
Corinaldesi C, Crevatin E, Del Negro PF, Marini M, Russo A, Fonda-Umani S, Danovaro R. Large-scale spatial distribution of virioplankton in the Adriatic Sea: testing the trophic state control hypothesis. Appl Environ Microbiol. 2003;69:2664–73.
Krstulović N, Šolić M, Šantić D, Maršić-Lučić J, Ordulj M, Šestanović S. Microbial community structure in two anchialine caves on Mljet Island (Adriatic Sea). Acta Adriat. 2013;54(2):183–98.
Ordulj M, Krstulović N, Šantić D, Jozić S, Šolić M. Distribution of marine viruses in the Central and South Adriatic Sea. Medit Mar Sci. 2015;16(1):65–72.
Karuza A, Del Negro P, Crevatin E, Fonda Umani S. Viral production in the Gulf of Trieste (Northern Adriatic Sea): preliminary results using different methodological approaches. J Exp Mar Biol Ecol. 2010;383:96–104.
Karuza A, Umani SF, Del Negro P. The (un) coupling between viruses and prokaryotes in the Gulf of Trieste. Estuar Coast Shelf Sci. 2012;115:87–97.
Šantić D, Krstulović N, Šolić M, Ordulj M, Kušpilić G. Dynamics of prokaryotic picoplankton community in the central and southern Adriatic Sea (Croatia). Helgol Mar Res. 2012;67(3):471.
Marie D, Partensky F, Jacquet S, Vaulot D. Enumeration and cell cycle analysis of natural populations of marine picoplankton by flow cytometry using the nucleic acid stain SYBR green I. Appl Environ Microbiol. 1997;63:186–93.
Fuhrman JA, Azam F. Bacterioplankton secondary production estimates for coastal waters of British Columbia, Antarctica, and California. Appl Environ Microbiol. 1980;39:1085–95.
Riemann B, Bjmnsen PK, Newell S, Fallon R. Calculation of cell production of coastal marine bacteria based on measured incorporation of [3H]-thymidine. Limnol Oceanogr. 1987;32:471–6.
Noble RT, Fuhrman JA. Use of SYBR Green I for rapid epifluorescence counts of marine viruses and bacteria. Aquat Microb Ecol. 1998;14:113–8.
Wen K, Ortmann AC, Suttle CA. Accurate estimation of viral abundance by epifluorescence microscopy. Appl Environ Microbiol. 2004;70(7):3862–7.
Wilhelm SW, Brigden SM, Suttle CA. A dilution technique for the direct measurement of viral production: a comparison in stratified and tidally mixed coastal waters. Microb Ecol. 2002;43:168–73.
Luef B, Luef F, Peduzzi P. Online program “VIPCAL” for calculating lytic viral production and lysogenic cell based on viral reduction approach. Environ Microbiol Rep. 2009;1:78–85.
Weinbauer MG, Rowe JM, Wilhelm SW. Determining rates of virus production in aquatic systems by the virus reduction approach. MAVE. 2010;1:1–8.
Wells LE, Deming JW. Significance of bacterivory and viral lysis in bottom waters of Franklin Bay, Canadian Arctic, during winter. Aquat Microb Ecol. 2006;43(3):209–21.
Middelboe M, Lyck PG. Regeneration of dissolved organic matter by viral lysis in marine microbial communities. Aquat Microb Ecol. 2002;27(2):187–94.
Boras A, Montserrat Sala M, Vasquez-Dominugues E, Weinbauer MG, Vaque D. Annual changes of bacterial mortality due to viruses and protest in an oligotrophic coastal environment (NW Mediterranean). Environ Microbiol. 2009;11(5):1181–93.
Noble RT, Fuhrman JA. Virus decay and it’s causes in coastal waters. Appl Environ Microbiol. 1997;63:77–83.
Ortmann AC, Lawrence JE, Suttle CA. Lysogeny and lytic viral production during a bloom of the cyanobacterium Synechococcus spp. Microb Ecol. 2002;43(2):225–31.
Lee S, Fuhrman JA. Relationships between biovolume and biomass of naturally derived marine bacterioplankton. Appl Environ Microbiol. 1987;53:1298–303.
Simon M, Azam F. Protein content and protein synthesis rates of planktonic bacteria. Mar Ecol Prog Ser. 1989;51:201–13.
Fukuda R, Ogawa H, Nagata T, Koike I. Direct determination of carbon and nitrogen contents of natural bacterial assemblage in marine environments. Appl Environ Microbiol. 1998;64:3352–8.
Gundersen K, Heldal M, Norland S, Purdie D, Knap A, Elemental C. N, and P cell content of individual bacteria collected at the Bermuda Atlantic time-series study (BATS) site. Limnol Oceanogr. 2002;47:525–1530.
Cole JJ, Pace ML, Caraco NF, Steinhart GS. Bacterial biomass and cell size distributions in lakes: more and larger cells in anoxic waters. Limnol Oceanogr. 1993;38:1627–32.
Krstulović N, Šolić M, Marasović I. Relationship between bacteria, phytoplankton and heterotrophic nanoflagellates along the trophic gradient. Helgöland Meeresuntersuch. 1997;51:433–43.
Šolić M, Krstulović N, Vilibić I, Bojanić N, Kušpilić G, Šestanović S, Šantić D, Ordulj M. Variability in the bottom-up and top-down control of bacteria on trophic and temporal scale in the middle Adriatic Sea. Aquat Microb Ecol. 2009;58:15–29.
Šantić D, Šestanović S, Šolić M, Krstulović N, Kušpilić G, Ordulj M, Ninčević Gladan Ž. Dynamics of the picoplankton community from coastal waters to the open sea in the central Adriatic. Medit Mar Sci. 2014;15(1):179–88.
Stopar D, Cerne A, Zeigman M, Poljsak-Prijatelj M, Turk V. Viral abundance and a high proportion of lysogens suggest that viruses are important members of the microbial community in the Gulf of Trieste. Microb Ecol. 2004;47:1–8.
Parikka KJ, Le Romancer M, Wauters N, Jacquet S. Deciphering the virus-to-prokaryote ratio (VPR): insights into virus–host relationships in a variety of ecosystems. Biol Rev. 2017;92(2):1081–100.
Jacquet S, Bratbak G. Effects of ultraviolet radiation on marine virus-phytoplankton interactions. FEMS Microbiol Ecol. 2003;44(3):279–89.
Wilhelm SW, Jeffrey WH, Dean AL, Meador J, Pakulski JD, Mitchell DL. UV radiation induced DNA damage in marine viruses along a latitudinal gradient in the southeastern Pacific Ocean. Aquat Microb Ecol. 2003;31(1):1–8.
Mojica KD, Brussaard CP. Factors affecting virus dynamics and microbial host–virus interactions in marine environments. FEMS Microbiol Ecol. 2014;89(3):495–515.
Siokou-Frangou I, Christaki U, Mazzocchi MG, Montresor M, Ribera d’Alcalá M, Vaqué D, Zingone A. Plankton in the open Mediterranean Sea: a review. Biogeosciences. 2010;7(5):1543–86.
Degobbis D, Precali R, Ivančić I, Smodlaka N, Fuks D, Kveder S. Long-term changes in the northern Adriatic ecosystem related to anthropogenic eutrophication. Int J Environ Pollut. 2000;13:495–533.
Alonso MC, Jimenez-Gomez F, Rodriguez J, Borrego JJ. Distribution of virus-like particles in an oligotrophic marine environment (Alboran Sea, Western Mediterranean). Microb Ecol. 2001;42:407–15.
Abedon ST, Herschler TD, Stopar D. Bacteriophage latent-period evolution as a response to resource availability. Appl Environ Microbiol. 2001;67:4233–41.
Šantić D, Krstulović N, Šolić M, Kušpilić G. Distribution of Synechococcus and Prochlorococcus in the central Adriatic Sea. Acta Adriat. 2011;52(1):101–14.
George T, George P, Costas D, Theodorou A. Assessing marine ecosystem response to nutrient inputs. Mar Pollut Bull. 2001;43(7):175–86.
Wilhelm SW, Weinbauer MG, Suttle CA, Jeffrey WH. The role of sunlight in the removal and repair of viruses in the sea. Limnol Oceanogr. 1998;43(4):586–92.
Guixa-Boixereu N, Vaqué D, Gasol JM, Sánchez-Cámara J, Pedrós-Alió C. Viral distribution and activity in Antarctic waters. Deep Sea Res Pt II. 2002;49(4):827–45.
Šolić M, Krstulović N. Role of predation in controlling bacterial and heterotrophic nanoflagellate standing stocks in the coastal Adriatic Sea: seasonal patterns. Mar Ecol Prog Ser. 1994;114:219–35.
Šolić M, Krstulović N. Bacterial carbon flux through the microbial loop in Kaštela Bay (Adriatic Sea). Ophelia. 1995;41:345–60.
Šolić M, Krstulović N, Kušpilić G, Ninčević Gladan Ž, Bojanić N, Šestanović S, Šantić D, Ordulj M. Changes in microbial food web structure in response to changed environmental trophic status: a case study of the Vranjic Basin (Adriatic Sea). Mar Environ Res. 2010;70:239–49.
Wilcox RM, Fuhrman JA. Bacterial viruses in coastal seawater: lytic rather than lysogenic production. Mar Ecol Prog Ser. 1994;114:35–45.
Ackermann HW, DuBow MS. Viruses of Prokaryotes. Vol. I. General properties of bacteriophages. Boca Raton: CRC Press Inc; 1987.
Riemann L, Holmfeldt K, Titelman J. Importance of viral lysis and dissolved DNA for bacterioplankton activity in a P-limited estuary, northern Baltic Sea. Microb Ecol. 2009;57:286–94.
Vukadin I, Stojanoski L. Phosphorus versus nitrogen limitation in the middle Adriatic Sea. Rapport du 36e Congrès de la Commission Internationale pour l’Exploration Scientifique de la mer Mediterranee. 2001;36:174.
Šolić M, Krstulović N, Šantić D, Šestanović S, Ordulj M, Bojanić N, Kušpilić G. Structure of microbial communities in phosphorus-limited estuaries along the eastern Adriatic coast. J Mar Biol Assoc UK. 2015;95:1565–78.
MO and NK designed this study. MO conducted the experiments, sampling, and processed the samples. MO wrote the paper with great help of MŠ and SJ and other authors. DŠ performed the bacterial counts on Flow cytometer. All authors read and approved the manuscript.
We would like to thank prof. Peter Peduzzi for his valuable advice and help with the manuscript. We would also like to express our gratitude to Mr. Mate Pavlović and the crew of R/V Bios Dva for their assistance. This research was supported by the Croatian Science Foundation as part of the research project: IP-2014-09-4143 “Marine microbial food web processes in global warming perspective” (MICROGLOB).
The authors declare that they have no competing interests.
Availability of data and materials
The datasets generated and analysed during the current study are not publicly available due to copyright restrictions, but may be available from the corresponding author on reasonable request.
Consent for publication
All authors have read and approved the manuscript.
Ethics approval and consent to participate
The research involved no human subjects, human material or human data. Sampling was done on-board the BIOS Dva research vessel.
This research was supported by the Croatian Science Foundation as part of the research project: IP-2014-09-4143 “Marine microbial food web processes in global warming perspective” (MICROGLOB).
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.