Anaplasma and Theileria Pathogens in Cattle of Lambwe Valley, Kenya: A Case for Pro-Active Surveillance in the Wildlife–Livestock Interface

Tick-borne pathogens (TBPs) are major constraints to livestock production and a threat to public health in Africa. This cross-sectional study investigated the risk of infection with TBPs in cattle of Lambwe Valley, Kenya. Blood samples of 680 zebu cattle from 95 herds in six geospatial clusters within 5 km of Ruma National Park were screened for bacterial and protozoan TBPs by high-resolution melting analysis and sequencing of PCR products. We detected Anaplasma bovis (17.4%), Anaplasma platys (16.9%), Anaplasma marginale (0.6%), Theileria velifera (40%), and Theileria mutans (25.7%), as well as an Anaplasma sp. (11.6%) that matched recently reported Anaplasma sp. sequences from Ethiopia. Babesia, Rickettsia, and Ehrlichia spp. were not detected. The animal and herd-level prevalences for TBPs were 78.5% (95% confidence intervals (CI): 75.3, 81.5) and 95.8% (95% CI: 91.8, 99.8), respectively. About 31.6% of cattle were co-infected with 13 combinations of TBPs. The prevalence of TBPs differed between clusters and age, but the risk of infection was not associated with sex, herd size, or the distance of homesteads from Ruma. This study adds insight into the epidemiology of TBPs around Ruma and highlights the need for proactive surveillance of TBPs in livestock–wildlife interfaces.


Introduction
Wildlife-livestock interfaces, or areas where game and domesticated animals co-habit [1], support a plethora of interactions that can be of epidemiological concern [2]. Such interfaces are a haven for ticks and many other arthropod vectors of diseases maintained by wildlife [1]. Global trends suggest that emerging infectious diseases are on the rise, with up to 80% of animal pathogens and 71.8% of emerging zoonoses in some regions predicted to have a wildlife component [3,4]. Over the past two decades, there has been a heightened incidence of ticks and tick-borne pathogens (TBPs) across a wide geographical spread [3] and new tick-borne disease agents often discovered [5]. Many of these infectious diseases transmitted by ticks between domestic and wild animals represent emerging and re-emerging burdens to global public health, economies, and the conservation of biodiversity [2,6]. Astonishingly, there is limited information on the diversity and prevalence of TBPs in many areas, including the wildlife-livestock interfaces. Where available, data are dated and might not be beneficial for the strategic management and control of TBPs [6].
In sub-Saharan Africa, theileriosis, ehrlichiosis, babesiosis, and anaplasmosis spread by ticks are the leading causes of losses in livestock production [7]. Several TBPs have been identified in wildlife [8,9], ticks [9][10][11], and livestock, including small ruminants kept by nomads [12]. Coexistence with wildlife can increase the risk of TBPs in livestock and man [13]. The most important tick-borne disease in sub-Saharan Africa, the corridor disease caused by Theileria parva, is believed to have co-evolved with the African cape buffalo, a host immune to the pathogen, before "jumping hosts" to cattle [14]. More recently, studies have reported diverse TBPs in Kenya, including a novel Rickettsia sp. identified in ticks in the coastal Shimba Hills National Reserve [5]. A prevalence of up to 5.5% spotted fever group rickettsioses was found in wildlife in Laikipia and the Maasai Mara National Reserve [8,9]. Tortoises in Baringo County and monitor lizards in Homa Bay County have been implicated as reservoirs of Ehrlichia ruminantium and Ehrlichia canis, respectively [10]. Similarly, Bunyamwera and West Nile viruses were identified from ticks collected from livestock and wildlife in Ijara District, northern Kenya [15]. These findings demonstrate the importance of wildlife in the circulation of diseases.
In Kenya, national parks, reserves, and conservancies account for about 11% of the country's landmass [16]. This is the direct result of wildlife conservation initiatives during the 1970s that led to the conversion of large portions of arid and semi-arid ecosystems, then considered to have little agricultural potential, into protected areas [1]. Ruma National Park, which occupies a third of Lambwe Valley in western Kenya, is a wildlife protected area where significant changes in land use have occurred over the last two decades [17]. In the last half century, intensive multiagency tsetse fly control activities in Lambwe Valley resulted in a marked decline in the incidence of trypanosomiasis [18] and contributed to the local elimination of human African sleeping sickness. The effective interventions impacted land use so that cultivated land increased by 22% over the last three decades [17]. Human populations on the periphery of the protected area have increased significantly [18], leading to an equivalent increase in livestock [19]. Such trends often increase the interactions between livestock, humans, and wildlife and may create a complex disease environment that leads to spill-over of diseases from wildlife to livestock and humans, and vice versa [20]. Despite the Ruma wildlife-livestock interface being a unique epidemiological niche of concern, empirical data on the epidemiology of TBPs circulating in livestock in the valley is limited.
Recently, smallholder farmers who reside within Lambwe Valley complained of episodes of outbreaks in livestock, causing significant fatality. Such reports also identified the presence of ticks and clinical signs compatible with haemoparasite infection in cattle. We hypothesised that the cattle were infected with TBPs and that cattle populations at the wildlife-livestock interface were at increased risk of TBP infection. We conducted a molecular epidemiological study to investigate the prevalence of TBPs and identify putative risk factors associated with exposure in cattle reared under a traditional system in the Lambwe Valley ecosystem.

