# Plant diversity drives global patterns of insect invasions

Andrew M. Liebhold, Takehiko Yamanaka, Alain Roques, Sylvie Augustin,
Steven L. Chown, Eckehard G. Brockerhoff, Petr Pyšek

<figure><img src="image_1.png" /></figure>

## To cite this version:

Andrew M. Liebhold, Takehiko Yamanaka, Alain Roques, Sylvie Augustin, Steven L. Chown, et al. Plant diversity drives global patterns of insect invasions. Scientific Reports, 2018, 8, 10.1038/s41598-018-30605-4. doi:10.1038/s41598-018-30605-4

HAL Id: hal-02623059

https://hal.inrae.fr/hal-02623059

Submitted on 26 May 2020

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_2.png" /></figure>
SCIENTIFIC REPORTS

<figure><img src="image_3.png" /></figure>

OPEN

# Plant diversity drives global patterns of insect invasions

Andrew M. Liebhold<fnref n="1" />,², Takehiko Yamanaka³, Alain Roques⁴, Sylvie Augustin⁴, Steven L. Chown⁵, Eckehard G. Brockerhoff⁶ & Petr Pyšek⁷,⁸,⁹

During the last two centuries, thousands of insect species have been transported (largely inadvertently) and established outside of their native ranges worldwide, some with catastrophic ecological and economic impacts. Global variation in numbers of invading species depends on geographic variation in propagule pressure and heterogeneity of environmental resistance to invasions. Elton's diversity-invasibility hypothesis, proposed over sixty years ago, has been widely explored for plants but little is known on how biodiversity affects insect invasions. Here we use species inventories from 44 land areas, ranging from small oceanic islands to entire continents in various world regions, to show that numbers of established insect species are primarily driven by diversity of plants, with both native and non-native plant species richness being the strongest predictor of insect invasions. We find that at large spatial scales, plant diversity directly explains variation in non-native insect species richness among world regions, while geographic factors such as land area, climate and insularity largely affect insect invasions indirectly via their effects on local plant richness.

Received: 12 February 2018

Accepted: 2 August 2018

Published online: 14 August 2018

Insects display an enormous diversity of life histories and exist in virtually every terrestrial environment across all world regions. Given their immense diversity, it comes as no surprise that insects outnumber all other taxa of non-native animals worldwide<fnref n="1" />. Many non-native insect species are notorious, with catastrophic impacts on agriculture, human health and natural ecosystems².

Even though most world biomes are affected by biological invasions, certain regions appear to be more prone to invasions than others. Geographic variation observed in numbers of established non-native species³,⁴ can be attributed both to variation in historical propagule pressure and to habitat characteristics that make certain areas more prone to invasions, i.e. "invasibility"⁵. The concept of biotic resistance describes the effects of community characteristics to promote or inhibit establishment of non-native species.

Both theoretical and experimental studies of invasive plants support the diversity resistance hypothesis, which posits that diverse communities are highly competitive and readily resist invasion⁶. But there is a scale dependency in such relationships; at larger spatial scales (areas > 10 km²) there tends to be an inverse relationship between plant diversity and resistance to plant invasions⁷. Less is known about effects of local diversity on resistance to insect invasions. At small spatial scales, most studies have found that insect abundance and diversity are negatively related to plant diversity though a few have found the opposite⁸⁻¹⁰. At larger spatial scales, geographic variation in numbers of insect invasions can be substantial<fnref n="1" /> but relatively little information exists about how invasibility is related to plant diversity¹¹.

Here we investigate global variation in numbers of invading insect species, seeking explanations for observed patterns. In particular, we explore the role of plant diversity in explaining historical insect invasions. We assemble an extensive data set on numbers of native and non-native insect species in 44 global regions, both mainland and island, and apply structural equation modeling to tease apart the roles of propagule pressure and habitat invasibility. We find that the dominant driver of both native and non-native insect diversity is plant diversity.

