Genetic rescue often leads to higher fitness as a result of increased heterozygosity across animal taxa.

Biodiversity loss has reached critical levels partly due to anthropogenic habitat loss and degradation. These landscape changes are damaging as they can fragment species distributions into small, isolated populations, resulting in limited gene flow, population declines and reduced adaptive potential. Genetic rescue, the translocation of individuals to increase genetic diversity and ultimately fitness, has produced promising results for fragmented populations but remains underutilized due to a lack of long-term data and monitoring. To promote a better understanding of genetic rescue and its potential risks and benefits over the short-term, we reviewed and analysed published genetic rescue attempts to identify whether genetic diversity increases following translocation, and if this change is associated with increased fitness. Our review identified 19 studies that provided genetic and fitness data from before and after the translocation; the majority of these were on mammals, and included experimental, natural and conservation-motivated translocations. Using a Bayesian meta-analytical approach, we found that on average, genetic diversity and fitness increased in populations post translocations, although there were some exceptions to this trend. Overall, genetic diversity was a positive predictor of population fitness, and in some cases this relationship extended three generations post-rescue. These data suggest a single translocation can have lasting fitness benefits, and support translocation as another tool to facilitate conservation success. Given the limited number of studies with long-term data, we echo the need for genetic monitoring of populations post-translocation to understand whether genetic rescue can also limit the loss of adaptive potential in the long-term.

A cost-effective RNA extraction and RT-qPCR approach to detect California serogroup viruses from pooled mosquito samples.

Mosquito-borne diseases pose ongoing global health concerns, demanding more cost-efficient methods to detect pathogens to support enhanced surveillance efforts. This study introduces an adapted TRIzol-based high-throughput RNA extraction protocol, tailored for the detection of California serogroup viruses in pooled mosquito samples in a rapid and cost-effective manner. This approach provided consistent RNA yields and sensitive viral detection relative to two commercial extraction kits (QIAGEN RNeasy Mini Kit and MACHEREY-NAGEL NucleoSpin RNA Plus Kit). The incorporation of a user-friendly and non-spiking-based RT-qPCR internal control designed for the 18S rRNA gene in mosquitoes minimizes potential false positives/negatives, improving the fidelity of viral detection outcomes. Effective RNA yields, purity, and successful target amplification across 25 mosquito species and varied pool sizes (1-50 mosquitoes per tube) affirm the reliability of our approach. The extraction method is cost-effective, with an incurred cost of $0.58 CAD per sample, in contrast to the $5.25 CAD cost per sample of the two kits, rendering it promising for mosquito-borne disease surveillance initiatives.

Relationships between water quality and mosquito presence and abundance: a systematic review and meta-analysis.

Mosquito-borne diseases (MBDs) are emerging in response to climate and land use changes. As mosquito (Diptera: Culicidae) habitat selection is often contingent on water availability for egg and larval development, studies have recognized water quality also influences larval habitats. However, underlying species-, genera-, and mosquito level preferences for water quality conditions are varied. This systematic review and meta-analysis aimed to identify, characterize, appraise, and synthesize available global data on the relationships between water quality and mosquito presence and abundance (MPA); with the goal to further our understanding of the geographic expansion of MBD risks. A systematic review was conducted to identify studies investigating the relationships between water quality properties and MPA. Where appropriate, random-effects meta-analyses were conducted to provide pooled estimates for the association between the most reported water quality properties and MPA. The most reported water quality parameters were pH (87%), nitrogen concentrations (56%), turbidity (56%), electrical conductivity (54%), dissolved oxygen (43%), phosphorus concentrations (30%), and alkalinity (10%). Overall, pH (P = 0.05), turbidity (P < 0.0001), electrical conductivity (P = 0.005), dissolved oxygen (P < 0.0001), nitrogen (P < 0.0001), and phosphorus (P < 0.0001) showed significantly positive pooled correlations with MPA, while alkalinity showed a nonsignificant null pooled correlation (P = 0.85). We observed high heterogeneity in most meta-analyses, and climate zonation was shown to influence the pooled estimates. Linkages between MPA and water quality properties will enhance our capacity to predict MBD risks under changing environmental and land use changes.

A roadmap to robust discriminant analysis of principal components.

