Next Article in Journal
Consistent Site-Specific Foraging Behaviours of Yellow-eyed Penguins/Hoiho Breeding on Stewart Island, New Zealand
Previous Article in Journal
The ABC Model of Happiness—Neurobiological Aspects of Motivation and Positive Mood, and Their Dynamic Changes through Practice, the Course of Life
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Pelagic Bacteria and Viruses in a High Arctic Region: Environmental Control in the Autumn Period

by
Vladimir G. Dvoretsky
*,
Marina P. Venger
,
Anastasya V. Vashchenko
,
Tatyana M. Maksimovskaya
,
Tatyana G. Ishkulova
and
Veronika V. Vodopianova
Murmansk Marine Biological Institute of the Russian Academy of Sciences (MMBI RAS), 183010 Murmansk, Russia
*
Author to whom correspondence should be addressed.
Biology 2022, 11(6), 845; https://doi.org/10.3390/biology11060845
Submission received: 21 April 2022 / Revised: 28 May 2022 / Accepted: 30 May 2022 / Published: 31 May 2022
(This article belongs to the Section Ecology)

Abstract

:

Simple Summary

Microbial plankton represent a pivotal component of aquatic ecosystems worldwide. Despite their importance in pelagic ecosystems, microbial assemblages are less studied, in many Arctic regions. The study was concerned with the abundance of marine bacteria and viruses, in relation to a set of environmental variables, after the main productive season. We found a mosaic horizontal distribution of the microbial plankton, while the number of bacteria and viruses decreased with depth. Nutrients and zooplankton carbon were significant drivers of microbial abundance. Bacterioplankton abundance was positively correlated with counts of viruses, indicating close relations between these groups of microbes. Comparisons with previous studies revealed strong seasonal variations in the total abundance of marine bacteria, with the maximum values being in the summer and spring periods, while virioplankton counts were comparable in the autumn and spring–summer seasons. Our study is the first report regarding microbial plankton in the northeastern Barents Sea, in the autumn period, and provides baseline information, expanding current knowledge on the structure of pelagic Arctic ecosystems.

Abstract

In the marine environment, bacteria and viruses play a significant role in carbon fluxes, remineralization processes, and the infection of various organisms. We performed a survey in the northeastern Barents Sea, a region adjacent to the Arctic Ocean, to investigate spatial patterns of microbial plankton, after the main productive period, in October 2020. Two main water masses occurred in the study region—colder Arctic Water and warmer Barents Sea Water, representing transformed Atlantic Water. Multivariate analyses detected patchiness in the horizontal distribution of bacteria and viruses, and their abundances showed no clear association with the water masses. There was an obvious vertical pattern in microbial concentration, with the highest estimates in the upper layers. Surface viral and bacterial abundance varied in a wide range (2.20 × 105–10.7 × 105 cells·mL−1 and 0.86 × 106–14.98 × 106 particles·mL−1, respectively) and were correlated with each other. Bacterioplankton was dominated by small-sized cells (<2 μm, 0.04–0.06 µm3), and the average volume of bacterial cells tended to increase toward the seafloor. The ratio of viral to bacterial abundance (VBR) was 11 ± 1 and did not differ between water masses and depth layers. VBR were higher, compared to summer values, suggesting a strong impact of viruses on bacterioplankton, after the main productive season. Redundancy and correlation analyses showed that inorganic nutrients (nitrate and phosphate) and organic carbon from zooplankton were most responsible for the total variability in the microbial parameters. Water temperature and salinity, also, had a measurable impact, but their influence was lower. Bacterial abundance was lower than in other seasons and regions of the Barents Sea, while viral abundance was comparable, suggesting a stronger viral impact on Arctic marine bacteria in the autumn season.

1. Introduction

Arctic ecosystems are very sensitive to changes in physical conditions and climatic variability; here, it is proposed to be more crucial, relative to lower latitudes [1,2]. There is a set of physical factors making the Arctic marine environment unique, namely low temperatures, ice cover, the presence of large shallow continental shelves, clear seasonal fluctuations in the light regime, and a high supply of freshwater discharge, originating from rivers and melting ice [2,3]. Strong climatic variations have been documented in the Arctic, since the end of the 20th century [3,4,5]. Climate forcing was reported to affect all trophic levels, from microbes to seabirds and marine mammals [5,6,7]. Studying ecological responses to climatic variation has a great significance in the sustainability of human activities and management of Arctic ecosystems.
The Barents Sea may be considered the most productive Arctic region [1,2]. The southern and southwestern areas represent the largest cod-fishery areas, located in the Arctic [1,2]. The Barents Sea is a transition zone, for Atlantic water entering the Arctic Ocean from the Norwegian Sea.
The pelagic ecosystem function and structure in the Barents Sea are, largely, determined by the hydrological regime, local circulation, frontal zones, seafloor topography, and general currents [1]. Recent climatic fluctuations had major impacts on marine ecosystems in the region, including northward expansion of boreal species, variability in the total productivity, and changes in the pelagic community structure [3,5]. Plankton assemblages play a key role in carbon fluxes, in ecosystems of the Barents Sea. Being sensitive to environmental influence, plankton communities are good indicators of climatic changes in Arctic seas [3,5,6,7].
In the pelagic ecosystems, bacterioplankton biomass is known to be a substantial fraction and, in some regions with low chlorophyll a concentrations, it can be higher than phytoplankton stocks [8,9]. Primary production and subsequent availability of dissolved organic matter promote heterotrophic bacterial growth and abundance [10]. Phytoplankton biomass is grazed upon by herbivorous zooplankton communities and, then, dead plankton as well as various excretes and metabolites are utilized by marine bacteria, serving microzooplankton (ciliates, tintinnids, and heterotrophic flagellates) as food resources [11]. Bacterioplankton assemblages are responsible for exchanging active metabolites (vitamin B12, antibiotics, and growth promoters) that may regulate algal–bacterial interactions [10,11]. Nutrient release to the seawater, due to decomposition of animals and plants, primarily by bacterial action as well as bacterial metabolites, is, thus, of great significance in allowing further photosynthetic cycles in the Barents Sea [1].
Viruses are abundant and active components of marine ecosystems, which, primarily, infect the numerically dominant bacterioplankton and eukaryotic marine microbes [12]. Being host-specific pathogens, they may affect bacterioplankton abundance and composition [13]. Viral lysis of host cells provides substrates for bacterioplankton and can modify energy and matter fluxes in the microbial plankton [12,13]. Therefore, the activity of marine viruses may be considered a significant factor, determining the productivity of higher trophic levels and the efficiency of organic carbon transfer, in the marine environment worldwide [14]. The abundance and distribution of marine microbes and their contributions to pelagic-food-web dynamics in polar waters have been, intensively, investigated during the past decades [15,16,17,18,19,20]. Information regarding the standing stocks of microbial plankton in the Barents Sea is, relatively, sparse. There are few studies dealing with bacterioplankton and virioplankton ecology in the region [21,22,23,24,25,26]. However, previous studies were focused on the western, central, and southern regions, in the spring and summer periods. The northern part of the Barents Sea is less studied and there are no data on the microbial plankton in the autumn period, after the main productive season.
The aims of the present study were to (1) describe the horizontal and vertical distribution of virioplankton and bacterioplankton, (2) test a hypothesis that microbial plankton is associated with different water masses, and (3) reveal which factors are responsible for spatial variations in microbial abundance and biomass, in the northeastern Barents Sea in autumn.

2. Materials and Methods

2.1. Sampling

Water samples were collected at nine stations at two transects located in the northeastern Barents Sea, during a research cruise of the R/V Dalnie Zelentsy, between 8 and 11 November 2020 (Figure 1, Table 1). Vertical profiles of water temperature and salinity were recorded with a CTD Sealogger (SeaBird SeaCat SBE-19, CTD) [27]. Water masses were distinguished, based on hydrological criteria (Figure S1) [2,4]. Water samples were collected at 5–6 depths (0 m, 10 m, 25 m, 50 m, 100 m and near the bottom) at each station, using 10 L Niskin water bottles, attached to a rosette. Subsamples were taken from each bottle, for the determination of nitrate (NO3), phosphate (PO43−), and dissolved oxygen, as well as bacterial and viral abundance, for each sampling depth. Chlorophyll a concentrations were determined, for the upper 0 m and 10 m layers. Considering the significance of plankton animals, as a potential organic substrate for microbial growth and development, we collected zooplankton in the upper 100 m layer, using a Juday net (180 µm mesh size and 0.11 m2 mouth opening). Zooplankton sampling was performed, in accordance with our previous studies [28,29,30,31].

2.2. Processing

