# Efficient Transovarial Transmission of Babesia Spp. in Rhipicephalus microplus Ticks Fed on Water Buffalo (Bubalus bubalis)

Dasiel Obregon, Belkis Corona-Gonzalez, Adrian Alberto Díaz Sánchez, Yasmani Armas, Eugenio Roque, Márcia Cristina de Sena Oliveira, Alejandro Cabezas Cruz

► To cite this version:

Dasiel Obregon, Belkis Corona-Gonzalez, Adrian Alberto Díaz Sánchez, Yasmani Armas, Eugenio Roque, et al.. Efficient Transovarial Transmission of Babesia Spp. in Rhipicephalus microplus Ticks Fed on Water Buffalo (Bubalus bubalis). Pathogens 2020, 2020, 9, 10.3390/pathogens9040280. hal-02951798

HAL Id: hal-02951798

https://hal.inrae.fr/hal-02951798v1

Submitted on 28 Sep 2020

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons CC BY 4.0 - Attribution - International License
Article

# Efficient Transovarial Transmission of *Babesia* Spp. in *Rhipicephalus microplus* Ticks Fed on Water Buffalo (*Bubalus bubalis*)

Dasiel Obregón ¹,²*, Belkis Corona-González ³, Adrian Alberto Díaz-Sánchez ³,⁴, Yasmani Armas ⁵, Eugenio Roque ⁵, Márcia Cristina de Sena Oliveira ⁶ and Alejandro Cabezas-Cruz ⁷,*

¹ School of Environmental Sciences, University of Guelph, Guelph, ON N1G 2W1, Canada

² Center for Nuclear Energy in Agriculture, University of São Paulo, Piracicaba, São Paulo 13400-970, Brazil

³ Centro Nacional de Sanidad Agropecuaria (CENSA), San José de Las Lajas, Mayabeque 32700, Cuba; bcorona@censa.edu.cu (B.C.-G.); adiasanz88@gmail.com (A.A.D.-S.)

⁴ Department of Biology, University of Saskatchewan, Saskatoon, SK S7N 5A2, Canada

⁵ Faculty of Veterinary Medicine, Universidad Agraria de La Habana, San José de Las Lajas, Mayabeque 32700, Cuba; yarmas2012@gmail.com (Y.A.); roque@unah.edu.cu (E.R.)

⁶ Embrapa Pecuária Sudeste, São Carlos, São Paulo 13560-970, Brazil; marcia.sena-oliveira@embrapa.br

⁷ UMR BIPAR, INRAE, ANSES, Ecole Nationale Vétérinaire d'Alfort, Université Paris-Est, 94700 Maisons-Alfort, France

* Correspondence: dasieloa@uoguelph.ca (D.O.); alejandro.cabezas@vet-alfort.fr (A.C.-C.)

Received: 10 March 2020; Accepted: 10 April 2020; Published: 11 April 2020

**Abstract:** Water buffaloes can be infected by tick-borne pathogens (TBPs) in endemic areas where cattle and buffalo coexist. Among TBPs affecting buffaloes is the Apicomplexan hemoparasites *Babesia bovis* and *B. bigemina*, transmitted by *Rhipicephalus microplus* ticks. However, little empirical evidence exists on whether buffalo can support TBPs' infection and transmission. A cohort study was designed to measure the infestation levels of *R. microplus* in buffaloes as well as the ability of buffalo-fed ticks to transmit *B. bovis* and *B. bigemina* to their offspring. Tick infestation of different life stages was quantified in cattle and buffalo kept in field conditions in western Cuba. Engorged adult female ticks were allowed to lay eggs in controlled conditions of humidity and temperature, and reproductive parameters were measured and analyzed. Hosts and tick larvae were tested for the presence of *Babesia* spp. using species-specific qPCR assays. Tick infestation was not observed in adult buffaloes. However, buffalo and cattle calves were equally infested, although the larval survival rate was higher in cattle calves than in buffalo calves. All larval pools (31) obtained from the adult female ticks were positive for *B. bovis*, whereas only 68% (21/31) was positive for *B. bigemina*. Among the 10 larval pools negative for *B. bigemina*, three proceeded from adult females fed on *Babesia*-negative buffaloes. The other seven pools were from *Babesia*-positive animals, three from cattle and four from buffalo calves. *Babesia* infection levels in tick larvae, quantified by qPCR, were similar in female ticks fed on buffalo and bovine calves. We conclude that water buffalo can sustain tick vector populations and support *Babesia* infection in levels high enough as to be infective for ticks. Our results also validated the hypothesis that adult female ticks fed on buffalo can transmit the pathogens *B. bovis* and *B. bigemina* to their offspring. Nevertheless, further laboratory studies are needed to address the question of whether the transovarial transmission of *Babesia* occurs in the following settings: (1) When adult females are infected previous to the feeding on the buffalo or/and (2) when the adult females acquire the infection while feeding on the buffalo.

**Keywords:** water buffalo; cattle; ticks; reproductive efficiency; *Babesia*; qPCR
# 1. Introduction

Bovine babesiosis, caused by the protozoan parasites *Babesia bovis*, *B. bigemina*, and *B. divergens* (Apicomplexa: Babesiidae), is a tick-borne disease (TBD) posing a constraint for livestock production in tropical and subtropical regions of the world. Specifically, *B. divergens* has been reported only in northeast Europe [1,2], whereas *B. bovis* and *B. bigemina* are widespread in several countries in Latin America, Africa, Australia, and Asia. The one-host ixodid tick *Rhipicephalus microplus* is the vector of these pathogens in the Caribbean and South American countries [3,4]. These piroplasms replicate within tick tissues and have transstadial and transovarial transmission [5,6]. This adaptation ensures the long-lasting persistence of infection in endemic areas, since ticks remain infected and infective for several generations without the need of reinfection [6]. In the vertebrate hosts, *Babesia* spp. infect the erythrocytes, producing a disease characterized for hemolytic anemia. Bovines (cattle) are the maintenance hosts of *B. bovis* and *B. bigemina*. Other ruminants including *Odocoileus virginianus* (white-tailed deer) [7], *Syncerus caffer* (African buffalo) [8], and *Bubalus bubalis* (water buffalo) [9] can be carriers of these pathogens. However, whether carrier animals can act as a reservoir of these pathogens remains to be established [9].

