Adaptive latitudinal cline of photoperiodic diapause induction in the parasitoid Nasonia vitripennis in Europe

Living in seasonally changing environments requires adaptation to seasonal cycles. Many insects use the change in day length as a reliable cue for upcoming winter and respond to shortened photoperiod through diapause. In this study, we report the clinal variation in photoperiodic diapause induction in populations of the parasitoid wasp Nasonia vitripennis collected along a latitudinal gradient in Europe. In this species, diapause occurs in the larval stage and is maternally induced. Adult Nasonia females were exposed to different photoperiodic cycles and lifetime production of diapausing offspring was scored. Females switched to the production of diapausing offspring after exposure to a threshold number of photoperiodic cycles. A latitudinal cline was found in the proportion of diapausing offspring, the switch point for diapause induction measured as the maternal age at which the female starts to produce diapausing larvae, and the critical photoperiod for diapause induction. Populations at northern latitudes show an earlier switch point, higher proportions of diapausing individuals and longer critical photoperiods. Since the photoperiodic response was measured under the same laboratory conditions, the observed differences between populations most likely reflect genetic differences in sensitivity to photoperiodic cues, resulting from local adaptation to environmental cycles. The observed variability in diapause response combined with the availability of genomic tools for N. vitripennis represent a good opportunity to further investigate the genetic basis of this adaptive trait.     

Local adaptation along an environmental cline in a species with an inversion polymorphism

Polymorphic inversions are ubiquitous across the animal kingdom and are frequently associated with clines in inversion frequencies across environmental gradients. Such clines are thought to result from selection favouring local adaptation; however, empirical tests are scarce. The seaweed fly Coelopa frigida has an α/β inversion polymorphism, and previous work demonstrated that the α inversion frequency declines from the North Sea to the Baltic Sea and is correlated with changes in tidal range, salinity, algal composition and wrackbed stability. Here, we explicitly test the hypothesis that populations of C. frigida along this cline are locally adapted by conducting a reciprocal transplant experiment of four populations along this cline to quantify survival. We found that survival varied significantly across treatments and detected a significant Location x Substrate interaction, indicating local adaptation. Survival models showed that flies from locations at both extremes had highest survival on their native substrates, demonstrating that local adaptation is present at the extremes of the cline. Survival at the two intermediate locations was, however, not elevated at the native substrates, suggesting that gene flow in intermediate habitats may override selection. Together, our results support the notion that population extremes of species with polymorphic inversions are often locally adapted, even when spatially close, consistent with the growing view that inversions can have direct and strong effects on the fitness of species.     

Footprints of selection in wild populations of Bicyclus anynana along a latitudinal cline

One of the major questions in ecology and evolutionary biology is how variation in the genome enables species to adapt to divergent environments. Here, we study footprints of thermal selection in candidate genes in six wild populations of the afrotropical butterfly Bicyclus anynana sampled along a c. 3000 km latitudinal cline. We sequenced coding regions of 31 selected genes with known functions in metabolism, pigment production, development and heat shock responses. These include genes for which we expect a priori a role in thermal adaptation and, thus, varying selection pressures along a latitudinal cline, and genes we do not expect to vary clinally and can be used as controls. We identified amino acid substitution polymorphisms in 13 genes and tested these for clinal variation by correlation analysis of allele frequencies with latitude. In addition, we used two F_ST‐based outlier methods to identify loci with higher population differentiation than expected under neutral evolution, while accounting for potentially confounding effects of population structure and demographic history. Two metabolic enzymes of the glycolytic pathway, UGP and Treh, showed clinal variation. The same loci showed elevated population differentiation and were identified as significant outliers. We found no evidence of clines in the pigmentation genes, heat shock proteins and developmental genes. However, we identified outlier loci in more localized parts of the range in the pigmentation genes yellow and black. We discuss that the observed clinal variation and elevated population divergence in UGP and Treh may reflect adaptation to a geographic thermal gradient.     

The influence of immune activation on thermal tolerance along a latitudinal cline

