Fitness-dependent mutation rates in finite populations

Replicator-mutator dynamics of the rock-paper-scissors game: Learning through mistakes

We generalize the Bush-Mosteller learning, the Roth-Erev learning, and the social learning to include mistakes, such that the nonlinear replicator-mutator equation with either additive or multiplicative mutation is generated in an asymptotic limit. Subsequently, we exhaustively investigate the ubiquitous rock-paper-scissors game for some analytically tractable motifs of mutation pattern for which the replicator-mutator flow is seen to exhibit rich dynamics that include limit cycles and chaotic orbits. The main result of this paper is that in both symmetric and asymmetric game interactions, mistakes can sometimes help the players learn; in fact, mistakes can even control chaos to lead to rational Nash-equilibrium outcomes. Furthermore, we report a hitherto-unknown Hamiltonian structure of the replicator-mutator equation.

Global coordination of the mutation and growth rates across the genetic and nutritional variety in Escherichia coli

Fitness and mutability are the primary traits of living organisms for adaptation and evolution. However, their quantitative linkage remained largely deficient. Whether there is any general relationship between the two features and how genetic and environmental variables influence them remained unclear and were addressed here. The mutation and growth rates of an assortment of Escherichia coli strain collections, including the wild-type strains and the genetically disturbed strains of either reduced genomes or deletion of the genes involved in the DNA replication fidelity, were evaluated in various media. The contribution of media to the mutation and growth rates was differentiated depending on the types of genetic disturbance. Nevertheless, the negative correlation between the mutation and growth rates was observed across the genotypes and was common in all media. It indicated the comprehensive association of the correlated mutation and growth rates with the genetic and medium variation. Multiple linear regression and support vector machine successfully predicted the mutation and growth rates and the categories of genotypes and media, respectively. Taken together, the study provided a quantitative dataset linking the mutation and growth rates, genotype, and medium and presented a simple and successful example of predicting bacterial growth and mutability by data-driven approaches.

Mutations Beget More Mutations-Rapid Evolution of Mutation Rate in Response to the Risk of Runaway Accumulation

The rapidity with which the mutation rate evolves could greatly impact evolutionary patterns. Nevertheless, most studies simply assume a constant rate in the time scale of interest (Kimura 1983; Drake 1991; Kumar 2005; Li 2007; Lynch 2010). In contrast, recent studies of somatic mutations suggest that the mutation rate may vary by several orders of magnitude within a lifetime (Kandoth et al. 2013; Lawrence et al. 2013). To resolve the discrepancy, we now propose a runaway model, applicable to both the germline and soma, whereby mutator mutations form a positive-feedback loop. In this loop, any mutator mutation would increase the rate of acquiring the next mutator, thus triggering a runaway escalation in mutation rate. The process can be initiated more readily if there are many weak mutators than a few strong ones. Interestingly, even a small increase in the mutation rate at birth could trigger the runaway process, resulting in unfit progeny. In slowly reproducing species, the need to minimize the risk of this uncontrolled accumulation would thus favor setting the mutation rate low. In comparison, species that starts and ends reproduction sooner do not face the risk and may set the baseline mutation rate higher. The mutation rate would evolve in response to the risk of runaway mutation, in particular, when the generation time changes. A rapidly evolving mutation rate may shed new lights on many evolutionary phenomena (Elango et al. 2006; Thomas et al. 2010, 2018; Langergraber et al. 2012; Besenbacher et al. 2019).

Evolution of Stress-Induced Mutagenesis in the Presence of Horizontal Gene Transfer

Stress-induced mutagenesis has been observed in multiple species of bacteria and yeast. It has been suggested that in asexual populations, a mutator allele that increases the mutation rate during stress can sweep to fixation with the beneficial mutations it generates. However, even asexual microbes can undergo horizontal gene transfer and rare recombination, which typically interfere with the spread of mutator alleles. Here we examine the effect of horizontal gene transfer on the evolutionary advantage of stress-induced mutator alleles. Our results demonstrate that stress-induced mutator alleles are favored by selection even in the presence of horizontal gene transfer and more so when the mutator alleles also increase the rate of horizontal gene transfer. We suggest that when regulated by stress, mutation and horizontal gene transfer can be complementary rather than competing adaptive strategies and that stress-induced mutagenesis has important implications for evolutionary biology, ecology, and epidemiology, even in the presence of horizontal gene transfer and rare recombination.

