Molecular Survey of Vector-Borne Pathogens in Ticks, Sheep Keds, and Domestic Animals from Ngawa, Southwest China

Abstract

Vector-borne pathogens are mainly transmitted by blood-feeding arthropods such as ticks, mosquitoes, fleas, lice, mites, etc. They pose a significant threat to animal and human health due to their worldwide distribution. Although much work has been performed on these pathogens, some neglected areas and undiscovered pathogens are still to be further researched. In this study, ticks (Haemaphysalis qinghaiensis), sheep keds (Melophagus ovinus), and blood samples from yaks and goats were collected in Ngawa Tibetan and Qiang Autonomous Prefecture located on the eastern edge of the Qinghai–Tibet Plateau, Southwest China. Several vector-borne bacterial pathogens were screened and studied. Anaplasma bovis strains representing novel genotypes were detected in ticks (8.83%, 37/419), yak blood samples (45.71%, 64/140), and goat blood samples (58.93%, 33/56). Two spotted fever group (SFG) Rickettsiae, Candidatus Rickettsia jingxinensis, and a novel Rickettsia species named Candidatus Rickettsia hongyuanensis were identified in ticks. Another Rickettsia species closely related to the Rickettsia endosymbiont of Polydesmus complanatus was also detected in ticks. Furthermore, a Coxiella species was detected in ticks (3.34%, 14/419), keds (1.89%, 2/106), and yak blood (0.71%, 1/140). Interestingly, another Coxiella species and a Coxiella-like bacterium were detected in a tick and a goat blood sample, respectively. These results indicate the remarkable diversity of vector-borne pathogens circulating in this area. Further investigations on their pathogenicity to humans and domestic animals are still needed.

1. Introduction

Vector-borne diseases (VBDs) are diseases transmitted by arthropods, mainly by hematophagous arthropods such as ticks, mosquitoes, fleas, sand flies, lice, and mites. Vector-borne pathogens include bacteria (Anaplasma, Rickettsia, Coxiella, Borrelia, etc.), viruses (dengue virus, Crimean–Congo hemorrhagic fever virus, etc.), protozoans (Plasmodium vivax, Leishmania donovani, etc.), and even parasitic helminths (Wuchereria bancrofti, etc.) [1,2,3]. Due to the wide geographical distribution of hematophagous arthropods, globally emerging VBDs infect millions of people annually, especially in developing countries, and hence represent a significant threat to human health. Currently, vector-borne bacterial pathogens of most concern in China include Rickettsia spp., Orientia spp. (family Rickettsiaceae), Anaplasma spp. (family Anaplasmataceae), and Coxiella spp. (family Coxiellaceae) [4,5]. Anaplasma phagocytophilum, A. bovis, and spotted fever group Rickettsia have been frequently detected in ticks and animals from multiple provinces in China [6,7]. Cases of human infection with these pathogens have occasionally been reported [4,8]. Notably, domestic ungulates are often infected with Anaplasma spp., causing weight loss, abortion, reduced milk production, or even death, leading to numerous economic losses in the livestock industry annually [9]. Meanwhile, many infected domestic ungulates are permanently and asymptomatically infected. They may act as reservoirs and sources of pathogens for vectors in their life cycles. Meanwhile, Coxiella burnetii, the etiologic agent of Q fever, is the most notorious pathogen of the genus Coxiella [10]. As a worldwide distributed pathogen, it has been detected in multiple animal species, including goats, cattle, horses, deer, rodents, and even birds [11,12,13]. Coxiella burnetii infection in domestic animals is usually asymptomatic but may lead to infertility, abortion, and neonatal death [14], thus resulting in tremendous economic losses. Furthermore, human Q fever caused by C. burnetii is a worldwide zoonosis that has long been under-reported and underdiagnosed due to its non-specific symptoms.

Ngawa Tibetan and Qiang Autonomous Prefecture has an area of 84,242 km². It is located in the northwest of Sichuan Province and on the eastern edge of the Qinghai–Tibet Plateau, China. The high altitude (982–6250 m, mostly 2500–4000 m) and the abundant grassland resources of Ngawa make animal husbandry the major income source for rural households. Yaks, goats, horses, and pigs are the common livestock animals in Ngawa. The frequent contact between herdsmen and livestock animals makes it possible for zoonotic and vector-borne pathogens to be easily transmitted from animals or their ectoparasites to humans. In recent years, some studies on vector-borne pathogens have been performed in some areas of the eastern Qinghai–Tibet Plateau [15,16,17]. In 2016, Yang et al. reported a novel Anaplasma species closely related to A. capra in H. qinghaiensis ticks from Gannan Tibetan Autonomous Prefecture located in the east Qinghai–Tibet Plateau [15]. In 2020, diverse Bartonella species were detected in H. qinghaiensis and Dermacentor everestianus ticks collected in Shiqu County, located on the Qinghai–Tibetan Plateau in Ganze Prefecture [16]. However, except for one study on Theileria detected in H. qinghaiensis, D. everestianus, Rhipicephalus microplus, and Melophagus ovinus reported in 2020 [17], vector-borne pathogens in Ngawa remain understudied.

The objective of our study was to update the epidemiological profile of circulating bacterial VBDs in Ngawa. In this study, we collected ticks (H. qinghaiensis), sheep keds (Melophagus ovinus), and the blood of yaks and goats in Hongyuan County of Ngawa. The presence, prevalence, and genetic characteristics of several vector-borne bacterial pathogens were studied.

2. Results

2.1. Sample Collection

From February to March 2021, a total of 419 ticks (adults and nymphs) were collected from yaks (Bos mutus), and 106 keds were collected from goats (Capra hircus) in Hongyuan County of Ngawa, China (Figure 1). All ticks were morphologically determined to be H. qinghaiensis and further confirmed by analysis of the COI gene (accession numbers: ON358162-ON358175). All keds were morphologically identified as M. ovinus and molecularly confirmed (accession numbers: ON358236-ON358247). From March to May 2021, 140 and 56 blood samples were collected from apparently healthy yaks and goats, respectively.

2.2. PCR Detection

The PCR products of the 16S rDNA gene that met the expected lengths (Rickettsia: 900 bp; Anaplasma: 900 bp; Coxiella: 600 bp) were selected and subjected to sequencing. Anaplasma bovis was detected in 37 tick samples (8.83%, 37/419), 64 yak blood samples (45.71%, 64/140), and as many as 33 goat blood samples (58.93%, 33/56) (Table 1). All sequences of the 16S rDNA gene have the highest (99.21–100%) identity to A. bovis isolate Zhouzhi-goat-29 (MH255938.1) identified from goats in Shaanxi Province, China.

In total, three Rickettsia species, including Candidatus R. jingxinensis (0.48%, 2/419), SFG Rickettsia strain tick61 (0.24%, 1/419), and Rickettsia sp. strain tick14 (0.24%, 1/419), were identified in ticks (Table 1). Their 16S genes have the highest identities of 99.92% to Candidatus Rickettsia jingxinensis strains, 100% to uncultured Rickettsia sp. Hja 192 (LC379493.1) from Ha. japonica in Japan, and 99.83% to uncultured Rickettsiales bacterium clone LG127_1 (MK616427.1) from Pompholyxophrys punicea, respectively. Three Coxiella species were detected in these samples. A Coxiella strain closely related to Coxiella sp. isolate XinXian-HL9 (MG906671.1) identified from H. longicornis in China (with a sequence identity of 99.24–99.32% for the 16S gene) was detected in 14 of the tick samples (3.34%, 14/419), 2 of the ked samples (1.89%, 2/106), and 1 of the yak blood samples (0.71%, 1/140) (Table 1). Coxiella strain tick8, closely related to Coxiella endosymbionts, was detected in a tick sample with a 99.32% identity to the Coxiella endosymbiont of Rhipicephalus turanicus (JQ480818.1) and 99.24% identity to the Coxiella endosymbiont of Rhipicephalus microplus strain Porto Alegre (KY026064.1) and the Coxiella endosymbiont of Rhipicephalus australis isolate Rhaus2 (KP994830.1). Interestingly, Coxiella-like bacterium strain goat12 was identified from one goat blood sample, and the 16S rDNA shares the highest (96.77%) identity with C. burnetii strains (strain CB_25, strain CB_30, etc.).

2.3. Amplification of Key Genes and Genetic and Phylogenetic Analysis

For further characterization of the detected bacterial strains, sequences of groEL (769 bp) and gltA (826 bp) genes were recovered from randomly selected A. bovis strains from goats, yaks, and ticks. The gltA sequences have 92.23–92.74% identities to uncultured Anaplasma sp. clone 499 (JN588561.1) from raccoons in Japan. The groEL sequences have 95.54–96.49% similarities to A. bovis isolate Zhouzhi-goat-29 and 94.93–95.32% similarities to uncultured Anaplasma sp. clone 499 (JN588562.1) (Table S1). These results suggest that these strains are divergent from previously reported genotypes of A. bovis. As shown in the phylogenetic trees, the gltA and groEL genes are divided into two clusters (Figure 2). The A. bovis strains tick103, goat6, and goat9 belong to group 1, while strains goat8, goat10, yak30, and yak33 belong to group 2. Interestingly, the gltA gene of strain goat67 belongs to group 1, whereas its groEL gene is clustered in group 2.

A longer 16S fragment (1201 bp) and the groEL (656–1053 bp), gltA (1006 bp), ompA (705 bp), ompB (376 bp), and htrA (411 bp) genes were successfully recovered for the Rickettsia strains. The groEL, gltA, ompA, and htrA genes of the Candidatus R. jingxinensis strains have the highest identities of 99.05%, 100.00%, 100%, and 99.76% to other strains (Candidatus Rickettsia jingxinensis isolate Meixian-Hl-107, isolate tick-XA172, etc.). For the SFG Rickettsia strain tick61, although the 16S gene is 100% identical to Rickettsia sp. Hja 192 (LC379493.1) [18], its groEL, gltA, ompA, and htrA genes show 98.77%, 99.90%, 99.83%, 99.76% nucleotide similarity to R. peacockii str. Rustic (CP001227.1) from Dermacentor andersoni [19], Candidatus Rickettsia principis str. douglasi 061 (AY578115.1) from Ha. japonica [20], uncultured Rickettsia sp. clone Y27-1 (KT921894.1) from H. qinghaiensis [21], and uncultured Rickettsia sp. Hf_01 (LC656405) from Ha. flava in Japan, respectively (Table S2). As to the ompB gene, it has only as low as 97.07% (365/376 nt) nucleotide identity to Candidatus Rickettsia tasmanensis strain T152 (GQ223393) identified from Ixodes tasmani in Australia [22]. In the phylogenetic trees, its groEL, gltA, ompA, ompB, and htrA genes are apparently divergent from other SFG Rickettsia species (Figure 3). Notably, its ompA gene is in a basal location in the phylogenetic tree but far from SFG Rickettsia species. When only compared to validated Rickettsia species, its 16S, gltA, groEL, ompA, and ompB genes have the highest identities (99.58%, 98.81%, 98.77%, 81.87%, and 96.81%) to R. conorii strain Malish7 (NR074480) [23], R. raoultii strain Khabarovsk (CP010969) [24], R. peacockii str. Rustic (CP001227.1) [19], R. japonica strain YH_M (AP017602) [25], and R. conorii subsp. caspia A-167 (AF123708) [26], respectively. According to the criteria established by Fournier et al. [27], these results suggest that it represents a novel Rickettsia species. Based on the location where it was first characterized, we named it Candidatus R. hongyuanensis. For the Rickettsia sp. strain tick14, its groEL sequence has 99.54% (648/651 nt) nucleotide identity to the Rickettsia endosymbiont of Polydesmus complanatus strain Bel-53 (MH618395.1) identified in millipede [28] and 80.54% nucleotide identity to the Rickettsia endosymbiont of Bemisia tabaci (EU435143.1) identified in whitefly [29].

RpoB (408–526 bp) and groEL (545 bp) genes were obtained from the Coxiella strains. For the Coxiella strain closely related to Coxiella sp. isolate XinXian-HL9, its rpoB sequence has 86.96% identity (400/460, 87% coverage) to Coxiella sp. strain CoxAsp (MT000811.1) from Amblyomma sp. and 83.88% identity (385/459, 86% coverage) to the Coxiella endosymbiont of Dermacentor marginatus isolate DmarOD1 (MK248730.1). Meanwhile, its groEL gene has the highest (91.56%) similarity (499/545 nt) to Coxiella sp. isolate DR275 (MG860511.1) from Dermacentor reticulatus in Slovakia (Table S3). As shown in the phylogenetic trees, its genes form distinct clusters from any other Coxiella species (Figure 4). For the Coxiella strain tick8, closely related to Coxiella endosymbionts, only the rpoB gene was obtained. The rpoB gene shows 99.51% similarity (403/405 nt) to the Coxiella endosymbiont of Rhipicephalus annulatus isolate BF12 (KY678165.1) and 98.77% similarity (400/405 nt) to the Coxiella endosymbiont of Rhipicephalus microplus strain Porto Alegre (KY026067.1). The rpoB and groEL genes were not recovered from the Coxiella-like bacterium detected in the goat sample.

All of these sequences have been submitted to Genbank (Shown in Table S4).

3. Discussion

In this study, most pathogens were detected in H. qinghaiensis ticks. Meanwhile, some pathogens detected in ticks were also detected in goats, yaks, and keds, indicating the possible role of H. qinghaiensis in the maintenance and transmission of these pathogens in this area. Haemaphysalis qinghaiensis is an endemic tick species widely distributed in the western plateau of China [15]. As a three-host tick, it often infests domestic animals, including goats, cattle, yaks, horses, etc. Occasionally, it also bites humans [15]. In previous studies, it has been demonstrated that H. qinghaiensis transmits a lot of human and animal pathogens, such as Theileria sp., Babesia sp., Borrelia burgdorferi, and A. phagocytophilum [21]. The potential role of H. qinghaiensis ticks in the transmission of human and animal diseases in this area should be further considered.

Anaplasma bovis has been well recognized as a tick-borne pathogen infecting cattle, buffalos, yaks, goats, and rodents, as well as humans [4,30,31,32]. It is one of the etiologic agents of bovine anaplasmosis, characterized by fluctuating fever, depression, and death [33]. Occasionally, it also infects humans and causes fever, headache, myalgia, rash, eschar, etc. [4]. In this study, we observed remarkably high positivity rates of A. bovis in apparently healthy goats (58.93%) and yaks (45.71%), indicating that A. bovis is circulating in domestic animals in this region. Genetic and phylogenetic analyses both indicate that these strains are different from previously reported A. bovis strains, and they may represent more than one novel genotype. Phylogenetic analysis based on groEL and gltA genes indicated that they are divided into two clusters, and genetic recombination between the two clusters may exist, suggesting the long-history evolution of A. bovis in this area. The pathogenicity of these strains to humans is still to be determined. Ticks also tested positive for this A. bovis strain, with a positivity rate of 8.83%. As reported previously, the detection rate of Anaplasma in ticks was likely skewed due to Anaplasma infection in livestock [33]. Although the possibility cannot be ruled out that the A. bovis DNA may come from the blood meal of ticks, we suspect that H. qinghaiensis ticks might play a role in the transmission of A. bovis in domestic animals.

Three Rickettsia species, including Rickettsia sp. strain tick14, Candidatus R. jingxinensis, and a novel SFG Rickettsia species named Candidatus R. hongyuanensis, were identified in ticks. Candidatus R. jingxinensis, belonging to SFG Rickettsia, was first identified in Ha. longicornis from northeast China in 2016 [34]. In recent years, it has been detected in Rhipicephalus microplus, Ha. turturis, and Ixodes persulcatus from multiple locations in China, Korea, Thailand, and India, indicating that it is a widely distributed Rickettsia circulating in Asia [35,36,37,38,39]. Notably, a gltA sequence (KU853023) from a patient in China was submitted to GenBank (unpublished), which shows 99.91% nucleotide identity to Ca. R. jingxinensis, suggesting its potential human pathogenicity. Furthermore, a novel SFG Rickettsia named Candidatus R. hongyuanensis was identified in H. qinghaiensis ticks. Its 16S, gltA, groEL, ompA, and ompB genes have the highest identities of 99.58%, 98.81%, 98.77%, 81.87%, and 96.81% to validated species, supporting that it represents a novel Rickettsia species based on criteria established by Fournier et al. [27]. Genetic analyses indicate that it is closely related to human pathogenic R. raoultii and R. conorii. Due to the fact that H. qinghaiensis is a three-host tick, it can easily transmit pathogens between different hosts during its life cycle [15]. Therefore, the role of H. qinghaiensis in the transmission of SFG Rickettsia should be further considered, and investigations on the human pathogenicity of these Rickettsia species are still needed. Furthermore, Rickettsia sp. strain tick14, closely related to the Rickettsia endosymbiont of Polydesmus complanatus Bel-53 [28], was identified, and it is located in a basic position of the phylogenetic trees based on 16S and groEL sequences.

Similar to Rickettsiales bacteria, the genus Coxiella, belonging to Coxiellaceae, is a group of intracellular bacteria [40]. This genus includes two recognized species (C. burnetii and C. cheraxi), one Candidatus species (Candidatus Coxiella mudrowiae), and some unclassified Coxiella symbionts [41,42]. In this study, three Coxiella species were detected, with two of them infecting domestic animals. This may be one of the few recent studies in which Coxiella species were found to infect mammals. One Coxiella species (strain yak17, strain tick103, and strain tick166) was detected in ticks (3.34%), sheep keds (1.89%), and one yak blood sample (0.71%). Interestingly, sheep keds also tested positive for this Coxiella. There might be three possibilities: (1). Sheep keds are infected, and they act as the vector of this Coxiella. (2). The Coxiella DNA is from the blood meal of sheep keds. That is, goats may be infected by this Coxiella. (3). The Coxiella strains could be just Coxiella-like endosymbionts that co-exist with hosts/vectors. Each possibility is intriguing, and further studies are still needed. Furthermore, it is noteworthy that a Coxiella-like bacterium was identified in a goat blood sample. Its 16S gene shared the highest (96.77%) identity with the human and animal pathogen C. burnetii strain CB_25 (LC464975.1), strain CB_30 (LC464974.1), etc. (unpublished). However, sequences of other key genes are not available for this strain, possibly due to the low identity to known strains.

In summary, we identified seven bacterial species in arthropods and domestic animals in Ngawa. Although the human pathogenicity for most of these bacteria is still to be determined, more attention should be paid to the risk of human infection and the possible circulation of these pathogens in local people.

4. Materials and Methods

4.1. Sample Collection and Processing

Ngawa Tibetan and Qiang Autonomous Prefecture is located in the northwest of Sichuan Province and on the eastern edge of Qinghai–Tibet Plateau, China. From February to March 2021, 419 ticks and 106 sheep keds were respectively collected from yaks and goats in Sedi town (102.98°E, 30.01°N, with an altitude of 3528 m), Hongyuan County of Ngawa, Sichuan Province, China. The ticks and keds were stored at −40 °C before DNA extraction. From March to May, blood samples were collected from asymptomatic domestic animals, including yaks and goats, in Hongyuan County. Whole blood samples were stored in tubes containing EDTA as anticoagulant and then stored at −40 °C before DNA extraction.

4.2. DNA Extraction

The cryopreserved tick and sheep ked samples were washed three times with phosphate-buffered saline (PBS) and then were individually placed into 2 mL sterile tubes with 0.5 mL of PBS solution. After homogenization using Mixer Mill MM 400 (Retsch, Hann, Germany), the arthropod suspension was directly subjected to DNA extraction with a Mollusc DNA Extraction Kit (Omega Bio-Tek, Norcross, GA, USA). The kit was operated according to the manufacturer’s instructions, and the genomic DNA was extracted and eluted with 60 μL of Elution Buffer twice.

According to the protocol of Blood & Tissue DNA Kit (Omega, Norcross, GA, USA), DNA was extracted from 200 μL of the animal blood samples. After extraction, the DNA samples were diluted into 60 μL of Elution Buffer repeatedly.

4.3. Species Identification of Ticks and Keds

Ticks and keds were morphologically identified by an arthropod taxonomist using microscopy and then stored at −40 °C until DNA isolation. The tick species was morphologically determined by observing the basis scapituli and palp, while the ked species was determined by checking the mouthparts, bristles on the body surface, and legs tipped with pointed claws [43]. Tick and ked species were further confirmed by PCR amplification and analysis of the COI gene (cytochrome C oxidase I) (primers shown in Table S5). PCR amplification was performed using a Sensoquest PCR System LabCycler (Göttingen, Germany). The cycling conditions were: 94 °C for 3 min, 35 cycles of 94 °C for 30 s, 55 °C for 1 min, 72 °C for 1 min, and a final extension of 72 °C for 5 min [44]. After electrophoresis, the amplicons were subjected to sequencing and analyzed by BLASTN alignment.

4.4. PCR Detection and Amplification of Key Genes

To determine the presence and prevalence of Rickettsia spp., Anaplasma spp., and Coxiella spp., semi-nested or nested PCR targeting a conserved domain of 16S rDNA gene was performed on all DNA samples. The primers are shown in Table S5 and previous reports [43,45]. PCR amplification was performed using a Sensoquest PCR System LabCycler (Germany). For further analysis of the molecular characteristics and phylogenetic positions of these bacterial strains, the sequences of two other genes (the gltA gene encoding citrate synthase and the groEL gene encoding the heat shock protein) were obtained for Rickettsia (gltA: 1006 bp; groEL: 656–1053 bp) and Anaplasma (gltA: 826 bp; groEL: 769 bp) positive samples (primers and details shown in Table S5). Additionally, a longer 16S fragment (1201 bp), ompA (705 bp), ompB (376 bp), and htrA (411 bp) genes were recovered from Rickettsia strains. The primers used were previously described [43,46] (Table S5). For Coxiella strains, the sequences of groEL and rpoB genes were obtained using the primers shown in Table S5 [45]. The details of amplification are provided in Table S5. After electrophoresis in 1.0% agarose gel, the PCR amplicons were subjected to sequencing by Sangon Biotech Company (Shanghai, China), and the sequencing results were analyzed by BLASTN alignment.

4.5. Genetic and Phylogenetic Analyses

Protein-coding sequences from Rickettsia spp., Anaplasma spp., and Coxiella spp. obtained in this study, as well as the reference sequences retrieved from GenBank, were aligned using ClustalW (“protein-coding genes” strategy) in the MEGA software, version 7.0 [47]. Nucleotide sequence identities were calculated by MegAlign program in the DNASTAR Lasergene package (DNASTAR, Inc., Madison, WI, USA). Phylogenetic trees of the 16S and other key genes based on the maximum likelihood method were reconstructed by PhyML v3.2 [48]. The substitution model test was performed to determine the best-fit phylogenetic model. The confidence values for individual branches of the tree were determined by bootstrap analysis with 100 repetitions.

Author Contributions

Conceptualization, K.L., H.Z. and H.J.; methodology, M.L.; validation, M.L. and W.W.; formal analysis, W.W. and X.Q.; resources, J.T. and H.J.; data curation, K.L.; writing—review and editing, K.L., H.Z. and H.J.; supervision, K.L.; project administration, K.L.; funding acquisition, M.L., J.T., H.Z., H.J., W.W., X.Q. and K.L. All authors have read and agreed to the published version of the manuscript.

Institutional Review Board Statement

This work was approved by the ethics committee of the National Institute for Communicable Disease Control and Prevention, the Chinese Center for Disease Control and Prevention (approval numbers: SYXK2021-011 and SYXK2017-0021), and Chengde Medical University (approval number: 2020004).

Informed Consent Statement

Not applicable.

Data Availability Statement

All sequence files are available from the NCBI database (shown in Table S4).

Acknowledgments

We sincerely thank Wenping Guo, Duoji, and Xiaojing Jin for their warmhearted help.

Conflicts of Interest

The authors declare no conflict of interest.

Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Figures and Table

Figure 1. A map showing the location of Ngawa Tibetan and Qiang Autonomous Prefecture, Sichuan Province, China.

Figure 2. Phylogenetic trees constructed by the PhyML 3.0 software based on the nucleotide sequences of 16S rRNA (881 bp), gltA (826 bp), and groEL (769 bp) genes of Anaplasma bovis strains. (A): 16S rRNA, (B): gltA, (C): groEL.

