Phylogeography of hepatitis B: The role of Portugal in the early dissemination of HBV worldwide

In Portugal, the genetic diversity, origin of HBV, and the Portuguese role in the dissemination of HBV worldwide were never investigated. In this work, we studied the epidemic history and transmission dynamics of HBV genotypes that are endemic in Portugal. HBV pol gene was sequenced from 130 patients followed in Lisbon. HBV genotype A (HBV/A) was the most prevalent (n=54, 41.5%), followed by D [HBV/D; (n=44, 33.8%)], and E [HBV/E; (n=32, 24.6%)]. Spatio-temporal evolutionary dynamics was reconstructed in BEAST using a Bayesian Markov Chain Monte Carlo method, with a GTR nucleotide substitution model, an uncorrelated lognormal relaxed molecular clock model, a Bayesian skyline plot, and a continuous diffusion model. HBV/D4 was the first subgenotype to be introduced in Portugal around 1857 (HPD 95% 1699-1931) followed by HBV/D3 and A2 a few decades later. HBV/E and HBV/A1 were introduced in Portugal later, almost simultaneously. Our results also indicate a very important role of Portugal in the exportation of HBV/D4 and A2 to Brazil and Cape Verde, respectively, at the beginning of the XX century. This work clarifies the epidemiological history of HBV in Portugal and shows that Portugal had an important role in the global spread of this virus.

Genetic diversity in enhancer II region of HBV genotype D and its association with advanced liver diseases

Background Hepatitis B Virus (HBV) is one of the most common human infectious agents, and the mutations in its genome may cause chronic hepatitis (CH), liver cirrhosis (LC), and hepatocellular carcinoma (HCC). This study was designed to characterize the enhancer-II (Enh-II) region of X gene in HBV positive patients to assess the association of such mutations with CH, LC, and HCC. Methods HBV positive samples (N = 200) with patients’ demographic and clinical data were collected from different regions of Khyber Pakhtunkhwa (KP), Pakistan. The Enh-II region of the HBx gene was sequenced and zanalyzed for polymorphism associated with advanced liver disease. Univariate and logistic regression analyses were performed to evaluate potent mutations associated with a risk for LC and HCC. Results HBV Enh-II region sequences analysis revealed 25 different mutations. The highest frequency of mutations S101F (62.2%), A102V/R/G/I (56.25%), M103L/A (68.75%)were found in HCC, followed in LC and CH patients as 57.1%, 42.8%, 28.52% 16%, 15.2% and 18.4% respectively. H94 deletion in the α-box of the Enh-II region, associated with a high risk of HCC was found in half of the HCC patients. This deletion was present in 28.5% of LC and 6.5% of CH patients. Importantly, the high frequency of some notable mutations such as E109A/Y, A110S/K, Y111D/E, and F112L was first time reported in the entire study population. The frequencies of these mutations were high in HCC (43.75%, 37.5%, 50% and 43.75% respectively) as compared to LC (14.28%, 14.28%, 28.2% and 42.8%) and CH patients (12.8%, 15.2%, 16.8% and 16% respectively). Conclusion Mutations associated with LC and HCC are prevalent in the Enh-II region in Pakistani HBV isolates. The mutations found are alarming in CH patients as these may progress to LC and HCC in a large number of patients.

Discrepant Results of Hepatitis B Virus Genotype Determination by PCR and DNA Sequencing

Currently, multiplex-PCR with genotype-specific primers is widely used for preliminary screening of hepatitis B virus (HBV) genotyping, despite its relatively lower accuracy compared with whole genome sequencing. Here, we present the discrepant results of HBV genotyping by PCR and full-length sequencing. HBV DNA was isolated from chronic hepatitis B serum and the HBV genotype was detected by PCR using genotype-specific primers and full-length genome sequencing. As a result, the determination of genotype B by the PCR method was consistent with the DNA sequencing results; however, PCR revealed that genotype C exhibited a mixed genotype of B and C in the current study. In conclusion, the PCR-based genotyping method may not provide accurate information of the HBV genotype and whole genome sequencing remains the “gold standard” method for HBV genotyping.

Canonical and Divergent N-Terminal HBx Isoform Proteins Unveiled: Characteristics and Roles during HBV Replication

