Multiyear analysis uncovers coordinated seasonality in stocks and composition of the planktonic food web in the Baltic Sea proper
Scientific Reports volume 13, Article number: 11865 (2023) Cite this article
499 Accesses
5 Altmetric
Metrics details
The planktonic realm from bacteria to zooplankton provides the baseline for pelagic aquatic food webs. However, multiple trophic levels are seldomly included in time series studies, hampering a holistic understanding of the influence of seasonal dynamics and species interactions on food web structure and biogeochemical cycles. Here, we investigated plankton community composition, focusing on bacterio-, phyto- and large mesozooplankton, and how biotic and abiotic factors correlate at the Linnaeus Microbial Observatory (LMO) station in the Baltic Sea from 2011 to 2018. Plankton communities structures showed pronounced dynamic shifts with recurring patterns. Summarizing the parts of the planktonic microbial food web studied here to total carbon, a picture emerges with phytoplankton consistently contributing > 39% while bacterio- and large mesozooplankton contributed ~ 30% and ~ 7%, respectively, during summer. Cyanophyceae, Actinobacteria, Bacteroidetes, and Proteobacteria were important groups among the prokaryotes. Importantly, Dinophyceae, and not Bacillariophyceae, dominated the autotrophic spring bloom whereas Litostomatea (ciliates) and Appendicularia contributed significantly to the consumer entities together with the more traditionally observed mesozooplankton, Copepoda and Cladocera. Our findings of seasonality in both plankton composition and carbon stocks emphasize the importance of time series analyses of food web structure for characterizing the regulation of biogeochemical cycles and appropriately constraining ecosystem models.
The Baltic Sea is one of the world’s largest brackish systems1. It offers multiple socio-economic services, like fisheries, recreational boating, swimming and overall enjoyment of the seascape2. The Baltic Sea has been, and still is, heavily influenced by multiple anthropogenic stressors and over exploitation of ecosystem services with severe effects on the ecosystem, e.g., eutrophication and pollution3,4. The production of plankton communities is transferred to higher trophic levels. However, the mechanisms regulating the transfer efficiency are still not well defined and additionally, may be affected by anthropogenic stressors5. As such, the Baltic Sea can shed light on how ecosystem services, e.g., wild/farmed seafood, removing nutrients and pollutants and tourism, are affected by anthropogenic stressors and thereby offer an expanded view on ecosystem-based management to ensure a sustainable utilization of the Baltic Sea.
Primary producers and consumers like bacterio-, phyto- and zooplankton communities constitute the base of pelagic food webs, together with fungi and archaea, and their metabolic activity has major impacts on global elemental cycling6,7,8,9,10,11,12,13. The total flow of energy and matter to higher trophic levels in the food web is highly dependent on microbial community structure and functioning since heterotrophic microorganisms recycle nutrients and dissipate energy produced by primary producers14,15,16. Organisms at the base of the aquatic food web produce biomass and compounds which are important for higher trophic levels throughout the food web. For instance, both hetero- and autotrophic plankton play a key role in the synthesis of e.g. vitamins that are universally needed by higher organisms17. Moreover, fish affect the plankton communities, either directly by predation or indirectly by predator-release18,19. Therefore, regional long-term investigations of the seasonal dynamics of environmental parameters and bacterio-, phyto- and zooplankton have the potential to provide novel knowledge on the structure and functioning of the food web.
Bacterioplankton play a pivotal role in the aquatic food web, attributable to high abundances, number of species, functional properties, and numerous interactions with other microorganisms as well as higher trophic levels. In aquatic ecosystems, the major primary producers are phytoplankton20,21, which only make up 0.2% of the global primary producer biomass, but account for 50% of the net production, e.g.21,22. Phytoplankton, in turn are grazed by zooplankton and it is estimated that zooplankton can consume 10–40% of the primary production through grazing daily, depending on the productivity of the ecosystem23. Importantly, even small changes in community- or size-structures, and production can have complex and large-scale effects on the entire ecosystem including fish production4,24,25. Nevertheless, relatively few studies include multiple trophic levels within the plankton realm, and seldomly cover total biomass estimates along with measures of community composition and abiotic parameters8,26,27,28,29,30,31,32.
Temporal similarity in community composition can provide insights on the community development, or stability, as shown by Yeh and Fuhrman32 in the San Pedro Time Series. Here, the prokaryotic community was found to be quite constant over time whilst the protist community was substantially altered over time32. Microbial community composition is affected by interactions with other compartments of the food web (e.g., grazing, viral lysis, allelopathy, symbiosis, competition, and parasitism) as well as abiotic factors (e.g., temperature, water mixing, and nutrient availability). Studies which combine field observations of community composition and biotic and abiotic factors with a correlative approach can serve as a steppingstone in designing mechanistic experimental studies33. Here we present long-term and high-frequency data of biotic and abiotic parameters retrieved from the Western Baltic Sea Proper—the Linnaeus Microbial Observatory (LMO)—covering the period 2011–2018. Short-term subsets of some of these data have previously been used to analyze microbial dynamics and seasonal variability as well as micronutrient dynamics, e.g.,34,35,36,37,38,39,40,41. In this study, we substantially extend the analysis to bacterio-, phyto- and large mesozooplankton in the planktonic food web over a period of eight years including community dynamics at several trophic levels. The aim of the study is to provide one of the first planktonic food web descriptions of this part of the Baltic Sea and to determine if the seasonal dynamics of the different trophic levels are correlated and to what extent abiotic and/or biotic factors affect the plankton food web.
Samples were collected at the Linnaeus Microbial Observatory (LMO; N56°55.8540´, E17°3.6420), located approximately 11 km (6 nautical miles) off the NE coast of the island Öland with a depth of 40 m34,35,40. In short, sampling started in March 2011 and was performed with various frequency over the years, from twice weekly to monthly samplings. In the present study, data from samplings performed from 2011-03-25 to 2018-11-27 are included, covering 270 sampling cruises. Water was collected from 2 m depth using a 3 or 5 L Ruttner water sampler at 08.00–10.00 am local time. Temperature and salinity were measured on-site. Water was collected in 10 L, acid-washed, polycarbonate bottles and transported to the laboratory within 1 h. Large mesozooplankton were collected by oblique hauls from the top 30 m using a plankton net (Ø50 cm, 200-μm mesh size) with a fitted flowmeter and stored in a cooling box. Mesozooplankton sampling was included in the sampling regime in November 2013 and sampled monthly. This sampling procedure mainly catch larger sized mesozooplankton and has lower retention of smaller zooplankton such as rotifers, small cladocera, nauplii and other microzooplankton. In addition, starting in 2013, the water column was profiled using a CTD probe (AAQ 1186-H, Alec Electronics, Japan) for temperature, salinity, and light intensity. In the laboratory, samples for nutrients, chlorophyll a (Chl a), dissolved organic carbon (DOC), colored dissolved organic matter (cDOM), as well as abundance, biomass, and community composition for bacterio-, phyto- and large mesozooplankton were collected and analyzed, see below.
Inorganic nutrients; nitrate and nitrite (NO3− and NO2−, together presented as nitrate), ammonium (NH4+), phosphate (PO43−), and silicate (SiO44−) were collected and frozen until analysis using colorimetric methods (UV-1600 Spectrometer, VWR and DR 5000, Hach Lange) according to Valderama42. Detection limits for inorganic nutrients were 0.02 µM for NO3− and PO43− and 0.1 µM for NH4+ and SiO44−. Chl a was sampled by filtration of ~ 500 mL water in duplicates on A/E filters (Pall Laboratory). Filters were extracted overnight in 96% ethanol in the dark and later analyzed according to Jespersen and Christoffersen43, using a Trilogy flourometer (Turner Designs, USA). DOC was sampled in duplicates; 20 mL of water was filtered through precombusted (475 °C, 2 h) GF/C glass fiber filters (nominal pore size ~ 1.2 µm) via gravity, acidified (200 µL 2 M HCl) and stored in precombusted glass vials (475 °C, 3 h) with acid washed lids at 6-8 °C until analysis. DOC content was measured via high temperature catalytic oxidation (HTCO) using a Shimadzu TOC-V analyzer coupled to a TNM-1 unit as well as determined with a Shimadzu TOC 5000 analyzer and a Shimadzu TOC-L Total Organic Carbon Analyzer using an acetanilide dilution series as a standard44,45. Detection limit for DOC measurement was ~ 10 µM. cDOM was analyzed by filtering (0.2 µm, Sterivex cartridge filters) a sample into an amber glass bottle, and the absorption coefficients were measured in a 10 cm quartz cuvette over the 420–780 nm range with 2 nm increments in a UV-1600 Spectrometer, VWR46. cDOM corresponds to the absorption coefficient at 440 nm.
Samples for bacterioplankton abundance were preserved with formaldehyde (3.7% final concentration) and kept at − 80 °C until analysis. Samples for total bacterioplankton abundance were enumerated with a flow cytometer (BD FACs Calibur (2011–2012) and Partec Cube8 (2013–2018)) with protocols adapted from Gasol and Morán47. Bacterioplankton abundance was converted to carbon biomass using a factor of 20 fg C cell-148. Nucleic acids from bacterioplankton for amplicon analysis were collected without prefiltration by filtering up to 10 L seawater on 0.2 µm Sterivex cartridge filters (Millipore). Filters were stored frozen at − 80 °C in TE buffer until DNA was extracted using a phenol–chloroform extraction as described by Boström, et al.49 and modified after Bunse, et al.50. The V3V4 region of the 16S rRNA gene was amplified using PCRs with the primer pair 341F-805R51 as described and validated in Hugerth, et al.52. DNA concentrations were analyzed using Qubit 2.0 Fluorometer (Life Technologies) and subsequent gel electrophoresis confirmed amplicon specificity. Sequencing was carried out at the Science for Life Laboratory, Sweden on a MiSeq platform (Illumina, Inc.), producing 2 × 300 bp paired-end reads. For bioinformatical processing, we used the Ampliseq pipeline (https://github.com/nf-core/ampliseq)53 using DADA2 to infer amplicon sequence variants (ASVs)54. The versions of the software used to produce these results were: nf-core/ampliseq (v1.2.0dev); Nextflow (v20.10.0); FastQC (v0.11.8); MultiQC (v1.9); Cutadapt (v2.8); QIIME2 (v2019.10.0). ASVs were taxonomically annotated against the SILVA database (v132). The resulting amplicon sequence variant (ASV) table was subsequently analyzed in R (see below). Relative abundances were normalized per genome size, using genome “16S rRNA/genome” published by Vetrovsky and Baldrian55 and relative abundances represent percentages per taxon. Note that cyanobacteria were deleted from the ASV-generated taxonomy and is only presented as cell counts within the phytoplankton fraction. This is in order to not have this group in both the bacteria and phytoplankton category when analyzing community composition.
Separate samples for phytoplankton and large mesozooplankton abundances were preserved with acidic Lugol’s solution (2% final concentration) and stored in the dark until analyses. Phytoplankton community composition was analyzed in sedimentation chambers according to standard techniques56,57 and counted using an inverted microscope (Nikon TMS). Phytoplankton were identified to the genus or species level and cell measurements were used to calculate biovolume and biomass according to Edler58 and Olenina, et al.59. For the period 2011–2014, size class was not recorded for the phytoplankton counts, therefore the median biovolume and biomass estimates of the size classes was used for the calculations. Starting in 2015, individual size classes were used for the phytoplankton counts. Mesozooplankton were identified to the genus or species level using a stereomicroscope (see table S1) and weights were calculated according to Hernroth60, HELCOM61 and Sprung62. Biomass was estimated from individual wet weight and an assumed carbon content of 5% of the wet weight63. As the mesh size of the plankton net was 200 µm, mainly large mesozooplankton are included in the analysis and these were classified into the groups: Copepods, Cladocera, Appendicularia, Rotifers and Other zooplankton (including Thecostraca, Bivalvia, Malacostraca, Gastropoda, Clitellata, Ostracoda, Polychaeta, Thaliacea, Oligotrichea as well as specimens of the phyla Chaetognatha, Echinodermata and unidentified), .
All data handling, statistical analyses and graphics were performed using R version 3.6.064 and the “ggplot2” package11 in combination with the “gridExtra” package65, and the “tidyverse” package66. Spearman correlations were performed and visualized using the “psych”67 and “corrplot” packages68. Spearman's rank correlation coefficient is reported as ρ (rho). If not stated otherwise, values correspond to average values per seasons. Average values correspond to rolling averages covering 15 julian days, spanning all years. If sampling was carried out the same julian day in different years, a mean was calculated prior to calculating the rolling average. The aim of this was to dampen the impact of singular samplings and to acquire an overall pattern. Seasons in the Baltic Proper were defined according to HELCOM as follows; Spring: March–May, Summer: June–September, Autumn: October–December, and Winter: January–February69. Bray–Curtis similarity (1- dissimilarity) over time was analyzed using the “vegan” package70, using “genus” level in all trophic levels. Non-metric multidimensional scaling (NMDS) analysis was performed using the “metaMDS()” function in the “vegan” package70 using median community composition to calculate Bray–Curtis dissimilarity to visualize differences in plankton community composition over the study period. Nutrient availability was calculated as described by Fleming and Kaitala71, with the equation \(\sqrt[3]{\left[{NO}_{3}^{-}+ {NO}_{2}^{-}\right]\times \left[{PO}_{4}^{3-}\right]\times \left[{SiO}_{4}^{4-}\right]}\) where [NO3− + NO2−], [PO43−] and [SiO44−] are the concentrations of the respective nutrients, to define the combined nutrient level. Mixed layer depth was calculated using the “oce” package72, which uses in situ pressure, temperature and salinity to calculate the density (sigma-theta, σθ). Mixed layer depth was defined as the depth where Δσθ > 0.125 kg m−3, compared to the surface density73. Over the sampling time span, 13 sampling occasions had a mixed layer depth between 0 and 1 m, which could be caused by the stabilization time being too short. However, all samplings are presented as this was considered to have limited effects on the interpretation of our results. In addition, STRÅNG data (sunshine duration and PAR) were extracted from the Swedish Meteorological and Hydrological Institute (SMHI) and produced with support from the Swedish Radiation Protection Authority and the Swedish Environmental Agency. Spearman correlations to investigate relationships among biotic and abiotic parameters used the complete dataset with individual samplings. To investigate relationships among bacterio-, phyto- and the large mesozooplankton community composition, Spearman correlations for relative biomasses of bacterio-, phyto- and large mesozooplankton were used and in addition to the five taxonomical groups of large mesozooplankton presented in Fig. 1. Litostomatea (principally Mesodinium rubrum) is a mixotroph74,75 but was here considered a predator in the analysis. The output from the analysis was visualized using the “circlize” package76 in combination with the “ComplexHeatmap” package77. ‘Other zooplankton’ included specimens from the classes Thecostraca, Bivalvia, Malacostraca, Gastropoda, Clitellata, Ostracoda, Polychaeta, Thaliacea, Oligotrichea as well as specimens of the phyla Chaetognatha, Echinodermata and unidentified. Contribution to the relative biomass was always < 2% for each separate class/phyla included in ‘Other zooplankton’.
Annual variation of sea-surface temperature (a), sunshine duration (b) and mixed layer depth (sigma theta 0.125) (c), salinity (d), nutrient availability (e), chlorophyll a concentration (f), levels of dissolved organic carbon, DOC (g), and colored dissolved organic matter, cDOM (h). Grey-shaded area illustrate moving standard deviation and the black line shows the moving average covering 15 julian days, combining all the years. Points illustrate separate samplings over the years with a different color for each year. Vertical dashed lines illustrate breaks for seasons.
The extensive sampling at LMO revealed strong seasonality for all investigated parameters. In general, surface (2 m) water temperature was stable around 4ºC during winter (January–February, days 0–59), followed by a steady increase from mid-spring (March–May, days 60–151), peaking at ~ 18ºC during summer (June–September, days 152–273), and a decline during autumn (October-December, 274–365) (Fig. 1a). Temperature displayed some interannual variability, e.g. in the range of 4–6 °C among years during summer, but with similar seasonality patterns among the years (Fig. 1a; Fig. S1a). Number of sunshine hours ranged from one to three hours per day in winter to ten hours of sunshine per day in summer during 2011–2019 (Fig. 1b). During winter and autumn, the interannual variation was lower, with ± 1 h, whereas summer values had a larger range among the years with ± 2.5 h (Fig. 1b; Fig. S1b). The water-column was well mixed during winter, until late spring (May) when it started to stratify by temperature which was maintained during summer, followed by well mixed conditions again during autumn (Fig. 1c). Mixed layer depth ranged from 20 to 25 m during winter and autumn but decreased during summer to ~ 10 m. There was a high variability within and among the years for the different seasons, with mixed layer depths ranging from 1 to 38 m during winter, whilst summer values were less variable ranging from 1 to 27 m. In autumn, mixed layer depth varied between 9 and 39 m over the study period (Fig. 1c; Fig. S1c). Salinity was less variable among seasons but was highest during winter and spring and lowest during summer (Fig. 1d) and exhibited some interannual variability (Fig. 1d; Fig. S1d).
The combined nutrient availability (\(\sqrt[3]{\left[{NO}_{3}^{-}+ {NO}_{2}^{-}\right]\times \left[{PO}_{4}^{3-}\right]\times \left[{SiO}_{4}^{4-}\right]}\)) was high during winter, with a rapid decrease during April (days 91–120) and then remained at a constant level from mid-spring to mid-autumn (October, days 274–304) (Fig. 1e). Nutrient availability showed some interannual variability, with recurring peaks autumn/winter and valleys during summer (Fig. 1e; Fig. S1e). Nitrate displayed the most pronounced seasonality pattern (Fig. S1a) ranging from ~ 2.5 to 0.2 µM. Phosphate showed similar seasonality, ranging from concentrations of 0.8 to 0.1 µM (Fig.S2c). Silicate concentrations ranged from 16 to 8 µM and tended to be higher during winter (Fig. S2d). Ammonium varied between 0.4 and 1.2 with no systematic patterns over time (Fig. S2b). Interannual variability for all nutrients is available in the supplementary files (Fig. S2; Fig. S3). DOC concentrations were at a relatively constant level during autumn, winter, and spring (~ 360 µM) but accumulated after mid-summer (July, days 182–212) reaching up to 410 µM. Towards the end of summer (September, days 244–273), DOC concentrations decreased to levels prior to mid-summer at around 360 µM (Fig. 1g). There was an increase of DOC concentrations in the years 2017 and 2018 which increased the interannual variability (Fig. 1g, Fig. S1g). cDOM was lowest during winter, ~ 0.25 followed by an increase by mid-spring (April, days 91–120) at around 0.30 followed by a decline. A larger and more pronounced peak (~ 0.33) developed around July and lasted until the end of August (days 182–243), by which cDOM decreased to winter values (Fig. 1h). Except for two samplings (spring 2018 and summer 2013), cDOM displayed a low interannual variability (Fig. 1h; Fig. S1h).
Bacterioplankton biomass generally started to increase during early spring, reaching a plateau around 40 mg C m−3 until July (days 182–212) when it increased further to reach around 90 mg C m−3 towards the middle of the summer (Fig. 2a). During autumn, bacterioplankton biomass decreased steadily down to winter values at around 20 mg C m−3. The trend was present for all the years, but the bacterioplankton biomass was greater for the period 2011–2012 compared to 2013–2018. However, within the two separate periods, the interannual variability was less pronounced (Fig. 2a; Fig. S4a). Chl a peaked in April (days 91–120) during the spring bloom (~ 5 µg L−1), after which it decreased to close to pre-bloom levels (~ 1 µg L−1). By mid-summer (July, days 182–212)), a summer peak (~ 3 µg L−1) was followed by a slow decrease during autumn to winter levels (~ 0.5 µg L−1) (Fig. 1f). Interannual variability was low, with only two samplings (spring 2011 and 2018) having notably high values > 10 µg L−1 (Fig. 1f; Fig. S1f.). Phytoplankton biomass displayed a similar pattern compared to Chl a, with a spring bloom in April (days 91–120), reaching ~ 250 mg C m−3. During late spring, phytoplankton biomass had decreased to around 100 mg C m−3 (Fig. 2c) and a second peak was reached by mid-summer (July, 182–212) at ~ 200 mg C m−3. Subsequently, biomass stabilized around 100 mg C m−3 and decreased to winter values in December-January (25–50 mg C m−3). Phytoplankton were quantified using slightly different analytical methods during 2011–2014 compared to 2015–2018 (see material and methods) indicating that separate analysis of the two time periods is necessary. Biomass did not display a large interannual variability during the two periods 2011–2014 and 2015–2018 when analyzed separately. However, biomass was higher throughout the first period compared to the later (Fig. S4b). This result was probably related to the differences in the analytical procedure for the two time periods. Total biomass of large mesozooplankton displayed a seasonal trend with winter biomass of ~ 1 mg C m−3, followed by a rapid increase in mid-spring (May, days 121–151) peaking during early summer (June, days 152–181)) with large mesozooplankton biomass of ~ 10 mg C m−3 being stable for about a month (Fig. 2e). Towards the end of July (day 212), large mesozooplankton started to decrease and by autumn reached ~ 5 mg C m−3. The biomass of large mesozooplankton displayed low interannual variability, with similar seasonality among the years (Fig. 2e; Fig. S4c).
Annual trend of bacterioplankton (a), phytoplankton (c), large mesozooplankton biomass (e), and relative bacterioplankton (b) phytoplankton (d) and large mesozooplankton community composition (f). For (b, d and f), the ridge plot shows the moving average covering 15 julian days, combining all the years whilst bars show discrete days where the coverage is partial. Points illustrate separate samplings over the years with a different color for each year and the black lines show the moving average covering 15 julian days, combining all years. Vertical dashed lines illustrate breaks for seasons. Note that cyanobacteria community composition is included in the phytoplankton category (d) and not in the bacterioplankton category (b) in this figure (see methods).
Analysis of 16S rRNA gene amplicons showed that Actinobacteria, Bacteroidetes, Proteobacteria were abundant along with Planctomycetes and Verrucomicrobia, and several of these phyla exhibited pronounced seasonal changes in relative abundance (Fig. 2b, Fig. 6). For example, Actinobacteria exhibited two periods of elevated relative abundance (in spring and autumn), while Verrucomicrobia peaked in summer. Bacterial community composition did not vary largely interannually but displayed recurring patterns (Fig. S5) that were similar over the study period (Fig. S8a). Microscopy analysis showed that the phytoplankton and mixotrophic communities during spring were dominated by Dinophyceae (47%), Litostomatea (ciliates, 36%), Bacillariophyceae (12%) and Flagellates (11%) (Fig. 2d, Fig. 6). During summer, Cyanophyceae (38%) and Litostomatea (35%) made up the largest proportion of the community, followed by Flagellates (23%) and Dinophyceae (11%). During autumn, Litostomatea dominated in terms of biomass, accounting for 59% of the phytoplankton community. Bacillariophyceae (20%), as well as Flagellates (16%) comprised a larger part of the community. The winter phytoplankton community had a large contribution of Litostomatea (41%), but Dinophyceae (27%), Cyanophyceae (18%) and Flagellates (11%) were also present. When investigating the periods 2011–2014 and 2015–2018 separately, phytoplankton community composition did not show a large interannual variability. However, comparing the two periods show a clear large difference (Fig. S4b; Fig. S6; Fig. S8b), related to the relative contribution of Litostomatea, Flagellates, Cryptophyceae and “Other phytoplankton” and this was probably mostly due to differences in analytical procedure and should not be interpreted as an ecological effect. In the large mesozooplankton fraction, Copepoda (63%), Cladocera (20%) and Appendicularia (16%) contributed the most to the total biomass during spring (Fig. 2f, Fig. 6, Fig. S7). Copepoda was the largest group in the community during summer (57%), when also Cladocera peaked (40%). Similarly, during autumn, Copepoda was the main group (71%), followed by Cladocera (28%) and Appendicularia (9%). Winter community’s main contributors were Copepoda (74%), Appendicularia (27%) and ‘Other zooplankton’ (21%) (Fig. 2f). The large mesozooplankton community composition showed similar recurring patterns on an annual scale (Fig. S7) and was similar over the study period (Fig. S8c). The most common bacterio-, phyto- and large mesozooplankton taxa are specified in Table S1.
Temporal analysis of Bray–Curtis similarity values of the respective community compositions between samplings (indicative of temporal community similarity/dissimilarity) displayed a sinusoidal annual pattern for all trophic levels (Fig. 3). The seasonality pattern was more pronounced for the large mesozooplankton. Moreover, for the period when large mesozooplankton composition data was available (5 years), the dynamics in similarity values were highly coordinated with distinct peaks within only a month’s time. However, the data for large mesozooplankton was only collected once per month reducing the resolution of the seasonal cycle in these taxa. Bacterioplankton and large mesozooplankton Bray–Curtis similarity mostly varied around 0.25, whereas similarity values for the phytoplankton decreased from around 0.45 to 0.1 over periods of up to 3 years after which it stabilized (Fig. 3).
Bray–Curtis similarity plot for bacterioplankton, phytoplankton and large mesozooplankton. Dots mark monthly averages of pairwise comparisons of samplings “n” months apart. Therefore, the first point shows the average similarity of all pairs of samples collected within one month of each other, the second point for all samples collected within 2 months of each other. Vertical dashed lines illustrate breaks for years (365 days). Note that cyanobacteria community composition is only included in the phytoplankton category and not in the bacterioplankton category in this figure (see methods).
Bacterioplankton biomass displayed significant correlations with most of the investigated variables, except for DOC and sunshine duration (Fig. 4). Bacterioplankton biomass showed positive correlations with phytoplankton biomass, large mesozooplankton biomass, temperature, Chl a, cDOM and solar irradiance whilst the relationship was significantly negative with salinity, nutrient availability, and mixed layer depth. Phytoplankton biomass only displayed positive correlations with bacterioplankton biomass, Chl a, cDOM, and solar irradiance. Mesozooplankton biomass was positively correlated with bacterioplankton biomass and temperature and showed negative correlation with nutrient availability. Moreover, as could be expected, the abiotic variables showed multiple correlative relationships, among which temperature and nutrient availability displayed the largest number of correlations. Both temperature and nutrient availability displayed significant relationships with all investigated variables, except for phytoplankton biomass (Fig. 4).
Spearman correlation matrix for biotic and abiotic parameters. Color intensity and size of circle corresponds to correlation value. Values in bold italics in diagonal show the total number of datapoints available for each parameter. Values below diagonal display number of pairwise comparisons for the separate correlations. Blank squares show when the correlation significance level was p > 0.05.
Relative biomasses of bacterio-, phyto- and large mesozooplankton were used to infer potential interactions in the planktonic food web (Fig. 5), although one should note that correlations could also appear without any direct or indirect interaction. Here, we considered large mesozooplankton and Litostomatea (predominantly the mixotrophic ciliate Mesodinium rubrum) as predators whilst bacterioplankton and phytoplankton groups were prey items. The relative biomass of Copepoda was significantly correlated with the relative biomasses for five out of the six phytoplankton groups included, as well as four out of the 11 bacterioplankton groups (Fig. 5a). Strongest correlations were obtained with Cyanophyceae and ‘Other phytoplankton’ as well as Epsilonbacteraeota. Cladocera displayed significant correlations with four phytoplankton groups, and six bacterioplankton groups. The relative biomass of Rotifera was negatively correlated with two phytoplankton groups, and four bacterioplankton groups (Fig. 5b). Appendicularia had the highest number of significant correlations, with all phytoplankton groups, and nine bacterioplankton groups. Litostomatea (ciliates) displayed significant correlations with four phytoplankton groups, and three bacterioplankton groups (Fig. 5c). ‘Other zooplankton’ only displayed one significant correlation (Cryptophyceae). From the bacterioplankton community, relative biomasses of Euryarchaeota did not show any significant correlations with predators. Note that some care should be applied when interpreting the correlations with phytoplankton since they were estimated using slightly different analytical methods during 2011–2014 compared to 2015–2018. Furthermore, Mesodinium rubrum also has an autotropic life style74,75 and the correlations found here should be evaluated in light of the mixotrohic capacity of this species.
Chord diagram visualizing Spearman correlations of the relative biomass among predators (large mesozooplankton and Litostomatea) and prey (phyto- and bacterioplankton). Copepoda and Cladocera (a), Appendicularia and Rotifera (b), Litostomatea and ‘Other zooplankton’ (c) versus prey item as well as correlations between phytoplankton and bacterioplankton (d). Width and color of links are proportional to correlation coefficient. Only significant (p < 0.05) correlations are displayed. Links among predators and phytoplankton are presented with black border for distinction.
Next, we determined potential interactions between organisms in the microbial food web that compete for resources such as nutrients and carbon sources and in some cases light for energy production. Main phytoplankton groups in spring, such as Dinophyceae, showed a strong positive correlation (co-occurrence) with the bacteria Bacteroidetes and negative correlations (different seasonal pattern) with Verrucomicrobia and Planctomycetes (Fig. 5d). Additionally, the most common phytoplankton in summer, Cyanophyceae, had positive correlations with Verrucomicrobia and Planctomycetes and Firmicutes but negative correlations with the remaining groups excluding Euryarchaeota and ‘Minor phyla’.
When quantifying the carbon content in various plankton groups in different seasons, both inter- and intra-seasonal dynamics were present. Bacterioplankton carbon was highest during summer, but as the total carbon content also peaked during summer, the relative contribution from bacterioplankton during summer was 30% of seston carbon content (Fig. 6; Table S2). The relative contribution from bacterioplankton peaked during winter (34%), even though the amount of carbon was the lowest. Bacterioplankton relative contribution to the total carbon in the plankton food web was lowest during spring, and during autumn bacterioplankton accounted for 25% of the total carbon (Fig. 6; Table S2). Phytoplankton consistently accounted for > 39% of the total carbon in the plankton food web. During spring, phytoplankton carbon and its relative contribution was highest (65%), whilst reaching minimal values in winter (39%). During summer and autumn, phytoplankton carbon was at a similar level, considering both amount and relative contribution (Fig. 6; Table S2). Litostomatea (exclusively the ciliate M. rubrum) made a noteworthy contribution to total carbon in the plankton communities, comprising > 20% of the total carbon throughout all seasons. During winter, Litostomatea had the lowest carbon levels, whilst the highest levels were found during summer. For most of the seasons, large mesozooplankton carbon and relative contribution was low (1–3%), except during summer when values peaked (7%). During winter and spring, large mesozooplankton carbon was at its lowest. Autumn large mesozooplankton carbon was two times higher than winter and spring levels (Fig. 6; Table S2).
Seasonal carbon content in various plankton groups. Only samplings when all trophic levels were available are included (n = 46). The size of circles corresponds to carbon content (mg C m−3), and relative contribution to the total carbon pool is presented inside or above each circle. Values in bold italics are seasonal total carbon content for the different compartments of the planktonic food web studied here (i.e. bacterioplankton, phytoplankton, Litosomatea (at 2 m) and large mesozooplankton from water column from 30 m up to surface)). Drawings illustrate the major taxa in each plankton group and its size illustrate relative contribution and correspond to drawings in Fig. 5. Data presented in Table S2.
This multi trophic-level and multi-year study showed seasonal patterns in all trophic levels and a strong connectedness between abiotic parameters, especially with temperature and nutrient availability. Furthermore, in the plankton communities, potential interactions were identified among groups in all trophic levels. These data can be used in a continued work to understand food web structure and interactions of bacterio-, phyto- and large mesozooplankton communities and to provide new knowledge on the interdependencies of trophic levels at the crossroads of the microbial and traditional food webs. In our study, focusing on bacterio-, phyto- and large mesozooplankton, combining biomass, community composition and co-occurrences, an extra noteworthy observation was that phytoplankton made up most of the plankton carbon biomass throughout all seasons. However, it should be emphasized that there are also other compartments of the microbial food web which are not included here which can contribute a significant amount of carbon to the system, e.g. picophytoplankton, viruses, fungi, archaea and smaller zooplankton. When investigating temporal community similarity some striking features were revealed; as expected, but rarely shown, Bray–Curtis similarity was coordinated over time and among all plankton groups studied. This is in line with other studies on specific taxa of plankton suggesting similar patterns of recurring community similarity8,9,10,78,79.This implies that whilst the community composition changes over the seasons, it returns to a stable state over a year. This is a prerequisite to create biomass which was, relatively stable over the seasons, except in winter when the planktonic carbon pool was at its lowest.
Temperature has been reported to alter bacterial community composition80,81, as well as having a positive relationship with bacterial biomass and production34,40,82,83. In the present study, bacterioplankton biomass peaked just after sea surface temperature had reached its maximum and the two parameters were strongly correlated. Phytoplankton biomass did not display any correlation with temperature, in contrast to the unimodal relationship described by Legrand, et al.34 as well as the positive relationship with filamentous cyanobacteria84 and other phytoplankton phyla40. Zooplankton biomass shows a strong positive correlation with temperature85. In a recent study, mixotrophic and autotrophic biomass was found to correlate positively with temperature, whereas the biomass of heterotrophs did not show this relationship86. This illustrates that while being a largely influential factor, temperature is not the sole controlling parameter for plankton communities in the Baltic Sea. The interannual and within-year variation in temperature, mixed layer depth, salinity, and nutrient availability also illustrates that the sampled station is dynamic regarding the abiotic conditions.
Nutrient availability was inversely correlated to several abiotic variables as well as to bacterio- and large mesozooplankton biomass. Phytoplankton biomass was not correlated to nutrient availability, which is opposite to previous findings, e.g.87. This discrepancy could be due to the bimodal phytoplankton dynamics in the Baltic Proper, having one spring bloom comprised of Dinophyceae, Litostomatea and Bacillariophyceae followed by a summer bloom dominated by Cyanophyceae and Litostomatea. For example, decoupling between nutrient levels and cyanobacterial abundances has previously been reported84,88. Filamentous, heterocystous N2 fixers dominated the summer bloom, and these could fuel the system with bioavailable nitrogen for hetero- and autotrophic picoplankton89, as well as grazers through ingestion and other phytoplankton and microbes by exudates90.
Bacterioplankton community composition displayed a seasonal pattern, and the Bray–Curtis similarity was stable over the years. This pattern indicate that the composition of the communities follow a seasonal pattern and that the communities are more similar during the same season (even when several years apart), but the community composition is most different when comparing opposite seasons (~ 6 months apart)91. A high Bray–Curtis similarity over time implies temporal stability, concurrently a low similarity implies changes in the communities91. The dominance of filamentous Aphanizomenon flos-aquae and Nodularia spumigena in the microscopy counts during summer was consistent with previous reports from the Baltic Proper34,69,84,92,93. The phytoplankton and mixotrophic community showed a pronounced seasonality, with Litostomatea, Dinophyceae, Bacillariophyceae, and Cyanophyceae all being major contributors at different time points on an annual scale, which has been shown partially previously28,69. A change in the spring community composition from a pelagic system where diatoms dominated spring and autumn blooms, and cyanobacteria and green algae prevailed during summer blooms, to a revised picture with Dinophyceae, rather than Bacillariophyceae, being the dominant phytoplankton in spring surface waters has been reported previously34,69,93,94 and was also true for the present study. Although the large mesozooplankton community was sampled at a lower frequency than bacterio- and phytoplankton, data covering five years still allowed for exploring the community composition. The common paradigm in zooplankton ecology is that Copepoda is the dominating large mesozooplankton taxa in the Baltic Sea35,95. Copepoda were indeed the major taxa in the large mesozooplankton community during all seasons, consistently accounting for > 55% of the community (average per season). Yet, other taxa such as Appendicularia appear to have a significant role in the Baltic Sea in relation to copepods and cladocerans. Small-sized zooplankton were not included in this study, due to the sampling technique, making inferences about them in terms of abundance and biomass in relation to large mesozooplankton impossible. However, other studies have shown that small-sized zooplankton, such as rotifers and small cladocerans, play a critical role, both in terms of biomass and abundance, in the Baltic food web8,10. Furthermore, ciliates (Litostomatea; M. rubrum) also contributed to a large fraction of the carbon pool throughout the year, although some studies suggest that their life-style is dominated by autotrophy74,75, questioning their importance as predators. This suggests that a diverse and dynamic plankton community should be considered when assessing the food web structure, which could help explain the interactions in the planktonic realm.
Bacterioplankton community composition was consistent over the study period and displayed recurring seasonality throughout the 8 years. Bacterioplankton communities were more similar when comparing seasons among years, than when comparing different seasons within a year. The largest differences in seasonal community composition were between summer and winter and was related to the relative contribution of Nitrospinae, Chloroflexi, Epsilonbacteraeota and Minor phyla during winter. In contrast, phytoplankton community composition displayed two periods of varying community composition. Analytical methods for microscopic phytoplankton counts were different in 2011–2014 and 2015–2018, respectively. Hence, patterns over time are not relevant to study. Within the separate periods (2011–2014 and 2015–2018), phytoplankton community showed a low variability, where relative contribution of Cyanophyceae and Dinophyceae had the largest impact on the seasonal distinction of the community. Mesozooplankton community composition was similar over the years of the study and was dominated by Copepoda throughout. Relative increases of Cladocera and Rotifera during summer in conjunction with relative increases of Appendicularia during winter had an impact on the community composition over seasons, but not among years.
The strong seasonality observed in the phyto- and large mesozooplankton communities indicated seasonally reoccurring patterns in community structures and that the turnover was similar for all investigated plankton communities. Recently, Yeh and Fuhrman32 found different diversity patterns for protists and bacterioplankton at the Californian coast with more pronounced changes in eukaryotic communities, using amplicon sequencing. In contrast, our microscopic counts indicated a strong seasonality in both phytoplankton and large mesozooplankton communities compared to bacterioplankton. When combining all trophic levels, carbon biomass was highest during summer, followed by spring, autumn, and lowest during winter, similar to findings from lakes and river plumes in temperate regions26,27. Also, the relative contribution from various plankton groups matched previous reports26,27, where phytoplankton contributed most to the planktonic total carbon pool, followed by bacterioplankton and Litostomatea, whilst large mesozooplankton only had a minor contribution to the carbon pool.
Interactions between organisms include e.g., grazing, symbiosis, competition for resources, and parasitism. As such, a negative correlation can be interpreted in several ways, including a strong top-down effect (predation) or occurrence in different seasons (e.g., spring or summer adapted taxa), alternatively organisms may occupy the habitat during different periods of the year. If predators and prey are influenced by similar environmental drivers without strong top-down effects, they will co-occur (positive correlation). Grazing might result in the release of substrates (sloppy feeding) which provides new resources for organisms, and bulk plankton analyses might also capture pathogenic microbiomes, which will lead to positive correlations between species. Also, trophic effects are possible e.g., when macro- and mesozooplankton graze on microzooplankton and thus release grazing pressure on other planktonic groups (positive and negative correlations of biomass between separate groups). However, one should always keep in mind that significant correlations could arise without any interaction between the organisms whatsoever. Hence, correlation analyses cannot offer a definitive answer on effects and causation, but it enables investigations for potential links. A general caveat in this study is also that small sized zooplankton are not included in the sampling which may underestimate the correlations between primary producers and consumers.
In our study, Copepoda correlated mainly negatively with various phytoplankton taxa, possibly due to grazing. However, phytoplankton which are not vulnerable to grazing (for instance due to morphology) could be negatively correlated to Copepoda not only due to grazing but due to increased competition over resources from other phytoplankton or bacterial species which are positively correlated with Copepoda, e.g., Cyanophyceae and Flagellates. Copepoda correlated mainly positively with various bacterioplankton taxa, which could be the result of sloppy feeding. However, it should be noted that zooplankton were sampled in the top 30 m of the sampling station whereas bacterio- and phytoplankton were sampled at two meters. Zooplankton are known to have both seasonal and diel habitat partitioning. For example, many cladocera are mostly abundant in summer in warm surface water, whereas many appendicularians predominantly avoid the warm surface water8,9,10. Hence, some zooplankton taxa in our data set may have a different niche compared to the bacterio- and phytoplankton we sampled in surface water (2 m), and the correlations observed may not be a function of true direct or indirect interactions.
Furthermore, zooplankton have been shown to interact with bacterioplankton in more ways than predation, both by attracting and enabling growth in the zoosphere, as well as farming them inside or on the outside of the body97. In the present study, relative biomass of Litostomatea had no significant correlations with other predators. Ciliates have been reported to forage successfully on phytoplankton, occasionally at rates more than double that of zooplankton98. Litostomatea was negatively correlated with Dinophyceae (predation and/or competition) and positively correlated with Cryptophyceae, Bacillariophyceae and ‘Other phytoplankton’ and showed mainly positive relationships with bacterioplankton, suggesting that Litostomatea does not exert a substantial feeding pressure on phytoplankton and bacterioplankton.
Relationships among phyto- and bacterioplankton have been studied in several systems and experiments, but studies covering longer time periods and with information on community composition are sparse99,100,101,102,103. However, Yeh and Fuhrman32 recently described how prokaryotic and protist communities differ in diversity over both depths and time. Considering the current findings and previous knowledge, plankton biomass and community composition is associated to shifts in e.g., temperature and nutrient availability. In the present study, several positive and negative correlations among phyto- and bacterioplankton were present, further suggesting functional couplings, some potentially being top-down control and some bottom-up control. Previous studies have found that Alphaproteobacteria and Bacteroidetes were associated to phytoplankton spring blooms, whereas the summer bloom was more associated with Alpha- and Gammaproteobacteria and Bacteroidetes104,105,106. Also, filamentous cyanobacteria had tight couplings to Bacteroidetes, Gammaproteobacteria and Verrucomicrobia107, which is in line with the finding in our study. Furthermore, species richness and diversity of particle-attached and free-living bacteria have been shown to have varying responses to phytoplankton blooms108. This shows the connectivity among trophic levels but also the complexity when interpreting multi-level associations. In the broader context there have been large changes occurring in the Baltic Sea during the twentieth century. Anthropogenic pressures such as climate change, eutrophication, and over-exploitation, together with climate variability, have led to regime shifts4,109. Recent studies suggest that the Central Baltic Sea is in a pelagic dominated, high productivity state with, e.g., low abundances of cod and copepods such as Pseudocalanus and high phytoplankton biomass during summer110. Hence, this dataset should be interpreted in this broader context and collection of time series data is important to detect future possible regime shifts.
This study provides insights into the structure and seasonal dynamics of different trophic levels of the pelagic microbial food web and how the trophic levels correlate with each other and environmental drivers. At present, studies covering several years, seasons, trophic levels, and environmental drivers are scarce. As such, our study offers a rare window into seasonal succession, food web structure, and interactions of bacterio-, phyto-, and large mesozooplankton communities. Our results showed a clear and stable seasonal succession of the investigated plankton communities. For phyto- and large mesozooplankton there was a pronounced reoccurring pattern with highest similarity every 12 months. Also, we found a strong interconnectivity between bacterioplankton and other parts of the aquatic environment. Temperature and nutrient availability correlated with all variables except for phytoplankton biomass, illustrating that these two environmental drivers are crucial for shaping the pelagic food web. Plankton community composition displayed a low interannual variability and the different seasons within a year were more dissimilar than matching seasons among years.
This study contributes a baseline of the food web structure in the plankton realm and here Dinophyceae, not Bacillariophyceae, dominated the spring phytoplankton bloom in the Baltic Proper between 2011 and 2018. Another striking result is that Litostomatea (ciliates) and Appendicularia contribute to the food web in such a large extent. Important members in the phytoplankton and bacterioplankton community were Cyanophyceae, Actinobacteria, Bacteroidetes, and Proteobacteria. Considering the recognition of the importance of interspecies interactions for community dynamics and biogeochemical processes111,112,113, time series research should include all trophic levels for a holistic understanding of the various food web interactions. Future studies should continue to determine the nature of the interactions using more mechanistic studies to determine the nature of direct predation, competition and mutualistic effects as well as including understudied groups such as picophytoplankton, viruses, fungi, archaea and microzooplankton.
Data will be available upon request. 16S data are deposited in the EMBL-EBI European Nucleotide Archive repository (https://www.ebi.ac.uk/ena), accession numbers PRJEB42455, SRP048666, PRJEB52855, PRJEB52782, PRJEB52780, PRJEB52772, PRJEB52627, PRJEB52496, PRJEB52828, PRJEB52837 and PRJEB52854.
Snoeijs-Lejonmalm, P., Schubert, H. & Radziejewska, T. Biological Oceanography of the Baltic Sea (Springer, 2017).
Book Google Scholar
Hyytiäinen, K. et al. Provision of aquatic ecosystem services as a consequence of societal changes: The case of the Baltic Sea. Popul. Ecol. 63, 61–74. https://doi.org/10.1002/1438-390X.12033 (2021).
Article Google Scholar
Casini, M. et al. Multi-level trophic cascades in a heavily exploited open marine ecosystem. Proc. R. Soc. B 275, 1793–1801. https://doi.org/10.1098/rspb.2007.1752 (2008).
Article PubMed PubMed Central Google Scholar
Möllmann, C. & Diekmann, R. in Global Change in Multispecies Systems: Part II Vol. 47 Advances in Ecological Research (eds G. Woodward, J. Ute, & E. J. O´Gorman) Ch. 4, 303–347 (Academic Press, 2012).
Eddy, T. D. et al. Energy flow through marine ecosystems: confronting transfer efficiency. Trends. Ecol. Evol. 36, 76–86. https://doi.org/10.1016/j.tree.2020.09.006 (2021).
Article PubMed Google Scholar
Azam, F. Microbial control of oceanic carbon flux: The plot thickens. Science 280, 694–696. https://doi.org/10.1126/science.280.5364.694 (1998).
Article CAS Google Scholar
Fuhrman, J. A., Sleeter, T. D., Carlson, C. A. & Proctor, L. M. Dominance of bacterial biomass in the Sargasso Sea and its ecological implications. Mar. Ecol. Prog. Ser. 57, 207–217 (1989).
Article ADS Google Scholar
Hernroth, L., Ackefors, H. & Havsfiskelaboratoriet. The Zooplankton of the Baltic Proper: A Long-term Investigation of the Fauna, Its Biology and Ecology. (Institute of Marine Research, 1977)
Schulz, J., Möllmann, C. & Hirche, H.-J. Vertical zonation of the zooplankton community in the Central Baltic Sea in relation to hydrographic stratification as revealed by multivariate discriminant function and canonical analysis. J Mar. Syst 67, 47–58. https://doi.org/10.1016/j.jmarsys.2006.09.004 (2007).
Article Google Scholar
Schulz, J. et al. Spatial and temporal habitat partitioning by zooplankton in the Bornholm Basin (central Baltic Sea). Progress Oceanogr. 107, 3–30. https://doi.org/10.1016/j.pocean.2012.07.002 (2012).
Article ADS Google Scholar
Wickham, H. ggplot2: Elegant Graphics for Data Analysis (Springer-Verlag, 2009).
Book MATH Google Scholar
Bar-On, Y. M., Phillips, R. & Milo, R. The biomass distribution on earth. Proc. Natl. Acad. Sci. U.S.A. 115, 6506–6511. https://doi.org/10.1073/pnas.1711842115 (2018).
Article ADS CAS PubMed PubMed Central Google Scholar
Whitman, W. B., Coleman, D. C. & Wiebe, W. J. Prokaryotes: The unseen majority. Proc. Natl. Acad. Sci. U.S.A. 95, 6578–6583. https://doi.org/10.1073/pnas.95.12.6578 (1998).
Article ADS CAS PubMed PubMed Central Google Scholar
Azam, F. & Hodson, R. E. Size distribution and activity of marine microheterotrophs. Limnol. Oceanogr. 22, 492–501. https://doi.org/10.4319/lo.1977.22.3.0492 (1977).
Article ADS CAS Google Scholar
Fuhrman, J. A., Sleeter, T. D., Carlson, C. A. & Proctor, L. M. Dominance of bacterial biomass in the Sargasso Sea and its ecological implications. Mar. Ecol. Prog. Ser. 57, 207–217. https://doi.org/10.3354/meps057207 (1989).
Article ADS Google Scholar
Pomeroy LR (2001) Caught in the food web complexity made simple. Sci Mar 65: 31–40
Sañudo-Wilhelmy, S. A., Gómez-Consarnau, L., Suffridge, C. & Webb, E. A. The role of B vitamins in marine biogeochemistry. Annu. Rev. Mar. Sci. 6, 339–367. https://doi.org/10.1146/annurev-marine-120710-100912 (2014).
Article ADS Google Scholar
Novotny, A., Jan, K. M. G., Dierking, J. & Winder, M. Niche partitioning between planktivorous fish in the pelagic Baltic Sea assessed by DNA metabarcoding, qPCR and microscopy. Sci. Rep. 12, 10952. https://doi.org/10.1038/s41598-022-15116-7 (2022).
Article ADS CAS PubMed PubMed Central Google Scholar
Reissig, M., Trochine, C., Queimaliños, C., Balseiro, E. & Modenutti, B. Impact of fish introduction on planktonic food webs in lakes of the patagonian plateau. Biol. Conserv. 132, 437–447. https://doi.org/10.1016/j.biocon.2006.04.036 (2006).
Article Google Scholar
Behrenfeld, M. J. et al. Climate-driven trends in contemporary ocean productivity. Nature 444, 752–755. https://doi.org/10.1038/nature05317 (2006).
Article ADS CAS PubMed Google Scholar
Field, C. B., Behrenfeld, M. J., Randerson, J. T. & Falkowski, P. Primary production of the biosphere: Integrating terrestrial and oceanic components. Science 281, 237–240. https://doi.org/10.1126/science.281.5374.237 (1998).
Article ADS CAS PubMed Google Scholar
Ducklow, H. W. in Microbial Ecology of the Oceans (ed D. L. Kirchman) 85–120 (John Wiley and Sons, 2000).
Calbet, A. Mesozooplankton grazing effect on primary production: A global comparative analysisin marine ecosystems. Limnol. Oceanogr. 46, 1824–1830. https://doi.org/10.4319/lo.2001.46.7.1824 (2001).
Article ADS Google Scholar
Ejsmond, M. J. et al. Modeling vitamin B1 transfer to consumers in the aquatic food web. Sci. Rep. 9, 10045. https://doi.org/10.1038/s41598-019-46422-2 (2019).
Article ADS CAS PubMed PubMed Central Google Scholar
Chust, G. et al. Biomass changes and trophic amplification of plankton in a warmer ocean. Glob. Chang. Biol. 20, 2124–2139. https://doi.org/10.1111/gcb.12562 (2014).
Article ADS PubMed Google Scholar
Bode, A. et al. Seasonal variability of plankton blooms in the Ria de Ferrol (NW Spain): II. Plankton abundance, composition and biomass. Estuarine Coast. Shelf Sci. 63, 285–300. https://doi.org/10.1016/j.ecss.2004.11.021 (2005).
Article ADS Google Scholar
Hochstädter, S. Seasonal changes of C: P ratios of seston, bacteria, phytoplankton and zooplankton in a deep, mesotrophic lake. Freshwater Biol. 44, 453–463 (2000).
Article Google Scholar
Novotny, A., Zamora-Terol, S. & Winder, M. DNA metabarcoding reveals trophic niche diversity of micro and mesozooplankton species. Proc. R. Soc. B 288, 20210908. https://doi.org/10.1098/rspb.2021.0908 (2021).
Article CAS PubMed PubMed Central Google Scholar
Tiselius, P., Belgrano, A., Andersson, L. & Lindahl, O. Primary productivity in a coastal ecosystem: A trophic perspective on a long-term time series. J. Plankton. Res. 38, 1092–1102. https://doi.org/10.1093/plankt/fbv094 (2016).
Article CAS Google Scholar
Andersson, A., Hajdu, S., Haecky, P., Kuparinen, J. & Wikner, J. Succession and growth limitation of phytoplankton in the Gulf of Bothnia (Baltic Sea). Mar. Biol. 126, 791–801 (1996).
Article Google Scholar
Mayer, J. et al. Seasonal successions and trophic relations between phytoplankton, zooplankton, ciliate and bacteria in a hypertrophic shallow lake in Vienna, Austria. Hydrobiologia 343(343), 165–174 (1997).
Article Google Scholar
Yeh, Y.-C. & Fuhrman, J. A. Contrasting diversity patterns of prokaryotes and protists over time and depth at the San-Pedro ocean time series. ISME Commun. https://doi.org/10.1038/s43705-022-00121-8 (2022).
Article PubMed Central Google Scholar
Worden, A. Z. et al. Rethinking the marine carbon cycle: Factoring in the multifarious lifestyles of microbes. Science 347, 1257594. https://doi.org/10.1126/science.1257594 (2015).
Article CAS PubMed Google Scholar
Legrand, C. et al. Interannual variability of phyto-bacterioplankton biomass and production in coastal and offshore waters of the Baltic Sea. AMBIO 44, 427–438. https://doi.org/10.1007/s13280-015-0662-8 (2015).
Article CAS PubMed PubMed Central Google Scholar
Fridolfsson, E. et al. Seasonal variation and species-specific concentrations of the essential vitamin B1 (thiamin) in zooplankton and seston. Mar. Biol. 166, 70. https://doi.org/10.1007/s00227-019-3520-6 (2019).
Article CAS Google Scholar
Nilsson, E. et al. Genomic and seasonal variations among aquatic phages infecting the baltic sea gammaproteobacteria Rheinheimera sp. BAL341. Appl. Environ. Microbiol. https://doi.org/10.1128/AEM.01003-19 (2019).
Article PubMed PubMed Central Google Scholar
Paerl, R. W. et al. Prevalent reliance of bacterioplankton on exogenous vitamin B1 and precursor availability. PNAS 115, E10447–E10456. https://doi.org/10.1073/pnas.1806425115 (2018).
Article ADS CAS PubMed PubMed Central Google Scholar
Lindh, M. V. et al. Disentangling seasonal bacterioplankton population dynamics by high-frequency sampling. Environ. Microbiol. 17, 2459–2476. https://doi.org/10.1111/1462-2920.12720 (2014).
Article Google Scholar
Berner, C., Bertos-Fortis, M., Pinhassi, J. & Legrand, C. Response of microbial communities to changing climate conditions during summer cyanobacterial blooms in the baltic sea. Front. Microbiol. 9, 1–16. https://doi.org/10.3389/fmicb.2018.01562 (2018).
Article CAS Google Scholar
Bunse, C. et al. High frequency multi-year variability in baltic sea microbial plankton stocks and activities. Front. Microbiol. 9, 3296. https://doi.org/10.3389/fmicb.2018.03296 (2019).
Article PubMed PubMed Central Google Scholar
Martínez-Garcia, S. et al. Seasonal dynamics in carbon cycling of marine bacterioplankton are lifestyle dependent. Front. Microbiol. https://doi.org/10.3389/fmicb.2022.834675 (2022).
Article PubMed PubMed Central Google Scholar
Valderama, J. C. in Manual on Harmful Marine Microalgae. IOC Manuals and Guides, vol. 33, (eds G. M. Hallegraeff, D. M. Anderson, & A. D. Cembella) 251–268 (UNESCO, 1995).
Jespersen, A. M. & Christoffersen, K. Measurements of chlorophyll-a from phytoplankton using ethanol as extraction solvent. Arch. Hydrobiol. 109, 445–454 (1987).
Article CAS Google Scholar
Cauwet, G. in Methods of Seawater Analysis (ed K. Grashoff, Kremling, K., Ehrhardt, M. ) 407 – 420 (Wiley-VCH, 1999).
Traving, S. J. et al. The effect of increased loads of dissolved organic matter on estuarine microbial community composition and function. Front. Microbiol. https://doi.org/10.3389/fmicb.2017.00351 (2017).
Article PubMed PubMed Central Google Scholar
Stedmon, C. A., Markager, S. & Kaas, H. Optical properties and signatures of chromophoric dissolved organic matter (CDOM) in danish coastal waters. Estuarine Coast. Shelf Sci. 51, 267–278. https://doi.org/10.1006/ecss.2000.0645 (2000).
Article ADS CAS Google Scholar
Gasol, J. M. & Morán, X. A. G. in Hydrocarbon and Lipid Microbiology Protocols: Single-Cell and Single-Molecule Methods (eds Terry J. McGenity, Kenneth N. Timmis, & Balbina Nogales) 159–187 (Springer, 2016).
Lee, S. & Fuhrman, J. A. Relationships between biovolume and biomass of naturally derived marine bacterioplankton. Appl. Environ. Microbiol. 53, 1298–1303 (1987).
Article ADS CAS PubMed PubMed Central Google Scholar
Boström, K. H., Simu, K., Hagström, Å. & Riemann, L. Optimization of DNA extraction for quantitative marine bacterioplankton community analysis. Limnol. Oceanogr. Methods 2, 365–373. https://doi.org/10.4319/lom.2004.2.365 (2004).
Article Google Scholar
Bunse, C. et al. Response of marine bacterioplankton pH homeostasis gene expression to elevated CO2. Nat. Clim. Chang. 6, 483–487. https://doi.org/10.1038/nclimate2914 (2016).
Article ADS CAS Google Scholar
Herlemann, D. P. et al. Transitions in bacterial communities along the 2000 km salinity gradient of the Baltic Sea. ISME J. 5, 1571–1579. https://doi.org/10.1038/ismej.2011.41 (2011).
Article CAS PubMed PubMed Central Google Scholar
Hugerth, L. W. et al. DegePrime, a program for degenerate primer design for broad-taxonomic-range PCR in microbial ecology studies. Appl. Environ. Microbiol. 80, 5116–5123. https://doi.org/10.1128/AEM.01403-14 (2014).
Article ADS CAS PubMed PubMed Central Google Scholar
Straub, D. et al. Interpretations of environmental microbial community studies are biased by the selected 16S rRNA (Gene) amplicon sequencing pipeline. Front. Microbiol. https://doi.org/10.3389/fmicb.2020.550420 (2020).
Article PubMed PubMed Central Google Scholar
Callahan, B. J. et al. DADA2: High-resolution sample inference from Illumina amplicon data. Nat. Methods 13, 581–583. https://doi.org/10.1038/nmeth.3869 (2016).
Article CAS PubMed PubMed Central Google Scholar
Vetrovsky, T. & Baldrian, P. The variability of the 16S rRNA gene in bacterial genomes and its consequences for bacterial community analyses. PloS one 8, e57923. https://doi.org/10.1371/journal.pone.0057923 (2013).
Article ADS CAS PubMed PubMed Central Google Scholar
Utermöhl, H. Zur vervollkommnung der quantitativen phytoplankton-methodik. Mitt. Int. Ver. Theor. Angew. Limnol. 9, 1–38 (1958).
Google Scholar
Utermöhl, H. Neue wege in der quantitativen erfassung des planktons (mit besonderer berucksichtigung des ultraplanktons). Verh. Int. Ver. theor. Angew. Limnnol. 5, 567–596. https://doi.org/10.1080/03680770.1931.11898492 (1931).
Article Google Scholar
Edler, L. Recommendations on methods for marine biological studies in the Baltic sea: Phytoplankton and chlorophyll. Balt. Mar. Biol. Publ. 5, 1–38 (1979).
Google Scholar
Olenina, I. et al. Biovolumes and size-classes of phytoplankton in the Baltic Sea. HELCOM Balt. Sea Environ. Proc. 106, 144 (2006).
Google Scholar
Hernroth, L. Recommendations on methods for marine biological studies in the Baltic Sea: Mesozooplankton biomass assessment. Vol. 14 (Baltic Marine Environment Protection Commission, 1985).
ZEN-ZIIM zooplankton database (2021).
Sprung, M. Physiological energetics of mussel larvae (Mytilus edulis). I. Shell growth and biomass. Mar. Ecol. Prog. Ser. 17, 283–293 (1984).
Article ADS Google Scholar
Omori, M. & Ikeda, T. Methods in Marine Zooplankton Ecology Vol. 60 (Wiley Interscience, 1984).
Google Scholar
R: A language and environment for statistical computing. R Foundation for Statistical Computing (Vienna, Austria., 2019).
gridExtra: Miscellaneous Functions for "Grid" Graphics (2017).
Wickham, H. et al. Welcome to the tidyverse. J. Open Sour. Softw. 4, 1686. https://doi.org/10.21105/joss.01686 (2019).
Article ADS Google Scholar
psych: Procedures for Personality and Psychological Research v. 2.1.6 (Northwestern University, Evanston, Illinois, USA, 2021).
R package "corrplot": Visualization of a Correlation Matrix (Version 0.89) (2021).
Wasmund, N., Tuimala, J., Suikkanen, S., Vandepitte, L. & Kraberg, A. Long-term trends in phytoplankton composition in the western and central Baltic Sea. J. Mar. Syst. 87, 145–159. https://doi.org/10.1016/j.jmarsys.2011.03.010 (2011).
Article Google Scholar
vegan: Community Ecology Package v. R package version 2.5–7 (2020).
Fleming, V. & Kaitala, S. Phytoplankton spring bloom intensity index for the Baltic Sea estimated for the years 1992 to 2004. Hydrobiologia 554, 57–65. https://doi.org/10.1007/s10750-005-1006-7 (2006).
Article Google Scholar
oce: Analysis of Oceanographic Data (2019).
Monterey, G. & Levitus, S. Seasonal variability of mixed layer depth for the world ocean. 96 (NOAA, Washington D.C., 1997).
Altenburger, A. et al. Limits to the cellular control of sequestered cryptophyte prey in the marine ciliate Mesodinium rubrum. ISME J. 15, 1056–1072. https://doi.org/10.1038/s41396-020-00830-9 (2021).
Article CAS PubMed Google Scholar
Kim, M., Drumm, K., Daugbjerg, N. & Hansen, P. J. Dynamics of sequestered cryptophyte nuclei in mesodinium rubrum during starvation and refeeding. Front. Microbiol. 8, 423. https://doi.org/10.3389/fmicb.2017.00423 (2017).
Article PubMed PubMed Central Google Scholar
Gu, Z., Gu, L., Eils, R., Schlesner, M. & Brors, B. Circlize implements and enhances circular visualization in R. Bioinformatics 30, 2811–2812. https://doi.org/10.1093/bioinformatics/btu393 (2014).
Article CAS PubMed Google Scholar
Gu, Z., Eils, R. & Schlesner, M. Complex heatmaps reveal patterns and correlations in multidimensional genomic data. (2016).
Smetacek, V. Diatoms and the ocean carbon cycle. Protist 150, 25–32 (1999).
Article CAS PubMed Google Scholar
Suikkanen, S., Laamanen, M. & Huttunen, M. Long-term changes in summer phytoplankton communities of the open northern Baltic Sea. Estuarine Coast Shelf Sci. 71, 580–592. https://doi.org/10.1016/j.ecss.2006.09.004 (2007).
Article ADS Google Scholar
Lindh, M. V. et al. Consequences of increased temperature and acidification on bacterioplankton community composition during a mesocosm spring bloom in the Baltic Sea. Environ. Microbiol. Rep. 5, 252–262. https://doi.org/10.1111/1758-2229.12009 (2013).
Article CAS PubMed Google Scholar
Lambert, S. et al. Rhythmicity of coastal marine picoeukaryotes, bacteria and archaea despite irregular environmental perturbations. ISME J. 13, 388–401. https://doi.org/10.1038/s41396-018-0281-z (2019).
Article PubMed Google Scholar
Andersson, A. et al. Projected future climate change and Baltic Sea ecosystem management. AMBIO 44, 345–356. https://doi.org/10.1007/s13280-015-0654-8 (2015).
Article CAS PubMed PubMed Central Google Scholar
Morán, X. A. G. et al. More, smaller bacteria in response to ocean’s warming?. Proc. R. Soc. B. 282, 20150371. https://doi.org/10.1098/rspb.2015.0371 (2015).
Article CAS PubMed PubMed Central Google Scholar
Bertos-Fortis, M. et al. Unscrambling cyanobacteria community dynamics related to environmental factors. Front. Microbiol. 7, 625. https://doi.org/10.3389/fmicb.2016.00625 (2016).
Article PubMed PubMed Central Google Scholar
Richardson, A. J. In hot water: Zooplankton and climate change. ICES J. Mar. Sci. 65, 279–295 (2008).
Article Google Scholar
Romano, F., Pitta, P. & Dolan, J. Relationships of pelagic ciliates with the microbial food web components at a coastal station in the oligotrophic Eastern Mediterranean Sea: Temporal and vertical variability. J. Plankton. Res. https://doi.org/10.1093/plankt/fbab053 (2021).
Article PubMed PubMed Central Google Scholar
Mousing, E. A., Richardson, K. & Ellegaard, M. Global patterns in phytoplankton biomass and community size structure in relation to macronutrients in the open ocean. Limnol. Oceanogr. 63, 1298–1312. https://doi.org/10.1002/lno.10772 (2018).
Article ADS CAS Google Scholar
Munkes, B., Löptien, U. & Dietze, H. Cyanobacteria blooms in the Baltic Sea: A review of models and facts. Biogeosciences 18, 2347–2378. https://doi.org/10.5194/bg-18-2347-2021 (2021).
Article ADS CAS Google Scholar
Ploug, H. et al. Carbon, nitrogen and O2 fluxes associated with the cyanobacterium Nodularia spumigena in the Baltic sea. AMBIO 5, 1549–1558. https://doi.org/10.1038/ismej.2011.20 (2011).
Article CAS Google Scholar
Karlson, A. M. et al. Nitrogen fixation by cyanobacteria stimulates production in Baltic food webs. AMBIO 44(S3), 413–426. https://doi.org/10.1007/s13280-015-0660-x (2015).
Article CAS PubMed PubMed Central Google Scholar
Fuhrman, J. A., Cram, J. A. & Needham, D. M. Marine microbial community dynamics and their ecological interpretation. Nat. Rev. Microbiol. 13, 133–146. https://doi.org/10.1038/nrmicro3417 (2015).
Article CAS PubMed Google Scholar
Zufia, J. A., Legrand, C. & Farnelid, H. Seasonal dynamics in picocyanobacterial abundance and clade composition at coastal and offshore stations in the Baltic Sea. Sci. Rep. 12(1), 14330. https://doi.org/10.1038/s41598-022-18454-8 (2022).
Article ADS CAS PubMed PubMed Central Google Scholar
Wasmund, N. & Uhlig, S. Phytoplankton trends in the Baltic Sea. ICES J. Mar. Sci. 60, 177–186. https://doi.org/10.1016/s1054-3139(02)00280-1 (2003).
Article Google Scholar
Hjerne, O., Hajdu, S., Larsson, U., Downing, A. S. & Winder, M. Climate driven changes in timing, Composition and magnitude of the Baltic sea phytoplankton spring bloom. Front. Mar. Sci. https://doi.org/10.3389/fmars.2019.00482 (2019).
Article Google Scholar
Diaz-Gil, C. et al. Spatio-temporal composition and dynamics of zooplankton in the Kalmar Sound (western Baltic Sea) in 2009–2010. Boreal Environ. Res. 19, 323–335 (2014).
Google Scholar
Smetacek, V. The annual cycle of protozooplankton in the Kiel Bight. Mar. Biol. 63, 1–11. https://doi.org/10.1007/bf00394657 (1981).
Article Google Scholar
Shoemaker, K. M., Duhamel, S. & Moisander, P. H. Copepods promote bacterial community changes in surrounding seawater through farming and nutrient enrichment. Environ. Microbiol. 21, 3737–3750. https://doi.org/10.1111/1462-2920.14723 (2019).
Article CAS PubMed Google Scholar
Johansson, M., Gorokhova, E. & Larsson, U. Annual variability in ciliate community structure potential prey and predators in the open northern Baltic Sea proper. J. Plankton Res. 26, 67–80. https://doi.org/10.1093/plankt/fbg115 (2004).
Article Google Scholar
Steinberg, D. K. et al. Overview of the US JGOFS bermuda atlantic time-series study (BATS): A decade-scale look at ocean biology and biogeochemistry. Deep Sea Res. Part II Top. Stud. Oceanogr. 48, 1405–1447. https://doi.org/10.1016/S0967-0645(00)00148-X (2001).
Article ADS CAS Google Scholar
Karl, D. M. & Church, M. J. Microbial oceanography and the Hawaii ocean time-series programme. Nat. Rev. Microbiol. 12, 699–713. https://doi.org/10.1038/nrmicro3333 (2014).
Article CAS PubMed Google Scholar
Mestre, M., Höfer, J., Sala, M. M. & Gasol, J. M. Seasonal variation of bacterial diversity along the marine particulate matter continuum. Front. Microbiol. https://doi.org/10.3389/fmicb.2020.01590 (2020).
Article PubMed PubMed Central Google Scholar
Cram, J. A. et al. Seasonal and interannual variability of the marine bacterioplankton community throughout the water column over ten years. ISME J. 9, 563–580. https://doi.org/10.1038/ismej.2014.153 (2015).
Article PubMed Google Scholar
Lambert, S., Lozano, J. C., Bouget, F. Y. & Galand, P. E. Seasonal marine microorganisms change neighbours under contrasting environmental conditions. Environ. Microbiol. 23, 2592–2604. https://doi.org/10.1111/1462-2920.15482 (2021).
Article PubMed Google Scholar
Pinhassi, J. & Hagström, Å. Seasonal succession in marine bacterioplankton. Aquat. Microb. Ecol. 21, 245–256 (2000).
Article Google Scholar
Buchan, A., LeCleir, G. R., Gulvik, C. A. & Gonzalez, J. M. Master recyclers: features and functions of bacteria associated with phytoplankton blooms. Nat. Rev. Microbiol. 12, 686–698. https://doi.org/10.1038/nrmicro3326 (2014).
Article CAS PubMed Google Scholar
Teeling, H. et al. Recurring patterns in bacterioplankton dynamics during coastal spring algae blooms. eLife https://doi.org/10.7554/eLife.11888 (2016).
Article PubMed PubMed Central Google Scholar
Andersson, A. F., Riemann, L. & Bertilsson, S. Pyrosequencing reveals contrasting seasonal dynamics of taxa within Baltic Sea bacterioplankton communities. ISME J. 4, 171–181. https://doi.org/10.1038/ismej.2009.108 (2010).
Article PubMed Google Scholar
Xu, H. et al. Distinct successional patterns and processes of free-living and particle-attached bacterial communities throughout a phytoplankton bloom. Freshwater Biol. 65, 1363–1375. https://doi.org/10.1111/fwb.13505 (2020).
Article Google Scholar
Casini, M. et al. Trophic cascades promote threshold-like shifts in pelagic marine ecosystems. PNAS 106, 197–202 (2009).
Article ADS CAS PubMed Google Scholar
Tomczak, M. T. et al. Reference state, structure, regime shifts, and regulatory drivers in a coastal sea over the last century: The Central Baltic Sea case. Limnol. Oceanogr. https://doi.org/10.1002/lno.11975 (2022).
Article Google Scholar
Needham, D. M. & Fuhrman, J. A. Pronounced daily succession of phytoplankton, archaea and bacteria following a spring bloom. Nat. Microbiol. 1, 16005. https://doi.org/10.1038/nmicrobiol.2016.5 (2016).
Article CAS PubMed Google Scholar
Martin-Platero, A. M. et al. High resolution time series reveals cohesive but short-lived communities in coastal plankton. Nat. Commun. 9, 266. https://doi.org/10.1038/s41467-017-02571-4 (2018).
Article ADS CAS PubMed PubMed Central Google Scholar
Sörenson, E., Farnelid, H., Lindehoff, E. & Legrand, C. Resource partitioning between phytoplankton and bacteria in the coastal Baltic sea. Front. Mar. Sci. 7, 608244. https://doi.org/10.3389/fmars.2020.608244 (2020).
Article Google Scholar
Download references
During the years, many people have been involved in the LMO sampling. We would like to thank everyone that has contributed to the sampling, especially Arnautovic S. and Månsson A. Furthermore, we acknowledge the help and assistance from Northern Offshore Services (NOS), M/V Provider crew, E. ON and RWE on the samplings. We would also like to thank Baltar F, Berglöf K, Bertos-Fortis M, Carlsson E, Helin A, Israelsson S, Jessen E, Karlsson C, Lautin J, Lindh MV, Lundin D, Martínez-García S, Mattsson L, Mollica T, Muthusamy S, Nham Q, Nilsson E, Fast M, Osbeck C, Oskarsson T, Pérez-Martínez C, Sjöstedt J for help during sampling and laboratory work. The research was supported by the Swedish Research Council FORMAS Strong Research environment EcoChange (Ecosystem dynamics in the Baltic Sea in a changing climate) to CL and JP, and by the Linnaeus University Center for Ecology and Evolution in Microbial model Systems (EEMiS). We acknowledge the support from Science for Life Laboratory, the National Genomics Infrastructure, NGI, and Uppmax (compute project SNIC 2017/7-419 and storage project SNIC 2020/16- 76), Sweden, for providing assistance in massive parallel sequencing and computational infrastructure.
Open access funding provided by Linnaeus University.
Centre for Ecology and Evolution in Microbial Model Systems - EEMiS, Linnaeus University, 39182, Kalmar, Sweden
Emil Fridolfsson, Carina Bunse, Elin Lindehoff, Hanna Farnelid, Benjamin Pontiller, Kristofer Bergström, Jarone Pinhassi, Catherine Legrand & Samuel Hylander
Department of Marine Sciences, University of Gothenburg, 40530, Gothenburg, Sweden
Carina Bunse
GEOMAR Helmholtz Centre for Ocean Research Kiel, E24105, Kiel, Germany
Benjamin Pontiller
School of Business, Innovation and Sustainability, Halmstad University, 30118, Halmstad, Sweden
Catherine Legrand
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
C.L. and J.P. designed the original sampling schemes. E.F., C.B., E.L., H.F., J.P., C.L., S.H. contributed to developing the sampling plan to its current state. E.F., C.B., B.P., and K.B. took part in the field and laboratory samplings. E.F. and C.B. also took part in the analyses and compiled all the variables. E.F., C.B., E.L., H.F., B.P., K.B., J.P., C.L., S.H. took part in the interpretation of the data. E.F. wrote the first draft and all authors contributed to finalizing the manuscript and approved the submitted version.
Correspondence to Jarone Pinhassi, Catherine Legrand or Samuel Hylander.
The authors declare no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
Reprints and Permissions
Fridolfsson, E., Bunse, C., Lindehoff, E. et al. Multiyear analysis uncovers coordinated seasonality in stocks and composition of the planktonic food web in the Baltic Sea proper. Sci Rep 13, 11865 (2023). https://doi.org/10.1038/s41598-023-38816-0
Download citation
Received: 12 October 2022
Accepted: 15 July 2023
Published: 22 July 2023
DOI: https://doi.org/10.1038/s41598-023-38816-0
Anyone you share the following link with will be able to read this content:
Sorry, a shareable link is not currently available for this article.
Provided by the Springer Nature SharedIt content-sharing initiative
By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.