Next Article in Journal
Molecular Correlation between Larval, Deutonymph and Adult Stages of the Water Mite Arrenurus (Micruracarus) Novus
Previous Article in Journal
Pharmacokinetic Interactions between Herbal Medicines and Drugs: Their Mechanisms and Clinical Relevance
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Assessing Viral Abundance and Community Composition in Four Contrasting Regions of the Southern Ocean

1
Institut de Ciències del Mar-Consejo Superior de Investigaciones Científicas (ICM-CSIC), Passeig Marítim de la Barceloneta, 37-49, 08003 Barcelona, Catalonia, Spain
2
Instituto Español de Oceanografía (IEO), Centro Oceanográfico de Málaga, 29640 Fuengirola, Spain
3
C/Condi de la Cañada s/n, 10857 Acebo (Cáceres), Spain
4
Laboratoire d’Océanographie Microbienne, (CNRS UMR7621), 66650 Banyuls sur Mer, France
5
Red Sea Research Center (RSRC), King Abdullah University of Science and Technology (KAUST), Thuwal 23955, Saudi Arabia
6
Institute of Oceanography and Global Change (IOCAG), University of Las Palmas de Gran Canaria (ULPGC), 35214 Telde, Spain
*
Authors to whom correspondence should be addressed.
Life 2020, 10(7), 107; https://doi.org/10.3390/life10070107
Submission received: 4 May 2020 / Revised: 1 July 2020 / Accepted: 3 July 2020 / Published: 5 July 2020
(This article belongs to the Special Issue Microbial Community Structure and Metabolic Networks in Polar Areas)

Abstract

:
We explored how changes of viral abundance and community composition among four contrasting regions in the Southern Ocean relied on physicochemical and microbiological traits. During January–February 2015, we visited areas north and south of the South Orkney Islands (NSO and SSO) characterized by low temperature and salinity and high inorganic nutrient concentration, north of South Georgia Island (NSG) and west of Anvers Island (WA), which have relatively higher temperatures and lower inorganic nutrient concentrations. Surface viral abundance (VA) was highest in NSG (21.50 ± 10.70 × 106 viruses mL−1) and lowest in SSO (2.96 ± 1.48 × 106 viruses mL−1). VA was positively correlated with temperature, prokaryote abundance and prokaryotic heterotrophic production, chlorophyll a, diatoms, haptophytes, fluorescent organic matter, and isoprene concentration, and was negatively correlated with inorganic nutrients (NO3−, SiO42−, PO43−), and dimethyl sulfide (DMS) concentrations. Viral communities determined by randomly amplified polymorphic DNA–polymerase chain reaction (RAPD-PCR) were grouped according to the sampling location, being more similar within them than among regions. The first two axes of a canonical correspondence analysis, including physicochemical (temperature, salinity, inorganic nutrients—NO3−, SiO42−, and dimethyl sulfoniopropionate -DMSP- and isoprene concentrations) and microbiological (chlorophyll a, haptophytes and diatom, and prokaryote abundance and prokaryotic heterotrophic production) factors accounted for 62.9% of the variance. The first axis, temperature-related, accounted for 33.8%; the second one, salinity-related, accounted for 29.1%. Thus, different environmental situations likely select different hosts for viruses, leading to distinct viral communities.

1. Introduction

Marine viruses, the smallest and most abundant biological entities in the sea, play key roles in biogeochemical cycles, shape microbial communities, and are the largest reservoir of diversity throughout the water column from the tropics to polar systems [1,2]. The Southern Ocean (SO), which surrounds Antarctica, plays an essential role in regulating the world’s climate by contributing to the global water circulation system. One of its main traits is the Antarctic Circumpolar Current (ACC), which flows clockwise around Antarctica. The SO is characterized by low temperature and relatively low salinity waters and comprises open sea regions in which high inorganic nutrient availability coexists with low chlorophyll concentrations (HNLC regions), primarily due to a lack of iron supply for phytoplankton to perform photosynthesis [3] and the limited light during most of year. The phenomenon of lower algal biomass than that expected for the high concentrations of inorganic nutrients present in the system is known as the Antarctic paradox [4]. However, during the Austral summer, notable phytoplankton blooms occur throughout the area [5,6]. These are followed by the proliferation of prokaryotes, heterotrophic protists and viruses, which can reach high levels of abundance, activity, and diversity [7]. Processes like grazing by protists and/or zooplankton (through sloppy feeding), as well as viral lysis of prokaryotes and phytoplankton, break microbial cells [8] and promote the leaching of organic matter and micronutrients, such as iron-rich organic compounds, which become available for the growth of prokaryotic and eukaryotic phytoplankton [9,10,11,12].
Furthermore, viruses, phytoplankton, and prokaryotes are presumably key agents involved in making secondary metabolites, including volatiles that escape to the atmosphere and eventually may evolve into marine secondary aerosols, crucial in the creation of cloud condensation nuclei and, therefore, having consequent effects on climate [13,14]. A conspicuous secondary metabolite is dimethyl sulfoniopropionate (DMSP), an algal osmolyte that is produced in high intracellular concentrations by many phytoplankton taxa [15,16,17]. Two of the major aerosol-forming volatiles are isoprene (2-methyl-1,3-butadiene), which is a by-product of algal photosynthesis [14,18], and dimethyl sulfide (DMS), which derives from DMSP through the action of enzymatic lyase activity [15,16,19,20]. Generally, haptophytes, cryptophytic, and small dinoflagellates are the groups recognized to be greater DMSP producers [16,21]. DMSP is released from the cell mainly through senescence or exudation in late phases of blooms [22,23], but most importantly through grazing [20,24] and viral attack [25,26,27] and is partly transformed into DMS. Isoprene is more related to photosynthetic activity as its concentration peaks match those of phytoplankton activity [28,29,30,31], and chlorophyll-normalized isoprene production rates are known for their different phytoplankton taxa (see Booge [32] for a review). However, the relationship between phytoplankton, viruses, and isoprene remains poorly explored.
The dominating groups of phytoplankton during austral summer are diatoms, dinoflagellates, haptophytes, cryptophytes, and prasinophytes [33,34,35,36]. Phytoplankton proliferations are followed by specific prokaryotic taxa [37] that covary with changes in the phytoplankton community [38,39,40]. Both prokaryotes and phytoplankton are subjected to grazing by protists and lysis by viruses [41,42,43,44]. Environmental fluctuations could influence the abundance and activity of potential hosts, which will imply changes in the abundance and composition of the viral community. This was observed in several Antarctic lakes ([45] and refs therein), where changes from a single-stranded DNA to a double-stranded DNA–virus-dominated assemblage appeared to respond to a seasonal shift in host organisms [46]. At the same time, virus mortality processes are constrained by environmental factors [47,48]. For example, it has been reported that virus-host interactions are affected by changes in temperature and salinity [49,50], as well as by inorganic and organic nutrient concentrations, as viruses have a high demand for nitrogen and phosphorous during replication [51]. However, for Antarctic marine waters, there is still little information known on the sources of viral abundance and community variation, except for the study carried out by Miranda [52] who described seasonal changes of ssRNA viral communities in coastal Antarctic waters, in relation to variations of phytoplankton communities in spring and summer.
In this study, our goal was to explore the relationships between physicochemical and biological factors, and the abundance and composition of marine viral communities from the surface and the DCM (deep chlorophyll maximum). We studied these relationships in four regions of the Southern Ocean, characterized by distinct hydrographic conditions and dominated by different phytoplankton groups. Specifically, we tested whether (i) viral abundance and community composition were different among the four regions, and (ii) to which degree these dissimilarities were linked to the variability of potential hosts or environmental factors. To achieve our objective, we assessed the abundance and viral community composition in these four regions, and examined their potential relation to an array of (1) physicochemical (temperature, salinity, inorganic nutrients, fluorescent dissolved organic matter -FDOM- as an indicator of organic matter, and DMSP, DMS, and isoprene concentrations), and (2) biological (prokaryote abundance and prokaryotic heterotrophic production, chlorophyll a concentration, and phytoplankton taxonomic composition and biomass) variables. Additionally, we discuss the potential role of the distinct viral communities detected at each site in the production of aerosol precursor compounds.

2. Materials and Methods

2.1. Sampling Regions and Strategy

The PEGASO cruise was carried out between January and early February 2015 on board of the R/V BIO-Hespérides in the Southern Ocean near the Antarctic Peninsula, and the Sub-Antarctic Ocean (Figure 1A). Four zones were chosen for a several-day study following a Lagrangian approach: north of the South Orkney Islands (NSO), southeast of the South Orkney Island (SSO), northwest of South Georgia Island (NSG), and west of Anvers Island (WA). The position of the main hydrographic fronts along the cruise was determined following the scheme of Orsi et al. (1995) [53] (see Nunes et al. [36] for more details), with reference to the continuous records of temperature and salinity (thermosalinograph SBE 21 SeaCAT). Current velocity and direction were measured with the Shipboard Acoustic Doppler Current Profiler (SADCP) “Ocean Surveyor” at 75 khz, and the synoptic modeling data were obtained from the Global Real-Time Forecast System (Global RTOFS) [54]. In addition, we used 8-day average satellite images of chlorophyll-a concentration and sea surface temperature obtained from the Visible and Infrared Scanner (VIRS), NASA. In three of the zones (NSO, NSG, and WA), the studied water bodies were marked using WOCE (World Ocean Circulation Experiment) standard drifters provided with the Iridium communication system; in SSO, icebergs were used as Lagrangian markers. Seawater samples were collected at least once a day, at around 8:30 solar time (local time), with 12 L PVC Niskin bottles attached to a rosette sampling system. Samples were collected at several depths, from the surface (4 m) to 200 m; however, for the present study we used the 4 m and the DCM depth samples. A total of 50 stations were sampled for our purposes (16 at NSO, 14 at SSO, 11 at NSG, and 9 at WA, Figure 1A, and Tables S1 and S2).

2.2. Physicochemical Variables

