The distinctive footprints of local hitchhiking in a varied environment and global hitchhiking in a subdivided population

Detecting the footprints of divergent selection in oaks with linked markers

Genome scans are increasingly used to study ecological speciation, providing a useful genome-wide perspective on divergent selection in the presence of gene flow. Here, we compare current approaches to detect footprints of divergent selection in closely related species. We analyzed 192 individuals from two interfertile European temperate oak species using 30 nuclear microsatellites from eight linkage groups. These markers present little intraspecific differentiation and can be used in combination to assign individual genotypes to species. We first show that different outlier detection tests give somewhat different results, possibly due to model constraints. Second, using linkage information for these markers, we further characterize the signature of divergent selection in the presence of gene flow. In particular, we show that recombination estimates for regions with outlier markers are lower than those for a control region, in line with a prediction from ecological speciation theory. Most importantly, we show that analyses at the haplotype level can distinguish between truly divergent (bi-directional) selection and positive selection in one of the two species, offering a new and improved method for characterizing the speciation process.

Gene flow at major transitional areas in sea bass (Dicentrarchus labrax) and the possible emergence of a hybrid swarm

The population genetic structure of sea bass (Dicentrarchus labrax) along a transect from the Atlantic Ocean (AO) to the Eastern Mediterranean (EM) Sea differs from that of most other marine taxa in this area. Three populations (AO, Western Mediterranean [WM], EM) are recognized today, which were originally two allopatric populations. How two ancestral genetic units have evolved into three distinct units has not been addressed yet. Therefore, to investigate mechanisms that lead to the emergence of the central WM population, its current status, and its connectivity with the two parental populations, we applied 20 nuclear loci that were either gene associated or gene independent. Results confirmed the existence of three distinct gene pools, with higher differentiation at two transitional areas, the Almeria-Oran Front (AOF) and of the Siculo-Tunisian Strait (STS), than within any population. Significant linkage disequilibrium and heterozygote excess indicated that the STS is probably another tension zone, as already described for the AOF. Neutrality tests fail to reveal marker loci that could be driven by selection within or among metapopulations, except for locus DLA0068. Collectively, results support that the central WM population arose by trapping two tensions zones at distinct geographic locations of limited connectivity. Population assignment further revealed that WM individuals were more introgressed than individuals from the other two metapopulations. This suggests that this population might result from hybrid swarming, and was or is still seeded by genes received through the filter of each tension zone.

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.

A framework for estimating the fixation time of an advantageous allele in stepping-stone models

Determining how population subdivision increases the fixation time of an advantageous allele is an important problem in evolutionary genetics as this influences many processes. Here, I lay out a framework for calculating the fixation time of a positively selected allele in a subdivided population, as a function of the number of demes present, the migration rate between them and the manner in which they are connected. Using this framework, it becomes clear that a beneficial allele's fixation time is significantly reduced through migration continuously introducing copies of the allele into a newly colonized subpopulation, increasing its frequency within these demes. The effect that migration has on allele frequency needs to be explicitly taken into account to produce a realistic estimate of fixation time. This behaviour is most prominent when demes are arranged on a two-dimensional torus, in comparison with populations where demes are arranged in a circle. This is because each subpopulation is connected to several neighbours over a torus, so that there are multiple paths that an allele can take in order to fix. As a consequence, some demes experience a greater influx and efflux of migrants than others. Analytical results are found to be very accurate when compared to stochastic simulations, and are generally robust if there are a large number of demes, or if the allele is weakly selected for.

AFLP genome scans suggest divergent selection on colour patterning in allopatric colour morphs of a cichlid fish

Genome scan-based tests for selection are directly applicable to natural populations to study the genetic and evolutionary mechanisms behind phenotypic differentiation. We conducted AFLP genome scans in three distinct geographic colour morphs of the cichlid fish Tropheus moorii to assess whether the extant, allopatric colour pattern differentiation can be explained by drift and to identify markers mapping to genomic regions possibly involved in colour patterning. The tested morphs occupy adjacent shore sections in southern Lake Tanganyika and are separated from each other by major habitat barriers. The genome scans revealed significant genetic structure between morphs, but a very low proportion of loci fixed for alternative AFLP alleles in different morphs. This high level of polymorphism within morphs suggested that colour pattern differentiation did not result exclusively from neutral processes. Outlier detection methods identified six loci with excess differentiation in the comparison between a bluish and a yellow-blotch morph and five different outlier loci in comparisons of each of these morphs with a red morph. As population expansions and the genetic structure of Tropheus make the outlier approach prone to false-positive signals of selection, we examined the correlation between outlier locus alleles and colour phenotypes in a genetic and phenotypic cline between two morphs. Distributions of allele frequencies at one outlier locus were indeed consistent with linkage to a colour locus. Despite the challenges posed by population structure and demography, our results encourage the cautious application of genome scans to studies of divergent selection in subdivided and recently expanded populations.

Genome-wide patterns of divergence during speciation: the lake whitefish case study

