Molecular identification of blood meal sources of ticks (Acari, Ixodidae) using cytochrome b gene as a genetic marker

Ernieenor Faraliana Che Lah; Salmah Yaakop; Mariana Ahamad; Shukor Md Nor

doi:10.3897/zookeys.478.8037

Abstract

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 (Leopoldamys sabanus, Rattus tiomanicus and Sundamys muelleri), shrews (Tupaia glis) and mammals (Tapirus indicus and Prionailurus bengalensis), 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.

Keywords

Ticks, Blood meal, Vector control, hosts, Cytochrome b

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 Ixodes ricinus 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 Ixodes ricinus. 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 methods

Study 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.

Figure 1.

Map of the study sites in Peninsular Malaysia: 1 Hulu Langat, Selangor 2 Labu, Negeri Sembilan 3 Janda Baik, Pahang 4 Gunung Tebu, Terengganu.

Tick Morphological identification

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).

Table 1.

Download as

CSV

XLSX

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, Apodemus sylvaticus (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 (Leopoldamys sabanus, Sundamys muelleri and Rattus tiomanicus); a shrew, Tupaia glis and two mammals identified as Tapirus indicus and Prionailurus bengalensis, were identified. Mixed stages of on-host ticks were collected from the animals and identified into four genera: Ixodes, Amblyomma, Haemaphysalis and Dermacentor.

Table 2.

Download as

CSV

XLSX

Details of ticks samples used in analysis.

Code samples	Tick species	Locality	Ecology
K1a	Amblyomma testudinarium	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	Ixodes granulatus	Labu, Negeri Sembilan	Scrubs
SBN12_1	Dermacentor sp.	Labu, Negeri Sembilan	Scrubs
SBN23_1	Ixodes granulatus	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	Ixodes granulatus	Hulu Langat, Selangor	Pristine tropical rainforest
HL02	Ixodes granulatus	Hulu Langat, Selangor	Pristine tropical rainforest
HL03_2	Ixodes granulatus	Hulu Langat, Selangor	Pristine tropical rainforest
HL03	Ixodes granulatus	Hulu Langat, Selangor	Pristine tropical rainforest
HL04_15	Ixodes granulatus	Hulu Langat, Selangor	Pristine tropical rainforest
HL02_4	Ixodes granulatus	Hulu Langat, Selangor	Pristine tropical rainforest
HL02_5	Ixodes granulatus	Hulu Langat, Selangor	Pristine tropical rainforest
HL02_3	Ixodes granulatus	Hulu Langat, Selangor	Pristine tropical rainforest
HL07	Ixodes granulatus	Hulu Langat, Selangor	Pristine tropical rainforest
HL02_1	Ixodes granulatus	Hulu Langat, Selangor	Pristine tropical rainforest
GT23_2	Ixodes granulatus	Gunung Tebu, Terengganu	Secondary growth
GT23_3	Ixodes granulatus	Gunung Tebu, Terengganu	Secondary growth
GT23_5	Ixodes granulatus	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.

Figure 2.

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., Amblyomma testudinarium, Haemaphysalis sp., Haemaphysalis sp.,) and Lanes M1, M2: 100 bp DNA ladder (Bioron, Germany).

Table 3.

Download as

CSV

XLSX

Blasting results against available sequences in GenBank.

