Composition and Interactions of Hepatitis B Virus Quasispecies Defined the Virological Response During Telbivudine Therapy

Early dynamic changes of quasispecies in the reverse transcriptase region of hepatitis B virus in telbivudine treatment

BACKGROUND

Telbivudine (LdT) - a synthetic thymidine β-L-nucleoside analogue (NA) - is an effective inhibitor for hepatitis B virus (HBV) replication. The quasispecies spectra in the reverse transcriptase (RT) region of the HBV genome and their dynamic changes associated with LdT treatment remains largely unknown.

METHODS

We prospectively recruited a total of 21 treatment-naive patients with chronic HBV infection and collected sequential serum samples at five time points (baseline, weeks 1, 3, 12, and 24 after LdT treatment). The HBV RT region was amplified and shotgun-sequenced by the Ion Torrent Personal Genome Machine (PGM)® system. We reconstructed full-length haplotypes of the RT region using an integrated bioinformatics framework, including de novo contig assembly and full-length haplotype reconstruction. In addition, we investigated the quasispecies' dynamic changes and evolution history and characterized potential NAs resistant mutations over the treatment course.

RESULTS

Viral quasispecies differed obviously between patients with complete (n = 8) and incomplete/no response (n = 13) at 12 weeks after LdT treatment. A reduced d_N/d_S ratio in quasispecies demonstrated a selective constraint resulting from antiviral therapy. The temporal clustering of sequential quasispecies showed different patterns along with a 24-week observation, although its statistic did not differ significantly. Several patients harboring pre-existing resistant mutations showed different clinical responses, while NAs resistant mutations were rare within a short-term treatment.

CONCLUSION

A complete profile of quasispecies reconstructed from in-depth shotgun sequencing may has important implications for enhancing clinical decision in adjusting antiviral therapy timely.

High prevalence of tryptophan-truncated S quasispecies in treatment-naïve chronic hepatitis B patients

Hepatitis B virus surface antigen (HBsAg) encoded by the S gene is highly expressed during the replication cycle of hepatitis B virus (HBV). However, the frequent usage of tryptophan in HBsAg, which leads to a high cost of biosynthesis, is inconsistent with the high expression level of this protein. Tryptophan-truncated mutation of HBsAg, that is, a tryptophan to stop codon mutation resulting in truncated HBsAg, might help to maintain its high expression with lower biosynthetic cost. We aimed to investigate the prevalence of tryptophan-truncated S quasispecies in treatment-naïve patients with chronic hepatitis B (CHB) by applying CirSeq as well as a site-by-site algorithm developed by us to identify variants at extremely low frequencies in the carboxyl terminus of HBsAg. A total of 730 mutations were identified in 27 patients with CHB, varying from seven to 56 mutations per sample. The number of synonymous mutations was much higher than that of nonsynonymous mutations in the reverse transcriptase (RT) coding region and vice versa in the S coding region, implying that the evolutionary constraints on the RT and S genes might be different. We showed that 25 (92.6 %) of 27 patients had at least one S-truncated mutation, most of which were derived from tryptophan, indicating a high prevalence of tryptophan-truncated S mutations in treatment-naïve patients with CHB. In terms of the RT gene, 21 (77.8 %) patients had pre-existing drug-resistant mutations, while no truncated mutations were detected. Our findings that tryptophan-truncated S quasispecies and drug-resistant RT mutants were highly prevalent in treatment-naïve patients with CHB provide new insights into the composition of the HBV population, which might help optimize the treatment and management of patients with CHB.

Methodology-dependent performance of serum HBV RNA in predicting treatment outcomes in chronic hepatitis B patients

BACKGROUND

Whether different serum HBV RNA detection assays can consistently predict treatment outcomes in patients with chronic hepatitis B remains controversial.

METHODS

We enrolled 188 patients who had stopped nucleos(t)ide analogues (NAs) (STOP cohort-1, -2) and 78 receiving entecavir (ETV) therapy (ETV cohort) and used double-target (targeting both 5' and 3' ends of the HBV pregenome RNA [DT-RNA]) and three single-target (targeting the S-region [S-RNA], X-region [X-RNA], and poly-A tail of HBV RNA [PolyA-RNA]) assays to predict treatment outcomes.

