Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

DNA barcoding of coastal ray-finned fishes in Vietnam

  • Pham The Thu ,

    Contributed equally to this work with: Pham The Thu, Wen-Chien Huang

    Roles Conceptualization, Data curation, Formal analysis, Funding acquisition, Investigation, Methodology, Project administration, Resources, Supervision, Validation, Writing – original draft, Writing – review & editing

    Affiliation Institute of Marine Environment and Resources, Vietnam Academy of Science and Technology, Hai Phong, Vietnam

  • Wen-Chien Huang ,

    Contributed equally to this work with: Pham The Thu, Wen-Chien Huang

    Roles Conceptualization, Data curation, Formal analysis, Investigation, Methodology, Project administration, Software, Supervision, Validation, Writing – original draft, Writing – review & editing

    Affiliations Doctoral Degree Program in Marine Biotechnology, National Sun Yat-sen University, Kaohsiung, Taiwan, Doctoral Degree Program in Marine Biotechnology, Academia Sinica, Taipei, Taiwan, Department of Oceanography, National Sun Yat-sen University, Kaohsiung, Taiwan

  • Tak-Kei Chou,

    Roles Data curation, Formal analysis, Investigation, Methodology, Software, Validation

    Affiliation Department of Oceanography, National Sun Yat-sen University, Kaohsiung, Taiwan

  • Nguyen Van Quan,

    Roles Conceptualization, Funding acquisition, Project administration, Resources

    Affiliation Institute of Marine Environment and Resources, Vietnam Academy of Science and Technology, Hai Phong, Vietnam

  • Pham Van Chien,

    Roles Investigation, Resources

    Affiliation Institute of Marine Environment and Resources, Vietnam Academy of Science and Technology, Hai Phong, Vietnam

  • Fan Li,

    Roles Data curation, Investigation

    Affiliation Department of Oceanography, National Sun Yat-sen University, Kaohsiung, Taiwan

  • Kwang-Tsao Shao,

    Roles Conceptualization, Funding acquisition

    Affiliations Biodiversity Research Center, Academia Sinica, Taipei, Taiwan, Institute of Marine Biology, National Taiwan Ocean University, Keelung, Taiwan

  • Te-Yu Liao

    Roles Conceptualization, Data curation, Funding acquisition, Investigation, Methodology, Project administration, Resources, Supervision, Validation, Writing – original draft, Writing – review & editing

    swp0117@gmail.com

    Affiliation Department of Oceanography, National Sun Yat-sen University, Kaohsiung, Taiwan

Abstract

DNA barcoding based on a fragment of the cytochrome c oxidase subunit I (COI) gene is widely applied in species identification and biodiversity studies. The aim of this study was to establish a comprehensive barcoding database of coastal ray-finned fishes in Vietnam. A total of 3,638 specimens were collected from fish landing sites in northern, central and southern Vietnam. Seven hundred and sixty-five COI sequences of ray-finned fishes were generated, belonging to 458 species, 273 genera, 113 families and 43 orders. A total of 59 species were newly recorded in Vietnam and sequences of six species were new to the Genbank and BOLD online databases. Only 32 species cannot be annotated to species level because difficulty in morphological identifications and their Kimura-2-Parameter (K2P) genetic distances to most similar sequences were more than 2%. Moreover, intra-specific genetic distances in some species are also higher than 2%, implying the existence of putative cryptic species. The mean K2P genetic distances within species, genera, families, orders and classes were 0.34%, 12.14%, 17.39%, 21.42%, and 24.80, respectively. Species compositions are quite different with only 16 common species among northern, central and southern Vietnam. This may attribute to multiple habitats and environmental factors across the 3,260 km Vietnamese coastline. Our results confirmed that DNA barcoding is an efficient and reliable tool for coastal fish identification in Vietnam, and also established a reliable DNA barcode reference library for these fishes. DNA barcodes will contribute to future efforts to achieve better monitoring, conservation, and management of fisheries in Vietnam.

Introduction

The DNA-based barcoding method has been proven to be a valuable molecular tool for species identification and it is accessible to non-specialists [13]. Currently, DNA barcoding has been employed to a large variety of organisms ranging from yeasts to humans with different standard region of genes for various organisms, including the gene of COI (cytochrome c oxidase 1) for animals [4] and it has been shown effective in discriminating fish species [57]. FISH-BOL [8] is also one of fish barcode database that is now well established and aims at DNA barcoding all the fishes of the world [8], with a total of approximately 7,800 fish species barcoded, among which 2019 species are distributed in Southeast Asia [9, 10]. When the reference sequence library is in place, new specimens and products can be identified by comparing their DNA barcode sequences against this barcode reference library. Hence, DNA barcoding has now been broadly applied to various fields such as preserving natural resources, protecting endangered species, controlling agriculture pests, identifying disease vectors, monitoring water quality, authentication of natural health products, identification of medicinal plants, food traceability and commercial-fraud prevention [4, 11].

