Skip to main content

Phylogenetic analysis of the Australasian paralysis ticks and their relatives (Ixodidae: Ixodes: Sternalixodes)

Abstract

Background

The Australasian paralysis ticks and their relatives, Ixodes Latrielle, subgenus Sternalixodes Schulze, are some of the most important ticks in the region. However, very little is known about their phylogenetic relationships. The aim of this study was to elucidate the evolutionary relationships of members of the subgenus Sternalixodes by undertaking phylogenetic analyses of morphological and molecular datasets.

Methods

Adult females (nā€‰=ā€‰64) of Sternalixodes, including Ixodes anatis Chilton, 1904, Ixodes confusus Roberts, 1960, Ixodes cornuatus Roberts, 1960, Ixodes cordifer Neumann, 1908, Ixodes dendrolagi Wilson, 1967, Ixodes hirsti Hassall, 1931, Ixodes holocyclus Neumann, 1899, Ixodes myrmecobii Roberts, 1962 and Ixodes trichosuri Roberts, 1960, were examined morphologically. Subsequently, these Ixodes spp. were genetically characterised using cytochrome c oxidase subunit 1 (cox1) gene and the internal transcribed spacer 2 (ITS-2) of the rRNA. Both morphological and molecular datasets were analysed using various phylogenetic methods to assess the evolutionary relationship of various members of the subgenus Sternalixodes.

Results

Phylogenetic analyses of the cox1 sequences and morphological characters datasets revealed that the Australian and Papuan Sternalixodes formed a distinct clade with the New Zealand member of the group I. anatis positioned basally, in a separate clade. Ixodes holocyclus, I. cornuatus and I. myrmecobii formed a distinctive clade in both the cox1 and morphological phylogenies. However, based on phylogenetic analysis of the ITS-2 data, I. holocyclus formed a separate clade whereas I. cornuatus and I. myrmecobii grouped in a different clade.

Conclusions

The cox1 and morphological data suggest that the subgenus Sternalixodes is paraphyletic, and I. anatis is not a sternalixodid tick; hence, it should not be included in the subgenus. Based on the phylogenetic analyses of cox1 and ITS-2 sequences, it appears that I. myrmecobii and I. cornuatus are not subspecies of I. holocyclus. Although this study provided better insights into the taxonomic status of the subgenus Sternalixodes, a complete morphological and molecular (using multiple markers) phylogenetic analysis including all members of the subgenus would be required to more accurately elucidate the evolutionary relationships within the subgenus.

Background

Ticks (Arachnida: Ixodidae) are important ectoparasites of humans and animals and can cause direct (e.g., paralysis, anaemia) as well as indirect (e.g., transmission of pathogens) effects on their hosts [1ā€“3]. Ticks can be found on all continents and are known to feed on all types of terrestrial vertebrates, including mammals, birds, reptiles and amphibians [1, 3]. The life-cycle of ixodid ticks consists of four developmental stages, the egg, and three active parasitic stages, larva, nymph, and adult (male and female). Depending on the type (hard, Ixodidae or soft, Argasidae) and species of ticks, their life-cycle can vary significantly.

To date, 70 species (56 and 14 members of families Ixodidae and Argasidae, respectively) of ticks have been recorded from a variety of hosts (humans and domestic animalsā€‰=ā€‰16; mammals, reptiles and birdsā€‰=ā€‰54) from Australia [4]. Among different genera of hard ticks prevalent in Australia, Ixodes is arguably the most important and its members are known to transmit and/or harbour pathogens, including Rickettsia australis (the causative agent of Queensland tick typhus) in humans [5], and flaviviruses, bunyaviruses and Cercopithifilaria johnstoni Mackerras (Nematoda: Filarioidea) in wildlife [6ā€“8]. Some Ixodes species such as I. holocyclus and I. cornuatus can also cause paralysis in humans, domestic animals, and wildlife [9]. Ixodes holocyclus is known to have an immunoeffectory action on humans, causing tick bite anaphylaxis [10].

Members of the genus Ixodes have not been intensively examined (e.g., phylogeny using combined morphological and molecular datasets, biology, life-cycle) within Australia. Few life-cycles have been elucidated and the bionomics of very few species are understood. Of the subgenera of Ixodes, Sternalixodes has received most attention. The subgenus comprises nine members, including Ixodes anatis, I. confusus, I. cornuatus, I. cordifer, I. dendrolagi, I. hirsti, I. holocyclus, I. myrmecobii and I. trichosuri [11]. However, the majority of studies have focussed on I. holocyclus and I. cornuatus, a species morphologically similar to I. holocyclus, aiming to determine their distribution [1, 12], morphological and molecular identification [1, 13] and phylogenetic relationships [14]. A number of questions therefore remain to be answered regarding the members of the subgenus Sternalixodes. For instance, the status of I. myrmecobii as subspecies of I. holocyclus as proposed by Roberts [1] needs to be tested. To date, Australian paralysis ticks and their relatives have not been analysed using morphological as well as molecular phylogenetics. Being an important subgenus, Sternalixodes requires systematic investigations to address a number of taxonomic questions regarding the validity of its members. Therefore, this study was designed to elucidate the evolutionary relationships of members of the subgenus Sternalixodes by undertaking phylogenetic analyses of morphological and molecular datasets.

