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.     

Strong phenotypic divergence in spite of low genetic structure in the endemic Mangrove Warbler subspecies (Setophaga petechia xanthotera) of Costa Rica

Despite the enormous advances in genetics, links between phenotypes and genotypes have been made for only a few nonmodel organisms. However, such links can be essential to understand mechanisms of ecological speciation. The Costa Rican endemic Mangrove Warbler subspecies provides an excellent subject to study differentiation with gene flow, as it is distributed along a strong precipitation gradient on the Pacific coast with no strong geographic barriers to isolate populations. Mangrove Warbler populations could be subject to divergent selection driven by precipitation, which influences soil salinity levels, which in turn influences forest structure and food resources. We used single nucleotide polymorphisms (SNPs) and morphological traits to examine the balance between neutral genetic and phenotypic divergence to determine whether selection has acted on traits and genes with functions related to specific environmental variables. We present evidence showing: (a) associations between environmental variables and SNPs, identifying candidate genes related to bill morphology (BMP) and osmoregulation, (b) absence of population genetic structure in neutrally evolving markers, (c) divergence in bill size across the precipitation gradient, and (d) strong phenotypic differentiation (P_ST) which largely exceeds neutral genetic differentiation (F_ST) in bill size. Our results indicate an important role for salinity, forest structure, and resource availability in maintaining phenotypic divergence of Mangrove Warblers through natural selection. Our findings add to the growing body of literature identifying the processes involved in phenotypic differentiation along environmental gradients in the face of gene flow.     

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.     

Local adaptation with high gene flow: temperature parameters drive adaptation to altitude in the common frog (Rana temporaria)

Both environmental and genetic influences can result in phenotypic variation. Quantifying the relative contributions of local adaptation and phenotypic plasticity to phenotypes is key to understanding the effect of environmental variation on populations. Identifying the selective pressures that drive divergence is an important, but often lacking, next step. High gene flow between high‐ and low‐altitude common frog (Rana temporaria) breeding sites has previously been demonstrated in Scotland. The aim of this study was to assess whether local adaptation occurs in the face of high gene flow and to identify potential environmental selection pressures that drive adaptation. Phenotypic variation in larval traits was quantified in R. temporaria from paired high‐ and low‐altitude sites using three common temperature treatments. Local adaptation was assessed using Q_ST–F_ST analyses, and quantitative phenotypic divergence was related to environmental parameters using Mantel tests. Although evidence of local adaptation was found for all traits measured, only variation in larval period and growth rate was consistent with adaptation to altitude. Moreover, this was only evident in the three mountains with the highest high‐altitude sites. This variation was correlated with mean summer and winter temperatures, suggesting that temperature parameters are potentially strong selective pressures maintaining local adaptation, despite high gene flow.     

Divergent selection and phenotypic plasticity during incipient speciation in Lake Victoria cichlid fish

Divergent selection acting on several different traits that cause multidimensional shifts are supposed to promote speciation, but the outcome of this process is highly dependent on the balance between the strength of selection vs. gene flow. Here, we studied a pair of sister species of Lake Victoria cichlids at a location where they hybridize and tested the hypothesis that divergent selection acting on several traits can maintain phenotypic differentiation despite gene flow. To explore the possible role of selection we tested for correlations between phenotypes and environment and compared phenotypic divergence (P_ST) with that based on neutral markers (F_ST). We found indications for disruptive selection acting on male breeding colour and divergent selection acting on several morphological traits. By performing common garden experiments we also separated the environmental and heritable components of divergence and found evidence for phenotypic plasticity in some morphological traits contributing to species differences.     

Remarkable life history polymorphism may be evolving under divergent selection in the silverleaf sunflower

