# Bottom-up and top-down effects of tree species diversity on leaf insect herbivory

Bastien Castagneyrol, Damien Bonal, Maxime Damien, Hervé Jactel, Céline Meredieu, Evalyne W. Muiruri, Luc Barbaro

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

## To cite this version:

Bastien Castagneyrol, Damien Bonal, Maxime Damien, Hervé Jactel, Céline Meredieu, et al., Bottom-up and top-down effects of tree species diversity on leaf insect herbivory. Ecology and Evolution, 2017, 7 (10), pp.3520-3531. 10.1002/ece3.2950. hal-01559963

HAL Id: hal-01559963

https://univ-rennes.hal.science/hal-01559963

Submitted on 21 Jun 2018

HAL is a multi-disciplinary open access archive for the deposit and dissemination of scientific research documents, whether they are published 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>

Distributed under a Creative Commons Attribution 4.0 International License
# Bottom-up and top-down effects of tree species diversity on leaf insect herbivory

Bastien Castagneyrol<fnref n="1" /> ID | Damien Bonal<fnref n="2" /> | Maxime Damien<fnref n="3" /> | Hervé Jactel<fnref n="1" /> |
Céline Meredieu<fnref n="1" /> | Evalyne W. Muiruri<fnref n="4" />,<fnref n="5" /> | Luc Barbaro<fnref n="1" />,<fnref n="6" />

**Correspondence**

Bastien Castagneyrol, INRA - UMR 1202 BIOGECO, Cestas, France.

Email: bastien.castagneyrol@inra.fr

**Funding information**

French Ministry of Ecology; French National Research Agency (ANR), Grant/Award

Number: ANR-11-LABX-0002-01

**Abstract**

The diversity of plant neighbors commonly results in direct, bottom-up effects on herbivore ability to locate their host, and in indirect effects on herbivores involving changes in plant traits and a top-down control by their enemies. Yet, the relative contribution of bottom-up and top-down forces remains poorly understood. We also lack knowledge on the effect of abiotic constraints such as summer drought on the strength and direction of these effects. We measured leaf damage on pedunculate oak (*Quercus robur*), alone or associated with birch, pine or both in a long-term tree diversity experiment (ORPHEE), where half of the plots were irrigated while the other half remained without irrigation and received only rainfall. We tested three mechanisms likely to explain the effects of oak neighbors on herbivory: (1) Direct bottom-up effects of heterospecific neighbors on oak accessibility to herbivores, (2) indirect bottom-up effects of neighbors on the expression of leaf traits, and (3) top-down control of herbivores by predators. Insect herbivory increased during the growth season but was independent of neighbor identity and irrigation. Specific leaf area, leaf toughness, and thickness varied with neighbor identity while leaf dry matter content or C:N ratio did not. When summarized in a principal component analysis (PCA), neighbor identity explained 87% of variability in leaf traits. PCA axes partially predicted herbivory. Despite greater rates of attack on dummy caterpillars in irrigated plots, avian predation, and insect herbivory remained unrelated. Our study suggests that neighbor identity can indirectly influence insect herbivory in mixed forests by modifying leaf traits. However, we found only partial evidence for these trait-mediated effects and suggest that more attention should be paid to some unmeasured plant traits such as secondary metabolites, including volatile organic compounds, to better anticipate the effects of climate change on plant-insect interactions in the future.

**KEYWORDS**

associational effects, biodiversity, drought, dummy caterpillars, leaf traits, predation, *Quercus robur*

<footnote n="1">¹BIOGECO, INRA, Univ. Bordeaux, 33610 Cestas, France</footnote>

<footnote n="2">²EEF, INRA, Université de Lorraine, 54280 Champenoux, France</footnote>

<footnote n="3">³ECOBIO, UMR CNRS 6553, Université de Rennes, 35042 Rennes, France</footnote>

<footnote n="4">⁴School of Biological Sciences, Royal Holloway University of London, Egham, Surrey TW20 0EX, UK</footnote>

<footnote n="5">⁵Department of Biosciences, Durham University, Stockton Road, Durham, DH1 3LE, UK</footnote>

<footnote n="6">⁶Dynafor, INPT, EI Purpan, INRA, Université de Toulouse, 31320 Auzeville, France</footnote>
# 1 | INTRODUCTION

Plant diversity is a key driver of insect herbivory in grassland, agricultural, and forest ecosystems (Allan et al., 2013; Broad, Schellhorn, Lisson, & Mendham, 2008; Castagneyrol et al., 2014), because heterospecific neighbors may either decrease (associational resistance) or increase (associational susceptibility) the likelihood of focal plants being attacked by insect herbivores (Barbosa et al., 2009; Root, 1973; Tahvanainen & Root, 1972; White & Whitham, 2000). Several reviews and meta-analyses suggest that associational resistance is the most common pattern (Andow, 1991; Barbosa et al., 2009), particularly in forest ecosystems (Castagneyrol et al., 2014; Jactel & Brockerhoff, 2007). However, despite the large body of evidence on associational effects (Moreira et al., 2016; Underwood, Inouye, & Hambäck, 2014), gaps in our understanding of underlying mechanisms hinder the development of a predictive framework of herbivore responses to tree diversity. In addition, due to the rising threats to forests that result from climate change, like more severe or frequent outbreaks of pest insects (Logan, Régnière, & Powell, 2003), it is critical to better evaluate the relative contribution of factors responsible for associational resistance or susceptibility under contrasting abiotic constraints.

Associational effects depend on the identity of heterospecific neighbors, both directly through bottom-up effects on focal host accessibility (e.g., Castagneyrol, Giffard, Péré, & Jactel, 2013) and indirectly through changes in focal host quality (e.g., Kos, Bukovinszky, Mulder, & Bezemer, 2015; Kostenko, Mulder, Courbois, & Bezemer, 2016; Mraja, Unsicker, Reichelt, Gershenzon, & Roscher, 2011) and top-down control of herbivore populations by predators (e.g., Muiruri et al. 2015). However, the relative importance of top-down and bottom-up effects remains unclear, in particular because few studies to date addressed them simultaneously (Abdala-Roberts et al., 2016; Moreira, Mooney, Zas, & Sampedro, 2012).

