Putative chromosomal rearrangements are associated primarily with ecotype divergence rather than geographic separation in an intertidal,poorly dispersing snail

Littorina saxatilis is becoming a model system for understanding the genomic basis of ecological speciation. The parallel formation of crab‐adapted ecotypes that exhibit partial reproductive isolation from wave‐adapted ecotypes has enabled genomic investigation of conspicuous shell traits. Recent genomic studies suggest that chromosomal rearrangements may enable ecotype divergence by reducing gene flow. However, the genomic architecture of traits that are divergent between ecotypes remains poorly understood. Here, we use 11,504 single nucleotide polymorphism (SNP) markers called using the recently released L. saxatilis genome to genotype 462 crab ecotype, wave ecotype and phenotypically intermediate Spanish L. saxatilis individuals with scored phenotypes. We used redundancy analysis to study the genetic architecture of loci associated with shell shape, shape corrected for size, shell size and shell ornamentation, and to compare levels of co‐association among different traits. We discovered 341 SNPs associated with shell traits. Loci associated with trait divergence between ecotypes were often located inside putative chromosomal rearrangements recently characterized in Swedish L. saxatilis. In contrast, we found that shell shape corrected for size varied primarily by geographic site rather than by ecotype and showed little association with these putative rearrangements. We conclude that genomic regions of elevated divergence inside putative rearrangements were associated with divergence of L. saxatilis ecotypes along steep environmental axes—consistent with models of adaptation with gene flow—but were not associated with divergence among the three geographical sites. Our findings support predictions from models indicating the importance of genomic regions of reduced recombination allowing co‐association of loci during ecological speciation with ongoing gene flow.     

Genomic divergence between Spanish Littorina saxatilis ecotypes unravels limited admixture and extensive parallelism associated with population history

The rough periwinkle, Littorina saxatilis, is a model system for studying parallel ecological speciation in microparapatry. Phenotypically parallel wave‐adapted and crab‐adapted ecotypes that hybridize within the middle shore are replicated along the northwestern coast of Spain and have likely arisen from two separate glacial refugia. We tested whether greater geographic separation corresponding to reduced opportunity for contemporary or historical gene flow between parallel ecotypes resulted in less parallel genomic divergence. We sequenced double‐digested restriction‐associated DNA (ddRAD) libraries from individual snails from upper, mid, and low intertidal levels of three separate sites colonized from two separate refugia. Outlier analysis of 4256 SNP markers identified 34.4% sharing of divergent loci between two geographically close sites; however, these sites each shared only 9.9%–15.1% of their divergent loci with a third more‐distant site. STRUCTURE analysis revealed that genotypes from only three of 166 phenotypically intermediate mid‐shore individuals appeared to result from recent hybridization, suggesting that hybrids cannot be reliably identified using shell traits. Hierarchical AMOVA indicated that the primary source of genomic differentiation was geographic separation, but also revealed greater similarity of the same ecotype across the two geographically close sites than previously estimated with dominant markers. These results from a model system for ecological speciation suggest that genomic parallelism is affected by the opportunity for historical or contemporary gene flow between populations.     

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.     

Genomics of Rapid Incipient Speciation in Sympatric Threespine Stickleback

Ecological speciation is the process by which reproductively isolated populations emerge as a consequence of divergent natural or ecologically-mediated sexual selection. Most genomic studies of ecological speciation have investigated allopatric populations, making it difficult to infer reproductive isolation. The few studies on sympatric ecotypes have focused on advanced stages of the speciation process after thousands of generations of divergence. As a consequence, we still do not know what genomic signatures of the early onset of ecological speciation look like. Here, we examined genomic differentiation among migratory lake and resident stream ecotypes of threespine stickleback reproducing in sympatry in one stream, and in parapatry in another stream. Importantly, these ecotypes started diverging less than 150 years ago. We obtained 34,756 SNPs with restriction-site associated DNA sequencing and identified genomic islands of differentiation using a Hidden Markov Model approach. Consistent with incipient ecological speciation, we found significant genomic differentiation between ecotypes both in sympatry and parapatry. Of 19 islands of differentiation resisting gene flow in sympatry, all were also differentiated in parapatry and were thus likely driven by divergent selection among habitats. These islands clustered in quantitative trait loci controlling divergent traits among the ecotypes, many of them concentrated in one region with low to intermediate recombination. Our findings suggest that adaptive genomic differentiation at many genetic loci can arise and persist in sympatry at the very early stage of ecotype divergence, and that the genomic architecture of adaptation may facilitate this.     