Methods

Tick collection and morphological identification

Female ticks (nā€‰=ā€‰74) used in this study were either available from The University of Melbourne (Ian Beveridge and Abdul Jabbar) or museums in Australia (South Australian Museum, Western Australian Museum, and Australian National Insect Collection), New Zealand (A. Heath, AgResearch, New Zealand), Papua New Guinea (Ifor L. Owen, National Veterinary Laboratory, Papua New Guinea) and South America (A. Guglielmone, Instituto Nacional de Technologia Agropecuaria, Argentina) (Fig.Ā 1; TableĀ 1). Following collection, each tick specimen was stored in 70% ethanol until used. Developmental stages of all species of Sternalixodes could not be examined as many are not yet described.

Fig. 1
figure 1

Collection sites for Ixodes species from Argentina, Australia, New Zealand and Uruguay, used in this study. Information linked to each unique number on the map is provided in TableĀ 1

Table 1 Specimens of Ixodes spp. used in molecular studies, and their sources and collection sites

For morphological identification, each tick was examined using a dissecting microscope (Olympus, Japan). In addition, electron micrographs were taken using a Hitachi TM3030 Tabletop Scanning Electron Microscope, Germany. All Australian and Papuan tick species were identified using keys by Roberts [1]; whereas I. anatis specimens were identified following Hardwick [15] and I. auritulus Neumann, 1904 specimens were identified by A. Heath and A. Guglielmone.

One or two legs were removed from each specimen using flame sterilized forceps and stored in 70% ethanol for molecular work.

DNA extraction, PCR amplification and DNA sequencing

Prior to DNA extraction, ethanol was removed and leg(s) of individual ticks were washed three times (30Ā min) in distilled H2O, and then ground using a plastic mortar. DNA was extracted using a DNeasy Blood and Tissue Kit (Qiagen, Hilden, Germany) following the protocol provided by the manufacturer.

Two loci (one mitochondrial and one nuclear ribosomal DNA) were PCR-amplified separately from each individual genomic DNA sample. The first locus (partial cox1 gene, ~850Ā bp) was amplified using the primers HCO2064 (5ā€²-GGT GGG CTC ATA CAA TAA ATC C-3ā€²) and HCOX1215 (5ā€²-GCC ATT TTA CCG CGA TGA-3ā€²); the second locus (partial second internal transcribed spacer, ITS-2; ~760Ā bp) was amplified employing primers ITS865 (5ā€²-CTC GCC TGA TCG TGA GGT CG-3ā€²) and ITS105 (5ā€²-GGT CGA ATT GCC CCT CCT CG-3ā€²) [14]. All PCRs were performed in a final volume of 50Ā Ī¼l, containing 10Ā mM Trisā€“HCl (pHĀ 8.4), 50Ā mM KCl, 3.5Ā mM of MgCl2, 200Ā Ī¼M of deoxynucleotide triphosphate, 100 pmol of each primer and 1 U of GoTaq polymerase (Promega, Madison, WI, USA) under the following cycling conditions: 94Ā Ā°C for 5Ā min (initial denaturation); 35Ā cycles of 94Ā Ā°C for 30Ā s (extension), 48Ā Ā°C (cox1) or 50Ā Ā°C (ITS-2) for 30Ā s (annealing) and 72Ā Ā°C for 50Ā s (extension), followed by final extension at 72Ā Ā°C for 5Ā min. For each set of PCRs, negative (no-DNA) and positive (I. holocyclus DNA) controls were included. No amplification was detected in any of the negative control reactions at any time during the study. Amplicons (5Ā Ī¼l) were examined on 1.5% agarose gels stained with ethidium bromide. Gels were examined using transillumination and were photographed using a GelDoc system (BioRad, Hercules, CA, USA). If amplicons were not detected on agarose gel, then semi-nested PCRs were used as follows: HCOX1240 (5ā€²-CCA CAA ATC ATA AAG ACA TTG G-3ā€²) was used in conjunction with HCO2064 to amplify cox1 and ITS130 (5ā€²-AGT TGT ACA TTG G-3ā€²) in conjunction with ITS865 was used to amplify ITS-2. PCR cycling conditions for semi-nested PCRs were same as used above.

