Int.J.Curr.Microbiol.App.Sci (2018) 7(7): 3351-3358

International Journal of Current Microbiology and Applied Sciences ISSN: 2319-7706 Volume 7 Number 07 (2018) Journal homepage: http://www.ijcmas.com

Original Research Article https://doi.org/10.20546/ijcmas.2018.707.390

Molecular Characterization of β-Tubulin Isotype-1 Gene of trigonocephalum

Ravi Kumar Khare1, A. Dixit3, G. Das4, A. Kumar1, K. Rinesh3, D.S. Khare4, D. Bhinsara1, Mohar Singh2, B.C. Parthasarathi2, P. Dipali2, M. Shakya5, J. Jayraw5, D. Chandra2 and M. Sankar1*

1Division of Temperate Husbandry, ICAR- IVRI, Mukteswar, India 2IVRI, Izatnagar, India 3College of Veterinary Science and A.H., Rewa, India 4College of Veterinary Sciences and A.H., Jabalpur, India 5College of Veterinary Sciences and A.H., Mhow, India

*Corresponding author

ABSTRACT

The mechanism of benzimidazoles resistance is linked to single nucleotide polymorphisms (SNPs) on beta -tubulin isotype-1 gene. The three known SNPs responsible for BZ K e yw or ds resistance are F200Y, F167Y and E198A on the beta-tubulin isotype-1. The present study

was aimed to characterize beta-tubulin isotype-1 gene of Bunostomum trigonocephalum, Benzimidazole for identifying variations on possible mutation sites. The adult parasites were collected resistance, Beta from Mukteswar, Uttarakhand. The parasites were thoroughly examined morphologically tubulin, and male parasites were subjected for RNA isolation. Complementary DNA (cDNA) was Bunostomum synthesised from total RNA using OdT. The PCR was performed using cDNA and self trigonocephalum, Small designed degenerative primers. The purified PCR amplicons were cloned into pGEMT easy vector and custom sequenced. The obtained sequences were analysed using DNA Article Info STAR, MEGA7.0 and Gene tool software. The deduced amino acid sequence showed 99> % homology with published B. trigonocephalum as well as closely related Accepted: Ancylostomum caninum and Strongyles of equines. It is also showed 98-98.7% identity 24 June 2018 Available Online: with Trichostrongylus species and 91.8-93.6% homology with other helminths like P.equorum, A. galli and F.hepatica. The information obtained from current study enlighten 10 July 2018 to investigate further related to benzimidazole resistance in B. trigonocephalum and

formulate effective control strategy as this parasite is one of the most pathogenic strongyles of small ruminants.

Introduction ruminants and the infection, bunostomiasis is characterised by anaemia due to blood sucking Bunostomum trigonocephalum (Order: of worm and dermatitis by larval penetration, , Family: ) is particularly lower part of infected host commonly known as hookworm of small (Soulsby, 1982). Few hundreds of worm can

3351

