60urn:lsid:arphahub.com:pub:056820A7-C438-5162-B00B-FC18BD6C2AA0urn:lsid:zoobank.org:pub:4B0FB9C5-4BE9-4A41-8BA6-2C2FD3522FC1ZoologiaZoologia1984-4689Pensoft Publishers10.3897/zoologia.38.e5933259332Short CommunicationActinopterygiiAnimaliaChordataOsteichthyesPiscesVertebrataBiodiversity & ConservationBiogeographyEcology & Environmental sciencesHabitats Ecosystems & Natural SpacesPopulation Management Harvesting and MonitoringPopulations & CommunitiesDNA barcoding of a tropical anguillid eel, Anguillabicolor (Actinopterygii: Anguilliformes), in Indo-Pacific region and notes on its population structureNorarfanAhmad Fathi1Azreena MoktiSiti Shazwani1TahaHusseinhttps://orcid.org/0000-0003-3272-43131AminMuhamadhttps://orcid.org/0000-0002-7633-83882AliMuhamad3AraiTakaomitakaomi.arai@ubd.edu.bnhttps://orcid.org/0000-0001-9440-79331Environmental and Life Sciences Programme, Faculty of Science, Universiti Brunei Darussalam. Jalan Tungku Link, Gadong, BE 1410, Brunei Darussalam.Universiti Brunei DarussalamBandar Seri BegawanBruneiDepartment of Aquaculture, Faculty of Fisheries and Marine, Universitas Airlangga. Surabaya, Indonesia.Universitas AirlanggaSurabayaIndonesiaLaboratory of Microbiology and Biotechnology, Faculty of Animal Science. Universitas Mataram. West-Nusa, Tenggara, Indonesia.Universitas MataramWest-Nusa TenggaraIndonesia
202124032021381780C949A9-333D-5574-88FD-173BC1F9EF0C9F3E23D8-C236-4629-A270-03551CD55E1246552590510202009122020Ahmad Fathi Norarfan, Siti Shazwani Azreena Mokti, Hussein Taha, Muhamad Amin, Muhamad Ali, Takaomi AraiThis 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/9F3E23D8-C236-4629-A270-03551CD55E12
The tropical anguillid eel, Anguillabicolor McCelland, 1844, includes two subspecies, Anguillabicolorbicolor McCelland, 1844 and Anguillabicolorpacifica Schmidt, 1928, and is distributed across the Indo-Pacific region. Although A.bicolor is widely distributed and recognized as an important fish resource in the Indo-Pacific region, few studies have been conducted on its genetic variation and population structure. DNA barcoding of A.bicolor specimens collected in the Indo-Pacific region was carried out in this study using mitochondrial cytochrome c oxidase subunit I. Anguillabicolor was found to diverge genetically, which supported its classification into two different subspecies. In addition, our study showed that A.bicolorbicolor had two genetically distinct populations/groups, and these different populations co-occur geographically in Indonesia and Malaysia in the eastern Indian Ocean. Our findings suggest that the eel larvae might be transported from at least two geographically different spawning grounds in the Indian Ocean, and then recruited to and settled in the same habitats in Indonesian and Malaysian waters. The molecular evidence calls for further research on the life history, stock assessment and protection of the populations of A.bicolorbicolor in Indonesia and Malaysia.
Anguillid eelhaplotype networkmolecular phylogenytropical species
Sixteen species of the genus Anguilla Schrank, 1798 are widely distributed across the world. Eleven species are found in tropical Indo-Pacific regions (Ege 1939, Arai 2016, 2020). Of the eleven species, six species, i.e., Anguillamarmorata Quoy & Gaimard, 1824, Anguillabengalensis Gray, 1831, Anguillacelebesensis Kaup, 1856, Anguillaborneensis Popta, 1924, Anguillainterioris Whitley, 1938, and Anguillabicolor McCelland, 1844, are distributed in the tropical waters of Indonesia and Malaysia. Anguillid eels have a catadromous life history, migrating between offshore spawning grounds and freshwater growth habitats (Arai 2020). Many tropical anguillid species share similar geographic distributions and spawning grounds (Tesch 2003), which raises the question of whether the different species hybridize naturally.
The short-finned tropical eel, A.bicolor, is distributed globally across the Indo-Pacific region, from East Africa to Papua New Guinea, and the Philippines through Indonesia (Ege 1939). Due to the wide geographical distribution, offshore spawning migration, passive larval transportation and the absence of significant physical barriers in the Indo-Pacific Ocean, a weak or an absence of population structure was expected for A.bicolor (Minegishi et al. 2012, Fahmi et al. 2015), as observed in temperate species such as Anguillaanguilla Linnaeus, 1758 (Palm et al. 2009) and Anguillajaponica Temminck & Schlegel, 1846 (Han et al. 2008). However, based on its allopatric distribution and slight differences in morphological characteristics, A.bicolor has been subdivided into A.bicolorbicolor (Indian Ocean) and A.bicolorpacifica (western Pacific Ocean) (Ege 1939, Watanabe et al. 2014).
Anguillabicolor was previously sampled across its entire geographic range and analyzed using microsatellites and mitochondrial control region sequences (Minegishi et al. 2012). They found that A.bicolor was divided into (i) a homogeneous eastern and western side of the Indian Ocean (A.bicolorbicolor), and (ii) a western side of the Pacific Ocean (A.bicolorpacifica). Analysis of mitochondrial cytochrome b sequences and microsatellite markers further confirmed the differentiation into two subspecies with the microsatellite loci, suggesting a moderate differentiation between the subspecies (Fahmi et al. 2015). Although A.bicolorbicolor in the Indian Ocean was found to have two mitochondrial sublineages with random occurrence among and within sampling locations (Minegishi et al. 2012, Fahmi et al. 2015), little work has been done on the population structure of A.bicolorbicolor in the Indian Ocean.
Knowledge on the life history of A.bicolor has been gradually accumulating (e.g., Arai et al. 1999a, 1999b, 2020, Robinet et al. 2003, Arai and Chino 2013, 2019, Arai and Abdul Kadir 2017a, 2017b, Chai and Arai 2018, Arai et al. 2020). Hatching and recruitment are estimated to be almost year-round and at a constant age, respectively (Arai et al. 2001). Tropical glass eels are recruited to a river mouth throughout the year (Arai et al. 1999c, Sugeha et al. 2001), which corresponds to a year-round spawning (Arai et al. 2016, Arai and Abdul Kadir 2017a). In contrast, temperate eels do not show year-round spawning and recruitment. The broad ranges of otolith strontium/calcium (Sr:Ca) ratios found in A.bicolorbicolor and A.bicolorpacifica suggest that these subspecies have flexible migratory behaviors in ambient waters (Chino and Arai 2010a, 2010b, Arai et al. 2013, 2020, Arai and Chino 2019).
The aim of the present study was to investigate the genetic variation through DNA barcoding of A.bicolor in Indonesian, Malaysian and Philippine waters. All eel specimens were examined for their morphological characteristics and cytochrome c oxidase subunit I (COI) gene sequence. We also discuss the potential dispersion mechanism and importance of the population structure of A.bicolor in these regions.
A total of 18 anguillid eels were collected from Southeast Asia: nine specimens in Indonesia (Cilacap and Yogyakarta, central Java), four specimens in Malaysia (Penang, northwestern Peninsular Malaysia) and five specimens in the Philippines (Mindanao) through traps, hooks and lines from 2009 to 2019 (Fig. 1, Table 1). The external morphometric characteristics were determined (Table 1). The fin difference index (FDI), which is the relative ratio of the ano-dorsal (the origin of the dorsal fin to the anus) length (Z) and the total length (LT) (Ege 1939), was calculated as follows: FDI = 100 Z LT-1.
Morphometric characteristics of tropical anguillid eels collected in Southeast Asia.
Map of sampling locations of Anguillabicolorbicolor from Indonesia and Malaysia and A.bicolorpacifica from the Philippines. Sampling locations are shown in red. Locations of Banda Aceh in northern Sumatra of Indonesia, Sukabumi in western Java of Indonesia, Andhra Pradesh in southeast India, Langkawi and Penang islands of northwest Peninsular Malaysia, Cebu in the Philippines, and Yangon of Myanmar, which were used in molecular phylogenetic and haplotype network analyses, are shown in blue. Base map is downloaded from http://viewer.nationalmap.gov/viewer (USGS 2020).
https://binary.pensoft.net/fig/523889
DNeasy Blood & Tissue Kit (QIAGEN, Germany) was used to extract genomic DNA from a dorsal fin clip. Approximately 600 bp region of the mitochondrial COI gene was amplified using different combinations of universal primers: FishF1 (5’TCA ACC AAC CAC AAA GAC ATT GGC AC3’), FishF2 (5’TCG ACT AAT CAT AAA GAT ATC GGC AC3’), FishR1 (5’TAG ACT TCT GGG TGG CCA AAG AAT CA3’), and FishR2 (5’ACT TCA GGG TGA CCG AAG AAT CAG AA3’) (Ward et al. 2005). The samples that could not be amplified by these primers were subsequently amplified with the primers, AngF (5’TCA CCC GTT GAT TCT TTT CT3’) and AngR (5’CCG ATA GCC ATT ATT GCT CAG3’), which were designed in this study to amplify 835 bp of COI (starting from position 8 until 842) from Anguilla eels. The primers were designed by targeting the conserved sites of the COI genes from different species of sequenced Anguilla eels (Minegishi et al. 2005). Each PCR reaction contained 2 µl of DNA sample, 2.5 µl of each 10 µM primer, 25 µl of 2x Taq PCR Master Mix (QIAGEN, Germany) and 18 µl of distilled water. The PCR conditions were initially 95°C for 2 mins, then 35 cycles of 94°C for 30 s, 50°C for 30 s and 72°C for 60 s, and finally 72°C for 10 minutes. PCR amplicons were purified using QIAquick Gel Extraction Kit (QIAGEN, Germany), and sequenced bi-directionally with the same primers.
MEGA version X (Kumar et al. 2018) was used to edit and assemble DNA sequence trace files. To validate the species identity, the contig sequences were compared for percent similarity with the reference sequences in the GenBank database (Benson et al. 2013) by using BLAST search (Altschul et al. 1990). The contig sequences of eels from Indonesia, Malaysia and Philippines were submitted to the GenBank database with accession numbers MT107677 – MT107685, MT107673 – MT107676 and MT107686 – MT107690, respectively. In addition to these eighteen sequences, the COI sequences of thirty A.bicolorbicolor specimens from other localities that were found in the GenBank were also included in the analysis: Banda Aceh in northern Sumatra of Indonesia (22 specimens; KY618768 – KY618795), Sukabumi in western Java of Indonesia (one specimen; KU692247), Andhra Pradesh in southeast India (two specimens; MG675613, KP979655), Langkawi (one specimen; KF182304) and Penang islands of northwest Peninsular Malaysia (3 specimens; KM875503 – KM875505), and Yangon of Myanmar (one specimen; AP007236). The COI sequences of A.bicolorpacifica from Cebu in the Philippines (one specimen; AP007237) and Anguillamegastoma Kaup, 1856 (as outgroup; AP007243) were also included (Fig. 1).
GUIDANCE version 2 (Sela et al. 2015) was used to create multiple sequence alignment via its inbuilt Clustal W version 2 (Larkin et al. 2007), which also assigned the alignment with a robust score of 1. The multiple sequence alignment was trimmed at both ends, resulting in all sequences having similar length of 445 bp. No internal gap was observed in the alignment. MEGA was used to construct three phylogenetic trees via Neighbour Joining (NJ), Maximum Parsimony (MP) and Maximum Likelihood (ML) algorithms. The NJ tree was constructed using the Kimura 2-parameter (K2P) model (incorporating transitions and transversions) with uniform rates. The MP tree was searched using Tree-Bisection-Reconnection with search level 1. The initial tree was obtained by 1000 random addition replicates. A total of 100 most parsimonious trees were retained. The ML tree was constructed with the best-fit nucleotide substitution model, K2P with discrete Gamma distribution (+G; 5 categories), which was determined using MEGA. A heuristic search starting with the initial tree (based on NJ and BioNJ algorithms) was conducted using the Nearest-Neighbour-Interchange method with no branch swap filter. Bootstrap test with 1000 replicates was also carried out for the NJ, MP and ML analyses. Haplotype analysis was carried out using DnaSP version 6 (Rozas et al 2017), and a median joining haplotype network was constructed using Network version 10 (www.fluxus-engineering.com). Statistical analysis was not carried out because the number of samples in some localities was too low to produce meaningful results.
Eighteen eel samples had short dorsal fin and skin without variegated markings. All the eel samples had FDI values of -0.3 to 4.7 % (Table 1). Based on the non-variegated skin and short dorsal fin, the eel specimens could be any species in the following list: A.bicolorbicolor, A.bicolorpacifica, Anguillaobscura Günther, 1872, Anguillaaustralisaustralis Richardson, 1841, Anguillaaustralisschmidtii Phillipps, 1925 (Ege 1939, Watanabe et al. 2004, Arai 2016). To further determine which species from the list, diagnostic geographical distribution (Ege 1939) was used. Thirteen specimens from Indonesia and Malaysia and five specimens from the Philippines were identified as A.bicolorbicolor and A.bicolorpacifica, respectively (Table 1).
DNA barcoding using the COI gene confirmed that the 13 specimens from Indonesia and Malaysia and five specimens from the Philippines were A.bicolorbicolor and A.bicolorpacifica, respectively (Table 1), with 99–100 % maximum identity matches with the reference sequences in the GenBank database. Apart from the universal primers used, this study showed that DNA barcoding could also be carried out using novel primers that were designed for Anguilla eels.
From our phylogenetic analyses (Fig. 2), A.bicolorbicolor samples were observed to form a clade (bootstrap proportion = >60 %) that was distinct from the clade (bootstrap proportion ≥ 90 %) formed by A.bicolorpacifica samples. This confirms the differentiation of A.bicolor into two subspecies. The MP tree strongly supports the presence of two subgroups within the A.bicolorbicolor clade (bootstrap proportion: subgroup 1 = 80%, subgroup 2 = 66%). These two subgroups are also found in the NJ tree, though weakly supported by bootstrap proportion (subgroup 1 = 50%, subgroup 2 = 58%). However, these two subgroups were not clearly observed in the ML tree.
Haplotype analysis revealed a total of 18 different haplotypes in A.bicolorbicolor (H1 – H18 in Fig. 3) and 2 different haplotypes in A.bicolorpacifica (H19–H20 in Fig. 3). The constructed haplotype network (Fig. 3) further confirmed the presence of two mitochondrial sublineages in A.bicolorbicolor, as the haloptypes could be grouped into two haplogroups or populations.
Maximum Parsimony tree using the COI gene sequences of Anguillabicolorbicolor and Anguillabicolorpacifica and additional sequences obtained from the GenBank database with indicated accession numbers. Anguillamegastoma was used as outgroup. The consistency index is 0.73, the retention index is 0.93, and the composite index is 0.82 for all sites and 0.68 for parsimony-informative sites. The percentages of the bootstrap test for Maximum Parsimony/Neighbour Joining/Maximum Likelihood are shown next to the branches.
Haplotype network constructed with AnguillabicolorCOI gene sequences. Each colour represents a sample site. The size of the circle is proportional to the number of samples that belong to each haplotype. White circle refers to median vector which is a hypothesized or missing haplotype. Each dash, which appears on the line that connects two haplotypes together, symbolizes one mutational step.
https://binary.pensoft.net/fig/523891
Anguillabicolor from the Indian Ocean (A.bicolorbicolor) and Pacific Ocean (A.bicolorpacifica) were found to diverge from each other, and the outcomes of our COI gene sequence were similar to those previous studies using microsatellites, mitochondrial control region and cytochrome b sequences (Minegishi et al. 2012, Fahmi et al. 2015). The previous studies suggest that A.bicolor has diverged into two subspecies in the two oceans. In the eastern and western sides of the Indian Ocean, the genetic profiles of A.bicolorbicolor were found to be similar. However, the genetic profile of either the eastern or western side of the Indian Ocean was observed to have two divergent mitochondrial sublineages (Minegishi et al. 2012, Fahmi et al. 2015). These two sublineages coexist in the Indian Ocean, occurring randomly among and within sampling locations. The present results also supported the occurrence of two mitochondrial sublineages, as there were two distinctive groups in the eastern Indian ocean (Figs 2, 3).
Based on the COI sequence analysis, the results suggest that different A.bicolorbicolor populations could co-occur geographically in Indonesian and Malaysian waters. Spawning areas of A.bicolorbicolor are thought to be located in the waters off west Sumatra (Jespersen 1942), and off east Madagascar (Jubb 1961, Robinet et al. 2003). These two geographically different spawning areas might explain the occurrence of two divergent populations of A.bicolorbicolor in the Indian Ocean. Thus, the present study supports the hypothetical occurrence of the two spawning grounds of A.bicolorbicolor in the Indian Ocean. Furthermore, the ages of recruitment of this species between the east (148–202 days in Java of Indonesia) and west coasts (68–96 days in Mascarene Islands) of the Indian Ocean definitely differed (Arai et al. 1999a, Robinet et al. 2003). This suggests that A.bicolorbicolor could be reasonably assumed to have different populations in the east and west sides of the Indian Ocean.
Anguillabicolorbicolor in the eastern Indian Ocean has a longer leptocephalus period (119–171 days) (Arai et al. 1999a) and is spawning throughout the year (Arai et al. 2016, Arai and Abdul Kadir 2017a). These life history strategies might enable A.bicolorbicolor to disperse across the Indian Ocean. Seasonal changes of the South Equatorial Current and the long larval migration period in the Indian Ocean, which was estimated to be more than one year in anguillid eels (Pous et al. 2010), might be capable of mixing A.bicolorbicolor larvae from different spawning areas. Further studies are needed to determine the exact spawning grounds and larval transportation mechanism and routes in A.bicolorbicolor. Nevertheless, these particular life history characteristics of A.bicolorbicolor and the oceanographic features of the Indian Ocean could explain the sympatric occurrence of different populations.
The specific choice of spawning area might be imprinted in anguillid eels including A.bicolorbicolor, although there is no physical barrier to migration to the different spawning areas. The occurrence of two populations in the eastern Indian Ocean suggests that population divergence in anguillid eels could still occur even in the absence of physical barriers during the spawning migration. Minegishi et al. (2012) presented a different hypothesis, as this previous study found that A.bicolorbicolor diverged into two sublineages based on the mitochondrial control region but not when using the microsatellite markers. The study hypothesizes that there could be an incomplete lineage sorting or a possible secondary contact between the two groups after an initial population splitting.
This study focused only in the eastern side but did not include the western side of Indian Ocean. The sample sizes were also small in all localities. Therefore, more comprehensive research is needed to understand the population structure of A.bicolorbicolor across the Indian Ocean. This study has revealed a unique population structure in A.bicolorbicolor, revealing a clear genetic divergence in the eastern Indian Ocean. The molecular evidence calls for further research on the life history, stock assessment and protection of the populations of A.bicolorbicolor in Indonesia and Malaysia.
Acknowledgements
This study was financially supported by University Brunei Darussalam under the Competitive Research Grant Scheme (UBD/OVACRI/CRGWG(003)) and under the Faculty/Institute/Center Research Grant (UBD/RSCH/1.4/FICBF(b)/2019/021 and UBD/RSCH/1.4/FICBF(b)/2020/029).
Literature citedAltschulSFGishWMillerWMyersEWLipmanDJ (1990) Basic local alignment search tool.Journal of Molecular Biology215: 403–410. https://doi.org/10.1016/S0022-2836(05)80360-2AraiT (2016) Taxonomy and distribution. In: AraiT (Ed.) Biology and ecology of anguillid eels.Boca Raton, CRC Press, 1–20.AraiT (2020) Ecology and evolution of migration in the freshwater eels of the genus Anguilla Schrank, 1798. Heliyon 6: e05176. https://doi.org/10.1016/j.heliyon.2020.e05176AraiTAbdulKadir SR (2017a) Opportunistic spawning of tropical anguillid eels Anguillabicolorbicolor and A.bengalensisbengalensis. Scientific Reports 7: 41649. https://doi.org/10.1038/srep41649AraiTAbdulKadir SR (2017b) Diversity, distribution and different habitat use among the tropical freshwater eels of genus Anguilla. Scientific Reports 7: 7593. https://doi.org/10.1038/s41598-017-07837-xAraiTAbdulKadir SRChinoN (2016) Year-round spawning by a tropical catadromous eel Anguillabicolorbicolor. Marine Biology 163: 37. https://doi.org/10.1007/s00227-015-2792-8AraiTAoyamaJLimbongDTsukamotoK (1999c) Species composition and inshore migration of the tropical eels Anguilla spp. recruiting to the estuary of the Poigar River, Sulawesi Island.Marine Ecology Progress Series188: 299–303. https://doi.org/10.3354/meps188299AraiTChaiIJIizukaYChangCW (2020) Habitat segregation and migration in tropical anguillid eels, Anguillabengalensisbengalensis and A.bicolorbicolor. Scientific Reports 10: 16890. https://doi.org/10.1038/s41598-020-72788-9AraiTChinoN (2013) Timing of maturation of a tropical eel, Anguillabicolorbicolor in Malaysia.Journal of Applied Ichthyology29: 271–273. https://doi.org/10.1111/jai.12040AraiTChinoN (2019) Variations in the migratory history of the tropical catadromous eels Anguillabicolorbicolor and A.bicolorpacifica in south-east Asian waters.Journal of Fish Biology94: 752–758. https://doi.org/10.1111/jfb.13952AraiTChinoNLeDQ (2013) Migration and habitat use of the tropical eels Anguillamarmorata and A.bicolorpacifica in Vietnam.Aquatic Ecology47: 57–65. https://doi.org/10.1007/s10452-012-9424-xAraiTLimbongDOtakeTTsukamotoK (1999a) Metamorphosis and inshore migration of tropical eels Anguilla spp. in the Indo-Pacific.Marine Ecology Progress Series182: 283–293. https://doi.org/10.3354/meps182283AraiTLimbongDOtakeTTsukamotoK (2001) Recruitment mechanisms of tropical eels, Anguilla spp., and implications for the evolution of oceanic migration in the genus Anguilla.Marine Ecology Progress Series216: 253–264. https://doi.org/10.3354/meps216253AraiTOtakeTLimbongDTsukamotoK (1999b) Early life history and recruitment of the tropical eel Anguillabicolorpacifica, as revealed by otolith microstructure and microchemistry.Marine Biology133: 319–326. https://doi.org/10.1007/s002270050470BensonDACavanaughMClarkKKarsch-MizrachiILipmanDJOstellJSayersEW (2013) GenBank. Nucleic Acids Research 41(D1): D36–D42. https://doi.org/10.1093/nar/gks1195ChaiIJAraiT (2018) Ages at maturation of tropical freshwater eels, Anguillabicolorbicolor and A.bengalensisbengalensis.Journal of Applied Animal Research46: 1108–1113. https://doi.org/10.1080/09712119.2018.1470090ChinoNAraiT (2010a) Occurrence of marine resident tropical eel Anguillabicolorbicolor in Indonesia.Marine Biology157: 1075–1081. https://doi.org/10.1007/s00227-009-1388-6ChinoNAraiT (2010b) Habitat use and habitat transitions in the tropical eel, Anguillabicolorbicolor.Environmental Biology of Fishes89: 571–578. https://doi.org/10.1007/s10641-010-9677-yEgeV (1939) A revision of the genus Anguilla Shaw.Dana Report16: 8–256.FahmiMRSolihinDDShaoZPouyaudLBerrebiP (2015) Population genetic structure of the tropical eel Anguillabicolor in Indonesian waters based on microsatellite markers.Journal of Vertebrate Biology64: 87–96. https://doi.org/10.25225/fozo.v64.i2.a1.2015HanYSSunYLLiaoYFShenKNLiaoICTzengWN (2008) Temporal analysis of population genetic composition in the overexploited Japanese eel Anguillajaponica.Marine Biology155: 613–621. https://doi.org/10.1007/s00227-008-1057-1JespersenP (1942) Indo-Pacific leptochephalids of the genus Anguilla.Dana Report22: 1–128.JubbRA (1961) The freshwater eels (Anguilla spp.) of southern Africa: an introduction to their identification and biology.Annuals of the Cape Provincial Museum1: 15–48.KumarSStecherGLiMKnyazCTamuraK (2018) MEGA X: Molecular Evolutionary Genetics Analysis across computing platforms.Molecular Biology and Evolution35: 1547–1549. https://doi.org/10.1093/molbev/msy096LarkinMABlackshieldsGBrownNPChennaRMcGettiganPAMcWilliamHValentinFWallaceIMWilmALopezRThompsonJDGibsonTJHigginsDG (2007) Clustal W and Clustal X version 2.0.Bioinformatics23(21): 2947–2948. https://doi.org/10.1093/bioinformatics/btm404MinegishiYAoyamaJInoueJGMiyaMNishidaMTsukamotoK (2005) Molecular phylogeny and evolution of the freshwater eels genus Anguilla based on the whole mitochondrial genome sequences.Molecular Phylogenetics and Evolution34: 134–146. https://doi.org/10.1016/j.ympev.2004.09.003MinegishiYGagnairePAAoyamaJBoscPFeunteunETsukamotoKBerrebiP (2012) Present and past genetic connectivity of the Indo-Pacific tropical eel Anguillabicolor.Journal of Biogeography39: 408–420. https://doi.org/10.1111/j.1365-2699.2011.02603.xPalmSDannewitzJPrestegaardTWickströmH (2009) Panmixia in European eel revisited: no genetic difference between maturing adults from southern and northern Europe.Heredity103: 82–89. https://doi.org/10.1038/hdy.2009.51PousSFeunteunEEllienC (2010) Investigation of tropical eel spawning area in the South-Western Indian Ocean: influence of the oceanic circulation.Progress in Oceanography86: 396–413. https://doi.org/10.1016/j.pocean.2010.06.002RobinetTLecomte-FinigerREscoubeyrouKFeunteunE (2003) Tropical eels Anguilla spp. recruiting to Reunion Island in the Indian Oceans: taxonomy, patterns of recruitment and early life histories.Marine Ecology Progress Series259: 263–272. https://doi.org/10.3354/meps259263RozasJFerrer-MataASánchez-DelBarrioJCGuirao-RicoSLibradoPRamos-OnsinsSESánchez-GraciaA (2017) DnaSP 6: DNA sequence polymorphism analysis of large datasets.Molecular Biology and Evolution34: 3299–3302. https://doi.org/10.1093/molbev/msx248SelaIAshkenazyHKatohKPupkoT (2015) GUIDANCE2: accurate detection of unreliable alignment regions accounting for the uncertainty of multiple parameters. Nucleic Acids Research 43: W7–W14. https://doi.org/10.1093/nar/gkv318SugehaHYAraiTMillerMJLimbongDTsukamotoK (2001) Inshore migration of the tropical eels Anguilla spp. Recruiting to the Poigar River estuary on North Sulawesi Island.Marine Ecology Progress Series221: 233–243. https://doi.org/10.3354/meps2212TeschFW (2003) .The eel.Blackwell Publishing, London, 3rd ed., 408 pp.USGS (2020) National map viewer. USA Department of the Interior, United States Geological Survey. http://viewer.nationalmap.gov/viewer [Accessed: 05/10/2020]WardRDZemlakTSInnesBHLastPRHerbertPDN (2005) DNA barcoding Australia’s fish species.Philosophical Transactions of the Royal Society B: Biological Sciences360: 1847–1857. https://doi.org/10.1098/rstb.2005.1716WatanabeSAoyamaJTsukamotoK (2004) Reexaminationof Ege’s (1939) use of taxonomic characters of the genus Anguilla.Bulletin of Marine Science74: 337–351.WatanabeSMillerMJAoyamaJTsukamotoK (2014) Evaluation of the population structure of Anguillabicolor and A.bengalensis using total number of vertebrae and consideration of the subspecies concept for the genus Anguilla.Ecology of Freshwater Fish23: 77–85. https://doi.org/10.1111/eff.12076