Phylogenomics reveals patterns of ancient hybridization and differential diversification that contribute to phylogenetic conflict in willows, poplars, and close relatives.

Despite the economic, ecological, and scientific importance of the genera Salix L. (willows) and Populus L. (poplars, cottonwoods, and aspens) Salicaceae, we know little about the sources of differences in species diversity between the genera and of the phylogenetic conflict that often confounds estimating phylogenetic trees. Salix subgenera and sections, in particular, have been difficult to classify, with one recent attempt termed a "spectacular failure" due to a speculated radiation of the subgenera Vetrix and Chamaetia. Here, we use targeted sequence capture to understand the evolutionary history of this portion of the Salicaceae plant family. Our phylogenetic hypothesis was based on 787 gene regions and identified extensive phylogenetic conflict among genes. Our analysis supported some previously described subgeneric relationships and confirmed the polyphyly of others. Using an fbranch analysis, we identified several cases of hybridization in deep branches of the phylogeny, which likely contributed to discordance among gene trees. In addition, we identified a rapid increase in diversification rate near the origination of the Vetrix-Chamaetia clade in Salix. This region of the tree coincided with several nodes that lacked strong statistical support, indicating a possible increase in incomplete lineage sorting due to rapid diversification. The extraordinary level of both recent and ancient hybridization in both Salix and Populus have played important roles in the diversification and diversity in these two genera.

Growth-defense trade-offs shape population genetic composition in an iconic forest tree species.

All organisms experience fundamental conflicts between divergent metabolic processes. In plants, a pivotal conflict occurs between allocation to growth, which accelerates resource acquisition, and to defense, which protects existing tissue against herbivory. Trade-offs between growth and defense traits are not universally observed, and a central prediction of plant evolutionary ecology is that context-dependence of these trade-offs contributes to the maintenance of intraspecific variation in defense [Züst and Agrawal, Annu. Rev. Plant Biol., 68, 513-534 (2017)]. This prediction has rarely been tested, however, and the evolutionary consequences of growth-defense trade-offs in different environments are poorly understood, especially in long-lived species [Cipollini et al., Annual Plant Reviews (Wiley, 2014), pp. 263-307]. Here we show that intraspecific trait trade-offs, even when fixed across divergent environments, interact with competition to drive natural selection of tree genotypes corresponding to their growth-defense phenotypes. Our results show that a functional trait trade-off, when coupled with environmental variation, causes real-time divergence in the genetic architecture of tree populations in an experimental setting. Specifically, competitive selection for faster growth resulted in dominance by fast-growing tree genotypes that were poorly defended against natural enemies. This outcome is a signature example of eco-evolutionary dynamics: Competitive interactions affected microevolutionary trajectories on a timescale relevant to subsequent ecological interactions [Brunner et al., Funct. Ecol. 33, 7-12 (2019)]. Eco-evolutionary drivers of tree growth and defense are thus critical to stand-level trait variation, which structures communities and ecosystems over expansive spatiotemporal scales.

Chemical Composition of Essential Oils from Natural Populations of Artemisia scoparia Collected at Different Altitudes: Antibacterial, Mosquito Repellent, and Larvicidal Effects.

The current study aimed to evaluate the presence of chemical variations in essential oils (EOs) extracted from Artemisia scoparia growing at different altitudes and to reveal their antibacterial, mosquito larvicidal, and repellent activity. The gas chromatographic-mass spectrometric analysis of A. scoparia EOs revealed that the major compounds were capillene (9.6-31.8%), methyleugenol (0.2-26.6%), ß-myrcene (1.9-21.4%), γ-terpinene (1.5-19.4%), trans-ß-caryophyllene (0.8-12.4%), and eugenol (0.1-9.1%). The EO of A. scoparia collected from the city of Attock at low elevation was the most active against Staphylococcus aureus, Bacillus subtilis, Escherichia coli, and Pseudomonas aeruginosa bacteria (minimum inhibitory concentration of 156-1250 µg/mL) and showed the best mosquito larvicidal activity (LC50, 55.3 mg/L). The EOs of A. scoparia collected from the high-altitude areas of Abbottabad and Swat were the most repellent for females of Ae. aegypti and exhibited repellency for 120 min and 165 min, respectively. The results of the study reveal that different climatic conditions and altitudes have significant effects on the chemical compositions and the biological activity of essential oils extracted from the same species.

