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

Quasispecies theory and emerging viruses: challenges and applications

Quasispecies theory revolutionized our understanding of viral evolution by describing viruses as dynamic populations of genetically diverse variants constantly adapting. This article explores the theory's role in virus-host interactions, including immune evasion, drug resistance, and viral emergence. We review the original model, recent advances, and key virus dynamics needing incorporation into quasispecies theory. We introduce the ultracube concept as a more realistic multidimensional sequence space to investigate virus evolutionary dynamics.

Current sampling and sequencing biases of Lassa mammarenavirus limit inference from phylogeography and molecular epidemiology in Lassa fever endemic regions

Lassa fever (LF) is a potentially lethal viral haemorrhagic infection of humans caused by Lassa mammarenavirus (LASV). It is an important endemic zoonotic disease in West Africa with growing evidence for increasing frequency and sizes of outbreaks. Phylogeographic and molecular epidemiology methods have projected expansion of the Lassa fever endemic zone in the context of future global change. The Natal multimammate mouse (Mastomys natalensis) is the predominant LASV reservoir, with few studies investigating the role of other animal species. To explore host sequencing biases, all LASV nucleotide sequences and associated metadata available on GenBank (n = 2,298) were retrieved. Most data originated from Nigeria (54%), Guinea (20%) and Sierra Leone (14%). Data from non-human hosts (n = 703) were limited and only 69 sequences encompassed complete genes. We found a strong positive correlation between the number of confirmed human cases and sequences at the country level (r = 0.93 (95% Confidence Interval = 0.71-0.98), p < 0.001) but no correlation exists between confirmed cases and the number of available rodent sequences (r = -0.019 (95% C.I. -0.71-0.69), p = 0.96). Spatial modelling of sequencing effort highlighted current biases in locations of available sequences, with increased sequencing effort observed in Southern Guinea and Southern Nigeria. Phylogenetic analyses showed geographic clustering of LASV lineages, suggestive of isolated events of human-to-rodent transmission and the emergence of currently circulating strains of LASV from the year 1498 in Nigeria. Overall, the current study highlights significant geographic limitations in LASV surveillance, particularly, in non-human hosts. Further investigation of the non-human reservoir of LASV, alongside expanded surveillance, are required for precise characterisation of the emergence and dispersal of LASV. Accurate surveillance of LASV circulation in non-human hosts is vital to guide early detection and initiation of public health interventions for future Lassa fever outbreaks.

Mammarenavirus Genetic Diversity and Its Biological Implications

Members of the family Arenaviridae are classified into four genera: Antennavirus, Hartmanivirus, Mammarenavirus, and Reptarenavirus. Reptarenaviruses and hartmaniviruses infect (captive) snakes and have been shown to cause boid inclusion body disease (BIBD). Antennaviruses have genomes consisting of 3, rather than 2, segments, and were discovered in actinopterygian fish by next-generation sequencing but no biological isolate has been reported yet. The hosts of mammarenaviruses are mainly rodents and infections are generally asymptomatic. Current knowledge about the biology of reptarenaviruses, hartmaniviruses, and antennaviruses is very limited and their zoonotic potential is unknown. In contrast, some mammarenaviruses are associated with zoonotic events that pose a threat to human health. This review will focus on mammarenavirus genetic diversity and its biological implications. Some mammarenaviruses including lymphocytic choriomeningitis virus (LCMV) are excellent experimental model systems for the investigation of acute and persistent viral infections, whereas others including Lassa (LASV) and Junin (JUNV) viruses, the causative agents of Lassa fever (LF) and Argentine hemorrhagic fever (AHF), respectively, are important human pathogens. Mammarenaviruses were thought to have high degree of intra-and inter-species amino acid sequence identities, but recent evidence has revealed a high degree of mammarenavirus genetic diversity in the field. Moreover, closely related mammarenavirus can display dramatic phenotypic differences in vivo. These findings support a role of genetic variability in mammarenavirus adaptability and pathogenesis. Here, we will review the molecular biology of mammarenaviruses, phylogeny, and evolution, as well as the quasispecies dynamics of mammarenavirus populations and their biological implications.

Genetic diversity through social heterosis can increase virulence in RNA viral infections and cancer progression

In viral infections and cancer tumours, negative health outcomes often correlate with increasing genetic diversity. Possible evolutionary processes for such relationships include mutant lineages escaping host control or diversity, per se, creating too many immune system targets. Another possibility is social heterosis where mutations and replicative errors create clonal lineages varying in intrinsic capability for successful dispersal; improved environmental buffering; resource extraction or effective defence against immune systems. Rather than these capabilities existing in one genome, social heterosis proposes complementary synergies occur across lineages in close proximity. Diverse groups overcome host defences as interacting 'social genomes' with group genetic tool kits exceeding limited individual plasticity. To assess the possibility of social heterosis in viral infections and cancer progression, we conducted extensive literature searches for examples consistent with general and specific predictions from the social heterosis hypothesis. Numerous studies found supportive patterns in cancers across multiple tissues and in several families of RNA viruses. In viruses, social heterosis mechanisms probably result from long coevolutionary histories of competition between pathogen and host. Conversely, in cancers, social heterosis is a by-product of recent mutations. Investigating how social genomes arise and function in viral quasi-species swarms and cancer tumours may lead to new therapeutic approaches.

Quasispecies characteristic in "a" determinant region is a potential predictor for the risk of immunoprophylaxis failure of mother-to-child-transmission of sub-genotype C2 hepatitis B virus: a prospective nested case-control study

OBJECTIVE

This study was performed to explore the correlation between the characteristics of hepatitis B virus (HBV) quasispecies in HBV-infected pregnant women and the risk of immunoprophylaxis failure for their infants.

DESIGN

In this prospective nested case-control study, the characteristics of HBV quasispecies in mothers whose infants were immunoprophylaxis success (control group) and those whose infants were immunoprophylaxis failure (case group) were analysed by the clone-based sequencing of full-length HBV genome and next-generation sequencing (NGS) of "a" determinant region, and were compared between the two groups.

RESULTS

The quasispecies characteristics including mutant frequency, Shannon entropy and mean genetic distance at amino acid level of "a" determinant region were significantly lower in case group than that in control group, using the full-length HBV genome clone-based sequencing assay. These results were confirmed by NGS assay. Notably, we discovered that the differences were also significant at nucleotide level by NGS assay. Furthermore, the risk of immunoprophylaxis failure could be predicted by analysing the three HBV quasispecies characteristics either at nucleotide level or at amino acid level of "a" determinant region, and the corresponding predictive values were tentatively set up.

CONCLUSIONS

HBV quasispecies with a more complex mutant spectrum in "a" determinant region might be more vulnerable to extinct through mother-to-child-transmission (MTCT). More importantly, analysing HBV quasispecies characteristics in pregnant women with high HBV DNA load might be helpful to predict the high-risk population of immunoprophylaxis failure, and consequently provide accurate intervention against MTCT of HBV.

The increasing impact of lethal mutagenesis of viruses

Selection of viral mutants resistant to compounds used in therapy is a major determinant of treatment failure, a problem akin to antibiotic resistance in bacteria. In this scenario, mutagenic base and nucleoside analogs have entered the picture because they increase the mutation rate of viral populations to levels incompatible with their survival. This antiviral strategy is termed lethal mutagenesis. It has found a major impulse with the observation that some antiviral agents, which initially were considered only inhibitors of virus multiplication, may in effect exert part of their antiviral activity through mutagenesis. Here, we review the conceptual basis of lethal mutagenesis, the evidence of virus extinction through mutagenic nucleotide analogs and prospects for application in antiviral designs.

Contribution of a Multifunctional Polymerase Region of Foot-and-Mouth Disease Virus to Lethal Mutagenesis

