2.1 Collection of blood and ectoparasites from bats
A total of 123 northern bats (Eptesicus nilssonii) were captured in Yakumo City (42°29`N, 140°18`E), Hokkaido Prefecture, located in the northern part of Japan. Bat samples were separately collected in August 2017 (N = 49), June 2018 (N = 38), and August 2018 (N = 36). The bats were living within the exterior walls of an abandoned building in dense colonies as shown in the supplementary video. Before performing the present study, we obtained individual permissions to capture bats from the Oshima general sub-prefectural bureau, Hokkaido government (license #: Oshima 27, 28, 29, 122, 123, and 124).
Blood samples were aseptically collected from the bats via heart puncture using the same procedures as previously reported [21]. The blood samples and the carcasses were immediately sent to the Laboratory of Veterinary Public Health at the Department of Veterinary Medicine, College of Bioresource Sciences, Nihon University under frozen conditions with dry ice. All of the samples were stored at − 70°C until they were examined and tested.
A total of 174 bat fleas and two batbugs were obtained from the surfaces of the carcasses in the laboratory. The ectoparasites were morphologically identified at the species level under a stereo microscope SZX16 (Olympus, Tokyo, Japan) with taxonomic keys [26, 27]. Rates of parasitism by species were compared using chi-square tests with P values < 0.01 and considered to be statistically significant.
2.2 Isolation of Bartonella bacteria from bats
Isolation of Bartonella bacteria from the bat blood samples exactly followed a previously described method [21]. Bartonella bacteria were tentatively identified by colony morphology (small, white, round shape) and three colonies per sample were sub-cultured on a fresh blood agar plate using the same conditions as the primary culture.
2.3 Detection of Bartonella DNA from ectoparasites
To avoid bacterial and fungal contamination on the surface of the ectoparasites, each sample was sterilized by immersing in 500 µl of 70% ethanol containing 0.35% povidone-iodine (Shionogi & Co., Ltd, Osaka, Japan) for 10 minutes. Then, the samples were once washed with 1 ml of 0.01 M PBS with 0.5% FBS and transferred into shatter-resistant 2.0 ml tubes (SSIbio, Lodi CA, USA). After adding 400 µl of sucrose phosphate glutamate (SPG; 10 mM sodium phosphate, 220 mM sucrose, and 0.5 mM L-glutamic acid), each sample was homogenized using a bead crusher µT-12 (TAITEC corp., Saitama, Japan) at 3,000 rpm for 1 min. Genomic DNA was extracted from the homogenate aliquot (100 µl) by InstaGene Matrix (Bio-Rad Laboratories, Inc. Hercules, CA, USA).
Real-time PCR targeting the transfer-mRNA (ssrA) gene of the genus Bartonella was used for screening Bartonella DNA [28]. Real-time PCR assays were performed in 20 µl reaction mixtures containing 10 µl of TB Green Premix Ex Taq II (Takara Bio Inc., Shiga, Japan), 2 µL of DNA sample (approximately 10 to 20 ng/µl), 1 µl of each primer (10 nM), and 6 µl of nuclease-free water. The PCR conditions were as follows: 95°C for 30 seconds, followed by 40 cycles at 95°C for 5 seconds and 60°C for 30 seconds. The genomic DNA from the northern bat isolate (EN2-1) was used as a positive control and nuclease-free water as a negative control in this study. The targeted DNA was amplified with a Thermal Cycler Dice Real Time System II (Takara Bio Inc). Melting curve analysis was applied to the results with a Ct value lower than 35 cycles. When a sample showed 80 ± 1°C melting temperature, the sample was defined as positive by the assay. Each positive sample was additionally analyzed by conventional PCR targeting the citrate synthase gene (gltA) [29]. The conventional PCR products were separated on 3% agarose gels by electrophoresis and were visualized by staining with ethidium bromide under UV light. When samples tested positive for both PCRs, the samples were defined as being infected with Bartonella bacteria. A band showing the expected size for gltA (approximately 380 bp) was purified by using the Wizard SV Gel and PCR Clean-Up System (Promega, Madison, WI, USA) following the manufacturer’s instructions.
The purified PCR products of the gltA gene were directly sequenced using a BigDye Terminator Cycle Sequencing Ready Reaction kit and a Genetic Analyzer model 3130 (Applied Biosystems). When double peaks were detected in the DNA sequencing chromatograms, the PCR products were cloned using the plasmid pGEM-T Easy vector system (Promega) and resequenced. The plasmids were purified from the transformed cell lysates with a plasmid purification kit (PureYield Plasmid Miniprep System; Promega) and then were sequenced using the primers SP6 (5’-CAAGCTATTTAGGTGACACTATAG-3’) and T7 (5’-TAATACGACTCACTATAGGG-3’).
2.4 Isolation of Bartonella bacteria from ectoparasites
The homogenates of Bartonella DNA-positive ectoparasites were subsequently used for the isolation of Bartonella bacteria following the previously described method [13].
2.5 PCR amplification of the gltA and rpoB genes from Bartonella isolates
Genomic DNA was extracted from each sub-cultured colony by using InstaGene Matrix (Bio-Rad Laboratories) and subjected to genus-specific PCR targeting the gltA [29] and RNA polymerase beta-subunit-encoding (rpoB) genes [30]. The PCR products were separated on 2% agarose gels by electrophoresis and visualized by staining with ethidium bromide under UV light. Any samples showing the expected band sizes for gltA and rpoB (approximately 850 bp) were considered as members of the genus Bartonella. The PCR products were then purified using the Wizard SV Gel and PCR Clean-Up System (Promega) and directly sequenced as previously described [21].
2.6 Genotyping.
Genotypes of the gltA genes were determined when unique sequence variants with ≥ 1 nucleotide difference were found in each isolate by comparing the gltA sequences using Genetyx software Ver 12 (Genetyx Corporation, Tokyo, Japan). Representative isolates from each genotype were further analyzed by sequencing their rpoB genes.
2.7 Sequence homology analysis
The gltA and rpoB sequences of representative strains were compared with genomic sequences of prokaryotes registered in the GenBank/EMBL/DDBJ database using the BLAST program.
2.8 Phylogenetic analysis
Based on evolutionary model selection using JModelTest2 [31] with Akaike’s information criterion corrected for finite sample sizes (AICc) [32], the generalized time-reversible substitution model with four gamma-distributed categories and a proportion of invariant sites (GTR + G + I) model was the best available model for the phylogenetic analyses based on the gltA sequences.
A phylogenetic tree of the gltA sequences was constructed using the maximum-likelihood method based on the GTR + G + I model in MEGA 7 [33]. Known Bartonella species (N = 40), bat-associated Bartonella strains (N = 420) derived from 36 countries, and Brucella melitensis 16M as an outgroup were included in this analysis. Strain names, host types, host scientific names, countries where the bats and ectoparasites were collected are summarized together with references and the gltA accession numbers in supplementary Table 1.