Identification of population structure is a common goal for a variety of applications, including conservation, wildlife management, and medical genetics. The outcome of these analyses can have far reaching implications; therefore, it is important to ensure an understanding of the strengths and weaknesses of the methodologies used. Increasing in popularity, the discriminant analysis of principal components (DAPC) method incorporates combinations of genetic variables (alleles) into a model that differentiates individuals into genetic clusters. However, users may not have a full understanding of how to best parameterize the model. In this issue of Thia (Molecular Ecology Resources, 2022) looks under the hood of the DAPC. Using simulated data, he demonstrates the importance of careful parameter selection in developing a DAPC model, what the implications are for over-fitting the model, and finally, how best to evaluate the results of DAPC models. This work highlights the issues that can arise when over-parameterizing the DAPC model when gene flow is high among clusters and provides important guidelines to ensure researchers are making conclusions that are biologically relevant.

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.

Extent and direction of introgressive hybridization of mule and white-tailed deer in western Canada.

Hybridization of mule deer (Odocoileus hemionus) and white-tailed deer (O. virginianus) appears to be a semi-regular occurrence in western North America. Previous studies confirmed the presence of hybrids in a variety of sympatric habitats, but their developing molecular resources limited identification to the earliest, most admixed generations. For this reason, estimates of hybrid production in wild populations often rely on anecdotal reports. As well, white-tailed deer populations' continued encroachment into historically mule deer-occupied habitats due to changes in land use, habitat homogenization, and a warming climate may increase opportunities for interspecific encounters. We sought to quantify the prevalence and extent of hybrid deer in the prairies of western Canada using a SNP assay with enhanced discriminating power. By updating the available molecular resources, we sought to identify and characterize previously cryptic introgression. We also investigated the influence of various parameters on hybridity by way of logistic regression. We observed overall hybridization rates of ~1.0%, slightly lower than that reported by previous studies, and found white-tailed-like hybrids to be more common than their mule deer-like counterparts. Here, we build upon past studies of hybridization in North American deer by increasing hybrid detection power, expanding sample sizes, demonstrating a new molecular resource applicable to future research and observing asymmetrical directionality of introgression.

Genetic association with boldness and maternal performance in a free-ranging population of grey seals (Halichoerus grypus).

Individual variation in quantitative traits clearly influence many ecological and evolutionary processes. Moderate to high heritability estimates of personality and life-history traits suggest some level of genetic control over these traits. Yet, we know very little of the underlying genetic architecture of phenotypic variation in the wild. In this study, we used a candidate gene approach to investigate the association of genetic variants with repeated measures of boldness and maternal performance traits (weaning mass and lactation duration) collected over an 11- and 28-year period, respectively, in a free-ranging population of grey seals on Sable Island National Park Reserve, Canada. We isolated and re-sequenced five genes: dopamine receptor D4 (DRD4), serotonin transporter (SERT), oxytocin receptor (OXTR), and melanocortin receptors 1 (MC1R) and 5 (MC5R). We discovered single nucleotide polymorphisms (SNPs) in each gene; and, after accounting for loci in linkage disequilibrium and filtering due to missing data, we were able to test for genotype-phenotype relationships at seven loci in three genes (DRD4, SERT, and MC1R). We tested for association between these loci and traits of 180 females having extreme shy-bold phenotypes using mixed-effects models. One locus within SERT was significantly associated with boldness (effect size = 0.189) and a second locus within DRD4 with weaning mass (effect size = 0.232). Altogether, genotypes explained 6.52-13.66% of total trait variation. Our study substantiates SERT and DRD4 as important determinants of behaviour, and provides unique insight into the molecular mechanisms underlying maternal performance variation in a marine predator.

Large-scale prion protein genotyping in Canadian caribou populations and potential impact on chronic wasting disease susceptibility.

