Hepatitis C Virus (HCV) NS3 sequence diversity and antiviral resistance-associated variant frequency in HCV/HIV coinfection

Dynamic changes of HCV genomes under selective pressure from DAAs therapy in relapsed patients

Currently, direct-acting antiviral drugs (DAAs) are widely used as therapeutic methods for hepatitis C virus (HCV)-positive patients, however, patients may experience treatment failure, and the dynamic changes of HCV genomes in these patients are unknown. In this study, three real-world DAAs cohorts were enrolled to observe clinical efficacy. In addition, serum samples from treatment failure patients at baseline and relapse were used to analyze changes of the HCV genomes at near full-length genome level, including resistance-associated variants (RAVs), viral quasispecies diversity and selection analysis. Next-generation sequencing was used as the detection method. The overall sustained virological response at 12 w after the end of treatment was achieved in 91.9% (57/62) of HCV patients, and 3 paired samples obtained from relapsed patients. All the 3 patients harbored baseline NS5A RAVs, the frequency of NS5A RAVs increased in 2 patients and a new NS5A RAV emerged in 1 patient at relapse, and almost all the viral strains existed with NS5A RAVs at relapse. The results of the viral quasispecies diversity analysis revealed that viral quasispecies diversity decreased at relapse compared to baseline, and the results of selection analysis indicated that the virus population experienced a bottleneck phenomenon, recent selective sweep and population expansion or was under purification selection after DAAs treatment. This study indicated that the clinical efficacy was excellent in real-world DAAs cohorts, and the viral strains existed at relapse were selective by DAAs therapy.

Baseline Hepatitis C Virus NS5A Resistance-Associated Polymorphisms in Patients With and Without Human Immunodeficiency Virus Coinfection in Mexico

Objective: We aimed to evaluate the frequency and associated factors of baseline NS5A resistance-associated substitutions (RASs) in patients coinfected with human immunodeficiency virus (HIV) and hepatitis C virus (HCV) monoinfection with genotype 1b (GT1b) or genotype 1a (GT1a). Moreover, we performed a phylogenetic analysis to evaluate the pattern of clustering among samples of patients with RASs. Results: Fifty-five patients were infected with GT1a, of whom 44 (80%) were HIV-infected patients. RAS prevalence in GT1a was 14% (6/44) and distributed as follows: 5 (11%) harbored M28V and 1 (2%) A92T. Twenty-four patients were infected with HCV GT1b, of whom only 5 (21%) were HIV coinfected; RASs were found in 17/24 (71%) patients, as follows: Y93H+F37L+Q54H (1/24), Y93H+F37L (1/24), P58S (1/24), L31F+F37L (1/24), F37L+H/Q54H (3/24), and F37L (10/24). Only GT1b was significantly associated with RASs (adjusted odds ratio 16.37; 95% confidence interval 2.74-97.48; p = 0.002) in the multivariate analysis. A cluster of sequences from HIV/HCV GT1a patients was found; however, we did not find phylogenetic relationships among sequences with NS5A RASs. Conclusions: In our population of HCV-infected patients, the frequency of NS5A RASs at baseline was somewhat similar to the previously reported worldwide rate. HCV GT1b showed the most significant association with harboring of NS5A RASs. Of note, despite there being clusters among sequences of HIV-coinfected patients, NS5A RASs were not transmitted.

Naturally occurring NS5A and NS5B resistant associated substitutions in HCV and HCV/HIV patients in iranian population

BACKGROUND

The introduction of direct acting antivirals (DAAs) for hepatitis C virus (HCV) treatment promises shorter treatment duration, higher cure rates and fewer side effects. Naturally, occurring Resistance Associated Substitutions (RASs) are major challenge to the success of the HCV antiviral therapy.

AIM

To determine the naturally occurring NS5A and NS5B RASs in Iranian HCV and HCV/human immunodeficiency virus (HIV) patients.

METHODS

A total of 209 DAA-naïve chronic HCV patients including 104 HCV mono-infected and 105 HCV/HIV co-infected cases were enrolled. Amplification and Sanger population sequencing of NS5A and NS5B regions of HCV genome were carried out. The amino acid sequence diversity of the NS5A and NS5B regions were analyzed using geno2pheno HCV.

RESULTS

NS5A RASs were detected in 25.5% of HCV and 16.9% of HCV/HIV subjects. In HCV cases, clinically relevant RASs were L28M followed by M28Vand Q30H and Y93H/N. In HCV/HIV subjects, clinically relevant RASs were Y93H/N followed by L28M and P58T and M28V/T and Q30R. NS5B RASs were observed in 11.8% of HCV and 5.9% of HCV/HIV subjects. Clinically relevant substitutions were included V321A/I, C316Y, S282R and L159F. The major S282T mutation was not observed.

CONCLUSION

