Spatial genetic structure of the mountain pine beetle (Dendroctonus ponderosae) outbreak in western Canada: historical patterns and contemporary dispersal

Rearranging and completing the puzzle: Phylogenomic analysis of bark beetles Dendroctonus reveals new hypotheses about genus diversification

Studies carried out on bark beetles within Dendroctonus have been extensive and revealed diverse information in different areas of their natural history, taxonomy, evolution, and interactions, among others. Despite these efforts, phylogenetic hypotheses have remained obscured mainly due to limited information analyzed (taxonomic, gene sampling, or both) in studies focused on obtaining evolutionary hypotheses for this genus. With the aim of filling these gaps in the evolutionary history for Dendroctonus, we analyzed ∼1800 loci mapped to a reference genome obtained for 20 of the 21 species recognized to date, minimizing the impact of missing information and improving the assumption of orthology in a phylogenomic framework. We obtained congruent phylogenetic topologies from two phylogenomic inference strategies: loci concatenation (ML framework) and a multispecies coalescent model (MSC) through the analysis of site pattern frequencies (SNPs). Dendroctonus is composed of two major clades (A and B), each containing five and four subclades, respectively. According to our divergence dating analysis, the MRCA for Dendroctonus dates back to the early Eocene, while the MRCA for each major clade diverged in the mid-Eocene. Interestingly, most of the speciation events of extant species occurred during the Miocene, which could be correlated with the diversification of pine trees (Pinus). The MRCA for Dendroctonus inhabited large regions of North America, with all ancestors and descendants of clade A having diversified within this region. The Mexican Transition Zone is important in the diversification processes for the majority of clade A species. For clade B, we identified two important colonization events to the Old World from America: the first in the early Oligocene from the Arctic to Asia (via Beringia), and the second during the Miocene from the Arctic-Western-Alleghany region to Europe and Siberia (also via Beringia). Our genomic analyses also supported the existence of hidden structured lineages within the frontalis complex, and also that D. beckeri represent a lineage independent from D. valens, as previously suggested. The information presented here updates the knowledge concerning the diversification of a genus with remarkable ecological and economic importance.

Chromosome-level genome assembly reveals genomic architecture of northern range expansion in the mountain pine beetle, Dendroctonus ponderosae Hopkins (Coleoptera: Curculionidae)

Genome sequencing methods and assembly tools have improved dramatically since the 2013 publication of draft genome assemblies for the mountain pine beetle, Dendroctonus ponderosae Hopkins (Coleoptera: Curculionidae). We conducted proximity ligation library sequencing and scaffolding to improve contiguity, and then used linkage mapping and recent bioinformatic tools for correction and further improvement. The new assemblies have dramatically improved contiguity and gaps compared to the originals: N50 values increased 26- to 36-fold, and the number of gaps were reduced by half. Ninety per cent of the content of the assemblies is now contained in 12 and 11 scaffolds for the female and male assemblies, respectively. Based on linkage mapping information, the 12 largest scaffolds in both assemblies represent all 11 autosomal chromosomes and the neo-X chromosome. These assemblies now have nearly chromosome-sized scaffolds and will be instrumental for studying genomic architecture, chromosome evolution, population genomics, functional genomics, and adaptation in this and other pest insects. We also identified regions in two chromosomes, including the ancestral-X portion of the neo-X chromosome, with elevated differentiation between northern and southern Canadian populations.

Champions of winter survival: cold acclimation and molecular regulation of cold hardiness in evergreen conifers

Evergreen conifers are champions of winter survival, based on their remarkable ability to acclimate to cold and develop cold hardiness. Counterintuitively, autumn cold acclimation is triggered not only by exposure to low temperature, but also by a combination of decreasing temperature, decreasing photoperiod and changes in light quality. These environmental cues control a network of signaling pathways that coordinate cold acclimation and cold hardiness in overwintering conifers, leading to cessation of growth, bud dormancy, freezing tolerance and changes in energy metabolism. Advances in genomic, transcriptomic and metabolomic tools for conifers have improved our understanding of how trees sense and respond to changes in temperature and light during cold acclimation and the development of cold hardiness, but there remain considerable gaps deserving further research in conifers. In the first section of this review, we focus on the physiological mechanisms used by evergreen conifers to adjust metabolism seasonally and to protect overwintering tissues against winter stresses. In the second section, we review how perception of low temperature and photoperiod regulate the induction of cold acclimation. Finally, we explore the evolutionary context of cold acclimation in conifers and evaluate challenges imposed on them by changing climate and discuss emerging areas of research in the field.

