Phylogenetic analysis of Puumala virus strains from Central Europe highlights the need for a full-genome perspective on hantavirus evolution

Genetic features of the Puumala virus (Hantaviridae: Orthohantavirus) identified in the Moscow region

INTRODUCTION

Puumala virus (family Hantaviridae, genus Orthohantavirus) is distributed in most regions of the European part of Russia. However, information about its genetic variants circulating on the territory of the Central Federal District is extremely scarce.

MATERIALS AND METHODS

Rodents' tissue samples were tested after reverse transcription by PCR for the presence of hantaviral RNA. The amplified fragments of the L segment were sequenced by the Sanger method. For two samples, sequences of all three segments were obtained using the NGS method. Phylogenetic trees were built in the MEGA-X software.

RESULTS

Puumala virus was found in six samples. Based on the phylogenetic analysis of sequences of three segments, the obtained genetic variants belong to the sublineage previously designated as W-RUS.

CONCLUSION

A genetic variant of the Puumala virus, belonging to the subline W-RUS, circulates on the territory of the Volokolamsk district of Moscow region.

Genetic Diversity of Puumala orthohantavirus in Rodents and Human Patients in Austria, 2012-2019

Puumala orthohantavirus (PUUV) has a wide distribution throughout Europe. Distinctive temporal patterns of spillover into the human population are related to population dynamics of the reservoir host, the bank vole (Clethrionomys glareolus). As the rodent host is tied to specific habitats with small individual ranges, PUUV genetic diversity is also highly correlated with geographic distance. Using sequenced portions of viral S and M segments, we determined whether geographic clusters were supported. Human cases of PUUV infections are concentrated in southeastern Austria. We detected four distinct genotypes: two genotypes of the Alpe-Adria (ALAD) lineage typically associated with southeast Europe, and two sublineages of the Central Europe (CE) lineage. One cluster of CE genotypes represents a phylogenetically distinct sublineage compared to previously reported CE clades, and extends the boundary of the CE lineage further south than previously reported.

Genetic diversity and evolution of Hantaan virus in China and its neighbors

BACKGROUND

Hantaan virus (HTNV; family Hantaviridae, order Bunyavirales) causes hemorrhagic fever with renal syndrome (HFRS), which has raised serious concerns in Eurasia, especially in China, Russia, and South Korea. Previous studies reported genetic diversity and phylogenetic features of HTNV in different parts of China, but the analyses from the holistic perspective are rare.

METHODOLOGY AND PRINCIPAL FINDINGS

To better understand HTNV genetic diversity and gene evolution, we analyzed all available complete sequences derived from the small (S) and medium (M) segments with bioinformatic tools. Eleven phylogenetic groups were defined and showed geographic clustering; 42 significant amino acid variant sites were found, and 19 of them were located in immune epitopes; nine recombinant events and eight reassortments with highly divergent sequences were found and analyzed. We found that sequences from Guizhou showed high genetic divergence, contributing to multiple lineages of the phylogenetic tree and also to the recombination and reassortment events. Bayesian stochastic search variable selection analysis revealed that Heilongjiang, Shaanxi, and Guizhou played important roles in HTNV evolution and migration; the virus may originate from Zhejiang Province in the eastern part of China; and the virus population size expanded from the 1980s to 1990s.

CONCLUSIONS/SIGNIFICANCE

These findings revealed the original and evolutionary features of HTNV, which will help to illustrate hantavirus epidemic trends, thus aiding in disease control and prevention.

Assessing Genome-Wide Diversity in European Hantaviruses through Sequence Capture from Natural Host Samples

Research on the ecology and evolution of viruses is often hampered by the limitation of sequence information to short parts of the genomes or single genomes derived from cultures. In this study, we use hybrid sequence capture enrichment in combination with high-throughput sequencing to provide efficient access to full genomes of European hantaviruses from rodent samples obtained in the field. We applied this methodology to Tula (TULV) and Puumala (PUUV) orthohantaviruses for which analyses from natural host samples are typically restricted to partial sequences of their tri-segmented RNA genome. We assembled a total of ten novel hantavirus genomes de novo with very high coverage (on average >99%) and sequencing depth (average >247×). A comparison with partial Sanger sequences indicated an accuracy of >99.9% for the assemblies. An analysis of two common vole (Microtus arvalis) samples infected with two TULV strains each allowed for the de novo assembly of all four TULV genomes. Combining the novel sequences with all available TULV and PUUV genomes revealed very similar patterns of sequence diversity along the genomes, except for remarkably higher diversity in the non-coding region of the S-segment in PUUV. The genomic distribution of polymorphisms in the coding sequence was similar between the species, but differed between the segments with the highest sequence divergence of 0.274 for the M-segment, 0.265 for the S-segment, and 0.248 for the L-segment (overall 0.258). Phylogenetic analyses showed the clustering of genome sequences consistent with their geographic distribution within each species. Genome-wide data yielded extremely high node support values, despite the impact of strong mutational saturation that is expected for hantavirus sequences obtained over large spatial distances. We conclude that genome sequencing based on capture enrichment protocols provides an efficient means for ecological and evolutionary investigations of hantaviruses at an unprecedented completeness and depth.

