Autumn to spring microbial community in the northern Baltic Sea: temporal variability in bacterial, viral and nanoflagellate abundance during the cold-water season

Marine microbial communities undergo drastic changes during the seasonal cycle in high latitude seas. Despite the dominance of microbial biomass in the oceans, comprehensive studies on the seasonal changes of microbial plankton during the complete winter period are lacking. To study the seasonal variation in abundance of the microbial community, water samples were collected weekly in the Northern Baltic Sea from October to May. During ice cover from mid-January to April, samples from the sea ice and the underlying water were taken in addition to the water column samples. Abundances of bacteria, virus-like particles, nanoflagellates, and chlorophyll a concentrations were measured from sea ice, under-ice water, and the water column, and examined in relation to environmental conditions. All studied organisms had clear seasonal changes in abundance, and the sea-ice microbial community had an independent wintertime development compared to the water column. Bacteria were observed to have a key role in the biotic interactions in both ice and the water column, and the dormant period during the cold-water months (October–May) was limited to before ice formation. Our results provide the first insights into the temporal dynamics of bacteria and viruses during the whole cold-water season (October–May) in coastal high latitude seas, and demonstrate that changes in the environmental conditions are likely to affect bacterial dynamics and have implications on trophic interactions.


Introduction
Microbial communities undergo changes in both abundance and community composition during a seasonal cycle in high latitude seas. During winter months when autotrophic production is limited, the importance of heterotrophic production by microbes is emphasized (Moreau et al. 2010;Enberg et al. 2018). The extreme seasonality in light availability in high latitudes limits primary production and results in only a short growth season that provides biomass for higher trophic levels for the entire year (Sakshaug 2004). Together with changes in the physicochemical properties of the water column, notably the annual ice cover, planktonic communities also exhibit pronounced seasonal succession patterns (Kuosa 1991;Pinhassi and Hagström 2000;Enberg et al. 2018). To understand how communities respond to the changing environmental conditions, it is crucial to examine the seasonal variation in the microbial community.
Given its subpolar location, seasonality is one of the key features of the Baltic Sea. A yearly ice cover extends on approximately 44% of the surface of the sea with a significant inter-annual variation in the extent and duration of the sea ice (Vihma and Haapala 2009). Baltic Sea ice is structurally similar to sea ice found in polar areas with interconnected brine pockets (Vihma and Haapala 2009) and hosts a diverse range of organisms living in the brine channels (Arrigo 2014;Thomas et al. 2017). Lower salinities in the Baltic Sea result in smaller brine volumes than in Polar sea ice, further limiting the available space (Weissenberger et al. 1992). Dissolved nutrients are concentrated within the brine pockets similarly to the salts present in the parent water, making the brine pockets an extremely nutrient concentrated habitat, with lower temperatures than in the underlying water column (Petrich and Eicken 2010). In addition to the brine channels, the under-ice water hosts a sub-ice habitat with organisms attached to or in close contact with the bottom ice (Boetius et al. 2015).
The initial microbial communities in sea ice are a result of the organisms incorporating in the brine pockets from the water column during ice formation (Gradinger and Ikävalko 1998;Rózańska et al. 2009;Eronen-Rasimus et al. 2015). Organisms may further be advected into the brine channels from the under-ice water through active brine movement (Hunt Jr et al. 2016). Given the limited space in the brine channels, Baltic Sea ice is dominated by smaller protists than Polar sea ice. The food web structure in ice thus differs from that of the water column, mainly as a result of absence of larger grazers (Krembs et al. 2000;Kaartokallio 2004;Caron et al. 2017). Heterotrophic bacteria are one of the most abundant organisms in the water column and sea ice both in terms of cell numbers and biomass (Lizotte 2003;Mock and Thomas 2005) and undergo seasonal changes in abundance and community composition (Mock et al. 1997;Pinhassi and Hagström 2000;Kaartokallio et al. 2008). Together with heterotrophic eukaryotes (Hassett et al. 2019), bacteria have a key role in nutrient recycling in ice, degrading particulate and dissolved organic matter and regulating nutrient availability (Zhou et al. 2014). Bacteria are grazed by a diverse group of nanoflagellates and other protists (Kuuppo-Leinikki 1990;Tophøj et al. 2018), which are an important link between bacteria and higher trophic levels (Azam et al. 1983).
Viral lysis is another important factor controlling bacterial abundance. High abundances of viruses have been observed in both in sea ice and the underlying water, and the number of viruses may be over tenfold to that of prokaryotic organisms (Whitman et al. 1998;Karner et al. 2001;Suttle 2005), with highest ratios in newly formed sea ice (Collins and Deming 2011). By lysing cells, viruses prevent the carbon flow to higher trophic levels, returning it to the microbial loop and controlling the nutrient and energy recycling in the water column (Wilhelm and Suttle 1999). While protistan grazing of bacteria transfers carbon to higher trophic levels, viral lysis releases carbon as dissolved organic matter, which can in some cases increase the activity of the remaining bacteria. Viruses are often very host-specific and can thus notably change the bacterial community composition. The relative importance of grazing and lysis in bacterial mortality further varies with depth and environmental conditions (Umani et al. 2010;Tsai et al. 2016). Given the key roles of microbial communities in biogeochemical cycling transferring energy and nutrients to higher trophic levels (Ejsmond et al. 2019), viral infections and shifts in community composition may have major ecological consequences (Proctor and Fuhrman 1990;Zhang et al. 2007).
Biogeochemical processes in the sea ice are important in mediating ecosystem processes also in the water column (Yakubov et al. 2019). Timing of sea-ice formation and its duration have been shown to be important for primary producers, and the vernal ice melt has an essential role for the benthic communities as organic matter sinks to the sea floor (Leu et al. 2015). Warming temperatures may have complex impacts on the phenologies of primary producers (Tedesco et al. 2019) and impact notably the size distribution of plankton communities (Daufresne et al. 2009). As increasing temperatures are expected to reduce seasonal ice cover in the Baltic Sea (Jylhä et al. 2008), the changing seasonality (Kahru et al. 2016) may therefore also affect the microbial community phenologies (Sydeman and Bograd 2009).
Given the dominance of microbial biomass in the oceans, predicting future changes in high latitude marine ecosystems requires emphasis on microbes and their temporal dynamics. Yet, proportionally few studies have addressed the succession patterns of microbial marine communities during a whole seasonal cycle, often focusing on specific taxonomic groups (Andersson et al. 1996;Pinhassi and Hagström 2000;Schauer et al. 2003;Mary et al. 2006;Iversen and Seuthe 2011;Hassett and Gradinger 2016). Comprehensive time series data on several important microbial variables have only been collected recently (Iversen and Seuthe 2011;Marquardt et al. 2016;Bunse et al. 2018), and especially viruses have been largely neglected (Sandaa et al. 2018). There is thus a pertinent need for more time series approaches to establish a baseline to observe potential shifts in marine microbial communities. Here, we study the seasonal changes in abundance of major groups of microorganisms in the water column, sea ice, and under-ice water. The focus of our investigation is on the abundance of heterotrophic bacteria, virus-like particles, and nanoflagellates from autumn to early summer. The term "microbial community" is used in this paper to describe the studied organism groups. We use the term bacteria to describe prokaryotic cells, although the observations may also include archaea.

