Genotyping hepatitis B virus dual infections using population-based sequence data

An Oxford Nanopore Technology-Based Hepatitis B Virus Sequencing Protocol Suitable For Genomic Surveillance Within Clinical Diagnostic Settings

Chronic hepatitis B virus (HBV) infection remains a significant public health concern, particularly in Africa, where there is a substantial burden. HBV is an enveloped virus, with isolates being classified into ten phylogenetically distinct genotypes (A - J) determined based on full-genome sequence data or reverse hybridization-based diagnostic tests. In practice, limitations are noted in that diagnostic sequencing, generally using Sanger sequencing, tends to focus only on the S-gene, yielding little or no information on intra-patient HBV genetic diversity with very low-frequency variants and reverse hybridization detects only known genotype-specific mutations. To resolve these limitations, we developed an Oxford Nanopore Technology (ONT)-based HBV genotyping protocol suitable for clinical virology, yielding complete HBV genome sequences and extensive data on intra-patient HBV diversity. Specifically, the protocol involves tiling-based PCR amplification of HBV sequences, library preparation using the ONT Rapid Barcoding Kit, ONT GridION sequencing, genotyping using Genome Detective software, recombination analysis using jpHMM and RDP5 software, and drug resistance profiling using Geno2pheno software. We prove the utility of our protocol by efficiently generating and characterizing high-quality near full-length HBV genomes from 148 left-over diagnostic Hepatitis B patient samples obtained in the Western Cape province of South Africa, providing valuable insights into the genetic diversity and epidemiology of HBV in this region of the world.

Nanopore sequencing from extraction-free direct PCR of dried serum spots for portable hepatitis B virus drug-resistance typing

BACKGROUND

Effective drug regimens for the treatment of hepatitis B virus (HBV) infections are essential to achieve the World Health Organisation commitment to eliminate viral hepatitis by 2030. Lamivudine (3TC) is widely used in countries with high levels of chronic HBV, however resistance has been shown to occur in up to 50 % of individuals receiving continuous monotherapy for 4 years. Telbivudine (LdT) is now more commonly used in place of lamivudine but is ineffective against 3TC-resistant HBV. Genotyping and identification of resistanceassociated substitutions (RAS) is not practical in many locations.

OBJECTIVES

A novel assay was designed to enable HBV genotyping and characterisation of resistance mutations directly from serum samples stored on filter paper, using Sanger and MinION sequencing.

STUDY DESIGN

The assay was applied to a cohort of 30 samples stored on filter paper for several years with HBV viral loads ranging from 8.2 × 10⁸ to 635 IU/mL. A set of 6 high-titre samples were used in a proof-of-principle study using the MinION sequencer.

RESULTS

The assay allowed determination of HBV genotype and elucidation of RAS down to 600 IU/mL using a 550bp amplicon. Sequencing of a 1.2 kb amplicon using a MinION sequencer gave results consistent with Sanger sequencing and allowed the identification of minor populations of variants.

CONCLUSIONS

We present two approaches for reliable HBV sequencing and RAS identification using methods suitable for resource-limited environments. This is the first demonstration of extraction-free DNA sequencing direct from DSS using MinION and these workflows are adaptable to the investigation of other DNA viruses.

A meta-analysis on genetic variability of RT/HBsAg overlapping region of hepatitis B virus (HBV) isolates of Bangladesh

Background and aim

Hepatitis B caused by HBV is a serious public health hazard prevalent worldwide including Bangladesh. Few scattered molecular studies of HBV have been reported in Bangladesh. This study aimed to analyze the genetic variability of RT/HBsAg overlapping region of HBV isolates of Bangladesh and determination of correlation among the genotype/serotype and HBsAg escape and/or drug-resistant mutations.

Methods

A total of 97 complete HBsAg sequences of Bangladeshi HBV isolates from 2005 to 2017 from NCBI GenBank were extracted and analyzed using several HBV bioinformatics tools such as Geno2pheno-HBV, HBV Serotyper, HIV-Grade:HBV-Tool, and CLC sequence viewer.

Results

