Evolutionary Relationships of Ljungan Virus Variants Circulating in Multi-Host Systems across Europe.

The picornavirus named 'Ljungan virus' (LV, species Parechovirus B) has been detected in a dozen small mammal species from across Europe, but detailed information on its genetic diversity and host specificity is lacking. Here, we analyze the evolutionary relationships of LV variants circulating in free-living mammal populations by comparing the phylogenetics of the VP1 region (encoding the capsid protein and associated with LV serotype) and the 3Dpol region (encoding the RNA polymerase) from 24 LV RNA-positive animals and a fragment of the 5' untranslated region (UTR) sequence (used for defining strains) in sympatric small mammals. We define three new VP1 genotypes: two in bank voles (Myodes glareolus) (genotype 8 from Finland, Sweden, France, and Italy, and genotype 9 from France and Italy) and one in field voles (Microtus arvalis) (genotype 7 from Finland). There are several other indications that LV variants are host-specific, at least in parts of their range. Our results suggest that LV evolution is rapid, ongoing and affected by genetic drift, purifying selection, spillover and host evolutionary history. Although recent studies suggest that LV does not have zoonotic potential, its widespread geographical and host distribution in natural populations of well-characterized small mammals could make it useful as a model for studying RNA virus evolution and transmission.

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.

Transmission ecology of taeniid larval cestodes in rodents in Sweden, a low endemic area for Echinococcus multilocularis.

Although local prevalence of Echinococcus multilocularis may be high, this zoonotic parasite has an overall low prevalence in foxes and rodents in Sweden. To better understand opportunities for E. multilocularis transmission in the Swedish environment, the aim of this study was to investigate other taeniid cestodes and to relate observed patterns to E. multilocularis. Cestode parasites were examined in fox feces and rodents caught in different habitats from four regions of Sweden. Arvicola amphibius and Microtus agrestis were parasitized with Versteria mustelae, Hydatigera taeniaeformis s. l., and E. multilocularis, whereas Myodes glareolus and Apodemus spp. were parasitized with V. mustelae, Taenia polyacantha, H. taeniaeformis s.l., and Mesocestoides spp. Rodents caught in field habitat (Ar. amphibius, Mi. agrestis) were more likely (OR 10, 95% CI 5-19) to be parasitized than rodents caught in forest habitat (My. glareolus, Apodemus spp.). The parasite preference for each rodent species was present regardless of the type of background contamination from fox feces. These results further support the importance of both ecological barriers and individual species susceptibility in parasite transmission, and indicate that future monitoring for E. multilocularis in the Swedish environment should focus in field habitats where Mi. agrestis and Ar. amphibius are abundant.

DISTRIBUTION AND SEASONAL VARIATION OF LJUNGAN VIRUS IN BANK VOLES (MYODES GLAREOLUS) IN FENNOSCANDIA.

Ljungan virus (LV) is a picornavirus originally isolated from Swedish bank voles ( Myodes glareolus ) in 1998. The association of LV with human disease has been debated ever since, but fundamental data on the ecology of the virus are still lacking. Here we present results of the first intensive study on the prevalence of LV in bank voles trapped in Fennoscandia (Sweden and Finland) from 2009-12 as determined by PCR. Using an LV-specific real-time reverse transcriptase PCR, LV was detected in the liver of 73 out of 452 (16.2%) individuals and in 13 out of 17 sampling sites across Sweden and Finland (mean per site prevalence 16%, SE 3%, range 0-50%). We found more infected animals in autumn compared to spring, and lighter and heavier individuals had a higher prevalence than those with intermediate body masses. The result that LV prevalence is also lower in heavier (i.e., older) animals suggests for the first time that LV infection is not persistent in rodents.

Islet autoantibodies present in association with Ljungan virus infection in bank voles (Myodes glareolus) in northern Sweden.

