Dating the origin of hepatitis B virus reveals higher substitution rate and adaptation on the branch leading to F/H genotypes

The distinct spatiotemporal evolutionary landscape of HBV and HDV largely determines the unique epidemic features of HDV globally

Chronic infection of hepatitis B virus (HBV) and hepatitis D virus (HDV) causes the most severe form of viral hepatitis. Due to the dependence on HBV, HDV was deemed to co-evolve and co-migrate with HBV. However, we previously found that the naturally occurred HDV/HBV combinations do not always reflect the most efficient virological adaptation (Wang et al., 2021). Moreover, regions with heavy HBV burden do not always correlate with high HDV prevalence (e.g., East Asia), and vice versa (e.g., Central Asia). Herein, we systematically elucidated the spatiotemporal evolutionary landscape of HDV to understand the unique epidemic features of HDV. We found that the MRCA of HDV was from South America around the late 13th century, was globally dispersed mainly via Central Asia, and evolved into eight genotypes from the 19th to 20th century. In contrast, the MRCA of HBV was from Europe ∼23.7 thousand years ago (Kya), globally dispersed mainly via Africa and East Asia, and evolved into eight genotypes ∼1100 years ago. When HDV stepped in, all present-day HBV genotypes had already formed and its global genotypic distribution had stayed stable geographically. Nevertheless, regionalized HDV adapted to local HBV genotypes and human lineages, contributing to the global geographical separation of HDV genotypes. Additionally, a sharp increase in HDV infections was observed after the 20th century. In conclusion, HDV exhibited a distinct spatiotemporal distribution path compared with HBV. This unique evolutionary relationship largely fostered the unique epidemic features we observe nowadays. Moreover, HDV infections may continue to ramp up globally, thus more efforts are urgently needed to combat this disease.

High Prevalence of Hepatitis B Virus Drug Resistance Mutations to Lamivudine among People with HIV/HBV Coinfection in Rural and Peri-Urban Communities in Botswana

(1) Background: We aimed to determine the prevalence of hepatitis B virus (HBV) resistance-associated mutations (RAMs) in people with HBV and human immunodeficiency virus (HBV/HIV) in Botswana. (2) Methods: We sequenced HBV deoxyribonucleic acid (DNA) from participants with HBV/HIV from the Botswana Combination Prevention Project study (2013-2018) using the Oxford Nanopore GridION platform. Consensus sequences were analyzed for genotypic and mutational profiles. (3) Results: Overall, 98 HBV sequences had evaluable reverse transcriptase region coverage. The median participant age was 43 years (IQR: 37, 49) and 66/98 (67.4%) were female. Most participants, i.e., 86/98 (87.8%) had suppressed HIV viral load (VL). HBV RAMs were identified in 61/98 (62.2%) participants. Most RAMs were in positions 204 (60.3%), 180 (50.5%), and 173 (33.3%), mostly associated with lamivudine resistance. The triple mutations rtM204V/L180M/V173L were the most predominant (17/61 [27.9%]). Most participants (96.7%) with RAMs were on antiretroviral therapy for a median duration of 7.5 years (IQR: 4.8, 10.5). Approximately 27.9% (17/61) of participants with RAMs had undetectable HBV VL, 50.8% (31/61) had VL < 2000 IU/mL, and 13/61 (21.3%) had VL ≥ 2000 IU/mL. (4) Conclusions: The high prevalence of lamivudine RAMs discourages the use of ART regimens with 3TC as the only HBV-active drug in people with HIV/HBV.

The Hepatitis B Virus Genotypes E to J: The Overlooked Genotypes

Hepatitis B virus (HBV) genotypes E to J are understudied genotypes. Genotype E is found almost exclusively in West Africa. Genotypes F and H are found in America and are rare in other parts of the world. The distribution of genotype G is not completely known. Genotypes I and J are found in Asia and probably result from recombination events with other genotypes. The number of reported sequences for HBV genotypes E to J is small compared to other genotypes, which could impact phylogenetic and pairwise distance analyses. Genotype F is the most divergent of the HBV genotypes and is subdivided into six subgenotypes F1 to F6. Genotype E may be a recent genotype circulating almost exclusively in sub-Saharan Africa. Genotype J is a putative genotype originating from a single Japanese patient. The paucity of data from sub-Saharan Africa and Latin America is due to the under-representation of these regions in clinical and research cohorts. The purpose of this review is to highlight the need for further research on HBV genotypes E to J, which appear to be overlooked genotypes.

