Original Article
Iranian J Publ Health, Vol. 43, No.10, Oct 2014, pp.1364-1371
Genotyping and Phylogenetic Analysis of Fasciola Spp. Isolated from Sheep and Cattle Using PCR-RFLP in Ardabil Province, Northwestern Iran Mojgan ARYAEIPOUR 1, Soheila ROUHANI 2, Mojgan BANDEHPOUR 2, Hadi MIRAHMADI 1, *Bahram KAZEMI 2, *Mohammad Bagher ROKNI 3, 4 1. Dept. of Medical Parasitology and Mycology, International Branch, Shahid Beheshti University of Medical Science, Tehran, Iran 2. Dept. Of Biotechnology, School of Medicine, Shahid Beheshti University of Medical Science, Tehran, Iran 3. Dept. of Medical Parasitology and Mycology, School of Public Health, Tehran University of Medical Sciences, Tehran, Iran 4. Center for Research of Endemic Parasites of Iran (CREPI), Tehran University of Medical Sciences, Tehran, Iran *Corresponding Authors: E-mail: [email protected], [email protected] (Received 17 Jun 2014; accepted 21 Sep 2014)
Abstract Background: The aim of this study was to detect the genotype of Fasciola spp. in Meshkin-Shahr, Ardabil Province, northwestern Iran in different hosts using PCR-RFLP. Methods: The parasite hosts included cattle, and sheep. Overall, 70 adult flukes from livers of slaughtered animals were collected from the abattoirs of aforementioned area. The included 35 samples from infected sheep and 35 samples from 35 infected cattle. PCR–RFLP and sequence analysis of the first nuclear ribosomal internal transcribed spacer (ITS 1) region from Fasciola species were used to conduct the study. Results: The fragment of approximately 700bp in all of the Fasciola samples was amplified. PCR products of ITS 1 were subjected for digestion by restriction enzyme. RsaI restriction enzyme was selected for RFLP method that caused the separation specifically of Fasciola species. Amplicons with the sequences of F. hepatica had a pattern of about 360, 100, and 60 bp band size, whereas F. gigantica worms had a profile of 360, 170, and 60 bp in size, respectively. Results based on PCR-RFLP analysis were confirmed by sequence analysis of representative ITS 1 amplicons. No hybrid forms were detected in the present study. All sheep were infected with F. hepatica but cattle were infected with both species. Conclusion: Both species of Fasciola are present in Ardabil. The method described here can be valuable for identification of Fasciola species in endemic parts for fasciolosis, regions with intermediate species and in that overlapping distribution area. Keywords: Fasciola, Genotyping, Fasciolosis, PCR, Iran
Introduction Fascioliasis, caused by Fasciola spp., is one of the most zoonotic diseases with a global extension and of a significant concern in areas with suitable situations (1).” About 2.4-17 million cases of human fascioliasis are estimated in the world and 180 million exposed at risk of fascioliasis that demonstrates the importance of the disease (2)”. 1364
Two great epidemics of human fascioliasis with about 10,000 people infected in each case occurred in the north of Iran (3-5). According to WHO, fascioliasis is one of the most important diseases and parasitic infections transmitted by food (6). Considering the importance of fascioliasis in Iran, various studies especially in order to Available at:
http://ijph.tums.ac.ir
Aryaeipour et al.: Genotyping and Phylogenetic Analysis of Fasciola Spp. …
identify and genotyping of Fasciola seems indispensable. Morphological features of adult worms and eggs are influenced by the different factors such as host type, age of parasite, fixation of the samples and intensity of the infection (7(.Distribution overlap of F. hepatica and F. gigantica in Asia and simultaneous infection with these two parasites in cattle in northern Iran and bug reporting interface is also introduced in other provinces of Iran (8). It is usually difficult to accurate differentiation between Fasciola species based on phenotype criteria (9-11). Today, molecular methods based on DNA analysis is a great impact in many fields, including identifying systematic parasitology, parasite infection, parasite epidemic, analysis of the genetic structure of parasite populations, such as genetic differences between the sexes, gene expression and drug resistance of and the vaccine.Therefore, it seems that the use of the accurate and reliable method for identification and differentiation of Fasciola species is necessary. For this purpose, PCR-RFLP based on 28s rRNA, 18s rRNA, ITS 1 or ITS 2 has been used (12-13). In the present study, region between 18S and 28S (ITS 1, 5.8s,) of ribosomal DNA was used by PCR-RFLP method. Since, ITS 1 sequence was suitable genetic markers for genotyping, interaspecific variations, and phylogenetic studies of parasites (14); we used only ITS 1 for molecular identification of Fasciola species.The adjacent of Ardabil Province to endemic areas of fasciolosis persuaded us to conduct the present study (15).
Materials and Methods Parasite
Overall, 70 Fasciola isolates were studied from two hosts in Ardabil Province, northwestern Iran (Fig. 1). These parasites were isolated from 35 infected sheep and 35 infected cattle. Individual worms were washed extensively in PBS buffer (37°C) and subsequently fixed in 20% ethanol and maintained at room temperature for several weeks until extraction of genomic DNA (Table 1).
Available at: http://ijph.tums.ac.ir
Fig. 1: Map of Iran (a) Ardabil Province, northwestern Iran (b)
DNA extraction
Approximately 10 mg tissue was removed from a portion of the lateral zone of adult flukes and crushed. The ethanol in each of these samples was allowed to evaporate for a few minutes and washed in distilled water for 3 times and squashing before the genomic DNA was extracted with (DNGTM-PLUS) kit was used for genomic DNA extraction of Fasciola Parasites. The extracted DNA, split into aliquots, and stored at -20°C until used in the PCR.
Primer designing
We designed forward (5’-ACCGGTGCTGAGAAGACG-3’) and Reverse (5’-CGACGTACGTGCAGTCCA3’) primers using Gene Runner and BLAST
(http://www.ncbi.nlm.nih.gov/blast) software. The primers were synthesized by Bioneer Company (Korea).
PCR
To amplify a 700 bp region of the ITS 1 sequence, PCR was performed using a pair primer. Total volume of reaction was 15µl containing 1.5µl DNA template, 5µl distilled water, 10pmol of each primers (Forward and Reverse), and 7.5µl master mix (amplicon). PCR amplification was performed in Eppendorf Mastercycler Gradient thermocycler . The temperature profile was as follows: an initiation of 95ºC for 5 min, followed by 30 cycles of 94ºC for 30s (denaturation), 60ºC for 30s (annealing), 72ºC for 30s (extension) and a final extension of 72ºC for 5min followed by cooling at 4ºC. 1365
Iranian J Publ Health, Vol. 43, No.10, Oct 2014, pp. 1364-1371
Table 1: Position of variable bases and differences in sequences of ITS 1 region of Fasciola spp Position Sample
ITS 1 Haplotype
F. hepatica F. gigantica F. gigantica
(H1) (G2) (G1)
Accession Number KF982049 KF982046 KF982047
59
68
186
276
370
448
468
T C C
C T C
C T T
A T T
C T T
T A A
C T T
PCR products were electrophoresed on 1.5% (for ITS 1 PCR) agarose gel in TAE buffer at 60 V for 40 min and visualized by UV illumination (UVITEC) after ethidium bromide staining. To estimate the size of the amplicons, a 100bp DNA ladder (Fermentas) was used in gels.
Sequencing and analysis
To select a suitable restriction enzyme we sequenced PCR products of one Fasciola samples from sheep and one Fasciola samples from cattle. The sequences were analyzed and aligned using Blast (http://blast.ncbi.nlm.nih.gov/Blast) and multalin (http://multalin.toulose.inra.fr/ multalin ) software in comparison with sequences of Fasciola ITS 1, previously published from other countries in GenBank. Nucleotides of 700bp from fasciolid species were subjected to comparison of restriction sites and selection the appropriate enzyme by Webcutter 2 using the following website: (http://bio.lundberg.gu.se/cutter2) software.
Restriction Fragment Length Polymorphism (RFLP)
Rsa I enzyme (Fermentas) was selected for RFLP method that caused the separation specifically of Fasciola species. To performance RFLP, 5 µl of Fasciola ITS 1 PCR product, 2.5 µl of supplied restriction enzyme buffer, 5 µl of restriction enzyme diluted, and D.W up to 25 µl were provided. According to the manufacturer instruction, the tubes were incubated at 37 °C for 7 h, to ensure full cutting of fragments. For analyzing the digestion products, 15 μl of each product in addition to 2 μl of loading buffer were electrophoresed on either 3% agarose gel in TAE buffer at 60 V for 40 min and visualized by UV illumination (UVITEC) af1366
ter ethidium bromide staining. To estimate the size of the separation specifically of Fasciola species, a 100bp DNA ladder (Fermentas) was used in gels.
DNA Sequencing and Phylogenic Analysis
PCR products of ITS 1, of 20 isolates from two hosts (cattle and sheep) were sequenced .DNA Analyzers Sequencing, using the same primers which were used in the PCR.The sequences were aligned and compared with those of existing sequences from the region, related to Fasciola spp. available in the GenBank, using the BLASTprogram of NCBI GenBank. Multiple alignments were performed with data related to Fasciola spp. from Iran and other countries deposited in GenBank. A maximum likelihood tree was constructed using Blast (http://blast.ncbi.nlm.nih.gov/Blast).
Results Genomic DNA was extracted from 70 isolates, which could amplify a fragment approximately 700bp in all of the samples (Fig. 2). Negative control did not produce any band on the gels.All amplified products of Fasciola were digested with the Rsa I restriction enzyme. The band patterns of the fragments digested with Rsa I were accurately distinguished among the two forms of Fasciola.
PCR-RFLP Analysis
Based on RFLP pattern, from 70 Fasciola isolates, 4 isolates (11.4%) from cattle had a RFLP pattern corresponding to F. hepatica, 31 isolates (88.6%) from cattle had a RFLP pattern corresponding to F. gigantica.35 isolates (100%) from sheep had a RFLP pattern corresponding to F. hepatica.
Available at:
http://ijph.tums.ac.ir
Aryaeipour et al.: Genotyping and Phylogenetic Analysis of Fasciola Spp. …
Fig. 2: PCR pattern of Fasciola. Lane M: 100bp DNA ladder, Lanes 6and 17: Fasciola from cattle, Lane 13 and Lane 16: Fasciola from sheep, Lanes N: Control negative
RFLP pattern from F. hepatica, which had 5 cutting sites, were predicted to be separated into fragments of 367, 104, 68,59, 54, and 28 bp in the amplicons. while 5 fragments were produced by 4 cutting sites from F. gigantica including 367, 172, 59, 54, and 28 bp in the amplicons (Fig. 3, 4). The band of 60 bp was thought to include the three fragments of 68 bp (only F. hepatica), 59 and 54 bp, and the small fragment of 28 bp was not detected on agarose gels. These findings showed that the band patterns between the amplicons with the sequences of F. hepatica and F. gigantica were evidently different.
Fig. 3: PCR-RFLP pattern of Fasciola after digestion with Rsa I restriction enzyme. Lane M: 100bp DNA ladder, Lanes 6: F. gigantica from cattle, Lane 17 F. hepatica from cattle, Lane 3: F. hepatica from sheep, Lanes16: F. hepatica from sheep after and before digestion with Rsa I restriction enzyme
Available at: http://ijph.tums.ac.ir
Fig. 4: PCR-RFLP pattern of Fasciola after digestion with Rsa I restriction enzyme. Lane M: 100bp DNA ladder, Lanes 7 F. gigantic, from cattle, Lane 9: PCR amplicons (upper) of F. gigantic and their restriction fragments (lower) produced by the PCR–RFLP method.from cattle, Lane 11 and 12: F. hepatica from sheep, Lane 14: F. gigantica from cattle before digestion with Rsa I restriction enzyme
The sequences of 700 bp ITS 1 of the flukes were analyzed and aligned with those of available sequences in GenBank using multalin (http://multalin.toulose.inra.fr/ multalin) software in comparison (Fig. 5, 6). Alignment of the sequences of ITS 1 showed DNA variable sites in which nucleotides at the position of 59, 68, 186, 276, 370,448, and 468 were single-base substituted resulting in segregation of the specimens (Table 1). Some sequences of F. hepatica and F. gigantica were analyzed and deposited in GenBank (Table 2). Ninety-hundred percent similarities were obtained in comparison of these sequences with all available data of Fasciola spp. in GenBank. Restriction sites of the nucleotides in the Fasciola species were studied by computer software for selection a suitable enzyme. Accordingly, Rsa I was selected as one of the best restriction enzyme for differentiation between F. hepatica and F. gigantica which showed effect on GTAC nucleotides.
Phylogenic tree
Phylogenic Tree of ITS 1 showed that flukes were scattered as pure F. hepatica and F. gigantica clades, suggesting that two genotypes of Fasciola are able to infect animals and probably human in Northwestern of Iran (Fig. 7). 1367
Iranian J Publ Health, Vol. 43, No.10, Oct 2014, pp. 1364-1371
Fig. 5: Sequence alignment of the ITSI, region from Fasciola hepatica (KF982049) and Fasciola gigantic (KF982046 & KF982047) of Ardabil Province, Iran
Fig. 6: Alignment of ITSI sequence of F. hepatica (GQ925431.1& AB514861.1) and F. gigantica (JF295001.1) deposited in GenBank with F. hepatica (KF982049) and F. gigantica (KF982046& KF982047) of Ardabil Province, Iran
1368
Available at:
http://ijph.tums.ac.ir
Aryaeipour et al.: Genotyping and Phylogenetic Analysis of Fasciola Spp. …
Table 2: Profile of Fasciola Spp. used in this study Specimen code ADS1 ADS8 ADS11 ADS16 ADS25 ADS18 ADC45 ADC47 ADC51 ADC53 ADC57
species F.hepatica F.hepatica F.hepatica F.hepatica F.hepatica F.hepatica F.gigantica F.hepatica F.gigantica F.gigantica F.gigantica
Haplotyp H1 H1 H1 H1 H1 H1 G1 H1 G2 G1 G1
Host Sheep Sheep Sheep Sheep Sheep Sheep Cattle Cattle Cattle Cattle Cattle
Sequence analysis F.hepatica F.hepatica F.hepatica F.hepatica F.hepatica F.hepatica F.gigantica F.hepatica F.gigantica F.gigantica F.gigantica
ITS 1 types Accession no. KF982049 KF982045 KJ728732 KJ728734 KJ728735 KJ728736 KF982047 KF982048 KF982046 KJ728737 KJ728738
PCR-RFLP Analysis F.hepatica F.hepatica F.hepatica F.hepatica F.hepatica F.hepatica F.gigantica F.hepatica F.gigantica F.gigantica F.gigantica
Fig. 7: Phylogenetic relationship of ITS 1 sequences of isolates of Fasciola hepatica and F. gigantica from Iran using Blast (http://blast.ncbi.nlm.nih.gov/Blast). Fasciolides magna (AN: EF051080) was used as the out group
Available at: http://ijph.tums.ac.ir
1369
Iranian J Publ Health, Vol. 43, No.10, Oct 2014, pp. 1364-1371
Discussion In this study, we used 700bp region contained ITS1, 18s rRNA for identification and differentiation of Fasciola isolates by PCR-RFLP method. Moreover, both F. hepatica and F. gigantic samples were digested with Rsa I restriction enzyme, which showed different RFLP patterns. PCRRFLP has been used in different studies in Iran to differentiate Fasciola spp. (8-16-18). To detect the rate of similarity or difference of Fasciola samples, altogether 11 samples of both F. hepatica and F. gigantica were sequenced. Accordingly, it was showed that there was about 1.5% difference among some nucleotides. Karimi et al. showed that using Bfrl enzyme there was similar pattern for both aforementioned parasites, while using Dral different profiles are obvious (8). Rokni et al. using TasI for ITS1 region showed differentiation between both species (19). Saki et al., using AvaII and DraII for differentiation between both species (16). In 2013, Shafiei et al.in PCR-RFLP method was used to specifically distinguish F. hepatica from F. gigantica in ITS1 with Rsa I enzyme and in IST2 with MSPI and KpnI (20). In their study no intermediate forms were seen. In 2009, Ghavami et al. reported a digested pattern of 330, 340 and 241 bp in F. hepatica using BamHI and PagI enzymes on ITS 2 gene (17). As for studies in other countries, Shalabi et al. used RAPD-PCR to differentiate F. hepatica and F. gigantica among imported sheep to Saudi Arabia from Sudan. They reported one intermediate case as well (21). In China, Huang et al. using ITS 2 and HSP9211 plus RcaI enzymes has differentiated both species besides intermediate cases. They recommended the first enzyme as better choice (14). In 2010, Madoka Ichikawa used PCR-RFLP method to specifically distinguish F. hepatica from F. gigantica in ITS1 with Rsa I enzyme (22). Ashrafi et al. used ITS2 nucleotide sequencing and showed two species of F. hepatica and F. gigantica in Guilan, northern Iran (23). Accordingly F. hepatica in that region was similar to Europe F. hepatica and F. gigatinca was similar to Burkina Faso isolate in terms of above mentioned method. The authors believed a kind of intermediate isolate in Guilan (23).
1370
Using phylogenetic tree, we could detect a close relationship of our isolates with those of other countries. As mentioned in tables and figures already, we detected two haplotypes of F. gigantica as G1 and G2. It is worth mentioning that recently a hypothesis has been proposed that this F. gigantica might be considered as separate species (24). Our findings showed that all isolates belonging to F. hepatica were of common H1 haplotype. Figure 6 shows that this haplotype has been reported already from Yasuj and Guilan from Iran, as well as from Egypt, Ireland, China and Australia as well. Besides, haplotype G1 of F. gigantica has been reported from Egypt, Vietnam, Indonesia, China, and Yasuj plus Guilan from Iran. Haplotype G2 is reported from Iran and Egypt as well as Saudi Arabia (All references in Fig. 7). As mentioned already no intermediate isolates were detected in our study. The limitation of our study is that more samples should have been collected to be sure of the deduction as much as possible. The strength might be told that this is for the first time that this study has been conducted in this area in its kind.
Conclusion Comparing F. hepatica and F. gigantica sequences showed 90-100% homology between our samples with those originated from other regions available in GenBank, but comparing F. hepatica and F. gigantica with each other showed a few nucleotide differences.
Ethical considerations Ethical issues (Including plagiarism, Informed Consent, misconduct, data fabrication and/or falsification, double publication and/or submission, redundancy, etc) have been completely observed by the authors.
Acknowledgments The study was supported financially by a grant from Tehran University of Medical Sciences No. 23177. The authors declare that there is no conflict of interests.
Available at:
http://ijph.tums.ac.ir
Aryaeipour et al.: Genotyping and Phylogenetic Analysis of Fasciola Spp. …
References 1. Chang Ei, Choi Ho, Park Y (1991). Subcutaneuse Fascioliasis: A case Report. The Korean J Parasitol, 29(4):403-405. 2. Mas-Coma S, Bargues MD, Esteban JG (1999). Human fascioliasis. In: Dalton JP. Editor. Fascioliasis. Wallingford, CAB International, pp.: 411–434. 3. Rokni MB (2008). The present status of human helminthic diseases in Iran. Ann Trop Med Parasitol, 102(4):283–295. 4. Asmar M, Milani A, Amirkhani A et al.(1991). Seroepidemiological investigation of fascioliasis in northern Iran. Med J Islamic Repub Iran, 5:23–27. 5. World Health Organization (1995). Control of Food-borne Trematode Infections.No.849 Technical Report Series. 6. Mahami-Oskouei M, Dalimi A, ForouzandehMoghadam M, Rokni MB (2011). Molecular Identification and Differentiation of Fasciola Isolates Using PCR- RFLP Method Based on Internal Transcribed Spacer (ITS1, 5.8S rDNA, ITS2). Iran J Parasitol, 6(3):35-42. 7. Mas-Coma S, Bargues MD, Valero MA (2005). Fasciolosis and other plant-borne trematodezoonoses. Int J Parasitol, 35:1255–1278. 8. Karimi A (2008). Genetic diagnosis of Fasciola species based on 18S ribosomal DNA sequences. J Biol Sci, 8:1166–1173. 9. Mas-Coma S, Bargues MD (2008). Human liver flukes: a review. Res Rev Parasitol, 57:145–218. 10. Itagaki T, Tsutsumi KI (1998). Triploid form of Fasciola in Japan: genetic relationships between Fasciola hepatica and Fasciola gigantica determined by ITS-2 sequence of nuclear rDNA. Int J Parasitol, 28:777–781. 11. Esteban JC, Bargues MD, Mas-Coma S (1998). Geographical distribution, diagnosis and treatment of human fascioliasis: a review. Res Rev Parasitol, 58:13–42. 12. Marcilla A, Bargues MD, Mas-Coma S (2002). A PCR-RFLP assay for the distinction between Fasciola hepatica and Fasciola gigantica. Mol Cell Probes, 16:327–333. 13. Prasad PK, Biswal DK, Goswami LM, Chatterjee A (2008). Molecular identification of the Indian liver fluke, Fasciola (Trematoda: Fasciolidae) based
Available at: http://ijph.tums.ac.ir
14.
15.
16. 17. 18.
19.
20.
21.
22. 23.
24.
on the ribosomal internal transcribed spacer regions. Parasitol Res, 103:1247–1255. Huang WY, He B, Wang CR, Zhu XQ (2004). Characterisation of Fasciola species from Mainland China by ITS-2 ribosomal DNA sequence. Vet Parasitol, 120:75–83. Asadian S, Mohebali M, Mahmoudi M et al.(2013). Seroprevalence of Human Fascioliasis in Meshkin-Shahr District, Ardabil Province, Northwestern Iran in 2012. Iran J Parasitol, 8(4):516-521. Saki J, Khademvatan S, Yousefi E (2011). Molecular identification of animal Fasciola isolates in Southwest of Iran. Aust J Basic Appl Sci, 5:1878-1883. Ghavami MB, Rahimi P, Haniloo A, Mosavinasab SN (2009). Genotypic and phenotypic analysis of Fasciola Isolates. Iran J Parasitol, 4:61-670. Hossein M, Mirhendi SH, Brandão J, Mirdashti R, Rosado L (2011). Comparison of enzymatic method rapid yeast plus system with RFLP-PCR for identification of isolated yeast from vulvovaginal candidiasis. Iran J Basic Med Sci, 14: 44350. Rokni MB, Mirhendi H, Behnia M (2010). Molecular Characterization of Fasciola hepatica Isolates by RAPD-PCR and Ribosomal ITS1 Sequencing. Iran Red Cresc Med J, 12; (1): 27-32. Shafiei R, Sarkari B, Moshfe A (2013). A Consistent PCR-RFLP Assay Based on ITS-2 Ribosomal DNA for Differentiation of Fasciola Species. Iran J Basic Med Sci, 16(12):1266-9. Shalaby I, Gherbawy Y, Banaja A (2013). Molecular characterization of Fasciola species isolated from imported sheep in Taif region (Saudi Arabia). Trop Biomed, 30(1):15-26. Ichikawa M, Itagaki T (2010). Discrimination of the ITS1 types of Fasciola spp. based on a PCR– RFLP method. Parasitol Res, 106:757–761. Ashrafi K. A survey on human and animal fascioliasis and genotypic and phenotypic characteristics of fasciolids and their relations-hip with lymnaeid snails in Gilan province, northern Iran (PhD Thesis). School of Public Health, Tehran University of Medical Sciences, Iran; 2004. Walker SM, Prodohl PA, Hoey EM et al. (2012). Substantial genetic divergence between morphologically indistinguishable populations of Fasciola suggests the possibility of cryptic speciation. Int J Parasitol, 42; 13-14: 1193–1199.
1371
Iranian J Publ Health, Vol. 43, No.10, Oct 2014, pp.1364-1371
Genotyping and Phylogenetic Analysis of Fasciola Spp. Isolated from Sheep and Cattle Using PCR-RFLP in Ardabil Province, Northwestern Iran Mojgan ARYAEIPOUR 1, Soheila ROUHANI 2, Mojgan BANDEHPOUR 2, Hadi MIRAHMADI 1, *Bahram KAZEMI 2, *Mohammad Bagher ROKNI 3, 4 1. Dept. of Medical Parasitology and Mycology, International Branch, Shahid Beheshti University of Medical Science, Tehran, Iran 2. Dept. Of Biotechnology, School of Medicine, Shahid Beheshti University of Medical Science, Tehran, Iran 3. Dept. of Medical Parasitology and Mycology, School of Public Health, Tehran University of Medical Sciences, Tehran, Iran 4. Center for Research of Endemic Parasites of Iran (CREPI), Tehran University of Medical Sciences, Tehran, Iran *Corresponding Authors: E-mail: [email protected], [email protected] (Received 17 Jun 2014; accepted 21 Sep 2014)
Abstract Background: The aim of this study was to detect the genotype of Fasciola spp. in Meshkin-Shahr, Ardabil Province, northwestern Iran in different hosts using PCR-RFLP. Methods: The parasite hosts included cattle, and sheep. Overall, 70 adult flukes from livers of slaughtered animals were collected from the abattoirs of aforementioned area. The included 35 samples from infected sheep and 35 samples from 35 infected cattle. PCR–RFLP and sequence analysis of the first nuclear ribosomal internal transcribed spacer (ITS 1) region from Fasciola species were used to conduct the study. Results: The fragment of approximately 700bp in all of the Fasciola samples was amplified. PCR products of ITS 1 were subjected for digestion by restriction enzyme. RsaI restriction enzyme was selected for RFLP method that caused the separation specifically of Fasciola species. Amplicons with the sequences of F. hepatica had a pattern of about 360, 100, and 60 bp band size, whereas F. gigantica worms had a profile of 360, 170, and 60 bp in size, respectively. Results based on PCR-RFLP analysis were confirmed by sequence analysis of representative ITS 1 amplicons. No hybrid forms were detected in the present study. All sheep were infected with F. hepatica but cattle were infected with both species. Conclusion: Both species of Fasciola are present in Ardabil. The method described here can be valuable for identification of Fasciola species in endemic parts for fasciolosis, regions with intermediate species and in that overlapping distribution area. Keywords: Fasciola, Genotyping, Fasciolosis, PCR, Iran
Introduction Fascioliasis, caused by Fasciola spp., is one of the most zoonotic diseases with a global extension and of a significant concern in areas with suitable situations (1).” About 2.4-17 million cases of human fascioliasis are estimated in the world and 180 million exposed at risk of fascioliasis that demonstrates the importance of the disease (2)”. 1364
Two great epidemics of human fascioliasis with about 10,000 people infected in each case occurred in the north of Iran (3-5). According to WHO, fascioliasis is one of the most important diseases and parasitic infections transmitted by food (6). Considering the importance of fascioliasis in Iran, various studies especially in order to Available at:
http://ijph.tums.ac.ir
Aryaeipour et al.: Genotyping and Phylogenetic Analysis of Fasciola Spp. …
identify and genotyping of Fasciola seems indispensable. Morphological features of adult worms and eggs are influenced by the different factors such as host type, age of parasite, fixation of the samples and intensity of the infection (7(.Distribution overlap of F. hepatica and F. gigantica in Asia and simultaneous infection with these two parasites in cattle in northern Iran and bug reporting interface is also introduced in other provinces of Iran (8). It is usually difficult to accurate differentiation between Fasciola species based on phenotype criteria (9-11). Today, molecular methods based on DNA analysis is a great impact in many fields, including identifying systematic parasitology, parasite infection, parasite epidemic, analysis of the genetic structure of parasite populations, such as genetic differences between the sexes, gene expression and drug resistance of and the vaccine.Therefore, it seems that the use of the accurate and reliable method for identification and differentiation of Fasciola species is necessary. For this purpose, PCR-RFLP based on 28s rRNA, 18s rRNA, ITS 1 or ITS 2 has been used (12-13). In the present study, region between 18S and 28S (ITS 1, 5.8s,) of ribosomal DNA was used by PCR-RFLP method. Since, ITS 1 sequence was suitable genetic markers for genotyping, interaspecific variations, and phylogenetic studies of parasites (14); we used only ITS 1 for molecular identification of Fasciola species.The adjacent of Ardabil Province to endemic areas of fasciolosis persuaded us to conduct the present study (15).
Materials and Methods Parasite
Overall, 70 Fasciola isolates were studied from two hosts in Ardabil Province, northwestern Iran (Fig. 1). These parasites were isolated from 35 infected sheep and 35 infected cattle. Individual worms were washed extensively in PBS buffer (37°C) and subsequently fixed in 20% ethanol and maintained at room temperature for several weeks until extraction of genomic DNA (Table 1).
Available at: http://ijph.tums.ac.ir
Fig. 1: Map of Iran (a) Ardabil Province, northwestern Iran (b)
DNA extraction
Approximately 10 mg tissue was removed from a portion of the lateral zone of adult flukes and crushed. The ethanol in each of these samples was allowed to evaporate for a few minutes and washed in distilled water for 3 times and squashing before the genomic DNA was extracted with (DNGTM-PLUS) kit was used for genomic DNA extraction of Fasciola Parasites. The extracted DNA, split into aliquots, and stored at -20°C until used in the PCR.
Primer designing
We designed forward (5’-ACCGGTGCTGAGAAGACG-3’) and Reverse (5’-CGACGTACGTGCAGTCCA3’) primers using Gene Runner and BLAST
(http://www.ncbi.nlm.nih.gov/blast) software. The primers were synthesized by Bioneer Company (Korea).
PCR
To amplify a 700 bp region of the ITS 1 sequence, PCR was performed using a pair primer. Total volume of reaction was 15µl containing 1.5µl DNA template, 5µl distilled water, 10pmol of each primers (Forward and Reverse), and 7.5µl master mix (amplicon). PCR amplification was performed in Eppendorf Mastercycler Gradient thermocycler . The temperature profile was as follows: an initiation of 95ºC for 5 min, followed by 30 cycles of 94ºC for 30s (denaturation), 60ºC for 30s (annealing), 72ºC for 30s (extension) and a final extension of 72ºC for 5min followed by cooling at 4ºC. 1365
Iranian J Publ Health, Vol. 43, No.10, Oct 2014, pp. 1364-1371
Table 1: Position of variable bases and differences in sequences of ITS 1 region of Fasciola spp Position Sample
ITS 1 Haplotype
F. hepatica F. gigantica F. gigantica
(H1) (G2) (G1)
Accession Number KF982049 KF982046 KF982047
59
68
186
276
370
448
468
T C C
C T C
C T T
A T T
C T T
T A A
C T T
PCR products were electrophoresed on 1.5% (for ITS 1 PCR) agarose gel in TAE buffer at 60 V for 40 min and visualized by UV illumination (UVITEC) after ethidium bromide staining. To estimate the size of the amplicons, a 100bp DNA ladder (Fermentas) was used in gels.
Sequencing and analysis
To select a suitable restriction enzyme we sequenced PCR products of one Fasciola samples from sheep and one Fasciola samples from cattle. The sequences were analyzed and aligned using Blast (http://blast.ncbi.nlm.nih.gov/Blast) and multalin (http://multalin.toulose.inra.fr/ multalin ) software in comparison with sequences of Fasciola ITS 1, previously published from other countries in GenBank. Nucleotides of 700bp from fasciolid species were subjected to comparison of restriction sites and selection the appropriate enzyme by Webcutter 2 using the following website: (http://bio.lundberg.gu.se/cutter2) software.
Restriction Fragment Length Polymorphism (RFLP)
Rsa I enzyme (Fermentas) was selected for RFLP method that caused the separation specifically of Fasciola species. To performance RFLP, 5 µl of Fasciola ITS 1 PCR product, 2.5 µl of supplied restriction enzyme buffer, 5 µl of restriction enzyme diluted, and D.W up to 25 µl were provided. According to the manufacturer instruction, the tubes were incubated at 37 °C for 7 h, to ensure full cutting of fragments. For analyzing the digestion products, 15 μl of each product in addition to 2 μl of loading buffer were electrophoresed on either 3% agarose gel in TAE buffer at 60 V for 40 min and visualized by UV illumination (UVITEC) af1366
ter ethidium bromide staining. To estimate the size of the separation specifically of Fasciola species, a 100bp DNA ladder (Fermentas) was used in gels.
DNA Sequencing and Phylogenic Analysis
PCR products of ITS 1, of 20 isolates from two hosts (cattle and sheep) were sequenced .DNA Analyzers Sequencing, using the same primers which were used in the PCR.The sequences were aligned and compared with those of existing sequences from the region, related to Fasciola spp. available in the GenBank, using the BLASTprogram of NCBI GenBank. Multiple alignments were performed with data related to Fasciola spp. from Iran and other countries deposited in GenBank. A maximum likelihood tree was constructed using Blast (http://blast.ncbi.nlm.nih.gov/Blast).
Results Genomic DNA was extracted from 70 isolates, which could amplify a fragment approximately 700bp in all of the samples (Fig. 2). Negative control did not produce any band on the gels.All amplified products of Fasciola were digested with the Rsa I restriction enzyme. The band patterns of the fragments digested with Rsa I were accurately distinguished among the two forms of Fasciola.
PCR-RFLP Analysis
Based on RFLP pattern, from 70 Fasciola isolates, 4 isolates (11.4%) from cattle had a RFLP pattern corresponding to F. hepatica, 31 isolates (88.6%) from cattle had a RFLP pattern corresponding to F. gigantica.35 isolates (100%) from sheep had a RFLP pattern corresponding to F. hepatica.
Available at:
http://ijph.tums.ac.ir
Aryaeipour et al.: Genotyping and Phylogenetic Analysis of Fasciola Spp. …
Fig. 2: PCR pattern of Fasciola. Lane M: 100bp DNA ladder, Lanes 6and 17: Fasciola from cattle, Lane 13 and Lane 16: Fasciola from sheep, Lanes N: Control negative
RFLP pattern from F. hepatica, which had 5 cutting sites, were predicted to be separated into fragments of 367, 104, 68,59, 54, and 28 bp in the amplicons. while 5 fragments were produced by 4 cutting sites from F. gigantica including 367, 172, 59, 54, and 28 bp in the amplicons (Fig. 3, 4). The band of 60 bp was thought to include the three fragments of 68 bp (only F. hepatica), 59 and 54 bp, and the small fragment of 28 bp was not detected on agarose gels. These findings showed that the band patterns between the amplicons with the sequences of F. hepatica and F. gigantica were evidently different.
Fig. 3: PCR-RFLP pattern of Fasciola after digestion with Rsa I restriction enzyme. Lane M: 100bp DNA ladder, Lanes 6: F. gigantica from cattle, Lane 17 F. hepatica from cattle, Lane 3: F. hepatica from sheep, Lanes16: F. hepatica from sheep after and before digestion with Rsa I restriction enzyme
Available at: http://ijph.tums.ac.ir
Fig. 4: PCR-RFLP pattern of Fasciola after digestion with Rsa I restriction enzyme. Lane M: 100bp DNA ladder, Lanes 7 F. gigantic, from cattle, Lane 9: PCR amplicons (upper) of F. gigantic and their restriction fragments (lower) produced by the PCR–RFLP method.from cattle, Lane 11 and 12: F. hepatica from sheep, Lane 14: F. gigantica from cattle before digestion with Rsa I restriction enzyme
The sequences of 700 bp ITS 1 of the flukes were analyzed and aligned with those of available sequences in GenBank using multalin (http://multalin.toulose.inra.fr/ multalin) software in comparison (Fig. 5, 6). Alignment of the sequences of ITS 1 showed DNA variable sites in which nucleotides at the position of 59, 68, 186, 276, 370,448, and 468 were single-base substituted resulting in segregation of the specimens (Table 1). Some sequences of F. hepatica and F. gigantica were analyzed and deposited in GenBank (Table 2). Ninety-hundred percent similarities were obtained in comparison of these sequences with all available data of Fasciola spp. in GenBank. Restriction sites of the nucleotides in the Fasciola species were studied by computer software for selection a suitable enzyme. Accordingly, Rsa I was selected as one of the best restriction enzyme for differentiation between F. hepatica and F. gigantica which showed effect on GTAC nucleotides.
Phylogenic tree
Phylogenic Tree of ITS 1 showed that flukes were scattered as pure F. hepatica and F. gigantica clades, suggesting that two genotypes of Fasciola are able to infect animals and probably human in Northwestern of Iran (Fig. 7). 1367
Iranian J Publ Health, Vol. 43, No.10, Oct 2014, pp. 1364-1371
Fig. 5: Sequence alignment of the ITSI, region from Fasciola hepatica (KF982049) and Fasciola gigantic (KF982046 & KF982047) of Ardabil Province, Iran
Fig. 6: Alignment of ITSI sequence of F. hepatica (GQ925431.1& AB514861.1) and F. gigantica (JF295001.1) deposited in GenBank with F. hepatica (KF982049) and F. gigantica (KF982046& KF982047) of Ardabil Province, Iran
1368
Available at:
http://ijph.tums.ac.ir
Aryaeipour et al.: Genotyping and Phylogenetic Analysis of Fasciola Spp. …
Table 2: Profile of Fasciola Spp. used in this study Specimen code ADS1 ADS8 ADS11 ADS16 ADS25 ADS18 ADC45 ADC47 ADC51 ADC53 ADC57
species F.hepatica F.hepatica F.hepatica F.hepatica F.hepatica F.hepatica F.gigantica F.hepatica F.gigantica F.gigantica F.gigantica
Haplotyp H1 H1 H1 H1 H1 H1 G1 H1 G2 G1 G1
Host Sheep Sheep Sheep Sheep Sheep Sheep Cattle Cattle Cattle Cattle Cattle
Sequence analysis F.hepatica F.hepatica F.hepatica F.hepatica F.hepatica F.hepatica F.gigantica F.hepatica F.gigantica F.gigantica F.gigantica
ITS 1 types Accession no. KF982049 KF982045 KJ728732 KJ728734 KJ728735 KJ728736 KF982047 KF982048 KF982046 KJ728737 KJ728738
PCR-RFLP Analysis F.hepatica F.hepatica F.hepatica F.hepatica F.hepatica F.hepatica F.gigantica F.hepatica F.gigantica F.gigantica F.gigantica
Fig. 7: Phylogenetic relationship of ITS 1 sequences of isolates of Fasciola hepatica and F. gigantica from Iran using Blast (http://blast.ncbi.nlm.nih.gov/Blast). Fasciolides magna (AN: EF051080) was used as the out group
Available at: http://ijph.tums.ac.ir
1369
Iranian J Publ Health, Vol. 43, No.10, Oct 2014, pp. 1364-1371
Discussion In this study, we used 700bp region contained ITS1, 18s rRNA for identification and differentiation of Fasciola isolates by PCR-RFLP method. Moreover, both F. hepatica and F. gigantic samples were digested with Rsa I restriction enzyme, which showed different RFLP patterns. PCRRFLP has been used in different studies in Iran to differentiate Fasciola spp. (8-16-18). To detect the rate of similarity or difference of Fasciola samples, altogether 11 samples of both F. hepatica and F. gigantica were sequenced. Accordingly, it was showed that there was about 1.5% difference among some nucleotides. Karimi et al. showed that using Bfrl enzyme there was similar pattern for both aforementioned parasites, while using Dral different profiles are obvious (8). Rokni et al. using TasI for ITS1 region showed differentiation between both species (19). Saki et al., using AvaII and DraII for differentiation between both species (16). In 2013, Shafiei et al.in PCR-RFLP method was used to specifically distinguish F. hepatica from F. gigantica in ITS1 with Rsa I enzyme and in IST2 with MSPI and KpnI (20). In their study no intermediate forms were seen. In 2009, Ghavami et al. reported a digested pattern of 330, 340 and 241 bp in F. hepatica using BamHI and PagI enzymes on ITS 2 gene (17). As for studies in other countries, Shalabi et al. used RAPD-PCR to differentiate F. hepatica and F. gigantica among imported sheep to Saudi Arabia from Sudan. They reported one intermediate case as well (21). In China, Huang et al. using ITS 2 and HSP9211 plus RcaI enzymes has differentiated both species besides intermediate cases. They recommended the first enzyme as better choice (14). In 2010, Madoka Ichikawa used PCR-RFLP method to specifically distinguish F. hepatica from F. gigantica in ITS1 with Rsa I enzyme (22). Ashrafi et al. used ITS2 nucleotide sequencing and showed two species of F. hepatica and F. gigantica in Guilan, northern Iran (23). Accordingly F. hepatica in that region was similar to Europe F. hepatica and F. gigatinca was similar to Burkina Faso isolate in terms of above mentioned method. The authors believed a kind of intermediate isolate in Guilan (23).
1370
Using phylogenetic tree, we could detect a close relationship of our isolates with those of other countries. As mentioned in tables and figures already, we detected two haplotypes of F. gigantica as G1 and G2. It is worth mentioning that recently a hypothesis has been proposed that this F. gigantica might be considered as separate species (24). Our findings showed that all isolates belonging to F. hepatica were of common H1 haplotype. Figure 6 shows that this haplotype has been reported already from Yasuj and Guilan from Iran, as well as from Egypt, Ireland, China and Australia as well. Besides, haplotype G1 of F. gigantica has been reported from Egypt, Vietnam, Indonesia, China, and Yasuj plus Guilan from Iran. Haplotype G2 is reported from Iran and Egypt as well as Saudi Arabia (All references in Fig. 7). As mentioned already no intermediate isolates were detected in our study. The limitation of our study is that more samples should have been collected to be sure of the deduction as much as possible. The strength might be told that this is for the first time that this study has been conducted in this area in its kind.
Conclusion Comparing F. hepatica and F. gigantica sequences showed 90-100% homology between our samples with those originated from other regions available in GenBank, but comparing F. hepatica and F. gigantica with each other showed a few nucleotide differences.
Ethical considerations Ethical issues (Including plagiarism, Informed Consent, misconduct, data fabrication and/or falsification, double publication and/or submission, redundancy, etc) have been completely observed by the authors.
Acknowledgments The study was supported financially by a grant from Tehran University of Medical Sciences No. 23177. The authors declare that there is no conflict of interests.
Available at:
http://ijph.tums.ac.ir
Aryaeipour et al.: Genotyping and Phylogenetic Analysis of Fasciola Spp. …
References 1. Chang Ei, Choi Ho, Park Y (1991). Subcutaneuse Fascioliasis: A case Report. The Korean J Parasitol, 29(4):403-405. 2. Mas-Coma S, Bargues MD, Esteban JG (1999). Human fascioliasis. In: Dalton JP. Editor. Fascioliasis. Wallingford, CAB International, pp.: 411–434. 3. Rokni MB (2008). The present status of human helminthic diseases in Iran. Ann Trop Med Parasitol, 102(4):283–295. 4. Asmar M, Milani A, Amirkhani A et al.(1991). Seroepidemiological investigation of fascioliasis in northern Iran. Med J Islamic Repub Iran, 5:23–27. 5. World Health Organization (1995). Control of Food-borne Trematode Infections.No.849 Technical Report Series. 6. Mahami-Oskouei M, Dalimi A, ForouzandehMoghadam M, Rokni MB (2011). Molecular Identification and Differentiation of Fasciola Isolates Using PCR- RFLP Method Based on Internal Transcribed Spacer (ITS1, 5.8S rDNA, ITS2). Iran J Parasitol, 6(3):35-42. 7. Mas-Coma S, Bargues MD, Valero MA (2005). Fasciolosis and other plant-borne trematodezoonoses. Int J Parasitol, 35:1255–1278. 8. Karimi A (2008). Genetic diagnosis of Fasciola species based on 18S ribosomal DNA sequences. J Biol Sci, 8:1166–1173. 9. Mas-Coma S, Bargues MD (2008). Human liver flukes: a review. Res Rev Parasitol, 57:145–218. 10. Itagaki T, Tsutsumi KI (1998). Triploid form of Fasciola in Japan: genetic relationships between Fasciola hepatica and Fasciola gigantica determined by ITS-2 sequence of nuclear rDNA. Int J Parasitol, 28:777–781. 11. Esteban JC, Bargues MD, Mas-Coma S (1998). Geographical distribution, diagnosis and treatment of human fascioliasis: a review. Res Rev Parasitol, 58:13–42. 12. Marcilla A, Bargues MD, Mas-Coma S (2002). A PCR-RFLP assay for the distinction between Fasciola hepatica and Fasciola gigantica. Mol Cell Probes, 16:327–333. 13. Prasad PK, Biswal DK, Goswami LM, Chatterjee A (2008). Molecular identification of the Indian liver fluke, Fasciola (Trematoda: Fasciolidae) based
Available at: http://ijph.tums.ac.ir
14.
15.
16. 17. 18.
19.
20.
21.
22. 23.
24.
on the ribosomal internal transcribed spacer regions. Parasitol Res, 103:1247–1255. Huang WY, He B, Wang CR, Zhu XQ (2004). Characterisation of Fasciola species from Mainland China by ITS-2 ribosomal DNA sequence. Vet Parasitol, 120:75–83. Asadian S, Mohebali M, Mahmoudi M et al.(2013). Seroprevalence of Human Fascioliasis in Meshkin-Shahr District, Ardabil Province, Northwestern Iran in 2012. Iran J Parasitol, 8(4):516-521. Saki J, Khademvatan S, Yousefi E (2011). Molecular identification of animal Fasciola isolates in Southwest of Iran. Aust J Basic Appl Sci, 5:1878-1883. Ghavami MB, Rahimi P, Haniloo A, Mosavinasab SN (2009). Genotypic and phenotypic analysis of Fasciola Isolates. Iran J Parasitol, 4:61-670. Hossein M, Mirhendi SH, Brandão J, Mirdashti R, Rosado L (2011). Comparison of enzymatic method rapid yeast plus system with RFLP-PCR for identification of isolated yeast from vulvovaginal candidiasis. Iran J Basic Med Sci, 14: 44350. Rokni MB, Mirhendi H, Behnia M (2010). Molecular Characterization of Fasciola hepatica Isolates by RAPD-PCR and Ribosomal ITS1 Sequencing. Iran Red Cresc Med J, 12; (1): 27-32. Shafiei R, Sarkari B, Moshfe A (2013). A Consistent PCR-RFLP Assay Based on ITS-2 Ribosomal DNA for Differentiation of Fasciola Species. Iran J Basic Med Sci, 16(12):1266-9. Shalaby I, Gherbawy Y, Banaja A (2013). Molecular characterization of Fasciola species isolated from imported sheep in Taif region (Saudi Arabia). Trop Biomed, 30(1):15-26. Ichikawa M, Itagaki T (2010). Discrimination of the ITS1 types of Fasciola spp. based on a PCR– RFLP method. Parasitol Res, 106:757–761. Ashrafi K. A survey on human and animal fascioliasis and genotypic and phenotypic characteristics of fasciolids and their relations-hip with lymnaeid snails in Gilan province, northern Iran (PhD Thesis). School of Public Health, Tehran University of Medical Sciences, Iran; 2004. Walker SM, Prodohl PA, Hoey EM et al. (2012). Substantial genetic divergence between morphologically indistinguishable populations of Fasciola suggests the possibility of cryptic speciation. Int J Parasitol, 42; 13-14: 1193–1199.
1371
