Genetic analysis of differentiation among breeding ponds reveals a candidate gene for local adaptation in Rana arvalis

One of the main questions in evolutionary and conservation biology is how geographical and environmental features of the landscape shape neutral and adaptive genetic variation in natural populations. The identification of genomic polymorphisms that account for adaptive variation can aid in finding candidate loci for local adaptation. Consequently, a comparison of spatial patterns in neutral markers and loci under selection may help disentangle the effects of gene flow, genetic drift and selection at the landscape scale. Many amphibians breed in wetlands, which differ in environmental conditions and in the degree of isolation, enhancing the potential for local adaptation. We used microsatellite markers to measure genetic differentiation among 17 local populations of Rana arvalis breeding in a network of wetlands. We found that locus RC08604 deviated from neutral expectations, suggesting that it is a good candidate for directional selection. We used a genetic network analysis to show that the allele distribution in this locus is correlated with habitat characteristics, whereas this was not the case at neutral markers that displayed a different allele distribution and population network in the study area. The graph approach illustrated the genomic heterogeneity (neutral loci vs. the candidate locus for directional selection) of gene exchange and genetic divergence among populations under directional selection. Limited gene flow between wetlands was only observed at the candidate genomic region under directional selection. RC08604 is partially located inside an up‐regulated thyroid‐hormone receptor (TRβ) gene coordinating the expression of other genes during metamorphosis and appears to be linked with variation in larval life‐history traits found among R. arvalis populations. We suggest that directional selection on genes coding larval life‐history traits is strong enough to maintain the divergence in these genomic regions, reducing the effective recombination of locally adapted alleles but not in other regions of the genome. Integrating this knowledge into conservation plans at the landscape scale will improve the design of management strategies to preserve adaptive genetic diversity in wetland networks.     

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.     

Chromosomal Fusions Facilitate Adaptation to Divergent Environments in Threespine Stickleback

Chromosomal fusions are hypothesized to facilitate adaptation to divergent environments, both by bringing together previously unlinked adaptive alleles and by creating regions of low recombination that facilitate the linkage of adaptive alleles; but, there is little empirical evidence to support this hypothesis. Here, we address this knowledge gap by studying threespine stickleback (Gasterosteus aculeatus), in which ancestral marine fish have repeatedly adapted to freshwater across the northern hemisphere. By comparing the threespine and ninespine stickleback (Pungitius pungitius) genomes to a de novo assembly of the fourspine stickleback (Apeltes quadracus) and an outgroup species, we find two chromosomal fusion events involving the same chromosomes have occurred independently in the threespine and ninespine stickleback lineages. On the fused chromosomes in threespine stickleback, we find an enrichment of quantitative trait loci underlying traits that contribute to marine versus freshwater adaptation. By comparing whole-genome sequences of freshwater and marine threespine stickleback populations, we also find an enrichment of regions under divergent selection on these two fused chromosomes. There is elevated genetic diversity within regions under selection in the freshwater population, consistent with a simulation study showing that gene flow can increase diversity in genomic regions associated with local adaptation and our demographic models showing gene flow between the marine and freshwater populations. Integrating our results with previous studies, we propose that these fusions created regions of low recombination that enabled the formation of adaptative clusters, thereby facilitating freshwater adaptation in the face of recurrent gene flow between marine and freshwater threespine sticklebacks.     

Widespread yet heterogeneous genomic divergence

Genetic differentiation during adaptive divergence and speciation is heterogeneous among genomic regions. Some regions can be highly differentiated between populations, for example, because they harbour genes under divergent selection or those causing reproductive isolation and thus are resistant to gene flow. Other regions might be homogenized by gene flow and thus weakly differentiated. Debates persist about the number of differentiated regions expected under divergence with gene flow, and their causes, size, and genomic distribution. In this issue of Molecular Ecology, a study of freshwater stickleback used next-generation sequencing to shed novel insight into these issues (Roesti et al. 2012). Many genomic regions distributed across the genome were strongly differentiated, indicating divergence with gene flow can involve a greater number of loci than often thought. Nonetheless, differentiation of some regions, such as those near the centre of chromosomes where recombination is reduced, was strongly accentuated over others. Thus, divergence was widespread yet highly heterogeneous across the genome. Moreover, different population pairs varied in patterns of differentiation, illustrating how genomic divergence builds up across stages of the speciation process. The study demonstrates how variation in different evolutionary processes, such as selection and recombination rate, can combine to result in similar genomic patterns. Future work could focus on teasing apart the contributions of different processes for causing differentiation, a task facilitated by experimental manipulations.     