*Rhipicephalus microplus* is the most important tick infesting cattle and represents 57% of the ixodid ticks' population in Cuba [10,11]. The climatic conditions that prevail in Cuba allow the development of tick life cycle on grassland during all months of the year [12,13]. Buffaloes are also parasitized by *R. microplus* in Cuba, although with a low infestation load [14]. Benitez et al. [15] showed experimentally that *R. microplus* can complete its life cycle in buffaloes. These authors also observed that only 5% of the tick larvae that fed on buffaloes completed their life cycle, while in cattle, 12% of the larvae completed the cycle. The elimination of ticks at different life stages can influence the dynamics of host infection by these piroplasms since only larvae transmit *B. bovis* while *B. bigemina* is transmitted by nymphs and adults [5].

Buffalo herds are increasing in Cuba and often coexist with cattle herds [16]. These animals may be carriers of *Babesia* spp. [9,17,18]. Therefore, they are potential reservoirs for these pathogens and can affect the eco-epidemiology of babesiosis in cattle [19–22]. However, evidence of the vectorial capacity of *R. microplus* that feed on carrier buffalo is scarce. As far as we know, no study has measured the transovarial transmission of *Babesia* spp. on ticks that feed on buffaloes. Thus, in this study, we verified the levels of *R. microplus* infestation in calves and adult buffaloes and compared to cattle, when raised together or in distant farms. We also quantified the levels of infection by *B. bovis* and *B. bigemina* in the progeny of ticks collected on these animals.

# 2. Results and Discussion

In this work, no tick infestation was found in the adult buffaloes of either the neighboring or control herds (Table 1). In contrast, buffalo calves were infested in both groups, with prevalence similar to that of adult cattle and cattle calves ($F = 1.52, p > 0.05$). The highest parasitic load among all the groups was observed in buffalo calves in immediate proximity to cattle herds ($F = 2.49, p = 0.04$) (Table 1). Similar results were obtained in Thailand, where Nithikathkul et al. [23] observed that adult buffalos were less infested by *R. microplus* than cattle in endemic areas. In addition, a previous study in the same region in Cuba reported that, unlike adult buffalo, buffalo calves were susceptible to tick infestation, presumably because they had thinner and hairier skin [14]. This result suggests a natural resistance of adult buffaloes to ticks. The resistance of adult buffaloes to tick infestation seems to be acquired since calves are not resistant to this ectoparasite. Another factor to consider may be that buffalo calves usually have less access to wallows (water and mud) than adult buffalo in the study region, which constitutes a biological control for ticks [15,19].

The nymphs/larvae ratio was higher in bovine calves when compared to buffalo calves ($F = 40.9, p < 0.05$, Table 1). These results are consistent with a study conducted in field conditions in Brazil in which the proportion of nymphs (21.3%) observed in buffaloes was much lower than that of larvae (74.7%). In that work, they also observed that larvae attack in buffalo calves happened mainly during the firsts months of life, usually causing intense inflammatory reactions at the skin sites more bitten
by ticks [24,25]. Furthermore, our results are comparable with the findings of Benitez et al. [15], who carried out an experimental infestation of cattle and buffalo calves in Argentina using 10,000 *R. microplus* larvae and observed that the larval survival was lower in buffalo than in cattle [15]. These authors also observed that only 5.4 % and 12% of larvae survived until the adult stages in buffalo and cattle, respectively. Benitez et al. [15] also described an intense inflammatory response on buffalo at skin sites where a large number of larvae were attached. The inflammatory reaction was not observed in cattle [15], suggesting that the immune system of buffalo may be more reactive to tick saliva components than that of cattle. Accordingly, in this work, inflammatory reactions were observed in the skin of buffaloes infested by ticks (Figure S2), but were not observed in cattle. These results suggest that buffalo acquire immune resistance to tick infestation, whereas cattle remain tolerant or susceptible to ticks.

<table><caption>**Table 1.** *R. microplus* infestation on buffalo and cattle from different ages and closeness groups.</caption><thead><tr><th rowspan="2">Indicators *</th><th rowspan="2">Age Group</th><th colspan="2">Neighboring-herds</th><th colspan="2">Control-herds</th></tr><tr><th>Cattle</th><th>Buffalo</th><th>Cattle</th><th>Buffalo</th></tr></thead><tbody><tr><td rowspan="2">(TIP) Tick infestation prevalence</td><td>Adults</td><td>80%</td><td>0%</td><td>60%</td><td>0%</td></tr><tr><td>Calves</td><td>50%</td><td>100%</td><td>70%</td><td>60%</td></tr><tr><td rowspan="2">(TIR) Tick infestation rate (x̄ ± SD)</td><td>Adults</td><td>7 ± 3<sup>ab</sup></td><td>0</td><td>5 ± 2<sup>b</sup></td><td>0</td></tr><tr><td>Calves</td><td>7 ± 4<sup>ab</sup></td><td>12 ± 10<sup>a</sup></td><td>5 ± 3<sup>b</sup></td><td>5 ± 4<sup>b</sup></td></tr><tr><td rowspan="2">Nymphs/larvae ratio (x̄ ± SE)</td><td>Adults</td><td>-</td><td>-</td><td>-</td><td>-</td></tr><tr><td>Calves</td><td>36 ± 0.009%<sup>c</sup></td><td>30 ± 0.009%<sup>b</sup></td><td>35 ± 0.009%<sup>c</sup></td><td>24 ± 0.008%<sup>a</sup></td></tr></tbody></table>

The immunological mechanisms underlying the natural resistance of buffalo to tick infestations remain to be identified. In cattle, host genetic factors define the immune resistance to ticks [26], and tick resistance is an inherited trait in *Bos taurus indicus* that can be passed to *B. t. taurus* by crossbreeding of these two subspecies [27]. The precise mechanisms underlying tick resistance in cattle remains poorly understood, but several host genetic factors associated with tick resistance and susceptibility have been identified [26,28]. Despite a combination of humoral and cellular immune response mechanisms seeming to be involved [26], the evidence for a correlation between the level of tick infestation and the level of immunoglobulin G (IgG) to tick antigens is contradictory [26]. However, there is consensus in that different subsets of inflammatory cells play an important role in tick resistance in bovines [26,29,30]. For example, higher numbers of immune cells involved in the histamine response of the skin, eosinophils, mast cells, and basophils, were found in the skin of resistant cattle [26,28].

