Barrier-Independent, Fitness-Associated Differences in Sofosbuvir Efficacy against Hepatitis C Virus

Synergism between remdesivir and ribavirin leads to SARS-CoV-2 extinction in cell culture

BACKGROUND AND PURPOSE

There is a need for effective anti-COVID-19 treatments, mainly for individuals at risk of severe disease such as the elderly and the immunosuppressed. Drug repositioning has proved effective in identifying drugs that can find a new application for the control of coronavirus disease, in particular COVID-19. The purpose of the present study was to find synergistic antiviral combinations for COVID-19 based on lethal mutagenesis.

EXPERIMENTAL APPROACH

The effect of combinations of remdesivir and ribavirin on the infectivity of SARS-CoV-2 in cell culture has been tested. Viral populations were monitored by ultra-deep sequencing, and the decrease of infectivity as a result of the treatment was measured.

KEY RESULTS

Remdesivir and ribavirin exerted a synergistic inhibitory activity against SARS-CoV-2, quantified both by CompuSyn (Chou-Talalay method) and Synergy Finder (ZIP-score model). In serial passage experiments, virus extinction was readily achieved with remdesivir-ribavirin combinations at concentrations well below their cytotoxic 50 value, but not with the drugs used individually. Deep sequencing of treated viral populations showed that remdesivir, ribavirin, and their combinations evoked significant increases of the number of viral mutations and haplotypes, as well as modification of diversity indices that characterize viral quasi-species.

CONCLUSION AND IMPLICATIONS

SARS-CoV-2 extinction can be achieved by synergistic combination treatments based on lethal mutagenesis. In addition, the results offer prospects of triple drug treatments for effective SARS-CoV-2 suppression.

Fitness-Dependent, Mild Mutagenic Activity of Sofosbuvir for Hepatitis C Virus

The concept of a mild mutagen was coined to describe a minor mutagenic activity exhibited by some nucleoside analogues that potentiated their efficacy as antiretroviral agents. In the present study, we report the mild mutagen activity of sofosbuvir (SOF) for hepatitis C virus (HCV). Serial passages of HCV in human hepatoma cells, in the presence of SOF at a concentration well below its cytotoxic concentration 50 (CC50) led to pre-extinction populations whose mutant spectra exhibited a significant increase of C→U transitions, relative to populations passaged in the absence of SOF. This was reflected in an increase in several diversity indices that were used to characterize viral quasispecies. The mild mutagenic activity of SOF was largely absent when it was tested with isogenic HCV populations that displayed high replicative fitness. Thus, SOF can act as a mild mutagen for HCV, depending on HCV fitness. Possible mechanisms by which the SOF mutagenic activity may contribute to its antiviral efficacy are discussed.

Hepatitis C virus fitness can influence the extent of infection-mediated epigenetic modifications in the host cells

Introduction

Cellular epigenetic modifications occur in the course of viral infections. We previously documented that hepatitis C virus (HCV) infection of human hepatoma Huh-7.5 cells results in a core protein-mediated decrease of Aurora kinase B (AURKB) activity and phosphorylation of Serine 10 in histone H3 (H3Ser10ph) levels, with an affectation of inflammatory pathways. The possible role of HCV fitness in infection-derived cellular epigenetic modifications is not known.

Methods

Here we approach this question using HCV populations that display a 2.3-fold increase in general fitness (infectious progeny production), and up to 45-fold increase of the exponential phase of intracellular viral growth rate, relative to the parental HCV population.

Results

We show that infection resulted in a HCV fitness-dependent, average decrease of the levels of H3Ser10ph, AURKB, and histone H4 tri-methylated at Lysine 20 (H4K20m3) in the infected cell population. Remarkably, the decrease of H4K20m3, which is a hallmark of cellular transformation, was significant upon infection with high fitness HCV but not upon infection with basal fitness virus.

Discussion

Here we propose two mechanisms ─which are not mutually exclusive─ to explain the effect of high viral fitness: an early advance in the number of infected cells, or larger number of replicating RNA molecules per cell. The implications of introducing HCV fitness as an influence in virus-host interactions, and for the course of liver disease, are warranted. Emphasis is made in the possibility that HCV-mediated hepatocellular carcinoma may be favoured by prolonged HCV infection of a human liver, a situation in which viral fitness is likely to increase.

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.

Quasispecies Fitness Partition to Characterize the Molecular Status of a Viral Population. Negative Effect of Early Ribavirin Discontinuation in a Chronically Infected HEV Patient