Localizing FST outliers on a QTL map reveals evidence for large genomic regions of reduced gene exchange during speciation‐with‐gene‐flow

Populations that maintain phenotypic divergence in sympatry typically show a mosaic pattern of genomic divergence, requiring a corresponding mosaic of genomic isolation (reduced gene flow). However, mechanisms that could produce the genomic isolation required for divergence‐with‐gene‐flow have barely been explored, apart from the traditional localized effects of selection and reduced recombination near centromeres or inversions. By localizing F_ST outliers from a genome scan of wild pea aphid host races on a Quantitative Trait Locus (QTL) map of key traits, we test the hypothesis that between‐population recombination and gene exchange are reduced over large ‘divergence hitchhiking’ (DH) regions. As expected under divergence hitchhiking, our map confirms that QTL and divergent markers cluster together in multiple large genomic regions. Under divergence hitchhiking, the nonoutlier markers within these regions should show signs of reduced gene exchange relative to nonoutlier markers in genomic regions where ongoing gene flow is expected. We use this predicted difference among nonoutliers to perform a critical test of divergence hitchhiking. Results show that nonoutlier markers within clusters of F_ST outliers and QTL resolve the genetic population structure of the two host races nearly as well as the outliers themselves, while nonoutliers outside DH regions reveal no population structure, as expected if they experience more gene flow. These results provide clear evidence for divergence hitchhiking, a mechanism that may dramatically facilitate the process of speciation‐with‐gene‐flow. They also show the power of integrating genome scans with genetic analyses of the phenotypic traits involved in local adaptation and population divergence.     

A universal mechanism generating clusters of differentiated loci during divergence‐with‐migration

Genome‐wide patterns of genetic divergence reveal mechanisms of adaptation under gene flow. Empirical data show that divergence is mostly concentrated in narrow genomic regions. This pattern may arise because differentiated loci protect nearby mutations from gene flow, but recent theory suggests this mechanism is insufficient to explain the emergence of concentrated differentiation during biologically realistic timescales. Critically, earlier theory neglects an inevitable consequence of genetic drift: stochastic loss of local genomic divergence. Here, we demonstrate that the rate of stochastic loss of weak local differentiation increases with recombination distance to a strongly diverged locus and, above a critical recombination distance, local loss is faster than local “gain” of new differentiation. Under high migration and weak selection, this critical recombination distance is much smaller than the total recombination distance of the genomic region under selection. Consequently, divergence between populations increases by net gain of new differentiation within the critical recombination distance, resulting in tightly linked clusters of divergence. The mechanism responsible is the balance between stochastic loss and gain of weak local differentiation, a mechanism acting universally throughout the genome. Our results will help to explain empirical observations and lead to novel predictions regarding changes in genomic architectures during adaptive divergence.     

Genetic Barriers to Historical Gene Flow between Cryptic Species of Alpine Bumblebees Revealed by Comparative Population Genomics

Evidence is accumulating that gene flow commonly occurs between recently diverged species, despite the existence of barriers to gene flow in their genomes. However, we still know little about what regions of the genome become barriers to gene flow and how such barriers form. Here, we compare genetic differentiation across the genomes of bumblebee species living in sympatry and allopatry to reveal the potential impact of gene flow during species divergence and uncover genetic barrier loci. We first compared the genomes of the alpine bumblebee Bombus sylvicola and a previously unidentified sister species living in sympatry in the Rocky Mountains, revealing prominent islands of elevated genetic divergence in the genome that colocalize with centromeres and regions of low recombination. This same pattern is observed between the genomes of another pair of closely related species living in allopatry (B. bifarius and B. vancouverensis). Strikingly however, the genomic islands exhibit significantly elevated absolute divergence (d_XY) in the sympatric, but not the allopatric, comparison indicating that they contain loci that have acted as barriers to historical gene flow in sympatry. Our results suggest that intrinsic barriers to gene flow between species may often accumulate in regions of low recombination and near centromeres through processes such as genetic hitchhiking, and that divergence in these regions is accentuated in the presence of gene flow.     