The prevalence of genotypes A, C, and D are 18, 46 and 35% which correspond to serotype adw, adr, and ayw, respectively. The prevalence of HBsAg escape mutations is 51% and most of which (62%) are found in the genotype D followed by 32% in genotype C and 6% in genotype A. Interestingly most (24/36) of the sequences of HBsAg escape mutations contained 128 V mutant which all belongs to only serotype ayw3 (Genotype D). Prevalence of drug-resistant mutations is ~ 11%, most of which are from genotype C (63.64%) and D (36.36%). Lamivudine resistant mutations were found in ~ 11% of sequences followed by Telbivudine 10% and Adefovir 3% where Tenofovir showed susceptibility to all 97 sequences. Moreover, 7 among of 97 sequences showed both HBsAg and drugs resistant mutations and none of them are found due to the same nucleotide substitutions.

Conclusion

There is a strong correlation among the genotype/serotype and HBsAg escape and/or drug-resistant mutations. This meta-analytical review will be helpful for genotype-serotype prediction by PCR-based diagnosis and development of vaccine and/or diagnostic kits, and the treatment against HBV infection in the future.

VirGenA: a reference-based assembler for variable viral genomes

Characterization of the within-host genetic diversity of viral pathogens is required for selection of effective treatment of some important viral infections, e.g. HIV, HBV and HCV. Despite the technical ability of detection, there are conflicting data regarding the clinical significance of low-frequency variants, partially because of the difficulty of their distinguishing from experimental artifacts. The issue of cross-contamination is relevant for all highly sensitive techniques, including deep sequencing: even trace contamination leads to a significant increase of false positives in identified SNVs. Determination of infections by multiple genotypes of some viruses, the incidence of which can be considerable, especially in risk groups, is also clinically significant in some cases. We developed a new viral reference-guided assembler, VirGenA, that can separate mixtures of strains of different intraspecies genetic groups (genotypes, subtypes, clades, etc.) and assemble a separate consensus sequence for each group in a mixture. It produced long assemblies for mixture components of extremely low frequencies (<1%) allowing detection of cross-contamination of samples by divergent genotypes. We tested VirGenA on both clinical and simulated data. On both types of data, VirGenA shows better or similar results than the existing de novo assemblers. Cross-platform implementation (including source code) is freely available at https://github.com/gFedonin/VirGenA/releases.

Genetic variability in coding regions of the surface antigen and reverse transcriptase domain of hepatitis B virus polymerase, Colombia, 2002-2014

Introduction: Despite the availability of an effective vaccine and treatment to reduce the viral load and progressive hepatocellular injury, approximately 240 million people worldwide are chronically infected with the hepatitis B virus (HBV). In Colombia, the circulation of different viral genotypes has been confirmed. Mutations in the genome have been associated to antiviral therapy resistance, viral escape to neutralizing antibodies, occult infection and progression to hepatocellular carcinoma. Objective: To identify the genotypes and the presence of mutations in the coding region of the surface (S) antigen and the reverse transcriptase (RT) domain of the polymerase of HBV obtained from serum samples for hepatitis B diagnosis received by the Instituto Nacional de Salud during the period 2002-2014. Materials and methods: A total of 495 serum samples with previous HBsAg reactive result were used for molecular detection. A fragment of 1,591 nucleotides was sequenced, and the corresponding phylogenetic analysis was performed. Results: We detected the viral genome of HBV in 66 samples and 28 were successfully sequenced. The phylogenetic analysis allowed the identification of subgenotypes F3 and A2. The L180M and M204V resistance mutations were simultaneously identified in one sample, while the I169L resistance mutation was identified in another one. A single escape mutation, P120Q, was identified in one more. Two samples showed a deletion of 105 nucleotides in the preS1-preS2 region. Conclusions: The circulation of genotypes/subgenotypes F3 and A2 of HBV in Colombia was corroborated, as well as the presence of some resistance and escape mutations. The present study constitutes a contribution to the molecular epidemiology of HBV in Colombia.

geno2pheno[ngs-freq]: a genotypic interpretation system for identifying viral drug resistance using next-generation sequencing data

Identifying resistance to antiretroviral drugs is crucial for ensuring the successful treatment of patients infected with viruses such as human immunodeficiency virus (HIV) or hepatitis C virus (HCV). In contrast to Sanger sequencing, next-generation sequencing (NGS) can detect resistance mutations in minority populations. Thus, genotypic resistance testing based on NGS data can offer novel, treatment-relevant insights. Since existing web services for analyzing resistance in NGS samples are subject to long processing times and follow strictly rules-based approaches, we developed geno2pheno[ngs-freq], a web service for rapidly identifying drug resistance in HIV-1 and HCV samples. By relying on frequency files that provide the read counts of nucleotides or codons along a viral genome, the time-intensive step of processing raw NGS data is eliminated. Once a frequency file has been uploaded, consensus sequences are generated for a set of user-defined prevalence cutoffs, such that the constructed sequences contain only those nucleotides whose codon prevalence exceeds a given cutoff. After locally aligning the sequences to a set of references, resistance is predicted using the well-established approaches of geno2pheno[resistance] and geno2pheno[hcv]. geno2pheno[ngs-freq] can assist clinical decision making by enabling users to explore resistance in viral populations with different abundances and is freely available at http://ngs.geno2pheno.org.

