Xosé López-Goldar1,2,3,4*
Caterina Villari2,5
Pierluigi Bonello2
Anna Karin Borg-Karlson3
Delphine Grivet4,6
Rafael Zas1
Luís Sampedro1- 1Misión Biológica de Galicia, Consejo Superior de Investigaciones Científicas, Pontevedra, Spain
- 2Department of Plant Pathology, The Ohio State University, Columbus, OH, United States
- 3Ecological Chemistry Group, Department of Chemistry, Royal Institute of Technology, Stockholm, Sweden
- 4Department of Forest Ecology and Genetics, Forest Research Centre, INIA, Madrid, Spain
- 5Daniel B. Warnell School of Forestry and Natural Resources, University of Georgia, Athens, GA, United States
- 6Sustainable Forest Management Research Institute, INIA-University of Valladolid, Palencia, Spain
Resistance to herbivores and pathogens is considered a key plant trait with strong adaptive value in trees, usually involving high concentrations of a diverse array of plant secondary metabolites (PSM). Intraspecific genetic variation and plasticity of PSM are widely known. However, their ecology and evolution are unclear, and even the implication of PSM as traits that provide direct effective resistance against herbivores is currently questioned. We used control and methyl jasmonate (MJ) induced clonal copies of genotypes within families from ten populations of the main distribution range of maritime pine to exhaustively characterize the constitutive and induced profile and concentration of PSM in the stem phloem, and to measure insect herbivory damage as a proxy of resistance. Then, we explored whether genetic variation in resistance to herbivory may be predicted by the constitutive concentration of PSM, and the role of its inducibility to predict the increase in resistance once the plant is induced. We found large and structured genetic variation among populations but not between families within populations in resistance to herbivory. The MJ-induction treatment strongly increased resistance to the weevil in the species, and the genetic variation in the inducibility of resistance was significantly structured among populations, with greater inducibility in the Atlantic populations. Genetic variation in resistance was largely explained by the multivariate concentration and profile of PSM at the genotypic level, rather than by bivariate correlations with individual PSM, after accounting for genetic relatedness among genotypes. While the constitutive concentration of the PSM blend did not show a clear pattern of resistance to herbivory, specific changes in the chemical profile and the increase in concentration of the PSM blend after MJ induction were related to increased resistance. To date, this is the first example of a comprehensive and rigorous approach in which inducibility of PSM in trees and its implication in resistance was analyzed excluding spurious associations due to genetic relatedness, often overlooked in intraspecific studies. Here we provide evidences that multivariate analyses of PSM, rather than bivariate correlations, provide more realistic information about the potentially causal relationships between PSM and resistance to herbivory in pine trees.
Introduction
Resistance to herbivores and pathogens has been widely recognized as a key plant trait with strong adaptive value (Futuyma and Agrawal, 2009), particularly at early stages of development (Goodger et al., 2013). In long-lived plants with high apparency, like trees, resistance against pests and pathogens usually relies on high concentrations of a diverse array of plant secondary metabolites (PSM) (Feeny, 1976; Wiggins et al., 2016). Such chemical cocktail generally produces dose dependent effects, where higher concentrations in specific plant tissues targeted by herbivores results in reduced enemy performance and/or plant damage (Zhao et al., 2011b). These chemical defenses are expensive to produce and maintain (Gershenzon, 1994; Sampedro et al., 2011a). All these features are certainly present in conifers, whose higher concentrations of terpenes and phenolics are commonly associated with increased resistance against insect herbivores (Phillips and Croteau, 1999; Witzell and Martín, 2008; Schiebe et al., 2012).
Within conifer species, variation in defensive PSM is usually difficult to explain based exclusively on an arms race between trees and their biotic aggressors because many other factors often affect plant–organism interactions (Moore et al., 2013). For example, heterogeneity in the biotic and abiotic environment (Lombardero et al., 2000; Sampedro et al., 2010), the large, repetitive, and idiosyncratic structure of conifer genomes (Theis and Lerdau, 2003), trade-offs with other life functions (Sampedro et al., 2011a; Moreira et al., 2014) and other evolutionary forces, such as gene flow and genetic drift (Kremer et al., 2012), can all be involved in maintaining large intraspecific variation in PSM as well as in resistance. Phylogenetic analyses has been widely used in recent years (Rasmann and Agrawal, 2011; Pellissier et al., 2012; Moreira et al., 2014; Carrillo-Gavilán et al., 2015) to control for the statistical dependence among related species in macroecology (Freckleton et al., 2002), leading to the identification of imprinted historical events in contemporary observed traits (Pagel, 1999). At the intraspecific level, demographic processes, habitat fragmentation, gene flow and local adaptation can also greatly influence the variation in a species’ traits among and within populations (Savolainen et al., 2007; De-Lucas et al., 2009), making necessary to account for the genetic structure of populations and the relative kinship among genotypes within families when exploring trait-trait or trait-environment associations (Yu et al., 2006).
Evidence is accumulating that PSM are highly plastic, as most plant species respond with large changes in the concentration and profile of PSM to the perception of a biotic threat, such as a pathogen or a herbivore (Karban et al., 1999; Heil, 2009). This sort of biochemical plasticity, referred to as inducibility, may be a very important ecological and evolutionary plant trait, as it is a source of functional phenotypic variation in resistance and thus, potentially, a trait under selection (Agrawal et al., 2015). Considering a single plant genotype, inducibility is usually estimated as the difference between the induced vs. the constitutive concentration of PSM. In large plants, estimating inducibility of a PSM may require a repeated measures design within the same subject, i.e., before and after challenge with real insect herbivory (Morris et al., 2006). However, a common undesired side-effect is that sampling plant tissues may itself be an inducer of plant defenses, which would alter the estimation of inducibility (see Moreira et al., 2012, for example). Even more challenging can be the establishment of any relationship between inducibility of PSM and the effective resistance against real herbivory. This requires knowing the concentration of PSM before and after exposure of both constitutive and induced phenotypes to real herbivory when measuring plant damage or insect performance. For small plants and to evaluate plant traits that require destructive sampling, this problem can only be solved using clonal replicates of the same genotype.
At the intraspecific level, while genetic variation and plasticity of PSM are widely recognized, the ecology and evolution of PSM are far from being clear, and the role itself of PSM as traits providing direct effective resistance against herbivores is currently under debate (Carmona et al., 2011; Cárdenas et al., 2014). Despite the common finding of positive associations between concentration of PSM and resistance to specific herbivores (Borg-Karlson et al., 2006; Erbilgin et al., 2006; Kannaste et al., 2009; Zhao et al., 2011b), cause-effect studies demonstrating the direct role of particular PSM or blends of PSM on resistance are much less frequent (Agrawal, 2005; Iason et al., 2011).
In this study we tested the hypothesis that intraspecific variation in effective resistance of maritime pine (Pinus pinaster) to a key insect herbivore, the pine weevil (Hylobius abietis), may be predicted by variation in inducibility of PSM. We used clonally replicated genotypes with known family structure from 10 populations representative of the main distribution range of maritime pine. We performed an exhaustive characterization of constitutive concentrations and profiles of PSM in the stem phloem and their associated induced changes after application of methyl jasmonate (MJ), a plant hormone involved in signaling of chewing herbivore damage (Sampedro et al., 2011b; Schiebe et al., 2012; Zas et al., 2014). Both constitutive and MJ-induced clonal copies were exposed to insect herbivory to measure plant damage as a proxy of resistance (i.e., the greater the damage, the lower the resistance). Then we explored to what extent resistance to herbivory may be predicted by the constitutive concentration of PSM in the stem, and the role of its inducibility in predicting the increase in resistance once the plant is induced. Since the non-independence due to relatedness among genotypes, families and populations may affect the relationships between PSM and resistance, we accounted for the population structure and the relative kinship among individuals within populations.
Materials and Methods
Study System
We used P. pinaster as a model tree system, a Mediterranean conifer species with a fragmented distribution displaying large and structured among-population genetic variation (Vendramin et al., 1998; Bucci et al., 2007; Jaramillo-Correa et al., 2015). Populations largely differ in life history traits (Tapias et al., 2004; Grivet et al., 2013), including those putatively related to resistance to biotic stress and PSM (Sampedro et al., 2011b; Zas et al., 2015).
The pine weevil H. abietis is a bark-chewing insect herbivore with Palaearctic distribution, and is specialized in young conifers, feeding on several species (Leather et al., 1999). It harbors low neutral genetic variation across its large distribution range (Conord et al., 2006). The pine weevil is a harmful forest pest that causes severe damage on young conifer seedlings, and has a strong impact on plant fitness and early survival (Langström and Day, 2004; Zas et al., 2006, 2014). Host resistance has been showed to be related to plant survival and to be heritable (Zas et al., 2005). The pine weevil became a model insect for studies of resistance in pine trees (see for instance Nordlander et al., 2003; Wainhouse et al., 2005; Björklund, 2007; Nordlander et al., 2011; Sampedro et al., 2011b; Zas et al., 2014, 2017). The success in the regeneration of conifer forests is highly contingent on this insect (Langström and Day, 2004), which causes important economic losses in large part of Europe (Leather et al., 1999). Larval stages feed on stump roots of recently cut trees (Nordenhem and Nordlander, 1994), and newly emerged adults migrate in massive numbers to fresh clear-cuts in spring (Solbreck and Gyldberg, 1979; Tan et al., 2011). Generation times may vary from 1 to 3 years depending on the location in Europe (Day et al., 2004; Wainhouse et al., 2014). The longer generation times observed at higher latitudes, and the continuous masses of conifer forests, particularly in North Europe, strongly facilitate pine weevil dissemination, establishment, and subsequent damage among areas with clear-cuts for at least two up to four seasons after harvesting (Örlander et al., 1997). Apart from control and management measures to reduce the impact of the insect on plant fitness and survival (Leather et al., 1999; Nordlander et al., 2011), it is important to know the defensive potential in conifers to deal with this pernicious pest.
Plant and Herbivore Material
Ten maritime pine populations belonging to the collection ‘CLONAPIN Bank 1’ were selected along a latitudinal gradient across the natural distribution range of the species (Figure 1), from the Atlantic France to North Africa. Briefly, open-pollinated seeds were collected from five mother trees randomly selected in natural stands within each of the ten populations, resulting in fifty half-sib families (i.e., known ‘mother’ and unknown ‘father’). A separation of at least 50 m between selected mother trees was performed to minimize inbreeding (González-Martínez et al., 2006). Five seeds from each mother tree were sown and cultured in the facilities of SERIDA (Asturias, Spain), producing the complete clonal ortet collection, which comprised 250 genotypes (10 populations × 5 families × 5 half-sib genotypes). These genotypes were then clonally propagated by recurrent early cutting as described in Majada et al. (2011). A total of 205 out of the 250 original genotypes (belonging to 48 families out of the original 50) were eventually included due to low germination rate and/or rooting efficiency.
FIGURE 1. Natural distribution range of maritime pine (shaded area) and location of the ten studied populations (colored circles). Circle color represents the genetic membership of the populations based on STRUCTURE (see Supplementary Methods SM2). Populations were coded with four capital letters (in brackets). Modified from EUFORGEN (2009).
The cuttings were grown at SERIDA Agricultural Station (Grado, Asturias, Spain) in 200 mL containers containing peat and perlite (70:30 v/v), then transported to Norfor Forest Nursery (Figueirido, Pontevedra, Spain) and transplanted to 2 L pots. In September 2012, when pine cuttings were around 2-years old, the plants were transferred to an environmentally controlled greenhouse in the MBG-CSIC (Pontevedra, Spain). Plants were fertilized with a slow release fertilizer (Granum, Soaga SL, Vilanova de Arousa, Spain, NPK 11-22-9), supplemented with a foliar fertilizer (Fertimón, Progando SL, Galicia, Spain), and watered twice a week. Pines were allowed to acclimatize until the start of the experiment by December 2012. During this period, mean temperatures and relative humidity were recorded during the day (25.6 ± 0.06°C and 57.2 ± 0.08%, respectively) and during the night (18.2 ± 0.05°C and 69.0 ± 0.05%, respectively).
Adult weevils were collected at a clear-felled pine forest in Monte Castrove (Pontevedra, Spain, 42°2702300N, 8°4301900W) using Nordlander traps (Nordlander, 1987) and maintained in culture chambers at 10°C several weeks until the start of the herbivory bioassays (see below). During that period, the pine weevils were fed weekly with fresh twigs of maritime pine.
Experimental Design
Plants were assigned to treatments following a factorial split-plot design replicated in 5 blocks, with population as the whole plot factor (ten populations) and the factorial combination of family (3–5 families per population) and the MJ induction treatment (2 levels: control and induced) as the split factor. Two clonal replicates of each genotype from each family were included in each of the two MJ-induction treatments. One of these copies was used for chemical analyses and the other for evaluation of the effective resistance against the pine weevil. The experiment comprised 814 plants from 10 populations, with 3–5 families per population, 5 half-sibs (genotypes) per family, and 4 clonal replicates per genotype (two of them, control and MJ, used for chemical analysis; and two, control and MJ, used in the weevil bioassays). A detailed graphical scheme of the experimental design can be found in the Supplementary Figure SM1.
Induction Treatment
Plasticity of PSM in response to biotic stress has been extensively studied in maritime pine. Large changes in the concentration and chemical profile of total or specific PSM using either the exogenous application of MJ (Moreira et al., 2009), real herbivory (Moreira et al., 2013a), or both (Sampedro et al., 2011b; Moreira et al., 2012) have been reported in this and other pine species (Moreira et al., 2014; Carrillo-Gavilán et al., 2015). In December 2012, plants allotted to induction were transferred to a separate chamber prior to treatment application to avoid cross-contamination with control plants. Then, MJ-induced plants were sprayed with a suspension of 25 mM of 95% MJ (Sigma–Aldrich, #39270-7) in deionized water with 2.5% ethanol (v/v). Control plants were treated just with the carrier solution. Both solutions were sprayed over the aboveground part of each plant to runoff (∼10 mL per plant). MJ concentration and time elapsed between MJ application and sampling (1 month, see below) were based on previous studies (Moreira et al., 2013b; Zas et al., 2014). This lapse of time allowed plants to develop induced and long-lasting anatomical structures (i.e., resin ducts and polyphenolic parenchyma cells) that are involved in the production of induced defensive PSM (Moreira et al., 2015).
Plant Sampling, Measurements and Analysis of Plant Secondary Metabolites
One month after MJ-application, two of the clonal copies of each genotype (one treated with MJ and one control) were destructively harvested by blocks, cutting the stem aboveground. Needles were separated from the stem and two 2-cm-long stem subsamples were immediately collected from its middle part, flash-frozen in liquid nitrogen, and stored at -80°C. We then analyzed the composition and concentration of PSM in 3 out of the 5 blocks (260 plants). Briefly, terpenoids were extracted from one of the frozen stem subsamples with HPLC-grade hexane (HiPersolv Chromanorm #83992.320) and separated, identified and quantified by gas chromatography-mass spectrometry (GC-MS) and GC-flame ionization detection (GC-FID). Phenolics were extracted from the other stem subsample with HPLC-grade methanol (HiPersolv Chromanorm #152506X) and separated, identified and quantified by ultra high performance liquid chromatography-mass spectrometry (UHPLC-MS) and UHPLC-diode array detection (UHPLC-DAD). A detailed description of the analytical procedures and data processing prior to statistical analyses can be found in the Supplementary Methods SM1).
We identified 118 PSM that were classified into eight chemical groups [monoterpenes, sesquiterpenes, diterpenes, flavonoids, hydroxycinnamic acids (HCAs), lignans, eugenols, and fatty acids]. For the purpose of the current study, we only considered those compounds present in at least 5 individuals of the control or MJ-induced plants for each population. Ninety-eight chemical compounds were finally included in the eight chemical groups (Supplementary Tables SM1, SM2).
Herbivory Bioassays
At the time plants not subjected to herbivory were harvested for chemical analyses, the remaining two clonal replicates (one control and one treated with MJ) were used for evaluating the resistance of each genotype to the pine weevil, following López-Goldar et al. (2016), with modifications. Briefly, a transparent plastic tube was fitted around each plant and a randomly selected pair of weevils (one male and one female) were put on each plant and allowed to feed along the entire plant for 48 h. A fine-mesh was fitted at the top of the tube to prevent weevil escape. Weevils were starved in a Petri dish with moist filter paper (48 h at 15°C in dark) and pre-weighed prior the herbivory treatment. After 48 h of exposure to herbivory along the entire plant, weevils were removed and the damage caused by the weevils in the stem of each plant was measured at the nearest mm2 using calibrated area templates.
Genetic Variation in Resistance to Herbivory
Variation in the damage caused by H. abietis in the stem was analyzed by fitting a mixed model using the PROC MIXED procedure in SAS v9.4 (Littell et al., 2006), with the induction treatment (MJ), population (P), family within population [F(P)] and MJ × P and MJ × F(P) interactions as fixed factors. The MJ × P and MJ × F(P) interactions are informative of the variation in the responses to the induction treatment among populations and families within populations, respectively. Block (B), the B × P interaction and the genotype [MJ × F(P) × B] were considered random factors. The B × P interaction was considered a random factor in order to analyze the main effect of P with the appropriate error term in the split-plot design (Littell et al., 2006). The genotype was also included as a random factor to account for the dependence between the two clonal replicates subjected to herbivory within each block. Damage was square root-transformed to achieve normality. Heterogeneous residual variance models across P were fitted. Weevil weight, as a proxy of weevil size, and stem diameter of the plant, both known to affect feeding rate (Wainhouse et al., 2004; Wainhouse et al., 2005), were included as covariates because they improved model fit. The same model was fit to explore the genetic variation among genetic groups of populations (GP) and within GP (families within each genetic group). Genetic groups of populations were identified and characterized elsewhere using molecular markers (Jaramillo-Correa et al., 2015; Serra-Varela et al., 2015), and corresponded to: Atlantic France (PLEU and MIMI), Atlantic Iberia (PTOV, CDVO, ARMY and SCRI), Central Iberia (COCA, ASPE), Southwestern Iberia (ORIA) and Morocco (TAMR) (Figure 1).
In order to quantify the relative contribution of each factor to the genetic variation among populations and families within populations in resistance to the pine weevil, the previous mixed model was run again assuming all factors as random factors. Variance components of each factor were estimated by restricted maximum likelihood.
Relationships Between PSM and Resistance to Herbivory
To explore the causal relationships between resistance and PSM in the stem of P. pinaster, we performed pairwise correlations and multiple regression analyses between the weevil damage and PSM concentration at the genotypic level. Exploring the implication of PSM in resistance to herbivory in individuals within populations across the main distribution range will allow us to predict which PSM or blend of PSM within a genotype might be the optimal to determine the highest resistance to herbivory in the entire species. The clonal structure of our pine collection allowed for this genotypic-based analysis.
Because genotypes belong to a hierarchical collection of families within populations and to different populations in maritime pine, we accounted for the non-independence among genotypes within families and populations to avoid genetic effects due to relatedness in the correlation analyses (see Supplementary Methods SM2). For that, we included the population structure (Q) and kinship (K) matrices at the genotypic level using a mixed-model approach as described in Yu et al. (2006). Briefly, Q and K matrices were constructed based on 126 SNPs: for the Q matrix, we performed a Bayesian cluster analysis in STRUCTURE v 2.3.4 (Pritchard et al., 2000), whereas for the K matrix we used SPAGeDi (Hardy and Vekemans, 2002) based upon Loiselle et al. (1995) kinship coefficients. The Q and K matrices were then incorporated to the mixed models in order to account for the multiple level of relatedness among individuals (see Supplementary Methods SM2).
In addition, we fitted a mixed model for constitutive (C) and for the inducibility (MJ – C) of each total and individual PSM and weevil damage. For each defensive mode, population (P) and family within population [F(P)] were considered as fixed factors. Block (B) and the B × P interaction were considered random factors as above. Q and K matrices were incorporated into each mixed model to account for the genetic relatedness among genotypes. The predicted values of the variables from each plant in each mixed model were then used for the pairwise correlations and multiple regression analyses. Because PSM data were available in 3 out of 5 blocks, pairwise correlations and multiple regression analyses with weevil damage were performed only in 102 genotypes.Given that we were not able to extract the predicted values for all 102 genotypes after implementing Q and K matrices, and also not all clonal copies of the same genotype were available, sample sizes were slightly lower than the original 102 genotypes for each defensive state.
Two types of evaluations were performed to explore the relationships between PSM and weevil damage in the stem of P. pinaster. First, to observe the basal resistance without interference of the induction treatment, we explored the relationships between constitutive (C) PSM and damage by H. abietis. Second, we explored the relationships between the estimated inducibility of each total and individual PSM (MJ – C) with the inducibility of resistance to the pine weevil after MJ application. For each type of analysis, Pearson correlations between the concentration of each total and individual PSM and weevil damage were performed using PROC CORR in SAS v9.4. To estimate all the pairwise correlations without inflating Type I error due to multiple testing, p-value adjustments were performed using the False Discovery Rate (FDR) for p < 0.05 (Benjamini and Hochberg, 1995). Also, for the same two evaluations (constitutive and inducibility), multiple regression analyses were performed with PROC REG in SAS v9.4 to explore to what extent the PSM blend predicted herbivory resistance using the stepwise model selection method on the 98 individual PSM. We applied a more stringent threshold than default parameters (α = 0.05; default, α = 0.15) when including and excluding variables during the stepwise method to reduce type I error due to the large number of PSM used. Predictor variables were previously standardized (mean = 0, variance = 1) and influential observations (measured with Cook’s D) that did not meet the requirements were excluded from the analysis. Normality of the residuals was achieved using square-root transformation of the response variable (weevil damage) when necessary.
Results
Sources of Variation in Resistance to Herbivory by Hylobius abietis in Pinus pinaster
Resistance to weevil herbivory largely and significantly varied among populations but not between families within populations (companion table in Figure 2). The highest constitutive resistance was observed for the Spanish Atlantic population PTOV, whereas the most constitutively susceptible was the French Atlantic population MIMI (Figure 3). Large genetic variation was also found among genetic groups of populations (F4,140 = 4.73, P = 0.001), revealing similar patterns of constitutive resistance among populations of the same genetic group. The highest constitutive resistance was found for the Atlantic Spanish group (61.5 ± 7.3 mm2 of damage), whereas the Atlantic French group was the most susceptible (94.9 ± 12.3 mm2 of damage).
FIGURE 2. Summary of the mixed model testing the effects of the induction treatment with methyl-jasmonate (MJ), population (P), family within population [F (P)] and the MJ × P and MJ × F(P) interactions on the early resistance to herbivory by the pine weevil on plants from ten maritime pine populations. Variance components for each effect and the residual are shown in the companion pie chart. Weevil weight and plant diameter were included as covariates. Significant p-values (p < 0.05) are highlighted in bold. A total sample size of 205 genotypes was used.
FIGURE 3. Intraspecific genetic variation in constitutive resistance (white bars) and MJ-induced resistance (gray bars) to the pine weevil across 10 maritime pine populations representing its main distribution range, grouped by the main genetic groups in the species. The inset panel shows the overall effect of the MJ-induction on the resistance to weevil herbivory in the species. Bars represent the least square mean ± SE (N = 8–25 plants for each population; N = 197–205 plants for each induction treatment in the small panel). A total sample size of 205 genotypes was used.
The MJ-induction treatment strongly increased resistance to herbivory in the stem of all populations (companion table in Figure 2), reducing the overall damage caused by the pine weevil to less than half compared to non-induced plants (Figure 3, inset panel). The French population PLEU showed the highest induced resistance, whereas COCA emerged as the less resistant after MJ application (Figure 3). The inducibility of resistance (i.e., the plasticity to MJ application) was similar for all populations and families [no significant MJ × P and MJ × F(P) interactions, companion table in Figure 2], but differed among genetic groups (MJ × GP interaction, F4,140 = 3.36, P = 0.012). Interestingly, greater inducibility of resistance after induction was found in both Atlantic genetic groups (71 and 67% in Atlantic France and Atlantic Spain, respectively) than in the other genetic groups (Central Spain, 37%; Southeastern Spain, 34%; and Morocco, 45%).
The main effects in our experimental design, namely MJ-induction, population, family and their interactions explained 36% of the total variance in resistance. MJ-induction was by far the most relevant factor, absorbing most of the explained variation in resistance among the studied effects (>80% of the explained variance), with the contribution of all remaining factors being comparatively much lower (Figure 2).
Concentration of PSM in the Stem as Predictor of Damage by the Pine Weevil at the Genotypic Level
We performed pairwise correlations in order to explore the relationship between concentration of total and individual PSM in the stem with weevil damage, while accounting for population structure and genetic relatedness among individuals. Total and individual PSM did not show any significant negative relationship with weevil damage at the constitutive level (Table 1 and Supplementary Table SM3). When we explored the inducibility, total flavonoids and total lignans showed significant negative relationships with weevil damage (Table 1). Consistently, when considering individual PSM, significant negative correlations were found for the inducibility of two lignan hexoside derivatives, and marginally significant for three lignan xyloside derivatives and taxifolin derivative 1, after FDR correction (Supplementary Table SM3). As an exception, eugenol was the only PSM that was positively related with weevil damage at the constitutive level (Supplementary Table SM3).
TABLE 1. Pearson correlations between damage by the pine weevil, as a proxy of resistance, and the concentration of total plant secondary metabolites (PSM) in the stem phloem of 102 genotypes from 10 populations, representing the main distribution range of maritime pine.
We used multiple regression analyses looking for the relative predictive value of individual PSM in the relationship between PSM and damage by H. abietis in the stem, while accounting for population structure and relative kinship among genotypes in P. pinaster. At the constitutive level, 18 individual PSM (10 terpenes and 8 phenolics; Table 2) significantly predicted 70% of the total variation in constitutive weevil damage (F14,66 = 18.97, P < 0.001). Six PSM (4 terpenes and 2 phenolics) were negatively related with weevil damage, as shown by their negative standardized regression coefficients (Table 2). With regard to inducibility, 12 PSM (8 terpenes, 3 phenolics and 1 fatty acid; Table 2) significantly predicted 68% of the total variance of weevil damage (F11,83 = 17.32, P < 0.001). Eight PSM (6 terpenes, 2 phenolics) showed standardized regression coefficients that were negatively related to the variation in weevil damage after MJ-induction (Table 2).
TABLE 2. Summary of the multiple regression analysis explaining the resistance to the pine weevil using the concentration of individual PSM as predictor variables in constitutive (upper part of the table) and in inducibility (lower part of the table) defensive modes in 102 genotypes from 10 maritime pine populations.
Discussion
In our study, large and structured genetic variation in resistance to herbivory by H. abietis in the stem was found among P. pinaster populations, as well as significant differences among genetic groups of populations in resistance after MJ-induction. In addition, the resistance patterns were strongly predicted by the multivariate concentration and profile of PSM of genotypes after accounting for the genetic relatedness, rather than by bivariate correlations with weevil damage. Furthermore, differences in the resistance patterns between the constitutive and inducibility strategies were observed. While a high proportion of the variance for weevil damage was explained by the constitutive PSM, the directionality of the pattern (i.e., resistance or susceptibility) was unclear. Conversely, inducibility of PSM, besides explaining also a large proportion of the variance in weevil damage, predicted greater inducibility of resistance mainly due to relevant plastic and specific changes in the profile and increases in the concentration of the chemical compounds after induction.
Inter- but Not Intrapopulation Genetic Variation and Increased Resistance to Herbivory After Induction
Our results show that P. pinaster populations are variable in resistance against the bark chewing herbivore H. abietis, but we found no evidence of variation between families within populations. Using the same clonal collection, Elvira-Recuenco et al. (2014) reported similar results for intraspecific variation in constitutive resistance against the fungal pathogen Fusarium circinatum. Atlantic populations showed greater constitutive resistance (with MIMI population as the exception in our study) than those from the Mediterranean Basin in both studies. This was, perhaps, an unexpected outcome, given the wide taxonomic dissimilarity of the two enemies. F. circinatum is considered an alien invasive pathogen in Europe (Wingfield et al., 2008) where, to date, it has been detected in nature in only one adult P. pinaster plantation in the northern Iberian Peninsula (Iturritxa et al., 2013). These results suggest that similar functional defensive traits may be effective against both biotic enemies, and that populations from growth-prone environments are better defended against these biotic agents at early stages. Similarly, we found extensive genetic variation among genetic groups of populations. This suggests that the differences in resistance levels among individual populations might be due to greater similarity among genetically closely related populations than distantly related ones. We did not find, however, genetic variation among families within populations for resistance, which is contrasting with previous results using the same pine species (Zas et al., 2005; Elvira-Recuenco et al., 2014). This finding might be due to the strong effect of the induction treatment, which may have overwhelmed most of the variation between families, or likely to sample size, that might be too low for detecting family variation.
Boosting MJ signaling pathways by exogenous application of MJ to produce a phenotype with increased expression of chemical and physical defenses is a well-known experimental tool in conifers (Heijari et al., 2005; Zhao et al., 2011a; Moreira et al., 2013b; Zas et al., 2014) and was confirmed in this study, with MJ being by far the most relevant factor (>80% of the total explained variance in resistance). The lack of a significant MJ by population interaction suggests that the patterns of resistance to herbivory elicited by MJ are highly regulated and conserved in P. pinaster, and probably derived from common ancestral defensive mechanisms shared in the conifer lineage (Hudgins et al., 2004). Results from previous studies that explored the genetic variation in the inducibility of resistance to insect herbivores after MJ application reported contrasting results. MJ induction reduced herbivory damage in P. radiata (Moreira et al., 2013b) and P. pinaster (Sampedro et al., 2011b), but family genetic variation in the inducibility of resistance was only detected for the folivore Thaumetopoea pityocampa and not for H. abietis (Sampedro et al., 2011b; Moreira et al., 2013b), suggesting that resistance mechanisms would be differentially regulated within populations depending on the herbivore identity and/or the target tissue. Despite that we did not find a MJ × population interaction, we did find it out among genetic groups of populations, indicating that the inducibility of resistance after MJ application is genetically structured in the species, i.e., resistance plasticity is more similar among genetically closely related populations than among distantly related ones. Our results highlight the evolutionary relevance and adaptive value of induced defenses and inducibility as putative key traits explaining variation in resistance. In addition, our findings are in agreement with the predictions regarding the intraspecific patterns of variation in plant defense (Hahn and Maron, 2016), although more profound analyses are required to appropriately address the postulates of the theory, which are out of scope in this study.
Accounting for Biases in Relationships Between PSM and Resistance Across and Within Species
While associations between concentration of PSM and resistance have been widely reported in conifers (Borg-Karlson et al., 2006; Erbilgin et al., 2006; Kannaste et al., 2009), cause-effect studies directly exploring the particular role of individual PSM or PSM blends with resistance have been little explored (Agrawal, 2005; Iason et al., 2011). Moreover, recent studies addressing the implication of PSM in resistance to herbivory across and within species show contrasting results (Moles et al., 2011; Rasmann and Agrawal, 2011; Abdala-Roberts et al., 2016; Moreira et al., 2018), putting the defensive role of PSM under strong criticism (Carmona et al., 2011; Cárdenas et al., 2014). Such controversy might have arisen due to, among other factors, a lack of control for genetic relatedness among species, populations and individuals, and the consideration of broad groups of PSM (e.g., all terpenoids or all phenolic) instead of small chemical groups or individual PSM. First, in order to explore the relationships between PSM and resistance, the genetic relatedness across or within species must be accounted for to avoid spurious trait-trait or trait-environment relationships. In this sense, phylogenetic correction allows to account for the variation in traits due to common ancestry across related species (Freckleton et al., 2002), whereas the implementation of population structure and relative kinship among individuals allow to account for the variation due genetic relatedness among populations and among individuals within families in a single species (Yu et al., 2006). Second, the role of PSM can be better and precisely described when explored in small chemical groups or at the individual scale, because exploring relationships only among broad classes of compounds with resistance might overstate the implication of some individual PSM from those large groups that are not adaptively configured for defense (Hamilton et al., 2001). Furthermore, previous studies that criticize the adaptive value of PSM against resistance to herbivory (Carmona et al., 2011; Cárdenas et al., 2014) explored not only the effect of a small fraction of the whole diversity of chemical compounds (mostly few types and broad groups of phenolics), but also overlooked the inducibility of PSM as a relevant trait associated with resistance. In this sense, calling into question the role of PSM in resistance to herbivory based only on the effect of concrete and broad groups of PSM and when their plasticity was not measured, may obscure the adaptive value of PSM due to biased overgeneralizations. An exhaustive characterization of the constitutive chemical profile and its inducibility of the individuals would be required to maximize the range of metabolites putatively related to resistance and to identify those that might accurately predict the variation in resistance across and within species.
PSM as Resistance Markers Against Herbivory
Plant secondary metabolites have been found to play an essential role in plant resistance to herbivory in many and phylogenetically distant plant species (for instance, see Agrawal et al., 2002; Zhao et al., 2010; Rasmann et al., 2014; Mason et al., 2016); thus, pointing to their adaptive value (Futuyma and Agrawal, 2009; Agrawal, 2010). Developing novel and heritable defensive traits that allow to resist herbivore attacks is of particular relevance in long lived plants like trees (Futuyma and Agrawal, 2009).
In our study, i) we maximized the range of phenotypic variation in the concentration of PSM potentially related to resistance against herbivores by exploring the constitutive allocation to PSM and their inducibility in the stem, ii) minimized possible genetic biases by using clonal replicates, iii) and accounted for relatedness among genotypes within families and populations of P. pinaster during the analyses. Our results showed that inducibility of flavonoids and lignans were strongly related to the inducibility of resistance to the pine weevil when exploring pairwise relationships. Previous studies have suggested that flavonoids and lignans are involved in plant defense (see reviews in Treutter, 2006; Witzell and Martín, 2008), acting as direct defenses against biotic agents (Wallis et al., 2010; Villari et al., 2016) as precursors in cell wall lignification for physical defense (Ralph et al., 2006), and being involved in protection against oxidative stress (Treutter, 2006). Only the concentration of eugenol at the constitutive level was found to be positively related with H. abietis damage (i.e., negatively related to resistance), which contrasts with previous findings where it was found to be an antifeedant for the pine weevil (Borg-Karlson et al., 2006; Eriksson et al., 2008). It is noteworthy that those previous studies were performed under laboratory conditions using purified extracts of eugenol, thus overlooking the interactive effects with other chemical constituents in the PSM blend that would affect insect feeding behavior. An alternative explanation for this finding is that the effect of eugenol was observed to be short-lasting, with decreased anti-feedant properties after several hours of exposure (Borg-Karlson et al., 2006), maybe due to the high volatility of this PSM. It is possible that the very low constitutive concentration of eugenol had an attractant effect rather than repellent in higher concentrations. Surprisingly, we did not detect any relationship between the constitutive concentration of terpene subgroups or individual terpenes, and their inducibility, with resistance to the pine weevil, despite previous results supporting this link (Danielsson et al., 2008; Kannaste et al., 2009; Schiebe et al., 2012; Lundborg et al., 2016). Noteworthy, however, when we explored the whole phenotypic range, including both constitutive and MJ-induced clonal replicates, all total terpene subgroups, total eugenols and total lignans, and many individual terpenes and phenolics were strongly and positively related to resistance against herbivory in the stem (data not shown) as previously found in multiple studies (Porth et al., 2011; Sampedro et al., 2011b; Schiebe et al., 2012; Zas et al., 2014; Lundborg et al., 2016). This result, together with the lack of significant pairwise associations between constitutive PSM (and inducibility) and resistance, suggest that MJ may have artificially forced the pairwise relationships when exploring the whole phenotypic variation. Previous results reporting PSM-resistance relationships upon control and induced plants should thus be interpreted with care. Based on our findings, we suggest that the strong effect of MJ in the reduction of weevil damage is product of a large multivariate induction rather than smaller and punctual individual effects on specific PSM.
In fact, when we explored the predictive power of the whole PSM blend in the stem, comprising 98 chemical compounds of different chemical nature, on resistance to H. abietis damage, we found strong evidence that resistance is largely mediated by the concentration of a specific blend of PSM and also by changes in the chemical profile after induction. On the one hand, constitutive resistance was explained at the multivariate level by terpene volatiles (mono- and sesquiterpenes) and many simple phenolics or precursors of complex molecules (for example, eugenols; Dudareva et al., 2013). This indicates that undamaged plants might initially invest more in simple carbon-based compounds in the stem as a cost-saving strategy when the biotic risk of attack is expected to be low (Karban, 2011). Among these PSM, β–phellandrene + limonene and lignan hexoside 1 have been showed to be strong and effective herbivore repellent and antifeedant, respectively, in conifers and other tree species, with increased resistance found in plants with greater concentrations of these PSM (Schiebe et al., 2012; Villari et al., 2016). While β–phellandrene and limonene coeluted during the chromatographic run in our study, limonene was previously observed to be significantly induced by MJ, whereas β–phellandrene was not affected (Zhao et al., 2010) and seedlings with greater concentrations of limonene in the stem were less attacked by H. abietis (Danielsson et al., 2008). On the other hand, plasticity of PSM revealed an important allocation toward more complex PSM, and changes in the chemical profile of the PSM blend after MJ application. Induction by MJ strongly increased the investment in larger molecules in the stem, such as diterpene resin acids, which are known to be toxic and function as deterrents (Larsson et al., 1986; Tomlin et al., 1996), but also contribute to the sealing of wounded tissue and mechanically impede feeding activity (Ericsson et al., 1988; Phillips and Croteau, 1999). Two phenolics previously found in the pairwise correlations, Unk P5 and lignan hexoside derivative 2, predicted almost half of the variance in increased resistance relative to the whole inducible PSM blend, supporting prior work (Ralph et al., 2006; Villari et al., 2016), and suggesting that these compounds are probably key in resistance, needing thus further investigation.
Despite their general ecological function, the role of fatty acids in conifer resistance is not well understood. We found that the inducibility of oleic acid is associated with susceptibility to pine weevil in maritime pine, confirming that fatty acid composition may affect host suitability for insect pests (Ishangulyyeva et al., 2016). However, due to the important investment in resources that induced defenses usually imply (Cipollini and Heil, 2010), it is possible that the production of induced PSM comes, in part, from the utilization of fatty acids after remobilization and digestion of lipid bodies in the plastids (Bernard-Dagan, 1988).
Interestingly, abietic acid, one of the most well-known diterpene resin acids associated with resistance against insect herbivores (Wagner et al., 1983; Tomlin et al., 1996; Powell and Raffa, 1999), showed an opposite pattern in our study. Recent studies found, however, that gut microbiota of stem chewers and borers can help overcome host defenses, particularly diterpene resin acids, when they are present at relatively low concentrations (Boone et al., 2013; Berasategui et al., 2017). Indeed, such microbiota-mediated processes may confer advantages to such insect guilds (Fäldt, 2000 and references therein), a phenomenon that may well apply to H. abietis and thus needs to be examined further.
Conclusion
Our study supports a large body of evidence accumulated over decades that host chemical responses are involved in resistance to herbivory at the intraspecific level. However, our study provides the first example of a comprehensive approach in which PSM inducibility in trees, and its role in resistance, are analyzed in a rigorous genetic framework that allows the exclusion of spurious associations due to among-subject relatedness. It follows that multivariate analyses of PSM, conducted in this way, provide more realistic information about the potentially causal relationships between PSM and resistance to herbivory than evaluating single pairwise relationships or using broad groups of compounds.
Data Accessibility
The SNP dataset analyzed in this study is publicly available in the Zenodo repository (doi: 10.5281/zenodo.1445313).
Author Contributions
LS and RZ designed the experiment, provided the reagents, performed the sampling, helped with the statistical and GC analyses and the interpretation of the results and improved the different versions of the manuscript. PB provided the infrastructure and know-how for the UHPLC analyses of phenolic compounds and helped with the interpretation of the associated data. CV helped with the UHPLC analyses of phenolic compounds and with the interpretation of the associated data. AB-K provided the infrastructure and background for GC-MS analyses of terpenoid compounds and helped with the interpretation of the associated data. DG developed the SNPs database, constructed the population structure and kinship matrices, and provided expertise on controlling for genetic relatedness. XL-G performed the sampling, all the chemical analyses, most of the statistical analyses, produced the results, wrote the first draft with RZ and LS, and lead on the improvement of subsequent versions of the manuscript. PB, CV, AB-K, and DG contributed equally to the final version of the manuscript.
Funding
This research was supported by the Spanish Government via the MINECO/FEDER grants FENOPIN (AGL2012-40151-C03-01), FUTURPIN (AGL2015-68274-C03-02-R) and AdapCon (CGL2011-30182-C02-01/02), by the European Union’s Horizon 2020 Programme for Research & Innovation (grant id 773383), and by the Xunta de Galicia – GAIN grant IN607A-2016/13. XL-G received financial support from Fundación Pedro Barrié de la Maza, and from the FPI and Short-Term Stay Grant programs (MINECO). We acknowledge support of the publication fee by the CSIC Open Access Publication Support Initiative through its Unit of Information Resources for Research (URICI).
Conflict of Interest Statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Acknowledgments
We thank Dr. Juan Majada (SERIDA, Asturias) and Ricardo Alía (CIFOR-INIA, Madrid) and their corresponding research groups for conceiving, growing and providing the plant material; José Carlos Rodríguez and Santiago Rodal (Viveros Figueirido, NORFOR, Pontevedra, Spain) for plant maintenance and culture, César Cendán and the field personnel of Misión Biológica de Galicia (MBG-CSIC) for their help during plant sampling, and María Luz Pato and Helena Pazó for the help with the plant sampling, measurements and chemical analyses. Also thanks to Jean-Christophe Cocuron for the support provided with the UHPLC-DAD-MS analyses, and Santiago C. González-Martínez, Carmen García Barriga, and Diana Barba Egido for the help with SNP selection, DNA extraction, and SNP editing, respectively.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpls.2018.01651/full#supplementary-material
References
Abdala-Roberts, L., Rasmann, S., Berny-Mier, Y. T. J. C., Covelo, F., Glauser, G., and Moreira, X. (2016). Biotic and abiotic factors associated with altitudinal variation in plant traits and herbivory in a dominant oak species. Am. J. Bot. 103, 2070–2078. doi: 10.3732/ajb.1600310
Agrawal, A. A. (2005). Natural selection on common milkweed by a community of specialized insect herbivores. Evol. Ecol. Res. 7, 651–667.
Agrawal, A. A. (2010). Current trends in the evolutionary ecology of plant defence. Funct. Ecol. 25, 420–432. doi: 10.1111/j.1365-2435.2010.01796.x
Agrawal, A. A., Conner, J. K., Johnson, M. T. J., and Wallsgrove, R. (2002). Ecological genetics of an induced plant defense against herbivores: additive genetic variance and costs of phenotypic plasticity. Evolution 56, 2206–2213. doi: 10.1111/j.0014-3820.2002.tb00145.x
Agrawal, A. A., Hastings, A. P., Bradburd, G. S., Woods, E. C., Zust, T., Harvey, J. A., et al. (2015). Evolution of plant growth and defense in a continental introduction. Am. Nat. 186, E1–E15. doi: 10.1086/681622
Benjamini, Y., and Hochberg, Y. (1995). Controlling the False Discovery Rate: a practical and powerful approach to multiple testing. J. R. Stat. Soc. Ser. B 57, 289–300.
Berasategui, A., Salem, H., Paetz, C., Santoro, M., Gershenzon, J., Kaltenpoth, M., et al. (2017). Gut microbiota of the pine weevil degrades conifer diterpenes and increases insect fitness. Mol. Ecol. 26, 4099–4110. doi: 10.1111/mec.14186
Bernard-Dagan, C. (1988). “Seasonal variations in energy sources and biosynthesis,” in Mechanisms of Woody Plant Defenses Against Insects: Search for Pattern, eds W. J. Mattson, J. Levieux, and C. Bernard-Dagan (New York, NY: Springer-Verlag).
Björklund, N. (2007). Cues for shelter use in a phytophagous insect. J. Insect Behav. 21, 9–23. doi: 10.1007/s10905-007-9100-5
Boone, C. K., Keefover-Ring, K., Mapes, A. C., Adams, A. S., Bohlmann, J., and Raffa, K. F. (2013). Bacteria associated with a tree-killing insect reduce concentrations of plant defense compounds. J. Chem. Ecol. 39, 1003–1006. doi: 10.1007/s10886-013-0313-0
Borg-Karlson, A. K., Nordlander, G., Mudalige, A., Nordenhem, H., and Unelius, C. R. (2006). Antifeedants in the feces of the pine weevil Hylobius abietis: identification and biological activity. J. Chem. Ecol. 32, 943–957. doi: 10.1007/s10886-006-9050-y
Bucci, G., Gonzalez-Martinez, S. C., Le Provost, G., Plomion, C., Ribeiro, M. M., Sebastiani, F., et al. (2007). Range-wide phylogeography and gene zones in Pinus pinaster Ait. revealed by chloroplast microsatellite markers. Mol. Ecol. 16, 2137–2153. doi: 10.1111/j.1365-294X.2007.03275.x
Cárdenas, R. E., Valencia, R., Kraft, N. J. B., Argoti, A., Dangles, O., and Swenson, N. (2014). Plant traits predict inter- and intraspecific variation in susceptibility to herbivory in a hyperdiverse Neotropical rain forest tree community. J. Ecol. 102, 939–952. doi: 10.1111/1365-2745.12255
Carmona, D., Lajeunesse, M. J., and Johnson, M. T. J. (2011). Plant traits that predict resistance to herbivores. Funct. Ecol. 25, 358–367. doi: 10.1111/j.1365-2435.2010.01794.x
Carrillo-Gavilán, A., Moreira, X., Zas, R., Gonzalez-Voyer, A., Vilà, M., Sampedro, L., et al. (2015). Phylogenetic and biogeographical patterns in defensive strategies and quantitative allocation to chemical defences in Palaearctic and Nearctic pine trees. J. Biogeogr. 42, 684–693. doi: 10.1111/jbi.12444
Cipollini, D., and Heil, M. (2010). Costs and benefits of induced resistance to herbivores and pathogens in plants. CAB Rev. Perspect. Agric. Vet. Sci. Nutr. Nat. Resour. 5, 1–25. doi: 10.1111/j.1469-8137.2009.03066.x
Conord, C., Lemperiere, G., Taberlet, P., and Despres, L. (2006). Genetic structure of the forest pest Hylobius abietis on conifer plantations at different spatial scales in Europe. Heredity 97, 46–55. doi: 10.1038/sj.hdy.6800837
Danielsson, M., Kännaste, A., Lindström, A., Hellqvist, C., Stattin, E., Långström, B., et al. (2008). Mini-seedlings of Picea abies are less attacked by Hylobius abietis than conventional ones: is plant chemistry the explanation? Scand. J. For. Res. 23, 299–306. doi: 10.1080/02827580802203560
Day, K. R., Nordlander, G., Kenis, M., and Halldorson, G. (2004). “General biology and life cycles of bark weevils,” in Bark and Wood Boring Insects in Living Trees in Europe, a Synthesis, eds F. Lieutier, K. R. Day, A. Battisti, J. C. Grégoire, and H. F. Evans (Dordrecht: Kluwer Academic Publishers), 331–349. doi: 10.1007/978-1-4020-2241-8_14
De-Lucas, A. I., González-Martínez, S. C., Vendramin, G. G., Hidalgo, E., and Heuertz, M. (2009). Spatial genetic structure in continuous and fragmented populations of Pinus pinaster Aiton. Mol. Ecol. 18, 4564–4576. doi: 10.1111/j.1365-294X.2009.04372.x
Dudareva, N., Klempien, A., Muhlemann, J. K., and Kaplan, I. (2013). Biosynthesis, function and metabolic engineering of plant volatile organic compounds. New Phytol. 198, 16–32. doi: 10.1111/nph.12145
Elvira-Recuenco, M., Iturritxa, E., Majada, J., Alía, R., and Raposo, R. (2014). Adaptive potential of maritime pine (Pinus pinaster) populations to the emerging pitch canker pathogen, Fusarium circinatum. PLoS One 9:e114971. doi: 10.1371/journal.pone.0114971
Erbilgin, N., Krokene, P., Christiansen, E., Zeneli, G., and Gershenzon, J. (2006). Exogenous application of methyl jasmonate elicits defenses in Norway spruce (Picea abies) and reduces host colonization by the bark beetle Ips typographus. Oecologia 148, 426–436. doi: 10.1007/s00442-006-0394-3
Ericsson, A., Gref, R., Hellqvist, C., and Lagstrom, B. (1988). “Wound response of living bark of scots pine seedlings and its influence on feeding by the weevil, Hylobius abietis,” in Mechanisms of Woody plant Defenses Against Insects, ed. W. J. Mattson (New York, NY: Springer-Verlag),227–235.
Eriksson, C., Mansson, P. E., Sjodin, K., and Schlyter, F. (2008). Antifeedants and feeding stimulants in bark extracts of ten woody non-host species of the pine weevil, Hylobius abietis. J. Chem. Ecol. 34, 1290–1297. doi: 10.1007/s10886-008-9525-0
EUFORGEN (2009). Distribution Map of Maritime pine (Pinus pinaster) [Online]. Available at: www.euforgen.org [accessed June 14, 2018].
Fäldt, J. (2000). Volatile Constituents in Conifers and Conifer-Related Wood-Decaying Fungi. Ph.D. thesis, Stockholm, KTH Royal Institute of Technology.
Feeny, P. (1976). “Plant apparency and chemical defense,” in Biochemical Interaction between Plants and Insects, 1st Edn, eds J. W. Wallace and R. L. Mansell (New York, NY: Plenum Press), 1–40.
Freckleton, R. P., Harvey, P. H., and Pagel, M. (2002). Phylogenetic analysis and comparative data: a test and review of evidence. Am. Nat. 160, 712–726. doi: 10.1086/343873
Futuyma, D. J., and Agrawal, A. A. (2009). Macroevolution and the biological diversity of plants and herbivores. Proc. Natl. Acad. Sci. U.S.A. 106, 18054–18061. doi: 10.1073/pnas.0904106106
Gershenzon, J. (1994). Metabolic costs of terpenoid accumulation in higher plants. J. Chem. Ecol. 20, 1281–1328. doi: 10.1007/BF02059810
González-Martínez, S. C., Burczyk, J., Nathan, R., Nanos, N., Gil, L., and Alía, R. (2006). Effective gene dispersal and female reproductive success in Mediterranean maritime pine (Pinus pinaster Aiton). Mol. Ecol. 15, 4577–4588. doi: 10.1111/j.1365-294X.2006.03118.x
Goodger, J. Q., Heskes, A. M., and Woodrow, I. E. (2013). Contrasting ontogenetic trajectories for phenolic and terpenoid defences in Eucalyptus froggattii. Ann. Bot. 112, 651–659. doi: 10.1093/aob/mct010
Grivet, D., Climent, J., Zabal-Aguirre, M., Neale, D. B., Vendramin, G. G., and Gonzalez-Martinez, S. C. (2013). Adaptive evolution of Mediterranean pines. Mol. Phylogenet. Evol. 68, 555–566. doi: 10.1016/j.ympev.2013.03.032
Hahn, P. G., and Maron, J. L. (2016). A framework for predicting intraspecific variation in plant defense. Trends Ecol. Evol. 31, 646–656. doi: 10.1016/j.tree.2016.05.007
Hamilton, J. G., Zangerl, A. R., DeLucia, E. H., and Berenbaum, M. R. (2001). The carbon-nutrient balance hypothesis: its rise and fall. Ecol. Lett. 4, 86–95. doi: 10.1046/j.1461-0248.2001.00192.x
Hardy, O., and Vekemans, X. (2002). SPAGeDi: a versatile computer program to analyse spatial genetic structure at the individual or population levels. Mol. Ecol. Notes 2, 618–620. doi: 10.1046/j.1471-8286.2002.00305.x
Heijari, J., Nerg, A., Kainulainen, P., Viiri, H., Vuorinen, M., and Holopainen, J. K. (2005). Application of methyl jasmonate reduces growth but increases chemical defence and resistance againt Hylobius abietis in Scots pine seedlings. Entomol. Exp. Appl. 115, 117–124. doi: 10.1111/j.1570-7458.2005.00263.x
Heil, M. (2009). Plastic defence expression in plants. Evol. Ecol. 24, 555–569. doi: 10.1007/s10682-009-9348-7
Hudgins, J. W., Christiansen, E., and Franceschi, V. R. (2004). Induction of anatomically based defense responses in stems of diverse conifers by methyl jasmonate: a phylogenetic perspective. Tree Physiol. 24, 251–264. doi: 10.1093/treephys/24.3.251
Iason, G. R., O’Reilly-Wapstra, J. M., Brewer, M. J., Summers, R. W., and Moore, B. D. (2011). Do multiple herbivores maintain chemical diversity of Scots pine monoterpenes? Philos. Trans. R. Soc. Lond. B Biol. Sci. 366, 1337–1345. doi: 10.1098/rstb.2010.0236
Ishangulyyeva, G., Najar, A., Curtis, J. M., and Erbilgin, N. (2016). Fatty acid composition of novel host jack pine do not prevent host acceptance and colonization by the invasive mountain pine beetle and its symbiotic fungus. PLoS One 11:e0162046. doi: 10.1371/journal.pone.0162046
Iturritxa, E., Ganley, R. J., Raposo, R., García-Serna, I., Mesanza, N., Kirkpatrick, S. C., et al. (2013). Resistance levels of Spanish conifers against Fusarium circinatum and Diplodia pinea. For. Pathol. 43, 488–495. doi: 10.1111/efp.12061
Jaramillo-Correa, J. P., Rodriguez-Quilon, I., Grivet, D., Lepoittevin, C., Sebastiani, F., Heuertz, M., et al. (2015). Molecular proxies for climate maladaptation in a long-lived tree (Pinus pinaster Aiton, Pinaceae). Genetics 199, 793–807. doi: 10.1534/genetics.114.173252
Kannaste, A., Nordenhem, H., Nordlander, G., and Borg-Karlson, A. K. (2009). Volatiles from a mite-infested spruce clone and their effects on pine weevil behavior. J. Chem. Ecol. 35, 1262–1271. doi: 10.1007/s10886-009-9708-3
Karban, R. (2011). The ecology and evolution of induced resistance against herbivores. Funct. Ecol. 25, 339–347. doi: 10.1111/j.1365-2435.2010.01789.x
Karban, R., Agrawal, A. A., Thaler, J. S., and Adler, L. S. (1999). Induced plant responses and information content about risk of herbivory. Trends Ecol. Evol. 14, 443–447. doi: 10.1016/S0169-5347(99)01678-X
Kremer, A., Ronce, O., Robledo-Arnuncio, J. J., Guillaume, F., Bohrer, G., Nathan, R., et al. (2012). Long-distance gene flow and adaptation of forest trees to rapid climate change. Ecol. Lett. 15, 378–392. doi: 10.1111/j.1461-0248.2012.01746.x
Langström, B., and Day, K. R. (2004). “Damage, control and management of weevil pests, especially Hylobius abietis,” in Bark and Wood Boring Insects in Living Trees in Europe, A Synthesis, eds F. Lieutier, K. R. Day, A. Battisti, J. C. Grégoire, and H. F. Evans (Berlin: Springer), 415–444.
Larsson, S., Björkman, C., and Gref, R. (1986). Responses of Neodiprion sertifer (Hym., Diprionidae) larvae to variation in needle resin acid concentration in Scots pine. Oecologia 70, 77–84. doi: 10.1007/BF00377113
Leather, S. R., Day, K. R., and Salisbury, A. N. (1999). The biology and ecology of the large pine weevil, Hylobius abietis (Coleoptera: Curculionidae): a problem of dispersal? Bull. Entomol. Res. 89, 3–16. doi: 10.1017/s0007485399000024
Littell, R. C., Milliken, G. A., Stroup, W. W., Wolfinger, R., and Schabenberger, O. (2006). SAS System for Mixed Models. Cary, NC: SAS Institute Inc.
Loiselle, B., Sork, V., Nason, J., and Graham, C. (1995). Spatial genetic structure of a tropical understory shrub, Psychotria officinalis (Rubiaceae). Am. J. Bot. 82, 1420–1425. doi: 10.1002/j.1537-2197.1995.tb12679.x
Lombardero, M. J., Ayres, M. P., Lorio, P. L., and Ruel, J. (2000). Environmental effects on constitutive and inducible resin defences of Pinus taeda. Ecol. Lett. 3, 329–339. doi: 10.1046/j.1461-0248.2000.00163.x
López-Goldar, X., Sampedro, L., and Zas, R. (2016). Carbon starvation by light deprivation does not constrain the ability of young pines to produce induced chemical defences in response to a bark-chewing herbivore. Environ. Exp. Bot. 130, 141–150. doi: 10.1016/j.envexpbot.2016.05.010
Lundborg, L., Fedderwitz, F., Bjorklund, N., Nordlander, G., and Borg-Karlson, A. K. (2016). Induced defenses change the chemical composition of pine seedlings and influence meal properties of the pine weevil Hylobius abietis. Phytochemistry 130, 99–105. doi: 10.1016/j.phytochem.2016.06.002
Majada, J., Martínez-Alonso, C., Feito, I., Kidelman, A., Aranda, I., and Alía, R. (2011). Minicuttings: an effective technique for the propagation of Pinus pinaster Ait. New For. 41, 399–412. doi: 10.1007/s11056-010-9232-x
Mason, C. M., Bowsher, A. W., Crowell, B. L., Celoy, R. M., Tsai, C.-J., and Donovan, L. A. (2016). Macroevolution of leaf defenses and secondary metabolites across the genus Helianthus. New Phytol. 209, 1720–1733. doi: 10.1111/nph.13749
Moles, A. T., Wallis, I. R., Foley, W. J., Warton, D. I., Stegen, J. C., Bisigato, A. J., et al. (2011). Putting plant resistance traits on the map: a test of the idea that plants are better defended at lower latitudes. New Phytol. 191, 777–788. doi: 10.1111/j.1469-8137.2011.03732.x
Moore, B. D., Andrew, R. L., Kulheim, C., and Foley, W. J. (2013). Explaining intraspecific diversity in plant secondary metabolites in an ecological context. New Phytol. 201, 733–750. doi: 10.1111/nph.12526
Moreira, X., Castagneyrol, B., Abdala-Roberts, L., Berny-Mier y Teran, J. C., Timmermans, B. G. H., Bruun, H. H., et al. (2018). Latitudinal variation in plant chemical defences drives latitudinal patterns of leaf herbivory. Ecography 41, 1124–1134. doi: 10.1111/ecog.03326
Moreira, X., Lundborg, L., Zas, R., Carrillo-Gavilan, A., Borg-Karlson, A. K., and Sampedro, L. (2013a). Inducibility of chemical defences by two chewing insect herbivores in pine trees is specific to targeted plant tissue, particular herbivore and defensive trait. Phytochemistry 94, 113–122. doi: 10.1016/j.phytochem.2013.05.008
Moreira, X., Zas, R., and Sampedro, L. (2013b). Additive genetic variation in resistance traits of an exotic pine species: little evidence for constraints on evolution of resistance against native herbivores. Heredity 110, 449–456. doi: 10.1038/hdy.2012.108
Moreira, X., Mooney, K. A., Rasmann, S., Petry, W. K., Carrillo-Gavilán, A., Zas, R., et al. (2014). Trade-offs between constitutive and induced defences drive geographical and climatic clines in pine chemical defences. Ecol. Lett. 17, 537–546. doi: 10.1111/ele.12253
Moreira, X., Sampedro, L., and Zas, R. (2009). Defensive responses of Pinus pinaster seedlings to exogenous application of methyl jasmonate: concentration effect and systemic response. Environ. Exp. Bot. 67, 94–100. doi: 10.1016/j.envexpbot.2009.05.015
Moreira, X., Zas, R., and Sampedro, L. (2012). Quantitative comparison of chemical, biological and mechanical induction of secondary compounds in Pinus pinaster seedlings. Trees 26, 677–683. doi: 10.1007/s00468-011-0602-6
Moreira, X., Zas, R., Solla, A., and Sampedro, L. (2015). Differentiation of persistent anatomical defensive structures is costly and determined by nutrient availability and genetic growth-defence constraints. Tree Physiol. 35, 112–123. doi: 10.1093/treephys/tpu106
Morris, W. F., Traw, M. B., and Bergelson, J. (2006). On testing for a tradeoff between constitutive and induced resistance. Oikos 112, 102–110. doi: 10.1111/j.0030-1299.2006.14253.x
Nordenhem, H., and Nordlander, G. (1994). Olfactory oriented migration through soil by root-living Hylobius abietis (L.) larvae (Col., Curculionidae). J. Appl. Entomol. 117, 457–462. doi: 10.1111/j.1439-0418.1994.tb00762.x
Nordlander, G. (1987). A method for trapping Hylobius abietis(L.) with a standardized bait and its potential for forecasting seedling damage. Scand. J. For. Res. 2, 199–213. doi: 10.1080/02827588709382458
Nordlander, G., Bylund, H., Örlander, G., and Wallertz, K. (2003). Pine weevil population density and damage to coniferous seedlings in a regeneration area with and without shelterwood. Scand. J. For. Res. 18, 438–448. doi: 10.1080/02827580310001634
Nordlander, G., Hellqvist, C., Johansson, K., and Nordenhem, H. (2011). Regeneration of European boreal forests: effectiveness of measures against seedling mortality caused by the pine weevil Hylobius abietis. For. Ecol. Manag. 262, 2354–2363. doi: 10.1016/j.foreco.2011.08.033
Örlander, G., Nilsson, U., and Nordlander, G. (1997). Pine weevil abundance on clear-cuttings of different ages: a 6-year study using pitfall traps. Scand. J. For. Res. 12, 225–240. doi: 10.1080/02827589709355405
Pagel, M. (1999). Inferring the historical patterns of biological evolution. Nature 401, 877–884. doi: 10.1038/44766
Pellissier, L., Fiedler, K., Ndribe, C., Dubuis, A., Pradervand, J. N., Guisan, A., et al. (2012). Shifts in species richness, herbivore specialization, and plant resistance along elevation gradients. Ecol. Evol. 2, 1818–1825. doi: 10.1002/ece3.296
Phillips, M. A., and Croteau, R. B. (1999). Resin-based defenses in conifers. Trends Plant Sci. 4, 184–190. doi: 10.1016/S1360-1385(99)01401-6
Porth, I., Hamberger, B., White, R., and Ritland, K. (2011). Defense mechanisms against herbivory in Picea: sequence evolution and expression regulation of gene family members in the phenylpropanoid pathway. BMC Genomics 12:608. doi: 10.1186/1471-2164-12-608
Powell, J. S., and Raffa, K. E. (1999). Effects of selected Larix laricina terpenoids on Lymantria dispar development and behavior. Environ. Entomol. 28, 148–154. doi: 10.1093/ee/28.2.148
Pritchard, J., Stephens, M., and Donnelly, P. (2000). Inference of population structure using multilocus genotype data. Genetics 155, 945–959.
Ralph, S. G., Yueh, H., Friedmann, M., Aeschliman, D., Zeznik, J. A., Nelson, C. C., et al. (2006). Conifer defence against insects: microarray gene expression profiling of Sitka spruce (Picea sitchensis) induced by mechanical wounding or feeding by spruce budworms (Choristoneura occidentalis) or white pine weevils (Pissodes strobi) reveals large-scale changes of the host transcriptome. Plant Cell Environ. 29, 1545–1570. doi: 10.1111/j.1365-3040.2006.01532.x
Rasmann, S., and Agrawal, A. A. (2011). Latitudinal patterns in plant defense: evolution of cardenolides, their toxicity and induction following herbivory. Ecol. Lett. 14, 476–483. doi: 10.1111/j.1461-0248.2011.01609.x
Rasmann, S., Buri, A., Gallot-Lavallée, M., Joaquim, J., Purcell, J., Pellissier, L., et al. (2014). Differential allocation and deployment of direct and indirect defences by Vicia sepium along elevation gradients. J. Ecol. 102, 930–938. doi: 10.1111/1365-2745.12253
Sampedro, L., Moreira, X., Llusia, J., Penuelas, J., and Zas, R. (2010). Genetics, phosphorus availability, and herbivore-derived induction as sources of phenotypic variation of leaf volatile terpenes in a pine species. J. Exp. Bot. 61, 4437–4447. doi: 10.1093/jxb/erq246
Sampedro, L., Moreira, X., and Zas, R. (2011a). Costs of constitutive and herbivore-induced chemical defences in pine trees emerge only under low nutrient availability. J. Ecol. 99, 818–827. doi: 10.1111/j.1365-2745.2011.01814.x
Sampedro, L., Moreira, X., and Zas, R. (2011b). Resistance and response of Pinus pinaster seedlings to Hylobius abietis after induction with methyl jasmonate. Plant Ecol. 212, 397–401. doi: 10.1007/s11258-010-9830-x
Savolainen, O., Pyhäjärvi, T., and Knürr, T. (2007). Gene flow and local adaptation in trees. Annu. Rev. Ecol. Evol. Syst. 38, 595–619. doi: 10.1146/annurev.ecolsys.38.091206.095646
Schiebe, C., Hammerbacher, A., Birgersson, G., Witzell, J., Brodelius, P. E., Gershenzon, J., et al. (2012). Inducibility of chemical defenses in Norway spruce bark is correlated with unsuccessful mass attacks by the spruce bark beetle. Oecologia 170, 183–198. doi: 10.1007/s00442-012-2298-8
Serra-Varela, M. J., Grivet, D., Vincenot, L., Broennimann, O., Gonzalo-Jiménez, J., and Zimmermann, N. E. (2015). Does phylogeographical structure relate to climatic niche divergence? A test using maritime pine (Pinus pinaster Ait.). Glob. Ecol. Biogeogr. 24, 1302–1313. doi: 10.1111/geb.12369
Solbreck, C., and Gyldberg, B. (1979). Temporal flight pattern of the large pine weevil Hylobius abietis L. (Coleoptera, Curculionidae), with special reference to the influence of weather. J. Appl. Entomol. 88, 532–536. doi: 10.1111/j.1439-0418.1979.tb02532.x
Tan, J. Y., Wainhouse, D., Morgan, G., and Day, K. R. (2011). Interaction between flight, reproductive development and oviposition in the pine weevil Hylobius abietis. Agric. For. Entomol. 13, 149–156. doi: 10.1111/j.1461-9563.2010.00503.x
Tapias, R., Climent, J., Pardos, J. A., and Gil, L. (2004). Life histories of Mediterranean pines. Plant Ecol. 171, 53–68. doi: 10.1023/B:VEGE.0000029383.72609.f0
Theis, N., and Lerdau, M. (2003). The evolution of function in plant secondary metabolites. Int. J. Plant Sci. 164(3 Suppl.), S93–S102. doi: 10.1086/374190
Tomlin, E. S., Borden, J. H., and Pierce, H. D. Jr. (1996). Relationship between cortical resin acids and resistance of Sitka spruce to the white pine weevil. Can. J. Bot. 74, 599–606. doi: 10.1139/b96-076
Treutter, D. (2006). Significance of flavonoids in plant resistance: a review. Environ. Chem. Lett. 4, 147–157. doi: 10.1007/s10311-006-0068-8
Vendramin, G. G., Anzidei, M., Madaghiele, A., and Bucci, G. (1998). Distribution of genetic diversity in Pinus pinaster Ait. as revealed by chloroplast microsatellites. Theor. Appl. Genet. 97, 456–463. doi: 10.1007/s001220050917
Villari, C., Herms, D. A., Whitehill, J. G., Cipollini, D., and Bonello, P. (2016). Progress and gaps in understanding mechanisms of ash tree resistance to emerald ash borer, a model for wood-boring insects that kill angiosperms. New Phytol. 209, 63–79. doi: 10.1111/nph.13604
Wagner, M. R., Benjamin, D. M., Clancy, K. M., and Schuh, B. A. (1983). Influence of diterpene resin acids on feeding and growth of larch sawfly, Pristiphora erichsonii (Hartig). J. Chem. Ecol. 9, 119–127. doi: 10.1007/BF00987776
Wainhouse, D., Boswell, R., and Ashburner, R. (2004). Maturation feeding and reproductive development in adult pine weevil, Hylobius abietis (Coleoptera: Curculionidae). Bull. Entomol. Res. 94, 81–87. doi: 10.1079/ber2003283
Wainhouse, D., Inward, D. J. G., and Morgan, G. (2014). Modelling geographical variation in voltinism of Hylobius abietis under climate change and implications for management. Agric. For. Entomol. 16, 136–146. doi: 10.1111/afe.12043
Wainhouse, D., Staley, J., Johnston, J., and Boswell, R. (2005). The effect of environmentally induced changes in the bark of young conifers on feeding behaviour and reproductive development of adult Hylobius abietis (Coleoptera: Curculionidae). Bull. Entomol. Res. 95, 151–159. doi: 10.1079/ber2004344
Wallis, C. M., Reich, R. W., Lewis, K. J., and Huber, D. P. W. (2010). Lodgepole pine provenances differ in chemical defense capacities against foliage and stem diseases. Can. J. For. Res. 40, 2333–2344. doi: 10.1139/x10-178
Wiggins, N. L., Forrister, D. L., Endara, M. J., Coley, P. D., and Kursar, T. A. (2016). Quantitative and qualitative shifts in defensive metabolites define chemical defense investment during leaf development in Inga, a genus of tropical trees. Ecol. Evol. 6, 478–492. doi: 10.1002/ece3.1896
Wingfield, M. J., Hammerbacher, A., Ganley, R., Steenkamp, E. T., Gordon, T. R., Wingfield, B. D., et al. (2008). Pitch canker caused by Fusarium circinatum - a growing threat to pine plantations and forests worldwide. Aust. Plant Pathol. 37, 319–334. doi: 10.1071/AP08036
Witzell, J., and Martín, J. A. (2008). Phenolic metabolites in the resistance of northern forest trees to pathogens — past experiences and future prospects. Can. J. For. Res. 38, 2711–2727. doi: 10.1139/x08-112
Yu, J., Pressoir, G., Briggs, W. H., Vroh Bi, I., Yamasaki, M., Doebley, J. F., et al. (2006). A unified mixed-model method for association mapping that accounts for multiple levels of relatedness. Nat. Genet. 38, 203–208. doi: 10.1038/ng1702
Zas, R., Bjorklund, N., Nordlander, G., Cendan, C., Hellqvist, C., and Sampedro, L. (2014). Exploiting jasmonate-induced responses for field protection of conifer seedlings against a major forest pest. Hylobius abietis. For. Ecol. Manage. 313, 212–223. doi: 10.1016/j.foreco.2013.11.014
Zas, R., Björklund, N., Sampedro, L., Hellqvist, C., Karlsson, B., Jansson, S., et al. (2017). Genetic variation in resistance of Norway spruce seedlings to damage by the pine weevil Hylobius abietis. Tree Genet. Genomes 13:111. doi: 10.1007/s11295-017-1193-1
Zas, R., Moreira, X., Ramos, M., Lima, M. R. M., Nunes da Silva, M., Solla, A., et al. (2015). Intraspecific variation of anatomical and chemical defensive traits in Maritime pine (Pinus pinaster) as factors in susceptibility to the pinewood nematode (Bursaphelenchus xylophilus). Trees 29, 663–673. doi: 10.1007/s00468-014-1143-6
Zas, R., Sampedro, L., Prada, E., and Fernandez-Lopez, J. (2005). Genetic variation of Pinus pinaster Ait. seedlings in susceptibility to the pine weevil Hylobius abietis L. Ann. For. Sci. 62, 681–688. doi: 10.1051/forest:2005064
Zas, R., Sampedro, L., Prada, E., Lombardero, M. J., and Fernandez-Lopez, J. (2006). Fertilization increases Hylobius abietis L. damage in Pinus pinaster Ait. seedlings. For. Ecol. Manag. 222, 137–144. doi: 10.1016/j.foreco.2005.10.008
Zhao, T., Borg-Karlson, A. K., Erbilgin, N., and Krokene, P. (2011a). Host resistance elicited by methyl jasmonate reduces emission of aggregation pheromones by the spruce bark beetle, Ips typographus. Oecologia 167, 691–699. doi: 10.1007/s00442-011-2017-x
Zhao, T., Krokene, P., Hu, J., Christiansen, E., Björklund, N., Långström, B., et al. (2011b). Induced terpene accumulation in Norway spruce inhibits bark beetle colonization in a dose-dependent manner. PLoS One 6:e26649. doi: 10.1371/journal.pone.0026649
Zhao, T., Krokene, P., Bjorklund, N., Langstrom, B., Solheim, H., Christiansen, E., et al. (2010). The influence of Ceratocystis polonica inoculation and methyl jasmonate application on terpene chemistry of Norway spruce, Picea abies. Phytochemistry 71, 1332–1341. doi: 10.1016/j.phytochem.2010.05.017
Keywords: genetic variation, herbivory, inducibility, maritime pine, plant secondary metabolites (PSM), phenolics, resistance, terpenes
Citation: López-Goldar X, Villari C, Bonello P, Borg-Karlson AK, Grivet D, Zas R and Sampedro L (2018) Inducibility of Plant Secondary Metabolites in the Stem Predicts Genetic Variation in Resistance Against a Key Insect Herbivore in Maritime Pine. Front. Plant Sci. 9:1651. doi: 10.3389/fpls.2018.01651
Received: 07 August 2018; Accepted: 24 October 2018;
Published: 21 November 2018.
Edited by:
Judy Simon, Universität Konstanz, GermanyReviewed by:
Axel Schmidt, Max-Planck-Institut für Chemische Ökologie, GermanyStefanie Wienkoop, Universität Wien, Austria
Copyright © 2018 López-Goldar, Villari, Bonello, Borg-Karlson, Grivet, Zas and Sampedro. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Xosé López-Goldar, xlgoldar@mbg.csic.es