Divergent selection and heterogeneous genomic divergence

Levels of genetic differentiation between populations can be highly variable across the genome, with divergent selection contributing to such heterogeneous genomic divergence. For example, loci under divergent selection and those tightly physically linked to them may exhibit stronger differentiation than neutral regions with weak or no linkage to such loci. Divergent selection can also increase genome‐wide neutral differentiation by reducing gene flow (e.g. by causing ecological speciation), thus promoting divergence via the stochastic effects of genetic drift. These consequences of divergent selection are being reported in recently accumulating studies that identify: (i) ‘outlier loci’ with higher levels of divergence than expected under neutrality, and (ii) a positive association between the degree of adaptive phenotypic divergence and levels of molecular genetic differentiation across population pairs [‘isolation by adaptation’ (IBA)]. The latter pattern arises because as adaptive divergence increases, gene flow is reduced (thereby promoting drift) and genetic hitchhiking increased. Here, we review and integrate these previously disconnected concepts and literatures. We find that studies generally report 5–10% of loci to be outliers. These selected regions were often dispersed across the genome, commonly exhibited replicated divergence across different population pairs, and could sometimes be associated with specific ecological variables. IBA was not infrequently observed, even at neutral loci putatively unlinked to those under divergent selection. Overall, we conclude that divergent selection makes diverse contributions to heterogeneous genomic divergence. Nonetheless, the number, size, and distribution of genomic regions affected by selection varied substantially among studies, leading us to discuss the potential role of divergent selection in the growth of regions of differentiation (i.e. genomic islands of divergence), a topic in need of future investigation.     

Adaptive divergence in the monkey flower Mimulus guttatus is maintained by a chromosomal inversion

Organisms exhibit an incredible diversity of life history strategies as adaptive responses to environmental variation. The establishment of novel life history strategies involves multilocus polymorphisms, which will be challenging to establish in the face of gene flow and recombination. Theory predicts that adaptive allelic combinations may be maintained and spread if they occur in genomic regions of reduced recombination, such as chromosomal inversion polymorphisms, yet empirical support for this prediction is lacking. Here, we use genomic data to investigate the evolution of divergent adaptive ecotypes of the yellow monkey flower Mimulus guttatus. We show that a large chromosomal inversion polymorphism is the major region of divergence between geographically widespread annual and perennial ecotypes. In contrast, ～40,000 single nucleotide polymorphisms in collinear regions of the genome show no signal of life history, revealing genomic patterns of diversity have been shaped by localized homogenizing gene flow and large‐scale Pleistocene range expansion. Our results provide evidence for an inversion capturing and protecting loci involved in local adaptation, while also explaining how adaptive divergence can occur with gene flow.     

Adaptive chromosomal divergence driven by mixed geographic mode of evolution

Chromosomal inversions are ubiquitous in nature and of great significance for understanding adaptation and speciation. Inversions were the first markers used to investigate the genetic structure of natural populations, leading to the concept of coadapted gene complexes and theories concerning founder effects and genetic drift in small populations. However, we still lack elements of a general theory accounting for the origins and distribution of inversions in nature. Here, we use computer simulations to show that a "mixed geographic mode" of evolution involving allopatric separation of populations followed by secondary contact and gene flow generates chromosomal divergence by natural selection under wider conditions than previous hypotheses. This occurs because inversions arising in allopatry contain a full complement of locally adapted genes. Once gene flow ensues, reduced recombination within inversions keeps these favorable genotypic combinations intact, resulting in inverted genomic regions being favored over collinear regions. This process allows inversions to establish to high frequencies. Our model can account for several classic patterns in the geographic distribution of inversions and highlights how selection on standing genetic variation allows rapid chromosomal evolution without the waiting time for new mutations. As inversion differences often separate closely related taxa, mixed modes of divergence could be common.     

