Search

We report that ribavirin exerts an inhibitory and mutagenic activity on SARS-CoV-2-infecting Vero cells, with a therapeutic index higher than 10. Deep sequencing analysis of the mutant spectrum of SARS-CoV-2 replicating in the absence or presence of ribavirin indicated an increase in the number of mutations, but not in deletions, and modification of diversity indices, expected from a mutagenic activity. Notably, the major mutation types enhanced by replication in the presence of ribavirin were A→G and U→C transitions, a pattern which is opposite to the dominance of G→A and C→U transitions previously described for most RNA viruses. Implications of the inhibitory activity of ribavirin, and the atypical mutational bias produced on SARS-CoV-2, for the search for synergistic anti-COVID-19 lethal mutagen combinations are discussed.

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.

Hepatitis C virus (HCV) core is a highly conserved and multifunctional protein that forms the viral capsid, making it an attractive target for HCV detection and inhibition. Aptamers are in vitro selected, single-stranded nucleic acids (RNA or ssDNA) with growing applicability in viral diagnostics and therapy. We have carried out DNA and RNA in vitro selection against six different variants of HCV core protein: two versions of the full-length protein of genotype 1, and the hydrophilic domain of genotypes 1 to 4. The aptamer populations obtained were analyzed by means of Ultra-Deep Sequencing (UDS), the most abundant sequences were identified and a number of highly represented sequence motifs were unveiled. Affinity (measured as the dissociation constant, Kd) of the most abundant DNA and RNA aptamers were quantified using Enzyme-Linked OligoNucleotide Assay (ELONA)-based methods. Some aptamers with nanomolar or subnanomolar Kd values (as low as 0.4 nM) were the common outcome of DNA and RNA selections against different HCV core variants. They were tested in sandwich and competitive biosensor assays, reaching a limit of detection for HCV core of 2 pM. Additionally, the two most prevalent and high affinity aptamers were assayed in Huh-7.5 reporter cell lines infected with HCV, where they decreased both the viral progeny titer and the extracellular viral RNA level, while increasing the amount of intracellular viral RNA. Our results suggest that these aptamers inhibit HCV capsid assembly and virion formation, thus making them good candidate molecules for the design of novel therapeutic approaches for hepatitis C., The work at CAB was funded by the Spanish Ministerio de Economía y Competitividad (MINECO) and EU FEDER program grants no. RTC-2014-1986-1 and BIO2016-79618-R, the Spanish Ministerio de Ciencia e Innovación (MICINN) grant no. PID2019-104903RB-I00, and the Spanish Agencia Estatal de Investigación (AEI) Project no. MDM-2017-0737 - Unidad de Excelencia “María de Maeztu”. The work at CAB also benefits from the interdisciplinary framework provided by CSIC through “LifeHUB.CSIC” initiative (PIE 202120E047-Conexiones-Life). B.T-V. is supported by the predoctoral fellowship TS17/16 from INTA and by the CSIC “Garantía Juvenil” contract CAM19_PRE_CAB_001 funded by Comunidad Autónoma de Madrid (CAM). The work at CBMSO and IIS-FJD was supported by grants SAF2017-87846-R and BFU2017-91384-EXP from Ministerio de Ciencia, Innovación y Universidades (MCIU), PI18/00210 from Instituto de Salud Carlos III (ISCIII), and S2018/BAA-4370 (PLATESA2) from CAM/FEDER. C.G.-C. is supported by the predoctoral contract PRE2018-083422 from MCIU. C.P. is supported by the Miguel Servet program of ISCIII (CPII19/00001), cofinanced by the European Regional Development Fund (ERDF). CIBERehd is funded by ISCIII. Institutional grants from the Fundación Ramón Areces and Banco Santander to the CBMSO are also acknowledged