Spatiotemporal analyses suggest the role of glacial history and the ice-free corridor in shaping American badger population genetic variation

Recurring glacial cycles through the Quaternary period drastically altered the size and distribution of natural populations of North American flora and fauna. The "southerly refugia model" has been the longstanding framework for testing the effects of glaciation on contemporary genetic patterns; however, insights from ancient DNA have contributed to the reconstruction of more complex histories for some species. The American badger, Taxidea taxus, provides an interesting species for exploring the genetic legacy of glacial history, having been hypothesized to have postglacially emerged from a single, southerly refugium to recolonize northern latitudes. However, previous studies have lacked genetic sampling from areas where distinct glacial refugia have been hypothesized, including the Pacific Northwest and American Far North (Yukon, Alaska). In order to further investigate the phylogeographic history of American badgers, we collected mitochondrial DNA sequence data from ancient subfossil material collected within the historical range (Alaska, Yukon) and combined them with new and previously published data from across the species' contemporary distribution (n = 1,207). We reconstructed a mostly unresolved phylogenetic tree and star-like haplotype network indicative of emergence from a largely panmictic glacial refugium and recent population expansion, the latter further punctuated by significantly negative Tajima's D and Fu's Fs values. Although directionality of migration cannot be unequivocally inferred, the moderate to high levels of genetic variation exhibited by American badgers, alongside the low frequency of haplotypes with indels in the Midwest, suggest a potential recolonization into central North America after the hypothesized ice-free corridor reopened ~13,000 years ago. Overall, the expanded reconstruction of phylogeographic history of American badgers offers a broader understanding of contemporary range-wide patterns and identifies unique genetic units that can likely be used to inform conservation of at-risk populations at the northern periphery.

Spatial genetic structure of Rocky Mountain bighorn sheep (Ovis canadensis canadensis) at the northern limit of their native range

The Canadian Rocky Mountains are one of the few places on Earth where the spatial genetic structure of wide-ranging species has been relatively unaffected by anthropogenic disturbance. We characterized the spatial genetic structure of Rocky Mountain bighorn sheep (Ovis canadensis canadensis Shaw, 1804) in the northern portion of their range. Using microsatellites from 1495 individuals and mitochondrial DNA sequences from 188 individuals, we examined both broad- and fine-scale spatial genetic structure, assessed sex-biased gene flow within the northern portion of the species range, and identified geographic patterns of genetic diversity. We found that broad-scale spatial genetic structure was consistent with barriers to movement created by major river valleys. The fine-scale spatial genetic structure was characterized by a strong isolation-by-distance pattern, and analysis of neighborhood size using spatial autocorrelation indicated gene flow frequently occurred over distances of up to 100 km. However, analysis of sex-specific spatial autocorrelation and analysis of mitochondrial haplotype distributions failed to detect any evidence of sex-biased gene flow. Finally, our analyses reveal decreasing genetic diversity with increasing latitude, consistent with patterns of post-glacial recolonization of the Rocky Mountains.

Repurposing population genetics data to discern genomic architecture: A case study of linkage cohort detection in mountain pine beetle (Dendroctonus ponderosae)

Genetic surveys of the population structure of species can be used as resources for exploring their genomic architecture. By adjusting filtering assumptions, genome-wide single-nucleotide polymorphism (SNP) datasets can be reused to give new insights into the genetic basis of divergence and speciation without targeted resampling of specimens. Filtering only for missing data and minor allele frequency, we used a combination of principal components analysis and linkage disequilibrium network analysis to distinguish three cohorts of variable SNPs in the mountain pine beetle in western Canada, including one that was sex-linked and one that was geographically associated. These marker cohorts indicate genomically localized differentiation, and their detection demonstrates an accessible and intuitive method for discovering potential islands of genomic divergence without a priori knowledge of a species' genomic architecture. Thus, this method has utility for directly addressing the genomic architecture of species and generating new hypotheses for functional research.

New microsatellites revealed strong gene flow among populations of a new outbreak pest, Athetis lepigone (Möschler)

