Bioinformatic curation and alignment of genotyped hepatitis B virus (HBV) sequence data from the GenBank public database

Hepatitis B virus (HBV) viremia despite tenofovir disoproxil fumarate-containing antiretroviral therapy in persons with HBV/HIV coinfection

BACKGROUND

The goal of treatment of hepatitis B virus (HBV) and human immunodeficiency virus (HIV) coinfection is suppression of both viruses; yet incomplete HBV suppression on tenofovir (TFV) disoproxil fumarate (TDF)-based antiretroviral therapy (ART) is common. This study investigated TFV resistance-associated mutations (RAMs) in individuals with HBV/HIV coinfection with viremia on TDF/lamivudine (3TC)-containing ART.

METHODS

Samples from individuals with HBV DNA levels ≥20 IU/mL in a cross-sectional study of 138 persons with HBV/HIV coinfection in Ghana were analyzed in the present study. HBV was sequenced for RAM analysis. TFV-diphosphate (TFV-DP) concentration in peripheral blood mononuclear cells (PBMCs) was used to assess ART adherence level.

RESULTS

Nine of 138 participants (6.5 %) had detectable HBV DNA levels ≥20 IU/mL while on ART. Seven of the nine participants had TFV-DP concentrations commensurate with 7 doses per week, and six had suppressed HIV RNA. Phylogenetic analysis revealed that eight sequences were HBV genotype E, with one genotype E/A recombinant. Ten previously-reported TFV RAMs were present in the study samples; eight were wild-type for HBV genotype E. The non-genotype-E-wild-type point mutations M267L and K333Q were found in two and one patients, respectively. No 3TC RAMs were found.

CONCLUSION

HBV viremia despite high adherence to TDF/3TC-based ART may be associated with the presence of TFV RAMs. These findings highlight the need for enhanced resistance monitoring and further research to examine the clinical significance of reported TFV RAMs. Individuals with HBV/HIV coinfection and TFV resistance on TDF-based ART may need alternative treatment strategies.

The Complex Role of HBeAg and Its Precursors in the Pathway to Hepatocellular Carcinoma

Hepatitis B virus (HBV) is one of the seven known human oncogenic viruses and has adapted to coexist with a single host for prolonged periods, requiring continuous manipulation of immunity and cell fate decisions. The persistence of HBV infection is associated with the pathogenesis of hepatocellular carcinoma, and various HBV proteins have been implicated in promoting this persistence. The precursor of hepatitis e antigen (HBeAg), is translated from the precore/core region and is post-translationally modified to yield HBeAg, which is secreted in the serum. HBeAg is a non-particulate protein of HBV and can act as both a tolerogen and an immunogen. HBeAg can protect hepatocytes from apoptosis by interfering with host signalling pathways and acting as a decoy to the immune response. By evading the immune response and interfering with apoptosis, HBeAg has the potential to contribute to the hepatocarcinogenic potential of HBV. In particular, this review summarises the various signalling pathways through which HBeAg and its precursors can promote hepatocarcinogenesis via the various hallmarks of cancer.

HexSE: Simulating evolution in overlapping reading frames

Gene overlap occurs when two or more genes are encoded by the same nucleotides. This phenomenon is found in all taxonomic domains, but is particularly common in viruses, where it may provide a mechanism to increase the information content of compact genomes. The presence of overlapping reading frames (OvRFs) can skew estimates of selection based on the rates of non-synonymous and synonymous substitutions, since a substitution that is synonymous in one reading frame may be non-synonymous in another and vice versa. To understand the impact of OvRFs on molecular evolution, we implemented a versatile simulation model of nucleotide sequence evolution along a phylogeny with any distribution of open reading frames in linear or circular genomes. We use a custom data structure to track the substitution rates at every nucleotide site, which is determined by the stationary nucleotide frequencies, transition bias and the distribution of selection biases (dN/dS) in the respective reading frames. Our simulation model is implemented in the Python scripting language. All source code is released under the GNU General Public License version 3 and are available at https://github.com/PoonLab/HexSE.

