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.

Environmental clines in organismal defensive traits are usually attributed to stronger selection by enemies at lower latitudes or near the host’s range center. Nonetheless, little functional evidence has supported this hypothesis, especially for coevolving plants and herbivores. We quantified cardenolide toxins in seeds of 24 populations of common milkweed ( Asclepias syriaca ) across 13 degrees of latitude, revealing a pattern of increasing cardenolide concentrations toward the host's range center. The unusual nitrogen-containing cardenolide labriformin was an exception and peaked at higher latitudes. Milkweed seeds are eaten by specialist lygaeid bugs that are even more tolerant of cardenolides than the monarch butterfly, concentrating most cardenolides (but not labriformin) from seeds into their bodies. Accordingly, whether cardenolides defend seeds against these specialist bugs is unclear. We demonstrate that Oncopeltus fasciatus (Lygaeidae) metabolized two major compounds (glycosylated aspecioside and labriformin) into distinct products that were sequestered without impairing growth. We next tested several isolated cardenolides in vitro on the physiological target of cardenolides (Na ⁺ /K ⁺ -ATPase); there was little variation among compounds in inhibition of an unadapted Na ⁺ /K ⁺ -ATPase, but tremendous variation in impacts on that of monarchs and Oncopeltu s. Labriformin was the most inhibitive compound tested for both insects, but Oncopeltus had the greater advantage over monarchs in tolerating labriformin compared to other compounds. Three metabolized (and stored) cardenolides were less toxic than their parent compounds found in seeds. Our results suggest that a potent plant defense is evolving by natural selection along a geographical cline and targets specialist herbivores, but is met by insect tolerance, detoxification, and sequestration.

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.