In the course of experiments aimed at deciphering the inhibition mechanism of mycophenolic acid and ribavirin in hepatitis C virus (HCV) infection, we observed an inhibitory effect of the nucleoside guanosine (Gua). Here, we report that Gua, and not the other standard nucleosides, inhibits HCV replication in human hepatoma cells. Gua did not directly inhibit the in vitro polymerase activity of NS5B, but it modified the intracellular levels of nucleoside di- and tri-phosphates (NDPs and NTPs), leading to deficient HCV RNA replication and reduction of infectious progeny virus production. Changes in the concentrations of NTPs or NDPs modified NS5B RNA polymerase activity in vitro, in particular de novo RNA synthesis and template switching. Furthermore, the Gua-mediated changes were associated with a significant increase in the number of indels in viral RNA, which may account for the reduction of the specific infectivity of the viral progeny, suggesting the presence of defective genomes. Thus, a proper NTP:NDP balance appears to be critical to ensure HCV polymerase fidelity and minimal production of defective genomes., CP is supported by the Miguel Servet program (grants CP14/00121 and CP19/00001) of the Instituto de Salud Carlos III cofinanced by FEDER. CP has received funding from Ministerio de Ciencia, Innovacio´n y Universidades (grant BFU2017-91384-EXP), from Instituto de Salud Carlos III (grants PI18/00210 and PI21/00139), from Fundación La Marato´ (grant 525/C/2021), and from CSIC (grant CSIC-COV19-014). ED has received funding from CIBERehd (Centro de Investigación en Red de Enfermedades Hepáticas y Digestivas, Instituto de Salud Carlos III), from Ministerio de Economía y Competitividad (grants SAF2014-52400-R, SAF2017-87846-R, and PID2020-113888RB-I00), and from Comunidad de Madrid/FEDER (grants S2013/ABI-2906 PLATESA, and S2018/BAA-4370 PLATESA2). AM has received funding from Ministerio de Economía y Competitividad (grants SAF2016-80451-P, PID2019-106068GB-I00, EQC2018-004420-P, and EQC2018-004631-P), and Plan Propio of Universidad de Castilla-La Mancha. A.G-P and L.DM have received funding from Consejería de Transformación Económica, Industria, Conocimiento y Universidades, Junta de Andalucía, cofinanced by FEDER and FSE (grants BIO-264, and P10-CVI-6561) and Plan Propio of Universidad de Málaga.

RNA viruses replicate as complex mutant spectra termed viral quasispecies. The frequency of each individual genome in a mutant spectrum depends on its rate of generation and its relative fitness in the replicating population ensemble. The advent of deep sequencing methodologies allows for the first-time quantification of haplotype abundances within mutant spectra. There is no information on the haplotype profile of the resident genomes and how the landscape evolves when a virus replicates in a controlled cell culture environment. Here, we report the construction of intramutant spectrum haplotype landscapes of three amplicons of the NS5A-NS5B coding region of hepatitis C virus (HCV). Two-dimensional (2D) neural networks were constructed for 44 related HCV populations derived from a common clonal ancestor that was passaged up to 210 times in human hepatoma Huh-7.5 cells in the absence of external selective pressures. The haplotype profiles consisted of an extended dense basal platform, from which a lower number of protruding higher peaks emerged. As HCV increased its adaptation to the cells, the number of haplotype peaks within each mutant spectrum expanded, and their distribution shifted in the 2D network. The results show that extensive HCV replication in a monotonous cell culture environment does not limit HCV exploration of sequence space through haplotype peak movements. The landscapes reflect dynamic variation in the intramutant spectrum haplotype profile and may serve as a reference to interpret the modifications produced by external selective pressures or to compare with the landscapes of mutant spectra in complex in vivo environments. IMPORTANCE The study provides for the first time the haplotype profile and its variation in the course of virus adaptation to a cell culture environment in the absence of external selective constraints. The deep sequencing-based self-organized maps document a two-layer haplotype distribution with an ample basal platform and a lower number of protruding peaks. The results suggest an inferred intramutant spectrum fitness landscape structure that offers potential benefits for virus resilience to mutational inputs.