The emergence of RASs is a growing issue in the setting of current treatment with DAAs. Although currently, screening of RASs is recommended before specific DAA regimens, it should be consider in patients with therapeutic failure and in the cases of retreatment.

Favorable outcome with direct-acting antiviral treatment in hepatitis C patients coinfected with HIV

AIM

To investigate the efficacy and safety of all-oral direct-acting antiviral treatments in patients coinfected with hepatitis C virus (HCV) and HIV.

METHODS

In all, 35 patients with HCV/HIV coinfection (22 patients with HCV genotype 1 infection, 6 with genotype 2, and 7 with genotype 3) were treated with sofosbuvir and ledipasvir (for genotype 1 patients) or sofosbuvir and ribavirin (for genotypes 2 and 3). Sustained virological response (SVR) at 24 weeks after end of treatment and adverse events were assessed.

RESULTS

The overall SVR rate was 91.4% (32/35). One patient with genotype 1 infection discontinued treatment on day 2 due to severe headache, which subsided after the cessation of medication; all other patients completed their treatment without severe adverse events. Two patients who had a relapse of HCV were infected with a genotype 3 strain. We observed hyperbilirubinemia in a patient with genotype 3, who was under antiretroviral therapy including atazanavir. He completed the treatment and achieved SVR.

CONCLUSION

Direct-acting antiviral treatment for patients coinfected with HCV/HIV is as effective as in patients infected only with HCV. It was generally well tolerated, except in one patient who discontinued the treatment due to severe headache.

Evolutionary pathways to NS5A inhibitor resistance in genotype 1 hepatitis C virus

Direct-acting antivirals (DAAs) targeting NS5A are broadly effective against hepatitis C virus (HCV) infections, but sustained virological response rates are generally lower in patients infected with genotype (gt)-1a than gt-1b viruses. The explanation for this remains uncertain. Here, we adopted a highly accurate, ultra-deep primer ID sequencing approach to intensively study serial changes in the NS5A-coding region of HCV in gt-1a- and gt-1b-infected subjects receiving a short course of monotherapy with the NS5A inhibitor, elbasvir. Low or undetectable levels of viremia precluded on-treatment analysis in gt-1b-infected subjects, but variants with the resistance-associated substitution (RAS) Y93H in NS5A dominated rebounding virus populations following cessation of treatment. These variants persisted until the end of the study, two months later. In contrast, while Y93H emerged in multiple lineages and became dominant in subjects with gt-1a virus, these haplotypes rapidly decreased in frequency off therapy. Substitutions at Q30 and L31 emerged in distinctly independent lineages at later time points, ultimately coming to dominate the virus population off therapy. Consistent with this, cell culture studies with gt-1a and gt-1b reporter viruses and replicons demonstrated that Y93H confers a much greater loss of replicative fitness in gt-1a than gt-1b virus, and that L31M/V both compensates for the loss of fitness associated with Q30R (but not Y93H) and also boosts drug resistance. These observations show how differences in the impact of RASs on drug resistance and replicative fitness influence the evolution of gt-1a and gt-1b viruses during monotherapy with an antiviral targeting NS5A.

Prevalence and impact of baseline resistance-associated substitutions on the efficacy of ledipasvir/sofosbuvir or simeprevir/sofosbuvir against GT1 HCV infection

Baseline resistance-associated substitutions (RASs) have variable impacts in clinical trials but their prevalence and impact in real-world patients remains unclear. We performed baseline resistance testing using a commercial assay (10% cutoff) for 486 patients treated with LDV/SOF or SMV/SOF, with or without ribavirin, in the multi-center, observational HCV-TARGET cohort. Linkage of RASs was evaluated in selected samples using a novel quantitative single variant sequencing assay. Our results showed that the prevalence of NS3, NS5A, NS5B RASs was 45%, 13%, and 8%, respectively, and 10% of patients harbored RASs in 2 or more drug classes. Baseline LDV RASs in GT1a, TE, and cirrhosis LDV/SOF subgroup was associated with 2-4% lower SVR12 rates. SMV RASs was associated with lower SVR12 rates in GT1a, treatment-experienced, cirrhotics SMV/SOF subgroup. Pooled analysis of all patients with baseline RASs revealed that SVR12 was 100% (19/19) in patients treated for longer than 98 days but was 87% (81/93) in patients treated for shorter than 98 days. These results demonstrate that RASs prevalence and their impact in real world practice are in general agreement with registration trials, and suggest that longer treatment duration may overcome the negative impact of baseline RASs on SVR12 rates in clinical practice.

Effect of HIV/HCV Co-Infection on the Protease Evolution of HIV-1B: A Pilot Study in a Pediatric Population

