Taking the beat of the Arctic: are lemming population cycles changing due to winter climate?

Reports of fading vole and lemming population cycles and persisting low populations in some parts of the Arctic have raised concerns about the spread of these fundamental changes to tundra food web dynamics. By compiling 24 unique time series of lemming population fluctuations across the circumpolar region, we show that virtually all populations displayed alternating periods of cyclic/non-cyclic fluctuations over the past four decades. Cyclic patterns were detected 55% of the time (n = 649 years pooled across sites) with a median periodicity of 3.7 years, and non-cyclic periods were not more frequent in recent years. Overall, there was an indication for a negative effect of warm spells occurring during the snow onset period of the preceding year on lemming abundance. However, winter duration or early winter climatic conditions did not differ on average between cyclic and non-cyclic periods. Analysis of the time series shows that there is presently no Arctic-wide collapse of lemming cycles, even though cycles have been sporadic at most sites during the last decades. Although non-stationary dynamics appears a common feature of lemming populations also in the past, continued warming in early winter may decrease the frequency of periodic irruptions with negative consequences for tundra ecosystems.

Trophic fate and biomagnification of organic micropollutants from staple food to a specialized predator.

The environmental burden of organic micropollutants has been shown in aquatic ecosystems, while trophic fate of many compounds in terrestrial food chains remains highly elusive. We therefore studied concentrations of 108 organic micropollutants in a common European mammal, the bank vole (Clethrionomys glareolus), and 82 of the compounds in a specialized predator, Tengmalm's owl (Aegolius funereus) relying to >90 % on voles as its prey. We studied compounds in whole voles (n = 19), pools of 4-8 bank voles (npools = 4), owl blood (n = 10) and in owl eggs (n = 10) in two regions in Sweden. For comparison, we also included previously published data on 23 PFAS (per- and polyfluoroalkyl substances) in bank vole liver (npools = 4) from the same regions. In voles, concentrations of the organic micropollutants caffeine (maxIndividual 220 ng/g ww) and DEET (N,N-diethyl-m-toluamide) (maxPool 150 ng/g ww) were 2-200 times higher in voles relative to owl blood and eggs. Conversely, concentrations of nicotine, oxazepam, salicylic acid, and tributyl citrate acetate were 1.3-440 times higher in owls. Several PFAS showed biomagnification in owls as revealed by maximum biomagnification factors (BMFs); PFNA (perfluorononanoate) BMF = 5.6, PFTeDA (perfluorotetradecanoic acid) BMF = 5.9, and PFOS (perfluorooctane sulfonate) BMF = 6.1. Concentrations of organic micropollutants, alongside calculated BMFs, and Tengmalm's owl's heavy reliance on bank vole as staple food, suggest, despite small sample size and potential spatio-temporal mismatch, accumulation of PFAS (especially PFNA, PFTeDA, and PFOS) in owls and biomagnification along the food chain. Concentrations of PFAS in owl eggs (e.g., 21 ng/g ww PFOS) highlight the likely pivotal role of maternal transfer in contaminant exposure for avian embryos. These concentrations are also of concern considering that certain predators frequently consume owl eggs, potentially leading to additional biomagnification of PFAS with yet undetermined consequences for ecosystem health.

Urban forest soils harbour distinct and more diverse communities of bacteria and fungi compared to less disturbed forest soils.

Anthropogenic changes to land use drive concomitant changes in biodiversity, including that of the soil microbiota. However, it is not clear how increasing intensity of human disturbance is reflected in the soil microbial communities. To address this issue, we used amplicon sequencing to quantify the microbiota (bacteria and fungi) in the soil of forests (n = 312) experiencing four different land uses, national parks (set aside for nature conservation), managed (for forestry purposes), suburban (on the border of an urban area) and urban (fully within a town or city), which broadly represent a gradient of anthropogenic disturbance. Alpha diversity of bacteria and fungi increased with increasing levels of anthropogenic disturbance, and was thus highest in urban forest soils and lowest in the national parks. The forest soil microbial communities were structured according to the level of anthropogenic disturbance, with a clear urban signature evident in both bacteria and fungi. Despite notable differences in community composition, there was little change in the predicted functional traits of urban bacteria. By contrast, urban soils exhibited a marked loss of ectomycorrhizal fungi. Soil pH was positively correlated with the level of disturbance, and thus was the strongest predictor of variation in alpha and beta diversity of forest soil communities, indicating a role of soil alkalinity in structuring urban soil microbial communities. Hence, our study shows how the properties of urban forest soils promote an increase in microbial diversity and a change in forest soil microbiota composition.