Integrated hepatitis B virus DNA maintains surface antigen production during antiviral treatment

The focus of hepatitis B functional cure, defined as sustained loss of hepatitis B virus (HBV) surface antigen (HBsAg) and HBV DNA from blood, is on eliminating or silencing the intranuclear template for HBV replication, covalently closed circular DNA (cccDNA). However, HBsAg also derives from HBV DNA integrated into the host genome (iDNA). Little is known about the contribution of iDNA to circulating HBsAg with current therapeutics. We applied a multiplex droplet digital PCR assay to demonstrate that iDNA is responsible for maintaining HBsAg quantities in some individuals. Using paired bulk liver tissue from 16 HIV/HBV-coinfected persons on nucleos(t)ide analog (NUC) therapy, we demonstrate that people with larger HBsAg declines between biopsies derive HBsAg from cccDNA, whereas people with stable HBsAg levels derive predominantly from iDNA. We applied our assay to individual hepatocytes in paired tissues from 3 people and demonstrated that the individual with significant HBsAg decline had a commensurate loss of infected cells with transcriptionally active cccDNA, while individuals without HBsAg decline had stable or increasing numbers of cells producing HBsAg from iDNA. We demonstrate that while NUC therapy may be effective at controlling cccDNA replication and transcription, innovative treatments are required to address iDNA transcription that sustains HBsAg production.

In silico analysis of mutations associated with occult hepatitis B virus (HBV) infection in South Africa

Occult hepatitis B virus (OBI) infection is defined by the presence of viral DNA in the liver and/or serum in absence of hepatitis B surface antigen (HBsAg). While multiple studies have identified mutations that are associated with OBI, only a small portion of these mutations have been functionally characterized in vitro. Using complementary in silico approaches, the effects of OBI-associated amino acid mutations on HBV protein function in HBV/HIV-positive ART-naïve South Africans were evaluated. Two OBI-associated mutations in the PreS1 region, one in the PreS2 region, and seven in the surface region of subgenotype A1 sequences were identified as deleterious. In subgenotype A2 sequences, 11 OBI-associated mutations in the PreS1 region, seven in the PreS2 region, and 31 in the surface region were identified as deleterious. In the polymerase region, 14 OBI-associated mutations in subgenotype A1 and 71 OBI-associated mutations in subgenotype A2 were identified as deleterious. This study utilized in silico approaches to characterize the likely impact of OBI-associated mutations on viral function, thereby identifying and prioritizing candidates and reducing the significant cost associated with functional studies that are essential for mechanistic studies of the OBI phenotype.

Identification of mutations in the S gene of hepatitis B virus in HIV positive Mexican patients with occult hepatitis B virus infection

INTRODUCTION AND AIM

Occult hepatitis B virus infection (OBI) is characterized by the presence of replication-competent hepatitis B virus (HBV) DNA in the liver and/or serum of patients with undetectable levels of the HBV surface antigen (HBsAg). Due to the shared infection routes HIV positive patients are at higher risk of developing OBI, thus, the aim of this study was to determine the frequency of OBI in Mexican HIV-infected patients and to identify mutations in the HBV S gene that could be associated to the development of OBI.

MATERIALS AND METHODS

Plasma samples from 50 HIV-infected patients with undetectable levels of the HBsAg were obtained and analyzed. The Core, PreS and S genes were amplified by nested PCR and sequenced by the Sanger method. To analyze HBV diversity in the OBI-positive patients, ten sequences of 762bp from the HBV S gene were selected, cloned, and subsequently sequenced for mutational analyses.

RESULTS

OBI infection was found with a frequency of 36% (18/50). All the HBV sequences corresponded to the H genotype. The most common mutations were: C19Y, Q129H, E164D, and I195M, with a frequency of 44%, 36%, 39% and 48% respectively.

CONCLUSIONS

