Evidence that Ophiostomatoid Fungal Symbionts of Mountain Pine Beetle Do Not Play a Role in Overcoming Lodgepole Pine Defenses During Mass Attack.

Mountain pine beetle (MPB; Dendroctonus ponderosae Hopkins) is a devastating forest insect pest that has killed millions of hectares of pines in western North America over the past two decades. Like other bark beetles, MPB vectors ophiostomatoid fungal species, some of which are pathogenic to host pine species. The phytopathogenicity of these fungal symbionts has sparked considerable debate regarding their role in facilitating MPB attack success. We tested the hypothesis that MPB ophiostomatoid fungal associates like Grosmannia clavigera (Robinson-Jeffrey and Davidson) Zipfel, de Beer and Wingfield contribute to overwhelming host defenses during MPB mass attack. We compared responses of mature lodgepole pine (Pinus contorta Dougl. ex Loud. var. latifolia Engelm.) trees growing in natural stands that were mass attacked by MPB with those inoculated with G. clavigera by examining host defense hormones, secondary metabolites, and gene expression profiles. The jasmonate and ethylene signatures of necrotrophic pathogen-triggered response were identified in G. clavigera-inoculated trees, but only the jasmonate signature of a herbivore-triggered response was measured in MPB-attacked trees. Several G. clavigera-induced changes in pine phenolic metabolite profiles and phenolic biosynthesis gene expression patterns were absent in MPB-attacked pines. These findings indicate that ophiostomatoid fungi like G. clavigera are not a major factor in overwhelming host defenses during MPB mass attack. Instead, fungal pathogenicity likely is more important in aiding MPB colonization and development within the host tree. Phenolics appear to play a larger role in the host response to G. clavigera than to MPB, although phenolics may also influence MPB feeding and behavior. [Formula: see text] Copyright © 2024 The Author(s). This is an open access article distributed under the CC BY 4.0 International license.

Evidence of Coevolution Between Cronartium harknessii Lineages and Their Corresponding Hosts, Lodgepole Pine and Jack Pine.

Variation in rate of infection and susceptibility of Pinus spp. to the fungus Cronartium harknessii (syn. Endocronartium harknessii), the causative agent of western gall rust, has been well documented. To test the hypothesis that there is a coevolutionary relationship between C. harknessii and its hosts, we examined genetic structure and virulence of C. harknessii associated with lodgepole pine (P. contorta var. latifolia), jack pine (P. banksiana), and their hybrids. A secondary objective was to improve assessment and diagnosis of infection in hosts. Using 18 microsatellites, we assessed genetic structure of C. harknessii from 90 sites within the ranges of lodgepole pine and jack pine. We identified two lineages (East and West, FST = 0.677) associated with host genetic structure (r = 0.81, P = 0.001), with East comprising three sublineages. In parallel, we conducted a factorial experiment in which lodgepole pine, jack pine, and hybrid seedlings were inoculated with spores from the two primary genetic lineages. With this experiment, we refined the phenotypic categories associated with infection and demonstrated that stem width can be used as a quantitative measure of host response to infection. Overall, each host responded differentially to the fungal lineages, with jack pine exhibiting more resiliency to infection than lodgepole pine and hybrids exhibiting intermediate resiliency. Taken together, the shared genetic structure between fungus and host species, and the differential interaction of the fungal species with the hosts, supports a coevolutionary relationship between host and pathogen.[Formula: see text] Copyright © 2022 The Author(s). This is an open access article distributed under the CC BY-NC-ND 4.0 International license.

Linking genotype to phenotype to identify genetic variation relating to host susceptibility in the mountain pine beetle system.