The changes occurring in viral quasispecies populations during infection have been monitored using diversity indices, nucleotide diversity, and several other indices to summarize the quasispecies structure in a single value. In this study, we present a method to partition quasispecies haplotypes into four fractions according to their fitness: the master haplotype, rare haplotypes at two levels (those present at <0.1%, and those at 0.1−1%), and a fourth fraction that we term emerging haplotypes, present at frequencies >1%, but less than that of the master haplotype. We propose that by determining the changes occurring in the volume of the four quasispecies fitness fractions together with those of the Hill number profile we will be able to visualize and analyze the molecular changes in the composition of a quasispecies with time. To develop this concept, we used three data sets: a technical clone of the complete SARS-CoV-2 spike gene, a subset of data previously used in a study of rare haplotypes, and data from a clinical follow-up study of a patient chronically infected with HEV and treated with ribavirin. The viral response to ribavirin mutagenic treatment was selection of a rich set of synonymous haplotypes. The mutation spectrum was very complex at the nucleotide level, but at the protein (phenotypic/functional) level the pattern differed, showing a highly prevalent master phenotype. We discuss the putative implications of this observation in relation to mutagenic antiviral treatment.

The apparent interferon resistance of transmitted HIV-1 is possibly a consequence of enhanced replicative fitness

HIV-1 transmission via sexual exposure is an inefficient process. When transmission does occur, newly infected individuals are colonized by the descendants of either a single virion or a very small number of establishing virions. These transmitted founder (TF) viruses are more interferon (IFN)-resistant than chronic control (CC) viruses present 6 months after transmission. To identify the specific molecular defences that make CC viruses more susceptible to the IFN-induced 'antiviral state', we established a single pair of fluorescent TF and CC viruses and used arrayed interferon-stimulated gene (ISG) expression screening to identify candidate antiviral effectors. However, we observed a relatively uniform ISG resistance of transmitted HIV-1, and this directed us to investigate possible underlying mechanisms. Simple simulations, where we varied a single parameter, illustrated that reduced growth rate could possibly underly apparent interferon sensitivity. To examine this possibility, we closely monitored in vitro propagation of a model TF/CC pair (closely matched in replicative fitness) over a targeted range of IFN concentrations. Fitting standard four-parameter logistic growth models, in which experimental variables were regressed against growth rate and carrying capacity, to our in vitro growth curves, further highlighted that small differences in replicative growth rates could recapitulate our in vitro observations. We reasoned that if growth rate underlies apparent interferon resistance, transmitted HIV-1 would be similarly resistant to any growth rate inhibitor. Accordingly, we show that two transmitted founder HIV-1 viruses are relatively resistant to antiretroviral drugs, while their matched chronic control viruses were more sensitive. We propose that, when present, the apparent IFN resistance of transmitted HIV-1 could possibly be explained by enhanced replicative fitness, as opposed to specific resistance to individual IFN-induced defences. However, further work is required to establish how generalisable this mechanism of relative IFN resistance might be.

Efficacy decrease of antiviral agents when administered to ongoing hepatitis C virus infections in cell culture

We report a quantification of the decrease of effectiveness of antiviral agents directed to hepatitis C virus, when the agents are added during an ongoing infection in cell culture vs. when they are added at the beginning of the infection. Major determinants of the decrease of inhibitory activity are the time post-infection of inhibitor administration and viral replicative fitness. The efficacy decrease has been documented with antiviral assays involving the combination of the direct-acting antiviral agents, daclatasvir and sofosbuvir, and with the combination of the lethal mutagens, favipiravir and ribavirin. The results suggest that strict antiviral effectiveness assays in preclinical trials may involve the use of high fitness viral populations and the delayed administration of the agents, relative to infection onset.

SARS-CoV-2 Mutant Spectra at Different Depth Levels Reveal an Overwhelming Abundance of Low Frequency Mutations

Populations of RNA viruses are composed of complex and dynamic mixtures of variant genomes that are termed mutant spectra or mutant clouds. This applies also to SARS-CoV-2, and mutations that are detected at low frequency in an infected individual can be dominant (represented in the consensus sequence) in subsequent variants of interest or variants of concern. Here we briefly review the main conclusions of our work on mutant spectrum characterization of hepatitis C virus (HCV) and SARS-CoV-2 at the nucleotide and amino acid levels and address the following two new questions derived from previous results: (i) how is the SARS-CoV-2 mutant and deletion spectrum composition in diagnostic samples, when examined at progressively lower cut-off mutant frequency values in ultra-deep sequencing; (ii) how the frequency distribution of minority amino acid substitutions in SARS-CoV-2 compares with that of HCV sampled also from infected patients. The main conclusions are the following: (i) the number of different mutations found at low frequency in SARS-CoV-2 mutant spectra increases dramatically (50- to 100-fold) as the cut-off frequency for mutation detection is lowered from 0.5% to 0.1%, and (ii) that, contrary to HCV, SARS-CoV-2 mutant spectra exhibit a deficit of intermediate frequency amino acid substitutions. The possible origin and implications of mutant spectrum differences among RNA viruses are discussed.

