Contrasting patterns of body shape and neutral genetic divergence in marine and lake populations of threespine sticklebacks

Comparisons of neutral marker and quantitative trait divergence can provide important insights into the relative roles of natural selection and neutral genetic drift in population differentiation. We investigated phenotypic and genetic differentiation among Fennoscandian threespine stickleback (Gasterosteus aculeatus) populations, and found that the highest degree of differentiation occurred between sea and freshwater habitats. Within habitats, morphological divergence was highest among the different freshwater populations. Pairwise phenotypic and neutral genetic distances among populations were positively correlated, suggesting that genetic drift may have contributed to the morphological differentiation among habitats. On the other hand, the degree of phenotypic differentiation (PST) clearly surpassed the neutral expectation set by FST, suggesting a predominant role for natural selection over genetic drift as an explanation for the observed differentiation. However, separate PST/FST comparisons by habitats revealed that body shape divergence between lake and marine populations, and even among marine populations, can be strongly influenced by natural selection. On the other hand, genetic drift can play an important role in the differentiation among lake populations.     

Does natural selection promote population divergence? A comparative analysis of population structure using amplified fragment length polymorphism markers and quantitative traits

Divergent natural selection is considered an important force in plant evolution leading to phenotypic differentiation between populations exploiting different environments. Extending an earlier greenhouse study of population differentiation in the selfing annual plant Senecio vulgaris, we estimated the degree of population divergence in several quantitative traits related to growth and life history and compared these estimates with those based on presumably neutral molecular markers (amplified fragment length polymorphisms; AFLPs). This approach allowed us to disentangle the effects of divergent selection from that of other evolutionary forces (e.g. genetic drift). Five populations were examined from each of two habitat types (ruderal and agricultural habitats). We found a high proportion of total genetic variance to be among populations, both for AFLP markers (phiST = 0.49) and for quantitative traits (range of QST: 0.26-0.77). There was a strong correlation between molecular and quantitative genetic differentiation between pairs of populations (Mantel's r = 0.59). However, estimates of population differentiation in several quantitative traits exceeded the neutral expectation (estimated from AFLP data), suggesting that divergent selection contributed to phenotypic differentiation, especially between populations from ruderal and agricultural habitats. Estimates of within-population variation in AFLP markers and quantitative genetic were poorly correlated, indicating that molecular marker data may be of limited value to predict the evolutionary potential of populations of S. vulgaris.     

Adaptive differentiation in floral traits in the presence of high gene flow in scarlet gilia (Ipomopsis aggregata)

Plant–pollinator interactions are thought to be major drivers of floral trait diversity. However, the relative importance of divergent pollinator‐mediated selection vs. neutral processes in floral character evolution has rarely been explored. We tested for adaptive floral trait evolution by comparing differentiation at neutral genetic loci to differentiation at quantitative floral traits in a putative Ipomopsis aggregata hybrid zone. Typical I. aggregata subsp. candida displays slender white tubular flowers that are typical of flowers pollinated by hawkmoths, and subsp. collina displays robust red tubular flowers typical of flowers pollinated by hummingbirds; yet, hybrid flower morphs are abundant across the East Slope of the Colorado Rockies. We estimated genetic differentiation (F_ST) for nuclear and chloroplast microsatellite loci and used a half‐sib design to calculate quantitative trait divergence (Q_ST) from collection sites across the morphological hybrid zone. We found little evidence for population structure and estimated mean F_ST to be 0.032. Q_ST values for several floral traits including corolla tube length and width, colour, and nectar volume were large and significantly greater than mean F_ST. We performed multivariate comparisons of neutral loci to genetic correlations within and between populations and found a strong signal for divergent selection, suggesting that specific combinations of floral display and reward traits may be the targets of selection. Our results show little support for historical subspecies categories, yet floral traits are more diverged than expected due to drift alone. Non‐neutral divergence for multivariate quantitative traits suggests that selection by pollinators is maintaining a correlation between display and reward traits.     

Selection on life‐history traits and genetic population divergence in rotifers

A combination of founder effects and local adaptation – the Monopolization hypothesis – has been proposed to reconcile the strong population differentiation of zooplankton dwelling in ponds and lakes and their high dispersal abilities. The role genetic drift plays in genetic differentiation of zooplankton is well documented, but the impact of natural selection has received less attention. Here, we compare differentiation in neutral genetic markers (F_ST) and in quantitative traits (Q_ST) in six natural populations of the rotifer Brachionus plicatilis to assess the importance of natural selection in explaining genetic differentiation of life‐history traits. Five life‐history traits were measured in four temperature × salinity combinations in common‐garden experiments. Population differentiation for neutral genetic markers – 11 microsatellite loci – was very high (F_ST = 0.482). Differentiation in life‐history traits was higher in traits related to sexual reproduction than in those related to asexual reproduction. Q_ST values for diapausing egg production (a trait related to sexual reproduction) were higher than their corresponding F_ST in some pairs of populations. Our results indicate the importance of divergent natural selection in these populations and suggest local adaptation to the unpredictability of B. plicatilis habitats.     

Genetic isolation of fragmented populations is exacerbated by drift and selection

Reduced genetic variation at marker loci in small populations has been well documented, whereas the relationship between quantitative genetic variation and population size has attracted little empirical investigation. Here we demonstrate that both neutral and quantitative genetic variation are reduced in small populations of a fragmented plant metapopulation, and that both drift and selective change are enhanced in small populations. Measures of neutral genetic differentiation (F(ST)) and quantitative genetic differentiation (Q(ST)) in two traits were higher among small demes, and Q(ST) between small populations exceeded that expected from drift alone. This suggests that fragmented populations experience both enhanced genetic drift and divergent selection on phenotypic traits, and that drift affects variation in both neutral markers and quantitative traits. These results highlight the need to integrate natural selection into conservation genetic theory, and suggests that small populations may represent reservoirs of genetic variation adaptive within a wide range of environments.     

Divergent selection along climatic gradients in a rare central European endemic species,Saxifraga sponhemica

Background and Aims The effects of habitat fragmentation on quantitative genetic variation in plant populations are still poorly known. Saxifraga sponhemica is a rare endemic of Central Europe with a disjunct distribution, and a stable and specialized habitat of treeless screes and cliffs. This study therefore used S. sponhemica as a model species to compare quantitative and molecular variation in order to explore (1) the relative importance of drift and selection in shaping the distribution of quantitative genetic variation along climatic gradients; (2) the relationship between plant fitness, quantitative genetic variation, molecular genetic variation and population size; and (3) the relationship between the differentiation of a trait among populations and its evolvability.Methods Genetic variation within and among 22 populations from the whole distribution area of S. sponhemica was studied using RAPD (random amplified polymorphic DNA) markers, and climatic variables were obtained for each site. Seeds were collected from each population and germinated, and seedlings were transplanted into a common garden for determination of variation in plant traits.Key Results In contrast to previous results from rare plant species, strong evidence was found for divergent selection. Most population trait means of S. sponhemica were significantly related to climate gradients, indicating adaptation. Quantitative genetic differentiation increased with geographical distance, even when neutral molecular divergence was controlled for, and Q_ST exceeded F_ST for some traits. The evolvability of traits was negatively correlated with the degree of differentiation among populations (Q_ST), i.e. traits under strong selection showed little genetic variation within populations. The evolutionary potential of a population was not related to its size, the performance of the population or its neutral genetic diversity. However, performance in the common garden was lower for plants from populations with reduced molecular genetic variation, suggesting inbreeding depression due to genetic erosion.Conclusions The findings suggest that studies of molecular and quantitative genetic variation may provide complementary insights important for the conservation of rare species. The strong differentiation of quantitative traits among populations shows that selection can be an important force for structuring variation in evolutionarily important traits even for rare endemic species restricted to very specific habitats.     

Patterns of differentiation in a colour polymorphism and in neutral markers reveal rapid genetic changes in natural damselfly populations

The existence and mode of selection operating on heritable adaptive traits can be inferred by comparing population differentiation in neutral genetic variation between populations (often using F(ST) values) with the corresponding estimates for adaptive traits. Such comparisons indicate if selection acts in a diversifying way between populations, in which case differentiation in selected traits is expected to exceed differentiation in neutral markers [F(ST )(selected) > F(ST )(neutral)], or if negative frequency-dependent selection maintains genetic polymorphisms and pulls populations towards a common stable equilibrium [F(ST) (selected) < F(ST) (neutral)]. Here, we compared F(ST) values for putatively neutral data (obtained using amplified fragment length polymorphism) with estimates of differentiation in morph frequencies in the colour-polymorphic damselfly Ischnura elegans. We found that in the first year (2000), population differentiation in morph frequencies was significantly greater than differentiation in neutral loci, while in 2002 (only 2 years and 2 generations later), population differentiation in morph frequencies had decreased to a level significantly lower than differentiation in neutral loci. Genetic drift as an explanation for population differentiation in morph frequencies could thus be rejected in both years. These results indicate that the type and/or strength of selection on morph frequencies in this system can change substantially between years. We suggest that an approach to a common equilibrium morph frequency across all populations, driven by negative frequency-dependent selection, is the cause of these temporal changes. We conclude that inferences about selection obtained by comparing F(ST) values from neutral and adaptive genetic variation are most useful when spatial and temporal data are available from several populations and time points and when such information is combined with other ecological sources of data.     

Does habitat fragmentation reduce fitness and adaptability? A case study of the common frog (Rana temporaria)

Studies examining the effects of anthropogenic habitat fragmentation on both neutral and adaptive genetic variability are still scarce. We compared tadpole fitness-related traits (viz. survival probability and body size) among populations of the common frog (Rana temporaria) from fragmented (F) and continuous (C) habitats that differed significantly in population sizes (C > F) and genetic diversity (C > F) in neutral genetic markers. Using data from common garden experiments, we found a significant positive relationship between the mean values of the fitness related traits and the amount of microsatellite variation in a given population. While genetic differentiation in neutral marker loci (F(ST)) tended to be more pronounced in the fragmented than in the continuous habitat, genetic differentiation in quantitative traits (Q(ST)) exceeded that in neutral marker traits in the continuous habitat (i.e. Q(ST) > F(ST)), but not in the fragmented habitat (i.e. Q(ST) approximately F(ST)). These results suggest that the impact of random genetic drift relative to natural selection was higher in the fragmented landscape where populations were small, and had lower genetic diversity and fitness as compared to populations in the more continuous landscape. The findings highlight the potential importance of habitat fragmentation in impairing future adaptive potential of natural populations.     

The effects of selection,drift and genetic variation on life‐history trait divergence among insular populations of natterjack toad,Bufo calamita

Although loss of genetic variation is frequently assumed to be associated with loss of adaptive potential, only few studies have examined adaptation in populations with little genetic variation. On the Swedish west coast, the northern fringe populations of the natterjack toad Bufo calamita inhabit an atypical habitat consisting of offshore rock islands. There are strong among‐population differences in the amount of neutral genetic variation, making this system suitable for studies on mechanisms of trait divergence along a gradient of within‐population genetic variation. In this study, we examined the mechanisms of population divergence using Q_ST–F_ST comparisons and correlations between quantitative and neutral genetic variation. Our results suggest drift or weak stabilizing selection across the six populations included in this study, as indicated by low Q_ST–F_ST values, lack of significant population × temperature interactions and lack of significant differences among the islands in breeding pond size. The six populations included in this study differed in both neutral and quantitative genetic variation. Also, the correlations between neutral and quantitative genetic variation tended to be positive, however, the relatively small number of populations prevents any strong conclusions based on these correlations. Contrary to the majority of Q_ST–F_ST comparisons, our results suggest drift or weak stabilizing selection across the examined populations. Furthermore, the low heritability of fitness‐related traits may limit evolutionary responses in some of the populations.     

Comparison of genetic differentiation at marker loci and quantitative traits

The comparison of the degree of differentiation in neutral marker loci and genes coding quantitative traits with standardized and equivalent measures of genetic differentiation (F_ST and Q_ST, respectively) can provide insights into two important but seldom explored questions in evolutionary genetics: (i) what is the relative importance of random genetic drift and directional natural selection as causes of population differentiation in quantitative traits, and (ii) does the degree of divergence in neutral marker loci predict the degree of divergence in genes coding quantitative traits? Examination of data from 18 independent studies of plants and animals using both standard statistical and meta‐analytical methods revealed a number of interesting points. First, the degree of differentiation in quantitative traits (Q_ST) typically exceeds that observed in neutral marker genes (F_ST), suggesting a prominent role for natural selection in accounting for patterns of quantitative trait differentiation among contemporary populations. Second, the F_ST – Q_ST difference is more pronounced for allozyme markers and morphological traits, than for other kinds of molecular markers and life‐history traits. Third, very few studies reveal situations were Q_ST < F_ST, suggesting that selection pressures, and hence optimal phenotypes, in different populations of the same species are unlikely to be often similar. Fourth, there is a strong correlation between Q_ST and F_ST indices across the different studies for allozyme (r=0.81), microsatellite (r=0.87) and combined (r=0.75) marker data, suggesting that the degree of genetic differentiation in neutral marker loci is closely predictive of the degree of differentiation in loci coding quantitative traits. However, these interpretations are subject to a number of assumptions about the data and methods used to derive the estimates of population differentiation in the two sets of traits.     

Local selection modifies phenotypic divergence among Rana temporaria populations in the presence of gene flow

In ectotherms, variation in life history traits among populations is common and suggests local adaptation. However, geographic variation itself is not a proof for local adaptation, as genetic drift and gene flow may also shape patterns of quantitative variation. We studied local and regional variation in means and phenotypic plasticity of larval life history traits in the common frog Rana temporaria using six populations from central Sweden, breeding in either open‐canopy or partially closed‐canopy ponds. To separate local adaptation from genetic drift, we compared differentiation in quantitative genetic traits (Q_ST) obtained from a common garden experiment with differentiation in presumably neutral microsatellite markers (F_ST). We found that R. temporaria populations differ in means and plasticities of life history traits in different temperatures at local, and in F_ST at regional scale. Comparisons of differentiation in quantitative traits and in molecular markers suggested that natural selection was responsible for the divergence in growth and development rates as well as in temperature‐induced plasticity, indicating local adaptation. However, at low temperature, the role of genetic drift could not be separated from selection. Phenotypes were correlated with forest canopy closure, but not with geographical or genetic distance. These results indicate that local adaptation can evolve in the presence of ongoing gene flow among the populations, and that natural selection is strong in this system.     

Clonal erosion and genetic drift in cyclical parthenogens – the interplay between neutral and selective processes

The occurrence of alternating phases of clonal and sexual reproduction may strongly impact the interplay between neutral and selective genetic variation in populations. Using a physiologically structured model of the life history of Daphnia, we investigated to what extent clonal erosion associated with selection during the clonal phase affects the genetic structure as observed by neutral markers. Incorporating conservative levels of quantitative genetic variation at 11 physiological and life history traits induces strong clonal erosion, reducing clonal diversity (CD) near the end of the simulations (1000 days) to a level between 1 and 5, even in habitats with high initial CD (10⁸ clones). This strong clonal erosion caused by selection can result in reduced genetic diversity, significant excess of heterozygotes and significant genetic differentiation between populations as observed by neutral markers. Our results indicate that, especially in relatively small habitats, clonal selection may strongly impact the genetic structure and may contribute to the often observed high level of neutral genetic differentiation among natural populations of cyclical parthenogens.     

Recent local adaptation of sockeye salmon to glacial spawning habitats

Salmonids spawn in highly diverse habitats, exhibit strong genetic population structuring, and can quickly colonize newly created habitats with few founders. Spawning traits often differ among populations, but it is largely unknown if these differences are adaptive or due to genetic drift. To test if sockeye salmon (Oncorhynchus nerka) populations are adapted to glacial, beach, and tributary spawning habitats, we examined variation in heritable phenotypic traits associated with spawning in 13 populations of wild sockeye salmon in Lake Clark, Alaska. These populations were commonly founded between 100 and 400 hundred sockeye salmon generations ago and exhibit low genetic divergence at 11 microsatellite loci (F _ST < 0.024) that is uncorrelated with spawning habitat type. We found that mean P _ST (phenotypic divergence among populations) exceeded neutral F _ST for most phenotypic traits measured, indicating that phenotypic differences among populations could not be explained by genetic drift alone. Phenotypic divergence among populations was associated with spawning habitat differences, but not with neutral genetic divergence. For example, female body color was lighter and egg color was darker in glacial than non-glacial habitats. This may be due to reduced sexual selection for red spawning color in glacial habitats and an apparent trade-off in carotenoid allocation to body and egg color in females. Phenotypic plasticity is an unlikely source of phenotypic differences because Lake Clark sockeye salmon spend nearly all their lives in a common environment. Our data suggest that Lake Clark sockeye salmon populations are adapted to spawning in glacial, beach and tributary habitats and provide the first evidence of a glacial spawning ecotype in salmonids. Glacial spawning habitats are often young (i.e., <200 years old) and ephemeral. Thus, local adaptation of sockeye salmon to glacial habitats appears to have occurred recently.     

Contrasting Levels of Variation in Neutral and Quantitative Genetic Loci on Island Populations of Moor Frogs (<Emphasis Type="Italic">Rana arvalis</Emphasis>)

Reduced levels of genetic variability and a prominent differentiation in both neutral marker genes and phenotypic traits are typical for many island populations as compared to their mainland conspecifics. However, whether genetic diversity in neutral marker genes reflects genetic variability in quantitative traits, and thus, their evolutionary potential, remains typically unclear. Moreover, the phenotypic differentiation on islands could be attributable to phenotypic plasticity, selection or drift; something which seldom has been tested. Using eight polymorphic microsatellite loci and quantitative genetic breeding experiments we conducted a detailed comparison on genetic variability and differentiation between Nordic islands (viz. Gotland, Öland and Læsø) and neighbouring mainland populations of moor frogs (Rana arvalis). As expected, the neutral variation was generally lower in island than in mainland populations. But as opposed to this, higher levels of additive genetic variation (V _A) in body size and tibia length were found on the island of Gotland as compared to the mainland population. When comparing the differentiation seen in neutral marker genes (F _ST) with the differentiation in genes coding quantitative traits (Q _ST) two different evolutionary scenarios were found: while selection might explain a smaller size of moor frogs on Gotland, the differentiation seen in tibia length could be explained by genetic drift. These results highlight the limited utility of microsatellite loci alone in inferring the causes behind an observed phenotypic differentiation, or in predicting the amount of genetic variation in ecologically important quantitative traits.     

Potential adaptive divergence between subspecies and populations of snapdragon plants inferred from QST–FST comparisons

Phenotypic divergence among natural populations can be explained by natural selection or by neutral processes such as drift. Many examples in the literature compare putatively neutral (F_ST) and quantitative genetic (Q_ST) differentiation in multiple populations to assess their evolutionary signature and identify candidate traits involved with local adaptation. Investigating these signatures in closely related or recently diversified species has the potential to shed light on the divergence processes acting at the interspecific level. Here, we conducted this comparison in two subspecies of snapdragon plants (eight populations of Antirrhinum majus pseudomajus and five populations of A. m. striatum) in a common garden experiment. We also tested whether altitude was involved with population phenotypic divergence. Our results identified candidate phenological and morphological traits involved with local adaptation. Most of these traits were identified in one subspecies but not the other. Phenotypic divergence increased with altitude for a few biomass‐related traits, but only in A. m. striatum. These traits therefore potentially reflect A. m. striatum adaptation to altitude. Our findings imply that adaptive processes potentially differ at the scale of A. majus subspecies.     

Molecular and quantitative trait variation within and among small fragmented populations of the endangered plant species Psilopeganum sinense

Background and Aims

Natural selection and genetic drift are important evolutionary forces in determining genetic and phenotypic differentiation in plant populations. The extent to which these two distinct evolutionary forces affect locally adaptive quantitative traits has been well studied in common plant and animal species. However, we know less about how quantitative traits respond to selection pressures and drift in endangered species that have small population sizes and fragmented distributions. To address this question, this study assessed the relative strengths of selection and genetic drift in shaping population differentiation of phenotypic traits in Psilopeganum sinense, a naturally rare and recently endangered plant species.

Methods

Population differentiation at five quantitative traits (Q_ST) obtained from a common garden experiment was compared with differentiation at putatively neutral microsatellite markers (F_ST) in seven populations of P. sinense. Q_ST estimates were derived using a Bayesian hierarchical variance component method.

Key Results

Trait-specific Q_ST values were equal to or lower than F_ST. Neutral genetic diversity was not correlated with quantitative genetic variation within the populations of P. sinense.

Conclusions

Despite the prevalent empirical evidence for Q_ST > F_ST, the results instead suggest a definitive role of stabilizing selection and drift leading to phenotypic differentiation among small populations. Three traits exhibited a significantly lower Q_ST relative to F_ST, suggesting that populations of P. sinense might have experienced stabilizing selection for the same optimal phenotypes despite large geographical distances between populations and habitat fragmentation. For the other two traits, Q_ST estimates were of the same magnitude as F_ST, indicating that divergence in these traits could have been achieved by genetic drift alone. The lack of correlation between molecular marker and quantitative genetic variation suggests that sophisticated considerations are required for the inference of conservation measures of P. sinense from neutral genetic markers.     

Selection and drift influence genetic differentiation of insular Canada lynx (Lynx canadensis) on Newfoundland and Cape Breton Island

Island populations have long been important for understanding the dynamics and mechanisms of evolution in natural systems. While genetic drift is often strong on islands due to founder events and population bottlenecks, the strength of selection can also be strong enough to counteract the effects of drift. Here, we used several analyses to identify the roles of genetic drift and selection on genetic differentiation and diversity of Canada lynx (Lynx canadensis) across eastern Canada, including the islands of Cape Breton and Newfoundland. Specifically, we assessed whether we could identify a genetic component to the observed morphological differentiation that has been reported across insular and mainland lynx. We used a dinucleotide repeat within the promoter region of a functional gene that has been linked to mammalian body size, insulin‐like growth factor‐1 (IGF‐1). We found high genetic differentiation at neutral molecular markers but convergence of allele frequencies at the IGF‐1 locus. Thus, we showed that while genetic drift has influenced the observed genetic structure of lynx at neutral molecular markers, natural selection has also played a role in the observed patterns of genetic diversity at the IGF‐1 locus of insular lynx.     

LOCAL ADAPTATION MAINTAINS CLINAL VARIATION IN MELANIN‐BASED COLORATION OF EUROPEAN BARN OWLS (TYTO ALBA)

Ecological parameters vary in space, and the resulting heterogeneity of selective forces can drive adaptive population divergence. Clinal variation represents a classical model to study the interplay of gene flow and selection in the dynamics of this local adaptation process. Although geographic variation in phenotypic traits in discrete populations could be remainders of past adaptation, maintenance of adaptive clinal variation requires recurrent selection. Clinal variation in genetically determined traits is generally attributed to adaptation of different genotypes to local conditions along an environmental gradient, although it can as well arise from neutral processes. Here, we investigated whether selection accounts for the strong clinal variation observed in a highly heritable pheomelanin‐based color trait in the European barn owl by comparing spatial differentiation of color and of neutral genes among populations. Barn owl's coloration varies continuously from white in southwestern Europe to reddish‐brown in northeastern Europe. A very low differentiation at neutral genetic markers suggests that substantial gene flow occurs among populations. The persistence of pronounced color differentiation despite this strong gene flow is consistent with the hypothesis that selection is the primary force maintaining color variation among European populations. Therefore, the color cline is most likely the result of local adaptation.     

driftsel: an R package for detecting signals of natural selection in quantitative traits

Approaches and tools to differentiate between natural selection and genetic drift as causes of population differentiation are of frequent demand in evolutionary biology. Based on the approach of Ovaskainen et al. (2011), we have developed an R package (driftsel ) that can be used to differentiate between stabilizing selection, diversifying selection and random genetic drift as causes of population differentiation in quantitative traits when neutral marker and quantitative genetic data are available. Apart from illustrating the use of this method and the interpretation of results using simulated data, we apply the package on data from three‐spined sticklebacks (Gasterosteus aculeatus) to highlight its virtues. driftsel can also be used to perform usual quantitative genetic analyses in common‐garden study designs.     

Genetic analysis of differentiation among breeding ponds reveals a candidate gene for local adaptation in Rana arvalis

One of the main questions in evolutionary and conservation biology is how geographical and environmental features of the landscape shape neutral and adaptive genetic variation in natural populations. The identification of genomic polymorphisms that account for adaptive variation can aid in finding candidate loci for local adaptation. Consequently, a comparison of spatial patterns in neutral markers and loci under selection may help disentangle the effects of gene flow, genetic drift and selection at the landscape scale. Many amphibians breed in wetlands, which differ in environmental conditions and in the degree of isolation, enhancing the potential for local adaptation. We used microsatellite markers to measure genetic differentiation among 17 local populations of Rana arvalis breeding in a network of wetlands. We found that locus RC08604 deviated from neutral expectations, suggesting that it is a good candidate for directional selection. We used a genetic network analysis to show that the allele distribution in this locus is correlated with habitat characteristics, whereas this was not the case at neutral markers that displayed a different allele distribution and population network in the study area. The graph approach illustrated the genomic heterogeneity (neutral loci vs. the candidate locus for directional selection) of gene exchange and genetic divergence among populations under directional selection. Limited gene flow between wetlands was only observed at the candidate genomic region under directional selection. RC08604 is partially located inside an up‐regulated thyroid‐hormone receptor (TRβ) gene coordinating the expression of other genes during metamorphosis and appears to be linked with variation in larval life‐history traits found among R. arvalis populations. We suggest that directional selection on genes coding larval life‐history traits is strong enough to maintain the divergence in these genomic regions, reducing the effective recombination of locally adapted alleles but not in other regions of the genome. Integrating this knowledge into conservation plans at the landscape scale will improve the design of management strategies to preserve adaptive genetic diversity in wetland networks.