# Benefits of jasmonate-dependent defenses against vertebrate herbivores in nature

Ricardo Ar Machado, Mark Mcclure, Maxime Hervé, Ian T Baldwin, Matthias Erb

► To cite this version:

Ricardo Ar Machado, Mark Mcclure, Maxime Hervé, Ian T Baldwin, Matthias Erb. Benefits of jasmonate-dependent defenses against vertebrate herbivores in nature. eLife, 2016, 5, 10.7554/eLife.13720. hal-02638800

HAL Id: hal-02638800

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

Submitted on 28 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 Attribution 4.0 International License
<figure><img src="image_2.png" /></figure>

# Benefits of jasmonate-dependent defenses against vertebrate herbivores in nature

Ricardo AR Machado¹,²,³, Mark McClure⁴, Maxime R Hervé²,⁵, Ian T Baldwin³, Matthias Erb¹,²*

¹Root-Herbivore Interactions Group, Department of Molecular Ecology, Max Planck Institute for Chemical Ecology, Jena, Germany; ²Institute of Plant Sciences, University of Bern, Bern, Switzerland; ³Department of Molecular Ecology, Max-Planck Institute for Chemical Ecology, Jena, Germany; ⁴School of the Environment, Washington State University, Washington, United States; ⁵Institut de Génétique, Environment et Protection des Plantes, Le Rheu, France

**Abstract** Endogenous jasmonates are important regulators of plant defenses. If and how they enable plants to maintain their reproductive output when facing community-level herbivory under natural conditions, however, remains unknown. We demonstrate that jasmonate-deficient *Nicotiana attenuata* plants suffer more damage by arthropod and vertebrate herbivores than jasmonate-producing plants in nature. However, only damage by vertebrate herbivores translates into a significant reduction in flower production. Vertebrate stem peeling has the strongest negative impact on plant flower production. Stems are defended by jasmonate-dependent nicotine, and the native cottontail rabbit *Sylvilagus nuttallii* avoids jasmonate-producing *N. attenuata* shoots because of their high levels of nicotine. Thus, endogenous jasmonates enable plants to resist different types of herbivores in nature, and jasmonate-dependent defenses are important for plants to maintain their reproductive potential when facing vertebrate herbivory. Ecological and evolutionary models on plant defense signaling should aim at integrating arthropod and vertebrate herbivory at the community level.

DOI: 10.7554/eLife.13720.001

## Introduction

In nature, plants are attacked by a multitude of herbivore species. By removing plant tissues and interrupting essential physiological processes, herbivory can negatively affect plant fitness, impose selection pressure, and thus drive the evolution of plant defenses (Züst et al., 2012; Agrawal et al., 2012; Huber et al., 2016a; 2016b). Herbivory does not necessarily reduce plant survival and reproduction, however (Bruelheide and Scheidel, 1999; Morris et al., 2007; Corbett et al., 2011). For example, minor insect herbivory that occurs early in the growing season can induce plant defenses, and thereby ensure enhanced protection during times of greater herbivore abundance (Agrawal, 1998; Kessler and Baldwin, 2004; Poelman et al., 2008). Likewise, herbivory that reduces competitor abundance may benefit a given individual (Agrawal et al., 2012). Accordingly, measuring the impact of herbivory on plant reproductive fitness in an ecologically relevant setting (e.g., one that includes the full suite of native herbivores), is essential for identifying the selective forces underlying the evolution of plant defenses.

Most plant defenses are regulated by jasmonates (Howe and Jander, 2008). Interrupting the biosynthesis of these hormones results in decreased constitutive and induced production of volatile and non-volatile secondary metabolites, defensive proteins and structural defenses (Thaler et al., 2002;

<footnote>*For correspondence: matthias.erb@ips.unibe.ch</footnote>

<footnote>Competing interest: See page 15</footnote>

<footnote>Funding: See page 15</footnote>

<footnote>Received: 10 December 2015</footnote>

<footnote>Accepted: 10 May 2016</footnote>

<footnote>Published: 29 June 2016</footnote>

<footnote>Reviewing editor: Merijn R Kant, University of Amsterdam, Netherlands</footnote>

<footnote>© Copyright Machado et al. This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.</footnote>
**eLife digest** Plants are attacked by many different herbivores, including insects and mammals, and often produce toxins in response to protect themselves. Toxin production is regulated by plant hormones called jasmonates. It is commonly assumed this ability helps plants to survive and reproduce in nature. However, proof that a plant's own jasmonates (also known as "endogenous jasmonates") can increase a plant's fitness in the wild is lacking, especially in the context of attack by mammals.

Machado et al. have now asked whether endogenous jasmonates increase the fitness of coyote tobacco plants that were under attack by herbivores in their natural habitats in Southwestern Utah. Plants that lacked jasmonates were attacked more strongly by various herbivores, yet unexpectedly only the damage by mammals – including gophers, deer and rabbits – caused the plants to produce fewer flowers. Since plants with more flowers tend to produce more offspring, the number of flowers is a measure of a plant's fitness. Damage by insects, which are often seen as major enemies of plants, did not result in a significant impact on the number of flowers.

Laboratory experiments then revealed that damaging plants in a similar way to mammalian herbivores strongly reduced the plants' fitness. However mimicking insect damage did not have such a large effect. Finally, feeding experiments with cottontail rabbits revealed that jasmonate-producing plants are protected by higher levels of the nicotine toxin, which can explain why these plants fare better when attacked by mammals in nature.

Jasmonates are well known to regulate plant defenses and provide protection against a wide variety of herbivores. However, these new findings show that this only translates into fitness benefits for the plants against a subset of herbivores. A major challenge in the future will be to study how diverse communities of herbivores shape the evolution of plant defense signaling. Including larger herbivores, like mammals, into such experiments will be challenging but necessary to understand how plants survive in nature.

DOI: 10.7554/eLife.13720.002

Li et al., 2004; Paschold et al., 2007; Howe and Jander, 2008; Zhou et al., 2009). Surprisingly, however, despite the availability of many different jasmonate deficient mutants (McConn et al., 1997; Thaler et al., 2002; Zhou et al., 2009; Kallenbach et al., 2012), fitness benefits of endogenous jasmonate production in herbivore-attacked plants in natural environments have not been demonstrated so far. Exogenous applications of jasmonic acid (JA) or methyl jasmonate (MeJA) are known to enhance plant defenses and fitness (Baldwin, 1998; Thaler, 1999; Heil et al., 2001), but these treatments remain difficult to interpret because they may not simulate dosages, localizations, and timings of endogenous jasmonate responses. Apart from inducing defenses, jasmonates can also reduce plant growth (Baldwin, 1998; Hanley, 1998; Redman et al., 2001; Zavala and Baldwin, 2006; Bodenhausen and Reymond, 2007; Machado et al., 2013), by interfering, for example, with the gibberellin and auxin signaling cascades (Onkokesung et al., 2010; Yang et al., 2012). This suggests that the evolution of jasmonate defense signaling is accompanied by growth suppression that may reduce fitness benefits under low herbivore pressure.

Our understanding of jasmonates in plant-herbivore interactions is also wanting because most research has focused on leaf-feeding arthropods (Mafli et al., 2012; Falk et al., 2014), and much less is known with respect to vertebrate herbivores, even though vertebrates are often the primary consumers in plant communities (Paige and Whitham, 1987; Hodgson and Illius, 1996). What we do know is that vertebrate herbivores can exert strong selective pressure on plants (Collins et al., 1998; Becerra, 2015) by influencing growth (Paige and Whitham, 1987; Bergman, 2002; Liu et al., 2012; Ishihama et al., 2014), structural defenses (Abrahamson, 1975; White, 1988; Young and Okello, 1998; Takada et al., 2001; Wilson and Kerley, 2003; Young et al., 2003; Kato et al., 2008), reproductive timing (Zamora et al., 2001) and mortality (Veblen et al., 1989; Gill, 1992; Vila and Guibal, 2001; Saint-Andrieux et al., 2009). We also know that vertebrates tend to avoid plants that are rich in secondary metabolites, including condensed tannins and phenolics (Cooper and Owen-Smith, 1985; Owen Smith, 1993; Furstenburg and van Hoven, 1994; O'reilly-Wapstra et al., 2004; Jansen et al., 2007; DeGabriel et al., 2009; Rosenthal and Berenbaum,
2012). Furthermore, it has been shown that silencing the production of the nervous toxin nicotine can engender increased leaf damage by vertebrate browsers in the field (Steppuhn et al., 2008). Despite these advances, it remains unclear to what extent endogenous jasmonates actually help plants to maintain their fitness when facing herbivore communities that include vertebrates.

To assess the importance of jasmonates in protecting plants from native arthropod and verte-
brate herbivores, we studied three experimental *N. attenuata* populations in their native environ-
ment in the Great Basin Desert (United States). Each population consisted of a mix of jasmonate
deficient and wild type plants. For each population, we characterized the damage that was caused
by vertebrates and arthropods and then correlated damage patterns with plant flower production as
a measure of the plant's reproductive potential. We then assessed the impact of simulated herbivory
on jasmonate-dependent flower production and defoliation tolerance under glasshouse conditions.
As part of this glasshouse work, we quantified primary and secondary metabolites in the specific
plant parts that experienced herbivory and damage from the different kinds of herbivores in the
field. Based on our findings we conducted a controlled feeding experiment to assess the impact of
jasmonate deficiency on consumption rates by cottontail rabbits that resided in our study area.
Finally, we conducted a complementation experiment with the same rabbits to assess whether nico-
tine was the main jasmonate-dependent deterrent affecting consumption rates.

Results

Jasmonate-deficiency increases arthropod and vertebrate damage and decreases flower production

