4.1. August Planktonic Food Resources of Juvenile Sticklebacks
In general, the zooplankton of the studied sites was represented by a boreal warm-watermass related taxa assemblage typical of the coastal areas of the White Sea in summer [
29]. Copepoda was the dominant group, reaching 90-100% of total biomass and density of the whole zooplankton community; with the addition of taxa belonging to Cladocera. These results correspond well to major studies of zooplankton in the region that describe Copepoda as the major contributing group to the White Sea zooplankton [
30,
31].
Absolute zooplankton biomass varied between 38.5 and 5040 mg/m
3 with the general mean of 779 mg/m
3. When Koliushkovaya lagoon is excluded from the analysis, due to exclusive dominance of
A. longiremis, total biomass of zooplankton in ‘open’ marine sites was 281 mg/m
3. These values correspond well to the within-decadal variation of zooplankton biomass in Kandalaksha Bay between 140 and 257 mg/m
3 [
32]. These studies also indicated that the zooplankton biomass is generally higher in coastal areas compared to the open sea [
29].
Among the nursery grounds studied, the highest total biomass of zooplankton was observed at eelgrass (
Zostera marina) beds. As dense eelgrass thickets significantly reduce the water currents, this community located on the edge of subtidal and intertidal zones form a unique habitat for both planktonic and benthic species [
33,
34]. Locally calm hydrodynamic conditions provide optimal habitat for zooplankton [
35,
36] and thus support the extensive food resources for juvenile sticklebacks for which eelgrass beds are also a most important nursery habitat.
4.4. Stomach Content in Juvenile G. aculeatus
In total, 31 prey items were recorded in stomachs of juvenile sticklebacks, varying between seven and 17 at various nursery habitats. According to the studies conducted in 1960s [
42] the number of taxa in stomach contents varied between five and 25, indicating the same level of diversity of prey objects in stickleback diet throughout a five-decade period between studies. In general, juvenile sticklebacks were found to be planktivorous, as planktonic taxa prevailed in stomach content by abundance and biomass. Among the diversity of prey items, nine benthic taxa were found in stomachs, rarely contributing up to 46% by biomass to the stomach content. Krugliash was the only site where benthic taxa were important for the stickleback diet during both early (36%) and late (46%) August.
Prevalence of planktonic diet in juvenile stickleback diet has previously been reported from the White Sea [
21,
37], the Baltic Sea [
43] and freshwater lakes at British Isles [
44]. In adult White Sea three-spined sticklebacks, diversity of prey items is also comparable [
18], and similar values of diversity of prey items have been observed in nine-spined stickleback
P. pungitius from lacustrine environments despite the fact that adult fishes were mostly considered in the study [
45].
D-index and number of taxa in individual stomachs varied significantly between dates and sites and did not correlate to FI. Fish diet varied between early and late August and no relationship between the size of fishes (TL) and taxonomic richness of prey items in stomachs was observed. No clear differentiation between the taxonomic richness of prey items consumed by stickleback juveniles related to the particular type of marine nursery grounds was found. Nevertheless, within each sampling period, the diversity of prey items observed in stomachs was significantly different between most sites with minor exceptions indicating differences in the availability of prey items.
The highest number of effective taxa was observed in Letnaya and Podpakhta (4.6 ± 0.16 and 4.5 ± 0.16), contrasting with marine lagoon Koliushkovaya where the lowest number of effective taxa (2.52 ± 0.14) was found in fish diet. At the eelgrass site Seldianaya the diversity of prey items was high in early August, and on average in late August. In the stomachs of individual fishes, 2-5 prey items were observed, indicating that despite the broad taxonomic range of potential prey items found in stomachs at each site, these <5 prey items prevail in fish diet. According to our data, 70% of fish diet was formed by ‘effective’, and other 30% were randomly consumed. Across all the sites studied, two taxa were the most important contributors to the diet of juvenile sticklebacks, namely
T. longicornis and Orthocladiinae, similar to the data from Seldianaya obtained earlier [
21], with the only exception of the marine lagoon, where
A.longiremis overwhelmingly prevailed in the zooplankton. In late August samples, the proportion of
T. longicornis was higher, probably indicating that with the growth of individual fishes, their ability to catch this crustacean increases. Consumption of chironomid larvae (Orthocladiinae) by stickleback larvae probably indicates the patchy distribution of this prey item. Chironomids are common and an important element of fish diet in marine and freshwater ecosystems [
4,
44,
45,
46]. In White Sea juvenile sticklebacks, when Orthocladiinae are present in stomach contents, they represent the majority of the biomass of food consumed, while
T. longicornis occurs evenly across the individual fishes within the sample. On average, three individuals of Orthicladiinae were observed per stomach of
G. aculeatus, with a maximum of eight larvae. No differences in consumption rate of Orthocladiinae were observed between sites and dates of sample collection. Comparison of consumption of other effective taxa by stickleback juveniles revealed that consumption of
P. leuckartii varied significantly among sites only, while consumption of
T. longicornis,
M. norvegica, Copepoda juv. and
A. longiremis varied between both sampling dates and sites. Observed spatial and temporal variability in diet of juvenile sticklebacks correspond well to the previous findings [
26,
42,
43,
47].
4.5. Prey Selectivity by Juvenile Sticklebacks in the White Sea
Prey selectivity has previously been studied in both adult and juvenile three-spined sticklebacks, yet these studies were to a great extent implemented using aquarium experiments or in freshwater environments and more often were focused on adult fishes. Yet, these studies contributed significantly to an understanding of the role of local hydrology in trophic preferences [
48], profitability and visual characteristics of prey items [
49], and competition with other species within the same habitat [
50,
51], with a special focus on ecosystems where three-spined stickleback is a non-native species [
52,
53]. As the body size in sticklebacks is comparable to many macroinvertebrate predators (e.g., insect larvae), potential overlaps of their dietary preferences were also considered [
54]. One of the most detailed descriptions of dietary shifts of sympatric species of sticklebacks, including
G. aculeatus,
P. pungitius,
G. wheatlandi and
Apeltes quadractus was implemented at tidal saltmarsh spawning grounds in Canada [
55]. This study considered seasonal dynamics in selective feeding of sticklebacks from early summer when adults come to the spawning grounds in late August when both adults and juveniles were present in samples [
55].
In our study, in early August, T. longicornis was selectively preferred by juvenile sticklebacks (E values up to 1) at all sites studied, based on E index averaged per sample. Considering individually calculated E index values, similar results or neutral selectivity were obtained. Selectivity for T. longicornis (mean body size 0.56 ± 0.003 mm) varied among sites in late August when both approaches of calculating E index were applied. While E index values for average fish varied from - 0.2 to + 0.9 at various sites, individual data yielded values between - 0.42 to + 0.34, both indicating neutral or even avoiding a pattern of selectivity for this species. Similar seasonal variation in selectivity by sticklebacks was observed for P. leuckartii. At most sites, this species was selected positively during early August. In late August, positive selectivity for this species was observed at Keret (reaching + 0.64 by individuals and + 0.8 by average fish). On the contrary, P. leuckartii was avoided by juveniles at Seldianaya in late August, according to both methods of calculating E index.
Two prey taxa (
M. norvegica and Copepoda juv.), being the smallest organisms in stickleback diet considered, demonstrated opposed patterns of selectivity over the two periods of observation. Mean size of
M. norvegica was 0.41 ± 0.001 mm), and Copepoda juv. were mainly represented by nauplii Copepoda (0.11 ± 0.001 mm). These taxa were among the prey items of mainly smaller juvenile sticklebacks. Both prey items were negatively or neutrally selected in early August (
E from -0.75 to +0.06 for
M. norvegica and from - 0.8 to - 0.16 for Copepoda juv.), and positively or neutrally selected in late August (
E from + 0.23 to + 1 for
M. norvegica and from - 0.19 to + 1 for Copepoda juv.). A plausible explanation for the change in dietary preferences of these fishes might be that this shift is linked to changes in abundance of prey items. In early August, both
M. norvegica and Copepoda juv. each contributed about 15-20% to the total zooplankton biomass (
Figure 4). Later, their abundance in the zooplankton decreased, but small sizes of juvenile sticklebacks were still present at the nursery grounds due to the extended spawning period (
Figure 2). The observed decrease in abundance of
M. norvegica and Copepoda juv. in the zooplankton might be related either to seasonal trends or be a consequence of selective feeding by juvenile sticklebacks.
Several dominant planktonic taxa were mostly avoided by sticklebacks. In both periods of observation,
O. similis was avoided at all sites where this species was present in the zooplankton (see
Figure 3 and
Figure 4 for details). As the size of individual
O. similis (0.4 ± 0.005 mm) is similar to
M. norvegica, differential selectivity patterns for these two species can be linked to several factors. First,
M. norvegica (as in many other Harpacticoida) is bright-coloured (red to orange). As the same part of the spectrum is recognizable by at least female sticklebacks during the mating season [
56], this can potentially increase visibility of this species for juvenile sticklebacks. Second, these species demonstrate very different swimming behaviour. Despite spending most of the time motionless, cyclopoids belonging to genus
Oithona (e.g.
O. plumifera) are able to jump fast for up to 20 body lengths (speed 22.1 mm/s), and almost exclusively upwards [
57]. This behaviour can be beneficial for
O. similis to avoid attacks by sticklebacks. On the contrary,
M. norvegica demonstrate very slow swimming velocity (0.6 mm/s) [
58], and even
T. longicornis,almost always selected by sticklebacks, has a speed of 5.9 mm/s in adults, and lower values during early life stages [
59].
Another species selectively avoided (
E index values mostly around -1) by juvenile sticklebacks was
A. longiremis at nearly all sites and over both periods. In zooplankton communities this species reached about 60% of total biomass, and in the marine lagoon Koliushkovaya reached 100% of total biomass. Considering that FI of stickleback juveniles at this site was generally low (
Figure 3), it is possible to assume that this is a case of forced feeding on the only available prey item despite high prey availability. The possible evidence for forced dietary preferences in the lagoon is that selectivity for the second prey item group available - Copepoda juv. – was highly positive despite remarkably lower abundance in the zooplankton. This probably indicates that when individual fish could find larvae of copepods (remarkably this group can also include early stages of
A. longiremis to some extent), it preferred to feed on them, and if not, the only available food were adult
A. longiremis.
While comparison of two methods of calculating
E index mostly yielded similar results, three cases can be considered: (1)
E index values calculated per individual and per average fish both simultaneously indicated either very high or very low positive selectivity, (2) each index indicated values within the ‘neutral’ selectivity zone between - 0.25 and + 0.25, (3) averaged fish index indicated positive selectivity while individual-based index indicated slightly negative or neutral selectivity. Individual variation of prey selectivity for
P. leuckartii by
G. aculeatus at various sites in early August is a good example of these three patterns (
Figure 6). Thus, the conventional way of calculating
E index, i.e., the values calculated per average fish at the nursery ground (i.e., local population), mostly indicated a general trend of prey selectivity by individual stickleback population at each nursery ground. On the contrary, individual
E values provided more accurate estimates and indicated sites and dates where various fishes within the sample had different consumption of prey item by individuals within the local population.