Study site and sampling
Sampling was performed weekly in the northwest coast of Gulf of Finland (Fig. 1). The study site Storfjärden (59° 51.250′ N, 23° 15.815′ E) is approximately 30 m deep, characterized by salinity of 4-6, relative exposure to winds and thus prone to sudden sea-ice break-ups. Samples were collected from a boat or on the ice from October to May, with the exceptions of 24 and 31 Dec 2012, 7 Jan 2013, 11 to 18 Feb 2013 and 22 Apr 2013 when samples could not be collected due to poor sea-ice conditions. Water samples were collected using a hose sampler with an internal diameter of 6 cm from 0 to 15 m depth (Majaneva et al. 2009). Water temperature and salinity were measured with a Falmouth Scientific NXIC CTD.
During ice cover from January to April, samples from the sea ice and the underlying water were collected in addition to the water column samples. Ice sampling was performed using a motorized CRREL-type ice-coring auger (9 cm internal diameter; Kovacs Enterprises, Roseburg, OR, USA). Ice samples were rapidly melted by crushing without filtered sea water additions as described in Rintala et al. (2014). Prior to ice sampling, snow cover on ice was measured from three random replicate samples with 1 cm precision using a motorized CRREL-type ice-coring auger (9 cm internal diameter; Kovacs Enterprises). Sea-ice temperatures were measured with a Testo 110 thermometer, and the bulk salinities of the melted ice samples were measured with an YSI 63 m (Yellow Springs Instruments). Under-ice water samples were obtained by submerging a bottle in a hole drilled in ice. After ice melted, we continued to sample the surface water until the end of the study period from a similar depth than the under-ice water samples. Concentrations of inorganic nutrients (PO 4 , NH 4 , NO 2 + NO 3 , SiO 4 ) were determined with an interval of 2 to 3 weeks from all sampled habitats using a Hitachi U-110 Spectrophotometer (Hitachi High-Technologies) according to standard seawater protocols described by Hansen and Koroleff (1999). Nutrient concentrations in ice were normalized to the average bulk salinities of melted sea ice to adjust for salinity-related variations in the nutrient concentrations.
Chlorophyll a concentration was determined from two 100 ml aliquots of sample water (ice and water) filtered on 25 mm Whatman GF/F filters and extracted in 94% ethanol. After being incubated in dark in room temperature for a minimum of 24 h, the extract was filtered again through a Whatman GF/F filter and measured with a Cary Varian Eclipse spectrofluorometer. Chlorophyll a concentrations were calculated according to HELCOM (1988).
Abundances of virus-like particles (VLP) and bacterial cells were determined using a CyFlow© Cube 8 flow cytometer as per instructions of Manual of Viral Ecology  (Brussaard et al. 2010). Subsamples of 1.5 ml were obtained from all of the water and sea-ice samples and placed into plastic cryovials after which they were preserved with paraformaldehyde (final concentration 1%). The samples were left to settle for 30 min before flashfreezing with liquid nitrogen and were stored at − 80 °C (Marie et al. 1999). Prior to enumeration, samples were thawed in room temperature, diluted in 1:10 in sterile TEbuffer (AppliChem GmbH, Darmstadt, Germany, pH 8) and stained for 15 min in the dark with SYBR Green I nucleic acid stain (final concentration 0.0036%, Molecular Probes). After staining, 0.5 µm fluorescent microspheres (Molecular Probes) were added into the solution for size reference. Samples were analyzed using a flow rate of 12 μl min −1 . The data were acquired on a dot plot displaying green fluorescence (488 nm) versus side scatter signal, both on a logarithmic scale. The detection trigger was depicted as a green fluorescence. The obtained data were analyzed with a FCS Express 4-software (De Novo Software).
Bacterial production was measured using the 3H-thymidine incorporation method Azam 1980, 1982) to examine bacterial net biomass production based on DNA synthesis. Analyses were carried out at 1 to 3 week intervals on both ice and water samples. The ice samples were prepared as described in Kaartokallio (2004). Duplicate 10 ml subsamples and a formalin-killed blank were incubated with [methyl-3H]-thymidine at a saturation level (25 nM, specific activity 20 Ci mmol −1 ) in dark at − 0.2° for 2-19 h according to the predicted level of activity. Incubation was stopped by adding formaldehyde and samples were processed with the standard cold-TCA extraction method using 0.2 μm mixed cellulose ester filters (Advantec MFS). The filters were dissolved in InstaGel scintillation cocktail (PerkinElmer) and the incorporated radioactive thymidine was assayed with a Wallac WinSpectral 1414 liquid scintillation counter (PerkinElmer).
For nanoflagellate enumeration, a subsample of 20 ml was taken from each water and ice sample and preserved with 1 ml particle free 25% glutaraldehyde. Samples were stored at + 4 °C. Nanoflagellates (< 10 µm) were enumerated with a Leitz Aristoplan epifluorescence microscope under blue light using a 100× objective magnification together with 12.5 × ocular magnifications. Five to ten ml of sample was filtered onto 0.2 µm black track etched membrane filters (Whatman) and stained with 50 to 100 µl 0.033% proflavine solution (Haas 1982). A minimum of 30 cells or 50 fields of vision were counted from each sample with the accuracy of ± 28% at a 95% confidence interval (Lund et al. 1958). Separation of autotrophic and heterotrophic nanoflagellates was not performed systematically.

