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.

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.

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.

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.

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.

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.

Spatial patterns of neutral and functional genetic variations reveal patterns of local adaptation in raccoon (Procyon lotor) populations exposed to raccoon rabies.

Local adaptation is necessary for population survival and depends on the interplay between responses to selective forces and demographic processes that introduce or retain adaptive and maladaptive attributes. Host-parasite systems are dynamic, varying in space and time, where both host and parasites must adapt to their ever-changing environment in order to survive. We investigated patterns of local adaptation in raccoon populations with varying temporal exposure to the raccoon rabies virus (RRV). RRV infects approximately 85% of the population when epizootic and has been presumed to be completely lethal once contracted; however, disease challenge experiments and varying spatial patterns of RRV spread suggest some level of immunity may exist. We first assessed patterns of local adaptation in raccoon populations along the eastern seaboard of North America by contrasting spatial patterns of neutral (microsatellite loci) and functional, major histocompatibility complex (MHC) genetic diversity and structure. We explored variation of MHC allele frequencies in the light of temporal population exposure to RRV (0-60 years) and specific RRV strains in infected raccoons. Our results revealed high levels of MHC variation (66 DRB exon 2 alleles) and pronounced genetic structure relative to neutral microsatellite loci, indicative of local adaptation. We found a positive association linking MHC genetic diversity and temporal RRV exposure, but no association with susceptibility and resistance to RRV strains. These results have implications for landscape epidemiology studies seeking to predict the spread of RRV and present an example of how population demographics influence the degree to which populations adapt to local selective pressures.

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 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.

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.

Temporal analysis of genetic structure to assess population dynamics of reintroduced swift foxes.

Reintroductions are increasingly used to reestablish species, but a paucity of long-term postrelease monitoring has limited understanding of whether and when viable populations subsequently persist. We conducted temporal genetic analyses of reintroduced populations of swift foxes (Vulpes velox) in Canada (Alberta and Saskatchewan) and the United States (Montana). We used samples collected 4 years apart, 17 years from the initiation of the reintroduction, and 3 years after the conclusion of releases. To assess program success, we genotyped 304 hair samples, subsampled from the known range in 2000 and 2001, and 2005 and 2006, at 7 microsatellite loci. We compared diversity, effective population size, and genetic connectivity over time in each population. Diversity remained stable over time and there was evidence of increasing effective population size. We determined population structure in both periods after correcting for differences in sample sizes. The geographic distribution of these populations roughly corresponded with the original release locations, which suggests the release sites had residual effects on the population structure. However, given that both reintroduction sites had similar source populations, habitat fragmentation, due to cropland, may be associated with the population structure we found. Although our results indicate growing, stable populations, future connectivity analyses are warranted to ensure both populations are not subject to negative small-population effects. Our results demonstrate the importance of multiple sampling years to fully capture population dynamics of reintroduced populations. Análisis Temporal de la Estructura Genética para Evaluar la Dinámica Poblacional de Zorros (Vulpes velox) Reintroducidos.

Mountain pine beetle host-range expansion threatens the boreal forest.

The current epidemic of the mountain pine beetle (MPB), an indigenous pest of western North American pine, has resulted in significant losses of lodgepole pine. The leading edge has reached Alberta where forest composition shifts from lodgepole to jack pine through a hybrid zone. The susceptibility of jack pine to MPB is a major concern, but there has been no evidence of host-range expansion, in part due to the difficulty in distinguishing the parentals and their hybrids. We tested the utility of a panel of microsatellite loci optimized for both species to classify lodgepole pine, jack pine and their hybrids using simulated data. We were able to accurately classify simulated individuals, and hence applied these markers to identify the ancestry of attacked trees. Here we show for the first time successful MPB attack in natural jack pine stands at the leading edge of the epidemic. This once unsuitable habitat is now a novel environment for MPB to exploit, a potential risk which could be exacerbated by further climate change. The consequences of host-range expansion for the vast boreal ecosystem could be significant.

Of glaciers and refugia: a decade of study sheds new light on the phylogeography of northwestern North America.