Three alternative mechanisms have been proposed to explain associational resistance. First, associational resistance has been proposed to result primarily from a reduced ability of herbivores to locate and reach their host plants among heterospecific neighbors. This may be due to patches with greater plant diversity being less attractive than monospecific patches (i.e., the resource concentration hypothesis, Root, 1973; Andersson, Löfstedt, & Hambäck, 2013) or from heterospecific neighbors reducing the physical (Castagneyrol et al., 2013; Damien et al., 2016; Floater & Zalucki, 2000) or chemical (Jactel, Birgersson, Andersson, & Schlyter, 2011; Zhang & Schlyter, 2004) apparenty of host plants. Alternatively, the attractant-decoy hypothesis predicts that herbivores can be diverted from a given plant and aggregate on more apparent, more attractive, or more palatable neighbors (Atsatt & O'Dowd, 1976; Hahn & Orrock, 2016).

Second, heterospecific neighbors may affect insect herbivory on a focal plant through changes in host food quality, although this is less well documented (but see, e.g., Walter et al., 2011; Kos et al., 2015). Because the amount of insect herbivory depends on plant traits such as leaf toughness, water content, C:N ratio, and secondary metabolites (Loranger et al., 2013; Moreira, Abdala-Roberts, Parra-Tabla, & Mooney, 2014; Pearse, 2011), any biotic or abiotic factor affecting the expression

of these plant traits may also affect herbivores. For instance, competition or facilitation among the focal plant and its neighbors may change the nutritional value of the host plant tissues (Kos et al., 2015; Schädler, Brandl, & Haase, 2007; Walter et al., 2011) or its production of secondary metabolites to defend against herbivores (Moreira et al., 2014).

Third, the enemies' hypothesis posits that predators and parasitoids of herbivores are more abundant and diverse in species-rich plant communities (Elton, 1958; Root, 1973; Schuldt et al., 2011; Straub et al., 2014). Associational resistance may therefore also result from an enhanced top-down control of herbivores by predation or parasitism in species-rich plant communities (Riihimäki, Kaitaniemi, Koricheva, & Vehviläinen, 2004). Besides its effect on abundance and richness of enemies, plant diversity may indirectly influence the strength of herbivory suppression by predators. For instance, greater predation is expected when herbivores are more exposed to their enemies as they spend more time foraging for less accessible resource (Straub et al., 2014). Plant diversity may also modify the magnitude of top-down effects by altering the proportion of generalist versus specialist herbivores because herbivore communities dominated by generalist herbivores are more sensitive to predation (Singer et al., 2014).

The strength and direction of plant-herbivores-enemies (predators or parasitoids) interactions is expected to change along environmental gradients (Bauerfeind & Fischer, 2013; Péré, Jactel, & Kenis, 2013; Rodríguez-Castañeda, 2013; Walter et al., 2011). Rooted in the plant stress hypothesis (White, 1974) that predicts increasing plant susceptibility to herbivores with higher water stress, several studies suggest that drought favors leaf-feeding herbivores and reduces damage caused by sap-feeders (Huberty & Denno, 2004; Jactel et al., 2012). In particular, drought might affect the palatability of plant tissues, notably through change in carbohydrate content and C:N ratio (Jactel et al., 2012; Walter et al., 2011). In addition to changes in leaf traits interfering with plant quality, drought may also indirectly affect herbivory by modifying predator attraction (Aslam, Johnson, & Karley, 2013; Staley et al., 2006; Weldegergis, Zhu, Poelman, & Dicke, 2015). Yet, although underlying mechanisms remain unclear, it is increasingly acknowledged that neighbors can also modify plant's response to drought (Forrester, Theiveyanathan, Collopy, & Marcar, 2010; Grossiord et al., 2014). Insect herbivory may therefore be affected by both plant diversity and water stress, as well as their potential interacting effects (Walter et al., 2011).

Experimental evidence that plant neighborhood mediates the effect of water stress on herbivores and their enemies—or the opposite—is still lacking given the difficulties to control for both the composition of tree neighborhood and climatic variables. In this study, we searched for mechanisms responsible for associational effects under contrasting abiotic conditions. We measured insect herbivory as the percentage of defoliation by chewing herbivores on pedunculate oak (*Quercus robur*) in a tree diversity experiment with a factorial design crossing the identity of their neighbors (homospecific vs. heterospecific) with an irrigation treatment (rainfall only vs. rainfall plus irrigation). We specifically addressed mechanisms underlying neighbors and drought effects on insect herbivory. We predicted that: (1) herbivory is greater in tree monocultures where individual oaks are more apparent than in tree mixtures where they are protected by taller nonoak neighbors, (2)
tree neighbors indirectly mediate the expression of traits involved in herbivore-oak interaction; (3) heterospecific neighbors increase top-down predation by natural enemies, and (4) water stress changes the strength of top-down and bottom-up processes. Overall, this study addresses how drought indirectly mediates tree neighbor effects on insect herbivory via changes in leaf traits and top-down control by herbivores' enemies. We provide partial support to the hypothesis of indirect, bottom-up, effect of tree species diversity on herbivory through changes in leaf traits, regardless of drought.

# 2 | MATERIALS AND METHODS

## 2.1 | Ethics statement

This study did not involve manipulations of humans or animals. No specific permissions were required for our field work. The study did not involve endangered or protected species.

## 2.2 | The ORPHEE experiment

This study was carried out in the ORPHEE experiment in SW France (44°440 N, 00°460 W). The design has been fully described in Castagneyrol et al. (2013), and we will thus only provide here a brief overview. In 2008, eight blocks were established covering 12 ha, with 32 plots in every block corresponding to the 31 possible combinations of one to five tree species (*Betula pendula*, *Quercus robur*, *Q. pyrenaica*, *Q. ilex*, and *Pinus pinaster*) with an additional replicate of the five-species plot. Each plot contained 10 rows of 10 trees planted 2 m apart (100 trees). Tree species mixtures were established according to a substitutive design, keeping tree density and the identity of tree neighbors equal across plots. Within plots, individual trees from different species were planted in a regular alternate pattern, such that a tree from a given species had at least one neighbor from each of the other species within a 2-m radius (Fig. S1 and www.facebook.com/orpheeexperiment).

In this study, we focused on a subsample of plots to reconstruct a tree diversity gradient spanning from the monoculture of pedunculate oak, the two-two-species mixtures associating oak with pine or birch, and the three-species mixtures of oak, birch, and pine, for a total of four compositions. These plots were chosen in order (1) to span a gradient of oak apparency by increasing the number of neighbors taller than oaks (i.e., high apparency in monocultures, intermediate apparency in oak-birch and oak-pine mixtures and low apparency in oak-birch-pine mixtures, Castagneyrol et al., 2013) and (2) to contrast functional diversity by associating oak with a broadleaved or a conifer species. Within the selected 32 plots (4 plots × 8 blocks), we sampled at random six individual trees of the 36 innermost trees of each plot, for a total of 192 sampled oaks.

In 2015, half of the blocks were irrigated (Fig. S2). Irrigation consisted in sprinkling *ca* 42 m³ per night and per block from early May to late September, corresponding to *ca* 3 mm/day per plot. This volume was calculated based on regional climatic data (evapotranspiration) and was assumed to avoid any soil water deficit in the irrigated blocks during the entire growing season.

## 2.3 | Insect herbivory

Insect herbivory was assessed twice to test for a season effect, in early June and early August 2015 by visual inspection of 30 leaves per sampled oak (Johnson, Bertrand, & Turcotte, 2016). Three branches were selected at the top, middle, and bottom of each tree. Five leaves were then randomly chosen at the tip, and five at the basis of each branch. The percentage of leaf area removed (LAR) by insect defoliators (chewers and skeletonizers) was estimated on each sampled leaf by a unique observer all along the experiment (BC) using seven classes (0%, 1%-5%, 6%-15%, 16%-25%, 26%-50%, 51%-75% and >75% LAR) and then averaged per sampled tree using the midpoint of each damage class.

## 2.4 | Leaf traits

We measured five leaf traits known to significantly influence insect herbivory: specific leaf area (SLA), leaf dry matter content (LDMC), leaf toughness and thickness, and leaf C:N ratio (e.g., Loranger et al., 2013; Pearse, 2011).

Leaf traits may vary in response to the identity of oak neighbors, herbivory, or both. To avoid confounding the effect of tree species neighbors and insect herbivory on leaf traits, half of the sampled trees were sprayed every two weeks with a broad action spectrum insecticide (5% λ-cyhalothrin, KARATE®, Syngenta, diluted at 15 g/hl). Sprayed and unsprayed oaks were at least 4 m apart to reduce accidental drift of insecticide onto control oaks.

In June 2015, we randomly collected six sunlit, mature, and fully expanded leaves from the top part of the canopy on each target oak (n = 6 per plot, for a total of 192 trees), stored them in zipped plastic bags and immediately put them into a cool box. In order to standardize trait measurements, leaves were rehydrated for 48 hr (Pérez-Harguindeguy et al., 2013) and then scanned and weighted individually. Leaf toughness was assessed as the resistance to piercing, with six measurements per leaf. Measurements were conducted with a dial tension gauge model with peak hold (Mitutoyo Messgeräte Leonberg GmbH, Leonberg, Germany). Leaves were then dried until constant weight (at least 48 hr at 60°C) and weighted again to obtain LDMC (mg/g) and SLA (cm²/g). Leaf thickness was derived from SLA and LDMC (Pérez-Harguindeguy et al., 2013). In early August 2015, we again collected six leaves per tree as explained above. Leaves were dried and then bulked for each tree and finely grinded. Leaf carbon isotope composition (δ¹³C, per mil), foliar carbon (C%, %), and nitrogen (N%, %) contents were measured on these bulked samples. Isotope and elementary analyses were performed at the INRA Nancy Technical Facility of Functional Ecology (OC 081) with an EA/GA-IRMS (Carlo Erba, Elementar, Finnigan, Isoprime, Bremen, Germany). C% and N% were used to calculate C:N ratios. Together with the punctual measurement of leaf water potential (see below), δ¹³C was used to confirm that irrigation actually alleviated water stress.

## 2.5 | Tree water status

To assess the effectiveness of irrigation at alleviating water stress, we compared tree water status of oak trees by measuring predawn leaf
water potential ($\psi_w$) in early August 2015 on 1–3 oaks per plot. $\psi_w$ is used as a proxy of the spatially integrated water potential of the soil explored by roots. We used a Scholander-type pressure chamber (model 1000; PMS Instruments, Corvalis, OR) to measure $\psi_w$ on one leaf per oak tree. Leaves were sampled between 4:30 and 6:30 a.m. local time (before sunrise) with clippers.

## 2.6 | Insectivory

Total and avian insect predation rates were assessed using a standardized method (e.g., Barbaro et al., 2014; Mäntylä et al., 2008; Muiruri, Rainio, & Koricheva, 2016). We modeled 2-cm-long dummy caterpillars made with green modeling clay to mimic larvae of locally abundant moth species, notably Geometridae and Noctuidae. Dummy caterpillars were secured on oak branches with a thin metallic wire. Three to six oaks were randomly selected at the center of each sampled plot. We installed three caterpillars per oak, at the tip of three different branches at ca 1.5 m from the ground. Insectivory was assessed twice, after 3 and 6 weeks of exposure (June 9–10 and July 2) by counting the number of lures showing obvious beak (birds), teeth (micro-mammals), and mandibles (arthropods) marks per individual tree.

## 2.7 | Analyses

Irrigation being applied at the block level (i.e., whole plot), the ORPHEE experiment is a split-plot experiment which requires adapting the calculation of degrees of freedom and mean sum of squares of residuals. This was achieved using linear mixed effect models (LMM), with Block and Block × Irrigation as random factors (1|Block:Irrigation in R syntax). Herbivory and leaf trait data were analyzed at the tree level. Nonindependence of individual oaks within plots was accounted for by defining plot as a random factor, nested within blocks (Schielzeth & Nakagawa, 2013). For each test, we first built a full model including all fixed effects and their interactions (see details for each response variable below). The full model was then simplified by sequentially dropping nonsignificant terms, starting with highest-order interactions. Model parameters were finally estimated on the most parsimonious model obtained using restricted maximum likelihood (REML).

## 2.8 | Leaf traits

Leaf traits were measured on all trees, insect damage might trigger induced defenses modifying leaf traits. To test the direct, independent effect of oak neighbor identity (here after "*Neighbors*") on leaf traits, we limited the analysis of tree diversity on leaf traits to insecticide sprayed trees. The insecticide treatment reduced herbivory by ca 33% in spring (Kruskal-Wallis test: $\chi^2 = 40, p < .0001$) and 45% in summer (Kruskal-Wallis test: $\chi^2 = 50, p < .0001$). Traits were analyzed individually, using the same LMM modeling approach, with *Irrigation*, *Neighbors*, and their interaction as fixed effects.

Then, we tested whether oak leaf traits could account for variability in herbivory. We conducted a principal component analysis (PCA) on traits measured on unsprayed trees and extracted coordinates on the

two-first PCA axes (PC1 and PC2, respectively) that together explained 88% of variability in leaf traits. We used LMM (see above) to test whether PC1 and PC2 were explained by *Irrigation*, *Neighbors*, or their interaction.

## 2.9 | Insect herbivory

Only unsprayed trees were included in this set of analyses. We first tested the effect of oak neighbors and irrigation on insect herbivory using the same approach as for traits. Fixed effects were *Season*, *Neighbors*, and *Irrigation*, and all interactions. Because herbivory was assessed twice on the same individual oaks, we added tree identity as a random factor, nested within *Block* and *Plot* factors, to account for repeated measurements.

We then tested whether the effect of neighbors could be accounted for through changes in leaf traits measured on the same trees or predation pressure at the plot scale (i.e., proportion of attacked caterpillars pooled across the two surveys). *Traits and Predation* were included into the model as fixed effects, together with *Irrigation* and *Neighbor* factors and *Neighbors × Irrigation* interactions. Because trait values may have changed during the course of the season, and because predation was only assessed in early summer, we tested the effects of SLA, LDMC, leaf toughness, thickness, and predation on early season herbivory, and the effect of C:N on late season herbivory, separately. Given that both leaf traits and predation rates might have been influenced by irrigation and neighbor identity, this approach allows testing whether some residual variance can still be explained by plot-level factors (i.e., *Irrigation* and *Neighbors*) when accounting for traits and predation, and conversely. We expected that, should herbivory be better explained by leaf traits or predation than by irrigation or neighbor identity, significant effects of any of these factors (or their interaction) would become nonsignificant once traits and predation are accounted for. Conversely, should irrigation or neighbor identity effects remain significant after including leaf traits and predation, this would suggest that irrigation and neighbor identity influence herbivory through other unmeasured plant traits or unknown natural enemies (e.g., parasitoids).

Finally, we summarized the univariate effects of all individual traits using coordinates on the two PCA axes and followed the same approach to model insect herbivory as a function of PC1, PC2, *Irrigation*, *Neighbors*, PC1 × PC2, and *Irrigation × Neighbors* interactions. We then applied the same simplification procedure as explained above.

## 2.10 | Insectivory

We analyzed predation rates using generalized LMM (GLMM), with binomial error and logit link. The response variable was the number of attacked versus nonattacked dummy caterpillars per tree. We modeled separately the response of total predation (i.e., all predators combined, including birds, small rodents, and arthropods) and avian predation only to neighbor identity and irrigation.

For every model, conditions of application were visually checked and the response variable was log-transformed whenever necessary. All analyses were conducted in R (R Core Team, 2016) using the *ImerT-est* and *ade4* libraries.
# 3 | RESULTS

## 3.1 | Effects of irrigation and neighbors on water stress in oaks

Predawn leaf water potential of oaks ($\psi_w$) was significantly lower ($F = 35.13, p < .0001$) in blocks receiving only rainfall (-0.37 MPa) than in irrigated blocks (-0.19 MPa), confirming the effectiveness of the irrigation treatment (Fig. S1). There was a significant effect of plot composition on $\psi_w$ in oaks ($F = 4.14, p = .017$), with lower $\psi_w$ in nonirrigated plots containing pines (Fig. S1). Consistently, $\delta^{13}C$ was greater in the mixture with pine as compared to the two other mixtures and monocultures, thus confirming the greater water stress of oaks in two- and three-species mixtures including pine (Figure 2).

## 3.2 | Effects of irrigation and oak neighbor identity on insect herbivory

On average (±SE), insect herbivory on unsprayed trees was twice as high in summer compared to spring (16.7 ± 1.0 in summer and 8.7 ± 0.6% leaf area damaged in spring, respectively, $F_{(1,92.4)} = 78.72, p < .001$, Figure 1a).

Insect herbivory did not vary with the identity of oak neighbors ($F_{(3, 85.1)} = 0.60, p = .617$, Figure 1b), regardless of the season (Neighbors × Season: $F_{(3,89.8)} = 1.58, p = .201$).

The effect of irrigation was not significant ($F_{(1,6.0)} = 0.01, p = .913$), neither separately nor in interaction with the season ($F_{(1,88.8)} = 0.50, p = .479$) nor with oak neighbor identity ($F_{(3,81.9)} = 0.72, p = .543$). The simplified model only retained the season as predictor and explained 21% of variability in herbivory ($R^2_m = 0.21, R^2_c = 0.49$).

## 3.3 | Effects of irrigation and neighbor identity on oak leaf traits

In sprayed oak trees, SLA, leaf toughness, and leaf thickness significantly varied with oak neighbor identity (Figure 2), but none was influenced by irrigation or *Irrigation × Neighbor* interaction (Table 1). SLA was higher in plots where oak was associated with pine (two- and three-species mixtures). Toughness was lower in plots where pine was

present (two- and three-species mixtures). Thickness was higher in monocultures than in mixed plots. There was no significant effect of oak neighbor identity on LDMC or C:N.

The first and second axes of the principal component analysis explained 64.3% (PC1) and 23.5% (PC2) of variability in leaf traits, respectively. PC1 was driven by SLA, thickness, and toughness, positive values being associated with large, thin, and soft leaves (Figure 3a). PC2 was driven by water content, positive values being associated with high LDMC, that is, low water content (Figure 3a). There was a significant effect of oak neighbors on PC1, which was independent of irrigation (Table 1). Oak leaf traits were not significantly different between monocultures and oak-birch mixtures, but with two- and three-species mixtures containing pine (Figure 3b). There was no effect of neighbors or irrigation on PC2.

## 3.4 | Effects of irrigation and oak neighbors on predation

Overall, the predation rate of dummy caterpillars averaged across predator identity, plots, and dates was 12.7 ± 1.7%. Among the 110 lures displaying obvious predator marks (of 1,359 observations), 70.9% were attributed to birds (Figure 4). Attacks by small rodents (16.4%) and arthropods (12.7%) were too scarce and unevenly distributed to be analyzed separately. Insectivory was thus analyzed first for all predators and then for birds only. Overall, total insectivory was higher in irrigated than in rainfall only plots ($\chi^2 = 4.62, p = .032$, Figure 4). Differences in predation rates among plots with different tree species composition were not significant ($\chi^2 = 1.89, p = .596$, Figure 4). There was no difference between the two surveys ($\chi^2 = 0.12, p = .728$), and results were identical when only bird insectivory was considered. There was no correlation between the proportion of predated dummy caterpillars and the percentage of leaf area removed by herbivores, neither at the tree (Pearson's *r* = -0.04, *p* = .814) nor at the plot scale (Pearson's *r* = 0.13, *p* = .100).

## 3.5 | Effects of oak leaf traits on insect herbivory

When considered individually, none of the leaf traits measured in early season (i.e., SLA, thickness, toughness, LDMC) on unsprayed oak trees

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

**FIGURE 1** Effects of season (a) and tree species composition (b) on herbivory (% Leaf Area Removed—%LAR). Rainfall only and irrigated plots are combined. Boxes indicate the lower and upper quartiles. Thick horizontal lines and dots represent the median and the mean, respectively. Different letters above boxes indicate significant differences. *n* indicates sample size. Qr, *Quercus robur*; QrBp, *Q. robur* + *Betula pendula*; QrPp, *Q. robur* + *Pinus pinaster*; QrBpPp, *Q. robur* + *B. pendula* + *P. pinaster*
<figure><img src="image_4.png" /><figcaption>**FIGURE 2** Effects of tree species composition on oak leaf traits. (a) Specific leaf area, (b) leaf dry matter content, (c) leaf toughness, (d) leaf thickness, (e) leaf δ¹³C (x-1), (f) Leaf C:N. Data were summarized across all blocks, regardless of irrigation. Traits were compared in trees sprayed with insecticide such that differences in leaf traits are independent of any herbivory effect. Boxes indicate the lower and upper quartiles. Thick horizontal lines and dots represent the median and the mean, respectively. Same letters above bars indicate nonsignificant differences. Number within brackets below boxes indicates sample size. Qr, Quercus robur; QrBp, Q. robur + Betula pendula; QrPp, Q. robur + Pinus pinaster; QrBpPp, Q. robur + B. pendula + P. pinaster</figcaption></figure>

had a significant effect on insect herbivory (Table S1). Insectivory had no effect on herbivory neither (Table S1). While individually, the measured traits failed to explain variability in insect herbivory, there was a significant effect of PC1 × PC2 interaction on early season herbivory ($F_{(1,57,6)} = 5.30, p = .025$). Although main effects of PC1 and PC2 were

not significant, model parameter estimates were positive for PC1 (0.08 ± 0.09) and negative for PC2 (-0.11 ± .16). Furthermore, the significant interaction between PC1 and PC2 had a positive parameter estimate (model parameter estimate ± SE: 0.20 ± 0.09) which indicates that herbivory tended to increase more along PC1 axis (i.e., with larger,
<table><caption>**TABLE 1** Summary of LMM testing the effects of oak neighbor identity and block irrigation on oak leaf traits</caption><thead><tr><th>Trait</th><th>Predictors</th><th>F value (df)</th><th>p value</th><th>R<sup>2</sup>m (R<sup>2</sup>c)</th></tr></thead><tbody><tr><td rowspan="3">SLA</td><td>Neighbors</td><td>4.20 (3, 20.73)</td><td><strong>.018</strong></td><td>0.12 (0.29)</td></tr><tr><td>Irrigation</td><td>0.01 (1, 7.36)</td><td>.91</td><td></td></tr><tr><td>Neighbors × irrigation</td><td>1.37 (3, 18.04)</td><td>.283</td><td></td></tr><tr><td rowspan="3">LDMC</td><td>Neighbors</td><td>1.77 (3, 20.26)</td><td>.185</td><td>- (0.30)</td></tr><tr><td>Irrigation</td><td>0.01 (1, 6.49)</td><td>.916</td><td></td></tr><tr><td>Neighbors × irrigation</td><td>2.81(3, 17.21)</td><td>.07</td><td></td></tr><tr><td rowspan="3">Toughness</td><td>Neighbors</td><td>5.08 (3, 81.66)</td><td><strong>.003</strong></td><td>0.05 (0.70)</td></tr><tr><td>Irrigation</td><td>0.19 (1, 8.00)</td><td>.672</td><td></td></tr><tr><td>Neighbors × irrigation</td><td>0.45 (3, 78.50)</td><td>.715</td><td></td></tr><tr><td rowspan="3">Thickness</td><td>Neighbors</td><td>13.21 (3, 81.91)</td><td><strong>&lt;.0001</strong></td><td>0.22 (0.56)</td></tr><tr><td>Irrigation</td><td>0.02 (1, 7.44)</td><td>.901</td><td></td></tr><tr><td>Neighbors × irrigation</td><td>0.46 (3, 78.63)</td><td>.711</td><td></td></tr><tr><td rowspan="3">δ13C</td><td>Neighbors</td><td>14.97 (3, 18.71)</td><td><strong>&lt;.0001</strong></td><td>0.35 (0.62)</td></tr><tr><td>Irrigation</td><td>4.65 (1, 5.93)</td><td>.075</td><td></td></tr><tr><td>Neighbors × Irrigation</td><td>0.64 (3, 15.46)</td><td>.6</td><td></td></tr><tr><td rowspan="3">C:N</td><td>Neighbors</td><td>0.49 (3, 19.16)</td><td>.695</td><td>- (0.34)</td></tr><tr><td>Irrigation</td><td>0.02 (1, 6.18)</td><td>.884</td><td></td></tr><tr><td>Neighbors × irrigation</td><td>0.65 (3, 17.04)</td><td>.591</td><td></td></tr><tr><td rowspan="3">PC1</td><td>Neighbors</td><td>13.17 (3, 19.32)</td><td><strong>.0001</strong></td><td>0.34 (0.65)</td></tr><tr><td>Irrigation</td><td>0.11 (1, 5.86)</td><td>.748</td><td></td></tr><tr><td>Neighbors × irrigation</td><td>1.49 (3,15.30)</td><td>.256</td><td></td></tr><tr><td rowspan="3">PC2</td><td>Neighbors</td><td>0.90 (3, 57.80)</td><td>.0447</td><td>- (0.38)</td></tr><tr><td>Irrigation</td><td>0.07 (1, 6.00)</td><td>.794</td><td></td></tr><tr><td>Neighbors × irrigation</td><td>0.07 (3, 54.70)</td><td>.978</td><td></td></tr></tbody></table>

Marginal ($R^2_m$) and conditional ($R^2_c$, $R^2$) were calculated for the simplified model. They correspond to the variance explained by fixed and fixed plus random factors, respectively. Bold values indicate significant effects ($P < 0.05$).

<figure><img src="image_5.png" /><figcaption>**FIGURE 3** Principal Components Analysis of leaf traits. (a) Correlation circle showing correlations among leaf traits (SLA: specific leaf area; LDMC: leaf dry matter content) and between leaf traits and PCA axes. (b-c) Projections of individual trees on PCA axes according to (b) tree species composition and (c) irrigation treatment. Qr, Quercus robur; QrBp, Q. robur + Betula pendula; QrPp, Q. robur + Pinus pinaster; QrBpPp, Q. robur + B. pendula + P. pinaster</figcaption></figure>

thinner, and smoother leaves) when PC2 values were high (i.e., leaves with a low water content) and conversely (Figure 5a). These results indicate that altogether, leaf trait variation can significantly contribute to explain insect herbivory, although the final model explained only a

small fraction (8%) of variability in early season herbivory ($R^2_m = 0.08$, $R^2_c = 0.25$).

In late season, insect herbivory significantly decreased with the C:N ratio of oak leaves ($F_{(1,66.5)} = 4.53$, $p = .037$, Figure 5b), irrespective
<figure><img src="image_6.png" /><figcaption>**FIGURE 4** Effects of block irrigation (a) and tree species composition (b) on predation (i.e., attacks of dummy caterpillars). (*) indicates significant difference between treatments. Qr, *Quercus robur*; QrBp, *Q. robur* + *Betula pendula*; QrPp, *Q. robur* + *Pinus pinaster*; QrBpPp, *Q. robur* + *B. pendula* + *P. pinaster*</figcaption></figure>

of irrigation (Irrigation: $F_{(1,29,4)} = 1.38, p = .249$; Irrigation × C:N: $F_{(1,62,6)} = 0.70, p = .405$). C:N ratio, however, only explained 6% of variability in late season herbivory ($R_m^2 = 0.06, R_c^2 = 0.19$).

Underwood, 2014; Kos et al., 2015). Our results partially support this hypothesis.

