Population genomics of secondary contact

Hybrid fitness effects modify fixation probabilities of introgressed alleles

Hybridization is a common occurrence in natural populations, and introgression is a major source of genetic variation. Despite the evolutionary importance of adaptive introgression, classical population genetics theory does not take into account hybrid fitness effects. Specifically, heterosis (i.e. hybrid vigor) and Dobzhansky–Muller incompatibilities influence the fates of introgressed alleles. Here, we explicitly account for polygenic, unlinked hybrid fitness effects when tracking a rare introgressed marker allele. These hybrid fitness effects quickly decay over time due to repeated backcrossing, enabling a separation-of-timescales approach. Using diffusion and branching process theory in combination with computer simulations, we formalize the intuition behind how hybrid fitness effects affect introgressed alleles. We find that hybrid fitness effects can significantly hinder or boost the fixation probability of introgressed alleles, depending on the relative strength of heterosis and Dobzhansky–Muller incompatibilities effects. We show that the inclusion of a correction factor (α, representing the compounded effects of hybrid fitness effects over time) into classic population genetics theory yields accurate fixation probabilities. Despite having a strong impact on the probability of fixation, hybrid fitness effects only subtly change the distribution of fitness effects of introgressed alleles that reach fixation. Although strong Dobzhansky–Muller incompatibility effects may expedite the loss of introgressed alleles, fixation times are largely unchanged by hybrid fitness effects.

Gene flow as a simple cause for an excess of high-frequency-derived alleles

Most human populations exhibit an excess of high-frequency variants, leading to a U-shaped site-frequency spectrum (uSFS). This pattern has been generally interpreted as a signature of ongoing episodes of positive selection, or as evidence for a mis-assignment of ancestral/derived allelic states, but uSFS has also been observed in populations receiving gene flow from a ghost population, in structured populations, or after range expansions. In order to better explain the prevalence of high-frequency variants in humans and other populations, we describe here which patterns of gene flow and population demography can lead to uSFS by using extensive coalescent simulations. We find that uSFS can often be observed in a population if gene flow brings a few ancestral alleles from a well-differentiated population. Gene flow can either consist in single pulses of admixture or continuous immigration, but different demographic conditions are necessary to observe uSFS in these two scenarios. Indeed, an extremely low and recent gene flow is required in the case of single admixture events, while with continuous immigration, uSFS occurs only if gene flow started recently at a high rate or if it lasted for a long time at a low rate. Overall, we find that a neutral uSFS occurs under more restrictive conditions in populations having received single pulses of gene flow than in populations exposed to continuous gene flow. We also show that the uSFS observed in human populations from the 1000 Genomes Project can easily be explained by gene flow from surrounding populations without requiring past episodes of positive selection. These results imply that uSFS should be common in non-isolated populations, such as most wild or domesticated plants and animals.

The genomic landscape at a late stage of stickleback speciation: High genomic divergence interspersed by small localized regions of introgression

Speciation is a continuous process and analysis of species pairs at different stages of divergence provides insight into how it unfolds. Previous genomic studies on young species pairs have revealed peaks of divergence and heterogeneous genomic differentiation. Yet less known is how localised peaks of differentiation progress to genome-wide divergence during the later stages of speciation in the presence of persistent gene flow. Spanning the speciation continuum, stickleback species pairs are ideal for investigating how genomic divergence builds up during speciation. However, attention has largely focused on young postglacial species pairs, with little knowledge of the genomic signatures of divergence and introgression in older stickleback systems. The Japanese stickleback species pair, composed of the Pacific Ocean three-spined stickleback (Gasterosteus aculeatus) and the Japan Sea stickleback (G. nipponicus), which co-occur in the Japanese islands, is at a late stage of speciation. Divergence likely started well before the end of the last glacial period and crosses between Japan Sea females and Pacific Ocean males result in hybrid male sterility. Here we use coalescent analyses and Approximate Bayesian Computation to show that the two species split approximately 0.68-1 million years ago but that they have continued to exchange genes at a low rate throughout divergence. Population genomic data revealed that, despite gene flow, a high level of genomic differentiation is maintained across the majority of the genome. However, we identified multiple, small regions of introgression, occurring mainly in areas of low recombination rate. Our results demonstrate that a high level of genome-wide divergence can establish in the face of persistent introgression and that gene flow can be localized to small genomic regions at the later stages of speciation with gene flow.

A new method to scan genomes for introgression in a secondary contact model

Secondary contact between divergent populations or incipient species may result in the exchange and introgression of genomic material. We develop a simple DNA sequence measure, called G_min, which is designed to identify genomic regions experiencing introgression in a secondary contact model. G_min is defined as the ratio of the minimum between-population number of nucleotide differences in a genomic window to the average number of between-population differences. Although it is conceptually simple, one advantage of G_min is that it is computationally inexpensive relative to model-based methods for detecting gene flow and it scales easily to the level of whole-genome analysis. We compare the sensitivity and specificity of G_min to those of the widely used index of population differentiation, F_ST, and suggest a simple statistical test for identifying genomic outliers. Extensive computer simulations demonstrate that G_min has both greater sensitivity and specificity for detecting recent introgression than does F_ST. Furthermore, we find that the sensitivity of G_min is robust with respect to both the population mutation and recombination rates. Finally, a scan of G_min across the X chromosome of Drosophila melanogaster identifies candidate regions of introgression between sub-Saharan African and cosmopolitan populations that were previously missed by other methods. These results show that G_min is a biologically straightforward, yet powerful, alternative to F_ST, as well as to more computationally intensive model-based methods for detecting gene flow.

Hybrid fitness, adaptation and evolutionary diversification: lessons learned from Louisiana Irises

Estimates of hybrid fitness have been used as either a platform for testing the potential role of natural hybridization in the evolution of species and species complexes or, alternatively, as a rationale for dismissing hybridization events as being of any evolutionary significance. From the time of Darwin's publication of The Origin, through the neo-Darwinian synthesis, to the present day, the observation of variability in hybrid fitness has remained a challenge for some models of speciation. Yet, Darwin and others have reported the elevated fitness of hybrid genotypes under certain environmental conditions. In modern scientific terminology, this observation reflects the fact that hybrid genotypes can demonstrate genotype × environment interactions. In the current review, we illustrate the development of one plant species complex, namely the Louisiana Irises, into a 'model system' for investigating hybrid fitness and the role of genetic exchange in adaptive evolution and diversification. In particular, we will argue that a multitude of approaches, involving both experimental and natural environments, and incorporating both manipulative analyses and surveys of natural populations, are necessary to adequately test for the evolutionary significance of introgressive hybridization. An appreciation of the variability of hybrid fitness leads to the conclusion that certain genetic signatures reflect adaptive evolution. Furthermore, tests of the frequency of allopatric versus sympatric/parapatric divergence (that is, divergence with ongoing gene flow) support hybrid genotypes as a mechanism of evolutionary diversification in numerous species complexes.