Genomic islands of divergence linked to ecotypic variation in sockeye salmon

Regions of the genome displaying elevated differentiation (genomic islands of divergence) are thought to play an important role in local adaptation, especially in populations experiencing high gene flow. However, the characteristics of these islands as well as the functional significance of genes located within them remain largely unknown. Here, we used data from thousands of SNPs aligned to a linkage map to investigate genomic islands of divergence in three ecotypes of sockeye salmon (Oncorhynchus nerka) from a single drainage in southwestern Alaska. We found ten islands displaying high differentiation among ecotypes. Conversely, neutral structure observed throughout the rest of the genome was low and not partitioned by ecotype. One island on linkage group So13 was particularly large and contained six SNPs with F_ST > 0.14 (average F_ST of neutral SNPs = 0.01). Functional annotation revealed that the peak of this island contained a nonsynonymous mutation in a gene involved in growth in other species (TULP4). The islands that we discovered were relatively small (80–402 Kb), loci found in islands did not show reduced levels of diversity, and loci in islands displayed slightly elevated linkage disequilibrium. These attributes suggest that the islands discovered here were likely generated by divergence hitchhiking; however, we cannot rule out the possibility that other mechanisms may have produced them. Our results suggest that islands of divergence serve an important role in local adaptation with gene flow and represent a significant advance towards understanding the genetic basis of ecotypic differentiation.     

Single-nucleotide polymorphism in populations of sockeye salmon Oncorhynchus nerka from Kamchatka Peninsula

The genetic variability of 45 single-nucleotide polymorphism loci was examined in the four largest wild populations of sockeye salmon Oncorhynchus nerka from drainages of the Asian coast of the Pacific Ocean (Eastern and Western Kamchatka). It was demonstrated that sockeye salmon from the Palana River were considerably different from all other populations examined. The most probable explanation of the observed differences is the suggestion on possible demographic events in the history of this population associated with the decrease in its effective number. To study the origin, colonization patterns, and evolution of Asian sockeye salmon, as well as to resolve some of the applied tasks, like population assignment and genetic identification, a differential approach to SNP-marker selection was suggested. Adaptively important loci that evolve under the pressure of balancing (stabilizing) selection were identified, owing to this fact the number of loci that provide the baseline classification error rates in the population assignment tests was reduced to 30. It was demonstrated that SNPs located in the MHC2 and GPH genes were affected by diversifying selection. Procedures for selecting single-nucleotide polymorphisms for phylogenetic studies of Asian sockeye salmon were suggested. Using principal-component analysis, 17 loci that adequately reproduce genetic differentiation within and among the regions of the origin of Kamchatka sockeye salmon were selected.     

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.     

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.     

The adaptive genomic landscape of beak morphology in Darwin's finches

Beak shape in Darwin's ground finches (Geospiza) is emblematic of natural selection and adaptive radiation, yet our understanding of the genetic basis of beak shape variation, and thus the genetic target of natural selection, is still evolving. Here we reveal the genomic architecture of beak shape variation using genomewide comparisons of four closely related and hybridizing species across 13 islands subject to parallel natural selection. Pairwise contrasts among species were used to identify a large number of genomic loci that are consistently related to species differences across a complex landscape. These loci are associated with hundreds of genes that have enriched GO categories significantly associated with development. One genomic region of particular interest is a section of Chromosome 1A with many candidate genes and increased linkage. The distinct, pointed beak shape of the cactus finch is linked to an excess of intermediate frequency alleles and increased heterozygosity in significant SNPs, but not across the rest of the genome. Alleles associated with pointier beaks among species were associated with pointier‐beaked populations within each species, thus establishing a common basis for natural selection, species divergence and adaptive radiation. The adaptive genomic landscape for Darwin's finches mirrors theoretical expectations based on morphological variation. The implication that a large number of genes are actively maintained to facilitate beak variation across parallel populations with documented interspecies admixture challenges our understanding of evolutionary processes in the wild.     

