Virus mutators and antimutators: roles in evolution, pathogenesis and emergence

Molecular Biology and Diversification of Human Retroviruses

Studies of retroviruses have led to many extraordinary discoveries that have advanced our understanding of not only human diseases, but also molecular biology as a whole. The most recognizable human retrovirus, human immunodeficiency virus type 1 (HIV-1), is the causative agent of the global AIDS epidemic and has been extensively studied. Other human retroviruses, such as human immunodeficiency virus type 2 (HIV-2) and human T-cell leukemia virus type 1 (HTLV-1), have received less attention, and many of the assumptions about the replication and biology of these viruses are based on knowledge of HIV-1. Existing comparative studies on human retroviruses, however, have revealed that key differences between these viruses exist that affect evolution, diversification, and potentially pathogenicity. In this review, we examine current insights on disparities in the replication of pathogenic human retroviruses, with a particular focus on the determinants of structural and genetic diversity amongst HIVs and HTLV.

An MHV-68 Mutator Phenotype Mutant Virus, Confirmed by CRISPR/Cas9-Mediated Gene Editing of the Viral DNA Polymerase Gene, Shows Reduced Viral Fitness

Drug resistance studies on human γ-herpesviruses are hampered by the absence of an in vitro system that allows efficient lytic viral replication. Therefore, we employed murine γ-herpesvirus-68 (MHV-68) that efficiently replicates in vitro as a model to study the antiviral resistance of γ-herpesviruses. In this study, we investigated the mechanism of resistance to nucleoside (ganciclovir (GCV)), nucleotide (cidofovir (CDV), HPMP-5azaC, HPMPO-DAPy) and pyrophosphate (foscarnet (PFA)) analogues and the impact of these drug resistance mutations on viral fitness. Viral fitness was determined by dual infection competition assays, where MHV-68 drug-resistant viral clones competed with the wild-type virus in the absence and presence of antivirals. Using next-generation sequencing, the composition of the viral populations was determined at the time of infection and after 5 days of growth. Antiviral drug resistance selection resulted in clones harboring mutations in the viral DNA polymerase (DP), denoted Y383S^GCV, Q827R^HPMP-5azaC, G302W^PFA, K442T^PFA, G302W+K442T^PFA, C297W^HPMPO-DAPy and C981Y^CDV. Without antiviral pressure, viral clones Q827R^HPMP-5azaC, G302W^PFA, K442T^PFA and G302W+K442T^PFA grew equal to the wild-type virus. However, in the presence of antivirals, these mutants had a growth advantage over the wild-type virus that was moderately to very strongly correlated with antiviral resistance. The Y383S^GCV mutant was more fit than the wild-type virus with and without antivirals, except in the presence of brivudin. The C297W and C981Y changes were associated with a mutator phenotype and had a severely impaired viral fitness in the absence and presence of antivirals. The mutator phenotype caused by C297W in MHV-68 DP was validated by using a CRISPR/Cas9 genome editing approach.

A proofreading-impaired herpesvirus generates populations with quasispecies-like structure

RNA virus populations are composed of highly diverse individuals that form a cloud of related sequences commonly referred to as a 'quasispecies'^1-3. This diversity arises as a consequence of low-fidelity genome replication^4,5. By contrast, DNA virus populations contain more uniform individuals with similar fitness⁶. Genome diversity is often correlated with increased fitness in RNA viruses, while DNA viruses are thought to require more faithful genome replication. During DNA replication, erroneously incorporated bases are removed by a 3'-5' exonuclease, a highly conserved enzymatic function of replicative DNA but not RNA polymerases. This proofreading process enhances replication fidelity and ensures the genome integrity of DNA organisms, including large DNA viruses⁷. Here, we show that a herpesvirus can tolerate impaired exonucleolytic proofreading, resulting in DNA virus populations, which, as in RNA viruses⁸, are composed of highly diverse genotypes of variable individual fitness. This indicates that herpesvirus mutant diversity may compensate for individual fitness loss. Notably, in vivo infection with diverse virus populations results in a marked increase in virulence compared to genetically homogenous parental virus. While we cannot exclude that the increase in virulence is caused by selection of and/or interactions between individual genotypes, our findings are consistent with quasispecies dynamics. Our results contrast with traditional views of DNA virus replication and evolution, and indicate that a substantial increase in population diversity can lead to higher virulence.

