Worldwide, ticks are considered second to mosquitoes in their ability to transmit disease causing pathogens, and are the most relevant disease vectors for domestic and wild animals [1] . Among the tick borne zoonotic pathogens (TBZ) vectored by ticks are tick-borne encephalitis virus (TBEV), Crimean Congo hemorrhagic fever virus (CCHFV), Babesia spp., Rickettsia, Anaplasma, Francisela, Borrelia, Bartonella, Coxiella and Ehrlichia [2] , which they can transmit transstadially and in some cases transovarially [3] . In addition, an individual tick can harbor multiple pathogens which could be co-transmitted during feeding [4] . Other microbes co-exist in ticks as endosymbionts or commensals, and have been reported to influence vector survival and pathogen transmission fitness [5] . Until the advent of next generation sequencing (NGS), methods for detection targeted known pathogens [6] , making routine pathogen discovery efforts harder. NGS is pathogen agnostic and allows unbiased detection of sequences in a sample that can then be analyzed bioinformatically for matches in microbial gene repositories.

For this study, 16S rRNA gene, which is universally present in all bacteria, was used for identification of bacteria communities in ticks. The gene comprises 9 hyper-variable regions (V1 to V9) that are interspaced by conserved regions [7] [8] . Although the variable regions have considerable sequence diversity, no single region can differentiate all bacteria and a combination of 2 or more regions is recommended [7] [9] [10] . A combination of V3-V4 region that covers 460 base-pairs (bp) increases the accuracy of taxonomic classification to genus level compared to other variable regions [10] [11] [12] . Using these regions, the Illumina 300 cycles sequencing kit can be used to process 460 bp paired-ends reads [12] .

The 16S rRNA approach has been used for detection of pathogenic bacteria and symbionts in ticks [13] [14] [15] [16] [17] . These studies have reported up to 80% dominance of Proteobacteria followed by Actinobacteria, Bacteroidetes or Firmicutes [13] [14] . Potentially pathogenic genera that have been identified in these studies include Anaplasma, Coxiella, Ehrlichia, Rickettsia, Borrelia, Acinetobacter, Burkholderia and Staphylococcus among others [14] [15] [16] . In addition to harboring bacteria of medical relevance, ticks also carry endosymbionts such as Midichloria mitochloria, Wolbachia spp., Neoehrlichia spp., Rickettsiella. Endosymbionts very similar to tick-transmitted pathogens, including Coxiella-like, Rickettsia-like, or Francisella-like [13] [17] [18] that are said to provide additional nutrients lacking in blood meals have also been identified [19] [20] . Thus, 16S rRNA metagenomics is a useful approach for analyzing the complexity of bacteria community in ticks, as well as identifying potential TBZ.

This study used 16S rRNA metagenomics to examine bacterial community in ticks that were collected from livestock. We highlight the complexity of tick microbiome that they can acquire from the environment, livestock, wildlife or humans, some of them with potential to cause zoonoses.

The 16S rRNA MiSeq sequences output were de-multiplexed and adapters trimmed using the MiSeq reporter software version 2.6.3 (Illumina, CA USA). Mothur pipeline (version 1.35) was used for paired end reads contig assembly, sequence quality filtering, chimera removal and taxonomic assignment [24] . In brief, contigs containing ambiguous bases, and those with lengths < 350 bp (bp) or greater than 466 bp were discarded. Sequences were then aligned to a customized V3-V4 region on the SILVA database [25] , followed by merging sequences that were not more than 2 bp different from each other using the pre-cluster command in Mothur [24] . The merged sequences were then filtered for chimeras before taxonomic assignment against the customized V3-V4 SILVA database [24] using Bayesian classifier with 80% confidence [24] [25] . Unassigned operational taxonomic unit (OTU) and those assigned to Chloroplast, Mitochondria, Archaea, and Eukaryote were discarded. Samples with less than 1000 (n = 40) sequences were excluded from downstream analysis as small library sizes often conceal biologically meaningful results [26] . Taxa detected in the non-template control were censored from the tick sample dataset [27] .

