The Bacterial Community in Questing Ticks From Khao Yai National Park in Thailand

Ticks are known vectors for a variety of pathogens including bacteria, viruses, fungi, and parasites. In this study, bacterial communities were investigated in active life stages of three tick genera (Haemaphysalis, Dermacentor, and Amblyomma) collected from Khao Yai National Park in Thailand. Four hundred and thirty-three questing ticks were selected for pathogen detection individually using real-time PCR assays, and 58 of these were subjected to further metagenomics analysis. A total of 62 ticks were found to be infected with pathogenic bacteria, for a 14.3% prevalence rate, with Amblyomma spp. exhibiting the highest infection rate (20.5%), followed by Haemaphysalis spp. (14.5%) and Dermacentor spp. (8.6%). Rickettsia spp. were the most prevalent bacteria (7.9%) found, followed by Ehrlichia spp. (3.2%), and Anaplasma spp. and Borrelia spp. each with a similar prevalence of 1.6%. Co-infection between pathogenic bacteria was only detected in three Haemaphysalis females, and all co-infections were between Rickettsia spp. and Anaplasmataceae (Ehrlichia spp. or Anaplasma spp.), accounting for 4.6% of infected ticks or 0.7% of all examined questing ticks. The prevalence of the Coxiella-like endosymbiont was also investigated. Of ticks tested, 65.8% were positive for the Coxiella-like endosymbiont, with the highest infection rate in nymphs (86.7%), followed by females (83.4%). Among tick genera, Haemaphysalis exhibited the highest prevalence of infection with the Coxiella-like endosymbiont. Ticks harboring the Coxiella-like endosymbiont were more likely to be infected with Ehrlichia spp. or Rickettsia spp. than those without, with statistical significance for Ehrlichia spp. infection in particular (p-values = 0.003 and 0.917 for Ehrlichia spp. and Rickettsia spp., respectively). Profiling the bacterial community in ticks using metagenomics revealed distinct, predominant bacterial taxa in tick genera. Alpha and beta diversities analyses showed that the bacterial community diversity and composition in Haemaphysalis spp. was significantly different from Amblyomma spp. However, when examining bacterial diversity among tick life stages (larva, nymph, and adult) in Haemaphysalis spp., no significant difference among life stages was detected. These results provide valuable information on the bacterial community composition and co-infection rates in questing ticks in Thailand, with implications for animal and human health.


INTRODUCTION
Ticks are recognized as a medically important group of arthropods that transmit a number of diseases to humans (1). Different tick species favor different habitats, which ultimately defines their geographical distribution and thus the risk areas for human or animal infections. Several pathogens are known to be carried by either hard ticks (Ixodidae) or soft ticks (Argasidae), including a wide range of viruses, bacteria, fungi, and protozoa (2)(3)(4)(5)(6).
Ticks also harbor many non-pathogenic organisms. Bacterial endosymbionts have been recognized as important microorganisms required for tick fitness, especially with regard to regulating host reproduction and immunity (18,19). Some studies suggest that these symbionts have a potential role in providing key vitamins absent in the bloodmeal (20,21). Others reported that they might have an important role in facilitating pathogen colonization in the gut as in the case in Ixodes scapularis ticks, in which alteration of the symbiont abundance resulted in decreased Borrelia burgdorferi colonization (22). Likewise, the level of Anaplasma marginale acquisition was lower in Dermacentor andersoni when their microbiome was altered by antibiotic treatment leading to an increase in proportion and quantity of Rickettsia bellii in the microbiome (23). Moreover, the prevalence and transovarial transmission of bacterial endosymbionts occurs at a high rate, suggesting that they might have an obligate relationship with the host (24,25). In addition to the bacterial microbiota in ticks, other microorganisms are as important and abundant. For example, several studies reported that two groups of bunyaviruses (South Bay virus and Phlebovirus) were commonly associated with Ixodes spp. ticks in America and Europe (26)(27)(28) and were more abundant than bacteria or eukaryote counterparts (29).
It is thought that because ticks harbor a diverse range of microorganisms, co-infection or co-occurrence of bacteria, parasites, and viruses in ticks is possible (30). Co-infection may lead to increased disease severity as it complicates disease diagnosis and treatment (31)(32)(33)(34). Co-infections between B. burgdorferi and Anaplasma phagocytophilum are widely recognized and reported (35). Other co-infections among different types of microorganisms have also been reported, such as between B. burgdorferi and South Bay virus (SBV), Babesia microti, and B. burgdorferi, as well as between a novel filarial worm (Onchocercidae sp. ex. Ixodes scapularis) and Wolbachia spp (29).
In this study, we used both conventional and high-throughput sequencing methods to study bacterial pathogen co-infections and pathogen association with bacterial endosymbionts in questing ticks collected in Khao Yai National Park. Metagenomics were used to determine the bacterial communities in individual ticks and conventional methods (real-time PCR, PCR, and Sanger sequencing) were used to detect pathogenic bacteria: Rickettsia spp., Anaplasma spp., Ehrlichia spp., Borrelia spp., Coxiella burnetii, and pathogenic Francisella spp. in Amblyomma spp., Dermacentor spp., and Haemaphysalis spp. ticks individually. All positive samples underwent DNA sequencing to identify pathogens to the species level. Coinfections between bacteria and the association between pathogenic bacteria and the Coxiella-like endosymbiont were investigated.