Desiccation-induced changes in recombination rate and crossover interference in Drosophila melanogaster: evidence for fitness-dependent plasticity

Meiotic recombination is evolutionarily ambiguous, as being associated with both benefits and costs to its bearers, with the resultant dependent on a variety of conditions. While existing theoretical models explain the emergence and maintenance of recombination, some of its essential features remain underexplored. Here we focus on one such feature, recombination plasticity, and test whether recombination response to stress is fitness-dependent. We compare desiccation stress effects on recombination rate and crossover interference in chromosome 3 between desiccation-sensitive and desiccation-tolerant Drosophila lines. We show that relative to desiccation-tolerant genotypes, desiccation-sensitive genotypes exhibit a significant segment-specific increase in single- and double-crossover frequencies across the pericentromeric region of chromosome 3. Significant changes (relaxation) in crossover interference were found for the interval pairs flanking the centromere and extending to the left arm of the chromosome. These results indicate that desiccation is a recombinogenic factor and that desiccation-induced changes in both recombination rate and crossover interference are fitness-dependent, with a tendency of less fitted individuals to produce more variable progeny. Such dependence may play an important role in the regulation of genetic variation in populations experiencing environmental challenges.

Evolution of the Mutational Process under Relaxed Selection in Caenorhabditis elegans

The mutational process varies at many levels, from within genomes to among taxa. Many mechanisms have been linked to variation in mutation, but understanding of the evolution of the mutational process is rudimentary. Physiological condition is often implicated as a source of variation in microbial mutation rate and may contribute to mutation rate variation in multicellular organisms.Deleterious mutations are an ubiquitous source of variation in condition. We test the hypothesis that the mutational process depends on the underlying mutation load in two groups of Caenorhabditis elegans mutation accumulation (MA) lines that differ in their starting mutation loads. "First-order MA" (O1MA) lines maintained under minimal selection for ∼250 generations were divided into high-fitness and low-fitness groups and sets of "second-order MA" (O2MA) lines derived from each O1MA line were maintained for ∼150 additional generations. Genomes of 48 O2MA lines and their progenitors were sequenced. There is significant variation among O2MA lines in base-substitution rate (µbs), but no effect of initial fitness; the indel rate is greater in high-fitness O2MA lines. Overall, µbs is positively correlated with recombination and proximity to short tandem repeats and negatively correlated with 10 bp and 1 kb GC content. However, probability of mutation is sufficiently predicted by the three-nucleotide motif alone. Approximately 90% of the variance in standing nucleotide variation is explained by mutability. Total mutation rate increased in the O2MA lines, as predicted by the "drift barrier" model of mutation rate evolution. These data, combined with experimental estimates of fitness, suggest that epistasis is synergistic.

Temperature effects on life-history trade-offs, germline maintenance and mutation rate under simulated climate warming

Mutation has a fundamental influence over evolutionary processes, but how evolutionary processes shape mutation rate remains less clear. In asexual unicellular organism, increased mutation rates have been observed in stressful environments and the reigning paradigm ascribes this increase to selection for evolvability. However, this explanation does not apply in sexually reproducing species, where little is known about how the environment affects mutation rate. Here we challenged experimental lines of seed beetle, evolved at ancestral temperature or under simulated climate warming, to repair induced mutations at ancestral and stressful temperature. Results show that temperature stress causes individuals to pass on a greater mutation load to their grand-offspring. This suggests that stress-induced mutation rates, in unicellular and multicellular organisms alike, can result from compromised germline DNA repair in low condition individuals. Moreover, lines adapted to simulated climate warming had evolved increased longevity at the cost of reproduction, and this allocation decision improved germline repair. These results suggest that mutation rates can be modulated by resource allocation trade-offs encompassing life-history traits and the germline and have important implications for rates of adaptation and extinction as well as our understanding of genetic diversity in multicellular organisms.