Code samples	Tick species	Host species (morphological)	% similarity with GenBank (species)
K1a	Amblyomma testudinarium	Tapirus indicus	99 (Tapirus indicus)
K2a	Dermacentor sp.	Tapirus indicus	99 (Tapirus indicus)
K3c	Dermacentor sp.	Tapirus indicus	99 (Tapirus indicus)
K4b	Amblyomma sp.	Tapirus indicus	99 (Tapirus indicus)
K5b	Amblyomma sp.	Tapirus indicus	99 (Tapirus indicus)
SBN 01	Ixodes granulatus	Rattus tiomanicus	99 (Rattus tiomanicus)
SBN12_1	Dermacentor sp	Rattus tiomanicus	99 (Rattus tiomanicus)
SBN23_1	Ixodes granulatus	Rattus tiomanicus	98 (Rattus tiomanicus)
SBN 14	Ixodes sp.	Rattus tiomanicus	99 (Prionailurus bengalensis)
JBB01_1	Dermacentor sp.	Tupaia glis	99 (Tupaia glis)
JBB03_1	Haemaphysalis sp.	Sundamys muelleri	91 (Sundamys muelleri)
JBB03_2	Haemaphysalis sp.	Sundamys muelleri	91 (Sundamys muelleri)
JBB03_3	Haemaphysalis sp.	Sundamys muelleri	91 (Sundamys muelleri)
JBB03_4	Haemaphysalis sp.	Sundamys muelleri	91 (Sundamys muelleri)
HL01	Ixodes granulatus	Leopoldamys sabanus	99 (Leopoldamys sabanus)
HL02	Ixodes granulatus	Leopoldamys sabanus	99 (Leopoldamys sabanus)
HL03_2	Ixodes granulatus	Leopoldamys sabanus	99 (Leopoldamys sabanus)
HL03	Ixodes granulatus	Leopoldamys sabanus	97 (Leopoldamys sabanus)
HL04_15	Ixodes granulatus	Sundamys muelleri	98 (Sundamys muelleri)
HL02_4	Ixodes granulatus	Leopoldamys sabanus	99 (Leopoldamys sabanus)
HL02_5	Ixodes granulatus	Leopoldamys sabanus	99 (Leopoldamys sabanus)
HL02_3	Ixodes granulatus	Leopoldamys sabanus	99 (Leopoldamys sabanus)
HL07	Ixodes granulatus	Leopoldamys sabanus	99 (Leopoldamys sabanus)
HL02_1	Ixodes granulatus	Leopoldamys sabanus	99 (Leopoldamys sabanus)
GT23_2	Ixodes granulatus	Leopoldamys sabanus	99 (Leopoldamys sabanus)
GT23_3	Ixodes granulatus	Leopoldamys sabanus	99 (Leopoldamys sabanus)
GT23_5	Ixodes granulatus	Leopoldamys sabanus	99 (Leopoldamys sabanus)

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.

Figure 3.

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.

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.

Figure 4.

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.

Discussion

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 Ixodes granulatus 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 Leopoldamys sabanus, Rattus diardii and Niviventer rapit 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.

References

Allan BF, Goessling LS, Storch GA, Thach RE (2010) Blood meal analysis to identify reservoir hosts for Amblyomma americanum ticks. Emerging Infectious Diseases 16(3): 433–440. doi: 10.3201/eid1603.090911

Altschul SF, Gish W, Miller W, Myers EW, Lipman DJ (1990) Basic local alignment search tool. Journal of Molecular Biology 215(3): 403–410. doi: 10.1016/S0022-2836(05)80360-2

Boakye DA, Tang J, Truc P, Merriweather A, Unnasch TR (1999) Identification of blood meals in haematophagous Diptera by cytochrome B heteroduplex analysis. Medical and Veterinary Entomology 13: 282–287. doi: 10.1046/j.1365-2915.1999.00193.x

Briciu VT, Titilincu A, Tatulescu DF, Carstina D, Lefkaditis M, Mihalca AD (2011) First survey on hard ticks (Ixodidae) collected from humans in Romania: possible risks for tick-borne diseases. Experimental and Applied Acarology 54(2): 199–204. doi: 10.1007/s10493-010-9418-0

Bursali A, Tekin S, Orhan M, Keskin A, Ozkan M (2010) Ixodid ticks (Acari: Ixodidae) infesting human in Tokat Province of Turkey: species diversity and seasonal activity. Journal Vector Ecology 35: 180–186. doi: 10.1111/j.1948-7134.2010.00075.x

Cadenas FM, Rais O, Humair PF, Douet V, Moret J, Gern L (2007) Identification of host bloodmeal source and Borrelia burgdorferi sensu-lato in field-collected Ixodes ricinus ticks in Chaumont (Switzerland). Journal of Medical Entomology 44(6): 1109–1117. doi: 10.1603/0022-2585(2007)44[1109:IOHBSA]2.0.CO;2

