Genome‐wide single nucleotide polymorphisms reveal population history and adaptive divergence in wild guppies

Adaptation of guppies (Poecilia reticulata) to contrasting upland and lowland habitats has been extensively studied with respect to behaviour, morphology and life history traits. Yet population history has not been studied at the whole‐genome level. Although single nucleotide polymorphisms (SNPs) are the most abundant form of variation in many genomes and consequently very informative for a genome‐wide picture of standing natural variation in populations, genome‐wide SNP data are rarely available for wild vertebrates. Here we use genetically mapped SNP markers to comprehensively survey genetic variation within and among naturally occurring guppy populations from a wide geographic range in Trinidad and Venezuela. Results from three different clustering methods, Neighbor‐net, principal component analysis (PCA) and Bayesian analysis show that the population substructure agrees with geographic separation and largely with previously hypothesized patterns of historical colonization. Within major drainages (Caroni, Oropouche and Northern), populations are genetically similar, but those in different geographic regions are highly divergent from one another, with some indications of ancient shared polymorphisms. Clear genomic signatures of a previous introduction experiment were seen, and we detected additional potential admixture events. Headwater populations were significantly less heterozygous than downstream populations. Pairwise F_ST values revealed marked differences in allele frequencies among populations from different regions, and also among populations within the same region. F_ST outlier methods indicated some regions of the genome as being under directional selection. Overall, this study demonstrates the power of a genome‐wide SNP data set to inform for studies on natural variation, adaptation and evolution of wild populations     

Analysis of F ST outliers at allozyme loci in Pacific salmon: implications for natural selection

Natural selection has been invoked to explain the observed geographic distribution of allozyme allele frequencies for a number of teleost species. The effects of selection on allozyme loci in three species of Pacific salmon were tested. A simulation-based approach to estimate the null distribution of population differentiation (F _ST) and test for F _ST outliers was used. This approach showed that a majority of allozyme loci conform to neutral expectations predicted by the simulation model, with relatively few F _ST outliers found. No consistent F _ST outlier loci were found across species. Analysis of population sub-groups based on geography and genetic identity reduced the number of outlier loci for some species, indicating that large geographic groups may include genetically divergent populations and/or that there is geographic heterogeneity in selection pressure upon allozyme loci. Two outlier allozyme loci found in this analysis, lactate dehydrogenase-B and malic enzyme, have been shown to be influenced by selection in other teleost species. This approach is also useful in identifying allozyme loci (or other genetic markers) that meet assumptions for population genetic study.     

Adaptive divergence despite strong genetic drift: genomic analysis of the evolutionary mechanisms causing genetic differentiation in the island fox (Urocyon littoralis)

The evolutionary mechanisms generating the tremendous biodiversity of islands have long fascinated evolutionary biologists. Genetic drift and divergent selection are predicted to be strong on islands and both could drive population divergence and speciation. Alternatively, strong genetic drift may preclude adaptation. We conducted a genomic analysis to test the roles of genetic drift and divergent selection in causing genetic differentiation among populations of the island fox (Urocyon littoralis). This species consists of six subspecies, each of which occupies a different California Channel Island. Analysis of 5293 SNP loci generated using Restriction‐site Associated DNA (RAD) sequencing found support for genetic drift as the dominant evolutionary mechanism driving population divergence among island fox populations. In particular, populations had exceptionally low genetic variation, small N_e (range = 2.1–89.7; median = 19.4), and significant genetic signatures of bottlenecks. Moreover, islands with the lowest genetic variation (and, by inference, the strongest historical genetic drift) were most genetically differentiated from mainland grey foxes, and vice versa, indicating genetic drift drives genome‐wide divergence. Nonetheless, outlier tests identified 3.6–6.6% of loci as high F_ST outliers, suggesting that despite strong genetic drift, divergent selection contributes to population divergence. Patterns of similarity among populations based on high F_ST outliers mirrored patterns based on morphology, providing additional evidence that outliers reflect adaptive divergence. Extremely low genetic variation and small N_e in some island fox populations, particularly on San Nicolas Island, suggest that they may be vulnerable to fixation of deleterious alleles, decreased fitness and reduced adaptive potential.     

High degree of genetic differentiation in marine three‐spined sticklebacks (Gasterosteus aculeatus)

