Research Article |
Corresponding author: Jonatan Rodríguez ( jonatan@uvigo.es ) Academic editor: Montserrat Vilà
© 2020 Jonatan Rodríguez, Ana Novoa, Adolfo Cordero-Rivera, David M. Richardson, Luís González.
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:
Rodríguez J, Novoa A, Cordero-Rivera A, Richardson DM, González L (2020) Biogeographical comparison of terrestrial invertebrates and trophic feeding guilds in the native and invasive ranges of Carpobrotus edulis. NeoBiota 56: 49-72. https://doi.org/10.3897/neobiota.56.49087
|
Plant invasions impact on biodiversity by altering the composition of native communities by disrupting taxonomic and functional diversity. Non-native plants are often released from their natural enemies, which might result in a reduction of the attack of primary consumers. However, they can also be exposed to the attack of new herbivores that they might not be able to tolerate. Hence, invertebrate communities can be influenced by invasive non-native plants, which in turn modify interactions and change environmental conditions. In this study, we examined the compositional and trophic diversity of invertebrate species, comparing ecosystems with and without the plant species Carpobrotus edulis in coastal areas in its native (South Africa) and introduced (Iberian Peninsula) ranges. Results show that C. edulis has a clear impact on invertebrate communities in its non-native range, reducing their abundance in invaded areas, and particularly affecting certain trophic groups. Invasive C. edulis also alters the invertebrate diversity by not only reducing abundance but also by altering species composition. Overall, the physical dominance of C. edulis modifies the co-occurrence of invertebrate assemblages, reducing the number of trophic groups and leading to substantial effects on primary consumers. Results suggest that the lack of natural enemies might be an important driver of the expansion of C. edulis in its introduced range. Further work is needed to examine long-term changes caused by non-native plants on invertebrate assemblages and the subsequent modification of biological interactions.
Alien species, beta-diversity, biological invasions, enemy release hypothesis, insects, invasion ecology, plant-animal interactions, species richness
Drivers of change such as globalization, habitat fragmentation, and climate change facilitate the arrival, establishment and proliferation of invasive non-native species (
Non-native plants are often released from their natural enemies (i.e. Enemy Release Hypothesis;
Plant invasions are a major concern in Europe (
We investigated the composition and trophic feeding guilds of invertebrate species by comparing coastal areas with and without C. edulis in the native (South African) and introduced (Iberian Peninsula) range of the species. We hypothesized that (i) areas with C. edulis will show higher differences in diversity compared with areas where the species is absent in its non-native range, the Iberian Peninsula. We also hypothesized that (ii) the composition of invertebrate species will be altered through reduced beta-diversity in invaded areas due to the lower replacement of species. We also expected (iii) a higher diversity of trophic groups in well-preserved areas in the native range of C. edulis.
The study was carried out at ten locations, five of them within the native range of C. edulis in the Cape Floristic Region (Western Cape, South Africa) and five in the introduced range in the north-western Iberian Peninsula (Galicia, Spain, and Região Norte, Portugal) (Suppl. material
We established two different sampling areas at each location: plant communities with and without C. edulis, with an area of 200 m2 each (20 m × 10 m) and at least 100 m apart (n = 5). To compare the invertebrate occurrences between areas by visual spot-checks (i.e. hand collection by inspecting each plant), we randomly dropped a standard quadrat grid (0.5 m × 0.5 m) at 3 m intervals ten times within each plot. To identify and quantify the number of aerial invertebrates we used a total of 100 yellow sticky traps (0.2 m × 0.1 m; two faces) in late spring (June in the Iberian Peninsula and December in South Africa). Five traps were evenly spaced within each plot (between 0.2 m to 0.5 m above ground) (see Fig.
Schematic representation of the experimental design for assessing differences in terrestrial invertebrates and trophic feeding guilds in the native (South Africa) and invasive (Iberian Peninsula) ranges of Carpobrotus edulis. Within the plot, standard quadrats (black squares) and yellow sticky traps (yellow rectangles) are represented.
We estimated the invertebrate species richness between the coastal areas of South Africa and the Iberian Peninsula that differ in plant community (with and without C. edulis) by using the specpool function in the ‘vegan’ R package v.2.5–2 (
To evaluate the differences in the invertebrate community, we used Generalized Linear Mixed Models (GLMMs) to test the effect of the region (South Africa and Iberian Peninsula) and plant community (with and without C. edulis) for differences in invertebrate species abundance, richness and diversity indices (Margalef, Shannon, Simpson and Pielou evenness). We compared invertebrate trophic and taxonomic diversity between areas using GLMMs to test the effect of the region and plant community for differences on abundance and species richness. To establish statistical comparisons between models, we used the glmer function in ‘lme4’ R package v.1.1–19 (
To examine beta-diversity differences, we used a multivariate test for homogeneity of dispersion analyses (PERMDISP) of differences in invertebrate composition among studied areas. We assessed the species turnover (replacement of one species by others among locations of the same plant community), nestedness (species richness gain or loss among locations where species lists vary across different plant communities) and total beta-diversity (
Lastly, the species composition was analysed using a non-metric multidimensional scaling (NMDS) ordination based on a Bray-Curtis similarity matrix of standardised and log(x+1) transformed data. Differences were tested using a permutational multivariate analysis of variance (PERMANOVA) with the adonis function (strata = location) in ‘vegan’ R package. PERMANOVA was used to test for differences among coastal areas of South Africa (n = 5) and the Iberian Peninsula (n = 5) for the effect of plant community (with and without C. edulis), distance to the sea, vegetation cover, and their interaction. All statistical analyses were performed using the software programme R (R Development Core Team, 2019; v.3.6.1).
A total of 13,785 invertebrate individuals were identified by visual spot-checks and yellow sticky traps; these were assigned to 90 families and 170 morphospecies (Suppl. material
Relative abundances of taxa (order) within each studied area. The abundance of each taxon was calculated as the percentage of sequences per location for a given invertebrate group. The group ‘Other taxa’ encompasses grouped orders with lower abundance.
Species accumulation curves and the estimated number of species in relation to the species observed indicated that the sampling effort was adequate to capture the majority of the species (Suppl. material
Sample-size-based rarefaction (solid line segment) sampling curves with 95% confidence intervals (shaded areas) for the invertebrate species richness of areas with and without Carpobrotus edulis for both regions, South Africa (SA) and Iberian Peninsula (IP), separated by diversity order: species richness (left panel), Shannon diversity (central panel) and Simpson diversity (right panel).
Number of collected invertebrate species (Sobs) and individuals as well as the estimated richness of species (calculated by Chao2, Jack 1, Jack 2 and bootstrap species estimators) from South Africa and Iberian Peninsula at locations that differ in plant community (with and without Carpobrotus edulis).
Area | Sobs | Individuals | Chao2 (±SD) | Jack1 (±SD) | Jack2 | Bootstrap (±SD) |
---|---|---|---|---|---|---|
Overall (all combined) | 171 | 13785 | 210.86 (17.47) | 207.63 (7.35) | 227.40 | 187.87 (4.06) |
Overall South Africa | 104 | 8283 | 123.64 (11.32) | 124.58 (5.82) | 134.39 | 113.65 (3.29) |
Carpobrotus edulis | 89 | 4684 | 105.28 (9.34) | 109.16 (7.08) | 117.02 | 98.66 (3.95) |
No Carpobrotus edulis | 91 | 3599 | 110.53 (10.73) | 113.08 (6.20) | 122.78 | 101.26 (3.29) |
Overall Iberian Peninsula | 94 | 5502 | 124.11 (15.87) | 119.48 (6.06) | 134.09 | 105.30 (3.23) |
Carpobrotus edulis | 75 | 2283 | 100.13 (13.70) | 98.04 (6.86) | 110.42 | 85.32 (3.66) |
No Carpobrotus edulis | 72 | 3219 | 118.46 (29.51) | 93.12 (6.28) | 108.95 | 80.87 (3.26) |
GLMMs indicated that the levels of abundance and invertebrate diversity indices (species richness, Margalef, Shannon, Simpson, Pielou evenness) tended to be significantly higher in South Africa than in the Iberian Peninsula (Fig.
Effect of the plant community (with and without Carpobrotus edulis) for differences in A abundance B species richness C Margalef D Shannon E Simpson and F Pielou Evenness recorded among regions of South Africa and Iberian Peninsula. Model-adjusted least square means values ± SE are shown. Different letters indicate statistical significance at p ≤ 0.05 level using Generalized Linear Mixed Models.
Results from the Generalized Linear Mixed Models (GLMM) to test the effect of the region (South Africa and Iberian Peninsula) and plant community (with and without Carpobrotus edulis) for differences in abundance, species richness and diversity indices (Margalef, Shannon, Simpson and Pielou evenness) between the invertebrate samples collected. Models were carried with species nested within plant community using Wald Chi-square Method and restricted maximum likelihood (REML). Abundance and species richness were fitted by maximum likelihood (Laplace Approximation) using Poisson distribution and link function = log.
Effect | df | χ2 | Pr(>Chisq) |
---|---|---|---|
Abundance | |||
Region (R) | 1, 100 | 4.717 | 0.029* |
Plant community (PC) | 1, 100 | 0.903 | 0.342 |
R × PC | 1, 100 | 377.912 | 0.001*** |
Species richness | |||
Region (R) | 1, 100 | 18.851 | 0.001*** |
Plant community (PC) | 1, 100 | 4.021 | 0.044* |
R × PC | 1, 100 | 7.661 | 0.005** |
Margalef index | |||
Region (R) | 1, 100 | 14.824 | 0.001*** |
Plant community (PC) | 1, 100 | 11.956 | 0.001*** |
R × PC | 1, 100 | 6.085 | 0.013* |
Shannon index (H’) | |||
Region (R) | 1, 100 | 15.014 | 0.001*** |
Plant community (PC) | 1, 100 | 13.294 | 0.001*** |
R × PC | 1, 100 | 5.587 | 0.018* |
Simpson index (D) | |||
Region (R) | 1, 100 | 7.618 | 0.006** |
Plant community (PC) | 1, 100 | 18.651 | 0.001*** |
R × PC | 1, 100 | 17.924 | 0.001*** |
Pielou evenness (J) | |||
Region (R) | 1, 100 | 3.959 | 0.046* |
Plant community (PC) | 1, 100 | 13.053 | 0.001*** |
R × PC | 1, 100 | 32.060 | 0.001*** |
PERMDISP analyses revealed no dissimilarity for species turnover, nestedness and total beta-diversity, indicating that beta-diversities are very similar among areas with and without C. edulis (Fig.
Non-metric multi-dimensional scaling based on a Bray–Curtis matrix of dissimilarities in A South Africa and in B Iberian Peninsula, showing distances between areas for the species composition and differences in the plant communities (with and without Carpobrotus edulis) and distances to the sea.
Results from the multivariate permutational analysis (PERMANOVA) of differences for South Africa and Iberian Peninsula at locations that differ in plant community (with and without Carpobrotus edulis), distance to the sea and vegetation cover. The PERMANOVA with the adonis function (strata = location) in the ‘vegan’ package in R was based on a Bray-Curtis similarity matrix of standardised on log(x + 1) transformed data.
PERMANOVA | df | SS | MS | Pseudo-F | R2 | Pr(>F) |
---|---|---|---|---|---|---|
(A) South Africa | ||||||
Plant community (PC) | 1 | 0.095 | 0.095 | 0.767 | 0.014 | 0.413 |
Distance to the sea (DS) | 1 | 0.388 | 0.388 | 3.121 | 0.056 | 0.001*** |
Vegetation cover (VC) | 1 | 0.227 | 0.227 | 1.824 | 0.032 | 0.403 |
PC × DS | 1 | 0.418 | 0.418 | 3.365 | 0.060 | 0.001*** |
PC × VC | 1 | 0.140 | 0.140 | 1.124 | 0.020 | 0.446 |
DS × VC | 1 | 0.254 | 0.254 | 2.040 | 0.036 | 0.131 |
PC × DS × VC | 1 | 0.198 | 0.198 | 1.593 | 0.028 | 0.190 |
Residuals | 42 | 5.224 | 0.114 | 0.752 | ||
Total | 49 | 6.945 | 1.000 | |||
(B) Iberian Peninsula | ||||||
Plant community (PC) | 1 | 0.300 | 0.2997 | 1.744 | 0.031 | 0.032* |
Distance to the sea (DS) | 1 | 0.248 | 0.2485 | 1.445 | 0.026 | 0.002** |
Vegetation cover (VC) | 1 | 0.284 | 0.2843 | 1.654 | 0.030 | 0.336 |
PC × DS | 1 | 0.758 | 0.7577 | 4.408 | 0.080 | 0.001*** |
PC × VC | 1 | 0.109 | 0.1096 | 0.638 | 0.011 | 0.864 |
DS × VC | 1 | 0.249 | 0.2491 | 1.449 | 0.026 | 0.343 |
PC × DS × VC | 1 | 0.338 | 0.3386 | 1.970 | 0.035 | 0.027* |
Residuals | 42 | 7.219 | 0.1792 | 0.759 | ||
Total | 49 | 9.507 | 1.000 |
GLMMs revealed significant differences for almost all the groups (except detritivores), while for species richness significant differences between plant communities and regions were only found for herbivores, nectar feeders, parasites and predators (Suppl. material
Within taxonomic groupings, results revealed a significantly higher abundance of almost all groups in areas with and without C. edulis in South Africa (except Diptera), while the areas without C. edulis in the Iberian Peninsula had the lowest values, except for Araneae, Diptera and Other Taxa (Suppl. material
The Cape Floristic Region of South Africa is recognised as a global biodiversity hotspot (
Biodiversity is not always reduced following invasion of non-native plants. Invasive plants can induce the replacement of the native biota by species with similar characteristics (
Plant invasions frequently alter the structure of native communities by disrupting their trophic interactions and changing the ecological processes of the invaded areas, which often creates new environmental scenarios (
We expected a higher diversity of trophic groups in well-preserved areas in the native range of C. edulis. Accordingly, in the Iberian Peninsula, the abundance and species richness of trophic groups were significantly lower than in South Africa. Moreover, our results show that, in the Iberian Peninsula, the abundance and species richness of feeding guilds tended to be significantly reduced by the presence of C. edulis. These results agree with previous research which found that introduced plants can alter the trophic diversity in invaded areas (
At the taxonomic level, our results show differences mainly between regions, but also the invasive C. edulis in the Iberian Peninsula significantly influences the abundance of Diptera, Hemiptera and other taxa groups. In areas invaded by C. edulis, the reduction of secondary consumers as predators could imply a reduction in predation pressure on herbivorous invertebrates (members of Hemiptera and Coleoptera) (
Ecosystem processes and species functional diversity at a local scale can be influenced by the invasion of non-native plants (
Overall, coastal areas in South Africa are more diverse and had a greater abundance of invertebrates compared to the Iberian Peninsula. Invertebrate species composition was affected by C. edulis in the Iberian Peninsula, while the results did not indicate the same in South Africa (see Fig.
Theoretical diagram describing the main results obtained comparing ecosystems with and without the plant species Carpobrotus edulis in coastal areas in its native and introduced ranges. On the left side, the native range is represented (Western Cape, South Africa). On the right, the invasion of C. edulis causes changes in the invaded areas (Galicia, Spain).
This work was funded by Xunta de Galicia, Spain (CITACA Strategic Partnership, Reference: ED431E 2018/07) and carried out within the framework of the project “Retos en la gestión de la planta invasora Carpobrotus edulis. Variabilidad fenotípica y cambios en la relación suelo-planta durante el proceso de invasion” (in Spanish), reference CGL2013-48885-C2-1-R, funded by the Ministry of Economy and Competence (Spanish Government). JR was supported by a research contract (GRC2015/012) from the “Xunta de Galicia/FEDER, Consellería de Educación y Ordenación Universitaria” and a research contract from “Plan de mellora do Centro de Investigacións Agroalimentarias CIA3 do Campus de Ourense, Universidade de Vigo”. JR also acknowledges a short-term research grant awarded by University of Vigo to visit the DSI-NRF Centre of Excellence for Invasion Biology. JR and AN acknowledge funding from the Czech Science Foundation (Project No. 19-13142S and EXPRO grant no. 19-28807X); and the Czech Academy of Sciences (long-term research development project No. RVO 67985939). DMR and AN acknowledge support from the DSI-NRF Centre of Excellence for Invasion Biology and the National Research Foundation of South Africa (grant 85417). AN also acknowledges funding from the South African Department of Environment, Forestry, and Fisheries (DEFF), noting that this publication does not necessarily represent the views or opinions of DEFF or its employees. We thank Javier Puig Ochoa for helpful assistance in the identification of species, and Beatriz Rodríguez-Salvador for helping with data analysis. We are most grateful for the valuable and constructive comments from Subject Editor Montserrat Vilà and the anonymous reviewer that have substantially improved our manuscript.
Biogeographical comparison of terrestrial invertebrates and trophic feeding guilds in the native and invasive ranges of Carpobrotus edulis
Data type: Studied areas, characteristics, occurrences, tables, images