Comparative Metagenomic Profiling of Symbiotic Bacterial Communities Associated with Ixodes persulcatus, Ixodes pavlovskyi and Dermacentor reticulatus Ticks

Ixodes persulcatus, Ixodes pavlovskyi, and Dermacentor reticulatus ticks inhabiting Western Siberia are responsible for the transmission of a number of etiological agents that cause human and animal tick-borne diseases. Because these ticks are abundant in the suburbs of large cities, agricultural areas, and popular tourist sites and frequently attack people and livestock, data regarding the microbiomes of these organisms are required. Using metagenomic 16S profiling, we evaluate bacterial communities associated with I. persulcatus, I. pavlovskyi, and D. reticulatus ticks collected from the Novosibirsk region of Russia. A total of 1214 ticks were used for this study. DNA extracted from the ticks was pooled according to tick species and sex. Sequencing of the V3-V5 domains of 16S rRNA genes was performed using the Illumina Miseq platform. The following bacterial genera were prevalent in the examined communities: Acinetobacter (all three tick species), Rickettsia (I. persulcatus and D. reticulatus) and Francisella (D. reticulatus). B. burgdorferi sensu lato and B. miyamotoi sequences were detected in I. persulcatus and I. pavlovskyi but not in D. reticulatus ticks. The pooled samples of all tick species studied contained bacteria from the Anaplasmataceae family, although their occurrence was low. DNA from A. phagocytophilum and Candidatus Neoehrlichia mikurensis was first observed in I. pavlovskyi ticks. Significant inter-species differences in the number of bacterial taxa as well as intra-species diversity related to tick sex were observed. The bacterial communities associated with the I. pavlovskyi ticks displayed a higher biodiversity compared with those of the I. persulcatus and D. reticulatus ticks. Bacterial community structure was also diverse across the studied tick species, as shown by permutational analysis of variance using the Bray-Curtis dissimilarity metric (p = 0.002). Between-sex variation was confirmed by PERMANOVA testing in I. persulcatus (p = 0.042) and I. pavlovskyi (p = 0.042) ticks. Our study indicated that 16S metagenomic profiling could be used for rapid assessment of the occurrence of medically important bacteria in tick populations inhabiting different natural biotopes and therefore the epidemic danger of studied foci.