The impact of quasispecies dynamics on the use of therapeutics

Molecular mechanism of rhinovirus escape from the Pyrazolo[3,4-d]pyrimidine capsid-binding inhibitor OBR-5-340 via mutations distant from the binding pocket: Derivatives that brake resistance

Rhinoviruses (RVs) cause the common cold. Attempts at discovering small molecule inhibitors have mainly concentrated on compounds supplanting the medium chain fatty acids residing in the sixty icosahedral symmetry-related hydrophobic pockets of the viral capsid of the Rhinovirus-A and -B species. High-affinity binding to these pockets stabilizes the capsid against structural changes necessary for the release of the ss(+) RNA genome into the cytosol of the host cell. However, single-point mutations may abolish this binding. RV-B5 is one of several RVs that are naturally resistant against the well-established antiviral agent pleconaril. However, RV-B5 is strongly inhibited by the pyrazolopyrimidine OBR-5-340. Here, we report on isolation and characterization of RV-B5 mutants escaping OBR-5-340 inhibition and show that substitution of amino acid residues not only within the binding pocket but also remote from the binding pocket hamper inhibition. Molecular dynamics network analysis revealed that strong inhibition occurs when an ensemble of several sequence stretches of the capsid proteins enveloping OBR-5-340 move together with OBR-5-340. Mutations abrogating this dynamic, regardless of whether being localized within the binding pocket or distant from it result in escape from inhibition. Pyrazolo [3,4-d]pyrimidine derivatives overcoming OBR-5-340 escape of various RV-B5 mutants were identified. Our work contributes to the understanding of the properties of capsid-binding inhibitors necessary for potent and broad-spectrum inhibition of RVs.

Quasispecies Nature of RNA Viruses: Lessons from the Past

Viral quasispecies are distinct but closely related mutants formed by the disparity in viral genomes due to recombination, mutations, competition, and selection pressure. Theoretical derivation for the origin of a quasispecies is owed to the error-prone replication by polymerase and mutants of RNA replicators. Here, we briefly addressed the theoretical and mathematical origin of quasispecies and their dynamics. The impact of quasispecies for major salient human pathogens is reviewed. In the current global scenario, rapid changes in geographical landscapes favor the origin and selection of mutants. It comes as no surprise that a cauldron of mutants poses a significant risk to public health, capable of causing pandemics. Mutation rates in RNA viruses are magnitudes higher than in DNA organisms, explaining their enhanced virulence and evolvability. RNA viruses cause the most devastating pandemics; for example, members of the Orthomyxoviridae family caused the great influenza pandemic (1918 flu or Spanish flu), the SARS (severe acute respiratory syndrome) and MERS (Middle East respiratory syndrome) outbreak, and the human immunodeficiency viruses (HIV), lentiviruses of the Retroviridae family, caused worldwide devastation. Rapidly evolving RNA virus populations are a daunting challenge for the designing of effective control measures like vaccines. Developing awareness of the evolutionary dispositions of RNA viral mutant spectra and what influences their adaptation and virulence will help curtail outbreaks of past and future pathogens.

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.

Viral Fitness, Population Complexity, Host Interactions, and Resistance to Antiviral Agents

Fitness of viruses has become a standard parameter to quantify their adaptation to a biological environment. Fitness determinations for RNA viruses (and some highly variable DNA viruses) meet with several uncertainties. Of particular interest are those that arise from mutant spectrum complexity, absence of population equilibrium, and internal interactions among components of a mutant spectrum. Here, concepts, fitness measurements, limitations, and current views on experimental viral fitness landscapes are discussed. The effect of viral fitness on resistance to antiviral agents is covered in some detail since it constitutes a widespread problem in antiviral pharmacology, and a challenge for the design of effective antiviral treatments. Recent evidence with hepatitis C virus suggests the operation of mechanisms of antiviral resistance additional to the standard selection of drug-escape mutants. The possibility that high replicative fitness may be the driver of such alternative mechanisms is considered. New broad-spectrum antiviral designs that target viral fitness may curtail the impact of drug-escape mutants in treatment failures. We consider to what extent fitness-related concepts apply to coronaviruses and how they may affect strategies for COVID-19 prevention and treatment.

The Intrapopulation Genetic Diversity of RNA Virus May Influence the Sensitivity of Chlorine Disinfection

