The many-demes limit for selection and drift in a subdivided population

A diffusion approximation is obtained for the frequency of a selected allele in a population comprised of many subpopulations or demes. The form of the diffusion is equivalent to that for an unstructured population, except that it occurs on a longer time scale when migration among demes is restricted. This many-demes diffusion limit relies on the collection of demes always being in statistical equilibrium with respect to migration and drift for a given allele frequency in the total population. Selection is assumed to be weak, in inverse proportion to the number of demes, and the results hold for any deme sizes and migration rates greater than zero. The distribution of allele frequencies among demes is also described.     

Fixation probability for a beneficial allele and a mutant strategy in a linear game under weak selection in a finite island model

The effect of population structure on the probability of fixation of a newly introduced mutant under weak selection is studied using a coalescent approach. Wright's island model in a framework of a finite number of demes is assumed and two selection regimes are considered: a beneficial allele model and a linear game among offspring. A first-order approximation of the fixation probability for a single mutant with respect to the intensity of selection is deduced. The approximation requires the calculation of expected coalescence times, under neutrality, for lineages starting from two or three sampled individuals. The results are obtained in a general setting without assumptions on the number of demes, the deme size or the migration rate, which allows for simultaneous coalescence or migration events in the genealogy of the sampled individuals. Comparisons are made with limit cases as the deme size or the number of demes goes to infinity or the migration rate goes to zero for which a diffusion approximation approach is possible. Conditions for selection to favor a mutant strategy replacing a resident strategy in the context of a linear game in a finite island population are addressed.     

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.     

Genealogy of Neutral Genes and Spreading of Selected Mutations in a Geographically Structured Population

In a geographically structured population, the interplay among gene migration, genetic drift and natural selection raises intriguing evolutionary problems, but the rigorous mathematical treatment is often very difficult. Therefore several approximate formulas were developed concerning the coalescence process of neutral genes and the fixation process of selected mutations in an island model, and their accuracy was examined by computer simulation. When migration is limited, the coalescence (or divergence) time for sampled neutral genes can be described by the convolution of exponential functions, as in a panmictic population, but it is determined mainly by migration rate and the number of demes from which the sample is taken. This time can be much longer than that in a panmictic population with the same number of breeding individuals. For a selected mutation, the spreading over the entire population was formulated as a birth and death process, in which the fixation probability within a deme plays a key role. With limited amounts of migration, even advantageous mutations take a large number of generations to spread. Furthermore, it is likely that these mutations which are temporarily fixed in some demes may be swamped out again by non-mutant immigrants from other demes unless selection is strong enough. These results are potentially useful for testing quantitatively various hypotheses that have been proposed for the origin of modern human populations.     

Polymorphism and divergence for island-model species

Estimates of the scaled selection coefficient, gamma of Sawyer and Hartl, are shown to be remarkably robust to population subdivision. Estimates of mutation parameters and divergence times, in contrast, are very sensitive to subdivision. These results follow from an analysis of natural selection and genetic drift in the island model of subdivision in the limit of a very large number of subpopulations, or demes. In particular, a diffusion process is shown to hold for the average allele frequency among demes in which the level of subdivision sets the timescale of drift and selection and determines the dynamic equilibrium of allele frequencies among demes. This provides a framework for inference about mutation, selection, divergence, and migration when data are available from a number of unlinked nucleotide sites. The effects of subdivision on parameter estimates depend on the distribution of samples among demes. If samples are taken singly from different demes, the only effect of subdivision is in the rescaling of mutation and divergence-time parameters. If multiple samples are taken from one or more demes, high levels of within-deme relatedness lead to low levels of intraspecies polymorphism and increase the number of fixed differences between samples from two species. If subdivision is ignored, mutation parameters are underestimated and the species divergence time is overestimated, sometimes quite drastically. Estimates of the strength of selection are much less strongly affected and always in a conservative direction.     

Multilocus selection in subdivided populations I. Convergence properties for weak or strong migration

