Molecular survey and genetic diversity of piroplasmids in equids from Midwestern Brazil

Levantamento molecular e diversidade genética de piroplasmídeos em equídeos do Centro-Oeste do Brasil

Abstract

We evaluated the distribution of piroplasmids in equids from the Mato Grosso state in Midwestern Brazil using molecular methods and the interspecific genetic diversity. For this, 1,624 blood samples of equids from 973 farms were examined by PCR, using primer pairs that amplify a fragment of the genes rap-1 and ema-1 of Babesia caballi and Theileria equi, respectively. For molecular characterization and phylogenetic studies, 13 and 60 sequences of the rap-1 and ema-1 genes, respectively, were used to build a dendogram using maximum parsimony. B. caballi and T. equi were detected in 4.11% and 28.16% of the farms, respectively, and molecular prevalence was 2.74% for B. caballi and 25.91% for T. equi. The location of the farms and animals raised in the Pantanal ecoregion influence the probability of equids testing positive for B. caballi and T. equi . Moreover, age and herd purpose were variables significantly associated with T . equi infection. The sequences of B. caballi presented 1.95% intraspecific variability, contrasting with 2.99% in T. equi. Dendrograms for both species demonstrated the presence of subgroups with high values of support of branches. However, it is not possible to associate these groups with geographic origin and/or ecoregion.

Keywords:
Babesia caballi, Theileria equi; equine babesiosis; equine piroplasmosis; phylogenetic analysis

Resumo

Foi avaliada a distribuição de piroplasmídeos em equídeos do Estado de Mato Grosso, no Centro-Oeste do Brasil, utilizando-se métodos moleculares e a diversidade genética interespecífica. Para isso, 1.624 amostras de sangue de equídeos de 973 fazendas foram examinadas pela PCR, usando pares de oligonucleotídeos que amplificam um fragmento dos genes rap-1and ema-1 de Babesia caballi e Theileria equi, respectivamente. Para caracterização molecular e estudos filogenéticos, foram utilizadas 13 e 60 sequências dos genes rap-1 e ema-1, respectivamente, para construção de um dendograma utilizando máxima parcimônia. B. caballi e T . equi foram detectados em 4,11% e 28,16% das fazendas, respectivamente, e a prevalência molecular foi de 2,74% para B. caballi e 25,91% para T. equi. A localização das fazendas e animais criados na ecorregião do Pantanal influenciam a probabilidade de equídeos serem positivos para B. caballi e T. equi. Além disso, idade e propósito do rebanho foram variáveis, significativamente, associadas à infecção por T. equi. As sequências de B . caballi apresentaram variabilidade intraespecífica de 1,95%, contrastando com 2,99% em T. equi. Dendrogramas para ambas as espécies demonstraram a presença de subgrupos com altos valores de sustentação dos ramos. No entanto, não é possível associar esses grupos com origem geográfica e/ou ecorregião.

Palavras-chave:
Babesia caballi; Theileria equi; babesiose equina; piroplasmose equina; análise filogenética

Introduction

Equine piroplasmosis is a tick-borne disease of horses and other Equidae (donkeys, mules, and zebras) caused by the hemoprotozoan parasites Theileria equi and Babesia caballi ( WISE et al., 2013 ; WISE et al., 2014 ). These two parasites are classified within the phylum Apicomplexa, and despite their biological differences, they cause a similar pathology and have similar life cycles and vector relationships. Acute disease is characterized by fever, malaise and reduced appetite, increased pulse rate and respiration, anorexia, constipation followed by diarrhea, tachycardia, petechiae, splenomegaly, thrombocytopenia, and hemolytic anemia leading to hemoglobinuria and icterus ( MEHLHORN & SCHEIN, 1998 ; UILENBERG, 2006 ; ROTHSCHILD, 2013 ; WISE et al., 2013 ; HABIBI et al., 2016 ).

Horses that recover from acute disease remain chronically-infected carriers without any overt signs of disease and can be reservoirs for the transmission of these protozoan pathogens with ticks as a major vector. Parasitemia is often too low to detect on blood smears, but infected animals can be identified by serology or polymerase chain reaction (PCR) ( WISE et al., 2013 ; WISE et al., 2014 ; SCOLES & UETI, 2015 ). Ticks are the definitive hosts for both T. equi and B. caballi, and in Brazil, these two parasites are mainly transmitted by Rhipicephalus (Boophilus) microplus and Dermacentor nitens, respectively ( SCOLES & UETI 2015 ). Furthermore, equine piroplasmosis has an economic importance because of the potential effect of the international movement of these horses for trade and for participation in international sporting events ( WISE et al., 2013 ; WISE et al., 2014 ; SCOLES & UETI, 2015 ).

Piroplasmosis occurs in most countries around the world and infection is maintained within equine populations as long as competent vectors are present. According to the World Organization for Animal Health (OIE), Central and South America, Cuba, Africa, Asia, the Middle East, and Southern Europe are considered to be endemic regions. Within South America, the disease is readily identified in all regions, except for the southernmost areas of Chile and Argentina ( WISE et al., 2013 ). In Brazil, previous studies using serological methods have shown that equine piroplasmosis is endemic, with prevalence rates varying from 19.4% to 100% ( KERBER et al., 1999 ; HEIM et al., 2007 ; KERBER et al., 2009 ; BALDANI et al., 2010 ; SANTOS et al., 2011 ; MACHADO et al., 2012 ; TORRES et al., 2012 ; VIEIRA et al., 2013 ; PROCHNO et al., 2014 ; FERREIRA et al., 2016 ; GUIMARÃES et al., 2016 ; BRAGA et al., 2017 ).