RNA virus populations are not clonal; rather, they comprise a mutant swarm in which sequences are slightly different from the master sequence. Genetic diversity within a population (intrapopulation genetic diversity) is critical for RNA viruses to survive under environmental stresses. Disinfection has become an important practice in the control of pathogenic viruses; however, the impact of intrapopulation genetic diversity on the sensitivity of disinfection, defined as -log10 (postdisinfected infectious titer/predisinfected titer), has not been elucidated. In this study, we serially passaged populations of rhesus rotavirus. We demonstrated that populations with reduced chlorine sensitivity emerged at random and independently of chlorine exposure. Sequencing analysis revealed that compared with sensitive populations, less-sensitive ones had higher non-synonymous genetic diversity of the outer capsid protein gene, suggesting that changes in the amino acid sequences of the outer capsid protein were the main factors influencing chlorine sensitivity. No common mutations were found among less-sensitive populations, indicating that rather than specific mutations, the diversity of the outer capsid protein itself was associated with the disinfection sensitivity and that the disinfection sensitivity changed stochastically. Simulation results suggest that the disinfection sensitivity of a genetically diverse population is destabilized if cooperative viral clusters including multiple sequences are formed. These results advocate that any prevention measures leading to low intrapopulation genetic diversity are important to prevent the spread and evolution of pathogenic RNA viruses in society.

Population Disequilibrium as Promoter of Adaptive Explorations in Hepatitis C Virus

Replication of RNA viruses is characterized by exploration of sequence space which facilitates their adaptation to changing environments. It is generally accepted that such exploration takes place mainly in response to positive selection, and that further diversification is boosted by modifications of virus population size, particularly bottleneck events. Our recent results with hepatitis C virus (HCV) have shown that the expansion in sequence space of a viral clone continues despite prolonged replication in a stable cell culture environment. Diagnosis of the expansion was based on the quantification of diversity indices, the occurrence of intra-population mutational waves (variations in mutant frequencies), and greater individual residue variations in mutant spectra than those anticipated from sequence alignments in data banks. In the present report, we review our previous results, and show additionally that mutational waves in amplicons from the NS5A-NS5B-coding region are equally prominent during HCV passage in the absence or presence of the mutagenic nucleotide analogues favipiravir or ribavirin. In addition, by extending our previous analysis to amplicons of the NS3- and NS5A-coding region, we provide further evidence of the incongruence between amino acid conservation scores in mutant spectra from infected patients and in the Los Alamos National Laboratory HCV data banks. We hypothesize that these observations have as a common origin a permanent state of HCV population disequilibrium even upon extensive viral replication in the absence of external selective constraints or changes in population size. Such a persistent disequilibrium—revealed by the changing composition of the mutant spectrum—may facilitate finding alternative mutational pathways for HCV antiviral resistance. The possible significance of our model for other genetically variable viruses is discussed.

Rare haplotype load as marker for lethal mutagenesis

RNA viruses replicate with a template-copying fidelity, which lies close to an extinction threshold. Increases of mutation rate by nucleotide analogues can drive viruses towards extinction. This transition is the basis of an antiviral strategy termed lethal mutagenesis. We have introduced a new diversity index, the rare haplotype load (RHL), to describe NS5B (polymerase) mutant spectra of hepatitis C virus (HCV) populations passaged in absence or presence of the mutagenic agents favipiravir or ribavirin. The increase in RHL is more prominent in mutant spectra whose expansions were due to nucleotide analogues than to multiple passages in absence of mutagens. Statistical tests for paired mutagenized versus non-mutagenized samples with 14 diversity indices show that RHL provides consistently the highest standardized effect of mutagenic treatment difference for ribavirin and favipiravir. The results indicate that the enrichment of viral quasispecies in very low frequency minority genomes can serve as a robust marker for lethal mutagenesis. The diagnostic value of RHL from deep sequencing data is relevant to experimental studies on enhanced mutagenesis of viruses, and to pharmacological evaluations of inhibitors suspected to have a mutagenic activity.

Differential Shape of Geminivirus Mutant Spectra Across Cultivated and Wild Hosts With Invariant Viral Consensus Sequences