The dynamics and equilibrium structure of a deterministic population-genetic model of migration and selection acting on multiple multiallelic loci is studied. A large population of diploid individuals is distributed over finitely many demes connected by migration. Generations are discrete and nonoverlapping, migration is irreducible and aperiodic, all pairwise recombination rates are positive, and selection may vary across demes. It is proved that, in the absence of selection, all trajectories converge at a geometric rate to a manifold on which global linkage equilibrium holds and allele frequencies are identical across demes. Various limiting cases are derived in which one or more of the three evolutionary forces, selection, migration, and recombination, are weak relative to the others. Two are particularly interesting. If migration and recombination are strong relative to selection, the dynamics can be conceived as a perturbation of the so-called weak-selection limit, a simple dynamical system for suitably averaged allele frequencies. Under nondegeneracy assumptions on this weak-selection limit which are generic, every equilibrium of the full dynamics is a perturbation of an equilibrium of the weak-selection limit and has the same stability properties. The number of equilibria is the same in both systems, equilibria in the full (perturbed) system are in quasi-linkage equilibrium, and differences among allele frequencies across demes are small. If migration is weak relative to recombination and epistasis is also weak, then every equilibrium is a perturbation of an equilibrium of the corresponding system without migration, has the same stability properties, and is in quasi-linkage equilibrium. In both cases, every trajectory converges to an equilibrium, thus no cycling or complicated dynamics can occur.      

Consequences of population structure on genes under balancing selection

This paper describes a new approach to modeling population structure for genes under strong balancing selection of the type seen in plant self-incompatibility systems and the major histocompatibility complex (MHC) system of vertebrates. Simple analytic solutions for the number of alleles maintained at equilibrium and the expected proportion of alleles shared between demes at various levels are derived and checked against simulation results. The theory accurately captures the dynamics of allele number in a subdivided population and identifies important values of m (migration rate) at which allele number and distribution change qualitatively. Starting from a panmictic population, as migration among demes decreases a qualitative change in dynamics is seen at approximately m(crit) approximately equal to the square root of(s/4piNT) where NT is the total population size and s is a measure of the strength of selection. At this point, demes can no longer maintain their panmictic allele number, due to increasing isolation from the total population. Another qualitative change occurs at a migration rate on the same order of magnitude as the mutation rate, mu. At this point, the demes are highly differentiated for allele complement, and the total number of alleles in the population is increased. Because in general u < m<(crit) at intermediate migration rates slightly fewer alleles may be maintained in the total population than are maintained at panmixia. Within this range, total allele number may not be the best indicator of whether a population is effectively panmictic, and some caution should be used when interpreting samples from such populations. The theory presented here can help to analyze data from genes under balancing selection in subdivided populations.     

Stochasticity, complex spatial structure, and the feasibility of the shifting balance theory

Sewall Wright's shifting balance theory of evolution posits a mechanism by which a structured population may escape local fitness optima and find a global optimum. We examine a one-locus, two-allele model of underdominance in populations with differing spatial arrangements of demes, both analytically and with Monte Carlo simulations. We find that inclusion of variance in interpatch connectivities can significantly reduce the number of generations required for fixation of the more favorable allele relative to island and stepping-stone models. Although time to fixation increases with migration rate in all cases, the presence of one or two relatively isolated demes may reduce the number of generations by 80% or more. These results suggest that the shifting balance process may operate under less restrictive conditions than those found with a simple spatial arrangement of demes.     

Quantifying the Role of Population Subdivision in Evolution on Rugged Fitness Landscapes