Glacial cycles have played a dominant role in shaping the genetic structure and distribution of biota in northwestern North America. The two major ice age refugia of Beringia and the Pacific Northwest were connected by major mountain chains and bordered by the Pacific Ocean. As a result, numerous refugial options were available for the regions taxa during glacial advances. We reviewed the importance of glaciations and refugia in shaping northwestern North America's phylogeographic history. We also tested whether ecological variables were associated with refugial history. The recurrent phylogeographic patterns that emerged were the following: (i) additional complexity, i.e. refugia within refugia, in both Beringia and the Pacific Northwest; and (ii) strong evidence for cryptic refugia in the Alexander Archipelago and Haida Gwaii, the Canadian Arctic and within the ice-sheets. Species with contemporary ranges that covered multiple refugia, or those with high dispersal ability, were significantly more likely to have resided in multiple refugia. Most of the shared phylogeographic patterns can be attributed to multiple refugial locales during the last glacial maximum or major physiographic barriers like rivers and glaciers. However, some of the observed patterns are much older and appear connected to the orogeny of the Cascade-Sierra chain or allopatric differentiation during historic glacial advances. The emergent patterns from this review suggest we should refine the classic Beringian-southern refugial paradigm for northwestern North American biota and highlight the ecological and evolutionary consequences of colonization from multiple refugia.

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.

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.

Differential permeability of rivers to raccoon gene flow corresponds to rabies incidence in Ontario, Canada.

The correlation of landscape features with genetic discontinuities reveals barriers to dispersal that can contribute to understanding present and future spread of wildlife diseases. This knowledge can then be used for targeting control efforts. The impact of natural barriers on raccoon dispersal was assessed through genetic analysis of samples from two regions, Niagara (N = 666) and St. Lawrence (N = 802). These areas are transected by major rivers and are at the northern front of a raccoon rabies epizootic. Genetic clusters were identified in each region using Bayesian clustering algorithms. In the Niagara region, two clusters were identified corresponding to either side of the Niagara River. For the St. Lawrence region, spatially congruent clusters were not identified, despite the presence of the intervening St. Lawrence River. These genetic data are consistent with raccoon rabies incidence data where rabies has been detected across the St. Lawrence River in Ontario while no cases have been detected in Ontario across the Niagara River. This is despite expectations of rabies incidence in Niagara before the St. Lawrence based on the progression of rabies from New York. The results from the two regions suggest different permeabilities to raccoons between New York and Ontario that may be attributed to the rivers. However, other factors have also been explored that could contribute to this difference between these study sites including the shape of the landscape and resource distribution.

Assessing a landscape barrier using genetic simulation modelling: implications for raccoon rabies management.

Landscape barriers influence movement patterns of animals, which in turn, affect spatio-temporal spread of infectious wildlife disease. We compare genetic data from computer simulations to those acquired from field samples to measure the effect of a landscape barrier on raccoon (Procyon lotor) movement, enabling risk assessment of raccoon rabies disease spread across the Niagara River from New York State into Ontario, an area currently uninfected by rabies. An individual-based spatially explicit model is used to simulate the expansion of a raccoon population to cross the Niagara River, for different permeabilities of the river to raccoon crossings. Since the model records individual raccoon genetics, the genetic population structure of neutral mitochondrial DNA haplotypes are characterised in the expanding population, every 25 years, using a genetic distance measure, phi ST, Mantel tests and a gene diversity measure. The river barrier effect is assessed by comparing genetic measures computed from model outputs to those calculated from 166 raccoons recently sampled from the same landscape. The "best fit" between modelled scenarios and field data indicate the river prevents 50% of attempts to cross the river. Founder effects dominated the colonizing genetic population structure, and, as the river barrier effect increased, its genetic diversity decreased. Using gene flow to calibrate the effect of the river as a barrier to movement provides an estimate of the effect of a river in reducing the likelihood of cross-river infection. Including individual genetic markers in simulation modelling benefits investigations of disease spread and control.

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.