Despite DNA barcoding has been successfully identified marine ichthyofauna and provided the wealth of DNA barcode information in many places, such as Portugal [12], Canada [6, 13], Australia [5], China [1417], India [18], and Taiwan [19]. However, in Southeast Asian waters, an area exceptionally rich in fish biodiversity with more than 3,400 marine and 1000 freshwater species [20, 21], DNA barcoding studies have just been concentrated recently, such as 18 species in Laguna de Bay [22] and 27 commercially important species of the Serranidae [23] in the Philippines; eight species of the Ariidae [24] and 36 species of the Carangidae in the coast of Malaysia [25]; some species of genera Amphiprion and Premnas in Thailand [26] and 26 species of the Gobiidae in Vietnam [27]. Zhang and Hanner [15] conducted a relatively comprehensive survey of fish biodiversity in the South China Sea (SCS), with 242 species successfully identified by DNA barcoding from 1336 specimens. However, Zhang and Hanner’s survey mainly focused on the southern coast of China and the middle SCS. Peripheral countries such as Vietnam with long coastlines and various coastal landforms, which may contribute massive biodiversity, are not included in their study.

Vietnam’s coast has been drawing much attention due to its unique geographical location and marine environment characteristics which show physical, environmental and current gradients along its extraordinarily long coastline [2831]. The various gradients could be potential influences on biodiversity and structure of fish community in northern, central and southern areas, and variations in fish community could be the biological basis for fishery management in Vietnam. At least 2083 marine fish species belonging to 717 genera and 198 families have been recorded in Vietnam [32, 33], and among which 130 fish species are considered economically important species [34]. Fishery resources play an important role in the economy of Vietnam [35], and the fisheries are a major source of food and have helped coastal communities to maintain their livelihoods and community structure. However, the coastal ichthyofaunae of Vietnam remain unclear due to lack of taxonomists [36] whose training is time-consuming and species delimitation requires painstaking labor. Hence, application of suitable methods for accurate and quick identification of fish will be crucial to assist in managing fisheries for long-term sustainability, and will improve ecosystem research and conservation.

DNA barcoding would be a powerful tool to improve the knowledge of fish fauna in Vietnam. However, there have been only a few DNA barcoding studies of fishes in Vietnam and most of them were focused on a small taxonomic scale and/or samples collected from a limited area [3740]. A comprehensive DNA barcoding study on the coastal fishes of Vietnam is urgently needed, especially before its’ growing fishery industry threatening the diversity. The present study is aimed (1) to DNA barcode fishes collected from Vietnam in order to provide a species list confirmed by COI sequences and (2) to preliminarily examine the faunae among fishes collected from northern, central and southern Vietnam.

Material and methods

Sample collection

Through 2016 to 2017, fish specimens were collected from fish markets at five sampling sites that represented three geographic areas along the coast of Vietnam, including Do Son and Ha Long of Hai Phong–Quang Ninh Provinces (Northern area, abbreviated as N, 14 days of sampling effort), Phan Rang of Ninh Thuan Province (Central area, C, 34 days), and Tac Cau and Ha Tien of Kien Giang Province (Southern area, S, 15 days) (Fig 1). Each specimen was photographed, and fin clips or muscle tissues were cut and stored in 95% ethanol and frozen at − 20°C before DNA extraction. For those large specimens and/or expensive fishes, only tissue samples were preserved. Voucher specimens were fixed with 10% formalin and then transferred to 70% ethanol solution for preservation. Voucher specimens and tissue samples were transported to Taiwan and deposited in the Department of Oceanography, National Sun Yat-sen University, Kaohsiung (DOS).

thumbnail
Fig 1. Three sampling areas along the coast of Vietnam.

N comprises Do Son (20.719536, 106.793336) and Ha Long (20.948599, 107.083191) of Hai Phong–Quang Ninh Provinces; C comprises Phan Rang (11.591625, 109.046675; 11.546292, 109.025067) of Ninh Thuan Province; S comprises Tac Cau (9.876155, 105.124151) and Ha Tien (10.378664, 104.485162) of Kien Giang Province. GPS coordinates are for the major fish markets visited every day.

https://doi.org/10.1371/journal.pone.0222631.g001

DNA sequencing

An initial species checklist was built based on morphological identification of photographs and specimens. One specimen of each species from each geographic area depending on morphological-based identification was selected to extract DNA for further molecular analysis. Crude DNA was extracted from the tissue samples using a GeneMark DNA Purification Kit (GMbiolab, Taichung, Taiwan). Polymerase chain reaction (PCR) was performed to amplify mitochondrial cytochrome c oxidase subunit I (COI) using primers designated by Ward, Zemlak (5]: FishF1 (5’-TCAACCAACCACAAAGACATTGGCAC-3’); FishF2 (5’-TCGACTAATCATAAAGATATCGGCAC-3’) and FishR1 (5’- TAGACTTCTGGGTGGCCAAAGAATCA-3’); FishR2 (5’-ACTTCAGGGTGACCGAAGAATCAGAA-3’). PCRs were run with a reaction volume of 25 μL containing 3 μL of template DNA (50 ng μL-1), 3 μL of 10x buffer, 2 μL of dNTPs (2.5 mM), 1.2 μL of each primer (10 μM), 0.13 μL of ProTaq Plus polymerase (Protech, Taipei) and 14.47 μL of deionized water. PCR was performed with an initial denaturation at 94°C for 5 min, followed by 35 cycles of amplification (denaturing at 94°C for 30 s, annealing at 50°C for 30 s, and extension at 70°C for 1 min) and a final extension at 72°C for 8 min. PCR products were checked by electrophoresis with 1.5% agarose gel and then purified by the SAP-Exo purification kit (Jena Bioscience, Jena) according to the manufacturer’s protocols. Sequences were generated by an ABI 3730 automated sequencer and edited manually using MEGA version X [41]. All sequences were translated into amino acids to confirm the effectiveness of the sequences and to detect the presence of nuclear DNA pseudogenes, insertions, deletions, or stop codons. All sequences were submitted to GenBank, and their accession numbers were listed in S1 Table while catalog numbers of tissue samples were included in the description of each sequence in the NCBI webpage. Voucher specimens and tissue samples have identical catalog numbers.

