Research Article |
Corresponding author: Joanna Grabowska ( joanna.grabowska@biol.uni.lodz.pl ) Academic editor: Emili García-Berthou
© 2024 Joanna Grabowska, Mateusz Płóciennik, Michał Grabowski.
This is an open access article distributed under the terms of the Creative Commons Attribution License (CC BY 4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Citation:
Grabowska J, Płóciennik M, Grabowski M (2024) Detailed analysis of prey taxonomic composition indicates feeding habitat partitioning amongst co-occurring invasive gobies and native European perch. NeoBiota 92: 1-23. https://doi.org/10.3897/neobiota.92.116033
|
One of the negative impacts of non-native invasive species on trophic interactions in an invaded ecosystem occurs via increased interspecific competition for food resources between the invader and local species of the same food niche. In freshwaters, there are usually several fish species that feed on similar food resources. Ponto-Caspian gobies are amongst the most successful and widespread invaders colonising European waterways. They have a wide food niche and an opportunistic feeding strategy, with a focus on benthic invertebrates and piscivory occurring occasionally mainly in the case of large individuals. Competition with native percids for food resources is predicted on the basis of high dietary overlap. However, studies published so far provide no unequivocal answer. In order to resolve this question, we conducted a comparative taxonomic analysis of gut content, with an emphasis on chironomids and amphipods, of the invasive monkey goby (Neogobius fluviatilis), racer goby (Babka gymnotrachelus) and the native Eurasian perch (Perca fluviatilis) occurring sympatrically in a large lowland European river, the Bug River in Poland. We found that each species forages in slightly different habitats, as indicated by the different composition of prey species in the gut content. This suggests feeding niche partitioning between the studied species facilitating their co-existence and reduction or avoidance of competition for food resources. Resource partitioning regarding prey types and foraging habitats is a mechanism for permitting the co-existence of closely-related alien gobies with similar food preferences in the invaded waters and co-occurrence with local species. This mechanism can contribute to their invasion success, as observed in European waters during the recent decades. We also demonstrate that precise prey identification to the lowest possible taxon is crucial to reveal the dietary overlap between co-occurring fish species and to predict the impact of alien invaders on native species through interspecific competition, as well as to recommend such an approach in studies upon fish foraging strategies.
Feeding niche, fish diet, non-native species, resource partitioning, trophic interactions
Freshwater ecosystems, together with their biodiversity, are amongst the most threatened and altered environments on the planet, due to the intensive human exploitation of water resources. Widespread invasions of introduced non-native species are amongst the five main threats for such ecosystems (
As fish species change their trophic status over the course of their lives, displaying ontogenetic niche shifts, many European freshwater fish do not fall into discrete trophic categories (
The racer goby Babka gymnotrachelus (Kessler, 1857) and monkey goby Neogobius fluviatilis (Pallas, 1814) are the first two invasive gobies that arrived in the Vistula River system, almost at the same time, in the mid-1990s and soon spread there rapidly (
The diet of racer and monkey gobies is similar and mainly comprises benthic macroinvertebrates, though chironomid larvae, other insects larvae, amphipods, molluscs and occasionally also small fish fry, predominate in their diet (
We hypothesised that the three fish species, although potentially feeding on the same type of prey, slightly vary the composition of their diet, for example, by exploring different foraging habitats to minimise interspecific competition when co-occurring in the same section of a river. We verified this prediction by comparative analysis of fish diet, based on the detailed taxonomic identification of selected prey taxa, focusing on chironomids and amphipods, which can differ in terms of the occupied microhabitats. We achieved this by analysing the gut content of racer goby, monkey goby and European perch occurring sympatrically in a large lowland river flowing through the East European Plain, which constitutes a crucial part of the Central Invasion Corridor for westward expansion of the Ponto-Caspian fauna, as defined by
All the three studied species, monkey goby, racer goby and perch, were sampled from three sites (Fig.
Location of sampling sites Z, R, B in the Western Bug River and places of first record of racer goby (red dot) and monkey goby (green dot) in Poland.
Morphometry of sampling sites: Z, R, B in the Western Bug River. Explanations: a) – absent, + very little/few, ++ common, +++ abundant; b) % of bed cover; c) % of bank overgrown; d) pa- pastures, ł- meadows, cr- cropland, bl- buildings.
Site | Z | R | B | |
---|---|---|---|---|
1. | Distance from the mouth [km] | 176.7 | 108.2 | 46.0 |
2. | Mean width [m] | 120 | 109 | 114 |
3. | Mean depth in current [m] | 1.8 | 0.8 | 1.5 |
Maximal depth in the current [m] | > 3.0 | 1.5 | > 3.0 | |
4.a) | Pools/riffles | +++ | + | + |
5. | Mud cover 0–100 [%] | 0 | 5 | 0 |
Bottom substrate: sand/gravel/stone [%] | 70/20/10 | 50/30/20 | 60/20/20 | |
6.b) | Submerged plants | 10 | 20 | 10 |
7.c) | Emerged plants | 30 | 50 | 30 |
8.a) | Trees along banks | +++ | +++ | + |
9.d) | Adjacent area | pa, ł | cr, pa, bl | bl, ł, pa |
10. | Water temperature [oC] | 19.0 | 19.1 | 19.2 |
11. | pH | 8.62 | 8.75 | 8.89 |
12. | Dissolved oxygen [mg O2 dm-3] | 13.4 | 5.9 | 9.26 |
13 | Conductivity [μS cm-1] | 509 | 502 | 476 |
Fish were sampled at depths from 0.5 m to 1.7 m, along the riverbank by electrofishing with a battery-powered unit, 350 V, 20–100 Hz, wading ca.100 m upstream along the bank and from the boat drifting 500 m downstream.
Fish were anaesthetised (MS-222) and preserved in 4% formaldehyde. In the laboratory, the fish were measured for total length (LT; to the nearest 1 mm) and weighed (with 0.01 g accuracy). Their alimentary tracts were dissected. The gut contents (in each fish, the same section of alimentary tract, i.e. stomach and first half of intestine) were weighed (to 0.0001 g accuracy) and prey items were identified under a stereomicroscope. Chironomidae larval stages can be identified to genera or groups of closely-related species, but only rarely to the species level. Their remains from the fish gut lack many features that are necessary for precise identification and, thus, following the main key used in this study (
Animal prey remains were identified to the lowest readily recognisable taxon, counted and the proportional weight was estimated. The percentage contribution by weight of each food category to the biomass of total stomach content was estimated visually (
The
Dietary overlap between each pair of fish species was calculated using Schoener’s index (Wallace 1981): α = 1–0.5 [Ʃn = 1 (pij – pik)], where pij is the proportion of the ith resource used by species j and pik is the proportion of the ith resources used by species k; overlap values exceeding 0.6 were regarded as high or biologically significant (Wallace 1981). As the proportion in Schoener’s index calculation, we used %N proportion: numbers of given prey type to the total number of prey found in fish gut. To show how the accuracy of prey identification influences the evaluation of dietary overlap, we estimated it based on protocol 1 – considering main food categories, i.e. prey pooled into taxonomic groups usually applied in fish diet studies, for example, Amphipoda, Chironomidae (called later Schoener’s index 1) and protocol 2 – considering detailed food categories, i.e. prey identified to the lowest possible taxonomic level (called later Schoener’s index 2).
To compare the taxonomic composition of the diet between fish species overall (all sites pooled) and at each sampling site, one-way permutation analysis of similarity (ANOSIM, Bray-Curtis similarity coefficient) was used, based on prey. ANOSIM is analogous to an ANOVA procedure, with non-parametric permutation applied to rank similarity matrix of samples. The similarity percentage procedure (SIMPER) was used to identify which prey taxa were most likely responsible for the patterns detected by ANOSIM. SIMPER provided the average dissimilarities between the species and identified which prey taxa made the greatest contribution to any dissimilarities between analysed categories (
Dietary niche width was calculated as a Simpson diversity index: 1 – D = 1 – Σpi 2 and Shannon diversity index: H = – Ʃpi log2 pi, where pi is the proportion of different prey in the diet (
The fish species recorded from the sampling sites were mainly bleak Alburnus alburnus and roach Rutilus rutilus. These two species constituted 45%–64% of all fish caught at the study sites and dominated in abundance along the whole middle and lower river course. The other species that occurred at all three sites were common bream Abramis brama, white bream Blicca bjoerkna, pike Esox lucius, chub Squalius cephalus, ide Leuciscus idus, common rudd Scardinus erythrophtalmus, spined loach Cobitis taenia and bitterling Rhodeus sericeus (Suppl. material
In total, 63 individuals of racer goby, 77 of monkey goby and 62 of perch were caught in three sampling sites. In four out of 202 dissected individuals, the alimentary tracts were empty and not considered in further analysis.
In all three fish species, prey belonging to Amphipoda, Chironomidae larvae and pupae, Gastropoda, Trichoptera larvae, Coleoptera larvae, Oligochaeta and Hirudinea were found in the diet (Suppl. material
Feeding strategy displayed using the
Diet composition of the European perch, racer goby and monkey goby (mean, minimum and maximum total length of fish – TL) expressed as relative abundance (%N) and relative biomass (%B) of main food categories in gut content at the three studied sites (Z, R, B).
Species | Racer goby | Monkey goby | European perch | |||||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
Site | Z | R | B | Z | R | B | Z | R | ||||||||
N of specimens | 34 | 20 | 27 | 20 | 32 | 24 | 31 | 30 | ||||||||
Mean TL [mm] | 70.24 (±14.50) | 60.70 (±15.53) | 74.65 (±14.50) | 96.90 (±27.36) | 92.36 (±14.68) | 85.86 (±16.03) | 95.24 (±29.95) | 98.154 (±7.95) | ||||||||
Min-max TL [mm] | 51–101 | 42–96 | 52–103 | 52–124 | 61–120 | 54–112 | 53–150 | 84–115 | ||||||||
Food categories | %N | %B | %N | %B | %N | %B | %N | %B | %N | %B | %N | %B | %N | %B | %N | %B |
Amphipoda | 27.5 | 48 | 19.6 | 65 | 13.8 | 38 | 16.7 | 22 | 6.1 | 19 | 54.9 | 53 | 54.2 | 57 | ||
Chironomidae larvae | 43.7 | 12 | 67.4 | 29 | 78.8 | 17 | 39.7 | 8 | 76.2 | 44 | 76.1 | 17 | 20.4 | 6 | 23.8 | 5 |
Chironomidae pupas | 4.8 | 4 | 4.3 | 1 | 3.8 | 3 | 5.4 | 7 | 1.0 | 2 | 2.6 | 2 | 7.0 | 3 | ||
Trichoptera larvae | 0.6 | > 1 | 12.8 | 4 | 2.2 | 3 | 18.1 | 18 | 0.4 | > 1 | 0.5 | 2 | ||||
Odonata larvae | 0.6 | > 1 | 0.5 | 6 | 1.0 | 1.7 | > 1 | |||||||||
Coleoptera larvae | 1.8 | 2 | 2.5 | 1 | 1.0 | 4 | 0.9 | > 1 | 1.9 | 2 | ||||||
Bivalvia | 14.4 | 19 | 2.2 | 1 | 11.5 | 26 | 7.9 | 5 | 1.9 | 6 | ||||||
Gastropoda | 5.4 | 4 | 6.5 | 4 | 2.5 | 10 | 12.8 | 27 | 1.8 | 2 | 1.4 | 2 | ||||
Mollusca not ident. | 6 | 14 | ||||||||||||||
Oligochaeta | 1.3 | 1 | 1.0 | 6 | 0.9 | 6 | ||||||||||
Hirudinea | 0.6 | > 1 | 2.6 | 10 | ||||||||||||
Pisces | 0.6 | > 1 | 1.3 | 21 | 19.1 | 34 | 10.3 | 23 | ||||||||
Detritus (plant) | 12 | 23 | 2 | |||||||||||||
Sand | 12 | |||||||||||||||
Fish eggs | 12 | 12 |
The plot of prey specific abundance (%Nps) and frequency of occurrence (%F) of the main components of the diet showed that chironomid larvae were the prey of higher importance for gobies, while, for European perch, amphipods were more important (Fig.
Diet overlap, as calculated for the 11 main food categories (Schoener’s index 1), occurred amongst all three species if data for all sites were pooled, which indicated that their prey spectrum was very similar (Table
The dietary overlap estimated, based on two protocols: calculated for general (Schoener’s index 1) and detailed food identification (Schoener’s index 2) categories. Pairwise comparisons (Bonferroni test) of fish diet following one-way ANOSIM and SIMPER analysis based on detailed identified food categories.
Comparisons | Schoener’s index 1 | Schoener’s index 2 | ANOSIM | SIMPER | |
---|---|---|---|---|---|
R | p | dissimilarity | |||
Fish species (site pooled) | |||||
perch vs racer goby vs. monkey goby | 0.1893 | 0.0001 | 91.08 | ||
perch vs. racer goby | 0.826 | 0.424 | 0.1333 | 0.0060 | 87.82 |
perch vs. monkey goby | 0.721 | 0.220 | 0.3520 | 0.0030 | 94.95 |
racer goby vs. monkey goby | 0.818 | 0.480 | 0.1266 | 0.0030 | 91.22 |
Site Z | |||||
perch vs. racer goby vs. monkey goby | 0.2953 | 0.0001 | 91.07 | ||
perch vs. racer goby | 0.530 | 0.346 | 0.3543 | 0.0003 | 88.47 |
perch vs. monkey goby | 0.400 | 0.175 | 0.0353 | 0.9140 | 97.26 |
racer goby vs. monkey goby | 0.784 | 0.572 | 0.2656 | 0.0470 | 92.53 |
Site R | |||||
perch vs. racer goby vs. monkey goby | 0.1883 | 0.0010 | 91.47 | ||
perch vs. racer goby | 0.491 | 0.274 | 0.1393 | 0.0003 | 84.30 |
perch vs. monkey goby | 0.370 | 0.175 | 0.3250 | 1.0000 | 94.57 |
racer goby vs. monkey goby | 0.816 | 0.556 | -0.0870 | 1.0000 | 90.53 |
Site B | |||||
racer goby vs. monkey goby | 0.790 | 0.389 | 0.3243 | 0.0130 | 88.83 |
Up to 56 taxa (including 28 chironomids and 4 amphipods) dominated the food categories shared by the three studied fish species. For the analysis, we rejected taxa that were found in only one fish, which reduced the number of prey taxa to 42 (including 24 chironomids). The values of Schoener’s index 2 indicated that there was no dietary overlap between gobies and perch at any site, but there was also no dietary overlap between racer goby and monkey goby at site B or it was moderate (ca. 0.5) at the other two sites, Z and R (Table
Results of SIMPER analysis identifying prey categories with the highest contribution to the overall dissimilarity amongst fish species diets and their mean relative abundance (%N) in diets of perch (PF), racer goby (BG) and monkey goby (NF).
Taxon | Contribution % | Cumulative % | Mean PF | Mean BG | Mean NF |
---|---|---|---|---|---|
Glyptotendipes cf. pallens | 15.1 | 15.1 | 17.0 | 28.0 | 11.9 |
Dikerogammarus villosus | 14.1 | 29.2 | 29.1 | 13.3 | 0.1 |
Chironomus cf. riparius | 9.6 | 38.8 | 0.0 | 3.8 | 22.8 |
Pisces | 7.8 | 46.6 | 21.0 | 1.3 | 0.0 |
Pontogammarus robustoides | 6.4 | 53.0 | 14.6 | 2.3 | 1.3 |
Polypedilum cf. nubeculosum | 4.8 | 57.8 | 2.8 | 5.1 | 6.1 |
Sphaeriidae | 4.4 | 62.3 | 0.0 | 7.5 | 3.8 |
Microtendipes cf. pedellus | 4.3 | 66.5 | 0.0 | 4.8 | 6.8 |
Gastropoda not identified | 4.2 | 70.7 | 0.6 | 2.4 | 8.0 |
Hydropsyche sp. larvae | 3.6 | 74.3 | 0.5 | 0.0 | 8.9 |
Dikerogammarus not identified | 2.4 | 76.7 | 0.3 | 2.7 | 3.4 |
Rheocricotopus cf. chalybeatus | 2.1 | 78.8 | 0.0 | 3.6 | 1.6 |
Gammaridae not identified | 1.9 | 80.7 | 1.3 | 2.4 | 1.5 |
Dicrotendipes cf. nervosus | 1.9 | 82.6 | 0.6 | 3.1 | 1.2 |
Glyptotendipes cauliginellus pupae | 1.6 | 84.2 | 2.4 | 1.3 | 0.9 |
Chaetogammarus ischnus | 1.5 | 85.7 | 1.8 | 2.2 | 0.0 |
Cryptochironomus sp. | 1.4 | 87.1 | 0.0 | 0.0 | 3.6 |
Dikerogammarus haemobaphes | 1.3 | 88.4 | 0.4 | 1.5 | 1.4 |
Rheotanytarsus sp. | 1.2 | 89.6 | 0.0 | 1.3 | 1.8 |
Coleoptera (Gyrinus sp.) | 1.1 | 90.7 | 1.7 | 1.2 | 0.0 |
Lipiniella moderata | 1.1 | 91.8 | 1.2 | 0.3 | 1.5 |
Tanytarsini not identified | 1.1 | 92.9 | 0.0 | 0.3 | 2.4 |
Zygoptera larvae | 1.0 | 93.8 | 2.4 | 0.3 | 0.0 |
Trichoptera larvae not identified | 0.9 | 94.7 | 0.0 | 0.0 | 2.3 |
Coleoptera larvae not identified | 0.7 | 95.5 | 0.5 | 0.9 | 0.4 |
Polypedilum cf. sordens | 0.5 | 95.9 | 0.0 | 1.2 | 0.0 |
Ten taxa, i.e. Glyptotendipes cf. pallens, Dikerogammarus villosus, Chironomus cf. riparius, Pisces, Pontogammarus robustoides, Polypedilum cf. nubeculosum, Microtendipes pedellus-type, gastropods, caseless larvae (Hydropsyche sp.) of Trichoptera and Rheocricotopus cf. chalybeatus, out of 56 analysed, contributed to 80% of the overall dissimilarity amongst the diets of perch, racer and monkey gobies, though the mean abundance of particular prey varied between sites (Fig.
Average relative abundance of prey taxa (%N) in the gut contents of PF – European perch, BG – racer goby, NF – monkey goby, which, in total, contributed to 95% of dissimilarity (SIMPER) amongst fish species diets at A site Z B site R and C site B. n.d. – means not identified.
The perch mainly fed on amphipods, i.e. P. robustoides at site Z, D. villosus at site R and fish at sites Z and R. Amongst the Chironomidae larvae, the G. cf. pallens contributed the most to the perch diet at each site. This chironomid was the most abundant in the diet of racer goby at sites Z and B. D. villosus was also an important food item of racer goby at site Z. Sphaeridae were not recorded in the diet of perch, but contributed to the diet of both goby species, especially at site Z, where monkey goby fed also on gastropods. Contrary to the other two co-occurring fish species, monkey goby consumed many caseless trichopteran larvae, as well as the chironomids: C. cf. riparius at all sites and P. cf. nubeculosum at site B.
The prey diversity was lower for perch than for gobies. Concerning the latter, prey diversity tended to be higher for racer goby than for monkey goby at sites Z and R; however, it was equal at site B (Fig.
The studied fish species, native European perch as well as non-native racer goby and monkey goby, fed on similar prey taxa, which suggests a high dietary overlap. Nevertheless, more detailed identification of taxa in the most abundant food categories, i.e. chironomid larvae and amphipods, revealed that they foraged on different prey at sites where they co-occurred. Thus, although the majority of prey taxa were recorded in guts of all the three studied fish species, their contribution to the diet at a given site was different. This supports the hypothesis of resource partitioning to avoid competition for food between native and non-native species.
Although several experimental studies showed the higher competitive ability of invader versus native species and the greater potential of the former to utilise resources (
Moreover, our findings proved that accuracy in taxonomic identification of prey taxa is essential to provide reliable data for dietary overlap or resource partitioning assessment. It is especially crucial in the case of fish species, for example, racer goby and monkey goby, feeding on the same type of prey that is very diverse considering its body size and occupied microhabitats. Identification of prey to the lowest possible taxon also allows us to determine the habitat preferences of fish species based on the knowledge of their prey microhabitat preferences. Our results showed that, in the case of gobies, resource partitioning is realised by utilisation of different habitats.
Both goby species and European perch fed on the macrozoobenthos. The perch is known to shift toward piscivory with its ontogenetic development (
In general, the diets were more similar between the goby species than between either of the gobies and perch. However, the diet of perch was more similar to that of the racer goby than to that of the monkey goby. Both the Eurasian perch and the racer goby fed on prey that indicated their association with macrophytes. Macrophyte patches are refuges for small fish, as well as hiding places or substrate for several macrozoobenthic groups, such as amphipods, insect larvae, for example, Diptera and Zygoptera larvae or gastropods, that are attractive food for many fish species (
Amphipods were especially important food items for perch and racer goby in the Bug River. Depending on the site, these fish mainly ate Pontogammarus robustoides or Dikerogammarus villosus and less D. haemobaphes. Field observations have shown that all three species are rather eurytopic (
Both goby species and European perch feed on small gastropods, such as Bithynia sp., Valvata sp. and Potammopyrgus antipodarum, which are also associated with submerged macrophytes (
In summary, alien gobies, in particular the racer goby, and European perch possibly used similar habitats for foraging, i.e. macrophyte patches in areas of more stagnant water and muddy bottom. However, perch with a body length similar to that of co-occurring gobies, was more piscivorous. The dietary overlap between perch and gobies usually comprised prey items that are very common in the riverine environment, such as amphipods and large chironomid larvae (
Racer goby and monkey goby had similar diets. They fed mainly on Chironomidae larvae, on the basis of the relative abundance and frequency of this prey in the fish gut content. The detailed identification of taxa within this food category showed that, in fact, the gobies foraged in different microhabitats, even at the same sites and their mode of foraging was also slightly different. Our study shows that several taxa of chironomids contributed to 60% dissimilarity between the diet of the studied goby species. Chironomids are a prevalent group in the freshwater macrozoobenthos, often standing out in their abundance and species and functional diversity, which makes them key elements of freshwater food webs (
Glyptotendipes cf. pallens dominated amongst chironomid larvae in the diet of racer goby. In summer, this taxon can be found in silty tubes built on macrophytes, mining in their decaying parts, but also on other firm surfaces, such as decaying wood or stones. In large rivers, in particular, the larvae of this species are more numerous on stones than on plants. This taxon avoids fast running waters and prefers more stagnant parts of the river channel (
To conclude, we show that detailed prey identification to the lowest possible taxon is crucial to properly justify the diet overlap between co-occurring fish species and to verify the suggested impact of alien invaders on native species through interspecific competition. Resource partitioning considering prey types and foraging habitats is one of the ways of allowing the co-existence of closely-related alien gobies with similar food preferences in the invaded waters and their co-occurrence with local fish species. Together with an opportunistic feeding strategy, it is likely to be a major factor behind their invasion success observed in European waters in the last decades. We therefore recommend that, in order to gain more detailed insights into the foraging strategy of fish, in future studies, researchers should not limit their dietary analysis only to the identification of higher taxa, but should identify prey down to the lowest possible level, especially in taxonomic groups consisting of species that differ in the microhabitats they occupy.
We would like to thank Małgorzata Dukowska and Joanna Leszczyńska, University of Lodz, for their comments on the species-specific habitats of chironomid larvae. We are also grateful to Stephen Venn, University of Lodz, for language editing and all the comments that helped us to improve the quality of our manuscript. Finally, we would like to acknowledge all the colleagues that helped with sampling: Dariusz Pietraszewski, Szymon Tybulczuk, Lidia Marszał, Grzegorz Zięba, Łukasz Kapusta and Andrzej Kruk.
The authors have declared that no competing interests exist.
No ethical statement was reported.
No funding was reported.
Conceptualization: JG. Formal analysis: JG, MP. Investigation: MG, MP, JG. Resources: JG. Visualization: JG. Writing – original draft: MG, JG. Writing – review and editing: MP.
Joanna Grabowska https://orcid.org/0000-0001-9924-0650
Mateusz Płóciennik https://orcid.org/0000-0003-1487-6698
Michał Grabowski https://orcid.org/0000-0002-4551-3454
All of the data that support the findings of this study are available in the main text or Supplementary Information.
Relative abundance of species (%N) in fish assemblages found at sites Z, R, B in the Western Bug River in August 2007 (
Data type: xlsx
Relative abundance of prey categories (%N) (number of given prey category in relation to total number of prey) identified in fish guts at sites Z, R, B in the Western Bug River in August 2007
Data type: xlsx