In this study, we report the presence of OBI in a cohort of Mexican HIV-infected patients with an overall prevalence of 36%. Mutational analyses revealed that four non-silent mutations were frequent in different regions of the HBsAg gene, suggesting that they might be associated to the development of OBI in this population, nevertheless, further studies are required to determine their role in the pathogenesis of OBI.

The distribution of hepatitis B virus surface antigen polymorphisms at positions associated with vaccine escape

Hepatitis B virus (HBV) infects over 250 million people worldwide. Vaccination is effective at preventing infection, although several mutations within the "a" determinant region of the HBV surface antigen (HBsAg) are associated with vaccine escape. We evaluated the frequency, genotype, and global distribution of polymorphisms at sites associated with vaccine escape in 4244 unique full-length HBV genomes. The "a" determinant within the Surface gene was inspected for polymorphisms at sites identified previously associated with vaccine escape. Nearly, 268 (6.3%) sequences from 36 countries contained a polymorphism at a site associated with vaccine escape including 22 genotype A, 99 genotype B, 93 genotype C, 32 genotype D, 14 genotype E, 3 genotype F, 2 genotype G, and 3 genotype I. In genotype A, the most common polymorphism occurred at M133. In genotype B, Q129 and M133 occurred 45 and 51 times, respectively, accounting for 94% of polymorphisms. Polymorphisms at G145 were most frequent in genotype C, while P120 was most common in genotype D. Among all genotypes, polymorphisms at M133 were the most common and accounted for 30.9% of polymorphisms. Polymorphisms at T116, P120, F134, K141, and P142 occurred in geographically diverse locations, whereas polymorphisms at Q129, M133, D144, and G145 were concentrated in East Asia. While the sample size is large, this approach relied on convenience sampling within each country, and many countries have no data available, thereby highlighting the need for additional routine surveillance of surface antigen mutations associated with vaccine escape.

Global and regional dispersal patterns of hepatitis B virus genotype E from and in Africa: A full-genome molecular analysis

Description of the spatial characteristics of viral dispersal is important in understanding the history of infections. Nine hepatitis B virus (HBV) genotypes (A-I), and a putative 10th genotype (J), with distinct geographical distribution, are recognized. In sub-Saharan Africa (sub)-genotypes A1, D3 and E circulate, with E predominating in western Africa (WA), where HBV is hyperendemic. The low genetic diversity of genotype E (HBV/E) suggests its recent emergence. Our aim was to study the dispersal of HBV/E using full-length, non-redundant and non-recombinant sequences available in public databases. HBV/E was confirmed, and the phylogeny reconstruction performed using maximum likelihood (ML) with bootstrapping. Phylogeographic analysis was conducted by reconstruction of ancestral states using the criterion of parsimony on the estimated ML phylogeny. 46.5% of HBV/E sequences were found within monophyletic clusters. Country-wise analysis revealed the existence of 50 regional clusters. Sequences from WA were located close to the root of the tree, indicating this region as the most probable origin of the HBV/E epidemic and expanded to other geographical regions, within and outside of Africa. A localized dispersal was observed with sequences from Nigeria and Guinea as compared to other WA countries. Based on the sequences available in the databases, the phylogenetic results suggest that European strains originated primarily from WA whereas a majority of American strains originated in Western Central Africa. The differences in regional dispersal patterns of HBV/E suggest limited cross-border transmissions because of restricted population movements.

Cross-Protection of Hepatitis B Vaccination among Different Genotypes

Hepatitis B (HB) vaccination is the most effective method for preventing HB virus (HBV) infection. Universal HB vaccination containing recombinant HB surface antigens (HBsAg) is recommended. Our data revealed that human monoclonal HB surface antibody (anti-HBs) from individuals inoculated with genotype C-based HB vaccine induced cross-protection against HBV genotype A infection. An in vitro infection model demonstrated anti-HBs-positive sera from individuals inoculated with genotype A- or C-based HB vaccine harbored polyclonal anti-HBs that could bind to non-vaccinated genotype HBV. However, because there were low titers of anti-HBs specific for HBsAg of non-vaccinated genotype, high anti-HBs titers would be required to prevent non-vaccinated genotype HBV infection. Clinically, the 2015 Centers for Disease Control and Prevention guidelines state that periodic monitoring of anti-HBs levels after routine HB vaccination is not needed and that booster doses of HB vaccine are not recommended. However, the American Red Cross suggests that HB-vaccine-induced immune memory might be limited; although HB vaccination can prevent clinical liver injury (hepatitis), subclinical HBV infections of non-vaccinated genotypes resulting in detectable HB core antibody could not be completely prevented. Therefore, monitoring anti-HBs levels after routine vaccination might be necessary for certain subjects in high-risk groups.

