Genotype- and Subtype-Independent Full-Genome Sequencing Assay for Hepatitis C Virus

Determination of hepatitis C virus subtype prevalent in Sindh, Pakistan: a phylogenetic analysis

Hepatitis is a major public health issue, affecting 10-17 million people worldwide, with its prevalence continuously increasing. The Hepatitis C virus (HCV) is responsible for liver related diseases, which include liver cirrhosis, hepatocellular carcinoma, and chronic hepatitis. Pakistan is experiencing a serious rise in HCV cases. This study aimed to assess the prevalence and distribution of HCV genotypes in Sindh, Pakistan. Serum samples from HCV-positive patients were collected from various local hospitals in Sindh. These samples were first screened for HCV antibodies using ELISA. Samples that tested positive for HCV RNA underwent further genotyping through sequencing using the standard Sanger method. The genotypes were identified by comparing the sequences with those available in the National Center for Biotechnology Information (NCBI) database, and a phylogenetic tree was constructed. The phylogenetic analysis showed that all isolates in this study were clustered with genotypes 3a and 3b, except for one sequence that was clustered with genotype 1a. No isolates were found to be clustered with reference genomes of genotypes 2, 4, 5, 6, and 7 suggesting that genotype 3a is endemic in this region. The analyzed sequences demonstrated a 98% similarity with reference and isolated sequences. In summary, sequencing of the HCV 5' UTR essential for identifying the predominant genotype of HCV RNA in the Sindh region Further research on the distribution of HCV genotypes in other regions of Pakistan could aid in improving screening processes, identifying more effective treatment options, and developing suitable prevention strategies.

Sequence characterization of extracellular HBV RNA in patient plasma

Antiviral nucleos(t)ide analogue therapies inhibit HBV replication and suppress the HBV DNA levels in patients with chronic HBV infection. Since HBV RNAs are expressed from cccDNA or HBV integrated sequences, independently of viral genome replication, levels of HBV RNAs in plasma may remain high following treatment with nucleos(t)ide analogue. Thus, HBV RNAs have been proposed to be used as a viral biomarker for treatment outcome and disease progression. Recent investigations of plasma HBV RNAs described the presence of full length as well as subgenomic forms of RNA. To support the usage of plasma HBV RNAs as a viral biomarker, further understanding of HBV RNA composition in clinical samples is needed. Here, sequence of extracellular HBV RNAs was characterized in plasma samples of patients with chronic HBV infection using two independent RNA amplification methods that do not use HBV-specific primers for amplification: total RNA (NuGEN RNAseq) and mRNA (TruSeq RNAseq). Sequencing coverage was obtained across the full length of HBV genome for both methods, confirming the presence of full-length HBV RNA in plasma. The sequence of HBV RNA was nearly identical to plasma HBV DNA sequence in each sample with only 0-14 (median 4) mismatches over 3 kb. Thus, sequence of HBV RNA plasma reflects the intrahepatic viral reservoir and can be used for monitoring of sequence variants such as resistance in clinical trials. Additionally, RNA splice forms, different polyA tails start positions and presence of HBV-human chimeric transcript were identified.

Safety and Efficacy of Limited Laboratory Monitoring for Hepatitis C Treatment: A Blinded Clinical Trial in Rwanda

Direct-acting antivirals for hepatitis C virus (HCV) are highly effective and well-tolerated. However, only a small percentage of HCV-infected individuals globally have received therapy. Reducing the complexity of monitoring during HCV therapy, if shown to be safe, could facilitate greater access to HCV services, particularly in resource-limited settings such as sub-Saharan Africa. We enrolled a total of 300 patients who were chronically infected with genotype 4 HCV in Rwanda and treated them with fixed-dose ledispasvir/sofosbuvir for 12 weeks. For 60 consecutive participants enrolled, we blinded the study clinician to on-treatment laboratory results. We compared the efficacy, safety, and tolerability in those with blinded laboratory results to those with standard laboratory monitoring. Baseline characteristics among those with blinded laboratory values were comparable to those with standard monitoring. Among both groups, the median age was 63 years, and the median HCV viral load was 5.9 log (versus 64 years and 6.0 log, respectively). Sustained virologic response rates at 12 weeks after treatment completion were similar in those with blinded laboratories (87%) compared to those with standard laboratory monitoring (87%). There was no increase in adverse events in those with blinded laboratory results, and no participants discontinued the study medication because of an adverse event. Conclusion: On-treatment laboratory monitoring did not improve patient outcomes in those treated with ledispasvir/sofosbuvir. Eliminating this monitoring in treatment programs in resource-limited settings may facilitate and accelerate scale-up of HCV therapy.

