Heterogeneous recombination among Hepatitis B virus genotypes

Molecular Evolution of SARS-CoV-2 during the COVID-19 Pandemic

The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) produced diverse molecular variants during its recent expansion in humans that caused different transmissibility and severity of the associated disease as well as resistance to monoclonal antibodies and polyclonal sera, among other treatments. In order to understand the causes and consequences of the observed SARS-CoV-2 molecular diversity, a variety of recent studies investigated the molecular evolution of this virus during its expansion in humans. In general, this virus evolves with a moderate rate of evolution, in the order of 10-3-10-4 substitutions per site and per year, which presents continuous fluctuations over time. Despite its origin being frequently associated with recombination events between related coronaviruses, little evidence of recombination was detected, and it was mostly located in the spike coding region. Molecular adaptation is heterogeneous among SARS-CoV-2 genes. Although most of the genes evolved under purifying selection, several genes showed genetic signatures of diversifying selection, including a number of positively selected sites that affect proteins relevant for the virus replication. Here, we review current knowledge about the molecular evolution of SARS-CoV-2 in humans, including the emergence and establishment of variants of concern. We also clarify relationships between the nomenclatures of SARS-CoV-2 lineages. We conclude that the molecular evolution of this virus should be monitored over time for predicting relevant phenotypic consequences and designing future efficient treatments.

Consequences of Genetic Recombination on Protein Folding Stability

Genetic recombination is a common evolutionary mechanism that produces molecular diversity. However, its consequences on protein folding stability have not attracted the same attention as in the case of point mutations. Here, we studied the effects of homologous recombination on the computationally predicted protein folding stability for several protein families, finding less detrimental effects than we previously expected. Although recombination can affect multiple protein sites, we found that the fraction of recombined proteins that are eliminated by negative selection because of insufficient stability is not significantly larger than the corresponding fraction of proteins produced by mutation events. Indeed, although recombination disrupts epistatic interactions, the mean stability of recombinant proteins is not lower than that of their parents. On the other hand, the difference of stability between recombined proteins is amplified with respect to the parents, promoting phenotypic diversity. As a result, at least one third of recombined proteins present stability between those of their parents, and a substantial fraction have higher or lower stability than those of both parents. As expected, we found that parents with similar sequences tend to produce recombined proteins with stability close to that of the parents. Finally, the simulation of protein evolution along the ancestral recombination graph with empirical substitution models commonly used in phylogenetics, which ignore constraints on protein folding stability, showed that recombination favors the decrease of folding stability, supporting the convenience of adopting structurally constrained models when possible for inferences of protein evolutionary histories with recombination.

Methodologies for Microbial Ancestral Sequence Reconstruction

The reconstruction of genetic material of ancestral organisms constitutes a powerful application of evolutionary biology. A fundamental step in this inference is the ancestral sequence reconstruction (ASR), which can be performed with diverse methodologies implemented in computer frameworks. However, most of these methodologies ignore evolutionary properties frequently observed in microbes, such as genetic recombination and complex selection processes, that can bias the traditional ASR. From a practical perspective, here I review methodologies for the reconstruction of ancestral DNA and protein sequences, with particular focus on microbes, and including biases, recommendations, and software implementations. I conclude that microbial ASR is a complex analysis that should be carefully performed and that there is a need for methods to infer more realistic ancestral microbial sequences.

Molecular characterization and phylogenetic analyses of full-length viral genomes from Iranian patients with chronic hepatitis B virus

Aim: Chronic hepatitis B infection is the main cause of liver complications such as hepatic failure, liver cirrhosis and hepatocellular carcinoma (HCC). In this study, we attempted to evaluate molecular characterization and phylogenetic analyses of full-length viral genomes from chronic hepatitis B virus (HBV)-infected patients. Methods: The full-length genomic sequence of the five HBV isolates from Ahvaz (city of Iran) patients was amplified, cloned in pTZ57R/T vector, sequenced and examined. Results: Phylogenetic analyses showed that all isolates belonged to genotype D (D1/D3). Serotyper tool identified ayw2 serotype in all HBV isolates. YMDE mutation was detected in an HBV isolate in the reverse transcriptase domain. Conclusion: In the present study, the analyses of full-length sequence of genome revealed that the HBV genotype D, sub-genotype D1/D3, and subtype ayw2 were predominant among Ahvaz HBV strains. As HBV genome replicates and is mediated via reverse transcription process, periodic investigations of full HBV genome are needed.

