- Research article
- Open Access
Amblyomma cajennense (Fabricius, 1787) (Acari: Ixodidae), the Cayenne tick: phylogeography and evidence for allopatric speciation
BMC Evolutionary Biology volume 13, Article number: 267 (2013)
Amblyomma cajennense F. is one of the best known and studied ticks in the New World because of its very wide distribution, its economical importance as pest of domestic ungulates, and its association with a variety of animal and human pathogens. Recent observations, however, have challenged the taxonomic status of this tick and indicated that intraspecific cryptic speciation might be occurring. In the present study, we investigate the evolutionary and demographic history of this tick and examine its genetic structure based on the analyses of three mitochondrial (12SrDNA, d-loop, and COII) and one nuclear (ITS2) genes. Because A. cajennense is characterized by a typical trans-Amazonian distribution, lineage divergence dating is also performed to establish whether genetic diversity can be linked to dated vicariant events which shaped the topology of the Neotropics.
Total evidence analyses of the concatenated mtDNA and nuclear + mtDNA datasets resulted in well-resolved and fully congruent reconstructions of the relationships within A. cajennense. The phylogenetic analyses consistently found A. cajennense to be monophyletic and to be separated into six genetic units defined by mutually exclusive haplotype compositions and habitat associations. Also, genetic divergence values showed that these lineages are as distinct from each other as recognized separate species of the same genus. The six clades are deeply split and node dating indicates that they started diverging in the middle-late Miocene.
Behavioral differences and the results of laboratory cross-breeding experiments had already indicated that A. cajennense might be a complex of distinct taxonomic units. The combined and congruent mitochondrial and nuclear genetic evidence from this study reveals that A. cajennense is an assembly of six distinct species which have evolved separately from each other since at least 13.2 million years ago (Mya) in the earliest and 3.3 Mya in the latest lineages. The temporal and spatial diversification modes of the six lineages overlap the phylogeographical history of other organisms with similar extant trans-Amazonian distributions and are consistent with the present prevailing hypothesis that Neotropical diversity often finds its origins in the Miocene, after the Andean uplift changed the topology and consequently the climate and ecology of the Neotropics.
Amblyomma cajennense Fabricius (Figure 1) is one of the most widely distributed tick species in the New World. Its range extends from northern Argentina, to the Caribbean and the southernmost part of the U.S. (from 27° N to 29° S). Throughout its distribution, this tick has adapted to widely different ecological conditions, including ecosystems as different as semi-arid grasslands and subtropical secondary forests . The geographical area occupied by this tick is interspersed with major geographical barriers: the Andes, the Gulf of Mexico, and large rivers [1, 2] (Figure 2A).
Early morphological studies resulted in contradicting opinions about the taxonomy of A. cajennense. Some authors considered phenotypic differences (presence and number of ventral plates, proportional length of festoons, ornamentation, scutal size and shape, punctation, and shape of porose areas) to identify different species within A. cajennense[3–5], whereas other authors ascribed morphological differences to mere intraspecific polymorphism [6–9]. The latter point of view prevailed and, since 1953 , all related species were relegated to synonyms of A. cajennense. Lack of variation in the cuticular hydrocarbons reported from twelve geographically distinct populations of A. cajennense supported the synonymy .
Nevertheless, after observing different distinct developmental features (lengths of feeding, pre-molting, pre-ovoposition, and egg incubation periods) in laboratory colonies of A. cajennense from different Neotropical areas (Argentina, Brazil, Cuba, U.S., and Trinidad), Guglielmone et al.  suggested that A. cajennense might be a complex of sibling species. More recently, biological crossbreeding incompatibility was observed between colonies of A. cajennense from different geographical areas [12, 13]. Unsuccessful crossmating experiments constitute a strong indication of the occurrence of different species. Nevertheless, laboratory breeding conditions may not fully represent natural settings and may artificially impede mating processes that would normally occur in nature.
In order to further investigate these apparently contradicting findings it was therefore necessary to use molecular methods as additional tools for developing better informed taxonomic opinions. Molecular sequence analyses were also used to evaluate the phylogeographical evolutionary history of the taxon as its present distribution, similar to that of other Neotropical organisms [14–22], can contribute to better understand the respective influence of Pleistocene versus Miocene biogeographical events in generating diversity in the Neotropics.
The present distribution of taxa and their diversity is the result of past demographic events such as colonization, expansion, and extinction, usually driven by habitat modifications. The evolution and population differentiation in ticks, obligate blood-sucking ectoparasites of vertebrates, are shaped by two main driving forces: biogeographical vicariance and host association [23–26].
In terms of vicariance, in the Neotropical area, the richness in species diversity had initially been attributed to habitat fragmentation due to the fluctuations in levels of dryness  or temperature  during the Pleistocene. By degrees, the scientific community realized that the diversification of organisms in the Neotropics could not be solely be attributed to the effect of Quaternary glaciation on Neotropical climate [29–34], but was likely to have started earlier under the effect of the Andean orogenesis . More recently, particular attention has been given to vicariant trans-Amazonian taxa distribution. The geographical range of such organisms often overlaps the remaining nuclei of the Neotropical Seasonally Dry Tropical Forest (SDTF) and of the adjoining dry savannas of the Chaco and Cerrado regions, which are thought to have been widespread during the Neogene [15, 29–31, 36–38]. The distribution of the Cayenne tick  roughly coincides with a typical trans-Amazonian distribution and significantly overlaps SDTF, Chaco and Cerrado refugial foci (Figure 2B) [29, 39]. Ticks collected in deeper Amazonian areas (i. e. Rondonia) usually inhabit drier corridors and are found at the interface between forest and savanna and not in the thick of the rain forest . One can, therefore, hypothesize, that like many other organisms with a trans-Amazonian distribution, the Cayenne tick has undergone allopatric speciation in the different ecological regions which resulted from the fragmentation of SDTF and of environments drier and more markedly seasonal than rain forest habitats. However, if we consider host association as the main driving mechanism underlying tick diversification, there is no doubt that A. cajennense is very eclectic in its feeding preferences. Adult specimens mainly infest large wild and domestic mammals, ungulates in particular ; immature stages are less specific in their host choice and can feed on a broader range of vertebrates, including marsupials, small and large mammals, humans, and sometimes birds [42–45]. Therefore, this tick not only can find suitable hosts practically everywhere, but it can easily be carried over large distances on wild ungulates, cattle displaced by humans, or on birds. Hillburn and Sattler  postulated that populations of livestock ticks, when not host-specific and in the presence of abundant host fauna, should closely approach panmixia.
Hence, two opposite evolutionary hypotheses can let us expect either that isolation mechanisms have resulted in the subdivision of A. cajennense into genetically distinct populations, or that persistent gene flow, maintained by dispersal on hosts, has created a genetic gradient of progressive differentiation throughout the tick distribution range. In the present study, we investigated the phylogeny and population structure of A. cajennense throughout its distribution range by analyzing three mitochondrial and one nuclear gene sequences. Data were used to determine the extent of genetic structure within the taxon and its compatibility with the hypothesized occurrence of distinct cryptic species, to evaluate the demographic and phylogeographical history of the tick, and to tentatively date lineage radiation events.
Our samples included wild caught ticks and ticks from colonies (three Brazilian and one Colombian) (Table 1). The Brazilian colonies originated from São Paulo, Rio de Janeiro, and Rondonia and the Colombian from Cundinamarca. As colony ticks can evolve and differentiate very rapidly from wild caught specimens , all the colony ticks used in this study were first generation specimens that had not yet experienced the effects of inbreeding. Fresh wild samples were collected in Mexico, Argentina, Costa Rica, Perú, French Guiana, Colombia, and Ecuador. Additional ticks were selected from the holdings in the U.S. National Tick Collection (USNTC, Georgia Southern University), the Collection of the Instituto Nacional de Tecnología Agropecuaria (INTA, Rafaela, Argentina), the Coleção Nacional de Carrapatos (CNC, Faculty of Veterinary Medicine and Zootechnoloy, Univeristy of São Paulo, Brazil), and the Coleção Acarológica do Instituto Butantan (IBSP, Butantan Institute, São Paulo, Brazil). Sequences of Amblyomma imitator Kohls, 1958, a North American species often confused with A. cajennense, and Amblyomma americanum L. were used for outgroup comparisons.
DNA extraction, PCR, cloning, and sequencing
DNA was extracted and the cuticles of the ticks were preserved for further morphological analysis following previously published protocols [48, 49]. This involved cutting off a small portion of the postero-lateral idiosoma of each tick with a disposable scalpel, before an overnight incubation in 180 μl Qiagen ATL lysis buffer (Qiagen, Valencia, CA) and 40 μl of 14.3 mg/ml proteinase K (Roche Applied Sciences, Indianapolis, IN). After complete lysis of the tick tissues and repeated vortexing, the cuticle was stored in 70% ethanol and kept as a voucher specimen. The lysed tissues were further processed as previously described [48, 49]. Two mitochondrial gene sequences, 12SrDNA and the control region or d-loop (DL) were amplified as previously reported [48, 49]. Primers F2LITS2 = 5’ - tgagggtcggatcayatatca -3’ and McLn = 5’ - gtgaattctatgcttaaattcagggggt - 3’  were used to amplify approx. 950 bp long fragments of the Internal Transcribed Spacer 2 (ITS2) . Primers for the COII gene fragment, COIIF (5’-tca gaa cay wcy tty aat caa aat -3’) and COIIR2 (5’-cca caa att tct gaa cat tgw cca-3’), were selected within the cytochrome c oxidase subunit I and at the end of the COII sequence respectively, by comparing the complete mitochondrial genomes of Rhipicephlaus sanguineus (GenBank accession number: NC002074.1) and Ixodes hexagonus (GenBank accession number: NC002010.1). A touch-down annealing sequence from 55 to 46°C was used to amplfy COII sequences. PCR was performed by using a MasterTaq kit (5-Prime, Gaithersburg, MD). Each reaction contained 2.5 μl of tick DNA, 2.5 μl of 10 × Taq buffer, 5 μl of 5 × TaqMaster PCR Enhancer, 1.5 μl of MgAc (25 mM), 0.5 μl dNTP mix (10 mM each), 0.1 μl of Taq polymerase (5U/ μl), 1.25 μl of each primer from a 10 pmoles/ μl stock solution (Invitrogen, Life Technologies Corporation, Grand Island, NY), and 14.6 μl molecular biology grade H2O. The two DNA strands of each amplicon were purified and sequenced at the High-Throughput Genomics Unit (HTGU, University of Washington, Seattle, WA) and were assembled with Sequencer 4.5 (Gene Codes Corporation, Ann Arbor, MI). Because the ITS2 region occurs in multiple copies in the same organism, it was necessary to compare its variability within single tick specimens, within populations, and between populations to ensure that, in A. cajennense, its evolution reflected concerted evolution. For a subset of geographical regions, 2 specimens were randomly selected and the corresponding ITS2 PCR product cloned using TOPO-TA (Invitrogen). For each cloning, eight colonies containing the tick DNA insert were grown overnight in LB medium supplemented with 10% glycerol and ampicillin (Teknova, Hollister, CA), frozen, and sent for sequencing to the HTGU.
Genetic structure and demography
Relationships between sequences were first investigated by generating unrooted networks with the statistical parsimony method implemented in TCS 1.13  and a confidence interval of 95%. Haplotypes were considered to be distinct when they differed by at least one bp. Individual sequence alignments were imported into DnaSP (DNA Sequence Polymorphism)  version 5.10.01 for analysis of nucleotide polymorphism. The program was used to calculate haplotype and nucleotide diversity (Hd ±S D and π±S D) and to perform Ramos-Onsin and Roza’s R2  and Fu’s Fs [54, 55] neutrality tests, considered to be the most robust methods for studying the effect of demographical events on DNA sequence data . For each dataset, genetic differentiation (F ST ) between the clades identified through TCS and phylogenetic analyses were estimated by also using DnaSP. The statistical significances of all tests performed with DnaSP were estimated by the coalescent method with 95% confidence interval and 10,000 permutations. The significance of fixation indices was determined by comparison to a null distribution of these values based on 10,000 random permutations .
Sequences were manually aligned with McClade 4.07 OSX . Secondary structure was considered in aligning 12SrDNA  and DL sequences, and codon organization in aligning the COII data set. The homogeneity of base frequencies across our sample was evaluated with a χ2 goodness-of-fit test using PAUP 4.0b10  prior to all phylogenetic analyses. Substitution saturation was evaluated with DAMBE 18.104.22.168 [59–61]. Each data set was analyzed by maximum parsimony (MP) and maximum likelihood (ML) with PAUP, and Bayesian analysis (BA) with MrBayes 2.01 . Branch support was assessed by bootstrap analysis (1000 replica) with PAUP for MP, with PHYML (100 replica)  in Phylogeny.fr  for ML, and by posterior probability with MrBayes. MP heuristic searches were performed by branch-swapping using the tree bisection-reconnection (TBR) algorithm, ACCTRAN character optimization, with all substitutions given equal weight, and with 10 random sequence addition replicates. Gaps were treated either as a 5th (in 12srDNA and D-loop analyses) or as a missing character. Maximum likelihood heuristic searches were run after the nucleotide substitution model best fitting the data was selected by Modeltest v3.7 . The MP tree with the best ML score was used as the starting tree for ML searches. Two runs, with four chains each, were run simultaneously for BA analyses (1,000,000 generations). Trees were sampled every 100 iteration. Trees saved before the average standard deviation of split fragments converged to a value < 0.01 were discarded from the final sample. When necessary, the number of generations was increased so that the number of discarded samples would not exceed 25% of the total sampled trees. The 50% majority-rule consensus tree of the remaining trees was inferred and posterior probabilities recorded for each branch. The four data matrices were compared for congruence by using the partition-homogeneity test, with 100 replicates and significance threshold value P < 0.05 as implemented in PAUP . Congruent data sets were combined for total evidence analyses. Two concatenated data sets, one including only mitochondrial (mtDNA) and one including both, mitochondrial a nuclear sequences (n+mtDNA), were analyzed following the same procedure outlined for the separate analyses. BA concatenated analyses were partitioned by gene and codon position (for the COII portion of the dataset).
Molecular clock and divergence dates
In order to test substitution rate variation among lineages, for each dataset (four separate genes, mtDNA, and n+mtDNA), relative-rate tests were applied to all sister clades by using DAMBE . In addition, the molecular clock hypothesis was tested by likelihood ratio-test also in DAMBE . Tentative estimates of divergence time were performed by using the n+mtDNA concatenated data matrix with the relaxed molecular clock model implemented in BEAST v.1.7.4 [67, 68]. BEAST analyses were run analyzing the four-gene set, but allowing for independent model parameters for each partition. The ingroup was defined as being monophyletic in agreement with the phylogenetic analysis results. The tree priors was set to Yule Process, and the molecular clock set to uncorrelated lognormal distribution. Chain lengths were set to 1,500,000 and data were sampled every 750 iterations with a random starting tree. Two independent runs were combined with LogCombiner v.1.7.4 in order to reach stable posterior distributions in Tracer v.1.5. After deleting 10% of the generated trees, the remaining trees were summarized as a combined maximum clade credibility tree by using TreeAnnotator v.1.7.4. FigTree v.1.3.1 was used to visualize tree structure, with mean divergence times. Because dated fossil records are missing for this group of ticks, node dating was based on known vicariant events which shaped the territory occupied by A. cajennense. Different and independent node dating strategies were applied in BEAST analysis and the results compared. In order to account for the uncertainty of these dates, normally distributed tree priors were used. Assuming a vicariance model, dating was first attempted by calibrating single nodes which is allowed when applying relaxed molecular clocks. As the basal lineages always included the PE branch, which is located in a SDTF hot spot supposedly isolated for at least 10 Mya in the Inter-Andean Valley, we surmised that the differentiation within A. cajennense initiated before that, most likely during a period coinciding with the beginning of the Andean uplift. Therefore, the node between out- and ingroup was set at 20 Mya ±5 Mya; the PE - AR node was set to a median of 10 Mya ±2 Mya; the origin of the NE-CO-NW clade was set to a median of 6 Mya ±1 Mya which corresponds to the end of the uplift of the eastern Cordillera in its northernmost portion; the origin of the NW clades was set at 3 ± 0.5 Mya which would represent the closure of the Isthmus of Panama. Next, we applied two different double node calibrations: the first with outgroup - ingroup at 20 Mya ±5 Mya and PE-AR at 10 Mya ±2 Mya, the second with PE-AR at 10 Mya ±2 Mya and NW-NE-CO at 6 Mya ±1 Mya.
Results and discussion
The tick sample (Table 1) contained specimens from most of the New World biomes  inhabited by A. cajennense. The localities included areas with temperate grasslands, tropical and subtropical grasslands, tropical and subtropical moist broadleaf forests, tropical and subtropical dry broadleaf forests, and desert and xeric shrubland. Our sampling extended the known geographical distribution of the species to the montane shrubland of the inter-Andean Valley of Perú and the coastal mixed forest of Ecuador [1, 2]. Overall, the analysis included specimens from 19 localities and 11 countries (Table 1). Collection sites or colony origins are represented in Figure 2A. For clarity, following abbreviations are used throughout the paper and in the illustrations: TX = Texas, MX = Mexico, CR = Costa Rica, EC = Ecuador, VE = Venezuela, FG = French Guiana, CHO = Chaco Occidental (Argentina), CHS = Chaco Serrano (Argentina), YU = Yungas (Argentina), AF = Atlantic Forest (Brazil), RO = Rondonia (Brazil), MG = Mato Grosso (Brazil), PA = Paraguay, PE = Perú, CO = Colombia. The phylogenetic clades were named after the general geographical range they occupy within the distribution of A. cajennense: NW = north-western clade (Texas, Mexico, Costa Rica, Ecuador); CO = Colombia; NE = north-eastern clade (French Guiana, Brazil-Rondonia); PE = Perú (Inter-Andean valley); EA = eastern clade (Brazil - Atlantic Forest - Mato Grosso and Argentina - Yungas); AR = Argentine clade (Argentine and Paraguayan Chaco).
DNA extraction, PCR, cloning, and sequencing
The number of sequences obtained for each geographical area was variable due to the fact that some tick samples had been freshly collected, whereas older ethanol-preserved samples from the USNTC sometimes yielded little or no DNA at all. The amplification success also varied depending on primer sets, with 12SrDNA being the easiest gene fragment to amplify. Haplotype distributions for each gene are listed in Additional files 1, 2, 3 and 4. GenBank accession numbers are listed in the Availability of supporting data section.
Sequence diversity and TCS analyses
The 123 12SrDNA sequences were represented by 33 distinct haplotypes, the 60 COII sequences by 26 haplotypes, and the 110 DL sequences by 31 haplotypes (Additional files 1, 2 and 3). The initial 444 bp long DL alignment needed to be reduced to 394 characters, because it included a hypervariable region alignable within clades, but not between clades. This portion of the alignment was, therefore, eliminated from further data analysis. The TCS analyses (95% parsimony cut-off level) separated mitochondrial haplotypes in distinct networks that could not be joined even after reducing the cut-off level to 70–80% (Additional files 5, 6 and 7). Independently on the gene analyzed, the geographical composition of the subnetworks (labeled NW, CO, NE, EA, AR, and PE) was the same. NW consistently included samples from TX, MX, CR and EC; NE samples from FG and RO; CO samples from Colombia; EA samples from YU, MG, AF; AR samples from CHO, CHS, and PA; and PE samples from the inter-Andean Valley of Perú. The network subdivision indicated that mitochondrial lineages of A. cajennense were deeply split and separated by a large number of parsimonious steps. It also showed that the three mitochondrial markers were congruent in splitting A. cajennense in the exact same subunits. In general, haplotypes were mostly unique to their area of origin, with the exception of haplotypes consistently shared by EC, MX, and CR, and by AF, YU, and MG. In addition, when available, the PA haplotypes were identical or at least very similar (1 bp difference) to the CHO, rather than to the CHS haplotypes, which is geographically easily explained. Also, when available, the VE and FG haplotypes clustered together. The 78 ITS2 gene sequences were collapsed into 20 distinct genotypes (Additional file 4). The ITS2 sequences were more conserved than the mitochondrial gene sequences with, for instance, French Guiana and Rondonia sharing identical genotypes. Similarly, although, separated by long branches (10 steps), CO and NE were included in the same network as were NW and EA (Additional file 8). By increasing the cut-off level to 97% (approx. 24 bp/825 total bp), TCS resulted in the exact same split recorded for the mitochondrial genes. The individual networks were structured, often including long branches and loops, and lacked the star patterns identifying recent sudden radiations. The populations they represented, therefore, appeared to have diversified over time. These findings indicated that the analysis of both mitochondrial and nuclear markers revealed congruent marked divergence between samples from different areas. The lack of connection between subnetworks did not allow us to determine the possible geographical pathways followed by A. cajennense throughout its radiation around the Amazon basin. The ITS2 network, however, indicated that the NW ticks evolved from the EA cluster, and that NE and CO were more closely related to each other than to the other groups. Measures of ML sequence divergences (Tables 2 and 3) showed that diversity within clades was very low varying from 0 to 4% (reaching 4% only with COII, the most diverse gene). Values between clades, between outgroups and ingroup, and between outgroups were comparable and typically about 8–10 times higher than intra-clade values, with a single exception: CO and NE only differed by 1.15% in ITS2, which is significantly lower than the other observed inter-clade values. This was the only significant discrepancy observed between mtDNA and ITS2 data, indicating that the separation between CO and NE may not be as clearly defined as it is between other clades when using a nuclear recombining rather than a mitochondrial marker. In general, however, our data showed that the diversity between clades was similar or higher than that observed between the two very different outgroups or between the outgroups and the ingroup. This strongly suggested that the differences between the six clades identified through TCS and phylogenetic analyses (see below) were compatible with the occurrence of six separate species.
Tree properties and scores are shown in Table 4 and substitution models selected by the Akaike Information Criterion in Modeltest  for each alignment in Table 5. Phylogenetic reconstructions obtained with the separate datasets are shown in Additional files 5, 6, 7 and 8, whereas the mtDNA and the n+mtDNA trees are shown in Figures 3 and 4, respectively.
With very few exceptions, for each gene, the three phylogenetic methods (MP, ML, and BA) generated trees with overall similar topologies and lineage support. Although phylogenetic reconstructions based on the three mitochondrial genes were not topologically fully congruent (Additional files 5, 6 and 7), they consistently identified the same monophyletic units which corresponded to the subnetworks identified through TCS analysis. In general, 12SrDNA sequences appeared to provide better overall resolution than COII or DL, confirming the good level of information offered by small ribosomal subunit genes when investigating relationships among closely related taxa, but also within a single species [48, 49, 70, 71]. DAMBE did not reveal significant saturation in the 12SrDNA dataset. In the 12SrDNA reconstruction, PE was the basal well-supported lineage followed, in order of divergence, by the strongly supported AR and EA clades. The most recently evolving lineage grouped CO-NE-NW, which were each monophyletic. NE and CO were sister clades. For the COII dataset, saturation levels were tested for the whole data matrix, for each codon position separately, and for the two first codon positions together. None of these sets showed significant nucleotide saturation based on transition and transversion rates. Therefore, third codon positions were included in all further analyses. The topology of the COII tree differed from the 12SrDNA tree in few main aspects: PE-AR clustered in a monophyletic lineage, the lineages grouping the remaining clades were only supported in BA, and the position of CO was unresolved. The DL is one of the most variable regions in the mitochondrial genome . Although it’s usefulness for the study of the evolutionary history of arthropods has not yet been thoroughly investigated, this gene has proved to be informative at the intraspecific level for some arthropods [73, 74] and for two tick species, Ixodes ricinus and A. variegatum[49, 70]. In this study, the overall DL tree structure also lacked in support, particularly at the base of the otherwise well-supported CO-NE-EA-NW split. Although substitution saturation was initially suspected to be at the origin of the lack of resolution in DL and COII, all tests failed to find significant saturation in the mitochondrial markers.
Mitochondrial and nuclear genes do not always portray the same evolutionary history due to their distinct inheriting mechanisms. The non-recombining mitochondrial genes usually provide excellent phylogeographical information. Although they are increasingly being used for the delimitation of species (DNA barcoding), the identification of mitochondrial deep divergent splits may merely reflect past biogeographical events that do not always imply speciation. On that account, and in order to verify whether or not the analysis of a nuclear gene would result in the same clear cut subdivision of our samples, sequences of the noncoding rapidly evolving ITS2 regions were also sequenced. However, nuclear ribosomal DNA (rDNA) is known to occur in multiple copies in the genome. Different copies within the same individual usually evolve as a single-copy gene through concerted evolution [75, 76]. Nevertheless, exceptions to this rule have been observed in arthropods including ticks, particularly in the noncoding rDNA ITS regions which separate the transcribed genes [77–82]. If ribosomal DNA copies within a specimen are more diverse than between specimens, phylogenetic analyses using one representative sequence for each sampled specimen may not accurately represent relationships by descent. Phylogenetic and population genetics studies in ticks have, however, paradoxically often been based on analyses of ITS sequences [50, 83–88]. The ITS2 gene proved to be an informative marker at the intrageneric level in Neotropical Amblyomma species , and within populations of Amblyomma americanum, but not within Amblyomma variegatum, probably a more recently evolving species. After sequencing randomly chosen cloned ITS2 sequences from single specimens, the variability within (0.00–0.15%) was slightly lower than that between specimens from the same clade (0.00–0.88%). Moreover, when the intra-specimen substitutions were visually inspected, it became clear that they occurred as parsimony uninformative singletons and never involved informative segregating sites. Consequently, we concluded that ITS2 gene fragments were suitable markers for evolutionary studies in A. cajennense. The ITS2 tree topology (Additional file 8) was strongly supported at each hierarchical level and identified the same monophyletic clades revealed by the analysis of the mitochondrial genes. PE and AR were clustered in the basal monophyletic clade, a sister group to the well-supported EA-NW-NE-CO lineage. Within the latter, EA-NW, and NE-CO were sister branches.
The four data matrices proved to be phylogenetically informative with relatively little homoplasy (Table 4), with homogeneous base frequencies, and non significant substitution saturation. Partition homogeneity tests revealed that 12SrDNA-DL and COII-DL were congruent datasets, with p = 0.86 and p = 0.74 respectively, while 12SrDNA-COII were not (p=0.04). Because the latter significance level was low, and there was no significant conflict between DL and both other mitochondrial genes, the three datasets were concatenated. There was also no significant conflict’s signal between mtDNA and ITS2 (p=0.28) and, consequently, the four datasets were combined for total evidence analysis. The topology of the two concatenated trees (Figures 3 and 4) was identical, with strongly supported nodes (> 90%) at all levels, if we exclude the AR-PE lineage which was less supported in the MP mtDNA tree. The six main clades were all characterized by 100% bootstrap support. In order of divergence from basal to crown, the AR-PE clade was followed by the EA, the NW, and the NE-CO clades. The nodes giving rise to the six clades were set deep in the phylogenies, while the terminal branches were often polytomic (with the exception of the more diversified EA clade).
Historical demography and population genetics
We used diversity statistics and neutrality tests to look for signals of population growth in each of the clades (Table 6). Under a model of population expansion, because newer haplotypes are retained in the population and are recently derived, Hd is expected to be high whereas nucleotide diversity is expected to be low. It must be noted that, because indels are excluded from the DnaSP analyses, the number of haplotypes detected by DnaSP is lower than the number of haplotypes listed in Additional files 1, 2, 3 and 4. Although not all our samples included sufficient specimens or different haplotypes for robust analysis, in the total sample, observed haplotype diversity was significantly higher than expected (Table 6). Nucleotide diversity varied among mitochondrial genes, but was not significantly lower than expected in the mtDNA and the nuclear datasets within clades and in the total sample. This is an indication that the six clades constitute stable distinct populations and are not the result of recent population expansion or growth. Fu’s F and Ramos-Onsins R2 statistics performed on the nuclear and mtDNA datasets were not significant (Table 6), corroborating the idea that the six clades did not experience any sudden demographic event and that they have been evolving in stable conditions for a long period. The important level of genetic diversity observed in the clades may indicate that the ancestral taxon was genetically diverse before it became fragmented and that this genetic diversity was maintained in isolated allopatric refugia because they provided stable environmental conditions allowing for genetic accumulation. Relatively important genetic diversity can also be the result of limited dispersal and recent evidence  strongly corroborates the hypothesis  that many ticks are more strictly dependent on suitable environmental conditions than on the availability of specific hosts. Therefore, even if the tick theoretically could feed on a number of vertebrate hosts with high dispersal ability, it appears that they would not venture too far from the ecological niche they prefer. Whether some of this diversity is maintained through secondary contact along the suture zones between the areas occupied by the six clades will have to be determined through intensive sampling efforts and additional molecular analyses. In terms of population structure, F ST values between clades for each gene separately and for the concatenated datasets were very high (i.e. 0.91 – 0.99 for 12SrDNA, 0.92 – 1.00 for DL, 0.96 – 1.00 for COII, 0.80 – 1.00 for ITS2, and 0.95 – 0.99 for mtDNA) and were all highly significant with p values < 0.001 and < 0.01 confirming, once more, that the six lineages are genetically very distinct form each other.
Molecular clock and divergence dates
With all data sets, relative rate tests did not reveal significant differences between sister taxa at all evolutionary hierarchic levels within each tree. Rates were also not statistically different between ingroup and outgroup sequences, indicating that rate variation among lineages would not be the cause of misleading divergence date estimations. The molecular clock hypothesis was tested for each gene by the least-square method and the likelihood ratio test implemented in DAMBE . The molecular clock hypothesis could not be rejected in any of them with the exception of the LRT in ITS2 (p=0.03). As this was the only significant finding and its level of significance was low, we deduced that, overall, the datasets did not depart significantly form clockwise evolution. Divergence dates obtained by the different calibration criteria on the total evidence dataset (one node vs. two nodes) resulted in similar average radiation times, although the confidence intervals observed with the one-node calibrations were approximately 10% more important. The average dates obtained with the two two-node calibrations are shown in Figure 5 and the corresponding confidence intervals are listed in Table 7. If the confidence intervals were quite large at the base of the tree, particularly between outgroups and ingroup, they became, however, less prominent towards the crown bifurcations. Nevertheless, the average dates inferred with the two calibrating criteria were well within confidence intervals of both analyses. Therefore, independently on the calibration strategy applied to the analysis and within the hypothesized timeframe, the timing of the diverging events appeared to be reasonably consistent. Naturally, node dating based solely on what we believe are the biogeographical events involved in shaping the topology of the area occupied by A. cajennense may be misleading as the geological records are also estimates. Nevertheless, for lack of fossil records and based on these tentative evaluations, we can formulate some hypotheses about the temporal phylogeographical sequence of events based on an allopatric vicariant model of divergence (Figure 6). We can assume that the geographical distribution of the ancestor of A. cajennense was larger and covered the northern half of South-America during the first half of the Miocene, in environments that corresponded to the overall present ecological requirement of the tick . Reasons for PE and AR being isolated in the same clade are difficult to establish. One possible explanation is that the ancestors of the AR-PE and the ancestors of the EA-NW-NC-CO clades became isolated in the south-west and north-east, respectively, of the Oceanic introgression called Paranean Sea in the middle - late Miocene  (Figure 6). The PE clade, like other Andean lineages [91, 92], is separated by deep divergence from the related lowland taxa. We assume that PE became trapped between the progressively rising Andean Cordilleras because the area where our Peruvian samples were collected is thought to be a SDFT refugium isolated for at least 10-5 Mya since the rapid final phase of the Central Andean uplift [31, 93, 94]. After the retreat of the Paranean Sea from the end of the Miocene to the end of the Pliocene (10-3 Mya)  the AR and the EA clades may have reached each other again within a secondary contact zone corresponding to north-western Argentina. The NE-CO-NW clade split from the EA around 8-9 Mya when lacustrine ecosystems and swamps covered large parts of what would become Amazonia, creating an environment unsuitable for A. cajennense. Therefore, we can surmise that the tick populations were progressively being shifted to the periphery of the rain forest [35, 95, 96]. EA became established in the Atlantic Forest along the eastern coast of modern Brazil and in part of the Cerrado which is one of the elements of the so-called dry diagonal separating the Amazonian rain forest from the Atlantic Forest. A number of species are known to occur concomitantly in the Atlantic Forest and the Cerrado biome . The fact that EA ticks are also found in forests of the Yungas in north-western Argentina is not surprising as the Amazon forest and the Atlantic Forest were reportedly linked by continuous forest [97, 98]. The phylogeographical history of the NE-NW-CO clade is difficult to unravel. Nevertheless, we can speculate that it became separated from the EA clade in the late Miocene and that it adapted to ecosystems at the periphery of the Amazon basin along the western edge of the the pre-existing Cerrado-Caatinga formations. Diversification times for Cerrado woody lineages date back to 9.8 Mya . The divergence between NW and NE-CO might have happened after their common ancestor had established itself along the northern coast of South America. The origin of the NW, a mostly Central and North American clade with an extension along the Pacific coast of Ecuador, seems to predate the NE-CO split (approx. 6-8 Mya) and may coincide with the end of the orogenesis of the Central Northern Cordillera which could have confined this population along its western slopes. In addition, diversification could also have been driven by the presence of the Pebas aquatic system which isolated the north-western edge of South America from the Guiana shield  in the early-middle Miocene (24-11 Mya). The dating, however, would not coincide with the split between NW and NE-CO, which is more recent. From there, the ticks probably reached Central and North America after the closure of the Isthmus of Panama at the beginning of the Pliocene , although earlier dispersal on birds over water cannot be dismissed. The NE group includes samples collected in French Guiana and Rondonia and, in order to link these disjunct areas, the most likely hypothesis would suggest that the populations belonging to this clade occupied the northern coast of South American and, following the western side of the Cerrado-Caatinga corridor, exploited relatively drier environments within the rain forest. These ecosystems developed and have persisted since the early Miocene  as an entity distinct from the newly formed Amazon basin. The closest relative of the NE lineage is the CO and the most parsimonious explanation for its present location would involve vicariant separation from NE in the late Miocene, when the northern part of the Eastern Cordillera completed its orogenesis. More specifically, the eastern Cordillera started developing 25 Mya but only reached appreciable elevations in the late Miocene, as did the Merida Andes in Venezuela [103, 104]. The Sabana de Bogotà, where the Magdalena Valley is located and where the Colombian samples were collected, progressively increased in elevation between 15 and 3 Mya. The split between the NE and the CO lineages appears to coincide with the end of the rise of the Sabana de Bogotà. During approximately 3 to 9 Mya the main A. cajennense lineages did not appear to diversify indicating that they might have experienced a long period of genetic stability probably maintained by non-fluctuating environmental conditions, an hypothesis also confirmed by our demographic data. Nevertheless, we cannot exclude that extinction events during the late Miocene reduced the number of lineages to the extant ones. Diversification resumed in the late Pliocene (EA) and during the Pleistocene probably in response to the well documented climatic variability of the Quaternary age. Alternatively, as the basal lineages are PE and AR, we could imagine that the tick was well established in what corresponds now to its south-western range in the early Miocene and dispersed from there into acceptable ecological niches that developed around the newly formed Amazon Basin. The deep lineage split between clades would, however, rather suggest that the tick already occupied its present overall distribution area before the Amazon basin completed its development (6 Mya - ) and that its present allopatric distribution is the result of habitat fragmentation and not of a progressive northward dispersal.
Phylogenetic analyses and taxonomic considerations
Our data demonstrate that A. cajennense is subdivided into six genetically distinct groups. The question arises, therefore, as to whether or not these reciprocally monophyletic units represent different species. Regardless of the gene used, the six groups are characterized by unique fixed characters and mutually exclusive sequences, which is in agreement with the phylogenetic species definition . The level of genetic divergence between the clades is comparable to, or much higher than, genetic diversity recorded between different Amblyomma species represented here by the outgroup taxa. Significant very high among-clade differentiation estimates also confirm that there is significantly little genetic exchange between these entities. These results prompted a thorough reassessment of the morphology of A. cajennense, which resulted in the identification of previously ignored fixed phenotypic characters, differentiating the six groups . Therefore, if allopatric speciation between populations of A. cajennense has occurred, it is certainly not “cryptic”. While this study was completed, some of us [12, 13] carried out cross-breeding experiments with colonies established for some of the clades. The experiments showed that colonies of AR, EA, NE and CO were not compatible, at least in laboratory conditions. Therefore, morphological and biological studies further strengthen our results and collectively argue for the occurrence of six species, indicating that Koch  and Tonelli-Rondelli [4, 5] were mostly correct in their species delimitations. Formal descriptions of these species will be provided separately .
Phylogeographical patterns and tentative dating of the major lineages splits
The presence of Amblyomma fossils in Dominican amber deposits, dated from the late Oligocene-early Miocene (approx. 25 Mya), indicates that the genus was already well established in the New World at the end of the Oligocene. In present times, adults A. cajennense appear to favor ungulate hosts, such as horses, cows, deer, and pigs, but can also be found feeding on carnivores, marsupials and Edentata. Therefore, although the present preferred hosts invaded the Neotropical region from the Northern Hemisphere after the establishment of the Panama Isthmus, we can surmise that suitable hosts were available earlier. The geographical area occupied by A. cajennense significantly overlaps, or is closely adjacent to, the vestigial refugia of SDTF (Figure 2A-B) which are disjunct geographical areas with similar ecological conditions, identified mostly through the study of hundreds of plant lineages with coincidental distribution patterns [14, 15]. In addition, the tick distribution includes also areas of Chaco and Cerrado which were excluded by Pennington from the definition of SDTF  based on soil and vegetation compositions. Nevertheless Chaco and Cerrado are characterized by seasonality and long dry seasons like SDTF. Ticks, unlike plants, are more likely to depend on climatic conditions than on soil composition. Although the ecology of the ancestral SDTF was probably very similar, their long isolation from each other also coincided with gradual ecological modifications to which the ticks responded through increasing allopatric adaptation. For instance, while some clades may share similar coastal ecological conditions (EA, NE in French Guiana, and NW) there is no doubt that AR and PE occur in very different environments, the arid Chaco and the montane dry Inter-Andean valleys. The EA and AR ticks can be found 50 km apart, but they do not appear to venture out from their respective endemic areas. Nevertheless, a thorough exploration of the possible secondary contact zones between clades has yet to be undertaken and may shed some light on the real level of incompatibility in the natural environment between the identified species. Our attempt at node dating is in part speculative, as fossil records are largely missing for ticks, and because the dating of biogeographical topological modifications are also sometimes conjectural. Nevertheless, the phylogeography of A. cajennese is, sometimes partly but sometimes extensively congruent with the phylogeographical estimates generated for organisms with similar trans-Amazonian distributions, which include reptiles, rodents, and birds [16–20]. In some cases, divergence dating in these studies were supported by dated fossil records. Not only are the radiation patterns similar, but their dating often also matches our results, indicating that A. cajennense is only one of many organisms with an early Miocene origin and with affinity for trans-Amazonian regions with marked seasonality and long dry seasons [16–20]. Interestingly, the most extensive phylogeographical similarities are found in unrelated taxa, such as caviomorph rodents  and Crotalus spp. snakes . To the best of our knowledge, this is the first study dealing with an invertebrate with such a large and typical trans-Amazonian distribution range. The identification of these new species is not of solely taxonomic interest. It also has important implications for public health issues, because these ticks are vectors of important human pathogens. Accordingly, the subdivision revealed by our study also appears to correlate, in some cases, with distinct tick-pathogen associations [107–111].
Availability of supporting data
GenBank accession numbersfor 12SrDNA sequences are: EU791583-97, EU791599-600, EU791603-609, EU791611-615, AY342288, and JX987796-890. GenBank accession numbers for COII sequences are: KF787572-787631 and FJ860250-251; GenBank accession numbers for d-loop sequences are:KF527299-408. GenBank accession numbers for ITS2 sequences are: JN866835-JN866905, JN866908-910, and KF527286-298. Sequence alignments (concatenated mtDNA and mtDNA + ITS2) and the corresponding ML trees are accessible through TreeBase (URL: http://purl.org/phylo/treebase/phylows/study/TB2:S14894).
Estrada-Peña A, Guglielmone AA, Mangold AJ: The distribution and ecological ’preferences’ of the tick Amblyomma cajennense (Acari: Ixodidae), an ectoparasite of humans and other mammals in the Americas. Ann Trop Med Parasitol. 2004, 98 (3): 283-292.
Walker JB, Olwage A: The tick vectors of Cowdria ruminantium (Ixodoidea, Ixodidae, genus Amblyomma) and their distribution. Onderstepoort J Vet Res. 1987, 54 (3): 353-379.
Koch CL: Systematische Übersicht über die Ordnung der Zecken. Arch Naturgesch. 1844, 10: 217-239.
Tonelli-Rondelli M: Ixodoidea. Parte I. - Amblyomma ovale Koch, Amblyomma cajennense Fabricius e le specie a loro affini nuove o poco note. Riv Parassitol. 1937, 1: 273-299.
Tonelli-Rondelli M: Ixodoidea. Parte II. Contributo alla conoscenza della fauna ixodologica sud-americana. Spedizione del Prof. Beccari nella Gujana inglese ed a Trinidad. Riv Parassitol. 1939, 3: 39-55.
Neumann LG: Révision de la famille des Ixodidés (3ème mémoire). Mém Soc Zool France. 1899, 12: 107-294.
Neumann LG: Ixodidae. Das Tierreich, Volume 26 Lief. 1911, Berlin: R. Friedlander und Sohn, 169-169.
Aragão HB, da Fonseca F: Notas de ixodologia. V. A propósito da validade de algumas espécies do gênero Amblyomma do continente Americano (Acari: Ixodidae). Mem Inst Oswaldo Cruz (Rio de Janeiro). 1953, 51: 485-492.
Kohls GM: Amblyomma imitator, a new species of tick from Texas and Mexico, and remarks on the synonymy of A. cajennense (Fabricius) (Acarina-Ixodidae). J Parasitol. 1958, 44 (4): 430-433. Section 1
Estrada-Peña A, Guglielmone AA, Mangold AJ, Castellà J: Patterns of cuticular hydrocarbon variation and genetic similarity between natural populations of Amblyomma cajennense (Acari: Ixodoidea). Acta Trop. 1993, 55: 61-78.
Guglielmone AA, Mangold AJ, Oyola BC: Ciclo de vida de Amblyomma cajennense (Fabricius, 1787) (Acari: Ixodidae) en condiciones de laboratorio. Rev Med Vet (Buenos Aires). 1992, 73: 184-187.
Mastropaolo M, Nava S, Guglielmone AA, Mangold AJ: Biological differences between two allopatric populations of Amblyomma cajennense (Acari: Ixodidae) in Argentina. Exp Appl Acarol. 2011, 53 (4): 371-375.
Labruna MB, Soares JF, Martins TF, Soares HS, Cabrera RR: Cross-mating experiments with geographically different populations of Amblyomma cajennense (Acari: Ixodidae). Exp Appl Acarol. 2011, 54: 41-49.
Prado DE, Gibbs PE: Patterns of species distributions in the dry seasonal forests of South America. Ann Missouri Bot Gard. 1993, 80: 902-927.
Pennington RT, Prado DE, Pendry CA: Neotropical seasonally dry forests and Quaternary vegetation changes. J Biogeogr. 2000, 27: 261-273.
Solomon S, Bacci M, Martins J, Gonçalves-Vinha G, Mueller U: Paleodistributions and comparative molecular phylogeography of leafcutter ants (Atta spp.) provide new insight into the origins of Amazonian diversity. PLoS One. 2008, 3 (7): e2738-
Werneck FP, Gamble T, Colli GR, Rodrigues MT, Sites JW Jr: Deep diversification and long-term persistence in the South American “dry diagonal”: integrating continent-wide phylogeography and distribution modeling of geckos. Evolution. 2012, 66 (10): 3014-3034.
Martins FM, Templeton AR, Pavan ACO, Kohlbach BC, Morgante JS: Phylogeography of the common vampire bat (Desmodus rotundus): marked population structure, Neotropical Pleistocene vicariance and incongruence between nuclear and mtDNA markers. BMC Evol Biol. 2009, 9: 294-
Salazar-Bravo J, Dragoo JW, Tinnin DS, Yates TL: Phylogeny and evolution of the neotropical rodent genus Calomys: inferences from mitochondrial DNA sequence data. Mol Phylogenet Evol. 2001, 20 (2): 173-184.
Ribas CC, Miyaki CY, Cracraft J: Phylogenetic relationships, diversification and biogeography in Neotropical Brotogeris parakeets. J Biogeogr. 2009, 36 (9): 1712-1729.
Upham NS, Patterson BD: Diversification and biogeography of the Neotropical caviomorph lineage Octodontoidea (Rodentia: Hystricognathi). Mol Phylogenet Evol. 2012, 63 (2): 417-429.
Wüster W, Ferguson JE, Quijada-Mascareñas JA, Pook CE, Salomão MdG, Thorpe RS: Tracing an invasion: landbridges, refugia, and the phylogeography of the Neotropical rattlesnake (Serpentes: Viperidae: Crotalus durissus). Mol Ecol. 2005, 14 (4): 1095-1108.
Hoogstraal H, Aeschlimann A: Tick-host specificity. Mitt Schweiz Entomol Ges. 1982, 55: 5-32.
Klompen S, Black WC4th, Keirans JE, Oliver JH Jr: Evolution of ticks. Annu Rev Entomol. 1996, 41: 141-161.
McCoy KD, Boulinier T, Tirard C, Michalakis Y: Host-dependent genetic structure of parasite populations: differential dispersal of seabird tick host races. Evolution. 2003, 57 (2): 288-296.
McCoy KD, Chapuis E, Tirard C, Boulinier T, Michalakis Y, Bohec CL, Maho YL, Gauthier-Clerc M: Recurrent evolution of host-specialized races in a globally distributed parasite. P Roy Soc Lond B Bio. 2005, 272 (1579): 2389-2395.
Haffer J: Speciation in Amazonian Forest Birds. Science. 1969, 165: 131-137.
Colinvaux P: Pleistocene biogeography and diversity in tropical forests of South America. Biological Relationships Between Africa and South America. Edited by: Goldblatt P. 1993, New Haven: Yale University Press, 473-499.
Pennington RT, Lavin M, Prado DE, Pendry CA, Pell SK, Butterworth CA: Historical climate change and speciation: neotropical seasonally dry forest plants show patterns of both tertiary and quaternary diversification. Philos Trans R Soc Lond B Biol Sci. 2004, 359 (1443): 515-537.
Caetano S, Prado D, Pennington RT, Beck S, Oliveira-Filho A, Spichiger R, Naciri Y: The history of seasonally dry tropical forests in eastern South America: inferences from the genetic structure of the tree Astronium urundeuva (Anacardiaceae). Mol Ecol. 2008, 17 (13): 3147-3159.
Pennington RT, Lavin M, Särkinen T, Lewis GP, Klitgaard BB, Hughes CE: Contrasting plant diversification histories within the Andean biodiversity hotspot. Proc Natl Acad Sci USA. 2010, 107 (31): 13783-13787.
Werneck FP, Costa GC, Colli GR, Prado DE, Sites JW Jr: Revisiting the historical distribution of Seasonally Dry Tropical Forests: new insights based on paleodistribution modelling and palynological evidence. Global Ecol Biogeogr. 2011, 20: 272-288.
Patton JL, Costa P: Molecular phylogeography and species limits in rainforest didelphid marsupials of South America. Predators with Pouches: The Biology of Carnivorous Marsupials. Edited by: Jones M, Dickman C, Archer M. 2003, Collingwood: CSIRO Publishing, 63-81.
Antonelli A, Quijada-Mascareñas A, Crawford AJ, Bates JM, Velazco PM, Wüster W, Wesselingh FP: Molecular studies and phylogeography of Amazonian tetrapods and their relation to geological and climatic models. Amazonia, Landscape and Species Evolution: A Look into the Past. Edited by: Hoorn C. 2010, Hoboken, NJ: Wiley-Blackwell, 386-404.
Hoorn C, Wesselingh FP, ter Steege, Bermudez MA, Mora A, Sevink J, Sanmartín I, Sanchez-Meseguer A, Anderson CL, Figueiredo JP, Jaramillo C, Riff D, Negri FR, Hooghiemstra H, Lundberg J, Stadler T, Särkinen T, Antonelli A: Amazonia through time: Andean uplift, climate change, landscape evolution, and biodiversity. Science. 2010, 330 (6006): 927-931.
Ireland H, Kite G, Veitch N, Chase M, Schrire B, Lavin M, Linares J, Pennington RT: Biogeographical, ecological and morphological structure in a phylogenetic analysis of Ateleia (Swartzieae, Fabaceae) derived from combined molecular, morphological and chemical data. Bot J Linn Soc. 2010, 162: 39-53.
Simon MF, Grether R, de Queiroz LP, Skema C, Pennington RT, Hughes CE: Recent assembly of the Cerrado, a neotropical plant diversity hotspot, by in situ evolution of adaptations to fire. Proc Natl Acad Sci USA. 2009, 106 (48): 20359-20364.
Muellner AN, Pennington TD, Koecke AV, Renner SS: Biogeography of Cedrela (Meliaceae, Sapindales) in central and South america. Am J Bot. 2010, 97 (3): 511-518.
Särkinen T, Iganci JRV, Linares-Palomino R, Simon MF, Prado DE: Forgotten forests - issues and prospects in biome mapping using seasonally dry tropical forests as a case study. BMC Ecol. 2011, 11: 27-
Labruna MB, Camargo MA, Terrassini FA, Ferreira F, Schumaker TTS, Camargo EP: Ticks (Acari: Ixodidae) from the state of Rondônia, western Amazon, Brazil. Syst Appl Acarol. 2005, 10: 17-32.
Castagnolli KC, de Figueiredo LB, Santana DA, de Castro MB, Romano MA, Szabó MPJ: Acquired resistance of horses to Amblyomma cajennense (Fabricius, 1787) ticks. Vet Parasitol. 2003, 117 (4): 271-283.
Aragão HB: Ixodidas brasileiros e de alguns paizes limitrophes. Mem Inst Oswaldo Cruz (Rio de Janeiro). 1936, 31: 759-843.
Robinson LE: The genus Amblyomma. Ticks: A Monograph of the Ixodoidea. Part IV. Edited by: Nuttall GHF, Warburton C, Robinson LE. 1926, London: Cambridge University Press, 302-302.
Floch H, Fauran P: Ixodidés de la Guyane et des Antilles Françaises. Arch Inst Pasteur Guyane et Terr Inini. 1958, 4: 1-31.
Rojas R, Marini MA, Zanatta Coutinho AT: Wild birds as hosts of Amblyomma cajennense (Fabricius, 1787) (Acari: Ixodidae). Mem Inst Oswaldo Cruz (Rio de Janeiro). 1999, 94 (3): 315-322.
Hilburn LR, Sattler PW: Are tick populations really less variable and should they be?. Heredity. 1986, 57: 1133-1117.
Reichard MV, Kocan AA, Van Den Bussche RA, Barker RW, Wyckoff JH3rd, Ewing SA: Sequence variation of the ribosomal DNA second internal transcribed spacer region in two spatially-distinct populations of Amblyomma americanum (L.) (Acari: Ixodidae). J Parasitol. 2005, 91 (2): 260-263.
Beati L, Keirans JE: Analysis of the systematic relationships among ticks of the genera Rhipicephalus and Boophilus (Acari: Ixodidae) based on mitochondrial 12S ribosomal DNA gene sequences and morphological characters. J Parasitol. 2001, 87: 32-48.
Beati L, Patel J, Lucas-Williams H, Adakal H, Kanduma EG, Tembo-Mwase E, Krecek R, Mertins JW, Alfred JT, Kelly S, Kelly P: Phylogeography and demographic history of Amblyomma variegatum (Fabricius) (Acari: Ixodidae), the tropical bont tick. Vector Borne Zoonotic Dis. 2012, 12 (6): 514-525.
McLain DK, Wesson DM, Collins FH, Oliver JH Jr: Variation in ribosomal DNA internal transcribed spacers 1 among eastern populations of Ixodes scapularis (Acari: Ixodidae). J Med Entomol. 1995, 32 (3): 353-360.
Clement M, Posada D, Crandall KA: TCS: a computer program to estimate gene genealogies. Mol Ecol. 2000, 9 (10): 1657-1659.
Librado P, Rozas J: DnaSP v5: a software for comprehensive analysis of DNA polymorphism data. Bioinformatics. 2009, 25: 1451-1452.
Ramos-Onsins SE, Rozas J: Statistical properties of new neutrality tests against population growth. Mol Biol Evol. 2002, 19 (12): 2092-2100.
Fu YX, Li WH: Statistical tests of neutrality of mutations. Genetics. 1993, 133: 693-709.
Fu XY: Statistical tests of neutrality of mutations against population growth, hitchhiking, and background selection. Genetics. 1997, 147: 915-925.
Schneider S, Kueffer JM, Roessli D, Excoffier L: Arlequin: A Software for Population Genetic Data Analysis. Ver 2.000. 2000, Geneva: Genetics and Biometry Laboratory, University of Geneva
Maddison DR, Maddison WP: MacClade 4: Analysis of Phylogeny and Character Evolution. 2000, Sunderland, MA (CD-ROM): Sinauer Associates
Swofford DL: PAUP*: Phylogenetic Analysis Using Parsimony (and Other Methods) 4.0 Beta. 2000, Sunderland, MA: Sinauer Associates
Xia X, Xie Z: DAMBE: software package for data analysis in molecular biology and evolution. J Hered. 2001, 92 (4): 371-373.
Xia X, Xie Z, Salemi M, Chen L, Wang Y: An index of substitution saturation and its application. Mol Phylogenet Evol. 2003, 26: 1-7.
Xia X: Assessing substitution saturation with DAMBE. The Phylogenetic Handbook. Edited by: Lemey P, Salemi M, Vandamme AM. 2009, London: Cambridge University Press, 723-723.
Huelsenbeck JP, Ronquist F: MRBAYES: Bayesian inference of phylogenetic trees. Bioinformatics. 2001, 17 (8): 754-755.
Guindon S, Gascuel O: A simple, fast, and accurate algorithm to estimate large phylogenies by maximum likelihood. Syst Biol. 2003, 52 (5): 696-704.
Dereeper A, Guignon V, Blanc G, Audic S, Buffet S, Chevenet F, Dufayard JF, Guindon S, Lefort V, Lescot M, Claverie JM, Gascuel O: Phylogeny.fr: robust phylogenetic analysis for the non-specialist. Nucleic Acids Res. 2008, 36 (Web Server issue): W465-W469.
Posada D, Crandall KA: MODELTEST: testing the model of DNA substitution. Bioinformatics. 1998, 14: 817-818.
Bull JJ, Huelsenbeck JP, Cunningham CW, Swofford DL, Waddell PJ: Partitioning and combining data in phylogenetic analysis. Syst Biol. 1993, 42 (3): 384-397.
Drummond AJ, Ho S, Phillips MJ, Rambaut A: Relaxed phylogenetics and dating with confidence. PLoS Biol. 2006, 4: 699-710.
Drummond AJ, Rambaut A: BEAST: Bayesian evolutionary analysis by sampling trees. BMC Evol Biol. 2007, 7: 214-
Olson DM, Dinerstein E, Wikramanayake ED, Burgess ND, Powell GVN, Underwood EC, D’amico JA, Itoua I, Strand HE, Morrison JC, Loucks CJ, Allnutt TF, Ricketts TH, Kura Y, Lamoreux JF, Wettengel WW, Hedao P, Kassem KR: Terrestrial ecoregions of the world: a new map of life on Earth. BioScience. 2001, 51 (11): 933-938.
Casati S, Bernasconi MV, Gern L, Piffaretti JC: Assessment of intraspecific mtDNA variability of European Ixodes ricinus sensu stricto (Acari: Ixodidae). Infect Genet Evol. 2008, 8 (2): 152-158.
Norris DE, Klompen JS, Keirans JE: Black WC 4th: Population genetics of Ixodes scapularis (Acari: Ixodidae) based on mitochondrial 16S and 12S genes. J Med Entomol. 1996, 33: 78-89.
Zhang DX, Hewitt G: Insect mitochondrial control region: a review of its structure, evolution and usefulness in evolutionary studies. Biochem Syst Ecol. 1997, 25 (2): 99-120.
Schultheis AS, Weigt LA, Hendricks AC: Arrangement and structural conservation of the mitochondrial control region of two species of Plecoptera: utility of tandem repeat-containing regions in studies of population genetics and evolutionary history. Insect Mol Biol. 2002, 11 (6): 605-610.
Lessinger AC, Azeredo-Espin AM: Evolution and structural organisation of mitochondrial DNA control region of myiasis-causing flies. Med Vet Entomol. 2000, 14: 71-80.
Turner BJ, Elder JF Jr: Concerted evolution of repetitive DNA sequences in eukaryotes. Q Rev Biol. 1995, 70 (3): 297-320.
Graur D, Li WH: Fundamentals of Molecular Evolution. 2000, Cambridge: Sinauer Associates
Crease TJ: Ribosomal DNA evolution at the population level: nucleotide variation in the intergenic spacer arrays of Daphnia pulex. Genetics. 1995, 141: 1327-1337.
Fenton B, Malloch G, Moxey E: Analysis of eriophyid mite rDNA internal transcribed spacer sequences reveals variable simple sequence repeats. Insect Mol Biol. 1997, 6: 23-32.
Rich SM, Rosenthal BM, Telford SR3rd, Spielman A, Hartl DL, Ayala FJ: Heterogeneity of the internal transcribed spacer (ITS-2) region within individual deer ticks. Insect Mol Biol. 1997, 6 (2): 123-129.
Harris DJ, Crandall KA: Intragenomic variation within ITS1 and ITS2 of freshwater crayfishes (Decapoda: Cambaridae): implications for phylogenetic and microsatellite studies. Mol Biol Evol. 2000, 17 (2): 284-291.
Parkin EJ, Butlin RK: Within- and between-individual sequence variation among ITS1 copies in the meadow grasshopper Chorthippus parallelus indicates frequent intrachromosomal gene conversion. Mol Biol Evol. 2004, 21 (8): 1595-1601.
Keller I, Chintauan-Marquier IC, Veltsos P, Nichols RA: Ribosomal DNA in the grasshopper Podisma pedestris: escape from concerted evolution. Genetics. 2006, 174 (2): 863-874.
McLain DK, Wesson DM, Collins FH, Oliver JH Jr: Evolution of the rDNA spacer, ITS 2, in the ticks Ixodes scapularis and I. pacificus (Acari: Ixodidae). Heredity. 1995, 75 (3): 303-319.
Zahler M, Filippova NA, Morel PC, Gothe R, Rinder H: Relationships between species of the Rhipicephalus sanguineus group: a molecular approach. J Parasitol. 1997, 83 (2): 302-306.
Hlinka O, Murrell A, Barker SC: Evolution of the secondary structure of the rRNA internal transcribed spacer 2 (ITS2) in hard ticks (Ixodidae, Arthropoda). Heredity. 2002, 88 (4): 275-279.
Shaw M, Murrell A, Barker SC: Low intraspecific variation in the rRNA internal transcribed spacer 2 (ITS2) of the Australian paralysis tick, Ixodes holocyclus. Parasitol Res. 2002, 88 (3): 247-252.
Marrelli MT, Souza LF, Marques RC, Labruna MB, Matioli SR, Tonon AP, Ribolla PEM, Marinotti O, Schumaker TTS: Taxonomic and phylogenetic relationships between neotropical species of ticks from genus Amblyomma (Acari: Ixodidae) inferred from second internal transcribed spacer sequences of rDNA. J Med Entomol. 2007, 44 (2): 222-228.
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.
Nava S, Guglielmone AA: A meta-analysis of host specificity in Neotropical hard ticks (Acari: Ixodidae). Bull Entomol Res. 2013, 103 (2): 216-224.
Donato M, Posadas P, Miranda-Esquivel R, Ortiz Jaureguizar E, Cladera G: Historical biogeography of the Andean region: evidence from Listroderina (Coleoptera: Curculionidae: Rhytirrhinini) in the context of the South American geobiotic scenario. Biol J Linnean Soc. 2003, 80: 339-352.
Ron SR: Biogeographic area relationships of lowland Neotropical rainforest based on raw distributions of vertebrate groups. Biol J Linnean Soc. 2000, 71: 379-402.
Brumfield RT, Edwards SV: Evolution into and out of the Andes: a bayesian analysis of historical diversification in Thamnophilus antshrikes. Evolution. 2007, 61: 346-367.
Gregory-Wodzicki KM: Uplift history of the Central and Northern Andes: a review. Bull Geol Soc Am. 2000, 112: 1091-1105.
Antonelli A, Nylander JAA, Persson C, Sanmartín I: Tracing the impact of the Andean uplift on Neotropical plant evolution. Proc Natl Acad Sci USA. 2009, doi: 10.1073/pnas.0811421106
Hoorn C: The birth of the mighty Amazon. Sci Am. 2006, 294 (5): 52-59.
Hoorn C, Wesselingh FP, Hovikoski J, Guerrero J: The development of the Amazonian mega-wetland (Miocene; Brazil, Colombia, Peru, Bolivia). Amazonia, Landscape and Species Evolution: A Look into the Past. Edited by: Hoorn C, Wesselingh F. 2010, Hoboken, NJ: Wiley-Blackwell, 123-142.
Costa LP: The historical bridge between the Amazon and the Atlantic Forest of Brazil: a study of molecular phylogeography with small mammals. J Biogeogr. 2003, 30: 71-86.
Martins FM: Historical biogeography of the Brazilian Atlantic forest and the Carnaval-Moritz model of Pleistocene refugia: what do phylogeographical studies tell us?. Biol J Linnean Soc. 2011, 104: 499-509.
Werneck FP, Nogueira C, Colli GR, Costa GC, Sites JW Jr: Climatic stability in the Brazilian Cerrado: implications for the biogeographical connections of South American savannas, species richness and conservation in a biodiversity hotspot. J Biogeogr. 2012, 39: 1695-1706.
Hovikoski J, Wesselingh FP, Räsänen M, Gingras M, Vonhof HB: Marine influence in Amazonia: evidence from geological record. Amazonia, Landscape and Species Evolution: a Look into the Past. Edited by: Hoorn C, Wesselingh F. 2010, Hoboken, NJ: Wiley-Blackwell, 143-161.
Smith BT, Klicka J: The profound influence of the Late Pliocene Panamanian uplift on the exchange, diversification, and distribution of New World birds. Ecography. 2010, 33: 333-342.
Nogueira C, Ribeiro S, Costa GC, Colli GR: Vicariance and endemism in a Neotropical savanna hotspot: distribution patterns of Cerrado squamate reptiles. J Biogeogr. 2011, 38: 1907-1922.
Mora A, Baby P, Roddaz M, Parra M, Brusset S, Hermoza W, Espurt N: Tectonic history of the Andes and sub-Andean zones: implications for the development of the Amazon basin. Amazonia, Landscape and Species Evolution: A Look into the Past. Edited by: Hoorn C, Wesselingh FP. 2010, Hoboken, NJ: Wiley-Blackwell, 38-60.
Albert JS, Lovejoy NR, Crampton WGR: Miocene tectonism and the separation of cis- and trans-Andean river basins: evidence from the Neotropics. J South Am Earth Sci. 2006, 21: 14-27.
Davis JI, Nixon KC: Populations, genetic variation, and species delimitation of phylogenetic species. Syst Biol. 1992, 41 (4): 421-435.
Nava S, Beati L, Labruna MB, Caceres AG, Mangold AJ, Guglielmone AA: Reassessment of the taxonomic status of Amblyomma cajennense (Fabricius, 1787) with the description of three new species, Amblyomma tonelliae n. sp.,Amblyomma interandinum n. sp. and Amblyomma patinoi n. sp., and resurrection of Amblyomma mixtum Koch, 1844 and Amblyomma sculptum Berlese, 1888 (Ixodida: Ixodidae). Ticks and Tick-Borne Dis. 2013,, in press
Galvão MAM, Mafra CL, Moron C, Anaya E, Walker DH: Rickettsiosis of the genus Rickettsia in South America. Ann N Y Acad Sci. 2003, 990: 57-61.
Guedes E, Leite RC, Prata MCA, Pacheco RC, Walker DH, Labruna MB: Detection of Rickettsia rickettsii in the tick Amblyomma cajennense in a new Brazilian spotted fever-endemic area in the state of Minas Gerais. Mem Inst Oswaldo Cruz (Rio de Janeiro). 2005, 100 (8): 841-845.
Billings AN, Yu XJ, Teel PD, Walker DH: Detection of a spotted fever group rickettsia in Amblyomma cajennense (Acari: Ixodidae) in south Texas. J Med Entomol. 1998, 35 (4): 474-478.
Labruna MB, Whitworth T, Bouyer DH, McBride J, Camargo LMA, Camargo EP, Popov V, Walker DH: Rickettsia bellii and Rickettsia amblyommii in Amblyomma ticks from the State of Rondônia, Western Amazon, Brazil. J Med Entomol. 2004, 41 (6): 1073-1081.
Tomassone L, Nuñez P, Ceballos L, Gürtler R, Kitron U, Farber M: Detection of “Candidatus Rickettsia sp. strain Argentina” and Rickettsia bellii in Amblyomma ticks (Acari: Ixodidae) from Northern Argentina. Exp Appl Acarol. 2011, 52: 93-100.
This study was funded by NSF grant #1026146 to L. Beati and L.A. Durden. Our gratitude also goes to Jeff Alfred (USDA, Ames, ID) who provided us with samples from Texas. ICPS article number 1.
The authors declare that they have no competing interests.
LB and EJB designed the study, generated, analyzed, and interpreted data. Part of this study was the object of EJB’s master thesis project. LB, DBB, MBL, AG, SN, and LD participated in the conception of the study, acquisition and interpretation of data. AGC, CGC, RL, and JLHF were instrumental in data acquisition. All authors read and approved the final manuscript.
Electronic supplementary material
Additional file 1: 12SrDNA haplotype distribution. Haplotype distribution for 12SrDNA sequences. TX = Texas (U.S.), MX = Mexico, CR = Costa Rica, EC = Ecuador, CO = Colombia, FG = French Guiana, RO = Rondonia (Brazil, AF-MG = Atlantic Forest and Mato Grosso (Brazil), Yu = Yungas (Argentina), PA = Paraguay, CHS- CHO = Chaco Serrano and Chaco Occidental (Argentina), PE = Perú. Gray areas delimit haplotypes found clustered in the same TCS networks and phylogenetic clades. (XLSX 45 KB)
Additional file 2: COII haplotype distribution. Haplotype distribution for COII sequences. TX = Texas (U.S.), MX = Mexico, CR = Costa Rica, EC = Ecuador, CO = Colombia, FG = French Guiana, RO = Rondonia (Brazil, AF-MG = Atlantic Forest and Mato Grosso (Brazil), Yu = Yungas (Argentina), PA = Paraguay, CHS- CHO = Chaco Serrano and Chaco Occidental (Argentina), PE = Perú. Gray areas delimit haplotypes found clustered in the same TCS networks and phylogenetic clades. (XLSX 11 KB)
Additional file 3: DL haplotype distribution. Haplotype distribution for DL sequences. TX = Texas (U.S.), MX = Mexico, CR = Costa Rica, EC = Ecuador, CO = Colombia, FG = French Guiana, RO = Rondonia (Brazil, AF-MG = Atlantic Forest and Mato Grosso (Brazil), Yu = Yungas (Argentina), PA = Paraguay, CHS- CHO = Chaco Serrano and Chaco Occidental (Argentina), PE = Perú. Gray areas delimit haplotypes found clustered in the same TCS networks and phylogenetic clades. (XLSX 12 KB)
Additional file 4: ITS2 genotype distribution. Genotype distribution for ITS2 sequences. LTX = Texas (U.S.), MX = Mexico, CR = Costa Rica, EC = Ecuador, CO = Colombia, FG = French Guiana, RO = Rondonia (Brazil, AF-MG = Atlantic Forest and Mato Grosso (Brazil), Yu = Yungas (Argentina), PA = Paraguay, CHS- CHO = Chaco Serrano and Chaco Occidental (Argentina), PE = Perú. Gray areas delimit haplotypes found clustered in the same TCS networks and phylogenetic clades. (XLSX 43 KB)
Additional file 5: 12SrDNA maximum likelihood tree. Tree representing the relationships between A. cajennense inferred by ML analysis of 12SrDNA gene sequences. NW = Texas, Mexico, Cost Rica, Ecuador clade, NE = French Guiana and Rondonia (Brazil) clade, CO = Colombia, EA = Yungas Argentina + Atlantic Forest of Brazil, AR = Chaco (Argentina and Paraguay), PE = inter-Andean Valley of Perú. Numbers over the branches represent MP bootstrap values (1000 replicates), ML bootstrap values (100 replicates), and BA posterior probabilities respectively. (B) Unrooted TCS Network (95% parsimony cut-off). Same colors in A and B represent the same samples. (PPTX 86 KB)
Additional file 6: COII maximum likelihood tree. Tree representing the relationships between A. cajennense inferred by ML analysis of COII gene sequences. NW = Texas, Mexico, Cost Rica, Ecuador clade, NE = French Guiana and Rondonia (Brazil) clade, CO = Colombia, EA = Yungas Argentina + Atlantic Forest of Brazil, AR = Chaco (Argentina and Paraguay), PE = inter-Andean Valley of Perú. Numbers over the branches represent MP bootstrap values (1000 replicates), ML bootstrap values (100 replicates), and BA posterior probabilities respectively. (B) Unrooted TCS Network (95% parsimony cut-off). Same colors in A and B represent the same samples. (PPTX 76 KB)
Additional file 7: DL maximum likelihood tree. Tree representing the relationships between A. cajennense inferred by ML analysis of DL sequences. NW = Texas, Mexico, Cost Rica, Ecuador clade, NE = French Guiana and Rondonia (Brazil) clade, CO = Colombia, EA = Yungas Argentina + Atlantic Forest of Brazil, AR = Chaco (Argentina and Paraguay), PE = inter-Andean Valley of Perú. Numbers over the branches represent MP bootstrap values (1000 replicates), ML bootstrap values (100 replicates), and BA posterior probabilities respectively. (B) Unrooted TCS Network (95% parsimony cut-off). Same colors in A and B represent the same samples. (PPTX 82 KB)
Additional file 8: ITS2 maximum likelihood tree. Tree representing the relationships between A. cajennense inferred by ML analysis of ITS2 sequences. NW = Texas, Mexico, Cost Rica, Ecuador clade, NE = French Guiana and Rondonia (Brazil) clade, CO = Colombia, EA = Yungas Argentina + Atlantic Forest of Brazil, AR = Chaco (Argentina and Paraguay), PE = inter-Andean Valley of Perú. Numbers over the branches represent MP bootstrap values (1000 replicates), ML bootstrap values (100 replicates), and BA posterior probabilities respectively. (B) Unrooted TCS Network (95% parsimony cut-off). Same colors in A and B represent the same samples. (PPTX 76 KB)
Authors’ original submitted files for images
About this article
Cite this article
Beati, L., Nava, S., Burkman, E.J. et al. Amblyomma cajennense (Fabricius, 1787) (Acari: Ixodidae), the Cayenne tick: phylogeography and evidence for allopatric speciation. BMC Evol Biol 13, 267 (2013) doi:10.1186/1471-2148-13-267
- COII Sequence
- Eastern Cordillera
- Substitution Saturation
- Phylogeographical History
- Molecular Clock Hypothesis