BATVI: Fast, sensitive and accurate detection of virus integrations

Background

The study of virus integrations in human genome is important since virus integrations were shown to be associated with diseases. In the literature, few methods have been proposed that predict virus integrations using next generation sequencing datasets. Although they work, they are slow and are not very sensitive.

Results and discussion

This paper introduces a new method BatVI to predict viral integrations. Our method uses a fast screening method to filter out chimeric reads containing possible viral integrations. Next, sensitive alignments of these candidate chimeric reads are called by BLAST. Chimeric reads that are co-localized in the human genome are clustered. Finally, by assembling the chimeric reads in each cluster, high confident virus integration sites are extracted.

Conclusion

We compared the performance of BatVI with existing methods VirusFinder and VirusSeq using both simulated and real-life datasets of liver cancer patients. BatVI ran an order of magnitude faster and was able to predict almost twice the number of true positives compared to other methods while maintaining a false positive rate less than 1%. For the liver cancer datasets, BatVI uncovered novel integrations to two important genes TERT and MLL4, which were missed by previous studies. Through gene expression data, we verified the correctness of these additional integrations.

BatVI can be downloaded from http://biogpu.ddns.comp.nus.edu.sg/~ksung/batvi/index.html.

Widespread hepatitis B virus genotype G (HBV-G) infection during the early years of the HIV epidemic in the Netherlands among men who have sex with men

Background

Hepatitis B virus (HBV) variants belong to different genotypes, A-J, whose worldwide distribution is linked with geography, probably because viral spread was associated with ancient human migrations. HBV genotype G (HBV-G) is an aberrant genotype with little sequence divergence, suggesting a recent origin. HBV-G is strongly associated with certain risk groups such as intravenous drug users (IDUs) and men who have sex with men (MSM), but hardly with geography. The origin and epidemiology of HBV-G remain unresolved, as is the disease association.

Methods

To estimate the prevalence and possible time of introduction of HBV-G into the MSM community in Amsterdam, the Netherlands, we have retrospectively analysed 226 blood serum samples from HBsAg positive MSM enrolled in the Amsterdam Cohort Studies (ACS) on HIV infection and AIDS dating from 1984 to 1999 using genotype-specific PCR assays.

Results

Of the 226 HBsAg-positive samples, 149 were HBV DNA positive. Of those, 104 were positive for HBV genotype A (HBV-A) and five for HBV-G, and 40 showed a dual infection with both HBV-A and HBV-G. Being HIV-infected was significantly associated with a reduced HBV DNA viral load in blood, but not with the prevalence of HBV-G. Early virus already contained stop codons in the precore region and a 36 bp insertion in the core gene which are the characteristics of HBV-G.

Conclusions

HBV-G was introduced before 1985 into the Amsterdam MSM community. Early isolates show very limited sequence variation, confirming a low evolutionary rate. HBV-G acquisition was independent of HIV infection, but being HIV-infected was significantly associated with a reduced HBV viral load in blood, indicating a beneficial effect of early HIV infection in controlling HBV replication.

Electronic supplementary material

The online version of this article (doi:10.1186/s12879-016-1599-7) contains supplementary material, which is available to authorized users.

Overview of hepatitis B viral replication and genetic variability

Chronic infection with hepatitis B virus (HBV) greatly increases the risk for liver cirrhosis and hepatocellular carcinoma (HCC). HBV isolates worldwide can be divided into ten genotypes. Moreover, the immune clearance phase selects for mutations in different parts of the viral genome. The outcome of HBV infection is shaped by the complex interplay of the mode of transmission, host genetic factors, viral genotype and adaptive mutations, as well as environmental factors. Core promoter mutations and mutations abolishing hepatitis B e antigen (HBeAg) expression have been implicated in acute liver failure, while genotypes B, C, subgenotype A1, core promoter mutations, preS deletions, C-terminal truncation of envelope proteins, and spliced pregenomic RNA are associated with HCC development. Our efforts to treat and prevent HBV infection are hampered by the emergence of drug resistant mutants and vaccine escape mutants. This paper provides an overview of the HBV life cycle, followed by review of HBV genotypes and mutants in terms of their biological properties and clinical significance.

