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

Genome-wide support for incipient Tula hantavirus species within a single rodent host lineage

Evolutionary divergence of viruses is most commonly driven by co-divergence with their hosts or through isolation of transmission after host shifts. It remains mostly unknown, however, whether divergent phylogenetic clades within named virus species represent functionally equivalent byproducts of high evolutionary rates or rather incipient virus species. Here, we test these alternatives with genomic data from two widespread phylogenetic clades in Tula orthohantavirus (TULV) within a single evolutionary lineage of their natural rodent host, the common vole Microtus arvalis. We examined voles from forty-two locations in the contact region between clades for TULV infection by reverse transcription (RT)-PCR. Sequencing yielded twenty-three TULV Central North and twenty-one TULV Central South genomes, which differed by 14.9-18.5 per cent at the nucleotide and 2.2-3.7 per cent at the amino acid (AA) level without evidence of recombination or reassortment between clades. Geographic cline analyses demonstrated an abrupt (<1 km wide) transition between the parapatric TULV clades in continuous landscape. This transition was located within the Central mitochondrial lineage of M. arvalis, and genomic single nucleotide polymorphisms showed gradual mixing of host populations across it. Genomic differentiation of hosts was much weaker across the TULV Central North to South transition than across the nearby hybrid zone between two evolutionary lineages in the host. We suggest that these parapatric TULV clades represent functionally distinct, incipient species, which are likely differently affected by genetic polymorphisms in the host. This highlights the potential of natural viral contact zones as systems for investigating the genetic and evolutionary factors enabling or restricting the transmission of RNA viruses.

Cocirculation of Leptospira spp. and multiple orthohantaviruses in rodents, Lithuania, Northern Europe

In Europe, zoonotic Leptospira spp. and orthohantaviruses are mainly associated with specific rodent hosts. These pathogens cause febrile human diseases with similar symptoms and disease progression. In Lithuania, the presence of Dobrava-Belgrade orthohantavirus (DOBV), Tula orthohantavirus (TULV) and Leptospira spp. in rodent reservoirs is still unknown, and Puumala orthohantavirus (PUUV) was detected in bank voles (Clethrionomys glareolus) at only one site. Therefore, we collected and screened 1617 rodents and insectivores from Lithuania for zoonotic (re-)emerging Leptospira and orthohantaviruses. We detected Leptospira DNA in six rodent species, namely striped field mouse (Apodemus agrarius), yellow-necked mouse (Apodemus flavicollis), bank vole, common vole (Microtus arvalis), field vole (Microtus agrestis) and root vole (Microtus oeconomus). Leptospira DNA was detected with an overall mean prevalence of 4.4% (range 3.7%-7.9% per rodent species). We detected DOBV RNA in 5.6% of the striped field mice, PUUV RNA in 1% of bank voles and TULV RNA in 4.6% of common voles, but no Leptospira DNA in shrews and no hantavirus-Leptospira coinfections in rodents. Based on the complete coding sequences of the three genome segments, two distant DOBV phylogenetic lineages in striped field mice, one PUUV strain in bank voles and two TULV strains in common voles were identified. The Leptospira prevalence for striped field mice and yellow-necked mice indicated a significant negative effect of the distance to water points. The detection of (re-)emerging human pathogenic Leptospira and three orthohantaviruses in rodent reservoirs in Lithuania calls for increased awareness of public health institutions and allows the improvement of molecular diagnostics for pathogen identification.

Molecular Characterisation and Phylogeny of Tula Virus in Kazakhstan

Orthohantaviruses are zoonotic pathogens that play a significant role in public health. These viruses can cause haemorrhagic fever with renal syndrome in Eurasia. In the Republic of Kazakhstan, the first human cases were registered in the year 2000 in the West Kazakhstan region. Small mammals can be reservoirs of orthohantaviruses. Previous studies showed orthohantavirus antigens in wild-living small mammals in four districts of West Kazakhstan. Clinical studies suggested that there might be further regions with human orthohantavirus infections in Kazakhstan, but genetic data of orthohantaviruses in natural foci are limited. The aim of this study was to investigate small mammals for the presence of orthohantaviruses by molecular biological methods and to provide a phylogenetic characterization of the circulating strains in Kazakhstan. Small mammals were trapped at 19 sites in West Kazakhstan, four in Almaty region and at seven sites around Almaty city during all seasons of 2018 and 2019. Lung tissues of small mammals were homogenized and RNA was extracted. Orthohantavirus RT-PCR assays were applied for detection of partial S and L segment sequences. Results were compared to published fragments. In total, 621 small mammals from 11 species were analysed. Among the collected small mammals, 2.4% tested positive for orthohantavirus RNA, one sample from West Kazakhstan and 14 samples from Almaty region. None of the rodents caught in Almaty city were infected. Sequencing parts of the small (S) and large (L) segments specified Tula virus (TULV) in these two regions. Our data show that geographical distribution of TULV is more extended as previously thought. The detected sequences were found to be split in two distinct genetic clusters of TULV in West Kazakhstan and Almaty region. TULV was detected in the common vole (Microtus arvalis) and for the first time in two individuals of the forest dormouse (Dryomys nitedula), interpreted as a spill-over infection in Kazakhstan.

