Molecular detection of Bartonella spp. and Rickettsia spp. in bat ectoparasites in Brazil

The family Streblidae comprises a monophyletic group of Hippoboscoidea, hematophagous dipterans that parasitize bats. Bartonella spp. and Rickettsia spp. have been reported in bats sampled in Europe, Africa, Asia, North, Central and South America. However, there are few reports on the Bartonella and Rickettsia bacteria infecting Hippoboscoidea flies and mites. While Spinturnicidae mites are ectoparasites found only in bats, those belonging to the family Macronyssidae comprise mites that also parasitize other mammal species. This study investigates the occurrence and assesses the phylogenetic positioning of Bartonella spp. and Rickettsia spp. found in Streblidae flies and Spinturnicidae and Macronyssidae mites collected from bats captured in Brazil. From May 2011 to April 2012 and September 2013 to December 2014, 400 Streblidae flies, 100 Macronyssidaes, and 100 Spinturnicidae mites were collected from bats captured in two sites in northeastern Nova Iguaçu, Rio de Janeiro, southeastern Brazil. Forty (19.8%) out of 202 Streblidae flies were positive for Bartonella spp. in qPCR assays based on the nuoG gene. Among the flies positive for the bacterium, six (18%) were Paratrichobius longicrus, seven (29%) Strebla guajiro, two (40%) Aspidoptera phyllostomatis, five (11%) Aspidoptera falcata, one (10%) Trichobius anducei, one (25%) Megistopoda aranea, and 18 (32%) Trichobius joblingi, and collected from bats of the following species: Artibeus lituratus, Carollia perspicillata, Artibeus planirostris, Sturnira lilium, and Artibeus obscurus. Six sequences were obtained for Bartonella (nuoG [n = 2], gltA [n = 2], rpoB [n = 1], ribC = 1]). The phylogenetic analysis based on gltA (750pb) gene showed that the Bartonella sequences clustered with Bartonella genotypes detected in bats and ectoparasites previously sampled in Latin America, including Brazil. Only one sample (0.49%) of the species Trichobius joblingi collected from a specimen of Carollia perspicillata was positive for Rickettsia sp. in cPCR based on the gltA gene (401bp). This sequence was clustered with a ‘Candidatus Rickettsia andaenae" genotype detected in an Amblyomma parvum tick collected from a rodent in the southern region of Brazilian Pantanal. The sampled Macronyssidae and Spinturnicidae mites were negative for Bartonella spp. and Rickettsia spp. This study demonstrated the first occurrence of Bartonella spp. and Rickettsia spp. DNA in Streblidae flies collected from bats in Brazil.