Global Patterns of Recombination across Human Viruses

Viral recombination is a major evolutionary mechanism driving adaptation processes, such as the ability of host-switching. Understanding global patterns of recombination could help to identify underlying mechanisms and to evaluate the potential risks of rapid adaptation. Conventional approaches (e.g., those based on linkage disequilibrium) are computationally demanding or even intractable when sequence alignments include hundreds of sequences, common in viral data sets. We present a comprehensive analysis of recombination across 30 genomic alignments from viruses infecting humans. In order to scale the analysis and avoid the computational limitations of conventional approaches, we apply newly developed topological data analysis methods able to infer recombination rates for large data sets. We show that viruses, such as ZEBOV and MARV, consistently displayed low levels of recombination, whereas high levels of recombination were observed in Sarbecoviruses, HBV, HEV, Rhinovirus A, and HIV. We observe that recombination is more common in positive single-stranded RNA viruses than in negatively single-stranded RNA ones. Interestingly, the comparison across multiple viruses suggests an inverse correlation between genome length and recombination rate. Positional analyses of recombination breakpoints along viral genomes, combined with our approach, detected at least 39 nonuniform patterns of recombination (i.e., cold or hotspots) in 18 viral groups. Among these, noteworthy hotspots are found in MERS-CoV and Sarbecoviruses (at spike, Nucleocapsid and ORF8). In summary, we have developed a fast pipeline to measure recombination that, combined with other approaches, has allowed us to find both common and lineage-specific patterns of recombination among viruses with potential relevance in viral adaptation.

Genetic diversity of hepatitis B virus in Yunnan, China: identification of novel subgenotype C17, an intergenotypic B/I recombinant, and B/C recombinants

Yunnan is considered to be a geographical hotspot for the introduction, mutation and recombination of several viruses in China. However, there are limited data regarding the genotypic profiles of hepatitis B virus (HBV) in this region. In this study, we characterized 206 HBV strains isolated from chronic hepatitis B patients in Yunnan, China. Initial genotyping based on 1.5 kb sequences revealed that genotype C was the most prevalent at 52.4 % (108/206), followed by genotype B at 30.6 % (63/206) and unclassified genotypes at 17.0 % (35/206). To characterize the 35 unclassified strains, 32 complete HBV genomes were amplified and analysed; 17 isolates were classified within a known subgenotype, 8 were classified as B/C recombinants, 1 was classified as a B/I recombinant and 6 constituted a potentially novel C subgenotype that we designated as C17, based on the characteristics of a monophyletic cluster, >4 % genetic distances, no significant evidence of recombination and no epidemiological link among individuals. Thus, multiple subgenotypes - namely B1, B2, B4, C1, C2, C3, C4, C8 and C17 - and two distinct intergenotypic recombinants exist in Yunnan, China, highlighting the complex and diverse distribution pattern of HBV genotypic profiles.

HBV evolution and genetic variability: Impact on prevention, treatment and development of antivirals