Viral RNA-dependent RNA polymerases (RdRps) are major determinants of high mutation rates and generation of mutant spectra that mediate RNA virus adaptability. The RdRp of the picornavirus foot-and-mouth disease virus (FMDV), termed 3D, is a multifunctional protein that includes a nuclear localization signal (NLS) in its N-terminal region. Previous studies documented that some amino acid substitutions within the NLS altered nucleotide recognition and enhanced the incorporation of the mutagenic purine analogue ribavirin in viral RNA, but the mutants tested were not viable and their response to lethal mutagenesis could not be studied. Here we demonstrate that NLS amino acid substitution M16A of FMDV serotype C does not affect infectious virus production but accelerates ribavirin-mediated virus extinction. The mutant 3D displays polymerase activity, RNA binding, and copying processivity that are similar to those of the wild-type enzyme but shows increased ribavirin-triphosphate incorporation. Crystal structures of the mutant 3D in the apo and RNA-bound forms reveal an expansion of the template entry channel due to the replacement of the bulky Met by Ala. This is a major difference with other 3D mutants with altered nucleotide analogue recognition. Remarkably, two distinct loop β9-α11 conformations distinguish 3Ds that exhibit higher or lower ribavirin incorporation than the wild-type enzyme. This difference identifies a specific molecular determinant of ribavirin sensitivity of FMDV. Comparison of several polymerase mutants indicates that different domains of the molecule can modify nucleotide recognition and response to lethal mutagenesis. The connection of this observation with current views on quasispecies adaptability is discussed.IMPORTANCE The nuclear localization signal (NLS) of the foot-and-mouth disease virus (FMDV) polymerase includes residues that modulate the sensitivity to mutagenic agents. Here we have described a viable NLS mutant with an amino acid replacement that facilitates virus extinction by ribavirin. The corresponding polymerase shows increased incorporation of ribavirin triphosphate and local structural modifications that implicate the template entry channel. Specifically, comparison of the structures of ribavirin-sensitive and ribavirin-resistant FMDV polymerases has identified loop β9-α11 conformation as a determinant of sensitivity to ribavirin mutagenesis.

Viral replication modes in single-peak fitness landscapes: A dynamical systems analysis

Positive-sense, single-stranded RNA viruses are important pathogens infecting almost all types of organisms. Experimental evidence from distributions of mutations and from viral RNA amplification suggest that these pathogens may follow different RNA replication modes, ranging from the stamping machine replication (SMR) to the geometric replication (GR) mode. Although previous theoretical work has focused on the evolutionary dynamics of RNA viruses amplifying their genomes with different strategies, little is known in terms of the bifurcations and transitions involving the so-called error threshold (mutation-induced dominance of mutants) and lethal mutagenesis (extinction of all sequences due to mutation accumulation and demographic stochasticity). Here we analyze a dynamical system describing the intracellular amplification of viral RNA genomes evolving on a single-peak fitness landscape focusing on three cases considering neutral, deleterious, and lethal mutants. We analytically derive the critical mutation rates causing lethal mutagenesis and error threshold, governed by transcritical bifurcations that depend on parameters α (parameter introducing the mode of replication), replicative fitness of mutants (k₁), and on the spontaneous degradation rates of the sequences (ϵ). Our results relate the error catastrophe with lethal mutagenesis in a model with continuous populations of viral genomes. The former case involves dominance of the mutant sequences, while the latter, a deterministic extinction of the viral RNAs during replication due to increased mutation. For the lethal case the critical mutation rate involving lethal mutagenesis is μ_c=1-ɛ/α. Here, the SMR involves lower critical mutation rates, being the system more robust to lethal mutagenesis replicating closer to the GR mode. This result is also found for the neutral and deleterious cases, but for these later cases lethal mutagenesis can shift to the error threshold once the replication mode surpasses a threshold given by α=ϵ/k₁.

Extinction of Zika Virus and Usutu Virus by Lethal Mutagenesis Reveals Different Patterns of Sensitivity to Three Mutagenic Drugs

Flaviviruses constitute an increasing source of public health concern, with growing numbers of pathogens causing disease and geographic spread to temperate climates. Despite a large body of evidence supporting mutagenesis as a conceivable antiviral strategy, there are currently no data on the sensitivity to increased mutagenesis for Zika virus (ZIKV) and Usutu virus (USUV), two emerging flaviviral threats.

Flaviviruses constitute an increasing source of public health concern, with growing numbers of pathogens causing disease and geographic spread to temperate climates. Despite a large body of evidence supporting mutagenesis as a conceivable antiviral strategy, there are currently no data on the sensitivity to increased mutagenesis for Zika virus (ZIKV) and Usutu virus (USUV), two emerging flaviviral threats. In this study, we demonstrate that both viruses are sensitive to three ribonucleosides, favipiravir, ribavirin, and 5-fluorouracil, that have shown mutagenic activity against other RNA viruses while remaining unaffected by a mutagenic deoxyribonucleoside. Serial cell culture passages of ZIKV in the presence of these compounds resulted in the rapid extinction of infectivity, suggesting elevated sensitivity to mutagenesis. USUV extinction was achieved when a 10-fold dilution was applied between every passage, but not in experiments involving undiluted virus, indicating an overall lower susceptibility than ZIKV. Although the two viruses are inhibited by the same three drugs, ZIKV is relatively more susceptive to serial passage in the presence of purine analogues (favipiravir and ribavirin), while USUV replication is suppressed more efficiently by 5-fluorouracil. These differences in sensitivity typically correlate with the increases in the mutation frequencies observed in each nucleoside treatment. These results are relevant to the development of efficient therapies based on lethal mutagenesis and support the rational selection of different mutagenic nucleosides for each pathogen. We will discuss the implications of these results to the fidelity of flavivirus replication and the design of antiviral therapies based on lethal mutagenesis.

Resistance of high fitness hepatitis C virus to lethal mutagenesis

Viral fitness quantifies the degree of virus adaptation to a given environment. How viral fitness can influence the mutant spectrum complexity of a viral quasispecies subjected to lethal mutagenesis has not been investigated. Here we document that two high fitness hepatitis C virus populations display higher resistance to the mutagenic nucleoside analogues favipiravir and ribavirin than their parental, low fitness HCV. All populations, however, exhibited a mutation transition bias indicative of active mutagenesis. Resistance to the analogues was associated with a limited expansion of mutant spectrum complexity, as evidenced by several diversity indices used to characterize mutant spectra. The results are consistent with a replicative site-drug competition mechanism that was previously proposed for HCV fitness-associated resistance to non-mutagenic inhibitors. Other alternative, non-mutually exclusive mechanisms are considered. The results introduce viral fitness as a relevant parameter to evaluate the response of viruses to lethal mutagenesis, with implications for antiviral designs.

The fundamental theorem of natural selection with mutations

The mutation-selection process is the most fundamental mechanism of evolution. In 1935, R. A. Fisher proved his fundamental theorem of natural selection, providing a model in which the rate of change of mean fitness is equal to the genetic variance of a species. Fisher did not include mutations in his model, but believed that mutations would provide a continual supply of variance resulting in perpetual increase in mean fitness, thus providing a foundation for neo-Darwinian theory. In this paper we re-examine Fisher's Theorem, showing that because it disregards mutations, and because it is invalid beyond one instant in time, it has limited biological relevance. We build a differential equations model from Fisher's first principles with mutations added, and prove a revised theorem showing the rate of change in mean fitness is equal to genetic variance plus a mutational effects term. We refer to our revised theorem as the fundamental theorem of natural selection with mutations. Our expanded theorem, and our associated analyses (analytic computation, numerical simulation, and visualization), provide a clearer understanding of the mutation-selection process, and allow application of biologically realistic parameters such as mutational effects. The expanded theorem has biological implications significantly different from what Fisher had envisioned.

The Janus Face of Follicular T Helper Cells in Chronic Viral Infections

