Research Article |
Corresponding author: Lars Pelikan ( lars.pelikan@web.de ) Academic editor: Jaimie T.A. Dick
© 2024 Lars Pelikan, Eglė Šidagytė-Copilas, Andrius Garbaras, Jonas Jourdan, Denis Copilaș-Ciocianu.
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:
Pelikan L, Šidagytė-Copilas E, Garbaras A, Jourdan J, Copilaș-Ciocianu D (2024) Competitive interaction in headwaters: slow upstream migration leads to trophic competition between native and non-native amphipods. NeoBiota 90: 193-216. https://doi.org/10.3897/neobiota.90.112383
|
The spread of non-native species is one of the outcomes of global change, threatening many native communities through predation and competition. Freshwater ecosystems are particularly affected by species turnover with non-native species. One species that has been established in Central Europe for many decades – or even a few centuries – is the amphipod crustacean Gammarus roeselii. Although G. roeselii is nowadays widespread in major river systems, there have been recent reports of its spread into smaller streams that are typically inhabited by the native species Gammarus fossarum. Due to their leaf shredding ability, G. fossarum takes up a key position in headwater streams. This raises the important question, to what extent G. roeselii can equivalently take over this function. To answer this question, we collected both species from nine different sites in a mid-mountain river system (Kinzig catchment, Hesse, Germany) and investigated their functional similarity using a combination of stable isotope analysis, gut content and functional morphology. The species hardly differed in morphological characteristics, only females showed differences in some traits. Gut content analysis indicated a broad dietary overlap, while stable isotopes showed a higher trophic position of G. roeselii. The observed functional overlap could intensify interspecific competition and allow the larger and more predaceous G. roeselii to replace G. fossarum in the future as a headwater keystone species. However, the differentiation in the stable isotopes also shows that co-existence can occur by occupying different trophic niches. Moreover, the wide range of inhabited sites and exploited resources demonstrate the omnivorous lifestyle of G. roeselii, which is likely to help the species succeed in rapidly changing environments.
Freshwater ecosystems, functional morphology, Gammarus, gut content analysis, stable isotope analysis, trophic niche
The introduction and spread of non-native species pose a threat to native communities globally (
A key group that is currently enormously affected by species-turnover is the taxonomic group of amphipods (
Investigating the trophic ecology of invasive species is necessary to gain a better understanding of the community-wide effects of invasions (
The ecological niche of an organism is connected with its functional morphology (
Stable isotope analyses are useful for answering general questions about trophic structure (
Stable isotope analysis is a common method used to reveal the trophic position of organisms in the field (
We hypothesise that, in the focal study area, G. fossarum and G. roeselii share a similar trophic niche. This equivalence is suggested by a laboratory experiment where G. roeselii showed the same leaf consumption rate as G. fossarum (
Gammarus fossarum and G. roeselii were collected with the kick-sampling method on 10 and 11 August 2021 at nine sampling sites in the Kinzig catchment in Hesse, Germany (Table
The nine sampling sites within the Kinzig catchment. The site ID (according to
Site ID | Stream | Species | GPS |
---|---|---|---|
1 | Gründau | G. fossarum + G. roeselii | 50°14.93'N, 9°9.33'E |
5 | Bracht | G. fossarum + G. roeselii | 50°22.62'N, 9°16.22'E |
6 | Bracht | G. fossarum + G. roeselii | 50°26.22'N, 9°16.43'E |
7 | Salz | G. roeselii | 50°25.00'N, 9°21.80'E |
10 | Ulmbach | G. fossarum + G. roeselii | 50°20.18'N, 9°25.70'E |
20 | Schwarzbach | G. fossarum + G. roeselii | 50°21.53'N, 9°33.11'E |
33 | Kinzig headwater | G. fossarum + G. roeselii | 50°18.87'N, 9°36.98'E |
37 | Haselsbach | G. fossarum | 50°13.87'N, 9°22.21'E |
105 | Riedbach, Kinzig | G. fossarum | 50°22.62'N, 9°31.58'E |
The methodology largely followed
Two gammarids within each 2 mm size class were used for the analysis. Amphipods with a body size lower than 7 mm were considered juvenile. In total, 26 juveniles, 18 females and 24 males of G. fossarum and 19 juveniles, 36 females and 28 males of G. roeselii were used. Three leaf replicates per site were used as baseline for trophic position estimates. All samples were sorted, washed with distilled water and dried for 48 h at 60 °C. Afterwards, they were ground to a fine powder with a pestle and mortar. The powder for each amphipod sample was aimed to be 1 mg. In case juveniles did not have a high enough body mass, a composite sample of more individuals of the same species, sex and size class was used. The aimed weight for the leaves was between 3 mg and 4 mg powder for each sample. Carbon and nitrogen stable isotope ratios were measured at the Isotopic Research Laboratory of the Centre for Physical Sciences and Technology in Vilnius, Lithuania. Here, an elemental analyser coupled to the isotope ratio mass spectrometer (EA-IRMS, Flash EA1112–Thermo V Advantage) via the ConFlo III interface was used for the measurement.
In our amphipod samples, the C:N mass ratio varied over 3.5 and in our leaf samples the C exceeded 40%; thus, we corrected the δ13C values for lipid content using the relevant formulae for aquatic animals and plants from
To correct for body size, the morphological measurements were first regressed against the body length and the residuals across all sampling sites were pooled into four species by sex groups (female G. roeselii, male G. roeselii, female G. fossarum, and male G. fossarum). The mean value of residuals was used in case of missing values. The gut content data of all sampling sites were grouped in the same manner. Subsequently, principal component analyses (PCA) were conducted in PAST 4 (version 4.08;
We analysed the difference in trophic position between species by building linear mixed-effects models (LMEMs) for each of the referenced metrics (ΔC and TL). In these models, we considered the interacting fixed effects of size, species and the syntopy (or co-occurrence) factor, while site was included as a random factor. The continuous size variable was centred around the global mean of 8.25 mm for more proper effect testing (but back-transformations were applied for the provided visuals). The effects were tested using type III analysis of variance with Satterthwaite’s approximation for denominator degrees of freedom. These analyses were conducted by employing the R packages lme4 v. 1.1-32 and lmerTest v. 3.1-3 (
In the bivariate stable isotope space, we approximated the population isotopic niches as ellipses containing 95% of the data with their area estimates (BEA95%). We also estimated the overlaps between the species in the six syntopic sites and standardised them as proportions of the sum of the non-overlapping ellipse areas (0 – no overlap, 1 – complete overlap). For this, we used the Bayesian estimation available in R package SIBER v. 1.2.7 (
Results from omnibus PERMANOVA testing for all traits (F = 12.7, p = 0.001) revealed a significant morphological differentiation between sexes within G. roeselii and within G. fossarum and between females of both species, but not males (see Suppl. material
The PCA of the gut content indicated that the most important differentiation between specimens was amongst the detritus, sand vs. plant axis, explaining 29.9% of variation (Fig.
The results of stable isotope analysis revealed pronounced niche differentiation between G. fossarum and G. roeselii, with G. roeselii generally occupying a higher trophic level (Fig.
Referenced stable isotope biplot of studied amphipods showing trophic niches by site (means ± SD). Point size reflects animal body size. Green labels at the means correspond to site IDs from Table
Both LMEMs of ΔC and TL (Table
Effects of amphipod size, species and their syntopic occurrence on isotopic metrics of trophic position within the linear mixed-effects models of a, b referenced δ13C (ΔC) and c, d of trophic level (TL) by a, c non-syntopic vs. b, d syntopic sites. See Table
Results of analysis of variance (type III decomposition) from the linear mixed-effects models of isotopic metrics of trophic position – referenced δ13C (ΔC) and trophic level (TL) – testing for the interacting effects of amphipod size, species (Gammarus fossarum vs. G. roeselii) and their syntopic occurrence. See Fig.
Tested term | df | ΔC model | TL model | ||||
---|---|---|---|---|---|---|---|
df denominator | F | p | df denominator | F | p | ||
Size | 1 | 142.6 | 35.2 | < 0.001 | 142.1 | 43.4 | < 0.001 |
Species | 1 | 10.6 | 125.1 | < 0.001 | 9.4 | 7.0 | 0.026 |
Syntopy | 1 | 9.7 | 114.8 | < 0.001 | 9.1 | 0.6 | 0.459 |
Size : Species | 1 | 142.7 | 6.3 | 0.013 | 142.1 | 5.0 | 0.027 |
Size : Syntopy | 1 | 142.6 | 14.7 | < 0.001 | 142.1 | 0.6 | 0.444 |
Species : Syntopy | 1 | 10.6 | 87.8 | < 0.001 | 9.4 | 1.9 | 0.202 |
Size : Species : Syntopy | 1 | 142.7 | 7.1 | 0.008 | 142.1 | 4.3 | 0.039 |
Regarding the effect of sex across the adult dataset, the stepwise procedure removed all the effects apart from sex from the linear model of ΔC and indicated a model without predictors for TL (although sex was removed last). Thus, we ended up applying simple t-tests using only the sex factor. These indicated a marginally higher female ΔC (t104 = 1.9, p = 0.061), but no effect of sex on TL (t104 = 0.4, p = 0.69).
Isotopic niche widths of the populations and the overlaps between species are provided in Table
Population isotopic niche widths as ellipse areas (BEA95%) and their absolute and relative overlaps by study site. The Bayesian estimates are provided as modes and 95% credible intervals.
Site ID | G. fossarum BEA95% | G. roeselii BEA95% | Overlap | %Overlap |
---|---|---|---|---|
37 | 3.47 (1.56–6.59) | - | - | - |
105 | 5.20 (3.01–9.21) | - | - | - |
6 | 3.51 (2.10–6.79) | 3.05 (1.96–5.37) | 0.02 (0.00–1.44) | 0.00 (0.00–0.20) |
10 | 2.10 (0.99–4.07) | 4.76 (3.09–9.63) | 1.61 (0.46–2.88) | 0.26 (0.08–0.49) |
20 | 1.31 (0.66–2.94) | 1.93 (1.34–3.67) | 0.70 (0.00–1.47) | 0.27 (0.00–0.55) |
33 | 3.53 (2.09–8.38) | 4.85 (3.10–8.40) | 1.73 (0.15–3.17) | 0.24 (0.05–0.39) |
1 | 3.63 (1.56–9.25) | 4.78 (1.79–11.19) | 0.04 (0.00–3.30) | 0.00 (0.00–0.32) |
5 | 3.84 (1.94–10.98) | 2.39 (1.42–4.85) | 1.19 (0.00–2.50) | 0.02 (0.00–0.40) |
7 | - | 6.75 (3.66–12.41) | - | - |
Our study revealed a strong overlap in morphology and gut content between the native G. fossarum and non-native G. roeselii in headwater streams. However, stable isotopes indicated a stronger dietary differentiation between the two species when occurring alone and a more similar trophic niche when occurring together, with G. roeselii, however, generally occupying a higher trophic level. In addition, the gut content analysis confirmed that G. fossarum appears to have a different diet when occurring alone. This indicates that, despite their apparent functional morphological equivalence, the two species exploit different food resources. Below, we expand on the significance of these findings.
Overall, morphological differentiation occurs between the sexes rather than between species when looking at the combined data, pooled over all sampling sites. Only in some traits of the females (gnathopods, pereiopods and coxae) could we see differences between the species. This differentiation between females of both species could be explained by different reproduction characteristics of the species, such as different thermal optima for maximum fecundity (
The gut content analysis showed a strong dietary overlap between the species when pooled over all sampling sites. Our results thus indicate that the foraging on the same food sources might lead to competition between the two species. In principle, a strategy adopted by different amphipods to reduce competition for limited resources could be to utilise resources in different ways, in different microhabitats or at different times (
The stable isotope analysis revealed that both species had more similar trophic niches when they co-occurred. Specifically, G. fossarum underwent a trophic level increase while G. roeselii a decrease which was also accompanied by a shift in the carbon source. Although this result was unexpected and non-intuitive at a first glance, it could be explained by reciprocal predation on juveniles or recently moulted individuals of the other species, as observed in another native and non-native amphipod species pair (
Nevertheless, our stable isotope analysis revealed that, even though both species experience a niche shift, G. roeselii still has generally a higher trophic position compared to G. fossarum when both species occur together, contradicting our first hypothesis of a shared trophic niche. One reason why we have not found this more carnivorous lifestyle in the gut content could be that animal material can be digested more quickly (
Overall, our results indicate that G. roeselii exhibits a broader trophic niche than G. fossarum confirming our second hypothesis. It has been shown that, when G. roeselii co-occurs with G. fossarum, it significantly affects their micro-distribution (
The omnivorous diet of both amphipods, the higher trophic niche of G. roeselii seen in our stable isotope analysis and the observed microhabitat partitioning in the field (
Our study revealed that the non-native G. roeselii is morphologically similar to the native G. fossarum in headwater streams. We also found similar food items in the gut content, which showed the generally omnivorous lifestyle of both species. However, stable isotopes indicated that the trophic niches of both species differ substantially, with G. roeselii being more predaceous and generally having a broader niche. This indicates that, despite their shared morphological characteristics and omnivorous tendencies, there is a noticeable niche differentiation in G. roeselii, consequently influencing the headwater food web. In situations of ample resource availability, co-existence between both species may be possible. However, in cases of resource scarcity, we expect G. roeselii to be competitively superior, particularly given its ability to exploit a broader range of food resources, regardless of their quality.
We thank Erik Aschenbrenner and Jana Kabus for fieldwork assistance. We thank an anonymous reviewer for constructive comments that helped to improve the manuscript. Lars Pelikan thanks the Laboratory of Evolutionary Ecology of Hydrobionts at the Nature Research Centre in Vilnius, Lithuania, for their kind hospitality during his stay.
Lars Pelikan was funded by an Erasmus+ scholarship. The open access publication of this article was funded by the Open Access Publication Fund of Goethe University Frankfurt am Main.
Supplementary information
Data type: docx
Explanation note: table S1. Results (p-values) of the PERMANOVA for the morphological traits of GR-F (female G. roeselii), GR-M (male G. roeselii), GF-F (female G. fossarum) and GF-M (male G. fossarum). PERMANOVA was performed with 9999 permutations and Euclidean similarity index. Bonferroni correction was applied for multiple comparisons between group pairs. Significant p-values (< 0.05) are marked in bold. table S2. Trophic position metrics of different amphipod size groups by study site (derived from stable isotope analysis relative to tree-leaf detritus). Provided values are means with standard deviations.