Despite the endemicity of piroplasmosis among equids in Brazil, studies involving the molecular detection of the parasite in carrier animals ( HEIM et al., 2007 ; BALDANI et al., 2010 ; PECKLE et al., 2013 ; MACHADO et al., 2012 ; BARROS et al., 2015 ; FERREIRA et al., 2016 ; BRAGA et al., 2017 ) that have also examined genetic diversity ( BRAGA et al., 2017 ; PECKLE et al., 2018 ) are scarce. The aim of this study is to survey the distribution of piroplasmid (T . equi and B. caballi) parasites using molecular methods in equids from the Mato Grosso state and evaluate the interspecific genetic diversity.

Materials and Methods

Study area

With nearly one thousand square kilometers, the state of Mato Grosso is located in Midwestern Brazil, encompassing biotic elements from Amazonia, Cerrado (Brazilian savannah), and Pantanal ecoregions. Southern Amazon forests correspond to nearly 53% of the state and are located in its northern portion. The Cerrado – a huge tropical savannah that extends across the Brazilian Central Plateau – occupies nearly 40% of the state, in its central portion; and the Pantanal – the largest tropical wetland in the world – occupies another 7% of the state area, in its southwestern portion ( BRASIL, 2014 ).

Blood samples were collected during the course of a previous study on the epidemiology of equine infectious anemia (EIA) within the three ecoregions in the state of Mato Grosso. This study was conducted from September to December 2014 ( BARROS, 2017 ). The stored samples were made available for this study after the conclusion of the previous study.

Concomitant with sampling, each farmer was asked to fill out a questionnaire for subsequent epidemiological data analysis.

Sampling procedures

According to the Mato Grosso Institute for Agricultural Defense ( INDEA/MT, 2014 ), there are currently 306,931 equids on 53,710 farms distributed through three ecoregions (Cerrado, Pantanal and Amazonia) in the state, as observed in table 1 .

Table 1
Columns left to right: (1) Ecoregion in Mato Grosso state, Brazil. (2) Total number of farms in each ecoregion. (3) Number of farms sampled in the current survey. (4) Total number of equids in each ecoregion. (5) Number of equids sampled in the current survey. All the data were collected between September and December 2014.

Ecoregion1	No. of farms2	No. of sampled farms3	No. of equids4	No. of sampled equids5
Amazonia	32,069	394	165,413	661
Cerrado	18,114	312	110,656	498
Pantanal	3,527	267	30,862	465
Total	53,710	973	306,931	1,624

The sample size was determined by the statistical formula $n_0=[Np(1-p)/(N-1){\left(d/Z_{\alpha /2}\right)}^2+p(1-p)].deff$ ( SCHAEFFER et al., 2011 ). This formula considered population size (306,931 animals); 50% estimated prevalence (due the absence of prior studies and the maximum number of occurrences stay around 50% on the normal distribution); 3% maximum error (d); 95% confidence interval (CI); and a design effect of 1.45 (deff). The sample size was adjusted from 1,542 to 1,624 animals (1,389 equines, 227 donkeys and eight asinines) for the possibility of loss, and then distributed among the three ecoregions ( Table 1 ). Considering that 77% of the farms had no more than five equids, the following criteria were used: only one animal was sampled on farms with five or fewer equids, while five animals were randomly sampled on farms that had more than five equids.

Sampling was performed in two stages: selection of farms and selection of Equidae. Selection procedures were as follows: a pre-established number of random farms (primary sampling units) was selected, which resulted in the sampling of 973 farms ( Table 1 ). Next, a random selection of equids on each farm (secondary sampling units) was performed. Figure 1 shows the distribution of the farms sampled for the screening of B. caballi and T. equi infection in the three ecoregions (Amazonia, Pantanal and Cerrado) of the state of Mato Grosso.

Figure 1
Distribution of farms sampled for the detection of Babesia caballi and Theileria equi in equids from the three ecoregions (Pantanal, Cerrado and Amazonia) of the Mato Grosso state, Brazil, between September and December 2014.

DNA extraction and PCR

DNA was extracted from 60 μg of each blood clot sample by the phenol-chloroform- isoamyl method as described by Sambrook & Russel (2001) . In order to verify the presence of amplifiable DNA, 163 randomly selected samples were submitted to internal control PCR assays targeting fragments of mammalian glyceraldehyde-3-phosphate dehydrogenase (GAPDH), as previously described ( BIRKENHEUER et al., 2003 ). Each sample of extracted DNA was submitted for PCR analysis, using species-specific primer pairs based on the genes rap-1 for B. caballi and ema-1 for T. equi, which amplified DNA fragments of 570 and 750 base pairs (bp), respectively. This PCR analysis procedure was adapted from Alhassan et al. (2005) . The reaction mixture consisted of 1.5 µl template DNA (B. caballi ) or 2 µl (T. equi), 2.5 µl 10× PCR buffer, 4.0 µl MdNTP mix, 0.15 µl Taq DNA polymerase (Sigma-Aldrich), 0.5 µl each of 10 pmol B. caballi and T. equi specific primers, and 15.35 µl of distilled water (B. caballi) or 15.85 µl ( T. equi), for a total volume of 25 µl. The amplification conditions for B. caballi and T. equi were as follows: 40 cycles with an initial enzyme activation at 96 °C for 10 min, denaturation at 96°C for 1 min, primer annealing at 60.5 °C for 1 min, extension at 72 °C for 1 min and final extension at 72 °C for 10 min. For all PCR assay, B. caballi and T . equi DNA samples obtained from naturally infected horse ( BARROS et al., 2015 ) and ultra-pure sterile water were used as positive and negative controls, respectively. PCR product sizes were analyzed on 1.5% agarose gel stained with GelRed™ Nucleic Acid Gel Stain, and visualized in a ChemiDoc XRS system.

