Genome architecture enables local adaptation of Atlantic cod despite high connectivity

Adaptation to local conditions is a fundamental process in evolution; however, mechanisms maintaining local adaptation despite high gene flow are still poorly understood. Marine ecosystems provide a wide array of diverse habitats that frequently promote ecological adaptation even in species characterized by strong levels of gene flow. As one example, populations of the marine fish Atlantic cod (Gadus morhua) are highly connected due to immense dispersal capabilities but nevertheless show local adaptation in several key traits. By combining population genomic analyses based on 12K single nucleotide polymorphisms with larval dispersal patterns inferred using a biophysical ocean model, we show that Atlantic cod individuals residing in sheltered estuarine habitats of Scandinavian fjords mainly belong to offshore oceanic populations with considerable connectivity between these diverse ecosystems. Nevertheless, we also find evidence for discrete fjord populations that are genetically differentiated from offshore populations, indicative of local adaptation, the degree of which appears to be influenced by connectivity. Analyses of the genomic architecture reveal a significant overrepresentation of a large ~5 Mb chromosomal rearrangement in fjord cod, previously proposed to comprise genes critical for the survival at low salinities. This suggests that despite considerable connectivity with offshore populations, local adaptation to fjord environments may be enabled by suppression of recombination in the rearranged region. Our study provides new insights into the potential of local adaptation in high gene flow species within fine geographical scales and highlights the importance of genome architecture in analyses of ecological adaptation.     

Genomic Architecture of Rapid Parallel Adaptation to Fresh Water in a Wild Fish

Rapid adaptation to novel environments may drive changes in genomic regions through natural selection. However, the genetic architecture underlying these adaptive changes is still poorly understood. Using population genomic approaches, we investigated the genomic architecture that underlies rapid parallel adaptation of Coilia nasus to fresh water by comparing four freshwater-resident populations with their ancestral anadromous population. Linkage disequilibrium network analysis and population genetic analyses revealed two putative large chromosome inversions on LG6 and LG22, which were enriched for outlier loci and exhibited parallel association with freshwater adaptation. Drastic frequency shifts and elevated genetic differentiation were observed for the two chromosome inversions among populations, suggesting that both inversions would undergo divergent selection between anadromous and resident ecotypes. Enrichment analysis of genes within chromosome inversions showed significant enrichment of genes involved in metabolic process, immunoregulation, growth, maturation, osmoregulation, and so forth, which probably underlay differences in morphology, physiology and behavior between the anadromous and freshwater-resident forms. The availability of beneficial standing genetic variation, large optimum shift between marine and freshwater habitats, and high efficiency of selection with large population size could lead to the observed rapid parallel adaptive genomic change. We propose that chromosomal inversions might have played an important role during the evolution of rapid parallel ecological divergence in the face of environmental heterogeneity in C. nasus. Our study provides insights into the genomic basis of rapid adaptation of complex traits in novel habitats and highlights the importance of structural genomic variants in analyses of ecological adaptation.     

Detecting adaptive evolution based on association with ecological gradients: Orientation matters!

Population genetic signatures of local adaptation are frequently investigated by identifying loci with allele frequencies that exhibit high correlation with ecological variables. One difficulty with this approach is that ecological associations might be confounded by geographic variation at selectively neutral loci. Here, we consider populations that underwent spatial expansion from their original range, and for which geographical variation of adaptive allele frequency coincides with habitat gradients. Using range expansion simulations, we asked whether our ability to detect genomic regions involved in adaptation could be impacted by the orientation of the ecological gradients. For three ecological association methods tested, we found, counter-intuitively, fewer false-positive associations when ecological gradients aligned along the main axis of expansion than when they aligned along any other direction. This result has important consequences for the analysis of genomic data under non-equilibrium population genetic models. Alignment of gradients with expansion axes is likely to be common in scenarios in which expanding species track their ecological niche during climate change while adapting to changing environments at their rear edge.     

Ancient chromosomal rearrangement associated with local adaptation of a postglacially colonized population of Atlantic Cod in the northwest Atlantic