To evaluate the impact of jasmonates on herbivory-dependent plant fitness, we established three experimental *N. attenuata* populations (henceforth called 'Lytle', 'Poplar' and 'Snow') across the field station of the Lytle Ranch Preserve (St. George, UT, USA; Figure 1—figure supplement 1). In each plot, at least 12 jasmonate-deficient inverted repeat allene-oxide cyclase (irAOC) plants and empty vector controls (EV, 'wild type') were planted. The irAOC line has been characterized previously (Kallenbach et al., 2012). Its herbivory-induced jasmonate levels are reduced by more than 95% (Kallenbach et al., 2012; Machado et al., 2013) while flower production is similar to WT plants in the absence of herbivore attack (Machado et al., 2013). Five to seven weeks after the establishment of the populations, we recorded herbivore damage and counted the number of flowers on each plant as a strong predictor of Darwinian fitness (Van Dam and Baldwin, 2001; Glawe et al., 2003; Baldwin, 2003). Across all three plots, we observed four main herbivore damage types: leaf removal, stem peeling, leaf chewing and leaf sucking/piercing. The different damage types are characteristic for different herbivores including deer, rabbits, wood rats, gophers, caterpillars, ants, mirid bugs and leafhoppers (Meldau et al., 2009; Stitz et al., 2011; Kallenbach et al., 2012; Schuman et al., 2012; Đinh et al., 2013; Schäfer et al., 2013). Across all plots, jasmonate deficiency significantly increased herbivore damage and decreased flower production (Figure 1). In the Lytle plot, jasmonate-deficient plants suffered more stem peeling and leaf removal by vertebrates, but similar arthropod damage compared to wild type plants. Jasmonate-deficient plants also produced fewer flowers. Likelihood ratio tests based on Generalized Linear Models (GLMs) showed that this effect was associated with the higher occurrence of vertebrate stem peeling in irAOC plants. In the Snow plot, irAOC plants suffered more vertebrate leaf removal and damage from leaf chewing and leaf sucking/piercing insects than the wild type plants, but only leaf removal by vertebrates was associated with a reduction in flower production. In the Poplar plot, no vertebrate damage was observed, and jasmonate-deficiency increased damage by leaf sucking/piercing arthropods, but did not decrease flower production (Figure 1, Figure 1—figure supplement 2). Overall, jasmonate deficiency increased vertebrate damage more strongly than arthropod damage, and only jasmonate-dependent changes in vertebrate damage translated into a decrease in flower production.

Mimicking vertebrate damage reduces flower production more strongly than mimicking arthropod attack

To understand the impact of herbivore damage patterns and jasmonate-deficiency on plant flower
production in more detail, we mimicked the different types of damage that we observed in the field
in a controlled glasshouse experiment and quantified flower production over the entire flowering
<figure><img src="image_3.png" /><figcaption>**Figure 1.** Jasmonate-deficiency reduces flower production by increasing vertebrate damage in nature. Effects of jasmonate deficiency on vertebrate and invertebrate damage by damage type and effect of damage type on flower production across three experimental plots ('Lytle', 'Poplar' and 'Snow') are shown (n = 12–19). Solid lines indicate significant effects of jasmonate-deficiency on herbivore damage patterns and flower production. The herbivores responsible for the different damage types were identified based on field observations and characteristic feeding patterns. Jasmonate-deficiency increases damage by vertebrate and arthropod herbivores, but only vertebrate damage leads to a reduction of flower production as a strong predictor of plant reproductive potential.</figcaption></figure>

DOI: 10.7554/eLife.13720.003

Figure 1 continued on next page
Figure 1 continued

The following source data and figure supplements are available for figure 1:

**Source data 1.** Parameters used to determine the presence of different herbivores in the three experimental populations.

DOI: 10.7554/eLife.13720.004

**Source data 2.** Quantification and herbivore association of the different damage types observed in the three experimental populations.

DOI: 10.7554/eLife.13720.005

**Source data 3.** Leaf-herbivore damage screen and fitness measurements.

DOI: 10.7554/eLife.13720.006

**Figure supplement 1.** Overview of the different experimental *N. attenuata* populations used in the present study.

DOI: 10.7554/eLife.13720.007

**Figure supplement 2.** Detailed results of the analysis of the overall effect of jasmonate signaling and herbivore damage on *N. attenuata* flower production in the field.

DOI: 10.7554/eLife.13720.008

**Figure supplement 3.** Photographic evidence of herbivore damage and the associated herbivores in the Lytle plot.

DOI: 10.7554/eLife.13720.009

**Figure supplement 4.** Photographic evidence of herbivore damage and the associated herbivores in the Snow plot.

DOI: 10.7554/eLife.13720.010

period (Figure 2). Attack by chewing invertebrates was mimicked by wounding (W) the plants with a pattern wheel and treating the wounds with *Manduca sexta* oral secretions (OS) (W+OS treatments) (Machado et al., 2013). Vertebrate damage was mimicked by removing either the rosette leaves (rosette defoliation), the rosette and stem leaves (full defoliation) or the rosette and stem leaves plus the stem bark (stem peeling). The artificially peeled stems looked similar to the vertebrate-damaged stems in the field, and histological staining revealed that stem peeling resulted in the removal of the epidermis and cortex from the stems (Figure 2—figure supplement 1). In wild type plants, W+OS induction did not significantly reduce peak flower production (Figure 2a). All other treatments led to significant reductions in flower production: rosette defoliation reduced flower production by 24%, rosette and stem leaf removal reduced flower production by 53%, and stem peeling resulted in a marked 78% reduction in flower numbers. Full defoliation and stem peeling also delayed the onset of flowering by more than two weeks. Overall, the same flowering patterns were observed in EV and irAOC plants (Figure 2b). However, the total number of flowers was significantly higher in defoliated and stem-peeled irAOC plants than in similarly treated EV plants.

## Jasmonates increase secondary metabolites and decrease carbohydrates in leaves and stems

As a first step to understand the mechanism behind the significant differences in vertebrate damage between EV and irAOC plants in the field, we profiled primary and secondary metabolites in the stems and leaves of glasshouse-grown plants. Redundancy analyses (RDAs) showed that simulated herbivory led to dramatic changes in both leaf- and stem carbohydrate profiles and secondary metabolite levels in a jasmonate-dependent manner (permutation test for the effect of the genotype x treatment interaction: p<0.001 in both tissues) (Figure 3 and Figure 3—figure supplement 1). The RDAs further revealed that diterpene glycosides (DTGs), rutin, nicotine, glucose and fructose in the leaves and nicotine, glucose and fructose in the stems explained most of the metabolic differences between induced EV and irAOC plants (Figure 3). Closer inspection of these metabolites showed that under control conditions, irAOC plants have lower nicotine and DTG concentrations in the leaves and lower nicotine levels in the stems than EV plants. Upon induction, these differences became even more pronounced (Figure 3). Constitutive rutin, glucose and fructose concentrations on the other hand did not differ between genotypes, but were more strongly depleted in induced EV plants (Figure 3).

## Jasmonate-dependent nicotine determines feeding preferences of cottontail rabbits

Feeding preferences by vertebrate herbivores are influenced by plant chemistry. Toxic secondary metabolites in particular can reduce food uptake (Bryant et al., 1991; Iason, 2005; Rosenthal and Berenbaum, 2012; Camp et al., 2015). To determine whether the nicotine deficiency observed in
<figure><img src="image_4.png" /><figcaption>**Figure 2.** Simulated vertebrate attack strongly reduces plant flower production. Average (±SE) flower production of wild type (a) and jasmonate deficient irAOC plants (b) following different simulated herbivory treatments

*Figure 2 continued on next page*</figcaption></figure>
Figure 2 continued

(n = 10–12). Different letters indicate significant differences between treatments within genotypes (p<0.05). Asterisks indicate significant differences between genotypes within treatments (*p<0.05). W+OS: Wounding + application of *Manduca sexta* oral secretions.

DOI: 10.7554/eLife.13720.011

The following source data and figure supplement are available for figure 2:

**Source data 1.** Fitness costs of induction and defoliation in the glasshouse.

DOI: 10.7554/eLife.13720.012

**Figure supplement 1.** Toluene-blue staining of stem-peeled EV plants at the early bolting stage.

DOI: 10.7554/eLife.13720.013

irAOC plants could explain the higher stem and leaf-removal by vertebrates in the field, we developed a biotechnology-driven *in vitro* feeding assay for cottontail rabbits (*Sylvilagus nuttallii*) (*Figure 4—figure supplement 1*). Cottontail rabbits can feed on *N. attenuata* in nature (*Karban, 1997; Baldwin, 2003*) and were present in our field experiment (*Figure 1—figure supplement 3*). In a first step, we produced food pellets by drying *N. attenuata* shoots of wild type and irAOC plants and pressing them into pellets. To determine whether drying affects *N. attenuata* defenses, we measured secondary metabolites in plants before and after drying (*Figure 4—figure supplement 2*). The drying process did not affect the water-loss corrected abundance of DTGs. The abundances of nicotine, chlorogenic acid and rutin were reduced, but jasmonate-dependent differences were conserved. The most abundant phenol amides dicaffeoylspermidine, dicoumarylspermidine, ferolylputrescine and caffeoylputrescine were degraded during the drying process, while less abundant phenol amides increased in concentration. From these data, we concluded that the dried plant material could be used to assess the impact of jasmonate-dependent alkaloids, DTGs and phenolic compounds, but not phenol amides, on vertebrate feeding preferences. In a first set of experiments, we presented eight individual cottontail rabbits with three types of food pellets: one consisting of pure dried alfalfa leaves as a rabbit maintenance food, and two others consisting of 1) a mixture of dried alfalfa leaves and EV or 2) dried alfalfa leaves and irAOC plant material at a ratio of 5:1 to account for typical food mixing by the rabbits. We then monitored food consumption as a response variable for feeding preference over 24 hr (*Figure 4*). Comparisons of adjusted confidence intervals of asymptotes of logistic models fitted for each treatment showed that cottontails rejected EV pellets and strongly preferred to feed from both irAOC pellets and pure alfalfa (*Figure 4a, Figure 4—figure supplement 3*). Similar feeding preferences were observed in two independent experiments (*Figure 4a, Figure 4—figure supplement 4a*). From our metabolite profiling, we had identified nicotine as a prominent candidate that might drive feeding preference of vertebrates. We therefore tested the hypothesis that nicotine may be responsible for the observed feeding patterns by complementing nicotine levels of irAOC pellets to match those observed in EV pellets (*Figure 4b*, inset). Strikingly, rabbits refused both nicotine complemented-irAOC and EV pellets and preferentially fed on alfalfa pellets (*Figure 4b*). Similar results were obtained in two independent experiments (*Figure 4b, Figure 4—figure supplement 4b*).