## 4 | DISCUSSION

Contrary to our expectations (i.e., associational effects in presence of heterospecific neighbors), insect herbivory on oak was not influenced by oak neighbor identity. Insect herbivory was, however, partially influenced by leaf traits, which responded to the identity of oak neighbors. Despite higher predation in irrigated plots, water stress had no direct or indirect effect on insect herbivory.

### 4.1 | Oak neighbor identity had weak indirect trait-mediated effects on insect herbivory

Contrary to previous studies reporting an effect of tree composition or diversity on insect herbivores associated with oaks (Alalouni, Brandl, Auge, & Schädler, 2014; Castagneyrol et al., 2013; Setiawan et al., 2014), we observed no difference in insect herbivory among different neighbor identities. Two main mechanisms have been proposed to explain why herbivory should vary with the composition of forests: (1) The probability of a tree being colonized by herbivores (i.e., its apparency, frequency, and concentration) may depend on the diversity and identity of its neighbors (Barbosa et al., 2009; Castagneyrol et al., 2013; Hambäck, Inouye, Andersson, & Underwood, 2014; Setiawan et al., 2016) and (2) the abundance and activity of herbivore enemies may change with forest composition and structure (Muiruri et al., 2016; Riihimäki et al., 2004; Schuldt et al., 2011). Our results support neither of these. Oak apparency, frequency, and concentration (i.e., oak accessibility) were consistently low in all mixtures as compared to monocultures, and yet, they did not receive lower herbivory in the presence of taller birches or pines (Castagneyrol et al., 2013). It is unlikely that the effects of greater apparency on herbivore recruitment were compensated by higher predation rates in oak monocultures, as predation remained unaffected by oak neighbor identity.