Global patterns in the metacommunity structuring of lake macrophytes: regional variations and driving factors.

We studied community-environment relationships of lake macrophytes at two metacommunity scales using data from 16 regions across the world. More specifically, we examined (a) whether the lake macrophyte communities respond similar to key local environmental factors, major climate variables and lake spatial locations in each of the regions (i.e., within-region approach) and (b) how well can explained variability in the community-environment relationships across multiple lake macrophyte metacommunities be accounted for by elevation range, spatial extent, latitude, longitude, and age of the oldest lake within each metacommunity (i.e., across-region approach). In the within-region approach, we employed partial redundancy analyses together with variation partitioning to investigate the relative importance of local variables, climate variables, and spatial location on lake macrophytes among the study regions. In the across-region approach, we used adjusted R2 values of the variation partitioning to model the community-environment relationships across multiple metacommunities using linear regression and commonality analysis. We found that niche filtering related to local lake-level environmental conditions was the dominant force structuring macrophytes within metacommunities. However, our results also revealed that elevation range associated with climate (increasing temperature amplitude affecting macrophytes) and spatial location (likely due to dispersal limitation) was important for macrophytes based on the findings of the across-metacommunities analysis. These findings suggest that different determinants influence macrophyte metacommunities within different regions, thus showing context dependency. Moreover, our study emphasized that the use of a single metacommunity scale gives incomplete information on the environmental features explaining variation in macrophyte communities.

Spatial prediction and validation of zoonotic hazard through micro-habitat properties: where does Puumala hantavirus hole - up?

BACKGROUND: To predict the risk of infectious diseases originating in wildlife, it is important to identify habitats that allow the co-occurrence of pathogens and their hosts. Puumala hantavirus (PUUV) is a directly-transmitted RNA virus that causes hemorrhagic fever in humans, and is carried and transmitted by the bank vole (Myodes glareolus). In northern Sweden, bank voles undergo 3-4 year population cycles, during which their spatial distribution varies greatly. METHODS: We used boosted regression trees; a technique inspired by machine learning, on a 10 - year time-series (fall 2003-2013) to develop a spatial predictive model assessing seasonal PUUV hazard using micro-habitat variables in a landscape heavily modified by forestry. We validated the models in an independent study area approx. 200 km away by predicting seasonal presence of infected bank voles in a five-year-period (2007-2010 and 2015). RESULTS: The distribution of PUUV-infected voles varied seasonally and inter-annually. In spring, micro-habitat variables related to cover and food availability in forests predicted both bank vole and infected bank vole presence. In fall, the presence of PUUV-infected voles was generally restricted to spruce forests where cover was abundant, despite the broad landscape distribution of bank voles in general. We hypothesize that the discrepancy in distribution between infected and uninfected hosts in fall, was related to higher survival of PUUV and/or PUUV-infected voles in the environment, especially where cover is plentiful. CONCLUSIONS: Moist and mesic old spruce forests, with abundant cover such as large holes and bilberry shrubs, also providing food, were most likely to harbor infected bank voles. The models developed using long-term and spatially extensive data can be extrapolated to other areas in northern Fennoscandia. To predict the hazard of directly transmitted zoonoses in areas with unknown risk status, models based on micro-habitat variables and developed through machine learning techniques in well-studied systems, could be used.

Sublethal Lead Exposure Alters Movement Behavior in Free-Ranging Golden Eagles.