Samples for nutrient determination were collected into acid-cleaned polyethylene tubes, after thorough rinsing and filtration, through a 45 µm polycarbonate filter (FMPA, Vladisart, Vladimir, Russia) inserted in a filter holder. Samples for nutrient analysis were collected in 20 mL vials and, immediately, frozen in liquid nitrogen and, then, stored at −20 °C to further analysis. Determination of dissolved-oxygen concentration was carried out using the Winkler method with an automatic-endpoint-detection burette Digital Burette VITLAB (63762 Grossostheim, Bavaria, Germany). Analyses of dissolved inorganic nutrients were carried out at the Murmansk Marine Biological Institute laboratory, according to standard methods [32].
Water samples for chlorophyll a determination with a volume of 5 L for each sampling layer were filtered, immediately, through 0.6 µm Vladiopore filters. The filters were kept frozen until the onboard analysis. The filters were, subsequently, extracted in 90% acetone, placed in a freezer at 4 °C for 24 h, and centrifuged, following the standard procedure [33]. Chlorophyll a concentrations were measured, using a Nicolett Evolution 500 spectrophotometer (Spectronic Unicam, Scotia, NY, USA), which had been calibrated with commercially purified chlorophyll a preparations.
In the laboratory, the zooplankton organisms were identified and counted under an MBS-10 stereomicroscope. Abundance was expressed as individuals·m−3. Calculations of biomass of particular taxa were made using published mean individual wet, dry, or carbon weights as well as length/weight relationships [34,35,36]. All values were computed as mg carbon mass (DM) per cubic meter, using the following relationship: 1 mg wet mass = 0.2 mg dry weight = 0.1 mg C [37]. The conversion of wet to dry weight, for Ctenophora and medusae, assumed 1 mg wet mass = 0.04 mg dry weight [37].
Water samples for marine bacteria were preserved with particle-free 40% formaldehyde (final concentration of 2%) and, immediately, stored in slide boxes at −20 °C, for further microscopic analysis in the laboratory. Sub-samples (20 to 100 mL) were stained with DAPI for 8–10 min and filtered onto 0.2 µm Nuclepore filters using a low vacuum [38]. The blackened Nuclepore filters were supported by GF/F filters, during the filtration process, to facilitate a homogenous distribution of cells. Filtered samples were mounted on slides in immersion oil and stored at −30 °C, in preparation for microscopic examination. Microbial cells were enumerated, using an OlympusBX 53 epifluorescence microscope, with a barrier filter for UV-light excitation. Approximately 400 to 500 bacterial cells were counted per filter, using the 1000× objective. Microbial cells were fractionated into size classes: small single cells (<2 µm), ultra-small single cells (cell volume <0.04 µm3), large cells (>2 µm), and chain cells (cell length-width ratio ≥10). The dimensions of at least 30 small cells in each sample were measured, and the average value (ABV–average bacterial cell volume, µm3) was calculated. All cells >2 µm were counted, directly, and their biovolume was calculated, according to shape. The wet bacterial biomass was converted to carbon biomass (C, fg·cell−1), using the following equation: C = 120 × V0.72, where V–average cell volume, µm3 [39]. Virus-like particles were enumerated, after staining the samples with SYBR Green I fluorochrome (Molecular Probes, Eugene, OR, USA) [40]. Depending on the number of viruses, between 0.5–1.0 mL of seawater was filtered through 0.02-μm pore-size Anodisc aluminum oxide membrane filters (Whatman). Two filters were analyzed for each water sample, and a minimum of 400 viruses were counted, on each filter.
Contact rates (R) between viruses and bacteria at each sample were determined from the equation: R = (Sh × 2π × w × Dv) × V × P [41], where S is the Sherwood number (dimensionless), w is the diameter of a marine bacterium (0.45 × 10−4 cm; Ref. [42]), D is the diffusivity of viruses (3.456 × 103 cm2·d−1), and V and P are the viral and bacterial abundances, respectively. Additionally, we calculated the ratio of viral to bacterial abundance (VBR).

2.3. Statistical Analyses

Differences in environmental parameters (temperature, salinity, nutrient and chlorophyll a concentrations, zooplankton abundance, and biomass) between different water masses and depth layers were tested using one-way ANOVA or the Kruskal–Wallis test and a post hoc Tukey test (p < 0.05). Assumptions of ANOVA were checked using a Kolmogorov–Smirnov test for normality and Levene’s test for homogeneity of variances. The Kruskal–Wallis test is the non-parametric analog of a one-way ANOVA, which does not make assumptions about normality. It is performed on the ranks of the measurement observations [43]. Non-parametric multivariate techniques were applied to microbial data using the PRIMER v5 software package [44]. Hierarchical cluster analysis was performed on the bacterioplankton and virioplankton abundance data. Prior to the analyses, the data were lg(x + 1)-transformed, to ensure the normality and homogeneity of variances. Similarities between the samples were assessed, using the Bray–Curtis index [45]. The relative contribution of each microbial cell size class and total virioplankton number to the dissimilarities between stations was assessed, by applying a one-way similarity percentage analysis (SIMPER; [44]). To determine the spatial variations, non-metric multidimensional scaling (nMDS) ordinations were computed and Analysis of Similarities (ANOSIM) was used, to detect significant (p = 0.001) differences between water masses and depth layers [45], with respect to microbial abundance. ANOSIM is sensitive to differences in multivariate dispersion. Therefore, to test the hypothesis of equality in group dispersions, we applied a multivariate analog to Leneve’s test. Homogeneity of dispersion among groups was calculated, as an average distance (±SE) of group members (samples) from the group’s centroid [46]. Comparisons of bacterial and viral abundances (ANOSIM), initially, showed no significant differences between the 0 m, 10 m, and 25 m layers and between the 50–100 m layers, although significant differences were found, between the 0–25 m and 50–100 m layers. Therefore, environmental and biological data were grouped into three groups, corresponding to the surface (0–25 m), intermediate (25–100 m), and bottom (192–365 m) layers.
To study the relationship between environmental variables (water temperature, salinity, sampling depth, nutrient concentrations, chlorophyll a content, zooplankton abundance, and biomass) and microbial data (bacterial abundance and biomass, viral abundance, ABV, R, and VBR), constrained ordination procedures were performed in the CANOCO for Windows v4.5 software [47]. The linear model (redundancy analysis, RDA) was chosen, after test runs of detrended correspondence analysis, which showed that the gradient length of the first axis was 0.235, suggesting that microbial data had linear microbial data–environment responses [47]. RDA represents a linear model, i.e., by linear combinations of environmental variables. Before including the selected environmental variables in the RDA, we tested our data for collinearity and, then, removed the variables that had a variance inflation factor higher than 10 (VIF > 10).
Ordination techniques represent multivariate methods that arrange sites along axes on the basis of biological data (species composition, abundance, biomass, etc.). The aim of this approach is to summarize multivariate data, in a convenient way, in scatter diagrams [47]. RDA is a canonical version of principal component analysis (PCA), of dependent variables Y (microbial parameters) constrained by independent variables X (environmental factors) [48,49]. The goal of RDA is to find the combination of independent variables that best explains the variation (or dispersion) of dependent variables. As in regression, each biological variable constitutes a response variable but, in ordination, these response variables are analyzed, simultaneously [47]. Ordination techniques are, widely, used in various studies to assess possible relationships between environmental factors (explanatory variables) and biological data [50]. In most studies, the RDA is used to reveal variations in biological communities, across a range of different environmental conditions. However, this method may be applied to the integrated parameters or several “traits” (e.g., total abundance, total biomass, and some others) [47]. In all cases, the primary data set contains records on a collection of observations–samples (sampling units) [47]. Each sample comprises measures for multiple species or, less often, the other kinds of descriptors (integrated community parameters) [47]. The primary data can be represented, as a rectangular matrix, where the rows refer to individual samples, and the columns indicate individual variables [47]. After log-transformation, the environmental and biological variables were, automatically, centered and standardized by the CANOCO program. Standardizing variables is the centering of the variables, around a mean of zero. Thus, the variables were converted to a z-score, so that the means were equal to 0 and the variance to 1 [49,50].
The forward procedure of non-collinear environmental variables in CANOCO was used, with a significance level of 0.05. The statistical significance of each environmental variable was evaluated, by Monte Carlo permutation tests [45]. RDA was supplemented, with the Pearson correlation analysis, to find possible significant correlations between biotic and environmental variables. The Benjamini–Hochberg procedure was used to correct p-values, for multiple testing. In each case, a false discovery rate was set at 0.1. Means are presented with standard errors (±SE).

3. Results

3.1. Environmental Parameters

Water temperature varied in a wide range (−1.65 °C to 2.51 °C) and the mean values tended to decrease, from the southeastern to the northwestern part of the study area (Table 1, Figures S1 and S2). The upper 10 m layer, at all stations, was characterized by positive temperatures, while below 50–100 m there were colder waters (<0 °C). The bottom layers were occupied by low temperature waters (−0.9 °C to −1.3 °C) (Figure 2). Only small changes in salinity were recorded over the study period. Salinity ranged from 33.49 to 34.84, with the lowest values in the surface (0–10 m) layers (Figure 2). There was an increase in salinity with depth. There was a weak pycnocline located at 30–60 m. Two main types of water masses were recorded in the study area (Figure S1). The Barents Sea Water (BSW), a water mass formed mainly from AW, as a result of the heat loss, was identified at two stations (St. 7 and 8) located near Novaya Zemlya. The Arctic Water (ArW) was present at the rest stations. There were significant differences in water temperature and salinity between the water masses and depth layers (one-way ANOVA or Kruskal–Wallis test, p < 0.05), except for the bottom layers (Table 2). The sea was ice-free during the study period.
Nutrient concentrations decreased from the southern to the northern part of the study area (Table 2). This trend was, particularly, evident for phosphate and nitrate. Phosphate concentrations ranged from <0.01 µM to 0.89 µM, while nitrate concentrations varied, from 0.1 µM to 18.65 µM. Vertical distributions of these nutrients showed an increase in the deepwater layers (Table 2). Mean oxygen concentrations demonstrated low horizontal variability (6.81–8.20 mL·L−1), although there was a clear vertical trend, with the lowest values being recorded below 50–100 m (Table 2). One-way ANOVA revealed significant differences (p < 0.001) in phosphate and oxygen concentrations, between water masses and depth layers, while nitrate concentrations were similar in both water masses but differed regarding sampling layers (p < 0.001). Phosphate and nitrate concentrations were significantly positively correlated with salinity and negatively correlated with water temperature (Pearson’s correlation, Benjamini–Hochberg p < 0.1). In contrast, oxygen concentrations tended to increase with temperature and decrease with salinity (Pearson’s correlation, Benjamini–Hochberg p < 0.1).
Surface chlorophyll a concentrations varied between 0.18 and 0.36 mg·m−3 and were, relatively, homogeneous in the study area (Table 2). No clear relationships (Pearson’s correlation, Benjamini–Hochberg p < 0.1) were detected between temperature, salinity, nutrients, and chlorophyll a concentrations.
Zooplankton density ranged between 752 ind.·m−3 (st. 10) and 1278 ind.·m−3 (st. 13), averaging 1012 ± 87 ind.·m−3 (Table 2, Figure 3). The total zooplankton biomass varied from 0.8 to 11.8 mgC·m−3, with a mean value of 4.9 ± 1.7 mgC·m−3 (Table 2, Figure 3). Both parameters increased along Transects 1 and 2 northwestward. Significant spatial differences were found, regarding the total zooplankton biomass, with the highest estimates in ArW (one-way ANOVA, p < 0.001) (Table 2, Figure 3). Copepods were the most diverse and numerous group, accounting for 94–99% of the total abundance and 40–92% of the total biomass (Figure 3). The zooplankton assemblage of BSW was characterized by a high proportion of small copepods (Microcalanus pygmaeus, Microsetella norvegica, Oithona similis, Triconia borealis, Pseudocalanus spp.), in the total zooplankton biomass (Figure 3), while larger species (Calanus spp. and Metridia longa) amounted to >50% of the total biomass in ArW (Figure 3). The total zooplankton biomass was found to be positively correlated with surface chlorophyll a (Pearson’s correlation, Benjamini–Hochberg p < 0.1).