Athetis lepigone (Möschler) (Lepidoptera: Noctuidae) is a new outbreak pest in China. Consequently, it is unclear whether the emergence and spread of the outbreak of this pest are triggered by rapid in situ population size increases in each outbreak area, or by immigrants from a potential source area in China. In order to explore the outbreak process of this pest through a population genetics approach, we developed ten novel polymorphic expressed sequence tags (EST)-derived microsatellites. These new microsatellites had moderately high levels of polymorphism in the tested population. The number of alleles per locus ranged from 3 to 19, with an average of 8.6, and the expected heterozygosity ranged from 0.269 to 0.783. A preliminary population genetic analysis using these new microsatellites revealed a lack of population genetic structure in natural populations of A. lepigone. The estimates of recent migration rate revealed strong gene flow among populations. In conclusion, our study developed the first set of EST-microsatellite markers and shed a new light on the population genetic structure of this pest in China.

Selection of the sex-linked inhibitor of apoptosis in mountain pine beetle (Dendroctonus ponderosae) driven by enhanced expression during early overwintering

The mountain pine beetle (Dendroctonus ponderosae) is an insect native to western North America; however, its geographical range has recently expanded north in BC and east into Alberta. To understand the population structure in the areas of expansion, 16 gene-linked microsatellites were screened and compared to neutral microsatellites using outlier analyses of F_st and F_ct values. One sex-linked gene, inhibitor of apoptosis (IAP), showed a strong signature of positive selection for neo-X alleles and was analyzed for evidence of adaptive variation. Alleles of IAP were sequenced, and differences between the neo-X and neo-Y alleles were consistent with neutral evolution suggesting that the neo-Y allele may not be under functional constraints. Neo-Y alleles were amplified from gDNA, but not effectively from cDNA, suggesting that there was little IAP expression from neo-Y alleles. There were no differences in overall IAP expression between males and females with the common northern neo-X allele suggesting that the neo-X allele in males compensates for the reduced expression of neo-Y alleles. However, males lacking the most common northern neo-X allele thought to be selected for in northern populations had reduced overall IAP expression in early October-at a time when beetles are preparing for overwintering. This suggests that the most common allele may have more rapid upregulation. The reduced function of neo-Y alleles of IAP suggested by both sequence differences and lower levels of expression may foster a highly selective environment for neo-X alleles such as the common northern allele with more efficient upregulation.

Cold Tolerance of Mountain Pine Beetle (Coleoptera: Curculionidae) Eggs From the Historic and Expanded Ranges

Winter mortality is expected to be a key factor determining the ability of mountain pine beetle, Dendroctonus ponderosae Hopkins (Coleoptera: Curculionidae), to expand its range in Canada. We determined the mortality rate and supercooling points of eggs from the beetle's historic range in southern British Columbia as well as the recently expanded range in north-central Alberta and tested if eggs require an extended period of chilling to reach their maximum cold tolerance. We found no effect of population source or acclimation time on egg cold tolerance. Although 50% of eggs can survive brief exposure to -20.5 °C (LT50), storage at 0.3 °C and -7.5 °C for 59 d resulted in 50% and 100% mortality, respectively. Our results indicate that eggs suffer significant prefreeze mortality and are not well-adapted to overwintering: eggs are unlikely to survive winter throughout much of the beetle's range. Our results provide information that can be used to help model the climatic suitability of mountain pine beetle, including how changes in seasonality associated with new or changing climates may affect winter survival. In addition to lower lethal temperatures, it is critical that the duration of exposure to sublethal cold temperatures are considered in a comprehensive index of cold tolerance and incorporated into survival and population models.

Confirmation of independent introductions of an exotic plant pathogen of Cornus species, Discula destructiva, on the east and west coasts of North America

Cornus florida (flowering dogwood) and C. nuttallii (Pacific dogwood) are North American native tree species that belong to the big-bracted group of dogwoods. Cornus species are highly valued for their ornamental characteristics, and have fruits that contain high fat content for animals. Also, they are an important understory tree in natural forests. Dogwood anthracnose, caused by Discula destructiva, was observed in the late 1970s on the east and west coasts of the United States and by 1991 had quickly spread throughout most of the native ranges of C. florida and C. nuttalli. We investigated the genetic diversity and population structure of 93 D. destructiva isolates using 47 microsatellite loci developed from the sequenced genome of the type strain of D. destructiva. Clone-corrected data indicated low genetic diversity and the presence of four genetic clusters that corresponded to two major geographic areas, the eastern United States and the Pacific Northwest, and to the two collection time periods when the isolates were collected (pre- and post-1993). Linkage disequilibrium was present in five out of six subpopulations, suggesting that the fungus only reproduced asexually. Evidence of population bottlenecks was indicated across four identified genetic clusters, and was probably the result of the limited number of founding individuals on both coasts. These results support the hypothesis that D. destructiva is an exotic pathogen with independent introductions on the east and west coasts of North America. We also tested the cross-amplification of these microsatellite primers to other Discula species. Genomic DNA from 17 isolates of four other Discula species and two isolates of Juglanconis species (formerly Melanconis species) were amplified by 17 of 47 primer pairs. These primers may be useful for investigating the genetic diversity and population structure of these Discula species.

