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

Contemporary Insights into Hepatitis C Virus: A Comprehensive Review

Hepatitis C virus (HCV) remains a significant global health challenge. Approximately 50 million people were living with chronic hepatitis C based on the World Health Organization as of 2024, contributing extensively to global morbidity and mortality. The advent and approval of several direct-acting antiviral (DAA) regimens significantly improved HCV treatment, offering potentially high rates of cure for chronic hepatitis C. However, the promising aim of eventual HCV eradication remains challenging. Key challenges include the variability in DAA access across different regions, slightly variable response rates to DAAs across diverse patient populations and HCV genotypes/subtypes, and the emergence of resistance-associated substitutions (RASs), potentially conferring resistance to DAAs. Therefore, periodic reassessment of current HCV knowledge is needed. An up-to-date review on HCV is also necessitated based on the observed shifts in HCV epidemiological trends, continuous development and approval of therapeutic strategies, and changes in public health policies. Thus, the current comprehensive review aimed to integrate the latest knowledge on the epidemiology, pathophysiology, diagnostic approaches, treatment options and preventive strategies for HCV, with a particular focus on the current challenges associated with RASs and ongoing efforts in vaccine development. This review sought to provide healthcare professionals, researchers, and policymakers with the necessary insights to address the HCV burden more effectively. We aimed to highlight the progress made in managing and preventing HCV infection and to highlight the persistent barriers challenging the prevention of HCV infection. The overarching goal was to align with global health objectives towards reducing the burden of chronic hepatitis, aiming for its eventual elimination as a public health threat by 2030.

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.

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.

Composition and Dynamics of H1N1 and H7N9 Influenza A Virus Quasispecies in a Co-infected Patient Analyzed by Single Molecule Sequencing Technology

A human co-infected with H1N1 and H7N9 subtypes influenza A virus (IAV) causes a complex infectious disease. The identification of molecular-level variations in composition and dynamics of IAV quasispecies will help to understand the pathogenesis and provide guidance for precision medicine treatment. In this study, using single-molecule real-time sequencing (SMRT) technology, we successfully acquired full-length IAV genomic sequences and quantified their genotypes abundance in serial samples from an 81-year-old male co-infected with H1N1 and H7N9 subtypes IAV. A total of 26 high diversity nucleotide loci was detected, in which the A-G base transversion was the most abundant substitution type (67 and 64%, in H1N1 and H7N9, respectively). Seven significant amino acid variations were detected, such as NA:H275Y and HA: R222K in H1N1 as well as PB2:E627K and NA: K432E in H7N9, which are related to viral drug-resistance or mammalian adaptation. Furtherly, we retrieved 25 H1N1 and 22 H7N9 genomic segment haplotypes from the eight samples based on combining high-diversity nucleotide loci, which provided a more concise overview of viral quasispecies composition and dynamics. Our approach promotes the popularization of viral quasispecies analysis in a complex infectious disease, which will boost the understanding of viral infections, pathogenesis, evolution, and precision medicine.

Study of Quasispecies Complexity and Liver Damage Progression after Liver Transplantation in Hepatitis C Virus Infected Patients

Cirrhosis derived from chronic hepatitis C virus (HCV) infection is still a common indication for liver transplantation (LT). Reinfection of the engrafted liver is universal in patients with detectable viral RNA at the time of transplant and causes fast progression of cirrhosis (within 5 years) in around one-third of these patients. To prevent damage to the liver graft, effective direct-acting antiviral (DAA) therapy is required as soon as possible. However, because of post-LT clinical instability, it is difficult to determine the optimal time to start DAAs with a low risk of complications. Evaluate changes in quasispecies complexity following LT and seek a predictive index of fast liver damage progression to determine the timing of DAA initiation. HCV genomes isolated from pre-LT and 15-day post-LT serum samples of ten patients, who underwent orthotopic LT, were quantified and sequenced using a next-generation sequencing platform. Sequence alignments, phylogenetic trees, quasispecies complexity measures, biostatistics analyses, adjusted R2 values, and analysis of variance (ANOVA) were carried out. Three different patterns of reinfection were observed (viral bottlenecking, conserved pre-LT population, and mixed populations), suggesting that bottlenecking or homogenization of the viral population is not a generalized effect after liver graft reinfection. None of the quasispecies complexity measures predicted the future degree of liver damage. Higher and more uniform viral load (VL) values were observed in all pre-LT samples, but values were more dispersed in post-LT samples. However, VL increased significantly from the pre-LT to 15-day post-LT samples in patients with advanced fibrosis at 1-year post-LT, suggesting that a VL increase on day 15 may be a predictor of fast liver fibrosis progression. HCV kinetics after LT differ between patients and are not fibrosis-dependent. Higher VL at day 15 post-LT versus pre-LT samples may predict fast liver fibrosis progression.

Late Relapse and Reinfection in HCV Patients Treated with Direct-Acting Antiviral (DAA) Drugs

The risk of hepatitis C virus (HCV) recurrence after direct-acting antiviral (DAA) treatment is <0.5%. However, the distinction between HCV RNA late relapse and reinfection still represents a challenge in virological diagnostics. The aim of this study was to employ next-generation sequencing (NGS) to investigate HCV RNA recurrence in patients achieving a sustained virologic response (SVR) at least six months post-treatment. NGS was performed on plasma samples from six HCV-positive patients (Pt1–6) treated with DAA. NGS of HCV NS5B was analyzed before treatment (T0), after HCV RNA rebound (T1), and, for Pt3, after a second rebound (T2). Reinfection was confirmed for Pt5, and for the first rebound observed in Pt3. Conversely, viral relapse was observed when comparing T0 and T1 for Pt6 and T1 and T2 for Pt3. Z-scores were calculated and used to predict whether HCV-positive patient samples at different time points belonged to the same quasispecies population. A low Z-score of <2.58 confirmed that viral quasispecies detected at T0 and T1 were closely related for both Pt1 and Pt2, while the Z-score for Pt4 was suggestive of possible reinfection. NGS data analyses indicate that the Z-score may be a useful parameter for distinguishing late relapse from reinfection.

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.

