Molecular Characterization and Phylogeny Analysis Based on Sequences of Cytochrome Oxidase gene From Hemiscorpius lepturus of Iran

Document Type : Genetics - Immunology

Author

Department of Basic Sciences, Biochemistry and Molecular Biology Section, Faculty of Veterinary Medicine, Shahid Chamran University of Ahvaz, Ahvaz, Iran

Abstract

Abstract:
Background: Hemiscorpius lepturus is a medically important scorpion found along the Iranian borders, especially near to Khuzestan Province in the south-west of Iran. This is the only non-buthid scorpion which is potentially lethal in southern Iran and is responsible for severe dermonecrotic scorpionism.
OBJECTIVES: In this study, DNA fragment of the mitochondrial cytochrome c oxidase subunit I (COXI) gene of H. Lepturus for the molecular phylogenetic analysis was amplified.
METHODS: We amplified  a 624 bp gene fragment of cytochrome C oxidase subunit I (COXI) from H. lepturus collected from Khuzestan, Ahvaz by PCR. After sequencing of the PCR products, the phylogenetic analysis was performed using the neighbor-joining method with 1000 replicates of bootstrapping using the MEGA7 software.
RESULTS: The results of phylogenetic analysis revealed four distinct clusters (A1, A2, B and C) belonging to the family of Hemiscorpionidae that were grouped together with bootstrap score between 77-96%. The gene fragment of Hl-Kh formed a cluster relative to the only scorpion of H. Lepturus (Hl) from Izeh, Iran with a good bootstrap score of 96. These two samples isolated from the Khuzestan province of Iran are closely related, as they clustered together as Cluster B. The genetic distances of Hl-Kh among the cluster A1, A2 and C ranged from 16 to19% and the lowest interspecific distance was in cluster B between Hl-Kh and Hl (6%).
CONCLUSIONS: Hl-Kh isolated from the Khuzestan province of Iran are closely related to Hl, as they showed the greatest interspecific variation observed in Hemisccorpius genus in this study. Although the sample size in Cluster B is not large enough to draw a final conclusion, the percentage of sequence divergence was high enough for interspecific comparisons to provide separation of species.

Keywords


Article Title [Persian]

شناسایی ملکولی و آنالیز فیلوژنی عقرب همیسکورپیوس لپچروس ایران بر اساس توالی ژن سیتوکروم C اکسیداز

Author [Persian]

  • عباس جلودار
گروه علوم پایه بخش بیوشیمی و بیولوژی مولکولی ، دانشکده دامپزشکی، دانشگاه شهید چمران اهواز، اهواز، ایران
Abstract [Persian]

زمینه مطالعه: همیسکورپیوس لپچروس یک عقرب با اهمیت از نظر پزشکی است که بموازات مرز ایران بویژه در استان خوزستان واقع در جنوی غرب ایران یافت می شود. این تنها گونه عقرب از خانواده غیر بوتیده است که در جنوب ایران علاوه بر اینکه بطور بالقوه‌ای کشنده است، عامل اصلی درماتیت شدید ناشی از گزش عقربها نیز می باشد.
هدف: هدف از این مطالعه شناسایی ملکولی و آنالیز فیلوژنیک عقرب همیسکورپیوس لپچروس با استفاده از ژن سیتوکروم اکسیداز C بود.
روش کار: قطعه ژن 624 نوکلئوتیدی سیتوکروم اکسیداز C زیر واحد 1 (COXI) با استفاده از آغازگرهای اختصاصی بوسیله واکنش پلیمراز زنجیره ای تکثیر و سپس توالی یابی گردید.
نتایج: پس از توالی یابی محصول PCR، مطابق نتایج فیلوژنیک، توالی های مشابه از خانواده همیسکورپیونیده در چهار گروه مشخص (A1 ,A2 ,B and C) با در صد مشابهت در محدوده بین 96-77 گروه بندی شدند. در این گروه بندی، قطعه ژن تکثیر شده Hl-Kh با ژن همولوگ عقرب همیسکورپیوس لپچروس جدا شده از ایذه در استان خوزستان ایران که در گروه B قرار گرفتند نزدیکی نسبتاً بالایی با یکدیکر نشان می‌دادند (96). فاصله ژنتیکی قطعه ژن Hl-Kh با گروه‌های A1، A2 و C در محدوده 19-16 درصد قرار داشت، در حالیکه کمترین فاصله ژنتیکی این قطعه ژن با قطعه جدا شده Hl در گروه B بمیزان 6 درصد گزارش گردید.
نتیجه گیری نهایی: هرچند که با توجه به کوچکی نمونه‌ها در گروه B، نمی‌توان به یک نتیجه گیری قطعی رسید، ولی درصد اختلاف ژنتیکی ابن دو قطعه به اندازه کافی زیاد هست که بتوان Hl-Kh را یک تفاوت داخل گونه ای انگاشت.

