Dynamics of defective hepatitis C virus clones in reinfected liver grafts in liver transplant recipients: ultradeep sequencing analysis

Single-molecule sequencing of the whole HCV genome revealed envelope deletions in decompensated cirrhosis associated with NS2 and NS5A mutations

BACKGROUND

Defective hepatitis C virus (HCV) genomes with deletion of the envelope region have been occasionally reported by short-read sequencing analyses. However, the clinical and virological details of such deletion HCV have not been fully elucidated.

METHODS

We developed a highly accurate single-molecule sequencing system for full-length HCV genes by combining the third-generation nanopore sequencing with rolling circle amplification (RCA) and investigated the characteristics of deletion HCV through the analysis of 21 patients chronically infected with genotype-1b HCV.

RESULT

In 5 of the 21 patients, a defective HCV genome with approximately 2000 bp deletion from the E1 to NS2 region was detected, with the read frequencies of 34-77%, suggesting the trans-complementation of the co-infecting complete HCV. Deletion HCV was found exclusively in decompensated cirrhosis (5/12 patients), and no deletion HCV was observed in nine compensated patients. Comparing the amino acid substitutions between the deletion and complete HCV (DAS, deletion-associated substitutions), the deletion HCV showed higher amino acid mutations in the ISDR (interferon sensitivity-determining region) in NS5A, and also in the TMS (transmembrane segment) 3 to H (helix) 2 region of NS2.

CONCLUSIONS

Defective HCV genome with deletion of envelope genes is associated with decompensated cirrhosis. The deletion HCV seems susceptible to innate immunity, such as endogenous interferon with NS5A mutations, escaping from acquired immunity with deletion of envelope proteins with potential modulation of replication capabilities with NS2 mutations. The relationship between these mutations and liver damage caused by HCV deletion is worth investigating.

Analysis of direct-acting antiviral-resistant hepatitis C virus haplotype diversity by single-molecule and long-read sequencing

The method of analyzing individual resistant hepatitis C virus (HCV) by a combination of haplotyping and resistance-associated substitution (RAS) has not been fully elucidated because conventional sequencing has only yielded short and fragmented viral genomes. We performed haplotype analysis of HCV mutations in 12 asunaprevir/daclatasvir treatment-failure cases using the Oxford Nanopore sequencer. This enabled single-molecule long-read sequencing using rolling circle amplification (RCA) for correction of the sequencing error. RCA of the circularized reverse-transcription polymerase chain reaction products successfully produced DNA longer than 30 kilobase pairs (kb) containing multiple tandem repeats of a target 3 kb HCV genome. The long-read sequencing of these RCA products could determine the original sequence of the target single molecule as the consensus nucleotide sequence of the tandem repeats and revealed the presence of multiple viral haplotypes with the combination of various mutations in each host. In addition to already known signature RASs, such as NS3-D168 and NS5A-L31/Y93, there were various RASs specific to a different haplotype after treatment failure. The distribution of viral haplotype changed over time; some haplotypes disappeared without acquiring resistant mutations, and other haplotypes, which were not observed before treatment, appeared after treatment. Conclusion: The combination of various mutations other than the known signature RAS was suggested to influence the kinetics of individual HCV quasispecies in the direct-acting antiviral treatment. HCV haplotype dynamic analysis will provide novel information on the role of HCV diversity within the host, which will be useful for elucidating the pathological mechanism of HCV-related diseases.

Single-molecular real-time deep sequencing reveals the dynamics of multi-drug resistant haplotypes and structural variations in the hepatitis C virus genome

While direct-acting antivirals (DAAs) for hepatitis C virus (HCV) have dramatically progressed, patients still suffer from treatment failures. For the radical eradication of HCV, a deeper understanding of multiple resistance-associated substitutions (RASs) at the single-clone level is essential. To understand HCV quasispecies and their dynamics during DAA treatment, we applied single-molecule real-time (SMRT) deep sequencing on sera from 12 patients with genotype-1b HCV infections with DAA treatment failures, both pre- and post-treatment. We identified >3.2 kbp sequences between NS3 and NS5A genes of 187,539 clones in total, classifying into haplotype codes based on the linkage of seven RAS loci. The number of haplotype codes during the treatment, per sample, significantly decreased from 14.67 ± 9.12 to 6.58 ± 7.1, while the number of nonsynonymous codons on the seven RAS loci, per clone, significantly increased from 1.50 ± 0.92 to 3.64 ± 0.75. In five cases, the minority multi-drug resistant haplotypes at pre-treatment were identical to the major haplotypes at relapse. Moreover, various structural variations (SVs) were detected and their dynamics analysed. These results suggest that SMRT deep sequencing is useful for detecting minority haplotypes and SVs, and to evaluate the dynamics of viral genomes at the single-clone level.