Introduction
The order Chiroptera is the second largest group of mammals in the world, comprising approximately 20% of mammals and more than 1200 species present in all continents, except Antarctica [1]. In Brazil, about 47% of species diversity is found in urban areas [2].
Among all the ectoparasites of bats, Streblidae flies are the ones most frequently reported in the Neotropics. These ectoparasites are often found associated with bat species of the Phyllostomidae and Noctilionidae families [3,4,5]. The family Streblidae comprises a monophyletic group of Hippoboscoidea dipterans [6,7,8], characterized by adenotrophic viviparity, consisting of not only winged but also brachypterous and apterous species that are obligate hematophagous ectoparasites of bats [5,9].
The mites of the Suborder Mesostigmata belong to four families: Macronyssidae, Laelapidae, Spelaeorhynchidae, and Spinturnicidae. The species belonging to the families Macronyssidae and Laelapidae can parasitize several mammals species, including bats, whereas those of the families Spelaeorhynchidae and Spinturnicidae are known to parasitize Chiroptera exclusively [10,11].
Bartonella species includes Gram-negative facultative intracellular alpha-proteobacteria belonging to the order Rhizobiales [12]. These reemerging agents parasitize the erythrocytes and endothelial cells of mammals, being associated with diseases in humans and animals [13,14]. Bartonella [29]. Additionally, the role of bats as carriers of Bartonella species and genotypes with zoonotic potential has been investigated. For instance, "Candidatus Bartonella mayotimonensis", an agent associated with cases of endocarditis in humans in Iowa, USA [33], has been detected in bats in Finland [34], France, Spain [35], USA [36], and in ectoparasites (flies and fleas) in Finland [33].
The genus Rickettsia includes obligatory intracellular Gram-negative bacteria belonging to the Phylum Proteobacteria, Class Alphaproteobacteria, Order Rickettsiales and Family Rickettsiaceae. The pathogenic Rickettsia species, which causes the group of diseases known as rickettsioses, are divided into two groups: Typhus, which comprises species mainly transmitted by fleas, and the Spotted Fever, that include Rickettsia species transmitted mostly by ticks [37]. Rickettsia spp. has been already detected in bats sampled in the United States [38], Saint Kitts islands, Galapagos [39], South Africa, Swaziland [29], and Argentina [40]. In Brazil, serological evidence of exposure to Rickettsia spp. (9.5% to R. rickettsii, 9.5% to R. parkeri, 7.8% to R. amblyommii, and 1.1% to R. rhipicephali) has been reported among bats sampled in São Paulo state [41].
Although Bartonella spp. and Rickettsia spp. have not been detected in Macronyssidae and Spinturnicidae mites parasitizing bats so far, these agents were detected in Ornithonyssus bacoti, a Macronyssidae mite species found parasitizing rodents in Egypt [42].
Furthermore, Rickettsia spp. has been molecularly detected in ticks collected from bats in the United States [38,43], France [44], French Guiana [45], and Poland [46]. Additionally, Rickettsia sp. of the Spotted Fever Group was detected in flies collected from bats sampled in the USA [47] and Malaysia [48].
All DNA samples were initially tested in duplicates. All duplicates whose Cq difference was higher than 0.5 were re-tested in triplicate. Amplification efficiency (E) was calculated from the slope of the standard curve in each run using the following formula (E = 10 −1/slope ). The standard curves generated by 10-fold dilutions were used to determine the amount of DNA that could be detected with 95% of sensitivity [66].
The products of all cPCR assays were separated by electrophoresis on a 1% agarose gel stained with ethidium bromide (Life Technologies™, Carlsbad, CA, USA) under 100V/150mA for 50 minutes. The gels were imaged under ultraviolet light (ChemiDoc MP Imaging System, Bio Rad™) using the Image Lab Software Version 4.1.