Intraspecific diversity is central to the management and conservation of exploited species, yet knowledge of how this diversity is distributed and maintained in the genome of many marine species is lacking. Recent advances in genomic analyses allow for genome‐wide surveys of intraspecific diversity and offer new opportunities for exploring genomic patterns of divergence. Here, we analysed genome‐wide polymorphisms to measure genetic differentiation between an offshore migratory and a nonmigratory population and to define conservation units of Atlantic Cod (Gadus morhua) in coastal Labrador. A total of 141 individuals, collected from offshore sites and from a coastal site within Gilbert Bay, Labrador, were genotyped using an ~11k single nucleotide polymorphism array. Analyses of population structure revealed strong genetic differentiation between migratory offshore cod and nonmigratory Gilbert Bay cod. Genetic differentiation was elevated for loci within a chromosomal rearrangement found on linkage group 1 (LG1) that coincides with a previously found double inversion associated with migratory and nonmigratory ecotype divergence of cod in the northeast Atlantic. This inverted region includes several genes potentially associated with adaptation to differences in salinity and temperature, as well as influencing migratory behaviour. Our work provides evidence that a chromosomal rearrangement on LG1 is associated with parallel patterns of divergence between migratory and nonmigratory ecotypes on both sides of the Atlantic Ocean.     

Testing for local adaptation and evolutionary potential along altitudinal gradients in rainforest Drosophila: beyond laboratory estimates

Predicting how species will respond to the rapid climatic changes predicted this century is an urgent task. Species distribution models (SDMs) use the current relationship between environmental variation and species’ abundances to predict the effect of future environmental change on their distributions. However, two common assumptions of SDMs are likely to be violated in many cases: (i) that the relationship of environment with abundance or fitness is constant throughout a species’ range and will remain so in future and (ii) that abiotic factors (e.g. temperature, humidity) determine species’ distributions. We test these assumptions by relating field abundance of the rainforest fruit fly Drosophila birchii to ecological change across gradients that include its low and high altitudinal limits. We then test how such ecological variation affects the fitness of 35 D. birchii families transplanted in 591 cages to sites along two altitudinal gradients, to determine whether genetic variation in fitness responses could facilitate future adaptation to environmental change. Overall, field abundance was highest at cooler, high‐altitude sites, and declined towards warmer, low‐altitude sites. By contrast, cage fitness (productivity) increased towards warmer, lower‐altitude sites, suggesting that biotic interactions (absent from cages) drive ecological limits at warmer margins. In addition, the relationship between environmental variation and abundance varied significantly among gradients, indicating divergence in ecological niche across the species’ range. However, there was no evidence for local adaptation within gradients, despite greater productivity of high‐altitude than low‐altitude populations when families were reared under laboratory conditions. Families also responded similarly to transplantation along gradients, providing no evidence for fitness trade‐offs that would favour local adaptation. These findings highlight the importance of (i) measuring genetic variation in key traits under ecologically relevant conditions, and (ii) considering the effect of biotic interactions when predicting species’ responses to environmental change.     

Chromosome polymorphisms track trans‐Atlantic divergence and secondary contact in Atlantic salmon

Pleistocene glaciations drove repeated range contractions and expansions shaping contemporary intraspecific diversity. Atlantic salmon (Salmo salar) in the western and eastern Atlantic diverged >600,000 years before present, with the two lineages isolated in different southern refugia during glacial maxima, driving trans‐Atlantic genomic and karyotypic divergence. Here, we investigate the genomic consequences of glacial isolation and trans‐Atlantic secondary contact using 108,870 single nucleotide polymorphisms genotyped in 80 North American and European populations. Throughout North America, we identified extensive interindividual variation and discrete linkage blocks within and between chromosomes with known trans‐Atlantic differences in rearrangements: Ssa01/Ssa23 translocation and Ssa08/Ssa29 fusion. Spatial genetic analyses suggest independence of rearrangements, with Ssa01/Ssa23 showing high European introgression (>50%) in northern populations indicative of post‐glacial trans‐Atlantic secondary contact, contrasting with low European ancestry genome‐wide (3%). Ssa08/Ssa29 showed greater intrapopulation diversity, suggesting a derived chromosome fusion polymorphism that evolved within North America. Evidence of potential selection on both genomic regions suggests that the adaptive role of rearrangements warrants further investigation in Atlantic salmon. Our study highlights how Pleistocene glaciations can influence large‐scale intraspecific variation in genomic architecture of northern species.     

Ecological genomics meets community‐level modelling of biodiversity: mapping the genomic landscape of current and future environmental adaptation