HCV Defective Genomes Promote Persistent Infection by Modulating the Viral Life Cycle

Defective interfering (DI) RNAs have been detected in several human viruses. HCV in-frame deletions mutants (IFDMs), missing mainly the envelope proteins, have been found in patient sera and liver tissues. IFDMs replicate independently and can be trans-packaged into infectious virions in the presence of full length viral genome. So far, their biological role is unclear. In this study, we have isolated and cloned IFDMs from sera samples and liver tissues of patients infected with HCV genotypes 1b, 2a, and 3a. IFDMs were present in up to 26% of samples tested. Using the in vitro HCV cell culture system, co-expression of the wild type (wt) HCV replicon with HCV IFDMs RNA resulted in increased HCV replication. Additionally, co-transfection of the HCV full length genome RNA and a defective mutant missing the envelope region led to increased viral release, collectively suggesting an important biological role for IFDMs in the virus life cycle. Recently, exosomes, masters of intercellular communication, have been implicated in the transport of HCV viral genomes. We report for the first time that exosomal RNA isolated from HCV sera samples contains HCV defective genomes. We also demonstrate that inhibition of exosomal biogenesis and release influences HCV viral replication. Overall, we provide evidence that the presence of HCV IFDMs affects both viral replication and release. IFDMs exploit exosomes as means of transport, a way to evade the immune system, to spread more efficiently and possibly maintain persistent infection.

DI-tector: defective interfering viral genomes' detector for next-generation sequencing data

Defective interfering (DI) genomes, or defective viral genomes (DVGs), are truncated viral genomes generated during replication of most viruses, including live viral vaccines. Among these, "panhandle" or copy-back (cb) and "hairpin" or snap-back (sb) DI genomes are generated during RNA virus replication. 5' cb/sb DI genomes are highly relevant for viral pathogenesis since they harbor immunostimulatory properties that increase virus recognition by the innate immune system of the host. We have developed DI-tector, a user-friendly and freely available program that identifies and characterizes cb/sb genomes from next-generation sequencing (NGS) data. DI-tector confirmed the presence of 5' cb genomes in cells infected with measles virus (MV). DI-tector also identified a novel 5' cb genome, as well as a variety of 3' cb/sb genomes whose existence had not previously been detected by conventional approaches in MV-infected cells. The presence of these novel cb/sb genomes was confirmed by RT-qPCR and RT-PCR, validating the ability of DI-tector to reveal the landscape of DI genome population in infected cell samples. Performance assessment using different experimental and simulated data sets revealed the robust specificity and sensitivity of DI-tector. We propose DI-tector as a universal tool for the unbiased detection of DI viral genomes, including 5' cb/sb DI genomes, in NGS data.

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.

Persistent RNA virus infections: do PAMPS drive chronic disease?

Chronic disease associated with persistent RNA virus infections represents a key public health concern. While human immunodeficiency virus-1 and hepatitis C virus are perhaps the most well-known examples of persistent RNA viruses that cause chronic disease, evidence suggests that many other RNA viruses, including re-emerging viruses such as chikungunya virus, Ebola virus and Zika virus, establish persistent infections. The mechanisms by which RNA viruses drive chronic disease are poorly understood. Here, we discuss how the persistence of viral RNA may drive chronic disease manifestations via the activation of RNA sensing pathways.

Understanding the complex evolution of rapidly mutating viruses with deep sequencing: Beyond the analysis of viral diversity

With the advent of affordable deep sequencing technologies, detection of low frequency variants within genetically diverse viral populations can now be achieved with unprecedented depth and efficiency. The high-resolution data provided by next generation sequencing technologies is currently recognised as the gold standard in estimation of viral diversity. In the analysis of rapidly mutating viruses, longitudinal deep sequencing datasets from viral genomes during individual infection episodes, as well as at the epidemiological level during outbreaks, now allow for more sophisticated analyses such as statistical estimates of the impact of complex mutation patterns on the evolution of the viral populations both within and between hosts. These analyses are revealing more accurate descriptions of the evolutionary dynamics that underpin the rapid adaptation of these viruses to the host response, and to drug therapies. This review assesses recent developments in methods and provide informative research examples using deep sequencing data generated from rapidly mutating viruses infecting humans, particularly hepatitis C virus (HCV), human immunodeficiency virus (HIV), Ebola virus and influenza virus, to understand the evolution of viral genomes and to explore the relationship between viral mutations and the host adaptive immune response. Finally, we discuss limitations in current technologies, and future directions that take advantage of publically available large deep sequencing datasets.

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.

Quasispecies dynamics in hepatitis C liver transplant recipients receiving grafts from hepatitis C virus infected donors