Geminiviruses (family Geminiviridae) possess single-stranded circular DNA genomes that are replicated by cellular polymerases in plant host cell nuclei. In their hosts, geminivirus populations behave as ensembles of mutant and recombinant genomes, known as viral quasispecies. This favors the emergence of new geminiviruses with altered host range, facilitating new or more severe diseases or overcoming resistance traits. In warm and temperate areas several whitefly-transmitted geminiviruses of the genus Begomovirus cause the tomato yellow leaf curl disease (TYLCD) with significant economic consequences. TYLCD is frequently controlled in commercial tomatoes by using the dominant Ty-1 resistance gene. Over a 45 day period we have studied the diversification of three begomoviruses causing TYLCD: tomato yellow leaf curl virus (TYLCV), tomato yellow leaf curl Sardinia virus (TYLCSV) and tomato yellow leaf curl Malaga virus (TYLCMaV, a natural recombinant between TYLCV and TYLCSV). Viral quasispecies resulting from inoculation of geminivirus infectious clones were examined in plants of susceptible tomato (ty-1/ty-1), heterozygous resistant tomato (Ty-1/ty-1), common bean, and the wild reservoir Solanum nigrum. Differences in virus fitness across hosts were observed while viral consensus sequences remained invariant. However, the complexity and heterogeneity of the quasispecies were high, especially in common bean and the wild host. Interestingly, the presence or absence of the Ty-1 allele in tomato did not lead to differences in begomovirus mutant spectra. However, the fitness decrease of TYLCSV and TYLCV in tomato at 45 dpi might be related to an increase in CP (Coat protein) mutation frequency. In Solanum nigrum the recombinant TYLCMaV, which showed lower fitness than TYLCSV, at 45 dpi actively explored Rep (Replication associated protein) ORF but not the overlapping C4. Our results underline the importance of begomovirus mutant spectra during infections. This is especially relevant in the wild reservoir of the viruses, which has the potential to maintain highly diverse mutant spectra without modifying their consensus sequences.

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.

Extinction of West Nile Virus by Favipiravir through Lethal Mutagenesis

Favipiravir is an antiviral agent effective against several RNA viruses. The drug has been shown to protect mice against experimental infection with a lethal dose of West Nile virus (WNV), a mosquito-borne flavivirus responsible for outbreaks of meningitis and encephalitis for which no antiviral therapy has been licensed; however, the mechanism of action of the drug is still not well understood. Here, we describe the potent in vitro antiviral activity of favipiravir against WNV, showing that it decreases virus-specific infectivity and drives the virus to extinction. Two passages of WNV in the presence of 1 mM favipiravir-a concentration that is more than 10-fold lower than its 50% cytotoxic concentration (CC₅₀)-resulted in a significant increase in mutation frequency in the mutant spectrum and in a bias toward A→G and G→A transitions relative to the population passaged in the absence of the drug. These data, together with the fact that the drug is already licensed in Japan against influenza virus and in a clinical trial against Ebola virus, point to favipiravir as a promising antiviral agent to fight medically relevant flaviviral infections, such as that caused by WNV.

Quasispecies of Hepatitis C Virus Participate in Cell-Specific Infectivity

It is well documented that a variety of viral quasispecies are found in the patients with chronic infection of hepatitis C virus (HCV). However, the significance of quasispecies in the specific infectivity to individual cell types remains unknown. In the present study, we analyzed the role of quasispecies of the genotype 2a clone, JFH1 (HCVcc), in specific infectivity to the hepatic cell lines, Huh7.5.1 and Hep3B. HCV RNA was electroporated into Huh7.5.1 cells and Hep3B/miR-122 cells expressing miR-122 at a high level. Then, we adapted the viruses to Huh7 and Hep3B/miR-122 cells by serial passages and termed the resulting viruses HCVcc/Huh7 and HCVcc/Hep3B, respectively. Interestingly, a higher viral load was obtained in the homologous combination of HCVcc/Huh7 in Huh7.5.1 cells or HCVcc/Hep3B in Hep3B/miR-122 cells compared with the heterologous combination. By using a reverse genetics system and deep sequence analysis, we identified several adaptive mutations involved in the high affinity for each cell line, suggesting that quasispecies of HCV participate in cell-specific infectivity.

Collective Infectious Units in Viruses

Increasing evidence indicates that viruses do not simply propagate as independent virions among cells, organs, and hosts. Instead, viral spread is often mediated by structures that simultaneously transport groups of viral genomes, such as polyploid virions, aggregates of virions, virion-containing proteinaceous structures, secreted lipid vesicles, and virus-induced cell-cell contacts. These structures increase the multiplicity of infection, independently of viral population density and transmission bottlenecks. Collective infectious units may contribute to the maintenance of viral genetic diversity, and could have implications for the evolution of social-like virus-virus interactions. These may include various forms of cooperation such as immunity evasion, genetic complementation, division of labor, and relaxation of fitness trade-offs, but also noncooperative interactions such as negative dominance and interference, potentially leading to conflict.