Chronic infections with non-cytopathic viruses constitutively expose virus-specific adaptive immune cells to cognate antigen, requiring their numeric and functional adaptation. Virus-specific CD8 T cells are compromised by various means in their effector functions, collectively termed T cell exhaustion. Alike CD8 T cells, virus-specific CD4 Th1 cell responses are gradually downregulated but instead, follicular T helper (T_FH) cell differentiation and maintenance is strongly promoted during chronic infection. Thereby, the immune system promotes antibody responses, which bear less immune-pathological risk compared to cytotoxic and pro-inflammatory T cell responses. This emphasis on T_FH cells contributes to tolerance of the chronic infection and is pivotal for the continued maturation and adaptation of the antibody response, leading eventually to the emergence of virus-neutralizing antibodies, which possess the potential to control the established chronic infection. However, sustained high levels of T_FH cells can also result in a less stringent B cell selection process in active germinal center reactions, leading to the activation of virus-unspecific B cells, including self-reactive B cells, and to hypergammaglobulinemia. This dispersal of B cell help comes at the expense of a stringently selected virus-specific antibody response, thereby contributing to its delayed maturation. Here, we discuss these opposing facets of T_FH cells in chronic viral infections.

Lethal mutagenesis of an RNA plant virus via lethal defection

Lethal mutagenesis is an antiviral therapy that relies on increasing the viral mutation rate with mutagenic nucleoside or base analogues. Currently, the molecular mechanisms that lead to virus extinction through enhanced mutagenesis are not fully understood. Increasing experimental evidence supports the lethal defection model of lethal mutagenesis of RNA viruses, where replication-competent-defectors drive infective virus towards extinction. Here, we address lethal mutagenesis in vivo using 5-fluorouracil (5-FU) during the establishment of tobacco mosaic virus (TMV) systemic infections in N. tabacum. The results show that 5-FU decreased the infectivity of TMV without affecting its viral load. Analysis of molecular clones spanning two genomic regions showed an increase of the FU-related base transitions A → G and U → C. Although the mutation frequency or the number of mutations per molecule did not increase, the complexity of the mutant spectra and the distribution of the mutations were altered. Overall, our results suggest that 5-FU antiviral effect on TMV is associated with the perturbation of the mutation-selection balance in the genomic region of the RNA-dependent RNA polymerase (RdRp). Our work supports the lethal defection model for lethal mutagenesis in vivo in a plant RNA virus and opens the way to study lethal mutagens in plant-virus systems.

The Combined Effect of Oseltamivir and Favipiravir on Influenza A Virus Evolution

Influenza virus inflicts a heavy death toll annually and resistance to existing antiviral drugs has generated interest in the development of agents with novel mechanisms of action. Favipiravir is an antiviral drug that acts by increasing the genome-wide mutation rate of influenza A virus (IAV). Potential synergistic benefits of combining oseltamivir and favipiravir have been demonstrated in animal models of influenza, but the population-level effects of combining the drugs are unknown. In order to elucidate the underlying evolutionary processes at play, we performed genome-wide sequencing of IAV experimental populations subjected to serial passaging in vitro under a combined protocol of oseltamivir and favipiravir. We describe the interplay between mutation, selection, and genetic drift that ultimately culminates in population extinction. In particular, selective sweeps around oseltamivir resistance mutations reduce genome-wide variation while deleterious mutations hitchhike to fixation given the increased mutational load generated by favipiravir. This latter effect reduces viral fitness and accelerates extinction compared with IAV populations treated with favipiravir alone, but risks spreading both established and newly emerging mutations, including possible drug resistance mutations, if transmission occurs before the viral populations are eradicated.

Molecular and Functional Bases of Selection against a Mutation Bias in an RNA Virus

The selective pressures acting on viruses that replicate under enhanced mutation rates are largely unknown. Here, we describe resistance of foot-and-mouth disease virus to the mutagen 5-fluorouracil (FU) through a single polymerase substitution that prevents an excess of A to G and U to C transitions evoked by FU on the wild-type foot-and-mouth disease virus, while maintaining the same level of mutant spectrum complexity. The polymerase substitution inflicts upon the virus a fitness loss during replication in absence of FU but confers a fitness gain in presence of FU. The compensation of mutational bias was documented by in vitro nucleotide incorporation assays, and it was associated with structural modifications at the N-terminal region and motif B of the viral polymerase. Predictions of the effect of mutations that increase the frequency of G and C in the viral genome and encoded polymerase suggest multiple points in the virus life cycle where the mutational bias in favor of G and C may be detrimental. Application of predictive algorithms suggests adverse effects of the FU-directed mutational bias on protein stability. The results reinforce modulation of nucleotide incorporation as a lethal mutagenesis-escape mechanism (that permits eluding virus extinction despite replication in the presence of a mutagenic agent) and suggest that mutational bias can be a target of selection during virus replication.

Epistatic Interactions within the Influenza A Virus Polymerase Complex Mediate Mutagen Resistance and Replication Fidelity

RNA viruses exist as genetically diverse populations. This standing genetic diversity gives them the potential to adapt rapidly, evolve resistance to antiviral therapeutics, and evade immune responses. Viral mutants with altered mutation rates or mutational tolerance have provided insights into how genetic diversity arises and how it affects the behavior of RNA viruses. To this end, we identified variants within the polymerase complex of influenza virus that are able tolerate drug-mediated increases in viral mutation rates. We find that drug resistance is highly dependent on interactions among mutations in the polymerase complex. In contrast to other viruses, influenza virus counters the effect of higher mutation rates primarily by maintaining high levels of genome replication. These findings suggest the importance of maintaining large population sizes for viruses with high mutation rates and show that multiple proteins can affect both mutation rate and genome synthesis.

Lethal mutagenesis is a broad-spectrum antiviral strategy that employs mutagenic nucleoside analogs to exploit the high mutation rate and low mutational tolerance of many RNA viruses. Studies of mutagen-resistant viruses have identified determinants of replicative fidelity and the importance of mutation rate to viral population dynamics. We have previously demonstrated the effective lethal mutagenesis of influenza A virus using three nucleoside analogs as well as the virus’s high genetic barrier to mutagen resistance. Here, we investigate the mutagen-resistant phenotypes of mutations that were enriched in drug-treated populations. We find that PB1 T123A has higher replicative fitness than the wild type, PR8, and maintains its level of genome production during 5-fluorouracil (2,4-dihydroxy-5-fluoropyrimidine) treatment. Surprisingly, this mutagen-resistant variant also has an increased baseline rate of C-to-U and G-to-A mutations. A second drug-selected mutation, PA T97I, interacts epistatically with PB1 T123A to mediate high-level mutagen resistance, predominantly by limiting the inhibitory effect of nucleosides on polymerase activity. Consistent with the importance of epistatic interactions in the influenza virus polymerase, our data suggest that nucleoside analog resistance and replication fidelity are strain dependent. Two previously identified ribavirin {1-[(2R,3R,4S,5R)-3,4-dihydroxy-5-(hydroxymethyl)oxolan-2-yl]-1H-1,2,4-triazole-3-carboxamide} resistance mutations, PB1 V43I and PB1 D27N, do not confer drug resistance in the PR8 background, and the PR8-PB1 V43I polymerase exhibits a normal baseline mutation rate. Our results highlight the genetic complexity of the influenza A virus polymerase and demonstrate that increased replicative capacity is a mechanism by which an RNA virus can counter the negative effects of elevated mutation rates.

IMPORTANCE RNA viruses exist as genetically diverse populations. This standing genetic diversity gives them the potential to adapt rapidly, evolve resistance to antiviral therapeutics, and evade immune responses. Viral mutants with altered mutation rates or mutational tolerance have provided insights into how genetic diversity arises and how it affects the behavior of RNA viruses. To this end, we identified variants within the polymerase complex of influenza virus that are able to tolerate drug-mediated increases in viral mutation rates. We find that drug resistance is highly dependent on interactions among mutations in the polymerase complex. In contrast to other viruses, influenza virus counters the effect of higher mutation rates primarily by maintaining high levels of genome replication. These findings suggest the importance of maintaining large population sizes for viruses with high mutation rates and show that multiple proteins can affect both mutation rate and genome synthesis.

