# The impact of tree diversity on different aspects of insect herbivory along a global temperature gradient - A meta-analysis.

Stephan Kambach, Ingolf Kuhn, Bastien Castagneyrol, Helge Bruelheide

► To cite this version:

Stephan Kambach, Ingolf Kuhn, Bastien Castagneyrol, Helge Bruelheide. The impact of tree diversity on different aspects of insect herbivory along a global temperature gradient - A meta-analysis. PLoS ONE, 2016, 11 (11), 14 p. 10.1371/journal.pone.0165815. [https://doi.org/10.1371/journal.pone.0165815](https://doi.org/10.1371/journal.pone.0165815).

HAL Id: hal-02630854

https://hal.inrae.fr/hal-02630854v1

Submitted on 27 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_1.png" /></figure>

Distributed under a Creative Commons CC BY 4.0 - Attribution - International License
<figure><img src="image_2.png" /></figure>

CrossMark

click for updates

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

OPEN ACCESS

**Citation:** Kambach S, Kühn I, Castagneyrol B, Bruelheide H (2016) The Impact of Tree Diversity on Different Aspects of Insect Herbivory along a Global Temperature Gradient - A Meta-Analysis. PLoS ONE 11(11): e0165815. doi:10.1371/journal.pone.0165815

**Editor:** Hideyuki Doi, University of Hyogo, JAPAN

**Received:** June 17, 2016

**Accepted:** October 18, 2016

**Published:** November 11, 2016

**Copyright:** © 2016 Kambach et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

**Data Availability Statement:** All relevant data are within the paper and its Supporting Information files.

**Funding:** The work of SK was funded by the Helmholtz Centre for Environmental Research - UFZ, the Helmholtz Research School for Ecosystem Services under Changing Land-use and Climate (ESCALATE) and the Martin Luther University Halle-Wittenberg. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

RESEARCH ARTICLE

# The Impact of Tree Diversity on Different Aspects of Insect Herbivory along a Global Temperature Gradient - A Meta-Analysis

Stephan Kambach¹,²*, Ingolf Kühn¹,²,³, Bastien Castagneyrol⁴, Helge Bruelheide²,³

1 Department of Community Ecology, Helmholtz-Centre for Environmental Research, Halle (Saale), Germany, 2 Institute of Biology/Geobotany and Botanical Garden, Martin Luther University Halle-Wittenberg, Halle (Saale), Germany, 3 German Centre for Integrative Biodiversity Research (iDiv) Halle-Jena-Leipzig, Leipzig, Germany, 4 BIOGECO, INRA, Univ. Bordeaux, 33610 Cestas, France

* stephan.kambach@ufz.de

## Abstract

Forests with higher tree diversity are often assumed to be more resistant to insect herbivores but whether this effect depends on climatic conditions is so far poorly understood. In particular, a forest's resistance to herbivory may depend on mean annual temperature (MAT) as a key driver of plant and insect phenology. We carried out a global meta-analysis on regression coefficients between tree diversity and four aspects of insect herbivory, namely herbivore damage, abundance, incidence rate and species richness. To test for a potential shift of tree diversity effects along a global gradient of MAT we applied mixed-effects models and estimated grand mean effect sizes and the influence of MAT, experimental vs. observational studies and herbivores diet breadth. There was no overall effect of tree diversity on the pooled effect sizes of insect herbivore damage, abundance and incidence rate. However, when analysed separately, we found positive grand mean effect sizes for herbivore abundance and species richness. For herbivore damage and incidence rate we found a significant but opposing shift along a gradient of MAT indicating that with increasing MAT diversity effects on herbivore damage tend towards associational resistance whereas diversity effects on incidence rates tend towards associational susceptibility. Our results contradict previous meta-analyses reporting overall associational resistance to insect herbivores in mixed forests. Instead, we report that tree diversity effects on insect herbivores can follow a biogeographic pattern calling for further in-depth studies in this field.

## Introduction

Insect herbivores can compromise the functioning of forest ecosystems [1]. Insect herbivory is controlled by top-down mechanisms involving natural enemies of herbivores [2], bottom-up mechanisms including tree defences [3] and associational effects provided by tree diversity [4, 5]. So-called associational effects occur if herbivory on individual trees is influenced by the identity and density of neighbouring trees [6] and are assumed to be key regulators of herbivory [6, 7]. An increase in tree diversity is often reported to decrease herbivore pressure
**Competing Interests:** The authors have declared that no competing interests exist.

(associational resistance-AR; [7-9]), but examples for an opposite relationship can also be found (associational susceptibility-AS; [7, 10]).

Associational resistance, on the one hand, occurs if a higher diversity of tree species reduces the availability of host trees for specialist herbivores because of reduced resource availability and reduced encounter rates of herbivores and hosts (resource dilution effects; [11, 12]). Mixtures of host and non-host tree species can also reduce the success of herbivores to detect suitable resources, both visually [13] and olfactorily [14].

Associational susceptibility, on the other hand, occurs if a population of herbivores builds up on a more preferred host species and then 'spills over' to admixed less palatable species [10, 15]. Generalist herbivores might furthermore benefit from a more diversified nutrition [16, 17] as suggested by the dietary mixing hypothesis [18, 19].

When accounting for the whole community of herbivores that are able to attack a particular tree, any change in the proportion of generalist vs. specialist herbivores may then shift the balance between AS and AR. Yet, the proportion of specialized insect herbivores was reported to increases towards lower latitudes [20]. As such, it is likely that the direction of associational effects changes along latitudinal gradients.

Several other key features of plant-herbivore interactions are known to vary with latitude. Numerous studies also found that global gradients exist for plant species richness [21, 22], herbivore species richness [23-25], plant defences [26, 27], herbivore pressure [25, 28, 29], leaf herbivory [30] and trophic interactions [31]. Since temperature is a key driver of herbivore development and abundance [32], consumption rates [33] and host-plant choices [34] it is likely to influence these plant-herbivore interactions [33, 34]. Yet, and surprisingly, it still remains to be tested whether the direction and strength of associational effects also change along a global gradient of mean annual temperature (MAT).

In addition to the overlooked effect of MAT on associational effects, the current understanding of AR and AS suffers from several methodological and conceptual biases, including a lack of considering diversity gradients, and confusion between functional (i.e. consumption or proportion of attacked tissue) and quantitative (i.e. abundance and species richness) responses of herbivores.

Previous meta-analyses on associational effects in tree stands mainly compared single species (i.e. monocultures) with mixed stands, irrespective of species richness and species evenness [8, 9]. Yet, the relative share of tree species within tree stands may critically change the way herbivores perceive stand quality [8, 35]. To our knowledge, no study so far has investigated whether the strength and direction of associational effects depend on the metric of tree diversity applied. We therefore compiled studies that reported on either the Shannon or Simpson diversity of tree stands reflecting the relative contribution of host concentration and relative frequency (see [6] and [36] for a discussion about host concentration vs. frequency). If species proportions were not available we applied a gradient of tree species richness or, if species richness was not reported, conducted a comparison between single species and mixed tree stands.

Beside the simplification to monoculture-mixture comparisons, the meta-analyses mentioned above [8, 9] also simplified the response of herbivores by pooling studies that reported on the actual damage inflicted and studies that reported on the abundance of insect herbivores, assuming that more herbivores always cause more damage. This assumption mainly holds for some herbivores such as bark-beetles, leaf miners or galls, but is more controversial for most of defoliators. For instance, herbivore abundance and the associated damage are not necessarily correlated [37, 38] and these two aspects have been reported to sometimes respond differently to plant diversity [7]. Beside, how the species richness of insect herbivores depends on tree diversity has not been summarized in previous meta-analyses. Thus, we separately
analysed the response of insect herbivores to increasing tree diversity into the response of her-
bivore damage, abundance, incidence rate (i.e. the proportion of attacked plant tissue) and the
species richness of herbivores.

In the present study, we went beyond previous meta-analyses by testing whether the direc-
tion and strength of associational effects on forest trees can partly be attributed to the MAT at
the study location. Moreover, we asked whether the focus on different aspects of herbivory (i.e.
the amount of insect herbivore damage, abundance, incidence rate and species richness) as
well as on different metrics of tree diversity (i.e. the Shannon/Simpson diversity, the species
richness and the comparison of single species vs. mixed tree stands) can lead to diverging asso-
ciational effects.

Specifically, we hypothesized that

1. insect herbivory (pooled over herbivore damage, abundance and incidence rate) is negatively related to tree species diversity;

2. insect herbivore damage, abundance, incidence rate and species richness differ in their relationship to tree species diversity;

3. the relationships between the four aspects of insect herbivory (i.e. herbivore damage, abundance, incidence rate and species richness) and tree diversity change with mean annual temperature.

**Material and Methods**

**2.1. Search strategy and inclusion criteria**

We performed a scoping search with combinations of relevant search terms including the fol-
lowing key words: forest, tree, diversity, richness, herbivores, pest, damage, monoculture, mix-
ture and plantation. The following literature databases were queried in August 2016: Thomson
Reuters Web of Science, Google Scholar and Cab Direct. These search queries were partly con-
ducted with the help of the program *Publish or Perish 4* [39]. The PRISMA checklist of items to
include when reporting a systematic review is shown in S1 PRISMA checklist [40] and the
search protocol used for each database is fully described in S1 Text. After additionally review-
ing the literature cited in relevant articles, our search initially yielded 3,707 articles. The search
and subsequent selection process is depicted in the flow chart Fig 1.

After the initial screening, we retained 1,234 articles for further in-depth examination. To
be included in the present meta-analysis, an article had to report on i) either the Shannon
diversity, Simpson diversity, species richness or a monoculture-mixture comparison of tree
stands and ii) either the amount of herbivory sustained by trees, the abundance or incidence
rate (e.g. the proportion of individuals or parts that were attacked) or species richness of insect
herbivores. All causes leading to dismissing articles are fully described in S1 Text.