For each locus, amplicon(s) representing each Ixodes species were purified using shrimp alkaline phosphatase and exonuclease 1 [16] prior to automated DNA sequencing (ABI3730XL automatic sequencer at Macrogen Cooperation, South Korea). Sequencing of the cox1 and ITS-2 region was conducted using the primers HCO2064 and HCOX1215 or HCO2064 and HCOX1240 (cox1) and ITS865 and ITS105 or ITS865 and ITS130 (ITS-2), in separate reactions. The quality of each sequence obtained was appraised using the program Geneious Pro 6.5 (Biomatters Ltd., Auckland, New Zealand) [17]. Partial cox1 sequences were identified by local alignment comparison (set reading frames) using amino acid sequences conceptually translated using an online tool http://www.ebi.ac.uk/Tools/st/emboss_transeq/ from the respective loci of the reference sequence of I. holocyclus are available from GenBank.

Phylogenetic analyses

For morphological phylogenetics, the character matrix was based on adult female specimens. All characters are morphological, collected by examining specimens using light and/or scanning electron microscopy. Characters that could not be scored with complete accuracy in some taxa were excluded from the analysis. Morphological data were analysed employing Maximum Parsimony (MP) in TnT [18], gaps were treated as missing characters, and bootstrap replicates and maximum trees were set at 10,000. In addition, data were analysed using Bayesian Inference (BI) by employing the Markov K model in MrBayes 3.2.6 [19ā€“21]. Lset rates were set to gamma and coding was set to variable. Four simultaneous tree-building chains were used to calculate posterior probabilities (pp) for 2,000,000 generations, saving every 100th tree produced. Based on the final 75% of trees generated, a consensus tree was constructed. Ixodes tasmani Neumann, 1899 was used as the outgroup.

For molecular phylogenetics, nucleotide sequences were aligned using the MUSCLE V 3.8.31 program [22] and adjusted manually employing the program Mesquite V 3.03 [23]. Based on pairwise comparisons, sequence differences were calculated using the program MEGA 6.0. [24]. Two separate datasets representing cox1 and ITS-2 were compiled, together with reference sequences from GenBank [14, 25, 26]; I. tasmani and I. uriae White, 1852 were used as the outgroups, respectively. Both cox1 and ITS-2 sequences were aligned over 519 and 610Ā bp, respectively, and adjusted manually as described above. Phylogenetic analyses were performed on individual cox1 and ITS-2 datasets using Maximum Likelihood (ML), Neighbour-Joining (NJ) and BI methods. The ML and NJ analyses were performed using MEGA 6.0. and the nodes were tested for robustness with 10,000 bootstrap replicates. The data format was set to DNA and gaps were treated as missing data (10,000 bootstrap replicates, Max. trees was set at 10,000). The likelihood parameters for the BI (TIM2+I+G for pCXO1 and TVM+G for ITS-2) and ML (Tamura 3-parameter model for both cox1 and ITS-2) analyses were selected based on the Akaike Information Criterion (AIC) test in jModeltest v2.1.5 [27]. The BI was conducted, using Monte Carlo Markov Chain (MCMC) analysis in MrBayes 3.1.2. Four simultaneous tree-building chains were used to calculate posterior probabilities (pp) for 2,000,000 generations, saving every 100th tree produced. Based on the final 75% of trees generated, a consensus tree was constructed.

The phylogenetic trees produced for both morphological and molecular datasets were visually compared separately for concordance in their topologies.

Results

Morphological characterisation

Out of 74 individual specimens of female Ixodes examined, 64 belonged to the subgenus Sternalixodes, including I. anatis (nā€‰=ā€‰7), I. dendrolagi (nā€‰=ā€‰3), I. cordifer (nā€‰=ā€‰5), I. cornuatus (nā€‰=ā€‰5), I. hirsti (nā€‰=ā€‰10), I. holocyclus (nā€‰=ā€‰13), I. myrmecobii (nā€‰=ā€‰17) and I. trichosuri (nā€‰=ā€‰4); whereas, remaining 10 belonged to two subgenera Endopalpiger Schulze (I. tasmani; nā€‰=ā€‰6) and Multidentatus Neumann (I. auritulus; nā€‰=ā€‰4).

Character states are presented in TableĀ 2, and the morphological data matrix is provided in TableĀ 3. In addition, principal features of the capitulum used as characters are shown in Fig.Ā 2.

Table 2 List of morphological characters (character numbers, name of character, character states)
Table 3 Morphological character matrix of character states for each taxon used to construct morphological phylogeny
Fig. 2
figure 2

Principal features of the capitulum used as characters in TableĀ 2. Capituli of a, Ixodes cornuatus (ventral view); b, Ixodes hirsti (ventral view); c, Ixodes tasmani (ventral view); d, Ixodes tasmani (dorsal view), and e, Ixodes auritulus (dorsal view)

Molecular characterisation