Global change is shifting both temperature patterns and the geographic distribution of pathogens, and infection has already been shown to substantially reduce host thermal performance, potentially placing populations at greater risk that previously thought. But what about individuals that are able to successfully clear an infection? Whilst the direct damage a pathogen causes will likely lead to reductions in host's thermal tolerance, the response to infection often shares many underlying pathways with the general stress response, potentially acting as a buffer against subsequent thermal stress. Here, by exposing Drosophila melanogaster to heat‐killed bacterial pathogens, we investigate how activation of a host's immune system can modify any response to both heat and cold temperature stress. In a single focal population, we find that immune activation can improve a host's knockdown times during heat shock, potentially offsetting some of the damage that would subsequently arise as an infection progresses. Conversely, immune activation had a detrimental effect on CT_max and did not influence lower thermal tolerance as measured by chill‐coma recovery time. However, we also find that the influence of immune activation on heat knockdown times is not generalizable across an entire cline of locally adapted populations. Instead, immune activation led to signals of local adaptation to temperature being lost, erasing the previous advantage that populations in warmer regions had when challenged with heat stress. Our results suggest that activation of the immune system may help buffer individuals against the detrimental impact of infection on thermal tolerance; however, any response will be population specific and potentially not easily predicted across larger geographic scales, and dependent on the form of thermal stress faced by a host.     

Physiological and regulatory underpinnings of geographic variation in reptilian cold tolerance across a latitudinal cline

Understanding the mechanisms that produce variation in thermal performance is a key component to investigating climatic effects on evolution and adaptation. However, disentangling the effects of local adaptation and phenotypic plasticity in shaping patterns of geographic variation in natural populations can prove challenging. Additionally, the physiological mechanisms that cause organismal dysfunction at extreme temperatures are still largely under debate. Using the green anole, Anolis carolinensis, we integrate measures of cold tolerance (CT_min), standard metabolic rate, heart size, blood lactate concentration and RNAseq data from liver tissue to investigate geographic variation in cold tolerance and its underlying mechanisms along a latitudinal cline. We found significant effects of thermal acclimation and latitude of origin on variation in cold tolerance. Increased cold tolerance correlates with decreased rates of oxygen consumption and blood lactate concentration (a proxy for oxygen limitation), suggesting elevated performance is associated with improved oxygen economy during cold exposure. Consistent with these results, co‐expression modules associated with blood lactate concentration are enriched for functions associated with blood circulation, coagulation and clotting. Expression of these modules correlates with thermal acclimation and latitude of origin. Our findings support the oxygen and capacity‐limited thermal tolerance hypothesis as a potential contributor to variation in reptilian cold tolerance. Moreover, differences in gene expression suggest regulation of the blood coagulation cascade may play an important role in reptilian cold tolerance and may be the target of natural selection in populations inhabiting colder environments.     

Bud phenology and growth are subject to divergent selection across a latitudinal gradient in Populus angustifolia and impact adaptation across the distributional range and associated arthropods

Temperate forest tree species that span large geographical areas and climatic gradients often have high levels of genetic variation. Such species are ideal for testing how neutral demographic factors and climate‐driven selection structure genetic variation within species, and how this genetic variation can affect ecological communities. Here, we quantified genetic variation in vegetative phenology and growth traits in narrowleaf cottonwood, Populus angustifolia, using three common gardens planted with genotypes originating from source populations spanning the species' range along the Rocky Mountains of North America (ca. 1700 km). We present three main findings. First, we found strong evidence of divergent selection (Q_ST > F_ST) on fall phenology (bud set) with adaptive consequences for frost avoidance. We also found evidence for selection on bud flush duration, tree height, and basal diameter, resulting in population differentiation. Second, we found strong associations with climate variables that were strongly correlated with latitude of origin. More strongly differentiated traits also showed stronger climate correlations, which emphasizes the role that climate has played in divergent selection throughout the range. We found population × garden interaction effects; for some traits, this accounted for more of the variance than either factor alone. Tree height was influenced by the difference in climate of the source and garden locations and declined with increasing transfer distance. Third, growth traits were correlated with dependent arthropod community diversity metrics. Synthesis. Overall, we conclude that climate has influenced genetic variation and structure in phenology and growth traits and leads to local adaptation in P. angustifolia, which can then impact dependent arthropod species. Importantly, relocation of genotypes far northward or southward often resulted in poor growth, likely due to a phenological mismatch with photoperiod, the proximate cue for fall growth cessation. Genotypes moved too far southward suffer from early growth cessation, whereas those moved too far northward are prone to fall frost and winter dieback. In the face of current and forecasted climate change, habitat restoration, forestry, and tree breeding efforts should utilize these findings to better match latitudinal and climatic source environments with management locations for optimal future outcomes.     

Polymorphism in the neurofibromin gene,Nf1, is associated with antagonistic selection on wing size and development time in Drosophila melanogaster