Substantial intraspecific variation in life history is rare and potentially a signal of incipient ecological speciation, if variation is driven by geographically heterogenous natural selection. We present the first report of extensive life history polymorphism in Helianthus argophyllus, the silverleaf sunflower, and examine evidence for its evolution by divergent selection. In 18 populations sampled from across the species range and grown in a common garden, most quantitative traits covaried such that individuals could be assigned to two distinct life history syndromes: tall and late flowering with small initial flowerheads, or short and early flowering with larger initial flowerheads. Helianthus argophyllus exhibits regional genetic structure, but this population structure does not closely correspond with patterns of phenotypic variation. The early‐flowering syndrome is primarily observed in populations from coastal barrier islands, while populations from the nearby mainland coast, although geographically and genetically close, are primarily late flowering. Additionally, several traits are more differentiated among regions than expected based on neutral genetic divergence (Q_ST > F_ST), including the first principal component score corresponding with life history syndrome. This discordance between patterns of phenotypic and genetic variation suggests that divergent selection is driving genetic differences in life history across the species range. If so, the silverleaf sunflower may be in early stages of ecological speciation.     

Female mate choice of male signals is unlikely to promote ecological adaptation in Enchenopa treehoppers (Hemiptera: Membracidae)

A key question in speciation research is how ecological and sexual divergence arise and interact. We tested the hypothesis that mate choice causes local adaptation and ecological divergence using the rationale that the performance~signal trait relationship should parallel the attractiveness~signal trait relationship. We used female fecundity as a measure of ecological performance. We used a species in the Enchenopa binotata treehopper complex, wherein speciation involves adaptation to novel environments and divergence in sexual communication. We used a full‐sibling, split‐family rearing design to estimate genetic correlations (r_G) between fecundity and signal traits, and compared those relationships against population‐level mate preferences for the signal traits. Animal model estimates for r_G between female fecundity and male signal traits overlapped zero—rejecting the hypothesis—but could reflect sample size limitations. The magnitude of r_G correlated with the strength of the mate preferences for the corresponding signal traits, especially for signal frequency, which has the strongest mate preference and the most divergence in the complex. However, signal frequencies favored by the population‐level mate preference are not associated with high fecundity. Therefore, mate preferences do not appear to have been selected to favor high‐performance genotypes. Our findings suggest that ecological and sexual divergence may arise separately, but reinforce each other, during speciation.     

The genomic bases of morphological divergence and reproductive isolation driven by ecological speciation in Senecio (Asteraceae)

Ecological speciation, driven by adaptation to contrasting environments, provides an attractive opportunity to study the formation of distinct species, and the role of selection and genomic divergence in this process. Here, we focus on a particularly clear‐cut case of ecological speciation to reveal the genomic bases of reproductive isolation and morphological differences between closely related Senecio species, whose recent divergence within the last ~200 000 years was likely driven by the uplift of Mt. Etna (Sicily). These species form a hybrid zone, yet remain morphologically and ecologically distinct, despite active gene exchange. Here, we report a high‐density genetic map of the Senecio genome and map hybrid breakdown to one large and several small quantitative trait loci (QTL). Loci under diversifying selection cluster in three 5 cM regions which are characterized by a significant increase in relative (F_ST), but not absolute (d_XY), interspecific differentiation. They also correspond to some of the regions of greatest marker density, possibly corresponding to ‘cold‐spots’ of recombination, such as centromeres or chromosomal inversions. Morphological QTL for leaf and floral traits overlap these clusters. We also detected three genomic regions with significant transmission ratio distortion (TRD), possibly indicating accumulation of intrinsic genetic incompatibilities between these recently diverged species. One of the TRD regions overlapped with a cluster of high species differentiation, and another overlaps the large QTL for hybrid breakdown, indicating that divergence of these species may have occurred due to a complex interplay of ecological divergence and accumulation of intrinsic genetic incompatibilities.     

Isolation by distance,not incipient ecological speciation,explains genetic differentiation in an Andean songbird (Aves: Furnariidae: Cranioleuca antisiensis,Line‐cheeked Spinetail) despite near threefold body size change across an environmental gradient