Natural selection drives populations towards higher fitness, but crossing fitness valleys or plateaus may facilitate progress up a rugged fitness landscape involving epistasis. We investigate quantitatively the effect of subdividing an asexual population on the time it takes to cross a fitness valley or plateau. We focus on a generic and minimal model that includes only population subdivision into equivalent demes connected by global migration, and does not require significant size changes of the demes, environmental heterogeneity or specific geographic structure. We determine the optimal speedup of valley or plateau crossing that can be gained by subdivision, if the process is driven by the deme that crosses fastest. We show that isolated demes have to be in the sequential fixation regime for subdivision to significantly accelerate crossing. Using Markov chain theory, we obtain analytical expressions for the conditions under which optimal speedup is achieved: valley or plateau crossing by the subdivided population is then as fast as that of its fastest deme. We verify our analytical predictions through stochastic simulations. We demonstrate that subdivision can substantially accelerate the crossing of fitness valleys and plateaus in a wide range of parameters extending beyond the optimal window. We study the effect of varying the degree of subdivision of a population, and investigate the trade-off between the magnitude of the optimal speedup and the width of the parameter range over which it occurs. Our results, obtained for fitness valleys and plateaus, also hold for weakly beneficial intermediate mutations. Finally, we extend our work to the case of a population connected by migration to one or several smaller islands. Our results demonstrate that subdivision with migration alone can significantly accelerate the crossing of fitness valleys and plateaus, and shed light onto the quantitative conditions necessary for this to occur.     

Separation of time scales, fixation probabilities and convergence to evolutionarily stable states under isolation by distance

To a first order of approximation, selection is frequency independent in a wide range of family structured models and in populations following an island model of dispersal, provided the number of families or demes is large and the population is haploid or diploid but allelic effects on phenotype are semidominant. This result underlies the way the evolutionary stability of traits is computed in games with continuous strategy sets. In this paper similar results are derived under isolation by distance. The first-order effect on expected change in allele frequency is given in terms of a measure of local genetic diversity, and of measures of genetic structure which are almost independent of allele frequency in the total population when the number of demes is large. Hence, when the number of demes increases the response to selection becomes of constant sign. This result holds because the relevant neutral measures of population structure converge to equilibrium at a rate faster than the rate of allele frequency changes in the total population. In the same conditions and in the absence of demographic fluctuations, the results also provide a simple way to compute the fixation probability of mutants affecting various ecological traits, such as sex ratio, dispersal, life-history, or cooperation, under isolation by distance. This result is illustrated and tested against simulations for mutants affecting the dispersal probability under a stepping-stone model.     

Separation of time scales, fixation probabilities and convergence to evolutionary stable states under isolation by distance

To a first order of approximation, selection is frequency independent in a wide range of family structured models and in populations following an island model of dispersal, provided the number of families or demes is large and the population is haploid or diploid but allelic effects on phenotype are semidominant. This result underlies the way the evolutionary stability of traits is computed in games with continuous strategy sets. In this paper similar results are derived under isolation by distance. The first-order effect on expected change in allele frequency is given in terms of a measure of local genetic diversity, and of measures of genetic structure which are almost independent of allele frequency in the total population when the number of demes is large. Hence, when the number of demes increases the response to selection becomes of constant sign. This result holds because the relevant neutral measures of population structure converge to equilibrium at a rate faster than the rate of allele frequency changes in the total population. In the same conditions and in the absence of demographic fluctuations, the results also provide a simple way to compute the fixation probability of mutants affecting various ecological traits, such as sex ratio, dispersal, life-history, or cooperation, under isolation by distance. This result is illustrated and tested against simulations for mutants affecting the dispersal probability under a stepping-stone model.     

The Effects of Background and Interference Selection on Patterns of Genetic Variation in Subdivided Populations

It is well known that most new mutations that affect fitness exert deleterious effects and that natural populations are often composed of subpopulations (demes) connected by gene flow. To gain a better understanding of the joint effects of purifying selection and population structure, we focus on a scenario where an ancestral population splits into multiple demes and study neutral diversity patterns in regions linked to selected sites. In the background selection regime of strong selection, we first derive analytic equations for pairwise coalescent times and F_ST as a function of time after the ancestral population splits into two demes and then construct a flexible coalescent simulator that can generate samples under complex models such as those involving multiple demes or nonconservative migration. We have carried out extensive forward simulations to show that the new methods can accurately predict diversity patterns both in the nonequilibrium phase following the split of the ancestral population and in the equilibrium between mutation, migration, drift, and selection. In the interference selection regime of many tightly linked selected sites, forward simulations provide evidence that neutral diversity patterns obtained from both the nonequilibrium and equilibrium phases may be virtually indistinguishable for models that have identical variance in fitness, but are nonetheless different with respect to the number of selected sites and the strength of purifying selection. This equivalence in neutral diversity patterns suggests that data collected from subdivided populations may have limited power for differentiating among the selective pressures to which closely linked selected sites are subject.     