Genetic and morphological divergence between Littorina fabalis ecotypes in Northern Europe

Low dispersal marine intertidal species facing strong divergent selective pressures associated with steep environmental gradients have a great potential to inform us about local adaptation and reproductive isolation. Among these, gastropods of the genus Littorina offer a unique system to study parallel phenotypic divergence resulting from adaptation to different habitats related with wave exposure. In this study, we focused on two Littorina fabalis ecotypes from Northern European shores and compared patterns of habitat‐related phenotypic and genetic divergence across three different geographic levels (local, regional and global). Geometric morphometric analyses revealed that individuals from habitats moderately exposed to waves usually present a larger shell size with a wider aperture than those from sheltered habitats. The phenotypic clustering of L. fabalis by habitat across most locations (mainly in terms of shell size) support an important role of ecology in morphological divergence. A genome scan based on amplified fragment length polymorphisms (AFLPs) revealed a heterogeneous pattern of differentiation across the genome between populations from the two different habitats, suggesting ecotype divergence in the presence of gene flow. The contrasting patterns of genetic structure between nonoutlier and outlier loci, and the decreased sharing of outlier loci with geographic distance among locations are compatible with parallel evolution of phenotypic divergence, with an important contribution of gene flow and/or ancestral variation. In the future, model‐based inference studies based on sequence data across the entire genome will help unravelling these evolutionary hypotheses, improving our knowledge about adaptation and its influence on diversification within the marine realm.     

The role of local ecology during hybridization at the initial stages of ecological speciation in a marine snail

Hybrid zones of ecologically divergent populations are ideal systems to study the interaction between natural selection and gene flow during the initial stages of speciation. Here, we perform an amplified fragment length polymorphism (AFLP) genome scan in parallel hybrid zones between divergent ecotypes of the marine snail Littorina saxatilis, which is considered a model case for the study of ecological speciation. Ridged‐Banded (RB) and Smooth‐Unbanded (SU) ecotypes are adapted to different shore levels and microhabitats, although they present a sympatric distribution at the mid‐shore where they meet and mate (partially assortatively). We used shell morphology, outlier and nonoutlier AFLP loci from RB, SU and hybrid specimens captured in sympatry to determine the level of phenotypic and genetic introgression. We found different levels of introgression at parallel hybrid zones and nonoutlier loci showed more gene flow with greater phenotypic introgression. These results were independent from the phylogeography of the studied populations, but not from the local ecological conditions. Genetic variation at outlier loci was highly correlated with phenotypic variation. In addition, we used the relationship between genetic and phenotypic variation to estimate the heritability of morphological traits and to identify potential Quantitative Trait Loci to be confirmed in future crosses. These results suggest that ecology (exogenous selection) plays an important role in this hybrid zone. Thus, ecologically based divergent natural selection is responsible, simultaneously, for both ecotype divergence and hybridization. On the other hand, genetic introgression occurs only at neutral loci (nonoutliers). In the future, genome‐wide studies and controlled crosses would give more valuable information about this process of speciation in the face of gene flow.     

Genome‐wide association analyses reveal polygenic genomic architecture underlying divergent shell morphology in Spanish Littorina saxatilis ecotypes

Gene flow between diverging populations experiencing dissimilar ecological conditions can theoretically constrain adaptive evolution. To minimize the effect of gene flow, alleles underlying traits essential for local adaptation are predicted to be located in linked genome regions with reduced recombination. Local reduction in gene flow caused by selection is expected to produce elevated divergence in these regions. The highly divergent crab‐adapted and wave‐adapted ecotypes of the marine snail Littorina saxatilis present a model system to test these predictions. We used genome‐wide association (GWA) analysis of geometric morphometric shell traits associated with microgeographic divergence between the two L. saxatilis ecotypes within three separate sampling sites. A total of 477 snails that had individual geometric morphometric data and individual genotypes at 4,066 single nucleotide polymorphisms (SNPs) were analyzed using GWA methods that corrected for population structure among the three sites. This approach allowed dissection of the genomic architecture of shell shape divergence between ecotypes across a wide geographic range, spanning two glacial lineages. GWA revealed 216 quantitative trait loci (QTL) with shell size or shape differences between ecotypes, with most loci explaining a small proportion of phenotypic variation. We found that QTL were evenly distributed across 17 linkage groups, and exhibited elevated interchromosomal linkage, suggesting a genome‐wide response to divergent selection on shell shape between the two ecotypes. Shell shape trait‐associated loci showed partial overlap with previously identified outlier loci under divergent selection between the two ecotypes, supporting the hypothesis of diversifying selection on these genomic regions. These results suggest that divergence in shell shape between the crab‐adapted and wave‐adapted ecotypes is produced predominantly by a polygenic genomic architecture with positive linkage disequilibrium among loci of small effect.     

