Research Article |
Corresponding author: Fabio Ercoli ( fabio.ercoli@emu.ee ) Academic editor: Emili García-Berthou
© 2024 Fabio Ercoli, Mikko Kiljunen, Paul Teesalu, Arvo Tuvikene, Meelis Tambets, Einar Kärgenberg, Tiina Nõges.
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:
Ercoli F, Kiljunen M, Teesalu P, Tuvikene A, Tambets M, Kärgenberg E, Nõges T (2024) Niche partitioning of invasive Amur sleeper (Perccottus glenii) amongst native fish communities in three different freshwater ecosystems. NeoBiota 95: 181-198. https://doi.org/10.3897/neobiota.95.116327
|
The invasive fish, Amur sleeper, poses a significant and growing threat to Central European freshwater ecosystems. Despite its rapid spread, the ecological implications of its invasion have been poorly explored. Recent findings confirm its presence in various Estonian freshwater systems, raising concerns about its imminent expansion into larger lakes. To better understand its potential ecological impacts, we explored the isotopic niche of the Amur sleeper in comparison with native fish species co-existing in three Estonian freshwater ecosystems. We employed carbon and nitrogen stable isotope analyses alongside gut content analyses. Our findings show that the Amur sleeper’s diet in newly-invaded Estonian water bodies predominantly comprises benthic macroinvertebrates, although it may also include fish, confirming its role as a predator in the local food web. Notably, Amur sleeper populations exhibited clear isotopic niche partitioning in three invaded ecosystems. A logistic regression model, based on stomach content analyses, revealed an ontogenetic diet shift from benthivorous to piscivorous feeding habits from small to large specimens. Amur sleeper exhibits voracious, non-selective feeding habits, which can negatively impact native freshwater communities. The ability to occupy a distinct isotopic niche, with minimal overlap with native fish populations, may reduce interspecific competition, facilitating the spread and establishment of Amur sleeper in newly-invaded habitats. Managing the spread of this invasive species thus becomes even more critical to safeguard the integrity of native aquatic ecosystems.
gut content, invasive species, isotopic niche, native fish community, ontogenetic shift, stable isotopes
The spread of invasive alien species is recognised as one of the greatest threats to global biodiversity (
Amur sleeper invasion could alter trophic relationships in aquatic communities by generalist feeding on different trophic levels and prey items mainly associated with submerged vegetation (
Amur sleeper has already been reported in Estonian freshwater ecosystems (
In summer 2008, northern pike (Esox lucius Linnaeus, 1758) were introduced in the pond to eradicate Amur sleeper. Intensive capture and predation by pike and perch reduced the population biomass of Amur sleeper by 80% (
Due to its feeding habits, ability to exploit a wide spectrum of food resources and rapid invasion rate, it is crucial to better understand the potential ecological effects of the Amur sleeper on native communities for impact assessment and risk management. In this study, we aimed to: (1) evaluate the diet plasticity of invasive Amur sleeper, (2) assess if the species displays a similar isotopic niche in each of the three different invaded ecosystems and (3) quantify its potential overlap with native fish species. We hypothesised that: 1). Amur sleeper has a wider isotopic niche than benthic native fish and 2). Amur sleeper is a strong competitor for local predatory-benthivorous fish species (e.g. perch) due to niche overlap and a similar diet.
Samples were collected from the Narva Reservoir (59°21.02'N, 28°10.79'E), the cooling water channel of the Baltic Thermal Power Plant (BPP) (59°18.94'N, 28°5.12'E) and a small pond located near, but not connected, to the reservoir (59°21.19'N, 28°11.17'E) (Fig.
Samples of Amur sleeper and native fish species were collected in September 2017 from the reservoir, pond and channel using electrofishing and gillnets. Nordic-type, multi-section benthic and pelagic gillnets (5–55 mm from knot to knot) were supplemented with larger mesh sized nets (65 mm from knot to knot). Gillnets were set up overnight in each sampling site and retrieved approximately 15 hours later the next day. Electrofishing (mean area 618 m2) was performed in shallow water with abundant vegetation near the areas where gillnets were deployed, at the same time or after their removal. The pond was also sampled prior to 2017 and in 2018, using only electrofishing. In September 2017, three replicate samples for macroinvertebrates were collected using a kick-net with a mesh size of 0.5 mm (0.5 to 1.5 m water depth), along littoral areas at sites near where gillnets were set and removed and electrofishing was performed. Additionally, Ekman grab was used to sample macroinvertebrates from the bottom of the pond.
After sampling, fish and macroinvertebrates were identified to species and at the lowest taxonomic levels, respectively. For stable isotopes analysis (SIA), a piece of white dorsal muscle was collected from each specimen, while macroinvertebrate samples were prepared from whole organisms to represent macroinvertebrate food sources for Amur sleeper (Table
Number (n) of macroinvertebrates species sampled in pond, Narva Reservoir and channel and their respective mean values (± standard deviation) of carbon and nitrogen stable isotopes.
Ecosystem | Species | n | d13C‰ | d15N‰ |
---|---|---|---|---|
Pond | Gastropoda | |||
Bithynia tentaculata | 3 | -29.58 ± 0.66 | 3.82 ± 0.46 | |
Crustacea | ||||
Asellus aquaticus | 5 | -25.26 ± 0.63 | 3.35 ± 0.34 | |
Ephemeroptera | ||||
Ephemerellidae | 5 | -26.92 ± 0.79 | 2.2 ± 0.27 | |
Coleoptera | ||||
Dytiscidae | 6 | -25.55 ± 2.86 | 4.81 ± 0.65 | |
Bivalvia | ||||
Sphaerium sp | 9 | -32.22 ± 0.31 | 3.98 ± 0.36 | |
Narva Reservoir | Gastropoda | |||
Radix balthica | 4 | -27.67 ± 1.00 | 8.95 ± 0.44 | |
Crustacea | ||||
Asellus aquaticus | 6 | -26.70 ± 0.27 | 7.73 ± 1.43 | |
Gmelinoides fasciatus | 16 | -25.88 ± 0.75 | 8.43 ± 1.17 | |
Bivalvia | ||||
Dreissena polymorpha | 4 | -30.75 ± 0.23 | 10.16 ± 0.08 | |
Odonata | ||||
Coenagrionidae | 3 | -29.25 ± 0.12 | 11.13 ± 0.07 | |
Trichoptera | ||||
Phryganeidae | 3 | -25.40 ± 3.04 | 10.69 ± 0.47 | |
BPP channel | Gastropoda | |||
Lymnaea stagnalis | 3 | -21.99 ± 0.71 | 8.18 ± 0.16 | |
Crustacea | ||||
Procambarus virginalis | 6 | -29.09 ± 1.29 | 8.28 ± 0.62 | |
Ephemeroptera | ||||
Caenis horaria | 2 | -31.02 ± 1.01 | 6.94 ± 0.12 | |
Bivalvia | ||||
Dreissena polymorpha | 3 | -27.70 ± 0.73 | 9.3 ± 0.31 | |
Odonata | ||||
Aeshnidae | 3 | -26.66 ± 0.12 | 9.57 ± 0.02 | |
Trichoptera | ||||
Phryganeidae | 3 | -28.38 ± 0.66 | 8.84 ± 0.19 |
All SIA samples were freeze-dried for 48 hours to constant weight, ground to a fine, homogenous powder, weighed into tin cups (~ 0.6 mg of material) and encapsulated. If needed, small macroinvertebrates of the same taxa were pooled to achieve enough material for SIA. Analyses of stable carbon and nitrogen isotopes were conducted using a Thermo Finnigan DELTAplus Advantage continuous-flow isotope ratio mass spectrometer coupled to a Flash EA 1112 elemental analyser (Thermo Scientific, Bremen, Germany) at Jyväskylä University (Finland). SI values for carbon and nitrogen are expressed as parts per thousand (‰) delta values (δ13C, δ15N) relative to international standards:
X = (Rsample / Rstandard -1) × 1000
where X is either carbon or nitrogen SI value and R is the ratio of heavy to light SI of carbon or nitrogen in samples and standards.
Reference materials were used as internal standards, with known relationships to the international standards of The Vienna Pee Dee Belemnite (VPDB) for carbon SI and atmospheric nitrogen (Air-N2) for nitrogen SI. White muscle tissue of northern pike (Esox lucius) and birch leaves (Betula pendula L.), with known isotopic compositions, were used as internal working standards to ensure analytical precision. One internal standard was run repeatedly after every five samples in each sequence. Standard deviations within reference samples in each analytical run were always less than 0.1‰ for carbon and 0.2‰ for nitrogen in pike and birch leaf samples. Sample analysis also yielded percentage carbon and nitrogen from which C:N ratios (by weight) were derived.
Stable isotope values were used for identifying isotopic niches and to evaluate and compare isotopic niche similarity of Amur sleeper and native fish species in each of the studied freshwater ecosystems. Isotopic niches and isotopic niches overlapping calculations were performed using the SIBER (Stable Isotope Bayesian Ellipses in R) package in R (
Amur sleeper is a predatory, omnivorous fish and its diet was assessed by stomach content analyses to support SIA interpretations. Stomach contents were identified and divided into the following categories: macrophytes, zooplankton, macroinvertebrate parts, gastropods, bivalves, Chironomidae, Trichoptera, Ephemeroptera, Asellus aquaticus, Coleoptera, Odonata, fish, Oligochaeta, Gammaridae, Micronecta and detritus. Stomach fullness was estimated, ranging from 0% (empty stomach) to 100% (full stomach), based on the sum of percentages of each food category (
where Ji is the number of Amur sleeper individuals containing prey i and P is the number of Amur sleeper individuals with food in their stomach (
To assess the probability of Amur sleeper undergoing an ontogeny diet shift, a logistic regression model was used to estimate the length when individuals could switch their diet from benthivorous to piscivorous feeding habits
y = [e(α+βx)] [1 + e(α+βx)]−1
where y indicates the occurrence of fish in Amur sleeper stomach, x is the length of Amur sleeper specimens and α and β are the coefficients estimated by the model.
From stomach content analyses, individuals with fish in their stomachs were given the value of 1, whereas those without fish in their stomach were given the value of 0. The logistic regression model was performed in R (
A total of 66 Amur sleeper individuals were caught: 17 from the pond (4 large and 13 small), 36 from the Narva Reservoir (8 large and 28 small) and 13 from the BPP channel (1 large and 12 small) (Table
Number of Amur sleeper and other fish species individuals used for stable isotope analyses and mean values of their length (cm) and carbon and nitrogen stable isotope values (± standard deviation). SEAc indicates the Standard Ellipse Area corrected of each fish population isotopic niche in the pond, Narva Reservoir, and BPP channel. SEAc of Amur sleeper is calculated including large and small specimens. SEAc overlapping % represents the degree of isotopic niche similarity between Amur sleeper and other fish species in each ecosystem.
Ecosystem | Fish species | n | Length (cm) | d13C‰ | d15N‰ | SEAc‰2 | SEAc overlapping% |
---|---|---|---|---|---|---|---|
Pond | Amur sleeper (Perccottus glenii) large | 4 | 14 ± 5.2 | -38.19 ± 0.59 | 9.58 ± 1.82 | 77.5 | – |
Amur sleeper (Perccottus glenii) small | 13 | 5.6 ± 1.0 | -38.16 ± 10.61 | 6.79 ± 2.48 | – | – | |
Roach (Rutilus rutilus) | 10 | 20.0 ± 5.1 | -38.17 ± 0.95 | 7.88 ± 0.23 | 0.43 | 0 | |
Rudd (Scardinius erythrophthalmus) | 5 | 10.9 ± 2.8 | -26.47 ± 1.12 | 8.10 ± 0.10 | 0.49 | 0 | |
Perch (Perca fluviatilis) | 9 | 15.9 ± 2.9 | -26.06 ± 1.12 | 8.80 ± 0.37 | 1.51 | 0 | |
Pike (Esox lucius) | 4 | 17.2 ± 3.0 | -26.33 ± 0.14 | 8.08 ± 0.46 | 0.27 | 0 | |
Narva Reservoir | Amur sleeper (Perccottus glenii) large | 8 | 11.5 ± 2.8 | -27.04 ± 0.28 | 12.86 ± 0.27 | 0.88 | – |
Amur sleeper (Perccottus glenii) small | 28 | 5.7 ± 2.0 | -27.13 ± 0.52 | 11.83 ± 0.38 | – | – | |
Roach (Rutilus rutilus) | 11 | 20.8 ± 2.0 | -26.46 ± 1.68 | 13.27 ± 0.75 | 4.81 | 14 | |
Rudd (Scardinius erythrophthalmus) | 8 | 22.1 ± 2.6 | -24.10 ± 1.86 | 11.69 ± 1.53 | 10.7 | 0 | |
Perch (Perca fluviatilis) | 12 | 23.4 ± 3.4 | -26.62 ± 0.66 | 14.89 ± 0.67 | 1.53 | 0 | |
Ruffe (Gymnocephalus cernua) | 8 | 6.8 ± 2.5 | -30.02 ± 2.36 | 11.78 ± 1.45 | 5.19 | 0 | |
Tench (Tinca tinca) | 4 | 31.3 ± 8.7 | -27.19 ± 0.81 | 12.87 ± 0.62 | 2.09 | 35 | |
Bleak (Alburnus alburnus) | 5 | 11.3 ± 0.8 | -29.31 ± 1.28 | 12.61 ± 0.50 | 1.66 | 0 | |
Spined loach (Cobitis taenia) | 5 | 9.4 ± 0.3 | -27.78 ± 0.42 | 12.37 ± 0.28 | 0.50 | 32 | |
White bream (Blicca bjoerkna) | 6 | 12.6 ± 1.6 | -26.08 ± 0.64 | 14.01 ± 0.57 | 1.71 | 0 | |
BPP channel | Amur sleeper (Perccottus glenii) large | 1 | 9.2 | -27.64 | 11.35 | 0.59 | – |
Amur sleeper (Perccottus glenii) small | 12 | 4.1 ± 1.6 | -29.0 ± 0.45 | 10.52 ± 0.48 | – | – | |
Roach (Rutilus rutilus) | 5 | 12.1 ± 1.4 | -27.04 ± 1.13 | 12.65 ± 0.91 | 4.95 | 0 | |
Rudd (Scardinius erythrophthalmus) | 5 | 23.1 ± 3.7 | -23.61 ± 1.33 | 11.74 ± 0.84 | 5.93 | 0 | |
Perch (Perca fluviatilis) | 6 | 19.3 ± 3.5 | -26.46 ± 1.39 | 13.28 ± 0.17 | 0.93 | 0 | |
Spined loach (Cobitis taenia) | 5 | 7.4 ± 1.1 | -28.74 ± 0.71 | 12.24 ± 0.54 | 1.46 | 0 |
Catch Per Unit Effort (CPUE) indicates the number of fish per 100 m2 of sampling area in electrofishing and the number of fish per net in gillnet sampling, caught in pond, Narva Reservoir and BPP channel.
Place and sampling type | Amur sleeper | Perch | Pike | Roach | Gibel carp | Rudd | Ruffe | Spined loach | Bleak | Tench | White bream |
---|---|---|---|---|---|---|---|---|---|---|---|
Pond | |||||||||||
Electrofishing | 0.2 | 0 | 0.3 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 |
Benthic Nordic | – | 18.3 | 0.33 | 4.3 | 0 | 3 | 0 | 0 | 0 | 0 | 0 |
Pelagic Nordic | – | 11 | 0 | 6 | 0 | 6 | 0 | 0 | 0 | 0 | 0 |
Benthic 65 mm (knot to knot) | – | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 |
Narva Reservoir | |||||||||||
Electrofishing | 6 | 0 | 0.1 | 0 | 0 | 0 | 0 | 0.3 | 0 | 0.5 | 0 |
Benthic Nordic | – | 10.3 | 0 | 6.3 | 0 | 0.7 | 6.3 | 0.7 | 0 | 0.3 | 0.3 |
Pelagic Nordic | – | 0 | 0 | 20 | 0 | 1 | 0 | 0 | 7 | 0 | 0 |
Benthic 65 mm (knot to knot) | – | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 0 | 2 | 0 |
BPP channel | |||||||||||
Electrofishing | 43 | 0 | 0.1 | 0 | 0 | 0 | 0 | 13 | 0.4 | 0 | 0 |
Benthic Nordic | – | 10.3 | 0 | 11.3 | 0 | 10.7 | 0 | 0.3 | 0.3 | 0 | 1 |
Benthic 65 mm (knot to knot) | – | 0 | 0 | 0 | 4 | 0 | 0 | 0 | 0 | 0 | 0 |
SIA results show that values of nitrogen stable isotope of Amur sleeper individuals differ between the pond and reservoir and between the pond and channel, being highest in the reservoir and lowest in the pond (Fig.
Carbon and nitrogen stable isotopes mean values (± standard deviation) of Amur sleeper large and small individuals (open and filled black circles), native fish species and macroinvertebrates, in the (A) pond, (B) Narva Reservoir and (C) BPP channel.
Isotopic niches of each fish species represented by ellipses and fish individuals represented by open circles in the (A) pond, (B) Narva Reservoir and (C) BPP channel.
SIBER model results indicated that Amur sleeper isotopic niches did not overlap with those of co-occurring native fish populations, except in Narva Reservoir, where it overlapped with tench (35%), spined loach (32%) and roach (14%). The isotopic niche (SEAc) of Amur sleeper was very wide in the pond (77.5‰2) due to low carbon isotope values, while the niches were smaller in the channel and Reservoir (0.59‰2 and 0.88‰2, respectively), compared to other fish species within each system. Only spined loach in the reservoir occupied a smaller niche area compared to Amur sleeper (Fig.
The logistic regression curve indicated that Amur sleeper undergoes an ontogenetic shift from < 50% benthivorous to > 50% piscivorous diet at an individual total length (TL) of 8.7 cm (p = 0.004; 95% CI, Fig.
Logistic regression curve was fitted by the following equation y = [e(−13·226 + 0·426x)][1+ e(−13·226 + 0·426x)]−1 showing ontogenetic diet shift from benthivory (0) to piscivory (1). Shadow area represents 95% Confidence Interval.
Stomachs of all sampled Amur sleeper individuals (66) were analysed, of which 65 were 100% full, regardless of fish size and one was 90% full. The amount of the analysed stomachs of large and small individuals were respectively 4 and 13 from pond; 8 and 28 from Narva Reservoir and 1 and 12 from the channel. In all the invaded ecosystems, the Amur sleeper diet comprised mainly macroinvertebrates and, to some extent, fish. In the channel, only one individual represented the large group (TL > 8.7 cm), its stomach being 100% full of macrophytes and detritus. In all studied ecosystems, the diet of small (TL < 8.7 cm) specimens consisted mostly of macroinvertebrates (fragments), ranging from 17% in the channel to 46% in the pond, with Chironomidae, Trichoptera, Ephemeroptera and Odonata being the most common (Fig.
Invasive Amur sleeper was well-established in each of the studied waterbodies and occupied the isotopic niche between predatory-omnivorous fish species and macroinvertebrates. In the reservoir and channel, the nitrogen isotope values of Amur sleeper, both large and small, show its intermediate role between bottom-up and top-down energy fluxes of the food web, while in the pond, large individuals had similar isotope values with native fish species. Yet, in the pond, the isotopic niche of Amur sleeper was much wider than those of co-occurring native fish. Substantial variations in nitrogen and carbon isotope values and much lower values for Amur sleeper compared to other fish, contributed to the wide isotopic niche in the pond. Low carbon SI values may indicate anoxic conditions at the pond bottom, favouring methanogenesis and conversion of methane to microbial biomass by methane-oxidising bacteria (MOB) (
Chironomidae and Oligochaeta feeding on MOB would be the best prey candidates for low carbon SI values. Unfortunately, we were not able to analyse SI values from these organisms, since the bottom of the pond was fully covered by macrophytes, preventing sampling of profundal benthic macroinvertebrates. However, amongst the sites and both size groups, highest proportions of Oligochaeta were found from stomachs of small Amur sleeper in the pond, suggesting that Oligochaeta was an important prey item and potential source for methane-derived carbon.
Invasive Amur sleeper is a predatory-omnivorous fish species mainly feeding on macroinvertebrates, but can be also piscivorous (
In accordance with prior research (
Pike were introduced in the pond to help control the Amur sleeper population. Similar to the findings of
Amongst the prey items found in Amur sleeper stomachs, spined loach (Cobitis taenia) was identified. Spined loach is a protected fish species in Estonia (
Interestingly, amongst species consumed by Amur sleeper, there were also different invasive species present, such as zebra mussels (Dreissena polymorpha), amphipod (Gmelinoides fasciatus) and marbled crayfish (Procambarus virginalis), which are all well established in the studied waterbodies. All these invasive species occupy different trophic positions and could represent food sources for Amur sleeper, possibly facilitating Amur sleeper invasion and establishment. These results support the invasion meltdown hypothesis (
In our studied ecosystems, Amur sleeper can affect native communities at different trophic levels, especially native macroinvertebrate communities. These feeding characteristics exert both top-down and bottom-up effects, disrupting energy transfer to higher trophic levels. In a recent study,
Our findings showed that Amur sleeper may affect native fish communities in the studied ecosystems, serving as a warning for stakeholders and resource managers of the potential threats that this invasive fish species can pose on invaded and nearby, connected aquatic habitats.
We thank Dr Mark McCarthy for proofreading the manuscript.
The authors have declared that no competing interests exist.
No ethical statement was reported.
This project has received funding from the European Union’s Horizon 2020 Research and Innovation Programme under grant agreement No 951963. The project was financed by the Estonian University of Life Sciences (Project Number P190254PKKH to FE) and by the Estonian Research Council (Mobilitas Pluss project MOBJD29 to FE). The Estonian Research Council grant PRG709, the Climate Change Mitigation and Adaptation Programme were financed by the European Economic Area Financial Mechanism. Additional funding came from the Environmental Investment Centre projects “Implementation of the action plan for the use and protection of crayfish stocks” and “Eradication of aquatic invasive species in Estonian freshwaters” No 4-17/16674.
FE and MK- Conceptualisation. MT, EK and FE - Investigation. FE, MK, PT and AT- Formal analysis, data curation, methodology, software and visualisation. FE, MK, PT and TN - Writing original draft. FE, MK, PT, EK, MT, AT and TN - Writing, review and editing.
Fabio Ercoli https://orcid.org/0000-0003-2433-6222
Mikko Kiljunen https://orcid.org/0000-0002-7411-1331
Paul Teesalu https://orcid.org/0000-0003-1713-4876
Arvo Tuvikene https://orcid.org/0000-0001-5520-1166
Meelis Tambets https://orcid.org/0000-0003-1620-5396
Einar Kärgenberg https://orcid.org/0000-0003-4391-1368
Tiina Nõges https://orcid.org/0000-0002-4558-7373
All of the data that support the findings of this study are available in the main text.