Bank voles are known reservoirs for Puumala hantavirus and probably also for Ljungan virus (LV), a suggested candidate parechovirus in type 1 diabetes etiology and pathogenesis. The aim of this study was to determine whether wild bank voles had been exposed to LV and if exposure associated to autoantibodies against insulin (IAA), glutamic acid decarboxylase 65 (GADA), or islet autoantigen-2 (IA-2A). Serum samples from bank voles (Myodes glareolus) captured in early summer or early winter of 1997 and 1998, respectively, were analyzed in radio binding assays for antibodies against Ljungan virus (LVA) and Puumala virus (PUUVA) as well as for IAA, GADA, and IA-2A. LVA was found in 25% (189/752), IAA in 2.5% (18/723), GADA in 2.6% (15/615), and IA-2A in 2.5% (11/461) of available bank vole samples. LVA correlated with both IAA (P = 0.007) and GADA (P < 0.001), but not with IA-2A (P = 0.999). There were no correlations with PUUVA, detected in 17% of the bank voles. Compared to LVA negative bank voles, LVA positive animals had higher levels of both IAA (P = 0.002) and GADA (P < 0.001), but not of IA-2A (P = 0.205). Levels of LVA as well as IAA and GADA were higher in samples from bank voles captured in early summer. In conclusion, LVA was detected in bank voles and correlated with both IAA and GADA but not with IA-2A. These observations suggest that exposure to LV may be associated with islet autoimmunity. It remains to be determined if islet autoantibody positive bank voles may develop diabetes in the wild. J. Med. Virol. 89:24-31, 2017. © 2016 Wiley Periodicals, Inc.

Support for targeted sampling of red fox (Vulpes vulpes) feces in Sweden: a method to improve the probability of finding Echinococcus multilocularis.

BACKGROUND: Localized concentrations of Echinococcus multilocularis eggs from feces of infected red fox (Vulpes vulpes) can create areas of higher transmission risk for rodent hosts and possibly also for humans; therefore, identification of these areas is important. However, in a low prevalence environment, such as Sweden, these areas could be easily overlooked. As part of a project investigating the role of different rodents in the epidemiology of E. multilocularis in Sweden, fox feces were collected seasonally from rodent trapping sites in two regions with known parasite status and in two regions with unknown parasite status, 2013-2015. The aim was to evaluate background contamination in rodent trapping sites from parasite eggs in these regions. To maximize the likelihood of finding fox feces positive for the parasite, fecal collection was focused in habitats with the assumed presence of suitable rodent intermediate hosts (i.e. targeted sampling). Parasite eggs were isolated from feces through sieving-flotation, and parasite species were then confirmed using PCR and sequencing. RESULTS: Most samples were collected in the late winter/early spring and in open fields where both Arvicola amphibius and Microtus agrestis were captured. Fox feces positive for E. multilocularis (41/714) were found within 1-3 field collection sites within each of the four regions. The overall proportion of positive samples was low (≤5.4%) in three regions, but was significantly higher in one region (22.5%, P < 0.001). There was not a significant difference between seasons or years. Compared to previous national screenings, our sampling strategy identified multiple E. multilocularis positive feces in all four regions, including the two regions with previously unknown parasite status. CONCLUSIONS: These results further suggest that the distribution of E. multilocularis is highly aggregated in the environment and provide support for further development of a targeted sampling strategy. Our results show that it was possible to identify new areas of high contamination in low endemic environments. After further elaboration, such a strategy may be particularly useful for countries designing surveillance to document freedom from disease.

First identification of Echinococcus multilocularis in rodent intermediate hosts in Sweden.

