==== Front BMC Vet Res BMC Vet Res BMC Veterinary Research 1746-6148 BioMed Central London 2702 10.1186/s12917-020-02702-y Research Article Active surveillance for antibodies confirms circulation of lyssaviruses in Palearctic bats Seidlova Veronika seidlovav@vfu.cz 1 Zukal Jan 23 Brichta Jiri 1 Anisimov Nikolay 4 Apoznański Grzegorz 5 Bandouchova Hana 1 Bartonička Tomáš 3 Berková Hana 2 Botvinkin Alexander D. 6 Heger Tomas 1 Dundarova Heliana 7 Kokurewicz Tomasz 5 Linhart Petr 1 Orlov Oleg L. 48 Piacek Vladimir 1 Presetnik Primož 9 Shumkina Alexandra P. 10 Tiunov Mikhail P. 11 Treml Frantisek 12 Pikula Jiri 113 1 grid.412968.00000 0001 1009 2154Department of Ecology and Diseases of Game, Fish and Bees, University of Veterinary and Pharmaceutical Sciences Brno, Palackého tř. 1946/1, 612 42 Brno, Czech Republic 2 grid.418095.10000 0001 1015 3316Institute of Vertebrate Biology, Czech Academy of Sciences, Květná 8, 603 65 Brno, Czech Republic 3 grid.10267.320000 0001 2194 0956Department of Botany and Zoology, Masaryk University, Kotlářská 267/2, 611 37 Brno, Czech Republic 4 grid.446209.d0000 0000 9203 3563Land Use and Biodiversity, International Complex Research Laboratory for Study of Climate Change, Tyumen State University, Volodarckogo 6, 625003 Tyumen, Russia 5 grid.411200.60000 0001 0694 6014Institute of Biology, Department of Vertebrate Ecology and Palaeontology, Wrocław University of Environmental and Life Sciences, Wrocław, Poland 6 grid.446313.70000 0001 0451 2298Irkutsk State Medical University, Krasnogo Vosstania street 1, 664003 Irkutsk, Russian Federation 7 grid.424727.00000 0004 0582 9037Department of Ecosystem Research, Environment Risk Assessment and Conservation Biology, Institute of Biodiversity and Ecosystem Research, Tsar Osvoboditel 1, 1000 Sofia, Bulgaria 8 grid.467075.70000 0004 0480 6706Department of Biochemistry, Ural State Medical University, Repina 3, 620014 Ekaterinburg, Russia 9 grid.433381.eCentre for Cartography of Fauna and Flora, Antoličičeva 1, SI-2204 Miklavž na Dravskem polju, Slovenia 10 Western Baikal protected areas, Federal State Budgetary Institution “Zapovednoe Pribaikalye”, Baikalskaya st. 291B, 664050 Irkutsk, Russia 11 grid.417808.20000 0001 1393 1398Institute of Biology and Soil Science, Far East Branch of the Russian Academy of Sciences, Pr- t 100-letiya Vladivostoka 159, 690022 Vladivostok, Russia 12 grid.412968.00000 0001 1009 2154Department of Infectious Diseases and Microbiology, University of Veterinary and Pharmaceutical Sciences Brno, Palackého tř. 1946/1, 612 42 Brno, Czech Republic 13 grid.412968.00000 0001 1009 2154CEITEC - Central European Institute of Technology, University of Veterinary and Pharmaceutical Sciences Brno, Brno, Czech Republic 10 12 2020 10 12 2020 2020 16 48227 11 2019 2 12 2020 © The Author(s) 2020Open AccessThis article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.Background Palearctic bats host a diversity of lyssaviruses, though not the classical rabies virus (RABV). As surveillance for bat rabies over the Palearctic area covering Central and Eastern Europe and Siberian regions of Russia has been irregular, we lack data on geographic and seasonal patterns of the infection. Results To address this, we undertook serological testing, using non-lethally sampled blood, on 1027 bats of 25 species in Bulgaria, the Czech Republic, Poland, Russia and Slovenia between 2014 and 2018. The indirect enzyme-linked immunosorbent assay (ELISA) detected rabies virus anti-glycoprotein antibodies in 33 bats, giving an overall seroprevalence of 3.2%. Bat species exceeding the seroconversion threshold included Myotis blythii, Myotis gracilis, Myotis petax, Myotis myotis, Murina hilgendorfi, Rhinolophus ferrumequinum and Vespertilio murinus. While Myotis species (84.8%) and adult females (48.5%) dominated in seropositive bats, juveniles of both sexes showed no difference in seroprevalence. Higher numbers tested positive when sampled during the active season (10.5%), as compared with the hibernation period (0.9%). Bat rabies seroprevalence was significantly higher in natural habitats (4.0%) compared with synanthropic roosts (1.2%). Importantly, in 2018, we recorded 73.1% seroprevalence in a cave containing a M. blythii maternity colony in the Altai Krai of Russia. Conclusions Identification of such “hotspots” of non-RABV lyssavirus circulation not only provides important information for public health protection, it can also guide research activities aimed at more in-depth bat rabies studies. Keywords Chiropterarabiesblood samplesseroprevalenceEuropeSiberiahttp://dx.doi.org/10.13039/501100001824Grantová Agentura České Republiky17-20286Shttp://dx.doi.org/10.13039/501100003335Ministry of Education, Youth and ScienceGrant No. KП-06-H31/14issue-copyright-statement© The Author(s) 2020 ==== Body Background Lyssaviruses are zoonotic agents of rabies that cause fatal encephalomyelitis in mammals. Different bat species act as principal reservoirs for most lyssaviruses, though carnivores only host the type species rabies virus (RABV), which is responsible for the majority of human rabies cases [1]. While dog-bite mediated rabies can be eliminated by control measures such as obligatory animal vaccination and pre- and post-exposure prophylactic treatment, an estimated 59 000 people still die annually from rabies in underdeveloped countries [2]. In both Europe and North America, however, large-scale oral rabies vaccination campaigns have been successful in eliminating the risk of exposure to rabid wild carnivores [3–6]. Unlike New World insectivorous bats, Palearctic bats only host non-RABV lyssaviruses. Some of these are only known from single or few isolates and have been associated with either no or only sporadic human rabies cases contracted via bat bite [7–16]. Mass vaccination of reservoir populations is presently unfeasible as a control and elimination strategy for bat rabies; not only as culling of rabies-positive bat colonies runs counter to present international legislation regarding conservation of threatened species, especially EC Directive 92/43/EEC of 21 May 1992 on the Conservation of Natural Habitats and of Wild Fauna and Flora, but also as this approach may have the opposite effect on rabies epidemiology by stimulating bat dispersal [17]. Accepting that widespread endemicity and persistence of bat lyssaviruses cannot be prevented, surveillance may prove the best option for risk evaluation and public health protection [16] [18]. Two protocols presently exist for screening bat lyssavirus infection. The first is passive surveillance, which involves testing dead bats at roosting sites or close to human habitation. Laboratory submissions also include diseased bats suspected of having rabies, those that die in rescue centres or those that have injured humans [15] [19] [20]. This protocol relies on the vigilance of both the public and bat specialists in order for wildlife casualties to be reported and presented for examination. The second protocol involves active surveillance of live bats that are captured and sampled using non-lethal methods [10] [15]. Serological screening is the main method of active surveillance for bat rabies. Bat rabies surveillance activities in different countries of the Palearctic region and bat species are irregular, with decreasing intensity from west to east [10] [15] [21–26]. Bats sampled using these protocols fall into two non-overlapping groups, i.e. diseased and healthy bats, though bats from both groups may be possible carriers of lyssaviruses. Possible sources of bias, however, include the fact that readily encountered synanthropic species tend to prevail among species presented for examination and that both protocols target bats in the active season of their annual life cycle. For several years, our group has been collecting blood for the study of host-pathogen interactions between hibernating Palearctic bats and the white-nose syndrome fungus Pseudogymnoascus destructans [27–37]. As such, we have hundreds of stored blood samples that provide an opportunity for screening rabies virus anti-glycoprotein antibodies. Given the varying intensity of surveillance for bat rabies over an extensive area of the Palearctic, we hypothesise that (i) circulation of lyssaviruses has previously gone undetected and, consequently, infection is underreported in some countries, and (ii) seroprevalence between summer- and winter-sampled bats will differ. Results Bats showed no clinical signs suggestive of rabies encephalomyelitis during the short period (hours) of sampling in the field or weeks of observation in the rescue centres. With 33 positive individuals, the overall rabies seroprevalence in our Palearctic collection was 3.2% (Tables 1 and 2). Bat species exceeding the seroconversion threshold of 0.123 EU/ml (derived from a calibration curve) included R. ferrumequinum in Bulgaria, M. myotis in the Czech Republic and Poland and N. noctula in the Czech Republic. The highest diversity of seropositive bat species, including M. blythii, M. gracilis, M. petax, M. hilgendorfi and V. murinus, was found in the Altai Krai of Russia. Antibody titres were mostly in the range of insufficient seroconversion, with only three bats (two M. myotis and one M. blythii) testing close to the threshold of sufficient seroconversion. A further M. blythii female was recorded over this threshold and had the highest measured titre of 0.805 EU/ml. Aside from the Bat’s Cave (Peshchera letuchikh myshey, Chineta reserve) in Altai Krai, Russia, where 26 bats of three species were examined and 19 individuals of two species were confirmed positive (M. blythii and M. petax) with a seroprevalence of 73.1%, seropositive bats tended to occur only sporadically. Antibody detection revealed three bat species newly identified as a seropositive, i.e. M. hilgendorfi, M. petax and M. gracilis. Myotis species (84.8%) and females (48.5%) were dominant among seropositive adult bats (χ2 = 10.07; p = 0.002), while the percentage of seropositive juvenile bats was the same for both sexes (χ2 = 0.01; p = 0.943). A higher percentage of bats tested positive during the active season (10.5%) than the hibernation period (0.9%), and rabies antibody prevalence was significantly higher in natural habitats (4.0%) than in synanthropic roosts (1.2%; p = 0.034). Table 1 Summary of geographic sampling details and rabies infection status (seropositivity). BG = Bulgaria, CZ = Czech Republic, PL = Poland, RU = Russia, SLO = Slovenia, N.A. = not available Country Sample site Habitat Number examined Species examined Gender data (F/M) Age data (AD/SAD/JUV) Seropositive bats BG Rakitovo (Lepenitza cave) natural 51 4 12/39 51/0/0 0 BG Stoychovtsi (Andaka cave) natural 26 1 16/10 26/0/0 1 CZ city of Brno (building) synanthropic 10 1 8/1 (1 N.A.) N.A. 0 CZ Ivanovice na Hané (building) synanthropic 32 1 22/10 32/0/0 0 CZ Jinačovice (rescue centre) synanthropic 10 3 0/4 (6 N.A.) 0/4 (6 N.A.) 0 CZ Mořina (Malá Amerika mines) natural 41 2 11/30 32/5/4 0 CZ Malá Morávka (Šimon and Juda mines) natural 99 2 31/66 (2 N.A.) 71/20/6 (2 N.A.) 1 CZ Adamov (Býčí skála cave) natural 40 1 14/26 31/3/6 0 CZ Suchdol (Kateřinská cave) natural 68 7 12/49 (7 N.A.) 51/4/5 (8 N.A.) 0 CZ Sloup (Sloupsko-Šošůvske caves) natural 52 1 16/36 41/8/3 0 CZ city of Plzeň (reconstructed building) synanthropic 68 1 53/15 68/0/0 0 CZ Poděbrady (rescue centre) synanthropic 5 1 4/1 5/0/0 1 CZ Vranov nad Dyjí (Ledové Sluje caves) natural 149 12 36/84 (29 N.A.) 8/2 (139 N.A.) 1 PL Nietoperek (bunker) synanthropic 110 3 51/59 100/9/1 2 RU Tigirek (Cave Jashchur) natural 11 5 6/5 8/3 1 RU Tigirek (Strashnaya cave) natural 1 1 1/0 1/0/0 0 RU Tigirek (netting) natural 9 3 1/8 7/0/2 4 RU Sakhyurta (Mechta cave) natural 30 2 8/22 30/0/0 0 RU Yelansy (Vologodskovo cave) natural 9 1 1/8 9/0/0 0 RU Ust-Pustynka (Peshchera letuchikh myshey) natural 26 3 18/8 12/0/14 19 RU Kozevnikovo (steppe lakes) natural 1 1 1/0 1/0/0 1 RU Spirino (Ob river) natural 3 2 2/1 1/0/1 2 RU Arakaevo (Arakaevskaya cave) natural 10 2 5/5 10/0/0 0 RU Pokrovskoje (Smolinskaya cave) natural 21 3 4/17 21/0/0 0 RU Severouralsk (Dačnaya cave) natural 1 1 0/1 1/0/0 0 RU Severouralsk (Komsomolskaya cave) natural 17 3 9/7 (1 N.A.) 17/0/0 0 RU Aramashevo (Sharkanskaya cave) natural 25 4 12/13 25/0/0 0 RU Severouralsk (Partizanskaya cave) natural 14 2 12/2 14/0/0 0 RU Cheremuchovo (Chertovo gorodishte) natural 19 3 7/12 9/10/0 0 RU Sergeyevka (Primorskiy Velikan cave) natural 32 3 15/17 19/13/0 0 RU Barsukovo (Barsukovskaya cave) natural 20 2 7/13 14/6/0 0 SLO Lovrenc na Pohorju (church Devica Maria) synanthropic 17 1 10/7 17/0/0 0 Total 1027 26 405/576 (46 N.A.) 737/83/41 (165) 33 Table 2 Characteristics of rabies seropositive bats. BG = Bulgaria, CZ = Czech Republic, PL = Poland, RU = Russia, SLO = Slovenia, N.A. = not available Bat species Country Sample site Sampling date Gender Age Titre (EU/ml) Rhinolophus ferrumequinum BG Rakitovo (Lepenitza cave) 01.04.2017 female adult 0.141 Myotis myotis CZ Malá Morávka (Šimon and Juda mines) 15.04.2014 N.A. N.A. 0.499 Myotis myotis CZ Vranov nad Dyjí (Ledové Sluje caves) 10.04.2015 male adult 0.147 Nyctalus noctula CZ Poděbrady (rescue centre) 07.12.2015 female adult 0.191 Myotis myotis PL Nietoperek (bunker) 19.03.2016 female sub-adult 0.125 Myotis myotis PL Nietoperek (bunker) 20.03.2016 female adult 0.473 Murina hilgendorfi RU Tigirek (Cave Jashchur) 15.04.2018 male adult 0.440 Vespertilio murinus RU Kozevnikovo (steppe lakes) 05.08.2018 female adult 0.235 Vespertilio murinus RU Spirino (Ob river) 07.08.2018 female juvenile 0.156 Myotis gracilis RU Spirino (Ob river) 07.08.2018 male adult 0.271 Myotis blythii RU Ust-Pustynka (Peshchera letuchikh myshey) 13.08.2018 female adult 0.136 Myotis blythii RU Ust-Pustynka (Peshchera letuchikh myshey) 13.08.2018 female adult 0.123 Myotis blythii RU Ust-Pustynka (Peshchera letuchikh myshey) 13.08.2018 male juvenile 0.321 Myotis blythii RU Ust-Pustynka (Peshchera letuchikh myshey) 13.08.2018 female juvenile 0.176 Myotis blythii RU Ust-Pustynka (Peshchera letuchikh myshey) 13.08.2018 female adult 0.370 Myotis blythii RU Ust-Pustynka (Peshchera letuchikh myshey) 13.08.2018 male juvenile 0.174 Myotis blythii RU Ust-Pustynka (Peshchera letuchikh myshey) 13.08.2018 female adult 0.805 Myotis blythii RU Ust-Pustynka (Peshchera letuchikh myshey) 13.08.2018 female adult 0.163 Myotis blythii RU Ust-Pustynka (Peshchera letuchikh myshey) 13.08.2018 female adult 0.233 Myotis blythii RU Ust-Pustynka (Peshchera letuchikh myshey) 13.08.2018 female adult 0.198 Myotis blythii RU Ust-Pustynka (Peshchera letuchikh myshey) 13.08.2018 female adult 0.451 Myotis blythii RU Ust-Pustynka (Peshchera letuchikh myshey) 13.08.2018 female juvenile 0.169 Myotis blythii RU Ust-Pustynka (Peshchera letuchikh myshey) 13.08.2018 female adult 0.328 Myotis petax RU Tigirek (netting) 10.08.2018 male adult 0.145 Myotis petax RU Tigirek (netting) 10.08.2018 male adult 0.139 Myotis petax RU Tigirek (netting) 12.08.2018 male juvenile 0.235 Myotis petax RU Tigirek (netting) 12.08.2018 male juvenile 0.163 Myotis petax RU Ust-Pustynka (Peshchera letuchikh myshey) 13.08.2018 male juvenile 0.323 Myotis petax RU Ust-Pustynka (Peshchera letuchikh myshey) 13.08.2018 female adult 0.229 Myotis petax RU Ust-Pustynka (Peshchera letuchikh myshey) 13.08.2018 male juvenile 0.139 Myotis petax RU Ust-Pustynka (Peshchera letuchikh myshey) 13.08.2018 female adult 0.187 Myotis petax RU Ust-Pustynka (Peshchera letuchikh myshey) 13.08.2018 male juvenile 0.277 Myotis petax RU Ust-Pustynka (Peshchera letuchikh myshey) 13.08.2018 female adult 0.180 Discussion In this study, we were able to extend active surveillance of bat rabies to selected countries of Central and Eastern Europe and eastern Palearctic (Siberian and Far East) regions of Russia. Despite the fact that the bats were primarily sampled for a different purpose, the seroprevalence data obtained suggests both ongoing exposure of Palearctic bats to lyssavirus circulation and spatial and temporal variation in bat rabies infection. Rabies invariably proves fatal for mammals once clinical signs develop. On the other hand, much is still poorly understood regarding bats infected with non-RABV lyssaviruses, e.g. persistence in natural hosts, pathogenesis and epidemiology [15] [38]. For example, some bats survive natural exposure to non-RABV lyssaviruses and go on to produce post-infection antibodies [15] [21–24] [39], and yet little is known about antiviral functions of the chiropteran immune system [40–42]. In vitro experiments with brain cell lines derived from M. myotis, however, suggest that high up-regulation of pattern recognition receptors and a vigorous interferon-mediated response are invoked to control lyssavirus infection [43]. Lyssaviruses identified to circulate in bat species of the studied geographic region include European bat lyssavirus type 1 and 2, Bokeloh bat lyssavirus, Lleida bat lyssavirus, West Caucasian bat lyssavirus, and lyssaviruses Aravan, Khujand and Irkut [7] [9] [13] [14] [16]. The commercial ELISA kit used in this study employes RABV glycoprotein to detect anti-rabies antibodies. It is a multi-species assay using protein A (a surface molecule of Staphylococcus aureus binding to the Fc portion of immunoglobulins from a large number of animal species) coupled with peroxidase as an enzymatic conjugate. Detection of cross-reactivity between sera against divergent members of the lyssavirus genus was recorded [1] [44] [45] [46]. Though it is for the classical RABV detection, using this kit has proved advantageous when testing samples from mixed bat populations collected over an extensive area with several possible non-RABV lyssaviruses circulating. However, bats exposed to phylogenetically distant lyssaviruses may not have been detected; meaning that rabies prevalence is probably underestimated in our study. ELISA including RABV as one of antigens has already been used to screen bats for rabies in Old World bats [47]. Interestingly, approximately 11% of 789 bats sampled in Northern Vietnam had neutralizing antibodies against RABV [48]. As ELISA detects binding antibodies to rabies glycoprotein G, it may differ in sensitivity and susceptibility when compared with standard virus neutralization assays, i.e. the fluorescent antibody virus neutralization (FAVN) or the rapid fluorescent focus inhibition (RFFIT) [47]. The body mass of Palearctic bat species limits the blood volume available for non-lethal collection. It is therefore very difficult to yield enough blood for screening a range of different lyssaviruses using FAVN or RFFIT reactions. Though there are no validated standards for bat sera [47] and the study by Wasniewski et al. (2014) [49] suggested decreasing the cut-off threshold (ELISA tests) to increase the test sensitivity for domestic carnivores and wildlife, our data show that ELISA could be used as a quick screening tool for bat exposure to lyssaviruses. However, further validation and effectiveness testing of the method for bat samples in combination with reference neutralization tests will be necessary. Problems associated with diagnostic sensitivity, specificity, sample volume requirements and detection of multiple lyssaviruses could be resolved by development of a bead-based multiplex immunoassay [50] for bat rabies screening. To our knowledge, this is the first active surveillance study to report M. hilgendorfi, M. petax and M. gracilis as rabies seropositive, though M. hilgendorfi has previously been associated with rabies as M. leucogaster, recorded as positive in the 2002 Irkut lyssavirus isolation [7], is now considered to be M. hilgendorfi in the Baikal Lake region. In addition, seropositivity of M. petax and M. gracilis, sibling species to M. daubentonii and M. brandtii, respectively, matches the documented association between Myotis bat species and European bat lyssavirus type 2 (EBLV-2) [10] [15]. The overall rabies seroprevalence in our Palearctic collection was 3.2% and when we remove seropositive bat colony (Altai Krai, Russia), seroprevalence in the Europe decreases (2.0%). Myotis species represented the majority (84.8%) of seropositive records in this study (Table 2), though one seropositive sampling site stood out from all others with a seroprevalence level of 73.1%. While the two species confirmed positive (M. blythii and M. petax) occupied different roosts in the Bat’s cave (Peshchera letuchikh myshey, Chineta reserve, Altai Krai, Russia), seropositivity indicated possible interspecies lyssavirus circulation. Large, tightly clustered maternity colonies of M. blythii (Fig. 1), allow for high contact intensity and may be the driver of lyssavirus circulation at this site. An analogy can be drawn between this site and serological status associated with circulation of European bat lyssavirus type 1 (EBLV-1) in maternity colonies of Eptesicus serotinus in France [24] and M. myotis in Spain [21] [39] Myotis petax and Myotis gracilis are examples of cryptic bat species paired with Palearctic bat rabies positive species, i.e. Myotis daubentonii and Myotis brandtii [51]. Two boreal species, European Myotis brandtii and Asian Myotis gracilis which are closely related, belong to the New World clade. These morphologically almost indistinguishable species were long considered as synonyms of Myotis mystacinus. On the contrary, morphologically very similar Palearctic species, Asian Myotis petax and European Myotis daubentonii, appear distantly related within Old World Myotis. They were long considered a single species with the Western Palearctic M. daubentonii, but multiple lines of evidence [52] [53] clearly show that they are not even sister-species [51]. Even though the evolutionary history of these pairs is different, we may expect that the potential for harbouring pathogens will be similar, including the cross-species positivity. Fig. 1 Bats roosting in the Peshchera letuchikh myshey (Chineta reserve, Altai Krai, Russia). This Myotis blythii maternity colony comprises approximately 1400 bats, including females and young-of-the-year. As shown, individual bats are tightly clustered during the nursery period, allowing for high contact intensity Detection of cross-reactive antibodies using ELISA provides no clue to the identity of circulating lyssaviruses, however, other than it belongs to phylogroup I [1]. Interestingly, in 1991, Aravan lyssavirus (ARAV) was isolated from a male M. blythii in Kyrgyzstan (Central Asia) [9], suggesting that bats roosting in this cave need to be sampled for identification of circulating viruses. The higher percentage of seropositive adult females (48.5%) in this study corresponds with the findings of Picard-Meyer et al. (2011) and Robardet et al. (2017), suggesting that females may act as a suitable indicator of lyssavirus circulation in bat populations. In comparison, the more solitary behaviour of male bats results in a lower probability of exposure to lyssaviruses. Interestingly, male vampire bats were shown to contribute disproportionately to rabies spatial spread [54]. Our finding of no difference in seropositive prevalence between the sexes in juvenile bats (ca. two months) is most likely due to an equal chance of exposure during the nursery period or the passing of maternally acquired antibodies to both sexes. Conclusions Our results suggest that active surveillance of bat rabies is more efficient (i.e. providing higher rate of positive findings) when samples are collected during the bat’s active season, rather than during the hibernation period. While both innate and adaptive immune functions are suppressed during hibernation [41] [42] [55], the rate of infection transmission contacts is also reduced. Colony size, species richness, high mobility and the social behaviour of bats during the active season when summer colonies are formed and swarming occurs provide higher chances for pathogen transmission between individuals [18] [21] [56]. Likewise, during the active season the bats´ immune system better responds to antigenic stimulation increasing antibody production and chances to detect seropositive animals in a population in which lyssaviruses circulate. Further, it would appear that human-bat contact in natural habitats poses a higher exposure risk, which is of importance as regards public health. Obviously, there is a considerable difference in risk of exposure for the general public and bat specialists, researchers and wildlife rehabilitators. The risks, however, cannot be generalised, meaning that geographic, species diversity, population density and annual life cycle issues have to be considered separately in each case. Methods Material collection Blood samples were taken from 1027 bats of 25 different species in Bulgaria (n = 77), the Czech Republic (n = 574), Poland (n = 110), Russia (n = 249) and Slovenia (n = 17) between 2014 and 2018. Every bat was sampled just one time. Using morphological traits and/or sequencing of the mitochondrial gene for cytochrome b (mtcyb), bat species were identified as Barbastella barbastellus (n = 11), Eptesicus nilssonii (n = 26), Eptesicus serotinus (n = 3), Miniopterus schreibersii (n = 44), Murina hilgendorfi (n = 42), Myotis alcathoe (n = 2), Myotis bechsteinii (n = 18), Myotis blythii (n = 20), Myotis bombinus (n = 2), Myotis brandtii (n = 33), Myotis dasycneme (n = 52), Myotis daubentonii (n = 86), Myotis emarginatus (n = 20), Myotis gracilis (n = 36), Myotis myotis (n = 328), Myotis mystacinus (n = 3), Myotis naterreri (n = 47), Myotis petax (n = 27), Nyctalus noctula (n = 118), Pipistrellus pipistrellus (n = 3), Plecotus auritus (n = 50), Plecotus ognevi (n = 4), Rhinolophus euryale (n = 13), Rhinolophus ferrumequinum (n = 35) and Vespertilio murinus (n = 6) [31] [34] The external genitalia of each bat was inspected for gender determination and age was estimated based on epiphyseal ossification of the thoracic limb fingers and teeth abrasion [57]. Bats were sampled at underground hibernacula, summer colonies, swarming sites and rescue centres [31] [32] [58], with bats captured during the annual life cycle active season and the hibernation period numbering 247 and 780 specimens, respectively. For further detail on sampling site, habitat type, number of specimens examined, species and gender, see Tables 1 and 2. Every attempt was made to minimise the impact of disturbance, handling stress and sampling procedure duration. Bats emerging from hibernacula, roosts or swarming sites were captured using mist-nets, with active bats sampled immediately following capture and released as soon as possible thereafter. For those bats in torpor, we allowed a re-warming period of about 60 minutes to allow for efficient blood collection. In all cases, the skin surface was first disinfected with alcohol, the uropatagial vessel was then lanced using a sterile needle and a 15 µl blood sample was taken using a pipette tip. A drop of tissue surgical glue (Surgibond, SMI AG, Belgium) was then applied to seal the lacerated skin and prevent further bleeding. Prior to release, all bats were provided with 5% glucose for rapid replenishment of energy and fluids. Detection of rabies virus anti-glycoprotein antibodies Blood samples were examined using the indirect enzyme-linked immunosorbent assay (ELISA) (PLATELIA™ RABIES II kit, BIO-RAD, France), which allows for qualitative detection and titration of IgG antibodies against the rabies virus glycoprotein G in animal serum. Detailed data concerning the antigen, reagents and procedures used in the Platelia kit can be found in [59]. In general, we followed the manufacturer’s recommendations and used the supplied positive controls and quantification standards; the only modification to the recommended protocol being mixing of the whole blood sample with sample diluent immediately after blood collection and subsequent centrifugation to remove formed elements. Samples were then kept frozen (-20 °C) until serological analysis. Prior to testing bat samples we performed in-house validation by comparison of results obtained using the commercial PLATELIA RABIES II kit and the reference FAVN test. Positive (post-vaccination) and negative paired whole blood and serum samples of domestic and wild carnivores (n = 30) were used for this purpose, providing 100% specificity for both sample types and insignificant differences in sensitivity (88 and 94%, respectively). A set of these carnivore samples was included later in each run with bat samples and, together with positive and negative controls supplied with the kit, used to construct a calibration curve for conversion of measured optical density values into titres as well as to derive the classification thresholds. The following threshold levels were applied for titre interpretation: (a) negative results representing undetectable seroconversion (< 0.123 Equivalent Units/ml), (b) positive results with insufficient seroconversion (0.123 to 0.5 EU/ml), and (c) positive results with sufficient seroconversion (> 0.5 EU/ml). Percentage bat seropositivity in different habitat types was compared by testing the difference between two proportions. The Chi-square test with Yate’s correction was used to detect seropositivity patterns by age and sex. Abbreviations EBLV-1European bat lyssavirus type 1 EBLV-2European bat lyssavirus type 2 ELISAEnzyme-linked immunosorbent assay RABVRabies virus Publisher’s Note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations. Acknowledgements We would like to thank Alena Hulcova, Anna Kompauerova, Natália Martínková, Margita Prikrylova, Alexander Osintsev, Alexander Stashkin and Julia Kosheleva for their help with research. We are grateful to Dr. Kevin Roche for correction and improvement of the English text. Authors' contributions VS, JZ and JP conceived and designed the study; JZ organised the field trips and collected material, with support from VS, JB, HBan, TB, HBer, ADB, TH, HD, TK, GA, PL, NA, OLO, VP, PP, MPT, APS, FT and JP; VS and JB performed the laboratory measurements; VS, JZ and JP analysed the data and drafted the manuscript, to which all authors contributed with critical comments. All authors have read and approved the manuscript. Funding This study was supported by the Czech Science Foundation (Grant No. 17-20286S) and Nacional Science Fund, Ministry of Education and Sciences, Bulgaria (Grant No. KП-06-H31/14). The funders had no role in the study design, data analysis, and decision to publish, or preparation of the manuscript. Availability of data and materials All data needed to evaluate the conclusions are present in the paper. The raw datasets used and/or analysed during the current study available from the corresponding author on reasonable request. Ethics approval and consent to participate Team members were authorised to handle wild bats according to the Czech Certificate of Competency (No. CZ01341; § 17, Act No. 246/1992 Coll.). Experimental procedures were approved by the Ethical Committee of the Czech Academy of Sciences (Document No. 169/2011). In the Czech Republic, capturing of bats for sample collection complied with Czech Law No. 114/1992 on Nature and Landscape Protection and was based on permits issued by the Agency for Nature Conservation and Landscape Protection of the Czech Republic (01662/MK/2012S/00775/MK/2012, 866/JS/2012 and 00356/KK/2008/AOPK) and Podyjí National Park (NPP 0312/2018). The II Local Ethical Commission in Wrocław approved sampling at the Natura 2000 site “Nietoperek” in Poland (No. 45/2015). Sampling in Bulgaria, Poland, Russia and Slovenia was approved by the Bulgarian Ministry of Environment and Water (No. 645/13.08.2015 and No. 683/04.07.2016), the Regional Directorate for Environmental Protection in Gorzów Wielkopolski, Poland (No. WPN-I-6205.10.2015.AI), the Institute of Plant and Animal Ecology - Ural Division of the Russian Academy of Sciences (No. 16353 − 2115/325) and the Ministry of Environment and Spatial Planning of the Republic of Slovenia, Slovenian Environment Agency (No. 35601-35/2010-6). Consent for publication Not applicable. Competing interests The authors declare that they have no competing interests. ==== Refs References 1. Malerczyk C Freuling C Gniel D Giesen A Selhorst T Cross-neutralization of antibodies induced by vaccination with purified chick embryo cell vaccine (PCECV) against different Lyssavirus species Hum Vaccin Immunother 2014 10 2799 804 10.4161/21645515.2014.972741 25483634 2. Hampson K Coudeville L Lembo T Sambo M Kieffer A Attlan M Estimating the global burden of endemic canine rabies PLoS Negl Trop Dis 2015 9 e0003709 10.1371/journal.pntd.0003709 25881058 3. Slate D Rupprecht CC Rooney JA Donovan D Lein DH Chipman RH Status of oral rabies vaccination in wild carnivores in the United States Virus Res 2005 111 68 76 10.1016/j.virusres.2005.03.012 15896404 4. Slate D Algeo TP Nelson KM Chipman RB Donovan D Oral Rabies Vaccination in North America: Opportunities, Complexities, and Challenges PLOS Negl Trop Dis 2009 3 12 e549 10.1371/journal.pntd.0000549 20027214 5. Wasniewski M Almeid I Baur A Bedekovic T Boncea D Chaves LB First international collaborative study to evaluate rabies antibody detection method for use in monitoring the effectiveness of oral vaccination programmes in fox and raccoon dog in Europe J Virol Methods 2016 238 77 85 10.1016/j.jviromet.2016.10.006 27751949 6. Robardet E Bosnjak D Englund L Demetriou P Martín PR Cliquet F Zero Endemic Cases of Wildlife Rabies (Classical Rabies Virus, RABV) in the European Union by 2020: An Achievable Goal Trop Med Infect Dis 2019 4 124 10.3390/tropicalmed4040124 6958318 7. Botvinkin AD Poleschuk EM Kuzmin IV Borisova TI Gazaryan SV Yager P Novel lyssaviruses isolated from bats in Russia Emerg Infect Dis 2003 9 1623 5 10.3201/eid0912.030374 14720408 8. Fooks AR Brookes S Johnson N McElhinney LM Hutson AM European bat lyssaviruses: an emerging zoonosis Epidemiol Infect. 2003 131 1029 1039 14959767 9. Kuzmin IV Orciari LA Arai YT Smith JS Hanlon CA Kameoka Y Bat lyssaviruses (Aravan and Khujand) from Central Asia: phylogenetic relationships according to N, P and G gene sequences Virus Res 2003 97 65 79 10.1016/S0168-1702(03)00217-X 14602198 10. Vos A Kaipf I Denzinger A Fooks AR Johnson N Müller T European bat lyssaviruses - an ecological enigma Acta Chiropterol 2007 9 283 96 10.3161/1733-5329(2007)9[283:EBLAEE]2.0.CO;2 11. Johnson N Selden D Parsons G Healy D Brookes SM McElhinney LM Hutson AM Fooks AR Isolation of a European Bat Lyssavirus type 2 from a Daubentons bat in the United Kingdom Vet Rec 2003 152 383 38 10.1136/vr.152.13.383 12696703 12. Johnson N Vos A Freuling C Tordo N Fooks AR Müller T Human rabies due to lyssavirus infection of bat origin Vet Microbiol 2010 142 151 9 10.1016/j.vetmic.2010.02.001 20188498 13. Ceballos NA Morón SV Berciano JM Nicolás O López CA Juste J Novel lyssavirus in bat, Spain Emerg Infect Dis 2013 19 793 5 10.3201/eid1905.121071 23648051 14. Freuling CM Beer M Conraths FJ Finke S Hoffmann B Keller B Novel lyssavirus in Natterer’s bat, Germany Emerg Infect Dis 2011 17 1519 22 21801640 15. Schatz J Fooks AR Mcelhinney L Horton D Echevarria J Vázquez-Moron S Bat Rabies Surveillance in Europe Zoonoses Public HLTH 2013 60 22 34 10.1111/zph.12002 16. Banyard AC Evans JS Luo TR Fooks AR Lyssaviruses and bats: emergence and zoonotic threat Viruses 2014 6 8 2974 90 10.3390/v6082974 25093425 17. Blackwood JC Streicker DG Altizer S Rohani P Resolving the roles of immunity, pathogenesis, and immigration for rabies persistence in vampire bats PNAS 2013 110 20837 10.1073/pnas.1308817110 24297874 18. Colombi D Serra-Cobo J Métras R Apolloni A Poletto Ch Mechanisms for lyssavirus persistence in non-synanthropic bats in Europe: insights from a modeling study Sci Rep 2019 9 537 10.1038/s41598-018-36485-y 30679459 19. Van Der Poel WHM Van Der Heide R Verstraten E Takumi K Lina PHC Kramps JA European bat lyssaviruses, the Netherlands Emerg Infect Dis 2005 11 1854 9 10.3201/eid1112.041200 16485470 20. Hajkova P Pikula J Veterinary treatment of evening bats (Vespertilionidae) in the Czech Republic Vet Rec 2007 161 139 40 10.1136/vr.161.4.139 17660469 21. Amengual B Bourhy H López-Roig M Serra-Cobo J Temporal dynamics of European bat lyssavirus type 1 and survival of Myotis myotis bats in natural colonies PLoS ONE 2007 2 e566 10.1371/journal.pone.0000566 17593965 22. Picard-Meyer E Dubourg-Savage MJ Arthur L Barataud M Bécu D Bracco S Active surveillance of bat rabies in France: A 5-year study (2004–2009) Vet Microbiol 2011 51 390 5 10.1016/j.vetmic.2011.03.034 23. Nokireki T Huovilainen A Lilley T Kyheröinen EM Ek-Kommonen C Sihvonen L Bat rabies surveillance in Finland BMC Vet Res 2013 9 174 10.1186/1746-6148-9-174 24011337 24. Robardet E Borel C Moinet M Jouan D Wasniewski M Barrat J Longitudinal survey of two serotine bat (Eptesicus serotinus ) maternity colonies exposed to EBLV-1 (European Bat Lyssavirus type 1): Assessment of survival and serological status variations using capture-recapture models PLoS Negl Trop Dis 2017 11 e0006048 10.1371/journal.pntd.0006048 29149215 25. Leopardi S, Priori P, Zecchin B, Poglayen G, Trevisiol K, Lelli D, et al. Active and passive surveillance for bat lyssaviruses in Italy revealed serological evidence for their circulation in three bat species. Epidemiol Infect. 2018;1–6. 26. Šimić I Lojkić I Krešić N Cliquet F Picard-Meyer E Wasniewski M Molecular and serological survey of lyssaviruses in Croatian bat populations BMC Vet Res 2018 14 1 274 10.1186/s12917-018-1592-z 30189884 27. Pikula J Bandouchova H Novotny L Meteyer CU Zukal J Irwin N Histopathology confirms white-nose syndrome in bats in Europe J Wildl Dis 2012 48 207 11 10.7589/0090-3558-48.1.207 22247393 28. Zukal J Bandouchova H Bartonička T Berkova H Brack V Brichta J White-nose syndrome fungus: A generalist pathogen of hibernating bats PLoS ONE 2014 9 e97224 10.1371/journal.pone.0097224 24820101 29. Bandouchova H Bartonicka T Berkova H Brichta J Cerny J Kovacova V Pseudogymnoascus destructans : Evidence of virulent skin invasion for bats under natural conditions, Europe Transbound Emerg Dis 2015 62 1 5 10.1111/tbed.12282 25268034 30. Flieger M Bandouchova H Cerny J Chudickova M Kolarik M Kovacova V Vitamin B2 as a virulence factor in Pseudogymnoascus destructans skin infection Sci Rep 2016 6 33200 10.1038/srep33200 27620349 31. Zukal J Bandouchova H Brichta J Cmokova A Jaroň KS Kolařík M White-nose syndrome without borders: Pseudogymnoascus destructans infection tolerated in Europe and Palearctic Asia but not in North America Sci Rep 2016 6 19829 10.1038/srep19829 26821755 32. Pikula J Bandouchova H Kovacova V Linhart P Piacek V Zukal J Reproduction of rescued Vespertilionid bats (Nyctalus noctula ) in captivity: Veterinary and physiologic aspects. Veterinary Clinics of North America Exot Anim Pract 2017 20 665 77 10.1016/j.cvex.2016.11.013 33. Bandouchova H Bartonička T Berkova H Brichta J Kokurewicz T Kovacova V Alterations in the health of hibernating bats under pathogen pressure Sci Rep 2018 8 6067 10.1038/s41598-018-24461-5 29666436 34. Kovacova V Zukal J Bandouchova H Botvinkin AD Harazim M Martínková N White-nose syndrome detected in bats over an extensive area of Russia BMC Vet Res 2018 14 192 10.1186/s12917-018-1521-1 29914485 35. Martínková N Pikula J Zukal J Kovacova V Bandouchova H Bartonička T Hibernation temperature-dependent Pseudogymnoascus destructans infection intensity in Palearctic bats Virulence 2018 9 1 1734 50 10.1080/21505594.2018.1548685 36. Zahradníková A Kovacova V Martínková N Orlova MV Orlov OL Historic and geographic surveillance of Pseudogymnoascus destructans possible from collections of bat parasites Transbound Emerg Dis 2018 65 303 8 10.1111/tbed.12773 29181887 37. Davy CM Donaldson ME Bandouchova H Breit AM Dorville N Dzal YA Transcriptional host-pathogen responses of Pseudogymnoascus destructans and three species of bats with white-nose syndrome Virulence 2020 11 1 781 94 10.1080/21505594.2020.1768018 32552222 38. Schatz J Teifke JP Mettenleiter TC Aue A Stiefel D Müller T Lyssavirus distribution in naturally infected bats from Germany Vet Microbiol 2014 169 33 41 10.1016/j.vetmic.2013.12.004 24440375 39. Serra-Cobo J Amengual B Abellán C Bourhy H European bat lyssavirus infection in Spanish bat populations Emerg Infect Dis 2002 8 413 20 10.3201/eid0804.010263 11971777 40. Baker M Schountz LT Wang LF Antiviral immune responses of bats: A review Zoonoses Public HLTH 2013 60 104 16 10.1111/j.1863-2378.2012.01528.x 41. Pikula J Heger T Bandouchova H Kovacova V Nemcova M Phagocyte activity reflects mammalian homeo- and hetero-thermic physiological states BMC Vet Res 2020 16 1 1 13 10.1186/s12917-020-02450-z 31900161 42. Heger T Zukal J Seidlová V Němcová M Nečas D Measurement of phagocyte activity in heterotherms Acta Veterinaria Brno 2020 89 1 79 87 10.2754/avb202089010079 43. He X Korytář T Zhu Y Pikula J Bandouchova H Zukal J Establishment of Myotis myotis cell lines - model for investigation of host-pathogen interaction in a natural host for emerging viruses PLoS ONE 2014 9 e109795 10.1371/journal.pone.0109795 25295526 44. Brookes SM Parsons G Johnson N McElhinney LM Fooks AR Rabies human diploid cell vaccine elicits cross-neutralising and cross-protecting immune responses against European and Australian bat lyssavirusses Vaccine 2005 23 32 4101 9 10.1016/j.vaccine.2005.03.037 15964478 45. Brookes SM Healy DM Fooks AR Ability of rabies vaccine strains to elicit crossneutralising antibodies Dev Biol (Basel) 2006 125 185 93 16878476 46. Malerczyk C Selhorst T Tordo N Moore S Müller T Antibodies induced by vaccination with purified chick embryo cell culture vaccine (PCECV) cross-neutralize non-classical bat lyssavirus strains Vaccine 2009 27 5320 5 10.1016/j.vaccine.2009.06.095 19615958 47. Reynes J-M Molia S Audry L Hout S Ngin S Walston J Serologic evidence of lyssavirus infection in bats, Cambodia Emerg Infect Dis 2004 10 2231 4 10.3201/eid1012.040459 15663870 48. Nguyen AT Nguyen TT Noguchi A Bat lyssaviruses, northern Vietnam Emerg Infect Dis 2014 20 1 161 3 10.3201/eid2001.130813 24377728 49. Wasniewski M Labbe A Tribout L Rieder J Labadie A Schereffer JL Evaluation of a rabies ELISA as an alternative method to seroneutralisation tests in the context of international trade of domestic carnivores J Virol Methods 2014 195 211 20 10.1016/j.jviromet.2013.10.021 24161815 50. Smits GP van Gageldonk PG Schouls LM van der Klis FR Berbers GAM Development of a Bead-Based Multiplex Immunoassay for Simultaneous Quantitative Detection of IgG Serum Antibodies against Measles, Mumps, Rubella, and Varicella-Zoster Virus Clin Vaccine Immunol 2012 19 396 10.1128/CVI.05537-11 22237896 51. Ruedi M Stadelmann B Gager Y Douzery EJP Francis CM Lin LK Guillen-Servent A Cibois A Molecular phylogenetic reconstructions identify East Asia as the cradle for the evolution of the cosmopolitan genus Myotis (Mammalia, Chiroptera) Mol Phylogenet Evol 2013 69 437 49 10.1016/j.ympev.2013.08.011 23988307 52. Datzmann T Dolch D Batsaikhan N Kiefer A Helbig-Bonitz M Cryptic diversity in Mongolian vespertilionid bats (Vespertilionidae, Chiroptera, Mammalia). Results of the Mongolian–German biological expeditions since 1962, No. 299 Acta Chiropterol 2012 14 243 64 10.3161/150811012X661594 53. Matveev VA Kruskop SV Kramerov DA Revalidation of Myotis petax Hollister, 1912 and its new status in connection with M. daubentonii (Kuhl, 1817) (Vespertilionidae, Chiroptera) Acta Chiropterol 2005 7 23 37 10.3161/1733-5329(2005)7[23:ROMPHA]2.0.CO;2 54. Streicker DG Winternitz JC Satterfield DA Condori-Condori RE Broos A Tello C Host-pathogen evolutionary signatures reveal dynamics and future invasions of vampire bat rabies PNAS 2016 113 39 10926 10.1073/pnas.1606587113 27621441 55. Bouma HR Carey HV Kroese FGM Hibernation: the immune system at rest? J Leukoc Biol 2010 88 619 24 10.1189/jlb.0310174 20519639 56. Serra-Cobo J López-Roig M Seguí M Sánchez LP Nadal J Ecological Factors Associated with European Bat Lyssavirus Seroprevalence in Spanish Bats PLoS ONE. 2013 8 e64467 10.1371/journal.pone.0064467 23700480 57. Brunet-Rossinni A, Kwilkinson GS. Methods for age estimation and the study of senescence in bats.2009; In: Ecological and Behavioral Methods for the Study of Bats (eds Kunz, T. H. and Parsons, S.) 315–325 (The Johns Hopkins University Press, 2009). 58. Pikula J Amelon SK Bandouchova H Bartonicka T Berkova H Brichta J White-nose syndrome pathology grading in Nearctic and Palearctic bats PLoS ONE 2017 12 e0180435 10.1371/journal.pone.0180435 28767673 59. Servat A Feyssaguet M Blanchard I Morize JL Schereffer JL Boue F A quantitative indirect ELISA to monitor the effectiveness of rabies vaccination in domestic and wild carnivores J Immunol Methods 2007 318 1 10 10.1016/j.jim.2006.07.026 17166510