Enhanced fitness of hepatitis C virus increases resistance to direct-acting antivirals

Drug resistance mutations of hepatitis C virus (HCV) negatively impact viral replicative fitness. RNA viruses are known to change their replication behaviour when subjected to suboptimal selection pressure. Here, we assess whether mutation supply in HCV is sufficiently large to allow the selection of its variants during dual or triple direct-acting antiviral (DAA) treatment associated with augmented virus fitness or impairment. We engineered randomly mutagenized full-genome libraries to create a highly diverse population of replication-competent HCV variants in cell culture. These variants exhibited escape when treated with NS5A/NS5B inhibitors (daclatasvir/sofosbuvir), and relapse on treatment with a combination of NS3/NS5A/NS5B inhibitors (simeprevir or paritaprevir/daclatasvir/sofosbuvir). Analysis of the relationship between virus fitness and drug resistance of JFH1-derived NS5A-5B variants showed a significant positive correlation (P=0.003). At the earliest time points, intracellular RNA levels remain unchanged in both the subgenomic replicon and infection assays, whereas extracellular RNA levels increased upto ten-fold compared to wild-type JFH1. Beneficial substitutions hyperstimulated phosphatidylinositol 4-phosphate during DAA treatment, and showed decreased dependence on cyclophilins during cyclosporine A treatment, indicating an interplay of virus-host molecular mechanisms in beneficial substitution selection that may necessitate infectious virus production. This comprehensive study demonstrates a possible role for HCV fitness of overcoming drug-mediated selection pressure.

Guanosine inhibits hepatitis C virus replication and increases indel frequencies, associated with altered intracellular nucleotide pools

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.

A Two-Level, Intramutant Spectrum Haplotype Profile of Hepatitis C Virus Revealed by Self-Organized Maps

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.

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.

Sofosbuvir Selects for Drug-Resistant Amino Acid Variants in the Zika Virus RNA-Dependent RNA-Polymerase Complex In Vitro

The nucleotide analog sofosbuvir, licensed for the treatment of hepatitis C, recently revealed activity against the Zika virus (ZIKV) in vitro and in animal models. However, the ZIKV genetic barrier to sofosbuvir has not yet been characterized. In this study, in vitro selection experiments were performed in infected human hepatoma cell lines. Increasing drug pressure significantly delayed viral breakthrough (p = 0.029). A double mutant in the NS5 gene (V360L/V607I) emerged in 3 independent experiments at 40–80 µM sofosbuvir resulting in a 3.9 ± 0.9-fold half- maximal inhibitory concentration (IC₅₀) shift with respect to the wild type (WT) virus. A triple mutant (C269Y/V360L/V607I), detected in one experiment at 80 µM, conferred a 6.8-fold IC₅₀ shift with respect to the WT. Molecular dynamics simulations confirmed that the double mutant V360L/V607I impacts the binding mode of sofosbuvir, supporting its role in sofosbuvir resistance. Due to the distance from the catalytic site and to the lack of reliable structural data, the contribution of C269Y was not investigated in silico. By a combination of sequence analysis, phenotypic susceptibility testing, and molecular modeling, we characterized a double ZIKV NS5 mutant with decreased sofosbuvir susceptibility. These data add important information to the profile of sofosbuvir as a possible lead for anti-ZIKV drug development.

Prevalence of Naturally-Occurring NS5A and NS5B Resistance-Associated Substitutions in Iranian Patients With Chronic Hepatitis C Infection

BACKGROUND

Hepatitis C virus (HCV), non-structural 5A (NS5A), and non-structural 5B (NS5B) resistance-associated substitutions (RASs) are the main causes of failure to direct-acting antiviral agents (DAAs). NS5A and NS5B RASs can occur in patients with HCV infection naturally and before exposure to DAAs.

OBJECTIVES

This study aimed to evaluate naturally-occurring NS5A and NS5B RASs in Iranian patients with HCV genotype 1a (HCV-1a) and -3a infections.