Local adaptation is a central feature of most species occupying spatially heterogeneous environments, and may factor critically in responses to environmental change. However, most efforts to model the response of species to climate change ignore intraspecific variation due to local adaptation. Here, we present a new perspective on spatial modelling of organism–environment relationships that combines genomic data and community‐level modelling to develop scenarios regarding the geographic distribution of genomic variation in response to environmental change. Rather than modelling species within communities, we use these techniques to model large numbers of loci across genomes. Using balsam poplar (Populus balsamifera) as a case study, we demonstrate how our framework can accommodate nonlinear responses of loci to environmental gradients. We identify a threshold response to temperature in the circadian clock gene GIGANTEA‐5 (GI5), suggesting that this gene has experienced strong local adaptation to temperature. We also demonstrate how these methods can map ecological adaptation from genomic data, including the identification of predicted differences in the genetic composition of populations under current and future climates. Community‐level modelling of genomic variation represents an important advance in landscape genomics and spatial modelling of biodiversity that moves beyond species‐level assessments of climate change vulnerability.     

Landscape genomics and a common garden trial reveal adaptive differentiation to temperature across Europe in the tree species Alnus glutinosa

The adaptive potential of tree species to cope with climate change has important ecological and economic implications. Many temperate tree species experience a wide range of environmental conditions, suggesting high adaptability to new environmental conditions. We investigated adaptation to regional climate in the drought‐sensitive tree species Alnus glutinosa (Black alder), using a complementary approach that integrates genomic, phenotypic and landscape data. A total of 24 European populations were studied in a common garden and through landscape genomic approaches. Genotyping‐by‐sequencing was used to identify SNPs across the genome, resulting in 1990 SNPs. Although a relatively low percentage of putative adaptive SNPs was detected (2.86% outlier SNPs), we observed clear associations among outlier allele frequencies, temperature and plant traits. In line with the typical drought avoiding nature of A. glutinosa, leaf size varied according to a temperature gradient and significant associations with multiple outlier loci were observed, corroborating the ecological relevance of the observed outlier SNPs. Moreover, the lack of isolation by distance, the very low genetic differentiation among populations and the high intrapopulation genetic variation all support the notion that high gene exchange combined with strong environmental selection promotes adaptation to environmental cues.     

Environmental gradients of selection for an alpine-obligate bird,the white-tailed ptarmigan (Lagopus leucura)

The warming climate will expose alpine species adapted to a highly seasonal, harsh environment to novel environmental conditions. A species can shift their distribution, acclimate, or adapt in response to a new climate. Alpine species have little suitable habitat to shift their distribution, and the limits of acclimation will likely be tested by climate change in the long-term. Adaptive genetic variation may provide the raw material for species to adapt to this changing environment. Here, we use a genomic approach to describe adaptive divergence in an alpine-obligate species, the white-tailed ptarmigan (Lagopus leucura), a species distributed from Alaska to New Mexico, across an environmentally variable geographic range. Previous work has identified genetic structure and morphological, behavioral, and physiological differences across the species’ range; however, those studies were unable to determine the degree to which adaptive divergence is correlated with local variation in environmental conditions. We used a genome-wide dataset generated from 95 white-tailed ptarmigan distributed throughout the species’ range and genotype–environment association analyses to identify the genetic signature and environmental drivers of local adaptation. We detected associations between multiple environmental gradients and candidate adaptive loci, suggesting ptarmigan populations may be locally adapted to the plant community composition, elevation, local climate, and to the seasonality of the environment. Overall, our results suggest there may be groups within the species’ range with genetic variation that could be essential for adapting to a changing climate and helpful in guiding conservation action.Subject terms: Ecological genetics, Evolutionary ecology     

Signatures of polygenic adaptation associated with climate across the range of a threatened fish species with high genetic connectivity

Adaptive differences across species’ ranges can have important implications for population persistence and conservation management decisions. Despite advances in genomic technologies, detecting adaptive variation in natural populations remains challenging. Key challenges in gene–environment association studies involve distinguishing the effects of drift from those of selection and identifying subtle signatures of polygenic adaptation. We used paired‐end restriction site‐associated DNA sequencing data (6,605 biallelic single nucleotide polymorphisms; SNPs) to examine population structure and test for signatures of adaptation across the geographic range of an iconic Australian endemic freshwater fish species, the Murray cod Maccullochella peelii. Two univariate gene–association methods identified 61 genomic regions associated with climate variation. We also tested for subtle signatures of polygenic adaptation using a multivariate method (redundancy analysis; RDA). The RDA analysis suggested that climate (temperature‐ and precipitation‐related variables) and geography had similar magnitudes of effect in shaping the distribution of SNP genotypes across the sampled range of Murray cod. Although there was poor agreement among the candidate SNPs identified by the univariate methods, the top 5% of SNPs contributing to significant RDA axes included 67% of the SNPs identified by univariate methods. We discuss the potential implications of our findings for the management of Murray cod and other species generally, particularly in relation to informing conservation actions such as translocations to improve evolutionary resilience of natural populations. Our results highlight the value of using a combination of different approaches, including polygenic methods, when testing for signatures of adaptation in landscape genomic studies.     

