Assessing larval fish diversity and conservation needs in the Luzon strait using DNA barcoding

Introduction

Fishes account for approximately 50% of marine vertebrate species, serving as a key component of biodiversity in global marine ecosystems. Larval fish can provide the fundamental source for the replenishment and sustainable utilization of fish stocks (Panprommin et al., 2020; Panprommin et al., 2021; Huang et al., 2022). Accurate identification and classification of larval fish species are crucial for effective fisheries management, as it enables the estimation of parent stocks and future recruitment, as well as the prediction of important spawning and nursing grounds (Baumgartner et al., 2004; Xing et al., 2022). Traditional morphological identification methods based on morphological features such as spination, body shape, fin placement, and photophores often yield misidentification due to distinct characteristics between larval and adult fish and developmental variations (Azmir et al., 2017; Sheraliev and Peng, 2021). Furthermore, limited diagnostic characteristics hinder precise identification, especially for related taxa (Valdez-Moreno et al., 2010; Kwun, 2018). Moreover, the lack of larval fish taxonomists profoundly influences accurate identifications based on morphology (Azmir et al., 2017). Therefore, the inherent limitations of morphology-based identification systems underscore the necessity for accurate and efficient approaches to identify larval fish species (Zhang and Hanner, 2012).

Molecular methods, such as DNA barcoding, have become effective tools for accurate and rapid species determination, including fish (Zhan et al., 2008; Hubert et al., 2010; Hubert et al., 2015; Panprommin et al., 2021; Huang et al., 2022). DNA barcoding utilizes a standardized DNA region, such as the mitochondrial cytochrome c oxidase subunit I (COI) gene, to identify species based on their unique genetic sequences. It has been successfully applied to identify incomplete or immature fish samples, such as fish eggs (Hou et al., 2022), larval fish (Azmir et al., 2017), and seafood products (Shehata et al., 2018). Compared to traditional morphological identification methods for larval fish species, DNA barcoding provides advantages in terms of efficiency, repeatability, and cost-effectiveness (Azmir et al., 2017). Additionally, online DNA barcode databases, such as the Barcode of Life Data System (BOLD) and the National Center for Biotechnology Information (NCBI), offer valuable resources for species annotation. These databases provide extensive reference sequences from diverse fish species, encompassing different taxonomic groups and geographic regions, facilitating reliable comparison and identification of unknown samples. These databases also provide additional information on geographical distribution, ecological traits, and morphological characteristics, which would enhance our understanding of the biology and ecology of the identified species. Environmental DNA (eDNA) metabarcoding has recently emerged as a revolutionary tool for monitoring fish biodiversity (Yao et al., 2022; Zhang et al., 2022). Nonetheless, a critical challenge arises when utilizing public databases in fish eDNA metabarcoding, leading to increased potential for ambiguous species identification or false positives. In contrast, the use of local databases can significantly enhance the efficiency of species identification, ensuring identification at the species level (Yao et al., 2022). Thus, conducting DNA barcoding analysis of larval fish in specific sea areas provides valuable insights into fish biodiversity, population dynamics, and ecological interactions in marine ecosystems, contributing to the exploration, management, and conservation of fisheries resources (Hubert et al., 2015; Yi et al., 2017).

The Luzon Strait, located between the Luzon Island in the Philippines and Taiwan, serves as a vital channel connecting the South China Sea to Western Pacific Ocean. Recognized as a hotspot for marine biodiversity, this strait holds significant ecological importance. It encompasses diverse habitats such as open ocean waters, coral reefs, seagrass beds, and mangrove forests (Gu et al., 2019). The Luzon Strait is home to a variety of coral reef-associated fish species, including snappers (e.g., Lutjanus erythropterus and L. madras), wrasses (e.g., Thalassoma amblycephalum and T. quinquevittatum), and groupers (e.g., Epinephelus lanceolatus and E. fasciatus) (Huang et al., 2015). It also supports commercially important fish species, such as skipjack tuna (Katsuwonus pelamis), yellowfin tuna (Thunnus albacares), bigeye tuna (Thunnus obesus), as well as various groupers and mackerels, highlighting its significance for fisheries (Zhang and Hanner, 2012; Xu et al., 2021). However, limited information is available regarding the larval fish diversity in this region. The Luzon Strait experiences high oceanographic variability influenced by multiple oceanic currents, eddies, and fronts. Particularly, the Kuroshio current originating from the Western North Pacific brings warm, salty, and nutrient-rich waters into the strait, significantly affecting its hydrological dynamics and fish diversity (Gu et al., 2019; Zhang et al., 2022). However, the region is facing ecological and environmental challenges due to frequent anthropogenic disturbances such as overfishing and pollution (Chang et al., 2017; Hou et al., 2018). These activities have altered fish community structures and resulted in the decline of fish resources (Olden et al., 2007; Chang et al., 2017). The destruction of habitats, including fish spawning and nursing grounds, further threatens the sustainability of marine fish populations and ecosystems (Olden et al., 2007). Therefore, monitoring the composition, diversity, and dynamics of larval fishes in the Luzon Strait and its surrounding waters is crucial for identifying potential threats and implementing effective conservation strategies.

In this study, we utilized DNA barcoding to identify larval fish species collected from stations in the Luzon Strait and adjacent waters. The aim of this study is to assess the species composition, diversity, and geographical distribution of the larval fish community in this region. Additionally, we investigated the habitat type, human use, and IUCN conservation status of each larval fish species. Our findings are expected to complement the current fish database of the Luzon Strait and surrounding waters, and further contribute to identifying important spawning and nursing grounds of fish populations, thereby supporting sustainable fisheries resource management as well as implementing conservation programs.

