The effect of population bottlenecks on mutation rate evolution in asexual populations

In the absence of recombination, a mutator allele can spread through a population by hitchhiking with beneficial mutations that appear in its genetic background. Theoretical studies over the past decade have shown that the survival and fixation probability of beneficial mutations can be severely reduced by population size bottlenecks. Here, we use computational modelling and evolution experiments with the yeast S. cerevisiae to examine whether population bottlenecks can affect mutator dynamics in adapting asexual populations. In simulation, we show that population bottlenecks can inhibit mutator hitchhiking with beneficial mutations and are most effective at lower beneficial mutation supply rates. We then subjected experimental populations of yeast propagated at the same effective population size to three different bottleneck regimes and observed that the speed of mutator hitchhiking was significantly slower at smaller bottlenecks, consistent with our theoretical expectations. Our results, thus, suggest that bottlenecks can be an important factor in mutation rate evolution and can in certain circumstances act to stabilize or, at least, delay the progressive elevation of mutation rates in asexual populations. Additionally, our findings provide the first experimental support for the theoretically postulated effect of population bottlenecks on beneficial mutations and demonstrate the usefulness of studying mutator frequency dynamics for understanding the underlying dynamics of fitness‐affecting mutations.     

Competition between high- and higher-mutating strains of Escherichia coli

Experimental studies have shown that a mutator allele can readily hitchhike to fixation with beneficial mutations in an asexual population having a low, wild-type mutation rate. Here, we show that a genotype bearing two mutator alleles can supplant a population already fixed for one mutator allele. Our results provide experimental support for recent theory predicting that mutator alleles will tend to accumulate in asexual populations by hitchhiking with beneficial mutations, causing an ever-higher genomic mutation rate.     

Limits to adaptation in asexual populations

In asexual populations, the rate of adaptation is basically limited by the frequency and properties of spontaneous beneficial mutations. Hence, knowledge of these mutational properties and how they are affected by particular evolutionary conditions is a precondition for understanding the process of adaptation. Here, we address how the rate of adaptation of asexual populations is limited by its population size and mutation rate, as well as by two factors affecting the fraction of mutations that confer a benefit, i.e. the initial adaptedness of the population and the variability of the environment. These factors both influence which mutations are likely to occur, as well as the probability that they will ultimately contribute to adaptation. We attempt to separate the consequences of these basic population features in terms of their effect on the rate of adaptation by using results from evolution experiments with microorganisms.     

Adaptive evolution of asexual populations under Muller's ratchet

We study the population genetics of adaptation in nonequilibrium haploid asexual populations. We find that the accumulation of deleterious mutations, due to the operation of Muller's ratchet, can considerably reduce the rate of fixation of advantageous alleles. Such reduction can be approximated reasonably well by a reduction in the effective population size. In the absence of Muller's ratchet, a beneficial mutation can only become fixed if it creates the best possible genotype; if Muller's ratchet operates, however, mutations initially arising in a nonoptimal genotype can also become fixed in the population, since the loss of the least-loaded class implies that an initially nonoptimal background can become optimal. We show that, while the rate at which adaptive mutations become fixed is reduced, the rate of fixation of deleterious mutations due to the ratchet is not changed by the presence of beneficial mutations as long as the rate of their occurrence is low and the deleterious effects of mutations (s(d)) are higher than the beneficial effects (s(a)). When s(a) > s(d), the advantage of a beneficial mutation can outweigh the deleterious effects of associated mutations. Under these conditions, a beneficial allele can drag to fixation deleterious mutations initially associated with it at a higher rate than in the absence of advantageous alleles. We propose analytical approximations for the rates of accumulation of deleterious and beneficial mutations. Furthermore, when allowing for the possible occurrence of interference between beneficial alleles, we find that the presence of deleterious mutations of either very weak or very strong effect can marginally increase the rate of accumulation of beneficial mutations over that observed in the absence of such deleterious mutations.     

Variance effective population size is affected by census size in sub-structured populations