Int.J.Curr.Microbiol.App.Sci (2018) 7(7): 3351-3358 kill an animal (Soulsby, 1982) and characterize β –tubulin gene of Bunostomiasis is reported to be affecting all B.trigonocephalum for analyzing and age groups, mainly of growing young ones (5- predicting mutation pattern with respect to BZ 8 months aged).The infection is more resistance. With this aim, the present study prevalent in warm and humid regions (Tariq et was planned to characterize beta-tubulin al., 2008, 2010), and is also reported as major isotype 1 gene of B.trigonocephalum of cause of economic losses in the livestock Mukteswar. industry in temperate areas (Stancampiano, 2007). The prevalence of Bunostomum spp. is Materials and Methods restricted in few pockets of India, mainly from central(Singh et al., 2016; Rajpoot et al., Study area and collection of parasites 2017), north east (Yadav and Tondon, 1989; Bandyopadhyay et al., 2010) and Kashmir Adult Bunostomum trigonocephalum isolate valley (Tariq et al., 2008, 2010). The were collected from gastrointestinal tract of prevalence of B. trigonocephalumis very high, goats slaughtered at local abattoir at especially in Kashmir valley, where Mukteswar (29°28’N and 79°39’E, 7500 feet prevalence in sheep and goat were 37.7% and above mean sea level), Uttarakhand. Parasites 30.1%, respectively (Tariq et al., 2008, 2010). were washed thoroughly in PBS (pH 7.4) and Control of gastrointestinal nematodiasis identified as per the morphological keys includingbunostomiasis has achieved by using (Johnson, 1965; Soulsby, 1982). The adult broad spectrum chemotherapeutic agents like male worms were used for extraction of total benzimidazoles (BZs), imidazothiazoles, RNA. tetrahydropyrimidines and macrocyclic lactones. The excessive and frequent use of Isolation of total RNA and cDNA synthesis anthelmintics has resulted in substantial and widespread emergence of anthelminthic Total RNA was isolated from adult male B. resistance (AR), particularly against BZin trigonbocephalum using RNeasy minikit populations (Kaplan et al., 2004; (Qiagen, Germany) as per manufacturer’s Garg and Yadav, 2009; Chandra et al., 2014, instructions. The complementary DNA 2015; Dixit et al., 2017). Maximum reports of (cDNA) was synthesized from the total RNA BZ resistance are restricted in three main of adult male B.trigonocehalum using oligo gastrointestinal nematodes Haemonchus dT primer and by using RevertAid reverse contortus. Trichostrongylus colubriformis and transcriptase enzyme (Thermo scientific, Teladorsagia circumcincta (Kwa et al., 1993, USA). 1994, 1995; Silvestre and Humbert, 2000; Ghisi et al., 2007; Rufener et al., 2009; Garg Polymerase chain reaction for amplification and Yadav, 2009; Chandra et al., 2014, 2015; of β-tubulin isotype-1 gene Dixit et al., 2017). BZ resistance is primarily linked to a point mutation at amino acid 200 PCR was standardized to amplify the β-tubulin (Phe to Tyr) (Kwa et al., 1993, 1994, isotype-1 gene of B.trigonocephalum from 1995),167 (Phe to Tyr) (Ghisi et al., 2007) and cDNA. The open reading frame of truncated 198 (Glu to Ala) (Rufener et al., 2009, β-tubulin gene was amplified using the self- Chaudhary et al., 2015) of β-tubulin isotope-1 designed degenerative primers (forward gene. However, works on Bunostomum genus, primer 5’GCC GGW CAR TGC GGH AAC particularly on the B.trigonocephalum is CAG 3’ and reverse primer 5’GTG AAY meager. Therefore, it is necessary to TCC ATT TCG TCC ATA C 3’) and were

3352

Int.J.Curr.Microbiol.App.Sci (2018) 7(7): 3351-3358 designed to amplify all the expected mutations B.trigonocephalum and other related for BZ resistance present in the gene such as Strongylus species. The β-tubulin gene 167th, 198th and 200th position. The PCR sequences of other Strongylus and other mixture consisted of cDNA as template 1.0µl, helminthes were retrieved from NCBI TerraTM PCR Direct Red Dye Premix 12.5 µl, database and used for comparative analysis 10 pmoles of each primer and adding nuclease purpose. water to make final volume 25 µl. The reaction was standardized with annealing Results and Discussion temperature at 60ºC. The amplicons were confirmed in 1.2% agarose gel The anterior end of B.trigonocephalum is bent electrophoresis. in a dorsal direction; therefore the parasite is looks like hooks. The buccal capsule is Molecular cloning, sequencing and triangular funnel shape opens anterio-dorsally; characterization of β-tubulin gene of it has a large dorsal tooth and two short B.trigonocephalum ventral teeth. There are two sub-ventral cutting plates and small pair of dorsal plates The amplicons were gel purified using near moth opening. Large dorsalcone is Qiaquick Gel extraction kit (Qiagen, characteristics of this species, which projects Germany) and ligated with 50ng of pGEM®-T in to the ⅔ of buccal cavity (Fig.1). The bursa easy TA cloning vector (Promega) in 1:3 of male B.trigonocephalumis well developed ratios (Vector: Purified amplicons). The with small asymmetrical dorsal lobe, which is recombinant plasmid was transformed in to not well demarcated from lateral lobes. The E.coli Top10 competent cells by heat-shock spicules are spirally twisted and united method at 420C for 90 sec. The transformed posteriorly. The spicules are 0.6-0.64 mm culture was plated over the freshly prepared long, slender and alate. The gubernaculums is LB Amp+ X gal+ IPTG+ plates and incubated absent (Fig.2). The morphological features of overnight at 37°C. The positive colonies were B.trigonocephalum are documented elsewhere selected using blue-white screening method (Johnson, 1965; Soulsby, 1982; Suresh Singh, (α-complementation) and further confirmed by 2003) colony PCR and release of desired products from vector using EcoRI enzyme by restriction Amplification of β-tubulin gene sequence of enzyme digestion. Subsequently, the positive B.trigonocephalum clone was inoculated in LB stab culture and custom sequenced. The PCR was amplified approximately 1178 bp size fragment of β-tubulin isotype-1 gene Genetic characterization in agarose gel electrophoresis (Fig.3). The PCR product was purified and the Stab cultures of positive clones harbouring the concentration of purified β-tubulin gene was desired β-tubulin gene was sent for custom 32ng/µl. DNA sequencing to Department of Biochemistry, Delhi University, South The ligated amplicon with pGEM®-T easy campus. The sequence information received TA cloning vector was successfully was analyzed by using ClustalW pair distance transformed as evidenced by appearance of method (DNA Star) and phylogenetic tree was desired white colonies in the LB Amp+ X gal+ constructed using maximum likelihood IPTG+ plates and by colony PCR (Fig.4). The method (MEGA version 7.0) with published presence of insert was further confirmed by beta tubulin isotype 1 gene of restriction enzyme analysis (Fig.5). 3353