In many invertebrates, body size shows genetically based clines, with size increasing in colder climates. Large body size is typically associated with prolonged development times. We consider variation in the CNS‐specific gene neurofibromin 1 (Nf1) and its association with body size and development time. We identified two major Nf1 haplotypes in natural populations, Nf1‐insertion‐A and Nf1‐deletion‐G. These haplotypes are characterized by a 45‐base insertion/deletion (INDEL) in Nf1 intron 2 and an A/G synonymous substitution (locus L17277). Linkage disequilibrium (LD) between the INDEL and adjacent sites is high but appears to be restricted within the Nf1 gene interval. In Australia, the frequency of the Nf1‐insertion‐A haplotype increases with latitude where wing size is larger, independent of the chromosomal inversion In(3R)Payne. Unexpectedly, the Nf1‐insertion‐A haplotype is negatively associated with wing size. We found that the Nf1‐insertion‐A haplotype is enriched in females with shorter development time. This suggests that the Nf1 haplotype cline may be driven by selection for development time rather than size; females from southern (higher latitude) D. melanogaster populations maintain a rapid development time despite being relatively larger, and the higher incidence of Nf1‐insertion‐A in Southern Australia may contribute to this pattern, whereas the effects of the Nf1 haplotypes on size may be countered by other loci with antagonistic effects on size and development time. Our results point to the potential complexity involved in identifying selection on genetic variants exhibiting pleiotropic effects when studies are based on spatial patterns or association studies.     

Genotype‐environment interactions rule the response of a widespread butterfly to temperature variation

Understanding how organisms adapt to complex environments is a central goal of evolutionary biology and ecology. This issue is of special interest in the current era of rapidly changing climatic conditions. Here, we investigate clinal variation and plastic responses in life history, morphology and physiology in the butterfly Pieris napi along a pan‐European gradient by exposing butterflies raised in captivity to different temperatures. We found clinal variation in body size, growth rates and concomitant development time, wing aspect ratio, wing melanization and heat tolerance. Individuals from warmer environments were more heat‐tolerant and had less melanised wings and a shorter development, but still they were larger than individuals from cooler environments. These findings suggest selection for rapid growth in the warmth and for wing melanization in the cold, and thus fine‐tuned genetic adaptation to local climates. Irrespective of the origin of butterflies, the effects of higher developmental temperature were largely as expected, speeding up development; reducing body size, potential metabolic activity and wing melanization; while increasing heat tolerance. At least in part, these patterns likely reflect adaptive phenotypic plasticity. In summary, our study revealed pronounced plastic and genetic responses, which may indicate high adaptive capacities in our study organism. Whether this may help such species, though, to deal with current climate change needs further investigation, as clinal patterns have typically evolved over long periods.     

Genomic studies on the nature of species: adaptation and speciation in Mimulus

Evolutionary biology is in an exciting era, in which powerful genomic tools make the answers accessible to long‐standing questions about variation, adaptation and speciation. The availability of a suite of genomic resources, a shared knowledge base and a long history of study have made the phenotypically diverse plant genus Mimulus an important system for understanding ecological and evolutionary processes. An international Mimulus Research Meeting was held at Duke University in June 2014 to discuss developments in ecological and evolutionary genetic studies in Mimulus. Here, we report major recent discoveries presented at the meeting that use genomic approaches to advance our understanding of three major themes: the parallel genetic basis of adaptation; the ecological genomics of speciation; and the evolutionary significance of structural genetic variation. We also suggest future research directions for studies of Mimulus and highlight challenges faced when developing new ecological and evolutionary model systems.     

Adaptation,ancestral variation and gene flow in a ‘Sky Island’ Drosophila species

Over time, populations of species can expand, contract, fragment and become isolated, creating subpopulations that must adapt to local conditions. Understanding how species maintain variation after divergence as well as adapt to these changes in the face of gene flow is of great interest, especially as the current climate crisis has caused range shifts and frequent migrations for many species. Here, we characterize how a mycophageous fly species, Drosophila innubila, came to inhabit and adapt to its current range which includes mountain forests in south‐western USA separated by large expanses of desert. Using population genomic data from more than 300 wild‐caught individuals, we examine four populations to determine their population history in these mountain forests, looking for signatures of local adaptation. In this first extensive study, establishing D. innubila as a key genomic "Sky Island" model, we find D. innubila spread northwards during the previous glaciation period (30–100 KYA) and have recently expanded even further (0.2–2 KYA). D. innubila shows little evidence of population structure, consistent with a recent establishment and genetic variation maintained since before geographic stratification. We also find some signatures of recent selective sweeps in chorion proteins and population differentiation in antifungal immune genes suggesting differences in the environments to which flies are adapting. However, we find little support for long‐term recurrent selection in these genes. In contrast, we find evidence of long‐term recurrent positive selection in immune pathways such as the Toll signalling system and the Toll‐regulated antimicrobial peptides.     