Identifying genetic variants responsible for phenotypic variation under selective pressure has the potential to enable productive gains in natural resource conservation and management. Despite this potential, identifying adaptive candidate loci is not trivial, and linking genotype to phenotype is a major challenge in contemporary genetics. Many of the population genetic approaches commonly used to identify adaptive candidates will simultaneously detect false positives, particularly in nonmodel species, where experimental evidence is seldom provided for putative roles of the adaptive candidates identified by outlier approaches. In this study, we use outcomes from population genetics, phenotype association, and gene expression analyses as multiple lines of evidence to validate candidate genes. Using lodgepole and jack pine as our nonmodel study species, we analyzed 17 adaptive candidate loci together with 78 putatively neutral loci at 58 locations across Canada (N > 800) to determine whether relationships could be established between these candidate loci and phenotype related to mountain pine beetle susceptibility. We identified two candidate loci that were significant across all population genetic tests, and demonstrated significant changes in transcript abundance in trees subjected to wounding or inoculation with the mountain pine beetle fungal associate Grosmannia clavigera. Both candidates are involved in central physiological processes that are likely to be invoked in a trees response to stress. One of these two candidate loci showed a significant association with mountain pine beetle attack status in lodgepole pine. The spatial distribution of the attack-associated allele further coincides with other indicators of susceptibility in lodgepole pine. These analyses, in which population genetics was combined with laboratory and field experimental validation approaches, represent first steps toward linking genetic variation to the phenotype of mountain pine beetle susceptibility in lodgepole and jack pine, and provide a roadmap for more comprehensive analyses.

A novel application of RNase H2-dependent quantitative PCR for detection and quantification of Grosmannia clavigera, a mountain pine beetle fungal symbiont, in environmental samples.

Mountain pine beetle (Dendroctonus ponderosae Hopkins; MPB) is an economically and ecologically important pest of pine species in western North America. Mountain pine beetles form complex multipartite relationships with microbial partners, including the ophiostomoid fungi Grosmannia clavigera (Robinson-Jeffrey and Davidson) Zipfel, de Beer and Wingfield, Ophiostoma montium (Rumbold) von Arx, Grosmannia aurea (Robinson-Jeffrey and Davidson) Zipfel, de Beer and Wingfield, Leptographium longiclavatum (Lee, Kim, and Breuil) and Leptographium terebrantis (Barras and Perry). These fungi are vectored by MPB to new pine hosts, where the fungi overcome host defenses to grow into the sapwood. A tree's relative susceptibility to these fungi is conventionally assessed by measuring lesions that develop in response to fungal inoculation. However, these lesions represent a symptom of infection, representing both fungal growth and tree defense capacity. In order to more objectively assess fungal virulence and host tree susceptibility in studies of host-pathogen interactions, a reliable, consistent, sensitive method is required to accurately identify and quantify MPB-associated fungal symbionts in planta. We have adapted RNase H2-dependent PCR, a technique originally designed for rare allele discrimination, to develop a novel RNase H2-dependent quantitative PCR (rh-qPCR) assay that shows greater specificity and sensitivity than previously published PCR-based methods to quantify MPB fungal symbionts in pine xylem and MPB whole beetles. Two sets of assay probes were designed: one that amplifies a broad range of ophiostomoid species, and a second that amplifies G. clavigera but not other MPB-associated ophiostomoid species. Using these primers to quantify G. clavigera in pine stems, we provide evidence that lesion length does not accurately reflect the extent of fungal colonization along the stem nor the quantity of fungal growth within this colonized portion of stem. The sensitivity, specificity, reproducibility, cost effectiveness and high-throughput potential of the rh-qPCR assay makes the technology suitable for identification and quantification of a wide array of pathogenic and beneficial microbes that form associations with plants and other organisms, even when the microbial partner is present in low abundance.

Genetic and genomic evidence of niche partitioning and adaptive radiation in mountain pine beetle fungal symbionts.

Bark beetles form multipartite symbiotic associations with blue stain fungi (Ophiostomatales, Ascomycota). These fungal symbionts play an important role during the beetle's life cycle by providing nutritional supplementation, overcoming tree defences and modifying host tissues to favour brood development. The maintenance of stable multipartite symbioses with seemingly less competitive symbionts in similar habitats is of fundamental interest to ecology and evolution. We tested the hypothesis that the coexistence of three fungal species associated with the mountain pine beetle is the result of niche partitioning and adaptive radiation using SNP genotyping coupled with genotype-environment association analysis and phenotypic characterization of growth rate under different temperatures. We found that genetic variation and population structure within each species is best explained by distinct spatial and environmental variables. We observed both common (temperature seasonality and the host species) and distinct (drought, cold stress, precipitation) environmental and spatial factors that shaped the genomes of these fungi resulting in contrasting outcomes. Phenotypic intraspecific variations in Grosmannia clavigera and Leptographium longiclavatum, together with high heritability, suggest potential for adaptive selection in these species. By contrast, Ophiostoma montium displayed narrower intraspecific variation but greater tolerance to extreme high temperatures. Our study highlights unique phenotypic and genotypic characteristics in these symbionts that are consistent with our hypothesis. By maintaining this multipartite relationship, the bark beetles have a greater likelihood of obtaining the benefits afforded by the fungi and reduce the risk of being left aposymbiotic. Complementarity among species could facilitate colonization of new habitats and survival under adverse conditions.