Int.J.Curr.Microbiol.App.Sci (2018) 7(7): 3351-3358

Sequencing and genetic characterization of Phylogenetic tree was constructed with β-tubulin isotype-1 gene of B. Maximum Likelihood (ML) method using trigonocephalum Tamura-Nei model with 500 bootstrap replications to confirm the authenticity of the The positive clones harboring β-tubulin gene taxa analysed for each node. Input file was of B. trigonocephalum were custom sequenced obtained by applying the BioNJ method to a and analysis result revealed that the amplicon matrix of pairwise distances estimated using size is 1178 bp. The deduced amino acid of the Maximum Composite Likelihood (MCL) sequenced region of the gene consisted of all approach. the possible and reported mutation sites for BZ resistance i.e. F167Y, E198A and F200Y. A discrete Gamma distribution was used to Mutation at amino acid 200 of the beta tubulin model evolutionary rate differences and isotype-1 (Phe to Tyr) is mostly responsible number of substitutions per site (Fig.7). For for BZ resistance in H.contortus (Kwa et al., understanding the real situation of BZ 1994; Rufener et al., 2009). However, resistance at field level, characterization of mutations at 167 (Phe to Tyr) and 198 (Glu to beta tubulin gene is pre-requisite. Ala) are also reported to be associated with Characterization studies on target genes resistance in some isolates of H. contortus enable to identify their polymorphism, if any, (Prichard, 2001; Ghisi et al., 2007; Rufener et which may provide great platform in al., 2009). formulating effective control strategies. The phylogenetic analysis of deduced amino acid The characterization of beta tubulin isotype-1 sequences revealed all the isolates of gene of B. trigonocephalum revealed that the B.trigonocephalum clustered in one clade with organism is susceptible to benzimidazole boot value more than 500and other strongyles resistance as amino acid on 167 and 200th were in another clade.. position is phenylalanine, and 198th position is glutamic acid. The sequenced genes were As expected B. trigonocephalum are much aligned and analysed with available β- tubulin closed associated with strongyles of same gene sequence of B.trigonocephalum and super family such as A. canimum and also other Strongylus sequences (Fig.6). The beta- with horse strongyles β-tubulin isotype 1 tubulin isotype-1 gene of B. trigonocephalum sequences. The analysis also suggested that β- has >99% with published sequence of same tubulin isotype 1 gene sequences of B. organism and much closed related parasites trigonocephalum isolates have closely related A.caninum and Cyathostomes. with Trichostrongyles and other helminthes clustered separately. The identity was 98-98.7% with Trichostrongylus species such as H.contortus, The study is concluded that beta tubulin T.colubriformis, Cooperia pectinata and isotype-1 gene of B.trigonocephalum and C.oncophora. Further, 91.8-93.6% homology other trichostrongyles are highly conserved. with other species of helminths like The information of beta tubulin isotype-1 gene P.equorum, A. galli and F.hepatica. Since the of B.trigonocephalum provided an idea for β-tubulin is one of the framework proteins of development of molecular tools to diagnosis the cell function, many of the amino acid will of benzimidazole resistance at the early stage be conserved. Therefore, single mutation in the country. creates considerable functional consequences.