Calicivirus RNA-Dependent RNA Polymerases: Evolution, Structure, Protein Dynamics, and Function

The Caliciviridae are viruses with a positive-sense, single-stranded RNA genome that is packaged into an icosahedral, environmentally stable protein capsid. The family contains five genera (Norovirus, Nebovirus, Sapovirus, Lagovirus, and Vesivirus) that infect vertebrates including amphibians, reptiles, birds, and mammals. The RNA-dependent RNA polymerase (RdRp) replicates the genome of RNA viruses and can speed up evolution due to its error-prone nature. Studying calicivirus RdRps in the context of genuine virus replication is often hampered by a lack of suitable model systems. Enteric caliciviruses and RHDV in particular are notoriously difficult to propagate in cell culture; therefore, molecular studies of replication mechanisms are challenging. Nevertheless, research on recombinant proteins has revealed several unexpected characteristics of calicivirus RdRps. For example, the RdRps of RHDV and related lagoviruses possess the ability to expose a hydrophobic motif, to rearrange Golgi membranes, and to copy RNA at unusually high temperatures. This review is focused on the structural dynamics, biochemical properties, kinetics, and putative interaction partners of these RdRps. In addition, we discuss the possible existence of a conserved but as yet undescribed structural element that is shared amongst the RdRps of all caliciviruses.

Complexities of Viral Mutation Rates

Many viruses evolve rapidly. This is due, in part, to their high mutation rates. Mutation rate estimates for over 25 viruses are currently available. Here, we review the population genetics of virus mutation rates. We specifically cover the topics of mutation rate estimation, the forces that drive the evolution of mutation rates, and how the optimal mutation rate can be context-dependent.

Antigenic Drift Defines a New D4 Subgenotype of Measles Virus

The measles virus hemagglutinin (MeV-H) protein is the main target of protective neutralizing antibodies. Using a panel of monoclonal antibodies (MAbs) that recognize known major antigenic sites in MeV-H, we identified a D4 genotype variant that escapes neutralization by MAbs targeting the neutralizing epitope (NE) antigenic site. By site-directed mutagenesis, L249P was identified as the critical mutation disrupting the NE in this genotype D4 variant. Forty-two available D4 genotype gene sequences were subsequently analyzed and divided into 2 groups according to the presence or absence of the L249P MeV-H mutation. Further analysis of the MeV-N gene sequences of these 2 groups confirmed that they represent clearly definable, sequence-divergent D4 subgenotypes, which we named subgenotypes D4.1 and D4.2. The subgenotype D4.1 MeVs were isolated predominantly in Kenya and Ethiopia, whereas the MAb-resistant subgenotype D4.2 MeVs were isolated predominantly in France and Great Britain, countries with higher vaccine coverage rates. Interestingly, D4.2 subgenotype viruses showed a trend toward diminished susceptibility to neutralization by human sera pooled from approximately 60 to 80 North American donors. Escape from MAb neutralization may be a powerful epidemiological surveillance tool to monitor the evolution of new MeV subgenotypes.IMPORTANCE Measles virus is a paradigmatic RNA virus, as the antigenic composition of the vaccination has not needed to be updated since its discovery. The vaccine confers protection by inducing neutralizing antibodies that interfere with the function of the hemagglutinin protein. Viral strains are indistinguishable serologically, although characteristic nucleotide sequences differentiate 24 genotypes. In this work, we describe a distant evolutionary branch within genotype D4. Designated subgenotype D4.2, this virus is distinguishable by neutralization with vaccine-induced monoclonal antibodies that target the neutralizing epitope (NE). The subgenotype D4.2 viruses have a higher predominance in countries with intermediary levels of vaccine coverage. Our studies demonstrate that subgenotype D4.2 lacks epitopes associated with half of the known antigenic sites, which significantly impacts our understanding of measles virus evolution.

Population Heterogeneity in Mutation Rate Increases the Frequency of Higher-Order Mutants and Reduces Long-Term Mutational Load