What drives the evolution of condition-dependent recombination in diploids? Some insights from simulation modelling

While the evolutionary advantages of non-zero recombination rates have prompted diverse theoretical explanations, the evolution of essential recombination features remains underexplored. We focused on one such feature, the condition dependence of recombination, viewed as the variation in within-generation sensitivity of recombination to external (environment) and/or internal (genotype) conditions. Limited empirical evidence for its existence comes mainly from diploids, whereas theoretical models show that it only easily evolves in haploids. The evolution of condition-dependent recombination can be explained by its advantage for the selected system (indirect effect), or by benefits to modifier alleles, ensuring this strategy regardless of effects on the selected system (direct effect). We considered infinite panmictic populations of diploids exposed to a cyclical two-state environment. Each organism had three selected loci. Examining allele dynamics at a fourth, selectively neutral recombination modifier locus, we frequently observed that a modifier allele conferring condition-dependent recombination between the selected loci displaced the allele conferring the optimal constant recombination rate. Our simulations also confirm the results of theoretical studies showing that condition-dependent recombination cannot evolve in diploids on the basis of direct fitness-dependent effects alone. Therefore, the evolution of condition-dependent recombination in diploids can be driven by indirect effects alone, i.e. by modifier effects on the selected system.This article is part of the themed issue 'Evolutionary causes and consequences of recombination rate variation in sexual organisms'.

Errors in mutagenesis and the benefit of cell-to-cell signalling in the evolution of stress-induced mutagenesis

Stress-induced mutagenesis is a widely observed phenomenon. Theoretical models have shown that stress-induced mutagenesis can be favoured by natural selection due to the beneficial mutations it generates. These models, however, assumed an error-free regulation of mutation rate in response to stress. Here, we explore the effects of errors in the regulation of mutagenesis on the evolution of stress-induced mutagenesis, and consider the role of cell-to-cell signalling. Using theoretical models, we show (i) that stress-induced mutagenesis can be disadvantageous if errors are common; and (ii) that cell-to-cell signalling can allow stress-induced mutagenesis to be favoured by selection even when error rates are high. We conclude that cell-to-cell signalling can facilitate the evolution of stress-induced mutagenesis in microbes through second-order selection.

Coordinated Changes in Mutation and Growth Rates Induced by Genome Reduction

Genome size is determined during evolution, but it can also be altered by genetic engineering in laboratories. The systematic characterization of reduced genomes provides valuable insights into the cellular properties that are quantitatively described by the global parameters related to the dynamics of growth and mutation. In the present study, we analyzed a small collection of W3110 Escherichia coli derivatives containing either the wild-type genome or reduced genomes of various lengths to examine whether the mutation rate, a global parameter representing genomic plasticity, was affected by genome reduction. We found that the mutation rates of these cells increased with genome reduction. The correlation between genome length and mutation rate, which has been reported for the evolution of bacteria, was also identified, intriguingly, for genome reduction. Gene function enrichment analysis indicated that the deletion of many of the genes encoding membrane and transport proteins play a role in the mutation rate changes mediated by genome reduction. Furthermore, the increase in the mutation rate with genome reduction was highly associated with a decrease in the growth rate in a nutrition-dependent manner; thus, poorer media showed a larger change that was of higher significance. This negative correlation was strongly supported by experimental evidence that the serial transfer of the reduced genome improved the growth rate and reduced the mutation rate to a large extent. Taken together, the global parameters corresponding to the genome, growth, and mutation showed a coordinated relationship, which might be an essential working principle for balancing the cellular dynamics appropriate to the environment.

