Clonal Interference and the Evolution of RNA Viruses

See allHide authors and affiliations

Science  10 Sep 1999:
Vol. 285, Issue 5434, pp. 1745-1747
DOI: 10.1126/science.285.5434.1745


In asexual populations, beneficial mutations that occur in different lineages compete with one another. This phenomenon, known as clonal interference, ensures that those beneficial mutations that do achieve fixation are of large effect. Clonal interference also increases the time between fixations, thereby slowing the adaptation of asexual populations. The effects of clonal interference were measured in the asexual RNA virus vesicular stomatitis virus; rates and average effects of beneficial mutations were quantified.

Populations adapt through the appearance and subsequent fixation of beneficial mutations. In a large population, beneficial mutations may arise frequently enough that two or more are simultaneously present in independent lineages. Once beneficial mutations have arisen, there is a certain probability of losing them by drift while their frequency is low. However, after this period dominated by drift, they reach a substantial frequency in the population. For a sexual system, these beneficial mutations will eventually recombine, ensuring their fixation together. If the system is asexual, the lineages created by these beneficial mutations will compete; only the mutation with largest effect will be fixed. Thus, asexual populations must fix beneficial mutations sequentially (1, 2). The possibility of simultaneous fixation of beneficial mutations in sexual populations is often contrasted with the sequential fixation in asexual populations as an argument for the evolutionary advantage of sex (2).

The idea that beneficial mutations must compete in asexual populations was originally proposed by Muller (3), and it has been developed theoretically (4), as well as experimentally demonstrated to be important in determining the rate of adaptation of the bacterium Escherichia coli (5). Gerrish and Lenski (4) modeled the fate of beneficial mutations by considering clonal interference among them as a major factor. The main conclusions of their model were as follows: (i) The probability of fixation of a given beneficial mutation decreases with both population size and mutation rate. (ii) As population size or mutation rate increase, adaptive substitutions result in larger fitness increases. (iii) The rate of adaptation is an increasing, but decelerating, function of both population size and mutation rate. (iv) Beneficial mutations that become transiently common but do not achieve fixation because of interfering beneficial mutations are relatively abundant. (v) Transient polymorphisms may give rise to a “leapfrog” effect, where the most common genotype at a given moment might be less closely related to the immediately preceding one than with an earlier genotype.

RNA viruses show the highest mutation rates in nature (6). This, together with their potentially large effective population sizes and the fact that their reproduction is not obligately sexual, suggests that clonal interference may play an important role in their adaptive evolution. Our goal here is to infer the presence of clonal interference acting on viral populations. Following (4), for increasing population sizes, we predicted that (i) the fitness effect associated with fixed beneficial mutation will tend to be larger and (ii) the rate of adaptation will tend toward a limit.

To detect experimentally the fixation of a beneficial mutation in a viral population, we mixed, at equal proportions, two variants of vesicular stomatitis virus (VSV) that differ only in their ability to grow in the presence of a monoclonal antibody (7). The two variants were selectively equivalent in the absence of monoclonal antibody (7), implying that they should stably coexist until a successful beneficial mutation appears in one of them. Seven different evolutionary regimes were designed, each one differing from the others in effective population size (N e). As shown in Table 1, N e ranged in these seven regimes between ∼100 and ∼108 viral particles (8). Each regime was independently replicated five times, for a total of 35 experimental lines. Each mixture was kept under the appropriate batch transfer conditions (9) until one of the two variants became fixed. Then, the winner variant that carried a beneficial mutation was isolated. This variant was then placed in head-to-head competition with its nonevolved counterpart (10) to estimate the fitness effect (W) associated with the beneficial mutation that drove it to fixation.

Table 1

Parameters describing the fixation of beneficial mutations under each N e. The number of lines that increased fitness is reported in the second column. The third column shows the number of generations elapsed until fixation of a beneficial mutation (8). A significant correlation between log N e and the time to fixation has been observed (ρS = 0.75, n = 5, one-tailP = 0.0261). The last column shows the probability of fixation by random genetic drift (15).

View this table:

For the smallest population sizes, one can expect genetic drift to play a role in fixing neutral, or even deleterious, mutations. However, previous results have shown that for MARM C clone, the smallest N e used here did not have a considerable deleterious effect (11). Thus, we can safely assume that the fixation of deleterious mutations during our experiment will be minimized by purifying selection.