Mutation rate is a crucial evolutionary parameter that has typically been treated as a constant in population genetic analyses. However, the propensity to mutate is likely to vary among co-existing individuals within a population, due to genetic polymorphisms, heterogeneous environmental influences, and random physiological fluctuations. We review the evidence for mutation rate heterogeneity and explore its consequences by extending classic population genetic models to allow an arbitrary distribution of mutation rate among individuals, either with or without inheritance. With this general new framework, we rigorously establish the effects of heterogeneity at various evolutionary timescales. In a single generation, variation of mutation rate about the mean increases the probability of producing zero or many simultaneous mutations on a genome. Over multiple generations of mutation and selection, heterogeneity accelerates the appearance of both deleterious and beneficial multi-point mutants. At mutation-selection balance, higher-order mutant frequencies are likewise boosted, while lower-order mutants exhibit subtler effects; nonetheless, population mean fitness is always enhanced. We quantify the dependencies on moments of the mutation rate distribution and selection coefficients, and clarify the role of mutation rate inheritance. While typical methods of estimating mutation rate will recover only the population mean, analyses assuming mutation rate is fixed to this mean could underestimate the potential for multi-locus adaptation, including medically relevant evolution in pathogenic and cancerous populations. We discuss the potential to empirically parameterize mutation rate distributions, which have to date hardly been quantified.

Mechanisms of viral mutation

The remarkable capacity of some viruses to adapt to new hosts and environments is highly dependent on their ability to generate de novo diversity in a short period of time. Rates of spontaneous mutation vary amply among viruses. RNA viruses mutate faster than DNA viruses, single-stranded viruses mutate faster than double-strand virus, and genome size appears to correlate negatively with mutation rate. Viral mutation rates are modulated at different levels, including polymerase fidelity, sequence context, template secondary structure, cellular microenvironment, replication mechanisms, proofreading, and access to post-replicative repair. Additionally, massive numbers of mutations can be introduced by some virus-encoded diversity-generating elements, as well as by host-encoded cytidine/adenine deaminases. Our current knowledge of viral mutation rates indicates that viral genetic diversity is determined by multiple virus- and host-dependent processes, and that viral mutation rates can evolve in response to specific selective pressures.

Effect of mismatch repair on the mutation rate of bacteriophage ϕX174

Viral mutation rates vary widely in nature, yet the mechanistic and evolutionary determinants of this variability remain unclear. Small DNA viruses mutate orders of magnitude faster than their hosts despite using host-encoded polymerases for replication, which suggests these viruses may avoid post-replicative repair. Supporting this, the genome of bacteriophage ϕX174 is completely devoid of GATC sequence motifs, which are required for methyl-directed mismatch repair in Escherichia coli. Here, we show that restoration of the randomly expected number of GATC sites leads to an eightfold reduction in the rate of spontaneous mutation of the phage, without severely impairing its replicative capacity over the short term. However, the efficacy of mismatch repair in the presence of GATC sites is limited by inefficient methylation of the viral DNA. Therefore, both GATC avoidance and DNA under-methylation elevate the mutation rate of the phage relative to that of the host. We also found that the effects of GATC sites on the phage mutation rate vary extensively depending on their specific location within the phage genome. Finally, the mutation rate reduction afforded by GATC sites is fully reverted under stress conditions, which up-regulate repair pathways and expression of error-prone host polymerases such as heat and treatment with the base analog 5-fluorouracil, suggesting that access to repair renders the phage sensitive to stress-induced mutagenesis.

Expected Effect of Deleterious Mutations on Within-Host Adaptation of Pathogens

Adaptation is a common theme in both pathogen emergence, for example, in zoonotic cross-species transmission, and pathogen control, where adaptation might limit the effect of the immune response and antiviral treatment. When such evolution requires deleterious intermediate mutations, fitness ridges and valleys arise in the pathogen's fitness landscape. The effect of deleterious intermediate mutations on within-host pathogen adaptation is examined with deterministic calculations, appropriate for pathogens replicating in large populations with high error rates. The effect of deleterious intermediate mutations on pathogen adaptation is smaller than their name might suggest: when two mutations are required and each individual single mutation is fully deleterious, the pathogen can jump across the fitness valley by obtaining two mutations at once, leading to a proportion of adapted mutants that is 20-fold lower than that in the situation where the fitness of all mutants is neutral. The negative effects of deleterious intermediates are typically substantially smaller and outweighed by the fitness advantages of the adapted mutant. Moreover, requiring a specific mutation order has a substantially smaller effect on pathogen adaptation than the effect of all intermediates being deleterious. These results can be rationalized when the number of routes of mutation available to the pathogen is calculated, providing a simple approach to estimate the effect of deleterious mutations. The calculations discussed here are applicable when the effect of deleterious mutations on the within-host adaptation of pathogens is assessed, for example, in the context of zoonotic emergence, antigenic escape, and drug resistance.