3354

Int.J.Curr.Microbiol.App.Sci (2018) 7(7): 3351-3358

Fig.1 Anterior end of B.trigonocephalum Fig.2 Posterior end of B.trigonocephalum showing well showing spirally twisted

Fig.3 PCR amplification of β-tubulin gene of Fig.4 Colony PCR amplification of β-tubulin B.trigonocephalum Lane M: 1kb DNA ladder gene of B.trigonocephalum Lane M: 1kb DNA Lane 1&2: β- tubulin gene ladder Lane 1&2: β- tubulin gene 1 2 M M 1 2

Fig.5 Release of insert from recombinant clone by EcoRI enzyme Lane M: 1kb DNA ladder Lane C1-C3: β- tubulin gene insert C1 C2 C3 M

3355

Int.J.Curr.Microbiol.App.Sci (2018) 7(7): 3351-3358

Fig.6 Pairwise relationship between B. trigonocephalum with other helminths

Fig.7 Molecular phylogenetic analysis by Maximum Likelihood method

100 B.trigonocephalum/Mukteswar/Ind 82 B.trigonocephalum/Ind

51 A.caninum/USA Cylicocyclus_nassatus/U.K 100 100 C.longibursatus/U.K 91 C.goldi/U.K 82 O.columbianum/Mukteswar/Ind T.colubriformis/France C.oncophora/Canada 99 100 46 H.contortus/Mukteswar/Ind 100 H.contortus/Switz Cylicocyclus)_nassatus/Scot Parascaris_equorum/Sweden 98 A.galli/Sweden Fasciola_hepatica/U.K

0.1

3356

Int.J.Curr.Microbiol.App.Sci (2018) 7(7): 3351-3358

Acknowledgements 1478. Garg, and Yadav. 2009. Genotyping of The authors are highly thankful to Indian benzimidazole susceptible and Council of Agricultural Research (ICAR) resistant alleles in different New Delhi, India for funding through network populations of Haemonchus contortus programme on gastrointestinal , from Himalayan and sub-Himalayan Project Coordinator of Network Programme regions of North West India. Trop. on GI Parasitism and Director of Indian Anim. Health Prod. 41 (7): 1127– Veterinary Research Institute, Izatnagar, India 1131. for providing facilities for research. Ghisi M., Kaminsky R. and Mäser P. 2007. Phenotyping and genotyping of References Haemonchus contortus isolates reveals a new putative candidate mutation for Bandyopadhyay, S., Devi, P., Bera, A., benzimidazole resistance in Bandyopadhyay, S. and Bhattacharya, nematodes. Vet. Parasitol. 144(3): D., 2010. Prevalence of 313–320. gastrointestinal parasite in goats in Johnson, S. 1965. Observations on the Shillong, Meghalaya, India. systematic position of the genus Webmedcentral.1(9):WMC00777. Gaigeria Railliet & Henry, 1910 Borgsteede, F.H.M., Tibben, J., Cornelissen, (subfamily Necatorinae Lane, 1907: J.B., Agneessens, J. and Gaasenbeek, family Ancylostomidae Lane, 1907: C.P., 2000. Nematode parasites of Nematoda). Biol. Trop., 13(2): 213- adult dairy cattle in the Netherlands. 217. Vet. Parasitol. 89: 287-296. Kaplan, R.M., Burke, J.M., Terrill, T.H., Chandra, S., Prasad, A., Sankar, M., Yadav, Miller, J.E., Getz, W.R., Mobini, S., N. and Dalal, S. 2014. Molecular Valencia, E., Williams, M.J., diagnosis of benzimidazole Williamson, L.H., Larsen, M., and resistance in in sheep from different Vatta, A.F., 2004. Validation of the geographic regions of North India FAMACHA eye color chart for Haemonchus contortus. Veterinary detecting clinical anemia in sheep and World. 7: 337-341. goats on farms in the southern United Chandra, S., Prasad, A., Yadav, N., States. Vet. Parasitol. 123:105-120. Latchumikanthan, A., Rakesh, R.L., Kwa, M. S., Kooyman, F. N. J., Boersema, J. Praveen, K., Khobra, V, Subramani, H., and Roos, M. H. 1993. Effect of K.V., Misri, J. and Sankar, M. 2015. selection for benzimidazole resistance Status of benzimidazole resistance in in Haemonchus contortus in β-tubulin Haemonchus contortus of goats from isotype-1 and isotype-2 genes. different geographic regions of Uttar Biochem. Biophys. Res. Commun. Pradesh, India. Vet Parasitol. 208 191: 413-419. (34):263-267. Kwa, M. S., Veenstra, J. G. and Roos, M.H. Dixit, A.K., Das, G., Dixit, P., Singh, A.P., 1994. Benzimidazole resistance in Kumbhkar, N.K., Sankar, M. and Haemonchus contortus is correlated Sharma, R.L. 2017. An assessment of with a conserved mutation at amino benzimidazole resistance against acid 200 in β-tubulin isotype–1. Mol. caprine nematodes in Central India. Biochem. Parasitol. 63: 299-303. Trop Anim Health Prod. 49: 1471- Kwa, M. S., Veenstra, J. G., and Roos, M. H.

