Molecular Confirmation of Rhipicephalus Haemaphysaloides Infesting Ruminants in Wayanad, Kerala, India
Total Page:16
File Type:pdf, Size:1020Kb
Int.J.Curr.Microbiol.App.Sci (2017) 6(8): 2304-2309 International Journal of Current Microbiology and Applied Sciences ISSN: 2319-7706 Volume 6 Number 8 (2017) pp. 2304-2309 Journal homepage: http://www.ijcmas.com Original Research Article https://doi.org/10.20546/ijcmas.2017.608.271 Molecular Confirmation of Rhipicephalus haemaphysaloides Infesting Ruminants in Wayanad, Kerala, India Murikoli Nimisha1, Reghu Ravindran1*, Pradeep R. Kariyappa1, B. Mallappa Amrutha1, Prashant Somalingappa Kurbet1, Karapparambu Gopalan Ajith Kumar1, Anju Varghese1, Chundayil Kalarikkal Deepa1 and Sanis Juliet2 1Department of Veterinary Parasitology, College of Veterinary and Animal Sciences, Pookode, Lakkidi, P.O., Wayanad-673576, Kerala, India 2Department of Veterinary Pharmacology and Toxicology, College of Veterinary and Animal Sciences, Pookode, Lakkidi, P.O., Wayanad-673576, Kerala, India *Corresponding author ABSTRACT Rhipicephalus haemaphysaloides is a prevalent multi host tick species with K e yw or ds high vector potential in south India. There are several reports based on morphological identification of the tick species from various parts of the Rhipicephalus country. However, there were limited attempts for molecular confirmation haemaphysaloides , PCR, Phylogeny, of this tick species. In the present study, the presence of R. haemaphysaloides in Wayanad, Kerala was confirmed by polymerase chain Wayanad, Kerala. reaction based on amplification of mitochondrial 16S rRNA gene and V4 Article Info region of 18S rRNA gene. Two tick samples, one engorged on a bovine Accepted: male calf of Pookode, Wayanad and other engorged on a sambar deer in 21 June 2017 Available Online: Meppadi were used for the study. Phylogeny based on 16S rRNA and 18S 10 August 2017 rRNA of R. haemaphysaloides revealed genetic relatedness with Chinese isolates of the tick species. Introduction Parasitic diseases are one of the major La Fuente et al., 2008). Tick infestation obstacles for the health of the animals, results in direct damages like reduction in thereby causing severe economic constraints milk yield and live weight and indirect effects globally Ticks are the obligate due to pathogen transmission and toxicosis haematophagousecto parasites of both resulting in paralysis, irritation and allergy livestock and humans. They were rank second (Aktas, 2014). after mosquitoes as vectors of human diseases. Ticks play an important role in the The direct injuries due to ixodid ticks are transmission of various pathogens like severe in tropical climate (Ghosh and Nagar, bacteria, viruses, protozoa and helminths, 2014). There is a recent surge in studies many of which have zoonotic significance (de regarding tick and tick borne diseases due to 2304 Int.J.Curr.Microbiol.App.Sci (2017) 6(8): 2304-2309 the changes in tick distribution mainly template for the polymerase chain reaction contributed by the changes in climate (PCR). Molecular confirmation of tick species (Jaensonet al., 2009). According to Burger et was done by PCR amplification of al., (2014), 904 tick species were identified mitochondrial 16S rRNA gene and V4 region globally. Almost 80 per cent of the tick of 18S rRNA gene (Crampton et al., 1996; species comes under the family Ixodidae or Kumar et al., 2011). The primer sequences for hard ticks (705 species). Family Argasidae or PCR amplification are mentioned in table 1. soft ticks contain 198 species and one under the family Nuttalliellidae. Amplified PCR products were analyzed by gel electrophoresis using 1.5 per cent agarose India is home for approximately 109 tick gel. PCR amplified products were sequenced species, under 12 genera (Ghosh et al., at Sci. Genome Pvt. Ltd., Cochin and the 2007).The climatic conditions in India are sequence data were BLAST analyzed. highly favorable for the propagation of ticks. Most predominant tick species identified in Phylogenetic analysis was performed for both India belong to the genus Rhipicephalus and 16S rRNA and 18S rRNA genesby Neighbor- Hyalomma (Ghosh et al., 2007). joining tree method using Mega 7. Nucleotide sequences were aligned by Clustal W for The prevalence of Rhipicephalus analysis. Akaike information criterion haemaphysaloides based on morphological implemented in MEGA 7.0 was used to identification was recorded from various parts determine the best fitting models. The of India including Kerala (Geevarghese and evolutionary distances were computed using Dhanda, 1995; Rajendran and Hafeez, 2003; Tamura 3-parameter method and Jukes- Prakasan and Ramani, 2007; Soundararajan et Cantor method for 16S rRNA and 18S rRNA al., 2014). However, they were no studies respectively. regarding the molecular confirmation of the tick species. In the present study, molecular Results and Discussion confirmation of R. haemaphysaloides was performed using ticks collected from Tick isolates from Pookode and Meppadi Wayanad, a northern district of Kerala. were morphologically identified as R. haemaphysaloides. Molecular analysis of tick Materials and Methods mitochondrial16S rRNA gene of Pookode and Meppadi isolates revealed 95 per cent and 93 Tick samples were collected froma bovine per cent identity respectively with R. male Calffrom College of Veterinary and haemaphysaloides (AY972534) from China. Animal Sciences, Pookode, Wayanad. Another isolate was collected from a sambar BLAST analysis of 18S rRNA sequence of deer from Meppadi, Wayanad. Ticks were Pookode and Meppadi isolate showed 96 per morphologically identified (Arthur, 1960) and cent and 100 per cent identity respectively stored in -80°C until RNA extractionwas with R. haemaphysaloides (DQ839552) from performed. RNA extraction was done using China. The sequences were submitted to Gen RNeasy Mini Kit (Qiagen, Netherlands) based Bankto obtain accession numbers. Pookode on the manufacture’s protocols. Synthesis of isolate was assigned with accession numbers, cDNA was performed by Revert Aid H minus KU895511 and KU895510 for 16S rRNA and First strand cDNA synthesis kit (Thermo 18S rRNA respectively. Similarly, Meppadi Scientific, USA). This cDNA was used as isolate was assigned with accession numbers 2305 Int.J.Curr.Microbiol.App.Sci (2017) 6(8): 2304-2309 MF351994 and MF351848 for 16S rRNA and haemaphysaloidesis usually prevalent in 18S rRNA respectively. Oriental, Australasian and Palearctic zoogeographic regions with wide host range. Phylogenetic analysis of mitochondrial16S The usual hosts of these ticks are birds and rRNA revealed that both Pookode and mammals including humans (Guglielmone et Meppadi isolates of ticks were genetically al., 2014). R. haemaphysaloides is one of the closely related. Also, both shared identity most abundant cattle ticks in Sri Lanka (Diyes with Chinese and Thailand isolates of R. and Rajakaruna, 2015). Rajendran and Hafeez haemaphysaloides. Based on evolutionary (2003) reported a prevalence of 3.29 per cent analysis of 18S rRNA, Pookode isolate of R. for R. haemaphysaloides among crossbred haemaphysaloides was more ancestral to both cattle of Andhra Pradesh. Soundararajan et Meppadi and Chinese isolates (Figs. 1 and 2). al., (2014) recorded a prevalence of 3.13 per cent in goats of Tamil Nadu. Prakasan and In the present study, the presence of R. Ramani (2007) reported R. haemaphysaloides haemaphysaloides in Wayanad was from cattle, buffalo, goat and pigs in Kerala. confirmed by molecular techniques. R. Fig.1 Molecular phylogeny based on 16S r RNA. The evolutionary history was inferred using the Neighbor-Joining method. The optimal tree with the sum of branch length = 0.56035004 is shown. The percentage of replicate trees in which the associated taxa clustered together in the bootstrap test (1000 replicates) is shown next to the branches. The tree is drawn to scale, with branch lengths in the same units as those of the evolutionary distances used to infer the phylogenetic tree. The evolutionary distances were computed using the Tamura 3-parameter method and are in the units of the number of base substitutions per site. The rate variation among sites was modeled with a gamma distribution (shape parameter = 1). The analysis involved 14 nucleotide sequences 100 KU895511 R. haemaphysaloides 16S rRNA Pookode isolate India 93 MF351994 R. haemaphysaloides 16S rRNA Meppadi isolate India KU664522 R. haemaphysaloides 16S rRNA China 97 100 KC170743 R. haemaphysaloides 16S rRNA Thailand KC170744 R. sanguineus 16S rRNA Thailand 42 KY069269 R. turanicus 16S rRNA China 100 KU183524 R. turanicus 16S rRNA China 100 KP210048 H. anatolicum excavatum 16S rRNA Haryana India KP210065 H. detritum detritum 16S rRNA Haryana India 49 KP210050 R. sanguineus 16S rRNA Haryana India 100 HM536970 R. microplus 16S rRNA Punjab India 52 GU323288 R. microplus 16S rRNA UP India KY458970 H. bispinosa 16S rRNA India 100 KC853420 H. bispinosa 16S rRNA India 2306 Int.J.Curr.Microbiol.App.Sci (2017) 6(8): 2304-2309 Fig.2 Molecular phylogeny based on 18S rRNA. The evolutionary history was inferred using the Neighbor-Joining method. The optimal tree with the sum of branch length = 1.67598566 is shown. The percentage of replicate trees in which the associated taxa clustered together in the bootstrap test (1000 replicates) is shown next to the branches. The tree is drawn to scale, with branch lengths in the same units as those of the evolutionary distances used to infer the phylogenetic tree. The evolutionary distances were computed using the Jukes-Cantor method and are in the units of the number of base substitutions per site. The analysis involved 11 nucleotide sequences