Evolutional transition of HBV genome during the persistent infection determined by single-molecule real-time sequencing.

BACKGROUND: Although HBV infection is a serious health issue worldwide, the landscape of HBV genome dynamics in the host has not yet been clarified. This study aimed to determine the continuous genome sequence of each HBV clone using a single-molecule real-time sequencing platform, and clarify the dynamics of structural abnormalities during persistent HBV infection without antiviral therapy. PATIENTS AND METHODS: Twenty-five serum specimens were collected from 10 untreated HBV-infected patients. Continuous whole-genome sequencing of each clone was performed using a PacBio Sequel sequencer; the relationship between genomic variations and clinical information was analyzed. The diversity and phylogeny of the viral clones with structural variations were also analyzed. RESULTS: The whole-genome sequences of 797,352 HBV clones were determined. The deletion was the most common structural abnormality and concentrated in the preS/S and C regions. Hepatitis B e antibody (anti-HBe)-negative samples or samples with high alanine aminotransferase levels have significantly diverse deletions than anti-HBe-positive samples or samples with low alanine aminotransferase levels. Phylogenetic analysis demonstrated that various defective and full-length clones evolve independently and form diverse viral populations. CONCLUSIONS: Single-molecule real-time long-read sequencing revealed the dynamics of genomic quasispecies during the natural course of chronic HBV infections. Defective viral clones are prone to emerge under the condition of active hepatitis, and several types of defective variants can evolve independently of the viral clones with the full-length genome.

Mutational spectrum of hepatitis C virus in patients with chronic hepatitis C determined by single molecule real-time sequencing.

The emergence of hepatitis C virus (HCV) with resistance-associated substitution (RAS), produced by mutations in the HCV genome, is a major problem in direct acting antivirals (DAA) treatment. This study aimed to clarify the mutational spectrum in HCV-RNA and the substitution pattern for the emergence of RASs in patients with chronic HCV infection. HCV-RNA from two HCV replicon cell lines and the serum HCV-RNA of four non-liver transplant and four post-liver transplant patients with unsuccessful DAA treatment were analyzed using high-accuracy single-molecule real-time long-read sequencing. Transition substitutions, especially A>G and U>C, occurred prominently under DAAs in both non-transplant and post-transplant patients, with a mutational bias identical to that occurring in HCV replicon cell lines during 10-year culturing. These mutational biases were reproduced in natural courses after DAA treatment. RASs emerged via both transition and transversion substitutions. NS3-D168 and NS5A-L31 RASs resulted from transversion mutations, while NS5A-Y93 RASs was caused by transition substitutions. The fidelity of the RNA-dependent RNA polymerase, HCV-NS5B, produces mutational bias in the HCV genome, characterized by dominant transition mutations, notably A>G and U>C substitutions. However, RASs are acquired by both transition and transversion substitutions, and the RASs-positive HCV clones are selected and proliferated under DAA treatment pressure.

Association of HLA-DRB1*04:05 allele with drug-induced interstitial lung disease in Japanese population.

Drug-induced interstitial lung disease (DILD) is a life-threatening adverse reaction. The Japanese population is more susceptible to DILD as compared with other populations, suggesting its pathogenesis could vary depending on ethnic genetic background. We conducted case-control studies to elucidate the association between DILD and HLA alleles in the Japanese. The 177 clinically diagnosed DILD patients and 3002 healthy controls for exploration and 55 DILD patients and 201 healthy controls for validation were genotyped for four HLA genes. HLA-DRB1*04:05 was significantly associated with DILD (corrected p = 0.014); this was also validated in the other set of patients/controls. Chemical drugs other than protein therapeutics showed this association (p = 1.7 × 10-4) . The Japanese population showed a higher HLA-DRB1*04:05 frequency than most other populations. In conclusion, HLA-DRB1*04:05 could be associated with DILD susceptibility in Japanese individuals, and its high general frequency may explain the high reported incidence of DILD in Japanese.

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.

Exploring the hepatitis C virus genome using single molecule real-time sequencing.

Single molecular real-time (SMRT) sequencing, also called third-generation sequencing, is a novel sequencing technique capable of generating extremely long contiguous sequence reads. While conventional short-read sequencing cannot evaluate the linkage of nucleotide substitutions distant from one another, SMRT sequencing can directly demonstrate linkage of nucleotide changes over a span of more than 20 kbp, and thus can be applied to directly examine the haplotypes of viruses or bacteria whose genome structures are changing in real time. In addition, an error correction method (circular consensus sequencing) has been established and repeated sequencing of a single-molecule DNA template can result in extremely high accuracy. The advantages of long read sequencing enable accurate determination of the haplotypes of individual viral clones. SMRT sequencing has been applied in various studies of viral genomes including determination of the full-length contiguous genome sequence of hepatitis C virus (HCV), targeted deep sequencing of the HCV NS5A gene, and assessment of heterogeneity among viral populations. Recently, the emergence of multi-drug resistant HCV viruses has become a significant clinical issue and has been also demonstrated using SMRT sequencing. In this review, we introduce the novel third-generation PacBio RSII/Sequel systems, compare them with conventional next-generation sequencers, and summarize previous studies in which SMRT sequencing technology has been applied for HCV genome analysis. We also refer to another long-read sequencing platform, nanopore sequencing technology, and discuss the advantages, limitations and future perspectives in using these third-generation sequencers for HCV genome analysis.