The weight of fully engorged adult females that fed on buffalo calves (0.150 ± 0.028 g) was similar to that of ticks that fed on cattle (0.151 ± 0.033 g). Similar results were reported by Benitez et al. [15], who found no differences in the body weight of adult female ticks that fed on cattle and buffaloes, or in the times of preoviposition, oviposition, and incubation. The reproductive performance of adult female ticks is shown in Table 2. No significant difference was found in the reproductive parameters CE (F = 0.93, p = 0.47), HR (F = 0.34, p = 0.88), and REI (F = 0.77; p = 0.57) of adult females fed on bovine or buffalo calves. After observing similar results, Benitez et al. [15] concluded that *R. microplus* adult females can complete their life cycles in buffalo calves and cattle with comparable reproductive efficiency.

The qPCR analysis revealed that the 31 groups of tick larvae (Figure S3) and the 31 hosts tested were positive for *B. bovis*. However, only 68% (21/31) of the larval groups was positive for *B. bigemina*. Of the 10 larval groups negative for *B. bigemina*, three were from qPCR-negative buffalo calves, and the other seven were from adult females fed on qPCR-positive animals (four buffalo calves and three adult cattle). The infection levels of *B. bovis* were similar in larvae from females fed on buffalo and bovine calves (Table 3), while larvae from adult cattle harbor the highest values of infection (F = 4.90, p = 0.05). In the case of *B. bigemina* there were no differences in the infection levels in tick larvae relative to host type (F = 3.06, p = 0.06). Besides, we found higher values in the count of parasitic infection by *B. bovis* than *B. bigemina* (t = 7.08, p < 0.05).

<footnote>* Different letters in the same indicator indicate values significantly different between groups (p ≤ 0.05).</footnote>
<table><caption>**Table 2.** Reproductive performance of *R. microplus* engorged females fed on buffaloes and cattle from different age groups and distances between herds.</caption><thead><tr><th rowspan="2">Indicators</th><th rowspan="2">Age Group</th><th colspan="2">Neighboring-herds</th><th colspan="2">Control-herds</th></tr><tr><th>Cattle</th><th>Buffalo</th><th>Cattle</th><th>Buffalo</th></tr></thead><tbody><tr><td rowspan="2">(CE) Conversion efficiency (x̄ ± SE)</td><td>Adults</td><td>43 ± 0.04%</td><td>-</td><td>42 ± 0.07%</td><td>-</td></tr><tr><td>Calves</td><td>38 ± 0.05%</td><td>49 ± 0.03%</td><td>41 ± 0.05%</td><td>50 ± 0.07%</td></tr><tr><td rowspan="2">(IR) Incubation rate (x̄ ± SE)</td><td>Adults</td><td>98 ± 0.02%</td><td>-</td><td>99 ± 0.04%</td><td>-</td></tr><tr><td>Calves</td><td>98 ± 0.02%</td><td>95 ± 0.01%</td><td>96 ± 0.03%</td><td>99 ± 0.04%</td></tr><tr><td rowspan="2">(REI) Reprod. efficiency index (x̄ ± SD)</td><td>Adults</td><td>8.4 ± 2.0 × 10<sup>5</sup></td><td>-</td><td>8.4 ± 5.2 × 10<sup>5</sup></td><td>-</td></tr><tr><td>Calves</td><td>7.4 ± 3.0 × 10<sup>5</sup></td><td>9.5 ± 2.2 × 10<sup>5</sup></td><td>7.9 ± 1.6 × 10<sup>5</sup></td><td>9.8 ± 2.0 × 10<sup>5</sup></td></tr></tbody></table>

<table><caption>**Table 3.** Molecular quantification (qPCR) of *B. bovis* and *B. bigemina* infection in the offspring of adult females of *R. microplus* fed on buffalo and cattle.</caption><thead><tr><th rowspan="2">Hemoparasite *</th><th rowspan="2">Age Group</th><th colspan="2">Neighboring-herds</th><th colspan="2">Control-herds</th></tr><tr><th>Cattle</th><th>Buffalo</th><th>Cattle</th><th>Buffalo</th></tr></thead><tbody><tr><td rowspan="2">B. bovis (LogCN/µL)</td><td>Adults</td><td>8.5 ± 1.3<sup>b</sup></td><td>-</td><td>8.6 ± 1.4<sup>b</sup></td><td>-</td></tr><tr><td>Calves</td><td>7.8 ± 0.7<sup>b</sup></td><td>7.3 ± 0.6<sup>a</sup></td><td>7.7 ± 1.0<sup>b</sup></td><td>7.2 ± 0.5<sup>a</sup></td></tr><tr><td rowspan="2">B. bigemina (LogCN/µL)</td><td>Adults</td><td>5.3 ± 1.0<sup>ab</sup></td><td>-</td><td>4.8 ± 1.4<sup>a</sup></td><td>-</td></tr><tr><td>Calves</td><td>7.2 ± 1.3<sup>b</sup></td><td>5.6 ± 0.7<sup>a</sup></td><td>6.0 ± 2.0<sup>ab</sup></td><td>5.4 ± 0.4<sup>a</sup></td></tr></tbody></table>

The higher frequency of *B. bovis* compared to *B. bigemina* in the offspring of adult females of *R. microplus* fed on cattle is a distinctive epidemiological feature of our study, since most studies report that *B. bigemina* infection is more frequent than *B. bovis* in ticks fed on cattle [5,31]. Specifically, research conducted in Australia revealed that only ~0.04% of larval ticks were infected with *B. bovis* on farms with *B. t. taurus* cattle, whereas the infection with *B. bigemina* was usually higher (i.e., ~0.23%) [32]. As a consequence, a higher transmission rate of *B. bigemina* than of *B. bovis* could be expected in regions where both *Babesia* species are present, therefore, causing the highest prevalence of *B. bigemina* in cattle herds [5,33]. However, our results are consistent with those of an epidemiological study previously conducted in the same region as our study, in which a high molecular prevalence of both *B. bovis* and *B. bigemina* was found in buffalo [22]. Eco-epidemiological factors favoring the presence of *B. bovis* in this region of Cuba should be studied, as *B. bovis* could cause severe disease in susceptible cattle, particularly in South America and Australia [5,34].