Impact of increased mutagenesis on adaptation to high temperature in bacteriophage Qβ

RNA viruses replicate with very high error rates, which makes them more sensitive to additional increases in this parameter. This fact has inspired an antiviral strategy named lethal mutagenesis, which is based on the artificial increase of the error rate above a threshold incompatible with virus infectivity. A relevant issue concerning lethal mutagenesis is whether incomplete treatments might enhance the adaptive possibilities of viruses. We have addressed this question by subjecting an RNA virus, the bacteriophage Qβ, to different transmission regimes in the presence or the absence of sublethal concentrations of the mutagenic nucleoside analogue 5-azacytidine (AZC). Populations obtained were subsequently exposed to a non-optimal temperature and analyzed to determine their consensus sequences. Our results show that previously mutagenized populations rapidly fixed a specific set of mutations upon propagation at the new temperature, suggesting that the expansion of the mutant spectrum caused by AZC has an influence on later evolutionary behavior.

Viral quasispecies complexity measures

Mutant spectrum dynamics (changes in the related mutants that compose viral populations) has a decisive impact on virus behavior. The several platforms of next generation sequencing (NGS) to study viral quasispecies offer a magnifying glass to study viral quasispecies complexity. Several parameters are available to quantify the complexity of mutant spectra, but they have limitations. Here we critically evaluate the information provided by several population diversity indices, and we propose the introduction of some new ones used in ecology. In particular we make a distinction between incidence, abundance and function measures of viral quasispecies composition. We suggest a multidimensional approach (complementary information contributed by adequately chosen indices), propose some guidelines, and illustrate the use of indices with a simple example. We apply the indices to three clinical samples of hepatitis C virus that display different population heterogeneity. Areas of virus biology in which population complexity plays a role are discussed.

Arenavirus Quasispecies and Their Biological Implications

The family Arenaviridae currently comprises over 20 viral species, each of them associated with a main rodent species as the natural reservoir and in one case possibly phyllostomid bats. Moreover, recent findings have documented a divergent group of arenaviruses in captive alethinophidian snakes. Human infections occur through mucosal exposure to aerosols or by direct contact of abraded skin with infectious materials. Arenaviruses merit interest both as highly tractable experimental model systems to study acute and persistent infections and as clinically important human pathogens including Lassa (LASV) and Junin (JUNV) viruses, the causative agents of Lassa and Argentine hemorrhagic fevers (AHFs), respectively, for which there are no FDA-licensed vaccines, and current therapy is limited to an off-label use of ribavirin (Rib) that has significant limitations. Arenaviruses are enveloped viruses with a bi-segmented negative strand (NS) RNA genome. Each genome segment, L (ca 7.3 kb) and S (ca 3.5 kb), uses an ambisense coding strategy to direct the synthesis of two polypeptides in opposite orientation, separated by a noncoding intergenic region (IGR). The S genomic RNA encodes the virus nucleoprotein (NP) and the precursor (GPC) of the virus surface glycoprotein that mediates virus receptor recognition and cell entry via endocytosis. The L genome RNA encodes the viral RNA-dependent RNA polymerase (RdRp, or L polymerase) and the small (ca 11 kDa) RING finger protein Z that has functions of a bona fide matrix protein including directing virus budding. Arenaviruses were thought to be relatively stable genetically with intra- and interspecies amino acid sequence identities of 90-95 % and 44-63 %, respectively. However, recent evidence has documented extensive arenavirus genetic variability in the field. Moreover, dramatic phenotypic differences have been documented among closely related LCMV isolates. These data provide strong evidence of viral quasispecies involvement in arenavirus adaptability and pathogenesis. Here, we will review several aspects of the molecular biology of arenaviruses, phylogeny and evolution, and quasispecies dynamics of arenavirus populations for a better understanding of arenavirus pathogenesis, as well as for the development of novel antiviral strategies to combat arenavirus infections.

Trends in Antiviral Strategies

Viral populations are true moving targets regarding the genomic sequences to be targeted in antiviral designs. Experts from different fields have expressed the need of new paradigms for antiviral interventions and viral disease control. This chapter reviews several strategies that aim at counteracting the adaptive capacity of viral quasispecies. The proposed designs are based on combinations of different antiviral drugs and immune modulators, or in the administration of virus-specific mutagenic agents, in an approach termed lethal mutagenesis of viruses. It consists of decreasing viral fitness by an excess of mutations that render viral proteins sub-optimal or non-functional. Viral extinction by lethal mutagenesis involves several sequential, overlapping steps that recapitulate the major concepts of intra-population interactions and genetic information stability discussed in preceding chapters. Despite the magnitude of the challenge, the chapter closes with some optimistic prospects for an effective control of viruses displaying error-prone replication, based on the combined targeting of replication fidelity and the induction of the innate immune response.

Deep sequencing and phylogenetic analysis of variants resistant to interferon-based protease inhibitor therapy in chronic hepatitis induced by genotype 1b hepatitis C virus

Because of recent advances in deep sequencing technology, detailed analysis of hepatitis C virus (HCV) quasispecies and their dynamic changes in response to direct antiviral agents (DAAs) became possible, although the role of quasispecies is not fully understood. In this study, to clarify the evolution of viral quasispecies and the origin of drug-resistant mutations induced by interferon (IFN)-based protease inhibitor therapy, the nonstructural-3 (NS3) region of genotype 1b HCV in 34 chronic hepatitis patients treated with telaprevir (TVR)/pegylated interferon (PEG-IFN)/ribavirin (RBV) was subjected to a deep sequencing study coupled with phylogenetic analysis. Twenty-six patients (76.5%) achieved a sustained viral response (SVR), while 8 patients did not (non-SVR; 23.5%). When the complexity of the quasispecies was expressed as the mutation frequency or Shannon entropy value, a significant decrease in the IFNL3 (rs8099917) TT group and a marginal decrease in the SVR group were found soon (12 h) after the introduction of treatment, whereas there was no decrease in the non-SVR group and no significant decrease in mutation frequency in the IFNL3 TG/GG group. In the analysis of viral quasispecies composition in non-SVR patients, major populations greatly changed, accompanied by the appearance of resistance, and the compositions were unlikely to return to the pretreatment composition even after the end of therapy. Clinically TVR-resistant variants were observed in 5 non-SVR patients (5/8, 62.5%), all of which were suspected to have acquired resistance by mutations through phylogenetic analysis. In conclusion, results of the study have important implications for treatment response and outcome in interferon-based protease inhibitor therapy.

IMPORTANCE

In the host, hepatitis C virus (HCV) consists of a variety of populations (quasispecies), and it is supposed that dynamic changes in quasispecies are closely related to pathogenesis, although this is poorly understood. In this study, recently developed deep sequencing technology was introduced, and changes in quasispecies associated with telaprevir (TVR)/pegylated interferon (PEG-IFN)/ribavirin (RBV) triple therapy and their clinical significance were investigated extensively by phylogenetic tree analysis. Through this study, the associations among treatment response, changes in viral quasispecies complexity in the early stage of treatment, changes in the quasispecies composition, and origin of TVR-resistant variant HCV were elucidated.

Deep sequencing analysis of viral infection and evolution allows rapid and detailed characterization of viral mutant spectrum

Motivation: The study of RNA virus populations is a challenging task. Each population of RNA virus is composed of a collection of different, yet related genomes often referred to as mutant spectra or quasispecies. Virologists using deep sequencing technologies face major obstacles when studying virus population dynamics, both experimentally and in natural settings due to the relatively high error rates of these technologies and the lack of high performance pipelines. In order to overcome these hurdles we developed a computational pipeline, termed ViVan (Viral Variance Analysis). ViVan is a complete pipeline facilitating the identification, characterization and comparison of sequence variance in deep sequenced virus populations.

