Research Article |
Corresponding author: Piotr Kłosiński ( klosinski.piotr@doktorant.umk.pl ) Academic editor: Nicola Smith
© 2025 Piotr Kłosiński, Jarosław Kobak, Tomasz Kakareko.
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:
Kłosiński P, Kobak J, Kakareko T (2025) Competitive interactions for food resources between invasive Ponto-Caspian gobies and their native competitors in the context of global warming. NeoBiota 97: 91-119. https://doi.org/10.3897/neobiota.97.134566
|
Climate warming can modify the process of biological invasions by affecting the outcomes of competition between alien species and their native counterparts in invaded environments. Inland freshwaters are particularly vulnerable to the intensification of such phenomena due to the accumulation of invaders, including thermophilic species that may benefit from warming. We intended to check whether an elevated summer temperature (25 vs. 17 °C) affects the abilities of the Ponto-Caspian gobies to compete for food. These fish are considered temperature-tolerant, highly invasive freshwater fish in Europe. In laboratory experiments, we tested single- and two-species pairs of juvenile specimens of two goby species and their native counterparts from the same ecological guild (the racer goby Babka gymnotrachelus versus European bullhead Cottus gobio, and monkey goby Neogobius fluviatilis versus native gudgeon Gobio gobio). The fish competed for food (live chironomidae larvae provided at rates below satiation) for 1 hour at night. We analysed behaviours associated with direct interactions (aggression acts) and foraging activity (time to enter the feeder and the time spent in the feeder). We found that although the gobies did not show higher aggression than the natives, they more actively accessed food compared to the latter, irrespective of temperature. Our results suggest that, in the wild, the invasive fish have a competitive advantage over the native ones due to better resource allocation (gaining food without incurring the costs of aggression) and will maintain this advantage as water warming continues.
Aggressive behaviour, biological invasions, climate change, food competition, freshwater species, nonnative species
Nowadays, biological invasions constitute a leading threat to global biodiversity (
Invasions of alien species can be aided by changes in environmental conditions, such as climate warming (
The Ponto-Caspian region constitutes the major donor of alien taxa for European waters (
The success of the Ponto-Caspian gobies is often linked to their effective competition (
Ongoing global warming can reconfigure interspecific interactions between invasive species and their native counterparts (
The tolerance to elevated temperature and competitive efficiency are separate issues which likely interact with each other, but their interacting effects are unknown. Therefore, we aimed to study these interactions experimentally. First, we assumed that interspecific competition between alien gobies and their native counterparts belonging to the same ecological guild is an effect of overlapping food niches (
We tested two goby species of Ponto-Caspian origin, the racer goby and monkey goby, paired with their coexisting native competitors: the European bullhead and gudgeon Gobio gobio (Linnaeus, 1758), respectively. These two pairs of species were chosen as they co-occur in the same habitats of European freshwater environments sharing similar biology and ecology (
Fish were transferred from the stock tanks to 85-L acclimation tanks in groups of 10–12 individuals, at an initial temperature of 17 °C (as in the stock tanks). The acclimation tanks were filled with conditioned (24 h aged, aerated) tap water and furnished in the same way as the experimental tanks (see below, “Experimental setup”). The progressive adjustment of a temperature up to 25 °C was reached within 8 days using aquarium heaters with an accuracy of 0.25 °C (AQUAEL Ultra Heater 150 W; Suwałki, Poland). During acclimation, fish were fed ad libitum once a day with unfrozen chironomid larvae and uneaten prey were removed from the acclimation tanks. Food was delivered with a small amount of water to the acclimation tank on the Petri dish (a feeder placed on the bottom) through the PVC hose and the transparent glass tube. The fish to be tested at 17 °C were transferred from the stock tanks to the acclimation tanks for the same amount of time, but not subject to other temperature alterations. After 8 days in the acclimation tanks, when the temperature reached 25 °C, the fish were transferred to the experimental tanks.
Experiments were carried out in 27-L tanks (30 × 30 × 30 cm) filled with aged (24 h), aerated tap water. To reduce the effects of handling and visual disturbance on the test fish, the experimental tanks were isolated on all sides by Styrofoam screens. Each tank was furnished with an aerator, two shelters, aquarium heater (between the shelters), and feeder (Suppl. materials
We took the fish for the research randomly, firstly from the field and then from the stock tanks. The total length of the fish was measured from digital photographs taken during tests using ImageJ 1.53k (freeware by W.S. Rasband, U.S. National Institutes of Health, Bethesda, Maryland, USA: https://imagej.net/ij/). Mean (± SD) total lengths (TL) were: 4.68 ± 0.58 cm, 4.87 ± 0.67 cm, 4.70 ± 0.69 cm and 5.31 ± 0.52 cm for the racer goby, bullhead, monkey goby and gudgeon, respectively. Within each species pair, the fish were tested in dyads of similar TL (average difference in TL of 0.10–0.26 cm). Mean TLs of fish in pairs were not significantly different between the species (Student’s t tests for dependent samples; see Suppl. material
The last feeding took place 40 h before the beginning of the experiment. Two fish (depending on the treatment) were selected from the acclimation tanks and placed in the experimental tank 16 h (at 15:00) before the start of the trial to get familiar with the experimental arena (the adaptation period) (Suppl. material
Analysis of all the video recordings of fish behaviour was carried out manually, always by the same person, to avoid bias due to differences in the interpretation of fish behaviour. We noted one variable related to aggression and two variables related to foraging: (1) the number of aggressive actions directed towards the opponent, when one fish moved quickly towards the other, which ended in a physical contact between the individuals, such as hitting or pushing (so, the opponents had to touch each other at some moment of the interaction to count the event as aggression). This allowed us to establish clear, strict and objective criteria of aggressiveness, which did not raise any doubts about their correct assessment by the observer; (2) the time to enter the feeder for the first time by each individual; (3) the percentage of time spent by the fish directly in the feeder, which was used as a proxy for food consumption, as it was challenging to observe it directly in darkness. We assume this as a good proxy for foraging, especially in the initial period of the exposure, directly after the food application, when the food was present in the feeder for sure. The animal needed to be present inside the feeder at this moment to have access to the food. In the one-species treatments, because of the visual similarity of the individuals, it was not possible to track them without mistaking particular individuals on video frames. Instead, the two individuals of the same species were tracked together and the final response consisted of summed up and averaged responses of these individuals.
We conducted the following types of statistical analyses: (i) comparison between the species within each pair in their single-species treatments (to test differences between the species); (ii) comparison between the species within each pair in the mixed-species treatment (to check which species has an advantage over the other when they are confronted in the same area); (iii) comparison of the behaviour of each species between the mixed vs. single species treatments (to test the impact of one species on the other). Dependent variables tested in the analyses were as follows: (i) the number of aggression events determined in six consecutive 10-min periods during the exposure (analysed using a General Linear Mixed Model; the use of a Generalized Linear Model designed for count data was not possible due to non-integer data points averaged for single species pairs); (ii) the time spent in the feeder (analysed using a General Linear Mixed Model); (iii) the time to enter the feeder (analysed using a Cox proportional hazard regression to account for the individuals that did not enter the feeder at all). Independent variables were as follows: (i) species (in the comparisons between the species); (ii) treatment (in the comparison between the mixed vs. single species treatments for each species); (iii) temperature (17 and 25 °C); (iv) exposure time counted since the food introduction to the feeder (for the models testing the number of aggression events and time spent in the feeder, a continuous covariate: 10, 20, 30, 40, 50, 60 min), (v) individual pair ID as a random factor (to group repeated measurements for each pair of individuals). Species was a within-subject factor when the species tested in mixed-species treatments were compared to each other. The summary of all the models used in the study is shown in Suppl. material
In all treatments, the number of aggression events exhibited by the racer goby and European bullhead decreased with time (a significant effect of exposure time), but was independent of temperature (Table
Numbers of aggression acts per 10 min (a single observation period) shown by the racer goby and European bullhead kept in separate single-species treatments (A) or together in the mixed-species treatment (B). Panels C, D present comparisons of the behaviour of the racer goby and European bullhead, respectively, between the single- and mixed-species treatments. Symbols represent raw data (means ± 95%CI) for each species, temperature and period. Lines are predicted by the models (with 95%CI as shaded areas). Common slopes were predicted for groups of data that did not differ significantly from each other in the models.
General Linear Mixed Models to test the impact of treatment, temperature, exposure time and species on the number of aggressive events shown by the racer goby and European bullhead. Non-significant higher order interactions were removed from the models in a simplification procedure.
Analysis | Effect | df | F | P |
---|---|---|---|---|
Racer goby vs European bullhead from single-species treatments | Species | 1, 38 | 8.43 | 0.006* |
Temperature | 1, 38 | 0.43 | 0.519 | |
Exposure time C | 1, 204 | 7.85 | 0.006* | |
Racer goby vs European bullhead from the mixed-species treatment | SpeciesWS | 1, 295 | 0.33 | 0.565 |
Temperature | 1, 25 | 0.75 | 0.395 | |
Exposure timeC | 1, 295 | 13.78 | <0.001* | |
Racer goby from mixed- vs single-species treatments | Treatment | 1, 46 | 0.001 | 0.976 |
Temperature | 1, 46 | 0.09 | 0.766 | |
Exposure timeC | 1, 244 | 13.40 | <0.001* | |
European bullhead from mixed- vs single-species treatments | Treatment | 1, 43 | 4.81 | 0.034* |
Temperature | 1, 43 | 3.81 | 0.057 | |
Exposure timeC | 1, 229 | 8.69 | 0.004* |
The numbers of interspecific aggression events displayed by the racer goby and European bullhead in the mixed-species treatment (Fig.
The racer goby showed similar levels of intra- and interspecific aggression (Fig.
In the single-species treatments (Fig.
Numbers of aggression acts per 10 min (a single observation period) shown by the monkey goby and gudgeon kept in separate single-species treatments (A) or together in the mixed-species treatment (B). Panels C, D present comparisons of the behaviour of the monkey goby and gudgeon, respectively, between the single- and mixed-species treatments. Symbols represent raw data (means ± 95%CI) for each species, temperature and period. Lines are predicted by the models (with 95%CI as shaded areas). Common slopes were predicted for groups of data that did not differ significantly from each other in the models. Horizontal lines indicate non-significant slopes.
General Linear Mixed Models to test the impact of treatment, temperature, exposure time and species on the number of aggressive events shown by the monkey goby and gudgeon. Non-significant higher order interactions were removed from the models in a simplification procedure.
Analysis | Effect | df | F | P |
---|---|---|---|---|
Monkey goby vs gudgeon from single-species treatments | Species | 1, 47 | 0.11 | 0.746 |
Temperature | 1, 47 | 0.06 | 0.806 | |
Exposure timeC | 1, 249 | 13.73 | <0.001* | |
Monkey goby vs gudgeon from the mixed-species treatment | SpeciesWS (Spec.) | 1, 305 | 17.58 | <0.001* |
Temperature | 1, 26 | 1.06 | 0.313 | |
Exposure timeC (Time) | 1, 305 | 17.28 | <0.001* | |
Spec.WS*Time | 1, 305 | 12.18 | 0.001* | |
Monkey goby from mixed- vs single-species treatments | Treatment | 1, 50 | 1.88 | 0.177 |
Temperature | 1, 50 | 1.41 | 0.241 | |
Exposure timeC | 1, 264 | 9.05 | 0.003* | |
Gudgeon from mixed- vs single-species treatments | Treatment (Treat.) | 1, 279 | 9.42 | 0.002* |
Temperature (Temp.) | 1, 279 | 1.56 | 0.212 | |
Exposure time (Time)C | 1, 261 | 19.98 | <0.001* | |
Treat.*Temp. | 1, 279 | 6.06 | 0.014* | |
Treat.*Time | 1, 261 | 6.36 | 0.012* | |
Temp.*Time | 1, 261 | 1.28 | 0.259 | |
Treat.*Temp.*Time | 1, 261 | 5.12 | 0.024* |
The number of interspecific aggression events displayed by these fish in the mixed-species treatment (Fig.
The monkey goby displayed similar levels of intra- and interspecific aggression (single vs. mixed-species treatments) irrespective of temperature (Fig.
On the other hand, the number of aggression events shown by the gudgeon (Fig.
In the single-species treatments (Fig.
Times to enter the feeder by the racer goby and European bullhead kept in separate single-species treatments (A) or together in the mixed-species treatment (B). Arrows indicate groups significantly differing from each other.
Cox proportional hazard regression models to test the effect of treatment, temperature and species on the time to enter the feeder by the racer goby and European bullhead.
Analysis | Effect | df | χ2 | P |
---|---|---|---|---|
Racer goby vs European bullhead from single-species treatments | Species | 1 | 12.92 | <0.001* |
Temperature | 1 | 1.16 | 0.282 | |
Racer goby vs European bullhead from the mixed-species treatment | Species | 1 | 17.54 | <0.001* |
Temperature | 1 | 6.99 | 0.008* | |
Racer goby from mixed- vs single-species treatments | Treatment | 1 | 4.96 | 0.026* |
Temperature | 1 | 4.45 | 0.035* | |
European bullhead from mixed- vs single-species treatments | Treatment | 1 | 10.58 | <0.001* |
Temperature | 1 | 0.63 | 0.429 |
In the mixed-species treatment (Fig.
The racer goby reached the feeder earlier in the presence of conspecifics than with the European bullhead (Table
In the single-species treatments (Fig.
Times to enter the feeder by the monkey goby and gudgeon kept in separate single-species treatments (A) or together in the mixed-species treatment (B). Arrows indicate groups significantly differing from each other.
Cox proportional hazard regression models to test the effect of treatment, temperature and species on the time to enter the feeder by monkey goby and gudgeon.
Analysis | Effect | df | χ2 | P |
---|---|---|---|---|
Monkey goby vs gudgeon from single-species treatments | Species | 1 | 5.33 | 0.021* |
Temperature | 1 | 6.42 | 0.011* | |
Monkey goby vs gudgeon from the mixed-species treatment | Species | 1 | 4.74 | 0.029* |
Temperature | 1 | 0.18 | 0.671 | |
Monkey goby from mixed- vs single-species treatments | Treatment | 1 | 0.56 | 0.454 |
Temperature | 1 | 0.10 | 0.754 | |
Gudgeon from mixed- vs single-species treatments | Treatment | 1 | 1.51 | 0.219 |
Temperature | 1 | 7.84 | 0.005* |
In the mixed-species treatment (Fig.
The entry time to the feeder shown by the monkey goby and gudgeon was independent of the species identity of the other individual in the pair (Table
In the single-species treatments (Fig.
Times spent in the feeder (as percentage of the total exposure time) by the racer goby and European bullhead kept in separate single-species treatments (A). Panels B, C present comparisons of the behaviour of the racer goby and European bullhead, respectively, between the single- and mixed-species treatments. Symbols represent raw data (means ± 95%CI) for each species, temperature and period. Lines are predicted by the models (with 95%CI as shaded areas). Common slopes were predicted for groups of data that did not differ significantly from each other in the models. Horizontal lines indicate non-significant slopes.
General Linear Mixed Models to test the impact of treatment, temperature, exposure time and species on the feeder occupancy time shown by the racer goby and European bullhead. Non-significant higher order interactions were removed from the models in a simplification procedure.
Analysis | Effect | df | F | P |
---|---|---|---|---|
Racer goby vs European bullhead from single-species treatments | Species (Spec.) | 1, 220 | 14.44 | <0.001* |
Temperature | 1, 38 | 6.74 | 0.013* | |
Exposure time (Time)C | 1, 203 | 55.08 | <0.001* | |
Spec.*Time | 1, 203 | 13.30 | <0.001* | |
Racer goby from mixed- vs single-species treatments | Treatment | 1, 46 | 2.43 | 0.126 |
Temperature | 1, 46 | 2.95 | 0.093 | |
Exposure timeC | 1, 244 | 103.98 | <0.001* | |
European bullhead from mixed- vs single-species treatments | Treatment | 1, 43 | 1.32 | 0.258 |
Temperature (Temp.) | 1, 246 | 9.10 | 0.003* | |
Exposure time (Time)C | 1, 228 | 3.19 | 0.075 | |
Temp.*Time | 1, 228 | 7.59 | 0.006* |
The presence of heterospecifics in the mixed-species treatment (Fig.
In the single-species treatments (Fig.
Times spent in the feeder (as percentage of the total exposure time) by the monkey goby and gudgeon kept in separate single-species treatments (A). Panels B, C present comparisons of the behaviour of the monkey goby and gudgeon, respectively, between the single- and mixed-species treatments. Symbols represent raw data (means ± 95%CI) for each species, temperature and period. Lines are predicted by the models (with 95%CI as shaded areas). Common slopes were predicted for groups of data that did not differ significantly from each other in the models. Horizontal lines indicate non-significant slopes.
General Linear Mixed Models to test the impact of treatment, temperature, exposure time and species on the feeder occupancy time shown by the monkey goby and gudgeon. Non-significant higher order interactions were removed from the models in a simplification procedure.
Analysis | Effect | df | F | P |
---|---|---|---|---|
Monkey goby vs gudgeon from single-species treatments | Species (Spec.) | 1, 264 | 27.41 | <0.001* |
Temperature (Temp.) | 1, 264 | 20.27 | <0.001* | |
Exposure Time (Time)C | 1, 247 | 91.76 | <0.001* | |
Spec.*Time | 1, 247 | 24.71 | <0.001* | |
Temp.*Time | 1, 247 | 17.78 | <0.001* | |
Monkey goby from mixed- vs single-species treatments | Treatment | 1, 50 | 8.28 | 0.006* |
Temperature | 1, 50 | 0.03 | 0.871 | |
Exposure timeC | 1, 264 | 108.14 | <0.001* | |
Gudgeon from mixed- vs single-species treatments | Treatment (Treat.) | 1, 290 | 0.06 | 0.803 |
Temperature (Temp.) | 1, 290 | 8.10 | 0.005* | |
Exposure time (Time)C | 1, 261 | 14.47 | <0.001* | |
Treat.*Temp. | 1, 290 | 3.58 | 0.060 | |
Treat.*Time | 1, 261 | 0.02 | 0.889 | |
Temp.*Time | 1, 261 | 6.28 | 0.013* | |
Treat.*Temp.*Time | 1, 261 | 4.00 | 0.046* |
The feeder was occupied for a longer time by the monkey goby in the presence of the gudgeon in the mixed-species treatment (6.9% of the total exposure time, Fig.
Whereas the gudgeon spent more time in the presence of conspecifics than with the monkey goby, but only at the beginning of exposure at 25 °C (6.7 vs. 4.0% of time during the first 20 min of the exposure, Fig.
Present work supported the first hypothesis that the non-native gobies are more successful food competitors than their native counterparts. Although the invaders did not consistently outperform the natives in terms of higher aggression, they revealed faster and longer food access compared to the natives. However, the second hypothesis was not confirmed. The effect of an elevated temperature on interspecific competition did not translate into a more apparent dominance of the gobies over the native fish.
In single-species treatments, the racer goby was more aggressive than the European bullhead. In contrast, the gudgeon and monkey goby did not differ in level of aggression in the second pair of co-existing species. This finding suggests that aggressive behaviour is of primary importance for establishing intra-species dominance in the racer goby, while it is not so in the other species tested. Interestingly, relatively higher aggression was revealed in inter-species interactions in these species in our study (
The relationship between temperature and aggression acts can be variable. Elevated temperature can either increase (
Regarding the first pair of fish studied, we found that the aggression shown by the European bullhead and racer goby was independent of temperature. This indicates the potential of the European bullhead to survive in warming waters, assuming they have access to abundant food resources and meet increased energy needs under such conditions. According to
We posit that the time to enter the feeder and the time spent in the feeder should be considered together. These two behaviours are likely to act together in the same direction to enhance the probability of success in food resource competition. Both gobies tended to reach the feeder before their native counterparts. This was likely to limit foraging of their native competitors and provided the invasives with better access to the richest food resources (directly after the food application), which has also been shown for larger (adult) European bullhead and racer goby (
An elevated temperature delayed the time to enter the feeder by the invasive racer goby, despite the fact that this species originates from a warmer climate than that in its invaded range. Hence, increased temperature may have a disruptive effect not only on natives, but also on invasive species. However, the native species, being less adapted to elevated temperatures, can use even more energy or even limit their foraging (thus causing difficulties with obtaining energy) at 25 °C. Therefore, in the longer term, indirect (exploitation) competition (
Our study has shown that, regardless of summer temperatures (normal or elevated) that occur in Central European rivers, the juvenile invasive gobies are more effective than their native counterparts in competing for access to limited food resources. This finding broadens the knowledge of the threat posed by the Ponto-Caspian gobies towards native European freshwater fishes (see a review by
The authors have declared that no competing interests exist.
We sampled and used the fish by the permit of the Local Committee for Ethics in Animal Research in Bydgoszcz, Poland, statement no. 30/2022 from 22 June 2022. In addition, the capture and use of the European bullhead, which is protected by law in Poland, was accepted by the Regional Directorate of Environmental Protection in Bydgoszcz, Poland (approval number: WOP.6401.4.52.2022.MO). The procedures conducted within the study met the ASAB/ABS guidelines for the use of animals in research (ASAB Ethical Committee and ABS Animal Care Committee 2019). The housing conditions guaranteed a high level of animal welfare, which was manifested by the overall activity and food intake of the fish throughout the research period. The fish had no external signs of stress or disease. Each specimen was used only once. After the experiments, the European bullhead and gudgeon were released into the wild from which they were caught. In turn, racer goby and monkey goby, as invasive species were euthanized by an overdose of buffered Tricaine Methanesulfonate (MS-222) and disposed of according to the Regulation of the Polish Minister of the Environment from 9 September 2011 (Journal of Laws No. 210, item 1260). Killing was conducted by a qualified, certified person (certificate No. 2355/2015 issued by the Polish Laboratory Animal Science Association).
This research was supported by “The Excellence Initiative – Research University” programme at Nicolaus Copernicus University in Toruń, Poland.
PK: Conceptualisation, Resources, Methodology, Investigation, Formal analysis, Data interpretation, Visualisation, Writing-Original draft preparation. JK: Conceptualisation, Resources, Formal analysis, Data interpretation, Visualisation, Writing-Review & Editing. TK: Conceptualisation, Resources, Methodology, Data interpretation, Writing-Review & Editing, Supervision.
Piotr Kłosiński https://orcid.org/0000-0002-6359-9781
Jarosław Kobak https://orcid.org/0000-0001-7660-9240
Tomasz Kakareko https://orcid.org/0000-0001-5054-7416
The data that support the findings of this study are available from the corresponding author, PK, upon reasonable request.
Experimental setup
Data type: tif
Explanation note: Experimental setup (all the dimensions are given in mm).
Experimental procedure
Data type: tif
Explanation note: Experimental procedure.
Numbers of replicates and individual size differences in particular experimental treatments
Data type: docx
Explanation note: Numbers of replicates (n) and individual size (total length, TL) differences in particular experimental treatments (pairs of fish in particular species configurations tested at specific temperatures). Individual sizes were compared between the species within each species pair in each experimental treatment using t-tests for dependent samples.
The summary of all the models used in the study
Data type: docx
Explanation note: The summary of all the models used in the study.
Tests of slope significance for particular levels of categorical factors interacting with exposure time in their effects on the number of aggression events shown by the monkey goby and gudgeon (see Table
Data type: docx
Explanation note: Tests of slope significance for particular levels of categorical factors interacting with exposure time in their effects on the number of aggression events shown by the monkey goby and gudgeon (see Table
Tests of slope significance for particular levels of categorical factors interacting with exposure time in their effects on the time spent in the feeder by the racer goby and European bullhead (see Table
Data type: docx
Explanation note: Tests of slope significance for particular levels of categorical factors interacting with exposure time in their effects on the time spent in the feeder by the racer goby and European bullhead (see Table