Fixation probability in spatially changing environments.

The fixation probability of a mutant in a subdivided population with spatially varying environments is investigated using a finite island model. This probability is different from that in a panmictic population if selection is intermediate to strong and migration is weak. An approximation is used to compute the fixation probability when migration among subpopulations is very weak. By numerically solving the two-dimensional partial differential equation for the fixation probability in the two subpopulation case, the approximation was shown to give fairly accurate values. With this approximation, we show in the case of two subpopulations that the fixation probability in subdivided populations is greater than that in panmictic populations mostly. The increase is most pronounced when the mutant is selected for in one subpopulation and is selected against in the other subpopulation. Also it is shown that when there are two types of environments, further subdivision of subpopulations does not cause much change of the fixation probability in the no dominance case unless the product of the selection coefficient and the local population size is less than one. With dominance, the effect of subdivision becomes more complex.     

Migration promotes mutator alleles in subdivided populations

Mutator alleles that elevate the genomic mutation rate may invade nonrecombining populations by hitchhiking with beneficial mutations. Mutators have been repeatedly observed to take over adapting laboratory populations and have been found at high frequencies in both microbial pathogen and cancer populations in nature. Recently, we have shown that mutators are only favored by selection in sufficiently large populations and transition to being disfavored as population size decreases. This population size‐dependent sign inversion in selective effect suggests that population structure may also be an important determinant of mutation rate evolution. Although large populations may favor mutators, subdividing such populations into sufficiently small subpopulations (demes) might effectively inhibit them. On the other hand, migration between small demes that otherwise inhibit hitchhiking may promote mutator fixation in the whole metapopulation. Here, we use stochastic, agent‐based simulations and evolution experiments with the yeast Saccharomyces cerevisiae to show that mutators can, indeed, be favored by selection in subdivided metapopulations composed of small demes connected by sufficient migration. In fact, we show that population structure plays a previously unsuspected role in promoting mutator success in subdivided metapopulations when migration is rare.     

The maintenance of sexual reproduction in a structured population

We present the results of a computer simulation model in which a sexual population produces an asexual mutant. We estimate the probability that the new asexual lineage will go extinct. We find that whenever the asexual lineage does not go extinct the sexual population is out-competed, and only asexual individuals remain after a sufficiently long period of time has elapsed. We call this type of outcome an asexual takeover. Our results suggest that, given repeated mutations to asexuality, asexual takeover is likely in an unstructured environment. However, if the environment is subdivided into demes that are connected by migration, then asexual takeover becomes less likely. The probability of asexual takeover declines towards zero as the number of demes increases and as the rate of migration decreases. The reason for this is that asexuality leads to a greater loss of fitness due to mutation and genetic drift, in comparison to what occurs under sexual reproduction. Population subdivision slows the spread of asexual lineages, which allows more time for the genetic degeneration caused by asexuality to take place.     

The effects of subdivision on the genetic divergence of populations and species

Abstract. An island model of migration is used to study the effects of subdivision within populations and species on sample genealogies and on between-population or between-species measures of genetic variation. The model assumes that the number of demes within each population or species is large. When populations (or species), connected either by gene flow or historical association, are themselves subdivided into demes, changes in the migration rate among demes alter both the structure of genealogies and the time scale of the coalescent process. The time scale of the coalescent is related to the effective size of the population, which depends on the migration rate among demes. When the migration rate among demes within populations is low, isolation (or speciation) events seem more recent and migration rates among populations seem higher because the effective size of each population is increased. This affects the probability of reciprocal monophyly of two samples, the chance that a gene tree of a sample matches the species tree, and relative likelihoods of different types of polymorphic sites. It can also have a profound effect on the estimation of divergence times.     

Linked selected and neutral loci in heterogeneous environments