Phylogenetic surveillance of travel-related Zika virus infections through whole-genome sequencing methods

In 2018, the World Health Organization identified the Zika virus (ZIKV) as a pathogen that should be prioritized for public health research due to its epidemic potential. In this study, whole-genome sequencing (WGS) of travel-acquired ZIKV infections was used to examine the limitations of phylogenetic analysis. WGS and phylogenetic analysis were performed to investigate geographic clustering of samples from five Canadians with travel-acquired ZIKV infections and to assess the limitations of phylogenetic analysis of ZIKV sequences using a phylogenetic cluster approach. Genomic variability of ZIKV samples was assessed and for context, compared with hepatitis C virus (HCV) samples. Phylogenetic analysis confirmed the suspected region of ZIKV infection for one of five samples and one sample failed to cluster with sequences from its suspected country of infection. Travel-acquired ZIKV samples depicted low genomic variability relative to HCV samples. A floating patristic distance threshold classified all pre-2000 ZIKV sequences into separate clusters, while only Cambodian, Peruvian, Malaysian, and South Korean sequences were similarly classifiable. While phylogenetic analysis of ZIKV data can identify the broad geographical region of ZIKV infection, ZIKV’s low genomic variability is likely to limit precise interpretations of phylogenetic analysis of the origins of travel-related cases.

2'-C-methylated nucleotides terminate virus RNA synthesis by preventing active site closure of the viral RNA-dependent RNA polymerase

The 2'-C-methyl ribonucleosides are nucleoside analogs representing an important class of antiviral agents, especially against positive-strand RNA viruses. Their value is highlighted by the highly successful anti-hepatitis C drug sofosbuvir. When appropriately phosphorylated, these nucleotides are successfully incorporated into RNA by the virally encoded RNA-dependent RNA polymerase (RdRp). This activity prevents further RNA extension, but the mechanism is poorly characterized. Previously, we had identified NMR signatures characteristic of formation of RdRp-RNA binary and RdRp-RNA-NTP ternary complexes for the poliovirus RdRp, including an open-to-closed conformational change necessary to prepare the active site for catalysis of phosphoryl transfer. Here we used these observations as a framework for interpreting the effects of 2'-C-methyl adenosine analogs on RNA chain extension in solution-state NMR spectroscopy experiments, enabling us to gain additional mechanistic insights into 2'-C-methyl ribonucleoside-mediated RNA chain termination. Contrary to what has been proposed previously, poliovirus RdRp that was bound to RNA with an incorporated 2'-C-methyl nucleotide could still bind to the next incoming NTP. Our results also indicated that incorporation of the 2'-C-methyl nucleotide does not disrupt RdRp-RNA interactions and does not prevent translocation. Instead, incorporation of the 2'-C-methyl nucleotide blocked closure of the RdRp active site upon binding of the next correct incoming NTP, which prevented further nucleotide addition. We propose that other nucleotide analogs that act as nonobligate chain terminators may operate through a similar mechanism.

Hepatitis C virus genotype 1 and 2 recombinant genomes and the phylogeographic history of the 2k/1b lineage

Recombination is an important driver of genetic diversity, though it is relatively uncommon in hepatitis C virus (HCV). Recent investigation of sequence data acquired from HCV clinical trials produced twenty-one full-genome recombinant viruses belonging to three putative inter-subtype forms 2b/1a, 2b/1b, and 2k/1b. The 2k/1b chimera is the only known HCV circulating recombinant form (CRF), provoking interest in its genetic structure and origin. Discovered in Russia in 1999, 2k/1b cases have since been detected throughout the former Soviet Union, Western Europe, and North America. Although 2k/1b prevalence is highest in the Caucasus mountain region (i.e., Armenia, Azerbaijan, and Georgia), the origin and migration patterns of CRF 2k/1b have remained obscure due to a paucity of available sequences. We assembled an alignment which spans the entire coding region of the HCV genome containing all available 2k/1b sequences (>500 nucleotides; n = 109) sampled in ninteen countries from public databases (102 individuals), additional newly sequenced genomic regions (from 48 of these 102 individuals), unpublished isolates with newly sequenced regions (5 additional individuals), and novel complete genomes (2 additional individuals) generated in this study. Analysis of this expanded dataset reconfirmed the monophyletic origin of 2k/1b with a recombination breakpoint at position 3,187 (95% confidence interval: 3,172–3,202; HCV GT1a reference strain H77). Phylogeography is a valuable tool used to reveal viral migration dynamics. Inference of the timed history of spread in a Bayesian framework identified Russia as the ancestral source of the CRF 2k/1b clade. Further, we found evidence for migration routes leading out of Russia to other former Soviet Republics or countries under the Soviet sphere of influence. These findings suggest an interplay between geopolitics and the historical spread of CRF 2k/1b.