A limitation in this work was that the qPCR quantification of *Babesia* spp. was carried out on tick larval pools (~100 larvae from egg batches laid by 10 engorged females). This approach could have influenced the high frequency of positive samples detected in our study, thus overestimating the level of transovarial transmission of *Babesia* spp., despite low transovarial transmissibility of these apicomplexan parasites [35]. Future studies including water buffaloes should quantify the infection level in each larva and establish the proportion of larvae that acquired the *Babesia* transovarially. This aspect is essential to determine the likelihood of host exposure to these hemoparasites when buffaloes are present in endemic areas with a potential impact on the epidemiology of bovine babesiosis (e.g., water buffalo could have a dilution effect [36] on the dynamics of *Babesia* spp. transmission). In this regard, studies carried out in Australia compared the transmission rates of *B. bovis* in *B. t. taurus* and *B. t. indicus* and concluded that the breeding of tick-resistant cattle might lead to the disappearance of ticks in environments not suitable for tick survival [33]. This situation could lead to low infection rates and lack of immunity on cattle herds, therefore, posing a threat to endemic stability [5,33].

In this work, the molecular quantification of the pathogens in vertebrate hosts revealed that the infection levels of *B. bovis* were higher in cattle than in buffalo when both species coexisted (neighboring-herds) (F = 7.98, p < 0.05), however, the infection levels were similar to those observed

<footnote>*Different letters in the same hemoparasite indicate values significantly different between groups (p ≤ 0.05).</footnote>
in the control groups. No differences were observed between age groups from the same host species (Table 4). In the case of *B. bigemina*, the highest infection levels were found in bovine calves ($F = 17.38, p < 0.05$), while the values were similar among buffalo calves, adult buffalo, and adult cattle. No correlation was found between the parasitic loads of *B. bovis* (CN/µL of DNA) in tick larvae and mammalian hosts ($r = 0.2, p = 0.27$). In contrast, there was a correlation between the parasite loads in tick and vertebrate hosts ($r = 0.4, p = 0.03$) of *B. bigemina*. There are sparse references on this subject, although Giglioti et al. [35] found a high correlation ($p < 0.01$) on the infection levels of *B. bovis* ($r = 0.58$) and *B. bigemina* ($r = 0.6$) between carrier cattle and ticks. These authors suggest the existence of a cause-effect relationship between parasitic loads in cattle and ticks; additional studies should address this issue.

<table><caption>**Table 4.** Molecular quantification (qPCR) of *B. bovis* and *B. bigemina* infection on carrier cattle and buffaloes.</caption><thead><tr><th rowspan="2">Hemoparasite *</th><th rowspan="2">Age Group</th><th colspan="2">Neighboring-herds</th><th colspan="2">Control-herds</th></tr><tr><th>Cattle</th><th>Buffalo</th><th>Cattle</th><th>Buffalo</th></tr></thead><tbody><tr><td rowspan="2">B. bovis (LogCN/µL)</td><td>Adults</td><td>8.1 ± 0.3<sup>c</sup></td><td>6.6 ± 1.4<sup>ab</sup></td><td>8.0 ± 0.0<sup>bc</sup></td><td>7.5 ± 1.0<sup>bc</sup></td></tr><tr><td>Calves</td><td>8.3 ± 0.3<sup>c</sup></td><td>6.7 ± 0.9<sup>a</sup></td><td>8.1 ± 0.3<sup>c</sup></td><td>7.5 ± 0.1<sup>abc</sup></td></tr><tr><td rowspan="2">B. bigemina (LogCN/µL)</td><td>Adults</td><td>6.1 ± 0.8<sup>b</sup></td><td>5.0 ± 1.0<sup>ab</sup></td><td>6.3 ± 1.1<sup>abc</sup></td><td>4.8 ± 1.8<sup>ab</sup></td></tr><tr><td>Calves</td><td>9.0 ± 0.7<sup>d</sup></td><td>6.5 ± 0.7<sup>ab</sup></td><td>8.4 ± 0.7<sup>cd</sup></td><td>5.0 ± 0.1<sup>a</sup></td></tr></tbody></table>

In general, our results showed that adult females of *R. microplus* fed on water buffalo can produce offspring infected by *B. bovis* and *B. bigemina*. This fact contributes as a piece of evidence of the participation of water buffalo as a potential reservoir of *B. bovis* and *B. bigemina*. The potential impact of these results for the epidemiology of babesiosis in western Cuba, and perhaps in other endemic regions where cattle and buffalo are bred together, remains to be fully elucidated. Further investigation is needed to assess: (1) Whether uninfected ticks can acquire *Babesia* when feeding on *Babesia*-infected buffaloes and (2) whether transovarial transmission of *Babesia* is equally efficient in uninfected ticks that acquire the pathogen while feeding on the buffalo as in ticks carrying the infection from previous generations. Another important question that remains to be answered is the relative contribution of horizontal transmission and vertical transmission in the epidemiological cycles of the system *R. microplus*-buffalo-*Babesia*.

# 3. Materials and Methods

## 3.1. Study Site and Samples' Selection

A cohort study was conducted in livestock of the municipality of San Jose, Mayabeque, Cuba. The climate of the region is tropical, seasonally humid, with an annual average of temperature between 22 and 28 °C, and relative humidity of 80% [37]. *R. microplus* is the only tick species identified in cattle and buffaloes of the region [38]. This tick species persists throughout the year with peaks of high-abundance population during the dry season (December–March) [38]. *B. bovis*, *B. bigemina*, and *A. marginale* are endemic in cattle populations of this region [4]. The prevalence of these pathogens in buffalo of the region is high [22,39]. For this study, two farms were selected where cattle and buffalo herds were raised together, either in immediate proximity separated by wire fences (hereafter the "Neighboring herd" group) or separated by a distance of 10 km (hereafter the "Control herd" group). Each group, neighboring and control, consisted of 20 cattle and 20 buffalo randomly selected. Each cohort consisted of 10 calves (0–1 year) and 10 adults (>5 years). The experimental design is shown in Figure 1.