Repeated and predictable patterns of ecotypic differentiation during a biological invasion: lake–stream divergence in parapatric Swiss stickleback

The relative importance of ecological selection and geographical isolation in promoting and constraining genetic and phenotypic differentiation among populations is not always obvious. Interacting with divergent selection, restricted opportunity for gene flow may in some cases be as much a cause as a consequence of adaptation, with the latter being a hallmark of ecological speciation. Ecological speciation is well studied in parts of the native range of the three‐spined stickleback. Here, we study this process in a recently invaded part of its range. Switzerland was colonized within the past 140 years from at least three different colonization events involving different stickleback lineages. They now occupy diverse habitats, ranging from small streams to the pelagic zone of large lakes. We use replicated systems of parapatric lake and stream populations, some of which trace their origins to different invasive lineages, to ask (i) whether phenotypic divergence occurred among populations inhabiting distinct habitats, (ii) whether trajectories of phenotypic divergence follow predictable parallel patterns and (iii) whether gene flow constrains divergent adaptation or vice versa. We find consistent phenotypic divergence between populations occupying distinct habitats. This involves parallel evolution in several traits with known ecological relevance in independent evolutionary lineages. Adaptive divergence supersedes homogenizing gene flow even at a small spatial scale. We find evidence that adaptive phenotypic divergence places constraints on gene flow over and above that imposed by geographical distance, signalling the early onset of ecological speciation.     

The role of structural genomic variants in population differentiation and ecotype formation in Timema cristinae walking sticks

Theory predicts that structural genomic variants such as inversions can promote adaptive diversification and speciation. Despite increasing empirical evidence that adaptive divergence can be triggered by one or a few large inversions, the degree to which widespread genomic regions under divergent selection are associated with structural variants remains unclear. Here we test for an association between structural variants and genomic regions that underlie parallel host‐plant‐associated ecotype formation in Timema cristinae stick insects. Using mate‐pair resequencing of 20 new whole genomes we find that moderately sized structural variants such as inversions, deletions and duplications are widespread across the genome, being retained as standing variation within and among populations. Using 160 previously published, standard‐orientation whole genome sequences we find little to no evidence that the DNA sequences within inversions exhibit accentuated differentiation between ecotypes. In contrast, a formerly described large region of reduced recombination that harbours genes controlling colour‐pattern exhibits evidence for accentuated differentiation between ecotypes, which is consistent with differences in the frequency of colour‐pattern morphs between host‐associated ecotypes. Our results suggest that some types of structural variants (e.g., large inversions) are more likely to underlie adaptive divergence than others, and that structural variants are not required for subtle yet genome‐wide genetic differentiation with gene flow.     

DIVERGENCE IS FOCUSED ON FEW GENOMIC REGIONS EARLY IN SPECIATION: INCIPIENT SPECIATION OF SUNFLOWER ECOTYPES

Early in speciation, as populations undergo the transition from local adaptation to incipient species, is when a number of transient, but potentially important, processes appear to be most easily detected. These include signatures of selective sweeps that can point to asymmetry in selection between habitats, divergence hitchhiking, and associations of adaptive genes with environments. In a genomic comparison of ecotypes of the prairie sunflower, Helianthus petiolaris, occurring at Great Sand Dunes National Park and Preserve (Colorado), we found that selective sweeps were mainly restricted to the dune ecotype and that there was variation across the genome in whether proximity to the nondune population constrained or promoted divergence. The major regions of divergence were few and large between ecotypes, in contrast with an interspecific comparison between H. petiolaris and a sympatric congener, Helianthus annuus. In general, the large regions of divergence observed in the ecotypic comparison swamped locus‐specific associations with environmental variables. In both comparisons, regions of high divergence occurred in portions of the genetic map with high marker density, probably reflecting regions of low recombination. The difference in genomic distributions of highly divergent regions between ecotypic and interspecific comparisons highlights the value of studies spanning the spectrum of speciation in related taxa.     