This pilot study evaluates in pediatric patients the impact of HIV/HCV coinfection in the molecular evolution of the HIV-1 subtype B protease (HIV-1BPR). For this study, HIV-1B/HCV coinfected (15) and HIV-1B monoinfected (56) patients with available HIV-1B pol sequences were enrolled. Both groups of patients had comparable gender frequencies and average age, time of infection, antiretroviral treatment (ART) exposure and time under ART. Prevalence of drug resistance mutations (DRM), genetic diversity, number of synonymous (d_S) and non-synonymous (d_N) mutations per site and selection pressures (d_N - d_S) in the HIV-1BPR were estimated and compared between mono- and coinfected patients. Both HIV-1B populations presented similar genetic diversity (0.050 ± 0.02 vs. 0.045 ± 0.01) and d_S (0.074 ± 0.03 vs. 0.078 ± 0.04). In turn, in coinfected patients the HIV-1BPR had higher d_N (0.045 ± 0.01 vs. 0.024 ± 0.01) and d_N-d_S (-0.026 ± 0.02 vs. -0.048 ± 0.04) values, and less amino acid sites under purifying selection (4.2% vs. 42.1%) than in monoinfected patients. Accordingly, in co-infection with HCV, the HIV-1BPR sites 50, 53, 82, 84 and 88 - associated with resistance to PIs - were under neutral evolution, whereas these sites were under purifying selection in monoinfected patients. This pilot study suggests that HIV-1B may evolve differently in the presence than in the absence of HCV.

Consequence of HIV and HCV co-infection on host immune response, persistence and current treatment options

Hepatitis C virus (HCV) is a common opportunistic pathogen especially among Human immunodeficiency virus (HIV) infected patients. Due to incongruous studies, the pathological effect of HCV on HIV induced disease are still not fully understood. While some studies have showed no effect of HCV on HIV infection, others reported a defined role of HCV in aggravating the rates of AIDS-related illnesses and mortality. The explanation of such variances may be due to the host immune response, viral genotypes, sub-type and quasi-species distribution. The factors that complicate the management of HIV/HCV patients are: (1) reduced HCV antibody production, (2) drug interactions, (3) liver disease and (4) different epidemiologic characteristics. However, it is abundantly clear that the morbidity and mortality caused by HCV have increased since the introduction of highly active antiretroviral therapy (HAART) against HIV. In this review, the consequence of HIV/HCV co-infection on host immune response, viral replication, disease progression, mortality and morbidity, viral load, persistence and current treatment options have been discussed. Based on the clinical studies, it is necessary to evaluate the effect of HCV therapy on HIV progression and to provide a fully active HCV treatment for patients receiving HIV treatment. In conclusion, it is recommended to provide fully active HAART therapy in combination with a known HCV therapy.

Characterization of HCV NS3 Protease Variants in HCV/HIV-Coinfected Patients by Ultra-Deep Sequence Analysis: Relationship with Hepatic Fibrosis

BACKGROUND

Treatment of HCV/HIV coinfection is now largely based on utilization of direct acting agents. Pretreatment viral resistant-associated variants (RAVs) and host liver condition may affect the sustained virological response. In this study, we explored relative prevalence of protease resistance-associated mutations, the evolution of those RAVs after 12 weeks of pegylated interferon alfa exposure, and the role hepatic fibrosis might have on RAV display.

METHODS

Thirty nonresponder HCV/HIV-coinfected subjects were evaluated before and after 12 weeks of PegIFN treatment. Ultra-deep sequence analysis of NS3 RAVs was performed. Hepatic fibrosis was determined by sensitive computer-assisted histomorphometry determination.

RESULTS

At baseline, protease inhibitor RAVs were present in 73.3% of patients and expanded to 83.3% of patients after 12 weeks of PegIFN exposure. Q80K showed the highest prevalence before and after treatment at 46.7% and 56.7%, respectively. The presence of Q80K is positively correlated with percent collagen content of the liver tissue.

CONCLUSIONS

Key RAVs for HCV protease inhibitors are present in a major portion of the HCV/HIV-coinfected population before therapy. Some variants get selected after exposure. Correlation of Q80K with collagen content of the liver suggests that compartmentalization within the liver may contribute to persistence of mutations less fit than wildtype.

The hepatitis C virus nonstructural protein 3 Q80K polymorphism is frequently detected and transmitted among HIV-infected MSM in the Netherlands

OBJECTIVES

The Q80K polymorphism is a naturally occurring resistance-associated variant in the hepatitis C virus (HCV) nonstructural protein 3 (NS3) region and is likely transmissible between hosts. This study describes the Q80K origin and prevalence among HCV risk groups in the Netherlands and examines whether Q80K is linked to specific transmission networks.

DESIGN AND METHODS

Stored blood samples from HCV genotype 1a-infected patients were used for PCR and sequencing to reconstruct the NS3 maximum likelihood phylogeny. The most recent common ancestor was estimated with a coalescent-based model within a Bayesian statistical framework.

RESULTS

