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.     

Environmental effects on molecular and phenotypic variation in populations of Eruca sativa across a steep climatic gradient

In Israel Eruca sativa has a geographically narrow distribution across a steep climatic gradient that ranges from mesic Mediterranean to hot desert environments. These conditions offer an opportunity to study the influence of the environment on intraspecific genetic variation. For this, we combined an analysis of neutral genetic markers with a phenotypic evaluation in common‐garden experiments, and environmental characterization of populations that included climatic and edaphic parameters, as well as geographic distribution. A Bayesian clustering of individuals from nine representative populations based on amplified fragment length polymorphism (AFLP) divided the populations into a southern and a northern geographic cluster, with one admixed population at the geographic border between them. Linear mixed models, with cluster added as a grouping factor, revealed no clear effects of environment or geography on genetic distances, but this may be due to a strong association of geography and environment with genetic clusters. However, environmental factors accounted for part of the phenotypic variation observed in the common‐garden experiments. In addition, candidate loci for selection were identified by association with environmental parameters and by two outlier methods. One locus, identified by all three methods, also showed an association with trichome density and herbivore damage, in net‐house and field experiments, respectively. Accordingly, we propose that because trichomes are directly linked to defense against both herbivores and excess radiation, they could potentially be related to adaptive variation in these populations. These results demonstrate the value of combining environmental and phenotypic data with a detailed genetic survey when studying adaptation in plant populations.     

Ecotypes of an ecologically dominant prairie grass (Andropogon gerardii) exhibit genetic divergence across the U.S. Midwest grasslands' environmental gradient

Big bluestem (Andropogon gerardii) is an ecologically dominant grass with wide distribution across the environmental gradient of U.S. Midwest grasslands. This system offers an ideal natural laboratory to study population divergence and adaptation in spatially varying climates. Objectives were to: (i) characterize neutral genetic diversity and structure within and among three regional ecotypes derived from 11 prairies across the U.S. Midwest environmental gradient, (ii) distinguish between the relative roles of isolation by distance (IBD) vs. isolation by environment (IBE) on ecotype divergence, (iii) identify outlier loci under selection and (iv) assess the association between outlier loci and climate. Using two primer sets, we genotyped 378 plants at 384 polymorphic AFLP loci across regional ecotypes from central and eastern Kansas and Illinois. Neighbour‐joining tree and PCoA revealed strong genetic differentiation between Kansas and Illinois ecotypes, which was better explained by IBE than IBD. We found high genetic variability within prairies (80%) and even fragmented Illinois prairies, surprisingly, contained high within‐prairie genetic diversity (92%). Using Bayenv 2, 14 top‐ranked outlier loci among ecotypes were associated with temperature and precipitation variables. Six of seven BayeScan F_ST outliers were in common with Bayenv 2 outliers. High genetic diversity may enable big bluestem populations to better withstand changing climates; however, population divergence supports the use of local ecotypes in grassland restoration. Knowledge of genetic variation in this ecological dominant and other grassland species will be critical to understanding grassland response and restoration challenges in the face of a changing climate.     

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

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

Outlier loci highlight the direction of introgression in oaks

Loci considered to be under selection are generally avoided in attempts to infer past demographic processes as they do not fit neutral model assumptions. However, opportunities to better reconstruct some aspects of past demography might thus be missed. Here we examined genetic differentiation between two sympatric European oak species with contrasting ecological dynamics (Quercus robur and Quercus petraea) with both outlier (i.e. loci possibly affected by divergent selection between species or by hitchhiking effects with genomic regions under selection) and nonoutlier loci. We sampled 855 individuals in six mixed forests in France and genotyped them with a set of 262 SNPs enriched with markers showing high interspecific differentiation, resulting in accurate species delimitation. We identified between 13 and 74 interspecific outlier loci, depending on the coalescent simulation models and parameters used. Greater genetic diversity was predicted in Q. petraea (a late‐successional species) than in Q. robur (an early successional species) as introgression should theoretically occur predominantly from the resident species to the invading species. Remarkably, this prediction was verified with outlier loci but not with nonoutlier loci. We suggest that the lower effective interspecific gene flow at loci showing high interspecific divergence has better preserved the signal of past asymmetric introgression towards Q. petraea caused by the species' contrasting dynamics. Using markers under selection to reconstruct past demographic processes could therefore have broader potential than generally recognized.     

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.     

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

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

