Combining deep sequencing and conventional molecular approaches reveals broad diversity and distribution of fleas and Bartonella in rodents and shrews from Arctic and Subarctic ecosystems

Abstract Background Bartonella are intracellular bacteria that are transmitted via animal scratches, bites and hematophagous arthropods. Rodents and their associated fleas play a key role in the maintenance of Bartonella worldwide, with > 22 species identified in rodent hosts. No studies have addressed the occurrence and diversity of Bartonella species and vectors for small mammals in Arctic and Subarctic ecosystems, which are increasingly impacted by invasive species and climate change. Methods In this study, we characterized the diversity of rodent fleas using conventional PCR targeting the mitochondrial cytochrome c oxidase II gene (COII) and Bartonella species in rodents and shrews (n = 505) from northern Canada using conventional PCR targeting the ITS (intergenic transcribed spacer) region and gltA (citrate synthase) gene. Metagenomic sequencing of a portion of the gltA gene was completed on a subset of 42 rodents and four rodent flea pools. Results Year, total summer precipitation the year prior to sampling, average minimum spring temperature and small mammal species were significant factors in predicting Bartonella positivity. Occurrence based on the ITS region was more than double that of the gltA gene and was 34% (n = 349) in northern red-backed voles, 35% (n = 20) in meadow voles, 37% (n = 68) in deer mice and 31% (n = 59) in shrews. Six species of Bartonella were identified with the ITS region, including B. grahamii, B. elizabethae, B. washoensis, Candidatus B. rudakovii, B. doshiae, B. vinsonii subsp. berkhoffii and subsp. arupensis. In addition, 47% (n = 49/105) of ITS amplicons had