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.     

Kin groups and trait groups: population structure and epidemic disease selection

A Monte Carlo simulation based on the population structure of a small-scale human population, the Semai Senoi of Malaysia, has been developed to study the combined effects of group, kin, and individual selection. The population structure resembles D.S. Wilson's structured deme model in that local breeding populations (Semai settlements) are subdivided into trait groups (hamlets) that may be kin-structured and are not themselves demes. Additionally, settlement breeding populations are connected by two-dimensional stepping-stone migration approaching 30% per generation. Group and kin-structured group selection occur among hamlets the survivors of which then disperse to breed within the settlement population. Genetic drift is modeled by the process of hamlet formation; individual selection as a deterministic process, and stepping-stone migration as either random or kin-structured migrant groups. The mechanism for group selection is epidemics of infectious disease that can wipe out small hamlets particularly if most adults become sick and social life collapses. Genetic resistance to a disease is an individual attribute; however, hamlet groups with several resistant adults are less likely to disintegrate and experience high social mortality. A specific human gene, hemoglobin E, which confers resistance to malaria, is studied as an example of the process. The results of the simulations show that high genetic variance among hamlet groups may be generated by moderate degrees of kin-structuring. This strong microdifferentiation provides the potential for group selection. The effect of group selection in this case is rapid increase in gene frequencies among the total set of populations. In fact, group selection in concert with individual selection produced a faster rate of gene frequency increase among a set of 25 populations than the rate within a single unstructured population subject to deterministic individual selection. Such rapid evolution with plausible rates of extinction, individual selection, and migration and a population structure realistic in its general form, has implications for specific human polymorphisms such as hemoglobin variants and for the more general problem of the tempo of evolution as well.     

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.     

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.     

Population structure promotes the evolution of costly sex in artificial gene networks

We build on previous observations that Hill–Robertson interference generates an advantage of sex that, in structured populations, can be large enough to explain the evolutionary maintenance of costly sex. We employed a gene network model that explicitly incorporates interactions between genes. Mutations in the gene networks have variable effects that depend on the genetic background in which they appear. Consequently, our simulations include two costs of sex—recombination and migration loads—that were missing from previous studies of the evolution of costly sex. Our results suggest a critical role for population structure that lies in its ability to align the long‐ and short‐term advantages of sex. We show that the addition of population structure favored the evolution of sex by disproportionately decreasing the equilibrium mean fitness of asexual populations, primarily by increasing the strength of Muller's Ratchet. Population structure also increased the ability of the short‐term advantage of sex to counter the primary limit to the evolution of sex in the gene network model—recombination load. On the other hand, highly structured populations experienced migration load in the form of Dobzhansky–Muller incompatibilities, decreasing the effective rate of migration between demes and, consequently, accelerating the accumulation of drift load in the sexual populations.     

Migration and recovery of the genetic diversity during the increasing density phase in cyclic vole populations

In cyclic populations, high genetic diversity is currently reported despite the periodic low numbers experienced by the populations during the low phases. Here, we report spatio-temporal monitoring at a very fine scale of cyclic populations of the fossorial water vole (Arvicola terrestris) during the increasing density phase. This phase marks the transition from a patchy structure (demes) during low density to a continuous population in high density. We found that the genetic diversity was effectively high but also that it displayed a local increase within demes over the increasing phase. The genetic diversity remained relatively constant when considering all demes together. The increase in vole abundance was also correlated with a decrease of genetic differentiation among demes. Such results suggest that at the end of the low phase, demes are affected by genetic drift as the result of being small and geographically isolated. This leads to a loss of local genetic diversity and a spatial differentiation among demes. This situation is counterbalanced during the increasing phase by the spatial expansion of demes and the increase of the effective migration among differentiated demes. We provide evidences that in cyclic populations of the fossorial water voles, the relative influence of drift operating during low density populations and migration occurring principally while population size increases interacts closely to maintain high genetic diversity.     

Muller''s Ratchet under Epistatic Selection