3357

Int.J.Curr.Microbiol.App.Sci (2018) 7(7): 3351-3358

1993. Molecular characterization of β- parasites. Exp. Parasitol. 95: 271-276. tubulin genes present in Singh, A.K., Das, G., Roy, B., Nath, S., and benzimidazole resistant populations of Naresh, R., 2016. Prevalence of Haemonchus contortus. Mol. Strongyle infections in goat of Maha Biochem. Parasitol. 60: 133-144. Koushal region, Madhya Pradesh, Kwa, M. S., Veenstra, J. G., Dijk, M. D. and India. J Parasit Dis.40 (2):289-291. Roos, M. H. 1995. β-tubulin genes Soulsby, E. J. L. 1982. Helminths, Arthropods from the parasite nematode and Protozoa of domesticated . Haemonchus contortus modulate drug 7th edn. Billiere resistance in Caenorhabditis elegans. Stancampiano, L. 2007. Parasites of the J. Mol. Biol. 246: 500-510. digestive tract in beef cattle imported Prichard, R. K. 2001. Genetic variability from France to Italy. Parassitologia following selection of Haemonchus 49(1-2):101-106. contortus with anthelmintics. Trends Suresh Singh, K.R. 2003. Veterinary in Parasitol. 17: 445-453. Helminthology, Book Published by Rajpoot, J., Shukla, S., Jatav, G.P., Garg, Indian Council for Agricultural U.K. and Agrawal, V., 2017. Research, New Delhi.pp379-400. Coproculture study of strongyle Tariq, K.A., Chishti, M.Z., Ahmad, F. and infection of goats from Malwa region Shawl, A.S. 2008. Epidemiology of of Madhya Pradesh.J.Entomol. gastrointestinal nematodes of sheep Zoology studies.5 (5): 876-878. managed under traditional husbandry Rufener L., Kaminsky R. and Maser P. 2009. system in Kashmir valley. Vet. In vitro selection of Haemonchus Parasitol. 158:138–143. contortus for benzimidazole resistance Tariq, K.A., Chishti, M.Z. and Ahmad, F., reveals a mutation at amino acid 198 2010. Gastro-intestinal nematode of beta-tubulin. Mol. Biochem. infections in goats relative to season, Parasitol. 168(1): 120–122. host sex and age from the Kashmir Silvestre, A. and Humbert, J.F. 2000. A valley, India. J. Helminthol. 84:93–97. molecular tool for species Yadav, A.K., Tandon, V., 1989. identification and benzimidazole Gastrointestinal nematode infection of resistance diagnosis in larval goat in a subtropical and humid zone communities of small of India. Vet.Parasitol. 33(2):135-142.

How to cite this article:

Ravi Kumar Khare, A. Dixit, G. Das, A. Kumar, K. Rinesh, D.S. Khare, D. Bhinsara, Mohar Singh, B.C. Parthasarathi, P. Dipali, M. Shakya, J. Jayraw, D. Chandra and Sankar, M. 2018. Molecular Characterization of β-Tubulin Isotype-1 Gene of Bunostomum trigonocephalum. Int.J.Curr.Microbiol.App.Sci. 7(07): 3351-3358. doi: https://doi.org/10.20546/ijcmas.2018.707.390

3358