Materials and methods

Ethics statement

The requirement of ethical approval was waived by the Animal Experiment Administration Committee of South China Sea Ecological Center, Ministry of Natural Resources of the People’s Republic of China. The studies were conducted in accordance with the local legislation and institutional requirements.

Sample collection

In May 2016, we collected larval fish specimens from 15 stations in the Luzon Strait and adjacent waters (Figure 1; Table S1). To capture a wide variety of fish species, we employed both vertical and horizontal trawls to collect larval fish individuals by using a zooplankton net with a length of 2.8 m, bore diameter of 80 cm, and mesh size of 0.505 mm. The nets were flushed repeatedly from top to bottom to collect the larval fish individuals. For vertical trawls, nets were lowered down to approximately the seabed and dragged with a speed of 0.5-0.8 m/s. For horizontal trawls, nets were towed for 15 min with a speed of 1.0-2.0 knots/h at each station. After sampling, all specimens were immediately preserved in 95% ethanol and transported to the laboratory for subsequent analyses.

FIGURE 1

Figure 1 The map of 15 stations in the Luzon Strait and its surrounding areas where larval fish samples were collected. The path of the Kuroshio Current was approximately shown by light blue shading regions. The colored stations were corresponding to their assigned clusters as shown in Figure 6.

DNA extraction, PCR amplification and sequencing

The genomic DNA of each larval fish individual was extracted from the muscle tissues using the DNeasy Blood and Tissue Kit (QIAGEN), and the quantity and quality of extracted DNA were subsequently measured using a NanoDrop 2000 spectrophotometer (Thermo Fisher Scientific, USA). We chose the cytochrome c oxidase subunit I (COI) mitochondrial gene for fish DNA barcoding. Sequences of approximately 650 bp within the 5’ region of the COI gene were amplified using the universal primers FishF1 [5′-TCAACCAACCACAAAGACATTGGCAC-3′] and FishR1 [5′-TAGACTTCTGGGTGGCCAAAGAATCA-3′] (Ward et al., 2005). Polymerase chain reaction (PCR) amplifications were performed in a volume of 25 μL, containing 2.5 μL 10 X PCR buffer⁺, 2.0 μL of dNTPs (10 mM each), 1.0 μL of each primer (0.01 mM), 0.5 μL of ExTaq polymerase, 1.0 μL of template DNA (approximately 100 ng), and 18 μL ultrapure water. The PCR amplification procedures consisted of an initial denaturation step at 95 °C for 5 min, followed by 35 cycles of 95°C for 30 s, 52°C for 30 s, and 72°C for 60 s, with a final elongation at 72°C for 10 min. All PCR products were visualized on a 1.2% agarose gel, and those successfully amplified PCR products were sequenced bidirectionally with BigDye Terminator Cycle Sequencing chemistry on an Applied Biosystems 3730XL automated DNA analyzer (Applied Biosystems, USA).

Data analysis