In a finite asexual population mean fitness may decrease by a process known as Muller's ratchet, which proceeds if all individuals with the minimum number of deleterious alleles are randomly lost. If these alleles have independent effects on fitness, previous analysis suggested that the rate of this decrease either remains constant or, if accumulation of mutations leads to the decline of the population size, grows. Here I show that this conclusion is quite sensitive to the assumption of independence. If deleterious alleles have synergistic fitness effects, then, as the ratchet advances, the frequency of the best available genotype will necessarily increase, making its loss less and less probable. As a result, sufficiently strong synergistic epistasis can effectively halt the action of Muller's ratchet. Instead of being driven extinct, a finite asexual population could then survive practically indefinitely, although with lower mean fitness than without random drift.     

Muller''s Ratchet, Epistasis and Mutation Effects

In this study, computer simulation is used to show that despite synergistic epistasis for fitness, Muller's ratchet can lead to lethal fitness loss in a population of asexuals through the accumulation of deleterious mutations. This result contradicts previous work that indicated that epistasis will halt the ratchet. The present results show that epistasis will not halt the ratchet provided that rather than a single deleterious mutation effect, there is a distribution of deleterious mutation effects with sufficient density near zero. In addition to epistasis and mutation distribution, the ability of Muller's ratchet to lead to the extinction of an asexual population under epistasis for fitness depends strongly on the expected number of offspring that survive to reproductive age. This strong dependence is not present in the nonepistatic model and suggests that interpreting the population growth parameter as fecundity is inadequate. Because a continuous distribution of mutation effects is used in this model, an emphasis is placed on the dynamics of the mutation effect distribution rather than on the dynamics of the number of least mutation loaded individuals. This perspective suggests that current models of gene interaction are too simple to apply directly to long-term prediction for populations undergoing the ratchet.     

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.     

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.     

Convergence to the island-model coalescent process in populations with restricted migration

In this article we apply some graph-theoretic results to the study of coalescence in a structured population with migration. The graph is the pattern of migration among subpopulations, or demes, and we use the theory of random walks on graphs to characterize the ease with which ancestral lineages can traverse the habitat in a series of migration events. We identify conditions under which the coalescent process in populations with restricted migration, such that individuals cannot traverse the habitat freely in a single migration event, nonetheless becomes identical to the coalescent process in the island migration model in the limit as the number of demes tends to infinity. Specifically, we first note that a sequence of symmetric graphs with Diaconis-Stroock constant bounded above has an unstructured Kingman-type coalescent in the limit for a sample of size two from two different demes. We then show that circular and toroidal models with long-range but restricted migration have an upper bound on this constant and so have an unstructured-migration coalescent in the limit. We investigate the rate of convergence to this limit using simulations.     

Detecting Isolation by Distance Using Phylogenies of Genes

We introduce a method for analyzing phylogenies of genes sampled from a geographically structured population. A parsimony method can be used to compute s, the minimum number of migration events between pairs of populations sampled, and the value of s can be used to estimate the effective migration rate M, the value of Nm in an island model with local populations of size N and a migration rate m that would yield the same value of s. Extensive simulations show that there is a simple relationship between M and the geographic distance between pairs of samples in one- and two-dimensional models of isolation by distance. Both stepping-stone and lattice models were simulated. If two demes k steps apart are sampled, then, s, the average value of s, is a function only of k/(Nm) in a one-dimensional model and is a function only of k/(Nm)2 in a two-dimensional model. Furthermore, log(M) is approximately a linear function of log(k). In a one-dimensional model, the regression coefficient is approximately -1 and in a two-dimensional model the regression coefficient is approximately -0.5. Using data from several locations, the regression of log(M) on log(distance) may indicate whether there is isolation by distance in a population at equilibrium and may allow an estimate of the effective migration rate between adjacent sampling locations. Alternative methods for analyzing DNA sequence data from a geographically structured population are discussed. An application of our method to the data of R. L. Cann, M. Stoneking and A. C. Wilson on human mitochondrial DNA is presented.     