Tick Collection
Questing ticks were collected by dragging a 1-m 2 cotton cloth over vegetation in four tourist attraction sites in Khoa Yai National Park (14 • 26'19.5"N 101 • 22'20.1"E). Sampling was conducted by six people at each site for 1 h. Ticks were collected in one trip in November 2020 and immediately preserved in 90% ethanol before transporting to the laboratory for further processing. All sampling procedures and experimental manipulations were reviewed and approved as part of the animal collection protocol entitled "Surveillance of Tickand Flea-Borne Diseases of Public Health Importance" (PN# 21-01). The project was also approved by Mahidol University -Institute animal care and use committee (MU-IACUC 2019/3). Research was conducted in compliance with the Animal Welfare Act and other federal statutes and regulations related to animals and experiments involving animals, and adhered to principles outlined in the "Guide for the Care and Use of Laboratory Animals, " NRC Publication, 2011 edition.

Tick Surface Sterilization and DNA Extraction
Questing ticks (unfed ticks) were morphologically identified using taxonomic keys (36,37). Each tick was vortexed for 1 min in 3% sodium hypochlorite and then transferred to a new tube and vortexed for 1 min in 70% alcohol followed by three washes in sterile PBS in the same manner. Ticks were air dried for 10 min on Whatman R filter paper (Sigma-Aldrich, Saint Louis, MO, USA) before DNA extraction. Whole ticks in 250 µl of ATL buffer (lysis buffer, component of the extraction kit) were punctured with a fine tip under a stereomicroscope to release the tissue from the hard chitin exoskeleton prior to adding 2 mg/ml of Proteinase K solution. Samples were then incubated at 55 • C overnight. A total volume of 250 µl of homogenized solution was then used for DNA extraction on the QIAsymphony R SP instrument with QIAsymphony R DSP DNA Mini Kit using Tissue LC 200 DSP protocol (Qiagen, Hombrechtikon, Switzerland). The DNA was eluted in 50 µl of ATE buffer (elution buffer, component of extraction kit) and stored at−20 • C until use. Ultrapure DNA/RNA-free distilled water as well as PBS buffer used for tick surface sterilization were also included as an extraction control.

Amplification of Bacterial 16S rDNA
Following DNA extraction, the bacterial-specific 16S rDNA (V3-V4, a 550-bp fragment) was amplified as previously described (15). Negative control PCR reactions were included in every experimental run using Ultrapure DNA/RNA-free distilled water in place of DNA template. PCR reactions were also performed with eluates from mock DNA extractions as well as from PBS buffer used for tick surface sterilization. PCR product was cleaned using AMPure magnetic beadbased purification system (Beckman Coulter, UK) following the manufacturer's instructions. Purified PCR products were eluted and quantified using the Qubit dsDNA HS Assay Kit (Invitrogen Life Technologies, MA).

Library Preparation and High-Throughput Sequencing
The library was prepared with dual indices and Illumina sequencing adapters attached to purified PCR products using the Nextera XT Index Kit following the manufacturer's protocol (Illumina Inc., San Diego, CA). For index control reaction, a combination of index primers that were not used with samples was also included with PCR grade water as template. The number of reads recovered from these particular index combinations were used to filter the cross-contaminations between indexed PCR primers and to identify errors in an Illumina sample sheet. Libraries were cleaned using Agencourt AMPure XP beads. The purity of the libraries was checked on the QIAxcel Advanced System (Qiagen) with a QIAxcel DNA High Resolution Cartridge. Purified amplicon libraries were quantified using the Qubit dsDNA HS Assay Kit (Invitrogen). DNA concentration was calculated and normalized to reach 4.0 nM for each library. Five microliters of DNA from each library was pooled for a NGS run. Pooled libraries were denatured and diluted to a final concentration of 8 pM with a 10% PhiX (Illumina) control. Sequencing was performed using the MiSeq Reagent Kit V3 on the Illumina MiSeq System.

Data Analysis for Metagenomics and Diversity Estimates
The sequence reads generated by the 16S rRNA on MiSeq sequencers were processed on CLC Genomics workbench v 12.0.3 (Qiagen, Aarhus A/S, http://www.clcbio.com). Highthroughput sequences were imported into CLC Genomics Workbench according to quality scores of Illumina pipeline 1.8. In order to achieve the highest quality sequences for clustering, paired reads were merged in CLC microbial genomics module v 4.8 using default settings (mismatch cost = 1; minimum score = 40; gap cost = 4 and maximum unaligned end mismatch = 5). Primer sequences were trimmed from merged reads using parameters (trim using quality scores = 0.01, trim ambiguous nucleotides = 2, and discard read length shorter than 150 bp). Samples were removed from analysis if the number of reads was <100 or <25% from the median (the median number of reads across all samples). Chimeric sequences were detected and removed. Only filtered and merged sequences were clustered into operational taxonomic units (OTUs) according to a threshold of 97% sequence identity. All such processes were performed using CLC microbial genomics module v 4.8. Reference OTU data used in the present study were downloaded from the Greengenes database V13.8 (38). OTUs with combined abundance <10 reads were removed from the downstream analysis. Alpha diversity estimates (Observed OTUs, Simpson's index, and Shannon entropy) were analyzed on the quality-filtered OTU table at the genus level using CLC Microbial Genomics Module v 4.8. Beta diversity analysis for microbiome compositional difference between groups was calculated using a distancebased non-parametric test, the generalized UniFrac distances (39). The statistical significance of microbiome compositional difference between groups (tick genera, tick developmental stages, Francisella persica infection, Coxiella-like endosymbiont infection, and Rickettsia spp. infection statuses) was then compared using PERMANOVA using 99,999 permutations of the distance values. All analyses mentioned here were performed with CLC Microbial Genomics Module v 4.8.

Bacterial Pathogen Detection by Real-Time PCR and Conventional PCR
Real-time PCR and PCR assays were performed on 433 individual ticks for the detection of bacteria (Rickettsia spp., Borrelia spp., Anaplasma spp., Ehrlichia spp., and Coxiella-like endosymbiont) and the taxonomic species assignment. Other potential pathogenic bacteria (Francisella spp. and Coxiella spp.) detected by NGS analysis (read count > 1) were also confirmed by real-time PCR and PCR assays. Detailed methods for assays and target gene(s) for selected pathogens are provided as online Supplementary data (Supplementary Table 1). For all real-time PCR, the reaction consisted of 1X Platinum quantitative PCR SuperMix-UDG (Invitrogen) using standard

DNA Sequencing
PCR amplicons were purified using the QIAquick R PCR Purification Kit (Qiagen) according to the manufacturer's instructions. The PCR products were cycle-sequenced using an ABI BigDye TM Terminator v3.1 Cycle Sequencing Kit, ethanol precipitated, and run on a SeqStudio Genetic Analyzer (Applied Biosystems ThermoFisher, Thailand). Sequences of each sample and pathogen were assembled using Sequencher TM ver. 5.4.6 (GeneCodes Corp., Ann Arbor, MI). The pathogen sequences were aligned with reference sequences retrieved from the GenBank database using the MUSCLE codon alignment algorithm (40). A maximum likelihood phylogenetic tree was then constructed from bacterial target genes (Supplementary Table 1) using the best fit model of nucleotide substitution with bootstrapping (1,000 replicates) in MEGA 6 (41).

Statistical Analysis
Differences in alpha diversity indices of the bacterial community composition, based on metagenomics data, were determined by Mann-Whitney U test (between two groups) or Kruskal-Wallis test (across all groups) and the critical range (p < 0.05) was determined. Statistical analyses (two-way ANOVA tests and mean, 95% confidence interval) and scatter plots were performed using GraphPad Prism version 5.04 for Windows (GraphPad Software, San Diego, CA. www.graphpad.com). Some graphical illustrations presented in this study as well as the Chi-square independence test and Fisher's exact test were performed in the R environment for statistical computing (42,43). A nucleotide distance matrix was generated using "Compute Pairwise Distance" in MEGA 6 (41).

Co-infection of Bacterial Pathogens
Co-infection between pathogenic bacteria examined in this study only occurred in a small number of ticks, primarily female Haemaphysalis spp. (n = 3). All co-infections occurred with Rickettsia spp. and Anaplasmataceae (Ehrlichia spp. or Anaplasma spp.) (Figures 2A,B) and accounted for 4.8% of infected ticks and 0.7% of all ticks examined. All three co-infected ticks were also positive for the Coxiella-like endosymbiont.

Bacterial Endosymbionts and Co-occurrence With Pathogenic Bacteria
The Coxiella-like endosymbiont was detected in 65.8% of all ticks examined, with the highest infection rate in nymphs (86.7%) and females (83.4%) ( Table 4). Of all genera, Haemaphysalis spp. exhibited the highest rate of infection with the Coxiella-like endosymbiont, with an 81.0% infection rate. The majority of Haemaphysalis spp. infected were females (95.1%), followed by nymphs (88.9%), males (87.0%), and, to a lesser extent, larvae (24.0%). The second highest infection       Table 1). The classification of OTUs from each sample was made against the Greengenes reference database and the similarity threshold was set at 0.97 in the CLC microbial genomics module. There were nine recorded phyla found among all tick samples studied: Proteobacteria (91%), Bacteroidetes (4%), Firmicutes (3%), Cyanobacteria (1%), and Actinobacteria (1%) (Figure 3A). There were five major phyla found in controls for DNA extraction, PCR, and Indexing; however, the majority of Proteobacteria phylum (73%) detected was genus Enterobacteria (89%), and very small amount of genus Rickettsia (0.0092%) (Supplementary Figure 2). Only 1% of reads could not be classified using Greengenes reference and were removed before performing downstream process of abundance analysis. Figure 3B shows the abundance of bacterial taxa in tick genera at the family level. There were one or two predominant bacterial taxa in each tick genus. Rickettsiaceae (83%) was the predominant bacterial taxon in Amblyomma spp., while Francisellaceae (50%) and Methylobacteriaceae (26%) were the major taxa found in Dermacentor spp. Haemaphysalis spp. harbored a more diverse bacterial spectrum than the other two tick genera. Coxiellaceae (43%) and Methylobacteriaceae (24%) were the predominant bacterial taxa found, with a lesser abundance of Rickettsiaceae (4.6%). Beta analysis using distance-based non-parametric test indicated that microbiome composition differed significantly across tick genera by the generalized UniFrac distance (df = 2, pseudo-F = 5.812, p-value = 0.00001). Alpha diversity estimates of bacterial communities in ticks and statistical significance (p-value) for group comparison are summarized in Table 5. The alpha diversity analyses showed that the bacterial community of Haemaphysalis spp. (n = 40) was significantly different from Amblyomma spp. (n = 8) with p-value < 0.05 for all tests (Observed OTUs, Simpson's index, Shannon entropy) when measured at the genus level as shown in Figure 4A. While there were few ticks positive for Francisella persica (n = 4), alpha diversity estimates indicated that their bacterial profiles were quite different from ticks that were not infected with F. persica with statistical significance (p-value < 0.001) ( Figure 4B).
As three Haemaphysalis spp. life stages (Larva = 9, Nymph = 15, Male = 6, Female = 10) were included in NGS, this genus was selected for testing the bacterial community difference among tick life stages (Figures 5A,B). Alpha diversity estimates (Simpson's index and Shannon entropy) showed that while the overall bacterial profiles among life stages were not significantly different with regard to bacterial taxa present, the number of observed OTUs among life stages was significantly different (Table 5; Figure 5B). Likewise, the beta diversity analysis using distance-based non-parametric test showed no significant difference across Haemaphysalis spp. life stages by the generalized UniFrac distance (df = 3, pseudo-F = 0.905, p-value = 0.5638). Comparing Coxiellaceae abundance across life stages, this group comprised 57-59% abundance of the total bacteria community in male and female ticks, 40% in nymphs, and 28% in larvae ( Figure 5A). The statistically significant differences between groups are indicated (*p-value < 0.05, **p-value < 0.001). The solid lines show mean and 95% confidence interval (CI) of alpha diversity for each group. Scatter plots and Mann-Whitney U test were created and tested using GraphPad Prism version 5.04.

DISCUSSION
In this study, we examined the bacterial community profile and co-infection of questing ticks collected by dragging from Khao Yai National Park using metagenomics. The majority of ticks collected were in the larval stage; however, adults and nymphs were also collected and identified as Haemaphysalis spp., Amblyomma spp., and Dermacentor spp. The greater number of Haemaphysalis spp. found in comparison to other genera may be due to the collection method used (dragging). The prevalence of Rickettsia spp. in ticks collected in Khao Yai National Park in this study was lower than in a previous report that found 30% of Amblyomma testudinarium positive for Rickettsia spp. and approximately 17% positive in Haemaphysalis spp. (44). However, they did not detect other bacteria such as Borrelia spp., Francisella spp., or the common symbiont Wolbachia spp. In this study, we detected Rickettsia DNA in 18% and 8% in Amblyomma spp. (n = 44) and Haemaphysalis spp. (n = 331) respectively, along with other bacteria such as Anaplasma spp., Ehrlichia spp., and Borrelia spp. at lower prevalence rates. We did not detect Wolbachia spp. in any of the ticks tested but found other bacterial endosymbionts, namely Rickettsia spp., Francisella spp., and Coxiella spp., using metagenomics and real-time PCR assay. Sequence analyses of Coxiella spp. and Francisella spp. detected by NGS confirmed the detection of these bacterial endosymbionts in examined ticks (Supplementary Figure 1). Rickettsia spp. that were detected by NGS but were not confirmed by subsequent qPCR and PCR assays are suspected to be a Rickettsia endosymbiont.
Co-infection of pathogenic bacteria occurred in a small number of ticks (n = 3), with the majority of co-infections between Rickettsia spp. and the Anaplasmataceae family. High prevalence of Rickettsia spp. in domestic animals and their ectoparasites, especially fleas, in Thailand has been recognized (15,16,45). Less than 10% of Rickettsia spp. was found in ticks collected from animals (45). However, some studies found infection rates as high as 24% for Rickettsia spp. and 32% for Anaplasma spp. in adult ticks collected from under leaves along animal trails across the country (46). The finding of co-infections between Rickettsia spp. and other bacteria in female ticks in this study was likely due to the ability of the bacteria to be maintained in ticks through transovarial and transstadial transmission (47)(48)(49). As females have already taken two bloodmeals during their development from larva to adult, the probability for Rickettsia-infected or Anaplasmataceae-infected ticks acquiring additional bacteria during feeding would be high (50)(51)(52). Our study also revealed the increased infection rate in adult ticks compared with larvae or nymphs. Similar findings were reported for an increasing A. phagocytophilum infection in adult ticks, especially in females, in Hanover, Germany (53). The study also reported that the co-infection between A. phagocytophilum and Rickettsia spp. was higher in females (5.2%) than in males (2.4%) or nymphs (1.6%).
A similar study was conducted by Nooroong et al. (46), in which ticks were collected under leaves along animal trails across Thailand were screened for bacterial pathogens. They found co-infections between Rickettsia spp. and Anaplasma spp. This co-infection was also observed in adult A. testudinarium ticks carrying the Coxiella-like endosymbiont in Nakhon Nayok, the same province where Khao Yai National Park is located. In Nooroong et al. (46), only adult ticks were collected and a 5.97% co-infection rate between Rickettsia spp. and Anaplasma spp., as well as a total infection rate of 35.8% for Rickettsia, was reported in the location near our study. This higher infection rate is likely due to the fact that Nooroong et al. (46) focused on the adult stage and did not include data from immature ticks. While all tick life stages collected by dragging were included in our study, the majority of ticks collected were larvae. Several previous studies in Thailand have primarily focused on ticks collected from animals; therefore, true co-infection status is inconclusive as pathogens may have been from the host animal. Our study focuses on coinfections in questing ticks actively seeking hosts, which directly represents the human risk of encountering ticks that can transmit multiple pathogens simultaneously. The co-infections investigated in previous studies were mostly among pathogenic bacteria and/or with bacterial endosymbionts (54), or among bacteria and protozoa such as Coxiella and Babesia spp. in H. bispinosa (55), or even co-infection with two to four pathogens of Anaplasmataceae, Babesia spp., and Hepatozoon spp. in Rhipicephalus sanguineus sensu lato collected from dogs in Bangkok, Thailand (56). Co-infection with more than three microorganisms was also reported in questing Ixodes scapularis ticks in Wisconsin, USA (29). Additionally, another study reported by Moutailler et al. (57) found up to 45% of questing I. ricinus ticks were co-infected with five to eight different pathogens including bacteria (Borrelia spp.), parasites, viruses, and endosymbionts. No significant interactions between endosymbionts and pathogens were found but there was a significant association between two Borrelia species: B. garinii and B. afzelii. Our results show that Ehrlichia infection in ticks harboring Coxiella-like endosymbiont was significantly higher than those without, suggesting that there might be some kind of relationship between Ehrlichia spp. and the Coxiellalike endosymbiont. However, we found no other association between Anaplasmataceae (Ehrlichia spp. and Anaplasma spp.) or Rickettsia spp. infection with the Coxiella-like endosymbiont. The data observed in this study implied that there was no negative or competitive effect of the Coxiella-like endosymbiont on the maintenance or existence of other bacteria in tick hosts. Endosymbionts relate to their host and pathogenic bacteria in many ways, such as providing nutrition lacking in blood meal, and sometimes possess an obligate relationship with the host (58,59). However, other endosymbionts may interfere with the transmission to vertebrate hosts as was shown in R. rickettsii and R. peacockii in D. andersoni (60,61). Another example of transmission interference can be seen in the salivary glands of Amblyomma spp., where a Coxiella-related symbiont impairs the transmission of E. chaffeensis (62). Alteration of the bacterial microbiome can also interfere with the colonization of pathogenic bacteria as seen in B. burgdorferi by modulation of the host immune response (63), which indicates that gut microbiota in ticks also plays a role in pathogen colonization in the gut lumen (22,63).
Bacterial profiles using metagenomics used in this study showed a few dominant bacterial taxa harbored by questing ticks collected in Khao Yai National Park. Three main bacteria taxa were found in each tick genus-Coxiella spp. in Haemaphysalis spp., Rickettsia spp. in Amblyomma spp., and Francisella spp. in Dermacentor spp., making up the majority of bacterial taxa in ticks with relative abundance ranged from 40 to 80%. Our findings were consistent with other studies reporting that hard tick microbiomes are dominated by a small number of bacterial species, most of which are endosymbionts (20,64). For example, Coxiella spp. was the main taxon found with a relative prevalence of 89.5% in Rhipicephalus turanicus (65), 89-100% for A. americanum (66)(67)(68), 98.2% in the female ovaries of Rhipicephalus (Boophilus) microplus (69), and 39.2% in Haemaphysalis spp. collected from animals in Malaysia (64). Coxiella spp. is a ubiquitous bacterium found in many tick species and maintained through transovarial transmission, as it was shown to pass on to their eggs and larvae from adult laboratory-reared R. sanguineus ticks (65,70). Other bacterial endosymbionts detected in our study were Rickettsia spp., and depending on the tick species, Rickettsia spp. could represent up to 83% of the bacterial community, as observed in Amblyomma ticks in our study, while it was found to be less prevalent in the other two genera (Haemaphysalis and Dermacentor). A Francisella-like endosymbiont was another taxon we detected in high prevalence along with Methylobacteria in Dermacentor spp. The same finding was previously reported by several groups in which two Dermacentor species were dominated by three core bacterial taxa: Francisella, Sphinogomonas, and Methylobacterium (61,71,72). In addition to Dermacentor spp., A. maculatum was reported to harbor a Francisella-like endosymbiont at high abundance as well (73,74). Our study also found that bacterial community composition varied significantly among tick genera, especially between Amblyomma and Haemaphysalis ticks. However, the species richness and diversity slightly decreased during development in Haemaphysalis spp. Other studies found that microbiome richness and diversity significantly decreased during development and varied greatly among species (72,75). In one study, the obvious difference was on core OTUs of endosymbiont bacteria among tick species where Dermacentor spp. was dominated by Francisella spp., while H. leporispalustris and I. pacificus had Coxiella spp. and Rickettsia spp. as dominant bacterial species, respectively (72). Differences in reports might be from the techniques used to study the microbiome, geography, and natural vs. laboratory-reared tick populations (22,(76)(77)(78). However, the striking similarity among these studies is the difference in bacterial endosymbionts among tick genera, suggesting that there could be a competition among symbionts within tick genera. As it was reported in R. turanicus, Coxiella was the primary symbiont and Rickettsia was the secondary, having lesser relative abundance (65). Other studies discovered competition between Rickettsia species in Dermacentor spp. in which one species prevented another from transovarial transmission (79,80).
Understanding the microbiome composition in ticks of different species, their vector capacities, as well as the role of bacterial endosymbionts on tick physiology, including their influence on pathogen transmission, may provide insight into vector control to prevent human infection and the emergence of tick-borne diseases in the future. This study provides important information on bacterial community composition and coinfection rates in questing ticks in Thailand with implications for animal and human health.

DATA AVAILABILITY STATEMENT
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found below: GenBank accession number(s) OK161258-61, OK161265-68, OK161272-79, OK161279-84, and OK180561-611. The SRA BioProject ID PRJNA766341.

ETHICS STATEMENT
The animal study was reviewed and approved by AFRIMS