Lead poisoning of animals due to ingestion of fragments from lead-based ammunition in carcasses and offal of shot wildlife is acknowledged globally and raises great concerns about potential behavioral effects leading to increased mortality risks. Lead levels in blood were correlated with progress of the moose hunting season. Based on analyses of tracking data, we found that even sublethal lead concentrations in blood (25 ppb, wet weight), can likely negatively affect movement behavior (flight height and movement rate) of free-ranging scavenging Golden Eagles (Aquila chrysaetos). Lead levels in liver of recovered post-mortem analyzed eagles suggested that sublethal exposure increases the risk of mortality in eagles. Such adverse effects on animals are probably common worldwide and across species, where game hunting with lead-based ammunition is widespread. Our study highlights lead exposure as a considerably more serious threat to wildlife conservation than previously realized and suggests implementation of bans of lead ammunition for hunting.

Average niche breadths of species in lake macrophyte communities respond to ecological gradients variably in four regions on two continents.

Different species' niche breadths in relation to ecological gradients are infrequently examined within the same study and, moreover, species niche breadths have rarely been averaged to account for variation in entire ecological communities. We investigated how average environmental niche breadths (climate, water quality and climate-water quality niches) in aquatic macrophyte communities are related to ecological gradients (latitude, longitude, altitude, species richness and lake area) among four distinct regions (Finland, Sweden and US states of Minnesota and Wisconsin) on two continents. We found that correlations between the three different measures of average niche breadths and ecological gradients varied considerably among the study regions, with average climate and average water quality niche breadth models often showing opposite trends. However, consistent patterns were also found, such as widening of average climate niche breadths and narrowing of average water quality niche breadths of aquatic macrophytes along increasing latitudinal and altitudinal gradients. This result suggests that macrophyte species are generalists in relation to temperature variations at higher latitudes and altitudes, whereas species in southern, lowland lakes are more specialised. In contrast, aquatic macrophytes growing in more southern nutrient-rich lakes were generalists in relation to water quality, while specialist species are adapted to low-productivity conditions and are found in highland lakes. Our results emphasise that species niche breadths should not be studied using only coarse-scale data of species distributions and corresponding environmental conditions, but that investigations on different kinds of niche breadths (e.g., climate vs. local niches) also require finer resolution data at broad spatial extents.

Selective predation on hantavirus-infected voles by owls and confounding effects from landscape properties.

It has been suggested that predators may protect human health through reducing disease-host densities or selectively preying on infected individuals from the population. However, this has not been tested empirically. We hypothesized that Tengmalm's owl (Aegolius funereus) selectively preys on hantavirus-infected individuals of its staple prey, the bank vole (Myodes glareolus). Bank voles are hosts of Puumala hantavirus, which causes a form of hemorrhagic fever in humans. Selective predation by owls on infected voles may reduce human disease risk. We compared the prevalence of anti-Puumala hantavirus antibodies (seroprevalence), in bank voles cached by owls in nest boxes to seroprevalence in voles trapped in closed-canopy forest around each nest box. We found no general difference in seroprevalence. Forest landscape structure could partly account for the observed patterns in seroprevalence. Only in more connected forest patches was seroprevalence in bank voles cached in nest boxes higher than seroprevalence in trapped voles. This effect disappeared with increasing forest patch isolation, as seroprevalence in trapped voles increased with forest patch isolation, but did not in cached voles. Our results suggest a complex relationship between zoonotic disease prevalence in hosts, their predators, and landscape structure. Some mechanisms that may have caused the seroprevalence patterns in our results include higher bank vole density in isolated forest patches. This study offers future research potential to shed further light on the contribution of predators and landscape properties to human health.

Impact of Beaver Pond Colonization History on Methylmercury Concentrations in Surface Water.