Genomic divergence across ecological gradients in the Central African rainforest songbird (Andropadus virens)

The little greenbul, a common rainforest passerine from sub‐Saharan Africa, has been the subject of long‐term evolutionary studies to understand the mechanisms leading to rainforest speciation. Previous research found morphological and behavioural divergence across rainforest–savannah transition zones (ecotones), and a pattern of divergence with gene flow suggesting divergent natural selection has contributed to adaptive divergence and ecotones could be important areas for rainforests speciation. Recent advances in genomics and environmental modelling make it possible to examine patterns of genetic divergence in a more comprehensive fashion. To assess the extent to which natural selection may drive patterns of differentiation, here we investigate patterns of genomic differentiation among populations across environmental gradients and regions. We find compelling evidence that individuals form discrete genetic clusters corresponding to distinctive environmental characteristics and habitat types. Pairwise F_ST between populations in different habitats is significantly higher than within habitats, and this differentiation is greater than what is expected from geographic distance alone. Moreover, we identified 140 SNPs that showed extreme differentiation among populations through a genomewide selection scan. These outliers were significantly enriched in exonic and coding regions, suggesting their functional importance. Environmental association analysis of SNP variation indicates that several environmental variables, including temperature and elevation, play important roles in driving the pattern of genomic diversification. Results lend important new genomic evidence for environmental gradients being important in population differentiation.     

Genome divergence during evolutionary diversification as revealed in replicate lake-stream stickleback population pairs

Evolutionary diversification is often initiated by adaptive divergence between populations occupying ecologically distinct environments while still exchanging genes. The genetic foundations of this divergence process are largely unknown and are here explored through genome scans in multiple independent lake-stream population pairs of threespine stickleback. We find that across the pairs, overall genomic divergence is associated with the magnitude of divergence in phenotypes known to be under divergent selection. Along this same axis of increasing diversification, genomic divergence becomes increasingly biased towards the centre of chromosomes as opposed to the peripheries. We explain this pattern by within-chromosome variation in the physical extent of hitchhiking, as recombination is greatly reduced in chromosome centres. Correcting for this effect suggests that a great number of genes distributed widely across the genome are involved in the divergence into lake vs. stream habitats. Analyzing additional allopatric population pairs, however, reveals that strong divergence in some genomic regions has been driven by selection unrelated to lake-stream ecology. Our study highlights a major contribution of large-scale variation in recombination rate to generating heterogeneous genomic divergence and indicates that elucidating the genetic basis of adaptive divergence might be more challenging than currently recognized.     

THE EFFICACY OF DIVERGENCE HITCHHIKING IN GENERATING GENOMIC ISLANDS DURING ECOLOGICAL SPECIATION

Genes under divergent selection flow less readily between populations than other loci. This observation has led to verbal “divergence hitchhiking” models of speciation in which decreased interpopulation gene flow surrounding loci under divergent selection can generate large regions of differentiation within the genome (genomic islands). The efficacy of this model in promoting speciation depends on the size of the region affected by divergence hitchhiking. Empirical evidence is mixed, with examples of both large and small genomic islands. To address these empirical discrepancies and to formalize the theory, we present mathematical models of divergence hitchhiking, which examine neutral differentiation around selected sites. For a single locus under selection, regions of differentiation do not extend far along a chromosome away from a selected site unless both effective population sizes and migration rates are low. When multiple loci are considered, regions of differentiation can be larger. However, with many loci under selection, genome‐wide divergence occurs and genomic islands are erased. The results show that divergence hitchhiking can generate large regions of differentiation, but that the conditions under which this occurs are limited. Thus, speciation may often require multifarious selection acting on many, isolated and physically unlinked genes. How hitchhiking promotes further adaptive divergence warrants consideration.     