An Efficient Microarray-Based Genotyping Platform for the Identification of Drug-Resistance Mutations in Majority and Minority Subpopulations of HIV-1 Quasispecies

The response of human immunodeficiency virus type 1 (HIV-1) quasispecies to antiretroviral therapy is influenced by the ensemble of mutants that composes the evolving population. Low-abundance subpopulations within HIV-1 quasispecies may determine the viral response to the administered drug combinations. However, routine sequencing assays available to clinical laboratories do not recognize HIV-1 minority variants representing less than 25% of the population. Although several alternative and more sensitive genotyping techniques have been developed, including next-generation sequencing (NGS) methods, they are usually very time consuming, expensive and require highly trained personnel, thus becoming unrealistic approaches in daily clinical practice. Here we describe the development and testing of a HIV-1 genotyping DNA microarray that detects and quantifies, in majority and minority viral subpopulations, relevant mutations and amino acid insertions in 42 codons of the pol gene associated with drug- and multidrug-resistance to protease (PR) and reverse transcriptase (RT) inhibitors. A customized bioinformatics protocol has been implemented to analyze the microarray hybridization data by including a new normalization procedure and a stepwise filtering algorithm, which resulted in the highly accurate (96.33%) detection of positive/negative signals. This microarray has been tested with 57 subtype B HIV-1 clinical samples extracted from multi-treated patients, showing an overall identification of 95.53% and 89.24% of the queried PR and RT codons, respectively, and enough sensitivity to detect minority subpopulations representing as low as 5–10% of the total quasispecies. The developed genotyping platform represents an efficient diagnostic and prognostic tool useful to personalize antiviral treatments in clinical practice.

Cooperative Interaction Within RNA Virus Mutant Spectra

RNA viruses usually consist of mutant spectra because of high error rates of viral RNA polymerases. Growth competition occurs among different viral variants, and the fittest clones predominate under given conditions. Individual variants, however, may not be entirely independent of each other, and internal interactions within mutant spectra can occur. Examples of cooperative and interfering interactions that exert enhancing and suppressing effects on replication of the wild-type virus, respectively, have been described, but their underlying mechanisms have not been well defined. It was recently found that the cooperation between wild-type and variant measles virus genomes produces a new phenotype through the heterooligomer formation of a viral protein. This observation provides a molecular mechanism underlying cooperative interactions within mutant spectra. Careful attention to individual sequences, in addition to consensus sequences, may disclose further examples of internal interactions within mutant spectra.

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.

Resistance of Hepatitis C Virus to Inhibitors: Complexity and Clinical Implications

Selection of inhibitor-resistant viral mutants is universal for viruses that display quasi-species dynamics, and hepatitis C virus (HCV) is no exception. Here we review recent results on drug resistance in HCV, with emphasis on resistance to the newly-developed, directly-acting antiviral agents, as they are increasingly employed in the clinic. We put the experimental observations in the context of quasi-species dynamics, in particular what the genetic and phenotypic barriers to resistance mean in terms of exploration of sequence space while HCV replicates in the liver of infected patients or in cell culture. Strategies to diminish the probability of viral breakthrough during treatment are briefly outlined.

Increased intrahepatic quasispecies heterogeneity correlates with off-treatment sustained response to nucleos(t)ide analogues in e antigen-positive chronic hepatitis B patients

Finite treatment with nucleos(t)ide analogues (NAs) remains a great challenge for chronic hepatitis B in the clinic. This study aimed to investigate the relationship between intrahepatic quasispecies heterogeneity and the NAs off-treatment outcomes in a prospective cohort. Eighteen HBeAg-positive patients with chronic hepatitis B who achieved the cessation criteria underwent liver biopsy, and stopped treatment thereafter. Patients were followed up prospectively for 1 year. The reverse transcriptase (RT) gene of intrahepatic hepatitis B virus (HBV) was cloned and sequenced. Intrahepatic quasispecies heterogeneity and specific gene mutations were analysed using bioinformatic methods. Ten patients achieved sustained response, and eight patients developed viral relapse. The intrahepatic quasispecies Shannon entropy and nucleotide diversity within either RT or the surface (S) region of patients with sustained response were significantly higher (p < 0.05) than those of patients who had a viral relapse. Intrahepatic quasispecies Shannon entropy at the nucleotide level predicted the sustained off-treatment response (area under receiver operating characteristics curve 0.925; 95% CI 0.807-1.000; p 0.003). More positive selection sites and N-glycosylation mutations within the S region were found in patients with sustained response than in the patients with viral relapse (p < 0.01). Most of the positive selection sites in patients with sustained response were located in reported HLA-I-restricted or HLA-II-restricted epitopes. Intrahepatic quasispecies heterogeneity at the end of treatment was correlated with off-treatment outcomes in HBeAg-positive patients with chronic hepatitis B. More immune escape mutations were found within the S region in patients with sustained response. The higher intrahepatic quasispecies heterogeneity indicated a more robust immune control over HBV, which in turn maintained a sustained response after withdrawal of NAs.