Influence of finite-sites mutation, population subdivision and sampling schemes on patterns of nucleotide polymorphism for species with molecular hyperdiversity

Molecular hyperdiversity has been documented in viruses, prokaryotes and eukaryotes. Such organisms undermine the assumptions of the infinite-sites mutational model, because multiple mutational events at a site comprise a non-negligible portion of polymorphisms. Moreover, different sampling schemes of individuals from species with subdivided populations can profoundly influence resulting patterns and interpretations of molecular variation. Inspired by molecular hyperdiversity in the nematode Caenorhabditis sp. 5, which exhibits average pairwise differences among synonymous sites of >5% as well as modest population structure, we investigated via coalescent simulation the joint effects of a finite-sites mutation (FSM) process and population subdivision on the variant frequency spectrum. From many demes interconnected through a stepping-stone migration model, we constructed local samples from a single deme, pooled samples from several demes and scattered samples of a single individual from numerous demes. Compared with a single panmictic population at equilibrium, we find that high population mutation rates induce a deficit of rare variants (positive Tajima's D) under a FSM model. Population structure also induces such a skew for local samples when migration is high and for pooled samples when migration is low. Contrasts of sampling schemes for C. sp. 5 imply high mutational input coupled with high migration. We propose that joint analysis of local, pooled and scattered samples for species with subdivided populations provides a means of improving inference of demographic history, by virtue of the partially distinct patterns of polymorphism that manifest when sequences are analyzed according to differing sampling schemes.     

The temporal dynamics of processes underlying Y chromosome degeneration

Y chromosomes originate from ordinary autosomes and degenerate by accumulating deleterious mutations. This accumulation results from a lack of recombination on the Y and is driven by interference among deleterious mutations (Muller's ratchet and background selection) and the fixation of beneficial alleles (genetic hitchhiking). Here I show that the relative importance of these processes is expected to vary over the course of Y chromosome evolution due to changes in the number of active genes. The dominant mode of degeneration on a newly formed gene-rich Y chromosome is expected to be Muller's ratchet and/or background selection due to the large numbers of deleterious mutations arising in active genes. However, the relative importance of these modes of degeneration declines rapidly as active genes are lost. In contrast, the rate of degeneration due to hitchhiking is predicted to be highest on Y chromosomes containing an intermediate number of active genes. The temporal dynamics of these processes imply that a gradual restriction of recombination, as inferred in mammals, will increase the importance of genetic hitchhiking relative to Muller's ratchet and background selection.     

The structured ancestral selection graph and the many-demes limit

We show that the unstructured ancestral selection graph applies to part of the history of a sample from a population structured by restricted migration among subpopulations, or demes. The result holds in the limit as the number of demes tends to infinity with proportionately weak selection, and we have also made the assumptions of island-type migration and that demes are equivalent in size. After an instantaneous sample-size adjustment, this structured ancestral selection graph converges to an unstructured ancestral selection graph with a mutation parameter that depends inversely on the migration rate. In contrast, the selection parameter for the population is independent of the migration rate and is identical to the selection parameter in an unstructured population. We show analytically that estimators of the migration rate, based on pairwise sequence differences, derived under the assumption of neutrality should perform equally well in the presence of weak selection. We also modify an algorithm for simulating genealogies conditional on the frequencies of two selected alleles in a sample. This permits efficient simulation of stronger selection than was previously possible. Using this new algorithm, we simulate gene genealogies under the many-demes ancestral selection graph and identify some situations in which migration has a strong effect on the time to the most recent common ancestor of the sample. We find that a similar effect also increases the sensitivity of the genealogy to selection.     

Theory of the effects of population structure and sampling on patterns of linkage disequilibrium applied to genomic data from humans

We develop predictions for the correlation of heterozygosity and for linkage disequilibrium between two loci using a simple model of population structure that includes migration among local populations, or demes. We compare the results for a sample of size two from the same deme (a single-deme sample) to those for a sample of size two from two different demes (a scattered sample). The correlation in heterozygosity for a scattered sample is surprisingly insensitive to both the migration rate and the number of demes. In contrast, the correlation in heterozygosity for a single-deme sample is sensitive to both, and the effect of an increase in the number of demes is qualitatively similar to that of a decrease in the migration rate: both increase the correlation in heterozygosity. These same conclusions hold for a commonly used measure of linkage disequilibrium (r(2)). We compare the predictions of the theory to genomic data from humans and show that subdivision might account for a substantial portion of the genetic associations observed within the human genome, even though migration rates among local populations of humans are relatively large. Because correlations due to subdivision rather than to physical linkage can be large even in a single-deme sample, then if long-term migration has been important in shaping patterns of human polymorphism, the common practice of disease mapping using linkage disequilibrium in "isolated" local populations may be subject to error.     

Temporal variation in selection accelerates mutational decay by Muller's ratchet

Asexual species accumulate deleterious mutations through an irreversible process known as Muller's ratchet. Attempts to quantify the rate of the ratchet have ignored the role of temporal environmental heterogeneity even though it is common in nature and has the potential to affect overall ratchet rate. Here we examine Muller's ratchet in the context of conditional neutrality (i.e., mutations that are deleterious in some environmental conditions but neutral in others) as well as more subtle changes in the strength (but not sign) of selection. We find that temporal variation increases the rate of the ratchet (mutation accumulation) and the rate of fitness decline over that of populations experiencing constant selection of equivalent average strength. Temporal autocorrelation magnifies the effects of temporal heterogeneity and can allow the ratchet to operate at large population sizes in which it would be halted under constant selection. Classic studies of Muller's ratchet show that the rate of fitness decline is maximized when selection is of a low but intermediate strength. This relationship changes quantitatively with all forms of temporal heterogeneity studied and changes qualitatively when there is temporal autocorrelation in selection. In particular, the rate of fitness decline can increase indefinitely with the strength of selection with some forms of temporal heterogeneity. Our finding that temporal autocorrelation in selection dramatically increases ratchet rate and rate of fitness decline may help to explain the paucity of asexual taxa.     

Muller's ratchet and the degeneration of Y chromosomes: a simulation study

A typical pattern in sex chromosome evolution is that Y chromosomes are small and have lost many of their genes. One mechanism that might explain the degeneration of Y chromosomes is Muller's ratchet, the perpetual stochastic loss of linkage groups carrying the fewest number of deleterious mutations. This process has been investigated theoretically mainly for asexual, haploid populations. Here, I construct a model of a sexual population where deleterious mutations arise on both X and Y chromosomes. Simulation results of this model demonstrate that mutations on the X chromosome can considerably slow down the ratchet. On the other hand, a lower mutation rate in females than in males, background selection, and the emergence of dosage compensation are expected to accelerate the process.     

Dynamic mutation-selection balance as an evolutionary attractor

The vast majority of mutations are deleterious and are eliminated by purifying selection. Yet in finite asexual populations, purifying selection cannot completely prevent the accumulation of deleterious mutations due to Muller's ratchet: once lost by stochastic drift, the most-fit class of genotypes is lost forever. If deleterious mutations are weakly selected, Muller's ratchet can lead to a rapid degradation of population fitness. Evidently, the long-term stability of an asexual population requires an influx of beneficial mutations that continuously compensate for the accumulation of the weakly deleterious ones. Hence any stable evolutionary state of a population in a static environment must involve a dynamic mutation-selection balance, where accumulation of deleterious mutations is on average offset by the influx of beneficial mutations. We argue that such a state can exist for any population size N and mutation rate U and calculate the fraction of beneficial mutations, ε, that maintains the balanced state. We find that a surprisingly low ε suffices to achieve stability, even in small populations in the face of high mutation rates and weak selection, maintaining a well-adapted population in spite of Muller's ratchet. This may explain the maintenance of mitochondria and other asexual genomes.