Christine Z (2011) Development of genetic barcodes for hosts of the blacklegged tick (Ixodes scapularis) in Southern New York. Dissertation, New York, United State: Fordham University.

Chul-Min K, Ying-Hua Y, Do-Hyeon Y et al. (2006) Tick-borne rickettsial pathogens in ticks and small mammals in Korea. Applied & Environmental Microbiology 72: 5766–76. doi: 10.1128/AEM.00431-06

Dalgic A, Kandogan T, Kavak H, Ari A, Erkan N, Ozuer MZ (2010) Ticks in the external auditory canal. Hong Kong Journal of Emergency Medicine 17(2): 190–192.

Elias SL, Mariana SF, Luiz Claudio MDO, Jeronimo A, Carlos BM (2010) Blood meal identification of selected mosquitoes in Rio De Janeiro, Brazil. Journal of the American Mosquito Control Association 26(1): 18–23. doi: 10.2987/09-5914.1

Ernieenor FCL, Mariana A, Mohd Subail H, Ho TM (2012) Establishment of a molecular tool for blood meal identification in Malaysia. Asian Pacific Journal of Tropical Biomedicine 2(3): 223–227. doi: 10.1016/S2221-1691(12)60046-X

Francis CM (2008) Field guide to the mammals of South-east Asia. Princeton Press, London, 392 pp.

Gariepy TD, Lindsay R, Ogden N, Gregory TR (2012) Identifying the last supper utility of the DNA barcode library for blood meal identification in ticks. Molecular Resources Ecology 12: 646–652. doi: 10.1111/j.1755-0998.2012.03140.x