<footnote n="1">¹US Forest Service Northern Research Station, Morgantown, WV, 26505, USA. ²Czech University of Life Sciences Prague, Faculty of Forestry and Wood Sciences, Praha 6 - Suchdol, CZ 165 21, Czech Republic. ³Division of Informatics and Inventory, Institute for Agro-Environmental Sciences, NARO, Ibaraki, Japan. ⁴INRA UR0633, Zoologie Forestière, 45075, Orléans, France. ⁵School of Biological Sciences, Monash University, Victoria, 3800, Australia. ⁶Scion (New Zealand Forest Research Institute), Christchurch, 8540, New Zealand. ⁷The Czech Academy of Sciences, Institute of Botany, CZ 25243, Průhonice, Czech Republic. ⁸Department of Ecology, Faculty of Science, Charles University, Viničná 7, CZ 12844, Prague 2, Czech Republic. ⁹Centre for Invasion Biology, Department of Botany & Zoology, Stellenbosch University, Stellenbosch, South Africa. Correspondence and requests for materials should be addressed to A.M.L. (email: aliebhold@fs.fed.us)</footnote>
<figure><img src="image_4.png" /><figcaption>**Figure 1.** Fit of the reduced structural equation model predicting native and non-native plant and insect species richness. Regression parameter estimates are shown next to arrows; black arrows indicate positive estimates, blue arrows indicate negative estimates and weight of each arrow is proportional to the estimated value. Dashed arrows correspond to non-significant relationships. Distance represents insularity and is measured by distance to the mainland.</figcaption></figure>

Since propagule pressure cannot be directly measured, we evaluated the following potential proxies: occupants (i.e., resident human population), gross domestic product (GDP) and distance to mainland (for continental regions this was set at zero). Potential proxies of invasibility considered were latitude, land area, yearly mean temperature, yearly mean humidity, net primary productivity, number of native plant species and number of naturalized alien plant species.

## Results

Pairwise correlations indicate that native and non-native insect species richness are mutually correlated with each other and with most proxies for propagule pressure and invasibility (Supplementary Fig. S1). Given, the high collinearity of these explanatory factors and possibility for complex networks of causality, we applied structural equation modeling (SEM) to investigate these relationships¹².

In the full SEM model (Supplementary Fig. S2), the strongest determinant of non-native insect species richness was native plant richness, though non-native plant richness had a nearly equal effect. Number of human occupants did not have a significant effect on non-native insect richness but both occupants and GDP had strong influences on non-native plant species. Covariances among land area, distance, GDP, occupants and net productivity were generally high which caused an overall low concordance of the model with the data and obscured our ability to differentiate among the effects of these factors.

Following an iterative procedure (see Materials and Methods), a reduced model (Fig. 1) with good concordance with the data was identified. In the reduced model, where occupants, GDP, temperature and productivity are removed, the strongest relationships were the positive effects of land area on non-native and native plant richness. The next strongest relationship was the positive effect of native plant richness on native insects. Again, non-native insect richness was primarily determined by non-native plant richness, followed closely by native plants. Though land area and humidity had strong positive effects on both native and non-native plant richness, they did not have substantial direct effects on either native or non-native insects.

## Discussion

Overall, these results indicate the dominance of plant diversity as a driver of insect invasions. While direct competition plays a key role in plant community assembly¹³, it appears to play a lesser role in insects¹⁴. Herbivory is the dominant life history among insect species, including non-native species¹⁵ and most non-herbivorous species use herbivores as hosts. Furthermore, most insect herbivores are either monophagous or oligophagous¹⁶. Consequently, insect community assembly is typically structured around host plant phylogeny¹⁷. The ability of phytophagous insects to colonize new regions is likely dependent upon locating hosts, and consequently regions supporting more diverse plant communities offer greater opportunities for herbivore colonization.

Land area, latitude, climate, and insularity are all well known drivers of global biodiversity¹⁸,¹⁹. However, our analysis indicates that none of these factors directly affect non-native insect diversity (Fig. 1). These factors more strongly affect plant diversity and thus it appears that their effects on insect diversity mostly operate indirectly via native and non-native plant diversity. Several different mechanisms have been proposed to explain latitudinal gradients of diversity²⁰.
The role of non-native plant diversity as a driver of insect invasions may be related to the phenomenon of “invasive meltdown” in which initial invasions promote subsequent invasions²¹. While invasive meltdown is often considered to occur as a result of invasion-driven disturbance, plant invasions may simply create additional niches for non-native insects to exploit. For example, in Europe 46% of non-native herbivores are limited to non-native host plants²². Many insect herbivores that feed on non-native hosts also utilize native hosts and thus can adversely affect native flora. An important applied implication of these results is that biosecurity efforts aimed at limiting plant invasions may have added benefits in reducing potentially damaging insect invasions as well. Some caution should be used in interpreting the non-native plant/insect invasion association since part of the relationship observed here could have arisen from economic factors (e.g. trade) functioning as drivers of both types of invasions.