Therapy Implications of Hepatitis C Virus Genetic Diversity

Hepatitis C virus (HCV) is an important human pathogen with a high chronicity rate. An estimated 71 million people worldwide are living with chronic hepatitis C (CHC) infection, which carries the risk of progression to hepatic fibrosis, cirrhosis, and hepatocellular carcinoma (HCC). Similar to other RNA viruses, HCV has a high rate of genetic variability generated by its high mutation rate and the actions of evolutionary forces over time. There are two levels of HCV genetic variability: intra-host variability, characterized by the distribution of HCV mutant genomes present in an infected individual, and inter-host variability, represented by the globally circulating viruses that give rise to different HCV genotypes and subtypes. HCV genetic diversity has important implications for virus persistence, pathogenesis, immune responses, transmission, and the development of successful vaccines and antiviral strategies. Here we will discuss how HCV genetic heterogeneity impacts viral spread and therapeutic control.

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 quasispecies

Viral quasispecies refers to a population structure that consists of extremely large numbers of variant genomes, termed mutant spectra, mutant swarms or mutant clouds. Fueled by high mutation rates, mutants arise continually, and they change in relative frequency as viral replication proceeds. The term quasispecies was adopted from a theory of the origin of life in which primitive replicons) consisted of mutant distributions, as found experimentally with present day RNA viruses. The theory provided a new definition of wild type, and a conceptual framework for the interpretation of the adaptive potential of RNA viruses that contrasted with classical studies based on consensus sequences. Standard clonal analyses and deep sequencing methodologies have confirmed the presence of myriads of mutant genomes in viral populations, and their participation in adaptive processes. The quasispecies concept applies to any biological entity, but its impact is more evident when the genome size is limited and the mutation rate is high. This is the case of the RNA viruses, ubiquitous in our biosphere, and that comprise many important pathogens. In virology, quasispecies are defined as complex distributions of closely related variant genomes subjected to genetic variation, competition and selection, and that may act as a unit of selection. Despite being an integral part of their replication, high mutation rates have an upper limit compatible with inheritable information. Crossing such a limit leads to RNA virus extinction, a transition that is the basis of an antiviral design termed lethal mutagenesis.

Subtype-Specific Prevalence of Hepatitis C Virus NS5A Resistance Associated Substitutions in Mainland China

Resistance associated substitutions (RASs) can reduce the efficacy of direct-acting antiviral agents (DAAs) targeting hepatitis C virus (HCV) and lead to treatment failure. Clinical data of HCV NS5A RASs prevalence are limited in China and need to be investigated. A total of 878 unique patient samples with different genotypes (GT) (1b: n = 489, 2a: n = 203, 3a: n = 60, 3b: n = 78, 6a: n = 48) were collected from around mainland China by KingMed Laboratory and analyzed for NS5A RASs distribution by Sanger sequencing. Phylogeographic analyses based on NS5A domain 1 sequences indicated circulation of both locally and nationally epidemic strains. Relatively high frequency of Y93H (14.1%) was only detected in GT1b but not in other subtypes. High frequency of L31M was found in both GT2a (95.6%) and GT3b (98.7%) sequences. Due to the overlapping incidence of A30K, 96% of GT3b isolates had NS5A RASs combination A30K + L31M, which confers high levels of resistance to most NS5A inhibitors. No RASs were detected in GT6a strains. Meanwhile, baseline NS5A RASs fingerprints were also evaluated in 185 DAA treatment-naive GT1b patients with next generation sequencing method. Patients presenting with Y93H had statistically higher entropy of HCV NS5A sequences. Taken together, subtype-specific distribution patterns of NS5A RASs were observed. GT1b patients with higher HCV complexity tend to have a greater chance of Y93H presence, while GT3b patients are naturally resistant to current NS5A inhibitors and their treatment may pose a challenge to real-world DAA application.

Prevalence and Factors Related to Natural Resistance-Associated Substitutions to Direct-Acting Antivirals in Patients with Genotype 1 Hepatitis C Virus Infection

This study aimed to assess the prevalence of natural resistance-associated substitutions (RASs) to NS3, NS5A and NS5B inhibitors in 86 genotype 1 Hepatitis C Virus (HCV)-infected patients from Buenos Aires, Argentina, and to determine their effect on therapy outcome. Additionally, virological, clinical and host genetic factors were explored as predictors of the presence of baseline RASs. NS3 RASs (39.2%) were more prevalent than NS5A RASs (25%) and NS5B RASs (8.9%). In the three regions, the frequencies of RASs were significantly higher in HCV-1b than in HCV-1a. The prevalence of Y93H, L159F and Q80K were 1.3%, 6.3% and 2.5%, respectively. IFNL3 CC genotype was identified as an independent predictor of the presence of baseline RASs in NS5A and NS3 genes (p = 0.0005 and p = 0.01, respectively). Sustained virologic response was achieved by 93.3% of the patients after receiving direct-acting antivirals (DAAs), although 48.7% of them showed baseline RASs related to the DAA-regimen. Notably, the prevalence of clinically relevant RASs in the three genes was lower than that observed around the world. The baseline presence of RASs in both subtypes did not appear to affect therapy outcome. These results support the need to evaluate resistance patterns in each particular country since RASs´ prevalence significantly vary worldwide.