During the process of ecological speciation, reproductive isolation results from divergent natural selection and leads to a positive correlation between genetic divergence and adaptive phenotypic divergence, that is, isolation by adaptation (IBA). In natural populations, phenotypic differentiation is often autocorrelated with geographic distance, making IBA difficult to distinguish from the neutral expectation of isolation by distance (IBD). We examined these two alternatives in a dramatic case of clinal phenotypic variation in an Andean songbird, the Line‐cheeked Spinetail (Cranioleuca antisiensis). At its geographic extremes, this species shows a near threefold difference in body mass (11.5 to 31.0 g) with marked plumage differences. We analysed phenotypic, environmental and genetic data (5,154 SNPs) from 172 individuals and 19 populations sampled along its linear distribution in the Andes. We found that body mass was tightly correlated with environmental temperature, consistent with local adaptation as per Bergmann's rule. Using a P_ST–F_ST analysis, we found additional support for natural selection driving body mass differentiation, but these results could also be explained by environment‐mediated phenotypic plasticity. When we assessed the relative support for patterns of IBA and IBD using variance partitioning, we found that IBD was the best explanation for genetic differentiation along the cline. Adaptive phenotypic or environmental divergence can reduce gene flow, a pattern interpreted as evidence of ecological speciation's role in diversification. Our results provide a counterexample to this interpretation. Despite conditions conducive to ecological speciation, our results suggest that dramatic size and environmental differentiation within C. antisiensis are not limiting gene flow.     

Population genetic structure and a possible role for selection in driving phenotypic divergence in a rainbowfish (Melanotaeniidae)

Quantitative genetic divergence may be driven by drift or selection. The rainbowfish Melanotaenia australis exhibits phenotypic divergence among populations in Western Australia, although the mechanisms driving this divergence are unknown. We used microsatellites to assess neutral genetic divergence (F_ST), and found a hierarchical pattern of subdivision consistent with low divergence between upstream and downstream populations (within drainages), moderate divergence between drainages (within regions), and high divergence between regions. Using a common‐garden approach, we measured quantitative genetic divergence in phenotypic traits (Q_ST). By comparing this to expectations from neutral processes (F_ST), we concluded that the effect of selection varies depending on the spatial scale considered. Within drainages, selection may be causing divergence between upstream and downstream phenotypes but, between regions, selection appears to homogenize phenotypes. This highlights the importance of spatial scale in studies of this kind, and suggests that, because variance in selection pressures can drive speciation, polymorphism in M. australis may represent speciation in action. © 2010 The Linnean Society of London, Biological Journal of the Linnean Society, 2011, 102 , 144–160.     

Diversification across a heterogeneous landscape

Adaptation to contrasting environments across a heterogeneous landscape favors the formation of ecotypes by promoting ecological divergence. Patterns of fitness variation in the field can show whether natural selection drives local adaptation and ecotype formation. However, to demonstrate a link between ecological divergence and speciation, local adaptation must have consequences for reproductive isolation. Using contrasting ecotypes of an Australian wildflower, Senecio lautus in common garden experiments, hybridization experiments, and reciprocal transplants, we assessed how the environment shapes patterns of adaptation and the consequences of adaptive divergence for reproductive isolation. Local adaptation was strong between ecotypes, but weaker between populations of the same ecotype. F1 hybrids exhibited heterosis, but crosses involving one native parent performed better than those with two foreign parents. In a common garden experiment, F2 hybrids exhibited reduced fitness compared to parentals and F1 hybrids, suggesting that few genetic incompatibilities have accumulated between populations adapted to contrasting environments. Our results show how ecological differences across the landscape have created complex patterns of local adaptation and reproductive isolation, suggesting that divergent natural selection has played a fundamental role in the early stages of species diversification.     

Global patterns of population genetic differentiation in seed plants

