Gene genealogies in a metapopulation

A simple genealogical process is found for samples from a metapopulation, which is a population that is subdivided into a large number of demes, each of which is subject to extinction and recolonization and receives migrants from other demes. As in the migration-only models studied previously, the genealogy of any sample includes two phases: a brief sample-size adjustment followed by a coalescent process that dominates the history. This result will hold for metapopulations that are composed of a large number of demes. It is robust to the details of population structure, as long as the number of possible source demes of migrants and colonists for each deme is large. Analytic predictions about levels of genetic variation are possible, and results for average numbers of pairwise differences within and between demes are given. Further analysis of the expected number of segregating sites in a sample from a single deme illustrates some previously known differences between migration and extinction/recolonization. The ancestral process is also amenable to computer simulation. Simulation results show that migration and extinction/recolonization have very different effects on the site-frequency distribution in a sample from a single deme. Migration can cause a U-shaped site-frequency distribution, which is qualitatively similar to the pattern reported recently for positive selection. Extinction and recolonization, in contrast, can produce a mode in the site-frequency distribution at intermediate frequencies, even in a sample from a single deme.     

The evolutionary dynamics of evolvability in a gene network model

Evolvability, the ability of populations to adapt, has recently emerged as a major unifying concept in biology. Although the study of evolvability offers new insights into many important biological questions, the conceptual bases of evolvability, and the mechanisms of its evolution, remain controversial. We used simulated evolution of a model of gene network dynamics to test the contentious hypothesis that natural selection can favour high evolvability, in particular in sexual populations. Our results conclusively demonstrate that fluctuating natural selection can increase the capacity of model gene networks to adapt to new environments. Detailed studies of the evolutionary dynamics of these networks establish a broad range of validity for this result and quantify the evolutionary forces responsible for changes in evolvability. Analysis of the genotype–phenotype map of these networks also reveals mechanisms connecting evolvability, genetic architecture and robustness. Our results suggest that the evolution of evolvability can have a pervasive influence on many aspects of organisms.     

Populational heritability: empirical studies of evolution in metapopulations

Using demes from experimental metapopulations of the flour beetle, Tribolium castaneum, we investigated phase 3 of Wright's shifting balance process. Using parent demes of high, intermediate, and low mean fitness, we experimentally modeled migration of varying amounts from demes of high mean fitness into demes of lower mean fitness (like phase 3) as well as the reciprocal (the opposite of phase 3). In natural populations, some migration among demes occurs independently of deme fitness. In this case, demes of high mean fitness are likely to receive migrants from demes of lower mean fitness; these effects might limit the effectiveness of phase 3 but have not been studied experimentally. We estimated the populational heritability of mean fitness by the regression of offspring deme means on the weighted parental means and found moderate levels of demic heritability one (0.641-0.690) and two (0.518-0.552) generations after migration. We discuss our findings in relation to the role of interdemic migration in "adaptive peak shifts" in metapopulations and the controversies over group selection and the units of inheritance.     

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.      

Intra-deme molecular diversity in spatially expanding populations

We report here a simulation study examining the effect of a recent spatial expansion on the pattern of molecular diversity within a deme. We first simulate a range expansion in a virtual world consisting in a two-dimensional array of demes exchanging a given proportion of migrants (m) with their neighbors. The recorded demographic and migration histories are then used under a coalescent approach to generate the genetic diversity in a sample of genes. We find that the shape of the gene genealogies and the overall pattern of diversity within demes depend not only on the age of the expansion but also on the level of gene flow between neighboring demes, as measured by the product Nm, where N is the size of a deme. For small Nm values (< approximately 20 migrants sent outwards per generation), a substantial proportion of coalescent events occur early in the genealogy, whereas with larger levels of gene flow, most coalescent events occur around the time of the onset of the spatial expansion. Gene genealogies are star shaped, and mismatch distributions are unimodal after a range expansion for large Nm values. In contrast, gene genealogies present a mixture of both very short and very long branch lengths, and mismatch distributions are multimodal for small Nm values. It follows that statistics used in tests of selective neutrality like Tajima's D statistic or Fu's F(S) statistic will show very significant negative values after a spatial expansion only in demes with high Nm values. In the context of human evolution, this difference could explain very simply the fact that analyses of samples of mitochondrial DNA sequences reveal multimodal mismatch distributions in hunter-gatherers and unimodal distributions in post-Neolithic populations. Indeed, the current simulations show that a recent increase in deme size (resulting in a larger Nm value) is sufficient to prevent recent coalescent events and thus lead to unimodal mismatch distributions, even if deme sizes (and therefore Nm values) were previously much smaller. The fact that molecular diversity within deme is so dependent on recent levels of gene flow suggests that it should be possible to estimate Nm values from samples drawn from a single deme.     

