Abstract
Hydatid cyst caused by Echinococcus granulosus is one of the most important parasitic diseases around the world and many countries in Asia, including Iran, are involved with this infection. This disease can cause high mortality in humans as well as economic losses in livestock. To date, several molecular methods have been used to determine the genetic diversity of E. granulosus. So far, identification of E. granulosus using real-time PCR fluorescence-based quantitative assays has not been studied worldwide, also in Iran. Therefore, the aim of this study was to investigate the genetic diversity of E. granulosus from center of Iran using real-time PCR method. A total of 71 hydatid cysts were collected from infected sheep, goat, and cattle slaughtered in Isfahan, Iran during 2013. DNA was extracted from protoscolices and/or germinal layers from each individual cyst and used as template to amplify the mitochondrial cytochrome c oxidase subunit 1 gene (cox1) (420 bp). Five cattle isolates out of 71 isolates were sterile and excluded from further investigation. Overall, of 66 isolates, partial sequences of the cox1 gene of E. granulosus indicated the presence of genotypes G1 in 49 isolates (74.2%), G3 in 15 isolates (22.7%), and G6 in 2 isolates (3.0%) in infected intermediate hosts. Sixteen sequences of G1 genotype had microgenetic variants, and they were compared to the original sequence of cox1. However, isolates identified as G3 and G6 genotypes were completely consistent with original sequences. G1 genotype in livestock was the dominant genotype in Isfahan region, Iran.
-
Key words: Echinococcus granulosus, cox1 gene, genotyping
INTRODUCTION
Hydatid cyst, caused by
Echinococcus granulosus is one of the most important parasitic diseases worldwide and distribute in many Asia countries, including Iran. This disease can cause high mortality in humans as well as economic losses in livestock [
1,
2,
3]. Nowadays, genetic studies on a partial sequence of mitochondrial cytochrome
c oxidase subunit 1 (
cox1) and nuclear genetic markers suggest that
E. granulosus has 10 distinct genotypes termed G1-G10 [
4,
5,
6,
7,
8,
9]. These different genotypes include G1 and G2 as sheep strains, G3 and G5 as bovine strains, G4 and G6 as horse and camel strains, G7 as a pig strain, G8 and G10 as cervid strains [
10]. Recent taxonomic studies have suggested
E. granulosus in 4 distinct groups, including
E. granulosus sensu stricto (G1-G3 genotypes),
Echinococcus equinus (G4),
Echinococcus ortleppi (G5), and
Echinococcus canadensis (G6-G10) [
11]. The identification of
E. granulosus strains has been carried out in different laboratories using various analytical methods (morphology, physiology, biochemistry, and molecular genetics), and all of them were proved to be useful, particularly when used together [
12,
13,
14,
15,
16]. The possible susceptibility factors for genetic variation that occur in
E. granulosus may influence many different phenotypic features such as the life cycle patterns, development rate, host specificity, geographical distribution, transmission dynamics, infectivity, antigenicity, control of diseases, and chemotherapeutic agents.
PCR-based methods have been extensively used to characterize the different strains of
E. granulosus [
12]. On the other hand, various techniques such as PCR-RFLP, nested PCR, and real-time PCR were used to determine genetic variations within
E. granulosus. However, fluorescence-based quantitative PCR assays has not fully developed yet [
17,
18]. Real-time PCR is a combination of PCR chemistry with fluorescent probe for detection of amplicon which shows results with curve [
19]. Specifically, EVA Green™-based real-time PCR, in association with melting curve analysis, can be useful for detection of nucleotide variation such as single nucleotide polymorphism (SNP). Different PCR products often have a different melting temperature (T
m), and this T
m is related to the percentage of guanine-cytosine (GC), length, and nucleotide composition [
20]. In Iran, especially in Isfahan, sheep is the most common and important intermediate host for
E. granulosus. According to previous studies, the prevalence of infection and cyst fertility rates were higher in sheep than in other slaughtered animals [
21,
22,
23]. None of the previous studies in Iran was carried out using real-time PCR method for detection of
E. granulosus. Therefore, the aim of this study was to investigate the genetic diversity of
E. granulosus from center of Iran using real-time PCR method.
MATERIALS AND METHODS
Sample collection
Seventy-one hydatid cysts were collected from infected slaughtered animals in different regions of Isfahan, including Fasaran, Khomeinishahr, and Najafabad from January 2013 until March 2013. Separation of hydatid cysts was performed under sterile condition, and then protoscolices and/or the germinal layer were collected from an individual hydatid cyst. Protoscolices and the germinal layer were preserved in 70% (v/v) ethanol and stored at -20℃ until DNA extraction.
DNA extraction
Before DNA extraction, ethanol was removed by 3 times washing with sterile distilled water. Total genomic DNA was extracted using a genomic DNA extraction kit (Bioneer, Daejeon, Korea) according to the manufacturer's instructions with some modifications. The concentration of extracted DNA was evaluated by Nano Drop, and then the samples were kept at -20℃ for further analysis.
PCR and DNA sequencing
Genomic DNA was amplified using specific primers for
cox1 gene, i.e., forward JB3 primer (5-TTT TTT GGG CAT CCT GAG GTT TAT-3) and reverse JB4.5 primer (5-TAA AGA AAG AAC ATA ATG AAA ATG-3'] as described by Bowles et al. [
5]. PCR was done in a 25 µl final volume containing 10 µl master mixes (Type-it HRM PCR Kit, Qiagen, Hilden, Germany), 9 µl distilled water, 1 µl each primer, and 4 µl template DNA.
Enzymatic reaction was performed under the following conditions: an initial denaturation step at 95℃ for 10 min, followed by 40 cycles at 95℃ for 10 sec, 55℃ for 30 sec, 72℃ for 27 sec, and a final extension step at 72℃ for 5 min.
The master mix provided with the kit contains the novel double-stranded DNA-binding fluorescent dye, Eva Green, and includes an optimized HRM buffer, Hot StarTaq plus DNA polymerase, Q-Solution, and dNTPs. Together, these components provided reliable PCR specificity and acceptable results, even with hard genomic loci. Real-time PCR was carried out in a Mini Opticon real-time PCR detection system (Rotor-Gene 6,000, Hilden, Germany) under the following conditions: after a pre-incubation step (95℃ for 3 min) to activate the polymerase enzyme, 40 cycles of amplification were performed, each one including 95℃ for 10 sec, 55℃ for 30 sec, and 72℃ for 27 sec, followed by a final extension at 72℃ for 5 min. The fluorescence signal was measured at each cycle after the extension step. The melting experiment was performed from 55℃ to 95℃ at 1℃/sec with continuous fluorescence monitoring. Melting peaks were visualized by plotting the first derivate against the melting temperature; the Tm was defined as the peak of the curve. The quantitative detection and Tm were obtained using the Rotor-gene 6,000 series software version 1.7 (Corbett, Hilden, Germany).
Two replicates of each sample were analyzed and Tm analysis was repeated 3 times in each run to confirm the repeatability of the Tm assay by estimating the Tm variation within PCR amplification (intra-assay) and between PCR amplification (inter-assay), respectively. The coefficient of variation (CV) was calculated dividing the standard deviation by the arithmetic mean of the measured values of Tm (CV=standard deviation [SD]/mean value).
In addition, to check the uniformity of temperature of the cycler block, indicating the reliability of results, replicates of an already genotyped sample were amplified in different positions of the cycler block within the same amplification cycle. The intra-assay CVs represent the mean CVs for the results obtained from the replicates of all the E. granulosus genotypes in the same run. The inter-assay CVs represent the mean CVs for the results obtained from 3 separate runs.
To confirm the identified genotypes, 22 samples were randomly sequenced for cox1 mtDNA. Three samples that have been already sequenced for cox1 and identified as G1, G3, and G6 (accession nos. HM563013, HM563017, and HM563018, respectively) were included in each PCR set as positive controls.
Sequence homology and phylogenetic analysis
The obtained continuous sequences were compared with previously published sequences of the mitochondrial
cox1 for
E. granulosus genotypes in NCBI using BLAST system. Phylogenetic analyses of the sequence data were inferred with maximum likelihood and multiple alignments and done using MEGA 5 of the computer program CLC main workbench software (version 5.2.1, 2013) [
24,
25,
26].
RESULTS
Out of 71 isolates, 5 cattle isolates were sterile and excluded from further investigation. Molecular survey and PCR amplification were successfully performed on all 66 isolates. Of 66 samples, 51 isolates were collected from sheep (35 livers and 16 lungs), 8 from goat (4 livers and 4 lungs), and 7 cattle origin (lungs). Real-time PCR using the partial sequences of cox1 gene of 66 E. granulosus isolates showed that 49, 15, and 2 isolates were identified as G1, G3, and G6 genotypes, respectively. Sheep, goat, and cattle isolates were genotypically categorized as G1 (77%, 63%, and 72%), G3 (21%, 25%, and 28%), and G6 (2%, 12%, and 0%), respectively (data not shown).
As shown in
Table 1, real-time PCR melting curve results indicated that the mean T
m of G1, G3, and G6 genotypes were 79.3, 78.8, and 77.3, respectively. Assessment of intra and inter-assay variability showed low and acceptable CVs (
Table 1). Melting curve and T
m of the analyzed hydatid cysts identification by sequencing are shown in
Fig. 1.
Sequencing results confirmed that all isolates were correctly differentiated by real-time PCR. All the isolates identified by real-time PCR were clustered along with the corresponding reference genotypes as shown in
Fig. 2. The phylogenetic tree was divided into 2 main clades, the first clade contained 2 subclades corresponding to the G1 and G3 genotypes (
E. granulosus sensu stricto), and the second one corresponding to the isolates identified as G6 (
E. canadensis) along with the individual reference genotypes. Forty-nine samples with T
m between 79.1℃ and 79.5℃ had identical sequences to the G1genotype, while the 15 samples with T
m between 78.6℃ and 79.0℃ had identical sequences to the G3 genotype, and finally 2 samples with T
m between 77.1℃ and 77.5℃ had identical sequences to the G6 genotype.
Phylogenetic analysis also showed that 16 isolates of G1 genotype had micro-variants, compared to the original sequences of cox1 gene (accession no. M84661). Four samples identified as G3 and 2 samples identified as G6 were completely identical to the G3 and G6 reference sequences (accession nos. M84663 and M84666, respectively).
DISCUSSION
In the present study, Tm analysis was used for detection of E. granulosus strains (G1, G3, and G6 genotypes) in Isfahan Province, Iran. Due to high sensitivity, specificity, accuracy, and repeatability of this novel method, we confirmed that this technique is useful for detection of Iranian genotypes. Because of 0.5℃ mean difference of Tm between G1 and G3 genotypes, the Tm ranges did not overlap, therefore we could differentiate these 2 genotypes easily and rapidly. On the other hand, G6 genotype could be easily distinguished from the other ones (G1 and G3), because of 1.5℃ different Tm. The accuracy of the Tm assay was confirmed by comparing these results with direct sequencing.
The advantage of using real-time PCR is that this technique is more cost-effective and less time-consuming than conventional PCR. Actually, with this method, a large number of samples could be tested quickly and repeatable results could be achieved with significant saving in time and materials [
17,
27,
28]. On the other hand, this method has some disadvantages, such as existing instruments with limited multiplexing capabilities, incompatibility of several systems with some fluorogenic chemistries, and necessity of high technical skills for its launching [
27].
Epidemiological studies showed that G1 genotype is the predominant genotype in both human and animals [
29]. Maurelli et al. [
17] used real-time PCR for identification of G1 and G2/G3 genotypes of
E. granulosus in Italy. They found that T
m analysis can be used for discriminating between G1 and G2/G3 genotypes of
E. granulosus sensu stricto. Our results were compatible with their findings.
Recently, Boubaker et al. [
30] in Switzerland developed a single-tube multiplex PCR (mPCR) for discrimination of the 3 levels of the
E. granulosus complex genotyping. They found
E. granulosus sensu stricto (G1/G2/G3),
E. granulosus equinus (G4), and
E. canadensis (G7) as predominant genotypes in their study. This study showed that this method could provide the unique opportunity to address directly speciation and genotyping within the framework of large-scale studies.
Recent studies conducted analysis of various strains in the Middle East. For example, a study by Adwan et al. [
12] in Palestine showed the presence of G1-G3 genotypes in this region using
cox1 loci. In another study, Simsek et al. [
31] reported a G1/G3 cluster (
E. granulosus sensu stricto) in Turkey.
In recent years, several studies reported different genotypes of
E. granulosus from intermediate and definitive hosts, including the sheep, goat, cattle, camel, buffalo, human, and dog in different parts of Iran [
4,
22,
32,
33]. The results of this study were similar to previous studies in Iran showing that genotypes G1, G3, and G6 are present in intermediate hosts [
29]. In the few studies, genotypes G1, G3, and G6 were simultaneously reported [
4,
34]. Recently, Rostami et al. [
29] used the post real-time PCR with HRM method and differentiated G1, G3, and G6 genotypes of sensu lato in Iran. Rajabloo et al. [
35] reported G1 and G6 genotypes from goats in Iran. A similar study in Ilam Province on hydatid cysts isolated from animals and humans showed
E. granulosus sensu stricto (
G1-
G3) in this region [
36]. A study on 120 hydatid cysts isolated from goats in Mazandaran Province, Northern Iran, by Youssefi et al. [
37] showed the G1 genotype as the dominant strain. Among various studies, Rostami et al. [
29] reported G1, G3, and G6 genotypes in the intermediate hosts that are similar to our results.
In conclusion, we believe that real-time PCR that has been developed in the present study, provides powerful means for molecular and epidemiological studies on human and animal E. granulosus infections. Therefore, we suggest that the extensive use of this technique with high sample size can help to find undetected genotypes using different molecular targets.
Isfahan University of Medical Sciences
Notes
-
The authors declare no conflict of interest related to this work.
ACKNOWLEDGMENTS
This study was carried out as a part of M.Sc. thesis done by Ahmad Hosseini Safa and was financially supported by the vice-chancellor for Research, Isfahan University of Medical Sciences, Iran.
References
- 1. Wang Z, Xiaoming W, Xiaoqing L. Echinococcosis in China, a review of the epidemiology of Echinococcus spp. Ecohealth 2008;5:115-126.
- 2. Sadjjadi SM. Present situation of echinococcosis in the Middle East and Arabic North Africa. Parasitol Int 2006;55:S197-S202.
- 3. Torgerson PR, Oguljahan B, Muminov AE, Karaeva RR, Kuttubaev OT, Aminjanov M, Shailkenov B. Present situation of cystic echinococcosis in Central Asia. Parasitol Int 2006;55:S207-S212.
- 4. Sharbatkhori M, Harandi MF, Mirhendi H, Hajialilo E, Kia EB. Sequence analysis of cox1 and nad1 genes in Echinococcus granulosus G3 genotype in camels (Camelus dromedarius) from central Iran. Parasitol Res 2011;108:521-527.
- 5. Bowles J, Blair D, McManus DP. Genetic variants within the genus Echinococcus identified by mitochondrial DNA sequencing. Mol Biochem Parasitol 1992;54:165-173.
- 6. 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.
- 7. Scott JC, Stefaniak J, Pawlowski Z, McManus D. Molecular genetic analysis of human cystic hydatid cases from Poland: identification of a new genotypic group (G9). Parasitology 1997;114:37-43.
- 8. Thompson RC, McManus DP. Towards a taxonomic revision of the genus Echinococcus. Trends Parasitol 2002;18:452-457.
- 9. Lavikainen A, Lehtinen M, Meri T, Hirvela-Koski V, Meri S. Molecular genetic characterization of the Fennoscandian cervid strain, a new genotypic group (G10) of Echinococcus granulosus. Parasitology 2003;127:207-215.
- 10. Grosso G, Gruttadauria S, Biondi A, Marventano S, Mistretta A. Worldwide epidemiology of liver hydatidosis including the Mediterranean area. World J Gastroenterol 2012;18:1425.
- 11. Nakao M, Yanagida T, Okamoto M, Knapp J, Nkouawa A, Sako Y, Ito A. State-of-the-art Echinococcus and Taenia: Phylogenetic taxonomy of human-pathogenic tapeworms and its application to molecular diagnosis. Infect Genet Evol 2010;10:444-452.
- 12. Adwan G, Adwan KM, Bdir S, Abuseir S. Molecular characterization of Echinococcus granulosus isolated from sheep in Palestine. Exp parasitol 2013;134:195-199.
- 13. Singh BB, Sharma JK, Tuli A, Sharma R, Bal MS, Aulakh RS, Gill JPS. Prevalence and morphological characterisation of Echinococcus granulosus from north India. J Parasit Dis 2014;38:36-40.
- 14. Cui SJ, Xu LL, Zhang T, Xu M, Yao J, Fang CY, Feng Z, Yanga PY, Hu W, Liu F. Proteomic characterization of larval and adult developmental stages in Echinococcus granulosus reveals novel insight into host-parasite interactions. J Proteomics 2013;84:158-175.
- 15. Shield JM. Dipylidium caninum, Echinococcus granulosus and Hydatigera taeniaeformis: Histochemical identification of cholinesterases. Exp Parasitol 1969;25:217-231.
- 16. McManus DP. Immunodiagnosis of sheep infections with Echinococcus granulosus: in 35 years where have we come? Parasite Immunol 2014;36:125-130.
- 17. Maurelli MP, Rinaldi L, Capuano F, Perugini AG, Cringoli G. Development of a real-time PCR for the differentiation of the G1 and G2/G3 genotypes of Echinococcus granulosus. Parasitol Res 2009;105:255-259.
- 18. Parija SC. Hydatid fluid as a clinical specimen for the etiological diagnosis of a suspected hydatid cyst. J Parasit Dis 2004;28:64-68.
- 19. Espy MJ, Uhl J, Sloan L, Buckwalter S, Jones M, Vetter EA, Yao JD, Wengenack NL, Rosenblatt JE, Cockerill FR 3rd, Smith TF. Real-time PCR in clinical microbiology: applications for routine laboratory testing. Clin Microbiol Rev 2006;19:165-256.
- 20. Herrmann MG, Durtschi JD, Bromley LK, Wittwer CT, Voelkerding KV. Amplicon DNA melting analysis for mutation scanning and genotyping: cross-platform comparison of instruments and dyes. Clin Chem 2006;52:494-503.
- 21. Rostami Nejad M, Taghipour N, Nochi Z, Mojarad EN, Mohebbi S, Harandi MF, Zali MR. Molecular identification of animal isolates of Echinococcus granulosus from Iran using four mitochondrial genes. J Helminthol 2012;86:485-492.
- 22. Shahnazi M, Hejazi H, Salehi M, Andalib AR. Molecular characterization of human and animal Echinococcus granulosus isolates in Isfahan, Iran. Acta Trop (Basel) 2011;117:47-50.
- 23. Zhang L, Eslami A, Hosseini S, McManus D. Indication of the presence of two distinct strains of Echinococcus granulosus in Iran by mitochondrial DNA markers. Am J Trop Med Hyg 1998;59:171-174.
- 24. Kimura M. A simple method for estimating evolutionary rates of base substitutions through comparative studies of nucleotide sequences. J Mol Evol 1980;16:111-120.
- 25. Kumar S, Nei M, Dudley J, Tamura K. MEGA: a biologist-centric software for evolutionary analysis of DNA and protein sequences. Brief Bioinform 2008;9:299-306.
- 26. 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.
- 27. Morsczeck C, Langendörfer D, Schierholz JM. A quantitative real-time PCR assay for the detection of tetR of Tn 10 in Escherichia coli using SYBR Green and the Opticon. J Biochem Biophys Methods 2004;59:217-227.
- 28. Harasawa R, Mizusawa H, Fujii M, Yamamoto J, Mukai H, Uemori T, Asada K, Kato I. Rapid detection and differentiation of the major mycoplasma contaminants in cell cultures using real-time PCR with SYBR Green I and melting curve analysis. Microbiol Immunol 2005;49:859-863.
- 29. McManus DP, Thompson R. Molecular epidemiology of cystic echinococcosis. Parasitology 2003;127:S37-S51.
- 30. Boubaker G, Macchiaroli N, Prada L, Cucher MA, Rosenzvit MC, Ziadinov I, Deplazes P, Saarma U, Babba H, Gottstein B, Spiliotis M. A Multiplex PCR for the Simultaneous Detection and Genotyping of the Echinococcus granulosus Complex. PLoS Negl Trop Dis 2013;7:e2017.
- 31. Simsek S, Kaplan M, Ozercan IH. A comprehensive molecular survey of Echinococcus granulosus in formalin-fixed paraffin-embedded tissues in human isolates in Turkey. Parasitol Res 2011;109:411-416.
- 32. Pour AA, Hosseini SH, Shayan P. Comparative genotyping of Echinococcus granulosus infecting buffalo in Iran using cox1 gene. Parasitol Res 2011;108:1229-1234.
- 33. Parsa F, Harandi MF, Rostami S, Sharbatkhori M. Genotyping Echinococcus granulosus from dogs from Western Iran. Exp Parasitol 2012;132:308-312.
- 34. Rostami S, Talebi S, Babaei Z, Sharbatkhori M, Ziaali N, Rostami H, Harandi MF. High resolution melting technique for molecular epidemiological studies of cystic echinococcosis: differentiating G1, G3, and G6 genotypes of Echinococcus granulosus sensu lato. Parasitol Res 2013;112:3441-3447.
- 35. Rajabloo M, Hosseini SH, Jalousian F. Morphological and molecular characterisation of Echinococcus granulosus from goat isolates in Iran. Acta Trop 2012;123:67-71.
- 36. Dousti M, Abdi J, Bakhtiyari S, Mohebali M, Mirhendi S, Rokni M. Genotyping of hydatid cyst isolated from human and domestic animals in Ilam Province, Western Iran using PCR-RFLP. Iran J Parasitol 2013;8:47-52.
- 37. 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:362-366.
- 38. Nakao M, McManus D, Schantz P, Craig P, Ito A. A molecular phylogeny of the genus Echinococcus inferred from complete mitochondrial genomes. Parasitology 2007;134:713-722.
Fig. 1Tm of the analyzed hydatid cysts. (A) E. granulosus G1 identified by sequencing, (B) E. granulosus G3 identified by sequencing, and (C) E. granulosus G6 identification by sequencing.
Fig. 2Molecular phylogenetic tree of 22
E. granulosus isolates of sheep, goat, and cow along with reference isolates based on
CO1 gene sequence. The evolutionary history was inferred by using the Maximum Likelihood method based on the Kimura 2 parameter model [
19]. The tree with the highest log likelihood (-814.2241) is shown. Reference accession nos.: G1, M64661; G3, M64663; G6, M84666.
Table 1.Mean Tm (℃) and SD calculated and intra- and inter-assay CVs of the G1, G3, and G6 genotypes of E. granulosus
Table 1.
|
E. granulosus genotype |
Mean Tm (°C) |
SD |
Intra-assay CV (%) |
Inter-assay CV (%) |
|
G1 |
79.3 |
0.15 |
0.07 |
0.12 |
|
G3 |
78.8 |
0.14 |
0.09 |
0.14 |
|
G6 |
77.3 |
0.12 |
0.05 |
0.12 |
