Nonproductive Hepatitis B Virus Covalently Closed Circular DNA Generates HBx-Related Transcripts from the HBx/Enhancer I Region and Acquires Reactivation by Superinfection in Single Cells

Hepatitis B virus (HBV) infection remains a public health problem worldwide. Persistent HBV infection relies on active transcription of the covalently closed circular DNA (cccDNA) in hepatocytes, which is less understood at the single-cell level. In this study, we isolated primary human hepatocytes from liver-humanized FRG mice infected with HBV and examined cccDNA transcripts in single cells based on 5' end sequencing. Our 5' transcriptome sequencing (RNA-seq) analysis unambiguously assigns different viral transcripts with overlapping 3' sequences and quantitatively measures viral transcripts for structural genes (3.5 kb, 2.4 kb, and 2.1 kb) and the nonstructural X gene (0.7 kb and related) in single cells. We found that an infected cell either can generate all viral transcripts, signifying active transcription, or presents only transcripts from the X gene and its associated enhancer I domain and no structural gene transcripts. Results from cell infection assays with recombinant HBV show that nonproductive transcription of cccDNA can be activated by incoming virus through superinfection. Moreover, upon HBV infection, cccDNA apparently can be transcribed in the absence of HBx and produces HBx, needed for productive transcription of other viral genes. These results shed new light on cccDNA transcription at the single-cell level and provide insights useful for improving the treatment strategy against chronic HBV infection. IMPORTANCE Hepatitis B virus (HBV) infection can be effectively suppressed but rarely cured by available drugs. Chronic HBV infection is based on persistence of covalently closed circular DNA (cccDNA) and continuous infection and reinfection with HBV in the liver. Understanding transcriptional regulation of cccDNA will help to achieve permanent transcriptional silencing, i.e., functional cure of HBV. In our study, we found that an infected cell either can generate all viral transcripts, signifying active transcription, or presents only transcripts from the X gene and its associated enhancer I domain and no structural gene transcripts. The nonproductive transcription of cccDNA can be activated by incoming virus through superinfection. Upon an infection, cccDNA apparently can be transcribed in the absence of HBx to produce HBx, necessary for subsequent transcription of other HBV genes. Our studies shed new light on the mechanism of HBV infection and may have implications for a functional cure regimen for HBV.

Polymorphisms predicting phylogeny in hepatitis B virus

Hepatitis B viruses (HBVs) are compact viruses with circular genomes of ∼3.2 kb in length. Four genes (HBx, Core, Surface, and Polymerase) generating seven products are encoded on overlapping reading frames. Ten HBV genotypes have been characterised (A-J), which may account for differences in transmission, outcomes of infection, and treatment response. However, HBV genotyping is rarely undertaken, and sequencing remains inaccessible in many settings. We set out to assess which amino acid (aa) sites in the HBV genome are most informative for determining genotype, using a machine learning approach based on random forest algorithms (RFA). We downloaded 5,496 genome-length HBV sequences from a public database, excluding recombinant sequences, regions with conserved indels, and genotypes I and J. Each gene was separately translated into aa, and the proteins concatenated into a single sequence (length 1,614 aa). Using RFA, we searched for aa sites predictive of genotype and assessed covariation among the sites with a mutual information-based method. We were able to discriminate confidently between genotypes A-H using ten aa sites. Half of these sites (5/10) sites were identified in Polymerase (Pol), of which 4/5 were in the spacer domain and one in reverse transcriptase. A further 4/10 sites were located in Surface protein and a single site in HBx. There were no informative sites in Core. Properties of the aa were generally not conserved between genotypes at informative sites. Among the highest co-varying pairs of sites, there were fifty-five pairs that included one of these 'top ten' sites. Overall, we have shown that RFA analysis is a powerful tool for identifying aa sites that predict the HBV lineage, with an unexpectedly high number of such sites in the spacer domain, which has conventionally been viewed as unimportant for structure or function. Our results improve ease of genotype prediction from limited regions of HBV sequences and may have future applications in understanding HBV evolution.

Spacer Domain in Hepatitis B Virus Polymerase: Plugging a Hole or Performing a Role?

Hepatitis B virus (HBV) polymerase is divided into terminal protein, spacer, reverse transcriptase, and RNase domains. Spacer has previously been considered dispensable, merely acting as a tether between other domains or providing plasticity to accommodate deletions and mutations. We explore evidence for the role of spacer sequence, structure, and function in HBV evolution and lineage, consider its associations with escape from drugs, vaccines, and immune responses, and review its potential impacts on disease outcomes.

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

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

Origin, Evolution and Stability of Overlapping Genes in Viruses: A Systematic Review