Patterns of selection and allele diversity of class I and class II major histocompatibility loci across the species range of sockeye salmon (Oncorhynchus nerka)

The major histocompatibility complex (MHC), an important component of the vertebrate immune system, provides an important suite of genes to examine the role of genetic diversity at non‐neutral loci for population persistence. We contrasted patterns of diversity at the two classical MHC loci in sockeye salmon (Oncorhynchus nerka), MHC class I (UBA) and MHC class II (DAB), and neutral microsatellite loci across 70 populations spanning the species range from Washington State to Japan. There was no correlation in allelic richness or heterozygosity between MHC loci or between MHC loci and microsatellites. The two unlinked MHC loci may be responding to different selective pressures; the distribution of F_ST values for the two loci was uncorrelated, and evidence for both balancing and directional selection on alleles and lineages of DAB and UBA was observed in populations throughout the species range but rarely on both loci within a population. These results suggest that fluctuating selection has resulted in the divergence of MHC loci in contemporary populations.     

Signatures of natural selection among lineages and habitats in Oncorhynchus mykiss

Recent advances in molecular interrogation techniques now allow unprecedented genomic inference about the role of adaptive genetic divergence in wild populations. We used high-throughput genotyping to screen a genome-wide panel of 276 single nucleotide polymorphisms (SNPs) for the economically and culturally important salmonid Oncorhynchus mykiss. Samples included 805 individuals from 11 anadromous and resident populations from the northwestern United States and British Columbia, and represented two major lineages including paired populations of each life history within single drainages of each lineage. Overall patterns of variation affirmed clear distinctions between lineages and in most instances, isolation by distance within them. Evidence for divergent selection at eight candidate loci included significant landscape correlations, particularly with temperature. High diversity of two nonsynonymous mutations within the peptide-binding region of the major histocompatibility complex (MHC) class II (DAB) gene provided signatures of balancing selection. Weak signals for potential selection between sympatric resident and anadromous populations were revealed from genome scans and allele frequency comparisons. Our results suggest an important adaptive role for immune-related functions and present a large genomic resource for future studies.     

Persistent inter‐ and intraspecific gene exchange within a parallel radiation of caterpillar hunter beetles (Calosoma sp.) from the Galápagos

When environmental gradients are repeated on different islands within an archipelago, similar selection pressures may act within each island, resulting in the repeated occurrence of ecologically similar species on each island. The evolution of ecotypes within such radiations may either result from dispersal, that is each ecotype evolved once and dispersed to different islands where it colonized its habitat, or through repeated and parallel speciation within each island. However, it remains poorly understood how gene flow during the divergence process may shape such patterns. In the Galápagos islands, three phenotypically similar species of the beetle genus Calosoma occur at higher elevations of different islands, while lowlands are occupied by a fourth species. By genotyping all major populations within this radiation for two nuclear and three mitochondrial gene fragments and seven microsatellite markers, we found strong support that the oldest divergence separates the highland species of the oldest island from the remaining species. Despite their morphological distinctness, highland species of the remaining islands were genetically closely related to the lowland population on each island and within the same magnitude as lowland populations sampled at different islands. Repeated evolution of highland ecotypes out of the lowland species appears the most likely scenario and estimates of geneflow rates revealed extensive admixture among ecotypes within islands, as well as between islands. These findings indicate that gene exchange among the different populations and species may have shaped the phylogenetic relationships and the repeated evolution of these ecotypes.     

Selection from parasites favours immunogenetic diversity but not divergence among locally adapted host populations