Echinococcus multilocularis is a zoonotic tapeworm with a sylvatic lifecycle and an expanding range in Europe. Monitoring efforts following its first identification in 2011 in Sweden have focused on the parasite's definitive host, the red fox (Vulpes vulpes). However, identifying rodent intermediate hosts is important to recognize opportunities for parasite transmission. During 2013-2015, livers from a total of 1566 rodents from four regions in Sweden were examined for E. multilocularis metacestode lesions. Species identity of suspect parasite lesions was confirmed by PCR and sequencing. E. multilocularis positive lesions >6 mm in diameter were also examined histologically. One Microtus agrestis out of 187 (0.5%, 95%CI: 0-2.9%), 8/439 (1.8%, 95%CI: 0.8-3.6%) Arvicola amphibius, 0/655 (0%, 95%CI: 0-0.6%) Myodes glareolus, and 0/285 (0%, 95%CI: 0-1.3%) Apodemus spp. contained E. multilocularis metacestode lesions. Presence of protoscoleces was confirmed in the infected M. agrestis and in three of eight infected A. amphibius. Six of the nine positive rodents were captured from the same field. This is the first report of E. multilocularis in intermediate hosts in Sweden. The cluster of positive rodents in one field shows that local parasite prevalence can be high in Sweden despite overall low national prevalence in foxes (<0.1%). The presence of protoscoleces in infected M. agrestis and A. amphibius indicate these species can serve as competent intermediate hosts in Sweden. However, their relative importance for E. multilocularis transmission in the Swedish environment is not yet possible to assess. In contrast, the negative findings in all M. glareolus and Apodemus spp. suggest that these species are of no importance.

Reappraisal of Hydatigera taeniaeformis (Batsch, 1786) (Cestoda: Taeniidae) sensu lato with description of Hydatigera kamiyai n. sp.

The common cat tapeworm Hydatigera taeniaeformis is a complex of three morphologically cryptic entities, which can be differentiated genetically. To clarify the biogeography and the host spectrum of the cryptic lineages, 150 specimens of H. taeniaeformis in various definitive and intermediate hosts from Eurasia, Africa and Australia were identified with DNA barcoding using partial mitochondrial cytochrome c oxidase subunit 1 gene sequences and compared with previously published data. Additional phylogenetic analyses of selected isolates were performed using nuclear DNA and mitochondrial genome sequences. Based on molecular data and morphological analysis, Hydatigera kamiyai n. sp. Iwaki is proposed for a cryptic lineage, which is predominantly northern Eurasian and uses mainly arvicoline rodents (voles) and mice of the genus Apodemus as intermediate hosts. Hydatigera taeniaeformis sensu stricto (s.s.) is restricted to murine rodents (rats and mice) as intermediate hosts. It probably originates from Asia but has spread worldwide. Despite remarkable genetic divergence between H. taeniaeformis s.s. and H. kamiyai, interspecific morphological differences are evident only in dimensions of rostellar hooks. The third cryptic lineage is closely related to H. kamiyai, but its taxonomic status remains unresolved due to limited morphological, molecular, biogeographical and ecological data. This Hydatigera sp. is confined to the Mediterranean and its intermediate hosts are unknown. Further studies are needed to classify Hydatigera sp. either as a distinct species or a variant of H. kamiyai. According to previously published limited data, all three entities occur in the Americas, probably due to human-mediated introductions.

Shaping zoonosis risk: landscape ecology vs. landscape attractiveness for people, the case of tick-borne encephalitis in Sweden.

BACKGROUND: In this paper, the hazard and exposure concepts from risk assessment are applied in an innovative approach to understand zoonotic disease risk. Hazard is here related to the landscape ecology determining where the hosts, vectors and pathogens are and, exposure is defined as the attractiveness and accessibility to hazardous areas. Tick-borne encephalitis in Sweden was used as a case study. METHODS: Three boosted regression tree models are compared: a hazard model, an exposure model and a global model which combines the two approaches. RESULTS: The global model offers the best predictive power and the most accurate modelling. The highest probabilities were found in easy-to-reach places with high landscape diversity, holiday houses, waterbodies and, well-connected forests of oak, birch or pine, with open-area in their ecotones, a complex shape, numerous clear-cuts and, a variation in tree height. CONCLUSION: While conditions for access and use of hazardous areas are quite specific to Scandinavia, this study offers promising perspectives to improve our understanding of the distribution of zoonotic and vector-borne diseases in diverse contexts.