Host genetic factors associated with the range limit of a European hantavirus

The natural host ranges of many viruses are restricted to very specific taxa. Little is known about the molecular barriers between species that lead to the establishment of this restriction or generally prevent virus emergence in new hosts. Here, we identify genomic polymorphisms in a natural rodent host associated with a strong genetic barrier to the transmission of European Tula orthohantavirus (TULV). We analysed the very abrupt spatial transition between two major phylogenetic clades in TULV across the comparatively much wider natural hybrid zone between evolutionary lineages of their reservoir host, the common vole (Microtus arvalis). Genomic scans of 79,225 single nucleotide polymorphisms (SNPs) in 323 TULV‐infected host individuals detected 30 SNPs that were consistently associated with the TULV clades CEN.S or EST.S in two replicate sampling transects. Focusing the analysis on 199 voles with evidence of genomic admixture at the individual level (0.1–0.9) supported statistical significance for all 30 loci. Host genomic variation at these SNPs explained up to 37.6% of clade‐specific TULV infections. Genes in the vicinity of associated SNPs include SAHH, ITCH and two members of the Syngr gene family, which are involved in functions related to immune response or membrane transport. This study demonstrates the relevance of natural hybrid zones as systems not only for studying processes of evolutionary divergence and speciation, but also for the detection of evolving genetic barriers for specialized parasites.

FREQUENT LEPTOSPIRA SPP. DETECTION BUT ABSENCE OF TULA ORTHOHANTAVIRUS IN MICROTUS SPP. VOLES, NORTHWESTERN SPAIN

The common vole (Microtus arvalis) is a major agricultural pest in Europe and is a reservoir for several zoonotic agents, such as Leptospira spp. and Tula orthohantavirus (TULV). However, little is known about the occurrence of those pathogens in voles from Spain, where the species has largely expanded its distribution range in the past decades, causing agricultural pests and zoonotic diseases. For a molecular survey, 580 common voles and six Lusitanian pine voles (Microtus lusitanicus) were collected in 26 localities from four provinces of northwestern Spain. We assessed the presence of Leptospira spp. DNA in kidney tissue by PCR targeting the lipL32 gene, detecting a prevalence of 7.9% (95% confidence interval, 5.9-10.4) for common voles and of 33.3% (95% confidence interval, 4.3-77.7) for Lusitanian pine voles. We identified Leptospira kirschneri in 24 animals and Leptospira borgpetersenii in two animals, using secY gene-specific PCR. We analyzed environmental and demographic factors (such as age class, weight, and sex) and population dynamics data for their potential effect on the Leptospira spp. prevalence in those voles. The Leptospira spp. DNA detection rate in common voles increased significantly with maximum air temperature, vole weight, and amount of accumulated rainfall during the 90 d before capture and within the peak phase of the population cycle. We assessed the presence of TULV in lung tissue of 389 voles by reverse-transcription PCR, with no positive results. The absence of TULV might be explained by the evolutionary isolation of the common vole in Spain. The detection of two Leptospira genomospecies underlines the necessity for further typing efforts to understand the epidemiology of leptospiral infection in the common vole and the potential risk for human health in Spain.

Geographical Distribution and Genetic Diversity of Bank Vole Hepaciviruses in Europe

The development of new diagnostic methods resulted in the discovery of novel hepaciviruses in wild populations of the bank vole (Myodes glareolus, syn. Clethrionomys glareolus). The naturally infected voles demonstrate signs of hepatitis similar to those induced by hepatitis C virus (HCV) in humans. The aim of the present research was to investigate the geographical distribution of bank vole-associated hepaciviruses (BvHVs) and their genetic diversity in Europe. Real-time reverse transcription polymerase chain reaction (RT-qPCR) screening revealed BvHV RNA in 442 out of 1838 (24.0%) bank voles from nine European countries and in one of seven northern red-backed voles (Myodes rutilus, syn. Clethrionomys rutilus). BvHV RNA was not found in any other small mammal species (n = 23) tested here. Phylogenetic and isolation-by-distance analyses confirmed the occurrence of both BvHV species (Hepacivirus F and Hepacivirus J) and their sympatric occurrence at several trapping sites in two countries. The broad geographical distribution of BvHVs across Europe was associated with their presence in bank voles of different evolutionary lineages. The extensive geographical distribution and high levels of genetic diversity of BvHVs, as well as the high population fluctuations of bank voles and occasional commensalism in some parts of Europe warrant future studies on the zoonotic potential of BvHVs.

