Selection in a subdivided population with local extinction and recolonization

In a subdivided population, local extinction and subsequent recolonization affect the fate of alleles. Of particular interest is the interaction of this force with natural selection. The effect of selection can be weakened by this additional source of stochastic change in allele frequency. The behavior of a selected allele in such a population is shown to be equivalent to that of an allele with a different selection coefficient in an unstructured population with a different size. This equivalence allows use of established results for panmictic populations to predict such quantities as fixation probabilities and mean times to fixation. The magnitude of the quantity N(e)s(e), which determines fixation probability, is decreased by extinction and recolonization. Thus deleterious alleles are more likely to fix, and advantageous alleles less likely to do so, in the presence of extinction and recolonization. Computer simulations confirm that the theoretical predictions of both fixation probabilities and mean times to fixation are good approximations.     

Effect of dominance on heterozygosity and the fixation probability in a subdivided population

Under the assumptions of a subdivided population and the presence of dominance for fitness, the expected sum of heterozygosity in the total population during the lifetime of mutant was investigated. It was shown analytically and by computer simulations that in the island model the effect of dominance on the expected sum of heterozygosity decreases as the migration rate decreases and is lost almost completely when the migration rate is very low. In addition to the expected sum of heterozygosity, the fixation probability of mutant was also investigated. The effect of dominance on the fixation probability also decreases as the migration rate decreases but is not completely lost when the migration rate is very low.     

Gene frequency patterns in the levene subdivided population model

In the context of a multideme population structure subject to selection, migration, and mating forces, it is desired to ascertain the stable evolutionary outcomes for various classes of selection regimes. These include selection patterns as (i) a mosaic of local directed selection effects, (ii) overdominance throughout with varying intensities of the local heterozygote advantage, (iii) varying degrees of underdominance throughout, (iv) a mixed underdominant-overdominant regime. An accounting of the nature of the equilibrium configurations in the Levene population subdivision model was done with respect to the above classes of selection regimes. In particular, it is established that multiple polymorphic equilibria do not arise for selection structures (i) and (ii), while for the mixed underdominant-overdominant selection form (iv) with appropriate parameter ranges there can exist two stable internal equilibria. A surprising finding is that the number and/or character of the equilibria is not changed by increased population division beyond that of two habitats, while there is a significant difference in the equilibrium possibilities for a one- as against a two-deme population. These results appear to be a special limiting feature of the Levene population subdivision formulation.     

Selection and drift in subdivided populations: a straightforward method for deriving diffusion approximations and applications involving dominance, selfing and local extinctions

Population structure affects the relative influence of selection and drift on the change in allele frequencies. Several models have been proposed recently, using diffusion approximations to calculate fixation probabilities, fixation times, and equilibrium properties of subdivided populations. We propose here a simple method to construct diffusion approximations in structured populations; it relies on general expressions for the expectation and variance in allele frequency change over one generation, in terms of partial derivatives of a "fitness function" and probabilities of genetic identity evaluated in a neutral model. In the limit of a very large number of demes, these probabilities can be expressed as functions of average allele frequencies in the metapopulation, provided that coalescence occurs on two different timescales, which is the case in the island model. We then use the method to derive expressions for the probability of fixation of new mutations, as a function of their dominance coefficient, the rate of partial selfing, and the rate of deme extinction. We obtain more precise approximations than those derived by recent work, in particular (but not only) when deme sizes are small. Comparisons with simulations show that the method gives good results as long as migration is stronger than selection.     

Evolution of coadaptation in a subdivided population

The interplay between population subdivision and epistasis is investigated by studying the fixation probability of a coadapted haplotype in a subdivided population. Analytical and simulation models are developed to study the evolutionary fate of two conditionally neutral mutations that interact epistatically to enhance fitness. We find that the fixation probability of a coadapted haplotype shows a marked increase when the population is genetically subdivided and subpopulations are loosely connected by migration. Moderate migration and isolation allow the propagation of the mutant alleles across subpopulations, while at the same time preserving the favorable allelic combination established within each subpopulation. Together they create the condition most favorable for the ultimate fixation of the coadapted haplotype. On the basis of the analytical and simulation results, we discuss the fundamental role of population subdivision and restricted gene flow in promoting the evolution of functionally integrated systems, with some implications for the shifting-balance theory of evolution.     

