Clinical application of estimating hepatitis B virus quasispecies complexity by massive sequencing: correlation between natural evolution and on-treatment evolution

Inspecting the Ribozyme Region of Hepatitis Delta Virus Genotype 1: Conservation and Variability

The hepatitis delta virus (HDV) genome has an autocatalytic region called the ribozyme, which is essential for viral replication. The aim of this study was to use next-generation sequencing (NGS) to analyze the ribozyme quasispecies (QS) in order to study its evolution and identify highly conserved regions potentially suitable for a gene-silencing strategy. HDV RNA was extracted from 2 longitudinal samples of chronic HDV patients and the ribozyme (nucleotide, nt 688-771) was analyzed using NGS. QS conservation, variability and genetic distance were analyzed. Mutations were identified by aligning sequences with their specific genotype consensus. The main relevant mutations were tested in vitro. The ribozyme was conserved overall, with a hyper-conserved region between nt 715-745. No difference in QS was observed over time. The most variable region was between nt 739-769. Thirteen mutations were observed, with three showing a higher frequency: T23C, T69C and C64 deletion. This last strongly reduced HDV replication by more than 1 log in vitro. HDV Ribozyme QS was generally highly conserved and was maintained during follow-up. The most conserved portion may be a valuable target for a gene-silencing strategy. The presence of the C64 deletion may strongly impair viral replication, as it is a potential mechanism of viral persistence.

Analysis of entire hepatitis B virus genomes reveals reversion of mutations to wild type in natural infection, a 15 year follow-up study

It has been reported that some mutations in the genome of hepatitis B virus (HBV) may predict the outcome of the virus infection. However, evolutionary data derived from long-term longitudinal analysis of entire HBV genomes using next generation sequencing (NGS) remain rare. In this study, serum samples were collected from asymptomatic hepatitis B surface antigen (HBsAg) carriers from a long-term prospective cohort. The entire HBV genome was amplified by polymerase chain reaction (PCR) and sequenced using NGS. Twenty-eight time series serum samples from nine subjects were successfully analysed. The Shannon entropy (Sn) ranged from 0 to 0.89, with a median value of 0.76, and the genetic diversity (D) ranged from 0 to 0.013, with a median value of 0.004. Intrahost HBV viral evolutionary rates ranged from 2.39E-04 to 3.11E-03. Double mutations at nt1762(A → T) and 1764(G → A) and a stop mutation at nt1896(G → A) were seen in all sequences from subject BO129 in 2007. However, in 2019, most sequences were wild type at these positions. Deletions between nt 2920-3040 were seen in all sequences from subject TS115 in 2007 and 2013 but these were not present in 2004 or 2019. Some sequences from subject CC246 had predicted escape substitutions (T123N, G145R) in the surface protein in 2004, 2013 and 2019 but none of the sequences from 2007 had these changes. In conclusion, HBV mutations may revert to wild type in natural infection. Clinicians should be wary of predicting long-term prognoses on the basis of the presence of mutations.

Next-generation sequencing for the diagnosis of hepatitis B: current status and future prospects

INTRODUCTION

Hepatitis B virus (HBV) causes a complex and persistent infection with a major impact on patients health. Viral-genome sequencing can provide valuable information for characterizing virus genotype, infection dynamics and drug and vaccine resistance.

AREAS COVERED

This article reviews the current literature to describe the next-generation sequencing progress that facilitated a more comprehensive study of HBV quasispecies in diagnosis and clinical monitoring.

EXPERT OPINION

HBV variability plays a key role in liver disease progression and treatment efficacy. Second-generation sequencing improved the sensitivity for detecting and quantifying mutations, mixed genotypes and viral recombination. Third-generation sequencing enables the analysis of the entire HBV genome, although the high error rate limits its use in clinical practice.

The evolution and clinical impact of hepatitis B virus genome diversity