Hepatitis B Virus Genotype D: An Overview of Molecular Epidemiology, Evolutionary History, and Clinical Characteristics

The hepatitis B virus (HBV) genotype D (HBV/D) is the most extensively distributed genotype worldwide with distinct molecular and epidemiological features. This report provides an up-to-date review on the history of HBV/D subgenotyping and misclassifications, along with large-scale analysis of over 1000 HBV/D complete genome sequences, with the aim of gaining a thorough understanding of the global prevalence and geographic distribution of HBV/D subgenotypes. We have additionally explored recent paleogenomic findings, which facilitated the detection of HBV/D genomes dating back to the late Iron Age and provided new perspectives on the origins of modern HBV/D strains. Finally, reports on distinct disease outcomes and responses to antiviral therapy among HBV/D subgenotypes are discussed, further highlighting the complexity of this genotype and the importance of HBV subgenotyping in the management and treatment of hepatitis B.

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.

Women in the European Virus Bioinformatics Center

Viruses are the cause of a considerable burden to human, animal and plant health, while on the other hand playing an important role in regulating entire ecosystems. The power of new sequencing technologies combined with new tools for processing "Big Data" offers unprecedented opportunities to answer fundamental questions in virology. Virologists have an urgent need for virus-specific bioinformatics tools. These developments have led to the formation of the European Virus Bioinformatics Center, a network of experts in virology and bioinformatics who are joining forces to enable extensive exchange and collaboration between these research areas. The EVBC strives to provide talented researchers with a supportive environment free of gender bias, but the gender gap in science, especially in math-intensive fields such as computer science, persists. To bring more talented women into research and keep them there, we need to highlight role models to spark their interest, and we need to ensure that female scientists are not kept at lower levels but are given the opportunity to lead the field. Here we showcase the work of the EVBC and highlight the achievements of some outstanding women experts in virology and viral bioinformatics.

Detection of Ancient Viruses and Long-Term Viral Evolution

The COVID-19 outbreak has reminded us of the importance of viral evolutionary studies as regards comprehending complex viral evolution and preventing future pandemics. A unique approach to understanding viral evolution is the use of ancient viral genomes. Ancient viruses are detectable in various archaeological remains, including ancient people's skeletons and mummified tissues. Those specimens have preserved ancient viral DNA and RNA, which have been vigorously analyzed in the last few decades thanks to the development of sequencing technologies. Reconstructed ancient pathogenic viral genomes have been utilized to estimate the past pandemics of pathogenic viruses within the ancient human population and long-term evolutionary events. Recent studies revealed the existence of non-pathogenic viral genomes in ancient people's bodies. These ancient non-pathogenic viruses might be informative for inferring their relationships with ancient people's diets and lifestyles. Here, we reviewed the past and ongoing studies on ancient pathogenic and non-pathogenic viruses and the usage of ancient viral genomes to understand their long-term viral evolution.

The unexpected high prevalence of HBV subgenotype D4 in patients with chronic hepatitis B in Galicia, a northwestern Spanish region, reflects strong links with Latin America

BACKGROUND

Hepatitis B virus (HBV) comprises 9 genotypes and multiple subgenotypes that depict differences in geographic distribution, clinical outcome and response to antiviral therapy. However, the molecular epidemiology of HBV geno/subgenotypes is globally scarce. In Spain, HBV genotype D seems to be more prevalent in the northwestern regions compared to the rest of the country for unclear reasons.

METHODS

HBV genotyping was performed using geno2pheno on a S gene fragment amplified from plasma collected from all chronic hepatitis B individuals attended at one reference hospital in Santiago de Compostela, the Galicia's capital town. Phylogenetic and phylogeographic analyses using a fragment of 345 bp were performed in all viremic specimens. To avoid misleading allocation as consequence of short fragment analysis, several bioinformatic controls were used.

RESULTS

A total of 320 individuals with persistent serum HBsAg+ and detectable HBV-DNA were seen between 2000 and 2016 (male 68.4%; median age, 52 years-old; native Spaniards 83.8%). HBV genotype distribution was as follows: A 15.3%; B 1.6%; C 2.5%; D 71.6%; E 3.1%; F 2.2%; G 3.1%; and H 0.6%. HBV genotype D was mostly represented by D4 and D2 subgenotypes (33.4% and 15% of total, respectively). Compared to chronic hepatitis B patients with genotypes B, C, E and G, HBV-D4 carriers tended to be older (54.2% had >50 years-old) and HBeAg-negative (85%). Moreover, 43% were female, 4.7% had cirrhosis, 10.2% hepatitis C and 6.4% HIV coinfection. Phylogenetic analyses could be performed on 82 HBV-D4 specimens; and 79 were confirmed as HBV-D4 using PhyML. Phylogeography using FasTree suggested at least two distinct introductions of HBV-D4 in Galicia, one from the Caribbean and South America, and another from India.