Keywords [Persian]

  • همیسکورپیوس لپچروس
  • عقرب
  • فیلوژنتیک
  • سیتوکروم اکسیداز C زیرواحد 1
  • زهر

Hemiscorpius lepturus (Scorpions: Hemiscorpionidae) is a medically important scorpion found along the Iranian borders, especially adjacent to Khuzestan province in the south-west of Iran. H. lepturus venom is known to have extremely toxic effects on the central nervous system and cardiovascular system (Radmanesh, 1990, 1998). This is the only dangerous and potentially deadly non-Buthid scorpion which is responsible for severe ulceration of the skin and also hemolysis of blood cells during the warm months of the year. Although only 10-12% of total scorpion stings during the hot season in Khuzestan province are due to this species, it is responsible for more than 95% of deaths (Manouchehrifar et al., 2013). H. lepturus seems to be restricted to Iraq and the western and south western regions of Iran. However, the distribution area of this species also includes the west of Pakistan (Fet, 2000). Three species belong to the family Hemiscorpiidae that are known to exist in in Iran are H. gaillardi, H. lepturus and H. persicus. These species are morphologically very close to each other and difficult to distinguish for a non-skilled person, therefore we assume that H. lepturus is probably not the only species responsible for significant scorpionism problems in the southern provinces of Iran.

DNA sequences of the mitochondrial cytochrome oxidase I (COXI) gene have been one of the most widely used molecular markers for identifying all kinds of animals. This is a key enzyme in aerobic metabolism. In fact, the evolution of this gene is rapid enough to allow the discrimination of not only closely related species, but also phylogenetic variations inside a particular species (Hebert et al., 2003). Phylogenetic analysis using COXI gene sequences were extensively carried out by several workers in different groups of organisms like the genus members of Puntius belongs to Cyprinidae family (Garg et al., 2017), 42 species of Culicoides (Diptera, Ceratopogonidae) from three continents (Augot et al., 2017 and Linto et al., 2002), white backed plant hopper Sogatella furcifera (Sreejith  and Sebastian, 2014) and Commercially Available Pangasius hypophthalmus in Italy (Bellagamba et al., 2015). This gene is widely accepted as a DNA barcode for the accurate and easy identification of species.

In this study, DNA sequence of the mitochondrial cytochrome c oxidase subunit I (COXI) gene of H. Lepturus for phylogenetic analysis was amplified.

Materials and Methods

Scorpion samples: Iranian scorpions H. lepturus were collected from Khuzestan Province in Iran. The specimens were obtained from the reference laboratory of the Razi Vaccine and Serum Research Institute in 95% ethanol and were frozen at -20 °C. An individual scorpion of Mesobuthus eupeus served as outgroup for phylogenetic analyses.

Genomic DNA extraction: Scorpion tissue (0.5-1.0 g) was homogenized in 10 ml of RSB buffer (10 mM Tris–HCl pH 7.4, 10 mM NaCl/25 mM EDTA, SDS 1%). The DNA was extracted with an equal volume of phenol/chloroform, and the aqueous phase was then extracted with chloroform.  Precipitated Nucleic Acid is pelleted out by centrifuging at 12000 rpm for 15min before the pelleted DNA was washed with 70% ethanol, air dried and resuspended in 50μl TE buffer. DNA concentration and purity were determined by calculating the absorbance ratio A260/280.