Study participants (n = 150) were either MSM (39%), people who inject drugs (17%), or patients with other (15%) or unknown/unreported (29%) risk behavior. Overall 45% was coinfected with HIV. Q80K was present in 36% (95% confidence interval 28-44%) of patients throughout the sample collection period (2000-2015) and was most prevalent in MSM (52%, 95% confidence interval 38-65%). Five MSM-specific transmission clusters were identified, of which three exclusively contained sequences with Q80K. The HCV-1a most recent common ancestor in the Netherlands was estimated in 1914 (95% higher posterior density 1879-1944) and Q80K originated in 1957 (95% higher posterior density 1942-1970) within HCV-1a clade I. All Q80K lineages could be traced back to this single origin.

CONCLUSION

Q80K is a highly stable and transmissible resistance-associated variant and was present in a large part of Dutch HIV-coinfected MSM. The introduction and expansion of Q80K variants in this key population suggest a founder effect, potentially jeopardizing future treatment with simeprevir.

NS3 Resistance-Associated Variants (RAVs) in Patients Infected with HCV Genotype 1a in Spain

Background

Resistance-associated variants have been related to treatment failure of hepatitis C virus (HCV) therapy with direct-acting antiviral drugs. The aim of our study was to analyze the prevalence of clinically relevant resistance-associated variants within NS3 in patients infected with HCV genotype 1a (GT1a) in Spain.

Methods