Genome reduction is a powerful approach for investigating the fundamental rules for living systems. Whether genetically disturbed genomes have any specific properties that are different from or similar to those of natively evolved genomes has been under investigation. In the present study, we found that Escherichia coli cells with reduced genomes showed accelerated nucleotide substitution errors (mutation rates), although these cells retained the normal DNA mismatch repair systems. Intriguingly, this finding of correlation between reduced genome size and a higher mutation rate was consistent with the reported evolution of mutation rates. Furthermore, the increased mutation rate was quantitatively associated with a decreased growth rate, indicating that the global parameters related to the genome, growth, and mutation, which represent the amount of genetic information, the efficiency of propagation, and the fidelity of replication, respectively, are dynamically coordinated.

Hsp90 and Physiological Stress Are Linked to Autonomous Transposon Mobility and Heritable Genetic Change in Nematodes

Transposable elements (TEs) have been recognized as potentially powerful drivers of genomic evolutionary change, but factors affecting their mobility and regulation remain poorly understood. Chaperones such as Hsp90 buffer environmental perturbations by regulating protein conformation, but are also part of the PIWI-interacting RNA pathway, which regulates genomic instability arising from mobile TEs in the germline. Stress-induced mutagenesis from TE movement could thus arise from functional trade-offs in the dual roles of Hsp90. We examined the functional constraints of Hsp90 and its role as a regulator of TE mobility by exposing nematodes (Caenorhabditis elegans and Caenorhabditis briggsae) to environmental stress, with and without RNAi-induced silencing of Hsp90. TE excision frequency increased with environmental stress intensity at multiple loci in several strains of each species. These effects were compounded by RNAi-induced knockdown of Hsp90. Mutation frequencies at the unc-22 marker gene in the offspring of animals exposed to environmental stress and Hsp90 RNAi mirrored excision frequency in response to these treatments. Our results support a role for Hsp90 in the suppression of TE mobility, and demonstrate that that the Hsp90 regulatory pathway can be overwhelmed with moderate environmental stress. By compromising genomic stability in germline cells, environmentally induced mutations arising from TE mobility and insertion can have permanent and heritable effects on both the phenotype and the genotype of subsequent generations.

Low Genetic Quality Alters Key Dimensions of the Mutational Spectrum

Mutations affect individual health, population persistence, adaptation, diversification, and genome evolution. There is evidence that the mutation rate varies among genotypes, but the causes of this variation are poorly understood. Here, we link differences in genetic quality with variation in spontaneous mutation in a Drosophila mutation accumulation experiment. We find that chromosomes maintained in low-quality genetic backgrounds experience a higher rate of indel mutation and a lower rate of gene conversion in a manner consistent with condition-based differences in the mechanisms used to repair DNA double strand breaks. These aspects of the mutational spectrum were also associated with body mass, suggesting that the effect of genetic quality on DNA repair was mediated by overall condition, and providing a mechanistic explanation for the differences in mutational fitness decline among these genotypes. The rate and spectrum of substitutions was unaffected by genetic quality, but we find variation in the probability of substitutions and indels with respect to several aspects of local sequence context, particularly GC content, with implications for models of molecular evolution and genome scans for signs of selection. Our finding that the chances of mutation depend on genetic context and overall condition has important implications for how sequences evolve, the risk of extinction, and human health.

Stress-induced mutagenesis and complex adaptation

Because mutations are mostly deleterious, mutation rates should be reduced by natural selection. However, mutations also provide the raw material for adaptation. Therefore, evolutionary theory suggests that the mutation rate must balance between adaptability-the ability to adapt-and adaptedness-the ability to remain adapted. We model an asexual population crossing a fitness valley and analyse the rate of complex adaptation with and without stress-induced mutagenesis (SIM)-the increase of mutation rates in response to stress or maladaptation. We show that SIM increases the rate of complex adaptation without reducing the population mean fitness, thus breaking the evolutionary trade-off between adaptability and adaptedness. Our theoretical results support the hypothesis that SIM promotes adaptation and provide quantitative predictions of the rate of complex adaptation with different mutational strategies.