Investigating the extent of parallelism in morphological and genomic divergence among lake trout ecotypes in Lake Superior

Understanding the emergence of species through the process of ecological speciation is a central question in evolutionary biology which also has implications for conservation and management. Lake trout (Salvelinus namaycush) is renowned for the occurrence of different ecotypes linked to resource and habitat use throughout North America. We aimed to unravel the fine genetic structure of the four lake trout ecotypes in Lake Superior. A total of 486 individuals from four sites were genotyped at 6822 filtered SNPs using RADseq technology. Our results revealed different extent of morphological and genetic differentiation within the different sites. Overall, genetic differentiation was weak but significant and was on average three times higher between sites (mean F_ST = 0.016) than between ecotypes within sites (mean F_ST = 0.005) indicating higher level of gene flow or a more recent shared ancestor between ecotypes within each site than between populations of the same ecotype. Evidence of divergent selection was also found between ecotypes and/or in association with morphological variation. Outlier loci found in genes related to lipid metabolism and visual acuity were of particular interest in this context of ecotypic divergence. However, we did not find clear indication of parallelism at the genomic level, despite the presence of phenotypic parallelism among some ecotypes from different sampling sites. Overall, the occurrence of different levels of both genomic and phenotypic differentiation between ecotypes within each site with several differentiated loci linked to relevant biological functions supports the presence of a continuum of divergence in lake trout.     

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

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

Regional environmental pressure influences population differentiation in turbot (Scophthalmus maximus)

Unravelling the factors shaping the genetic structure of mobile marine species is challenging due to the high potential for gene flow. However, genetic inference can be greatly enhanced by increasing the genomic, geographical or environmental resolution of population genetic studies. Here, we investigated the population structure of turbot (Scophthalmus maximus) by screening 17 random and gene‐linked markers in 999 individuals at 290 geographical locations throughout the northeast Atlantic Ocean. A seascape genetics approach with the inclusion of high‐resolution oceanographical data was used to quantify the association of genetic variation with spatial, temporal and environmental parameters. Neutral loci identified three subgroups: an Atlantic group, a Baltic Sea group and one on the Irish Shelf. The inclusion of loci putatively under selection suggested an additional break in the North Sea, subdividing southern from northern Atlantic individuals. Environmental and spatial seascape variables correlated marginally with neutral genetic variation, but explained significant proportions (respectively, 8.7% and 10.3%) of adaptive genetic variation. Environmental variables associated with outlier allele frequencies included salinity, temperature, bottom shear stress, dissolved oxygen concentration and depth of the pycnocline. Furthermore, levels of explained adaptive genetic variation differed markedly between basins (3% vs. 12% in the North and Baltic Sea, respectively). We suggest that stable environmental selection pressure contributes to relatively strong local adaptation in the Baltic Sea. Our seascape genetic approach using a large number of sampling locations and associated oceanographical data proved useful for the identification of population units as the basis of management decisions.     

Natural Selection and Neutral Evolution Jointly Drive Population Divergence between Alpine and Lowland Ecotypes of the Allopolyploid Plant Anemone multifida (Ranunculaceae)

Population differentiation can be driven in large part by natural selection, but selectively neutral evolution can play a prominent role in shaping patters of population divergence. The decomposition of the evolutionary history of populations into the relative effects of natural selection and selectively neutral evolution enables an understanding of the causes of population divergence and adaptation. In this study, we examined heterogeneous genomic divergence between alpine and lowland ecotypes of the allopolyploid plant, Anemone multifida. Using peak height and dominant AFLP data, we quantified population differentiation at non-outlier (neutral) and outlier loci to determine the potential contribution of natural selection and selectively neutral evolution to population divergence. We found 13 candidate loci, corresponding to 2.7% of loci, with signatures of divergent natural selection between alpine and lowland populations and between alpine populations (Fst  = 0.074–0.445 at outlier loci), but neutral population differentiation was also evident between alpine populations (F_ST  = 0.041–0.095 at neutral loci). By examining population structure at both neutral and outlier loci, we determined that the combined effects of selection and neutral evolution are associated with the divergence of alpine populations, which may be linked to extreme abiotic conditions and isolation between alpine sites. The presence of outlier levels of genetic variation in structured populations underscores the importance of separately analyzing neutral and outlier loci to infer the relative role of divergent natural selection and neutral evolution in population divergence.     

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.     