The global burden of hepatitis B virus (HBV) is enormous, with 257 million persons chronically infected, resulting in more than 880,000 deaths per year worldwide. HBV exists as nine different genotypes, which differ in disease progression, natural history and response to therapy. HBV is an ancient virus, with the latest reports greatly expanding the host range of the Hepadnaviridae (to include fish and reptiles) and casting new light on the origins and evolution of this viral family. Although there is an effective preventive vaccine, there is no cure for chronic hepatitis B, largely owing to the persistence of a viral minichromosome that is not targeted by current therapies. HBV persistence is also facilitated through aberrant host immune responses, possibly due to the diverse intra-host viral populations that can respond to host-mounted and therapeutic selection pressures. This Review summarizes current knowledge on the influence of HBV diversity on disease progression and treatment response and the potential effect on new HBV therapies in the pipeline. The mechanisms by which HBV diversity can occur both within the individual host and at a population level are also discussed.

Kinetics of hepatitis B surface antigen quasispecies during REP 2139-Ca therapy in HBeAg-positive chronic HBV infection

Chronic HBV infection results in various clinical manifestations due to different levels of immune response. In recent years, hepatitis B treatment has improved by long-term administration of nucleos(t)ide analogues (NUCs) and peg-interferon. Nucleic acid polymers (NAPs; REP 2139-Ca and REP 2139-Mg) are new antiviral drugs that block the assembly of subviral particles, thus preventing the release of HBsAg and allowing its clearance and restoration of functional control of infection when combined with various immunotherapies. In the REP 102 study (NCT02646189), 9 of 12 patients showed substantial reduction of HBsAg and seroconversion to anti-HBs in response to REP 2139-Ca, whereas 3 of 12 patients did not show responses (>1 log reduction of HBsAg and HBV DNA from baseline). We characterized the dynamic changes of HBV quasispecies (QS) within the major hydrophilic region (MHR) of the 'pre-S/S' open reading frame including the 'a' determinant in responders and nonresponders of the REP 102 study and four untreated matched controls. HBV QS complexity at baseline varied slightly between responders and nonresponders (P = .28). However, these responders showed significant decline in viral complexity (P = .001) as REP 2139-Ca therapy progressed but no significant change in complexity was observed among the nonresponders (P = .99). The MHR mutations were more frequently observed in responders than in nonresponders and matched controls. No mutations were observed in 'a' determinant of major QS population which may interfere with the detection of HBsAg by diagnostic assays. No specific mutations were found within the MHR which could explain patients' poor HBsAg response during REP 2139-Ca therapy.

Next generation sequencing identifies baseline viral mutants associated with treatment response to pegylated interferon in HBeAg-positive chronic hepatitis B

Current data of hepatitis B virus (HBV) variants associated with treatment outcome identified by next generation sequencing (NGS) are limited. This study was aimed at determining the role of baseline sequence variations in the enhancer II (EnhII), basal core promotor (BCP) and pre-core (PC) regions of HBV genotype C in patients treated with pegylated interferon (PEG-IFN). Patients with HBeAg-positive chronic hepatitis B (CHB) treated with 48-week PEG-IFN were enrolled. Combined response (CR) at week 96 was defined by HBeAg seroconversion plus HBV DNA < 2000 IU/mL and HBsAg < 1000 IU/mL. Pre-treatment viral mutations were characterized by Sanger sequencing and NGS (Miseq Illumina platform). Among 47 patients (32 male, mean age 32.4 years), CR was achieved in 12 (25.5%) individuals. Overall, NGS was superior to Sanger sequencing in detecting mutations (61.7% vs. 38.3%, P < 0.001). Based on NGS, the prevalence of T1753V (T1753C/A/G) and A1762T/G1764A variants were significantly lower in responders compared to non-responders (8.3% vs. 51.4%, P = 0.009 and 33.3% vs. 68.6%, P = 0.032, respectively). No significant difference between groups was found regarding C1653T and G1896A mutants. The absence of T1753V and A1762T/G1764A mutations were factors associated with CR (OR 11.65, 95%CI 1.36-100.16, P = 0.025, and OR 4.36, 95%CI 1.08-17.63, P = 0.039, respectively). The existence of pre-treatment T1753V, A1762T/G1764A mutations and their combination yielded negative predictive values of 94.7%, 85.7% and 93.8%, respectively. The presence of HBV mutants in the BCP region determined by NGS at baseline was associated with poor treatment outcome in patients with HBeAg-positive CHB receiving PEG-IFN.