Adaptive genetic structure in phytophagous insect populations

Genetic variation in insect populations is frequently structured into discrete groups, or demes, that form in response to stochastic forces or natural selection. Because host-plant populations are often highly heterogeneous, phytophagous insects may form demes that are adapted to the unique traits of individual plants. Recent field experiments indicate that selection pressures imposed by host-plants can promote rapid adaptive evolution in natural insect populations at very fine spatial scales. Adaptive deme formation may be more common among endophagous insects, which feed and reside within plant tissue, than for externally feeding insects, because internal feeders experience stronger plant-mediated selection pressures.     

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.     

The coalescent in an island model of population subdivision with variation among demes

A simple genealogical structure is found for a general finite island model of population subdivision. The model allows for variation in the sizes of demes, in contributions to the migrant pool, and in the fraction of each deme that is replaced by migrants every generation. The ancestry of a sample of non-recombining DNA sequences has a simple structure when the sample size is much smaller than the total number of demes in the population. This allows an expression for the probability distribution of the number of segregating sites in the sample to be derived under the infinite-sites mutation model. It also yields easily computed estimators of the migration parameter for each deme in a multi-deme sample. The genealogical process is such that the lineages ancestral to the sample tend to accumulate in demes with low migration rates and/or which contribute disproportionately to the migrant pool. In addition, common ancestor or coalescent events tend to occur in demes of small size. This provides a framework for understanding the determinants of the effective size of the population, and leads to an expression for the probability that the root of a genealogy occurs in a particular geographic region, or among a particular set of demes.     

Dominance and the maintenance of polymorphism in multiallelic migration-selection models with two demes

The maintenance of genetic variation in a spatially heterogeneous environment has been one of the main research themes in theoretical population genetics. Despite considerable progress in understanding the consequences of spatially structured environments on genetic variation, many problems remain unsolved. One of them concerns the relationship between the number of demes, the degree of dominance, and the maximum number of alleles that can be maintained by selection in a subdivided population. In this work, we study the potential of maintaining genetic variation in a two-deme model with deme-independent degree of intermediate dominance, which includes absence of G×E interaction as a special case. We present a thorough numerical analysis of a two-deme three-allele model, which allows us to identify dominance and selection patterns that harbor the potential for stable triallelic equilibria. The information gained by this approach is then used to construct an example in which existence and asymptotic stability of a fully polymorphic equilibrium can be proved analytically. Noteworthy, in this example the parameter range in which three alleles can coexist is maximized for intermediate migration rates. Our results can be interpreted in a specialist-generalist context and (among others) show when two specialists can coexist with a generalist in two demes if the degree of dominance is deme independent and intermediate. The dominance relation between the generalist allele and the specialist alleles play a decisive role. We also discuss linear selection on a quantitative trait and show that G×E interaction is not necessary for the maintenance of more than two alleles in two demes.     

Spatial heterogeneity in the strength of selection against deleterious alleles and the mutation load

According to current estimates of genomic deleterious mutation rates (which are often of the order 0.1-1) the mutation load (defined as a reduction in the average fitness of a population due to the presence of deleterious alleles) may be important in many populations. In this paper, I use multilocus simulations to explore the effect of spatial heterogeneity in the strength of selection against deleterious alleles on the mutation load (for example, it has been suggested that stressful environments may increase the strength of selection). These simulations show contrasted results: in some situations, spatial heterogeneity may greatly reduce the mutation load, due to the fact that migrants coming from demes under stronger selection carry relatively few deleterious alleles, and benefit from a strong advantage within demes under weaker selection (where individuals carry many more deleterious alleles); in other situations, however, deleterious alleles accumulate within demes under stronger selection, due to migration pressure from demes under weaker selection, leading to fitness erosion within those demes. This second situation is more frequent when the productivity of the different demes is proportional to their mean fitness. The effect of spatial heterogeneity is greatly reduced, however, when the response to environmental differences is inconsistent across loci.     

F(st) in a Hierarchical Island Model

It is shown that in a hierarchical island model, in which demes within a neighborhood exchange migrants at a much higher rate than do demes in different neighborhoods, hierarchical F statistics introduced by S. Wright can indicate the extent of gene flow within and between neighborhoods. At equilibrium, the within-neighborhood inbreeding coefficient, FSN, is approximately 1/(1 + 4Nm1) where N is the deme size and m1 is the migration rate among demes in the same neighborhood. The between-neighborhood inbreeding coefficient, FNT, is approximately 1/(1 + 4Ndm2) where d is the number of demes in a neighborhood and m2 is the migration rate among demes in different neighborhoods.     