PCR amplification and sequencing: PCR amplification was carried out on aliquots of the H. Lepturus genomic DNA as template with initial denaturation for 5 min at 95 °C, followed by 35 cycles of 45 s at 94 °C, 1 min at 48 °C, and 1 min at 72 °C and, finally, 7 min of incubation at 72 °  C in a final reaction volume of 25 µl containing 50 mM KCl, 10mM Tris-HCl (pH 8.3), 20 pmol of each primer, deoxynucleotides (each at 220 µM), 1.5 mM MgCl2, 100 ng genomic DNA, and 1 U Taq polymerase. The primers were COXI-F 5’-GGTCAACAAATCATAAAGATATTGG and COXI-R 5’- TAAACTTCAGGGTGACCAAAAAATCA (Folmer et al., 1994). The Amplified PCR products were electrophoresed through a 1% agarose gel, and stained with DNA Safe Stain (Sinaclon, Iran) before detection by UV transillumination. The amplified DNA fragments were extracted from agarose gel prior to performing DNA sequencing according to the dideoxy termination method using an automated Applied Biosystems 373 DNA sequencer. The sequence was determined for both strands by using overlapping fragments. 

DNA analysis: The DNA sequences comparisons and putative functions of the DNA fragments were done using the blastn and blastx algorithms programs in the Nartional Center for Biotechnology Information (NCBI) GenBank database (Altschul et al., 1997). Similar sequences are downloaded from the NCBI database. The sequence obtained after sequencing was translated and its secondary structure of protein was predicted using the Prediction of Transmembrane Regions and Orientation server (embnet.vital-it.ch/software/TMPRED_form.html). The phylogenetic analysis was performed using the neighbor-joining method with 1000 replicates of bootstrapping using the MEGA7 software. A COXI sequence of Mesobuthus eupeus was used as outgroup.

Results

Amplification and sequence analysis. Starting with 0.05 g of fresh tissue, 5 µg of genomic DNA was obtained. The results obtained by PCR analysis confirmed the presence of COXI in H. lepturus since the expected 624 bp cDNA fragment was amplified (Fig. 1). No amplification product was obtained in the control samples where template DNA was excluded. According to the sequencing results (Fig. 2), the DNA contained a single open reading frame of 208 amino acids. The Blastx search indicated that the sequenced segment is closely related to different scorpion species from the order chelicerata. Domain analysis of Hl-Kh fragment showed e-value of 6.81e-93 with conserved domain of cytochrome c oxidase subunit I (accession number MTH00142) with good query coverage. This multiple- domain contains Heme-copper oxidase subunit I Superfamily with e-value 3.20e-50. The PSIPRED protein sequence analysis method suggested the secondary structure of the protein. Each amino acid is allocated values for alpha helix, beta sheet and coils using a window of 7 amino acids. The likelihood of a given amino acid calculated, and the secondary structure with the largest conformation probability is determined by the amino acid. We obtained a prediction of 60% alpha-helical for the coding sequence of COXI from H. lepturus (Fig. 3). The secondary structure data (Table 1) were used to enable a consensus confirmation of predicted transmembrane (TM) helices.