CONCLUSIONS

HBV subgenotype D4 is the most prevalent HBV variant in chronic hepatitis B patients living in the northwest of Spain, representing 33.4% (107/320) of all chronic hepatitis B infections. This rate of HBV-D4 is among the highest reported worldwide. Epidemiological and phylogenetic analyses suggest a strong association with historical migrant exchanges with Latin America, and especially with the Caribbean basin.

Dating the origin and dispersal of global hepatitis B virus genotype C in humans

Hepatitis B virus genotype C (HBV/C) is one of the most prevalent HBV strains worldwide, especially in the Western Pacific and the South-East Asia. However, the origin and evolutionary timescale of HBV/C remains largely unresolved. We analyzed the evolutionary rate and molecular clock phylogeny of 101 full-genome HBV/C sequences sampled globally using a Bayesian Markov Chain Monte Carlo (MCMC) approach. We inferred the spatiotemporal dynamics of the HBV/C worldwide by the Bayesian Stochastic Search Variable Selection (BSSVS). We found that the estimated mean evolution rate of the HBV/C genotype full-genome was 4.32 × 10^-5 subs/site/year (95% highest posterior density 3.02 × 10^-6 - 8.97 × 10^-5). Phylogeographic reconstruction was able to identify a single location for the origin of the global HBV/C in Australia around A.D. 715. The subgenotype C4 diverged earliest and mainly circulated in Australia, C1 mainly in Southeast Asia, C2 mainly in East Asia and C3 in Remote Oceania. The effective number of HBV infection presented a rapid exponential increase between the 1760s and 1860s followed by a maintained high level until now. Our study, for the first time, provides an estimated timescale for the HBV/C epidemic, and brings new insight to the dispersal of HBV/C in humans globally. Based on the continuous presence of a highly effective viral population, this study provides further evidence of the challenge from a population-based molecular level to eliminate HBV by 2030, and calls for a concerted effort from policy makers, health providers, and society in the globalized world.

Hepatitis B Virus Genotype G: The Odd Cousin of the Family

With a widespread distribution but low prevalence worldwide, the hepatitis B virus (HBV) genotype G (HBV/G) is a recently described genotype for which the origin and biology are poorly understood. Some unique features make HBV/G the most peculiar of all genotypes. In this review, we reflect on the major milestones in HBV/G research, highlighting the main aspects of its discovery, molecular epidemiology, and virological and clinical characteristics. We also illustrate common pitfalls in the routine detection, which may lead to underestimated rates of HBV/G infection. Large-scale analysis of data from dozens of articles was further performed, with the aim of gaining comprehensive insights into the epidemiological aspects of HBV/G. Finally, we point out recent findings on HBV/G origins and discuss new perspectives regarding the evolutionary history of HBV/G and the plausibility of an African geographic re-emergence of this genotype.

Hepatitis B Virus Integration into Transcriptionally Active Loci and HBV-Associated Hepatocellular Carcinoma

Hepatitis B Virus (HBV) DNA integrations into the human genome are considered major causative factors to HBV-associated hepatocellular carcinoma development. In the present study, we investigated whether HBV preferentially integrates parts of its genome in specific genes and evaluated the contribution of the integrations in HCC development per gene. We applied dedicated in-house developed pipelines on all of the available HBV DNA integration data and performed a statistical analysis to identify genes that could be characterized as hotspots of integrations, along with the evaluation of their association with HBV-HCC. Our results suggest that 15 genes are recurrently affected by HBV integrations and they are significantly associated with HBV-HCC. Further studies that focus on HBV integrations disrupting these genes are mandatory in order to understand the role of HBV integrations in clonal advantage gain and oncogenesis promotion, as well as to determine whether inhibition of the HBV-disrupted genes can provide a therapy strategy for HBV-HCC.

Molecular archeology of human viruses

The evolution of human-virus associations is usually reconstructed from contemporary patterns of genomic diversity. An intriguing, though still rarely implemented, alternative is to search for the genetic material of viruses in archeological and medical archive specimens to document evolution as it happened. In this chapter, we present lessons from ancient DNA research and incorporate insights from virology to explore the potential range of applications and likely limitations of archeovirological approaches. We also highlight the numerous questions archeovirology will hopefully allow us to tackle in the near future, and the main expected roadblocks to these avenues of research.

