Analysis of potential antiviral resistance mutation profiles within the HBV reverse transcriptase in untreated chronic hepatitis B patients using an ultra-deep pyrosequencing method

Whole genome analysis of hepatitis B virus before and during long-term therapy in chronic infected patients: Molecular characterization, impact on treatment and liver disease progression

Hepatitis B virus (HBV) infection remains a serious public health concern worldwide despite the availability of an efficient vaccine and the major improvements in antiviral treatments. The aim of the present study is to analyze the mutational profile of the HBV whole genome in ETV non-responder chronic HBV patients, in order to investigate antiviral drug resistance, immune escape, and liver disease progression to Liver Cirrhosis (LC) or Hepatocellular Carcinoma (HCC). Blood samples were collected from five chronic hepatitis B patients. For each patient, two plasma samples were collected, before and during the treatment. Whole genome sequencing was performed using Sanger technology. Phylogenetic analysis comparing the studied sequences with reference ones was used for genotyping. The mutational profile was analyzed by comparison with the reference sequence M32138. Genotyping showed that the studied strains belong to subgenotypes D1, D7, and D8. The mutational analysis showed high genetic variability. In the RT region of the polymerase gene, 28 amino acid (aa) mutations were detected. The most significant mutations were the pattern rtL180M + rtS202G + rtM204V, which confer treatment resistance. In the S gene, 35 mutations were detected namely sP120T, sT126S, sG130R, sY134F, sS193L, sI195M, and sL216stop were previously described to lead to vaccine, immunotherapy, and/or diagnosis escape. In the C gene, 34 mutations were found. In particular, cG1764A, cC1766G/T, cT1768A, and cC1773T in the BCP; cG1896A and cG1899A in the precore region and cT12S, cE64D, cA80T, and cP130Q in the core region were associated with disease progression to LC and/or HCC. Other mutations were associated with viral replication increase including cT1753V, cG1764A/T, cC1766G/T, cT1768A, and cC1788G in the BCP as well as cG1896A and cG1899A in the precore region. In the X gene, 30 aa substitutions were detected, of which substitutions xT36D, xP46S, xA47T, xI88F, xA102V, xI127T, xK130M, xV131I, and xF132Y were previously described to lead to LC and/or HCC disease progression. In conclusion, our results show high genetic variability in the long-term treatment of chronic HBV patients causing several effects. This could contribute to guiding national efforts to optimize relevant HBV treatment management in order to achieve the global hepatitis elimination goal by 2030.

Characterization of mutations in the reverse transcriptase region of hepatitis B virus in treated and untreated chronic hepatitis B patients

BACKGROUND

The reverse transcriptase (RT) region of the hepatitis B virus (HBV) is the target of antiviral treatment. However, the discrepancy in RT mutations between nucleos(t)ide analogue (NA)-treated and -untreated chronic hepatitis B (CHB) patients is un clear.

METHODS

Serum samples were collected from 119 NA-treated and 135 NA-untreated patients. The sampling time was decided by the clinician. Full-length HBV RT regions were amplified using nest polymerase chain reaction. The mutations within the RT region were analysed by direct sequencing.

RESULTS

The incidence of RT mutations in treated patients was higher than that in untreated patients (p<0.05). The classic drug-resistant mutations were detected in 44.5% (53/119) of treated patients, which was significantly higher than in untreated patients (6.7% [9/135]) (p<0.05). The non-classical mutations showed their complexity and diversity in both patient groups. Multiple mutations (three or more) were more frequent in treated patients than in untreated patients (p<0.05). Several novel mutations might be related to NA resistance.

CONCLUSIONS

The selection pressures of NAs accelerated the development of RT mutations, especially within the functional domain. Mutations in the RT region occurred not only at classical sites, but also at other non-classical sites, which might be related to drug resistance and/or viral replication. The biological function and fitness of HBV isolates harbouring these novel mutations need further in vitro and in vivo verification experiments.

Prevalence of Potential Resistance Related Variants Among Chinese Chronic Hepatitis B Patients Not Receiving Nucleos(T)ide Analogues

Background and Aims

Potential drug resistance (DR) related variants in the hepatitis B virus (HBV) reverse transcriptase (RT) region may be associated with the effectiveness of antiviral drugs and disease progression. The aim of this study was to investigate the prevalence and clinical characteristics of potential DR-related variants in Chinese CHB patients not receiving nucleos(t)ide analogues (NAs).

