- Open Access
Automatic barcode gap discovery reveals diverse clades of Rhipicephalus spp. and Haemaphysalis spp. ticks from small mammals in 'Asir, Saudi Arabia
Parasites & Vectors volume 14, Article number: 541 (2021)
The ixodid tick genera Rhipicephalus and Haemaphysalis contain several species of medical and/or veterinary importance, but their diversity in some regions of the world remains under-explored. For instance, very few modern studies have been performed on the taxonomy of these genera on the Arabian Peninsula.
In this study, we trapped small mammals in the 'Asir Mountains of south-western Saudi Arabia and collected tick specimens for morphological examination and molecular barcoding, targeting three mitochondrial loci: cox1, 16S rRNA and 12S rRNA.
We obtained a total of 733 ticks (608 Haemaphysalis spp. and 125 Rhipicephalus spp.) from 75 small mammal hosts belonging to six species. All tick specimens were immature except for nine adults recovered from a hedgehog (Paraechinus aethiopicus). Morphologically, the Rhipicephalus ticks resembled R. camicasi, but the Haemaphysalis ticks showed differences in palp morphology compared with species previously described from Saudi Arabia. Phylogenetic analysis and automatic barcode gap discovery identified a novel clade of Rhipicephalus sp. representing most of the nymphs. This was most closely related to R. leporis, R. guilhoni and R. linnaei. The adult ticks and a small proportion of nymphs clustered with R. camicasi sequences from a previous study. Finally, the Haemaphysalis nymphs formed two distinct clades that were clearly separated from all reference sequences but closest to some African species.
This apparent high level of tick diversity observed in a single study site of only ~ 170 km2, on a relatively small number of hosts, highlights the potential for the discovery of new tick species on the Arabian Peninsula.
The Ixodidae (hard ticks) is by far the most speciose family of ticks, with over 700 validly described species . Until comparatively recently, our understanding of the relationships between tick species was founded almost exclusively on analysis of morphological features. Due to their large and complex genomes, whole nuclear genome data for ticks remain sparse  compared with insects of medical and/or veterinary importance, and investigations of possible species complexes within morphologically similar tick groups have proceeded slowly. However, molecular confirmation of tick species identity using mitochondrial barcodes and phylogenetic analyses based on concatenated mitochondrial loci, or more recently, nucleotide and amino acid datasets from whole mitogenomes, have begun to revolutionise both the taxonomic status of closely related species and the higher-level relationships between tick genera and families [3,4,5,6,7].
There have been increasing reports of discordance between morphological features and genetic characteristics within ixodid taxa, including Ixodes and Rhipicephalus; two of the most intensely studied genera of medical and veterinary importance. For instance, a recent study showed that certain Australian Ixodes spp. specimens were highly divergent genetically but morphologically indistinguishable, whereas other specimens were morphologically distinct but poorly resolved genetically . Moreover, two of the most important Rhipicephalus spp. globally, the Asian blue tick, R. microplus, and the brown dog tick, R. sanguineus, are each now known to be formed of several distinct lineages, which are becoming recognised as distinct species [9,10,11,12,13,14,15]. The highly diverse genus Haemaphysalis has been the subject of far fewer molecular studies, although substantial discrepancies between morphology-based classification and molecular characteristics have recently been noted for this taxon too [4, 7, 16]. One generic approach to resolving species diversity using objective molecular criteria is automatic barcode gap discovery (ABGD), which is founded on the principle that the genetic divergence should be smaller within species than between species . This allows a confidence limit to be assigned to intraspecific divergence, thus partitioning gene sequences into bins or operational taxonomic units (OTUs). The ABGD approach and related methods are gaining in popularity in molecular studies of ticks worldwide [18,19,20].
One geographic region in which the diversity of Ixodidae is under-explored is the Arabian Peninsula. A key to the ticks of Yemen was published by Hoogstraal and Kaiser  and for Saudi Arabia by Hoogstraal et al. . Recent reports of ticks from the region have focused primarily on identification of species collected from domestic animals and pathogen screening [23,24,25], with a smaller number of studies on tick specimens obtained from wild hosts [26,27,28]. Importantly, to the best of our knowledge, no molecular data from ticks collected from wildlife in Saudi Arabia have been published to date. Here, we identify a novel clade of Rhipicephalus spp. ticks feeding on rodents in the 'Asir Mountains of south-western Saudi Arabia, which is molecularly distinct from sympatric specimens that cluster with R. camicasi. We also present preliminary evidence for two novel clades of Haemaphysalis spp. ticks infesting the same hosts.
Field site and small mammal trapping
Details of the study site and small mammal collection have been published previously . Briefly, small mammals were trapped overnight in the summers of 2016 and 2017 near three villages (Al Ous', Alogl and Wosanib) on the upper escarpment of the 'Asir Mountains in south-western Saudi Arabia, between the towns of Abha and Muhayil Asir. An additional brief excursion to the same area was undertaken in October 2020. Rodents were identified morphologically with reference to the work of Harrison and Bates . Molecular confirmation was performed by amplification of a cytochrome b gene barcode using conventional polymerase chain reaction (PCR) with primers L14841 and H15149 . These rodent sequences were submitted to the Barcode of Life Data Systems (BOLD) (http://www.boldsystems.org) under project code SSS.
Morphological examination of ticks
Mammal carcasses were examined for ticks with the naked eye and then under a dissecting microscope. Ticks were removed with fine forceps, fixed in 70% ethanol and maintained at 4 °C prior to enumeration. Approximately 5% of specimens from each host were selected for morphological or molecular analysis, prioritising nymphs over larvae due to the low DNA yields and problems in identification associated with the latter. Semi-engorged immature stages selected for morphological examination were placed in distilled water for 10 min, transferred to a macerating solution (10% potassium hydroxide) and incubated at 37 °C for up to 10 min until the cuticle had cleared sufficiently to visualise key morphological features. The specimens were again placed in distilled water for 10 min and then dehydrated serially using 50%, 70% and 100% ethanol (10 min at each concentration). Finally, specimens were transferred to a glass slide with a drop of DPX mountant (VWR International), covered and examined using an Axio Imager M2 microscope with ZEN 2011 imaging software (Zeiss). Adult ticks (males only in this study, as females were fully engorged) were examined directly from 70% ethanol under a dissecting microscope without further processing. Morphological features of the ticks were compared with those described in keys and other taxonomic reference works for ticks, focusing on the Middle East, Southern Europe and North Africa [14, 21, 22, 32,33,34,35,36,37,38].
DNA extraction, PCR and sequencing
DNA extraction was performed with a DNeasy Blood & Tissue Kit (Qiagen) according to the manufacturer’s instructions. In the case of immature ticks, DNA was extracted from the whole specimen, whereas for adults, DNA extraction was performed on the anterior portion only to reduce carryover of the blood meal in engorged specimens. The amplification of fragments of three mitochondrial loci (cox1, 12S rRNA and 16S rRNA) was attempted for each specimen (Table 1) using previously published primers from Low and Prakash , Beati and Keirans , and Black and Piesman , respectively. Expected product sizes were 550 bp for cox1, 336 bp for 12S rRNA and 460 bp for 16S rRNA. The PCR assays were performed on a T1 Thermoblock thermocycler (Biometra) using BioMix Red reaction mix (Meridian Bioscience) in 20-μl volumes containing 5 μl DNA template. Following agarose gel electrophoresis, PCR products were purified using a QIAquick PCR Purification Kit (Qiagen) and sequenced in both directions by Eurofins Genomics.
Phylogenetic analysis and automatic barcode gap discovery
Reference sequences were selected based on Low and Prakash , Chandra et al.  and Kanduma et al.  for Rhipicephalus spp., and Hornok et al.  for Haemaphysalis spp. Additional closest reference nucleotide sequences displayed in the Basic Local Sequence Alignment Tool  were also included for phylogenetic tree construction. All sequences were preliminarily aligned using CLUSTAL X  and edited using BioEdit . Phylogenetic relationships were inferred using the neighbour-joining method in MEGA X . The neighbour-joining bootstrap values were estimated using 1000 replicates with Kimura’s two-parameter model of substitution (K2P distance), and bootstrap proportions of > 70% were considered well supported . Gaps and missing data were eliminated. Statistical congruence was calculated using a partition homogeneity test implemented in PAUP 4.0b10 . No significant differences were found among separate gene regions (P = 0.800); hence, cox1, 12S and 16S sequences were concatenated for further analyses. To assess the genetic divergence of taxa, uncorrected (p) pairwise genetic distances among species were estimated using PAUP 4.0b10 . The species boundary among tick taxa was assigned by automatic barcode gap discovery (ABGD) analysis performed on the web server using the Kimura (K80) TS/TV model. Entity recognition was based on the suggested partition at P = 0.01 .
Material obtained and examined
We obtained 75 small mammal hosts across the three sites, which belonged to six species (Table 2): the eastern spiny mouse (Acomys dimidiatus), king jird (Meriones rex), Yemeni mouse (Myomyscus yemeni), black rat (Rattus rattus), house mouse (Mus musculus) and desert hedgehog (Paraechinus aethiopicus). These were infested with a total of 733 ticks (608 Haemaphysalis spp. and 125 Rhipicephalus spp.), all of which were immature except for nine adults (seven males and two females) recovered from the hedgehog. The overall prevalence of tick infestation was 70.7%, with a mean abundance per host of 9.8. Most subsampled specimens (80–95%) from each host were prioritised for molecular analysis, including pathogen screening (to be reported separately), and we focused primarily on Rhipicephalus spp. due to its greater potential regional importance as a disease vector. All specimens subjected to PCR [Rhipicephalus spp. nymphs (n = 33), one pool of Rhipicephalus spp. larvae, and eight Haemaphysalis spp. nymphs] generated at least one mitochondrial gene sequence (Table 1). At least two specimens per life cycle stage of each tick genus were examined morphologically.
Rhipicephalus spp. nymphs displayed variation in the length and shape of the palps as well as the appearance of the scutum, which slightly overlapped coxa III in some individuals only (Fig. 1c, d). Nymphs exhibited a highly reduced external spur on coxa I, and the internal spur appeared vestigial (Fig. 1d). According to the works of Pegram et al. [35, 36] on the R. sanguineus group, these features of the spurs together with the ratio of length to width of the capitulum would position these specimens closer in morphology to R. camicasi than to R. turanicus or R. sanguineus sensu lato (s.l.). In addition, the adanal plates of the adult males (Fig. 1e) lacked the distinctly concave shape proximal to the anus reported by Nava et al.  in their re-description of R. sanguineus sensu stricto (s.s.).
The Haemaphysalis spp. nymphs displayed palps that were flared posteriorly (Fig. 2), which according to Hoogstraal et al.  is a feature of H. erinacei that distinguishes it from H. sulcata. However, the ventral spur on palp segment I (Fig. 2b) had a triangular profile unlike that of H. erinacei. Since Hoogstraal and Kaiser  and Hoogstraal et al.  also reported H. leachi from the Arabian Peninsula, we consulted the descriptions and re-descriptions of this species and the closely related H. elliptica from Africa [32, 37]. The posterior margin of the basis capituli in both of these species is convex, but in some of the specimens from 'Asir, it is straight (compare Fig. 2b, c).
Sequence analysis of Rhipicephalus spp.
At least one mitochondrial gene sequence was amplified and sequenced successfully from a total of 33 Rhipicephalus spp. adult or nymphal tick specimens and one pool of larvae, obtained from two villages and four species of small mammal host (Table 1). The Rhipicephalus spp. phylogeny based on cox1 indicated that the vast majority of nymphal specimens belonged to a single, novel clade with 90% bootstrap support; this was distinct from all other Rhipicephalus spp. included in the analysis (Fig. 3). The novel clade exhibited closest relationships with R. leporis, R. guilhoni and R. linnaei. In contrast, a single nymph (R25 from host A. dimidiatus in Wosanib) clustered with an adult specimen from the current study (H1_2 from host P. aethiopicus, also from Wosanib) and previously published sequences from “R. cf. camicasi” from Riyadh Province. The novel lineage was separated from other species by a minimum genetic distance of 2.2% (for R. leporis) to a maximum of 15.4% (for R. simus) (Additional file 1: Table S1). The ABGD analysis delimited 18 operational taxonomic units (OTUs) and supported the novel clade comprising most nymph specimens (OTU 1) as a distinct taxon (Fig. 3).
For 12S rRNA, the novel lineage was also resolved for all nymphs except R25 (> 90% bootstrap support). The clade differed from other members of the genus with lower genetic distances of 1.8% (for R. leporis) to 11.1% for the R. simus complex (including an unidentified Rhipicephalus sp. from Kenya; Additional file 1: Table S2). The ABGD analysis identified 16 OTUs, and although lower interspecific genetic distances were observed, the delimitation analysis demonstrated the novel lineage as a distinct OTU (Fig. 4). The “R. cf. camicasi” specimens from the previous study in Riyadh Province (obtained from camels and a dog) were split into three distinct OTUs, suggesting cryptic diversity in this species. One of these (from a camel) clustered with nymph specimen R25. Interestingly, the pool of six larvae (R29 from Wosanib) was placed in a unique OTU separated from all nymph specimens (Fig. 4). This was most closely related to members of the R. simus complex from Africa, especially R. praetextatus; indeed, the larval pool was not differentiated from the R. simus complex in the PAUP analysis (Additional file 1: Table S2).
In the case of 16S rRNA, the novel lineage (> 70% bootstrap support) was also distantly separated from other members of the genus with genetic distances ranging from 4.1% (for R. guilhoni) to 12.3% (for R. muhsamae) (Additional file 1: Table S3). A total of 15 OTUs were delimited, one of which was associated with the novel lineage (Fig. 5). Rhipicephalus cf. camicasi comprised two OTUs, populated by adult specimens from P. aethiopicus, four nymph specimens and the previously published sequences from specimens collected from camels in Riyadh Province. An incongruence was noted for one of the tick samples, nymph R9_7 from Alogl, which was classified in the novel lineage by cox1 and 12S rRNA genes but clustered with R. cf. camicasi OTU 4 by 16S rRNA (Fig. 5). The pool of larvae (R29) formed its own OTU (#12 in Fig. 5) that was most closely related to a sequence (OTU 13) from an unidentified Rhipicephalus sp. collected from a dog in Kenya (GenBank: MN266945).
Sufficient sequence data were obtained from 10 nymph specimens for a concatenated analysis of cox1, 12S rRNA and 16S rRNA genes alongside references for R. sanguineus s.s., R. linnaei, R. cf. camicasi, R. turanicus and R. simus. The novel clade comprised eight specimens and was distinct from all references (bootstrap support 86%), demonstrating closest affinity with R. linnaei (Fig. 6). In concordance with the single-gene trees, specimen R25 clustered with one of two R. cf. camicasi OTUs, whereas the incongruent specimen R9_7 formed its own OTU in proximity to R. linnaei (Fig. 6). As only short sequences (~ 200 bases) for 12S rRNA could be obtained from the two adult ticks from P. aethiopicus, they were excluded in the concatenated analysis. However, these short sequences exhibited 100% identity with the previously published R. cf. camicasi sequences from Riyadh Province (GenBank MH094506 and MH094507 from camel hosts).
Sequence analysis of Haemaphysalis spp.
The Haemaphysalis nymph samples collected in this study were resolved robustly into two lineages (100% bootstrap support) in the 16S rRNA phylogenetic tree (Fig. 7). While OTU 1 demonstrated a sister relationship with H. spinulosa from South Africa (genetic distance, 7.4%), OTU 4 showed closer relationships with H. muhsamae and H. elliptica, also from sub-Saharan Africa, with genetic distances of 6.8% and 8.2%, respectively (Additional file 1: Table S4). The species delimitation analysis split the Saudi specimens and references into a total of 15 OTUs, with the Saudi nymphs distinctly separated from all other species included in the analysis (Fig. 7). Notably, these two novel OTUs did not segregate by geographic location (Table 1), with OTU 1 containing specimens from both Alogl (M. musculus as host) and Al Ous' (A. dimidiatus as hosts).
In this study, we found widespread tick infestations represented by two genera feeding on small mammals in a relatively small region (approximately 170 km2) in the 'Asir Mountains. The tick abundance per host (9.8) was moderately high compared with previous studies of ticks on rodents in Saudi Arabia; for instance, in Riyadh Province, the mean abundance on gerbils was < 1, and in Ta'if (Makkah Province) it was 0.6–6.2 on gerbils and 1.4 on A. dimidiatus . However, in Ha'il Province in the north, a mean tick abundance of ~ 20 was recorded on R. rattus in the most heavily infested geographical sites, although the average abundance on A. dimidiatus was only 3.2 . Remarkably, the Rhipicephalus and Haemaphysalis ticks recovered from 'Asir not only were genetically diverse, comprising four and two OTUs, respectively, but all but one (R. camicasi) of these OTUs appeared to be novel. The strongest evidence for a previously unrecognised taxon was for Rhipicephalus OTU 1, which formed a distinct clade in the cox1, 12S rRNA, 16S rRNA and concatenated analyses. This clade was found on three species of rodent hosts trapped in agricultural areas surrounding the villages of Alogl and Wosanib. It was most closely related to R. leporis, R. guilhoni and the “tropical lineage” of R. sanguineus s.l. (recently identified as R. linnaei ). Due to the limited number of sequenced mitochondrial markers available for R. leporis and R. guilhoni, we were only able to include R. linnaei and more distantly related Rhipicephalus spp. in the concatenated phylogeny, but this analysis clearly separated the novel OTU 1 from R. linnaei.
Prior phylogenetic analyses have sometimes assigned R. leporis and R. guilhoni to the same clade as R. linnaei, along with R. camicasi, depending on the loci included [14, 15, 47,48,49]. The taxonomy and biogeography of the R. sanguineus group are notoriously complex due to their morphological similarity and the tendency for different species or clades to be spread worldwide on domestic hosts. Estrada-Peña et al.  consider R. guilhoni and R. camicasi as tropical species that have invaded Palearctic regions, whereas R. leporis appears to be a Palearctic species that has been introduced into sub-Saharan Africa . There are few molecular data available for R. camicasi, but the sequences provided by Chandra et al.  for “R. cf. camicasi” from Riyadh Province are clearly distinct from available references for other Rhipicephalus spp. and clustered with a small proportion of our nymph specimens from rodents. To add further to the complexity, R. camicasi from Saudi Arabia did not form a single OTU in our analyses, including in the concatenated phylogeny.
Rhipicephalus camicasi was originally described from Northeast Africa in 1976 . It was not included in the tick fauna of Saudi Arabia by Hoogstraal et al. , who listed only two native Rhipicephalus spp. (R. sanguineus s.l. and R. turanicus), excluding the subgenus Boophilus. However, they noted the presence of unidentified Rhipicephalus spp. on numerous mammalian hosts, including A. dimidiatus, M. rex and M. musculus. Subsequently, Pegram et al.  stated that R. camicasi could be found on livestock (ruminants, camels and donkeys) in Yemen and Saudi Arabia without details of specific locations. More recently, R. camicasi has been reported from sheep in Makkah Province  and from camels and dogs in Riyadh Province [25, 26], as well as from A. dimidiatus (as nymphs and larvae) in Ta'if . To the best of our knowledge, R. camicasi has not been reported from a hedgehog host previously worldwide . Our incidental finding of R. camicasi on a single P. aethiopicus in this study should be followed by a targeted survey to determine whether this common and widespread host acts as a vehicle or reservoir to maintain R. camicasi populations nationwide.
Very few studies have attempted to identify ticks from small mammal hosts from Saudi Arabia or Yemen previously. However, the classic wild mammal survey of Yemen (which borders 'Asir) by Sanborn and Hoogstraal  reported R. simus, R. sanguineus s.l. and Ornithodoros sp. from M. musculus; H. leachi and R. simus from A. dimidiatus; and R. simus and “Ixodes sp. nov.” from M. rex, among a wide range of other hosts examined. Similar host–ectoparasite relationships were recorded by Hoogstraal et al.  for Saudi Arabia, with the addition of immature Hyalomma spp. observed on all three rodent species. Our finding of Rhipicephalus larvae on M. rex that appeared to be closely related to the R. simus complex supports these early observations of Hoogstraal regarding the introduction of African Rhipicephalus spp. into the Arabian Peninsula. Asiry and Fetoh  described R. turanicus infestations on A. dimidiatus, alongside R. sanguineus s.l. and R. turanicus feeding on R. rattus, from Ha'il Province. Notably, the most recent prior survey by Harrison et al.  echoed the work of Hoogstraal et al.  in reporting the presence of an unidentified immature Rhipicephalus sp. on rodents in Riyadh and Ta'if. It was most common on M. rex in Ta'if, but was also found on Meriones lybicus in Riyadh and in smaller numbers on Gerbillus nanus in both locations. Only a single specimen was found on A. dimidiatus (in Ta'if), a host species on which it was apparently outcompeted by R. camicasi (see above). However, no morphological description (in particular, how the specimens were differentiated from R. camicasi) or molecular barcode was provided for this unidentified Rhipicephalus sp. Overall, these studies from Arabia highlight distinct differences compared with the wider Middle East, as a recent systematic review reported that Hyalomma rhipicephaloides and Ixodes eldaricus were the most prevalent ticks found on rodents in the whole region, representing 69.7% and 15.7% of ticks identified, respectively .
The only native Haemaphysalis spp. recorded from Saudi Arabia in Hoogstraal et al.  were H. erinacei and H. sulcata; while in Yemen, H. leachi (presumably introduced from Africa) was reported on A. dimidiatus [21, 51]. Prior to the emergence of severe fever with thrombocytopenia syndrome virus and the global spread of its vector, H. longicornis, molecular analyses of the genus Haemaphysalis had been relatively limited . However, sufficient data are available to conclude that neither the morphology nor the 16S rRNA sequences of our Haemaphysalis spp. specimens are fully compatible with species previously recorded from Arabia. The two distinct OTUs we identified exhibited closest relationships with African Haemaphysalis spp. (H. spinulosa, H. muhsamae and H. elliptica) that primarily parasitise carnivores or erinaceids in the adult stage and rodents as immature stages [37, 54]. The previous surveys of rodent ticks conducted in Arabia (see above) suggest that Haemaphysalis spp. are restricted (or at least more abundant) in Yemen and southern Saudi Arabia compared with more northern regions. Whether the novel Haemaphysalis OTUs represent undescribed species native to the southern Arabian Peninsula will require further investigations, including locating adult specimens for comprehensive morphological and molecular analyses.
This first molecular analysis of ticks collected from rodents in the Arabian Peninsula raises many questions about the evolution and distribution of Rhipicephalus spp. and Haemaphysalis spp. in this understudied region. For instance, the taxonomic status and native geographical range of R. camicasi is still poorly defined, especially with respect to its relationship with R. linnaei. As highlighted by Hekimoglu et al. , Asia Minor and the Middle East constitute a bridge between Europe and Africa in the evolutionary history of the R. sanguineus group, in which the role of R. camicasi remains enigmatic. A limitation of our study was that in order to maximise DNA yields for multiple PCR assays, a portion of each specimen was not retained as a voucher  prior to DNA extraction. Hence, it is not clear whether R. camicasi and Rhipicephalus OTU 1 are morphologically distinct in the immature stages, which would be suggested by the work of Harrison et al. —if OTU 1 is indeed the species they recorded from Riyadh and Ta'if.
The strongest evidence that OTU 1 constitutes a distinct species is that the ABGD analysis consistently binned known Rhipicephalus spp. at each locus into separate OTUs, while OTU 1 was delimited with moderate to strong phylogenetic support (bootstrap values ≥ 70%; ) for each locus but especially in the concatenated tree. However, the genetic distance of OTU 1 from R. leporis was modest (~ 2%) at the cox1 and 12S rRNA loci; while as noted above, a lack of 16S rRNA data for this species prevented its inclusion in the concatenated analysis. In a previous tick barcoding study across multiple genera using neighbour-joining analysis only, optimal species delimitation boundaries were considered to be 5.3% for 16S rRNA and 6.1% for cox1, although few Rhipicephalus spp. were included in this work . The formal characterisation of Rhipicephalus OTU 1 will require sampling of adult stages from the environment [or host(s)—which of course remain unknown currently], followed by detailed morphological and molecular comparisons with closely related species and (ideally) laboratory experiments to determine reproductive compatibility. This would provide substantive evidence that OTU 1 is a novel species or a subspecies, rather than just a divergent mitochondrial haplotype. Indeed, although mitochondrial versus nuclear marker-based phylogenies for ticks are generally congruent [3, 20], nuclear–mitochondrial discordance has been observed within tick species previously . Moreover, the incongruent results between mitochondrial loci for specimen R9_7 could indicate hybridisation between Rhipicephalus OTU 1 and R. camicasi.
In a small region of the 'Asir Mountains in south-western Saudi Arabia, small mammals were found to be infested with Rhipicephalus spp. and Haemaphysalis spp. ticks that formed four and two clades, respectively, by the ABGD method. In addition to two clades of R. camicasi-like adult and nymphal ticks and one clade of R. simus-like larvae, a novel OTU composed of Rhipicephalus nymphs was found infesting three species of rodent hosts. It was related to, but distinct from, R. leporis, R. guilhoni and R. linnaei. Taken together, our findings indicate that a hotspot of tick diversity may exist in the 'Asir Mountains that deserves further faunistic, ecological and genetic investigations. This conclusion is supported by the fact that most previous studies of tick genetic diversity have identified population structuring at the level of whole countries or continents rather than at a highly localised level [58,59,60,61,62]. However, it remains provisional until additional molecular studies are conducted on ticks derived from small mammals in this part of the world. Prior ectoparasite sampling from rodents trapped in other regions of Saudi Arabia suggest that Rhipicephalus OTU 1 might constitute a widespread novel species, and future studies should focus on locating adult specimens to permit a formal description of the taxon.
Availability of data and materials
Sequence data are available from GenBank with identifiers MW742686-MW742711 for cox1, MW756110-MW756125 for 12S rRNA, and MW763030-MW763059 for 16S rRNA.
Automatic barcode gap discovery
Operational taxonomic unit
Guglielmone A, Robbins RG, Apanaskevich DA, Petney TN, Estrada-Pena A, Horak IG, et al. The Argasidae, Ixodidae and Nuttalliellidae (Acari: Ixodida) of the world: a list of valid species names. Zootaxa. 2010;2528:1–28.
Murgia MV, Bell-Sakyi L, de la Fuente J, Kurtti TJ, Makepeace BL, Mans B, et al. Meeting the challenge of tick-borne disease control: A proposal for 1000 Ixodes genomes. Ticks Tick Borne Dis. 2019;10:213–8.
Mans BJ, Featherston J, Kvas M, Pillay KA, de Klerk DG, Pienaar R, et al. Argasid and ixodid systematics: Implications for soft tick evolution and systematics, with a new argasid species list. Ticks Tick Borne Dis. 2019;10:219–40.
Kelava S, Mans BJ, Shao R, Moustafa MA, Matsuno K, Takano A, et al. Phylogenies from mitochondrial genomes of 120 species of ticks: Insights into the evolution of the families of ticks and of the genus Amblyomma. Ticks Tick Borne Dis. 2021;12:101577.
Kwak ML, Beveridge I, Koehler AV, Malipatil M, Gasser RB, Jabbar A. Phylogenetic analysis of the Australasian paralysis ticks and their relatives (Ixodidae: Ixodes: Sternalixodes). Parasit Vectors. 2017;10:122.
Black WC, Piesman J. Phylogeny of hard- and soft-tick taxa (Acari: Ixodida) based on mitochondrial 16S rDNA sequences. Proc Natl Acad Sci USA. 1994;91:10034–8.
Hornok S, Wang Y, Otranto D, Keskin A, Lia RP, Kontschan J, et al. Phylogenetic analysis of Haemaphysalis erinacei Pavesi, 1884 (Acari: Ixodidae) from China, Turkey. Italy and Romania Parasit Vectors. 2016;9:643.
McCann KM, Grant WN, Spratt DM, Hedtke SM. Cryptic species diversity in ticks that transmit disease in Australia. Int J Parasitol Parasites Wildl. 2019;10:125–31.
Burger TD, Shao R, Barker SC. Phylogenetic analysis of mitochondrial genome sequences indicates that the cattle tick, Rhipicephalus (Boophilus) microplus, contains a cryptic species. Mol Phylogenet Evol. 2014;76:241–53.
Low VL, Prakash BK. First genetic characterization of the brown dog tick Rhipicephalus sanguineus sensu lato in Peninsular Malaysia. Exp Appl Acarol. 2018;75:299–307.
Roy BC, Estrada-Pena A, Krucken J, Rehman A, Nijhof AM. Morphological and phylogenetic analyses of Rhipicephalus microplus ticks from Bangladesh, Pakistan and Myanmar. Ticks Tick Borne Dis. 2018;9:1069–79.
Low VL, Tay ST, Kho KL, Koh FX, Tan TK, Lim YA, et al. Molecular characterisation of the tick Rhipicephalus microplus in Malaysia: new insights into the cryptic diversity and distinct genetic assemblages throughout the world. Parasit Vectors. 2015;8:341.
Dantas-Torres F, Latrofa MS, Annoscia G, Giannelli A, Parisi A, Otranto D. Morphological and genetic diversity of Rhipicephalus sanguineus sensu lato from the New and Old Worlds. Parasit Vectors. 2013;6:213.
Nava S, Beati L, Venzal JM, Labruna MB, Szabo MPJ, Petney T, et al. Rhipicephalus sanguineus (Latreille, 1806): Neotype designation, morphological re-description of all parasitic stages and molecular characterization. Ticks Tick Borne Dis. 2018;9:1573–85.
Chitimia-Dobler L, Langguth J, Pfeffer M, Kattner S, Kupper T, Friese D, et al. Genetic analysis of Rhipicephalus sanguineus sensu lato ticks parasites of dogs in Africa north of the Sahara based on mitochondrial DNA sequences. Vet Parasitol. 2017;239:1–6.
Burger TD, Shao R, Barker SC. Phylogenetic analysis of the mitochondrial genomes and nuclear rRNA genes of ticks reveals a deep phylogenetic structure within the genus Haemaphysalis and further elucidates the polyphyly of the genus Amblyomma with respect to Amblyomma sphenodonti and Amblyomma elaphense. Ticks Tick Borne Dis. 2013;4:265–74.
Puillandre N, Lambert A, Brouillet S, Achaz G. ABGD, Automatic Barcode Gap Discovery for primary species delimitation. Mol Ecol. 2012;21:1864–77.
Lv J, Wu S, Zhang Y, Zhang T, Feng C, Jia G, et al. Development of a DNA barcoding system for the Ixodida (Acari: Ixodida). Mitochondrial DNA. 2014;25:142–9.
Evans ML, Egan S, Irwin PJ, Oskam CL. Automatic Barcode Gap Discovery reveals large COI intraspecific divergence in Australian Ixodidae. Zootaxa. 2019;4656:393–6.
Kanduma EG, Bishop RP, Githaka NW, Skilton RA, Heyne H, Mwacharo JM. Mitochondrial and nuclear multilocus phylogeny of Rhipicephalus ticks from Kenya. Mol Phylogenet Evol. 2019;140:106579.
Hoogstraal H, Kaiser MN. Ticks (Ixodoidea) of Arabia, with special reference to the Yemen. v.39:no.28 (1959). [Chicago]: Chicago Natural History Museum.
Hoogstraal H, Wassef HY, Büttiker W. Ticks (Acarina) of Saudi Arabia, fam Argasidae, Ixodidae. Fauna Saudi Arabia. 1981;3:25–110.
Alanazi AD, Alouffi AS, Alshahrani MY, Alyousif MS, Abdullah H, Allam AM, et al. A report on tick burden and molecular detection of tick-borne pathogens in cattle blood samples collected from four regions in Saudi Arabia. Ticks Tick Borne Dis. 2021;12:101652.
Alanazi AD, Nguyen VL, Alyousif MS, Manoj RRS, Alouffi AS, Donato R, et al. Ticks and associated pathogens in camels (Camelus dromedarius) from Riyadh Province, Saudi Arabia. Parasit Vectors. 2020;13:110.
Chandra S, Smith K, Alanazi AD, Alyousif MS, Emery D, Slapeta J. Rhipicephalus sanguineus sensu lato from dogs and dromedary camels in Riyadh, Saudi Arabia: low prevalence of vector-borne pathogens in dogs detected using multiplexed tandem PCR panel. Folia Parasitol (Praha). 2019;66:007.
Alanazi AD, Al-Mohammed HI, Alyousif MS, Said AE, Salim B, Abdel-Shafy S, et al. Species diversity and seasonal distribution of hard ticks (Acari: Ixodidae) infesting mammalian hosts in various districts of Riyadh Province, Saudi Arabia. J Med Entomol. 2019;56:1027–32.
Harrison A, Robb GN, Alagaili AN, Hastriter MW, Apanaskevich DA, Ueckermann EA, et al. Ectoparasite fauna of rodents collected from two wildlife research centres in Saudi Arabia with discussion on the implications for disease transmission. Acta Trop. 2015;147:1–5.
Asiry KA, Fetoh BS. Occurrence of ectoparasitic arthropods associated with rodents in Hail region northern Saudi Arabia. Environ Sci Pollut Res Int. 2014;21:10120–8.
Stekolnikov AA, Al-Ghamdi SQ, Alagaili AN, Makepeace BL. First data on chigger mites (Acariformes: Trombiculidae) of Saudi Arabia, with a description of four new species. Syst Appl Acarol. 2019;24:1937–63.
Harrison DL, Bates PJJ. The mammals of Arabia. 2nd ed. Sevenoaks, Kent: Harrison Zoological Museum; 1991.
Kocher TD, Thomas WK, Meyer A, Edwards SV, Paabo S, Villablanca FX, et al. Dynamics of mitochondrial DNA evolution in animals: amplification and sequencing with conserved primers. Proc Natl Acad Sci USA. 1989;86:6196–200.
Hoogstraal H. Notes on African Haemaphysalis ticks. IV. Description of Egyptian populations of the yellow dog-tick, H. leachii leachii (Audouin, 1827) (Ixodoidea, Ixodidae). J Parasitol. 1958;44:548–58.
Morel PC, Mouchet J, Rodhain F. Description of Rhipicephalus camicasi n. sp. (Acaridae, Ixodida) of subdesert steppes of the plain of Afar. Rev Elev Med Vet Pays Trop. 1976;29:337–40.
Pegram RG, Walker JB, Clifford CM, Keirans JE. Comparison of populations of the Rhipicephalus simus group: R. simus, R. praetextatus, and R. muhsamae (Acari: Ixodidae). J Med Entomol. 1987;24:666–82.
Pegram RG, Clifford CM, Walker JB, Keirans JE. Clarification of the Rhipicephalus sanguineus group (Acari, Ixodoidea, Ixodidae). 1. Rhipicephalus sulcatus Neumann, 1908 and Rhipicephalus turanicus Pomerantsev, 1936. Syst Parasitol. 1987;10:3–26.
Pegram RG, Keirans JE, Clifford CM, Walker JB. Clarification of the Rhipicephalus sanguineus group (Acari, Ixodoidea, Ixodidae). 2. Rhipicephalus sanguineus (Latreille, 1806) and related species. Syst Parasitol. 1987;10:27–44.
Apanaskevich DA, Horak IG, Camicas JL. Redescription of Haemaphysalis (Rhipistoma) elliptica (Koch, 1844), an old taxon of the Haemaphysalis (Rhipistoma) leachi group from East and southern Africa, and of Haemaphysalis (Rhipistoma) leachi (Audouin, 1826) (Ixodida, Ixodidae). Onderstepoort J Vet Res. 2007;74:181–208.
Estrada-Pena A, Mihalca AD, Petney TN. Ticks of Europe and North Africa: a guide to species identification. Cham: Springer Nature; 2017.
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.
Altschul SF, Gish W, Miller W, Myers EW, Lipman DJ. Basic local alignment search tool. J Mol Biol. 1990;215:403–10.
Thompson JD, Gibson TJ, Plewniak F, Jeanmougin F, Higgins DG. The CLUSTAL_X windows interface: flexible strategies for multiple sequence alignment aided by quality analysis tools. Nucleic Acids Res. 1997;25:4876–82.
Hall TA. BioEdit: a user-friendly biological sequence alignment editor and analysis program for Windows 95/98/NT. Nucleic Acids Symp Ser. 1999;41:95–8.
Kumar S, Stecher G, Li M, Knyaz C, Tamura K. MEGA X: Molecular Evolutionary Genetics Analysis across computing platforms. Mol Biol Evol. 2018;35:1547–9.
Hillis DM, Bull JJ. An empirical test of bootstrapping as a method for assessing confidence in phylogenetic analysis. Syst Biol. 1993;42:182–92.
Swofford D. PAUP*. Phylogenetic Analysis Using Parsimony (*and Other Methods). Version 4.0b10. vol. Version 4.0; 2002.
Slapeta J, Chandra S, Halliday B. The, “tropical lineage” of the brown dog tick Rhipicephalus sanguineus sensu lato identified as Rhipicephalus linnaei (Audouin, 1826). Int J Parasitol. 2021;51:431–6.
Hornok S, Sandor AD, Tomanovic S, Beck R, D’Amico G, Kontschan J, et al. East and west separation of Rhipicephalus sanguineus mitochondrial lineages in the Mediterranean Basin. Parasit Vectors. 2017;10:39.
Zemtsova GE, Apanaskevich DA, Reeves WK, Hahn M, Snellgrove A, Levin ML. Phylogeography of Rhipicephalus sanguineus sensu lato and its relationships with climatic factors. Exp Appl Acarol. 2016;69:191–203.
Bakkes DK, Chitimia-Dobler L, Matloa D, Oosthuysen M, Mumcuoglu KY, Mans BJ, et al. Integrative taxonomy and species delimitation of Rhipicephalus turanicus (Acari: Ixodida: Ixodidae). Int J Parasitol. 2020;50:577–94.
El-Azazy OM, Scrimgeour EM. Crimean-Congo haemorrhagic fever virus infection in the western province of Saudi Arabia. Trans R Soc Trop Med Hyg. 1997;91:275–8.
Sanborn CC, Hoogstraal H. Some mammals of Yemen and their ectoparasites. Chicago: Chicago Natural History Museum; 1953.
Islam MM, Farag E, Eltom K, Hassan MM, Bansal D, Schaffner F, et al. Rodent ectoparasites in the Middle East: A systematic review and meta-analysis. Pathogens. 2021;10:139.
Thompson AT, Dominguez K, Cleveland CA, Dergousoff SJ, Doi K, Falco RC, et al. Molecular characterization of Haemaphysalis species and a molecular genetic key for the identification of Haemaphysalis of North America. Front Vet Sci. 2020;7:141.
Tomlinson JA, Horak IG, Apanaskevich DA. Identity of Haemaphysalis (Rhipistoma) muhsamae Santos Dias, 1954 (Acari: Ixodidae) and H. (R.) subterra Hoogstraal, El Kammah & Camicas, 1992, parasites of carnivores and rodents in eastern and southern Africa. Syst Parasitol. 2018;2018(95):673–91.
Hekimoglu O, Saglam IK, Ozer N, Estrada-Pena A. New molecular data shed light on the global phylogeny and species limits of the Rhipicephalus sanguineus complex. Ticks Tick Borne Dis. 2016;7:798–807.
Scott JD, Foley JE, Young MR, Durden LA. First report of a blacklegged tick, Ixodes scapularis Say (Acari: Ixodidae), parasitizing a raptor in Canada. Syst Appl Acarol. 2017;22:208–16.
Leo SS, Pybus MJ, Sperling FA. Deep mitochondrial DNA lineage divergences within Alberta populations of Dermacentor albipictus (Acari: Ixodidae) do not indicate distinct species. J Med Entomol. 2010;47:565–74.
Al-Khafaji AM, Clegg SR, Pinder AC, Luu L, Hansford KM, Seelig F, et al. Multi-locus sequence typing of Ixodes ricinus and its symbiont Candidatus Midichloria mitochondrii across Europe reveals evidence of local co-cladogenesis in Scotland. Ticks Tick Borne Dis. 2019;10:52–62.
Jia N, Wang J, Shi W, Du L, Sun Y, Zhan W, et al. Large-scale comparative analyses of tick genomes elucidate their genetic diversity and vector capacities. Cell. 2020;182:1328–40.
Poli P, Lenoir J, Plantard O, Ehrmann S, Roed KH, Leinaas HP, et al. Strong genetic structure among populations of the tick Ixodes ricinus across its range. Ticks Tick Borne Dis. 2020;11:101509.
Xu G, Wielstra B, Rich SM. Northern and southern blacklegged (deer) ticks are genetically distinct with different histories and Lyme spirochete infection rates. Sci Rep. 2020;10:10289.
Lado P, Smith ML, Carstens BC, Klompen H. Population genetic structure and demographic history of the lone star tick, Amblyomma americanum (Ixodida: Ixodidae): New evidence supporting old records. Mol Ecol. 2020;29:2810–23.
Underwood W, Anthony R, Cartner S, Corey D, Grandin T, Greenacre CB, et al. AVMA guidelines for the euthanasia of animals. American Veterinary Medical Association: Schaumburg; 2013.
Charbonneau R, Niel L, Olfert E, von Keyserlingk M, Griffin C. CCAC guidelines on: Euthanasia of animals used in science. Ottawa: Canadian Council on Animal Care; 2010.
This work was funded by doctoral scholarships from Al Baha University (to SQA) and from Imam Muhammad Ibn Saud Islamic University (to HAA). The Vice Deanship of Research Chairs, Deanship of Scientific Research of the King Saud University, Kingdom of Saudi Arabia, provided support to ANA. The funders had no roles in the design of the study; collection, analysis, and interpretation of data; or in writing the manuscript.
Ethics approval and consent to participate
Permission for rodent trapping and euthanasia was received from the Saudi Wildlife Authority by co-author ANA before fieldwork commenced. Ethical approval from the University of Liverpool’s Animal Welfare and Ethics Review Board was also obtained. Rodents were euthanised by inhaled anaesthetic overdose according to guidelines published by the American Veterinary Medical Association Council on Research  and the Canadian Council on Animal Care . A single individual of the desert hedgehog (Paraechinus aethiopicus) was released live after removal of ticks.
Consent for publication
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
cox1 interspecific genetic distances among Rhipicephalus taxa. Table S2. 12S rRNA interspecific genetic distances among Rhipicephalus taxa. Table S3. 16S rRNA interspecific genetic distances among Rhipicephalus taxa. Table S4. 16S rRNA interspecific genetic distances among Haemaphysalis taxa.
About this article
Cite this article
Alghamdi, S.Q., Low, V.L., Alkathiry, H.A. et al. Automatic barcode gap discovery reveals diverse clades of Rhipicephalus spp. and Haemaphysalis spp. ticks from small mammals in 'Asir, Saudi Arabia. Parasites Vectors 14, 541 (2021). https://doi.org/10.1186/s13071-021-05049-x
- Molecular barcoding
- Meriones rex
- Acomys dimidiatus
- Brown dog tick