Abstract
The Lesser Sunda Seascape (LSS) is a critical habitat for various megafauna, including the spinetail devil ray (Mobula mobular). For decades, bycatch from drift gillnet fisheries has increased in the LSS, posing a greater threat to their populations. Migratory species require population-based management to effectively address their conservation needs and ensure their long-term survival. Fifty-five samples were collected from three locations within the LSS, i.e., East Flores (19), West Manggarai (15), and East Java (21). This study aims to understand the population genetic structure of M. mobular within the three different locations around the LSS using eight microsatellite markers and mitochondrial DNA (mtDNA) NADH-dehydrogenase subunit 5. The analysis of molecular variances result indicated the significant result between the three locations based on microsatellite data (0.0876, p < 0.05 for eight loci and 0.1360, p < 0.05 for four loci) and mtDNA (0.04820, p < 0.05). The pairwise genetic differentiation index (Fst) revealed that East Flores versus West Manggarai has a nonsignificant value (p > 0.05). Meanwhile, East Flores or West Manggarai versus East Java showed a significant Fst value (p < 0.05). In conclusion, microsatellite and mtDNA results revealed that in the LSS, spinetail devil rays are divided into two distinct subpopulations: the Savu Sea (East Flores and West Manggarai) subpopulation and the East Java subpopulation. Even though there are two potential subpopulations in the LSS, co-management of this species should be prioritized in the LSS due to the indication of shared haplotypes between those locations.
Similar content being viewed by others
Data availability
The raw mitochondrial DNA for ND5 gen of this study has been submitted to National Center for Biotechnology Information (NCBI) GenBank database with accession numbers OR248469–OR248523. Then, the raw microsatellite DNA score result of this study is available from the corresponding author upon reasonable request.
References
Armstrong AJ, Dudgeon CL, Bustamante C, Bennett MB, Ovenden JR (2019) Development and characterization of 17 polymorphic microsatellite markers for the reef manta ray (Mobula alfredi). BMC Res Notes 12:1–5. https://doi.org/10.1186/s13104-019-4270-8
Astanei I, Gosling E, Wilson JIM, Powell E (2005) Genetic variability and phylogeography of the invasive zebra mussel, Dreissena polymorpha (Pallas). Mol Ecol 14(6):1655–1666. https://doi.org/10.1111/j.1365-294X.2005.02530.x
Ayre DJ, Hughes TP (2000) Genotypic diversity and gene flow in brooding and spawning corals along the Great Barrier Reef, Australia. Evoluation 54(5):1590–1605. https://doi.org/10.1111/j.0014-3820.2000.tb00704.x
Bernard AM, Feldheim KA, Heithaus MR, Wintner SP, Wetherbee BM, Shivji MS (2016) Global population genetic dynamics of a highly migratory, apex predator shark. Mol Ecol 25(21):5312–5329. https://doi.org/10.1111/mec.13845
Billington N (2003) Mitochondrial DNA. In: Hallerman EM (ed) Population genetics: principle and application for fisheries scientists. American Fisheries Society, Bethesda, pp 59–100
Bograd SJ, Block BA, Costa DP, Godley BJ (2010) Biologging technologies: new tools for conservation. Introduction. Endanger Species Res 10:1–7. https://doi.org/10.3354/esr00269
Canese S, Cardinali A, Romeo T, Giusti M, Salvati E, Angiolillo M, Greco S (2011) Diving behavior of the giant devil ray in the Mediterranean Sea. Endanger Species Res 14(2):171–176. https://doi.org/10.3354/esr00349
Chapuis MP, Estoup A (2007) Microsatellite null alleles and estimation of population differentiation. Mol Biol Evol 24(3):621–631. https://doi.org/10.1093/molbev/msl191
Couturier LIE, Marshall AD, Jaine FRA, Kashiwagi T, Pierce SJ, Townsend KA, Weeks SJ, Bennett MB, Richardson AJ (2012) Biology, ecology and conservation of the Mobulidae. J Fish Biol 80(5):1075–1119. https://doi.org/10.1111/j.1095-8649.2012.03264.x
Croll DA, Dewar H, Dulvy NK, Fernando D, Francis MP, Galván-Magaña F, Hall M, Heinrichs S, Marshall A, McCauley D, Newton KM, Notarbartolo-Di-Sciara G, O’Malley M, O’Sullivan J, Poortvliet M, Roman M, Stevens G, Tershy BR, White WT (2016) Vulnerabilities and fisheries impacts: the uncertain future of manta and devil rays. Aquat Conserv 26(3):562–575. https://doi.org/10.1002/aqc.2591
Dakin EE, Avise JC (2004) Microsatellite null alleles in parentage analysis. Heredity 93(5):504–509. https://doi.org/10.1038/sj.hdy.6800545
DeBoer TS, Naguit MRA, Erdmann MV, Ablan-Lagman MCA, Carpenter KE, Toha AHA, Barber PH (2014) Concordant phylogenetic patterns inferred from mitochondrial and microsatellite DNA in the giant clam Tridacna crocea. Bull Mar Sci 90(1):301–329. https://doi.org/10.5343/bms.2013.1002
Dewar H (2002) Preliminary report: Manta Harvest in Lamakera. Report from the PflegerInstitue of Environmental Research and the Nature Conservancy, P. Oceanside, USA, vol 3
Dharmadi, Fahmi, Satria F (2015) Fisheries management and conservation of sharks in Indonesia. Afr J Mar Sci 37(2):249–258. https://doi.org/10.2989/1814232X.2015.1045431
di Sciara GN, Stevens G, Fernando D (2020) The giant devil ray Mobula mobular (Bonnaterre, 1788) is not giant, but it is the only spinetail devil ray. Mar Biodivers Rec 13(1):1–5. https://doi.org/10.1186/s41200-020-00187-0
Dulvy NK, Sadovy Y, Reynolds JD (2003) Extinction vulnerability in marine populations. Fish Fish 4(1):25–64. https://doi.org/10.1046/j.1467-2979.2003.00105.x
Evanno G, Regnaut S, Goudet J (2005) Detecting the number of clusters of individuals using the software structure: a simulation study. Mol Ecol 14(8):2611–2620. https://doi.org/10.1111/j.1365-294X.2005.02553.x
Excoffier L, Lischer HE (2010) Arlequin suite ver 3.5 a new series of programs to perform population genetics analyses under Linux and Windows. Mol Ecol Resour 10(3):564–567. https://doi.org/10.1111/j.1755-0998.2010.02847.x
Francis MP, Jones EG (2016) Movement, depth distribution and survival of spinetail devilrays (Mobula japanica) tagged and released from purse-seine catches in New Zealand. Aquat Conserv Mar Freshw Ecosyst 27:219–236. https://doi.org/10.1002/aqc.2641
Garza JC, Williamson EG (2001) Detection of reduction in population size using data from microsatellite loci. Mol Ecol 10(2):305–318. https://doi.org/10.1046/j.1365-294X.2001.01190.x
Gillett R (2011) Bycatch in small-scale tuna fisheries: a global study. FAO Fisheries and Aquaculture Technical Paper, (560), I
Gilmour J (2002) Substantial asexual recruitment of mushroom corals contributes little to population genetics of adults in conditions of chronic sedimentation. Mar Ecol Prog Ser 235:81–91. https://doi.org/10.3354/meps235081
Hale ML, Burg TM, Steeves TE (2012) Sampling for microsatellite-based population genetic studies: 25 to 30 individuals per population is enough to accurately estimate allele frequencies. PLoS One 7:e45170. https://doi.org/10.1371/journal.pone.0045170
Hays GC, Mortimer JA, Ierodiaconou D, Esteban N (2014) Use of long-distance migration patterns of an endangered species to inform conservation planning for the world’s largest marine protected area. Conserv Biol 28(6):1636–1644. https://doi.org/10.1111/cobi.12325
Hedgecock D, Li G, Hubert S, Bucklin K, Ribes V (2004) Widespread null alleles and poor cross-species amplification of microsatellite DNA loci cloned from the Pacific oyster, Crassostrea gigas. J Shellfish Res 23:379–386
Hoffman J, Dasmahapatra K, Amos W, Phillips C, Gelatt TS, Bickham JW (2009) Contrasting patterns of genetic diversity at three different genetic markers in a marine mammal metapopulation. Mol Ecol 18:2961–2978. https://doi.org/10.1111/j.1365-294X.2009.04246.x
Hui M, Nuryanto A, Kochzius M (2017) Concordance of microsatellite and mitochondrial DNA markers in detecting genetic population structure in the boring giant clam Tridacna crocea across the Indo-Malay Archipelago. Mar Ecol 38(1):e12389. https://doi.org/10.1111/maec.12389
Jombart T (2008) ADEGENET: a R package for the multivariate analysis of genetic markers. Bioinformatics 24(11):1403–1405. https://doi.org/10.1093/bioinformatics/btn129
Kashiwagi T, Marshall AD, Bennett MB, Ovenden JR (2012) The genetic signature of recent speciation in manta rays (Manta alfredi and M. birostris). Mol Phylogenet Evol 64(1):212–218. https://doi.org/10.1016/j.ympev.2012.03.020
Kimura M (1980) A simple method for estimating evolutionary rates of base substitutions through comparative studies of nucleotide sequences. J Mol Evol 16:111–120. https://doi.org/10.1007/BF01731581
Kopelman NM, Mayzel J, Jakobsson M, Rosenberg NA, Mayrose I (2015) CLUMPAK: a program for identifying clustering modes and packaging population structure inferences across K. Mol Ecol Resour 15(5):1179–1191. https://doi.org/10.1111/1755-0998.12387
Krützen M, Sherwin WB, Berggren P, Gales N (2004) Population structure in an inshore cetacean revealed by microsatellite and mtDNA analysis: bottlenose dolphins (Tursiops sp.) in Shark Bay, Western Australia. Mar Mamm Sci 20(1):28–47. https://doi.org/10.1111/j.1748-7692.2004.tb01139.x
Kumar S, Stecher G, Li M, Knyaz C, Tamura K (2018) MEGA X: molecular evolutionary genetics analysis across computing platforms. Mol Biol 35(6):1547–1549. https://doi.org/10.1093/molbev/msy096
Lassauce H, Dudgeon CL, Armstrong AJ, Wantiez L, Carroll EL (2022) Evidence of fine-scale genetic structure for reef manta rays Mobula alfredi in New Caledonia. Endanger Species Res 47:249–264. https://doi.org/10.3354/esr01178
Launey S, Hedgecock D (2001) High genetic load in the Pacific oyster Crassostrea gigas. Genetics 159:255–265. https://doi.org/10.1093/genetics/159.1.255
Lawson JM, Fordham SV, O’Malley MP, Davidson LNK, Walls RHL, Heupel MR, Stevens G, Fernando D, Budziak A, Simpfendorfer CA, Ender I, Francis MP, Notarbartolo di Sciara G, Dulvy NK (2017) Sympathy for the devil: a conservation strategy for devil and manta rays. PeerJ 5:e3027. https://doi.org/10.7717/peerj.3027
Leigh JW, Bryant D (2015) POPART: full-feature software for haplotype network construction. Methods Ecol Evol 6(9):1110–1116. https://doi.org/10.1111/2041-210X.12410
Lewis SA, Setiasih N, Fahmi F, Dharmadi D, O’Malley MP, Campbell SJ, Sianipar AB (2015) Assessing Indonesian manta and devil ray populations through historical landings and fishing community interviews. PeerJ Prepr. https://doi.org/10.7287/peerj.preprints.1334v1
Lezama-Ochoa N, Hall MA, Pennino MG, Stewart JD, López J, Murua H (2019) Environmental characteristics associated with the presence of the Spinetail devil ray (Mobula mobular) in the eastern tropical Pacific. PLoS ONE 14(8):e0220854. https://doi.org/10.1371/journal.pone.0220854
Li G, Hubert S, Bucklin K, Ribes V, Hedgecock D (2003) Characterization of 79 microsatellite DNA markers in the Pacific oyster Crassostrea gigas. Mol Ecol Notes 3:228–232. https://doi.org/10.1046/j.1471-8286.2003.00406.x
Mackenzie JB, Munday PL, Willis BL, Miller DJ, van Oppen MJH (2004) Unexpected patterns of genetic structuring among locations but not colour morphs in Acropora nasuta (Cnidaria; Scleractinia) ). Mol Ecol 13:9–20. https://doi.org/10.1046/j.1365-294X.2003.02019.x
Maier E, Tollrain R, Rinkevich B, Nurnberger B (2005) Isolation by distance in the scleractinian coral Seriatopora hystrix from the red sea. Mar Biol 147:1109–1120. https://doi.org/10.1007/s00227-005-0013-6
Makhawi AM, Liu XB, Yang SR, Liu QY (2013) Genetic variations of ND5 gene of mtDNA in populations of Anopheles sinensis (Diptera: Culicidae) malaria vector in China. Parasites Vectors 6:1–11. https://doi.org/10.1186/1756-3305-6-290
Malik MDA, Putra MIH, Topan E, Pertiwi NPD, Artiningsih EY, Sari SK, Lewis S, Prabuning D, Sembiring A (2022) Genetic variation of oceanic manta ray (Mobula birostris) based on mtDNA data in the Savu Sea, Indonesia. Biodiversitas. https://doi.org/10.13057/biodiv/d230362
Marshall A, Barreto R, Carlson J, Fernando D, Fordham S, Francis MP, Herman K, Jabado RW, Liu KM, Rigby CL, Romanov E (2022) Mobula mobular (amended version of 2020 assessment). The IUCN Red List of Threatened Species 2022: e.T110847130A214381504. https://doi.org/10.2305/IUCN.UK.2022-1.RLTS.T110847130A214381504.en. Accessed 18 June 2023
Masangcay SIG, Metillo EB, Hayashizaki KI, Tamada S, Nishida S (2018) Feeding habits of Mobula japanica (Chondrichthyes, Mobulidae) in Butuan Bay, Mindanao Island, Philippines. Sci Diliman 30(1):24–44
Mayer C, Schiegg K, Pasinelli G (2009) Patchy population structure in a short-distance migrant: evidence from genetic and demographic data. Mol Ecol 18(11):2353–2364. https://doi.org/10.1111/j.1365-294X.2009.04200.x
Meglecz E, Petenian F, Danchin E, D’Acier AC, Rasplus JY, Faure E (2004) High similarity between flanking regions of different microsatellites detected within each of two species of Lepidoptera: Parnassius apollo and Euphydryas aurinia. Mol Ecol 13(6):1693–1700. https://doi.org/10.1111/j.1365-294X.2004.02163.x
Meirmans PG (2006) Using the AMOVA framework to estimate a standardized genetic differentiation measure. Evolution 60(11):2399–2402. https://doi.org/10.1111/j.0014-3820.2006.tb01874.x
Musick JA, Ellis JK, Hamlett W (2005) Reproductive evolution of chondrichthyans. In: Reproductive biology and phylogeny of chondrichthyes: sharks, batoids and chimaeras, vol 3, pp 45–80
Nishikawa A, Sakai K (2005) Genetic connectivity of the scleractinian coral Goniastrea aspera around Okinawa Islands. Coral Reefs 24:318–323. https://doi.org/10.1007/s00338-005-0484-4
Peakall ROD, Smouse PE (2006) GENALEX 6: genetic analysis in Excel. Population genetic software for teaching and research. Mol Ecol Notes 6(1):288–295. https://doi.org/10.1111/j.1471-8286.2005.01155.x
Pichler FB, Baker CS (2000) Loss of genetic diversity in the endemic Hector’s dolphin due to fisheries-related mortality. Proc R Soc B: Biol Sci 267(1438):97–102. https://doi.org/10.1098/rspb.2000.0972
Piry S, Luikart G, Cornuet JM (1999) Computer note. BOTTLENECK: a computer program for detecting recent reductions in the effective population size using allele frequency data. Heredity 90(4):502–503. https://doi.org/10.1093/jhered/90.4.502
Poortvliet M, Galván-Magana F, Bernardi G, Croll DA, Olsen JL (2011) Isolation and characterization of twelve microsatellite loci for the Japanese Devilray (Mobula japanica). Conserv Genet Resour 3(4):733–735. https://doi.org/10.1007/s12686-011-9445-0
Pritchard JK, Stephens M, Donnelly P (2000) Inference of population structure using multilocus genotype data. Genetics 155:945–959. https://doi.org/10.1093/genetics/155.2.945
Pujolar JM, Schiavina M, Di Franco A, Melià P, Guidetti P, Gatto M, De Leo GA, Zane L (2013) Understanding the effectiveness of marine protected areas using genetic connectivity patterns and Lagrangian simulations. Divers Distrib 19(12):1531–1542. https://doi.org/10.1111/ddi.12114
Putra MIH, Mustika PLK (2020) Incorporating in situ prey distribution into foraging habitat modelling for marine megafauna in the Solor waters of the Savu Sea, Indonesia. Aquat Conserv 30(12):2384–2401. https://doi.org/10.1002/aqc.3379
Putra MIH, Setyawan E, Laglbauer BJ, Lewis S, Dharmadi D, Sianipar A, Ender I (2020) Predicting mobulid ray distribution in coastal areas of Lesser Sunda Seascape: implication for spatial and fisheries management. Ocean Coast Manag 198:105328. https://doi.org/10.1016/j.ocecoaman.2020.105328
Raymond M, Rousset F (1995) GENEPOP (version 1.2): population genetics software for exact tests and ecumenicism. J Hered 86:248–249. https://doi.org/10.1093/oxfordjournals.jhered.a111573
Rice WR (1989) Analyzing tables of statistical tests. Evolution 43:223–225. https://doi.org/10.2307/2409177
Rousset F (2008) Genepop’007: a complete reimplementation of the Genepop software for Windows and Linux. Mol Ecol Resour 8:103–106. https://doi.org/10.1111/j.1471-8286.2007.01931.x
Rozas J, Ferrer-Mata A, Sánchez-DelBarrio JC, Guirao-Rico S, Librado P, Ramos-Onsins SE, Sánchez-Gracia A (2017) DnaSP 6: DNA sequence polymorphism analysis of large data sets. Mol Biol Evol 34(12):3299–3302. https://doi.org/10.1093/molbev/msx248
Sahri A, Putra MIH, Mustika PLK, Kreb D, Murk AJ (2021) Cetacean habitat modelling to inform conservation management, marine spatial planning, and as a basis for anthropogenic threat mitigation in Indonesia. Ocean Coast Manag 205:105555. https://doi.org/10.1016/j.ocecoaman.2021.105555
Sahri A, Jak C, Putra MIH, Murk AJ, Andrews-Goff V, Double MC, Van Lammeren RJ (2022) Telemetry-based home range and habitat modelling reveals that the majority of areas important for pygmy blue whales are currently unprotected. Biol Conserv 272:109594. https://doi.org/10.1016/j.biocon.2022.109594
Saitou N, Nei M (1987) The neighbor-joining method: a new method for reconstructing phylogenetic trees. Mol Biol Evol 4(4):406–425. https://doi.org/10.1093/oxfordjournals.molbev.a040454
Smith PJ, Francis RICC, McVeagh M (1991) Loss of genetic diversity due to fishing pressure. Fish Res 10(3–4):309–316. https://doi.org/10.1016/0165-7836(91)90082-Q
Spaet JL, Jabado RW, Henderson AC, Moore AB, Berumen ML (2015) Population genetics of four heavily exploited shark species around the Arabian Peninsula. Ecol Evol 5(12):2317–2332. https://doi.org/10.1002/ece3.1515
Underwood JN, Smith LD, van Oppen MJH, Gilmour JP (2007) Multiple scales of genetic connectivity in a brooding coral on isolated reefs following catastrophic bleaching. Mol Ecol 16:771–784. https://doi.org/10.1111/j.1365-294X.2006.03187.x
Van Oosterhout C, Hutchinson WF, Wills DP, Shipley P (2004) MICRO-CHECKER: software for identifying and correcting genotyping errors in microsatellite data. Mol Ecol Notes 4(3):535–538. https://doi.org/10.1111/j.1471-8286.2004.00684.x
Venables SK, Marshall AD, Armstrong AJ, Tomkins JL, Kennington WJ (2021) Genome-wide SNPs detect no evidence of genetic population structure for reef manta rays (Mobula alfredi) in southern Mozambique. Heredity 126(2):308–319. https://doi.org/10.1038/s41437-020-00373-x
Wardana ED, Putra ING, Malik MDA, Yusmalinda NLA, Niningsih EY, Pertiwi NPD, Salim MG, Kanedi MH, Putra MIH, Sembiring A (2023) Population genetic structure of the bentfin devil ray (Mobula thurstoni) in the South Indonesia Sea with limited sample based on ND5 gene. Biodiversitas 24(7):1. https://doi.org/10.13057/biodiv/d240711
Ward-Paige CA, David B, Worm B (2013) Global population trends and human use patterns of Manta and Mobula rays. PLoS ONE 8(9):e74835. https://doi.org/10.1371/journal.pone.0074835
Whitaker K (2004) Non-random mating and population genetic subdivision of two broadcast spawning corals at Ningaloo reef, Western Australia. Mar Biol 144:593–603. https://doi.org/10.1007/s00227-003-1220-7
Wright S (1965) The interpretation of population structure by F-statistics with special regard to systems of mating. Evolution 19(3):395–420. https://doi.org/10.2307/2406450
Yeh FC, Boyle T (1997) POPGENE version 1.2, Microsoft Windows-based software for population genetics analysis. University of Alberta, Alberta, pp 180–189
You EM, Chiu TS, Liu KF, Tassanakajon A, Klinbunga S, Triwitayakorn K, de la Pena YL, Yu HT (2008) Microsatellite and mitochondrial haplotype diversity reveals population differentiation in the tiger shrimp (Penaeus monodon) in the Indo-Pacific region. Anim Genet 39(3):267–277. https://doi.org/10.1111/j.1365-2052.2008.01724.x
Acknowledgements
We thank our funders Vulcan Inc., MacArthur Foundation, Shark Conservation Fund, and WildAid Canada Society, for supporting our research program in the Savu Sea as a part of Lamakera Manta Conservation Project. This field activity and genetic analysis also has been funded and supported by the Rufford Foundation (small grant refs: 30112-1 and 34956-2), the IDEA WILD (MALIINDO0122), National Academy of Sciences (NAS), and USAID under the USAID Prime Award Number AID-OAA-A-11-00012, and the Pew Fellows Program in Marine Conservation at The Pew Charitable Trusts. In addition, we extend our thanks to the Secretary of DJPRL KKP for their support of this research through MOU 02/PRL/KKP/PKS/II/2021. Any opinions, findings, conclusions, or recommendations expressed in this material are those of the authors and do not necessarily reflect the views of the funders.
Author information
Authors and Affiliations
Contributions
MDAM: conceptualization, data curation, methodology, formal analysis, resources, visualization, investigation, writing—original draft, funding acquisition. MIHP: conceptualization, data curation, investigation, resources, visualization, writing—review and editing. ET: data curation, investigation. NLAY: project administration, validation, formal analysis. NPDP: writing—review and editing. YFS: validation, writing—review and editing. NKDC: writing—review and editing, supervision. EYA: investigation. SL: funding acquisition. LNLT: writing—review and editing. MGS: data curation, writing—review and editing. FT: data curation. FAA: data curation, writing—review and editing. DP: funding acquisition, writing—review and editing. AS: writing—review and editing, funding acquisition, supervision.
Corresponding author
Ethics declarations
Conflict of interest
The authors declare no conflict of interest.
Additional information
Publisher’s Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary Information
Below is the link to the electronic supplementary material.
Rights and permissions
Springer Nature or its licensor (e.g. a society or other partner) holds exclusive rights to this article under a publishing agreement with the author(s) or other rightsholder(s); author self-archiving of the accepted manuscript version of this article is solely governed by the terms of such publishing agreement and applicable law.
About this article
Cite this article
Al Malik, M.D., Putra, M.I.H., Topan, E. et al. Population structure of endangered spinetail devil ray (Mobula mobular) in the Lesser Sunda Seascape, Indonesia, revealed using microsatellite and mitochondrial DNA. Aquat Sci 86, 6 (2024). https://doi.org/10.1007/s00027-023-01020-3
Received:
Accepted:
Published:
DOI: https://doi.org/10.1007/s00027-023-01020-3