Next-generation sequencing for the clinical management of hepatitis C virus infections: does one test fits all purposes?

While the prospect of viral cure is higher than ever for individuals infected with the hepatitis C virus (HCV) due to ground-breaking progress in antiviral treatment, success rates are still negatively influenced by HCV's high genetic variability. This genetic diversity is represented in the circulation of various genotypes and subtypes, mixed infections, recombinant forms and the presence of numerous drug resistant variants among infected individuals. Common misclassifications by commercial genotyping assays in combination with the limitations of currently used targeted population sequencing approaches have encouraged researchers to exploit alternative methods for the clinical management of HCV infections. Next-generation sequencing (NGS), a revolutionary and powerful tool with a variety of applications in clinical virology, can characterize viral diversity and depict viral dynamics in an ultra-wide and ultra-deep manner. The level of detail it provides makes it the method of choice for the diagnosis and clinical assessment of HCV infections. The sequence library provided by NGS is of a higher magnitude and sensitivity than data generated by conventional methods. Therefore, these technologies are helpful to guide clinical practice and at the same time highly valuable for epidemiological studies. The decreasing costs of NGS to determine genotypes, mixed infections, recombinant strains and drug resistant variants will soon make it feasible to employ NGS in clinical laboratories, to assist in the daily care of patients with HCV.

Identification of 19 Novel Hepatitis C Virus Subtypes-Further Expanding HCV Classification

Background

Hepatitis C virus (HCV) is currently classified into 8 genotypes and 86 subtypes. The objective of this study was to characterize novel HCV subtypes and to investigate the impact of subtypes on treatment outcome.

Methods

Full-genome sequencing was performed on HCV plasma samples with <85% sequence homology of NS3, NS5A, and/or NS5B to HCV genotype (GT) 1–8 reference strains.

Results

A total of 14 653 patients with GT1–6 HCV infection were enrolled in clinical studies of sofosbuvir-based regimens. For the majority of the patients, a specific subtype could be assigned based on a close genetic relationship to previously described subtypes. However, for 19 patients, novel subtypes were identified with <85% homology compared with previously described subtypes. These novel subtypes had the following genotypes: 9 in GT2, 5 in GT4, 2 in GT6, and 1 each in GT1, GT3, and GT5. Despite the presence of polymorphisms at resistance-associated substitution positions, 18 of the 19 patients treated with sofosbuvir-containing therapy achieved SVR12.

Conclusions

Nineteen novel HCV subtypes were identified, suggesting an even greater genetic diversity of HCV subtypes than previously recognized.

Treatment of chronic hepatitis C virus infection in Rwanda with ledipasvir-sofosbuvir (SHARED): a single-arm trial

BACKGROUND

Limited treatment data are available for hepatitis C virus (HCV) in sub-Saharan Africa, especially for genotype 4. Our objective was to establish the safety and efficacy of ledipasvir-sofosbuvir for chronic HCV genotype 1 or 4 infection in adults in Rwanda.

METHODS

We did a single-arm trial to evaluate the safety and efficacy of ledipasvir-sofosbuvir in Rwandan adults with chronic HCV infection at a single study site (Rwanda Military Hospital, Kigali, Rwanda). We enrolled individuals aged 18 years or older with HCV genotype 1 or 4 infection and a plasma HCV RNA concentration of more than 1000 IU/mL at screening. All participants were given ledipasvir (90 mg) and sofosbuvir (400 mg) in a single combination tablet once daily for 12 weeks. We established HCV genotype using an Abbott platform, and HCV subtype with PCR amplification. The primary endpoint was the proportion of participants with a sustained virological response 12 weeks after therapy (SVR12). All patients enrolled in the study were included in the primary endpoint analyses. This study is registered with ClinicalTrials.gov, number NCT02964091.

