Research Article |
Corresponding author: Dangen Gu ( gudangen@163.com ) Academic editor: Eric Larson
© 2024 Hui Wei, Lorenzo Vilizzi, Mingsi Zhang, Ying Jiang, Meng Xu, Miao Fang, Fandong Yu, Lu Shu, Xuejie Wang, Dangen Gu.
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:
Wei H, Vilizzi L, Zhang M, Jiang Y, Xu M, Fang M, Yu F, Shu L, Wang X, Gu D (2024) Community-level trophic characteristics and interactions between native and non-native fish: The example of the Lower Pearl River Basin of China. NeoBiota 96: 129-149. https://doi.org/10.3897/neobiota.96.129121
|
Aquatic ecosystems can harbour more than one non-native fish species and this can represent a threat due to trophic interactions with native fishes. However, research on interactions amongst multiple co-occurring native and non-native fish remains scarce. In this study, 551 organisms from 44 native fish, 11 non-native fish, 35 macroinvertebrates (of which one was non-native), together with 162 samples of basal resources were collected from six rivers of the Lower Pearl River Basin of China. Nitrogen and carbon stable isotope analysis was used to calculate community-wide trophic metrics and the degree of trophic overlap between native and non-native fish at both the community and functional feeding group level, together with diet composition. At the community level, there was a high degree of trophic niche overlap between native and non-native fish as a result of similarities in trophic characteristics. At the functional feeding group level, both native and non-native functional feeding groups demonstrated the capacity to occupy the niche space of each other. A significant trophic niche overlap, exceeding 50%, was found between non-native detritivorous and omnivorous fish, suggesting competition. The difference in diet composition between some native and non-native fish depended on the category of diet source across the rivers, suggesting dietary segregation. Albeit limited, the present findings suggest that trophic interaction between native and non-native fish is likely to reach a dynamic equilibrium status in the community owing to trophic segregation of fish species and the antagonistic effects amongst non-native fish.
Diet, impact, multiple invasions, stable isotope analysis, trophic interactions
The management and control of non-native species has become a priority for biodiversity conservation, as invasive species are recognised as one of the major drivers of global environmental change (
The co-existence of species within communities, driven by species interactions, predominantly stems from niche differences (e.g. disparities in environmental requirements), facilitating resource partitioning and reduced interspecific competition (
The Lower Pearl River Basin of China faces a heightened risk of non-native fish invasions due to intensive aquaculture and the ornamental fish trade (
This study employed nitrogen (δ15N) and carbon (δ13C) stable isotope analysis to achieve three objectives: (i) elucidate differences in trophic characteristics between native and non-native fish at the community level; (ii) quantify the degree of niche overlap between native and non-native fish at the community and functional feeding group levels; and (iii) determine whether native and non-native fish have similar diet composition. The outcomes of this study are expected to provide a deeper understanding of the interactions between multiple native and non-native fish, which is crucial for identifying the consequences of multiple species invasions.
Sampling in the Lower Pearl River Basin included the rivers Beijiang, Dongjiang, Liuxihe, Xijiang, Xizhijiang and Zengjiang (Fig.
Map of the six rivers of the Lower Pearl River Basin of China sampled for native and non-native fish, macroinvertebrates and basal resources.
To predict the diet composition of native and non-native fish, potential food sources were collected from the six rivers where fish were sampled. Aquatic insects were captured using a D-shaped net (edge length of 30 cm and pore size of 500 μm:
All samples were dried at 60 °C for 48 h and then ground to powder using a mortar and pestle. For small-bodied macroinvertebrates, two to three individuals were homogenised into one sample. Basal resources including CPOM, FPOM, periphyton, plants and seston were also dried, ground and acidified to remove inorganic carbon. The powdered samples were loaded into individual tin capsules and weighed. Up to three replicates of each animal species and basal resources per river were combusted. The N and C content and isotope ratios were determined using a continuous-flow carrier-gas system (Conflo) equipped with a stable isotope mass spectrometer (Delta V Advantage, Thermo Finnigan, Germany) and an elemental analyser (Thermo Fisher, USA) at Wuhan Botanical Garden, Chinese Academy of Sciences. Replicates of isotopic standard samples (urea) were processed to calibrate for any potential drift (
The δ15N and δ13C were corrected to compare differences amongst rivers. Accordingly, δ15N was converted to trophic position (TP) as per
TP = 2 + (δ15Nfish – δ15Nprey) / 3.4,
where δ15Nfish is the nitrogen stable isotope ratio of each individual fish and δ15Nprey is the average δ15N value of the macroinvertebrate prey resources from each river. In this study, the bivalves of the primary consumers and Angulyagra polyzonata were used in the calculation. Bivalves and A. polyzonata, which were widely distributed across the six rivers, are long-lived filter-feeders and have relatively simple diet sources. Note that the macroinvertebrates were enriched in δ15N due to the pollutions from domestic sewage in the Lower Pearl River (e.g.
δ13Ccorr = (δ13Cfish-δ13CmeanMI) / CRMI,
where δ13Ccorr is the corrected δ13C, δ13Cfish is the δ13C of each individual fish from each river, δ13CmeanMI is the average δ13C of the bivalves and A. polyzonata from each river and CRMI is the δ13C range (δ13Cmax−δ13Cmin) of the bivalves and A. polyzonata (
The effect of fish origin (Origin: native, non-native) and FFG on δ13Ccorr and TP was analysed by linear mixed modelling using the R package lme4. In the models, δ13Ccorr and TP were the response variables, Origin and FFG the fixed effects, River and SL the random effects, as follows:
δ13Ccorr ~ Origin + FFG + Origin:FFG +1 / River + 1/ SL
TP ~ Origin + FFG + Origin:FFG + 1 / River + 1/ SL
Layman metrics were computed for the native and non-native fish in each river using the R package SIBER (
The probability of overlap between the isospace of native and non-native fish was estimated by Bayesian analysis using the R package nicheROVER, which is not sensitive to sample size, with a normal-independent-inverse-Wishart prior to simulate the posterior distribution of the models’ parameters (
The relative proportions of potential sources contributing to the diet of native and non-native fish in each river were analysed using a Bayesian mixing model with R package MixSIAR (
Permutational multivariate analysis of variance was conducted with medium values using the adonis function from the vegan R package to investigate the impacts of River, FFG, Origin and their interaction terms on the potential diet contribution to the fish. A Bray-Curtis dissimilarity measure was employed, with 9999 unrestricted permutations of the raw data and with statistical effects evaluated at α = 0.05. Differences in diet composition amongst fish species were determined by their 95% credible intervals, with overlapping credible intervals indicating no differences amongst fish species (
The TP and δ13Ccorr of native fish were higher, though not statistically different, than that of non-native fish (Fig.
Fixed and random effect coefficients for a linear mixed model describing the effect of Origin (native, non-native), functional feeding group (FFG) and their interaction on trophic position (TP) and corrected δ13C (δ13Ccorr) for fish in six rivers of the Lower Pearl River Basin of China. SL = standard length. See also Fig. S2. Significant results indicated as *** (p < 0.001), ** (p < 0.01), * (p < 0.05) and † (p < 0.1).
Effect | Source | TP | δ13Ccorr |
---|---|---|---|
Random | SL | 0.11 | 0.03 |
River | 0.79 | 0.82 | |
Fixed | Intercept | 1.64 | 0.76 |
Origin (Native) | 0.61 | 0.56 | |
FFG (Herbivore) | -0.33 | 3.24** | |
FFG (Invertivore) | 3.41*** | -0.36 | |
FFG (Omnivore) | 1.52 | 2.10* | |
FFG (Piscivore) | -1.67† | 1.10 | |
FFG (Planktivore) | 0.03 | -1.84† | |
Origin (Native) * FFG (Omnivore) | -0.87 | -0.22 | |
Origin (Native) * FFG (Piscivore) | 2.34* | -0.43 |
Differences between native and non-native fish in trophic position (TP) and corrected δ13C and amongst functional feeding groups (FFG). Significant effects are indicated as *** p < 0.001, ** p < 0.01 and † p < 0.1.
Layman metrics did not differ between native and non-native fish for the δ15N and δ13C ranges, nor for TA and CD (respectively: F = 2.25, p = 0.33; F = 3.53, p = 0.17; F = 55.57, p = 0.18; F = 0.39, p = 0.58). However, the δ15N and δ13C ranges of native fish were wider than those of non-native fish (Fig.
Layman metrics for native and non-native fish for: (A) δ15N range, (B) δ13C range, (C) total area (TA), (D) mean distance to centroid (CD), (E) standard deviation of nearest neighbour distance (SDNND), (F) mean nearest neighbour distance (MNND). Insets show the Layman metrics pooled across the six rivers. Significant effects at p ≤ 0.05 are marked with an asterisk, at p ≤ 0.01 marked with two asterisks.
At the community level, the probability that a non-native fish would be found in the niche of a native fish was higher compared to the opposite (Fig.
Bi-dimensional projections (with 95% CI) of the δ15N and δ13C isotopic niche region for native (red) and non-native (blue) fish.
At the FFG level, the niche space of non-native detritivores and omnivores significantly overlapped with that of the other native FFGs (Fig.
Niche overlap estimates and niche size of native and non-native functional feeding groups (FFG). Only estimates > 50% are shown. (A) Probability of non-native fish FFGs (light red) overlapping with native fish FFGs (light blue) indicated by the direction of the arrows. Line thickness indicates the degree of niche overlap with rivers labelled above the lines. (B) Probability of native fish FFGs overlapping with non-native fish FFGs. (C-H) Niche size of overlapping functional feeding groups in the six rivers. The differences in niche sizes among FFGs can be found in Table S2. Fish FFG abbreviations: Detr_nat = detritivorous native; Detr_nonnat = detritivorous non-native; Herb_nat = herbivorous native; Inve_nat = invertivorous native; Omni_nat = omnivorous native; Omni_nonnat = omnivorous non-native; Pisc_nat = piscivorous native; Plan_nat = planktivorous native. River abbreviations: BJ = Beijiang; DJ = Dongjiang; LXH = Liuxihe; XJ = Xijiang; XZJ = Xizhijiang; ZJ = Zengjiang.
In Dongjiang, non-native detritivores occupied 67.73% of the niche space of native omnivores, while native detritivores and planktivores occupied 51.00% and 72.76% of non-native detritivores. The niche sizes of these non-native fish were similar, except for the niche size of native detritivores which was smaller than that of native omnivores (Fig.
In Xijiang, non-native detritivores occupied 72.57% of the niche space of native omnivores. On the other hand, native planktivores, piscivores and omnivores occupied 75.30% 85.98% and 88.94% of the niche space of non-native detritivores, respectively, while native planktivores, omnivores and piscivores occupied 59.74%, 69.94% and 82.82% of non-native omnivores, respectively. Non-native detritivores occupied 74.35% of the niche space of non-native omnivores, while non-native omnivores occupied 66.31% of non-native detritivores. The niche sizes of these non-native fish were similar, except for the niche size of native planktivores which was smaller than that of non-native detritivores and native omnivores (Fig.
River significantly affected the potential diet composition of fish, but not the interactions with Origin and FFG (Table
Permutational multivariate analysis of variance for the effect of River, functional feeding group (FFG), Origin (native, non-native) and their interaction terms on the diet composition (i.e. median diet proportions) of fish. Significant effects (α = 0.05) in bold.
Source | df | R 2 | F | p |
---|---|---|---|---|
River | 5 | 0.330 | 8.88 | < 0.01 |
FFG | 5 | 0.050 | 1.42 | 0.12 |
Origin | 1 | 0.001 | 0.20 | 0.93 |
River*FFG | 18 | 0.110 | 0.86 | 0.77 |
River*Origin | 5 | 0.050 | 1.37 | 0.14 |
FFG*Origin | 2 | 0.010 | 0.87 | 0.53 |
River*FFG*Origin | 3 | 0.004 | 0.18 | 1.00 |
Residual | 91 | 0.390 |
Although previous studies have indicated that the invasion of non-native fish may lead to a reduction in the trophic position of native fish species as a result of interspecific competition (
No discernible differences in trophic structure and diet composition were observed between native and non-native fish, leading to a high degree of overlap. These findings suggest that both native and non-native fish may face intense competition. Although some native fish exhibited distinct dietary composition from non-native fish across the study rivers, native fish species may confront heightened competition not only from non-native fish, but also amongst themselves. This is attributable to the increased trophic redundancy within the native fish community, with native fish depending on a limited subset of resources (
Interspecific interactions between native and non-native fish have been shown to result mainly in negative or neutral impacts (
The use of FFG level analysis clarifies interactions between native and non-native fish in the Lower Pearl River Basin. Non-native detritivorous fish had a high probability to occupy the niche space of native FFGs. The non-native detritivorous fish C. mrigala, O. niloticus and Pterygoplichthys spp. are the most abundant non-native fish species in the Pearl River Basin (
The niche space of native omnivorous fish was relatively broader than that of non-native detritivorous fish in Zengjiang, which might mitigate the competitive pressure exerted by non-native fish. Native omnivorous, piscivorous and planktivorous fish also had high probabilities of occupying the niche space of non-native detritivorous and omnivorous fish. These results suggest that non-native fish might also face intense competition from native fish. Overall, this study has revealed that non-native fish have established novel trophic interactions with native fish within the community. However, observations also indicated that some native FFGs demonstrated minimal interactions with non-native FFGs, suggesting variability in the degree of ecological integration (see Suppl. material
Co-occurring non-native fish species have the potential to alter community structure through direct or indirect interactions, either facilitating or suppressing one or both invaders (
Understanding the consequence of interactions between multiple non-native fish and native aquatic organisms is a major challenge to manage multiple invasions (e.g.
We thank Dr Wen Zhou and Fujun Gong from Wuhan Botanical Garden, Chinese Academy of Science, for their suggestions to improve data analysis and writing. We also thank Fangcan Chen from Guangzhou Qianjiang Water Ecology Technology Limited Company and all volunteers for their assistance in processing the samples. We also appreciate two anonymous reviewers and the handling editor Dr Eric R. Larson for their comments to improve the quality of this manuscript.
The authors have declared that no competing interests exist.
No ethical statement was reported.
This study was funded by the National Natural Science Foundation of China (Grant No. 32371746 and 31700473), the Central Public-Interest Scientific Institution Basal Research Fund, CAFS (Grant No. 2023TD17 and 2020GH04) and the China Agriculture Research System of MOF and MARA (Grant No. CARS-45).
H. Wei and D. Gu-Conceptualization; Investigation; Data curation; Methodology; Writing-original draft; Funding acquisition. L. Villizi- Data curation; Writing- review & editing. M. Zhang, Y. Jiang, M. Xu, M. Fang, F. Yu, L. Su and X. Wang- Investigation, Writing- review & editing.
Hui Wei https://orcid.org/0000-0001-9641-8214
Lorenzo Vilizzi https://orcid.org/0000-0001-8103-885X
Meng Xu https://orcid.org/0000-0002-2906-6471
Miao Fang https://orcid.org/0000-0002-0924-8323
Dangen Gu https://orcid.org/0000-0002-1310-789X
The datasets and R code generated in this study are available from the first author upon request.
Additional information
Data type: docx
Explanation note: This file includes three supplementary figures and eight tables. The titles of the figures and tables are listed as following: fig. S1. Biplot of δ15N vs δ13C for food resources. fig. S2. Differences in trophic position (left) and corrected δ13C‰ (right) between native and non-native detritivore, omnivore and piscivore fish. fig. S3. Niche overlap estimates of native and non-native functional feeding groups in six rivers of the Lower Pearl River Basin of China. table S1. Fish and macroinvertebrates, with indication of feeding group and origin, sampled from the Lower Pearl River Basin of China. table S2. The mean, standard error and 95% credible intervals of the niche size for native and non-native functional feeding group in six rivers. tables S3–S8. MixSIAR results for the probability of diet contribution for native and non-native fish in Beijiang, Dongjiang, Liuxihe, Xijiang, Xizhijiang and Zengjiang. Species name abbreviations in table S1.