Viral quasispecies are dynamic distributions of nonidentical but closely related mutant and recombinant viral genomes subjected to a continuous process of genetic variation, competition, and selection that may act as a unit of selection. The quasispecies concept owes its theoretical origins to a model for the origin of life as a collection of mutant RNA replicators. Independently, experimental evidence for the quasispecies concept was obtained from sampling of bacteriophage clones, which revealed that the viral populations consisted of many mutant genomes whose frequency varied with time of replication. Similar findings were made in animal and plant RNA viruses. Quasispecies became a theoretical framework to understand viral population dynamics and adaptability. The evidence came at a time when mutations were considered rare events in genetics, a perception that was to change dramatically in subsequent decades. Indeed, viral quasispecies was the conceptual forefront of a remarkable degree of biological diversity, now evident for cell populations and organisms, not only for viruses. Quasispecies dynamics unveiled complexities in the behavior of viral populations,with consequences for disease mechanisms and control strategies. This review addresses the origin of the quasispecies concept, its major implications on both viral evolution and antiviral strategies, and current and future prospects., SAF2014-52400-R from Ministerio de Economía y Competitividad (MINECO), SAF2017-87846-R and BFU2017-91384-EXP from Ministerio de Ciencia, Innovación y Universidades (MCIU), PI18/00210 from Instituto de Salud Carlos III, S2013/ABI-2906 (PLATESA from Comunidad de Madrid/FEDER), and S2018/BAA-4370 (PLATESA2 from Comunidad de Madrid/FEDER). C.P. is supported by the Miguel Servet program of the Instituto de Salud Carlos III (CP14/00121 and CPII19/00001), cofinanced by the European Regional Development Fund (ERDF). CIBERehd (Centro de Investigación en Red de Enfermedades Hepáticas y Digestivas) is funded by Instituto de Salud Carlos III. Institutional grants from the Fundación Ramón Areces and Banco Santander

The error rate displayed during template copying to produce viral RNA progeny is a biologically relevant parameter of the replication complexes of viruses. It has consequences for virus–host interactions, and it represents the first step in the diversification of viruses in nature. Measurements during infections and with purified viral polymerases indicate that mutation rates for RNA viruses are in the range of 10 to 10 copying errors per nucleotide incorporated into the nascent RNA product. Although viruses are thought to exploit high error rates for adaptation to changing environments, some of them possess misincorporation correcting activities. One of them is a proofreading-repair 3 to 5 exonuclease present in coronaviruses that may decrease the error rate during replication. Here we review experimental evidence and models of information maintenance that explain why elevated mutation rates have been preserved during the evolution of RNA (and some DNA) viruses. The models also offer an interpretation of why error correction mechanisms have evolved to maintain the stability of genetic information carried out by large viral RNA genomes such as the coronaviruses., Ministerio de Economía y Competitividad (MINECO), grant number SAF2014-52400-R, SAF2017-87846-R and BFU2017-91384-EXP from Ministerio de Ciencia, Innovación y Universidades (MCIU), PI18/00210 from Instituto de Salud Carlos III, S2013/ABI-2906 (PLATESA from Comunidad de Madrid/FEDER), and S2018/BAA-4370 (PLATESA2 from Comunidad de Madrid/FEDER). C.P. is supported by the Miguel Servet program of the Instituto de Salud Carlos III (CPII19/00001), cofinanced by the European Regional Development Fund (ERDF). CIBERehd (Centro de Investigación en Red de Enfermedades Hepáticas y Digestivas) is funded by Instituto de Salud Carlos III. Institutional grants from the Fundación Ramón Areces and Banco Santander to the CBMSO are also acknowledged