Modelling zoonotic diseases in humans: comparison of methods for hantavirus in Sweden.

Because their distribution usually depends on the presence of more than one species, modelling zoonotic diseases in humans differs from modelling individual species distribution even though the data are similar in nature. Three approaches can be used to model spatial distributions recorded by points: based on presence/absence, presence/available or presence data. Here, we compared one or two of several existing methods for each of these approaches. Human cases of hantavirus infection reported by place of infection between 1991 and 1998 in Sweden were used as a case study. Puumala virus (PUUV), the most common hantavirus in Europe, circulates among bank voles (Myodes glareolus). In northern Sweden, it causes nephropathia epidemica (NE) in humans, a mild form of hemorrhagic fever with renal syndrome.Logistic binomial regression and boosted regression trees were used to model presence and absence data. Presence and available sites (where the disease may occur) were modelled using cross-validated logistic regression. Finally, the ecological niche model MaxEnt, based on presence-only data, was used.In our study, logistic regression had the best predictive power, followed by boosted regression trees, MaxEnt and cross-validated logistic regression. It is also the most statistically reliable but requires absence data. The cross-validated method partly avoids the issue of absence data but requires fastidious calculations. MaxEnt accounts for non-linear responses but the estimators can be complex. The advantages and disadvantages of each method are reviewed.

Hantaviruses and their hosts in Europe: reservoirs here and there, but not everywhere?

Five hantaviruses are known to circulate among rodents in Europe, and at least two among insectivores. Four (Dobrava, Saaremaa, Seoul, and Puumala [PUUV] viruses) are clearly associated with hemorrhagic fever with renal syndrome (HFRS). PUUV, the most common etiological agent of HFRS in Europe, is carried by the bank vole (Myodes glareolus), one of the most widespread and abundant mammal species in Europe. This host-virus system is among hantaviruses also the most studied one in Europe. However, HFRS incidence varies throughout the continent. The spatial as well as temporal variation in the occurrence of HFRS is linked to geographic differences in the population dynamics of the reservoir rodents in different biomes of Europe. While rodent abundance may follow mast seeding events in many parts of temperate Europe, in northern (N) Europe multiannual cycles in population density exist as the result of the interaction between rodent populations and specialist predator populations in a delayed density-dependent manner. The spatial distribution of hantaviruses further depends on parameters such as forest patch size and connectivity of the most suitable rodent habitats, and the conditions for the survival of the virus outside the host, as well as historical distribution patterns (phylogeographies) of hosts and viruses. In multiannually fluctuating populations of rodents, with population increases of great amplitude, one should expect a simultaneous build-up of recently hantavirus-infected (shedding) rodents. The increasing number of infectious, virus-shedding rodents leads to a rapid transmission of hantavirus across the rodent population, and to humans. Our review discusses these aspects for PUUV, the only European hantavirus for which there is a reasonable, yet still far from complete, ecological continental-wide understanding. We discuss how this information could translate to other European hantavirus-host systems, and where the most important questions lie for further research.

Tnf-α expression and promoter sequences reflect the balance of tolerance/resistance to Puumala hantavirus infection in European bank vole populations.