Phylogenetic analysis: By comparing the sequence of Hl-Kh with the NCBI database, a total of 13 sequences including Hemiscorpius lepturus cytochrome oxidase subunit I (COI) gene (Hl; KU341987), Diplocentrus tehuacanus voucher AMCC:LP2044 cytochrome oxidase subunit I (COI) gene  (Dt; KM514661), Brachistosternus ferrugineus isolate LP_7632 cytochrome c oxidase subunit I (COI) gene (Bf; KT446962), Scorpio maurus fuliginosus isolate Sc51 cytochrome oxidase subunit I (COI) gene (Smf; FJ198060), Brachistosternus paposo isolate S027 cytochrome c oxidase subunit 1 (COI) gene (Bp; KX517218), Brachistosternus roigalsinai isolate S021 cytochrome c oxidase subunit 1 (COI) gene (Br; KX517215), Pandinus cavimanus cytochrome oxidase subunit I (COI) gene (Pc; AY156580), Brachistosternus simoneae isolate LP_10708 cytochrome c oxidase subunit I (COI) gene (Bs; KT447023), Brachistosternus ehrenbergii isolate LP_9279 cytochrome c oxidase subunit I (COI) gene (Be; KT446957), Diplocentrus bereai voucher AMCC:LP6532 cytochrome oxidase subunit I (COI) gene (Db; KM514638), Diplocentrus reddelli voucher AMCC:LP10981 cytochrome oxidase subunit I (COI) gene (Dr; KM514657),  Diplocentrus cozumel voucher AMCC:LP4102 cytochrome oxidase subunit I (COI) gene (Dc; KM514642), Scorpio fuscus isolate 12228H cytochrome oxidase subunit I (COI) gene (Sf; KT188243) and Mesobuthus eupeus cytochrome c oxidase subunit I gene (Me-Kh) were retrieved. The nucleotide sequence of Hl-Kh with the retrieved sequence data was subjected to multiple sequence alignment and phylogenetic analysis using MEGA7 software (Fig. 4). Almost all the internal branches showed high bootstrap value bigger than 80. Branches corresponding to partitions reproduced in less than 50% bootstrap replicates are collapsed. The tree is drawn to scale, with branch lengths in the same units as those of the evolutionary distances used to conclude the phylogenetic tree. M. eupeus (Me) used as an out group was clearly branched separately, which signified the credibility of constructed phylogram. This phylogeny allows addressing several important taxonomic issues relating to the Scorpionidae family. The constructed phylogram showed four distinct clusters (A1, A2, B and C). The phylogenetic tree shows that all the species belonging to Cluster A are joined together with a reasonable bootstrap score of 77. This major cluster was further subdivided into two sub-clusters with a high bootstrap value of 98 (Cluster A1; Brachistosternus) and 88 (Cluster A2; Diplocentrus).  The species belonging to the genus Hemiscorpius were grouped together with a good bootstrap score of 96. The Khuzestan’s sample Hl-Kh formed a cluster (Cluster B) relative to the only retrieved sample of Genbank H. Lepturus (Hl) from Iran. Within this cluster, there is a strong support for species Hl and Hl-Kh (96%). The similarity of the species belonging to Pandinus cavimanus (Pc) and Scorpio fuscus (Sf) formed a single cluster with a score of 83 (Cluster C). As it was expected Hl-Kh used in this study was perfectly arranged in cluster B showing higher genetic similarity with Hl and showed higher genetic distance particularly with species in Cluster A.

Genetic distances: The genetic distances between Hemiscorpionidae family scorpions (Cluster A1, A2, B and C) and the outgroup species M. eupeus were calculated (Table 2). The genetic variations of Hl-Kh with cluster A1, A2 and C ranged from 16 to19%. The highest distance was observed with cluster C (18-19%) and the lowest distance was calculated between Hl-Kh and Hl (6%). Finally, as expected, the distances between Hl-Kh and the outgroup species M. eupeus was 25%.

 

 

 

Table 1. The amino acid sequences predicted to form transmemberane helises (Hl-Kh). Only scores above  500 are considered significant.

Inside to outside

Outside to inside

From

To

Score

From

To

Score

29

51

2650

29

51

2153

73

90

1426

75

91

1478

123

144

1448

123

142

1840

163

183

1040

157

178

1037

 

 

Table 2. Genetic distance table of the cytochrome oxidase subunit 1 (Hl-Kh). 

 

1

2

3

4

5

6

7

8

9

10

11

12

13

14

15

1. Hl-Kh

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2. Hl

0.06

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3. Dt

0.16

0.15

 

 

 

 

 

 

 

 

 

 

 

 

 

4. Bf

0.18

0.17

0.18

 

 

 

 

 

 

 

 

 

 

 

 

5. Smf

0.17

0.15

0.16

0.18

 

 

 

 

 

 

 

 

 

 

 

5. Bp

0.18

0.17

0.14

0.13

0.18

 

 

 

 

 

 

 

 

 

 

7. Br

0.19

0.16

0.16