Statistical analysis
A repeated measures analysis of variance (ANOVA) was applied to test for differences between organism abundances between sampled habitats at specific sampling days. Normality checks and Levene's test were carried out and the assumptions were met. Bivariate correlations were used to examine relations between the studied organisms and environmental parameters. To account for the temporal autocorrelation in the data and to specifically remove seasonal variation in the time series, residuals were first calculated by subtracting monthly mean values from the corresponding cell abundance values and environmental data values (Chatfield 2016). Pearson's product-moment correlation analysis was then applied on the residuals to examine linear relations between the seasonally adjusted organism abundances and environmental parameters (Tsai et al. 2013). Stationarity of the residuals and absence of serial correlation were verified with autocorrelation and partial autocorrelation function plots provided by R 3.6.1 (R Core Team 2019) base package prior to correlation tests.
Due to the potentially changing direction of interactions between organisms during the seasonal cycle, we also examine the seasonal interactions between organisms in the water column while maintaining the time series structure of the data. In a time series, the value of an observation at a given point is not truly random or independent, but depends on the values of preceding time points (Borcard et al. 2018). Due to this temporal dependency and strong autocorrelation of the observations, traditional multivariate analysis techniques are not statistically adequate for studying these relations. To account for these issues, we applied a first order vector autoregressive model (VAR(1)) which allows the variables to be modeled jointly over present and past time periods. An autoregressive model using the vars package in R (Pfaff 2008a, b) was applied only on water sample data given the relatively small number of samples in ice, using a lag of 1 week. To examine the interactions proposed by the model, we tested how well values of one time series can be used to predict the values of another (Faust et al. 2015) by performing Granger causality tests between groups for which the model suggested significant interactions (Granger 1969). All analyses were performed with R 3.6.1 (R Core Team 2019).