Divergent Macroparasite Infections in Parapatric Swiss Lake-Stream Pairs of Threespine Stickleback (Gasterosteus aculeatus)

Spatial heterogeneity in diversity and intensity of parasitism is a typical feature of most host-parasite interactions, but understanding of the evolutionary implications of such variation is limited. One possible outcome of infection heterogeneities is parasite-mediated divergent selection between host populations, ecotypes or species which may facilitate the process of ecological speciation. However, very few studies have described infections in population-pairs along the speciation continuum from low to moderate or high degree of genetic differentiation that would address the possibility of parasite-mediated divergent selection in the early stages of the speciation process. Here we provide an example of divergent parasitism in freshwater fish ecotypes by examining macroparasite infections in threespine stickleback (Gasterosteus aculeatus) of four Swiss lake systems each harbouring parapatric lake-stream ecotype pairs. We demonstrate significant differences in infections within and between the pairs that are driven particularly by the parasite taxa transmitted to fish from benthic invertebrates. The magnitude of the differences tended to correlate positively with the extent of neutral genetic differentiation between the parapatric lake and stream populations of stickleback, whereas no such correlation was found among allopatric populations from similar or contrasting habitats. This suggests that genetic differentiation is unrelated to the magnitude of parasite infection contrasts when gene flow is constrained by geographical barriers while in the absence of physical barriers, genetic differentiation and the magnitude of differences in infections tend to be positively correlated.     

Pollinator shifts between Ophrys sphegodes populations: might adaptation to different pollinators drive population divergence?

Local adaptation to different pollinators is considered one of the possible initial stages of ecological speciation as reproductive isolation is a by‐product of the divergence in pollination systems. However, pollinator‐mediated divergent selection will not necessarily result in complete reproductive isolation, because incipient speciation is often overcome by gene flow. We investigated the potential of pollinator shift in the sexually deceptive orchids Ophrys sphegodes and Ophrys exaltata and compared the levels of floral isolation vs. genetic distance among populations with contrasting predominant pollinators. We analysed floral hydrocarbons as a proxy for floral divergence between populations. Floral adoption of pollinators and their fidelity was tested using pollinator choice experiments. Interpopulation gene flow and population differentiation levels were estimated using AFLP markers. The Tyrrhenian O. sphegodes population preferentially attracted the pollinator bee Andrena bimaculata, whereas the Adriatic O. sphegodes population exclusively attracted A. nigroaenea. Significant differences in scent component proportions were identified in O. sphegodes populations that attracted different preferred pollinators. High interpopulation gene flow was detected, but populations were genetically structured at species level. The high interpopulation gene flow levels independent of preferred pollinators suggest that local adaptation to different pollinators has not (yet) generated detectable genome‐wide separation. Alternatively, despite extensive gene flow, few genes underlying floral isolation remain differentiated as a consequence of divergent selection. Different pollination ecotypes in O. sphegodes might represent a local selective response imposed by temporal variation in a geographical mosaic of pollinators as a consequence of the frequent disturbance regimes typical of Ophrys habitats.     

Immigrant and extrinsic hybrid seed inviability contribute to reproductive isolation between forest and dune ecotypes of Epipactis helleborine (Orchidaceae)

Reproductive isolation caused by divergent natural selection arising from differences between ecological environments or ecological interactions represents a major mechanism contributing to speciation, but its relative importance is poorly known. In this study, controlled reciprocal crossings and seed germination experiments were combined with genetic and morphometric analyses to test the hypothesis that previously described differences in mycorrhizal communities between forest and dune ecotypes of Epipactis helleborine were sufficiently strong to create complete reproductive isolation between the two ecotypes. Molecular analyses using 770 SNP markers showed that the two ecotypes were genetically distinct and that populations of the dune ecotype were genetically impoverished compared to populations of the forest ecotype. Morphologically, the two ecotypes were also significantly different, with plants of the dune ecotype generally being smaller than plants of the forest ecotype. The results further showed that immigrant seeds had a significantly lower probability of protocorm formation than native seeds, indicating strong immigrant inviability. Although both ecotypes were able to cross easily and to produce a large number of viable seeds, hybrid seeds showed significantly lower protocorm formation than pure seeds, further contributing to reproductive isolation. Overall, these results indicate that interfertile populations of a widespread orchid adapting to contrasting environments diverge as a consequence of concurrent selection acting against immigrants and hybrids and suggest that mycorrhizal fungi can play a role in the early stages of plant speciation.     

