https://doi.org/10.5812/archcid.14168

Discrimination of Mixed Infections of Echinococcus Species Based on in Silico Sequence Analysis: A New Way of Reflecting Overlapped Strains in Indigenous Areas

authors:

avatar Maryam Karimi 1 , avatar Reza Ghasemikhah 2 , avatar Hadi Mirahmadi ORCID 3 , avatar Adel Spotin ORCID 4 , avatar Soheila Rouhani 1 , avatar Seyyed Javad Seyyed Tabaei 1 , *

Department of Parasitology and Mycology, Shahid Beheshti University of Medical Sciences, Tehran, Iran
Department of Parasitology and Mycology, Arak University of Medical Sciences, Arak, Iran
Department of Parasitology and Mycology, Zahedan University of Medical Sciences, Zahedan, Iran
Immunology Research Center, Tabriz University of Medical Sciences, Tabriz, Iran
Archives of Clinical Infectious Diseases: Vol.12, issue 4; e14168
published online: September 30, 2017
article type: Research Article
received: August 17, 2016
revised: May 28, 2017
accepted: September 19, 2017

how to cite: Karimi M, Ghasemikhah R, Mirahmadi H, Spotin A, Rouhani S, et al. Discrimination of Mixed Infections of Echinococcus Species Based on in Silico Sequence Analysis: A New Way of Reflecting Overlapped Strains in Indigenous Areas. Arch Clin Infect Dis. 2017;12(4):e14168. https://doi.org/10.5812/archcid.14168.

Abstract

Objectives:

Hydatidosis, as a cyclo-zoonotic indigenous disease, is marked by various strains and multiple infections around the world where different overlapping cycles of Echinococcus spp. are explicitly circulating among intermediate and definitive hosts. One of the current challenges in Echinococcus taxonomy is associated with the imprecise identification of dual infections in endemic regions. The aim of this study was to genotype Echinococcus strains, particularly those ignored during mixed infections, by employing phylogenetic analysis in order to determine the accurate status of the parasite in endemic areas.

Methods:

60 isolates were collected from the liver and lungs of twenty humans, fifteen sheep, fifteen goat, and ten cattle during 2012 - 14 in Markazi province, Iran. Subsequently, the extraction of DNA was performed, followed by amplification and identification of the DNA through both PCR for internal transcribed spacer 1 (ITS1-rDNA) of ribosomal DNA and restriction fragment length polymorphism (RFLP) using two endonuclease ‎enzymes: Rsa1 and HpaII. Then, the isolates were directly sequenced in order to accomplish the definite identification of strains/haplotypes. In addition, the mixed infections were analyzed by mixed sequences reader (MSR).

Results:

The RFLP, phylogenetic analyses of ITS1-rDNA sequences and MSA finding based on analyzed chromatograms strongly revealed the sheep strain (n = 58) and two mixed infection G1 and G6 genotypes (n = 2). The five specific haplotypes that represent moderate diversity indices were identified. The haplotype diversity was estimated to be 0.522, with a nucleotide diversity of 0.02.

Conclusions:

The G1 was reported to be the dominant strain in the central regions of Iran; however, the first dual discrimination of G1 and G6 genotypes in these regions showed that a sheep-camel/dog life cycle is circulating in central Iran. The current findings can facilitate the parasite taxonomic classification based on in silico analysis particularly about those that are neglected in dual infections. Additional studies are highly needed with a greater sample size and in various regions of Iran.

1. Background

Cystic Echinococcosis (CE), a life-threatening zoonotic infection, is caused by the metacestode, a larval life stage of a dog tapeworm, of Echinococcus granulosus. Currently, ten genotypes (G1-G10) of genus Echinococcus sensu lato have been characterized by using mitochondrial markers of 16 species and 13 subspecies (1-3). Interestingly, the following taxonomic positions of these genotypes are defined based on an extensive revision: E. granulosus sensu stricto (G1 - G3), E. equines (G4) that are classified as non-zoonotic strains, and E. ortleppi (G5), E. canadensis (G6-G10) and E. felidis (‘lion strain’) as zoonotic strains (4-7). However, some genotypes have remained unknown in the consequence of overlapping cycles among various intermediate - definitive hosts (8, 9).

