Kick-starting the ratchet: the fate of mutators in an asexual population

Endosymbiotic ratchet accelerates divergence after organelle origin: The Paulinella model for plastid evolution: The Paulinella model for plastid evolution

We hypothesize that as one of the most consequential events in evolution, primary endosymbiosis accelerates lineage divergence, a process we refer to as the endosymbiotic ratchet. Our proposal is supported by recent work on the photosynthetic amoeba, Paulinella, that underwent primary plastid endosymbiosis about 124 Mya. This amoeba model allows us to explore the early impacts of photosynthetic organelle (plastid) origin on the host lineage. The current data point to a central role for effective population size (N_e ) in accelerating divergence post-endosymbiosis due to limits to dispersal and reproductive isolation that reduce N_e , leading to local adaptation. We posit that isolated populations exploit different strategies and behaviors and assort themselves in non-overlapping niches to minimize competition during the early, rapid evolutionary phase of organelle integration. The endosymbiotic ratchet provides a general framework for interpreting post-endosymbiosis lineage evolution that is driven by disruptive selection and demographic and population shifts. Also see the video abstract here: https://youtu.be/gYXrFM6Zz6Q.

Evolutionary dynamics and genomic features of the Elizabethkingia anophelis 2015 to 2016 Wisconsin outbreak strain

An atypically large outbreak of Elizabethkingia anophelis infections occurred in Wisconsin. Here we show that it was caused by a single strain with thirteen characteristic genomic regions. Strikingly, the outbreak isolates show an accelerated evolutionary rate and an atypical mutational spectrum. Six phylogenetic sub-clusters with distinctive temporal and geographic dynamics are revealed, and their last common ancestor existed approximately one year before the first recognized human infection. Unlike other E. anophelis, the outbreak strain had a disrupted DNA repair mutY gene caused by insertion of an integrative and conjugative element. This genomic change probably contributed to the high evolutionary rate of the outbreak strain and may have increased its adaptability, as many mutations in protein-coding genes occurred during the outbreak. This unique discovery of an outbreak caused by a naturally occurring mutator bacterial pathogen provides a dramatic example of the potential impact of pathogen evolutionary dynamics on infectious disease epidemiology.

Elizabethkingia anophelis is an emerging pathogen of high antimicrobial resistance. Perrin and colleagues sequenced isolates of a 2015/2016 E. anophelis outbreak in Wisconsin and found substantial genetic diversity, accelerated evolutionary rate and a disruptive mutation in the DNA repair gene mutY.

Dynamics and Fate of Beneficial Mutations Under Lineage Contamination by Linked Deleterious Mutations

Beneficial mutations drive adaptive evolution, yet their selective advantage does not ensure their fixation. Haldane's application of single-type branching process theory showed that genetic drift alone could cause the extinction of newly arising beneficial mutations with high probability. With linkage, deleterious mutations will affect the dynamics of beneficial mutations and might further increase their extinction probability. Here, we model the lineage dynamics of a newly arising beneficial mutation as a multitype branching process. Our approach accounts for the combined effects of drift and the stochastic accumulation of linked deleterious mutations, which we call lineage contamination We first study the lineage-contamination phenomenon in isolation, deriving dynamics and survival probabilities (the complement of extinction probabilities) of beneficial lineages. We find that survival probability is zero when [Formula: see text] where U is deleterious mutation rate and [Formula: see text] is the selective advantage of the beneficial mutation in question, and is otherwise depressed below classical predictions by a factor bounded from below by [Formula: see text] We then put the lineage contamination phenomenon into the context of an evolving population by incorporating the effects of background selection. We find that, under the combined effects of lineage contamination and background selection, ensemble survival probability is never zero but is depressed below classical predictions by a factor bounded from below by [Formula: see text] where [Formula: see text] is mean selective advantage of beneficial mutations, and [Formula: see text] This factor, and other bounds derived from it, are independent of the fitness effects of deleterious mutations. At high enough mutation rates, lineage contamination can depress fixation probabilities to values that approach zero. This fact suggests that high mutation rates can, perhaps paradoxically, (1) alleviate competition among beneficial mutations, or (2) potentially even shut down the adaptive process. We derive critical mutation rates above which these two events become likely.

Role of epistasis on the fixation probability of a non-mutator in an adapted asexual population

The mutation rate of a well adapted population is prone to reduction so as to have a lower mutational load. We aim to understand the role of epistatic interactions between the fitness affecting mutations in this process. Using a multitype branching process, the fixation probability of a single non-mutator emerging in a large asexual mutator population is analytically calculated here. The mutator population undergoes deleterious mutations at constant, but at a much higher rate than that of the non-mutator. We find that antagonistic epistasis lowers the chances of mutation rate reduction, while synergistic epistasis enhances it. Below a critical value of epistasis, the fixation probability behaves non-monotonically with variation in the mutation rate of the background population. Moreover, the variation of this critical value of the epistasis parameter with the strength of the mutator is discussed in the appendix. For synergistic epistasis, when selection is varied, the fixation probability reduces overall, with damped oscillations.

The rich phase structure of a mutator model