Expansion Under Climate Change: The Genetic Consequences

Range expansion and range shifts are crucial population responses to climate change. Genetic consequences are not well understood but are clearly coupled to ecological dynamics that, in turn, are driven by shifting climate conditions. We model a population with a deterministic reaction-diffusion model coupled to a heterogeneous environment that develops in time due to climate change. We decompose the resulting travelling wave solution into neutral genetic components to analyse the spatio-temporal dynamics of its genetic structure. Our analysis shows that range expansions and range shifts under slow climate change preserve genetic diversity. This is because slow climate change creates range boundaries that promote spatial mixing of genetic components. Mathematically, the mixing leads to so-called pushed travelling wave solutions. This mixing phenomenon is not seen in spatially homogeneous environments, where range expansion reduces genetic diversity through gene surfing arising from pulled travelling wave solutions. However, the preservation of diversity is diminished when climate change occurs too quickly. Using diversity indices, we show that fast expansions and range shifts erode genetic diversity more than slow range expansions and range shifts. Our study provides analytical insight into the dynamics of travelling wave solutions in heterogeneous environments.

Adaptive and neutral markers both show continent-wide population structure of mountain pine beetle (Dendroctonus ponderosae)

Assessments of population genetic structure and demographic history have traditionally been based on neutral markers while explicitly excluding adaptive markers. In this study, we compared the utility of putatively adaptive and neutral single‐nucleotide polymorphisms (SNPs) for inferring mountain pine beetle population structure across its geographic range. Both adaptive and neutral SNPs, and their combination, allowed range‐wide structure to be distinguished and delimited a population that has recently undergone range expansion across northern British Columbia and Alberta. Using an equal number of both adaptive and neutral SNPs revealed that adaptive SNPs resulted in a stronger correlation between sampled populations and inferred clustering. Our results suggest that adaptive SNPs should not be excluded prior to analysis from neutral SNPs as a combination of both marker sets resulted in better resolution of genetic differentiation between populations than either marker set alone. These results demonstrate the utility of adaptive loci for resolving population genetic structure in a nonmodel organism.

Modeling Landscape-Level Spatial Variation in Sex Ratio Skew in the Mountain Pine Beetle (Coleoptera: Curculionidae)

Through their influence on effective population sizes, sex ratio skew affects population dynamics. We examined spatial variation in female-biased sex ratios in the mountain pine beetle (Dendroctonus ponderosae Hopkins) outbreak in western Canada to better understand how environmental context affects sex ratio skew. Our specific objectives were to: 1) characterize spatial variation in mountain pine beetle sex ratio; 2) test previously asserted hypotheses that beetle sex ratio varies with tree diameter and year in outbreak; and 3) develop predictive models of sex ratio skew for larval and adult populations. Using logistic regression, we modeled the probability that an individual beetle (n = 2,369) was female as a function of multiple environmental variables across 34 stands in British Columbia and Alberta, Canada. We identified a consistent female-biased sex ratio with significantly greater skew in adults (2:1, n = 713) than in larvae (1.76:1, n = 1,643). We found that the proportion of larval females increased with decreasing tree size and with outbreak age. However, adults did not respond to tree size and larvae did not respond to outbreak age. Predictive models differed between larvae and adults. All identified models perform well and included predictors related to weather, tree diameter, and year in outbreak. Female-biased sex ratios appear to originate from differential male mortality during development rather than from sex-biased oviposition, suggesting sex ratio skew is not the cause of outbreaks, but rather a consequence.

Polygamy and an absence of fine-scale structure in Dendroctonus ponderosae (Hopk.) (Coleoptera: Curcilionidae) confirmed using molecular markers