Hepatitis B virus (HBV) X protein (HBx) is a viral regulatory and multifunctional protein. It is well-known that the canonical HBx reading frame bears two phylogenetically conserved internal in-frame translational initiation codons at Met2 and Met3, thus possibly generating divergent N-terminal smaller isoforms during translation. Here, we demonstrate that the three distinct HBx isoforms are generated from the ectopically expressed HBV HBx gene, named XF (full-length), XM (medium-length), and XS (short-length); they display different subcellular localizations when expressed individually in cultured hepatoma cells. Particularly, the smallest HBx isoform, XS, displayed a predominantly cytoplasmic localization. To study HBx proteins during viral replication, we performed site-directed mutagenesis to target the individual or combinatorial expression of the HBx isoforms within the HBV viral backbone (full viral genome). Our results indicate that of all HBx isoforms, only the smallest HBx isoform, XS, can restore WT levels of HBV replication, and bind to the viral mini chromosome, thereby establishing an active chromatin state, highlighting its crucial activities during HBV replication. Intriguingly, we found that sequences of HBV HBx genotype H are devoid of the conserved Met3 position, and therefore HBV genotype H infection is naturally silent for the expression of the HBx XS isoform. Finally, we found that the HBx XM (medium-length) isoform shares significant sequence similarity with the N-terminus domain of the COMMD8 protein, a member of the copper metabolism MURR1 domain-containing (COMMD) protein family. This novel finding might facilitate studies on the phylogenetic origin of the HBV X protein. The identification and functional characterization of its isoforms will shift the paradigm by changing the concept of HBx from being a unique, canonical, and multifunctional protein toward the occurrence of different HBx isoforms, carrying out different overlapping functions at different subcellular localizations during HBV genome replication. Significantly, our current work unveils new crucial HBV targets to study for potential antiviral research, and human virus pathogenesis.

Variability in the response of HBV D-subgenotypes to antiviral therapy: designing pan D-subgenotypic reverse transcriptase inhibitors

Nucleos(t)ide analogues entecavir (ETV ) and tenofovir disoproxil fumarate (TDF ) are recommended as first - line monotherapies for chronic hepatitis B (CHB). Multiple HBV genotypes/subgenotypes have been described, but their impact on treatment response remains largely elusive. We investigated the effectiveness of ETV/TDF on HBV/D-subgenotypes, D1/D2/D3/D5, studied the structural/functional differences in subgenotype-specific reverse transcriptase (RT) domains of viral polymerase and identified novel molecules with robust inhibitory activity on various D-subgenotypes. Transfection of Huh7 cells with full-length D1/D2/D3/D5 and in vitro TDF/ETV susceptibility assays demonstrated that D1/D2 had greater susceptibility to TDF/ETV while D3/D5 exhibited poorer response. Additionally, HBV load was substantially reduced in TDF-treated CHB patients carrying D1/D2 but minimally reduced in D3/D5-infected patients. Comparison of RT sequences of D-subgenotypes led to identification of unique subgenotype-specific residues and molecular modeling/docking/simulation studies depicted differential bindings of TDF/ETV to the active site of their respective RTs. Replacement of signature residues in D3/D5 HBV clones with corresponding amino acids seen in D1/D2 improved their susceptibility to TDF/ETV. Using high throughput virtual screening, we identified N(9)-[3 - fluoro - 2-(phosphonomethoxy)propyl] (FPMP) derivatives of purine bases, including N 6 -substituted ( S )-FPMP derivative of 2,6-diaminopurine (DAP) (OB-123-VK), as potential binders of RT of different D-subgenotypes. We synthesized ( S )-FPMPG prodrugs (FK-381-FEE/FK-381-SEE/FK-382) and tested their effectiveness along with OB-123-VK. Both OB-123-VK and FK-381-FEE exerted similar antiviral activities against all D-subgenotypes, although FK-381-FEE was more potent. Our study highlighted the natural variation in therapeutic response of D1/D2/D3/D5 and emphasized the need for HBV subgenotype determination before treatment. Novel molecules described here could benefit future design/discovery of pan-D-subgenotypic inhibitors. Importance: Current treatment of chronic hepatitis B relies heavily on nucleotide/nucleoside analogs in particular, tenofovir disoproxil fumarate (TDF) and entecavir (ETV) to keep HBV replication under control and prevent end-stage liver diseases. However, it was unclear whether the therapeutic effects of TDF/ETV differ among patients infected with different HBV genotypes and subgenotypes. HBV genotype D is the most widespread of all HBV genotypes and multiple D-subgenotypes have been described. We here report that different subgenotypes of HBV genotype-D exhibit variable response towards TDF and ETV and this could be attributed to naturally occurring amino acid changes in the reverse transcriptase domain of the subgenotype-specific polymerase. Further, we identified novel molecules and also synthesized prodrugs that are equally effective on different D-subgenotypes and could facilitate management of HBV/D-infected patients irrespective of D-subgenotype.

Understanding the Genetics of Viral Drug Resistance by Integrating Clinical Data and Mining of the Scientific Literature