Study Area
The study was carried out in Lambwe Valley, in village clusters around the Ruma National Park (Figure 1). Lambwe Valley (latitude 0 • 38 35.52 S, longitude 34 • 16 48 E) is 1200-1600 m above sea level and is infested by tsetse flies, the major vector for African trypanosomiasis. The terrain consists of rolling grasslands with open woodland and thickets dominated by acacia and a variety of grass species. The soil is predominantly black cotton, and the climate is hot and humid with an average annual air temperature of 25 • C. Rainfall is bimodal, peaking in April-June and September-November. The annual dry and hot period is from January to March [21]. The community in Lambwe Valley practice subsistence farming, fishing, and animal husbandry. In the park, the main grazing and browsing wild ruminant populations consist of roan antelope (Hippotragus equinus langheldi), Jackson's hartebeest (Alcelaphus buselaphus Jacksonii), oribi (Ourebia ourebi), and Rothschild giraffes (Giraffa camelopardalis rothschildi). The indigenous zebu cattle breed is the dominant domestic species followed by small ruminants. Livestock graze in the open fields extending to the park's fence, thus creating an interface for interaction between wildlife and livestock.

Study Design and Sample Size Determination
We used a cross-sectional design with stratified one-stage cluster sampling. All zebu cattle within a 5 km radius from the park's fence were defined as the study population. Construction of a sampling frame was not possible due to lack of reliable cattle demographic data. Therefore, we defined clusters based on the coordinates of a map grid (3 km 2 ), which later was named after the nearest village ( Figure 1). Data were collected between December 2018 and February 2019. Indigenous zebu cattle that were at least one year old and managed in smallholder farms in the wildlife-livestock interface of the Lambwe Valley were considered eligible for the study. This age category is likely to interact with other herds or wildlife at watering points and/or during grazing and is thus amenable to TBP surveillance in the area. Sample sizes were determined based on one-stage cluster sampling using previously described methods [22][23][24].
The sample size to estimate the prevalence with a specified precision is given by: where N is the sample size, P is the prevalence, D is the design effect, SE is the standard error of an estimated proportion P, g is the average number of animal sampled per cluster, and c the number of clusters sampled [22]. The design effect was given by the formula, We restricted the number of sampled villages (clusters) to 20% (c = 6) of the total 35 clusters, which were selected randomly, but accounting for accessibility. The intra-cluster correlation coefficient (ICC) is a measure of homogeneity of clustered data. We used the ICC value of 0.04 for Anaplasma marginale from previous studies [25]. Considering the possibility to collect a maximum of 100 samples by a team of six people per day per cluster, D was calculated as 5 (Eqution (2)). Prevalence information for TBPs is limited in the study area. Therefore, the study assumed a prevalence (P) of 50% for the pathogens. With an expected prevalence of 50%, a cluster size of 6 and the desired precision of 5%, a total sample size of 500 animals was estimated. The study was announced to farmers for sensitisation through local radio, and we accounted for a 25% non-response rate. To consider the sample size for the non-response rate: The above equations gave an estimated sample size of 667, and the final dataset consisted of a total of 680 cattle.
Epidemiological data on potential risk factors were collected from each sampled household using a questionnaire that captured sex, herd size, herd husbandry practices, grazing system, commonly used acaricides, and the frequency of application.

