Abbreviations
CO1: Cytochrome oxidase I
PCR: Polymerase chain reaction
mtCO1: mitochondrial cytochrome oxidase subunit I
MBST: Molecular Biological System Transfer
Introduction
Oestrosis is nasal myiasis, resulting from infestation with the larvae of flies belonging to the genus Oestrus (Diptera: Oestridae). It is regarded as a severe parasitosis in goats and sheep and sometimes in other animal species [ 1 ]. Nasal discharge and sneezing are the most common clinical signs in animals [ 2 ]. Parasitic larvae seriously affect the host's well-being and cause the loss of milk, meat, and wool [ 3 ]. The sheep bot fly (Oestrus ovis L.) causes severe cavitary myiasis in all areas worldwide where goats and sheep are reared [ 4 ]. O. ovis can be diagnosed by the morphological characterization of the larvae (e.g., the slits of the posterior spiracles located on the posterior spiracular plates), clinical presentation, and infrequent adult fly identification [ 5 ]. A few years ago, molecular techniques were extensively used to achieve insights into the systematics and taxonomy of various insects. Particularly, hypervariable regions within ribosomal and mitochondrial DNA are robust genetic markers to tackle taxonomic diagnostic issues for some flies that cause myiasis and belong to the Oestridae family [ 6 , 8 ]. The mtCO1 gene was chosen to characterize the Oestridae larval relationships. It has previously been revealed to be necessary for various molecular phylogenetic targets due to its comparatively enormous size and the inclusion of both highly conserved and variable regions with different mutational rates [ 9 ]. The CO1 658 bp region is a standard and broadly considered animal taxa universal marker [ 10 ]. Over the past decade, the barcoding of DNA has provided an effective instrument for phylogenetic characterization using CO1 sequences and molecular identification of various insect species comprising the Oestridae family [ 11 , 14 ]. Therefore, the current study aimed to utilize PCR and partial mtCO1 gene sequencing as a diagnostic tool to identify and demonstrate the phylogenies of larvae in different climatic regions of Iran. The molecular detection of O. ovis is substantial to address its epidemiology and control.
Results
Morphological identification
The identification of the mature larvae of O. ovis was based on morphological characteristics. The third instar larva is up to 20 mm long, yellowish‐white, and tapering anteriorly. Every fragment has an obscure lateral band on the dorsal side. They have a big black oral hook linked to the internal cephalopharyngeal skeleton. The ventral surface is lined with small spines. The results corresponded with what is mentioned in the identification key [ 15 ].
PCR amplification and sequencing of mtCOI gene
Figure 1 shows that the isolates in this study had a similar sequence alignment and they all are placed in the same consensus modification group. The 610 bp fragment of mtCO1 from all isolates was successfully amplified (Figure 1). The 610 bp CO1 region covered by the consensus sequences from all reads was effectively acquired with high-quality chromatogram scores. The sequence of the amplified fragments revealed a genome arrangement typical of the O. ovis sequence in GenBank. Phylogenetic tree elucidated phylogenetic relationships between other members of the family Oestridae published in GenBank and O. ovis larvae amplicons as shown in Figure 3. The mean reported interspecific similarity value of CO1s in the Oestrinae subfamily was 86.7% [ 16 ]. However, the sequence indicated a 99% identity using the CO1 sequences of O. ovis available at GenBank and < 91% identity using the CO1 sequences of other Oestrus species (Oestrus sp.) [ 17 ]. Furthermore, this sequence identity using the CO1 sequences of other Oestrid species was < 90%. These findings established that the examined larvae belonged to O. ovis.
Sequencing results of Iranian O. ovis isolates from six sheep sinuses in this study showed no difference in partial CO1 genes among six Iranian isolates, and they had identical gene sequences (Figure 3). The CO1 gene sequences from the other mitochondrial genomes available in GenBank were chosen and aligned. The obtained sequences in this study (MZ972997 - MZ973002) were compared with GenBank sequences.
Figure 2 demonstrates the comparison of the isolated sequences in this study with sequences in the gene bank, which shows the differences and similarities in the alignment process. Based on multiple nucleotide sequence alignments (Figure 2) and sequence identity matrix algorithms, the percent identities among O. ovis isolates in this study were 100% with O. ovis from Iraq (MW555828), Bosnia and Herzegovina (MG755264), and Croatia (MN845130), and 99% with O. ovis from Iran (KX268655), Spain (MW145179), and Brazil (KR820703). The nucleotide differences with other isolates are exhibited in Figure 3. As a result, a phylogenetic tree was inferred with CO1 gene sequences from this study and previously published correlated species of mitochondrial sequences. The phylogenetic reconstruction showed two main branches (A and B) and a sub-branch (A1) that contained our isolate (MZ972997 - MZ973002) and homolog sequences with high similarity from Iraq, Bosnia and Herzegovina, and Croatia (Figure 3).
About 100% similarity of Iranian sheep isolates with isolates from Iraq, Bosnia and Herzegovina, and Croatia, shows a very close relationship, and homology between the nucleotide sequences of O. ovis isolates from sheep sinuses mtCO1 gene in these geographical areas. Figure 4 shows the alignment process between the isolates of this study and the isolates of the same genus and species identified from the GenBank.
Discussion
This study analyzed the genetic diversity of O. ovis based on mtCO1 gene nucleotide sequences isolated from the sinuses of Iranian sheep. The mtDNA provides more accurate identifications of the species to which larvae belong, and even closely related species can have different developmental characteristics [ 18 , 19 ]. Phylogenetic analysis demonstrated that none of the sheep sequences obtained in this study had intraspecific variations and high similarity with the sequences of O. ovis isolates from sheep mtCO1 gene from Iraq, Bosnia and Herzegovina, Croatia, Brazil, and Spain (Figure 3).
The results of the nucleotide alignment of MZ972997 - MZ973002 with other sequences revealed that variation in the A chain of the phylogenetic tree between O. ovis sequences was in position 258 (T > C) for KR820703, MW145179, KX268655, and MT124626, 372 (G > A) for MT124626, 435 (T > A) for MW555828, 456 (C > A) for MT124626, and 459 (T > C) for MW145179, MN845130, MG755264, and MT124626 (Figure 4). The genes encoding proteins in the mtDNA are more diverse than the ribosomal genes [ 20 ]. Based on molecular alignment and phylogenetic tree analysis in this study, the close relationship of Iranian isolates with Iraq, Bosnia and Herzegovina, and Croatia from the same hosts was maybe correlated with the same geographical conditions and the high rate of trade between countries [ 21 , 22 ]. Furthermore, the results of this study showed that all O. ovis obtained in this research and from GenBank were grouped into a monophyletic group according to the results of the previous work [ 23 , 25 ]. Molecular studies are being used to investigate taxonomy and perform substantial status for identifying Oestridae larvae that lead to obligate myiasis [ 26 ].
The physical traits of the larvae and morphological identification are unreliable [ 27 ]. Therefore, the molecular identification of myiasis-causing fly species may be considered a credible substitutional to morphological identification as it is challenging to recognize larvae to the genus level [ 27 , 29 ]. In this study, the high similarity between the isolates from Iran and other countries isolated from the same host and interestingly between other species from other hosts were observed. The results revealed that the genetic variability of O. ovis isolated from different geographical areas of Iran was not clear based on a 610 bp fragment of the mtCO1 gene. Although a low rate of intraspecific variation existed in O. ovis, especially in different hosts, more O. ovis isolates should be collected for further investigations, and complete gene sequences from isolates would be required. The differences between O. ovis isolates may result from variations in the animal breed, geographical area, and climate [ 22 , 24 , 25 , 31 , 32 ].
Bosly, in 2018, reported a molecular identification of O. ovis larvae, which was recognized as a reliable replacement for morphological identification [ 25 ]. This research gives the initial molecular depiction data on O. ovis larvae in sheep in Iran using the barcode CO1 sequences. Our results showed a close association between O. ovis from different countries, even though limited sequences were included in the dataset. The barcode CO1 sequence successfully analyses the species of the subfamilies in Oestriade, and this result shows the usefulness of the CO1 as a suitable diagnostic molecular marker for identifying and characterizing bot fly species.
Collection and identification of larvae
O. ovis larvae were collected from the sinuses of infected sheep in the slaughterhouse of three climatic regions of Iran, including the coastal areas of the Caspian Sea, mountain plateau areas, and flat plains around the Persian Gulf. Two provinces were selected from each region, and ten samples were collected from each of these provinces. Collected larvae were conserved in 70% C2H5OH for morphometric identification and molecular survey. The larvae were identified based on their morphological traits and the three slits at the posterior larvae’s spiracle (Figure 5) [ 15 , 33 ].
Morphological identification and microscopic examination
The mature larvae of O. ovis were identified based on the morphological characteristics of the larval stages (Figure 6). The third instar larva is 20 mm long and yellow when young, changing to light brown later, and in the maturing stage finally, bands are appear on the dorsal of the third-stage larva. It shows wide transverse blackish bands dorsally. The second segment has a variable number of small teeth on the dorsal side with the rest of the segments being bare. Ventrally, the segments bear rows of strong spines (Figures 6A and B). The findings corresponded with what was mentioned in the identification key [ 15 ].
DNA extraction
DNA was extracted from larvae using a DNA isolation kit (MBST, Tehran, Iran) according to the manufacturer's instructions.
PCR procedures and sequencing
Extracted DNA from each sample was used for amplifying a 610 bp fragment of the mtCOI gene barcoding region using the FFCOI forward (5'-GGAGCATTAATYGGRGAYG-3') (34) and R-HCO reverse (5'- TAAACTTCAGGGTGACCAAAAATCA-3') primers (35). The PCR reaction was a total of 25 µl containing 1 µl of DNA template, 0.25 mM dNTPs, 1x enzyme buffer, 0.25 mM MgCl2, 2.5 µM primers (FFCOI / R-HCO), and 0.5 U of Taq DNA polymerase. The PCR conditions were 40 cycles of initial denaturation at 94°C for 30 sec, followed by denaturation at 94°C for 30 sec, annealing at 52°C for 30 sec, and extension at 72°C for 30 sec. At the final extension step, the reaction terminated at 72°C for 10 min. During PCR, control reactions were performed without adding DNA to the reaction to exclude contamination. To assess the amplified fragment, the amplicons from the PCR were electrophoresed on 1.5% agarose gel and were compared with a marker of standard molecular weight. All amplicons were sequenced using an ABI-3730XL capillary machine (Macrogen Inc., South Korea). O. ovis COI nucleotide sequences were trimmed for any error in nucleotides and the primer sequences were removed using DNA SISMax 3.0 software and then, were submitted to GenBank (Table 1). Based on nucleotide sequences, Bioedit 7.0 software was used for multiple sequence alignments with the Clustal W method. Using the Neighbor-Joining method, phylogenetic trees, and statistical analyses were conducted. Moreover, the evolutionary distances were computed utilizing the p-distance method with the MEGA X software. The analysis involved 36 nucleotide sequences with several 1000 bootstrap replications. Sequences were compared with the identities of isolates confirmed using the NCBI BLAST database (). The nucleotide sequences were retrieved from GenBank and were utilized to paradigm a phylogenetic tree.
Isolate | Country-Provenience | Host | Nucleotide accession number |
---|---|---|---|
A-COI | Iran- Rasht | Sheep | MZ972997 |
B-COI | Iran-Urmia | Sheep | MZ972998 |
C-COI | Iran-Shiraz | Sheep | MZ972999 |
D-COI | Iran-Kermanshah | Sheep | MZ973000 |
E-COI | Iran-Sari | Sheep | MZ973001 |
F-COI | Iran-Ahvaz | Sheep | MZ973002 |
Authors' Contributions
M.T. designed the study. B.E. helped for designing the study and writing of manuscript. S.R. sampling and laboratory examination. All authors reviewed the manuscript.
Acknowledgements
The authors acknowledge the financial support provided by the Vice-Chancellor for Research and Technology, Urmia University. We are also very grateful to Dr. Mohammad Ali Bayatzadeh and Dr. Awat Samiei for their technical assistance.
Competing Interests
The authors declare that they have no competing interests.
References
- Horak I. Parasites of domestic and wild animals in South Africa. XLVI 1-Oestrid fly larvae of sheep, goats, springbok, and black wildebeest in the Eastern Cape Province. Onderstepoort J Vet Res. 2005; 72(4): 315-20. DOI
- Hanan BA. Seasonal prevalence of Oestrus ovis L.(Diptera: Oestridae) larvae in infested sheep in Jazan Region, Saudi Arabia. J Parasitol Vector Biol. 2013; 5(5): 66-71. DOI
- Angulo-Valadez CE, Scholl PJ, Cepeda-Palacios R, Jacquiet P, Dorchies P. Nasal bots a fascinating world! Vet Parasitol. 2010; 174(1-2):19-25. DOI
- Dorchies PH, Tabouret G, Hoste H, Jacquiet PH, Colwell DD, Hall MJR, Scholl PJ. Oestrinae host-parasite interactions. (éd.): The Oestrid flies: biology, host-parasite relationships, impact and management, CABI, UK. 2006; 191-200. DOI
- Ogo NI, Onovoh E, Okubanjo OO, Galindo RC, Lastra JMP. Molecular identification of Cordylobia anthropophaga Blanchard (Diptera: Calliphoridae) larvae collected from dogs (Canis familiaris) in Jos South, Plateau State, Nigeria. Onderstepoort Vet Res. 2012; 79(1):1-4. DOI
- Otranto D, Traversa D. Molecular evidence indicating that Przhevalskiana silenus, Przhevalskiana aegagri, and Przhevalskiana crossii (Diptera, Oestridae) are one species. Acta Parasitol. 2004; 2(49): 173-6.
- Otranto D, Colwell DD, Traversa D, Stevens JR. Species identification of Hypoderma (Diptera, Oestridae) affecting domestic and wild ruminants by Scanning Electron Microscopy (SEM) and molecular characterization of the cytochrome oxidase I (COI) gene. Med Vet Entomol. 2003; 17;: 316-25. DOI
- Otranto D, Traversa D, Colwell DD, Guan G, Giangaspero A, Boulard C, Yin H. A third species of Hypoderma (Diptera, Oestridae) affecting cattle and yaks in China: molecular and morphological evidence. J Parasitol. 2004; 90(5): 958-65. DOI
- Lunt DH, Zhang DX, Szymura JM, Hewitt GM. The insect cytochrome oxidase I gene: evolutionary patterns and conserved primers for phylogenetic studies. Insect Mol Biol. 1996; 5(3): 153-65. DOI
- Hebert PD, Stoeckle MY, Zemlak TS, Francis CM. Identification of birds through DNA Barcodes. PLoS Biol. 2004; 2(10):e312. DOI
- Moreno V, Romero-Fernández I, Marchal JA, Beltrán M, Granados JE, Habela MA, Tamadon A, Rakhshandehroo E, Sarasa M, Pérez JM, Sánchez A. Molecular characterization of bot flies, Oestrus spp (Diptera, Oestridae), from domestic and wild Bovidae hosts. Vet Parasitol. 2015; 212(3-4): 473-7. DOI
- Cavallero S, Pombi M, Perrone V, Milardi GL, D’Amelio S, Giuliani C, Gabrielli S. Gasterophilus intestinalis (Diptera: Oestridae) in the diaphragmatic muscle: An unusual finding. Vet Parasitol. 2017; 237: 117-21. DOI
- Rakhshandehroo E, Razavi SM, Farzaneh R, Esmailnejad A, Asadpour M, Shams S. Phylogenetic analysis of goat warble fly (Przhevalskiana silenus) based on mitochondrial COI gene. J Parasit Dis. 2019; 43(2): 304-7. DOI
- Li XY, Pape T, Zhang D. Gasterophilus flavipes (Oestridae: Gasterophilinae) A horse stomach bot fly brought back from oblivion with morphological and molecular evidence. PLoS One. 2019; 14(8):e0220820. DOI
- Zumpt F. Myiasis in Man and Animals in the Old World: A Textbook for Physicians, Veterinarians, and Zoologists; 1965.
- Otranto D, Milillo P, Traversa D, Colwell DD. Morphological variability and genetic identity in Rhinoestrus spp. Causing horse nasal myiasis. Med Vet Entomol. 2005; 19(1): 96-100. DOI
- Moreno V, Romero-Fernández I, Marchal JA, Beltrán M, Granados JE, Habela MA, Tamadon A, Rakhshandehroo E, Sarasa M, Pérez JM, Sánchez A. Molecular characterization of bot flies, Oestrus spp., (Diptera, Oestridae) from domestic and wild Bovidae hosts. Vet Parsitol. 2015; 212(3-4): 473-77. DOI
- Pape T. Phylogeny of Oestridae (Insecta: Diptera). Syst Entomol. 2001; 26(2): 133-171. DOI
- Guo Y, Cai J, Chang Y, Li X, Liu Q, Wang X, Wang X, Zhong M, Wen J, Wang J. Identification of forensically important sarcophagid flies (Diptera: Sarcophagidae) in China, based on COI and 16S rDNA gene sequences. J Forensic Sci. 2011; 56(6): 1534-40. DOI
- Tan X, Lu H, Sun Y, Chen X, Wang D, Jia F, Zhang K. Expanding the materials space of DNA via organic-phase ring-opening metathesis polymerization. Chem. 2019; 5(6): 1584-96. DOI
- Banafshi O, Hanafi-Bojd AA, Karimi M, Faghihi F, Beik-Mohammadi M, Gholami S, Javaherizadeh S, Edalat H, Vatandoost H, Telmadarraiy Z. Tick ectoparasites of animals in borderline of Iran-Iraq and their role on disease transmission. J Arthropod Borne Dis. 2018; 12(3): 252. DOI
- Kettle PR. A study on the sheep botfly, Oestrus ovis (Diptera: Oestridae) in New Zealand. N Z Entomol. 1973; 5(2):185-191. DOI
- Metwally DM, Albasyouni SA, Barakat IAH, Al-Turaiki IM, Almuhanna AM, Bashir MA, Baghdadi HB, El-Khadragy MF, Alajmi RA. Prevalence Rate and Molecular Characteristics of Oestrus ovis L.(Diptera, Oestridae) in Sheep and goats from Riyadh, Saudi Arabia. Animals. 2021; 11(3): 689. DOI
- Negm-Eldin MM, Elmadawy RS, Hanan GM. Oestrus ovis larval infestation among sheep and goats of Green Mountain areas in Libya. J Adv Vet Anim Res. 2015; 2(4): 382-7. DOI
- Bosly AH. Molecular identification studies on Oestrus ovis L. (Diptera: Oestridae) larvae-infested sheep in Jazan region, Saudi Arabia. Indian J Anim Res. 2018; 52(1): 105-10. DOI
- Otranto D, Stevens JR. Molecular approaches to the study of myiasis-causing larvae. Int J Parasitol. 2002; 32(11): 1345-60. DOI
- Ogo NI, Onovoh E, Okubanjo OO, Galindo RC, Lastra JMP. Molecular identification of Cordylobia anthropophaga Blanchard (Diptera: Calliphoridae) larvae collected from dogs (Canis familiaris) in Jos South, Plateau State, Nigeria. Onderstepoort J Vet Res. 2012; 79(1): 01-04. DOI
- Ames C, Turner B, Daniel B. The use of mitochondrial cytochrome oxidase I gene (COI) to differentiate two UK blowfly species–Calliphora vicina and Calliphora vomitoria. Forensic Sci Int. 2006; 164(2): 179-82. DOI
- Hall M, Adams Z, Wyatt N, Testa J, Edge W, Nikolausz M, Farkas R, Ready P. Morphological and mitochondrial DNA characters for identification and phylogenetic analysis of the myiasis causing flesh fly Wohlfahrtia magnifica and its relatives, with a description of Wohlfahrtia monegrosensis sp. n. Wyatt & Hall. Hall Med Vet Entomol. 2009; 23(s1): 59-71. DOI
- Wang Y, Liu GH, Li JY, Xu MJ, Ye YG, Zhou DH, Song HQ, Lin RQ, Zhu XQ. Genetic variability among Trichuris ovis isolates from different hosts in Guangdong Province, China revealed by sequences of three mitochondrial genes. mtDNA. 2013; 24(1): 50-54. DOI
- Aleix-Mata G, Peréz JM, Sánchez A. The complete mitochondrial genome of Oestrus ovis (Linnaeus, 1758)(Diptera: Oestridae). MitDNA Part B. 2021; 6(7):1847-48. DOI
- Hall MJ, Smith KG. Diptera causing myiasis in man. Medical insects and arachnids, Springer. 1993; 429-69. DOI
- Blacket MJ, Semeraro L, Malipatil MB. Barcoding Queensland fruit flies (Bactrocera tryoni): impediments and improvements. Mol Ecol Resour. 2012; 12(3): 428-36. DOI
- Folmer O, Black M, Hoeh W, Lutz R, Vrijenhoek R. DNA primers for amplification of mitochondrial cytochrome c oxidase Subunit 1 from diverse metazoan invertebrates. Mol Marine Biol Biotechnol. 1994; 3(5): 294-9. PubMed
Send comment about this article