IMPORTANCE

Adaptation is critical for pathogens after zoonotic transmission into a new host species or to achieve antigenic immune escape and drug resistance. Using a deterministic approach, the effects of deleterious intermediate mutations on pathogen adaptation were calculated while avoiding commonly made simplifications that do not apply to large pathogen populations replicating with high mutation rates. Perhaps unexpectedly, pathogen adaptation does not halt when the intermediate mutations are fully deleterious. The negative effects of deleterious mutations are substantially outweighed by the fitness gains of adaptation. To gain an understanding of the effect of deleterious mutations on pathogen adaptation, a simple approach that counts the number of routes available to the pathogen with and without deleterious intermediate mutations is introduced. This methodology enables a straightforward calculation of the proportion of the pathogen population that will cross a fitness valley or traverse a fitness ridge, without reverting to more complicated mathematical models.

Alphavirus mutator variants present host-specific defects and attenuation in mammalian and insect models

Arboviruses cycle through both vertebrates and invertebrates, which requires them to adapt to disparate hosts while maintaining genetic integrity during genome replication. To study the genetic mechanisms and determinants of these processes, we use chikungunya virus (CHIKV), a re-emerging human pathogen transmitted by the Aedes mosquito. We previously isolated a high fidelity (or antimutator) polymerase variant, C483Y, which had decreased fitness in both mammalian and mosquito hosts, suggesting this residue may be a key molecular determinant. To further investigate effects of position 483 on RNA-dependent RNA-polymerase (RdRp) fidelity, we substituted every amino acid at this position. We isolated novel mutators with decreased replication fidelity and higher mutation frequencies, allowing us to examine the fitness of error-prone arbovirus variants. Although CHIKV mutators displayed no major replication defects in mammalian cell culture, they had reduced specific infectivity and were attenuated in vivo. Unexpectedly, mutator phenotypes were suppressed in mosquito cells and the variants exhibited significant defects in RNA synthesis. Consequently, these replication defects resulted in strong selection for reversion during infection of mosquitoes. Since residue 483 is conserved among alphaviruses, we examined the analogous mutations in Sindbis virus (SINV), which also reduced polymerase fidelity and generated replication defects in mosquito cells. However, replication defects were mosquito cell-specific and were not observed in Drosophila S2 cells, allowing us to evaluate the potential attenuation of mutators in insect models where pressure for reversion was absent. Indeed, the SINV mutator variant was attenuated in fruit flies. These findings confirm that residue 483 is a determinant regulating alphavirus polymerase fidelity and demonstrate proof of principle that arboviruses can be attenuated in mammalian and insect hosts by reducing fidelity.

Chikungunya (CHIKV) is a re-emerging mosquito-borne virus that constitutes a major and growing human health burden. Like all RNA viruses, during viral replication CHIKV copies its genome using a polymerase that makes an average of one mistake per replication cycle. Therefore, a single virus generates millions of viral progeny that carry a multitude of distinct mutations in their genomes. In this study, we isolated CHIKV mutators (strains that make more errors than the wildtype virus), to study how higher mutation rates affect fitness in arthropod-borne viruses (arboviruses). CHIKV mutators have reduced virulence in mice and severe replication defects in Aedes mosquito cells. However, these replication defects result in selective pressure for reversion of mutators to a wildtype polymerase in mosquito hosts. To examine how mutators would behave in an insect model in absence of this genetic instability, we isolated mutators of a related virus, Sindbis virus (SINV). SINV mutators had no replication defect in fruit fly (Drosophila) cells, and a SINV mutator strain was stable and attenuated in fruit flies. This work shows proof of principle that arbovirus mutators can exhibit attenuation in both mammalian and insect hosts, and may remain a viable vaccine strategy.

Persistence versus reversion of 3TC resistance in HIV-1 determine the rate of emergence of NVP resistance