We performed a cross-sectional study on 2568 patients from 115 hospitals throughout Spain (2014–2015). The viral NS3 protease gene was amplified by nested polymerase chain reaction and sequenced by Sanger sequencing using an ABI PRISM 377 DNA sequencer. Additionally, clade information for genotype 1a was obtained by using the software geno2pheno (http://hcv.geno2pheno.org/).

Results

In total, 875 out of 2568 samples were from human immunodeficiency virus (HIV)/HCV-coinfected patients. Q80K was the main RAV found in our patients (11.1%) and the rest of the resistance-associated variants had a lower frequency, including S122G (6.23%), T54S (3.47%), V55A (2.61%), and V55I (2.15%), which were among the most frequent after Q80K. Overall, 286 samples had the Q80K polymorphism (11.1%) and 614 (23.9%) were GT1a clade I. HIV/HCV-coinfected patients had a higher frequency of Q80K and GT1a clade I than HCV-monoinfected patients (12.9% vs. 9.6% [p = 0.012] and 28.5% vs. 21.4% [p<0.001], respectively). Both the prevalence of Q80K and GT1a clade I were not uniform throughout the country (p<0.001), which ranged from 7.3%-22.2% and 15.7%-42.5%, respectively. The frequency of the Q80K polymorphism was far higher in patients infected with GT1a clade I than in patients infected with GT1a clade II (41.5% vs. 1.6%; p<0.001).

Conclusions

The prevalence of most resistance-associated variants in NS3 was low in patients infected with HCV GT1a in Spain, except for Q80K (11.1%), which was also notably higher in HIV/HCV-coinfected patients. The vast majority of Q80K polymorphisms were detected in GT1a clade I.

Natural Polymorphisms Conferring Resistance to HCV Protease and Polymerase Inhibitors in Treatment-Naïve HIV/HCV Co-Infected Patients in China

BACKGROUND

The advent of direct-acting agents (DAAs) has improved treatment of HCV in HIV co-infection, but may be limited by primary drug resistance. This study reports the prevalence of natural polymorphisms conferring resistance to NS3/4A protease inhibitors and NS5B polymerase inhibitors in treatment-naïve HIV/HCV co-infected individuals in China.

METHODS

Population based NS3/4A sequencing was completed for 778 treatment-naïve HIV/HCV co-infected patients from twelve provinces. NS3 sequences were amplified by nested PCR using in-house primers for genotypes 1-6. NS5B sequencing was completed for genotyping in 350 sequences. Resistance-associated variants (RAVs) were identified in positions associated with HCV resistance.

RESULTS

Overall, 72.8% (566/778) of all HCV sequences had at least one RAV associated with HCV NS3/4A protease inhibitor resistance. Variants were found in 3.6% (7/193) of genotype 1, 100% (23/23) of genotype 2, 100% (237/237) of genotype 3 and 92% (299/325) of genotype 6 sequences. The Q80K variant was present in 98.4% of genotype 6a sequences. High-level RAVs were rare, occurring in only 0.8% of patients. 93% (64/69) patients with genotype 1b also carried the C316N variant associated with NS5B low-level resistance.

CONCLUSIONS

The low frequency of high-level RAVs associated with primary HCV DAA resistance among all genotypes in HIV/HCV co-infected patients is encouraging. Further phenotypic studies and clinical research are needed.

Antiviral drug resistance as an adaptive process

Antiviral drug resistance is a matter of great clinical importance that, historically, has been investigated mostly from a virological perspective. Although the proximate mechanisms of resistance can be readily uncovered using these methods, larger evolutionary trends often remain elusive. Recent interest by population geneticists in studies of antiviral resistance has spurred new metrics for evaluating mutation and recombination rates, demographic histories of transmission and compartmentalization, and selective forces incurred during viral adaptation to antiviral drug treatment. We present up-to-date summaries on antiviral resistance for a range of drugs and viral types, and review recent advances for studying their evolutionary histories. We conclude that information imparted by demographic and selective histories, as revealed through population genomic inference, is integral to assessing the evolution of antiviral resistance as it pertains to human health.

Naturally occurring hepatitis C virus protease inhibitors resistance-associated mutations among chronic hepatitis C genotype 1b patients with or without HIV co-infection

The aim of this study was to measure the frequency of natural mutations in hepatitis C virus (HCV) mono-infected and HIV/HCV co-infected protease inhibitor (PI)-naive patients.

METHODS

Population sequence of the non-structural (NS)3 protease gene was evaluated in 90 HCV mono-infected and 96 HIV/HCV co-infected PI treatment-naive patients. The natural prevalence of PI resistance mutations in both groups was compared.

RESULTS

Complete HCV genotype 1b NS3 sequence information was obtained for 152 (81.72%) samples. Seven sequences (8.33%) of the 84 HCV mono-infected patients and 21 sequences (30.88%) of the 68 HIV/HCV co-infected patients showed amino acid substitutions associated with HCV PI resistance. There was a significant difference in the natural prevalence of PI resistance mutations between these two groups (P = 0.000). The mutations T54S, R117H and N174F were observed in 1.19%, 5.95% and 1.19% of HCV mono-infected patients. The mutations F43S, T54S, Q80K/R, R155K, A156G/V, D168A/E/G and V170A were found in 1.47%, 4.41%, 1.47%/1.47%, 2.94%, 23.53%/1.47%, 1.47%/1.47%/1.47% and 1.47% of HIV/HCV co-infected patients, respectively. In addition, the combination mutations in the NS3 region were detected only in HIV/HCV genotype 1b co-infected patients.

CONCLUSION

Naturally occurring HCV PI resistance mutations existed in HCV mono-infected and HIV/HCV co-infected genotype 1b PI-naive patients. HIV co-infection was associated with a greater frequency of PI resistance mutations. The impact of HIV infection on baseline HCV PI resistance mutations and treatment outcome in chronic hepatitis C (CHC) patients should be further analyzed.

CD4+ T-cell recovery with suppressive ART-induced rapid sequence evolution in hepatitis C virus envelope but not NS3

OBJECTIVES

CD4 T-cell depletion from HIV infection leads to a global decline in anti-hepatitis C virus (HCV) envelope neutralizing antibody (nAb) response, which may play a role in accelerating liver fibrosis. An increase in anti-HCV nAb titers has been reported during antiretroviral therapy (ART) but its impact on HCV remains poorly understood. The objective of this study is to determine the effects of ART on long-term HCV evolution.

DESIGN AND METHODS

We examined HCV quasispecies structure and long-term evolution in HIV/HCV coinfected patients with ART-induced CD4 T-cell recovery, and compared with patients with CD4 T-cell depletion from delayed ART. We applied a single-variant sequencing (SVS) method to construct authentic viral quasispecies and compared sequence evolution in HCV envelope, the primary target for humoral immune responses, and NS3, a target for cellular immunity, between the two cohorts.

RESULTS

The SVS method corrected biases known to skew the proportions of viral variants, revealing authentic HCV quasispeices structures. We observed higher rates of HCV envelope sequence evolution in patients with ART-induced CD4 T-cell recovery, compared with patients with CD4 T-cell depletion from delayed ART (P = 0.03). Evolutionary rates for NS3 were considerably lower than the rates for envelope (P < 0.01), with no significant difference observed between the two groups.

CONCLUSION

ART-induced CD4 T-cell recovery results in rapid sequence evolution in HCV envelope, but not in NS3. These results suggest that suppressive ART disproportionally enhances HCV-specific humoral responses more than cellular responses, resulting in rapid sequence evolution in HCV envelope but not NS3.

Evolution and function of the HCV NS3 protease in patients with acute hepatitis C and HIV coinfection

Little is known about the importance of the hepatitis C virus (HCV) NS3 protease in acute hepatitis C. In this prospective study, 82 consecutive patients with acute hepatitis C and human immunodeficiency virus (HIV) coinfection were enrolled. Individuals were infected with highly related HCV strains and the baseline NS3 quasispecies diversity and complexity was higher compared to a chronic hepatitis C control group (P<0.0001). Both parameters were comparable in patients with spontaneous clearance (n=6) versus treatment-induced SVR (n=5) or development of chronic hepatitis C (n=9). Longitudinal NS3 quasispecies kinetics showed a trend to a decreasing diversity and complexity (P<0.05) within 4 weeks in patients with spontaneous clearance compared to the other groups. The innate immune signalling protein CARDIF was cleaved to a similar extent independent of the outcome. Together with a more pronounced viral load decline (P<0.05), an early decreasing NS3 quasispecies evolution indicates spontaneous clearance of acute hepatitis C.

Similarities between Human Immunodeficiency Virus Type 1 and Hepatitis C Virus Genetic and Phenotypic Protease Quasispecies Diversity

Human immunodeficiency virus type 1 (HIV-1) and hepatitis C virus (HCV) are two highly variable RNA viruses that cause chronic infections in humans. Although HCV likely preceded the AIDS epidemic by some decades, the global spread of both viruses is a relatively recent event. Nevertheless, HCV global diversity is higher than that of HIV-1. To identify differences in mutant diversity, we compared the HIV-1 protease and HCV NS3 protease quasispecies. Three protease gene quasispecies samples per virus, isolated from a total of six infected patients, were genetically and phenotypically analyzed at high resolution (HIV-1, 308 individual clones; HCV, 299 clones). Single-nucleotide variant frequency did not differ between quasispecies from the two viruses (HIV-1, 2.4 × 10(-3) ± 0.4 × 10(-3); HCV, 2.1 × 10(-3) ± 0.5 × 10(-3)) (P = 0.1680). The proportion of synonymous substitutions to potential synonymous sites was similar (3.667 ± 0.6667 and 2.183 ± 0.9048, respectively) (P = 0.2573), and Shannon's entropy values did not differ between HIV-1 and HCV (0.84 ± 0.02 and 0.83 ± 0.12, respectively) (P = 0.9408). Of note, 65% (HIV-1) and 67% (HCV) of the analyzed enzymes displayed detectable protease activity, suggesting that both proteases have a similar mutational robustness. In both viruses, there was a rugged protease enzymatic activity landscape characterized by a sharp peak, representing the master sequence, surrounded by a collection of diverse variants present at lower frequencies. These results indicate that nucleotide quasispecies diversification during chronic infection is not responsible for the higher worldwide genetic diversity observed in HCV.

IMPORTANCE

HCV global diversity is higher than that of HIV-1. We asked whether HCV genetic diversification during infection is responsible for the higher worldwide genetic diversity observed in HCV. To this end, we analyzed and compared the genotype and enzymatic activities of HIV-1 and HCV protease quasispecies existing in infected individuals. Our results indicate that HIV-1 and HCV protease quasispecies have very similar genetic diversity and comparable rugged enzymatic activity landscapes. Therapy for HCV has expanded, with new therapeutic agents such as the direct-acting antivirals (DAAs). DAAs, which target HCV NS3 protease and other virus proteins, have improved cure rates. However, major questions remain to be elucidated regarding the virologic correlates of HCV eradication. The findings shown here may help our understanding of the different therapeutic responses observed during chronic HCV infection.

Development and Validation of Two Screening Assays for the Hepatitis C Virus NS3 Q80K Polymorphism Associated with Reduced Response to Combination Treatment Regimens Containing Simeprevir

Persons with hepatitis C virus (HCV) genotype 1a (GT1a) infections harboring a baseline Q80K polymorphism in nonstructural protein 3 (NS3) have a reduced virologic response to simeprevir in combination with pegylated interferon-alfa and ribavirin. We aimed to develop, validate, and freely disseminate an NS3 clinical sequencing assay to detect the Q80K polymorphism and potentially other HCV NS3 drug resistance mutations. HCV RNA was extracted from frozen plasma using a NucliSENS easyMAG automated nucleic acid extractor, amplified by nested reverse transcription-PCR, and sequenced using Sanger and/or next-generation (MiSeq) methods. Sanger chromatograms were analyzed using in-house software (RECall), and nucleotide mixtures were called automatically. MiSeq reads were iteratively mapped to the H77 reference genome, and consensus NS3 sequences were generated with nucleotides present at >20% called as mixtures. The accuracy, precision, and sensitivity for detecting the Q80K polymorphism were assessed in 70 samples previously sequenced by an external laboratory. A comparison of the sequences generated by the Sanger and MiSeq methods with those determined by an external lab revealed >98.5% nucleotide sequence concordance and zero discordant calls of the Q80K polymorphism. The results were both highly repeatable and reproducible (>99.7% nucleotide concordance and 100% Q80K concordance). The limits of detection (>2 and ∼5 log10 IU/ml for the Sanger and MiSeq assays, respectively) are sufficiently low to allow genotyping in nearly all chronically infected treatment-naive persons. No systematic bias in the under- or overamplification of minority variants was observed. Coinfection with other viruses (e.g., HIV and hepatitis B virus [HBV]) did not affect the assay results. The two independent HCV NS3 sequencing assays with the automated analysis procedures described here are useful tools to screen for the Q80K polymorphism and other HCV protease inhibitor drug resistance mutations.

Previous failure of interferon-based therapy does not alter the frequency of HCV NS3 protease or NS5B polymerase inhibitor resistance-associated variants: longitudinal analysis in HCV/HIV co-infected patients

Since 2011, treatment of chronic hepatitis C virus (HCV) includes direct-acting antivirals (DAAs) in addition to pegylated interferon-α (peg-IFN) and ribavirin (RBV). IFN-based treatment induces strong cytotoxic T-lymphocyte activity directed to the protease- and polymerase-derived epitopes. This enhanced immunological pressure could favour the emergence of viral epitope variants able to evade immune surveillance and, when resistance-associated variants (RAVs) are implicated, could also be co-selected as a hitchhiking effect. This study analysed the dynamics of the frequency of protease and polymerase inhibitor RAVs that could affect future HCV treatment in human immunodeficiency virus (HIV) co-infected patients on stable antiretroviral therapy with previous IFN-based treatment failure. HCV genotype 1a RNA was extracted from plasma samples of 18 patients prior to and during (24h and 4, 12, 24 and 48 weeks) therapy with peg-IFN+RBV. Next-generation sequencing was performed on HCV-RNA populations using NS3 and NS5B PCR-amplified coding regions. Two measures of genetic diversity were used to compare virus populations: average pairwise nucleotide diversity (π) and Tajima's D statistic. Several protease and polymerase RAVs were detected in all subjects at very low frequencies (<5%), and in most cases their presence was not constant during follow-up. Only samples from two patients for each region exhibited Q80R/K/L and A421V as highly predominant variants. No significant differences were observed among sampling times for either π or D values. In conclusion, previous therapy and failure of peg-IFN+RBV were not associated with an increase in DAA-targeting NS3 or NS5B RAVs that naturally exist in HIV co-infected subjects.

Primer ID Validates Template Sampling Depth and Greatly Reduces the Error Rate of Next-Generation Sequencing of HIV-1 Genomic RNA Populations

Validating the sampling depth and reducing sequencing errors are critical for studies of viral populations using next-generation sequencing (NGS). We previously described the use of Primer ID to tag each viral RNA template with a block of degenerate nucleotides in the cDNA primer. We now show that low-abundance Primer IDs (offspring Primer IDs) are generated due to PCR/sequencing errors. These artifactual Primer IDs can be removed using a cutoff model for the number of reads required to make a template consensus sequence. We have modeled the fraction of sequences lost due to Primer ID resampling. For a typical sequencing run, less than 10% of the raw reads are lost to offspring Primer ID filtering and resampling. The remaining raw reads are used to correct for PCR resampling and sequencing errors. We also demonstrate that Primer ID reveals bias intrinsic to PCR, especially at low template input or utilization. cDNA synthesis and PCR convert ca. 20% of RNA templates into recoverable sequences, and 30-fold sequence coverage recovers most of these template sequences. We have directly measured the residual error rate to be around 1 in 10,000 nucleotides. We use this error rate and the Poisson distribution to define the cutoff to identify preexisting drug resistance mutations at low abundance in an HIV-infected subject. Collectively, these studies show that >90% of the raw sequence reads can be used to validate template sampling depth and to dramatically reduce the error rate in assessing a genetically diverse viral population using NGS.

IMPORTANCE

Although next-generation sequencing (NGS) has revolutionized sequencing strategies, it suffers from serious limitations in defining sequence heterogeneity in a genetically diverse population, such as HIV-1 due to PCR resampling and PCR/sequencing errors. The Primer ID approach reveals the true sampling depth and greatly reduces errors. Knowing the sampling depth allows the construction of a model of how to maximize the recovery of sequences from input templates and to reduce resampling of the Primer ID so that appropriate multiplexing can be included in the experimental design. With the defined sampling depth and measured error rate, we are able to assign cutoffs for the accurate detection of minority variants in viral populations. This approach allows the power of NGS to be realized without having to guess about sampling depth or to ignore the problem of PCR resampling, while also being able to correct most of the errors in the data set.

A molecular phylogenetics-based approach for identifying recent hepatitis C virus transmission events

Improved surveillance methods are needed to better understand the current hepatitis C virus (HCV) disease burden and to monitor the impact of prevention and treatment interventions on HCV transmission dynamics. Sanger sequencing (HCV NS5B, HVR1 and Core-E1-HVR1) and phylogenetics were applied to samples from individuals diagnosed with HCV in British Columbia, Canada in 2011. This included individuals with two or three sequential samples collected <1 year apart. Patristic distances between sequential samples were used to set cutoffs to identify recent transmission clusters. Factors associated with transmission clustering were analyzed using logistic regression. From 618 individuals, 646 sequences were obtained. Depending on the cutoff used, 63 (10%) to 92 (15%) unique individuals were identified within transmission clusters of predicted recent origin. Clustered individuals were more likely to be <40 years old (Adjusted Odds Ratio (AOR) 2.12, 95% CI 1.21-3.73), infected with genotype 1a (AOR 6.60, 95% CI 1.98-41.0), and to be seroconverters with estimated infection duration of <1 year (AOR 3.13, 95% CI 1.29-7.36) or >1 year (AOR 2.19, 95% CI 1.22-3.97).

CONCLUSION

Systematic application of molecular phylogenetics may be used to enhance traditional surveillance methods through identification of recent transmission clusters.

Prevalence of HCV NS3 pre-treatment resistance associated amino acid variants within a Scottish cohort

BACKGROUND

Protease inhibitors (PI) including boceprevir, telaprevir and simeprevir have revolutionised HCV genotype 1 treatment since their introduction. A number of pre-treatment resistance associated amino acid variants (RAVs) and polymorphisms have been associated with reduced response to treatment.

OBJECTIVES

We measured the prevalence of RAVs/polymorphisms in a PI treatment-naïve HCV genotype 1 Scottish cohort using Sanger sequencing.

STUDY DESIGN

Chronically infected, treatment-naïve, HCV genotype 1 patients (n = 146) attending NHS Greater Glasgow and Clyde clinics were investigated for RAVs/polymorphisms to the PIs boceprevir, telaprevir and simeprevir. The NS3/4A region was amplified by nested polymerase chain reaction. The 1.4 kb amplified product was sequenced using an ABI 3710XL DNA sequencer. Sequence analysis was performed using web-based ReCall (beta 2.10). Amino acid positions 36, 41, 43, 54, 55, 80, 109, 122, 155, 156, 168 and 170 were analysed for RAVs/polymorphisms.

RESULTS

Overall, 23.29% (34/146) of patients had an RAV or polymorphism detected. Overall, 13.69% (20/146) of patients had HCV virus that contained the Q8 K polymorphism. Other RAVs detected were: V36 M 0.70% (1/146), V36L 0.70% (1/146), T54S 6.85% (10/146), V55A 3.42% (5/146) and V/I170A 0.68% (1/146). Four patients had dual combinations of mutations (T54S+V36L; T54S+V55A and 2 patients with T54S+Q80K).

CONCLUSIONS

Q80K was the most prevalent baseline polymorphism detected in the Scottish cohort. Simeprevir treatment is not recommended in patients infected with the Q80K genotype 1a variant. This highlights the need for baseline sequencing prior to administration of this drug in this population.

Deconvoluting the composition of low-frequency hepatitis C viral quasispecies: comparison of genotypes and NS3 resistance-associated variants between HCV/HIV coinfected hemophiliacs and HCV monoinfected patients in Japan

Pre-existing low-frequency resistance-associated variants (RAVs) may jeopardize successful sustained virological responses (SVR) to HCV treatment with direct-acting antivirals (DAAs). However, the potential impact of low-frequency (∼0.1%) mutations, concatenated mutations (haplotypes), and their association with genotypes (Gts) on the treatment outcome has not yet been elucidated, most probably owing to the difficulty in detecting pre-existing minor haplotypes with sufficient length and accuracy. Herein, we characterize a methodological framework based on Illumina MiSeq next-generation sequencing (NGS) coupled with bioinformatics of quasispecies reconstruction (QSR) to realize highly accurate variant calling and genotype-haplotype detection. The core-to-NS3 protease coding sequences in 10 HCV monoinfected patients, 5 of whom had a history of blood transfusion, and 11 HCV/HIV coinfected patients with hemophilia, were studied. Simulation experiments showed that, for minor variants constituting more than 1%, our framework achieved a positive predictive value (PPV) of 100% and sensitivities of 91.7–100% for genotyping and 80.6% for RAV screening. Genotyping analysis indicated the prevalence of dominant Gt1a infection in coinfected patients (6/11 vs 0/10, p = 0.01). For clinical samples, minor genotype overlapping infection was prevalent in HCV/HIV coinfected hemophiliacs (10/11) and patients who experienced whole-blood transfusion (4/5) but none in patients without exposure to blood (0/5). As for RAV screening, the Q80K/R and S122K/R variants were particularly prevalent among minor RAVs observed, detected in 12/21 and 6/21 cases, respectively. Q80K was detected only in coinfected patients, whereas Q80R was predominantly detected in monoinfected patients (1/11 vs 7/10, p < 0.01). Multivariate interdependence analysis revealed the previously unrecognized prevalence of Gt1b-Q80K, in HCV/HIV coinfected hemophiliacs [Odds ratio = 13.4 (3.48–51.9), p < 0.01]. Our study revealed the distinct characteristics of viral quasispecies between the subgroups specified above and the feasibility of NGS and QSR-based genetic deconvolution of pre-existing minor Gts, RAVs, and their interrelationships.