The tracer files of sequencing samples were initially checked and edited, and the raw DNA sequences were assembled using the Lasergene package in the DNAStar software. The high-quality sequences were aligned and edited using BioEdit v.7.0.5.3 (Hall, 1999) and MEGA v.7.0 (Kumar et al., 2016). The obtained sequences were searched against two databases, Barcode of Life Data System (BOLD, http://www.boldsystems.org/) and GenBank (https://www.ncbi.nlm.nih.gov/) for species identification. We set a sequence similarity of 98% as a threshold for the species level annotation (Hou et al., 2021). In addition, the sequence with an identity from 95% to 98% matched to the database sequence was classified at the genus level, and a lower similarity than 95% was considered unidentified. The sequences were trimmed to equal length and formatted into a single fasta format using BioEdit, and subsequently arranged using ClustalW in MEGA. Intraspecific and interspecific genetic distances were calculated using the Kimura 2-parameter (K2P) distance model, and the phylogenetic tree was subsequently constructed by the neighbor-joining (NJ) method using the K2P model. The K2P distances and NJ trees were estimated with 1,000 bootstrap replicates in MEGA (Kumar et al., 2016). The NJ tree was visualized, edited, and annotated using interactive tree of life (iTOL, https://itol.embl.de/) (Letunic et al., 2021).

Diversity and geographical distribution of larval fish

We first estimated taxonomic diversity of larval fish using the species richness (i.e., the number of fish), and calculated the occurrence frequency of each larval fish species at different stations using the R package BiodiversityR (Kindt and Coe, 2005). To further measure the differences in larval fish community compositions among different sites, we conducted the non-metric multidimensional scaling (NMDS) analysis based on the Bray-Curtis similarity (Xiong et al., 2017; Xiong et al., 2019; Yang et al., 2019b). In addition, the hydrodynamic process of the Luzon Strait is greatly controlled by the Kuroshio Current (Gu et al., 2019), which would be likely to affect the diversity and community composition of larval fish. We thus divided 15 stations into two clusters based on the path of the Kuroshio current. Cluster I included eight western stations (LS1, LS2, LS3, LS4, LS5, LS13, LS14 and LS15), while cluster II included seven eastern stations (LS6, LS7, LS8, LS9, LS10, LS11 and LS12). The significance of dissimilarity in species composition between two clusters was further assessed using the analysis of similarity (ANOSIM) test (Xiong et al., 2019; Yang et al., 2019a). Both the NMDS analysis and ANOSIM test were conducted using online Omicshare CloudTools (http://www.omicshare.com/tools).

The habitat, human uses, and conservation status

Each larval fish species was searched for the habitat as well as human uses based on the FishBase webpage (http://www.fishbase.org/), which is a global biodiversity information system on fishes. We downloaded the habitat (pelagic, demersal, bathypelagic, reef-associated, etc.) and human use (fisheries, aquaculture, aquariums, etc.) data for each target species. Furthermore, the conservation status of each species was determined according to the International Union for Conservation of Nature (IUCN) webpage (https://www.iucnredlist.org/).

Results

Species identification

A total of 385 larval fish individuals were collected from 15 stations in the Luzon Strait and surrounding waters, and 361 individuals (93.77%) were successfully barcoded. The remaining 24 individuals were excluded from subsequent analyses, as they failed to produce any PCR targeted product or high-quality sequence. After conducting sequence alignment, trimming, and removing noisy sites, we obtained 623 bp DNA barcodes for the 361 fish individuals. Through blast searches against both BOLD and GenBank databases, 347 larval sequences were successfully assigned to 131 taxa to the species level, using a similarity threshold of 98%. Sixteen larval sequences were determined to the genus level with similarities ranging from 95% to 98%, while the remaining seven individuals were unidentified with similarities lower than 95%. Overall, 354 larvae with high similarities belonged to 147 fish taxa, which were categorized into 22 orders, 44 families, and 93 genera (Table S2). Among the identified larval fish, the order Myctophiformes accounted for the highest quantity (31.07%) and variety (26.53%), followed by Stomiiformes (21.19% in quantity and 11.56% in variety) and Perciformes (13.56% in quantity and 19.05% in variety) (Figure 2). Additionally, six orders including Beryciformes, Centrarchiformes, Gadiformes, Gobiiformes, Lophiiformes, and Spariformes were represented by a single species (Figure 2). The family Myctophidae (110 individuals) and the genus Vinciguerria (54 individuals) were the most abundant, accounting for approximately 31% and 15% of all fish individuals, respectively (Table S2).

FIGURE 2

Figure 2 The quantity and variety of species in each order of larval fish identified in this study.

Based on the Kimura 2-parameter (K2P) distance, intraspecific genetic distances ranged from 0 to 3.98%, with an average value of 0.33% (Table S3). In contrast, interspecific genetic distances varied from 5.22% to 39.18%, with an average value of 18.03% (Table S3). These findings highlighted the increased genetic divergence observed between different species compared to within species. To further validate the accuracy of larval fish species identification, we constructed a phylogenetic NJ tree based on DNA barcode sequences. The tree revealed 147 independent lineages (Figure 3), providing additional support for the reliability of species identification achieved through the utilization of the BOLD and NCBI databases.

FIGURE 3

Figure 3 The neighbor-joining (NJ) tree of 147 larval fish using COI sequences.

Species diversity and geographical distribution

The species richness varied across different stations, with the highest richness observed at the station LS12 (55 species), followed by LS8 (29 species) and LS10 (24 species). In contrast, LS1 exhibited the lowest number of species, comprising only two species (Figure 4). Vinciguerria nimbaria emerged as the most frequently occurring larval fish species, which was detected at 11 out of 15 stations. Furthermore, the relative abundance of the dominant orders displayed significant variability across different stations. The order Myctophiformes dominated over half of the stations, including LS2, LS3, LS4, LS5, LS6, LS7, LS11, LS12, LS13, LS14, and LS15. In addition, LS1, LS8, LS9, and LS10 were dominated by the Labriformes (~ 85%), Syngnathiformes (~ 34%), Clupeiformes (~ 58%), and Stomiiformes (~ 33%), respectively (Figure 5).

FIGURE 4

Figure 4 Species richness of larval fish at each station.

FIGURE 5

Figure 5 Taxonomic composition of larval fish at the order level at each station.

Based on the NMDS analysis, we observed the patterns of dissimilarity in larval fish community compositions among different stations. Notably, within the second cluster (II cluster) in the eastern regions, larval fish community compositions displayed higher similarities, as shown by relatively higher Bray-Curtis similarities in the NMDS plot (Figure 6). In contrast, the fish community compositions at stations in the western regions largely varied, with lower similarities and higher dissimilarities. There was no significant difference between the community compositions of the two clusters as indicated by the results of ANOSIM (R = 0.010, P = 0.450).

FIGURE 6

Figure 6 Nonmetric multidimensional scaling ordination (NMDS) based on the Bray-Curtis similarity of larval fish communities at each sampling station.

The habitat, human use, and IUCN status

In terms of the habitat of the identified fish species, they were categorized as follows: 36.64% were bathypelagic, 35.88% were reef-associated, 12.21% were pelagic-oceanic, 7.63% were demersal, 5.34% were benthopelagic, 1.53% were pelagic-neritic, and 0.76% were pelagic (Table S4; Figure 7). When assessing for human use, nearly half of the evaluated fish species were utilized in fisheries (49.62%), while 46.56% were not assessed (NA). Other types of human uses included aquariums (19.08%), gamefish (9.16%), bait (3.82%), and aquaculture (3.05%) (Figure 7). The remaining species (12.98%) were considered of no interest in fisheries. Among the recorded species, 35 were found to serve for multiple purposes, and Seriola dumerili was used for four purposes, including fisheries, aquaculture, gamefish, and aquariums (Table S4). The IUCN status of most species was determined to be least concern (LC, 84.73%), followed by not assessed (NA, 11.45%) and data deficient (DD, 2.29%). In addition, Istiophorus platypterus was designated as vulnerable (VU), and Xiphias gladius was categorized as near threatened (NT).

FIGURE 7

Figure 7 The percentage of fish species according to the habitat, human use and the International Union for Conservation of Nature (IUCN) status, respectively. The information of each species was obtained from the FishBase and IUCN webpages.

Discussion

In this study, we employed DNA barcoding to identify larval fish species collected from the Luzon Strait and adjacent waters. This region is home to diverse fish species, and accurate species identification at larval stages is crucial for understanding their ecology and population dynamics. The accurate determination of larval fish species solely depending on traditional morphology-based methods has proven to be challenging and time-consuming (Neira et al., 2015; Panprommin et al., 2021). However, recent advancements in molecular methods, such as DNA (meta)barcoding, have provided a powerful tool for species identification (Zhan and MacIsaac, 2015; Xiong et al., 2020). The COI gene used for DNA barcoding produced sequences with an average length of approximately 650 bp, which has been widely recognized for its efficacy in identifying larval fish species (Chang et al., 2017; Panprommin et al., 2021). Numerous studies have demonstrated the reliability and applicability of DNA barcoding for identifying larval fish species across various geographic regions, including the marine fish larval community in the South China Sea (Hou et al., 2021), larval fish in mangrove areas from Peninsular Malaysia (Azmir et al., 2017), and larval fish in tropical estuarine waters of Malacca Straits in Malaysia (Chu et al., 2019). The construction of a comprehensive local species reference database is of vital importance to ensure precise identification and complete coverage of fish diversity (Yao et al., 2022). In this study, we identified a set of larval fish in the Luzon Strait and adjacent waters, including 22 orders, 44 families, 93 genera, and 147 species. DNA barcoding allowed to delimit 96.12% of all larval fish sequences to the species level with a 98% similarity threshold. The intraspecific K2P divergence observed in this study was found to be approximately 0.33%, while the interspecific K2P divergence was approximately 18.03%. This represents an approximately 55-fold increase in genetic distance between species when compared to the average intraspecific genetic distances, which is consistent with the DNA barcoding principle that the interspecific divergence largely outscores the intraspecific divergence (Hebert et al., 2003). The constructed NJ tree based on the COI gene sequences also demonstrated that each taxonomic-level species was assigned to an independent lineage. These findings highlight the feasibility of using COI gene sequences as DNA barcodes for discriminating larval fish species in the Luzon Strait and adjacent waters.

We found that 39 of identified larval fish species (26.53%) were categorized into the order of Myctophiformes, which also emerged as the dominant order at more than half of the sampling stations. These findings indicate that this order is likely to be the dominant fish taxon in the Luzon Strait and adjacent waters. While Myctophiformes, specifically the lanternfishes, do not have high fisheries values, they are a significant component of deep-water fauna, which constitute a substantial and play a crucial role in the diet of high economic value species. Despite the ecological significance of Myctophiformes, there remains a lack of molecular studies focusing on their identification and classification (Poulsen et al., 2013; Pappalardo et al., 2015). The first attempt at lanternfish larvae identification using DNA barcoding was conducted in the Sicilian Channel of the Mediterranean Sea, where five Myctophidae species were identified from 20 lanternfish larvae fish individuals (Pappalardo et al., 2015). Similarly, a recent survey of larval fish in the South China Sea implemented DNA barcoding and revealed that the Myctophidae was one of the most abundant families (Hou et al., 2021).

Our study revealed notable variations in fish species diversity and community composition across different sampling stations. Cluster II exhibited higher species richness, with an average of 22 species, when compared with cluster I (average 11 species). Importantly, we identified several economically valuable species at these stations, including Abudefduf vaigiensis, Elagatis bipinnulata, Parupeneus barberinoides, Parupeneus heptacanthus, and Istiophorus platypterus. These findings indicate the ecological significance of these stations as potential spawning or nursing grounds. The presence of economically important species underscores the significance of protecting these areas to ensure the sustainability of fish populations and support fisheries management initiatives. The NMDS analysis in this study demonstrated that the larval fish community compositions in the western regions were more diverse and dissimilar compared to the eastern regions, as indicated by the higher dispersion of stations in the western region. Importantly, the two distinct patterns of larval fish communities observed in this study were consistent with the grouping results based on the Kuroshio current. Specifically, the cluster I, located in the western regions, was situated away from the path of the Kuroshio current. While the cluster II, located in the eastern regions, was situated within the path of the Kuroshio current. Understanding the impacts of oceanic currents can provides valuable insights into predicting and managing the fish diversity and community structures in the region, benefiting both conservation efforts and commercial fisheries (Zhan et al., 2009; Almany et al., 2017). Oceanic currents, such as the Kuroshio current, play a vital role in the transport of fish eggs and larvae, aiding their dispersal until they settle and develop into sedentary adult individuals (Chang et al., 2017; Xu et al., 2021). The Kuroshio current, originating from the western North Pacific, can bring warm, salty, and nutrient-rich waters into the strait, which is expected to cause a significant impact on the fish diversity and community structures in the Luzon Strait and surrounding waters. This inflow of nutrients sustains rich planktonic and benthic communities, providing food sources for a variety of fish species (Xu et al., 2021). The water temperature was different between northern and southern Taiwan, owing to three main ocean currents, including Kuroshio, China Coast and South China Sea (Chang et al., 2017). Consequently, the strait has become an important fishing ground for commercial fisheries. Furthermore, the Kuroshio current plays a significant role in creating different environmental conditions and habitats, which may support specific fish species with adaptive capacities to warm and high-salinity environments (Gu et al., 2019). The presence of certain species in cluster II indicates that larval fish community compositions, particularly in this cluster, are likely influenced by the Kuroshio current. For instance, Ceratoscopelus warmingii, also known as the warming’s lantern fish, is a dominant species among mesopelagic fishes and is known to thrive in warm and high-salinity habitats in tropical and subtropical waters (Kong et al., 2016; Zhang et al., 2021). Its consistent presence in approximately all stations of the cluster II (6 out of 7 stations) indicates that these stations are likely influenced by the Kuroshio current. The current provides warm, high-oxygenation, and high-salinity water conditions that are suitable for this species (Gu et al., 2019). In contrast, the most frequently occurring larval fish in stations of the cluster I was Vinciguerria nimbaria, an oceanic lightfish. This species was reported with rich biological production in low temperature environments, as they were supported by abundant zooplankton (Kim et al., 2023). The current can act as a transport mechanism for fish eggs and larvae, shape the diverse marine ecosystem by providing nutrients, and create distinct environmental conditions and habitats. Overall, these findings support that larval fish community compositions, especially in the cluster II, were likely influenced by the Kuroshio current.

The differences in fish habitat preferences observed in this study highlight the unique characteristics of the Luzon Strait and surrounding waters. In this study, the majority of fish species were found to prefer bathypelagic (36.64%) and reef-associated habitats (35.88%), while demersal-related species were less dominant, accounting for only 7.63% of all detected species. These findings contrast with a previous survey conducted in 230 sites in the South China Sea, where demersal-associated species were identified as the most dominant (37.62%), followed by reef-associated species (32.67%). The divergent habitat preferences observed could be attributed to several factors, including differences in oceanographic conditions, hydrological features, and local environmental factors specific to the Luzon Strait and surrounding waters (Xu et al., 2021; Chen et al., 2022). Overall, the fish habitat preferences observed between this study emphasize the unique ecological characteristics of the Luzon Strait and surrounding waters, as well as the significance of comprehending local factors that shape fish community compositions.

The fish species assessed in this study have significant human uses, particularly in fisheries. We identified a number of fishes with high commercial value, such as Thunnus albacares, Lutjanus madras, Lutjanus erythropterus, and Decapterus maruadsi. However, it is important to note that several of these commercially valuable fish species are facing conservation concerns. Overfishing, along with other anthropogenic impacts such as habitat degradation, biological invasion, and climate change, pose significant threats to fish populations (Chang et al., 2017; Hou et al., 2018). The decreasing population trends reported in the two pelagic-oceanic species, Swordfish (Xiphias gladius) and Sailfish (Istiophorus platypterus), are concerning and highlight the effects of overfishing on their populations. X. gladius, a highly important food and gamefish species, was found only at one station in this survey, indicating a localized distribution. The International Union for Conservation of Nature (IUCN) lists this species as near threatened (NT) due to a significant decline of at least 22% over the past 20 years (Collette et al., 2022a). This decline is consistent with the global trend of decreasing Swordfish populations, particularly in regions such as the Southwest Pacific (58% decrease), South Atlantic (45% decrease), and Indian Ocean (40% decrease), which are heavily impacted by overfishing (Collette et al., 2022b). Similarly, another pelagic-oceanic species, Sailfish (Istiophorus platypterus), was found at three stations in this study (LS1, LS10, and LS12). This species is classified as Vulnerable on the IUCN Red List. Overfishing is also the primary driver of population decline, which are commonly targeted for their value as food and gamefish (Collette et al., 2022a). Globally, Sailfish populations have declined by approximately 30% in the past 12 years (Collette et al., 2022b). The Indian Ocean stock, representing the majority (87%) of the global population, has experienced an estimated 27% to 33% decline due to overfishing. Therefore, it is crucial to implement effective conservation measures, including DNA barcoding, to ensure the survival and recovery of endangered and vulnerable fish species.

Conclusion

Our study demonstrates the efficacy of DNA barcoding as a highly effective tool for accurate species identification of larval fish in the Luzon Strait and adjacent waters. Among the 385 larval fish individuals, we successfully barcoded and identified 354 individuals, revealing 22 orders, 44 families, 93 genera, and 147 species. This comprehensive identification has provided valuable insights into the species composition, diversity, and geographical distribution of larval fish communities in the region. Furthermore, our research has shed light on the significance of the Luzon Strait and surrounding waters as crucial spawning and nursery grounds for various commercially important fish species, such as tunas, groupers, and mackerels. The presence of these species highlights the importance of protecting and managing this ecosystem to safeguard both conservation and sustainable fishing practices. While we acknowledge the challenges of traditional DNA barcoding techniques when dealing with a large sample size of fish larvae, recent advancements in biotechnology, such as next-generation sequencing and environmental DNA (eDNA) metabarcoding, have significantly improved the efficiency of assessing fish diversity. These advancements have facilitated the generation of larger genomic datasets, enabling more accurate and comprehensive species identification in larval fish. Additionally, the species information generated through our study, including habitat types, human uses, and IUCN conservation status, has the potential to support fisheries management efforts, conservation initiatives, and the implementation of effective strategies to protect and preserve the diverse marine ecosystems in the region. In conclusion, this study contributes to our understanding of larval fish diversity and their recruitment dynamics in the Luzon Strait and surrounding waters. By employing DNA barcoding techniques, we have provided valuable insights that can inform fisheries management efforts and conservation actions, ultimately ensuring the sustainable utilization and preservation of this important marine ecosystem.

Data availability statement

The original contributions presented in the study are publicly available. This data can be found here: https://figshare.com/articles/dataset/Assessing_Larval_Fish_Diversity_and_Conservation_Needs_in_the_Luzon_Strait_using_DNA_Barcoding/23949270.

Ethics statement

The animal study was approved by the Animal Experiment Administration Committee of South China Sea Ecological Center, Ministry of Natural Resources of the People’s Republic of China. The study was conducted in accordance with the local legislation and institutional requirements.

Author contributions

HL: Conceptualization, Data curation, Funding acquisition, Investigation, Methodology, Writing – original draft, Writing – review & editing. YC: Conceptualization, Data curation, Methodology, Software, Supervision, Writing – original draft, Writing – review & editing. XL: Methodology, Software, Visualization, Writing – review & editing. PZ: Data curation, Investigation, Writing – review & editing. XX: Data curation, Methodology, Writing – review & editing.

Funding

The author(s) declare financial support was received for the research, authorship, and/or publication of this article. This work was supported by the National Natural Science Foundation of China (41606191) and the National Key Research and Development Program of China (2022YFC3106302).

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmars.2023.1268399/full#supplementary-material

References

Almany G. R., Planes S., Thorrold S. R., Berumen M. L., Bode M., Saenz-agudelo P., et al. (2017). Larval fish dispersal in a coral-reef seascape. Nat. Ecol. Evol. 1, 1–7. doi: 10.1038/s41559-017-0148

PubMed Abstract | CrossRef Full Text | Google Scholar

Azmir I. A., Esa Y., Amin S. M. N., Md Yasin I. S., MD Yusof F. Z. (2017). Identification of larval fish in mangrove areas of Peninsular Malaysia using morphology and DNA barcoding methods. J. Appl. Ichthyol 33, 998–1006. doi: 10.1111/jai.13425

CrossRef Full Text | Google Scholar

Baumgartner G., Nakatani K., Gomes L. C., Bialetzki A., Sanches P. V., Makrakis M. C. (2004). Identification of spawning sites and natural nurseries of fishes in the upper Paraná River, Brazil. Environ. Biol. Fishes 71, 115–125. doi: 10.1007/s10641-004-0098-z

CrossRef Full Text | Google Scholar

Chang C. H., Shao K. T., Lin H. Y., Chiu Y. C., Lee M. Y., Liu S. H., et al. (2017). DNA barcodes of the native ray-finned fishes in Taiwan. Mol. Ecol. Resour. 17, 796–805. doi: 10.1111/1755-0998.12601

PubMed Abstract | CrossRef Full Text | Google Scholar

Chen K. S., Chen H.-S., Chen C. Y., Su Y. L., Meng P. J., Chen M. H. (2022). Multivariate analysis of the spatial species diversity of demersal fish assemblages in relation to habitat characteristics in a subtropical national park, Taiwan. Mar. Biodivers. 52, 1–18. doi: 10.1007/s12526-021-01225-5

CrossRef Full Text | Google Scholar

Chu C., Loh K. H., Ng C. C., Ooi A. L., Konishi Y., Huang S. P., et al. (2019). Using DNA barcodes to aid the identification of larval fishes in tropical estuarine waters (Malacca straits, Malaysia). Zool. Stud. 58, 1–15. doi: 10.6620/ZS.2019.58-30

CrossRef Full Text | Google Scholar

Collette B. B., Di Natale A., Fox W., Graves J., Juan Jorda M., Pohlot B., et al. (2022a). Istiophorus platypterus. IUCN Red List Threatened Species. doi: 10.2305/IUCN.UK.2022-1.RLTS.T170338A46649664.en

CrossRef Full Text | Google Scholar

Collette B. B., Di Natale A., Fox W., Graves J., Juan Jorda M., Pohlot B., et al. (2022b). Xiphias gladius. IUCN Red List Threatened Species. doi: 10.2305/IUCN.UK.2022-1.RLTS.T23148A46625751.en

CrossRef Full Text | Google Scholar

Gu H., Zhao F., Yu W., Li R., Zhao L., Li D. (2019). Distributions of 224, 223Ra and 137Cs in the Luzon Strait and its adjacent waters and their response to Typhoon Rainbow. Acta Oceanol Sin. 38, 1–7. doi: 10.1007/s13131-019-1508-1

CrossRef Full Text | Google Scholar

Hall T. A. (1999). BioEdit: a user-friendly biological sequence alignment editor and analysis program for Windows 95/98/NT. Nucleic Acids Symp. Ser. 41, 95–98. doi: 10.14601/Phytopathol_Mediterr-14998u1.29

CrossRef Full Text | Google Scholar

Hebert P. D. N., Cywinska A., Ball S. L., Jeremy R., B P. R. S. L. (2003). Biological identifications through DNA barcodes. Proc. R. Soc Lond. Ser. B. 270, 313–321. doi: 10.1098/rspb.2002.2218

CrossRef Full Text | Google Scholar

Hou G., Chen W. T., Lu H. S., Cheng F., Xie S. G. (2018). Developing a DNA barcode library for perciform fishes in the South China Sea: Species identification, accuracy and cryptic diversity. Mol. Ecol. Resour. 18, 137–146. doi: 10.1111/1755-0998.12718

PubMed Abstract | CrossRef Full Text | Google Scholar

Hou G., Chen Y., Wang J., Pan C., Lin J., Feng B., et al. (2022). Molecular identification of species diversity using pelagic fish eggs in spring and late autumn-winter in the eastern beibu gulf, China. Front. Mar. Sci. 8, 806208. doi: 10.3389/fmars.2021.806208

CrossRef Full Text | Google Scholar

Hou G., Chen Y., Wang S., Wang J., Chen W., Zhang H. (2021). Formalin-fixed fish larvae could be effectively identified by DNA barcodes: A case study on thousands of specimens in south China sea. Front. Mar. Sci. 8. doi: 10.3389/fmars.2021.634575

CrossRef Full Text | Google Scholar

Huang D., Chen J., Xu L., Wang X., Ning J., Li Y., et al. (2022). Larval fish assemblages and distribution patterns in the zhongsha atoll (Macclesfield bank, south China sea). Front. Mar. Sci. 8. doi: 10.3389/fmars.2021.787765

CrossRef Full Text | Google Scholar

Huang D., Licuanan W. Y., Hoeksema B. W., Chen C. A., Ang P. O., Huang H., et al. (2015). Extraordinary diversity of reef corals in the South China Sea. Mar. Biodivers. 45, 157–168. doi: 10.1007/s12526-014-0236-1

CrossRef Full Text | Google Scholar

Hubert N., Delrieu-Trottin E., Irisson J. O., Meyer C., Planes S. (2010). Identifying coral reef fish larvae through DNA barcoding: A test case with the families Acanthuridae and Holocentridae. Mol. Phylogenet. Evol. 55, 1195–1203. doi: 10.1016/j.ympev.2010.02.023

PubMed Abstract | CrossRef Full Text | Google Scholar

Hubert N., Espiau B., Meyer C., Planes S. (2015). Identifying the ichthyoplankton of a coral reef using DNA barcodes. Mol. Ecol. Resour. 15, 57–67. doi: 10.1111/1755-0998.12293

PubMed Abstract | CrossRef Full Text | Google Scholar

Kim S., Chin B. S., Wang S. Y. (2023). Evidence of the Intrusion of the Oceanic Lightfish (Vinciguerria nimbaria) into Korean Waters Based on High-Throughput Sequencing of Mixed Fish Eggs. J. Mar. Sci. Eng 11, 257. doi: 10.3390/jmse11020257

CrossRef Full Text | Google Scholar

Kindt R., Coe R. (2005). Tree diversity analysis: a manual and software for common statistical methods for ecological and biodiversity studies. Nairobi: World Agroforestry Centre.

PubMed Abstract | Google Scholar

Kong X., Jiang Y., Gong Y., Chen Z., Zhang J., Fang J. (2016). A preliminary study on fishery biology of Ceratoscopelus warmingii in the central and northern South China Sea. South China Fisheries Sci. 12, 117–124. doi: 10.3969/j.issn.2095-0780.2016.04.015

CrossRef Full Text | Google Scholar

Kumar S., Stecher G., Tamura K. (2016). MEGA7: molecular evolutionary genetics analysis version 7.0 for bigger datasets brief communication. Mol. Biol. Evol. 33, 1870–1874. doi: 10.1093/molbev/msw054

PubMed Abstract | CrossRef Full Text | Google Scholar

Kwun H. J. (2018). Species identification of juvenile fishes of the genus Pseudoblennius using mitochondrial DNA barcoding. Mitochondrial DNA Part B Resour. 3, 405–408. doi: 10.1080/23802359.2018.1456982

CrossRef Full Text | Google Scholar

Letunic I., Bork P., Gmbh B. S. (2021). Interactive Tree Of Life (iTOL) v5: an online tool for phylogenetic tree display and annotation. Nucleic Acids Res. 49, 1–4. doi: 10.1093/nar/gkab301

PubMed Abstract | CrossRef Full Text | Google Scholar

Neira F. J., Perry R. A., Burridge C. P., Lyle J. M., Keane J. P. (2015). Molecular discrimination of shelf-spawned eggs of two co-occurring Trachurus spp. (Carangidae) in southeastern Australia: a key step to future egg-based biomass estimates Francisco. ICES J. Mar. Sci. 72, 614–624. doi: 10.1093/icesjms/fsu151

CrossRef Full Text | Google Scholar

Olden J. D., Hogan Z. S., Zanden Vander M. J. (2007). Small size-biased extinction risk of the world ‘ s freshwater and marine fishes. Glob. Ecol. Biogeogr. 16, 694–701. doi: 10.1111/j.1466-8238.2007.00337.x

CrossRef Full Text | Google Scholar

Panprommin D., Iamchuen N., Soontornprasit K., Tuncharoen S. (2020). The utility of DNA barcoding for the species identification of larval fish in the lower ing river, Thailand. Turkish J. Fish Aquat Sci. 20, 671–679. doi: 10.4194/1303-2712-v20_9_02

CrossRef Full Text | Google Scholar

Panprommin D., Soontornprasit K., Tuncharoen S., Iamchuen N. (2021). Efficacy of DNA barcoding for the identification of larval fish species in the Upper and Middle Ing River, Thailand. Gene Rep. 23, 101057. doi: 10.1016/j.genrep.2021.101057

CrossRef Full Text | Google Scholar

Pappalardo A. M., Cuttitta A., Sardella A., Musco M., Maggio T., Patti B., et al. (2015). DNA barcoding and COI sequence variation in Mediterranean lanternfishes larvae. Hydrobiologia 749, 155–167. doi: 10.1007/s10750-014-2161-5

CrossRef Full Text | Google Scholar

Poulsen J. Y., Byrkjedal I., Willassen E., Rees D., Takeshima H., Satoh T. P. (2013). Mitogenomic sequences and evidence from unique gene rearrangements corroborate evolutionary relationships of myctophiformes (Neoteleostei). BMC Evol. Biol. 13, 1–22. doi: 10.1186/1471-2148-13-111

PubMed Abstract | CrossRef Full Text | Google Scholar

Shehata H. R., Naaum A. M., Gardu R. A. (2018). DNA barcoding as a regulatory tool for seafood authentication in Canada. Food Control 92, 147–153. doi: 10.1016/j.foodcont.2018.04.045

CrossRef Full Text | Google Scholar

Sheraliev B., Peng Z. (2021). Molecular diversity of Uzbekistan’s fishes assessed with DNA barcoding. Sci. Rep. 11, 1–12. doi: 10.1038/s41598-021-96487-1

PubMed Abstract | CrossRef Full Text | Google Scholar

Valdez-Moreno M., Vásquez-Yeomans L., Elías-Gutírrez M., Ivanova N. V., Hebert P. D. N. (2010). Using DNA barcodes to connect adults and early life stages of marine fishes from the Yucatan Peninsula, Mexico: Potential in fisheries management. Mar. Freshw. Res. 61, 665–671. doi: 10.1071/MF09222

CrossRef Full Text | Google Scholar

Ward R. D., Zemlak T. S., Innes B. H., Last P. R., Hebert P. D. N. (2005). DNA barcoding Australia's fish species. Phil. Trans. R. Soc. B. 360, 1847–1857. doi: 10.1098/rstb.2005.1716

CrossRef Full Text | Google Scholar

Xing B., Wang C., Wu Q., Wang Y., Chen X., Sun R., et al. (2022). Species identification of larval fish in hawaiian waters using DNA barcoding. Front. Mar. Sci. 9. doi: 10.3389/fmars.2022.825395

PubMed Abstract | CrossRef Full Text | Google Scholar

Xiong W., Huang X., Chen Y., Fu R., Du X., Chen X., et al. (2020). Zooplankton biodiversity monitoring in polluted freshwater ecosystems: a technical review. Environ. Sci. Ecotechnol 1, 100008. doi: 10.1016/j.ese.2019.100008

CrossRef Full Text | Google Scholar

Xiong W., Ni P., Chen Y., Gao Y., Li S., Zhan A. (2019). Biological consequences of environmental pollution in running water ecosystems: A case study in zooplankton. Environ. pollut. 252, 1483–1490. doi: 10.1016/j.envpol.2019.06.055

PubMed Abstract | CrossRef Full Text | Google Scholar

Xiong W., Ni P., Chen Y., Gao Y., Shan B., Zhan A. (2017). Zooplankton community structure along a pollution gradient at fine geographical scales in river ecosystems: The importance of species sorting over dispersal. Mol. Ecol. 26, 4351–4360. doi: 10.1111/mec.14199

PubMed Abstract | CrossRef Full Text | Google Scholar

Xu L., Wang X., Van Damme K., Huang D., Li Y., Wang L., et al. (2021). Assessment of fish diversity in the South China Sea using DNA taxonomy. Fish Res. 233, 105771. doi: 10.1016/j.fishres.2020.105771

CrossRef Full Text | Google Scholar

Yang Y., Gao Y., Chen Y., Li S., Zhan A. (2019a). Interactome-based abiotic and biotic impacts on biodiversity of plankton communities in disturbed wetlands. Divers. Distrib. 25, 1416–1428. doi: 10.1111/ddi.12949

CrossRef Full Text | Google Scholar

Yang Y., Li S., Gao Y., Chen Y., Zhan A. (2019b). Environment-driven geographical distribution of bacterial communities and identification of indicator taxa in Songhua River. Ecol. Indic 101, 62–70. doi: 10.1016/j.ecolind.2018.12.047

CrossRef Full Text | Google Scholar

Yao M., Zhang S., Lu Q., Chen X., Zhang S.-Y., Kong Y., et al. (2022). Fishing for fish environmental DNA: Ecological applications, methodological considerations, surveying designs, and ways forward. Mol. Ecol. 31, 5132–5164. doi: 10.1111/mec.16659

PubMed Abstract | CrossRef Full Text | Google Scholar

Yi S., Zhong J., Huang S., Wang S., Wang W. (2017). Morphological comparison and DNA barcoding of four closely related species in the genera Misgurnus and Paramisgurnus (Cypriniformes: Cobitidae). Biochem. Syst. Ecol. 70, 50–59. doi: 10.1016/j.bse.2016.10.019

CrossRef Full Text | Google Scholar

Zhan A., Hu J., Hu X., Lu W., Wang W., Peng W., et al. (2008). Fast identification of scallop adductor muscles using species-specific microsatellite markers. Eur. Food Res. Technol. 227, 353–359. doi: 10.1007/s00217-007-0728-3

CrossRef Full Text | Google Scholar

Zhan A., Hu J., Hu X., Zhou Z., Hui M., Wang S., et al. (2009). Fine-scale population genetic structure of Zhikong scallop (Chlamys farreri): Do local marine currents drive geographical differentiation? Mar. Biotechnol. 11, 223–235. doi: 10.1007/s10126-008-9138-1

CrossRef Full Text | Google Scholar

Zhan A., MacIsaac H. J. (2015). Rare biosphere exploration using high-throughput sequencing: research progress and perspectives. Conserv. Genet. 16, 513–522. doi: 10.1007/s10592-014-0678-9

CrossRef Full Text | Google Scholar

Zhang J., Hanner R. (2012). Molecular approach to the identification of fish in the South China Sea. PloS One 7, 1–9. doi: 10.1371/journal.pone.0030621

CrossRef Full Text | Google Scholar

Zhang J., Wang Y., Chen Z., Jiang Y., Xu S. (2021). Age and growth of Ceratoscopelus warmingii (Myctophidae) in the South China Sea based on sagittal otolith microstructure. Mar. Biol. Res. 17, 733–743. doi: 10.1080/17451000.2021.2015390

CrossRef Full Text | Google Scholar

Zhang S., Zheng Y., Zhan A., Dong C., Zhao J., Yao M. (2022). Environmental DNA captures native and non-native fish community variations across the lentic and lotic systems of a megacity. Sci. Adv. 8 (6), eabk0097. doi: 10.1175/JPO-D-21-0284.1

PubMed Abstract | CrossRef Full Text | Google Scholar