Genetic component of flammability variation in a Mediterranean shrub

Recurrent fires impose a strong selection pressure in many ecosystems worldwide. In such ecosystems, plant flammability is of paramount importance because it enhances population persistence, particularly in non‐resprouting species. Indeed, there is evidence of phenotypic divergence of flammability under different fire regimes. Our general hypothesis is that flammability‐enhancing traits are adaptive; here, we test whether they have a genetic component. To test this hypothesis, we used the postfire obligate seeder Ulex parviflorus from sites historically exposed to different fire recurrence. We associated molecular variation in potentially adaptive loci detected with a genomic scan (using AFLP markers) with individual phenotypic variability in flammability across fire regimes. We found that at least 42% of the phenotypic variation in flammability was explained by the genetic divergence in a subset of AFLP loci. In spite of generalized gene flow, the genetic variability was structured by differences in fire recurrence. Our results provide the first field evidence supporting that traits enhancing plant flammability have a genetic component and thus can be responding to natural selection driven by fire. These results highlight the importance of flammability as an adaptive trait in fire‐prone ecosystems.     

Narrow hybrid zones in spite of very low population differentiation in neutral markers in an island bird species complex

Patterns of phenotypic and genic frequencies across hybrid zones provide insight into the origin and evolution of reproductive isolation. The Reunion grey white‐eye, Zosterops borbonicus, exhibits parapatrically distributed plumage colour forms across the lowlands of the small volcanic island of Reunion (Mascarene archipelago). These forms meet and hybridize in regions that are natural barriers to dispersal (rivers, lava fields). Here, we investigated the relationship among patterns of differentiation at neutral genetic (microsatellite) markers, phenotypic traits (morphology and plumage colour) and niche characteristics across three independent hybrid zones. Patterns of phenotypic divergence revealed that these hybrid zones are among the narrowest ever documented in birds. However, the levels of phenotypic divergence stand in stark contrast to the lack of clear population neutral genetic structure between forms. The position of the hybrid zones coincides with different natural physical barriers, yet is not associated with steep changes in vegetation and related climatic variables, and major habitat transitions are shifted from these locations by at least 18 km. This suggests that the hybrid zones are stabilized over natural dispersal barriers, independently of environmental boundaries, and are not associated with niche divergence. A striking feature of these hybrid zones is the very low levels of genetic differentiation in neutral markers between forms, suggesting that phenotypic divergence has a narrow genetic basis and may reflect recent divergence at a few linked genes under strong selection, with a possible role for assortative mating in keeping these forms apart.     

Evidence for nonallopatric speciation among closely related sympatric Heliotropium species in the Atacama Desert

The genetic structure of populations of closely related, sympatric species may hold the signature of the geographical mode of the speciation process. In fully allopatric speciation, it is expected that genetic differentiation between species is homogeneously distributed across the genome. In nonallopatric speciation, the genomes may remain undifferentiated to a large extent. In this article, we analyzed the genetic structure of five sympatric species from the plant genus Heliotropium in the Atacama Desert. We used amplified fragment length polymorphisms (AFLPs) to characterize the genetic structure of these species and evaluate their genetic differentiation as well as the number of loci subject to positive selection using divergence outlier analysis (DOA). The five species form distinguishable groups in the genetic space, with zones of overlap, indicating that they are possibly not completely isolated. Among‐species differentiation accounts for 35% of the total genetic differentiation (F_ST = 0.35), and F_ST between species pairs is positively correlated with phylogenetic distance. DOA suggests that few loci are subject to positive selection, which is in line with a scenario of nonallopatric speciation. These results support the idea that sympatric species of Heliotropium sect. Cochranea are under an ongoing speciation process, characterized by a fluctuation of population ranges in response to pulses of arid and humid periods during Quaternary times.     