Data analyses

DNA barcoding was conducted using COI sequences blasted in the GenBank and BOLD online databases with a threshold of similarity higher than 98% [42]. Specimens would be re-examined if any conflict exists between morphological and molecular identification. The final species list was compared with four literatures of Vietnamese fishes as well as the Fishbase online database to determine the number of new record species [4346]. After blast in online database, COI sequences were aligned and trimmed to the same length using MEGA X for subsequent analysis. The Kimura-2-Parameter (K2P) genetic distances [47] were calculated at different taxonomic levels, including intraspecific, interspecific within the same genus, inter-generic within the same family, inter-family within the same order, and inter-order in the same class. For intraspecific level, species that only contain one sequence were excluded from the distance analysis; for intra-generic level, genera that contain only one species were excluded, and this criterion was applied for the higher levels of genetic distance analysis. Genetic distances at all levels were calculated disregarding the source of the specimens. The values may contain intra- and/or inter-regional divergences. According to Ward [48], species that shows intraspecific variations above 2% were treated as different species. Higher taxonomy, order and family, in the present study followed Betancur-R et al. [49]. A neighbour-joining (NJ) tree of all analyzed COI sequences was built using Tamura-Nei + Γ model recommended by MEGA X [50] with 1,000 bootstrap replications [51]. Both model test and tree construction were conducted by MEGA X.

Results

A total of 3,638 tissue samples were collected from Vietnam, among which 765 COI sequences were obtained (Table 1). Base on morphological and molecular identifications, these samples represented 458 species of 273 genera, 113 families and 43 orders (S1 Table). Among the 458 species, 59 fishes were new records in Vietnam, sequences of six species were new to the Genbank and BOLD online databases, and only 32 species cannot be annotated to species level because difficulty in morphological identifications and their genetic distances to most similar sequences were more than 2% (S1 Table). The length of all barcode sequences ranged from 313 to 683 bp, with an average of 588 bp and 81% of sequences were longer than 550 bp. No stop codon, insertion, or deletion was observed in any of the obtained sequences.

thumbnail
Table 1. Summary of specimens, sequences and species number in each area in this study.

https://doi.org/10.1371/journal.pone.0222631.t001

Thirty specimens and their COI sequences from Nguyen et al. [27] were included and re-identified in this study. In total, there were 795 COI sequences of 478 species, 284 genera, 113 families and 43 orders for the tree reconstruction (S1 Table and Fig 2). After alignment, all sequences were trimmed to 500 bp and those shorted than 500 bp were abandoned. In total, 707 sequences, including those not able to be determined to species level, were remained for subsequent analyses (444 species of 264 genera, 105 families and 42 orders).

thumbnail
Fig 2. Variation of taxonomic levels in each area.

Data include records from [27].

https://doi.org/10.1371/journal.pone.0222631.g002

The K2P genetic distances within each taxonomic level are summarized in Table 2. The average intraspecific K2P genetic distance was 0.34 ± 0.03%. Intra-generic, -families, -orders, and -class genetic distances were 12.14 ± 0.61%, 17.39 ± 0.63%, 21.42 ± 0.72% and 24.80 ± 1.70%, respectively (Table 2). An increase in genetic variation at higher taxonomic levels was observed, but the rate of increase declined with the higher taxonomic categories (Fig 3). K2P genetic distances of 10 species (Platycephalus indicus, Fistularia petimba, Plotosus lineatus, Uroconger lepturus, Trypauchen vagina, Scolopsis vosmeri, Scatophagus argus, Seriolina nigrofasciata, Ablennes hians, Decapterus maruadsi) were greater than 2% (Table 3) and their sequences were divided into two clades. However, despite the high genetic distance, all their sequences were identified as one same species in BLAST results. In order to validate the sequences of these 10 species, additional COI sequences of each species were obtained from BOLD and analysed together with the sequences generated in the present study. The NJ trees with selected outgroups (S1 Fig) showed that each of those species can be separated into different clusters.

thumbnail
Fig 3. Distribution of K2P distances (percent) for mitochondrial COI (707 sequences, 500 bp) at different taxonomic levels.

Data include records from [27].