3.2. Spatial Variations of Microbial Plankton

The bacterial abundance varied from 2.20 × 105 to 10.68 × 105 cells·mL−1 (Figure 4a), with an average of 4.97 ± 0.27 × 105 cells·mL−1. The surface and bottom counts increased from the Novaya Zemlya coast to Franz Joseph Land (Figure 4a). Mean values detected in BSW were comparable with those found in ArW (Table 2), although there were significant differences in the total bacterial abundance between the 0–100 m and bottom layers (Kruskal–Wallis Test, p < 0.001). The abundance showed two peaks at the intermediate layer (50 m) as well as stations 10 and 14 (Figure 4a). The bacterial biomass varied from 2.70 mgC·m−3 to 12.48 mgC·m−3 (mean 5.90 ± 0.33 mgC·m−3) and demonstrated mosaic distribution within the study area, although maximum records were found at stations affected by ArW (Figure 4b). Small and ultra-small bacterial cells were the most numerous in the total bacterioplankton, comprising 98.7–99.7% (Figure 5a,b). Large bacterial cells accounted for 0.1–5.1% in the bacterioplankton and tended to be concentrated in the surface layer (Figure 5c). Chain-bacterial cells were, rarely, found in the samples and their proportion, in the total bacterial density, did not exceed 0.3% (Figure 5d). They reached maximal abundance in the intermediate layer. ABV ranged from 0.020 μm3 to 0.136 μm3 and tended to be slightly higher in ArW, relative to BSW, although the differences were non-significant (one-way ANOVA, p > 0.05) (Table 2).
Viral abundance ranged from 0.86 × 106 particles·mL−1 to 14.98 × 106 particles·mL−1, averaging 6.10 ± 0.45 × 106 particles·mL−1 (Figure 4c). There were no significant differences between BSW and ArW, in the total number of viruses (one-way ANOVA, p > 0.05), but the surface and intermediate layers had greater values (Kruskal–Wallis Test, p < 0.001). The viral abundance was positively correlated with the bacterial abundance and biomass (Pearson’s correlation, Benjamini–Hochberg p < 0.1). VBR was estimated to be 1–28 (11 ± 1) (Table 2). This variable did not differ between water masses and depth layers (one-way ANOVA or Kruskal–Wallis Test, p > 0.05). Contact rates (R) between viruses and bacteria were 0.3–9.6 (3.5 ± 0.3) (Table 2) and significantly decreased with depth (Kruskal–Wallis Test, p < 0.001). Pearson correlation analysis found significant positive correlations of R, with viral and bacterial abundances (Pearson’s correlation, Benjamini–Hochberg p < 0.1), while R was negatively correlated with ABV (Pearson’s correlation, Benjamini–Hochberg p < 0.1).
Multivariate analyses (hierarchical cluster analysis and NMDS ordination), based on the abundance of microbial plankton, revealed a high degree of similarity between the two types of water masses in the region (Figure 6). The ANOSIM showed significant differences in the bacterioplankton and virioplankton abundances, between different sampling layers (global R: 0.195, p = 0.013), although there were no significant differences in the microbial density averaged for water masses (global R: −0.161, p = 0.974). Based on the SIMPER analysis, dissimilarities ranged from 19.9% to 44.4%. The highest dissimilarity occurred for the intermediate and bottom layers. The SIMPER analysis, also, indicated that viral abundance contributed, mostly (92–93%), to the total dissimilarities between sampling layers.

3.3. Environmental Control of Microbial Plankton

Possible relationships, between the microbial plankton and the environmental dataset, were analyzed using RDA analysis and supplemented with the Pearson correlation analysis (Figure 7, Table 3 and Table 4). The first two RDA axes (eigenvalues 0.480 and 0.127) explained 48.0% and 12.7% of the total variability in the microbial parameters, during the study period. The first RDA axis was significantly positively correlated with nitrate content and negatively correlated with water temperature and oxygen concentration (Figure 7).
The second axis was significantly positively correlated with zooplankton biomass and oxygen concentration and negatively correlated with salinity (Figure 7). The forward selection of environmental factors, using a Monte Carlo permutation test (999 permutations), showed that zooplankton biomass, nitrate concentration, and salinity were the factors that contribute significantly to the observed variability in microbial-plankton variations (Figures S2–S11 and Figure 7, Table 3). The three environmental variables together explained 45% of the total variance in microbial parameters (Table 3). The high bacterial abundances and biomasses were associated with low nitrate and phosphate contents as well as high oxygen concentrations (Figure 7, Table 3). The total viral abundances were connected with high temperature and low zooplankton biomass (Figure 7, Table 3 and Table 4). VBR tended to decrease with increasing zooplankton biomass (Figure 7, Table 3 and Table 4). ABV was positively correlated with nitrates (Figure 7, Table 3 and Table 4). High R was associated with low phosphate and nitrate as well as high oxygen concentrations (Figure 7, Table 3 and Table 4).

4. Discussion

The present study aimed to fill the gap in our knowledge of microbial plankton, in a less studied region strongly affected by cold waters from the Arctic Ocean, after the main productive season. Our data is the first report describing spatial patterns in bacterioplankton and virioplankton, in the northeastern Barents Sea. In general, we revealed that microbial abundance was strongly affected by a set of environmental variables (nutrient concentrations, zooplankton carbon biomass, and hydrological conditions).

4.1. Environmental Parameters

A clear warming trend has been detected in the Barents Sea, since the 2000s, and the mean water temperature was found to have risen [1,2,3,4]. In particular, in most regions of the Barents Sea, during the period of 1998–2017, the highest values of the mean annual surface temperature were registered in regions affected by warm AW, while the lowest estimates were in the northeastern Barents Sea [3]. Remote sensing data suggest that a 1.0 °C increase over the 20-year period was evident for the region [3]. We revealed that water temperature in the upper 0–50 m layer was somewhat higher, relative to the mean values recorded in October (1.0–1.2 °C vs. −1–0 °C), based on observations from 1950 to 1998 [51]. Moreover, in temperate and cold years, the region of Franz Joseph Land is covered by ice, but there was no sea ice during our cruise. This finding is in line with a general trend of ice loss, in the Arctic Ocean and adjacent areas [1,3,4]. Surface salinity, recorded in October 2020, was slightly higher, compared to the mean values obtained in 1950–1998 [51], suggesting a lower impact of cold ArW.
Nutrient concentrations demonstrate strong seasonal and spatial fluctuations in the Barents Sea [1,2,52]. Total depletion of phosphate and nitrate may occur in the upper layers during the spring–summer period, due to phytoplankton bloom in the Arctic [1,2,53,54]. Spring bloom leads to a decreased stock of nutrients, and in ArW of the Marginal Ice Zone, the outburst may occur earlier than in AW [54]. For instance, nitrate and phosphate concentrations in ArW (central Barents Sea), during spring bloom, were 2–3 μM and 0.2–0.3 μMm, respectively, while these in AW were 5–6 μM and 0.5–0.6 μM, respectively [55]. In the eastern Barents Sea, nitrate and phosphate concentrations in spring were 0.27–0.30 µM and 1.34–3.81 µM, respectively [53]. In the western Barents Sea, a clear seasonal cycle of nitrate has been documented, with a maximum of 9–10 µM in March and a minimum of ~1 µM in August [56]. In more northern regions of the Barents Sea, a peak of blooming is, typically, recorded in July–August [1]. In ArW, there is a rise in nutrient supply from late summer to winter, when nitrate and phosphate reach their maximum values [1].
In our study, we recorded a clear vertical pattern in the nutrient concentrations, with the richest waters being located in deeper layers. This pattern can be explained by intense vertical mixing processes, in the upper 50 m layer, and convection, as revealed in earlier studies [1,2]. Enhanced bottom concentrations of nutrients seemed to be associated with decreased water exchange, at deepwater sites (e.g., St. Anna Through), relative to intermediate and surface waters. Phytoplankton, also, play an important role in the regulation of overall nutrient stocks in the Arctic. The autumn stage of phytoplankton succession in the high Arctic regions may be specified with decreased abundance and biomass of pelagic microalgae [1,54]. Autumn–winter chlorophyll a concentrations are low and do not exceed 0.1–0.5 mg·m−3 [1,26,54,56]. We, also, revealed low surface chlorophyll a values, which confirmed this general pattern, for phytoplankton density to decrease, after the main productive season in the northern Barents Sea.
Our study found that zooplankton assemblages were dominated by copepods. The zooplankton composition was typical for ArW and BSW with Calanus spp., Pseudocalanus spp., Microcalanus pygmaeus, Oithona similis, and Metridia longa being common taxa [2,28,30,57]. Moreover, these species are among the most frequent and abundant zooplankters in other Arctic regions, including Svalbard waters, the northern Kara Sea, and Greenland waters [1,2,28,29,30,31,58]. We observed an increase in the total mesozooplankton biomass in ArW, compared to BSW, and this pattern was associated with the prevalence of small taxa in BSW. Larger copepods have a higher relative abundance in ArW, and the mean biomass here is 3.7 times higher than in BSW. Zooplankton abundance and biomass estimates in the present study were lower, when compared to the values recorded in ArW in the Kola Transect (33°30′ E), in the autumn–winter periods of 2011–2012 [30], and these differences can be explained by sampling methods and climatic conditions. Our study was focused on the upper 100 m layer, while the data obtained in 2011–2012 considered the whole water layer, plus the period of 2020 was colder, relative to 2011–2012 [3,46]. Local variability of hydrological conditions and water circulation might, also, be responsible for differences in the zooplankton density [1,2,29,30]. We, also, established a positive correlation between the zooplankton biomass and surface chlorophyll a. This result is expected because most species in our study were found to be herbivorous or omnivorous, and phytoplankton, apparently, played an important role in their feeding in the Barents Sea and other Arctic regions [1,2,59].

