2urn:lsid:arphahub.com:pub:45048D35-BB1D-5CE8-9668-537E44BD4C7Eurn:lsid:zoobank.org:pub:91BD42D4-90F1-4B45-9350-EEF175B1727AZooKeysZK1313-29891313-2970Pensoft Publishers10.3897/zookeys.478.80374599Research ArticleArachnidaMolecular biologyMalaysiaMolecular identification of blood meal sources of ticks (Acari, Ixodidae) using cytochrome b gene as a genetic markerChe LahErnieenor Faraliana12erniee@imr.gov.myYaakopSalmah2AhamadMariana1Md NorShukor2Acarology Unit, Infectious Diseases Research Centre, Institute for Medical Research, Jalan Pahang, 50588, Kuala Lumpur, Malaysia
School of Environmental and Natural Resource Sciences, Faculty of Sciences and Technology, University Kebangsaan Malaysia, 43600 Bangi, Selangor, Malaysia
Corresponding author: Ernieenor Faraliana Che Lah (erniee@imr.gov.my)
Academic editor: D. Apanaskevich
2015280120154782743FFE0FFD1-6A0E-FFE5-B230-FF83FFA01451A6589F9C-293D-4EBF-B419-6FF2E36779225787090506201430122014Ernieenor Faraliana Che Lah, Salmah Yaakop, Mariana Ahamad, Shukor Md NorThis is an open access article distributed under the terms of the Creative Commons Attribution License (CC BY 4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.http://zoobank.org/A6589F9C-293D-4EBF-B419-6FF2E3677922
Blood meal analysis (BMA) from ticks allows for the identification of natural hosts of ticks (Acari: Ixodidae). The aim of this study is to identify the blood meal sources of field collected on-host ticks using PCR analysis. DNA of four genera of ticks was isolated and their cytochrome b (Cyt b) gene was amplified to identify host blood meals. A phylogenetic tree was constructed based on data of Cyt b sequences using Neighbor Joining (NJ) and Maximum Parsimony (MP) analysis using MEGA 5.05 for the clustering of hosts of tick species. Twenty out of 27 samples showed maximum similarity (99%) with GenBank sequences through a Basic Local Alignment Search Tool (BLAST) while 7 samples only showed a similarity range of between 91–98%. The phylogenetic trees showed that the blood meal samples were derived from small rodents (Leopoldamyssabanus, Rattustiomanicus and Sundamysmuelleri), shrews (Tupaiaglis) and mammals (Tapirusindicus and Prionailurusbengalensis), supported by 82–88% bootstrap values. In this study, Cyt b gene as a molecular target produced reliable results and was very significant for the effective identification of ticks’ blood meal. The assay can be used as a tool for identifying unknown blood meals of field collected on-host ticks.
TicksBlood mealVector controlhostsCytochrome bNational Institute of Health, Ministry of Health Malaysia and University Kebangsaan MalaysiaCitation
Che Lah EF, Yaakop S, Ahamad M, Md Nor S (2015) Molecular identification of blood meal sources of ticks (Acari, Ixodidae) using cytochrome b gene as a genetic marker. ZooKeys 478: 27–43. doi: 10.3897/zookeys.478.8037
Introduction
The tick is a member of the class Arachnida that belongs to the sub-class Acari. It relies heavily on other animals as hosts to complete their life cycle (Mihalca et al. 2012). The tick is the second important vector after mosquitoes and is able to transmit disease in both man and livestock (Dalgic et al. 2010, Smith et al. 2010). Ticks may transmit a range of disease agents that are of medical and veterinary importance when taking a blood meal. Ticks that live in close proximity to humans, especially those on rodents, play a significant role in the transmission of several diseases (Chul-Min et al. 2006). Besides their role as hosts, rodents also serve as reservoirs of tick-borne pathogens (Paulauskas et al. 2009, Harrison et al. 2010, Paziewska et al. 2010). The close association between rodents and human is a risk factor for the transmission of diseases such as rickettsiosis, babesiosis and Lyme borreliosis (Killilea et al. 2008, Kia et al. 2009, Mihalca et al. 2012).
Research on the identification of a natural host for ticks is valuable and is the main goal of a blood meal analysis (BMA). The analysis is a tool used to identify the hosts of blood feeding arthropods (Gariepy et al. 2012). Research into the composition of blood meals of arthropod vectors has been shown to provide evidence of link vector species with their specific hosts that can be used to determine possible disease reservoirs (Scott et al. 2012, Onder et al. 2013). Knowledge of vector host identification of ticks is critical in understanding the transmission cycle of vector-borne diseases (Christine 2011) and thus can provide important information to public health stakeholders for the development of more effective control strategies. Initially, attempts to identify blood meals of arthropods were made using serological techniques such as precipitin and latex agglutination tests (Walker and Davies 1971, Nevil and Anderson 1972). An improvement on the sensitivity of blood meal origin analysis was later studied using enzyme-linked immunosorbent assays (ELISA) (Elias et al. 2010, Maria Stella et al. 2012). While those techniques continue to provide valuable and insightful data, the identification of many blood meal sources is limited to just the order or family level and species could only be identified when antibodies were available. The advert of the polymerase chain reaction (PCR) and the availability of DNA sequence data of various vertebrates quickly led to a shift towards DNA-based identification (Mukabana et al. 2002, Humair et al. 2007, Kent 2009, Garros et al. 2011).
The progressive development of molecular analysis has resulted in the use of markers that were originally designed to study phylogenetic relationships among vertebrates being applied to determine the blood meal origin of haematophagous arthropods (Boakye et al. 1999, Ninio et al. 2011). Tobolewski et al. (1992) reported that DNA techniques are much faster and more reliable to determine host species compared to other techniques that depend on antibodies. Humair et al. (2007) used 12S rDNA as a genetic marker for identification of blood meal sources in Ixodesricinus ticks while Pichon et al. (2003) used a set of universal primers to amplify part of the vertebrate 18S rRNA gene followed by reverse line blot hybridization for host identification of Ixodesricinus. Understanding the origin of blood meals will give some information on the species involved in disease transmission cycles, as well as the disease risk posed to human.
Analysis of blood meals of local ticks to determine their natural hosts is still poorly understood. There is a critical necessity to document information of potential hosts for local ticks as part of the nation’s preparedness for emerging and re-emerging infections. Thus, the aim of this study was to identify the blood meal sources of ticks in Peninsular Malaysia based on the cytochrome b (Cyt b) gene. The obtained Cyt b genes were then used to construct a phylogenetic tree for the clustering of hosts of tick species.
Materials and methodsStudy area and ticks sampling
On-host ticks were collected from four different states of Peninsular Malaysia (Figure 1). The study sites were in Janda Baik, Pahang; Labu, Negeri Sembilan; Hulu Langat, Selangor; and Gunung Tebu, Terengganu. The collections were conducted from March to July 2012. The localities were chosen based on the records from available previous data of high numbers of tick infestation on small animals. The ecology for all localities was similar; consisting of mainly pristine tropical lowland rainforest, secondary growths, scrubs and riverine vegetation. One hundred wire traps were used to capture wild rodents and tree-shrews in each study site. Traps were placed on the ground and on tree branches along existing trails at approximately five meter intervals. Traps were baited with bananas, oil palm fruits, tapioca or potatoes and checked once daily for 5 consecutive days of trapping. Caught animals were placed in cloth bags and brought back to the Institute for Medical Research (IMR). The animals were anesthetized with chloroform before screening and the ticks were collected in the laboratory. The epidemiology data such as locality and ecology were recorded.
The species of animals were identified by their morphological traits following Medway (1983), Francis (2008) and Payne et al. (2005). A total of 27 on-host ticks were collected either using sterile soft forceps or sharpened wooden applicator sticks. The ticks were then kept individually in vials containing 70% ethanol prior to identification using specific illustrated morphological taxonomic keys (Kohl 1957, Walker et al. 2003).
DNA extraction of blood meals in ticks
Prior to DNA extraction, each engorged tick was individually washed 3 times with sterile distilled water. Extraction of DNA using QIAamp Mini Kit (Qiagen, Germany) was performed according to the manufacturer’s protocols. DNA of ticks was extracted by adding 80 µl of PBS buffer and 100 µl of ATL buffer into the sample. The ticks were then macerated using sterile tips for 5 minutes before adding of 20 µl of proteinase K. The samples were incubated at 56 °C (6 hours) for complete lyses. The following steps were the same as those in the manufacture’s protocols. The DNA was then used for subsequent PCR.
PCR amplification of the mtDNA Cyt b gene
A portion of mitochondrial DNA, Cyt b gene was amplified by Polymerase Chain Reaction (PCR) with a set of vertebrate-universal primers and reaction conditions as described by Kent and Norris (2005). The primers used were UNFOR403 and UNREV1025 (Table 1) that preferentially amplified a 623-bp region of the Cyt b gene from the mitochondrial DNA of vertebrates. The PCR reactions were conducted in 50 µl reaction tubes with the following reagents: 25 µl Taq PCR master 2X, 2.5 µl of 0.5 µM of each primer, 10 µl of nuclease free water and 10 µl of DNA template. The amplification program consists of a total of 35 cycles, denaturing at 94 °C for 3 min, annealing at 52 °C for 1 min, and extension at 72 °C for 1 min, with an initial denaturation at 94 °C for 1 min. PCR was carried out using an Eppendorf Master Cycler Personal machine (Eppendorf, Germany). For each PCR reaction, a negative control containing the DNA of unfed ticks and a positive control containing vertebrate DNA was included. The amplicons were visualized in 1.5% agarose gels stained with ethidium bromide and viewed under an ultraviolet trans-illuminator (wavelength 254 nm).
List of primers used for PCR amplification
Gene
Primer name
Sequences (5’–3’)
Cyt b
UNFOR403
5'-TGA GGA CAA ATA TCA TTC TGA GG-3'
UNREV1025
5'-GGT TGT CCT CCA ATT CAT GTT A-3'
Sequencing analysis and Alignment
The PCR product was excised with a sterile gel cutter and purified using 5 Prime PCR Agarose Gel Extract Mini Kit (Hamburg, Germany) according to the manufacturer’s protocols. The purified product was then sent to the sequencing service company, Medigene Sdn. Bhd. in Petaling Jaya, Selangor. The sequencing was bi-directional for all specimens and the primer combination for this step was the same as that used in the PCR amplification. Sequencing results were exported as FASTA sequence files. The Cyt b gene sequences of samples were aligned using ClustalW multiple alignment of BioEdit (Thompson et al. 1994, Hall 2005).
BLAST analysis
The obtained sequences were then compared with available sequences in the GenBank database using the Basic Local Alignment Search Tool search (NCBI website, http://www.ncbi.nlm.nih.gov/BLAST/) for the identification of the host species. This approach was reported to be simple and robust for rapid comparison of query sequences to database sequences leading to species identification (Altschul et al. 1990, Mitler et al. 2010). The approach enabled the similarity of sequences to be measured depending on several criterias such as expected value, maximum identical, query coverage and maximum score.
Clustering analysis
The clustering analysis for all sequences of hosts was carried out by performing phylogenetic analysis using MEGA software (version 5.05). For distance analysis, a neighbor-joining (NJ) tree was generated from a Kimura two-parameter distance matrix. Maximum-parsimony (MP) analysis was performed with Tree-Bisection-Reconnection (TBR) heuristic algorithm to reconstruct a character based phylogenetic tree. Internal branches of both trees were statistically supported by bootstrapping with 1,000 replications. In this study, Apodemussylvaticus (GenBank Accession no. AJ298599.1) was selected as an outgroup for Cyt b gene.
Results
A total of 27 engorged ticks collected from four different localities (Table 2) were amplified from Cyt b using PCR analysis. Six species of hosts comprising three rodents (Leopoldamyssabanus, Sundamysmuelleri and Rattustiomanicus); a shrew, Tupaiaglis and two mammals identified as Tapirusindicus and Prionailurusbengalensis, were identified. Mixed stages of on-host ticks were collected from the animals and identified into four genera: Ixodes, Amblyomma, Haemaphysalis and Dermacentor.
Details of ticks samples used in analysis.
Code samples
Tick species
Locality
Ecology
K1a
Amblyommatestudinarium
Krau Wildlife Reserved, Pahang
Pristine tropical rainforest
K2a
Dermacentor sp.
Krau Wildlife Reserved, Pahang
Pristine tropical rainforest
K3c
Dermacentor sp.
Krau Wildlife Reserved, Pahang
Pristine tropical rainforest
K4b
Amblyomma sp.
Krau Wildlife Reserved, Pahang
Pristine tropical rainforest
K5b
Amblyomma sp.
Krau Wildlife Reserved, Pahang
Pristine tropical rainforest
SBN 01
Ixodesgranulatus
Labu, Negeri Sembilan
Scrubs
SBN12_1
Dermacentor sp.
Labu, Negeri Sembilan
Scrubs
SBN23_1
Ixodesgranulatus
Labu, Negeri Sembilan
Scrubs
SBN 14
Ixodes sp.
Labu, Negeri Sembilan
Scrubs
JBB01_1
Dermacentor sp.
Janda Baik, Pahang
Riverine vegetation
JBB03_1
Haemaphysalis sp.
Janda Baik, Pahang
Riverine vegetation
JBB03_2
Haemaphysalis sp.
Janda Baik, Pahang
Riverine vegetation
JBB03_3
Haemaphysalis sp.
Janda Baik, Pahang
Riverine vegetation
JBB03_4
Haemaphysalis sp.
Janda Baik, Pahang
Riverine vegetation
HL01
Ixodesgranulatus
Hulu Langat, Selangor
Pristine tropical rainforest
HL02
Ixodesgranulatus
Hulu Langat, Selangor
Pristine tropical rainforest
HL03_2
Ixodesgranulatus
Hulu Langat, Selangor
Pristine tropical rainforest
HL03
Ixodesgranulatus
Hulu Langat, Selangor
Pristine tropical rainforest
HL04_15
Ixodesgranulatus
Hulu Langat, Selangor
Pristine tropical rainforest
HL02_4
Ixodesgranulatus
Hulu Langat, Selangor
Pristine tropical rainforest
HL02_5
Ixodesgranulatus
Hulu Langat, Selangor
Pristine tropical rainforest
HL02_3
Ixodesgranulatus
Hulu Langat, Selangor
Pristine tropical rainforest
HL07
Ixodesgranulatus
Hulu Langat, Selangor
Pristine tropical rainforest
HL02_1
Ixodesgranulatus
Hulu Langat, Selangor
Pristine tropical rainforest
GT23_2
Ixodesgranulatus
Gunung Tebu, Terengganu
Secondary growth
GT23_3
Ixodesgranulatus
Gunung Tebu, Terengganu
Secondary growth
GT23_5
Ixodesgranulatus
Gunung Tebu, Terengganu
Secondary growth
PCR products and BLAST analysis
Vertebrate DNA was successfully amplified from 27 engorged ticks. The amplification of a single fragment encoding a 623 bp sequence of the Cyt b gene yielded the expected amplification products (Figure 2). Negative control (unfed ticks) yielded no PCR product implying that only host DNA patterns were detected in the amplifying specimens. All blood meals were identified down to species level. In general, the percent similarity between queried (unknown) sequences and the closest match in GenBank was between 91–100% (Table 3). Twenty out of 27 sequences (74.1%) showed maximum homology (99%) while 2 sequences (HL04_15 and SBN23_1) displayed 98% similarity with the online database. One sequence (HL03) showed maximum homology to L.sabanus with a 97% similarity while four sequences from Janda Baik matched to S.muelleri with a 91%. Some of the differences were probably due to intraspecific variation and in some cases, it could be due to poor sequence quality, particularly when unassigned nucleotides (n’s) were present in the sequence.
Amplification of Cyt b gene produced 623 bp of PCR products from ticks species. Lane 1: unfed ticks (negative control); Lane 2: vertebrate DNA (positive control); Lanes 3–8: DNA of field-collected ticks (Ixodes sp., Dermacentor sp., Amblyomma sp., Amblyommatestudinarium, Haemaphysalis sp., Haemaphysalis sp.,) and Lanes M1, M2: 100 bp DNA ladder (Bioron, Germany).
https://binary.pensoft.net/fig/37668
Blasting results against available sequences in GenBank.
Code samples
Tick species
Host species (morphological)
% similarity with GenBank (species)
K1a
Amblyommatestudinarium
Tapirusindicus
99 (Tapirusindicus)
K2a
Dermacentor sp.
Tapirusindicus
99 (Tapirusindicus)
K3c
Dermacentor sp.
Tapirusindicus
99 (Tapirusindicus)
K4b
Amblyomma sp.
Tapirusindicus
99 (Tapirusindicus)
K5b
Amblyomma sp.
Tapirusindicus
99 (Tapirusindicus)
SBN 01
Ixodesgranulatus
Rattustiomanicus
99 (Rattustiomanicus)
SBN12_1
Dermacentor sp
Rattustiomanicus
99 (Rattustiomanicus)
SBN23_1
Ixodesgranulatus
Rattustiomanicus
98 (Rattustiomanicus)
SBN 14
Ixodes sp.
Rattustiomanicus
99 (Prionailurusbengalensis)
JBB01_1
Dermacentor sp.
Tupaiaglis
99 (Tupaiaglis)
JBB03_1
Haemaphysalis sp.
Sundamysmuelleri
91 (Sundamysmuelleri)
JBB03_2
Haemaphysalis sp.
Sundamysmuelleri
91 (Sundamysmuelleri)
JBB03_3
Haemaphysalis sp.
Sundamysmuelleri
91 (Sundamysmuelleri)
JBB03_4
Haemaphysalis sp.
Sundamysmuelleri
91 (Sundamysmuelleri)
HL01
Ixodesgranulatus
Leopoldamyssabanus
99 (Leopoldamyssabanus)
HL02
Ixodesgranulatus
Leopoldamyssabanus
99 (Leopoldamyssabanus)
HL03_2
Ixodesgranulatus
Leopoldamyssabanus
99 (Leopoldamyssabanus)
HL03
Ixodesgranulatus
Leopoldamyssabanus
97 (Leopoldamyssabanus)
HL04_15
Ixodesgranulatus
Sundamysmuelleri
98 (Sundamysmuelleri)
HL02_4
Ixodesgranulatus
Leopoldamyssabanus
99 (Leopoldamyssabanus)
HL02_5
Ixodesgranulatus
Leopoldamyssabanus
99 (Leopoldamyssabanus)
HL02_3
Ixodesgranulatus
Leopoldamyssabanus
99 (Leopoldamyssabanus)
HL07
Ixodesgranulatus
Leopoldamyssabanus
99 (Leopoldamyssabanus)
HL02_1
Ixodesgranulatus
Leopoldamyssabanus
99 (Leopoldamyssabanus)
GT23_2
Ixodesgranulatus
Leopoldamyssabanus
99 (Leopoldamyssabanus)
GT23_3
Ixodesgranulatus
Leopoldamyssabanus
99 (Leopoldamyssabanus)
GT23_5
Ixodesgranulatus
Leopoldamyssabanus
99 (Leopoldamyssabanus)
Clustering inferens
From the 27 aligned DNA sequences, a total of 575-bp portion of the Cyt b gene was used for analysis. Out of 575 characters from Cyt b fragments, 252 variable sites were detected, among which 203 (35.3%) variable characters were parsimony-informative while 49 (8.5%) characters were parsimony-uninformative. Additionally, the conserved sites were constituted by 323 (56.1%) characters showing that Cyt b gene is a very conserved gene in the mtDNA.
The mitochondrial Cyt b gene sequences of host species were grouped into 2 major clades by NJ and MP analysis; one group of small rodents and another consisting of shrews and large mammals. The tree topology showed that all 27 host sequences examined fell into two distinct genetic lineages: Clade A (consists of L.sabanus, R.tiomanicus and S.muelleri) and Clade B (consists of T.glis, P.bengalensis and T.indicus). NJ tree topology revealed a distinction with 95% bootstrap value for Clade A but a lower bootstrap value of 85% for Clade B (Figure 3). Significant grouping of L.sabanus, R.tiomanicus and S.muelleri in each independent monophyletic subclade was obtained with a 100% bootstrap value. The NJ tree also showed the infestation of Ixodes ticks on 3 species of different hosts that is L.sabanus, R.tiomanicus and P.bengalensis while the Haemaphysalis ticks infested only on a single host, S.muelleri.
Neighbor-joining tree constructed from 28 sequences (including one outgroup sequence) of the Cyt b gene. The numbers at the branches stand for bootstrap values 70% and above of 1000 replications. Genera of ticks represented by blue for Ixodes sp., orange for Dermacentor sp., violet for Haemaphysalis sp. and green for Amblyomma sp.
https://binary.pensoft.net/fig/37669
Seven parsimonious trees were produced by the MP analysis using equally weighted TBR. The best tree had 497 steps (Figure 4) with a consistency index of 0.67114, a homoplasy index of 0.70422 and a retention index of 0.88804. The MP tree is concordant with the NJ tree as it shares a relatively similar tree topology of phylogenetic relationships. Two distinct monophyletic clades are clearly shown on the phylogenetic tree of potential hosts of ticks. These small rodents and large mammals were fully supported by high bootstrap values of 82% and 88%, respectively.
Maximum parsimony tree constructed from 28 sequences (including one outgroup sequence) of the Cyt b gene. The numbers at the branches stand for bootstrap values 70% and above of 1000 replications. Genera of ticks represented by blue for Ixodes sp., orange for Dermacentor sp., violet for Haemaphysalis sp. and green for Amblyomma sp.
https://binary.pensoft.net/fig/37670Discussion
The study produced the first analysis of host identification for local ticks that identifies animals to the species level by detecting their Cyt b gene in engorged on-host ticks. The usefulness of blood meal analysis in determining the identity of the host species by PCR has been demonstrated by a recent study (Ernieenor et al. 2012). The identification of the vertebrate hosts using molecular detection of the Cyt b gene displayed an exact match and this adds more confidence to the use of this analysis. In the present case, all the samples were conclusively identified to species level. The results of this study also indicated that local ticks are not host specific because they have a wide host range from small rodents to mammals.
In this study, the Cyt b gene was successfully proven as a discriminatory molecular marker for the identification of host DNA in ticks. Cyt b was selected as the target gene because of its track record in blood meal identification assays and its utilization in developing mammalian phylogeny (Kirstein and Gray 1996, Pierce et al. 2009). Moreover, the Cyt b gene has been widely used to make reliable blood meal identification of arthropods inclusive ticks due to its high copy numbers as mitochondrial genes and sufficient genetic variations at the primary sequence level among vertebrate taxa (Ngo and Kramer 2003, Maleki-Ravasan et al. 2009).
The most significant aspect of the method is its sensitivity to detect minuscule amounts of host DNA. In this study, host DNA could be detected in the 27 engorged ticks as determined by electrophoresis on agarose gel. The successful identification presumably because ticks sample was in the freshly engorged state, sufficient intact host DNA was present in the midgut. Adult stage of the ticks was used in this analysis may also eventually give positive identification of host. Cadenas et al. (2007) had similar complications, reporting that the ability to detect host DNA in adult ticks was significantly higher than in nymphs, probably due to the higher quantity of blood ingested by ticks during their previous bloodmeals as nymphs and larvae.
This study shows that Ixodesgranulatus may probably pose greater problems than most other ticks because of its capability to infest various hosts such as rodents and larger animals. This finding is in accordance with studies that reported I.granulatus as the most notable tick infesting rodents (Mariana et al. 2008, Nadchatram 2008, Paulauskas et al. 2011, Mihalca et al. 2012, Mihalca and Sandor 2013). The species is a common acarine ectoparasite of rodents in Malaysia with wide distribution extends from Southeast Asia to eastern India and China (Nadchatram 2008). Ixodes are the most common genus of ticks that feed on humans (Briciu et al. 2011) may also infest rodents that can serve as reservoir hosts for human diseases (Marchette 1966, Schorn et al. 2011). Several genera of ticks found feeding on one individual of host is a common feeding behavior in wild rodents (Mariana et al. 2005). The sharing of host by different genera or species of ticks is important mainly for the bridging of microbial pathogens through reservoir hosts (Bursali et al. 2010). The ecological importance of reservoir hosts is greater if the animals are also common hosts for competent ticks (Mihalca et al. 2012).
An interesting finding was observed in this study where Ixodes ticks (SBN 14) collected from R.tiomanicus gave a high similarity (99%) of blood meal to P.bengalensis. It was probably due to incomplete or an interruption of feeding that occurred before the tick attached itself to another host (Allan et al. 2010). Another possible reason is the short duration of feeding on the first host before the partially fed ticks dropped off-host and attached to a new host (P.bengalensis) for a blood meal. A similar finding known as co-feeding transmission was also reported by Gern and Rais (1996) and Shih and Spielman (1993).
The finding of T.indicus as one of the hosts for Amblyomma ticks (K1a, K4b and K5b) was rare because the ticks were generally host specific on wild reptiles and amphibians (Pandit et al. 2011). That appears to indicate that the population of reptiles and amphibians surrounding the locality are scarcely found and the ticks have no choice to feed on any other animals available in the area. Rodents such as Leopoldamyssabanus, Rattusdiardii and Niviventerrapit were also reported as alternative hosts for Amblyomma ticks (Nadchatram et al. 1966, Paramasvaran et al. 2009).
Partial Cyt b mtDNA gene sequences used in this study seems to be effective in identifying the phylogenetic relationships between ticks and their hosts. This is because both NJ and MP tree topology showed very clear distinction between groups of small rodents, shrews and large mammals with a highly supported monophyletic clade. In small rodents, the node is solved in the NJ and MP tree as two monophyletic clades, representing the species of L.sabanus, R.tiomanicus and S.muelleri. The tree topologies also show that T.glis formed their own distinct monophyletic clade separated from other larger mammals which consists of P.bengalensis and T.indicus. The information obtained has further corroborated the morphological identification and classification of T.glis which comes under the group of shrews. Such knowledge gained through host preference studies is essential to understanding the relationship of host and vector and their roles in the enzootic transmission cycle (Ngo and Kramer 2003). Therefore, analyses of ticks bloodmeals will create a better understanding of epidemiologically important vectors-their hosts that can be lead to the design of more effective control strategies (Kirstein and Gray 1996).
Of particular interest was the high numbers of ticks infesting small rodents compared to larger mammals. This finding is in agreeable with a previous study which reported that the abundance of small rodents compensate even if the intensity of tick parasitism on them is smaller than the larger hosts (Randolph 2004). The capability of rodents as host for ticks has been extensively investigated until to date mainly because they can be easily captured in large numbers and are easier to handle or maintain in the laboratory (Gern and Humair 2002). This is in line with our catch which was represented mainly by rodents as they were abundantly found in all localities.
For further studies, it is recommended that ticks are collected from sites frequently visited by animals including the wallows, wildlife main trails and river banks to include a more diverse group of animals in order to generate better results (Pichon et al. 2005, Pierce et al 2009). A bigger sample size with heterogeneous host species collection needs to be examined in parallel with different molecular markers to enable researchers to draw reliable conclusions and provide alternative views on the relationship between ticks and their natural hosts in Malaysia.
Conclusion
The PCR direct sequencing system using vertebrate Cyt b gene is a potential screening tool for the identification of ticks’ blood meal. The Cyt b gene was selected for this study based on their higher nucleotide variations for the effective identification of hosts. Blood meal identification of field collected ticks by molecular methods offer a direct and efficient approach for understanding the contributions of both competent and incompetent hosts for the transmission dynamics of tick-borne diseases. Furthermore, this valuable information can confirm a strong association between hard ticks and hosts (especially rodents) and this will assist public health officials with efforts to outline an effective tick-borne diseases control program.
Acknowledgements
The authors wish to thank the Director-General of Health, Malaysia, for permission to publish this paper. We wish to thank staff of the Acarology Unit, IMR for their assistance in the field. For the revision of English of the manuscript, we thank Serina Abdul Rahman. The study was supported by National Institute of Health grant (Code: JPP-IMR 11-010) from the Ministry of Health, Malaysia and ERGS/1/2011/STWN/UKM/03/9.
ReferencesAllanBFGoesslingLSStorchGAThachRE (2010) Blood meal analysis to identify reservoir hosts for Amblyommaamericanum ticks.16(3): 433–440. doi: 10.3201/eid1603.090911AltschulSFGishWMillerWMyersEWLipmanDJ (1990) Basic local alignment search tool.215(3): 403–410. doi: 10.1016/S0022-2836(05)80360-2BoakyeDATangJTrucPMerriweatherAUnnaschTR (1999) Identification of blood meals in haematophagous Diptera by cytochrome B heteroduplex analysis.13: 282–287. doi: 10.1046/j.1365-2915.1999.00193.xBriciuVTTitilincuATatulescuDFCarstinaDLefkaditisMMihalcaAD (2011) First survey on hard ticks (Ixodidae) collected from humans in Romania: possible risks for tick-borne diseases.54(2): 199–204. doi: 10.1007/s10493-010-9418-0BursaliATekinSOrhanMKeskinAOzkanM (2010) Ixodid ticks (Acari: Ixodidae) infesting human in Tokat Province of Turkey: species diversity and seasonal activity.35: 180–186. doi: 10.1111/j.1948-7134.2010.00075.xCadenasFMRaisOHumairPFDouetVMoretJGernL (2007) Identification of host bloodmeal source and Borreliaburgdorferi sensu-lato in field-collected Ixodesricinus ticks in Chaumont (Switzerland).44(6): 1109–1117. doi: 10.1603/0022-2585(2007)44[1109:IOHBSA]2.0.CO;2ChristineZ (2011) Development of genetic barcodes for hosts of the blacklegged tick (Ixodesscapularis) in Southern New York. Dissertation, New York, United State: Fordham University.Chul-MinKYing-HuaYDo-HyeonY et al. (2006) Tick-borne rickettsial pathogens in ticks and small mammals in Korea.72: 5766–76. doi: 10.1128/AEM.00431-06DalgicAKandoganTKavakHAriAErkanNOzuerMZ (2010) Ticks in the external auditory canal.17(2): 190–192.EliasSLMarianaSFLuiz ClaudioMDOJeronimoACarlosBM (2010) Blood meal identification of selected mosquitoes in Rio De Janeiro, Brazil.26(1): 18–23. doi: 10.2987/09-5914.1ErnieenorFCLMarianaAMohd SubailHHoTM (2012) Establishment of a molecular tool for blood meal identification in Malaysia.2(3): 223–227. doi: 10.1016/S2221-1691(12)60046-XFrancisCM (2008) Princeton Press, London, 392 pp.GariepyTDLindsayROgdenNGregoryTR (2012) Identifying the last supper utility of the DNA barcode library for blood meal identification in ticks.12: 646–652. doi: 10.1111/j.1755-0998.2012.03140.xGarrosCGardesLAlleneXRakotoarivonyIViennetERossiSBalenghienT (2011) Development of an allele-specific multiplex PCR assay for the identification of blood meal source in Culicoides (Ceratopogonidae: Diptera). Applications on Palaeartic bitting midges species.11: 1103–1110. doi: 10.1016/j.meegid.2011.04.002GernLHumairPF (2002) Ecology of Borreliaburgdorferi sensu lato in Europe. In: GrayJSKahlOLaneRSStanekG (Eds) . CABI, 149–174. doi: 10.1079/9780851996325.0149GernLRaisO (1996) Efficient transmission of Borreliaburgdorferi between cofeeding Ixodesricinus ticks (Acari: Ixodidae).33: 189–192.HallT (2005) Bioedit version 7.0.4. Department of Microbiology, North Carolina State University.HarrisonAScantleburyMMontgomeryWI (2010) Body mass and sex-biased parasitism in wood mice Apodemussylvaticus.119(7): 1099–1104. doi: 10.1111/j.1600-0706.2009.18072.xHumairPFDouetVCadenasFMSchoulsLMVan De PolIGemL (2007) Molecular identification of blood meal source in Ixodesricinus ticks using 12S rDNA as a genetic marker.44: 869–880. doi: 10.1603/0022-2585(2007)44[869:MIOBSI]2.0.CO;2KentRJ (2009) Molecular methods for arthropod bloodmeal identification and applications to ecological and vector-borne disease studies.9: 4–18. doi: 10.1111/j.1755-0998.2008.02469.xKentRJNorrisD (2005) Identification of mammalian blood meals in mosquito by a multiplexed Polymerase Chain Reaction targeting cytochrome b.73(2): 336–342.KiaEBMoghddas-SaniHHassanpoorH et al. (2009) Ectoparasites of rodents captured in Bandar Abbas, Southern Iran.3: 44–9.KillileaMESweiALaneRSBriggsCJOstfeldRS (2008) Spatial dynamics of Lyme disease: A review.5(2): 167–195. doi: 10.1007/s10393-008-0171-3KirsteinFGrayJS (1996) A molecular marker for the identification of the zoonotic reservoirs of Lyme borreliosis by analysis of the blood meal in its European vector Ixodesricinus.62: 4060–4065.KohlsGM (1957) Ticks (Ixodidae) of Borneo and Malaya.28: 65–94.Maleki RavasanNOshaghiMAJavadianERassiYSadraeiJMohtaramiF (2009) Blood meal identification in field-captured sand flies: Comparison of PCR-RFLP and ELISA assays.3(1): 8–18.MarchetteNJ (1966) Rickettsiosis (Tick typhus, Q-fever, Urban-typhus) in Malaya.4: 393–71.Maria StellaBBAndrea de BarrosPVNatalia GoesB (2012) Standardization of Enzyme Linked Immunosorbent Assay (ELISA) for the identification of blood meal in black flies (Diptera: Simuliidae).14(1): 12–16.MarianaAZuraidawatiZHoTMMohd KulaimiBSalehIShukorMNShahrul-AnuarMS (2005) A survey of ectoparasites in Gunung Stong Forest Reserve, Kelantan, Malaysia.36(5): 1125–1131.MarianaAZuraidawatiZHoTm et al. (2008) Ticks (Ixodidae) and other ectoparasites in Ulu Muda Forest Reserve, Kedah, Malaysia.39: 496–506.MedwayL (1983) . Oxford University Press, Kuala Lumpur.MihalcaADSandorAD (2013) The role of rodents in the ecology of Ixodesricinus and associated pathogens in Central and Eastern Europe.3(56): 1–3.MihalcaADDumitracheMOSandorADMagdasCOlteanM et al. (2012) Tick parasites of rodents in Romania: host preferences, community structure and geographical distribution.5(266): 1–7.MittlerTLevyMChadFKarenS (2010) MULTBLAST: A web application for multiple BLAST searches.5(5): 224–226. doi: 10.6026/97320630005224MukabanaWRTakkenWSedaPKillenGFHawleyWAKnolsBGJ (2002) Analysis of arthropod blood meals using molecular genetic markers.18(11): 505–509. doi: 10.1016/S1471-4922(02)02364-4NadchatramM (2008) The beneficial rain forest ecosystem with environmental effects on zoonoses involving ticks and mites (acari), a Malaysian perspective and review.25: 1–92.NadchatramMDamrowRNgCK (1966) Parasitic acarina of the mammals.31: 129–140.NevilEMAndersonD (1972) Host preference of Culicoides midges (Diptera: Ceratopogonidae) in South Africa as determined by precipitin test and light trap catches.39: 147–152.NgoKAKramerLD (2003) Identification of mosquito blood meals using polymerase chain reaction (PCR) with order-specific primers.40: 215–222. doi: 10.1603/0022-2585-40.2.215NinioCDAugotJCDelecolleBDufourBDepaquitJ (2011) Contribution to the knowledge of Culicoides (Diptera: Ceratopogonidae) host preference in France.108: 657–663. doi: 10.1007/s00436-010-2110-9OnderOShaoWKempsBDLamHBrissonD (2013) Identifying sources of blood meals using unidentified tandem mass spectral libraries.4: 1746. doi: 10.1038/ncomms2730PanditPBandivdekarRGeevargheseGPandeSMandkaO (2011) Tick infestation on wild snakes in Northern part of Western Ghats of India.48(3): 504–507. doi: 10.1603/ME10164ParamasvaranSSaniRAHassanLKrishnasamyMJefferyJOothumanP et al. (2009) Ectoparasite fauna of rodents and shrews from four habitats in Kuala Lumpur and the states of Selangor and Negeri Sembilan, Malaysia and its public health significance.26(3): 303–311.PaulauskasARadzijevskajaJRosefOTurcinavicieneJAmbrasieneD (2009) Infestation of mice and voles with Ixodesricinus ticks in Lithuania and Norway.58(2): 112–125. doi: 10.3176/eco.2009.2.05PaulauskasARadzijevskajaJRosefO (2011) Rodents as carriers of tick-borne zoonotic diseases and their eological impact.432: 182–183.PayneJFrancisCMPhillipsK (2005) Sabah Society, Kota Kinabalu, 332 pp.PaziewskaAZwolinskaLHarrisPDBajerASinskiE (2010) Utilisation of rodent species by larvae and nymphs of hard ticks (Ixodidae) in two habitats in NE Poland.50(1): 79–91. doi: 10.1007/s10493-009-9269-8PichonBEganDRogersMGrayJ (2003) Detection and identification of pathogens and host DNA in unfed host-seeking Ixodesricinus L. (Acari: Ixodidae).40: 723–731. doi: 10.1603/0022-2585-40.5.723PichonBRogerMEganDGrayJ (2005) Blood meal analysis for the identification of reservoir hosts of tick-borne pathogen in Ireland.5(2): 172–180. doi: 10.1089/vbz.2005.5.172PierceKAPaddockCDSumnerJWNicholsonWL (2009) Pathogen prevalence and blood meal identification in Amblyomma ticks as a mean of reservoir host determination for ehrlichial pathogens.15 (Suppl. 2): 37–38. doi: 10.1111/j.1469-0691.2008.02166.xRandolphSE (2004) Tick ecology: processes and patterns behind the epidemiological risk posed by ixodid ticks as vectors.129: S37–S65. doi: 10.1017/S0031182004004925SchornSPfisterKReulenHMahlingMSilaghiC (2011) Occurrence of Babesia spp., Rickettsia spp., and Bartonella spp. In: Ixodes ricinus in Bavarian public parks, Germany.4: 135. doi: 10.1186/1756-3305-4-135ScottMCHarmonJrTsaoJIJonesCJHicklingG (2012) Reverse line blot probe design and polymerase chain reaction optimization for blood meal analysis of ticks from the eastern United States.49(3): 697–709. doi: 10.1603/ME11162ShihCMSpielmanA (1993) Accelerated transmission of Lyme disease spirochetes by partially fed vector ticks.31: 2878–2881.ThompsonJDHigginsDGGibsonTJ (1994) CLUSTAL W: Improving the sensitivity of progressive multiple sequence alignment through sequence weighting, position-specific gap penalties and weight matrix choice.22(22): 4673–4680. doi: 10.1093/nar/22.22.4673TobolewskiJKaliszewskiMJColwellRKOliverJH (1992) Detection and identification of mammalian DNA from the gut of museum specimens of ticks.29: 1049–1051.WalkerARDaviesFG (1971) A preliminary survey of the epidemiology of bluetongue in Kenya.69(1): 47–60. doi: 10.1017/S0022172400021239WalkerARBouattourACamicasJLEstrada-PenaAHorakIGLatifAPegramRGPrestonPM (2003) Biosciences Reports, United Kingdom, 157 pp.