Virus Replication as a Phenotypic Version of Polynucleotide Evolution

In this paper we revisit and adapt to viral evolution an approach based on the theory of branching process advanced by Demetrius, Schuster and Sigmund ("Polynucleotide evolution and branching processes", Bull. Math. Biol. 46 (1985) 239-262), in their study of polynucleotide evolution. By taking into account beneficial effects we obtain a non-trivial multivariate generalization of their single-type branching process model. Perturbative techniques allows us to obtain analytical asymptotic expressions for the main global parameters of the model which lead to the following rigorous results: (i) a new criterion for "no sure extinction", (ii) a generalization and proof, for this particular class of models, of the lethal mutagenesis criterion proposed by Bull, Sanju\'an and Wilke ("Theory of lethal mutagenesis for viruses", J. Virology 18 (2007) 2930-2939), (iii) a new proposal for the notion of relaxation time with a quantitative prescription for its evaluation, (iv) the quantitative description of the evolution of the expected values in in four distinct "stages": extinction threshold, lethal mutagenesis, stationary "equilibrium" and transient. Finally, based on these quantitative results we are able to draw some qualitative conclusions.Comment: 23 pages, 1 figure, 2 tables. arXiv admin note: substantial text overlap with arXiv:1110.336

A Multi-Step Process of Viral Adaptation to a Mutagenic Nucleoside Analogue by Modulation of Transition Types Leads to Extinction-Escape

Resistance of viruses to mutagenic agents is an important problem for the development of lethal mutagenesis as an antiviral strategy. Previous studies with RNA viruses have documented that resistance to the mutagenic nucleoside analogue ribavirin (1-β-D-ribofuranosyl-1-H-1,2,4-triazole-3-carboxamide) is mediated by amino acid substitutions in the viral polymerase that either increase the general template copying fidelity of the enzyme or decrease the incorporation of ribavirin into RNA. Here we describe experiments that show that replication of the important picornavirus pathogen foot-and-mouth disease virus (FMDV) in the presence of increasing concentrations of ribavirin results in the sequential incorporation of three amino acid substitutions (M296I, P44S and P169S) in the viral polymerase (3D). The main biological effect of these substitutions is to attenuate the consequences of the mutagenic activity of ribavirin —by avoiding the biased repertoire of transition mutations produced by this purine analogue—and to maintain the replicative fitness of the virus which is able to escape extinction by ribavirin. This is achieved through alteration of the pairing behavior of ribavirin-triphosphate (RTP), as evidenced by in vitro polymerization assays with purified mutant 3Ds. Comparison of the three-dimensional structure of wild type and mutant polymerases suggests that the amino acid substitutions alter the position of the template RNA in the entry channel of the enzyme, thereby affecting nucleotide recognition. The results provide evidence of a new mechanism of resistance to a mutagenic nucleoside analogue which allows the virus to maintain a balance among mutation types introduced into progeny genomes during replication under strong mutagenic pressure

Does Mutational Robustness Inhibit Extinction by Lethal Mutagenesis in Viral Populations?

Lethal mutagenesis is a promising new antiviral therapy that kills a virus by raising its mutation rate. One potential shortcoming of lethal mutagenesis is that viruses may resist the treatment by evolving genomes with increased robustness to mutations. Here, we investigate to what extent mutational robustness can inhibit extinction by lethal mutagenesis in viruses, using both simple toy models and more biophysically realistic models based on RNA secondary-structure folding. We show that although the evolution of greater robustness may be promoted by increasing the mutation rate of a viral population, such evolution is unlikely to greatly increase the mutation rate required for certain extinction. Using an analytic multi-type branching process model, we investigate whether the evolution of robustness can be relevant on the time scales on which extinction takes place. We find that the evolution of robustness matters only when initial viral population sizes are small and deleterious mutation rates are only slightly above the level at which extinction can occur. The stochastic calculations are in good agreement with simulations of self-replicating RNA sequences that have to fold into a specific secondary structure to reproduce. We conclude that the evolution of mutational robustness is in most cases unlikely to prevent the extinction of viruses by lethal mutagenesis

Mutagenesis-Mediated Virus Extinction: Virus-Dependent Effect of Viral Load on Sensitivity to Lethal Defection