Match probabilities in a finite, subdivided population

We generalize a recently introduced graphical framework to compute the probability that haplotypes or genotypes of two individuals drawn from a finite, subdivided population match. As in the previous work, we assume an infinite-alleles model. We focus on the case of a population divided into two subpopulations, but the underlying framework can be applied to a general model of population subdivision. We examine the effect of population subdivision on the match probabilities and the accuracy of the product rule which approximates multi-locus match probabilities as a product of one-locus match probabilities. We quantify the deviation from predictions of the product rule by R, the ratio of the multi-locus match probability to the product of the one-locus match probabilities. We carry out the computation for two loci and find that ignoring subdivision can lead to underestimation of the match probabilities if the population under consideration actually has subdivision structure and the individuals originate from the same subpopulation. On the other hand, under a given model of population subdivision, we find that the ratio R for two loci is only slightly greater than 1 for a large range of symmetric and asymmetric migration rates. Keeping in mind that the infinite-alleles model is not the appropriate mutation model for STR loci, we conclude that, for two loci and biologically reasonable parameter values, population subdivision may lead to results that disfavor innocent suspects because of an increase in identity-by-descent in finite populations. On the other hand, for the same range of parameters, population subdivision does not lead to a substantial increase in linkage disequilibrium between loci. Those results are consistent with established practice.     

Inference from gene trees in a subdivided population

This paper studies gene trees in subdivided populations which are constructed as perfect phylogenies from the pattern of mutations in a sample of DNA sequences and presents a new recursion for the probability distribution of such gene trees. The underlying evolutionary model is the coalescent process in a subdivided population. The infinitely-many-sites model of mutation is assumed. Ancestral inference questions that are discussed are maximum likelihood estimation of migration and mutation rates; detection of population growth by likelihood techniques; determining the distribution of the time to the most recent common ancestor of a sample of sequences; determining the distribution of the age of the mutations on the gene tree; determining in which subpopulation the most recent common ancestor of all the sequences was; determining subpopulation ancestors, where they were, and times to them; and determining in which subpopulations mutations occurred. A computational technique of Griffiths and Tavaré used is a computer intensive Markov chain simulation, which simulates gene trees conditional on their topology implied by the mutation pattern in the sample of DNA sequences. The software GENETREE, which implements these ancestral inference techniques, is available.     

Variation after a selective sweep in a subdivided population

The effect of genetic hitchhiking on neutral variation is analyzed in subdivided populations with differentiated demes. After fixation of a favorable mutation, the consequences on particular subpopulations can be radically different. In the subpopulation where the mutation first appeared by mutation, variation at linked neutral loci is expected to be reduced, as predicted by the classical theory. However, the effect in the other subpopulations, where the mutation is introduced by migration, can be the opposite. This effect depends on the level of genetic differentiation of the subpopulations, the selective advantage of the mutation, the recombination frequency, and the population size, as stated by analytical derivations and computer simulations. The characteristic outcomes of the effect are three. First, the genomic region of reduced variation around the selected locus is smaller than that predicted in a panmictic population. Second, for more distant neutral loci, the amount of variation increases over the level they had before the hitchhiking event. Third, for these loci, the spectrum of gene frequencies is dominated by an excess of alleles at intermediate frequencies when compared with the neutral theory. At these loci, hitchhiking works like a system that takes variation from the between-subpopulation component and introduces it into the subpopulations. The mechanism can also operate in other systems in which the genetic variation is distributed in clusters with limited exchange of variation, such as chromosome arrangements or genomic regions closely linked to targets of balancing selection.     

The propagation of a cultural or biological trait by neutral genetic drift in a subdivided population

