GENETIC HITCHHIKING AND THE DYNAMIC BUILDUP OF GENOMIC DIVERGENCE DURING SPECIATION WITH GENE FLOW

A major issue in evolutionary biology is explaining patterns of differentiation observed in population genomic data, as divergence can be due to both direct selection on a locus and genetic hitchhiking. “Divergence hitchhiking” (DH) theory postulates that divergent selection on a locus reduces gene flow at physically linked sites, facilitating the formation of localized clusters of tightly linked, diverged loci. “Genome hitchhiking” (GH) theory emphasizes genome‐wide effects of divergent selection. Past theoretical investigations of DH and GH focused on static snapshots of divergence. Here, we used simulations assessing a variety of strengths of selection, migration rates, population sizes, and mutation rates to investigate the relative importance of direct selection, GH, and DH in facilitating the dynamic buildup of genomic divergence as speciation proceeds through time. When divergently selected mutations were limiting, GH promoted divergence, but DH had little measurable effect. When populations were small and divergently selected mutations were common, DH enhanced the accumulation of weakly selected mutations, but this contributed little to reproductive isolation. In general, GH promoted reproductive isolation by reducing effective migration rates below that due to direct selection alone, and was important for genome‐wide “congealing” or “coupling” of differentiation (F_ST) across loci as speciation progressed.     

Spatially explicit models of divergence and genome hitchhiking

Strong barriers to genetic exchange can exist at divergently selected loci, whereas alleles at neutral loci flow more readily between populations, thus impeding divergence and speciation in the face of gene flow. However, ‘divergence hitchhiking’ theory posits that divergent selection can generate large regions of differentiation around selected loci. ‘Genome hitchhiking’ theory suggests that selection can also cause reductions in average genome‐wide rates of gene flow, resulting in widespread genomic divergence (rather than divergence only around specific selected loci). Spatial heterogeneity is ubiquitous in nature, yet previous models of genetic barriers to gene flow have explored limited combinations of spatial and selective scenarios. Using simulations of secondary contact of populations, we explore barriers to gene flow in various selective and spatial contexts in continuous, two‐dimensional, spatially explicit environments. In general, the effects of hitchhiking are strongest in environments with regular spatial patterning of starkly divergent habitat types. When divergent selection is very strong, the absence of intermediate habitat types increases the effects of hitchhiking. However, when selection is moderate or weak, regular (vs. random) spatial arrangement of habitat types becomes more important than the presence of intermediate habitats per se. We also document counterintuitive processes arising from the stochastic interplay between selection, gene flow and drift. Our results indicate that generalization of results from two‐deme models requires caution and increase understanding of the genomic and geographic basis of population divergence.     

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.     

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.     

Isolation by adaptation in Neochlamisus leaf beetles: host-related selection promotes neutral genomic divergence

This study tests how divergent natural selection promotes genomic differentiation during ecological speciation. Specifically, we use adaptive ecological divergence (here, population divergence in host plant use and preference) as a proxy for selection strength and evaluate the correlation between levels of adaptive and genetic differentiation across pairwise population comparisons. Positive correlations would reveal the pattern predicted by our hypothesis, that of 'isolation by adaptation' (IBA). Notably, IBA is predicted not only for selected loci but also for neutral loci. This may reflect the effects of divergent selection on neutral loci that are 'loosely linked' to divergently selected loci or on geneflow restriction that facilitates genetic drift at all loci, including neutral loci that are completely unlinked to those evolving under divergent selection. Here, we evaluate IBA in maple- and willow-associated populations of Neochlamisus bebbianae leaf beetles. To do so, we collected host preference data to construct adaptive divergence indices and used AFLPs (amplified fragment length polymorphisms) and mitochondrial sequences to quantify genetic differentiation. Partial Mantel tests showed significant IBA in 'pooled' analyses of putatively neutral and of putatively selected ('outlier') AFLP loci. This pattern was also recovered in 12% of 'locus-specific' analyses that separately evaluated genetic differentiation at individual neutral loci. These results provided evidence for widespread effects of selection on neutral genomic divergence. Our collective findings indicate that host-related selection may play important roles in the population genomic differentiation of both neutral and selected gene regions in herbivorous insects.     

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.     

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.     