<footnote>*Different letters in the same hemoparasite indicate values significantly different between groups ($p \le 0.05$).</footnote>
<figure><img src="image_2.png" /><figcaption>**Figure 1.** Diagram of experimental design. A cohort study was carried out, with two study groups, consisting of 20 buffaloes and 20 cattle each. In one group, bovines and buffalo cohabited on adjoining ranches, with grazing areas divided by a perimeter fence (Neighboring herds), while the other group was the reference panel (Control herds), consisting of two groups of cattle and buffalo from distant herds (~10 km between them). Besides, each subgroup consisted of 10 calves (0–1 year) and 10 adult animals (>5 years).</figcaption></figure>

## 3.2. Tick Infestation

Adult female ticks (≥4.5 mm) [40] were counted on the left side of each animal, and the resulting value was multiplied by two to determine the parasite burden. This data were used to calculate the prevalence of tick infestation (TIP) (i.e., % of infested animals) and the intensity of tick infestation (TIR) (i.e., adult females/infested animals) in each group [41]. Tick counting was done between 7:00 to 8:00 a.m. [42]. The nymph/larvae ratio was compared between buffalo and bovine calves to determine the larval survival rate. For this purpose, the tick life stages were counted as described by Veríssimo et al. [43]; briefly, a 10 cm diameter area of the perineal region of each animal was scraped with a razor blade (Figure S1). The scraping material was collected and placed in tubes with a screw cap (15 mL) containing 3 mL of distilled water. The material was transferred to Petri dishes, and up to 100 ticks were counted on each plate using a Stemi DV4 Stereo microscope (Carl Zeiss, Oberkochen, Germany). In the counts, the numbers of larvae or nymphs were discriminated, and adult ticks were not considered. The nymph/larvae ratio was calculated as: Nymphs (%) = (nymphs/nymphs + larvae) × 100.

## 3.3. Analysis of the Reproductive Parameters of Adult Female Ticks

The reproductive performance was evaluated in adult female ticks collected from 31 hosts (12 buffalo calves, nine bovine calves, and 10 adult cattle). Ticks were manually detached at the end of engorgement (≥4.5 mm, with dorsal grooves not visible, [42]) and placed in 15-mL conical tubes, covered with cotton. In the laboratory, groups of up to 10 adult females were placed in Petri dishes and incubated at 27 ± 1 °C and 85-86% relative humidity for 15 days for oviposition to occur, as described by Oliveira et al. [44]. The egg mass of each group of adult females was weighed, placed in closed syringes with cotton, and incubated for 21 days. The larval hatching rate (HR) was determined as previously described by Giglioti et al. [45]. Three estimators of the reproductive performance of female ticks were calculated: (1) Conversion efficiency (CE): $CE\% = (WE / WT) \times 100$; (2) hatching rate (HR): $HR\% = (NL/NL + NE) \times 100$; and (3) reproductive efficiency index (REI): $REI = (WE \times IR \times 20,000)/WT$, where: WE is the weight of eggs, WT is the weight of female, NL is the number of larvae, NE is the number of eggs not incubated, and 20,000 is an estimated number of larvae that are produced from 1 g of eggs.
### 3.4. DNA Extraction and Quantification of Babesia Infection Levels

A blood sample was collected from each animal in 4 mL vacutainer EDTA tubes (BD Vacutainer Blood Collection, Franklin Lakes, NJ, USA), and kept frozen at -80 °C in cryovials until processing. The genomic DNA was extracted from 300 µL of blood by using the Wizard® Genomic DNA Purification Kit (Promega, Madison, WI, USA) according to the manufacturer's instructions. For DNA extraction from ticks, approximately 100 larvae were taken from each sample, poured into a sterile 1.5 mL Eppendorf tube, flash-frozen in liquid nitrogen and immediately crushed with a plastic macerator, and homogenized in a Tissue Lyser® (Qiagen, Hilden, Germany). DNA extraction was performed using DNeasy® Blood and Tissue DNA Purification Kit (Qiagen, Hilden, Germany). DNA quantity and quality were evaluated using a Nanodrop Spectrophotometer 1000 v.3.5 (Thermo Fisher Scientific, Waltham, MA, USA) and stored at -20 °C. The SYBR green-based qPCR was used for high sensitivity detection of *B. bovis* and *B. bigemina*. The primers used were those described by Buling et al. [46], which amplify 88 base pairs (bp) fragments of the mitochondrial cytochrome b gene (*mt-Cytb*). The standardized thermocycling conditions, as well as the quantification procedures, were as described previously by Obregon et al. [47].

### 3.5. Statistical Analysis

The TIP was compared between groups with the Fisher Exact test and Duncan test (95% confidence interval (CI)), while the ratio nymphs/larvae and the TIR were analyzed using the ANOVA test and Duncan's multiple range test (95% CI). The ANOVA and Duncan's multiple range tests were also used to compare the reproductive efficiency indicators (CE, HR, and REI) between groups, as well as the infection levels of hemoparasites (measured as DNA copy number (CN), CN/µL) in the DNA samples from tick larvae and mammalian hosts. Pearson's correlation coefficient (r) was used to explore the association on the infection levels of *B. bovis* and *B. bigemina* between tick larvae and the hosts wherein the adult females fed. The analyses were performed using the statistical software package Statgraphics Centurion v. 16.1.03 (StatPoint technologies Inc, Warrenton, VA, USA).

### 3.6. Ethics Statement

All procedures in this work were according to the principles established by the International Guiding Principles for Biomedical Research Involving Animals (2012). The committee on ethics and animal welfare at Centro Nacional de Sanidad Agropecuaria (CENSA) approved this research. No animal was sacrificed for the purposes of the study and the field study did not involve endangered or protected species. The epidemiological survey did not involve harm or cruelty to animals. The border veterinary service of the Institute of Veterinary Medicine (IMV) from Cuba authorized the exportation of DNA samples used in this work (certification number R.S.1522010).

## 4. Conclusions

