Molecular classification of echinococcus granulosus strains from livestock animals in Libya.

Link t full russian text
Download PDF format

Russian Journal of Parasitology, 2016, V.38, Iss.4
DOI: …………………
Received 01.09.2016
Accepted 28.11.2016


MOLECULAR CLASSIFICATION OF ECHINOCOCCUS GRANULOSUS STRAINS FROM LIVESTOCK ANIMALS IN LIBYA

Elmajdoub O. L1, W. A. Rahman2, M. F.F. Wajdii2, M.N. Siti Azizah2
1 Zoology Department, Science College, Misurata University, Libya
elmajdoublayla@sci.misuratau.edu.ly
2School of Biological Sciences, Universiti Sains Malaysia, 11800 Minden, Penang

Abstract

This study demonstrates the distribution of various E. granulosus strains in sheep, cattle and camel. The common sheep strain G1 is mainly found in sheep and cattle, but also parasitized camels.  In contrast, the camel strain G6 is found mainly in camels and rarely in sheep and cattle. However, the study also revealed the possible presence of cryptic species that are closely related to both genotypes in livestock of Libya as evident by high mutations in several specimens. Based on the occurrence of overlapping hosts of E. granulosus in Libya, more research on the transmission cycles and genotypes of E. granulosus in Libya is required. In addition, it is suggested that surveys on potential intermediate hosts, including in humans with dogs as the major final host in Libya using higher resolution molecular tools such as microsatellite markers is recommended. 

Keywords: Echinococcus granulosus, liver, lung, common sheep strain G1

Introduction

Through the past five decades, significant phenotypic and genetic variabilities have been recognized and identified in various strains of E. granulosus isolated from different regions (Van Herwerden et al., 2000; Thompson and McManus 2001, 2002; Pearson et al., 2002 and Huttner et al., 2008) . These studies have revealed that the different strains of E. granulosus consist of heterogeneous groups of genetic variants (McManus, 2002). Thompson (1995) illustrated that different strains may display variations in morphology, host specificity, development rate, pathogenicity and geographical distributions. Moreover, many studies have been conducted to determine the host and geographic ranges of these strains, and whether genetic variations were characteristic and specific to the different endemic areas throughout the world (Jenkins and Thompson, 2005). Thompson and Kumaralilake (1982) observed that some strains of E. granulosus share similar morphological characters but showed epidemiological differences; thus, this parasite showed high diversity. Identification of strain types of E. granulosus is very important in strategizing and implementing an Echinococcosis control and management programmer. 

Until now, 10 strains or genotypes, namely G1- G10 have been recognized and described in E. granulosus based on mitochondrial and nuclear gene analyses (Bowles et al., 1992, 1994; Scott et al., 1997; Lavikainen et al., 2003; Thompson et al., 2008; Saarma et al., 2009). According to mitochondrial data, E. granulosus has been traditionally assigned to the various taxonomic species in relation to the G1- G10 genotypes; E. granulosus sensu stricto (G1, G2, and G3), E. equinus (G4), E. ortleppi (G5), and E. canadensis (G6-G10) (Nakao et al., 2007; Moks et al., 2008). However, recent studies based on nuclear data categorized the genotypes (G6- G10) into two strains; cervid genotypes G8 and G10 belonging to E. canadensis, whilst camel and pig genotypes G6, G7 belonged to E. intermedius (Saarma et al., 2009).

The aim of this study was to determine E. granulosus genotypes present in Libyan livestock by molecular genetic strain typing. For the analysis, the mitochondrial (ATP6) and nuclear (Act II) genes were utilized.

Materials and Methods

Parasite samples

