Mitochondrial genome analysis across different populations reveals the intraspecific variation and phylogeography of the Caucasian soft tick relapsing fever vector, Ornithodoros (Pavlovskyella) verrucosus (Ixodida: Argasidae)

Territories in southern parts of Eastern Europe and in the Caucasus are endemic for tick-borne relapsing fever (TBRF), caused by Borrelia caucasica. This spirochete is transmitted exclusively by the bites of Ornithodoros verrucosus; however, the distribution and genetic diversity of the tick vector have not been explored. To address this, we performed a phylogeographic study of O. verrucosus specimens collected across a large geographic distribution. We sequenced and analyzed complete mitochondrial genomes of 54 individual O. verrucosus ticks representing 23 geographically diverse populations from Ukraine, Georgia, and Azerbaijan. We detected 47 unique haplotypes, with every collection site exhibiting distinct polymorphisms. This, along with other population genetic indices, suggests little evidence of gene flow between populations. The Bayesian coalescent analysis revealed the presence of four lineages that diverged in the Middle Pleistocene (770-126 kya). Two lineages were widespread and present in all study regions, while the other two were restricted to the southern foothills of the Lesser Caucasus mountain range. The sympatry of these ancient lineages suggests that isolation by environment, in addition to geographic distance, may play a role in the intraspecific divergence of tick populations. Using a phylogeographic approach, we provide a snapshot of genetic diversity in O. verrucosus and discuss the evolutionary history of the tick vector.