With this investigation, we concluded that *R. microplus* infestation, although rare in adult buffaloes, is similar in frequency and intensity in cattle and buffalo calves living in endemic areas in western Cuba. We found that the reproductive performance of female ticks that fed on buffalo calves was similar to that of female ticks that fed on cattle. More importantly, the offspring of adult females fed on buffalo calves carried *B. bovis* and *B. bigemina*, suggesting transovarial transmission of *Babesia* spp. in ticks fed on buffalo calves. Furthermore, the levels of *Babesia* DNA were similar in the offspring of adult females that fed on buffalo and cattle. Future studies are needed to establish the transmission rate of *Babesia* spp. in regions where the life cycle of *R. microplus* is completed in buffaloes. Other components of the possible reservoir competence of water buffalo for *B. bovis* and *B. bigemina*, such as the infectiousness of the strains transmitted by ticks fed on buffalos, should also be analyzed.
**Supplementary Materials:** The following are available online at www.mdpi.com/2076-0817/9/4/280/s1, Figure S1: Tick sampling procedure for tick nymph/larval ratio analysis. The figures represent the main four stages of sampling: (1) First, a circular area (10 cm in diameter) was delimited in the perineal region of each bovine and buffalo calf, (2) then the skin was moistened with soapy water, (3) and was shaved with a razor and sterile blade, and (4) the scraping was placed in capped tubes (15 mL) containing 3 mL of distilled water. Figure S2: The skin of buffalo calves in anatomic areas infested with *R. microplus* ticks. The arrow's point indicates lesions due to inflammatory reactions in response to bites and tick fixation. Figure S3: Real-time qPCR analysis of tick larvae for babesial hemoparasites. (a) Amplification signal from tick DNA samples on the qPCR assay for *B. bovis*. (b) Melting curve analysis of resulting *B. bovis* amplicons (melting peaks at 79.5 °C. (c) Amplification signal from tick DNA samples on the qPCR assay for *B. bigemina*. (d) Melting curve analysis of resulting *B. bigemina* amplicons (melting peaks at 78.5 °C. Each sample was analyzed in duplicate, NTC, no template control. Graphics automatically generated by the software Bio-Rad CFX96 Manager v. 3.1 (BioRad, USA).

**Author Contributions:** Conceptualization, D.O., B.C.-G., and E.R.; formal analysis, D.O., A.A.D.-S., and Y.A.; funding acquisition, B.C.-G.; investigation, D.O., A.A.D.-S., and Y.A.; methodology, D.O., A.A.D.-S., and M.C.S.O.; resources, B.C.-G., M.C.S.O., and A.C.-C.; supervision, E.R.; writing—original draft, D.O. and A.C.-C.; writing—review and editing, D.O., B.C.-G., A.A.D.-S., E.R., Y.A., and M.C.S.O. All authors have read and agreed to the published version of the manuscript.

**Funding:** This work was supported by funding from the Brazilian National Council for Scientific and Technological Development (CNPq-474648/210-9), Brazil; Empresa Brasileira de Pesquisa Agropecuária (Embrapa Pecuária Sudeste); and the National Priority Program for Animal and Plant Health, Ministry of Agriculture (MINAGRIP131LH003007).

**Acknowledgments:** Cuba. Dasiel Obregon would like to acknowledge the CNPq and the Three World Academy of Science (TWAS) for his sandwich Ph.D. grant (TWAS-CNPq 2013, Process 302999/2013-2). We are grateful to Rosangela Z. Machado (Universidade Estadual Paulista-Campus de Jaboticabal) for her recommendations for this work.

**Conflicts of Interest:** The authors declare that they have no commercial or financial relationships that could be construed as a potential conflict of interest.

**Data Availability Statement:** All data generated or analyzed during this study are included in this published article.

**References**

1. Suarez, C.E.; Noh, S. Emerging perspectives in the research of bovine babesiosis and anaplasmosis. *Vet. Parasitol.* **2011**, *180*, 109–125.

2. OIE. Babesiosis bovina. In *Manual of Diagnostic Tests and Vaccines for Terrestrial Animals* 2012; World Organization for Animal Health: Paris, France, 2012; pp. 1–16.

3. Nari, A. Strategies for the control of one-host ticks and relationship with tick-borne diseases in South America. *Vet. Parasitol.* **1995**, *57*, 153–165.

4. Alonso, M.; Arellano-Sota, C.; Cereser, V.H.; Cordoves, C.O.; Kessler, A.A.; Solari, M.A.; Vega, C.A. Epidemiology of bovine anaplasmosis and babesiosis in Latin America and the Caribbean. *Int. J. Epidemiol.* **1992**, *11*, 713–733.

5. Bock, R.; Jackson, L.; Vos, A.; Jorgensen, W. Babesiosis of cattle. *Parasitology* **2004**, *129*, S247–S269.

6. Chauvin, A.C.; Oreau, E.M.; Onnet, S.B.; Lantard, O.P.; Alandrin, L.M.; 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.

7. Cantu, A.; Ortega-S, J.A.; Mosqueda, J.; Garcia-Vazquez, Z.; Henke, S.E.; George, J.E. Immunologic and molecular identification of *Babesia bovis* and *Babesia bigemina* in free-ranging white-tailed deer in northern Mexico. *J. Wildl. Dis.* **2007**, *43*, 504–507.

8. Iseki, H.; Alhassan, A.; Ohta, N.; Thekisoe, O.M.M.; Yokoyama, N.; Inoue, N.; Nambota, A.; Yasuda, J.; Igarashi, I. Development of a multiplex loop-mediated isothermal amplification (mLAMP) method for the simultaneous detection of bovine *Babesia* parasites. *J. Microbiol. Methods* **2007**, *71*, 281–287.

9. Obregón, D.; Oliveira, M.C.S.; Tizioto, P.; Funnes, E.M.; Martínez, S.; Roque, E.; Fonseca, A.H.; Corona, B. Diagnóstico de *Babesia bovis* en búfalos de la región occidental de Cuba a través de un ensayo de nPCR. *Rev. Salud Anim.* **2012**, *34*, 101–108.
10. Barros-Battesti, D.; Reyes, M.H.; Onofrio, V.C.; Beati, L.; Famada, K.; Guglielmone, A.A. The ixodid ticks (Acari: Ixodidae) of Cuba. *Syst. Appl. Acarol.* **2009**, *14*, 101–128.