Size selection by a gape‐limited predator of a marine snail: Insights into magic traits for speciation

The intertidal snail Littorina saxatilis has repeatedly evolved two parallel ecotypes assumed to be wave adapted and predatory shore crab adapted, but the magnitude and targets of predator‐driven selection are unknown. In Spain, a small, wave ecotype with a large aperture from the lower shore and a large, thick‐shelled crab ecotype from the upper shore meet in the mid‐shore and show partial size‐assortative mating. We performed complementary field tethering and laboratory predation experiments; the first set compared the survival of two different size‐classes of the crab ecotype while the second compared the same size‐class of the two ecotypes. In the first set, the large size‐class of the crab ecotype survived significantly better than the small size‐class both on the upper shore and in the laboratory. In the second set, the small size‐class of the crab ecotype survived substantially better than that of the wave ecotype both on the upper shore and in the laboratory. Shell‐breaking predation on tethered snails was almost absent within the lower shore. In the laboratory shore crabs (Pachygrapsus marmoratus) with larger claw heights selected most strongly against the small size‐class of the crab ecotype, whereas those with medium claw heights selected most strongly against the thin‐shelled wave ecotype. Sexual maturity occurred at a much larger size in the crab ecotype than in the wave ecotype. Our results showed that selection on the upper shore for rapid attainment of a size refuge from this gape‐limited predator favors large size, thick shells, and late maturity. Model parameterization showed that size‐selective predation restricted to the upper shore resulted in the evolution of the crab ecotype despite gene flow from the wave ecotype snails living on the lower shore. These results on gape‐limited predation and previous ones showing size‐assortative mating between ecotypes suggest that size may represent a magic trait for the thick‐shelled ecotype.     

Postglacial ecotype formation under outcrossing and self‐fertilization in Arabidopsis lyrata

The formation of ecotypes has been invoked as an important driver of postglacial biodiversity, because many species colonized heterogeneous habitats and experienced divergent selection. Ecotype formation has been predominantly studied in outcrossing taxa, while far less attention has been paid to the implications of mating system shifts. Here, we addressed whether substrate‐related ecotypes exist in selfing and outcrossing populations of Arabidopsis lyrata subsp. lyrata and whether the genomic footprint differs between mating systems. The North American subspecies colonized both rocky and sandy habitats during postglacial range expansion and shifted the mating system from predominantly outcrossing to predominantly selfing in a number of regions. We performed an association study on pooled whole‐genome sequence data of 20 selfing or outcrossing populations, which suggested genes involved in adaptation to substrate. Motivated by enriched gene ontology terms, we compared root growth between plants from the two substrates in a common environment and found that plants originating from sand grew roots faster and produced more side roots, independent of mating system. Furthermore, single nucleotide polymorphisms associated with substrate‐related ecotypes were more clustered among selfing populations. Our study provides evidence for substrate‐related ecotypes in A. lyrata and divergence in the genomic footprint between mating systems. The latter is the likely result of selfing populations having experienced divergent selection on larger genomic regions due to higher genome‐wide linkage disequilibrium.     

VARIABLE PROGRESS TOWARD ECOLOGICAL SPECIATION IN PARAPATRY: STICKLEBACK ACROSS EIGHT LAKE-STREAM TRANSITIONS

