Carla Vázquez-González
Xosé López-Goldar
Rafael Zas
Luis Sampedro- Misión Biológica de Galicia, Consejo Superior de Investigaciones Científicas (CSIC), Pontevedra, Spain
Resin ducts are important anatomical defensive traits related to biotic resistance in conifers. Previous studies have reported intraspecific genetic variation in resin duct characteristics. However, little is currently known about the micro-evolutionary patterns and adaptive value of these defensive structures. Here, we quantified inter-population genetic variation in resin duct features and their inducibility in Pinus pinaster and assessed whether such variation was associated with climate gradients. To that end, we characterized the resin duct system of 2-year-old saplings from 10 populations across the species’ distribution range. We measured axial resin duct features (density, mean size, and percentage conductive area of resin ducts) and their inducibility in response to methyl jasmonate. Genotyping of single nucleotide polymorphisms allowed to account for the population genetic structure in our models in order to avoid spurious correlations between resin duct characteristics and climate. We found large inter-population variation in resin duct density and conductive area, but not in their inducibility. Our results suggest that population variation in the percentage conductive area of resin ducts likely arise from adaptation to local climate conditions. This study highlights the adaptive relevance of resin ducts and helps to shed light on the micro-evolutionary patterns of resin-based defenses in conifers.
Introduction
Conifers have developed a broad array of physical and chemical defensive traits and strategies to reduce the damage caused by a large and diverse community of antagonist organisms (Franceschi et al., 2005; Mason et al., 2019; Whitehill et al., 2019). Among them, resin ducts are the structures that produce and exude the oleoresin (Fahn, 1979), an effective terpenoid-based chemical defense against insects and pathogens (Keeling and Bohlmann, 2006). Resin duct characteristics (e.g. area or abundance) are considered important proxies of defensive investment in conifer trees because they have been positively associated with enhanced resistance against insects (Kane and Kolb, 2010; Ferrenberg et al., 2014; Hood et al., 2015) and pathogens (Christiansen et al., 1999; Krokene, et al., 2003). Larger resin ducts may enhance resistance because small increases in resin duct lumen area can considerably increase resin flow (Schopmeyer et al., 1954). Accordingly, several studies have reported a positive correlation between resin duct investment and resin flow (Lombardero et al., 2000; Rodriguez-Garcia et al., 2014; Hood and Sala, 2015). Besides their role in oleoresin accumulation, resin ducts act as an anatomical barrier, and their spatial distribution may be crucial in determining the success of an attack. Supporting this idea, Boucher et al. (2001) reported that a denser network of smaller resin ducts was positively correlated with resistance to Pissodes strobi Peck because it leaves fewer undefended gaps for the entrance of organisms. In contrast, bigger and fewer resin ducts would lead to larger gaps, favoring insect accessibility to feeding tissues (Reed et al., 1986).
Resin ducts are present constitutively (i.e., in the absence of any biotic stress), and are also induced in response to damage signaling. Several external stimuli have been reported to induce the formation of "traumatic" resin ducts, including insect attacks, pathogen infestation, or wounding alone (Franceschi et al., 2000; Nagy et al., 2000; Chen et al., 2018). Resin ducts can also be induced by the application of methyl jasmonate (MJ) (Franceschi et al., 2002; Martin et al., 2002; Krokene et al., 2008; Moreira et al., 2015b), a hormone implicated in damage signaling and defensive responses (Hudgins and Franceschi, 2004). Several studies have reported increased resistance against insects and pathogens by means of resin duct induction with different elicitors. For instance, pre-inoculations of Endoconidiophora polonica (previously Ceratocystis polonica (Siemaszko) C. Moreau) or MJ application in Picea abies (L.) Karst trees enhanced traumatic resin duct production and biotic resistance to bark beetles and their associated pathogens in later infections (Christiansen et al., 1999; Franceschi et al., 2002; Krokene et al., 2003; Erbilgin et al., 2006; Krokene et al., 2008). The development of traumatic or induced resin ducts implies, however, a cellular differentiation process at the cambium that can take several weeks (Nagy et al., 2000). Thus, induced resin ducts may not fully develop in time to avoid the effects of an initial attack.
Research on tree defenses has become a priority given the rising biotic stress caused by the progressive expansion of native and newly arrived alien antagonist organisms (Raffa et al., 2008; Santini et al., 2013; Early et al., 2016). Moreover, abiotic stress driven by climate change can also affect biotic resistance in tree species, because more limiting conditions can affect resource allocation to defenses (Herms and Mattson, 1992; Sampedro et al., 2011; Moreira et al., 2015b). Phenotypic plasticity is therefore an important mechanism for dealing with increasing environmental stress in long-lived plant species. However, the speed and magnitude of global change impose serious questions about whether phenotypic plasticity alone will be enough to counter the negative effects of global change (Alberto et al., 2013). Adaptation of tree populations, which ultimately depends on intraspecific genetic variation in functional traits, may be the only way to cope with rapid environmental change (Aitken et al., 2008; Kremer et al., 2012; Alberto et al., 2013). While intraspecific genetic variation in conifers has been widely studied for traits related to growth and abiotic resistance (de la Mata et al., 2012; Gaspar et al., 2013; Grivet et al., 2013), much less is known about variation in anatomical defensive traits such as resin ducts.
Some studies have shown that resin duct features in conifers harbor high intraspecific variation between populations (O’Neill et al., 2002; Zas et al., 2015) as well as within-populations (King et al., 2011; Moreira et al., 2012; Moreira et al., 2015b). However little is yet known about the evolutionary origin of such variation or the adaptive value of this anatomical trait. Provenance trials can be useful to better understand the environmental drivers of genetic variation among populations. For instance O’Neill et al. (2002) reported that climate conditions and geographic variables at the origin of Picea sitchensis (Bong.) Carr. populations shaped genetic variation in resin duct characteristics. However, genetic variation includes not only adaptive variation resulting from diversifying selection in fitness-related traits, but also neutral variation due to processes such as genetic drift or gene flow (González-Martínez et al., 2006; Holderegger et al., 2006). It is noteworthy that most studies aiming to assess the adaptive value of quantitative traits by testing the existence of environmental gradients have failed to differentiate neutral from adaptive variation. Accounting for the neutral population genetic structure can, however, considerably change the interpretation of clinal patterns (Xia et al., 2018; López‐Goldar et al., 2019).
Here, we explored inter-population variation in constitutive resin duct characteristics (density, mean size, and percentage of conductive area covered by resin ducts) and their inducibility in Pinus pinaster. To determine the adaptive origin of such variation, we assessed the extent to which inter-population differentiation was explained by the climate conditions at the origin. To differentiate neutral from potentially adaptive genetic variation, and to avoid spurious correlations between traits and climate variables, we accounted for the population structure derived from neutral molecular markers in our models. To our knowledge, this is the first study that aims to explain inter-population variation in resin duct characteristics as a function of adaptation to local climate conditions. This knowledge will help to better understand the microevolution patterns of anatomical defenses in conifer species and their relevance in adaptive forest management under global change scenarios.
Materials and Methods
Study System
Maritime pine (Pinus pinaster Ait.) is a Mediterranean conifer species of high economic and ecological value in Southwest Europe and North Africa, inhabiting a wide range of heterogeneous environments, from wet-coastal to regions with prolonged summer drought (Figure 1A). P. pinaster populations are strongly differentiated for traits related to growth, abiotic stress tolerance, reproduction, and chemical defenses (Tapias et al., 2004; Arrabal et al., 2005; Correia et al., 2008; Santos-del-Blanco et al., 2012). Moreover, the demographic history of this species, with multiple isolated postglacial refugia as sources of migration routes, has led to high and geographically structured neutral genetic variation among populations (Bucci et al., 2007; Jaramillo-Correa et al., 2015). Resin duct features in this species have shown large intraspecific genetic variation among populations (Zas et al., 2015) and families within populations (Moreira et al., 2015b). However, the evolutionary drivers of such variation remain unexplored.
Figure 1 (A) Distribution map of Pinus pinaster (from EUFORGEN 2009, www.euforgen.org) showing the location of the 10 studied population (See Table S1). Colors indicate populations from different genetic clusters based on 1745 SNPs genotyped in 153 randomly sampled and unrelated individuals (from Jaramillo-Correa et al., 2015). (B) Principal component analysis of climate variables and plot of the 10 populations in the principal components ordination space. The climate variables included were annual mean temperature (AMT), maximum temperature of the warmest month (MTWM), minimum temperature of the coldest month (MTCM), temperature seasonality (TS), annual precipitation (AP), and potential evapotranspiration (ETP). Higher values of PC1 were indicative of Atlantic climates with higher precipitation and lower thermic oscillation, while lower values were related to Mediterranean climates. PC2 was positively associated with annual mean temperature (AMT) and potential evapotranspiration (ETP).
Plant Material, Experimental Design, and Sampling
To explore inter-population genetic variation in resin duct characteristics and their inducibility we used 10 populations of P. pinaster from CLONAPIN bank 1 collection (Majada et al., 2011). The populations were spread across the natural distribution range of the species, covering 10 latitude degrees from the southern to northern range ends, representing a wide environmental gradient. At each native population, five mother trees were selected and five open-pollinated seeds were collected from each of them. A total of 250 half-sib seeds with known mother and unknown father were grown under a common environment and comprised the ortet collection (10 population × 5 families × 5 genotypes) (López-Goldar et al., 2018). For the current experiment, two genotypes from 2 to 5 families within each of 10 the population were randomly selected and vegetatively propagated by minicuttings as described in Majada et al. (2011) obtaining two clonal copies, one exposed to methyl jasmonate treatment (MJ) and the other used as a control (CT). When cuttings were 2-years old they were relocated to a glasshouse in the MBG-CSIC (Pontevedra, Spain) and distributed following a split-plot design replicated in two blocks with population as the whole plot factor, and the factorial combination of family and MJ as the split factor. A total of N = 130 plants were included in the experiment (10 populations × 2–5 families × 2 blocks × 2 clonally propagated genotypes, each one subjected to either induction or control treatment). Resin duct induction was performed on plants allocated to the MJ treatment and sprayed with a solution of 25 mM of MJ (Sigma-Aldrich #39270-7) in deionized water with 2.5% ethanol (v/v). MJ was sprayed over the foliage to runoff (~10 ml per plant). Control plants were sprayed with the carrier solution in the same way (López-Goldar et al., 2018). Plants allotted to MJ treatment were transferred to a separate chamber for 6 days to avoid indirect induction effects on control plants by volatile emissions. All plants were watered weekly during the experiment and maintained in the same irrigation conditions to avoid any effects of differential water availability among plants.
One month after the induction treatment, all plants were harvested by cutting the stem aboveground. Then, a short section (ca. 2 cm long) of the main stem at mid height was sampled and stored in FAA (5% formaldehyde, 90% ethanol, 5% acetic acid) and then transferred to 75% ethanol (v/v) as previously reported in Moreira et al. (2015b). Thin transversal sections (40–60 µm) were obtained using a rotatory microtome (Leica RD 2235), stained with 0.01% safranin, and left overnight. Transversal sections were cleaned with progressive ethanol solutions (30, 50, 75, and 90%) for dehydration and to remove the excess of safranin. Finally, the stem sections were mounted on microscope slides and preserved in glycerol.
Resin Duct Measurements
Images from the transversal sections were obtained with a camera (VisiScope SCM154) coupled to a binocular stereoscope and analyzed using ImageJ 1.51k image analysis software (Abràmoff et al., 2004). After calibration with a precision slide (0.01 mm), we measured the area of phloem and xylem and separately counted the number of axial resin ducts in each tissue and measured their lumen area (Figure 2). Measurements were taken in two randomly selected quarters of the whole section. In each tissue we calculated the density of resin ducts as the number of resin ducts per square mm of tissue (hereafter RD density – number/mm2), the mean area of resin ducts (hereafter RD mean size – mm2), and the percentage of the total tissue area occupied by the lumen of resin ducts (hereafter RD conductive area – %) (Table 1). Resin duct measurements relative to tissue area represent an estimation of resin duct production in relation to plant growth.
Table 1 Description of resin duct characteristics measured in the xylem and the phloem of Pinus pinaster plants from 10 populations.
Figure 2 Examples of stem transversal sections of experimental plants showing axial resin ducts in the phloem (A) and xylem (B) in two year old Pinus pinaster saplings.
Statistical Analysis
To have an initial picture of the variation and covariation of the different analyzed resin duct metrics, main descriptive statistics of each variable and Pearson’s correlation coefficient among them were calculated. To explore inter-population genetic variation in resin ducts characteristics (RD density, RD mean size and RD conductive area in the xylem and phloem) and their response to MJ treatment (inducibility), we applied linear mixed models using the lmer function from lme4 package in R (Bates et al., 2014). Population (Pop), methyl jasmonate treatment (MJ), and their interaction (Pop × MJ) were included as fixed factors. Family within population (F[Pop]), the interaction between families and MJ (F[Pop] × MJ), and the genotype were included as random factors to reduce the model error rather than to explore genetic variation among families and genotypes. All resin duct variables were log-transformed to meet normality assumptions.
We obtained average climate data for the period 1950–2000 at the origin of populations from regional (Gonzalo, 2007) and global (Worldclim database - Hijmans et al., 2005) climate models. Both models compute gridded monthly climate data with a spatial resolution of 1 km2. The model of Gonzalo (2007) was used for Spanish populations due to its higher accuracy modeling climate variables in this region. Climate variables were summarized using principal component analysis as described in López‐Goldar et al. (2019). We used annual mean temperature (AMT), maximum temperature of the warmest month (MTWM), minimum temperature of the coldest month (MTCM), temperature seasonality (TS), annual precipitation (AP), and potential evapotranspiration (ETP). The first principal component (PC1) explained 53.4% of the total variance, with higher values associated with lower thermic oscillation and higher precipitation (Atlantic climates). In contrast, lower values of this index were associated with higher temperature extremes and lower precipitation (Mediterranean climates). The second principal component (PC2) explained 36% of the total variance and was positively correlated with mean annual temperature and potential evapotranspiration (see Figure 1B).
To test whether resin duct differentiation among populations was associated with climate, we performed correlation analysis between those resin duct characteristics that showed significant inter-population variation and the climate indices (PC1 and PC2). To that end, population least square means for each resin duct trait were obtained from the afore-mentioned model using the LSmeans function from the homonym package in R (Lenth, 2016).
Given the genetic population structure of the species, the non-independence between populations can generate spurious correlations (Desdevises et al., 2003). To account for this source of bias, we retrieved the population structure data (Q matrix) based on the genotyping of 126 SNPs from an independent set of clonal replicates (see more details in López-Goldar et al., 2018) publicly available in the Zenodo repository (http://doi.org/10.5281/zenodo.1445313). The Q matrix assigns to each genotype a genetic membership coefficient for each of the five genetic clusters (see Figure 1A) and was included as a covariate in the afore-mentioned model, as described in Yu et al. (2006). Correlation analyses with climate variables were performed once again, this time using the corrected population least square means for resin duct characteristics, adjusted for the neutral population structure.
Results
RD density in control and induced plants was higher in the xylem than in the phloem (2.60 ± 0.006 and 0.58 ± 0.002 RD/mm2 respectively; mean ± SE., N = 130). However, RD mean size was a 10-fold higher in the phloem than in the xylem (0.01 ± 0.00006 and 0.001 ± 0.000003 mm2 respectively; mean ± SE., N = 130). Consequently, RD conductive area was higher in the phloem than in the xylem (0.67 ± 0.003 and 0.25 ± 0.000007% respectively; mean ± SE, N = 130) (see also Table S2).
RD density was negatively correlated to RD mean size both in the xylem (r = ‑0.41, p < 0.001, N = 130) and in the phloem (r = ‑0.22, p < 0.05, N = 130). Moreover, RD density in the phloem was negatively correlated with RD mean size in the xylem (r = ‑0.3, p < 0.001, N = 130) and vice versa, RD density in the xylem was negatively correlated with RD mean size in the phloem (r = ‑0.24, p < 0.05, N = 130) (see Table S3).
Inter-Population Variation in Resin Duct Characteristics and Their Inducibility
RD density in the xylem and in the phloem showed large inter-population genetic variation (Table 2). Constitutive RD density in the xylem varied between 1.73 and 2.94 RD/mm2 between populations. TAMR was the population with the highest RD density in the xylem whereas PLEU showed the lowest value for this trait (Figure 3A). In the phloem, constitutive RD density ranged between 0.31 and 0.73 RD/mm2. PTOV and TAMR were the populations that had the highest and lowest value respectively for this trait (Figure 3B). Although no significant inter-population variation was found for RD mean size in the xylem and the phloem, RD conductive area was highly variable in both tissues (Table 2). Constitutive RD conductive area in the xylem ranged from 0.17 and 0.34 percent between populations, with TAMR and MIMI being the populations with the highest and lowest values respectively (Figure 3C). Constitutive RD conductive area in the phloem ranged between 0.25 and 0.82 percent between populations and TAMR showed the lowest value, whereas CDVO had the highest value for this trait (Figure 3D).
Table 2 Summary of the mixed model testing the effects of Population (Pop), induction by methyl jasmonate (MJ), and their interaction (Pop × MJ) on resin duct density (RD density—RD/mm2), mean size (RD mean size—mm2) and relative conductive area (RD conductive area—%) in each tissue (xylem and phloem) of Pinus pinaster saplings.
Figure 3 Intraspecific genetic variation in constitutive resin duct density (RD density, RD/mm2, panels A and B) and relative conductive area (RD conductive area, %, panels C and D) in the xylem (A, C) and phloem (B, D) among 10 Pinus pinaster populations covering a south–north latitudinal gradient. Bars represent population least square means ± SE extracted from the mixed model (Control plants; N = 63).
MJ significantly increased RD density in the xylem and decreased RD mean size in the phloem (Table 2; Figure 4), but had no effect on other resin duct characteristics (Table 2, Figure S1). We did not find inter-population variation in the inducibility of resin ducts for any of the analyzed traits (non-significant Pop × MJ interaction; Table 2).
Figure 4 Resin duct density (RD density, RD/mm2) in the xylem (A) and resin duct mean size (RD mean size, mm2) in the phloem (B) in MJ-induced (N = 67) and control plants (N = 63) of Pinus pinaster from 10 populations. Bars represent least square means ± SE extracted from the mixed model. Significant p-values computed from the corresponding mixed model are indicated with asterisks (p < .01**; p < .05*).
Inter-population variation in constitutive RD conductive area in the phloem was significantly and negatively correlated with that in the xylem (Figure 5A). A similar trend was observed for RD density although the correlation was not significant (Figure 5B).
Figure 5 Relationship between constitutive resin duct features in the xylem and phloem across 10 Pinus pinaster populations. (A) Resin duct conductive area (RD conductive area, %) and (B) resin duct density (RD density, RD/mm2). Each point represents the population least square means ± SE extracted from the mixed model. Pearson correlation coefficient (r) and associated p-values (p) are shown.
Clinal Variation With Climate of Resin Duct Characteristics
Population differentiation in several resin duct characteristics was associated with climate across their distribution range. Before accounting for the population genetic structure, constitutive RD density and RD conductive area in the phloem were positively correlated with the PC1. However, RD conductive area in the xylem was negatively correlated with this index (Table 3—Uncorrected). Thus, Atlantic populations had lower RD conductive area in the xylem (Figure 6A—Uncorrected) and higher RD density and RD conductive area in the phloem (Figures 6B, C—Uncorrected) than Mediterranean populations. We did not find significant correlations between resin duct characteristics and the PC2.
Table 3 Pearson correlation coefficient (r) between constitutive resin duct characteristics—resin duct density (RD density—RD/mm2) and relative conductive area (RD conductive area—%)—in the xylem and phloem and climate at the origin of 10 Pinus pinaster populations.
Figure 6 Relationship between climate indices (PC1 and PC2) and the constitutive resin duct characteristics—Resin duct density (RD density, RD/mm2, panel B) and/or conductive area (RD conductive area, %, panels A, C and D)—in the xylem (A) and phloem (B, C, D) across 10 Pinus pinaster populations. White circles represent uncorrected population least squared means, whereas dark circles represent corrected data after accounting for the population genetic structure in the corresponding mixed models. Higher values of PC1 are indicative of Atlantic climates with higher precipitation and lower thermic oscillation, whereas lower values are associated with Mediterranean climates. PC2 is positively associated with temperature and potential evapotranspiration. Pearson correlation coefficients (r) and associated p-values (p) are shown for uncorrected (gray) and corrected (black) data. Significant correlations are also represented with trend lines (dashed—uncorrected; solid—corrected).
After accounting for the population genetic structure only RD conductive area in the xylem remained negatively correlated with the PC1 (Table 3; Figure 6A—Corrected). RD density and RD conductive area in the phloem were no longer correlated with this climate index (Table 3; Figures 6B, C—Corrected). Although we didn’t find any significant correlations between resin duct characteristics and the PC2 before accounting for the population structure of the species, a significant positive correlation between RD conductive area in the phloem and the PC2 emerged after including the Q matrix in our models (Table 3; Figure 6D—Corrected).
Discussion
We found large inter-population variation in constitutive resin ducts, but no evidence of inter-population genetic variation in their inducibility. Of all the resin duct characteristics examined, only RD conductive area both in the xylem and in the phloem were correlated with climate at the origin of the populations after accounting for the neutral demographic history of the species. These results suggest that adaptive diversifying selection associated with heterogeneous climates may be contributing to differentiation between populations for these traits only.
Inter-Population Variation in Resin Duct Characteristics and Their Inducibility
Inter-population variation in defensive traits is expected because allocation to defenses is costly and resource availability varies across the distribution range of the species. Thus, plants must balance resource investment between different life functions (e.g. growth, reproduction, maintenance) (Herms and Mattson, 1992; Endara and Coley, 2011). Abiotic environmental conditions may select for divergent allocation priorities and thus, heterogeneous environments would favor population differentiation in terms of defensive investment. Supporting this prediction, we found that P. pinaster harbors large inter-population variation in the characteristics of their resin duct system, particularly for those metrics relative to secondary growth (RD density and RD conductive area). These results also agree with previous findings. In particular, Zas et al. (2015) found that Iberian populations of this species differ greatly in xylem RD density and cortical RD mean size. Likewise, O’Neill et al. (2002) found that constitutive resin duct characteristics (area and number) had high intraspecific genetic variation in Picea sp. Variation among Pinus sylvestris L. and P. nigra J.F. Arnold populations in resin duct characteristics has also been reported (Martin et al., 2010; Esteban et al., 2012), although these studies did not use common gardens and thus the observed variation cannot be exclusively linked to genetic (population) effects.
The observed variation in resin ducts among P. pinaster populations suggests important implications in potential biotic resistance in this species. It has been suggested that resource allocation to resin ducts provides enhanced resistance to different antagonists. For instance, in North American pine species, trees that survived or avoided bark beetle attacks had more or larger resin ducts than trees that died (Kane and Kolb, 2010; Ferrenberg et al., 2014; Gaylord et al., 2015; Hood et al., 2015; Zhao and Erbilgin, 2019). Furthermore, Picea sitchensis provenances previously characterized as resistant to the white pine weevil (Pissodes strobi) exhibited a denser network of resin ducts (King et al., 2011). Enhanced resistance to other antagonistic organisms has also been associated with higher production of resin ducts (Christiansen et al., 1999; Krokene et al., 2003), but not always (Zas et al., 2015). Higher resin duct production may lead to enhanced resistance by means of increasing terpene production, one of the major chemical defensive compounds in conifers (Raffa, 2014). Indeed, we found that the density of phloem resin ducts was positively correlated with the concentration of diterpenes in the same populations (r = 0.69; p-value < 0.05; N = 10), which were analyzed in previous work (López-Goldar et al., 2018; López‐Goldar et al., 2019). The association between resin duct production and chemical defenses highlights the relevance of resin ducts as an important defensive trait in conifers.
The induction of resin ducts can also be crucial for conifer resistance to insects and pathogens (Christiansen et al., 1999; Krokene et al., 2003). Here we found that MJ increased RD density in the xylem, probably by means of de novo production of traumatic resin duct, as expected according to previous studies in Pinus (Vivas et al., 2012; Moreira et al., 2015b) and Picea species (Franceschi et al., 2002; Martin et al., 2002; Erbilgin et al., 2006; Krokene et al., 2008). In contrast, the exogenous application of MJ reduced RD mean size in the phloem compared to control untreated plants. The observed reduction in RD mean size in the phloem was surprising and may be explained as a side effect of MJ caused by enhanced resin production in the epithelial cells. It is known that MJ increases the synthesis of resin defenses (Sampedro et al., 2011; López-Goldar et al., 2018) and resin flow (Franceschi et al., 2002). Since resin is synthesized in the epithelial cells surrounding the resin duct lumen (Fahn, 1979) increased resin production and accumulation after induction could lead to the swelling of the epithelial cells reducing the lumen area of resin ducts. Further research should test this idea, exploring the effects of MJ on the activation of epithelial cells in the cortical resin ducts. Moreover, the state of hydration of resin duct epithelial cells could be affecting their lumen area. However, as previously stated in the methods section, all plants were well watered and remained in the same conditions during the experiment. Thus, this factor was unlikely contributing to the observed variation between groups.
Given the impact of induced resin defenses on resistance (Krokene et al., 2003; López-Goldar et al., 2018), the costs of their production (Sampedro et al., 2011) and the heterogeneous abiotic and biotic conditions across the distribution range of the species, variation among populations in the plastic responses to biotic stimuli was also expected. Indeed, intraspecific genetic variation among families in the inducibility of resin ducts and chemical defenses was previously reported for P. pinaster (Sampedro et al., 2011; Moreira et al., 2015b). In contraposition, a parallel study using the same plant material as our study (López‐Goldar et al., 2019) found scarce intraspecific genetic variation among and within populations in the inducibility of plant secondary metabolites. Similarly, we did not find a significant interaction between pine populations and the induction treatment, indicating that all populations responded similarly to MJ. Our results may suggest that stabilizing selection has been operating on the inducibility of resin ducts leading to a highly conserved response among populations.
Constitutive RD conductive area in the xylem and the phloem were negatively correlated among populations pointing to a possible trade-off between the development of resin ducts in each tissue. While resin ducts in the xylem are the main resin reservoir during mature stages, resin ducts in the primary phloem are only present during the first years of development and are an important defensive mechanism in conifer seedlings (Franceschi et al., 2005). Theory suggests that redundant and costly defensive traits may trade-off among each other, although the evidences supporting this premise are scarce (Koricheva et al., 2004; Agrawal et al., 2010). Higher allocation priorities to phloem resin ducts could therefore lead to under defend the xylem, and vice versa. Differences in local environmental conditions could explain divergent defensive allocation priorities in each tissue among populations. For instance, populations exposed to high biotic stress may produce more resin ducts in the primary phloem in order to ensure tree survival at early developmental stages. However, our study lacks data on biotic pressures at the origin of populations and therefore this hypothesis cannot be tested. Further studies should explore the broad existence of this pattern in this and other species and assess the possible causes underlying this trade-off.
Adaptive Significance of Climate Gradients in Resin Duct Characteristics
Before accounting for the population genetic structure of the species, we found several significant correlations between resin duct characteristics in the phloem and xylem and the climate index PC1 (indicative of Atlantic climate). Phenotypic differences in resin duct characteristics among P. nigra provenances have also been associated with climate conditions in the Iberian Peninsula, with trees from colder environments producing more and larger radial resin ducts (Esteban et al., 2012). Likewise, differences in constitutive and traumatic resin duct characteristics among populations from the Picea sitchensis × Picea glauca (Moench) Voss introgression zone were found to be positively correlated with aridity and continentality at the origin of the populations (O’Neill et al., 2002).
According to the results obtained by this and other studies, it may be assumed that differential climate conditions act as a selective force contributing to the adaptive differentiation among populations (Alberto et al., 2013). However, inter-population variation does not only include adaptive but also neutral variation associated with the migratory and demographic history of the species (González‐Martínez et al., 2006; Holderegger et al., 2006). P. pinaster, like other European tree species, experienced severe bottlenecks during the ice ages (Gómez et al., 2005). This fact, together with its singular postglacial expansion from several relict emplacements, as well as its fragmented distribution, has led maritime pine to harbor large and geographically structured neutral genetic variation (Jaramillo-Correa et al., 2015; Serra‐Varela et al., 2015; López‐Goldar et al., 2019). As a result, the observed differentiation among populations in many phenotypic traits may not be related to enhanced tree fitness (i.e., adaptation). That seems to be the case for most plant secondary metabolites in P. pinaster, whose genetic differentiation among populations (QST) was not different from that expected by genetic drift alone (FST) (López‐Goldar et al., 2019).
When making inferences in comparative ecology, not accounting for the genetic relatedness can lead to the misinterpretation of potential clinal patterns due to the non-independence of the populations (Desdevises et al., 2003). After accounting for the population genetic structure we found that RD density in the phloem was no longer significantly correlated with the PC1 (Atlanticity). This finding indicates that the correlation between RD density in phloem and climate was probably not due to adaptive processes but rather caused by geographically structured patterns of neutral variation. Similar results were obtained for chemical defenses in this species, for which climate clines where reported to be due to neutral demographic processes in most cases (López‐Goldar et al., 2019).
The opposite situation would be those resin duct characteristics that were significantly correlated with climate even after accounting for the demographic history of the species. Such is the case of RD conductive area. In the xylem, this trait remained significantly and negatively correlated with the PC1 of Atlanticity after accounting for the neutral variation, suggesting that climate heterogeneity is contributing to inter-population differentiation of resin ducts in this tissue. In addition, a positive significant correlation between RD conductive area in the phloem and the PC2 emerged after accounting for the neutral variation. In both cases, adaptation to local climate conditions, alongside the demographic history of the species, seems to help explain the differentiation of populations in these traits. Current hypothesis on growth-defense trade-offs state that resource allocation to defenses is promoted in limiting environments to the detriment of growth capacity (Coley et al., 1985; Herms and Mattson, 1992; Endara and Coley, 2011). Our results somewhat support such a hypothesis, as populations from Atlantic regions characterized by less limiting conditions for tree growth invest less in resin duct production in the xylem than Mediterranean populations. Similarly, differentiation of populations towards higher allocation to phloem resin ducts is promoted in environments with higher temperature and potential evapotranspiration indicative of higher water stress. Moreover, we should acknowledge that resource allocation to defenses may also trade-off with investment in reproduction (Herms and Mattson, 1991; Redmond et al., 2019). Therefore, environmental conditions favoring higher reproductive investment could also limit defense production. RD conductive area is a more informative metric than RD density, as it accounts for both the density and the area of resin ducts in a tissue. Indeed, previous studies reported that a combination of smaller but denser resin ducts enhanced biotic resistance (Reed et al., 1986; Boucher et al., 2001). Accordingly, our results suggest that the relative conductive area of resin ducts has a greater adaptive and evolutionary significance than their relative number or their average size.
We must consider that separating neutral from adaptive variation in most European tree species may be challenging, as the postglacial migration history of these species parallels current environmental gradients (Serra‐Varela et al., 2015; López‐Goldar et al., 2019). Further studies in this and other conifer species should include comparative analysis of neutral versus quantitative variation (FST -QST analysis) in resin duct characteristics to broadly understand the evolutionary origin of this anatomical trait. Moreover, assessment of adaptive genetic variation in resin ducts should also be tested in mature trees, as previous studies showed that ontogeny plays an important role modulating resin-based defenses in conifers (Erbilgin and Colgan, 2012). We should also consider that maternal effects may be contributing to the observed differences among populations. Seeds were collected on natural populations growing in highly contrasted environments, which can lead to differential transgenerational plastic effects through mechanisms of non-genetic heritance (Yakovlev et al., 2012). However, the seeds collected in each population were initially planted in a common environment from which the cuttings used here were obtained. Although, common environmental conditions have likely helped to reduce such epigenetic effects, carryover and transgenerational effects could still be occurring. Moreover, heterogeneity in the biotic environment among populations (i.e., differential biotic pressure at the origin of populations) which is likely associated to climate gradients (Moreira et al., 2015a), could also contribute to explain patterns of genetic variation in this species. Unfortunately, we lack historical information on herbivory/pathogen pressures at the origin of the populations and thus this hypothesis could not be tested. We are aware that the sample size within populations in our study may limit our capacity to draw general assumptions. However, the strong relationship found between genetic differentiation in resin duct traits and climate still has great relevance in the theoretical framework of plant defense theory.
Altogether our findings show that even after accounting for the neutral variation in our models, differentiation between populations in resin duct characteristics was associated with climate gradients. Our results therefore are evidence that variation in the quantitative characteristics of this anatomical defensive trait in P. pinaster likely has an adaptive origin and suggests that resin duct characteristics have differentially evolved across heterogeneous climates. Although the inducibility of resin ducts is supposed to be an important mechanism conferring enhanced resistance, here we found no evidence of inter-population variation in their response to MJ. Thus, the magnitude of the biotic plasticity may not have been subject to diversifying selection, but rather widely conserved across the distribution range of P. pinaster. We aim to highlight the importance of accounting for the population genetic structure when inferring climate clines in phenotypic traits in species with large spatially structured neutral variation.
Data Availability Statement
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, performed the sampling, helped with the statistical analyses and the interpretation of the results, and improved the different versions of the manuscript. XL-G performed the sampling, produced the Q matrix, helped with the statistical analysis and the interpretation of the results, and contributed to improve the manuscript. CV-G processed the samples and obtained the images, measured resin duct features, performed the statistical analyses, produced the results, wrote the first draft with RZ, LS, and XL-G, and led the improvement of the manuscript along the peer-review process.
Funding
This research was supported by the Spanish government (MINECO/FEDER) grants FUTURPIN (AGL2015-68274-C03-02-R), AdapCon (CGL2011-30182-C02-01/02) and RESILPINE (RTI2018-094691-B-C33), and by the Xunta de Galicia-GAIN grant (IN607A2016/013). CVG received financial support from the MINECO FPI Grant program (MINECO-Spain BES-2016-076624). 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
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, Spain) and Dr. Ricardo Alía (INIA-CIFOR, Madrid, Spain) for conceiving and providing the plant material; José Carlos Rodríguez and Santiago Rodal (Viveros Figueirido, NORFOR, Pontevedra, Spain) for plant culture.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpls.2019.01613/full#supplementary-material
References
Abràmoff, M. D., Magalhães, P. J., Ram, S. J. (2004). Image processing with ImageJ. Biophotonics Int. 11, 36–42.
Agrawal, A. A., Conner, J. K., Rasmann, S. (2010). Tradeoffs and negative correlations in evolutionary ecology. In Evolution Since Darwin: The First 150 Years. (Sunderland Massachusetts: Sinauer Assoc. Inc.), pp. 243–268.
Aitken, S. N., Yeaman, S., Holliday, J. A., Wang, T., Curtis-McLane, S. (2008). Adaptation, migration or extirpation: climate change outcomes for tree populations. Evol. App. 1, 95–111. doi: 10.1111/j.1752-4571.2007.00013.x
Alberto, F. J., Aitken, S. N., Alía, R., González-Martínez, S. C., Hänninen, H., Kremer, A., et al. (2013). Potential for evolutionary responses to climate change–evidence from tree populations. Global Change Biol. 19, 1645–1661. doi: 10.1111/gcb.12181
Arrabal, C., Cortijo, M., de Simón, B. F., Vallejo, M. C. G., Cadahía, E. (2005). Differentiation among five Spanish Pinus pinaster provenances based on its oleoresin terpenic composition. Biochem. Syst. Ecol. 33, 1007–1016. doi: 10.1016/j.bse.2005.03.003
Bates, D., Mächler, M., Bolker, B., Walker, S. (2014). Fitting linear mixed-effects models using lme4. J. Stat. Software 67, 1–48. doi: 10.18637/jss.v067.i01
Boucher, D., Lavallee, R., Mauffette, Y. (2001). Biological performance of the white pine weevil in relation to the anatomy of the resin canal system of four different host species. Can. J. For. Res. 31, 2035–2041. doi: 10.1139/x01-135
Bucci, G., González-Martínez, 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 Chloroplast Microsatellite Markers. Mol. Ecol. 16, 2137–2153. doi: 10.1111/j.1365-294X.2007.03275.x
Chen, R., He, X., Chen, J., Gu, T., Liu, P., Xu, T., et al. (2018). Traumatic resin duct development, terpenoid formation, and related synthase gene expression in Pinus massoniana under feeding pressure of Monochamus alternatus. J. Plant Growth Regul. 38(3), 897–908. doi: 10.1007/s00344-018-9900-1
Christiansen, E., Krokene, P., Berryman, A. A., Franceschi, V. R., Krekling, T., Lieutier, F., et al. (1999). Mechanical injury and fungal infection induce acquired resistance in Norway spruce. Tree Physiol. 19, 399–403. doi: 10.1093/treephys/19.6.399
Coley, P. D., Bryant, J. P., Chapin, F. S. (1985). Resource availability and plant antiherbivore defense. Science 230, 895–899. doi: 10.1126/science.230.4728.895
Correia, I., Almeida, M. H., Aguiar, A., Alía, R., David, T. S., Pereira, J. S. (2008). Variations in growth, survival and carbon isotope composition (δ13C) among Pinus pinaster populations of different geographic origins. Tree Physiol. 28, 1545–1552. doi: 10.1093/treephys/28.10.1545
de la Mata, R., Voltas, J., Zas, R. (2012). Phenotypic plasticity and climatic adaptation in an Atlantic maritime pine breeding population. Ann. For. Sci. 69, 477–487. doi: 10.1007/s13595-011-0173-0
Desdevises, Y., Legendre, P., Azouzi, L., Morand, S. (2003). Quantifying phylogenetically structured environmental variation. Evolution 57, 2647–2652. doi: 10.1554/02-695
Early, R., Bradley, B. A., Dukes, J. S., Lawler, J. J., Olden, J. D., Blumenthal, D. M., et al. (2016). Global threats from invasive alien species in the twenty-first century and national response capacities. Nat. Commun. 7, 12485. doi: 10.1038/ncomms12485
Endara, M. J., Coley, P. D. (2011). The resource availability hypothesis revisited: a meta-analysis. Funct. Ecol. 25, 389–398. doi: 10.1111/j.1365-2435.2010.01803.x
Erbilgin, N., Colgan, L. J. (2012). Differential effects of plant ontogeny and damage type on phloem and foliage monoterpenes in jack pine (Pinus banksiana). Tree Physiol. 32, 946–957. doi: 10.1093/treephys/tps047
Erbilgin, N., Krokene, P., Christiansen, E., Zeneli, G., 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
Esteban, L. G., Martin, J. A., de Palacios, P., Fernandez, F. G. (2012). Influence of region of provenance and climate factors on wood anatomical traits of Pinus nigra Arn. Subsp. Salzmannii. Eur. J. For. Res. 131, 633–645.
Ferrenberg, S., Kane, J. M., Mitton, J. B. (2014). Resin duct characteristics associated with tree resistance to bark beetles across lodgepole and limber pines. Oecologia 174, 1283–1292. doi: 10.1007/s00442-013-2841-2
Franceschi, V. R., Krokene, P., Krekling, T., Christiansen, E. (2000). Phloem parenchyma cells are involved in local and distant defense responses to fungal inoculation or bark-beetle attack in Norway spruce (Pinaceae). Am. J. Bot. 87, 314–326
Franceschi, V. R., Krekling, T., Christiansen, E. (2002). Application of methyl jasmonate on Picea abies (Pinaceae) stems induces defense-related responses in phloem and xylem. Am. J. Bot. 89, 578–586. doi: 10.3732/ajb.89.4.578
Franceschi, V. R., Krokene, P., Christiansen, E., Krekling, T. (2005). Anatomical and chemical defenses of conifer bark against bark beetles and other pests. New Phytol. 167, 353–376. doi: 10.1111/j.1469-8137.2005.01436.x
Gómez, A., Vendramin, G., González-Martínez, S., Alía, R. (2005). Genetic diversity and differentiation of two Mediterranean pines (Pinus halepensis Mill. and Pinus pinaster Ait.) along a latitudinal cline using chloroplast microsatellite markers. Diversity Distrib. 11, 257–263. doi: 10.1111/j.1366-9516.2005.00152.x
Gaspar, M. J., Velasco, T., Feito, I., Alía, R., Majada, J. (2013). Genetic variation of drought tolerance in Pinus pinaster at three hierarchical levels: a comparison of induced osmotic stress and field testing. PloS One 8, e79094. doi: 10.1371/journal.pone.0079094
Gaylord, M. L., Kolb, T. E., McDowell, N. G. (2015). Mechanisms of pinon pine mortality after severe drought: A retrospective study of mature trees. Tree Physiol. 35, 806–816. doi: 10.1093/treephys/tpv038
González-Martínez, S. C., Krutovsky, K. V., Neale, D. B. (2006). Forest-tree population genomics and adaptive evolution. New Phytol. 170, 227–238. doi: 10.1111/j.1469-8137.2006.01686.x
Gonzalo, J. (2007). Phytoclimatic analysis of the Spanish Peninsula: update and geostatistical analysis. PhD Thesis, University of Valladolid, Palencia, Spain.
Grivet, D., Climent, J., Zabal-Aguirre, M., Neale, D. B., Vendramin, G. G., González-Martínez, S. C. (2013). Adaptive evolution of Mediterranean pines. Mol. Phylogenet. Evol. 68, 555–566. doi: 10.1016/j.ympev.2013.03.032
Herms, D. A., Mattson, W. J. (1991). “Does reproduction compromise defense in woody plants?” in Forest Insect Guilds: Patterns of Interaction with Host Trees. Eds. Baranchikov, Y. N., Mattson, W. J., Hain, F. P., Payne, T. L.; 1989 August 13-17; Abakan, Siberia, USSR Gen. Tech. Rep. NE-153. Radnor, PA: US Department of Agriculture, Forest Service, Northeastern Forest Experiment Station. 35–46. 153:35–46.
Herms, D. A., Mattson, W. J. (1992). The dilemma of plants: to grow or defend. Q. Rev. Biol. 67, 283–335.
Hijmans, R. J., Cameron, S. E., Parra, J. L., Jones, P. G., Jarvis, A. (2005). Very high resolution interpolated climate surfaces for global land areas. Int. J. Climatol. 25, 1965–1978.
Holderegger, R., Kamm, U., Gugerli, F. (2006). Adaptive vs. neutral genetic diversity: Implications for landscape genetics. Landscape Ecol. 21, 797–807. doi: 10.1007/s10980-005-5245-9
Hood, S., Sala, A. (2015). Ponderosa pine resin defenses and growth: metrics matter. Tree Physiol. 35, 1223–1235. doi: 10.1093/treephys/tpv098
Hood, S., Sala, A., Heyerdahl, E. K., Boutin, M. (2015). Low-severity fire increases tree defense against bark beetle attacks. Ecology 96, 1846–1855. doi: 10.1890/14-0487.1
Hudgins, J. W., Franceschi, V. R. (2004). Methyl jasmonate-induced ethylene production is responsible for conifer phloem defense responses and reprogramming of stem cambial zone for traumatic resin duct formation. Plant Physiol. 135, 2134–2149. doi: 10.1104/pp.103.037929
Jaramillo-Correa, J.-P., Rodríguez-Quilón, 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
Kane, J. M., Kolb, T. E. (2010). Importance of resin ducts in reducing ponderosa pine mortality from bark beetle attack. Oecologia 164, 601–609. doi: 10.1007/s00442-010-1683-4
Keeling, C. I., Bohlmann, J. (2006). Genes, enzymes and chemicals of terpenoid diversity in the constitutive and induced defence of conifers against insects and pathogens. New Phytol. 170, 657–675. doi: 10.1111/j.1469-8137.2006.01716.x
King, J. N., Alfaro, R. I., Lopez, M. G., Van Akker, L. (2011). Resistance of Sitka spruce (Picea sitchensis (Bong.) Carr.) to white pine weevil (Pissodes strobi Peck): characterizing the bark defence mechanisms of resistant populations. Forestry 84, 83–91. doi: 10.1093/forestry/cpq047
Koricheva, J., Nykänen, H., Gianoli, E. (2004). Meta-analysis of trade-offs among plant antiherbivore defenses: are plants jacks-of-all-trades, masters of all? Am. Natural. 163, E64–E75. doi: 10.1086/382601
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
Krokene, P., Solheim, H., Krekling, T., Christiansen, E. (2003). Inducible anatomical defense responses in Norway spruce stems and their possible role in induced resistance. Tree Physiol. 23, 191–197. doi: 10.1093/treephys/23.3.191
Krokene, P., Nagy, N. E., Solheim, H. (2008). Methyl jasmonate and oxalic acid treatment of Norway spruce: anatomically based defense responses and increased resistance against fungal infection. Tree Physiol. 28, 29–35. doi: 10.1093/treephys/28.1.29
López-Goldar, X., Villari, C., Bonello, P. E., Borg-Karlson, A. K., Grivet, D., Zas, R., et al. (2018). Inducibility of plant secondary metabolites predicts genetic variation in resistance against a key insect herbivore in maritime pine. Front. In Plant Sci. 9, 1651. doi: 10.3389/fpls.2018.01651
López-Goldar, X., Villari, C., Bonello, P., Borg-Karlson, A. K., Grivet, D., Sampedro, L., et al. (2019). Genetic variation in the constitutive defensive metabolome and its inducibility are geographically structured and largely determined by demographic processes in maritime pine. J. Ecol. 107, 2464–2477. doi: 10.1111/1365-2745.13159
Lenth, R. V. (2016). Least-squares means: the R package lsmeans. J. Stat. Software 69, 1–33. doi: 10.18637/jss.v069.i01
Lombardero, M., Ayres, M. P., Lorio, P. L., Jr., Ruel, J. 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
Majada, J., Martínez-Alonso, C., Feito, I., Kidelman, A., Aranda, I., Alía, R. (2011). Mini-cuttings: an effective technique for the propagation of Pinus pinaster Ait. New Forests 41, 399–412. doi: 10.1007/s11056-010-9232-x
Martin, D. M., Tholl, D., Gershenzon, J., Bohlmann, J. (2002). Methyl jasmonate induces traumatic resin ducts, terpenoid resin biosynthesis, and terpenoid accumulation in developing xylem of Norway spruce stems. Plant Physiol. 129, 1003–1018. doi: 10.1104/pp.011001
Martin, J. A., Esteban, L. G., de Palacios, P., Fernandez, F. G. (2010). Variation in wood anatomical traits of Pinus sylvestris L. between Spanish regions of provenance. Trees - Struct. Funct. 24, 1017–1028. doi: 10.1007/s00468-010-0471-4
Mason, C. J., Keefover-Ring, K., Villari, C., Klutsch, J. G., Cook, S., Bonello, P., et al. (2019). Anatomical defenses against bark beetles relate to degree of historical exposure between species and are allocated independently of chemical defenses within trees. Plant Cell Environ. 42, 633–646. doi: 10.1111/pce.13449
Moreira, X., Alfaro, R. I., King, J. N. (2012). Constitutive defenses and damage in Sitka spruce progeny obtained from crosses between white pine weevil resistant and susceptible parents. Forestry 85, 87–97. doi: 10.1093/forestry/cpr060
Moreira, X., Abdala-Roberts, L., Parra-Tabla, V., Mooney, K. A. (2015a). Latitudinal variation in herbivory: influences of climatic drivers, herbivore identity and natural enemies. Oikos 124, 1444–1452. doi: 10.1111/oik.02040
Moreira, X., Zas, R., Solla, A., Sampedro, L. (2015b). 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
Nagy, N. E., Franceschi, V. R., Solheim, H., Krekling, T., Christiansen, E. (2000). Wound-induced traumatic resin duct development in stems of Norway spruce (Pinaceae): Anatomy and cytochemical traits. Am. J. Bot. 87, 302–313.
O’Neill, G. A., Aitken, S. N., King, J. N., Alfaro, R. I. (2002). Geographic variation in resin canal defenses in seedlings from the Sitka spruce x white spruce introgression zone. Can. J. For. Res. 32, 390–400. doi: 10.1139/x01-206
Raffa, K. F., Aukema, B. H., Bentz, B. J., Carroll, A. L., Hicke, J. A., Turner, M. G., et al. (2008). Cross-scale drivers of natural disturbances prone to anthropogenic amplification: the dynamics of bark beetle eruptions. BioScience 58, 501–517. doi: 10.1641/B580607
Raffa, K. F. (2014). Terpenes tell different tales at different scales: Glimpses into the chemical ecology of conifer-bark beetle-microbial interactions. J. Chem. Ecol. 40, 1–20. doi: 10.1007/s10886-013-0368-y
Redmond, M. D., Davis, T. S., Ferrenberg, S. M., Wion, A. P. (2019). Resource allocation trade-offs in a mast-seeding conifer: Piñon pine prioritizes reproduction over defense. AoB Plants. 11, plz070. doi: 10.1093/aobpla/plz070
Reed, A. N., Hanover, J. W., Furniss, M. M. (1986). Douglas-fir and western larch: Chemical and physical properties in relation to Douglas-fir bark beetle attack. Tree Physiol. 1, 277–287. doi: 10.1093/treephys/1.3.277
Rodriguez-Garcia, A., Lopez, R., Martin, J. A., Pinillos, F., Gil, L. (2014). Resin yield in Pinus pinaster is related to tree dendrometry, stand density and tapping-induced systemic changes in xylem anatomy. For. Ecol. Manag. 313, 47–54. doi: 10.1016/j.foreco.2013.10.038
Sampedro, L., Moreira, X., Zas, R. (2011). 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
Santini, A., Ghelardini, L., De Pace, C., Desprez-Loustau, M.-L., Capretti, P., Chandelier, A., et al. (2013). Biogeographical patterns and determinants of invasion by forest pathogens in Europe. New Phytol. 197, 238–250. doi: 10.1111/j.1469-8137.2012.04364.x
Santos-del-Blanco, L., Climent, J., González-Martínez, S., Pannell, J. (2012). Genetic differentiation for size at first reproduction through male versus female functions in the widespread Mediterranean tree Pinus pinaster. Ann. Bot. 110, 1449–1460. doi: 10.1093/aob/mcs210
Schopmeyer, C., Mergen, F., Evans, T. C. (1954). Applicability of Poiseuille’s law to exudation of oleoresin from wounds on slash pine. Plant Physiol. 29, 82. doi: 10.1104/pp.29.1.82
Serra-Varela, M., Grivet, D., Vincenot, L., Broennimann, O., Gonzalo-Jiménez, J., Zimmermann, N. (2015). Does phylogeographical structure relate to climatic niche divergence? A test using maritime pine (Pinus pinaster Ait.). Global Ecol. Biogeogr. 24, 1302–1313. doi: 10.1111/geb.12369
Tapias, R., Climent, J., Pardos, J. A., Gil, L. (2004). Life histories of Mediterranean pines. Plant Ecol. 171, 53–68. doi: 10.1023/B:VEGE.0000029383.72609.f0
Vivas, M., Martin, J. A., Gil, L., Solla, A. (2012). Evaluating methyl jasmonate for induction of resistance to Fusarium oxysporum, F. circinatum and Ophiostoma novo-ulmi. For. Syst. 21, 289–299. doi: 10.5424/fs/2012212-02172
Whitehill, J. G., Yuen, M. M., Henderson, H., Madilao, L., Kshatriya, K., Bryan, J., et al. (2019). Functions of stone cells and oleoresin terpenes in the conifer defense syndrome. New Phytol. 221, 1503–1517. doi: 10.1111/nph.15477
Xia, H., Wang, B., Zhao, W., Pan, J., Mao, J. F., Wang, X. R. (2018). Combining mitochondrial and nuclear genome analyses to dissect the effects of colonization, environment, and geography on population structure in Pinus tabuliformis. Evol. App. 11, 1931–1945. doi: 10.1111/eva.12697
Yakovlev, I., Fossdal, C. G., Skrøppa, T., Olsen, J. E., Jahren, A. H., Johnsen, Ø. (2012). An adaptive epigenetic memory in conifers with important implications for seed production. Seed Sci. Res. 22, 63–76. doi: 10.1017/S0960258511000535
Yu, J., Pressoir, G., Briggs, W. H., Bi, I. V., 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. doi: 10.1038/ng1702
Zas, R., Moreira, X., Ramos, M., Lima, M. R. M., da Silva, M. N., 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 - Struct. Funct. 29, 663–673. doi: 10.1007/s00468-014-1143-6
Keywords: adaptive variation, climate gradients, conifer defences, genetic variation, neutral variation, Pinus pinaster, resin canals
Citation: Vázquez-González C, López-Goldar X, Zas R and Sampedro L (2019) Neutral and Climate-Driven Adaptive Processes Contribute to Explain Population Variation in Resin Duct Traits in a Mediterranean Pine Species. Front. Plant Sci. 10:1613. doi: 10.3389/fpls.2019.01613
Received: 14 September 2019; Accepted: 15 November 2019;
Published: 13 December 2019.
Edited by:
Henry D. Adams, Oklahoma State University, United StatesReviewed by:
Teemu Hölttä, University of Helsinki, FinlandRaimund Nagel, University of Leipzig, Germany
Miranda D. Redmond, Colorado State University, United States
Copyright © 2019 Vázquez-González, López-Goldar, 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: Carla Vázquez-González, Y3ZhenF1ZXpAbWJnLmNzaWMuZXM=