Findings reported here of positive effects of plant richness on insect invasions stand in contrast to analyses at smaller spatial scales indicating plant invasions lead to decreased insect diversity<sup>23,24</sup>. Our results also differ from previous findings that insect abundance and diversity tend to be negatively related to richness of all plant species at small spatial scales<sup>9,10</sup>. These previous studies are inconsistent with the positive effect of plant diversity on insect invasions observed here but this inconsistency may be analogous to the “invasion paradox” observed in plant invasions in which there is scale dependency in the effect of diversity on biotic resistance to invasions; at large spatial scales invasibility to plants is positively associated with native plant richness, but at smaller scales invasibility generally decreases with native plant richness<sup>7,25</sup>.

Similar to our report, Hawkins and Porter²⁶ found that native insect herbivore richness was positively associated with plant species richness at large spatial scales. However, they concluded that this did not reflect a direct causal relationship; instead they hypothesized that both plant and insect richness were independently driven by primary productivity and moisture availability. However our analysis indicates the opposite, both native and alien insect diversity are directly determined by plant diversity and the effects of climate and productivity are indirect, operating via plant diversity. Just why our study and theirs have reached contrary conclusions is not yet clear, though differences in the spatial scale of analysis is likely to be a contributing factor, suggesting that investigations of scale effects would be a useful further avenue to explore.

At large spatial scales, plant diversity may increase the number of available niches for insect herbivores and thus promote invasions; this phenomenon is referred to as the “facilitation effect” in the infectious disease literature<sup>27</sup>. But at smaller spatial scales, plant diversity may impair the ability of insects to locate hosts; this is referred to as the “dilution effect”<sup>28</sup> and may cause a negative association between host richness and invasion success. Such opposing effects of facilitation and dilution may explain scale-dependent influences of plant diversity on habitat invasibility to insect invasions.

While the 44 regions analyzed here comprise a diverse sample of land areas from many parts of the world, they are not a random sample and we acknowledge the possibility that unknown bias associated with the areas we have used may affect our conclusions. This analysis utilized all known exhaustive inventories of native and non-native insect species; however, these are most readily available from economically developed countries in the northern hemisphere and selected islands in the southern hemisphere<sup>29</sup>. Future work, especially on scaling effects, would benefit from a broader data set that includes inventories from less well-investigated regions. The development of a global register of introduced and invasive species is a useful step to promote the availability of such information³⁰.

## Materials and Methods

Total numbers of established native and non-native of insect species, as well as numbers of native and non-native naturalized vascular plant species, were derived from a series of species inventories obtained from 44 land areas ranging from small islands to continents (Supplementary Methods online, Supplementary Table S1). For these same regions, we also assembled proxy data on habitat invasibility (land area, latitude, mean annual temperature, mean annual humidity and net primary productivity) and propagule pressure (distance from mainland, occupants and GDP). All data can be viewed in Supplementary Table S2. Numerous pathways facilitate insect invasions so we acknowledge that these proxies do not capture geographical variation in all forms of propagule pressure.

Relationships among propagule pressure proxies, invasibility proxies, native and non-native naturalized plant species richness and native and non-native insect species richness in each of the 44 land areas were quantified using structural equation modeling (SEM). The SEM approach can reveal causal relationships among multiple variables in complex networks (particularly when some of these variables are highly collinear)¹². An *a priori* structural equation model network was identified to represent hypothesized dependencies between predictor and response variables based on logical causal relationships. The full model (Supplementary Fig. 2) included all logical causal connections between invasibility and propagule pressure proxies with native and non-native plant and insect species richness. The model was evaluated using standard SEM procedures¹² implemented with the lavaan package in the R language. All variables (species richness and environmental data) were log-transformed and standardized to stabilize the analyses except latitude was not log-transformed but transformed as absolute values. A maximum likelihood procedure in the lavaan package provides a chi-square statistic that can be used to test the hypothesis of model-data consistency. A value of *P* (chi-square test) below the standard critical value of 0.05 indicates a significant deviation between observed and model implied covariances. Larger *P*-values indicate greater model concordance with data and values of *P* > 0.05 are indicative of overall model adequacy¹²,³¹.

