- Research article
- Open Access
Allopatric speciation in ticks: genetic and reproductive divergence between geographic strains of Rhipicephalus (Boophilus) microplus
BMC Evolutionary Biologyvolume 9, Article number: 46 (2009)
The cattle tick, Rhipicephalus (Boophilus) microplus, economically impact cattle industry in tropical and subtropical regions of the world. The morphological and genetic differences among R. microplus strains have been documented in the literature, suggesting that biogeographical and ecological separation may have resulted in boophilid ticks from America/Africa and those from Australia being different species. To test the hypothesis of the presence of different boophilid species, herein we performed a series of experiments to characterize the reproductive performance of crosses between R. microplus from Australia, Africa and America and the genetic diversity of strains from Australia, Asia, Africa and America.
The results showed that the crosses between Australian and Argentinean or Mozambican strains of boophilid ticks are infertile while crosses between Argentinean and Mozambican strains are fertile. These results showed that tick strains from Africa (Mozambique) and America (Argentina) are the same species, while ticks from Australia may actually represent a separate species. The genetic analysis of mitochondrial 12S and 16S rDNA and microsatellite loci were not conclusive when taken separately, but provided evidence that Australian tick strains were genetically different from Asian, African and American strains.
The results reported herein support the hypothesis that at least two different species share the name R. microplus. These species could be redefined as R. microplus (Canestrini, 1887) (for American and African strains) and probably the old R. australis Fuller, 1899 (for Australian strains), which needs to be redescribed. However, experiments with a larger number of tick strains from different geographic locations are needed to corroborate these results.
The cattle tick, Rhipicephalus (Boophilus) microplus, is distributed in tropical and subtropical regions of the world . Infestations with R. microplus economically impact cattle industry by reducing weight gain and milk production, and by transmitting pathogens that cause babesiosis (Babesia bovis and B. bigemina) and anaplasmosis (Anaplasma marginale) [2, 3]. The morphological and genetic differences among R. microplus strains have been documented in the literature [4, 5]. Although Londt and Arthur  argued that morphological differences between these strains are too slight to warrant species status, other results provided opposing evidence . These authors demonstrated that hybrids resulting from crossing R. (B.) microplus strains from South Africa and Australia were sterile. Thus, R. microplus from South Africa and Australia might be regarded as different species .
Contradictory results have been reported after the analysis of different gene loci in R. microplus strains collected at different locations. It has been  found that the analysis of ribosomal internal transcribed spacer 2 (ITS2) sequences resulted in similar divergence between Australian and South African R. microplus strains and among tick strains from Australia, Brazil, Kenya and South Africa. However, a clear divergence in the acetylcholinesterase sequence between Australian and African/Mexican strains of R. microplus has been detected . The mitochondrial 12S rDNA divergence between American and African R. microplus strains is lower than the divergence between these strains and the R. microplus from Australia and Nepal .
These results suggested that biogeographical and ecological separation between boophilid lineages may have resulted in different species of ticks in America/Africa and those from Australia, collectively known as R. microplus, but being actually separate species . To test the hypothesis of the presence of different species in R. microplus ticks, herein we performed a series of experiments to characterize the the reproductive performance of crosses between ticks from Australia, Africa and America and the genetic diversity of R. microplus strains from Australia, Asia, Africa and America.
Tick crosses and fertility analysis
Tick crosses were conducted using strains originally determined as R. microplus from Argentina (ARG), Mozambique (MOZ) and Australia (AUS, Yeerongpilly strain) (Table 1). In experiment I, calves were infested with larvae of ARG and AUS, each strain separately infesting a cotton sleeve (20 cm diameter) glued (Kamar heat detector adhesive, Kamar, Steamboat Springs, CO, USA) to the shaved dorsum skin of two tick-naïve Holstein calves. After 12–15 days of infestation, sleeves were opened and engorged nymphs attached to the skin were manually removed from each sleeve and held in an incubator at 34°C and 95% RH for molting. Emerged adults, less than 24 hours from their ecdysis were then sorted to form the infestation male × female crosses into new sleeves glued to the dorsum of the same calves from where engorged nymphs were collected. Crosses were made with adult ticks from the same strain (AUS × AUS, and ARG × ARG; homologous crosses), and with adult ticks from different strains (AUS × ARG, and ARG × AUS; heterologous crosses) (Table 2). Each cross consisted of 15–20 males and 20–30 females per sleeve. Two additional sleeves, each containing only 20–30 unfed females from each strain were prepared to be the control of virgin females for each strain.
All engorged females recovered from each sleeve (each separate cross) were individually weighed and left in an incubator at 25°C and 85% RH. The total egg mass produced by each female was weighed and its hatchability was determined as described . The egg production efficiency (EPE) was determined as (weight of eggs/weight of the engorged female) × 100 .
In experiment II, tick-naïve New Zealand rabbits were infested with larvae from ARG, AUS, and MOZ strains inside a cotton sleeve (20 cm diameter) glued to the shaved dorsum skin of the rabbit. At 12–15 days post-infestation, sleeves were opened and engorged nymphs attached to the skin were manually removed from each sleeve and held in an incubator at 34°C and 95% RH for molting. Emerged adults, less than 24 hours from their ecdysis, were sexed to form the infestation male × female crosses in sleeves glued to the dorsum of a calf. A total of 12 sleeves were mounted on the calf to encompass the nine possible crosses between the three strains, as well as three additional sleeves for virgin control females for each strain (Table 3). Each cross consisted of 20–25 males and 25–30 females per sleeve, and procedures for recovery and analysis of engorged females recovered from each sleeve were as described for experiment I.
In order to verify successful reproductive compatibility between strains, unfed larvae (20–30 days old) obtained from the crosses (designated as F1 larvae) were used to infest a calf for experiment I-F1 ticks, and another calf for experiment II-F1 ticks (one sleeve used for infestation with larvae generated from each separate cross). The resulting engorged females were processed as previously described.
Feeding periods of ticks from each cross were compared by the non-parametric Mann-Whitney test, whereas engorged female weights and reproductive parameters (which show normal distribution) were compared by the Student t-test. Variables were considered significantly different if P < 0.05.
Phylogenetic analysis of tick mitochondrial 12S and 16S rDNA
Specimens of sixteen strains of R. microplus from America, Africa, Asia and Oceania and one strain of each R. annulatus and R. decoloratus from Africa were used for DNA extraction and sequencing of mitochondrial 12S and 16S rDNA (Table 1). One to three specimens from each strain were analyzed. DNA was extracted from alcohol-preserved specimens and polymerase chain reaction (PCR) amplification of a fragment of the mitochondrial 16S rDNA was conducted as described . The PCR conditions for 12S rDNA amplification were as described . Amplified PCR products were purified using Wizard PCR Preps DNA Purification System (Promega Corporation, Madison, Wisconsin, USA). The purified DNA was directly sequenced at the IMyZA (Instituto de Microbiología y Zoología Agrícola, INTA, Castelar, Buenos Aires, Argentina). Both DNA strands were sequenced and assembled using BioEdit 7.05.3 software .
For phylogenetic analysis, the following tick mitochondrial 16S and 12S rDNA sequences available in the Genbank were also included:
R. annulatus, Spain (Z97877), R. annulatus, USA (L34311), R. microplus, USA (L34310), R. microplus, strain 9 from Taiwan (AY974232), Dermacentor andersoni, USA (L34299), D. marginatus, Spain (Z97879), D. variabilis, USA (L34300), Haemaphysalis juxtakochi, Uruguay (AY762323), H. leporispalustris, USA (L34309), H. punctata, Spain (Z97880), Hyalomma dromedarii, UK (L34306), H. lusitanicum, Spain (Z97881), Rhipicephalus appendiculatus, UK (L34301), R. bursa, Spain (Z97878), R. pusillus, Spain (Z97883) and R. sanguineus, Spain (Z97884).
R. annulatus (U95866), R. annulatus, Italy (AM410573), R. annulatus, Israel (AF133058), R. microplus, Nepal (AF150042), R. microplus, Australia (AF031847), R. decoloratus, Zimbabwe (AF150044), R. decoloratus, Kenya (AF031846), R. kohlsi, Jordan (AF150043), R. kohlsi, origin unknown (AY008686), Dermacentor albipictus, USA (AF150041), D. andersoni, USA (AF150040), D. reticulatus, France (AF150038), Haemaphysalis leachi, Zimbabwe (AF150035), H. punctata, Switzerland (AF150032), Hyalomma dromedarii, Morocco (AF150036), H. truncatum, Zimbabwe (AF150031), Rhipicephalus appendiculatus, Uganda (AF150028), R. bursa, Spain (AF150053), R. evertsi evertsi, Zimbabwe (AF150052), R. pusillus, France (AF150022), R. sanguineus, France (AF150020), and R. turanicus, France (AF150018).
Multiple sequence alignments were done using Clustal W . Phylogenetic and molecular evolutionary analyses were conducted using MEGA version 4 . Phylogenetic relationships between sequences were assessed by neighbored-joining (NJ) method . The NJ topologies were examined using Tamura-Nei distances  and relative support for the internal nodes was tested by bootstrapping over 1,000 replications . Tamura-Nei distance measures were used because this model corrects for multiple hits, taking into account the differences in substitution rate between nucleotides and the inequality of nucleotide frequencies. It distinguishes between transitional substitution rates between purines and transversional substitution rates between pyrymidines. All positions containing alignment gaps were eliminated only in pairwise sequence comparisons. Ticks species of the genus Haemaphysalis were treated as out-group.
Analysis of microsatellite polymorphism
Microsatellite polymorphism was analyzed in tick strains from Argentina, Australia, Mozambique, India and New Caledonia (Table 1). For Argentina, Australia and Mozambique strains, femalestrain 1 × malestrain 2 and femalestrain 2 × malestrain 1 homologous and heterologous crosses were used. For each cross, egg masses from two separate crosses were used for genotyping. Tick DNA was obtained from egg batches (ARG, AUS and MOZ strains), pooled larvae (New Caledonia strain) or pooled whole ticks (India and Indonesia strains). Tick tissues were homogenized in liquid N or with a 1 ml tuberculin syringe with a 25-G needle to extract DNA with Tri Reagent (Sigma, St. Louis, MO, USA) following manufacturer's recommendations. The purified DNA was dissolved in distilled water and the concentration determined using the NanoDrop 1000 (Thermo Fisher Scientific, Wilmington, DE, USA).
Polymorphism was analyzed at microsatellite loci BmA12 (Genbank accession number DQ001904), BmA06 (DQ001905), BmB12 (DQ001906), BmC03 (DQ001907), BmC07 (DQ001909) and BmD12 (DQ001911) using oligonucleotide primers and PCR cycling conditions described . The PCR was done with labeled forward primers (BmA12 and BmD12, 6-FAM; BmA06 and BmC07, NED; BmB12, VIC; BmC03, PET) in a 50-μl volume (1.5 mM MgSO4, 1X avian myeloblastosis virus (AMV) RT/Thermus flavus (Tfl) reaction buffer, 0.2 mM each deoxynucleoside triphosphate (dNTP), 5 u Tfl DNA polymerase, 0.2 μM of each oligonucleotide primer) employing the Access RT-PCR system (Promega, Madison, WI, USA). Reactions were performed in an automated DNA thermal cycler (Techne model TC-512, Cambridge, England, UK). The reaction was terminated after a final extension at 68°C for 5 min. Control reactions were done using the same procedures, but without DNA added to check contamination of the PCR reaction. PCR products were electrophoresed on 1% agarose gels to check the size of amplified fragments by comparison to a DNA molecular weight marker (1 Kb Plus DNA Ladder, Promega). Fragments were separated using an ABI 3730 automated DNA sequencer (Applied Biosystems, Inc. Foster City, CA, USA) and sized relative to a ROX-labeled internal size standard (GeneScan-500LIZ, Applera, Norwalk, CT, USA). The data were analyzed using program Peak Scanner (Applied Biosystems).
For statistical analysis, a binary matrix was constructed by scoring the different alleles from each microsatellite locus as presence (1) and absence (0) of the PCR bands for each cross. Based on the binary matrix, similarity matrixes were calculated using Jaccard's and Dice's coefficients [23, 24]. The similarity matrix was subjected to Sequential Agglomerative Hieratical Nested Clustering (SAHN) and dendograms were constructed employing the Unweighted Pair Group Method of Arithmetic Averages (UPGMA) of Sneath and Sokal  to group the progenies into clusters. All the analyses were performed using statistical software package NTSYSpc version 2.01e .
The crosses between Australia and Argentina or Mozambique strains are infertile
Tick crosses were conducted between AUS, MOZ and ARG R. microplus strains (Tables 2 and 3). Both homologous and heterologous crosses of ARG and MOZ produced fertile offspring. AUS ticks produced a viable, fertile offspring in homologous crosses. However, AUS females failed to yield viable larvae in heterologous crosses. Even though heterologous crosses of AUS males with ARG or MOZ females produced viable offspring (mean egg hatching varying from 15.6% to 53.6%), the resultant F1 hybrid ticks were infertile in all cases (Table 4). In contrast, F1 ticks from heterologous crosses between ARG and MOZ were as fertile as the F1 ticks from homologous crosses. This result suggested a lack of genetic compatibility between Australian ticks and strains from Africa and America. It was interesting to note the lower weight and feeding period of engorged AUS females, even in homologous crosses (Tables 2 and 4). Regardless of being fertile or not, females from all crosses were significantly heavier and showed feeding periods significantly shorter than their corresponding virgin female controls.
The Australian tick strain is genetically more divergent than Asian, African and American strains
Approximately 400 bp of the mitochondrial 16S rDNA gene were successfully sequenced from specimens of R. microplus representing 16 different strains from America, Africa, Asia and Oceania. The analysis of 16S rDNA sequences variation revealed very little genetic diversity (0.0 to 0.3%) between American and African strains (Fig. 1). In contrast, genetic diversity at this locus was 1.3–1.6% and 1.6–2.4% when these strains were compared to R. microplus from Oceania and Asia, respectively (Fig. 1). For 12S rDNA, approximately 350 bp were sequenced from R. microplus representing 15 strains from America, Africa, Asia and Oceania. No differences were observed in 12S rDNA sequences between American and African strains (Fig. 2). However, when these strains were compared to R. microplus from Oceania and Asia, the genetic diversity was 1.3–1.6%, and 1.6–2.5%, respectively.
The phylogenetic analysis of 16S sequences showed that sequences of R. microplus from America form a clade with African strains and the Taiwan (AY974232) sequence (98% bootstrap support) whereas the strains from Australia, New Caledonia and Indonesia grouped together in a different cluster (Fig. 3). Surprisingly, sequences of R. annulatus from USA and Spain formed a third cluster with R. microplus from India with a relatively high bootstrap support (76%) (Fig. 3). The phylogenetic analysis of 12S rDNA sequences was similar to the analysis of 16S rDNA (Fig. 4). This analysis provided strong support for a clade including R. microplus from America and Africa (98% bootstrap support). The strains from Australia, New Caledonia and Indonesia formed another cluster with a relatively high bootstrap support (78%). Finally, a strong support (99%) was found for a third clade containing R. microplus from India and Nepal (AF150043). In this analysis, the R. annulatus 12S rDNA sequences grouped in a separate cluster with a low bootstrap support (65%).
The analysis of microsatellite loci was conducted to provide additional evidences of genetic polymorphism between tick strains using nuclear genome markers. In these experiments, egg masses from two separate homologous and heterologous crosses were used and allele number and size was determined for each locus (Table 5 and Additional file 1). The UPGMA clustering analysis of microsatellite alleles obtained in the progenies of R. microplus crosses showed that homologous and heterologous crosses involving AUS had lower similarity coefficients when compared to other homologous and heterologous crosses between MOZ and ARG strains (Fig. 5A and 5B). Some microsatellite loci did not provide a clear genotype for the MOZ × MOZ cross and in ticks from India and New Caledonia (Additional file 1), probably due to the presence of null alleles in these strains.
This paper introduced a combined approach based on biological and molecular tools to provide additional support for the hypothesis of the existence of at least two different tick species currently considered as R. microplus. The most important evidences were provided by studies of crosses between tick strains which showed reproductive isolation of the Australian strain. Genetic analyses of mitochondrial and nuclear genome markers were not conclusive but provided support for the genetic divergence between the Australian and African/American strains.
Previous results have suggested that the Australian strains of R. microplus show biological differences when compared to other R. microplus strains from Africa and America . In these studies, sterile progenies were obtained from crosses between strains from Australia and South Africa. In addition, the mean number of eggs produced by engorged females of American strains is higher than that obtained from Australian strains [27–30]. Furthermore, the amount of the insect growth regulator fluazuron required for complete growth inhibition in Australian tick strains is lower than that required for the Argentinean strains . These results suggested that Australian R. microplus differ from strains in America and Africa and may have evolved into a separate species.
To test this hypothesis, we conducted experiments to characterize the reproductive performance and genetic diversity of R. microplus strains from Australia, Africa and America. The results clearly showed that tick strains from Africa (Mozambique) and America (Argentina) have no biological barrier for reproduction, being therefore the same species. All the reproductive parameters were significantly reduced when the Australian strain was involved in heterologous crosses suggesting the existence of a separate species. The genetic analysis of mitochondrial 12S and 16S rDNA as well as nuclear microsatellite loci sequences further suggested that strains from Africa and South America are conspecific. According to these analyses, ticks from Australia, Indonesia and New Caledonia would be also conspecific, different from the American-African strains. While clustered in a different clade, support for a separate species in India and Nepal was inconclusive. Interestingly, the results of the microsatellite analysis suggested the absence of transmission ratio distortion between paternal and maternal alleles even in the offspring of the AUS crosses with either ARG and/or MOZ strains, reducing the possibility of parental genetic effects to explain strain incompatibility.
It is well known that in the absence of males, Metastriata females increase their feeding period, have smaller repletion weights and oviposit none or only few fertile eggs when compared to females mated with conspecific males . In the present study, females from all crosses showed feeding periods significantly shorter and significantly higher repletion weights than the virgin female controls. This result indicated that all crosses resulted in copula regardless of being fertile or not and evidenced that AUS females were attractive to ARG and MOZ males and vice versa. Attractiveness and cross-mating between R. microplus and R. annulatus has been reported, resulting in sterile hybrids , similarly to our results for heterologous crosses with AUS ticks.
One possible drawback of this study is the use of the Australian Yeerongpilly strain, which has been kept in the laboratory for over 60 years . This strain has been also used in previous analyses of the reproductive success of heterologous crosses . However, while low reproductive performance could be ascribed to continuous maintenance of the strain in the laboratory, differences were not observed in 12S and 16S rDNA sequences between the Yeerongpilly strain, the other two Australian strains analyzed in this study (A and B; Tables 4 and 5) and those reported by Murrell et al. . Therefore, reproductive isolation seems to be a factual character associated to the strain and not a side effect of long-term laboratory maintenance.
Although mitochondrial 12S and 16S rDNA phylogenies were constructed with sequences from multiple Australian strains, the studies of the reproductive performance and microsatellite polymorphism were conducted with a single Australian (Yeerongpilly) strain, a fact that should be taken with caution when generalizing these results to Australian R. microplus in general. However, despite the limitations of the study, the results reported herein support the hypothesis that at least two species exist under the name R. microplus. These species could be redefined as R. microplus (Canestrini, 1887) (for American and African strains) and probably R. australis Fuller, 1899 (for Australian strains). After its description, R. australis was reported in Africa. However, these records were synonymized under R. annulatus. Interestingly, R. annulatus appeared within the monophyletic clade of R. microplus in the 16S mtDNA tree or as a paraphyletic branche in the 12S mtDNA tree suggesting that the denomination of R. annulatus may needs revision. Currently, R. australis is regarded as a synonym of R. microplus. Biological features reported here, including the lack of fertile F1 progeny in heterologous crosses and smaller, significantly different weight of engorged Australian (Yeerongpilly strain) females as well as molecular findings involving different Australian strains are clearly supportive of the hypothesis of a reproductive isolation between Australian and American/African ticks. We can not ascribe the Australian-Indonesian-Caledonian specimens to that species until the finding of its types and further morphological comparisons are conducted. It is also possible that other boophilid tick strains may inhabit small areas of the South Pacific. This issue warrants further studies by conducting detailed morphologic, genetic and reproductive studies with multiple tick strains collected at different geographic locations and cured from their endosymbionts which may affect tick reproductive performance.
It is important to note that the introduction of R. microplus into Africa took place with animals from India, after the epidemic of Rinderpest in the 19th century. Although it is hypothesized that R. microplus entered into South America along with Indian cattle, information is currently not available to support this hypothesis. However, the introduction of R. microplus in Australia took place with animals from Indonesia (most probably cattle from Java around 1870)  and then later into New Caledonia with Australian cattle (Barré, pers. comm.). This fact is in agreement with the reproductive and molecular findings of this study: ticks from Indonesia, Australia and New Caledonia may be a different, reproductively isolated species. As mentioned before, the strains from India and Nepal may constitute a different species, closer to African and American R. microplus, but adequate biological evidence is lacking and molecular evidence presented herein are not conclusive and/or may be the result of allopatric evolution. Detailed morphological, genetic and reproductive studies of tick strains from Africa, America, Australia, Indonesia, New Caledonia, India and Nepal would be necessary before confirming the redescription and possible resurrection of other boophilid tick species. Nevertheless, until further studies are conducted, caution should be taken in extrapolating results from studies on Australian ticks determined as R. microplus. These data have important implications in the field of animal health, because many studies involving ecology and control of cattle ticks performed on Australian tick strains considered as R. microplus may not be directly applicable to tick strains in America or Africa.
The results reported herein provided biological and genetic evidences of allopatric speciation in R. microplus and suggested a reproductive isolation for Australian-Indonesian-Caledonian strains resulting in a different species. These results support the hypothesis that at least two different species have evolved under the name R. microplus. These species must be redefined after careful inspection of type specimens and further analysis of local tick fauna, including studies with several tick strains from each geographic location.
Estrada-Peña A, Bouattour A, Camicas JL, Guglielmone A, Horak I, Jongejan F, Latif A, Pegram R, Walker AR: 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. 10.1007/s10493-006-0003-5.
Peter RJ, Bossche Van den P, Penzhorn BL, Sharp B: Tick, fly, and mosquito control-Lessons from the past, solutions for the future. Vet Parasitol. 2005, 132: 205-215. 10.1016/j.vetpar.2005.07.004.
de la Fuente J, Estrada-Peña A, Venzal JM, Kocan KM, Sonenshine DE: Overview: Ticks as vectors of pathogens that cause disease in humans and animals. Frontiers in Biosciences. 2008, 13: 6938-6946. 10.2741/3200.
García-García JC, González IL, González DM, Valdés M, Méndez L, Lamberti J, D'Agostino B, Citroni D, Fragoso H, Ortiz M, Rodríguez M, de la Fuente J: Sequence variations in the Boophilus microplus Bm86 locus and implications for immunoprotection in cattle vaccinated with this antigen. Exp Appl Acarol. 1999, 23: 883-895. 10.1023/A:1006270615158.
de la Fuente J, García-García JC, González DM, Izquierdo G, Ochagavia ME: Molecular analysis of Boophilus spp. (Acari: Ixodidae) tick strains. Vet Parasitol. 2000, 92: 209-222. 10.1016/S0304-4017(00)00308-3.
Londt JGH, Arthur D: The structure and parasitic life cycle of Boophilus microplus (Canestrini, 1888) in South Africa (Acarina: Ixodidae). J Entomol Soc S Afr. 1975, 38: 321-340.
Spickett AM, Malan AR: Genetic incompatibility between Boophilus decoloratus (Koch 1844) and Boophilus microplus (Canestrini, 1888) and hybrid sterility of Australian and South African Boophilus microplus (Acarina: Ixodidae). Onderstepoort J Vet Res. 1978, 45: 149-153.
Guglielmone AA, Estrada-Peña A, Keirans JE, Robbins RG: Ticks (Acari: Ixodidae) of the Neotropical Zoogeographic Region. 2003, Special Publication, International Consortium on Ticks and Tick-borne Diseases, Atalanta, Houten, The Netherlands, 173-
Barker SC: Distinguishing species and strains of Rhipicephaline ticks with ITS2 ribosomal RNA. J Parasitol. 1998, 84: 887-892. 10.2307/3284614.
Baxter GD, Barker SC: Comparison of acetylcholinesterase genes from cattle ticks. Int J Parasitol. 1999, 29: 1765-1774. 10.1016/S0020-7519(99)00136-8.
Nava S, Guglielmone AA, Mangold AJ: An overview of systematic and evolution of ticks. Front Biosciences. 2009, 14: 3012-3023.
Drummond RO, Ernst SE, Trevino JL, Gladney WJ, Graham OH: Boophilus annulatus and B. microplus : laboratory tests of insecticides. J Econ Entomol. 1973, 66: 130-133.
Bennett GF: Oviposition of Boophilus microplus (Canestrini) (Acarina: Ixodidae). Acarologia. 1974, 16: 1652-1661.
Mangold AJ, Bargues MD, Mas-Coma S: Mitochondrial 16S rDNA sequences and phylogenetic relationships of species of Rhipicephalus and other tick genera among Metastriata (Acari: Ixodidae). Parasitol Res. 1998, 84: 478-484. 10.1007/s004360050433.
Szabó MPJ, Mangold AJ, Joao CF, Bechara GH, Guglielmone AA: Biological and DNA evidence of two dissimilar populations of the Rhipicephalus sanguineus tick group (Acari: Ixodidae) in South America. Vet Parasitol. 2005, 130: 131-140. 10.1016/j.vetpar.2005.03.008.
Hall TA: BioEdit: a user-friendly biological sequence alignment editor and analysis program for Windows 95/98/NT. Nucl Acids Symp. 1999, 41: 95-98.
Thompson JD, Higgins DG, Gibson TJ: Improving the sensitivity of progressive multiple sequence alignment through sequence weighting, position-specific gap penalties and weight matrix choice. Nucl Acids Res. 1994, 22: 4673-4680. 10.1093/nar/22.22.4673.
Tamura K, Dudley J, Nei M, Kumar S: MEGA4: Molecular Evolutionary Genetics Analysis (MEGA) software version 4.0. Mol Biol Evol. 2007, 24: 1596-1599. 10.1093/molbev/msm092.
Saitou N, Nei M: The neighbor-joining method: A new method for reconstructing phylogenetic trees. Mol Biol Evol. 1987, 4: 406-425.
Tamura K, Nei M: Estimation of the number of nucleotide substitutions in the control region of mitochondrial DNA in humans and chimpanzees. Mol Biol Evol. 1993, 10: 512-526.
Felsenstein J: Confidence limits on phylogenies: An approach using the bootstrap. Evolution. 1985, 39: 783-791. 10.2307/2408678.
Koffi BB, Risterucci AM, Joulia D, Durand P, Barré N, de Meeús T, Chevillon C: Characterization of polymorphic microsatellite loci within a young Boophilus microplus metapopulation. Mol Ecol Notes. 2006, 6: 502-504. 10.1111/j.1471-8286.2006.01295.x.
Jaccard P: Nouvelles recherches sur la distribution florale. Bull Soc Vaud Sci Nat. 1908, 44: 223-270.
Dice LR: Measures of the amount of ecologic association between species. Ecology. 1945, 26: 297-302. 10.2307/1932409.
Sneath PHA, Sokal RR: Numerical taxonomy: the principles and practice of numerical classification. 1973, W.H. Freeman & Co., San Francisco, 573-
Rohlf FJ: Adaptive hierarchical clustering schemes. Systematic Zool. 1970, 19: 58-82. 10.2307/2412027.
Legg J: Some observations on the life history of the cattle tick (Boophilus australis). Proc Roy Soc Queensland. 1930, 41: 121-132.
Hitchcock LF: Studies on the parasitic stages of the cattle tick Boophilus microplus. Can Aust J Zool. 1954, 3: 145-155. 10.1071/ZO9550145.
Alvarado RU, González JC: A postura e a viabilidade do Boophilus microplus (Canestrini, 1888) (Acarina, Ixodidae) en condicioes de laboratorio. Rev Lat Amer Microbiol. 1979, 21: 31-36.
Guglielmone AA, Mangold AJ, Aguirre DH, Gaído AB, Olsen AA: The effect of infection by Babesia sp. on some biological parameters of engorged females of Boophilus microplus. Folia Parasitol. 1989, 36: 1-6.
Guglielmone AA, Volpogni MM, Anziani OS, Hereu C: Evaluación del fluazuron (Acatak) para el control de poblaciones naturales de Boophilus microplus (Acari: Ixodidae) en Santiago del Estero, Argentina. Vet Arg. 1998, 15: 710-714.
Oliver JH: Biology and systematics of ticks (Acari: Ixodida). Ann Rev Ecol Syst. 1989, 20: 397-430. 10.1146/annurev.es.20.110189.002145.
Graham OH, Price MA, Trevino JL: Cross-mating experiments with Boophilus annulatus and B. microplus (Acarina: Ixodidae). J Med Entomol. 1972, 9: 531-537.
Hitchcock LF, Roulston WJ: Arsenic resistance in a strain of cattle tick (Boophilus microplus (Canestrini)) from northern New South Wales. Aust J Agricult Res. 1955, 6: 666-671. 10.1071/AR9550666.
Murrell A, Campbell NJH, Barker SC: Mitochondrial 12S rDNA indicates that the Rhipicephalinae (Acari: Ixodida: Ixodidae) is paraphyletic. Mol Phylogen Evol. 1999, 12: 83-86. 10.1006/mpev.1998.0595.
Roberts FHS: Australian ticks. 1970, CSIRO, Melbourne, 267-
We thank Peter Willadsen (CSIRO Livestock Industries, Queensland, Australia) for the information about the Yeerongpilly tick strain, Martín N. García (Intituto de Recursos Naturales, INTA, Castelar, Argentina) for technical assistance and Wilson R. Fernandes (Faculty of Veterinary, University of São Paulo) for logistic support for tick infestations on cattle. S. Ghosh (Entomology Laboratory, Division of Parasitology, Indian Veterinary Research Institute, Izatnagar, India) is acknowledged for supplying the Indian tick strain. Members of the Integrated Consortium on Ticks and Tick-borne Diseases (ICTTD-3) are acknowledged for fruitful discussions and suggestions. This work was supported by ICTTD-3, financed by the International Cooperation Program of the European Union, coordination action project No. 510561, the Consejería de Educación y Ciencia, JCCM, Spain (project PAI06-0046-5285) (to JF), the Conselho Nacional de Desenvolvimento Científico e Tecnológico – CNPq (to MBL), the Consejo Nacional de Investigaciones Cientícas y Técnicas de Argentina (PIP 2058) and INTA Rafaela (TCP 426100) (to AJM, AAG and CT). V. Naranjo was funded by Junta de Comunidades de Castilla – La Mancha (JCCM), Spain.
MBL conducted the tick crosses. MBL and AE-P analyzed the results of tick crosses. VN and AJM obtained new tick sequences. AJM and AAG conducted phylogenetic analyses. VN, CT and JF obtained and analyzed tick microsatellites. MBL, AE-P, AAG, FJ and JF conceived and participated in study design. MBL, AE-P, AAG, AJM, FJ and JF helped to draft the manuscript. All authors read and approved the final manuscript.