Viral quasispecies

New generation sequencing is greatly expanding the capacity to examine the composition of mutant spectra of viral quasispecies in infected cells and host organisms. Here we review recent progress in the understanding of quasispecies dynamics, notably the occurrence of intra-mutant spectrum interactions, and implications of fitness landscapes for virus adaptation and de-adaptation. Complementation or interference can be established among components of the same mutant spectrum, dependent on the mutational status of the ensemble. Replicative fitness relates to an optimal mutant spectrum that provides the molecular basis for phenotypic flexibility, with implications for antiviral therapy. The biological impact of viral fitness renders particularly relevant the capacity of new generation sequencing to establish viral fitness landscapes. Progress with experimental model systems is becoming an important asset to understand virus behavior in the more complex environments faced during natural infections.

Modulation of flavivirus population diversity by RNA interference

To test the hypothesis that RNA interference (RNAi) imposes diversifying selection on RNA virus genomes, we quantified West Nile virus (WNV) quasispecies diversity after passage in Drosophila cells in which RNAi was left intact, depleted, or stimulated against WNV. As predicted, WNV diversity was significantly lower in RNAi-depleted cells and significantly greater in RNAi-stimulated cells relative to that in controls. These findings reveal that an innate immune defense can shape viral population structure.

Antiviral Strategies Based on Lethal Mutagenesis and Error Threshold

The concept of error threshold derived from quasispecies theory is at the basis of lethal mutagenesis, a new antiviral strategy based on the increase of virus mutation rate above an extinction threshold. Research on this strategy is justified by several inhibitor-escape routes that viruses utilize to ensure their survival. Successive steps in the transition from an organized viral quasispecies into loss of biologically meaningful genomic sequences are dissected. The possible connections between theoretical models and experimental observations on lethal mutagenesis are reviewed. The possibility of using combination of virus-specific mutagenic nucleotide analogues and broad-spectrum, non-mutagenic inhibitors is evaluated. We emphasize the power that quasispecies theory has had to stimulate exploration of new means to combat pathogenic viruses.

Theories of Lethal Mutagenesis: From Error Catastrophe to Lethal Defection

RNA viruses get extinct in a process called lethal mutagenesis when subjected to an increase in their mutation rate, for instance, by the action of mutagenic drugs. Several approaches have been proposed to understand this phenomenon. The extinction of RNA viruses by increased mutational pressure was inspired by the concept of the error threshold. The now classic quasispecies model predicts the existence of a limit to the mutation rate beyond which the genetic information of the wild type could not be efficiently transmitted to the next generation. This limit was called the error threshold, and for mutation rates larger than this threshold, the quasispecies was said to enter into error catastrophe. This transition has been assumed to foster the extinction of the whole population. Alternative explanations of lethal mutagenesis have been proposed recently. In the first place, a distinction is made between the error threshold and the extinction threshold, the mutation rate beyond which a population gets extinct. Extinction is explained from the effect the mutation rate has, throughout the mutational load, on the reproductive ability of the whole population. Secondly, lethal defection takes also into account the effect of interactions within mutant spectra, which have been shown to be determinant for the understanding the extinction of RNA virus due to an augmented mutational pressure. Nonetheless, some relevant issues concerning lethal mutagenesis are not completely understood yet, as so survival of the flattest, i.e. the development of resistance to lethal mutagenesis by evolving towards mutationally more robust regions of sequence space, or sublethal mutagenesis, i.e., the increase of the mutation rate below the extinction threshold which may boost the adaptability of RNA virus, increasing their ability to develop resistance to drugs (including mutagens). A better design of antiviral therapies will still require an improvement of our knowledge about lethal mutagenesis.

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.

Evaluation of anti-HIV-1 mutagenic nucleoside analogues