Profiles of salinity, temperature, and fluorescence were obtained with a CTD probe (SeaBird 911 Plus) attached to the rosette. The CTD cast was deployed down to 200 m deep at all stations. Inorganic nutrients, fluorescent dissolved organic matter (FDOM), dimethyl sulfoniopropionate (DMSP), dimethyl sulfide (DMS), and isoprene concentrations were sampled at all stations, but only assessed at the surface (4 m). For inorganic nutrient (nitrate (NO3), nitrite (NO2), and ammonium (NH4+), silicate (SiO42−), and phosphate, P (PO43−)) concentrations, samples of 10 mL were collected from all stations and kept frozen (−20 °C) until analyzed in the laboratory of the Institut de Ciències del Mar, CSIC (ICM-CSIC) in an auto-analyzer (Bran + Luebbe AA3), following standard spectrophotometric methods [55]. For FDOM, we collected 250 mL sample in acid-cleaned glass bottles. The water sample was filtered through 0.2 µm acid-cleaned polycarbonate filters and analyzed onboard using a Perkin Elmer LS55 luminescence spectrometer. Single measurements for all surface samples were carried out in a 1 cm acid-cleaned quartz fluorescence cell at a constant room temperature and rinsed with a filtered water sample before analyses. We used the excitation/emission (Ex/Em) of 280/350 nm wavelengths to determine FDOM peak-T as an indicator of protein-like substances and of labile DOM [56,57]. Fluorescence measurements were expressed in Raman Units (R.U.), following Lawaetz and Stedmon [58], by normalization to the integrated water Raman scattering band of Milli-Q water freshly generated onboard every day. For DMSP analysis, two pellets of NaOH were added to 30 mL surface samples in gas-tight glass vials for 24-h hydrolysis to volatile DMS. Once in the ICM-CSIC lab, aliquots of 0.1 to 1 mL were injected into a purge flask with 2 mL of high purity water and purged for 4–6 min with ultrapure Helium. Volatiles were trapped in a Teflon tube loop submerged in liquid nitrogen, from where they were revolatilized by dipping the loop in hot water and injecting it into a gas chromatography system (GC-14A, Shimadzu). Sulfur compounds were separated using a packed CarbopackH 60/80 mesh column (Sigma-Aldrich, St. Louis, MO, USA) maintained at 170 °C and detected with a flame photometric detector. Retention time for DMS was 0.9 min, and the detection limit was 3 pmol. Analytical precision was better than 5%. Calibration was performed by syringe injection into the purge vial of varying volumes of a gaseous mixture of He and DMS released by a weight-calibrated permeation tube (Dynacal, Valco Instruments Co. Inc., Houston, TX, USA). DMSP concentrations were calculated by subtraction of the endogenous DMS. Isoprene and DMS were determined, along with other volatiles, on a gas chromatography-mass spectrometry system (5975-T LTM GC/MS, Agilent Technologies, Santa Clara, CA, USA). Aliquots of 25 mL were drawn from the glass bottle with a glass syringe with a Teflon tube and filtered through a 25 mm glass fiber filter while being introduced into a purge and trap system (Stratum, Tekmar Teledyne, Mason, OH, USA). Volatiles were stripped by bubbling with 40 mL min−1 of ultrapure He for 12 min, trapped on solid adsorbent (VOCARB 3000) at room temperature, and thermally desorbed (250 °C) into the GC. DMS and isoprene, monitored as m/z 62 and 67, respectively, in selected ion monitoring mode, had retention times of 2.4 and 2.5 min in the LTM DB-VRX chromatographic column held at 35 °C. The detection limit was 0.01 pmol (0.5 pmol L−1 in the samples), and the median analytical precision was 5%. DMS calibration was performed with DMS solutions generated by dissolution and hydrolysis of solid DMSP (TCI, Tokyo, Japan) in high purity water. Calibration of isoprene was performed by injections of a gaseous mixture of isoprene in N2 (provided by the University of California Irvine).

2.3. Microbiological Variables

The biomass of phytoplankton, abundances of viruses and prokaryotes, and prokaryotic heterotrophic production, were assessed in all stations (surface and DCM), while viral communities’ compositions were estimated in 20 stations at the surface and 17 at the DCM (Table S2). Chlorophyll a concentration (Chl-a) was determined on board in 250 mL sub-samples that were filtered through glass-fiber filters (Whatman GFF). This filtration was followed by extraction with 90% acetone of at least 24 h in the dark at 4 °C. The fluorescence of the extracts was measured with a calibrated fluorometer (Turner Designs), and no phaeophytin corrections were applied [59]. Microalgae were identified and quantified by pigment composition determined by HPLC [60]. CHEMTAX (chemical taxonomy software, version 1.95) was used to derive the contribution of microalgal groups for the total Chl-a biomass from pigment data [61]. Seven pigmentary classes were quantified: chlorophytes, cryptophytes, diatoms, dinoflagellates, haptophytes, prasynophytes, and pelagophytes (see Nunes et al. [36] for details). For our purposes, we only used the three most representative groups (cryptophytes, haptophytes, and diatoms) that each contributed with more than 20% each to the total chlorophyll a concentration in the visited areas. Sub-samples (2 mL) for viral abundance (VA, dsDNA-double stranded DNA- viruses) were fixed with glutaraldehyde (0.5% final concentration), flash-frozen in liquid nitrogen, and stored at −80 °C. Once in the ICM-CSIC laboratory, these sub-samples were thawed and stained with SYBR Green I before counting [62]. Groups of viruses were determined in bivariate scatter plots of green fluorescence of stained nucleic acids versus side scatter [42]. Viruses were classified as low (V1), medium (V2), and high (V3) according to fluorescence signal. Presumably, the V1 and V2 fractions are mainly attributed to bacteriophages, although some eukaryotic algal viruses may display low fluorescence, and V3 to phytoplankton viruses [42]. Prokaryote abundance (PA) was determined in 2 mL sub-samples preserved with a 1% paraformaldehyde +0.05% glutaraldehyde solution, flash-frozen in liquid N2, and stored at −80 °C until analysis. Once in the ICM-CSIC laboratory, samples were stained with SYTO13 (SYTOTM13, ThermoFisher) and counted following the protocol of Gasol and Del Giorgio [63]. Counts of viruses and prokaryotes were carried out in a flow cytometer (FACSCalibur, Becton & Dickinson). Samples for prokaryotes were run with fluorescent yellow-green latex beads (0.92-μm) as an internal standard, and for viral counts using a flow rate of ~60 µL min−1. Prokaryotic heterotrophic production (PHP) was estimated by the radioactive 3H-leucine incorporation technique [64], according to the modifications established for the use of microcentrifuge [65]. Leucine incorporation was assessed on board in a scintillation counter (Beckman). Then,
PHP = Leu × CF (μgC L−1 d−1),
where Leu is the 3H-leucine incorporation (pmol L−1 h−1), and CF is the conversion factor used: 1.5 kg C mol Leu−1 (cf. [66]).

2.4. Viral Community Composition: Random Amplified Polymorphic DNA (RAPD)

The randomly amplified polymorphic DNA-polymerase chain reaction (RAPD-PCR) is a practical and efficient method for measuring dsDNA viral diversity, providing viral assemblage comparisons through banding patterns. By gathering the required DNA concentration, we ensure that the analyses will give us a representative image of the original richness of the area and that the evenness distribution between samples will only respond to the initial diversity of the sample. Seawater samples (5L) from representative selected stations of the different regions and from two depths—surface and DCM (Table S2)—were sequentially filtered by 0.8 µm and 0.2 µm. Then, the 0.2 µm filtered samples were concentrated up to 50 mL with tangential flow filtration by using a spiral-wound cartridge (30-KDa VIVAFlow 200). Concentrates were stored at 4 °C until analysis. Once in the laboratory of the ICM-CSIC, the 50 mL concentrates were prefiltered through a polycarbonate filter (0.2 µm nucleopore) to eliminate any remaining prokaryote that would add prokaryotic DNA to the viral banding pattern, and ultra-concentrate by centrifuging at 4000 rpm for 10 min in 30-kDa ultra centrifugal filter tubes (Amicon, Sigma) [67]. Then, the viral ultra-concentrates were placed into agarose plug-like molding pieces. Before amplification of the viral DNA, the absence of the 16S rRNA prokaryotic gene was checked via PCR to ensure that each sample was prokaryotic-DNA free and, if it was not the case, it was treated with DNAse and further rechecked. Then, proteinase K was added in order to release the viral DNA from the protein capsid, and afterwards, treated with a cleaning protocol (Pefabloc®, Sigma), to inactivate this proteinase, as it can act as a Taq–polymerase inhibitor. Randomly amplified polymorphic DNA-polymerase chain reaction (RAPD–PCR) was designed (10′ at 94 °C, 30 times 3′ at 35 °C, 1′ at 72 °C and 30 s at 94 °C, 3′ at 35 °C, and 10′ at 72 °C) to amplify random sequences of the viral DNA by using the CRA22 primer (5′-CCG CAG CCA A- 3′). This primer was chosen due to its better results during the optimization of the protocol [68]. Consecutive RAPD-PCR reactions were done until 450 ng of DNA were obtained from each sample. DNA quantification was performed using a QuBit fluorometer. After that, the viral community composition was determined by using a gel electrophoresis technique. Briefly, the gel characteristics and running settings were 1.8% AG-2 at 70 Volts during 120 min and DNA was stained with 7 μL of SYBR Safe (10000X concentrate in DMSO). Finally, a UV image was taken from every gel with the gel imaging system ChemiDocTM (Bio-Rad) with a time exposure of around 2.1 s. Gel band pattern analysis was performed using the Quantity One 4.5.2 software (Bio-Rad). To make the band-distribution comparable between two different gels, before to final results we conducted banding pattern calibrations. Moreover, in the final gels, we ran replicates of some samples between gels in order to be able to couple and adjust them according to the duplicates and the ladder. Comparison of bands (operational taxonomic units, OTU) between lanes was performed by hand. Information was converted into a binary matrix representing the presence-absence of viral OTU; only communities present at more than two stations were considered in the analyses. The banding pattern analysis, i.e., the diversity analysis, was conducted with all samples, surface, and DCM, pooled together, so that we could identify which viral populations were present in both the surface and DCM.