Dissecting the role of a large chromosomal inversion in life history divergence throughout the Mimulus guttatus species complex

Chromosomal inversions can play an important role in adaptation, but the mechanism of their action in many natural populations remains unclear. An inversion could suppress recombination between locally beneficial alleles, thereby preventing maladaptive reshuffling with less‐fit, migrant alleles. The recombination suppression hypothesis has gained much theoretical support but empirical tests are lacking. Here, we evaluated the evolutionary history and phenotypic effects of a chromosomal inversion which differentiates annual and perennial forms of Mimulus guttatus. We found that perennials likely possess the derived orientation of the inversion. In addition, this perennial orientation occurs in a second perennial species, M. decorus, where it is strongly associated with life history differences between co‐occurring M. decorus and annual M. guttatus. One prediction of the recombination suppression hypothesis is that loci contributing to local adaptation will predate the inversion. To test whether the loci influencing perenniality pre‐date this inversion, we mapped QTLs for life history traits that differ between annual M. guttatus and a more distantly related, collinear perennial species, M. tilingii. Consistent with the recombination suppression hypothesis, we found that this region is associated with life history in the absence of the inversion, and this association can be broken into at least two QTLs. However, the absolute phenotypic effect of the LG8 inversion region on life history is weaker in M. tilingii than in perennials which possess the inversion. Thus, while we find support for the recombination suppression hypothesis, the contribution of this inversion to life history divergence in this group is likely complex.     

Successful despite poor flight performance: range expansion is associated with enhanced exploratory behaviour and fast development

Anthropogenic interference forces species to respond to changing environmental conditions. One possible response is dispersal and concomitant range shifts, allowing individuals to escape unfavourable conditions or to track the shifting climate niche. Range expansions depend on both dispersal capacity and the ability to establish populations beyond the former range. We here compare well‐established core populations with recently established edge populations in the currently northward expanding butterfly Lycaena tityrus. Edge populations were characterized by shorter development times and smaller size, a higher sensitivity to high temperature and an enhanced exploratory behaviour. The differences between core and edge populations found suggest adaptation to local climates and an enhanced dispersal ability in edge populations. In particular, enhanced exploratory behaviour may be advantageous in all steps of the dispersal process and may have facilitated the current range expansion. This study describes differences associated with a current range expansion, knowledge which might be useful for a better understanding of species responses to environmental change. We further report on variation between males and females in morphology and flight behaviour, with males showing a longer flight endurance and more pronounced exploratory behaviour than females.     

The roles of pre‐ and post‐hatching growth rates in generating a latitudinal cline of body size in the eastern fence lizard (Sceloporus undulatus)