Characterization of hepatitis B virus X gene quasispecies complexity in mono-infection and hepatitis delta virus superinfection

BACKGROUND

Hepatitis delta virus (HDV) seems to strongly suppress hepatitis B virus (HBV) replication, although little is known about the mechanism of this interaction. Both these viruses show a dynamic distribution of mutants, resulting in viral quasispecies. Next-generation sequencing is a viable approach for analyzing the composition of these mutant spectra. As the regulatory hepatitis B X protein (HBx) is essential for HBV replication, determination of HBV X gene (HBX) quasispecies complexity in HBV/HDV infection compared to HBV mono-infection may provide information on the interactions between these two viruses.

AIM

To compare HBV quasispecies complexity in the HBX 5’ region between chronic hepatitis delta (CHD) and chronic HBV mono-infected patients.

METHODS

Twenty-four untreated patients were included: 7/24 (29.2%) with HBeAg-negative chronic HBV infection (CI, previously termed inactive carriers), 8/24 (33.3%) with HBeAg-negative chronic hepatitis B (CHB) and 9/24 (37.5%) with CHD. A serum sample from each patient was first tested for HBV DNA levels. The HBX 5’ region [nucleotides (nt) 1255-1611] was then PCR-amplified for subsequent next-generation sequencing (MiSeq, Illumina, United States). HBV quasispecies complexity in the region analyzed was evaluated using incidence-based indices (number of haplotypes and number of mutations), abundance-based indices (Hill numbers of order 1 and 2), and functional indices (mutation frequency and nucleotide diversity). We also evaluated the pattern of nucleotide changes to investigate which of them could be the cause of the quasispecies complexity.

RESULTS

CHB patients showed higher median HBV-DNA levels [5.4 logIU/mL, interquartile range (IQR) 3.5-7.9] than CHD (3.4 logIU/mL, IQR 3-7.6) (P = n.s.) or CI (3.2 logIU/mL, IQR 2.3-3.5) (P < 0.01) patients. The incidence and abundance indices indicated that HBV quasispecies complexity was significantly greater in CI than CHB. A similar trend was observed in CHD patients, although only Hill numbers of order 2 showed statistically significant differences (CHB 2.81, IQR 1.11-4.57 vs CHD 8.87, 6.56-11.18, P = 0.038). There were no significant differences in the functional indices, but CI and CHD patients also showed a trend towards greater complexity than CHB. No differences were found for any HBV quasispecies complexity indices between CHD and CI patients. G-to-A and C-to-T nucleotide changes, characteristic of APOBEC3G, were higher in CHD and CI than in CHB in genotype A haplotypes, but not in genotype D. The proportion of nt G-to-A vs A-to-G changes and C-to-T vs T-to-C changes in genotype A and D haplotypes in CHD patients showed no significant differences. In CHB and CI the results of these comparisons were dependent on HBV genotype.

CONCLUSION

The lower-replication CHD and CI groups show a trend to higher quasispecies complexity than the higher-replication CHB group. The mechanisms associated with this greater complexity require elucidation.

Insights From Deep Sequencing of the HBV Genome-Unique, Tiny, and Misunderstood

Hepatitis B virus (HBV) is a unique, tiny, partially double-stranded, reverse-transcribing DNA virus with proteins encoded by multiple overlapping reading frames. The substitution rate is surprisingly high for a DNA virus, but lower than that of other reverse transcribing organisms. More than 260 million people worldwide have chronic HBV infection, which causes 0.8 million deaths a year. Because of the high burden of disease, international health agencies have set the goal of eliminating HBV infection by 2030. Nonetheless, the intriguing HBV genome has not been well characterized. We summarize data on the HBV genome structure and replication cycle, explain and quantify diversity within and among infected individuals, and discuss advances that can be offered by application of next-generation sequencing technology. In-depth HBV genome analyses could increase our understanding of disease pathogenesis and allow us to better predict patient outcomes, optimize treatment, and develop new therapeutics.

Hepatitis B virus: The challenge of an ancient virus with multiple faces and a remarkable replication strategy