PCR amplification was successful for 27 (out of 64) genomic DNA samples extracted from individual tick specimens (TableĀ 1). Considerable variation in the size (~650 to 750Ā bp) of amplicons (nā€‰=ā€‰15) for the ITS-2 was detected on agarose gel, whereas the amplicon size (~700Ā bp) for cox1 (nā€‰=ā€‰27) did not differ. DNA sequencing of amplicons for both loci revealed 27 and 15 unique sequences for cox1 and ITS-2, respectively. Sequence length, G+C content, pairwise differences and GenBank accession numbers for cox1 (KY213767ā€“KY213793) and ITS-2 (KY213752ā€“KY213766) sequences are given in TableĀ 4. The length of cox1 sequences for each tick species was 674Ā bp, whereas that of ITS-2 ranged from 630 to 704Ā bp. Among various members of Sternalixodes, the highest genetic variation was detected in I. cornuatus (number of sequencesā€‰=ā€‰4; pairwise differences 0.2ā€“12%) followed by I. anatis (nā€‰=ā€‰2; 1.2%), I. myrmecobii (nā€‰=ā€‰6; 0.2ā€“1.1%), I. auritulus (nā€‰=ā€‰2; 1.1%), I. holocyclus (nā€‰=ā€‰5; 0.2ā€“0.9%), I. hirsti (nā€‰=ā€‰3; 0.2ā€“0.3%) and I. trichosuri (nā€‰=ā€‰2; 0.3%) (TableĀ 4). Based on ITS-2 sequences, multiple sequences for individual ticks were obtained only for I. holocyclus and I. myrmecobii and their pairwise differences were 1.4ā€“8.6% and 0.8ā€“3.6%, respectively (TableĀ 4).

Table 4 Characteristics of sequences of Ixodes species determined in this study. GenBank accession numbers, G+C content and length of each sequence, and pairwise differences for each species with more than one specimens are provided

Phylogenetic analyses

The topology of the phylogenetic trees generated for morphological data employing BI and MP methods were similar (data not shown); hence, the MP tree is presented here, with nodal support values given for both methods (Fig.Ā 3). The morphological phylogram showed six main clades, clade numbers including taxa of the preceding clade. Ixodes confusus and I. dendrolagi grouped together in clade 1, with moderate statistical support (posterior probability for BI: 0.99; bootstrap value for MP: 87%) (Fig.Ā 3). Ixodes hirsti, I. trichosuri, I. auritulus and I. anatis each formed a clade (2, 3, 5 and 6, respectively), with no to high statistical support (Fig.Ā 3). The common Australian paralysis tick, I. holocyclus, and I. cordifer, I. cornuatus and I. myrmecobii formed clade 4, with low to moderate statistical support (0.92, 79%).

Fig. 3
figure 3

Morphological phylogram of specimens of Sternalixodes from Argentina, Australia, New Zealand and Uruguay. The relationships were inferred based on 38 morphological characters using Maximum Parsimony (MP) and Bayesian Inference (BI) methods. Ixodes tasmani was used as the outgroup. There was a concordance in the topology between this MP tree and that produced using BI (not shown). Nodal support (from left to right) is given as a posterior probability for BI and bootstrap values for MP. For simplicity, each clade number includes taxa in the preceding clade. The scale-bar indicates the number of inferred substitutions per character

Molecular phylogenetic analyses revealed that the topology of trees generated from the cox1 (aligned over 519 positions) and ITS-2 (608 positions) sequence data were similar using BI, NJ and ML (data not shown); hence, only the NJ trees for both loci are presented here (Figs.Ā 4 and 5). The cox1 tree had three major clades (Fig.Ā 4) in which I. cornuatus, I. holocyclus and I. myrmecobii formed Clade 1, with mixed statistical support (posterior probability for BI: 0.90; bootstrap value for NJ and ML: 97 and 86%). Individually, five cox1 sequences of I. holocyclus determined herein (GenBank accession nos. KY213779ā€“KY213782) grouped together with those previously published from Australia, with strong statistical support (0.99, 100, 99%) (Fig.Ā 4). All six cox1 sequences of I. myrmecobii grouped together with strong statistical support (0.99, 100, 96%), whereas four cox1 sequences of I. cornuatus found in this study formed two sub-clades with strong statistical support (1.0, 100, 99%) in which two sequences from Tasmania (KY213792 and KY213793) grouped outside the other two sequences from this study (KY213790 and KY213791) as well as previously published sequences (Fig.Ā 4). Clade 2 contained I. dendrolagi, I. hirsti and I. trichosuri but without statistical support (0.68, 58, 51%; Fig.Ā 4). However, individual sequences of I. hirsti determined here (KY213773ā€“KY213775) formed a separate sub-clade compared with previously published sequences of this species. Ixodes auritulus and I. anatis formed Clade 3 with weak to moderate statistical support (0.98, 67, 70%) (Fig.Ā 4).