Heterogeneous genomic differentiation between walking-stick ecotypes: "isolation by adaptation" and multiple roles for divergent selection

Genetic differentiation can be highly variable across the genome. For example, loci under divergent selection and those tightly linked to them may exhibit elevated differentiation compared to neutral regions. These represent "outlier loci" whose differentiation exceeds neutral expectations. Adaptive divergence can also increase genome-wide differentiation by promoting general barriers to neutral gene flow, thereby facilitating genomic divergence via genetic drift. This latter process can yield a positive correlation between adaptive phenotypic divergence and neutral genetic differentiation (described here as "isolation-by-adaptation"). Here, we examine both these processes by combining an AFLP genome scan of two host plant ecotypes of Timema cristinae walking-sticks with existing data on adaptive phenotypic divergence and ecological speciation in these insects. We found that about 8% of loci are outliers in multiple population comparisons. Replicated comparisons between population-pairs using the same versus different host species revealed that 1-2% of loci are subject to host-related selection specifically. Locus-specific analyses revealed that up to 10% of putatively neutral (nonoutlier) AFLP loci exhibit significant isolation-by-adaptation. Our results suggest that selection may affect differentiation directly, via linkage, or by facilitating genetic drift. They thus illustrate the varied and sometimes nonintuitive contributions of selection to heterogeneous genomic differentiation.     

An ecological genetic delineation of local seed‐source provenance for ecological restoration

An increasingly important practical application of the analysis of spatial genetic structure within plant species is to help define the extent of local provenance seed collection zones that minimize negative impacts in ecological restoration programs. Here, we derive seed sourcing guidelines from a novel range‐wide assessment of spatial genetic structure of 24 populations of Banksia menziesii (Proteaceae), a widely distributed Western Australian tree of significance in local ecological restoration programs. An analysis of molecular variance (AMOVA) of 100 amplified fragment length polymorphism (AFLP) markers revealed significant genetic differentiation among populations (Φ_PT = 0.18). Pairwise population genetic dissimilarity was correlated with geographic distance, but not environmental distance derived from 15 climate variables, suggesting overall neutrality of these markers with regard to these climate variables. Nevertheless, Bayesian outlier analysis identified four markers potentially under selection, although these were not correlated with the climate variables. We calculated a global R‐statistic using analysis of similarities (ANOSIM) to test the statistical significance of population differentiation and to infer a threshold seed collection zone distance of ~60 km (all markers) and 100 km (outlier markers) when genetic distance was regressed against geographic distance. Population pairs separated by >60 km were, on average, twice as likely to be significantly genetically differentiated than population pairs separated by <60 km, suggesting that habitat‐matched sites within a 30‐km radius around a restoration site genetically defines a local provenance seed collection zone for B. menziesii. Our approach is a novel probability‐based practical solution for the delineation of a local seed collection zone to minimize negative genetic impacts in ecological restoration.     

Population genomics of parallel phenotypic evolution in stickleback across stream-lake ecological transitions

Understanding the genetics of adaptation is a central focus in evolutionary biology. Here, we use a population genomics approach to examine striking parallel morphological divergences of parapatric stream-lake ecotypes of threespine stickleback fish in three watersheds on the Haida Gwaii archipelago, western Canada. Genome-wide variation at greater than 1000 single nucleotide polymorphism loci indicate separate origin of giant lake and small-bodied stream fish within each watershed (mean F(ST) between watersheds = 0.244 and within = 0.114). Genome scans within watersheds identified a total of 21 genomic regions that are highly differentiated between ecotypes and are probably subject to directional selection. Most outliers were watershed-specific, but genomic regions undergoing parallel genetic changes in multiple watersheds were also identified. Interestingly, several of the stream-lake outlier regions match those previously identified in marine-freshwater and benthic-limnetic genome scans, indicating reuse of the same genetic loci in different adaptive scenarios. We also identified multiple new outlier loci, which may contribute to unique aspects of differentiation in stream-lake environments. Overall, our data emphasize the important role of ecological boundaries in driving both local and broadly occurring parallel genetic changes during adaptation.