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.     

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.     

Climate variables explain neutral and adaptive variation within salmonid metapopulations: the importance of replication in landscape genetics

Understanding how environmental variation influences population genetic structure is important for conservation management because it can reveal how human stressors influence population connectivity, genetic diversity and persistence. We used riverscape genetics modelling to assess whether climatic and habitat variables were related to neutral and adaptive patterns of genetic differentiation (population‐specific and pairwise F_ST) within five metapopulations (79 populations, 4583 individuals) of steelhead trout (Oncorhynchus mykiss) in the Columbia River Basin, USA. Using 151 putatively neutral and 29 candidate adaptive SNP loci, we found that climate‐related variables (winter precipitation, summer maximum temperature, winter highest 5% flow events and summer mean flow) best explained neutral and adaptive patterns of genetic differentiation within metapopulations, suggesting that climatic variation likely influences both demography (neutral variation) and local adaptation (adaptive variation). However, we did not observe consistent relationships between climate variables and F_ST across all metapopulations, underscoring the need for replication when extrapolating results from one scale to another (e.g. basin‐wide to the metapopulation scale). Sensitivity analysis (leave‐one‐population‐out) revealed consistent relationships between climate variables and F_ST within three metapopulations; however, these patterns were not consistent in two metapopulations likely due to small sample sizes (N = 10). These results provide correlative evidence that climatic variation has shaped the genetic structure of steelhead populations and highlight the need for replication and sensitivity analyses in land and riverscape genetics.     

Evidence for adaptive phenotypic differentiation in Baltic Sea sticklebacks

The evidence for adaptive phenotypic differentiation in mobile marine species remains scarce, partly due to the difficulty of obtaining quantitative genetic data to demonstrate the genetic basis of the observed phenotypic differentiation. Using a combination of phenotypic and molecular genetic approaches, we elucidated the relative roles of natural selection and genetic drift in explaining lateral plate number differentiation in threespine sticklebacks (Gasterosteus aculeatus) across the entire Baltic Sea basin (approximately 392 000 km²). We found that phenotypic differentiation (P_ST = 0.213) in plate number exceeded that in neutral markers (F_ST = 0.008), suggesting an adaptive basis for the observed differentiation. Because a close correspondence was found between plate phenotype and genotype at a quantitative trait loci (QTL; STN381) tightly linked to the gene (Ectodysplasin) underlying plate variation, the evidence for adaptive differentiation was confirmed by comparison of F_ST at the QTL (F_STQ = 0.089) with F_ST at neutral marker loci. Hence, the results provide a comprehensive demonstration of adaptive phenotypic differentiation in a high‐gene‐flow marine environment with direct, rather than inferred, verification for the genetic basis of this differentiation. In general, the results illustrate the utility of P_ST–F_ST–F_STQ comparisons in uncovering footprints of natural selection and evolution and add to the growing evidence for adaptive genetic differentiation in high‐gene‐flow marine environments, including that of the relatively young Baltic Sea.     

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.     

Population size is weakly related to quantitative genetic variation and trait differentiation in a stream fish

How population size influences quantitative genetic variation and differentiation among natural, fragmented populations remains unresolved. Small, isolated populations might occupy poor quality habitats and lose genetic variation more rapidly due to genetic drift than large populations. Genetic drift might furthermore overcome selection as population size decreases. Collectively, this might result in directional changes in additive genetic variation (V_A) and trait differentiation (Q_ST) from small to large population size. Alternatively, small populations might exhibit larger variation in V_A and Q_ST if habitat fragmentation increases variability in habitat types. We explored these alternatives by investigating V_A and Q_ST using nine fragmented populations of brook trout varying 50‐fold in census size N (179–8416) and 10‐fold in effective number of breeders, N_b (18–135). Across 15 traits, no evidence was found for consistent differences in V_A and Q_ST with population size and almost no evidence for increased variability of V_A or Q_ST estimates at small population size. This suggests that (i) small populations of some species may retain adaptive potential according to commonly adopted quantitative genetic measures and (ii) populations of varying sizes experience a variety of environmental conditions in nature, however extremely large studies are likely required before any firm conclusions can be made.     

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.     

A meta-analysis of the effects of plant traits and geographical scale on the magnitude of adaptive differentiation as measured by the difference between QST and FST