Fig. 4
figure 4

Genetic relationships of specimens of Sternalixodes from Argentina, Australia, New Zealand and Uruguay. The relationships were inferred based on phylogenetic analyses of the cytochrome c oxidase subunit 1 (cox1) sequence data determined herein (bold) using Bayesian Inference (BI), distance-based Neighbor Joining (NJ) and Maximum Likelihood (ML) methods. Previously published sequences of Ixodes species were obtained from GenBank (see accession numbers). Ixodes tasmani was used as the outgroup. There was a concordance in the topology between this NJ tree and those produced using BI and ML (not shown). Nodal support (from left to right) is given as a posterior probability for BI and bootstrap values for NJ and ML. The scale-bar indicates the number of inferred substitutions per nucleotide site

Fig. 5
figure 5

Genetic relationships of specimens of Sternalixodes from Argentina, Australia, New Zealand and Uruguay. The relationships were inferred based on phylogenetic analyses of the second internal transcribed spacer (ITS-2) sequence data determined herein (bold) using Bayesian Inference (BI), distance-based Neighbor Joining (NJ) and Maximum Likelihood (ML) methods. Previously published sequences of Ixodes species were obtained from GenBank (see accession numbers). Ixodes uriae was used as the outgroup. There was a concordance in the topology between this NJ tree and those produced using BI and ML (not shown). Nodal support (from left to right) is given as a posterior probability for BI and bootstrap values for NJ and ML. The scale-bar indicates the number of inferred substitutions per nucleotide site

The ITS-2 tree contained four major clades (Fig.Ā 5). However, the composition of some clades was different from that found in the cox1 tree. For example, in the ITS-2 tree, Clade 1 contained only I. holocyclus with strong statistical support (1.0, 100, 99%) whereas I. cornuatus and I. myrmecobii formed a separate clade (Clade 3) with strong statistical support (1.0, 100, 99%) (Fig.Ā 5). Similarly in the cox1 tree, I. dendrolagi, I. hirsti and I. trichosuri formed a separate clade (Clade 2) with strong statistical support (1.0, 99, 95%) whereas I. anatis formed a separate clade outside Australian and Papuan Ixodes spp. (Fig.Ā 5).

Discussion

This study addressed the evolutionary relationships amongst the species of Sternalixodes, using both morphological and molecular phylogenetic methods. The cox1 and morphological data suggest that the subgenus is paraphyletic with I. anatis which is congruent with a previous suggestion by Heath [28].

The topology was similar in the cox1 and morphological trees, with the Australian and Papuan Sternalixodes forming a distinct clade and the New Zealand member of the group I. anatis positioned basally, in a separate clade. Ixodes holocyclus, I. cornuatus and I. myrmecobii formed a distinctive clade in both the cox1 and morphological phylogenies. This pattern supports comments made by Roberts [1], who suggested that these three species were closely related and that I. myrmecobii and I. cornuatus may be subspecies of I. holocyclus. Ixodes hirsti, I. trichosuri and I. dendrolagi formed a separate clade distinct from the I. holocyclus species group. It is possible that the resulting tree would show I. dendrolagi forming a distinct clade with these northern species rather than grouping with I. hirsti. The cox1 phylogeny suggests that the outgroup I. tasmani contains a cryptic species based on the long branch lengths of the two I. tasmani samples (see Fig.Ā 4) as previously proposed by Roberts [1].

While the morphological and cox1 trees showed similar topologies, there were some differences. The position of I. auritulus was different between the cox1 tree and the morphological tree, while I. auritulus and I. anatis formed a distinct clade in the cox1 phylogram. However, in the morphological tree I. auritulus formed a clade with the Australian/Papuan Sternalixodes. It is likely that the groups are only distantly related as I. anatis, I. auritulus and the Australian Sternalixodes are all morphologically and molecularly distinct. More extensive morphological and molecular phylogenetic analyses are required to adequately illuminate the evolutionary relationships between the three groups. This would require more extensive morphological character sets as well as the use of other molecular markers such as 16S, 28S or complete mitochondrial genomes.

Topological differences were also present between the two molecularly derived trees (see Figs.Ā 4 and 5). The members of the I. holocyclus species group did not form a distinct clade in the ITS-2 tree as was seen in the cox1 tree, but instead formed two separate clades comprising I. cornuatus and I. myrmecobii in one and I. holocyclus in the other. The morphologically similar species, I. holocyclus and I. cornuatus grouped in highly divergent clades contrary to morphological evidence, as well as inferences from previous study of the interspecific relationship of these species by Song et al. [14]. This was likely due to the conserved nature of this marker in ticks. Song et al. [14] commented on the conserved nature of ITS-2 in Sternalixodes noting that intraspecific variation between I. holocyclus and I. cornuatus was as low as 0.19%. Despite suggestions that ITS-2 is suitable for inferring evolutionary relationships in ticks [14], it appears that it may not be suited for revealing the relationships between more distantly related species within subgenera.