Fitness landscapes reflect the adaptive potential of viruses. There is no information on how fitness peaks evolve when a virus replicates extensively in a controlled cell culture environment. Here we report the construction of Self-Organized Maps (SOMs), based on deep sequencing reads of three amplicons of the NS5A-NS5B-coding region of hepatitis C virus (HCV). A two-dimensional neural network was constructed and organized according to sequence relatedness. The third dimension of the fitness profile was given by the haplotype frequencies at each neuron. Fitness maps were derived for 44 HCV populations that share a common ancestor that was passaged up to 210 times in human hepatoma Huh-7.5 cells. As the virus increased its adaptation to the cells, the number of fitness peaks expanded, and their distribution shifted in sequence space. The landscape consisted of an extended basal platform, and a lower number of protruding higher fitness peaks. The function that relates fitness level and peak abundance corresponds a power law, a relationship observed with other complex natural phenomena. The dense basal platform may serve as spring-board to attain high fitness peaks. The study documents a highly dynamic, double-layer fitness landscape of HCV when evolving in a monotonous cell culture environment. This information may help interpreting HCV fitness landscapes in complex in vivo environments.IMPORTANCEThe study provides for the first time the fitness landscape of a virus in the course of its adaptation to a cell culture environment, in absence of external selective constraints. The deep sequencing-based self-organized maps document a two-layer fitness distribution with an ample basal platform, and a lower number of protruding, high fitness peaks. This landscape structure offers potential benefits for virus resilience to mutational inputs.

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., The work at CBMSO was supported by grants SAF2014-52400-R from Ministerio de Economía y Competitividad (MINECO), SAF2017-87846-R and BFU2017-91384-EXP from Ministerio de Ciencia, Innovación y Universidades (MCIU), PI18/00210 from Instituto de Salud Carlos III, S2013/ABI-2906 (PLATESA from Comunidad de Madrid/FEDER), and S2018/BAA-4370 (PLATESA2 from Comunidad de Madrid/FEDER). C.P. is supported by the Miguel Servet program of the Instituto de Salud Carlos III (CPII19/00001), cofinanced by the European Regional Development Fund (ERDF). CIBERehd (Centro de Investigación en Red de Enfermedades Hepáticas y Digestivas) is funded by Instituto de Salud Carlos III. Institutional grants from the Fundación Ramón Areces and Banco Santander to the CBMSO are also acknowledged. The team at CBMSO belongs to the Global Virus Network (GVN). The work in Barcelona was supported by Instituto de Salud Carlos III, cofinanced by the European Regional Development Fund (ERDF) Grant No. PI19/00301 and by the Centro para el Desarrollo Tecnológico Industrial (CDTI) from the MICIU, Grant No. IDI-20200297. Work at CAB was supported by MINECO grant BIO2016-79618R and PID2019-104903RB-I00 (funded by the EU under the FEDER program) and by the Spanish State research agency (AEI) through project number MDM-2017-0737 Unidad de Excelencia “María de Maeztu”-Centro de Astrobiología (CSIC-INTA). C.G.-C. is supported by predoctoral contract PRE2018-083422 from MCIU. B.M.-G. is supported by predoctoral contract PFIS FI19/00119 from Instituto de Salud Carlos III (Ministerio de Sanidad y Consumo), cofinanced by Fondo Social Europeo (FSE).

Despite the high virological response rates achieved with current directly-acting antiviral agents (DAAs) against hepatitis C virus (HCV), around 2% to 5% of treated patients do not achieve a sustained viral response. Identification of amino acid substitutions associated with treatment failure requires analytical designs, such as subtype-specific ultra-deep sequencing (UDS) methods for HCV characterization and patient management. Using this procedure, we have identified six highly represented amino acid substitutions (HRSs) in NS5A and NS5B of HCV from 220 patients who failed therapy, which are not bona fide resistance-associated substitutions (RAS). They were present frequently in basal and post-treatment virus of patients who failed therapy to different DAA-based therapies. Contrary to several RAS, HRSs belong to the acceptable subset of substitutions according to the PAM250 replacement matrix. Coherently, their mutant frequency, measured by the number of deep sequencing reads within the HCV quasispecies that encode the relevant substitutions, ranged between 90% and 100% in most cases. Also, they have limited predicted disruptive effects on the three-dimensional structures of the proteins harboring them. Possible mechanisms of HRS origin and dominance, as well as their potential predictive value of treatment response are discussed.