Ten millennia of hepatitis B virus evolution

[Figure: see text].

Hepatitis B Virus DNA Integration, Chronic Infections and Hepatocellular Carcinoma

Hepatitis B Virus (HBV) is an Old World virus with a high mutation rate, which puts its origins in Africa alongside the origins of Homo sapiens, and is a member of the Hepadnaviridae family that is characterized by a unique viral replication cycle. It targets human hepatocytes and can lead to chronic HBV infection either after acute infection via horizontal transmission usually during infancy or childhood or via maternal-fetal transmission. HBV has been found in ~85% of HBV-related Hepatocellular Carcinomas (HCC), and it can integrate the whole or part of its genome into the host genomic DNA. The molecular mechanisms involved in the HBV DNA integration is not yet clear; thus, multiple models have been described with respect to either the relaxed-circular DNA (rcDNA) or the double-stranded linear DNA (dslDNA) of HBV. Various genes have been found to be affected by HBV DNA integration, including cell-proliferation-related genes, oncogenes and long non-coding RNA genes (lincRNAs). The present review summarizes the advances in the research of HBV DNA integration, focusing on the evolutionary and molecular side of the integration events along with the arising clinical aspects in the light of WHO's commitment to eliminate HBV and viral hepatitis by 2030.

Ancient viral genomes reveal introduction of human pathogenic viruses into Mexico during the transatlantic slave trade

After the European colonization of the Americas, there was a dramatic population collapse of the Indigenous inhabitants caused in part by the introduction of new pathogens. Although there is much speculation on the etiology of the Colonial epidemics, direct evidence for the presence of specific viruses during the Colonial era is lacking. To uncover the diversity of viral pathogens during this period, we designed an enrichment assay targeting ancient DNA (aDNA) from viruses of clinical importance and applied it to DNA extracts from individuals found in a Colonial hospital and a Colonial chapel (16th-18th century) where records suggest that victims of epidemics were buried during important outbreaks in Mexico City. This allowed us to reconstruct three ancient human parvovirus B19 genomes and one ancient human hepatitis B virus genome from distinct individuals. The viral genomes are similar to African strains, consistent with the inferred morphological and genetic African ancestry of the hosts as well as with the isotopic analysis of the human remains, suggesting an origin on the African continent. This study provides direct molecular evidence of ancient viruses being transported to the Americas during the transatlantic slave trade and their subsequent introduction to New Spain. Altogether, our observations enrich the discussion about the etiology of infectious diseases during the Colonial period in Mexico.

Comprehensive Analysis of Clinically Significant Hepatitis B Virus Mutations in Relation to Genotype, Subgenotype and Geographic Region

Hepatitis B virus (HBV) is a highly variable DNA virus due to its unique life cycle, which involves an error-prone reverse transcriptase. The high substitution rate drives the evolution of HBV by generating genetic variants upon which selection operates. HBV mutants with clinical implications have been documented worldwide, indicating the potential for spreading and developing their own epidemiology. However, the prevalence of such mutants among the different HBV genotypes and subgenotypes has not been systematically analyzed. In the current study, we performed large-scale analysis of 6,479 full-length HBV genome sequences from genotypes A-H, with the aim of gaining comprehensive insights into the relationships of relevant mutations associated with immune escape, antiviral resistance and hepatocellular carcinoma (HCC) development with HBV (sub)genotypes and geographic regions. Immune escape mutations were detected in 10.7% of the sequences, the most common being I/T126S (1.8%), G145R (1.2%), M133T (1.2%), and Q129R (1.0%). HBV genotype B showed the highest rate of escape mutations (14.7%) while genotype H had no mutations (P < 0.001). HCC-associated mutations were detected in 33.7% of the sequences, with significantly higher frequency of C1653T, T1753V and A1762T/G1764A in genotype G than C (P < 0.001). The overall frequencies of lamivudine-, telbivudine-, adefovir-, and entecavir-resistant mutants were 7.3, 7.2, 0.5, and 0.2%, respectively, while only 0.05% showed reduced susceptibility to tenofovir. In particular, the highest frequency of lamivudine-resistant mutations was observed in genotype G and the lowest frequency in genotype E (32.5 and 0.3%; P < 0.001). The prevalence of HBV mutants was also biased by geographic location, with North America identified as one of the regions with the highest rates of immune escape, antiviral resistance, and HCC-associated mutants. The collective findings were discussed in light of natural selection and the known characteristics of HBV (sub)genotypes. Our data provide relevant information on the prevalence of clinically relevant HBV mutations, which may contribute to further improvement of diagnostic procedures, immunization programs, therapeutic protocols, and disease prognosis.