Stress-induced hypermutation as a physical property of life, a force of natural selection and its role in four thought experiments

The independence of genetic mutation rate from selection is central to neo-Darwinian evolutionary theory. However, it has been continuously challenged for more than 30 years by experimental evidence of genetic mutation rate transiently increasing in response to stress (stress-induced hypermutation, SIH). The prominent concept of evolved evolvability (EE) explains that natural selection for strategies more competitive at evolutionary adaptation itself gives rise to mechanisms dynamically adjusting mutation rates to environmental stress. Here, we theoretically investigate the alternative (not mutually exclusive) hypothesis that SIH is an inherent physical property of all genetically reproducing life. We define stress as any condition lowering the capability of utilizing metabolic resources for genome storage and replication. This thermodynamical analysis indicates stress-induced increases in the genetic mutation rate in genome storage and in genome replication as inherent physical properties of genetically reproducing life. Further integrating SIH into an overall organismic thermodynamic budget identifies SIH as a force of natural selection, alongside death rate, replication rate and constitutive mutation rate differences. We execute four thought experiments with a non-recombinant lesion mutant strain to predict experimental observations due to SIH in response to different stresses and stress combinations. We find (1) acceleration of adaptation over models without SIH, (2) possibility of adaptation at high stresses which are not explicable by mutation in genome replication alone and (3) different adaptive potential under high growth-inhibiting versus high lethal stresses. The predictions are directly comparable to culture experiments (colony size time courses, antibacterial resistance assay and occurrence of lesion-reversion mutant colonies) and genome sequence analysis. Considering suggestions of drug-mediated disruption of SIH and attempts to target mutation-associated sites with chemotherapeutic agents to prevent resistance, our findings seem to be relevant knowledge for resistance-averse drug development and administration.

Genomic background and generation time influence deleterious mutation rates in Daphnia

Understanding how genetic variation is generated and how selection shapes mutation rates over evolutionary time requires knowledge of the factors influencing mutation and its effects on quantitative traits. We explore the impact of two factors, genomic background and generation time, on deleterious mutation in Daphnia pulicaria, a cyclically parthenogenic aquatic microcrustacean, using parallel mutation-accumulation experiments. The deleterious mutational properties of life-history characters for individuals from two different populations, and for individuals maintained at two different generation times, were quantified and compared. Mutational properties varied between populations, especially for clutch size, suggesting that genomic background influences mutational properties for some characters. Generation time was found to have a greater effect on mutational properties, with higher per-generation deleterious mutation rates in lines with longer generation times. These results suggest that differences in genetic architecture among populations and species may be explained in part by demographic features that significantly influence generation time and therefore the rate of mutation.

Mutation Load: The Fitness of Individuals in Populations Where Deleterious Alleles Are Abundant

Many multicellular eukaryotes have reasonably high per-generation mutation rates. Consequently, most populations harbor an abundance of segregating deleterious alleles. These alleles, most of which are of small effect individually, collectively can reduce substantially the fitness of individuals relative to what it would be otherwise; this is mutation load. Mutation load can be lessened by any factor that causes more mutations to be removed per selective death, such as inbreeding, synergistic epistasis, population structure, or harsh environments. The ecological effects of load are not clear-cut because some conditions (such as selection early in life, sexual selection, reproductive compensation, and intraspecific competition) reduce the effects of load on population size and persistence, but other conditions (such as interspecific competition and load on resource use efficiency) can cause small amounts of load to have strong effects on the population, even extinction. We suggest a series of studies to improve our understanding of the effects of mutation load.

Invariance (?) of mutational parameters for relative fitness over 400 generations of mutation accumulation in Caenorhabditis elegans