Hepatitis B virus (HBV) poses a major global health burden with 260 million people being chronically infected and 890,000 dying annually from complications in the course of the infection. HBV is a small enveloped virus with a reverse-transcribed DNA genome that infects hepatocytes and can cause acute and chronic infections of the liver. HBV is endemic in humans and apes representing the prototype member of the viral family Hepadnaviridae and can be divided into 10 genotypes. Hepadnaviruses have been found in all vertebrate classes and constitute an ancient viral family that descended from non-enveloped progenitors more than 360 million years ago. The de novo emergence of the envelope protein gene was accompanied with the liver-tropism and resulted in a tight virus-host association. The oldest HBV genomes so far have been isolated from human remains of the Bronze Age and the Neolithic (~7000 years before present). Despite the remarkable stability of the hepadnaviral genome over geological eras, HBV is able to rapidly evolve within an infected individual under pressure of the immune response or during antiviral treatment. Treatment with currently available antivirals blocking intracellular replication of HBV allows controlling of high viremia and improving liver health during long-term therapy of patients with chronic hepatitis B (CHB), but they are not sufficient to cure the disease. New therapy options that cover all HBV genotypes and emerging viral variants will have to be developed soon. In addition to the antiviral treatment of chronically infected patients, continued efforts to expand the global coverage of the currently available HBV vaccine will be one of the key factors for controlling the rising global spread of HBV. Certain improvements of the vaccine (e.g. inclusion of PreS domains) could counteract known problems such as low or no responsiveness of certain risk groups and waning anti-HBs titers leading to occult infections, especially with HBV genotypes E or F. But even with an optimal vaccine and a cure for hepatitis B, global eradication of HBV would be difficult to achieve because of an existing viral reservoir in primates and bats carrying closely related hepadnaviruses with zoonotic potential.

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.

High-throughput sequencing (HTS) for the analysis of viral populations

The development of High-Throughput Sequencing (HTS) technologies is having a major impact on the genomic analysis of viral populations. Current HTS platforms can capture nucleic acid variation across millions of genes for both selected amplicons and full viral genomes. HTS has already facilitated the discovery of new viruses, hinted new taxonomic classifications and provided a deeper and broader understanding of their diversity, population and genetic structure. Hence, HTS has already replaced standard Sanger sequencing in basic and applied research fields, but the next step is its implementation as a routine technology for the analysis of viruses in clinical settings. The most likely application of this implementation will be the analysis of viral genomics, because the huge population sizes, high mutation rates and very fast replacement of viral populations have demonstrated the limited information obtained with Sanger technology. In this review, we describe new technologies and provide guidelines for the high-throughput sequencing and genetic and evolutionary analyses of viral populations and metaviromes, including software applications. With the development of new HTS technologies, new and refurbished molecular and bioinformatic tools are also constantly being developed to process and integrate HTS data. These allow assembling viral genomes and inferring viral population diversity and dynamics. Finally, we also present several applications of these approaches to the analysis of viral clinical samples including transmission clusters and outbreak characterization.

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.

Molecular epidemiology, phylogenetic analysis and genotype distribution of hepatitis B virus in Saudi Arabia: Predominance of genotype D1

Despite the implementation of various vaccination programs, hepatitis B virus (HBV) poses a considerable health problem in Saudi Arabia. Insight on HBV evolutionary history in the region is limited. We performed a comprehensive epidemiological and phylogenetic reconstruction based on a large cohort of HBV infected patients. Three hundred and nineteen HBV-infected patients with different clinical manifestations, including inactive and active chronic carriers and patients with cirrhosis and hepatocellular carcinoma (HCC), were enrolled in this study. The full-length large S gene was amplified and sequenced. Phylogenetic analysis was performed to determine the genotype and subgenotypes of the isolates. Phylogenetic tree analysis revealed that genotype D is the most dominant genotype among patients. Moreover, this analysis identified two strains with genotype E isolated from active carriers. Detailed phylogenetic analyses confirmed the presence of four HBV D subgenotypes, D1 (93%, n = 296), D2 (0.02%, n = 5), D3 (0.003%, n = 1), and D4 (0.003%, n = 1). In addition, six genotype D strains were not assigned to any existing HBV D subgenotype. The large S gene of eight strains showed signatures of genotype recombination between the genotypes D and A and between D and E. Several strains harbored medically important point mutations at the protein level. Along with the dominance of the HBV genotype D, isolation of the E genotype and several recombinant strains from patients with Saudi Arabian origin is an essential result for decisions involving therapeutic measures for patients. Development of vaccines and detection of diagnostic escape mutations at antigenic epitopes on the HBsAg will be valuable to public health authorities. Furthermore, the diversity at the nucleotide and amino acid levels and different proportions of dN/dS at the PreS1, PreS2, and HBsAg reveal the selective pressure trend from inactive status towards advanced liver diseases.