Measurement of allele frequency shifts between temporally spaced samples has long been used for assessment of effective population size (N_e), and this ‘temporal method’ provides estimates of N_e referred to as variance effective size (N_eV). We show that N_eV of a local population that belongs to a sub-structured population (a metapopulation) is determined not only by genetic drift and migration rate (m), but also by the census size (N_c). The realized N_eV of a local population can either increase or decrease with increasing m, depending on the relationship between N_e and N_c in isolation. This is shown by explicit mathematical expressions for the factors affecting N_eV derived for an island model of migration. We verify analytical results using high-resolution computer simulations, and show that the phenomenon is not restricted to the island model migration pattern. The effect of N_c on the realized N_eV of a local subpopulation is most pronounced at high migration rates. We show that N_c only affects local N_eV, whereas N_eV for the metapopulation as a whole, inbreeding (N_eI), and linkage disequilibrium (N_eLD) effective size are all independent of N_c. Our results provide a possible explanation to the large variation of N_e/N_c ratios reported in the literature, where N_e is frequently estimated by N_eV. They are also important for the interpretation of empirical N_e estimates in genetic management where local N_eV is often used as a substitute for inbreeding effective size, and we suggest an increased focus on metapopulation N_eV as a proxy for N_eI.     

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.     

Samples from subdivided populations yield biased estimates of effective size that overestimate the rate of loss of genetic variation

Many empirical studies estimating effective population size apply the temporal method that provides an estimate of the variance effective size through the amount of temporal allele frequency change under the assumption that the study population is completely isolated. This assumption is frequently violated, and the magnitude of the resulting bias is generally unknown. We studied how gene flow affects estimates of effective size obtained by the temporal method when sampling from a population system and provide analytical expressions for the expected estimate under an island model of migration. We show that the temporal method tends to systematically underestimate both local and global effective size when populations are connected by gene flow, and the bias is sometimes dramatic. The problem is particularly likely to occur when sampling from a subdivided population where high levels of gene flow obscure identification of subpopulation boundaries. In such situations, sampling in a manner that prevents biased estimates can be difficult. This phenomenon might partially explain the frequently reported unexpectedly low effective population sizes of marine populations that have raised concern regarding the genetic vulnerability of even exceptionally large populations.     

Allelic and genotypic diversity in long-term asexual populations of the pea aphid, Acyrthosiphon pisum in comparison with sexual populations

Many aphid species exhibit geographical variation in the mode of reproduction that ranges from cyclical parthenogenesis with a sexual phase to obligate parthenogenesis (asexual reproduction). Theoretical studies predict that organisms reproducing asexually should maintain higher allelic diversity per locus but lower genotypic diversity than organisms reproducing sexually. To corroborate this hypothesis, we evaluated genotypic and allelic diversities in the sexual and asexual populations of the pea aphid, Acyrthosiphon pisum (Harris). Microsatellite analysis revealed that populations in central Japan are asexual, whereas populations in northern Japan are obligatorily sexual. No mixed populations were detected in our study sites. Phylogenetic analysis using microsatellite data and mitochondrial cytochrome oxidase subunit I (COI) gene sequences revealed a long history of asexuality in central Japan and negated the possibility of the recent origin of the asexual populations from the sexual populations. Asexual populations exhibited much lower genotypic diversity but higher allelic richness per locus than did sexual populations. Asexual populations consisted of a few predominant clones that were considerably differentiated from one another. Sexual populations on alfalfa, an exotic plant in Japan, were most closely related to asexual populations associated with Vicia sativa L. The alfalfa-associated sexual populations harboured one COI haplotype that was included in the haplotype clade of the asexual populations. Available evidence suggests that the sexuality of the alfalfa-associated populations has recently been restored through the northward migration and colonization of alfalfa by V. sativa- associated lineages. Therefore, our results support the theoretical predictions and provide a new perspective on the origin of sexual populations.     

Adaptation,Clonal Interference,and Frequency-Dependent Interactions in a Long-Term Evolution Experiment with Escherichia coli

Twelve replicate populations of Escherichia coli have been evolving in the laboratory for >25 years and 60,000 generations. We analyzed bacteria from whole-population samples frozen every 500 generations through 20,000 generations for one well-studied population, called Ara−1. By tracking 42 known mutations in these samples, we reconstructed the history of this population’s genotypic evolution over this period. The evolutionary dynamics of Ara−1 show strong evidence of selective sweeps as well as clonal interference between competing lineages bearing different beneficial mutations. In some cases, sets of several mutations approached fixation simultaneously, often conveying no information about their order of origination; we present several possible explanations for the existence of these mutational cohorts. Against a backdrop of rapid selective sweeps both earlier and later, two genetically diverged clades coexisted for >6000 generations before one went extinct. In that time, many additional mutations arose in the clade that eventually prevailed. We show that the clades evolved a frequency-dependent interaction, which prevented the immediate competitive exclusion of either clade, but which collapsed as beneficial mutations accumulated in the clade that prevailed. Clonal interference and frequency dependence can occur even in the simplest microbial populations. Furthermore, frequency dependence may generate dynamics that extend the period of coexistence that would otherwise be sustained by clonal interference alone.     