Evaluating the factors that drive patterns of population differentiation in plants is critical for understanding several biological processes such as local adaptation and incipient speciation. Previous studies have given conflicting results regarding the significance of pollination mode, seed dispersal mode, mating system, growth form and latitudinal region in shaping patterns of genetic structure, as estimated by F_ST values, and no study to date has tested their relative importance together across a broad scale. Here, we assembled a 337‐species data set for seed plants from publications with data on F_ST from nuclear markers and species traits, including variables pertaining to the sampling scheme of each study. We used species traits, while accounting for sampling variables, to perform phylogenetic multiple regressions. Results demonstrated that F_ST values were higher for tropical, mixed‐mating, non‐woody species pollinated by small insects, indicating greater population differentiation, and lower for temperate, outcrossing trees pollinated by wind. Among the factors we tested, latitudinal region explained the largest portion of variance, followed by pollination mode, mating system and growth form, while seed dispersal mode did not significantly relate to F_ST. Our analyses provide the most robust and comprehensive evaluation to date of the main ecological factors predicted to drive population differentiation in seed plants, with important implications for understanding the basis of their genetic divergence. Our study supports previous findings showing greater population differentiation in tropical regions and is the first that we are aware of to robustly demonstrate greater population differentiation in species pollinated by small insects.     

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.     

MULTIVARIATE ADAPTATION BUT NO INCREASE IN COMPETITIVE ABILITY IN INVASIVE GERANIUM CAROLINIANUM L. (GERANIACEAE)

Adaptive evolution can affect the successful establishment of invasive species, but changes in selective pressures, loss of genetic variation in relevant traits, and/or altered trait correlations can make adaptation difficult to predict. We used a common‐garden experiment to assess trait correlations and patterns of adaptation in the invasive plant, Geranium carolinianum, sampled across 20 populations in its native (United States) and invasive (China) ranges. We used multivariate Q_ST – F_ST tests to determine if phenotypic differences between countries are attributable to adaptation. We also compared population‐level variation within each country to assess whether local adaptation resulted in similar multivariate phenotypes in the United States and China. Between countries, most phenotypic differences are indistinguishable from genetic drift, although we detected a signature of adaptation to the colder, drier winters in China. There was no evidence for increases in invasive traits in China. Within countries, strong multivariate adaptation appears to be driven by latitudinal climatic variation in the United States, but not in China. Additionally, adaptive trait combinations as well as their underlying correlations differ between the two countries, indicating that adaptation in invasive populations does not parallel patterns in native populations due to differences in selection pressures, genetic constraints, or both.     

Genomewide patterns of variation in genetic diversity are shared among populations,species and higher‐order taxa

Genomewide screens of genetic variation within and between populations can reveal signatures of selection implicated in adaptation and speciation. Genomic regions with low genetic diversity and elevated differentiation reflective of locally reduced effective population sizes (N_e) are candidates for barrier loci contributing to population divergence. Yet, such candidate genomic regions need not arise as a result of selection promoting adaptation or advancing reproductive isolation. Linked selection unrelated to lineage‐specific adaptation or population divergence can generate comparable signatures. It is challenging to distinguish between these processes, particularly when diverging populations share ancestral genetic variation. In this study, we took a comparative approach using population assemblages from distant clades assessing genomic parallelism of variation in N_e. Utilizing population‐level polymorphism data from 444 resequenced genomes of three avian clades spanning 50 million years of evolution, we tested whether population genetic summary statistics reflecting genomewide variation in N_e would covary among populations within clades, and importantly, also among clades where lineage sorting has been completed. All statistics including population‐scaled recombination rate (ρ), nucleotide diversity (π) and measures of genetic differentiation between populations (F_ST, PBS, d_xy) were significantly correlated across all phylogenetic distances. Moreover, genomic regions with elevated levels of genetic differentiation were associated with inferred pericentromeric and subtelomeric regions. The phylogenetic stability of diversity landscapes and stable association with genomic features support a role of linked selection not necessarily associated with adaptation and speciation in shaping patterns of genomewide heterogeneity in genetic diversity.     