Patients and Methods

Two hundred and six untreated CHB patients from Huzhou Central Hospital in eastern China were recruited for this study. The serum DNA was extracted and the HBV RT region was amplified using nest polymerase chain reaction (nest-PCR). The 42 potential DR-related variants were analyzed by direct sequencing.

Results

Among these CHB patients, HBV genotype B and genotype C were identified in 121 (58.7%) and 85 (41.3%) patients, respectively. Potential DR-related variants were detected in 42.7% (88/206) of patients. Primary and secondary DR variants were found in 7.3% (15/206) of patients, including rtL80I/V, rtI169T, rtV173L rtL180M, rtA181T/V, rtM204I/V, and rtN236T. The variants at rt53, rt82, rt221, rt233, rt237, and rt256 were specific for genotype B, and those at rt38, rt84, rt126, rt139, rt153, rt191, rt214, rt238, and rt242 were specific for genotype C. Moreover, the variation frequency in the A-B interdomain (3.96%) was significantly higher than that in the functional domains (1.17%) and non-A-B interdomains (1.11%). Multivariate logistic regression analysis showed that lower HBV-DNA load (<10⁶ IU/mL) was an independent factor associated with potential DR-related variants in untreated CHB patients (P <0.05).

Conclusion

Potential DR-related variants were frequent and complex in untreated Chinese CHB patients. Furthermore, the variants may contribute to decreased serum HBV-DNA loads. However, the effects of potential DR-related variants on the antiviral therapy and liver disease progression require further study.

Molecular and serological characterization of occult hepatitis B virus infection among patients with hemophilia

BACKGROUND

The occult hepatitis B virus infection (OBI) is a health concern among high-risk groups and immunosuppressed individuals. There is still a paucity of data regarding the occult hepatitis B virus infection among hemophilic patients. With this in mind, we aimed to evaluate the molecular prevalence of OBI among clients with hemophilia.

METHODS

Totally, 87 hemophilic patients were selected to be studied. To detect OBI, nested polymerase chain reaction test was used to amplify HBV-S, X, and Core regions. Viral load was determined using an in-house real-time PCR assay. Finally, sequence of S gene was used for genotyping and analysis of mutations.

RESULTS

The mean age of patients was 28.4 ± 5.3 years old, with 90.7% of whom were men. HBV-DNA was detected in eight subjects (9.3%). The rate of OBI was much higher in anti-HBs seronegative subjects than that in other patients (P = 0.019). All OBI cases had HBV genotype D, subgenotype D1. In addition, five out of eight cases (62.5%) showed detectable viral loads (a mean viral load of 4.5 × 10 ² copies/mL). sR73H, sI110L, sP120A, sP127T, sQ129H, sG130R, and sC137S were shown to be the most determinant escape mutation and OBI-relevant mutants.

CONCLUSION

The rate of OBI among the studied population of hemophilia seems to be remarkable. Therefore, screening for OBI must be a routine practice in patients with hemophilia and also patients undergoing immunosuppressive treatments. Amino acid substitutions were observed in the major hydrophilic region. However further investigations are needed for analysis of exact function.

Molecular Characterization of Near Full-Length Genomes of Hepatitis B Virus Isolated from Predominantly HIV Infected Individuals in Botswana

The World Health Organization plans to eliminate hepatitis B and C Infections by 2030. Therefore, there is a need to study and understand hepatitis B virus (HBV) epidemiology and viral evolution further, including evaluating occult (HBsAg-negative) HBV infection (OBI), given that such infections are frequently undiagnosed and rarely treated. We aimed to molecularly characterize HBV genomes from 108 individuals co-infected with human immunodeficiency virus (HIV) and chronic hepatitis B (CHB) or OBI identified from previous HIV studies conducted in Botswana from 2009 to 2012. Full-length (3.2 kb) and nearly full-length (~3 kb) genomes were amplified by nested polymerase chain reaction (PCR). Sequences from OBI participants were compared to sequences from CHB participants and GenBank references to identify OBI-unique mutations. HBV genomes from 50 (25 CHB and 25 OBI) individuals were successfully genotyped. Among OBI participants, subgenotype A1 was identified in 12 (48%), D3 in 12 (48%), and E in 1 (4%). A similar genotype distribution was observed in CHB participants. Whole HBV genome sequences from Botswana, representing OBI and CHB, were compared for the first time. There were 43 OBI-unique mutations, of which 26 were novel. Future studies using larger sample sizes and functional analysis of OBI-unique mutations are warranted.

