- Research article
- Open Access
First Report on Circulation of Echinococcus ortleppi in the one Humped Camel (Camelus dromedaries), Sudan
BMC Veterinary Researchvolume 9, Article number: 127 (2013)
Echinococcus granulosus (EG) complex, the cause of cystic echinococcosis (CE), infects humans and several other animal species worldwide and hence the disease is of public health importance. Ten genetic variants, or genotypes designated as (G1-G10), are distributed worldwide based on genetic diversity. The objective of this study was to provide some sequence data and phylogeny of EG isolates recovered from the Sudanese one-humped camel (Camelus dromedaries). Fifty samples of hydatid cysts were collected from the one- humped camels (Camelus dromedaries) at Taboul slaughter house, central Sudan. DNAs were extracted from protoscolices and/or associated germinal layers of hydatid cysts using a commercial kit. The mitochondrial NADH dehydrogenase subunit 1 (NADH1) gene and the cytochrome C oxidase subunit 1 (cox1) gene were used as targets for polymerase chain reaction (PCR) amplification. The PCR products were purified and partial sequences were generated. Sequences were further examined by sequence analysis and subsequent phylogeny to compare these sequences to those from known strains of EG circulating globally.
The identity of the PCR products were confirmed as NADH1 and cox1 nucleotide sequences using the Basic Local Alignment Search Tool (BLAST) of NCBI (National Center for Biotechnology Information, Bethesda, MD). The phylogenetic analysis showed that 98% (n = 49) of the isolates clustered with Echinococcus canadensis genotype 6 (G6), whereas only one isolate (2%) clustered with Echinococcus ortleppi (G5).
This investigation expands on the existing sequence data generated from EG isolates recovered from camel in the Sudan. The circulation of the cattle genotype (G5) in the one-humped camel is reported here for the first time.
The larval stage of Echinococcus granulosus complex (EG) causes cystic echinococcosis (CE) in domestic livestock and humans. Because of the involvement of the vital organs, CE in humans is considered a critical public health problem. In addition to human health concerns, infections in cattle breeding areas may result in economic losses . Moreover, CE represents one of the neglected tropical diseases, especially in the Sub- Saharan Africa . In the Sudan, many reports of cystic echinococcosis have been described in humans and animals [3–12]. In addition to its importance as a major public health problem in the country, CE is also considered as one of the major causes of condemnation of sheep carcasses during meat inspection . Currently, 10 distinct genotypes of EG designated as G1-G10 have been described worldwide on the basis of genetic diversity related to nucleotide sequences of the mitochondrial cytochrome C oxidase subunit 1 (cox1) and NADH dehydrogenase subunit 1 (NADH 1) genes. These different genotypes are associated with distinct intermediate hosts including sheep, pigs, cattle, horses, camels, goats and cervids [14–25]. Of the ten genotypes of EG, the sheep (G1), the cattle (G5) and the camel (G6) strains were reported in humans and livestock in the Sudan [8–10, 26–28]. Recent epidemiological studies indicated that the camel strain (G6) was the most prevalent strain in Sudan [10, 27, 28]. EG complex comprises a number of intra-specific variants, strains or genotypes at the genetic level [29–32]. It was suggested that the extensive intra-specific genetic variation of E. granulosus could be better understood within the context of variations in the life cycle pattern . In addition, different genotypes would probably exhibit different antigenicity, transmission profiles, and sensitivity to chemotherapeutic agents as well as different pathological consequences [29, 33, 34]. These biological variations should be considered in developing vaccines, diagnostic kits and pharmacological therapies for control of CE. Therefore, due to epidemiological implementation and control strategies, it is essential that circulating EG genotypes in a given area of endemicity should be clearly defined [30, 35]. In this study, hydatid cysts recovered from the one- humped camel (Camelus dromedaries) in an endemic area of Tamboul, Central Sudan, were defined by genetic studies and subsequent phylogenetic analysis. The molecular characterization was made possible by targeting fragments of the mitochondrial NADH 1 and cox1 subunit 1 genes to define the circulating genetic variants in the hyper endemic area of Tmboul, Central Sudan.
Collection of samples
Hydatid cysts (n = 50) were collected over a period of one month from camel at the slaughterhouse of Tamboul, a village located at the camel producing region of Central Sudan. This slaughterhouse represents one of the major abattoirs of camel in the Sudan. The hydatid cysts were obtained from camels soon after slaughtering and transferred in thermo- flasks to the Molecular Biology Laboratory at the Faculty of Veterinary Medicine, University of Khartoum, for processing and molecular characterization studies. The cysts were mainly located in the lung (n = 45) and the rest were found in the liver (n = 5). All cysts were fertile and measured 2–10 cm in diameter. Hydatid cysts containing protoscolices and associated germinal layers were aspirated with sterile needles. The aspirates were transferred to clean sterile 50 ml tubes to which 70% alcohol was added as preservative.
DNA Extraction from intact cysts
The suspensions containing protoscolices and/or associated germinal layers were washed in nucleic acid free water to remove excess alcohol. Extraction of DNA from hydatid cysts was made possible using a commercially available QIAamp tissue kit (QIAGEN, Hilden, Germany) according to the manufacturer’s instructions. Briefly, 200 μl of the suspended aspirate, 20 μl of proteinase K stock solution, and 200 μl of lysing buffer were pipetted into 1.5 ml eppendorf tube. The mixture was incubated at 37°C for 1 h and then at 70°C for 30 min before the addition of 200 μl of absolute alcohol and mixing by vortexing. The mixture was then transferred to the QIAamp spin column placed in a clean 2 ml collection tube and centrifuged at 8000 RPM in MiniSpin centrifuge (Eppendorf, Wesseling-Berzdorf, Germany) for 1 min at room temperature. The QIAamp spin column was washed twice with 500 μl of the washing buffers by spinning for 1 min. The QIAamp spin column was placed in a clean 1.5 ml eppendorf tube and the DNA was eluted with 200 μl of double distilled water preheated at 70°C. Maximum DNA yield was obtained by spinning at 12,000 RPM for 1 min at room temperature. From the suspended nucleic acid 5 μl was used in the PCR amplification.
Selection of primers
The primers were selected from mitochondrial NADH dehydrogenase subunit 1 (NADH 1) gene and cytochrome C oxidase subunit 1 gene (cox1) gene. The NADH 1 primers used in this study were basically described previously . For the first amplification step, a pair of outer primers JBl1: 5' AGATTCGTAAGGGGCCTAATA 3' and JB12: 5' ACCACTAACTAATTCACTTTC 3' were used to amplify a 530 bp PCR product from EG isolates. A pair of internal sequencing primers JB11.5: 5' TTATGGTAGATATTATAG 3' and JB12.5: 5' CACACACATAAAACAAGC 3' designed to conserved segments of the EG NADH 1 sequences, were used to generate a 471 bp PCR product. The cox1 primers used in this study include a forward primer CO1: 5' GAG GTT TAT TTT TTT GGG CAT CCT 3' and a reverse primer CO1: 5' TAA AGA AAG ATA ATG AAA ATG 3' as described before . CO1 and CO2 primers would produce a 415 bp specific PCR product.
Polymerase chain reaction
A stock buffered solution containing 150 μl 10x PCR buffer, 100 μl of 25 mM MgCl2, 12.5 μl of each dATP, dTTP, dGTP and dCTP at a concentration of 10 mM was prepared in 1.5 ml eppendorf tube. The primers were used at a concentration of 20 pg/μl, and double distilled water was added to bring the volume of the stock buffer solution to 1.5 ml. Each 0.5 ml PCR reaction tube contained 2 μl of the primers, 1 μl (5.0 U) of Taq DNA polymerase (QIAGEN), 5.0 μl of the target DNA and 42 μl of the stock buffered solution. For nested PCR, 2 μl of the primary PCR product was used as DNA template. The thermal cycling profiles were as follows: a 2 min initial incubation at 95°C, followed by 40 cycles of 95°C for 1 min, 55°C for 30 sec and 72°C for 45 sec, and a final incubation at 72°C for 10 min. Thermal profiles were performed on a Techne TC-412 thermal cycler (Techne, Staffordshire, UK). Following amplification, 15 μl from each PCR containing amplified products were loaded onto gels of 1.0% agarose and electrophoresed for 1 h. The gels were stained with ethidium bromide and the PCR products were easily identified following visualization under UV light.
Sequence analysis and construction of phylogenetic tree
The primary PCR products were purified using QIAquick PCR purification kit (QIAGEN) and sent for sequencing to commercial companies (Seqlab, Göttingen, Germany, and Macrogen, Seoul, Korea). Resulted sequences were edited and aligned using BioEdit software (Ibis Biosciences, Carlsbad, CA, USA). The Basic Local Alignment Search Tool (BLAST) of NCBI (National Center for Biotechnology Information, Bethesda, MD, USA) was used to confirm the identity of the generated sequences in relation to the GenBank nucleotide database. The sequences were then aligned with the corresponding regions of NADH 1 and cox1 subunit genes of known genotypes from other countries. The phylogenetic trees were constructed using Unweighted Pair Group Method with Arithmetic mean (UPGMA) implemented in MEGA software version 5.0 . The country of origin, the GenBank accession numbers and the genotype were given for each EG isolate when available. Bootstrap analysis of 1000 replicates was applied and values were given at relevant nodes of the constructed tree. Corresponding nucleotide sequences of NADH 1 and cox 1 of Taenia multiceps with GenBank accession numbers HM143887.1 and JX535576.1, respectively; were used as out groups in the constructed phylogenetic trees.
Polymerase chain reaction (PCR)
The PCR-based assay with primers specific for NADH1 and cox1 genes yielded amplification products from of all of the fifty hydatid cysts obtained from naturally infected camels. For NADH1, the outer pair of primers produced a primary 530 bp PCR product and the nested primers produced a 471 bp PCR product. Cox 1 primers produced 415 bp PCR products.
Sequence analysis and phylogenetic relationship
The sequences obtained from the PCR products were found to align with corresponding regions for NADH1 and cox1 genes in the GenBank confirming the cysts to contain the EG complex. Aligned with BioEdit, partial sequences for NADH 1 and cox1 showed 100% homology among 49 out of 50 EG isolates recovered in this study. To investigate for the relationship between these EG isolates and the other EG genotypes identified globally, phylogenetic trees were constructed (Figures 1 and 2). Forty nine (98%) of EG isolates (represented by one sequence in the tree) clustered with E. canadensis genotype 6 (G6), the camel genotype, of EG complex obtained from other parts of the world with a strong bootstrap (1000 replicates). However, only one EG isolate (2%) clustered with E. ortleppi, the cattle genotype (G5). The partial sequences of NADH 1 gene representing the camel genotype (G6) and of the cattle genotype (G5) were submitted to the GenBank under accession numbers JN637176, JN637177, respectively. The partial sequences generated from mitochondrial cytochrome C oxidase subunit 1 (cox1) for G6 and G5 strains were submitted to the GenBank under accession numbers KC417036, JX912709, respectively.
In Sudan, camels are owned by migratory pastoralists as a source of milk, meat, riding animals, as well as sign of wealth. Extensive research has been conducted to evaluate the role played by camels in transmission of parasitic infections with special emphasis on cystic echinococosis [9, 10, 37]. The camel strain (G6) was reported to be the most prevalent genotype of EG in the Sudan, and that camels seem to play an important role in the transmission cycle of the parasite and the epidemiology of the disease. In contrast, the majority of the Sudanese ecotypes of desert sheep and Nubian goats seem to harbor calcified or infertile cysts of EG-complex [9, 10, 38]. It was, therefore, suggested that sheep and goats have natural resistance to infection with EG-complex. However, this assumption requires further investigation. To advance beyond the current knowledge of the epidemiology of the disease in camel, attempts were made to better understand the life cycle and molecular epidemiology of EG isolates circulating in different parts of the Sudan. An initial step in controlling the life cycle of EG and minimizing infections is to determine the genotype(s) involved in the transmission cycle. In a previous study, different genotypes including G5 and G6 were confirmed in the Sudan by partial sequence analysis of the mitochondrial cytochrome oxidase 1 (cox1) and NADH dehydrogenase subunit 1 genes . However, no phylogenetic study was conducted to determine the phylogenetic relationship of the identified genotypes of EG in the Sudan. In this study, the phylogenetic analysis illustrated that Echinococcus canadensis genotype 6 (G6) is the most infectious and widespread genotype in the Sudan, confirming the results of our previous studies [8, 27]. Nevertheless, the present study indicated that Echinococcus ortleppi, the cattle genotype 5 (G5) should equally be considered as an infectious form of EG-complex in the one humped camels in Sudan. Sequence analysis and subsequent phylogenetic studies of fragments of the mitochondrial NADH subunit 1 gene indicated the presence of a single isolate of the cattle genotype (G5). However, the remaining 49 (98%) isolates belonged to E. canadensis (G6). Similar results were obtained when constructing phylogenetic trees using partial sequences derived from cox1 subunit gene, thus confirming the identity of the presence of cattle genotype (G5). This finding provides an alarming evidence for the circulation of the cattle genotype (G5) in the one-humped camels. It is, however, uncertain whether the G5-infected camel was maintained in Tamboul area, Central Sudan, or brought from a neighboring area for the purpose of slaughtering. The circulation of a major variant (G6) in Sudanese camels suggests that specific mechanisms are responsible for its persistence in the endemic area of Tamboul, Central Sudan. This is probably due to close relationships between dogs and camels in the study area . In rural communities with resource-poor settings, such as this study area, the practice of animal slaughtering is usually performed in the open space. Under these conditions, dogs would have free access to feed on livestock viscera, which may harbor hydatid cysts, the infective stage. Therefore, it is believed that this practice of livestock slaughtering could effectively contribute towards the persistence of the camel genotype (G6) in the study area.
The addition of EG strains sequences from Sudan enhances our understanding of the expansion and, to some degree, maintenance of the parasites in the intermediate hosts. Ongoing surveillance and EG strains characterization should also aid in determining the distribution of this cestode parasite in the country. As more sequencing data and prediction tools become more accurate and available, these data will provide the public health authorities an opportunity to anticipate and prepare for treatment and subsequent control programs for the disease.
In conclusion, the result of this study indicates the circulation of E. ortleppi, the cattle genotype (G5) in the one humped camel in the Sudan for the first time. Therefore, the G5 strain should be considered during epidemiological surveys for this important parasitic infection in Sudan. In addition, this investigation expands on the existing data on sequences generated from EG isolates recovered from the one humped camel in the endemic area of Tamboul, Central Sudan.
Thompson RC: The taxonomy, phylogeny and transmission of Echinococcus. Exp Parasitol. 2008, 119: 439-446. 10.1016/j.exppara.2008.04.016.
Wahlers K, Menezes CN, Wong ML, Zeyhle E, Ahmed ME, Ocaido M, Stijnis C, Romig T, Kern P, Grobusch MP: Cystic echinococcosis in sub-Saharan Africa. Lancet Infect Dis. 2012, 12: 871-880. 10.1016/S1473-3099(12)70155-X.
Abdel Malek E: Check-list of helminth parasites of domestic animals in the Sudan. Indian Vet J. 1959, 36: 28-286.
Eisa AM, Mustfa AA, Soliman KN: Preliminary report on cysticercosis and hydatidosis in the Southern Sudan. Sud J Vet Sc Anim Husb. 1962, 3: 97-102.
Eisa AM, El Khawat SE, Slepnev NK: A survey of parasites of dogs in Khartoum province. In Proceedings of the 8th Veterinary Conference.Khartoum- Sudan: Livestock and animal production in the Sudan; 1977.
El Khawat SE, Eisa AM, Slepnev NK, Saad MB: Hydatidosis of domestic animals in the central region of the Sudan. Bull Anim Hlth Prod Afr. 1979, 27: 249-251.
El Husein AM, Ali SM: A note on hydatidosis in camels at El Damer province, Northern State, Sudan. Sud J Vet Res. 1990, 10: 63-64.
Elmahdi IE, Ali QM, Magzoub MM, Ibrahim AM, Saad MB, Romig T: Cystic echinococcosis of livestock and humans in central Sudan. Ann Trop Med Parasitol. 2004, 98: 473-479. 10.1179/000349804225003578.
Omer RA, Dinkel A, Romig T, Mackenstedt U, Elamin M, Elnahas A, Aradaib I: Strains characterization of human hydatidosis in Sudan. Int Arch Hydatidosis. 2004, 35: 41.
Osman AM: Pathological and Biomolecular Study on Hydatid Disease in Camels. Cattle and Sheep. M.V.Sc: Thesis, University of Khartoum; 2007.
Ahmed ME: Surgical pulmonary hydatid disease in Khartoum, Sudan. Athens,Greece: XXII International Congress of Hydatidology; 2007.
Ahmed ME, Abdelrahim MI, Ahmed FM: Hydatid disease, a morbid drop needs awareness. Sudan Med J. 2011, 47: 4-8.
Musa NO, Eltom K, Awad S, Gameel AA: Causes of condemnation of sheep carcasses in abattoirs in Khartoum. In Tropentag Book of Abstracts:Ahmed et al. BMC Veterinary Research 2013, 9:127 Page 5 of 6 http://www.biomedcentral.com/1746-6148/9/127 Tropentag 19–21 Sept. 2012. Edited by Tielkes E. Göttingen, Germany:Cuvillier; 2012:54.
Bowles J, Blair D, McManus DP: Genetic variants within the genus Echinococcus identified by mitochondrial DNA sequencing. Mol Biochem Parasitol. 1992, 54: 165-173. 10.1016/0166-6851(92)90109-W.
Bowles J, Blair D, McManus DP: A molecular phylogeny of the genus Echinococcus. Parasitology. 1995, 110: 317-328. 10.1017/S0031182000080902.
Bowles J, McManus DP: Rapid discrimination of Echinococcus species and strains using a polymerase chain reaction-based RFLP method. Mol Biochem Parasitol. 1993, 57: 231-239. 10.1016/0166-6851(93)90199-8.
Bowles J, McManus DP: NADH dehydrogenase 1 gene sequences compared for species and strains of the genus Echinococcus. Int J Parasitol. 1993, 23: 969-972. 10.1016/0020-7519(93)90065-7.
Bowles J, McManus DP: Molecular variation in Echinococcus. Acta Trop. 1993, 53: 291-305. 10.1016/0001-706X(93)90035-A.
Scott JC, Stefaniak J, Pawlowski ZS, McManus DP: Molecular genetic analysis of human cystic hydatid cases from Poland: identification of a new genotypic group (G9) of Echinococcus granulosus. Parasitology. 1997, 114: 37-43. 10.1017/S0031182096008062.
Zhang LH, Chai JJ, Jiao W, Osman Y, McManus DP: Mitochondrial genomic markers confirm the presence of the camel strain (G6 genotype) of Echinococcus granulosus in north-western China. Parasitology. 1998, 116: 29-33. 10.1017/S0031182097001881.
McManus DP: The molecular epidemiology of Echinococcus granulosus and cystic hydatid disease. Trans R Soc Trop Med Hyg. 2002, 96: 51-157.
Lavikainen A, Lehtinen MJ: Meri, T, Hirvila-Koski V, Meri S: Molecular genetic characterization of the Fennoscandien cervid strain, a new genotypic group (G10) of Echinococcus granulosus. Parasitology. 2003, 127: 207-215. 10.1017/S0031182003003780.
Eckert J, Deplazes P: Biological, epidemiological, and clinical aspects of echinococcosis, a zoonosis of increasing concern. Am Soc Mic. 2004, 17: 107-135.
Busi M, Snabel V, Vercasia A, Garippa G, Perrone V, De Liberato C, D’Amelio S: Genetic variation within and between G1 and G3 genotypes of Echinococcus granulosus in Italy revealed by multilocus DNA sequencing. Vet Parasitol. 2007, 150: 75-83. 10.1016/j.vetpar.2007.09.003.
Moks E, Jõgisalu I, Valdmann H, Saarma U: First report of Echinococcus granulosus G8 in Eurasia and a reappraisal of the phylogenetic relationships of genotypes G5-G10. Parasitology. 2008, 135: 647-654.
Ahmed ME, Aradaib IE: Molecular characterization of Echinococcus granulosus using polymerase chain reaction. Int J Trop Med. 2006, 1: 30-32.
Omer RA, Dinkle A, Romig T, Mackenstedt U, Elnahas AA, Aradaib IE, Ahmed ME, Elmalik KH, Adam A: A molecular survey of cystic echinococcosis in Sudan. Vet Parasitol. 2010, 169: 340-346. 10.1016/j.vetpar.2010.01.004.
Ibrahim K, Romig T, Peter K, Omer RA: A molecular survey on cystic echinococcosis in Sinnar area, Blue Nile state (Sudan). Chinese Med J. 2011, 124: 2829-2833.
Thompson RC, McManus DP: Towards a taxonomic revision of the genus Echinococcus. Trends Parasitol. 2002, 18: 452-457. 10.1016/S1471-4922(02)02358-9.
McManus DP, Thompson RCA: Molecular epidemiology of cystic echinococcosis. Parasitology. 2003, 127: 37-51. 10.1017/S0031182003003275.
Nakao M, McManus DP, Schantz PM, Craig PS, Ito A: A molecular phylogeny of the genus Echinococcus inferred from complete mitochondrial genomes. Parasitology. 2007, 134: 713-722. 10.1017/S0031182006001934.
Nakao M, Li T, Han X, Ma X, Xiao N, Qiu J, Wang H, Yanagida T, Mamuti W, Wen H, Moro PL, Giraudoux P, Craig PS, Ito A: Genetic polymorphisms of Echinococcus tapeworms in China as determined by mitochondrial and nuclear DNA sequences. Int J Parasitol. 2010, 40: 7.
Thompson RC, Lymbery AJ: The nature, extent and significance of variation within the genus Echinococcus. Adv Parasitol. 1988, 27: 209-258.
Thompson RC, Lymbery AJ, Constantine CC: Variation in Echinococcus: towards a taxonomic revision of the genus. Adv Parasitol. 1995, 35: 145-176.
McManus DP, Bowles J: Molecular genetic approaches to parasite identification: their value in diagnostic parasitology and systematics. Int J Parasitol. 1996, 26: 687-704. 10.1016/0020-7519(96)82612-9.
Tamura K, Peterson D, Peterson N, Stecher G, Nei M, Kumar S: MEGA5: Molecular Evolutionary Genetics Analysis using Maximum Likelihood, Evolutionary Distance, and Maximum Parsimony Methods. Mol Biol Evol. 2011, 28: 2731-2739. 10.1093/molbev/msr121.
Saad MB, Magzoub M: Echinococcus granulosus infection in Tamboul, Sudan. J Helminthol. 1986, 60: 299-300. 10.1017/S0022149X0000852X.
Saad MB, Magzoub M: Hydatidosis in camels and cattle in the Sudan. Sud J Vet Sci Anim Husb. 1989, 28: 27-32.
Dinkel A, Njoroge E, Zimmerman A, Wälz M, Zeyhle E, Elmahdi I, Mackenstedt U, Romig T: A PCR system for detection of species and genotypes of the Echinococcus granulosus complex, with reference to the epidemiological situation in eastern Africa. Int J Parasitol. 2004, 34: 645-653. 10.1016/j.ijpara.2003.12.013.
This study was made possible by the invaluable assistance provided by Dr. Ibrahim Higazi, Veterinary Medical officer, Tamboul slaughterhouse, Tamboul, Gezira State, Central Sudan. The authors are very grateful to Mr. Abdall M. Fadl Elmoula for technical assistance. This study received financial support from the Deanship of Scientific Research Alneelain University, Sudan.
The authors declare that they have no competing interests.
MEA collected hydatid cyst samples, extracted the DNA and optimized the polymerase chain reaction-based detection assay; KAE analyzed the sequences and designed the study; IAA and FME designed the experiment; NOM, edited the sequences and helped with experimental design. IEA designed the experiment and prepared the final manuscript. All authors read and approved the final version of the manuscript.