Epidemiologically, CE is distributed in hyperendemic regions of the world, including the Middle East, Northern and Eastern Africa, Central Asia, Mediterranean, Eastern Europe, and some foci in the United Kingdom (3, 10) where different intermediate hosts are circulating sympatrically. Nevertheless, due to the cross transmissions of Echinococcus species among diverse hosts and host-parasite interactions, the occurrence of mixed (dual/triple) infections would not be unexpected. In biological sequence analysis, based on bioinformatics tools, the mixed infections are suspected due to multiple overlapping peaks (chromatograms) observed in some sequences (11, 12). The lack of precise diagnosis of mixed infections leads to an insufficient understanding of the pathogenesis, classification, and biology of the parasite in pertinent regions.

To date, characterization of the mixed (concomitant) infections has been practically assisted by cloning and/or multiplex Real-Time PCR strategies that employ specific probes-primers in some parasitic infections (12-14). Some of the challenges of implementing these methods are their high cost and the need for reference laboratories, as well as designing specific probe-primers and optimizing conditions. Notably, discrimination of these infections has not been supported by PCR-RFLP and Nested-PCR methods due to the formation of nonspecific digestions and unidentified bands (15-17). Therefore, utilizing an available in silico program to analyze DNA sequences with heterozygous bases can reflect the valid status of neglected Echinococcus spp. in endemic areas. Furthermore, the identification of distinct genotypes by phylogenetic analysis together with clarifying their potentially novel strains/haplotypes demonstrates the differences in morphology aspects, transmission dynamics, and development and pathogenicity rate (1, 3, 7).

So far, at least four genotypes (G1, G2, G3, and G6) of E. granulosus have been reported in different parts of Iran (9, 18-20) where according to the circulation of a variety of intermediate hosts, this seems to be underestimated. However, many Iranian researchers are focusing on genotyping of E. granulosus in different regions of Iran and among populations of sheep, cattle, goat, and camel (19-24). With regard to the various intermediate hosts, including sheep, buffalo, pig, camel, goat, and cervid in these regions, the isolation and characterization of mixed genotypes as well as discrimination of their heterogeneity traits warrant a more realistic image of the parasite status in the region.

The present study aimed to introduce an in silico program for analyzing DNA sequences with heterozygous bases via phylogenetic analyses. This could facilitate the identification of E. granulosus isolates and provide insights into the determination of neglected Echinococcus spp. taxonomic status in central Iran.

2. Methods

2.1. Sample Preparation

A total of sixty hydatid cysts were obtained from hospitals and the slaughterhouse of Arak (Markazi province, geographically located in central Iran). The cysts were taken separately from 20 humans, 15 sheep, 15 goats, and 10 cattle liver and lung tissues. Then, these samples were aseptically fixed using 80% (v/v) ethanol. In addition, the presence of possible protoscolices was examined by optical microscopes and on aspirated hydatid fluids. The viability rate of protoscolices was confirmed by 0.1% aqueous eosin staining, preceded by washing with phosphate buffered saline (PBS, pH 7.2) (25). Whilst the live protoscolices remained colorless, dead ones absorbed eosin and became reddish. For additional molecular modalities, storage of the washed protoscolices was performed at 4°C in 70% ethanol.

2.2. Processes of DNA Extraction and Polymerase Chain Reaction (PCR)