Migration-induced phenotypic divergence: the migration-selection balance of correlated traits

Genetically correlated traits are known to respond to indirect selection pressures caused by directional selection on other traits. It is however unclear how local adaptation in populations diverging along some phenotypic traits but not others is affected by the joint action of gene flow and genetic correlations among traits. This simulation study shows that although gene flow is a potent constraining mechanism of population adaptive divergence, it may induce phenotypic divergence in traits under homogeneous selection among habitats if they are genetically correlated with traits under divergent selection. This correlated phenotypic divergence is a nonmonotonous function of migration and increases with mutational correlation among traits. It also increases with the number of divergently selected traits provided their genetic autonomy relative to the uniformly selected trait is reduced by specific patterns of genetic covariances: populations with lower effective trait dimensionality are more likely to generate very large correlated divergence. The correlated divergence is likely to be picked up by Q(ST)-F(ST) analysis of population genetic differentiation and be erroneously ascribed to adaptive divergence under divergent selection. This study emphasizes the necessity to understand the interaction between selection and the genetic basis of adaptation in a multivariate rather than univariate context.     

When can ecological speciation be detected with neutral loci?

It is not yet clear under what conditions empirical studies can reliably detect progress toward ecological speciation through the analysis of allelic variation at neutral loci. We use a simulation approach to investigate the range of parameter space under which such detection is, and is not, likely. We specifically test for the conditions under which divergent natural selection can cause a ‘generalized barrier to gene flow’ that is present across the genome. Our individual‐based numerical simulations focus on how population divergence at neutral loci varies in relation to recombination rate with a selected locus, divergent selection on that locus, migration rate and population size. We specifically test whether genetic differences at neutral markers are greater between populations in different environments than between populations in similar environments. We find that this expected signature of ecological speciation can be detected under part of the parameter space, most consistently when divergent selection is strong and migration is intermediate. By contrast, the expected signature of ecological speciation is not reliably detected when divergent selection is weak or migration is low or high. These findings provide insights into the strengths and weaknesses of using neutral markers to infer ecological speciation in natural systems.     

The genomic signature of parallel adaptation from shared genetic variation

Parallel adaptation is common and may often occur from shared genetic variation, but the genomic consequences of this process remain poorly understood. We first use individual‐based simulations to demonstrate that comparisons between populations adapted in parallel to similar environments from shared variation reveal a characteristic genomic signature around a selected locus: a low‐divergence valley centred at the locus and flanked by twin peaks of high divergence. This signature is initiated by the hitchhiking of haplotype tracts differing between derived populations in the broader neighbourhood of the selected locus (driving the high‐divergence twin peaks) and shared haplotype tracts in the tight neighbourhood of the locus (driving the low‐divergence valley). This initial hitchhiking signature is reinforced over time because the selected locus acts as a barrier to gene flow from the source to the derived populations, thus promoting divergence by drift in its close neighbourhood. We next empirically confirm the peak‐valley‐peak signature by combining targeted and RAD sequence data at three candidate adaptation genes in multiple marine (source) and freshwater (derived) populations of threespine stickleback. Finally, we use a genome‐wide screen for the peak‐valley‐peak signature to discover additional genome regions involved in parallel marine‐freshwater divergence. Our findings offer a new explanation for heterogeneous genomic divergence and thus challenge the standard view that peaks in population divergence harbour divergently selected loci and that low‐divergence regions result from balancing selection or localized introgression. We anticipate that genome scans for peak‐valley‐peak divergence signatures will promote the discovery of adaptation genes in other organisms.     

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.     

Consequences of divergence and introgression for speciation in Andean cloud forest birds