Spatial and Temporal Dynamics and Molecular Evolution of Tula orthohantavirus in German Vole Populations

Tula orthohantavirus (TULV) is a rodent-borne hantavirus with broad geographical distribution in Europe. Its major reservoir is the common vole (Microtus arvalis), but TULV has also been detected in closely related vole species. Given the large distributional range and high amplitude population dynamics of common voles, this host-pathogen complex presents an ideal system to study the complex mechanisms of pathogen transmission in a wild rodent reservoir. We investigated the dynamics of TULV prevalence and the subsequent potential effects on the molecular evolution of TULV in common voles of the Central evolutionary lineage. Rodents were trapped for three years in four regions of Germany and samples were analyzed for the presence of TULV-reactive antibodies and TULV RNA with subsequent sequence determination. The results show that individual (sex) and population-level factors (abundance) of hosts were significant predictors of local TULV dynamics. At the large geographic scale, different phylogenetic TULV clades and an overall isolation-by-distance pattern in virus sequences were detected, while at the small scale (<4 km) this depended on the study area. In combination with an overall delayed density dependence, our results highlight that frequent, localized bottleneck events for the common vole and TULV do occur and can be offset by local recolonization dynamics.

Identification of a novel hantavirus strain in the root vole (Microtus oeconomus) in Lithuania, Eastern Europe

Hantaviruses are zoonotic pathogens that can cause subclinical to lethal infections in humans. In Europe, five orthohantaviruses are present in rodents: Myodes-associated Puumala orthohantavirus (PUUV), Microtus-associated Tula orthohantavirus, Traemmersee hantavirus (TRAV)/ Tatenale hantavirus (TATV)/ Kielder hantavirus, rat-borne Seoul orthohantavirus, and Apodemus-associated Dobrava-Belgrade orthohantavirus (DOBV). Human PUUV and DOBV infections were detected previously in Lithuania, but the presence of Microtus-associated hantaviruses is not known. For this study we screened 234 Microtus voles, including root voles (Microtus oeconomus), field voles (Microtus agrestis) and common voles (Microtus arvalis) from Lithuania for hantavirus infections. This initial screening was based on reverse transcription-polymerase chain reaction (RT-PCR) targeting the S segment and serological analysis. A novel hantavirus was detected in eight of 79 root voles tentatively named "Rusne virus" according to the capture location and complete genome sequences were determined. In the coding regions of all three genome segments, Rusne virus showed high sequence similarity to TRAV and TATV and clustered with Kielder hantavirus in phylogenetic analyses of partial S and L segment sequences. Pairwise evolutionary distance analysis confirmed Rusne virus as a strain of the species TRAV/TATV. Moreover, we synthesized the entire nucleocapsid (N) protein of Rusne virus in Saccharomyces cerevisiae. We observed cross-reactivity of antibodies raised against other hantaviruses, including PUUV, with this new N protein. ELISA investigation of all 234 voles detected Rusne virus-reactive antibodies exclusively in four of 79 root voles, all being also RNA positive, but not in any other vole species. In conclusion, the detection of Rusne virus RNA in multiple root voles at the same trapping site during three years and its absence in sympatric field voles suggests root voles as the reservoir host of this novel virus. Future investigations should evaluate host association of TRAV, TATV, Kielder virus and the novel Rusne virus and their evolutionary relationships.

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.

How Bank Vole-PUUV Interactions Influence the Eco-Evolutionary Processes Driving Nephropathia Epidemica Epidemiology-An Experimental and Genomic Approach

In Europe, Puumala virus (PUUV) is responsible for nephropathia epidemica (NE), a mild form of hemorrhagic fever with renal syndrome (HFRS). Despite the presence of its reservoir, the bank vole, on most of French territory, the geographic distribution of NE cases is heterogeneous and NE endemic and non-endemic areas have been reported. In this study we analyzed whether bank vole-PUUV interactions could partly shape these epidemiological differences. We performed crossed-experimental infections using wild bank voles from French endemic (Ardennes) and non-endemic (Loiret) areas and two French PUUV strains isolated from these areas. The serological response and dynamics of PUUV infection were compared between the four cross-infection combinations. Due to logistical constraints, this study was based on a small number of animals. Based on this experimental design, we saw a stronger serological response and presence of PUUV in excretory organs (bladder) in bank voles infected with the PUUV endemic strain. Moreover, the within-host viral diversity in excretory organs seemed to be higher than in other non-excretory organs for the NE endemic cross-infection but not for the NE non-endemic cross-infection. Despite the small number of rodents included, our results showed that genetically different PUUV strains and in a lesser extent their interaction with sympatric bank voles, could affect virus replication and diversity. This could impact PUUV excretion/transmission between rodents and to humans and in turn at least partly shape NE epidemiology in France.