Mutant spectra in virus behavior

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.

Assessing in vivo mutation frequencies and creating a high-resolution genome-wide map of fitness costs of Hepatitis C virus

Like many viruses, Hepatitis C Virus (HCV) has a high mutation rate, which helps the virus adapt quickly, but mutations come with fitness costs. Fitness costs can be studied by different approaches, such as experimental or frequency-based approaches. The frequency-based approach is particularly useful to estimate in vivo fitness costs, but this approach works best with deep sequencing data from many hosts are. In this study, we applied the frequency-based approach to a large dataset of 195 patients and estimated the fitness costs of mutations at 7957 sites along the HCV genome. We used beta regression and random forest models to better understand how different factors influenced fitness costs. Our results revealed that costs of nonsynonymous mutations were three times higher than those of synonymous mutations, and mutations at nucleotides A or T had higher costs than those at C or G. Genome location had a modest effect, with lower costs for mutations in HVR1 and higher costs for mutations in Core and NS5B. Resistance mutations were, on average, costlier than other mutations. Our results show that in vivo fitness costs of mutations can be site and virus specific, reinforcing the utility of constructing in vivo fitness cost maps of viral genomes.

Novel Synthesis and Phenotypic Analysis of Mutant Clouds for Hepatitis E Virus Genotype 1

Many RNA viruses exist as an ensemble of genetically diverse, replicating populations known as a mutant cloud. The genetic diversity (cloud size) and composition of this mutant cloud may influence several important phenotypic features of the virus, including its replication capacity. We applied a straightforward, bacterium-free approach using error-prone PCR coupled with reverse genetics to generate infectious mutant RNA clouds with various levels of genetic diversity from a genotype 1 strain of hepatitis E virus (HEV). Cloning and sequencing of a genomic fragment encompassing 70% of open reading frame 1 (ORF1) or of the full genome from variants in the resultant clouds showed the occurrence of nucleotide mutations at a frequency on the order of 10^-3 per nucleotide copied and the existence of marked genetic diversity, with a high normalized Shannon entropy value. The mutant clouds showed transient replication in cell culture, while wild-type HEV did not. Cross-sectional data from these cell cultures supported the existence of differential effects of clouds of various sizes and compositions on phenotypic characteristics, such as the replication level of (+)-RNA progeny, the amounts of double-stranded RNA (a surrogate for the rate of viral replication) and ORF1 protein, and the expression of interferon-stimulated genes. Since mutant cloud size and composition influenced the viral phenotypic properties, a better understanding of this relationship may help to provide further insights into virus evolution and prediction of emerging viral diseases.IMPORTANCE Several biological or practical limitations currently prevent the study of phenotypic behavior of a mutant cloud in vitro We developed a simple and rapid method for synthesizing mutant clouds of hepatitis E virus (HEV), a single-stranded (+)-RNA [ss(+) RNA] virus, with various and controllable levels of genetic diversity, which could then be used in a cell culture system to study the effects of cloud size and composition on viral phenotype. In a cross-sectional analysis, we demonstrated that a particular mutant cloud which had an extremely high genetic diversity had a replication rate exceeding that of wild-type HEV. This method should thus provide a useful model for understanding the phenotypic behavior of ss(+) RNA viruses.

Involvement of a joker mutation in a polymerase-independent lethal mutagenesis escape mechanism

We previously characterized a foot-and-mouth disease virus (FMDV) with three amino acid replacements in its polymerase (3D) that conferred resistance to the mutagenic nucleoside analogue ribavirin. Here we show that passage of this mutant in the presence of high ribavirin concentrations resulted in selection of viruses with the additional replacement I248T in 2C. This 2C substitution alone (even in the absence of replacements in 3D) increased FMDV fitness mainly in the presence of ribavirin, prevented an incorporation bias in favor of A and U associated with ribavirin mutagenesis, and conferred the ATPase activity of 2C decreased sensitivity to ribavirin-triphosphate. Since in previous studies we described that 2C with I248T was selected under different selective pressures, this replacement qualifies as a joker substitution in FMDV evolution. The results have identified a role of 2C in nucleotide incorporation, and have unveiled a new polymerase-independent mechanism of virus escape to lethal mutagenesis.

•

A replacement in FMDV protein 2C confers reduced sensitivity to the mutagen ribavirin.

•

The effect of the replacement is to prevent a mutational bias evoked by ribavirin.

•

2C has an effect in nucleotide incorporation by the FMDV polymerase.

•

We describe a new molecular mechanism of escape to ribavirin-mediated extinction.

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.

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.

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.

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.

Mutagenesis-mediated decrease of pathogenicity as a feature of the mutant spectrum of a viral population

BACKGROUND

RNA virus populations are heterogeneous ensembles of closely related genomes termed quasispecies. This highly complex distribution of variants confers important properties to RNA viruses and influences their pathogenic behavior. It has been hypothesized that increased mutagenesis of viral populations, by treatment with mutagenic agents, can induce alterations in the pathogenic potential of a virus population. In this work we investigate whether mutagenized foot-and-mouth disease virus (FMDV) populations display changes in their virulence in mice.