The full model (Supplementary Fig. 2) provided poor concordance with the data (minimum test statistic = 44.518, *df* = 15, *P* = 0.000). Subsequently, all non-significant connections (*P*(>|z|) > 0.5) were removed from the full model and the revised model was evaluated. Next, each connection that was removed in the previous step was re-evaluated using a modification index (*MI*), which represents the reduction of the chi-square value associated with removing the connection. Those connections with *MI* > 3.84 (95% percentile of the chi-square distribution with *df* = 1) were retained and checked again in the next step. These steps were repeated until the overall
P-value of the chi-square test was maximized. In cases where inclusion of a specific connection had a negligible effect on the model P-value, we also compared AIC values among candidate models and selected the model with the smallest AIC. The final reduced model shown in Fig. 1 provided good concordance with the data (minimum test statistic = 1.692, df = 6, P = 0.946).

**Data and materials availability.** The authors declare that all data supporting the findings of this study are available within the article and its Supplementary Information files.

**References**

1. Seebens, H. et al. No saturation in the accumulation of alien species worldwide. *Nat. Commun.* **8**, 14435 (2017).

2. Bradshaw, C. J. et al. Massive yet grossly underestimated global costs of invasive insects. *Nat. Commun.* **7**, 12986 (2016).

3. van Kleunen, M. et al. Global exchange and accumulation of non-native plants. *Nature* **525**, 100–103 (2015).

4. Dawson, W. et al. Global hotspots and correlates of alien species richness across taxonomic groups. *Nat. Ecol. Evol.* **1**, 0186 (2017).

5. Pyšek, P. et al. Disentangling the role of environmental and human pressures on biological invasions across Europe. *Proc. Nat. Acad. Sci.* **107**, 12157–12162 (2010).

6. Kennedy, T. A. et al. Biodiversity as a barrier to ecological invasion. *Nature* **417**, 636–638 (2002).

7. Fridley, J. D. et al. 2007. The invasion paradox: reconciling pattern and process in species invasions. *Ecology* **88**, 3–17 (2007).

8. Andow, D. A. Vegetational diversity and arthropod population response. *Ann. Rev. Entomol.* **36**, 561–586 (1991).

9. Knops, J. M. et al. Effects of plant species richness on invasion dynamics, disease outbreaks, insect abundances and diversity. *Ecol. Lett.* **2**, 286–293 (1999).

10. Jactel, H. & Brockerhoff, E. G. Tree diversity reduces herbivory by forest insects. *Ecol. Lett.* **10**, 835–848 (2007).

11. Liebhold, A. M. et al. A highly aggregated geographical distribution of forest pest invasions in the USA. *Divers. Distrib.* **19**, 1208–1216 (2013).

12. Grace, J. B. *Structural equation modeling and natural systems*. (Cambridge University Press, Cambridge, 2006).

13. Tilman, D. Niche tradeoffs, neutrality, and community structure: a stochastic theory of resource competition, invasion, and community assembly. *Proc. Nat. Acad. Sci.* **101**, 10854–10861 (2004).

14. Kaplan, I. & Denno, R. F. Interspecific interactions in phytophagous insects revisited: a quantitative assessment of competition theory. *Ecol. Lett.* **10**, 977–994 (2007).

15. Roques, A. et al. Temporal and interspecific variation in rates of spread for insect species invading Europe during the last 200 years. *Biol. Invas.* **18**, 907–920 (2016).

16. Forister, M. L. et al. The global distribution of diet breadth in insect herbivores. *Proc. Nat. Acad. Sci.* **112**, 442–447 (2015).

17. Strong, D. R., Lawton, J. H. & Southwood, S. R. *Insects on plants. Community patterns and mechanisms* (Blackwell Scientific Publications. 1984).

18. Gaston, K. J. Regional numbers of insect and plant species. *Functional Ecology* **6**, 243–247 (1992).

