Previous Article | Next Article 
Journal of Virology, December 2008, p. 12589-12590, Vol. 82, No. 24
0022-538X/08/$08.00+0 doi:10.1128/JVI.01394-08
Copyright © 2008, American Society for Microbiology. All Rights Reserved.
A Linear Relationship between Fitness and the Logarithm of the Critical Bottleneck Size in Vesicular Stomatitis Virus Populations
,
Isabel S. Novella,1*
Ranendra N. Dutta,1 and
Claus O. Wilke2
Department of Medical Microbiology and Immunology, College of Medicine, University of Toledo Health Science Campus, Toledo, Ohio 43614,1
Center for Computational Biology and Bioinformatics, Section of Integrative Biology, and Institute for Cell and Molecular Biology, University of Texas at Austin, Austin, Texas 787122
Received 3 July 2008/
Accepted 18 September 2008

ABSTRACT
We explored the relationship between fitness change and population
size during transmission in vesicular stomatitis populations
of very high fitness. The results show a linear correlation
between the logarithm of the critical bottleneck size (population
size at which there are no significant fitness changes after
20 passages) and the initial fitness of the population. In addition,
limits to fitness increases during large-population passages
of very-high-fitness strains were abolished by increasing the
population size during transmission, indicating that beneficial
variation is still available in these populations.

TEXT
Viral evolution is the result of selection and random drift
operating on the genetic variation that arises during replication.
Random drift dominates during bottleneck transmissions because
sampling effects lead to the loss of beneficial mutations or
the fixation of deleterious mutations. Bottlenecks are common
during natural infections, particularly among respiratory viruses.
In an earlier report, we found a correlation between the fitness
of vesicular stomatitis virus (VSV) strains and the critical
bottleneck size (CBS) in low-fitness populations subjected to
repeated severe bottlenecks (transmission sizes of 2 to 30 PFU)
(
8). We defined the CBS as the number of virus particles used
for transmission that resulted in no net fitness changes at
the end of 20 passages. Later, we and others reported that limits
to fitness gains (and fitness losses) during large-population
passages were consistent with a linear relationship between
initial fitness and the logarithm of the CBS (
9,
10). Miralles
et al. (
6) argued against this relationship and proposed that
limits to further fitness gains reflected the exhaustion of
beneficial variation. Their main criticism was that a linear
relationship observed at small bottleneck sizes could not be
extrapolated to bottleneck sizes that were several orders of
magnitude larger. To resolve this question, we tested whether
the relationship between the logarithm of the CBS and fitness
was still linear for strains with high fitness.
We generated several VSV populations with fitness levels between 5.5 and 15 in BHK-21 cells (Table 1) through 35 to 50 large-population passages (2 x 105 PFU/passage unless otherwise indicated) of monoclonal antibody-resistant mutant (MARM) U, a strain that has neutral fitness compared to the wild type (wt) and differs from the wt only in a genetic marker that provides resistance to monoclonal antibody I1 (5). Each strain was obtained from an independent replica of MARM U passages except strain Marilyn, which is the progeny of Bonnie after 20 additional passages at 2 x 104 PFU/passage, and Victoria, which is the progeny of Marilyn after 20 additional passages at 4 x 106 PFU/passage. All passages were done at a constant multiplicity of infection of 0.1 PFU/cell by using flasks of the appropriate size for each population size; for all passages, the number of cells infected equals the number of PFU used for infection. For each population, we carried out six replicas of 20 passages at the expected CBS, and for some of them, we also carried out passages at transmission sizes below or above the expected CBS (Table 1). We determined the fitness of the progeny for each replica in competitions between the test MARM virus and reference wt as previously described (5). For each competition, we mixed the test mutant and wt (2 x 105 PFU), and we infected a monolayer of BHK-21 cells at a multiplicity of infection of 0.1 and incubated the flask for 20 to 24 h at 37°C. We then determined the MARM/wt ratios before (R0) and after (R1) competition by using triplicate plaque assays in the presence and absence of monoclonal antibody I1, and we calculated fitness by dividing R1 by R0. Each fitness value was calculated as the average of three to six independent determinations.
Our results showed that the correlation between the logarithm
of the CBS and fitness was linear within the range of initial
fitness from previous work and for the strains reported here
(Fig.
1) (for a complete data set, see Table S1 in the supplementary
material). Furthermore, as predicted, passages at population
sizes below the CBS consistently led to fitness losses, with
two exceptions: the fitness loss in strain Marilyn was not significant
after correction for multiple testing (
P = 0.09), and strain
Bob gained fitness at and below the predicted CBS. These results
confirmed our suggestion that sampling effects influence virus
evolution at large-population transmissions if viral fitness
is sufficiently high. Furthermore, Marilyn had significant fitness
gains at a transmission size 20-fold larger than the CBS. Thus,
beneficial variation was still present even for strains such
as the high-fitness strain.
In this work, we defined the CBS based on the (lack of) fitness
change after 20 passages. We have previously shown that the
fitness change during the first 20 passages was dominated by
the sampling of preexisting variation (
4). Therefore, the fitness
increases we observed were conceivably due to beneficial mutations
carried over from the previous passages. However, it is also
possible that these mutations were generated de novo during
replication and sampled when the population size at transmission
was increased sufficiently. Finally, more-complex interactions
among the larger number of variants during replication at higher
transmission population sizes (TPS) may contribute to the overall
fitness increases.
We do not have a solid explanation for why strain Bob behaves differently from the rest of the strains. One possibility is that Bob is located on a different fitness peak than the other strains (3), and its beneficial mutation rate is increased. Bonnie, Marilyn, and Victoria are a continuation of the same replica of passages, and it is reasonable to propose that they are all in different sites of the same fitness peak. In contrast, of all the high-fitness populations, only Marco and Bob were generated independently. Nevertheless, the predictions of the model were generally good for all the strains in the data set, suggesting that displacement into different fitness peaks is probably infrequent.