https://doi.org/10.1371/journal.pone.0222631.g003

thumbnail
Table 2. Summary of genetic divergences (K2P percent) within various taxonomic levels based on 707 COI sequences.

https://doi.org/10.1371/journal.pone.0222631.t002

thumbnail
Table 3. Summary of K2P genetic divergences of 10 species with intraspecific distances more than 2%.

https://doi.org/10.1371/journal.pone.0222631.t003

The NJ tree based on 707 COI fragments showed that all sequences from the same species were monophyletic (10 aforementioned species with large intraspecific distances were treated as different species; S2 Fig). At the intra-generic level, however, 36 genera were not monophyletic, including Acanthocepola, Acentrogobius, Alepes, Ambassis, Arnoglossus, Callionymus, Carangoides, Chaetodon, Chelon, Cynoglossus, Evynnis, Gerres, Gymnothorax, Hyporhamphus, Leiognathus, Lethrinus, Lutjanus, Ostorhinchus, Paramonacanthus, Parupeneus, Pentapodus, Priacanthus, Pseudorhombus, Rogadius, Sardinella, Saurida, Scorpaenopsis, Secutor, Sillago, Sphyraena, Stolephorus, Strongylura, Synodus, Taeniamia, Terapon, Zebrias.

Only 16 out of 478 species were collected from all three areas (S1 Table). Numbers of common species are much higher between areas, with 45, 48 and 40 for N-C, C-S and N-S, respectively (Table 4). The numbers of common species show a trend with decreased distances.

Discussion

As in previous research [5, 6, 1219, 2226, 52], DNA barcoding analysis based on the COI gene accompanying morphological identifications was able to identify most fishes to the species level (446 out of 478) in the coast of Vietnam with limited similarity in fish fauna among our sampling areas (S1 Table). A reliable DNA barcode reference library for the coastal fish in Vietnam was established, and it would facilitate future studies in taxonomy, conservation, fishery management and so on by linking barcode sequences with carefully identified voucher specimens. However, since 32 species cannot be annotated to species level, it implies that more efforts are needed to sequence COI fragments in order to enrich the databases.

Although it is highly controversial [53], the distance-based technique remains as the standard approach in DNA barcoding [54]. In this study, the K2P model was used to ensure consistency and comparability with other barcoding studies. These results indicate that using COI gene sequences as DNA barcodes to discriminate fish species in Vietnam is feasible. Increasing average genetic distance values were obtained at higher taxonomic levels in this study (Fig 3). The average genetic distances between individuals within species, genera, families and orders were 0.34%, 12.14%, 17.39% and 21.42%, respectively, consistent with the patterns observed in other fish barcoding studies. For example, the K2P values within species, genera, families, and orders were 0.39%, 9.93%, 15.46%, and 22.18% in Australia [5]; 0.30%, 6.60%, 9.91%, and 16.00% in India [18]; 0.18%, 13.55%, 19.65%, and 24.05% in SCS [14]; 0.21%, 5.28%, 21.3%, and 23.63% in Rongcheng Bay, China [17]. The intraspecific divergence (about 0.34%) was much smaller than 2% in Vietnam, but the much lower intraspecific divergence is also observed in previous studies (0.3% in Lakra et al. [18]; 0.319% in Zhang, 2011 [55]; 0.18% in Zhang and Hanner [15]; 0.37% in Jaafar et al. [25]; 0.21% in Bingpeng et al. [16]). It is worth mentioning that sequences used for analyses in this study only with a length of 500 bp while all references mentioned above used the lengths between 605–655 bp, indicating the sufficient effectiveness even though using shorter sequences.

In a total of 478 fish species obtained from 795 COI sequences, there were 139, 285 and 171 species in northern, central and southern areas, respectively. In addition to species diversity, highest numbers of all taxonomic levels were observed in central Vietnam (Fig 2). The highest diversity may be due to the longer sampling period (34 vs. 14 and 15 days in the north and south, respectively), and sampling seasons may also be a concern (Dec., June, and Sep. for north, central and south, respectively). In addition, habitat diversity may also contribute to the fish fauna biodiversity. Most substrates in the northern and southern Vietnam are sandy and muddy while the central part of Vietnam contains various landforms of rocks and coral reefs in addition to sandy substrates. Habitat diversity may result in the higher fish biodiversity in central Vietnam.

There was a huge difference in species composition among geographic areas, only 16 out of 478 species were collected from all three areas (S1 Table). Numbers of common species are much higher between areas, with 45, 48 and 40 for N-C, C-S and N-S, respectively (Table 4). Since the geological distance of C-S is the shortest among three areas, the numbers of common species show a trend with decreased distances. However, variation in physical and environmental gradients may also be associated with the trend along the coastline over 3,260 km [28]. Our northern and southern sampling areas located at nearly the distal ends of the Vietnamese coastline, and it has been reported that the coast of Vietnam exhibits a remarkable environmental shift from north to south. The environmental gradients include average sea surface temperature (23–26–27°C or 25–28–30°C depending on studies); salinity (31–34–31 psu), and dissolved oxygen (4.5–4.0–3.8 ml L-1) [29, 30]. These physical factors are definitely important influences shaping the fish faunae. Moreover, the prevailing currents along the coast of Vietnam are also divided into three segments, approximately corresponding to our sampling sites [31]. The middle segment off the central Vietnamese coast may have a meandering northward coastal current and cyclonic gyre in different seasons while the north and south segments are rather stable in current directions [31]. Since currents may influence pelagic larval dispersal, complex currents might raise biodiversity of coastal species as shown in our central area.