## Discussion

Since their discovery 25 years ago (*Johnson et al., 1990*), jasmonates have emerged as key signals that govern plant responses and resistance to herbivores (*Howe and Jander, 2008*). Early experiments demonstrated that the application of jasmonates strongly reduced herbivore damage in the field and increased seed production under high herbivore pressure (*Baldwin, 1998; Agrawal, 1999*). Experiments with jasmonate biosynthesis and perception mutants subsequently revealed that a reduction in endogenous jasmonate signaling increased plant susceptibility to a large number of organisms including detritivorous crustaceans (*Farmer and Dubugnon, 2009*), caterpillars of noctuid moths (*Bodenhausen and Reymond, 2007*), spider mites (*Li et al., 2002*), beetles (*Kessler et al., 2004*) slugs (*Falk et al., 2014*) and tortoises (*Mafli et al., 2012*). However, despite the increasing number of available jasmonate mutants in different plant species, it has remained unclear whether endogenous jasmonate production actually benefits plant fitness in nature. Filling this gap of knowledge is critical to understand the evolution and natural variation of defense
<figure><img src="image_5.png" /><figcaption>**Figure 3.** Herbivory-induced jasmonates increase secondary metabolites and decrease sugars in leaves and stems. Results of a redundancy analysis of the metabolic profiles in the leaves (a) and in the stems (b) of wild type (EV) and jasmonate-deficient (irAOC) plants are shown (n = 3–5). Insets depict Figure 3 continued on next page</figcaption></figure>
Figure 3 continued

average (±SE) concentrations of metabolites explaining most of the variation between induced EV plants and the other treatment* genotype combinations (correlation coefficients > |0.8|). Letters indicate significant differences between treatments within genotypes (p<0.05). Asterisks indicate significant differences between genotypes within treatments (p<0.05). Control: intact plants; W+W: wounded and water-treated plants; W+OS: wounded and M. sexta oral secretions-treated plants; DTGs: diterpene glycosides.

DOI: 10.7554/eLife.13720.014

The following source data and figure supplement are available for figure 3:

**Source data 1.** Plant primary and secondary metabolite concentrations.

DOI: 10.7554/eLife.13720.015

**Figure supplement 1.** Secondary metabolite profiles in leaves and stems of *N. attenuata* in response to simulated herbivore attack.

DOI: 10.7554/eLife.13720.016

signaling (Machado et al., 2013; Li et al., 2015; Xu et al., 2015b; Xu et al., 2015a). Our experiments show that in the field, jasmonate-deficient plants were more strongly damaged by vertebrate and invertebrate herbivores. However, only in populations that experienced vertebrate herbivory did this effect translate into a reduction in flower production. Stem peeling in particular was associated with a jasmonate-dependent reduction in flower production. Mimicking vertebrate and arthropod herbivory in the greenhouse confirmed that the type of damage caused by vertebrates has a strong negative impact on the plant's reproductive potential.

Plant damage at later developmental stages was not assessed in the present study. Furthermore, the outcrossing success of the different genotypes was not determined due to safety precautions related to the use of transgenic plants. Earlier studies demonstrate however that leaf damage at later developmental stages has little impact on plant fitness (Zavala and Baldwin, 2006) and that flower production in the absence of future damage is a valid and robust approximation for plant fitness in the predominantly self-pollinating *N. attenuata* (Baldwin et al., 1997; Baldwin, 2003; Hettenhausen et al., 2012; Schuman et al., 2012). Our results therefore highlight the importance of jasmonates in safeguarding the plant's reproductive potential under vertebrate attack.

Vertebrate herbivores can have a profound impact on plant fitness and the distribution of plant species in nature (Hulme, 1994; 1996; Shimoda et al., 1994; Palmisano and Fox, 1997; Waller and Alverson, 1997; Gómez and Zamora, 2000; Sessions and Kelly, 2001; Warner and Cushman, 2002; Maron and Crone, 2006; Maron and Kauffman, 2006; Nisi et al., 2015). For instance, *Paliurus ramosissimus* trees that are intensively bark peeled by *sika deer* (Cervus nippon) will take nearly 30 years to recover (Ishihama et al., 2014). Similarly, strong growth delays and increased tree mortality have been observed in other species (Veblen et al., 1989; Gill, 1992; Palmisano and Fox, 1997; Vila and Guibal, 2001; Saint-Andrieux et al., 2009). Surprisingly however, vertebrate herbivores have rarely been considered in the context of plant defense signaling (Dyer, 1980; Smith et al., 1991; Liu et al., 2012). Instead, plant defense regulation has mostly been investigated in the context of arthropod herbivores and pathogens (Dixon et al., 1994; Chen et al., 1995; Yang et al., 1997; Rojo et al., 2003; Wu and Baldwin, 2009, Wu and Baldwin, 2010; Schmelz, 2015). Although arthropod damage in the field is generally low (Steppuhn et al., 2004; Rayapuram and Baldwin, 2007; Pandey and Baldwin, 2008; Rayapuram et al., 2008; Johnson et al., 2009; Meldau et al., 2009; Stitz et al., 2011; Fischer et al., 2012; Kallenbach et al., 2012; Schuman et al., 2012; Schäfer et al., 2013), local outbreaks can lead to substantial plant damage (Van Bael et al., 1999; Cease et al., 2012; DeRose and Long, 2012; Dinh et al., 2013), and several studies show that arthropods influence plant fitness (Baldwin, 1998; Maron, 1998; Kessler and Baldwin, 2004; Machado et al., 2013) and may drive the evolution of plant defenses (Huber et al., 2016; Poelman and Kessler, 2016; Züst and Agrawal, 2016) including plant defense signaling (Durrant et al., 2015; Xu et al., 2015a). Given the context-dependent potential of both arthropod and vertebrate herbivores to act as agents of natural selection, it seems crucial to include both groups when studying the ecology and evolution of plant defense regulation (Strauss, 1991; Hulme, 1994; 1996; Palmisano and Fox, 1997; Sessions and Kelly, 2001; Oduor et al., 2010).

Numerous studies have found that plant secondary metabolites can influence vertebrate foraging patterns (Freeland and Janzen, 1974; McArthur et al., 1995; Mithen et al., 1995; Dearing et al., 2005; Scogings et al., 2011; Mkhize et al., 2015). Vertebrates tend to avoid plants that are rich in
<figure><img src="image_6.png" /><figcaption>**Figure 4.** Jasmonate-dependent nicotine determines feeding preferences of a native vertebrate herbivore. Average amounts of food consumption by cottontail rabbits in a choice experiment. Choice experiments were *Figure 4 continued on next page*</figcaption></figure>
Figure 4 continued

carried out with EV, irAOC and alfalfa pellets (a) (n = 8) or with EV, irAOC+nicotine and alfalfa pellets (b) (n = 9).

Letters indicate significant differences in curve asymptotes tested based on asymptote confidence intervals of logistic models (p<0.001). Insets: nicotine contents of the different pellets. Letters indicate significant differences (p<0.05). N.D: not detected.

DOI: 10.7554/eLife.13720.017

The following source data and figure supplements are available for figure 4:

**Source data 1.** Impact of drying on secondary metabolite profiles.

DOI: 10.7554/eLife.13720.018

**Source data 2.** Nicotine concentrations in pellets.

DOI: 10.7554/eLife.13720.019

**Source data 3.** Results of choice experiment (repetition).

DOI: 10.7554/eLife.13720.020

**Source data 4.** Results of nicotine complementation experiment (repetition).

DOI: 10.7554/eLife.13720.021

**Figure supplement 1.** Experimental setup for feeding preference assays for cottontail rabbits (*Sylvilagus nuttallii*).

DOI: 10.7554/eLife.13720.022

**Figure supplement 2.** Secondary metabolite profiles in the leaves of *N. attenuata* before and after drying.

DOI: 10.7554/eLife.13720.023

**Figure supplement 3.** Asymptotes and confidence intervals of the feeding preference assays for cottontail rabbits (*Sylvilagus nuttallii*).

DOI: 10.7554/eLife.13720.024

**Figure supplement 4.** Results of the repetition of the feeding preference assays for cottontail rabbits (*Sylvilagus nuttallii*).

DOI: 10.7554/eLife.13720.025