Sequencing and phylogenetic analyses
Amplified products were purified using the Silica Bead DNA gel extraction kit (Thermo Fisher Scientific™, Waltham, MA, USA) and sequenced using the BigDye™ Terminator v3.1 Cycle Sequencing Kit (Thermo Fisher Scientific™, Waltham, MA, USA) and the ABI PRISM 310DNA Analyzer (Applied Biosystems™, Foster City, CA, USA) [78]. The primers used in the sequencing reactions have been previously described in PCR assays for Bartonella spp.
The sequences obtained from positive samples were first submitted to a screening test using Phred-Phrap software version 23 [79,80] to evaluate the electropherogram quality and to obtain consensus sequences from the alignment of the sense and antisense sequences. The BLAST program [81] was used to analyze the sequences of nucleotides (BLASTn), to browse and compare with sequences from an international database (GenBank) [82]. The consensus sequences obtained in this study and those retrieved from GenBank were aligned using the Clustal/W software [83] via Bioedit v. 7.0.5.3 [84]. Phylogenetic inference was based on Bayesian Inference (BI) and Maximum Likelihood (ML) methods. The Bayesian inference (BI) analysis was performed with MrBayes 3.1.2 [85] via CIPRES Science Gateway [86]. Markov Chain Monte Carlo (MCMC) simulations were run for 10 6 generations with a sampling frequency of every 100 generations and a burn-in of 25%. The Maximum-likelihood (ML) analysis was inferred with the W-IQ-Tree tool available online (http://iqtree.cibiv.univie.ac.at/) [87,88] using 1000 bootstrapping replicates. The best evolution model was selected by the program jModelTest2 (version 2.1.6) on XSEDE [89], under the Akaike Information Criterion (AIC) and Bayesian Information Criterion (BIC) [90]. All trees were examined in Treegraph 2.0.56-381 beta [91].
In PCR assays, one (0.49%) out of 202 flies samples was positive for Rickettsia spp. based on the gltA gene, being identified as 'Candidatus Rickettsia andeanae' after sequencing (Table 3). However, because this positive sample had low amount of rickettsial DNA, the subsequent PCR assays based on the ompA, ompB and htrA 17-kDa genes were negative, precluding additional phylogenetic inferences.

BLAST analysis and phylogenetic inference
Based on BLAST analysis, while one Bartonella nuoG sequence (GenBank accession number MG65471) was 93% identical to Bartonella sp. WD16.2 previously isolated from a deer sampled in Japan (GenBank accession number CP019781), the other GenBank accession number MG65472) was 84% identical to B. alsatica (GenBank accession number EF659935). Two Bartonella gltA sequences GenBank accession numbers MG551538) were 93-98% identical to  Table 3). The only detected gltA-Rickettsia sp. was 100% identical to 'Candidatus Rickettsia andeanae' (GenBank accession number MG65475), previously described in a tick sample of the species Amblyomma parvum, collected from a rodent found in the Pantanal Sul-Matogrossense, Brazil, with query coverage of 100% (Table 3).
The phylogenetic tree inferred by Bayesian analysis based on sequences of the Bartonella gltA gene formed two distinct clusters. The Bartonella sequence (GenBank accession number MG551538) detected in a Strebla guajiro specimen collected from Carollia perspicillata in Rio de Janeiro state was positioned alone in a branch but closely related to Bartonella genotypes previously detected in bats from South America, one genotype detected in a bat (Sturnira lillium) in Paraná state, Brazil (KY356753), and other genotypes detected in bats from Guatemala, Mexico, and Costa Rica, with 100% branch support. Additionally, such sequences were positioned in a larger clade related to Bartonella sequences detected in rodents sampled in Brazil and U.S.A., together with a Bartonella genotype detected in a Polygenis gwyni flea collected from a Sigmodon hispidus rodent in the U.S.A., with clade support value of 83% in BI analysis. In another cluster, the Bartonella sequence (GenBank accession number MG65470) obtained in a Strebla guajiro specimen collected from Carollia perspicillata sampled in Rio de Janeiro state was closely positioned to a genotype previously detected in a specimen of Trichobius sp. fly collected in the Dominican Republic (JX416249), together with genotypes detected in bats in Mexico (MF467776) and Costa Rica (KJ816683; KJ816678; KJ816672), with branch support value of 86% probability in BI analysis. In addition, a larger clade grouped a Bartonella genotype (KY356754) detected in a specimen of Glossophaga soricina sampled in Paraná, Brazil, and sequences previously detected in cervids and bovines, such as B. capreoli (AF293392), B. schoenbuchii (AJ278181) and B. chomelii (AY254308), with 55% clade support in the BI analysis (Fig 1).
The Bartonella rpoB sequence (GenBank accession number MG65474) obtained from a M. aranea specimen collected from C. perspicillata sampled in Rio de Janeiro state was closely related to a genotype detected in a bat (S. lillium) previously sampled in Paraná state, southern Brazil, with 100% of branch support. These two sequences were positioned in the same cluster formed by Bartonella taylorii (AF165995) and Bartonella genotypes detected in rodents (AB290276) and in a bat (Myodes rufocanus) (AB779537) from Asia, with a branch support of 93% of probability BI analysis (Fig 2).
The Bartonella nuoG sequences (GenBank accession numbers MG65471; MG65472) obtained from Strebla guajiro and Megistopoda aranea specimens collected from bats of the species Carollia perspicillata and Artibeus obscuros, respectively, in Rio de Janeiro state, were positioned in a single clade separated from the others described, with 92% clade support in BI analysis (Fig 4).
The Rickettsia gltA sequence obtained in a specimen of Trichobius joblingi collected from C. perspicillata sampled in Rio de Janeiro was closely related to a 'Candidatus Rickettsia andeanae' previously detected in A. parvum tick collected from a rodent in the wetlands of Pantanal, Brazil, with 99% with branch support of in ML analysis (Fig 5).

