Highly diversified shrew hepatitis B viruses corroborate ancient origins and divergent infection patterns of mammalian hepadnaviruses.

Shrews, insectivorous small mammals, pertain to an ancient mammalian order. We screened 693 European and African shrews for hepatitis B virus (HBV) homologs to elucidate the enigmatic genealogy of HBV. Shrews host HBVs at low prevalence (2.5%) across a broad geographic and host range. The phylogenetically divergent shrew HBVs comprise separate species termed crowned shrew HBV (CSHBV) and musk shrew HBV (MSHBV), each containing distinct genotypes. Recombination events across host orders, evolutionary reconstructions, and antigenic divergence of shrew HBVs corroborated ancient origins of mammalian HBVs dating back about 80 million years. Resurrected CSHBV replicated in human hepatoma cells, but human- and tupaia-derived primary hepatocytes were resistant to hepatitis D viruses pseudotyped with CSHBV surface proteins. Functional characterization of the shrew sodium taurocholate cotransporting polypeptide (Ntcp), CSHBV/MSHBV surface peptide binding patterns, and infection experiments revealed lack of Ntcp-mediated entry of shrew HBV. Contrastingly, HBV entry was enabled by the shrew Ntcp. Shrew HBVs universally showed mutations in their genomic preCore domains impeding hepatitis B e antigen (HBeAg) production and resembling those observed in HBeAg-negative human HBV. Deep sequencing and in situ hybridization suggest that HBeAg-negative shrew HBVs cause intense hepatotropic monoinfections and low within-host genomic heterogeneity. Geographical clustering and low MSHBV/CSHBV-specific seroprevalence suggest focal transmission and high virulence of shrew HBVs. HBeAg negativity is thus an ancient HBV infection pattern, whereas Ntcp usage for entry is not evolutionarily conserved. Shrew infection models relying on CSHBV/MSHBV revertants and human HBV will allow comparative assessments of HBeAg-mediated HBV pathogenesis, entry, and species barriers.

Molecular epidemiology and phylogenetics of human enteroviruses: Is there a forest behind the trees?

Enteroviruses are among the best studied small non-enveloped enteric RNA viruses. Most enteroviruses are easy to isolate in cell culture, and many non-polio enterovirus strains were archived worldwide as a byproduct of the WHO poliovirus surveillance system. Common outbreaks and epidemics, most prominently the epidemic of hand-foot-and-mouth disease with severe neurological complications in East and South-East Asia, justify practical interest of non-polio enteroviruses. As a result, there are over 50 000 enterovirus nucleotide sequences available in GenBank. Technical possibilities have been also improving, as Bayesian phylogenetic methods with an integrated molecular clock were introduced a decade ago and provided unprecedented opportunities for phylogenetic analysis. As a result, hundreds of papers were published on the molecular epidemiology of enteroviruses. This review covers the modern methodology, structure, and biases of the sequence dataset available in GenBank. The relevance of the subtype classification, findings of co-circulation of multiple genetic variants, previously unappreciated complexity of viral populations, and global evolutionary patterns are addressed. The most relevant conclusions and prospects for further studies on outbreak emergence mechanisms are discussed.

Molecular evolution of types in non-polio enteroviruses.

Non-polio enteroviruses are a ubiquitous and divergent group of non-enveloped RNA viruses. Novel types are reported regularly in addition to over 100 known types; however, mechanisms of emergence of novel types remain obscure. Here, the 33 most common types represented by 35-629 non-redundant partial VP1 sequences in GenBank were studied in parallel using Bayesian coalescent molecular clock analysis to investigate common evolutionary trends among enterovirus types. Inferred substitution rates were in the range of 0.41×10-2 to 3.07×10-2 substitutions per site per year. The most recent common ancestors of known isolates of each type presumably existed between 55 and 200 years ago. Phylogenetic analysis results suggested that global type populations underwent bottlenecks that could repeatedly reset the common ancestor dates. Nevertheless, species-level analysis suggested that the contemporary enterovirus types emerged within the last millennium. Analysis of 2657 complete VP1 sequences of the 24 most common types indicated that the type criterion based upon 75 % nucleotide sequence identity remains generally valid, despite exponential growth of the number of known sequences and a high rate of mutation fixation. However, in few types there was evidence that enteroviruses can drift slightly beyond the type threshold, up to 73 % identity, and both amino acid and nucleotide sequences should be considered for type identification. Analysis of sequence distances within types implied that sequence-identity-based identification of genotypes is rational within some, but not all, types and distinct genotype cut-offs (9-20 %) may be useful for different types.