Hydatid cyst samples were collected from the major slaughtered livestock namely sheep, camel and cattle from government abattoirs located in four regions in Libya from January to end of 2010. In total, 120 samples were used for molecular analysis, from liver, lung and other organs from the three major livestock. To obtain fertile cysts, the protoscoleces were collected from the sediment of hydatid sand by pipette under sterile conditions and rinsed three times with normal saline. Suspensions of protoscoleces were fixed in 90% (v/v) ethanol and then stored at - 20°C until DNA extraction. To obtain infertile cysts, the thin germinal layer from the wall cyst was cut into small pieces and washed with normal saline and then stored at - 20°C for further processing.

DNA extraction

The protoscoleces were washed several times in nucleic acid-free water to remove the alcohol preservation solution. Genomic DNA was extracted using a QIAamp DNeasy mini kit (Qiagen, Germany) according to the manufacturer`s protocol.

Electrophoresis was preceded on a 1% agarose gel with GelRed Nuclei Acid Gel Stain as marker, at 100 volt for 30 min to assess the success of DNA extraction. After that, the agarose gel was visualized in a gel documentation system (GENE Flash, Syngene Bio Imaging, USA) for the presence of the extracted DNA bands. To investigate the presence of sterile cyst (from cattle), DNA extraction was conducted on 25 mg of the cut pieces of infertile germinal layer and placed into a 1.5 ml Eppendrof tube. Then according to the manufacturer`s protocol.

Polymerase chain reaction (PCR) amplification

The partial fragments of the mitochondrial gene (ATP6) were PCR-amplified on the DNA extract using specific primers designed by Xiao et al. (2005).

ATP6 forward: 5`- GCA TCA ATT TGA AGA GTT GGG GAT AAC-3`

ATP6 reverse: 5`- CCA AAT AAT CTA TCA ACT ACA CAA CAC-3`

The PCR reaction contained 5.5 µL of 5X PCR buffer, 4 µL 25 mM MgCl2 solution, 0.3 µL of 5u/ µL Taq DNA polymerase, 0.7 µL of 10 mM dNTP (Promega, USA), 0.5 µL of each primer and 2.0 µL of the target DNA in a total volume of 25 µL. The PCR protocol consisted of an initial incubation at 94°C for 30 s, 35 cycles at 94 °C, 55 °C for 30 s and 72 °C for 1 min, 72 °C for 5 min in the final extension, The PCR amplification was conducted using the Mastercycler Gradient- Thermal cycler (Eppendorf, Germany).

The Act II fragment was PCR amplified using primers specifically designed by De Silva et al. (1993).

Act II forward: 5`- TCT TCC CCT CTA TCG TGG G-3`

Act II reverse: 5`- CTA ATG AAA TTA GTG CTT TGT GCG C-3`

The PCR was carried out in a 25 µL volume containing 2 µL target DNA, 5 µL of 5X PCR buffer, 5 µL 25 mM MgCl2 solution, 0.25 µL of 5u/ µL Taq DNA polymerase, 0.5 µL of 10 mM dNTP (Promega, USA), 0.5 µL of each primer. The PCR conditions were as follows: 94°C for 30 s, 40 cycles at 94 °C for 30 s, 60 °C for 30 s and 72 °C for 1.5 min, and then a final extension at 72 °C for 5 min, The PCR amplification used the Mastercycler Gradient- Thermal cycler (Eppendrof, Germany). 

After the PCR amplification, PCR products were assessed by electrophoresis in a 1.5 % agarose gel using GelRed Nuclei Acid Gel Stain as marker at 100 volt for 25 min. finally the purified DNA samples were sent to the service provider for sequencing procedure (First BASE Laboratories Sdn. Bhd. and Centre of Chemical Biology, CCB at Universiti Sains Malaysia).

Data analysis

The obtained sequences were edited using MEGA 5.05 program (Tamura et al., 2007). The genetic relationships between haplotypes were determined by constructing phylogeny trees based on Neighbour-Joining (NJ). GenBank sequences of E. granulosus were downloaded for comparisons with the current data for each gene. For ATP6 gene the following GenBank sequences were used- Acc. No AF 297617.1 sheep strain and AB208063.1 camel strain G6. For Act II gene the GenBank sequences AF 528499.1 sheep strain G1 and AF 528500.1 camel strain G6 were used to compare with the sequences in this study.        