METHODS

In this cross-sectional study, viral RNA was extracted from serum specimens. NS5A and NS5B regions were amplified using RT-PCR followed by DNA sequencing. The results of nucleotide sequences were aligned against reference sequences of HCV-1a and -3a and the amino acid substitutions were analyzed using geno2pheno [hcv] web application.

RESULTS

Among 135 patients with hepatitis C, NS5A amino acid substitutions/RASs were identified in 26.4% and 15.9% of patients with HCV-1a and -3a infections, respectively. The identified amino acid substitutions/RASs in the NS5A region of patients with HCV-1a infection were M28T/V/I 11.1%, Q30R/H 4.2%, L31M 1.4%, and H58Y/P/C/D/Q/S/T 16.7%. Y93H substitution was not found in HCV-1a sequences. In patients with HCV-3a infection, NS5A amino acid substitutions/RASs were A30T/K 9.5%, L31F 1.6%, P58S/T/C 3.2%, Y93H 3.2%, and Y93N 3.2%. No resistance substitutions were identified in NS5B sequences from patients with HCV-1a and -3a infections.

CONCLUSION

In this study, baseline amino acid substitutions/RASs were only identified in the NS5A region in Iranian patients with HCV-1a and -3a infections, and the prevalence of these amino acid substitutions/RASs were in accordance with similar studies. There were no RASs in the HCV-1a and -3a NS5B region.

Amino Acid Substitutions Associated with Treatment Failure for Hepatitis C Virus Infection

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.

Virological patterns of hepatitis C virus patients with failure to the current-generation direct-acting antivirals

There are few data on the virological characterisation of patients with failure to current-generation direct-acting antivirals (DAAs), namely elbasvir/grazoprevir, sofosbuvir/velpatasvir and glecaprevir/pibrentasvir. This study aimed to characterise virological patterns in patients with failure to current DAA regimens as well as the efficacy of re-treatment. All 61 consecutive hepatitis C virus (HCV) treatment-naïve patients with failure to current DAAs from January 2018 to February 2019 were enrolled. Sanger sequencing of NS3, NS5A and NS5B proteins was performed using homemade protocols. NS5A resistance-associated substitutions (RASs) were more frequent in the 17 patients treated with sofosbuvir/velpatasvir (89.5%) and 33 patients treated with elbasvir/grazoprevir (97%) compared with the 11 patients treated with glecaprevir/pibrentasvir (18.2%) (P = 0.002 and 0.000, respectively). NS3 RASs were more often detected in the 33 patients with failure to elbasvir/grazoprevir (30.3%) than in the 11 patients treated with glecaprevir/pibrentasvir (9.1%). NS3 RASs were also detected in 12% of sofosbuvir/velpatasvir-treated patients. NS5B RASs were infrequently identified. Of the glecaprevir/pibrentasvir-treated patients, 73% did not show RASs in any HCV regions, a prevalence higher than that observed in those treated with elbasvir/grazoprevir (0%; P < 0.05) or sofosbuvir/velpatasvir (12%; P < 0.05). Of the 61 patients, 21 (34.4%) were re-treated with sofosbuvir/velpatasvir and voxilaprevir. All patients achieved sustained virological response at 12 weeks (SVR12). To our knowledge, this is one of the first real-life studies describing patients who failed current-generation DAAs; the prevalence of RASs differed according to the DAA regimen used, and the efficacy of re-treatment was high.

Broad and Dynamic Diversification of Infectious Hepatitis C Virus in a Cell Culture Environment

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.

Viral fitness: history and relevance for viral pathogenesis and antiviral interventions

The quasispecies dynamics of viral populations (continuous generation of variant genomes and competition among them) has as one of its frequent consequences variations in overall multiplication capacity, a major component of viral fitness. This parameter has multiple implications for viral pathogenesis and viral disease control, some of them unveiled thanks to deep sequencing of viral populations. Darwinian fitness is an old concept whose quantification dates back to the early developments of population genetics. It was later applied to viruses (mainly to RNA viruses) to quantify relative multiplication capacities of individual mutant clones or complex populations. The present article reviews the fitness concept and its relevance for the understanding of the adaptive dynamics of viruses in constant and changing environments. Many studies have addressed the fitness cost of escape mutations (to antibodies, cytotoxic T cells or inhibitors) as an influence on the efficacy of antiviral interventions. Here, we summarize the evidence that the basal fitness level can be a determinant of inhibitor resistance.