Divergent selection between contrasting habitats can sometimes drive adaptive divergence and the evolution of reproductive isolation in the face of initially high gene flow. "Progress" along this ecological speciation pathway can range from minimal divergence to full speciation. We examine this variation for threespine stickleback fish that evolved independently across eight lake-stream habitat transitions. By quantifying stickleback diets, we show that lake-stream transitions usually coincide with limnetic-benthic ecotones. By measuring genetically based phenotypes, we show that these ecotones often generate adaptive divergence in foraging morphology. By analyzing neutral genetic markers (microsatellites), we show that adaptive divergence is often associated with the presence of two populations maintaining at least partial reproductive isolation in parapatry. Coalescent-based simulations further suggest that these populations have diverged with gene flow within a few thousand generations, although we cannot rule out the possibility of phases of allopatric divergence. Finally, we find striking variation among the eight lake-stream transitions in progress toward ecological speciation. This variation allows us to hypothesize that progress is generally promoted by strong divergent selection and limited dispersal across the habitat transitions. Our study thus makes a case for ecological speciation in a parapatric context, while also highlighting variation in the outcome.     

Quantitative genetic inheritance of morphological divergence in a lake-stream stickleback ecotype pair: implications for reproductive isolation

Ecological selection against hybrids between populations occupying different habitats might be an important component of reproductive isolation during the initial stages of speciation. The strength and directionality of this barrier to gene flow depends on the genetic architecture underlying divergence in ecologically relevant phenotypes. We here present line cross analyses of inheritance for two key foraging-related morphological traits involved in adaptive divergence between stickleback ecotypes residing parapatrically in lake and stream habitats within the Misty Lake watershed (Vancouver Island, Canada). One main finding is the striking genetic dominance of the lake phenotype for body depth. Selection associated with this phenotype against first- and later-generation hybrids should therefore be asymmetric, hindering introgression from the lake to the stream population but not vice versa. Another main finding is that divergence in gill raker number is inherited additively and should therefore contribute symmetrically to reproductive isolation. Our study suggests that traits involved in adaptation might contribute to reproductive isolation qualitatively differently, depending on their mode of inheritance.     

Genetic evidence for ecological divergence in kokanee salmon

The evolution of locally adapted phenotypes among populations that experience divergent selective pressures is a central mechanism for generating and maintaining biodiversity. Recently, the advent of high‐throughput DNA sequencing technology has provided tools for investigating the genetic basis of this process in natural populations of nonmodel organisms. Kokanee, the freshwater form of sockeye salmon (Oncorhynchus nerka), occurs as two reproductive ecotypes, which differ in spawning habitat (tributaries vs. shorelines); however, outside of the spawning season the two ecotypes co‐occur in many lakes and lack diagnostic morphological characteristics. We used restriction site‐associated DNA (RAD) sequencing to identify 6145 SNPs and genotype kokanee from multiple spawning sites in Okanagan Lake (British Columbia, Canada). Outlier tests revealed 18 loci putatively under divergent selection between ecotypes, all of which exhibited temporally stable allele frequencies within ecotypes. Six outliers were annotated to sequences in the NCBI database, two of which matched genes associated with early development. There was no evidence for neutral genetic differentiation; however, outlier loci demonstrated significant structure with respect to ecotype and had high assignment accuracy in mixed composition simulations. The absence of neutral structure combined with a small number of highly divergent outlier loci is consistent with theoretical predictions for the early stages of ecological divergence. These outlier loci were then applied to a realistic fisheries scenario in which additional RAD sequencing was used to genotype kokanee collected by trawl in Okanagan Lake, providing preliminary evidence that this approach may be an effective tool for conservation and management.     

Recurrent adaptation in a low‐dispersal trait

The study of natural populations from contrasting environments has greatly enhanced our understanding of ecological‐dependent selection, adaptation and speciation. Cases of parallel evolution in particular have facilitated the study of the molecular and genetic basis of adaptive variation. This includes the type and number of genes underlying adaptive traits, as well as the extent to which these genes are exchanged among populations and contribute repeatedly to parallel evolution. Yet, surprisingly few studies provide a comprehensive view on the evolutionary history of adaptive traits from mutation to widespread adaptation. When did key mutations arise, how did they increase in frequency, and how did they spread? In this issue of Molecular Ecology, Van Belleghem et al. ( 2015 ) reconstruct the evolutionary history of a gene associated with wing size in the salt marsh beetle Pogonus chalceus. Screening the entire distribution range of this species, they found a single origin for the allele associated with the short‐winged ecotype. This allele seemingly evolved in an isolated population and rapidly introgressed into other populations. These findings suggest that the adaptive genetic variation found in sympatric short‐ and long‐winged populations has an allopatric origin, confirming that allopatric phases may be important at early stages of speciation.