RESULTS

In STOP cohorts, DT-RNA, S-RNA and X-RNA at NAs cessation showed higher predictive powers for clinical relapse (time-dependent areas under the curve [AUCs] for years 1, 2, 3, and 4 ranged between 0.724 and 0.772 in cohort-1, and between 0.741 and 0.824 in cohort-2) than the PolyA-RNA (AUCs between 0.604 and 0.611 in cohort-1; and between 0.530 and 0.584 in cohort-2). The predictive power for 2-year HBeAg loss of the four targeted RNAs in the ETV cohort at 6 months were similar (AUCs, 0.848, 0.838, 0.825, and 0.801), and superior to that of the HBV DNA level at 6 months (AUC, 0.721).

CONCLUSION

The outcome prediction performance of serum HBV RNAs is methodology-dependent. PolyA-RNA detection was not recommended to predict off-treatment relapses.

Rapid Turnover of Hepatitis B Virus Covalently Closed Circular DNA Indicated by Monitoring Emergence and Reversion of Signature-Mutation in Treated Chronic Hepatitis B Patients

BACKGROUND AND AIMS

Hepatitis B virus (HBV) covalently closed circular DNA (cccDNA) plays a pivotal role in the establishment and persistence of HBV infection. Understanding the turnover time of preexisting cccDNA pools would be helpful in designing strategies to clear HBV by fully blocking the de novo generation of cccDNA.

APPROACH AND RESULTS

In this study, we retrospectively monitored the emergence and reversion of the rtM204I/V mutant, a signature lamivudine resistance (LAMR ) mutation serving as a biomarker of cccDNA turnover in liver biopsies and longitudinal serum samples from two clinical trials. Methodologies were optimized to differentially isolate and sequence HBV virion DNA, cccDNA, and HBV RNA from clinical samples. A strong correlation was observed between LAMR composition of cccDNA with that of serum and intrahepatic HBV RNA in paired liver and serum samples (r = 0.96 and 0.90, respectively), suggesting that serum HBV RNA can serve as a surrogate marker of cccDNA genetic composition when liver biopsies are unavailable. LAMR mutations emerged and increased from undetectable to 40%-90% within 16-28 weeks in serum HBV RNA from telbivudine-treated patients experiencing virological breakthrough. Similarly, in lamivudine-resistant patients who switched to interferon therapy, serum HBV-RNA population bearing 100% LAMR mutations fully reversed back to wild type within 24-48 weeks.

CONCLUSIONS

The genetic composition dynamics of serum HBV RNA and biopsy cccDNA in treated HBV patients indicates that cccDNA turnover occurs relatively rapidly (several months), offering a possibility of HBV cure with finite therapy through completely blocking cccDNA replenishment.

Factors associated with the biphasic kinetics of serum HBV RNA in patients with HBeAg-positive chronic hepatitis B treated with nucleos(t)ide analogues

BACKGROUND

Serum hepatitis B virus (HBV) RNA is a novel biomarker for evaluating treatment response. Detailed information regarding serum HBV RNA kinetics during treatment with nucleos(t)ide analogues (NAs) is limited.

AIMS

To ascertain serum HBV RNA kinetics during long-term NA treatment and identify associated factors.

METHODS

We enrolled 76 HBeAg-positive chronic hepatitis B patients receiving NA from randomised controlled trials. Laboratory assays were undertaken every 3 months. Factors associated with serum HBV RNA kinetics were identified by generalised estimating equations.

RESULTS