Results: Applying ViVan on deep sequenced data obtained from samples that were previously characterized by more classical approaches, we uncovered novel and potentially crucial aspects of virus populations. With our experimental work, we illustrate how ViVan can be used for studies ranging from the more practical, detection of resistant mutations and effects of antiviral treatments, to the more theoretical temporal characterization of the population in evolutionary studies.

Availability and implementation: Freely available on the web at http://www.vivanbioinfo.org

Contact: nshomron@post.tau.ac.il

Supplementary information:Supplementary data are available at Bioinformatics online.

Getting to Know Viral Evolutionary Strategies: Towards the Next Generation of Quasispecies Models

Viral populations are formed by complex ensembles of genomes with broad phenotypic diversity. The adaptive strategies deployed by these ensembles are multiple and often cannot be predicted a priori. Our understanding of viral dynamics is mostly based on two kinds of empirical approaches: one directed towards characterizing molecular changes underlying fitness changes and another focused on population-level responses. Simultaneously, theoretical efforts are directed towards developing a formal picture of viral evolution by means of more realistic fitness landscapes and reliable population dynamics models. New technologies, chiefly the use of next-generation sequencing and related tools, are opening avenues connecting the molecular and the population levels. In the near future, we hope to be witnesses of an integration of these still decoupled approaches, leading into more accurate and realistic quasispecies models able to capture robust generalities and endowed with a satisfactory predictive power.

Favipiravir elicits antiviral mutagenesis during virus replication in vivo

Lethal mutagenesis has emerged as a novel potential therapeutic approach to treat viral infections. Several studies have demonstrated that increases in the high mutation rates inherent to RNA viruses lead to viral extinction in cell culture, but evidence during infections in vivo is limited. In this study, we show that the broad-range antiviral nucleoside favipiravir reduces viral load in vivo by exerting antiviral mutagenesis in a mouse model for norovirus infection. Increased mutation frequencies were observed in samples from treated mice and were accompanied with lower or in some cases undetectable levels of infectious virus in faeces and tissues. Viral RNA isolated from treated animals showed reduced infectivity, a feature of populations approaching extinction during antiviral mutagenesis. These results suggest that favipiravir can induce norovirus mutagenesis in vivo, which in some cases leads to virus extinction, providing a proof-of-principle for the use of favipiravir derivatives or mutagenic nucleosides in the clinical treatment of noroviruses.

DOI:http://dx.doi.org/10.7554/eLife.03679.001

Viruses can infect, take control of and replicate themselves inside the living cells of other organisms. Some viral diseases can be treated with antiviral drugs, which stop viral infections either by making it more difficult for viruses to enter cells or by preventing the virus replicating once inside. As antiviral drugs are currently only available to treat a handful of viral infections, efforts are underway to develop and test experimental antiviral drugs.

One such experimental drug is called favipiravir, which is proving to be effective against several viruses that store their genetic information in the form of RNA molecules. These viruses include those that cause diseases such as influenza, gastroenteritis, and Ebola. Along with ongoing work determining how safe and effective favipiravir is for treating viral infections, researchers are also attempting to better understand how favipiravir works.

Whenever a strand of RNA is copied to allow a new virus to form, there is a risk that mistakes—or mutations—that could harm the virus are introduced into the genetic code. Previous experiments performed on cells grown in the laboratory suggested that favipiravir works against RNA viruses by increasing how often these mutations occur. RNA viruses naturally experience a large number of mutations and the ability to make mutations is in fact a benefit for viruses as it allows them to evolve rapidly and to escape immune responses. However, there is a limit to how many mutations can be tolerated in the viral genome before it can no longer replicate. Therefore, a slight increase in how often mutations occur—as thought to be caused by favipiravir—is able to stop the RNA virus replicating and halt the infection. However, favipiravir's mode of action had yet to be confirmed in living animals.

Using mice, Arias et al. tested favipiravir's ability to treat a persistent infection by norovirus—the most common cause of viral gastroenteritis in humans and also responsible for life-threatening chronic diarrhoea in immunodeficient patients. Treatment increased the number of mutations that occurred when the viral RNA replicated and could reduce the amount of virus in the mice to undetectable levels. In addition, favipiravir did not show toxicity in mice after 8 weeks of treatment. This suggests that favipiravir has the potential to be used safely and effectively to treat norovirus and other RNA viruses, although further studies are required before it can be developed into a clinical treatment.

DOI:http://dx.doi.org/10.7554/eLife.03679.002

Changes in protein domains outside the catalytic site of the bacteriophage Qβ replicase reduce the mutagenic effect of 5-azacytidine

The high genetic heterogeneity and great adaptability of RNA viruses are ultimately caused by the low replication fidelity of their polymerases. However, single amino acid substitutions that modify replication fidelity can evolve in response to mutagenic treatments with nucleoside analogues. Here, we investigated how two independent mutants of the bacteriophage Qβ replicase (Thr210Ala and Tyr410His) reduce sensitivity to the nucleoside analogue 5-azacytidine (AZC). Despite being located outside the catalytic site, both mutants reduced the mutation frequency in the presence of the drug. However, they did not modify the type of AZC-induced substitutions, which was mediated mainly by ambiguous base pairing of the analogue with purines. Furthermore, the Thr210Ala and Tyr410His substitutions had little or no effect on replication fidelity in untreated viruses. Also, both substitutions were costly in the absence of AZC or when the action of the drug was suppressed by adding an excess of natural pyrimidines (uridine or cytosine). Overall, the phenotypic properties of these two mutants were highly convergent, despite the mutations being located in different domains of the Qβ replicase. This suggests that treatment with a given nucleoside analogue tends to select for a unique functional response in the viral replicase.

IMPORTANCE

In the last years, artificial increase of the replication error rate has been proposed as an antiviral therapy. In this study, we investigated the mechanisms by which two substitutions in the Qβ replicase confer partial resistance to the mutagenic nucleoside analogue AZC. As opposed to previous work with animal viruses, where different mutations selected sequentially conferred nucleoside analogue resistance through different mechanisms, our results suggest that there are few or no alternative AZC resistance phenotypes in Qβ. Also, despite resistance mutations being highly costly in the absence of the drug, there was no sequential fixation of secondary mutations. Bacteriophage Qβ is the virus with the highest reported mutation rate, which should make it particularly sensitive to nucleoside analogue treatments, probably favoring resistance mutations even if they incur high costs. The results are also relevant for understanding the possible pathways by which fidelity of the replication machinery can be modified.

RNA virus evolution at variable error rate

ABSTRACT: 

RNA viruses replicate their genomes with very high error rates, which leads to the generation of a large genetic diversity that makes them highly adaptable to most environmental pressures, including antiviral drugs and immune responses. However, since most mutations are deleterious, an excess of errors can be very negative for RNA viruses, entailing that error rates must be finely regulated. Currently, the manipulation of the error rate is emerging as a promising antiviral therapy that could minimize the problem of virus adaptation to classical treatments. This review provides a detailed analysis of the different outcomes that can result from the variation of the error rate in RNA viruses, on the basis of the more relevant findings obtained in experimental studies.

Inference with viral quasispecies diversity indices: clonal and NGS approaches

Given the inherent dynamics of a viral quasispecies, we are often interested in the comparison of diversity indices of sequential samples of a patient, or in the comparison of diversity indices of virus in groups of patients in a treated versus control design. It is then important to make sure that the diversity measures from each sample may be compared with no bias and within a consistent statistical framework. In the present report, we review some indices often used as measures for viral quasispecies complexity and provide means for statistical inference, applying procedures taken from the ecology field. In particular, we examine the Shannon entropy and the mutation frequency, and we discuss the appropriateness of different normalization methods of the Shannon entropy found in the literature. By taking amplicons ultra-deep pyrosequencing (UDPS) raw data as a surrogate of a real hepatitis C virus viral population, we study through in-silico sampling the statistical properties of these indices under two methods of viral quasispecies sampling, classical cloning followed by Sanger sequencing (CCSS) and next-generation sequencing (NGS) such as UDPS. We propose solutions specific to each of the two sampling methods-CCSS and NGS-to guarantee statistically conforming conclusions as free of bias as possible.