Evidence is accumulating that individuals in poor physiologic condition may accumulate mutational damage faster than individuals in good condition. If poor condition results from pre-existing deleterious mutations, the result is "fitness-dependent mutation rate," which has interesting theoretical implications. Here we report a study in which 10 mutation accumulation (MA) lines of the nematode Caenorhabditis elegans that had previously accumulated mutations for 250 generations under relaxed selection were expanded into sets of "second-order" MA lines and allowed to accumulate mutations for an additional 150 generations. The 10 lines were chosen on the basis of the relative change in fitness over the first 250 generations of MA, five high-fitness lines and five low-fitness lines. On average, the mutational properties (per-generation change in mean relative fitness, mutational variance, and Bateman-Mukai estimates of genomic mutation rate and average mutational effect) of the high-fitness and low-fitness did not differ significantly, and averaged over all lines, the point estimates were extremely close to those of the first-order MA experiment after 200 generations of MA. However, several nonsignificant trends indicate that low-fitness lines may in fact be more likely to suffer mutational damage than high-fitness lines.

DNA repair pathway choice is influenced by the health of Drosophila melanogaster

In nature, individuals vary tremendously in condition and this may be an important source of variation in mutation rate. Condition is likely to affect cell state and thereby impact the amount of DNA damage sustained and/or the way it is repaired. Here, we focus on DNA repair. If low-condition individuals are less capable of devoting the same level of resources to accurate repair, they may suffer higher mutation rates. However, repair decisions are also governed by various aspects of cell physiology, which may render the prediction that "higher-condition individuals use better repair mechanisms" too simplistic. We use a larval diet manipulation in Drosophila melanogaster to create high- and low-condition individuals and then contrast their relative usage of three repair pathways [homologous recombination (HR), single-strand annealing (SSA), and nonhomologous end joining (NHEJ)] that differ in their mechanistic requirements and their mutational consequences. We find that low-condition flies are more likely than high-condition flies to use the most conservative of these three repair pathways, suggesting that physiological constraints on repair pathway usage may be more important than energetic costs. We also show that the repair differences between high- and low-condition flies resemble those between young and old flies, suggesting the underlying mechanisms may be similar. Finally, we observe that the effect of larval diet on adult repair increases as flies age, indicating that developmental differences early in life can have long-lasting consequences.

Evidence for elevated mutation rates in low-quality genotypes

The deleterious mutation rate plays a key role in a number of important topics in biology, from mating system evolution to human health. Despite this broad significance, the nature and causes of variation in mutation rate are poorly understood, especially in multicellular organisms. We test whether genetic quality, the presence or absence of deleterious alleles, affects the mutation rate in Drosophila melanogaster by using a modified mutation accumulation approach. We find evidence that genotypes constructed to carry deleterious "treatment" alleles on one chromosome during mutation accumulation experience an elevated mutation rate on a different chromosome. Further, this elevation is correlated with the effect of the treatment alleles on phenotypic condition, measured as body mass. Treatment alleles that reduce mass by 10% cause a doubling in the rate of mutational decline. Our results show that mutation rates are sensitive to genetic stress, such that individuals with low-quality genotypes will produce offspring of even lower genetic quality, in a mutational positive feedback loop. This type of variation in mutation rate is expected to alter a variety of predictions based on mutation load theory and accelerate adaptation to new environments. Positive mutational feedback could affect human health by increasing the rate of germline mutation, and possibly somatic mutation, in individuals of poor health because of genetic or environmental stress.

The evolution of stress-induced hypermutation in asexual populations

Numerous empirical studies show that stress of various kinds induces a state of hypermutation in bacteria via multiple mechanisms, but theoretical treatment of this intriguing phenomenon is lacking. We used deterministic and stochastic models to study the evolution of stress-induced hypermutation in infinite and finite-size populations of bacteria undergoing selection, mutation, and random genetic drift in constant environments and in changing ones. Our results suggest that if beneficial mutations occur, even rarely, then stress-induced hypermutation is advantageous for bacteria at both the individual and the population levels and that it is likely to evolve in populations of bacteria in a wide range of conditions because it is favored by selection. These results imply that mutations are not, as the current view holds, uniformly distributed in populations, but rather that mutations are more common in stressed individuals and populations. Because mutation is the raw material of evolution, these results have a profound impact on broad aspects of evolution and biology.