The impact of directional selection on specific trait types in plant species, and how a species’ life history mediates this response to selection remains understudied. Discovering such interactions is however crucial for understanding the interplay between ecological and genetic processes underlying local adaptation in plants, and to evaluate a species’ evolutionary potential with respect to changing environments. Furthermore, it remains unclear whether the degree of adaptive differentiation generally increases with the geographical distance between plant populations. Here, we present a weighted mixed model based meta-analysis aimed at unraveling the potential interactions between plant trait types, life history characteristics and Q_ST–F_ST comparisons, and assessing the effect of geographical scale on population differentiation. Based on 51 studies we found that Q_ST values exceeded their corresponding F_ST values in 71.74 % out of 401 cases. Furthermore, different trait types were found to be differently susceptible to natural selection and the magnitude of Q_ST–F_ST comparisons was mediated by a plant species’ life span. These findings may be closely related to the genetic architectures of trait types and life histories, with the proportion of large-effect genes likely shaping the response to natural selection. Q_ST–F_ST values also increased with increasing distance between populations, pinpointing the combined effects of environmental differentiation and isolation by distance on the magnitude of population divergence. Finally, our model showed an inverse relationship between F_ST and Q_ST–F_ST values, presumably resulting from isolation by distance, the exchange of advantageous alleles, or genetic correlations among traits.     

Phenotypic landscapes: phenological patterns in wild and cultivated barley

Phenotypic variation in natural populations is the outcome of the joint effects of environmentally induced adaptations and neutral processes on the genetic architecture of quantitative traits. In this study, we examined the role of adaptation in shaping wild barley phenotypic variation along different environmental gradients. Detailed phenotyping of 164 wild barley (Hordeum spontaneum) accessions from Israel (of the Barley1K collection) and 18 cultivated barley (H. vulgare) varieties was conducted in common garden field trials. Cluster analysis based on phenotypic data indicated that wild barley in this region can be differentiated into three ecotypes in accordance with their ecogeographical distribution: north, coast and desert. Population differentiation (Q_ST) for each trait was estimated using a hierarchical Bayesian model and compared to neutral differentiation (F_ST) based on 42 microsatellite markers. This analysis indicated that the three clusters diverged in morphological but not in reproductive characteristics. To address the issue of phenotypic variation along environmental gradients, climatic and soil gradients were compared with each of the measured traits given the geographical distance between sampling sites using a partial Mantel test. Flowering time and plant growth were found to be differentially correlated with climatic and soil characteristic gradients, respectively. The H. vulgare varieties were superior to the H. spontaneum accessions in yield components, yet resembled the Mediterranean types in vegetative characteristics and flowering time, which may indicate the geographical origin of domesticated barley.     

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.     

Phenotypic divergence of the common toad (Bufo bufo) along an altitudinal gradient: evidence for local adaptation

Variation in the environment can induce different patterns of genetic and phenotypic differentiation among populations. Both neutral processes and selection can influence phenotypic differentiation. Altitudinal phenotypic variation is of particular interest in disentangling the interplay between neutral processes and selection in the dynamics of local adaptation processes but remains little explored. We conducted a common garden experiment to study the phenotypic divergence in larval life-history traits among nine populations of the common toad (Bufo bufo) along an altitudinal gradient in France. We further used correlation among population pairwise estimates of quantitative trait (Q_ST) and neutral genetic divergence (F_ST from neutral microsatellite markers), as well as altitudinal difference, to estimate the relative role of divergent selection and neutral genetic processes in phenotypic divergence. We provided evidence for a neutral genetic differentiation resulting from both isolation by distance and difference in altitude. We found evidence for phenotypic divergence along the altitudinal gradient (faster development, lower growth rate and smaller metamorphic size). The correlation between pairwise Q_STs–F_STs and altitude differences suggested that this phenotypic differentiation was most likely driven by altitude-mediated selection rather than by neutral genetic processes. Moreover, we found different divergence patterns for larval traits, suggesting that different selective agents may act on these traits and/or selection on one trait may constrain the evolution on another through genetic correlation. Our study highlighted the need to design more integrative studies on the common toad to unravel the underlying processes of phenotypic divergence and its selective agents in the context of environmental clines.     

Comparison of quantitative and molecular variation in agroforestry populations of the shea tree (Vitellaria paradoxa C.F. Gaertn) in Mali