Baseline serum HBV RNA was 8.5 ± 1.0 log₁₀  copies/mL. Decline in serum HBV RNA during NA therapy was biphasic: the first phase (HBV DNA detectable) had a fast decrease (median slope, -0.207 log₁₀  copies/mL/month) and was followed by a second phase (HBV DNA undetectable) with slow decrease (median slope, -0.071 log₁₀  copies/mL/month). In the first phase, factors independently associated with lower initial serum HBV RNA were male sex (OR, 0.685, P = 0.044), low baseline HBsAg (OR, 0.525, P = 0.001) and rapid virological response (RVR) (OR, 0.624, P = 0.031). In the second phase, only RVR was independently associated with serum HBV RNA kinetics, including its lower initial level (OR, 0.694, P = 0.043) and greater decline (OR, 0.966, P = 0.002). Based on viral dynamics, time needed to achieve undetectable serum HBV RNA from baseline was 43.56 (IQR: 29.49-66.40) months.

CONCLUSION

RVR was a significant determinant for biphasic decline in serum HBV RNA during NA treatment, which significantly influenced the treatment duration required to achieve undetectable serum HBV RNA.

Viral quasispecies of hepatitis B virus in patients with YMDD mutation and lamivudine resistance may not predict the efficacy of lamivudine/adefovir rescue therapy

The association between hepatitis B virus (HBV) quasispecies (QS) and the efficacy of nucleos(t)ide analog therapy is currently not well defined, particularly in the case of lamivudine (LAM)/adefovir (ADV) combination rescue therapy for patients with chronic HBV infection (CHB) presenting with LAM resistance. In the present study, 16 CHB patients with the rtM204I/V mutation in the tyrosine-methionine-aspartate-aspartate motif of the C domain of the polymerase gene who switched to LAM/ADV treatment due to LAM resistance were assessed. HBV DNA was isolated from these patients and the reverse transcriptase (RT) region was sequenced. The QS heterogeneity and distribution was analyzed, the mutation sites were recorded and the phylogenetic trees were constructed. The results indicated that QS heterogeneity and distribution in the RT and S regions were not significantly different between responders (RS) and non-RS (NRS) at baseline (P>0.05), except for the higher frequency of a dominant strain in the RT region at the nucleotide level in the RS group (P=0.039). In addition, in NRS, no significant difference in QS heterogeneity or distribution in these regions was identified at six months vs. the baseline. Furthermore, although in the non-responder group the frequency of the LAM resistance-associated mutations (rtM204V/I) decreased at 6 months compared with the baseline, it did not disappear in any of the patients after six months of treatment. Analysis of individual patients did not indicate any consistent selection of specific HBV mutants during LAM/ADV rescue therapy. In conclusion, the baseline HBV QS within the RT and S regions may not be a valid predictor of the response to LAM/ADV rescue treatment in CHB patients with LAM resistance.

Deep sequencing of HBV pre-S region reveals high heterogeneity of HBV genotypes and associations of word pattern frequencies with HCC

Hepatitis B virus (HBV) infection is a common problem in the world, especially in China. More than 60–80% of hepatocellular carcinoma (HCC) cases can be attributed to HBV infection in high HBV prevalent regions. Although traditional Sanger sequencing has been extensively used to investigate HBV sequences, NGS is becoming more commonly used. Further, it is unknown whether word pattern frequencies of HBV reads by Next Generation Sequencing (NGS) can be used to investigate HBV genotypes and predict HCC status. In this study, we used NGS to sequence the pre-S region of the HBV sequence of 94 HCC patients and 45 chronic HBV (CHB) infected individuals. Word pattern frequencies among the sequence data of all individuals were calculated and compared using the Manhattan distance. The individuals were grouped using principal coordinate analysis (PCoA) and hierarchical clustering. Word pattern frequencies were also used to build prediction models for HCC status using both K-nearest neighbors (KNN) and support vector machine (SVM). We showed the extremely high power of analyzing HBV sequences using word patterns. Our key findings include that the first principal coordinate of the PCoA analysis was highly associated with the fraction of genotype B (or C) sequences and the second principal coordinate was significantly associated with the probability of having HCC. Hierarchical clustering first groups the individuals according to their major genotypes followed by their HCC status. Using cross-validation, high area under the receiver operational characteristic curve (AUC) of around 0.88 for KNN and 0.92 for SVM were obtained. In the independent data set of 46 HCC patients and 31 CHB individuals, a good AUC score of 0.77 was obtained using SVM. It was further shown that 3000 reads for each individual can yield stable prediction results for SVM. Thus, another key finding is that word patterns can be used to predict HCC status with high accuracy. Therefore, our study shows clearly that word pattern frequencies of HBV sequences contain much information about the composition of different HBV genotypes and the HCC status of an individual.