Statistical analysis and data visualization were done from OTU tables generated by Mothur on R software environment version 4, with Phyloseq, vegan and ggplot2 statistical adds-on [28] [29] [30] [31] . To down sample the data for alpha diversity calculation and account for unequal sequencing between samples, rarefaction was done using rarefy_even_depth command in phloseq with replacement [32] . The rarefied data was used to determine Shannon diversity index by first determining bacterial diversity for each tick species before computing the mean [33] .

This study presents data on bacterial communities associated with ticks of the genus Rhipicephalus, Amblyomma and Hyalomma that were collected from domestic animals (cattle, sheep and goats). Ticks belonging to these genera feed on multiple hosts, including wildlife and domestic livestock, in addition to being homophilic, vector multiple zoonotic pathogens [3] , and are therefore very important component of “One Health”. To the best of our knowledge, our study is the first in Kenya to report on the 16S rRNA approach in tick microbial diversity.

Four bacteria phyla carried nearly all the tick microbiome (Figure 1). The patterns of bacterial phyla is consistent with findings from other studies that reported over-abundance of Proteobacteria, Firmicutes and Actinobacteria, singly or in different combinations [14] [17] [18] [34] [35] . Of these four, Proteobacteria that comprises Gram negative bacteria, many of them pathogenic to humans and animals [36] [37] accounted for 33.3% of the OTUs. Firmicutes comprise Gram-positive bacteria that include notable pathogens that are found in different environments [38] and were the second abundant phylum at 23.4% of the OTUs. Actinobacteria, Gram-positive bacteria that comprise species that are plant and animal pathogens inhabits diverse environments [39] [40] was the third most abundant and contained 20.0% of the OTUs. Our finding is similar to other reports that indicate Actinobacteria as the third most common phylum after Proteobacteria and Firmicutes [41] . Bacteroidetes was the 4^th dominant phylum and constituted 11.6% of the OTUs. Bacteria in this phylum are Gram negative and inhabit many diverse habitats, where they play a beneficial role in the degradation of organic matter. Nevertheless, some members of the phylum are commonly isolated in human and animal clinical samples, including blood, urine, wounds, and feces [42] [43] . The remaining 23 phyla contributed 11.2% of the OTUs. Members of these phyla have also been reported in other tick studies [44] .

Amongst the 645 genera, three groups were discernible (Figure 2). The biggest group comprised commensals/symbionts (Figure 2, Panel A) that contributed 93.6% of the genera, but their individual sequence contribution was very low. Group two comprised genera that are known to contain potentially pathogenic species that were probably picked from the livestock during the ticks’ blood meals (Figure 2, Panel B). Of these, Coxiella was the most abundant and contributed 41.8%% of the total sequence reads. In studies conducted in Malaysia and China, Haemaphysalis ticks were found to contain overabundance of Coxiella [34] [45] . Coxiella also dominated in Amblyomma ticks collected in the United States [46] . Due to inability of the 16S rRNA V3-V4 approach to resolve genera to species, the Coxiella in ticks could contain pathogenic or symbionts species [11] . We opted to error on the side of caution and referred the genera as potentially pathogenic. Corynebacterium in the phylum Actinobacteria, is widely distributed in nature as part of animal and human microbiota, but some species such as C. diphtheria cause human infections [40] . This genus was the second most abundant at 13.9%. High abundance of Corynebacterium was found in Hyalomma ticks from United Arab Emirates (UAE) [47] . Other potentially pathogenic bacteria identified included Enterococcus, Staphylococcus, Klebsiella, and Acinetobacter (ESKAPE) that comprise species with high rates of antibiotic resistance [48] . Others included Acinetobacter (4.3%) and Staphylococcus (3.8%) and have been reported to occur in Ixodes, Amblyomma and Rhipicephalus ticks [13] [45] [49] . Similarly, Enterococcus and Klebsiella have been reported in low prevalence in ticks collected from dogs in China [13] . Although pathogenic, Burkhoderia, Escherichia-Shegella, Achromobacter, Haemophilus, Legionella, Campylobacter, Treponema, Elizabethkingia, Mycoplasma, Bordetella, Vibrio, and Brucella were present in low abundance. Multiple reports indicate presence of these genera in ticks [13] [44] [50] [51] . Group three comprised genera known to contain TBZ pathogens (Figure 2, Panel B). Rickettsiae, Anaplasma, Francisella, Erlichia, Bartonella and Borellia individually contributed <1% of the sequences (Figure 2 Panel C). Unlike other tick microbiome studies that have reported dominance of Rickettsia and Borrelia [14] [18] , in this study, these pathogens occurred at <1% relative abundance. As shown in Figure 2, panel C, Rickettsia was identified in all tick species accessed but notably in lower prevalence in Hyalomma compared to Amblyomma and Rhipeciphalus. Similar studies in the United Arab Emirates (UAE) and Saudi Arabia have reported absence of Rickettsia in Hyalomma tick collected from camels [47] [52] . Similar to this finding, Anaplasma and Ehrlichia have been reported in relatively low abundance in Tennessee, USA for example [44] . Francisella has also been shown to be dominant in Hyalomma ticks comprising up to 99.1% in some locations in the UAE [47] , but in our study, Francisella was present in much lower abundance (relative abundance of 0.005) and was absent in Hyalomma ticks. Similar to studies in Germany and China [13] [53] , Bartonella was detected in low abundance and only in Amblyomma spp and H. truncatum. As a caveat, we cannot say with certainty that the Rickettsia and Francisella identified in this study contained pathogenic species [5] [54] [55] . But, for Rickettsia, we know that Kenyan ticks carry a high prevalence of Rickettsia africae, R. aeschlimannii, R. mongolotimonae, R. conorii subsp. israelensis), Candidatus Rickettsia kulagini, and other unspeciated Rickettsiae [21] [56] .