Results

A total of 120 E. granulosus samples were successfully amplified by using the optimized PCR conditions. Length of partial fragment of ATP6 mtDNA was 513bp and partial fragment of Act II DNA fragment was 267bp in length.

Multiple sequence alignments of ATP6 and Act II genes were carried out for 102 individuals from livestock hydatid cyst (60 sheep, 31 camels and 11 cattle) and 89 individuals from livestock hydatid cyst (54 sheep, 25 camels and 10 cattle) respectively using MEGA 5 Software ClustalW (1.6) DNA weight matrix. Blastn search showed that all the aligned sequences belonged to ATP6 and Act II genes with high similarities (96-100%).

Neighbour Joining (NJ) analysis

The Neighbour-Joining (NJ) analysis of ATP6 gene was carried out based on Kimura 2- parameter as presented in Fig 1. The same GenBank sequences of E. granulosus from the previous sheep population study (AB 208063.1; AB235847.1; AB031283.1) were included. Two main clusters were formed; the first cluster consisted of mixed populations (different organs from different animals) with a low bootstrap confidence level of 51% and a second cluster with high support (99%). This cluster was divided into seven internal subclusters with low to moderate support but with no obvious genetic relationships to host or organ types. The GenBank (AB031283.1) representing sheep strain grouped in this cluster. The second monophyletic cluster with 99% support consisted of four subclusters.

Unlike Cluster 1, the second cluster was generally made up of parasites of camel origin with GenBank taxon (AB 208063.1) and (AB235847.1) representing camel strain grouping in this cluster. However, there were two exceptions, namely of parasites from host sheep, SMT1 (KF255896) and SMT2 (KF255899) grouping into this cluster. This provided further evidence of the previous population analysis that these two sheep parasite sequences in the mesentery are more closely related to camel strain. The four subclusters were low to moderately support. These internal subclusters did not show any obvious pattern to specific organs.

The NJ tree formed with 1000 replicates of Act II gene is presented in Fig 2. Sheep strain (gi 22653316) and camel strain (gi 22653318) GenBank haplotypes. There was no phylogenetic structuring into major clusters. However, a strongly supported cluster (99%) was observed consisting of several haplotypes (KP843657 with groups of KP843653, KP843652, KP843647, KP843639, KP843626) from mixed populations and the GenBank G6 camel strain haplotype. Four weakly supported and a single moderately supported terminal clusters consisting of mixed populations (different organs from different regions) were also observed.



Fig 1: Neighbour joining phylogenetic tree of E. granulosus haplotypes from Libyan combined sheep, cattle and camel livestock of ATP6 geneFig 1: Neighbour joining phylogenetic tree of E. granulosus haplotypes from Libyan combined sheep, cattle and camel livestock of ATP6 gene.


Fig 2: Neighbour joining phylogenetic tree of E. granulosus haplotypes from Libyan livestock of Act II gene.Fig 2: Neighbour joining phylogenetic tree of E. granulosus haplotypes from Libyan livestock of Act II gene.

Discussion

To date, six genotypes of E. granulosus complex have been identified in Africa; G1 (sheep strain), G2 (Tasmanian sheep strain), G4 (horse strain), G5 genotype (cattle strain) and the G6 (camel strain) genotypes (Thompson and McManus, 2002; Dinkel et al., 2004; Huttner et al., 2008 and Casulli et al., 2010). Limited studies in Libya, primarily conducted by Tashani et al. (2002) in eastern Libya reported that all livestock animals (sheep, cattle and camels) were infected with the same genotype of E. granulosus (G1, sheep strain). Another study by Abushhewa et al. (2010) from different areas in Libya recorded two groups; the first group belonging to G1-G3, and the second group belonging to G6-G10. Thus, the present study is a valuable contribution of data on the population genetics, diversity as well as complements existing taxonomic knowledge of this parasitic species from various organs in different hosts from Libya based on each of a mitochondrial and nuclear gene.  

             The DNA sequence variation of the partial mtDNA ATP6 of E. granulosus conducted in this investigation has indicated the transmission of two main strains in livestock animals (sheep, camels and cattle) in Libya. This study is the first documentation on ATP6 sequences of E. granulosus from livestock in Libya and has provided data on the common strain distributed in Libya of E. granulosus.