11. Gaínza, Y.A.; Martínez-Marrero, S.; Corona-González, B. Ticks of veterinary interest in Cuba, and its importance in the changing climatic conditions. *Rev. Electrónica Vet.* **2014**, *15*, 1–22.

12. Estrada-Peña, A.; Bouattour, A.; Camicas, J.-L.; Guglielmone, A.; Horak, I.; Jongejan, F.; Latif, A.; Pegram, R.; Walker, A. The known distribution and ecological preferences of the tick subgenus *Boophilus* (Acari: Ixodidae) in Africa and Latin America. *Exp. Appl. Acarol.* **2006**, *38*, 219–235.

13. Estrada-Peña, A.; Naranjo, V.; Acevedo-Whitehouse, K.; Mangold, A.J.; Kocan, K.M.; De la Fuente, J. Phylogeographic analysis reveals association of tick-borne pathogen, *Anaplasma marginale*, MSP1a sequences with ecological traits affecting tick vector performance. *BMC Biol.* **2009**, *7*, 1–13.

14. Obregón, D.; Rodríguez, J.D.; Roque, E.; Alemán, Y. *Rhipicephalus (Boophilus) microplus* (acari: Ixodidae) en búfalos (*Bubalus bubalis*), en Cuba. *Rev. Salud Anim.* **2010**, *32*, 132–134.

15. Benitez, D.; Cetrá, B.; Florin-Christensen, M. *Rhipicephalus (Boophilus) microplus* ticks can complete their life cycle on the water buffalo (*Bubalus bubalis*). *J. Buffalo Sci.* **2012**, *1*, 193–197.

16. Mitat, A.O.B. Búfalos de agua en Cuba: Origen y evolución. *Rev. ACPA* **2009**, *3*, 45-48.

17. Ferreri, L.; Benitez, D.; Dominguez, M.; Rodriguez, A.; Asenzo, G.; Mesplet, M.; Florin-Christensen, M.; Schnittger, L. Water Buffalos as carriers of *Babesia bovis* in Argentina. *Ann. N. Y. Acad. Sci.* **2008**, *1149*, 149–151.

18. Néo, T.A.; Giglioti, R.; Obregón, D.; Bilhassi, T.B.; Oliveira, H.N.; Machado, R.Z.; Aníbal, F.F.; Brito, L.G.; Malagó, W.; Bressani, F.A.; et al. Detection of *Babesia bovis* and *B. bigemina* in Water buffaloes (*Bubalus bubalis*) in endemic areas of São Paulo State, Brazil. *Open J. Vet. Med.* **2016**, *6*, 75–84.

19. Terkawi, M.A.; Huyen, N.X.; Shinuo, C.; Inpankaew, T.; Maklon, K.; Aboulaila, M.; Ueno, A.; Goo, Y.K.; Yokoyama, N.; Jittapalapong, S.; et al. Molecular and serological prevalence of *Babesia bovis* and *Babesia bigemina* in water buffaloes in the northeast region of Thailand. *Vet. Parasitol.* **2011**, *178*, 201–207.

20. Romero-Salas, D.; Mira, A.; Mosqueda, J.; García-Vázquez, Z.; Hidalgo-Ruiz, M.; Vela, N.A.O.; de León, A.A.P.; Florin-Christensen, M.; Schnittger, L. Molecular and serological detection of *Babesia bovis* and *Babesia bigemina*-infection in bovines and water buffaloes raised jointly in an endemic field. *Vet. Parasitol.* **2016**, *217*, 101–107.

21. Silva, J.B.; André, M.R.; da Fonseca, A.H.; de Albuquerque Lopes, C.T.; da Silva Lima, D.H.; de Andrade, S.J.T.; Oliveira, C.M.C.; Barbosa, J.D. Molecular and serological prevalence of *Babesia bovis* and *Babesia bigemina* in water buffaloes in the north region of Brazil. *Vet. Parasitol.* **2013**, *197*, 678–681.

22. Obregón, D.; Cabezas-cruz, A.; Armas, Y.; Silva, J.B.; Fonseca, A.H.; André, M.R.; Alfonso, P.; Oliveira, M.C.S.; Machado, R.Z.; Corona-González, B. High co-infection rates of *Babesia bovis*, *Babesia bigemina*, and *Anaplasma marginale* in water buffalo in Western Cuba. *Parasitol. Res.* **2019**, *118*, 1–13.

23. Nithikathkul, C.; Polseela, P.; Changsap, B.; Leemingsawat, S. Ixodid ticks on domestic animals in Samut Prakan Province, Thailand. *Southeast Asian J. Trop. Med. Public Health* **2002**, *33*, 41–44.

24. Gomes, R.A.; Machado, R.Z.; Starke-buzetti, W.A.; Bonesso, M.A. Immune-humoral response of water buffalo (*Bubalus bubalis*) against *Anaplasma marginale*. *Rev. Bras. Parasitol. Vet.* **2008**, *17*, 73–80.

25. Gomes, R.A. *Resposta Imune-Humoral de Búfalos (Bubalus Bubalis) Infectados Naturalmente por Babesia Bovis, B. Bigemina e Anaplasma Marginale*; UNESP-Campus de Jaboticabal: Jaboticabal, Brazil, 2007.

26. Tabor, A.E.; Ali, A.; Rehman, G.; Garcia, G.R.; Zangirolamo, A.F.; Malardo, T.; Jonsson, N.N. Cattle Tick *Rhipicephalus microplus*-host interface: A review of resistant and susceptible host responses. *Front. Cell. Infect. Microbiol.* **2017**, *7*, 1–18.

27. Porto Neto, L.R.; Jonsson, N.N.; D'Occhio, M.J.; Barendse, W. Molecular genetic approaches for identifying the basis of variation in resistance to tick infestation in cattle. *Vet. Parasitol.* **2011**, *180*, 165–172.

28. Robbertse, L.; Richards, S.A.; Maritz-Olivier, C. Bovine immune factors underlying tick resistance: Integration and future directions. *Front. Cell. Infect. Microbiol.* **2017**, *7*, 1–16.