CONTACT

josep.gregori@gmail.com Supplementary information: Supplementary data are available at Bioinformatics online.

Genetic variation in vitro and in vivo of an attenuated Lassa vaccine candidate

The attenuated Lassa vaccine candidate ML29 is a laboratory-produced reassortant between Lassa and Mopeia viruses, two Old World arenaviruses that differ by 40% in nucleic acid sequence. In our previous studies, ML29 elicited sterilizing immunity against Lassa virus challenge in guinea pigs and marmosets and virus-specific cell-mediated immunity in both simian immunodeficiency virus (SIV)-infected and uninfected rhesus macaques. Here, we show that ML29 is stable after 12 passages in vitro without losing its plaque morphology or its attenuated phenotype in suckling mice. Additionally, we used deep sequencing to characterize the viral population comprising the original stock of ML29, the stock of ML29 after 12 passages in Vero cells, and the ML29 isolates obtained from vaccinated animals. Twenty-seven isolates bore approximately 77 mutations that exceeded 20% of the single-nucleotide polymorphism (SNP) changes at any single locus. Of these 77 mutations, 5 appeared to be host specific, for example, appearing in mice but not in primates. None of these mutations were reversions of ML29 to the sequences of the parental Lassa and Mopeia viruses. The host-specific mutations indicate viral adaptations to virus-host interactions, and such interactions make reasonable targets for antiviral approaches. Variants capable of chronic infection did not emerge from any of the primate infections, even in immune-deficient animals, indicating that the ML29 reassortant is reasonably stable in vivo. In conclusion, the preclinical studies of ML29 as a Lassa virus vaccine candidate have been advanced, showing high levels of protection in nonhuman primates and acceptable stability both in vitro and in vivo.

Coronaviruses lacking exoribonuclease activity are susceptible to lethal mutagenesis: evidence for proofreading and potential therapeutics

No therapeutics or vaccines currently exist for human coronaviruses (HCoVs). The Severe Acute Respiratory Syndrome-associated coronavirus (SARS-CoV) epidemic in 2002–2003, and the recent emergence of Middle East Respiratory Syndrome coronavirus (MERS-CoV) in April 2012, emphasize the high probability of future zoonotic HCoV emergence causing severe and lethal human disease. Additionally, the resistance of SARS-CoV to ribavirin (RBV) demonstrates the need to define new targets for inhibition of CoV replication. CoVs express a 3′-to-5′ exoribonuclease in nonstructural protein 14 (nsp14-ExoN) that is required for high-fidelity replication and is conserved across the CoV family. All genetic and biochemical data support the hypothesis that nsp14-ExoN has an RNA proofreading function. Thus, we hypothesized that ExoN is responsible for CoV resistance to RNA mutagens. We demonstrate that while wild-type (ExoN+) CoVs were resistant to RBV and 5-fluorouracil (5-FU), CoVs lacking ExoN activity (ExoN−) were up to 300-fold more sensitive. While the primary antiviral activity of RBV against CoVs was not mutagenesis, ExoN− CoVs treated with 5-FU demonstrated both enhanced sensitivity during multi-cycle replication, as well as decreased specific infectivity, consistent with 5-FU functioning as a mutagen. Comparison of full-genome next-generation sequencing of 5-FU treated SARS-CoV populations revealed a 16-fold increase in the number of mutations within the ExoN− population as compared to ExoN+. Ninety percent of these mutations represented A:G and U:C transitions, consistent with 5-FU incorporation during RNA synthesis. Together our results constitute direct evidence that CoV ExoN activity provides a critical proofreading function during virus replication. Furthermore, these studies identify ExoN as the first viral protein distinct from the RdRp that determines the sensitivity of RNA viruses to mutagens. Finally, our results show the importance of ExoN as a target for inhibition, and suggest that small-molecule inhibitors of ExoN activity could be potential pan-CoV therapeutics in combination with RBV or RNA mutagens.

RNA viruses have high mutation rates (10⁻³ to 10⁻⁵ mutations/nucleotide/round of replication), allowing for rapid viral adaptation in response to selective pressure. While RNA viruses have long been considered unable to correct mistakes during replication, CoVs such as SARS-CoV and the recently emerged MERS-CoV are important exceptions to this paradigm. All CoVs encode an exoribonuclease activity in nonstructural protein 14 (nsp14-ExoN) that is proposed to prevent and/or remove misincorporated nucleotides. Because of the demonstrated resistance of SARS-CoV to the antiviral drug ribavirin (RBV), we hypothesized that ExoN is responsible for CoV resistance to RNA mutagens. Using RBV and the RNA mutagen 5-fluorouracil (5-FU), we show that CoVs lacking ExoN activity (ExoN−) are highly susceptible to RBV and 5-FU, in contrast to wild-type (ExoN+) CoVs. The inhibitory activity of 5-FU against ExoN− viruses resulted specifically from 5-FU incorporation during viral RNA synthesis that lead to extensive mutagenesis within the viral population, and was associated with a profound decrease in virus specific infectivity. These results demonstrate the proofreading activity of ExoN during virus replication and suggest that inhibitors of ExoN activity could be broadly useful inhibitors of CoV replication in combination with RBV or RNA mutagens.

Molecular dissection of a viral quasispecies under mutagenic treatment: positive correlation between fitness loss and mutational load

Low fidelity replication and the absence of error-repair activities in RNA viruses result in complex and adaptable ensembles of related genomes in the viral population, termed quasispecies, with important implications for natural infections. Theoretical predictions suggested that elevated replication error rates in RNA viruses might be near to a maximum compatible with viral viability. This fact encouraged the use of mutagenic nucleosides as a new antiviral strategy to induce viral extinction through increased replication error rates. Despite extensive evidence of lethal mutagenesis of RNA viruses by different mutagenic compounds, a detailed picture of the infectivity of individual genomes and its relationship with the mutations accumulated is lacking. Here, we report a molecular analysis of a foot-and-mouth disease virus population previously subjected to heavy mutagenesis to determine whether a correlation between increased mutagenesis and decreased fitness existed. Plaque-purified viruses isolated from a ribavirin-treated quasispecies presented decreases of up to 200-fold in infectivity relative to clones in the reference population, associated with an overall eightfold increase in the mutation frequency. This observation suggests that individual infectious genomes of a quasispecies subjected to increased mutagenesis lose infectivity by their continuous mutagenic ‘poisoning’. These results support the lethal defection model of virus extinction and the practical use of chemical mutagens as antiviral treatment. Even when extinction is not achieved, mutagenesis can decrease the infectivity of surviving virus, and facilitate their clearance by host immune responses or complementing antiviral approaches.

Arenavirus variations due to host-specific adaptation

Arenavirus particles are enveloped and contain two single-strand RNA genomic segments with ambisense coding. Genetic plasticity of the arenaviruses comes from transcription errors, segment reassortment, and permissive genomic packaging, and results in their remarkable ability, as a group, to infect a wide variety of hosts. In this review, we discuss some in vitro studies of virus genetic and phenotypic variation after exposure to selective pressures such as high viral dose, mutagens and antivirals. Additionally, we discuss the variation in vivo of selected isolates of Old World arenaviruses, particularly after infection of different animal species. We also discuss the recent emergence of new arenaviruses in the context of our observations of sequence variations that appear to be host-specific.

Arenaviruses and lethal mutagenesis. Prospects for new ribavirin-based interventions

Lymphocytic choriomeningitis virus (LCMV) has contributed to unveil some of the molecular mechanisms of lethal mutagenesis, or loss of virus infectivity due to increased mutation rates. Here we review these developments, and provide additional evidence that ribavirin displays a dual mutagenic and inhibitory activity on LCMV that can be relevant to treatment designs. Using 5-fluorouracil as mutagenic agent and ribavirin either as inhibitor or mutagen, we document an advantage of a sequential inhibitor-mutagen administration over the corresponding combination treatment to achieve a low LCMV load in cell culture. This advantage is accentuated in the concentration range in which ribavirin acts mainly as an inhibitor, rather than as mutagen. This observation reinforces previous theoretical and experimental studies in supporting a sequential inhibitor-mutagen administration as a possible antiviral design. Given recent progress in the development of new inhibitors of arenavirus replication, our results suggest new options of ribavirin-based anti-arenavirus treatments.