The hepatitis B virus (HBV) is the prototype member of the Hepadnaviridae, an ancient family of hepatotropic DNA viruses, which may have originated from 360 to 430 million years ago and with evidence of endogenization in reptilian genomes >200 million years ago. The virus is currently estimated to infect more than 250 million humans. The extremely successful spread of this pathogen among the human population is explained by its multiple particulate forms, effective transmission strategies (particularly perinatal transmission), long induction period and low associated mortality. These characteristics confer selective advantages, enabling the virus to persist in small, disperse populations and spread worldwide, with high prevalence rates in many countries. The HBV replication strategy is remarkably complex and includes a multiplicity of particulate structures. In addition to the common virions containing DNA in a relaxed circular (rcDNA) or double-stranded linear (dslDNA) forms, the viral population includes virion-like particles containing RNA or "empty" (viral envelopes and capsids without genomes), subviral particles (only an envelope) and even naked capsids. Consequently, several forms of the genome coexist in a single infection: (i) the "traveler" forms found in serum, including rcDNA and dslDNA, which originate from retrotranscription of a messenger RNA (the pregenomic RNA, another form of the viral genome itself) and (ii) forms confined to the host cell nucleus, including covalently closed circular DNA (cccDNA), which leads to a minichromosome form associated with histones and viral proteins, and double-stranded DNA integrated into the host genome. This complex composition lends HBV a kind of "multiple personality". Are these additional particles and genomic forms simple intermediaries/artifacts or do they play a role in the viral life cycle?

Host and viral factors associated with serum hepatitis B virus RNA levels among patients in need for treatment

Hepatitis B virus (HBV) RNA in serum is a novel biomarker for intrahepatic HBV replication and treatment response. For its proper use, it is essential to identify factors influencing serum HBV RNA level. Using a rapid amplification of complimentary DNA (cDNA) ends (RACE) PCR technique (lower limit of detection [LLD], 800 copies/mL [c/mL]), serum HBV RNA levels were measured in samples of 488 untreated individuals with chronic HBV infection who were eligible to treatment according to currently used recommendations. We explored the association of serum levels of HBV RNA with patient- and virus-associated factors. HBV genotype distribution was 21/10/20/46/3% for A/B/C/D/other. Mean HBV RNA serum level was 5.9 (1.6) log₁₀ c/mL (hepatitis B e antigen [HBeAg]-positive chronic hepatitis B [CHB], 6.5 [1.2] log c/mL; HBeAg-negative CHB, 4.1 [1.2] log c/mL; P < 0.001). By multivariable linear regression, factors associated with lower HBV RNA level were HBeAg negativity (β = -0.69; P < 0.001), HBV genotypes A (β = -0.13; P = 0.002), B (β = -0.07; P = 0.049), and C (β = -0.61; P < 0.001) in comparison to D, and presence of HBV basal core promoter mutation either alone (β = -0.14; P = 0.001) or in combination with precore mutation (β = -0.22; P < 0.001). Higher serum alanine aminotransferase (ALT) was associated with higher HBV RNA (β = 0.23; P < 0.001). HBV RNA correlated strongly with HBV DNA (HBeAg-pos, r = 0.72; P < 0.001; HBeAg-neg, r = 0.78; P < 0.001) and moderately with quantitative hepatitis B surface antigen (qHBsAg; HBeAg-pos, r = 0.54; P < 0.001; HBeAg-neg, r = 0.19; P = 0.04) and quantitative hepatitis B surface antigen (qHBeAg; r = 0.41; P < 0.001).

CONCLUSION

In this multiethnic cohort of 488 untreated individuals with CHB, factors associated with serum HBV RNA level were HBeAg status, serum ALT, HBV genotype, and presence of basal core promotor mutations. For the future use of serum HBV RNA as a clinical marker, it seems mandatory to take these factors into consideration. (Hepatology 2018).

Deep sequencing identifies hepatitis B virus core protein signatures in chronic hepatitis B patients

BACKGROUND

We aimed to identify HBc amino acid differences between subgroups of chronic hepatitis B (CHB) patients.

METHODS

Deep sequencing of HBc was performed in samples of 89 CHB patients (42 HBeAg positive, 47 HBeAg negative). Amino acid types were compared using Sequence Harmony to identify subgroup specific sites between HBeAg-positive and -negative patients, and between patients with combined response and non-response to peginterferon/adefovir combination therapy.