In this study we investigated the within- and between-population genetic variation using microsatellite markers and quantitative traits of the shea tree, Vitellaria paradoxa, an important agroforestry tree species of the Sudano–Sahelian region in Africa. Eleven populations were sampled across Mali and in northern Côte d’Ivoire. Leaf size and form and growth traits were measured in a progeny test at the nursery stage. Eight microsatellites were used to assess neutral genetic variation. Low levels of heterozygosity were recorded (1.6–3.0 alleles/locus; H_E = 0.25–0.42) and the fixation index (F_IS = −0.227–0.186) was not significantly different from zero suggesting that Hardy–Weinberg equilibrium is encountered in all populations sampled. Quantitative traits exhibited a strong genetic variation between populations and between families within populations. The degree of population differentiation of the quantitative traits (Q_ST = 0.055–0.283, Q_STmean = 0.189) strongly exceeds that in eight microsatellite loci (F_ST = −0.011–0.142, F_STmean = 0.047). Global and pairwise F_ST values were very low and not significantly different from zero suggesting agroforestry practices are amplifying gene flow (Nm = 5.07). The population means for quantitative traits and the rainfall variable were not correlated, showing variation was not linked with this climatic cline. It is suggested that this marked differentiation for quantitative traits, independent of environmental clines and despite a high gene flow, is a result of local adaptation and human selection of shea trees. This process has induced high linkage disequilibrium between underlying loci of polygenic characters.     

Environmental versus Anthropogenic Effects on Population Adaptive Divergence in the Freshwater Snail Lymnaea stagnalis

Repeated pesticide contaminations of lentic freshwater systems located within agricultural landscapes may affect population evolution in non-target organisms, especially in species with a fully aquatic life cycle and low dispersal ability. The issue of evolutionary impact of pollutants is therefore conceptually important for ecotoxicologists. The impact of historical exposure to pesticides on genetic divergence was investigated in the freshwater gastropod Lymnaea stagnalis, using a set of 14 populations from contrasted environments in terms of pesticide and other anthropogenic pressures. The hypothesis of population adaptive divergence was tested on 11 life-history traits, using Q _ST -F _ST comparisons. Despite strong neutral differentiation (mean F _ST = 0.291), five adult traits or parameters were found to be under divergent selection. Conversely, two early expressed traits showed a pattern consistent with uniform selection or trait canalization, and four adult traits appeared to evolve neutrally. Divergent selection patterns were mostly consistent with a habitat effect, opposing pond to ditch and channel populations. Comparatively, pesticide and other human pressures had little correspondence with evolutionary patterns, despite hatching rate impairment associated with global anthropogenic pressure. Globally, analyses revealed high genetic variation both at neutral markers and fitness-related traits in a species used as model in ecotoxicology, providing empirical support for the need to account for genetic and evolutionary components of population response in ecological risk assessment.     

Intraspecific divergence in sperm morphology of the green sea urchin, <Emphasis Type="Italic">Strongylocentrotus droebachiensis</Emphasis>: implications for selection in broadcast spawners

Background  

Sperm morphology can be highly variable among species, but less is known about patterns of population differentiation within species. Most studies of sperm morphometric variation are done in species with internal fertilization, where sexual selection can be mediated by complex mating behavior and the environment of the female reproductive tract. Far less is known about patterns of sperm evolution in broadcast spawners, where reproductive dynamics are largely carried out at the gametic level. We investigated variation in sperm morphology of a broadcast spawner, the green sea urchin (Strongylocentrotus droebachiensis), within and among spawnings of an individual, among individuals within a population, and among populations. We also examined population-level variation between two reproductive seasons for one population. We then compared among-population quantitative genetic divergence (Q _ST) for sperm characters to divergence at neutral microsatellite markers (F _ST).     

Post‐glacial colonization routes coincide with a life‐history breakpoint along a latitudinal gradient

Although adaptive divergence along environmental gradients has repeatedly been demonstrated, the role of post‐glacial colonization routes in determining phenotypic variation along gradients has received little attention. Here, we used a hierarchical Q_ST–F_ST approach to separate the roles of adaptive and neutral processes in shaping phenotypic variation in moor frog (Rana arvalis) larval life histories along a 1,700 km latitudinal gradient across northern Europe. This species has colonized Scandinavia via two routes with a contact zone in northern Sweden. By using neutral SNP and common garden phenotypic data from 13 populations at two temperatures, we showed that most of the variation along the gradient occurred between the two colonizing lineages. We found little phenotypic divergence within the lineages; however, all phenotypic traits were strongly diverged between the southern and northern colonization routes, with higher growth and development rates and larger body size in the north. The Q_ST estimates between the colonization routes were four times higher than F_ST, indicating a prominent role for natural selection. Q_ST within the colonization routes did not generally differ from F_ST, but we found temperature‐dependent adaptive divergence close to the contact zone. These results indicate that lineage‐specific variation can account for much of the adaptive divergence along a latitudinal gradient.     