The estimates obtained for W, under the seven different population-size regimes (12), are shown in Fig. 1. The first prediction we made on the basis of the clonal interference model is completely fulfilled: A significant correlation exists between log N eand the magnitude of the fitness effect (ρS = 0.8929, n = 7, one-tailed P = 0.0034). The larger the population size is (that is, the stronger the clonal interference), the larger the magnitude of the beneficial effect needed to fix a mutation is. As population size increases, there is a shorter waiting time between two consecutive events of beneficial mutation, and thereby more beneficial mutations coexist at a given time.

Figure 1

(left). Influence of effective population size on the magnitude of the fixed fitness effect. Error bars represent standard errors (n = 5). The solid line represents the fit to the theoretical model described in brief in (13) and in more depth in (4).

Each one of the W values used to generate Fig. 1 was transformed into rates of evolution by subtracting from them the fitness of the initial MARM C clone (7) and dividing by the approximate time it took each mutation to become fixed in the population (Table 1). Following (5), we then regressed these rates against N e using (i) a linear model, which implies that the rate of adaptation is proportional to the effective population size, and (ii) a hyperbolic model, which implies that clonal interference will impose a deceleration on the rate of adaptation. These data, as well as the fitting of both models, are shown in Fig. 2. The linear model gave a significant fit to the data [R 2 = 0.6092,F(1,6) = 24.7942, P = 0.0025], indicating that the rate of adaptation increases withN e. However, the hyperbolic model, despite losing a degree of freedom, provides a much better fit to the data than the linear model [partial F test:F(1,5) = 10.1493, P = 0.0244], showing a limit to the rate of adaptation of viral populations imposed by clonal interference. This finding confirms the second prediction we made on the basis of the clonal interference model (4).

Figure 2

(right). Influence of effective population size on the rate of evolution. Error bars represent standard errors (n = 5). Because the rate of adaptation in an asexual population of size zero must be null, the y intercept has been fixed at zero. The dashed line represents the fit of the experimental data to a linear model. The solid line represents the fit to a hyperbolic model [R 2 = 0.8710, F(2,5) = 36.3758,P = 0.0010]. Both curves appear to be exponential because of the common logarithmic scale in the xaxis.

The model of clonal interference (4) has the additional advantage of allowing us to estimate the beneficial mutation rate (μb) as well as the mean selective advantage,E(s), of all beneficial mutations produced in the population (not just those that are fixed). The probability density for the selective coefficients of successful beneficial mutations is given as a part of eq. 6 in (4). From this density, we computed the maximum likelihood estimates and associated variance for μb and E(s) = 1/α (where α is the parameter of an exponential distribution) from the data shown inFig. 1 (13). The solid line drawn in Fig. 1 represents the maximum likelihood fit of the model to the data. The estimated value for the beneficial mutation rate was μb = 6.387 × 10−8 beneficial mutations per genome and generation, with 95% confidence interval in the range 2.57 × 10−8 ≤ μb ≤ 1.58 × 10−7. This value of μb almost warrants that in all seven N e explored, by chance a beneficial mutation will arise, because all lines reached population size of ∼7 × 109 viruses at the end of the daily growth (Nf μb > 1). Also, this value supports our assumption that a single beneficial mutation fixed in each lineage was responsible for the fitness increase. [As an illustration, the probability of generating two beneficial mutations will be Nf μ2 b≪ 1.] Even at the smallest N e, a significant beneficial effect was detected in one of the five replicates (Table 1;W = 2.3692 ± 0.3520,t 2 = 3.8903, one-tail P = 0.0301). In contrast, only at the smallest N edid a line show a significant decline in fitness (W = 0.5286 ± 0.0312, t 2 = −15.1057, one-tail P = 0.0022). Drake et al.(6) estimated the total genomic mutation rate for VSV to be about 3.5 substitutions per genome and generation. Comparing this figure with our estimate of μb, we can infer that about one in 2 × 108 mutations produced in VSV can be considered beneficial. This number is two orders of magnitude smaller than that estimated for E. coli (4). This difference could result from the simpler genome of VSV when compared with E. coli: The more complex a genome is, the more room it has for improvement. Another possible explanation could be the difference in the degree of adaptation of each organisms to its experimental environment: The VSV clones used here have a history in the cell system, which can condition the number of possible beneficial mutations available (5).

The maximum likelihood estimate for E(s) was 0.3062 per day, with 95% confidence interval in the range 0.2371 ≤ E(s) ≤ 0.4309. In other words, the average fitness effect associated with the beneficial mutations produced (not necessarily fixed) is around 31% per day.