HBV infection can lead to many liver complications including hepatocellular carcinoma (HCC), one of the most common liver cancers in China. High-throughput sequencing technologies have recently been used to study the genotype sequence compositions of HBV infected individuals and to distinguish chronic HBV (CHB) infection from HCC. We used NGS to sequence the pre-S region of a large number of CHB and HCC individuals and designed novel word pattern based approaches to analyze the data. We have several surprising key findings. First, most HBV infected individuals contained mixtures of genotypes B and C sequences. Second, multi-dimensional scaling (MDS) analysis of the data showed that the first principal coordinate was closely associated with the fraction of genotype B (or C) sequences and the second principal coordinate was highly associated with the probability of HCC. Third, we also designed K-nearest neighbors (KNN) and support vector machine (SVM) based classifiers for CHB and HCC with high prediction accuracy. The results were validated in an independent data set.

Deep sequencing in the management of hepatitis virus infections

The hepatitis viruses represent a major public health problem worldwide. Procedures for characterization of the genomic composition of their populations, accurate diagnosis, identification of multiple infections, and information on inhibitor-escape mutants for treatment decisions are needed. Deep sequencing methodologies are extremely useful for these viruses since they replicate as complex and dynamic quasispecies swarms whose complexity and mutant composition are biologically relevant traits. Population complexity is a major challenge for disease prevention and control, but also an opportunity to distinguish among related but phenotypically distinct variants that might anticipate disease progression and treatment outcome. Detailed characterization of mutant spectra should permit choosing better treatment options, given the increasing number of new antiviral inhibitors available. In the present review we briefly summarize our experience on the use of deep sequencing for the management of hepatitis virus infections, particularly for hepatitis B and C viruses, and outline some possible new applications of deep sequencing for these important human pathogens.

Clinical features and viral quasispecies characteristics associated with infection by the hepatitis B virus G145R immune escape mutant

Coexistence of the hepatitis B surface antigen (HBsAg) and hepatitis B surface antibody (anti-HBs) is an uncommon phenomenon, and the underlying mechanisms remain largely unknown. Amino-acid (aa) substitution from glycine to arginine at aa 145 (G145R), in the major hydrophilic region, has been reported in patients with HBsAg and anti-HBs coexistence. However, there is limited knowledge about the clinical features and viral quasispecies characteristics associated with G145R mutant hepatitis B virus (HBV) infection. We herein describe the dynamic changes in the serological and virological markers in a case of hepatitis B with coexisting HBsAg and anti-HBs, caused by a G145R immune escape mutant (genotype C). Entecavir was administered during the 4th week after admission. Alanine aminotransferase peaked in the 16th week, while both the HBsAg and HBeAg declined rapidly. HBsAg clearance and hepatitis B e antigen (HBeAg)/hepatitis B e antibody (anti-HBe) seroconversion were achieved in the 36th week, and then entecavir was withdrawn. A follow-up of 96 weeks showed that HBV DNA remained undetectable and that anti-HBs was maintained above 100 mIU/mL. The quasispecies characteristics of the G145R mutant HBV were investigated via ultra-deep sequencing. The complexity and genetic distance of the S and RT regions were much higher in the 8th week than at baseline or in the 4th week. Moreover, the frequencies of mutations (L173P, Q181R and A184V) in cytotoxic T lymphocyte epitopes increased before entecavir treatment. These findings extend understanding of the evolution of HBV under host immune pressure and of the clinical outcomes of affected patients.

Clinical and biological insights from viral genome sequencing