FINDINGS

300 participants were enrolled between Feb 6, 2017, and Sept 18, 2017, and the follow-up period was completed on March 1, 2018. On genotyping, 248 (83%) participants were reported as having genotype 4, four (1%) genotype 1, and 48 (16%) both genotype 1 and genotype 4. Subsequent viral sequencing showed all participants actually had genotype 4 infection with subtype 4k (134 [45%]), subtype 4r (48 [16%]), subtype 4q (42 [14%]), and subtype 4v (24 [8%]) predominating. Overall, 261 (87%, 95% CI 83-91) participants achieved SVR12. In participants with genotype 4r, SVR12 was observed in 27 (56%, 95% CI 41-71) participants versus 234 (93%, 90-96) individuals with other subtypes. There were no drug-related treatment discontinuations due to ledipasvir-sofosbuvir. The most common adverse events were hypertension (97 [32%]), headache (78 [26%]), dizziness (61 [20%]), and fatigue (56 [19%]). There were six serious adverse events; none were assessed to be due to the study drug. 296 participants had data for pill counts at week 4 and 8; 271 (92%) had 100% adherence and only one (<1%) had an adherence of less than 90%.

INTERPRETATION

This is the first large-scale prospective study reporting direct-acting antiviral outcomes in sub-Saharan Africa. The high adherence and treatment success without intensive support measures or highly specialised clinical providers, and lack of treatment discontinuations due to adverse events support the feasibility of HCV treatment decentralisation and scale-up in sub-Saharan Africa. Genotype 4r is uniquely expressed in this region and associated with high rates of treatment failure, suggesting a need for rigorous test-of-cure in clinical practice and consideration of the use of newer pangenotypic direct-acting antiviral regimens in this region.

FUNDING

Gilead Sciences.

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.

Pipeline for specific subtype amplification and drug resistance detection in hepatitis C virus

BACKGROUND

Despite the high sustained virological response rates achieved with current directly-acting antiviral agents (DAAs) against hepatitis C virus (HCV), around 5-10% of treated patients do not respond to current antiviral therapies, and basal resistance to DAAs is increasingly detected among treatment-naïve infected individuals. Identification of amino acid substitutions (including those in minority variants) associated with treatment failure requires analytical designs that take into account the high diversification of HCV in more than 86 subtypes according to the ICTV website (June 2017).

METHODS

The methodology has involved five sequential steps: (i) to design 280 oligonucleotide primers (some including a maximum of three degenerate positions), and of which 120 were tested to amplify NS3, NS5A-, and NS5B-coding regions in a subtype-specific manner, (ii) to define a reference sequence for each subtype, (iii) to perform experimental controls to define a cut-off value for detection of minority amino acids, (iv) to establish bioinformatics' tools to quantify amino acid replacements, and (v) to validate the procedure with patient samples.

RESULTS

A robust ultra-deep sequencing procedure to analyze HCV circulating in serum samples from patients infected with virus that belongs to the ten most prevalent subtypes worldwide: 1a, 1b, 2a, 2b, 2c, 2j, 3a, 4d, 4e, 4f has been developed. Oligonucleotide primers are subtype-specific. A cut-off value of 1% mutant frequency has been established for individual mutations and haplotypes.

CONCLUSION

The methodological pipeline described here is adequate to characterize in-depth mutant spectra of HCV populations, and it provides a tool to understand HCV diversification and treatment failures. The pipeline can be periodically extended in the event of HCV diversification into new genotypes or subtypes, and provides a framework applicable to other RNA viral pathogens, with potential to couple detection of drug-resistant mutations with treatment planning.

A near full-length open reading frame next generation sequencing assay for genotyping and identification of resistance-associated variants in hepatitis C virus

BACKGROUND

The current treatment options for hepatitis C virus (HCV), based on direct acting antivirals (DAA), are dependent on virus genotype and previous treatment experience. Treatment failures have been associated with detection of resistance-associated substitutions (RASs) in the DAA targets of HCV, the NS3, NS5A and NS5 B proteins.