In Silico Prediction of Human Leukocytes Antigen (HLA) Class II Binding Hepatitis B Virus (HBV) Peptides in Botswana

Hepatitis B virus (HBV) is the primary cause of liver-related malignancies worldwide, and there is no effective cure for chronic HBV infection (CHB) currently. Strong immunological responses induced by T cells are associated with HBV clearance during acute infection; however, the repertoire of epitopes (epi) presented by major histocompatibility complexes (MHCs) to elicit these responses in various African populations is not well understood. In silico approaches were used to map and investigate 15-mers HBV peptides restricted to 9 HLA class II alleles with high population coverage in Botswana. Sequences from 44 HBV genotype A and 48 genotype D surface genes (PreS/S) from Botswana were used. Of the 1819 epi bindings predicted, 20.2% were strong binders (SB), and none of the putative epi bind to all the 9 alleles suggesting that multi-epitope, genotype-based, population-based vaccines will be more effective against HBV infections as opposed to previously proposed broad potency epitope-vaccines which were assumed to work for all alleles. In total, there were 297 unique epi predicted from the 3 proteins and amongst, S regions had the highest number of epi (n = 186). Epitope-densities (Depi) between genotypes A and D were similar. A number of mutations that hindered HLA-peptide binding were observed. We also identified antigenic and genotype-specific peptides with characteristics that are well suited for the development of sensitive diagnostic kits. This study identified candidate peptides that can be used for developing multi-epitope vaccines and highly sensitive diagnostic kits against HBV infection in an African population. Our results suggest that viral variability may hinder HBV peptide-MHC binding, required to initiate a cascade of immunological responses against infection.

Analysis of genomic-length HBV sequences to determine genotype and subgenotype reference sequences