Deep-sequencing reveals broad subtype-specific HCV resistance mutations associated with treatment failure

A percentage of hepatitis C virus (HCV)-infected patients fail direct acting antiviral (DAA)-based treatment regimens, often because of drug resistance-associated substitutions (RAS). The aim of this study was to characterize the resistance profile of a large cohort of patients failing DAA-based treatments, and investigate the relationship between HCV subtype and failure, as an aid to optimizing management of these patients. A new, standardized HCV-RAS testing protocol based on deep sequencing was designed and applied to 220 previously subtyped samples from patients failing DAA treatment, collected in 39 Spanish hospitals. The majority had received DAA-based interferon (IFN) α-free regimens; 79% had failed sofosbuvir-containing therapy. Genomic regions encoding the nonstructural protein (NS) 3, NS5A, and NS5B (DAA target regions) were analyzed using subtype-specific primers. Viral subtype distribution was as follows: genotype (G) 1, 62.7%; G3a, 21.4%; G4d, 12.3%; G2, 1.8%; and mixed infections 1.8%. Overall, 88.6% of patients carried at least 1 RAS, and 19% carried RAS at frequencies below 20% in the mutant spectrum. There were no differences in RAS selection between treatments with and without ribavirin. Regardless of the treatment received, each HCV subtype showed specific types of RAS. Of note, no RAS were detected in the target proteins of 18.6% of patients failing treatment, and 30.4% of patients had RAS in proteins that were not targets of the inhibitors they received. HCV patients failing DAA therapy showed a high diversity of RAS. Ribavirin use did not influence the type or number of RAS at failure. The subtype-specific pattern of RAS emergence underscores the importance of accurate HCV subtyping. The frequency of "extra-target" RAS suggests the need for RAS screening in all three DAA target regions.

Quasispecies dynamics in disease prevention and control

Medical interventions to prevent and treat viral disease constitute evolutionary forces that may modify the genetic composition of viral populations that replicate in an infected host and influence the genomic composition of those viruses that are transmitted and progress at the epidemiological level. Given the adaptive potential of viruses in general and the RNA viruses in particular, the selection of viral mutants that display some degree of resistance to inhibitors or vaccines is a tangible challenge. Mutant selection may jeopardize control of the viral disease. Strategies intended to minimize vaccination and treatment failures are proposed and justified based on fundamental features of viral dynamics explained in the preceding chapters. The recommended use of complex, multiepitopic vaccines, and combination therapies as early as possible after initiation of infection falls under the general concept that complexity cannot be combated with simplicity. It also follows that sociopolitical action to interrupt virus replication and spread as soon as possible is as important as scientifically sound treatment designs to control viral disease on a global scale.

Synergistic lethal mutagenesis of hepatitis C virus

Lethal mutagenesis is an antiviral approach that consists in 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 that renders the strategy potentially effective against emergent RNA viral infections. Here we describe 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 as implemented in the CompuSyn graphics, with average dose reduction index (DRI) above 1 (68.02±101.6 for favipiravir, and 5.83±6.07 for ribavirin), and average combination indices (CI) below 1 (0.52±0.28). Furthermore, analogue concentrations that individually did not extinguish high fitness HCV in ten serial infections, when used in combination they extinguished high fitness HCV in one to two passages. Although both analogues display 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. Prospects of synergy among mutagenic nucleotides as a strategy to confront emerging viral infections are discussed.

Worldwide prevalence of baseline resistance-associated polymorphisms and resistance mutations in HCV against current direct-acting antivirals

BACKGROUND

HCV infections can now be completely cured, thanks to the currently marketed direct-acting antivirals (DAAs). It is known that HCV patients carry viral populations with baseline polymorphisms and/or mutations that make them resistant against some of these DAAs, which can negatively impact the patient's treatment outcome. Using complete HCV coding sequences isolated from 1,306 treatment-naive patients of genotypes (GTs) 1, 2, 3, 4 and 6 from around the globe, we studied the prevalence of baseline resistance-associated polymorphisms (RAPs) and resistance mutations (RMs) against DAAs that are currently on the market or in clinical trials.

METHODS

The HCV genome sequences used in this study were retrieved from the NCBI database. RAPs and RMs, with reference to HCV GT1a, were identified using the HCV Geno2pheno web server.

RESULTS

Nearly 50% of the total amino acid positions (including NS3 protease, NS5A and NS5B) studied are baseline polymorphisms that differentiated one GT from the rest. A proportion of these baseline polymorphisms and baseline non-polymorphic RMs could confer a significant increase in resistance against DAAs.