Here, we additionally investigate a third, complementary hypothesis, where (3) oak neighbor identity could modify leaf traits involved in tree-herbivore interactions (Halpern, Bednar, Chisholm, &

#### 4.1.1 | Oak leaf traits varied with oak neighbor identity

The main differences in oak leaf traits were detected between oak monocultures and mixtures (Figure 2), in particular between monocultures and oak-pine mixtures. In the ORPHEE experiment, in 2015, oaks (mean height ± SE: 110.0 ± 0.8 cm) were much smaller than birches (510.0 ± 2.0 cm) and pines (563.0 ± 1.5 cm) (Damien et al., 2016). The total amount of light received by oaks was thus likely reduced in mixed plots where oaks were dominated by pines and birches. Although we did not quantify light interception by dominant species, this is consistent with the observed differences in SLA, toughness, δ¹³C, and thickness among treatments (Figure 2). Higher SLA, more negative δ¹³C values, and lower leaf thickness and toughness, as observed in mixtures, are indeed typical response patterns to more shaded conditions (Pérez-Harguindeguy et al., 2013). Furthermore, differences in leaf traits between oak monocultures and oak-pine mixtures were clearly shown on the first PCA axes (Figure 3) while leaf traits in oak-birch mixtures were more similar to leaf traits in monocultures. Given the canopy shape and foliar content of birch and pine trees, one can hypothesize that even though the two species have similar mean heights, solar radiation is more intercepted by pine trees (which are more opaque). Light could thus be one factor controlling leaf traits of oaks growing below the canopy of taller neighbors, but it was probably not the only one. Furthermore, our results do not support the hypothesis that variability in leaf traits among plots with different specific compositions was driven by differential response to water stress (Walter et al., 2011; Forey et al., 2016) as irrigation had no effect on leaf traits, neither in univariate nor in multivariate analyses. It must be acknowledged that we applied irrigation only few weeks before traits were measured, and it cannot be excluded that stronger effects will emerge in the future. Here, we can only speculate on the mechanisms responsible for such differences among mixtures, but it is likely that litter composition or understory vegetation differentially also affected oak traits in monocultures and oak-birch or oak-pine mixtures (Nickmans, Verheyen, Guiz, Jonard, & Ponette, 2015).
<figure><img src="image_7.png" /><figcaption>**FIGURE 5** Effect of leaf traits on herbivory. (a) Interactive effects of PC1 and PC2 on herbivory. Gray scale colors represent herbivory (% leaf area removed—%LAR) predicted by GLMM retaining only PC1, PC2, and their interaction as predictors. Individual trees are plotted as in Figure 3b. Herbivory increased strongly with PC1 for high values of PC2, but decreased slightly with PC1 for low values on PC2. Herbivory increased along PC2 axis when PC1 values were high, while it decreased slightly for low values of PC1. (b) Effects of C:N ratio on late season herbivory. Dashed and solid regression lines are for irrigated and rainfall only blocks, respectively. Shaded area represent model SE. Qr, *Quercus robur*; QrBp, *Q. robur* + *Betula pendula*; QrPp, *Q. robur* + *Pinus pinaster*; QrBpPp, *Q. robur* + *B. pendula* + *P. pinaster*</figcaption></figure>

### 4.1.2 | Leaf traits predicted only a limited fraction of insect herbivory

As oak neighbor identity affected variability in oak leaf traits, we expected that particular traits might be the functional link between oak

neighbor identity and herbivory. This would have been in accordance with recent work showing trait-mediated effects of heterospecific plant neighbors on insect herbivores and herbivory (Halpern et al., 2014; Ohgushi & Hambäck, 2015; Kos et al., 2015; Kostenko et al. 2017). Our results weakly support this prediction. Indeed, only leaf C:N ratio had an effect on herbivory (Figure 5), yet this trait did not vary with oak neighbor identity (Figure 2). Individually, none of the other traits measured in early season had a significant effect on herbivory. However, traits may have had a combined effect on herbivory as suggested by the multivariate analysis. This effect was, however, complex and explained only 8% of total variability in herbivory.

Although our results partially support the idea that tree species composition had indirect trait-mediated effects on overall insect herbivory, they should be considered with caution. Previous studies reported mixed evidence for such a relationship between plant diversity, leaf traits, and herbivory (Moreira et al., 2014; Abdala-Roberts, 2016; Kostenko et al. 2017). Here, because we measured total chewing damage, we cannot dismiss the fact that the overall pattern obscures species-specific variability in herbivore response to leaf traits (Heiermann & Schütz, 2008; Plath et al., 2011), which would deserve more investigation.

### 4.2 | Predation barely changed with oak neighbor identity

The enemies' hypothesis posits that more species-rich forest plots should shelter a greater diversity or abundance of predators than monocultures, enhancing potential biological control of herbivores. In particular, mixing conifer and deciduous tree species at both forest stand and landscape scales is expected to increase insectivorous predator abundance (Charbonnier et al., 2016; Oxbrough, García-Tejero, Spence, & O'Halloran, 2016), with cascading effects on pest control (Felton et al., 2016; Giffard, Barbaro, Jactel, & Corcket, 2013). Although the enemies' hypothesis has received some empirical support in grasslands and agro-ecosystems (Andow, 1991; Langellotto & Denno, 2004; Letourneau et al., 2011; Scherber et al., 2010), its validity for forest is still debated (Letourneau, Jedlicka, Bothwell, & Moreno, 2009; Muiruri et al., 2016; Riihimäki et al., 2004; Schuldt et al., 2011). In addition, it must be noticed that the use of dummy caterpillars underestimates the importance of top-down processes as it does not capture the effect of other natural enemies such as spiders or parasitoids which are known to, respectively, respond to structural and chemical complexity of their habitat (Kostenko et al., 2015; Langellotto & Denno, 2004).