2.5. Data Analyses

The map of the study regions, the temperature-salinity (T-S) diagram, and the temperature and salinity profiles were plotted using the software package Ocean Data View (ODV) [69]. Normality distribution of environmental and biological data (except for algal biomass) was checked with the Shapiro–Wilk test and, when necessary, continuous variables were logarithmically transformed prior to analyses to fulfill the requirements of statistical tests. The mean value of all variables was calculated considering all measurements performed at surface (n = 47–50) and at the DCM (n = 43–44). Spearman correlation analyses were used to determine the relationship between viral (total and V1-V3) abundances and biological and environmental variables. The Bonferroni–Holm correction was applied to each set of correlation coefficients. Multiple regression analysis was used to assess the predictor variables for the isoprene metabolite concentration. This was performed with the lm function from the R program. The best models are considered to be those that maximize the R2 and minimize the standard error of the estimate, with the minimum number of predictive variables. The distribution of abundances and concentrations of the biological variables and the analyses of the band pattern obtained with gel electrophoresis (see Section 2.4 for a detailed description of the technique) were plotted using the KaleidaGraph program. Distribution of viral community composition (i.e., similarity between band pattern among samples and sites) were plotted with non-metric multidimensional scaling plots (NMDS) using the metaMDS function in the R program. The dissimilarity matrix calculated using the Jaccard distance index and the stress factor is the statistical scale that indicates how good the 2-D representation of the original data is. Finally, the relationship between viral community composition and selected environmental variables was explored using constrained (or canonical) correspondence analyses (CCA), a multivariate ordination technique [70] in which the inertia can be interpreted as the explained variance. The group of environmental variables included in the CCA was chosen according to a step-wise selection, in a backward direction, with an AIC of 96.641. The CCA shows how the selected variables (represented as arrows) contribute to the structure of marine viral communities (represented as points colored according to the sampling zone). The lengths of the arrows are equivalent to their influence on the represented variables and exert in the viral communities (the longer the arrow, the stronger the influence). All statistical tests were performed using the vegan package from R [71].

3. Results

3.1. Characterization of the Study Area: Environmental and Biological Variables

3.1.1. Temperature and Salinity

The T–S diagram of the sampled stations revealed hydrographic differences among the visited regions (Figure 1B); this is also reflected in the vertical temperature and salinity profiles shown in Figure S1A,B. The main hydrographic fronts during the PEGASO cruise are described in Nunes et al. [36], following Orsi et al. [53]. The NSO and NSG regions were located within meanders of the southern boundary of the eastwards Antarctic circumpolar current (SBACC) and the polar front (PF). SSO was nearly 60 nautical miles north of the Weddell front (WF) and next to the marginal ice zone of the Weddell Sea in winter. WA was also placed on the southern boundary of the Antarctic circumpolar current (SBACC). NSG had the warmest waters, reaching 4.8 °C at the surface, and presented rather homogeneous salinity profiles (around 33.80, Figure S1A). NSO and SSO, in spite of their proximity, had different mean temperatures (seasonally warmed) and salinity (ice-melt water influenced), meaning SSO was significantly colder and less saline than NSO (Table 1, Figure S1). Finally, WA had also low salinity, probably was also affected by ice-melt water, and intermediate water temperatures (Figure S1B). The thermocline was appreciable at around 30 m (NSO), 16 m (SSO), 50 m (NSG), and 23 m (WA) (Figure S1).

3.1.2. Inorganic Nutrients, FDOM and Secondary Compounds

Mean concentrations of inorganic nutrients (NO3, NO2, NH4+, SiO42−, PO43−) and fluorescent dissolved organic matter (FDOM–peakT), as well as secondary metabolites (DMSP, DMS, and isoprene) measured at the surface, varied among sampling regions (Table 1). Nitrite and ammonia had almost negligible concentrations, compared to those of nitrate. Nitrate and phosphate concentration averages reached higher values in SSO and NSO (27.76–28.44 µM NO3 and 2.04–2.15 µM PO43−, respectively) than in NSG and WA (18.35–18.63 µM NO3−, and 1.42–1.81 µM PO43−, respectively). Then, silicate showed extremely low average values at NSG (mean: 2.02 ± 20.73 µM SiO42−, Table 1) and very high between 45.72 and 49.75 µM at the other regions. Although we did not measure micronutrients, there is evidence from prior studies that NSG is placed in an iron-sufficient zone [36,72,73]. The FDOM-peakT was remarkably higher at WA (40.8 ± 6.3 R.U, Table 1) than in the other areas. DMSP concentration was two-fold higher at NSO than in the rest of the regions (249.1 ± 101.2 nM, Table 1) and DMS was highest at SSO (7.8 ± 1.6 nM, Table 1). Isoprene reached the highest concentration at NSG (76.5 ± 8.3 nM) and the lowest at SSO (4.7 ± 1.6 nM, Table 1).

3.1.3. Chlorophyll a, Phytoplankton Taxa and Prokaryote Concentrations

Chl-a concentration was two-fold higher in NSG and WA than in NSO and was at nearly the limit of detection at SSO, both at the surface and at the DCM (Figure 2A, Table 1). The contribution of the dominant phytoplankton groups (cryptophytes, diatoms, and haptophytes) varied also among the sampling regions (Figure 2B,C). NSG was dominated by diatom biomass (83% of the total Chl-a-derived biomass of the three algal groups at the surface, and 85.6% at the DCM) (Table 1, Figure 2B,C). WA was dominated by cryptophytes (around 50%), followed by haptophytes (around 40%), at both depth layers (Table 1, Figure 2B,C). In NSO, the surface’s phytoplankton community showed co-dominance between cryptophytes (42.4%) and haptophytes (40.4%), with a lower presence of diatoms (17.2%), while at the DCM, diatoms and haptophytes co-dominated (44.1 and 36.7%, respectively). In the SSO, haptophytes were always dominant, reaching 78.7% of the total algal biomass (of the three groups) at the surface and 56.4% at the DCM (Table 1, Figure 2B,C). Prokaryote abundance (PA) showed the highest mean values at NSG, both at the surface and at the DCM (5.4 ± 0.9 × 105 prokaryote mL−1, and 5.2 ± 1.0 × 105 prokaryote mL−1, respectively (Table 1, Figure 2D)). PA differences between surface and DCM were only noticeable at SSO, where the DCM was enriched by a factor of 1.2 (Table 1, Figure 2D). Finally, PHP was highest in NSG and lowest in SSO (1.04 ± 0.44 µg C L−1 d−1 and 0.31 ± 0.18 µg C L−1 d−1, respectively (Table 1)).

3.2. Viral Abundance

Viral abundance (VA), like almost all the other variables, achieved the highest mean values at NSG (21.5 ± 10.7 × 106 viruses mL−1 at the surface and 17.9 ± 6.6 × 106 viruses mL−1 at the DCM (Figure 3A)) and the lowest at the surface of SSO (3.0 ± 1.5 × 106 viruses mL−1). Mean VA generally tended to be higher at the surface than at the DCM (Table 1), except for SSO (Table 1, Figure 3A). The contributions of low and medium fluorescence viral abundances (V1 and V2) to total viral abundance (VA) were similar (Figure 3B,C,E), except for NSG, where V2 was slightly higher than V1 (Table 1, Figure 3B,E). The V3 fraction (viruses with high fluorescence content, Figure 3 D) was less abundant at all regions, reaching on average ~7%, 6%, 10%, and 6% of the total VA at NSO, SSO, NSG, and WA, respectively (Figure 3E). Viral abundance was positively correlated with temperature, PA, Chl-a, PHP, diatom and haptophyte biomass, nitrite, FDOM, and isoprene concentrations (Table 2). Conversely, VA covaried negatively with nitrate, silicate, phosphate, and DMS concentrations (Table 2).

3.3. Viral Community Composition

We assessed viral diversity from 20 stations at the surface (8 in NSO, 6 in SSO, 3 in NSG, and 3 in WA) and from 17 at the DCM (7 in NSO, 5 in SSO, 2 in NSG, and 3 in WA). Considering all stations, we observed 22 distinct RAPD bands ranging in size from about 250 bp to 2000 bp (Figures S2 and S3). The banding pattern from each sample displayed in the agarose gels (Figure S2), corresponds to the most abundant populations of the viral community present at each station and gives an idea of the viral community composition. Thus, if we compare the band pattern of each lane (i.e., station), we can assess the viral diversity among and within regions and depths layers. NSG accounted for 76.7% of the detected bands, followed by NSO with 66.7%, WA with 60%, and SSO with 53.3%. From those, NSO was the region with a higher number of bands relative to the sampling effort (i.e., sampled stations per region). On average, half of the bands observed at the surface were also detected at the DCM of the same station (58% at NSO, 37% at SSO, 58% at NSG, and 61% at WA (Figure S3)). Several bands were only present at the NSO region (e.g., bands: 3, 9, 11, 14 (Figure S3)) suggesting some viral population endemism, while other bands (e.g., bands 15, 16, and 19 (Figure S3)) were found in different regions, pointing towards more cosmopolitan viral members. We grouped viral communities using a non-metric multidimensional scaling analysis (NMDS) (Figure 4), based on the RAPD’s band pattern, with a satisfactory stress factor of 0.171 (Figure 4 and Figure S4). Some viral communities were more similar within the same region than among regions, indicating that viruses differentiate regionally (Figure 4 and Figure S4). This can be especially appreciated with the viruses in the NSO, which form a compact group at both depth layers, revealing a highly homogeneous viral community, also appreciable for viruses from the surface at SSO (Figure 4A). In contrast, the viral communities of the WA and NSG regions were more widely dispersed, both at the surface and at the DCM (Figure 4A,B and Figure S4).
We used constrained correspondence analyses (CCA) (Figure 5) to explore the relationships between viral community structures and the resultant best explanatory variables according to step-wise selection. These variables were: temperature, salinity, nutrient concentrations (NO3−, SiO42−), PA, PHP, Chl-a, and the biomass of diatoms and haptophytes. In addition, secondary metabolites (DMSP and isoprene) were considered, in order to address the role of viral communities on the gas and particle exchange for aerosolization. The first two ordination axes of the CCA accounted for 62.9% of the variance in the viral community composition with a total constrained inertia of 77.8% (Figure 5). The first ordination axis of the CCA accounting for 33.8% of the variance was mainly related to temperature on the positive side and with nutrients on the negative one. Therefore, the stations of the more temperate regions (NSG and WA, Figure S1) appeared on the right side of the graph, while those of the colder regions were found on the left side (NSO and SSO (Figure S1)). The second axis accounted for a slightly smaller proportion of the variance (29.1%) and was positively related to salinity and algae biomass (most importantly in the WA region (Figure 5)); thus, most samples of SSO, the less saline region, occupied the negative side of this axis.