An understanding of mating systems and fine-scale spatial genetic structure is required to effectively manage forest pest species such as Dendroctonus ponderosae (mountain pine beetle). Here we used genome-wide single-nucleotide polymorphisms to assess the fine-scale genetic structure and mating system of D. ponderosae collected from a single stand in Alberta, Canada. Fine-scale spatial genetic structure was absent within the stand and the majority of genetic variation was best explained at the individual level. Relatedness estimates support previous reports of pre-emergence mating. Parentage assignment tests indicate that a polygamous mating system better explains the relationships among individuals within a gallery than the previously reported female monogamous/male polygynous system. Furthermore, there is some evidence to suggest that females may exploit the galleries of other females, at least under epidemic conditions. Our results suggest that current management models are likely to be effective across large geographic areas based on the absence of fine-scale genetic structure.

The genetic signature of rapid range expansions: How dispersal, growth and invasion speed impact heterozygosity and allele surfing

As researchers collect spatiotemporal population and genetic data in tandem, models that connect demography and dispersal to genetics are increasingly relevant. The dominant spatiotemporal model of invasion genetics is the stepping-stone model which represents a gradual range expansion in which individuals jump to uncolonized locations one step at a time. However, many range expansions occur quickly as individuals disperse far from currently colonized regions. For these types of expansion, stepping-stone models are inappropriate. To more accurately reflect wider dispersal in many organisms, we created kernel-based models of invasion genetics based on integrodifference equations. Classic theory relating to integrodifference equations suggests that the speed of range expansions is a function of population growth and dispersal. In our simulations, populations that expanded at the same speed but with spread rates driven by dispersal retained more heterozygosity along axes of expansion than range expansions with rates of spread that were driven primarily by population growth. To investigate surfing we introduced mutant alleles in wave fronts of simulated range expansions. In our models based on random mating, surfing alleles remained at relatively low frequencies and surfed less often compared to previous results based on stepping-stone simulations with asexual reproduction.

Population structure of mountain pine beetle symbiont Leptographium longiclavatum and the implication on the multipartite beetle-fungi relationships

Over 18 million ha of forests have been destroyed in the past decade in Canada by the mountain pine beetle (MPB) and its fungal symbionts. Understanding their population dynamics is critical to improving modeling of beetle epidemics and providing potential clues to predict population expansion. Leptographium longiclavatum and Grosmannia clavigera are fungal symbionts of MPB that aid the beetle to colonize and kill their pine hosts. We investigated the genetic structure and demographic expansion of L. longiclavatum in populations established within the historic distribution range and in the newly colonized regions. We identified three genetic clusters/populations that coincide with independent geographic locations. The genetic profiles of the recently established populations in northern British Columbia (BC) and Alberta suggest that they originated from central and southern BC. Approximate Bayesian Computation supports the scenario that this recent expansion represents an admixture of individuals originating from BC and the Rocky Mountains. Highly significant correlations were found among genetic distance matrices of L. longiclavatum, G. clavigera, and MPB. This highlights the concordance of demographic processes in these interacting organisms sharing a highly specialized niche and supports the hypothesis of long-term multipartite beetle-fungus co-evolutionary history and mutualistic relationships.

How the mountain pine beetle (Dendroctonus ponderosae) breached the Canadian Rocky Mountains

The mountain pine beetle (MPB; Dendroctonus ponderosae Hopkins), a major pine forest pest native to western North America, has extended its range north and eastward during an ongoing outbreak. Determining how the MPB has expanded its range to breach putative barriers, whether physical (nonforested prairie and high elevation of the Rocky Mountains) or climatic (extreme continental climate where temperatures can be below −40 °C), may contribute to our general understanding of range changes as well as management of the current epidemic. Here, we use a panel of 1,536 single nucleotide polymorphisms (SNPs) to assess population genetic structure, connectivity, and signals of selection within this MPB range expansion. Biallelic SNPs in MPB from southwestern Canada revealed higher genetic differentiation and lower genetic connectivity than in the northern part of its range. A total of 208 unique SNPs were identified using different outlier detection tests, of which 32 returned annotations for products with putative functions in cholesterol synthesis, actin filament contraction, and membrane transport. We suggest that MPB has been able to spread beyond its previous range by adjusting its cellular and metabolic functions, with genome scale differentiation enabling populations to better withstand cooler climates and facilitate longer dispersal distances. Our study is the first to assess landscape-wide selective adaptation in an insect. We have shown that interrogation of genomic resources can identify shifts in genetic diversity and putative adaptive signals in this forest pest species.