### 4.3 | Irrigation had no direct effect on herbivory but indirectly changed predation

Although predawn leaf water potentials ($\psi_w$) and carbon isotope values ($\delta^{13}C$) confirmed that irrigation alleviated water stress, neither leaf traits nor insect herbivory was influenced by irrigation. However, water stress varied among oaks with different neighbors as indicated by the more negative $\psi_w$ and less negative $\delta^{13}C$ in oak-pines mixtures, both indicating
stronger water stress. This could be due to pines being much taller than oaks and exploiting a greater amount of water. Such an effect of inter-specific competition for water may have been less pronounced with birches that lose leaves during summer, while pines kept their needles.

The lack of effect of irrigation on leaf traits and herbivory was unexpected, as the influence of drought on leaf functional traits has been known for decades (Chaves, Maroco, & Pereira, 2003). However, $\psi_w$ in the rainfall only plots did not reach very negative values during the 2015 growing season. They remained above -0.5 MPa and the difference in $\psi_w$ suggests that during the 2015 growing season, unirrigated trees were subject to low water stress. Yet, the response of herbivores may vary nonlinearly with stress intensity (Gutbrodt, Dorn, & Mody, 2011; Jactel et al., 2012). Therefore, we cannot exclude that the lack of effect of irrigation on leaf traits and herbivory was due to low drought intensity.