4. Discussion

Viral communities showed some segregation differences among the sampled regions, both at surface waters and at the DCM (Figure 4). However, while some viral populations were geographically constrained, others were widely distributed, as shown by the banding pattern of the RAPD (Figures S2 and S3). Indeed, several studies have reported a high global viral diversity and almost as high local diversity [74,75,76,77,78,79], as well as connectivity along the water column [80,81,82,83]. This might be based, on one hand, in the migration capability (i.e., transport by oceanic currents and sinking attached to particles) of marine viruses [84], and on the other hand, on environment selective pressure [15,17].
Authors such as Breitbart and Rohwer [75] suggested the seed-bank theory, where high local viral diversity relies on host availability and diversity. There, viruses would move from the seed-bank state to the active phase whenever their host would be “blooming”. Once the predominant host cells decay, the empty niche would be occupied by another hosting species, and the consequent displacement of viruses from the seed-bank state would occur [75]. These ecological dynamics may apply in the present study, explaining the obtained grouping of Antarctic viral community structures according to the sampled regions (in particular in the NSO and SSO regions (Figure 4)). In addition, Brum et al. [77] proposed temperature as a relevant environmental factor that influences changes in viral activity and diversity. Thus, during the spring–summer transition, they detected a reversion from lysogenic to lytic viruses that increased the viral and microbial diversity. All these processes will favor the geographical differentiation of viral assemblages.
In our study, temperature and salinity substantially varied among regions (Figure 1B and Figure S1), and the different degrees of stratification of the water column (with mixed layers from 16 m in SSO to 50 m in NSG) might have affected the light penetration and nutrient supply from deeper layers to the upper mixed layer, as is reported in Nunes et al. [36], for the studied area. In addition, it is expected that the development of different phytoplankton species [35] and prokaryotes [82,83] may influence the environmental conditions, modifying nutrient concentration and contributing to the production of secondary metabolites, such as DMSP, DMS, and isoprene [14,15,16,17,18,85,86]. At the same time, viruses could participate in these biogeochemical processes infecting phytoplankton and prokaryotes and releasing to the environment dissolved organic matter (DOM) from the cell-enclosed material [87,88]. It makes sense, then, that our results showed a strong significant correlation between fluorescent dissolved organic matter (FDOM, peak T, considered an indicator of proteinaceous compounds [89], and labile DOM [57]), and all viral abundance fractions (Table 2). Then, viruses would intervene in regulating microbial biomass and diversity, which in turn would be reflected in the viral community composition [25,26,27,84,88]. Viral abundance related positively with that of their potential hosts (prokaryotes and phytoplankton), but negatively with some inorganic nutrient concentrations (NO3, SiO42− and PO43− (Table 2)). This could be a consequence of inorganic nutrient uptake by algae at rates which vary among different phytoplankton taxa and depend on micronutrient availability [36]. Thus, in the iron-rich waters of NSG [36], where diatoms dominated the phytoplankton biomass (Figure 3), the inorganic nutrient concentrations, especially silicate, registered the lowest values (Table 1). In contrast, in the other regions (NSO, SSO, and WA), where iron was limited [3], inorganic nutrient concentrations were higher (Table 1).
The potential role of viruses in producing secondary metabolites as DMSP (DMS precursor) by lysing phytoplankton host cells, was experimentally reported by Hill et al. [26]. Then, we expected that viral abundance and DMSP concentration would correlate; however, that was not reflected in our results (Table 2). Bratbak et al. [25] also observed no correlation between viral abundance and DMSP concentration and posed that bacteria degradation represents the major sink for DMSP and DMS, which is also supported by Kiene and Service [90]. Furthermore, other DMSP-freeing pathways (i.e., grazing, senescence, apoptosis, etc.) may also play a role in increasing its release [20,22,23,91,92]. In the case of DMS, its concentration was negatively correlated with viruses. DMS has been shown to be strongly dependent on microbial community composition (e.g., presence of DMS consumers [15] or demethylating bacteria [17]) and to environmental conditions (e.g., wind-related ventilation and photolysis), which modify the structure of the microbial food web [15]. This may be the case of the area south of South Orkney (SSO) area, where the highest irradiance was registered, and thus, DMSP to DMS transformation could be enhanced. The formation of isoprene is mainly attributed to the activity and physiology of phytoplankton [28,29,30,31], while the effect of viral lysis is still unclear [14,18]. In our study, although isoprene concentration and viral abundance strongly covaried (Table 2), we believe that it is not a cause–effect. We are aware that diatom viruses are mainly ssRNA and ssDNA [93], and viral abundances reported here, by FCM, refer to dsDNA. Then, through multiple regression analyses, we obtained that the model that better predicted the isoprene concentration variability highlighted diatom biomass as the only explanatory factor and excluded the dsDNA viral abundance (log isoprene conc. = −3.53 (±1.43) + 0.37(±0.07) × log diatom conc.; n = 18, R2 = 0.78, p < 0.0001). Hence, further investigation is needed in studying the interaction between diatoms and their specific viruses in the formation of isoprene. Indeed, some studies with other photosynthetic microorganisms reported an increase of isoprene production after viral infection, as is the case of Prochlorococcus infected by phages [14,18].
In summary, the observed regional segregation of the different viral communities (Figure 4 and Figure 5) in this study, could be a result—to a greater or lesser extent—of the presence of different potential prokaryote and phytoplankton hosts, as observed by several authors [35,50,51,77,82,83], and to the variability of some physicochemical parameters [15,17]. Viral communities at the NSO and SSO regions were associated with the coldest waters and with high inorganic nutrient concentrations (Table 1, Figure 5). In particular, the group of viruses dominating the NSO region also coexisted with high DMSP concentration (Table 1) and a high proportion of haptophyte and cryptophyte biomass (Figure 2B,C), which according to Stefels et al. [16], are the main phytoplankton taxa responsible of the DMSP production. However, this was not reflected in the increase of DMS, perhaps due to low bacterial lyase activity, important in the DMSP degradation [91,92,94], or to environmental factors (i.e., temperature and sunlight etc. [15]). In the NSG and WA regions, although the number of samples was low, viral communities were more heterogeneous. Indeed, high temperatures and the availability of different dominant potential prokaryotic and phytoplanktonic hosts may influence the viral community composition. The high prokaryotic heterotrophic production and prokaryote abundance recorded at the NSG could be the results of a high growth rate of different successful prokaryote communities, derived from the carbon released by the various diatom taxa blooming during the sampling time [36]. High isoprene concentrations can be also associated with this bloom of diatoms, according to several studies [29,95,96]. Finally, regarding the variability of viral assemblages in the WA, relatively high temperatures could be associated with an advanced state of phytoplankton succession [36], which may have been driven by the high abundance and diversity of potential hosts such as haptophytes, cryptophytes, and prokaryotes [75].

5. Conclusions

In the present study, we found that the Antarctic and sub-Antarctic marine viral abundances are mainly associated with temperature and the abundance and biomass of their potential hosts. Furthermore, the viral community composition showed higher similarities within areas than among regions, particularly in surface waters. Presumably, the combination of physicochemical and biological factors led to distinct scenarios, which could drive variability in marine viral communities and simultaneously influence the various components of the microbial food web. However, these relationships must be interpreted in a dynamic way, because while viruses are causing changes in prokaryote and phytoplankton communities, at the same time it is the viruses that are changing in response to shifts in prokaryotic and phytoplankton community composition. While it is known that different phytoplankton taxa are associated with distinct prokaryote communities, further characterization of prokaryotic community composition would provide more insight on this potential diversification driver. Concerning secondary metabolites, we did not find a clear relationship between viruses and DMSP, DMS, and isoprene; thus, further experimental and field studies are needed in order to disentangle the complex relationships occurring in the marine food web that are involved in their release.

Supplementary Materials

The following are available online at https://www.mdpi.com/2075-1729/10/7/107/s1, Figure S1: Temperature and salinity profiles, Figure S2: Randomly amplified polymorphic DNA gel image, Figure S3: Viral community banding patter, Figure S4: NMDS of surface and DCM viral communities, Table S1: Physicochemical variables measured during the PEGASO cruise, Table S2: Biological variables measured during the PEGASO cruise.

Author Contributions

Conceptualization, R.S. and D.V.; data curation, A.S.-G., E.V.-D., and M.E. (Mikhail Emelianov); formal analysis, A.S.-G.; funding acquisition, M.D. and R.S.; investigation, A.S.-G., M.D., R.S., and D.V.; methodology, A.S.-G. and D.V.; project administration, R.S.; resources, P.C., Y.M.C., E.O.-R., and M.S.C.; software, A.S.-G.; supervision, D.V.; validation, I.F., M.E. (Marta Estrada) and M.S.; visualization, A.S.-G.; writing—original draft, A.S.-G. and D.V.; writing—review & editing, M.M.S., M.E. (Marta Estrada), C.M., S.N., and M.D. All authors have read and agreed to the published version of the manuscript.

Funding

The study was supported by the Spanish Ministry of Economy through projects PEGASO (CTM2012-37615) to R. Simó, and BIO-NUC (CGL2013–49020-R) to M. Dall’Osto.

Acknowledgments