The tumor necrosis factor-alpha (TNF-α) influences the ability to limit parasite infection but its over-production might result in inflammatory disorders. The level of Tnf-α gene expression could thus mediate a balance of tolerance/resistance to infections. This study focused on Puumala hantavirus (PUUV) infection in its rodent host, the bank vole (Myodes glareolus). In humans, PUUV is responsible of a mild form of hemorrhagic fever with renal syndrome, nephropathia epidemica (NE). The severity of NE is associated with an over-production of TNF-α. By contrast, PUUV infection in bank vole is chronic and asymptomatic. It is likely that different coevolutionary histories between PUUV and its hosts could lead to different balances of resistance/tolerance to PUUV infection, at least partly mediated by variable production levels of TNF-α. We investigated the hypothesis that bank voles from PUUV endemic areas should exhibit higher levels of tolerance, i.e. lower levels of TNF-α production, than bank voles from areas where PUUV prevalence is low. For this purpose, we analysed variations of Tnf-α gene expression and promoter sequences among European populations of bank voles. Our results revealed an absence of up-regulation of Tnf-α gene expression in PUUV infected bank voles and significant differences in Tnf-α gene expression level with regard to PUUV endemicity. These results corroborated the hypothesis of different balances of tolerance/resistance to PUUV. Two single-nucleotide polymorphism genotypes within the Tnf-α promoter (-302 GG/GG and -296 A/A) were associated with higher Tnf-α gene expression and were more frequent in non-endemic areas. This study emphasized the potential influence of selection acting on TNF-α production and mediating a tolerance/resistance balance to PUUV in bank voles. Further investigations, including the role of phenotypic plasticity and parasite communities on Tnf-α expression levels, should provide important keys to understand the prevalence of PUUV over Europe.

Associations between MHC genes and Puumala virus infection in Myodes glareolus are detected in wild populations, but not from experimental infection data.

We analysed the influence of MHC class II Dqa and Drb genes on Puumala virus (PUUV) infection in bank voles (Myodes glareolus). We considered voles sampled in five European localities or derived from a previous experiment that showed variable infection success of PUUV. The genetic variation observed in the Dqa and Drb genes was assessed by using single-strand conformation polymorphism and pyrosequencing methods, respectively. Patterns were compared with those obtained from 13 microsatellites. We revealed significant genetic differentiation between PUUV-seronegative and -seropositive bank voles sampled in wild populations, at the Drb gene only. The absence of genetic differentiation observed at neutral microsatellites confirmed the important role of selective pressures in shaping these Drb patterns. Also, we found no significant associations between infection success and MHC alleles among laboratory-colonized bank voles, which is explained by a loss of genetic variability that occurred during the captivity of these voles.

Puumala hantavirus and Myodes glareolus in northern Europe: no evidence of co-divergence between genetic lineages of virus and host.

The genus Hantavirus (family Bunyaviridae) includes negative-strand RNA viruses that are carried by persistently infected rodent and insectivore species. Puumala virus (PUUV), carried by bank voles (Myodes glareolus), is a pathogenic hantavirus that causes outbreaks of mild haemorrhagic fever with renal syndrome across Europe. In northern Europe, PUUV is represented by several genetic lineages that are maintained by distinct phylogroups of bank voles. The present study describes sequences of new PUUV strains recovered from northern and southern regions of Scandinavia and compares phylogenetic relationships between north-European PUUV strains and M. glareolus. This analysis revealed contradictions in phylogenetic clustering and remarkable differences in estimated divergence times between the lineages of PUUV and its host, suggesting that the established PUUV lineages did not co-diverge with the distinct phylogroups of M. glareolus that carry them at present.

Predicting high risk for human hantavirus infections, Sweden.

An increased risk for hemorrhagic fever with renal syndrome caused by Puumala hantavirus was forecast for Sweden in 2007. The forecast was based on a predicted increase in the number of Myodes glareolus rodents (reservoir hosts). Despite raised awareness and preparedness, the number of human cases during July 2007-June 2008 was 1,483, a new high.

Puumala hantavirus genetic variability in an endemic region (Northern Sweden).