Divergence with gene flow is well documented and reveals the influence of ecological adaptation on speciation. Yet, it remains intuitive that gene exchange inhibits speciation in many scenarios, particularly among ecologically similar populations. The influence of gene flow on the divergence of populations facing similar selection pressures has received less empirical attention than scenarios where differentiation is coupled with local environmental adaptation. I used a paired study design to test the influence of genomic divergence and introgression on plumage differentiation between ecologically similar allopatric replacements of Andean cloud forest birds. Through analyses of short‐read genome‐wide sequences from over 160 individuals in 16 codistributed lineages, I found that plumage divergence is associated with deep genetic divergence, implicating a prominent role of geographic isolation in speciation. By contrast, lineages that lack plumage divergence across the same geographic barrier are more recently isolated or exhibit a signature of secondary genetic introgression, indicating a negative relationship between gene flow and divergence in phenotypic traits important to speciation. My results suggest that the evolutionary outcomes of cycles of isolation and divergence in this important theatre of biotic diversification are sensitive to time spent in the absence of gene flow.     

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.     

Interpreting the genomic landscape of speciation: a road map for finding barriers to gene flow

Speciation, the evolution of reproductive isolation among populations, is continuous, complex, and involves multiple, interacting barriers. Until it is complete, the effects of this process vary along the genome and can lead to a heterogeneous genomic landscape with peaks and troughs of differentiation and divergence. When gene flow occurs during speciation, barriers restricting gene flow locally in the genome lead to patterns of heterogeneity. However, genomic heterogeneity can also be produced or modified by variation in factors such as background selection and selective sweeps, recombination and mutation rate variation, and heterogeneous gene density. Extracting the effects of gene flow, divergent selection and reproductive isolation from such modifying factors presents a major challenge to speciation genomics. We argue one of the principal aims of the field is to identify the barrier loci involved in limiting gene flow. We first summarize the expected signatures of selection at barrier loci, at the genomic regions linked to them and across the entire genome. We then discuss the modifying factors that complicate the interpretation of the observed genomic landscape. Finally, we end with a road map for future speciation research: a proposal for how to account for these modifying factors and to progress towards understanding the nature of barrier loci. Despite the difficulties of interpreting empirical data, we argue that the availability of promising technical and analytical methods will shed further light on the important roles that gene flow and divergent selection have in shaping the genomic landscape of speciation.     

Selection and genomic differentiation during ecological speciation: isolating the contributions of host association via a comparative genome scan of Neochlamisus bebbianae leaf beetles

This study uses a comparative genome scan to evaluate the contributions of host plant related divergent selection to genetic differentiation and ecological speciation in maple- and willow-associated populations of Neochlamisus bebbianae leaf beetles. For each of 15 pairwise population comparisons, we identified "outlier loci" whose strong differentiation putatively reflects divergent selection. Of 447 AFLP loci, 15% were outliers across multiple population comparisons, and low linkage disequilibrium indicated that these outliers derived from multiple regions of the genome. Outliers were further classified as "host-specific" if repeatedly observed in "different-host" population comparisons but never in "same-host" comparisons. Outliers exhibiting the opposite pattern were analogously classified as "host-independent." Host-specific outliers represented 5% of all loci and were more frequent than host-independent outliers, thus revealing a large role for host-adaptation in population genomic differentiation. Evidence that host-related selection can promote divergence despite gene flow was provided by population trees. These were structured by host-association when datasets included host-specific outliers, but not when based on neutral loci, which united sympatric populations. Lastly, three host-specific outliers were highly differentiated in all nine different-host comparisons. Because host-adaptation promotes reproductive isolation in these beetles, these loci provide promising candidate gene regions for future molecular studies of ecological speciation.     

What can patterns of differentiation across plant genomes tell us about adaptation and speciation?

Genome scans have become a common approach to identify genomic signatures of natural selection and reproductive isolation, as well as the genomic bases of ecologically relevant phenotypes, based on patterns of polymorphism and differentiation among populations or species. Here, we review the results of studies taking genome scan approaches in plants, consider the patterns of genomic differentiation documented and their possible causes, discuss the results in light of recent models of genomic differentiation during divergent adaptation and speciation, and consider assumptions and caveats in their interpretation. We find that genomic regions of high divergence generally appear quite small in comparisons of both closely and more distantly related populations, and for the most part, these differentiated regions are spread throughout the genome rather than strongly clustered. Thus, the genome scan approach appears well-suited for identifying genomic regions or even candidate genes that underlie adaptive divergence and/or reproductive barriers. We consider other methodologies that may be used in conjunction with genome scan approaches, and suggest further developments that would be valuable. These include broader use of sequence-based markers of known genomic location, greater attention to sampling strategies to make use of parallel environmental or phenotypic transitions, more integration with approaches such as quantitative trait loci mapping and measures of gene flow across the genome, and additional theoretical and simulation work on processes related to divergent adaptation and speciation.     