OBJECTIVE

To develop a next generation sequencing based method that provides genotype and detection of HCV NS3, NS5A, and NS5 B RASs without prior knowledge of sample genotype.

STUDY DESIGN

In total, 101 residual plasma samples from patients with HCV covering 10 different viral subtypes across 4 genotypes with viral loads of 3.84-7.61 Log IU/mL were included. All samples were de-identified and consequently prior treatment status for patients was unknown. Almost full open reading frame amplicons (∼ 9 kb) were generated using RT-PCR with a single primer set. The resulting amplicons were sequenced with high throughput sequencing and analysed using an in-house developed script for detecting RASs.

RESULTS

The method successfully amplified and sequenced 94% (95/101) of samples with an average coverage of 14,035; four of six failed samples were genotype 4a. Samples analysed twice yielded reproducible nucleotide frequencies across all sites. RASs were detected in 21/95 (22%) samples at a 15% threshold. The method identified one patient infected with two genotype 2b variants, and the presence of subgenomic deletion variants in 8 (8.4%) of 95 successfully sequenced samples.

CONCLUSIONS

The presented method may provide identification of HCV genotype, RASs detection, and detect multiple HCV infection without prior knowledge of sample genotype.

Identification of hepatitis C virus 2k/1b intergenotypic recombinants in Georgia

BACKGROUND AND AIMS

This study aimed to evaluate the prevalence of the hepatitis C virus intergenotype recombinant strain RF1_2k/1b in Georgia, confirm viral recombination by full genome sequencing, and determine a genetic relationship with previously described recombinant hepatitis C viruses.

METHODS

We retrospectively analysed data from 1421 Georgian patients with chronic hepatitis C. Genotyping was performed with the INNO-LiPA VERSANT HCV Genotype 2.0 Assay.

RESULTS

Virus isolates were assigned to nonspecific hepatitis C genotypes 2a/2c (n = 387) as performed by sequencing of core and NS5B genes. Subsequently, sequencing results classified the core region as genotype 2k and the NS5B region as genotype 1b for 72% (n = 280) of genotype 2 patients, corresponding to 19.7% of hepatitis C patients in Georgia. Eight samples were randomly selected for full genome sequencing which was successful in 7 of 8 samples. Analysis of the generated consensus sequences confirmed that all 7 viruses were 2k/1b recombinants, with the recombination breakpoint located within 73-77 amino acids before the NS2-NS3 junction, similar to the previously described RF1_2k/1b virus. Phylogenetic analysis revealed clustering of the Georgian 2k/1b viruses and RF1_2k/1b, suggesting that they are genetically related.

CONCLUSIONS

The 19.7% prevalence of RF1_2k/1b in Georgia patients is far higher than has generally been reported to date worldwide. Identification of recombinants in low income countries with a high prevalence of HCV infection might be reasonable for choosing the most cost-effective treatment regimens.

Global epidemiology of HCV subtypes and resistance-associated substitutions evaluated by sequencing-based subtype analyses

BACKGROUND & AIMS

HCV genotype, subtype, and presence of resistance-associated substitutions (RASs) are key determinants for the selection of direct-acting antiviral (DAA) treatment regimens. However, current HCV genotyping assays have limitations in differentiating between HCV subtypes, and RAS prevalence is largely undefined. The aim of this study was to investigate HCV epidemiology in 12,615 patient samples from 28 different countries across five geographic regions.

METHODS

We compared HCV genotype and subtypes using INNO-LiPA 2.0 vs. amplicon sequencing among 8,945 patients from phase II/III clinical trials of DAAs. Global HCV molecular epidemiology in 12,615 patients was investigated. Subtype RAS prevalence was determined by population or deep sequencing, and phylogenetic analyses investigating subtype diversity were performed.

RESULTS

Although there was high concordance between INNO-LiPA and sequencing for genotype determination, INNO-LiPA was insufficient for subtype determination for genotype 2, 3, 4, and 6. Sequencing provided subtype refinement for 42%, 10%, 81%, and 78% of genotype 2, 3, 4, and six patients, respectively. Genotype discordance (genotype 2-genotype 1) was observed in 28 of 950 (3%) genotype 2 patients, consistent with inter-genotype recombinants. Sequencing-based analyses demonstrated variations in regional subtype prevalence, notably within genotype 2, 4 and 6. RAS prevalence varied by subtype, with the clinically relevant NS3 RAS Q80K found in genotype 1a, 5a and 6a and the NS5A RAS Y93H in genotype 1b, 3a, 4b, 4r and 7.