CHROMOSOMAL INVERSIONS AND SPECIES DIFFERENCES: WHEN ARE GENES AFFECTING ADAPTIVE DIVERGENCE AND REPRODUCTIVE ISOLATION EXPECTED TO RESIDE WITHIN INVERSIONS?

Many factors can promote speciation, and one which has received much attention is chromosomal inversions. A number of models propose that the recombination suppressing effects of inversions facilitate the maintenance of differences between interbreeding populations in genes affecting adaptive divergence and reproductive isolation. These models predict that such genes will disproportionately reside within inversions, rather than in collinear regions. This hypothesis has received some support, but exceptions exist. Additionally, the effects of known low levels of recombination within inversions on these models are uninvestigated. Here, simulations are used to compare the maintenance of genetic differences between populations following secondary contact and hybridization in different inversion models. We compare regions with no recombination within them to regions with low recombination and to collinear regions with free recombination. Our most general finding is that the low levels of recombination within an inversion often result in the loss of accentuated divergence in inverted regions compared to collinear ones. We conclude that inversions can facilitate the maintenance of species differences under some conditions, but that large or qualitative differences between inverted and collinear regions need not occur. We also find that strong selection facilitates maintenance of divergence in a manner analogous to inversions.     

Extensive linkage disequilibrium and parallel adaptive divergence across threespine stickleback genomes

Population genomic studies are beginning to provide a more comprehensive view of dynamic genome-scale processes in evolution. Patterns of genomic architecture, such as genomic islands of increased divergence, may be important for adaptive population differentiation and speciation. We used next-generation sequencing data to examine the patterns of local and long-distance linkage disequilibrium (LD) across oceanic and freshwater populations of threespine stickleback, a useful model for studies of evolution and speciation. We looked for associations between LD and signatures of divergent selection, and assessed the role of recombination rate variation in generating LD patterns. As predicted under the traditional biogeographic model of unidirectional gene flow from ancestral oceanic to derived freshwater stickleback populations, we found extensive local and long-distance LD in fresh water. Surprisingly, oceanic populations showed similar patterns of elevated LD, notably between large genomic regions previously implicated in adaptation to fresh water. These results support an alternative biogeographic model for the stickleback radiation, one of a metapopulation with appreciable bi-directional gene flow combined with strong divergent selection between oceanic and freshwater populations. As predicted by theory, these processes can maintain LD within and among genomic islands of divergence. These findings suggest that the genomic architecture in oceanic stickleback populations may provide a mechanism for the rapid re-assembly and evolution of multi-locus genotypes in newly colonized freshwater habitats, and may help explain genetic mapping of parallel phenotypic variation to similar loci across independent freshwater populations.     

The genetic mosaic suggests a new role for hitchhiking in ecological speciation

Early in ecological speciation, the genomically localized effects of divergent selection cause heterogeneity among loci in divergence between incipient species. We call this pattern of genomic variability in divergence the 'genetic mosaic of speciation'. Previous studies have used F(ST) outliers as a way to identify divergently selected genomic regions, but the nature of the relationship between outlier loci and quantitative trait loci (QTL) involved in reproductive isolation has not yet been quantified. Here, we show that F(ST) outliers between a pair of incipient species are significantly clustered around QTL for traits that cause ecologically based reproductive isolation. Around these key QTL, extensive 'divergence hitchhiking' occurs because reduced inter-race mating and negative selection decrease the opportunity for recombination between chromosomes bearing different locally adapted QTL alleles. Divergence hitchhiking is likely to greatly increase the opportunity for speciation in populations that are sympatric, regardless of whether initial divergence was sympatric or allopatric. Early in ecological speciation, analyses of population structure, gene flow or phylogeography based on different random or arbitrarily chosen neutral markers should be expected to conflict--only markers in divergently selected genomic regions will reveal the evolutionary history of adaptive divergence and ecologically based reproductive isolation. Species retain mosaic genomes for a very long time, and gene exchange in hybrid zones can vary dramatically among loci. However, in hybridizing species, the genomic regions that affect ecologically based reproductive isolation are difficult to distinguish from regions that have diverged for other reasons.     

