New insights from molecular characterization of the tick Rhipicephalus (Boophilus) microplus in Brazil

Original Article Braz. J. Vet. Parasitol., Jaboticabal, v. 25, n. 3, p. 317-326, jul.-set. 2016 ISSN 0103-846X (Print) / ISSN 1984-2961 (Electronic) Doi: http://dx.doi.org/10.1590/S1984-29612016053 New insights from molecular characterization of the tick Rhipicephalus (Boophilus) microplus in Brazil Novos pontos de vista sobre a caracterização molecular em carrapatos Rhipicephalus (Boophilus) microplus, no Brasil Bárbara Guimarães Csordas1; Marcos Valério Garcia2§; Rodrigo Casquero Cunha3; Poliana Fernanda Giachetto4; Isabella Maiumi Zaidan Blecha5; Renato Andreotti2* 1 Pós-graduação em Doenças Infecciosas e Parasitárias, Faculdade de Medicina, Universidade Federal de Mato Grosso do Sul – UFMS, Campo Grande, MS, Brasil 2 Laboratório de Biologia Molecular do Carrapato, Departamento de Sanidade Animal, Embrapa Gado de Corte, Campo Grande, MS, Brasil 3 Programa de Pós-graduação em Biotecnologia, Universidade Federal de Pelotas – UFPel, Pelotas, RS, Brasil Laboratório Multiusuário de Bioinformática, Embrapa Informática Agropecuária, Campinas, SP, Brasil 4 5 Programa de Pós-graduação em Ciência Animal, Faculdade de Medicina Veterinária, Universidade Federal de Mato Grosso do Sul – UFMS, Campo Grande, MS, Brasil § Bolsista DCR Fundect/MS. Received April 11, 2016 Accepted July 20, 2016 Abstract The Rhipicephalus (Boophilus) microplus complex currently consists of five taxa, namely R. australis, R. annulatus, R. (B.) microplus clade A sensu, R. microplus clade B sensu, and R. (B.) microplus clade C sensu. Mitochondrial DNA-based methods help taxonomists when they are facing the morpho-taxonomic problem of distinguishing members of the R. (B.) microplus complex. The purpose of this study was to perform molecular characterization of ticks in all five regions of Brazil and infer their phylogenetic relationships. Molecular analysis characterized 10 haplotypes of the COX-1 gene. Molecular network analysis revealed that haplotype H-2 was the most dispersed of the studied populations (n = 11). Haplotype H-3 (n = 2) had the greatest genetic differentiation when compared to other Brazilian populations. A Bayesian phylogenetic tree of the COX-1 gene obtained strong support. In addition, it was observed that the population of R. (B.) microplus haplotype H-3 exhibited diverging branches among the other Brazilian populations in the study. The study concludes that the different regions of Brazil have R. (B.) microplus tick populations with distinct haplotypes. Keywords: Rhipicephalus (Boophilus) microplus, cattle tick, COX-I gene, ITS-2 gene, Brazil. Resumo Carrapatos do complexo R. (B.) microplus se distribuem em cinco taxa: R. australis, R. annulatus, R. (B.) microplus clado A sensu R. microplus clado B sensue e R. (B.) microplus clado C sensu. Métodos baseados no DNA mitocondrial podem auxiliar taxonomistas quando há dificuldades em estabelecer diferenças morfológicas para distinguir membros do complexo R. (B.) microplus. O objetivo deste estudo foi a caracterização molecular e a inferência de relações filogenéticas em carrapatos de todas as cinco regiões geográficas do Brasil. Para o gene COX-1, a análise molecular caracterizou 10 haplótipos. Na análise molecular em rede foi observado que o haplótipo H-2 é o mais disperso entre as populações (n=11). O haplótipo H-3 (n=2) foi o que obteve maior diferenciação genética ao ser comparado com outras populações brasileiras. A árvore filogenética Bayesiana de gene COX-1 gerou suporte robusto e foi observado que a população de R. (B.) microplus haplótipo H-3 apresentou ramificação com divergência entre as outras populações brasileiras apresentadas neste estudo. Conclui-se que as populações brasileiras possuem diversidade haplotípica com divergência entre as diversas populações de R. (B.) microplus no Brasil. Palavras-chave: Rhipicephalus (Boophilus) microplus, carrapato-do-boi, gene COX-I, gene ITS-2, Brasil. *Corresponding author: Renato Andreotti. Departamento de Sanidade Animal, Embrapa Gado de Corte, Avenida Rádio Maia, 830, Vila Popular, CEP 79106-550, Campo Grande, MS, Brasil. e-mail: renato.andreotti@embrapa.br www.cbpv.org.br/rbpv 318 Csordas, B.G. et al. Introduction The bovine tick Rhipicephalus (B.) microplus can be found in multiple tropical and subtropical regions worldwide (ESTRADA-PEÑA et al., 2006). This parasite may transmit pathogens that cause babesiosis (Babesia bovis and B. bigemina) and anaplasmosis (Anaplasma marginale) (PETER et al., 2005; de la FUENTE et al., 2008). Infestation of herds with R. (B.) microplus cause losses estimated at US$3.24 billion per year in Brazil (GRISI et al., 2014). In most cattle-producing countries, control methods against this tick are costly and require the use of acaricides. However, acaricide-resistant R. (B.) microplus populations have become a worldwide problem, and molecular ecology studies and new tick control technologies are now required to preserve cattle production (RODRÍGUEZ-VIVAS et al., 2007; ANDREOTTI et al., 2011; GUERRERO et al., 2014). Six years ago, Labruna et al. (2009) used 12S and 16S molecular markers to show a lack of divergence among populations in the Americas, including Brazil. Later, Burger et al. (2014) used the COX-1 molecular marker to show that the R. (B.) microplus populations of Brazil (GenBank: KC503261) are organized in clade A, along with populations of R. (B.) microplus from China (GenBank: HM193863) and Cambodia (GenBank: KC503260). Thus, molecular markers may resolve very close taxonomic relationships and provide new knowledge about population structure, allowing a deeper understanding of R. (B.) microplus populations (BURGER et al., 2014). Although the history of the dissemination of R. (B.) microplus is not well documented, the species is known to have originated in India (HOOGSTRAL, 1986). According to Barré & Uilenberg (2010), R. (B.) microplus originated in the southern and southeastern regions of Asia and was spread throughout the tropical and subtropical belts via cattle, arriving in Brazil between the 16th and 17th centuries. Mitochondrial and nuclear genome markers are increasingly being used to better understand phylogenetic relationships in order to elucidate whether evolution has occurred among populations of a species (KANDUMA et al., 2012). The DNA barcode has been proposed as a universal tool for identifying biological diversity, and can be used as a molecular marker (HEBERT et al., 2003). DNA barcoding is based on information gathered from a fragment of approximately 688 base pairs of mitochondrial DNA base sequences (mtDNA) from the cytochrome oxidase I gene (COX-I) of different species. A number of studies have shown that the DNA barcode is a universal code highly effective for the identification of species (HEBERT et al., 2004; BARRETT & HEBERT, 2005). In addition to the DNA barcode, ribosomal DNA (rDNA) may be used as a nuclear molecular marker. The internal transcribed spacer 2 (ITS-2) is located between the 5.8S and 28S ribosomal subunits of rDNA (CRUICKSHANK, 2002). These genes are arranged in repeated units known as ribosomal cistrons, which have repeated copies and have been used to study phylogenetic Braz. J. Vet. Parasitol. relationships (CAMPBELL et al., 1993; SONG et al., 2011; BURGER et al., 2014). The objective of this study was to infer the phylogenetic and phylogeographic relationships among R. (B.) microplus tick populations in Brazil based on COX-I mitochondrial DNA gene sequences and ITS-2 nuclear DNA, and compare them with R. (B.) microplus tick populations from other countries. Materials and Methods Tick collection Ticks were obtained from 22 locations in different geographical regions of Brazil (Figure 1, Table 1). A pool of larvae and 22 engorged females were collected from natural environments. The morphology and identification key was based on Barros-Battesti et al. (2006). Samples were stored in an ultra-freezer at –80 °C for subsequent analysis. DNA extraction and PCR amplification DNA extraction was performed on a pool of larvae and individual samples of engorged females using DNAzol Reagent (Invitrogen) according to the manufacturer’s recommendations. The DNA was quantified using a Nano Drop 1000 spectrophotometer (Thermo Scientific, Wilmington, DE, USA), and an A260nm/A280nm ratio above 1.6 was established as acceptable. The sample concentration was corrected to 50 ng µL–1. Next, the DNA was used in a polymerase chain reaction (PCR) for COX-I and ITS-2 region amplification with new primers designed by our group. PCR generated products of 643bp and 580bp for COX-1 and ITS-2, respectively. COX-I gene fragments were amplified using the COX-I.2F (5’-CTTCAGCCATTTTACCGCGA-3’) and COX-I.2R (5’-CTCCGCCTGAAGGGTCAAA-3’) starter oligonucleotides. ITS-2 fragments were amplified with the ITS-2F (5’-CGGATCACATATCAAGAGAG-3’) and ITS-2R (5’-CCCAACTGGAGTGGCCCAGTTT-3’) primers. The PCR was standardized for a final volume of 25 µL with 1X buffer, 1.5 mM MgCl, 25 mM dNTPs, 10 µM of each primer, 0.5 unit of Platinum Taq DNA Polymerase High Fidelity (Invitrogen by Life TechnologiesTM, Massachusetts, USA), and 1µL of DNA at 50 ngµL-1. PCR conditions were optimized for each reaction, and the annealing temperature was adjusted to suit the primers used. General PCR conditions were 94 °C for 2 min, followed by 35 cycles of 94 °C for 1 min, 56 °C/60 °C (COX-I/ITS-2, respectively) for 30s, 68 °C for 1 min, and a final extension of 68 °C for 3 min. PCR products were visualized in a 2% agarose gel stained with ethidium bromide. The PCR-amplified products were purified using a Purelinkkit (Invitrogen by Life TechnologiesTM, Massachusetts, USA). Fragments were then cloned with pGEM-T Easy Vector Systems cloning vector according to the manufacturer’s manual (Promega, Madison, WI-USA) and sequenced as described below. v. 25, n. 3, jul.-set. 2016 Molecular characterization of Rhipicephalus (Boophilus) microplus in Brazil 319 Figure 1. Geographic locations where Rhipicephalus (Boophilus) microplus ticks were collected in five regions of Brazil. Genetic analysis All analyses generated in this study were performed in silico. Sanger sequencing was performed in a 48 well plate with robotic instrumentation using T7 Transcription Start (5’-CTAATACGACTCACTATAGGG-3’) universal sequencing primer (Promega, Madison, WI-USA). Samples were sequenced using an Applied BiosystemsTM ABI 3730 DNA Analyzer and the conditions described by Song et al. (2011). The sequencing reactions were done using the Big Dye Terminator v3.1 cycle sequencing kit. The runs were performed in 36-cm capillaries using the POP7 polymer, and the sequences were generated by the Sequencing Analysis Software v5.3.1 through the Caller KB Phred program (EWING & GREEN, 1998; EWING et al., 1998). Plasmid sequences were identified and removed using an NCBI (VecScreen) tool (ALTSCHUL et al., 1997), and a hash search algorithm was used to remove the contaminant sequences. Consensus sequences were aligned with the GenBank data using the BLASTN program (ALTSCHUL et al., 1997). The partition homogeneity test was conducted using the ARLEQUIN program, version 3.5.1.2 (EXCOFFIER & LISCHER, 2010). The aligned COX-I and ITS-2 sequences comprised 681bpsand 649bps, respectively. Database sequence R. (B.) microplus MS- Brazil (GenBank: KC503261) clade A sensu Burger et al. (2014), was used for alignment with samples from this study. This sequence alignment revealed 10 haplotypes (H1-H10) based on COX-I genes alone (Figure 2). Distinct sequences of COX-I (KP226159-KP226180) and ITS-2 (KP226139-KP226158) were deposited in NCBI GenBank. A median-joining analysis implemented in the program Network Version 5.0 (BANDELT et al., 1999) was used for the intraspecific analysis of the evolutionary relationships among haplotypes. Uncorrected (p) pairwise genetic distances were calculated using ARLEQUIN program, version 3.5.1.2 Csordas, B.G. et al. 320 Braz. J. Vet. Parasitol. Table 1. Brazilian states, geographic coordinates related to the collection and GenBank accession number of this study. Region North Northeast Midwest Southeast South State Geographic coordinates AC AM PA RO RR TO MA CE RN PB PI BA MT GO MS ES MG RJ SP PR SC RS 9°59’30”S/67°48’36”W 3°6’0”S/60°1’0”W 1°17’49”S/47°55’19”W 8°45’43”S/63°54’14”W 2°49’10”S/60°40’17”W 10°11’4”S/48°20’01”W 4°51’32”S/43°21’21”W 3°43’6”S/38°32’34”W 5°47’42”S/3512’32”W 6°57’46”S/35°41’31”W 9°4’26”S/44°21’32”W 12°58’16”S/38°30’39”W 13°49’44”S/56°04’56”W 15°47’56”S/47°52’00”W 20°26’34”S/54°38’45”W 20°45’50”S/41°31’58”W 18°55’8”S/48°16’37”W 21°38’24”S/41°3’3”W 21°15’18”S/48°19’19”W 25°23’42”S/51°27’28”W 27°20’27”S/51°37’4”W 31°46’19”S/52°20’34”W (EXCOFFIER & LISCHER, 2010) to assess the genetic divergence of R. (B.) microplus in both COX-I and ITS-2 genes. To assess the level of genetic differentiation, we used the program DnaSP 5.0 (LIBRADO & ROZAS, 2009). Haplotype diversity (Hd), nucleotide diversity (pi), number of haplotypes (h), and genetic differentiation (FST) were determined. The levels of genetic differentiation were defined as FST > 0.25 (great differentiation), 0.15 to 0.25 (moderate differentiation) and FST< 0.05 (negligible differentiation) (WRIGHT, 1931). The levels of gene flow were defined as Nm> 1 (high gene flow), 0.25 to 0.99 (intermediate gene flow) and Nm< 0.25 (low gene flow). The distinct COX-I and 16S rRNA haplotypes identified in the Brazilian R. (B.) microplus ticks were aligned with all representative sequences of Rhipicephalus taxa available in the NCBI GenBank. R. annulatus was used as an outgroup for the construction of Bayesian phylogenetic trees based on Brazilian R. (B.) microplus ticks COX-I only; Dermacentor nitens was used as an outgroup for the construction of Bayesian phylogenetic trees based on COX-I and ITS-2 sequences. The Bayesian phylogenetic analysis performed using MrBayes 3.1 (HUELSENBECK & RONQUIST, 2001; RONQUIST & HUELSENBECK, 2003) were plotted using FigTree software version 1.4.2 (TREE BIO, 2016). The ARLEQUIN software revealed population structure by means of AMOVA, which uncovers the existence of population differentiation at both intra- and inter-population levels (EXCOFFIER & LISCHER, 2010). GenBank (COX-1) KP226159 KP226162 KP226163 KP226178 KP226161 KP226180 KP226160 KP226177 KP226172 KP226176 KP226174 KP226168 KP226167 KP226166 KP226164 KP226175 KP226171 KP226165 KP226170 KP226179 KP226173 KP226169 GenBank (ITS2) KP226139 KP226153 KP226142 KP226146 KP226157 KP226140 KP226145 KP226141 KP226148 KP226156 KP226143 KP226151 KP226155 KP226154 KP226150 KP226149 KP226158 KP226147 KP226144 KP226152 Results and Discussion Nucleotide analyses and Haplotype The phylogenetic and phylogeographic relations of the Brazilian R. (B.) microplus were inferred from the COX-I sequences (22 samples), R. (B.) microplus Brazil (GenBank: KC503261) and the ITS2 sequences (20 samples) obtained from all five Brazilian geographic regions (Figure 1). This is the only study that has compared a representative sample of R. (B.) microplus, which together represent the whole country and the authors were careful to represent each region by just having differentiated biogeography. Assays performed by Dantas-Torres (2015), can demonstrate that certain factors such as climate change and biodiversity boost the expansion of tick populations. The COX-I dataset showed haplotype diversity (0.767) greater than that revealed by the ITS-2 (0.567) gene. However, the ITS-2 showed greater nucleotide diversity (0.00277) than the COX-I (0.00432) dataset. Haplotype network analysis (COX-I) revealed ten distinct haplotype clusters among Brazilian populations, with no clear separation by tick populations/geographical areas, indicating that an overlap of the ten genetically divergent haplotypes of R. (B.) microplus (Figure 2A) exists. Haplotype H-2 was the most widespread haplotype (n = 11) in Northern, Midwest, Southeast and Southern populations (Figure 2B). Haplotypes H-5, H-6 and H-8 (Northern region) originates from H-4 and H-2 haplotypes, as shown in Figure 2B. v. 25, n. 3, jul.-set. 2016 Molecular characterization of Rhipicephalus (Boophilus) microplus in Brazil 321 Figure 2. (A) Bayesian phylogenetic tree recovered for cytochrome oxidase subunit 1 gene (COX-I) of 643 nucleotides. R. annulatus (GenBank: KC503256) was used as the outgroup (B) median joining network of R. (B.) microplus from five different regions populations in Brazil. A circle represents each haplotype. Relative areas of the circles indicate haplotype frequency of each sample. Csordas, B.G. et al. 322 Braz. J. Vet. Parasitol. regions can only be answered by analyzing a larger number of samples from various municipalities in the States of Rondônia and Paraíba. All haplotypes are derived from haplotype H-2. However, it is not possible to assert where haplotype H-2 has originated. Additionally, there is no haplotype H-2 in northeastern Brazil. The fact that the country has cattle production for imported and local breeds and that animals are transported between different regions do not allow us to infer the exact origin of haplotypes. The natural dispersal ability of ticks might occur in parallel with human-mediated dispersal, causing genetic admixture and high gene flow in Brazilian populations. However, how these haplotypes are associated in other regions is not clear. ITS-2 genotypes indicate a lack of clear separation of tick populations (data no shown). Genetic distance, genetic differentiation and gene flow In the uncorrected “p” distance matrix, the COX-I gene indicated stronger resolving power than ITS-2 (1.36–11.00) (Table 2). The genetic distance of the COX-I H-2 haplotype of R. (B.) microplus indicated 0.73% of divergence within this population (red value, Table 2). Overall, a high level of significant genetic differentiation was observed among all populations (P < 0.001). Great genetic differentiation (FST = 1.00), revealed by the COX-I gene, was found with in all haplotypes (Table 3). Additionally, genetic differentiation between H-2 haplotypes and others varied between 0.47-0.89 (Table 3). High levels (30.00) of gene flow (Nm) were observed among tick populations as shown by the COX-I (data no shown). Among the ten haplotypes studied, R. (B.) microplus haplotype H-3 from Rondônia (RO) and Paraíba (PB) were genetically diverse as supported by the highest haplotype and nucleotide diversity based on COX-I (Table 3). Higher inter-population levels (82.69%) and (FST = 0.8269) revealed the existence of population structure in all haplotypes analyzed by means of AMOVA (Table 4). Livestock transition in Brazil is very intense. Whether these regions actually have different populations not found in other Phylogenetic analyses Bayesian analyses produced phylogenetic trees with high valued probabilities (Figure 2 and Figure 3). The 10 haplotypes (H1-H10) generated from the COX-I sequences of Brazilian R. (B.) microplus were subjected to phylogenetic analyses. By the topology of the phylogenetic trees generated, the R. (B.) microplus Brazilian population comprises two clades: clade A (H3) and clade B (H1, H2 and H4 to H10) (Figure 2A). Different than results revealed by the haplotype network analysis (Figure 2B), according to Bayesian phylogenetic tree, haplotypes H5, H6 and H8 were associated with ticks from the North region of the country and are derived from haplotype H-4. The tree also indicates that all haplotypes derive from H-2. Table 2. Percentage of uncorrected “p” distance matrix among the ten haplotypes of COX-1 of R. (B.) microplus in Brazil. H_1 H_2 H_3 H_4 H_5 H_6 H_7 H_8 H_9 H_10 H_1 3.18 4.00 4.00 7.00 5.00 7.00 8.00 4.00 5.00 H_2 H_3 H_4 H_5 H_6 H_7 H_8 H_9 H_10 0.73 1.36 1.36 4.36 2.36 4.36 6.82 1.36 2.27 2.00 5.00 3.00 5.00 8.00 2.00 3.00 3.00 1.00 5.00 6.00 2.00 3.00 4.00 8.00 9.00 5.00 6.00 6.00 7.00 3.00 4.00 11.00 5.00 6.00 8.00 8.00 3.00 - Table 3. Genetic differentiation (FST) based on COX-1, in all the haplotypes found of R. (B.) microplus populations in Brazil. H_1 H_2 H_3 H_4 H_5 H_6 H_7 H_8 H_9 H_10 H_1 0.77143 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 H_2 H_3 H_4 H_5 H_6 H_7 H_8 H_9 H_10 0.5589 0.59878 0.83333 0.69231 0.83333 0.89333 0.46667 0.68000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 - Molecular characterization of Rhipicephalus (Boophilus) microplus in Brazil v. 25, n. 3, jul.-set. 2016 The ITS-2 Bayesian phylogenetic tree revealed one R. (B.) microplus genetic clade (Figure 4). The COX-I Bayesian phylogenetic tree revealed two genetic clades of R. (B.) microplus, Clade A (Brazil, Texas, Malaysia, R. australis) and Clade B (R. (B.) microplus China and R. annulatus). Rhipicephalus (B.) microplus clade A showed a sister relationship with R. australis. By contrast, R. (B.) microplus clade B is more closely related to R. annulatus, with which it forms a sister group relationship, corroborating the data by Burger et al. (2014) (Figure 3). Table 4. Analysis of molecular variance of molecular marker, gene COX-1. Source of variation d.f Sum of squares Among populations Within populations Total Fixation Index FST: 9 13 22 0.8269 Percentage of variation 25.059 Variance Components 1.34Va 3.636 0.28Vb 17.31 28.696 P < 0.001 1.62 82.69 323 The R. (B.) microplus populations in Brazil have internal nodes, indicating the occurrence of subpopulations (GenBank: KP226160, KP226167, KP226169, KP226172, KP226174, KP226177). The haplotype H-3 formed another smaller node (GenBank: KP226176, KP226178) (Figure 3). The intraspecific genetic distance or haplotype frequency, like that of haplotype H2 reported here, are notably higher than previously described for R. (B.) microplus sensu Low et al. (2015). Furthermore, the two, and perhaps three, genetic assemblages inferred from Brazilian R. (B.) microplus are more genetically diverse than those reported from other regions (Figure 3). This study sought to obtain high support for the topology of tree and the COX-I molecular marker was found to have the advantage of presenting greater mutation rates than the nuclear marker, corroborating the data presented by Lv et al. (2014). However, the ITS-2 molecular marker failed to resolve discrepancies between the populations of R. (B.) microplus and this study corroborates the findings of Burger et al. (2014) in which they showed that the ITS-2 molecular marker has insufficient power to distinguish between very close species. The comparative analysis in this study shows that R. (B.) microplus ticks represent at least two different populations. Moreover, this study also attempted to investigate whether haplotypes are Figure 3. Bayesian phylogenetic tree recovered for cytochrome oxidase subunit 1 gene (COX-I) of 643 nucleotides, for Rhipicephalus taxa samples. In blue and green are the smallest nodes formed the Brazilian populations of R. (B.) microplus. Dermacentor nitens (GenBank: KC503258) was used as the outgroup. 324 Csordas, B.G. et al. Braz. J. Vet. Parasitol. Figure 4. Bayesian Cladogram recovered for internal transcript spacer 2 (ITS-2) of 580 nucleotides, for sample Rhipicephalus taxa. Dermacentor nitens KC503275 was used as the outgroup. distributed according to the phylogeographical areas. It is unknown whether this indicates a cause-effect association of haplotypes and geographical areas. However, the prevalence of ticks has been typically associated with cattle breeds, the highly invasive and widespread movement of R. (B.) microplus may facilitate its occurrence (IBELLI et al., 2012; BUSCH et al., 2014). Conclusion This study provides new insights into the distinct genetic assemblages of the tick R. (B.) microplus in Brazil and reveals that the Brazilian R. (B.) microplus population consists of at least two different populations. In addition, future studies using a large number of samples from the Americas using COX-1 may v. 25, n. 3, jul.-set. 2016 Molecular characterization of Rhipicephalus (Boophilus) microplus in Brazil be able to show how these populations differ in their inter- and intra- relationships. 325 Acknowledgements Estrada-Peña A, Bouattour A, Camicas JL, Guglielmone A, Horak I, Jongejan F, et al. The known distribution and ecological preferences of the tick subgenus Boophilus (Acari: Ixodidae) in Africa and Latin America. Exp Appl Acarol 2006; 38(2-3): 219-235. http://dx.doi.org/10.1007/ s10493-006-0003-5. PMid:16596355. We thank FUNDECT - MS, the State of Mato Grosso do Sul, CAPES and CNPq for financial support. Ewing B, Green P. Base-calling of automated sequencer traces using phred. II. Error probabilities. Genome Res 1998; 8(3): 186-194. http:// dx.doi.org/10.1101/gr.8.3.186. PMid:9521922. References Ewing B, Hillier L, Wendl MC, Green P. Base-calling of automated sequencer traces using phred. I. Accuracy assessment. Genome Res 1998; 8(3): 175-185. http://dx.doi.org/10.1101/gr.8.3.175. PMid:9521921. Altschul SF, Madden TL, Schäffer AA, Zhang J, Zhang Z, Miller W, et al. Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. Nucleic Acids Res 1997; 25(17): 3389-3402. http:// dx.doi.org/10.1093/nar/25.17.3389. PMid:9254694. Excoffier L, Lischer HEL. Arlequin suite ver 3.5: a new series of programs to perform population genetics analyses under Linux and Windows. Mol Ecol Resour 2010; 10(3): 564-567. http://dx.doi.org/10.1111/j.17550998.2010.02847.x. PMid:21565059. Andreotti R, Guerrero FD, Soares MA, Barros JC, Miller RJ, Léon AP. Acaricide resistance of Rhipicephalus (Boophilus) microplus in State of Mato Grosso do Sul, Brazil. Rev Bras Parasitol Vet 2011; 20(2): 127-133. http:// dx.doi.org/10.1590/S1984-29612011000200007. PMid:21722487. Grisi L, Leite RC, Martins JR, Barros AT, Andreotti R, Cançado PH, et al. Reassessment of the potential economic impact of cattle parasites in Brazil. Rev Bras Parasitol Vet 2014; 23(2): 150-156. http://dx.doi.org/10.1590/ S1984-29612014042. PMid:25054492. Bandelt HJ, Forster P, Röhl A. Median-joining networks for inferring intraspecific phylogenies. Mol Biol Evol 1999; 16(1): 37-48. http:// dx.doi.org/10.1093/oxfordjournals.molbev.a026036. PMid:10331250. Guerrero FD, Andreotti R, Bendele KG, Cunha RC, Miller RJ, Yeater K, et al. Rhipicephalus (Boophilus) microplus aquaporin as an effective vaccine antigen to protect against cattle tick infestations. Parasit Vectors 2014; 7: 475. PMid:25306139. Barré N, Uilenberg G. Spread of parasites transported with their hosts: case study of two species of cattle tick. Rev Sci Tech 2010; 29(1): 149-160, 135-147. PMid:20617654. Hebert PDN, Cywinska A, Ball SL, deWaard JR. Biological identifications through DNA barcodes. Proc Biol Sci 2003; 270(1512): 313-321. Barros-Battesti DM, Arzua M, Bechara GH. Carrapatos de importância médico-veterinária da região neotropical: Um guia ilustrado para identificação de espécies. São Paulo: Instituto Butantan, 2006. 223 p. Barrett RDH, Hebert PDN. Identifying spiders through DNA barcodes. Can J Zool 2005; 83(3): 481-491. http://dx.doi.org/10.1139/z05-024. Busch JD, Stone NE, Nottingham R, Araya-Anchetta A, Lewis J, Hochhalter C, et al. Widespread movement of invasive cattle fever ticks (Rhipicephalus microplus) in southern Texas leads to shared local infestations on cattle and deer. Parasit Vectors 2014; 7(1): 188. http:// dx.doi.org/10.1186/1756-3305-7-188. PMid:24742041. Burger TD, Shao R, Barker SC. Phylogenetic analysis of mitochondrial genome sequences indicates that the cattle tick, Rhipicephalus (Boophilus) microplus, contains a cryptic species. Mol Phylogenet Evol 2014; 76: 241253. http://dx.doi.org/10.1016/j.ympev.2014.03.017. PMid:24685498. Campbell BC, Steffen-Campbell JD, Werren JH. Phylogeny of the Nasonia species complex (Hymenoptera: Pteromalidae) inferred from an internal transcribed spacer (ITS2) and 28S rDNA sequences. Insect Mol Biol 1993; 2(4): 225-237. http://dx.doi.org/10.1111/j.1365-2583.1994. tb00142.x. PMid:9087560. Cruickshank RH. Molecular markers for the phylogenetics of mites and ticks. Syst Appl Acarol 2002; 7(1): 3-14. http://dx.doi.org/10.11158/ saa.7.1.1. Dantas-Torres F. Climate change, biodiversity, ticks and tick-borne diseases: the butterfly effect. Int J Parasitol Parasites Wildl 2015; 4(3): 452-461. http://dx.doi.org/10.1016/j.ijppaw.2015.07.001. PMid:26835253. 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. Front Biosci 2008; 13(13): 6938-6946. http://dx.doi. org/10.2741/3200. PMid:18508706. Hebert PDN, Stoeckle MY, Zemlak TS, Francis CM. Identification of birds through DNA barcodes. PLoS Biol 2004; 2(10): e312. http://dx.doi. org/10.1371/journal.pbio.0020312. PMid:15455034. Hoogstral H. Theobald Smith: his scientific work and impact. Bull Entomol Soc Am 1986; 32(1): 22-35. http://dx.doi.org/10.1093/besa/32.1.22. Huelsenbeck JP, Ronquist F. MRBAYES: bayesian inference of phylogenetic trees. Bioinformatics 2001; 17(8): 754-755. http://dx.doi.org/10.1093/ bioinformatics/17.8.754. PMid:11524383. Ibelli AM, Ribeiro AR, Giglioti R, Regitano LC, Alencar MM, Chagas ACS, et al. Resistance of cattle of various genetic groups to the tick Rhipicephalus microplus and the relationship with coat traits. Vet Parasitol 2012; 186(3-4): 425-430. http://dx.doi.org/10.1016/j.vetpar.2011.11.019. PMid:22115946. Kanduma EG, Mwacharo JM, Sunter JD, Nzuki I, Mwaura S, Kinyanjui PW, et al. Micro-and minisatellite-expressed sequence tag (EST) markers discriminate between populations of Rhipicephalus appendiculatus. Ticks Tick Borne Dis 2012; 3(3): 128-136. http://dx.doi.org/10.1016/j. ttbdis.2012.05.001. PMid:22789728. Labruna MB, Naranjo V, Mangold AJ, Thompson C, Estrada-Peña A, Guglielmone AA, et al. Allopatric speciation in ticks: genetic and reproductive divergence between geographic strains of Rhipicephalus (Boophilus) microplus. BMC Evol Biol 2009; 9(1): 46. http://dx.doi. org/10.1186/1471-2148-9-46. PMid:19243585. Librado P, Rozas J. DnaSP v5: a software for comprehensive analysis of DNA polymorphism data. Bioinformatics 2009; 25(11): 1451-1452. http://dx.doi.org/10.1093/bioinformatics/btp187. PMid:19346325. Low VL, Tay ST, Kho KL, Koh FX, Tan TK, Lim YAL, et al. Molecular characterisation of the tick Rhipicephalus microplus in Malaysia: new insights into the cryptic diversity and distinct genetic assemblages throughout the world. Parasit Vectors 2015; 8(1): 341. http://dx.doi.org/10.1186/ s13071-015-0956-5. PMid:26104478. 326 Csordas, B.G. et al. Lv J, Wu S, Zhang Y, Chen Y, Feng C, Yuan X, et al. Assessment of four DNA fragments (COI, 16S rDNA, ITS2, 12S rDNA) for species identification of the Ixodida (Acari: Ixodida). Parasit Vectors 2014; 7(1): 93. http://dx.doi.org/10.1186/1756-3305-7-93. PMid:24589289. Peter RJ, Van den Bossche P, Penzhorn BL, Sharp B. Tick, fly, and mosquito Control-Lessons from the past, solutions for the future. Vet Parasitol 2005; 132(3-4): 205-215. http://dx.doi.org/10.1016/j. vetpar.2005.07.004. PMid:16099104. Rodríguez-Vivas RI, Rivas AL, Chowell G, Fragoso SH, Rosario CR, García Z, et al. Spatial distribution of acaricide profiles (Boophilus microplus strains susceptible or resistant to acaricides) in southeastern Mexico. Vet Parasitol 2007; 146(1-2): 158-169. http://dx.doi.org/10.1016/j. vetpar.2007.01.016. PMid:17349747. Braz. J. Vet. Parasitol. Ronquist F, Huelsenbeck JP. MRBAYES 3: bayesian phylogenetic inference under mixed models. Bioinformatics 2003; 19(12): 1572-1574. http:// dx.doi.org/10.1093/bioinformatics/btg180. PMid:12912839. Song S, Shao R, Atwell R, Barker S, Vankan D. Phylogenetic and phylogeographic relationships in Ixodes holocyclus and Ixodes cornuatus (Acari: Ixodidae) inferred from COX1 and ITS2 sequences. Int J Parasitol 2011; 41(8): 871-880. http://dx.doi.org/10.1016/j.ijpara.2011.03.008. PMid:21540032. Tree Bio. FigTree [online]. London; 2016 [cited 2016 Apr 11]. Available from: http://tree.bio.ed.ac.uk/ Wright S. Evolution in mendelian populations. Genetics 1931; 16(2): 97-159. PMid:17246615.