Genomics of Divergence along a Continuum of Parapatric Population Differentiation

The patterns of genomic divergence during ecological speciation are shaped by a combination of evolutionary forces. Processes such as genetic drift, local reduction of gene flow around genes causing reproductive isolation, hitchhiking around selected variants, variation in recombination and mutation rates are all factors that can contribute to the heterogeneity of genomic divergence. On the basis of 60 fully sequenced three-spined stickleback genomes, we explore these different mechanisms explaining the heterogeneity of genomic divergence across five parapatric lake and river population pairs varying in their degree of genetic differentiation. We find that divergent regions of the genome are mostly specific for each population pair, while their size and abundance are not correlated with the extent of genome-wide population differentiation. In each pair-wise comparison, an analysis of allele frequency spectra reveals that 25–55% of the divergent regions are consistent with a local restriction of gene flow. Another large proportion of divergent regions (38–75%) appears to be mainly shaped by hitchhiking effects around positively selected variants. We provide empirical evidence that alternative mechanisms determining the evolution of genomic patterns of divergence are not mutually exclusive, but rather act in concert to shape the genome during population differentiation, a first necessary step towards ecological speciation.     

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.     

Convergent evolution and divergent selection: lizards at the White Sands ecotone

Ecological transition zones, where organismal phenotypes result from a delicate balance between selection and migration, highlight the interplay of local adaptation and gene flow. Here, I study the response of an entire species assemblage to natural selection across a common ecotone. Three lizard species, distributed along a dramatic environmental gradient in substrate color, display convergent adaptation of blanched coloration on the gypsum dunes of White Sands National Monument. I investigate the role of gene flow in modulating phenotypic response to selection by quantifying color variation and genetic variation across the ecotone. I find species differences in degree of background matching and in genetic connectivity of populations across the ecotone. Differences among species in phenotypic response to selection scale precisely to levels of genetic isolation. Species with higher levels of gene flow across the ecotone exhibit less dramatic responses to selection. Results also reveal a strong signal of ecologically mediated divergence for White Sands lizards. For all species, phenotypic variation is better explained by habitat similarity than genetic similarity. Convergent evolution of blanched coloration at White Sands clearly reflects the action of strong divergent selection; however, adaptive response appears to be modulated by gene flow and demographic history and can be predicted by divergence-with-gene-flow models.     

Genome scans reveal candidate regions involved in the adaptation to host plant in the pea aphid complex

A major goal in evolutionary biology is to uncover the genetic basis of adaptation. Divergent selection exerted on ecological traits may result in adaptive population differentiation and reproductive isolation and affect differentially the level of genetic divergence along the genome. Genome‐wide scan of large sets of individuals from multiple populations is a powerful approach to identify loci or genomic regions under ecologically divergent selection. Here, we focused on the pea aphid, a species complex of divergent host races, to explore the organization of the genomic divergence associated with host plant adaptation and ecological speciation. We analysed 390 microsatellite markers located at variable distances from predicted genes in replicate samples of sympatric populations of the pea aphid collected on alfalfa, red clover and pea, which correspond to three common host‐adapted races reported in this species complex. Using a method that accounts for the hierarchical structure of our data set, we found a set of 11 outlier loci that show higher genetic differentiation between host races than expected under the null hypothesis of neutral evolution. Two of the outliers are close to olfactory receptor genes and three other nearby genes encoding salivary proteins. The remaining outliers are located in regions with genes of unknown functions, or which functions are unlikely to be involved in interactions with the host plant. This study reveals genetic signatures of divergent selection across the genome and provides an inventory of candidate genes responsible for plant specialization in the pea aphid, thereby setting the stage for future functional studies.