RESULTS

We identified 54 positions in HBc where the frequency of appearing amino acids was significantly different between HBeAg-positive and -negative patients. In HBeAg negative patients, 22 positions in HBc were identified which differed between patients with treatment response and those with non-response. The fraction non-consensus sequence on selected positions was significantly higher in HBeAg-negative patients, and was negatively correlated with HBV DNA and HBsAg levels.

CONCLUSIONS

Sequence Harmony identified a number of amino acid changes associated with HBeAg-status and response to peginterferon/adefovir combination therapy.

Detection of hyper-conserved regions in hepatitis B virus X gene potentially useful for gene therapy

AIM

To detect hyper-conserved regions in the hepatitis B virus (HBV) X gene (HBX) 5’ region that could be candidates for gene therapy.

METHODS

The study included 27 chronic hepatitis B treatment-naive patients in various clinical stages (from chronic infection to cirrhosis and hepatocellular carcinoma, both HBeAg-negative and HBeAg-positive), and infected with HBV genotypes A-F and H. In a serum sample from each patient with viremia > 3.5 log IU/mL, the HBX 5’ end region [nucleotide (nt) 1255-1611] was PCR-amplified and submitted to next-generation sequencing (NGS). We assessed genotype variants by phylogenetic analysis, and evaluated conservation of this region by calculating the information content of each nucleotide position in a multiple alignment of all unique sequences (haplotypes) obtained by NGS. Conservation at the HBx protein amino acid (aa) level was also analyzed.

RESULTS

NGS yielded 1333069 sequences from the 27 samples, with a median of 4578 sequences/sample (2487-9279, IQR 2817). In 14/27 patients (51.8%), phylogenetic analysis of viral nucleotide haplotypes showed a complex mixture of genotypic variants. Analysis of the information content in the haplotype multiple alignments detected 2 hyper-conserved nucleotide regions, one in the HBX upstream non-coding region (nt 1255-1286) and the other in the 5’ end coding region (nt 1519-1603). This last region coded for a conserved amino acid region (aa 63-76) that partially overlaps a Kunitz-like domain.

CONCLUSION

Two hyper-conserved regions detected in the HBX 5’ end may be of value for targeted gene therapy, regardless of the patients’ clinical stage or HBV genotype.

Quantitative characterization of hepatitis delta virus genome edition by next-generation sequencing

AIM

To determine the capacity of next-generation sequencing (NGS) for quantifying edited and unedited HDV populations, and to confirm if edition is a general phenomenon taking place along the entire HDV region analyzed, as we previously reported (Homs M et al. PLoS One 2016, 11, e0158557).

METHODS

Four serum samples from 4 patients with chronic HDV/HBV infection were included in the study. The region selected for analysis covered 360 nucleotides (nt), positions 910-1270 of the HDV genome, which included the HDAg ORF editing site (nt 1014 within codon 196). Quantification of edited and unedited genomes was performed by molecular cloning and Sanger sequencing and by NGS. To evaluate the reliability of the NGS values obtained, we combined a clone with an edited codon and one with an unedited codon in known percentages in a series of artificial mixtures, which were then analyzed by NGS. In addition, we determined the nt changes occurring over the complete amplified region after excluding the editing codon (196) to evaluate edition along it.

RESULTS

In total, 11,208 quality-filtered sequences were obtained in the 4 samples. The 95% confidence intervals for the proportions of unedited populations by molecular cloning and NGS were overlapping, and those of cloning were wider, indicating that they are comparable and that NGS is more precise than cloning. Unedited genomes predominated over edited ones in all 4 samples analyzed by NGS and in 3 of the 4 samples analyzed by molecular cloning. In total, 83,276 quality-filtered sequences were obtained from the artificial mixtures. Percentages of the two viral populations detected by NGS in these mixtures were comparable to the expected percentages. Evaluation of edition along the HDV coding region showed that transitions were more frequent than transversions, accounting for 63.09% and 36.91%, respectively. Interestingly, among the 4 possible transition-type changes, G:A and A:G accounted for 73.86% of the total.

CONCLUSION

Next-generation sequencing proved useful to quantify edited and unedited HDV genomes, and provided relevant information on the HDV quasispecies.