ACKNOWLEDGMENTS
We thank Arthur Chan for technical support and Douglas Lyles
(Wake Forest University) for the I1 hybridoma.
This work was supported by NIH grant R01 AI065960. C.O.W. was also supported by a Reeder Centennial Fellowship in systematic and evolutionary biology.

FOOTNOTES
* Corresponding author. Mailing address: Department of Microbiology and Immunology, University of Toledo Health Science Campus, Mail stop 1210, Toledo, OH 43614. Phone: (419) 383-6442. Fax: (419) 383-4221. E-mail:
isabel.novella{at}utoledo.edu 
Published ahead of print on 1 October 2008. 
Supplemental material for this article may be found at http://jvi.asm.org/. 

REFERENCES
1 - Benjamini, Y., and Y. Hochberg. 1995. Controlling the false-discovery rate: a practical and powerful approach to multiple testing. J. R. Stat. Soc. B 57:289-300.
2 - Benjamini, Y., and D. Yekutieli. 2001. The control of the false discovery rate in multiple testing under dependency. Ann. Stat. 29:1165-1188.[CrossRef]
3 - Burch, C. L., and L. Chao. 2000. Evolvability of an RNA virus is determined by its mutational neighbourhood. Nature 406:625-628.[CrossRef][Medline]
4 - Dutta, R. N., I. M. Rouzine, S. D. Smith, C. O. Wilke, and I. S. Novella. 2008. Rapid adaptive amplification of preexisting variation in an RNA virus. J. Virol. 82:4354-4362.[Abstract/Free Full Text]
5 - Holland, J. J., J. C. de la Torre, D. K. Clarke, and E. Duarte. 1991. Quantitation of relative fitness and great adaptability of clonal populations of RNA viruses. J. Virol. 65:2960-2967.[Abstract/Free Full Text]
6 - Miralles, R., A. Moya, and S. F. Elena. 2000. Diminishing returns of population size in the rate of RNA virus adaptation. J. Virol. 74:3566-3571.[Abstract/Free Full Text]
7 - Novella, I. S., and B. E. Ebendick-Corpus. 2004. Molecular basis of fitness loss and fitness recovery in vesicular stomatitis virus. J. Mol. Biol. 342:1423-1430.[CrossRef][Medline]
8 - Novella, I. S., S. F. Elena, A. Moya, E. Domingo, and J. J. Holland. 1995. Size of genetic bottlenecks leading to virus fitness loss is determined by mean initial population fitness. J. Virol. 69:2869-2872.[Abstract]
9 - Novella, I. S., J. Quer, E. Domingo, and J. J. Holland. 1999. Exponential fitness gains of RNA virus populations are limited by bottleneck effects. J. Virol. 73:1668-1671.[Abstract/Free Full Text]
10 - Rouzine, I. M., J. Wakeley, and J. M. Coffin. 2003. The solitary wave of asexual evolution. Proc. Natl. Acad. Sci. USA 100:587-592.[Abstract/Free Full Text]
Journal of Virology, December 2008, p. 12589-12590, Vol. 82, No. 24
0022-538X/08/$08.00+0 doi:10.1128/JVI.01394-08
Copyright © 2008, American Society for Microbiology. All Rights Reserved.