A General Model to Explain Repeated Turnovers of Sex Determination in the Salicaceae.

Dioecy, the presence of separate sexes on distinct individuals, has evolved repeatedly in multiple plant lineages. However, the specific mechanisms by which sex systems evolve and their commonalities among plant species remain poorly understood. With both XY and ZW sex systems, the family Salicaceae provides a system to uncover the evolutionary forces driving sex chromosome turnovers. In this study, we performed a genome-wide association study to characterize sex determination in two Populus species, P. euphratica and P. alba. Our results reveal an XY system of sex determination on chromosome 14 of P. euphratica, and a ZW system on chromosome 19 of P. alba. We further assembled the corresponding sex-determination regions, and found that their sex chromosome turnovers may be driven by the repeated translocations of a Helitron-like transposon. During the translocation, this factor may have captured partial or intact sequences that are orthologous to a type-A cytokinin response regulator gene. Based on results from this and other recently published studies, we hypothesize that this gene may act as a master regulator of sex determination for the entire family. We propose a general model to explain how the XY and ZW sex systems in this family can be determined by the same RR gene. Our study provides new insights into the diversification of incipient sex chromosomes in flowering plants by showing how transposition and rearrangement of a single gene can control sex in both XY and ZW systems.

Genetic mapping of sexually dimorphic volatile and non-volatile floral secondary chemistry of a dioecious willow.

Secondary chemistry often differs between sexes in dioecious plant species, a pattern attributed to its possible role in the evolution and/or maintenance of dioecy. We used GC-MS to measure floral volatiles emitted from, and LC-MS to quantitate non-volatile secondary compounds contained in, female and male Salix purpurea willow catkins from an F2 family. Using the abundance of these chemicals, we then performed quantitative trait locus (QTL) mapping to locate them on the genome, identified biosynthetic candidate genes in the QTL intervals, and examined expression patterns of candidate genes using RNA-seq. Male flowers emitted more total terpenoids than females, but females produced more benzenoids. Male tissue contained greater amounts of phenolic glycosides, but females had more chalcones and flavonoids. A flavonoid pigment and a spermidine derivative were found only in males. Male catkins were almost twice the mass of females. Forty-two QTL were mapped for 25 chemical traits and catkin mass across 16 of the 19 S. purpurea chromosomes. Several candidate genes were identified, including a chalcone isomerase associated with seven compounds. A better understanding of the genetic basis of the sexually dimorphic chemistry of a dioecious species may shed light on how chemically mediated ecological interactions may have helped in the evolution and maintenance of dioecy.

Population-specific responses of floral volatiles to abiotic factors in changing environments.

PREMISE: Shifts in abiotic factors can affect many plant traits, including floral volatiles. This study examined the response of floral volatiles to water availability and whether phenotypic plasticity to water availability differs among populations. It also investigated genetic differentiation in floral volatiles, determined the effect of temperature on phenotypic plasticity to water availability, and assessed temporal variation in floral scent emission between day and evening, since pollinator visitation differs at those times. METHODS: Rocky Mountain columbine plants (Aquilegia coerulea), started from seeds collected in three wild populations in Colorado, Utah, and Arizona, were grown under two water treatments in a greenhouse in Madison, Wisconsin, United States. One population was also grown under the two water treatments, at two temperatures. Air samples were collected from enclosed flowers using dynamic headspace methods and floral volatiles were identified and quantified by gas chromatography (GC) with mass spectrometry (MS). RESULTS: Emission of three floral volatiles increased in the wetter environment, indicating phenotypic plasticity. The response of six floral volatiles to water availability differed among populations, suggesting genetic differentiation in phenotypic plasticity. Five floral volatiles varied among populations, and emission of most floral volatiles was greater during the day. CONCLUSIONS: Phenotypic plasticity to water availability permits a quick response of floral volatiles in changing environments. The genetic differentiation in phenotypic plasticity suggests that phenotypic plasticity can evolve but complicates predictions of the effects of environmental changes on a plant and its pollinators.