Despite the high virological response rates achieved with current directly acting antiviral agents (DAAs) against hepatitis C virus (HCV), around 2% to 5% of treated patients do not achieve a sustained viral response. The identification of amino acid substitutions associated with treatment failure requires analytical designs, such as subtype-specific ultradeep sequencing (UDS) methods, for HCV characterization and patient management. Using this procedure, we have identified six highly represented amino acid substitutions (HRSs) in NS5A and NS5B of HCV, which are not bona fide resistance-associated substitutions (RAS), from 220 patients who failed therapy. They were present frequently in basal and posttreatment virus of patients who failed different DAA-based therapies. Contrary to several RAS, HRSs belong to the acceptable subset of substitutions according to the PAM250 replacement matrix. Their mutant frequency, measured by the number of deep sequencing reads within the HCV quasispecies that encode the relevant substitutions, ranged between 90% and 100% in most cases. They also have limited predicted disruptive effects on the three-dimensional structures of the proteins harboring them. Possible mechanisms of HRS origin and dominance, as well as their potential predictive value for treatment response, are discussed., The work at CBMSO was supported by grants SAF2014-52400-R from Ministerio de Economía y Competitividad (MINECO), SAF2017-87846-R and BFU2017-91384-EXP from Ministerio de Ciencia, Innovación y Universidades (MICIU), PI18/00210 from Instituto de Salud Carlos III, S2013/ABI-2906 (PLATESA from Comunidad de Madrid/FEDER), and S2018/BAA-4370 (PLATESA2 from Comunidad de Madrid/FEDER). C.P. is supported by the Miguel Servet program of the Instituto de Salud Carlos III (CP14/00121 and CPII19/00001), cofinanced by the European Regional Development Fund (ERDF). CIBERehd (Centro de Investigación en Red de Enfermedades Hepáticas y Digestivas) is funded by Instituto de Salud Carlos III. Institutional grants from the Fundación Ramón Areces and Banco Santander to the CBMSO are also acknowledged. The team at CBMSO belongs to the Global Virus Network (GVN). The work in Barcelona was supported by Instituto de Salud Carlos III, cofinanced by the European Regional Development Fund (ERDF) grant number PI19/00301 and by the Centro para el Desarrollo Tecnológico Industrial (CDTI) from the MICIU, grant number IDI-20151125. Work at CAB was supported by MINECO grant BIO2016-79618R and PID2019-104903RB-I00 (funded by the EU under the FEDER program) and by the Spanish State research agency (AEI) through project number MDM-2017-0737 Unidad de Excelencia “María de Maeztu”-Centro de Astrobiología (CSIC-INTA). Work at IBMB was supported by MICIN grant BIO2017-83906-P (funded by the EU under the FEDER program). C.G.-C. is supported by predoctoral contract PRE2018-083422 from MICIU. B.M.-G. is supported by predoctoral contract PFIS FI19/00119 from Instituto de Salud Carlos III (Ministerio de Sanidad y Consumo), cofinanced by Fondo Social Europeo (FSE).

Previous studies documented that long-term hepatitis C virus (HCV) replication in human hepatoma Huh-7.5 cells resulted in viral fitness gain, expansion of the mutant spectrum, and several phenotypic alterations. In the present work, we show that mutational waves (changes in frequency of individual mutations) occurred continuously and became more prominent as the virus gained fitness. They were accompanied by an increasing proportion of heterogeneous genomic sites that affected 1 position in the initial HCV population and 19 and 69 positions at passages 100 and 200, respectively. Analysis of biological clones of HCV showed that these dynamic events affected infectious genomes, since part of the fluctuating mutations became incorporated into viable genomes. While 17 mutations were scored in 3 biological clones isolated from the initial population, the number reached 72 in 3 biological clones from the population at passage 200. Biological clones differed in their responses to antiviral inhibitors, indicating a phenotypic impact of viral dynamics. Thus, HCV adaptation to a specific constant environment (cell culture without external influences) broadens the mutant repertoire and does not focus the population toward a limited number of dominant genomes. A retrospective examination of mutant spectra of foot-and-mouth disease virus passaged in cell cultures suggests a parallel behavior here described for HCV. We propose that virus diversification in a constant environment has its basis in the availability of multiple alternative mutational pathways for fitness gain. This mechanism of broad diversification should also apply to other replicative systems characterized by high mutation rates and large population sizes. IMPORTANCE The study shows that extensive replication of an RNA virus in a constant biological environment does not limit exploration of sequence space and adaptive options. There was no convergence toward a restricted set of adapted genomes. Mutational waves and mutant spectrum broadening affected infectious genomes. Therefore, profound modifications of mutant spectrum composition and consensus sequence diversification are not exclusively dependent on environmental alterations or the intervention of population bottlenecks., The work at CBMSO was supported by grants SAF2014-52400-R from Ministerio de Economía y Competitividad (MINECO); SAF2017-87846-R and BFU2017-91384-EXP from Ministerio de Ciencia, Innovación y Universidades (MICIU); and S2013/ABI-2906 (PLATESA from Comunidad de Madrid [CAM]/FEDER), S2018/BAA-4370 (PLATESA2 from CAM/FEDER), and PI18/00210 from Instituto de Salud Carlos III. C.P. is supported by the Miguel Servet program of the Instituto de Salud Carlos III (CP14/00121), cofinanced by the European Regional Development Fund (ERDF). CIBERehd is funded by Instituto de Salud Carlos III. Institutional grants from the Fundación Ramón Areces and Banco Santander to the CBMSO are also acknowledged. The team at CBMSO belongs to the Global Virus Network (GVN). The work in Barcelona was supported by Instituto de Salud Carlos III, cofinanced by the European Regional Development Fund (ERDF) (grant numbers PI16/00337) and by the Centro para el Desarrollo Tecnológico Industrial (CDTI) from the Spanish Ministry of Economy and Business (grant number IDI-20151125). Work at CAB was supported by MINECO grant BIO2016-79618-R (funded by the European Union under the FEDER program) and by the Spanish State Research Agency (AEI) through project number MDM-2017-0737 Unidad de Excelencia “María de Maeztu”-Centro de Astrobiologia (CSIC-INTA).

