During the years 2008–2010 I. ricinus and I. persulcatus ticks were collected from 64 sites in mainland Estonia and on the island Saaremaa. Presence of B. miyamotoi was found in 0.9% (23/2622) of ticks. The prevalence in I. persulcatus and I. ricinus ticks differed significantly, 2.7% (15/561) and 0.4% (8/2061), respectively. The highest prevalence rates were in found South-Eastern Estonia in an area of I. persulcatus and I. ricinus sympatry and varied from 1.4% (1/73) to 2.8% (5/178). Co-infections with B. burgdorferi s.l. group spirochetes and tick-borne encephalitis virus were also revealed. Genetic characterization of partial 16S rRNA, p66 and glpQ genes demonstrated that Estonian sequences belong to two types of B. miyamotoi and cluster with sequences from Europe and the European part of Russia, as well as with sequences from Siberia, Asia and Japan, here designated as European and Asian types, respectively. Estonian sequences of the European type were obtained from I. ricinus ticks only, whereas the Asian type of B. miyamotoi was shown for both tick species in the sympatric regions.

Ixodes and Dermacentor ticks harbor Borrelia, Anaplasma/Ehrlichia, Bartonella, and Babesia species.

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.