Environmental Noise,Genetic Diversity and the Evolution of Evolvability and Robustness in Model Gene Networks

The ability of organisms to adapt and persist in the face of environmental change is accepted as a fundamental feature of natural systems. More contentious is whether the capacity of organisms to adapt (or “evolvability”) can itself evolve and the mechanisms underlying such responses. Using model gene networks, I provide evidence that evolvability emerges more readily when populations experience positively autocorrelated environmental noise (red noise) compared to populations in stable or randomly varying (white noise) environments. Evolvability was correlated with increasing genetic robustness to effects on network viability and decreasing robustness to effects on phenotypic expression; populations whose networks displayed greater viability robustness and lower phenotypic robustness produced more additive genetic variation and adapted more rapidly in novel environments. Patterns of selection for robustness varied antagonistically with epistatic effects of mutations on viability and phenotypic expression, suggesting that trade-offs between these properties may constrain their evolutionary responses. Evolution of evolvability and robustness was stronger in sexual populations compared to asexual populations indicating that enhanced genetic variation under fluctuating selection combined with recombination load is a primary driver of the emergence of evolvability. These results provide insight into the mechanisms potentially underlying rapid adaptation as well as the environmental conditions that drive the evolution of genetic interactions.     

Selection for recombination in structured populations

In finite populations, linkage disequilibria generated by the interaction of drift and directional selection (Hill-Robertson effect) can select for sex and recombination, even in the absence of epistasis. Previous models of this process predict very little advantage to recombination in large panmictic populations. In this article we demonstrate that substantial levels of linkage disequilibria can accumulate by drift in the presence of selection in populations of any size, provided that the population is subdivided. We quantify (i) the linkage disequilibrium produced by the interaction of drift and selection during the selective sweep of beneficial alleles at two loci in a subdivided population and (ii) the selection for recombination generated by these disequilibria. We show that, in a population subdivided into n demes of large size N, both the disequilibrium and the selection for recombination are equivalent to that expected in a single population of a size intermediate between the size of each deme (N) and the total size (nN), depending on the rate of migration among demes, m. We also show by simulations that, with small demes, the selection for recombination is stronger than both that expected in an unstructured population (m = 1 - 1/n) and that expected in a set of isolated demes (m = 0). Indeed, migration maintains polymorphisms that would otherwise be lost rapidly from small demes, while population structure maintains enough local stochasticity to generate linkage disequilibria. These effects are also strong enough to overcome the twofold cost of sex under strong selection when sex is initially rare. Overall, our results show that the stochastic theories of the evolution of sex apply to a much broader range of conditions than previously expected.     

Invasion and fixation of sex-reversal genes

We simulated a meta-population with random dispersal among demes but local mating within demes to investigate conditions under which a dominant female-determining gene W, with no individual selection advantage, can invade and become fixed in females, changing the population from male to female heterogamety. Starting with one mutant W in a single deme, the interaction of sex ratio selection and random genetic drift causes W to be fixed among females more often than a comparable neutral mutation with no influence on sex determination, even when YY males have slightly reduced viability. Meta-population structure and interdeme selection can also favour the fixation of W. The reverse transition from female to male heterogamety can also occur with higher probability than for a comparable neutral mutation. These results help to explain the involvement of sex-determining genes in the evolution of sex chromosomes and in sexual selection and speciation.     

EVOLUTION OF ALTRUISM IN KIN-STRUCTURED AND RANDOM SUBDIVIDED POPULATIONS

A. population structure favorable to the evolution of an altruistic trait is studied by Monte Carlo simulation. The model is based on a small-scale nonindustrial human society but seems generalizable to other highly social mammals. Three hierarchical levels are recognized: 1) the ecologically isolated local group (hamlet) which may be composed of kin and/or unrelated individuals; 2) the deme (settlement) comprising several such groups which interbreed; and 3) the set of demes (metapopulation) among which gene flow occurs. The first two levels of the model are based on D. S. Wilson's structured deme concept; the third allows for gene flow among demes in the metapopulation and for the structured diffusion of alleles across a wider area than might be included within the scope of a single deme. The simulation models genetic drift by a process of hamlet formation which may be random, or variously kin-structured. Hamlets may then become extinct based on a probability function of their gene frequencies. Individual selection within settlements is modeled deterministically, and gene flow among settlements is modeled as two-dimensional steppingstone migration of random or kin-structured groups. Results of the simulations show that, with realistic values for group sizes, moderate extinction rate, and high rates of migration (m > 27%), disadvantageous alleles (s = 10% and 25%) may increase markedly due to differential hamlet extinction over the course of 50 generations. The greater the degree of kin-structuring of founder groups, the higher the variance among hamlets and the faster the rate of increase of the allele for altruism. Nonetheless, even in some randomly founded groups, a clear increase in the altruism gene frequency occurred. It is also notable that kin-structured group selection by hamlet extinction may be effective when the initial frequency of altruism genes is very low (average of one per deme) and among a relatively small number of demes (25). Thus the process of group extinction in a hierarchically structured population allows rapid increase of an allele for altruism under plausible demographic conditions.     