Patterns and Temporal Dynamics of Natural Recombination in Noroviruses.

Noroviruses infect a wide range of mammals and are the major cause of gastroenteritis in humans. Recombination at the junction of ORF1 encoding nonstructural proteins and ORF2 encoding major capsid protein VP1 is a well-known feature of noroviruses. Using all available complete norovirus sequences, we systematically analyzed patterns of natural recombination in the genus Norovirus both throughout the genome and across the genogroups. Recombination events between nonstructural (ORF1) and structural genomic regions (ORF2 and ORF3) were found in all analyzed genogroups of noroviruses, although recombination was most prominent between members of GII, the most common genogroup that infects humans. The half-life times of recombinant forms (clades without evidence of recombination) of human GI and GII noroviruses were 10.4 and 8.4-11.3 years, respectively. There was evidence of many recent recombination events, and most noroviruses that differed by more than 18% of nucleotide sequence were recombinant relative to each other. However, there were no distinct recombination events between viruses that differed by over 42% in ORF2/3, consistent with the absence of systematic recombination between different genogroups. The few inter-genogroup recombination events most likely occurred between ancient viruses before they diverged into contemporary genogroups. The recombination events within ORF1 or between ORF2/3 were generally rare. Thus, noroviruses routinely exchange full structural and nonstructural blocks of the genome, providing a modular evolution.

The International Virus Bioinformatics Meeting 2023.

The 2023 International Virus Bioinformatics Meeting was held in Valencia, Spain, from 24-26 May 2023, attracting approximately 180 participants worldwide. The primary objective of the conference was to establish a dynamic scientific environment conducive to discussion, collaboration, and the generation of novel research ideas. As the first in-person event following the SARS-CoV-2 pandemic, the meeting facilitated highly interactive exchanges among attendees. It served as a pivotal gathering for gaining insights into the current status of virus bioinformatics research and engaging with leading researchers and emerging scientists. The event comprised eight invited talks, 19 contributed talks, and 74 poster presentations across eleven sessions spanning three days. Topics covered included machine learning, bacteriophages, virus discovery, virus classification, virus visualization, viral infection, viromics, molecular epidemiology, phylodynamic analysis, RNA viruses, viral sequence analysis, viral surveillance, and metagenomics. This report provides rewritten abstracts of the presentations, a summary of the key research findings, and highlights shared during the meeting.

Cases of Acute Flaccid Paralysis Associated with Coxsackievirus A2: Findings of a 20-Year Surveillance in the Russian Federation.

Surveillance for acute flaccid paralysis syndrome (AFP) in children under 15 is the backbone of the Global Polio Eradication Initiative. Laboratory examination of stool samples from AFP cases allows the detection of, along with polioviruses, a variety of non-polio enteroviruses (NPEV). The etiological significance of these viruses in the occurrence of AFP cases has been definitively established only for enteroviruses A71 and D68. Enterovirus Coxsackie A2 (CVA2) is most often associated with vesicular pharyngitis and hand, foot and mouth disease. Among 7280 AFP cases registered in Russia over 20 years (2001-2020), CVA2 was isolated only from five cases. However, these included three children aged 3 to 4 years, without overt immune deficiency, immunized with 4-5 doses of poliovirus vaccine in accordance with the National Vaccination Schedule. The disease resulted in persistent residual paralysis. Clinical and laboratory data corresponded to poliomyelitis developing during poliovirus infection. These findings are compatible with CVA2 being the cause of AFP. Molecular analysis of CVA2 from these patients and a number of AFP cases in other countries did not reveal association with a specific phylogenetic group, suggesting that virus genetics is unlikely to explain the pathogenic profile. The overall results highlight the value of AFP surveillance not just for polio control but for studies of uncommon AFP agents.

Modular Evolution of Coronavirus Genomes.