Polymorphisms within the prion protein gene (Prnp) are an intrinsic factor that can modulate chronic wasting disease (CWD) pathogenesis in cervids. Although wild European reindeer (Rangifer tarandus tarandus) were infected with CWD, as yet there have been no reports of the disease in North American caribou (R. tarandus spp.). Previous Prnp genotyping studies on approximately 200 caribou revealed single nucleotide polymorphisms (SNPs) at codons 2 (V/M), 129 (G/S), 138 (S/N), 146 (N/n) and 169 (V/M). The impact of these polymorphisms on CWD transmission is mostly unknown, except for codon 138. Reindeer carrying at least one allele encoding for asparagine (138NN or 138SN) are less susceptible to clinical CWD upon infection by natural routes, with the majority of prions limited to extraneural tissues. We sequenced the Prnp coding region of two caribou subspecies (n = 986) from British Columbia, Saskatchewan, Yukon, Nunavut and the Northwest Territories, to identify SNPs and their frequencies. Genotype frequencies at codon 138 differed significantly between barren-ground (R. t. groenlandicus) and woodland (R. t. caribou) caribou when we excluded the Chinchaga herd (p < .05). We also found new variants at codons 153 (Y/F) and 242 (P/L). Our findings show that the 138N allele is rare among caribou in areas with higher risk of contact with CWD-infected species. As both subspecies are classified as Threatened and play significant roles in North American Indigenous culture, history, food security and the economy, determining frequencies of Prnp genotypes associated with susceptibility to CWD is important for future wildlife management measures.

The influence of a priori grouping on inference of genetic clusters: simulation study and literature review of the DAPC method.

Inference of genetic clusters is a key aim of population genetics, sparking development of numerous analytical methods. Within these, there is a conceptual divide between finding de novo structure versus assessment of a priori groups. Recently developed, Discriminant Analysis of Principal Components (DAPC), combines discriminant analysis (DA) with principal component (PC) analysis. When applying DAPC, the groups used in the DA (specified a priori or described de novo) need to be carefully assessed. While DAPC has rapidly become a core technique, the sensitivity of the method to misspecification of groups and how it is being empirically applied, are unknown. To address this, we conducted a simulation study examining the influence of a priori versus de novo group designations, and a literature review of how DAPC is being applied. We found that with a priori groupings, distance between genetic clusters reflected underlying FST. However, when migration rates were high and groups were described de novo there was considerable inaccuracy, both in terms of the number of genetic clusters suggested and placement of individuals into those clusters. Nearly all (90.1%) of 224 studies surveyed used DAPC to find de novo clusters, and for the majority (62.5%) the stated goal matched the results. However, most studies (52.3%) omit key run parameters, preventing repeatability and transparency. Therefore, we present recommendations for standard reporting of parameters used in DAPC analyses. The influence of groupings in genetic clustering is not unique to DAPC, and researchers need to consider their goal and which methods will be most appropriate.

Long-Term Incubation PrPCWD with Soils Affects Prion Recovery but Not Infectivity.

Chronic wasting disease (CWD) is a contagious prion disease of cervids. The infectious agent is shed from animals at the preclinical and clinical stages of disease where it persists in the environment as a reservoir of CWD infectivity. In this study, we demonstrate that long-term incubation of CWD prions (generated from tg-mice infected with deer or elk prions) with illite, montmorillonite (Mte) and whole soils results in decreased recovery of PrPCWD, suggesting that binding becomes more avid and irreversible with time. This continual decline of immunoblot PrPCWD detection did not correlate with prion infectivity levels. Bioassay showed no significant differences in incubation periods between mice inoculated with 1% CWD brain homogenate (BH) and with the CWD-BH pre-incubated with quartz or Luvisolic Ae horizon for 1 or 30 weeks. After 55 weeks incubation with Chernozem and Luvisol, bound PrPCWD was not detectable by immunoblotting but remained infectious. This study shows that although recovery of PrPCWD bound to soil minerals and whole soils with time become more difficult, prion infectivity is not significantly altered. Detection of prions in soil is, therefore, not only affected by soil type but also by length of time of the prion-soil interaction.

Predicting the spread-risk potential of chronic wasting disease to sympatric ungulate species.

Wildlife disease incidence is increasing, resulting in negative impacts on the economy, biodiversity, and potentially human health. Chronic wasting disease (CWD) is a fatal, transmissible spongiform encephalopathy of cervids (wild and captive) which continues to spread geographically resulting in exposure to potential new host species. The disease agent (PrPCWD) is a misfolded conformer of the cellular prion protein (PrPC). In Canada, the disease is endemic in Alberta and Saskatchewan, affecting mule and white-tail deer, with lesser impact on elk and moose. As the disease continues to expand, additional wild ungulate species including bison, bighorn sheep, mountain goat, and pronghorn antelope may be exposed. To better understand the species-barrier, we reviewed the current literature on taxa naturally or experimentally exposed to CWD to identify susceptible and resistant species. We created a phylogeny of these taxa using cytochrome B and found that CWD susceptibility followed the species phylogeny. Using this phylogeny we estimated the probability of CWD susceptibility for wild ungulate species. We then compared PrPC amino acid polymorphisms among these species to identify which sites segregated between susceptible and resistant species. We identified sites that were significantly associated with susceptibility, but they were not fully discriminating. Finally, we sequenced Prnp from 578 wild ungulates to further evaluate their potential susceptibility. Together, these data suggest the host-range for CWD will potentially include pronghorn, mountain goat and bighorn sheep, but bison are likely to be more resistant. These findings highlight the need for monitoring potentially susceptible species as CWD continues to expand.