Various mutations in different populations of sheep, camels and cattle were observed. Of these, 80.5% and 19.4% haplotypes belonged to common sheep strain (G1) and camel strain (G6) respectively. These findings corresponded to the results by Eryildiz and Sakru (2012) in a Turkish study, who recorded that most of their isolates, belonged to G1-G3. Only one isolate belonged to G6-G10. Moreover, Hailemariam et al. (2012) recorded similar findings in an Ethiopian study which observed 87% of livestock samples were identified as G1 and 13% as the G6. On the other hand, Abushhewa et al. (2010) in a Libyan study reported that all isolates from camels belonged to camel strain G6. This indicated the all camels in their study lived in the private pastures with no connection to sheep farms.

Phylogenetic analysis of E. granulosus populations in livestock

Phylogenetic analysis based on NJ method generated trees with similar topologies, which appeared less efficient based on the clustering and bootstrap approaches. Generally, in the present study, the topology of the NJ tree from combined sheep, camel and cattle populations, showed that all sheep populations referred to as G1 genotype sheep strain combined together, while only two individuals belonged to G6 genotype camel strain. However, all cattle individuals belonged to sheep strain, while several individuals from camels belonged to the sheep strain while the rest belonged to camel strain. Previous studies conducted in Libya by Tashani et al. (2002) recorded that all E. granulosus isolates from different livestock (sheep, cattle and camels) belonged to the common sheep strain (G1). In addition, the present study showed that all hydatid cysts from cattle were sterile; due to the infection by the common sheep strain (G1). This was also observed by Tashani et al. (2002). 

Abushhewa et al. (2010) in their study on Libya and Omer et al. (2011) in a study in Sudan reported that all E. granulosus isolates from camels (100%) belonged to the G6-G10 complex. Interestingly, 87% of cattle hydatid cysts investigated by Abushhewa et al. (2010) and 99% by Omer et al. (2011) belonged to the G6 genotype. More recent evidence by Abdel Aaty et al. (2012) and Omer et al. (2011) reported that all isolates from sheep, camels, pigs and cattle were identified as G6 camel strain in Egypt and Sudan. Thus their findings suggested that the camel strain play the major role in the transmission cycle of E. granulosus in Egypt and Sudan.

Considering the overall livestock, high substitution rates were observed in haplotypes closely related to G6.  Furthermore, a single sample (camel liver) from this group had unusually high nucleotide substitutions. This sample was closely related to GenBank G6 Acc No AB208063.1 (97%) and GenBank G7 Acc no. AB235847.1 G7 (96%). Farjallah et al. (2007) found G7 in slaughtered camels from Tunisia and Mauritania. This is the first record of a strain belonging to G7 or very closely related to it. Therefore, it is suggested that the prevalence reported in this study is due to the movement of livestock animals from neighboring countries such as Tunisia or due to the occurrence of high random mutation in this sample.

Furthermore, many related studies in North Africa correspond to the present findings; from Algeria, Bart et al. (2004) observed two distinct well supported clusters (G1 and G6) based on the mitochondrial (ND1, COX1) genes. The same situation was observed in Middle Africa, from Kenya and Sudan, where many genetic studies (Bowles et al., 1992; Wachira et al., 1993 and Dinkel et al., 2004) have demonstrated the importance of G1 and G6 in livestock. But Dinkel et al. (2004) also noted other strains in Kenya originating from pig and in Sudan from cattle (E. ortleppi).