Moreover, irrigation started only a few weeks before trait measurements. This period was enough to influence soil water availability in the rainfall only versus irrigation plots (as seen in $\psi_w$ values, Fig. S1), but it may not have been long enough, and the soil water deficit may not have been severe enough, to significantly influence leaf morphological and physiological properties. This later assumption is also supported by the fact that many temperate *Quercus* species have preformed growth of shoots in buds (Fontaine et al., 2000): The organogenesis period for the first growth unit ranges approximately from August to October of the previous year (Fontaine et al., 2000), whereas the elongation periods range from March to September of the current year. The influence of changing environmental conditions in spring on *Quercus* leaf traits may have thus been lower than for species in which leaf organogenesis is directly occurring in spring. This explanation is also consistent with the fact that changes in herbivory-related leaf traits involving cell division and differentiation (such as those we measured: SLA, LDMC, toughness) take longer than changes in leaf chemistry involving the production of secondary metabolites (Herms & Mattson, 1992; Moreira et al., 2014) and may contribute to explaining why we did not detect any significant effect of irrigation on leaf structural traits.

Predation was the only response variable that was influenced by drought, with greater predation in irrigated plots, especially by rodents and arthropods. This may result from understorey vegetation being denser and more stratified in irrigated plots (*personal observations*), potentially resulting in larger niche opportunities, higher prey availability, or lower top-down predation risk by apex predators (Figure 4a). Similar effects of structural complexity have been observed by Poch and Simonetti (2013) who detected higher predation rates in forest plantations with a more developed understorey. Further investigation will be needed to tease apart the potential direct effects of irrigation on herbivory and indirect effects mediated by predation.

5 | CONCLUSION

We did not find strong support to some of the most commonly proposed hypotheses to explain herbivory patterns, notably the apparenty hypothesis. We found no clear support for the enemies'

hypothesis. We found only partial support for the emerging idea that associational effects could be mediated by change in leaf traits relevant to insect herbivory. Although we could only speculate on the mechanisms linking tree species composition, variability in leaf traits, and herbivory, our results suggest that indirect trait-mediated effects (including secondary metabolites) will be a key aspect to consider in the future studies addressing associational effects of plant diversity on higher trophic levels.

ACKNOWLEDGMENTS

The research in ORPHEE was funded by the GIP-ECOFOR program from the French Ministry of Ecology, under the project BIOPICC ECOFOR-2014-15. The work was supported also by a grant overseen by the French National Research Agency (ANR) as part of the "Investissements d'Avenir" program (ANR-11-LABX-0002-01, Laboratory of Excellence ARBRE). We thank Rosalie Lefeu, Inge van Halder, Fabrice Vétillard, Bruno Garnier, Bastien Gravellier, and Mélanie Thiery for their help with field work. We thank Nicolas Angeli and Christian Hossann for the isotope and elementary analyses (INRA Nancy Technical Facility of Functional Ecology-OC 081). Thanks to the experiment unit UE Forêt-Pierroton (INRA), and especially to Bernard Issenhuth for yearly management of the ORPHEE experiment.

CONFLICT OF INTEREST

None declared.

REFERENCES

Abdala-Roberts, L., Hernandez-Cumplido, J., Chel-Guerrero, L., Betancur-Ancona, D., Benrey, B., & Moreira, X. (2016). Effects of plant intraspecific diversity across three trophic levels: Underlying mechanisms and plant traits. *American Journal of Botany*, *103*, 1810-1818.

Alalouni, U., Brandl, R., Auge, H., & Schädler, M. (2014). Does insect herbivory on oak depend on the diversity of tree stands? *Basic and Applied Ecology*, *15*, 685-692. doi:10.1016/j.baae.2014.08.013

Allan, E., Weisser, W. W., Fischer, M., Schulze, E.-D., Weigelt, A., Roscher, C., ... Schmid, B. (2013). A comparison of the strength of biodiversity effects across multiple functions. *Oecologia*, *173*, 223-237.

Andersson, P., Löfstedt, C., & Hambäck, P. A. (2013). Insect density-plant density relationships: A modified view of insect responses to resource concentrations. *Oecologia*, *173*, 1333-1344. doi:10.1007/s00442-013-2737-1

Andow, D. A. (1991). Vegetational diversity and arthropod population response. *Annual Review of Entomology*, *36*, 561-586.

Aslam, T. J., Johnson, S. N., & Karley, A. J. (2013). Plant-mediated effects of drought on aphid population structure and parasitoid attack. *Journal of Applied Entomology*, *137*, 136-145.

Atsatt, P. R., & O'Dowd, D. J. (1976). Plant defense guilds. *Science*, *193*, 24-29. doi:10.1126/science.193.4247.24

Barbaro, L., Giffard, B., Charbonnier, Y., van Halder, I., & Brockerhoff, E. G. (2014). Bird functional diversity enhances insectivory at forest edges: a transcontinental experiment. *Diversity and Distributions*, *20*, 149-159.

Barbosa, P., Hines, J., Kaplan, I., Martinson, H., Szczepaniec, A., & Szendrei, Z. (2009). Associational Resistance and Associational Susceptibility: Having Right or Wrong Neighbors. *The Annual Review of Ecology, Evolution, and Systematics*, *40*, 1-20.
Bauerfeind, S. S., & Fischer, K. (2013). Targeting the right trait: The relative suitability of a host plant depends on the herbivore trait considered and ambient temperature. *Basic and Applied Ecology, 14*, 555-564. doi:10.1016/j.baae.2013.08.010

Broad, S. T., Schellhorn, N. A., Lisson, S. N., & Mendham, N. J. (2008). Host location and oviposition of lepidopteran herbivores in diversified broccoli cropping systems. *Agricultural and Forest Entomology, 10*, 157-165. doi:10.1111/j.1461-9563.2008.00374.x

Castagneyrol, B., Giffard, B., Péré, C., & Jactel, H. (2013). Plant apparenty, an overlooked driver of associational resistance to insect herbivory. *Journal of Ecology, 101*, 418-429. doi:10.1111/1365-2745.12055

Castagneyrol, B., Jactel, H., Vacher, C., Brockerhoff, E. G., & Koricheva, J. (2014). Effects of plant phylogenetic diversity on herbivory depend on herbivore specialization. *Journal of Applied Ecology, 51*, 134-141.