4.2. Spatial Variation of Microbial Plankton

In the high Arctic, prolonged periods of low water temperature, ice coverage, and darkness in the autumn–winter season lead to strong seasonal differences in bacterial abundance and biomass as well as concentrations of viruses [15,18,19,60]. In late spring, ice melting and increasing light irradiance promote phytoplankton bloom, which causes a subsequent rise in the abundance of Arctic microbial plankton.
Many studies have reported marked spatial and seasonal variations, in microbial characteristics in the Barents Sea. Spring and summer patterns of bacterioplankton have been investigated in the Marginal Ice Zone and in the northern Barents Sea [21,22,23,24,25]. The total abundance of bacteria was recorded to be 0.41 × 106–4.1 × 106 cells·mL−1, during summer in the central Barents Sea, and these counts were lower than the estimates reported for the northern Barents Sea, during the spring bloom [22] and summer periods [23]. Other studies found that the number of marine bacteria tended to be higher in AW compared to ArW in the late summer period [24,25]. There is only one study that reported autumn values of bacterioplankton abundance and biomass [26]. The authors established that the average bacterioplankton abundance/biomass in the upper 50 m layer varied between 0.4 × 106 cells·mL−1/66 mg·m−3 (ArW) and 0.6 × 106 cells·mL−1/41 mg·m−3 (AW), in mid-November 2013 [26]. Our study revealed lower values of the total bacterial biomass in the upper 50 m layer, while the bacterial abundance was comparable with previous autumn estimates [26]. In the northern Kara Sea (St. Anna Through region), the autumn abundance of marine bacteria in the surface layer was reported to be much lower [61,62], in comparison to our data, but the total bacterial biomass was comparable [61]. Spatial and seasonal differences may be explained by climatic variations, local hydrology, nutrient dynamics, and abundance of phytoplankton and zooplankton. In particular, previous studies conducted in the Barents Sea were performed during a warmer period (2010–2016) [24,25,26]. Other studies were focused on the spring and summer seasons, which are considered to be the most productive periods in the year, when biomass of phytoplankton was the highest [21,22,23]. The lower abundance of autumn bacterioplankton, in the northern Kara Sea, relative to our study, may be explained by oligotrophic conditions [62].
Virioplankton abundance in the Barents Sea, also, demonstrated seasonal fluctuations, with the highest values occurring in the summer period. In the central Barents Sea, the number of virus-like particles ranged between 0.66 × 106 and 8.6 × 106 particles·mL−1 [21]. Similar counts were reported for AW, while in ArW, virioplankton abundance tended to be lower [24,25]. Autumn abundances of viruses in the upper 50 m layer in the Barents Sea [26] were similar to those found in the summer period [24,25]. Our surface values were in good accordance, with previously detected summer and autumn estimations. In the northern Kara Sea, the mean abundance of virioplankton in autumn [63] was much lower compared to our data, which may, also, be explained by oligotrophic conditions in the region. In other Arctic regions, namely, the Laptev Sea and the East Siberian Sea, autumn virioplankton counts were similar to the mean values recorded in our study [64,65].
We detected a positive significant correlation between the abundance of viruses and bacterial stock, and this pattern appears to be common in the Arctic marine environment [20]. The variability of viral abundance is, largely, influenced by bacteria, pelagic microalgae, protists, and zooplankton [1,12]. In the Arctic, phytoplankton assemblages, during the spring-bloom succession, actively consume nutrients to grow and develop, and they produce organic matter, which is utilized by marine bacteria, promoting subsequent viral infection and viral abundance in seawater [20,66]. This pattern can be detected by analyzing VBR. The maximum VBR levels are typical, for regions with higher bacterial abundance and biomass. Previous investigations in the Barents Sea have documented a tendency of summer VBR to be greater, in regions with an enhanced number of marine bacteria [24,25]. During autumn and winter periods, VBR decreased, due to a drop in the total microbial stock in the Barents and Kara seas [26,63]. Our VBR estimations were higher, compared to summer values [24,25], suggesting a higher impact of viruses on bacterioplankton, after the main productive season. We revealed that contact rates (R) between viruses and bacteria significantly decreased with depth, while an inverse relationship was detected in the case of AVB. Our finding suggests that viral particles were associated, mainly, with smaller bacterial cells occurring in the upper layers, although it might be expected that larger bacteria would be preferable for viruses, due to a larger surface for attachment. However, small bacterial cells seem, likely, to be more accessible, owing to their greater abundance. Larger cells were concentrated in the bottom layers and had lower densities.
Results of the present study showed patchiness in the horizontal distribution of bacteria and viruses in the Barents Sea. This is a well-documented phenomenon, reflecting the discontinuous distribution of individual organisms in aquatic environments [59]. Aggregation of plankton may be connected with frontal zones, local eddies, environmental gradients, food stocks, predator abundance, and biotic interactions [59]. Hierarchical cluster analysis indicated a low degree of dissimilarity between sampling stations, and NMDS ordination based on microbial parameters showed a high stress value. Moreover, ANOSIM indicated non-significant differences in bacterioplankton and virioplankton, between BSW and ArW. Therefore, our hypothesis regarding the association of microbial plankton to water masses must be rejected. Nevertheless, we revealed a clear vertical pattern in the distribution of marine microbes, with minimum parameters being registered in the bottom layer. AVB, also, increased with depth, and larger bacteria were accumulated near the seafloor. SIMPER indicated the importance of marine viruses, to distinguish depth layer by microbial parameters. Greater abundances of marine bacteria and viruses, in the upper and intermediate layers, are a common feature reported for various Arctic regions [60,61,62,63,64,65], including the Barents Sea [21,22,23,24,25,26]. This pattern can be explained by the nutrient distribution and local hydrological properties.

4.3. Environmental Control of Microbial Plankton

The RDA revealed that a set of environmental variables (temperature, salinity, nutrients, chlorophyll a, and zooplankton stock) accounted for 61% of the total fluctuation in microbial plankton. Three parameters were found to be significantly associated with bacterioplankton and virioplankton. The most important factor affecting the distribution of microbial plankton was nitrate concentration. Negative correlations between nutrient supply and microbial abundance are considered to be, indirectly, connected with phytoplankton. Chlorophyll a was found to be one of the main drivers, controlling both abundance and biomass of marine microbes in the Arctic, during the spring and summer seasons [21,22,62,64,65]. In our study, chlorophyll a values were low, suggesting that phytoplankton biomass was lower than in the most productive season. Nevertheless, microalgae seem to be important organisms, utilizing inorganic nutrients in the upper layer. On the other hand, there were great concentrations of nitrate and phosphate in the bottom layer, where bacterial abundance was significantly lower. This, also, might explain negative correlations between inorganic nutrients and microbial plankton. The enhanced supply of organic matter increases bacterial growth and leads to an increase in viral abundance in the marine environment [66,67]. Phytoplankton are the main source of organic carbon for microbial plankton, in the spring and summer periods, while other organisms may be important after the main productive period [11]. Microzooplankton and copepods can, rapidly, recycle and regenerate both dissolved organic and inorganic nutrients, through feeding, excretion, and releasing of fecal pellets in autumn [68]. Moreover, zooplankton, rather than phytoplankton, appear to be the key autochthonous source of bacterial substrates.
In contrast, virioplankton abundance and VBR negatively correlated with zooplankton biomass. This result is in contrast to other studies reporting an increase in viral abundance with carbon biomass, in the spring and summer periods [60,64]. It is more likely that these differences might be related to seasonal fluctuations in the distribution of plankton animals and their vertical migrations as well as the dominance of small-sized bacterioplankton in our study. Viruses are associated with suspended and particulate particles as well as small heterotrophic plankton [36]. Larger zooplankters, often, consume microzooplankton in the autumn and winter seasons [59,68], thus reducing the number of potential hosts for viruses. The observed pattern may, also, be explained by a direct influence of the smaller zooplankton, on potential hosts for marine viruses. Nauplii and omnivorous copepods might ingest detritus and particulate organic matter [68], a potential surface for marine viruses, and this might result in reduced virioplankton abundance. An alternative explanation is food selectivity of the larger zooplankton, which can consume different-sized fractions of microzooplankton [68]. They might prefer to ingest larger cells of protists, while the remaining nanoplankton and picoplankton seemed to be responsible for grazing marine bacteria and virus-like particles. Additionally, some modeling studies [69,70] have shown that zooplankton organisms are able to avoid virally infected phytoplankton and this, therefore, might be reflected as a negative correlation between the total zooplankton abundance/biomass and viral density. Our study and previous investigations dealing with autumn and winter Arctic zooplankton [29,30] highlighted the dominance of small copepods in the total zooplankton abundance, most of which are found to be omnivorous, preferring to ingest protozooplankton in the periods of low phytoplankton availability [68].
Salinity, together with water temperature, explained part of the variation in microbial abundance in the Barents Sea. Bacterioplankton and virioplankton abundances demonstrated a tendency to be lower with salinity, and an opposite pattern was revealed, in the case of water temperature. This result can be interpreted, considering the vertical distribution of microbes. Greater densities were located in warmer surface waters, while bottom high-saline and colder waters showed lower microbial counts. The activity of marine bacteria is known to be temperature-dependent [71] and, therefore, their growth would be slower at lower temperatures, near the bottom. Other investigations documented that temperature and salinity play a minor role in determining spatial variations, in the microbial abundance and biomass in the Arctic regions [9,22,72]. Thus, it is the concentration of organic and inorganic nutrients that are the most responsible factors affecting microbial plankton distribution in the Barents Sea, but sources of organic carbon, obviously, differ seasonally.