**2.2. Effect sizes**

Due to the heterogeneity of studies, the compiled measurements of herbivory vastly differed in
dimensions and order of magnitude. In the case of herbivore damage, studies reported the per-
centage of leaf area removed or percentage of crown volume damaged. In the case of herbivore
abundance studies reported on the number or density of individuals, galls and egg clusters.
Incidences rate were determined as percentage of trees, leaves, or branches attacked or with
herbivores present. Species richness was measured as the number of insect herbivore species
captured on trees or in traps.
<figure><img src="image_4.png" /><figcaption>**Fig 1. Flow chart of the study search and selection process.**</figcaption></figure>

doi:10.1371/journal.pone.0165815.g001

Our main objective was to test the correlation between an increase in tree diversity and the four aspects of tree herbivory (i.e. damage, abundance, incidence rate and species richness). Since Pearson's correlation coefficient *r* becomes skewed as it approaches ±1 we transformed *r* to Fisher's z-scores which range from -∞ to +∞, have the same sign as *r* values and are commonly applied in meta-analyses [41–43]. Uncertainty for each effect size was estimated by calculating the corresponding variance estimate (v = 1/(n-3), where n is the sample size [41].

When an article reported multiple measurements (e.g. for different taxa or different aspects of herbivory) we retained them as separate study cases within the same forest and accounted for this non-independence in the statistical analyses (see below).

When studies did not report any *r* value but provided tables or figures with information on tree diversity and any aspect of herbivory we extracted the raw data using the software *ImageJ* [44] and re-calculated the corresponding *r* values. This could only be achieved when there were at least four records for both tree diversity and insect response.
We extracted the mean annual temperature (MAT) for the approximate study locations from the *WorldClim*-global climate data [45]. In order to test for potentially confounding factors, we furthermore distinguished between experimental (i.e. plantations with purposefully manipulated tree diversity) and observational studies (in either semi-natural forests or forests that had not been planted to test for any diversity effects) and classified the diet breadth of herbivores as being either specialists that feed on a single species or genus of tree species or generalists which can utilize a wide range of tree species.

## 2.3. Statistical Analyses

We first tested the grand mean effect size corresponding to the overall effect of increasing tree diversity on the pooled effect sizes of herbivore damage, abundance and incidence rate. Using a random-effect meta-analysis model [43] with restricted maximum likelihood we calculated the model intercept, (i.e. averaged Fisher's z-scores) and the corresponding bootstrap confidence interval (CI). Effect sizes were weighted based on the inverse of their variance estimate. Since multiple effect sizes from the same forest location cannot be considered as fully independent we incorporated a hierarchical random-effect structure with the single effect sizes nested within forest locations. The grand mean effect size was considered statistically significant if the 95% CIs did not include zero. We also estimated the amount of residual heterogeneity ($\tau^2$; [43]) and tested whether effect sizes displayed significant between-study heterogeneity by applying the weighted Cochran's Q-test.

To test whether effect sizes depended on the aspect of herbivory we pooled all effect sizes on herbivore damage, abundance, incidence rate and species richness of herbivores, included this aspect of herbivory as a moderator in the previous random-effect model and calculated the test statistic for the omnibus test of model coefficients ($Q_M$). We then split the dataset based on the aspect of herbivory, applied the intercept only-model from above and calculated for each aspect separately the grand mean effect size, $\tau^2$ and Q-test statistic.

For each aspect of herbivory, we then tested the correlation between effect sizes and MAT, the type of study (experimental vs. observational) and the diet breadth of insect herbivores.

We thereby included these variables as moderators in the previous random-effect models and tested the significance of the obtained parameter estimates against a normal distribution. To exclude statistically insignificant moderators and to yield robust, most parsimonious models, we performed a backward selection based on error probabilities ($\alpha = 0.05$) calculated with maximum likelihood. The parameter estimates of the resulting four minimal adequate models were then calculated with restricted maximum likelihood. For each model we separately calculated $\tau^2$, $Q_M$ and the test statistic for the amount of residual (i.e. unexplained) heterogeneity ($Q_E$, weighted Cochran's Q-test). The amount of variation in effect sizes that predictors accounted for was finally calculated as $Q_M / (Q_M + Q_E)$.

To check for publication bias resulting from the omission of extreme or unlikely results, we visually inspected funnel plots and tested for their asymmetry by applying Egger's regression test [42, 46]. We furthermore tested whether the year of publication correlated with effect sizes or MAT at the study location by calculating Pearson correlation coefficients and the corresponding test statistics. We finally checked whether effect sizes correlated with the reported metric of tree diversity (i.e. Shannon/Simpson diversity, species richness and single vs. mixed species stands) by including this moderator in the previous random-effect model and calculating $Q_M$.

All statistical analyses were conducted in R [47] using the packages *metafor* for meta-analyses [48], *vegan* for calculating the Shannon diversity [49], *raster* for extracting the *WorldClim* data [50] and *ggplot2* for graphical representations [51].
## Results

The final dataset consisted of 60 studies with 173 study cases that reported on the correlation between tree diversity and insect herbivore damage (53 study cases), abundance (52 study cases), incidence rate (40 study cases) and species richness (28 study cases), respectively. Regarding the metric of tree diversity, 94 and 7 study cases reported the Shannon or Simpson diversity of tree species, respectively, whereas 44 study cases reported the richness of tree species. In the remaining 28 study cases, we could apply a comparison between tree monocultures and mixtures of tree species, without any further quantification of tree species diversity.

In 92 study cases, herbivory was measured on a focal tree species whereas in the remaining 81 cases herbivory and especially herbivore abundance and species richness were measured at the plot level, either by summarizing herbivory over all tree species or by reporting capture rates in insect traps that could not be related to certain tree species. All study cases and study characteristics are documented in S1 File.

The included study cases covered a latitudinal range of -36.7° to 62.8° and spanned a gradient of -3.3°C to 26.9°C mean annual temperature (MAT). Study sites were most frequent in Europe and North to Middle America and sparse in Asia, Africa and Australia (Fig 2).

Pooled over all study cases that reported herbivore damage, abundance and incidence rate, the grand mean effect size was not significant and the corresponding funnel plot was symmetric (Table 1, S2 Text, Egger's test: p = 0.58). There was a significant amount of residual heterogeneity indicating that heterogeneity in true effect-sizes could be accounted for by moderators.

Across the whole dataset (including effect sizes on herbivore species richness) the aspect of herbivory was a significant moderator of effect sizes as indicated by the significant omnibus test of model coefficients ($Q_M = 10, df = 3, p = 0.02$). When analysed separately, study cases on herbivore abundance and species richness yielded significantly positive grand mean effect sizes whereas herbivore damage and incidence rate showed no significant relationship with increasing tree diversity (Table 1, Fig 3). Funnels plots for these separate models were symmetric following Egger's regression test (damage: p = 0.14, abundance: p = 0.7, species richness: p = 0.25) except for herbivore incidence rate that displayed a positive correlation between effect sizes and the corresponding variance estimates (estimated effect sizes: $2.1 \pm 0.9$ standard error, p = 0.02).

<figure><img src="image_5.png" /><figcaption>Fig 2. Locations of the studies included in this meta-analysis reporting on the relationship between the four aspects of insect herbivory and tree Shannon/Simpson diversity (circles), tree species richness (squares) or the comparison between single vs. mixed stands (triangles). The colouring indicates the different aspects of insect herbivory. Made with Natural Earth. Free vector and raster map data (naturalearthdata.com).</figcaption></figure>

doi:10.1371/journal.pone.0165815.g002
<table><caption>Table 1. Intercept-only and most parsimonious mixed-effects meta-regression models.</caption><thead><tr><th>Model</th><th colspan="3">Model statistics</th><th>Parameters</th><th>Mean</th><th>Se</th><th>z-value</th><th>p-value</th></tr></thead><tbody><tr><td colspan="9"><strong>Grand mean all</strong></td></tr><tr><td>(= damage, abundance and incidence rate pooled)</td><td>τ<sup>2</sup> = 0.11<br>Q = 395.35</td><td>df = 144</td><td>p &lt; 0.001</td><td>Intercept</td><td>0.01</td><td>0.03</td><td>0.31</td><td>0.76</td></tr><tr><td colspan="9"><strong>Damage</strong></td></tr><tr><td>Grand mean</td><td>τ<sup>2</sup> = 0.07<br>Q = 92.24</td><td>AICc = 67.03<br>df = 52</td><td>p &lt; 0.001</td><td>Intercept</td><td>0.06</td><td>0.06</td><td>1.12</td><td>0.26</td></tr><tr><td>Most parsimonious model</td><td>τ<sup>2</sup> = 0.07<br>Q<sub>E</sub> = 77.99<br>Q<sub>M</sub> = 8.22</td><td>AICc = 65.63<br>df = 50</td><td>p = 0.01<br>p = 0.02</td><td>Intercept<br>(observational at zero MAT)</td><td>0.14</td><td>0.1</td><td>1.33</td><td>0.18<br><strong>MAT</strong><br><strong>Experimental vs. observational</strong></td></tr><tr><td colspan="9"><strong>Abundance</strong></td></tr><tr><td>Grand mean</td><td>τ<sup>2</sup> = 0.19<br>Q = 220.33</td><td>AICc = 93<br>df = 51</td><td>p &lt; 0.001</td><td>Intercept</td><td>0.08</td><td>0.04</td><td>2.13</td><td>0.03</td></tr><tr><td>Most parsimonious model</td><td>τ<sup>2</sup> = 0.15<br>Q<sub>E</sub> = 210.51<br>Q<sub>M</sub> = 3.64</td><td>AICc = 92.17<br>df = 50</td><td>p &lt; 0.001<br>p = 0.06</td><td>Intercept<br>(generalist)<br><strong>Specialists vs generalists</strong></td><td>0.15<br><strong>- 0.26</strong></td><td>0.08<br><strong>0.14</strong></td><td>1.94<br><strong>- 1.91</strong></td><td>0.05<br><strong>0.06</strong></td></tr><tr><td colspan="9"><strong>Incidence rate</strong></td></tr><tr><td>Grand mean</td><td>τ<sup>2</sup> = 0.06<br>Q = 73.18</td><td>AICc = 56.44<br>df = 39</td><td>p &lt; 0.001</td><td>Intercept</td><td>- 0.08</td><td>0.07</td><td>- 1.08</td><td>0.28</td></tr><tr><td>Most parsimonious model</td><td>τ<sup>2</sup> = 0.04<br>Q<sub>E</sub> = 57.75<br>Q<sub>M</sub> = 23.97</td><td>AICc = 43.54<br>df = 38</td><td>p = 0.02<br>p &lt; 0.001</td><td>Intercept<br><strong>MAT</strong></td><td>- 0.42<br><strong>0.03</strong></td><td>0.07<br><strong>0.01</strong></td><td>- 5.67<br><strong>4.9</strong></td><td><strong>&lt; 0.001</strong><br><strong>&lt; 0.001</strong></td></tr><tr><td colspan="9"><strong>Species richness</strong></td></tr><tr><td>Grand mean</td><td>τ<sup>2</sup> = 0.34<br>Q = 99.23</td><td>AICc = 57<br>df = 27</td><td>p &lt; 0.001</td><td>Intercept</td><td>0.36</td><td>0.15</td><td>2.35</td><td>0.02</td></tr></tbody></table>

Mixed-effects models tested the effect of mean annual temperature (MAT), study design (tree plantations or semi-natural forests) and herbivore specialization on the transformed correlation coefficient (Fisher's z-scores) between the diversity of tree species and the four different aspects of herbivory. In each model the intercept denotes the reference level of coefficient estimates, $τ^2$ denotes the variance between study cases, $Q/Q_E$ relate to Cochran's Q-test for residual heterogeneity and $Q_M$ denotes to the omnibus test of model coefficients. Significant parameter estimates are in bold.

doi:10.1371/journal.pone.0165815.t001

Effect sizes for herbivore damage showed a decrease with MAT, indicating a tendency toward associational resistance in warm regions, and associational susceptibility in cold regions (Fig 4), and were negative in experimental (estimated effect size: -0.02 ± 0.08, Table 1) and positive in observational studies (0.16 ± 0.11, Table 1). The abundance of generalist herbivores increased with tree diversity (estimated effect size: 0.15 ± 0.08, Table 1) while it decreases for specialist herbivores (-0.12 ± 0.1, Table 1). The included predictors explained 9.5 and 1.6 percent of variation for herbivore damage and abundance, respectively. As for herbivore incidence rate, effect sizes displayed a significantly positive relationship with MAT which accounted for 29.3 percent of variation. Effect sizes for herbivore species richness were not related to any of the proposed predictors (Table 1, Fig 4).

Addressing publication bias, most of funnel plots were symmetric, indicating that studies reporting positive and negative correlations, with low and high sample sizes were equally likely to be published (S2 Text). The omnibus test of model coefficients indicated that effect sizes did not depend on the metric of tree diversity applied ($Q_M$ = 5.66, df = 3, p = 0.129; S2 Text). We
<figure><img src="image_6.png" /><figcaption>**Fig 3.** Forest plots for the transformed correlation coefficients (Fisher's z-scores) between tree diversity and a) the damage sustained by, b) the abundance of, c) the incidence rate of, d) the species richness of insect herbivores. Each point represents the Fisher's z-score and the approximated confidence interval (= mean ± standard error × 1.96) for an individual study case. Negative values indicate associational resistance while positive values indicate associational susceptibility. Grand mean effect sizes, together with their 95% bootstrap confidence intervals, are shown in black at the bottom of each forest plot.</figcaption></figure>

doi:10.1371/journal.pone.0165815.g003

found a negative correlation between the year of publication and MAT (r = -0.18, p = 0.01), showing that studies in cooler climates were carried out later, but not between the year of publication and the actual effect sizes (r = -0.01, p = 0.87), indicating that there was no temporal shift in our understanding of associational effects.

## Discussion

Our meta-analysis does not support previous claims of pervasive, globally consistent associational resistance (AR) to insect herbivores in mixed forests [8, 9]. Importantly, we show that several sources of variation in the strength and direction of associational effects have been
<figure><img src="image_7.png" /><figcaption>**Fig 4.** Relationships between mean annual temperature (MAT) and the transformed correlation coefficients (Fisher's z-scores) between the diversity of tree species and a) the damage sustained by, b) the abundance of, c) the incidence rate of, d) the species richness of insect herbivores. Each point represents an individual study case for which negative values indicate associational resistance while positive values indicate associational susceptibility. The size of each point indicates its weight for estimating the regression slope (solid line) and the corresponding approximated credible interval (dotted lines) in a mixed-effects meta-regression model. Coefficient estimates for MAT are reported in Table 1.</figcaption></figure>

doi:10.1371/journal.pone.0165815.g004

overlooked, in particular mean annual temperature (MAT). We found no negative relation-
ship between tree diversity and herbivory, regardless if the dataset was split according to the
aspect of herbivory or the metric of tree diversity and, thus, have to reject our first hypothesis.

Instead, and in confirmation of our second hypothesis, we demonstrated that grand mean
effect sizes differed between the four aspects of herbivory. Notably, herbivore species richness,
an aspect not analysed in previous syntheses, increased with tree diversity. Importantly, we
detected a significant relationship between MAT and the response of herbivore damage and
incidence rate to increasing tree diversity. This finding partly confirms our third hypothesis
and strongly improves our understanding of heterogeneity among studies.
We found that both aspects of herbivory displayed quite opposite tendencies. With increas-
ing MAT, diversity effects on herbivore damage shifted from associational susceptibility (AS)
in cold regions towards AR in warmer regions, whereas the opposite pattern was encountered
for herbivore incidence rate, with a tendency towards AS in warmer regions.

As our study was based on correlation coefficients which become significant if any increase
in tree diversity is accompanied by a steady increase or decrease in herbivory, we cannot make
any statement on the change in the absolute amount of herbivory. However, we demonstrated
that future syntheses on the strength of associational effects in forest systems should differenti-
ate between the aspects of herbivory and consider the environmental context of the study site.
This is critical because different aspects of tree diversity may influence damage (i.e. actual con-
sumption) and incidence rate. For instance, incidence rate could reflect herbivore foraging
behaviour, while damage additionally could include tree diversity effects on food quality and
herbivore survival.

In search for factors that underlie the documented biogeographic gradient, we ask for more
research on how associational effects depend on and follow global gradients in i) regional spe-
cies diversity of herbivores and host plants, ii) herbivore density and pressure, especially dur-
ing times of outbreaks [25, 28–30, 52–56], iii) herbivore specialization [4, 5, 8, 9, 20, 57], iv) tree defences [25, 27, 29, 31, 52, 58] and v) abiotic factors affecting tree or herbivore develop-
ment (e.g. precipitation, climate stability, CO₂-concentration, UVB-radiation).

The literature on global gradients in herbivore pressure is controversial, providing exam-
ples of either increasing [25, 29, 52–54], unaffected [25, 29, 53, 55] or decreasing [28, 56] pres-
sure towards warmer climates. Regarding the level of background herbivory, e.g. the regular
loss of woody plant foliage, a recent global analysis even pointed out that background herbi-
vore pressure could show a nonlinear global pattern that is peaking in temperate regions [30].

The highest levels of herbivore pressure can be studied during insect outbreaks when herbi-
vores face strong intraspecific competition and resource depletion. If, during outbreaks, her-
bivorous individuals are forced to forage less selectively, this might shift associational effects
from AR to AS, regardless of the composition of tree species mixtures. Unfortunately, our
dataset included only 5 study cases (3 original studies) that measured herbivory during out-
breaks (S1 File) and we are not aware of any study that investigated global gradients in the fre-
quency and intensity of insect herbivore outbreaks in forest ecosystems.

Global patterns in the relative abundance of generalist vs. specialist insect herbivores, how-
ever, had already been addressed. Here, Novotny et al. [59] and Schuldt et al. [4] proposed that
tropical and subtropical forest are dominated by generalist herbivores, which, according to
previous meta-analyses [8, 9], are not affected or even benefit from the diversification of tree
stands. A higher proportion of generalist species could explain the shift in herbivore incidence
rates from AR to AS with increasing MAT. However, more comprehensive and global analyses
concluded that the proportion of generalist species actually decreases towards tropical regions
[20, 57].

Beside the difficulty to relate the documented shifts in associational effects with MAT to a
single, underlying factor, any global synthesis of plant-plant-herbivore interactions might fur-
thermore be confronted with gradients that are non-linear, non-additive or interrelated (i.e.
show interactions such as latitudinal changes in plant defences being counterbalanced by
higher herbivore pressure [29, 52]). In addition to the abiotic and biotic factors already men-
tioned, many decisions on the design of a study, such as the spatial scale of the investigated
plant neighbourhood [60], the age of the forest stand [61] and the sampling date [62] are likely
to determine the sampling success of insect herbivore communities and thus impact the con-
clusions on the direction and strength of associational effects.
Given the multitude of potentially confounding factors, it is hardly surprising that our study documented a high amount of unexplained heterogeneity, and thus, highlights the limits of our meta-analytical approach. Here, research co-operations, such as the globally distributed network of tree diversity experiments (www.treedivnet.ugent.be), can offer great future opportunities to experimentally study associational effects along replicated global abiotic and biotic gradients [63].

**Conclusion**

Our results indicate that studies of associational effects on herbivory might need to consider the biogeographical context in which plant-plant-herbivore interactions occur. Yet, along such gradients, joint impacts of insect herbivore diversity, pressure, specialization and abiotic factors on a global gradient of associational effects are difficult to disentangle.

We recommend that information on the damage, incidence, abundance and diversity, together with the identity of herbivores, be systematically recorded in observational and experimental tree diversity studies to provide a sounder understanding of mechanisms involved in AR and AS.

The inability to replicate major findings of previous meta-analyses, namely the significance of AR of mixed tree stands to insect herbivory [8, 9], could even indicate that associational effects are nonlinear, thus adding another layer of complexity.

A better understanding of the mechanisms at play will require addressing the identity, functional characteristics, density and diversity of both, insect herbivores and focal tree species, in a systematic way along large geographical gradients.

**Supporting Information**

**S1 PRISMA Checklist.**

(DOC)

**S1 File. Coding table for this analysis.** Including a list of publications and researchers that provided the data for this meta-analysis.

(XLSX)

**S1 Text. Search protocol.**

(DOCX)

**S2 Text. Funnel plot diagrams.**

(DOCX)

**Acknowledgments**

The work of S. Kambach was funded by the Helmholtz Centre for Environmental Research-UFZ, the Helmholtz Research School for Ecosystem Services under Changing Land-use and Climate (ESCALATE) and the Martin Luther University Halle-Wittenberg. We are indebted to the researchers who shared their raw data with us. Publications and researches that provided data for this meta-analysis are listened in S1 File.

**Author Contributions**

**Conceptualization:** SK HB.

**Data curation:** SK.
**Formal analysis:** SK BC HB IK.

**Funding acquisition:** HB IK.

**Investigation:** SK.

**Methodology:** SK BC HB IK.

**Software:** SK BC.

**Supervision:** HB IK.

**Visualization:** SK.

**Writing – original draft:** SK.

**Writing – review & editing:** BC HB IK.

## References

1. Boyd IL, Freer-Smith PH, Gilligan CA, Godfray HCJ. The consequence of tree pests and diseases for ecosystem services. Science. 2013; 342: 1235773. doi: 10.1126/science.1235773 PMID: 24233727

2. Price PW, Bouton CE, Gross P, McPheron BA, Thompson JN, Weis AE. Interactions among three trophic levels: influence of plants on interactions between insect herbivores and natural enemies. Annu Rev Ecol Syst. 1980; 11: 41–65. doi: 10.1146/annurev.es.11.110180.000353

3. Coley PD. Herbivory and defensive characteristics of tree species in a lowland tropical forest. Ecol Monogr. 1983; 53: 209–229. doi: 10.2307/1942495

4. Schuldt A, Baruffol M, Bohnke M, Bruelheide H, Hardtle W, Lang AC, et al. Tree diversity promotes insect herbivory in subtropical forests of south-east China. J Ecol. 2010; 98: 917–926. doi: 10.1111/j.1365-2745.2010.01659.x PMID: 20852667

5. Schuldt A, Bruelheide H, Härdtle W, Assmann T, Li Y, Ma K, et al. Early positive effects of tree species richness on herbivory in a large-scale forest biodiversity experiment influence tree growth. J Ecol. 2015; 103: 563–571. doi: 10.1111/1365-2745.12396 PMID: 26690688

6. Underwood N, Inouye BD, Hambäck PA. A conceptual framework for associational effects: when do neighbors matter and how would we know? Q Rev Biol. 2014; 89: 1–19. doi: 10.1086/674991 PMID: 24672901

7. Barbosa P, Hines J, Kaplan I, Martinson H, Szczepaniec A, Szendrei Z. Associational resistance and associational susceptibility: having right or wrong neighbors. Annu Rev Ecol Evol Syst. 2009; 40: 1–20. doi: 10.1146/annurev.ecolsys.110308.120242

8. Castagneyrol B, Jactel H, Vacher C, Brockerhoff EG, Koricheva J. Effects of plant phylogenetic diversity on herbivory depend on herbivore specialization. J Appl Ecol. 2014; 51: 134–141. doi: 10.1111/1365-2664.12175

9. Jactel H, Brockerhoff EG. Tree diversity reduces herbivory by forest insects. Ecol Lett. 2007; 10: 835–848. doi: 10.1111/j.1461-0248.2007.01073.x PMID: 17663717

10. White JA, Whitham TG. Associational susceptibility of cottonwood to a box elder herbivore. Ecology. 2000; 81: 1795–1803. doi: 10.1890/0012-9658(2000)081[1795:ASOCTA]2.0.CO;2

11. Otway SJ, Hector A, Lawton JH. Resource dilution effects on specialist insect herbivores in a grassland biodiversity experiment. J Anim Ecol; 2005; 74: 234–240. doi: 10.1111/j.1365-2656.2005.00913.x

12. Hambäck PA, Inouye BD, Andersson P, Underwood N. Effects of plant neighborhoods on plant–herbivore interactions: resource dilution and associational effects. Ecology. 2014; 95: 1370–1383. doi: 10.1890/13-0793.1 PMID: 25000768

13. Castagneyrol B, Giffard B, Péré C, Jactel H. Plant apparency, an overlooked driver of associational resistance to insect herbivory. J Ecol. 2013; doi: 10.1111/1365-2745.12055

14. Jactel H, Birgersson G, Andersson S, Schlyter F. Non-host volatiles mediate associational resistance to the pine processionary moth. Oecologia. 2011; 166: 703–711. doi: 10.1007/s00442-011-1918-z PMID: 21311911

15. Brown BJ, Ewel JJ. Herbivory in complex and simple tropical successional ecosystems. Ecology. 1987; 68: 108–116. doi: 10.2307/1938810
16. Unsicker SB, Oswald A, Köhler G, Weisser WW. Complementarity effects through dietary mixing enhance the performance of a generalist insect herbivore. Oecologia. 2008; 156: 313–324. doi: 10.1007/s00442-008-0973-6 PMID: 18301924

17. Raubenheimer D, Simpson SJ. Integrating nutrition: a geometrical approach. Proceedings of the 10th International Symposium on Insect-Plant Relationships. Springer Netherlands; 1999. pp. 67–82. 10.1007/978-94-017-1890-5_8

18. Bernays EA, Bright KL, Gonzalez N, Angel J. Dietary mixing in a generalist herbivore: tests of two hypotheses. Ecology. 1994; 75: 1997–2006. doi: 10.2307/1941604

19. Tahvanainen JO, Root RB. The influence of vegetational diversity on the population ecology of a specialized herbivore, Phyllotreta cruciferae (Coleoptera: Chrysomelidae). Oecologia. 1972; 10: 321–346. doi: 10.1007/BF00345736

20. Forister ML, Novotny V, Panorska AK, Baje L, Basset Y, Butterill PT, et al. The global distribution of diet breadth in insect herbivores. Proc Natl Acad Sci. 2015; 112: 442–447. doi: 10.1073/pnas.1423042112 PMID: 25548168

21. Kreft H, Jetz W. Global patterns and determinants of vascular plant diversity. Proc Natl Acad Sci. 2007; 104: 5925–5930. doi: 10.1073/pnas.0608361104 PMID: 17379667

22. Gillman LN, Wright SD, Cusens J, McBride PD, Malhi Y, Whittaker RJ. Latitude, productivity and species richness. Glob Ecol Biogeogr. 2015; 24: 107–117. doi: 10.1111/geb.12245

23. Price PW, Fernandes GW, Lara ACF, Brawn J, Barrios H, Wright MG, et al. Global patterns in local number of insect galling species. J Biogeogr. 1998; 25: 581–591. doi: 10.1046/j.1365-2699.1998.2530581.x

24. Hillebrand H. On the generality of the latitudinal diversity gradient. Am Nat. 2004; 163: 192–211. doi: 10.1086/381004 PMID: 14970922

25. Salazar D, Marquis RJ. Herbivore pressure increases toward the equator. Proc Natl Acad Sci. 2012; 109: 12616–12620. doi: 10.1073/pnas.1202907109 PMID: 22802664

26. Moles AT, Bonser SP, Poore AGB, Wallis IR, Foley WJ. Assessing the evidence for latitudinal gradients in plant defence and herbivory. Funct Ecol. 2011; 25: 380–388. doi: 10.1111/j.1365-2435.2010.01814.x

27. Pearse IS, Hipp AL. Global patterns of leaf defenses in oak species. Evolution. 2012; 66: 2272–2286. doi: 10.1111/j.1558-5646.2012.01591.x PMID: 22759301

28. Adams JM, Zhang Y. Is there more insect folivory in warmer temperate climates? A latitudinal comparison of insect folivory in eastern North America. J Ecol. 2009; 97: 933–940. doi: 10.1111/j.1365-2745.2009.01523.x

29. Pennings SC, Ho C-K, Salgado CS, Wleski K, Davé N, Kunza AE, et al. Latitudinal variation in herbivore pressure in Atlantic coast salt marshes. Ecology. 2009; 90: 183–195. doi: 10.1890/08-0222.1 PMID: 19294924

30. Kozlov MV, Lanta V, Zverev V, Zvereva EL. Global patterns in background losses of woody plant foliage to insects. Glob Ecol Biogeogr. 2015; doi: 10.1111/geb.12347

31. Rodríguez-Castañeda G. The world and its shades of green: a meta-analysis on trophic cascades across temperature and precipitation gradients. Glob Ecol Biogeogr. 2013; 22: 118–130. doi: 10.1111/j.1466-8238.2012.00795.x

32. Bale JS, Masters GJ, Hodkinson ID, Awmack C, Bezemer TM, Brown VK, et al. Herbivory in global climate change research: direct effects of rising temperature on insect herbivores. Glob Chang Biol. 2002; 8: 1–16. doi: 10.1046/j.1365-2486.2002.00451.x

33. Lemoine NP, Burkepile DE, Parker JD. Variable effects of temperature on insect herbivory. PeerJ. 2014; 2: e376. doi: 10.7717/peerj.376/table-1 PMID: 24860701

34. Clissold FJ, Coggan N, Simpson SJ. Insect herbivores can choose microclimates to achieve nutritional homeostasis. J Exp Biol. 2013; 216: 2089–2096. doi: 10.1242/jeb.078782 PMID: 23430995

35. Verschut TA, Becher PG, Anderson P, Hambäck PA. Disentangling associational effects: both resource density and resource frequency affect search behaviour in complex environments. Funct Ecol. 2016; 1365–2435. doi: 10.1111/1365-2435.12670

36. Kim TN, Underwood N. Plant neighborhood effects on herbivory: damage is both density and frequency dependent. Ecology. 2015; 96: 1431–1437. doi: 10.1890/14-1097.1 PMID: 26236855

37. Rhainds M, English-Loeb G. Testing the resource concentration hypothesis with tarnished plant bug on strawberry: density of hosts and patch size influence the interaction between abundance of nymphs and incidence of damage. Ecol Entomol. 2003; 28: 348–358. doi: 10.1046/j.1365-2311.2003.00508.x

38. Stoepler TM, Rehill B. Forest habitat, not leaf phenotype, predicts late-season folivory of Quercus alba saplings. Funct Ecol. 2012; 26: 1205–1213. doi: 10.1111/j.1365-2435.2012.02033.x

39. Harzing AW. Publish or Perish. 2014. Available: http://www.harzing.com.
40. Moher D, Liberati A, Tetzlaff J, Altman DG, PRISMA Group. Preferred reporting items for systematic reviews and meta-analyses: the PRISMA statement. Int J Surg. Am Coll Physicians; 2010; 8: 336–341. doi: 10.1016/j.ijsu.2010.02.007 PMID: 20171303

41. Fisher RA. On the "probable error" of a coefficient of correlation deduced from a small sample. Metron; 1921; 1(1): 3–32.

42. Borenstein M, Hedges LV, Higgins JPT, Rothstein HR. Introduction to meta-analysis. John Wiley & Sons; 2011.

43. Koricheva J, Gurevitch J, Mengersen K. Handbook of meta-analysis in ecology and evolution. Princeton University Press; 2013.

44. Rasband WS. ImageJ. 2014. Available: 2014. Available: http://imagej.nih.gov/ij/.

45. Hijmans RJ, Cameron S, Parra J, Jones P, Jarvis A, Richardson K. WorldClim. 2008. Available: http://www.worldclim.org.

46. Sterne JAC, Egger M. Regression methods to detect publication and other bias in meta-analysis. Publication bias in meta-analysis: Prevention, assessment and adjustments. John Wiley & Sons. 2005; 99–110. doi: 10.1002/0470870168.ch6

47. R Core Team. R: A Language and environment for statistical computing. 2013. Available: http://www.R-project.org.

48. Viechtbauer W. Conducting meta-analyses in R with the metafor package. J Stat Softw. 2010; 36: 1–48. doi: 10.18637/jss.v036.i03

49. Oksanen J, Blanchet FG, Kindt R, Legendre P, Minchin PR, O'Hara RB, et al. vegan: community ecology package. 2011. Available: http://vegan.r-forge.r-project.org.

50. Hijmans RJ. raster: Geographic data analysis and modeling. 2013. Available: http://CRAN.R-project.org/package=raster.

51. Wickham H. ggplot2: elegant graphics for data analysis. Springer New York; 2009. doi: 10.1111/j.1467-985X.2010.00676_9.x

52. Coley PD, Barone JA. Herbivory and plant defenses in tropical forests. Annu Rev Ecol Syst. 1996; 27: 305–335. doi: 10.1146/annurev.ecolsys.27.1.305

53. Kozlov MV. Losses of birch foliage due to insect herbivory along geographical gradients in Europe: a climate-driven pattern? Clim Change. 2008; 87: 107–117. doi: 10.1007/s10584-007-9348-y

54. Adams JM, Zhang Y, Basri M, Shukor N. Do tropical forest leaves suffer more insect herbivory? A comparison of tropical versus temperate herbivory, estimated from leaf litter. Ecol Res. 2009; 24: 1381–1392. doi: 10.1007/s11284-009-0623-4

55. Andrew NR, Hughes L. Herbivore damage along a latitudinal gradient: relative impacts of different feeding guilds. Oikos. 2005; 108: 176–182. doi: 10.1111/j.0030-1299.2005.13457.x

56. del-Val E, Armesto JJ. Seedling mortality and herbivory damage in subtropical and temperate populations: testing the hypothesis of higher herbivore pressure toward the tropics. Biotropica. 2010; 42: 174–179. doi: 10.1111/j.1744-7429.2009.00554.x

57. Dyer LA, Singer MS, Lill JT, Stireman JO, Gentry GL, Marquis RJ, et al. Host specificity of Lepidoptera in tropical and temperate forests. Nature. 2007; 448: 696–699. doi: 10.1038/nature05884 PMID: 17687325

58. Marquis RJ, Ricklefs RE, Abdala-Roberts L. Testing the low latitude/high defense hypothesis for broad-leaved tree species. Oecologia. 2012; 169: 811–820. doi: 10.1007/s00442-012-2249-4 PMID: 22271200

59. Novotny V, Basset Y, Miller SE, Weiblen GD, Bremer B, Cizek L, et al. Low host specificity of herbivorous insects in a tropical forest. Nature. 2002; 416: 841–844. doi: 10.1038/416841a PMID: 11976681

60. Champagne E, Tremblay J-P, Côté SD. Spatial extent of neighboring plants influences the strength of associational effects on mammal herbivory. Ecosphere. 2016; 7(6): e0137. doi: 10.1002/ecs2.1371

61. Jeffries JM, Marquis RJ, Forkner RE. Forest age influences oak insect herbivore community structure, richness, and density. Ecol Appl. 2006; 16: 901–912. doi: 10.1890/1051-0761(2006)016[0901:FAIOIH]2.0.CO;2 PMID: 16826990

62. van Asch M, Visser ME. Phenology of forest caterpillars and their host trees: the importance of synchrony. Annu Rev Entomol. 2007; 52: 37–55. doi: 10.1146/annurev.ento.52.110405.091418 PMID: 16842033

63. Verheyen K, Vanhellemont M, Auge H, Baeten L, Baraloto C, Barsoum N, et al. Contributions of a global network of tree diversity experiments to sustainable forest plantations. Ambio. 2015; doi: 10.1007/s13280-015-0685-1 PMID: 26264716