Search

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.

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.