Only 98 out of 478 species (20.50%) in this study have been repeatedly observed in Zhang and Hanner [15]. Numbers of common species between each area in this study and Zhang and Hanner [15] are 58 out of 139 (41.73%) in N, 46 out of 285 (16.14%) in C, and 46 out of 171 (26.90%) in S. The species composition in N is most similar to Zhang and Hanner’s study in SCS. Most sampling of Zhang and Hanner’s survey were carried out along the southern coast of China, and results in higher number of common species due to closer geological distance as shown in C-S. Besides, our sampling sites in each area were mainly focused on fish landing sites that most of fishes were obtained from the nearshore area. This may partially explain why species composition of N is rather similar to that of Zhang and Hanner [15] in SCS.

In the NJ tree of 707 DNA barcode sequences, monophyly occurred in every single species but not in each genus. This result indicates that the DNA barcode is reliable in diagnosis at species level but the resolution is not good enough for determining higher taxonomic levels (e.g. genus and family). There are 36 genera not monophyletic in the NJ tree. We used counterpart sequences in the same family from BOLD to reconstruct additional trees for each non-monophyletic genus (data not shown). Most non-monophyletic genera can be monophyletic in additional trees that only contain sequences from the same family except for eight genera (Acentrogobius, Chaetodon, Gymnothorax, Leiognathus, Lutjanus, Secutor, Sillago and Stolephorus), implying the needs of further taxonomic revisions. The genus Acentrogobius of the Gobiidae is paraphyletic nested by Parachaeturichthys. The paraphyly is consistent with the COI tree in Lavery Clements [56], and further studies may be necessary to clarify their relationship. Chaetodon of the Chaetodontidae is paraphyletic nested by Parachaetodon, which has been considered a synonym of Chaetodon [57]. The largest genus Gymnothorax in the Muraenidae is polyphyletic, consistent with the mitochondrial study using 12S and 16S genes [58]. Lutjanus of the Lutjanidae is a paraphyletic group nested by Caesio, Dipterygonotus and Pterocaesio, three genera used to be placed in the Caesionidae but now in the Lutjanidae [49]. The result is consistent with the relationships shown in Guo [59] based on complete mitochondrial genomes. Leiognathus and Secutor of the Leiognathidae are polyphyletic in our additional tree. Leiognathus has been revealed as a paraphyletic genus in Ikejima [60]. However, Secutor was monophyletic according to Zhang and Hanner [15]. The ambiguous relationship among genera in the Leiognathidae may be attributed to massive misidentifications of previously uploaded sequences [15]. In the Sillaginidae, Sillaginodes is nested in Sillago. However, this result is opposite to the morphological tree reconstructed by Kaga [61]. Similar to the Leiognathidae, sequences of sillaginid fishes may contain misidentifications due to their highly similar morphologies. Stolephorus of the Engraulidae is polyphyletic in this study but monophyletic in Qu [62] reconstructed using complete mitochondrial genomes. The 500 bp length sequences may not provide sufficient information to separate Stolephorus from its counterparts.

Conclusions

It is confirmed that DNA barcoding is efficient and reliable tools for identification of coastal fishes at the species level in Vietnam. In this study, 765 specimens from 458 different species of 273 genera, 113 families in 43 orders were DNA-barcoded. A total of 59 species were newly recorded in Vietnam. The average genetic distance using K2P analysis of individuals within species, genera, families, orders and classes were 0.34%, 12.14%, 17.39%, 21.42% and 24.80%, respectively. Further, establishment of a reliable COI barcode database of coastal fish fauna in Vietnam may serve as a reference library for accurate identification of fishes that may facilitate ichthyological research, including taxonomy, fishery, biodiversity management and so on. Presence of high degree of pair-wise intra-specific divergence among the individuals of some species was revealed, suggesting the presence of cryptic species and advocating the need for more comprehensive studies on fish fauna of Vietnam.

Supporting information

S1 Fig. Neighbor-joining trees based on COI sequences from 10 species (showed in Table 3) with intraspecific K2P distances greater than 2%.

Best chosen models were used with 1,000 bootstrap replications for tree reconstructions.

https://doi.org/10.1371/journal.pone.0222631.s001

(PDF)

S2 Fig. Neighbor-joining tree based on 707 COI sequences using Tamura-Nei + Γ model with 1,000 bootstrap replications.

*and** denote 10 species with intraspecific distances more than 2%.

https://doi.org/10.1371/journal.pone.0222631.s002

(PDF)

S1 Table. The list of fish species DNA barcoded in this study.