Disentangling structural genomic and behavioural barriers in a sea of connectivity

Genetic divergence among populations arises through natural selection or drift and is counteracted by connectivity and gene flow. In sympatric populations, isolating mechanisms are thus needed to limit the homogenizing effects of gene flow to allow for adaptation and speciation. Chromosomal inversions act as an important mechanism maintaining isolating barriers, yet their role in sympatric populations and divergence with gene flow is not entirely understood. Here, we revisit the question of whether inversions play a role in the divergence of connected populations of the marine fish Atlantic cod (Gadus morhua), by exploring a unique data set combining whole‐genome sequencing data and behavioural data obtained with acoustic telemetry. Within a confined fjord environment, we find three genetically differentiated Atlantic cod types belonging to the oceanic North Sea population, the western Baltic population and a local fjord‐type cod. Continuous behavioural tracking over 4 year revealed temporally stable sympatry of these types within the fjord. Despite overall weak genetic differentiation consistent with high levels of gene flow, we detected significant frequency shifts of three previously identified inversions, indicating an adaptive barrier to gene flow. In addition, behavioural data indicated that North Sea cod and individuals homozygous for the LG12 inversion had lower fitness in the fjord environment. However, North Sea and fjord‐type cod also occupy different depths, possibly contributing to prezygotic reproductive isolation and representing a behavioural barrier to gene flow. Our results provide the first insights into a complex interplay of genomic and behavioural isolating barriers in Atlantic cod and establish a new model system towards an understanding of the role of genomic structural variants in adaptation and diversification.     

Broad-scale adaptive genetic variation in alpine plants is driven by temperature and precipitation

Identifying adaptive genetic variation is a challenging task, in particular in non-model species for which genomic information is still limited or absent. Here, we studied distribution patterns of amplified fragment length polymorphisms (AFLPs) in response to environmental variation, in 13 alpine plant species consistently sampled across the entire European Alps. Multiple linear regressions were performed between AFLP allele frequencies per site as dependent variables and two categories of independent variables, namely Moran's eigenvector map MEM variables (to account for spatial and unaccounted environmental variation, and historical demographic processes) and environmental variables. These associations allowed the identification of 153 loci of ecological relevance. Univariate regressions between allele frequency and each environmental factor further showed that loci of ecological relevance were mainly correlated with MEM variables. We found that precipitation and temperature were the best environmental predictors, whereas topographic factors were rarely involved in environmental associations. Climatic factors, subject to rapid variation as a result of the current global warming, are known to strongly influence the fate of alpine plants. Our study shows, for the first time for a large number of species, that the same environmental variables are drivers of plant adaptation at the scale of a whole biome, here the European Alps.     

Locally Adaptive Inversions Modulate Genetic Variation at Different Geographic Scales in a Seaweed Fly

Across a species range, multiple sources of environmental heterogeneity, at both small and large scales, create complex landscapes of selection, which may challenge adaptation, particularly when gene flow is high. One key to multidimensional adaptation may reside in the heterogeneity of recombination along the genome. Structural variants, like chromosomal inversions, reduce recombination, increasing linkage disequilibrium among loci at a potentially massive scale. In this study, we examined how chromosomal inversions shape genetic variation across a species range and ask how their contribution to adaptation in the face of gene flow varies across geographic scales. We sampled the seaweed fly Coelopa frigida along a bioclimatic gradient stretching across 10° of latitude, a salinity gradient, and a range of heterogeneous, patchy habitats. We generated a chromosome-level genome assembly to analyze 1,446 low-coverage whole genomes collected along those gradients. We found several large nonrecombining genomic regions, including putative inversions. In contrast to the collinear regions, inversions and low-recombining regions differentiated populations more strongly, either along an ecogeographic cline or at a fine-grained scale. These genomic regions were associated with environmental factors and adaptive phenotypes, albeit with contrasting patterns. Altogether, our results highlight the importance of recombination in shaping adaptation to environmental heterogeneity at local and large scales.     