Ethical Approval
This study strictly adhered to the experimental guidelines and procedures approved by the Institutional Animal Care and Use Committee at International Centre of Insect Physiology and Ecology (icipe) and Kenya's animal welfare laws under the Veterinary Surgeons and Veterinary Para-professionals Act, 2011 (Cap. 366). Blood samples were collected only after receiving informed verbal consent from cattle keepers. Verbal rather than written consent was adopted as the pragmatic approach because all farmers were neither able to read nor write. Blood samples were collected by an experienced veterinarian with the aim of minimising pain and discomfort.

Blood Sample Collection and Processing
Blood samples were collected from the marginal ear vein to assess the anaemia status of each animal by measuring packed cell volume (PCV) using the micro-haematocrit method. Briefly, capillary blood was drawn into ethylenediaminetetraacetic acid (EDTA) microhematocrit tubes and centrifuged at 13,000× g for 10 min. PCV levels were calculated using the Hawksley haematocrit reader. For molecular detection of pathogens, about 4 mL of blood was collected from the jugular vein of each cattle using sterile vacutainer needles and EDTA vacutainers. Each tube was identified by a unique animal's ID number. Blood in the EDTA vacutainers were transferred into uniquely-labelled cryovials, stored in liquid nitrogen, and transported to icipe's Martin Lüscher-Emerging Infectious Diseases (ML-EID) Laboratory for molecular analysis.

PCR and High-Resolution Melting (HRM) Analysis
DNA was extracted from blood samples using the Bioline Isolate II genomic DNA kit (Meridian Life Science, Memphis, TN, USA) as described by the manufacturer and screened for the presence of Anaplasma, Theileria, Ehrlichia, Rickettsia, and Babesia spp. by PCR using single-plex genus-specific primers, as shown in Table 1. The 10 µL PCR volumes consisted of 2 µL of the extracted DNA template, 2 µL 5× HOT FIREPol ® EvaGreen HRM Mix (Solis BioDyne, Tartu, Estonia), 0.5 µL each of the 10 µM forward and reverse primers for the respective genus-specific reactions, and 5 µL PCR water. Each assay on different pathogens (Ehrlichia spp., Rickettsia spp., Anaplasma spp., Theileria spp., and Babesia spp.) had a set of known positive controls and a negative control (nuclease-free water) All PCR-HRM reactions were done in the Rotor-Gene Q 5-Plex HRM capable thermocycler (Qiagen, Hannover, Germany). The amplifications involved touchdown PCR conditions including an initial activation of the polymerase enzyme at 95 • C for 15 min, followed by 35 cycles of denaturation at a temperature of 9 • C for 20 sec. The 25 sec per cycle annealing temperature was reduced by 1 • C each cycle during the first ten cycles from 65 • C to 55 • C, and followed by extension at 72 • C for 25 s per cycle and a final extension step for 3 min at 72 • C. Then, the products were maintained for 3 min at 72 • C before proceeding with HRM analysis, in which the temperature was gradually increased from 75 • C to 95 • C with fluorescence acquisition after 2 s at 0.1 • C increments. Changes in fluorescence with change in temperature (dF/dT) were recorded and plotted. Melt curves were visualised and analysed on Rotor-Gene Q software version 2.3.1 (Qiagen, Hannover, Germany). The graphs were normalised between 100 and 0% fluorescence. The positive PCR samples were identified by comparing the different melting profiles with the positive controls using the HRM analysis software version 2.1.0 (Qiagen, Hannover, Germany). The PCR products were run in 2% agarose gels and purified using ExoSap-IT (USB Corporation, Cleveland, OH, USA) as per the manufacturers' instruction and then Sanger sequenced at Macrogen Inc. (Seoul, Korea). Gene sequences produced from this study were deposited in the GenBank database of the National Centre for Biotechnology Information. To amplify larger DNA fragments of Anaplasma, we performed nested PCR targeting family Anaplasmataceae 16S rRNA genes [27] in samples with Anaplasma amplicons that had unique HRM profiles. Amplification was done in a TM100 thermal cycler (BioRad, Hercules, CA, USA) in 10 µL PCR volumes comprised of 2 µL of 5X HOT FIREPol ® Blend Master Mix (Solis BioDyne, Tartu, Estonia), 0.5 µL each of the 10 µM forward and reverse primers, 2 µL template, and 5 µL PCR-grade water. Primary PCR amplifications using the EHR16SD and pH1522 primers (Table 1)