Within the cox1 phylogeny, I. hirsti was divided into two distinct subclades (see Fig.Ā 4). The GenBank sequences utilised were from ticks collected in South Australia [29], while the sequences obtained during this study were from Victorian specimens. This pattern in the cox1 sequences coupled with the geographic difference between the two groups suggests the group may be undergoing genetic differentiation.

The results of this study inform a number of historical questions and uncertainties concerning the subgenus Sternalixodes. Previously, Song et al. [14] used molecular techniques to assess the validity of the I. holocyclus species group. In the present study, based on the cox1 and ITS-2 trees, it appears that I. myrmecobii and I. cornuatus are not subspecies of I. holocyclus, each being a valid species. This contradicts the suggestion made by Roberts [1] that these two species may be subspecies of the widespread I. holocyclus. The results of this study are congruent with the results of Song et al. [14] and Jackson et al. [13]. However neither of these studies included the Western Australian species, I. myrmecobii.

The findings of this study also provided insights into the validity of Sternalixodes as a subgenus. The apparent paraphyletic status of the subgenus based on the position of I. anatis in the cox1 and morphological phylograms validates the suggestion by Heath [28] that the species should not be included in Sternalixodes. Heath [28] made the suggestion, citing the morphology of I. anatis as being incongruent with the morphological definition of Sternalixodes. Based on both morphological and cox1 data, it appears that I. anatis is not a sternalixodid tick and should not be included in the subgenus. This species does not meet the criteria of any of the subgenera of Ixodes defined by Clifford et al. [11]. However, Clifford et al. [11] noted that the classification of the subgenera of Ixodes was inaccurate in some situations, especially with regard to the subgenus Ixodes. An extensive examination of all existing subgenera using molecular and morphological data should be made in future to provide a more accurate hypothesis of the evolutionary relationships between the subgenera and the validity of the species within them. Based on the distinctive morphology of I. anatis and the fact that it does not meet the diagnostic requirements of any of the known subgenera of Ixodes, it may require the erection of a new subgenus. However more extensive molecular data should be accumulated and examined before this can occur.

Although seven of the nine species of Sternalixodes were examined in this study, molecular sequences and morphologically complete specimens were not located for I. confusus and I. cordifer. Although these species are most likely members of Sternalixodes, a complete molecular phylogenetic analysis including these species would be desirable to more accurately illuminate the evolutionary relationships within the subgenus.

As I. myrmecobii clusters within the I. holocyclus species group, a set of ticks known to cause paralysis, the question of its ability to also induce paralysis is raised. Tick induced paralysis has been extensively studied on the east coast of Australia; however, little information exists concerning ticks in Western Australia, let alone tick paralysis in Western Australia [9]. Studies of I. myrmecobii should be undertaken to determine if this species can induce paralysis. Roberts [30] and Kemp [31] noted that I. hirsti has been recorded to cause paralysis. Kemp [31] also proposed that all sternalixodid ticks may be capable of causing paralysis. As I. hirsti clustered with I. trichosuri and I. dendrolagi within the morphological and cox1 phylogeny, it is possible that these species may also be capable of inducing paralysis, however, this should be investigated.

Conclusion

In conclusion, the cox1 and morphological data suggest that the subgenus Sternalixodes is paraphyletic, and I. anatis should not be included in this subgenus. Based on the phylogenetic analyses of cox1 and ITS-2 sequences, it appears that I. myrmecobii and I. cornuatus are not subspecies of I. holocyclus, each being a valid species. Although this study has improved insights into the taxonomic status of the subgenus Sternalixodes, a complete morphological and molecular (using multiple markers) phylogenetic analysis including all nine species of the subgenus would be desirable to more accurately illuminate the evolutionary relationships within the subgenus.

Abbreviations

AIC:

Akaike information criterion

BI:

Bayesian inference

cox1:

Cytochrome c oxidase subunit 1 gene

ITS-2:

Internal transcribed spacer 2

MCMC:

Monte Carlo Markov Chain

ML:

Maximum likelihood

NJ:

Neighbour-joining

