Frequency and factors associated of potential zoonotic pathogens (Borrelia spp., Rickettsia spp., Leishmania spp., and Anaplasma phagocytophilum) in equids in the state of Bahia, Brazil
Parasites & Vectors volume 14, Article number: 275 (2021)
Currently, various zoonotic diseases are classified as emerging or reemerging. Because equids have a direct relationship with various vectors, they are possibly more frequently exposed to zoonotic agents than are humans. The undeniable importance of diseases such as human granulocytic anaplasmosis, spotted fever, and leishmaniasis for both public and animal health, as well as the possibility of equids acting as sources, reservoirs, or even sentinels for these pathogens, justifies the detection of their frequency and factors associated with infection in equids from northeastern Brazil.
Blood samples were collected from 569 equids (528 horses, 33 donkeys, and 8 mules), 516 from a rural area and 53 from an urban area. Pathogen detection was carried out as follows: Borrelia spp. and Rickettsia spp., serological analysis; Leishmania spp., serological analysis and polymerase chain reaction (PCR); Anaplasma phagocytophilum, PCR. Determination of associated factors was carried out through generalized linear models.
The frequencies of positivity for the pathogens observed in equids were as follows: Borrelia spp., 13.9% (79/569); Leishmania spp., 3.5% (20/569); Rickettsia spp. 33.4% (190/569). Regarding factors associated with infection, male sex was associated with protection against Borrelia spp.; donkeys and mules were associated with protection against Rickettsia spp., while a younger age was a risk factor. The infection of A. phagocytophilum was not detected in the sampled population. Co-infection was detected in 5.1% (29/569) of the animals.
Most of the studied pathogenic agents are present in the prospected area, indicating a possible risk for both human and animal health. This demonstrates that equids can be considered important sentinels in the assessment of pathogens with zoonotic potential in the region.
Equids can be infected with various pathogens that also infect humans [1,2,3]. The close contact between these two species, which spans for over 5000 years, potentially allowed some agents to develop the capacity to develop and/or multiply in both of them. In studies on populations, concepts like ecological fitting and fitness space, which explain the capacity for interactions among microbiota, flora, fauna, climate, and environment in the search for equilibrium, also seem to explain this phenomenon .
Currently, various zoonotic diseases are classified as emerging or reemerging . Climate change, increased possibility of dispersion, and interactions between vertebrate and invertebrate hosts, seem to be associated with the appearance, or even outbreaks, of certain diseases in different regions of the planet [4, 6, 7].
Because equids have a direct relationship with various vectors, which is associated with the environment and management practices they are exposed to, they can be more frequently exposed to hemoparasites and other zoonotic agents than are humans [8,9,10]. As such, they are able to play an important role as sentinels for certain diseases [11,12,13,14,15].
Diseases such as human granulocytic anaplasmosis, spotted fever, and leishmaniasis have an impact on public health. This causes significant direct and indirect economic losses, morbidity, and, in some cases, death in the human population. Similar to leishmaniasis , these diseases are highly prevalent in certain regions of the world. The causative agents of these illnesses also infect equids and invertebrate hosts present in the same environment, potentially increasing their exposure to infection. Horses are normally carriers, but clinical signs of borreliosis (Lyme disease) [17, 18], granulocytic anaplasmosis [19, 20], leishmaniasis [21, 22], and rickettsiosis [23, 24] have been described in previous reports.
The undeniable importance of the above-mentioned diseases for public and animal health, as well as the possibility of equids acting as sources, reservoirs, or even sentinels for these pathogens, justifies the detection of the frequency and factors associated with infections in equids from northeastern Brazil.
Study area and sample population
Data were collected between August 2013 and December 2014 in the microregion of Ilhéus-Itabuna, in the state of Bahia, northeastern Brazil. This geographic region is part of the mesoregion of southern Bahia and has an estimated equid population of 90,974 animals . The study area is located in the Atlantic Forest. The annual average rainfall is 1445 mm, with a relative humidity of 80% and a temperature of 24 °C . Five counties from this mainly rural microregion were selected for the study and were ranked according to the size of their equid population as follows: Itaju do Colônia (15° 08′ S 39° 43′ O), Itapé (14° 52′ S 39° 25′ O), Ibicaraí (14° 51′ S 39° 35′ O), Santa Cruz da Vitória (14º57'S 39° 48′ O), and Floresta Azul (14° 50′ S 39° 39′ O). The county of Itabuna (15° 8′ S 39° 43′ W), which is a mostly urban area, was also included in the study (Additional file 1: Figure S1).
Animals, farms, and counties were selected based on their convenience. The number of animals per county was proportional to their equid population. Blood samples were collected from 569 equids (528 horses, 33 donkeys, and 8 mules); 516 out of 569 equids were from 20 rural properties; 53 horses out of 569 equids were from urban areas and were used by mounted police, to draw coaches, or for horseback riding.
To evaluate potential risk factors, information regarding signalment (species, age, sex), farm characteristics, and management (animal kept in a stable, presence of ticks, contact with other animal species) was obtained through semi-structured interviews with handlers (staff) or owners. Interviews were always conducted by the same researcher.
The study was carried out according to the standards established by the Brazilian College of Ethics and Animal Welfare. The research proposal was approved by the Committee for Ethics in Research with Animals (protocol 002/2013) at the State University of Santa Cruz, Ilhéus, BA, Brazil.
Sample collection and processing
Blood (20 ml) was collected from each animal through jugular venipuncture using disposable needles (25 × 8 mm) connected to vacuum tubes with and without anticoagulant (EDTA). Tubes with anticoagulant were centrifuged for 10 min at 699×g. The supernatant plasma was discarded and then both the leukocyte layer and the packed red blood cells were removed, poured into DNase- and RNase-free plastic tubes, and frozen at − 20 °C for subsequent DNA extraction. To obtain the sera, the tubes without anticoagulant were centrifuged at 699×g for 10 min; the sera separated through aspiration were placed in plastic tubes and frozen at − 20 °C until the serological tests were performed.
Ticks were collected from equids during visits to farms. Taxonomic identification of the ticks was performed using a stereoscopic microscope, according to previously published morphological keys [27,28,29]. The ticks were stored in micro-tubes or Falcon tubes, depending on the quantity, and labeled with the name of the property from which they were collected. Subsequently, they were submerged in absolute alcohol and stored at − 20 °C for future morphological identification and DNA extraction.
For the detection of antibodies against the crude B. burgdorferi (s.l.) antigen G39/40 strain, indirect ELISA was performed using the antigen at a concentration of 15 µg/ml, serum in phosphate buffer saline (PBS Tween 20, 0.05%, pH 7.4) at a dilution of 1:800, and conjugated antibodies at a dilution of 1:5000 (anti-horse IgG, A6063, Sigma®) . The cut-off point for the assay was determined using the previously described methodology , which is 2.5 times the mean value of negative control (animal previously tested) absorbance. Optical density was measured at 405 nm. Four negative samples from previously tested animals were used as negative controls, and there were two positive controls from animals inoculated with crude B. burgdorferi antigen, G39/40 strain.
For the detection of antibodies against the soluble L. braziliensis antigen, indirect ELISA was performed using the antigen at a concentration of 10 ng/ml, serum in phosphate buffer saline (PBS Tween 80, 0.05%, pH 7.4) at a dilution of 1:200, and conjugated antibodies at a dilution of 1:30,000 (anti-horse IgG, A6063, Sigma®) . The cut-off point for the assay was determined using the previously described methodology , which is 2.5 times the mean value of negative control (animal previously tested) absorbance. Optical density was measured at 405 nm.
Rickettsia rickettsii and R. parkeri
The sera were tested using an indirect fluorescent antibody test (IFAT); a dilution of 1:64  was the cut-off point for antigens of R. rickettsii Taiaçu strain  and R. parkeri At24 strain . The slides had been previously sensitized with an antigen produced through the cultivation of R. rickettsii and R. parkeri in Vero cells. Conjugated Anti-Horse IgGs (F7759, Sigma-Aldrich ®) were used at a dilution of 1:80 . A microscope with an epifluorescence system (OLYMPUS, BX 51) was used for reading the slides. Reactions with complete fluorescence in the periphery of the agents were considered positive. Positive and negative controls were obtained from the Molecular Epidemiology Laboratory of the Fluminense Federal University (UFF). The positive samples were two-fold titrated.
Extraction of DNA from equine blood samples and molecular diagnosis
DNA was extracted from blood samples using a commercial QIAamp® DNA Blood Mini Kit (QIAGEN™), according to the manufacturer's recommendations. DNA samples were labeled with accession numbers and stored in a freezer at − 20 °C for subsequent polymerase chain reaction (PCR).
Specific RV1/RV2 primers (Table 1) for L. infantum  and B1/B2 primers (Table 1) for L. braziliensis were used . PCR was carried out using 10 × Taq DNA polymerase buffer, 0.2 mM of MgCl2, 2 mM of each dNTP, 10 pmol of each L. infantum primer, 1.25 U of Taq DNA polymerase (Invitrogen®), and 5 µl of genomic DNA. Similar conditions were adopted for L. braziliensis using 3 mM of MgCl2 and 2.5 U of Taq DNA polymerase. The thermocyclic conditions were as follows: initial denaturation, 94 °C for 5 min; 35 cycles of denaturation, 94 °C for 1 min, annealing of primers, 59 °C (L. infantum) or 65 ºC (L. braziliensis) for 1 min and extension, 72 ºC for 1 min; final extension, 72 °C for 5 min. Ultrapure water was used as a negative control. Positive controls of the Huhu strain were provided by the Gonçalves Muniz Research Center (FIOCRUZ-BA).
Reactions were carried out using primers described in a previous study  to increase the region of the 16S rRNA gene, gE3a/gE10R in the first reaction (932pb), and gE2/gE9f in the second reaction (546pb) (Table 1), using a final volume of 12.5 µl containing 2.5 µl of genomic DNA, 10 × reaction buffer, 1.5 mM of MgCl2, 0.2 mM of each dNTP, 0.4 µM of each primer, and 1.25 U of Taq polymerase; ultrapure water was added until the final volume was obtained. The thermocyclic conditions were as follows: initial denaturation, 94 °C for 5 min; 40 cycles of denaturation, 94 °C for 30 s, annealing, 55 °C for 1 min, extension, 72 °C for 1 min; and final extension, 72 °C for 5 min . For the nested-PCR, the same concentrations and final volume were used for the “mix,” that is, 0.5 µl of the first reaction’s product was used. The number of cycles was reduced to 30 while maintaining the thermocycler’s time and temperature conditions. Ultrapure water was used as a negative control.
PCR for GAPDH (glyceraldehyde-3-phosphate dehydrogenase) detection
To verify DNA integrity and the presence of potential inhibitors, negative samples were subjected to PCR for detection of the GAPDH gene using primers described by Birkenheuer et al.  (Table 1). In the PCR reactions, a final volume of 25 μl was used, composed of 5 μl of genomic DNA, 10 × reaction buffer, 2.0 mM of MgCl2, 0.2 mM of each dNTP, 0.4 μM of each primer, 1.25 U of Taq polymerase, and ultrapure water until the final volume was reached. The amplification protocol used consisted of an initial denaturation step at 95 °C for 5 min, followed by 40 cycles at 94 °C for 30 s for denaturation, annealing at 52 °C for 1 min, extension at 72 °C for 1 min, and final extension at 72 °C for 5 min .
Extraction of genomic DNA from the ticks and PCR for Rickettsia spp. detection
Total DNA was extracted from the ticks individually, using the phenol–chloroform method .
All the samples were subjected to a reaction using specific oligonucleotides for the 16S gene (Table 1) to verify DNA integrity and the possible presence of PCR inhibitors. The thermocyclical conditions used were as follows: initial denaturation, 94 °C for 2 min; 35 cycles of denaturation, 94 °C for 30 s, annealing, 55 °C for 30 s, extension, 72 °C for 45 s; the final extension, 72 °C for 7 min .
To determine the presence of Rickettsia spp. DNA, the pair of primers CS-239 and CS-1069 was used, which amplifies 834 pb of the gltA gene (Table 1). This gene is present in all known Rickettsia species. The thermocyclical conditions used were as follows: initial denaturation, 95 °C for 5 min; 40 cycles of denaturation, 95 °C for 20 s, annealing, 52 °C for 20 s, extension, 72 °C for 40 s; final extension, 72 °C for 5 min .
PCR products were detected using 2% agarose gel electrophoresis in a tris-acetate-EDTA (TAE) running buffer (40 mM Tris-acetate, 2 mM EDTA pH 8.0). The gel was run at 80 V, 180 mA for 30 min, and then stained with ethidium bromide (0.5 μg/ml). A DNA molecular weight standard control (1 Kb Plus DNA Ladder, Invitrogen®) was used to estimate the size of the amplified products. Amplified products were visualized under ultraviolet (UV) transilluminator (LPIX, Loccus Biotecnologia®) and photographed on a coupled image analyzer.
The variables were categorized for the purposes of statistical modeling as follows: species (horse or donkey + mule); age, in the form of age ranges, i.e. young (≤ 3 years), adult (> 3 and < 12 years), or senior (≥ 12 years); sex (male or female); equids kept in a stable (yes or no); equids who had contact with sheep (yes or no); equids who had contact with goats (yes or no); equids who had contact with poultry (yes or no); equids who had contact with cattle (yes or no); presence of rats on the farm (yes or no); presence of toxic plants in the pastures (yes or no); and inserted into the models according to the biological plausibility for each agent. The presence of Borrelia spp., Rickettsia spp., Leishmania spp. and Anaplasma phagocytophilum infection (yes/no) was considered as outcome variables. Donkeys and mules were inserted into a single category because of the low number of animals.
Blood samples from animals from urban areas were collected only in the county of Itabuna; a total of 53 horses were sampled, including those used by the mounted police, to draw coaches, or for horseback riding. These animals were excluded from modeling because of the different management practices applied to these horses. They were only used to compare frequency distributions of agents between animals living in rural and urban areas.
Generalized linear models with binomial distribution were used to perform bivariable and multivariable analyses. Because of the possibility of clusters, intraclass correlation coefficients (ICCs) were calculated. Null models  were estimated to obtain ICCs for municipality and farm variables because of the possibility of observations of animals from the same municipality and/or farm being correlated, forming clusters . In cases of cluster formation, the variables tested were considered random, and generalized linear mixed models were used for the analyses.
The modeling strategy used in the multivariable analyses was backward, that is, all variables were initially included in the model. Starting from this initial model, variables were selected at each step based on the Wald test until the most parsimonious model that best explained the outcome was obtained. The significance level for variables to remain in the final model was set at 5%. The Akaike information criterion (AIC) was used to evaluate the fit of the models. Frequency distributions among animals living in rural areas and animals living in urban areas were compared using the chi-square test.
Odds ratios (OR) and their respective 95% confidence intervals (CI) were calculated based on the regression coefficients that had been estimated through the models. Statistical calculations were performed using R software, version 3.2.5, for Windows  through the lmer4 package , version 1.1-12.
Regarding Rickettsia spp., 33.39% (190/569) of the equids showed positivity for at least one of the tested antigens. The proportion of animals that reacted serologically to only R. rickettsii antigens was 14.6% (83/569)—14.9% (77/516), 11.3% (6/53) of which belonged to the rural and urban populations, respectively (p > 0.05). The antibody titers varied from 1:64 to 1:1024 (Additional file 1: Table S1). Serologically reactive equids were found in all the studied municipalities, with the proportion of positivity varying from 6.6–44.7%. Among the 20 rural properties in the state, animals in five properties did not show positivity, while the proportion of positivity varied from 4.0–50.0% in the others. Positivity among the equids was distributed as follows: horses, 15.3% (81/528); donkeys and mules, 4.8% (2/41).
Regarding R. parkeri antigens, 15.1% (86/569) of the animals had antibodies against only these antigens—15.5% (80/516) and 11.3% (6/53) of which belong to the rural and urban populations, respectively (p > 0.05). The antibody titers varied from 1:64 to 1:512 (Additional file 1: Table S1). Positivity of equids was observed in all of the studied municipalities, with proportions varying from 5.6 to 33.7% (Additional file 1: Table S2). Among the 20 rural properties, animals in three properties did not show positivity. On the properties showing reactivity, the proportion of positivity varied from 3.3 to 31.5%. Among the evaluated species, positivity was distributed as follows: horses, 15.3% (81/528); donkeys and mules, 12.2% (5/41).
In 3.7% (21/569) of the animals, responses to both R. parkeri and R. rickettsii antigens were noted. Such animals constituted the “double-reaction” group. Among the animals presenting double reactivity, 66.6% (14/21) and 9.5% (2/21) had titers of 1:64 and 1:128, respectively, for both agents; however, five animals had divergent titers among the species.
The ICCs (intraclass correlation coefficients) for the municipalities (1.26%) and farms (2.96%) did not indicate the formation of clusters. Generalized linear models were therefore used (Additional file 1: Tables S3 and S4). The final model (Table 2) shows other species (donkey + mule) associated with protection and the young age group as a risk factor for contracting the infection.
The serological analysis for Borrelia spp. showed that 13.9% (79/569) of the animals tested positive, with 13.9% (72/516) and 13.2% (7/53) in the rural and urban zones (p > 0.05). Animals that tested positive were identified in all of the municipalities. On 85% (17/20) of the evaluated farms, at least one positive equid was found. Besides, 13.4% (71/528) of the horses and 19.5% (8/41) of the donkeys and mules tested positive for Borrelia spp. The ICCs for the counties (1.5%) and farms (4.5%) indicated clusters only at the farm level. Therefore, mixed regression models with a farm as the random effect were used for the analysis of possible risk factors for the occurrence of Borrelia spp. infection (Additional file 1: Table S5, S6) in the final model; this analysis indicated that male sex [OR 0.42(0.19–0.91); p = 0.03] was associated with protection against Borrelia spp.
For Leishmania spp., 3.5% (20/569) of animals tested seropositive in ELISA, all of which belonged to the rural zone; none of the animal tested positive as determined by PCR analysis. All the municipalities of the prospected rural area presented seropositive animals. At least one positive animal was observed on 40% (8/20) of the farms. Among the horses, 2.6% (14/528) tested positive for Leishmania spp., while this percentage was 14.6% (6/41) for donkeys and mules. Considering the low number of positive animals, which could generate erroneous results, performing a descriptive analysis of the results was considered appropriate. Thus, greater positivity was identified in donkeys and mules (Additional file 1: Table S7).
Distribution of co-infections among equids was as follows: Borrelia spp. and Rickettsia spp., 3.86% (22/569); Borrelia spp. and Leishmania spp., 0.7% (4/569); and Rickettsia spp. and Leishmania spp., 0.52% (3/569) (Fig. 1).
Tick control was reported at all collection locations; however, ticks were found at all the sampled properties. A total of 262 ticks were collected, which were identified as follows: Dermacentor nitens, 202; Amblyomma sculptum, 61; and Rhipicephalus (Boophilus) microplus, 23. Only the ticks of the A. sculptum species were subjected to PCR for the detection of Rickettsia spp.; however, all the reactions were negative.
Amplification of A. phagocytophilum DNA was not observed in the assessed animals. Amplification for GAPDH was noted in all the negative samples.
This study involved the largest number of equids so far for the detection of positivity for antigens of Rickettsia spp., Leishmania spp., and Anaplasma spp. in Brazil and the second largest sample for that of Borrelia spp. The study was conducted in a region characterized by a humid tropical climate, favoring the biological cycle of the vectors throughout the year. It is also in the Atlantic Forest biome, where there are various rodents of Brazilian fauna, such as the capybara (Hydrochaeris hydrochaeris), which is suspected to play an important role in the life cycles of Rickettsia spp.  and Borrelia spp. . Ticks were observed on all the sampled properties, justifying the ample spread of these pathogens.
The frequency of animals showing antigens related only to R. rickettsii (14.6%) or R. parkeri (15.1%) were higher in the areas included in this study than in a neighboring area (5.8% and 8.7%, respectively) . This fact can be explained by the small sample size when compared to our study, which covered the geographical area and a number of samples about ten times larger than those of the neighboring area. Ticks tested positive for R. parkeri in this region ; therefore, it is possible that animals with elevated anti-R. parkeri titers are in fact positive for this pathological agent.
Epidemiological surveys carried out on equids with respect to Rickettsia spp. rarely assessed risk factors. In the present study, however, the age range below 3 years was identified as a possible risk factor for infection by this pathological agent. This indicates early exposure of equids to the agent, with a possible decrease in circulating antibodies with time. Unfortunately, it was not possible to find data from other authors to promote discussion of the theme.
The results of the present study showed a lower frequency of Rickettsia spp. in donkeys and mules than in horses, which corroborates previous findings [50, 51]. This result may be associated with a greater resistance against Amblyomma spp. infestation in mules and donkeys than in horses ; however, positivity can be high in animals greatly challenged by vectors . Another interesting point was the absence of positivity for R. rickettsii antigens in donkeys. This may suggest that besides resistance to vectors, there may be differences in resistance against species of Rickettsia.
The absence of molecular detection of the gltA gene corroborates the findings of the state of Goiás  and those of Bahia , where not a single tick collected from horses tested positive for the gltA gene.
In the state of Bahia, between 2005 and 2016, there was at least one notification of infections caused by Rickettsia spp. per year in human beings, with no deaths . The importance of domestic animals as amplifying hosts of Rickettsia spp. is yet to be completely explained . These previously reported findings, together with the frequency observed in the present study, demonstrate that Rickettsia spp. is circulating in the state, considering that equids are indicated as sentinel animals by epidemiological studies and that they can be infected without presenting clinical symptomatology [56,57,58]. Furthermore, the positivity in animals of urban origin, as found in the present study, raises the risk of urbanization of the disease .
Borreliosis is a disease that is neglected in most of the world . The few studies on the disease in equids in Brazil indicate distinct seropositivity varying from 7.2 to 44.7% [3, 18, 60,61,62]. Even in studies carried out in the same state, large variations can occur; for example, among the studies carried out in the state of Pará, a seropositivity of 26.7%  was detected in one study while that of 7.2% was detected in another study . The frequency found in the present study was much higher than that found in dogs (1%) from the same microregion . One of the explanations may be the preference of A. sculptum (Amblyomma cajennense complex)  to parasitize equids, making these species better indicators of the presence of the pathological agent than dogs. Rhipicephalus microplus tick parasitism may be another important factor because of the possibility of cross-reaction with Borrelia theileri, a tick-borne spirochete known to infect cattle and other mammals like horses, sheep, and deer .
Similarities in the results of urban and rural animals may reflect the increase in the urban population of vectors (ticks) and reservoirs (capybaras) and demonstrate the risk of human exposition to the pathogen in the urban zone. No previous studies on this type of population were found.
Factors associated with Borrelia spp. infection were rarely cited [3, 8, 9, 18, 66], with contact and intensity of tick infestation [61, 67] and age  being the most common. In the present study, all the animals presented tick infestation; as such, this variable cannot be evaluated. However, the male sex was identified as a possible protective factor, which may be explained by the fact that most of the male equids/stallions receive differentiated management, which probably leads to lower infestation and exposure to transmitting agents. Following the logic that tick infestation intensity increases the risk of infection, researchers attribute a greater prevalence in a determined location and breed to the environment where they are found, which in reality is more conducive to maintenance of the vector .
As in the case of rickettsiosis, horses can be considered sentinels for borreliosis , as corroborated by the results of the present study. The detection of Borrelia spp. DNA in 43% of the ticks collected from horses  enables us to reach the conclusion that equids can be considered reservoirs of infection. Thus, the equids may be a multiplier and spreader of infected vectors to the peri-urban environment, as with Rickettsia spp. .
The low frequency of Leishmania spp. infection in the equids of the present study partly corroborates the data from the Ministry of Health’s National Disease Notification System, which, in 2015, registered only ten cases of human tegumentary leishmaniosis in the municipalities included in the study. There were no recorded cases of visceral leishmaniosis in the region , which has approximately 279,464 inhabitants .
It can be observed that despite its low frequency in the animals, this rate was much higher than the recorded cases of the disease in humans in the region. This result indicates that equids can also be considered sentinels for leishmaniosis , as well as sources of infection for the vectors [74, 75], because DNA of the parasite has already been detected in the blood of horses [15, 22].
Initially, leishmaniosis was described in rural areas because of the characteristics of its vectors, which maintain a sylvatic cycle . Currently, the disease can be identified in urban or peri-urban areas as a result of adaptation of the vector to this new environment and establishment of new reservoirs [15, 21, 22]. Unlike that of Borrelia spp. and Rickettsia spp., detection of Leishmania spp. in animals the present study occurred only in the rural environment. This may be indicative of absence or low prevalence of the vector or the agent from the studied urban environment or of the fact that equids are not good sentinels for this agent in areas with low prevalence of this pathological agent.
Despite the small number of animals that tested positive for Leishmania spp., which compromised the performance of statistical analysis, greater positivity was found in mules and donkeys than in horses. Researchers have identified donkeys as probable sources of infection introduction in a Venezuelan outbreak area, because the first skin lesions suggestive of the disease were reported in donkeys coming from the endemic region . Furthermore, it was observed that these animals were the main source of food for phlebotominae; when compared to the harassment of vectors between donkeys, humans, and dogs these animals were observed to be the main food sources for phlebotominae, unlike the vector hosting among donkeys, humans, and dogs . Nevertheless, epidemiological studies showing donkeys and mules as risk factors for leishmaniosis were not found in the literature.
The negative result of the PCR made it impossible to identify the species involved. However, given the absence of recorded cases of visceral leishmaniosis in humans and dogs in the region, it is probable that the evaluated equids are asymptomatic carriers of tegumentary leishmaniosis. This reinforces the importance of controlling the transportation of equids, especially donkeys, from endemic to non-endemic areas.
Despite having evaluated the largest sample of horses in Brazil, it was not possible to detect A. phagocytophilum DNA in the animals in the present study. The lack of detection, as well as the detection of low prevalences, as in the findings from eastern Europe (1.4% prevalence) , can be explained through the parasitemia of the pathological agent in horses, occurring for a short period (± 129 days) [78, 79]. This makes serology, the most commonly used method [20, 80], more sensitive than PCR. Performing serology would enable identification of chronic cases, confirmation of endemics of the region, and provision of conditions for evaluating the role of equids in the epidemiology of the disease; however, it was not possible to carry it out.
Few reports on co-infections among zoonotic organisms are available in the literature [81, 82]. The present study detected co-infections of Leishmania spp., Borrelia spp., and Rickettsia spp.; to our knowledge, this is the first report of this type of co-infection. It was not possible to determine how these interactions interfere with animal well-being and public health.
We were able to verify that infections caused by Leishmania spp., Borrelia spp., and Rickettsia spp. are present in the equine population of the studied area with positivity levels above those observed in other species. This may be a good indicator of the species as a sentinel for these infections. Despite the detected frequencies not being considered high, the participation of equids in the transmission cycle of these pathogens to their vectors or to the human population cannot be ruled out. The detection of animals that tested positive for Borrelia spp. and Rickettsia spp. in the evaluated urban zone is worth highlighting, which demonstrates the risk of these agents for public health in the region.
Availability of data and materials
http://www.biblioteca.uesc.br/biblioteca/bdtd/201370049T.pdf. Other data will be made available on request.
Butler CM, Houwers DJ, Jongejan F, van der Kolk JH. Borrelia burgdorferi infections with special reference to horses. A review. Vet Q. 2005;27:146–56.
Ishag MY, Majid AM, Abuobeida MM. Studies on Toxoplasmosis in horses in Khartoum State-Sudan. J Vet Med Anim Prod. 2014;5:29–35.
Galo KR, Fonseca AH, Madureira RC, Neto JDB. Frequência de anticorpos homólogos anti-Borrelia burgdorferi em equinos na mesorregião metropolitana de Belém. Estado do Pará Pesqui Vet Bras. 2009;29:229–32.
Brooks DR, Hoberg EP. How will global climate change affect parasite-host assemblages? Trends Parasitol. 2007;23:571–4.
Brasil. Leishmaniose visceral: recomendações clínicas para redução da letalidade. Ministério da Saúde. 2011;9:1–78.
Naucke TJ, Menn B, Massberg D, Lorentz S. Sandflies and leishmaniasis in Germany. Parasitol Res. 2008;103:65–8.
Rikihisa Y. Mechanisms of obligatory intracellular infection with Anaplasma phagocytophilum. Clin Microbiol Rev. 2011;24:469–89.
Hansen MGB, Christoffersen M, Thuesen LR, Petersen MR, Bojesen AM. Seroprevalence of Borrelia burgdorferi sensu lato and Anaplasma phagocytophilum in Danish horses. Acta Vet Scand. 2010;52:3. https://doi.org/10.1186/1751-0147-52-3.
Kiss T, Cadar D, Krupaci AF, Bordeanu A, Brudaşcă GF, Mihalca AD, et al. Serological reactivity to Borrelia burgdorferi sensu lato in dogs and horses from distinct areas in Romania. Vector Borne Zoonotic Dis [Internet]. 2011;11:1259–62. https://www.researchgate.net/profile/Andrei_Mihalca/publication/51165173_Serological_reactivity_to_Borrelia_burgdorferi_sensu_lato_in_dogs_and_horses_from_distinct_areas_in_Romania/links/0deec532317d0e0bf1000000.pdf
Maurizi L, Marie J, Aoun O, Courtin C, Gorsane S, Chal D, et al. Seroprevalence survey of Equine Lyme Borreliosis in France and in Sub-Saharan Africa. Vector-Borne Zoonotic Dis. 2010;10:7–10. https://doi.org/10.1089/vbz.2009.0083.
Brüning A. Equine piroplasmosis an update on diagnosis, treatment and prevention. Br Vet J. 1996;152:139–51.
Coiro CJ, Langoni H, da Silva RC. Epidemiological aspects in the Leptospira spp. and Toxoplasma gondii infection in horses from Botucatu, São Paulo Brazil. J Equine Vet Sci. 2012;32:620–3. https://doi.org/10.1016/j.jevs.2012.02.008.
Reuss SM. Review of equine cutaneous leishmaniasis : not just a foreign animal disease. Am Assoc Equine Pract. 2013;59:256–60.
Ribeiro AJ, Cardoso L, Maia JM, Coutinho T, Cotovio M. Prevalence of Theileria equi, Babesia caballi, and Anaplasma phagocytophilum in horses from the north of Portugal. Parasitol Res. 2013;112:2611–7.
Truppel JH, Otomura F, Teodoro U, Massafera R, Da Costa-Ribeiro MCV, Catarino CM, et al. Can equids be a reservoir of Leishmania braziliensis in endemic areas? PLoS ONE. 2014;9:1–6.
Brasil. Manual de Vigilância e Controle da Leishmaniose Visceral. Ministério da Saúde. 2006;1a edição:120.
DeVilbiss BA, Mohammed HO, Divers TJ. Perception of equine practitioners regarding the occurrence of selected equine neurologic diseases in the Northeast over a 10-year period. J Equine Vet Sci. 2009;29:237–46. https://doi.org/10.1016/j.jevs.2009.03.002.
Madureira RC, Rangel CP, Barbosa-Neto DJ, Fonseca AH. Sorologia para Borrelia Burgdorferi em equinos da ilha de Marajó e município de Castanhal, Pará. Brasil Rev Ciências da Vida. 2009;29:9–15.
Parra AC. Investigação diagnóstica de doença concomitante Babesiose e Anaplasmose em rebanho equino, por técnicas de Nested PCR, c – ELISA ou ELISA indireto. Universidade de São Paulo; 2009.
Salvagni CA, Dagnone AS, Gomes TS, Mota JS, Andrade GM, Baldani CD, et al. Serologic evidence of equine granulocytic anaplasmosis in horses from central West Brazil. Rev Bras Parasitol Vet. 2010;19:135–40.
Soares IR, Silva SO, Moreira FM, Prado LG, Fantini P, Maranhão R, et al. First evidence of autochthonous cases of Leishmania (Leishmania) infantum in horse (Equus caballus) in the Americas and mixed infection of Leishmania infantum and Leishmania (Viannia) braziliensis. Vet Parasitol. 2013;197:665–9. https://doi.org/10.1016/j.vetpar.2013.06.014.
Vedovello-Filho D, Jorge FA, Lonardoni MVC, Teodoro U, Silveira TGV. American cutaneous leishmaniasis in horses from endemic areas in the north-central Mesoregion of Paraná state, Brazil. Zoonoses Public Health. 2008;55:149–55.
Ricketts HT. Some aspects of Rocky Moutain spotted fever as shown by recent investigations. Med Rec. 1909;76:843–5.
Heinemann PG, Moore JJ. Experimental therapy of Rocky Moutain spotted fever: the preventive and curative action of a serum for spotted fever, and the inefficiency of sodium cacodylate as a curative agent for this disease in guinea pigs. J Infect Dis. 1912;10:294–304.
Brasil. Produção da Pecuária Municipal 2010. Inst Bras Geogr e Estatística -IBGE [Internet]. 2010;38:1–65. http://biblioteca.ibge.gov.br/visualizacao/livros/liv81830.pdf
Brasil. Normais Climatológicas (1961–1999), Versão Revisada e Ampliada [Internet]. Digital. INMET, editor. Ministério da Agricultura e Reforma Agrária; 2010. http://www.inmet.gov.br/portal/index.php?r=clima/normaisClimatologicas
Aragão HB, Fonseca F. Notas de Ixodologia. VIII. Lista e chave para os representantes da fauna ixodológica brasileira. Mem Inst Oswaldo Cruz. 1961;59:115–29.
Martins TF, Onofrio VC, Barros-Battesti DM, Labruna MB. Nymphs of the genus Amblyomma (Acari: Ixodidae) of Brazil: descriptions, redescriptions, and identification key. Ticks Tick Borne Dis. 2010;1:75–99. https://doi.org/10.1016/j.ttbdis.2010.03.002.
Barros-Battesti DM, Arzua M, Bechara GH. Carrapatos de importância médico-veterinária da região neotropical: um guia ilustrado para identificação de espécie. São Paulo: Vox/ICTTD-3/Butantan; 2006.
Cordeiro MD, De Meireles GS, Da Silva JB, De Souza MMS, Da Fonseca AH. Soroprevalência para Borrelia spp. em cães no município de Seropédica, estado do Rio de Janeiro*. Rev Bras Med Vet. 2012;34:251–6.
Machado RZ, Montassier HJ, Pinto AA, Lemos EG, Machado MRF, Valadão IFF, et al. An enzyme-linked immunosorbent assay (ELISA) for the detection of antibodies against Babesia bovis in cattle. Vet Parasitol. 1997;71:17–26.
Labruna MB, Horta MC, Aguiar DM, Cavalcante GT, Pinter A, Gennari SM, et al. Prevalence of Rickettsia infection in dogs from the urban and rural areas of Monte Negro Municipality, Western Amazon Brazil. Vector-Borne Zoonotic Dis. 2007;7:249–55.
Pinter A, Labruna MB. Isolation of Rickettsia rickettsii and Rickettsia bellii in cell culture from the tick Amblyomma aureolatum in Brazil. Ann N Y Acad Sci. 2006;1078:523–9.
Silveira I, Pacheco RC, Szabó MPJ, Ramos HGC, Labruna MB. Rickettsia parkeri in Brazil. Emerg Infect Dis. 2007;13:1111–3.
Lachaud L, Chabbert E, Dubessay P, Dereure J, Lamothe J, Dedet JP, et al. Value of two PCR methods for the diagnosis of canine visceral leishmaniasis and the detection of asymptomatic carriers. Parasitology. 2002;125:197–207.
de Bruijn MHL, Barker DC. Diagnosis of New World leishmaniasis: specific detection of species of the Leishmania braziliensis complex by amplification of kinetoplast DNA. Acta Trop. 1992;52:45–58.
Massung RF, Slater K, Owens JH, Nicholson WL, Mather TN, Solberg VB, et al. Nested PCR assay for detection of granulocytic Ehrlichiae. J Clin Microbiol. 1998;36:1090–5. http://jcm.asm.org/content/36/4/1090.short%5Cnpapers2://publication/uuid/30E23B41-1C9D-4E38-9078-ADB10C1BEAB4
Birkenheuer AJ, Levy MG, Breitschwerdt EB. Development and evaluation of a Seminested PCR for detection and differentiation of Babesia gibsoni (Asian Genotype) and B. canis DN. J Clin Microbiol. 2003;41:4172–7.
Lacerda L, Silva A, Freitas J, Cruz R, Said R, Munhoz A. Feline immunodeficiency virus and feline leukemia virus: frequency and associated factors in cats in northeastern Brazil. Genet Mol Res. 2017;16:8.
Santolin IDAC, Famadas KM, McIntosh D. Detection and identification of Rickettsia agents in ticks collected from wild birds in Brazil by polymerase chain reaction-restriction fragment length polymorphism (PCR-RFLP) analysis. Rev Bras Med Vet. 2013;35:68–73.
Mangold AJ, Bargues MD, Mas-Coma S. Mitochondrial 16S rDNA sequences and phylogenetic relationships of species of Rhipicephalus and other tick genera among Metastriata (Acari: Ixodidae). Parasitol Res. 1998;84:478–84.
McIntosh D, Bezerra RA, Luz HR, Faccini JLH, Gaiotto FA, Giné GAF, et al. Detection of Rickettsia bellii and Rickettsia amblyommii in Amblyomma longirostre (Acari: Ixodidae) from Bahia state, Northeast Brazil. Brazil J Microbiol. 2015;46:879–83.
Starkweather J. Linear Mixed Effects Modelling using R. [Internet]. 2010. http://researchsupport.unt.edu/class/Jon/Benchmarks/LinearMixedModels_JDS_Dec201
Dohoo IR, Martin S. Veterinary Epidemiologic Research. 1st ed. Charlottrtown: University Prince Edward Island; 2003.
R Development Core Team. R: a language and environment for statistical computing [Internet]. Vienna: R Foundation for Statistical Computing; 2016. http://www.r-project.org
Bates D, Mächler M, Bolker BM, Walker SC. Fitting linear mixed-effects models using lme4. J Stat Softw. 2015;67:1–48.
Krawczak FS, Nieri-Bastos FA, Nunes FP, Soares JF, Moraes-Filho J, Labruna MB. Rickettsial infection in Amblyomma cajennense ticks and capybaras (Hydrochoerus hydrochaeris) in a Brazilian spotted fever-endemic area. Parasit Vectors. 2014;7:7.
Basile RC, Yoshinari NH, Mantovani E, Bonoldi VN, da Macoris D, de Queiroz-Neto A. Brazilian borreliosis with special emphasis on humans and horses. Brazil J Microbiol. 2017;48:167–72.
de Oliveira PB, Harvey TV, Fehlberg HF, Rocha JM, Martins TF, da Acosta ICL, et al. Serologic and molecular survey of Rickettsia spp. in dogs, horses and ticks from the Atlantic rainforest of the state of Bahia Brazil. Exp Appl Acarol. 2019;78:431–42. https://doi.org/10.1007/s10493-019-00397-x.
Horta MC, Labruna MB, Sangione LA, Vianna MCB, Gennari SM, Galvão M a M, et al. Prevalence of antibodies to spotted fever group Rickettsiae in humans and domestic animals in a brazilian spotted fever-endemic area in the state of São Paulo, Brazil: serologic evidence for infection by Rickettsia Rickettsii and another spotted fever group Rickettsia. Am J Trop Med Hyg. 2004;71:93–7. http://www.ajtmh.org/docserver/fulltext/14761645/71/1/0700093.pdf?expires=1559505429&id=id&accname=guest&checksum=8D7A90D49597152263A71DA91CF742D2
Otomura FH, Sangioni LA, Pacheco RC, Labruna MB, Galhardo JA, Ribeiro MG, et al. Anticorpos anti-rickettsias do grupo da febre maculosa em equídeos e caninos no norte do Estado do Paraná, Brasil. Arq Bras Med Vet e Zootec. 2010;62:761–4.
Castagnolli KC, De Figueiredo LB, Santana DA, De Castro MB, Romano MA, Szabó MPJ. Acquired resistance of horses to Amblyomma cajennense (Fabricius, 1787) ticks. Vet Parasitol. 2003;117:271–83.
Lopes MG. Infecção por Rickettsia spp. em equídeos e carrapatos do Centro-Norte do Piauí. 2012;
Martins MEP, Brito WMED, Labruna MB, Moraes-Filho J. Identificação e pesquisa de Rickettsia spp. em carrapatos colhidos em cães e equinos de Quirinópolis, Goiás, Brasil. Multi-Science J. 2015;1:120–7.
BRASIL. DATASUS/FEBRE MACULOSA - Casos confirmados no Sistema de Informação de Agravos de Notificação [Internet]. SinanNet. 2018 [cited 2017 Dec 20]. http://tabnet.datasus.gov.br/cgi/deftohtm.exe?sinannet/cnv/febremaculosaba.def
Paddock CD, Brenner O, Vaid C, Boyd DB, Berg JM, Joseph RJ, et al. Short report: concurrent Rocky Mountain spotted fever in a dog and its owner. Am J Trop Med Hyg. 2002;66:197–9.
Sangioni LA, Horta MC, Vianna MCB, Gennari SM, Soares RM, Galvão MAM, et al. Rickettsial infection in animals and Brazilian spotted fever endemicity. Emerg Infect Dis. 2005;11:265–70.
Pinter A, Horta MC, Pacheco RC, Moraes-Filho J, Labruna MB. Serosurvey of Rickettsia spp. in dogs and humans from an endemic area for Brazilian spotted fever in State of São Paulo. Brazil Cad Saúde Pública. 2008;24:247–52.
Nasser JT, Lana RC, dos Silva CM, Lourenço RW, da Silva DC, Donalísio MR. Urbanization of Brazilian spotted fever in a municipality of the southeastern region: epidemiology and spatial distribution. Rev Bras Epidemiol. 2015;18:299–312.
Prado RFS, De CCHC, Cordeiro MD, Pires MS, Cepeda MB, Mafra C, et al. Seroprevalence of homologous anti-Borrelia burgdorferi antibodies in horses of military use in Brazil. Rev Bras Ciência Veterinária. 2017;24:72–6.
Salles RS, Fonseca AH, Madureira RC, Yoshinari NH. Sorologia para Borrelia burgdorferi latu sensu em equinos no estado do Rio de Janeiro. Hora Veterinária. 2002;127:46–9.
Madureira RC, Corrêa F do N, Cunha NC da, Guedes D da S, Fonseca AH da. Ocorrência de anticorpos homólogos anti-Borrelia burgdorferi em eqüinos de propriedades dos municípios de Três Rios e Vassouras, estado do Rio de Janeiro. Rev Bras Ciência Veterinária. 2007;14:43–6.
Carlos RSA, Muniz-Neta ES, Spagnol FH, Oliveira LLS, Brito RLLDE, Albuquerque GR. Frequência de anticorpos anti Erlichia canis, Borrelia burgdorferi e antígenos de Dirofilaria immitis em cães na microrregião Ilhéus-Itabuna, Bahia, Brasil. Rev Bras Parasitol Veterinária. 2007;16:117–20.
Yoshinari NH, Mantovani E, Bonoldi VLN, Marangoni RG, Gauditano G. Doença de Lyme-Símile Brasileira ou Síndrome Baggio-Yoshinari: Zoonose exótica e emergente transmitida por carrapatos. Rev Assoc Med Bras. 2010;56:363–9.
Cordeiro MD, Bahia M, Magalhães-Matos PC, Cepeda MB, Guterres A, Fonseca AH. Morphological, molecular and phylogenetic characterization of Borrelia theileri in Rhipicephalus microplus. Rev Bras Parasitol Vet. 2018;27:555–61.
Nascimento DA, da Vieira RF, Vieira TSWJ, dos Todelo R, Tamekuni K, dos Santos NJ, et al. Serosurvey of Borrelia in dogs, horses, and humans exposed to ticks in a rural settlement of southern Brazil. Brazil J Vet Parasitol. 2016;25:418–22.
Egenvall A, Franzén P, Gunnarsson A, Engvall EO, Vågsholm I, Wikström UB, et al. Cross-sectional study of the seroprevalence to Borrelia burgdorferi sensu lato and granulocytic Ehrlichia spp. and demographic, clinical and tick-exposure factors in Swedish horses. Prev Vet Med. 2001;49:191–208.
Funk RA, Pleasant RS, Witonsky SG, Reeder DS, Werre SR, Hodgson DR. Seroprevalence of Borrelia burgdorferi in horses presented for Coggins testing in Southwest Virginia and change in positive test results approximately 1 year later. J Vet Intern Med. 2016;30:1300–4.
Lee SH, Yun SH, Choi E, Park YS, Lee SE, Cho GJ, et al. Serological detection of Borrelia Burgdorferi among horses in Korea. Korean J Parasitol. 2016;54:97–101.
Passamonti F, Veronesi F, Cappelli K, Capomaccio S, Reginato A, Miglio A, et al. Polysynovitis in a horse due to Borrelia burgdorferi sensu lato infection – Case study. Ann Agric Environ Med. 2015;22:247–50.
Butler CM, Van O-O, Werners AH, Stout TAE. Borrelia burgdorferi and Anaplasma phagocytophilum in ticks and their equine hosts: a prospective clinical and diagnostic study of 47 horses following removal of a feeding tick. Pferdeheilkunde. 2016;32:335.
Brasil. Leishmaniose Tegumentar Americana - Casos confirmados notificados no sistema de informação de agravos de notificação. Ministério da Saúde. 2017;26:2–3.
Brasil. Produção da pecuária municipal 2015. Inst Bras Geogr e Estatística - IBGE. 2015;43:1–49.
Aguilar CM, Fernandez E, de Fernandez R, Deane LM. Study of an outbreak of cutaneous leishmaniasis in Venezuela. The role of domestic animals. Mem Inst Oswaldo Cruz. 1984;79:181–95.
Benassi JC, Benvenga GU, Ferreira HL, Soares RM, Silva DT, Pereira VF, et al. Molecular and serological detection of Leishmania spp. in horses from an endemic area for canine visceral leishmaniasis in southeastern Brazil. Pesqui Vet Bras. 2018;38:1058–63.
Follador I, Araujo C, Cardoso MA, Tavares-Neto J, Barral A, Miranda JC, et al. Surto de leishmaniose tegumentar americana em Canoa, Santo Amaro, Bahia. Brasil Rev Soc Bras Med Trop. 1999;32:497–503.
Slivinska K, Víchová B, Werszko J, Szewczyk T, Wróblewski Z, Peťko B, et al. Molecular surveillance of Theileria equi and Anaplasma phagocytophilum infections in horses from Ukraine, Poland and Slovakia. Vet Parasitol. 2016;215:35–7. https://doi.org/10.1016/j.vetpar.2015.10.025.
Dziegiel B, Adaszek Ł, Kalinowski M, Winiarczyk S. Equine granulocytic anaplasmosis. Res Vet Sci. 2013;95:316–20.
Franzén P, Aspan A, Egenvall A, Gunnarsson A, Åberg L, Pringle J. Acute clinical, hematologic, serologic, and polymerase chain reaction findings in Horses experimentally infected with a European strain of Anaplasma phagocytophilum. J Vet Intern Med. 2005;19:232–9. https://doi.org/10.1111/j.1939-1676.2005.tb02687.x.
Prado LG, Palhares MS, Bastos C V., da Silveira JAG, Ribeiro ÁAR, Miranda ALS, et al. Anaplasma phagocytophilum direct detection and exposure evidence in equines from two breeding farms from Minas Gerais, Brazil. Arq do Inst Biológico [Internet]. 2018;85:1–6. http://www.scielo.br/pdf/aib/v85/1808-1657-aib-85-e0902017.pdf
Lopes AP, Sousa S, Dubey JP, Ribeiro AJ, Silvestre R, Cotovio M, et al. Prevalence of antibodies to Leishmania infantum and Toxoplasma gondii in horses from the north of Portugal. Parasit Vectors. 2013;6:178.
Kouam MK, Diakou A, Kanzoura V, Papadopoulos E, Gajadhar AA, Theodoropoulos G. A seroepidemiological study of exposure to Toxoplasma, Leishmania, Echinococcus and Trichinella in equids in Greece and analysis of risk factors. Vet Parasitol. 2010;170:170–5. https://doi.org/10.1016/j.vetpar.2010.02.004.
We thank Fábio Rozemberg (Bio Rødberg) for making the Graphical abstract image.
To Fundação de Amparo à Pesquisa do Estado da Bahia–FAPESB and Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES) (no. 003/2014) for funding the airline tickets and accommodation and to Fundação de Amparo à Pesquisa do Estado da Bahia–FAPESB for the fellowship.
Ethics approval and consent to participate
The study was carried out according to the standards established by the Brazilian College of Ethics and Animal Welfare. The research proposal was approved by the Committee for Ethics in Research with Animals (protocol 002/2013) at the State University of Santa Cruz, Ilhéus, BA, Brazil.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Geographical indication of the municipalities participating in the study. A (Santa Cruz da Vitória); B (Floresta Azul); C (Ibicaraí); D (Itabuna); E (Itaju do Colônia), F (Itapé). Table S1. Distribution of the antibody titers observed against R. rickettsii and R. parkeri antigens in isolated responses in equids of the Ilhéus-Itabuna microregion, Bahia. Table S2. Distribution of infections by Rickettsia spp, Borrelia spp., and Leishmania spp. in naturally infected equids, according to municipality in the Ilhéus-Itabuna microregion, Bahia. Table S3. Generalized bivariate linear models for factors associated with Rickettsia spp. infection in naturally infected equids from the Ilhéus-Itabuna microregion, Bahia. Table S4. Generalized linear classic multivariate model for factors associated with Rickettsia spp. infection in naturally infected equids from the Ilhéus-Itabuna microregion, Bahia. Full model. Table S5. Generalized bivariate linear models for factors associated with Borrelia spp. infection in naturally infected equids from the Ilhéus-Itabuna microregion, in the state of Bahia. Table S6. Generalized linear mixed multivariate model for factors associated with Borrelia spp. infection in naturally infected equids from the Ilhéus-Itabuna microregion, Bahia. Full model. Table S7. Generalized bivariate linear models for factors associated with Leishmania spp. infection in naturally infected equids from the Ilhéus-Itabuna microregion, Bahia.
About this article
Cite this article
Costa, S.C.L., de Souza Freitas, J., Carvalho, F.S. et al. Frequency and factors associated of potential zoonotic pathogens (Borrelia spp., Rickettsia spp., Leishmania spp., and Anaplasma phagocytophilum) in equids in the state of Bahia, Brazil. Parasites Vectors 14, 275 (2021). https://doi.org/10.1186/s13071-021-04777-4