Sequencing and phylogenetic analysis

For molecular characterization and phylogenetic studies of B. caballi and T. equi isolates, fragments of 569 and 675 bp from rap-1 and ema-1 genes, respectively, were used. DNA from positive samples was purified and sequenced in an automated sequencer (ABI-PRISM 3500 Genetic Analyzer, Foster City, CA).

All the sequences obtained were aligned with sequences previously determined for B. caballi and T. equi, available in GenBank, using ClustalX ( THOMPSON et al., 1997 ). The sequences were adjusted manually using GeneDoc ( NICHOLAS et al., 1997 ).

The rap-1 and ema-1 sequences were used to build a dendogram using maximum parsimony, as implemented in PAUP version 4.0b10 ( SWOFFORD, 2002 ), with 500 bootstrap replicates, random stepwise addition starting trees (with random addition sequences), and TBR branch swapping.

Prevalence and statistical analysis

Piroplasmid prevalence on each farm was analyzed based on the logic that a farm would be considered positive if at least one animal tested positive. A multiple logistic regression model was created in several steps as follows: 1. For continuous variables, a linearity test was performed using graphs for log odds. When the visual evaluation indicated a nonlinear trend, the continuous variable was categorized. The optimization methodology of the software package CatPredi R was used ( R DEVELOPMENT CORE TEAM, 2015 ). 2. For the categorical variables, a chi-square or exact Fisher test was performed on the results and various predictor variables. All variables with a value of p <0.20 were incorporated into the model. 3. The selected variables were tested for collinearity including a multicollinearity analysis according to Vatcheva et al. (2016) , to ensure a mean variance inflation factor (VIF) <10 prior to applying them in the models. 4. The model was constructed by inserting all the variables selected in the previous steps. The less-significant variables (according to the Wald test) were then removed and the logistic regression analysis was then repeated. This process was repeated, and the adjusted model was compared to the previous one by a likelihood ratio test to check for confounding factors. When there was a change of more than 30% in the parameter estimates, the removed variable was considered to be a confounding factor and again included in the model.

The variables used to construct the herd model were ecoregion (biome), type of holding (farm, breeding facility or horseback riding facility), herd purpose (service, sport/recreation or reproduction), type of management (extensive, intensive or semi-extensive), needle/syringe sharing among animals, flooded area at the farm, veterinary assistance, and animal density (greater than 5 hectares per animal or less than 5 hectares per animal).