Populations of widespread marine organisms are typically characterized by a low degree of genetic differentiation in neutral genetic markers, but much less is known about differentiation in genes whose functional roles are associated with specific selection regimes. To uncover possible adaptive population divergence and heterogeneous genomic differentiation in marine three‐spined sticklebacks (Gasterosteus aculeatus), we used a candidate gene‐based genome‐scan approach to analyse variability in 138 microsatellite loci located within/close to (<6 kb) functionally important genes in samples collected from ten geographic locations. The degree of genetic differentiation in markers classified as neutral or under balancing selection—as determined with several outlier detection methods—was low (F_ST = 0.033 or 0.011, respectively), whereas average F_ST for directionally selected markers was significantly higher (F_ST = 0.097). Clustering analyses provided support for genomic and geographic heterogeneity in selection: six genetic clusters were identified based on allele frequency differences in the directionally selected loci, whereas four were identified with the neutral loci. Allelic variation in several loci exhibited significant associations with environmental variables, supporting the conjecture that temperature and salinity, but not optic conditions, are important drivers of adaptive divergence among populations. In general, these results suggest that in spite of the high degree of physical connectivity and gene flow as inferred from neutral marker genes, marine stickleback populations are strongly genetically structured in loci associated with functionally relevant genes.     

Evidence of divergent selection for drought and cold tolerance at landscape and local scales in Abies alba Mill. in the French Mediterranean Alps

Understanding local adaptation in forest trees is currently a key research and societal priority. Geographically and ecologically marginal populations provide ideal case studies, because environmental stress along with reduced gene flow can facilitate the establishment of locally adapted populations. We sampled European silver fir (Abies alba Mill.) trees in the French Mediterranean Alps, along the margin of its distribution range, from pairs of high‐ and low‐elevation plots on four different mountains situated along a 170‐km east–west transect. The analysis of 267 SNP loci from 175 candidate genes suggested a neutral pattern of east–west isolation by distance among mountain sites. F_ST outlier tests revealed 16 SNPs that showed patterns of divergent selection. Plot climate was characterized using both in situ measurements and gridded data that revealed marked differences between and within mountains with different trends depending on the season. Association between allelic frequencies and bioclimatic variables revealed eight genes that contained candidate SNPs, of which two were also detected using F_ST outlier methods. All SNPs were associated with winter drought, and one of them showed strong evidence of selection with respect to elevation. Q_ST–F_ST tests for fitness‐related traits measured in a common garden suggested adaptive divergence for the date of bud flush and for growth rate. Overall, our results suggest a complex adaptive picture for A. alba in the southern French Alps where, during the east‐to‐west Holocene recolonization, locally advantageous genetic variants established at both the landscape and local scales.     

Signatures of diversifying selection at EST‐SSR loci and association with climate in natural Eucalyptus populations

Understanding the environmental parameters that drive adaptation among populations is important in predicting how species may respond to global climatic changes and how gene pools might be managed to conserve adaptive genetic diversity. Here, we used Bayesian F_ST outlier tests and allele–climate association analyses to reveal two Eucalyptus EST‐SSR loci as strong candidates for diversifying selection in natural populations of a southwestern Australian forest tree, Eucalyptus gomphocephala (Myrtaceae). The Eucalyptus homolog of a CONSTANS‐like gene was an F_ST outlier, and allelic variation showed significant latitudinal clinal associations with annual and winter solar radiation, potential evaporation, summer precipitation and aridity. A second F_ST outlier locus, homologous to quinone oxidoreductase, was significantly associated with measures of temperature range, high summer temperature and summer solar radiation, with important implications for predicting the effect of temperature on natural populations in the context of climate change. We complemented these data with investigations into neutral population genetic structure and diversity throughout the species range. This study provides an investigation into selection signatures at gene‐homologous EST‐SSRs in natural Eucalyptus populations, and contributes to our understanding of the relationship between climate and adaptive genetic variation, informing the conservation of both putatively neutral and adaptive components of genetic diversity.     

Development and preliminary evaluation of a genomewide single nucleotide polymorphisms resource generated by RAD‐seq for the small yellow croaker (Larimichthys polyactis)