5. Conclusions

Microbial communities are, strongly, involved in biogeochemical cycles and represent a major link in the marine food webs and carbon fluxes. Plankton assemblages are considered good indicators of environmental impacts. The Arctic Ocean is the region, where clear climatic fluctuations have been reported. Studying plankton in the Arctic provides valuable data, to predict overall biotic responses to climatic changes. The present study was focused on a less-investigated region of the Arctic Ocean—the northeastern Barents Sea. Bacterioplankton abundance and biomass as well as concentration of virus-like particles demonstrated mosaic horizontal distribution in the autumn period. There were no clear associations of the microbial variables, with two main water masses (colder Arctic Water and warmer Barents Sea Water). A clear vertical pattern in microbial abundance was evident, with the highest concentrations in the upper layers. Small-sized cells were the most dominant group of the bacterioplankton, and average bacterial-cell volume tended to increase, near the seafloor. Bacterial abundance was lower than in other seasons and regions of the Barents Sea, while viral abundance was comparable. Bacterioplankton abundance was positively correlated with concentrations of viruses, indicating close relationships between these groups of microbes. Nutrients and organic carbon of zooplankton, together with hydrological variables, explained most of the spatial variability in microbial parameters.
Our study expands current knowledge on the structure of pelagic ecosystems in the Arctic and provides baseline information on marine microbes in the autumn season. This investigation may be used in future research, focusing on the ecology of Arctic plankton assemblages, during the period of global environmental changes.

Supplementary Materials

The following are available online at https://www.mdpi.com/article/10.3390/biology11060845/s1, Figure S1: T-S diagram of different water masses in the northeastern Barents Sea, in autumn 2020. BSW—Barents Sea Water, ArW—Arctic Water. Figure S2: Relationship between bacterial abundance (NB) and phosphate concentration in the northeastern Barents Sea, in autumn 2020. R2—determination coefficient. N = 30. Figure S3: Relationship between bacterial abundance (NB) and nitrate concentration in the northeastern Barents Sea, in autumn 2020. R2—determination coefficient. N = 30. Figure S4: Relationship between viral abundance (NV) and water temperature in the northeastern Barents Sea, in autumn 2020. R2—determination coefficient. N = 30. Figure S5: Relationship between viral abundance (NV) and nitrate concentration in the northeastern Barents Sea, in autumn 2020. R2—determination coefficient. N = 30. Figure S6: Relationship between viral abundance (NV) and zooplankton biomass in the northeastern Barents Sea, in autumn 2020. R2—determination coefficient. N = 30. Figure S7: Relationship between VBR (the ratio of viral to bacterial abundance) and zooplankton biomass in the northeastern Barents Sea, in autumn 2020. R2—determination coefficient. N = 30. Figure S8: Relationship between ABV (average bacterial cell volume) and nitrate concentration in the northeastern Barents Sea, in autumn 2020. R2—determination coefficient. N = 30. Figure S9: Relationship between R (contact rate between viruses and bacteria) and water temperature in the northeastern Barents Sea, in autumn 2020. R2—determination coefficient. N = 30. Figure S10: Relationship between R (contact rate between viruses and bacteria) and phosphate concentration in the northeastern Barents Sea, in autumn 2020. R2—determination coefficient. N = 30. Figure S11: Relationship between R (contact rate between viruses and bacteria) and nitrate concentration in the northeastern Barents Sea, in autumn 2020. R2—determination coefficient. N = 30.

Author Contributions

Conceptualization, V.G.D. and M.P.V.; methodology, V.G.D., M.P.V., A.V.V., T.M.M., T.G.I. and V.V.V.; software, V.G.D.; validation, V.G.D.; formal analysis, V.G.D. and M.P.V.; investigation, V.G.D., M.P.V., A.V.V., T.M.M. and T.G.I.; resources, V.G.D.; data curation, V.G.D.; writing—original draft preparation, V.G.D.; writing—review and editing, V.G.D., M.P.V., A.V.V. and T.M.M.; visualization, V.G.D., M.P.V., A.V.V. and T.M.M.; supervision, V.G.D.; project administration, V.G.D. All authors have read and agreed to the published version of the manuscript.

Funding

This study was funded by the Ministry of Science and Higher Education of the Russian Federation.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data are available on request, from the corresponding author.

Acknowledgments