Naturally occurring hepatitis B virus reverse transcriptase mutations related to potential antiviral drug resistance and liver disease progression

The annual number of deaths caused by hepatitis B virus (HBV)-related disease, including cirrhosis and hepatocellular carcinoma (HCC), is estimated as 887000. The reported prevalence of HBV reverse transcriptase (RT) mutation prior to treatment is varied and the impact of preexisting mutations on the treatment of naïve patients remains controversial, and primarily depends on geographic factors, HBV genotypes, HBeAg serostatus, HBV viral loads, disease progression, intergenotypic recombination and co-infection with HIV. Different sensitivity of detection methodology used could also affect their prevalence results. Several genotype-dependent HBV RT positions that can affect the emergence of drug resistance have also been reported. Eight mutations in RT (rtL80I, rtD134N, rtN139K/T/H, rtY141F, rtM204I/V, rtF221Y, rtI224V, and rtM309K) are significantly associated with HCC progression. HBeAg-negative status, low viral load, and genotype C infection are significantly related to a higher frequency and prevalence of preexisting RT mutations. Preexisting mutations are most frequently found in the A-B interdomain of RT which overlaps with the HBsAg “a” determinant region, mutations of which can lead to simultaneous viral immune escape. In conclusion, the presence of baseline RT mutations can affect drug treatment outcomes and disease progression in HBV-infected populations via modulation of viral fitness and host-immune responses.

Deep sequencing analysis of quasispecies in the HBV pre-S region and its association with hepatocellular carcinoma

BACKGROUND

The association between the evolution of hepatitis B virus (HBV) quasispecies and the development of hepatocellular carcinoma (HCC) is unknown.

METHODS

We used deep sequencing to examine the dynamics of HBV quasispecies and their relationship to HCC development. Thirty-two chronic hepatitis B (CHB) patients with HCC (HCC group) and 32 matched CHB patients without HCC (controls) were recruited. Fourteen patients from each group had serial sera available up to 9 years before the time of the present study. Deep sequencing of the HBV pre-S regions was performed. HBV quasispecies complexity, diversity, and intrapatient prevalence of pre-S deletions/mutations were analyzed.

RESULTS

Compared with control patients, HCC patients had a significant greater quasispecies complexity (p = 0.04 at the nucleotide level), greater diversity (p = 0.004 and 0.009 at the nucleotide level and the amino acid level respectively), and a trend of greater complexity at the amino acid level (p = 0.065). HCC patients had a higher intrapatient prevalence of pre-S deletions and point mutations (at codons 4, 27, and 167) compared with the control patients (all p < 0.05). Longitudinal observation in the sera of 14 HCC patients showed that quasispecies complexity (p = 0.027 and 0.024 at the nucleotide level and the amino acid level respectively) and diversity (p = 0.035 and 0.031 at the nucleotide level and the amino acid level respectively) increased as the disease progressed to HCC.

CONCLUSIONS

Increased HBV quasispecies complexity and diversity in the pre-S region, probably reflecting enhanced virus-host interplay, was associated with disease progression from CHB to HCC.

Selection of the highly replicative and partially multidrug resistant rtS78T HBV polymerase mutation during TDF-ETV combination therapy

BACKGROUND & AIMS

Patients chronically infected with the hepatitis B virus (HBV) and receiving long-term treatment with nucleoside or nucleotide analogues are at risk of selecting HBV strains with complex mutational patterns. We herein report two cases of HBV-infected patients with insufficient viral suppression, despite dual antiviral therapy with entecavir (ETV) and tenofovir (TDF). One patient died from aggressive hepatocellular carcinoma (HCC).

METHODS

Serum samples from the two patients at different time points were analyzed using ultra-deep pyrosequencing analysis. HBV mutations were identified and transiently transfected into hepatoma cells in vitro using replication-competent HBV vectors, and functionally analyzed. We assessed replication efficacy, resistance to antivirals and potential impact on HBV secretion (viral particles, exosomes).