Molecular epidemiology of hepatitis B virus infection in Switzerland: a retrospective cohort study

BACKGROUND

Chronic hepatitis B virus (HBV) infection affects up to 7% of the European population. Specific HBV genotypes are associated with rapid progression to end-stage liver disease and sub-optimal interferon treatment responses. Although the geographic distribution of HBV genotypes differs between regions, it has not been studied in Switzerland, which lies at the crossroads of Europe.

METHODS

In a retrospective analysis of 465 HBV samples collected between 2002 and 2013, we evaluated the HBV genotype distribution and phylogenetic determinants, as well as the prevalence of serological evidence of hepatitis delta, hepatitis C and HIV infections in Switzerland. Baseline characteristics of patients were compared across their region of origin using Fisher's exact test and ANOVA, and risk factors for HBeAg positivity were assessed using logistic regression.

RESULTS

The Swiss native population represented 15.7% of HBV-infected patients living in Switzerland. In the overall population, genotype D was most prevalent (58.3%), whereas genotype A (58.9%) was the predominant genotype among the Swiss native population. The prevalence of patients with anti-HDV antibodies was 4.4%. Patients of Swiss origin were most likely to be HBeAg-positive (38.1%). HBV genotypes of patients living in Switzerland but sharing the same original region of origin were consistent with their place of birth.

CONCLUSIONS

The molecular epidemiology of HBV infection in Switzerland is driven by migration patterns and not by the genotype distribution of the native population. The prevalence of positive anti-HDV antibodies in our cohort was very low.

HBV genotypes and response to tenofovir disoproxil fumarate in HIV/HBV-coinfected persons

Background

Hepatitis B virus (HBV) genotypes can influence treatment outcome in HBV-monoinfected and human immunodeficiency virus (HIV)/HBV-coinfected patients. Tenofovir disoproxil fumarate (TDF) plays a pivotal role in antiretroviral therapy (ART) of HIV/HBV-coinfected patients. The influence of HBV genotypes on the response to antiviral drugs, particularly TDF, is poorly understood.

Methods

HIV/HBV-co-infected participants with detectable HBV DNA prior to TDF therapy were selected from the Swiss HIV Cohort Study. HBV genotypes were identified and resistance testing was performed prior to antiviral therapy, and in patients with delayed treatment response (>6 months). The efficacy of TDF to suppress HBV (HBV DNA <20 IU/mL) and the influence of HBV genotypes were determined.

Results

143 HIV/HBV-coinfected participants with detectable HBV DNA were identified. The predominant HBV genotypes were A (82 patients, 57 %); and D (35 patients, 24 %); 20 patients (14 %) were infected with multiple genotypes (3 % A + D and 11 % A + G); and genotypes B, C and E were each present in two patients (1 %). TDF completely suppressed HBV DNA in 131 patients (92 %) within 6 months; and in 12 patients (8 %), HBV DNA suppression was delayed. No HBV resistance mutations to TDF were found in patients with delayed response, but all were infected with HBV genotype A (among these, 5 patients with genotype A + G), and all had previously been exposed to lamivudine.

Conclusion

In HIV/HBV-coinfected patients, infection with multiple HBV genotypes was more frequent than previously reported. The large majority of patients had an undetectable HBV viral load at six months of TDF-containing ART. In patients without viral suppression, no TDF-related resistance mutations were found. The role of specific genotypes and prior lamivudine treatment in the delayed response to TDF warrant further investigation.

Performance of HBsAg point-of-care tests for detection of diagnostic escape-variants in clinical samples

BACKGROUND

Hepatitis B viruses (HBV) harboring mutations in the a-determinant of the Hepatitis B surface antigen (HBsAg) are associated with reduced reactivity of HBsAg assays.

OBJECTIVES

To evaluate the sensitivity and specificity of three HBsAg point-of-care tests for the detection of HBsAg of viruses harboring HBsAg mutations.

STUDY DESIGN