Garros C, Gardes L, Allene X, Rakotoarivony I, Viennet E, Rossi S, Balenghien T (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. Infection Genetic & Evolution 11: 1103–1110. doi: 10.1016/j.meegid.2011.04.002

Gern L, Humair PF (2002) Ecology of Borrelia burgdorferi sensu lato in Europe. In: Gray JS, Kahl O, Lane RS, Stanek G (Eds) Lyme Borreliosis: Biology, Epidemiology and Control. CABI, 149–174. doi: 10.1079/9780851996325.0149

Gern L, Rais O (1996) Efficient transmission of Borrelia burgdorferi between cofeeding Ixodes ricinus ticks (Acari: Ixodidae). Journal of Medical Entomology 33: 189–192.

Hall T (2005) Bioedit version 7.0.4. Department of Microbiology, North Carolina State University.

Harrison A, Scantlebury M, Montgomery WI (2010) Body mass and sex-biased parasitism in wood mice Apodemus sylvaticus. Oikos 119(7): 1099–1104. doi: 10.1111/j.1600-0706.2009.18072.x

Humair PF, Douet V, Cadenas FM, Schouls LM, Van De Pol I, Gem L (2007) Molecular identification of blood meal source in Ixodes ricinus ticks using 12S rDNA as a genetic marker. Journal of Medical Entomology 44: 869–880. doi: 10.1603/0022-2585(2007)44[869:MIOBSI]2.0.CO;2

Kent RJ (2009) Molecular methods for arthropod bloodmeal identification and applications to ecological and vector-borne disease studies. Molecular Ecology Resources 9: 4–18. doi: 10.1111/j.1755-0998.2008.02469.x

Kent RJ, Norris D (2005) Identification of mammalian blood meals in mosquito by a multiplexed Polymerase Chain Reaction targeting cytochrome b. American Journal of Tropical Medicine Hygience 73(2): 336–342.

Kia EB, Moghddas-Sani H, Hassanpoor H et al. (2009) Ectoparasites of rodents captured in Bandar Abbas, Southern Iran. Iranian Journal of Arthropod-Borne Diseases 3: 44–9.

Killilea ME, Swei A, Lane RS, Briggs CJ, Ostfeld RS (2008) Spatial dynamics of Lyme disease: A review. Ecohealth 5(2): 167–195. doi: 10.1007/s10393-008-0171-3

Kirstein F, Gray JS (1996) A molecular marker for the identification of the zoonotic reservoirs of Lyme borreliosis by analysis of the blood meal in its European vector Ixodes ricinus. Appllied of Environmental Microbiology 62: 4060–4065.

Kohls GM (1957) Ticks (Ixodidae) of Borneo and Malaya. Study, Institute for Medical Research Malaya 28: 65–94.

Maleki Ravasan N, Oshaghi MA, Javadian E, Rassi Y, Sadraei J, Mohtarami F (2009) Blood meal identification in field-captured sand flies: Comparison of PCR-RFLP and ELISA assays. Iranian Journal of Arthropod-Borne Diseases 3(1): 8–18.

Marchette NJ (1966) Rickettsiosis (Tick typhus, Q-fever, Urban-typhus) in Malaya. Journal of Medical Entomology 4: 393–71.

Maria Stella BB, Andrea de Barros PV, Natalia Goes B (2012) Standardization of Enzyme Linked Immunosorbent Assay (ELISA) for the identification of blood meal in black flies (Diptera: Simuliidae). Revista Panamericana Infectologia 14(1): 12–16.

Mariana A, Zuraidawati Z, Ho TM, Mohd Kulaimi B, Saleh I, Shukor MN, Shahrul-Anuar MS (2005) A survey of ectoparasites in Gunung Stong Forest Reserve, Kelantan, Malaysia. Southeast Asian Journal of Tropical Medicine of Public Health 36(5): 1125–1131.

Mariana A, Zuraidawati Z, Ho Tm et al. (2008) Ticks (Ixodidae) and other ectoparasites in Ulu Muda Forest Reserve, Kedah, Malaysia. Southeast Asian Journal of Tropical Medicine of Public Health 39: 496–506.

Medway L (1983) The wild mammals of Malaya (Peninsular Malaysia and Singapore). Oxford University Press, Kuala Lumpur.

Mihalca AD, Sandor AD (2013) The role of rodents in the ecology of Ixodes ricinus and associated pathogens in Central and Eastern Europe. Frontiers in Cellular and Infections Microbiology 3(56): 1–3.

Mihalca AD, Dumitrache MO, Sandor AD, Magdas C, Oltean M et al. (2012) Tick parasites of rodents in Romania: host preferences, community structure and geographical distribution. Parasites & Vector 5(266): 1–7.

Mittler T, Levy M, Chad F, Karen S (2010) MULTBLAST: A web application for multiple BLAST searches. Bioinformation 5(5): 224–226. doi: 10.6026/97320630005224

Mukabana WR, Takken W, Seda P, Killen GF, Hawley WA, Knols BGJ (2002) Analysis of arthropod blood meals using molecular genetic markers. Trends in Parasitology 18(11): 505–509. doi: 10.1016/S1471-4922(02)02364-4

Nadchatram M (2008) The beneficial rain forest ecosystem with environmental effects on zoonoses involving ticks and mites (acari), a Malaysian perspective and review. Tropical Biomedicine 25: 1–92.

Nadchatram M, Damrow R, Ng CK (1966) Parasitic acarina of the mammals. Bulletin of National Musuem of Singapore 31: 129–140.

Nevil EM, Anderson D (1972) Host preference of Culicoides midges (Diptera: Ceratopogonidae) in South Africa as determined by precipitin test and light trap catches. Onderstepoort Journal Veterinary Research 39: 147–152.

Ngo KA, Kramer LD (2003) Identification of mosquito blood meals using polymerase chain reaction (PCR) with order-specific primers. Journal of Medical Entomology 40: 215–222. doi: 10.1603/0022-2585-40.2.215

Ninio CD, Augot JC, Delecolle B, Dufour B, Depaquit J (2011) Contribution to the knowledge of Culicoides (Diptera: Ceratopogonidae) host preference in France. Parasitology Research 108: 657–663. doi: 10.1007/s00436-010-2110-9

Onder O, Shao W, Kemps BD, Lam H, Brisson D (2013) Identifying sources of blood meals using unidentified tandem mass spectral libraries. Nature Communication 4: 1746. doi: 10.1038/ncomms2730

Pandit P, Bandivdekar R, Geevarghese G, Pande S, Mandka O (2011) Tick infestation on wild snakes in Northern part of Western Ghats of India. Journal of Medical Entomology 48(3): 504–507. doi: 10.1603/ME10164

Paramasvaran S, Sani RA, Hassan L, Krishnasamy M, Jeffery J, Oothuman P 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. Tropical Biomedicine 26(3): 303–311.

Paulauskas A, Radzijevskaja J, Rosef O, Turcinaviciene J, Ambrasiene D (2009) Infestation of mice and voles with Ixodes ricinus ticks in Lithuania and Norway. Estonian Journal of Ecology 58(2): 112–125. doi: 10.3176/eco.2009.2.05

Paulauskas A, Radzijevskaja J, Rosef O (2011) Rodents as carriers of tick-borne zoonotic diseases and their eological impact. 8th European Vertebrate Pest Management Conference Julius-Kuhn-Archiv 432: 182–183.

Payne J, Francis CM, Phillips K (2005) A field guide to the mammals of Borneo. Sabah Society, Kota Kinabalu, 332 pp.

Paziewska A, Zwolinska L, Harris PD, Bajer A, Sinski E (2010) Utilisation of rodent species by larvae and nymphs of hard ticks (Ixodidae) in two habitats in NE Poland. Experimental and Applied Acarology 50(1): 79–91. doi: 10.1007/s10493-009-9269-8

Pichon B, Egan D, Rogers M, Gray J (2003) Detection and identification of pathogens and host DNA in unfed host-seeking Ixodes ricinus L. (Acari: Ixodidae). Journal of Medical Entomology 40: 723–731. doi: 10.1603/0022-2585-40.5.723

Pichon B, Roger M, Egan D, Gray J (2005) Blood meal analysis for the identification of reservoir hosts of tick-borne pathogen in Ireland. Vector-Borne & Zoonotic Diseases 5(2): 172–180. doi: 10.1089/vbz.2005.5.172

Pierce KA, Paddock CD, Sumner JW, Nicholson WL (2009) Pathogen prevalence and blood meal identification in Amblyomma ticks as a mean of reservoir host determination for ehrlichial pathogens. Clinical Microbiology and Infectious Diseases 15 (Suppl. 2): 37–38. doi: 10.1111/j.1469-0691.2008.02166.x

Randolph SE (2004) Tick ecology: processes and patterns behind the epidemiological risk posed by ixodid ticks as vectors. Parasitology 129: S37–S65. doi: 10.1017/S0031182004004925

Schorn S, Pfister K, Reulen H, Mahling M, Silaghi C (2011) Occurrence of Babesia spp., Rickettsia spp., and Bartonella spp. In: Ixodes ricinus in Bavarian public parks, Germany.Parasitology Vectors 4: 135. doi: 10.1186/1756-3305-4-135

Scott MC, Harmon Jr, Tsao JI, Jones CJ, Hickling G (2012) Reverse line blot probe design and polymerase chain reaction optimization for blood meal analysis of ticks from the eastern United States. Journal of Medical Entomology 49(3): 697–709. doi: 10.1603/ME11162

Shih CM, Spielman A (1993) Accelerated transmission of Lyme disease spirochetes by partially fed vector ticks. Journal of Clinical Microbiology 31: 2878–2881.

Thompson JD, Higgins DG, Gibson TJ (1994) CLUSTAL W: Improving the sensitivity of progressive multiple sequence alignment through sequence weighting, position-specific gap penalties and weight matrix choice. Nucleic Acid Research 22(22): 4673–4680. doi: 10.1093/nar/22.22.4673

Tobolewski J, Kaliszewski MJ, Colwell RK, Oliver JH (1992) Detection and identification of mammalian DNA from the gut of museum specimens of ticks. Journal of Medical Entomology 29: 1049–1051.

Walker AR, Davies FG (1971) A preliminary survey of the epidemiology of bluetongue in Kenya. Journal of Hygience 69(1): 47–60. doi: 10.1017/S0022172400021239

Walker AR, Bouattour A, Camicas JL, Estrada-Pena A, Horak IG, Latif A, Pegram RG, Preston PM (2003) Ticks of domestic animals in Africa, A guide to identification of species. Biosciences Reports, United Kingdom, 157 pp.