The evolution and clinical impact of hepatitis B virus genome diversity

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

Long-term evolution of hepatitis B virus genotype F: Strong association between viral diversification and the prehistoric settlement of Central and South America

The genotype F (HBV-F) is an autochthonous Native American strain of the hepatitis B virus. In this study, we reconstruct the HBV-F long-term evolution under a hypothesis of co-divergence with humans in Central and South America, since their entry into the region 14.5-16 thousand years ago. The Bayesian phylogeographic reconstruction supported a virus-host co-expansion; however, two evolutionary scenarios would have been present. Whereas subgenotype F1 spreads along a Pacific coastal route and would have evolved associated with Central American and Andean cultures from the west of the continent, subgenotypes F2-F6 spread along the Atlantic coastline and inner pathways associated with communities inhabiting the tropical forest lowlands. Then, we propose a model for HBV-F evolution in which the selection of differential biological characteristics in these two main groups would be related to their evolution in host populations with different genetic backgrounds and dissimilar demographic conditions.

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.

Hepatitis B Virus: Alternative phylogenetic hypotheses and its impact on molecular evolution inferences

Characterizing molecular evolution patterns of the Hepatitis B Virus (HBV) is important for a better understanding of the natural history of this infection. However, several molecular evolution estimates are conditioned on tree topology. There is no consensus about the phylogenetic relationships of HBV genotypes, and different studies often find alternative topologies. While most studies consider HBV genotypes F and H as sister to all other human genotypes, a recent study suggested an alternative HBV phylogeny that indicates an accelerated substitution rate for HBV-F/H partially driven by positive selection. In this study, we evaluate the impact of alternative HBV topologies on inferences of HBV phylogeny, rate acceleration, and positive selection on the HBV-F/H branch. Our results indicate that under certain methodological approaches alternative HBV topologies are equally likely. Considering phylogenetic uncertainty, there is no evidence that HBV-F/H had an accelerated substitution rate, even though inferences of positive selection are robust to alternative background topologies. Our results further suggest that, under reasonable assumptions, HBV-F/H most likely represents the sister lineage to all other human/ape HBV genotypes. Understanding the full range of likely topologies will be crucial for elaborating, testing, and refining hypothesis about the evolutionary HBV origins in our species.

Mode and tempo of human hepatitis virus evolution

Human viral hepatitis, a major cause of morbidity and mortality worldwide, is caused by highly diverse viruses with different genetic, ecological, and pathogenetic features. Technological advances that allow throughput sequencing of viral genomes, as well as the development of computational tools to analyze such genome data, have largely expanded our knowledge on the host range and evolutionary history of human hepatitis viruses. Thus, with the exclusion of hepatitis D virus, close or distant relatives of these human pathogens were identified in a number of domestic and wild mammals. Also, sequences of human viral strains isolated from different geographic locations and over different time-spans have allowed the application of phylogeographic and molecular dating approaches to large viral phylogenies. In this review, we summarize the most recent insights into our understanding of the evolutionary events and ecological contexts that determined the origin and spread of human hepatitis viruses.

Evaluation of HBV-Like Circulation in Wild and Farm Animals from Brazil and Uruguay

The origin of the hepatitis B virus is a subject of wide deliberation among researchers. As a result, increasing academic interest has focused on the spread of the virus in different animal species. However, the sources of viral infection for many of these animals are unknown since transmission may occur from animal to animal, human to human, animal to human, and human to animal. The aim of this study was to evaluate hepadnavirus circulation in wild and farm animals (including animals raised under wild or free conditions) from different sites in Brazil and Uruguay using serological and molecular tools. A total of 487 domestic wild and farm animals were screened for hepatitis B virus (HBV) serological markers and tested via quantitative and qualitative polymerase chain reaction (PCR) to detect viral DNA. We report evidence of HBsAg (surface antigen of HBV) and total anti-HBc (HBV core antigen) markers as well as low-copy hepadnavirus DNA among domestic and wild animals. According to our results, which were confirmed by partial genome sequencing, as the proximity between humans and animals increases, the potential for pathogen dispersal also increases. A wider knowledge and understanding of reverse zoonoses should be sought for an effective One Health response.