Broad Concordance in the Spatial Distribution of Adaptive and Neutral Genetic Variation across an Elevational Gradient in Deer Mice

When species are continuously distributed across environmental gradients, the relative strength of selection and gene flow shape spatial patterns of genetic variation, potentially leading to variable levels of differentiation across loci. Determining whether adaptive genetic variation tends to be structured differently than neutral variation along environmental gradients is an open and important question in evolutionary genetics. We performed exome-wide population genomic analysis on deer mice sampled along an elevational gradient of nearly 4,000 m of vertical relief. Using a combination of selection scans, genotype−environment associations, and geographic cline analyses, we found that a large proportion of the exome has experienced a history of altitude-related selection. Elevational clines for nearly 30% of these putatively adaptive loci were shifted significantly up- or downslope of clines for loci that did not bear similar signatures of selection. Many of these selection targets can be plausibly linked to known phenotypic differences between highland and lowland deer mice, although the vast majority of these candidates have not been reported in other studies of highland taxa. Together, these results suggest new hypotheses about the genetic basis of physiological adaptation to high altitude, and the spatial distribution of adaptive genetic variation along environmental gradients.     

The behavioural diversity of Atlantic cod: insights into variability within and between individuals

Cod (Gadus morhua) are an iconic fish species of cultural, historical and economical significance across the Atlantic and adjacent seas. Among many scholarly investigations, this interest has prompted behavioural research, rendering cod one of the few commercially harvested marine fishes for which behaviour has been studied in a comprehensive manner. In our review of this behavioural work, we examine the variability in cod behaviour across five functional domains: foraging, predation, social interactions, migration and reproduction. Research to date suggests a high level of behavioural sophistication in cod that is underpinned by complex learning strategies and long-term memory. Cod also demonstrate substantial variability in how they respond to different ecological circumstances. Considerable variation is evident both within and between individuals, and in some instances, between populations. There are a number of pathways from which this variation appears to arise, such as asocial and social learning, environmental control of phenotypic plasticity and genetic control, but there are no known examples of behaviours that are purely the result of one of these mechanisms. Behavioural variation is therefore likely to result from a combination of these factors, underscoring the need for a quantitative, multivariate approach to understand behavioural variation in cod.     

Shared ancestral polymorphisms and chromosomal rearrangements as potential drivers of local adaptation in a marine fish

Gene flow has tremendous importance for local adaptation, by influencing the fate of de novo mutations, maintaining standing genetic variation and driving adaptive introgression. Furthermore, structural variation as chromosomal rearrangements may facilitate adaptation despite high gene flow. However, our understanding of the evolutionary mechanisms impending or favouring local adaptation in the presence of gene flow is still limited to a restricted number of study systems. In this study, we examined how demographic history, shared ancestral polymorphism, and gene flow among glacial lineages contribute to local adaptation to sea conditions in a marine fish, the capelin (Mallotus villosus). We first assembled a 490‐Mbp draft genome of M. villosus to map our RAD sequence reads. Then, we used a large data set of genome‐wide single nucleotide polymorphisms (25,904 filtered SNPs) genotyped in 1,310 individuals collected from 31 spawning sites in the northwest Atlantic. We reconstructed the history of divergence among three glacial lineages and showed that they probably diverged from 3.8 to 1.8 million years ago and experienced secondary contacts. Within each lineage, our analyses provided evidence for large N_e and high gene flow among spawning sites. Within the Northwest Atlantic lineage, we detected a polymorphic chromosomal rearrangement leading to the occurrence of three haplogroups. Genotype–environment associations revealed molecular signatures of local adaptation to environmental conditions prevailing at spawning sites. Our study also suggests that both shared polymorphisms among lineages, resulting from standing genetic variation or introgression, and chromosomal rearrangements may contribute to local adaptation in the presence of high gene flow.     

LANDSCAPE GENOMICS IN ATLANTIC SALMON (SALMO SALAR): SEARCHING FOR GENE–ENVIRONMENT INTERACTIONS DRIVING LOCAL ADAPTATION