Protoscolices taken from each sample were washed three times with phosphate-buffered saline (pH 7.2). Subsequently, they were sonicated to facilitate the breakdown of protoscolices. The extraction of total genomic DNA from either 50 μL of each samples’ protoscolices or germinal layers was performed based on the standard extraction procedure recommended by QIA DNeasy blood and tissue kit (QIAGEN, Hilden Germany) according to the manufacturer’s instructions. The extracted DNA concentration was assessed by NanoDrop system (thermo scientific Inc., Wilmington, DE). Furthermore, the extracted DNA was stored at -20°C until the molecular assessments. The forward (BD1:5’-GTC GTA ACA AGG TTT CCG TA-3’) and reverse primers (‎4S: 5’-TCT AGA TGC GTT CGA A(G/A)T GTC GAT G-3‎’) were used to amplify the 1000 bp ribosomal DNA (rDNA), also known as internal transcribed spacer 1 (ITS1). The following thermal cycling profile was exerted ‎for each PCR reaction (total volume: 30 μL): 1 μL of genomic ‎DNA (100 ng), 15 μL‏ ‏Master mix (2 ×), ‎1 μL‎ (20 ‎pmol) of primers (forward and reverse), and ‎12 μL‎ Double distilled water (DDW). The primary denaturation cycle was conducted at 94°C for 5 minutes. 35 cycles of PCR contained secondary denaturation (45 seconds at 94°C), annealing (35 seconds at ‎‎55°C), primary extension (‎50 seconds at ‎72°C‎), and secondary extension (72°C for 7 minutes). DDW ‎was used in all PCR assays representing the negative control. ‎Also, 4 µL of PCR products was added to 1% agarose gel in the process of Electrophoresis, followed by ethidium bromide staining (45 minutes at 100 V). The bands were detected through a gel documentation system (KODAK gel logic 200 ‎imaging system).‎

2.3. Restriction Fragment Length Polymorphism Polymerase Chain Reaction (PCR-RFLP) and Sequence Analysis