Reconstruction of the spatial and temporal dynamics of hepatitis B virus genotype D in the Americas

Hepatitis B virus (HBV) genotype D (HBV/D) is globally widespread, and ten subgenotypes (D1 to D10) showing distinct geographic distributions have been described to date. The evolutionary history of HBV/D and its subgenotypes, for which few complete genome sequences are available, in the Americas is not well understood. The main objective of the current study was to determine the full-length genomic sequences of HBV/D isolates from Brazil and frequency, origin and spread of HBV/D subgenotypes in the Americas. Complete HBV/D genomes isolated from 39 Brazilian patients infected with subgenotypes D1 (n = 1), D2 (n = 10), D3 (n = 27), and D4 (n = 1) were sequenced and analyzed together with reference sequences using the Bayesian coalescent and phylogeographic framework. A search for HBV/D sequences available in GenBank revealed 209 complete and 926 partial genomes from American countries (Argentina, Brazil, Canada, Chile, Colombia, Cuba, Haiti, Martinique, Mexico, USA and Venezuela), with the major circulating subgenotypes identified as D1 (26%), D2 (17%), D3 (36%), D4 (21%), and D7 (1%) within the continent. The detailed evolutionary history of HBV/D in the Americas was investigated by using different evolutionary time scales. Spatiotemporal reconstruction analyses using short-term substitution rates suggested times of the most recent common ancestor for the American HBV/D subgenotypes coincident with mass migratory movements to Americas during the 19^th and 20^th centuries. In particular, significant linkages between Argentina and Syria (D1), Brazil and Central/Eastern Europe (D2), USA and India (D2), and Brazil and Southern Europe (D3) were estimated, consistent with historical and epidemiological data.

Tracing Ancient Human Migrations into Sahul Using Hepatitis B Virus Genomes

The entry point and timing of ancient human migration into continental Sahul (the combined landmass of Australia, New Guinea, and Tasmania) are subject to debate. Unique strains of hepatitis B virus (HBV) are endemic among modern-day Australian Aboriginals (HBV/C4) and Indigenous Melanesians (HBV/C3). We postulated that HBV genomes could be used to infer human population movements because the main HBV transmission route in endemic populations is via mother-to-child for genotypes B and C infections. Phylogenetic and phylogeographic analyses of HBV genomes inferred the origin of HBV/C4 to be >59 thousand years ago (ka) (95% HPD: 34-85 ka), and most likely to have occurred on the Sunda Shelf (southeast extension of the continental shelf of Southeast Asia). Our analysis further suggested an age of >51 ka (95% Highest Posterior Density (HPD): 36-67 ka) for the most recent common ancestor of HBV/C4 in Australia, correlating with the arrival time of anatomically modern humans into Australia, with the entry point suggested along a southern route via Timor. While we also inferred the origin of HBC/C3 to be on the Sunda Shelf, our analyses suggested that it was carried into Melanesia by Indigenous Melanesians who migrated through New Guinea north of the highlands. These findings reveal that HBV genomes can be used to infer ancient human population movements.

Evolutionary biology of human hepatitis viruses

Hepatitis viruses are major threats to human health. During the last decade, highly diverse viruses related to human hepatitis viruses were found in animals other than primates. Herein, we describe both surprising conservation and striking differences of the unique biological properties and infection patterns of human hepatitis viruses and their animal homologues, including transmission routes, liver tropism, oncogenesis, chronicity, pathogenesis and envelopment. We discuss the potential for translation of newly discovered hepatitis viruses into preclinical animal models for drug testing, studies on pathogenesis and vaccine development. Finally, we re-evaluate the evolutionary origins of human hepatitis viruses and discuss the past and present zoonotic potential of their animal homologues.

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

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

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

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

Unravelling the history of hepatitis B virus genotypes A and D infection using a full-genome phylogenetic and phylogeographic approach

Hepatitis B virus (HBV) infection constitutes a global public health problem. In order to establish how HBV was disseminated across different geographic regions, we estimated the levels of regional clustering for genotypes D and A. We used 916 HBV-D and 493 HBV-A full-length sequences to reconstruct their global phylogeny. Phylogeographic analysis was conducted by the reconstruction of ancestral states using the criterion of parsimony. The putative origin of genotype D was in North Africa/Middle East. HBV-D sequences form low levels of regional clustering for the Middle East and Southern Europe. In contrast, HBV-A sequences form two major clusters, the first including sequences mostly from sub-Saharan Africa, and the second including sequences mostly from Western and Central Europe. Conclusion: We observed considerable differences in the global dissemination patterns of HBV-D and HBV-A and different levels of monophyletic clustering in relation to the regions of prevalence of each genotype.