NR–New recoded species in Vietnam; ND–Not yet determined; N, C, S–northern, central, and southern areas in Vietnam; *and** denote 10 species with intraspecific distances more than 2% (in Table 3); species in bold were distributed in all areas; ● species recorded in this study; ○ species recorded from Nguyen et al. [27]; numbers in parentheses represent total species excluding those recorded in Nguyen et al. [27].

https://doi.org/10.1371/journal.pone.0222631.s003

(DOC)

Acknowledgments

We are grateful to Claudio Oliveira and another anonymous referee for their constructive comments; to Jie-Cheng Chuang and Ming-hsien Chen of Sheng Long Biotech for their kind assistance during the survey in Ninh Thuan; to Zen-Wei Lin, Han-Ting Huang, Chih-Yen Tung, Po-Hao Chiu, Ching-Ju Chen, Jiehong Wei and You-Chen Liu for collecting specimens; to Wei-Cheng Jhuang, Sung-Yu Tsai and Wei-Che Hung for assistance in fish identification; and to Sih-Yu Chen for map drawing.

References

  1. 1. Hebert PD, Cywinska A, Ball SL. Biological identifications through DNA barcodes. Proceedings of the Royal Society of London B: Biological Sciences. 2003;270(1512):313–21.
  2. 2. Frézal L, Leblois R. Four years of DNA barcoding: current advances and prospects. Infection, Genetics and Evolution. 2008;8(5):727–36. pmid:18573351
  3. 3. Leray M, Knowlton N. DNA barcoding and metabarcoding of standardized samples reveal patterns of marine benthic diversity. Proceedings of the National Academy of Sciences. 2015;112(7):2076–81.
  4. 4. Kaur S. DNA Barcoding and Its Applications. International Journal of Engineering Research and General Science. 2015;3(2):602–4.
  5. 5. Ward RD, Zemlak TS, Innes BH, Last PR, Hebert PD. DNA barcoding Australia's fish species. Philosophical Transactions of the Royal Society of London B: Biological Sciences. 2005;360(1462):1847–57. pmid:16214743
  6. 6. Hubert N, Hanner R, Holm E, Mandrak NE, Taylor E, Burridge M, et al. Identifying Canadian freshwater fishes through DNA barcodes. PLoS one. 2008;3(6):e2490. pmid:22423312
  7. 7. Valdez-Moreno M, Ivanova N, Elías-Gutiérrez M, Contreras-Balderas S, Hebert P. Probing diversity in freshwater fishes from Mexico and Guatemala with DNA barcodes. Journal of Fish Biology. 2009;74(2):377–402. pmid:20735566
  8. 8. Ward RD, Hanner R, Hebert PD. The campaign to DNA barcode all fishes, FISH-BOL. Journal of fish biology. 2009;74(2):329–56. pmid:20735564
  9. 9. Becker S, Hanner R, Steinke D. Five years of FISH-BOL: brief status report. Mitochondrial DNA. 2011;22(sup1):3–9.
  10. 10. Mecklenburg CW, Møller PR, Steinke D. Biodiversity of arctic marine fishes: taxonomy and zoogeography. Marine Biodiversity. 2011;41(1):109–40.
  11. 11. Galimberti A, De Mattia F, Losa A, Bruni I, Federici S, Casiraghi M, et al. DNA barcoding as a new tool for food traceability. Food Research International. 2013;50(1):55–63.
  12. 12. Costa FO, Landi M, Martins R, Costa MH, Costa ME, Carneiro M, et al. A ranking system for reference libraries of DNA barcodes: application to marine fish species from Portugal. PLoS One. 2012;7(4):e35858. pmid:22558244
  13. 13. Steinke D, Zemlak TS, Boutillier JA, Hebert PDN. DNA barcoding of Pacific Canada’s fishes. Marine Biology. 2009;156(12):2641.
  14. 14. Zhang J. Species identification of marine fishes in China with DNA barcoding. Evidence-Based Complementary and Alternative Medicine. 2011;2011:10. http://dx.doi.org/10.1155/2011/978253.
  15. 15. Zhang J, Hanner R. Molecular approach to the identification of fish in the South China Sea. PLoS One. 2012;7(2):e30621. pmid:22363454
  16. 16. Bingpeng X, Heshan L, Zhilan Z, Chunguang W, Yanguo W, Jianjun W. DNA barcoding for identification of fish species in the Taiwan Strait. PloS one. 2018;13(6):e0198109. pmid:29856794
  17. 17. Wang L, Wu Z, Liu M, Liu W, Zhao W, Liu H, et al. DNA barcoding of marine fish species from Rongcheng Bay, China. PeerJ. 2018;6:e5013. pmid:29967722
  18. 18. Lakra W, Verma M, Goswami M, Lal KK, Mohindra V, Punia P, et al. DNA barcoding Indian marine fishes. Molecular Ecology Resources. 2011;11(1):60–71. pmid:21429101
  19. 19. Chang C-H, Shao K-T, Lin H-Y, Chiu Y-C, Lee M-Y, Liu S-H, et al. DNA barcodes of the native ray-finned fishes in Taiwan2017. 796–805 p.
  20. 20. Alien GR, Amaoka K, Anderson WD Jr, Bellwood DR, Bohlke EB, Bradbury MG, et al. A checklist of the fishes of the South China Sea. The Raffles Bulletin of Zoology. 2000;(8):569–667.
  21. 21. Zakaria-Ismail M. Zoogeography and biodiversity of the freshwater fishes of Southeast Asia. Ecology and conservation of Southeast Asian marine and freshwater environments including wetlands: Springer; 1994. p. 41–8.
  22. 22. Aquino LMG, Tango JM, Canoy RJC, Fontanilla IKC, Basiao ZU, Ong PS, et al. DNA barcoding of fishes of Laguna de Bay, Philippines. Mitochondrial DNA. 2011;22(4):143–53. pmid:22040082
  23. 23. Alcantara SG, Yambot AV. DNA barcoding of commercially important grouper species (Perciformes, Serranidae) in the Philippines. Mitochondrial DNA Part A. 2016;27(6):3837–45.
  24. 24. Abdullah ML, Nor SM, Naim DM. Use of DNA barcode in the identification of catfishes (Siluriformes: Ariidae) from Malaysia. Biodiversitas Journal of Biological Diversity. 2017;18(4):1358–66.
  25. 25. Jaafar TNAM, Taylor MI, Nor SAM, de Bruyn M, Carvalho GR. DNA barcoding reveals cryptic diversity within commercially exploited Indo-Malay Carangidae (Teleosteii: Perciformes). PLoS One. 2012;7(11):e49623. pmid:23209586
  26. 26. Ayudhaya PT, Muangmai N, Banjongsat N, Singchat W, Janekitkarn S, Peyachoknagul S, et al. Unveiling cryptic diversity of the anemonefish genera Amphiprion and Premnas (Perciformes: Pomacentridae) in Thailand with mitochondrial DNA barcodes. Agriculture and Natural Resources. 2017;51(3):198–205.
  27. 27. Nguyen ML, Nguyen VQ, Pham VC, Dao HL, Dinh VN, Dam TL. DNA barcoding application of mitochondrial COI gene to identify some fish of family Gobiidae in Vietnam. Vietnam Journal of Marine Science and Technology. 2018;18(4):443–51.
  28. 28. Pham TT, Furukawa M. Impact of Sea Level Rise on Coastal Zone of Vietnam. Bull Faculty Sci Univ Ryukyus. 2007;84.
  29. 29. Rojana-anawat P, Pradit S, Sukramongkol N, Siriraksophon S, editors. Temperature, salinity, dissolved oxygen and water masses of Vietnamese waters. Proceedings of the SEAFDEC Seminar on Fishery Resources in the South China Sea, Area IV: Vietnamese Waters; 2001; Thailand: Southeast Asian Fisheries Development Center.
  30. 30. Ngo-Duc T. Climate change in the coastal regions of Vietnam. Coastal Disasters and Climate Change in Vietnam: Elsevier; 2014. p. 175–98.
  31. 31. Fang G, Wang G, Fang Y, Fang W. A review on the South China Sea western boundary current. Acta Oceanologica Sinica. 2012;31(5):1–10.
  32. 32. Nguyen KH. Fish fauna of Vietnam. Collection of Marine Research Works 1995;IV:121–6.
  33. 33. Bui DC, Chu TV, Nguyen HD. Marine fisheries resources—Basic development of marine fisheries industry in Vietnam. In Proceedings of Marine Fisheries Research, Research Institute for Marine Fisheries, Vietnam Ministry of Fisheries, 560 p. 2001;II:200–10.
  34. 34. Dao MS, Pham T, editors. Management of coastal fisheries in Vietnam. Assessment, Management and Future Directions for Coastal Fisheries in Asian Countries: WorldFish Center Conference Proceedings; 2003.
  35. 35. Nguyen H, Hien P, Lebailly P. Vietnam’s fisheries and aquaculture development’s policy: are exports performance targets sustainable. 2017.
  36. 36. Lutaenko KA. Proceedings of the Workshop Coastal Marine Biodiversity and Bioresources of Vietnam and Adjacent Areas to the South China Sea, NhaTrang, Vietnam, November 24–25, 2011. 2011:Vladivostok-NhaTrang: Dalnauka, 2011. 123 pp.
  37. 37. Duong TY, Nguyen K, Bui SN, Nguyen VT, Nguyen BL, Tran DD. DNA barcodes and morphology of P. krempfi, P. mekongensis, and P. elongatus. Vietnam Journal of Biotechnology. 2016;14(1):29–37.
  38. 38. Vu DHQ, Le TTP, Tran TO, Le TT, Dang TB. Phylogenetic relationships of freshwater fish in Vietnamese Mekong. International Conference on Biological, Environment and Food Engineering (BEFE-2015) May 15–16, 2015 Singapore. 2015:58–61. http://dx.doi.org/10.15242/IICBE.C0515014.
  39. 39. Tran TO, Vu DHQ, Dang TB. Phylogenetic Relationships of Emperors (Lethrinidae) and Snappers (Lutjanidae) in Vietnam based on Mitochondrial DNA Sequences. International Conference on Biological, Environment and Food Engineering (BEFE-2015) May 15–16, 2015 Singapore. 2015:62–5. http://dx.doi.org/10.15242/IICBE.C0515016.
  40. 40. Tran TVT, Phan KL, Durand J-D. Diversity and distribution of cryptic species within the Mugil cephalus species complex in Vietnam. Mitochondrial DNA Part A. 2017;28(4):493–501.
  41. 41. Kumar S, Stecher G, Li M, Knyaz C, Tamura K. MEGA X: molecular evolutionary genetics analysis across computing platforms. Molecular biology and evolution. 2018;35(6):1547–9. pmid:29722887
  42. 42. Arroyave J, Stiassny ML. DNA barcoding reveals novel insights into pterygophagy and prey selection in distichodontid fishes (Characiformes: Distichodontidae). Ecology and evolution. 2014;4(23):4534–42. pmid:25512849
  43. 43. Tran D, Nguyen NT. List of marine fish in Vietnam. Marine research collection 1985:19–45.
  44. 44. Nguyen KH. Fauna of Vietnam: Hanoi Science and Technology publisher; 2001. 324 p.
  45. 45. Nguyen NT, Nguyen VQ. Biodiversity and living resources of the coral reef fish in Vietnam marine waters: Hanoi Science and Technology publisher; 2005. 120 p.
  46. 46. Seishi K, Hisashi I, Nguyen VQ, Pham TD. Fishes of Ha Long Bay: Fisheries Research Laboratory, Mie University; 2018. 314 p.
  47. 47. Kimura M. A simple method for estimating evolutionary rates of base substitutions through comparative studies of nucleotide sequences. Journal of molecular evolution. 1980;16(2):111–20. pmid:7463489
  48. 48. Ward RD. DNA barcode divergence among species and genera of birds and fishes. Molecular ecology resources. 2009;9(4):1077–85. pmid:21564845
  49. 49. Betancur-R R, Wiley EO, Arratia G, Acero A, Bailly N, Miya M, et al. Phylogenetic classification of bony fishes. BMC Evolutionary Biology. 2017;17(1):162. pmid:28683774
  50. 50. Tamura K, Nei M. Estimation of the number of nucleotide substitutions in the control region of mitochondrial DNA in humans and chimpanzees. Molecular biology and evolution. 1993;10(3):512–26. pmid:8336541
  51. 51. Felsenstein J. Confidence limits on phylogenies: an approach using the bootstrap. Evolution. 1985;39(4):783–91. pmid:28561359
  52. 52. Landi M, Dimech M, Arculeo M, Biondo G, Martins R, Carneiro M, et al. DNA barcoding for species assignment: the case of Mediterranean marine fishes. PLoS One. 2014;9(9):e106135. pmid:25222272
  53. 53. Srivathsan A, Meier R. On the inappropriate use of Kimura-2-parameter (K2P) divergences in the DNA-barcoding literature. Cladistics. 2012;28(2):190–4.
  54. 54. Reid B, Le M, McCord W, Iverson J, Georges A, Bergmann T, et al. Comparing and combining distance-based and character-based approaches for barcoding turtles. Molecular Ecology Resources. 2011;11(6):956–67. pmid:21635698
  55. 55. Zhang J. Species identification of marine fishes in China with DNA barcoding. Evidence-Based Complementary and Alternative Medicine. 2011;2011.
  56. 56. Lavery S, Clements KD, Hickey A. Molecular identification of cryptogenic/invasive gobies in New Zealand: Ministry of Fisheries; 2006.
  57. 57. Bellwood D, Klanten S, Cowman P, Pratchett M, Konow N, Van Herwerden L. Evolutionary history of the butterflyfishes (f: Chaetodontidae) and the rise of coral feeding fishes. Journal of evolutionary biology. 2010;23(2):335–49. pmid:20487131
  58. 58. Tang KL, Fielitz C. Phylogeny of moray eels (Anguilliformes: Muraenidae), with a revised classification of true eels (Teleostei: Elopomorpha: Anguilliformes). Mitochondrial DNA. 2013;24(1):55–66. pmid:22967094
  59. 59. Guo Y, Bai Q, Yan T, Wang Z, Liu C. Mitogenomes of genus Pristipomoides, Lutjanus and Pterocaesio confirm Caesionidae nests in Lutjanidae. Mitochondrial DNA Part A. 2016;27(3):2198–9.
  60. 60. Ikejima K, Ishiguro NB, Wada M, Kita-Tsukamoto K, Nishida M. Molecular phylogeny and possible scenario of ponyfish (Perciformes: Leiognathidae) evolution. Molecular Phylogenetics and Evolution. 2004;31(3):904–9. pmid:15120388
  61. 61. Kaga T. Phylogenetic systematics of the family Sillaginidae (Percomorpha: order Perciformes). Zootaxa. 2013;3642(1):1–105.
  62. 62. Qu L, Ma C, Ma H, Dong Y, Wang W, Ren G, et al. The complete mitochondrial genome of Stolephorus commersonii. Mitochondrial DNA Part B. 2017;2(2):573–4.