The viral family Coronaviridae comprises four genera, termed Alpha-, Beta-, Gamma-, and Deltacoronavirus. Recombination events have been described in many coronaviruses infecting humans and other animals. However, formal analysis of the recombination patterns, both in terms of the involved genome regions and the extent of genetic divergence between partners, are scarce. Common methods of recombination detection based on phylogenetic incongruences (e.g., a phylogenetic compatibility matrix) may fail in cases where too many events diminish the phylogenetic signal. Thus, an approach comparing genetic distances in distinct genome regions (pairwise distance deviation matrix) was set up. In alpha, beta, and delta-coronaviruses, a low incidence of recombination between closely related viruses was evident in all genome regions, but it was more extensive between the spike gene and other genome regions. In contrast, avian gammacoronaviruses recombined extensively and exist as a global cloud of genes with poorly corresponding genetic distances in different parts of the genome. Spike, but not other structural proteins, was most commonly exchanged between coronaviruses. Recombination patterns differed between coronavirus genera and corresponded to the modular structure of the spike: recombination traces were more pronounced between spike domains (N-terminal and C-terminal parts of S1 and S2) than within domains. The variability of possible recombination events and their uneven distribution over the genome suggest that compatibility of genes, rather than mechanistic or ecological limitations, shapes recombination patterns in coronaviruses.

Evaluating the Impact of Anthropogenic Factors on the Dissemination of Contemporary Cosmopolitan, Arctic, and Arctic-like Rabies Viruses.

Rabies is a globally prevalent viral zoonosis that causes 59,000 deaths per year and has important economic consequences. Most virus spread is associated with the migration of its primary hosts. Anthropogenic dissemination, mainly via the transportation of rabid dogs, shaped virus ecology a few hundred years ago and is responsible for several current outbreaks. A systematic analysis of aberrant long-distance events in the steppe and Arctic-like groups of rabies virus was performed using statistical (Bayesian) phylogeography and plots of genetic vs. geographic distances. The two approaches produced similar results but had some significant differences and complemented each other. No phylogeographic analysis could be performed for the Arctic group because polar foxes transfer the virus across the whole circumpolar region at high velocity, and there was no correlation between genetic and geographic distances in this virus group. In the Arctic-like group and the steppe subgroup of the cosmopolitan group, a significant number of known sequences (15-20%) was associated with rapid long-distance transfers, which mainly occurred within Eurasia. Some of these events have been described previously, while others have not been documented. Most of the recent long-distance transfers apparently did not result in establishing the introduced virus, but a few had important implications for the phylogeographic history of rabies. Thus, human-mediated long-distance transmission of the rabies virus remains a significant threat that needs to be addressed.

TBEV Subtyping in Terms of Genetic Distance.

Currently, the lowest formal taxon in virus classification is species; however, unofficial lower-level units are commonly used in everyday work. Tick-borne encephalitis virus (TBEV) is a species of mammalian tick-borne flaviviruses that may cause encephalitis. Many known representatives of TBEV are grouped into subtypes, mostly according to their phylogenetic relationship. However, the emergence of novel sequences could dissolve this phylogenetic grouping; in the absence of strict quantitative criterion, it may be hard to define the borders of the first TBEV taxonomic unit below the species level. In this study, the nucleotide/amino-acid space of all known TBEV sequences was analyzed. Amino-acid sequence p-distances could not reliably distinguish TBEV subtypes. Viruses that differed by less than 10% of nucleotides in the polyprotein-coding gene belonged to the same subtype. At the same time, more divergent viruses were representatives of different subtypes. According to this distance criterion, TBEV species may be divided into seven subtypes: TBEV-Eur, TBEV-Sib, TBEV-FE, TBEV-2871 (TBEV-Ob), TBEV-Him, TBEV-178-79 (TBEV-Bkl-1), and TBEV-886-84 (TBEV-Bkl-2).

Tick-Borne Encephalitis Virus: An Emerging Ancient Zoonosis?

Tick-borne encephalitis (TBE) is one of the most important viral zoonosis transmitted by the bite of infected ticks. In this study, all tick-borne encephalitis virus (TBEV) E gene sequences available in GenBank as of June 2019 with known date of isolation (n = 551) were analyzed. Simulation studies showed that a sample bias could significantly affect earlier studies, because small TBEV datasets (n = 50) produced non-overlapping intervals for evolutionary rate estimates. An apparent lack of a temporal signal in TBEV, in general, was found, precluding molecular clock analysis of all TBEV subtypes in one dataset. Within all subtypes and most of the smaller groups in these subtypes, there was evidence of many medium- and long-distance virus transfers. These multiple random events may play a key role in the virus spreading. For some groups, virus diversity within one territory was similar to diversity over the whole geographic range. This is best exemplified by the virus diversity observed in Switzerland or Czech Republic. These two countries yielded most of the known European subtype Eu3 subgroup sequences, and the diversity of viruses found within each of these small countries is comparable to that of the whole Eu3 subgroup, which is prevalent all over Central and Eastern Europe. Most of the deep tree nodes within all three established TBEV subtypes dated less than 300 years back. This could be explained by the recent emergence of most of the known TBEV diversity. Results of bioinformatics analysis presented here, together with multiple field findings, suggest that TBEV may be regarded as an emerging disease.