Nuclear gene Act II

To further investigate the strain identity of E. granulosus from sheep, camels and cattle in Libya, a 262 bp of Act II gene sequence data was analyzed. This is attributed to the lower mutation rates of nuclear markers compared to mitochondrial DNA. According to the phylogenetic tree two main genotypes, G1 (common sheep strain) and G6 (camel strain) with 98-100% homology with GenBank (AF528499 and AF528500 respectively) were observed. These findings were in agreement with previous studies by Gudewar et al. (2009) in India and Maillard et al. (2007) in Africa using the same Act II gene. Bart et al. (2006) used BG1/3 nuclear gene to identify E. granulosus among sheep, cattle and pigs, and identified two genotypes, a sheep strain and the pig strain.  

The present results showed that most of individuals from different sheep populations were placed in common sheep strain (G1), while only two samples (Misurata liver and mesentery) populations were placed in camel strain (G6). The results from mitochondrial and nuclear markers revealed that most of the sheep hydatid cysts from different organs in different areas belonged to the sheep strain, but were rarely infected by camel strain. This indicates that generally the camel strain was not effective in infecting the sheep host in Libya with several exceptions – in four sheep individuals, specifically 3 individuals from mesentery organ and another from liver organ infected with the camel strain. The G1 genotype is effective in infecting different organs in camels and cattle. In contrast, the camel strain is ineffective to infect the cattle host i.e. all cattle individuals were only infected with sheep strain.  For the nuclear marker only two samples from liver and lung cattle were infected by camel strain G6, while the sheep strain G1 was more suitable to infect cattle host than the G6. However, all hydatid cysts in different organs of cattle were sterile. This indicates that sheep and camel strains are not the adapted for propagation in cattle in Libya.  Furthermore, the G6 strain infection was also recorded in several individuals of sheep and cattle based on the nuclear gene, an observation also previously noted by Kamenetzky et al. (2002) in Argentina and Haag et al. (2004).  Thus, both the mitochondrial and nuclear gene analyses generated two main clades (G1 and G6) of E. granulosus which represented the common intermediate hosts (sheep, cattle and camels) from Libya as identified based to the Genbank databases.

Why are there incidences of G1 and G6 being present in the same host? An explanation that could be put forward is the interaction between the camel-dog and sheep-dog cycles. Dog is the final host of this parasite. Due to the co-existence of camels with sheep and cattle together in a close neighborhood, cross transmission of camel and sheep strains may occur in different livestock as a result of overlapping cycles. Clearly, more focused studies on the G1 and G6 in overlapping hosts of E. granulosus in Libya need to be conducted. In addition, there is a possibility of cross fertilization between genotypes of G1 and G6 in the dog which is the final host and thus, producing a hybrid genotype.

References

1.      Abushhewa, M. H., M. H. S. Abushhiwa, M. J. Nolan, A. R. Jex, B. E. Campbell, A. Jabbar and R. B. Gasser. 2011. Genetic classification of Echinococcus granulosus cysts from humans, cattle and camels in Libya using mutation scanning-based analysis of mitochondrial loci. Molecular and Cellular Probes 24(6):346-351.

2.      Bart, J. M., K. Bardonnet, M. C. B. Elfegoun, H. Dumon, L. Dia, D. A. Vuitton and R. Piarroux. 2004. Echinococcus granulosus strain typing in North Africa: Comparison of eight nuclear and mitochondrial DNA fragments. Parasitology 128(2):229-234.

3.      Bowles, J., D. Blair and D. P. McManus. 1992. Genetic variants within the genus Echinococcus identified by mitochondrial DNA sequencing. Molecular and Biochemical Parasitology 54(2):165-173.