We are grateful to the crew and staff of the R/V Dalnie Zelentsy, for their support during sampling. We thank the academic editor and the three anonymous reviewers, for the constructive and helpful comments.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Sakshaug, E.; Johnsen, G.; Kovacs, K. (Eds.) Ecosystem Barents Sea; Tapir Academic Press: Trondheim, Norway, 2009. [Google Scholar]
  2. Jakobsen, T.; Ozhigin, V.K. (Eds.) The Barents Sea: Ecosystem, Resources, Management: Half a Century of Russian-Norwegian Cooperation; Tapir Academic Press: Trondheim, Norway, 2011. [Google Scholar]
  3. ICES. Working Group on the Integrated Assessments of the Barents Sea (WGIBAR). ICES Sci. Rep. 2021, 3, 77, 1–236. [Google Scholar]
  4. Ozhigin, V.; Ivshin, V.; Trofimov, A.; Karsakov, A.L.; Antsiferov, M. The Barents Sea Water: Structure, Circulation, Variability; PINRO Press: Murmansk, Russia, 2016. [Google Scholar]
  5. Polyakov, I.V.; Alkire, M.B.; Bluhm, B.A.; Brown, K.A.; Carmack, E.C.; Chierici, M.; Danielson, S.L.; Ellingsen, I.; Ershova, E.A.; Gårdfeldt, K.; et al. Borealization of the Arctic Ocean in response to anomalous advection from sub-arctic seas. Front. Mar. Sci. 2020, 7, 491. [Google Scholar] [CrossRef]
  6. Ardyna, M.; Arrigo, K. Phytoplankton dynamics in a changing Arctic Ocean. Nat. Clim. Change 2020, 10, 892–903. [Google Scholar] [CrossRef]
  7. Dvoretsky, V.G.; Dvoretsky, A.G. Epiplankton in the Barents Sea: Summer variations of mesozooplankton biomass, community structure and diversity. Cont. Shelf Res. 2013, 52, 1–11. [Google Scholar] [CrossRef]
  8. Kirchman, D.L.; Morán, X.A.G.; Ducklow, H. Microbial growth in the polar oceans—Role of temperature and potential impact of climate change. Nat. Rev. Microbiol. 2009, 7, 451–459. [Google Scholar] [CrossRef] [PubMed]
  9. Boras, J.A.; Sala, M.M.; Arrieta, J.M.; Sá, E.L.; Felipe, J.; Agustí, S.; Duarte, C.M.; Vaqué, D. Effect of ice melting on bacterial carbon fluxes channeled by viruses and protists in the Arctic Ocean. Polar Biol. 2010, 33, 1695–1707. [Google Scholar] [CrossRef]
  10. Ducklow, H.W. Minireview: The bacterial content of the oceanic euphotic zone. FEMS Microbiol. Ecol. 1999, 30, 1–10. [Google Scholar] [CrossRef]
  11. Fuhrman, J.A. Microbial community structure and its functional implications. Nature 2009, 459, 193–199. [Google Scholar] [CrossRef]
  12. Wommack, K.E.; Colwell, R.R. Virioplankton: Viruses in aquatic ecosystems. Microbiol. Mol. Biol. Rev. 2000, 64, 69–114. [Google Scholar] [CrossRef] [Green Version]
  13. Suttle, C.A. Marine viruses—Major players in the global ecosystem. Nat. Rev. Microbiol. 2007, 5, 801–812. [Google Scholar] [CrossRef]
  14. Wilhelm, S.; Suttle, A.C. Viruses and Nutrient Cycles in the Sea. BioScience 1999, 49, 781–788. [Google Scholar] [CrossRef] [Green Version]
  15. Steward, G.F.; Fandino, L.B.; Hollibaugh, J.T.; Whitledge, T.E.; Azam, F. Microbial biomass and viral infections of heterotrophic prokaryotes in the sub-surface layer of the central Arctic Ocean. Deep-Sea Res. I 2007, 54, 1744–1757. [Google Scholar] [CrossRef]
  16. Steward, G.F.; Smith, D.C.; Azam, F. Abundance and production of bacteria and viruses in the Bering and Chukchi Seas. Mar. Ecol. Prog. Ser. 1996, 131, 287–300. [Google Scholar] [CrossRef]
  17. Sherr, B.F.; Sherr, E.B. Community respiration/production and bacterial activity in the upper water column of the central Arctic Ocean. Deep-Sea Res. I 2003, 50, 529–542. [Google Scholar] [CrossRef]
  18. Sherr, E.B.; Sherr, B.F.; Wheeler, P.A.; Thompson, K. Temporal and spatial variation in stocks of autotrophic and heterotrophic microbes in the upper water column of the central Arctic Ocean. Deep-Sea Res. I 2003, 50, 557–571. [Google Scholar] [CrossRef]
  19. Rokkan Iversen, K.; Seuthe, L. Seasonal microbial processes in a high-latitude fjord (Kongsfjorden, Svalbard): I. Heterotrophic bacteria, picoplankton and nanoflagellates. Polar Biol. 2011, 34, 731–749. [Google Scholar] [CrossRef] [Green Version]
  20. Yau, S.; Seth-Pasricha, M. Viruses of Polar Aquatic Environments. Viruses 2019, 11, 189. [Google Scholar] [CrossRef] [Green Version]
  21. Howard-Jones, M.H.; Ballard, V.D.; Allen, A.E.; Frischer, M.E.; Verity, P.G. Distribution of bacterial biomass and activity in the marginal ise zone of the central Barents Sea during summer. J. Mar. Syst. 2002, 38, 77–91. [Google Scholar] [CrossRef]
  22. Sturluson, M.; Nielsen, T.G.; Wassmann, P. Bacterial abundance, biomass and production during spring blooms in the northern Barents Sea. Deep-Sea Res. II 2008, 55, 2186–2198. [Google Scholar] [CrossRef]
  23. Tammert, H.; Olli, K.; Sturluson, M.; Hodal, H. Bacterial biomass and activity in the marginal ice zone of the northern Barents Sea. Deep-Sea Res. II 2008, 55, 2199–2209. [Google Scholar] [CrossRef]
  24. Venger, M.P.; Shirikolobova, T.I.; Makarevich, P.R.; Vodopyanova, V.V. Viruses in the pelagic zone of the Barents Sea. Dokl. Biol. Sci. 2012, 446, 345–349. [Google Scholar] [CrossRef]
  25. Venger, M.P.; Kopylov, A.I.; Zabotkina, E.A.; Makarevich, P.R. The influence of viruses on bacterioplankton of the offshore and coastal parts of the Barents Sea. Rus. J. Mar. Biol. 2016, 42, 26–35. [Google Scholar] [CrossRef]
  26. Shirikolobova, T.I.; Zhichkin, A.P.; Venger, M.P.; Vodopyanova, V.V.; Moiseev, D.V. Bacteria and viruses of the ice-free aquatic area of the Barents sea at the beginning of polar night. Dokl. Biol. Sci. 2016, 469, 182–186. [Google Scholar]
  27. Makarevich, P.R. (Ed.) Cruise Report of the Complex Expedition to the R/V Dalnie Zelentsy in the Barents Sea (October 01–October 15, 2020); MMBI: Murmansk, Russia, 2020. [Google Scholar]
  28. Dvoretsky, V.G.; Dvoretsky, A.G. Copepod communities off Franz Josef Land (northern Barents Sea) in late summer of 2006 and 2007. Polar Biol. 2011, 34, 1231–1238. [Google Scholar] [CrossRef]
  29. Dvoretsky, V.G.; Dvoretsky, A.G. Arctic marine mesozooplankton at the beginning of the polar night: A case study for southern and south-western Svalbard waters. Polar Biol. 2020, 43, 71–79. [Google Scholar] [CrossRef]
  30. Dvoretsky, V.G.; Dvoretsky, A.G. Mesozooplankton in the Kola Transect (Barents Sea): Autumn and winter structure. J. Sea Res. 2018, 142, 18–22. [Google Scholar] [CrossRef]
  31. Dvoretsky, V.G.; Dvoretsky, A.G. Coastal mesozooplankton assemblages during spring bloom in the eastern Barents Sea. Biology 2022, 11, 204. [Google Scholar] [CrossRef]
  32. Parsons, T.R.; Maita, Y.; Lalli, C.M. A Manual of Chemical and Biological Methods for Sea Water Analysis; Pergamon Press: New York, NY, USA, 1992. [Google Scholar]
  33. Aminot, A.; Rey, F. Standard Procedure for the Determination of Chlorophyll A by Spectroscopic Methods; International Council for the Exploration of the Sea: Copenhagen, Denmark, 2000. [Google Scholar]
  34. Chislenko, L.L. Nomogrammes to Determine Weights of Aquatic Organisms Based on the Size and Form of Their Bodies (Marine Mesobenthos and Plankton); Nauka Press: Leningrad, USSR, 1968. (In Russian) [Google Scholar]
  35. Hanssen, H. Mesozooplankton of the Laptev Sea and the adjacent eastern Nansen Basin—Distribution and community structure in late summer. Ber. Polarforsch. 1997, 229, 1–131. [Google Scholar]
  36. Richter, C. Regional and seasonal variability in the vertical distribution of mesozooplankton in the Greenland Sea. Ber. Polarforsch. 1994, 154, 1–90. [Google Scholar]
  37. Harris, R.P.; Wiebe, P.H.; Lenz, J.; Skjoldal, H.R.; Huntley, M. (Eds.) ICES Zooplankton Methodology Manual; Academic Press: San Diego, CA, USA, 2000. [Google Scholar]
  38. Porter, K.G.; Feig, Y.S. The use of DAPI for identifying and counting aquatic microflora. Limnol. Oceanogr. 1980, 25, 943–948. [Google Scholar] [CrossRef]
  39. Norland, S. The Relationships between Biomass and Volume of Bacteria. In Handbook of Methods in Aquatic Microbial Ecology; Kemp, P., Sherr, B., Sherr, E., Cole, J., Eds.; Lewis Publishing: Boca Raton, FL, USA, 1993; pp. 303–308. [Google Scholar]
  40. Noble, R.T.; Fuhrman, J.A. Virus decay and its causes in coastal waters. Appl. Environ. Microbiol. 1997, 63, 77–83. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  41. Murray, A.G.; Jackson, G.A. Viral dynamics: A model of the effects of size, shape, motion and abundance of single-celled planktonic organisms and other particles. Mar. Ecol. Prog. Ser. 1992, 89, 103–116. [Google Scholar] [CrossRef]
  42. Lee, S.; Fuhrman, J.A. Relationships between biovolume and biomass of naturally derived marine bacterioplankton. Appl. Environ. Microbiol. 1987, 53, 1298–1303. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  43. Kruskal, W.; Wallis, W. Use of ranks in one-criterion variance analysis. J. Am. Stat. Assoc. 1952, 47, 583–621. [Google Scholar] [CrossRef]
  44. Clarke, K.R.; Gorley, R.N. PRIMER v5: User Manual/Tutorial; PRIMER-E Ltd.: Plymouth, UK, 2001. [Google Scholar]
  45. Clarke, K.R.; Warwick, R.M. Changes in Marine Communities: An Approach to Statistical Analysis and Interpretation; Plymouth Marine Laboratory, NERC: Plymouth, UK, 1994. [Google Scholar]
  46. Anderson, M.J. Distance-based tests for homogeneity of multivariate dispersions. Biometrics 2006, 62, 245–253. [Google Scholar] [CrossRef] [PubMed]
  47. Ter Braak, C.J.F.; Smilauer, P. CANOCO Reference Manualand CanoDraw for Windows User’s Guide: Software for Canonical Community Ordination (version 4.5); Microcomputer Power: Ithaca, NY, USA, 2002. [Google Scholar]
  48. Rao, C.R. The use and interpretation of principal component analysis in applied research. Sankhy Na Ser. A 1964, 26, 329–358. [Google Scholar]
  49. Legendre, P.; Legendre, L. Numerical Ecology; Elsevier Science: Amsterdam, The Netherlands, 1998. [Google Scholar]
  50. Legendre, P.; Anderson, M.J. Distance-based redundancy analysis: Testing multispecies responses in multifactorial ecological experiments. Ecol. Monogr. 1999, 69, 1–24. [Google Scholar] [CrossRef]
  51. Matishov, G.; Zuyev, A.; Golubev, V.; Adrov, N.; Timofeev, S.; Karamusko, O.; Pavlova, L.; Fadyakin, O.; Buzan, A.; Braunstein, A.; et al. Climatic Atlas of the Arctic Seas 2004: Part I. Database of the Barents, Kara, Laptev, and White Seas—Oceanography and Marine Biology; NOAA Atlas NESDIS 58; U.S. Government Printing Office: Washington, DC, USA, 2004.
  52. Olsen, A.; Johannessen, T.; Rey, F. On the nature of the factors that control spring bloom development at the entrance to the Barents Sea and their interannual variability. Sarsia 2003, 88, 379–393. [Google Scholar] [CrossRef] [Green Version]
  53. Makarevich, P.R.; Vodopianova, V.V.; Bulavina, A.S.; Vashchenko, P.S.; Ishkulova, T.G. Features of the distribution of chlorophyll-a concentration along the western coast of the Novaya Zemlya Archipelago in spring. Water 2021, 13, 3648. [Google Scholar] [CrossRef]
  54. Makarevich, P.; Druzhkova, E.; Larionov, V. Primary producers of the Barents Sea. In Diversity of Ecosystems; Mahamane, A., Ed.; InTech: Rijeka, Croatia, 2012; pp. 367–392. [Google Scholar]
  55. Wassmann, P.; Ratkova, T.; Andreassen, I.; Vernet, M.; Pedersen, G.; Rey, F. Spring bloom development in the Marginal Ice Zone and the Central Barents Sea. Mar. Ecol. 1999, 20, 321–346. [Google Scholar] [CrossRef] [Green Version]
  56. Dalpadado, P.; Arrigo, K.R.; van Dijken, G.L.; Skjoldal, H.R.; Bagøien, E.; Dolgov, A.V.; Prokopchuk, I.P.; Sperfeld, E. Climate effects on temporal and spatial dynamics of phytoplankton and zooplankton in the Barents Sea. Progr. Oceanogr. 2020, 182, 102320. [Google Scholar] [CrossRef]
  57. Dvoretsky, V.G.; Dvoretsky, A.G. Checklist of fauna found in zooplankton samples from the Barents Sea. Polar Biol. 2010, 33, 991–1005. [Google Scholar] [CrossRef]
  58. Dvoretsky, V.G.; Dvoretsky, A.G. Regional differences of mesozooplankton communities in the Kara Sea. Cont. Shelf Res. 2015, 105, 26–41. [Google Scholar] [CrossRef]
  59. Raymont, J.E.G. Plankton and Productivity in the Oceans, 2nd ed.; Pergamon Press: Southampton, UK, 1983; Volume 2. [Google Scholar]
  60. Payet, J.P.; Suttle, C.A. Physical and biological correlates of virus dynamics in the southern Beaufort Sea and Amundsen Gulf. J. Mar. Syst. 2008, 74, 933–945. [Google Scholar] [CrossRef]
  61. Sazhin, A.F.; Romanova, N.D.; Mosharov, S.A. Bacterial and primary production in the pelagic zone of the Kara Sea. Oceanology 2010, 50, 759–765. [Google Scholar] [CrossRef]
  62. Mosharova, I.V.; Mosharov, S.A.; Ilinskiy, V.V. Distribution of bacterioplankton with active metabolism in waters of the St. Anna Trough, Kara Sea, in autumn 2011. Oceanology 2017, 57, 114–121. [Google Scholar] [CrossRef]
  63. Kopylov, A.I.; Sazhin, A.F.; Zabotkina, E.A.; Romanova, N.D. Virioplankton in the Kara Sea: The Impact of Viruses on Mortality of Heterotrophic Bacteria. Oceanology 2015, 55, 620–631. [Google Scholar] [CrossRef]
  64. Kopylov, A.I.; Zabotkina, E.A.; Romanenko, A.V.; Kosolapov, D.B.; Sazhin, A.F. Viruses in the water column and the sediment of the eastern part of the Laptev Sea. Estuar. Coast. Shelf Sci. 2020, 242, 106836. [Google Scholar] [CrossRef]
  65. Kopylov, A.I.; Kosolapov, D.B.; Zabotkina, E.A.; Romanenko, A.V.; Sazhin, A.F. Viruses and viral infection of heterotrophic prokaryotes in shelf waters of the western part of the East Siberian Sea. J. Mar. Syst. 2021, 218, 103544. [Google Scholar] [CrossRef]
  66. Hodges, L.R.; Bano, N.; Hollibaugh, J.T.; Yager, P. Illustrating the importance of particulate organic matter to pelagic microbial abundance and community structure—An Arctic case study. Aquat. Microb. Ecol. 2005, 40, 217–227. [Google Scholar] [CrossRef] [Green Version]
  67. Kirchman, D.L.; Elifantz, H.; Dittel, A.I.; Malmstrom, R.R.; Cottrell, M.T. Standing stocks and activity of Archaea and Bacteria in the western Arctic Ocean. Limnol. Oceanogr. 2007, 52, 495–507. [Google Scholar] [CrossRef]
  68. Turner, J.T. The importance of small planktonic copepods and their roles in pelagic marine food webs. Zool. Stud. 2004, 43, 255–266. [Google Scholar]
  69. Biswas, S.; Tiwari, P.K.; Bona, F.; Pal, S.; Venturino, E. Modeling the avoidance behavior of zooplankton on phytoplankton infected by free viruses. J. Biol. Phys. 2020, 46, 1–31. [Google Scholar] [CrossRef] [PubMed]
  70. Biswas, S.; Tiwari, P.K.; Kang, Y.; Pal, S. Effects of zooplankton selectivity on phytoplankton in an ecosystem affected by free-viruses and environmental toxins. Math. Biosci. Eng. 2020, 17, 1272–1317. [Google Scholar] [CrossRef]
  71. Rivkin, R.B.; Anderson, M.R.; Lajzerowicz, C. Microbial processes in cold oceans: Relationship between temperature and bacterial growth rate. Aquat. Microb. Ecol. 1996, 10, 243–254. [Google Scholar] [CrossRef] [Green Version]
  72. He, J.; Zhang, F.; Lin, L.; Ma, Y.; Chen, J. Bacterioplankton and picophytoplankton abundance, biomass, and distribution in the Western Canada Basin during summer 2008. Deep-Sea Res. II 2012, 81–84, 36–45. [Google Scholar] [CrossRef]