Confidently identifying the correct K value using the ΔK method: When does K = 2?

Populations delineated based on genetic data are commonly used for wildlife conservation and management. Many studies use the program structure combined with the ΔK method to identify the most probable number of populations (K). We recently found K = 2 was identified more often when studies used ΔK compared to studies that did not. We suggested two reasons for this: hierarchical population structure leads to underestimation, or the ΔK method does not evaluate K = 1 causing an overestimation. The present contribution aims to develop a better understanding of the limits of the method using one, two and three population simulations across migration scenarios. From these simulations we identified the "best K" using model likelihood and ΔK. Our findings show that mean probability plots and ΔK are unable to resolve the correct number of populations once migration rate exceeds 0.005. We also found a strong bias towards selecting K = 2 using the ΔK method. We used these data to identify the range of values where the ΔK statistic identifies a value of K that is not well supported. Finally, using the simulations and a review of empirical data, we found that the magnitude of ΔK corresponds to the level of divergence between populations. Based on our findings, we suggest researchers should use the ΔK method cautiously; they need to report all relevant data, including the magnitude of ΔK, and an estimate of connectivity for the research community to assess whether meaningful genetic structure exists within the context of management and conservation.

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.

Development of a Novel Mule Deer Genomic Assembly and Species-Diagnostic SNP Panel for Assessing Introgression in Mule Deer, White-Tailed Deer, and Their Interspecific Hybrids.

Mule deer (Odocoileus hemionus) are endemic to a wide variety of habitats in western North America, many of which are shared in sympatry with their closely related sister-species white-tailed deer (Odocoileus virginianus), whom they hybridize with in wild populations. Although mule deer meet many ideal conditions for a molecular ecological research species, such as high abundance, ecological importance, and broad dispersal and gene flow, conservation genetic studies have been limited by a relative lack of existing genomic resources and inherent difficulties caused by introgression with white-tailed deer. Many molecular tools currently available for the study of cervids were designed using reference assemblies of divergent model species, specifically cattle (Bos taurus). Bovidae and Cervidae diverged approximately 28 million years ago, therefore, we sought to ameliorate the available resources by contributing the first mule deer whole genome sequence draft assembly with an average genome-wide read depth of 25X, using the white-tailed genome assembly (Ovir.te_1.0) as a reference. Comparing the two assemblies, we identified ∼33 million single nucleotide polymorphisms (SNPs) and insertion/deletion variants. We then verified fixed SNP differences between the two species and developed a 40-loci SNP assay capable of identifying pure mule deer, white-tailed deer, and interspecific hybrids. Assignment capacity of the panel, which was tested on simulated datasets, is reliable up to and including the third backcross hybrid generation. Identification of post-F1 hybrids will be necessary for hybrid zone population studies going forward, and the new mule deer assembly will be a valuable resource for genetic and comparative genomics studies.

Soil humic acids degrade CWD prions and reduce infectivity.

Chronic wasting disease (CWD), an environmentally transmissible, fatal prion disease is endemic in North America, present in South Korea and has recently been confirmed in northern Europe. The expanding geographic range of this contagious disease of free-ranging deer, moose, elk and reindeer has resulted in increasing levels of prion infectivity in the environment. Soils are involved in CWD horizontal transmission, acting as an environmental reservoir, and soil mineral and organic compounds have the ability to bind prions. Upper horizons of soils are usually enriched with soil organic matter (SOM), however, the role of SOM in prion conservation and mobility remains unclear. In this study, we show that incubation of PrPCWD with humic acids (HA), a major SOM compound, affects both the molecular weight and recovery of PrPCWD. Detection of PrPCWD is reduced as HA concentration increases. Native HA extracted from pristine soils also reduces or entirely eliminates PrPCWD signal. Incubation of CWD prions with HA significantly increased incubation periods in tgElk mice demonstrating that HA can reduce CWD infectivity.