CONCLUSIONS

Together, these analyses provide an understanding of subtyping accuracy and RAS distribution that are crucial for the implementation of global HCV treatment strategies.

LAY SUMMARY

Hepatitis C virus (HCV) is highly variable, with seven genotypes and 67 subtypes characterized to date. The aim of this study was to i) compare two different methods of discriminating between genotypes; ii) investigate the prevalence of HCV subtypes for each genotype around the world; iii) find the prevalence of resistance-associated substitutions (RASs) in different subtypes. We found that both methods showed high concordance in genotype discrimination, but specific subtypes were not always identified accurately. Sequencing-based analyses demonstrated variations in regional subtype prevalence for some genotypes, notably within GT2, 4 and 6. RAS prevalence also varied by subtype. These variations could determine how successful different drugs are for treating HCV.

Sequencing of hepatitis C virus for detection of resistance to direct-acting antiviral therapy: A systematic review

The significance of the clinical impact of direct‐acting antiviral (DAA) resistance‐associated substitutions (RASs) in hepatitis C virus (HCV) on treatment failure is unclear. No standardized methods or guidelines for detection of DAA RASs in HCV exist. To facilitate further evaluations of the impact of DAA RASs in HCV, we conducted a systematic review of RAS sequencing protocols, compiled a comprehensive public library of sequencing primers, and provided expert guidance on the most appropriate methods to screen and identify RASs. The development of standardized RAS sequencing protocols is complicated due to a high genetic variability and the need for genotype‐ and subtype‐specific protocols for multiple regions. We have identified several limitations of the available methods and have highlighted areas requiring further research and development. The development, validation, and sharing of standardized methods for all genotypes and subtypes should be a priority. (Hepatology Communications 2017;1:379–390)

Treatment of a patient with genotype 7 hepatitis C virus infection with sofosbuvir and velpatasvir

During a phase 3 study evaluating the combination of sofosbuvir-velpatasvir for 12 weeks in patients with genotype 1, 2, 4, 5, and 6 hepatitis C virus (HCV) infection, we enrolled a patient who was subsequently found to be infected with genotype 7 HCV. This patient tolerated the study regimen well and achieved sustained virological response 12 weeks after treatment (SVR12). (Hepatology 2016;64:983-985).

Next-Generation Sequencing: a Diagnostic One-Stop Shop for Hepatitis C?