Clonal interference is alleviated by high mutation rates in large populations

When a beneficial mutation is fixed in a population that lacks recombination, the genetic background linked to that mutation is fixed. As a result, beneficial mutations on different backgrounds experience competition, or "clonal interference," that can cause asexual populations to evolve more slowly than their sexual counterparts. Factors such as a large population size (N) and high mutation rates (mu) increase the number of competing beneficial mutations, and hence are expected to increase the intensity of clonal interference. However, recent theory suggests that, with very large values of Nmu, the severity of clonal interference may instead decline. The reason is that, with large Nmu, genomes including both beneficial mutations are rapidly created by recurrent mutation, obviating the need for recombination. Here, we analyze data from experimentally evolved asexual populations of a bacteriophage and find that, in these nonrecombining populations with very large Nmu, recurrent mutation does appear to ameliorate this cost of asexuality.     

The role of migration in mutant dynamics in fragmented populations

Mutant dynamics in fragmented populations have been studied extensively in evolutionary biology. Yet, open questions remain, both experimentally and theoretically. Some of the fundamental properties predicted by models still need to be addressed experimentally. We contribute to this by using a combination of experiments and theory to investigate the role of migration in mutant distribution. In the case of neutral mutants, while the mean frequency of mutants is not influenced by migration, the probability distribution is. To address this empirically, we performed in vitro experiments, where mixtures of GFP-labelled (“mutant”) and non-labelled (“wid-type”) murine cells were grown in wells (demes), and migration was mimicked via cell transfer from well to well. In the presence of migration, we observed a change in the skewedness of the distribution of the mutant frequencies in the wells, consistent with previous and our own model predictions. In the presence of de novo mutant production, we used modelling to investigate the level at which disadvantageous mutants are predicted to exist, which has implications for the adaptive potential of the population in case of an environmental change. In panmictic populations, disadvantageous mutants can persist around a steady state, determined by the rate of mutant production and the selective disadvantage (selection-mutation balance). In a fragmented system that consists of demes connected by migration, a steady-state persistence of disadvantageous mutants is also observed, which, however, is fundamentally different from the mutation-selection balance and characterized by higher mutant levels. The increase in mutant frequencies above the selection-mutation balance can be maintained in small ($N<N_c$) demes as long as the migration rate is sufficiently small. The migration rate above which the mutants approach the selection-mutation balance decays exponentially with $N/N_c$. The observed increase in the mutant numbers is not explained by the change in the effective population size. Implications for evolutionary processes in diseases are discussed, where the pre-existence of disadvantageous drug-resistant mutant cells or pathogens drives the response of the disease to treatments.     

Effective population size,genetic diversity,and coalescence time in subdivided populations

A formula for the effective population size for the finite island model of subdivided populations is derived. The formula indicates that the effective size can be substantially greater than the actual number of individuals in the entire population when the migration rate among subpopulations is small. It is shown that the mean nucleotide diversity, coalescence time, and heterozygosity for genes sampled from the entire population can be predicted fairly well from the theory for randomly mating populations if the effective population size for the finite island model is used.     

Effective size of density‐dependent populations in fluctuating environments

Reliable estimates of effective population size are of central importance in population genetics and evolutionary biology. For populations that fluctuate in size, harmonic mean population size is commonly used as a proxy for (multi‐) generational effective size. This assumes no effects of density dependence on the ratio between effective and actual population size, which limits its potential application. Here, we introduce density dependence on vital rates in a demographic model of variance effective size. We derive an expression for the ratio in a density‐regulated population in a fluctuating environment. We show by simulations that yearly genetic drift is accurately predicted by our model, and not proportional to as assumed by the harmonic mean model, where N is the total population size of mature individuals. We find a negative relationship between and N. For a given N, the ratio depends on variance in reproductive success and the degree of resource limitation acting on the population growth rate. Finally, our model indicate that environmental stochasticity may affect not only through fluctuations in N, but also for a given N at a given time. Our results show that estimates of effective population size must include effects of density dependence and environmental stochasticity.     

Contrasting effects of large density changes on relative testes size in fluctuating populations of sympatric vole species

