Introduction
The emergence of zoonotic diseases represents a substantial threat to global biodiversity and human health [Reference Daszak, Cunningham and Hyatt1]. Due to a unique combination of life-history traits such as a long lifespan (average 10–20 years), high survival, great population densities, roosting behaviour and the ability to fly, bats are increasingly considered to be ideal reservoir hosts for emerging zoonotic pathogens [Reference Calisher2–Reference Allocati5]. Bats can roost in close contact with both domestic animals and humans, contaminating their local environment with guano and urine. They can indirectly transmit infectious agents to humans through intermediate hosts or arthropod vectors or directly through bat meat consumption, accidental bites and in some cases aerosol transmission [Reference Han6]. Most documented investigations involving bats and infectious agents are largely limited to viruses responsible for diseases in humans, such as rabies and other lyssaviruses, Ebola, Nipah and Marburg viruses, severe acute respiratory syndrome and Middle East respiratory syndrome coronaviruses (see reviews in 3 and 5).
Unlike viruses, bacteria in bats and their putative threat to human health remain poorly studied [Reference Allocati5]. Recently, bat ectoparasites (ticks, mites, fleas, cimicid bugs and bat flies) have been identified as potential vectors for pathogens such as Anaplasma, Bartonella, Borrelia, Rickettsia and Spiroplasma species [Reference Loftis7–Reference Reeves, Dowling and Dasch9]. High prevalence of the haemotrophic vector-borne bacteria Bartonella spp. were detected in the blood and organ tissues of several bat species, suggesting that bat ectoparasites are frequently implicated in transmitting bacteria [Reference Concannon10–Reference Lin13]. Veikkolainen et al. [Reference Veikkolainen14] also found Bartonella spp. DNA in faeces of vespertilionid bats, which might have originated either from bleeding into the bat intestine or from the insect prey of the bats, which include abundant blood-feeding bat ectoparasites. In the same study, the authors isolated B. mayotimonensis, identified as an aetiologic agent of endocarditis in humans, by culturing of peripheral blood samples of Daubenton's bats (Myotis daubentonii). Because Bartonella spp. are haemotrophic bacteria, it is unclear if viable bacteria can survive in bat faeces and be involved in the transmission to other hosts, including humans. Veikkolainen et al. [Reference Veikkolainen14] hypothesised that the faeces of blood-fed bat ectoparasites might transmit Bartonella spp. to human hosts through superficial scratching or trauma to tissue in the skin, which has yet to be confirmed experimentally. This does not exclude the possibility that humans might be exposed to pathogens from ectoparasite bites [Reference Gill15]. The presence of bacteria associated with human pathogens in bat guano has been confirmed in the Leschenault's Rousette (Rousettus leschenaultia) in India and in four insectivorous bat species in South Africa [Reference Banskar16, Reference Dietrich17]. While the presence of bacteria in bat faeces has become clearer [Reference Mühldorfer18], their prevalence in bat faeces is largely unknown, although this is important information to better understand the role of bats in the transmission of potentially pathogenic bacteria.
In the present study, we focused on a pathogenic bacterium commonly spread among wildlife in Europe [Reference Stuen19] and causing emerging diseases in humans: Anaplasma phagocytophilum which is transmitted through blood-feeding arthropod vectors. Anaplasma phagocytophilum is an obligate gram-negative intracellular bacterium that infects neutrophils. It can cause influenza-like symptoms in both animal and human hosts and on rare occasions is even a fatal condition in humans [Reference Ismail, Bloch and McBride20].
In this context, we investigated the faecal DNA prevalence of A. phagocytophilum in lesser horseshoe bat (Rhinolophus hipposideros) maternity colonies located in the Franche-Comté region (eastern France). The lesser horseshoe bat is a non-migratory species, with winter and maternity roosts located less than 30 kms apart [Reference Gaisler and Chytil21]. Females roost at the end of spring when they can form large colonies (generally 10–200 females) to give birth and raise their single offspring throughout the summer [Reference Gaisler22]. Maternity roosts are generally located in close proximity to humans, in old buildings such as churches or barns.
The aims of the present study were (i) to assess the prevalence of lesser horseshoe bats carrying A. phagocytophilum DNA in their faeces; and (ii) to explore whether bacterial DNA prevalence in bat faeces varied among colonies and was related to colony size. The tendency of bat females to share roosts, crowded together reaching large population densities, creates the potential for the frequent transmission of ectoparasites, which can transmit pathogenic microorganisms. In vespertilionid bats, adult females can host more ectoparasites than males, likely due to the clustering of many female bats in maternity colonies or to the physiological stresses of reproduction [Reference Zahn23–Reference Pearce and O'Shea25]. We thus hypothesised that roosts with larger lesser horseshoe bat maternity colonies would have higher vector-borne bacterial DNA faecal prevalence. Taking advantage of recent developments in non-invasive approaches in wildlife studies, the bats were studied exclusively through sampling and analysing 552 faecal samples from 23 maternity roosts.
Materials and methods
Ethics statement
Field sampling was carried out with the authorisation of the owners of the buildings hosting bat roosts. This study was based on droppings, no animal was captured or disturbed during the study period. No particular permit was required.
Faecal collection and storage
Faecal samples were collected from 23 lesser horseshoe bat maternity roosts located in three administrative departments of the Franche-Comté region, eastern France: Doubs (13 maternity roosts), Haute-Saône (three maternity roosts) and Jura (seven maternity roosts) (Fig. 1). All roosts were located in buildings (churches, barns) surrounded mainly by woodlands and shrublands, grasslands and arable lands [Reference Afonso26]. Sampling was conducted in June–July 2011, when individuals were abundant within the roost and before the birth of offspring. Plastic sheets were placed on the ground beneath the main clusters of bats and were left for 15 days. In each roost, 24 faeces samples were randomly collected directly from the plastic sheets, placed individually in 2 ml microtubes containing 90% alcohol and stored at room temperature.
Locations of the 23 lesser horseshoe bat maternity roosts sampled in the Franche-Comté region.
Because samples were collected without identification of the individual that excreted the faeces, each sample was genotyped to assign it to an individual. To check the expected sex structure of the colonies (i.e. essentially pregnant and lactating females and a few isolated males), all individuals were sexed using molecular typing. DNA extraction from bat faeces, genotyping and molecular sexing are detailed in [Reference Afonso26] and in the Technical Appendix 1.
Molecular detection of A. phagocytophilum in bat faeces
The faecal DNA prevalence of A. phagocytophilum was measured using molecular diagnosis with the DNA extracted from faeces. Only samples with a reliable genotype were tested for the presence of bacterial DNA. Because bacterial DNA in bat faeces was expected to be fragmented, we used polymerase chain reaction (PCR) primers amplifying short DNA sequences. The amplification of A. phagocytophilum was performed with the PCR primers Msp2 – 3F (5′ – CCAGCGTTTAGCAAGATAAGAG – 3′) and Msp2 – 3R (5′ – GCCCAGTAACAACATCATAAGC – 3′), amplifying a 334 bp fragment of the msp2 gene [Reference Massung and Slater27]. The PCR amplifications were conducted in a reaction mixture (25 µl) consisting of 2 µl of DNA extract, 1 × HotStarTaq Master Mix (Qiagen, Courtaboeuf, France), 0.2 µM of each primer and PCR-grade water. A negative control (PCR-grade water) was included for every 12 samples. The amplification cycling programme consisted of an activation step of 15 min at 95 °C followed by 40 cycles of denaturation at 94 °C for 30 s, annealing at 55 °C for 60 s and primer extension at 74 °C for 90 s. A final extension was performed at 72 °C for 10 min. The PCR products were separated and visualised using the QIAxcel device and a QIAxcel DNA high-resolution kit (Qiagen).
Samples PCR-positive for A. phagocytophilum were confirmed by sequencing. The amplified products were purified using the QIAquick PCR Purification Kit (Qiagen) according to the manufacturer's instructions. Direct sequencing of the PCR products was performed with an automated sequencer (Applied Biosystems 3130 Genetic Analyzer) using the primers employed for the PCR reactions. A homology search of the sequences generated in this study was then performed in February 2016 by conducting an online search with the NCBI Basic Local Alignment Search Tool for nucleotides (BLASTN) in the GenBank database (http://blast.ncbi.nlm.nih.gov/Blast.cgi).
Faecal DNA prevalence
Because several faecal samples could have been excreted by the same individual, each sample was assigned to one individual with the genotyping analyses. Among the 425 samples with a consensus genotype obtained for the eight microsatellite markers, 199 were unique and 79 had a genotype that was present at least twice in the dataset. Some of the individuals with at least two tests for A. phagocytophilum might have discordant diagnostics during the sampling period (15 days). The faecal DNA prevalence was then estimated by resampling to take into account that individuals could have both positive and negative PCR results. To estimate the prevalence of A. phagocytophilum in each bat colony, the number of DNA-positive individuals was divided by the total number of individuals sampled in the maternity roost. One sample was randomly selected for each individual and the faecal DNA prevalence was calculated. This procedure was repeated 100 times to estimate a mean individual faecal prevalence per maternity roost, i.e. the mean number of DNA-positive individuals among the 100 replicates divided by the number of individuals sampled in the maternity roost. This faecal DNA prevalence represents the mean proportion of individuals with bacterial DNA in their faeces. We estimated 95% binomial confidence intervals of the faecal DNA prevalence using a bootstrap procedure and the R package Hmisc.
Colony size
A statistical approach was used to evaluate whether colony size could be related to faecal bacterial DNA prevalence. Bat colony size is sometimes difficult to readily assess by visual estimation within maternity roosts. The number of bats in the roost varies daily and individuals are difficult to count when they are located in inaccessible places or when they are flying within the roost. The size of each colony was thus estimated via statistical methods, using the sequential Bayesian estimator method developed by Petit and Valière [Reference Petit and Valiere28] and implemented in the R 2.14.0 software (R Development Core Team. R: A Language and Environment for Statistical Computing, R Foundation for Statistical Computing, Vienna, Austria, 2011. http://www.R-project.org) using a script provided by Eric Petit (University Rennes 1). A bootstrap-method for correlations using a Spearman rank correlation test in the R software was used to investigate the relationships between the prevalence of A. phagocytophilum DNA and colony size.
Results
A total of 552 samples were collected from 23 maternity colonies. We obtained 425 (76%) consistent genetic profiles from 278 distinct individuals (9–20 individuals per colony, Table 1). Among the 278 individuals sampled, 213 were females and four were males (61 were undetermined due to a failure in PCR-control amplification). With colonies containing from 16 to 153 individuals (according to statistical estimates), a variable portion of the colony was sampled in each roost (the number of individuals sampled and the estimates of colony size are given Table 1).
Samples collected in 23 lesser horseshoe bat maternity roost during June–July 2011
Department, administrative French department; Roost, code name of the village hosting the maternity roost; Effective size, estimated number of individuals in the colony based on a sequential Bayesian estimator; Genotypes, number of reliable genotypes based on 24 collected faeces samples; PCR-positive, mean number of individuals positive after 100 resampling tests for A. phagocytophilum; Individuals, number of distinct individuals; Faecal DNA prevalence, mean number of individuals positive/number of distinct individuals, numbers in brackets are confidence intervals (95%) based on a bootstrap procedure.
Sequence analysis of DNA from samples PCR-positive for A. phagocytophilum revealed one unique DNA sequence (see Technical Appendix 2). The closest match in Genbank was A. phagocytophilum (GenBank accession number KX591651) with a similarity score of 100% with the sequence 11AL05 (292/292).
Anaplasma phagocytophilum DNA was detected in 63 of the 278 individuals (0.23, 95% CI 0.18–0.28). The faecal DNA prevalence of A. phagocytophilum was highly variable depending on the colony, with values between 0.00 and 0.82 (Table 1). The faecal DNA prevalence of A. phagocytophilum was significantly the highest in the Jura Department (Fig. 2). The faecal DNA prevalence was not related to the effective size of the colony estimated by a sequential Bayesian estimator and a bootstrap procedure using a Spearman's rank correlation test (P = 0.880).
Faecal DNA prevalence of Anaplasma phagocytophilum in lesser horseshoe bats sampled in three departments of the Franche-Comté region. Probabilities resulting from the pairwise comparison of proportional tests are given for each pair of departments. Significant differences in the faecal DNA prevalence of two departments are symbolised by * (P < 0.05).
Discussion
This study is the first investigation of the faecal DNA prevalence of A. phagocytophilum in a rhinolophid bat species and provides evidence that A. phagocytophilum DNA is prevalent in bat populations. DNA from Anaplasma spp. has been previously detected in mites (Spinturnix psi) collected off bats [Reference Reeves, Dowling and Dasch9]. In the present study, we report the faecal prevalence of A. phagocytophilum DNA in 23 lesser horseshoe bat colonies located in France based on the sampling of 278 individuals. The overall proportion of individuals with bacterial DNA in their faeces was high (22.7%). Such a high prevalence of bacterial DNA in lesser horseshoe bat colonies might suggest persistent infection of long-lived bats, as proposed for Bartonella sp. in other bat species [Reference Bai12]. This high prevalence could also reflect the presence of bacterial DNA in the insect preys of the bats (which can include mosquitoes), which is largely unknown.
Individual faecal DNA prevalence was highly variable among lesser horseshoe bat colonies (0–82%). One can assume that the high bacterial DNA prevalence reported in some roosts can indirectly reflect high exposure of bats to ectoparasites. Interestingly, faecal DNA prevalence was the highest in the Jura Department. According to the French National Institute for Health Survey, Jura Department is the department of the study area with the highest number of reported cases of tick bites in humans and the highest density of ticks (unpublished data).
The variability in prevalence among colonies was not related to the colony size. This result contradicts our hypothesis that roosts with larger lesser horseshoe bat colonies would have a higher prevalence of vector-borne bacteria. We speculated that if large colonies can be associated with high ectoparasite loads, as shown in the big brown bat, Eptesicus fuscus [Reference Pearce and O'Shea25], the chance that an arthropod vector carrying bacteria will parasitise a bat would be the highest in large colonies. Our finding could reflect the high complexity in the study of vector-borne microorganism transmission in bats. Furthermore, the faecal DNA prevalence of bacteria is not directly related to the proportion of infested hosts, as it may reflect both bat bleeding into the intestine or ingestion of ectoparasites carrying bacteria [Reference Veikkolainen14].
The high prevalence of A. phagocytophilum in faecal material in close proximity to humans (churches and barns in rural villages) raises one question: to what extent is high bacterial DNA prevalence in bat faeces a concern for humans? In the present study, all the maternity colonies were located in buildings, mainly in church bell towers. Sometimes, guano had spread onto the stairs leading to the bell tower or even into the church nave. The churchgoers regularly participate in the cleaning of the church, but the guano in the bell tower is removed more rarely, only when it reaches substantial quantities. Although we can assume that A. phagocytophilum cannot be transmitted by bat bites and likely not via direct contact with faeces, it is not excluded that bat ectoparasites can directly transmit bacteria to humans frequenting bell towers. However, even though anaplasmosis is an emerging disease in France, human cases remain rare [Reference Edouard29], contrasting with the high faecal DNA prevalence found in bat guano. However, in Europe, human prevalence studies on tick-borne diseases are mainly conducted on workers exposed to tick bites, such as forestry workers or farmers [Reference Chmielewska-Badora30, Reference Rigaud31]. Our study shows that surveys should also pay particular attention to churchgoers that frequent places of worship sheltering bat colonies to evaluate the risks of disease transmission by tick bites in such communities.
Even if bat guano is unlikely to directly transmit haemotrophic bacteria to other hosts, bats appear as epidemiologically connected to them. A future perspective of this work would be to investigate if blood-feeding arthropods parasitising bats would be good bioindicators of risk transmission to other hosts, including humans. Moreover, because bat guano can reflect both bat and insect preys A. phagocytophilum contamination, it could be easily used as a bioindicator of the spread of this bacterium in the environment. In the present study, we highlighted a high variability in faecal DNA prevalence of A. phagocytophilum among colonies and we assume that such data could serve as a basis for eco-epidemiological studies, for example by linking prevalence data with landscape components. Bats have been proposed as bioindicator taxa that show measurable responses to climate change, habitat loss and environmental pollution [Reference Jones32, Reference Zukal, Pikula and Bandouchova33] and we assume that they also could serve as a good indicator of bacteria circulation in ecosystems. Bats exhibit several life-history traits favourable to disease transmission: the colonial structure of their populations, close physical contact, heavy ectoparasite infestations and a long lifespan. Moreover, bat guano is an easy-to-collect material which analysis recently benefited from many technical optimisations.
To conclude, further research is needed to clearly understand the high faecal DNA prevalence of A. phagocytophilum we found in this study. The first point to elucidate is probably from which host (bat and/or blood-feeding arthropods) the DNA from A. phagocytophilum comes from. Bats with their guano positive for A. phagocytophilum DNA and showing a positive blood test and clinical signs might provide some answer elements. Because bat captures are difficult to proceed in maternity roosts for ethical reasons, one can hardly imagine collecting their body condition as well as the presence of bacteria in their blood and/or their ectoparasites. However, such a study might be considered on corpses sometimes found in roosts and could be envisaged with a large screening of arthropod-borne bacteria in both bats and their ectoparasites.
Supplementary material
The supplementary material for this article can be found at https://doi.org/10.1017/S0950268818001279
Acknowledgements
The authors sincerely thank Sébastien Roué and all the persons who participated in the field work and in genotyping analyses, especially Pierline Tournant and Pierre-Emmanuel Baurand. The authors also thank Eric Petit for providing an R script for the estimation of bat colony size and Kristin Muehldorfer and two anonymous reviewers for their helpful advices on previous versions of the manuscript. The study sites in Franche-Comté belong partly to the Long Term Ecological Research (LTER) Network ‘Zone Atelier Arc Jurassien’ (http://zaaj.univ-fcomte.fr/).
Financial support
This work was supported by the Région Franche-Comté, the Fondation pour la Recherche sur la Biodiversité (FRB) and the French Ministry of Ecology, Energy, Sustainable Development and the Sea (ITTECOP program) as part of the Graphab project managed by the USR 3124 MSHE Ledoux.
Conflict of interest
None.