The unprecedented polymorphism in the major histocompatibility complex (MHC) genes is thought to be maintained by balancing selection from parasites. However, do parasites also drive divergence at MHC loci between host populations, or do the effects of balancing selection maintain similarities among populations? We examined MHC variation in populations of the livebearing fish Poecilia mexicana and characterized their parasite communities. Poecilia mexicana populations in the Cueva del Azufre system are locally adapted to darkness and the presence of toxic hydrogen sulphide, representing highly divergent ecotypes or incipient species. Parasite communities differed significantly across populations, and populations with higher parasite loads had higher levels of diversity at class II MHC genes. However, despite different parasite communities, marked divergence in adaptive traits and in neutral genetic markers, we found MHC alleles to be remarkably similar among host populations. Our findings indicate that balancing selection from parasites maintains immunogenetic diversity of hosts, but this process does not promote MHC divergence in this system. On the contrary, we suggest that balancing selection on immunogenetic loci may outweigh divergent selection causing divergence, thereby hindering host divergence and speciation. Our findings support the hypothesis that balancing selection maintains MHC similarities among lineages during and after speciation (trans‐species evolution).     

Gene flow, effective population size and selection at major histocompatibility complex genes: brown trout in the Hardanger Fjord, Norway

Brown trout populations in the Hardanger Fjord, Norway, have declined drastically due to increased exposure to salmon lice from salmonid aquaculture. We studied contemporary samples from seven populations and historical samples (1972 and 1983) from the two largest populations, one of which has declined drastically whereas the other remains stable. We analysed 11 microsatellite loci, including one tightly linked to the UBA gene of the major histocompatibility class I complex (MHC) and another locus linked to the TAP2A gene, also associated with MHC. The results revealed asymmetric gene flow from the two largest populations to the other, smaller populations. This has important conservation implications, and we predict that possible future population recoveries will be mediated primarily by the remaining large population. Tests for selection suggested diversifying selection at UBA, whereas evidence was inconclusive for TAP2A. There was no evidence for temporally fluctuating selection. We assessed the distribution of adaptive divergence among populations. The results showed the most pronounced footprints of selection between the two largest populations subject to the least immigration. We suggest that asymmetric gene flow has an important influence on adaptive divergence and constrains local adaptive responses in the smaller populations. Even though UBA alleles may not affect salmon louse resistance, the results bear evidence of adaptive divergence among populations at immune system genes. This suggests that similar genetic differences could exist at salmon louse resistance loci, thus rendering it a realistic scenario that differential population declines could reflect differences in adaptive variation.     

Signals of heterogeneous selection at an MHC locus in geographically proximate ecotypes of sockeye salmon

The genes of the major histocompatibility complex (MHC) are an important component of the vertebrate immune system and can provide insights into the role of pathogen‐mediated selection in wild populations. Here, we examined variation at the MHC class II peptide‐binding region in 27 populations of sockeye salmon (Oncorhynchus nerka), distributed among three distinct spawning ecotypes, from a complex of interconnected rivers and lakes in south‐western Alaska. We also obtained genotypes from 90 putatively neutral single nucleotide polymorphisms for each population to compare the relative roles of demography and selection in shaping the observed MHC variation. We found that MHC divergence was generally partitioned by spawning ecotype (lake beaches, rivers and streams) and was 30 times greater than variation at neutral markers. Additionally, we observed substantial differences in modes of selection and diversity among ecotypes, with beach populations displaying higher levels of directional selection and lower MHC diversity than the other two ecotypes. Finally, the level of MHC differentiation in our study system was comparable to that observed over much larger geographic ranges, suggesting that MHC variation does not necessarily increase with increasing spatial scale and may instead be driven by fine‐scale differences in pathogen communities or pathogen virulence. The low levels of neutral structure and spatial proximity of populations in our study system indicate that MHC differentiation can be maintained through strong selective pressure even when ample opportunities for gene flow exist.     

Population structure of lake‐type and river‐type sockeye salmon in transboundary rivers of northern British Columbia