Recent advances in high‐throughput sequencing technologies have offered the possibility to generate genomewide sequence data to delineate previously unidentified genetic structure, obtain more accurate estimates of demographic parameters and to evaluate potential adaptive divergence. Here, we identified 27 556 single nucleotide polymorphisms for the small yellow croaker (Larimichthys polyactis) using restriction‐site‐associated DNA (RAD) sequencing of 24 individuals from two populations. Significant sources of genetic variation were identified, with an average nucleotide diversity (π) of 0.00105 ± 0.000425 across individuals, and long‐term effective population size was thus estimated to range between 26 172 and 261 716. According to the results, no differentiation between the two populations was detected based on the SNP data set of top quality score per contig or neutral loci. However, the two analysed populations were highly differentiated based on SNP data set of both top F_ST value per contig and the outlier SNPs. Moreover, local adaptation was highlighted by an F_ST‐based outlier tests implemented in LOSITAN and a total of 538 potentially locally selected SNPs were identified. blast2go annotation of contigs containing the outlier SNPs yielded hits for 37 (66%) of 56 significant blastx matches. Candidate genes for local adaptation constituted a wide array of biological functions, including cellular response to oxidative stress, actin filament binding, ion transmembrane transport and synapse assembly. The generated SNP resources in this study provided a valuable tool for future population genetics and genomics studies of L. polyactis.     

Targeted resequencing reveals geographical patterns of differentiation for loci implicated in parallel evolution

Parallel divergence and speciation provide evidence for the role of divergent selection in generating biological diversity. Recent studies indicate that parallel phenotypic divergence may not have the same genetic basis in different geographical locations – ‘outlier loci’ (loci potentially affected by divergent selection) are often not shared among parallel instances of phenotypic divergence. However, limited sharing may be due, in part, to technical issues if false‐positive outliers occur. Here, we test this idea in the marine snail Littorina saxatilis, which has evolved two partly isolated ecotypes (adapted to crab predation vs. wave action) in multiple locations independently. We argue that if the low extent of sharing observed in earlier studies in this system is due to sampling effects, we expect outliers not to show elevated F_ST when sequenced in new samples from the original locations and also not to follow predictable geographical patterns of elevated F_ST. Following a hierarchical sampling design (within vs. between country), we applied capture sequencing, targeting outliers from earlier studies and control loci. We found that outliers again showed elevated levels of F_ST in their original location, suggesting they were not generated by sampling effects. Outliers were also likely to show increased F_ST in geographically close locations, which may be explained by higher levels of gene flow or shared ancestral genetic variation compared with more distant locations. However, in contrast to earlier findings, we also found some outlier types to show elevated F_ST in geographically distant locations. We discuss possible explanations for this unexpected result.     

RAD sequencing reveals genomewide divergence between independent invasions of the European green crab (Carcinus maenas) in the Northwest Atlantic

Genomic studies of invasive species can reveal both invasive pathways and functional differences underpinning patterns of colonization success. The European green crab (Carcinus maenas) was initially introduced to eastern North America nearly 200 years ago where it expanded northwards to eastern Nova Scotia. A subsequent invasion to Nova Scotia from a northern European source allowed further range expansion, providing a unique opportunity to study the invasion genomics of a species with multiple invasions. Here, we use restriction‐site‐associated DNA sequencing‐derived SNPs to explore fine‐scale genomewide differentiation between these two invasions. We identified 9137 loci from green crab sampled from 11 locations along eastern North America and compared spatial variation to mitochondrial COI sequence variation used previously to characterize these invasions. Overall spatial divergence among invasions was high (pairwise F_ST ~0.001 to 0.15) and spread across many loci, with a mean F_ST ~0.052 and 52% of loci examined characterized by F_ST values >0.05. The majority of the most divergent loci (i.e., outliers, ~1.2%) displayed latitudinal clines in allele frequency highlighting extensive genomic divergence among the invasions. Discriminant analysis of principal components (both neutral and outlier loci) clearly resolved the two invasions spatially and was highly correlated with mitochondrial divergence. Our results reveal extensive cryptic intraspecific genomic diversity associated with differing patterns of colonization success and demonstrates clear utility for genomic approaches to delineating the distribution and colonization success of aquatic invasive species.     

Background selection and FST: Consequences for detecting local adaptation

Background selection is a process whereby recurrent deleterious mutations cause a decrease in the effective population size and genetic diversity at linked loci. Several authors have suggested that variation in the intensity of background selection could cause variation in F_ST across the genome, which could confound signals of local adaptation in genome scans. We performed realistic simulations of DNA sequences, using recombination maps from humans and sticklebacks, to investigate how variation in the intensity of background selection affects F_ST and other statistics of population differentiation in sexual, outcrossing species. We show that, in populations connected by gene flow, Weir and Cockerham's (1984; Evolution, 38 , 1358) estimator of F_ST is largely insensitive to locus‐to‐locus variation in the intensity of background selection. Unlike F_ST, however, d_XY is negatively correlated with background selection. Moreover, background selection does not greatly affect the false‐positive rate in F_ST outlier studies in populations connected by gene flow. Overall, our study indicates that background selection will not greatly interfere with finding the variants responsible for local adaptation.     