Prevalence of the Puumala orthohantavirus Strains in the Pre-Kama Area of the Republic of Tatarstan, Russia

Puumala orthohantavirus (PUUV) causes nephropathia epidemica (NE), a mild form of hemorrhagic fever with renal syndrome (HFRS) commonly diagnosed in Europe. The majority of HFRS cases in the European part of Russia are diagnosed in the Volga Federal District, which includes the Republic of Tatarstan (RT). The current study aims to analyze the genetic variability of PUUV in Pre-Kama region of the RT bounded by the Volga, Kama, and Vyatka rivers. In 2017, bank voles were caught in seven isolated forest traps in the Pre-Kama region and for the 26 PUUV-positive samples, the partial small (S), medium (M), and large (L) genome segment sequences were obtained and analyzed. It was determined that all identified PUUV strains belong to the Russian (RUS) genetic lineage; however, the genetic distance between strains is not directly correlated with the geographical distance between bank vole populations. One of the identified strains has S and L segments produced from one parental strain, while the M segment was supplied by another, suggesting that this strain could be the reassortant. We suggest that the revealed pattern of the PUUV strains distribution could be the result of a series of successive multidirectional migratory flows of the bank voles to the Pre-Kama region in the postglacial period.

Prediction of the Spatial Origin of Puumala Virus Infections Using L Segment Sequences Derived from a Generic Screening PCR

To screen diagnostic specimens for the presence of hantavirus genomes or to identify new hantaviruses in nature, the pan-hanta L-PCR assay, a broadly reactive nested reverse transcription polymerase chain reaction (RT-PCR) assay targeting the L segment, is highly preferred over other assays because of its universality and high sensitivity. In contrast, the geographic allocation of Puumala virus strains to defined outbreak regions in Germany was previously done based on S segment sequences. We show that the routinely generated partial L segment sequences resulting from the pan-hanta L-PCR assay provide sufficient phylogenetic signal to inform the molecular epidemiology of the Puumala virus. Consequently, an additional S segment analysis seems no longer necessary for the identification of the spatial origin of a virus strain.

Identifying the patterns and drivers of Puumala hantavirus enzootic dynamics using reservoir sampling

Hantaviruses are zoonotic hemorrhagic fever viruses for which prevention of human spillover remains the first priority in disease management. Tailored intervention measures require an understanding of the drivers of enzootic dynamics, commonly inferred from distorted human incidence data. Here, we use longitudinal sampling of approximately three decades of Puumala orthohantavirus (PUUV) evolution in isolated reservoir populations to estimate PUUV evolutionary rates, and apply these to study the impact of environmental factors on viral spread. We find that PUUV accumulates genetic changes at a rate of ∼10⁻⁴ substitutions per site per year and that land cover type defines the dispersal dynamics of PUUV, with forests facilitating and croplands impeding virus spread. By providing reliable short-term PUUV evolutionary rate estimates, this work facilitates the evaluation of spatial risk heterogeneity starting from timed phylogeographic reconstructions based on virus sampling in its animal reservoir, thereby side-stepping the need for difficult-to-collect human disease incidence data.

Reassortment events in the evolution of hantaviruses

Hantaviruses (order Bunyavirales, family Hantaviridae), known as important zoonotic human pathogens, possess the capacity to exchange genome segments via genetic reassortment due to their tri-segmented genome. Although not as frequent as in the arthropod-borne bunyaviruses, reports indicating reassortment events in the evolution of hantaviruses have been recently accumulating. The intra- and inter-lineage reassortment between closely related variants has been repeatedly reported for several hantaviruses including the rodent-borne human pathogens such as Sin Nombre virus, Puumala virus, Dobrava-Belgrade virus, or Hantaan virus as well as for the more recently recognized shrew-borne hantaviruses, Imjin and Seewis. Reassortment between more distantly related viruses was rarely found but seems to play a beneficial role in the process of crossing the host species barriers. Besides the findings based on phylogenetic studies of naturally occurring strains, hantavirus reassortants were generated also in in vitro studies. Interestingly, only reassortants with exchanged M segments could be generated suggesting that a high degree of genetic compatibility is required for the S and L segments while the exchange of M segment is better tolerated or is particularly beneficial. Altogether, the numerous reports on hantavirus reassortment, summarized in this review, clearly demonstrate that reassortment events play a significant role in hantavirus evolution and contributed to the currently recognized hantavirus diversity.