Elevated concentrations of methylmercury (MeHg) in freshwater ecosystems are of major environmental concern in large parts of the northern hemisphere. Beaver ponds have been identified as a potentially important source of MeHg. The role of beavers might be especially pronounced in large parts of Europe, where beaver populations have expanded rapidly following near-extirpation. This study evaluates the role of the age and colonization history (encompassing patterns of use and reuse) of ponds constructed by the Eurasian beaver Castor fiber in regulating MeHg concentrations in Swedish streams. In 12 beaver systems located in three regions, we quantified MeHg concentrations together with other relevant parameters on five occasions per year in 2012-2013. Five were pioneer systems, inundated for the first time since beaver extirpation, and seven were recolonized, with dams reconstructed by newly recolonizing beavers. MeHg concentrations in pioneer but not in recolonized beaver systems were up to 3.5 fold higher downstream than upstream of the ponds, and varied between seasons and years. Our results show that pioneer inundation by beavers can increase MeHg concentrations in streams, but that this effect is negligible when dams are reconstructed on previously used ponds. We therefore expect that the recovery and expansion of beavers in the boreal system will only have a transitional effect on MeHg in the environment.

Exploring the influence of host community composition on the outbreak potential of Anaplasma phagocytophilum and Borrelia burgdorferi s.l.

In large parts of the northern hemisphere, multiple deer species coexist, and management actions can strongly influence wild deer communities. Such changes may also indirectly influence other species in the community, such as small mammals and birds, because deer can have strong effects on their habitats and resources. Deer, small mammals and birds play an important role in the dynamics of tick-borne zoonotic diseases. It is, however, relatively underexplored how the abundance and composition of vertebrate communities may affect the outbreak potential, maintenance and circulation of tick-borne pathogens. In this study we focus on the outbreak potential by exploring how the basic reproduction number R0 for different tick-borne pathogens depends on host community composition. We used published data on co-varying roe deer (Capreolus capreolus) and fallow deer (Dama dama) densities following a hunting ban, and different small mammal and bird densities, to investigate how the change in host community influences the R0 of four tick-borne pathogens: one non-zoonotic, namely Anaplasma phagocytophilum ecotype 2, and three zoonotic, namely A. phagocytophilum ecotype 1, Borrelia afzelii and Borrelia garinii. We calculated R0 using a next generation matrix approach, and used elasticities to quantify the contributions to R0 of the different groups of host species. The value of R0 for A. phagocytophilum ecotype 1 was higher with high fallow deer density and low roe deer density, while it was the other way round for A. phagocytophilum ecotype 2. For B. afzelii, R0 was mostly related to the density of small mammals and for B. garinii it was mostly determined by bird density. Our results show that the effect of species composition is substantial in the outbreak potential of tick-borne pathogens. This implies that also management actions that change this composition, can (indirectly and unintentionally) affect the outbreak potential of tick-borne diseases.

Mystery of fatal 'staggering disease' unravelled: novel rustrela virus causes severe meningoencephalomyelitis in domestic cats.

'Staggering disease' is a neurological disease entity considered a threat to European domestic cats (Felis catus) for almost five decades. However, its aetiology has remained obscure. Rustrela virus (RusV), a relative of rubella virus, has recently been shown to be associated with encephalitis in a broad range of mammalian hosts. Here, we report the detection of RusV RNA and antigen by metagenomic sequencing, RT-qPCR, in-situ hybridization and immunohistochemistry in brain tissues of 27 out of 29 cats with non-suppurative meningoencephalomyelitis and clinical signs compatible with'staggering disease' from Sweden, Austria, and Germany, but not in non-affected control cats. Screening of possible reservoir hosts in Sweden revealed RusV infection in wood mice (Apodemus sylvaticus). Our work indicates that RusV is the long-sought cause of feline 'staggering disease'. Given its reported broad host spectrum and considerable geographic range, RusV may be the aetiological agent of neuropathologies in further mammals, possibly even including humans.

Puumala Orthohantavirus Infection Does Not Affect the Trapping Success of Its Reservoir Host.

Pathogens might affect behavior of infected reservoir hosts and hence their trappability, which could bias population estimates of pathogen prevalence. In this study, we used snap-trapping data on Puumala orthohantavirus (PUUV)-infected (n = 1619) and noninfected (n = 6940) bank voles (Myodes glareolus) from five vole cycles, normally representing increase, peak, and decline phase, to evaluate if infection status affected trapping success. If PUUV infection, as previously suggested, increases activity and/or mobility, we would expect a higher proportion of infected than noninfected specimens in the first trapping night. However, the proportion of PUUV-infected voles did not differ across the three trapping nights. We conclude that PUUV infection did not affect trapping success, confirming snap trapping as an appropriate trapping method for studies on PUUV prevalence and likely other orthohantaviruses.