Elucidation of the polyamine biosynthesis pathway during Brazilian pine (Araucaria angustifolia) seed development.

Polyamines (PAs), such as spermidine and spermine, as well as amino acids that are substrates for their biosynthesis, are known to be essential for plant development. However, little is known about the gene expression and metabolic switches associated with the ornithine/arginine and PA biosynthetic pathway during seed development in conifers. To understand these metabolic switches, the enzyme activity of arginine decarboxylase and ornithine decarboxylase, as well as the contents of PAs and amino acids were evaluated in three Araucaria angustifolia (Bertol. Kuntze) seed developmental stages in combination with expression profile analyses of genes associated with the ornithine/arginine and PA biosynthetic pathway. Twelve genes were selected for further analysis and it was shown that the expression profiles of AaADC and AaSAMDC were up-regulated during zygotic embryo development. Polyamines and amino acids were found to accumulate differently in embryos and megagametophytes, and the transition from the globular to the cotyledonary stage was marked by an increase in free and conjugated spermidine and spermine contents. Putrescine is made from arginine, which was present at low content at the late embryogenesis stage, when high content of citrulline was observed. Differences in amino acids, PAs and gene expression profiles of biosynthetic genes at specific seed stages and at each seed transition stage were investigated, providing insights into molecular and physiological aspects of conifer embryogenesis for use in future both basic and applied studies.

Poplar CBF1 functions specifically in an integrated cold regulatory network.

The C-repeat binding factors (CBFs), also termed dehydration-responsive element-binding protein 1 (DREB1) family members, play crucial roles in the acquisition of stress tolerance, but in trees, the underlying mechanisms of stress tolerance remain elusive. To gain insight into these mechanisms, we isolated five CBF1 orthologs from four poplar sections (Populus spp.) and assessed their expression under drought, cold, heat and salt stress conditions. Globally induced expression in response to cold suggested a correlation between poplar CBF1 expression and the acquisition of cold tolerance. Responses that varied between sections may reflect section-specific stress tolerance mechanisms, suggesting an effect of ecological context on the development of CBF1-mediated stress tolerance in poplar. We then used a genome-wide search strategy in Populus trichocarpa to predict 2263 putative CBF target genes; the identified genes participate in multiple biological processes and pathways. Almost all of the putative target genes contained multiple cis-acting elements that mediate responses to various environmental and endogenous signals, consistent with an important role of CBF1s in an integrated cold regulatory network. Finally, analysis of an association population of 528 individuals of Populus simonii identified six single-nucleotide polymorphisms (false discovery rate Q < 0.10) significantly (P < 0.005) associated with malondialdehyde production and electrolyte leakage, suggesting the potential importance of PsCBF1 in the regulation of some membrane-related functions. Our findings provide new insights into the function of PsCBF1 and shed light on the CBF-mediated regulatory network in poplar.

Advances in ecological genomics in forest trees and applications to genetic resources conservation and breeding.

Forest trees are an unparalleled group of organisms in their combined ecological, economic and societal importance. With widespread distributions, predominantly random mating systems and large population sizes, most tree species harbour extensive genetic variation both within and among populations. At the same time, demographic processes associated with Pleistocene climate oscillations and land-use change have affected contemporary range-wide diversity and may impinge on the potential for future adaptation. Understanding how these adaptive and neutral processes have shaped the genomes of trees species is therefore central to their management and conservation. As for many other taxa, the advent of high-throughput sequencing methods is expected to yield an understanding of the interplay between the genome and environment at a level of detail and depth not possible only a few years ago. An international conference entitled 'Genomics and Forest Tree Genetics' was held in May 2016, in Arcachon (France), and brought together forest geneticists with a wide range of research interests to disseminate recent efforts that leverage contemporary genomic tools to probe the population, quantitative and evolutionary genomics of trees. An important goal of the conference was to discuss how such data can be applied to both genome-enabled breeding and the conservation of forest genetic resources under land use and climate change. Here, we report discoveries presented at the meeting and discuss how the ecological genomic toolkit can be used to address both basic and applied questions in tree biology.