Associations between serum HBX quasispecies and their integration in hepatocellular carcinoma

HBV quasispecies are closely related to the course and outcome of liver disease. However, whether the complexity and diversity of HBX quasispecies affects its integration in the liver cell and thereby enhances the resultant carcinogenesis is still not clear. 15 HCC patients were recruited; genomic DNA and HBV DNA were extracted from liver cancer tissue and serum respectively. The integrated HBX fragment in liver cancer tissue was amplified by Alu repeat sequence-polymerase chain reaction (Alu-PCR) and sequenced. The serum HBX gene was amplified by nested PCR and sequenced. Quasispecies complexity and diversity, phylogenetic characteristics, lymphocyte count and survival time between HBX-integrated and HBX-unintegrated patients were evaluated. Results showed that the integrated HBX fragment was detected in the tumor tissue of nine patients, and the integration rate was 60.00% (9/15). Compared with the HBX-unintegrated patients, the HBX-integrated patients had a higher quasispecies complexity (P=0.028 and 0.004, at the nucleotide and amino acid levels, respectively). The HBX-integrated patients had a tendency of higher quasispecies diversity, lower lymphocyte count and the survival time. A total of 12 mutation sites were revealed in the HBX-integrated fragment after alignment with the reference sequence. In these, the HBX-integrated groups had significantly higher mutation frequencies at C1497T, A1630G, G1721A, A1762T/G1764A and A1774G. This study revealed influence factors of HBX integration both in virus and the host. The increased complexity and diversity of HBX quasispecies might destroy the host immune balance, and lead to HBX integration ultimately.

Deep sequencing shows that HBV basal core promoter and precore variants reduce the likelihood of HBsAg loss following tenofovir disoproxil fumarate therapy in HBeAg-positive chronic hepatitis B

OBJECTIVE

Hepatitis B e antigen (HBeAg) seroconversion and hepatitis B surface antigen (HBsAg) loss are important clinical outcomes for patients with chronic hepatitis B (CHB) treated with antiviral therapy. To date, there have been few studies that have evaluated viral sequence markers predicting serological response to nucleos(t)ide analogue (NA) treatment.

DESIGN

We used next-generation sequencing (NGS) and quantitative HBV serology (HBeAg and HBsAg) to identify viral sequence markers associated with serological response to long-term tenofovir disoproxil fumarate therapy among HBeAg-positive patients. In the GS-US-174-0103 study, approximately half the patients seroconverted to anti-HBe by week 192 and 11% of patients exhibited HBsAg loss, the closest outcome to functional cure. The frequency of HBV variants that have previously been associated with HBV clinical outcomes was evaluated. HBV viral diversity in baseline sequences generated by NGS was calculated using Shannon entropy.

RESULTS

NGS analysis of HBV sequences from 157 patients infected with genotypes A to D showed the frequency of variants in the basal core promoter (BCP) and precore (PC) regions varied by genotype and that these mutations were associated with the absence of HBsAg loss. This was the case even when mutations were present at frequencies below the threshold of detection by population sequencing. Increased viral diversity across the HBV genome as determined by NGS was also associated with reduced likelihood of HBsAg loss.

CONCLUSION

Patients with detectable BCP and/or PC variants and higher viral diversity have a lower probability of HBsAg loss during long-term NA therapy. Strategies to achieve functional cure of HBV infection through combination therapy should consider using NGS to stratify patients according to BCP/PC sequence. Consideration should also be given to earlier initiation of therapy prior to the emergence of BCP/PC variants.

TRIAL REGISTRATION NUMBER

NCT00116805; Post result.

Quantifying perinatal transmission of Hepatitis B viral quasispecies by tag linkage deep sequencing