The animals were analyzed following the same steps used to analyze the farms. However, a complex analysis, according to Korn & Graubard (1999) , was applied to the tests to correct the effect of the weighted values used for the selection of animals. The following formula was used to determine the weighted value of each animal: (# total animals/# total animals sampled) x (# total animals on the farm/# total animals sampled on the farm).

The model was built using the following variables: ecoregion (biome), animal species, sex, age (classified as greater than six months and less than 3 years, greater than 3 years and less than 10 years, or greater ten than years), and animal purpose (work or sports/recreation). Statistical analysis was performed using the R statistical software package ( R DEVELOPMENT CORE TEAM, 2015 ).

This study was approved by the Animal Research Ethics Committee of the Federal University of Mato Grosso, under reference no. 23108.016503/14-9.

Results

B. caballi and T. equi were detected in at least one animal on 40 and on 274 of 973 analyzed farms ( Figure 1 ), resulting in a herd prevalence of 4.11% (CI 95%, 2.86%; 5.36%) and 28.16% (CI 95%; 25.33%; 30.99%), respectively. We observed distribution of herd prevalence for B. caballi and T. equi, respectively, through the three ecoregions (biomes) of the Mato Grosso as follows: 5.33% (CI 95%; 3.11%; 7.55%) and 35.03% (CI 95%; 30.31%; 39.74%) in the Amazonia, 4.81% (CI 95%; 2.43%; 7.18%) and 27.88% (CI 95%; 22.91%; 32.86%) in the Cerrado, and 1.50% (CI 95%; 0.04%; 2.96%) and 18.35% (CI 95%; 13.71%; 23.00%) in the Pantanal ecoregion. All 163 selected DNA samples amplified the predicted product for GAPDH gene.

Moreover, molecular prevalence of the total population was 2.74% (CI 95%; 1.76%; 3.71%) and 25.91% (CI 95%; 21.67%; 30.15%) for B. caballi and T. equi, respectively, and mixed infection was detected in 1.97% (32/1624). Distribution of prevalence among the three ecoregions (biomes) are as follows for B. caballi and T. equi , respectively: 3.63% (CI 95%; 2.20%; 5.06%) and 30.26% (CI 95%; 26.76%; 33.76%) from the Amazonia, 3.41% (CI 95%; 1.82%; 5.01%) and 28.31% (CI 95%; 24.36%; 32.27%) from the Cerrado, and 0.86% (CI 95%; 0.02%; 1.70%) and 12.69% (CI 95%; 9.66%; 15.71%) from the Pantanal ecoregion.

In the multivariate analysis, after the model was applied, location (ecoregion) was the only variable found to be significantly associated with the prevalence of babesiosis on farms. Farms in the Pantanal ecoregion were 73% less likely to have equine infected with B. caballi . Table 2 describes the relationship between the positive farms and the variable of the logistic model. Furthermore, animals raised in the Pantanal ecoregion were 3% less likely to be infected with B. caballi ( Table 3 ).

Table 2
Results of the logistic regression model for risk factors associated with Babesia caballi herd prevalence on the 973 farms sampled in the state of Mato Grosso, Brazil, between September and December 2014.

Associated variable	Estimate	p-value	Odds Ratio	CI 95%
Location in the Pantanal	-1.30	0.01	0.27	0.07	0.71

Table 3
Results of the logistical regression model for risk factors associated with the prevalence of infection by Babesia caballi in individual equids from the state of Mato Grosso, Brazil, between September and December 2014.

Associated variable	Estimate	p-value	Odds Ratio	CI 95%
Location in the Pantanal	-0.02	<0.01	0.97	0.95	0.98

Two variables were found to be significantly associated with the prevalence of theileriosis on farms: location (ecoregion) in the state, and herd purpose. Table 4 describes the relationship between the positive farms and the variables of the logistic model. Farms in the Pantanal ecoregion were 57% less likely to have equine infected with T. equi, and farms that used equids to perform service were 84% less likely to have T. equi infected animals.

Table 4
Results of the logistic regression model for risk factors associated with Theileria equi herd prevalence on the 973 farms sampled in the state of Mato Grosso, Brazil, between September and December 2014.

Associated variable	Estimate	p-value	Odds Ratio	CI 95%
Location in the Pantanal	-0.29	< 0.01	0.43	0.29	0.63
Herd purpose (service, sport/recreation or reproduction)	-17.73	< 0.01	0.16	0.06	0.42

Age and location of equids were variables significantly associated with T . equi infection. Equids between six months and three years old were 1.2 times (20%), more likely to be infected with T. equi, and equines older than 10 years were 10% less likely to have T. equi infection. Furthermore, animals raised in the Pantanal ecoregion were 20% less likely to be infected with T. equi . Table 5 describes the relationship between the prevalence of infection by T. equi in the Equidae population and the variables of the logistic model of the 1,624 animals sampled.

Table 5
Results of the logistical regression model for risk factors associated with the prevalence of infection by Theileria equi in individual equids from the state of Mato Grosso, Brazil, between September and December 2014.

Associated variable	Estimate	p-value	Odds Ratio	CI 95%
Age* (greater than six months and less than three years)	0.18	<0.01	1.2	1.08	1.33
Age (greater than ten years)	-0.09	<0.01	0.9	0.85	0.95
Location in the Pantanal	-0.21	<0.01	0.8	0.73	0.87

• ^*
  Animal age, classified as greater than six months and less than 3 years, greater than 3 years and less than 10 years and greater than ten years.

Thirteen sequences of B. caballi were obtained from 13 farms and 60 sequences of T. equi were obtained from 52 properties. The sequences of B. caballi presented 1.95% intraspecific variability, contrasting with 2.99% in T. equi. Dendrograms based on maximum parsimony for B. caballi ( Figure 2 ) and T. equi ( Figure 3 ) demonstrated the presence of subgroups with high support values of branches ( Figure 3 ). The GenBank nucleotide sequence accession numbers generated in the present study are: MG906574- MG906586 for partial sequences of the rap-1 gene of B . caballi and MG906587– MG906646 for partial sequences of the ema-1 gene of T. equi.

Figure 2
Maximum parsimony tree inferred from rap-1 gene sequences of Babesia caballi. Numbers at nodes are the support values for the major branches (bootstrap over 500 replicates).

Figure 3
Maximum parsimony tree inferred from ema-1 gene sequences of Theileria equi. Numbers at nodes are the support values for the major branches (bootstrap over 500 replicates).

Discussion

In this study, we investigated the occurrence of piroplasmids (T. equi and B. caballi) using molecular methods in 1,624 equids (1,389 equine, 227 donkeys, eight asinines) from 973 farms from an extensive region of Midwestern Brazil, which covers three different ecoregions (Amazonia, Cerrado and Pantanal) of the Mato Grosso state. Furthermore, we evaluated variables (risk factor) thought to be associated with the occurrence of infection in animals and Equidae herds, as well as the interspecific genetic diversity of B. caballi and T. equi.

Despite the fact that B. caballi and T . equi have previously been detected in the Pantanal ecoregion ( BARROS et al., 2015 ), this is the first study, to our knowledge, to show the prevalence of piroplasmids in Equidae in the entire state. Herein, we confirm that infection is widespread in the studied region, with higher prevalence of T. equi (25.91%) than B. caballi (2.74%) among equids. These findings are similar to those found in other studies conducted in Brazil ( HEIM et al., 2007 ; VIEIRA et al., 2018 ) and other countries ( POSADA-GUZMÁN et al., 2015 ; GUVEN et al., 2017 ). This may reflect differences in the biology of the vectors, as well as differences in persistence of the parasite in the equine host ( SCOLES & UETI, 2015 ). In most cases, the horses develop a chronic infection and become asymptomatic carriers. While horses infected with B. caballi can undergo self-clearance of the parasite without treatment, this is a lifelong chronic infection for equids infected with T. equi ( WISE et al., 2013 ). Furthermore, the tick vector R. (B.) microplus is the only known biological vector responsible for the intrastadial and transstadial transmission of T. equi in South America, while D. nitens is the only transovarial vector for B . caballi in the Americas ( SCOLES & UETI, 2015 ).

Unexpectedly, when the purpose of the herd was for service, the likelihood of a farm to be positive for T. equi decreased by 83.02%, once horse's activities related to work have already been described as factors that affect the likelihood of infection by T. equi ( PECKLE et al., 2013 ; ABUTARBUSH et al., 2012 ). Furthermore, the activities that the horse exerts on the farm, mainly, cattle handling, allows move into large pasture perimeters and may increase their chances of infestation by potential tick vector ( PECKLE et al., 2013 ).

In the present study, farms located in the Pantanal are less likely to have animals infected with B. caballi and T. equi . Moreover, animals raised in the Pantanal ecoregion were less likely to be infected with both protozoans (B. caballi and T. equi). In fact, the lowest prevalence among farms and the Equidae population, both B. caballi and T. equi , were observed in the Pantanal ecoregion.

The Pantanal ecosystem is characterized by flood and drought cycles, where over 80% of the plains are submerged during the rainy period ( ALHO, 2008 ). Therefore, the total or partial flooding of pastures may eliminate ticks from the vegetation, and could also affect tick survival. Ramos et al. (2016) have studied ecological relationships between Nellore cattle (Bos indicus ) and ticks within the Brazilian Pantanal, and although R. (B. ) microplus (the known vector of T. equi ) was the most abundant tick species on cattle (primary tick-host), the overall highest tick prevalence on bovines in the dry season was that of Amblyomma sculptum nymphs. Additionally, during the entire study, R. (B.) microplus was not found in the environment.

According to Campbell & Glines (1979) , conditions of greater or lesser flooding of land have a great influence on the biology of Ixodidae, in addition to other environmental factors such as ambient temperature and humidity. Louzada & Daemon (2003) and Santos et al. (2012) have verified the effects of the immersion of engorged females of R. ( B.) microplus, and observed a negative effect on reproduction parameters, through a negative effect over the weight of the eggs, alterations in the egg-laying and in the hatching. Moreover, Paula et al. (2004) have observed that immersion in water is deleterious to the free-living phase of D . nitens larvae (the known vector of B. caballi). Thus, the differences between ecoregions and their edaphoclimatic conditions should influence the epidemiology of equine piroplasmosis.

T. equi seroprevalence rates increase with age, with the older group related to a chronic infectious status ( DE WAAL, 1992 ). Santos et al. (2011) have observed in 714 equids (702 horses and 12 mules) from the Rio de Janeiro state, that animals aged six to 12 months showed 87.5% of T. equi seropositivity. Furthermore, considering the absence of maternal antibodies in this age group, they have suggested that the first T. equi infection occurred between six and 12 months of age or possibly even younger since there is no passive immunity in animals over six months old. We also observed by molecular methods that equines older than six months and less than three years have a 1.2-fold increased risk of T. equi infection, however animals older than 10 years have a 10.08% decrease in the chance of being infected with T. equi . Our results are in agreement with Bartolomé Del Pino et al. (2016) , who also showed prevalence for T. equi decreasing with age as detected by PCR. However, they are contrary to some reports in literature that have highlighted a cumulative age-dependent increase in the PCR-prevalence of T. equi ( RÜEGG et al., 2008 ) and to others that have shown no association with age and infection rates ( PECKLE et al., 2013 ). Hence, the circumstances that influence positivity rates of T. equi during the lifetime of a horse, require further study due to the sensitivity of the various methods used in studies involving host–parasite interactions ( BARTOLOMÉ DEL PINO et al., 2016 ). Indeed, the associations of age and infection that are reported in the literature may be due to the small sample size of animals in some categories, which could possibly interfere with the analysis of the data ( SANTOS et al., 2011 ).

Characterizations of piroplasm sequences from Brazilian horses are scarce. Intraspecific analysis of the sequences obtained from equids of different ecoregions in the state of Mato Grosso showed 1.95% sequence divergence for B. caballi and 2.99% for T. equi. Studies based on ribosomal (18S) genes demonstrated variability between 0.5% and 2.1% for T. equi, allowing the segregation into different genotypes ( BRAGA et al., 2017 ; PECKLE et al., 2018 ).

The intraspecific divergence of T. equi found within the ema-1 gene was greater than that presented by the 18S gene. However, it was not possible to distinguish between subgroups/genotypes related to the different ecogeographic regions studied. The absence of clear patterns related to the geographic region and/or ecoregion may be due to the internal migrations of the animals, making it impossible to form the barriers or geographic isolation needed to identify clusters.

Further studies with the generation of new sequences of B. caballi and T. equi from different regions and Brazilian biomes could minimize the effects of animal traffic in the analysis of intraspecific variability.

Conclusions

Equine piroplasmosis occurs throughout the three ecoregions of the Mato Grosso state in Midwestern Brazil, with the Amazonia ecoregion containing the largest number of animals and farms infected by both agents (B. caballi and T. equi), while farms and animals raised in the Pantanal ecoregion were less likely to be infected with both protozoans. Animals between six months and two years and eleven months old are at an increased risk of T. equi infection. Interestingly, animals over the age of ten years were shown to be at reduced risk of T . equi infection. Intraspecific analysis of the sequences obtained for B. caballi and T. equi did not allow us to distinguish between subgroups/genotypes related to the different ecogeographic regions under area studied.

Acknowledgements

The authors gratefully acknowledge INDEA/MT for its technical support during the field work and for providing samples, as well as to FAPEMAT, Foundation for Research Support of the State of Mato Grosso (grant #223183/2015) for its financial support of this work, and CNPq (National Council for Scientific and Technological Development) for the awarding a research productivity grant to L. Nakazato, R.C. Pacheco and V. Dutra.

References

• Abutarbush SM, Alqawasmeh DM, Mukbel MR, Al-Majali AM. Equine Babesiosis: seroprevalence, risk factors and comparison of different diagnostic methods in Jordan. Transbound Emerg Dis 2012; 59(1): 72-78. http://dx.doi.org/10.1111/j.1865-1682.2011.01244.x. PMid:21771287.
  » http://dx.doi.org/10.1111/j.1865-1682.2011.01244.x
• Alhassan A, Pumidonming W, Okamura M, Hirata H, Battsetseg B, Fujisaki K, et al. Development of a single-round and multiplex PCR method for the simultaneous detection of Babesia caballi and Babesia equi in horse blood. Vet Parasitol 2005; 129(1-2): 43-49. http://dx.doi.org/10.1016/j.vetpar.2004.12.018. PMid:15817201.
  » http://dx.doi.org/10.1016/j.vetpar.2004.12.018
• Alho CJR. Biodiversity of the Pantanal: response to seasonal flooding regime and to environmental degradation. Braz J Biol 2008;68(4 Suppl): 957-966. http://dx.doi.org/10.1590/S1519-69842008000500005. PMid:19197468.
  » http://dx.doi.org/10.1590/S1519-69842008000500005
• Baldani CD, Nakaghi ACH, Machado RZ. Occurrence of Theileria equi in horses raised in the Jaboticabal microregion, São Paulo State, Brazil. Rev Bras Parasitol Vet 2010; 19(4): 228-232. http://dx.doi.org/10.1590/S1984-29612010000400007. PMid:21184699.
  » http://dx.doi.org/10.1590/S1984-29612010000400007
• Barros EM, Braga IA, Santos LGF, Ziliani TF, Melo ALT, Borges AMCM, et al. Detecção de Theileria equi e Babesia caballi e anticorpos anti- Ehrlichia spp. em equídeos do Pantanal Mato-Grossense, Brasil. Arq Bras Med Vet Zootec 2015; 67(3): 716-722. http://dx.doi.org/10.1590/1678-4162-7930.
  » http://dx.doi.org/10.1590/1678-4162-7930
• Barros ML. Prevalência, distribuição espacial de focos e possíveis fatores de risco associados à anemia infecciosa equina (aie) no estado de Mato Grosso [dissertation]. Cuiabá: Universidade Federal de Mato Grosso; 2017.
• Bartolomé Del Pino LE, Nardini R, Veneziano V, Iacoponi F, Cersini A, Autorino GL, et al. Babesia caballi and Theileria equi infections in horses in Central-Southern Italy: Sero-molecular survey and associated risk factors. Ticks Tick Borne Dis 2016; 7(3): 462-469. http://dx.doi.org/10.1016/j.ttbdis.2016.01.011. PMid:26847198.
  » http://dx.doi.org/10.1016/j.ttbdis.2016.01.011
• 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 DNA in canine blood samples. J Clin Microbiol 2003; 41(9): 4172-4177. http://dx.doi.org/10.1128/JCM.41.9.4172-4177.2003. PMid:12958243.
  » http://dx.doi.org/10.1128/JCM.41.9.4172-4177.2003
• Braga MDSCO, Costa FN, Gomes DRM, Xavier DR, André MR, Gonçalves LR, et al. Genetic diversity of piroplasmids species in equids from island of São Luís, northeastern Brazil. Rev Bras Parasitol Vet 2017; 26(3): 331-339. http://dx.doi.org/10.1590/S1984-29612017046. PMid:28977247.
  » http://dx.doi.org/10.1590/S1984-29612017046
• Brasil. Governo do Estado de Mato Grosso. Secretaria de Estado de Meio Ambiente – SEMA. Mapa dos Biomas Mato-Grossenses Cuiabá; 2014 [cited 2018 Apr 30]. Available from: http://www.sema.mt.gov.br/index.php?option=com_content&view=article&id=170&Itemid=107
  » http://www.sema.mt.gov.br/index.php?option=com_content&view=article&id=170&Itemid=107
• Campbell A, Glines MV. Development, survival, and oviposition of the rabbit tick, Haemaphysalis leporispalustris (Packard) (Acari: Ixodidae), at constant temperatures. J Parasitol 1979; 65(5): 777-781. http://dx.doi.org/10.2307/3280361. PMid:512768.
  » http://dx.doi.org/10.2307/3280361
• de Waal DT. Equine piroplasmosis: a review. Br Vet J 1992; 148(1): 6-14. http://dx.doi.org/10.1016/0007-1935(92)90061-5. PMid:1551016.
  » http://dx.doi.org/10.1016/0007-1935(92)90061-5
• Ferreira EP, Vidotto O, Almeida JC, Ribeiro LPS, Borges MV, Pequeno WHC, et al. Serological and molecular detection of Theileria equi in sport horses of northeastern Brazil. Comp Immunol Microbiol Infect Dis 2016; 47: 72-76. http://dx.doi.org/10.1016/j.cimid.2016.06.004. PMid:27477509.
  » http://dx.doi.org/10.1016/j.cimid.2016.06.004
• Guimarães AM, Bruhn FRP, Ribeiro MJM, Rosa MHF, Garcia AM, Rocha CMBM, et al. Prevalence and risk factors of Theileria equi infection in horses in Minas Gerais, Brazil. Vet Parasitol 2016; 3-4: 18-22. http://dx.doi.org/10.1016/j.vprsr.2016.05.005.
  » http://dx.doi.org/10.1016/j.vprsr.2016.05.005
• Guven E, Avcioglu H, Deniz A, Balkaya İ, Abay U, Yavuz Ş, et al. Prevalence and molecular characterization of Theileria equi and Babesia caballi in jereed horses in Erzurum, Turkey. Acta Parasitol 2017; 62(1): 207-213. http://dx.doi.org/10.1515/ap-2017-0025. PMid:28030350.
  » http://dx.doi.org/10.1515/ap-2017-0025
• Habibi G, Esmaeilnia K, Hablolvarid MH, Afshari A, Zamen M, Bozorgi S. Microscopic and molecular detection of Theileria (Babesia) equi infection in equids of Kurdistan province, Iran. Iran J Parasitol 2016; 11(1): 86-90. PMid:27095973.
• Heim A, Passos LF, Ribeiro MFB, Costa-Júnior LM, Bastos CV, Cabral DD, et al. Detection and molecular characterization of Babesia caballi and Theileria equi isolates from endemic areas of Brazil. Parasitol Res 2007; 102(1): 63-68. http://dx.doi.org/10.1007/s00436-007-0726-1. PMid:17828553.
  » http://dx.doi.org/10.1007/s00436-007-0726-1
• Instituto de Defesa Agropecuária do Estado de Mato Grosso – INDEA/MT. Relatório de equinos existentes [online]. Cuiabá: INDEA; 2014 [cited 2014 Aug 14]. Available from: http://www3.indea.mt.gov.br/defesa-sanitaria-animal/sanidade-equidea/
  » http://www3.indea.mt.gov.br/defesa-sanitaria-animal/sanidade-equidea/
• Kerber CE, Ferreira F, Pereira MC. Control of equine piroplasmosis in Brazil. Onderstepoort J Vet Res 1999; 66(2): 123-127. PMid:10486829.
• Kerber CE, Labruna MB, Ferreira F, de Waal DT, Knowles DP, Gennari SM. Prevalence of equine piroplasmosis and its association with tick infestation in the State of São Paulo, Brazil. Rev Bras Parasitol Vet 2009; 18(4): 1-8. http://dx.doi.org/10.4322/rbpv.01804001. PMid:20040201.
  » http://dx.doi.org/10.4322/rbpv.01804001
• Korn EL, Graubard BI. Analysis of health surveys New York: John Wiley and Sons; 1999. (Wiley Series in Probability and Statistics Survey Methodology Section).
• Louzada GL, Daemon E. Efeito da imersão de fêmeas ingurgitadas de Boophilus microplus (Canestrini, 1887) (Acari: Ixodidae) em água destilada sobre os parâmetros biológicos ligados à oviposição. Rev Bras Parasitol Vet 2003; 12(3): 115-120.
• Machado RZ, Toledo CZP, Teixeira MCA, André MR, Freschi CR, Sampaio P. Molecular and serological detection of Theileria equi and Babesia caballi in donkeys (Equus asinus) in Brazil. Vet Parasitol 2012; 186(3-4): 461-465. http://dx.doi.org/10.1016/j.vetpar.2011.11.069. PMid:22186194.
  » http://dx.doi.org/10.1016/j.vetpar.2011.11.069
• Mehlhorn H, Schein E. Redescription of Babesia equi Laveran, 1901 as Theileria equi Mehlhorn, Schein, 1998. Parasitol Res 1998; 84(6): 467-475. http://dx.doi.org/10.1007/s004360050431. PMid:9660136.
  » http://dx.doi.org/10.1007/s004360050431
• Nicholas KB, Nicholas HB Jr, Deerfield DW. GeneDoc: analysis and visualization of genetic variation. Embnew News 1997; 4: 14.
• Paula AR, Piranda EM, Faccini JLH, Daemon E. Efeito da imersão de larvas não alimentadas de Anocentor nitens (Neumann, 1897) e de Amblyomma cajennense (Fabricius, 1787) (Acari: Ixodidae) em água destilada. Rev Bras Parasitol Vet 2004; 13(1): 13-17.
• Peckle M, Pires MS, Santos TM, Roier ECR, Silva CB, Vilela JAR, et al. Molecular epidemiology of Theileria equi in horses and their association with possible tick vectors in the state of Rio de Janeiro, Brazil. Parasitol Res 2013; 112(5): 2017-2025. http://dx.doi.org/10.1007/s00436-013-3360-0. PMid:23474658.
  » http://dx.doi.org/10.1007/s00436-013-3360-0
• Peckle M, Pires MS, Silva CB, Costa RL, Vitari GLV, Senra MVX, et al. Molecular characterization of Theileria equi in horses from the state of Rio de Janeiro, Brazil. Ticks Tick Borne Dis 2018; 9(2): 349-353. http://dx.doi.org/10.1016/j.ttbdis.2017.11.011. PMid:29223587.
  » http://dx.doi.org/10.1016/j.ttbdis.2017.11.011
• Posada-Guzmán MF, Dolz G, Romero-Zúñiga JJ, Jiménez-Rocha AE. Detection of Babesia caballi and Theileria equi in blood from equines from four indigenous communities in Costa Rica. Vet Med Int 2015; 2015: 1-6. http://dx.doi.org/10.1155/2015/236278. PMid:26649225.
  » http://dx.doi.org/10.1155/2015/236278
• Prochno HC, Scorsin LM, Melo FR, Baldani CD, Falbo MK, Aquino LCT, et al. Seroprevalence rates of antibodies against Theileria equi in team roping horses from central-western region of Paraná. Rev Bras Parasitol Vet 2014; 23(1): 85-89. http://dx.doi.org/10.1590/S1984-29612014012. PMid:24728366.
  » http://dx.doi.org/10.1590/S1984-29612014012
• R Development Core Team. R: A language and environment for statistical computing . Vienna: R foundation for statistical computing; 2015.
• Ramos VN, Piovezan U, Franco AHA, Rodrigues VS, Nava S, Szabó MPJ. Nellore cattle ( Bos indicus) and ticks within the Brazilian Pantanal: ecological relationships. Exp Appl Acarol 2016; 68(2): 227-240. http://dx.doi.org/10.1007/s10493-015-9991-3. PMid:26613758.
  » http://dx.doi.org/10.1007/s10493-015-9991-3
• Rothschild CM. Equine piroplasmosis. J Equine Vet Sci 2013; 33(7): 497-508. http://dx.doi.org/10.1016/j.jevs.2013.03.189.
  » http://dx.doi.org/10.1016/j.jevs.2013.03.189
• Rüegg SR, Heinzmann D, Barbour AD, Torgerson PR. Estimation of the transmission dynamics of Theileria equi and Babesia caballi in horses. Parasitology 2008; 135(05): 555-565. http://dx.doi.org/10.1017/S0031182008004204. PMid:18302805.
  » http://dx.doi.org/10.1017/S0031182008004204
• Sambrook JJ, Russel DW. Molecular cloning: a laboratory manual New York: Cold Spring Harbor Laboratory Press; 2001. 999 p.
• Santos FCC, Vogel FSF, Monteiro SG, Roll VFB. Efeito da imersão em água destilada sobre os parâmetros reprodutivos de teleóginas de Rhipicephalus (Boophilus) microplus. Semina: Ciênc Agrár 2012; 33(3): 1127-1132. http://dx.doi.org/10.5433/1679-0359.2012v33n3p1127.
  » http://dx.doi.org/10.5433/1679-0359.2012v33n3p1127
• Santos TM, Roier EC, Santos HA, Pires MS, Vilela JA, Moraes LM, et al. Factors associated to Theileria equi in equids of two microregions from Rio de Janeiro, Brazil. Rev Bras Parasitol Vet 2011; 20(3): 235-241. http://dx.doi.org/10.1590/S1984-29612011000300011. PMid:21961755.
  » http://dx.doi.org/10.1590/S1984-29612011000300011
• Schaeffer RL, Mendenhall W, Ott L, Gerow KG. Elementary survey sampling . Boston: Cengage Learning Press; 2011.
• Scoles GA, Ueti MW. Vector ecology of equine piroplasmosis. Annu Rev Entomol 2015; 60(1): 561-580. http://dx.doi.org/10.1146/annurev-ento-010814-021110. PMid:25564746.
  » http://dx.doi.org/10.1146/annurev-ento-010814-021110
• Swofford DL. PAUP: Phylogenetic analysis using parsimony. Beta Version 4.0b10. Sunderland: Sinauer and Associates; 2002.
• Thompson JD, Gibson TJ, Plewniak F, Jeanmougin F, Higgins DG. The CLUSTAL_X windows interface: flexible strategies for multiple sequence alignment aided by quality analysis tools. Nucleic Acids Res 1997; 25(24): 4876-4882. http://dx.doi.org/10.1093/nar/25.24.4876. PMid:9396791.
  » http://dx.doi.org/10.1093/nar/25.24.4876
• Torres AJ, Finger IS, Farias NAR, Nizoli LQ, Silva SS, Nogueira CEW. Aspectos epidemiológicos da theileriose equina e sua relação com o carrapato Rhipicephalus (Boophilus) microplus em duas propriedades na região da campanha do Rio Grande do Sul - Brasil. Rev Ibero-Latinoam Parasitol 2012; 71(1): 70-77.
• Uilenberg G. Babesia: a historical overview. Vet Parasitol 2006; 138(1-2): 3-10. http://dx.doi.org/10.1016/j.vetpar.2006.01.035. PMid:16513280.
  » http://dx.doi.org/10.1016/j.vetpar.2006.01.035
• Vatcheva KP, Lee M, McCormick JB, Rahbar MH. Multicollinearity in Regression Analyses Conducted In Epidemiologic Studies. Epidemiology (Sunnyvale) 2016; 6(2): 227. http://dx.doi.org/10.4172/2161-1165.1000227. PMid:27274911.
  » http://dx.doi.org/10.4172/2161-1165.1000227
• Vieira MIB, Costa MM, Oliveira MT, Gonçalves LR, André MR, Machado RZ. Serological detection and molecular characterization of piroplasmids in equids in Brazil. Acta Trop 2018; 179: 81-87. http://dx.doi.org/10.1016/j.actatropica.2017.12.028. PMid:29291385.
  » http://dx.doi.org/10.1016/j.actatropica.2017.12.028
• Vieira TSWJ, Vieira RFC, Finger MAP, Nascimento DAG, Sicupira PML, Dutra LH, et al. Seroepidemiological survey of Theileria equi and Babesia caballi in horses from a rural and from urban areas of Paraná State, southern Brazil. Ticks Tick Borne Dis 2013; 4(6): 537-541. http://dx.doi.org/10.1016/j.ttbdis.2013.07.005. PMid:24209493.
  » http://dx.doi.org/10.1016/j.ttbdis.2013.07.005
• Wise LN, Kappmeyer LS, Mealey RH, Knowles DP. Review of equine piroplasmosis. J Vet Intern Med 2013; 27(6): 1334-1346. http://dx.doi.org/10.1111/jvim.12168. PMid:24033559.
  » http://dx.doi.org/10.1111/jvim.12168
• Wise LN, Pelzel-McCluskey AM, Mealey RH, Knowles DP. Equine piroplasmosis. Vet Clin North Am Equine Pract 2014; 30(3): 677-693. http://dx.doi.org/10.1016/j.cveq.2014.08.008. PMid:25300637.
  » http://dx.doi.org/10.1016/j.cveq.2014.08.008

Publication Dates

• Publication in this collection
  30 Aug 2018
• Date of issue
  Oct-Dec 2018

History

• Received
  04 May 2018
• Accepted
  06 June 2018