This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow Supplemental material
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Services
Right arrow Similar articles in this journal
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Right arrowReprints and Permissions
Right arrow Copyright Information
Right arrow Books from ASM Press
Right arrow MicrobeWorld
Citing Articles
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Novella, I. S.
Right arrow Articles by Wilke, C. O.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Novella, I. S.
Right arrow Articles by Wilke, C. O.

 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{triangledown} ,§

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


arrow
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.


arrow
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.


View this table:
[in this window]
[in a new window]

 
TABLE 1. TPS used for the different VSV strains included in this study

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.


Figure 1
View larger version (20K):
[in this window]
[in a new window]

 
FIG. 1. TPS versus initial fitness of the VSV strain. Symbols represent the fitness change after 20 passages: open squares, significant fitness loss; gray circles, no significant fitness change; filled diamonds, significant fitness gain. Data for Bob are moved slightly to the lower right because otherwise they overlap with data for Marilyn. TPS is the CBS for data represented as gray circles, because those correspond to regimens in which there was not a fitness change. We set our significance level to an {alpha} value of 0.05 and corrected for multiple testings using the false-discovery-rate correction (1, 2). The solid line is the best linear fit of initial fitness versus log transmission size for gray circles, and the dashed line represents the same for open squares.

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.


arrow
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.


arrow
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 Back

{triangledown} Published ahead of print on 1 October 2008. Back

§ Supplemental material for this article may be found at http://jvi.asm.org/. Back


arrow
REFERENCES
 
    1
  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. 2
  3. Benjamini, Y., and D. Yekutieli. 2001. The control of the false discovery rate in multiple testing under dependency. Ann. Stat. 29:1165-1188.[CrossRef]
  4. 3
  5. Burch, C. L., and L. Chao. 2000. Evolvability of an RNA virus is determined by its mutational neighbourhood. Nature 406:625-628.[CrossRef][Medline]
  6. 4
  7. 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]
  8. 5
  9. 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]
  10. 6
  11. 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]
  12. 7
  13. 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]
  14. 8
  15. 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]
  16. 9
  17. 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]
  18. 10
  19. 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.





This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow Supplemental material
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Services
Right arrow Similar articles in this journal
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Right arrowReprints and Permissions
Right arrow Copyright Information
Right arrow Books from ASM Press
Right arrow MicrobeWorld
Citing Articles
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Novella, I. S.
Right arrow Articles by Wilke, C. O.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Novella, I. S.
Right arrow Articles by Wilke, C. O.