Trait plasticity and trade-offs shape intra-specific variation in competitive response in a foundation tree species.

The ability to tolerate neighboring plants (i.e. degree of competitive response) is a key determinant of plant success in high-competition environments. Plant genotypes adjust their functional trait expression under high levels of competition, which may help explain intra-specific variation in competitive response. However, the relationships between traits and competitive response are not well understood, especially in trees. In this study, we investigated among-genotype associations between tree trait plasticity and competitive response. We manipulated competition intensity in experimental stands of trembling aspen (Populus tremuloides) to address the covariance between competition-induced changes in functional trait expression and aspects of competitive ability at the genotype level. Genotypic variation in the direction and magnitude of functional trait responses, especially those of crown foliar mass, phytochemistry, and leaf physiology, was associated with genotypic variation in competitive response. Traits exhibited distinct plastic responses to competition, with varying degrees of genotypic variation and covariance with other trait responses. The combination of genotypic diversity and covariance among functional traits led to tree responses to competition that were coordinated among traits yet variable among genotypes. Such relationships between tree traits and competitive success have the potential to shape stand-level trait distributions over space and time.

Combined drought and bark beetle attacks deplete non-structural carbohydrates and promote death of mature pine trees.

How carbohydrate reserves in conifers respond to drought and bark beetle attacks are poorly understood. We investigated changes in carbohydrate reserves and carbon-dependent diterpene defences in ponderosa pine trees that were experimentally subjected to two levels of drought stress (via root trenching) and two types of biotic challenge treatments (pheromone-induced bark beetle attacks or inoculations with crushed beetles that include beetle-associated fungi) for two consecutive years. Our results showed that trenching did not influence carbohydrates, whereas both biotic challenges reduced amounts of starch and sugars of trees. However, only the combined trenched-bark beetle attacked trees depleted carbohydrates and died during the first year of attacks. While live trees contained higher carbohydrates than dying trees, amounts of constitutive and induced diterpenes produced did not vary between live and beetle-attacked dying trees, respectively. Based on these results we propose that reallocation of carbohydrates to diterpenes during the early stages of beetle attacks is limited in drought-stricken trees, and that the combination of biotic and abiotic stress leads to tree death. The process of tree death is subsequently aggravated by beetle girdling of phloem, occlusion of vascular tissue by bark beetle-vectored fungi, and potential exploitation of host carbohydrates by bark beetle symbionts as nutrients.

Evolutionary history and ecology shape the diversity and abundance of phytochemical arsenals across monkeyflowers.

We examine the extent to which phylogenetic effects and ecology are associated with macroevolutionary patterns of phytochemical defence production across the Mimulus phylogeny. We grew plants from 21 species representing the five major sections of the Mimulus phylogeny in a common garden to assess how the arsenals (NMDS groupings) and abundances (concentrations) of a phytochemical defence, phenylpropanoid glycosides (PPGs), vary across the phylogeny. Very few PPGs are widespread across the genus, but many are common to multiple sections of the genus. Phytochemical arsenals cluster among sections in an NMDS and are not associated with total concentration of PPGs. There is a strong phylogenetic signal for phytochemical arsenal composition across the Mimulus genus, whereas ecological variables such as growing season length, latitude, and elevation do not significantly influence arsenal. In contrast, there is little phylogenetic signal for total PPG concentration, and this trait is significantly influenced by several ecological factors. Phytochemical arsenals and abundances are influenced by plant life history form. Both phylogenetic effects and ecology are related to phytochemical patterns across species, albeit in different ways. The independence of phytochemical defence concentrations from arsenal compositions indicates that these aspects of defence may continue to evolve independently of one another.