References

  1. Roberts FHS. Australian ticks. Victoria: CSIRO Publishing; 1970.

    Google ScholarĀ 

  2. Service M. Medical entomology for students. 3rd ed. Cambridge: Cambridge University Press; 2004.

    Google ScholarĀ 

  3. Sonenshine DE, Roe RM. Biology of ticks, vol 2. 2nd ed. New York: Oxford University Press; 2014.

    Google ScholarĀ 

  4. Barker SC, Walker AR, Campelo D. A list of the 70 species of Australian ticks; diagnostic guides to and species accounts of Ixodes holocyclus (paralysis tick), Ixodes cornuatus (southern paralysis tick) and Rhipicephalus australis (Australian cattle tick); and consideration of the place of Australia in the evolution of ticks with comments on four controversial ideas. Int J Parasitol. 2014;44:941ā€“53.

    ArticleĀ  PubMedĀ  Google ScholarĀ 

  5. Sexton DJ, Dwyer B, Kemp R, Graves S. Spotted fever group rickettsial infections in Australia. Rev Infect Dis. 1991;13:876ā€“86.

    ArticleĀ  CASĀ  PubMedĀ  Google ScholarĀ 

  6. St George TD, Doherty RL, Carley JG, Filippich C, Brescia A, Casals J, et al. The isolation of arboviruses including a new flavivirus and a new bunyavirus from Ixodes (Ceratixodes) uriae (Ixodoidea, Ixodidae) collected at Macquarie Island, Australia, 1975ā€“1979. Am J Trop Med Hyg. 1985;34:406ā€“12.

    CASĀ  PubMedĀ  Google ScholarĀ 

  7. St George TD, Standfast HA, Doherty RL, Carley JG, Fillipich C, Brandsma J. The isolation of Saumarez Reef virus, a new flavivirus, from bird ticks Ornithodoros capensis and Ixodes eudyptidis in Australia. Aust J Exp Biol Med Sci. 1977;55:493ā€“9.

    ArticleĀ  CASĀ  PubMedĀ  Google ScholarĀ 

  8. Spratt DM, Haycock P. Aspects of the life history of Cercopithifilaria johnstoni (Nematoda:Filarioidea). Int J Parasitol. 1988;18:1087ā€“92.

    ArticleĀ  CASĀ  PubMedĀ  Google ScholarĀ 

  9. Barker SC, Walker AR. Ticks of Australia. The species that infest domestic animals and humans. Zootaxa. 2014;3816:1ā€“144.

    ArticleĀ  Google ScholarĀ 

  10. Brown AF, Hamilton DL. Tick bite anaphylaxis in Australia. J Accid Emerg Med. 1998;15:111ā€“3.

    ArticleĀ  CASĀ  PubMedĀ  PubMed CentralĀ  Google ScholarĀ 

  11. Clifford CM, Sonenshine DE, Keirans JE, Kohls GM. Systematics of the subfamily Ixodinae (Acarina: Ixodidae). 1. The subgenera of Ixodes. Ann Entomol Soc Am. 1973;66:489ā€“500.

    ArticleĀ  Google ScholarĀ 

  12. Jackson J, Beveridge I, Chilton NB, Andrews RH. Distributions of the paralysis ticks Ixodes cornuatus and Ixodes holocyclus in south-eastern Australia. Aust Vet J. 2007;85:420ā€“4.

    ArticleĀ  CASĀ  PubMedĀ  Google ScholarĀ 

  13. Jackson J, Chilton NB, Beveridge I, Morris M, Andrews RH. Genetic variation within the ticks Ixodes holocyclus and Ixodes cornuatus from south-eastern Australia. Int J Parasitol. 2000;30:1159ā€“66.

    ArticleĀ  CASĀ  PubMedĀ  Google ScholarĀ 

  14. 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:871ā€“80.

    ArticleĀ  CASĀ  PubMedĀ  Google ScholarĀ 

  15. Hardwick S. Electronic polytomous and dichotomous keys to the genera and species of hard ticks (Acari: Ixodidae) present in New Zealand. Syst Appl Acarol. 2010;15:163ā€“83.

    ArticleĀ  Google ScholarĀ 

  16. Werle E, Schneider C, Renner M, Vƶlker M, Fiehn W. Convenient single-step, one tube purification of PCR products for direct sequencing. Nucleic Acids Res. 1994;22:4354.

    ArticleĀ  CASĀ  PubMedĀ  PubMed CentralĀ  Google ScholarĀ 

  17. Kearse M, Moir R, Wilson A, Stones-Havas S, Cheung M, Sturrock S, et al. Geneious Basic: an integrated and extendable desktop software platform for the organization and analysis of sequence data. Bioinformatics. 2012;28:1647ā€“9.

    ArticleĀ  PubMedĀ  PubMed CentralĀ  Google ScholarĀ 

  18. Goloboff PA, Farris JS, Nixon KC. TNT, a free program for phylogenetic analysis. Cladistics. 2008;24:774ā€“86.

    ArticleĀ  Google ScholarĀ 

  19. Lewis P. A likelihood approach to estimating phylogeny from discrete morphological character data. Syst Biol. 2001;50:913ā€“25.

    ArticleĀ  CASĀ  PubMedĀ  Google ScholarĀ 

  20. Huelsenbeck JP, Ronquist F. MRBAYES: Bayesian inference of phylogenetic trees. Bioinformatics. 2001;17:754ā€“5.

    ArticleĀ  CASĀ  PubMedĀ  Google ScholarĀ 

  21. Ronquist F, Huelsenbeck JP. MRBAYES 3: Bayesian phylogenetic inference under mixed models. Bioinformatics. 2003;19:1572ā€“4.

    ArticleĀ  CASĀ  PubMedĀ  Google ScholarĀ 

  22. Edgar RC. MUSCLE: multiple sequence alignment with high accuracy and high throughput. Nucleic Acids Res. 2004;32:1792ā€“7.

    ArticleĀ  CASĀ  PubMedĀ  PubMed CentralĀ  Google ScholarĀ 

  23. Maddison WP, Maddison DR. Mesquite: a modular system for evolutionary analysis. Version 3.10, 2016. http://mesquiteproject.org

  24. Tamura K, Stecher G, Peterson D, Filipski A, Kumar S. MEGA 6: molecular evolutionary genetics analysis version 6.0. Mol Biol Evol. 2013;30:2725ā€“9.

    ArticleĀ  CASĀ  PubMedĀ  PubMed CentralĀ  Google ScholarĀ 

  25. Shao R, Aoki Y, Mitani H, Tabuchi N, Barker SC, Fukunaga M. The mitochondrial genomes of soft ticks have an arrangement of genes that has remained unchanged for over 400 million years. Insect Mol Biol. 2004;13:219ā€“24.

    ArticleĀ  CASĀ  PubMedĀ  Google ScholarĀ 

  26. Hammer JF, Emery D, Bogema DR, Jenkins C. Detection of Theileria orientalis genotypes in Haemaphysalis longicornis ticks from southern Australia. Parasit Vectors. 2015;8:229.

    ArticleĀ  PubMedĀ  PubMed CentralĀ  Google ScholarĀ 

  27. Darriba D, Taboada GL, Doalla R, Posada D. jModelTest 2: more models, new heuristics and parallel computing. Nat Methods. 2012;9:772.

    ArticleĀ  CASĀ  PubMedĀ  PubMed CentralĀ  Google ScholarĀ 

  28. Heath ACG. Zoogeography of the New Zealand tick fauna. Tuatara. 1977;23:26ā€“38.

    Google ScholarĀ 

  29. Chapman TW, Marando L, Oorebeek M, Kleindorfe S. Genetic structure in ixodid ticks from Kangaroo Island and the South Australian mainland. Aust J Entomol. 2009;48:40ā€“6.

    ArticleĀ  Google ScholarĀ 

  30. Roberts FHS. Tick paralysis in South Australia. Aust Vet J. 1961;37:440.

    ArticleĀ  Google ScholarĀ 

  31. Kemp DH. Identity of Ixodes holocyclus and other paralysis ticks in Australia. Aust Adv Vet Sci. 1979;71:71.

    Google ScholarĀ 