0.14

0.19

0.07

 

 

 

 

 

 

 

 

 

8. Pc

0.19

0.17

0.17

0.18

0.16

0.21

0.22

 

 

 

 

 

 

 

 

9. Bs

0.19

0.18

0.17

0.06

0.19

0.15

0.15

0.19

 

 

 

 

 

 

 

10. Be

0.18

0.16

0.16

0.14

0.18

0.12

0.10

0.19

0.14

 

 

 

 

 

 

11. Db

0.16

0.15

0.11

0.18

0.16

0.14

0.14

0.16

0.17

0.15

 

 

 

 

 

12. Dr

0.18

0.16

0.11

0.18

0.18

0.14

0.15

0.18

0.18

0.14

0.10

 

 

 

 

13.Dc

0.18

0.16

0.12

0.17

0.18

0.14

0.15

0.19

0.17

0.17

0.12

0.05

 

 

 

14. Sf

0.18

0.18

0.18

0.19

0.09

0.17

0.17

0.17

0.19

0.18

0.17

0.17

0.17

 

 

15. Me-Kh

0.25

0.25

0.24

0.25

0.24

0.27

0.25

0.26

0.26

0.22

0.22

0.24

0.24

0.26

 

 

 

 

Discussion

In common with other protein coding genes, COXI appears to have a greater range of phylogenetic signal than any other mitochondrial gene. The third-position nucleotides of this gene show a high frequency rate of base alterations. This makes it about three times more valuable than that of 12S or 16S rDNA in molecular phylogeny (Knowlton and Weigt, 1998).

This study was undertaken with the aim to sequence the mitochondrial Cytochrome C Oxidase subunit 1 for the species H. lepturus. Comparison of 624 bp of the encoded COX I gene provided good resolution of the species that were examined in this study. TMs were found as shown in Fig. 3. The aim of this TM search step was to confirm our sequence database based hypothesis of the putative sequence being identified as a cytochrome C oxidase subunit I. Cythochrome C oxidases are known to have these TM regions and the identification of these regions served as a confirmatory step in determining the validity of the predicted structure (Mohamed et al., 2003). Although COXI gene is somehow conserved among species, both at the level of nucleotide and amino acid sequence, the presence of those TM extents predicted to form amphipathic alpha-helical secondary structures are relatively conserved (Saraste et al., 1994).

The Khuzestan’s Hemiscorpius scorpions in cluster B including Hl-Kh and Hl (isolated from Izeh in Khuzestan) are very divergent from the other Hemiscorpionidae family members (Cluster A1, A2, and C) which are isolated from North Africa, Middle east, South America and Brazil. From the phylogenetic studies we also understand that the Hl-Kh sequence was clustered with the only entry in the Genbank database corresponding to H. lepturus (Hl) with a relatively high distance (6%) among Hemiscorpius genus in Cluster B. Hl-Kh isolated from the Khuzestan province of Iran are closely related to Hl, as they showed the greatest interspecific variation observed in Hemisccorpius genus in this study. Although the sample size in this cluster is not large enough to draw a final conclusion, the percentage of sequence divergence was high enough for interspecific comparisons to provide separation of species. But, further studies need to be conducted to prove this concept by analyzing total mitochodrial DNA sequence.

Acknowledgments

This work was supported by a grant from the Shahid Chamran University of Ahvaz, Ahvaz, Iran.

Conflicts of interest

The author declared no conflict of interest.