Background: Lethal mutagenesis is a transition towards virus extinction mediated by enhanced mutation rates during viral genome replication, and it is currently under investigation as a potential new antiviral strategy. Viral load and virus fitness are known to influence virus extinction. Here we examine the effect or the multiplicity of infection (MOI) on progeny production of several RNA viruses under enhanced mutagenesis. Results: The effect of the mutagenic base analogue 5-fluorouracil (FU) on the replication of the arenavirus lymphocytic choriomeningitis virus (LCMV) can result either in inhibition of progeny production and virus extinction in infections carried out at low multiplicity of infection (MOI), or in a moderate titer decrease without extinction at high MOI. The effect of the MOI is similar for LCMV and vesicular stomatitis virus (VSV), but minimal or absent for the picornaviruses foot-and-mouth disease virus (FMDV) and encephalomyocarditis virus (EMCV). The increase in mutation frequency and Shannon entropy (mutant spectrum complexity) as a result of virus passage in the presence of FU was more accentuated at low MOI for LCMV and VSV, and at high MOI for FMDV and EMCV. We present an extension of the lethal defection model that agrees with the experimental results. Conclusions: (i) Low infecting load favoured the extinction of negative strand viruses, LCMV or VSV, with an increase of mutant spectrum complexity. (ii) This behaviour is not observed in RNA positive strand viruses, FMDV or EMCV. (iii) The accumulation of defector genomes may underlie the MOI-dependent behaviour. (iv) LCMV coinfections are allowed but superinfection is strongly restricted in BHK-21 cells. (v) The dissimilar effects of the MOI on the efficiency of mutagenic-based extinction of different RNA viruses can have implications for the design of antiviral protocols based on lethal mutagenesis, presently under development. © 2012 Moreno et al.Centro de Biología Molecular Severo Ochoa; Ministerio de Ciencia e Innovación (MICINN); Fundación Ramón ArecesPeer Reviewe

Using enhanced number and brightness to measure protein oligomerization dynamics in live cells

Protein dimerization and oligomerization are essential to most cellular functions, yet measurement of the size of these oligomers in live cells, especially when their size changes over time and space, remains a challenge. A commonly used approach for studying protein aggregates in cells is number and brightness (N&B), a fluorescence microscopy method that is capable of measuring the apparent average number of molecules and their oligomerization (brightness) in each pixel from a series of fluorescence microscopy images. We have recently expanded this approach in order to allow resampling of the raw data to resolve the statistical weighting of coexisting species within each pixel. This feature makes enhanced N&B (eN&B) optimal for capturing the temporal aspects of protein oligomerization when a distribution of oligomers shifts toward a larger central size over time. In this protocol, we demonstrate the application of eN&B by quantifying receptor clustering dynamics using electron-multiplying charge-coupled device (EMCCD)-based total internal reflection microscopy (TIRF) imaging. TIRF provides a superior signal-to-noise ratio, but we also provide guidelines for implementing eN&B in confocal microscopes. For each time point, eN&B requires the acquisition of 200 frames, and it takes a few seconds up to 2 min to complete a single time point. We provide an eN&B (and standard N&B) MATLAB software package amenable to any standard confocal or TIRF microscope. The software requires a high-RAM computer (64 Gb) to run and includes a photobleaching detrending algorithm, which allows extension of the live imaging for more than an hour

Composition and Interactions of Hepatitis B Virus Quasispecies Defined the Virological Response During Telbivudine Therapy

Reverse transcriptase (RT) mutations contribute to hepatitis B virus resistance during antiviral therapy with nucleos(t)ide analogs. However, the composition of the RT quasispecies and their interactions during antiviral treatment have not yet been thoroughly defined. In this report, 10 patients from each of 3 different virological response groups, i.e., complete virological response, partial virological response and virological breakthrough, were selected from a multicenter trial of Telbivudine treatment. Variations in the drug resistance-related critical RT regions in 107 serial serum samples from the 30 patients were examined by ultra-deep sequencing. A total of 496,577 sequence reads were obtained, with an average sequencing coverage of 4,641X per sample. The phylogenies of the quasispecies revealed the independent origins of two critical quasispecies, i.e., the rtA181T and rtM204I mutants. Data analyses and theoretical modeling showed a cooperative-competitive interplay among the quasispecies. In particular, rtM204I mutants compete against other quasispecies, which eventually leads to virological breakthrough. However, in the absence of rtM204I mutants, synergistic growth of the drug-resistant rtA181T mutants with the wild-type quasispecies could driveNational Science and Technology Major Project [2012ZX10002003, 2013ZX10004104-002]; key Clinical Specialty Discipline Construction Program, Shanghai Commission for Science and Technology [12ZR1421200]; National Natural Science Foundation of China [91331107, 91331109, 31171279]SCI(E)[email protected]; [email protected]; [email protected]