condensed tannins, phenolics and alkaloids (Cooper and Owen-Smith, 1985; Owen Smith, 1993; Furstenburg and van Hoven, 1994; O'reilly-Wapstra et al., 2004; Jansen et al., 2007; Steppuhn et al., 2008; DeGabriel et al., 2009; Rosenthal and Berenbaum, 2012): Quinolizidine alkaloids for instance seem to protect lupin (Lupinus spp) plants from rabbits (Oryctolagus spp.) (Wink, 1985; 1988; Fattorusso and Taglialatela-Scafati, 2008). However, apart from toxic secondary metabolites, primary metabolites such as carbohydrates can also affect mammalian feeding preferences (Mayland et al., 2000). In this study, leaves and stems of jasmonate-deficient plants consistently contained lower levels of nicotine and higher levels of glucose and fructose compared to wild type plants. While both processes may have influenced vertebrate feeding preferences, we hypothesized that nicotine as a strong nervous toxin may have played a dominant role in increasing vertebrate feeding in jasmonate-deficient plants. The experiments with captive cottontails support this inference: Cottontails refused to feed on pellets containing 20% wild type *N. attenuata* shoots, but readily fed on pellets containing jasmonate-deficient, and thus nicotine-deficient, irAOC shoots. When irAOC pellets were complemented with nicotine to match the concentration in wild-type plants, rabbits avoided the two pellet types and only consumed alfalfa. This finding suggests that the amount of nicotine produced via jasmonate signaling is sufficient to dictate cottontail feeding preferences. From the plant's perspective, these results illustrate that nicotine effectively protects *N. attenuata* against sympatric vertebrates in a similar manner as it does against leaf feeding arthropods (Steppuhn and Baldwin, 2007). They further suggest that the increased damage and reduced flower production in jasmonate-deficient plants in the field can be attributed to reduced nicotine concentrations.

Nicotine is a nervous system toxin that can bind to acetylcholine receptors of synaptic cells leading to neuromuscular blockade (Albuquerque et al., 2009; Prochaska and Benowitz, 2016). Its median lethal dose (LD₅₀) varies between 3–50 mg*kg⁻¹ for small rodents (Wink, 1993; Karáčonji, 2005). As nicotine is produced at concentrations of up to 3 mg*g FW⁻¹ in *N. attenuata* shoots, the consumption of a single wild-type plant can be lethal to a small vertebrate herbivore. Interestingly, certain vertebrates such as the desert woodrat (*Neotoma lepida*) seem to excise *N. attenuata* leaves, dry them and transfer them to their nests (Baldwin, 2003). It has been proposed that this behavior may allow woodrats to use nicotine as pharmaceutical means to reduce ectoparasite loads
(Baldwin, 2003). Whether nicotine can provide benefits to vertebrates in this way remains to be
tested.

In conclusion, our study demonstrates the importance of jasmonates as regulatory signals that boost nicotine production and thereby protect the leaves and stems of *N. attenuata* from a diverse range of herbivores. By studying the interaction between plants and their natural herbivore communities, we show that jasmonates can provide fitness benefits in nature, and that these fitness benefits are derived from protecting vital plant tissues from vertebrates.

**Material and methods**

**Establishment of experimental populations**

In spring 2012, we established three experimental *N. attenuata* populations ('Lytle', 'Poplar' and 'Snow') across the field station of the Lytle Ranch Preserve (St. George, UT, USA; **Figure 1—figure supplement 1**). Each population consisted of at least twelve jasmonate-deficient inverted repeat allene-oxide cyclase plants (irAOC, line 457) and empty vector controls (EV, line A-03-9-1, 'wild type') that were planted in alternation approximately 1 m apart. The irAOC line has been characterized previously (Kallenbach et al., 2012) and was chosen for its strongly reduced jasmonate levels (Kallenbach et al., 2012; Machado et al., 2013). Plants were germinated and planted as described (Krügel et al., 2002; Schuman et al., 2012; Machado et al., 2013). Seeds of the transformed *N. attenuata* lines were imported under APHIS notification number 07-341-101n and experiments were conducted under notification number 06-242-02r. In the Snow population, half of the plants were induced by wounding the leaves with a pattern wheel and applying *Manduca sexta* oral secretions (W+OS) as described (23). As the induction treatment had no significant effect on herbivore damage or flower production, the parameter was not included in the final statistical models.

**Herbivore damage screen and fitness measurements**

Wild type (EV) and jasmonate-deficient irAOC plants from all three populations were screened for herbivore damage five weeks after transplantation. The type of damage was classified as leaf removal, stem peeling, leaf chewing and leaf sucking/piercing and was attributed to different herbivores by using knowledge from previous studies (Baldwin, 2003; Kessler et al., 2004; Steppuhn et al., 2004; Meldau et al., 2009; Stitz et al., 2011; Kallenbach et al., 2012; Schuman et al., 2012; Đinh et al., 2013; Schäfer et al., 2013). The parameters used to determine the presence of herbivores and details about the quantification of the different damage types are summarized in Figure 1—source datas 1, 2. Photographic evidence of some herbivores is provided in Figure 1—figure supplements 3, 4. As a measure of reproductive potential, the numbers of flowers were counted for each plant in the different populations immediately following the herbivore impact assessments (Van Dam and Baldwin, 2001; Glawe et al., 2003; Sime and Baldwin, 2003). Flowers were removed after counting according to APHIS regulations for the use of transgenic plants in the field. In the plot surrounding the Snow population, we observed a spike in gopher activity in the week after the flower count. To be able to include this attack into our statistical model, we reassessed herbivore damage and counted the number of regrowing flowers two weeks later. This second flower count was then used for statistical analysis.

**Fitness costs of induction and defoliation in the glasshouse**

To determine the relative value of leaves and stem bark for plant reproduction, we evaluated reproductive output in plants that were subjected to four different regimes of simulated herbivory: simulated insect attack, stem peeling, and full and partial defoliation (n =10-12). Wild type (EV) and irAOC plants were grown as described elsewhere (Krügel et al., 2002). When plants reached the bolting stage, insect attack by chewing herbivores was mimicked by wounding 3 leaves and applying *M. sexta* oral secretions as described (W+OS) (Machado et al., 2015). In another batch of plants, vertebrate herbivory was simulated by removing either the rosette leaves, or by removing both rosette and stem-leaves (i.e., full defoliation). Finally, we mimicked the vertebrate damage patterns observed in the Lytle plot population by removing all leaves and peeling off the bark. To remove the bark, a small incision was made into the epidermis at the stem base. The epidermis and bark were then pulled off with forceps from the bottom to the top, leaving only xylem and pith of the stems
(See below 'Histological staining')*(Figure 2—figure supplement 1). Intact plants were used as controls. Following the different treatments, all plants were left to grow in a fully randomized setup, and the number of flowers was counted every 3-4 days until the end of the flowering period.

## Plant primary and secondary metabolite profiling

As a first step to understand the mechanism behind the significant differences in vertebrate grazing patterns between wild type (EV) and jasmonate-deficient irAOC plants in the field, we profiled primary and secondary metabolites in the stems and leaves for several plant treatments. Prolonged leaf-wounding and simulated insect attack treatments were carried out as described (Machado et al., 2013; 2015) (n = 3–5). Nicotine, chlorogenic acid, rutin and total 17-hydroxygeranyllinalool diterpene glycosides were determined by HPLC-DAD as described (Keinänen et al., 2001). Sugars and starch were measured as described (Velterop and Vos, 2001; Smith and Zee-man, 2006; Machado et al., 2013; 2015).

## Cottontail rabbit feeding preference

To determine which of the jasmonate-dependent changes in host plant chemistry may be responsible for the preference for irAOC plants observed in the field, we developed an in vitro feeding assay for cottontail rabbits (Sylvilagus nuttallii) as follows.

### Effect of drying process on N. attenuata secondary metabolite profiles

To test the effect of drying on N. attenuata secondary metabolite profiles, we measured secondary metabolite levels in rosette stage plants that were both harvested and immediately frozen in liquid nitrogen or were harvested and dried (50°C, 24 hr) prior to secondary metabolite measurements (n = 5). Secondary metabolites analyses were carried out as described (Ferrieri et al., 2015). To compare the concentrations between fresh and dried plant materials, we measured water loss during drying and corrected the measured concentrations of the dried materials accordingly.

### Food pellet preparation

Stems and leaves of early elongated EV and irAOC N. attenuata plants were harvested and dried (50°C, 24 hr). Dried plant material was then ground using a Thomas Scientific Wiley Mill equipped with a 2-mm sieve (Thomas Scientific, Swedesboro, NJ, USA). The ground material was then mixed (1:5) with Alfalfa Purina rabbit chow (Purina Mills, St. Louis, Missouri, USA) and pelletized by a Buskirk PM605 Lab Mill customized with 3-mm wide holes (Buskirk engineering, Ossian, IN, USA). Alfalfa Purina rabbit chow ('Alfalfa') was used as the rabbit's maintenance food. All three pellet types were divided into two batches. The first batch was left intact and used for the first set of experiments and the second one was used for the nicotine complementation experiments. To complement nicotine levels in irAOC pellets, we first measured nicotine content in the pellets (Keinänen et al., 2001) and then, by evenly spraying a solution of nicotine on irAOC pellets, we equalized the nicotine levels of irAOC to those observed in EV pellets (Figure 4a,b: insets). EV and alfalfa pellets were sprayed with pure water. Pellets were allowed to dry at room temperature and nicotine levels were measured again (n = 5).

### Feeding preference assay

In a first set of experiments, we tested whether our field observations could be reproduced by our in vitro feeding assay with cottontail rabbits. For this, we presented rabbits with 20 g (±0.1 g) of each pellet type (EV, irAOC and alfalfa) and allowed them to choose and feed freely for a 24 hr period (Figure 4—figure supplement 1). Eight rabbits were used in total. One animal was tested at a time in 24 hr-intervals and the full experiment was carried out twice for each rabbit. Food consumption was set as response variable for feeding preference. To quantify food consumption, pellets were placed separately in 9-cm wide and 6-cm deep circular bowls that sat on the top of a balance (A&D EJ3000) (Figure 4—figure supplement 1). The balances were adjacent to each other and enclosed in a plywood box (66 cm wide, 38 cm long, and 16 cm high). Automatic recordings of pellet weight were logged on a laptop computer and plotted against time. In a second set of experiments, we tested whether differences in nicotine content between EV and irAOC plants might be responsible for the feeding preferences observed. For this, we carried out a procedure as described above, but
we contrasted feeding preference of cottontails on irAOC+nicotine and EV pellets. The rabbit's maintenance food ('Alfalfa') was used as a control in both sets of experiments. All the choice experiments were carried out using the same rabbits. Trials were approved by the International Care and Use Committees (IACUC) via permit 04495-002. The feeding trials were carried out in the Washington State University's Small Mammal Research Facility, located 4 km east of Pullman, Washington.

## Histological staining

To investigate which tissues were removed from the stems by vertebrates, we simulated this type of herbivory by experimentally removing the epidermal tissue from bolting wild type plants as described above. Peeled and non-peeled stems where then stained with toluene blue for 15 s and photographed using a stereomicroscope equipped with a digital CCD camera (SteREO Discovery. V8, 14 Carl Zeiss Microimaging) and processed with AxioVision LE software (Carl Zeiss 15 Microimaging).

## Statistical analyses

All statistical analyses were performed using the R software (R Development Core Team, 2015) unless otherwise stated:

### Leaf-herbivore damage screen and fitness measurements

Since different herbivore communities were observed in the different plots, the effect of damage intensity on flower production was assessed for the three experimental populations separately. ANOVAs were used to compare damage intensity between plant genotypes, for each damage type. The influence of the different damage types on plant fitness (i.e. flower production) was assessed using likelihood ratio tests applied on Generalized Linear Models (GLM; family: negative binomial, link: log) (package 'MASS', function 'glm.nb' (Venables and Ripley, 2002)). Flower production was compared between plant genotypes using ANOVAs by considering the natural logarithm of the number of flowers as response variable.

### Fitness costs of induction and defoliation in the glasshouse

The effect of the different simulated herbivore treatments on flower production was analyzed separately for each plant genotype using a likelihood ratio test applied on a GLM (family: negative binomial, link: log) considering the number of flowers produced the day of flowering peak as response variable (EV: at day 22 for control, W+OS and rosette defoliation, and at day 26 for full defoliation and stem peeling. irAOC: at day 17 for control, W+OS and rosette defoliation, and at day 26 for full defoliation and stem peeling). Pairwise comparisons of Least Squares Means (LSMeans; package 'lsmeans', function 'lsmeans') (Lenth, 2015) were computed using the False Discovery Rate (FDR) (Benjamini and Hochberg, 1995) correction for *P*-values. Plant genotypes were compared separately for each treatment using the same procedure, and *P*-values of all tests were FDR-adjusted.

### Plant primary and secondary metabolite profiling

To assess jasmonate-dependent differences in stem and leaf metabolic profiles, a redundancy analysis (RDA), considering as constraints the plant genotype, the treatment and the interaction between these two factors (package 'vegan', function 'rda') was carried out for leaves and stems separately (Oksanen et al., 2015). Constraints explained 70.2% and 61.0% of total variance in leaves and stems, respectively. Significance of these three terms was assessed using a permutation test (with 9999 permutations) as described (Legendre and Legendre, 2012). The interaction term had a significant effect in both RDAs (p<0.001 in both cases). Metabolites explaining most of the variation between induced EV plants (EV:W+W and EV:W+OS) and the other treatment genotype combinations (EV:Control, irAOC:Control, irAOC:W+W and irAOC:W+OS) were identified by computing the Pearson correlation coefficient between chemical data and sample scores by using the directions in the two dimensional vector spaces which discriminate these treatments (axis 1 in the leaves and to the bottom right-to-top left diagonal in the stems). Metabolites with correlation coefficient values > |0.8| were then included as bar graphs in Figure 3. To test the effect of simulated herbivory on N. attenuata secondary metabolite profiles, two-way ANOVAs with treatment and genotype as factors were carried out with Sigma Plot 12.0 (Systat Software Inc., San Jose,
CA, USA). Levene's and Shapiro-Wilk tests were applied to determine error variance and normality. Holm-Sidak post hoc tests were used for multiple comparisons.

## Cottontail rabbit feeding preference

To test rabbit feeding preference, a logistic model was built for each pellet type separately, using the average amount of pellet consumed by 8 rabbits at each time point as response variable. Asymptote values of these models were considered as treatment effects. These asymptotes were compared using their Bonferroni-adjusted 95% confidence intervals. Non-overlapping confidence intervals show significantly different asymptotes. Concentration of nicotine in pellets was compared between treatments using pairwise t-tests, using the FDR correction for P-values.

## Effect of the drying process on N. attenuata secondary metabolite profiles

To test the effect of drying on N. attenuata secondary metabolite profiles, two-way ANOVAs with water status and genotype as factors were carried out with Sigma Plot 12.0 (Systat Software Inc., San Jose, CA, USA). Water-loss corrected values were used for all analyses. Levene's and Shapiro-Wilk tests were applied to determine error variance and normality. Holm-Sidak post hoc tests were used for multiple comparisons.

## Acknowledgements

Celia Diezel, Youngjoo Oh and Danny Kessler helped with the field experiments. Abigail Ferrieri and Rajarajan Ramakrishnan supported the metabolite analyses. Arne Weinhold provided photographic evidence for gopher grazing. Pavan Kumar provided photographic evidence for the presence of cot-tontail rabbits. Lucas Cortés Llorca assisted with the histological staining. The work of ME was supported by a Marie Curie Intra European Fellowship (grant no. 273107). This work was supported by the Max Planck Society and the Collaborative Research Centre ChemBioSys (CRC 1127 ChemBioSys) of the Deutsche Forschungsgemeinschaft (DFG).

## Additional information

### Competing interests

ITB: Senior Editor, eLife. The other authors declare that no competing interests exist.

### Funding

<table><thead><tr><td>Funder</td><td>Grant reference number</td><td>Author</td></tr></thead><tbody><tr><td>Marie Curie Intra European Scholarship</td><td>273107</td><td>Matthias Erb</td></tr><tr><td>Deutsche Forschungsgemeinschaft</td><td>CRC 1127 ChemBioSys</td><td>Ian T Baldwin</td></tr><tr><td>Institut National de la Recherche Agronomique</td><td></td><td>Maxime R Hervé</td></tr></tbody></table>

The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.

### Author contributions

RARM, MM, ME, Conception and design, Acquisition of data, Analysis and interpretation of data, Drafting or revising the article; MRH, Analysis and interpretation of data, Drafting or revising the article; ITB, Conception and design, Analysis and interpretation of data, Drafting or revising the article

### Author ORCIDs

Ian T Baldwin, http://orcid.org/0000-0001-5371-2974

Matthias Erb, http://orcid.org/0000-0002-4446-9834
## Ethics

Animal experimentation: Trials were approved by the International Care and Use Committees (IACUC) via permit 04495-002. The feeding trials were carried out in the Washington State University's Small Mammal Research Facility, located 4 km east of Pullman, Washington. Field experiments: Seeds of the transformed *N. attenuata* lines were imported under APHIS notification number 07-341-101n and experiments were conducted under notification number 06-242-02r.

## References

* Abrahamson WG. 1975. Reproduction of rubus hispidus L. in different habitats. *American Midland Naturalist* **93**: 471-478. doi: 10.2307/2424184

* Agrawal AA, Hastings AP, Johnson MTJ, Maron JL, Salminen J-P. 2012. Insect herbivores drive real-time ecological and evolutionary change in plant populations. *Science* **338**:113–116. doi: 10.1126/science.1225977

* Agrawal AA. 1998. Induced responses to herbivory and increased plant performance. *Science* **279**:1201–1202. doi: 10.1126/science.279.5354.1201

* Agrawal AA. 1999. Induced responses to herbivory in wild radish: Effects on several herbivores and plant fitness. *Ecology* **80**:1713–1723. doi: 10.1890/0012-9658(1999)080[1713:IRTHIW]2.0.CO;2

* Albuquerque EX, Pereira EF, Alkondon M, Rogers SW. 2009. Mammalian nicotinic acetylcholine receptors: from structure to function. *Physiological Reviews* **89**:73–120. doi: 10.1152/physrev.00015.2008

* Baldwin IT, Preston C, Euler M, Gorham D. 1997. Patterns and consequences of benzyl acetone floral emissions from nicotiana attenuata plants. *Journal of Chemical Ecology* **23**:2327–2343. doi: 10.1023/B:JOEC.0000006677.56380.cd

* Baldwin IT. 1998. Jasmonate-induced responses are costly but benefit plants under attack in native populations. *Proceedings of the National Academy of Sciences of the United States of America* **95**:8113–8118. doi: 10.1073/pnas.95.14.8113

* Baldwin IT. 2003. Curing of leaves by small mammals does not decrease nicotine contents. *Western North American Naturalist* **63**:114–117.

* Becerra JX. 2015. On the factors that promote the diversity of herbivorous insects and plants in tropical forests. *Proceedings of the National Academy of Sciences of the United States of America* **112**:6098–6103. doi: 10.1073/pnas.1418643112

* Benjamini Y, Hochberg Y. 1995. Controlling the false discovery rate: A practical and powerful approach to multiple testing. *Journal of the Royal Statistical Society. Series B*:289-300.

* Bergman M. 2002. Can saliva from moose, alces alces, affect growth responses in the sallow, *Oikos* **96**:164–168. doi: 10.1034/j.1600-0706.2002.960118.x

* Bodenhausen N, Reymond P. 2007. Signaling pathways controlling induced resistance to insect herbivores in arabidopsis. *Molecular Plant-Microbe Interactions* **20**:1406-1420. doi: 10.1094/MPMI-20-11-1406

* Bruelheide H, Scheidel U. 1999. Slug herbivory as a limiting factor for the geographical range of *arnica montana*. *Journal of Ecology* **87**:839–848. doi: 10.1046/j.1365-2745.1999.00403.x

* Bryant JP, Provenza FD, Pastor J, Reichardt PB, Clausen TP, du Toit JT. 1991. Interactions between woody plants and browsing mammals mediated by secondary metabolites. *Annual Review of Ecology and Systematics* **22**:431–446. doi: 10.1146/annurev.es.22.110191.002243

* Camp MJ, Shipley LA, Johnson TR, Forbey JS, Rachlow JL, Crowell MM. 2015. Modeling trade-offs between plant fiber and toxins: A framework for quantifying risks perceived by foraging herbivores. *Ecology* **96**:3292–3302. doi: 10.1890/14-2412.1

* Cease AJ, Elser JJ, Ford CF, Hao S, Kang L, Harrison JF. 2012. Heavy livestock grazing promotes locust outbreaks by lowering plant nitrogen content. *Science* **335**:467–469. doi: 10.1126/science.1214433

* Chen Z, Malamy J, Henning J, Conrath U, Sanchez-Casas P, Silva H, Ricigliano J, Klessig DK. 1995. Induction, modification, and transduction of the salicylic acid signal in plant defense responses. *Proceedings of the National Academy of Sciences of the United States of America* **92**:4134–4137. doi: 10.1073/pnas.92.10.4134

* Collins SL, Knapp AK, Briggs JM, Blair JM, Steinauer EM. 1998. Modulation of diversity by grazing and mowing in native tallgrass prairie. *Science* **280**:745–747. doi: 10.1126/science.280.5364.745

* Cooper SM, Owen-Smith N. 1985. Condensed tannins deter feeding by browsing ruminants in a south african savanna. *Oecologia* **67**:142–146. doi: 10.1007/BF00378466

* Corbett BP, Jia L, Sayler RJ, Arevalo-Soliz LM, Goggin F. 2011. The effects of root-knot nematode infection and mi-mediated nematode resistance in tomato on plant fitness. *Journal of Nematology* **43**:82–89.

* Dearing MD, Foley WJ, McLean S. 2005. The influence of plant secondary metabolites on the nutritional ecology of herbivorous terrestrial vertebrates. *Annual Review of Ecology, Evolution, and Systematics* **36**:169–189. doi: 10.1146/annurev.ecolsys.36.102003.152617

* DeGabriel JL, Moore BD, Foley WJ, Johnson CN. 2009. The effects of plant defensive chemistry on nutrient availability predict reproductive success in a mammal. *Ecology* **90**:711–719. doi: 10.1890/08-0940.1

* DeRose RJ, Long JN. 2012. Factors influencing the spatial and temporal dynamics of engelmann spruce mortality during a spruce beetle outbreak on the markagunt plateau, utah. *Forest Science* **58**:1–14. doi: 10.5849/forsci.10-079

* Dixon RA, Harrison MJ, Lamb CJ. 1994. Early events in the activation of plant defense responses. *Annual Review of Phytopathology* **32**:479–501. doi: 10.1146/annurev.py.32.090194.002403
Dinh ST, Gális I, Baldwin IT, Dinh S. 2013. UVB radiation and 17-hydroxygeranyllinalool diterpene glycosides provide durable resistance against mirid (tupiocoris notatus) attack in field-grown. *Plant, Cell & Environment* **36**:590–606. doi: 10.1111/j.1365-3040.2012.02598.x

Durrant M, Boyer J, Baldwin IT, Xu S. 2015. An evolutionary hourglass of herbivore-induced transcriptomic responses in Nicotiana attenuata. *bioRxiv*. doi: 10.1101/034603

Dyer MI. 1980. Mammalian epidermal growth factor promotes plant growth. *Proceedings of the National Academy of Sciences of the United States of America* **77**:4836–4837. doi: 10.1073/pnas.77.8.4836

Falk KL, Kästner J, Bodenhausen N, Schramm K, Paetz C, Vassão DG, Reichelt M, von Knorre D, Bergelson J, Erb M, Gershenzon J, Meldau S. 2014. The role of glucosinolates and the jasmonic acid pathway in resistance of arabidopsis thaliana against molluscan herbivores. *Molecular Ecology* **23**:1188–1203. doi: 10.1111/mec.12610

Farmer EE, Dubugnon L. 2009. Detritivorous crustaceans become herbivores on jasmonate-deficient plants. *Proceedings of the National Academy of Sciences of the United States of America* **106**:935–940. doi: 10.1073/pnas.0812182106

Fattorusso E, Taglialatela-Scafati O. 2008. *Modern Alkaloids: Structure, Isolation, Synthesis, and Biology*. John Wiley & Sons.

Ferrieri AP, Arce CC, Machado RA, Meza-Canales ID, Lima E, Baldwin IT, Erb M. 2015. A nicotiana attenuata cell wall invertase inhibitor (nacwii) reduces growth and increases secondary metabolite biosynthesis in herbivore-attacked plants. *The New Phytologist* **208**:519–530. doi: 10.1111/nph.13475

Fischer M, Weyand A, Rudmann-Maurer K, Stöcklin J. 2012. Omnipresence of leaf herbivory by invertebrates and leaf infection by fungal pathogens in agriculturally used grasslands of the swiss alps, but low plant damage. *Alpine Botany* **122**:95–107. doi: 10.1007/s00035-012-0108-x

Freeland WJ, Janzen DH. 1974. Strategies in herbivory by mammals: The role of plant secondary compounds. *The American Naturalist* **108**:269–289. doi: 10.1086/282907

Furstenburg D, van Hoven W. 1994. Condensed tannin as anti-defoliate agent against browsing by giraffe (giraffe camelopardalis) in the kruger national park. *Comparative Biochemistry and Physiology Part A: Physiology* **107**:425–431. doi: 10.1016/0300-9629(94)90402-2

Gill RM. 1992. A review of damage by mammals in north temperate forests: 3. impact on trees and forests. *Forestry* **65**:363–388. doi: 10.1093/forestry/65.4.363-a

Glawe GA, Zavala JA, Kessler A, Van Dam NM, Baldwin IT. 2003. ECOLOGICAL costs and benefits correlated with trypsin protease inhibitor production in nicotiana attenuata. *Ecology* **84**:79–90. doi: 10.1890/0012-9658 (2003)084[0079:ECABCW]2.0.CO;2

Gómez JM, Zamora R. 2000. Differential impact of vertebrate and invertebrate herbivores on the reproductive output of hormathophylla spinosa. *EcoScience* **7**:299–306.

Hanley M. 1998. Seedling herbivory, community composition and plant life history traits. *Perspectives in Plant Ecology, Evolution and Systematics* **1**:191–205. doi: 10.1078/1433-8319-00058

Heil M, Koch T, Hilpert A, Fiala B, Boland W, Linsenmair KE. 2001. Extrafloral nectar production of the ant-associated plant, macaranga tanarius, is an induced, indirect, defensive response elicited by jasmonic acid. *Proceedings of the National Academy of Sciences of the United States of America* **98**:1083–1088. doi: 10.1073/pnas.98.3.1083

Hettenhausen C, Baldwin IT, Wu J. 2012. Silencing MPK4 in nicotiana attenuata enhances photosynthesis and seed production but compromises abscisic acid-induced stomatal closure and guard cell-mediated resistance to pseudomonas syringae pv tomato DC3000. *Plant Physiology* **158**:759–776. doi: 10.1104/pp.111.190074

Hodgson J, Illius AW. 1996. *The Ecology and Management of Grazing Systems*. Oxford: CAB international Press.

Howe GA, Jander G. 2008. Plant immunity to insect herbivores. *Annual Review of Plant Biology* **59**:41–66. doi: 10.1146/annurev.arplant.59.032607.092825

Huber M, Bont Z, Fricke J, Brillatz T, Aziz Z, Gershenzon J, Erb M. 2016b. A below-ground herbivore shapes root defensive chemistry in natural plant populations. *Proceedings of the Royal Society B: Biological Sciences* **283**: 20160285. doi: 10.1098/rspb.2016.0285

Huber M, Epping J, Schulze Gronover C, Fricke J, Aziz Z, Brillatz T, Swyers M, Köllner TG, Vogel H, Hammerbacher A, Triebwasser-Freese D, Robert CA, Verhoeven K, Preite V, Gershenzon J, Erb M, Gronover CS. 2016a. A latex metabolite benefits plant fitness under root herbivore attack. *PLoS Biology* **14**:e1002332. doi: 10.1371/journal.pbio.1002332

Hulme PE. 1994. Seedling herbivory in grassland: Relative impact of vertebrate and invertebrate herbivores. *The Journal of Ecology* **82**:873–880. doi: 10.2307/2261451

Hulme PE. 1996. Herbivores and the performance of grassland plants: A comparison of arthropod, mollusc and rodent herbivory. *The Journal of Ecology* **84**:43–51. doi: 10.2307/2261698

Iason G. 2005. The role of plant secondary metabolites in mammalian herbivory: Ecological perspectives. *Proceedings of the Nutrition Society* **64**:123–131. doi: 10.1079/PNS2004415

Ishihama F, Fujii S, Yamamoto K, Takada T. 2014. Estimation of dieback process caused by herbivory in an endangered root-sprouting shrub species, paliurus ramosissimus (lour.) poir., using a shoot-dynamics matrix model. *Population Ecology* **56**:275–288. doi: 10.1007/s10144-013-0414-1

Jansen DAWAM, van Langevelde F, de Boer WF, Kirkman KP. 2007. Optimisation or satiation, testing diet selection rules in goats. *Small Ruminant Research* **73**:160–168. doi: 10.1016/j.smallrumres.2007.01.012

Johnson MTJ, Smith SD, Rausher MD. 2009. Plant sex and the evolution of plant defenses against herbivores. *Proceedings of the National Academy of Sciences of the United States of America* **106**:18079–18084. doi: 10.1073/pnas.0904695106
Johnson R, Lee JS, Ryan CA. 1990. Regulation of expression of a wound-inducible tomato inhibitor I gene in transgenic nightshade plants. *Plant Molecular Biology* **14**:349–356. doi: 10.1007/BF00028771

Kallenbach M, Bonaventure G, Gilardoni PA, Wissgott A, Baldwin IT. 2012. Empoasca leafhoppers attack wild tobacco plants in a jasmonate-dependent manner and identify jasmonate mutants in natural populations. *Proceedings of the National Academy of Sciences of the United States of America* **109**:E1548–E1557. doi: 10.1073/pnas.1200363109

Karačonjić IB. 2005. Facts about nicotine toxicity. *Arhiv Za Higijenu Rada I Toksikologiju* **56**:363–371.

Karban R. 1997. Neighbourhood affects a plant's risk of herbivory and subsequent success. *Ecological Entomology* **22**:433–439. doi: 10.1046/j.1365-2311.1997.00081.x

Kato T, Ishida K, Sato H. 2008. The evolution of nettle resistance to heavy deer browsing. *Ecological Research* **23**:339–345. doi: 10.1007/s11284-007-0387-7

Keinänen M, Oldham NJ, Baldwin IT. 2001. Rapid HPLC screening of jasmonate-induced increases in tobacco alkaloids, phenolics, and diterpene glycosides in nicotiana attenuata. *Journal of Agricultural and Food Chemistry* **49**:3553–3558. doi: 10.1021/jf010200+

Kessler A, Baldwin IT. 2004. Herbivore-induced plant vaccination. part I. the orchestration of plant defenses in nature and their fitness consequences in the wild tobacco nicotiana attenuata. *The Plant Journal* **38**:639–649. doi: 10.1111/j.1365-313X.2004.02076.x

Kessler A, Halitschke R, Baldwin IT. 2004. Silencing the jasmonate cascade: Induced plant defenses and insect populations. *Science* **305**:665–668. doi: 10.1126/science.1096931

Krügel T, Lim M, Gase K, Halitschke R, Baldwin IT. 2002. Agrobacterium-mediated transformation of nicotiana attenuata, a model ecological expression system. *Chemoecology* **12**:177–183. doi: 10.1007/PL00012666

Legendre P, Legendre L. 2012. *Numerical Ecology*. Amsterdam: Elsevier.

Lenth R. 2015. *Ismeans: Least-squares means*. R Package Version 2.20-23.

Li C, Williams MM, Loh YT, Lee GI, Howe GA. 2002. Resistance of cultivated tomato to cell content-feeding herbivores is regulated by the octadecanoid-signaling pathway. *Plant Physiology* **130**:494–503. doi: 10.1104/pp.005314

Li D, Baldwin IT, Gaquerel E. 2015. Navigating natural variation in herbivory-induced secondary metabolism in coyote tobacco populations using MS/MS structural analysis. *Proceedings of the National Academy of Sciences of the United States of America* **112**:E4147–E4155. doi: 10.1073/pnas.1503106112

Li L, Zhao Y, McCaig BC, Wingerd BA, Wang J, Whalon ME, Pichersky E, Howe GA. 2004. The tomato homolog of CORONATINE-INSENSITIVE1 is required for the maternal control of seed maturation, jasmonate-signaled defense responses, and glandular trichome development. *The Plant Cell* **16**:126–143. doi: 10.1105/tpc.017954

Liu J, Wang L, Wang D, Bonser SP, Sun F, Zhou Y, Gao Y, Teng X. 2012. Plants can benefit from herbivory: Stimulatory effects of sheep saliva on growth of leymus chinensis. *PloS One* **7**:e29259. doi: 10.1371/journal.pone.0029259

Machado RA, Arce CC, Ferrieri AP, Baldwin IT, Erb M. 2015. Jasmonate-dependent depletion of soluble sugars compromises plant resistance to manduca sexta. *The New Phytologist* **207**:91–105. doi: 10.1111/nph.13337

Machado RA, Ferrieri AP, Robert CA, Glauser G, Kallenbach M, Baldwin IT, Erb M, Am Robert C. 2013. Leaf-herbivore attack reduces carbon reserves and regrowth from the roots via jasmonate and auxin signaling. *The New Phytologist* **200**:1234–1246. doi: 10.1111/nph.12438

Mafli A, Goudet J, Farmer EE. 2012. Plants and tortoises: Mutations in the arabidopsis jasmonate pathway increase feeding in a vertebrate herbivore. *Molecular Ecology* **21**:2534–2541. doi: 10.1111/j.1365-294X.2012.05533.x

Maron JL, Crone E. 2006. Herbivory: Effects on plant abundance, distribution and population growth. *Proceedings. Biological Sciences / the Royal Society* **273**:2575–2584. doi: 10.1098/rspb.2006.3587

Maron JL, Kauffman MJ. 2006. Habitat-specific impacts of multiple consumers on plant population dynamics. *Ecology* **87**:113–124. doi: 10.1890/05-0434

Maron JL. 1998. INSECT herbivory above- and belowground: Individual and joint effects on plant fitness. *Ecology* **79**:1281–1293. doi: 10.1890/0012-9658(1998)079[1281:IIHAABI]2.0.CO;2

Mayland HF, Shewmaker GE, Harrison PA, Chatterton NJ. 2000. Nonstructural Carbohydrates in Tall Fescue Cultivars. *Agronomy Journal* **92**:1203–1206. doi: 10.2134/agronj2000.9261203x

McArthur C, Sanson GD, Beal AM. 1995. Salivary proline-rich proteins in mammals: Roles in oral homeostasis and counteracting dietary tannin. *Journal of Chemical Ecology* **21**:663–691. doi: 10.1007/BF02033455

McConn M, Creelman RA, Bell E, Mullet JE, Browse J. 1997. Jasmonate is essential for insect defense in arabidopsis. *Proceedings of the National Academy of Sciences of the United States of America* **94**:5473–5477. doi: 10.1073/pnas.94.10.5473

Meldau S, Wu J, Baldwin IT. 2009. Silencing two herbivory-activated MAP kinases, SIPK and WIPK, does not increase nicotiana attenuata's susceptibility to herbivores in the glasshouse and in nature. *The New Phytologist* **181**:161–173. doi: 10.1111/j.1469-8137.2008.02645.x

Mithen R, Raybould AF, Giamoustaris A. 1995. Divergent selection for secondary metabolites between wild populations of brassica oleracea and its implications for plant-herbivore interactions. *Heredity* **75**:472–484. doi: 10.1038/hdy.1995.164

Mkhize NR, Heitkönig IMA, Scogings PF, Dziba LE, Prins HHT, de Boer WF. 2015. Condensed tannins reduce browsing and increase grazing time of free-ranging goats in semi-arid savannas. *Applied Animal Behaviour Science* **169**:33–37. doi: 10.1016/j.applanim.2015.04.012
Morris WF, Hufbauer RA, Agrawal AA, Bever JD, Borowicz VA, Gilbert GS, Maron JL, Mitchell CE, Parker IM, Power AG, Torchin ME, Vázquez DP. 2007. Direct and interactive effects of enemies and mutualists on plant performance: A meta-analysis. *Ecology* **88**:1021–1029. doi: 10.1890/06-0442

Nisi AC, Hernández DL, English LP, Rogers ES. 2015. Patterns of selective herbivory on five prairie legume species. *The American Midland Naturalist* **173**:110–121. doi: 10.1674/0003-0031-173.1.110

O'Reilly-Wapstra JM, Mcarthur C, Potts BM. 2004. Linking plant genotype, plant defensive chemistry and mammal browsing in a eucalyptus species. *Functional Ecology* **18**:677–684. doi: 10.1111/j.0269-8463.2004.00887.x

Oduor AMO, Gómez JM, Strauss SY. 2010. Exotic vertebrate and invertebrate herbivores differ in their impacts on native and exotic plants: A meta-analysis. *Biological Invasions* **12**:407–419. doi: 10.1007/s10530-009-9622-1

Oksanen J, Blanchet G, Kindt R, Legendre P. 2015. *Vegan: Community ecology package*. R package version 2.3-1.

Onkokesung N, Gális I, von Dahl CC, Matsuoka K, Saluz HP, Baldwin IT. 2010. Jasmonic acid and ethylene modulate local responses to wounding and simulated herbivory in nicotiana attenuata leaves. *Plant Physiology* **153**:785–798. doi: 10.1104/pp.110.156232

Owen Smith N. 1993. Woody plants, browsers and tannins in southern african savannas. *South African Journal of Science* **89**:505–510.

Paige KN, Whitham TG. 1987. Overcompensation in response to mammalian herbivory: The advantage of being eaten. *The American Naturalist* **129**:407–416. doi: 10.1086/284645

Palmisano S, Fox LR. 1997. Effects of mammal and insect herbivory on population dynamics of a native californian thistle, cirsium occidentale. *Oecologia* **111**:413–421. doi: 10.1007/s004420050253

Pandey SP, Baldwin IT. 2008. Silencing rna-directed RNA polymerase 2 increases the susceptibility of nicotiana attenuata to UV in the field and in the glasshouse. *The Plant Journal* **54**:845–862. doi: 10.1111/j.1365-313X.2008.03450.x

Paschold A, Halitschke R, Baldwin IT. 2007. Co(i)-ordinating defenses: Nacoi1 mediates herbivore- induced resistance in nicotiana attenuata and reveals the role of herbivore movement in avoiding defenses. *The Plant Journal* **51**:79–91. doi: 10.1111/j.1365-313X.2007.03119.x

Poelman EH, Broekgaarden C, Van Loon JJ, Dicke M. 2008. Early season herbivore differentially affects plant defence responses to subsequently colonizing herbivores and their abundance in the field. *Molecular Ecology* **17**:3352–3365. doi: 10.1111/j.1365-294X.2008.03838.x

Poelman EH, Kessler A. 2016. Keystone Herbivores and the Evolution of Plant Defenses. *Trends in Plant Science* **21**:477–485. doi: 10.1016/j.tplants.2016.01.007

Prochaska JJ, Benowitz NL. 2016. The Past, Present, and Future of Nicotine Addiction Therapy. *Annual Review of Medicine* **67**:467–486. doi: 10.1146/annurev-med-111314-033712

Rayapuram C, Baldwin IT. 2007. Increased SA in npr1-silenced plants antagonizes JA and ja-dependent direct and indirect defenses in herbivore-attacked nicotiana attenuata in nature. *The Plant Journal* **52**:700–715. doi: 10.1111/j.1365-313X.2007.03267.x

Rayapuram C, Wu J, Haas C, Baldwin IT. 2008. PR-13/thionin but not PR-1 mediates bacterial resistance in nicotiana attenuata in nature, and neither influences herbivore resistance. *Molecular Plant-Microbe Interactions* **21**:988–1000. doi: 10.1094/MPMI-21-7-0988

R Development Core Team. 2015. *R: A language and environment for statistical computing*. Vienna, Austria: R foundation for statistical computing. http://www.r-project.org.

Redman AM, Jr. DFC, Schultz JC, Cipollini Jr DF. 2001. Fitness costs of jasmonic acid-induced defense in tomato, lycopersicon esculentum. *Oecologia* **126**:380–385. doi: 10.1007/s004420000522

Rojo E, Solano R, Sánchez-Serrano JJ. 2003. Interactions between signaling compounds involved in plant defense. *Journal of Plant Growth Regulation* **22**:82–98. doi: 10.1007/s00344-003-0027-6

Rosenthal GA, Berenbaum MR. 2012. *Herbivores: Their Interactions with Secondary Plant Metabolites: Ecological and Evolutionary Processes*. Academic Press.

Saint-Andrieux C, Bonenfant C, Toïgo C, Basille M, Klein F. 2009. Factors affecting beech tica bark stripping by red deer us elaphus in a mixed forest. *Wildlife Biology* **15**:187–196. doi: 10.2981/07-100

Schmelz EA. 2015. Impacts of insect oral secretions on defoliation-induced plant defense. *Current Opinion in Insect Science* **9**:7–15. doi: 10.1016/j.cois.2015.04.002

Schuman MC, Barthel K. Baldwin IT. 2012. Herbivory-induced volatiles function as defenses increasing fitness of the native plant nicotiana attenuata in nature. *eLife* **1**:e00007. doi: 10.7554/eLife.00007

Schäfer M, Brütting C, Gase K, Reichelt M, Baldwin I, Meldau S. 2013. 'Real time' genetic manipulation: A new tool for ecological field studies. *The Plant Journal* **76**:506–518. doi: 10.1111/tpj.12301

Scogings PF, Hjältén J, Skarpe C. 2011. Secondary metabolites and nutrients of woody plants in relation to browsing intensity in african savannas. *Oecologia* **167**:1063–1073. doi: 10.1007/s00442-011-2042-9

Sessions LA, Kelly D. 2001. Heterogeneity in vertebrate and invertebrate herbivory and its consequences for new zealand mistletoes. *Austral Ecology* **26**:571–581. doi: 10.1046/j.1442-9993.2001.01129.x

Shimoda K, Kimura K, Kanzaki M, Yoda K. 1994. The regeneration of pioneer tree species under browsing pressure of sika deer in an evergreen oak forest. *Ecological Research* **9**:85–92. doi: 10.1007/BF02347245

Sime KR, Baldwin IT. 2003. Opportunistic out-crossing in nicotiana attenuata (solanaceae), a predominantly self-fertilizing native tobacco. *BMC Ecology* **3**:6. doi: 10.1186/1472-6785-3-6

Smith AM, Zeeman SC. 2006. Quantification of starch in plant tissues. *Nature Protocols* **1**:1342–1345. doi: 10.1038/nprot.2006.232
Smith WHM, Griffiths DW, Boag B. 1991. Overwinter variation in glucosinolate content of green tissue of rape (brassica napus) in response to grazing by wild rabbit. *Journal of the Science of Food and Agriculture* **56**:511–521. doi: 10.1002/jsfa.2740560411

Steppuhn A, Gase K, Krock B, Halitschke R, Baldwin IT. 2004. Nicotine's defensive function in nature. *PLoS Biology* **2**:e217. doi: 10.1371/journal.pbio.0020217

Steppuhn A, Baldwin IT. 2007. Resistance management in a native plant: Nicotine prevents herbivores from compensating for plant protease inhibitors. *Ecology Letters* **10**:499–511. doi: 10.1111/j.1461-0248.2007.01045.x

Steppuhn A, Schuman MC, Baldwin IT. 2008. Silencing jasmonate signalling and jasmonate-mediated defences reveals different survival strategies between two nicotiana attenuata accessions. *Molecular Ecology* **17**:3717–3732. doi: 10.1111/j.1365-294X.2008.03862.x

Stitz M, Baldwin IT, Gaquerel E. 2011. Diverting the flux of the JA pathway in Nicotiana attenuata compromises the plant's defense metabolism and fitness in nature and glasshouse. *PLoS One* **6**:e25925. doi: 10.1371/journal.pone.0025925

Strauss SY. 1991. Direct, indirect, and cumulative effects of three native herbivores on a shared host plant. *Ecology* **72**:543–558. doi: 10.2307/2937195

Takada M, Asada M, Miyashita T. 2001. Regional differences in the morphology of a shrub damnacanthus indicus: An induced resistance to deer herbivory? *Ecological Research* **16**:809–813. doi: 10.1046/j.1440-1703.2001.00436.x

Thaler JS, Farag MA, Paré PW, Dicke M. 2002. Jasmonate-deficient plants have reduced direct and indirect defences against herbivores. *Ecology Letters* **5**:764–774. doi: 10.1046/j.1461-0248.2002.00388.x

Thaler JS. 1999. Induced resistance in agricultural crops: Effects of jasmonic acid on herbivory and yield in tomato plants. *Environmental Entomology* **28**:30–37. doi: 10.1093/ee/28.1.30

Van Bael SA, Aiello A, Valderrama A, Medianero E, Samaniego M, Wright SJ. 1999. General herbivore outbreak following an el nino-related drought in a lowland panamanian forest. *Journal of Tropical Ecology* **20**:625–633. doi: 10.1017/S0266467404001725

Van Dam NM, Baldwin IT. 2001. Competition mediates costs of jasmonate-induced defences, nitrogen acquisition and transgenerational plasticity in nicotiana attenuata. *Functional Ecology* **15**:406–415. doi: 10.1046/j.1365-2435.2001.00533.x

Veblen TT, Mermoz M, Martin C, Ramilo E. 1989. Effects of exotic deer on forest regeneration and composition in northern patagonia. *The Journal of Applied Ecology* **26**:711–724. doi: 10.2307/2404095

Velterop JS, Vos F. 2001. A rapid and inexpensive microplate assay for the enzymatic determination of glucose, fructose, sucrose, l-malate and citrate in tomato (lycopersicon esculentum) extracts and in orange juice. *Phytochemical Analysis* **12**:299–304. doi: 10.1002/pca.598

Venables WN, Ripley BD. 2002. *Modern Applied Statistics with S (4th ed)*. New York: Springer. doi: 10.1007/978-0-387-21706-2

Vila B, Guibal F. 2001. Assessing browsing influence in forest understory using dendrochronology on Haida Gwaii archipelago (British Columbia, Canada). *Dendrochronologia* **19**:139–151.

Waller DM, Alverson WS. 1997. The white-tailed deer: A keystone herbivore. *Wildlife Society Bulletin*:217–226.

Warner PJ, Cushman HJ. 2002. Influence of herbivores on a perennial plant: Variation with life history stage and herbivore species. *Oecologia* **132**:77–85. doi: 10.1007/s00442-002-0955-z

White PS. 1988. Prickle distribution in aralia spinosa (araliaceae). *American Journal of Botany* **75**:282–285. doi: 10.2307/2443895

Wilson SL, Kerley GIH. 2003. The effect of plant spinescence on the foraging efficiency of bushbuck and boergoats: Browsers of similar body size. *Journal of Arid Environments* **55**:150–158. doi: 10.1016/S0140-1963(02)00254-9

Wink M. 1985. Chemische verteidigung der lupinen: Zur biologischen bedeutung der chinolizidinalkaloide. *Plant Systematics and Evolution* **150**:65–81. doi: 10.1007/BF00985568

Wink M. 1988. Plant breeding: Importance of plant secondary metabolites for protection against pathogens and herbivores. *Theoretical and Applied Genetics* **75**:225–233. doi: 10.1007/BF00303957

Wink M. 1993. Allelochemical properties or the raison d'etre of alkaloids. In: *The Alkaloids: Chemistry and Pharmacology*. Elsevier. p 1–118. doi: 10.1016/S0099-9598(08)60134-0

Wu J, Baldwin IT. 2009. Herbivory-induced signalling in plants: Perception and action. *Plant, Cell & Environment* **32**:1161–1174. doi: 10.1111/j.1365-3040.2009.01943.x

Wu J, Baldwin IT. 2010. New insights into plant responses to the attack from insect herbivores. *Annual Review of Genetics* **44**:1–24. doi: 10.1146/annurev-genet-102209-163500

Xu S, Zhou W, Baldwin IT. 2015a. The rapidly evolving associations among herbivore associated elicitor-induced phytohormones in nicotiana. *Plant Signaling & Behavior* **10**:e1035850. doi: 10.1080/15592324.2015.1035850

Xu S, Zhou W, Pottinger S, Baldwin IT. 2015b. Herbivore associated elicitor-induced defences are highly specific among closely related nicotiana species. *BMC Plant Biology* **15**:2. doi: 10.1186/s12870-014-0406-0

Yang D-L, Yao J, Mei C-S, Tong X-H, Zeng L-J, Li Q, Xiao L-T, Sun T.-p., Li J, Deng X-W, Lee CM, Thomashow MF, Yang Y, He Z, He SY. 2012. Plant hormone jasmonate prioritizes defense over growth by interfering with gibberellin signaling cascade. *Proceedings of the National Academy of Sciences of the United States of America* **109**:E1192–E1200. doi: 10.1073/pnas.1201616109

Yang Y, Shah J, Klessig DF. 1997. Signal perception and transduction in plant defense responses. *Genes & Development* **11**:1621–1639. doi: 10.1101/gad.11.13.1621
Young TP, Okello BD. 1998. Relaxation of an induced defense after exclusion of herbivores: Spines on acacia drepanolobium. *Oecologia* **115**:508–513. doi: 10.1007/s004420050548

Young TP, Stanton ML, Christian CE. 2003. Effects of natural and simulated herbivory on spine lengths of acacia drepanolobium in kenya. *Oikos* **101**:171–179. doi: 10.1034/j.1600-0706.2003.12067.x

Zamora R, Gómez JM, Hódar JA, Castro J, García D. 2001. Effect of browsing by ungulates on sapling growth of scots pine in a mediterranean environment: Consequences for forest regeneration. *Forest Ecology and Management* **144**:33–42. doi: 10.1016/S0378-1127(00)00362-5

Zavala JA, Baldwin IT. 2006. Jasmonic acid signalling and herbivore resistance traits constrain regrowth after herbivore attack in nicotiana attenuata. *Plant, Cell & Environment* **29**:1751–1760. doi: 10.1111/j.1365-3040.2006.01551.x

Zhou G, Qi J, Ren N, Cheng J, Erb M, Mao B, Lou Y. 2009. Silencing oshi-lox makes rice more susceptible to chewing herbivores, but enhances resistance to a phloem feeder. *The Plant Journal* **60**:638–648. doi: 10.1111/j.1365-313X.2009.03988.x

Züst T, Heichinger C, Grossniklaus U, Harrington R, Kliebenstein DJ, Turnbull LA. 2012. Natural enemies drive geographic variation in plant defenses. *Science* **338**:116–119. doi: 10.1126/science.1226397

Züst T, Agrawal AA. 2016. Mechanisms and evolution of plant resistance to aphids. *Nature Plants* **2**:15206. doi: 10.1038/nplants.2015.206