Change in soil fungal community structure driven by a decline in ectomycorrhizal fungi following a mountain pine beetle (Dendroctonus ponderosae) outbreak.

Western North American landscapes are rapidly being transformed by forest die-off caused by mountain pine beetle (Dendroctonus ponderosae), with implications for plant and soil communities. The mechanisms that drive changes in soil community structure, particularly for the highly prevalent ectomycorrhizal fungi in pine forests, are complex and intertwined. Critical to enhancing understanding will be disentangling the relative importance of host tree mortality from changes in soil chemistry following tree death. Here, we used a recent bark beetle outbreak in lodgepole pine (Pinus contorta) forests of western Canada to test whether the effects of tree mortality altered the richness and composition of belowground fungal communities, including ectomycorrhizal and saprotrophic fungi. We also determined the effects of environmental factors (i.e. soil nutrients, moisture, and phenolics) and geographical distance, both of which can influence the richness and composition of soil fungi. The richness of both groups of soil fungi declined and the overall composition was altered by beetle-induced tree mortality. Soil nutrients, soil phenolics and geographical distance influenced the community structure of soil fungi; however, the relative importance of these factors differed between ectomycorrhizal and saprotrophic fungi. The independent effects of tree mortality, soil phenolics and geographical distance influenced the community composition of ectomycorrhizal fungi, while the community composition of saprotrophic fungi was weakly but significantly correlated with the geographical distance of plots. Taken together, our results indicate that both deterministic and stochastic processes structure soil fungal communities following landscape-scale insect outbreaks and reflect the independent roles tree mortality, soil chemistry and geographical distance play in regulating the community composition of soil fungi.

The joint influence of photoperiod and temperature during growth cessation and development of dormancy in white spruce (Picea glauca).

Timely responses to environmental cues enable the synchronization of phenological life-history transitions essential for the health and survival of north-temperate and boreal tree species. While photoperiodic cues will remain persistent under climate change, temperature cues may vary, contributing to possible asynchrony in signals influencing developmental and physiological transitions essential to forest health. Understanding the relative contribution of photoperiod and temperature as determinants of the transition from active growth to dormancy is important for informing adaptive forest management decisions that consider future climates. Using a combination of photoperiod (long = 20 h or short = 8 h day lengths) and temperature (warm = 22 °C/16 °C and cool = 8 °C/4 °C day/night, respectively) treatments, we used microscopy, physiology and modeling to comprehensively examine hallmark traits of the growth-dormancy transition-including bud formation, growth cessation, cold hardiness and gas exchange-within two provenances of white spruce [Picea glauca (Moench) Voss] spanning a broad latitude in Alberta, Canada. Following exposure to experimental treatments, seedlings were transferred to favorable conditions, and the depth of dormancy was assessed by determining the timing and ability of spruce seedlings to resume growth. Short photoperiods promoted bud development and growth cessation, whereas longer photoperiods extended the growing season through the induction of lammas growth. In contrast, cool temperatures under both photoperiodic conditions delayed bud development. Photoperiod strongly predicted the development of cold hardiness, whereas temperature predicted photosynthetic rates associated with active growth. White spruce was capable of attaining endodormancy, but its release was environmentally determined. Dormancy depth varied substantially across experimental treatments suggesting that environmental cues experienced within one season could affect growth in the following season, which is particularly important for a determinate species such as white spruce. The joint influence of these environmental cues points toward the importance of including local constant photoperiod and shifting temperature cues into predictive models that consider how climate change may affect northern forests.

Responses of hybrid aspen over-expressing a PIP2;5 aquaporin to low root temperature.

Aquaporins mediate the movement of water across cell membranes. Plasma membrane intrinsic protein 2;5 from Populus trichocarpa×deltoides (PtdPIP2;5) was previously demonstrated to be a functionally important water conducting aquaporin. To study the relevance of aquaporin-mediated root water transport at low temperatures, we generated transgenic Populus tremula×alba over-expressing PtdPIP2;5 under control of the maize ubiquitin promoter, and compared the physiological responses and water transport properties of the PtdPIP2;5 over-expressing lines (PtdPIP2;5ox) with wild-type plants. We hypothesized that over-expression of PtdPIP2;5 would reduce temperature sensitivity of root water transport and gas exchange. Decreasing root temperatures to 10 and 5°C significantly decreased hydraulic conductivities (Lp) in wild-type plants, but had no significant effect on Lp in PtdPIP2;5ox plants. Recovery of Lp in the transgenic lines returned to 20°C from 5°C was faster than in the wild-type plants. Low root temperature did not induce major changes in transcript levels for other PIPs. When roots were exposed to 5°C in solution culture and shoots were exposed to 20°C, wild-type plants had significantly lower net photosynthetic and transpiration rates compared to PtdPIP2;5ox plants. Taken together, our results demonstrate that over-expression of PtdPIP2;5 in P. tremula×alba was effective in alleviating the effects of low root temperature on Lp and gas exchange.