Alphacoronavirus in a Daubenton's Myotis Bat (Myotis daubentonii) in Sweden.

The ongoing COVID-19 pandemic has stimulated a search for reservoirs and species potentially involved in back and forth transmission. Studies have postulated bats as one of the key reservoirs of coronaviruses (CoVs), and different CoVs have been detected in bats. So far, CoVs have not been found in bats in Sweden and we therefore tested whether they carry CoVs. In summer 2020, we sampled a total of 77 adult bats comprising 74 Myotis daubentonii, 2 Pipistrellus pygmaeus, and 1 M. mystacinus bats in southern Sweden. Blood, saliva and feces were sampled, processed and subjected to a virus next-generation sequencing target enrichment protocol. An Alphacoronavirus was detected and sequenced from feces of a M. daubentonii adult female bat. Phylogenetic analysis of the almost complete virus genome revealed a close relationship with Finnish and Danish strains. This was the first finding of a CoV in bats in Sweden, and bats may play a role in the transmission cycle of CoVs in Sweden. Focused and targeted surveillance of CoVs in bats is warranted, with consideration of potential conflicts between public health and nature conservation required as many bat species in Europe are threatened and protected.

Climate change accelerates winter transmission of a zoonotic pathogen.

Many zoonotic diseases are weather sensitive, raising concern how their distribution and outbreaks will be affected by climate change. At northern high latitudes, the effect of global warming on especially winter conditions is strong. By using long term monitoring data (1980-1986 and 2003-2013) from Northern Europe on temperature, precipitation, an endemic zoonotic pathogen (Puumala orthohantavirus, PUUV) and its reservoir host (the bank vole, Myodes glareolus), we show that early winters have become increasingly wet, with a knock-on effect on pathogen transmission in its reservoir host population. Further, our study is the first to show a climate change effect on an endemic northern zoonosis, that is not induced by increased host abundance or distribution, demonstrating that climate change can also alter transmission intensity within host populations. Our results suggest that rainy early winters accelerate PUUV transmission in bank voles in winter, likely increasing the human zoonotic risk in the North.

Population fluctuations and synanthropy explain transmission risk in rodent-borne zoonoses.

Population fluctuations are widespread across the animal kingdom, especially in the order Rodentia, which includes many globally important reservoir species for zoonotic pathogens. The implications of these fluctuations for zoonotic spillover remain poorly understood. Here, we report a global empirical analysis of data describing the linkages between habitat use, population fluctuations and zoonotic reservoir status in rodents. Our quantitative synthesis is based on data collated from papers and databases. We show that the magnitude of population fluctuations combined with species' synanthropy and degree of human exploitation together distinguish most rodent reservoirs at a global scale, a result that was consistent across all pathogen types and pathogen transmission modes. Our spatial analyses identified hotspots of high transmission risk, including regions where reservoir species dominate the rodent community. Beyond rodents, these generalities inform our understanding of how natural and anthropogenic factors interact to increase the risk of zoonotic spillover in a rapidly changing world.

Wild ungulate species differ in their contribution to the transmission of Ixodes ricinus-borne pathogens.

