Research Article |
Academic editor: Jonathan M. Jeschke
© 2024 Marek Šmejkal, Kiran Thomas, Zuzana Šmejkalová, Yevdokiia Stepanyshyna, Daniel Bartoň, Sandip Tapkir, Travis Meador, Mojmír Vašek.
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:
Šmejkal M, Thomas K, Šmejkalová Z, Stepanyshyna Y, Bartoň D, Tapkir S, Meador T, Vašek M (2024) Isotopic niches reveal the impact of topmouth gudgeon and gibel carp on native crucian carp. NeoBiota 93: 203-224. https://doi.org/10.3897/neobiota.93.119274
|
Invasive species pose a major threat to natural ecosystems and directly outcompete many native species, placing them at imminent threat of extinction. The topmouth gudgeon (Pseudorasbora parva) is on the EU’s blacklist of invasive freshwater species and threatens biodiversity, especially in wetland and floodplain ecosystems, aquacultures and village ponds. The crucian carp (Carassius carassius) is native to Europe and its populations have declined in large part of its native range, with invasive gibel carp (C. gibelio) suspected as a major cause of its decline. Invasions by topmouth gudgeon have been implicated in the decline of crucian carp populations but this still needs to be verified. The aim of this study was to evaluate by the experimental approach the competitive interaction between the two species, topmouth gudgeon and crucian carp, focusing on isotopic niche sizes and their overlap in syntopy. A four-month mesocosm experiment was performed to determine the isotopic niche of crucian carp and topmouth gudgeon living alone and in syntopy. Additionally, stable isotope data were collected at the sites where the two species co-occur to compare niche sizes and overlaps. Experimental data showed that the isotopic space of topmouth gudgeon responded more flexibly (reducing niche size at syntopy) than that of the crucian carp and confirmed a high isotopic niche overlap between the species. Field studies have shown that topmouth gudgeon has invaded the isotopic niche of the crucian carp, especially when another invasive species, the gibel carp, lived in the community (25% at 40% ellipse area and 50% at 95% ellipse area). When only the topmouth gudgeon and crucian carp were present in the field, the overlap was lower (3% and 48%, respectively) and directional overlap modelling showed that the crucian carp was more likely to invade the isotopic niche of topmouth gudgeon than vice versa. The data indicated that competition between crucian carp and topmouth gudgeon is likely, especially in syntopy with other invasive species. This study shows that the feeding plasticity of topmouth gudgeon likely facilitates its establishment outside its native range and, due to high isotopic niche overlap, threatens native fish with similar feeding ecology with competitive displacement.
Biological invasions, freshwater conservation, interspecific competition, invasive species, SIA, stable isotope analysis
Wetland and floodplain ecosystems are amongst the most threatened freshwater ecosystems (
Fish from the Cyprinidae family were the first to be spread by humans outside their native ranges due to aquaculture, angling and ornamental fish breeding purposes and now pose a serious problem for biodiversity worldwide (
Stable isotope analysis has become a useful tool to assess the competitive relationship between invasive and native species (
There is a general trend in decline of native fish species of floodplain ecosystems in Europe (
The extent of isotopic niche overlap between the invasive topmouth gudgeon and the declining crucian carp was investigated in the mesocosm experiment, with an aim to estimate the shift in isotopic niche position in syntopy compared to allopatry. In addition, two populations of crucian carp and topmouth gudgeon living in syntopy were sampled in the field together with four populations, where the gibel carp was present as a second invasive species.
Specifically, we hypothesised that:
A mesocosm experiment was designed to test the size of the niches under conditions of syntopy and allopatry between the invasive topmouth gudgeon and the native crucian carp. The experiment was conducted using an array of 24 outdoor mesocosms at the Institute of Hydrobiology, Biology Centre of the Czech Academy of Sciences, České Budějovice, Czechia. To provide near-natural conditions, zooplankton and phytoplankton were collected from a nearby pond in České Budějovice (48.9671469°N, 14.4528125°E) and transferred to two separate 500 litre circular vats kept outdoors. The phytoplankton and zooplankton inoculum were kept in the vats filled with water from the fish ponds for two days. Each mesocosm (rectangular shape, 98 cm length × 91 cm height × 90 cm width, made of high-density polyethylene) received 800 litres of aged tap water and 10 kg of sand. Two litres of the well-mixed phytoplankton and zooplankton inoculum were then added to each tank on 19 April 2023. The set-up was maintained for one month before the fish were added to stabilise the system for phytoplankton and zooplankton production.
The fish for the mesocosm experiment were collected on 7 May 2023 in the Smilovice and Ujezdec ponds (49.8669103°N, 14.9854556°E and 49.9006147°N, 14.9486450°E; topmouth gudgeon and crucian carp, respectively). The fish were transported to České Budějovice and kept in 300 l aquaria under laboratory conditions at 20 °C and a photoperiod of 12L:12D for a 12-day acclimatisation and antiparasitic treatment (Costapur, Sera, Heinsberg, Germany). They were fed 1% of body weight in dry food (C-3 Carpe F, Skretting, Stavanger, Norway) per day.
To determine the overlap of isotopic niches under single species and syntopy conditions, three treatments (crucian carp, topmouth gudgeon and syntopy) were distributed in a randomised block design with eight replicates per treatment. This arrangement was chosen to obtain a sufficient number of replicated mesocosms and a sufficient number of individual fish in each treatment. Length and weight of individual fish (0.1 g accuracy, SI-132-3 balance, Excell, New Taipei City, Taiwan) were recorded at the beginning and end of the experiment. The initial fish biomass in each tank was roughly 8 g, depending on fish size variability (mean ± standard deviation (hereafter SD), crucian carp = 7.42 ± 0.62 g; topmouth gudgeon = 8.22 ± 0.71 g, SL = 88 ± 12 mm; syntopy = 8.22 ± 0.39 g). Crucian carp and topmouth gudgeon used for the experiment had an SL = 35 ± 6 mm and a weight of 1.29 ± 0.73 g and a SL = 44 ± 3 mm and a weight of 1.36 ± 0.24, respectively. Six fish were used in each mesocosm (6 crucian carp, 6 topmouth gudgeon or 3 and 3 in syntopy treatment, i.e. a total of 72 fish per species). The initial biomass of 8 g corresponds to a population density of 80 kg.ha-1, allowing for growth up to the carrying capacity of the mesocosm. The experiment ran for four months, from 19 May to 19 September 2023 and the mesocosms were monitored for transparency at the beginning, middle and end of the experiment using a Sneller tube (
At the end of the experiment, a total of 80 fin clips were collected (20 topmouth gudgeon and 20 crucian carp living in the syntopy treatment and 20 individuals of each species living in the single species treatment, 2–3 fish per species and mesocosm unit). As young individuals with a high growth rate were used, the stable isotope turnover during the warmest four months of the growing season should ensure a reasonable estimate of the individual isotopic niche position (
Fish were collected at six sites in the Czechia (small lentic waters with surface area between 0.01 and 0.20 ha, Table
Sampling sites involved in the study, selected on the basis of the presence of crucian carp (Carassius carassius) and topmouth gudgeon (Pseudorasbora parva) in the syntopy, with additional presence of gibel carp (C. gibelio) at four sites. CC = crucian carp, TMG = topmouth gudgeon, GC = gibel carp. The site ID corresponds to Fig.
Site ID | Sites | Macrophytes (%) | Latitude/ Longitude | Area (ha) | Sampling date | No. of individuals | Standard length (mm) | CPUE |
---|---|---|---|---|---|---|---|---|
1 | Soudkuv lom | 50 | 49.6280844°N, 15.3921969°E | 0.04 | 14/6/2023 | CC = 12 | CC = 105 ± 22 | CC = 0.16 |
TMG = 14 | TMG = 40 ± 3 | TMG = 0.41 | ||||||
GC = 15 | GC = 56 ± 25 | GC = 0.40 | ||||||
2 | Lipi | 100 | 48.9414672°N, 14.3579775°E | 0.01 | 21/6/2023 | CC = 15 | CC = 68 ± 10 | CC = 2.18 |
TMG = 15 | TMG = 56 ± 7 | TMG = 0.93 | ||||||
GC = 15 | GC = 48 ± 8 | GC = 1.59 | ||||||
3 | U Spacku | 100 | 48.9477278°N, 14.4882947°E | 0.11 | 23/6/2023 | CC = 15 | CC = 50 ± 6 | CC = 0.75 |
TMG = 12 | TMG = 48 ± 5 | TMG = 1.45 | ||||||
GC = 15 | GC = 56 ± 24 | GC = 0.95 | ||||||
4 | Vrabce | 30 | 48.9209244°N, 14.3631514°E | 0.06 | 28/6/2023 | CC = 15 | CC = 56 ± 9 | CC = 2.40 |
TMG = 11 | TMG = 32 ± 3 | TMG = 64.5 | ||||||
GC = 11 | GC = 109 ± 29 | GC = 0.10 | ||||||
5 | Smilovice | 10 | 49.8668344°N, 14.9854161°E | 0.13 | 8/5/2023 | CC = 15 | CC = 98 ± 10 | CC = 0.30 |
TMG = 15 | TMG = 46 ± 3 | TMG = 3.80 | ||||||
6 | Valcha | 55 | 49.7043211°N, 13.3039761°E | 0.20 | 14/7/2023 | CC = 6 | CC = 69 ± 9 | CC = 0.07 |
TMG = 15 | TMG = 48 ± 6 | TMG = 8.21 |
Map of sampling sites where specimens were collected for stable isotope analysis to estimate the trophic niche overlap of topmouth gudgeon (Pseudorasbora parva) and crucian carp (Carassius carassius). Species: Syn = topmouth gudgeon + crucian carp; syn+GC = topmouth gudgeon + crucian carp + gibel carp (C. gibelio). Site details are given in Table
Eighty fin clips from mesocosms from the experimental part of the study were subjected to stable isotope analysis (SIA). To compare the experimental results with the field data, scales from the field data were used, as scales are an alternative to dorsal muscle or fin clips for SIA and lead to comparable results (
A total of 216 sampled scales were selected from the six different sites (Table
To assess the position of the isotopic niches of the tested species in syntopy and allopatry, standard ellipse areas (40% EA) were used, which are a bivariate measure of the distribution of individuals in isotopic space, where the ellipse includes the mean 40% of the data as it represents the typical resource utilisation of the population (
A permutation analysis of variance with the adonis function (PERMANOVA, distance = Euclidean, nperm = 10,000) in the R package vegan was used to estimate the difference in syntopy between the isotopic niches of crucian carp and topmouth gudgeon (
Finally, following
Identical statistical and visual approaches as above were used to estimate the overlap between the isotopic niches of crucian carp and topmouth gudgeon at six syntopic sites, out of which the invasive gibel carp was also present at four of them. Statistical analyses and figures were generated using R version 4.2.2 (
On average, the fish in the mesocosm increased in biomass during the experiments, with the highest increase observed in the crucian carp single species and syntopy treatments (248.4 ± 44.5% and 221.0 ± 55.2% of the initial biomass, respectively) and lower growth in the topmouth gudgeon single-species treatments (106.4 ± 25.4% and 107.8 ± 32.2%, respectively). Seven out of 16 treatments with topmouth gudgeon (both single species and syntopy) produced offspring of topmouth gudgeon, while crucian carp produced no offspring. Five topmouth gudgeon died during the experimental period, two in syntopy and three in the single species treatment.
The syntopy treatment in the mesocosm showed that the mean isotopic values of topmouth gudgeon and crucian carp differed significantly (PERMANOVA: F1,38 = 10.56, p < 0.001; Fig.
Layman’s metrics for isotopic niches of crucian carp (CC, Carassius carassius), topmouth gudgeon (TMG, Pseudorasbora parva) and co-existing gibel carp (GC, C. gibelio) in the mesocosm experiment and field evidence. CR = δ13C range, NR = δ15N range, CD = mean distance to centroid, NND = mean nearest neighbour distance, SDNND = standard deviation of nearest neighbour distance.
Mesocosm data | ||||||
Treatment | Species | CR | NR | CD | NND | SDNND |
Single-species | CC | 4.2 | 3.8 | 4.5 | 0.6 | 0.3 |
Syntopy | CC | 2.9 | 7.0 | 1.5 | 0.7 | 0.6 |
Single-species | TMG | 2.4 | 9.4 | 2.3 | 0.7 | 0.4 |
Syntopy | TMG | 2.8 | 3.9 | 1.0 | 0.4 | 0.3 |
Field data | ||||||
Site | Species | CR | NR | CD | NND | SDNND |
Soudkuv lom | CC | 2.60 | 1.80 | 0.87 | 0.40 | 0.13 |
TMG | 1.50 | 1.20 | 0.45 | 0.17 | 0.15 | |
GC | 2.50 | 1.50 | 0.69 | 0.28 | 0.11 | |
Lipi | CC | 4.30 | 2.40 | 1.11 | 0.43 | 0.39 |
TMG | 4.00 | 1.70 | 1.10 | 0.35 | 0.23 | |
GC | 2.40 | 3.90 | 1.16 | 0.49 | 0.32 | |
U Spacku | CC | 4.00 | 2.90 | 1.15 | 0.55 | 0.28 |
TMG | 4.00 | 1.40 | 0.83 | 0.49 | 0.59 | |
GC | 3.30 | 1.90 | 0.97 | 0.40 | 0.25 | |
Vrabce | CC | 2.20 | 1.90 | 0.65 | 0.32 | 0.25 |
TMG | 1.40 | 1.70 | 0.63 | 0.29 | 0.18 | |
GC | 2.20 | 2.50 | 0.71 | 0.40 | 0.48 | |
Smilovice | CC | 1.80 | 2.60 | 0.64 | 0.30 | 0.29 |
TMG | 3.70 | 2.20 | 1.02 | 0.39 | 0.33 | |
Valcha | CC | 2.40 | 1.10 | 0.77 | 0.62 | 0.30 |
TMG | 3.80 | 1.60 | 0.72 | 0.48 | 0.40 |
Mesocosm experiment (A) with treatments of crucian carp (CC, Carassius carassius), syntopy (Syn) and topmouth gudgeon (TMG, Pseudorasbora parva) and field evidence from six sites (B). Isotopic niches shown as 40% standard ellipse area for three species at each site (darker shading) and 95% ellipse area (lighter shading). Solid lines = ellipse areas, symbols = individual values. CC = Crucian carp (Carassius carassius), GC = Gibel carp (C. gibelio); TMG = Topmouth gudgeon (Pseudorasbora parva).
Posterior probability distributions of estimated Bayesian 95% ellipse area for mesocosm experiment and field dataset, modelled with nicheROVER package (
The Bayesian overlap between the isotopic niches of crucian carp and topmouth gudgeon was 56% and 26% at 95% and 40% ellipse area, respectively (Fig.
The study sites differed significantly in their mean δ15N values (PERMANOVA: F5,210 = 232.4, p < 0.001; Fig.
Results of Analysis of Variance (ANOVA) of δ15N and δ13C values at six sites and post-hoc Tukey tests. The Valcha site showed a non-significant permutation analysis of variance for δ15N and δ13C; therefore, no ANOVAs were performed. CC = crucian carp (Carassius carassius), GC = gibel carp (C. gibelio); TMG = topmouth gudgeon (Pseudorasbora parva).
Site | ID | F | df | p | CC×TMG | CC×GC | TMG×GC | |
---|---|---|---|---|---|---|---|---|
Soudkuv lom | 1 | δ15N | 34.9 | 2 | < 0.001 | 0.95 | < 0.001 | < 0.001 |
δ13C | 22.9 | 2 | < 0.001 | 0.85 | < 0.001 | < 0.001 | ||
Lipi | 2 | δ15N | 4.8 | 2 | 0.013 | 1.00 | 0.03 | 0.02 |
δ13C | 30.6 | 2 | < 0.001 | < 0.001 | < 0.001 | 0.01 | ||
U Spacku | 3 | δ15N | 38.7 | 2 | < 0.001 | 0.1 | < 0.001 | < 0.001 |
δ13C | 2.0 | 2 | 0.152 | – | – | – | ||
Vrabce | 4 | δ15N | 4.5 | 2 | 0.018 | 0.172 | < 0.05 | 0.572 |
δ13C | 23.8 | 2 | < 0.001 | < 0.01 | < 0.01 | < 0.001 | ||
Smilovice | 5 | δ15N | 13.8 | 1 | < 0.001 | – | – | – |
δ13C | 28.5 | 1 | < 0.001 | – | – | – | ||
Valcha | 6 | δ15N | – | – | – | – | – | – |
δ13C | – | – | – | – | – | – |
The isotopic niche of crucian carp was generally larger than those of topmouth gudgeon and gibel carp at sites where all three species occurred, while the niche of topmouth gudgeon was larger than that of crucian carp at the Smilovice and Valcha sites, where only these two species occurred (Table
Ellipse areas (EA) of the isotopic niche encompassing 40% or 95% of the data. The ellipse areas and their isotopic niche overlaps were calculated with the R package Siber using Bayesian posterior probability distributions (
Site | ID | Species | 40% EA | 40% EA overlap | 95% EA | 95% EA overlap | P of ind. overlap at 95 % | ||||||
---|---|---|---|---|---|---|---|---|---|---|---|---|---|
CC | TMG | GC | CC | TMG | GC | CC | GC | TMG | |||||
Soudkuv lom | 1 | CC | 1.44 | – | 43.24 | 0 | 8.45 | – | 43.45 | 28.94 | NA | 22.53 | 45.68 |
TMG | 0.40 | – | – | 0 | 2.35 | – | – | 23.90 | 96.56 | 27.28 | NA | ||
GC | 0.86 | – | – | – | 5.09 | – | – | – | 27.61 | NA | 8.89 | ||
Lipi | 2 | CC | 2.22 | – | 12.51 | 0 | 13.04 | – | 48.95 | 10.16 | NA | 3.8 | 41.59 |
TMG | 1.35 | – | – | 9.06 | 7.95 | – | – | 34.56 | 63.80 | 42.5 | NA | ||
GC | 2.26 | – | – | – | 13.27 | – | – | – | 6.73 | NA | 37.24 | ||
U Spacku | 3 | CC | 2.21 | – | 36.49 | 0 | 12.96 | – | 52.56 | 20.04 | NA | 8.57 | 50.61 |
TMG | 1.61 | – | – | 0 | 9.43 | – | – | 31.48 | 73.1 | 20.79 | NA | ||
GC | 1.47 | – | – | – | 8.61 | – | – | – | 14.9 | NA | 20.58 | ||
Vrabce | 4 | CC | 0.99 | – | 6.37 | 0 | 5.81 | – | 55.97 | 42.54 | NA | 40.25 | 49.95 |
TMG | 0.86 | – | – | 0 | 5.05 | – | – | 24.25 | 63.04 | 16.85 | NA | ||
GC | 1.27 | – | – | – | 7.46 | – | – | – | 41.73 | NA | 13.19 | ||
Smilovice | 5 | CC | 0.69 | – | 0.08 | – | 4.09 | – | 42.01 | – | NA | NA | 71.82 |
TMG | 1.80 | – | – | – | 10.55 | – | – | – | 23.30 | NA | NA | ||
GC | – | – | – | – | – | – | – | – | NA | NA | NA | ||
Valcha | 6 | CC | 1.20 | – | 4.94 | – | 7.06 | – | 54.43 | – | NA | NA | 61.16 |
TMG | 1.30 | – | – | – | 7.65 | – | – | – | 37.72 | NA | NA | ||
GC | – | – | – | – | – | – | – | – | NA | NA | NA |
The data collected for this study indicate that the presence of the topmouth gudgeon leads to an overlap of isotopic niches with the native crucian carp, which confirms the first general hypothesis. Both in the field and in the mesocosm, there was evidence of strong overlap between the invasive topmouth gudgeon and the native crucian carp. A very important competitive interaction can occur between young-of-the-year crucian carp and the topmouth gudgeon, as the latter often reaches enormous densities and has a very strong influence on the availability of zooplankton (
At three sites, Lipi, Vrabce and Smilovice, crucian carp was, on average, more enriched in 13C than topmouth gudgeon. This suggests that crucian carp may utilise more aquatic littoral food sources than topmouth gudgeon at these sites, as the base of the aquatic littoral food web generally exhibits higher δ13C values compared to the base of the pelagic food web (
The relative δ15N values are an indicator of the trophic position of the species (
The topmouth gudgeon had a very large isotopic niche in the single species treatment, whereas the syntopic treatment showed a flexible shift in its niche. This was only partly consistent with our third hypothesis and other studies (
Regarding the overlap of isotopic niches, topmouth gudgeon overlapped more with the crucian carp than with the invasive gibel carp and the overlap with topmouth gudgeon was greatest in the presence of the second invasive species, the gibel carp, which confirms the fourth hypothesis. The crucian carp adapts its isotopic niche to the aquatic community (
In this study, we have only addressed part of the interactions between species, namely isotopic niche overlap, which was assessed via a stable isotope approach. However, an organism’s niche consists of many different aspects of resource sharing in space and time (
The topmouth gudgeon has also been found to harm other fish through parasitic behaviour (
Although we tried to study the competition between the two species, topmouth gudgeon and the crucian carp, it turned out that the fish communities very often also contained gibel carp and we were not able to obtain samples from the field where only the two former species co-occur. Gibel carp have very similar habitat requirements to crucian carp, with adaptation to anoxia conditions (
Crucian carp populations have declined considerably in Europe (
The six sites selected for this study were part of the verification of the citizen-science project “Save the Crucian Carp” (
The decline of the crucian carp has led to conservation actions to protect this species in England (
We sincerely thank all those who responded as part of our citizen-science project, which has helped to reveal sites that were involved in the study. We also thank the media for their continuous and comprehensive public coverage. We specifically thank to Vojtěch Kresl and Ondřej Remeš, who identified sites and helped during collection of fish samples.
The authors have declared that no competing interests exist.
The field sampling methods and experimental protocols used in this study were performed by the guidelines and permission from the Ministry of Environment of Czechia (ZN/MZP/2022/630/1428).
The collection of data resulting in this manuscript has been supported by the programme of Regional Cooperation of the Czech Academy of Sciences (R200962201) and the Research Programme Strategy AV21 Water for life for valuable support. Stable isotope measurements were supported by the European Regional Development Fund MEYS CZ.02.1.01/0.0/0.0/16_013/0001782, the Czech MEYS Large Infrastructure for Research MEYS LM2015075, NSF-OCE 1736656 (LIA).
Marek Šmejkal: Conceptualisation, Investigation, Data curation, Formal analysis, Writing – Original draft, Writing – Review & Editing, Visualisation. Kiran Thomas: Data curation, Investigation, Writing – Review & Editing. Zuzana Šmejkalová: Data curation, Writing – Review & Editing. Yevdokiia Stepanyshyna: Investigation, Writing – Review & Editing. Sandip Tapkir: Data curation, Writing – Review & Editing. Travis Meador: Investigation, Formal analysis, Writing – Review & Editing. Mojmír Vašek: Investigation, Formal analysis, Writing – Review & Editing.
Marek Šmejkal https://orcid.org/0000-0002-7887-6411
Kiran Thomas https://orcid.org/0000-0002-4318-1732
Daniel Bartoň https://orcid.org/0000-0001-8042-4564
Sandip Tapkir https://orcid.org/0000-0003-3456-8406
Travis Meador https://orcid.org/0000-0001-8255-5883
Mojmír Vašek https://orcid.org/0000-0001-6386-4015
All of the data that support the findings of this study are available in the main text.