Signatures of adaptation and genetic structure among the mainland populations of Pinus radiata (D. Don) inferred from SNP loci

Insights into the relative contributions of locus specific and genome-wide effects on population genetic diversity can be gained through separation of their resulting genetic signals. Here we explore patterns of adaptive and neutral genetic diversity in the disjunct natural populations of Pinus radiata (D. Don) from mainland California. A first-generation common garden of 447 individuals revealed significant differentiation of wood phenotypes among populations (P _ST), possibly reflecting local adaptation in response to environment. We subsequently screened all trees for genetic diversity at 149 candidate gene single nucleotide polymorphism (SNP) loci for signatures of adaptation. Ten loci were identified as being possible targets of diversifying selection following F _ST outlier tests. Multivariate canonical correlation performed on a data set of 444 individuals identified significant covariance between environment, adaptive phenotypes and outlier SNP diversity, lending support to the case for local adaptation suggested from F _ST and P _ST tests. Covariation among discrete sets of outlier SNPs and adaptive phenotypes (inferred from multivariate loadings) with environment are supported by existing studies of candidate gene function and genotype–phenotype association. Canonical analyses failed to detect significant correlations between environment and 139 non-outlier SNP loci, which were applied to estimate neutral patterns of genetic differentiation among populations (F _ST 4.3 %). Using this data set, significant hierarchical structure was detected, indicating three populations on the mainland. The hierarchical relationships based on neutral SNP markers (and SSR) were in contrast with those inferred from putatively adaptive loci, potentially highlighting the independent action of selection and demography in shaping genetic structure in this species.     

Is local selection so widespread in river organisms? Fractal geometry of river networks leads to high bias in outlier detection

Identifying local adaptation is crucial in conservation biology to define ecotypes and establish management guidelines. Local adaptation is often inferred from the detection of loci showing a high differentiation between populations, the so‐called F_ST outliers. Methods of detection of loci under selection are reputed to be robust in most spatial population models. However, using simulations we showed that F_ST outlier tests provided a high rate of false‐positives (up to 60%) in fractal environments such as river networks. Surprisingly, the number of sampled demes was correlated with parameters of population genetic structure, such as the variance of F_STs, and hence strongly influenced the rate of outliers. This unappreciated property of river networks therefore needs to be accounted for in genetic studies on adaptation and conservation of river organisms.     

Evidence of genomic adaptation to climate in Eucalyptus microcarpa: Implications for adaptive potential to projected climate change

Understanding whether populations can adapt in situ or whether interventions are required is of key importance for biodiversity management under climate change. Landscape genomics is becoming an increasingly important and powerful tool for rapid assessments of climate adaptation, especially in long‐lived species such as trees. We investigated climate adaptation in Eucalyptus microcarpa using the DArTseq genomic approach. A combination of F_ST outlier and environmental association analyses were performed using >4200 genomewide single nucleotide polymorphisms (SNPs) from 26 populations spanning climate gradients in southeastern Australia. Eighty‐one SNPs were identified as putatively adaptive, based on significance in F_ST outlier tests and significant associations with one or more climate variables related to temperature (70/81), aridity (37/81) or precipitation (35/81). Adaptive SNPs were located on all 11 chromosomes, with no particular region associated with individual climate variables. Climate adaptation appeared to be characterized by subtle shifts in allele frequencies, with no consistent fixed differences identified. Based on these associations, we predict adaptation under projected changes in climate will include a suite of shifts in allele frequencies. Whether this can occur sufficiently rapidly through natural selection within populations, or would benefit from assisted gene migration, requires further evaluation. In some populations, the absence or predicted increases to near fixation of particular adaptive alleles hint at potential limits to adaptive capacity. Together, these results reinforce the importance of standing genetic variation at the geographic level for maintaining species’ evolutionary potential.     

Positive selection of IL‐33 in adaptive immunity of domestic Chinese goats