Altschul, S. F., Madden, T.L., Schaffer, A.A., Zhang, J., Zhang, Z., Miller, W., Lipman D.J. (1997). Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. Nucleic Acids Res, 25, 3389-3402. https://doi.org/10.1093/nar/25.17.3389 PMID: 9254694
Augot, D., Mathieu, B., Hadj-Henni, L., Barriel, V., Zapata Mena, S., Smolis, S. Slama, D., Randrianambinintsoa, F.J., Trueba, G., Kaltenbach, M., Rahola, N., Depaquit, J. (2017). Molecular phylogeny of 42 species of Culicoides (Diptera, Ceratopogonidae) from three continents. Parasite, 24, 23. https://doi.org/10.1051/parasite/2017020 PMID: 28643630
Bellagamba, F., Velayutham, D., Cristina Cozzi, M., Caprino, F., Vasconi, M., Letizia Busetto, M., Bagnato, A., Moretti, V.M.  (2015). Cytochrome Oxidase-I Sequence Based Studies of Commercially Available Pangasius Hypophthalmus in Italy. J Anim Sci, 14, 3928.  https://doi.org/10.4081/ijas.2015.3928
Fet, V. (2000). Subfamily Hemiscorpiinae Pocock, 1893. In: Catalogue of the Scorpions of the World (1758-1998). Fet, V., Sissom, D. W., LOWE, G., Braunwalder, M. E. (eds.). The New York Entomological Society, New York, USA. p. 428-431.
Folmer, O., Black, M., Hoeh, W., Lutz, R., Vrijenhoek, R. (1994). DNA primers for amplification of mitochondrial cytochrome c oxidase subunit I from diverse metazoan invertebrates. Mol Mar Biol Biotech, 3, 294-299. PMID: 7881515
Garg, R. K., Dubey, N., Batav, N., Pandey, P., Singh R.K. (2017). Mitochondrial COI gene sequence analyses of puntius ticto and compared with seven species of genus puntius of family cyprinidae: A finding for phylogenetic positioning and DNA barcoding as model study for cryptic species. J Proteomics & Bioinformatics, 10(9), 214-221. https://doi.org/10.4172/jpb.1000445
Hebert, P.D.N., Sujeevan, R., Jeremy, R. (2003). Barcoding animal life: cytochrome c oxidase subunit 1 divergences among closely related species. Proc R Soc Lond, 270, 96-99. https://doi.org/10.1098/rsbl.2003.0025 PMID: 12952648
Knowlton, N., Weigt, L.A. (1998). New dates and new rates for divergence across the Isthmus of Panama. Proc R Soc Lond B Biol Sci, 265, 2257-2263. https://doi.org/10.1098/rspb.1998.0568
Linto, Y.M., Mordue Luntz, A.J., Cruickshank, R.H., Meiswinkel, R., Mellor, P. S., Dallas, J.F. (2002). Phylogenetic analysis of the mitochondrial cytochrome oxidase subunit I gene of five species of the Culicoides imicola species complex. Med Vet Entomol, 16(2), 139-46. https://doi.org/10.1046/j.1365-2915.2002.00356.x PMID: 12109707
Manouchehrifar, M., Khosravi, S.H., Khavanin, A., Derakhshandeh, N. (2013). Report of 267 cases of scorpion bite referring to an emergency department during one year. Emergency, 1(1), 24-6. https://doi.org/10.22037/emergency.v1i1.5263 PMID: 26495332
Mohamed, R., Raih, M.R., Sailan, A.T., Zamrod, Z., Embi, M.N. (2003). A predicted structure of the cytochrome c oxidase from Burkholderia Pseudomallei. J Boitec, 6, 17-28. https://doi.org/10.2225/vol6-issue1-fulltext-5
Radmanesh, M. (1990). Clinical study of Hemiscorpion lepturus in Iran. J Trop Med Hyg, 93, 327-332. PMID: 2231841
Radmanesh, M. (1998). Cutaneous manifestations of the Hemiscorpius lepturus sting: a clinical study. Int J Dermatol, 37, 500-507. https://doi.org/10.1046/j.1365-4362.1998.00386.x PMID: 9679690
Saraste, M., Castresana, J., Higgins, D.G., Lubben, M. (1994). Evolution of cytochrome oxidase, an enzyme older than atmospheric oxygen. EMBO J, 13, 2516-2525. https://doi.org/10.1002/j.1460-2075.1994.tb06541.x PMID: 8013452
Sreejith, K., Sebastian, C.D. (2014). Phylogenetic analysis and sequencing of the mitochondrial cytochrome oxidase subunit I (COI) of white backed plant hopper, Sogatella furcifera (Horvath). Int Res J Pharm, 5, 887-890. https://doi.org/10.7897/2230-8407.0512180