Hydraulic conductivity and aquaporin transcription in roots of trembling aspen (Populus tremuloides) seedlings colonized by Laccaria bicolor.

Ectomycorrhizal fungi have been reported to increase root hydraulic conductivity (L pr) by altering apoplastic and plasma membrane intrinsic protein (PIP)-mediated cell-to-cell water transport pathways in associated roots, or to have little effect on root water transport, depending on the interacting species and imposed stresses. In this study, we investigated the water transport properties and PIP transcription in roots of aspen (Populus tremuloides) seedlings colonized by the wild-type strain of Laccaria bicolor and by strains overexpressing a major fungal water-transporting aquaporin JQ585595. Inoculation of aspen seedlings with L. bicolor resulted in about 30 % colonization rate of root tips, which developed dense mantle and the Hartig net that was restricted in the modified root epidermis. Transcript abundance of the aspen aquaporins PIP1;2, PIP2;1, and PIP2;2 decreased in colonized root tips. Root colonization by JQ585595-overexpressing strains had no significant impact on seedling shoot water potentials, gas exchange, or dry mass; however, it led to further decrease in transcript abundance of PIP1;2 and PIP2;3 and the significantly lower L pr than in non-inoculated roots. These results, taken together with our previous study that showed enhanced root water hydraulics of L. bicolor-colonized white spruce (Picea glauca), suggest that the impact of L. bicolor on root hydraulics varies by the ectomycorrhiza-associated tree species.

Differences in defence responses of Pinus contorta and Pinus banksiana to the mountain pine beetle fungal associate Grosmannia clavigera are affected by water deficit.

We tested the hypotheses that responses to the mountain pine beetle fungal associate Grosmannia clavigera will differ between the evolutionarily co-evolved host lodgepole pine (Pinus contorta var. latifolia) and the naïve host jack pine (Pinus banksiana) and that these responses will be influenced by water availability. G. clavigera inoculation resulted in more rapid stem lesion development in lodgepole than in jack pine; water deficit delayed lesion development in both species. Decreased hydraulic conductivity was observed in inoculated lodgepole pine seedlings, likely because of tracheid occlusion by fungal hyphae and/or metabolite accumulation. Drought but not inoculation significantly impacted bark abscisic acid levels. Jasmonic and salicylic acid were implicated in local and systemic responses of both species to G. clavigera, with salicylic acid appearing to play a greater role in jack pine response to G. clavigera than lodgepole pine. Water deficit increased constitutive levels and/or attenuated induced responses to G. clavigera for several monoterpenes in lodgepole but not jack pine. Instead, inoculation of well-watered but not water deficit jack pine resulted in a greater number of xylem resin ducts. These findings reveal mechanisms underlying differences in G. clavigera-induced responses between lodgepole and jack pine hosts, and how water availability modulates these responses.

Transcript profiling of aquaporins during basidiocarp development in Laccaria bicolor ectomycorrhizal with Picea glauca.

Sporocarp formation is part of the reproductive stage in the life cycle of many mycorrhizal macrofungi. Sporocarp formation is accompanied by a transcriptomic switch and profound changes in regulation of the gene families that play crucial roles in the sporocarp initiation and maturation. Since sporocarp growth requires efficient water delivery, in the present study, we investigated changes in transcript abundance of six fungal aquaporin genes that could be cloned from the ectomycorrhizal fungus Laccaria bicolor strain UAMH8232, during the initiation and development of its basidiocarp. Aquaporins are intrinsic membrane proteins facilitating the transmembrane transport of water and other small neutral molecules. In controlled-environment experiments, we induced basidiocarp formation in L. bicolor, which formed ectomycorrhizal associations with white spruce (Picea glauca) seedlings. We profiled transcript abundance corresponding to six fungal aquaporin genes at six different developmental stages of basidiocarp growth and development. We also compared physiological parameters of non-inoculated to mycorrhizal seedlings with and without the presence of basidiocarps. Two L. bicolor aquaporins--JQ585592, a functional channel for CO2, NO and H2O2, and JQ585595, a functional water channel--showed the greatest degree of upregulation during development of the basidiocarp. Our findings point to the importance of aquaporin-mediated transmembrane water and CO2 transport during distinct stages of basidiocarp development.