The population structure of 'lake‐type' and 'river‐type' sockeye salmon Oncorhynchus nerka , primarily in transboundary rivers in northern British Columbia, was examined with a survey of microsatellite variation. Variation at 14 microsatellite loci was surveyed from c . 3000 lake‐type and 3200 river‐type sockeye salmon from 47 populations in six river drainages in British Columbia. The mean F _ST for the 14 microsatellite loci and 47 populations was 0·068, and 0·034 over all river‐type populations. River‐type sockeye salmon were more genetically diverse than lake‐type sockeye salmon, with expected heterozygosity of river‐type sockeye salmon 0·72 and with an average 12·7 alleles observed per locus, whereas expected heterozygosity of lake‐type sockeye salmon was 0·65 with and average 10·5 alleles observed per locus. River drainage of origin was a significant unit of population structure. There was clear evidence of genetic differentiation among river‐type populations of sockeye salmon from different drainages over a broad geographic range in British Columbia.     

Sequence analysis of MHC class I α2 from sockeye salmon (Oncorhynchus nerka)

Most studies assessing adaptive MHC diversity in salmon populations have focused on the classical class II DAB or DAA loci, as these have been most amenable to single PCR amplifications due to their relatively low level of sequence divergence. Herein, we report the characterization of the classical class I UBA α2 locus based on collections taken throughout the species range of sockeye salmon (Oncorhynchus nerka). Through use of multiple lineage-specific primer sets, denaturing gradient gel electrophoresis and sequencing, we identified thirty-four alleles from three highly divergent lineages. Sequence identity between lineages ranged from 30.0% to 56.8% but was relatively high within lineages. Allelic identity within the antigen recognition site (ARS) was greater than for the longer sequence. Global positive selection on UBA was seen at the sequence level (dN:dS = 1.012) with four codons under positive selection and 12 codons under negative selection.     

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.     

Site-specific genetic divergence in parallel hybrid zones suggests nonallopatric evolution of reproductive barriers

The evolution of reproductive isolation in the presence of gene flow is supported by theoretical models but rarely by data. Empirical support might be gained from studies of parallel hybrid zones between interbreeding taxa. We analysed gene flow over two hybrid zones separating ecotypes of Littorina saxatilis to test the expectation that neutral genetic markers will show site-specific differences if barriers have evolved in situ. Distinct ecotypes found in contrasting shore habitats are separated by divergent selection and poor dispersal, but hybrid zones appear between them. Swedish islands formed by postglacial uplift 5000 years ago provide opportunities to assess genetic structure in a recently evolved system. Each island houses a discrete population containing subpopulations of different ecotypes. Hybrid zones between ecotypes may be a product of ecological divergence occurring on each island or a consequence of secondary overlap of ecotypes of allopatric origin that have spread among the islands. We used six microsatellite loci to assess gene flow and genetic profiles of hybrid zones on two islands. We found reduced gene flow over both hybrid zones, indicating the presence of local reproductive barriers between ecotypes. Nevertheless, subpopulations of different ecotypes from the same island were genetically more similar to each other than were subpopulations of the same ecotype from different islands. Moreover, neutral genetic traits separating the two ecotypes across hybrid zones were site-specific. This supports a scenario of in situ origin of ecotypes by ecological divergence and nonallopatric evolution of reproductive barriers.     

Divergent and linked selection shape patterns of genomic differentiation between European and North American Atlantic salmon (Salmo salar)

As populations diverge many processes can shape genomic patterns of differentiation. Regions of high differentiation can arise due to divergent selection acting on selected loci, genetic hitchhiking of nearby loci, or through repeated selection against deleterious alleles (linked background selection); this divergence may then be further elevated in regions of reduced recombination. Atlantic salmon (Salmo salar) from Europe and North America diverged >600,000 years ago and despite some evidence of secondary contact, the majority of genetic data indicate substantial divergence between lineages. This deep divergence with potential gene flow provides an opportunity to investigate the role of different mechanisms that shape the genomic landscape during early speciation. Here, using 184,295 single nucleotide polymorphisms (SNPs) and 80 populations, we investigate the genomic landscape of differentiation across the Atlantic Ocean with a focus on highly differentiated regions and the processes shaping them. We found evidence of high (mean F_ST = 0.26) and heterogeneous genomic differentiation between continents. Genomic regions associated with high trans‐Atlantic differentiation ranged in size from single loci (SNPs) within important genes to large regions (1–3 Mbp ) on four chromosomes (Ssa06, Ssa13, Ssa16 and Ssa19). These regions showed signatures consistent with selection, including high linkage disequilibrium, despite no significant reduction in recombination. Genes and functional enrichment of processes associated with differentiated regions may highlight continental differences in ocean navigation and parasite resistance. Our results provide insight into potential mechanisms underlying differences between continents, and evidence of near‐fixed and potentially adaptive trans‐Atlantic differences concurrent with a background of high genome‐wide differentiation supports subspecies designation in Atlantic salmon.     