Despite full immunoprophylaxis, mother-to-child transmission (MTCT) of Hepatitis B Virus still occurs in approximately 2–5% of HBsAg positive mothers. Little is known about the bottleneck of HBV transmission and the evolution of viral quasispecies in the context of MTCT. Here we adopted a newly developed tag linkage deep sequencing method and analyzed the quasispecies of four MTCT pairs that broke through immunoprophylaxis. By assigning unique tags to individual viral sequences, we accurately reconstructed HBV haplotypes in a region of 836 bp, which contains the major immune epitopes and drug resistance mutations. The detection limit of minor viral haplotypes reached 0.1% for individual patient sample. Dominance of “a determinant” polymorphisms were observed in two children, which pre-existed as minor quasispecies in maternal samples. In all four pairs of MTCT samples, we consistently observed a significant overlap of viral haplotypes shared between mother and child. We also demonstrate that the data can be potentially useful to estimate the bottleneck effect during HBV MTCT, which provides information to optimize treatment for reducing the frequency of MTCT.

Viral minority variants in the core promoter and precore region identified by deep sequencing are associated with response to peginterferon and adefovir in HBeAg negative chronic hepatitis B patients

BACKGROUND AND AIM

Precore (PC) and basal core promoter (BCP) mutations are associated with responses to interferon-based treatment in HBeAg-positive chronic hepatitis B (CHB) patients. Here, we identify viral minority variants in these regions and assess association with response to peginterferon-alfa (Peg-IFN) and adefovir combination therapy.

PATIENTS AND METHODS

Ultra-deep pyrosequencing analysis of the BCP and PC region was performed for 89 CHB patients (42 HBeAg-positive; 47 HBeAg-negative), at baseline and during treatment. Specifically, associations of individual positions with the HBeAg-negative phenotype were studied, as well as the association of the most prevalent mutations with combined response in HBeAg-positive and -negative patients at week 72 (HBeAg negativity, HBV-DNA <2000 IU/mL and ALT normalization at 24 weeks of treatment-free follow-up).

RESULTS

The mutations most strongly correlated with the HBeAg-negative phenotype were at positions 1762/1764 and 1896/1899 in the BCP and PC region, respectively. No major changes in nucleotide composition of these positions were observed during treatment. In HBeAg-negative patients, a combined presence of 1764A and 1896A was correlated with lower ALT levels (p = 0.004), whereas the presence of 1899A was correlated with higher age (p = 0.030), lower HBV-DNA level (p = 0.036), and previous IFN therapy (p = 0.032). The presence of 1764A/1896A or the absence of 1899A at baseline, was associated with lower response rates, after adjustment for HBV genotype (p = 0.031 and p = 0.017) and HBsAg level (p = 0.035 and p = 0.022).

CONCLUSION

We identified novel correlations between common BCP and PC variants with response to Peg-IFN and adefovir in HBeAg-negative patients. Ultimately, this may guide the selection of those patients most likely to benefit from Peg-IFN-based treatment.

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.

Review of Laboratory Tests used in Monitoring Hepatitis B Response to Pegylated Interferon and Nucleos(t)ide Analog Therapy

There are only two currently approved classes of hepatitis B virus (HBV) antiviral agents, pegylated interferon (Peg-IFN), and nucleos(t)ide analogs (NAs) for chronic HBV infection. Although Peg-IFN is used for a finite 48-week duration and offers a greater chance of sustained off-treatment virological response, it is poorly tolerated and can only be offered to selected patients. The NAs are well tolerated but require prolonged therapy due to risk of relapse with treatment cessation. There is evolving data that novel virological assays (e.g., quantitative hepatitis B surface antigen, quantitative hepatitis B core antigen, quantitative antibody to core protein) in combination with hepatitis B genotype and more sensitive HBV DNA polymerase chain reaction (PCR) assays may be useful to predict response to IFN as well as off-treatment NA durability. Utilization of these clinical laboratory tests may be important given the development of novel anti-HBV therapies, hoping to achieve a cure for chronic hepatitis B infection.

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.

Evidence of an Exponential Decay Pattern of the Hepatitis Delta Virus Evolution Rate and Fluctuations in Quasispecies Complexity in Long-Term Studies of Chronic Delta Infection