METHODOLOGY AND PRINCIPAL FINDINGS

FMDV C-S8c1 was passaged in BHK cells in the presence of the mutagenic agent ribavirin. Decline in viral titer and viral RNA progeny was observed in the first passage, fluctuating around a constant value thereafter. Hence, the specific infectivity remained stable during the passages. The viral population harvested from passage 9 (P9 R) showed decreased virulence in mice, with a lethal dose 50 (LD(50)) >10(4) PFU, as compared with LD(50) of 50 PFU of the parental population FMDV C-S8c1. This decrease in virulence was associated to a 20-fold increase in the mutation frequency of the P9 R population with respect to C-S8c1. Interestingly, individual biological clones isolated from the attenuated population P9 R were as virulent as the parental virus C-S8c1. Furthermore, a mixed population of C-S8c1 and P9 R was inoculated into mice and showed decreased virulence as compared to C-S8c1, suggesting that population P9 R is able to suppress the virulent phenotype of C-S8c1.

CONCLUSION

Ribavirin-mediated mutagenesis of an FMDV population resulted in attenuation in vivo, albeit a large proportion of its biological clones displayed a highly virulent phenotype. These results, together with the suppression of C-S8c1 by mutagenized P9 R population, document a suppressive effect of mutagenized viral quasispecies in vivo, and suggest novel approaches to the treatment and prevention of viral diseases.

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.

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.

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.

Lethal mutagenesis of viruses

Lethal mutagenesis aims at extinguishing viruses by increased mutagenesis prompted by virus-specific mutagenic agents, mainly nucleoside analogues. It is derived from the error threshold relationship of quasispecies theory, and it is slowly finding its way towards a clinical application. We summarize the current situation of research in this field of antiviral therapy.

Quasispecies as a matter of fact: viruses and beyond

We review the origins of the quasispecies concept and its relevance for RNA virus evolution, viral pathogenesis and antiviral treatment strategies. We emphasize a critical point of quasispecies that refers to genome collectivities as the unit of selection, and establish parallels between RNA viruses and some cellular systems such as bacteria and tumor cells. We refer also to tantalizing new observations that suggest quasispecies behavior in prions, perhaps as a result of the same quantum-mechanical indeterminations that underlie protein conformation and error-prone replication in genetic systems. If substantiated, these observations with prions could lead to new research on the structure-function relationship of non-nucleic acid biological molecules.

Lethal mutagenesis of foot-and-mouth disease virus involves shifts in sequence space

Lethal mutagenesis or virus transition into error catastrophe is an antiviral strategy that aims at extinguishing a virus by increasing the viral mutation rates during replication. The molecular basis of lethal mutagenesis is largely unknown. Previous studies showed that a critical substitution in the foot-and-mouth disease virus (FMDV) polymerase was sufficient to allow the virus to escape extinction through modulation of the transition types induced by the purine nucleoside analogue ribavirin. This substitution was not detected in mutant spectra of FMDV populations that had not replicated in the presence of ribavirin, using standard molecular cloning and nucleotide sequencing. Here we selectively amplify and analyze low-melting-temperature cDNA duplexes copied from FMDV genome populations passaged in the absence or presence of ribovirin Hypermutated genomes with high frequencies of A and U were present in both ribavirin -treated and untreated populations, but the major effect of ribavirin mutagenesis was to accelerate the occurrence of AU-rich mutant clouds during the early replication rounds of the virus. The standard FMDV quasispecies passaged in the absence of ribavirin included the salient transition-modulating, ribavirin resistance mutation, whose frequency increased in populations treated with ribavirin. Thus, even nonmutagenized FMDV quasispecies include a deep, mutationally biased portion of sequence space, in support of the view that the virus replicates close to the error threshold for maintenance of genetic information.

Influence of mutagenesis and viral load on the sustained low-level replication of an RNA virus

Lethal mutagenesis is an antiviral strategy that aims to extinguish viruses as a consequence of enhanced mutation rates during virus replication. The molecular mechanisms that underlie virus extinction by mutagenic nucleoside analogues are not well understood. When mutagenic agents and antiviral inhibitors are administered sequentially or in combination, interconnected and often conflicting selective constraints can influence the fate of the virus either towards survival through selection of mutagen-escape or inhibitor-escape mutants or towards extinction. Here we report a study involving the mutagenesis of foot-and-mouth disease virus (FMDV) by the nucleoside analogue ribavirin (R) and the effect of R-mediated mutagenesis on the selection of FMDV mutants resistant to the inhibitor of RNA replication, guanidine hydrochloride (GU). The results show that under comparable (and low) viral load, an inhibitory activity by GU could not substitute for an equivalent inhibitory activity by R in driving FMDV to extinction. Both the prior history of R mutagenesis and the viral population size influenced the selection of GU-escape mutants. A sufficiently low viral load allowed continued viral replication without selection of inhibitor-escape mutants, irrespective of the history of mutagenesis. These observations imply that reductions of viral load as a result of a mutagenic treatment may provide an opportunity either for immune-mediated clearing of a virus or for an alternative antiviral intervention, even if extinction is not initially achieved.