During their long evolutionary history viruses generated many proteins de novo by a mechanism called “overprinting”. Overprinting is a process in which critical nucleotide substitutions in a pre-existing gene can induce the expression of a novel protein by translation of an alternative open reading frame (ORF). Overlapping genes represent an intriguing example of adaptive conflict, because they simultaneously encode two proteins whose freedom to change is constrained by each other. However, overlapping genes are also a source of genetic novelties, as the constraints under which alternative ORFs evolve can give rise to proteins with unusual sequence properties, most importantly the potential for novel functions. Starting with the discovery of overlapping genes in phages infecting Escherichia coli, this review covers a range of studies dealing with detection of overlapping genes in small eukaryotic viruses (genomic length below 30 kb) and recognition of their critical role in the evolution of pathogenicity. Origin of overlapping genes, what factors favor their birth and retention, and how they manage their inherent adaptive conflict are extensively reviewed. Special attention is paid to the assembly of overlapping genes into ad hoc databases, suitable for future studies, and to the development of statistical methods for exploring viral genome sequences in search of undiscovered overlaps.

Recent Advances in Hepatitis B Treatment

Hepatitis B virus infection affects over 250 million chronic carriers, causing more than 800,000 deaths annually, although a safe and effective vaccine is available. Currently used antiviral agents, pegylated interferon and nucleos(t)ide analogues, have major drawbacks and fail to completely eradicate the virus from infected cells. Thus, achieving a "functional cure" of the infection remains a real challenge. Recent findings concerning the viral replication cycle have led to development of novel therapeutic approaches including viral entry inhibitors, epigenetic control of cccDNA, immune modulators, RNA interference techniques, ribonuclease H inhibitors, and capsid assembly modulators. Promising preclinical results have been obtained, and the leading molecules under development have entered clinical evaluation. This review summarizes the key steps of the HBV life cycle, examines the currently approved anti-HBV drugs, and analyzes novel HBV treatment regimens.

The characteristic of the synonymous codon usage and phylogenetic analysis of hepatitis B virus

Background

Hepatitis B virus (HBV) infection is a crucial medical issue worldwide. The dependence of HBV replication on host cell machineries and their co-evolutionary interactions prompt the codon usage pattern of viral genes to translation selection and mutation pressure.

Objective

The evolutionary characteristics of HBV and the natural selection effects of the human genome on the codon usage characteristics were analyzed to provide a basis for medication development for HBV infection.

Methods

The codon usage pattern of sequences from different HBV genotypes of our isolates and reference HBV genome sequences downloaded from the National Center for Biotechnology Information (NCBI) database were analyzed by computing the relative synonymous codon usage (RSCU), nucleotide content, codon adaptation index (CAI) and the effective number of codons (ENC).

Results

The highest ENC values were observed in the C genotypes, followed by the B genotypes. The ENC values indicated a weak codon usage bias (CUB) in HBV genome. The number of codons differentially used between the three genotypes was markedly higher than that of similarly used codons. High CAI values indicated a good adaptability of HBV to its host. The ENC plot indicated the occurrence of mutational pressure in the three genotypes. The mean Ka/Ks ratios in the three genotypes were lower than 1, which indicated a negative selection pressure. The CAI and GC3% plot indicated the existence of CUB in the HBV genome.

Conclusions

Nucleotide composition, mutation bias, negative selection and mutational pressure are key factors influencing the CUB and phylogenetic diversity in HBV genotypes. The data provided here could be useful for developing drugs for HBV infection.

Electronic supplementary material

The online version of this article (10.1007/s13258-020-00932-w) contains supplementary material, which is available to authorized users.

New insights into the evolutionary features of viral overlapping genes by discriminant analysis

Overlapping genes originate by a mechanism of overprinting, in which nucleotide substitutions in a pre-existing frame induce the expression of a de novo protein from an alternative frame. In this study, I assembled a dataset of 319 viral overlapping genes, which included 82 overlaps whose expression is experimentally known and the respective 237 homologs. Principal component analysis revealed that overlapping genes have a common pattern of nucleotide and amino acid composition. Discriminant analysis separated overlapping from non-overlapping genes with an accuracy of 97%. When applied to overlapping genes with known genealogy, it separated ancestral from de novo frames with an accuracy close to 100%. This high discriminant power was crucial to computationally design variants of de novo viral proteins known to possess selective anticancer toxicity (apoptin) or protection against neurodegeneration (X protein), as well as to detect two new potential overlapping genes in the genome of the new coronavirus SARS-CoV-2.

Asymmetric evolution in viral overlapping genes is a source of selective protein adaptation

