Two tick species were identified in this study: H. qinghaiensis (in Maga only) and D. everestianus (in all four sites). D. everestianus was reported only in Northwestern China and Nepal [15] with an altitude of 2,600 - 4,700 m [16]. Larvae and nymphs of this tick specie often infest lagomorphs and rodents, while adult ticks usually utilize medium-large sized, modest and wild mammals as hosts, including hares, sheep, yaks, and horses [15, 16]. However, H. qinghaiensis is only reported in China [17-21], particularly prevalent in the western plateau, including the provinces of Qinghai, Gansu, Sichuan, and Tibet [21]. Its natural hosts include sheep, goats, yaks, cattle, and hares (Lepus oiostolus). All life stages of the tick could develop in sheep, goats, yaks, and cattle [21-27]. Contrary to ticks, H. qinghaiensis mostly performs its activity at low altitude. Arizha, Changxgma, and Derongma belong to the sub-frigid zone, whereas Maga village is located in the cold temperate zone. Due to the significant difference in altitude between Maga and the other three villages, H. qinghaiensis was only found in Maga.
All types of ticks were found to contain Bartonella DNA, although in varying percentages and locations. A survey regarding ticks from 16 states of the United States revealed that the overall prevalence of B. henselae in Ixodes ticks was 2.5% [28]. In Austria, Bartonella spp. (B. henselae, B. doshiae, and B. grahamii) were detected in 2.1% of I. ricinus, with the highest rate in ticks from Vienna (with an infection rate of 7.5%) and higher prevalence in adult ticks than other life stages [29]. Furthermore, a recent One Health perspective review on Bartonella indicated that the overall presence of Bartonella in ticks (combining evidence from multiple surveillance studies) was about 15% [30]. In our results, a total prevalence of 30.1% in ticks (especially in Maga, 76.8%) indicated the severity in Shiqu county.
B. melophagi, a human bacterial pathogen, was first isolated from sheep blood in 2007 [31], and the same bacteria were then isolated from the blood of two female patients with pericarditis and skin lesions in the United States of America [32]. As a result, the first report of DNA of B. melophagi detected in D. everestianus and H. qinghaiensis was obtained, which was the first molecular evidence of B. melophagi in Shiqu county. However, there is no current evidence supporting the ability of these ticks to transmit B. melophagi to livestock or human. To address this issue, experiments should be performed to assess the vector competency of D. everestianus and H. qinghaiensis to transmit B. melophagi in the future.
Bartonella infection has been mostly reported in Rodentia [33-42], with few cases reported in Lagomorpha. Until now, there has been only one report of Bartonella infection in plateau pika with a positive rate of 18.99% [43]. A total of 15 Bartonella strains were obtained, and most of them were closely related to B. taylorii and B. grahamii [43]. Based on our research, B. grahamii, a pathogenic strain in humans, was detected in all four villages, while B. queenslandensis was detected only in Maga. Nevertheless, similar to B. coopersplainsensis, the zoonotic potential of B. queenslandensis has not been reported. Additionally, for two unvalidated Bartonella species (Bartonella.sp* and Bartonella.sp**) found in Ariza and Changxgma respectively, sequences analysis showed that: 1) based on gltA gene, they were clustered with B. rochalimae and B. queenslandensis respectively and 2) based on rpoB, however, they were clustered with B. melophagi. There are possible explanations for why this conflicting result may occur. For instance, the potential presence of multiple Bartonella species in the sample although it is not commonly based on culturing. Secondly, different primer sets may also have amplification bias towards particular species based on the annealing affinity, which may cause the observed Bartonella diversity to differ, depending on the primer sets used for amplification. Lastly, homologous recombination, a specific form of LGT (lateral gene transfer) among Bartonella spp., have been reported in other studies based on sequencing multiple protein-coding loci [12, 44-48]. However, culturing, sequencing multiple loci (including 16S rRNA, ftsZ, gltA, groEL, ribC and rpoB, and ITS), cloning sequences into vectors before sequencing, or implementing deep sequencing approaches may discover a potentially novel Bartonella specie or subspecies and may differentiate these possible scenarios.
In Shiqu, plateau pika, the largest population of local small rodents, is in close contact with local people and livestock and can be infested with fleas and ticks, implicating them in transmission cycles of Bartonella spp. In China, Bartonella infections among humans have been mainly reported in the central plain area such as Jiangsu, Zhejiang, Anhui, and Hubei province. No cases or suspected cases have been reported in the Qinghai-Tibetan plateau. Therefore, the relationship of plateau pika and the transmission of Bartonella should be further studied. A thorough analysis with controlled experiments should be conducted to determine the exact routes of transmission between plateau pika, the transmission between plateau pika and their vectors, as well as the transmission from plateau pika to humans and livestock.