RESULTS

Sequencing analyses revealed the selection of the rtS78T HBV polymerase mutation in both cases that simultaneously creates a premature stop codon at sC69 and thereby deletes almost the entire small HBV surface protein. One of the patients had an additional 261bp deletion in the preS1/S2 region. Functional analyses of the mutations in vitro revealed that the rtS78T/sC69∗ mutation, but not the preS1/S2 deletion, significantly enhanced viral replication and conferred reduced susceptibility to ETV and TDF. The sC69∗ mutation caused truncation of HBs protein, leading to impaired detection by commercial HBsAg assay, without causing intracellular HBsAg retention or affecting HBV secretion.

CONCLUSIONS

The rtS78T/sC69∗ HBV mutation, associated with enhanced replication and insufficient response to antiviral treatment, may favor long-term persistence of these isolates. In addition to the increased production of HBV transcripts and the sustained secretion of viral particles in the absence of antigenic domains of S protein, this HBV mutation may predispose patients to carcinogenic effects.

LAY SUMMARY

Long-term treatment with antiviral drugs carries the risk of selecting mutations in the hepatitis B virus (HBV). We herein report two cases of patients with insufficient response to dual tenofovir and entecavir therapy. Molecular analyses identified a distinct mutation, rtS78T/sC69∗, that abolishes HBsAg detection, enhances replication, sustains exosome-mediated virion secretion and decreases susceptibility to antivirals, thereby representing a potentially high-risk mutation for HBV-infected individuals.

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.

Dynamics of Genotypic Mutations of the Hepatitis B Virus Associated With Long-Term Entecavir Treatment Determined With Ultradeep Pyrosequencing: A Retrospective Observational Study

The aim of the study is to explore the evolution of genotypic mutations within the reverse transcriptase region in partial virological responders (PVRs) receiving long-term entecavir (ETV) treatment. A total of 32 patients were classified as completely virological responders (CVRs) (n = 12) or PVRs (n = 20). Five partial responders were hepatitis B virus (HBV)-DNA positive after long-term therapy, which lasted for >3 years. A total of 71 serum samples from these 32 patients were assayed by ultra-deep pyrosequencing (UDPS): 32 samples were from all patients at baseline, and 39 were from PVRs with sequential inter-treatment. Approximately 84,708 sequences were generated per sample. At baseline, the quasispecies heterogeneity did not significantly differ between the 2 groups. The frequencies of substitutions indicating pre-existence of nucleos(t)ide analog resistant (NAr) mutants ranged from 0.10% to 6.70%, which did not statistically differ between groups either. However, the substitutions associated with the NAr mutants were significantly different from those associated with the non-NAr mutants in 13 patients; 6 of these patients were PVRs and the others were CVRs. Five patients were HBV DNA positive after regular ETV monotherapy for >3 years, and 4 of these patients underwent mild NAr substitution fluctuations (<20%). One patient developed virological breakthrough while bearing single, double, and triple (rtL180 M, rtM204 V, rtS202G) substitutions. In addition to the common substitutions, unknown amino acid substitutions, such as rtL145 M/S, rtF151Y/L, rtR153Q, rtI224 V, rtN248H, rtS223A, rtS256C, need to be further verified. NAr substitutions are observed at frequencies of 0.10% to 6.7% before therapy. Long-term ETV therapy generally results in virological responses, as long as the proportion of resistance mutations remains at a relatively low level. Genotypic resistance to ETV is detected in all PVRs receiving long-term ETV therapy.

Implementation of Next-Generation Sequencing for Hepatitis B Virus Resistance Testing and Genotyping in a Clinical Microbiology Laboratory