Because of their high mutation rates, RNA viruses and retroviruses replicate close to the threshold of viability. Their existence as quasi-species has pioneered the concept of "lethal mutagenesis" that prompted us to synthesize pyrimidine nucleoside analogues with antiviral activity in cell culture consistent with an accumulation of deleterious mutations in the HIV-1 genome. However, testing all potentially mutagenic compounds in cell-based assays is tedious and costly. Here, we describe two simple in vitro biophysical/biochemical assays that allow prediction of the mutagenic potential of deoxyribonucleoside analogues. The first assay compares the thermal stabilities of matched and mismatched base pairs in DNA duplexes containing or not the nucleoside analogues as follows. A promising candidate should display a small destabilization of the matched base pair compared with the natural nucleoside and the smallest gap possible between the stabilities of the matched and mismatched base pairs. From this assay, we predicted that two of our compounds, 5-hydroxymethyl-2'-deoxyuridine and 5-hydroxymethyl-2'-deoxycytidine, should be mutagenic. The second in vitro reverse transcription assay assesses DNA synthesis opposite nucleoside analogues inserted into a template strand and subsequent extension of the newly synthesized base pairs. Once again, only 5-hydroxymethyl-2'-deoxyuridine and 5-hydroxymethyl-2'-deoxycytidine are predicted to be efficient mutagens. The predictive potential of our fast and easy first line screens was confirmed by detailed analysis of the mutation spectrum induced by the compounds in cell culture because only compounds 5-hydroxymethyl-2'-deoxyuridine and 5-hydroxymethyl-2'-deoxycytidine were found to increase the mutation frequency by 3.1- and 3.4-fold, respectively.

Current perspectives on HIV-1 antiretroviral drug resistance

Current advancements in antiretroviral therapy (ART) have turned HIV-1 infection into a chronic and manageable disease. However, treatment is only effective until HIV-1 develops resistance against the administered drugs. The most recent antiretroviral drugs have become superior at delaying the evolution of acquired drug resistance. In this review, the viral fitness and its correlation to HIV-1 mutation rates and drug resistance are discussed while emphasizing the concept of lethal mutagenesis as an alternative therapy. The development of resistance to the different classes of approved drugs and the importance of monitoring antiretroviral drug resistance are also summarized briefly.

Hepatitis C virus: standard-of-care treatment

Hepatitis C virus (HCV) infection is curable by therapy. The antiviral treatment of chronic hepatitis C has been based for decades on the use of interferon (IFN)-α, combined with ribavirin. More recently, new therapeutic approaches that target essential components of the HCV life cycle have been developed, including direct-acting antiviral (DAA) and host-targeted agents (HTA). A new standard-of-care treatment has been approved in 2011 for patients infected with HCV genotype 1, based on a triple combination of pegylated IFN-α, ribavirin, and either telaprevir or boceprevir, two inhibitors of the HCV protease. New triple and quadruple combination therapies including pegylated IFN-α, ribavirin, and one or two DAAs/HTAs, respectively, are currently being evaluated in Phase II and III clinical trials. In addition, various options for all-oral, IFN-free regimens are currently being evaluated. This chapter describes the characteristics of the different drugs used in the treatment of chronic hepatitis C and those currently in development and provides an overview of the current and future standard-of-care treatments of chronic hepatitis C.

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.

Variability in mutational fitness effects prevents full lethal transitions in large quasispecies populations

The distribution of mutational fitness effects (DMFE) is crucial to the evolutionary fate of quasispecies. In this article we analyze the effect of the DMFE on the dynamics of a large quasispecies by means of a phenotypic version of the classic Eigen's model that incorporates beneficial, neutral, deleterious, and lethal mutations. By parameterizing the model with available experimental data on the DMFE of Vesicular stomatitis virus (VSV) and Tobacco etch virus (TEV), we found that increasing mutation does not totally push the entire viral quasispecies towards deleterious or lethal regions of the phenotypic sequence space. The probability of finding regions in the parameter space of the general model that results in a quasispecies only composed by lethal phenotypes is extremely small at equilibrium and in transient times. The implications of our findings can be extended to other scenarios, such as lethal mutagenesis or genomically unstable cancer, where increased mutagenesis has been suggested as a potential therapy.

Attenuation of human enterovirus 71 high-replication-fidelity variants in AG129 mice