Across species, there is usually a positive relationship between sperm competition level and male reproductive effort on ejaculates, typically measured using relative testes size (RTS). Within populations, demographic and ecological processes may drastically alter the level of sperm competition and thus, potentially affect the evolution of testes size. Here, we use longitudinal records (across 38 years) from wild sympatric Fennoscandian populations of five species of voles to investigate whether RTS responds to natural fluctuations in population density, i.e. variation in sperm competition risk. We show that for some species RTS increases with density. However, our results also show that this relationship can be reversed in populations with large-scale between-year differences in density. Multiple mechanisms are suggested to explain the negative RTS–density relationship, including testes size response to density-dependent species interactions, an evolutionary response to sperm competition levels that is lagged when density fluctuations are over a certain threshold, or differing investment in pre- and post-copulatory competition at different densities. The results emphasize that our understanding of sperm competition in fluctuating environments is still very limited.     

Three components of genetic drift in subdivided populations

Wright's metaphor of sampling is extended to consider three components of genetic drift: those occurring before, during, and after migration. To the extent that drift at each stage behaves like an independent random sample, the order of events does not matter. When sampling is not random, the order does matter, and the effect of population size is confounded with that of mobility. The widely cited result that genetic differentiation of local groups depends only on the product of group size and migration rate holds only when nonrandom sampling does not occur prior to migration in the life cycle.     

Adaptive evolution in a spatially structured asexual population

We study the process of adaptation in a spatially structured asexual haploid population. The model assumes a local competition for replication, where each organism interacts only with its nearest neighbors. We observe that the substitution rate of beneficial mutations is smaller for a spatially structured population than that seen for populations without structure. The difference between structured and unstructured populations increases as the adaptive mutation rate increases. Furthermore, the substitution rate decreases as the number of neighbors for local competition is reduced. We have also studied the impact of structure on the distribution of adaptive mutations that fix during adaptation.     

Drift load in populations of small size and low density

According to theory, drift load in randomly mating populations is determined by past population size, because enhanced genetic drift in small populations causes accumulation and fixation of recessive deleterious mutations of small effect. In contrast, segregating load due to mutations of low frequency should decline in smaller populations, at least when mutations are highly recessive and strongly deleterious. Strong local selection generally reduces both types of load. We tested these predictions in 13 isolated, outcrossing populations of Arabidopsis lyrata that varied in population size and plant density. Long-term size was estimated by expected heterozygosity at 20 microsatellite loci. Segregating load was assessed by comparing performance of offspring from selfings versus within-population crosses. Drift load was the heterosis effect created by interpopulation outbreeding. Results showed that segregating load was unrelated to long-term size. However, drift load was significantly higher in populations of small effective size and low density. Drift load was mostly expressed late in development, but started as early as germination and accumulated thereafter. The study largely confirms predictions of theory and illustrates that mutation accumulation can be a threat to natural populations.     

Parthenogenetic female populations in the brown alga Scytosiphon lomentaria (Scytosiphonaceae,Ectocarpales): decay of a sexual trait and acquisition of asexual traits

In isogamous brown algae, the sexuality of populations needs to be tested by laboratory crossing experiments, as the sexes of gametophytes are morphologically indistinguishable. In some cases, gamete fusion is not observed and the precise reproductive mode of the populations is unknown. In the isogamous brown alga Scytosiphon lomentaria in Japan, both asexual (gamete fusion is unobservable) and sexual populations (gamete fusion is observable) have been reported. In order to elucidate the reproductive mode of asexual populations in this species, we used PCR‐based sex markers to investigate the sex ratio of three asexual and two sexual field populations. The markers indicated that the asexual populations consisted only of female individuals, whereas sexual populations are composed of both males and females. In culture, female gametes of most strains from asexual populations were able to fuse with male gametes; however, they had little to no detectable sexual pheromones, significantly larger cell sizes, and more rapid parthenogenetic development compared to female/male gametes from sexual populations. Investigations of sporophytic stages in the field indicated that alternation of gametophytic and parthenosporophytic stages occur in an asexual population. These results indicate that the S. lomentaria asexual populations are female populations that lack sexual reproduction and reproduce parthenogenetically. It is likely that females in the asexual populations have reduced a sexual trait (pheromone production) and have acquired asexual traits (larger gamete sizes and rapid parthenogenetic development).     

The effect of spatial structure on adaptation in Escherichia coli

Populations of organisms are generally organized in a given spatial structure. However, the vast majority of population genetic studies are based on populations in which every individual competes globally. Here we use experimental evolution in Escherichia coli to directly test a recently made prediction that spatial structure slows down adaptation and that this cost increases with the mutation rate. This was studied by comparing populations of different mutation rates adapting to a liquid (unstructured) medium with populations that evolved in a Petri dish on solid (structured) medium. We find that mutators adapt faster to both environments and that adaptation is slower if there is spatial structure. We observed no significant difference in the cost of structure between mutator and wild-type populations, which suggests that clonal interference is intense in both genetic backgrounds.