A novel hepatitis B virus species discovered in capuchin monkeys sheds new light on the evolution of primate hepadnaviruses

BACKGROUND & AIMS

All known hepatitis B virus (HBV) genotypes occur in humans and hominoid Old World non-human primates (NHPs). The divergent woolly monkey HBV (WMHBV) forms another orthohepadnavirus species. The evolutionary origins of HBV are unclear.

METHODS

We analysed sera from 124 Brazilian monkeys collected during 2012-2016 for hepadnaviruses using molecular and serological tools, and conducted evolutionary analyses.

RESULTS

We identified a novel orthohepadnavirus species in capuchin monkeys (capuchin monkey hepatitis B virus [CMHBV]). We found CMHBV-specific antibodies in five animals and high CMHBV concentrations in one animal. Non-inflammatory, probably chronic infection was consistent with an intact preCore domain, low genetic variability, core deletions in deep sequencing, and no elevated liver enzymes. Cross-reactivity of antisera against surface antigens suggested antigenic relatedness of HBV, CMHBV, and WMHBV. Infection-determining CMHBV surface peptides bound to the human HBV receptor (human sodium taurocholate co-transporting polypeptide), but preferentially interacted with the capuchin monkey receptor homologue. CMHBV and WMHBV pseudotypes infected human hepatoma cells via the human sodium taurocholate co-transporting polypeptide, and were poorly neutralised by HBV vaccine-derived antibodies, suggesting that cross-species infections may be possible. Ancestral state reconstructions and sequence distance comparisons associated HBV with humans, whereas primate hepadnaviruses as a whole were projected to NHP ancestors. Co-phylogenetic analyses yielded evidence for co-speciation of hepadnaviruses and New World NHP. Bayesian hypothesis testing yielded strong support for an association of the HBV stem lineage with hominoid ancestors. Neither CMHBV nor WMHBV was likely the ancestor of the divergent human HBV genotypes F/H found in American natives.

CONCLUSIONS

Our data suggest ancestral co-speciation of hepadnaviruses and NHP, and an Old World origin of the divergent HBV genotypes F/H. The identification of a novel primate hepadnavirus offers new perspectives for urgently needed animal models of chronic hepatitis B.

LAY SUMMARY

The origins of HBV are unclear. The new orthohepadnavirus species from Brazilian capuchin monkeys resembled HBV in elicited infection patterns and could infect human liver cells using the same receptor as HBV. Evolutionary analyses suggested that primate HBV-related viruses might have emerged in African ancestors of New World monkeys millions of years ago. HBV was associated with hominoid primates, including humans and apes, suggesting evolutionary origins of HBV before the formation of modern humans. HBV genotypes found in American natives were divergent from those found in American monkeys, and likely introduced along prehistoric human migration. Our results elucidate the evolutionary origins and dispersal of primate HBV, identify a new orthohepadnavirus reservoir, and enable new perspectives for animal models of hepatitis B.

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.

The paradox of HBV evolution as revealed from a 16th century mummy

Hepatitis B virus (HBV) is a ubiquitous viral pathogen associated with large-scale morbidity and mortality in humans. However, there is considerable uncertainty over the time-scale of its origin and evolution. Initial shotgun data from a mid-16th century Italian child mummy, that was previously paleopathologically identified as having been infected with Variola virus (VARV, the agent of smallpox), showed no DNA reads for VARV yet did for hepatitis B virus (HBV). Previously, electron microscopy provided evidence for the presence of VARV in this sample, although similar analyses conducted here did not reveal any VARV particles. We attempted to enrich and sequence for both VARV and HBV DNA. Although we did not recover any reads identified as VARV, we were successful in reconstructing an HBV genome at 163.8X coverage. Strikingly, both the HBV sequence and that of the associated host mitochondrial DNA displayed a nearly identical cytosine deamination pattern near the termini of DNA fragments, characteristic of an ancient origin. In contrast, phylogenetic analyses revealed a close relationship between the putative ancient virus and contemporary HBV strains (of genotype D), at first suggesting contamination. In addressing this paradox we demonstrate that HBV evolution is characterized by a marked lack of temporal structure. This confounds attempts to use molecular clock-based methods to date the origin of this virus over the time-frame sampled so far, and means that phylogenetic measures alone cannot yet be used to determine HBV sequence authenticity. If genuine, this phylogenetic pattern indicates that the genotypes of HBV diversified long before the 16th century, and enables comparison of potential pathogenic similarities between modern and ancient HBV. These results have important implications for our understanding of the emergence and evolution of this common viral pathogen.