Figure 1. Location of sampling stations in the northeastern Barents Sea, in October 2020.
Figure 1. Location of sampling stations in the northeastern Barents Sea, in October 2020.
Biology 11 00845 g001
Figure 2. Vertical profiles of temperature, °C (a), and salinity (b) along the transects in the northeastern Barents Sea, in October 2020. Distance shows the intervals between sampling points.
Figure 2. Vertical profiles of temperature, °C (a), and salinity (b) along the transects in the northeastern Barents Sea, in October 2020. Distance shows the intervals between sampling points.
Biology 11 00845 g002
Figure 3. Zooplankton composition, abundance ((a)—individuals·m−3) and carbon biomass ((b)—mgC·m−3) in the 0–100 m layer in the northeastern Barents Sea, in October 2020. Water masses: BSW—Barents Sea Water, ArW—Arctic Water.
Figure 3. Zooplankton composition, abundance ((a)—individuals·m−3) and carbon biomass ((b)—mgC·m−3) in the 0–100 m layer in the northeastern Barents Sea, in October 2020. Water masses: BSW—Barents Sea Water, ArW—Arctic Water.
Biology 11 00845 g003
Figure 4. Vertical profiles of bacterial abundance ((a), 105 cells·mL−1), bacterial biomass ((b), mgC·m−3), and viral abundance ((c), 106 particles·mL−1), along the transects preformed in the northeastern Barents Sea, in October 2020.
Figure 4. Vertical profiles of bacterial abundance ((a), 105 cells·mL−1), bacterial biomass ((b), mgC·m−3), and viral abundance ((c), 106 particles·mL−1), along the transects preformed in the northeastern Barents Sea, in October 2020.
Biology 11 00845 g004
Figure 5. Abundance of bacteria ((a)—small, (b)—ultra small, (c)—large, (d)—chain cells) at each station in the northeastern Barents Sea, in October 2020.
Figure 5. Abundance of bacteria ((a)—small, (b)—ultra small, (c)—large, (d)—chain cells) at each station in the northeastern Barents Sea, in October 2020.
Biology 11 00845 g005
Figure 6. Cluster dendrogram (a) and NMDS ordination plot (b), showing relationships between stations, based on microbial parameters in the northeastern Barents Sea, in October 2020. Water masses: BSW—Barents Sea Water, ArW—Arctic Water.
Figure 6. Cluster dendrogram (a) and NMDS ordination plot (b), showing relationships between stations, based on microbial parameters in the northeastern Barents Sea, in October 2020. Water masses: BSW—Barents Sea Water, ArW—Arctic Water.
Biology 11 00845 g006
Figure 7. Redundancy analysis (RDA) on lg(x + 1)-transformed microbial abundance, indicating the ordination of microbial parameters and environmental variables in the northeastern Barents Sea, in October 2020. NB—bacterial abundance (105 cells·mL−1), BB—bacterial biomass (mgC·m−3), NV—viral abundance (106 particles·mL−1), VBR—the ratio of viral to bacterial abundance, ABV—average bacterial cell volume, R—contact rate between viruses and bacteria. T—temperature (°C), Sal—salinity, Chl a—surface chlorophyll a concentration (mg·m−3), PO4—phosphate concentration (µM), O2—dissolved oxygen concentration (mL·L−1), NO3—nitrate concentration (µM), ZA—zooplankton abundance in the upper 100 m layer (ind.·m−3), ZB—zooplankton biomass in the upper 100 m layer (mgC·m−3). Significant variables (forward selection procedure, Monte Carlo permutation test, p < 0.05) are marked with asterisks.
Figure 7. Redundancy analysis (RDA) on lg(x + 1)-transformed microbial abundance, indicating the ordination of microbial parameters and environmental variables in the northeastern Barents Sea, in October 2020. NB—bacterial abundance (105 cells·mL−1), BB—bacterial biomass (mgC·m−3), NV—viral abundance (106 particles·mL−1), VBR—the ratio of viral to bacterial abundance, ABV—average bacterial cell volume, R—contact rate between viruses and bacteria. T—temperature (°C), Sal—salinity, Chl a—surface chlorophyll a concentration (mg·m−3), PO4—phosphate concentration (µM), O2—dissolved oxygen concentration (mL·L−1), NO3—nitrate concentration (µM), ZA—zooplankton abundance in the upper 100 m layer (ind.·m−3), ZB—zooplankton biomass in the upper 100 m layer (mgC·m−3). Significant variables (forward selection procedure, Monte Carlo permutation test, p < 0.05) are marked with asterisks.
Biology 11 00845 g007
Table 1. Location and hydrology of sampling stations in the northeastern Barents Sea, in October 2020 [27].
Table 1. Location and hydrology of sampling stations in the northeastern Barents Sea, in October 2020 [27].
IDDateNEDepth, mTempSalWater Mass
708.10.202077°23′67°44′360−1.28–2.5133.95–34.89BSW
809.10.202077°50′66°46′330−0.98–1.5534.21–34.94BSW
909.10.202078°19′65°41′365−1.50–1.4334.08–34.88ArW
1009.10.202078°47′64°32′355−1.27–0.9833.92–34.84ArW
1109.10.202079°15′63°16′255−1.54–0.6133.49–34.83ArW
1210.10.202079°25′59°24′247−1.79–1.0733.72–34.83ArW
1310.10.202078°56′60°41′192−1.66–0.0933.57–34.81ArW
1410.10.202078°26′61°55′290−1.33–1.4133.88–34.84ArW
1510.10.202077°57′63°27′340−1.10–2.0534.02–34.85ArW
Note. Temp—temperature (°C), Sal—salinity. Water mass: BSW—Barents Sea Water, ArW—Arctic Water.
Table 2. Environmental parameters and microbial variables in different water masses and sampling layers in the northeastern Barents Sea, in autumn 2020.
Table 2. Environmental parameters and microbial variables in different water masses and sampling layers in the northeastern Barents Sea, in autumn 2020.
ParameterA
Surface (0–25 m)
B
Intermediate (50–100 m)
C
Bottom (192–365 m)
RangeMean ± SERangeMean ± SERangeMean ± SE
Arctic Water
NB, 105 cell·mL−13.33–9.255.67 ± 0.35 C2.96–10.686.12 ± 0.77 C2.27–4.093.20 ± 0.27 AB
BB, mgC·m−33.4–12.56.2 ± 0.6 C4.1–127.7 ± 0.8 C2.9–8.24.5 ± 0.7 AB
VBR1–2813 ± 27–1411 ± 14–149 ± 1
ABV, µm30.02–0.0740.037 ± 0.003 BC0.033–0.0850.049 ± 0.005 A0.031–0.1360.055 ± 0.014 A
R0.8–7.14.0 ± 0.53.0–9.65.7 ± 0.70.3–1.40.9 ± 0.2 AB
NV, 106 particles·mL−10.9–157.1 ± 1.0 C5.7–10.17.6 ± 0.5 C1.4–4.32.9 ± 0.5 AB
Temperature, °C−0.2–2.01.0 ± 0.2 BC−1.7–0.4−1.2 ± 0.2 A−1.3–−0.8−1.0 ± 0.1 A
Salinity33.49–34.0833.84 ± 0.06 BC34.06–34.7434.63 ± 0.0634.73–34.8434.81 ± 0.02
PO4, µM0.0–0.20.1 ± 0.0 BC0.0–0.70.4 ± 0.1 AC0.3–0.90.6 ± 0.1 AB
O2, mL·L−17.8–8.08.0 ± 0.0 C7.0–8.27.6 ± 0.1 C6.8–7.06.9 ± 0.0 AB
NO3, µM0.1–1.30.6 ± 0.1 BC0.4–16.18.4 ± 1.4 AC15–18.117.2 ± 0.4 AB
Chl a, mg·m−30.19–0.360.28 ± 0.02
ZA, ind.·m−3 (0–100 m)752–1278, 1067 ± 35
ZB, mgC·m−3 (0–100 m)0.8–11.8, 6.3 ± 0.6
Barents Sea Water
NB, 105 cell·mL−14.46–6.365.47 ± 0.27 C4.33–5.805.06 ± 0.73 C2.20–3.532.87 ± 0.67 AB
BB, mgC·m−34.6–9.66.1 ± 0.8 C5.6–6.15.8 ± 0.2 C2.7–5.54.1 ± 1.4 AB
VBR0–1511 ± 214–1515 ± 1 C10–1110 ± 0
ABV, µm30.026–0.0560.037 ± 0.005 C0.03–0.0510.041 ± 0.01 C0.042–0.0590.051 ± 0.008 AB
R2.2–6.14.3 ± 0.7 C2.8–4.73.7 ± 0.9 C0.5–1.40.9 ± 0.4 AB
NV, 106 particles·mL−14.4–9.17.0 ± 0.7 C6.0–8.97.5 ± 1.4 C2.4–3.63.0 ± 0.6 AB
Temperature, °C1.5–2.52.0 ± 0.2 BC−0.4–−0.2–0.3 ± 0.1 AC−1–−0.9–0.9 ± 0.0 AB
Salinity33.95–34.2734.14 ± 0.06 BC34.62–34.734.66 ± 0.04 A34.83–34.8334.83 ± 0.01 A
P-PO4, µM0.1–0.10.1 ± 0.0 BC0.3–0.50.4 ± 0.1 AC0.6–0.70.7 ± 0.0 AB
O2, mL·L−17.6–7.87.7 ± 0.0 C7.0–7.77.4 ± 0.3 C6.8–6.96.8 ± 0.0 AB
N-NO3, µM0.7–1.40.9 ± 0.1 BC2.5–3.93.2 ± 0.7 AC18.1–18.718.4 ± 0.3 BC
Chl a, mg·m−30.18–0.320.27 ± 0.03
ZA, ind.·m−3 (0–100 m)930–1008, 970 ± 13
ZB, mgC·m−3 (0–100 m)1.5–1.9, 1.7 ± 0.1
Note. NB—bacterial abundance, BB—bacterial biomass, VBR—the ratio of viral to bacterial abundance, ABV—average bacterial cell volume, R—contact rate between viruses and bacteria, NV—viral abundance. Chl a—surface chlorophyll a concentration (mg·m−3), PO4—phosphate concentration (µM), O2—dissolved oxygen concentration (mL·L−1), NO3—nitrate concentration (µM), ZA—zooplankton abundance in the upper 100 m layer (ind.·m−3), ZB—zooplankton biomass in the upper 100 m layer (mgC·m−3). The letters show significant differences between sampling layers (p < 0.05).
Table 3. Results of Monte Carlo permutation test in RDA: list of environmental variables affected microbial plankton in the northeastern Barents Sea, in autumn 2020.
Table 3. Results of Monte Carlo permutation test in RDA: list of environmental variables affected microbial plankton in the northeastern Barents Sea, in autumn 2020.
VariableLambdaAPFVariance Inflation Factor (VIF)Explained Variation, %
Nitrate0.270.00110.239.45527
Zooplankton biomass0.120.0055.444.95912
Salinity0.060.0453.015.5016
Temperature0.050.0692.775.3685
Dissolved oxygen0.050.0682.659.8455
Phosphate0.050.112.398.0115
Zooplankton abundance0.010.5140.72.7021
Chlorophyll a0.0010.9130.141.585<1
Table 4. Pearson’s correlation coefficients between environmental parameters and microbial variables in the northeastern Barents Sea, in autumn 2020.
Table 4. Pearson’s correlation coefficients between environmental parameters and microbial variables in the northeastern Barents Sea, in autumn 2020.
ParameterTemperatureSalinityPO4O2NO3Chl aZAZB
NB, 105 cell·mL−10.26−0.50−0.440.64−0.63−0.030.080.19
BB, mgC·m−3−0.02−0.28−0.230.34−0.340.150.250.28
NV, 106 particles·mL−10.40−0.17−0.290.46−0.52−0.26−0.21−0.37
ABV, µm30.290.13−0.040.10−0.18−0.26−0.29−0.54
VBR−0.350.270.24−0.350.380.310.320.15
R0.42−0.35−0.430.60−0.65−0.13−0.07−0.14
Note. NB—bacterial abundance, BB—bacterial biomass, VBR—the ratio of viral to bacterial abundance, ABV—average bacterial cell volume, R—contact rate between viruses and bacteria, NV—viral abundance. T—temperature (°C), Sal –salinity, Chl a—surface chlorophyll a concentration (mg·m−3), PO4—phosphate concentration (µM), O2—dissolved oxygen concentration (mL·L−1), NO3—nitrate concentration (µM), ZA—zooplankton abundance in the upper 100 m layer (ind.·m−3), ZB—zooplankton biomass in the upper 100 m layer (mgC·m−3). Bold font indicates significant p-values, confirmed by the Benjamini–Hochberg procedure.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Dvoretsky, V.G.; Venger, M.P.; Vashchenko, A.V.; Maksimovskaya, T.M.; Ishkulova, T.G.; Vodopianova, V.V. Pelagic Bacteria and Viruses in a High Arctic Region: Environmental Control in the Autumn Period. Biology 2022, 11, 845. https://doi.org/10.3390/biology11060845

AMA Style

Dvoretsky VG, Venger MP, Vashchenko AV, Maksimovskaya TM, Ishkulova TG, Vodopianova VV. Pelagic Bacteria and Viruses in a High Arctic Region: Environmental Control in the Autumn Period. Biology. 2022; 11(6):845. https://doi.org/10.3390/biology11060845

Chicago/Turabian Style

Dvoretsky, Vladimir G., Marina P. Venger, Anastasya V. Vashchenko, Tatyana M. Maksimovskaya, Tatyana G. Ishkulova, and Veronika V. Vodopianova. 2022. "Pelagic Bacteria and Viruses in a High Arctic Region: Environmental Control in the Autumn Period" Biology 11, no. 6: 845. https://doi.org/10.3390/biology11060845

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