1. Introduction
During the last two decades, three major outbreaks of coronavirus-caused diseases occurred in the world. The first epidemic (caused by the Severe Acute Respiratory Syndrome Coronavirus, SARS-CoV) happened in China in 2002 and caused 744 deaths among 8 096 patients [
1,
2]. After the outbreak of MERS-CoV disease in the Middle East in 2012 a total of 2 298 cases were confirmed including 811 deaths [
3]. When several cases of a strange pneumonia were reported in city of Wuhan (Hubei, China) in late 2019 [
4], barely anyone could imagine that this would lead to a pandemic causing 6,990,067 deaths among 773,119,173 patients (as of December 24
th 2023) [
5].
Coronaviruses (among which is the SARS-CoV-2, the betacoronavirus responsible for the latest pandemic) can be hosted by multiple mammal species like bats, rats, mice, dogs, camels, calves, turkeys, rabbits, pigs. Although they mainly cause respiratory disorders, some particular strains of coronaviruses may result in gastroenteritis, hepatitis, encephalitis or peritonitis in infected individual [
6,
7]. In humans, coronavirus infections may vary from asymptomatic to severe pneumonia accompanied by fever, cough, or gastrointestinal irritation [
8,
9]
Among animal hosts, bats are considered to be of major epidemiological importance. This is because they have been reported to be a reservoir and carriers for many viruses [
10,
11,
12]. Moreover, they form large colonies making it easier for the viruses to transfer between individuals. Finally, bats can travel vast distances spreading the diseases both within their own species and across other animal species [
13,
14]. Intermediate coronavirus hosts are crucial for the animal-human transmission as people rarely get infected from the bats themselves [
15]. Masked palm civets and pangolins are believed to be the direct source of SARS-CoV and SARS-CoV-2 infections for humans although they don’t play the role of animal reservoirs. It contrasts with dromedary camels which are both intermediate hosts and reservoirs for MERS-CoV [
16,
17]. Searching for animal hosts and reservoirs, developing the knowledge about relations between them and possible consequences of their coexistence is vital for better understanding the dynamics of the diseases caused by coronaviruses. In this research we focused on virological screening of bats residing in military bunkers in the Natura 2000 site PLH080003 "Nietoperek", western Poland.
2. Materials and Methods
2.1. Bat Sampling
The samples were collected in section 7 of the central sector of Międzyrzecz Fortified Front (MFF) (Festungsfront Oder-Warthe-Bogen) (52°25′N, 15°32′ E) situated in Western Poland (Map. 1). The fortification system was built in the period 1934–1944 and is composed of a large main underground corridor system (ca 32 km long, 20-30 m below ground) and one or two storey bunkers situated on the surface. Every year up to 40,000 bats of 13 species hibernate in both the main underground system of "Nietoperek" and the surrounding bunkers, making it one of the largest and most important bat hibernation sites in Europe [
18,
19,
20]. Tracking data from the Dresden bat ringing center (Fledermausmarkierungszentrale, Germany) shows bats migrating to "Nietoperek" from the expansive Central European Lowlands. The furthest travel distances recorded to "Nietoperek" are 257 km for
Myotis daubentonii, 226.7 km for
M. myotis, and 242.1 km for
M. brandtii. The minimal convex polygon (MCP) for large mouse-eared bats, based on recaptured individuals, spans at least 17,000 km² across significant parts of German territories: Brandenburg, Mecklemburg-Vorpommern, Sachsen-Anhalt, and western Poland (Lubuskie, Wielkopolskie, Zachodniopomorskie Voivodeships) [
21]. Current migration patterns in Nietoperek have been deduced solely from ringing data; future research using isotopes and DNA analysis may uncover a broader influence of this location. In November 2007, the underground system with the surrounding surface area of 7377.37 ha became protected as Natura 2000 site “Nietoperek” (area code PLH080003), with the bats as a target of protection. Section 7 is visited by tourists year-round, with potential human infection risk. Bats were individually collected from their colonies. Their species and sex were identified visually. We measured each bat’s body weight and arm length. For each bat, oral and anal swabs were taken using laboratory swabs. These swabs were then placed in Eppendorf tubes filled with a virus-inactivating medium and transported to the lab under cool conditions at +4°C for RNA extraction. After sample collection and measurements, the bats were carefully returned to their colony to minimize disturbance.
Map 1.
Lubrza “Nietoperek Reserve” is located in western Poland (52°23′25″N 15°31′02″E). Map Data, Google, 2023.
Map 1.
Lubrza “Nietoperek Reserve” is located in western Poland (52°23′25″N 15°31′02″E). Map Data, Google, 2023.
2.2. RNA Extraction
The RNA was extracted from oral and anal swab samples and used in PCR assays to detect coronaviruses. The extraction was done using the QIAamp 96 Virus QIAcube HT Kit and QIAcube HT system (QIAGEN, Hilden, Germany) following the manufacturer’s instructions. Briefly, the swabs were transferred to tubes containing 800 µL of ATL buffer and 30 µL of proteinase K. The tubes were incubated for 1 hour at 56 °C and homogenized by vortexing within intervals of 15 minutes. After this, 200 µL of lysate was transferred to a sample block and the extraction procedure was carried out with QIAcube, using onboard lysis. The final elution volume was set to 120 µL, and the extracted RNA was stored in -80 °C.
2.3. Pan-Coronavirus Molecular Screening Assay
Initially the RNA extract obtained from oral and anal swab were screened for the presence of coronaviruses using a broad range one-step real time RT-PCR assay. The reaction mix was prepared with qScript One-Step SYBR® Green RT-qPCR Kit (QUANTABIO, Beverly, MA, USA) and 600 nM of each primer, 11-FW (5’-TGATGATGSNGTTGTNTGYTAYAA-3’) and 13-RV (5’-GCATWGTRTGYTGNGARCARAATTC-3’) [
22]. The reactions included 5µL of template RNA and prepared to a final volume of 20 µL. The thermocycling conditions were as follows: first the reactions were incubated at 50 °C for 10 min, followed by initial denaturation at 95 °C for 5 min and then 45 cycles consisting of a denaturation step at 95 °C for 10 sec, an annealing step at 50 °C for 20 sec, and an elongation step at 72 °C for 30 sec. Fluorescence was collected after each cycle. At the end, a melt curve was obtained by increasing the temperature by 0.5°C /5 sec, from 55 to 95 °C. The amplification was carried out using the AriaMx Real-time PCR System and the fluorescence was analysed using AriaMx Software (Agilent Technologies, Santa Clara, CA, USA). Samples selected based on their temperature melting curve were taken for further investigation. RNA extracted from SARS CoV2 was used as positive control in the screening and conventional PCR.
2.4. Pan-Coronavirus Confirmatory PCR Assay
A broad range, semi-nested PCR assay was used to confirm the presence of coronaviruses in the previously selected samples. Complementary DNA (cDNA) was produced from the extracted RNA using Maxima Reverse Transcriptase (Thermo Scientific, Waltham, Massachusetts, USA) and used in the confirmatory PCR assay. The RT reactions were prepared to a final volume of 20 µL. Initially, a reaction mix containing 5 µL of template RNA, 1 µL of random hexamers (50 µM), 1 dNTP’s (10 mM each) and 7.5 µL of water was prepared and incubated at 65 °C for 5 min. Next, we added a second reaction mix containing RT Buffer to final concentration of 1X, 20 U of RiboLock RNase Inhibitor (Thermo Scientific™) and 200 U of Maxima Reverse Transcriptase. The final mix was initially incubated at 25 °C for 10 min, followed by 30 min at 50 °C and 5 min at 85°C, after which it was stored at -20 °C to be used as template in the confirmatory PCR reaction.
To confirm the presence of coronaviruses in the selected samples we used the primers targeting the RdRp published by Holbrook and collaborators [
23]. In this study we used the DreamTaq DNA Polymerase kit (Thermo Scientific) and the reactions were modified as follows. Initial and semi-nested reactions were prepared with 1X buffer and 0.2 mM of dNTP’s (each). In the first reaction, 2.5 µL of template cDNA was mixed with 2.5 U of polymerase, 4 % of DMSO and 1 µM of each primer – F1 (5’-GGTGGGAYTAYCCHAARTGYGA-3’), R1 (5’-CCRTCATCAGAHARWATCAT-3’) and R2 (5’-CCRTCATCACTHARWATCAT-3’), to a 25 µL final volume reaction. The second reaction was prepared with 2 µL of template (PCR product of the first reaction), 1.25 U of polymerase and 400 nM of each primer - F2 (5’-GAYTAYCCHAARTGTGAYAGA-3’), F3 (5’-GAYTAYCCHAARTGTGAYMGH-3’), R1, R2), to a final volume of 50 µL. The thermocycling conditions for both reactions were: an initial denaturation step at 94 °C for 3 min; 40 cycles of denaturation at 94 °C for 30 sec, annealing at 48 °C for 30 sec, and elongation at 72 °C for 40 sec; and a final elongation step at 72 °C for 5 min. The product of both reactions was visualized in 2% agarose gel stained with ethidium bromide.
2.5. SARS-CoV-2 Specific PCR Assay
In parallel to the confirmatory PCR, the presence of SARS-CoV-2 was investigated with the Luna® SARS-CoV-2 RT-qPCR Multiplex Assay Kit (New England Biolabs, Ipswich, Massachusetts, USA), a SARS-CoV-2 specific qPCR assay, as per manufacturer’s instructions. This assay was carried out on samples selected based on the initial screening assay.
2.6. Sequencing
The PCR products of reactions that amplified fragments of expected length were purified using the GeneJET PCR Purification Kit (Thermo Scientific) and sent for sequencing. The sequencing was carried out using Applied Biosystems BigDye Terminator v3.1 Cycle Sequencing Kit (Part No. 4336921) according to the manufacturer’s instructions and further analysed using Applied Biosystems ABI3130XL Genetic Analyzer (16-capillaries).
3. Results
We sampled 138 individuals belonging to six bat species: greater mouse-eared bat
Myotis myotis, Natterer’s bat
M. nattereri, western barbastelle
Barbastella barbastellus, Daubenton’s bat
M. daubentonii, Brandt’s bat
M. brandtii and whiskered bat
M. mystacinus (
Table 1).
As a first screen, samples were sorted for the likelihood of having coronaviruses using a pan coronavirus qPCR. The criteria for selecting samples at this initial screening stage were adjusted to include samples with very low viral load and virus genetic diversity. Therefore, sample with both a threshold cutting value (Ct) lower than those of the no-template controls (NTC) and negative controls, and a melting temperature close to that of the positive control, were used for follow-up assays. A total of 83 individuals out of 138 had one or both (oral and anal) samples selected for a pan Coronavirus PCR confirmatory assay targeting the RdRp. In this assay, a single positive sample was detected, an anal swab sample from catchment 5 (A5). The amplicon purified and prepared for sanger sequencing. Unfortunately, no readable sequence could be obtained.
Finally, the presence of SARS-CoV-2 was investigated using a commercial diagnostic test. No SARS-CoV-2 positive samples were found.
4. Discussion
The prevalence of coronaviruses in bats exhibits significant spatial and temporal variation, however the mechanistic factors aren’t completely clear [
24]. Factors like animal health, colony characteristics, and reproductive cycles can directly or indirectly influence prevalence. For instance, a study carried out in China, 2006, revealed an overall prevalence of 6.5%, with some colonies reaching rates as high as 55% [
25]. In Northern Germany, the overall prevalence was 9.8%, varying between 5.2% and 25.4% among different bat species [
26]. Additionally, a longitudinal study in Zimbabwe revealed varying prevalence of coronavirus in bats at the Chirundu farm site, ranging from 1.95% to 44.2% during pregnancy and weaning periods, respectively [
27]. In Magweto cave, during the pregnancy period, prevalence reached 10.5%, while for
Macronycteris gigas, residing in the same cave, it was recorded at 56.4%. The examples above illustrate significant variability in coronavirus prevalence in different regions, host species and under different periods. Even different species in the same area and under similar conditions might exhibit very different coronavirus prevalence.
In our study, out of 138 individuals, only one anal swab tested positive in conventional PCR, and no readable sequence was obtained. The number contrasts with findings in northern Germany, relatively closer to our study area, where bats caught outside the shelter during their activity period showed higher coronavirus prevalence (3). The low prevalence in the colony studied here might be due to unfavorable physiological and behavioral conditions during bats’ hibernation period, which might affect viral replication and reduce bats overall mobility and interactions, thereby decreasing the probability of virus transmission. Contrary to our results, Subudhi et al. (2017) found coronavirus persistence for up to 4 months in North American little brown bats during hibernation in laboratory settings [
28]. Nonetheless, this data is limited to a specific host and to laboratory conditions and can be different to other host species in the natural environment, as our studied specimens. The hibernation factor may explain our unexpected results, but a comprehensive understanding of the role of hibernation in coronavirus maintenance and transmission requires further research. At this stage, we believe the risk of coronavirus transmission from bats to tourists visiting underground sites is negligible.
5. Conclusions
We did not find hibernating bats to be infected with SARS-CoV-2 virus. And the very low overall prevalence of coronavirus (1/138) might be caused by the hibernation process. Hibernating bats are not the most optimal source for screening of zoonotic pathogens. On the one side they are easy to access, and material collection can be done with a relatively small disturbance to animals. On the other hand, hibernation cause significant decrease of metabolic processes in bat organisms. Therefore, hibernating bats produce very little amount of guano, and the bat metabolic changes might affect virus-host interactions. Considering past reports on SARS-CoV-2 origins, and susceptibility of bats to carry highly zoonotic viruses we believe that bat populations should be constantly screened against emerging and reemerging diseases to provide most accurate epidemiological image.
Author Contributions
The study was conceived and designed by MG & TK. Supervision of the biomonitoring by MG and TS. Sampling permits: TK. Samples were collected in the field by GA, TK, AR, AG and MG. Bats handling and measurements: GA, TK, AR. Molecular analysis and laboratory work was conducted by AG, JN, MG, VB, LD, ES, VV, RK. Data handling: LD, VB, MG. The manuscript was written by AG, LD, VB and MG in consultation with all co-authors. MG, LD and VB revised the manuscript. All authors accepted the final manuscript version.
Funding
This research was co-funded through the 2018–2019 BiodivERsA joint call for research proposals under the BiodivERsA3 ERA-Net COFUND program; the funding organisations ANR (France), DFG (Germany), EPA (Ireland), FWO (Belgium), and NCN (Poland). JN, MG and AG were supported by the National Science Centre, Poland, under the BiodivERsA3 programme (2019/31/Z/NZ8/04028).
Institutional Review Board Statement
This study was carried out in strict accordance with the recommendations in the Guidelines for the Care and Use of Laboratory Animals of the Polish National Ethics Committee for Animal Experimentation and according to the Polish national law for field work involving the trapping and culling of wild unprotected vertebrates for scientific purposes (Resolution No. 12/2022 of Polish National Ethics Committee for Animal Experiments, 11th March 2022). The study was performed according to the ARRIVE guidelines 2.0. Data collection was carried out under permission no. WPN-I.6205.24.2021.MG issued by Regional Director for Environmental Protection in Gorzów Wielkopolski (Poland)
Data Availability Statement
All data used in this study are provided within the manuscript.
Acknowledgments
MG thaks Ewa Zieliniewicz for her assistance in laboratory work.
Conflicts of Interest
The authors declare no conflict of interest.
References
- World Health Organisation Summary of Probable SARS Cases with Onset of Illness from 1 November 2002 to 31 July 2003; 2015. 1 November 2003.
- Xu, R.-H.; He, J.-F.; Evans, M.R.; Peng, G.-W.; Field, H.E.; Yu, D.-W.; Lee, C.-K.; Luo, H.-M.; Lin, W.-S.; Lin, P.; et al. Epidemiologic Clues to SARS Origin in China. Emerg Infect Dis 2004, 10, 1030–1037. [Google Scholar] [CrossRef]
- Al Awaidy, S.T.; Khamis, F. Middle East Respiratory Syndrome Coronavirus (MERS-CoV) in Oman: Current Situation and Going Forward. Oman Med J 2019, 34, 181–183. [Google Scholar] [CrossRef] [PubMed]
- World Health Organisation WHO-Convened Global Study of Origins of SARS-CoV-2: China Part; 2021.
- World Health Organisation WHO COVID-19 Dashboard. Available online: https://covid19.who.int/ (accessed on 4 January 2024).
- Sharma, A.; Ahmad Farouk, I.; Lal, S.K. COVID-19: A Review on the Novel Coronavirus Disease Evolution, Transmission, Detection, Control and Prevention. Viruses 2021, 13, 202. [Google Scholar] [CrossRef] [PubMed]
- Kahn, J.S.; McIntosh, K. History and Recent Advances in Coronavirus Discovery. Pediatric Infectious Disease Journal 2005, 24, S223–S227. [Google Scholar] [CrossRef] [PubMed]
- Wilder-Smith, A.; Teleman, M.D.; Heng, B.H.; Earnest, A.; Ling, A.E.; Leo, Y.S. Asymptomatic SARS Coronavirus Infection among Healthcare Workers, Singapore. Emerg Infect Dis 2005, 11, 1142–1145. [Google Scholar] [CrossRef]
- Chen, N.; Zhou, M.; Dong, X.; Qu, J.; Gong, F.; Han, Y.; Qiu, Y.; Wang, J.; Liu, Y.; Wei, Y.; et al. Epidemiological and Clinical Characteristics of 99 Cases of 2019 Novel Coronavirus Pneumonia in Wuhan, China: A Descriptive Study. The Lancet 2020, 395, 507–513. [Google Scholar] [CrossRef] [PubMed]
- Smith, I.; Wang, L.-F. Bats and Their Virome: An Important Source of Emerging Viruses Capable of Infecting Humans. Curr Opin Virol 2013, 3, 84–91. [Google Scholar] [CrossRef]
- Calisher, C.H.; Childs, J.E.; Field, H.E.; Holmes, K. V.; Schountz, T. Bats: Important Reservoir Hosts of Emerging Viruses. Clin Microbiol Rev 2006, 19, 531–545. [Google Scholar] [CrossRef]
- Hayman, D.T.S. Bats as Viral Reservoirs. Annu Rev Virol 2016, 3, 77–99. [Google Scholar] [CrossRef]
- Sharma, A.; Ahmad Farouk, I.; Lal, S.K. COVID-19: A Review on the Novel Coronavirus Disease Evolution, Transmission, Detection, Control and Prevention. Viruses 2021, 13, 202. [Google Scholar] [CrossRef]
- Ge, X.-Y.; Wang, N.; Zhang, W.; Hu, B.; Li, B.; Zhang, Y.-Z.; Zhou, J.-H.; Luo, C.-M.; Yang, X.-L.; Wu, L.-J.; et al. Coexistence of Multiple Coronaviruses in Several Bat Colonies in an Abandoned Mineshaft. Virol Sin 2016, 31, 31–40. [Google Scholar] [CrossRef] [PubMed]
- Xiao, K.; Zhai, J.; Feng, Y.; Zhou, N.; Zhang, X.; Zou, J.-J.; Li, N.; Guo, Y.; Li, X.; Shen, X.; et al. Author Correction: Isolation of SARS-CoV-2-Related Coronavirus from Malayan Pangolins. Nature 2021, 600, E8–E10. [Google Scholar] [CrossRef] [PubMed]
- de Wit, E.; van Doremalen, N.; Falzarano, D.; Munster, V.J. SARS and MERS: Recent Insights into Emerging Coronaviruses. Nat Rev Microbiol 2016, 14, 523–534. [Google Scholar] [CrossRef] [PubMed]
- Kadam, S.B.; Sukhramani, G.S.; Bishnoi, P.; Pable, A.A.; Barvkar, V.T. SARS-CoV-2, the Pandemic Coronavirus: Molecular and Structural Insights. J Basic Microbiol 2021, 61, 180–202. [Google Scholar] [CrossRef] [PubMed]
- Kokurewicz, T.; Apoznanski, G.; Gyselings, R.; Kirkpatrick, L.; De Bruyn, L.; Haddow, J.; Glover, A. 45 Years of Bat Study and Conservation in Nietoperek Bat Reserve (Western Poland). Nyctalus 2019, 19, 252–269. [Google Scholar]
- De Bruyn, L.; Gyselings, R.; Kirkpatrick, L.; Rachwald, A.; Apoznański, G.; Kokurewicz, T. Temperature Driven Hibernation Site Use in the Western Barbastelle Barbastella Barbastellus (Schreber, 1774). Sci Rep 2021, 11, 1464. [Google Scholar] [CrossRef] [PubMed]
- Apoznański, G.; Carr, A.; Gelang, M.; Kokurewicz, T.; Rachwald, A. Trophic Relationship between Salix Flowers, Orthosia Moths and the Western Barbastelle. Sci Rep 2023, 13, 7364. [Google Scholar] [CrossRef]
- Rogowska, K.; Kokurewicz, T. The Longest Migrations of Three Bat Species to the „Nietoperek” Bat Reserve (Western Poland). Berichte der Naturforschenden Gesellschaft der Oberlausitz 2007, 15, 53–60. [Google Scholar]
- Escutenaire, S.; Mohamed, N.; Isaksson, M.; Thorén, P.; Klingeborn, B.; Belák, S.; Berg, M.; Blomberg, J. SYBR Green Real-Time Reverse Transcription-Polymerase Chain Reaction Assay for the Generic Detection of Coronaviruses. Arch Virol 2007, 152, 41–58. [Google Scholar] [CrossRef]
- Holbrook, M.G.; Anthony, S.J.; Navarrete-Macias, I.; Bestebroer, T.; Munster, V.J.; van Doremalen, N. Updated and Validated Pan-Coronavirus PCR Assay to Detect All Coronavirus Genera. Viruses 2021, 13, 599. [Google Scholar] [CrossRef]
- Ruiz-Aravena, M.; McKee, C.; Gamble, A.; Lunn, T.; Morris, A.; Snedden, C.E.; Yinda, C.K.; Port, J.R.; Buchholz, D.W.; Yeo, Y.Y.; et al. Ecology, Evolution and Spillover of Coronaviruses from Bats. Nat Rev Microbiol 2022, 20, 299–314. [Google Scholar] [CrossRef]
- Tang, X.C.; Zhang, J.X.; Zhang, S.Y.; Wang, P.; Fan, X.H.; Li, L.F.; Li, G.; Dong, B.Q.; Liu, W.; Cheung, C.L.; et al. Prevalence and Genetic Diversity of Coronaviruses in Bats from China. J Virol 2006, 80, 7481–7490. [Google Scholar] [CrossRef]
- Gloza-Rausch, F.; Ipsen, A.; Seebens, A.; Göttsche, M.; Panning, M.; Drexler, J.F.; Petersen, N.; Annan, A.; Grywna, K.; Müller, M.; et al. Detection and Prevalence Patterns of Group I Coronaviruses in Bats, Northern Germany. Emerg Infect Dis 2008, 14, 626–631. [Google Scholar] [CrossRef]
- Chidoti, V.; De Nys, H.; Pinarello, V.; Mashura, G.; Missé, D.; Guerrini, L.; Pfukenyi, D.; Cappelle, J.; Chiweshe, N.; Ayouba, A.; et al. Longitudinal Survey of Coronavirus Circulation and Diversity in Insectivorous Bat Colonies in Zimbabwe. Viruses 2022, 14, 781. [Google Scholar] [CrossRef]
- Subudhi, S.; Rapin, N.; Bollinger, T.K.; Hill, J.E.; Donaldson, M.E.; Davy, C.M.; Warnecke, L.; Turner, J.M.; Kyle, C.J.; Willis, C.K.R.; et al. A Persistently Infecting Coronavirus in Hibernating Myotis Lucifugus, the North American Little Brown Bat. Journal of General Virology 2017, 98, 2297–2309. [Google Scholar] [CrossRef]
|
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content. |
© 2024 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/).