In order to successfully control tick-borne diseases, it is very important to develop good identification techniques, which enable successful separation of the different species of ticks. The RAPD-PCR is a molecular technique, which was used before the development of DNA sequencing. It is a good molecular technique, which is utilized for many years by taxonomists successfully and reliably   . However, with the advent of DNA sequencing new molecular identification techniques were developed by using segments of certain genes to provide good species identification. In a study  , researchers developed a DNA barcoding system for the identification of Ixodid ticks using segment of the 16S and COI genes. The identifications were successful, however it was mentioned that there were serious deficiencies in the information of 16S and COI of some species of ticks, and thus additional research was needed to resolve this problem. In another study, medically important ticks were identified using a 658 base pair DNA barcode sequence from the COI gene . Additionally, a genetic cluster analysis was performed on the barcode sequences and 26 Ixodid species were included. It was reported that DNA barcoding was a very good tool that clinicians could utilize to get the correct identification of tick specimens. To date, the use of a COI gene fragment is the standard marker for DNA barcoding. However, for COI, there are some limitations in its use with some species, and for some taxa . In order to study their relative effectiveness, the COI, 16S ribosomal DNA (rDNA), nuclear ribosomal internal transcribed spacer 2 (ITS2) and 12S rDNA were compared for tick species identification. The results showed that the intra-specific divergence of each marker was much lower than the inter-specific divergence. In addition, the aforementioned study indicated that COI was not significantly better than the other markers in terms of its rate of correct sequence identification. Nonetheless, COI should remain as the first choice for tick species identification. Moreover, a study reported that the DNA barcoding was successfully used in tick identification and that the K2P distances in most interspecific divergences exceeded 8%, while intraspecific distances were usually lower than 2% . However, the study indicated that more research was required to reveal ambiguous species delimitation in some problematic genera. But is should be mentioned that there is a chance for some misidentification when using DNA-based techniques. A study was conducted as a comparative test of identification of ticks occurring in Western Europe and Northern Africa. There was a sizable misidentification rate and thus the recommendation was that a combination of certain genes may adequately identify the target species of ticks .
Ticks are vectors to many pathogens and parasites, which can cause serious diseases to human, domestic animal and wildlife  . In different parts of the world such as North Africa, Southern Europe, Middle East, Central Asia and China Hyalomma (Acari: Ixodidae) ticks are common and may transmit Crimean-Congo hemorrhagic fever (CCHF), rickettsiae and tropical theileriosis  . One of the potentially important parasites to domestic animal health in the Middle East is Theileria annulata because it infects cattle and is transmitted by Hyalomma spp.   . Ornithodoros savignyi collected from camels in Egypt was reported to harbor Rickettsia-like microorganisms . Anaplasma marginale, Coxiella burnetii, Rickettsia aeschlim were reported from Egypt, Somalia and Nigeria from several Hyalomma spp. including H. dromedarii Koch, 1844  . Seventeen species of the genus Rickettsia are categorized within the Spotted Fever Group (SFG) rickettsiae  and several have been documented in the Middle East, although some reports are considered questionable . In the UAE, a spotted fever group Rickettsia sp. and Theileria annulata were recorded for the first time in H. dromedarii ticks in the country . Also, a Coxiella-like endosymbiont was reported in argasid ticks (Ornithodoros muesebecki) from a Socotra Cormorant colony in Umm Al Quwain, United Arab Emirates .
The aim of this study was to investigate the intraspecific genetic diversity of four H. dromedarii tick populations from the UAE.
2. Materials and Methods
2.1. Tick Collection
Camel ticks, H. dromedarii, Figure 1 were collected manually from the infested animals. Collected ticks were kept in 50-mL plastic tubes inside a −80˚C freezer. Ticks were collected from four locations: Falaj Al-Muaalla, Bayata, Al-Wagan, and Samha. Falaj Al-Mualla and Bayata are located in the Emirate of Umm Al-Qaiwain while Al-Wagan and Al-Samha are in the Emirate of Abu Dhabi. Twelve ticks were collected per location.
2.2. DNA Extraction and PCR
Genomic DNA of camel ticks was extracted from tick legs using a Maxwell 16 Tissue DNA purification kit (Promega, USA) following the manufacturer’s protocol. For the mitochondrial cytochrome c oxidase subunit I gene (COI) the following primer pair was used FishF1:
FishR1: 5’-TAGACTTCTGGGTGGCCAAAGAATCA-3’. The template DNA was amplified in a 25 μl reaction mixture containing 50 ng of mite DNA, 10 pmol of each primer pair, and 12.5 μl 2x PCR Master Mix (Promega, USA). Reaction mixtures were preheated at 95˚C for 5 min. Amplifications were carried out for 35 cycles (95˚C for 30 s, 58˚C for 30 s, and 72˚C for 1 min (final extension cycle at 72˚C for 5 min) in a Swift MaxPro thermocycler (ESCO, Singapore). For the RAPD-PCR assay, a set of 13 decamer primers (Operon Technologies, USA) were used. The sequences of these primers are presented in Table 1. These primers were selected as they generated consistent and reproducible RAPD banding patterns on the agarose gel. HotStarTaq™ Master Mix kit (Qiagen, USA) was used for PCR reactions. The reaction mixture included 20 ng template DNA, 12.5 μL HotStarTaq™ Master Mix, and 50 pM primer in a total volume of 25 μL. The program of thermal cycling was as follows: initial activation step at 95˚C for 5 min followed by 40 cycles of 1 min at 94˚C, 1 min at 36˚C and 2 min at 72˚C, with a final extension step at 72˚C for 5 min. To detect contamination, a negative
Figure 1. H. dromedarii adult female (a) and male (b). Ticks shriveled slightly after being stored in 70% ethanol.
Table 1. Random primers and number of generated bands of randomly amplified polymorphic DNA polymerase chain reaction (RAPD-PCR) in camel ticks.
control, in which all PCR components except DNA template were included, was run in every experiment. In order to insure reproducibility of the individual bands, PCR reactions on the extracted DNA were conducted at least three times. Amplification products were resolved by electrophoresis in a 1.5% agarose gel, which was stained with ethidium bromide and photographed by UVDI gel documentation system (Major Science, Taiwan).
2.3. PCR Amplicon Sequencing and Analysis
In the RAPD-PCR, the banding pattern produced by each primer on the agarose gel was scored as one or zero (one for the presence of a band and zero for the absence of the band). Band scores were analyzed using cluster analysis in SPSS statistical package (IBM SPSS Statistics for Windows). The PCR amplicons were purified by a PCR purification kit (Qiagen, Germany) following the manufacturer’s protocol. The purified PCR amplicons were sequenced (Sanger sequencing) by Source Bioscience (Cambridge, UK). Bases with too low confidence calling levels were trimmed from the beginning and end of sequence. Sequences were submitted in the GenBank and received accession numbers. They were analyzed by NCBI BLAST (http://blast.ncbi.nlm.nih.gov/Blast.cgi) analysis tool. In addition, sequence alignments and tree generation were conducted in MEGA6 . The percentage of replicate trees in which the associated taxa clustered together in the bootstrap test (1000 replicates) was shown next to the branches. Genetic distances were computed using the Kimura 2-parameter method and were in the units of the number of base substitutions per site.
The analysis of the banding patterns produced by the random primers revealed that the four H. dromedarii populations were grouped into four clusters Figure 2. In this grouping of cluster analysis, the Bayata, Al-Wagan, and Falaj Al-Mualla were closer to each other compared to Al-Samha.
3.2. COI Gene
The sequences of the PCR amplicons (≈600 bp) of the mitochondrial COI gene were submitted in the GenBanK and received accession numbers: Falaj Al-Mualla: MG188798, Bayata: MG188799, and Al-Samha: MG188800. The sequence of Al-Wagan KJ022631 was submitted already in the GenBank in a previous study .
The sequences of Falaj Al-Mualla, Al-Wagan, and Bayata were 100% similar to each other, while the one of Al-Samha was different from them in two nucleotides. The intra-specific divergence value was 0.002. In the maximum-likelihood Figure 3 and neighbor-joining Figure 4 trees the four sequences appeared in cluster of H. dromedarii sequences in which the sequence of Al-Samha showed up with some distance from the other three sequences. There was 99% - 100% similarity and 99% - 100% sequence coverage between the four sequences of the UAE and other H. dromedarii sequences in the GenBank such as KT920181.1, AJ437061.1, and KY548842.1. The inter-specific divergence value was 0.091.
The RAPD-PCR is a little bit an older molecular technique compared to Sanger sequencing, yet it is a good tool in studying genetic diversity. In a study  scientists reviewed more than 9000 papers that used RAPD and reported that
Figure 2. Dendogram of a cluster analysis showing three tick populations with higher DNA similarity versus a fourth population with less similarity (SPSS statistical package, IBM SPSS Statistics for Windows).
Figure 3. Maximum-likelihood tree constructed from nucleotide sequences (of the COI) of Hyalomma sp. based on Kimura 2-parameter. The Bootstrap branch support values were assessed with 1000 replicates. Falaj Al-Mualla (MG188798), Bayata (MG188799), Al-Samha (MG188800), and Al-Wagan (KJ022631) are four H. dromedarii sequences from the UAE.
Figure 4. Neighbor-Joinig tree constructed from nucleotide sequences (of the COI) of Hyalomma sp. based on Kimura 2-parameter. The Bootstrap branch support values were assessed with 1000 replicates. Falaj Al-Mualla (MG188798), Bayata (MG188799), Al-Samha (MG188800), and Al-Wagan (KJ022631) are four H. dromedarii sequences from the UAE.
with proper optimization, RAPD is a reliable, sensitive and reproducible assay and has the potential to detect a wide range of DNA damage as well as mutations. The current study supported this finding and showed that the results of RAPD-PCR were in agreement with the results of the mitochondrial cytochrome c oxidase subunit I (COI) gene. Accordingly, three tick populations (Falaj Al-Mualla, Bayata, and Al-Wagan) were very similar, while Al-Samha had a small genetic difference as represented by the dendogram.
The use of DNA sequencing of the COI gene is a very reliable identification tool in animal species. The current study did not show high levels of genetic diversity among the H. dromedarii ticks collected from the four locations in the UAE. In part, this could be explained by the fact that the sampling was limited to a few areas of UAE in a relatively short period of time. Therefore, future studies should investigate more locations and screen larger numbers of ticks in order to possibly find bigger genetic differences. In this study, the four tick populations showed no intraspecific sequence difference except in one population (Al-Samha) there were two different nucleotides. Moreover, there was 99% - 100% similarity between the four sequences of and other H. dromedarii sequences in the GenBank. The maximum-likelihood and neighbor-joining trees produced almost similar tree topologies and the branches were supported by high bootstrp values. The UAE tick sequences appeared in clusters of H. dromedarii on the trees. Other studies reported low genetic diversity when using COI gene. In a study on molecular identification of necrophagous Muscidae and Sarcophagidae fly species by COI, the data showed very low intraspecific sequence distances and species-level monophylies . The COI gene was used successfully in conducting an assessment of genetic variability in dog ticks . Studying genetic diversity in H. dromedarii tick populations is very important if some populations were different by harboring certain disease causing agents or by having acaricide resistance genes. Therefore, future studies should include as many geographical areas as possible in the UAE. This is very important in having successful disease and acaricide resistance management programs.
In conclusion, the current study demonstrated that H. dromedarii ticks in the UAE are very similar at the genetic level based on the four sampled locations. Also, the RAPD-PCR and the sequencing of the COI gene produced similar results. To our knowledge, the present study provides the first report of genetic diversity in H. dromedarii in the UAE.
We thank the Emirates Foundation in the UAE for providing the research grant from which this work was done.