Overexpression of Laccaria bicolor aquaporin JQ585595 alters root water transport properties in ectomycorrhizal white spruce (Picea glauca) seedlings.

The contribution of hyphae to water transport in ectomycorrhizal (ECM) white spruce (Picea glauca) seedlings was examined by altering expression of a major water-transporting aquaporin in Laccaria bicolor. Picea glauca was inoculated with wild-type (WT), mock transgenic or L. bicolor aquaporin JQ585595-overexpressing (OE) strains and exposed to root temperatures ranging from 5 to 20°C to examine the root water transport properties, physiological responses and plasma membrane intrinsic protein (PIP) expression in colonized plants. Mycorrhization increased shoot water potential, transpiration, net photosynthetic rates, root hydraulic conductivity and root cortical cell hydraulic conductivity in seedlings. At 20°C, OE plants had higher root hydraulic conductivity compared with WT plants and the increases were accompanied by higher expression of P. glauca PIP GQ03401_M18.1 in roots. In contrast to WT L. bicolor, the effects of OE fungi on root and root cortical cell hydraulic conductivities were abolished at 10 and 5°C in the absence of major changes in the examined transcript levels of P. glauca root PIPs. The results provide evidence for the importance of fungal aquaporins in root water transport of mycorrhizal plants. They also demonstrate links between hyphal water transport, root aquaporin expression and root water transport in ECM plants.

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.

Decline of ectomycorrhizal fungi following a mountain pine beetle epidemic.

Forest die-off caused by mountain pine beetle (MPB; Dendroctonus ponderosa) is rapidly transforming western North American landscapes. The rapid and widespread death of lodgepole pine (Pinus contorta) will likely have cascading effects on biodiversity. One group particularly prone to such declines associated with MPB are ectomycorrhizal fungi, symbiotic organisms that can depend on pine for their survival, and are critical for stand regeneration. We evaluated the indirect effects of MPB on above- (community composition of epigeous sporocarps) and belowground (hyphal abundance) occurrences of ectomycorrhizal fungi across 11 forest stands. Along a gradient of mortality (0-82% pine killed), macromycete community composition changed; this shift was driven by a decrease in the species richness of ectomycorrhizal fungi. Both the proportion of species that were ectomycorrhizal and hyphal length in the soil declined with increased MPB-caused pine mortality; < 10% of sporocarp species were ectomycorrhizal in stands with high pine mortality compared with > 70% in stands without MPB attacks. The rapid range expansion of a native insect results not only in the widespread mortality of an ecologically and economically important pine species, but the effect of MPB may also be exacerbated by the concomitant decline of fungi crucial for recovery of these forests.

Cross-species outlier detection reveals different evolutionary pressures between sister species.

Lodgepole pine (Pinus contorta var. latifolia) and jack pine (Pinus banksiana) hybridize in western Canada, an area of recent mountain pine beetle range expansion. Given the heterogeneity of the environment, and indications of local adaptation, there are many unknowns regarding the response of these forests to future outbreaks. To better understand this we aim to identify genetic regions that have adaptive potential. We used data collected on 472 single nucleotide polymorphism (SNP) loci from 576 tree samples collected across 13 lodgepole pine-dominated sites and four jack pine-dominated sites. We looked at the relationship of genetic diversity with the environment, and we identified candidate loci using both frequency-based (arlequin and bayescan) and correlation-based (matsam and bayenv) methods. We found contrasting relationships between environmental variation and genetic diversity for the species. While we identified a number of candidate outliers (34 in lodgepole pine, 25 in jack pine, and 43 interspecific loci), we did not find any loci in common between lodgepole and jack pine. Many of the outlier loci identified were correlated with environmental variation. Using rigorous criteria we have been able to identify potential outlier SNPs. We have also found evidence of contrasting environmental adaptations between lodgepole and jack pine which could have implications for beetle spread risk.