Molecular and Phylogenetic analysis revealed new genotypes of Theileria annulata parasites from India
- Neena George†1,
- Vasundhra Bhandari†1,
- D. Peddi Reddy1 and
- Paresh Sharma1Email authorView ORCID ID profile
© George et al. 2015
Received: 30 June 2015
Accepted: 4 September 2015
Published: 17 September 2015
Tick borne diseases impinge cattle worldwide causing mortality and resulting in huge economic losses. Theileriosis is one of the important tick borne diseases mainly caused by Theileria annulata and one of the commonly occurring infections among the livestock. T. annulata causes immense loss to the livestock industry and therefore, efficacious eradication and control strategies are needed for the control of the disease. Genetic diversity among T. annulata parasites is another important aspect which is overlooked in India. Thus, the present study aims to evaluate the prevalence along with genetic diversity and phylogeny of the prevailing T. annulata population of India.
Genomic DNA was extracted from cattle blood samples (n = 862) from different regions of Andhra Pradesh. Molecular diagnosis using T. annulata 18S rRNA based PCR was performed to detect parasites in cattle. Further, 18S rRNA gene was cloned and sequenced to determine similarity and diversity from the known T. annulata sequences.
We observed an overall prevalence rate of 32.40 % T. annulata infection in Andhra Pradesh based on PCR assay. The sequence analysis revealed novel genotypes among the T. annulata strains from India. Thirteen strains showed closed proximity with a strain from China whereas one Indian strain showed similarity with a South African strain [Theileria sp (buffalo)] based on phylogenetic analysis. Nucleotide heterogeneity of the 18S rRNA sequence among the strains examined varied from 0.1 to 8.6 % when compared with the published strains.
The present study provides us with the molecular prevalence of theileriosis, and will support the accomplishment of actions or in design of strategy to control theileriosis transmission to cattle. Additionally, it highlights the emergence of strains with novel genotypes from India.
Theileriosis is a tick borne disease which results in mortality and economic loss amounting to US$ 800million approximately in India . It is a lymphoproliferative disease caused by an apicomplexan intracellular parasite, Theileria annulata in cattle and water buffaloes. The clinical symptoms comprise of fever, anaemia and jaundice and the disease is fatal, if left untreated . Control of infection is mainly banking on an attenuated vaccine which is not effective and has problems like cold chain and short-term protection. An additional problem associated with the disease is high treatment cost not bearable by farmers, emerging drug resistance and lack of sufficient studies focussed on epidemiology and prevalence of theileriosis in India. During theileriosis recovery, primary infection may result in the development of a carrier state or relapse case and further, these carrier animals play a crucial role in the maintenance of the life cycle of the parasite by the ticks.
Diagnosis of the disease based on microscopic examination of Giemsa stained blood smears and lymph node biopsy is disadvantageous in case of low parasitemia and carrier state animals [3, 4]. However, the molecular based assays such as PCR have proven to be the most reliable tool for detecting Theileria parasites in clinical and subclinical cases of theileriosis especially detecting carrier state [3–7]. A number of control measures have been employed such as acaricides application, chemotherapy and vaccination, however, due to lack of knowledge on the prevalence of Theileria infection these measure are not effective [8, 9].
In addition to that, studies investigating the genotypic diversity among T. annulata parasites from India are scarce [10, 11]. The role of 18S rRNA gene in Theileria genetic diversity studies holds importance due to the presence of a conserved and a hypervariable region (V4) which are well studied and crucial for determining evolutionary patterns and similarity among the Theileria species [12, 13]. Moreover, it has been sequenced from a large number of strains belonging to different parts of the globe and database set is readily available for comparison [11, 12, 14–18]. In India, there is a need for studies focusing on the prevalence of theileriosis and to identify genetic diversity among T. annulata strains to make the control programs effective and successful. In the present study, we have randomly collected samples from the state of Andhra Pradesh (now divided into Seemandhra and Telangana) focusing on the epidemiology, molecular characterization of T. annulata parasite and to investigate genotypic variations among the strains.
The study involved drawing of ~ 5 ml blood from jugular vein aseptically from cattle with the consent of the farm owners. There is no specific law for blood sample collection and hence no approval was mandatory.
Sample collection and DNA isolation
Details of blood sample collected from different districts of Andhra Pradesh
Prevalence rate (%)
Culture of T annulata parasite and genomic DNA isolation
T. annulata infected macrophages were isolated from the clinically infected cattle and culture was maintained at 37 °C in RPMI-1640 medium supplemented with 25 mM HEPES, 4 mM L-glutamine, 100 μg/mL penicillin, 100 μg/mL streptomycin and 10 % Foetal bovine serum (FBS) with 5 % CO2 . Parasite enrichment was done from the T. annulata infected macrophages cell line using percoll gradient as per Langsley et al.  with slight modifications. DNA was isolated from the enriched parasites using QIAamp DNA Mini Kit as per manufacturer’s instructions (Qiagen).
PCR Amplification of 18S rRNA gene of T. annulata
PCR was used to amplify 125 bp fragment of T. annulata 18S rRNA gene [Genbank:EU083801.1] with species-specific primers, Forward (5′- ACGACTCCTTCAGCACCTTG-3′) and Reverse (5′-AAATTAAGCCGCAGCTCCAC-3′). PCR reactions were performed using Speed Star HS DNA polymerase (Takara) in an automatic DNA thermocycler (Biorad-T100 Thermal cycler). Reactions were performed in 25 μl volumes with 100 ng of purified genomic DNA. PCR conditions consisted of a first denaturation step of 1 min at 95 °C, followed by a second step of 35 cycles of 10 s at 95 °C, 20s at 61 °C and 10 s at 72 °C, and a final extension step of 1 min at 72 °C. DNA extracted from in vitro cultured T. annulata parasite was used as positive whereas distilled water and T. oreintalis positive DNA sample were used as negative controls, respectively in all PCR reactions. PCR products were separated by 2 % agarose gel electrophoresis and stained with ethidium bromide to assess the presence of specific bands indicative of T. annulata . Gels were photographed under gel documentation system (Syngene) and the size of each PCR product determined by reference to a DNA ladder (NEB).
Senstivity and specificity of 18S rRNA based PCR
Ten fold serial dilution of plasmid cloned 18S rRNA gene was prepared and further, utilized to evaluate the sensitivity of the PCR using conditions as described above. The highest dilution of the DNA showing a visible band was taken as the detection limit. Specificity of primer was also determined by using DNA samples positive for T. orientalis species infection.
Cloning and sequencing of 18S rRNA gene
One representative positive PCR sample (n = 14) from each district was selected for 18S rRNA gene cloning. The gene was amplified using primers (18S rRNA Forward: GGC CAG TAG TCA TAT GCT TGT and Reverse: TGA TCC TTC CGC AGG TTC ACC) and further, the amplified PCR products (base pair (bp) = 1728) were purified by Nucleospin Gel and PCR cleanup kit (Macherey Nagel, Germany) following the manufacturer’s instructions. After purification, 1728 bp PCR product was cloned into a TOPO cloning vector (Invitrogen, Life Technologies) followed by transformation into Top10 cells with ampicillin as a marker. Positive colonies were selected by PCR followed by restriction digestion with EcoR1 enzyme. Plasmid from positive colonies were isolated using Nucleospin plasmid kit (Macherey Nagel, Germany) and sequencing (five clones per sample) was performed using universal primers of M13 gene present on the backbone of the TOPO vector. Nucleotide sequence of the 18S rRNA gene was determined by DNA sequencing sent to a commercial laboratory by Sangers sequencing. Sequences were compared with 18S rRNA gene sequences of T. annulata parasites available in the NCBI database using the Basic Local Alignment Search Tool (BLAST) followed by multiple alignment using ClustalW2. The representative sequences obtained were registered in the NCBI GenBank nucleotide sequence database under the assigned accession numbers from KT367866 to KT367879 which corresponds to TA 1 – TA 14.
Phylogenetic tree and multiple alignment using 18S rRNA gene
Multiple alignment was done using ClustalW2 alignment method using 18S rRNA gene sequences. The evolutionary history was inferred by using the Maximum Likelihood method based on the Kimura 2-parameter model . The tree with the highest log likelihood (−25870.0411) is shown. The percentage of trees in which the associated taxa clustered together is shown next to the branches. Initial tree(s) for the heuristic search were obtained automatically by applying Neighbor-Join and BioNJ algorithms to a matrix of pairwise distances estimated using the Maximum Composite Likelihood (MCL) approach, and then selecting the topology with superior log likelihood value. The tree is drawn to scale, with branch lengths measured in the number of substitutions per site. The analysis involved 42 nucleotide sequences. All positions containing gaps and missing data were eliminated. There were a total of 965 positions in the final dataset. Evolutionary analyses were conducted in MEGA6 .
Divergence within the 18S rRNA gene sequenced from different districts was calculated in MEGA6 by comparing with the other available 18S rRNA sequences of T. annulata from NCBI. Divergence is calculated by comparing sequence pairs in relation to the phylogeny reconstructed by MEGA6. MEGA6 calculates the divergence between two branches by comparing sequence pairs in relation to the reconstructed phylogeny. Total Distance is the sum of all branch lengths. Distance (i,j) = sum (residue distances) + (gaps × gap penalty) + (gap residues × gap length penalty). When a direct pathway exists between two sequences on the tree, Distance (i,j) equals the sum of the branch lengths. Divergence (i,j) = 100 [Distance (i,j)]/Total Distance.
Confidence interval (CI) at 95 % was calculated using online software (https://www.mccallum-layton.co.uk/tools/statistic-calculators/confidence-interval-for-proportions-calculator).
Genetic diversity and phylogenetic analysis
There are number of hemo-protozoan parasites infections in cattle out of which T.annulata accounts for maximum infections . The disease causes economic loss due to high mortality, reduction in milk yield, emaciation as well as an increased burden due to treatment cost to the farmer . Further, impeding the problem is the lack of information available about the prevalence rate of the infection with the consequence of poor implementation of control policy.
In our study, we found high prevalence rate of 32.40 % for T. annulata infections in cattle at Andhra Pradesh by utilizing 18S rRNA based PCR. The samples were randomly collected from cattle, either cross breed or native breed from different districts of Andhra Pradesh. In India, there are few reports with less numbers showing the presence of T. annulata infection using molecular based assay for e.g., a study from Gujarat where 74 samples were detected positive out of 113 , in Bangalore 41 cases were positive out of 132  while a study from Punjab showed a prevalence rate of 14.65 % using microscopic examination . The current finding will help to shed light on the occurrence of T. annulata infection in Andhra Pradesh and bring it to the forefront of devising control policy to detect and eliminate the infection resulting in reducing the economic losses caused due to the disease.
Further, we have performed phylogenetic analysis and multiple alignment of the prevailing T. annulata strain population from Andhra Pradesh (Telangana and Seemandhra) using 18S rRNA gene. Multiple alignment revealed TA 4 & TA 5 strains with novel insertions (n = 9) and SNPs clustering with a strain from China [GenBank:KF559356] which may be due to evolutionary changes and might be important for adaptation. Further, maximum number of strains [TA 1-TA 13] showed proximity to China whereas TA14 to South African Theileria (buffalo) species [Genbank:HQ895982], however, forming a separate subclade on the evolutionary tree . The sequence similarity indicates that it might be due to bi-directional transfer or as a result of evolutionary changes. Similarity of TA 14 with South African strain raises concern since, South African strain is reported to be potentially lethal in cattle .
Reported diversity studies of T. annulata strains had considered the complete gene sequence or only the hypervariable region, V4 (~200 bp) of the 18S rRNA gene [9, 10, 12–15, 25]. In the current study, novel SNPs were found in all fourteen strains in the conserved as well as the hypervariable region indicating presence of new genotypes of T. annulata strain from India. TA 4 and TA 5 strains showed maximum changes in their nucleotide composition in comparison with the published sequence of Iran strain [GenBank:KF429795]. The study showed the presence of genotypes which cluster with the known T. annulata strains and other leukocyte transforming Theileria strains, however, forming separate sub-clades with 0.1 to a maximum 8.6 % nucleotide heterogeneity. Our data clearly show that Theileria genotypes detected are clearly distinct from earlier reported strains of India.
Overall, our study indicates the molecular prevalence of theileriosis and identification of novel genotypes of T. annulata strains in Andhra Pradesh, India.
The current study reveals a prevalence rate of about 32.40 % in the states of Andhra Pradesh and Telangana. It is a first step towards representing the occurrence of the T. annulata infection in this region. It will possibly raise the alarm for effective and urgent measures against the disease control as well as identification of novel genotypes of T. annulata in the prevailing population.
The study is part of intramural research funding from the National Institute of Animal Biotechnology, an autonomous institute of the Department of Biotechnology, New Delhi, India. We thank all the districts animal husbandry departments, veterinary doctors and local farming community for their great cooperation.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Brown CGD. Dynamics and impact of tick-borne diseases of cattle. Trop Anim Health Prod. 1997;29:1S–3.View ArticlePubMedGoogle Scholar
- Bishop R, Musoke A, Morzaria S, Gardner M, Nene V. Theileria: intracellular protozoan parasites of wild and domestic ruminants transmitted by ixodid ticks. Parasitol. 2004;129:S271–83.View ArticleGoogle Scholar
- Khatoon S, Kolte SW, Kurkure NV, Chopde NA, Jahan A. Detection of tropical bovine theileriosis by polymerase chain reaction in cattle. J Parasit Dis. 2015;39(1):53–6.View ArticlePubMedGoogle Scholar
- Kundave VR, Patel AK, Patel PV, Hasnani JJ, Joshi CG. Qualitative and quantitative assessment of Theileria annulata in cattle and buffaloes Polymerase Chain Reaction. Trop Biomed. 2014;31(4):728–35.PubMedGoogle Scholar
- Noaman. A molecular study on Theileria and Babesia in cattle from Isfahan province, Central Iran. J Parasit Dis. 2013;37(2):208–10.PubMed CentralView ArticlePubMedGoogle Scholar
- Aktas M, Dumanli N, Cetinkaya B, Cakmak A. Field evaluation of PCR in detecting Theileria annulata infection in cattle in eastern Turkey. Vet Rec. 2002;150:548–9.View ArticlePubMedGoogle Scholar
- Aktas M, Altay K, Dumanli N. A molecular survey of bovine Theileria parasites among apparently healthy cattle and with a note on the distribution of ticks in eastern Turkey. Vet Parasitol. 2006;138:179–85.View ArticlePubMedGoogle Scholar
- Beniwal RK, Nichani AK, Sharma RD, Rakha NK, Suri D. Responses in animals vaccinated with the Theileria annulata (Hisar) cell culture vaccine. Trop Anim Health Prod. 1997;29:109S–13.View ArticlePubMedGoogle Scholar
- Livestock Census, Animal Husbandry Department. http://dahd.nic.in/dahd/statistics/livestock-census.aspx. 2012.
- Manuja A, Malhotra DV, Sikka VK, Sangwan AK, Sharma R, Kumar B, et al. Isolates of Theileria annulata collected from different parts of India show phenotypic and genetic diversity. Vet Parasitol. 2006;137(3–4):242–52.View ArticlePubMedGoogle Scholar
- Aparna M, Vimalkumar MB, Varghese S, Senthilvel K, Ajithkumar KG, Raji K, et al. Phylogenetic analysis of bovine Theileria spp. isolated in south India. Trop. Biomed. 2013;30(2):281–90.Google Scholar
- Chae J, Lee J, Kwon O, Holman PJ, Waghela SD, Wagner GG. Nucleotide sequence heterogeneity in the small subunit ribosomal RNA gene variable (V4) region among and within geographic isolates of Theileria from cattle, elk and white-tailed deer. Vet Parasitol. 1998;75(1):41–52.View ArticlePubMedGoogle Scholar
- Sivakumar T, Hayashida K, Sugimoto C, Yokoyama N. Evolution and genetic diversity of Theileria. Infect Genet Evol. 2014;27:250–63.View ArticlePubMedGoogle Scholar
- Chansiri K, Kawazu S, Kamio T, Terada Y, Fujisaki K, Philippe H, et al. Molecular phylogenetic studies on Theileria parasites based on small subunit ribosomal RNA gene sequences. Vet Parasitol. 1999;83(2):99–105.View ArticlePubMedGoogle Scholar
- Aktas M, Bendele KG, Altay K, Dumanli N, Tsuji M, Holman PJ. Sequence polymorphism in the ribosomal DNA internal transcribed spacers differs among Theileria species. Vet Parasitol. 2007;147:221–30.View ArticlePubMedGoogle Scholar
- Liu Q, Zhou YQ, He GS, Oosthuizen MC, Zhou DN, Zhao JL. Molecular phylogenetic studies on Theileria spp. isolates (China) based on small subunit ribosomal RNA genesequences. Trop Anim Health Prod. 2010;42(1):109–14.View ArticlePubMedGoogle Scholar
- Habibi G. Phylogenetic analysis of Theileria annulata infected cell line S15 Iran vaccine strain. Iran J Parasitol. 2012;7(2):73–81.PubMed CentralPubMedGoogle Scholar
- Hooge J, Howe L, Ezenwa VO. Identification of novel Theileria genotypes from Grant’s gazelle. Int J Parasitol Parasites Wildl. 2015;4(2):239–43.PubMed CentralView ArticlePubMedGoogle Scholar
- Baumgartner M, Tardieux I, Ohayon H, Gounon P, Langsley G. The use of nocodazole in cell cycle analysis and parasite purification from Theileria parva-infected B cells. Microb Infect. 1999;1181–88.Google Scholar
- Kimura M. A simple method for estimating evolutionary rate of base substitutions through comparative studies of nucleotide sequences. J Mol Evol. 1980;16(2):111–20.View ArticlePubMedGoogle Scholar
- Tamura K, Stecher G, Peterson D, Filipski A, Kumar S. MEGA6: Molecular Evolutionary Genetics Analysis version 6.0. Mol Biol Evol. 2013;30:2725–9.PubMed CentralView ArticlePubMedGoogle Scholar
- Singh NK, Jyoti S, Haque M, Rath SS. Prevalence of parasitic infections in cattle of Ludhiana district, Punjab. J Parasit Dis. 2012;36(2):256–9.PubMed CentralView ArticlePubMedGoogle Scholar
- Ananda KJ, D’Souza PE, Putta lakshmamma GC. Prevalence of Haemoprotozoan diseases in crossbred cattle in Bangalore north. Vet Wor. 2009;2:5–16.Google Scholar
- Chaisi ME, Sibeko KP, Collins NE, Potgieter FT, Oosthuizen MC. Identification of Theileria parva and Theileria sp. (buffalo) 18S rRNA gene sequence variants in the African Buffalo (Syncerus caffer) in southern Africa. Vet Parasitol. 2011;182:150–62.View ArticlePubMedGoogle Scholar
- Khan MK, He L, Hussain A, Azam S, Zhang WJ, Wang LX, et al. Molecular epidemiology of Theileria annulata and identification of 18S rRNA gene and ITS regions sequences variants in apparently healthy buffaloes and cattle in Pakistan. Infect Genet Evol. 2013;13:124–32.View ArticlePubMedGoogle Scholar