Charbonnier, Y. M., Barbaro, L., Barnagaud, J.-Y., Ampoorter, E., Nezan, J., Verheyen, K., & Jactel, H. (2016). Bat and bird diversity along independent gradients of latitude and tree composition in European forests. *Oecologia, 182*(2), 529-537. doi:10.1007/s00442-016-3671-9

Chaves, M. M., Maroco, J. P., & Pereira, J. S. (2003). Understanding plant responses to drought—from genes to the whole plant. *Functional Plant Biology, 30*, 239. doi:10.1071/F02076

R Core Team (2016). *R: A language and environment for statistical computing*. Vienna, Austria: R Foundation for statistical computing.

Damien, M., Jactel, H., Meredieu, C., Régolini, M., van Halder, I., & Castagneyrol, B. (2016). Pest damage in mixed forests: Disentangling the effects of neighbor identity, host density and host apparenty at different spatial scales. *Forest Ecology and Management, 378*, 103-110.

Elton, C. S. (1958). *The ecology of invasions by animals and plants*. London: Methuen.

Felton, A., Nilsson, U., Sonesson, J., Felton, A. M., Roberge, J.-M., Ranius, T., ... Wallertz, K. (2016). Replacing monocultures with mixed-species stands: Ecosystem service implications of two production forest alternatives in Sweden. *Ambio, 45*, 124-139.

Floater, G. J., & Zalucki, M. P. (2000). Habitat structure and egg distributions in the processionary caterpillar *Ochrogaster lunifer*: Lessons for conservation and pest management. *Journal of Applied Ecology, 37*, 87-99. doi:10.1046/j.1365-2664.2000.00468.x

Fontaine, F., Chaar, H., Colin, F., Clément, C., Burrus, M., & Druelle, J.-L. (2000). Preformation and neoformation of growth units on 3-year-old seedlings of Quercus petraea. *Canadian Journal of Botany, 77*, 1623-1631.

Forey, E., Langlois, E., Lapa, G., Korboulewsky, N., Robson, T. M., & Aubert, M. (2016). Tree species richness induces strong intraspecific variability of beech (*Fagus sylvatica*) leaf traits and alleviates edaphic stress. *European Journal of Forest Research, 1*–11.

Forrester, D. I., Theiveyanathan, S., Collopy, J. J., & Marcar, N. E. (2010). Enhanced water use efficiency in a mixed *Eucalyptus globulus* and *Acacia mearnsii* plantation. *Forest Ecology and Management, 259*, 1761-1770. doi:10.1016/j.foreco.2009.07.036

Giffard, B., Barbaro, L., Jactel, H., & Corcket, E. (2013). Plant neighbours mediate bird predation effects on arthropod abundance and herbivory. *Ecological Entomology, 448*–455, doi:10.1111/een.12035

Grossiord, C., Granier, A., Gessler, A., Jucker, T., & Bonal, D. (2014). Does Drought Influence the Relationship between Biodiversity and Ecosystem Functioning in Boreal Forests? *Ecosystems, 17*, 394-404. doi:10.1007/s10021-013-9729-1

Gutbrodt, B., Dorn, S., & Mody, K. (2011). Drought stress affects constitutive but not induced herbivore resistance in apple plants. *Arthropod-Plant Interactions, 6*, 171-179.

Hahn, P. G., & Orrock, J. L. (2016). Neighbor palatability generates associational effects by altering herbivore foraging behavior. *Ecology, 97*, 2103-2111. doi:10.1002/ecy.1430

Halpern, S. L., Bednar, D., Chisholm, A., & Underwood, N. (2014). Plant-mediated effects of host plant density on a specialist herbivore of *Solanum carolinense*: Plant-mediated effects on Letourneau oviposition. *Ecological Entomology, 39*, 217-225. doi:10.1111/een.12088

Hambäck, P. A., Inouye, B. D., Andersson, P., & Underwood, N. (2014). Effects of plant neighborhoods on plant-herbivore interactions: Resource dilution and associational effects. *Ecology, 95*, 1370-1383. doi:10.1890/13-0793.1

Heiermann, J., & Schütz, S. (2008). The effect of the tree species ratio of European beech (*Fagus sylvatica* L.) and Norway spruce (*Picea abies* (L.) Karst.) on polyphagous and monophagous pest species - *Lymantria monacha* L. and *Calliteara pudibunda* L. (Lepidoptera: Lymantriidae) as an example. *Forest Ecology and Management, 255*, 1161-1166.

Herms, D. A., & Mattson, W. J. (1992). The dilemma of plants: To grow or defend. *The Quarterly Review of Biology, 67*, 283-335.

Huberty, A. F., & Denno, R. F. (2004). Plant water stress and its consequences for herbivorous insects: A new synthesis. *Ecology, 85*, 1383-1398. doi:10.1890/03-0352

Jactel, H., Birgersson, G., Andersson, S., & Schlyter, F. (2011). Non-host volatiles mediate associational resistance to the pine processionary moth. *Oecologia, 166*, 703-711.

Jactel, H., & Brockerhoff, E. G. (2007). Tree diversity reduces herbivory by forest insects. *Ecology Letters, 10*, 835-848. doi:10.1111/j.1461-0248.2007.01073.x

Jactel, H., Petit, J., Desprez-Loustau, M.-L., Delzon, S., Piou, D., Battisti, A., & Koricheva, J. (2012). Drought effects on damage by forest insects and pathogens: a meta-analysis. *Global Change Biology, 18*, 267-276.

Johnson, M. T. J., Bertrand, J. A., & Turcotte, M. M. (2016). Precision and accuracy in quantifying herbivory. *Ecological Entomology, 41*, 112-121. doi:10.1111/een.12280

Kos, M., Bukovinszky, T., Mulder, P. P. J., & Bezemer, T. M. (2015). Disentangling above- and belowground neighbor effects on the growth, chemistry, and arthropod community on a focal plant. *Ecology, 96*, 164-175. doi:10.1890/14-0563.1

Kostenko, O., Lammers, M., Grootemaat, S., Kroon, T., Harvey, J. A., Van Geem, M., & Bezemer, T. M. (2015). Effects of plant diversity and structural complexity on parasitoid behaviour in a field experiment. *Ecological Entomology, 40*, 748-758.

Kostenko, O., Mulder, P. P. J., Courbois, M., & Bezemer, T. M. (2017). Effects of Plant Diversity on the Concentration of Secondary Plant Metabolites and the Density of Arthropods on Focal Plants in the Field. *Journal of Ecology, 139*, 1-10. doi:10.1007/s00442-004-1497-3

Letourneau, D. K., Armbrecht, I., Rivera, B. S., Lerma, J. M., Carmona, E. J., Daza, M. C., ... Trujillo, A. R. (2011). Does plant diversity benefit agro-ecosystems? A synthetic review. *Ecological Applications, 21*, 9-21.

Letourneau, D. K., Jedlicka, J. A., Bothwell, S. G., & Moreno, C. R. (2009). Effects of natural enemy biodiversity on the suppression of arthropod herbivores in terrestrial ecosystems. *Annual Review of Ecology Evolution and Systematics, 40*, 573-592. doi:10.1146/annurev. ecolsys.110308.120320

Logan, J. A., Régnière, J., & Powell, J. A. (2003). Assessing the impacts of global warming on forest pest dynamics. *Frontiers in Ecology and the Environment, 1*, 130-137. doi:10.1890/1540-9295(2003)001[0130:ATIOGW]2.0.CO;2

Loranger, J., Meyer, S. T., Shipley, B., Kattge, J., Loranger, H., Roscher, C., ... Weisser, W. W. (2013). Predicting invertebrate herbivory from plant traits: Polycultures show strong nonadditive effects. *Ecology, 94*, 1499-1509.

Mäntylä, E., Alessio, G. A., Blande, J. D., Heijjari, J., Holopainen, J. K., Laaksonen, T., ... Klemola, T. (2008). From Plants to Birds: Higher Avian Predation Rates in Trees Responding to Insect Herbivory (ed S Rands). *PLoS One, 3*, e2832.