Mutation and recombination in pathogen evolution: Relevance, methods and controversies

Mutation and recombination drive the evolution of most pathogens by generating the genetic variants upon which selection operates. Those variants can, for example, confer resistance to host immune systems and drug therapies or lead to epidemic outbreaks. Given their importance, diverse evolutionary studies have investigated the abundance and consequences of mutation and recombination in pathogen populations. However, some controversies persist regarding the contribution of each evolutionary force to the development of particular phenotypic observations (e.g., drug resistance). In this study, we revise the importance of mutation and recombination in the evolution of pathogens at both intra-host and inter-host levels. We also describe state-of-the-art analytical methodologies to detect and quantify these two evolutionary forces, including biases that are often ignored in evolutionary studies. Finally, we present some of our former studies involving pathogenic taxa where mutation and recombination played crucial roles in the recovery of pathogenic fitness, the generation of interspecific genetic diversity, or the design of centralized vaccines. This review also illustrates several common controversies and pitfalls in the analysis and in the evaluation and interpretation of mutation and recombination outcomes.

Tracing hepatitis B virus (HBV) genotype B5 (formerly B6) evolutionary history in the circumpolar Arctic through phylogeographic modelling

Background

Indigenous populations of the circumpolar Arctic are considered to be endemically infected (>2% prevalence) with hepatitis B virus (HBV), with subgenotype B5 (formerly B6) unique to these populations. The distinctive properties of HBV/B5, including high nucleotide diversity yet no significant liver disease, suggest virus adaptation through long-term host-pathogen association.

Methods

To investigate the origin and evolutionary spread of HBV/B5 into the circumpolar Arctic, fifty-seven partial and full genome sequences from Alaska, Canada and Greenland, having known location and sampling dates spanning 40 years, were phylogeographically investigated by Bayesian analysis (BEAST 2) using a reversible-jump-based substitution model and a clock rate estimated at 4.1 × 10⁻⁵ substitutions/site/year.

Results

Following an initial divergence from an Asian viral ancestor approximately 1954 years before present (YBP; 95% highest probability density interval [1188, 2901]), HBV/B5 coalescence occurred almost 1000 years later. Surprisingly, the HBV/B5 ancestor appears to locate first to Greenland in a rapid coastal route progression based on the landscape aware geographic model, with subsequent B5 evolution and spread westward. Bayesian skyline plot analysis demonstrated an HBV/B5 population expansion occurring approximately 400 YBP, coinciding with the disruption of the Neo-Eskimo Thule culture into more heterogeneous and regionally distinct Inuit populations throughout the North American Arctic.

Discussion

HBV/B5 origin and spread appears to occur coincident with the movement of Neo-Eskimo (Inuit) populations within the past 1000 years, further supporting the hypothesis of HBV/host co-expansion, and illustrating the concept of host-pathogen adaptation and balance.

Heterogeneous recombination among Hepatitis B virus genotypes

The rapid evolution of Hepatitis B virus (HBV) through both evolutionary forces, mutation and recombination, allows this virus to generate a large variety of adapted variants at both intra and inter-host levels. It can, for instance, generate drug resistance or the diverse viral genotypes that currently exist in the HBV epidemics. Concerning the latter, it is known that recombination played a major role in the emergence and genetic diversification of novel genotypes. In this regard, the quantification of viral recombination in each genotype can provide relevant information to devise expectations about the evolutionary trends of the epidemic. Here we measured the amount of this evolutionary force by estimating global and local recombination rates in >4700 HBV complete genome sequences corresponding to nine (A to I) HBV genotypes. Counterintuitively, we found that genotype E presents extremely high levels of recombination, followed by genotypes B and C. On the other hand, genotype G presents the lowest level, where recombination is almost negligible. We discuss these findings in the light of known characteristics of these genotypes. Additionally, we present a phylogenetic network to depict the evolutionary history of the studied HBV genotypes. This network clearly classified all genotypes into specific groups and indicated that diverse pairs of genotypes are derived from a common ancestor (i.e., C-I, D-E and, F-H) although still the origin of this virus presented large uncertainty. Altogether we conclude that the amount of observed recombination is heterogeneous among HBV genotypes and that this heterogeneity can influence on the future expansion of the epidemic.

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.