The clonal interference model used to estimate μb andE(s) makes two important assumptions. The first is that selection coefficients of beneficial mutations are exponentially distributed. The general shape of this distribution was proposed by Fisher (1), and a statistical argument supporting the use of the exponential was given in (14). The second assumption of the model is to ignore the effect of deleterious mutations. In a small population, deleterious mutations may accumulate through fixation by drift and Muller's ratchet, thus reducing the population's fitness. As stated above, evidence of a reduction in fitness was seen in a line at the smallest N e. To eliminate any possible effect in the computation of μb andE(s), we corrected our data as follows. If the lowest winning fitness of the five replicates was less than one, we assumed that this population had fixed a beneficial mutation of arbitrarily small effect such that the fitness of that population was determined solely by the accumulation of deleterious mutations. The fitness loss through the accumulation of deleterious mutations was assumed to be similar in all five replicates. [This assumption is reasonable because the deleterious mutation rate is quite high (6) such that fitness loss is relatively deterministic.] Thus, the higher fitnesses of the other five replicates (although these may also be less than one) resulted from the fixation of beneficial mutations of significant effect. Although one may question the assumption that the smallest beneficial mutation fixed was of insignificant effect, we assert that the net effect of the assumption is negligible because (i) this correction (and hence this assumption) was only necessary in the smallestN e and (ii) in small N e, clonal interference does not insure the fixation of large-effect mutations; thus, mutations of very small effect may be fixed in these populations just by pure genetic drift. To check this last possibility, we computed the probability of random fixation of an allele in populations with N e equal to those used in our experiment (15). These computations showed that even for our smallest N e, the probability of chance fixation for a neutral allele was <5%. Thus, we can disregard the possibility that most of the mutations fixed during our experiment resulted from genetic drift.

Our results provide evidence that clonal interference occurs in viral populations. This evidence, along with models of clonal interference (1–4), allows certain properties of the adaptive evolution of RNA viruses to be inferred:

1) Adaptive substitutions appear as discrete, rare events, regardless of mutation rate or population size. They often do not occur simply as the result of a single mutation but instead represent the best of several competing mutations. This fact has consequences for the dynamics of drug resistance and the search for resistance mutations.

2) In medium to large populations, the rate of fitness increase is hardly affected by changes in either mutation rate or population size. Some have speculated that the high mutation rates of RNA viruses are maintained evolutionarily because of the great adaptive capacity they bestow (16). When clonal interference is present, however, this argument becomes questionable. When mutation rates are already high, changes in those rates have little effect on the adaptive capacity. Thus, a decrease in mutation rate could come about with little or no decrease in adaptive capacity. Furthermore, a decrease in mutation rate would benefit the population by slowing the accumulation of deleterious mutations. In this light, how is one to explain the high mutation rates of RNA viruses? When clonal interference is considered, it becomes more reasonable to speculate that lower mutation rates would be mechanistically costly or impossible (in terms of complex enzymatic systems required for error detection and correction) for RNA viruses than to suppose that their high mutation rates are maintained for the increased adaptive capacity they confer. If an increased production of beneficial mutations is of little or no advantage to a viral population, then mutator alleles (alleles that confer an elevated mutation rate) must be strictly deleterious because of their increased production of deleterious mutations.

3) Resident populations are protected from invaders simply because of their numerical advantage. A high-fitness VSV clone seeded at low frequency into a resident population of low-fitness variants was displaced by the low-fitness competitors (17). When its initial frequency was above a certain threshold, however, the high-fitness clone always outcompeted the low-fitness variants in the resident population (17). This observation can be easily explained in light of the clonal interference model: If the high-fitness clone is initially present at very low frequencies, it is probable that beneficial mutations arise in the most frequent genotypes, improving their fitness and interfering with the intruder, with the final result of eliminating it from the population. In contrast, when the initial frequency of the high-fitness clone is high, it increases in frequency in the population before the low-fitness variants have a chance to find the right beneficial mutation, resulting in fixation of the high-fitness invader. Such population dynamics in nature might prevent the emergence of dangerous viral pathogens: The existence of a frequency threshold for dominance imposes an element of uncertainty in virus sampling during outbreaks.

  • * Present address: Theoretical Biology and Biophysics Group, Los Alamos National Laboratory, Los Alamos, NM 87545, USA.

  • To whom correspondence should be addressed. E-mail: santiago.elena{at}


View Abstract

Navigate This Article