Assessment of exposure to piroplasms in sheep grazing in communal mountain pastures by using a multiplex DNA bead-based suspension array
© Ros-García et al.; licensee BioMed Central Ltd. 2013
Received: 2 July 2013
Accepted: 17 September 2013
Published: 24 September 2013
Piroplasms are tick-borne hemoprotozoans with a major impact on extensive management systems. Detection of sub-clinical low-level carriers, which can act as source of infection for vector ticks, is key to protect livestock trade and facilitate preventive control programs. The purpose of this study was to develop a method for the detection of ovine piroplasms and to use it in a field study aimed at investigating piroplasms infection in semi-extensive production systems in the Basque Country (northern Spain).
A DNA bead-based suspension array using the Luminex® xMAP technology that included a generic Theileria-Babesia control probe, 6 species-specific probes, and an internal control probe was developed to detect and identify piroplasms that infect sheep. To monitor piroplasm infection in clinically healthy sheep from 4 flocks that share communal mountain pastures, blood samples were collected during 2 grazing seasons.
Piroplasms were detected in 48% (214/446) of blood samples, nearly half of them (49.1%, 105/214) as mixed infections. Five different piroplasms were identified: Theileria sp. OT3 in 34.8% of the samples, Theileria ovis in 20.9%, and at lower prevalences Babesia motasi (12.3%), Theileria luwenshuni/OT1 (10.5%) and Babesia ovis (6.3%). Despite differences among flocks associated to differences in management, an increasing trend in the incidence of piroplasm infection with increasing age of animals after increased tick exposure was observed. This increment could be attributed to continued re-infection associated with re-exposure to ticks at grazing. Ticks were collected from animals (4 species) and vegetation (8 species), and associations between tick abundance seasonality and risk of infection with the different piroplasms were established.
The multiplex Luminex® xMAP procedure is a rapid and high throughput technique that provided highly specific and sensitive identification of single and mixed piroplasm infections in blood of sheep carriers. This study confirmed a situation of endemic stability for piroplasm infection in the region, where infection is present in the absence of clinical signs, and mountain grazing allows for sufficient inoculation rates to maintain such situation.
KeywordsBabesia Theileria Luminex Suspension microarray Sheep RLB Tick-borne pathogens Piroplasmosis
Theileria and Babesia species are tick-borne haemoprotozoan parasites that infect livestock and wildlife in tropical and sub-tropical regions of the world, including sheep. The impact of piroplasmosis is higher on management systems where animals spend long periods grazing in mountain pastures exposed to tick bites. In the Basque Country (northern Spain) there are approximately 324,000 sheep, of which 90% are Latxa breed, the native dairy sheep of the Basque Country, whose production characteristics have been reported elsewhere . Animal husbandry is semi-extensive; sheep are kept on farmland pastures from winter to early spring for lambing (one lambing per ewe per year) and milking, and on communal mountain pastures otherwise. Milking starts after the lamb(s) are slaughtered or weaned, and milked ewes are normally dried off at the beginning of the summer with total lactation length (from lambing to drying off) of about five months. Most of milk is processed to Idiazabal cheese. A previous study carried out in the region  showed a relatively high prevalence of sub-clinical infections in the sheep population and identified five different piroplasms: Babesia ovis, Babesia motasi, Theileria ovis, Theileria sp. OT1 and Theileria sp. OT3. While B. ovis and B. motasi are considered highly and moderately virulent species, respectively, the Theileria species found are considered to be less virulent or benign [3, 4]. Theileria sp. OT1 shares 99.6% similarity in the 18S rRNA gene with a Theileria sp. recently described in China (Theileria luwenshuni) considered highly virulent . Theileria sp. OT1 is widespread in the Basque Country among healthy sheep and is considered non-pathogenic .
Animals that survive acute infection generally become low-level carriers of the parasites, and could remain persistently infected for years without apparent clinical signs [3, 4]. Ticks feeding on these sub-clinical carriers can become infected. Ticks of the genera Rhipicephalus and Haemaphysalis are considered the vectors for T. ovis, B. ovis and B. motasi[3, 4]. For T. luwenshuni the principal vector in China seems to be Haemaphysalis qinghaiensis, whereas the vector of Theileria sp. OT3 is as yet unknown. Theileria lestoquardi, a highly virulent Theileria not reported in Spain, is transmitted by Hyalomma spp. , a tick genus very rarely found in the Basque Country. In Babesia spp., there is both transovarial and transstadial transmission, whereas in Theileria spp. only transstadial transmission occurs [3, 7]. Also, whereas most Babesia spp. inoculated by ticks invade erythrocytes of the vertebrate host directly, Theileria species first infect lymphocytes leading to clonal expansion of the infected cells before invading erythrocytes .
Carrier animals are known to exhibit fluctuating low parasitaemia which sometimes even escapes detection. This is particularly true for Babesia spp. infection where intraerythrocytic piroplasms in the circulating bloodstream rarely exceed a few percent even during acute infection [9, 10]. Thus, detection of sub-clinical carriers requires sensitive diagnostic tools able to detect and identify the different piroplasms that infect sheep. Reverse Line Blot (RLB) has been widely used for this purpose [2, 11, 12]. However, despite its multiplexing capacity, RLB is long and tedious. Here, we developed a DNA bead-based suspension array test based on the Luminex® xMAP technology to detect and differentiate ovine Babesia and Theileria species. This technique has already been shown to provide higher sensitivity and throughput than RLB . Once developed, we used it to monitor piroplasms infection among clinically healthy sheep grazing in a communal mountain pasture in a region where sub-clinical piroplasm infection was detected in the past . Ticks from animals and vegetation were also collected to study piroplasms in relation with tick abundance and distribution.
Animal selection and sample collection
The study was carried out in a communal mountain pasture (700–1,200 m above sea level) located in Alava, a province located in the Basque Country, northern Spain. Climate is of transitional type, resulting from the interaction of an Atlantic climate and a Continental Mediterranean climate, with cold winters and mild summers. Mean annual rainfall is ca. 900 mm. Vegetation is mainly composed of beech trees (Fagus sylvatica) with areas of pines (Pinus spp.) and small patches of oaks (Quercus spp.) interspersed with patches of natural grass.
Four flocks of Latxa breed (F1, F2, F3 and F4) that graze in these communal mountain pastures were monitored for piroplasms infections: flocks F1 & F2 during two grazing seasons (2011–2012), flock F3 solely in 2011, and flock F4 in 2012. There were some differences among flocks regarding milk production and management. Thus, in F1 and F2 management was more intensive and milk production was clearly higher than in F3 and F4. The use of communal mountain pastures was also different among flocks. In this way, in F1 and F2 only a few animals (non-milking ewes) grazed in these pastures from April to middle of July, whereas in F3 and F4 a higher number of animals (non-milking ewes, replacement lambs and their dams) were kept on mountains in this period of the year. From middle of July until November all animals of the four flocks grazed in mountain pastures. While grazing in the mountains, no feed supplementation was given and feeding relied only on pasture availability. No specific treatments against ticks were carried out in any of the flocks. Size of flocks ranged from 120 to 599 animals.
All flocks were sampled three times per year, before going to the communal pastures (spring), and twice during the grazing season (summer and autumn). Blood samples were randomly collected from ca. 25 animals per flock and sampling, divided into three age categories (A, <1 year; B, 1–2 years; C, >2 years). Animals were apparently healthy at sampling. Thus, a total of 446 blood samples were collected from clinically healthy sheep aged 2 months - 8 years. During the second year sampled animals were inspected for the presence of ticks, mainly on the ears, along the nape of the neck, perineum, udder and tail base, which were manually removed, counted and identified [14, 15].
Questing ticks sampling
To investigate tick abundance in the same mountain pastures where animals graze, ticks were monthly collected between May 2011 and April 2013 from the vegetation by blanket dragging a 1 m2 white cotton towel over vegetation transects of 100 m, stopping every 10 m to collect and count all attached adult and nymph ticks. For larvae counts, towels were examined at the laboratory and counted; exceptionally, when the amount of collected larvae was too large, an estimate was made after visually dividing the towel into 10 cm stripes and counting three of them (first, middle and last). Ticks were identified using taxonomic keys [14, 15] and stored at −20°C. To compare the abundance of each tick species and stages between samplings, tick count data were converted into tick abundance indexes (TAI) for larvae, nymphs and adults, which express the number of each tick species and stage collected in transects of 100 m2, as TAI = TR x 100 / a, where TR is the number of ticks recorded and a is the sampled area in square meters.
DNA was extracted from 200 μl of ovine blood using the QIAamp DNA Mini Kit (Qiagen, Hilden, Germany), including negative extraction controls every 10 samples. DNA yields were determined with a Nanodrop® ND-1000 Spectrophotometer (Nanodrop Technologies, DE, USA), and DNA was stored at − 20°C until subsequent analysis.
PCR was used to amplify the hypervariable V4 region of the 18S rRNA gene of the genera Babesia and Theileria using primers RLB-F2 (phosphorylated-5’-GACACAGGGAGGTAGTGACAAG-3’) and RLB-R2 (biotinylated-5’-CTAAGAATTTCACCTCTGACAGT-3’) (Sigma–Aldrich, MO, USA) as reported before . An internal amplification control (IAC), which is co-amplified with the same primers as the piroplasms DNA, was also included in each tube .
PCR amplification was performed in a final volume of 50 μl containing 50 ng of genomic DNA, 200 nM of each primer, 1 × PCR buffer, 1.5 mM MgCl2, 200 μM of each dNTP, 10 copies of the internal amplification control (IAC) plasmid and 1U of Taq Platinum Polymerase (Invitrogen, CA,USA). PCR conditions consisted of an enzyme activation step of 4 min at 94°C, and 40 cycles of 30 s at 94°C, 1 min at 51°C and 35 s at 72°C. Extraction controls and PCR negative (water) controls were included in each PCR reaction as negative controls.
Molecular detection and identification of ovine piroplasms by Luminex® xMAP technology
Luminex procedure description
Oligonucleotide probes used in the Luminex suspension array for the detection and species identification of ovine piroplasms
Sequence (5’- 3’)
Theileria sp. OT3
Theileria annulata / lestoquardi
The capture probes were bound to different polystyrene microsphere sets as described elsewhere . Biotinylated PCR products were hybridised to microspheres coupled with the probes mentioned above in 96-well plates in a suspension array. The optimal hybridisation temperature and incubation time proved to be 54°C and 30 min with shaking at 650 rpm in a Thermomix Comfort (Eppendorf). Detection was carried out by incubation at 54°C during 10 min with agitation at 650 rpm using streptavidin-phycoeritryn in 1x TMAC buffer. Finally, the signals produced for each bead were analyzed using the Luminex 200® system (Austin, TX) and expressed as mean fluorescent intensity (MFI) values. The cut-off value for a positive result was calculated for each assay and probe as described elsewhere .
Specificity, sensitivity and test performance
The analytical specificity of the probes was tested against 27 recombinant plasmids constructed as previously described  and containing inserts corresponding to the V4 variable region of the 18S rRNA gene of the following Theileria spp. and Babesia spp. from sheep, cattle, horse, dogs and wildlife: T. ovis, T. lestoquardi, Theileria sp. OT1, Theileria sp. OT3, B. ovis, B. motasi, Theileria annulata, Theileria buffeli, Theileria parva, Babesia bovis, Babesia divergens, Babesia bigemina, Babesia major, Babesia occultans, Theileria equi (genogroups A-D), Babesia caballi (genogroups A and B), Babesia canis, Babesia vogeli, Babesia gibsoni, Theileria annae, Babesia microti, Babesia sp. EU1, Theileria sp. 3185/02.
Analytical sensitivity of the hybridisation assay was assessed by processing serially diluted plasmids (103–1 copies) of the ovine Theileria spp. and Babesia spp. To assess the diagnostic test performance of the Luminex assay, 177 clinical samples from sheep were analysed in parallel by RLB hybridisation, which was performed as described elsewhere . To determine the detection limit of the Luminex assay, blood from an animal with clinical symptoms of piroplasmosis was collected. Haematological examination revealed low values in red blood cells (8.42 x 106 erythrocytes/mm3), low haematocrit (22%) and leukopenia (2.7 x 103 leukocytes/mm3). Giemsa-stained slides were examined under oil immersion using × 100 objective lens and 1.16% of the erythrocytes were parasitized by B. ovis. The number of erythrocytes and intracellular forms with morphology compatible with B. ovis present in 10 microscopic fields were counted. Mean parasitaemia was set at 0.0168 babesias/erythrocyte and parasitaemia estimated at 1.4 x 105 babesias/mm3. Blood was serially diluted (1/10 dilutions until 1.4 babesias/mm3, and 1/2 dilutions thereafter) and processed as described above. The V4 hypervariable region of the last positive dilution was sequenced using the ABI BigDye™ Terminator Cycle Sequencing Ready Reaction Kit and an ABI3130 Genetic Analyzer (Applied Biosystems, Foster City, CA, USA).
The level of agreement between methods (RLB and Luminex performed in parallel on the same samples) was tested by the Kappa (κ) index test at a 95% confidence interval using Win Episcope 2.0. A variable called complementary sensitivity (CSe), which measures the additional detection efficacy of method 1 over method 2 when both have similar specificity, was calculated as 100 × (no. samples positive by method 1 and negative by method 2/total no. samples positive by method 2) .
Differences in piroplasm infection or not (binomial) by each species regarding age [cat. A (<1 yr), cat. B (1–2 yr), cat. C (>2 yr)] (categorical; three levels), flock [F1–F4] (categorical; four levels) and sampling season (categorical; three levels), were assessed by analysis of variance and least square means comparison with the Tukey-Kramer adjustment for multiple comparisons in the GLM procedure of the SAS statistical package version 9.1 (SAS Institute Inc., Cary, NC, USA). P values less than 0.05 were considered statistically significant.
Samples were collected by clinical veterinarians as part of the usual screening scheme on farms and Spanish ethical guidelines and animal welfare regulations (RD 1201/2005) were strictly respected. All herd owners had given an informed consent prior to the study.
Optimization of the Luminex assay
Firstly, the panel of probes designed for ovine piroplasms identification by RLB hybridisation (B. ovis, B. motasi, T. ovis, T. luwenshuni/OT1, Theileria sp. OT3 and T. lestoquardi)  was tested in the Luminex system. Three of them (T. ovis, T. luwenshuni/OT1 and Theileria sp. OT3) showed good specificity and sensitivity, but those for B. ovis, B. motasi and T. lestoquardi produced very low MFI values (low sensitivity) and cross-reactions (lack of specificity). Therefore, new probes were specifically designed for Luminex for B. ovis and B. motasi. For T. lestoquardi, which shares 99.5% identity with T. annulata in the 18S rRNA gene sequence, the probe developed for the Luminex detection of T. annulata was used , and designated T. annulata/T. lestoquardi probe hereafter. The Theileria-Babesia conserved catch-all TB probe used in RLB hybridisation assays had already been tested in a Luminex suspension array .
Analytical sensitivity and specificity and diagnostic sensitivity
The detection limit of the Luminex assay was determined using serial dilutions of infected blood from an animal with a known level of B. ovis parasitaemia. The lowest parasitaemia that could be detected by the Luminex method corresponded to 0.35 parasites per μl of blood. Sequencing analysis of the amplification product obtained confirmed the presence of B. ovis DNA in this last dilution positive by Luminex.
Performance compared with RLB hybridisation
Luminex assay performance for the detection and identification of ovine piroplasms in sheep blood field samples compared with RLB hybridisation
CSe (%) Luminex
0.938 ± 0.027
1.000 ± 0.000
0.684 ± 0.120
0.962 ± 0.022
0.677 ± 0.113
Theileria sp. OT3
0.886 ± 0.035
Piroplasms distribution among clinically healthy sheep grazing in communal mountain pastures
Percentage of samples positive to piroplasms (single or mixed infection) and percentage of samples where each piroplasm species was identified
% piroplasm positive
Species identified (%)
T. luwenshuni/ OT1
Theileria sp. OT3
Ticks collected from the animals and the vegetation
During the second sampling season, a total of 192 ticks (125 adults, 60 nymphs and 7 larvae) were removed from 59 of the 222 animals examined. The mean number of ticks removed per infested animal was 3.5. However, most of the ticks (70%) were collected from a single flock (F4), where 30% of the animals were infested with an average of 6.4 ticks per infested animal. Four different tick species were found and identified as Ixodes ricinus (43.8%), Haemaphysalis punctata (34.9%), Rhipicephalus bursa (15.1%) and Dermacentor marginatus (6.3%). All D. marginatus (100%) and almost all H. punctata (95.5%) collected from the animals were adults, whereas the three tick stages of R. bursa and I. ricinus were collected. Still, R. bursa adults were most commonly collected (65.5%), all of them being collected in summer. Conversely, nymphs accounted for 59.5% of all I. ricinus, and were mainly collected in autumn. H. punctata was also most commonly collected in autumn than in any other season.
Little attention has been given to ovine piroplasmosis compared to bovine piroplasmosis despite its widespread distribution in southern Europe, the Middle East, China, and through tropical and subtropical areas. New techniques with higher sensitivity and large multiplexing capacity for the detection of carrier animals are necessary to better understand the epidemiology of ovine piroplasmosis. The Luminex® xMAP procedure presented here is a high throughput technique that combines large multiplexing capacity with high sensitivity, making it an ideal tool for the detection and identification of single and mixed piroplasm infections in sheep blood. As already shown for cattle , this study demonstrated the benefits of Luminex over RLB, another multiplex technique widely used for piroplasms detection and identification. Small modifications in probe design for B. ovis and B. motasi were needed to transfer the piroplasm RLB assay to the new hybridisation Luminex® xMAP platform. Thus, the new multiplex panel of ovine piroplasm probes included in the Luminex® xMAP assay allowed specific and sensitive detection of serial dilutions of a panel of recombinant plasmids (1 copy of target gene) for the different Theileria and Babesia species. Moreover, the Luminex® xMAP procedure was able to detected piroplasms in samples that were negative by RLB hybridisation. Finally, the inclusion of two control probes, the catch-all TB probe, that detects a wide spectrum of piroplasm species [2, 18], and the IAC probe, that monitors for the presence of PCR inhibitors [13, 19], assures that all new species and genotypes are detected at least at the group level and guarantees that false negatives do not occur.
Once validated, we used the Luminex® xMAP technology to monitor piroplasm infection throughout two years among clinically healthy sheep from four flocks grazing in the same communal mountain pasture in a region where subclinical piroplasm infection had been detected in the past . In this study, the same piroplasm species were detected in the Basque Country, but prevalences were different to those observed previously. While the percentage of animals positive for B. motasi increased considerably in this study, T. luwenshuni/OT1 showed a marked decrease compared to the study carried out in 2004. Then, samples were selected to represent the ovine population of the whole Basque Country (Study I – ); here, only four flocks restricted to a very specific location were sampled. Despite sharing the same mountain communal pastures, differences were found among flocks, probably due to differences in management. In fact, in flocks F1 and F2, milking animals and replacement lambs, which represented ca. 85% of the flock census, remained until July in farmland pastures with scarce presence of ticks, before going to communal mountain pastures; in flocks F3 and F4, 50% of the animals in the flock grazed in mountain pastures from April. These differences in tick exposure periods, along with a more intensive management system in F1 and F2 as compared to F3 and F4 where management is more traditional, would explain differences in piroplasms infection.
Taking into account that lambing in Latxa sheep in Alava occurs in January-February, three age categories (A, <1 year; B, 1–2 years; C, >2 years) were established in such manner that during the first annual sampling (spring), lambs (Cat. A) from the four flocks had not yet been exposed to ticks. Therefore, they were expected to be negative in spring and acquire the infection while the sampling period progressed and exposure to the tick vectors increased. Animals in Cat. B (hoggets) and C (ewes) had been exposed to ticks one or more grazing seasons, respectively, when sampled in spring. The analysis of piroplasm infection considering age at sampling and sampling season, showed increasing piroplasm prevalence with increasing age of animals and tick exposure. The increasing trend observed could be attributed to continued re-infection associated with re-exposure to ticks at grazing. Still, in the case of carrier animals infected with Babesia spp., where an intra-leukocytic phase does not occur, parasitaemia is much lower and sometimes barely perceptible . Hence, Babesia spp.-infected animals were more easily detected after flocks were re-exposed to infected ticks at grazing.
In order to evaluate the evolution of piroplasm infection in relation with tick abundance and distribution, ticks were removed from sampled animals and collected from the vegetation in the areas where sheep graze. The number of ticks collected from the animals was too low to draw many conclusions. Still, the natural vectors of ovine piroplasms, H. punctata and R. bursa, were collected from the animals at some point over the sampling year. The most abundantly collected tick species, both from the animals and the vegetation, was I. ricinus but abundance of questing I. ricinus in this area was lower than that reported in other parts of the Basque Country [20, 21]. This tick species has not been described as vector of any ovine piroplasm species, but DNA of T. ovis, Theileria sp. OT3 and B. ovis has been detected in I. ricinus questing ticks . On the other hand, H. punctata, which is vector of B. motasi and T. ovis, was the second most abundantly collected tick from the animals and showed higher abundance in vegetation than previously reported . Mainly adults were recovered from the animals and were mostly collected in autumn, the season when questing H. punctata adult ticks were also most active. Immature stages captured from the vegetation were more abundant in summer and autumn and decreased thereafter, until the activity of larvae disappeared in winter. In accordance with tick seasonality, summer and autumn would be the seasons of highest risk of infection with B. motasi. Accordingly, results showed that B. motasi infection was most prevalent in summer and autumn (data not shown), which could be a reason for concern considering the previously reported association between B. motasi and sheep abortion .
R. bursa transmits T. ovis and B. ovis[3, 7]. In this study, R. bursa accounted for 15.1% of the ticks collected directly from the animals (adults in summer and, nymphs and larvae in autumn), but only 1 adult was collected by blanket dragging. Other studies carried out in the Basque Country, detected larger proportions of questing R. bursa (2-10% of adult ticks collected) but, same as here, adult R. bursa ticks were all collected at the end of spring and early summer with peaks in June and July [21, 23]. Although B. ovis is described as the most virulent piroplasm species for sheep, here, animals positive to B. ovis did not show any clinical signs of disease at sampling. Incidence of B. ovis infection was particularly high in F4, as was infestation with R. bursa. In this flock, lambs are exposed to ticks for the first time at early spring, before maternal antibodies wane . Therefore, lambs bitten by B. ovis-infected R. bursa adult ticks at the end of spring and early summer would not develop disease but acquire protective immunity. Later, in autumn, when adults of R. bursa are not active, exposure to infected larvae and nymphs, would favour premunition . Due to differences in management, lambs in F1 and F2 start grazing in mountain pastures in July, therefore they have less probability of being bitten by R. bursa adult ticks. This could explain why they remained negative at sampling not only in summer but also in autumn. In the case of T. ovis, which as other Theileria spp., is detectable in blood for long periods, associations between the peak activity of vectors and host infection are more difficult to establish. Besides, T. ovis can be transmitted by both R. bursa and H. punctata. In any case, T. ovis is considered non-pathogenic and distinct clinical signs are rare . The principal vector of T. luwenshuni in China has been reported to be Haemaphysalis qinghaiensis, suggesting that another Haemaphysalis species could be involved in the transmission of T. luwenshuni/OT1 in Spain. In a previous study, Theileria sp. OT1 was detected in 3 different tick species, I. ricinus, H. punctata and D. reticulatus. The vector of Theileria sp. OT3 remains unknown, but a possible association with H. punctata was proposed .
In conclusion, this study confirmed a situation of endemic stability for piroplasm infection in the region, where infection is present in the absence of clinical signs, and mountain grazing allows for sufficient inoculation rates to maintain such a situation . Although mild or sporadic cases of clinical disease might go unnoticed during mountain grazing, when animals are not closely monitored, absence of clinical signs at sampling could be due to the fact that in endemic areas livestock can develop resistance to ticks and piroplasm infection [26, 27]. Future studies on the presence of piroplasms in the collected questing ticks will provide more information on the proportion of infection of the different combinations between piroplasm species and vector tick. For these studies, the Luminex® xMAP technology using probes for ovine piroplasms as described here, probes for piroplasm species that infect cattle as previously reported  and probes for horses currently under development, will facilitate rapid and abundant data in a multiplexing and highly sensitive format of high processability.
We thank Dr Dirk Geysen (Institute of Tropical Medicine, Belgium) for kindly providing T. parva, Prof. Dr. Ulrike Seitzer (Research Center Borstel, Germany) for T. lestoquardi, Prof. Dr. Tatjana Avsic-Zupanc (University of Ljubljana, Slovenia) for B. microti and Babesia sp. EU1 and, Dr. Sonia Almería (Universidad Autónoma de Barcelona, Spain) for T. ovis. The authors would also like to thank the farmers for their collaboration in sample collection. This work was conducted under financial support from the Spanish National Institute for Agricultural and Food Research and Technology (INIA, Project No. RTA2009-000-18-00-00), and the European Regional Development Fund (ERDF). AR is the recipient of a predoctoral fellowship from INIA.
- Gabiña D, Arrese F, Arranz J, Beltran de Heredia I: Average milk yields and environmental effects of Latxa sheep. J Dairy Sci. 1993, 76: 1191-1198. 10.3168/jds.S0022-0302(93)77448-2.View ArticleGoogle Scholar
- Nagore D, García-Sanmartín J, García-Pérez AL, Juste RA, Hurtado A: Identification, genetic diversity and prevalence of Theileria and Babesia species in a sheep population from Northern Spain. Int J Parasitol. 2004, 34: 1059-1067. 10.1016/j.ijpara.2004.05.008.View ArticlePubMedGoogle Scholar
- Preston PM: Theilerioses. Encyclopedia of arthropod-transmitted infections of man and domesticated animals. Edited by: Service M.W. 2001, Wallingford, UK: CABI Publishing, 487-502.Google Scholar
- Uilenberg G: Babesia-A historical overview. Vet Parasitol. 2006, 138: 3-10. 10.1016/j.vetpar.2006.01.035.View ArticlePubMedGoogle Scholar
- Yin H, Luo J, Schnittger L, Lu B, Beyer D, Ma M, Guan G, Bai Q, Lu C, Ahmed J: Phylogenetic analysis of Theileria species transmitted by Haemaphysalis qinghaiensis. Parasitol Res. 2004, 92: 36-42. 10.1007/s00436-003-0900-z.View ArticlePubMedGoogle Scholar
- Yin H, Luo J, Guan G, Lu B, Ma M, Zhang Q, Lu W, Lu C, Ahmed J: Experiments on transmission of an unidentified Theileria sp. to small ruminants with Haemaphysalis qinghaiensis and Hyalomma anatolicum anatolicum. Vet Parasitol. 2002, 108: 21-30. 10.1016/S0304-4017(02)00166-8.View ArticlePubMedGoogle Scholar
- Uilenberg G: Babesiosis. Encyclopedia of arthropod-transmitted infections of man and domesticated animals. Edited by: Service M.W. 2001, Wallingford, UK: CABI Publishing, 53-60.Google Scholar
- Sonenshine DE: Biology of ticks, Volume II. 1993, New York: Oxford University PressGoogle Scholar
- Chauvin A, Moreau E, Bonnet S, Plantard O, Malandrin L: Babesia and its hosts: adaptation to long-lasting interactions as a way to achieve efficient transmission. Vet Res. 2009, 40: 37-10.1051/vetres/2009020.PubMed CentralView ArticlePubMedGoogle Scholar
- Yeruham I, Hadani A, Galker F: Some epizootiological and clinical aspects of ovine babesiosis caused by Babesia ovis–a review. Vet Parasitol. 1998, 74: 153-163. 10.1016/S0304-4017(97)00143-X.View ArticlePubMedGoogle Scholar
- Liu AH, Yin H, Guan GQ, Schnittger L, Liu ZJ, Ma ML, Dang ZS, Liu JL, Ren QY, Bai Q, Ahmed JS, Luo JX: At least two genetically distinct large Babesia species infective to sheep and goats in China. Vet Parasitol. 2007, 147: 246-251. 10.1016/j.vetpar.2007.03.032.View ArticlePubMedGoogle Scholar
- Schnittger L, Yin H, Qi B, Gubbels JM, Beyer D, Niemann S, Jongejan E, Ahmed JS: Simultaneous detection and differentiation of Theileria and Babesia parasites infecting small ruminants by reverse line blotting. Parasitol Res. 2004, 92: 189-196. 10.1007/s00436-003-0980-9.View ArticlePubMedGoogle Scholar
- Ros-García A, Juste RA, Hurtado A: A highly sensitive DNA bead-based suspension array for the detection and species identification of bovine piroplasms. Int J Parasitol. 2012, 42: 207-214. 10.1016/j.ijpara.2011.12.001.View ArticlePubMedGoogle Scholar
- Gil-Collado J, Guillén-Llera JL, Zapatero-Ramos LM: Claves para la identificación de los Ixodoidea españoles (adultos). Rev Iber Parasitol. 1979, 39: 107-111.Google Scholar
- Manilla G: Fauna D’Italia: Acari, Ixodida. 1998, Bologna: Edizioni CalderiniGoogle Scholar
- Georges K, Loria GR, Riili S, Greco A, Caracappa S, Jongejan F, Sparagano O: Detection of haemoparasites in cattle by reverse line blot hybridisation with a note on the distribution of ticks in Sicily. Vet Parasitol. 2001, 99: 273-286. 10.1016/S0304-4017(01)00488-5.View ArticlePubMedGoogle Scholar
- Juste RA, Garrido JM, Geijo M, Elguezabal N, Aduriz G, Atxaerandio R, Sevilla I: Comparison of blood polymerase chain reaction and enzyme-linked immunosorbent assay for detection of Mycobacterium avium subsp. paratuberculosis infection in cattle and sheep. J Vet Diagn Invest. 2005, 17: 354-359. 10.1177/104063870501700409.View ArticlePubMedGoogle Scholar
- Gubbels JM, de Vos AP, van der Weide M, Viseras J, Schouls LM, de Vries E, Jongejan F: Simultaneous detection of bovine Theileria and Babesia species by reverse line blot hybridization. J Clin Microbiol. 1999, 37: 1782-1789.PubMed CentralPubMedGoogle Scholar
- Hoorfar J, Malorny B, Abdulmawjood A, Cook N, Wagner M, Fach P: Practical considerations in design of internal amplification controls for diagnostic PCR assays. J Clin Microbiol. 2004, 42: 1863-1868. 10.1128/JCM.42.5.1863-1868.2004.PubMed CentralView ArticlePubMedGoogle Scholar
- Barandika JF, Hurtado A, Juste RA, García-Pérez AL: Seasonal dynamics of Ixodes ricinus in a 3-year period in northern Spain: first survey on the presence of tick-borne encephalitis virus. Vector Borne Zoonotic Dis. 2010, 10: 1027-1035. 10.1089/vbz.2009.0148.View ArticlePubMedGoogle Scholar
- Barandika JF, Olmeda SA, Casado-Nistal MA, Hurtado A, Juste RA, Valcarcel F, Anda P, García-Pérez AL: Differences in questing tick species distribution between Atlantic and continental climate regions in Spain. J Med Entomol. 2011, 48: 13-19. 10.1603/ME10079.View ArticlePubMedGoogle Scholar
- García-Sanmartín J, Barandika JF, García-Pérez AL, Juste RA, Hurtado A: Distribution and molecular detection of Theileria and Babesia in questing ticks from Northern Spain. Med Vet Entomol. 2008, 22: 318-325. 10.1111/j.1365-2915.2008.00748.x.View ArticlePubMedGoogle Scholar
- Barandika JF, Berriatua E, Barral M, Juste RA, Anda P, García-Pérez AL: Risk factors associated with ixodid tick species distributions in the Basque region in Spain. Med Vet Entomol. 2006, 20: 177-188. 10.1111/j.1365-2915.2006.00619.x.View ArticlePubMedGoogle Scholar
- Yeruham I, Hadani A, Galker F, Rosen S: A study of an enzootic focus of sheep babesiosis (Babesia ovis, Babes, 1892). Vet Parasitol. 1995, 60: 349-354. 10.1016/0304-4017(95)00783-7.View ArticlePubMedGoogle Scholar
- Norval RAI, Perry BD, Young AS: The Epidemiology of Theileriosis in Africa. 1992, London, UK: Academic PressGoogle Scholar
- Dalgliesh RJ:Babesiosis. Immunology and molecular biology of parasitic infections. Edited by: Warren KS. 1993, Boston: Blackwell Scientific Publications, 3,Google Scholar
- Wikel SK: Tick modulation of host immunity: an important factor in pathogen transmission. Int J Parasitol. 1999, 29: 851-859. 10.1016/S0020-7519(99)00042-9.View ArticlePubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.