The identification of the candidate genes that play key role in phenotypic variation in livestock populations can provide new information about evolution and positive selection. IL‐33 (71954) (Interleukin) gene is associated with the increased nematode resistance in small ruminants; however, the role of IL‐33 for the genetic control of different diseases in Chinese goat breeds is poorly described in scientific literature. Therefore, the current investigation was performed for the better understanding of the molecular evolution and the positive selection of single‐nucleotide polymorphism in IL‐33 gene. Fixation Index (F_ST)‐based method was used for the outlier loci determination and found that IL‐33 was present in outlier area with the provisional combined allocation of mean heterozygosity and F_ST. Positively selected IL‐33 gene was significantly, that is, p(Simul F_ST < sample F_ST = 0.98*) present in corresponding positive selection area. Hence, our study provided novel information about the nucleotide variations in IL‐33 gene and found to be nonsynonymous which may helpful for the genetic control of diseases by enhancing the immune system in local Chinese goat breeds as well as in other analyzed vertebrate species.     

Differential introgression from a sister species explains high FST outlier loci within a mussel species

Scanning genomes for loci with high levels of population differentiation has become a standard of population genetics. F_ST outlier loci are most often interpreted as signatures of local selection, but outliers might arise for many other reasons too often left unexplored. Here, we tried to identify further the history and genetic basis underlying strong differentiation at F_ST outlier loci in a marine mussel. A genome scan of genetic differentiation has been conducted between Atlantic and Mediterranean populations of Mytilus galloprovincialis. The differentiation was low overall (F_ST = 0.03), but seven loci (2%) were strong F_ST outliers. We then analysed DNA sequence polymorphism at two outlier loci. The genetic structure proved to be the consequence of differential introgression of alleles from the sister‐hybridizing species Mytilus edulis. Surprisingly, the Mediterranean population was the most introgressed at these two loci, although the contact zone between the two species is nowadays localized along the Atlantic coasts of France and the British Isles. A historical contact between M. edulis and Mediterranean M. galloprovincialis should have happened during glacial periods. It proved difficult to disentangle two hypotheses: (i) introgression was adaptive, implying edulis alleles have been favoured in Mediterranean populations, or (ii) the genetic architecture of the barrier to edulis gene flow is different between the two M. galloprovincialis backgrounds. Five of the seven outliers between M. galloprovincialis populations were also outliers between M. edulis and Atlantic M. galloprovincialis, which would support the latter hypothesis. Differential introgression across semi‐permeable barriers to gene flow is a neglected scenario to interpret outlying loci that may prove more widespread than anticipated.     

Are QST–FST comparisons for natural populations meaningful?

Comparisons between putatively neutral genetic differentiation amongst populations, F_ST, and quantitative genetic variation, Q_ST, are increasingly being used to test for natural selection. However, we find that approximately half of the comparisons that use only data from wild populations confound phenotypic and genetic variation. We urge the use of a clear distinction between narrow‐sense Q_ST, which can be meaningfully compared with F_ST, and phenotypic divergence measured between populations, P_ST, which is inadequate for comparisons in the wild. We also point out that an unbiased estimate of Q_ST can be found using the so‐called ‘animal model’ of quantitative genetics.     

Intrapopulation genomics in a model mutualist: Population structure and candidate symbiosis genes under selection in Medicago truncatula

Bottom‐up evolutionary approaches, including geographically explicit population genomic analyses, have the power to reveal the mechanistic basis of adaptation. Here, we conduct a population genomic analysis in the model legume, Medicago truncatula, to characterize population genetic structure and identify symbiosis‐related genes showing evidence of spatially variable selection. Using RAD‐seq, we generated over 26,000 SNPs from 191 accessions from within three regions of the native range in Europe. Results from STRUCTURE analysis identify five distinct genetic clusters with divisions that separate east and west regions in the Mediterranean basin. Much of the genetic variation is maintained within sampling sites, and there is evidence for isolation by distance. Extensive linkage disequilibrium was identified, particularly within populations. We conducted genetic outlier analysis with F_ST‐based genome scans and a Bayesian modeling approach (PCAdapt). There were 70 core outlier loci shared between these distinct methods with one clear candidate symbiosis related gene, DMI1. This work sets that stage for functional experiments to determine the important phenotypes that selection has acted upon and complementary efforts in rhizobium populations.     

Population differentiation determined from putative neutral and divergent adaptive genetic markers in Eulachon (Thaleichthys pacificus,Osmeridae), an anadromous Pacific smelt

