- Research article
- Open Access
Topology of evolving, mutagenized viral populations: quasispecies expansion, compression, and operation of negative selection
© Ojosnegros et al; licensee BioMed Central Ltd. 2008
- Received: 29 February 2008
- Accepted: 17 July 2008
- Published: 17 July 2008
The molecular events and evolutionary forces underlying lethal mutagenesis of virus (or virus extinction through an excess of mutations) are not well understood. Here we apply for the first time phylogenetic methods and Partition Analysis of Quasispecies (PAQ) to monitor genetic distances and intra-population structures of mutant spectra of foot-and-mouth disease virus (FMDV) quasispecies subjected to mutagenesis by base and nucleoside analogues.
Phylogenetic and PAQ analyses have revealed a highly dynamic variation of intrapopulation diversity of FMDV quasispecies. The population diversity first suffers striking expansions in the presence of mutagens and then compressions either when the presence of the mutagenic analogue was discontinued or when a mutation that decreased sensitivity to a mutagen was selected. The pattern of mutations found in the populations was in agreement with the behavior of the corresponding nucleotide analogues with FMDV in vitro. Mutations accumulated at preferred genomic sites, and dn/ds ratios indicate the operation of negative (or purifying) selection in populations subjected to mutagenesis. No evidence of unusually elevated genetic distances has been obtained for FMDV populations approaching extinction.
Phylogenetic and PAQ analysis provide adequate procedures to describe the evolution of viral sequences subjected to lethal mutagenesis. These methods define the changes of intra-population structure more precisely than mutation frequencies and Shannon entropies. PAQ is very sensitive to variations of intrapopulation genetic distances. Strong negative (or purifying) selection operates in FMDV populations subjected to enhanced mutagenesis. The quantifications provide evidence that extinction does not imply unusual increases of intrapopulation complexity, in support of the lethal defection model of virus extinction.
- Shannon Entropy
- Viral Population
- Mutant Spectrum
- Equine Infectious Anemia Virus
- Average Genetic Distance
RNA viruses replicate as mutant distributions termed viral quasispecies. This is a consequence of high mutation rates operating during RNA genome copying, due to the absence of proofreading-repair activities in the relevant RNA-dependent RNA polymerases and RNA-dependent DNA polymerases [1, 2]. Most phylogenetic relationships among RNA viruses have been established using the consensus (or population) genomic sequences that represent a weighted average of multiple, closely related sequences present at each time point, in each virus sample obtained for analysis . Phylogenetic relationships established with consensus viral sequences have been instrumental to classify viruses and to determine origin of emergent viruses and rates of virus evolution [1, 4, 5]. For many purposes it is important to analyze phylogenetically the relationship among different genomes from the same mutant spectrum of a viral quasispecies. This type of analysis may reveal the existence of genome subpopulations within mutant spectra that might encode different phenotypic traits. Also, it allows the calculation of average genetic distances among individual components of the mutant spectrum, a parameter that can be a predictor of biological behaviour . As an example, a study with a poliovirus mutant which displays a -3-to 5-fold higher template-copying fidelity than the wild type documented that a narrow mutant spectrum impeded the virus to reach the brain of susceptible mice and produce neuropathology [6, 7]. An early study documented that complexity of the coronavirus murine hepatitis virus quasispecies influenced the pathogenic potential of this virus for mice . A broad hepatitis C (HCV) virus mutant spectrum was associated with poor response to treatment by ribavirin and interferon α , and rapid early evolution of the virus led to a chronic infection . Some studies have found an association between a reduction of mutant spectrum complexity of HCV at early stages of treatment and viral clearance (; reviewed in ). Recently, the role of the mutant spectrum in adaptation of West Nile virus has been documented [13, 14]. Therefore, there is a need to develop methods to describe the relationship among components of mutant spectra in viral populations.
An estimate of the complexity and internal relationships among genomes within a mutant spectrum can be obtained through application of distance-based phylogenetic methods, such as neighbour-joining (NJ)  or maximum likelihood  procedures. However, the reliability of the derived genome clusters may be questionable when the number of mutations distinguishing different components of a mutant spectrum is small. Small genetic distances among genomes of a mutant spectrum are found when a viral clone (single genome) has undergone a limited number of passages in cell culture . Despite this limitation, the general topology of a NJ tree may be robust enough to provide information about the evolutionary pattern undergone by the viral population.
An alternative method developed to group closely related sequences is Partition Analysis of Quasispecies (PAQ), which is considered a non-hierarchical clustering method . PAQ groups together the viral sequences separated by the shortest genetic distances. For this purpose a centre genotype that nucleates a group of sequences within a circle (cluster) with a previously set radius is selected. High compactness of a group is defined by having a larger number of variants near the centre of the circle than at its boundary. In this fashion, multiple, overlapping or non-overlapping groups can be defined within mutant spectra of viral quasispecies. PAQ has been successfully used to analyze subpopulations of equine infectious anemia virus , Ty1-copia retrotransposon sequences , regulatory sequences of bovine viral diarrhoea virus , and sequential isolates of hepatitis A virus . PAQ should find increasing application to describe viral quasispecies in view of its power to quantify relationships among related sequences, and of the critical role of mutant spectra in virus behaviour.
Our laboratory has studied quasispecies dynamics and evolution of foot-and-mouth disease virus (FMDV), a widespread picornavirus pathogen that causes the economically most important animal viral disease (reviews in [23, 24]). The molecular epidemiology of FMDV has been analyzed in detail by establishing an extensive database of genomic nucleotide sequences, as well as by defining phylogenetic relationships among and within serotypes and subtypes of the virus [25–27]. FMDV has been used as a model system to investigate lethal mutagenesis, a term coined to describe virus extinction through an increase of viral mutation rate during replication (; reviews in [29–31]). Studies on the interference that mutated FMDV and lymphocytic choriomeningitis virus (LCMV) populations exert on the infectivity of the corresponding standard genomes [32–35] led to the proposal of the lethal defection model of virus extinction . According to this model, extinction can occur with a modest average number of mutations per genome that generates a class of genomes termed defectors, that can replicate but interfere with replication of the standard genomes [34, 36]. Defector genomes differ from defective-interfering (DI) particles described for many viruses  in that DI particles are dependent on helper virus for replication while defectors are replication-competent [34, 35].
In the course of the studies on lethal mutagenesis of FMDV, biological clones of the virus were subjected to serial cytolytic passages in cell culture in the absence or presence of the mutagenic bases or nucleoside analogues 5-fluorouracil (FU), ribavirin (1-β-D ribofuranosyl-1, 2, 3- triazole- 3 carboxamide) (R), or 5-azacytidine (AZC) [32, 38–41]. These treatments have provided a number of viral populations differing in their mutant spectrum complexity, as measured by the average mutation frequency and Shannon entropy within mutant spectra [38–42]. However, possible alterations of phylogenetic relationships among components of the mutant spectra of mutagenized quasispecies have not been studied. To characterize possible changes among components of FMDV mutant spectra as a result of the mutagenic treatments, here we compare phylogenetic and PAQ analyses of several multiply-passaged FMDV populations, and their mutagenized counterparts. The results show that phylogenetic and clustering methods can provide a quantification of the mutagenic activity exerted on viral populations by reflecting changes in genetic distances among components of the mutant spectra of the viral quasispecies. PAQ analysis describes mutant distributions as spheres with size which is proportional to the genetic diversity, which varies as a result of enhanced mutagenesis. This type of representation has revealed that viral populations can respond rapidly to environmental changes, with striking switches between relaxation and compactness of the population diversity, that were not apparent from the comparison of mutation frequencies or Shannon entropies. Evidence is presented that identical or closely related consensus sequences may hide different subpopulations of genomes bearing distinct relationships among them.
Cells, viruses, and infections
The origin of BHK-21 cells, and procedures for cell growth, and infection of cell monolayers with foot-and-mouth disease virus (FMDV) have been previously described [39–41, 43, 44]. The following FMDV clonal populations have been used in the studies on lethal mutagenesis: i) FMDV C-S8c1, a plaque-purified derivative of the natural isolate C1 Santa-Pau Spain 70 , a representative of European serotype C FMDV , GenBank accession number NC_002554; ii) C-S8c1 multiply passaged at high multiplicity of infection (m.o.i.) in BHK-21 cells; the viral populations derived from each passage are indicated with a "p" before the passage number (i.e. C-S8p50 means C-S8c1 passaged 50 times in BHK-21 cells). iii) FMDV MARLS, a monoclonal antibody (mAb)-escape mutant, selected from population C-S8c1p213 , GeneBank accession number AF274010. iv) C-S8c1p260p3d, the viral population resulting from three diluted serial infections (mo.i. = 0.1) of C-S8c1p260 in BHK-21 cells ; gene bank accession number DQ409185; v) C10280, a clone derived from subjecting C-S8c1 to 280 plaque-to-plaque transfers in BHK-21 cells [47, 48].
Extraction of RNA, cDNA synthesis, PCR amplification and nucleotide sequencing
Procedures for extraction of RNA from the supernatants of infected cell cultures, reverse transcription (RT) of FMDV RNA and PCR amplification with high fidelity Pfu DNA polymerase have been previously described [32, 39, 40, 48]. Molecular cloning in plasmids pET-28a3Dpol or pMT-28, as well as controls to ensure that molecular clones reflect accurately the composition of the quasispecies under analysis, were described in the primary publications that produced the genomic sequences analyzed here [46, 49–52]. Nucleotide sequencing was performed using the Big Dye Terminator Cycle Sequencing Kit (Abi Prism; Applied Biosystems) and the automated sequencer ABI 373 and ABI 3730; all sequences were determined at least in duplicate, from independent sequencing reactions.
Mutagenic agents and treatments
Treatment with ribavirin
Treatment with fluorouracil, azacytidine and guanidine hydrochloride
5-fluorouracil (FU) and azacytidine (AZC) were used as mutagenic base analogues (Figure 1), while guanidine hydrochloride (G) was used as inhibitor of FMDV replication, as previously described [38, 39, 41]. FU medium contained 200 μg/ml FU, and FUG medium contained 200 μg/ml FU and 4 mM G. Solutions were sterilized by filtration and stored at 4°C for a maximum of 15 days before use. Infections in the presence of FUG were performed as previously described [38, 39, 41].
Multiple alignments of consensus sequences and molecular clones were carried out with the program CLUSTAL W , inserted into Bioedit package , using FMDV C-S8c1 and/or FMDV MARLS (GenBank accession numbers AJ133357 and AF274010, respectively) as the reference sequences. Pairwise distances matrix was generated by the program MEGA3.1 , using the Kimura-2 parameter model . Tree topology was inferred by three phylogenetic methods: (i) Neighbor-Joining (NJ)  using also the MEGA3.1 package ; bootstrap re-sampling (1000 data sets) of the multiple alignments was used to test the statistical robustness of the trees obtained by NJ . (ii) Maximum Parsimony (MP) (using the program DNAPARS from PHYLIP v3.5 package) . (iii) Maximum Likelihood (ML) trees were generated by the program PUZZLE  using the Tamura-Nei substitution model , and the Gamma distributed rates with eight parameters (TN-8Γ) as heterogeneity model (see additional file 1). Additionally, for the additional file 2, a ML analysis of mutagenized populations of FMDV was performed by means of the program Modelgenerator , useful to identify the optimal evolutionary model (Akaike Information Criteria and Hierarchical Likelihood Ratio Test indicated that the GTR model best fit the sequence data). Using this model, ML trees were constructed using software from the PhyML program , available at . As a measure of the robustness of each node, we used an approximate Likelihood Ratio Test (aLRT), which assesses that the branch being studied provides a significant likelihood gain, in comparison with the null hypothesis that involves collapsing the branch under study but leaving the rest of the tree topology unaffected .
Characterization of mutant spectra
in which p i is the proportion of each sequence of the mutant spectrum, and N is the total number of sequences compared.
in which n is the number of neighbors within the group with centre variant c, and D ic is the genetic distance (Hamming) between variants i and c. The populations are represented as spheres with the diameter proportional to the AD value. Standard errors have also been calculated (see additional file 3: Standard errors).
The nonsynonymous mutations corrected per nonsynonymous site (dn) and the synonymous mutations corrected per synonymous site (ds) were calculated using SNAP [66, 67]. Kn and Ks are the ratio of nonsynonymous and synonymous mutations respectively, per nucleotide (without any correction). The sliding-window software K-Estimator  was used to infer the Ks and Kn in the FMDV genome using 200 nucleotide windows and a shift of 100 nucleotides per step, to cover all the polyprotein-coding region. The software calculates the confidence interval using Monte Carlo simulations.
Software for phylogenetic and sequence analyses were retrieved from the corresponding references listed in the different sections of Methods.
One Way ANOVA, Tukey post hoc tests for non equal n, and standard errors were calculated using Statistica 6.0 software package (StatSoft 2001).
Phylogenetic evaluation of the mutagenesis-induced diversification of viral populations
The present study was aimed at analyzing retrospectively sets of consensus and clonal nucleotide sequences determined in our laboratory from FMDV populations subjected to serial cytolytic infections in BHK-21 cells, in the absence or presence of the mutagenic nucleotide analogues FU, R or AZC [32, 33, 39, 40] – (Figure 1).
Number and type of mutations scored, and resulting mutation frequencies, calculated for the mutant spectra of the FMDV populations analyzed in the present study
Number of clones
Total nt. sequenced
Minimum mutation frequency (×10-3)f
Maximum mutation frequency (×10-3)g
No Drug p1
[10 μg/ml] AZC
[200 μg/ml] FU
No drug p25
[200 μg/ml] FU
[800 μM] R, No drug
[500–800 μM] R
[800 μM] R
[800–5000 μM] R
Interestingly, the population that underwent 30 passages in the presence of R and then 5 passages in the absence of R, showed a more compact topology, and subclusters of clones were distinguished (RA0p35* in Figure 3B). It must be noted that, due to the close relatedness among the nucleotide sequences analyzed (see Background and Methods), the bootstrap values associated with the derived NJ trees suggest limited robustness of the derived clusterings. Nevertheless, the same general topological features are also observed when MP and ML were used (see additional file 1: neighbour-joining, maximum likelihood and maximum parsimony analysis of populations RAp35 and RA0p35) (see also statistical evaluation of PAQ, below).
The control population CAp35, obtained after 35 serial cytolytic infections of clone MARLS in cell culture in the absence of mutagen, shows a characteristic tree with evident lower diversity than the populations treated with mutagen (Figure 3), as documented by the very low mutation frequency value (12-fold smaller than population RAp35, see Table 1), and the average shortness of all the branches.
Partition analysis of quasispecies (PAQ) applied to the clonal analysis of mutagenized FMDV populations
PAQ analysis of the populations derived from clone C-S8c1 subjected to 25 passages either in the presence of FU or AZC or in the absence of mutagen, showed a 1.5-fold expansion of intrapopulation diversity as a result of virus propagation in the absence of mutagen, and 13.50-fold and 9.24-fold expansions as a result of passages in the presence of FU or AZC, respectively (Figure 4A). These values illustrate again the higher mutagenic potential of FU than AZC (Tukey post hoc test, p = 0.004) on FMDV replicating in BHK-21 cells .
A dramatic effect of R treatment was unveiled by PAQ through the comparison of MARLS populations passaged in the presence or absence of R (Figure 4B). Remarkably, multiple passages in the presence of up to 800 μM R, resulted in populations with significant differences in AD values by means of analysis of variance (One-way ANOVA: F = 85.5, df = 65, p < 0.001). RAp35 showed a 15.47-fold larger AD value (Tukey post hoc test, p < 0.001) than the population CAp35, the parallel control in the absence of the drug (Figure 4B).
The 3D of RAp35 and its derived populations harbours the substitution M296I which decreases its sensitivity to R . Further passages in the presence of 800 μM R led to a reduction of diversity of population RAp45 (1.90-fold reduction in AD value, Tukey post hoc test, p < 0.001; Table 1). Once the diversity decreased in RAp45, it did not further increase significantly after 15 additional passages under higher (up to 5000 μM) concentrations of R (RAp60) (1.10 fold, Tukey post hoc test, p = 0.9).
Population RAp30 (depicted in Figure 1B) was subjected in parallel to 5 further passages either in the absence (RA0p35) or the presence (RAp35) of R. The intrapopulation diversity in RA0p35 was 2.15-fold smaller than RAp35 (Tukey post hoc test, p < 0.001) (see Discussion).
The modification of the relationship among components of the mutant spectrum is not reflected in Shannon entropy (which is saturated in the highest value, 1, for all the populations from the RA lineage tested) (see Table 1). Also maximum mutation frequency varies within a narrow range of 1.7-fold to 2.4-fold when the different populations are compared with RAp35 (table 1). There are no statistically significant differences in the minimum/maximum mutation frequency of population RAp35 and RA0p35 (one-way ANOVA: F = 2.57, df = 1.33, p = 0.118). This ratio is larger (higher diversity) in populations RAp45 and RAp60, than in population RAp35 (Tukey post-hoc test, p = 0.0073 and p = 0.013, respectively), despite having lower AD values and, therefore lower intrapopulation diversity. These comparisons support that the AD values derived from PAQ are a more sensitive and realistic descriptor of intrapopulation diversity than mutation frequencies.
Mutational bias and rate of fixation of mutations
Evidence of purifying selection during enhanced mutagenesis
Corrected ratio of nonsynonymous to synonymous substitution in several FMDV populations
dn/ds (Molecular clones)b
dn/ds (Consensus sequences)c
Virus extinction without a large excess of mutations
Several remarkable modifications in the mutant spectrum occur when viral populations replicate in the presence of mutagenic agents such as base or nucleoside analogues. The initial experiments by J.J. Holland and colleagues documented very modest (1.1- to 2.8-fold) increases of mutation frequency at single base sites of poliovirus and vesicular stomatitis virus, associated with severe decrease of infectivity, as a consequence of enhanced mutagenesis by a variety of mutagenic agents . Subsequent work showed that replication of different viruses in the presence of mutagenic base or nucleoside analogues could lead to virus extinction accompanied of very modest (from barely measurable up to 6.4-fold) increases of mutation frequency and mutant spectrum complexity ([28, 39, 74]; reviews in [29, 30, 71]). In our previous studies with FMDV and lymphocytic choriomeningitis virus we have observed that in the course of the transition of the viruses towards extinction: i) there is a 102- to 103-fold decrease in specific infectivity (number of infectious units divided by the total amount of viral RNA); ii) low viral loads and low replication capacity (fitness) favour extinction; and iii) internal, interfering interactions among mutagenized components of the mutant spectra play an important role in the extinction [32–36, 38, 39, 75].
The critical participation of the mutant spectrum in viral extinction led to the proposal of "lethal defection" as a model of virus extinction by lethal mutagenesis [34, 35], and motivated the phylogenetic and partition analyses of mutant spectra reported in the present study. The results support our previous conclusions on the general increase of mutant spectrum complexity -here documented by average genetic distances in phylogenetic trees and quantified by the AD parameter in PAQ- as a result of the various mutagenic treatments. Moreover, new insights into the events associated with enhanced mutagenesis have been unveiled by the comparison of phylogenetic and PAQ analyses. The NJ tree topology of a population evolved under enhanced mutagenesis conditions consists of a basal consensus sequence, with the individual components of the population radiating from it (the more mutated the population the bigger the radiation of the branches) and with no apparent population structure in the form of subclusters. Trees constructed with FMDV populations whose 3D (polymerase) included a substitution that decreased the sensitivity to R (RAp45 and RAp60)  displayed a slight reduction of branch length, but maintained the same general topology. This is in sharp contrast with population RA0p35, which, after only 5 passages in the absence of R, presents an internal structure completely rearranged, as shown by the NJ tree, and confirmed with ML and MP methods (Figure 3B, and additional file 1: neighbour-joining, maximum likelihood and maximum parsimony analysis of populations RAp35 and RA0p35).
The average distance parameter AD of PAQ reflects the fine detail of the intra population structure of diversity better than the other estimators traditionally used, such as mutation frequencies or Shannon entropy. Mutation frequency does not capture the distribution of mutations among the individual components of a mutant spectrum, and the Shannon entropy reaches the maximum possible value of 1 when all the sequences under study are different, independently of the mutational load. This lack of resolution is circumvented by PAQ due to the pairwise comparison of all the clones with respect to the central sequence. In this sense, any subtle variation in the spatial distribution of mutation frequency may be amplified yielding a more sensitive, accurate and non-saturating description of the internal structure of the population. The present analysis has documented striking expansions of the average intrapopulation genetic distance associated with mutagenic treatments (Figure 4). Interestingly the analysis of the RA lineage revealed that, after an initial expansion of intrapopulation diversity, a remarkably fast contraction of the average distance occurred in two situations. First, a 2.15-fold reduction in AD value was observed when the drug treatment was discontinued for only 5 passages. Second, populations with FMDV including in its 3D (polymerase) replacement M296I -that decreases the sensitivity to R  – displayed a 1.90-fold reduction in AD value (RAp45 and RAp60 populations, Figure 4). The ratio of minimum over maximum mutation frequency yielded higher values (indicative of higher diversity), in populations RAp45 and RAp60 than in population RAp35 which is the most diverse. Therefore mutation frequencies (both maximum, minimum or their ratio) and Shannon entropy, are less definitory than PAQ (when analyzing AD) to characterize the internal genetic diversity of mutagenized viral populations.
A possible interpretation of intrapopulation contractions is that in the course of the mutagenic treatment subsets of genomes with better than average replication efficiency are produced. However, their potential replicative advantage can not be expressed due to the continuous mutational input due to the presence of R. When R pressure is removed, specific subsets of genomes showing higher than average replication capacity, replenish the population. A similar, albeit less pronounced, contraction would occur as a consequence of the dominance of viral mutants with decreased sensitivity to R.
The MARLS-derived populations that replicated in the presence of R (the RA lineage, Figure 1B), showed unexpected high resistance to extinction at passage 35 and subsequent passages, associated with a mutation that decreased the sensitivity to R . Despite such mutation, the number of total mutations fixed in the consensus sequence of population RAp35 was very high (1.13 mutations/passage, Figure 6) when compared with its ancestor (0.2 mutation/passage). Also the molecular clones derived from RAp35 presented a 12-fold increase in mutation frequency with respect to the control population, CAp35 (Table 1). The specific pattern of mutations of FU-treated, AZC-treated and R-treated populations (Figure 5), and the increase in mutant frequency with respect control populations (Table 1) suggest that the increase in mutation frequencies was effectively produced by the action of mutagens, in agreement with the mutational bias induced in the course of polymerization assays in the presence of R-triphosphate or FU-triphosphate using purified FMDV 3D (polymerase) (, Agudo et al submitted for publication).
Despite the high error frequencies induced by mutagen, mutations did not accumulate at random along the viral genome, but rather they accumulated at preferred genomic regions (Figure 7). The non-random accumulation of mutations and the excess of corrected synonymous mutations (Table 1) are consistent with purifying selection operating in the evolution of mutagenized populations. Viruses harbouring less deleterious mutations might have a selective advantage in a landscape of highly damaged viruses. In this view, the quasispecies would display a mutagenesis buffering activity, accepting those mutations that affect the more permissive regions of the viral genome. This dynamics is strikingly parallel to that observed in the course of bottleneck transfers carried out with the same FMDV clones C-S8c1 ([47, 48, 76–78] reviewed in ). In bottleneck transfers individual clones are selected for replication. In this situation Müller ratchet effect operates on the clones . In clones that resisted extinction, mutations accumulated preferentially at some genomic regions and, again, an excess of synonymous mutations was found [48, 77].
The present study supports the "lethal defection" model of viral extinction in that a limited number of mutations per genome can be sufficient to drive a virus to extinction [28, 34, 74]. Again, the comparison with FMDV clones subjected up to 409 bottleneck transfers is highly significant. In such clones extinction was not achieved despite the genomes reaching average mutation frequencies which are 5- to 37-fold higher than those associated with viral extinction by lethal mutagenesis [32, 39, 41, 48]. It has been proposed that in successive bottleneck transfers, the kinetics of mutation accumulation allows the fixation of compensatory beneficial mutations that counteract the Müller Ratchet effect [76–79].
Several theoretical models of lethal mutagenesis of viruses have been proposed, either as a direct consequence of quasispecies dynamics and its corollary concept of error catastrophe, or independently of error catastrophe [70, 81–84]. All models converge in that lethal mutagenesis is a feasible strategy to achieve virus extinction by mutagenic nucleotide analogues. What the models did not take into consideration is the key influence on extinction of internal interactions exerted among components of the mutant spectra. Lethal and interfering mutations impede a substantial "evaporation" (or "diffusion") of genomic sequences in sequence space , as it could not be otherwise, considering that we are dealing with loss of multiple biological functions compactly integrated in a viral genome . The phylogenetic and PAQ approaches used here should be extremely helpful to monitor in a quantitative fashion the evolution of mutagenized viral populations both in cell culture and in vivo, as they increase their mutational load and either succumb or escape extinction.
Phylogenetic and PAQ analyses have unveiled changes in the internal population structure of FMDV viral quasispecies subjected to mutagenesis by base and nucleoside analogues. Expansions and compressions of mutant spectra have been quantitated by comparing average genetic distances among components of mutant spectra. Comparisons of the types and distribution of mutations along the viral genomes have shown that negative (or purifying selection) acts in the course of enhanced mutagenesis. Virus extinction can be achieved with modest increases of population complexity. The average distance parameter (AD) reflects the intrapopulation structure better than the other estimators traditionally used, such as mutation frequencies or Shannon entropy.
The authors wish to thank O.Gordo for critical review of the manuscript, and C Escarmís for important contributions to this project. Work at Centro de Biología Molecular "Severo Ochoa" was supported by grants BFU2006-00863 from MEC, 36558/06 from FIPSE, and Fundación R.Areces. Work at Centro de Astrobiología was supported by INTA, MEC, CAM and UE. CIBERehd is funded by Instituto de Salud Carlos III. S.O. was supported by a predoctoral fellowship from the Ministerio de Educacion y Ciencia.
- Domingo E: Virus Evolution. Fields Virology. Edited by: Knipe DM, Howley PM. 2007, Philadelphia , Lappincott Williams & Wilkins, Vol.12: 389-421. 5thGoogle Scholar
- Domingo E, Biebricher C, Eigen M, Holland JJ: Quasispecies and RNA Virus Evolution: Principles and Consequences. 2001, Austin , Landes BioscienceGoogle Scholar
- Salemi M, Vandamme AM: The phylogeny Handbook. A practical approach to DNA and Protein Phylogeny. 2004, Cambridge , Cambridge University PressGoogle Scholar
- Grenfell BT, Pybus OG, Gog JR, Wood JL, Daly JM, Mumford JA, Holmes EC: Unifying the epidemiological and evolutionary dynamics of pathogens. Science. 2004, 303 (5656): 327-332.View ArticlePubMedGoogle Scholar
- Holmes EC: Comparative Studies of RNA Virus Evolution. Origins and Evolution of Viruses. Edited by: Domingo E, Parrish CR, Holland JJ. 2008, Oxford , Elsevier, 119-134. 2ndView ArticleGoogle Scholar
- Vignuzzi M, Stone JK, Arnold JJ, Cameron CE, Andino R: Quasispecies diversity determines pathogenesis through cooperative interactions in a viral population. Nature. 2006, 439: 344-348.PubMed CentralView ArticlePubMedGoogle Scholar
- Pfeiffer JK, Kirkegaard K: Increased fidelity reduces poliovirus fitness under selective pressure in mice. PLoS Pathog. 2005, 1: 102-110.View ArticleGoogle Scholar
- Rowe CL, Baker SC, Nathan MJ, Fleming JO: Evolution of mouse hepatitis virus: detection and characterization of spike deletion variants during persistent infection. J Virol. 1997, 71 (4): 2959-2969.PubMed CentralPubMedGoogle Scholar
- Pawlotsky JM: Mechanisms of antiviral treatment efficacy and failure in chronic hepatitis C. Antiviral Res. 2003, 59 (1): 1-11.View ArticlePubMedGoogle Scholar
- Farci P, Shimoda A, Coiana A, Diaz G, Peddis G, Melpolder JC, Strazzera A, Chien DY, Munoz SJ, Balestrieri A, Purcell RH, Alter HJ: The outcome of acute hepatitis C predicted by the evolution of the viral quasispecies. Science. 2000, 288 (5464): 339-344.View ArticlePubMedGoogle Scholar
- Farci P, Strazzera R, Alter HJ, Farci S, Degioannis D, Coiana A, Peddis G, Usai F, Serra G, Chessa L, Diaz G, Balestrieri A, Purcell RH: Early changes in hepatitis C viral quasispecies during interferon therapy predict the therapeutic outcome. Proc Natl Acad Sci USA. 2002, 99 (5): 3081-3086.PubMed CentralView ArticlePubMedGoogle Scholar
- Domingo E, Gomez J: Quasispecies and its impact on viral hepatitis. Virus Res. 2007Google Scholar
- Ciota AT, Ngo KA, Lovelace AO, Payne AF, Zhou Y, Shi PY, Kramer LD: Role of the mutant spectrum in adaptation and replication of West Nile virus. J Gen Virol. 2007, 88 (Pt 3): 865-874.PubMed CentralView ArticlePubMedGoogle Scholar
- Jerzak GV, Bernard K, Kramer LD, Shi PY, Ebel GD: The West Nile virus mutant spectrum is host-dependant and a determinant of mortality in mice. Virology. 2007, 360 (2): 469-476.PubMed CentralView ArticlePubMedGoogle Scholar
- Saitou N, Nei M: The neighbor-joining method: a new method for reconstructing phylogenetic trees. Mol Biol Evol. 1987, 4 (4): 406-425.PubMedGoogle Scholar
- Strimmer K, von Haeseler A: Likelihood-mapping: a simple method to visualize phylogenetic content of a sequence alignment. Proc Natl Acad Sci USA. 1997, 94 (13): 6815-6819.PubMed CentralView ArticlePubMedGoogle Scholar
- Fares MA, Barrio E, Becerra N, Escarmis C, Domingo E, Moya A: The foot-and-mouth disease RNA virus as a model in experimental phylogenetics. Int Microbiol. 1998, 1 (4): 311-318.PubMedGoogle Scholar
- Baccam P, Thompson RJ, Fedrigo O, Carpenter S, Cornette JL: PAQ: Partition Analysis of Quasispecies. Bioinformatics. 2001, 17 (1): 16-22.View ArticlePubMedGoogle Scholar
- Baccam P, Thompson RJ, Li Y, Sparks WO, Belshan M, Dorman KS, Wannemuehler Y, Oaks JL, Cornette JL, Carpenter S: Subpopulations of equine infectious anemia virus Rev coexist in vivo and differ in phenotype. J Virol. 2003, 77 (22): 12122-12131.PubMed CentralView ArticlePubMedGoogle Scholar
- Galindo LM, Gaitan-Solis E, Baccam P, Tohme J: Isolation and characterization of RNase LTR sequences of Ty1-copia retrotransposons in common bean (Phaseolus vulgaris L). Genome. 2004, 47 (1): 84-95.View ArticlePubMedGoogle Scholar
- Jones LR, Zandomeni R, Weber EL: Quasispecies in the 5' untranslated genomic region of bovine viral diarrhoea virus from a single individual. J Gen Virol. 2002, 83 (Pt 9): 2161-2168.View ArticlePubMedGoogle Scholar
- Costa-Mattioli M, Domingo E, Cristina J: Analysis of sequential hepatitis A virus strains reveals coexistence of distinct viral subpopulations. J Gen Virol. 2006, 87 (Pt 1): 115-118.View ArticlePubMedGoogle Scholar
- Rowlands DJ: Foot-and-mouth disease. Virus Res. 2003, 91: 1-161.View ArticleGoogle Scholar
- Sobrino F, Domingo E: Foot-and-Mouth Disease: Current Perspectives. Horizon Bioscience. 2004, Wymondham, England Google Scholar
- Picornavirus database. [http://www.picornaviridae.com/]
- Carrillo C, Tulman ER, Delhon G, Lu Z, Carreno A, Vagnozzi A, Kutish GF, Rock DL: Comparative genomics of foot-and-mouth disease virus. J Virol. 2005, 79 (10): 6487-6504.PubMed CentralView ArticlePubMedGoogle Scholar
- Knowles NJ, Samuel AR: Molecular epidemiology of foot-and-mouth disease virus. Virus Res. 2003, 91 (1): 65-80.View ArticlePubMedGoogle Scholar
- Loeb LA, Essigmann JM, Kazazi F, Zhang J, Rose KD, Mullins JI: Lethal mutagenesis of HIV with mutagenic nucleoside analogs. Proc Natl Acad Sci USA. 1999, 96: 1492-1497.PubMed CentralView ArticlePubMedGoogle Scholar
- Anderson JP, Daifuku R, Loeb LA: Viral Error Catastrophe by Mutagenic Nucleosides. Annu Rev Microbiol. 2004, 58: 183-205.View ArticlePubMedGoogle Scholar
- Domingo E: Virus entry into error catastrophe as a new antiviral strategy. Virus Res. 2005, 107: 115-228.View ArticleGoogle Scholar
- Graci JD, Cameron CE: Challenges for the development of ribonucleoside analogues as inducers of error catastrophe. Antivir Chem Chemother. 2004, 15 (1): 1-13.View ArticlePubMedGoogle Scholar
- Gonzalez-Lopez C, Gomez-Mariano G, Escarmis C, Domingo E: Invariant aphthovirus consensus nucleotide sequence in the transition to error catastrophe. Infect Genet Evol. 2005, 5 (4): 366-374.View ArticlePubMedGoogle Scholar
- Gonzalez-Lopez C, Arias A, Pariente N, Gomez-Mariano G, Domingo E: Preextinction viral RNA can interfere with infectivity. J Virol. 2004, 78 (7): 3319-3324.PubMed CentralView ArticlePubMedGoogle Scholar
- Grande-Perez A, Lazaro E, Lowenstein P, Domingo E, Manrubia SC: Suppression of viral infectivity through lethal defection. Proc Natl Acad Sci U S A. 2005, 102 (12): 4448-4452.PubMed CentralView ArticlePubMedGoogle Scholar
- Perales C, Mateo R, Mateu MG, Domingo E: Insights into RNA virus mutant spectrum and lethal mutagenesis events: replicative interference and complementation by multiple point mutants. J Mol Biol. 2007, 369 (4): 985-1000.View ArticlePubMedGoogle Scholar
- Grande-Perez A, Gomez-Mariano G, Lowenstein PR, Domingo E: Mutagenesis-induced, large fitness variations with an invariant arenavirus consensus genomic nucleotide sequence. J Virol. 2005, 79 (16): 10451-10459.PubMed CentralView ArticlePubMedGoogle Scholar
- Roux L, Simon AE, Holland JJ: Effects of defective interfering viruses on virus replication and pathogenesis in vitro and in vivo. Adv Virus Res. 1991, 40: 181-211.View ArticlePubMedGoogle Scholar
- Pariente N, Airaksinen A, Domingo E: Mutagenesis versus inhibition in the efficiency of extinction of foot-and-mouth disease virus. J Virol. 2003, 77 (12): 7131-7138.PubMed CentralView ArticlePubMedGoogle Scholar
- Sierra S, Davila M, Lowenstein PR, Domingo E: Response of foot-and-mouth disease virus to increased mutagenesis: influence of viral load and fitness in loss of infectivity. J Virol. 2000, 74 (18): 8316-8323.PubMed CentralView ArticlePubMedGoogle Scholar
- Sierra M, Airaksinen A, Gonzalez-Lopez C, Agudo R, Arias A, Domingo E: Foot-and-mouth disease virus mutant with decreased sensitivity to ribavirin: implications for error catastrophe. J Virol. 2007, 81 (4): 2012-2024.PubMed CentralView ArticlePubMedGoogle Scholar
- Pariente N, Sierra S, Lowenstein PR, Domingo E: Efficient virus extinction by combinations of a mutagen and antiviral inhibitors. J Virol. 2001, 75 (20): 9723-9730.PubMed CentralView ArticlePubMedGoogle Scholar
- Pariente N, Sierra S, Airaksinen A: Action of mutagenic agents and antiviral inhibitors on foot-and-mouth disease virus. Virus Res. 2005, 107 (2): 183-193.View ArticlePubMedGoogle Scholar
- Domingo E, Davila M, Ortin J: Nucleotide sequence heterogeneity of the RNA from a natural population of foot-and-mouth-disease virus. Gene. 1980, 11 (3-4): 333-346.View ArticlePubMedGoogle Scholar
- Sobrino F, Davila M, Ortin J, Domingo E: Multiple genetic variants arise in the course of replication of foot-and-mouth disease virus in cell culture. Virology. 1983, 128 (2): 310-318.View ArticlePubMedGoogle Scholar
- Charpentier N, Davila M, Domingo E, Escarmis C: Long-term, large-population passage of aphthovirus can generate and amplify defective noninterfering particles deleted in the leader protease gene. Virology. 1996, 223 (1): 10-18.View ArticlePubMedGoogle Scholar
- Garcia-Arriaza J, Manrubia SC, Toja M, Domingo E, Escarmis C: Evolutionary transition toward defective RNAs that are infectious by complementation. J Virol. 2004, 78 (21): 11678-11685.PubMed CentralView ArticlePubMedGoogle Scholar
- Escarmis C, Davila M, Charpentier N, Bracho A, Moya A, Domingo E: Genetic lesions associated with Muller's ratchet in an RNA virus. J Mol Biol. 1996, 264 (2): 255-267.View ArticlePubMedGoogle Scholar
- Escarmís C, Lázaro E, A. A, Domingo E: Repeated bottleneck transfers can lead to non-cytocidal forms of a cytopahic virus. Implications for viral extinction. J Mol Biol. 2008, 376 (2): 367-379.View ArticlePubMedGoogle Scholar
- Ferrer-Orta C, Arias A, Perez-Luque R, Escarmis C, Domingo E, Verdaguer N: Structure of foot-and-mouth disease virus RNA-dependent RNA polymerase and its complex with a template-primer RNA. J Biol Chem. 2004, 279 (45): 47212-47221.View ArticlePubMedGoogle Scholar
- Arias A, Agudo R, Ferrer-Orta C, Perez-Luque R, Airaksinen A, Brocchi E, Domingo E, Verdaguer N, Escarmis C: Mutant viral polymerase in the transition of virus to error catastrophe identifies a critical site for RNA binding. J Mol Biol. 2005, 353 (5): 1021-1032.View ArticlePubMedGoogle Scholar
- Toja M, Escarmis C, Domingo E: Genomic nucleotide sequence of a foot-and-mouth disease virus clone and its persistent derivatives. Implications for the evolution of viral quasispecies during a persistent infection. Virus Res. 1999, 64 (2): 161-171.View ArticlePubMedGoogle Scholar
- Sambrook J, Russell DW: Molecular Cloning: A Laboratory Manual. 2001, Cold Spring Harbor, New York , Cold Spring Harbor Laboratory Press, 3rdGoogle Scholar
- Airaksinen A, Pariente N, Menendez-Arias L, Domingo E: Curing of foot-and-mouth disease virus from persistently infected cells by ribavirin involves enhanced mutagenesis. Virology. 2003, 311 (2): 339-349.View ArticlePubMedGoogle Scholar
- Thompson JD, Higgins DG, Gibson TJ: CLUSTAL W: improving the sensitivity of progressive multiple sequence alignment through sequence weighting, position-specific gap penalties and weight matrix choice. Nucleic Acids Res. 1994, 22 (22): 4673-4680.PubMed CentralView ArticlePubMedGoogle Scholar
- Hall TA: BioEdit: a user-friendly biological sequence alignment editor and analysis program for Windows 95/98/NT. Nucl Acids SympSer. 1999, 41: 95-98.Google Scholar
- Kumar S, Tamura K, Nei M: MEGA3: Integrated software for Molecular Evolutionary Genetics Analysis and sequence alignment. Brief Bioinform. 2004, 5 (2): 150-163.View ArticlePubMedGoogle Scholar
- Kimura M: A simple method for estimating evolutionary rates of base substitutions through comparative studies of nucleotide sequences. J Mol Evol. 1980, 16 (2): 111-120.View ArticlePubMedGoogle Scholar
- Felsenstein J: Confidence limits on phylogenies: an approach using the bootstrap. Evolution. 1985, 39: 738-791.View ArticleGoogle Scholar
- Sober E: Reconstructing the past: Parsimony, Evolution and Inference. 1988, Cambridge, MA , MIT PressGoogle Scholar
- Tamura K, Nei M: Estimation of the number of nucleotide substitutions in the control region of mitochondrial DNA in humans and chimpanzees. Mol Biol Evol. 1993, 10 (3): 512-526.PubMedGoogle Scholar
- Keane TM, Creevey CJ, Pentony MM, Naughton TJ, McLnerney JO: Assessment of methods for amino acid matrix selection and their use on empirical data shows that ad hoc assumptions for choice of matrix are not justified. BMC Evol Biol. 2006, 6: 29-PubMed CentralView ArticlePubMedGoogle Scholar
- Guindon S, Lethiec F, Duroux P, Gascuel O: PHYML Online--a web server for fast maximum likelihood-based phylogenetic inference. Nucleic Acids Res. 2005, 33 (Web Server issue): W557-9.PubMed CentralView ArticlePubMedGoogle Scholar
- Phylogeny.fr. [http://www.phylogeny.fr/phylo_cgi/phyml]
- Anisimova M, Gascuel O: Approximate likelihood-ratio test for branches: A fast, accurate, and powerful alternative. Syst Biol. 2006, 55 (4): 539-552.View ArticlePubMedGoogle Scholar
- Volkenstein MV: Physical approaches to biological evolution. 1994, Berlin , Springer-VerlagView ArticleGoogle Scholar
- HCV sequence database. [http://hcv.lanl.gov/content/sequence/SNAP/perlsnap.html]
- Korber B: HIV sequence signatures and similarities. Computational and evolutionary analysis of HIV molecular sequences. Edited by: Rodrigo A, Learn GH. 2000, Dordrech: KluberGoogle Scholar
- Comeron JM: K-Estimator: calculation of the number of nucleotide substitutions per site and the confidence intervals. Bioinformatics. 1999, 15 (9): 763-764.View ArticlePubMedGoogle Scholar
- Felsenstein J: PHYLIP - Phylogeny inference package. Cladistics. 1989, 5: 164-166.Google Scholar
- Eigen M: Error catastrophe and antiviral strategy. Proc Natl Acad Sci USA. 2002, 99 (21): 13374-13376.PubMed CentralView ArticlePubMedGoogle Scholar
- Graci JD, Cameron CE: Quasispecies, error catastrophe, and the antiviral activity of ribavirin. Virology. 2002, 298 (2): 175-180.View ArticlePubMedGoogle Scholar
- Pincus SE, Wimmer E: Production of guanidine-resistant and -dependent poliovirus mutants from cloned cDNA: mutations in polypeptide 2C are directly responsible for altered guanidine sensitivity. J Virol. 1986, 60: 793-796.PubMed CentralPubMedGoogle Scholar
- Holland JJ, Domingo E, de la Torre JC, Steinhauer DA: Mutation frequencies at defined single codon sites in vesicular stomatitis virus and poliovirus can be increased only slightly by chemical mutagenesis. J Virol. 1990, 64: 3960-3962.PubMed CentralPubMedGoogle Scholar
- Tapia N, Fernandez G, Parera M, Gomez-Mariano G, Clotet B, Quinones-Mateu M, Domingo E, Martinez MA: Combination of a mutagenic agent with a reverse transcriptase inhibitor results in systematic inhibition of HIV-1 infection. Virology. 2005, 338 (1): 1-8.View ArticlePubMedGoogle Scholar
- Grande-Perez A, Sierra S, Castro MG, Domingo E, Lowenstein PR: Molecular indetermination in the transition to error catastrophe: systematic elimination of lymphocytic choriomeningitis virus through mutagenesis does not correlate linearly with large increases in mutant spectrum complexity. Proc Natl Acad Sci U S A. 2002, 99 (20): 12938-12943.PubMed CentralView ArticlePubMedGoogle Scholar
- Manrubia SC, Escarmis C, Domingo E, Lazaro E: High mutation rates, bottlenecks, and robustness of RNA viral quasispecies. Gene. 2005, 347 (2): 273-282.View ArticlePubMedGoogle Scholar
- Escarmis C, Gomez-Mariano G, Davila M, Lazaro E, Domingo E: Resistance to extinction of low fitness virus subjected to plaque-to-plaque transfers: diversification by mutation clustering. J Mol Biol. 2002, 315 (4): 647-661.View ArticlePubMedGoogle Scholar
- Escarmis C, Davila M, Domingo E: Multiple molecular pathways for fitness recovery of an RNA virus debilitated by operation of Muller's ratchet. J Mol Biol. 1999, 285 (2): 495-505.View ArticlePubMedGoogle Scholar
- Escarmís C, Lázaro E, Manrubia SC: Population bottlenecks in quasispecies dynamics. Curr Top Microbiol Immunol. 2006, 299: 141-170.PubMedGoogle Scholar
- Muller HJ: The relation of recombination to mutational advance. Mutat Res. 1964, 1: 2-9.View ArticleGoogle Scholar
- Biebricher CK, Eigen M: The error threshold. Virus Res. 2005, 107 (2): 117-127.View ArticlePubMedGoogle Scholar
- Bull JJ, Sanjuan R, Wilke CO: Theory of lethal mutagenesis for viruses. J Virol. 2007, 81 (6): 2930-2939.PubMed CentralView ArticlePubMedGoogle Scholar
- Bull JJ, Wilke CO: Lethal Mutagenesis. Origin and evolution of viruses. Edited by: Domingo E, Parrish C, Holland JJ. 2008, Amsterdan, San Diego , Elsevier (In Press), 2ndGoogle Scholar
- Schuster P, Stadler PF: Early Replication: Origin and Evolution. Origin and Evolution of Viruses. Edited by: Domingo E, Parrish C, Holland JJ. 2008, Amsterdan, San Diego , Elsevier (In press), 2ndGoogle Scholar
- Nei M, Gojobori T: Simple methods for estimating the numbers of synonymous and nonsynonymous nucleotide substitutions. Mol Biol Evol. 1986, 3 (5): 418-426.PubMedGoogle Scholar
- Korber B: HIV signature and sequence variation analysis. Comptutational Analysis of HIV molecular sequences. Edited by: Rodrigo AG, Learn GH. 2001, Dordrecht, Netherlands , Kluwer Academic Publishers, 4: 55-72.Google Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.