Rapid adaptive divergence between ecotypes of an aquatic isopod inferred from FST–QST analysis

Divergent natural selection is often thought to be the principal factor driving phenotypic differentiation between populations. We studied two ecotypes of the aquatic isopod Asellus aquaticus which have diverged in parallel in several Swedish lakes. In these lakes, isopods from reed belts along the shores colonized new stonewort stands in the centre of the lakes and rapid phenotypic changes in size and pigmentation followed after colonization. We investigated if selection was likely to be responsible for these observed phenotypic changes using indirect inferences of selection (F_ST–Q_ST analysis). Average Q_ST for seven quantitative traits were higher than the average F_ST between ecotypes for putatively neutral markers (AFLPs). This suggests that divergent natural selection has played an important role during this rapid diversification. In contrast, the average Q_ST between the different reed ecotype populations was not significantly different from the mean F_ST. Genetic drift could therefore not be excluded as an explanation for the minor differences between allopatric populations inhabiting the same source habitat. We complemented this traditional F_ST–Q_ST approach by comparing the F_ST distributions across all loci (n = 67–71) with the Q_ST for each of the seven traits. This analysis revealed that pigmentation traits had diverged to a greater extent and at higher evolutionary rates than size‐related morphological traits. In conclusion, this extended and detailed type of F_ST–Q_ST analysis provides a powerful method to infer adaptive phenotypic divergence between populations. However, indirect inferences about the operation of divergent selection should be analyzed on a per‐trait basis and complemented with detailed ecological information.     

ADAPTATION TO DESICCATION FAILS TO GENERATE PRE‐ AND POSTMATING ISOLATION IN REPLICATE DROSOPHILA MELANOGASTER LABORATORY POPULATIONS

Many laboratory speciation experiments have raised allopatric populations in different environments to determine whether reproductive isolation evolves as a byproduct of adaptation (a form of ecological speciation). Few, however, have addressed the evolution of both pre‐ and postmating isolation or investigated the conditions affecting the process. We present results of an evolution experiment in which 12 lines of Drosophila melanogaster were derived from a common population and then independently evolved for more than 57 generations under alternative selection treatments. Six “desiccation” populations had access to food and water removed during a period of their adult lives generating strong desiccation selection, and six “starvation” populations had access to food but not water removed for the same period, generating a mild starvation stress. Substantial divergence of cuticular hydrocarbons occurred between the desiccation and starvation populations, key traits that have been implicated in sexual isolation in Drosophila. Despite this divergence, there was no detectable premating isolation between desiccation and starvation populations and postmating isolation was asymmetrical: the fitness of F1 hybrids was reduced in the desiccation but not the starvation environment. This asymmetry was likely caused by the absence of divergent selection: adaptation to desiccation appears to have come at no cost to performance in the starvation environment. Novel environments are thought to promote the evolution of reproductive isolation. Understanding the conditions that favor or hamper this remains a key challenge for speciation research.     

Understanding the mechanisms of antitropical divergence in the seabird White‐faced Storm‐petrel (Procellariiformes: Pelagodroma marina) using a multilocus approach