A selection of 50 clinical plasma samples containing HBV with HBsAg mutations was used to evaluate the performance of three HBsAg point-of-care tests (Vikia(®), bioMérieux, Marcy-L'Étoile, France. Alere Determine HBsAg™, Iverness Biomedical Innovations, Köln, Germany. Quick Profile™, LumiQuick Diagnostics, California, USA) and compared to the ARCHITECT HBsAg Qualitative(®) assay (Abbott Laboratories, Sligo, Ireland).

RESULTS

The sensitivity of the point-of-care tests ranged from 98% to 100%. The only false-negative result occurred using the Quick Profile™ assay with a virus harboring a D144A mutation.

CONCLUSIONS

The evaluated point-of-care tests revealed an excellent sensitivity in detecting HBV samples harboring HBsAg mutations.

Re-appraisal of old and new diagnostic tools in the current management of chronic hepatitis B

Hepatitis B virus (HBV) is a very complex and intricate DNA structure associated with a particular genomic organization and replication cycle. However, many years of investigations allowed clarification of the real HBV natural history, through a deeper knowledge of the behavior of HBV antigens and viral structures. Several of the old diagnostic tools, such as HBV surface antigen (HBsAg) and HBV e antigen (HBeAg) determinations, gained prominence now, since the variation of both HBsAg and HBeAg plasma levels was shown to predict treatment response. In addition, the availability of more sensitive methods, such as HBV DNA detection by real-time PCR, has improved the current knowledge of the relationships between HBV replication levels and the natural history of the disease. It is now well established that some HBV genotypes are associated with a better response to treatment with pegylated interferon. Despite the widely accepted value of liver biopsy as a staging tool, transient elastography is being increasingly acknowledged as a non-invasive method to assess liver stiffness, chiefly for detection of advanced fibrosis. Current international guidelines for the management of chronic hepatitis B have provided several accurate biochemical and serological criteria for selecting patients for treatment, allowing a higher number of cases to be enrolled into antiviral therapy. This review describes the different serological markers used for the study of HBV and their clinical significance. It also deals with methods used for detection of genotypes and HBV DNA, emphasizing the effectiveness of such determinations for both patient selection and chronic hepatitis B therapy/monitoring.

Hepatitis B virus drug resistance tools: one sequence, two predictions

Drug resistance testing, genotype analysis, and the determination of immune and vaccine escape variants support personalized antiviral treatment for hepatitis B patients. As phenotypic drug resistance testing for hepatitis B virus (HBV) is especially labor-intensive, due to the lack of simple cell culture systems, genotypic methods play a very pronounced role. The genetic structure of HBV allows the simultaneous analysis of two different genes by examination of a single region in the genome of HBV. Nevertheless, the overlapping open reading frames of the hepatitis B surface antigen (HBsAg) and the reverse transcriptase (RT) have to be interpreted separately. In diagnostic procedures, standard Sanger type sequencing (mostly performed as a dye-dideoxynucleotide terminator system) is still the most commonly used method. This allows using established techniques for interpreting those types of genetic information. Besides reviewing the foundation of drug resistance interpretation for HBV, different interpretation systems, either commercially available (TRUGENE, Abbott, ViroScore) or free to use (geno2pheno[HBV], HIV-GRADE HBV tool), are compared in this article.

High prevalence of hepatitis B virus dual infection with genotypes A and G in HIV-1 infected men in Amsterdam, the Netherlands, during 2000-2011

Background

Hepatitis B virus (HBV) is divided into 8 definite (A-H) and 2 putative (I, J) genotypes that show a geographical distribution. HBV genotype G, however, is an aberrant genotype of unknown origin that demonstrates severe replication deficiencies and very little genetic variation. It is often found in co-infections with another HBV genotype and infection has been associated with certain risk groups such as intravenous drug users and men having sex with men (MSM). We aimed to estimate the prevalence of HBV-G in the Netherlands by analysing samples from HBV-positive patients visiting the Academic Medical Center in Amsterdam.

Methods

Ninety-six HBV-infected patients, genotyped as HBV-A or HBV-G infected, were retrieved from the clinical database. Blood plasma samples were analysed with a newly-developed real-time PCR assay that detects HBV-A and HBV-G. For three patients, the HBV plasma viral load (pVL) of both genotypes was followed longitudinally. In addition, three complete genomes of HBV-G were sequenced to determine their relationship to global HBV-G strains.

Results

Ten HBV-G infections were found in the selected Dutch patients. All concerned HIV-1 infected males with HBV-A co-infection. Dutch HBV-G strains were phylogenetically closely related to reference HBV-G strains.

Conclusions

In this study, HBV-G infection in the Netherlands is found exclusively in HIV-1 infected men as co-infection with HBV-A. A considerable percentage (37%) of men infected with HBV and HIV-1 are actually co- infected with two HBV genotypes.