Twelve eulachon (Thaleichthys pacificus, Osmeridae) populations ranging from Cook Inlet, Alaska and along the west coast of North America to the Columbia River were examined by restriction‐site‐associated DNA (RAD) sequencing to elucidate patterns of neutral and adaptive variation in this high geneflow species. A total of 4104 single‐nucleotide polymorphisms (SNPs) were discovered across the genome, with 193 putatively adaptive SNPs as determined by F_ST outlier tests. Estimates of population structure in eulachon with the putatively adaptive SNPs were similar, but provided greater resolution of stocks compared with a putatively neutral panel of 3911 SNPs or previous estimates with 14 microsatellites. A cline of increasing measures of genetic diversity from south to north was found in the adaptive panel, but not in the neutral markers (SNPs or microsatellites). This may indicate divergent selective pressures in differing freshwater and marine environments between regional eulachon populations and that these adaptive diversity patterns not seen with neutral markers could be a consideration when determining genetic boundaries for conservation purposes. Estimates of effective population size (N_e) were similar with the neutral SNP panel and microsatellites and may be utilized to monitor population status for eulachon where census sizes are difficult to obtain. Greater differentiation with the panel of putatively adaptive SNPs provided higher individual assignment accuracy compared to the neutral panel or microsatellites for stock identification purposes. This study presents the first SNPs that have been developed for eulachon, and analyses with these markers highlighted the importance of integrating genome‐wide neutral and adaptive genetic variation for the applications of conservation and management.     

Genetic divergence and isolation by thermal environment in geothermal populations of an aquatic invertebrate

Temperature is one of the most influential forces of natural selection impacting all biological levels. In the face of increasing global temperatures, studies over small geographic scales allowing investigations on the effects of gene flow are of great value for understanding thermal adaptation. Here, we investigated genetic population structure in the freshwater gastropod Radix balthica originating from contrasting thermal habitats in three areas of geothermal activity in Iceland. Snails from 32 sites were genotyped at 208 AFLP loci. Five AFLPs were identified as putatively under divergent selection in Lake Mývatn, a geothermal lake with an almost 20 °C difference in mean temperature across a distance of a few kilometres. In four of these loci, variation across all study populations was correlated with temperature. We found significant population structure in neutral markers both within and between the areas. Cluster analysis using neutral markers classified the sites mainly by geography, whereas analyses using markers under selection differentiated the sites based on temperature. Isolation by distance was stronger in the neutral than in the outlier loci. Pairwise differences based on outlier F_ST were significantly correlated with temperature at different spatial scales, even after correcting for geographic distance or neutral pairwise F_ST differences. In general, genetic variation decreased with increasing environmental temperature, possibly suggesting that natural selection had reduced the genetic diversity in the warm origin sites. Our results emphasize the influence of environmental temperature on the genetic structure of populations and suggest local thermal adaptation in these geothermal habitats.     

Candidate gene SNP variation in floodplain populations of pedunculate oak (Quercus robur L.) near the species' southern range margin: Weak differentiation yet distinct associations with water availability

Populations residing near species' low‐latitude range margins (LLMs) often occur in warmer and drier environments than those in the core range. Thus, their genetic composition could be shaped by climatic drivers that differ from those occurring at higher latitudes, resulting in potentially adaptive variants of conservation value. Such variants could facilitate the adaptation of populations from other portions of the geographical range to similar future conditions anticipated under ongoing climate change. However, very few studies have assessed standing genetic variation at potentially adaptive loci in natural LLM populations. We investigated standing genetic variation at single nucleotide polymorphisms (SNPs) located within 117 candidate genes and its links to putative climatic selection pressures across 19 pedunculate oak (Quercus robur L.) populations distributed along a regional climatic gradient near the species' southern range margin in southeastern Europe. These populations are restricted to floodplain forests along large lowland rivers, whose hydric regime is undergoing significant shifts under modern rapid climate change. The populations showed very weak geographical structure, suggesting extensive genetic connectivity and gene flow or shared ancestry. We identified eight (6.2%) positive F_ST‐outlier loci, and genotype–environment association analyses revealed consistent associations between SNP allele frequencies and several climatic variables linked to water availability. A total of 61 associations involving 37 SNPs (28.5%) from 35 annotated genes provided important insights into putative functional mechanisms in our system. Our findings provide empirical support for the role of LLM populations as sources of potentially adaptive variation that could enhance species’ resilience to climate change‐related pressures.