Overlapping genes represent an intriguing puzzle, as they encode two proteins whose ability to evolve is constrained by each other. Overlapping genes can undergo “symmetric evolution” (similar selection pressures on the two proteins) or “asymmetric evolution” (significantly different selection pressures on the two proteins). By sequence analysis of 75 pairs of homologous viral overlapping genes, I evaluated their accordance with one or the other model. Analysis of nucleotide and amino acid sequences revealed that half of overlaps undergo asymmetric evolution, as the protein from one frame shows a number of substitutions significantly higher than that of the protein from the other frame. Interestingly, the most variable protein (often known to interact with the host proteins) appeared to be encoded by the de novo frame in all cases examined. These findings suggest that overlapping genes, besides to increase the coding ability of viruses, are also a source of selective protein adaptation.

•

A dataset of 80 pairs of homologous overlapping genes from viruses is examined.

•

Its analysis reveals that half of overlapping genes undergo asymmetric evolution.

•

The most variable gene product is that encoded by the de novo overlapping gene.

•

Overlapping genes evolving asymmetrically are a source of selective protein adaptation.

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.

How viral genetic variants and genotypes influence disease and treatment outcome of chronic hepatitis B. Time for an individualised approach?

Chronic hepatitis B virus (HBV) infection remains a global problem. Several HBV genotypes exist with different biology and geographical prevalence. Whilst the future aim of HBV treatment remains viral eradication, current treatment strategies aim to suppress the virus and prevent the progression of liver disease. Current strategies also involve identification of patients for treatment, namely those at risk of progressive liver disease. Identification of HBV genotype, HBV mutants and other predictive factors allow for tailoured treatments, and risk-surveillance pathways, such as hepatocellular cancer screening. In the future, these factors may enable stratification not only of treatment decisions, but also of patients at risk of higher relapse rates when current therapies are discontinued. Newer technologies, such as next-generation sequencing, to assess drug-resistant or immune escape variants and quasi-species heterogeneity in patients, may allow for more information-based treatment decisions between the clinician and the patient. This article serves to discuss how HBV genotypes and genetic variants impact not only upon the disease course and outcomes, but also current treatment strategies. Adopting a personalised genotypic approach may play a role in future strategies to combat the disease. Herein, we discuss new technologies that may allow more informed decision-making for response guided therapy in the battle against HBV.

Deciphering the Origin and Evolution of Hepatitis B Viruses by Means of a Family of Non-enveloped Fish Viruses

Hepatitis B viruses (HBVs), which are enveloped viruses with reverse-transcribed DNA genomes, constitute the family Hepadnaviridae. An outstanding feature of HBVs is their streamlined genome organization with extensive gene overlap. Remarkably, the ∼1,100 bp open reading frame (ORF) encoding the envelope proteins is fully nested within the ORF of the viral replicase P. Here, we report the discovery of a diversified family of fish viruses, designated nackednaviruses, which lack the envelope protein gene, but otherwise exhibit key characteristics of HBVs including genome replication via protein-primed reverse-transcription and utilization of structurally related capsids. Phylogenetic reconstruction indicates that these two virus families separated more than 400 million years ago before the rise of tetrapods. We show that HBVs are of ancient origin, descending from non-enveloped progenitors in fishes. Their envelope protein gene emerged de novo, leading to a major transition in viral lifestyle, followed by co-evolution with their hosts over geologic eras.

•

Nackednaviruses are non-enveloped fish viruses related to hepadnaviruses

•

Both virus families separated from a common ancestor >400 million years ago

•

The envelope protein gene of hepadnaviruses emerged through two distinct processes

•

Hepadnaviruses mainly co-evolve with hosts while nackednaviruses jump between hosts

Impact of the biased nucleotide composition of viral RNA genomes on RNA structure and codon usage

We are interested in the influence of nucleotide composition on the fundamental characteristics of the virus RNA genome. Most RNA viruses have genomes with a distinct nucleotide composition, e.g. ranging from minimally 12.9 % to maximally 40.3 % (C- and U-count, respectively, in coronavirus HKU). We present a global analysis of diverse virus types, including plus-strand, minus-strand and double-strand RNA viruses, for the impact of this nucleotide preference on the predicted structure of the RNA genome that is packaged in virion particles and on the codon usage in the viral open reading frames. Several virus-specific features will be described, but also some general conclusions were drawn. Without exception, the virus-specific nucleotide bias was enriched in the unpaired, single-stranded regions of the RNA genome, thus creating an even more striking virus-specific signature. We present a simple mechanism that is based on elementary aspects of RNA structure folding to explain this general trend. In general, the nucleotide bias was the major determinant of the virus-specific codon usages, thus limiting a role for codon selection and translational control. We will discuss molecular and evolutionary scenarios that may be responsible for the diverse nucleotide biases of RNA viruses.