Drug resistance caused by mutations is a public health threat for existing and emerging viral diseases. A wealth of evidence about these mutations and their clinically-associated phenotypes is scattered across the literature, but a comprehensive perspective is usually lacking. This work aimed to produce a clinically-relevant view for the case of Hepatitis B virus (HBV) mutations by combining a chronic HBV clinical study with a compendium of genetic mutations systematically gathered from the scientific literature. We enriched clinical mutation data by systematically mining 2,472,725 scientific articles from PubMed Central in order to gather information about the HBV mutational landscape. By performing this analysis, we were able to identify mutational hotspots for each HBV genotype (A-E) and gene (C, X, P, S), as well as the location of disulfide bonds associated with these mutations. Through a modelling study, we also identified a mutational position common in both the clinical data and the literature that is located at the binding pocket for a known anti-HBV drug, namely entecavir. The results of this novel approach shows the potential of integrated analyses to assist in the development of new drugs for viral diseases that are more robust to resistance. Such analyses should be of particular interest due to the increasing importance of viral resistance in established and emerging viruses, such as for newly-developed drugs against SARS-CoV-2.

Complex structural variations in non-human primate hepatitis B virus

Background Recent genome sequence technology has revealed a novel type of genetic rearrangement referred to as complex structural variations (SVs). Previous studies have elucidated the complex SVs in human hepatitis B viruses (HBVs). In this study, we investigated the existence of complex SVs in HBVs from non-human primates (NHPs). Methods Searches for nucleotide sequences of NHP HBV were conducted using the PubMed, and genetic sequences were retrieved from databases. The candidate genetic sequences harboring complex SVs were analyzed using the CLUSTALW program and MAFFT. Additional bioinformatical analyses were performed to determine strains with complex SVs and to elucidate characteristics of NHP HBV strains. Results One hundred and fifty-four HBV strains from NHPs were identified from databases. SVs and complex SVs were observed in 11 (7.1%) strains. Three gibbon HBV (GiHBV) strains showed complex SVs consisting of an insertion and a deletion in the pre-S1 region. One GiHBV strain possessed a 6-nt insertion, which are normally specific to human HBV genotype A (HBV/A) in the Core region, and further analyses clarified that the 6-nt insertion was not caused by recombination, but rather by simple insertion. Another chimpanzee HBV strain showed complex SVs in the pre-S1 region, which were composed of human HBV/E, G-specific polymorphic SV, and an additional 6-nt insertion. Conclusions In this study, complex SVs were observed in HBV strains from NHPs, in addition to human HBV strains, as shown in previous studies. These data suggest that complex SVs could also be found in other members of hepadnaviruses, and may play a role in their genetic diversity.

Molecular Evolution of Hepatitis B Virus Genotype F in Latin America

Hepatitis B virus (HBV) genotype F evolution is not completely understood in Latin America. This study aims to evaluate the molecular evolution of HBV-F in Latin America by comparing 224 whole-genome sequences. Bayesian coalescent analysis was performed to estimate the time to the most recent common ancestor. Four main clades were formed dated back between 1245 and 1730. Also, four subclades were identified dated back between 1705 and 1801. HBV-F overall effective population size grew in the 18th century and showed an initial circulation of HBV-F from Venezuela to other countries from Latin America.

Virome in adult Aedes albopictus captured during different seasons in Guangzhou City, China

Background The mosquito Aedes albopictus is an important vector for many pathogens. Understanding the virome in Ae. albopictus is critical for assessing the risk of disease transmission, implementation of vector control measures, and health system strengthening. Methods In this study, viral metagenomic and PCR methods were used to reveal the virome in adult Ae. albopictus captured in different areas and during different seasons in Guangzhou, China. Results The viral composition of adult Ae. albopictus varied mainly between seasons. Over 50 viral families were found, which were specific to vertebrates, invertebrates, plants, fungi, bacteria, and protozoa. In rural areas, Siphoviridae (6.5%) was the most common viral family harbored by mosquitoes captured during winter and spring, while Luteoviridae (1.1%) was the most common viral family harbored by mosquitoes captured during summer and autumn. Myoviridae (7.0% and 1.3%) was the most common viral family in mosquitoes captured in urban areas during all seasons. Hepatitis B virus (HBV) was detected by PCR in a female mosquito pool. The first near full-length HBV genome from Ae. albopictus was amplified, which showed a high level of similarity with human HBV genotype B sequences. Human parechovirus (HPeV) was detected in male and female mosquito pools, and the sequences were clustered with HPeV 1 and 3 sequences. Conclusions Large numbers of viral species were found in adult Ae. albopictus, including viruses from vertebrates, insects, and plants. The viral composition in Ae. albopictus mainly varied between seasons. Herein, we are the first to report the detection of HPeV and HBV in mosquitoes. This study not only provides valuable information for the control and prevention of mosquito-borne diseases, but it also demonstrates the feasibility of xenosurveillance. Graphical Abstract