Environmental conditions
Winter 2012-2013 was characterized by long persisting sub-zero temperatures and ice cover from late December to mid-April. Ice formed at the sampling station between January 11th and January 14th, reaching a maximum depth of 52 cm on March 25th (Fig. 2). The ice cover melted partially between February 11th and February 18th, after which the ice was reformed and increased in depth until ice melt in April. Snow cover on the ice was relatively low throughout the study period, and varied between 0 and 3 cm (Fig. 2). Water temperatures decreased rapidly from October to December and remained between 0 and 1 °C from January to April (Fig. 3). A decline in salinity was observed in the under-ice water from January to April, with a minimum salinity of 2.4 pointing to a freshwater lens under the ice cover (see supplementary material in Enberg et al. 2018). The minimum salinity in the water column was 4.7 (Fig. 3). Concentrations of phosphate, silicate and mono-nitrogen oxides in the water column were highest from mid-December through March (Online Resource 1), with maximum concentrations of 0.9, (Phosphate), 26.9 (Si), 10.2 (nitrate) mmol m −3 . In ice, phosphate concentration decreased throughout the entire ice-covered period, with a minimum concentration 0.1 mmol m −3 (Online Resource 1).

Chlorophyll a
Chlorophyll a (chl a) concentration in the water column declined slightly in October and stayed below 3 µg l −1 until April, with the minimum of 0.312 µg l −1 and a maximum concentration of 10.87 µg chl a l −1 (Fig. 4). No clear spring bloom was observed during the study period. In sea ice, bulk chlorophyll a concentrations were higher than in the water column. The chl a in sea ice remained low until a sudden increase in February, with chl a concentration peaking at 21 µg chl a l −1 . This coincided with the melting and reformation of the ice field and corresponds to the maximum chlorophyll a concentration observed during the study period.

