Molecular diagnosis and phylogenetic analysis of Babesia bigemina and Babesia bovishemoparasites from cattle in South Africa
© Mtshali and Mtshali; licensee BioMed Central Ltd. 2013
Received: 15 April 2013
Accepted: 2 August 2013
Published: 8 August 2013
Babesia parasites, mainly Babesia bovis and B. bigemina, are tick-borne hemoparasites inducing bovine babesiosis in cattle globally. The clinical signs of the disease include, among others, anemia, fever and hemoglobinuria. Babesiosis is known to occur in tropical and subtropical regions of the world. In this study, we aim to provide information about the occurrence and phylogenetic relationship of B. bigemina and B. bovis species in cattle from different locations in nine provinces of South Africa.
A total of 430 blood samples were randomly collected from apparently healthy cattle. These samples were genetically tested for Babesia parasitic infections using nested PCR assays with species-specific primers.
Nested PCR assays with Group I primer sets revealed that the overall prevalence of B. bigemina and B. bovis in all bovine samples tested was 64.7% (95% CI = 60.0-69.0) and 35.1% (95% CI = 30.6-39.8), respectively. Only 117/430 (27.2%) animals had a mixed infection. The highest prevalence of 87.5% (95% CI = 77.2-93.5) for B. bigemina was recorded in the Free State province collection sites (Ficksburg, Philippolis and Botshabelo), while North West collection sites had the highest number of animals infected with B. bovis (65.5%; 95% CI = 52.7-76.4). Phylograms were inferred based on B. bigemina-specific gp45 and B. bovis-specific rap-1 nucleotide sequences obtained with Group II nested PCR primers. Phylogenetic analysis of gp45 sequences revealed significant differences in the genotypes of B. bigemina isolates investigated, including those of strains published in GenBank. On the other hand, a phylogeny based on B. bovis rap-1 sequences indicated a similar trend of clustering among the sequences of B. bovis isolates investigated in this study.
This study demonstrates the occurrence of Babesia parasites in cattle from different provinces of South Africa. It was also noted that the situation of Babesia parasitic infection in cattle from certain areas within the surveyed provinces had either reached endemic stability or was progressing towards stability.
KeywordsCattle Babesia bigemina Babesia bovis South Africa Nested PCR Phylogeny
Bovine babesiosis (also known as redwater) is one of the most economically important infectious diseases affecting cattle worldwide. This tick-borne disease is mainly induced by intraerythrocytic protozoan parasites of the genus Babesia, order Piroplasmida and phylum Apicomplexa . In southern Africa, Babesia bigemina and Babesia bovis are the two economically important species infecting cattle and have high prevalence in tropical and subtropical regions of the world . Clinical signs characterizing the disease caused by the two parasites are anemia, fever, hemoglobinuria, and in many cases death [3, 4].
Traditionally, the microscopic detection of Babesia parasites has always been considered as the gold standard for the diagnosis of acute babesiosis . However, the major drawback associated with microscopic examination of blood parasites is the low sensitivity offered by the technique, thus making it difficult to detect parasites in blood smears during low parasitemia in the case of carrier animals . Alternatively, PCR-based assays have been widely used for the detection of Babesia parasites owing to their high specificity and sensitivity [4, 6–9].
At present, there are still limited tools available for the prevention and eradication of bovine babesiosis. By far, the most recommended approach for controlling ticks and tick-borne diseases in South Africa is by integrating the strategic use of acaricides (i.e. the pesticides used for killing ticks) and application of vaccines . However, more than 90% of resource-poor farmers contend that the dip wash is not effective in killing the ticks. Consequently, these farmers complement the government dipping service with their own initiatives, which include spraying with conventional acaricides, using household disinfectants and manual removal [11, 12].
It is estimated that, in South Africa, approximately 18% of all cattle mortalities are due to tick-borne diseases including babesiosis, anaplasmosis and heartwater . These diseases have a considerable impact on the country’s economic security and also impact negatively on poor communities who are dependent on livestock production as their source of income and nutritional needs (meat and milk), and as labour for fieldwork and transport . As such, it is reported that the occurrence of B. bigemina and B. bovis in South African cattle population hampers the development of livestock industry, which accounts for up to 49% of the agricultural output .
Given the importance of livestock production in the South African economic landscape, in this study, we genetically investigated the occurrence and distribution of Babesia parasites, more specifically B. bigemina and B. bovis, in bovine samples randomly collected from different locations throughout South Africa. We also studied the phylogenetic relatedness amongst DNA sequences of randomly selected bovine samples.
Blood sample collection
Experimental collection of blood samples from cattle was approved by the NZG Ethics and Scientific Committee, National Zoological Gardens of South Africa. Between 2010 and 2012, a total of 430 blood samples were randomly collected from clinically healthy cattle occupying different locations in all nine provinces of South Africa: Mpumalanga province (Ehlanzeni South District; n = 48), KwaZulu-Natal province (Albert Falls, Ndaleni dip tank and Shallow Drift; n = 52), Limpopo province (Capricorn District; n = 47), North West province (Maubane and Mmatlhwaela; n = 58), Gauteng province (Bronkhorstspruit; n = 30), Free State province (Ficksburg, Philippolis and Botshabelo; n = 64), Eastern Cape province (Alice, Fort Beaufort and Adelaide; n = 60), Northern Cape province (Kuruman; n = 45) and Western Cape province (Boland area; n = 26). The total number of animals sampled in each collection site depended on the number of cattle present at the sampling stations (dip tanks/cattle farms). No information on the age groups, husbandry practices, vaccination histories and tick infestation status of the sampled cattle were available, given that other bovine blood samples were collected and provided by the farmers. Blood was collected from the coccygeal vein into EDTA-coated vacutainer tubes, transported to the laboratory on ice and stored at −20°C until further analysis.
Genomic DNA was extracted from 200 μl of blood using ZR Genomic DNA™-Tissue MiniPrep kit (Inqaba Biotechnical Industries, Pretoria, South Africa) according to the manufacturer’s instructions. Extracted DNA was eluted in 50 μl of DNA elution buffer and stored at −20°C until further analysis. DNA concentration was determined using a NanoDrop® ND-1000 (NanoDrop Technologies Inc., Wilmington, USA).
Primary and nested PCR primers used for PCR amplifications
Primer sequence (5′ → 3′)
Oligonucleotides used for screening (Group I) a
Oligonucleotides used for phylogenetic study (Group II) b
The new pairs of PCR and nested PCR primers were designed from conserved regions identified after performing multiple sequence alignments with CLUSTAL W algorithm  embedded in BioEdit software . The specificity of newly designed primers was tested against GenBank sequences using BLAST search . All primers were synthesized by Inqaba Biotechnical Industries.
Specificity of nested PCR
Purified DNA samples of B. bigemina, B. bovis, Anaplasma centrale, Theileria parva and Ehrlichia ruminantium were used to assess the specificity of Group I and Group II primers. These DNA samples were kindly provided by Dr Nicola Collins (Department of Veterinary Tropical Diseases, University of Pretoria, South Africa) and Dr Oriel Thekisoe (Parasitology Research Program, University of the Free State, South Africa). PCR and nested PCR mixtures were prepared and thermally cycled as described above using 2 μl of each of the purified DNA samples.
PCR detection assays
To detect the presence of B. bigemina and B. bovis from bovine samples using Group I primers, PCR was performed in a 25-μl reaction mixture containing 5 μl of the extracted DNA, 0.6 μM of each primer and 12.5 μl of DreamTaq Green PCR Master Mix (Inqaba Biotechnical Industries). Negative control reactions contained distilled water instead of template DNA. Reaction mixtures were subjected to PCR using BIO-RAD T100™ Thermal Cycler (Bio-Rad Laboratories, Johannesburg, South Africa). PCR amplifications (round 1) were performed at the following thermal conditions: 95°C for 3 min, followed by 35 cycles of 95°C for 30 sec, 55°C for 45 sec and 72°C for 1 min. Following the final extension step at 72°C for 10 min, 1 μl of each PCR product was added into the second (nested) PCR mixture comprising similar composition of reagents as the first round PCR, except that the external primers were replaced with the nested PCR primers. PCR mixtures were cycled as described above using annealing temperatures reflected in Table 1.
PCR-generated amplicons were analyzed by electrophoresis in 1.5% agarose gels stained with Biotium GelRed Acid Stain (Anatech Instruments, Johannesburg, South Africa) and visualized under UV illumination. GeneRuler™ 1 kb Plus DNA ladder (Inqaba Biotechnical Industries) was used as the standard molecular weight marker.
Sequencing and phylogenetic analysis
For phylogenetic study, genomic DNA of randomly selected field samples was amplified with Group II nested PCR primers (Table 1) targeting gp45 and rap-1 fragments of B. bigemina and B. bovis, respectively. Seven bovine samples were selected for B. bigemina and six for B. bovis. PCR mixtures were prepared and thermally cycled as described above employing the annealing temperatures reflected in Table 1. PCR-generated fragments of 853 bp (for B. bigemina) and 1009 bp (for B. bovis) were sent to Inqaba Biotechnical Industries for purification and sequencing in both directions using ABI 3130 XL Genetic Analyzer (Applied Biosystems, Johannesburg, South Africa). At least two individually amplified DNA fragments of each selected sample were sequenced.
Nucleotide sequences were assembled and aligned using BioEdit software program . The resulting consensus sequences were subsequently used to search for homologous sequences in GenBank. To construct phylogenetic trees, the consensus nucleotide sequences were trimmed manually to equivalent lengths. Phylogenies were inferred using the neighbour-joining algorithm of the MEGA v4.1 software . Bootstrapping analysis with 1000 replicates was used to estimate the robustness of individual branches .
Nucleotide sequence accession numbers
The determined B. bigemina gp45 and B. bovis rap-1 gene sequences were submitted to GenBank database under the accession numbers KC894392-KC894404.
The proportions for 95% confidence intervals (95% CI) were computed as CIs for proportions with binomial data employing no continuity correction.
Results and discussion
Nested PCR detection results of B. bigemina and B. bovis parasites in field blood samples
Total number of samples
The animals infected with B. bigemina and B. bovis were found in all the provinces surveyed, albeit there were differences observed in the distribution of these hemoparasites in cattle from certain provinces. The occurrence of B. bigemina and B. bovis in cattle from all surveyed areas could be attributed to the presence and distribution of tick vectors transmitting these parasites. Nested PCR assays detected mostly B. bigemina compared to B. bovis, and these results are similar to previous findings reported elsewhere [20–22]. However, a notable exception was evident for blood samples collected from Kuruman in the Northern Cape province where the occurrence of both Babesia species was below 5%. This low prevalence is not unusual for an area that is considered free of babesiosis tick vectors.
The overall prevalence of B. bigemina and B. bovis in all samples was 64.7% (95% CI = 60.0-69.0) and 35.1% (95% CI = 30.6-39.8), respectively. Samples from the Free State province had the highest frequency of B. bigemina (87.5%, 95% CI = 77.2-93.5), followed by KwaZulu-Natal province with as high as 84.6% (95% CI = 72.5-92.0) of samples tested positive for B. bigemina. These high prevalence values (>80%) were indicative of an endemically stable situation . By definition, an endemically stable situation occurs when 81–100% of the herd are infected with a particular Babesia species . The highest number of samples positive for B. bovis was recorded in the North West province collection sites (65.5%, 95% CI = 52.7-76.4), thus indicating the situation approaching endemic stability but with the potential occurrence of disease outbreaks . Of all the collection sites surveyed, Kuruman area in Northern Cape province had the lowest prevalence of 4.4% (2 out of 45; 95% CI = 1.2-14.8) for both B. bigemina and B. bovis. The cattle found to contain DNA of both B. bigemina and B. bovis (117/430, 27.2%) originated from all but one of the nine provinces from which blood samples were collected. None of the blood samples collected from Kuruman (Northern Cape province) possessed mixed infections.
In principle, the occurrence of B. bigemina and B. bovis in cattle is largely dependent on the distribution of tick vectors . The sole vector of B. bovis in South Africa is Rhiphicephalus (Boophilus) microplus, whereas B. bigemina is transmitted by R. (B.) microplus, Rhiphicephalus (Boophilus) decoloratus and Rhiphicephalus evertsi evertsi . Based on the geographical distribution of tick species in South Africa, R. (B.) microplus is more prevalent in KwaZulu-Natal, Eastern Cape and Limpopo provinces, while the high prevalence of R. (B.) decoloratus was recorded in Mpumalanga, Eastern Cape and KwaZulu-Natal provinces. On the other hand, R. evertsi evertsi is reported to be present throughout South Africa, except in arid areas of the Northern Cape and northern Western Cape provinces [25–29].
Therefore, it is not surprising that in the present study, there was a high percentage of bovine samples positive for B. bigemina and B. bovis infections. These results confirm the previous findings that Babesia infections are common in cattle in South Africa [14, 30, 31]. Using serology-based assays, a recent study  demonstrated the existence of B. bigemina and B. bovis parasites in cattle from eight South African provinces surveyed. Surprisingly, more than 20% of animals from Northern Cape province were found to possess B. bigemina and B. bovis when tested with IFAT and ELISA , and these values are considered higher for an area known to be free of R. (B.) microplus, R. (B.) decoloratus and R. evertsi evertsi tick vectors. According to the latter authors , the higher prevalence values could be explained by the possible outsourcing of animals from endemic areas in the neighbouring provinces.
To study the phylogenetic relationship between Babesia parasites of randomly selected bovine samples, new sets of nested PCR primers (Group II) were designed based on nucleotide sequences of the gp45 (B. bigemina-specific) and rap-1 (B. bovis-specific) genes. The specificity of both Group I and Group II primers was tested against purified DNA samples of B. bigemina, B. bovis, A. centrale, T. parva and E. ruminantium. As expected, the nested PCR assay with B. bigemina gene-specific primers only detected the DNA of B. bigemina. Likewise, B. bovis-specific nested PCR assay specifically identified B. bovis DNA sample, and no amplifications were observed for DNA samples derived from B. bigemina, A. marginale, T. parva and E. ruminantium.
The results presented in this study demonstrate a high incidence of cattle infection by B. bigemina and B. bovis in sampling areas from all provinces surveyed, with a notable exception of Kuruman area (Northern Cape province) that displayed a very low occurrence of Babesia parasites. The high prevalences (>80%) of B. bigemina in cattle from KwaZulu-Natal, Limpopo, Gauteng and Free State sampling areas suggest that the situation of B. bigemina infection has reached endemic stability. Conversely, the prevalence of B. bovis in bovine samples from KwaZulu-Natal, North West and Gauteng collection sites suggests that the situation in these areas is progressing towards endemic stability. Therefore, in order to attain endemic stability to Babesia parasites, a limited number of tick vectors should be allowed to infest cattle, particularly in the case of herds occupying areas with lower infection rates.
In addition, this study has expanded our current knowledge concerning the genetic diversity and phylogenetic relatedness among B. bigemina and B. bovis isolates of South African origin. Given the limited number of Babesia-specific gene sequences available in GenBank, particularly those of Babesia species originating from South African cattle, further studies incorporating sampling sites representative of each surveyed province would be required. Undoubtedly, this will enable a better understanding of the epidemiology of bovine babesiosis as well as the degree of genetic heterogeneities among B. bigemina and B. bovis isolates in South African cattle. Overall, the findings from this study will ultimately help farmers develop efficient control strategies to curtail cattle mortalities emanating from bovine babesiosis.
The authors are indebted to Dr Nicola Collins and Dr Oriel Thekisoe who kindly provided the purified DNA samples of B. bigemina, B. bovis, A. centrale, T. parva and E. ruminantium. The authors are also grateful to A. Tsotetsi, S. Ntanta, M. Legodi, A.M. Mutshembele, L. Seseni and I.D.B. Sekwati for their technical assistance. The farmers who participated in bovine sample collection are gratefully acknowledged. This work is based on the research supported by the National Research Foundation (NRF) and the National Zoological Gardens of South Africa (NZG). Any opinion, finding and conclusion or recommendation expressed in this material is that of the authors and the NRF does not accept any liability in this regard.
- McCosker PJ: The global importance of babesiosis. Babesiosis. Edited by: Ristic M, Kreier JP. New York: Academic; 1981:1-24.Google Scholar
- Bock R, Jackson L, de Vos A, Jorgensen W: Babesiosis of cattle. Parasitol. 2004, 129: 247-269. 10.1017/S0031182004005190.View ArticleGoogle Scholar
- de Vos AJ, Potgieter FT: Bovine babesiosis. Infectious diseases of livestock. Edited by: Coetzer JAW, Thomson GR, Tustin RC. 1994, Cape Town: Oxford University Press, 278-294. 1994Google Scholar
- Martins TM, Pedro OC, Caldeira RA, do Rosário VE, Neves L, Domingos A: Detection of bovine babesiosis in Mozambique by a novel seminested hot-start PCR method. Vet Parasitol. 2008, 153: 225-230. 10.1016/j.vetpar.2008.01.037.PubMedView ArticleGoogle Scholar
- Böse R, Jorgensen WK, Dalgliesh RJ, Friedhoff KT, de Vos AJ: Current state and future trends in the diagnosis of babesiosis. Vet Parasitol. 1995, 57: 61-74. 10.1016/0304-4017(94)03111-9.PubMedView ArticleGoogle Scholar
- Almeria S, Castella J, Ferrer D, Ortuno A, Estrada-Peña A, Gutierrez JF: Bovine piroplasms in Minorca (Balearic Islands, Spain): a comparison of PCR-based and light microscopy detection. Vet Parasitol. 2001, 99: 249-259. 10.1016/S0304-4017(01)00464-2.PubMedView ArticleGoogle Scholar
- Fahrimal Y, Goff WL, Jasmer DP: Detection of Babesia bovis carrier cattle by using polymerase chain reaction amplification of parasite DNA. J Clin Microbiol. 1992, 30: 1374-1379.PubMedPubMed CentralGoogle Scholar
- Figueroa JV, Chieves LP, Johnson GS, Buening GM: Multiplex polymerase chain reaction based assay for the detection of Babesia bigemina, Babesia bovis and Anaplasma marginale DNA in bovine blood. Vet Parasitol. 1993, 50: 69-81. 10.1016/0304-4017(93)90008-B.PubMedView ArticleGoogle Scholar
- AbouLaila M, Yokoyama N, Igarashi I: Development and evaluation of two nested PCR assays for the detection of Babesia bovis from cattle blood. Vet Parasitol. 2010, 172: 65-70. 10.1016/j.vetpar.2010.04.011.PubMedView ArticleGoogle Scholar
- De Waal DT: Anaplasmosis control and diagnosis in South Africa. Ann N Y Acad Sci. 2000, 916: 474-483.View ArticleGoogle Scholar
- Masika PJ, Sonandi A, van Averbeke W: Tick control by small-scale cattle farmers in the central Eastern Cape Province, South Africa. J S Afr Vet Assoc. 1997, 68: 45-48.PubMedGoogle Scholar
- Moyo B, Masika PJ: Tick control methods used by resource-limited farmers and the effect of ticks on cattle in rural areas of the Eastern Cape Province, South Africa. Trop Anim Health Prod. 2009, 41: 517-523. 10.1007/s11250-008-9216-4.PubMedView ArticleGoogle Scholar
- Perry BD, Sones KR: Poverty reduction through animal health. Science. 2007, 315: 333-334. 10.1126/science.1138614.PubMedView ArticleGoogle Scholar
- Terkawi MA, Thekisoe OMM, Katsande C, Latif AA, Mans BJ, Matthee O, Mkize N, Mabogoane N, Marais F, Yokoyama N, Xuan X, Igarashi I: Serological survey of Babesia bovis and Babesia bigemina in cattle in South Africa. Vet Parasitol. 2011, 182: 337-342. 10.1016/j.vetpar.2011.05.047.PubMedView ArticleGoogle Scholar
- Thompson JD, Higgins DG, Gibson TJ: CLUSTAL W: improving the sensitivity of progressive multiple sequence alignment through sequence weighting, position-specific gap penalties and weight matrix choice. Nucleic Acids Res. 1994, 22: 4673-4680. 10.1093/nar/22.22.4673.PubMedPubMed CentralView ArticleGoogle Scholar
- Hall TA: BioEdit: a user-friendly biological sequence alignment editor and analysis program for Windows 95/98/NT. Nucleic Acids Symp Ser. 1999, 41: 95-98.Google Scholar
- Altschul SF, Gish W, Miller W, Myers EW, Lipman DJ: Basic local alignment search tool. J Mol Biol. 1990, 215: 403-410.PubMedView ArticleGoogle Scholar
- Kumar S, Nei M, Dudley J, Tamura K: MEGA: a biologist-centric software for evolutionary analysis of DNA and protein sequences. Brief Bioinform. 2008, 9: 299-306. 10.1093/bib/bbn017.PubMedPubMed CentralView ArticleGoogle Scholar
- Felsenstein J: Confidence limits on phylogenies: an approach using the bootstrap. Evolu. 1985, 39: 783-791. 10.2307/2408678.View ArticleGoogle Scholar
- Applewhaite LM, Craig TM, Wagner GG: Serological prevalence of bovine babesiosis in Guyana. Trop Anim Health Prod. 1981, 13: 13-18. 10.1007/BF02237880.PubMedView ArticleGoogle Scholar
- Miller DK, Diall O, Craig TM, Wagner GG: Serological prevalence of bovine babesiosis in Mali. Trop Anim Health Prod. 1984, 16: 71-77. 10.1007/BF02239848.PubMedView ArticleGoogle Scholar
- Jongejan F, Perry BD, Moorhouse PD, Musisi FL, Pegram RG, Snacken M: Epidemiology of bovine babesiosis and anaplasmosis in Zambia. Trop Anim Health Prod. 1988, 20: 234-242. 10.1007/BF02239989.PubMedView ArticleGoogle Scholar
- Norval RAI, Fivaz BH, Lawrence JA, Dailecourt T: Epidemilogy of tick-borne diseases of cattle in Zimbabwe: I. Babesiosis. Trop Anim Health Prod. 1983, 15: 87-94. 10.1007/BF02239802.PubMedView ArticleGoogle Scholar
- de Vos AJ: Epidemiology of bovine babesiosis in South Africa. Onderstepoort J Vet Res. 1979, 50: 357-362.Google Scholar
- Baker MK, Ducasse FBW, Sutherst RW, Maywald GF: The seasonal tick populations on traditional and commercial cattle grazed at four altitudes in Natal. J S Afr Vet Assoc. 1989, 60: 95-101.PubMedGoogle Scholar
- Fourie LJ, Kok DJ, Heyne H: Adult ixodid ticks on two cattle breeds in the south-western Free State, and their seasonal dynamics. Onderstepoort J Vet Res. 1996, 63: 19-23.PubMedGoogle Scholar
- Bryson NR, Tice GA, Horak IG, Stewart CG, Du Plessis BJA: Ixodid ticks on cattle belonging to small-scale farmers at 4 communal grazing areas in South Africa. J S Afr Vet Assoc. 2002, 73: 98-103.PubMedGoogle Scholar
- Tønnesen MH, Penzhorn BL, Bryson NR, Stoltsz WH, Masibigiri T: Displacement of Boophilus decoloratus by Boophilus microplus in the Soutpansberg region, Limpopo Province, South Africa. Exp Appl Acarol. 2004, 32: 199-208.PubMedView ArticleGoogle Scholar
- Marufu MC, Chimonyo M, Mapiye C, Dzama K: Tick loads in cattle raised on sweet and sour rangelands in the low-input farming areas of South Africa. Trop Anim Health Prod. 2011, 43: 307-313. 10.1007/s11250-010-9690-3.PubMedView ArticleGoogle Scholar
- Dreyer K, Fourie LJ, Kok DJ: Epidemiology of tick-borne diseases of cattle in Botshabelo and Thaba Nchu in the Free State Province. Onderstepoort J Vet Res. 1998, 65: 285-289.PubMedGoogle Scholar
- Tønnesen MH, Penzhorn BL, Bryson NR, Stoltsz WH, Masibigiri T: Seroprevalence of Babesia bovis and Babesia bigemina in cattle in the Soutpansberg region, Limpopo Province, South Africa, associated with changes in vector-tick populations. J S Afr Vet Assoc. 2006, 77: 61-65.PubMedView ArticleGoogle Scholar
- Fisher TG, McElwain TF, Palmer GH: Molecular basis for variable expression of merozoite surface antigen gp45 among American isolates of Babesia bigemina. Infect Immun. 2001, 69: 3782-3790. 10.1128/IAI.69.6.3782-3790.2001.PubMedPubMed CentralView ArticleGoogle Scholar
- Silva MG, Henriques G, Sánchez C, Marques PX, Suarez CE, Oova A: First survey for Babesia bovis and Babesia bigemina infection in cattle from Central and Southern regions of Portugal using serological and DNA detection methods. Vet Parasitol. 2009, 166: 66-72. 10.1016/j.vetpar.2009.07.031.PubMedView ArticleGoogle Scholar
- Ramos CA, Araújo FR, Alves LC, de Souza II, Guedes DS, Soares CO: Genetic conservation of potentially immunogenic proteins among Brazilian isolates of Babesia bovis. Vet Parasitol. 2012, 187: 548-552. 10.1016/j.vetpar.2012.01.020.PubMedView ArticleGoogle Scholar