Hyperediting of human T-cell leukemia virus type 2 and simian T-cell leukemia virus type 3 by the dsRNA adenosine deaminase ADAR-1

RNA editing mediated by adenosine deaminases acting on RNA (ADARs) converts adenosine (A) to inosine (I) residues in dsRNA templates. While ADAR-1-mediated editing was essentially described for RNA viruses, the present work addresses the issue for two δ-retroviruses, human T-cell leukemia virus type 2 and simian T-cell leukemia virus type 3 (HTLV-2 and STLV-3). We examined whether ADAR-1 could edit HTLV-2 and STLV-3 virus genomes in cell culture and in vivo. Using a highly sensitive PCR-based method, referred to as 3DI-PCR, we showed that ADAR-1 could hypermutate adenosine residues in HTLV-2. STLV-3 hypermutation was obtained without using 3DI-PCR, suggesting a higher mutation frequency for this virus. Detailed analysis of the dinucleotide editing context showed preferences for 5' ArA and 5' UrA. In conclusion, the present observations demonstrate that ADAR-1 massively edits HTLV-2 and STLV-3 retroviruses in vitro, but probably remains a rare phenomenon in vivo.

The impact of quasispecies dynamics on the use of therapeutics

The application of quasispecies theory to viral populations has boosted our understanding of how endogenous and exogenous features condition their adaptation. Mounting empirical evidence demonstrates that internal interactions within mutant spectra may cause unexpected responses to antiviral treatments. In this scenario, increased mutagenesis could be efficient at low mutagen doses due to the lethal action of defective genomes, whereas sequential administration of antiviral drugs might be superior to combination therapies. Our ability to predict the outcome of a particular therapy takes advantage of the complementary use of in vivo observations, in vitro experiments, and mathematical models.

Stochastic simulations suggest that HIV-1 survives close to its error threshold

The use of mutagenic drugs to drive HIV-1 past its error threshold presents a novel intervention strategy, as suggested by the quasispecies theory, that may be less susceptible to failure via viral mutation-induced emergence of drug resistance than current strategies. The error threshold of HIV-1, , however, is not known. Application of the quasispecies theory to determine poses significant challenges: Whereas the quasispecies theory considers the asexual reproduction of an infinitely large population of haploid individuals, HIV-1 is diploid, undergoes recombination, and is estimated to have a small effective population size in vivo. We performed population genetics-based stochastic simulations of the within-host evolution of HIV-1 and estimated the structure of the HIV-1 quasispecies and . We found that with small mutation rates, the quasispecies was dominated by genomes with few mutations. Upon increasing the mutation rate, a sharp error catastrophe occurred where the quasispecies became delocalized in sequence space. Using parameter values that quantitatively captured data of viral diversification in HIV-1 patients, we estimated to be substitutions/site/replication, ∼2–6 fold higher than the natural mutation rate of HIV-1, suggesting that HIV-1 survives close to its error threshold and may be readily susceptible to mutagenic drugs. The latter estimate was weakly dependent on the within-host effective population size of HIV-1. With large population sizes and in the absence of recombination, our simulations converged to the quasispecies theory, bridging the gap between quasispecies theory and population genetics-based approaches to describing HIV-1 evolution. Further, increased with the recombination rate, rendering HIV-1 less susceptible to error catastrophe, thus elucidating an added benefit of recombination to HIV-1. Our estimate of may serve as a quantitative guideline for the use of mutagenic drugs against HIV-1.

Currently available antiretroviral drugs curtail HIV infection but fail to eradicate the virus. A strategy of intervention radically different from that employed by current drugs has been proposed by the molecular quasispecies theory. The theory predicts that increasing the viral mutation rate beyond a critical value, called the error threshold, would cause a severe loss of genetic information, potentially leading to viral clearance. Several chemical mutagens are now being developed that can increase the mutation rate of HIV-1. Their success depends on reliable estimates of the error threshold of HIV-1, which are currently lacking. The quasispecies theory cannot be applied directly to HIV-1: the theory considers an infinitely large population of asexually reproducing haploid individuals, whereas HIV-1 is diploid, undergoes recombination, and is estimated to have a small effective population size in vivo. We employed detailed stochastic simulations that overcome the limitations of the quasispecies theory and accurately mimic HIV-1 evolution in vivo. With these simulations, we estimated the error threshold of HIV-1 to be ∼2–6-fold higher than its natural mutation rate, suggesting that HIV-1 survives close to its error threshold and may be readily susceptible to mutagenic drugs.

Mutation induced extinction in finite populations: lethal mutagenesis and lethal isolation

Reproduction is inherently risky, in part because genomic replication can introduce new mutations that are usually deleterious toward fitness. This risk is especially severe for organisms whose genomes replicate "semi-conservatively," e.g. viruses and bacteria, where no master copy of the genome is preserved. Lethal mutagenesis refers to extinction of populations due to an unbearably high mutation rate (U), and is important both theoretically and clinically, where drugs can extinguish pathogens by increasing their mutation rate. Previous theoretical models of lethal mutagenesis assume infinite population size (N). However, in addition to high U, small N can accelerate extinction by strengthening genetic drift and relaxing selection. Here, we examine how the time until extinction depends jointly on N and U. We first analytically compute the mean time until extinction (τ) in a simplistic model where all mutations are either lethal or neutral. The solution motivates the definition of two distinct regimes: a survival phase and an extinction phase, which differ dramatically in both how τ scales with N and in the coefficient of variation in time until extinction. Next, we perform stochastic population-genetics simulations on a realistic fitness landscape that both (i) features an epistatic distribution of fitness effects that agrees with experimental data on viruses and (ii) is based on the biophysics of protein folding. More specifically, we assume that mutations inflict fitness penalties proportional to the extent that they unfold proteins. We find that decreasing N can cause phase transition-like behavior from survival to extinction, which motivates the concept of "lethal isolation." Furthermore, we find that lethal mutagenesis and lethal isolation interact synergistically, which may have clinical implications for treating infections. Broadly, we conclude that stably folded proteins are only possible in ecological settings that support sufficiently large populations.

Nucleoside analogue mutagenesis of a single-stranded DNA virus: evolution and resistance

It has been well established that chemical mutagenesis has adverse fitness effects in RNA viruses, often leading to population extinction. This is mainly a consequence of the high RNA virus spontaneous mutation rates, which situate them close to the extinction threshold. Single-stranded DNA viruses are the fastest-mutating DNA-based systems, with per-nucleotide mutation rates close to those of some RNA viruses, but chemical mutagenesis has been much less studied in this type of viruses. Here, we serially passaged bacteriophage X174 in the presence of the nucleoside analogue 5-fluorouracil (5-FU). We found that 5-FU was unable to trigger population extinction for the range of concentrations tested, but it negatively affected viral adaptability. The phage evolved partial drug resistance, and parallel nucleotide substitutions appearing in independently evolved lines were identified as candidate resistance mutations. Using site-directed mutagenesis, two single-nucleotide substitutions in the lysis protein E (T572C and A781G) were shown to be selectively advantageous in the presence of 5-FU. In RNA viruses, base analogue resistance is often mediated by changes in the viral polymerase, but this mechanism is not possible for X174 and other single-stranded DNA viruses because they do not encode their own polymerase. In addition to increasing mutation rates, 5-FU produces a wide variety of cytotoxic effects at the levels of replication, transcription, and translation. We found that substitutions T572C and A781G lost their ability to confer 5-FU resistance after cells were supplemented with deoxythymidine, suggesting that their mechanism of action is at the DNA level. We hypothesize that regulation of lysis time may allow the virus to optimize progeny size in cells showing defects in DNA synthesis.