Lethal mutagenesis is an antiviral approach that consists of extinguishing a virus by an excess of mutations acquired during replication in the presence of a mutagenic agent, often a nucleotide analogue. One of its advantages is its broad-spectrum nature, which renders the strategy potentially effective against emergent RNA viral infections. Here we describe the synergistic lethal mutagenesis of hepatitis C virus (HCV) by a combination of favipiravir (T-705) and ribavirin. Synergy has been documented over a broad range of analogue concentrations using the Chou-Talalay method implemented in CompuSyn graphics software, with the average dose reduction index (DRI) being above 1 (68.02 101.6 for favipiravir and 5.83 6.07 for ribavirin) and the average combination indices (CI) being below 1 (0.52 0.28). Furthermore, analogue concentrations that individually did not extinguish high-fitness HCV in 10 serial infections extinguished high-fitness HCV in 1 to 2 passages when used in combination. Although both analogues displayed a preference for G ¡ A and C ¡ U transitions, deep sequencing analysis of mutant spectra indicated a different preference of the two analogues for the mutation sites, thus unveiling a new possible synergy mechanism in lethal mutagenesis. The prospects for synergy among mutagenic nucleotides as a strategy to confront emerging viral infections are discussed, C.P. is supported by the Miguel Servet program of the Instituto de Salud Carlos III (grant CP14/00121), cofinanced by the European Regional Development Fund (ERDF). The Centro de Investigación en Red de Enfermedades Hepáticas y Digestivas (CIBERehd) is funded by the Instituto de Salud Carlos III. The work in Madrid, Spain, was supported by grant SAF2014-52400-R from the Ministerio de Economía y Competitividad, grants SAF2017-87846-R, BFU2017-91384-EXP, and PI18/00210 from the Ministerio de Ciencia, Innovación y Universidades, grant S2013/ABI-2906 from PLATESA from the Comunidad de Madrid/FEDER, and grant P2018/BAA-4370 from PLATESA2 from the Comunidad de Madrid/FEDER. The work in Barcelona, Spain, was also funded by the Instituto de Salud Carlos III; by grant PI16/00337, cofinanced by the European Regional Development Fund (ERDF); and by the Centro para el Desarrollo Tecnológico Industrial (CDTI), Spanish Ministry of Economics and Competitiveness (MINECO) (grant IDI-20151125). Institutional grants from the Fundación Ramón Areces and Banco Santander to CBMSO are also acknowledged.