When HIV-1 is exposed to lamivudine (3TC) at inhibitory concentrations, resistant variants carrying the reverse transcriptase (RT) substitution M184V emerge rapidly. This substitution confers high-level 3TC resistance and increased RT fidelity. We established a novel in vitro system to study the effect of starting nevirapine (NVP) in 3TC-resistant/NNRTI-naïve clinical isolates, and the impact of maintaining versus dropping 3TC pressure in this setting. Because M184V mutant HIV-1 seems hypersusceptible to adefovir (ADV), we also tested the effect of ADV pressure on the same isolates. We draw four conclusions from our experiments simulating combination therapy in vitro. (1) The presence of low-dose (1 μM) 3TC prevented reversal to wild-type from an M184V mutant background. (2) Adding low-dose 3TC in the presence of NVP delayed the selection of NVP-associated mutations. (3) The presence of ADV, in addition to NVP, led to more rapid reversal to wild-type at position 184 than NVP alone. (4) ADV plus NVP selected for greater numbers of mutations than NVP alone. Inference about the "selection of mutation" is based on two statistical models, one at the viral level, more telling, and the other at the level of predominance of mutation within a population. Multidrug pressure experiments lend understanding to mechanisms of HIV resistance as they bear upon new treatment strategies.

Advances in plant virus evolution: translating evolutionary insights into better disease management

Recent studies in plant virus evolution are revealing that genetic structure and behavior of virus and viroid populations can explain important pathogenic properties of these agents, such as host resistance breakdown, disease severity, and host shifting, among others. Genetic variation is essential for the survival of organisms. The exploration of how these subcellular parasites generate and maintain a certain frequency of mutations at the intra- and inter-host levels is revealing novel molecular virus-plant interactions. They emphasize the role of host environment in the dynamic genetic composition of virus populations. Functional genomics has identified host factors that are transcriptionally altered after virus infections. The analyses of these data by means of systems biology approaches are uncovering critical plant genes specifically targeted by viruses during host adaptation. Also, a next-generation resequencing approach of a whole virus genome is opening new avenues to study virus recombination and the relationships between intra-host virus composition and pathogenesis. Altogether, the analyzed data indicate that systematic disruption of some specific parameters of evolving virus populations could lead to more efficient ways of disease prevention, eradication, or tolerable virus-plant coexistence.

Variable mutation rates as an adaptive strategy in replicator populations

For evolving populations of replicators, there is much evidence that the effect of mutations on fitness depends on the degree of adaptation to the selective pressures at play. In optimized populations, most mutations have deleterious effects, such that low mutation rates are favoured. In contrast to this, in populations thriving in changing environments a larger fraction of mutations have beneficial effects, providing the diversity necessary to adapt to new conditions. What is more, non-adapted populations occasionally benefit from an increase in the mutation rate. Therefore, there is no optimal universal value of the mutation rate and species attempt to adjust it to their momentary adaptive needs. In this work we have used stationary populations of RNA molecules evolving in silico to investigate the relationship between the degree of adaptation of an optimized population and the value of the mutation rate promoting maximal adaptation in a short time to a new selective pressure. Our results show that this value can significantly differ from the optimal value at mutation-selection equilibrium, being strongly influenced by the structure of the population when the adaptive process begins. In the short-term, highly optimized populations containing little variability respond better to environmental changes upon an increase of the mutation rate, whereas populations with a lower degree of optimization but higher variability benefit from reducing the mutation rate to adapt rapidly. These findings show a good agreement with the behaviour exhibited by actual organisms that replicate their genomes under broadly different mutation rates.

5-Azacytidine can induce lethal mutagenesis in human immunodeficiency virus type 1

Ribonucleosides inhibit human immunodeficiency virus type 1 (HIV-1) replication by mechanisms that have not been fully elucidated. Here, we report the antiviral mechanism for the ribonucleoside analog 5-azacytidine (5-AZC). We hypothesized that the anti-HIV-1 activity of 5-AZC was due to an increase in the HIV-1 mutation rate following its incorporation into viral RNA during transcription. However, we demonstrate that 5-AZC's primary antiviral activity can be attributed to its effect on the early phase of HIV-1 replication. Furthermore, the antiviral activity was associated with an increase in the frequency of viral mutants, suggesting that 5-AZC's primary target is reverse transcription. Sequencing analysis showed an enrichment in G-to-C transversion mutations and further supports the idea that reverse transcription is an antiviral target of 5-AZC. These results indicate that 5-AZC is incorporated into viral DNA following reduction to 5-aza-2'-deoxycytidine. Incorporation into the viral DNA leads to an increase in mutant frequency that is consistent with lethal mutagenesis during reverse transcription as the primary antiviral mechanism of 5-AZC. Antiviral activity and increased mutation frequency were also associated with the late phase of HIV-1 replication; however, 5-AZC's effect on the late phase was less robust. These results reveal that the primary antiviral mechanism of 5-AZC can be attributed to its ability to increase the HIV-1 mutation frequency through viral-DNA incorporation during reverse transcription. Our observations indicate that 5-AZC can affect two steps in HIV-1 replication (i.e., transcription and reverse transcription) but that its primary antiviral activity is due to incorporation during reverse transcription.

