Invasive species usually occupy a wide geographical range in their native area. Invasive species are typically characterized by a number of traits that favor the establishment and spread across new ecosystems, including a broad environmental tolerances, high genetic variability, rapid growth, early sexual maturity combined with a high reproductive rate, short generation time, broad diet, gregariousness, rapid dispersal, and they are often commensal (Ricciardi & Rasmussen, 1998). Of course, not all invasive species meet all these criteria (Lodge, 1993). For example, successful invaders do not necessarily exhibit a broad diet (Vazquez, 2006). Yet a large dietary niche breadth is frequently considered as a hallmark of an invasive taxon.
The dietary niche is a component of the Eltonian niche, defined as the position of an organism, exhibiting a null population growth rate, in the trophic relationships with others organisms of the ecosystem such as its nutrients, predators and competitors (Chase & Leibold, 2003). Another aspect of the ecological niche is the Grinnellian niche, defined as the set of all values of the abiotic parameters enabling the occupancy of an area by a species (Soberón, 2007). The Grinnellian niche has known a recent intensification of research with the development of species distribution models (e.g., Angetter, Lötters & Rödder, 2011; Guisan et al., 2013). Based on a set of occurrence records and predictor variables, these models determine a given species fundamental niche and facilitate the assessment of potential niche shifts when projected onto novel conditions. While these models are based on the assumption that species retain their ancestral traits over time (see Ackerly, 2003; Wiens & Graham, 2005) recent evidence (Broennimann et al., 2007; Mukherjee et al., 2012; Stiels et al., 2015) revealing shifts in realized Grinnellian niches on a macro-ecological scale call this concept into question. How far Eltonian niches operating on a population level are variable, under the assumption of niche conservatism, is less well studied. Whether a population maintains its characteristics, or shifts them in the course of the invasion process, will contribute determining its ecological impact. This information is therefore of crucial importance for conservation practitioners facing the threat posed by invasive species.
Native to southern Africa, the African clawed frog, Xenopus laevis, has been introduced in many countries on four continents, where accidentally and deliberately released individuals have established viable populations (see Measey et al., 2012). Despite its importance as a biological model organism (Cannatella & De Sa, 1993), and the abundance of invasive populations, few field studies have been undertaken in colonized ranges (for a review, see Measey et al., 2012). Xenopus laevis has been reported to negatively affect the invaded ecosystems, and as a consequence has been ranked as having the second greatest impact on native ecosystems by any amphibian (Measey et al., 2016). Shifts in the Grinnellian niche of X. laevis have been recently demonstrated (Rödder et al., in press), whereas studies on changes in the Eltonian niches have not yet been undertaken.
The diet of X. laevis has been studied in the species’ native range of South-Africa (Schoonbee, Prinsloo & Nxiweni, 1992), as well as in several introduced populations in the United States of America (USA) (McCoid & Fritts, 1980), Wales (Measey, 1998a), Chile (Lobos & Measey, 2002), Italy (Faraone et al., 2008), Portugal (Amaral & Rebelo, 2012), and France (Courant et al., 2014). In most studies, the majority of prey items are aquatic, with zooplankton and dipteran larvae being the most frequent. Measey (1998b) also noticed the importance of terrestrial prey. While stomach content analyses conducted in Portugal (Amaral & Rebelo, 2012) and in the USA (McCoid & Fritts, 1980) revealed that X. laevis consumed eggs of fishes and amphibians, no study has reported a direct impact linked to predation. A similar study, conducted in South African aquaculture ponds, revealed that farmed fish larvae constituted a large proportion (5–25% occurrence frequency according to fish size) of frog stomach contents (Schramm, 1987), while another study found X. laevis to consume large quantities of anuran eggs and larvae (Vogt et al., 2017).
Given the wide diversity of prey items, dietary studies on X. laevis usually suggest a generalist feeding behaviour, but only one study has thus far investigated prey electivity (Measey, 1998a). No previous studies have explicitly compared the feeding behavior of populations in different ecological contexts, and with different invasion histories. In this study, we compiled published data on the diet of X. laevis from the USA (McCoid & Fritts, 1980), Wales (Measey, 1998a), Chile (Lobos & Measey, 2002), and South Africa (Vogt et al., 2017) with data collected during recent field work in Portugal and France to test the hypotheses that (i) trophic niche breadth is wider in native populations, hence releaving the capacity of the species to readily adapt to novel environments; and (ii) the diet of invasive populations differs significantly between the invaded ranges depending on local prey availability, and thus resulting in a low degree of electivity and population-specific niche shifts.
Material & Methods
Our dataset comprised 1,458 individuals from six countries, across four continents (Table 1). In most areas (Chile, South Africa, Wales, France), frogs were caught using funnel traps. In Portugal, animals were captured using electrofishing because the colonized habitats, mainly fast flowing streams, prevented the use of traps. In South Africa, research permission was issued by CapeNature (AAA007-01867) and SANParks (RC/2014-2015/001–2009/V1), with ethics clearance from Stellenbosch University Research Ethics Committee: Animal Care & Use (SU-ACUD15-00011). Animals from the Portuguese invasive population were captured under permit no 570/2014/CAPT from Instituto da Conservação da Natureza e das Florestas, in the scope of the “Plano de erradicação de Xenopus laevis nas ribeiras do Concelho de Oeiras” (Eradication plan of Xenopus laevis in the streams of Oeiras Municipality). In France, a research permit was provided by the prefecture of the Deux-Sèvres department.
Stomach content samples were obtained either by stomach flushing or dissection, following euthanasia of individuals by lethal injection of sodium pentobarbital or immersion in MS222. We considered that analyzing and comparing data collected with both dissection and flushing methods were valid and did not induce any bias (Wu, Li & Wang, 2007).
|Period of capture||From 06/2014 To-09/2014||From 05/1995 To 08/1996||From 05/2014 To 10/2014||01/1998 03/2001||From 06/2014 To 08/2014||1975–1976|
|Latitude/Longitude NWb||S34°18′24″ E18°25′35″||N51°27′33″ W3°33′11″||N47°16′14″ W0°33′56″||S33°29′ W70°54′||N38°45′09″ W9°17′27″||See McCoid & Fritts (1980)|
|Latitude/Longitude SEb||S34°20′06″ E19°04′29″||NA||N46°53′41″ W0°31′11″||S33°37′ W70°39′||N38°42′35″ W9°16′25″||See McCoid & Fritts (1980)|
|Capture occasion/site||From 1 to 4||29||3||1||From 1 to 4||1|
|Number of sites||8||1||26||2||12||1|
|Number of individuals||164||375||438||48||352||81|
|Prey collection method||Flushing/ Dissection||Flushing||Dissection||Dissection||Dissection||Dissection|
|Prey frequency in stomachs||Yes||Yes||No||Yes||No||Yes|
Prey items retrieved from the stomach content samples were identified to the lowest taxonomic level possible. However, for analytical consistency, we retained the lowest common taxonomic level that could be identified for all prey items. As volume and mass of prey items were not available for most studies, analyses were performed using prey frequencies. Each taxonomic prey category was assigned to one of the following ecological traits: plankton, benthos, nekton and terrestrial. Some groups of invertebrates belong to different ecological trait categories depending on their life stage, e.g., aquatic in their larval stage and terrestrial in their adult stage (Diptera, Ephemeroptera, Trichoptera). Thus, adults and larvae were treated separately when assigned to different ecological traits even though they belong to the same taxonomic prey category.
The diet of populations was first compared by calculating the relative abundance of each prey category and ecological trait. The relative abundance of prey classes (aquatic invertebrates, terrestrial invertebrates, vertebrates) was also calculated. To assess variation in diet between populations, a Principal Component Analysis (PCA) was performed on reduced and centered relative abundances of each prey category.
The concept of niche breadth can be applied to comparative studies of diet between populations or species (Slatyer, Hirst & Sexton, 2013), even if it has not always been treated within this contextual vocabulary (e.g., Rehage, Barnett & Sih, 2005; Luiselli et al., 2007; Dalpadado & Mowbray, 2013). We calculated niche breadth for all populations using the evenness measure J’. This index is based on the Shannon–Wiener’s index H’ (Shannon & Weaver, 1964), as recommended by Colwell & Futuyma (1971): where pi is the proportion of the prey category i in the diet and n is the number of food categories.
To test whether J′ is affected by the number of study sites the relationship of both variables was assessed using a nonlinear regression.
Prey availability was quantified in habitats for three of the six populations (France, Wales and South Africa) included in this study. Following Measey (1998a), the same sampling method was applied to all populations. Prey selection was assessed using the Vanderploeg & Scavia’s (1979) relativized electivity index recommended by Lechowicz (1982): where ri is the relative abundance of prey category i in the diet and pi is the relative abundance of prey category i in the environment. The number of prey categories included in the analysis is represented by n.
Across all samples zooplankton was the most common prey type with a mean relative abundance of 56.21% (Standard Deviation = 32.80%), followed by ephemeropteran larvae (10.31% ± 23.40%), dipteran larvae (9.68% ± 7.23%), and gastropods (7.24% ± 6.07%). The fifth and sixth most represented prey items were amphibian eggs, excluding X. laevis (4.86% ± 11.08%) and X. laevis eggs (3.46% ± 7.55%) respectively. Three aquatic invertebrate orders (Coleoptera, Diptera, and Heteroptera) were detected in all study sites (Table 2), while most of the terrestrial categories were exclusively found in one or two sites. Cannibalism of larvae and/or eggs was recorded in every locality, except Chile.
|Cat no||South Africa||Wales||France||Chile||Portugal||USA|
|34.80% ± 34.17%|
|7.74% ± 13.87%|
|0.64% ± 0.46%|
|9.45% ± 13.54%|
|40||X. laevis (larvae)||0.04||0.00||0.08||0.00||0.01||0.03|
|41||X. laevis (egg)||0.00||0.41||18.81||0.00||0.00||1.52|
Aquatic invertebrates represented the most consumed prey item class, with a relative abundance ranging from 66% in South Africa to 99% in Wales. Terrestrial invertebrates were rarely consumed and consequently relative abundance ranged from 0.02% in France and the USA to 1.5% in Chile. Variability in relative abundance was highest for vertebrate prey, reaching a maximal relative abundance of 33% in South Africa while being absent in Chile.
According to the PCA (Fig. 1), the diet of the Portuguese and South African populations was respectively characterized by high relative abundances of ephemeropteran larvae, and eggs of native amphibians. The cluster representing Wales and Chile was characterized by a high occurrence of zooplankton and a very low occurrence of all vertebrate prey items. Except for the eggs of X. laevis, French and American samples share the three most abundant prey categories (i.e., Diptera, Gastropoda, and Zooplankton), which explains the short distance between these samples in the PC space.
Benthic taxa represented 8.54% of the prey categories in South Africa, 6.87% in Wales, 14.85% in Chile, 28.8% in the USA, 55.87% in France, and 93.89% in Portugal (Fig. 2). Nektonic taxa represented 39.14% of the prey items in South Africa, 0.14% in Wales, 1.18% in Chile, 1.94% in the USA, 4.20% in France, and 3.01% in Portugal.
Standardized evenness was 0.48 in South Africa, 0.10 in Wales, 0.27 in Chile, 0.42 in the USA, 0.55 in France, 0.41 in Portugal (Fig. 3A). The relationship between evenness (J′) and the number of sampled sites (Ns) followed a logarithmic function (J′ = 0.086∗ln(Ns) + 0.249; Fig. 3B). The slope of the curve decreased from one to ten sites, and stabilized at 25 sites. There was no correlation between trophic niche breadth estimations and the number of stomach contents analyzed per population.
Negative electivity values were observed for most prey items (Fig. 4). In contrast, zooplankton was preferred (positive electivity) in the Welsh, South African and French populations. The electivity values of the other prey categories were variable, indicating that they were either selected, avoided, or not represented in the sampling.
In studies focusing on invasive species, both niche conservatism and niche shifts are commonly reported in the literature (Tillberg et al., 2007; Caut, Angulo & Courchamp, 2008; Comte, Cucherousset & Olden, 2016). Here, we report strong modifications in the realized dietary niche between naturally occurring and introduced or invasive populations of X. laevis. While these differences mainly constituted contractions or expansions of the realized dietary niche, the diet of the Portuguese population represented a shift in the species’ diet. Individuals from this population were captured in fast flowing streams using electrofishing, while all others frogs were caught in lentic environments. The difference might therefore be attributed to the habitat characteristics, reinforcing the hypothesis that X. laevis is a generalist predator that modifies its diet according to available resources.
Our study provides the first analysis of prey availability for a large number of sites that include both natural and invaded areas. Our findings indicate that X. laevis may expand or shift some dimensions of its trophic niche in novel environments. This result is significant, and has clear consequences for evaluating conditions that favor the invasion of newly introduced populations. Suitable trophic conditions allowing a positive population growth rate may be found in many places where prey abundance is sufficient. A recent macroecological assessment revealed that large areas of suitable climatic conditions were available outside the species’ native range, and that X. laevis was likely to expand its range in Europe as a result of climate change (Ihlow et al., 2016). The broad global trophic niche of X. laevis and its ability to adapt its diet according to local conditions, contribute to the strong invasive potential of this species, and the high impact it may induce on its environment (see Kumschick et al., 2017).
The diet of X. laevis has been studied in its native range (Schoonbee, Prinsloo & Nxiweni, 1992; Vogt et al., 2017), and in different invasive populations (McCoid & Fritts, 1980; Measey, 1998a; Lobos & Measey, 2002; Faraone et al., 2008; Amaral & Rebelo, 2012; Courant et al., 2014). The first study carried out in the native range was performed in a fish farm (Schoonbee, Prinsloo & Nxiweni, 1992), which does not necessarily represent the typical diet of native populations. In most studies, including ours, X. laevis was found to predominantly consume relatively small prey items. In Portugal, X. laevis mostly inhabits streams where zooplankton is rare and the main prey items are benthic ephemeropteran larvae. Prey availability in streams was not studied in Portugal, but collecting such data may help understanding the dietary shift we observed in this population.
In all populations but Chile, predation on amphibians was observed, including X. laevis itself, which represented the most frequent vertebrate taxon in the collective sample. This cannibalistic behavior has been recorded in non-native populations, including in the USA (McCoid & Fritts, 1980), Wales (Measey, 1998a), Chile (Lobos & Jaksic, 2005), and Italy (Faraone et al., 2008) but the high frequencies observed in France and South Africa are unprecedented. In the Chilean population, autochthonous amphibians, as well as X. laevis eggs and larvae, were not observed during the sampling period (J Measey, pers. comm., 2016) which may explain their absence in the diet. The small number of stomach contents sampled and sites analyzed for this population may also provide a biased assessment of the putative absence of native amphibians in this area. Predation on amphibians was minor in most populations, except for the native range of X. laevis. This low occurrence of predation on amphibians in other localities may be related to the season during which studies were carried out, and possible changes in the behavior of native amphibians which have co-existed with X. laevis for decades. Our study does not provide any evidence supporting this idea, but ongoing studies should bring new insights into this question. The noteworthy anurophagy reported in this study corroborates the conclusions of Measey et al. (2016) regarding the occurrence of this behavior among pipids.
In this study, terrestrial invertebrates represented the least consumed prey class. However, cumulatively, there were as many taxonomic categories among terrestrial, as among aquatic invertebrates. The occurrence of terrestrial prey items did not vary much between populations, suggesting that there were no local specializations for the capture of terrestrial prey. In previous dietary studies of X. laevis, authors concluded that the high portion of terrestrial prey could not exclusively be explained by the capture of terrestrial invertebrates that had fallen into the water (Measey, 1998b).
From a methodological perspective, sampling effort in each locality was heterogeneous with respect to the sampling period, the number of prospected sites (range: 1–26), and the number of individuals analyzed per population. Consequently, niche breadth could be under-estimated in populations where only few sites were sampled (Wales, Chile, and USA) or where sites were inter-connected (streams in Portugal). According to our results, there is a positive relationship between diet breadth and the number of prospected sites. As no threshold was identified for an optimal number of study sites, we would recommend using as many spatial and temporal replicates as possible for studies aiming at comparing niche breadths. The diet of an individual may be influenced by its size (Schafer et al., 2002), age (Gales, 1982; Rutz, Whittingham & Newton, 2006), and sex independently of size (Gales, 1982; Göçmen et al., 2011; Van Ngo, Lee & Ngo, 2014). In our study, these data were not available for every population, preventing us from analyzing their effect on diet breadth. In other respects, the unique native population included in our study may not be representative of the diet in the native range. Estimations of niche breadth and prey selectivity may have differed had they been based on a larger sample of native populations. Therefore, we encourage future investigators to consider as many naturally occurring populations as invasive populations in studies aimed at understanding the feeding ecology of invasive animals.
Our results indicate that no prey categories are strongly selected, except for zooplankton in Wales. This suggests that X. laevis does not usually specializes its diet and hence does not develop a population specific dietary niche. This characteristic may enhance its capacity to establish and spread in novel environments. Potential perturbations of X. laevis in its environment, linked to the large predation on small prey items, are crucial elements of the trophic relationships in aquatic ecosystems that still need to be demonstrated. Our results reflect the diet of X. laevis in invasive populations after decades of colonization and do not necessarily reflect the diet of the species at the moment of introduction. Some invasive species modify their diet during the years, or decades following habitat colonization (e.g., Tillberg et al., 2007; Gkenas et al., 2016). Comparing the diet of individuals at the core and the edge of a newly colonized area may be an effective approach to investigate the change in dietary composition of X. laevis during the colonization process.
Individual diet data for the Chilean, French, Portuguese and South African populations of Xenopus laevis
Calculation of the electivity index for the French, South African and Welsh populations of Xenopus laevis
Electivity values were calculated following the Vanderploeg and Scavia’s formula (1979) for each prey category detected in the diet of Xenopus laevis and in its environment.