Seasonal changes in bacterial abundance and production
Bacterial abundance in the water column decreased from 5 × 10 6 cells ml −1 in October to 2.2 × 10 6 cells ml −1 in February (Fig. 5, panel a). The minimum abundance was 2-2.5 × 10 6 cells ml −1 in March-April. After the ice melt, bacterial abundances increased, and the maximum was observed during the last sampling in late May with bacterial numbers reaching 7 × 10 6 cells ml −1 . During winter, the overall bacterial abundance was 2-2.5 million cells both in the under-ice water and the water column, but the abundance at specific sampling dates during the entire sampling period differed significantly between the sampled habitats (repeated measures ANOVA, F 2 = 89.7, p < 0.0001). Bulk sea-ice bacteria increased in numbers slightly a week after ice formation, then decreasing until the last week of ice winter. During winter, abundances in under-ice water were similar to those observed in the water column. After ice melt, bacterial abundance in the water column increased rapidly until the end of sampling period.
Total thymidine incorporation (TTI) rates, a proxy for bacterial cell production, declined steeply compared to bacterial abundance during autumn (Fig. 5, panel b). The minimum productive period in the water column was observed from mid-November to mid-January, with the total thymidine incorporation on average 0.5 pmol l −1 h −1 . The highest bacterial production rates were observed in the under-ice water, with the production peak of 6.42 TTI pmol l −1 h −1 on January 21. Thymidine incorporation rates in ice increased from the initial value of 0.51 after ice formation to 0.68 pmol l −1 h −1 end of February, and later peaking at 1.94 pmol l −1 h −1 in late March. Over the whole study period, cell-specific production rates were the highest in under-ice water (Fig. 5, panel c), yet the maximum specific productivity was measured from ice 2 weeks before ice melt. While bacterial abundance declined through the study period, cell-specific production had an opposite evolution, with increasing bacterial activity throughout the season.

Seasonal changes in VLP abundance
Abundance of virus-like particles in the water column increased from 7 × 10 6 ml −1 in October to 1.4 × 10 7 ml −1 in December, rapidly decreasing to 5 × 10 6 ml −1 prior to ice formation in January (Fig. 6). Most notable temporal variation was observed in the under-ice water where the VLP abundances were 2-6 times higher than in the water column

Seasonal changes in nanoflagellate abundance
Nanoflagellate abundance was the highest in October and decreased abruptly in November (Fig. 7). Maximum abundance during the study period was observed in October with nanoflagellate abundance reaching 12,500 cells ml −1 . During winter, cell abundances were 0-2500 cells ml −1 . Abundances in sea ice increased throughout the ice-covered period, and the maximum abundance of 7500 nanoflagellate cells ml −1 was observed in April before melting of the sea ice. The highest vernal cell counts were found in sea ice and the under-ice water only days before ice melted. A more detailed description of the nanoflagellate succession is given in Enberg et al. (2018) (site B), where the algal biomass during the ice-covered period was dominated by small flagellates.

Relation of temporal abundance patterns to environmental conditions
Water temperature correlated positively with bacteria and VLP abundance, (Table 1) whereas in ice no significant correlation was observed between organisms and the ice temperature (Table 2). Similarly, the VAR(1) model results (Table 3) and Granger causality test show that water temperature was found to Granger-cause bacterial abundance (F 23,−1 = 5.9253, p = 0.0235). With the exception of VLP abundance correlated negatively with water column salinity (r = 0.458, n = 29, p = 0.0243, Table 1), no other significant correlations were observed between water or ice salinity and abundance of the studied organisms. Bacterial and VLP abundance correlated negatively with ice depth both in the water column and within the ice (Table 2). Bacteria were found negatively correlated with inorganic nutrient concentrations both in the water column and in ice (Tables 1, 2).

Microbial community interactions
Pearson's product-moment correlation suggests strong positive correlation between the seasonally adjusted bacterial abundance and the number of VLPs in the water column Fig. 3 Water column, under-ice water and bulk ice temperature and salinity in Tvärminne, Gulf of Finland (59° 51.250′ N, 23° 15.815′ E). The timeline from the beginning of January to mid-April represents the icecovered season (r = 0.733, n = 29, p < 0.0001) ( Table 1). A similar positive correlation (r = 0.756, n = 12, p = 0.004) was observed in ice (Table 2). In ice a significant, albeit relatively low, negative correlation (r = − 0.367, n = 12, p = 0.007) was detected between bacteria and nanoflagellates.
VAR(1) model diagnostics indicate that bacteria Granger-cause VLP, nanoflagellates, and chlorophyll a concentration (F 3,76 = 4.3609, p = 0.0069). This result can be interpreted so that changes in bacterial abundance can be used to predict variation in the abundance of the other organisms. Further, changes in chl a concentration were found to precede changes in nanoflagellate abundance concentration (F 23,-1 = 4.2781, p = 0.0506).
In terms of biological interactions, the statistically significant coefficients of the VAR(1) model and the corresponding tests of Granger causality suggest that bacterial abundance and chlorophyll a concentration at a given time will only be affected by the abundance of the previous sampling of the variable itself, thus the week before in this case (Table 3). While VLP lagged abundance was a significant predictor for bacterial abundance, the interaction did not pass the Granger causality test (F 23,−1 = 0.5072, p = 0.4838). The coupled results from the correlation tests and VAR model show that the abundances of organisms are affected by both the short-term interactions, and lagged values from the precedent week.

