Comparative population genetic structure of two ixodid tick species (Ixodes ovatus and Haemaphysalis flava) in Niigata Prefecture, Japan

Ixodid tick species such as Ixodes ovatus and Haemaphysalis flava are essential vectors of tick-borne diseases in Japan. In this study, we investigated the population genetic structures and gene flow of I. ovatus and H. flava as affected by the tick host mobility. We hypothesized that I. ovatus and H. flava may have differences in their genetic structure due to the low mobility of small rodent hosts of I. ovatus at the immature stage in contrast to the mediated dispersal of avian hosts for immature H. flava. We collected 307 adult I. ovatus and 220 adult H. flava from 29 and 17 locations across Niigata Prefecture, Japan. We investigated the genetic structure at two mitochondrial loci (cox1, 16S rRNA gene). For I. ovatus, pairwise FST and analysis of molecular variance (AMOVA) analyses of cox1 sequences indicated significant genetic variation among populations. Both cox1 and 16S rRNA markers showed non-significant genetic variation among locations for H. flava. The Bayesian tree and haplotype network of cox1 marker for I. ovatus samples in Niigata Prefecture found 3 genetic groups wherein most haplotypes in group 2 were distributed in low altitudinal areas. When we added cox1 sequences of I. ovatus from China to the phylogenetic analysis, three genetic groups (China 1, China 2, and Niigata and Hokkaido, Japan) were formed in the tree suggesting the potential for cryptic species in the genetic group in Japan. Our results support our hypothesis and suggest that the host preference of ticks at the immature stage may influence the genetic structure and gene flow of the ticks. This information is vital in understanding the tick-host interactions in the field to better understand the tick-borne disease transmission and in designing an effective tick control program.