CONCLUSIONS

In this study, we show the presence and prevalence of RAPs and RMs in DAA treatment-naive patients against currently used DAAs or DAAs in clinical trials. Our study suggests that RAPs and RMs profiling of HCV patients should be performed before the start of the therapy. Our results should be relevant especially in low- and middle-income countries, where the patients have a large variation of GTs and subtypes, and where the generic HCV treatment is now increasingly available.

Management of HCV patients in cases of direct-acting antiviral failure

Introduction: Over 70 million people are infected with hepatitis C virus (HCV), increasing the risk of cirrhosis and hepatocellular carcinoma. Areas covered: Since the approval of the first interferon-free direct-acting antiviral (DAA) therapy in 2011, a number of DAAs have been approved, and HCV is now considered curable. Until recently, however, there were no clear guidelines on how to re-treat patients who fail DAA therapy. Current protease inhibitors (PIs) are generally unaffected by earlier resistance-associated variants (RAVs), but many NS5A inhibitors continue to have overlapping resistance profiles, and NS5A RAVs can persist even in the absence of DAAs. Expert opinion: Fortunately, RAVs affecting NS5B polymerase inhibitors are rare, making sofosbuvir a safe choice as the backbone of re-treatment therapies. Recent re-treatment guidelines that take into account genotype, fibrosis, treatment history, and RAV suggest that >90% of patients with prior treatment failures can be successfully re-treated with sofosbuvir/velpatasvir, sofosbuvir/velpatasvir/voxilaprevir or glecaprevir/pibrentasvir.

The increasing impact of lethal mutagenesis of viruses

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

A proof-of-concept study in HCV-infected Huh7.5 cells for shortening the duration of DAA-based triple treatment regimens

With the development of more effective direct-acting antivirals (DAAs), dual- or triple-therapy regimens represent the major strategy used to cure chronic hepatitis C virus (HCV) infection. Thus, shorter treatment duration regimens with low burden, few adverse effects and good patient adherence are urgently needed. This study theoretically demonstrates a proof-of-concept approach for shortening therapy duration by examining HCV-infected Huh7.5 cells after treatment with a high or low fixed dose of three DAAs (simeprevir + daclatasvir + sofosbuvir) for 6-15 days. The results demonstrated that HCV-infected Huh7.5 cells achieved an ultrarapid virologic response with undetectable HCV RNA and protein and were cured after treatment with the triple-therapy regimen for 15 days. When the treatment duration was shortened, virologic relapse might occur after treatment with a low fixed dose of the three DAAs for 9 days and did occur after treatment with a low fixed dose for 6 days, although HCV was below detectable levels at the end of treatment. However, virologic relapse could be avoided with treatment of a high fixed dose of the three DAAs for 9 or 6 days. Although a virologic breakthrough occurred after an intermittent treatment regimen at the low fixed dose, the high fixed dose cured HCV-positive Huh7.5 cells with intermittent treatment. In conclusion, HCV is persistently present below detectable levels in HCV-infected Huh7.5 cells for a long time after treatment, and a shortened therapy duration is associated with an increased risk of virologic relapse, but virologic relapse or breakthrough might be avoided by treatment with a combination of more highly effective DAAs.

Toll-Like Receptor 7 (TLR-7) and TLR-9 Agonists Improve Hepatitis C Virus Replication and Infectivity Inhibition by Plasmacytoid Dendritic Cells

