Abstract
Ixodid ticks, such as Ixodes ovatus and Haemaphysalis flava, are important vectors of tick-borne diseases in Japan, such as Japanese spotted fever caused by Rickettsia japonica. This study describes the Rickettsia infection rates influenced by the population genetic structure of I.ovatus and H. flava along an altitudinal gradient. A total of 346 adult I. ovatus and 243 H. flava were analyzed for the presence of Rickettsia by nested PCR targeting the 17kDA, gltA, rOmpA, and rOmpB genes. The population genetic structure was analyzed utilizing the mitochondrial cytochrome oxidase 1 (cox1) marker. The Rickettsia infection rates were 13.26% in I. ovatus and 6.17% in H. flava. For I. ovatus, the global FST value revealed significant genetic differentiation among the different populations, whereas H. flava showed non-significant genetic differentiation. The cox1 I. ovatus cluster dendrogram showed two cluster groups, while the haplotype network and phylogenetic tree showed three genetic groups. A significant difference was observed in Rickettsia infection rates and mean altitude per group between the two cluster groups and the three genetic groups identified within I. ovatus. No significant differences were found in the mean altitude or Rickettsia infection rates of H. flava. Our results suggest a potential correlation between the low gene flow in I. ovatus populations and the spatially heterogeneous Rickettsia infection rates observed along the altitudinal gradient. This information can be used in understanding the relationship between the tick vector, its pathogen, and environmental factors, such as altitude, and for the control of tick-borne diseases in Japan.