Viral quasispecies evolution

Evolution of RNA viruses occurs through disequilibria of collections of closely related mutant spectra or mutant clouds termed viral quasispecies. Here we review the origin of the quasispecies concept and some biological implications of quasispecies dynamics. Two main aspects are addressed: (i) mutant clouds as reservoirs of phenotypic variants for virus adaptability and (ii) the internal interactions that are established within mutant spectra that render a virus ensemble the unit of selection. The understanding of viruses as quasispecies has led to new antiviral designs, such as lethal mutagenesis, whose aim is to drive viruses toward low fitness values with limited chances of fitness recovery. The impact of quasispecies for three salient human pathogens, human immunodeficiency virus and the hepatitis B and C viruses, is reviewed, with emphasis on antiviral treatment strategies. Finally, extensions of quasispecies to nonviral systems are briefly mentioned to emphasize the broad applicability of quasispecies theory.

Biomedical implications of viral mutation and evolution

Mutation rates vary hugely across viruses and strongly determine their evolution. In addition, viral mutation and evolution are biomedically relevant because they can determine pathogenesis, vaccine efficacy and antiviral resistance. We review experimental methods for estimating viral mutation rates and how these estimates vary across viral groups, paying special attention to the more general trends. Recent advances positing a direct association between viral mutation rates and virulence, or the use of high-fidelity variants as attenuated vaccines, are also discussed. Finally, we review the implications of viral mutation and evolution for the design of rational antiviral therapies and for efficient epidemiological surveillance.

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.

Progress in the experimental therapy of severe arenaviral infections

A number of viruses in the family Arenaviridae cause severe illness in humans. Lassa virus in West Africa and a number of agents in South America produce hemorrhagic fever in persons exposed to aerosolized excretions of the pathogens' rodent hosts. Because arenaviruses are not transmitted by arthropods, and person-to-person spread is rare, human infections occur singly and sporadically, and are usually not diagnosed until the patient is severely ill. Because the arenaviruses are naturally transmitted by the airborne route, they also pose a potential threat as aerosolized bioterror weapons. The broad-spectrum antiviral drug ribavirin was shown to reduce mortality from Lassa fever, and has been tested against Argentine hemorrhagic fever, but it is not an approved treatment for either disease. Human immune convalescent plasma was proven to be effective for Argentine hemorrhagic fever in a controlled trial. New treatments are needed to block viral replication without causing toxicity and to prevent the increased vascular permeability that is responsible for hypotension and shock. In this paper, we review current developments in the experimental therapy of severe arenaviral infections, focusing on drugs that have been tested in animal models, and provide a perspective on future research.

Mathematical Modeling of Viral Zoonoses in Wildlife

Zoonoses are a worldwide public health concern, accounting for approximately 75% of human infectious diseases. In addition, zoonoses adversely affect agricultural production and wildlife. We review some mathematical models developed for the study of viral zoonoses in wildlife and identify areas where further modeling efforts are needed.

Mutagen-mediated enhancement of HIV-1 replication in persistently infected cells

Lethal mutagenesis, a new antiviral strategy to extinguish virus through elevated mutation rates, was explored in H61-D cells an HIV-1 persistently infected lymphoid cell line. Three mutagenic agents: 5-hydroxy-2(')-deoxycytidine (5-OHdC), 5-fluorouracil (5-FU) and 2,2(')-difluoro-2(')-deoxycytidine (gemcitabine) were used. After 54 passages, treatments with 5-FU and gemcitabine reduced virus infectivity, p24 and RT activity. Treatment with the pyrimidine analog 5-OHdC resulted in increases of p24 production, RT activity and infectivity. Rise in viral replication by 5-OHdC during HIV-1 persistence is in contrast with its inhibitory effect in acute infections. Viral replication enhancement by 5-OHdC was associated with an increase in intracellular HIV-1 RNA mutations. Mechanisms of HIV-1 replication enhancement by 5-OHdC are unknown but some potential factors are discussed. Increase of HIV-1 replication by 5-OHdC cautions against the use, without previous analyses, of mutagenic nucleoside analogs for AIDS treatment.

ADARs: viruses and innate immunity

Double-stranded RNA (dsRNA) functions both as a substrate of ADARs and also as a molecular trigger of innate immune responses. ADARs, adenosine deaminases that act on RNA, catalyze the deamination of adenosine (A) to produce inosine (I) in dsRNA. ADARs thereby can destablize RNA structures, because the generated I:U mismatch pairs are less stable than A:U base pairs. Additionally, I is read as G instead of A by ribosomes during translation and by viral RNA-dependent RNA polymerases during RNA replication. Members of several virus families have the capacity to produce dsRNA during viral genome transcription and replication. Sequence changes (A-G, and U-C) characteristic of A-I editing can occur during virus growth and persistence. Foreign viral dsRNA also mediates both the induction and the action of interferons. In this chapter our current understanding of the role and significance of ADARs in the context of innate immunity, and as determinants of the outcome of viral infection, will be considered.

Quasispecies Dynamics of RNA Viruses

RNA viruses, such as human immunodeficiency virus, hepatitis C virus, influenza virus, and poliovirus replicate with very high mutation rates and exhibit very high genetic diversity. The extremely high genetic diversity of RNA virus populations originates that they replicate as complex mutant spectra known as viral quasispecies. The quasispecies dynamics of RNA viruses are closely related to viral pathogenesis and disease, and antiviral treatment strategies. Over the past several decades, the quasispecies concept has been expanded to provide an adequate framework to explain complex behavior of RNA virus populations. Recently, the quasispecies concept has been used to study other complex biological systems, such as tumor cells, bacteria, and prions. Here, we focus on some questions regarding viral and theoretical quasispecies concepts, as well as more practical aspects connected to pathogenesis and resistance to antiviral treatments. A better knowledge of virus diversification and evolution may be critical in preventing and treating the spread of pathogenic viruses.

5-fluorouracil in lethal mutagenesis of foot-and-mouth disease virus

5-fluorouracil (FU) is a pyrimidine analogue extensively used in cancer chemotherapy. FU can be metabolized into 5-fluorouridine-triphosphate, which can be used as substrate for viral RNA-dependent RNA polymerases. This results in the incorporation of mutations into viral RNA. Accumulation of mutations may lead to loss of virus infectivity, in a process known as lethal mutagenesis. RNA virus pathogens are particularly difficult to control because they are highly mutable, and mutants resistant to antiviral agents are readily selected. Here, we review the basic principles of lethal mutagenesis as an antiviral approach, and the participation of FU in its development. Recent studies with foot-and-mouth disease virus indicate that FU can act both as an inhibitor and as a mutagen during foot-and-mouth disease virus replication. This dual activity renders FU an adequate drug for lethal mutagenesis. We suggest that structural and biochemical studies can contribute to the lead to new design of base or nucleoside analogues targeted specifically to viral polymerases.

Ribavirin can be mutagenic for arenaviruses

Arenaviruses include several important human pathogens, and there are very limited options of preventive or therapeutic interventions to combat these viruses. An off-label use of the purine nucleoside analogue ribavirin (1-β-d-ribofuranosyl-1-H-1,2,4-triazole-3-carboxamide) is the only antiviral treatment currently available for arenavirus infections. However, the ribavirin antiviral mechanism action against arenaviruses remains unknown. Here we document that ribavirin is mutagenic for the prototypic arenavirus lymphocytic choriomeningitis virus (LCMV) in cell culture. The mutagenic activity of ribavirin on LCMV was observed under single- and multiple-passage regimes and could not be accounted for by a decrease of the intracellular GTP pool promoted by ribavirin-mediated inhibition of inosine monophosphate dehydrogenase (IMPDH). Our findings suggest that the antiviral activity of ribavirin on arenaviruses might be exerted, at least partially, by lethal mutagenesis. Implications for antiarenavirus therapy are discussed.