Sex determination through X-Y heterogamety in Salix nigra.

The development of non-recombining sex chromosomes has radical effects on the evolution of discrete sexes and sexual dimorphism. Although dioecy is rare in plants, sex chromosomes have evolved repeatedly throughout the diversification of angiosperms, and many of these sex chromosomes are relatively young compared to those found in vertebrates. In this study, we designed and used a sequence capture array to identify a novel sex-linked region (SLR) in Salix nigra, a basal species in the willow clade, and demonstrated that this species has XY heterogamety. We did not detect any genetic overlap with the previously characterized ZW SLRs in willows, which map to a different chromosome. The S. nigra SLR is characterized by strong recombination suppression across a 2 MB region and an excess of low-frequency alleles, resulting in a low Tajima's D compared to the remainder of the genome. We speculate that either a recent bottleneck in population size or factors related to positive or background selection generated this differential pattern of Tajima's D on the X and autosomes. This discovery provides insights into factors that may influence the evolution of sex chromosomes in plants and contributes to a large number of recent observations that underscore their dynamic nature.

Phylogenomics of the genus Populus reveals extensive interspecific gene flow and balancing selection.

Phylogenetic analysis is complicated by interspecific gene flow and the presence of shared ancestral polymorphisms, particularly those maintained by balancing selection. In this study, we aimed to examine the prevalence of these factors during the diversification of Populus, a model tree genus in the Northern Hemisphere. We constructed phylogenetic trees of 29 Populus taxa using 80 individuals based on re-sequenced genomes. Our species tree analyses recovered four main clades in the genus based on consensus nuclear phylogenies, but in conflict with the plastome phylogeny. A few interspecific relationships remained unresolved within the multiple-species clade because of inconsistent gene trees. Our results indicated that gene flow has been widespread within each clade and also occurred among the four clades during their early divergence. We identified 45 candidate genes with ancient polymorphisms maintained by balancing selection. These genes were mainly associated with mating compatibility, growth and stress resistance. Both gene flow and selection-mediated ancient polymorphisms are prevalent in the genus Populus. These are potentially important contributors to adaptive variation. Our results provide a framework for the diversification of model tree genus that will facilitate future comparative studies.

To compete or defend: linking functional trait variation with life-history tradeoffs in a foundation tree species.

Although chemical deterrents to herbivory often exact costs in terms of plant growth, the manner in which those costs arise, and their physiological relationship to other functional traits, remain unclear. In the absence of appreciable herbivory, we examined interrelationships among chemical defense levels and other foliar functional traits (e.g., light-saturated photosynthesis, specific leaf area, nitrogen concentration) as co-determinants of tree growth and, by extension, competitive ability in high-density populations comprising 16 genotypes of Populus tremuloides. Across genotypes, concentrations of chemical defenses were not significantly related to other leaf functional traits, but levels of the salicinoid phenolic glycosides (SPGs) salicin, salicortin and tremulacin were each negatively correlated with relative mass growth (RMG) of aboveground woody tissue (P ≤ 0.001). RMG, in turn, underpinned 77% of the genotypic variation in relative height growth (our index of competitive ability). RMG was also positively related to light-saturated photosynthesis (P ≤ 0.001), which, together with the three SPGs, explained 86% of genotypic RMG variation (P ≤ 0.001). Moreover, results of a carbon balance simulation indicated that costs of resource allocation to SPGs, reaching nearly a third of annual crown photosynthesis, were likely mediated by substantial metabolic turnover, particularly for salicin. The lack of discernible links between foliar defense allocation and other (measured) functional traits, and the illustrated potential of metabolic turnover to reconcile influences of SPG allocation on RMG, shed additional light on fundamental physiological mechanisms underlying evolutionary tradeoffs between chemical defense investment and competitive ability in a foundation tree species.