Whole-genome sequencing (WGS) of pathogens is becoming increasingly important not only for basic research but also for clinical science and practice. In virology, WGS is important for the development of novel treatments and vaccines, and for increasing the power of molecular epidemiology and evolutionary genomics. In this Opinion article, we suggest that WGS of viruses in a clinical setting will become increasingly important for patient care. We give an overview of different WGS methods that are used in virology and summarize their advantages and disadvantages. Although there are only partially addressed technical, financial and ethical issues in regard to the clinical application of viral WGS, this technique provides important insights into virus transmission, evolution and pathogenesis.

Hepatitis B virus resistance substitutions: long-term analysis by next-generation sequencing

HBV phylogenetics and resistance-associated mutations (RAMs) were surveyed by next-generation sequencing of 21 longitudinal samples from seven patients entering antiviral therapy. The virus populations were dominated by a few abundant lineages that coexisted with substantial numbers of low-frequency variants. A few low-frequency RAMs were observed before treatment, but new ones emerged, and their frequencies increased during therapy. Together, these results support the idea that chronic HBV infection is dominated by a few virus lineages and that an accompanying plethora of diverse, low-frequency variants may function as a reservoir that potentially contribute to viral genetic plasticity, potentially affecting patient outcome.

Small surface antigen variants of HBV associated with responses to telbivudine treatment in chronic hepatitis B patients

BACKGROUND

Nucleoside/nucleotide analogues are widely used to treat chronic HBV infection, but drug resistance is common. The role of HBV surface gene variants in drug resistance to nucleoside/nucleotide analogues is unknown. We are trying to investigate the dynamics of S gene mutations and how they relate to a patient's virological response in this study.

METHODS

Thirty patients with chronic hepatitis B were enrolled and serum samples were collected at multiple time points during treatment with telbivudine (LdT). The coding regions of the small surface antigen (S-HBsAg) were amplified and sequenced using the 454 GS FLX⁺ System.

RESULTS

Sequencing results revealed different dynamics of non-synonymous mutations, such as sL9P, sN40S, sG44E, sW172*, sW182* and sS187F, between patients with a complete virological response and those with a partial virological response. The viral population heterogeneity decreased at week 12 of LdT treatment in patients with a complete virological response, with a concomitant decline in non-synonymous mutations (from an average of 14 to 9.9 per sample) and an increase in the frequencies of major variants (from 14.3% to 40.4%).

CONCLUSIONS

Our findings suggest that the decrease in viral population heterogeneity at an early stage of LdT treatment was associated with the subsequent optimal virological response, and the early appearance of some specific mutations, such as sG44E, sW172* and sW182*, is a potential indicator of a partial virological response in continuing therapy.

Different Variants in Reverse Transcriptase Domain Determined by Ultra-deep Sequencing in Treatment-naïve and Treated Indonesian Patients Infected with Hepatitis B Virus

A nucleos(t)ide analog (NA) is the common antiviral drug available for directly treating hepatitis B virus (HBV) infection. However, its application has led to the emergence of NA-resistant mutations mostly in a conserved region of the reverse transcriptase domain of HBV polymerase. Harboring NA-resistant mutations decreases drug effectiveness and increases the frequency of end-stage liver disease. The invention of next-generation sequencing that can generate thousands of sequences from viral complex mixtures provides opportunities to detect minor changes and early viral evolution under drug stress. The present study used ultra-deep sequencing to evaluate discrepant quasispecies in the reverse transcriptase domain of HBV including NA-resistant hotspots between seven treatment-naïve Indonesian patients infected with HBV and five at the early phase of treatment. The most common sub-genotype was HBV B3 (83.34%). The substitution rate of variants determined among amino acids with a ratio of ≥ 1% changes was higher among the population in conserved regions (23.19% vs. 4.59%, P = 0.001) and in the inter-reverse transcriptase domain (23.95% vs. 2.94%, P = 0.002) in treatment naïve, than in treated patients. Nine hotspots of antiviral resistance were identified in both groups, and the mean frequency of changes in all patients was < 1%. The known rtM204I mutation was the most frequent in both groups. The lower rate of variants in HBV quasispecies in patients undergoing treatment could be associated with virus elimination and the extinction of sensitive species by NA therapy. The present findings imply that HBV quasispecies dynamically change during treatment.