4.      Bowles, J., D. Blair and D. P. McManus. 1994. Molecular genetic characterization of the cervid strain ('northern form') of Echinococcus granulosus. Parasitology 109(2):215-221.

5.      Da Silva, C. M. D., H. B. Ferreira, M. Picon, N. Gorfinkiel, R. Ehrlich and A. Zaha. 1993. Molecular cloning and characterization of actin genes from Echinococcus granulosus. Molecular and Biochemical Parasitology 60(2):209-220.

6.      Farjallah, S., M. Busi, M. O. Mahjoub, B. B. Slimane, K. Said and S. D'Amelio. 2007. .Molecular characterization of Echinococcus granulosus in Tunisia and Mauritania by mitochondrial rrnS gene sequencing. Parassitologia 49(4):239-246.

7.      Gudewar, J., D. Pan, A. K. Bera, S. K. Das, A. Konar, J. R. Rao, A. K. Tiwari and D. Bhattacharya. 2009. Molecular characterization of Echinococcus granulosus of Indian animal isolates on the basis of nuclear and mitochondrial genotype. Molecular Biology Reports 36(6):1381-1385.

8.      Haag, K. L., F. J. Ayala, L. Kamenetzky, A. M. Gutierrez and M. Rosenzvit. 2004. Livestock trade history, geography, and parasite strains: The mitochondrial genetic structure of Echinococcus granulosus in Argentina. Journal of Parasitology 90(2):234-239.

9.      Hailemariam, Z., M. Nakao, S. Menkir, A. Lavikainen, T. Yanagida, M. Okamoto and A. Ito. 2012 Molecular identification of unilocular hydatid cysts from domestic ungulates in Ethiopia: Implications for human infections. Parasitology International 61(2):375-377.

10. Huttner, M., M. Nakao, T. Wassermann, L. Siefert, J. D. F. Boomker, A. Dinkel, Y. Sako, U. Mackenstedt, T. Romig and A. Ito. 2008. Genetic characterization and phylogenetic position of Echinococcus felidis Ortlepp, 1937 (Cestoda: Taeniidae) from the African lion International Journal for Parasitology 38(7):861-868.

11. Jenkins, D. J., T. Romig and R. C. A. Thompson. 2005. Emergence/re-emergence of Echinococcus spp. - A global update. International Journal for Parasitology 35(11-12):1205-1219.

12. Kamenetzky, L., A. M. Gutierrez, S. G. Canova, K. L. Haag, E. A. Guarnera, A. Parra, G. E. García and M. C. Rosenzvit. 2002. Several strains of Echinococcus granulosus infect livestock and humans in Argentina. Infection, Genetics and Evolution 2(2):129-136.

13. Lavikainen, A., M. J. Lehtinen, T. Meri, V. HirvelÃ-Koski and S. Meri. 2003. Molecular genetic characterization of the Fennoscandian cervid strain, a new genotypic group (G10) of Echinococcus granulosus. Parasitology 127(3):207-215.

14. Maillard, S., M. C. Benchikh-Elfegoun, J. Knapp, J. M. Bart, P. Koskei, B. Gottstein and R. Piarroux. 2007. Taxonomic position and geographical distribution of the common sheep G1 and camel G6 strains of Echinococcus granulosus in three African countries. Parasitology Research 100(3):495-503.

15. McManus, D. P. 2002. The molecular epidemiology of Echinococcus granulosus and cystic hydatid disease. Transactions of the Royal Society of Tropical Medicine and Hygiene 96(SUPPL. 1):S1/151-S151/157.

16. Omer, R. A., A. Dinkel, T. Romig, U. Mackenstedt, A. A. Elnahas, I. E. Aradaib, M. E. Ahmed, K. H. Elmalik and A. Adam. 2011. A molecular survey of cystic echinococcosis in Sudan. Veterinary Parasitology 169(3-4):340-346.