Single nucleotide polymorphisms unravel hierarchical divergence and signatures of selection among Alaskan sockeye salmon (Oncorhynchus nerka) populations

Background  

Disentangling the roles of geography and ecology driving population divergence and distinguishing adaptive from neutral evolution at the molecular level have been common goals among evolutionary and conservation biologists. Using single nucleotide polymorphism (SNP) multilocus genotypes for 31 sockeye salmon (Oncorhynchus nerka) populations from the Kvichak River, Alaska, we assessed the relative roles of geography (discrete boundaries or continuous distance) and ecology (spawning habitat and timing) driving genetic divergence in this species at varying spatial scales within the drainage. We also evaluated two outlier detection methods to characterize candidate SNPs responding to environmental selection, emphasizing which mechanism(s) may maintain the genetic variation of outlier loci.     

Recent local adaptation of sockeye salmon to glacial spawning habitats

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

LANDSCAPE GENOMICS IN ATLANTIC SALMON (SALMO SALAR): SEARCHING FOR GENE–ENVIRONMENT INTERACTIONS DRIVING LOCAL ADAPTATION

A growing number of studies are examining the factors driving historical and contemporary evolution in wild populations. By combining surveys of genomic variation with a comprehensive assessment of environmental parameters, such studies can increase our understanding of the genomic and geographical extent of local adaptation in wild populations. We used a large‐scale landscape genomics approach to examine adaptive and neutral differentiation across 54 North American populations of Atlantic salmon representing seven previously defined genetically distinct regional groups. Over 5500 genome‐wide single nucleotide polymorphisms were genotyped in 641 individuals and 28 bulk assays of 25 pooled individuals each. Genome scans, linkage map, and 49 environmental variables were combined to conduct an innovative landscape genomic analysis. Our results provide valuable insight into the links between environmental variation and both neutral and potentially adaptive genetic divergence. In particular, we identified markers potentially under divergent selection, as well as associated selective environmental factors and biological functions with the observed adaptive divergence. Multivariate landscape genetic analysis revealed strong associations of both genetic and environmental structures. We found an enrichment of growth‐related functions among outlier markers. Climate (temperature–precipitation) and geological characteristics were significantly associated with both potentially adaptive and neutral genetic divergence and should be considered as candidate loci involved in adaptation at the regional scale in Atlantic salmon. Hence, this study significantly contributes to the improvement of tools used in modern conservation and management schemes of Atlantic salmon wild populations.     

Integrating QTL mapping and genome scans towards the characterization of candidate loci under parallel selection in the lake whitefish (Coregonus clupeaformis)

As natural selection must act on underlying genetic variation, discovering the number and location of loci under the influence of selection is imperative towards understanding adaptive divergence in evolving populations. Studies employing genome scans have hypothesized that the action of divergent selection should reduce gene flow at the genomic locations implicated in adaptation and speciation among natural populations, yet once 'outlier' patterns of variation have been identified the function and role of such loci needs to be confirmed. We integrated adaptive QTL mapping and genomic scans among diverging sympatric pairs of the lake whitefish (Coregonus clupeaformis) species complex in order to test the hypothesis that differentiation between dwarf and normal ecotypes at growth-associated QTL was maintained by directional selection. We found evidence of significantly high levels of molecular divergence among eight growth QTL where two of the strongest candidate loci under the influence of directional selection exhibited parallel reductions of gene flow over multiple populations.