We analyze a system of ordinary differential equations modeling haplotype frequencies at a physically linked pair of loci, one selected and one neutral, in a population consisting of two demes with divergent selection regimes. The system is singularly perturbed, with the migration rate m between the demes serving as a small parameter. We use geometric singular perturbation theory to show that when m is sufficiently small, each solution not initially fixed for the same selected allele in both demes approaches one of a 1-dimensional continuum of equilibria. We then obtain asymptotic expansions of the solutions and show their validity on arbitrarily long finite time intervals. From these expansions we obtain formulas for the transient dynamics of F _ST (a measure of population structure) at both loci, as well as for the rate of genotyping error if the allelic state at the selected locus is inferred from that at the neutral (marker) locus. We examine two cases in detail, one modeling two populations in secondary contact after a period of evolution in allopatry, and the other modeling the origination and spread of a resistance allele.Electronic supplementary material Supplementary material is available in the online version of this article at http://dx.doi.org/10.1007/s00285-006-0038-6 and is accessible for authorized users.     

The effect of hitch-hiking on genes linked to a balanced polymorphism in a subdivided population

The effect of multi-allelic balancing selection on nucleotide diversity at linked neutral sites was investigated by simulations of subdivided populations. The motivation is to understand the behaviour of self-recognition systems such as the MHC and plant self-incompatibility. For neutral sites, two types of subdivision are present: (1) into demes (connected by migration), and (2) into classes defined by different functional alleles at the selected locus (connected by recombination). Previous theoretical studies of each type of subdivision separately have shown that each increases diversity, and decreases the relative frequencies of low-frequency variants, at neutral sites or loci. We show here that the two types of subdivision act non-additively when sampling is at the whole population level, and that subdivision produces some non-intuitive results. For instance, in highly subdivided populations, genetic diversity at neutral sites may decrease with tighter linkage to a selected locus or site. Another conclusion is that, if there is population subdivision, balancing selection leads to decreased expected FST values for neutral sites linked to the selected locus. Finally, we show that the ability to detect balancing selection by its effects on linked variation, using tests such as Tajima's D, is reduced when genes in a subdivided population are sampled from the total population, rather than within demes.     

Multilocus selection in subdivided populations II. Maintenance of polymorphism under weak or strong migration

The potential of maintaining multilocus polymorphism by migration-selection balance is studied. A large population of diploid individuals is distributed over finitely many demes connected by migration. Generations are discrete and nonoverlapping, selection may vary across demes, and loci are multiallelic. It is shown that if migration and recombination are strong relative to selection, then with weak or no epistasis and intermediate dominance at every locus and in every deme, arbitrarily many alleles can be maintained at arbitrarily many loci at a stable equilibrium. If migration is weak relative to selection and recombination, then with weak or no epistasis and intermediate dominance at every locus and in every deme, as many alleles as there are demes can be maintained at arbitrarily many loci at equilibrium. In both cases open sets of such parameter combinations are constructed, thus the results are robust with respect to small, but arbitrary, perturbations in the parameters. For weak migration, the number of demes is, in fact, a generic upper bound to the number of alleles that can be maintained at any locus. Thus, several scenarios are identified under which multilocus polymorphism can be maintained by migration-selection balance when this is impossible in a panmictic population.      

Polymorphism in multiallelic migration–selection models with dominance

The evolution of the gene frequencies at a single multiallelic locus under the joint action of migration and viability selection with dominance is investigated. The monoecious, diploid population is subdivided into finitely many panmictic colonies that exchange adult migrants independently of genotype. Underdominance and overdominance are excluded. If the degree of dominance is deme independent for every pair of alleles, then under the Levene model, the qualitative evolution of the gene frequencies (i.e., the existence and stability of the equilibria) is the same as without dominance. In particular: (i) the number of demes is a generic upper bound on the number of alleles present at equilibrium; (ii) there exists exactly one stable equilibrium, and it is globally attracting; and (iii) if there exists an internal equilibrium, it is globally asymptotically stable. Analytic examples demonstrate that if either the Levene model does not apply or the degree of dominance is deme dependent, then the above results can fail. A complete global analysis of weak migration and weak selection on a recessive allele in two demes is presented.