We study fixation probabilities and times as a consequence of neutral genetic drift in subdivided populations, motivated by a model of the cultural evolutionary process of language change that is described by the same mathematics as the biological process. We focus on the growth of fixation times with the number of subpopulations, and variation of fixation probabilities and times with initial distributions of mutants. A general formula for the fixation probability for arbitrary initial condition is derived by extending a duality relation between forwards- and backwards-time properties of the model from a panmictic to a subdivided population. From this we obtain new formulae(formally exact in the limit of extremely weak migration) for the mean fixation time from an arbitrary initial condition for Wright's island model, presenting two cases as examples. For more general models of population subdivision, formulae are introduced for an arbitrary number of mutants that are randomly located, and a single mutant whose position is known. These formulae contain parameters that typically have to be obtained numerically, a procedure we follow for two contrasting clustered models. These data suggest that variation of fixation time with the initial condition is slight, but depends strongly on the nature of subdivision. In particular, we demonstrate conditions under which the fixation time remains finite even in the limit of an infinite number of demes. In many cases-except this last where fixation in a finite time is seen--the time to fixation is shown to be in precise agreement with predictions from formulae for the asymptotic effective population size.     

Persistence of an infectious disease in a subdivided population

The transmission dynamics of a communicable disease in a subdivided population where the spread among groups follows the proportionate mixing model while the within-group transmission can correspond to preferred mixing, proportionate mixing among subgroups, or mixing between social and nonsocial subgroups, is analyzed. It is shown that the threshold condition for the disease to persist is that either (i) the disease can persist within at least one group through intragroup contacts, or--if (i) does not hold--(ii) the intergroup transmission is sufficiently high. The among-group transmission is computed as an average where each subgroup's reproductive number is weighted according to its intragroup activity level squared and the total number of cases that one infectious individual will cause through intragroup contacts. The model thus allows for a study of the relative importance of communitywide disease transmission and of disease transmission within geographically or socially separate groups.     

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.     

The evolution of social communication systems in a subdivided population

Evolution of social communication systems is modeled with a quantitative genetic model. The mathematical model describes the coevolutionary process of a social signal (a social character) and responsiveness (a social preference) to the signal. The responsiveness is postulated to influence fitness of senders of the signal. Considerations are extended to subdivided population structure by combining the social selection model with a group selection model. The numerical results derived from the models indicate that the evolutionary rate of social communication systems depends largely on genetic correlation between the signal and the responsiveness. Group selection can reinforce the evolutionary rate and relax its dependence on genetic correlation. The origin of genetic correlation is discussed in relation to group selection.     

Coalescence time for two genes from a subdivided population

In this paper a new form of the solution for the Laplace transform and moments of the distribution of the waiting time for two genes to coalescence is presented. The two genes are sampled from a subdivided population where migration rates between populations are constant in time. Equal subpopulation size is not assumed. For the special case of an island model with equal migration rates between islands, the Laplace transform of the coalescence time and the first and second moments are found explicitly. The new form of the solutions allows numerical calculation. The connection of how the results relate to a panmictic population when migration rates are large is illustrated using strong-migration-limit theory. Received: 19 April 1999 / Revised version: 22 March 2001 / Published online: 19 September 2001     

The expected number of heterozygous sites in a subdivided population.

A simple, exact formula is derived for the expected number of heterozygous sites per individual at equilibrium in a subdivided population. The model of infinitely many neutral sites is posited; the linkage map is arbitrary. The monoecious, diploid population is subdivided into a finite number of panmictic colonies that exchange gametes. The backward migration matrix is arbitrary, but time independent and ergodic (i.e., irreducible and aperiodic). With suitable weighting, the expected number of heterozygous sites is 4Neu, where Ne denotes the migration effective population number and u designates the total mutation rate per gene (or DNA sequence). For diploid migration, this formula is a good approximation if Ne >> 1.     

A diffusion approximation for selection and drift in a subdivided population

The population-genetic consequences of population structure are of great interest and have been studied extensively. An area of particular interest is the interaction among population structure, natural selection, and genetic drift. At first glance, different results in this area give very different impressions of the effect of population subdivision on effective population size (N(e)), suggesting that no single value of N(e) can completely characterize a structured population. Results presented here show that a population conforming to Wright's island model of subdivision with genic selection can be related to an idealized panmictic population (a Wright-Fisher population). This equivalent panmictic population has a larger size than the actual population; i.e., N(e) is larger than the actual population size, as expected from many results for this type of population structure. The selection coefficient in the equivalent panmictic population, referred to here as the effective selection coefficient (s(e)), is smaller than the actual selection coefficient (s). This explains how the fixation probability of a selected allele can be unaffected by population subdivision despite the fact that subdivision increases N(e), for the product N(e)s(e) is not altered by subdivision.     

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.