Rates of evolutionary change in viruses: patterns and determinants

Understanding the factors that determine the rate at which genomes generate and fix mutations provides important insights into key evolutionary mechanisms. We review our current knowledge of the rates of mutation and substitution, as well as their determinants, in RNA viruses, DNA viruses and retroviruses. We show that the high rate of nucleotide substitution in RNA viruses is matched by some DNA viruses, suggesting that evolutionary rates in viruses are explained by diverse aspects of viral biology, such as genomic architecture and replication speed, and not simply by polymerase fidelity.

The evolution of mutation rate in finite asexual populations

In this article, we model analytically the evolution of mutation rate in asexual organisms. Three selective forces are present. First, everything else being equal, individuals with higher mutation rate have a larger fitness, thanks to the energy and time saved by not replicating DNA accurately. Second, as a flip side, the genome of these individuals is replicated with errors that may negatively affect fitness. Third, and conversely, replication errors have a potential benefit if beneficial mutations are to be generated. Our model describes the fate of modifiers of mutation rate under the three forces and allows us to predict the long-term evolutionary trajectory of mutation rate. We obtain three major results. First, in asexuals, the needs for both adaptation and genome preservation are not evolutionary forces that can stabilize mutation rate at an intermediate optimum. When adaptation has a significant role, it primarily destabilizes mutation rate and yields the emergence of strong-effect mutators. Second, in contrast to what is usually believed, the appearance of modifiers with large mutation rate is more likely when the fitness cost of each deleterious mutation is weak, because the cost of replication errors is then paid after a delay. Third, in small populations, and even if adaptations are needed, mutation rate is always blocked at the minimum attainable level, because the rate of adaptation is too slow to play a significant role. Only populations whose size is above a critical mass see their mutation rate affected by the need for adaptation.

Mutagenic outcome of combined antiviral drug treatment during human immunodeficiency virus type 1 replication

The development of antiviral drug resistance is an important problem in the treatment of human immunodeficiency virus type 1 (HIV-1) infection. Potent antiretroviral therapy (ART) is currently used for treatment, and typically consists of at least two reverse transcriptase (RT) inhibitors. To assess the impact of combination therapy on HIV-1 mutagenesis, the mutagenic outcome of combined drug treatment was determined with several different RT drug combinations. Significant increases in HIV-1 mutant frequencies were observed with combinations of nucleoside RT inhibitors as well as in combinations where nucleoside inhibitors were used along with hydroxyurea, a drug known to deplete nucleotide pools in cells. This indicates that combinations of RT drugs can act together to further increase HIV-1 mutant frequencies, which could have important implications for virus population dynamics and could compromise drug therapy regimens.

Mutator dynamics in fluctuating environments

Populations with high mutation rates (mutator clones) are being found in increasing numbers of species, and a clear link is being established between the presence of mutator clones and drug resistance. Mutator clones exist despite the fact that in a constant environment most mutations are deleterious, with the spontaneous mutation rate generally held at a low value. This implies that mutator clones have an important role in the adaptation of organisms to changing environments. Our study examines how mutator dynamics vary according to the frequency of environmental fluctuations. Although recent studies have considered a single environmental switch, here we investigate mutator dynamics in a regularly varying environment, seeking to mimic conditions present, for example, under certain drug or pesticide regimes. Our model provides four significant new insights. First, the results demonstrate that mutators are most prevalent under intermediate rates of environmental change. When the environment oscillates more rapidly, mutators are unable to provide sufficient adaptability to keep pace with the frequent changes in selection pressure and, instead, a population of 'environmental generalists' dominates. Second, our findings reveal that mutator dynamics may be complex, exhibiting limit cycles and chaos. Third, we demonstrate that when each beneficial mutation provides a greater gain in fitness, mutators achieve higher densities in more rapidly fluctuating environments. Fourth, we find that mutators of intermediate strength reach higher densities than very weak or strong mutators.