2urn:lsid:arphahub.com:pub:45048D35-BB1D-5CE8-9668-537E44BD4C7Eurn:lsid:zoobank.org:pub:91BD42D4-90F1-4B45-9350-EEF175B1727AZooKeysZK1313-29891313-2970Pensoft Publishers10.3897/zookeys.1158.101347101347Research ArticleElapidaeSystematicsAsiaMiddle EastGenetic diversity among sea snakes of the genus Hydrophis (Elapidae, Reptilia) in the Persian Gulf and Gulf of OmanHosseinian YousefkhaniSeyyed Saeedhttps://orcid.org/0000-0002-4871-62061YasserAmaal23NaserMurtada23Rezaie-AtagholipourMohsen4Askari HesniMajid5YousefabadiFariba6Rastegar PouyaniEskandarrastegarpouyani45@gmail.com7Department of Animal Science, School of Biology, Damghan University, Damghan, IranMarine Science Centre, University of Basrah, Basra, IraqAustralian Rivers Institute, Griffith University, 170 Kessels Road, Nathan, Queensland 4111, AustraliaSchool of Environment and Science, Griffith University, 170 Kessels Road, Nathan, Queensland, 4111, AustraliaQeshm Environmental Conservation Institute (QECI), Qeshm Island, Hormozgan Province, IranDepartment of Biology, Faculty of Sciences, Shahid Bahonar University of Kerman, Kerman Province, IranDepartment of Biology, Faculty of Science, Ferdowsi University of Mashhad, Mashhad, IranDepartment of Biology, Hakim Sabzevari University, Sabzevar, Iran
202320042023115812113167BBE695-F872-53BE-94FD-731B8136AAC66E071A9E-7978-4690-B8AE-75C58F44B1D10102202301042023Seyyed Saeed Hosseinian Yousefkhani, Amaal Yasser, Murtada Naser, Mohsen Rezaie-Atagholipour, Majid Askari Hesni, Fariba Yousefabadi, Eskandar Rastegar PouyaniThis is an open access article distributed under the terms of the Creative Commons Attribution License (CC BY 4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.http://zoobank.org/6E071A9E-7978-4690-B8AE-75C58F44B1D1
Sea snakes of the genus Hydrophis are important components of animal diversity in Iranian waters of the Persian Gulf and Gulf of Oman. Ten species of Hydrophis have been identified from the these waters and, in this study, genetic structure of seven species was compared with other populations in the eastern Indian Ocean and the West Pacific. We found that six species (H.platurus, H.cyanocinctus, H.spiralis, H.schistosus, H.gracilis, and H.lapemiodes) show high genetic similarity with conspecific populations in the Indian Ocean and Australia. However, H.curtus from southern Iran shows a high level of genetic differentiation from conspecific populations in Sri Lanka and Indonesia (0.6% and 6% genetic distance from Sri Lankan samples for 16S and COI gene fragments, respectively). Variation between Iranian and Southeast Asian populations may reflect new genetic lineages and suggest the need of further morphological evaluations to re-evaluate their taxonomic position.
DispersalIndian Oceanmicrohabitat adaptationtrue sea snakesvariationIranian National Science foundationCitation
Hosseinian Yousefkhani SS, Yasser A, Naser M, Rezaie-Atagholipour M, Askari Hesni M, Yousefabadi F, Rastegar Pouyani E (2023) Genetic diversity among sea snakes of the genus Hydrophis (Elapidae, Reptilia) in the Persian Gulf and Gulf of Oman. ZooKeys 1158: 121–131. https://doi.org/10.3897/zookeys.1158.101347
Introduction
Sea snakes are a well-known part of marine ecosystems and play an irreplaceable role in the food chain as both predators of small marine fauna and prey for larger predators within their geographical distribution range, which comprises the tropical Indo-Pacific (Voris 1972; Udyawer et al. 2018). True sea snakes of the tribe Hydrophiini share a common ancestor approximately six million years ago, but the major speciation events and radiation within the group and the genus Hydrophis have taken place over the last 3.5 million years (Sanders et al. 2013). While studies on different subjects (e.g., Ukuwela et al. 2012, 2016) have been conducted on sea snakes on the eastern half of their geographic range (i.e., Southeast Asia, India, Indonesia, and northern Australia), their populations remain relatively understudied in the western Indian Ocean (Rezaie-Atagholipour et al. 2016).
Iranian coastal waters of the Persian Gulf and Gulf of Oman are known for their high biological diversity (Owfi et al. 2016) and share a marine herpetofauna with the Indian Ocean (Gasperetti 1988; Heatwole 1999). Ten species of Hydrophis have been recorded in these gulfs: H.platurus, H.schistosus, H.curtus, H.viperinus, H.spiralis, H.cyanocinctus, H.ornatus, H.lapemiodes, H.gracilis, and H.cantoris (Rezaie-Atagholipour et al. 2016; Rajabizadeh 2019). Among these, H.platurus has a wider geographic range toward eastern Africa (Heatwole 1999), but the westernmost ranges of the other species extend to the Persian Gulf (Rezaie-Atagholipour et al. 2016). The semi-enclosed basin of the Persian Gulf (i.e., connected to the Gulf of Oman and western Indian Ocean via the narrow Strait of Hormuz) and its unique extreme environment with high salinity and temperatures (i.e., sea surface temperatures can exceed 34 °C and salinity more than 39 ppt in most areas; Sheppard et al. 2010) may pose ecological and geographical barriers and limit gene flow between the Gulf and Indian Ocean marine communities (e.g., Natoli et al. 2017). Therefore, it seems important to compare the genetic structure of sea snakes, which are known to have undergone rapid speciation (see Sanders et al. 2013), between the Persian Gulf and the rest of the Indo-Pacific.
In this study, available genetic markers (16S rRNA, COI, and G1888) of species of sea snakes of the genus Hydrophis from the Persian Gulf and the Gulf of Oman were compared to populations in other areas of the Indian Ocean, Southeast Asia, and Indonesia. Examining the genetic structure of the species in the Persian Gulf and comparing them with other populations can illustrate the degree of genetic connection between these populations.
Materials and methodsTissue sampling and DNA extraction
Tissue samples were obtained from the specimens collected in the southern coastal regions of Iran (Persian Gulf and Gulf of Oman) during fieldwork in 2013 (Rezaie-Atagholipour et al. 2016) and were preserved in 96% ethanol. A total of 31 samples used in this study belong to seven species from four localities in southern Iran (Fig. 1, Suppl. material 1) and 27 comparative sequences were downloaded from GenBank. Sequences of Ephalophisgreyae, used as the outgroup, were also downloaded from GenBank (Suppl. material 1). Total genomic DNA was extracted from muscle tissue using the standard proteinase K – salt method (Kabir et al. 2006), and the quality and concentration of extracted DNA were measured using Nanodrop 1000.
A map of the world and the selected region as study area B localities in South and Southeast Asia from which samples were used in this study C localities in southern Iran where samples were collected: 1 = Bushehr; 2 = Larak Island; 3 = Jask; 4 = Beris and Pasabandar.
https://binary.pensoft.net/fig/840536Mitochondrial and nuclear fragment sequences
Two mitochondrial sequences, 16S rRNA (16S) (Kocher et al. 1989) and cytochrome oxidase subunit 1 (COI), and one anonymous nuclear marker (G1888) (Bertozzi et al. 2012) were used to reconstruct the molecular phylogenetic relationships between Iranian and Southeast Asian species of sea snakes. Two of these markers were obtained from previous successful studies (Lukoschek and Keogh 2006; Sanders et al. 2013), but COI is a new marker generated for sea snakes in the present study. The following primers were used for each gene fragment: 16S: 16SL 5'-CGCCTGTTT ATCAAAAACAT-3'/ 16SH 5'-CCGGTCTGAACTCAGATCACG-3'; COI: RepCOIF 5'-TNTTMTCAACNAACCACAAAGA-3'/ RepCOIR 5'-ACTTCTGGRTGKCCAAARAATCA-3'; G1888: G1888F 5'- CAGGGCCTTGCCTTGTGCCA-3'/ G1888R 5'-ACCTCTGCGCACTATGACTCTTGA-3' (Bertozzi et al. 2012). All DNA sequences were amplified using a standard PCR protocol with denaturation at 94 °C for 5 min, annealing temperatures of 49 °C for mitochondrial fragments and 52 °C for anonymous nuclear marker for 45 s, and extension at 72 °C for 70 s (36 cycles), and the final extension for 8 min. Sequencing of the PCR products was performed by the Kodon Genetic Group in Tehran, Iran. Sequences of the same species and the same genetic fragments from Southeast Asia and the outgroup were downloaded from GenBank and added to the dataset (Suppl. material 1). All sequences were aligned using the ClustalW algorithm implemented within Bioedit v. 7.0.9.0 (Hall 1999), and the protein coding genes were translated to amino acid sequences using MEGA v. 6.0 (Tamura et al. 2013) to check for internal stop codons and to determine the correct reading frame. The sequences generated in this study were deposited in GenBank and will be added into Suppl. material 1. Uncorrected genetic distances (p-distances) were calculated (for 16S and COI) and number of variable sites (V) and parsimony informative sites (Pi) were obtained using MEGA v. 6.0 for all markers. The sequences of the COI gene fragment were few and we only calculate genetic distance among some clades that revealed relatively high values.
Molecular phylogenetic analyses and haplotype network
The one mitochondrial and the nuclear fragment were concatenated (total length: 1015 bp; 16S: 523 bp; G1888: 492 bp) and COI gene fragment (708 bp) separately were considered to reconstruction phylogenetic trees using Bayesian inference (BI) and maximum-likelihood (ML) methods. The best substitution models for each alignment were: TrNef+I, F81+I, TIM+G, GTR+I+G and TVM+I+G for the first, second, and third codon positions of COI and 16S and G1888, respectively. RaxML v. 7.4.2 (Stamatakis 2006) as implemented in RaxmlGUI v. 1.3 (Silvestro and Michalak 2012) was employed for the ML analyses. The substitution model was set to GTR+G+I in RaxMl software for the molecular phylogenetic reconstruction. Bootstrapping was set to 1000 replicates to estimate nodal support, and the analysis was run as a heuristic search (Felsenstein 1985). MrBayes v. 3.2.1 (Ronquist et al. 2012) was used to run the BI analysis and the number of generations were set as 107 with a sample frequency of every 1000 generations, four Markov chains, and a burn-in with each 1000 samples.
The relationships among lineages and subclades of Hydrophis were assessed with the mitochondrial 16S rRNA gene, because the number of sequences of 16S covered all species. Aligned sequences were entered in DNAsp v. 5.0. (Librado and Rozas 2009) and a *.rdf extension file was created. The haplotype network was then constructed using the median-joining method in Network 1.2.1.
Results
Our dataset includes 59 samples containing two mitochondrial gene fragments 16S (495 bp; 77 V; 40 Pi) and COI (661 bp; 158 V; 112 Pi) and one nuclear anonymous fragment G1888 (375 bp; 265 V; 250 Pi) totaling 1531 bp. Both the ML and BI trees show similar topology and therefore we present only the BI tree (Fig. 2). Our results show that all species of Hydrophis in the Persian Gulf and Gulf of Oman similar with their corresponding species in Southeast Asia except H.curtus and H.lapemiodes. Hydrophiscurtus in the western part of its geographic range (Persian Gulf and Gulf of Oman) shows clear divergence from Indian Ocean (Sri Lanka) and Indonesian samples (Ukuwela et al. 2014). Hydrophisornatus also shows minor variation between Iranian and Indian Ocean populations, but not as extensive.
Molecular phylogenetic tree of the genus Hydrophis. Red and green clades indicate distinct populations of H.curtus and H.ornatus from the Persian Gulf relative to other populations in the Indian Ocean. Data from samples with an asterisk (*) were downloaded from GenBank.
https://binary.pensoft.net/fig/840537
Tree topology was used to group the sequences and calculate the uncorrected genetic distance (p-distance) among the lineages (Table 1). Mean genetic distance in 16S rRNA is relatively low (average 3.25% among all lineages) in Hydrophis. Iranian lineages of H.curtus are differentiated from those in Southeast Asia by 0.6% and 6% in 16S and COI gene fragments, respectively (Tables 1, 2). Two other distinct species, H.platurus and H.viperinus, show the genetic differentiation of 0.6% in 16S and 5.1% in COI (Tables 1, 2). The haplotype network calculated for the 16S rRNA fragment confirmed the relationship among lineages in the concatenated phylogenetic tree (Fig. 3). The Iranian clade of H.curtus had a separate haplotype (haplotype number 2), but the other conspecific samples (Pacific Ocean) had a different distinct haplotype (haplotype number 16).
Haplotype network of the genus Hydrophis for the 16S gene fragment. The haplotype colors correspond with the those in the molecular phylogenetic tree in Fig. 2.
https://binary.pensoft.net/fig/840538Discussion
Sea snakes are one of the most interesting reptiles in biogeographic and phylogeographic studies (Ukuwela et al. 2022). A biogeographical study of the genus Hydrophis indicated that the diversification of true sea snakes has been the result of sea level changes during the last 2.5 million years (Ukuwela et al. 2016). Our molecular results indicate that the sea snakes in the Persian Gulf and Gulf of Oman share a genetic structure with Indian Ocean and Southeast Asian taxa (Fig. 2). Among them, H.curtus has differentiated relatively more than others, so that it may represent a genetically distinct species. Ukuwela et al. (2014) had the same findings for H.curtus. Genetic variation among Persian Gulf populations of H.curtus and other populations in southeast Asia (0.6% for 16S and 6% for COI) indicate that the former may representative a new genetic lineage in Southwest Asia (Tables 1, 2). This variation among this taxon indicates that the threshold of genetic variation may reveal distinct populations at the species level, although this variation has not yet been established on the basis of morphology (Rezaie-Atagholipour et al. 2016).
A recent study on the genetic structure of H.curtus in the Indo-West Pacific revealed the presence of high variation among its populations. The phylogeny indicated that this species was affected by climatic fluctuations in late Pliocene and early Pleistocene (Ukuwela et al. 2014, 2022). Adding the Iranian dataset confirms that H.curtus populations in the Persian Gulf in the western part of its range have a close relationship with those in Sri Lanka and exhibit genetic variability greater than those populations in Southeast Asia (Fig. 2). Other species of Hydrophis have not yet been studied, but our study indicates that their differentiation requires conducting a comprehensive study. However, the variation among populations of these species implies a pattern of variation in the Indian Ocean and West Pacific regions that may occur either in Sri Lanka or in the Persian Gulf as well (Lambeck 1996; Voris 2000; Lambeck et al. 2002).
Other species of the genus, H.platurus, H.cyanocinctus, H.spiralis, H.ornatus, H.schistosus and H.gracilis, included in the analyses, clearly placed within their conspecific Southeast Asian clades, and this may indicate high gene flow and recent (i.e. 10,000 years) dispersal from their source populations (Fig. 2). We assume that different local adaptations among H.curtus and the other six other taxa make them more variable.
The total diversity of the true sea snakes of the genus Hydrophis in the Persian Gulf and Gulf of Oman is 10 species (Rezaie-Atagholipour et al. 2016). According to the literature (Ukuwela et al. 2014; Rezaie-Atagholipour et al. 2016; Rajabizadeh 2019), H.schistosus, H.viperinus, H.lapemiodes, H.cantoris, H.gracilis, H.spiralis, and H.platurus are abundant in Gulf of Oman as opposed to the Persian Gulf, which means the species prefer to not inhabit the Gulf. But H.curtus and H.ornatus are mostly distributed in the Persian Gulf and in the Gulf of Oman. Our results revealed the variation between Persian Gulf and Southeast Asia population of H.curtus (6% genetic distance in the COI gene fragment) and reveals a new distinct genetic lineage. The genetic difference between DNA sequences of H.curtus populations in the Persian Gulf and Pacific Ocean is very likely due to genetic drift. In general, biodiversity conservation requires information (demographic and microhabitat information) from local populations. Only then we can encourage conservationists to reassess the status of the species in the Persian Gulf and the Gulf of Oman, because the latest International Union for the Conservation of Nature (IUCN) status of most species is Least Concern and one is Data Deficient (IUCN 2022).
Acknowledgements
We thank L. Lee Grismer, Ann Paterson, Aaron Bauer, and Kamelia Algiers who edited the early draft of the manuscript. This study was partially supported by Iran National Science Foundation (INSF) with project number of 99025967. Tissue samples were collected during field surveys financially supported by the Iranian UNDP/GEF/SGP (grant number: IRA/SGP/OP5/Y3/STAR/BD/13/12-188).
ReferencesBertozziTSandersKLSistromMJGardnerMG (2012) Anonymous nuclear loci in non-model organisms: Making the most of high-throughput genome surveys.28(14): 1807–1810. https://doi.org/10.1093/bioinformatics/bts284FelsensteinJ (1985) Confidence limits on phylogenies: An approach using the bootstrap.39(4): 783–791. https://doi.org/10.2307/2408678GasperettiJ (1988) The snakes of Saudi Arabia. In: BüttikerWKruppF (Eds) Fauna of Saudi Arabia (Vol., 169–450.HallTA (1999) BioEdit: A user-friendly biological sequence alignment editor and analysis program for Windows 95/98/NT.41: 95–98.HeatwoleH (1999) Sea Snakes. University of New South Wales Press, Sydney.IUCN (2022) The IUCN Red List of Threatened Species. Version 2022-2. [Accessed on 2023-01-10]KabirSShahriarMKabirAHUddinMG (2006) High salt SDS-based method for the direct extraction of genomic DNA from three different gram-negative organisms.1: 57–64.KocherTDThomasWKMeyerAEdwardsSVPääboSVillablancaFXWilsonAC (1989) Dynamics of mitochondrial DNA evolution in animals: Amplification and sequencing with conserved primers.86(16): 6196–6200. https://doi.org/10.1073/pnas.86.16.6196LambeckK (1996) Shoreline reconstructions for the Persian Gulf since the last glacial maximum.142(1–2): 43–57. https://doi.org/10.1016/0012-821X(96)00069-6LambeckKEsatTMPotterEK (2002) Links between climate and sea levels for the past three million years.419(6903): 199–206. https://doi.org/10.1038/nature01089LibradoPRozasJ (2009) DnaSP v5: A software for comprehensive analysis of DNA polymorphism data.25(11): 1451–1452. https://doi.org/10.1093/bioinformatics/btp187LukoschekVKeoghJS (2006) Molecular phylogeny of sea snakes reveals a rapidly diverged adaptive radiation. Biological Journal of the Linnean Society.89(3): 523–539. https://doi.org/10.1111/j.1095-8312.2006.00691.xNatoliAPhillipsKPRichardsonDSJabadoRW (2017) Low genetic diversity after a bottleneck in a population of a critically endangered migratory marine turtle species.491: 9–18. https://doi.org/10.1016/j.jembe.2017.01.009OwfiFBraulikGTRabbanihaM (2016) Species diversity and distribution pattern of marine mammals of the Persian Gulf and Gulf of Oman-Iranian waters.15: 927–944.RajabizadehM (2019) The Snakes of Iran. Iranshenasi Press, Tehran.Rezaie-AtagholipourMGhezellouPHesniMADakhtehSMHAhmadianHVidalN (2016) Sea snakes (Elapidae, Hydrophiinae) in their westernmost extent: An updated and illustrated checklist and key to the species in the Persian Gulf and Gulf of Oman.622: 129–164. https://doi.org/10.3897/zookeys.622.9939RonquistFTeslenkoMVan Der MarkPAyresDLDarlingAHöhnaSHuelsenbeckJP (2012) MrBayes 3.2: Efficient Bayesian phylogenetic inference and model choice across a large model space.61(3): 539–542. https://doi.org/10.1093/sysbio/sys029SandersKLLeeMSBertozziTRasmussenAR (2013) Multilocus phylogeny and recent rapid radiation of the viviparous sea snakes (Elapidae: Hydrophiinae).66(3): 575–591. https://doi.org/10.1016/j.ympev.2012.09.021SheppardCRCAl-HusianiMAl-JamaliFAl-YamaniFBaldwinRBishopJBenzoniFDutrieuxEDulvyNKDurvasulaSRVJonesDALoughlandRMedioDNithyanandanMPillingGMPolikarpovIPriceARGPurkisSRieglBSaburovaMNaminKSTaylorOWilsonSZainalK (2010) The Gulf: A young sea in decline.60(1): 13–38. https://doi.org/10.1016/j.marpolbul.2009.10.017SilvestroDMichalakI (2012) raxmlGUI: A graphical front-end for RAxML.12(4): 335–337. https://doi.org/10.1007/s13127-011-0056-0StamatakisA (2006) RAxML-VI-HPC: Maximum likelihood-based phylogenetic analyses with thousands of taxa and mixed models.22(21): 2688–2690. https://doi.org/10.1093/bioinformatics/btl446TamuraKStecherGPetersonDFilipskiAKumarS (2013) MEGA6: Molecular evolutionary genetics analysis version 6.0.30(12): 2725–2729. https://doi.org/10.1093/molbev/mst197UdyawerVBarnesPBonnetXBrischouxFCrowe-RiddellJMD’anastasiBFryBGGillettAGoiranCGuineaMLHeatwoleHHeupelMRHourstonMKangasMKendrickAKoefoedILillywhiteHBLoboASLukoschekWMcAuleyRNitschkeCRasmussenARSandersKLSheehyC IIIShineRSomaweeraRSweetSSVorisHK (2018) Future directions in the research and management of marine snakes. Frontiers in Marine Science 5: e00399. https://doi.org/10.3389/fmars.2018.00399UkuwelaKDSandersKLFryBG (2012) Hydrophisdonaldi (Elapidae, Hydrophiinae), a highly distinctive new species of sea snake from northern Australia.3201(1): 45–57. https://doi.org/10.11646/zootaxa.3201.1.3UkuwelaKDde SilvaAFryBGSandersKL (2014) Multilocus phylogeography of the sea snake Hydrophiscurtus reveals historical vicariance and cryptic lineage diversity.43(5): 472–484. https://doi.org/10.1111/zsc.12070UkuwelaKDLeeMSRasmussenARDe SilvaAFryBGGhezellouPSandersKL (2016) Evaluating the drivers of Indo‐Pacific biodiversity: speciation and dispersal of sea snakes (Elapidae: Hydrophiinae).43(2): 243–255. https://doi.org/10.1111/jbi.12636UkuwelaKDBde SilvaASivarubanASandersKL (2022) Diversity, distribution, and natural history of the marine snakes of Sri Lanka.52(2): 1–24. https://doi.org/10.1007/s12526-022-01259-3VorisHK (1972) The role of sea snakes (Hydrophiidae) in the trophic structure of coastal ocean communities.14: 429–442.VorisHK (2000) Maps of Pleistocene sea levels in Southeast Asia: Shorelines, river systems and time durations.27(5): 1153–1167. https://doi.org/10.1046/j.1365-2699.2000.00489.xSupplementary materials10.3897/zookeys.1158.101347.suppl1EB827EC3-2869-5CA5-A9AA-70BB78097688
Information of sequences used in this study and those obtained from GenBank
occurences, accession numbers
Information of sequences used in this study and those obtained from GenBank (sequences with an asterisk, *). Accession numbers for the sequences added to the table. Accession numbers in red related to the new sequences of this study.
https://binary.pensoft.net/file/840539This dataset is made available under the Open Database License (http://opendatacommons.org/licenses/odbl/1.0/). The Open Database License (ODbL) is a license agreement intended to allow users to freely share, modify, and use this Dataset while maintaining this same freedom for others, provided that the original source and author(s) are credited.Seyyed Saeed Hosseinian Yousefkhani, Amaal Yasser, Murtada Naser, Mohsen Rezaie-Atagholipour, Majid Askari Hesni, Fariba Yousefabadi, Eskandar Rastegar Pouyani10.3897/zookeys.1158.101347.suppl27DE08688-CD4A-5F7C-A2F7-C4F180352FBD
Bayesian tree of COI gene fragment that clearly shows variation in Hydrophiscurtus clade
figure
https://binary.pensoft.net/file/840540This dataset is made available under the Open Database License (http://opendatacommons.org/licenses/odbl/1.0/). The Open Database License (ODbL) is a license agreement intended to allow users to freely share, modify, and use this Dataset while maintaining this same freedom for others, provided that the original source and author(s) are credited.Seyyed Saeed Hosseinian Yousefkhani, Amaal Yasser, Murtada Naser, Mohsen Rezaie-Atagholipour, Majid Askari Hesni, Fariba Yousefabadi, Eskandar Rastegar Pouyani