Anatomical defences against bark beetles relate to degree of historical exposure between species and are allocated independently of chemical defences within trees.

Conifers possess chemical and anatomical defences against tree-killing bark beetles that feed in their phloem. Resins accumulating at attack sites can delay and entomb beetles while toxins reach lethal levels. Trees with high concentrations of metabolites active against bark beetle-microbial complexes, and more extensive resin ducts, achieve greater survival. It is unknown if and how conifers integrate chemical and anatomical components of defence or how these capabilities vary with historical exposure. We compared linkages between phloem chemistry and tree ring anatomy of two mountain pine beetle hosts. Lodgepole pine, a mid-elevation species, has had extensive, continual contact with this herbivore, whereas high-elevation whitebark pines have historically had intermittent exposure that is increasing with warming climate. Lodgepole pine had more and larger resin ducts. In both species, anatomical defences were positively related to tree growth and nutrients. Within-tree constitutive and induced concentrations of compounds bioactive against bark beetles and symbionts were largely unrelated to resin duct abundance and size. Fewer anatomical defences in the semi-naïve compared with the continually exposed host concurs with directional differences in chemical defences. Partially uncoupling chemical and morphological antiherbivore traits may enable trees to confront beetles with more diverse defence permutations that interact to resist attack.

Drought-Mediated Changes in Tree Physiological Processes Weaken Tree Defenses to Bark Beetle Attack.

Interactions between water stress and induced defenses and their role in tree mortality due to bark beetles are poorly understood. We performed a factorial experiment on 48 mature ponderosa pines (Pinus ponderosa) in northern Arizona over three years that manipulated a) tree water stress by cutting roots and removing snow; b) bark beetle attacks by using pheromone lures; and c) phloem exposure to biota vectored by bark beetles by inoculating with dead beetles. Tree responses included resin flow from stem wounds, phloem composition of mono- and sesqui-terpenes, xylem water potential, leaf gas exchange, and survival. Phloem contained 21 mono- and sesqui-terpenes, which were dominated by (+)-α-pinene, (-)-limonene, and δ-3-carene. Bark beetle attacks (mostly Dendroctonus brevicomis) and biota carried by beetles induced a general increase in concentration of phloem mono- and sesqui-terpenes, whereas water stress did not. Bark beetle attacks induced an increase in resin flow for unstressed trees but not water-stressed trees. Mortality was highest for beetle-attacked water-stressed trees. Death of beetle-attacked trees was preceded by low resin flow, symptoms of water stress (low xylem water potential, leaf gas exchange), and an ephemeral increase in concentrations of mono- and sesqui-terpenes compared to surviving trees. These results show a) that ponderosa pine can undergo induction of both resin flow and phloem terpenes in response to bark beetle attack, and that the former is more constrained by water stress; b) experimental evidence that water stress predisposes ponderosa pines to mortality from bark beetles.

Genetic variation in aspen phytochemical patterns structures windows of opportunity for gypsy moth larvae.

Empirical studies indicate that host-tree bud break will likely advance faster than spring-folivore egg hatch in response to predicted increases in temperature. How these phenological shifts will affect herbivory will depend on temporal patterns of foliar traits that occur during leaf expansion, and their effects on folivore performance. Through fine-scale time series sampling of newly flushed trembling aspen (Populus tremuloides) foliage, we observed a previously unknown peak in phenolic glycoside concentrations that coincides with the emergence of sensitive neonates of gypsy moths and rapidly declines soon after bud break. The magnitude and duration of the initial post-bud break peak in phenolic glycosides varied substantially among genotypes. In contrast, foliar nitrogen concentrations declined at a more uniform rate among genotypes throughout leaf expansion. In addition, leaf toughness remained uniformly low throughout these periods of phytochemical change, and did not rise or vary substantially among genotypes until after anticipated windows of climate change-induced shifts between bud break and egg hatch had elapsed. Controlled manipulation of intervals between gypsy moth egg hatch and aspen bud break generated differences in larval performance among hatch cohorts and host genotypes that corresponded with changes in foliar phenolic glycoside and nitrogen concentrations. These findings indicate that the effects of climate change-induced phenological shifts on herbivory will differ among host plant genotypes, and that genetic variation in foliar chemical patterns will strongly influence this heterogeneity.