19. Field, R. et al. Spatial species-richness gradients across scales: a meta-analysis. *J. Biog.* **36**, 132–147 (2009).

20. Mittelbach, G. C. et al. Evolution and the latitudinal diversity gradient: speciation, extinction and biogeography. *Ecol. Lett.* **10**, 315–331 (2010).

21. Simberloff, D. & Von Holle, B. Positive interactions of nonindigenous species: invasional meltdown? *Biol. Invas.* **1**, 21–32 (1999).

22. Roques, A. Alien forest insects in a warmer world and a globalized economy: Impacts of changes in trade, tourism and climate on forest biosecurity. *New Zeal. J. For* suppl. **40**, 77–94 (2010).

23. Burghardt, K. T. & Tallamy, D. W. Plant origin asymmetrically impacts feeding guilds and life stages driving community structure of herbivorous arthropods. *Divers. Distrib.* **19**, 1553–1565 (2013).

24. Hengstum, T., Hooftman, D. A., Oostermeijer, J. G. B. & Tienderen, P. H. Impact of plant invasions on local arthropod communities: a meta-analysis. *J. Ecol.* **102**, 4–11 (2014).

25. Levine, J. M. & D'Antonio, C. M. Elton revisited: a review of evidence linking diversity and invasibility. *Oikos* **87**, 15–26 (1999).

26. Hawkins, B. A. & Porter, E. E. Does herbivore diversity depend on plant diversity? The case of California butterflies. *Amer. Nat.* **161**, 40–49 (2002).

27. Civitello, D. J. et al. Biodiversity inhibits parasites: broad evidence for the dilution effect. *Proc. Nat. Acad. Sci.* **112**, 8667–8671 (2015).

28. Huang, Z. Y. X., Van Langevelde, F., Estrada-Peña, A., Suzán, G. & De Boer, W. F. The diversity-disease relationship: evidence for and criticisms of the dilution effect. *Parasitology* **143**, 1075–1086 (2016).

29. Chown, S. L., Gremmen, N. J. M. & Gaston, K. J. Ecological biogeography of southern ocean islands: species-area relationships, human impacts, and conservation. *Amer. Nat.* **152**, 562–575 (1998).

30. McGeoch, M. A. et al. Prioritizing species, pathways, and sites to achieve conservation targets for biological invasion. *Biol. Invas.* **18**, 299–314 (2016).

31. Grace, J. B., Anderson, T. M., Olff, H. & Scheiner, S. M. On the specification of structural equation models for ecological systems. *Ecol. Monog.* **80**, 67–87 (2010).

**Acknowledgements**

We thank Laura Blackburn and Eugene Luzader for assistance compiling data. AL was supported by grants from the National Science Foundation Macrosystems Biology Program (grant numbers 1241932, 1638702) and grant EVA4.0, No. CZ.02.1.01/0.0/0.0/16_019/0000803 financed by OP RDE. AR and SA were supported by the European COST projects Alien Challenge (TD1209) and Global Warning (TD1401). PP was supported by long-term research development project RVO 67985939 (The Czech Academy of Sciences) and project no. 14-36079G Centre of Excellence PLADIAS (Czech Science Foundation). SLC was supported by Australian Antarctic Science Program Project 4307. EGB was supported by MBIE core funding to Scion and the Better Border Biosecurity collaboration (www.b3nz.org).

**Author Contributions**

A.M.L., T.Y., A.R., S.A., S.L.C., E.G.B. and P.P. contributed to the conception of the study, assembly of data and drafting of the manuscript. A.M.L. and T.Y. performed statistical analyses.

**Additional Information**

**Supplementary information** accompanies this paper at https://doi.org/10.1038/s41598-018-30605-4.

**Competing Interests:** The authors declare no competing interests.
**Publisher's note:** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

<figure><img src="image_5.png" /></figure>

**Open Access** This article is licensed under a Creative Commons Attribution 4.0 International

License, which permits use, sharing, adaptation, distribution and reproduction in any medium or
format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Cre-
ative Commons license, and indicate if changes were made. The images or other third party material in this
article are included in the article's Creative Commons license, unless indicated otherwise in a credit line to the
material. If material is not included in the article's Creative Commons license and your intended use is not per-
mitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the
copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.

© The Author(s) 2018