Chronic HDV infection can cause a severe form of viral hepatitis for which there is no specific treatment. Characterization of the hepatitis B or C viral quasispecies has provided insight into treatment failure and disease recurrence following liver transplantation, has proven useful to understand hepatitis B e antigen seroconversion, and has helped to predict whether hepatitis C infection will resolve or become chronic. It is likely that characterization of the hepatitis delta virus (HDV) quasispecies will ultimately have similar value for the management of this infection. This study sought to determine the RNA evolution rates in serum of chronic hepatitis delta (CHD) treatment-naïve patients, using next-generation sequencing methods. The region selected for study encompassed nucleotide positions 910 to 1270 of the genome and included the amber/W codon. Amber/W is a substrate of the editing process by the ADAR1 host enzyme and is essential for encoding the 2 delta antigens (HDAg). The amber codon encodes the small (unedited) HDAg form and the W codon the large (edited) HDAg form. The evolution rate was analyzed taking into account the time elapsed between samples, the percentage of unedited and edited genomes, and the complexity of the viral population. The longitudinal studies included 29 sequential samples from CHD patients followed up for a mean of 11.5 years. In total, 121,116 sequences were analyzed. The HDV evolution rate ranged from 9.5x10-3 to 1.2x10-3 substitutions/site/year and showed a negative correlation with the time elapsed between samples (p<0.05). An accumulation of transition-type changes was found to be responsible for higher evolution rates. The percentages of unedited and edited genomes and the quasispecies complexity showed no relationships with the evolution rate, but the fluctuations in the percentages of genomes and in complexity suggest continuous adaptation of HDV to the host conditions.

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.

Detection of Low-Level Mixed-Population Drug Resistance in Mycobacterium tuberculosis Using High Fidelity Amplicon Sequencing

Undetected and untreated, low-levels of drug resistant (DR) subpopulations in clinical Mycobacterium tuberculosis (Mtb) infections may lead to development of DR-tuberculosis, potentially resulting in treatment failure. Current phenotypic DR susceptibility testing has a theoretical potential for 1% sensitivity, is not quantitative, and requires several weeks to complete. The use of "single molecule-overlapping reads" (SMOR) analysis with next generation DNA sequencing for determination of ultra-rare target alleles in complex mixtures provides increased sensitivity over standard DNA sequencing. Ligation free amplicon sequencing with SMOR analysis enables the detection of resistant allele subpopulations at ≥0.1% of the total Mtb population in near real-time analysis. We describe the method using standardized mixtures of DNA from resistant and susceptible Mtb isolates and the assay's performance for detecting ultra-rare DR subpopulations in DNA extracted directly from clinical sputum samples. SMOR analysis enables rapid near real-time detection and tracking of previously undetectable DR sub-populations in clinical samples allowing for the evaluation of the clinical relevance of low-level DR subpopulations. This will provide insights into interventions aimed at suppressing minor DR subpopulations before they become clinically significant.

Complex Genotype Mixtures Analyzed by Deep Sequencing in Two Different Regions of Hepatitis B Virus

This study assesses the presence and outcome of genotype mixtures in the polymerase/surface and X/preCore regions of the HBV genome in patients with chronic hepatitis B virus (HBV) infection. Thirty samples from ten chronic hepatitis B patients were included. The polymerase/surface and X/preCore regions were analyzed by deep sequencing (UDPS) in the first available sample at diagnosis, a pre-treatment sample, and a sample while under treatment. HBV genotype was determined by phylogenesis. Quasispecies complexity was evaluated by mutation frequency and nucleotide diversity. The polymerase/surface and X/preCore regions were validated for genotyping from 113 GenBank reference sequences. UDPS yielded a median of 10,960 sequences per sample (IQR 16,645) in the polymerase/surface region and 11,595 sequences per sample (IQR 14,682) in X/preCore. Genotype mixtures were more common in X/preCore (90%) than in polymerase/surface (30%) (p<0.001). On X/preCore genotyping, all samples were genotype A, whereas polymerase/surface yielded genotypes A (80%), D (16.7%), and F (3.3%) (p = 0.036). Genotype changes in polymerase/surface were observed in four patients during natural quasispecies dynamics and in two patients during treatment. There were no genotype changes in X/preCore. Quasispecies complexity was higher in X/preCore than in polymerase/surface (p = 0.004). The results provide evidence of genotype mixtures and differential genotype proportions in the polymerase/surface and X/preCore regions. The genotype dynamics in HBV infection and the different patterns of quasispecies complexity in the HBV genome suggest a new paradigm for HBV genotype classification.