Phylogenetic Analysis
All sequences were edited and aligned, using the MAFFT [29] plugin in Geneious Prime version 2019.0.4 software (created by Biomatters, Auckland, New Zealand). They were subsequently queried in the GenBank nr database (http://www.ncbi.nlm.nih.gov/) using the Basic Local Alignment Search Tool (BLAST) [30]. We constructed maximum likelihood phylogenetic relationships of the study's gene sequences to those in GenBank using PhyML version 3.0. [31], employing the Akaike information criterion for automatic model selection. Tree topologies were estimated using nearest neighbor interchange improvements over 1000 bootstrap replicates. Trees were visualised in FIGTREE 1.4.2 [32].

Data Management and Analysis
A database consisting of questionnaire and molecular data was established in MS Excel (Microsoft ® Excel, Redmond, WA, USA). Individual-and herd-level information, including animal sex, age, animal weight, anaemia status, herd size, village name, distance of homestead from the park's fence, acaricide use (type and frequency), and infection status, were recorded. The distances of homesteads from the park's fence were measured using a handheld GPS device. Cattle ages were categorised as ≤2 and >2 years age groups and small (<10 head of cattle) or medium (>10 head of cattle) herd sizes of smallholder farmers. Anaemia status were categorized into anaemic (PCV ≤ 24) and normal (PCV > 24). After checking and variable coding, data were transferred to SPSS version 25 (IBM Corp, Armonk, NY, USA) for statistical analysis the individual animal and herd levels.
The prevalence of each TBP was calculated by dividing the number of positive samples by the total number of animals sampled during study period. Herd level prevalences were obtained by dividing the number of infected herds by the total number of herds tested. A herd was declared as positive for a TBP if one or more animals within the herd tested positive by PCR-HRM analysis.
Associations between prevalence and putative risk factors of TBPs were analysed using univariable logistic regression models. Risk factors associated with infection with specific TBPs were identified by using a multivariable logistic regression model and the strength of their association was assessed using adjusted odds ratios (OR). Variables with a p-value ≤ 0.25 in univariable analysis were included in a multivariable logistic model. A backward elimination procedure was used for further selection of variables. The variables were tested for interaction effect using cross-product terms and for multiple-collinearity using the collinearity matrix index before building the final model. The model validity and predictive ability were assessed using the Hosmer-Lemeshow test and receiver operating characteristic (ROC) curve. The confidence level was set at 95% with α = 0.05.
The datasets supporting the conclusion of this article are included within the article and Supplementary Materials files. Nucleotide sequences reported in this article are available via GenBank.

Pathogen Diversity and Prevalence
Blood from 680 cattle (370 females and 310 males), mainly zebu (Bos indicus), from 95 herds were screened for presence of TBPs. The average herd size was seven (ranging from 1-21 cattle per herd) and mainly open grazed as their major form of feeding.

Risk Factors Associated with Anaplasma and Theileria Infections
Results of the univariable logistic regression are summarised in Tables S2 and S3 and Figure 1. Each of the identified TBPs (n = 6) was used as an outcome in independent logistic regression analyses. Our study considered village, herd size, sex, age and distance from park as the risk factors. Small and medium herds were the herd size found in Lambwe Valley, while the distance considered was ≤1 km from the park and >1 km from the park. During the univariate analysis, the first level of each independent variable was used as a reference category. The results indicate that herd size, animal sex, and distance from the park were not significantly associated with the presence of any of the Anaplasma and Theileria species. However, there was a significant association between the identified Anaplasma and Theileria species and geospatially clustered villages in which the odds of infection in Kamato, Odelo, and Ruma Pap were significantly higher (Tables S2 and S3).
With the exception of A. marginale, all the Anaplasma species were detected in all the study villages. The presence of A. marginale was confirmed in only two villages, Kamato and Ruma Pap. The odds of A. platys clade, A. bovis, and Anaplasma sp. Lambwe-1 infection were highest in Kamato, Odelo, and Ruma pap (Table S2). The age of the animal was found to be significantly associated with Theileria infection; adults had higher odds of infection (OR = 1.8, p = 0.014) compared to calves and weaners. In contrast, the odds of T. velifera infection were lower in adult cattle (OR = 0.5, p = 0.005) compared to calves and weaners.

Discussion
This study provides insight into the diversity and complex co-occurrence of TBPs in cattle from smallholder farms at the wildlife-livestock interface of Ruma National Park, a protected wildlife area in Kenya. The findings show that cattle are infected with at least four Anaplasma species and two Theileria species. They further establish that known pathogens such as A. marginale, A. bovis, T. mutans, and T. velifera are prevalent in cattle from the study area. Notably, the study provides molecular evidence of A. platys and an Anaplasma spp. with unknown pathogenicity to livestock or humans.
The overall high prevalence of TBPs found in this study mirrors those recorded in Kenya and other countries in the region within the last five years. We found a TBP prevalence of 78.5% around Ruma. Regionally, reported TBP prevalence in cattle have ranged from 64.5% in Tanzania [33], 75.6% in central Uganda [34] and 89.6% in western Kenya [12], to 96.9% in Ethiopia [35]. The high prevalence of TBPs in cattle around Ruma and the region shows that TBPs could be a significant constraint to animal productivity in the majority of herds.
Three pathogens reported here, A. bovis, A. platys, and A. marginale, are known causative agents of bovine anaplasmosis. The clinical manifestation of infections with these pathogens include anaemia, fever, abortion, weight loss, lymphadenopathy, and death [7]. These symptoms are consistent with those described by local farmers from the study area and suggest that these TBPs could have contributed to the outbreak experienced in the area before the study.
This study reports infections of cattle with A. bovis and A. platys. Of these, A. bovis, was most prevalent (17.4%) and present in all sampled villages. However, this is lower than previously reported in calves (39.9%) in western Kenya [12]. Several tick species, including Rhipicephalus appendiculatus, Amblyomma variegatum, Hyalomma sp., Rhipicephalus sanguineus, and Haemaphysalis spp. transmit A. bovis [36]. Indeed, a study in western Baringo, Kenya, found that about 3% of the Amblyomma ticks collected from livestock had A. bovis [37]. Anaplasma platys commonly infects dog platelets and also causes illness in humans [38]. Matei et al. [39] confirmed that dogs roaming rural communities in islands off the Kenyan coast were infected with A. platys. This study, however, found A. platys-like DNA in 16.9% of cattle in the study area, a prevalence notably higher than reported in cattle in China (4.35%) [40], Tunisia (3.5%) [41], Brazil (4.75%) [42], and Kenya [43]. However, as there was considerable sequence variation among A. platys clade sequences, our findings may represent an unknown diversity of A. platys-like pathogens circulating in the region as found in Tunisia [41].
This study recorded a very low prevalence of A. marginale (0.6%). In a cross-sectional study conducted in two semi-intensively managed dairy farms in Machakos and Ngong districts, Adjou Moumouni et al. [44] reported A. marginale prevalence of up to 32.5%. Similarly, higher A. marginale prevalences were reported in Uganda (57%) [45] and Ethiopia (14.5%) [35]. Cattle breeds are differently susceptible to bovine anaplasmosis, with crossbreeds known to have relatively higher prevalences than indigenous cattle breeds [35,45]. The animals in this study were all local breeds.
Although Omondi et al. [37] reported diverse Anaplasma, Ehrlichia, and Rickettsia spp. in Amblyomma and Rhipicephalus ticks sampled from livestock in adjacent areas of Homa Bay County, A. marginale was not identified in the ticks. This complements our data to suggest that A. marginale is not prevalent in the region. Similarly, no A. marginale, Ehrlichia spp., nor Babesia spp. were identified in ticks sampled from another wildlife-livestock interface, the Maasai Mara National Reserve [11]. However, A. marginale, as well as T. parva, R. africae, and Babesia bigemina have recently been reported in ticks sampled at livestock markets in Busia County [46], approximately 175 km north of Homa Bay County, which suggests that there may be significant regional variation in TBP diversity associated with disparate agro-ecological contexts.
In 11.6% of the cattle sampled, we also identified a distinct Anaplasma sp. Lambwe-1 16S rRNA gene sequence that is identical to a recently identified Anaplasma sp. (clones "Saso" and "Hadesa") in Ethiopia [35]. The pathogenicity of this novel Anaplasma in cattle and other species and its zoonotic potential remains unknown. Since all of the sampled cattle in this study were apparently healthy, the present study did not evaluate the effect of the identified Anaplasma sp. However, the presence of Anaplasma sp. Lambwe-1 in the blood was significantly correlated with reduced haematocrite PCV.
We found a higher infection rate with Theileria spp. than with Anaplasma spp., with T. mutans at 25.7% and T. velifera at 40%. These pathogens, transmitted by several Amblyomma tick species, are less pathogenic and commonly found in cattle [5,47]. Our findings are consistent with most from the horn of Africa and southern Africa, which found higher prevalences of T. mutans and T. velifera in cattle compared to other TBPs [12,48]. Infections with T. velifera in cattle are mostly asymptomatic with animals infected at younger ages having lifelong infections [35,48,49]. Herd prevalences of Anaplasma spp. (78.9%) and Theileria spp. (88.4%) were higher than their respective individual prevalences of 45.7% and 63.5%. Higher herd prevalence may indicate that the pathogens are more broadly distributed than can be assessed by individual prevalence alone. High herd infection prevalence poses a concern on the herd health of cattle in Lambwe Valley.
Our study provides the first detailed analysis of co-infections in bovine species in Lambwe Valley. The high prevalence of co-infection was mainly due to double infection, the highest rates of which were for A. bovis and T. velifera (33%), and for Anaplasma sp. Lambwe-1 and T. mutans (31%) co-infections. The overall prevalence of TBP co-infection was, however, lower than previous studies in western Kenya [12]. Njiiri et al. [12] showed that 87.1% of cattle had mixed infections and that T. velifera, Theileria sp. clone "sable", and T. mutans were more prominent. In Ethiopia, Hailemariam et al. [35] reported a higher rate (57.9%) of relatively more diverse co-infection ranging from double to sextuple pathogen coinfections. Co-infections can augment transmission of pathogens and/or promote the severity of the diseases [50]. In co-infected animals, the clinical symptoms may not be similar to those in hosts with single infections. Woolhouse et al. [51] demonstrated that mild Theileria species reduced the severity of East Coast fever infection in Kenya. This study does not investigate the implication of co-infections to the burden of disease in Lambwe Valley. It however confirms that co-infections do occur in the wildlife-livestock interface and could lead to misdiagnosis and mismanagement of livestock diseases. In addition to co-infections, we found Wolbachia, an obligate endobacterial symbiont of bovine Onchocerca species that serves as a source of energy for filarial nematodes [52]. Its detection in bovine blood suggests an occult infection of filarial nematodes.
Our risk factor analysis showed that livestock from different villages were at different risks of infection with TBPs. No significant risks were associated with the sex of the animals, the sizes of herds they came from, or the distances of respective homesteads from Ruma National Park. Three village clusters north-east of the park had Anaplasma species prevalences > 75%, while three in the west and south-western side had a combined prevalence of 21% for the same pathogens. The type of grazing, tick control interventions applied in different villages, and the composition and abundance of vectors for the pathogens could explain disparities in prevalence between villages. Unlike this study, a previous study found higher prevalence of A. marginale in smaller herd sizes compared to medium and larger herds in Peninsular Malaysia [53]. However, the Lambwe Valley communities practice open grazing and cattle share communal watering points. We had hypothesised that cattle populations nearer to the park share more space and potential contact with wildlife than cattle populations further away from the park. However, our analysis showed that the odds of TBPs infection in cattle that reside in close proximity to the park (< 1 km) are not different from those that reside further away. The extensive grazing and livestock management system in Lambwe Valley means that cattle move and the approach of using fixed distances for herds may be inappropriate. This study did not sample wildlife to determine their exposure status with regard to the targeted pathogens; we were therefore unable to confirm a role of wildlife on the observed molecular prevalences.
In conclusion, this study provides an update on the epidemiological status of TBPs circulating in cattle of Lambwe Valley, Kenya. It highlights the high prevalence of infections in the wildlife-livestock interface and reports the existence of highly pathogenic A. marginale and A. bovis, as well as other novel Anaplasma spp., T. mutans, and T. velifera, with considerable co-infection rates. The results presented offer neither proof of inter-species pathogen transmission nor the incrimination of vectors responsible for the transmission of different pathogens. However, the findings add emphasis on the need to proactively investigate the prevalence and composition of TBPs in wildlife-livestock interfaces.
Supplementary Materials: The following are available online at http://www.mdpi.com/2076-2607/8/11/1830/s1, Summary result of the univariable logistic analysis of risk factors with dependent genus Anaplasma and genus Theileria infection (Table S1), Anaplasma species (Table S2) and Theileria species (Table S3)   The funders had no role in study design, data collection and analysis, decision to design, data collection and analysis, decision to publish, or preparation of the manuscript. The views expressed in this publication are those of the authors and not necessarily those of any of the funding agencies.