All PCR amplicons were digested for 2 - 3 hours at 37°C using restriction endonucleases RsaI and HpaII and buffers advised by the manufacturer (Fermentas, Vilnius, Lithuania). Then, ITS1-rDNA gene underwent endonuclease reaction in a volume of 30 µL. This reaction consisted of 1 µL of RsaI with cut site GT↓AC and 1 µL of HpaII (cut site C↓CGG) 10 µL of PCR products, 3 µL of 10× buffer, and 16 µL of distilled water. Electrophoresis of restriction fragments of amplicons was performed using a 2% (w/v) agarose gel at 100 V for 60 minutes. Indeed, the identification of Echinococcus spp. through RFLP method was determined based on the recently reported standardized pattern (9). Furthermore, direct sequencing of 21 amplicons by targeting ITS1-rDNA gene in both directions via the BD1 and 4S primers by ABIPRISMTM 3130 genetic analyzer automated sequencer led to the re-confirmation of RFLP results (applied biosystem, USA). The standard IUPAC codes used for combinations of two or more bases contributed to the coding of ambiguous sites. Contigs from all samples were aligned, justified, and edited in consensus positions compared to GenBank sequences of all regional species, attempting to determine the similarity of new haplotypes by employing the Sequencher Tmv.4.1.4 Software for PC (gene codes corporation) (http://multalin.toulouse.inra.fr/multalin).

2.4. Mixed Sequence Reader

The MSR program (freely available at http://MSR.cs.nthu.edu.tw/) has been developed not only to identify heterogeneity in the chromatographic traces, but also to establish the physical positions of the detected variants in the mixed infections. The algorithm used in MSR was modified from that of Indelligent, but MSR is designed to use reference database alignment. The analytic steps used in MSR were described by Chang et al. in 2012. In order to import the DNA sequence chromatographs, the imported files were chromatography traces in the ABI format. The base peak positions, quality values, and four channel signals (A, C, G, and T) recorded in the .abi format file were extracted and analyzed to identify the major and the minor signals at each base location.

2.5. Phylogenetic and Haplotype Network Assessments

A wide range of haplotypes based on the sequences of ITS1-rDNA using the statistical parsimony method was also applied by TCS 1.2 software (26). In addition, the possible limit of network estimation was 95%. Confidence limits with a 95% confidence interval were obtained for infection rates. A Neighbor-Net network built in splits tree 4.0 served to analyze the phylogenetic information provided by ITS1-rDNA sequences (27) regarding the genetic distances measured as well as the Kimura-2 parameter model of nucleotide substitutions. Finally, the heterogeneity of haplotype, Hd, and Nucleotide (π) was estimated by DnaSP version 5.10 software.

3. Results

For all of the samples, a 1000 bp fragment was successfully amplified within ITS1-rDNA gene. Based on our in silico patterns, G1 (n = 58) was digested by restriction enzymes HpaII and RsaI into 300 bp, 700 bp and 345 bp, 655 bp fragments, respectively (Figure 1A and Figure 1B). However, no other expected genotype (such as G5) was distinguished in the cattle isolates. Furthermore, the digested fragments of mixed infections (n = 2) showed undefined patterns and were not supported by previously designed profiles (9). 21 out of 60 E. granulosus isolates (10 human, four sheep, four cattle, and three goat samples) were sequenced directly from ITS1-rDNA gene and determined firmly as corresponding to the 19 sheep strains (G1), two mixed infections. In addition, single-nucleotide variations (transition or transversion mutation models) were identified, including five unique haplotypes with a moderate range of diversity (Haplotype variability of 0.522 vs. Nucleotide variability of 0.02) in the constructed haplotype network. Two human (haplotype Mar 14), one sheep (haplotype Mar 11), one goat (haplotype Mar 07), and one cattle (haplotype Mar 09) isolates had seven nucleotide differences at positions 73, 74, 119, 160, and 312 bp. In addition, in an attempt to obtain the most likely major and minor variances, the BLAST analysis of overlapped chromatograms against a set of reference sequences showed that the isolated sequences of sheep strains belonged to infections with G1 and G6 genotypes (Figure 2A). The median mix ratio (major /minor) was 1.26. In addition, the LRi (log ratio of intensity) value for each sequence as the log ratio of the two intensities of the combined signal peaks was calculated with a cut-off of 2.69 (Figure 2B). Furthermore, the common haplotype Mar 08 (n = 16, frequency = 76%) was included without a notable heterogeneity in a consensus position. To identify a genealogical link within the haplotypes, we designed a statistical parsimony network. Figure 3 demonstrates all GenBank accession numbers for the inferred haplotypes of this study and for the reference genotypes/species used in the phylogenetic analysis.

The Obtained Patterns of Restriction Fragment Length Polymorphism (RFLP) Polymerase Chain Reaction (PCR) Digested by (A) RsaI and (B) HpaII Enzymes
M: DNA ladder marker (100 bp), -Ve: negative control; Restriction enzyme without PCR product, +Ve: positive control, understudied isolates: human (H), goat (G), cattle (Ct), and sheep (Sh).

M: DNA ladder marker (100 bp), -Ve: negative control; Restriction enzyme without PCR product, +Ve: positive control, understudied isolates: human (H), goat (G), cattle (Ct), and sheep (Sh).

Figure 1.

The Obtained Patterns of Restriction Fragment Length Polymorphism (RFLP) Polymerase Chain Reaction (PCR) Digested by (A) RsaI and (B) HpaII Enzymes

A, the scheme of mixed infection of G1 and G6 in overlapped chromatograms containing major and minor sequences; B, the LRi cut-off (log ratio of intensity) value for each sequence as the log ratio of the two intensities of the combined signal peaks.
A, the scheme of mixed infection of G1 and G6 in overlapped chromatograms containing major and minor sequences; B, the LRi cut-off (log ratio of intensity) value for each sequence as the log ratio of the two intensities of the combined signal peaks.
Figure 2.

A, the scheme of mixed infection of G1 and G6 in overlapped chromatograms containing major and minor sequences; B, the LRi cut-off (log ratio of intensity) value for each sequence as the log ratio of the two intensities of the combined signal peaks.

Neighbor-Net Graph Regarding the Kimura-2 Parameter Model and Sequences of Echinococcus spp. ITS1-rDNA Gene
Haplotypes Mar 07, 09, 11 and 14 in the G1 complex are marked by an asterisk (*).

Haplotypes Mar 07, 09, 11 and 14 in the G1 complex are marked by an asterisk (*).

Figure 3.

Neighbor-Net Graph Regarding the Kimura-2 Parameter Model and Sequences of Echinococcus spp. ITS1-rDNA Gene

4. Discussion

To date, one of the major difficulties in Echinococcus spp. taxonomic classification is related to the lake of accurate diagnosis of mixed infections in hyperendemic regions where different life cycles of Echinococcus spp. are sympatrically circulating with various heterogeneity traits and morphological aspects (21, 28). However, few studies have been carried out on the identification of Echinococcus spp. global concomitant infections (8, 29). In this study, an MSR program was used to directly analyze heterozygous base-calling chromatographs and subsequently detect the overlapped sequences and multiple structural variations, including SNPs (single nucleotide polymorphisms). The G1 and G6 were identified by MSR analysis showing that sheep and camel strains with novel haplotypes unequivocally exist in the region. These genotypes have been also reported in other central provinces of Iran (Semnan and Isfahan) (9); however, the type of mixed infections has not been determined precisely (30). Although G1 is the dominant strain in central Iran, the identification of both G1 and G6 genotypes in the region shows the presence of sheep-camel/dog life cycles in central Iran.

The liver and lung tissues of individual animals have shown the presence of the mixed infection caused by different haplotypes of E. granulosus spp. (8, 29). This can be arisen either from a single infection due to a definitive host concurrently harboring adult worms of the two haplotypes or genotypes, or from consecutive infections of the intermediate host. This study on 21 sequences of ITS1-rDNA gene led to the identification of merely five novel haplotypes in human, goat, cattle, and sheep isolates, while we expected to detect more haplotypes. In fact, this could be described by failing to use mitochondrial semi-conserved markers including cox 1 and nad 1. The ability of these markers in detecting more unknown haplotypes than nuclear-conserved gene (ITS1-rDNA) has already been proven (31, 32). The existence of novel haplotypes in G1 and mixed genotypes is related to the pathogenicity range of Echinococcus spp., the creation of re/emergent strains in the region, and the resistance of metacestodes to host innate immunity responses (33-35). Additionally, genotyping analyses of all cattle isolates exhibited the G1 strain. Nevertheless, we expected to detect other strains of E. granulosus such as E. ortoleppi since the cross-transmission of sheep and cattle genotypes of E. granulosus may occur in a wide range of hosts due to the overlapping cycles.

Several studies have attempted to genotype hydatid cysts in human, sheep, goat, cattle, buffalo, and camel populations in Iran on the basis of mitochondrial (cox1 and nad1) and nucleus markers, using RAPD-PCR, PCR-RFLP and sequencing strategies (19, 20, 22, 28, 30, 36, 37). However, none of them could provide remarkable applications in mixed genotypes. This issue authenticates the importance of discriminating overlapped E. granulosus genotypes existing in intermediate hosts to characterize the life cycle, host specificity, control programs, and classification of E. granulosus strains.

In summary, MSA, sequencing, and phylogenetic analyses showed that the common sheep strain (G1) and its concomitant infections with molecular diversity are unambiguously circulating in humans and livestock isolates in central Iran. The current results should be considered in improving the determination of the real epidemiological status of neglected Echinococcus spp. and launching disease control programs against echinococcosis previously understood life-cycle patterns. To appraise the multiple infection scenarios, further studies are needed focusing on cloning and multiplex real-time PCR in order to characterize the real taxonomic status of Echinococcus strains in wider areas of Iran and neighboring countries.

Acknowledgements

References

  • 1.

    Thompson RC, Jenkins DJ. Echinococcus as a model system: biology and epidemiology. Int J Parasitol. 2014;44(12):865-77. [PubMed ID: 25123067]. https://doi.org/10.1016/j.ijpara.2014.07.005.

  • 2.

    Thompson RC. Parasite zoonoses and wildlife: One Health, spillover and human activity. Int J Parasitol. 2013;43(12-13):1079-88. [PubMed ID: 23892130]. https://doi.org/10.1016/j.ijpara.2013.06.007.

  • 3.

    Thompson RC. The taxonomy, phylogeny and transmission of Echinococcus. Exp Parasitol. 2008;119(4):439-46. [PubMed ID: 18539274]. https://doi.org/10.1016/j.exppara.2008.04.016.

  • 4.

    Huttner M, Nakao M, Wassermann T, Siefert L, Boomker JD, Dinkel A, et al. Genetic characterization and phylogenetic position of Echinococcus felidis (Cestoda: Taeniidae) from the African lion. Int J Parasitol. 2008;38(7):861-8. [PubMed ID: 18062973]. https://doi.org/10.1016/j.ijpara.2007.10.013.

  • 5.

    Moks E, Jogisalu 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(5):647-54. [PubMed ID: 18261256]. https://doi.org/10.1017/S0031182008004198.

  • 6.

    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(Pt 5):713-22. [PubMed ID: 17156584]. https://doi.org/10.1017/S0031182006001934.

  • 7.

    Nakao M, Lavikainen A, Yanagida T, Ito A. Phylogenetic systematics of the genus Echinococcus (Cestoda: Taeniidae). Int J Parasitol. 2013;43(12-13):1017-29. [PubMed ID: 23872521]. https://doi.org/10.1016/j.ijpara.2013.06.002.

  • 8.

    Casulli A, Manfredi MT, La Rosa G, Cerbo AR, Genchi C, Pozio E. Echinococcus ortleppi and E. granulosus G1, G2 and G3 genotypes in Italian bovines. Vet Parasitol. 2008;155(1-2):168-72. [PubMed ID: 18514422]. https://doi.org/10.1016/j.vetpar.2008.04.004.

  • 9.

    Spotin A, Gholami S, Nasab AN, Fallah E, Oskouei MM, Semnani V, et al. Designing and conducting in silico analysis for identifying of Echinococcus spp. with discrimination of novel haplotypes: an approach to better understanding of parasite taxonomic. Parasitol Res. 2015;114(4):1503-9. [PubMed ID: 25645007]. https://doi.org/10.1007/s00436-015-4334-1.

  • 10.

    Dakkak A. Echinococcosis/hydatidosis: a severe threat in Mediterranean countries. Vet Parasitol. 2010;174(1-2):2-11. [PubMed ID: 20888694]. https://doi.org/10.1016/j.vetpar.2010.08.009.

  • 11.

    Babiker AM, Ravagnan S, Fusaro A, Hassan MM, Bakheit SM, Mukhtar MM, et al. Concomitant infection with Leishmania donovani and L. majorin single ulcers of cutaneous leishmaniasis patients from Sudan. J Trop Med. 2014;2014:1-8. https://doi.org/10.1155/2014/170859.

  • 12.

    Chang CT, Tsai CN, Tang CY, Chen CH, Lian JH, Hu CY, et al. Mixed sequence reader: a program for analyzing DNA sequences with heterozygous base calling. ScientificWorldJournal. 2012;2012:365104. [PubMed ID: 22778697]. https://doi.org/10.1100/2012/365104.

  • 13.

    Qvarnstrom Y, Visvesvara GS, Sriram R, da Silva AJ. Multiplex real-time PCR assay for simultaneous detection of Acanthamoeba spp., Balamuthia mandrillaris, and Naegleria fowleri. J Clin Microbiol. 2006;44(10):3589-95. [PubMed ID: 17021087]. https://doi.org/10.1128/JCM.00875-06.

  • 14.

    Hamzah Z, Petmitr S, Mungthin M, Leelayoova S, Chavalitshewinkoon-Petmitr P. Development of multiplex real-time polymerase chain reaction for detection of Entamoeba histolytica, Entamoeba dispar, and Entamoeba moshkovskii in clinical specimens. Am J Trop Med Hyg. 2010;83(4):909-13. [PubMed ID: 20889890]. https://doi.org/10.4269/ajtmh.2010.10-0050.

  • 15.

    Bowles J, McManus DP. Rapid discrimination of Echinococcus species and strains using a polymerase chain reaction-based RFLP method. Mol Biochem Parasitol. 1993;57(2):231-9. [PubMed ID: 8094539].

  • 16.

    Bowles J, Blair D, McManus DP. Genetic variants within the genus Echinococcus identified by mitochondrial DNA sequencing. Mol Biochem Parasitol. 1992;54(2):165-73. [PubMed ID: 1435857].

  • 17.

    Spotin A, Rouhani S, Parvizi P. The associations of Leishmania major and Leishmania tropica aspects by focusing their morphological and molecular features on clinical appearances in Khuzestan province, Iran. Biomed Res Int. 2014;2014:1-13. https://doi.org/10.1155/2014/913510.

  • 18.

    Mahami Oskouei M, Ghabouli Mehrabani N, Miahipour A, Fallah E. Molecular characterization and sequence analysis of Echinococcus granulosus from sheep isolates in East Azerbaijan province, northwest of Iran. J Parasit Dis. 2016;40(3):785-90. [PubMed ID: 27605785]. https://doi.org/10.1007/s12639-014-0579-3.

  • 19.

    Parsa F, Fasihi Harandi M, Rostami S, Sharbatkhori M. Genotyping Echinococcus granulosus from dogs from Western Iran. Exp Parasitol. 2012;132(2):308-12. [PubMed ID: 22884512]. https://doi.org/10.1016/j.exppara.2012.07.010.

  • 20.

    Pezeshki A, Akhlaghi L, Sharbatkhori M, Razmjou E, Oormazdi H, Mohebali M, et al. Genotyping of Echinococcus granulosus from domestic animals and humans from Ardabil Province, northwest Iran. J Helminthol. 2013;87(4):387-91. [PubMed ID: 23046636]. https://doi.org/10.1017/S0022149X1200051X.

  • 21.

    Ahmadi N, Dalimi A. Characterization of Echinococcus granulosus isolates from human, sheep and camel in Iran. Infect Genet Evol. 2006;6(2):85-90. [PubMed ID: 16503509]. https://doi.org/10.1016/j.meegid.2005.01.005.

  • 22.

    Jamali R, Ghazanchaei A, Asgharzadeh M. Identification and characterization of Echinococcus granulosus by PCR RFLP technique in Tabriz district. J Parasitic Dis. 2004;28(2):69-72.

  • 23.

    Karimi A, Dianatpour R. Genotypic and phenotypic characterization of Echinococcus granulosus of Iran. Biotechnol. 2008;7(4):757-62. https://doi.org/10.3923/biotech.2008.757.762.

  • 24.

    Rajabloo M, Hosseini SH, Jalousian F. Morphological and molecular characterisation of Echinococcus granulosus from goat isolates in Iran. Acta Trop. 2012;123(2):67-71. [PubMed ID: 22484229]. https://doi.org/10.1016/j.actatropica.2012.03.006.

  • 25.

    Smyth JD, McManus D, Barrett NJ, Bryceson A, Cowie AGA. In vitro culture of human hydatid material. Lancet. 1980;315(8161):202-3. https://doi.org/10.1016/s0140-6736(80)90685-6.

  • 26.

    Clement M, Posada D, Crandall KA. TCS, a computer program to estimate gene genealogies. Mol Ecol. 2000;9(10):1657-9. [PubMed ID: 11050560].

  • 27.

    Huson DH, Bryant D. Application of phylogenetic networks in evolutionary studies. Mol Biol Evol. 2006;23(2):254-67. [PubMed ID: 16221896]. https://doi.org/10.1093/molbev/msj030.

  • 28.

    Youssefi MR, Tabaripour R, Omrani VF, Spotin A, Esfandiari B. Genotypic characterization of echinococcus granulosus in Iranian goats. Asian Pac J Trop Dis. 2013;3(5):362-6. https://doi.org/10.1016/s2222-1808(13)60085-7.

  • 29.

    Umhang G, Richomme C, Boucher JM, Hormaz V, Boue F. Prevalence survey and first molecular characterization of Echinococcus granulosus in France. Parasitol Res. 2013;112(4):1809-12. [PubMed ID: 23271567]. https://doi.org/10.1007/s00436-012-3245-7.

  • 30.

    Shahnazi M, Hejazi H, Salehi M, Andalib AR. Molecular characterization of human and animal Echinococcus granulosus isolates in Isfahan, Iran. Acta Trop. 2011;117(1):47-50. [PubMed ID: 20858453]. https://doi.org/10.1016/j.actatropica.2010.09.002.

  • 31.

    Esposti MD, De Vries S, Crimi M, Ghelli A, Patarnello T, Meyer A. Mitochondrial cytochrome b: evolution and structure of the protein. Biochim Biophys Acta. 1993;1143(3):243-71. [PubMed ID: 8329437].

  • 32.

    Howell N, Gilbert K. Mutational analysis of the mouse mitochondrial cytochrome b gene. J Mol Biol. 1988;203(3):607-17. https://doi.org/10.1016/0022-2836(88)90195-7.

  • 33.

    Mokhtari Amirmajdi M, Sankian M, Eftekharzadeh Mashhadi I, Varasteh A, Vahedi F, Sadrizadeh A, et al. Apoptosis of human lymphocytes after exposure to hydatid fluid, Mashhad university of medical sciences. Iran J Parasitol. 2011;6(2):9-16.

  • 34.

    Spotin A, Majdi MM, Sankian M, Varasteh A. The study of apoptotic bifunctional effects in relationship between host and parasite in cystic echinococcosis: a new approach to suppression and survival of hydatid cyst. Parasitol Res. 2012;110(5):1979-84. [PubMed ID: 22167369]. https://doi.org/10.1007/s00436-011-2726-4.

  • 35.

    Kamenetzky L, Canova SG, Guarnera EA, Rosenzvit MC. Echinococcus granulosus: DNA extraction from germinal layers allows strain determination in fertile and nonfertile hydatid cysts. Exp Parasitol. 2000;95(2):122-7. [PubMed ID: 10910713]. https://doi.org/10.1006/expr.2000.4518.

  • 36.

    Hajialilo E, Harandi MF, Sharbatkhori M, Mirhendi H, Rostami S. Genetic characterization of Echinococcus granulosus in camels, cattle and sheep from the south-east of Iran indicates the presence of the G3 genotype. J Helminthol. 2012;86(3):263-70. [PubMed ID: 21749740]. https://doi.org/10.1017/S0022149X11000320.

  • 37.

    Shariatzadeh SA, Spotin A, Gholami S, Fallah E, Hazratian T, Mahami Oskouei M, et al. The first morphometric and phylogenetic perspective on molecular epidemiology of Echinococcus granulosus sensu lato in stray dogs in a hyperendemic Middle East focus, northwestern Iran. Parasit Vectors. 2015;8:409. [PubMed ID: 26242729]. https://doi.org/10.1186/s13071-015-1025-9.