The K = 2 conundrum.

Assessments of population genetic structure have become an increasing focus as they can provide valuable insight into patterns of migration and gene flow. structure, the most highly cited of several clustering-based methods, was developed to provide robust estimates without the need for populations to be determined a priori. structure introduces the problem of selecting the optimal number of clusters, and as a result, the ΔK method was proposed to assist in the identification of the "true" number of clusters. In our review of 1,264 studies using structure to explore population subdivision, studies that used ΔK were more likely to identify K = 2 (54%, 443/822) than studies that did not use ΔK (21%, 82/386). A troubling finding was that very few studies performed the hierarchical analysis recommended by the authors of both ΔK and structure to fully explore population subdivision. Furthermore, extensions of earlier simulations indicate that, with a representative number of markers, ΔK frequently identifies K = 2 as the top level of hierarchical structure, even when more subpopulations are present. This review suggests that many studies may have been over- or underestimating population genetic structure; both scenarios have serious consequences, particularly with respect to conservation and management. We recommend publication standards for population structure results so that readers can assess the implications of the results given their own understanding of the species biology.

Mule deer spatial association patterns and potential implications for transmission of an epizootic disease.

Animal social behaviour can have important effects on the long-term dynamics of diseases. In particular, preferential spatial relationships between individuals can lead to differences in the rates of disease spread within a population. We examined the concurrent influence of genetic relatedness, sex, age, home range overlap, time of year, and prion disease status on proximal associations of adult Rocky Mountain mule deer (Odocoileus hemionus hemionus) in a chronic wasting disease endemic area. We also quantified the temporal stability of these associations across different sex, age, and disease status classes. We used three years of high frequency telemetry data from 74 individuals to record encounters within 25 m of each other, and to calculate seasonal home range overlap measured by volume of intersection (VI). The strength of pairwise spatial association between adult mule deer was independent of genetic relatedness, age and disease status. Seasonal variation in association strength was not consistent across years, perhaps due to annual changes in weather conditions. The influence of home range overlap on association strength varied seasonally, whereby associations were stronger in pre-rut and fawning than in the rest of the seasons. The sexes of individuals also interacted with both VI and season. At increasing levels of VI, associations were stronger between females than between males and between females and males. The strongest associations in pre-rut were between males, while the strongest in rut were between females and males. The temporal stability of associations was markedly dependant on the sex and the diagnosis of the associating pair. Our findings highlight the importance of considering concurrent effects of biological and environmental factors when seeking to understand the role of social preference in behavioural ecology and disease spread. Applying this knowledge in epidemiological modelling will shed light on the dynamics of disease transmission among mule deer.

Circumpolar Genetic Structure and Recent Gene Flow of Polar Bears: A Reanalysis.

Recently, an extensive study of 2,748 polar bears (Ursus maritimus) from across their circumpolar range was published in PLOS ONE, which used microsatellites and mitochondrial haplotypes to apparently show altered population structure and a dramatic change in directional gene flow towards the Canadian Archipelago-an area believed to be a future refugium for polar bears as their southernmost habitats decline under climate change. Although this study represents a major international collaborative effort and promised to be a baseline for future genetics work, methodological shortcomings and errors of interpretation undermine some of the study's main conclusions. Here, we present a reanalysis of this data in which we address some of these issues, including: (1) highly unbalanced sample sizes and large amounts of systematically missing data; (2) incorrect calculation of FST and of significance levels; (3) misleading estimates of recent gene flow resulting from non-convergence of the program BayesAss. In contrast to the original findings, in our reanalysis we find six genetic clusters of polar bears worldwide: the Hudson Bay Complex, the Western and Eastern Canadian Arctic Archipelago, the Western and Eastern Polar Basin, and-importantly-we reconfirm the presence of a unique and possibly endangered cluster of bears in Norwegian Bay near Canada's expected last sea-ice refugium. Although polar bears' abundance, distribution, and population structure will certainly be negatively affected by ongoing-and increasingly rapid-loss of Arctic sea ice, these genetic data provide no evidence of strong directional gene flow in response to recent climate change.

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.