BACKGROUND: Several ungulate species are feeding and propagation hosts for the tick Ixodes ricinus as well as hosts to a wide range of zoonotic pathogens. Here, we focus on Anaplasma phagocytophilum and Borrelia burgdorferi (s.l.), two important pathogens for which ungulates are amplifying and dilution hosts, respectively. Ungulate management is one of the main tools to mitigate human health risks associated with these tick-borne pathogens. Across Europe, different species of ungulates are expanding their ranges and increasing in numbers. It is currently unclear if and how the relative contribution to the life-cycle of I. ricinus and the transmission cycles of tick-borne pathogens differ among these species. In this study, we aimed to identify these relative contributions for five European ungulate species. METHODS: We quantified the tick load and collected ticks and spleen samples from hunted fallow deer (Dama dama, n = 131), moose (Alces alces, n = 15), red deer (Cervus elaphus, n = 61), roe deer (Capreolus capreolus, n = 30) and wild boar (Sus scrofa, n = 87) in south-central Sweden. We investigated the presence of tick-borne pathogens in ticks and spleen samples using real-time PCR. We determined if ungulate species differed in tick load (prevalence and intensity) and in infection prevalence in their tissue as well as in the ticks feeding on them. RESULTS: Wild boar hosted fewer adult female ticks than any of the deer species, indicating that deer are more important as propagation hosts. Among the deer species, moose had the lowest number of female ticks, while there was no difference among the other deer species. Given the low number of infected nymphs, the relative contribution of all ungulate species to the transmission of B. burgdorferi (s.l.) was low. Fallow deer, red deer and roe deer contributed more to the transmission of A. phagocytophilum than wild boar. CONCLUSIONS: The ungulate species clearly differed in their role as a propagation host and in the transmission of B. burgdorferi and A. phagocytophilum. This study provides crucial information for ungulate management as a tool to mitigate zoonotic disease risk and argues for adapting management approaches to the local ungulate species composition and the pathogen(s) of concern.

Low-level environmental metal pollution is associated with altered gut microbiota of a wild rodent, the bank vole (Myodes glareolus).

Mining and related industries are a major source of metal pollution. In contrast to the well-studied effects of exposure to metals on animal physiology and health, the impacts of environmental metal pollution on the gut microbiota of wild animals are virtually unknown. As the gut microbiota is a key component of host health, it is important to understand whether metal pollution can alter wild animal gut microbiota composition. Using a combination of 16S rRNA amplicon sequencing and quantification of metal levels in kidneys, we assessed whether multi-metal exposure (the sum of normalized levels of fifteen metals) was associated with changes in gut microbiota of wild bank voles (Myodes glareolus) from two locations in Finland. Exposure to increased metal load was associated with higher gut microbiota species diversity (α-diversity) and altered community composition (ß-diversity), but not dispersion. Multi-metal exposure and increased levels of several metals (Cd, Hg, Pb and Se) were associated with differences in the abundance of microbial taxa, especially those within the families Clostridiales vadinBB60 group, Desulfovibrionaceae, Lachnospiraceae, Muribaculaceae and Ruminococcaceae. Our data indicate that even low-level metal pollution can affect the diversity of microbiota and be associated with deterministic differences in composition of host gut microbiota in wild animal populations. These findings highlight the need to study a broader range of metals and their cocktails that are more representative of the types of environmental exposure experienced by wild animals.

Prevalence of Orthohantavirus-Reactive Antibodies in Humans and Peri-Domestic Rodents in Northern Ethiopia.

In 2012, Tigray orthohantavirus was discovered in Ethiopia, but its seasonal infection in small mammals, and whether it poses a risk to humans was unknown. The occurrence of small mammals, rodents and shrews, in human inhabitations in northern Ethiopia is affected by season and presence of stone bunds. We sampled small mammals in two seasons from low- and high-density stone bund fields adjacent to houses and community-protected semi-natural habitats in Atsbi and Hagere Selam, where Tigray orthohantavirus was first discovered. We collected blood samples from both small mammals and residents using filter paper. The presence of orthohantavirus-reactive antibodies in blood was then analyzed using immunofluorescence assay (human samples) and enzyme linked immunosorbent assays (small mammal samples) with Puumala orthohantavirus as antigen. Viral RNA was detected by RT-PCR using small mammal blood samples. Total orthohantavirus prevalence (antibodies or virus RNA) in the small mammals was 3.37%. The positive animals were three Stenocephalemys albipes rats (prevalence in this species = 13.04%). The low prevalence made it impossible to determine whether season and stone bunds were associated with orthohantavirus prevalence in the small mammals. In humans, we report the first detection of orthohantavirus-reactive IgG antibodies in Ethiopia (seroprevalence = 5.26%). S. albipes lives in close proximity to humans, likely increasing the risk of zoonotic transmission.

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.