Plasmacytoid dendritic cells (pDCs) are innate immune cells with high antiviral activity triggered by Toll-like receptor 7 (TLR-7) and TLR-9 stimulation. Moreover, they are important mediators between innate and adaptive immunity. Although nowadays there is available an effective therapeutic arsenal against hepatitis C virus (HCV), a protective vaccine is not available. We have analyzed the pDCs' response to HCV infection in a hepatitis C virus (HCV)-Huh7.5 virus-cell system, which allows completion of the virus infectious cycle. pDCs were cocultured following human immunodeficiency virus (HIV) aldrithiol-2 (AT-2 [TLR-7 agonist]) inactivation and CpG (TLR-9 agonist) stimulation. We employed three virus derivatives-wild-type Jc1, interferon (IFN)-resistant virus IR, and high-replicative-fitness virus P100-in order to explore additional IFN-α-related virus inhibition mechanisms. pDCs inhibited HCV infectivity and replication and produced IFN-α. After TLR-7 and TLR-9 stimulation, inhibition of infectivity and IFN-α production by pDCs were enhanced. TLR-7 stimulation drove higher TNF-related apoptosis-inducing ligand (TRAIL) expression in pDCs. Additionally, TLR-7- and TLR-9-stimulated pDCs exhibited a mature phenotype, improving the antigen presentation and lymph node homing-related markers. In conclusion, pDCs could serve as a drug target against HCV in order to improve antiviral activity and as an enhancer of viral immunization.IMPORTANCE We implemented a coculture system of pDCs with HCV-infected hepatoma cell line, Huh7.5. We used three HCV derivatives in order to gain insight into pDCs' behavior against HCV and associated antiviral mechanisms. The results with this cell coculture system support the capacity of pDCs to inhibit HCV replication and infectivity mainly via IFN-α, but also through additional mechanisms associated with pDC maturation. We provided evidence that TLR agonists can enhance antiviral pDCs' function and can induce phenotypic changes that may facilitate the interplay with other immune cells. These findings suggest the possibility of including TLR agonists in the strategies of HCV vaccine development.

Hepatitis C virus early kinetics and resistance-associated substitution dynamics during antiviral therapy with direct-acting antivirals

The emergence of resistance-associated substitutions (RASs) can compromise the high efficacy of direct-acting antivirals (DAAs). Little is known about RASs selection at very early time points during DAA treatment. Therefore, we analyzed the potential emergence of RASs immediately after therapy initiation. Samples of 71 patients treated with different DAAs were collected at baseline, during therapy (hours 4 and 8; days 1-7; weeks 2-4) or until target not detected. HCV-RNA levels were determined by qPCR, and RASs were detected by deep sequencing. Sixty-three (89%) patients achieved a sustained virological response (SVR), 7 (10%) relapsed, and 1 (1%) experienced a breakthrough. Almost all non-SVR (7/8, 88%) showed RASs either at baseline or relapse. High-frequency RASs detected at baseline (Y93H and L159F+C316N) remained detectable at early time points during therapy and reappeared as most prevalent substitutions at relapse. Conversely, emergent RASs at relapse (Q80R, D168E/V, R155K and L31V) were not observed during the first hours-days, before HCV-RNA became undetectable. HCV-RNA decay and genetic evolution of the quasispecies followed a similar pattern during the first hours of therapy in SVR and non-SVR patients. In conclusion, the absence of early RASs selection and the similar dynamics of HCV kinetics and quasispecies in SVR and non-SVR patients after therapy initiation suggest that RASs selection may occur at later stages in the remaining reservoir, where viral populations persist hidden at very low replication levels. Nevertheless, we cannot completely exclude very early selection, when RASs are present below the sensitivity limit of deep sequencing.

Quasispecies dynamics and clinical significance of hepatitis C virus (HCV) antiviral resistance

Hepatitis C virus (HCV) follows quasispecies dynamics in infected hosts and this influences its biology, how the virus diversifies into several genotypes and many subtypes, and how viral populations respond to antiviral therapies. Despite current antiviral combinations being able to cure a great percentage of HCV-infected patients, the presence of resistance-associated substitutions (RASs) diminishes the success of antiviral therapies, which is a main concern in the re-treatment of patients treated with direct-acting antiviral agents. Current methodologies such as ultra deep sequencing are ideal tools to obtain a detailed representation of the mutant spectrum composition circulating in infected patients. Such knowledge should allow optimisation of rescue treatments. A new mechanism of antiviral resistance not based on the selection of RASs but on high viral fitness 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.

Resistance of high fitness hepatitis C virus to lethal mutagenesis

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

Quasispecies and virus

Quasispecies theory has been instrumental in the understanding of RNA virus population dynamics because it considered for the first time mutation as an integral part of the replication process. The key influences of quasispecies theory on experimental virology have been: (1) to disclose the mutant spectrum nature of viral populations and to evaluate its consequences; (2) to unveil collective properties of genome ensembles that can render a mutant spectrum a unit of selection; and (3) to identify new vulnerability points of pathogenic RNA viruses on three fronts: the need to apply multiple selective constraints (in the form of drug combinations) to minimize selection of treatment-escape variants, to translate the error threshold concept into antiviral designs, and to construct attenuated vaccine viruses through alterations of viral polymerase copying fidelity or through displacements of viral genomes towards unfavorable regions of sequence space. These three major influences on the understanding of viral pathogens preceded extensions of quasispecies to non-viral systems such as bacterial and tumor cell collectivities and prions. These developments are summarized here.