Download references

Acknowledgements

We are grateful to Tony Chiovitti and other staff from the Gene Technology Access Centre (GTAC), Melbourne for access to the scanning electron microscope. We are also thankful to Allen Heath (Ag Research, New Zealand) and Alberto Guglielmone (Instituto Nacional de Technologia Agropecuaria, Argentina), for both identifying and providing tick specimens for both morphological and molecular examination. In addition, we are indented to Mark Harvey (Western Australian Museum), Ifor Owen (National Veterinary Laboratory, PNG) and Bruce Halliday (Australian National Insect Collection) for the provision of tick specimens.

Funding

This project was funded by the Early Career Researcher grant from The University of Melbourne, Australia (AJ).

Availability of data and materials

The data supporting the conclusions of this article are included within the article. The nucleotide sequences generated during this study are deposited in GenBank under the accession numbers KY213752ā€“KY213793.

Authorsā€™ contributions

MK, IB and AJ conceived the project and participated in the study design. MK carried out the laboratory work, data analyses, interpretation of data, and also drafted the manuscript, with guidance from co-authors. MK, AVK and AJ undertook phylogenetic analyses. IB and AJ participated in data interpretation and provided critical inputs on the draft manuscript. All authors read and approved the final manuscript.

Competing interests

The authors declare that they have no competing interests.

Consent for publication

Not applicable.

Ethics approval and consent to participate

Not applicable.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Abdul Jabbar.

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Kwak, M.L., Beveridge, I., Koehler, A.V. et al. Phylogenetic analysis of the Australasian paralysis ticks and their relatives (Ixodidae: Ixodes: Sternalixodes). Parasites Vectors 10, 122 (2017). https://doi.org/10.1186/s13071-017-2045-4

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s13071-017-2045-4

Keywords