Modifying effects of phenotypic plasticity on interactions among natural selection, adaptation and gene flow

Divergent natural selection, adaptive divergence and gene flow may interact in a number of ways. Recent studies have focused on the balance between selection and gene flow in natural populations, and empirical work has shown that gene flow can constrain adaptive divergence, and that divergent selection can constrain gene flow. A caveat is that phenotypic diversification may be under the direct influence of environmental factors (i.e. it may be due to phenotypic plasticity), in addition to partial genetic influence. In this case, phenotypic divergence may occur between populations despite high gene flow that imposes a constraint on genetic divergence. Plasticity may dampen the effects of natural selection by allowing individuals to rapidly adapt phenotypically to new conditions, thus slowing adaptive genetic divergence. On the other hand, plasticity may promote future adaptive divergence by allowing populations to persist in novel environments. Plasticity may promote gene flow between selective regimes by allowing dispersers to adapt to alternate conditions, or high gene flow may result in the selection for increased plasticity. Here I expand frameworks for understanding relationships among selection, adaptation and gene flow to include the effects of phenotypic plasticity in natural populations, and highlight its importance in evolutionary diversification.     

Genomics of adaptive divergence with chromosome‐scale heterogeneity in crossover rate

Genetic differentiation between divergent populations is often greater in chromosome centres than peripheries. Commonly overlooked, this broadscale differentiation pattern is sometimes ascribed to heterogeneity in crossover rate and hence linked selection within chromosomes, but the underlying mechanisms remain incompletely understood. A literature survey across 46 organisms reveals that most eukaryotes indeed exhibit a reduced crossover rate in chromosome centres relative to the peripheries. Using simulations of populations diverging into ecologically different habitats through sorting of standing genetic variation, we demonstrate that such chromosome‐scale heterogeneity in crossover rate, combined with polygenic divergent selection, causes stronger hitchhiking and especially barriers to gene flow across chromosome centres. Without requiring selection on new mutations, this rapidly leads to elevated population differentiation in the low‐crossover centres relative to the high‐crossover peripheries of chromosomes (“Chromosome Centre‐Biased Differentiation”, CCBD). Using simulated and empirical data, we then show that strong CCBD between populations can provide evidence of polygenic adaptive divergence with a phase of gene flow. We further demonstrate that chromosome‐scale heterogeneity in crossover rate impacts analyses beyond that of population differentiation, including the inference of phylogenies and parallel adaptive evolution among populations, the detection of genetic loci under selection, and the interpretation of the strength of selection on genomic regions. Overall, our results call for a greater appreciation of chromosome‐scale heterogeneity in crossover rate in evolutionary genomics.     

Evolution of putative barrier loci at an intermediate stage of speciation with gene flow in campions (Silene)

Understanding the origin of new species is a central goal in evolutionary biology. Diverging lineages often evolve highly heterogeneous patterns of genetic differentiation; however, the underlying mechanisms are not well understood. We investigated evolutionary processes governing genetic differentiation between the hybridizing campions Silene dioica (L.) Clairv. and S. latifolia Poiret. Demographic modelling indicated that the two species diverged with gene flow. The best‐supported scenario with heterogeneity in both migration rate and effective population size suggested that a small proportion of the loci evolved without gene flow. Differentiation (F _ST) and sequence divergence (d _XY) were correlated and both tended to peak in the middle of most linkage groups, consistent with reduced gene flow at highly differentiated loci. Highly differentiated loci further exhibited signatures of selection. In between‐species population pairs, isolation by distance was stronger for genomic regions with low between‐species differentiation than for highly differentiated regions that may contain barrier loci. Moreover, differentiation landscapes within and between species were only weakly correlated, suggesting that linked selection due to shared recombination and gene density landscapes is not the dominant determinant of genetic differentiation in these lineages. Instead, our results suggest that divergent selection shaped the genomic landscape of differentiation between the two Silene species, consistent with predictions for 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.