Natural selection and neutral evolutionary processes contribute to genetic divergence in leaf traits across a precipitation gradient in the tropical oak Quercus oleoides

The impacts of drought are expanding worldwide as a consequence of climate change. However, there is still little knowledge of how species respond to long‐term selection in seasonally dry ecosystems. In this study, we used Q_ST‐F_ST comparisons to investigate (i) the role of natural selection on population genetic differentiation for a set of functional traits related to drought resistance in the seasonally dry tropical oak Quercus oleoides and (ii) the influence of water availability at the site of population origin and in experimental treatments on patterns of trait divergence. We conducted a thorough phenotypic characterization of 1912 seedlings from ten populations growing in field and greenhouse common gardens under replicated watering treatments. We also genotyped 218 individuals from the same set of populations using eleven nuclear microsatellites. Q_ST distributions for leaf lamina area, specific leaf area, leaf thickness and stomatal pore index were higher than F_ST distribution. Results were consistent across growth environments. Genetic differentiation among populations for these functional traits was associated with the index of moisture at the origin of the populations. Together, our results suggest that drought is an important selective agent for Q. oleoides and that differences in length and severity of the dry season have driven the evolution of genetic differences in functional traits.     

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.     

Natural selection contributes to geographic patterns of thermal plasticity in Plantago lanceolata

A long‐standing debate in evolutionary biology concerns the relative importance of different evolutionary forces in explaining phenotypic diversification at large geographic scales. For example, natural selection is typically assumed to underlie divergence along environmental gradients. However, neutral evolutionary processes can produce similar patterns. We collected molecular genetic data from 14 European populations of Plantago lanceolata to test the contributions of natural selection versus neutral evolution to population divergence in temperature‐sensitive phenotypic plasticity of floral reflectance. In P. lanceolata, reflectance plasticity is positively correlated with latitude/altitude. We used population pairwise comparisons between neutral genetic differentiation (F_ST and Jost's D) and phenotypic differentiation (P_ST) to assess the contributions of geographic distance and environmental parameters of the reproductive season in driving population divergence. Data are consistent with selection having shaped large‐scale geographic patterns in thermal plasticity. The aggregate pattern of P_ST versus F_ST was consistent with divergent selection. F_ST explained thermal plasticity differences only when geographic distance was not included in the model. Differences in the extent of cool reproductive season temperatures, and not overall temperature variation, explained plasticity differences independent of distance. Results are consistent with the hypothesis that thermal plasticity is adaptive where growing seasons are shorter and cooler, that is, at high latitude/altitude.     

Geographic clines in wing morphology relate to colonization history in New World but not Old World populations of yellow dung flies

Geographic clines offer insights about putative targets and agents of natural selection as well as tempo and mode of adaptation. However, demographic processes can lead to clines that are indistinguishable from adaptive divergence. Using the widespread yellow dung fly Scathophaga stercoraria (Diptera: Scathophagidae), we examine quantitative genetic differentiation (Q_ST) of wing shape across North America, Europe, and Japan, and compare this differentiation with that of ten microsatellites (F_ST). Morphometric analyses of 28 populations reared at three temperatures revealed significant thermal plasticity, sexual dimorphism, and geographic differentiation in wing shape. In North America morphological differentiation followed the decline in microsatellite variability along the presumed route of recent colonization from the southeast to the northwest. Across Europe, where S. stercoraria presumably existed for much longer time and where no molecular pattern of isolation by distance was evident, clinal variation was less pronounced despite significant morphological differentiation (Q_ST>F_ST). Shape vector comparisons further indicate that thermal plasticity (hot‐to‐cold) does not mirror patterns of latitudinal divergence (south‐to‐north), as might have been expected under a scenario with temperature as the major agent of selection. Our findings illustrate the importance of detailed phylogeographic information when interpreting geographic clines of dispersal traits in an adaptive evolutionary framework.