Sanger sequencing or DNA hybridization have been the primary modalities for hepatitis B (HBV) resistance testing and genotyping; however, there are limitations, such as low sensitivity and the inability to detect novel mutations. Next-generation sequencing (NGS) for HBV can overcome these limitations, but there is limited guidance for clinical microbiology laboratories to validate this novel technology. In this study, we describe an approach to implementing deep pyrosequencing for HBV resistance testing and genotyping in a clinical virology laboratory. A nested PCR targeting the pol region of HBV (codons 143 to 281) was developed, and the PCR product was sequenced by the 454 Junior (Roche). Interpretation was performed by ABL TherapyEdge based on European Association for the Study of the Liver (EASL) guidelines. Previously characterized HBV samples by INNO-LiPA (LiPA) were compared to NGS with discordant results arbitrated by Sanger sequencing. Genotyping of 105 distinct samples revealed a concordance of 95.2% (100/105), with Sanger sequencing confirming the NGS result. Resistance testing by NGS was concordant with LiPA in 85% (68/80) of previously characterized samples. Additional mutations were found in 8 samples, which related to the identification of low-level mutant subpopulations present at <10% (6/8). To balance the costs of testing for the validation study, reproducibility of the NGS was investigated through an analysis of sequence variants at loci not associated with resistance in a single patient sample. Our validation approach attempts to balance costs with efficient data acquisition.

Pre-existing mutations in reverse transcriptase of hepatitis B virus in treatment-naive Chinese patients with chronic hepatitis B

High rate of viral replication and lacking of proofreading activity in hepatitis B virus (HBV) polymerase lead to the generation of mutations in HBV virus. Mutations in the reverse transcriptase (RT) region of HBV polymerase are demonstrated to be strongly associated with drug resistance during antiviral treatment. However, the presence of mutations as well as its clinical significance in treatment-naïve hepatitis patients (defined as pre-existing mutations) need to be further investigated. In the present study, a total of 168 serum samples from treatment-naive chronic hepatitis B (CHB) patients were collected, and the RT region of HBV polymerase was sequenced. The results showed that pre-existing mutations in the RT region of HBV polymerase were detected in 43 of 168 (25.6%) treatment-naive CHB patients within which there were no well-characterized primary nucleotide analogs (NAs) resistance sites. Three dominant sites at rt191, rt207 and rt226 were found mutant in 7(16.28%), 8(18.60%), and 14(32.56%) samples respectively among these 43 patients. No significant correlation was found between pre-existing mutations and gender, age, HBV genotype, ALT, HBeAg or HBV DNA loads. However, patients with pre-existing RT mutations under HBeAg sero-negative status exhibited decreased HBV DNA loads, which contributed to the decreased HBV DNA loads in the total HBeAg sero-negative patients. The above investigation indicated that there was a prevalence of pre-existing mutations in RT region of HBV polymerase which might affect the serum HBV DNA level in treatment-naive CHB patients. Its effects on the occurrence of NAs resistance and the prognosis after treatment need to be further investigated.

Different mechanism of selection of adefovir-resistant mutant viruses during adefovir monotherapy in patients with lamivudine-resistant chronic hepatitis B

BACKGROUND

Adefovir (ADV) resistance is more frequent in lamivudine (LMV)-resistant chronic hepatitis B (CHB) patients than in nucleos(t)ide analogue-naïve patients. The majority of LMV-resistant mutants harbor the rtM204V/I mutation, while a minor fraction harbor the rtA181V/T mutation. We aimed to elucidate the mechanism of the high rate of ADV resistance in LMV-resistant patients during ADV therapy.

METHODS

We performed a clonal analysis of HBV reverse transcriptase in treatment-naïve (n = 3) and LMV-resistant patients before ADV therapy (n = 14). Dynamic changes in the viral population (n = 9) during ADV therapy were also analyzed.

RESULTS

Before ADV therapy, rtA181V/T was observed in 30 of 680 clones (4.4%) from 7 patients with LMV resistance under dominant rt204V/I mutation and in one of 150 clones in treatment-naïve patients. The rtA181V/T mutation was more frequently found in clones from LMV-resistant patients than in treatment-naïve patients (p = 0.029). The rtN236T mutation was not observed in any clone. During ADV therapy, most rtM204V/I mutants were replaced by wild type in all 3 patients without the rtA181V/T mutation and in one patient with the rtA181V/T mutation. Subsequently, wild type was replaced by the rtN236T and/or rtA181V/T mutant. In patients with the rtA181V/T mutation (n = 6), the rtA181V/T mutant overtook the rtM204V/I mutant in 3 of 4 patients with ADV resistance. In 2 patients without ADV resistance, most of the viral population was replaced by wild type by the last follow-up.

CONCLUSION

The high rate of ADV resistance in patients with LMV-resistance might be attributable to preexisting rtA181V/T mutant virus.