In a screen for ribavirin resistance, a novel high-fidelity variant of human enterovirus 71 (EV71) with the single amino acid change L123F in its RNA-dependent RNA polymerase (RdRp or 3D) was identified. Based on the crystal structure of EV71 RdRp, L123 locates at the entrance of the RNA template binding channel, which might form a fidelity checkpoint. EV71 RdRp-L123F variants generated less progeny in a guanidine resistance assay and virus populations with lower mutation frequencies in cell culture passage due to their higher replication fidelity. However, compared with wild-type viruses, they did not show growth defects. In vivo infections further revealed that high-fidelity mutations L123F and G64R (previously reported) negatively impacted EV71 fitness and greatly reduced viral pathogenicity alone or together in AG129 mice. Interestingly, a variant with double mutations, RG/B4-G64R/L123F (where RG/B4 is an EV71 genotype B4 virus constructed by reverse genetics [RG])showed higher fidelity in vitro and less virulence in vivo than any one of the above two single mutants. The 50% lethal dose (LD50) of the double mutant increased more than 500 times compared with the LD50 of wild-type RG/B4 in mice. The results indicated that these high-fidelity variants exhibited an attenuated pathogenic phenotype in vivo and offer promise as a live attenuated EV71 vaccine.

IMPORTANCE

The error-prone nature of the RNA-dependent RNA polymerase (RdRp) of RNA viruses during replication results in quasispecies and aids survival of virus populations under a wide range of selective pressures. Virus variants with higher replication fidelity exhibit lower genetic diversity and attenuated pathogenicity in vivo. Here, we identified a novel high-fidelity mutation L123F in the RdRp of human enterovirus 71 (EV71). We further elucidated that EV71 variants with the RdRp-L123F mutation and/or the previously identified high-fidelity mutation RdRp-G64R were attenuated in an AG129 mouse model. As EV71 has emerged as a serious worldwide health threat, especially in developing countries in the Asia-Pacific region, we urgently need EV71 vaccines. Learning from the poliovirus vaccination, we prefer live attenuated EV71 vaccines to inactivated EV71 vaccines in order to effectively control EV71 outbreaks at low cost. Our results imply a new means of attenuating EV71 and reducing its mutation rate at the same time.

Characteristics of HIV-1 natural drug resistance-associated mutations in former paid blood donors in Henan Province, China

Background

Natural drug resistance is a major cause of antiviral treatment failure. The characteristics of HIV-1 natural drug resistance-associated mutations in former paid blood donors in Henan Province remain unclear.

Methods

One hundred and fifty HIV-1-positive plasma samples were collected. Plasma viral RNA was extracted for pol gene amplification and sequencing. The sequencing results were submitted to the HIV-1 drug resistance database for drug-resistance analysis.

Results

The rates of natural drug resistance and resistance-associated mutations were 17.7% (19/107) and 40.2% (43/107), respectively. The rates of PI major, PI minor, NRTI, and NNRTI mutations were: 0, 30.8% (33/107), 10.3% (11/107), and 18.7% (20/107), respectively. Nine cases (8.4%) had both NRTI and NNRTI resistance-associated mutations. Seven cases (6.5%) had PI minor, NRTI and NNRTI resistance-associated mutations. NNRTI resistance was the most serious, followed by NRTI resistance and PI resistance. Polymorphism mutation sites with mutation rates in the protease region higher than 60.0% were: L63A/P/S/T 89.7%, V77I 82.2%, I72E/M/K/T/V 80.4%, I93L 75.7%, and E35D 72.9%. Polymorphism mutation sites with mutation rates in the RT region higher than 60.0% were: I135A/L/M/R/T/V 93.5%, T200A/E/I/P/V 89.7%, Q278E/K/N/T 88.8%, S162C/Y 82.2%, and K277R/S 66.4%. The distribution of 107 gene sequences was scattered, with some drug-resistant strains grouped in the same cluster.

Conclusion

The natural drug resistance mutation rate of HIV-1 in former paid blood donors in Henan Province was 17.7%, with NNRTI resistance the most serious. The distribution of drug-resistant strains was scattered, with some correlations found in certain resistance loci.

Extinction of hepatitis C virus by ribavirin in hepatoma cells involves lethal mutagenesis