We are very grateful to our colleagues P. Rodríguez-Ros, E. Borrull, C. Antequera, V. Balagué, C. Cardelús, and M. Zamanillo, who helped with sampling, and the crew of the R/V Bio-Hespérides for making possible the development of the field study, for their hospitality, and help with the logistics. “This research is part of POLARCSIC activities”.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Suttle, C.A. Viruses in the sea. Nature. 2005, 437, 356–361. [Google Scholar] [CrossRef] [PubMed]
  2. Danovaro, R.; Corinaldesi, C.; Dell’Anno, A.; Fuhrman, J.A.; Middelburg, J.J. Marine viruses and global climate change. FEMS Microbiol. Rev. 2011, 35, 993–1034. [Google Scholar] [CrossRef]
  3. Behrenfeld, M.J.; Kolber, Z.S. Widespread iron limitation of phytoplankton in the South Pacific Ocean. Science 1999, 283, 840–843. [Google Scholar] [CrossRef]
  4. Priddle, J.; Nedwell, D.B.; Whitehouse, M.; Reay, D.S.; Savidge, G.; Gilpin, L.; Murphy, E.; Ellis-Evansi, J.C. Re-examining the Antarctic Paradox: Speculation on the Southern Ocean as a nutrient-limited system. Ann. Glac. 1998, 27, 661–668. [Google Scholar] [CrossRef] [Green Version]
  5. Tréguer, P.; Jacques, G. Review Dynamics of nutrients and phytoplankton, and fluxes of carbon, nitrogen and silicon in the Antarctic Ocean. In Weddell Sea Ecology; Hempel, G., Ed.; Springer: Berlin/Heidelberg, Germany, 1992. [Google Scholar]
  6. Davison, B.; O’dowd, C.; Hewitt, C.N.; Smith, M.H.; Harrison, R.M.; Peel, D.A.; Wolf, E.; Mulvaney, R.; Schwikowski, M.; Baltenspergert, U. Dimethyl sulfide and its oxidation products in the atmosphere of the Atlantic and Southern Oceans. Atmosph. Environ. 1996, 30, 1895–1906. [Google Scholar] [CrossRef]
  7. Gowing, M.M.; Garrison, D.L.; Gibson, A.H.; Krupp, J.M.; Jeffries, M.O.; Fritsen, C.H. Bacterial and viral abundance in Ross Sea summer pack ice communities. Mar. Ecol. Prog. Ser. 2004, 279, 3–12. [Google Scholar] [CrossRef]
  8. Fuhrman, J.A. Marine viruses and their biogeochemical and ecological effects. Nature 1999, 399, 541–548. [Google Scholar] [CrossRef] [PubMed]
  9. Gobler, C.; Hutchins, A.; Fisher, N.S.; Cosper, E.M.; Sañudo-Wilhelmy, S.S. Release and bioavailability of C, N, P, Se, and Fe following viral lysis of a marine chrysophyte. Limnol. Oceanogr. 1997, 42, 1492–1504. [Google Scholar] [CrossRef] [Green Version]
  10. Poorvin, L.; Rinta-Kanto, J.M.; Hutchins, D.A.; Wilhelm, S. Viral release of Fe and its bioavailability to marine plankton. Limnol. Oceanog. 2004, 49, 1734–1741. [Google Scholar] [CrossRef] [Green Version]
  11. Mioni, C.E.; Poorvin, L.; Wilhelm, S.W. Virus and sidero- phore-mediated transfer of available Fe between hetero- trophic bacteria: Characterization using an Fe-specific bioreporter. Aquat. Microb. Ecol. 2005, 41, 233–245. [Google Scholar] [CrossRef] [Green Version]
  12. Evans, C.; Brussaard, C.P.D. Regional Variation in Lytic and Lysogenic Viral Infection in the Southern Ocean and Its Contribution to Biogeochemical Cycling. Appl. Environ. Microbiol. 2012, 78, 6741–6748. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Vogt, M.; Liss, P.S. Dimethylsulfide and climate. In Surface Ocean—Lower Atmosphere Processes; Le Quére, Ć., Saltzman, E.S., Eds.; AGU: Washington, DC, USA, 2009; Volume 187, pp. 197–232. [Google Scholar]
  14. Shaw, S.L.; Gantt, B.; Meskhidze, N. Production and emissions of marine isoprene and monoterpenes: A review. Adv. Meteor. 2010, 1–24. [Google Scholar] [CrossRef]
  15. Simó, R. Production of atmospheric sulfur by oceanic plankton: Biogeochemical, ecological and evolutionary links. Trends Ecol. Evol. 2001, 16, 287–294. [Google Scholar] [CrossRef]
  16. Stefels, J.; Steinke, M.; Turner, S.; Malin, G.; Belviso, S. Environmental constraints on the production and removal of the climatically active gas dimethylsulphide (DMS) and implications for ecosystem modelling. Biogeochemistry 2007, 83, 245–275. [Google Scholar] [CrossRef] [Green Version]
  17. Steiner, P.A.; Sintes, E.; Simó, R.; De Corte, D.; Pfannkuchen, D.M.; Ivancic, I.; Najdek, M.; Herndl, G.J. Seasonal dynamics of marine snow-associated and free-living demethylating bacterial communities in the coastal northern Adriatic Sea. Envion. Microb. Rep. 2019, 11, 699–707. [Google Scholar] [CrossRef]
  18. Shaw, S.L.; Chisholm, S.W.; Prinn, R.G. Isoprene production by Prochlorococcus, a marine cyanobacterium, and other phytoplankton. Mar. Chem. 2003, 80, 227–245. [Google Scholar] [CrossRef]
  19. Charlson, R.J.; Lovelock, J.E.; Andreae, M.O.; Warren, S.G. Oceanic phytoplankton, atmospheric sulphur, cloud albedo and climate. Nature 1987, 326, 655–661. [Google Scholar] [CrossRef]
  20. Archer, S.D.; Stelfox-Widdicombe, C.E.; Burkill, P.H.; Malin, G. A dilution approach to quantify the production of dissolved dimethylsulphoniopropionate and dimethyl sulphide due to microzooplankton herbivory. Aquat. Microb. Ecol. 2001, 23, 131–154. [Google Scholar] [CrossRef]
  21. Keller, M.D. Dimethyl sulfide production and marine phytoplankton: The importance of species composition and cell size. Biol. Oceanogr. 1989, 6, 375–382. [Google Scholar]
  22. Matrai, P.A.; Keller, M.D. Total organic sulfur and dimethylsulfoniopropionate in marine phytoplankton: Intracellular variations. Mar. Biol. 1994, 119, 61–68. [Google Scholar] [CrossRef]
  23. Laroche, D.; Vezina, A.F.; Levasseur, M.; Gosselin, M.; Stefels, J.; Keller, M.D.; Matrai, P.A.; Kwint, R.L.J. DMSP synthesis and exudation in phytoplankton: A modeling approach. Mar. Ecol. Prog. Ser. 1999, 180, 37–49. [Google Scholar] [CrossRef] [Green Version]
  24. Simó, R.; Saló, V.; Almeda, R.; Movilla, J.; Trepat, I.; Saiz, E.; Calbet, A. The quantitative role of microzooplankton grazing in dimethylsulfide (DMS) production in the NW Mediterranean. Biogeochem. 2018, 141, 125–142. [Google Scholar] [CrossRef] [Green Version]
  25. Bratbak, G.; Levasseur, M.; Michaud, S.; Cantin, G.; Fernández, E.; Heimdal, B.R.; Heldal, M. Viral activity in relation to Emiliania huxleyi blooms: A mechanism of DMSP release? Mar. Ecol. Prog. Ser. 1995, 128, 133–142. [Google Scholar] [CrossRef]
  26. Hill, R.W.; White, B.A.; Cottrell, M.T.; Dacey, J.W. Virus-mediated total release of dimethylsulfoniopropionate from marine phytoplankton: A potential climate process. Aquat. Microb. Ecol. 1998, 14, 1–6. [Google Scholar] [CrossRef]
  27. Malin, G.; Wilson, W.H.; Bratbak, G.; Liss, P.S.; Mann, N.H. Elevated production of dimethylsulfide resulting from viral infection of cultures of Phaeocystis pouchetii. Limnol. Oceanogr. 1998, 43, 1389–1393. [Google Scholar] [CrossRef]
  28. Bonsang, B.; Polle, C.; Lambert, G. Evidence for marine production of isoprene. Geoph. Res. Lett. 1992, 19, 1129–1132. [Google Scholar] [CrossRef]
  29. Milne, P.J.; Riemer, D.D.; Zika, R.G.; Brand, L.E. Measurement of vertical distribution of isoprene in surface seawater, its chemical fate, and its emission from several phytoplankton monocultures. Mar. Chem. 1995, 48, 237–244. [Google Scholar] [CrossRef]
  30. Broadgate, W.J.; Liss, P.S.; Penkett, S.A. Seasonal emissions of isoprene and other reactive hydrocarbon gases from the ocean. Geoph. Res. Lett. 1997, 24, 2675–2678. [Google Scholar] [CrossRef]
  31. Yokouchi, Y.; Li, H.J.; Machida, T.; Aoki, S.; Akimoto, H. Isoprene in the marine boundary layer (Southeast Asian Sea, eastern Indian Ocean, and Southern Ocean): Comparison with dimethyl sulfide and bromoform. J. Geophys. Res. Atm. 1999, 104, 8067–8076. [Google Scholar] [CrossRef]
  32. Booge, D.; Marandino, C.A.; Schlundt, C.; Palmer, P.I.; Schlundt, M.; Atlas, E.L.; Bracher, A.; Saltzman, E.S.; Wallace, D.W.R. Can simple models predict large scale surface ocean isoprene concentrations? Atmos. Chem. Phys. 2016, 16, 11807–11821. [Google Scholar] [CrossRef] [Green Version]
  33. Peloquin, J.A.; Smith, W.O. Phytoplankton blooms in the Ross Sea, Antarctica: Interannual variability in magnitude, temporal patterns, and composition. J. Geophys. Res. Ocean. 2007, 112, 1–12. [Google Scholar] [CrossRef] [Green Version]
  34. Liu, X.; Smith, W.O. Physiochemical controls on phytoplankton distributions in the Ross Sea, Antarctica. J. Mar. Syst. 2012, 94, 135–144. [Google Scholar] [CrossRef]
  35. Mangoni, O.; Saggiono, V.; Bolinesi, F.; Margiotta, F.; Budillon, G.; Cotroneo, Y.; Misic, C.; Rivaro, P.; Saggiorno, M. Phytoplankton blooms during austral summer in the Ross Sea, Antarctica: Driving factors and trophic implications. PLoS ONE 2017, 12, e0176033. [Google Scholar] [CrossRef] [Green Version]
  36. Nunes, S.; Latasa, M.; Delgado, M.; Simó, R.; Estrada, M. Phytoplankton community structure in contrasting ecosystems of the Southern Ocean and western Antarctic Peninsula. Deep-Sea-Res. I 2019, 151, 103059. [Google Scholar] [CrossRef] [Green Version]
  37. Buchan, A.; LeCleir, G.R.; Gulvik, C.A.; González, J.M. Master recyclers: Features and functions of bacteria associated with phytoplankton blooms. Nat. Rev. 2014, 12, 686–698. [Google Scholar] [CrossRef] [PubMed]
  38. Fandino, L.B.; Riemann, L.; Steward, G.F.; Long, R.A.; Azam, F. Variations in bacterial community structure during a dinoflagellate bloom analyzed by DGGE and 16S rDNA sequencing. Aquat. Microb. Ecol. 2001, 23, 119–130. [Google Scholar] [CrossRef]
  39. Fandino, L.B.; Riemann, L.; Steward, G.F.; Azam, F. Population dynamics of Cytophaga-Flavobacteria during marine phytoplankton blooms analyzed by real-time quantitative PCR. Aquat. Microb. Ecol. 2005, 40, 251–257. [Google Scholar] [CrossRef]
  40. Pinhassi, J.; Sala, M.M.; Havskum, H.; Peters, F.; Guadayol, O.; Malits, A.; Marrasé, C. Changes in bacterioplankton composition under different phytoplankton regimens. Appl. Environ. Microbiol. 2004, 70, 6753–6766. [Google Scholar] [CrossRef] [Green Version]
  41. Guixa-Boixereu, N.; Vaqué, D.; Gasol, J.M.; Sánchez-Cámara, J.; Pedrós-Alió, C. Viral distribution and activity in Antarctic waters. Deep Sea Res. II Trop. Stud. Ocean. 2002, 49, 827–845. [Google Scholar] [CrossRef]
  42. Evans, C.; Pearce, I.; Brussaard, C.P.D. Viral-mediated lysis of microbes and carbon release in the sub-Antarctic and Polar Frontal zones of the Australian Southern Ocean. Environ. Microbiol. 2009, 11, 2924–2934. [Google Scholar] [CrossRef]
  43. Christaki, U.; Lefèvre, D.; Georges, C.; Colombet, J.; Catala, P.; Courties, C.; Sime-Ngando, T.; Blain, S.; Obernosterer, I. Microbial food web dynamics during spring phytoplankton blooms in the naturally iron-fertilized Kerguelen area (Southern Ocean). Biogeosciences 2014, 11, 6739–6753. [Google Scholar] [CrossRef] [Green Version]
  44. Vaqué, D.; Boras, J.A.; Torrent-Llagostera, F.; Agustí, S.; Arrieta, J.M.; Lara, E.; Castillo, Y.N.; Duarte, C.M.; Sala, M.M. Viruses and protists induced-mortality of prokaryotes around the antarctic peninsula during the austral summer. Front. Microbiol. 2017, 8, 1–12. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  45. Wilkins, D.; Yau, S.; Williams, T.J.; Allen, M.A.; Brown, M.V.; DeMaere, M.Z.; Lauro, F.M.; Cavicchioli, R. Key microbial drivers in Antarctic aquatic environments. FEMS Microbial. Rev. 2013, 37, 303–333. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  46. López-Bueno, A.; Tamames, J.; Velázquez, D.; Moya, A.; Quesada, A.; Alcamí, A. High diversity of the viral community from an Antarctic lake. Science 2009, 326, 858–861. [Google Scholar] [CrossRef] [PubMed]
  47. Bouvy, M.; Bettarel, Y.; Bouvier, C.; Domaizon, I.; Jacquet, S.; Le Floc’h, E.; Montanié, H.; Mostajir, B.; Sime-Ngando, T.; Torréton, J.P.; et al. Trophic interactions between viruses, bacteria and nanoflagellates under various nutrient conditions and simulated climate change. Environ. Microb. 2011, 13, 1842–1857. [Google Scholar] [CrossRef]
  48. Maranger, R.; Vaqué, D.; Nguyen, D.; Hébert, M.P.; Lara, E. Pan-Arctic patterns of planktonic heterotrophic microbial abundance and processes: Controlling factors and potential impacts of warming. Prog. Oceanogr. 2015, 139, 221–232. [Google Scholar] [CrossRef]
  49. Finke, J.F.; Hunt, B.P.V.; Winter, C.; Carmack, E.C.; Suttle, C.A. Nutrients and other environmental factors influence virus abundances across oxic and hypoxic marine environments. Viruses 2017, 9, 152. [Google Scholar] [CrossRef]
  50. Kukkaro, P.; Bamford, D.H. Virus-host interactions in environments with a wide range of ionic strengths. Environ. Microbiol. Rep. 2009, 1, 71–77. [Google Scholar] [CrossRef]
  51. Jover, L.F.; Effler, T.C.; Buchan, A.; Wilhelm, S.W.; Weitz, J.S. The elemental composition of virus particles: Implications for marine biogeochemical cycles. Nat. Rev. Microbiol. 2014, 12, 519–528. [Google Scholar] [CrossRef]
  52. Miranda, J.A.; Culley, A.I.; Schvarcz, C.R.; Steward, G.F. RNA viruses as major contributors to Antarctic virioplankton. Environ. Microbiol. 2016, 18, 3714–3727. [Google Scholar] [CrossRef]
  53. Orsi, A.H.; Whitworth, T.; Nowlin, W.D., Jr. On the meridional extent and fronts of the Antarctic Circumpolar Current. Deep-Sea Res. I 1995, 42, 641–673. [Google Scholar] [CrossRef]
  54. Dall’Osto, M.; Ovadnevaite, J.; Paglione, M.; Beddows, D.C.O.; Ceburnis, D.; Cree, C.; Cortés, P.; Zamanillo, M.; Nunes, S.O.; Pérez, G.L.; et al. Antarctic sea ice region as a source of biogenic organic nitrogen in aerosols. Sci. Repp. 2017, 7, 1–10. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  55. Hansen, H.P.; Koroleff, F. Determination of nutrients. In Methods of Seawater Analysis, 3rd ed.; Grasshoff, K., Kremling, K., Ehrhardt, M., Eds.; John Wiley & Sons: Hoboken, NJ, USA, 1999. [Google Scholar]
  56. Coble, P.G. Characterization of marine and terrestrial DOM in seawater using excitation-emission matrix spectroscopy. Mar. Chem. 1996, 51, 325–346. [Google Scholar] [CrossRef]
  57. Yamashita, Y.; Tanoue, E. Chemical characterization of protein-like fluorophores in DOM in relation to aromatic amino acids. Mar. Chem. 2003, 82, 255–271. [Google Scholar] [CrossRef]
  58. Lawaetz, A.J.; Stedmon, C.A. Fluorescence Intensity Calibration Using the Raman Scatter Peak of Water. App. Spectr. 2009, 63, 936–940. [Google Scholar] [CrossRef]
  59. Strickland, J.D.H.; Parsons, T.R. A Practical Handbook of Seawater Analysis, 2nd ed.; Fisheries Research Board of Canada: Ottawa, ON, Canada, 1968. [Google Scholar]
  60. Latasa, M. A simple method to increase sensitivity for RP-HPLC phytoplankton pigment analysis. Limnol. Oceanogr. Methods 2014, 12, 46–53. [Google Scholar] [CrossRef] [Green Version]
  61. Mackey, M.D.; Mackey, D.J.; Higgins, H.W.; Wright, S.W. CHEMTAX—A program for estimating class abundances from chemical markers: Application to HPLC measurements of phytoplankton. Mar. Ecol. Prog. Ser. 1996, 144, 265–283. [Google Scholar] [CrossRef] [Green Version]
  62. Brussaard, C.P. Optimization of procedures for counting viruses by flow cytometry. Appl. Environ. Microbiol. 2004, 70, 1506–1513. [Google Scholar] [CrossRef] [Green Version]
  63. Gasol, J.M.; Del Giorgio, P.A. Using flow cytometry for counting natural planktonic bacteria and understanding the structure of planktonic bacterial communities. Sci. Mar. 2000, 64, 197–224. [Google Scholar] [CrossRef] [Green Version]
  64. Kirchman, D.; K’nees, E.; Hodson, R. Leucine incorporation and its potential as a measure of protein synthesis by bacteria in natural aquatic systems. Appl. Environ. Microbiol. 1985, 49, 599–607. [Google Scholar] [CrossRef] [Green Version]
  65. Smith, D.C.; Azam, F.A. Simple, economical method for measuring bacterial protein synthesis rates in seawater using 3H-leucine. Mar. Microb. Food Webs 1992, 6, 107–114. [Google Scholar]
  66. Kirchman, D.L. Incorporation of thymidine and leucine in the subarctic Pacific: Application to estimating bacterial production. Mar. Ecol. Prog. Ser. 1992, 1, 301–309. [Google Scholar] [CrossRef]
  67. Winget, D.M.; Wommack, K.E. Randomly amplified polymorphic DNA PCR as a tool for assessment of marine viral richness. Appl. Environ. Microbiol. 2008, 74, 2612–2618. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  68. Lara, E. Viruses in the Marine Environment: Community Dynamics, Phage-Host Interactions and Genomic Structure. Ph.D. Thesis, University of Las Palmas of Canary Island (ULPG), Las Palmas de Gran Canaria, Spain, 2014. [Google Scholar]
  69. Schlitzer, R. Ocean Data View. 2017. Available online: http://odv.awi.de (accessed on 15 January 2017).
  70. Borcard, D.; Gillet, F.; Legendre, P. Numerical Ecology with R; Springer: Heidelberg, Germany; Dublin, Ireland, 2018. [Google Scholar]
  71. R Core Team. R: A Language and Environment for Statistical Computing; R Foundation for Statistical Computing: Vienna, Austria, 2017; Available online: https://www.R-project.org/ (accessed on 20 January 2017).
  72. Martin, J.; Fitzwater, S.; Gordon, R. Iron Deficiency Limits Phytoplankton Growth in Antarctic Waters. Glob. Biogeochem. Cycles 1990, 4, 5–12. [Google Scholar] [CrossRef] [Green Version]
  73. Nielsdóttir, M.C.; Bibby, T.S.; Moore, C.M.; Hinz, D.J.; Sanders, R.; Whitehouse, M.; Korb, R.; Achterberg, E.P. Seasonal and spatial dynamics of iron availability in the Scotia Sea. Mar. Chem. 2012, 130, 62–72. [Google Scholar] [CrossRef]
  74. Short, S.M.; Suttle, C.A. Sequence analysis of marine virus communities reveals that groups of related algal viruses are widely distributed in nature. Appl. Environ. Microbiol. 2002, 68, 1290–1296. [Google Scholar] [CrossRef] [Green Version]
  75. Breitbart, M.; Rohwer, F. Here a virus, there a virus, everywhere the same virus? Trends Microbiol. 2005, 13, 278–284. [Google Scholar] [CrossRef]
  76. Angly, F.E.; Felts, B.; Breitbart, M.; Salamon, P.; Edwards, R.A.; Carlson, C.; Chan, A.M.; Haynes, M.; Kelley, S.; Liu, H.; et al. The marine viromes of four oceanic regions. PLoS Biol. 2006, 4, 11. [Google Scholar] [CrossRef]
  77. Brum, J.R.; Hurwitz, B.L.; Schofield, O.; Ducklow, H.W.; Sullivan, M.B. Seasonal time bombs: Dominant temperate viruses affect Southern Ocean microbial dynamics. ISME J. 2016, 10, 437–449. [Google Scholar] [CrossRef]
  78. Moreno-Pino, M.; De la Iglesia, R.; Valdivia, N.; Henríquez-Castilo, C.; Galán, A.; Díez, B.; Trefault, N. Variation in coastal Antarctic microbial community composition at sub-mesoscale: Spatial distance or environmental filtering? FEMS Microbiol. Ecol. 2016, 92. [Google Scholar] [CrossRef]
  79. Gregory, A.C.; Zayed, A.A.; Conceição-Neto, N.; Temperton, B.; Bolduc, B.; Alberti, A.; Ardyna, M.; Arkhipova, K.; Carmichael, M.; Cruaud, C.; et al. Marine DNA Viral Macro- and Microdiversity from Pole to Pole. Cell 2019, 177, 1109–1123. [Google Scholar] [CrossRef] [Green Version]
  80. De Corte, D.; Sintes, E.; Yokokawa, T.; Reinthaler, T.; Herndl, G.J. Links between viruses and prokaryotes throughout the water column along a North Atlantic latitudinal transect. ISME J. 2012, 6, 1566–1577. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  81. Frank, A.H.; Garcia, J.A.; Herndl, G.J.; Reinthaler, T. Connectivity between surface and deep waters determines prokaryotic diversity in the North Atlantic Deep Water. Envion. Microb. 2016, 18, 2052–2063. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  82. Sunagawa, S.; Coelho, L.P.; Chaffron, S.; Kultima, J.R.; Labadie, K.; Salazar, G.; Djahanschiri, B.; Zeller, G.; Mende, D.R.; Alberti, A.; et al. Structure and function of the global ocean microbiome. Science 2015, 348, 1261359. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  83. Hurwitz, B.L.; Brum, J.R.; Sullivan, M.B. Depth-stratified functional and taxonomic niche specialization in the ‘core’ and ‘flexible’ Pacific Ocean Virome. ISME J. 2015, 9, 472–484. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  84. Brum, J.R.; Ignacio-Espinoza, J.C.; Roux, S.; Doulcier, G.; Acinas, S.G.; Alberti, A.; Chaffron, S.; Cruaud, C.; de Vargas, C.; Gasol, J.M.; et al. Patterns and ecological drivers of ocean viral communities. Science 2015, 348, 1261498. [Google Scholar] [CrossRef] [Green Version]
  85. Ruiz-González, C.; Galí, M.; Gasol, J.M.; Simó, R. Sunlight effects on the DMSP-sulfur and leucine assimilation activities of polar heterotrophic bacterioplankton. Biogeochemistry 2012, 110, 57–74. [Google Scholar] [CrossRef]
  86. Johnson, W.M.; Kido-Soule, M.C.; Kujawinski, E.B. Evidence for quorum sensing and differential metabolite production by a marine bacterium in response to DMSP. ISME J. 2016, 10, 2304–2316. [Google Scholar] [CrossRef] [Green Version]
  87. Wilhelm, S.W.; Suttle, C.A. Viruses and Nutrient Cycles in the Sea—Viruses play critical roles in the structure and function of aquatic food webs. Bioscience 1999, 49, 781–788. [Google Scholar] [CrossRef] [Green Version]
  88. Suttle, C.A. Marine viruses—Major players in the global ecosystem. Nat. Rev. Microb. 2007, 5, 801–812. [Google Scholar] [CrossRef]
  89. Gabor, R.S.; Baker, A.; Mcknight, D.M.; Miller, M. Fluorescence indices and their interpretation. In Aquatic Organic Matter Fluorescence; Coble, P., Lead, J., Baker, A., Reynolds, D., Spencer, R., Eds.; Cambridge University Press: Cambridge, UK, 2014. [Google Scholar]
  90. Kiene, R.P.; Service, S.K. Decomposition of dissolved DMSP and DMS in estuarine waters: Dependence on temperature and substrate concentration. Mar. Ecol. Prog. Ser. 1991, 76, 1–11. [Google Scholar] [CrossRef]
  91. Belviso, S.; Kim, S.K.; Rassoulzadegan, F.; Krajka, B.; Nguyen, B.C.; Mihalopoulos, N.; Buat-Menard, P. Production of dimethylsulfonium propionate (DMSP) and dimethylsulfide (DMS) by a microbial food web. Limnol. Oceanogr. 1990, 35, 1810–1821. [Google Scholar] [CrossRef]
  92. Kiene, R.P.; Linn, L.J. Distribution and turnover of dissolved DMSP and its relationship with bacterial production and dimethylsulphide in the Gulf of Mexico. Limnol. Oceanogr. 2000, 45, 849–861. [Google Scholar] [CrossRef]
  93. Tomaru, Y.; Nagasaki, K. Flow Cytometric Detection and Enumeration of DNA and RNA Viruses Infecting Marine Eukaryotic Microalgae. J. Oceanogr. 2007, 63, 215–221. [Google Scholar] [CrossRef]
  94. Stefels, J.; van Leeuwe, M.A. Effects of iron and light stress on the biochemical composition of Antarctic Phaeocystis sp. (prymnesiophyceae). i. intracellular dmsp concentrations. J. Phycol. Phycol. 1998, 34, 486–495. [Google Scholar] [CrossRef]
  95. Moore, R.M.; Oram, D.E.; Penkett, S.A. Production of isoprene by marine phytoplankton cultures. Geophys. Res. Lett. 1994, 21, 2507–2510. [Google Scholar] [CrossRef]
  96. McKay, W.A.; Turner, M.F.; Jones, B.M.R.; Halliwell, C.M. Emissions of hydrocarbons from marine phytoplankton—Some results from controlled laboratory experiments. Atmos. Envion. 1996, 30, 2583–2593. [Google Scholar] [CrossRef]