Analytical methods that apply coalescent theory to multilocus data have improved inferences of demographic parameters that are critical to understanding population divergence and speciation. In particular, at the early stages of speciation, it is important to implement models that accommodate conflicting gene trees, and benefit from the presence of shared polymorphisms. Here, we employ eleven nuclear loci and the mitochondrial control region to investigate the phylogeography and historical demography of the pelagic seabird White‐faced Storm‐petrel (Pelagodroma marina) by sampling subspecies across its antitropical distribution. Groups are all highly differentiated: global mitochondrial Φ_ST = 0.89 (P < 0.01) and global nuclear Φ_ST varies between 0.22 and 0.83 (all P < 0.01). The complete lineage sorting of the mitochondrial locus between hemispheres is corroborated by approximately half of the nuclear genealogies, suggesting a long‐term antitropical divergence in isolation. Coalescent‐based estimates of demographic parameters suggest that hemispheric divergence of P. marina occurred approximately 840 000 ya (95% HPD 582 000–1 170 000), in the absence of gene flow, and divergence within the Southern Hemisphere occurred 190 000 ya (95% HPD 96 000–600 000), both probably associated with the profound palaeo‐oceanographic changes of the Pleistocene. A fledgling sampled in St Helena (tropical South Atlantic) suggests recent colonization from the Northern Hemisphere. Despite the great potential for long‐distance dispersal, P. marina antitropical groups have been evolving as independent, allopatric lineages, and divergence is probably maintained by philopatry coupled with asynchronous reproductive phenology and local adaptation.     

Role of selection and gene flow in population differentiation at the edge vs. interior of the species range differing in climatic conditions

Evaluating the relative importance of neutral and adaptive processes as determinants of population differentiation across environments is a central theme of evolutionary biology. We applied the Q_ST–F_ST comparison flanked by a direct test for local adaptation to infer the role of climate‐driven selection and gene flow in population differentiation of an annual grass Avena sterilis in two distinct parts of the species range, edge and interior, which represent two globally different climates, desert and Mediterranean. In a multiyear reciprocal transplant experiment, the plants of desert and Mediterranean origin demonstrated home advantage, and population differentiation in several phenotypic traits related to reproduction exceeded neutral predictions, as determined by comparisons of Q_ST values with theoretical F_ST distributions. Thus, variation in these traits likely resulted from local adaptation to desert and Mediterranean environments. The two separate common garden experiments conducted with different experimental design revealed that two population comparisons, in contrast to multi‐population comparisons, are likely to detect population differences in virtually every trait, but many of these differences reflect effects of local rather than regional environment. We detected a general reduction in neutral (SSR) genetic variation but not in adaptive quantitative trait variation in peripheral desert as compared with Mediterranean core populations. On the other hand, the molecular data indicated intensive gene flow from the Mediterranean core towards desert periphery. Although species range position in our study (edge vs. interior) was confounded with climate (desert vs. Mediterranean), the results suggest that the gene flow from the species core does not have negative consequences for either performance of the peripheral plants or their adaptive potential.     

Low genome‐wide divergence between two lizard populations with high adaptive phenotypic differentiation

Usually, adaptive phenotypic differentiation is paralleled by genetic divergence between locally adapted populations. However, adaptation can also happen in a scenario of nonsignificant genetic divergence due to intense gene flow and/or recent differentiation. While this phenomenon is rarely published, findings on incipient ecologically driven divergence or isolation by adaptation are relatively common, which could confound our understanding about the frequency at which they actually occur in nature. Here, we explore genome‐wide traces of divergence between two populations of the lacertid lizard Psammodromus algirus separated by a 600 m elevational gradient. These populations seem to be differentially adapted to their environments despite showing low levels of genetic differentiation (according to previously studies of mtDNA and microsatellite data). We performed a search for outliers (i.e., loci subject to selection) trying to identify specific loci with F_ST statistics significantly higher than those expected on the basis of overall, genome‐wide estimates of genetic divergence. We find that local phenotypic adaptation (in terms of a wide diversity of characters) was not accompanied by genome‐wide differentiation, even when we maximized the chances of unveiling such differentiation at particular loci with F_ST‐based outlier detection tests. Instead, our analyses confirmed the lack of genome‐wide differentiation on the basis of more than 70,000 SNPs, which is concordant with a scenario of local adaptation without isolation by environment. Our results add evidence to previous studies in which local adaptation does not lead to any kind of isolation (or early stages of ecological speciation), but maintains phenotypic divergence despite the lack of a differentiated genomic background.