Countergradient variation in norms of reaction can dampen the direct effects of environmental influences on phenotypic traits, allowing phenotypic similarity among populations despite exposure to different environmental conditions. Such norms of reaction may occur at any phase of the life‐history (e.g. growth rates during both embryonic and postembryonic stages may influence geographical variation in adult body size). We collected gravid female lizards (Sceloporus undulatus) from northern (Indiana), central (Mississippi), and southern (Florida) populations, spanning almost the full latitudinal range of the species. Adult females from the southern population were smaller. Intrinsic growth rates of hatchlings were higher for the central population than for the other two populations. This pattern does not parallel the countergradient variation previously found in embryonic developmental rates among these populations. Earlier hatching enhanced survival rates of juveniles to a similar degree among populations, although juvenile survival rates in the field generally increase with latitude in this species. Our data reveal geographical variation in the ways in which intrinsic developmental/growth rates and survival shift during ontogeny, and suggest that latitudinal patterns in adult body size (such as Bergmann's rule) can result from both faster growth, and longer periods of growth. © 2012 The Linnean Society of London, Biological Journal of the Linnean Society, 2012, 106 , 202–209.     

Adaptive divergence in body size overrides the effects of plasticity across natural habitats in the brown trout

The evolution of life-history traits is characterized by trade-offs between different selection pressures, as well as plasticity across environmental conditions. Yet, studies on local adaptation are often performed under artificial conditions, leaving two issues unexplored: (i) how consistent are laboratory inferred local adaptations under natural conditions and (ii) how much phenotypic variation is attributed to phenotypic plasticity and to adaptive evolution, respectively, across environmental conditions? We reared fish from six locally adapted (domesticated and wild) populations of anadromous brown trout (Salmo trutta) in one semi-natural and three natural streams and recorded a key life-history trait (body size at the end of first growth season). We found that population-specific reaction norms were close to parallel across different streams and Q_ST was similar – and larger than F_ST – within all streams, indicating a consistency of local adaptation in body size across natural environments. The amount of variation explained by population origin exceeded the variation across stream environments, indicating that genetic effects derived from adaptive processes have a stronger effect on phenotypic variation than plasticity induced by environmental conditions. These results suggest that plasticity does not “swamp” the phenotypic variation, and that selection may thus be efficient in generating genetic change.     

Adaptation to a seasonally varying environment: a strong latitudinal cline in reproductive diapause combined with high gene flow in Drosophila montana

Adaptation to seasonal changes in the northern hemisphere includes an ability to predict the forthcoming cold season from gradual changes in environmental cues early enough to prepare for the harsh winter conditions. The magnitude and speed of changes in these cues vary between the latitudes, which induces strong selection pressures for local adaptation.We studied adaptation to seasonal changes in Drosophila montana, a northern maltfly, by defining the photoperiodic conditions leading to adult reproductive diapause along a latitudinal cline in Finland and by measuring genetic differentiation and the amount of gene flow between the sampling sites with microsatellites. Our data revealed a clear correlation between the latitude and the critical day length (CDL), in which half of the females of different cline populations enter photoperiodic reproductive diapause. There was no sign of limited gene flow between the cline populations, even though these populations showed isolation by distance. Our results show that local adaptation may occur even in the presence of high gene flow, when selection for locally adaptive life-history traits is strong. A wide range of variation in the CDLs of the fly strains within and between the cline populations may be partly due to gene flow and partly due to the opposing selection pressures for fly reproduction and overwinter survival. This variation in the timing of diapause will enhance populations' survival over the years that differ in the severity of the winter and in the length of the warm period and may also help them respond to long-term changes in environmental conditions.     

Variation in soil aluminium tolerance genes is associated with local adaptation to soils at the Park Grass Experiment

Studies of the wild grass Anthoxanthum odoratum at the long‐term Park Grass Experiment (PGE, Harpenden, UK) document a well‐known example of rapid plant evolution in response to environmental change. Repeated fertilizer applications have acidified the soil in some experimental plots over the past 150+ years, and Anthoxanthum subpopulations have quickly become locally adapted. Early reciprocal transplants showed subpopulation differentiation specifically in response to soil aluminium (Al) toxicity across the experiment, even at small (30 m) spatial scales. Almost 40 years after its original measurement, we reassessed the degree of local adaptation to soil Al at the PGE using updated phenotyping methods and identified genes with variation linked to the tolerance trait. Root growth assays show that plants are locally adapted to soil Al at both the seedling and adult growth stages, but to a smaller extent than previously inferred. Among a large suite of candidate loci that were previously shown to have Al‐sensitive expression differences between sensitive and tolerant plants, three loci contained SNPs that are associated with both Al tolerance and soil acidity: an Al‐sensitive malate transporter (ALMT), a tonoplast intrinsic protein (TIP) and the putative homolog of the rice cell‐wall modification gene STAR1. Natural genetic variation at these loci is likely to have contributed to the recent rapid evolution at PGE. Continued study of Al tolerance variants in Anthoxanthum will allow us to test hypotheses about the nature and source of genetic variation that enables some species to adapt to soil acidification and other types of rapid environmental change.     

Symbiont community diversity is more variable in corals that respond poorly to stress

Coral reefs are declining globally as climate change and local water quality press environmental conditions beyond the physiological tolerances of holobionts—the collective of the host and its microbial symbionts. To assess the relationship between symbiont composition and holobiont stress tolerance, community diversity metrics were quantified for dinoflagellate endosymbionts (Family: Symbiodiniaceae) from eight Acropora millepora genets that thrived under or responded poorly to various stressors. These eight selected genets represent the upper and lower tails of the response distribution of 40 coral genets that were exposed to four stress treatments (and control conditions) in a 10‐day experiment. Specifically, four ‘best performer’ coral genets were analyzed at the end of the experiment because they survived high temperature, high pCO₂, bacterial exposure, or combined stressors, whereas four ‘worst performer’ genets were characterized because they experienced substantial mortality under these stressors. At the end of the experiment, seven of eight coral genets mainly hosted Cladocopium symbionts, whereas the eighth genet was dominated by both Cladocopium and Durusdinium symbionts. Symbiodiniaceae alpha and beta diversity were higher in worst performing genets than in best performing genets. Symbiont communities in worst performers also differed more after stress exposure relative to their controls (based on normalized proportional differences in beta diversity), than did best performers. A generalized joint attribute model estimated the influence of host genet and treatment on Symbiodiniaceae community composition and identified strong associations among particular symbionts and host genet performance, as well as weaker associations with treatment. Although dominant symbiont physiology and function contribute to host performance, these findings emphasize the importance of symbiont community diversity and stochasticity as components of host performance. Our findings also suggest that symbiont community diversity metrics may function as indicators of resilience and have potential applications in diverse disciplines from climate change adaptation to agriculture and medicine.     

Microgeographic morphological variation across larval wood frog populations associated with environment despite gene flow

Gene flow has historically been thought to constrain local adaptation; yet, recent research suggests that populations can diverge despite exchanging genes. Here I use a common garden experiment to assess the combined effects of gene flow and natural selection on morphological variation of 16 wood frog (Rana sylvatica) populations, a species known to experience divergent selection pressures in open‐ and closed‐canopy ponds across relatively small geographic scales. Wood frog tadpoles from different ponds showed significant morphological variation associated with canopy type with a trade‐off between tail length and body depth consistent with previous research. In contrast, neutral genetic differentiation of nine microsatellite loci as measured by Jost's D was not associated with canopy type, indicating no pattern of isolation by environment. Genetic structure analyses indicated some substructure across the 16 ponds (K = 4); however, three out of four assigned clusters included both open‐ and closed‐canopy ponds. Together, these results suggest that morphological divergence among these wood frog populations is occurring despite gene flow and that selection within these environments is strong. Furthermore, morphological variation among ponds differed across two sampling periods during larval development, demonstrating the importance of evaluating phenotypic divergence over multiple time periods and at a time relevant to the processes being studied.     

Progression in multiple sclerosis is associated with low endogenous NCAM

Multiple sclerosis (MS) is a CNS disorder characterized by demyelination and neurodegeneration. Although hallmarks of recovery (remyelination and repair) have been documented in early MS, the regenerative capacity of the adult CNS per se remains uncertain with the wide held belief that it is either limited or non‐existent. The neural cell adhesion molecule (NCAM) is a cell adhesion molecule that has been widely implicated in axonal outgrowth, guidance and fasciculation. Here, we used in vitro and in vivo of MS to investigate the role of NCAM in disease progression. We show that in health NCAM levels decrease over time, but this occurs acutely after demyelination and remains reduced in chronic disease. Our findings suggest that depletion of NCAM is one of the factors associated with or possibly responsible for disease progression in MS.     

Selection at behavioural,developmental and metabolic genes is associated with the northward expansion of a successful tropical colonizer

What makes a species able to colonize novel environments? This question is key to understand the dynamics of adaptive radiations and ecological niche shifts, but the mechanisms that underlie expansion into novel habitats remain poorly understood at a genomic scale. Lizards from the genus Anolis are typically tropical, and the green anole (Anolis carolinensis) constitutes an exception since it expanded into temperate North America from subtropical Florida. Thus, we used the green anole as a model to investigate signatures of selection associated with colonization of a new environment, namely temperate North America. To this end, we analysed 29 whole‐genome sequences, covering the entire native range of the species. We used a combination of recent methods to quantify both positive and balancing selection in northern populations, including F_ST outlier methods, machine learning and ancestral recombination graphs. We naively scanned for genes of interest and assessed the overlap between multiple tests. Strikingly, we identified many genes involved in behaviour, suggesting that the recent successful colonization of northern environments may have been linked to behavioural shifts as well as physiological adaptation. Using a candidate genes strategy, we determined that genes involved in response to cold or behaviour displayed more frequently signals of selection, while controlling for local recombination rate, gene clustering and gene length. In addition, we found signatures of balancing selection at immune genes in all investigated genetic groups, but also at genes involved in neuronal and anatomical development.