Figure 1. (A) Map of the location of the study regions: north of South Orkney, NSO; south of South Orkney, SSO; north of South Georgia, NSG; west of Anvers, WA. For each region, a zoomed window shows the position and number of the stations sampled. (B) Salinity–temperature diagram for the different sampled regions.
Figure 1. (A) Map of the location of the study regions: north of South Orkney, NSO; south of South Orkney, SSO; north of South Georgia, NSG; west of Anvers, WA. For each region, a zoomed window shows the position and number of the stations sampled. (B) Salinity–temperature diagram for the different sampled regions.
Life 10 00107 g001
Figure 2. Distribution of the microbiological variables at the surface and DCM in the four visited areas. (A) Chlorophyll a concentration. (B,C) Percentage of haptophyte, cryptophyte, and diatom biomass, with respect to the sum of the three phytoplankton groups biomass, at the surface and at DCM, respectively. (D) Prokaryote abundance. See explanation in Figure 1 for acronyms.
Figure 2. Distribution of the microbiological variables at the surface and DCM in the four visited areas. (A) Chlorophyll a concentration. (B,C) Percentage of haptophyte, cryptophyte, and diatom biomass, with respect to the sum of the three phytoplankton groups biomass, at the surface and at DCM, respectively. (D) Prokaryote abundance. See explanation in Figure 1 for acronyms.
Life 10 00107 g002
Figure 3. Distribution of viral abundance in the four visited areas, at the surface and DCM. (A) Total viral abundance. (BD) Abundance of different viral groups (V1, V2, and V3). (E) Percentage of viral groups in relation to total viral abundance. Presumably, fractions of V1 and V2 are mainly attributed to bacteriophages and V3 to viruses of phytoplankton. See explanation in Figure 1 for acronyms.
Figure 3. Distribution of viral abundance in the four visited areas, at the surface and DCM. (A) Total viral abundance. (BD) Abundance of different viral groups (V1, V2, and V3). (E) Percentage of viral groups in relation to total viral abundance. Presumably, fractions of V1 and V2 are mainly attributed to bacteriophages and V3 to viruses of phytoplankton. See explanation in Figure 1 for acronyms.
Life 10 00107 g003
Figure 4. Non-metric multidimensional scaling (NMDS) plots of the marine viral community composition. (A) Surface viral diversity. (B) DCM viral diversity. Plotted data arise from a matrix of similarity drawn by the banding pattern of randomly amplified polymorphic DNA (RAPD)-PCR products obtained for the pooled set of samples of the surface and DCM (Figure S2). Full dots represent surface communities and empty circles represent DCM samples. Numbers correspond to the station sampled. See explanation in Figure 1 for acronyms.
Figure 4. Non-metric multidimensional scaling (NMDS) plots of the marine viral community composition. (A) Surface viral diversity. (B) DCM viral diversity. Plotted data arise from a matrix of similarity drawn by the banding pattern of randomly amplified polymorphic DNA (RAPD)-PCR products obtained for the pooled set of samples of the surface and DCM (Figure S2). Full dots represent surface communities and empty circles represent DCM samples. Numbers correspond to the station sampled. See explanation in Figure 1 for acronyms.
Life 10 00107 g004
Figure 5. Ordination diagram displaying the canonical (or constrained) correspondence analyses (CCA). Surface marine viral communities from 20 sites, related to environmental (temperature, salinity, NO3−, and SiO42−) and biological (Chl-a, PA, PHP, diatoms, and haptophytes) variables, and to secondary metabolites (DMSP and isoprene concentrations). Variables are shown as arrows, and viral community data as dots. Longer arrows denote a higher significant relationship. Same dots coloring as used in Figure 4A. Blue numbers correspond to the stations sampled. See explanation in Figure 1 for acronyms.
Figure 5. Ordination diagram displaying the canonical (or constrained) correspondence analyses (CCA). Surface marine viral communities from 20 sites, related to environmental (temperature, salinity, NO3−, and SiO42−) and biological (Chl-a, PA, PHP, diatoms, and haptophytes) variables, and to secondary metabolites (DMSP and isoprene concentrations). Variables are shown as arrows, and viral community data as dots. Longer arrows denote a higher significant relationship. Same dots coloring as used in Figure 4A. Blue numbers correspond to the stations sampled. See explanation in Figure 1 for acronyms.
Life 10 00107 g005
Table 1. Mean ± standard deviation (SD) of environmental and biological variables measured at all stations. Mean of variables only available from surface (n = 47−50 stations): nutrients: nitrate (NO3), nitrite (NO2), and ammonia (NH4+), silicate (SiO42−), and phosphate (PO43-), fluorescent dissolved organic matter (FDOM-peakT), and secondary metabolites: dimethyl sulfoniopropionate (DMSP), dimethyl sulfide (DMS) and Isoprene at the different sampling areas. Concentrations are expressed in: * µM, ** nM, and *** pM. Mean of variables available from both surface (n = 47–50) and DCM (n = 43–44): temperature, salinity, chlorophyll a concentration (Chl-a, µg L−1), biomasses of cryptophytes, diatoms, and haptophytes (ng Chl-a L−1), prokaryote abundance (PA, 105 cells mL−1), prokaryotic heterotrophic production (PHP, µg C L−1 d−1), viral abundances (106 virus mL−1): total (VA), low and medium (V1, V2) and high fluorescence (V3). See explanation of Figure 1 for region acronyms.
Table 1. Mean ± standard deviation (SD) of environmental and biological variables measured at all stations. Mean of variables only available from surface (n = 47−50 stations): nutrients: nitrate (NO3), nitrite (NO2), and ammonia (NH4+), silicate (SiO42−), and phosphate (PO43-), fluorescent dissolved organic matter (FDOM-peakT), and secondary metabolites: dimethyl sulfoniopropionate (DMSP), dimethyl sulfide (DMS) and Isoprene at the different sampling areas. Concentrations are expressed in: * µM, ** nM, and *** pM. Mean of variables available from both surface (n = 47–50) and DCM (n = 43–44): temperature, salinity, chlorophyll a concentration (Chl-a, µg L−1), biomasses of cryptophytes, diatoms, and haptophytes (ng Chl-a L−1), prokaryote abundance (PA, 105 cells mL−1), prokaryotic heterotrophic production (PHP, µg C L−1 d−1), viral abundances (106 virus mL−1): total (VA), low and medium (V1, V2) and high fluorescence (V3). See explanation of Figure 1 for region acronyms.
AREAS NSOSSONSGWA
VARIABLESDepthsMeanSDMeanSDMeanSDMeanSD
NO3− *Surface27.761.8928.443.4218.355.4718.630.89
NO2− *Surface0.230.060.170.040.270.050.190.03
NH4+ *Surface1.113.511.031.030.782.022.481.87
SiO42−Surface47.214.0845.7211.962.0220.7349.753.66
PO43− *Surface2.040.212.150.281.420.431.810.16
FDOM (R.U.)Surface18.54.97.91.315.43.140.86.3
DMSP **Surface249.1101.2103.557.284.428.0115.714.61
DMS **Surface7.02.67.81.65.81.02.10.6
Isoprene ***Surface9.61.74.71.676.58.312.00.6
TemperatureSurface0.50.3−0.80.14.80.51.40.1
DCM0.00.3−1.50.44.30.41.70.6
SalinitySurface33.880.07333.140.06433.740.00733.410.002
DCM34.020.06834.200.14833.770.00433.420.053
Chl-aSurface1.820.330.310.055.051.884.450.46
DCM1.610.380.370.187.513.273.530.63
CryptophytesSurface288.1110.511.14.1628.2714.21898.7268.2
DCM117.8112.49.434.8034.4719.67676.5203.0
DiatomsSurface116.958.1710.94.941983.91451109.931.81
DCM272.2102.256.8250.4624021321.2139.433.33
HaptophytesSurface274.280.7580.7618.15373.5127.4629.5199.5
DCM225.750.4585.8541.23369.5150.8583.1151.1
PASurface2.40.52.90.65.40.93.40.9
DCM2.60.53.61.45.21.03.30.5
PHPSurface0.460.170.310.181.040.440.400.13
DCM0.300.100.270.131.240.410.290.05
VASurface4.52.03.01.521.510.713.54.1
DCM3.11.74.01.517.96.612.07.2
V1Surface2.41.11.60.88.74.16.82.1
DCM1.60.81.70.36.92.76.03.9
V2Surface1.80.81.20.710.85.96.12.0
DCM1.30.71.60.38.33.95.33.0
V3Surface0.30.10.20.12.21.20.70.3
DCM0.20.20.30.11.70.80.70.3
Table 2. Spearman correlation coefficients between total viral abundance (VA) and different viral abundance fractions (V1, V2, V3) with physicochemical and biological variables measured at surface (n = 47–50). Significant values after the Bonferroni–Holm correction (for a significance level of 0.05) are highlighted in bold. See Table 1 for units. Temp: temperature; Chl-a: chlorophyll a; PA: prokaryote abundance; PHP: prokaryotic heterotrophic production; Crypto: cryptophytes; Hapto: haptophytes; NO3−, NO2−, NH4+, SiO42−, PO4, FDOM: fluorescent dissolved organic matter -peakT; DMSP: dimethyl sulfoniopropionate, DMS: dimethyl sulfide.
Table 2. Spearman correlation coefficients between total viral abundance (VA) and different viral abundance fractions (V1, V2, V3) with physicochemical and biological variables measured at surface (n = 47–50). Significant values after the Bonferroni–Holm correction (for a significance level of 0.05) are highlighted in bold. See Table 1 for units. Temp: temperature; Chl-a: chlorophyll a; PA: prokaryote abundance; PHP: prokaryotic heterotrophic production; Crypto: cryptophytes; Hapto: haptophytes; NO3−, NO2−, NH4+, SiO42−, PO4, FDOM: fluorescent dissolved organic matter -peakT; DMSP: dimethyl sulfoniopropionate, DMS: dimethyl sulfide.
TempSalinityPAChl-aPHPCrypto.DiatomsHapto.NO3−NO2−NH4+SiO42−PO43−FDOMDMSPDMSIsoprene
VA0.780.140.500.780.530.430.780.65−0.820.38−0.05−0.50−0.700.47−0.18−0.530.78
V10.740.140.470.690.410.380.700.52−0.670.24−0.08−0.51−0.540.39−0.19−0.500.64
V20.780.140.540.780.500.400.770.64−0.800.35−0.07−0.49−0.700.44−0.21−0.510.74
V30.880.310.580.830.540.410.850.70−0.820.50−0.11−0.50−0.690.43−0.17−0.520.78