Bacterial community in different tick species varies, probably as a result of environmental factors, blood-meals sources, tick immunity and developmental stages [57] . Contamination at different steps of the DNA extraction, purification and amplification process has been identified as a contributor to the high bacteria diversity seen in ticks [58] . To offset this problem, negative controls are used to track contamination, and if identified, censored from the dataset. Using the “decontam” approached [27] , 6 OTUs were dropped from the dataset. Clearly, since the ticks were not surface sterilized before DNA extraction, exo-skeleton associated bacteria contributed to some of the community observed. But, as noted by Narasimhan et al., exo-skeleton associated microbes should be considered as being part of the ticks and they probably play a key role in maintaining a healthy physical immune barrier [59] .

As shown in Figure 3, the aggregated mean Shannon diversity index was low (ranged between 2.69 to 4.15), suggesting that the bacteria genera were shared among the tick species. Similar findings have been reported previously and similar to these studies, a few core bacteria taxa, likely endosymbionts dominate [60] [61] [62] .

16S rRNA was used to identify bacterial communities associated with different tick species that were collected from domestic livestock in Kenya. There was little difference in bacteria diversity between the ticks, probably because of similarity in the environment and the hosts that the ticks interact with. Amidst the hundreds of commensal taxa, six genera known to contain TBZ pathogens, namely Rickettsiae, Anaplasma, Francisella, Ehrlichia, Bartonella and Borrelia were identified and illustrate the central role ticks play in “one health” and their usefulness as bio-indicators of pathogens they likely transmit to humans, domestic livestock and wildlife. Future studies should, in addition to prokaryotic 16S rRNA, include eukaryotic 18S rRNA and viral discovery in order to comprehensively monitor microbial ecosystem in the ticks.

The tick samples were collected under protocol SSC#1248 that was reviewed and approved by the Animal Use Committee of the Kenya Medical Research Institute.

This work is published with the permission of the Director, Kenya Medical Research Institute. This work was funded by the Armed Forces Health Surveillance Division (AFHSD) and its Global Emerging Infections Surveillance and Research Branch (insert ProMIS ID P0136_19_KY, 2019 and P0095_21_KY, 2021.

Material has been reviewed by the Walter Reed Army Institute of Research. There is no objection to its publication. The opinions or assertions contained herein are the private views of the author, and they are not to be construed as official, or as reflecting true views of the Department of the Army or the Department of Defense.

The authors declare no conflicts of interest regarding the publication of this paper.

Mutai, B., Njaanake, K., Gathii, K., Estambale, B.B. and Waitumbi, J.N. (2022) Bacteriome in Ticks Collected from Domestic Livestock in Kenya. Advances in Microbiology, 12, 67-82. https://doi.org/10.4236/aim.2022.122006