The nature, size and distribution of the genomic regions underlying divergence and promoting reproductive isolation remain largely unknown. Here, we summarize ongoing efforts using young (12 000 yr BP) species pairs of lake whitefish (Coregonus clupeaformis) to expand our understanding of the initial genomic patterns of divergence observed during speciation. Our results confirmed the predictions that: (i) on average, phenotypic quantitative trait loci (pQTL) show higher F(ST) values and are more likely to be outliers (and therefore candidates for being targets of divergent selection) than non-pQTL markers; (ii) large islands of divergence rather than small independent regions under selection characterize the early stages of adaptive divergence of lake whitefish; and (iii) there is a general trend towards an increase in terms of numbers and size of genomic regions of divergence from the least (East L.) to the most differentiated species pair (Cliff L.). This is consistent with previous estimates of reproductive isolation between these species pairs being driven by the same selective forces responsible for environment specialization. Altogether, dwarf and normal whitefish species pairs represent a continuum of both morphological and genomic differentiation contributing to ecological speciation. Admittedly, much progress is still required to more finely map and circumscribe genomic islands of speciation. This will be achieved through the use of next generation sequencing data but also through a better quantification of phenotypic traits moulded by selection as organisms adapt to new environmental conditions.

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.

Dynamics of drift, gene flow, and selection during speciation in Silene

The mechanics of speciation with gene flow are still unclear. Disparity among genes in population differentiation (F(ST)) between diverging species is often interpreted as evidence for semipermeable species boundaries, with selection preventing "key" genes from introgressing despite ongoing gene flow. However, F(ST) can remain high before it reaches equilibrium between the lineage sorting of species divergence and the homogenizing effects of gene flow (via secondary contact). Thus, when interpreting F(ST), the dynamics of drift, gene flow, and selection need to be taken into account. We illustrate this view with a multigenic analyses of gene flow and selection in three closely related Silene species, S. latifolia, S. dioica, and S. diclinis. We report that although S. diclinis appears to have evolved in allopatry, isolation with (bidirectional) gene flow between S. latifolia and S. dioica is likely, perhaps as a result of parapatric speciation followed by more extensive sympatry. Interestingly, we detected the signatures of apparently independent instances of positive selection at the same locus in S. latifolia and S. dioica. Despite gene flow between the species, the adaptive alleles have not crossed the species boundary, suggesting that this gene has independently undergone species-specific (diversifying or parallel) selection.

Evidence for hitchhiking of deleterious mutations within the human genome

Deleterious mutations present a significant obstacle to adaptive evolution. Deleterious mutations can inhibit the spread of linked adaptive mutations through a population; conversely, adaptive substitutions can increase the frequency of linked deleterious mutations and even result in their fixation. To assess the impact of adaptive mutations on linked deleterious mutations, we examined the distribution of deleterious and neutral amino acid polymorphism in the human genome. Within genomic regions that show evidence of recent hitchhiking, we find fewer neutral but a similar number of deleterious SNPs compared to other genomic regions. The higher ratio of deleterious to neutral SNPs is consistent with simulated hitchhiking events and implies that positive selection eliminates some deleterious alleles and increases the frequency of others. The distribution of disease-associated alleles is also altered in hitchhiking regions. Disease alleles within hitchhiking regions have been associated with auto-immune disorders, metabolic diseases, cancers, and mental disorders. Our results suggest that positive selection has had a significant impact on deleterious polymorphism and may be partly responsible for the high frequency of certain human disease alleles.

Hitchhiking effect of a beneficial mutation spreading in a subdivided population

A central problem in population genetics is to detect and analyze positive natural selection by which beneficial mutations are driven to fixation. The hitchhiking effect of a rapidly spreading beneficial mutation, which results in local removal of standing genetic variation, allows such an analysis using DNA sequence polymorphism. However, the current mathematical theory that predicts the pattern of genetic hitchhiking relies on the assumption that a beneficial mutation increases to a high frequency in a single random-mating population, which is certainly violated in reality. Individuals in natural populations are distributed over a geographic space. The spread of a beneficial allele can be delayed by limited migration of individuals over the space and its hitchhiking effect can also be affected. To study this effect of geographic structure on genetic hitchhiking, we analyze a simple model of directional selection in a subdivided population. In contrast to previous studies on hitchhiking in subdivided populations, we mainly investigate the range of sufficiently high migration rates that would homogenize genetic variation at neutral loci. We provide a heuristic mathematical analysis that describes how the genealogical structure at a neutral locus linked to the locus under selection is expected to change in a population divided into two demes. Our results indicate that the overall strength of genetic hitchhiking--the degree to which expected heterozygosity decreases--is diminished by population subdivision, mainly because opportunity for the breakdown of hitchhiking by recombination increases as the spread of the beneficial mutation across demes is delayed when migration rate is much smaller than the strength of selection. Furthermore, the amount of genetic variation after a selective sweep is expected to be unequal over demes: a greater reduction in expected heterozygosity occurs in the subpopulation from which the beneficial mutation originates than in its neighboring subpopulations. This raises a possibility of detecting a "hidden" geographic structure of population by carefully analyzing the pattern of a selective sweep.