Share and Cite

MDPI and ACS Style

Sotomayor-Garcia, A.; Montserrat Sala, M.; Ferrera, I.; Estrada, M.; Vázquez-Domínguez, E.; Emelianov, M.; Cortés, P.; Marrasé, C.; Ortega-Retuerta, E.; Nunes, S.; et al. Assessing Viral Abundance and Community Composition in Four Contrasting Regions of the Southern Ocean. Life 2020, 10, 107. https://doi.org/10.3390/life10070107

AMA Style

Sotomayor-Garcia A, Montserrat Sala M, Ferrera I, Estrada M, Vázquez-Domínguez E, Emelianov M, Cortés P, Marrasé C, Ortega-Retuerta E, Nunes S, et al. Assessing Viral Abundance and Community Composition in Four Contrasting Regions of the Southern Ocean. Life. 2020; 10(7):107. https://doi.org/10.3390/life10070107

Chicago/Turabian Style

Sotomayor-Garcia, Ana, Maria Montserrat Sala, Isabel Ferrera, Marta Estrada, Evaristo Vázquez-Domínguez, Mikhail Emelianov, Pau Cortés, Cèlia Marrasé, Eva Ortega-Retuerta, Sdena Nunes, and et al. 2020. "Assessing Viral Abundance and Community Composition in Four Contrasting Regions of the Southern Ocean" Life 10, no. 7: 107. https://doi.org/10.3390/life10070107

APA Style

Sotomayor-Garcia, A., Montserrat Sala, M., Ferrera, I., Estrada, M., Vázquez-Domínguez, E., Emelianov, M., Cortés, P., Marrasé, C., Ortega-Retuerta, E., Nunes, S., M. Castillo, Y., Serrano Cuerva, M., Sebastián, M., Dall’Osto, M., Simó, R., & Vaqué, D. (2020). Assessing Viral Abundance and Community Composition in Four Contrasting Regions of the Southern Ocean. Life, 10(7), 107. https://doi.org/10.3390/life10070107

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop