Skip to main content

The bacterial patterns suggesting the dynamic features of tick-associated microorganisms in hard ticks

Abstract

Background

Ticks are blood-feeding significant arthropods that can harbour various microorganisms, including pathogens that pose health risks to humans and animals. Tick-symbiont microorganisms are believed to influence tick development, but the intricate interactions between these microbes and the relationships between different tick-borne microorganisms remain largely unexplored.

Results

Based on 111 tick pool samples presenting questing and engorged statuses including 752 questing tick and 1083 engorged tick from cattle and goats, which were collected in two types of geographic landscape (semi-desert and alpine meadow). We observed significant variations in the composition of tick-borne microorganisms across different environments and blood-engorgement statuses, with a pronounced divergence in symbionts compared to environmental bacteria. Metabolic predictions revealed over 90 differential pathways for tick-borne microorganisms in distinct environments and more than 80 metabolic variations in response to varying blood engorgement statuses. Interestingly, nine pathways were identified, particularly related to chorismate synthesis and carbohydrate metabolism. Moreover, microbial network relationships within tick-borne microorganism groups were highly distinct across different environments and blood-engorgement statuses. The microbial network relationships of symbionts involve some pathogenic and environmental microorganisms. Regression modelling highlighted positive correlations between the Coxiella symbiont and related pathogens, while some environmental bacteria showed strong negative correlations with Coxiella abundance. We also identified commensal bacteria/pathogens in bacterial cooccurrence patterns. Furthermore, we tested pathogenic microorganisms of each tick sample analysis revealed that 86.36% (1601/1855) of the tick samples carried one or more pathogenic microorganisms, The total carrier rate of bacterial pathogens was 43.77% ((812/1855). Most blood samples carried at least one pathogenic microorganism. The pathogens carried by the ticks have both genus and species diversity, and Rickettsia species are the most abundant pathogens among all pathogens.

Conclusion

Our findings underscore that the bacterial pattern of ticks is dynamic and unstable, which is influenced by the environment factors and tick developmental characteristics.

Peer Review reports

Introduction

Ticks are an important vector that transmit pathogens by biting humans, domestic animals, or wildlife [1, 2]. The compositions of the tick microbiome were reported to be affected by various factors, such as tick species, life stage, hyperaemia level, and habitats and forest structures [3,4,5,6,7]. With the diversity of the microorganisms, there are usually hundreds of thousands of microorganisms in the microbial community carried by ticks [8], which leads to diverse physiological interactions between the tick hosts and their microbes. Some of these bacterial microorganisms have been proven to affect the growth and development of ticks, the blood meal behavior of ticks, and the transmission of pathogens by ticks [9, 10].

Additionally, the relationship between the microorganisms carried by ticks has not been revealed. In nature, distinct relationships, such as coexistence, mutualism, symbiosis, antagonism, competition, parasitism, and predation, exist among different microorganisms [11]. The relationship between microorganisms has been used in human body. For example, with additional intervention, the intestinal microbiota related to metabolism could be changed to achieve the metabolic health of the human host [12, 13]. However, it is a difficult and arduous task to identify the relationship between the microorganisms carried by ticks. While Amplicon sequencing of variable region V3-V4 of 16S rRNA gene technology cannot elucidate complex interaction networks between microorganisms, it can provide insights into macro-level microbial abundance. Utilizing these abundance levels, it is possible to construct an abundance correlation network to unveil trends in microbial abundance changes [14]. Understanding the changes in the abundance of key microorganisms in the tick-borne microbial community, such as symbolic bacteria and tick-borne pathogens (TBP), may help achieve similar intervention to human intestinal, i.e. manually intervene the composition of the tick-borne microbial community. Unfortunate such study is rarely reported.

In this study, We identified and characterized of bacterial microorganisms carried by ticks collected from different locations and different blood-engorgement statuses. Additionally, the bacterial co-occurrence networks of these communities were analyzed to reveal potential interaction patterns of tick-borne bacterial microorganisms.

Materials and methods

Sample collection

During the period spanning from late April to mid-May each year in 2019–2020, when adult ticks exhibited their peak activity in the Xinjiang Uygur Autonomous Region (XUAR), northwestern China. These collections encompassed both questing ticks and engorged ticks. The collection sites were situated in Wenquan County (WQ) (81°08’-81°10’E-44°95’-44°98’N, alpine meadow) and Alataw City (ALSK) (82°48’51’’E-45°04’22’’N, semi-desert). The questing ticks were gathered using the artificial trapping method [15]. The fully engorged ticks were systematically sampled from entire body of each pastured cattle and goat at different intervals, and subsequently preserved in 70% ethanol for further analysis.

With the consent of their owners, we also obtained blood samples from 25 pastured cattle in Alataw City and 79 pastured sheep in WQ County. During the period, we captured 10 great gerbils (Rhombomys opimus) from Alataw City and 29 long-tailed ground squirrels (Spermophilus undulatus) from WQ County. During the necropsy process, blood samples were collected from each rodent and stored at -80 Â°C until DNA extraction. Furthermore, we collected 368 blood samples from local herdsmen in WQ County. Detailed information about all the specimens, including their geographical location, host, the number of ticks from animal’s body, and the date of collection, was meticulously recorded (shown in Table 1).

Table 1 Information of the samples include ticks and blood sample

DNA extraction and tick molecular species identification

Before DNA extraction, ticks were carefully cleaned by 70% ethanol and phosphate buffered saline (PBS 7.4) to prevent any potential DNA contamination from environmental microorganisms. Then, they were briefly frozen in liquid nitrogen for a minute, after which they were delicately sliced into small pieces using a sterile scalpel. These tick fragments were then meticulously homogenized with a sterile micropestle in 100 µL of lysis solution containing proteinase K. Genomic DNA was extracted from each individual whole tick and blood sample using a DNeasy Blood & Tissue Kit (Qiagen, Shanghai, China) according to the manufacturer’s instructions. The concentration of extracted DNA was quantified with a spectrophotometer (NanoDrop ND-2000; Thermo Scientific).

We amplified fragments of the 16S ribosomal (r) RNA gene (455 bp) and cytochrome oxidase subunit 1 (CO1) gene (732 bp) from the DNA of each tick. The primers used for PCR were 16S-F1 (5’-TTAAATTGCTGTRGTATT-3’) and 16S-R1 (5’-CCGGTCTGAACTCASAWC-3) [16] and COX1F (5’-GGAACAATATATTTAATTTTTGG-3’) and COX1R (5’-ATCTATCCCTACTGTAAATATATG-3’) [17]. The PCR product specificity was verified by electrophoresis and Sange sequencing. The sequences were adjusted manually and compared with similar sequences retrieved from the GenBank database using BLASTn (http://www.ncbi.nlm.nih.gov/BLAST/). Phylogenetic relationships were inferred with the maximum likelihood (ML) method in MEGA 7.0 software.

16S rRNA fragment sequencing for tick carrying microbiome

Prior to sequencing, each 2–20 tick samples were constructed into a tick pool according to tick species, engorged status, host and their geographical location [20]. A total of 111 tick pools were finally constructed in this study (details in Table 1). Preparation of extracted DNA for pyrosequencing followed the protocols described by Fierer et al., 2008 [18] and Bates et al., 2010 [19]. In brief, a region of the 16S rRNA gene was amplified with the 341F (5’-CCTAYGGGRBGCASCAG-3’) and 806R (5’-GGACTACNNGGGTATCTAAT-3’) primers. Triplicate PCR products were pooled together and normalized in equimolar amounts and then loaded on an Illumina MiSeq PE250 platform for 2 × 250 bp paired-end sequencing.

Raw sequence data generated from 16S rRNA sequencing were further processed. Sequences were trimmed to delete adapters using cutadapt v3.4 (https://cutadapt.readthedocs.io/) and low-quality sequence with vsearch v2.17.1 (https://github.com/torognes/vsearch). which discards sequences whose expected errors per base is greater than 1%. Trimmed sequences were classified with Kraken2 (the second version of the Kraken taxonomic sequence classification system) [20] against the Silva database ver. 138.1 [21]. The abundance of each genus is summarized in Bracken [22].

Diversity analysis

The microbiome diversity in 4 groups: ALSKqt (questing ticks from Alataw City), ALSKet (engorged ticks from Alataw City), WQqt (questing ticks from Wenquan County) and WQet (engorged ticks from Wenquan County), were analyzed with MicrobiomeAnalyst (https://www.microbiomeanalyst.ca/), an integrated platform to comprehend microbiome data. Alpha diversity was estimated with the Shannon index, which accounts for both richness and evenness [23], and beta diversity was assessed using Bray-Curtis distance metrics and visualized using NMDS (nonmetric multidimensional scaling) to detect differences between groups. Significant differences in alpha and beta diversity among groups were tested by ANOVA (analysis of variance) and ANOSIM (analysis of similarities), respectively. LEfSe (linear discriminant analysis effect size) was then applied to distinguish genera that had different abundances among groups.

Function prediction of microbial communities

Trimmed sequences in the sequence processing section were clustered into amplicon sequence variants (ASVs) with the UNOISE 3 algorithm [24]. A 97% identity threshold, which corresponds approximately to the taxonomic level of species for bacteria [25], was used to generate ASVs with abundances for functional analysis. The ASV abundance table was delivered to Picrust2 [26] that implements a full pipeline to predict functional abundances of microbial communities based only on marker gene sequences. Pathway abundances, generated through a predictive approach in Picrust2 relying on MinPath, were selected for differential analysis between various groups. To initiate the analysis, the ALSKet group was subdivided into two subsets that had the maximum disparity in the number of each tick species (weighted by proportion). This approach aimed to simulate different tick compositions, allowing us to uncover whether ticks with varying species compositions exhibit distinct microbial community functions. Next, we conducted comparisons between the ALSK and WQ, ALSKqt and WQqt, and ALSKet and WQet groups to investigate the environmental impact on these communities. This analysis helped us to uncover variations in microbiome functions between ticks in these two states. All comparisons were carried out with the statistical analysis software STAMP v2.1.3 [27] using two-sided Welch’s t test and multiple test correction following the Bonferroni method to reduce false positive findings.

Network inference

The microbial co-occurrence networks (CoNets) help to deduce potential associations among specific genera within ticks. The genus abundance table produced by Bracken was employed to perform pairwise correlation calculations of microorganisms using the Sparse Correlations for Compositional data (SparCC) method, utilizing its default parameters [28]. The threshold value for generating regional CoNets was determined with the random matrix theory method [29], which are 0.56, 0.608, 0.308, 0.332 for group ALSKqt, ALSKet, WQqt and WQet respectively. Next, nodes and edges were generated from the adjacency matrix of pairwise correlations with the R package igraph v1.2.6 [30] in R 4.1.1 [31]. Visualization of networks and analysis of features were performed with Gephi v0.9.2 [32].

Fitting of the multiple linear regression model

Multiple linear regression was applied for two dominant symbionts, Coxiella and Francisella, to interpret the abundance change of symbionts with the abundance of their correlated genera. First, correlated genera were obtained from the SparCC correlation matrix for each symbiont, which were both positively and negatively correlated with the symbiont. Next, abundance data of symbionts and their correlated genera were split to obtain a training dataset and a test dataset at a ratio of 7:3. Then, a multiple linear regression model of abundance between the symbiont and its correlated genera was fitted with the function lm in the stats package of R, followed by optimizing the independent variable with a function step to exclude genera that contribute to a limited explanation for the abundance change of the symbiont. The optimized model was finally evaluated with the function gvlma in the gvlma package [33]. The test dataset was used to check the accuracy of the model with the function predict_lm.

Pathogen examination

All collected samples were tested nine genera of tick-borne pathogens, including Anaplasma, Babesia, Bartonella, Borrelia, Coxiella, Ehrlichia, Francisella, Rickettsia and Wolbachia by conventional PCR. The sequences of the primers used for detecting these pathogens were adopted from previous reports [34,35,36,37,38,39,40,41,42] and are listed in Additional file 1, while the primers for detecting Francisella tularensis (rpoB Fran-F/Fran-R) and Ehrichia spp. (Ehrlichia 16S rRNA EBm52F/EBm52R) were self-designed. The optimization steps of the PCR program were adapted from literature as follows: (1) Different concentrations of plasmids carrying the target gene were used as templates to optimize the reaction conditions and test the sensitivity of the primers. (2) Based on the first step, plasmids with the concentration that was the most sensitive to the reaction were selected and mixed with 60 ng DNA of tetracycline-treated ticks to form a template for mimicking the positive infection samples. Template containing 60 ng of DNA from tetracycline-treated ticks and template containing no DNA were used as control groups. Then, all these templates were used to test the specificity and sensitivity of the primers in complex sample situations.

Through experimental optimization and evaluation of the results, the optimum reaction conditions were established for the 11 pairs of primers. The volume of the reaction system for the PCR test was 30 µL in total, containing 27 µL of 1.1×T3 PCR mix, 1 µL of forward primer, 1 µL of reverse primer, and 1 µL of DNA template. The reaction program was set as follows: 1 cycle at 98 °C for 3 min, 35 cycles of 98 °C for 10 s, 54–60 °C for 10 s, 72 °C for 10 ∼ 25 s, and 1 cycle at 72 °C for 2 min. All samples were tested with all primers. The samples that tested positive for Francisella spp. and Ehrlichia spp. Ehr 521/Ehr 747 were further tested with the primers Francisella tularensis_Fran-F/Fran-R and Ehrlichia spp. EBm52F/EBm52R, respectively. Samples that tested positive with any primers were subjected to Sanger sequencing to identify their species. In addition, the species of the Francisella tularensis-positive samples were further identified by qPCR [43].

Results

Sampling and tick species identification

A total of 1855 ticks were collected, including questing ticks (n = 752) and engorged ticks (n = 1083) from livestock cattle and goats. 20 ticks (WQrt) were collected from the long-tailed ground squirrels (Table 1). EDTA-anticoagulated blood specimens were collected from parasitized livestock cattle (ALSKcb, n = 25) and goats (WQgb, n = 79). We also collected blood samples from great gerbils (ALSKrb n = 10) and long-tailed ground squirrels (WQrb n = 29), and 368 (WQhb) local herdsman blood samples were collected.

Among the 1855 ticks examined, we identified three genera and six species. Notably, 1132 ticks, comprising both questing ticks (n = 600) and engorged ticks from livestock goats (n = 512) and the long-tailed ground squirrels (n = 20), were collected in the WQ area. These ticks were conclusively identified as Dermacentor nuttalli. In our study, we observed the presence of two or three different tick species on each livestock animal. Among the 571 engorged ticks examined, we found that 49.91% were Hyalomma asiaticum (n = 285), 40.28% were Dermacentor marginatus (n = 230), and 8.07% were Hyalomma detritum (n = 46), while only a small number of ticks belonged to Dermacentor nuttalli (n = 5) and Rhipicephalus sanguineus (n = 2). Three ticks could not be identified. Of the 152 questing ticks, 147 were identified as Hyalomma asiaticum, 4 were identified as Dermacentor nuttalli, and 1 was identified as Hyalomma marginatum (n = 1) (for details see Table 1).

16S partial region sequencing

Out of the 111 tick pools comprising a total of 1855 ticks, we obtained a total of 4,770,842 quality-filtered reads. On average, each sample yielded approximately 42,981 reads, with a standard deviation of 22,548 and a range spanning from 19,265 to 65,352. The quality of these reads, indicated by an average quality score of 36, was deemed satisfactory. It is noteworthy that all libraries generated from tick samples exhibited ample sequencing depth for subsequent analysis. This is evident from the rarefaction curves plotting the number of observed OTUs, which reached a plateau when considering a sequencing depth ranging from 5000 to 10,000 sequences (Additional file 2). This plateau suggests that our samples had achieved sufficient coverage, validating their suitability for further analysis.

Microbiome profile

At the phylum level, tick pools in Alataw City carried more diverse bacteria. that those in Wenquan County The shared bacteria between ticks from the two locations belong to 9 phyla: Acidobacteriota, Actinobacteriota, Bacteroidota, Deinococcota, Firmicutes, Gemmatimonadota, Proteobacteria, Fusobacteriota and Verrucomicrobiota, and there were 11 exclusive phyla in ALSK ticks but only one exclusive phylum in WQ ticks (Additional file 3). Regarding abundance differences among groups, ticks from the two environments showed bacterial differences in 14 phyla out of all 21 phyla (Additional file 4). Specifically, questing ticks from ALSK carried more Acidobacteriota (W = 325, P = 0.0117), Actinobacteriota (W = 447, P < 0.0001), Bacteroidota (W = 426.5, P = 0.0005), Cyanobacteria (W = 297.5, P = 0.0315), Gemmatimonadota (W = 297.5, P = 0.0315) and Patescibacteria (W = 420, P < 0.0001) than WQ ticks. In addition, engorged ones from ALSK carried more Actinobacteriota (W = 807.5, P < 0.0001), Bacteroidota (W = 775, P < 0.0001), Fusobacteriota (W = 516, P = 0.0367) and Patescibacteria (W = 577.5, P = 0.0016) but less Deinococcota (W = 188, P = 0.0004) and Proteobacteria (W = 269, P = 0.0215). Regarding questing ticks and engorged ticks, engorged tick pools carried more diverse bacteria than questing tick pools. There were 13 shared bacterial phyla and no exclusive phyla in questing ticks but 8 exclusive phyla in engorged ticks (Additional file 5). Furthermore, questing ticks from ALSK carried more Acidobacteriota (W = 363, P = 0.0291) and Deinococcota (W = 493, P = 0.0003) but less Fusobacteriota (W = 217.5, P = 0.0265) and Verrucomicrobiota (W = 232.5, P = 0.0494) than engorged ticks. The questing ticks from WQ carried more Deinococcota (W = 525, P = 0.0071) but less Bacteroidota (W = 232.5, P = 0.0059) and Firmicutes (W = 209, P = 0.0067) than engorged tick pools.

At the genus level, the abundance table consisted of 495 genera. Similar to the phylum level, tick pools in Alataw City carried more diverse bacteria than those in WQ County, and engorged tick pools carried more diverse bacteria than questing tick pools (Additional file 6–7). For the 4 groups, there were 203.5 ± 147.63 (mean ± standard deviation) genera for each. Thirty-nine (7.88%) were shared in all groups (Fig. 1A). The shared genera occupied 86.21% of the total abundance. The engorged tick pools in Alataw City had the most exclusive genera (238, 78.81% in all exclusive genera), and 17 among these (7.15%) were from exclusive phyla of this group. While the engorged tick pools in WQ County had the least exclusive genera (3, 0.99% in all exclusive genera), one of these belongs to the only exclusive phylum of this group. In general, the genera with the top 10 relative abundances were Coxiella (22.99%), Rickettsia (17.90%), Acinetobacter (12.23%), Francisella (9.65%), Anaplasma (5.94%), Psychrobacter (5.25%), Thermus (3.19%), Pseudomonas (2.92%), Anoxybacillus (2.39%) and Rickettsiales bacterium Ac37b (1.74%) in all 111 samples. These high abundance genera are included in 39 shared genera except the genus Anaplasma, and 8 in 10 belong to the phylum Proteobacteria, which showed that shared genera formed a major component of the tick-carrying bacterial microbiome. The diversity of the bacteria within individual samples (alpha diversity) shows richer in engorged tick pools in ALSK (P = 5.4e-6) (Fig. 1B). However, the alpha diversity of bacteria in WQ ticks displayed no difference between questing and engorged tick pools (P = 0.37). The beta diversity between groups was assessed using the Bray–Curtis distance metrics and visualized using NMDS (Fig. 1C). In beta diversity observations, the genera of the 4 groups were significantly distinct from each other (ANOSIM, p < 0.01). LEfSe determined genera that had different abundances between groups (Fig. 1D-E). It showed 7 differential genera between questing tick-carrying bacteria from two places and 22 between engorged tick-carrying bacteria (Fig. 1D). Among these, the top differential genera according to LDA score were intracellular bacteria. However, the most differential genera were environmental bacteria. Two common differential genera that occurred in both places were Coxiella and Thermus, which are more abundant in questing ticks. For different engorged statuses, there were 16 and 6 differential genera between questing ticks and engorged ticks from ALSK and WQ (Fig. 1E). Similarly, intracellular bacteria were the top differential genera, and the most differential genera were environmental bacteria.

Fig. 1
figure 1

Profiles of the bacterial microbiome from ticks at the genus level. (A) Venn diagram of bacterial genera for 4 groups of ticks: questing ticks from ALSK (ALSKqt), engorged ticks from ALSK (ALSKet), questing ticks from WQ (WQqt) and engorged ticks from WQ (WQet); (B) Alpha diversity of bacterial genera for the 4 groups of ticks above; (C) Beta diversity of bacterial genera for the 4 groups of ticks above; (D) Bacterial genera with different abundances between environments; (E) Bacterial genera with different abundances between engorgement statuses

The predicted function of microbial communities

Four groups were used to detect differences in the functions of microbial communities to inspect the relation of environment and parasite statuses on observed microbial communities. First, 40 samples of the ALSKet group were divided into two subsets, both of which contain 20 samples but have different tick compositions. In the two subsets, one consisted of 0.69% D. nuttalli, 61.51% Hy. asiaticum, 27.84% D. marginatum, 9.62% Hy. detritum and 0.34% R. sanguineus (n = 291), and the other subset consisted of 1.08% D. nuttalli, 38.27% Hy. asiaticum, 53.79% D. marginatum, 6.5% Hy. detritum and 0.36% R. sanguineus (n = 277). Statistical analysis showed no significant change in the functional abundance of microbial communities (Additional file 8), suggesting that tick species have a limited influence on the formation of microbial functions in ticks.

For ALSK and WQ tick pools, there were 222 differential function pathways in bacteria (Additional file 9). When engorged statuses were considered, there were 97 detected in questing ticks and 237 detected in engorged ticks (Additional file 10–11). This means that there are many metabolic differences among tick-carrying bacteria in different environments. By inspecting microbial function abundance between ALSKqt and ALSKet, WQqt and WQet, 106 and 83 differential function pathways were discovered (corrected Ps < 0.05, Additional file 12–13), respectively. Among these pathways, there were 9 shared differential function pathways (Fig. 2), including chorismate biosynthesis I, chorismate biosynthesis from 3-dehydroquinate, the superpathway of β-D-glucuronide and D-glucuronate degradation, glycogen degradation I (bacterial), GDP-mannose biosynthesis, D-fructuronate degradation, the superpathway of aromatic amino acid biosynthesis and tRNA processing, which were all less abundant in engorged ticks, and one pathway of biotin biosynthesis, which was more abundant in engorged ticks.

Fig. 2
figure 2

The 9 shared differential metabolic pathways between questing ticks and engorged ticks from two environments. Bar plot with extended error bar. Red and blue objects represent groups of questing ticks and engorged ticks, respectively

Difference between microbial co-occurrence networks (CoNets) of questing and engorged ticks

Microbial co-occurrence networks representing different groups were constructed, and the degree parameters of CoNets are displayed in Fig. 3. The vertices and average degrees of the CoNets of tick-carrying bacteria on questing and engorged ticks of ALSK are less than those of WQ. Similarly, the vertices and average degrees of the CoNets of tick-carrying bacteria on engorged ticks are less than those on questing ticks. These results suggest that the CoNets of tick-carrying bacteria were affected by the environments where ticks grow and in the engorgement status of the ticks.

Fig. 3
figure 3

Comparison between microbial co-occurrence networks for tick-carrying bacteria at different levels. (A) microbial co-occurrence networks for ticks from ALSK and WQ, left for group questing ticks, right for engorged ticks; (B) microbial co-occurrence networks for questing and engorged ticks, left for group ALSK, right for group WQ; (C) dominated genera which account for 50% total degree of net of each group, orange lines indicated shared key genera between two groups. Each ball represents a genus, and the size of the ball indicates the degree of this genus. Red and green balls indicate dominant genera of different groups, which account for 50% of the total degree of net. Orange balls indicate shared dominant genera for both groups. Links represent co-occurrence relationships between genera, and their colours show which group displayed this relationship. Values in the dashed box indicate the average degree of the color corresponding group in the cooccurrence network

All groups were composed of different taxonomic profiles and were dominated by the phylum Proteobacteria (Additional file 14). For the abundance of bacteria in CoNets, ticks from two places showed bacterial differences in 5 phyla out of all 10 phyla. Specifically, questing ticks from ALSK carried more Actinobacteriota (W = 417, P < 0.0001), Bacteroidota (W = 393, P = 0.0004), Firmicutes (W = 329, P = 0.0391), and Patescibacteria (W = 391, P < 0.0001) and less Proteobacteria (W = 134, P = 0.0176) than WQ ticks. In addition, engorged ones from ALSK carried more Actinobacteriota (W = 749.5, P < 0.0001), Bacteroidota (W = 709, P < 0.0001), and Patescibacteria (W = 540, P = 0.0018) and less Deinococcota (W = 172, P = 0.0005) and Proteobacteria (W = 112, P < 0.0001). Additionally, questing ticks from ALSK carried more Deinococcota (W = 458, P = 0.0002) and Proteobacteria (W = 418, P = 0.0028) than engorged ticks. Questing ticks from WQ carried more Deinococcota (W = 487, P = 0.0079) but less Bacteroidota (W = 463, P = 0.0224) and Firmicutes (W = 188, P = 0.0059) than engorged ticks. For dominant genera that account for 50% of the total degree of net, few of the same genera were detected between different places or different engorgement statuses. These results showed that the tick-carrying bacteria included in CoNets between the two environments had significant differences in taxonomy and abundance. Meanwhile, in same environment, no pattern for bacterial differences between tick-questing and engorged ticks was detected.

Pathogen-related genera correlated to symbionts Coxiella and Francisella

After extracting local correlations for pathogen-related genera and 2 high abundance symbionts (Coxiella and Francisella) from the SparCC correlation matrix, a correlation network was summarized, as shown in Fig. 4. There were 63 connected to Coxiella in a total of 85 genera and 28 connected to Francisella. Among these, 44 environmental genera, 9 pathogen-related genera, 6 gut-related genera, 2 symbiont genera and 2 undetermined genera were connected to Coxiella. In addition, 23 environmental genera, 2 pathogen-related genera, 1 gut-related genus and 2 symbiont genera were connected to Francisella. The remaining 21 genera connected to Coxiella were mediated by one genus, but Francisella still had 10 genera that were not connected to Coxiella. Pathogen-related genera that were positively correlated with Coxiella were Anaplasma, Mycoplasma, Rickettsia, Roseomanos, Spiroplasma and Ehrlichia. Arsenophonus, Moraxella, Samonella, and Trueperella were negatively correlated with Coxiella. In addition, Ehrlichia and Mycoplasma were directly correlated with Francisella, both positively, and 5 indirectly correlated with Coxiella. These results revealed that the main correlations of symbionts were environmental genera, and some important pathogen-related genera correlated to symbionts where these genera showed different correlations with each symbiont. Additionally, no stable pattern of relationships between the two symbionts exists in different environments or engorged statuses.

Fig. 4
figure 4

Genera correlated with endosymbionts Coxiella and Francisella. Each ball represents a genus, and the size of the ball indicates the degree of this genus. Orange, red, purple, blue and green balls represent endosymbionts, pathogen-related genera, environmental bacteria, gut bacteria and unknown genera, respectively. Names of endosymbionts and pathogen-related genera are displayed

Multiple linear regression model (MLRM) revealed a quantitative relation between Coxiella and its correlated genera

For the two symbionts, only the dataset of Coxiella passed the evaluation of MLRM in gvlma, which suggested that the abundance change between Francisella and its correlated genera did not satisfy a linear relation. The MLRM formula for the abundance of Coxiella comprised 7 intracellular parasites, 6 gut-related genera and 28 environmental genera (Fig. 5). The abundance of these 41 genera interpreted a 99.97% change in Coxiella abundance. For the test dataset, the R squared between the prediction of the model and observation was 0.7280 after the negative predicted value was corrected to 0. The top 3 coefficients with the most weights, which belong to the genera Phenylobacterium, Moraxellla and Mycoplasma, were all negative. The top 3 positive coefficients with the most weights were possessed by the genera Herminiimonas, UCG-009 and Cloacibacterium. The weight proportions of the 6 genera in the formula were 21.96%, 17.92%, 7.43%, 7.18%, 6.34% and 3.46%. These results revealed a quantitative relation between Coxiella and its correlated genera, and environmental genera were the most important factors influencing the abundance of symbionts of Coxiella.

Fig. 5
figure 5

The 41 genera and their coefficients that have abundance correlations with Coxiella revealed by a multiple linear regression model. Negative and positive coefficients indicate negative and positive correlations with Coxiella, respectively. The absolute value of the coefficient indicates the weight of the correlation with Coxiella

Pathogens and symbionts in ticks

A total of 8 microorganisms with different carrier rates were detected in samples from 10 sources, including pathogens and tick symbionts: Ehrilichia spp. (0.42%, 10/2336), Bartonella spp. (0.77%, 18/2336), Babesia spp. (5.77%, 135/2336), Borrelia spp. (7.49%, 175/2336), Francisella spp. (19.34%, 452/2336), Anaplasma spp. (21.48%, 502/2336), Rickettsia spp. (30.52%, 713/2336), and Coxiella spp. (52.48%, 1226/2336). In addition, 86.36% of the tick samples carried one or more of the detected bacterial microorganisms, and the carrying rates from low to high were 20% (4/20) for WQrt, 75.66% (432/571) for ALSKet from livestock cattle, 87.33% (524/600) for WQqt, 95.70% (490/512) for WQet from goats and 100% (152/152) for ALSKqt. Moreover, 43.77% (812/1855) of tick samples carried bacterial pathogens, and nearly 18% of blood samples carried one or two detected pathogens. More specifically, the bacterial pathogen carrying rates of the blood samples from low to high were 0% (0/29) for the WQrb, 0% (0/10) for the ALSKrb, 1.6% (6/368) for the WQhb, 46% (13/28) for the ALSKcb and 94% (74/79) for the WQgb (for details see Table 2).

Table 2 Detection the statistics of individual microbes within each sample type

Ehrilichia spp.

Sanger sequencing analysis showed that these Ehrilichia sequences were the same and had 100% homology with 16S rRNA sequence of Ehrilichia (JX402605.1) found in Hy. asiaticum ticks in Xinjiang Province (northwestern China), the Ehrilichia 16S rRNA sequence KX987325.1 carried by Boophilus microplus ticks in Hubei Province (central China) and the Ehrilichia 16S rRNA sequence KY046298.1 carried by Rhipicephalus microplus ticks in Malaysian [44]. By aligning the sequences with the Basic Local Alignment Search Tool (BLAST), it was found that this Ehrilichia sequence had 100% homology with MT875368.1 (Candidatus Ehrlichia Hainanensis), KJ513197.1 (Ehrlichia canis), MZ733621.1 (Candidatus Ehrlichia pampeana) and MT738235.1 (Ehrlichia ruminantium). These results suggest that this species may be a pathogen that is prevalent in humans, livestock, and rodents.

Bartonella spp.

Bartonella spp. were detected only in tick samples collected from Wenquan County. The results of Sanger sequencing showed that Bartonella spp. had 100% homology with Bartonella henselae. Bartonella henselae is a zoonotic pathogen causing neurological disease [45]. However, this pathogen was not detected in the human blood and goat blood samples collected from Wenquan County. Since the felines and canines in Wenquan County were not included in the scope of this investigation, the situation was still unknown.

Babesia spp.

Babesia spp. was detected in a total of 136 cases in ticks, livestock and human blood samples from Wenquan County and Alataw City. By sequencing and sequence alignment, it was found that the sequence of the Babesia spp. in questing ticks and parasitic ticks had 100% homology with LC553515.1, which is a Babesia strain closely related to the zoonotic Babesia gibsoni [46].

Borrelia spp.

Borrelia spp. were detected in a total of 175 ticks. The results of Sanger sequencing showed that all Borrelia spp. had the same sequence, which was 100% homology with the NR_170496.1 sequence of Borrelia maritima. B. maritima was a new species of the Borrelia burgdorferi sensu lato complex [47, 48].

Francisella spp.

Francisella spp. was detected in a total of 452 samples (5 cases in human blood samples and others in tick samples). The results of Sanger sequencing showed that 4 cases of Alataw City engorged ticks from livestock cattle carried Francisella tularensis subsp. (CP009653.1), and the others carried Francisella-like symbionts (KX852466), which is a tick symbiont that is transmitted vertically through the maternal line. These results suggest that Francisella-like symbionts are stable in the population of ticks in Alataw City.

Anaplasma spp.

Anaplasma spp. was detected in a total of 502 samples. Two species were identified by Sanger sequencing and sequence alignment, including Anaplasma ovis (NZ_CP015994.1), which was found in 416 cases, and Anaplasma phagocytophilum (NC_021880.1), which was found in 86 cases. Anaplasma ovis is the Anaplasma (Rickettsiales: Anaplasmataceae) of susceptible goat. The positive rate of goat blood samples collected from Wenquan County was 93.67% (74/79), and the positive rate of engorged ticks from goat, was 63.08% (323/512), suggesting that the goats that are raised in Wenquan County should be further tested to prevent the spread of pathogens. Moreover, Anaplasma phagocytophilum, an intracellular parasite that can cause human granulocytic anaplasmosis, was detected in tick samples from both places [49], suggesting that the prevention and control of A. phagocytophilum infection in both Wenquan County and Alataw City need to be reinforced.

Rickettsia spp.

Rickettsia spp. was detected in a total of 713 samples. Among those cases, only 1 case was from cattle blood samples in Alataw City, and the rest were from tick samples. Six species were identified by next-generation sequencing, and the sequence alignment results indicated that 5 of the six species were zoonotic pathogens, including Rickettsia raoultii (229/713) [50], Rickettsia sibirica (53/713) [50], Rickettsia sp. Hme_HirooL009 (4/713) (LC544134.1), uncultured Rickettsia sp. (KM587631.1; KM587632.1; KM587633.1; MT434895.1; KX591658.1) (111/713) [51] and Rickettsia sp. BJ-90 (2/713) [52]. Rickettsia sp. BJ-90 is a pathogen identified as belonging to the Rickettsia sibirica clade by evolutionary analysis. Taken together, the carrying rate of the Rickettsia pathogen in ticks was 55.96%. The Rickettsia species of the remaining 314 samples could not be determined, accounting for 44% of the total sequences.

Coxiella spp.

A total of 1226 Coxiella spp. infection cases were detected, and all these cases were ticks collected from the two counties. All Coxiella spp. had the same sequence, which was homologous to the sequence of uncultured Coxiella sp. clone Dx-56 and that of JX432012.1, indicating that they were tick Coxiella-like symbionts. The carrying rates of Coxiella spp. in ticks collected were 80.83% (485/600) for WQqt, 76.75% (393/512) for WQet from goats, 10% (2/20) for WQrt, 42.91% (245/570) for ALSKet from livestock cattle, and 66.45% (101/152) for ALSKqt. These results suggest that the Coxiella-like symbionts are widely distributed and stable in the tick populations in the two regions. According to the number of detected microorganisms, most of the pathogen-carrying blood samples only carried one detected microorganism, and only one case was detected in the Wenquan goat blood sample carrying both Ehrilichia spp. and Anaplasma ovis. However, tick samples often carry up to 6 pathogens (Fig. 6).

Fig. 6
figure 6

The proportion of the different coexist number of microorganisms in different sample types. Different color column represents the different coexist number of microorganisms in that sample type. WQqt: questing ticks from WQ; WQet: engorged ticks from WQ; WQgb: goat blood from WQ; WQst: squirrel’s ticks from WQ; WQsb: squirrel blood from WQ; WQhb: human blood from WQ; ALSKqt: questing ticks from ALSK; ALSKet: engorged ticks from ALSK; ALSKcb: cattle blood from ALSK; ALSKgb: gerbil blood from ALSK

Discussion

The microbiome of a variety of ticks has been described and characterized in previous reports [53, 54], and some of the components, such as Francisella/Coxiella symbionts, have been shown to be physiologically significant for ticks [10]. This study showed that the natural patterns of tick microbial communities is dynamic and unstable, the living environment and parasitic state have a significant impact on the microbiota composition and bacterial cooccurrence pattern of ticks. The microbial cooccurrence network takes environmental microorganisms as the central node, and no consistency is found between different groups(Special attention: In order to minimize the possibility of detecting contamination, tick samples are carefully cleaned with 70% ethanol and phosphate-buffered saline (PBS PH7.4) after collection to prevent any potential contamination of DNA by environmental microorganisms.). In addition, and explored a possible approach to control ticks or pathogens by manipulating the tick microbial community.

Metabolic pathways of tick-carrying bacteria

There is an intimate relationship between insect vectors and the microorganisms they carry [55]. The symbiotic bacteria of ticks have been shown to provide chemical molecules, such as chorismite [56]\5-HT [57]\B vitamins [9]\cbb3-cyt b oxidase [58]\ATP [58], to help ticks develop. In this study, the results indicated that the microorganisms carried by different tick species collected from same area has no preference. There were nine pathways involved in the parasitic state change, and three of them, the chorismate biosynthesis I, chorismate biosynthesis from 3-dehydroquinic acid, and aromatic amino acid synthesis pathways, have been shown to be associated with promoting blood meal appetite [57]. The results also suggest that changes in parasitic status did not just cause appetite changes. However, the biological significance of the pathway changes needs further study.

No stable pattern indicated by co-occurrence networks

Co-occurrence networks are widely used to study relationships between microbial communities [59]. Our results showed that the co-occurrence network of tick-carrying microorganisms were significantly affected by the grow location of ticks, their blood saturation level, and developmental status, which indicates that the tick-carrying microorganisms could not form a stable community in the host. Studies of laboratory-raised ticks have directly shown that ticks can acquire microorganisms from the environment and blood meal [8]. It can be reasoned that the acquisition of environmental microorganisms and the changes in the physiological environment of ticks before and after their blood meal could be the main reasons for the changes in the diversity and abundance of tick-borne microorganisms. In addition, the existence of mutualism, symbiosis, competition, and antagonism among microorganisms may result the differences in microbial composition in the ticks.

Coxiella symbionts

It has been reported that Coxiella symbionts are involved in the life activities of ticks in many ways [57, 59]. We found that the blood saturation level of ticks significantly affected the abundance of Coxiella symbionts and the abundances of 3 pathways related to blood meal behavior. We also found that Coxiella symbionts were associated with 10 pathogens. These results suggest that ticks carrying a high abundance of Coxiella may carry a higher abundance of related pathogens, although we don’t know if Coxiella influences pathogens, or vice versa. But a possible point is that inhibiting Coxiella abundance may control the tick borne pathogens. On the one hand, ticks with low Coxiella abundance have lower willingness to find hosts for blood meal. On the other hand, low Coxiella abundance may be not sufficient for related pathogen enrichment. Interestingly, Phenylobacterium was the negatively correlated bacterium with the Coxiella abundance. It has been reported that Phenylobacterium can be isolated from environmental samples and can only be enriched in media containing benzene ring drugs [60], which means that an environment with a high abundance of Phenylobacterium may not be beneficial to the enrichment of tick-carrying Coxiella. Thus, high concentrations of benzene-containing compounds may be able to control the abundance of tick-carrying Coxiella to reduce the blood-sucking appetite of ticks and control tick population.

Carrying and infection of different tick-borne pathogens

The process by which ticks transmit pathogens is affected by many factors [61, 62]. A high tick-borne pathogen carrying rate increases the probability of transmitting pathogens by ticks, but it will be affected by many factors, such as the number of pathogens carried, the suitability of hosts to the pathogens, the immune clearance ability of the hosts, and the immune evasion ability of the pathogens. Our research showed that Rickettsia had a low infection rate in animals and humans, which may be related to the content of the pathogen and the immunity of the host to the pathogen [63, 64]. In addition, Anaplasma ovis had a high infection rate in suitable host goat while the Anaplasma phagocytophilum could not be found in hosts. Thus, the suitable host, high tick carrying rate and high pathogen abundance cause a large-scale infection of Anaplasma ovis in goat while the zero infections occurring in non-suitable hosts could be related to immunity [65]. The infection rate of Babesia in livestock is higher than that in ticks in the Alataw City area, which may be due to the long existence of Babesia in host [66]. Moreover, F. tularensis is a highly virulent zoonotic pathogens [67], and inhalation of 10 CFU of F. tularensis can cause severe disease. Therefore, prevention and control of Francisella tularensis subsp transmission are still critical.

Limitation

Although we found a negative or positive correlation between the abundance of the Coxiella symbiont and some environmental bacteria or related pathogens by the abundance linear regression model, this finding needs to be further confirmed by more tests in the future.

Conclusions

This study shows that the species and quantities of bacteria carried by ticks are diverse as they are affected by the environment and parasitic states. The intracellular symbiotic bacteria, Coxiella/Francisella, shows the most variations in the abundance of the species and quantities. Meanwhile, Coxiella abundance is linearly correlated with the abundance of its associated microbes. Moreover, Phenylobacterium had a strong negative correlation with the symbiont Coxiella. In addition, tick-borne microbes endow ticks with several metabolic functions and the abundance of metabolic pathways of bacteria is affected by the environment and the ticks’ blood saturation levels. Thus, no stable coexistence relationship was determined in the tick-carrying microbiota coexistence network. The carrying rate of tick-borne pathogens was very high, and there were significant differences among different tick borne pathogens. Taken together, the microbes in the environment and the egg-borne symbionts play an important role in the composition of tick microbial communities. This study provides directions for further understanding and utilization of tick microbiomes, as well as control strategies for ticks and tick-borne pathogens.

Data availability

Sequence data that support the findings of this study have been deposited to to SRA with the BioProject ID PRJNA1048745(https://www.ncbi.nlm.nih.gov/bioproject/PRJNA1048745). The datasets supporting the conclusions of this article are included within the article and its additional files.

Abbreviations

WQ:

Wenquan county

ALSK:

Alataw City

NMDS:

Nonmetric multidimensional scaling

ANOVA:

Analysis of variance

ANOSIM:

Analysis of similarities

LEfSe:

Linear discriminant analysis effect size

ASV:

Amplicon sequence variant

CoNets:

Co-occurrence networks

SparCC:

Sparse Correlations for Compositional data

MLRM:

Multiple linear regression model

References

  1. de la Fuente J, Estrada-Pena A, Venzal JM, et al. Overview: ticks as vectors of pathogens that cause disease in humans and animals. Front Biosci. 2008;13:6938–46.

    Article  PubMed  Google Scholar 

  2. Narasimhan S, Swei A, Abouneameh S, et al. Grappling with the tick microbiome. Trends Parasitol. 2021;37(8):722–33.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  3. Chicana B, Couper LI, Kwan JY et al. Comparative microbiome profiles of Sympatric Tick species from the Far-Western United States. Insects. 2019;10(10).

  4. Moreno CX, Moy F, Daniels TJ, et al. Molecular analysis of microbial communities identified in different developmental stages of Ixodes scapularis ticks from Westchester and Dutchess Counties, New York. Environ Microbiol. 2006;8(5):761–72.

    Article  CAS  PubMed  Google Scholar 

  5. Swei A, Kwan JY. Tick microbiome and pathogen acquisition altered by host blood meal. ISME J. 2017;11(3):813–6.

    Article  PubMed  Google Scholar 

  6. Wielinga PR, Gaasenbeek C, Fonville M, et al. Longitudinal analysis of tick densities and Borrelia, Anaplasma, and Ehrlichia infections of Ixodes ricinus ticks in different habitat areas in the Netherlands. Appl Environ Microbiol. 2006;72(12):7594–601.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  7. van Overbeek L, Gassner F, van der Plas CL, et al. Diversity of Ixodes ricinus tick-associated bacterial communities from different forests. FEMS Microbiol Ecol. 2008;66(1):72–84.

    Article  PubMed  Google Scholar 

  8. Li SS, Zhang XY, Zhou XJ et al. Bacterial microbiota analysis demonstrates that ticks can acquire bacteria from habitat and host blood meal. Exp Appl Acarol. 2022 May 9.

  9. Duron O, Morel O, Noel V, et al. Tick-Bacteria Mutualism depends on B vitamin synthesis pathways. Curr Biol. 2018;28(12):1896–902. e5.

    Article  CAS  PubMed  Google Scholar 

  10. Ahantarig A, Trinachartvanit W, Baimai V, et al. Hard ticks and their bacterial endosymbionts (or would be pathogens). Folia Microbiol (Praha). 2013;58(5):419–28.

    Article  CAS  PubMed  Google Scholar 

  11. Ma B, Wang Y, Ye S, et al. Earth microbial co-occurrence network reveals interconnection pattern across microbiomes. Microbiome. 2020;8(1):82.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  12. Fan Y, Pedersen O. Gut microbiota in human metabolic health and disease. Nat Rev Microbiol. 2021;19(1):55–71.

    Article  CAS  PubMed  Google Scholar 

  13. Wargo JA. Modulating gut microbes. Science. 2020;369(6509):1302–3.

    Article  CAS  PubMed  Google Scholar 

  14. Rottjers L, Faust K. From hairballs to hypotheses-biological insights from microbial networks. FEMS Microbiol Rev. 2018;42(6):761–80.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  15. Tsunoda T, Kasuga S, Amano H. Estimation of the density of nymphs of the bush tick, Haemaphysalis longicornis (Acari: Ixodidae), by the catch effort method. Exp Appl Acarol. 2004;32(1–2):69–78.

    Article  CAS  PubMed  Google Scholar 

  16. Lv J, Wu S, Zhang Y, et al. Development of a DNA barcoding system for the Ixodida (Acari: Ixodida). Mitochondrial DNA. 2014;25(2):142–9.

    Article  CAS  PubMed  Google Scholar 

  17. Chitimia L, Lin RQ, Cosoroaba I, et al. Genetic characterization of ticks from southwestern Romania by sequences of mitochondrial cox1 and nad5 genes. Exp Appl Acarol. 2010;52(3):305–11.

    Article  PubMed  Google Scholar 

  18. Fierer N, Hamady M, Lauber CL, et al. The influence of sex, handedness, and washing on the diversity of hand surface bacteria. Proc Natl Acad Sci U S A. 2008;105(46):17994–9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  19. Bates ST, Berg-Lyons D, Caporaso JG, et al. Examining the global distribution of dominant archaeal populations in soil. ISME J. 2011;5(5):908–17.

    Article  CAS  PubMed  Google Scholar 

  20. Wood DE, Lu J, Langmead B. Improved metagenomic analysis with Kraken 2. Genome Biol. 2019;20(1):257.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  21. Quast C, Pruesse E, Yilmaz P, et al. The SILVA ribosomal RNA gene database project: improved data processing and web-based tools. Nucleic Acids Res. 2013;41(Database issue):D590–6.

    CAS  PubMed  Google Scholar 

  22. Lu J, Breitwieser FP, Thielen P et al. Bracken: estimating species abundance in metagenomics data. Peerj Comput Sci. 2017 Jan 2.

  23. Shannon CE, editor. The mathematical theory of communication1950.

  24. Edgar RC. UNOISE2: improved error-correction for Illumina 16S and ITS amplicon sequencing. bioRxiv. 2016:081257.

  25. Konstantinidis KT, Tiedje JM. Prokaryotic taxonomy and phylogeny in the genomic era: advancements and challenges ahead. Curr Opin Microbiol. 2007;10(5):504–9.

    Article  CAS  PubMed  Google Scholar 

  26. Douglas GM, Maffei VJ, Zaneveld JR, et al. PICRUSt2 for prediction of metagenome functions. Nat Biotechnol. 2020;38(6):685–8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  27. Parks DH, Tyson GW, Hugenholtz P, et al. STAMP: statistical analysis of taxonomic and functional profiles. Bioinformatics. 2014;30(21):3123–4.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  28. Friedman J, Alm EJ. Inferring correlation networks from genomic survey data. PLoS Comput Biol. 2012;8(9):e1002687.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  29. Luo F, Zhong J, Yang Y et al. Application of random matrix theory to biological networks. Phys Lett A. 2006 2006/09/25/;357(6):420–3.

  30. Csardi G, Nepusz T. The igraph software package for complex network research. Interjournal Complex Syst. 2006;1695.

  31. R Core Team. R: a language and environment for statistical computing. Foundation for Statistical Computing. Vienna, Austria.;; 2015.

  32. Bastian M, Heymann S, Jacomy M, Gephi. An Open Source Software for Exploring and Manipulating Networks. 2009. (2009).

  33. Pena EA, Slate EH. Global validation of linear model assumptions. J Am Stat Assoc. 2006;101(473):341.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  34. Pancholi P, Kolbert CP, Mitchell PD, et al. Ixodes dammini as a potential vector of human granulocytic ehrlichiosis. J Infect Dis. 1995;172(4):1007–12.

    Article  CAS  PubMed  Google Scholar 

  35. Zhan L, Chu CY, Zuo SQ, et al. Anaplasma phagocytophilum and Borrelia burgdorferi in rabbits from southeastern China. Vet Parasitol. 2009;162(3–4):354–6.

    Article  PubMed  Google Scholar 

  36. Masuzawa T, Sawaki K, Nagaoka H, et al. Identification of rickettsiae isolated in Japan as Coxiella burnetii by 16S rRNA sequencing. Int J Syst Bacteriol. 1997;47(3):883–4.

    Article  CAS  PubMed  Google Scholar 

  37. Rolain JM, Franc M, Davoust B, et al. Molecular detection of Bartonella quintana, B. Koehlerae, B. Henselae, B. clarridgeiae, Rickettsia felis, and Wolbachia pipientis in cat fleas, France. Emerg Infect Dis. 2003;9(3):338–42.

    Article  CAS  PubMed  Google Scholar 

  38. 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:1449–55.

    Article  CAS  PubMed  Google Scholar 

  39. Scoles GA. Phylogenetic analysis of the Francisella-like endosymbionts of Dermacentor ticks. J Med Entomol. 2004;41(3):277–86.

    Article  CAS  PubMed  Google Scholar 

  40. Werren JH, Zhang W, Guo LR. Evolution and phylogeny of Wolbachia: reproductive parasites of arthropods. Proc Biol Sci. 1995;261(1360):55–63.

    Article  CAS  PubMed  Google Scholar 

  41. Marconi RT, Garon CF. Development of polymerase chain reaction primer sets for diagnosis of Lyme disease and for species-specific identification of Lyme disease isolates by 16S rRNA signature nucleotide analysis. J Clin Microbiol. 1992;30(11):2830–4.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  42. Casati S, Sager H, Gern L, et al. Presence of potentially pathogenic Babesia sp. for human in Ixodes ricinus in Switzerland. Ann Agric Environ Med. 2006;13(1):65–70.

    CAS  PubMed  Google Scholar 

  43. Larson MA, Sayood K, Bartling AM et al. Differentiation of Francisella tularensis subspecies and subtypes. J Clin Microbiol. 2020;58(4).

  44. Lu M, Tian JH, Yu B, et al. Extensive diversity of rickettsiales bacteria in ticks from Wuhan, China. Ticks Tick Borne Dis. 2017;8(4):574–80.

    Article  PubMed  Google Scholar 

  45. Alsmark CM, Frank AC, Karlberg EO, et al. The louse-borne human pathogen Bartonella quintana is a genomic derivative of the zoonotic agent Bartonella henselae. Proc Natl Acad Sci U S A. 2004;101(26):9716–21.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  46. Chatanga E, Kainga H, Maganga E, et al. Molecular identification and genetic characterization of tick-borne pathogens in sheep and goats at two farms in the central and southern regions of Malawi. Ticks Tick Borne Dis. 2021;12(2):101629.

    Article  PubMed  Google Scholar 

  47. Margos G, Fedorova N, Becker NS, et al. Borrelia maritima sp. nov., a novel species of the Borrelia burgdorferi sensu lato complex, occupying a basal position to north American species. Int J Syst Evol Microbiol. 2020;70(2):849–56.

    Article  CAS  PubMed  Google Scholar 

  48. Pritt BS, Respicio-Kingry LB, Sloan LM, et al. Borrelia mayonii sp. nov., a member of the Borrelia burgdorferi sensu lato complex, detected in patients and ticks in the upper midwestern United States. Int J Syst Evol Microbiol. 2016;66(11):4878–80.

    Article  CAS  PubMed  Google Scholar 

  49. Stuen S, Granquist EG, Silaghi C. Anaplasma phagocytophilum–a widespread multi-host pathogen with highly adaptive strategies. Front Cell Infect Microbiol. 2013;3:31.

    Article  PubMed  PubMed Central  Google Scholar 

  50. Igolkina Y, Krasnova E, Rar V, et al. Detection of causative agents of tick-borne rickettsioses in Western Siberia, Russia: identification of Rickettsia raoultii and Rickettsia sibirica DNA in clinical samples. Clin Microbiol Infect. 2018;24(2):199. e9-199 e12.

    Article  Google Scholar 

  51. Ionita M, Silaghi C, Mitrea IL, et al. Molecular detection of Rickettsia conorii and other zoonotic spotted fever group rickettsiae in ticks, Romania. Ticks Tick Borne Dis. 2016;7(1):150–3.

    Article  PubMed  Google Scholar 

  52. Fournier PE, Dumler JS, Greub G, et al. Gene sequence-based criteria for identification of new rickettsia isolates and description of Rickettsia heilongjiangensis sp. nov. J Clin Microbiol. 2003;41(12):5456–65.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  53. Wu-Chuang A, Hodzic A, Mateos-Hernandez L, et al. Current debates and advances in tick microbiome research. Curr Res Parasitol Vector Borne Dis. 2021;1:100036.

    Article  PubMed  PubMed Central  Google Scholar 

  54. Narasimhan S, Fikrig E. Tick microbiome: the force within. Trends Parasitol. 2015;31(7):315–23.

    Article  PubMed  PubMed Central  Google Scholar 

  55. Thapa S, Zhang Y, Allen MS. Effects of temperature on bacterial microbiome composition in Ixodes scapularis ticks. Microbiologyopen. 2019;8(5):e00719.

    Article  PubMed  Google Scholar 

  56. 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(6):1779–96.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  57. Zhong Z, Zhong T, Peng Y, et al. Symbiont-regulated serotonin biosynthesis modulates tick feeding activity. Cell Host Microbe. 2021;29(10):1545–57. e4.

    Article  CAS  PubMed  Google Scholar 

  58. Sassera D, Lo N, Epis S, et al. Phylogenomic evidence for the presence of a flagellum and cbb(3) oxidase in the free-living mitochondrial ancestor. Mol Biol Evol. 2011;28(12):3285–96.

    Article  CAS  PubMed  Google Scholar 

  59. Guizzo MG, Tirloni L, Gonzalez SA, et al. Coxiella Endosymbiont of Rhipicephalus microplus modulates Tick Physiology with a major impact in blood feeding capacity. Front Microbiol. 2022;13:868575.

    Article  PubMed  PubMed Central  Google Scholar 

  60. Eberspächer J, Lingens F. The Genus Phenylobacterium. In: Dworkin M, Falkow S, Rosenberg E, et al. editors. The prokaryotes: volume 5: Proteobacteria: alpha and Beta subclasses. New York, NY: Springer New York; 2006. pp. 250–6.

    Google Scholar 

  61. Brites-Neto J, Duarte KM, Martins TF. Tick-borne infections in human and animal population worldwide. Vet World. 2015;8(3):301–15.

    Article  PubMed  PubMed Central  Google Scholar 

  62. Stich RW, Blagburn BL, Bowman DD, et al. Quantitative factors proposed to influence the prevalence of canine tick-borne disease agents in the United States. Parasit Vectors. 2014;7:417.

    Article  PubMed  PubMed Central  Google Scholar 

  63. Robinson MT, Satjanadumrong J, Hughes T, et al. Diagnosis of spotted fever group Rickettsia infections: the Asian perspective. Epidemiol Infect. 2019;147:e286.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  64. Laukaitis HJ, Macaluso KR. Unpacking the intricacies of Rickettsia-Vector interactions. Trends Parasitol. 2021;37(8):734–46.

    Article  PubMed  PubMed Central  Google Scholar 

  65. Willadsen P, Jongejan F. Immunology of the tick-host interaction and the control of ticks and tick-borne diseases. Parasitol Today. 1999;15(7):258–62.

    Article  CAS  PubMed  Google Scholar 

  66. Chauvin A, Moreau E, Bonnet S, et al. Babesia and its hosts: adaptation to long-lasting interactions as a way to achieve efficient transmission. Vet Res. 2009 Mar-Apr;40(2):37.

  67. Pechous RD, McCarthy TR, Zahrt TC. Working toward the future: insights into Francisella tularensis pathogenesis and vaccine development. Microbiol Mol Biol Rev. 2009;73(4):684–711.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

Download references

Acknowledgements

Not applicable.

Funding

This work was supported by National Science and Technology Major Program under Grant 2018ZX10101002-003; the National Key Research & Development Program of China under Grant 2022YFC2304004; the National Natural Science Foundation of China under Grant 82260399; the Science and Technology Leading Talent Team in Inner Mongolia Autonomous Region under Grant 2022LJRC0009; the Inner Mongolia Autonomous Region Science and Technology leading talent team: Zoonotic disease prevention and Control Technology innovation team under Grant 2022SLJRC0023; the Key Technology Project of Inner Mongolia Science and Technology Department under Grant 2021GG0171; the State Key Laboratory of Reproductive Regulation and Breeding of Grassland Livestock under Grant 2020ZD0008; the National Natural Science Foundation of China under Grant U22A20526; and the Guangxi University of Chinese Medicine Scientific Research Project under Grant 2020BS009.

Author information

Authors and Affiliations

Authors

Contributions

BX, WH, YW designed the study. BX, CS, JL, WBT, SWW. conducted all sampling, and extracted DNA. BX, QFW, CS, ZWZ, XLW performed the majority of experiments. MJG, QFW, BX, JWW analysed the data. BX, MJG, QFW wrote the paper. YZW, WH critically revised the manuscript. All authors discussed and approved the data and reviewed the manuscript.

Corresponding authors

Correspondence to Yuanzhi Wang or Wei Hu.

Ethics declarations

Ethics approval and consent to participate

All procedures performed in this study involving wild animals and human were in accordance with the ethical standards of the Animal Ethics Committee of Shihezi University (Approval No. A2018-143-01, KJ2020-062-01&A2020-113-01).

Consent for publication

Not applicable.

Competing interests

The authors declare no competing interests.

Disclosure statement

The authors report there are no competing interests to declare.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Electronic supplementary material

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Xu, B., Gu, M., Wu, Q. et al. The bacterial patterns suggesting the dynamic features of tick-associated microorganisms in hard ticks. BMC Microbiol 24, 179 (2024). https://doi.org/10.1186/s12866-024-03323-3

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12866-024-03323-3

Keywords