The influence of quasispecies dynamics on long-term virus diversification in nature is a largely unexplored question. Specifically, whether intra-host nucleotide and amino acid variation in quasispecies fit the variation observed in consensus sequences or data bank alignments is unknown. Genome conservation and dynamics simulations are used for the computational design of universal vaccines, therapeutic antibodies and pan-genomic antiviral agents. The expectation is that selection of escape mutants will be limited when mutations at conserved residues are required. This strategy assumes long-term (epidemiologically relevant) conservation but, critically, does not consider short-term (quasispecies-dictated) residue conservation. We calculated mutant frequencies of individual loci from mutant spectra of hepatitis C virus (HCV) populations passaged in cell culture and from infected patients. Nucleotide or amino acid conservation in consensus sequences of the same populations, or in the Los Alamos HCV data bank did not match residue conservation in mutant spectra. The results relativize the concept of sequence conservation in viral genetics and suggest that residue invariance in data banks is an insufficient basis for the design of universal viral ligands for clinical purposes. Our calculations suggest relaxed mutational restrictions during quasispecies dynamics, which may contribute to higher calculated short-term than long-term viral evolutionary rates. The research at Centro de Biologia Molecular Severo Ochoa (CSIC-UAM) (CBMSO) was funded by Ministerio de Economia y Competitividad (MINECO), grant numbers SAF2014-52400-R and SAF2017-87846-R, and by Ministerio de Ciencia, Innovacion y Universidades (MCIU), grant number BFU2017-91384-EXP, and by Instituto de Salud Carlos III, grant number PI18/00210, and by PLATESA from Comunidad de Madrid/FEDER, grant numbers S2013/ABI-2906 and S2018/BAA-4370. C.P. is supported by the Miguel Servet program of the Instituto de Salud Carlos III (CP14/00121 and CPII19/00001) cofinanced by the European Regional Development Fund (ERDF). CIBERehd (Centro de Investigacion en Red de Enfermedades Hepaticas y Digestivas) is funded by Instituto de Salud Carlos III. Institutional grants from the Fundacion Ramon Areces and Banco Santander to the CBMSO are also acknowledged. The research at Vall d'Hebron Institut de Recerca (VHIR) was funded by Instituto de Salud Carlos III, cofinanced by the European Regional Development Fund (ERDF), grant number PI19/00301 and by the Centro para el Desarrollo Tecnologico Industrial (CDTI) from the MCIU, grant number IDI-20151125. The research at Centro de Astrobiologia (CAB) was supported by MINECO, grant numbers BIO2016-79618R and PID2019-104903RB-I00 (funded by EU under the FEDER program) and by the Spanish State research agency (AEI), grant number MDM-2017-0737 (Unidad de Excelencia "Maria de Maeztu"-Centro de Astrobiologia (CSIC-INTA). C.G.-C. is supported by predoctoral contract PRE2018-083422 from MCIU. B.M.-G. is supported by predoctoral contract PFIS FI19/00119 from Instituto de Salud Carlos III (Ministerio de Sanidad y Consumo) cofinanced by Fondo Social Europeo (FSE); With funding from the Spanish government through the "María de Maeztu Unit of Excellence" accreditation (MDM-2017-0737) Peer review

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., We thank Dr. Charles M. Rice for the supply of plasmid Jc1FLAG2(p7-nsGluc2A) and helpful advice for the implement of HCV replicon in cell culture, and to Mercedes Guerrero for help with the statistics. Work in Barcelona was funded by Instituto de Salud Carlos III, by grants PI13/00456, PI15/00829, and PI16/00337 cofinanced by the European Regional Development Fund (ERDF), and by CDTI (Centro para el Desarrollo Tecnologico Industrial), Spanish Ministry of Economics and Competitiveness (MINECO), IDI-20151125. C.P. is supported by the Miguel Servet program of the Instituto de Salud Carlos III (CP14/00121) cofinanced by the European Regional Development Fund (ERDF). CIBERehd (Centro de Investigacion en Red de Enfermedades Hepaticas y Digestivas) is funded by Instituto de Salud Carlos III. The work in Madrid was supported by grants BFU-2011-23604, SAF2014-52400-R, SAF2017-87846-R by Spanish MINECO and S2013/ABI-2906 (PLATESA from Comunidad de Madrid/FEDER). Institutional grants from the Fundacion Ramon Areces and Banco Santander to the CBMSO are also acknowledged.