Deep sequencing in the management of hepatitis virus infections

The hepatitis viruses represent a major public health problem worldwide. Procedures for characterization of the genomic composition of their populations, accurate diagnosis, identification of multiple infections, and information on inhibitor-escape mutants for treatment decisions are needed. Deep sequencing methodologies are extremely useful for these viruses since they replicate as complex and dynamic quasispecies swarms whose complexity and mutant composition are biologically relevant traits. Population complexity is a major challenge for disease prevention and control, but also an opportunity to distinguish among related but phenotypically distinct variants that might anticipate disease progression and treatment outcome. Detailed characterization of mutant spectra should permit choosing better treatment options, given the increasing number of new antiviral inhibitors available. In the present review we briefly summarize our experience on the use of deep sequencing for the management of hepatitis virus infections, particularly for hepatitis B and C viruses, and outline some possible new applications of deep sequencing for these important human pathogens.

Internal Disequilibria and Phenotypic Diversification during Replication of Hepatitis C Virus in a Noncoevolving Cellular Environment

Viral quasispecies evolution upon long-term virus replication in a noncoevolving cellular environment raises relevant general issues, such as the attainment of population equilibrium, compliance with the molecular-clock hypothesis, or stability of the phenotypic profile. Here, we evaluate the adaptation, mutant spectrum dynamics, and phenotypic diversification of hepatitis C virus (HCV) in the course of 200 passages in human hepatoma cells in an experimental design that precluded coevolution of the cells with the virus. Adaptation to the cells was evidenced by increase in progeny production. The rate of accumulation of mutations in the genomic consensus sequence deviated slightly from linearity, and mutant spectrum analyses revealed a complex dynamic of mutational waves, which was sustained beyond passage 100. The virus underwent several phenotypic changes, some of which impacted the virus-host relationship, such as enhanced cell killing, a shift toward higher virion density, and increased shutoff of host cell protein synthesis. Fluctuations in progeny production and failure to reach population equilibrium at the genomic level suggest internal instabilities that anticipate an unpredictable HCV evolution in the complex liver environment.IMPORTANCE Long-term virus evolution in an unperturbed cellular environment can reveal features of virus evolution that cannot be explained by comparing natural viral isolates. In the present study, we investigate genetic and phenotypic changes that occur upon prolonged passage of hepatitis C virus (HCV) in human hepatoma cells in an experimental design in which host cell evolutionary change is prevented. Despite replication in a noncoevolving cellular environment, the virus exhibited internal population disequilibria that did not decline with increased adaptation to the host cells. The diversification of phenotypic traits suggests that disequilibria inherent to viral populations may provide a selective advantage to viruses that can be fully exploited in changing environments.

Lethal Mutagenesis of Hepatitis C Virus Induced by Favipiravir

Lethal mutagenesis is an antiviral approach that consists in extinguishing a virus by an excess of mutations acquired during replication in the presence of a mutagen. Here we show that favipiravir (T-705) is a potent mutagenic agent for hepatitis C virus (HCV) during its replication in human hepatoma cells. T-705 leads to an excess of G → A and C → U transitions in the mutant spectrum of preextinction HCV populations. Infectivity decreased significantly in the presence of concentrations of T-705 which are 2- to 8-fold lower than its cytotoxic concentration 50 (CC50). Passaging the virus five times in the presence of 400 μM T-705 resulted in virus extinction. Since T-705 has undergone advanced clinical trials for approval for human use, the results open a new approach based on lethal mutagenesis to treat hepatitis C virus infections. If proven effective for HCV in vivo, this new anti-HCV agent may be useful in patient groups that fail current therapeutic regimens.

My Cousin, My Enemy: quasispecies suppression of drug resistance

If a freshly minted genome contains a mutation that confers drug resistance, will it be selected in the presence of the drug? Not necessarily. During viral infections, newly synthesized viral genomes occupy the same cells as parent and other progeny genomes. If the antiviral target is chosen so that the drug-resistant progeny's growth is dominantly inhibited by the drug-susceptible members of its intracellular family, its outgrowth can be suppressed. Precedent for 'dominant drug targeting' as a deliberate approach to suppress the outgrowth of inhibitor-resistant viruses has been established for envelope variants of vesicular stomatitis virus and for capsid variants of poliovirus and dengue virus. Small molecules that stabilize oligomeric assemblages are a promising means to an unfit family to destroy the effectiveness of a newborn drug-resistant relative due to the co-assembly of drug-susceptible and drug-resistant monomers.