Emerging Concepts and Challenges in Rheumatoid Arthritis Gene Therapy.

Rheumatoid arthritis (RA) is a systemic inflammatory joint disease affecting about 1% of the population worldwide. Current treatment approaches do not ensure a cure for every patient. Moreover, classical regimens are based on nontargeted systemic immune suppression and have significant side effects. Biological treatment has advanced considerably but efficacy and specificity issues remain. Gene therapy is one of the potential future directions for RA therapy, which is rapidly developing. Several gene therapy trials done so far have been of moderate success, but experimental and genetics studies have yielded novel targets. As a result, the arsenal of gene therapy tools keeps growing. Currently, both viral and nonviral delivery systems are used for RA therapy. Herein, we review recent approaches for RA gene therapy.

Baltic Group Tick-Borne Encephalitis Virus Phylogeography: Systemic Inconsistency Pattern between Genetic and Geographic Distances.

Tick-Borne Encephalitis Virus (TBEV) is a dangerous arbovirus widely distributed in Northern Eurasia. The area of this pathogen changes over time. At the beginning of the 2000s, the Ixodes tick populations in Karelia increased. At the same time, the area of I. persulcatus, the main vector of the Siberian TBEV subtype, also expanded. Herein, we sequenced 10 viruses isolated from ticks collected in three locations from the Karelia region in 2008-2018. PCR positive samples were passaged in suckling mice or pig embryo kidney cells (PEK). After the second passage in suckling, mice viral RNA was isolated and E-gene fragment was sequenced. Viral sequences were expected to be similar or nearly identical. Instead, there was up to a 4.8% difference in nucleotide sequence, comparable with the most diverse viruses belonging to the Baltic subgroup in Siberian TBEV subtype (Baltic TBEV-Sib). To reveal whether this was systemic or incidental, a comprehensive phylogeographical analysis was conducted. Interestingly, viruses within each geographic region demonstrated comparable diversity to the whole Baltic TBEV-Sib. Moreover, Baltic TBEV-Sib has a distribution area limited by three ecological regions. This means that active virus mixing occurs in the vast geographic area forming one common virus pool. The most plausible explanation is the involvement of flying animals in the TBEV spread.

The Effect of Sample Bias and Experimental Artefacts on the Statistical Phylogenetic Analysis of Picornaviruses.

Statistical phylogenetic methods are a powerful tool for inferring the evolutionary history of viruses through time and space. The selection of mathematical models and analysis parameters has a major impact on the outcome, and has been relatively well-described in the literature. The preparation of a sequence dataset is less formalized, but its impact can be even more profound. This article used simulated datasets of enterovirus sequences to evaluate the effect of sample bias on picornavirus phylogenetic studies. Possible approaches to the reduction of large datasets and their potential for introducing additional artefacts were demonstrated. The most consistent results were obtained using "smart sampling", which reduced sequence subsets from large studies more than those from smaller ones in order to preserve the rare sequences in a dataset. The effect of sequences with technical or annotation errors in the Bayesian framework was also analyzed. Sequences with about 0.5% sequencing errors or incorrect isolation dates altered by just 5 years could be detected by various approaches, but the efficiency of identification depended upon sequence position in a phylogenetic tree. Even a single erroneous sequence could profoundly destabilize the whole analysis by increasing the variance of the inferred evolutionary parameters.

Using Statistical Phylogenetics for Investigation of Enterovirus 71 Genotype A Reintroduction into Circulation.

Neurovirulent enterovirus 71 (EV-A71) caused a massive epidemic in China in 2008-2011. While subgenotype C4 was the major causative agent, a few isolates were almost identical to the prototype EV-A71 strain and belonged to genotype A. This variant was allegedly extinct since 1970, and its identification in this epidemic suggests reintroduction of the archive virus. Regression analysis of genetic distances (TempEst software) was of moderate utility due to the low resolution of classical phylogenetic methods. Bayesian phylogenetic analysis (BEAST software) suggested artificial introduction event based on highly aberrant phylogenetic tree branch rates that differed by over three standard deviations from the mean substitution rate for EV71. Manual nucleotide-level analysis was used to further explore the virus spread pattern after introduction into circulation. Upon reintroduction, the virus accumulated up to seven substitutions in VP1, most of them non-synonymous and located within the capsid's canyon or at its rims, compatible with readaptation of a lab strain to natural circulation.