29. Constantinoiu, C.C.; Jackson, L. a; Jorgensen, W.K.; Lew-Tabor, a E.; Piper, E.K.; Mayer, D.G.; Venus, B.; Jonsson, N.N. Local immune response against larvae of *Rhipicephalus (Boophilus) microplus* in *Bos taurus indicus* and *Bos taurus taurus* cattle. *Int. J. Parasitol.* **2010**, *40*, 865–875.

30. Piper, E.K.; Jackson, L.A.; Bielefeldt-Ohmann, H.; Gondro, C.; Lew-Tabor, A.E.; Jonsson, N.N. Tick-susceptible *Bos taurus* cattle display an increased cellular response at the site of larval *Rhipicephalus (Boophilus) microplus* attachment, compared with tick-resistant *Bos indicus* cattle. *Int. J. Parasitol.* **2010**, *40*, 431–441.
31. Quintão-Silva, M.G.; Melo, M.N.; Ribeiro, M.F.B. Comparison of duplex PCR and microscopic techniques for the identification of *Babesia bigemina* and *Babesia bovis* in engorged female ticks of *Boophilus microplus*. *Zoonoses Public Health* **2007**, *54*, 147–151.

32. Mahoney, D.F.; Mirre, G.B. Bovine babesiasis: Estimation of infection rates in the tick vector *Boophilus microplus* (Canestrini). *Ann. Trop. Med. Parasitol.* **1971**, *65*, 309–317.

33. Jonsson, N.N.; Bock, R.E.; Jorgensen, W.K.; Morton, J.M.; Stear, M.J. Is endemic stability of tick-borne disease in cattle a useful concept? *Trends Parasitol.* **2012**, *28*, 85–89.

34. Gohil, S.; Kats, L.M.; Sturm, A.; Cooke, B.M. Recent insights into alteration of red blood cells by *Babesia bovis*: Moovin' forward. *Trends Parasitol.* **2010**, *26*, 591–599.

35. Giglioti, R.; de Oliveira, H.N.; Okino, C.H.; de Sena Oliveira, M.C. qPCR estimates of *Babesia bovis* and *Babesia bigemina* infection levels in beef cattle and *Rhipicephalus microplus* larvae. *Exp. Appl. Acarol.* **2018**, *75*, 235–240.

36. Huang, Z.Y.X.; Langevelde, F.; Estrada-Peña, A.; Suzán, G.; de Boer, F. The diversity-disease relationship: Evidence and criticisms of the dilution effect. *Parasitology* **2016**, *1*, 1–12.

37. INSMET. Instituto de Meteorología de la República de CUBA. Centro del Clima. El Clima de Cuba. Características Generales. Available online: http://www.met.inf.cu/asp/genesis.asp?TB0=PLANTILLAS&TB1=CLIMAC&TB2=/clima/ClimaCuba.htm (accessed on 1 January 2016).

38. LNP (National Parasitology Laboratory) *Serie Histórica de la Anaplasmosis y Babesiosis en Cuba*; Sistema de información estadística del IMV. Annual balance report: 2014; Ministry of Agriculture of the Republic of Cuba: Habana, Cuba, 2014.

39. Obregón, D.; Corona, B.G.; de la Fuente, J.; Cabezas-Cruz, A.; Gonçalves, L.R.; Matos, C.A.; Armas, Y.; Hinojosa, Y.; Alfonso, P.; Oliveira, M.C.S.; et al. Molecular evidence of the reservoir competence of water buffalo (*Bubalus bubalis*) for *Anaplasma marginale* in Cuba. *Vet. Parasitol. Reg. Stud. Reports* **2018**, *13*, 180–187.

40. De la Vega, R.; Camejo, A.; Dias, G. Dynamics of *Boophilus microplus* (Acari: Ixodidae) female feeding on bovine. *Rev. Salud Anim.* **2003**, *25*, 192–195.

41. De la Vega, R.; Moreno, A.; Diaz, G. Sampling method of the cattle tick (*Boophilus microplus*) in dairy cows. *Rev. Salud Anim.* **1984**, *6*, 397–406.

42. De la Vega, R.; Díaz, G.; Camejo, A.; García, I. Application of tick sampling on acaricidal bath frequency in *Rhipicephalus microplus* (acari: Ixodidae). *Rev. Salud Anim.* **2009**, *31*, 105–107.

43. Veríssimo, C.J.; da Silva, R.G.; Oliveira, A.D.; Ribeiro, W.R.; Rocha, U.F. *Boophilus micropius* tick instar counts on a crossbred herd. *B. Ind. Anim* **1997**, *54*, 19–24.

44. Oliveira, M.C.S.; Oliveira-Sequeira, T.C.G.; Araujo, J.P.; Amarante, a F. T.; Oliveira, H.N. *Babesia* spp. infection in *Boophilus microplus* engorged females and eggs in Sao Paulo State, Brazil. *Vet. Parasitol.* **2005**, *130*, 61–67.

45. Giglioti, R.; Forim, M.R.; Oliveira, H.N.; Chagas, A.C.S.; Ferrezini, J.; Brito, L.G.; Falcoski, T.O.R.S.; Albuquerque, L.G.; Oliveira, M.C.S. In vitro acaricidal activity of neem (*Azadirachta indica*) seed extracts with known azadirachtin concentrations against *Rhipicephalus microplus*. *Vet. Parasitol.* **2011**, *181*, 309–315.

46. Buling, A.; Criado-Fornelio, A.; Asenzo, G.; Benitez, D.; Barba-Carretero, J.C.; Florin-Christensen, M. A quantitative PCR assay for the detection and quantification of *Babesia bovis* and *B. bigemina*. *Vet. Parasitol.* **2007**, *147*, 16–25.

47. Obregón, D.; Rabelo, M.D.; Giglioti, R.; Bilhassi, T.B.; Néo, T.A.; Corona, B.; Alfonso, P.; Machado, R.Z.; Oliveira, M.C.S. Standardization of a SYBR green based real-time PCR system for detection and molecular quantification of *Babesia bovis* and *B. bigemina* in water buffaloes (*Bubalus bubalis*). *J. Buffalo Sci.* **2016**, *5*, 44–52.

<figure><img src="image_3.png" /></figure>