Moreira, X., Abdala-Roberts, L., Parra-Tabla, V., & Mooney, K. A. (2014). Positive effects of plant genotypic and species diversity on anti-herbivore defenses in a tropical tree species. *PLoS ONE, 9*, e105438. doi:10.1371/journal.pone.0105438
Moreira, X., Abdala-Roberts, L., Rasmann, S., Castagneyrol, B., & Mooney, K. A. (2016). Plant diversity effects on insect herbivores and their natural enemies: current thinking, recent findings, and future directions. *Current Opinion in Insect Science*, **14**, 1–7.

Moreira, X., Mooney, K. A., Zas, R., & Sampedro, L. (2012). Bottom-up effects of host-plant species diversity and top-down effects of ants interactively increase plant performance. *Proceedings of the Royal Society B: Biological Sciences*, **279**, 4464–4472.

Mraja, A., Unsicker, S. B., Reichelt, M., Gershenzon, J., & Roscher, C. (2011). Plant community diversity influences allocation to direct chemical defence in *Plantago lanceolata*. *PLoS ONE*, **6**, e28055.

Muiruri, E. W., Rainio, K., & Koricheva, J. (2016). Do birds see the forest for the trees? Scale-dependent effects of tree diversity on avian predation of artificial larvae. *Oecologia*, **180**, 619–630.

Nickmans, H., Verheyen, K., Guiz, J., Jonard, M., & Ponette, Q. (2015). Effects of neighbourhood identity and diversity on the foliar nutrition of sessile oak and beech. *Forest Ecology and Management*, **335**, 108–117.

Ohgushi, T., & Hambäck, P. A. (2015). Toward a spatial perspective of plant-based indirect interaction webs: Scaling up trait-mediated indirect interactions. *Perspectives in Plant Ecology, Evolution and Systematics*, **17**, 500–509. doi:10.1016/j.ppees.2015.09.006

Oxbrough, A., García-Tejero, S., Spence, J., & O'Halloran, J. (2016). Can mixed stands of native and non-native tree species enhance diversity of epigaeic arthropods in plantation forests? *Forest Ecology and Management*, **367**, 21–29. doi:10.1016/j.foreco.2016.02.023

Pearse, I. S. (2011). The role of leaf defensive traits in oaks on the preference and performance of a polyphagous herbivore, *Orygia vetusta*. *Ecological Entomology*, **36**, 635–642. doi:10.1111/j.1365-2311.2011.01308.x

Péré, C., Jactel, H., & Kenis, M. (2013). Response of insect parasitism to elevation depends on host and parasitoid life-history strategies. *Biology Letters*, **9**, 20130028. doi:10.1098/rsbl.2013.0028

Pérez-Harguindeguy, N., Díaz, S., Garnier, E., Lavorel, S., Poorter, H., Jaureguiberry, P., ... Cornelissen, J. H. C. (2013). New handbook for standardised measurement of plant functional traits worldwide. *Australian Journal of Botany*, **61**, 167–234.

Plath, M., Dorn, S., Riedel, J., Barrios, H., & Mody, K. (2011). Associational resistance and associational susceptibility: specialist herbivores show contrasting responses to tree stand diversification. *Oecologia*, **169**, 477–487.

Poch, T. J., & Simonetti, J. A. (2013). Ecosystem services in human-dominated landscapes: Insectivory in agroforestry systems. *Agroforestry Systems*, **87**, 871–879. doi:10.1007/s10457-013-9603-3

Riihimäki, J., Kaitaniemi, P., Koricheva, J., & Vehviläinen, H. (2004). Testing the enemies hypothesis in forest stands: The important role of tree species composition. *Oecologia*, **142**, 90–97. doi:10.1007/s00442-004-1696-y

Rodríguez-Castañeda, G. (2013). The world and its shades of green: A meta-analysis on trophic cascades across temperature and precipitation gradients. *Global Ecology and Biogeography*, **22**, 118–130. doi:10.1111/j.1466-8238.2012.00795.x

Root, R. B. (1973). Organization of a plant-arthropod association in simple and diverse habitats: The fauna of collards (*Brassica Oleracea*). *Ecological Monographs*, **43**, 95. doi:10.2307/1942161

Schädler, M., Brandl, R., & Haase, J. (2007). Antagonistic interactions between plant competition and insect herbivory. *Ecology*, **88**, 1490–1498. doi:10.1890/06-0647

Scherber, C., Eisenhauer, N., Weisser, W. W., Schmid, B., Voigt, W., Fischer, M., ... Tscharntke, T. (2010). Bottom-up effects of plant diversity on multitrophic interactions in a biodiversity experiment. *Nature*, **468**, 553–556.

Schielzeth, H., & Nakagawa, S. (2013). Nested by design: Model fitting and interpretation in a mixed model era. *Methods in Ecology and Evolution*, **4**, 14–24. doi:10.1111/j.2041-210x.2012.00251.x

Schuldt, A., Both, S., Bruelheide, H., Hardtle, W., Schmid, B., Zhou, H. Z., & Assmann, T. (2011). Predator Diversity and Abundance Provide Little Support for the Enemies Hypothesis in Forests of High Tree Diversity. *PLoS One*, **6**(7), e22905.

Setiawan, N. N., Vanhellemont, M., Baeten, L., Dillen, M., & Verheyen, K. (2014). The effects of local neighbourhood diversity on pest and disease damage of trees in a young experimental forest. *Forest Ecology and Management*, **334**, 1–9.

Setiawan, N. N., Vanhellemont, M., Baeten, L., Gobin, R., De Smedt, P., Proesmans, W., ... Verheyen, K. (2016). Does neighbourhood tree diversity affect the crown arthropod community in saplings? *Biodiversity and Conservation*, **25**, 169–185.

Singer, M. S., Lichter-Marck, I. H., Farkas, T. E., Aaron, E., Whitney, K. D., & Mooney, K. A. (2014). Herbivore diet breadth mediates the cascading effects of carnivores in food webs. *Proceedings of the National Academy of Sciences*, **111**, 9521–9526.

Staley, J. T., Mortimer, S. R., Masters, G. J., Morecroft, M. D., Brown, V. K., & Taylor, M. E. (2006). Drought stress differentially affects leaf-mining species. *Ecological Entomology*, **31**, 460–469.

Straub, C. S., Simasek, N. P., Dohm, R., Gapinski, M. R., Aikens, E. O., & Nagy, C. (2014). Plant diversity increases herbivore movement and vulnerability to predation. *Basic and Applied Ecology*, **15**, 50–58.

Tahvanainen, J. O., & Root, R. B. (1972). The influence of vegetational diversity on the population ecology of a specialized herbivore, *Phyllotreta cruciferae* (Coleoptera: Chrysomelidae). *Oecologia*, **10**, 321–346. doi:10.1007/BF00345736

Underwood, N., Inouye, B. D., & Hambäck, P. A. (2014). A conceptual framework for associational effects: When do neighbors matter and how would we know? *Quarterly Review of Biology*, **89**, 1–19.

Walter, J., Hein, R., Auge, H., Beierkuhnlein, C., Löffler, S., Reifenrath, K., ... Jentsch, A. (2011). How do extreme drought and plant community composition affect host plant metabolites and herbivore performance? *Arthropod-Plant Interactions*, **6**, 15–25.

Weldegergis, B. T., Zhu, F., Poelman, E. H., & Dicke, M. (2015). Drought stress affects plant metabolites and herbivore preference but not host location by its parasitoids. *Oecologia*, **177**, 701–713.

White, T. C. (1974). A hypothesis to explain outbreaks of looper caterpillars, with special reference to populations of *Selidosema suavis* in a plantation of *Pinus radiata* in New Zealand. *Oecologia*, **16**, 279–301.

White, J. A., & Whitham, T. G. (2000). Associational susceptibility of cottonwood to a box elder herbivore. *Ecology*, **81**, 1795–1803.

Zhang, Q.-H., & Schlyter, F. (2004). Olfactory recognition and behavioural avoidance of angiosperm nonhost volatiles by conifer-inhabiting bark beetles. *Agricultural and Forest Entomology*, **6**, 1–20.

## SUPPORTING INFORMATION

Additional Supporting Information may be found online in the supporting information tab for this article.

**How to cite this article:** Castagneyrol B, Bonal D, Damien M, et al. Bottom-up and top-down effects of tree species diversity on leaf insect herbivory. *Ecol Evol.* 2017;7: 3520–3531. https://doi.org/10.1002/ece3.2950