Discussion
In this study, we identified distinctive changes in bacterial, VLP, and nanoflagellate abundance in sea ice, the underice water, and the water column from October to May in the northern Baltic Sea. The most prominent differences between the sampled habitats were observed between ice and water column samples, and the sea-ice microbial community had an independent wintertime development from that of the water column. Bacteria, nanoflagellates, and chlorophyll a concentrations in ice decreased or increased more evidently than in the water during the ice-covered period.
During the same seasonal cycle, Enberg et al. (2018) described the succession of the eukaryotic communities Fig. 4 Chlorophyll a concentration in the water column, bulk ice and under-ice water in Tvärminne, Gulf of Finland (59° 51.250′ N, 23° 15.815′ E). Error bars represent standard deviations from replicate samples. The timeline from the beginning of January to mid-April represents the ice-covered season in ice and water column using microscopy and the 18S rRNA gene sequence analysis. Together with Majaneva et al. (2019), they show that different species of flagellates characterize sea ice, under-ice water, and the water column during the sea-ice-covered season. Coupled with results of, not to this study, the distinct communities in the different sea-ice habitats also seem to extend to the abundances of the heterotrophic microbial community.
Abundance patterns of organisms are a result of both environmental factors and biological interactions. Previous time series studies have found bacterial abundance to be connected to both abiotic and biotic factors (Bunse et al. 2018). Here, bacterial increase in the spring coincided with rapidly increasing water temperatures. Our results thus suggest that environmental factors and the history of bacterial abundance affect bacterial abundance more than biological interactions. By examining both the lagged and instantaneous response of organisms, our results show a strong connection between bacterial abundance and temperature. Bacterial growth is also controlled by the production of labile dissolved carbon extracted by algae (Kuparinen et al. 2007), and in our study, the sudden increase in bacterial abundance in the spring seemed to be more connected to the melting of ice than phytoplankton growth. Bacterial community development was thus likely associated with the availability of the algalderived substrate (Eronen-Rasimus et al. 2017). Bacterial . Error bars represent standard deviations from replicate samples. The timeline from the beginning of January to mid-April represents the ice-covered season abundance was not directly correlated to chl a concentration, and the carbon incorporated into the algal and nanoflagellate biomass in ice was therefore likely to promote the growth of bacteria through dissolved organic matter transport after ice melted (Kuparinen 1988;Sime-Ngando et al. 1999;Sala et al. 2010). A similar lack of correlation between bacteria and chl a has been shown in other mid-winter studies (Stewart and Fritsen 2004), likely resulting from the relatively low chl a concentrations in ice. This supports previous findings that combined effect of temperature and substrate availability would be most efficient in increasing bacterial abundance (Autio 1998;Kirchman et al. 2005), highlighting the importance of the timing of ice melt for microorganisms.
The bacterial production results of this study show that the dormant period during the winter months is limited to before ice formation, and the heterotrophic production then increases throughout the winter and spring months. Bacterial activity in newly formed sea ice is often suppressed before consolidation of sea ice (Grossmann and Gleitz 1993;Kaartokallio et al. 2008), yet we found that the cell-specific production rates in ice were restored quickly after ice formation. While the low snow cover in our study may have exposed the sea-ice communities to higher solar radiation in the spring, the effects of UVA on bacterial production are limited to surface layers of ice (Piiparinen and Kuosa 2011), and the cell-specific production rates increased throughout the study period. Bacterial production was highest in the under-ice water and within the same ranges in ice than measured by Kaartokallio et al. (2008). While our results show high bacterial abundance during the winter minima, they are in line with earlier coastal studies where bacterial abundances were found to remain relatively low until the melting of the ice (Kuuppo-Leinikki 1990;Heinänen and Kuparinen 1991;Tuomi 1997;Kaartokallio 2004).
The highest nanoflagellate abundances were observed at the end of the growth season when also bacterial abundances Fig. 6 VLP abundance in the water column, bulk ice and under-ice water in Tvärminne, Gulf of Finland (59° 51.250′ N,23° 15.815′ E). Error bars represent standard deviations from replicate samples. The timeline from the beginning of January to mid-April represents the ice-covered season were the highest, potentially providing sufficient resources for nanoflagellate growth (Kuosa and Kivi 1989;Kivi et al. 1993). Although traditionally autotrophic and heterotrophic flagellates are treated separately, even chlorophyll containing nanoflagellates have been observed to consume bacteria (Estep et al. 1986). Mixotrophic behaviour, which combines both autotrophic and heterotrophic modes of energy acquisition, is widely spread among flagellates (Sanders and Gast 2012). Following the eukaryote succession during the same seasonal cycle, Enberg et al. (2018) showed that cryptophytes and haptophytes were some of the most abundant flagellate taxa at the end of the ice-covered period, some of which are known capable of mixotrophy. The autoregressive model supports this finding by indicating that changes in nanoflagellate abundance can be explained by both bacterial abundance and chl a concentration. Both the increasing light intensities and higher bacterial biomass thus likely contributed to the development of the nanoflagellate biomass.
As bacteria are the most abundant hosts of viruses, viral production is often linked to bacterial production and abundance (Steward et al. 1996;Sandaa et al. 2018). In this study, VLP appeared in large quantities in the autumn at the beginning of the cold-water period. This phenomenon may be a result of induction of lysogenic viruses during the temperature decrease or viral control of the bacterial population during this period (Collins and Deming 2011). High virus abundances together with low bacterial abundance and activity have also been detected in newly formed sea ice (Collins and Deming 2011;Luhtanen et al. 2018). VLP abundances showed both great irregular temporal variation and distinctive abundance patterns in the sampled habitats, especially the under-ice water. Weekly averages of total virioplankton can have an annual variation up to one order of magnitude (Li and Dickie 2001), and the virus-to-host cell ratio decreases with microbial cell density (Wigington et al. 2016). The results from the autoregressive model support Fig. 7 Nanoflagellate abundance in the water column, bulk ice and under-ice water in Tvärminne, Gulf of Finland (59° 51.250′ N, 23° 15.815′ E). The timeline from the beginning of January to mid-April represents the ice-covered season this by showing that VLP and bacterial abundance are interconnected, and the changes in bacterial abundance (either increasing or decreasing) mediate VLP abundance.
While there is considerable natural variation in viral abundance, the observed temporal changes in VLP abundances in this study may in part be due to methodological limitations. We hypothesize that the high VLP abundances in the under-ice water during the winter-spring transition are a result of induction of lysogenic viruses in warming ice, and particle flushing from the enlarging brine channels. Although virioplankton abundances are typically higher in the surface layers of the ocean (Liang et al. 2014), we cannot exclude that the observed particles may be other organic particles with decayed nucleic acid fragments. However, compared to the water column and the brine channels in sea ice, the under-ice water is a more dynamic environment as it is exposed to convection at the ice water interface (Hudier et al. 1995). A shorter sampling interval would therefore be adequate to characterize the temporal viral dynamics (Sandaa et al. 2018) and interactions between other organisms in the marine environment. Furthermore, unraveling the complex interactions between viruses, their hosts, and environmental factors would require significantly longer time series to increase the number of data points.
Bacteria have a key role in the microbial food web, and we found fluctuations in bacterial abundance to precede changes in the abundances of VLP and nanoflagellates. As interactions between organisms may change within a seasonal cycle, assigning only unidirectional interactions and assessing their significance bears limited importance to understanding community dynamics within seasonal changes. Here, the use of an autoregressive model in addition to correlation analysis provided a more complete overview on the interactions between organisms. However, changes of biota in the ice and in the water column coupled with drastic variations in temperature and biotic interactions make the system quite complicated, and the independent effects of biological interaction cannot be fully depicted by the results of this study. While our study system is notably lower salinity than many other sub-Arctic marine areas, the changes observed here throughout the low-light season may reflect future conditions in other coastal high latitude systems (Carmack et al. 2016). Because of the importance of ice cover for timing of biomass peaks in the microbial community, changing ice conditions are likely to affect microbial community phenologies and interactions. Our results provide first insights into the temporal dynamics of bacteria and viruses during the whole cold-water season, and demonstrate that changes in the environmental conditions are likely to affect bacterial dynamics and alter trophic interactions.