J Arthropod-Borne Dis, December 2020, 14(4): 416–424 F Talebzadeh et al.: Molecular Species Identification of … 416 http://jad.tums.ac.ir Published Online: December 31, 2020 Original Article Molecular Species Identification of Six Forensically Important Iranian Flesh Flies (Diptera) Fahimeh Talebzadeh1; *Mohammad Ali Oshaghi1; *Kamran Akbarzadeh1; Somayeh Panahi- Moghadam2 1Department of Medical Entomology and Vector Control, School of Public Health, Tehran University of Medical Sciences, Tehran, Iran 2Department of Cell and Molecular Biology, Faculty of Life Sciences and Biotechnology, Shahid Beheshti University, Tehran, Iran *Corresponding authors: Dr Mohammad Ali Oshaghi, E-mail: moshaghi@tums.ac.ir, Dr Kamran Akbarzadeh, E-mail: kakbarzadeh@tums.ac.ir (Received 21 June 2020; accepted 22 Dec 2020) Abstract Background: Flesh flies (Diptera: Sarcophagidae) are considered as myiasis agents and important evidences in forensic investigations. However, their use has been restricted because, at all larval stages and female adults, morphological species identification is difficult or very challenging. This study investigated to test utility of mitochondrial cytochrome oxidase subunit I (mt-COI) sequences for differentiation of six forensically important Iranian flesh flies namely, Sar- cophaga crassipalpis, S. flagellifera, S. hirtipes, S. aegyptica, S. africa and S. argyrostoma. Methods: Male specimens were morphologically identified to species level and then the genomic DNA of the flies were extracted and subjected to polymerase chain reaction (PCR) against mt-COI gene. The PCR products were sequenced and the obtained sequences were analyzed for the species specific restriction fragment length polymorphisms (RFLPs). Results: Rate of genetic variation between species was 6–10% which was enough to find restriction enzymes (RE) that were able to produce species-specific RFLP profiles. Combinations of three REs: BsrFI, RsaI and HinfI, provided diag- nostic bands for identification of the six Sarcophaga species. Conclusion: The results of this study showed that molecular markers such as RFLPs enhancing the use of evidence from flesh flies in forensic investigation. However, lack proper restriction sites in the COI region inhibited introduction of a single restriction enzyme for easy species identification. It is recommended to apply larger part of DNA such as combination of COI and COII genes to provide better RFLP markers for species identification of flesh flies. Keywords: Sarcophagidae; Forensic entomology; Flesh fly, Mitochondrial cytochrome oxidase subunit I (mt-COI) Introduction A range of insects present on a corpse can be used as evidence in forensic investigations to estimate the post-mortem interval (PMI). Since rate of development in insects is spe- cies-specific, commonly, estimation of PMI using insect evidence requires accurate spe- cies identification, with subsequent examina- tion of thermos biological profiles to deter- mine insect age (1-4). When larvae or pupae in various stages of development are collected from the site of investigation and the growth rates of samples are known, an approximate time of oviposition or larviposition can be es- timated (5). For accuracy, forensic entomolo- gists preferentially use evidence from initial corpse colonizers, such as flesh flies (Diptera: Sarcophagidae) and blow flies (Diptera: Cal- liphoridae) (2). It can be argued that sarcoph- agids have the potential to provide a more ac- curate PMI estimate than calliphorids, based on differences in their life cycles (6). Sarcoph- agids lay live larvae (viviparous), providing developed immatures ready to start feeding immediately. However, most calliphorids lay eggs (oviparous), which will only hatch into larvae when the correct environmental condi- Copyright © 2020 The Authors. Published by Tehran University of Medical Sciences. This work is licensed under a Creative Commons Attribution-NonCommercial 4.0 International license (https://creativecommons.org/licenses/by- nc/4.0/). Non-commercial uses of the work are permitted, provided the original work is properly cited. http://jad.tums.ac.ir/ mailto:moshaghi@tums.ac.ir mailto:kakbarzadeh@tums.ac.ir https://creativecommons.org/licenses/by-nc/4.0/ https://creativecommons.org/licenses/by-nc/4.0/ J Arthropod-Borne Dis, December 2020, 14(4): 416–424 F Talebzadeh et al.: Molecular Species Identification of … 417 http://jad.tums.ac.ir Published Online: December 31, 2020 tions are met (3). Despite the prospective use of sarcophagids in forensic investigations, their use to date has been overshadowed by calli- phorids. This is due to the difficulties of mor- phological species-level identification at any life stage of flesh flies, and a lack of documented thermos-biological profiles of these insects. Sarcophagid flies comprise more than 2,500 species in more than 100 genera, of which approximately 800 species belong to the ge- nus Sarcophaga Meigen (7), widespread, and often among the first insects to colonize a corpse (6). The genus is notorious for having many species with a highly similar morpho- logical appearance, females are difficult to iden- tify, and immature stages are often unknown or lack species-specific diagnostic characters (6, 8). Adult sarcophagids can be easily identified at the family level, as most species share the characteristic features of longitudinal stripes on the thorax and a tessellated/chequered abdominal pattern. Identification of the species of the ge- nus Sarcophaga is traditionally accomplished on males using characters of the phallus (9, 10). However, species-level identification is difficult and requires examination of subtle morphological variation of bristle placement and length, hair coloration, body pigmentation and genitalic structure of adults (7). This is un- fortunate because the PMI is sometimes esti- mated using larvae and pupae. Rearing larvae and pupae to the adult level may allow iden- tification but is time consuming, often not pos- sible for all species, and may result in a num- ber of unidentifiable females. Considering this, molecular-based approaches for species iden- tifications have been proposed to eliminate is- sues with identifications based exclusively on taxonomy (11, 12). Species identification is essential for de- termining growth rates, as these rates are spe- cies-specific (2). Therefore, species identifica- tion is a key step in estimating the PMI from entomological evidence. The traditional spe- cies identification method is dependent on the morphological features of insects and is not easily applicable to immature samples such as eggs, larvae, and pupae (4-9). Moreover, only a few expert taxonomists specialize in foren- sically important insect species. DNA-based approaches have been devel- oped in an effort to improve accessibility to methods of species identification. The mito- chondrial encoded cytochrome c oxidase sub- unit I (COI) gene has been shown to be a ma- jor candidate gene to identify forensically im- portant insects (12-14). Numerous studies have evaluated the effectiveness of COI gene, with the approach shown to be unreliable for some Diptera (15, 16) but also proven successful for many groups of invertebrates, such as mos- quitoes, cockroaches, and black flies (17-22). Although the COI gene has been used widely for species identification of forensic insects (23-29), however, only a few previously report- ed studies have used the COI gene for a few fo- rensically important species in Iran (30). There- fore, there has been little effort to characterize the COI haplotypes of Iranian Sarcophagidae fly species. This study examined about 711 nu- cleotide sequences of the 5’ end of COI gene of six Sarcophagidae fly species from the ge- nus Sarcophaga collected in Iran. It is hoped that the results from this study will assist with the implementation of Iranian Sarcophagidae in forensic investigations. Materials and Methods Specimen collection Trapping at decayed meat baits (sheep’s or chicken’s liver and fish carcasses) in plas- tic bottle fly traps, and hand netting were em- ployed to collect adult sarcophagid specimens in Tehran provinces at the Laleh Park in cen- ter of Tehran as well as the livestock shopping center and close vicinity of slaughter houses in the east and south of Tehran, Iran. The cap- tured flies in entomological nets immediately transferred into glass jars. The jars or baits har- boring flies were then transported in a poly- styrene icebox to the laboratory of Medical http://jad.tums.ac.ir/ J Arthropod-Borne Dis, December 2020, 14(4): 416–424 F Talebzadeh et al.: Molecular Species Identification of … 418 http://jad.tums.ac.ir Published Online: December 31, 2020 Entomology, School of Public Health, Tehran University of Medical Sciences (SPH-TUMS). Live adult were transferred individually into a bottle trap including sand and meat to lay eggs at 28 °C±1, 40%±5 relative humidity and 12h photoperiodicity, protected with an external net curtain to avoid the entry of other insect spe- cies. After laying eggs, the dead specimens were identified morphologically by using the taxo- nomic keys (7, 30, 31). To identify the speci- mens to species level, the male specimens were dissected and the genitalia of each male spec- imen were examined. All specimens then were kept in 70% ethanol and stored at -20 °C for further molecular investigation in the labora- tory of SPH-TUMS. DNA Extraction, Amplification and Se- quencing One to three legs from each known adult male sarcophagid specimen were used as tis- sue for total genomic DNA extractions using the Qiagen DNA extraction kit, following the protocols recommended by the manufacturer (Qiagen, Hilden, Germany). The DNA was re- suspended in 50–100µl of fresh TE solution (1 mM Tris–HCl (pH 8), 0.1mM EDTA) and subsequently stored at 4° C. The COI region of approximately 711bp was amplified and se- quenced using the primers of COI-F 5’-GG TCAACAAATCATAAAGATATTGG-3’ and COI-R 5’-TAAACTTCAGGGTGACCAAA AAATCA-3’ and thermal cycling conditions ex- plained previously (18, 19, 33). The PCR prod- ucts were visualized on a UV transilluminator following electrophoresis on a 1.2% agarose gel containing ethidium bromide. The PCR prod- ucts were purified from seized gels according to the supplier's guidelines (Amersham Phar- macia Biotech, USA) and used directly for DNA sequencing in an automatic sequencer (Seqlab, Goä ttingen, Germany). The sequences were generated in both directions, and the resultant chromatograms were edited using ChromasPro Version 1.33 (available online at www.tech- nelysium.com.au/ChromasPro.html) and aligned using ClustalW2 (34). Homology of the se- quences with GenBank sequences was assessed using FASTA search (https://www.ebi.ac.uk/Tools/sss/). Consensus sequences for the COI region was deposited in the GenBank database. Multiple alignments of the nucleotide sequences were performed using ClustalW2 program (http://www.ebi.ac.uk/Tools/msa/clustalw2). Following sequencing, polymerase chain reaction-restriction fragment length polymor- phism (PCR-RFLP) assay was used to distin- guish the six species. Based on sequence vari- ation among the six species, three restriction enzymes were selected by manufacturer’s NEBcutter V2.0 program (http://nc2.neb.com/NEBcutter2/) and a phys- ical map provided (Fig. 1). For restriction fragment assays, the 711bp COI PCR prod- ucts were digested separately with restriction endonuclease BsrFI, RsaI, and HinfI (Vi- vantis) according to the manufacturer’s rec- ommendation. Then, the digested sample was resolved by electrophoresis in 2.5% TBE aga- rose gel, stained with ethidium bromide and viewed under UV-illumination. The restriction patterns were photographed and analyzed. Results Morphological Species Identifications In the present study, 273 specimens were collected from Tehran district and 86 speci- mens comprising six species of Sarcophagidae were identified using the taxonomic literature. These specimens collectively represent six iden- tifiable taxa: S. argyrostoma Robineau-Desvoidy 1830, S. africa Wiedemann 1824, S. aegypti- ca/dux Salem 1935/Thomson 1869, S. crassi- palpis Macquart 1839, S. flagellifera Grunin 1964, and S. hirtipes Wiedemann 1830. These identifications were based on male morpho- logical characters, as these are the only spec- imens that can be reliably identified. Difficul- ties were encountered in accurately identifying some of the 273 specimens using the available taxonomic keys. To further assist with identi- fications, each ‘unknown’ sequence was com- http://jad.tums.ac.ir/ http://www.technelysium.com.au/ChromasPro.html http://www.technelysium.com.au/ChromasPro.html http://www.ebi.ac.uk/Tools/msa/clustalw2 J Arthropod-Borne Dis, December 2020, 14(4): 416–424 F Talebzadeh et al.: Molecular Species Identification of … 419 http://jad.tums.ac.ir Published Online: December 31, 2020 pared with the available data in Genbank, how- ever no conclusive matches were obtained. We cannot confidently associate the unknown Sarcophaga species with a particular subgenus, and it is possible that these could represent new species, given that no broad work on the Iranian Sarcophaga fauna has been carried out in the past. It is notable that there are 68 re- ported Sarcophagidae species in Iran (35). Mitochondrial cytochrome oxidase subunit I (COI) Analysis The Qiagen DNeasy Blood and Tissue kit (Qiagen) was used to extract genomic DNA from one to three legs of each specimen in each extraction to retain the integrity of the speci- mens. Known specific mitochondrion COI pri- mers were used to minimize the possibility of amplifying nuclear pseudogenes. About 711 bp products of 5’ end of COI gene were able to be amplified using the primer combination of COI-F and COI-R from the DNA extracted from the specimen legs. Direct PCR-sequencing were used and only strongly amplified products were sequenced. Further evidence that the ob- tained sequences were of mitochondrial origin came from BLAST search followed by the ob- servation that they did not contain base ambi- guities or premature stop codons upon transla- tion. This mtDNA region was observed to have a strong AT bias (average 67%), which is typ- ical of insect mitochondrial DNA (36). The ambiguous parts from the beginning and end of sequences for the six species were removed and finally 621bp reliable data were analyzed. Sequence comparison among species showed a maximum of 0.2% intra-specific variation for the species whereas the interspecific varia- tion ranged 6.23–10%. The most divergence rate (10%) was observed between S. argyro- stoma and S. flagellifera and the least rate (6.23%) was seen between S. hirtipes and S. aegyptica. Based on nucleotide sequence analysis of the mtDNA-COI gene, a PCR-RFLP assay was developed for the diagnosis and differentiation of the six Iranian Sarcophaga species. It was accomplished at first stage by digesting the PCR products (711bp) with BsrFI restriction enzyme. This enzyme cuts the PCR–amplified mtDNA-COI fragment of S. aegyptica into two bands of 511 and 200bp whereas it does not cut the amplicons of other five species and the amplicons remain intact (Fig. 2). At next step, RsaI enzyme was used for the five remaining species. This enzyme produces three profiles, 1) a unique profile for S. hirtipes with 485 and 226bp, 2) a specific profile for S. felagellifera by producing 643 and 68bp bands, and 3) a profile with two bands of 580 and 130bp for three species of S. africa, S. crassipalpis, and S. argyrostoma (Fig. 2). At third stage, HinfI enzyme was used to separate 1) S. africa with three bands of 514, 172, and 25bp, 2) S. argy- rostoma with two bands of 550 and 80bp, and 3) S. crassipalpis with an intact 711bp band (Fig. 2). Fig. 1. Physical map of 711bp of 3’ end of mtDNA-COI for six Iranian Sarcophaga species. B: BsrFI, R: RsaI, H: Hin- fI restriction enzymes http://jad.tums.ac.ir/ J Arthropod-Borne Dis, December 2020, 14(4): 416–424 F Talebzadeh et al.: Molecular Species Identification of … 420 http://jad.tums.ac.ir Published Online: December 31, 2020 Fig. 2. PCR-RFLP assay of 711bp of COI region digested with BsrFI, RsaI, and HinfI. Lane 1 is a 100bp plus ladder (Sinaclon, Iran), 2: Sarcophaga hirtipes, 3: S. africa, 4: S. crassipalpis, 5: S. flagellifera, 6: S. aegyptica, 7: S. argyrostoma Discussion This study was aimed at carefully evaluat- ing the COI gene for species identification of six forensically important Iranian Sarcopha- gidae. Calculation of the percentage diver- gence between sequences is used to quantita- tively evaluate the success of COI DNA. For successful species-level resolution using the COI approach, genetic variation between spe- cies (interspecific) exceeds that of within spe- cies (intraspecific) in which intraspecific vari- ation should be less than 3% and interspecific variation should be more than 3% (37, 38). This was the case for all of the six Iranian Sarcoph- agidae examined in this study, where the per- centage divergences lie between 6.23 to 10 percent through the 621bp compared. Similar result to this was obtained when different part of COI gene were previously evaluated for identification of Sarcophagidae (23, 24, 26, 27). Indeed, various fragments of the COI gene (for example: mini-barcode 127bp, 272–278bp, standard barcode 658bp, 700bp, and entire gene 1,535bp) that show low sequence divergence within species but high divergences among spe- cies. Hence, this gene can be employed as tax- on “barcodes” and unknown samples can be placed accurately into species groups simply by calculating their pairwise genetic distances with reference sequences of a “barcode library” (37, 39). As such, an online COI identification system will provide a cost-effective and widely accessible identification tool for flies of fo- rensic interest, even for immature or damaged specimens. At present, such practice are being in- creased for species of the genus Sarcophaga, since COI sequence database has been devel- oped for 12 species from the south of the Czech Republic (12), for seven species within Canada and the USA (6), for 17 Malaysian, two Indonesian, one Japanese species (40), for 16 Australian species (26), for 56 species in West Europe (25), for several Chinese species (23, 24, 28), for six Korean species (27), ten Iranian species (41) and recently for ten Indi- an species (29). However, a database for other regions in the world is lacking. Since the above mentioned studies have shown that intraspe- cific sequence divergence in Sarcophaga is usually well below 1.5%, whereas interspecif- ic divergence is more than 2–2.5%, it can be indicated that the COI gene has a large poten- tial as an identification tool in this forensical- ly interesting taxon. In this study, we introduce PCR-RFLP markers using 711bp of 3’ end of COI gene to identify six species of the genus Sarcophaga that occur in Iran. This is the first mtDNA COI PCR-RFLP assay for species identifica- tion of Sarcophagidae in the literature. From the PCR-RFLP analysis using combination of three restriction enzymes, each species showed http://jad.tums.ac.ir/ J Arthropod-Borne Dis, December 2020, 14(4): 416–424 F Talebzadeh et al.: Molecular Species Identification of … 421 http://jad.tums.ac.ir Published Online: December 31, 2020 a distinct and unique restriction pattern com- pared to their closely related species. These six species are synanthropic and commonly present in the same collection site, and show high similarities in morphology and behavior. These results suggest that PCR-RFLP assay are useful to distinguish these species. PCR- RFLP technique is a very cheap and rapid meth- od for species identification of many organ- isms. This molecular assay was used to identi- fy many species of different group of organ- isms such as rodents (42), pathogens (43-46), Anopheles spp (17, 18, 42) and to determine blood type within insect guts (47, 48). This will provide a cost-effective solution that so many specimens can be studied with no need to provide DNA sequencing. In this method, using a set of conserved primers, a tiny amount of DNA can be amplified and then be identi- fied by restriction enzyme that overcomes lack of ample amount of flesh fly body parts. Conclusion The COI sequences of each species were unique and distinguishable from each other, although they showed high homology. Spe- cies identification from immature dipteran by the DNA sequences was simple and time-sav- ing because there was no need to wait for adult emergence or knowledge of morphological keys. Nevertheless, given that sarcophagid lar- vae are likely to be found on a human corpse and cannot be identified based on anatomy, and that the identification of females is very difficult and requires strong taxonomic exper- tise, the use of the COI gene and or PCR-RFLP as a species identifier will greatly enhance rap- id and correct species identification. Acknowledgements This research was supported by Tehran University of Medical Sciences (TUMS), Iran, grant No. 29257. The authors declare that there is no conflict of interests. References 1. Amendt J, Campobasso CP, Gaudry E, Reiter C, LeBlanc HN, Hall MJ (2007) Best practice in forensic entomology--stand- ards and guidelines. Int J Legal Med. 121(2): 90–104. 2. Amendt J, Krettek R, Zehner R (2004) Fo- rensic entomology. Die Naturwissen- schaften. 91(2): 51–65. 3. Byrd JH, Castner JL (2001) Forensic En- tomology the Utility of Arthropods in Legal Investigations. CRC Press, New York. 4. Goff ML (1993) Estimation of postmortem interval using arthropod development and successional patterns. Forensic Sci Rev. 5(2): 81–94. 5. Kamal AS (1958) Comparative study of thirteen species of sarcosaprophagous Cal- liphoridae and Sarcophagidae (Diptera) I. Bionomics. Ann Entomol Soc Am. 51 (3): 261–271. 6. Wells JD, Pape T, Sperling FA (2001) DNA- based identification and molecular sys- tematics of forensically important Sar- cophagidae (Diptera). J Forensic Sci. 46 (5): 1098–1102. 7. Pape T (1996) Catalogue of the Sarcopha- gidae of the World (Insecta: Diptera). Associated Publishers, Gainesville, Flor- ida. 8. Saigusa K, Takamiya M, Aoki Y (2005) Spe- cies identification of the forensically im- portant flies in Iwate prefecture, Japan based on mitochondrial cytochrome oxi- dase gene subunit I (COI) sequences. Leg Med (Tokyo). 7(3): 175–178. 9. Dahlem GA DW (1996) Revision of the ge- nus Boetcheria in America north of Mex- ico (Diptera: Sarcophagidae). Insecta Mundi. 10: 77–103. 10. Povolný DVY (1997) The flesh flies of Central Europe. Spixiana Zeit Zool Suppl. http://jad.tums.ac.ir/ J Arthropod-Borne Dis, December 2020, 14(4): 416–424 F Talebzadeh et al.: Molecular Species Identification of … 422 http://jad.tums.ac.ir Published Online: December 31, 2020 24: 1–263. 11. Wells JD, Stevens JR (2008) Application of DNA-based methods in forensic en- tomology. Annu Rev Entomol. 53: 103– 120. 12. Zehner R, Amendt J, Schutt S, Sauer J, Krettek R, Povolny D (2004) Genetic identification of forensically important flesh flies (Diptera: Sarcophagidae). Int J Legal Med. 118(4): 245–247. 13. Sperling FA, Anderson GS, Hickey DA (1994) A DNA-based approach to the identification of insect species used for postmortem interval estimation. J Fo- rensic Sci. 39(2): 418–427. 14. Harvey AG, Bryant RA, Tarrier N (2003) Cognitive behaviour therapy for post- traumatic stress disorder. Clin Psychol Rev. 23(3): 501–522. 15. Meier R, Zhang G, Ali F (2008) The use of mean instead of smallest interspecific distances exaggerates the size of the “Barcoding Gap” and leads to misidenti- fication. Syst Biol. 57(5): 809–813. 16. Whitworth TL, Dawson RD, Magalon H, Baudry E (2007) DNA barcoding cannot reliably identify species of the blowfly genus Protocalliphora (Diptera: Calli- phoridae). Proc Biol Sci. 274(1619): 1731–1739. 17. Oshaghi MA, Yaghoobi F, Abai MR (2006) Pattern of mitochondrial DNA variation between and within Anopheles stephensi (Diptera: Culicidae) biological forms sug- gests extensive gene flow. Acta Trop. 99(2–3): 226–233. 18. Oshaghi MA, Shemshad K, Yaghobi-Ershadi MR, Pedram M, Vatandoost H, Abai MR, Akbarzadeh K, Mohtarami F (2007) Ge- netic structure of the malaria vector Anopheles superpictus in Iran using mi- tochondrial cytochrome oxidase (COI and COII) and morphologic markers: a new species complex? Acta Trop. 101(3): 241– 248. 19. Karimian F, Oshaghi MA, Sedaghat MM, Waterhouse RM, Vatandoost H, Hanafi- Bojd AA, Ravasan NM, Chavshin AR (2014) Phylogenetic analysis of the Ori- ental-Palearctic-Afrotropical members of Anopheles (Culicidae: Diptera) based on nuclear rDNA and mitochondrial DNA characteristics. Jpn J Infect Dis. 67(5): 361–367. 20. Hashemi-Aghdam SS, Oshaghi MA (2015) A checklist of Iranian cockroaches (Blat- todea) with description of Polyphaga sp as a new species in Iran. J Arthropod Borne Dis. 9(2): 161–175. 21. Hashemi-Aghdam SS, Rafie G, Akbari S, Oshaghi MA (2017) Utility of mtDNA- COI barcode region for phylogenetic re- lationship and diagnosis of five common pest cockroaches. J Arthropod Borne Dis. 11(2): 182–193. 22. Khanzadeh F, Khaghaninia S, Maleki-Rav- asan N, Oshaghi MA, Adler PH (2020) Black flies (Diptera: Simuliidae) of the Aras River Basin: species composition and floral visitation. Acta Trop. 209: 105536. 23. Ratcliffe ST, Webb DW, Weinzievr RA, Robertson HM (2003) PCR-RFLP iden- tification of Diptera (Calliphoridae, Mus- cidae and Sarcophagidae)--a generally applicable method. J Forensic Sci. 48 (4): 783–785. 24. Guo YD, Cai JF, Meng FM, Chang YF, Gu Y, Lan LM, Liang L, Wen JF (2012) Identification of forensically important flesh flies based on a shorter fragment of the cytochrome oxidase subunit I gene in China. Med Vet Entomol. 26(3): 307– 313. 25. Guo Y, Zha L, Yan W, Li P, Cai J, Wu L (2014) Identification of forensically im- portant sarcophagid flies (Diptera: Sar- cophagidae) in China based on COI and period gene. Int J Legal Med. 128(1): 221– 228. 26. Jordaens K, Sonet G, Richet R, Dupont E, Braet Y, Desmyter S (2013) Identification http://jad.tums.ac.ir/ J Arthropod-Borne Dis, December 2020, 14(4): 416–424 F Talebzadeh et al.: Molecular Species Identification of … 423 http://jad.tums.ac.ir Published Online: December 31, 2020 of forensically important Sarcophaga spe- cies (Diptera: Sarcophagidae) using the mitochondrial COI gene. Int J Legal Med. 127(2): 491–504. 27. Meiklejohn KA, Wallman JF, Dowton M (2011) DNA-based identification of fo- rensically important Australian Sarcoph- agidae (Diptera). Int J Legal Med. 125 (1): 27–32. 28. Kim YH, Shin SE, Ham CS, Kim SY, Ko KS, Jo TH (2014) Molecular identifica- tion of necrophagous muscidae and sar- cophagidae fly species collected in Ko- rea by mitochondrial cytochrome C oxi- dase subunit I nucleotide sequences. Sci- entificWorldJournal. 2014: 275085. 29. Zhang C, Fu X, Xie K, Yan W, Guo Y (2015) MtDNA analysis for genetic iden- tification of forensically important sar- cophagid flies (Diptera: Sarcophagidae) in China. J Med Entomol. 52(6): 1225– 1233. 30. Sharma M, Singh D, Sharma AK (2015) Mitochondrial DNA based identification of forensically important Indian flesh flies (Diptera: Sarcophagidae). Forensic Sci Int. 247: 1–6. 31. Talebzadeh F (2017) Insect fauna of hu- man cadavers in Tehran d istrict. [MSc Thesis]. School of Public Health, Tehran University of Medical Sciences, Iran. 32. Zumpt F (1965) Myiasis in Man and An- imals in the Old World, a textbook for physicians, veterinarians and zoologist. Butterwort, London. 33. McAlpine JF (1981) Morphology and ter- minology-adults. In: McAlpine JF, Pe- terson BV, Shewell GE, Teskey HJ, Vock- eroth JR, Wood DM (Eds): Manual of Nearctic Diptera. Vol. 1. Research Branch, Agriculture Canada, Monograph No. 27, pp. 9–63. 34. Folmer O, Black M, Hoeh W, Lutz R, Vri- jenhoek R (1994) DNA primers for am- plification of mitochondrial cytochrome c oxidase subunit I from diverse meta- zoan invertebrates. Mol Mar Biol Bio- technol. 3(5): 294–299. 35. Larkin MA BG, Brown NP, Chenna R, McGettigan PA, McWilliam H, Valentin F, Wallace IM, Wilm A, Lopez R, Thompson JD, Gibson TJ, Higgins DG (2007) Clustal W and Clustal X version 2.0. Bioinformatics. 23(21): 2947–2948. 36. Akbarzadeh K, Saghafipour A, Jesri N, Karami-Jooshin M, Arzamani K, Haz- ratian T, Kordshouli RS, Afshar AA (2018) Spatial distribution of necro- phagous flies of infraorder Muscomor- pha in Iran using geographical infor- mation system. J Med Entomol. 55(5): 1071–1085. 37. Crozier RH, Crozier YC (1993) The mito- chondrial genome of the honeybee Apis mellifera: complete sequence and genome organization. Genetics. 133(1): 97–117. 38. Hebert PD, Cywinska A, Ball SL, deWaard JR (2003) Biological identifications through DNA barcodes. Proc Biol Sci. 270(1512): 313–321. 39. Hebert PD, Ratnasingham S, deWaard JR (2003) Barcoding animal life: cytochrome c oxidase subunit 1 divergences among closely related species. Proc Biol Sci. 270(Suppl 1): S96–9. 40. Ball SL HP, Burian SK, Webb JM (2005) Biological identifications of mayflies (Ephemeroptera) using DNA barcodes. J North Am Benthol Soc. 24: 508–524. 41. Tan M, Alhajj R, Polat F (2010) Scalable approach for effective control of gene regulatory networks. Artif Intell Med. 48 (1): 51–59. 42. Jafari S, Oshaghi MA, Akbarzadeh K, Abai MR, Koosha M, Mohtarami F (2019) Identification of forensically important flesh flies using the cytochrome c oxi- dase subunits I and II genes. J Med En- tomol. 56(5): 1253–1259. 43. Bakhshi H, Oshaghi MA, Abai MR, Rassi Y, Akhavan AA, Mohebali M (2013) MtDNA CytB structure of Rhombomys http://jad.tums.ac.ir/ J Arthropod-Borne Dis, December 2020, 14(4): 416–424 F Talebzadeh et al.: Molecular Species Identification of … 424 http://jad.tums.ac.ir Published Online: December 31, 2020 opimus (Rodentia: Gerbellidae), the main reservoir of cutaneous leishmaniasis in the borderline of Iran-Turkmenistan. J Arthropod Borne Dis. 7(2): 173–184. 44. Oshaghi MA, Rafinejad J, Choubdar N, Piazak N, Vatandoost H, Telmadarraiy Z (2011) Discrimination of relapsing fe- ver Borrelia persica and Borrelia mi- crotti by diagnostic species-specific pri- mers and polymerase chain reaction-re- striction fragment length polymorphism. Vector Borne Zoonotic Dis. 11(3): 201– 207. 45. Oshaghi MA, Ravasan NM, Hide M, Ja- vadian EA, Rassi Y, Sedaghat MM (2009) Development of species-specific PCR and PCR-restriction fragment length poly- morphism assays for L. infantum/L. do- novani discrimination. Exp Parasitol. 122 (1): 61–65. 46. Oshaghi MA, Rasolian M, Shirzadi MR, Mohtarami F, Doosti S (2010) First re- port on isolation of Leishmania tropica from sandflies of a classical urban cuta- neous leishmaniasis focus in southern Iran. Exp Parasitol. 126(4): 445–450. 47. Rassi Y, Oshaghi MA, Azani SM, Abai MR, Rafizadeh S, Mohebai M (2011) Molecular detection of Leishmania in- fection due to Leishmania major and Leishmania turanica in the vectors and reservoir host in Iran. Vector Borne Zo- onotic Dis. 11(2): 145–150. 48. Oshaghi MA, Chavshin AR, Vatandoost H (2006) Analysis of mosquito blood- meals using RFLP markers. Exp Parasi- tol. 114(4): 259–264. 49. Anjomruz M, Oshaghi MA, Sedaghat MM, Pourfatollah AA, Raeisi A, Vatandoost H (2014) ABO blood groups of residents and the ABO host choice of malaria vec- tors in southern Iran. Exp Parasitol. 136: 63–67. http://jad.tums.ac.ir/