Lethal mutagenesis, or virus extinction produced by enhanced mutation rates, is under investigation as an antiviral strategy that aims at counteracting the adaptive capacity of viral quasispecies, and avoiding selection of antiviral-escape mutants. To explore lethal mutagenesis of hepatitis C virus (HCV), it is important to establish whether ribavirin, the purine nucleoside analogue used in anti-HCV therapy, acts as a mutagenic agent during virus replication in cell culture. Here we report the effect of ribavirin during serial passages of HCV in human hepatoma Huh-7.5 cells, regarding viral progeny production and complexity of mutant spectra. Ribavirin produced an increase of mutant spectrum complexity and of the transition types associated with ribavirin mutagenesis, resulting in HCV extinction. Ribavirin-mediated depletion of intracellular GTP was not the major contributory factor to mutagenesis since mycophenolic acid evoked a similar decrease in GTP without an increase in mutant spectrum complexity. The intracellular concentration of the other nucleoside-triphosphates was elevated as a result of ribavirin treatment. Mycophenolic acid extinguished HCV without an intervening mutagenic activity. Ribavirin-mediated, but not mycophenolic acid-mediated, extinction of HCV occurred via a decrease of specific infectivity, a feature typical of lethal mutagenesis. We discuss some possibilities to explain disparate results on ribavirin mutagenesis of HCV.

Response of hepatitis C virus to long-term passage in the presence of alpha interferon: multiple mutations and a common phenotype

Cell culture-produced hepatitis C virus (HCV) has been subjected to up to 100 serial passages in human hepatoma cells in the absence or presence of different doses of alpha interferon (IFN-α). Virus survival, genetic changes, fitness levels, and phenotypic traits have been examined. While high initial IFN-α doses (increasing from 1 to 4 IU/ml) did not allow HCV survival beyond passage 40, a gradual exposure (from 0.25 to 10 IU/ml) allowed the virus to survive for at least 100 passages. The virus passaged in the presence of IFN-α acquired IFN-α resistance as evidenced by enhanced progeny production and viral protein expression in an IFN-α environment. A partial IFN-α resistance was also noted in populations passaged in the absence of IFN-α. All lineages acquired adaptative mutations, and multiple, nonsynonymous mutations scattered throughout the genome were present in IFN-α-selected populations. Comparison of consensus sequences indicates a dominance of synonymous versus nonsynonymous substitutions. IFN-α-resistant populations displayed decreased sensitivity to a combination of IFN-α and ribavirin. A phenotypic trait common to all assayed viral populations is the ability to increase shutoff host cell protein synthesis, accentuated in infections with IFN-α-selected populations carried out in the presence of IFN-α. The trait was associated with enhanced phosphorylation of protein kinase R (PKR) and eIF2α, although other contributing factors are likely. The results suggest that multiple, independent mutational pathways can confer IFN-α resistance to HCV and might explain why no unified picture has been obtained regarding IFN-α resistance in vivo.

Back to the future: revisiting HIV-1 lethal mutagenesis

The concept of eliminating HIV-1 infectivity by elevating the viral mutation rate was first proposed over a decade ago, even though the general concept had been conceived earlier for RNA viruses. Lethal mutagenesis was originally viewed as a novel chemotherapeutic approach for treating HIV-1 infection in which use of a viral mutagen would over multiple rounds of replication lead to the lethal accumulation of mutations, rendering the virus population noninfectious - known as the slow mutation accumulation model. There have been limitations in obtaining good efficacy data with drug leads, leaving some doubt on clinical translation. More recent studies of the apolipoprotein B mRNA editing complex 3 (APOBEC3) proteins as well as new progress in the use of nucleoside analogs for inducing lethal mutagenesis have helped to refocus attention on rapid induction of HIV-1 lethal mutagenesis in a single or limited number of replication cycles leading to a rapid mutation accumulation model.

Running loose or getting lost: how HIV-1 counters and capitalizes on APOBEC3-induced mutagenesis through its Vif protein

Human immunodeficiency virus-1 (HIV-1) dynamics reflect an intricate balance within the viruses’ host. The virus relies on host replication factors, but must escape or counter its host’s antiviral restriction factors. The interaction between the HIV-1 protein Vif and many cellular restriction factors from the APOBEC3 protein family is a prominent example of this evolutionary arms race. The viral infectivity factor (Vif) protein largely neutralizes APOBEC3 proteins, which can induce in vivo hypermutations in HIV-1 to the extent of lethal mutagenesis, and ensures the production of viable virus particles. HIV-1 also uses the APOBEC3-Vif interaction to modulate its own mutation rate in harsh or variable environments, and it is a model of adaptation in a coevolutionary setting. Both experimental evidence and the substantiation of the underlying dynamics through coevolutionary models are presented as complementary views of a coevolutionary arms race.

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.