Defence syndromes in lodgepole - whitebark pine ecosystems relate to degree of historical exposure to mountain pine beetles.

Warming climate is allowing tree-killing bark beetles to expand their ranges and access naïve and semi-naïve conifers. Conifers respond to attack using complex mixtures of chemical defences that can impede beetle success, but beetles exploit some compounds for host location and communication. Outcomes of changing relationships will depend on concentrations and compositions of multiple host compounds, which are largely unknown. We analysed constitutive and induced chemistries of Dendroctonus ponderosae's primary historical host, Pinus contorta, and Pinus albicaulis, a high-elevation species whose encounters with this beetle are transitioning from intermittent to continuous. We quantified multiple classes of terpenes, phenolics, carbohydrates and minerals. Pinus contorta had higher constitutive allocation to, and generally stronger inducibility of, compounds that resist these beetle-fungal complexes. Pinus albicaulis contained higher proportions of specific monoterpenes that enhance pheromone communication, and lower induction of pheromone inhibitors. Induced P. contorta increased insecticidal and fungicidal compounds simultaneously, whereas P. albicaulis responses against these agents were inverse. Induced terpene accumulation was accompanied by decreased non-structural carbohydrates, primarily sugars, in P. contorta, but not P. albicaulis, which contained primarily starches. These results show some host species with continuous exposure to bark beetles have more thoroughly integrated defence syndromes than less-continuously exposed host species.

Rapid modulation of ultraviolet shielding in plants is influenced by solar ultraviolet radiation and linked to alterations in flavonoids.

The accumulation of ultraviolet (UV)-absorbing compounds (flavonoids and related phenylpropanoids) and the resultant decrease in epidermal UV transmittance (TUV ) are primary protective mechanisms employed by plants against potentially damaging solar UV radiation and are critical components of the overall acclimation response of plants to changing solar UV environments. Whether plants can adjust this UV sunscreen protection in response to rapid changes in UV, as occurs on a diurnal basis, is largely unexplored. Here, we use a combination of approaches to demonstrate that plants can modulate their UV-screening properties within minutes to hours, and these changes are driven, in part, by UV radiation. For the cultivated species Abelmoschus esculentus, large (30-50%) and reversible changes in TUV occurred on a diurnal basis, and these adjustments were associated with changes in the concentrations of whole-leaf UV-absorbing compounds and several quercetin glycosides. Similar results were found for two other species (Vicia faba and Solanum lycopersicum), but no such changes were detected in Zea mays. These findings reveal a much more dynamic UV-protection mechanism than previously recognized, raise important questions concerning the costs and benefits of UV-protection strategies in plants and have practical implications for employing UV to enhance crop vigor and quality in controlled environments.

Rapid Induction of Multiple Terpenoid Groups by Ponderosa Pine in Response to Bark Beetle-Associated Fungi.