Evolution of evolvability in a developmental model

Evolvability, the ability of populations to adapt, can evolve through changes in the mechanisms determining genetic variation and in the processes of development. Here we construct and evolve a simple developmental model in which the pleiotropic effects of genes can evolve. We demonstrate that selection in a changing environment favors a specific pattern of variability, and that this favored pattern maximizes evolvability. Our analysis shows that mutant genotypes with higher evolvability are more likely to increase to fixation. We also show that populations of highly evolvable genotypes are much less likely to be invaded by mutants with lower evolvability, and that this dynamic primarily shapes evolvability. We examine several theoretical objections to the evolution of evolvability in light of this result. We also show that this result is robust to the presence or absence of recombination, and explore how nonrandom environmental change can select for a modular pattern of variability.     

Genetic load,inbreeding depression and heterosis in an age‐structured metapopulation

In a metapopulation, the process of recurrent local extinction and recolonization gives rise to an age structure among demes. Recently established demes will tend to differ from older demes in terms of the levels of genetic diversity found within them and the way this diversity is distributed among demes in the same and different ages. The effects of population turnover on average levels of genetic diversity among demes in a metapopulation have been the focus of much attention, both for neutral and nonneutral loci, but much less is known about the distribution of nonneutral genetic diversity among demes of different ages. In this paper, we used computer simulations to study the distribution of genetic load, inbreeding depression and heterosis in an age‐structured metapopulation. We found that, for mildly deleterious mutations, within‐deme inbreeding depression increased, whereas heterosis and genetic load decreased with deme age following severe colonization bottlenecks. In contrast, recessive lethal alleles tended to be purged during colonization, with older populations showing higher genetic load and higher within‐deme inbreeding depression. Heterosis caused by recessive lethal alleles and resulting from gene flow among different demes tended to be greatest for young demes, because the mutations responsible tended to be purged in the first few generations after colonization, but its effects increased again as populations grow older as a result of immigration. Our results point to a need for estimates of genetic diversity, genetic load, within‐deme inbreeding depression and heterosis in demes of different age classes separately.     

The effects of stress and sex on selection,genetic covariance,and the evolutionary response

The capacity of a population to adapt to selection (evolvability) depends on whether the structure of genetic variation permits the evolution of fitter trait combinations. Selection, genetic variance and genetic covariance can change under environmental stress, and males and females are not genetically independent, yet the combined effects of stress and dioecy on evolvability are not well understood. Here, we estimate selection, genetic (co)variance and evolvability in both sexes of Tribolium castaneum flour beetles under stressful and benign conditions, using a half‐sib breeding design. Although stress uncovered substantial latent heritability, stress also affected genetic covariance, such that evolvability remained low under stress. Sexual selection on males and natural selection on females favoured a similar phenotype, and there was positive intersex genetic covariance. Consequently, sexual selection on males augmented adaptation in females, and intralocus sexual conflict was weak or absent. This study highlights that increased heritability does not necessarily increase evolvability, suggests that selection can deplete genetic variance for multivariate trait combinations with strong effects on fitness, and tests the recent hypothesis that sexual conflict is weaker in stressful or novel environments.     

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.     

THE PROBABILITY OF FIXATION OF A NEW KARYOTYPE IN A CONTINUOUS POPULATION

We investigate the probability of fixation of a chromosome rearrangement in a subdivided population, concentrating on the limit where migration is so large relative to selection (m ? s) that the population can be thought of as being continuously distributed. We study two demes, and one- and two-dimensional populations. For two demes, the probability of fixation in the limit of high migration approximates that of a population with twice the size of a single deme: migration therefore greatly reduces the fixation probability. However, this behavior does not extend to a large array of demes. Then, the fixation probability depends primarily on neighborhood size (Nb), and may be appreciable even with strong selection and free gene flow (≈exp(-B ≈ Nb√s) in one dimension, ≈exp(-B ≈ Nb) in two dimensions). Our results are close to those for the more tractable case of a polygenic character under disruptive selection.