Effects of introgression on the genetic population structure of two ecologically and economically important conifer species: lodgepole pine (Pinus contorta var. latifolia) and jack pine (Pinus banksiana)

Forest trees exhibit a remarkable range of adaptations to their environment, but as a result of frequent and long-distance gene flow, populations are often only weakly differentiated. Lodgepole and jack pine hybridize in western Canada, which adds the opportunity for introgression through hybridization to contribute to population structure and (or) adaptive variation. Access to large sample size, high density SNP datasets for these species would improve our ability to resolve population structure, parameterize introgression, and separate the influence of demography from adaptation. To accomplish this, 454 transcriptome reads for lodgepole and jack pine were assembled using Newbler and MIRA, the assemblies mined for SNPs, and 1536 SNPs were selected for typing on lodgepole pine, jack pine, and their hybrids (N = 536). We identified population structure using both Bayesian clustering and discriminate analysis of principle components. Introgressed SNP loci were identified and their influence on observed population structure was assessed. We found that introgressed loci resulted in increased differentiation both within lodgepole and jack pine populations. These findings are timely given the recent mountain pine beetle population expansion in the hybrid zone, and will facilitate future studies of adaptive traits in these ecologically important species.

Draft genome of the mountain pine beetle, Dendroctonus ponderosae Hopkins, a major forest pest

BACKGROUND

The mountain pine beetle, Dendroctonus ponderosae Hopkins, is the most serious insect pest of western North American pine forests. A recent outbreak destroyed more than 15 million hectares of pine forests, with major environmental effects on forest health, and economic effects on the forest industry. The outbreak has in part been driven by climate change, and will contribute to increased carbon emissions through decaying forests.

RESULTS

We developed a genome sequence resource for the mountain pine beetle to better understand the unique aspects of this insect's biology. A draft de novo genome sequence was assembled from paired-end, short-read sequences from an individual field-collected male pupa, and scaffolded using mate-paired, short-read genomic sequences from pooled field-collected pupae, paired-end short-insert whole-transcriptome shotgun sequencing reads of mRNA from adult beetle tissues, and paired-end Sanger EST sequences from various life stages. We describe the cytochrome P450, glutathione S-transferase, and plant cell wall-degrading enzyme gene families important to the survival of the mountain pine beetle in its harsh and nutrient-poor host environment, and examine genome-wide single-nucleotide polymorphism variation. A horizontally transferred bacterial sucrose-6-phosphate hydrolase was evident in the genome, and its tissue-specific transcription suggests a functional role for this beetle.

CONCLUSIONS

Despite Coleoptera being the largest insect order with over 400,000 described species, including many agricultural and forest pest species, this is only the second genome sequence reported in Coleoptera, and will provide an important resource for the Curculionoidea and other insects.

Temperature-driven range expansion of an irruptive insect heightened by weakly coevolved plant defenses

Warming climate has increased access of native bark beetles to high-elevation pines that historically received only intermittent exposure to these tree-killing herbivores. Here we show that a dominant, relatively naïve, high-elevation species, whitebark pine, has inferior defenses against mountain pine beetle compared with its historical lower-elevation host, lodgepole pine. Lodgepole pines respond by exuding more resin and accumulating higher concentrations of toxic monoterpenes than whitebark pine, where they co-occur. Furthermore, the chemical composition of whitebark pine appears less able to inhibit the pheromonal communication beetles use to jointly overcome tree defenses. Despite whitebark pine's inferior defenses, beetles were more likely to attack their historical host in mixed stands. This finding suggests there has been insufficient sustained contact for beetles to alter their complex behavioral mechanisms driving host preference. In no-choice assays, however, beetles readily entered and tunneled in both hosts equally, and in stands containing less lodgepole pine, attacks on whitebark pines increased. High-elevation trees in pure stands may thus be particularly vulnerable to temperature-driven range expansions. Predators and competitors were more attracted to volatiles from herbivores attacking their historical host, further increasing risk in less coevolved systems. Our results suggest cold temperatures provided a sufficient barrier against herbivores for high-elevation trees to allocate resources to other physiological processes besides defense. Changing climate may reduce the viability of that evolutionary strategy, and the life histories of high-elevation trees seem unlikely to foster rapid counter adaptation. Consequences extend from reduced food supplies for endangered grizzly bears to altered landscape and hydrological processes.