We propose a modification of the Crow-Kimura and Eigen models of biological molecular evolution to include a mutator gene that causes both an increase in the mutation rate and a change in the fitness landscape. This mutator effect relates to a wide range of biomedical problems. There are three possible phases: mutator phase, mixed phase and non-selective phase. We calculate the phase structure, the mean fitness and the fraction of the mutator allele in the population, which can be applied to describe cancer development and RNA viruses. We find that depending on the genome length, either the normal or the mutator allele dominates in the mixed phase. We analytically solve the model for a general fitness function. We conclude that the random fitness landscape is an appropriate choice for describing the observed mutator phenomenon in the case of a small fraction of mutators. It is shown that the increase in the mutation rates in the regular and the mutator parts of the genome should be set independently; only some combinations of these increases can push the complex biomedical system to the non-selective phase, potentially related to the eradication of tumors.

Evolution of Mutation Rates in Rapidly Adapting Asexual Populations

Mutator and antimutator alleles often arise and spread in both natural microbial populations and laboratory evolution experiments. The evolutionary dynamics of these mutation rate modifiers are determined by indirect selection on linked beneficial and deleterious mutations. These indirect selection pressures have been the focus of much earlier theoretical and empirical work, but we still have a limited analytical understanding of how the interplay between hitchhiking and deleterious load influences the fates of modifier alleles. Our understanding is particularly limited when clonal interference is common, which is the regime of primary interest in laboratory microbial evolution experiments. Here, we calculate the fixation probability of a mutator or antimutator allele in a rapidly adapting asexual population, and we show how this quantity depends on the population size, the beneficial and deleterious mutation rates, and the strength of a typical driver mutation. In the absence of deleterious mutations, we find that clonal interference enhances the fixation probability of mutators, even as they provide a diminishing benefit to the overall rate of adaptation. When deleterious mutations are included, natural selection pushes the population toward a stable mutation rate that can be suboptimal for the adaptation of the population as a whole. The approach to this stable mutation rate is not necessarily monotonic: even in the absence of epistasis, selection can favor mutator and antimutator alleles that "overshoot" the stable mutation rate by substantial amounts.

Fixation of mutators in asexual populations: the role of genetic drift and epistasis

We study the evolutionary dynamics of an asexual population of nonmutators and mutators on a class of epistatic fitness landscapes. We consider the situation in which all mutations are deleterious and mutators are produced from nonmutators continually at a constant rate. We find that in an infinitely large population, a minimum nonmutator-to-mutator conversion rate is required to fix the mutators but an arbitrarily small conversion rate results in the fixation of mutators in a finite population. We calculate analytical expressions for the mutator fraction at mutation-selection balance and fixation time for mutators in a finite population when the difference between the mutation rate for mutator and nonmutator is smaller (regime I) and larger (regime II) than the selection coefficient. Our main result is that in regime I, the mutator fraction and the fixation time are independent of epistasis but in regime II, mutators are rarer and take longer to fix when the decrease in fitness with the number of deleterious mutations occurs at an accelerating rate (synergistic epistasis) than at a diminishing rate (antagonistic epistasis). Our analytical results are compared with numerics and their implications are discussed.

The balance between mutators and nonmutators in asexual populations

Mutator alleles, which elevate an individual's mutation rate from 10 to 10,000-fold, have been found at high frequencies in many natural and experimental populations. Mutators are continually produced from nonmutators, often due to mutations in mismatch-repair genes. These mutators gradually accumulate deleterious mutations, limiting their spread. However, they can occasionally hitchhike to high frequencies with beneficial mutations. We study the interplay between these effects. We first analyze the dynamics of the balance between the production of mutator alleles and their elimination due to deleterious mutations. We find that when deleterious mutation rates are high in mutators, there will often be many "young," recently produced mutators in the population, and the fact that deleterious mutations only gradually eliminate individuals from a population is important. We then consider how this mutator-nonmutator balance can be disrupted by beneficial mutations and analyze the circumstances in which fixation of mutator alleles is likely. We find that dynamics is crucial: even in situations where selection on average acts against mutators, so they cannot stably invade, the mutators can still occasionally generate beneficial mutations and hence be important to the evolution of the population.

Fitness-dependent mutation rates in finite populations

Mutation rate may be condition dependent, whereby individuals in poor condition, perhaps from high mutation load, have higher mutation rates than individuals in good condition. Agrawal (J. Evol. Biol.15, 2002, 1004) explored the basic properties of fitness-dependent mutation rate (FDMR) in infinite populations and reported some heuristic results for finite populations. The key parameter governing how infinite populations evolve under FDMR is the curvature (k) of the relationship between fitness and mutation rate. We extend Agrawal's analysis to finite populations and consider dominance and epistasis. In finite populations, the probability of long-term existence depends on k. In sexual populations, positive curvature leads to low equilibrium mutation rate, whereas negative curvature results in high mutation rate. In asexual populations, negative curvature results in rapid extinction via 'mutational meltdown', whereas positive curvature sometimes allows persistence. We speculate that fitness-dependent mutation rate may provide the conditions for genetic architecture to diverge between sexual and asexual taxa.