Figure 3. Phylogenetic trees constructed by the PhyML 3.0 software based on the nucleotide sequences of 16S rRNA (1201 bp), gltA (1006 bp), groEL (656–1053 bp), ompA (687–699 bp), ompB (376 bp), and htrA (411 bp) genes of Rickettsia strains. (A): 16S rRNA, (B): gltA, (C: groEL, (D): ompA, (E): ompB, (F): htrA.

Figure 4. Phylogenetic trees constructed by the PhyML 3.0 software based on the nucleotide sequences of 16S rRNA (1180–1181 bp), groEL (545 bp), and rpoB (408–526 bp) genes of Coxiella strains. (A): 16S rRNA, (B): groEL, (C): rpoB.

Table 1

The positivity rates of Anaplasma, Rickettsia, and Coxiella pathogens in different hosts (ticks, keds, yaks, and goats).

Bacterial Species	Haemaphysalis qinghaiensis	Melophagus ovinus	Bos mutus	Capra hircus
Anaplasma bovis	8.83% (37/419)	0.00% (0/106)	45.71% (64/140)	58.93% (33/56)
Candidatus Rickettsia jingxinensis	0.48% (2/419)	0.00% (0/106)	0.00% (0/140)	0.00% (0/56)
Candidatus Rickettsia hongyuanensis	0.24% (1/419)	0.00% (0/106)	0.00% (0/140)	0.00% (0/56)
Rickettsia sp. tick14	0.24% (1/419)	0.00% (0/106)	0.00% (0/140)	0.00% (0/56)
Coxiella spp.	3.34% (14/419)	1.89% (2/106)	0.71% (1/140)	0.00% (0/56)
Coxiella sp. tick8	0.24% (1/419)	0.00% (0/106)	0.00% (0/140)	0.00% (0/56)
Coxiella-like bacterium	0.00% (0/419)	0.00% (0/106)	0.00% (0/140)	1.79% (1/56)
Total	13.37% (56/419)	1.89% (2/106)	46.43% (65/140)	60.71% (34/56)

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/pathogens11050606/s1, Table S1. Nucleotide identity of 16S, gltA, and groEL genes of Anaplasma bovis strains to reported strains in the Genbank database. Table S2. Nucleotide identity of 16S, gltA, groEL, ompA, ompB, and htrA genes of Rickettsia strains to reported strains in the Genbank database. Table S3. Nucleotide identity of 16S, groEL, and rpoB genes of Coxiella strains to reported strains in the Genbank database. Table S4. Genbank numbers of Rickettsia spp., Anaplasma spp., and Coxiella spp. sequences obtained in this study. Table S5. The primers used for amplification of 16S, gltA, groEL, htrA, ompA, ompB, and rpoB genes from Rickettsia spp., Anaplasma spp., and Coxiella spp. by nested PCR or semi-nested PCR. Supporting Data. Nucleotide sequences of Candidatus Rickettsia hongyuanensis.

References

1. Smetanová, K.; Schwarzová, K.; Kocianová, E. Detection of Anaplasma phagocytophilum, Coxiella burnetii, Rickettsia spp., and Borrelia burgdorferi s. l. in ticks, and wild-living animals in western and middle Slovakia. Ann. N. Y. Acad. Sci.; 2006; 1078, pp. 312-315. [DOI: https://dx.doi.org/10.1196/annals.1374.058] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/17114728]

2. Todd, C.S.; Mansoor, G.F.; Buhler, C.; Rahimi, H.; Zekria, R.; Fernandez, S.; Mikhail, A.F.; Scott, P.T.; Yingst, S.L. Prevalence of zoonotic and vector-borne infections among Afghan National Army Recruits in Afghanistan. Vector Borne Zoonotic Dis.; 2016; 16, pp. 501-506. [DOI: https://dx.doi.org/10.1089/vbz.2015.1921] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/27304051]

3. Minetti, C.; Pilotte, N.; Zulch, M.; Canelas, T.; Tettevi, E.J.; Veriegh, F.B.D.; Osei-Atweneboana, M.Y.; Williams, S.A.; Reimer, L.J. Field evaluation of DNA detection of human filarial and malaria parasites using mosquito excreta/feces. PLoS Negl. Trop. Dis.; 2020; 14, e0008175. [DOI: https://dx.doi.org/10.1371/journal.pntd.0008175] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/32267840]

4. Lu, M.; Li, F.; Liao, Y.; Shen, J.J.; Xu, J.M.; Chen, Y.Z.; Li, J.H.; Holmes, E.C.; Zhang, Y.Z. Epidemiology and diversity of Rickettsiales bacteria in humans and animals in Jiangsu and Jiangxi provinces, China. Sci. Rep.; 2019; 9, 13176. [DOI: https://dx.doi.org/10.1038/s41598-019-49059-3] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/31511528]

5. Ni, J.; Lin, H.; Xu, X.; Ren, Q.; Aizezi, M.; Luo, J.; Luo, Y.; Ma, Z.; Chen, Z.; Tan, Y. et al. Coxiella burnetii is widespread in ticks (Ixodidae) in the Xinjiang areas of China. BMC Vet. Res.; 2020; 16, 317. [DOI: https://dx.doi.org/10.1186/s12917-020-02538-6]

6. Ge, Y.; Yin, H.; Rikihisa, Y.; Pan, W.; Yin, H. Molecular Detection of Tick-Borne Rickettsiales in Goats and Sheep from Southeastern China. Vector Borne Zoonotic Dis.; 2016; 16, pp. 309-316. [DOI: https://dx.doi.org/10.1089/vbz.2015.1884]

7. Wang, H.; Yang, J.; Mukhtar, M.U.; Liu, Z.; Zhang, M.; Wang, X. Molecular detection and identification of tick-borne bacteria and protozoans in goats and wild Siberian roe deer (Capreolus pygargus) from Heilongjiang Province, northeastern China. Parasit Vectors; 2019; 12, 296. [DOI: https://dx.doi.org/10.1186/s13071-019-3553-1]

8. Tsai, K.H.; Chung, L.H.; Chien, C.H.; Tung, Y.J.; Wei, H.Y.; Yen, T.Y.; Shu, P.Y.; Wang, H.C. Human granulocytic anaplasmosis in Kinmen, an offshore island of Taiwan. PLoS Negl. Trop. Dis.; 2019; 13, e0007728. [DOI: https://dx.doi.org/10.1371/journal.pntd.0007728]

9. Woldehiwet, Z. Immune evasion and immunosuppression by Anaplasma phagocytophilum, the causative agent of tick-borne fever of ruminants and human granulocytic anaplasmosis. Vet. J.; 2008; 175, pp. 37-44. [DOI: https://dx.doi.org/10.1016/j.tvjl.2006.11.019]

10. Kazar, J. Coxiella burnetii infection. Ann. N. Y. Acad. Sci.; 2005; 1063, pp. 105-114. [DOI: https://dx.doi.org/10.1196/annals.1355.018]

11. Cho, H.C.; Hwang, S.; Kim, E.M.; Park, Y.J.; Shin, S.U.; Jang, D.H.; Chae, J.S.; Choi, K.S. Prevalence and molecular characterization of Coxiella burnetii in cattle, goats, and horses in the Republic of Korea. Vector Borne Zoonotic Dis.; 2021; 21, pp. 502-508. [DOI: https://dx.doi.org/10.1089/vbz.2020.2764]

12. González-Barrio, D.; Hagen, F.; Tilburg, J.J.; Ruiz-Fons, F. Coxiella burnetii genotypes in Iberian wildlife. Microb. Ecol.; 2016; 72, pp. 890-897. [DOI: https://dx.doi.org/10.1007/s00248-016-0786-9]

13. Tokarevich, N.K.; Panferova, Y.A.; Freylikhman, O.A.; Blinova, O.V.; Medvedev, S.G.; Mironov, S.V.; Grigoryeva, L.A.; Tretyakov, K.A.; Dimova, T.; Zaharieva, M.M. et al. Coxiella burnetii in ticks and wild birds. Ticks Tick-Borne Dis.; 2019; 10, pp. 377-385. [DOI: https://dx.doi.org/10.1016/j.ttbdis.2018.11.020]

14. Marenzoni, M.L.; Stefanetti, V.; Papa, P.; Casagrande Proietti, P.; Bietta, A.; Coletti, M.; Passamonti, F.; Henning, K. Is the horse a reservoir or an indicator of Coxiella burnetii infection? Systematic review and biomolecular investigation. Vet. Microbiol.; 2013; 167, pp. 662-669. [DOI: https://dx.doi.org/10.1016/j.vetmic.2013.09.027]

15. Yang, J.; Liu, Z.; Niu, Q.; Liu, J.; Han, R.; Liu, G.; Shi, Y.; Luo, J.; Yin, H. Molecular survey and characterization of a novel Anaplasma species closely related to Anaplasma capra in ticks, northwestern China. Parasit Vectors; 2016; 9, 603. [DOI: https://dx.doi.org/10.1186/s13071-016-1886-6]

16. Hao, L.; Yuan, D.; Guo, L.; Hou, W.; Mo, X.; Yin, J.; Yang, A.; Li, R. Molecular detection of Bartonella in ixodid ticks collected from yaks and plateau pikas (Ochotona curzoniae) in Shiqu County, China. BMC Vet. Res.; 2020; 16, 235. [DOI: https://dx.doi.org/10.1186/s12917-020-02452-x]

17. Hao, L.; Yuan, D.; Li, S.; Jia, T.; Guo, L.; Hou, W.; Lu, Z.; Mo, X.; Yin, J.; Yang, A. et al. Detection of Theileria spp. in ticks, sheep keds (Melophagus ovinus), and livestock in the eastern Tibetan Plateau, China. Parasitol Res.; 2020; 119, pp. 2641-2648. [DOI: https://dx.doi.org/10.1007/s00436-020-06757-6]

18. Thu, M.J.; Qiu, Y.; Matsuno, K.; Kajihara, M.; Mori-Kajihara, A.; Omori, R.; Monma, N.; Chiba, K.; Seto, J.; Gokuden, M. et al. Diversity of spotted fever group rickettsiae and their association with host ticks in Japan. Sci. Rep.; 2019; 9, 1500. [DOI: https://dx.doi.org/10.1038/s41598-018-37836-5]

19. Felsheim, R.F.; Kurtti, T.J.; Munderloh, U.G. Genome sequence of the endosymbiont Rickettsia peacockii and comparison with virulent Rickettsia rickettsii: Identification of virulence factors. PLoS ONE; 2009; 4, e8361. [DOI: https://dx.doi.org/10.1371/journal.pone.0008361]

20. Mediannikov, O.; Sidelnikov, Y.; Ivanov, L.; Fournier, P.E.; Tarasevich, I.; Raoult, D. Far eastern tick-borne rickettsiosis: Identification of two new cases and tick vector. Ann. N. Y. Acad. Sci.; 2006; 1078, pp. 80-88. [DOI: https://dx.doi.org/10.1196/annals.1374.010]

21. Yang, J.; Tian., Z.; Liu, Z.; Niu, Q.; Han, R.; Li, Y.; Guan, G.; Liu, J.; Liu, G.; Luo, J. et al. Novel spotted fever group rickettsiae in Haemaphysalis qinghaiensis ticks from Gansu, Northwest China. Parasit Vectors; 2016; 9, 146. [DOI: https://dx.doi.org/10.1186/s13071-016-1423-7]

22. Izzard, L.; Graves, S.; Cox, E.; Fenwick, S.; Unsworth, N.; Stenos, J. Novel rickettsia in ticks, Tasmania, Australia. Emerg Infect Dis.; 2009; 15, pp. 1654-1656. [DOI: https://dx.doi.org/10.3201/eid1510.090799]

23. Ogata, H.; Audic, S.; Renesto-Audiffren, P.; Fournier, P.E.; Barbe, V.; Samson, D.; Roux, V.; Cossart, P.; Weissenbach, J.; Claverie, J.M. et al. Mechanisms of evolution in Rickettsia conorii and R. prowazekii. Science; 2001; 293, pp. 2093-2098. [DOI: https://dx.doi.org/10.1126/science.1061471]

24. El Karkouri, K.; Mediannikov, O.; Robert, C.; Raoult, D.; Fournier, P.E. Genome Sequence of the Tick-Borne Pathogen Rickettsia raoultii. Genome Announc.; 2016; 4, e00157-16. [DOI: https://dx.doi.org/10.1128/genomeA.00157-16]

25. Akter, A.; Ooka, T.; Gotoh, Y.; Yamamoto, S.; Fujita, H.; Terasoma, F.; Kida, K.; Taira, M.; Nakadouzono, F.; Gokuden, M. et al. Extremely Low Genomic Diversity of Rickettsia japonica Distributed in Japan. Genome Biol. Evol.; 2017; 9, pp. 124-133. [DOI: https://dx.doi.org/10.1093/gbe/evw304]

26. Roux, V.; Raoult, D. Phylogenetic analysis of members of the genus Rickettsia using the gene encoding the outer-membrane protein rOmpB (ompB). Int. J. Syst. Evol. Microbiol.; 2000; 50, Pt 4, pp. 1449-1455. [DOI: https://dx.doi.org/10.1099/00207713-50-4-1449]

27. Fournier, P.E.; Dumler, J.S.; Greub, G.; Zhang, J.; Wu, Y.; Raoult, D. Gene sequence-based criteria for identification of new rickettsia isolates and description of Rickettsia heilongjiangensis sp. nov. J. Clin. Microbiol.; 2003; 41, pp. 5456-5465. [DOI: https://dx.doi.org/10.1128/JCM.41.12.5456-5465.2003]

28. Li, K.; Stanojević, M.; Stamenković, G.; Ilić, B.; Paunović, M.; Lu, M.; Pešić, B.; Đurić Maslovara, I.; Siljic, M.; Cirkovic, V. et al. Insight into diversity of bacteria belonging to the order Rickettsiales in 9 arthropods species collected in Serbia. Sci. Rep.; 2019; 9, 18680. [DOI: https://dx.doi.org/10.1038/s41598-019-55077-y] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/31822714]

29. Gottlieb, Y.; Zchori-Fein, E.; Mozes-Daube, N.; Kontsedalov, S.; Skaljac, M.; Brumin, M.; Sobol, I.; Czosnek, H.; Vavre, F.; Fleury, F. et al. The transmission efficiency of tomato yellow leaf curl virus by the whitefly Bemisia tabaci is correlated with the presence of a specific symbiotic bacterium species. J. Virol.; 2010; 84, pp. 9310-9317. [DOI: https://dx.doi.org/10.1128/JVI.00423-10] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/20631135]

30. Guo, W.P.; Tie, W.F.; Meng, S.; Li, D.; Wang, J.L.; Du, L.Y.; Xie, G.C. Extensive genetic diversity of Anaplasma bovis in ruminants in Xi’an, China. Ticks Tick Borne Dis.; 2020; 11, 101477. [DOI: https://dx.doi.org/10.1016/j.ttbdis.2020.101477] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/32723632]

31. Panthawong, A.; Grieco, J.P.; Ngoen-Klan, R.; Chao, C.C.; Chareonviriyaphap, T. Detection of Anaplasma spp. and Bartonella spp. from wild-caught rodents and their ectoparasites in Nakhon Ratchasima Province, Thailand. J. Vector Ecol.; 2020; 45, pp. 241-253. [DOI: https://dx.doi.org/10.1111/jvec.12395]

32. Yang, J.; Liu, Z.; Niu, Q.; Liu, J.; Guan, G.; Xie, J.; Luo, J.; Wang, S.; Wang, S.; Yin, H. First molecular survey and identification of Anaplasma spp. in white yaks (Bos grunniens) in China. Parasitology; 2016; 143, pp. 686-691. [DOI: https://dx.doi.org/10.1017/S003118201600041X]

33. Von Fricken, M.E.; Qurollo, B.A.; Boldbaatar, B.; Wang, Y.W.; Jiang, R.R.; Lkhagvatseren, S.; Koehler, J.W.; Moore, T.C.; Nymadawa, P.; Anderson, B.D. et al. Genetic diversity of Anaplasma and Ehrlichia bacteria found in Dermacentor and Ixodes ticks in Mongolia. Ticks Tick Borne Dis.; 2020; 11, 101316. [DOI: https://dx.doi.org/10.1016/j.ttbdis.2019.101316]

34. Liu, H.; Li, Q.; Zhang, X.; Li, Z.; Wang, Z.; Song, M.; Wei, F.; Wang, S.; Liu, Q. Characterization of rickettsiae in ticks in northeastern China. Parasit Vectors; 2016; 9, 498. [DOI: https://dx.doi.org/10.1186/s13071-016-1764-2]

35. Bang, M.S.; Kim, C.M.; Pyun, S.H.; Kim, D.M.; Yun, N.R. Molecular investigation of tick-borne pathogens in ticks removed from tick-bitten humans in the southwestern region of the Republic of Korea. PLoS ONE; 2021; 16, e0252992. [DOI: https://dx.doi.org/10.1371/journal.pone.0252992]

36. Guo, W.P.; Wang, Y.H.; Lu, Q.; Xu, G.; Luo, Y.; Ni, X.; Zhou, E.M. Molecular detection of spotted fever group rickettsiae in hard ticks, northern China. Transbound Emerg. Dis.; 2019; 66, pp. 1587-1596. [DOI: https://dx.doi.org/10.1111/tbed.13184]

37. Liu, H.; Liang, X.; Wang, H.; Sun, X.; Bai, X.; Hu, B.; Shi, N.; Wang, N.; Zhang, X.; Huang, L. et al. Molecular evidence of the spotted fever group Rickettsiae in ticks from Yunnan Province, Southwest China. Exp. Appl. Acarol.; 2020; 80, pp. 339-348. [DOI: https://dx.doi.org/10.1007/s10493-020-00467-5]

38. Takhampunya, R.; Korkusol, A.; Pongpichit, C.; Yodin, K.; Rungrojn, A.; Chanarat, N.; Promsathaporn, S.; Monkanna, T.; Thaloengsok, S.; Tippayachai, B. et al. Metagenomic approach to characterizing disease epidemiology in a disease-endemic environment in Northern Thailand. Front. Microbiol.; 2019; 10, 319. [DOI: https://dx.doi.org/10.3389/fmicb.2019.00319]

39. Wang, Q.; Guo, W.B.; Pan, Y.S.; Jiang, B.G.; Du, C.H.; Que, T.C.; Zhan, L.; Wu, J.H.; Yu, M.H.; Cui, X.M. et al. Detection of novel spotted fever group Rickettsiae (Rickettsiales: Rickettsiaceae) in ticks (Acari: Ixodidae) in Southwestern China. J. Med. Entomol.; 2021; 58, pp. 1363-1369. [DOI: https://dx.doi.org/10.1093/jme/tjaa294]

40. Brenner, A.E.; Muñoz-Leal, S.; Sachan, M.; Labruna, M.B.; Raghavan, R. Coxiella burnetii and related tick endosymbionts evolved from pathogenic ancestors. Genome Biol. Evol.; 2021; 13, evab108. [DOI: https://dx.doi.org/10.1093/gbe/evab108]

41. Gottlieb, Y.; Lalzar, I.; Klasson, L. Distinctive genome reduction rates revealed by genomic analyses of two Coxiella-like endosymbionts in ticks. Genome Biol. Evol.; 2015; 7, pp. 1779-1796. [DOI: https://dx.doi.org/10.1093/gbe/evv108]

42. Tan, C.K.; Owens, L. Infectivity, transmission and 16S rRNA sequencing of a rickettsia, Coxiella cheraxi sp. nov., from the freshwater crayfish Cherax quadricarinatus. Dis. Aquat. Organ.; 2000; 41, pp. 115-122. [DOI: https://dx.doi.org/10.3354/dao041115]

43. Huang, B.; Shen, J. Classific Atlas of Parasites for Livestock and Poultry in China; China Agricultural Science and Technology Press: Beijing, China, 2006; (In Chinese)

44. Lu, M.; Tian, J.; Pan, X.; Qin, X.; Wang, W.; Chen, J.; Guo, W.; Li, K. Identification of Rickettsia spp., Anaplasma spp., and an Ehrlichia canis-like agent in Rhipicephalus microplus from Southwest and South-Central China. Ticks Tick Borne Dis.; 2022; 13, 101884. [DOI: https://dx.doi.org/10.1016/j.ttbdis.2021.101884]

45. Duron, O.; Noël, V.; McCoy, K.D.; Bonazzi, M.; Sidi-Boumedine, K.; Morel, O.; Vavre, F.; Zenner, L.; Jourdain, E.; Durand, P. et al. The recent evolution of a maternally-inherited endosymbiont of ticks led to the emergence of the Q fever pathogen, Coxiella burnetii. PLoS Pathog.; 2015; 11, e1004892. [DOI: https://dx.doi.org/10.1371/journal.ppat.1004892]

46. Borsoi, A.B.P.; Bitencourth, K.; de Oliveira, S.V.; Amorim, M.; Gazêta, G.S. Human Parasitism by Amblyomma parkeri Ticks Infected with Candidatus Rickettsia paranaensis, Brazil. Emerg. Infect Dis.; 2019; 25, pp. 2339-2341. [DOI: https://dx.doi.org/10.3201/eid2512.190988]

47. Kumar, S.; Stecher, G.; Tamura, K. MEGA7: Molecular Evolutionary Genetics Analysis Version 7.0 for Bigger Datasets. Mol. Biol. Evol.; 2016; 33, pp. 1870-1874. [DOI: https://dx.doi.org/10.1093/molbev/msw054]

48. Guindon, S.; Delsuc, F.; Dufayard, J.F.; Gascuel, O. Estimating maximum likelihood phylogenies with PhyML. Methods Mol. Biol.; 2009; 537, pp. 113-137.