A growing number of studies are examining the factors driving historical and contemporary evolution in wild populations. By combining surveys of genomic variation with a comprehensive assessment of environmental parameters, such studies can increase our understanding of the genomic and geographical extent of local adaptation in wild populations. We used a large‐scale landscape genomics approach to examine adaptive and neutral differentiation across 54 North American populations of Atlantic salmon representing seven previously defined genetically distinct regional groups. Over 5500 genome‐wide single nucleotide polymorphisms were genotyped in 641 individuals and 28 bulk assays of 25 pooled individuals each. Genome scans, linkage map, and 49 environmental variables were combined to conduct an innovative landscape genomic analysis. Our results provide valuable insight into the links between environmental variation and both neutral and potentially adaptive genetic divergence. In particular, we identified markers potentially under divergent selection, as well as associated selective environmental factors and biological functions with the observed adaptive divergence. Multivariate landscape genetic analysis revealed strong associations of both genetic and environmental structures. We found an enrichment of growth‐related functions among outlier markers. Climate (temperature–precipitation) and geological characteristics were significantly associated with both potentially adaptive and neutral genetic divergence and should be considered as candidate loci involved in adaptation at the regional scale in Atlantic salmon. Hence, this study significantly contributes to the improvement of tools used in modern conservation and management schemes of Atlantic salmon wild populations.     

Multiscale landscape genomic models to detect signatures of selection in the alpine plant Biscutella laevigata

Plant species are known to adapt locally to their environment, particularly in mountainous areas where conditions can vary drastically over short distances. The climate of such landscapes being largely influenced by topography, using fine‐scale models to evaluate environmental heterogeneity may help detecting adaptation to micro‐habitats. Here, we applied a multiscale landscape genomic approach to detect evidence of local adaptation in the alpine plant Biscutella laevigata. The two gene pools identified, experiencing limited gene flow along a 1‐km ridge, were different in regard to several habitat features derived from a very high resolution (VHR) digital elevation model (DEM). A correlative approach detected signatures of selection along environmental gradients such as altitude, wind exposure, and solar radiation, indicating adaptive pressures likely driven by fine‐scale topography. Using a large panel of DEM‐derived variables as ecologically relevant proxies, our results highlighted the critical role of spatial resolution. These high‐resolution multiscale variables indeed indicate that the robustness of associations between genetic loci and environmental features depends on spatial parameters that are poorly documented. We argue that the scale issue is critical in landscape genomics and that multiscale ecological variables are key to improve our understanding of local adaptation in highly heterogeneous landscapes.     

Habitat variability correlates with duplicate content of Drosophila genomes

The factors limiting the habitat range of species are crucial in understanding their biodiversity and response to environmental change. Yet the genetic and genomic architectures that produce genetic variation to enable environmental adaptation have remained poorly understood. Here we show that the proportion of duplicated genes (P(D)) in the whole genomes of fully sequenced Drosophila species is significantly correlated with environmental variability within the habitats measured by the climatic envelope and habitat diversity. Furthermore, species with a low P(D) tend to lose the duplicated genes owing to their faster evolution. These results indicate that the rapid relaxation of functional constraints on duplicated genes resulted in a low P(D) for species with lower habitat diversity, and suggest that the maintenance of duplicated genes gives organisms an ecological advantage during evolution. We therefore propose that the P(D) in a genome is related to adaptation to environmental variation.     

Characterisation of Adaptive Genetic Diversity in Environmentally Contrasted Populations of Eucalyptus camaldulensis Dehnh. (River Red Gum)

As an increasing number of ecosystems face departures from long standing environmental conditions under climate change, our understanding of the capacity of species to adapt will become important for directing conservation and management of biodiversity. Insights into the potential for genetic adaptation might be gained by assessing genomic signatures of adaptation to historic or prevailing environmental conditions. The river red gum (Eucalyptus camaldulensis Dehnh.) is a widespread Australian eucalypt inhabiting riverine and floodplain habitats which spans strong environmental gradients. We investigated the effects of adaptation to environment on population level genetic diversity of E. camaldulensis, examining SNP variation in candidate gene loci sampled across 20 climatically diverse populations approximating the species natural distribution. Genetic differentiation among populations was high (F_ST = 17%), exceeding previous estimates based on neutral markers. Complementary statistical approaches identified 6 SNP loci in four genes (COMT, Dehydrin, ERECTA and PIP2) which, after accounting for demographic effects, exhibited higher than expected levels of genetic differentiation among populations and whose allelic variation was associated with local environment. While this study employs but a small proportion of available diversity in the eucalyptus genome, it draws our attention to the potential for application of wide spread eucalypt species to test adaptive hypotheses.