Hepatitis B virus (HBV) is a diverse, partially double-stranded DNA virus, with 9 genotypes (A-I), and a putative 10th genotype (J), characterized thus far. Given the broadening interest in HBV sequencing, there is an increasing requirement for a consistent, unified approach to HBV genotype and subgenotype classification. We set out to generate an updated resource of reference sequences using the diversity of all genomic-length HBV sequences available in public databases. We collated and aligned genomic-length HBV sequences from public databases and used maximum-likelihood phylogenetic analysis to identify genotype clusters. Within each genotype, we examined the phylogenetic support for currently defined subgenotypes, as well as identifying well-supported clades and deriving reference sequences for them. Based on the phylogenies generated, we present a comprehensive set of HBV reference sequences at the genotype and subgenotype level. All of the generated data, including the alignments, phylogenies and chosen reference sequences, are available online (https://doi.org/10.6084/m9.figshare.8851946) as a simple open-access resource.

Identification of hepatitis B virus genotype A/E recombinants in Ghana

Hepatitis B virus (HBV) exhibits a high degree of heterogeneity with at least 10 genotypes (A-J) identified to date. Intergenotypic recombination is relatively common. Previously, we investigated HBV drug resistance in HIV/HBV co-infected individuals in Ghana. After identifying multiple circulating genotypes and a novel D/E recombinant, we sought to determine if additional individuals were also infected with recombinant HBV. Partial genome sequences from three individuals were initially identified as genotype A4. Full-length HBV genomes were obtained using rolling circle amplification followed by PCR and shown to cluster with known A/E recombinant viruses. Similar recombination breakpoints were observed in these three individuals suggesting local spread of this novel recombinant HBV in Ghana.

Identification of hepatitis B virus genotype I in Thailand

The rare hepatitis B virus genotype I (HBV-I) classification includes complex A/G/C/U recombinants identified amongst the individuals from China, India, Laos, and Vietnam. Herein we report the first HBV-I specimen from Thailand, with detectable HBsAg despite a 10-amino-acid truncation. This HBV-I genome has a similar recombinant pattern to reference strains, including a C region that branches basal to references, suggesting a premodern era recombination event gave rise to HBV-I. With an average sequence divergence from other genotypes ranging from 7.66% (standard deviation [SD], 0.42%; C) to 14.27% (SD, 0.31%; H), this new genome supports the HBV-I classification as a unique genotype.

A novel hepatitis B virus recombinant genotype D4/E identified in a South African population

BACKGROUND

Genetic diversity is a characteristic trait of the hepatitis B virus (HBV) and has been associated with different clinical outcomes. In South Africa, HBV infection is a major public health concern. Most HBV infections are caused by genotype A strains. However rare cases of infection with HBV genotype D have been reported. The purpose of this study was to investigate the molecular characteristics of a rare HBV subgenotype D4 isolate.

METHODS

The full-length genome of isolate ZADGM6964 was amplified in a one-step polymerase chain reaction. The amplified product was purified and cloned into a pGEM®-T Easy Vector System to investigate the genetic diversity of the viral quasi-populations. The primary isolate and clones were then directly sequenced and analysed using an array of bioinformatics software.

RESULTS

Phylogenetic analysis showed that the primary isolate and cloned sequences formed a monophyletic cluster away from subgenotype D4 reference strains. Further recombination analysis revealed that isolate ZADGM6964 was in fact a D4/E recombinant strain with breakpoints identified within the X and overlapping pre-Core/Core open reading frames with a >70% bootstrap confidence level. The recombinant genotype D4/E was found to be unique from other D/E strains archived in the genetic database, GenBank.

CONCLUSION

This study represents the first ever report on the isolation and molecular characterization of an HBV D4/E recombinant strain in South Africa. The findings provide evidence of further HBV genetic diversity in South Africa than has been previously reported.

CCT: a coordinate conversion tool for hepatitis B virus

The hepatitis B virus (HBV) genome has ∼3 200 nucleotides, coding for seven proteins in four overlapping open reading frames (ORFs). Comparison of genomic coordinates between different samples and/or published literature requires manual conversion. An online tool is presented to convert nucleotide or amino acid positions between ORFs, regions and domains of the HBV genome. The user enters a position into an interactive web page, which then shows this position in all other applicable ORFs, regions or domains and plots it on a diagrammatic representation of the HBV genome. This tool assists researchers to convert coordinates, thereby facilitating comparisons between samples.

Immunomodulatory Function of HBeAg Related to Short-Sighted Evolution, Transmissibility, and Clinical Manifestation of Hepatitis B Virus

Hepatitis B virus (HBV) infection, a global public health problem can be asymptomatic, acute or chronic and can lead to serious consequences of infection, including cirrhosis, and hepatocellular carcinoma. HBV, a partially double stranded DNA virus, belongs to the family Hepadnaviridae, and replicates via reverse transcription of an RNA intermediate. This reverse transcription is catalyzed by a virus-encoded polymerase that lacks proof reading ability, which leads to sequence heterogeneity. HBV is classified into nine genotypes and at least 35 subgenotypes, which may be characterized by distinct geographical distributions. This HBV diversification and distinct geographical distribution has been proposed to be the result of the co-expansion of HBV with modern humans, after their out-of-Africa migration. HBeAg is a non-particulate protein of HBV that has immunomodulatory properties as a tolerogen that allows the virus to establish HBV infection in vivo. During the natural course of infection, there is seroconversion from a HBeAg-positive phase to a HBeAg-negative, anti-HBe-positive phase. During this seroconversion, there is loss of tolerance to infection and immune escape-HBeAg-negative mutants can be selected in response to the host immune response. The different genotypes and, in some cases, subgenotypes develop different mutations that can affect HBeAg expression at the transcriptional, translational and post-translational levels. The ability to develop mutations, affecting HBeAg expression, can influence the length of the HBeAg-positive phase, which is important in determining both the mode of transmission and the clinical course of HBV infection. Thus, the different genotypes/subgenotypes have evolved in such a way that they exhibit different modes of transmission and clinical manifestation of infection. Loss of HBeAg may be a sign of short-sighted evolution because there is loss of tolerogenic ability of HBeAg and HBeAg-negative virions are less transmissible. Depending on their ability to lead to HBeAg seroconversion, the genotype/subgenotypes exhibit varying degrees of short-sighted evolution. The “arms race” between HBV and the immune response to HBeAg is multifaceted and its elucidation intricate, with transmissibility and persistence being important for the survival of the virus. We attempt to shed some light on this complex interplay between host and virus.

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.

In Silico Analysis of Hepatitis B Virus Occult Associated Mutations in Botswana Using a Novel Algorithm

Occult hepatitis B infections (OBI) represent a reservoir of undiagnosed and untreated hepatitis B virus (HBV), hence the need to identify mutations that lead to this phenotype. Functionally characterizing these mutations by in vitro studies is time-consuming and expensive. To bridge this gap, in silico approaches, which predict the effect of amino acid (aa) variants on HBV protein function, are necessary. We developed an algorithm for determining the relevance of OBI-associated mutations using in silico approaches. A 3 kb fragment of subgenotypes A1 and D3 from 24 chronic HBV-infected (CHB) and 24 OBI participants was analyzed. To develop and validate the algorithm, the effects of 68 previously characterized occult-associated mutations were determined using three computational tools: PolyPhen2, SNAP2, and PROVEAN. The percentage of deleterious mutations (with impact on protein function) predicted were 52 (76.5%) by PolyPhen2, 55 (80.9%) by SNAP2, and 65 (95.6%) by PROVEAN. At least two tools correctly predicted 59 (86.8%) mutations as deleterious. To identify OBI-associated mutations exclusive to Botswana, study sequences were compared to CHB sequences from GenBank. Of the 43 OBI-associated mutations identified, 26 (60.5%) were predicted by at least two tools to have an impact on protein function. To our knowledge, this is the first study to use in silico approaches to determine the impact of OBI-associated mutations, thereby identifying potential candidates for functional analysis to facilitate mechanistic studies of the OBI phenotype.

A systematic review of hepatitis B virus (HBV) drug and vaccine escape mutations in Africa: A call for urgent action

International sustainable development goals for the elimination of viral hepatitis as a public health problem by 2030 highlight the pressing need to optimize strategies for prevention, diagnosis and treatment. Selected or transmitted resistance associated mutations (RAMs) and vaccine escape mutations (VEMs) in hepatitis B virus (HBV) may reduce the success of existing treatment and prevention strategies. These issues are particularly pertinent for many settings in Africa where there is high HBV prevalence and co-endemic HIV infection, but lack of robust epidemiological data and limited education, diagnostics and clinical care. The prevalence, distribution and impact of RAMs and VEMs in these populations are neglected in the current literature. We therefore set out to assimilate data for sub-Saharan Africa through a systematic literature review and analysis of published sequence data, and present these in an on-line database (https://livedataoxford.shinyapps.io/1510659619-3Xkoe2NKkKJ7Drg/). The majority of the data were from HIV/HBV coinfected cohorts. The commonest RAM was rtM204I/V, either alone or in combination with associated mutations, and identified in both reportedly treatment-naïve and treatment-experienced adults. We also identified the suite of mutations rtM204V/I + rtL180M + rtV173L, that has been associated with vaccine escape, in over 1/3 of cohorts. Although tenofovir has a high genetic barrier to resistance, it is of concern that emerging data suggest polymorphisms that may be associated with resistance, although the precise clinical impact of these is unknown. Overall, there is an urgent need for improved diagnostic screening, enhanced laboratory assessment of HBV before and during therapy, and sustained roll out of tenofovir in preference to lamivudine alone. Further data are needed in order to inform population and individual approaches to HBV diagnosis, monitoring and therapy in these highly vulnerable settings.

Applications of next-generation sequencing analysis for the detection of hepatocellular carcinoma-associated hepatitis B virus mutations

BACKGROUND

Next-generation sequencing (NGS) is a powerful and high-throughput method for the detection of viral mutations. This article provides a brief overview about optimization of NGS analysis for hepatocellular carcinoma (HCC)-associated hepatitis B virus (HBV) mutations, and hepatocarcinogenesis of relevant mutations.

MAIN BODY

For the application of NGS analysis in the genome of HBV, four noteworthy steps were discovered in testing. First, a sample-specific reference sequence was the most effective mapping reference for NGS. Second, elongating the end of reference sequence improved mapping performance at the end of the genome. Third, resetting the origin of mapping reference sequence could probed deletion mutations and variants at a certain location with common mutations. Fourth, using a platform-specific cut-off value to distinguish authentic minority variants from technical artifacts was found to be highly effective. One hundred and sixty-seven HBV single nucleotide variants (SNVs) were found to be studied previously through a systematic literature review, and 12 SNVs were determined to be associated with HCC by meta-analysis. From comprehensive research using a HBV genome-wide NGS analysis, 60 NGS-defined HCC-associated SNVs with their pathogenic frequencies were identified, with 19 reported previously. All the 12 HCC-associated SNVs proved by meta-analysis were confirmed by NGS analysis, except for C1766T and T1768A which were mainly expressed in genotypes A and D, but including the subgroup analysis of A1762T. In the 41 novel NGS-defined HCC-associated SNVs, 31.7% (13/41) had cut-off values of SNV frequency lower than 20%. This showed that NGS could be used to detect HCC-associated SNVs with low SNV frequency. Most SNV II (the minor strains in the majority of non-HCC patients) had either low (< 20%) or high (> 80%) SNV frequencies in HCC patients, a characteristic U-shaped distribution pattern. The cut-off values of SNV frequency for HCC-associated SNVs represent their pathogenic frequencies. The pathogenic frequencies of HCC-associated SNV II also showed a U-shaped distribution. Hepatocarcinogenesis induced by HBV mutated proteins through cellular pathways was reviewed.

CONCLUSION

NGS analysis is useful to discover novel HCC-associated HBV SNVs, especially those with low SNV frequency. The hepatocarcinogenetic mechanisms of novel HCC-associated HBV SNVs defined by NGS analysis deserve further investigation.

Ancient hepatitis B viruses from the Bronze Age to the Medieval period

Hepatitis B virus (HBV) is a major cause of human hepatitis. There is considerable uncertainty about the timescale of its evolution and its association with humans. Here we present 12 full or partial ancient HBV genomes that are between approximately 0.8 and 4.5 thousand years old. The ancient sequences group either within or in a sister relationship with extant human or other ape HBV clades. Generally, the genome properties follow those of modern HBV. The root of the HBV tree is projected to between 8.6 and 20.9 thousand years ago, and we estimate a substitution rate of 8.04 × 10^-6-1.51 × 10^-5 nucleotide substitutions per site per year. In several cases, the geographical locations of the ancient genotypes do not match present-day distributions. Genotypes that today are typical of Africa and Asia, and a subgenotype from India, are shown to have an early Eurasian presence. The geographical and temporal patterns that we observe in ancient and modern HBV genotypes are compatible with well-documented human migrations during the Bronze and Iron Ages^1,2. We provide evidence for the creation of HBV genotype A via recombination, and for a long-term association of modern HBV genotypes with humans, including the discovery of a human genotype that is now extinct. These data expose a complexity of HBV evolution that is not evident when considering modern sequences alone.

Correlation of oxidative stress in patients with HBV-induced liver disease with HBV genotypes and drug resistance mutations

OBJECTIVE

This study aims to explore the correlation of oxidative stress (OxS) in patients with chronic hepatitis B (CHB) and the disease severity with HBV genotypes and drug resistance mutations.

METHODS

A total of 296 patients with CHB were enrolled into the study. PCR-reverse dot-blot hybridization was used to detect the HBV genotypes (B, C, and D) and the drug resistance-causing HBV mutant genes. In addition, the total oxidative stress (TOS) and total antioxidant status (TAS) were determined, and oxidative stress index (OSI) was calculated and compared.

RESULTS

Serum levels of TOS and OSI, the B/C ratio, and drug resistance mutation rate were increased along with the elevated disease severity degree (CHB

CONCLUSION

There is oxidative damage in patients with HBV-induced liver disease, and the damage degree is correlated with the HBV genotype and drug resistance mutation. Oxidative stress might be a useful indicator of the progression of HBV-induced liver disease in patients.

Collapse

Key Words

• Drug resistance-causing mutation
• Genotype
• HBV infection
• Liver disease severity
• Oxidative stress

Collapse

MESH Headings

• Adult
• Drug Resistance, Viral
• Female
• Genotype
• Hepatitis B/blood
• Hepatitis B/drug therapy
• Hepatitis B/genetics
• Hepatitis B virus/genetics
• Humans
• Male
• Middle Aged
• Mutation
• Oxidative Stress

Collapse

Grants Collapse

Affiliation(s)

Jianbo Xianyu Department of General Surgery, The Affiliated Mianyang Central Hospital of Southwest Medical University, 621000, Sichuan, China
Jiafu Feng Department of Laboratory Medicine, The Affiliated Mianyang Central Hospital of Southwest Medical University, 621000, Sichuan, China.
Yuwei Yang Department of Laboratory Medicine, The Affiliated Mianyang Central Hospital of Southwest Medical University, 621000, Sichuan, China
Jie Tang Department of Laboratory Medicine, The Affiliated Mianyang Central Hospital of Southwest Medical University, 621000, Sichuan, China
Gang Xie Department of Clinical Pathology, The Affiliated Mianyang Central Hospital of Southwest Medical University, 621000, Sichuan, China
Lingying Fan Department of Laboratory Medicine, The Affiliated Mianyang Central Hospital of Southwest Medical University, 621000, Sichuan, China; Laboratory Medicine, Southwest Medical University, 646000 Luzhou, China

Collapse

A new hepatitis B virus e antigen-negative strain gene used as a reference sequence in an animal model

Infection with hepatitis B virus (HBV) e-antigen (HBeAg)-negative strains is increasingly prevalent. Currently, detailed information of the obtained natural HBV strain is not available except for the B genotype and HBeAg-negative. The aim of the present study was to characterize the natural genetic variation of the HBeAg-negative strain and investigate its function. The genic sequence was determined using Sanger sequencing, and compared to related sequences using alignment and phylogenetic analysis. In vivo, virus-specific serum markers were investigated in CBA/CaJ mice. The sequence had a full genome length of 3215 nucleotides. Sites 122, 125, 127, and 160 in S regions were identified as lysine, threonine, proline, and lysine respectively. The main four point variants including A1762T, G1764A, G1896A, and G1899A were detected in the full-length genome. The genotype of the sequence was B, with sub-genotype B2 and serological subtype adw2. The characterize of the natural genetic variation strain showed no reported drug-resistant variant in P region and no reported immune escape site in S region. The strain will increase viral replication and infection for mutations A1762T and G1764A in the basal core promoter region, and mutations G1896A and G1899A in the pre-core region. The G1896A variant resulted in a premature stop codon and abolished HBeAg expression. HBsAg persisted for 26 weeks and HBeAg was still negative in CBA/CaJ mice. The present sequence is representative of the HBeAg-negative genome and may serve as a valuable reference for studying HBeAg-negative strains. The present findings were successfully verified in CBA/CaJ mice, demonstrating good applicability of the sequence.