Puumala hantavirus (PUUV), naturally harboured and shed by bank voles (Myodes [Clethrionomys] glareolus), is the etiological agent to nephropathia epidemica (NE), a mild haemorrhagic fever with renal syndrome. Both host and virus are found throughout much of the European continent and in northern Sweden NE is the second most prevalent serious febrile viral infection after influenza. The reliability of diagnostics by PCR depends on genetic variability for the detection of viral nucleic acids in unknown samples. In the present study we evaluated the genetic variability of PUUV isolated from bank voles in an area of northern Sweden highly endemic for NE. Genetic variability among bank voles was also investigated to evaluate co-evolutionary patterns. We found that the viral sequence appeared stable across the 80km study region, with the exception of the southernmost sampling site, which differed from its nearest neighbour by 7%, despite a geographical separation of only 10km. The southernmost sampling site demonstrated a higher degree of genetic similarity to PUUV previously isolated 100km south thereof; two locations appear to constitute a separate PUUV phylogenetic branch. In contrast to the viral genome, no phylogenetic variance was observed in the bank vole mtDNA in this study. Previous studies have shown that as a result of terrestrial mammals' postglacial re-colonization routes, bank voles and associated PUUV of a southern and a northern lineage established a dichotomous contact zone across the Scandinavian peninsula approximately 100-150km south of the present study sites. Our observations reveal evolutionary divergence of PUUV that has led to dissimilarities within the restricted geographical scale of the northern host re-colonization route as well. These results suggest either a static situation in which PUUV strains are regionally well adapted, or an ongoing process in which strains of PUUV circulate on a geographical scale not yet reliably described.

Puumala hantavirus viremia diagnosed by real-time reverse transcriptase PCR using samples from patients with hemorrhagic fever and renal syndrome.

Puumala virus (PUUV) is the endemic hantavirus in northern Sweden and causes nephropathia epidemica (NE), a milder form of hemorrhagic fever with renal syndrome. There is a need for fast and reliable diagnostics to differentiate the disease from other infections. By aligning virus RNA sequences isolated from 11 different bank voles and one human patient, we designed a real-time reverse transcriptase (RT) PCR method for detection of PUUV RNA. The real-time RT-PCR assay showed linearity from 20 to 2 x 10(6) virus copies with a correlation coefficient above 0.98 to 0.99 for all experiments. The detection threshold for PUUV cDNA was two copies per reaction. A two-step qualitative RT-PCR to detect PUUV RNA showed 100% concordance with the real-time RT-PCR assay. PUUV RNA viremia was detected in 33 of 34 PUUV immunoglobulin M (IgM)-positive patients with typical clinical NE disease from the region of endemicity. One PUUV IgM-negative sample had PUUV RNA, and 4 days later, the patient was IgM positive. Of samples with indeterminate IgM, 43% were PUUV RNA positive. The kinetics of antibody titers and PUUV viremia were studied, and five of six NE patients displayed a decrease in PUUV viremia a few days after disease outbreak coupled with an increase in PUUV IgM and IgG. In one patient with continuously high PUUV RNA levels but low IgM and no IgG response, the infection was lethal. These findings demonstrated that real-time RT-PCR is a useful method for diagnosis of PUUV viremia and for detecting PUUV RNA at early time points, before the appearance of IgM antibodies.

Habitat factors associated with bank voles (Clethrionomys glareolus) and concomitant hantavirus in northern Sweden.

Puumala virus (PUUV), genus hantavirus, causes nephropathia epidemica, a mild form of hemorrhagic fever with renal syndrome in humans. In this study, bank voles, the natural reservoir of PUUV, were captured at locations of previous human PUUV exposure and paired controls within a region of high incidence in northern Sweden. The aim of the study was to evaluate the influence of environmental factors on the abundance of bank voles and the occurrence of PUUV. The total number of voles and the number of PUUV-infected voles did not differ between locations of previous human PUUV exposure and paired controls. The number of bank voles expressing antibodies to PUUV infection increased linearly with total bank vole abundance implying density independent transmission. Using principal component and partial correlation analysis, we found that particular environmental characteristics associated with old-growth moist forests (i.e., those dominated by Alectoria spp., Picea abies, fallen wood, and Vaccinium myrtillus) were also associated with increased abundance of bank vole and hence the number of PUUV-infected bank voles, whereas there were no correlations with factors associated with dry environments (i.e., Pinus sylvestris and V. vitis-idea). This suggests that circulation and persistence of PUUV within bank vole populations was influenced by habitat factors. Future modeling of risk of exposure to hantavirus and transmission of PUUV within vole populations should include the influence of these factors.