Discussion
Studies on bats and diseases caused by bacteria have increased worldwide due to the role of this mammal group as reservoirs, hosts, and sources of infection of several pathogens [92]. The present work reports the occurrence and molecular characterization of Bartonella spp. and Rickettsia spp. in Streblidae flies parasites of bats in two localities in Rio de Janeiro state, southeastern Brazil. Streblidae flies are strictly hematophagous ectoparasites of bats, with usually high specificity for hosts. The Streblidae fly species usually parasites a single bat species or some closely related species [3,93,94], such parasitism pattern was also observed in the present study, in which a certain Streblidae species was found parasitizing no more than two different bat species belonging to the same genus or family.
Although the real role of Streblidae flies in the transmission of Bartonella spp. has not yet been confirmed, previous studies suggest that these dipterans may play an important role as invertebrate hosts for this group of pathogens, harboring a large diversity of Bartonella genotypes [11]. In the present study, two different Bartonella genotypes were observed in two Strebla guajiro specimens of the same Streblidae species, both collected from bats of the species Carollia perspicillata, which were placed in different clades in the phylogeny based on the gltA gene.
The Bartonella spp. was observed especially in Streblide flies collected from bats of the species Carollia perspicillata and Sturnira lilium, which have been previously recognized as hosts of the new Bartonella genotypes in Brazil [30]. However, a previous study with bats in Brazil reported the occurrence of Bartonella spp. (5.28%) lower than that found in this study (19.8%). This result corroborates the hypothesis that hemoconcentration occurs in the digestive tract of arthropods, which could improve the molecular diagnosis sensitivity of Bartonella. Thus, molecular assays performed on arthropods collected from hosts could reflect a more sensitive epidemiological model [95].
Recent studies with species of Nycteribiidae flies collected from bats in Madagascar aimed to relate bacterial ecology, transmission routes and host-vector specificity [48]. According to Wilkinson et al. [48], certain Bartonella genotypes and Nycteribiidae fly species may form mutualistic interactions, which may lead to host specificity. In the aforementioned study, although the found Bartonella genotypes were allocated in five different groups, an interchange of Bartonella genotypes was observed between Cyclopodia dubia and Basilia sp., Nycteribiidae flies that did not share the same bat species as hosts. The authors have suggested the existence of direct or indirect mechanisms among the vertebrate hosts that could lead to the intra-specific diversity of Bartonella observed in this family of ectoparasites [48]. Similarly, in this study, the phylogeny based on gltA gene showed that the Bartonella genotype detected in the Strebla guajiro specimen collected from Carollia perspicillata was closely related to a Bartonella genotype detected in a specimen of Trichobius sp. collected from Phyllonycteris poeyi, a bat species restrictedly distributed in Central America [96]. Therefore, Streblidae flies could act as interchangers of different Bartonella genotypes among their vertebrate hosts, leading to intraspecific diversity.
In conclusion, the phylogenetic inference based on gltA sequences also demonstrated that one of the Bartonella genotypes detected in a S. guajiro specimen collected from C. perspicillata was closely related to Bartonella genotypes previously detected in bats from Latin America. Additionally, this same Bartonella genotype also clustered with sequences previously detected in rodents sampled in the USA [97] and Brazil [98]. Similarly, the phylogenetic inference based on rpoB sequences also demonstrated that one Bartonella genotype detected in a S. guajiro specimen collected from C. perspicillata was closely related to Bartonella genotypes previously detected in bats from Brazil and Japan, and to Bartonella sp. detected in rodents from  Bartonella genotypes found in bats in South Africa and Swaziland also clustered with those detected in rodents sampled in Africa, with a low clade support.
The phylogenetic inference based on the gltA gene also showed a relationship between Bartonella genotypes detected in Strebidae and bats in Latin America with Bartonella species found in ruminants. Previously, a Bartonella genotype detected in a Carollia perspicilata specimen was closely related to a clade containing sequences of B. chomelli and B. schoenbuchensis, also isolated from ruminants, in a phylogenetic analysis based on the ftsZ gene [30].
Although with low clade support, Rickettsia genotypes detected in bats in South Africa and Swaziland were previously grouped with Rickettsia conorii [28], the causative agent of Mediterranean spotted fever [99] that has recently been detected in Rhipicephalus sanguineus ticks collected from rodents in Nigeria [100]. Similarly, the Rickettsia genotype detected in a Trichobius joblingi specimen collected from C. perspicillata in this study was phylogenetically related to the 'Candidatus Rickettsia andenae' detected in a Amblyomma parvum tick found parasitizing a rodent trapped in the Brazilian Pantanal [101]. A. parvum is a tick species that parasitizes several mammal species during its life cycle while the adult tick parasitizes mainly medium and large mammals (ruminants, equids, and carnivores), the larva and nymph are frequently found in small animals [102]. 'Candidatus Rickettsia andeanae´, whose zoonotic potential remains unknown [103], has been reported infecting ticks in Peru (Amblyomma maculatum and Ixodes bolivensis), Argentina (A. parvum) [104], and Paraguay (A. parvum) [105]. This agent was also detected in A. parvum collected from horses in the Pantanal biome in Brazil [103], in A. parvum and Amblyomma auricularium collected from horses and Turdus amaurochalinus in Northeast Brazil [103], and in Amblyomma sculptum collected from a wild animal in Mato Grosso, in central-western Brazil [106; 107]. More recently, 'Candidatus Rickettsia andeanae' was detected in A. parvum ticks collected from rodents in the wetlands of Pantanal, Brazil [101].
Although the occurrence of Bartonella and Rickettsia has not been previously reported in mites of the family Spinturnicidae, Bartonella spp. and Rickettsia have been molecularly detected in Macronyssidae mites collected from rodents in Egypt [42]. In the aforementioned study, BLAST analysis showed 81% identity with Bartonella sp. SE-BartB detected in a flea in Egypt. Regarding Rickettsia, the genotypes obtained in the study showed 100% identity with those previously detected in fleas in the U.S. and Egypt [108; 109]. However, the possible role of Macronyssidae and Spinturnicidae mites acting as reservoirs and vectors of Bartonella spp. and Rickettsia spp. among bats is still unknown.
Even though bat ectoparasites (flies, fleas, and mites) have not been found parasitizing rodents so far, the hypothesis of ticks parasitizing both mammal groups due to their low specificity in relation to their hosts, cannot be ruled out considering the high specificity between these arthropods and bats [3; 4; 93]. Ornithodoros mimon, an Argasid tick species described parasitizing bats in South America [110; 111; 112], has already been found in rodents in Brazil [113]. Landulfo et al. [114] simulated the life cycle of this tick species in laboratory conditions, using rabbits and rodents as hosts. The authors found a feeding pattern of O. mimon larval stage similar to that found in bats, demonstrating that this tick species can parasitize both rodents and bats. In addition, the occurrence of rodent ectoparasites in bats cannot be ruled out. This fact could explain the phylogenetic association between Bartonella and Rickettsia genotypes found in bats and rodents. The parasitism of bats by immature stages of A. parvum infected with 'Candidatus Rickettsia andeanae', for example, could explain the occurrence of this Rickettsia species in Streblidae flies. Mutual association between bats and rodents in the same habitat, such as caves, could provide ecological opportunities for exposure and sharing various ectoparasites and pathogens [48].
Finally, it is highlighted that the Bartonella genotypes detected in bat ectoparasites in this study were closely related to those previously detected in rodents and bats in Brazil [29; 96]; additionally, the detected Rickettsia genotype was shown to be closely related to 'Candidatus Rickettsia andenae' detected in a tick collected from a rodent in Brazil. In addition, further studies on the vector capacity of Streblidae dipterans in the transmission of Bartonella and Rickettsia among bats are needed, since 75% of emerging infectious diseases comprises zoonosis, and most of them are transmitted by arthropod vectors [115]. The increase of ecotourism in caves in Brazil associated with the fact that 47% of bat species diversity is found in urban areas [2] emphasize the need of further studies on bacterial zoonotic agents circulating in bats and ectoparasites.
Finally, the results of this study raise an interesting question about the phylogenetic relationship between the Bartonella spp. genotypes found in Streblidae flies according to the criteria defining a cut-off point for the Bartonella species based on the sequence identity of five gene regions (16S rRNA, gltA, groEL, rpoB, ftsZ, ribC) and of 16S-23S Intergenic spacer (ITS) previously established by La Scola et al. [116]. The low identity of the sequences with others previously described in GenBank allows suggesting that the genotypes found belong to a new Bartonella species circulating in bat ectoparasites, but phylogenetically close to those found in bats, rodents and ruminants.
To the best of authors' knowledge, the present work presents the first evidence of Bartonella and Rickettsia DNA in Streblidae flies collected from bats in Brazil. Future studies to evaluate the role of Streblidae flies as vectors for bacterial zoonotic agents in bats are desirable.