Novel subgenotype D11 of hepatitis B virus in NaPo County, Guangxi, bordering Vietnam

Hepatitis B virus has been classified into 10 genotypes and 48 subgenotypes worldwide. We found previously, through polymerase chain reaction (PCR) amplification of a sample collected in 2011, that an HBsAg carrier was infected with two genotypes (B and D) of HBV. We carried out cloning, sequencing and phylogenetic analysis of the complete genomes and, for confirmation, analysed a sample collected from the same individual in 2018. Fifteen complete sequences were obtained from each sample. The carrier was infected in 2011 by genotypes B and D and by various recombinants, but only genotype D was present in 2018. The major and minor parents of the recombinants are genotypes B and D, respectively, although the recombination breakpoints vary among them. All 23 genotype D isolates form a cluster, branching out from other subgenotype D sequences and supported by a 100 % bootstrap value. Based on complete genome sequences, almost all of the estimated intragroup nucleotide divergence values between our isolates and HBV subgenotypes D1-D10 exceed 4 %. Compared to the other subgenotypes (D1-D10), 35 unique amino acids were present in our isolates. Our data provide evidence for a novel subgenotype, provisionally designated HBV subgenotype D11.

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.

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.

Recombination rates along the entire Epstein Barr virus genome display a highly heterogeneous landscape

Epstein Barr virus (EBV) has a large DNA genome assumed to be stable, but also subject to mutational processes such as nucleotide substitution and recombination, the latter explored to a lesser extent. Moreover, differences in the extent of recombination events across herpes sub-families were recently reported. Given the relevance of recombination in viral evolution and its possible impact in pathogenesis, we aimed to fully characterize and quantify its extension in all available EBV complete genome by assessing global and local recombination rate values (⍴/bp). Our results provide the first EBV recombination map based on recombination rates assessment, both at a global and gene by gene level, where the mean value for the entire genome was 0.035 (HPDI 0.020-0.062) ⍴/bp. We quantified how this evolutionary process changes along the EBV genome, and proved it to be non-homogeneous, since regulatory regions depicted the lowest recombination rate values while repetitive regions the highest signal. Moreover, GC content rich regions seem not to be linked to high recombination rates as previously reported. At an intragenic level, four genes (EBNA3C, EBNA3B, BRRF2 and BBLF2-BBLF3) presented a recombination rate above genome average. We specifically quantified the signal strength among different recombination-initiators previously described features and concluded that those which elicited the greatest amount of changes in ⍴/bp, TGGAG and CCCAG, were two well characterized recombination inducing motifs in eukaryotic cells. Strikingly, although TGGAG was not the most frequently detected DNA motif across the EBV genome (697 hits), it still induced a significantly greater proportion of initiation events (0.025 events/hits) than other more represented motifs, p-value = 0.04; one tailed proportion test. Present results support the idea that diversity and evolution of herpesviruses are impacted by mechanisms, such as recombination, which extends beyond the usual consideration of point mutations.

Neolithic and medieval virus genomes reveal complex evolution of hepatitis B

The hepatitis B virus (HBV) is one of the most widespread human pathogens known today, yet its origin and evolutionary history are still unclear and controversial. Here, we report the analysis of three ancient HBV genomes recovered from human skeletons found at three different archaeological sites in Germany. We reconstructed two Neolithic and one medieval HBV genome by de novo assembly from shotgun DNA sequencing data. Additionally, we observed HBV-specific peptides using paleo-proteomics. Our results demonstrated that HBV has circulated in the European population for at least 7000 years. The Neolithic HBV genomes show a high genomic similarity to each other. In a phylogenetic network, they do not group with any human-associated HBV genome and are most closely related to those infecting African non-human primates. The ancient viruses appear to represent distinct lineages that have no close relatives today and possibly went extinct. Our results reveal the great potential of ancient DNA from human skeletons in order to study the long-time evolution of blood borne viruses.

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.