Research Article |
|
Corresponding author: Gabriel Lobos ( galobos@ug.uchile.cl ) Academic editor: Jaclyn Hill
© 2025 Gabriel Lobos, Gianina Tapia, Alejandra Alzamora, Nicolás Rebolledo, Hugo Salinas, Juan Carlos Trujillo, Juan Sánchez, Victoria Gómez-Aburto, Cristóbal Galbán-Malagón.
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:
Lobos G, Tapia G, Alzamora A, Rebolledo N, Salinas H, Trujillo JC, Sánchez J, Gómez-Aburto V, Galbán-Malagón C (2025) Revealing the trophic role of the invasive African clawed frog Xenopus laevis through combined analysis of stable isotopes and heavy metals in a Mediterranean stream from central Chile. Aquatic Invasions 20(4): 495-511. https://doi.org/10.3391/ai.2025.20.4.164197
|
The African clawed frog Xenopus laevis is invasive on four continents, and is recognized as one of the invasive amphibians that generates the greatest impacts in the ecosystems it invades. Although its diet has been studied in its native habitat and invaded areas, its trophic role is still unclear, especially in the communities it invades. We studied the diet of X. laevis, and looked at its stable isotope signatures and its bioaccumulation of heavy metals, to gain a better understanding of its trophic role. The diet was found to consist mainly of aquatic invertebrates, with some consumption of the native fish Cheirodon pisciculus. The isotope analysis revealed that the assimilation of prey by X. laevis is unrelated to the most-consumed item. Xenopus laevis occupied a high trophic position in its own stream and was segregated from fish in by its use of trophic resources. Despite its high trophic position, only biomagnification of copper and zinc was found in relation to some prey, but not manganese or arsenic.
Bioaccumulation, biomagnification, diet, 13C/12C, 15N/14N
Invasive herpetofauna can cause significant trophic disruptions in the ecosystems they invade, and predation is one of their main mechanisms, as it reduces the abundance of native species, causing population extinctions or the change in nutrient cycles (
The African clawed frog Xenopus laevis (Daudin, 1802) has invaded an extensive geographic area on several continents, with established populations in the USA (
Although there are extensive data on the diet of X. laevis in the literature (e.g.,
Heavy metals are typically found with minerals (such as carbonates and sulfates) and organic substances. They accumulate in aquatic sediments, which serves as their primary reservoir (
Xenopus laevis has been categorized as an aquatic top predator due to its large size, generalist diet, and fully aquatic life history (
The study was performed in El Cobre stream, located in the Valparaiso administrative region of Chile and forms part of the Río Aconcagua watershed (33°34.027'S, -70°37.896'W; 623 m elevation). The stream has an intermediate flow, depending on winter precipitation, although some pools persist during the dry season. The high part of the system is drained by the El Gallo and El Sauce streams, in whose watersheds are located “El Soldado” and “Los Navíos” copper mines. The vegetation is dominated by spiny and sclerophyllous shrubs typical of the semi-arid Mediterranean region of Chile (
Ten individuals of X. laevis (all from site CO-1) were collected and used for diet analysis (stomach contents); animals were euthanized using immersion in a 0.2% tricaine methanesulfonate (MS-222) bath. Stomach contents were identified under a stereoscopic microscope to the minimum possible taxonomic level (Fernández and Domínguez 2001). The length and width (mm) of each prey item were measured to estimate the volume and percentage contribution of each prey taxon (
Two samples of M. aquaticum were used to analyze the isotopes of macrophytes. For invertebrates, a pooled compound sample of each taxon was used due to the small size of these organisms after removing hard tissues (mollusk shells, exoskeletons). For fish and amphibians (adults only), a 4 mg sample of muscle tissue (dorsal muscle in fish and thigh muscle in anurans) was obtained. Samples were ground, dried, and lipids were extracted in petroleum ether using a Soxhlet extractor (
Equation 1
R sample and Rstandard are the corresponding ratios of heavy to light isotopes (13C/12C and 15N/14N) in the sample and the standard, respectively. Typical precision of the analyses was ± 0.23‰ for δ15N and 0.25‰ for δ13C.
To perform niche comparisons between X. laevis and the collected fishes, analyses were conducted using R (SIBER package;
The trophic level of organisms was estimated as follows:
Equation 2 (
Where NT is the trophic position n of an organism, δ15N consumer is the value of the organism evaluated, and δ15N base is the nitrogen value of a basal organism (in this case the invertebrate with the lowest position; Baetidae). It was assumed that the NT base in this case occupies position 2 (over the producers of the trophic web). For the trophic discrimination factor (TDF), the reference value 2.3‰ was used (
Cu, Mn, As, and Zn levels were evaluated in sediments, water, and aquatic biota. For the aquatic biota, a variable number of replicates for each species collected was obtained (Appendix
To integrate the information on isotopes and metals, the biomagnification factor (BMF) for X. laevis and the bioaccumulation factor (BAF) for the community (
Equation (3)
Equation (4)
Where CPredator and CPrey are the concentrations of metals in the predator and the prey, respectively, in mg K-1, and the δ15NPredator and δ15NPrey are the ‰ of the 15N content, respectively. Values of BMF greater than 1 are interpreted as a signal of a possible biomagnification process of the metal. In contrast, values lower than 1 are considered as biodilution of the metal in the predator. In the case of BAF, CBiota and CSediment represent the metal concentrations in biota and sediment in mg K-1. BAF values greater than 1 could be interpreted as a signal of bioaccumulation (
Three fish species (Cnesterodon decenmaculatus invasive; Basilichthys microlepidotus and Cheirodon pisciculus native) and one invasive amphibian (Xenopus laevis) were collected from the two sampling sites (Table
Fish and amphibians by sampling site (n = number of individuals). Size (cm), mass (g) are presented as mean ± SD. F indicates number of historical records for the stream between 2018–2022 (15 sampling events).
| Sites | Species | n | Size (cm) | Mass (g) | F |
|---|---|---|---|---|---|
| CO-1-Fishes | Cheirodon pisciculus | 7 | 4.44 ± 0.70 | 1.77 ± 0.59 | 14/15 |
| Cnesterodon decemmaculatus | 53 | 3.18 ± 0.39 | 0.29 ± 0.10 | 14/15 | |
| CO-1-Anurans | Xenopus laevis | 10 | 7.54 ± 0.91 | 41.31 ± 16.94 | 3/15 |
| CO-2-Fishes | Basilichthys microlepidotus | 5 | 8.92 ± 2.17 | 5.38 ± 2.95 | 8/15 |
| Cheirodon pisciculus | 1 | 3.7 | 0.46 | ||
| CO-2-Anurans | Xenopus laevis | 1 | 7 | 38 | |
| Other species seen (but not collected) during the 2018–2022 sampling events | |||||
| Fishes | Gambusia holbrooki (invasive, absent since 2021) | 7/15 | |||
| Trichomycterus areolatus (native) | 5/15 | ||||
| Anurans | Rhinella arunco (terrestrial, tadpoles with short aquatic development) | 10/15 | |||
Density of aquatic invertebrates (individuals/m2 and percentage of total invertebrates) by sampling site (CO-1 and CO-2) collected during Surber net sampling.
| Class | Order | Family | CO-1 | CO-2 | ||
|---|---|---|---|---|---|---|
| Density | % | Density | % | |||
| Gastropoda | Hygrophila | Physidae | 114.81 | 14.90 | 0.00 | 0.00 |
| Clitellata | Oligochaeta | 37.04 | 4.81 | 5.56 | 0.61 | |
| Arachnida | Acariformes | Hydrachnidia | 14.81 | 1.92 | 0.00 | 0.00 |
| Insecta | Coleoptera | Dytiscidae | 25.93 | 3.37 | 44.44 | 4.88 |
| Hydrophilidae | 11.11 | 1.44 | 5.56 | 0.61 | ||
| Gyrinidae | 40.74 | 5.29 | 0.00 | 0.00 | ||
| Diptera | Ceratopogonidae | 3.70 | 0.48 | 0.00 | 0.00 | |
| Chironomidae | 359.26 | 46.63 | 838.89 | 92.07 | ||
| Simuliidae | 3.70 | 0.48 | 0.00 | 0.00 | ||
| Ephemeroptera | Baetidae | 3.70 | 0.48 | 0.00 | 0.00 | |
| Leptophlebiidae | 3.70 | 0.48 | 0.00 | 0.00 | ||
| Hemiptera | Belostomatidae | 0.00 | 0.00 | 5.56 | 0.61 | |
| Notonectidae | 0.00 | 0.00 | 5.56 | 0.61 | ||
| Odonata | Aeschnidae | 29.63 | 3.85 | 5.56 | 0.61 | |
| Coenagrionidae | 118.52 | 15.38 | 0.00 | 0.00 | ||
| Protoneuridae | 3.70 | 0.48 | 0.00 | 0.00 | ||
| Total density (ind/m²) | 770.37 | 911.11 | ||||
| Total number of families | 14 | 7 | ||||
The diet of X. laevis was composed of 41 prey items from seven taxa, including invertebrates (insects and mollusks) and fish (Table
Diet of Xenopus laevis determined via stomach content analysis. Ni = number of prey, %N = percentage of total prey, Vi = volume of prey (mm3), %Vi = percentage of total volume, Oi = number of animals that consumed this prey and %RII = relative importance index, in percentages. The three highest RII values are shown in bold.
| Order | Family | Ni | %N | Vi | %V | Oi | % RII |
|---|---|---|---|---|---|---|---|
| Odonata | Coenagrionidae | 6 | 14.63 | 77.35 | 21.54 | 4 | 19.06 |
| Diptera | Chironomidae | 1 | 2.44 | 0.09 | 0.02 | 1 | 0.32 |
| Coleoptera | Dytiscidae | 5 | 12.20 | 73.11 | 20.36 | 4 | 17.16 |
| Hydrophilidae | 3 | 7.32 | 45.84 | 12.76 | 3 | 7.94 | |
| Ephemeroptera | Baetidae | 3 | 7.32 | 0.27 | 0.07 | 2 | 1.95 |
| Gastropoda | Physidae | 21 | 51.22 | 60.16 | 16.75 | 5 | 44.78 |
| Fish | C. pisciculus | 2 | 4.88 | 102.36 | 28.50 | 2 | 8.80 |
| Total | 41 | 100.00 | 359.17 | 100.00 | 10 | 100.00 | |
The isotope comparison of the community analyzed is shown in Fig.
δ 13C and δ 15N values of food web components in El Cobre stream. SD = standard deviation. Trophic groups are H herbivores, D detritivores and C carnivores. Composite samples of invertebrates were used due to their small body mass (only soft tissues).
| Taxon | Samples | δ15N | δ13C | DS N | DS C | Trophic position | Trophic groups |
|---|---|---|---|---|---|---|---|
| Macrophyte | |||||||
| M. aquaticum | 2 | 2.60 | -31.64 | 0.07 | 3.97 | 0.7 | H |
| Invertebrates | |||||||
| Baetidae | 1 | 2.96 | -28.47 | 0 | 0 | 1.0 | H |
| Physidae | 1 | 4.24 | -26.55 | 0 | 0 | 1.4 | H |
| Dytiscidae | 2 | 5.83 | -29.87 | 1.17 | 1.76 | 2.0 | D |
| Hydrophilidae | 1 | 6.57 | -32.45 | 0 | 0 | 2.4 | D |
| Coenagrionidae | 2 | 8.45 | -29.03 | 2.47 | 3.19 | 3.5 | C |
| Vertebrates | |||||||
| X. laevis | 8 | 11.85 | -23.96 | 1.22 | 1.45 | 4.7 | C |
| C. decenmaculatus | 8 | 12.17 | -29.96 | 0.38 | 0.96 | 4.9 | C |
| C. pisciculus | 8 | 13.08 | -29.37 | 0.79 | 1.55 | 5.25 | C |
| B. microlepidotus | 4 | 13.52 | -16.36 | 0.24 | 1.58 | 5.5 | C |
The nitrogen range in the community was 1.86‰ (Cis: 1.84–1.88), and the carbon range was 13.64‰ (Cis: 13.62–13.65). The trophic diversity, measured as the mean distance to the centroid CD, was 4.87‰ (Cis: 4.86–4.87), and the trophic evenness SDNND was 3.27‰ (Cis: 3.25–3.28). SEAc was different between species evaluated (SEAc C. decenmaculatus = 1.02‰2, SEAc B. microlepidotus = 0.79‰2, SEAc C. pisciculus = 4.49‰2, SEAc X. laevis = 6.10‰2), and position in the δ 13C-δ 15N biplot (Fig.
Stable isotope composition from El Cobre stream. Solid lines enclose the standard ellipse areas (SEAc) containing c. 40% of the data, showing the core isotope niche of each species (3 fish, 1 frog). Dotted lines are the convex hull areas, which are the areas encompassed by each species in the δ 13C-δ 15N plot.
The concentrations of heavy metals varied among the taxa analyzed (Appendix
Concentrations of copper (Cu), manganese (Mn), arsenic (As) and zinc (Zn) (mg/K) in relation to the trophic position. Baet Baetidae (mayfly), Macr macrophyte, Phys Physidae (snail), Dyti Dytiscidae (beetle), Hydr Hydrophylidae (beetle), Coen Coenagrionidae (dragonfly), Xela Xenopus laevis, Cnde Cnesterodon decenmaculatus, Chpi Cheirodon pisciculus, Bami Basilichthys microlepidotus. Green, blue, and red circles represent herbivores, detritivores and carnivores, respectively.
Biomagnification was found for Cu and Zn in X. laevis (Appendix
Invasive species may affect ecosystem functioning by modifying the trophic structure and dynamics (
Xenopus laevis was positioned high in the foodweb (δ 15N = 11.8±1.2), along with the native and invasive fish. However, despite the presence of fish in the stomachs of X. laevis, the frogs occupied a similar trophic position to the fish. The fish present in the stream have been primarily described as consumers of aquatic invertebrates (see
The δ15N values are higher than those reported for other native amphibians such as Eusophus insularis (δ 15N = 6.6±0.8,
The concentrations of heavy metals present in the stream sediments exceeded the international ISQG and PEL standards, aligning with significant contamination reported for some Chilean rivers such as the Aconcagua, attributed to industrial activities, with high levels of Cu, Cd, Zn, Hg, and As (
This study was financed by a fund for freshwater studies of Anglo American. The funders had no role in study design, data collection and analysis, decision to publish or preparation of the manuscript.
GL, CG and JS conceived the idea; GT, AA, JT, NR, HS and GL collected most data; GT, CG and VG organized the dataset and analyses; GL created the figures; GL and CG wrote the first draft.
This project was performed with number permit 448 Chilean Fishing and Aquiculture Subsecretary. All research pertaining to this article did not require any ethics committee approval.
The datasets generated during and/or analyzed during the current study are available from the corresponding author on reasonable request.
We thank to the national fishery and aquaculture service of Chile for the capture permits and anonymous reviewers of this manuscript. CG Acknowledges the funding through ANID/FONDECYT/Regular 1210946 and ANID-PIA-INACH-ACT192057.
Concentration of heavy metals in mg/K (except stream water which is in mg/L) in El Cobre stream (mean ± SD). Samples represent the n used for each analysis. Reference norms corresponded to the Interim Sediment Quality Guidelines (ISQG) and Probable effect levels (PEL).
| Samples | Cu | Mn | As | Zn | |
|---|---|---|---|---|---|
| Sediments | 2 | 3024±530 | 2208±858 | 26±7.6 | 323±57 |
| Stream water | 2 | <0.01 | <0.01 | <0.01 | <0.01 |
| Macrophytes | 2 | 189±66 | 278±169 | 2.13±1.20 | 56.3±9.40 |
| Physidae | 1 | 261.55 | 176.89 | 2.61 | 26.36 |
| Coenagrionidae | 3 | 124±63 | 122±38 | 1.48±1.17 | 99±20.30 |
| Baetidae | 1 | 458.32 | 964.53 | 4.88 | 309.69 |
| Dytiscidae | 1 | 116.47 | 63.74 | 2.06 | 109.87 |
| Hydrophilidae | 1 | 160.54 | 102.49 | 2.50 | 153.70 |
| Cheirodon pisciculus | 2 | 19.8±3.30 | 30.9±5.70 | 1.7±0.20 | 228±112 |
| Xenopus laevis | 5 | 80.9±25.30 | 18.42±4.80 | 0.53±0.20 | 135.9±29 |
| C. decenmaculatus | 1 | 66.60 | 58.30 | 2.06 | 110 |
| B. microlepidotus | 1 | 39.88 | 28.50 | 1.46 | 155.80 |
| ISQG | 35.70 | 500 | 5.90 | 123 | |
| PEL | 197 | 17 | 315 |
Biomagnification factor BMF of heavy metals in Xenopus laevis in relation to its prey. Significant biomagnification is highlighted in bold.
| BMF | Cu | Mn | As | Zn |
|---|---|---|---|---|
| Physidae | 0.11 | 0.04 | 0.07 | 1.84 |
| Coenagrionidae | 0.50 | 0.12 | 0.27 | 1.05 |
| Baetidae | 0.14 | 0.01 | 0.08 | 0.34 |
| Dytiscidae | 0.32 | 0.13 | 0.12 | 0.57 |
| Hydrophilidae | 0.23 | 0.08 | 0.10 | 0.41 |
| Cheirodon pisciculus | 4.50 | 0.66 | 0.34 | 0.66 |
| BAF | Cu | Mn | As | Zn |
|---|---|---|---|---|
| Physidae | 0.09 | 0.08 | 0.10 | 0.08 |
| Coenagrionidae | 0.04 | 0.06 | 0.06 | 0.31 |
| Baetidae | 0.15 | 0.46 | 0.19 | 0.96 |
| Dytiscidae | 0.04 | 0.03 | 0.08 | 0.34 |
| Hydrophilidae | 0.05 | 0.05 | 0.10 | 0.47 |
| Cheirodon pisciculus | 0.01 | 0.01 | 0.07 | 0.70 |
| Xenopus laevis | 0.03 | 0.01 | 0.02 | 0.42 |
| C. decenmaculatus | 0.02 | 0.03 | 0.08 | 0.34 |
| B. microlepidotus | 0.01 | 0.01 | 0.06 | 0.48 |