Ponderosa pine (Pinus ponderosa) is a major and widely distributed component of conifer biomes in western North America and provides substantial ecological and economic benefits. This tree is exposed to several tree-killing bark beetle-microbial complexes, including the mountain pine beetle (Dendroctonus ponderosae) and the phytopathogenic fungus Grosmannia clavigera that it vectors, which are among the most important. Induced responses play a crucial role in conifer defenses, yet these have not been reported in ponderosa pine. We compared concentrations of terpenes and a phenylpropanoid, two phytochemical classes with strong effects against bark beetles and their symbionts, in constitutive phloem tissue and in tissue following mechanical wounding or simulated D. ponderosae attack (mechanical wounding plus inoculation with G. clavigera). We also tested whether potential induced responses were localized or systemic. Ponderosa pines showed pronounced induced defenses to inoculation, increasing their total phloem concentrations of monoterpenes 22.3-fold, sesquiterpenes 56.7-fold, and diterpenes 34.8-fold within 17 days. In contrast, responses to mechanical wounding alone were only 5.2, 11.3, and 7.7-fold, respectively. Likewise, the phenylpropanoid estragole (4-allyanisole) rose to 19.1-fold constitutive levels after simulated attack but only 4.4-fold after mechanical wounding. Overall, we found no evidence of systemic induction after 17 days, which spans most of this herbivore's narrow peak attack period, as significant quantitative and compositional changes within and between terpenoid groups were localized to the wound site. Implications to the less frequent exploitation of ponderosa than lodgepole pine by D. ponderosae, and potential advantages of rapid localized over long-term systemic responses in this system, are discussed.

Patterns of phytochemical variation in Mimulus guttatus (yellow monkeyflower).

The search for general patterns in the production and allocation of plant defense traits will be facilitated by characterizing multivariate suites of defense, as well as by studying additional plant taxa, particularly those with available genomic resources. Here, we investigated patterns of genetic variation in phytochemical defenses (phenylpropanoid glycosides, PPGs) in Mimulus guttatus (yellow monkeyflower). We grew plants derived from several natural populations, consisting of multiple full-sibling families within each population, in a common greenhouse environment. We found substantial variation in the constitutive multivariate PPG phenotype and in constitutive levels of individual phytochemicals within plants (among leaves of different ages), within populations (among full-sibling families), and among populations. Populations consisting of annual plants generally, but not always, had lower concentrations of phytochemicals than did populations of perennial plants. Populations differed in their plastic response to artificial herbivory, both in the overall multivariate PPG phenotype and in the individual phytochemicals. The relationship between phytochemistry and another defense trait, trichomes, differed among populations. Finally, we demonstrated that one of the PPGs, verbascoside, acts as a feeding stimulant rather than a feeding deterrent for a specialist herbivore of M. guttatus, the buckeye caterpillar (Junonia coenia Nymphalidae). Given its available genetic resources, numerous, easily accessible natural populations, and patterns of genetic variation highlighted in this research, M. guttatus provides an ideal model system in which to test ecological and evolutionary theories of plant-herbivore interactions.

Bacteria associated with a tree-killing insect reduce concentrations of plant defense compounds.

Bark beetles encounter a diverse array of constitutive and rapidly induced terpenes when attempting to colonize living conifers. Concentrations of these compounds at entry sites can rapidly reach levels toxic to beetles, their brood, and fungal symbionts. Large numbers of beetles can overwhelm tree defenses via pheromone-mediated mass attacks, but the mechanisms are poorly understood. We show that bacteria associated with mountain pine beetles can metabolize monoterpenes and diterpene acids. The abilities of different symbionts to reduce concentrations of different terpenes appear complementary. Serratia reduced concentrations of all monoterpenes applied to media by 55-75 %, except for α-pinene. Beetle-associated Rahnella reduced (-)- and (+)-α-pinene by 40 % and 45 %, respectively. Serratia and Brevundimonas reduced diterpene abietic acid levels by 100 % at low concentrations. However, high concentrations exhausted this ability, suggesting that opposing rates of bacterial metabolism and plant induction of terpenes are critical. The two major fungal symbionts of mountain pine beetle, Grosmannia clavigera and Ophiostoma montium were highly susceptible to abietic acid. Grosmannia clavigera did not reduce total monoterpene concentrations in lodgepole pine turpentine. We propose the ability of bark beetles to exert landscape-scale impacts may arise partly from micro-scale processes driven by bacterial symbionts.