17. Pearson, M., T.H.Le., L.H. Zhang, D. Blair, T.H.N. Dai and D.P. MeManus 2002. Molecular taxonomy and strain analysis in Echinococcus. In: craig, P. and Z. Pawlowski (Eds). Cestode Echinococcosis and cysticercosis an emergent and global problem IOs Press, Amsterdam, the Netherlands, pp: 205-219.

18. Scott, J. C., J. Stefaniak, Z. S. Pawlowski and D. P. McManus. 1997. Molecular genetic analysis of human cystic hydatid cases from Poland: Identification of a new genotypic group (G9) of Echinococcus granulosus. Parasitology 114(1):37-43.

19. Tamura, K., J. Dudley, M. Nei and S. Kumar. 2007. MEGA4: Molecular Evolutionary Genetics Analysis (MEGA) software version 4.0. Molecular Biology and Evolution 24(8):1596-1599.

20. Tashani, O. A., L. H. Zhang, B. Boufana, A. Jegi and D. P. McManus. 2002. Epidemiology and strain characteristics of Echinococcus granulosus in the Benghazi area of eastern Libya. Annals of Tropical Medicine and Parasitology 96(4):369-381.

21. Thompson, R. C. A and D. P. McManus 2001. Aetiology: parasites and life cycles. In: Eckert, J., Gemmell, M.A., Meslin, F.X. and Pawlowski, Z.S. (eds) WHO/ OIE Manual on Echinococcosis in humans and animals a public health problem of global concern. World Organisation for Animal Health, Paris, pp: 1-19.

22. Thompson, R. C. A. 1995. Biology and systematic of Echinococcus and Hydatid disease. Thompson R.C.A and Laymbery, A.J. (eds). CAB International Wallingford, Oxen, 1-37.

23. Thompson, R. C. A. 2008. The taxonomy, phylogeny and transmission of Echinococcus. Experimental Parasitology 119(4):439-446.

24. Thompson, R. C. A. and D. P. McManus. 2002. Towards a taxonomic revision of the genus Echinococcus. Trends in Parasitology 18(10):452-457.

25. Thompson, R. C. A. and L. M. Kumaratilake. 1982. Intraspecific variation in Echinococcus granulosus: the Australian situation and perspectives for the future. Transactions of the Royal Society of Tropical Medicine and Hygiene 76(1):13-16.

26. Van Herwerden, L., R. B. Gasser and D. Blair. 2000. ITS-1 ribosomal DNA sequence variants are maintained in different species and strains of Echinococcus.  International Journal for Parasitology 30(2):157-169.

27. Van Herwerden, L., R. B. Gasser and D. Blair. 2000. ITS-1 ribosomal DNA sequence variants are maintained in different species and strains of Echinococcus.  International Journal for Parasitology 30(2):157-169.

28. Wachira, T. M., J. Bowles, E. Zeyhle and D. P. McManus. 1993. Molecular examination of the sympatry and distribution of sheep and camel strains of Echinococcus granulosus in Kenya. American Journal of Tropical Medicine and Hygiene 48(4):473-479.

29. Xiao, N., J. Qiu, M. Nakao, K. Nakaya, H. Yamasaki, Y. Sako, W. Mamuti, P. M. Schantz, P. S. Craig and A. Ito. 2003. Short report: Identification of Echinococcus species from a yak in the Qinghai-Tibet plateau region of China. American Journal of Tropical Medicine and Hygiene 69(4):445-446.

© 2016 The Author(s). Published by All-Russian Scientific Research Institute of Fundamental and Applied Parasitology of Animals and Plants named after K.I. Skryabin. This is an open access article under the Agreement of 02.07.2014 (Russian Science Citation Index (RSCI)http://elibrary.ru/projects/citation/cit_index.asp) and the Agreement of 12.06.2014 (CA-BI.org/Human Sciences section: http://www.cabi.org/Uploads/CABI/publishing/fulltext-products/cabi-fulltext-material-from-journals-b...)