Before starting chronic hepatitis C treatment, the viral genotype/subtype has to be accurately determined and potentially coupled with drug resistance testing. Due to the high genetic variability of the hepatitis C virus, this can be a demanding task that can potentially be streamlined by viral whole-genome sequencing using next-generation sequencing as demonstrated by an article in this issue of the Journal of Clinical Microbiology by E. Thomson, C. L. C. Ip, A. Badhan, M. T. Christiansen, W. Adamson, et al. (J Clin Microbiol. 54:2455-2469, 2016, http://dx.doi.org/10.1128/JCM.00330-16).

Development and Validation of a Template-Independent Next-Generation Sequencing Assay for Detecting Low-Level Resistance-Associated Variants of Hepatitis C Virus

To develop hepatitis C virus (HCV) direct-acting antiviral (DAA) drugs that can treat most HCV genotypes and offer higher barriers for treatment-resistant mutations, it is important to study resistance-associated variants (RAVs). Current commercially available RAV detection assays rely on genotype- or subtype-specific template-dependent PCR amplification. These assays are limited to genotypes and subtypes that are often prevalent in developed countries because of availability of public sequence databases. To support global clinical trials of DAAs, we developed and validated a template-independent (TI) next-generation sequencing (NGS) assay for HCV whole genome sequencing that can perform HCV subtyping, detect HCV mixed genotype or subtype infection, and identify low-level RAVs at a 5% fraction of the viral population with sensitivity and positive predictive value ≥ 0.9. We compared TI-NGS with commercial genotype- or subtype-specific Sanger sequencing assays, and found that TI-NGS both confirmed most of variants called by Sanger sequencing and avoided biases likely caused by PCR primers used in Sanger sequencing. To confirm TI-NGS assay's variant calls at the discrepant positions with Sanger sequencing, we custom designed template-dependent NGS assays and obtained 100% concordance with the TI-NGS assay. The ability to reliably detect low-level RAVs in HCV samples of any subtype without PCR primer-related bias makes this TI-NGS assay an important tool in studying HCV DAA drug resistance.

Sequencing Analysis of NS3/4A, NS5A, and NS5B Genes from Patients Infected with Hepatitis C Virus Genotypes 5 and 6

Direct-acting antivirals (DAAs) with activity against multiple genotypes of the hepatitis C virus (HCV) were recently developed and approved for standard-of-care treatment. However, sequencing assays to support HCV genotype 5 and 6 analysis are not widely available. Here, we describe the development of a sequencing assay for the NS3/4A, NS5A, and NS5B genes from HCV genotype 5 and 6 patient isolates. Genotype- and subtype-specific primers were designed to target NS3/4A, NS5A, and NS5B for cDNA synthesis and nested PCR amplification. Amplification was successfully performed for a panel of 32 plasma samples from HCV-infected genotype 5 and 6 patients with sequencing data obtained for all attempted samples. LiPA 2.0 (Versant HCV genotype 2.0) is a reverse hybridization line probe assay that is commonly used for genotyping HCV-infected patients enrolled in clinical studies. Using NS3/4A, NS5A, and NS5B consensus sequences, HCV subtypes were determined that were not available for the initial LiPA 2.0 result for genotype 6 samples. Samples amplified here included the following HCV subtypes: 5a, 6a, 6e, 6f, 6j, 6i, 6l, 6n, 6o, and 6p. The sequencing data generated allowed for the determination of the presence of variants at amino acid positions previously characterized as associated with resistance to DAAs. The simple and robust sequencing assay for genotypes 5 and 6 presented here may lead to a better understanding of HCV genetic diversity and prevalence of resistance-associated variants.

Hepatitis C virus whole genome sequencing: Current methods/issues and future challenges

Therapy for hepatitis C is currently undergoing a revolution. The arrival of new antiviral agents targeting viral proteins reinforces the need for a better knowledge of the viral strains infecting each patient. Hepatitis C virus (HCV) whole genome sequencing provides essential information for precise typing, study of the viral natural history or identification of resistance-associated variants. First performed with Sanger sequencing, the arrival of next-generation sequencing (NGS) has simplified the technical process and provided more detailed data on the nature and evolution of viral quasi-species. We will review the different techniques used for HCV complete genome sequencing and their applications, both before and after the apparition of NGS. The progress brought by new and future technologies will also be discussed, as well as the remaining difficulties, largely due to the genomic variability.

A Complex Network of Interactions between S282 and G283 of Hepatitis C Virus Nonstructural Protein 5B and the Template Strand Affects Susceptibility to Sofosbuvir and Ribavirin

The hepatitis C virus (HCV) RNA-dependent RNA-polymerase NS5B is essentially required for viral replication and serves as a prominent drug target. Sofosbuvir is a prodrug of a nucleotide analog that interacts selectively with NS5B and has been approved for HCV treatment in combination with ribavirin. Although the emergence of resistance to sofosbuvir is rarely seen in the clinic, the S282T mutation was shown to decrease susceptibility to this drug. S282T was also shown to confer hypersusceptibility to ribavirin, which is of potential clinical benefit. Here we devised a biochemical approach to elucidate the underlying mechanisms. Recent crystallographic data revealed a hydrogen bond between S282 and the 2'-hydroxyl of the bound nucleotide, while the adjacent G283 forms a hydrogen bond with the 2'-hydroxyl of the residue of the template that base pairs with the nucleotide substrate. We show that DNA-like modifications of the template that disrupt hydrogen bonding with G283 cause enzyme pausing with natural nucleotides. However, the specifically introduced DNA residue of the template reestablishes binding and incorporation of sofosbuvir in the context of S282T. Moreover, the DNA-like modifications of the template prevent the incorporation of ribavirin in the context of the wild-type enzyme, whereas the S282T mutant enables the binding and incorporation of ribavirin under the same conditions. Together, these findings provide strong evidence to show that susceptibility to sofosbuvir and ribavirin depends crucially on a network of interdependent hydrogen bonds that involve the adjacent residues S282 and G283 and their interactions with the incoming nucleotide and complementary template residue, respectively.

A method for near full-length amplification and sequencing for six hepatitis C virus genotypes

BACKGROUND

Hepatitis C virus (HCV) is a rapidly evolving RNA virus that has been classified into seven genotypes. All HCV genotypes cause chronic hepatitis, which ultimately leads to liver diseases such as cirrhosis. The genotypes are unevenly distributed across the globe, with genotypes 1 and 3 being the most prevalent. Until recently, molecular epidemiological studies of HCV evolution within the host and at the population level have been limited to the analyses of partial viral genome segments, as it has been technically challenging to amplify and sequence the full-length of the 9.6 kb HCV genome. Although recent improvements have been made in full genome sequencing methodologies, these protocols are still either limited to a specific genotype or cost-inefficient.

RESULTS

In this study we describe a genotype-specific protocol for the amplification and sequencing of the near-full length genome of all six major HCV genotypes. We applied this protocol to 122 HCV positive clinical samples, and had a successful genome amplification rate of 90%, when the viral load was greater than 15,000 IU/ml. The assay was shown to have a detection limit of 1-3 cDNA copies per reaction. The method was tested with both Illumina and PacBio single molecule, real-time (SMRT) sequencing technologies. Illumina sequencing resulted in deep coverage and allowed detection of rare variants as well as HCV co-infection with multiple genotypes. The application of the method with PacBio RS resulted in sequence reads greater than 9 kb that covered the near full-length HCV amplicon in a single read and enabled analysis of the near full-length quasispecies.

CONCLUSIONS

The protocol described herein can be utilised for rapid amplification and sequencing of the near-full length HCV genome in a cost efficient manner suitable for a wide range of applications.

Amplification and pyrosequencing of near-full-length hepatitis C virus for typing and monitoring antiviral resistant strains

Directly acting antiviral drugs have contributed considerable progress to hepatitis C virus (HCV) treatment, but they show variable activity depending on virus genotypes and subtypes. Therefore, accurate genotyping including recombinant form detection is still of major importance, as is the detection of resistance-associated mutations in case of therapeutic failure. To meet these goals, an approach to amplify the HCV near-complete genome with a single long-range PCR and sequence it with Roche GS Junior was developed. After optimization, the overall amplification success rate was 73% for usual genotypes (i.e. HCV 1a, 1b, 3a and 4a, 16/22) and 45% for recombinant forms RF_2k/1b (5/11). After pyrosequencing and subsequent de novo assembly, a near-full-length genomic consensus sequence was obtained for 19 of 21 samples. The genotype and subtype were confirmed by phylogenetic analysis for every sample, including the suspected recombinant forms. Resistance-associated mutations were detected in seven of 13 samples at baseline, in the NS3 (n = 3) or NS5A (n = 4) region. Of these samples, the treatment of one patient included daclatasvir, and that patient experienced a relapse. Virus sequences from pre- and posttreatment samples of four patients who experienced relapse after sofosbuvir-based therapy were compared: the selected variants seem too far from the NS5B catalytic site to be held responsible. Although tested on a limited set of samples and with technical improvements still necessary, this assay has proven to be successful for both genotyping and resistance-associated variant detection on several HCV types.

ve-SEQ: Robust, unbiased enrichment for streamlined detection and whole-genome sequencing of HCV and other highly diverse pathogens

The routine availability of high-depth virus sequence data would allow the sensitive detection of resistance-associated variants that can jeopardize HIV or hepatitis C virus (HCV) treatment. We introduce ve-SEQ, a high-throughput method for sequence-specific enrichment and characterization of whole-virus genomes at up to 20% divergence from a reference sequence and 1,000-fold greater sensitivity than direct sequencing. The extreme genetic diversity of HCV led us to implement an algorithm for the efficient design of panels of oligonucleotide probes to capture any sequence among a defined set of targets without detectable bias. ve-SEQ enables efficient detection and sequencing of any HCV genome, including mixtures and intra-host variants, in a single experiment, with greater tolerance of sequence diversity than standard amplification methods and greater sensitivity than metagenomic sequencing, features that are directly applicable to other pathogens or arbitrary groups of target organisms, allowing the combination of sensitive detection with sequencing in many settings.