The allocation of liver grafts from hepatitis C virus (HCV)-positive donors in HCV-infected liver transplant (LT) recipients leads to infection with two different viral populations. In a previous study, we examined quasispecies dynamics during reinfection by clonal sequencing, which did not allow an accurate characterization of coexistence and competition events. To overcome this limitation, here we used deep-sequencing analysis of a fragment of the HCV NS5B gene in six HCV-infected LT recipients who received HCV-infected grafts. Successive expansions and contractions of quasispecies complexity were observed, evolving in all cases towards a more homogeneous population. The population that became dominant was the one displaying the highest mutant spectrum complexity. In four patients, coexistence of minority mutants, derived from the donor or the recipient, were detected. In conclusion, our study shows that, during reinfection with a different HCV strain in LT recipients, the viral population with the highest diversity always becomes dominant.

Hepatitis C Virus Deletion Mutants Are Found in Individuals Chronically Infected with Genotype 1 Hepatitis C Virus in Association with Age, High Viral Load and Liver Inflammatory Activity

Hepatitis C virus (HCV) variants characterized by genomic deletions in the structural protein region have been sporadically detected in liver and serum of hepatitis C patients. These defective genomes are capable of autonomous RNA replication and are packaged into infectious viral particles in cells co-infected with the wild-type virus. The prevalence of such forms in the chronically HCV-infected population and the impact on the severity of liver disease or treatment outcome are currently unknown. In order to determine the prevalence of HCV defective variants and to study their association with clinical characteristics, a screening campaign was performed on pre-therapy serum samples from a well-characterized cohort of previously untreated genotype 1 HCV-infected patients who received treatment with PEG-IFNα and RBV. 132 subjects were successfully analyzed for the presence of defective species exploiting a long-distance nested PCR assay. HCV forms with deletions predominantly affecting E1, E2 and p7 proteins were found in a surprising high fraction of the subjects (25/132, 19%). Their presence was associated with patient older age, higher viral load and increased necroinflammatory activity in the liver. While the presence of circulating HCV carrying deletions in the E1-p7 region did not appear to significantly influence sustained virological response rates to PEG-IFNα/RBV, our study indicates that the presence of these subgenomic HCV mutants could be associated with virological relapse in patients who did not have detectable viremia at the end of the treatment.

Spatiotemporal Reconstruction of the Introduction of Hepatitis C Virus into Scotland and Its Subsequent Regional Transmission

A more comprehensive understanding of hepatitis C virus (HCV) transmission dynamics could facilitate public health initiatives to reduce the prevalence of HCV in people who inject drugs. We aimed to determine how HCV sequences entered and spread throughout Scotland and to identify transmission hot spots. A Scottish data set with embedded demographic data was created by sequencing the NS5B of 125 genotype 1a (Gt1a) samples and 166 Gt3a samples and analyzed alongside sequences from public databases. Applying Bayesian inference methods, we reconstructed the global origin and local spatiotemporal dissemination of HCV in Scotland. Scottish sequences mainly formed discrete clusters interspersed between sequences from the rest of the world; the most recent common ancestors of these clusters dated to 1942 to 1952 (Gt1a) and 1926 to 1942 (Gt3a), coincident with global diversification and distribution. Extant Scottish sequences originated in Edinburgh (Gt1a) and Glasgow (Gt3a) in the 1970s, but both genotypes spread from Glasgow to other regions. The dominant Gt1a strain differed between Edinburgh (cluster 2 [C2]), Glasgow (C3), and Aberdeen (C4), whereas significant Gt3a strain specificity occurred only in Aberdeen. Specific clusters initially formed separate transmission zones in Glasgow that subsequently overlapped, occasioning city-wide cocirculation. Transmission hot spots were detected with 45% of samples from patients residing in just 9 of Glasgow's 57 postcode districts. HCV was introduced into Scotland in the 1940s, concomitant with its worldwide dispersal likely arising from global-scale historical events. Cluster-specific transmission hubs were identified in Glasgow, the key Scottish city implicated in HCV dissemination. This fine-scale spatiotemporal reconstruction improves understanding of HCV transmission dynamics in Scotland.

IMPORTANCE

HCV is a major health burden and the leading cause of hepatocellular carcinoma. Public health needle exchange and "treatment as prevention" strategies targeting HCV are designed to reduce prevalence of the virus in people who inject drugs (PWID), potentially mitigating the future burden of HCV-associated liver disease. Understanding HCV transmission dynamics could increase the effectiveness of such public health initiatives by identifying and targeting regions playing a central role in virus dispersal. In this study, we examined HCV transmission in Scotland by analyzing the genetic relatedness of strains from PWID alongside data inferring the year individuals became infected and residential information at a geographically finer-scale resolution than in previous studies. Clusters of Scotland-specific strains were identified with regional specificity, and mapping the spread of HCV allowed the identification of key areas central to HCV transmission in Scotland. This research provides a basis for identifying HCV transmission hot spots.