Next Article in Journal
Influence of Long-Term Wind Variability on the Storm Activity in the Caspian Sea
Next Article in Special Issue
Sustainability Perspective of Minjiang Estuary Coastal Fisheries Management—Estimation of Fish Richness
Previous Article in Journal
Climate Change and Hydrological Response in the Ranwu Lake Basin of Southeastern Tibet Plateau
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Water
Volume 15
Issue 11
10.3390/w15112123
water-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Fish Diversity Monitoring Using Environmental DNA Techniques in the Clarion–Clipperton Zone of the Pacific Ocean

by
Weiyi He
1,
Lei Wang
1,
Danyun Ou
1,
Weiwen Li
1,2,3,4,*,
Hao Huang
1,*,
Rimei Ou
1,
Jinli Qiu
1,
Lei Cai
1,
Lianghua Lin
1 and
Yanxu Zhang
1
1
Laboratory of Marine Biology and Ecology, Third Institute of Oceanography State Oceanic Administration, Xiamen 361005, China
2
National Distant-Water Fisheries Engineering Research Center, Shanghai 201306, China
3
Key Laboratory of Sustainable Exploitation of Oceanic Fisheries Resources, Ministry of Education, Shanghai Ocean University, Shanghai 201306, China
4
Key Laboratory of Oceanic Fisheries Exploration, Ministry of Agriculture and Rural Affairs, Shanghai 201306, China
*
Authors to whom correspondence should be addressed.
Water 2023, 15(11), 2123; https://0-doi-org.brum.beds.ac.uk/10.3390/w15112123
Submission received: 9 May 2023 / Revised: 29 May 2023 / Accepted: 31 May 2023 / Published: 2 June 2023
(This article belongs to the Special Issue Coastal Ecology and Fisheries Management)
Browse Figures
Versions Notes

Abstract

:
Marine fish populations have suffered the consequences of overfishing for a long time, leading to a loss in biodiversity. Traditional methods have been historically used to survey fish diversity but are limited to commercial species, particularly on the high seas. Environmental DNA (eDNA) has been successfully used to monitor biodiversity in aquatic environments. In this study, we monitored fish diversity in the Clarion–Clipperton Zone (CCZ) of the Eastern Pacific Ocean using eDNA metabarcoding. Our results identified 2 classes, 35 orders, 64 families, and 87 genera. The genera Mugil, Scomberomorus, and Scomber had high relative abundance in the mesopelagic and demersal zone. Fish diversity varied with sampling sites, and the greatest number of species was found at a depth of 2500 m. Environmental changes drove fish aggregation, and our results indicated that Chla was negatively correlated with fish communities, while DO was positively correlated with fish communities. This study released the fish diversity pattern and the effects of the environment in the CCZ, which would provide useful information for biodiversity management and an environmental baseline for the International Seabed Authority.

1. Introduction

The Clarion–Clipperton Zone (CCZ) is an area of 6 million square kilometers [1] and is located in the northeastern subtropical mid-Pacific Ocean between Mexico and Hawaii [2]. The CCZ is delimited by two fault zones, the Clarion and Clipperton, and encompasses an extensive range of habitats, including hills, seamounts, fault zones, and vast abyssal plains [3,4]. The CCZ contains many metal nodules rich in manganese, nickel, copper, cobalt, iron, and other rare earth elements, and is an important area for deep-sea manganese nodule mining [5,6].
Mining equipment generates noise pollution and impacts the environment by disturbing ecosystems both physically and chemically [7,8,9]. Deep-sea mining pumps sediment and metallic nodules to the surface, releasing sediment plumes back into the water column. The nutrients of sediment plumes influence pelagic food webs [9]. Mining not only affects the area where metallic nodules are removed, but also disrupts adjacent areas through the redeposition of sediment plumes, affecting wider areas of the seafloor than those directly affected by nodule removal. These changes are likely to persist for decades to centuries [3,10]. Mining of nodules requires appropriate monitoring and conservation strategies [11]. Fish diversity is a part of the environmental baseline.
However, there are few studies on fish in the CCZ. Marine fish diversity is crucial to a healthy ecosystem [12,13]. Understanding the link between fish diversity and ecosystem processes can aid in management and conservation decisions [14]. Fish are indicators for monitoring changing habitats due to their relatively long lifespan, mobility (predation across different trophic levels), and sensitivity to environmental disturbances [14,15]. The fish community has been affected by overfishing and environmental pollution for a long time [16,17,18]. So, it is very important to conduct rapid and continuous surveys of marine fish communities in changing environments. Traditional methods (i.e., nets) for fish monitoring may miss some information [15], such as some species being difficult to capture [19]. In addition, it is hard to observe cryptobenthic or elusive species [20], and it is difficult to conduct scientific surveys by trawling or gillnet in nearshore areas [21]. Net-based monitoring methods are costly, slow, and require survey equipment.
Novel methods for monitoring aquatic organisms have employed environmental DNA (eDNA) to mitigate issues related to conventional netting. eDNA-based methods detect the genetic material of target species in the aquatic environment, and have rapidly emerged as an effective tool to improve aquatic biodiversity monitoring [22,23]. eDNA refers to genetic material extracted from environmental samples (such as water, soil, or sediment) and is derived from mitochondrial or nuclear DNA [24,25]. Sources for eDNA include secretions, feces, urine, tissue, mucus, eggs, and sperm [26,27,28,29]. High-throughput parallel DNA sequencing, known as eDNA metabarcoding, has been increasingly used in eDNA research [30]. eDNA metabarcoding is an important tool for detecting and cataloging biodiversity in local communities and is widely used in marine habitats [31,32]. eDNA metabarcoding may avoid the shortcomings of traditional survey techniques, thereby providing a useful and repeatable method for assessing biodiversity [33]. This method of sampling is cost-effective, using high-throughput screening to survey the abundance and biomass of target species in a non-invasive manner [34,35,36]. In addition, the eDNA method is a valuable tool for detecting many species that are difficult to study by traditional methods (e.g., electrofishing, visual observations, and gillnets) [37]. Environmental factors have a significant effect on driving fish communities in aquatic systems [38]. The marine fish ecosystem is a typical complex system, which is highly dynamic and sensitive to external environmental factors such as sea temperature. Ocean warming is altering fish ecosystems, with profound effects on their ecology (behavior, biomass, range, abundance, etc.), and reducing the diversity and abundance of fish [39]. Chlorophyll concentration is related to primary production and trophic level production, which directly affects the food supply of fish [40]. Understanding the fish community and its spatial distribution is fundamental to assessing the impacts of mining.
eDNA metabarcoding has been used to study fish in different habitats, such as lakes, rivers, estuaries, and the high sea [29,41,42,43]. However, fish diversity studies based on eDNA metabarcoding in the CCZ are rarely reported. In our study, we used eDNA metabarcoding and 12S rRNA primers (MiFish-U and MiFish-E) to reveal mesopelagic and demersal fish diversity and community in the CCZ from two cruises that took place in 2017 and 2018. We monitored fish diversity through eDNA metabarcoding in seamounts and sea basins to provide a bioinformatic basis for enhancing biodiversity protection efforts in the CCZ, and to study biological information at the species and community levels.

2. Materials and Methods

2.1. Field Sampling

The Eastern Pacific Ocean polymetallic nodule region has been affected by the crossover effects of climate change, coastal upwelling, and El Niño-Southern Oscillation (ENSO) from high latitudes. In addition, mining may have caused a huge impact on the environment by disrupting ecosystems, causing physical disturbances, and changing chemical conditions in CCZ. We collected seawater samples from 10 sites in the CCZ of the Eastern Pacific Ocean. The sites were chosen because they contained habitats of seamounts and sea basins. A total of 22 samples were collected on board R/V XIANG YANG HONG 03 during the China Ocean 45 cruise in July 2017 and the China Ocean 45 cruise in August 2018, through WTS-LV Large Volume Water Transfer System (McLANE, Carrollton, TX, USA) for greater biodiversity coverage. The characteristics of the sampling sites are shown in Table 1 and Figure 1. Real-time data of environmental factors in different depths were measured by CTD on board, including temperature (T), salinity (S), turbidity (NTU), Chlorophyll a (Chla), and dissolved oxygen (DO).

2.2. eDNA Collection, Filtration, and Extraction

Seawater samples were collected from mesopelagic and demersal depths using WTS-LV Large Volume Water Transfer System and filtered immediately after collection. All samples and filtration equipment for seawater collection were washed with Milli-Q water before use. Samples were filtered through a glass-fiber membrane with a nominal pore size of 0.3 μm (GF-75, ADVANTEC, Tokyo, Japan). After filtration, filter membranes were placed in cell culture dishes (NEST, Wuxi, China). All samples were immediately frozen at −80 °C until eDNA extraction. eDNA was shredded and extracted using DNeasy PowerWater kit (Qiagen, Hilden, Germany) following the manufacturer’s protocol in the laboratory. eDNA samples were stored at −80 °C until further analysis.

2.3. Metabarcoding of eDNA Samples

eDNA metabarcoding using universal MiFish primer pairs has been shown to amplify short fragments of fish DNA in various taxa from environmental samples [44]. Our samples were analyzed using two universal primer pairs (MiFish-U, MiFish-E, Shengong, Shanghai, China) to amplify the V5 region of the mitochondrial 12S rRNA gene. The multiplex polymerase chain reaction (PCR) volume was 50 µL, including 20 µL of sterile distilled H2O, 25 µL of Taq 2× Master Mix (Vazyme, Nanjing, China), 1 µL of each primer (MiFish-U-F: 5′-GTCGGTAAAACTCGTGCCAGC-3′; MiFish-U-R: 3′-GTTTGACCCTAATCTATGGGGTGATAC-5′; MiFish-E-F: 5′-GTTGGTAAATCTCGTGCCAGC-3′; and MiFish-E-R: 3′-GTTTGATCCTAATC TATGGGGTGATAC-5′), and 1 µL of DNA solution. The thermal cycle PCR process included an initial 2 min denaturation at 94 °C, followed by 30 cycles of denaturation at 98 °C for 5 s each. It was then annealed at 50 °C for 10 s, extended at 72 °C for 10 s, and completed with a final extension at 72 °C for 5 min. Once the PCR was complete, equal amounts of 1× loading buffer (containing SYBR green) and PCR products were mixed and electrophoresed on 1% agarose gels. Samples with a bright main strip of 297 ± 25 bp were selected. The mixed PCR products were then purified with GeneJET Gel Extraction Kit (Thermo Scientific, Waltham, MA, USA).

2.4. Bioinformatics

The quality screening was performed on paired-end reads in the FASTQ format. To analyze the original double C-terminal sequencing data, the sliding window method was used. A window size of 10 bp was used. The analysis results showed that data began to move at 1 bp from the 5′ end of the first base position. A quality score of 20 (Q20) was required for 99% accuracy using FASTQ. The first value was lower than average quality as a result of a truncated sequence. The truncation ceased at 150 bp. Ambiguous bases (Ns) were not permissible.
Following the quality screening, Fast Length Adjustment of Short Reads (FLASH v1.2.7; http://0-ccb-jhu-edu.brum.beds.ac.uk/software/FLASH/ (19 October 2020)) [45] software was used to merge paired-end reads. FLASH software is able to extend short reads by overlapping paired-end reads with a base length of 10 or higher bp and with base mismatch numbers that had less than 10% overlapping base length.
Finally, using index information (i.e., barcode sequence, a short base sequence used to identify the sample), the indexed sequence was matched to the correct corresponding sample.

2.5. Statistical Analysis

Analysis of the sequence was performed using QIIME2 [46], according to the official tutorial (https://docs.qiime2.org/2019.4/tutorials/ (10 January 2021)). The raw data obtained via high-throughput sequencing were screened according to sequence quality, and high-quality sequences were used for subsequent analysis. The raw sequences that passed the quality screening were divided according to index and barcode information, and the barcode sequences were removed. Sequences were then quality filtered, denoised, merged, and chimera was removed using the DADA2 [47]. Deduplicated sequences generated by DADA2 quality control were considered ASVs (amplicon sequence variants) [47,48]. ASV is equivalent to OTU with 100% similarity clustering [49]. Statistics were performed on the length distribution of ASVs to check whether the lengths of these sequences were equivalent to the target fragments or sequences of abnormal lengths. Databases downloaded from NCBI (https://0-www-ncbi-nlm-nih-gov.brum.beds.ac.uk/ (28 February 2021)) and MitoFish (http://mitofish.aori.u-tokyo.ac.jp (28 February 2021)) were used for taxonomy.
Heatmap is plotted using heatmap tools in the Genescloud platform (https://www.genescloud.cn (20 May 2021)). The tool was developed from the heatmap package (V1.0.8), which was slightly modified to improve the layout style. The data were normalized by z-scores. The package uses popular clustering distances and methods implemented in dist and hclust functions in R. The list of distances includes correlation, Euclidean (default), maximum, Manhattan, Canberra, binary, and Minkowski. The clustering method in our analysis is average (UPGMA). Krona software (Brian Ondov edited this page on 5 May 2022, 25 revisions) (https://github.com/marbl/Krona/wiki (19 May 2022)) was used to display community taxonomic composition and its interaction [50]. The Krona figure represents seven taxonomic levels of domain, phylum, class, order, family, genus, and species from inside to outside. The size of the sector reflects the relative abundance of different taxa, and gives specific values. At each taxonomic level, taxa are distinguished by different colors. To compare the differences in species composition between samples and show the species abundance, a heatmap was used for species composition analysis. ASV-level alpha diversity indices, such as the Chao1 richness estimator [51], Observed species, Shannon diversity index [52], and Simpson index [53], were calculated using the ASV table in QIIME2. For the grouped samples, R script can be used to draw the data into boxplots to visually show the differences in alpha diversity among different groups. Kruskal–Wallis rank sum test and Dunn’s test can be used as post hoc tests. The significance of the difference was verified (Kruskal–Wallis test was equivalent to Wilcoxon test for two groups of samples). A principal coordinates analysis (PCoA) was performed to visualize the similarity among the fish communities in different samples using Bray–Curtis index. Redundancy analysis (RDA) was used to analyse the relationship between the fish community and environmental factors [54]. Temperature, salinity, turbidity, Chlorophyll a, and dissolved oxygen were analyzed the correlation to fish assemblage.

3. Results

3.1. eDNA Metabarcoding Sequencing Results

The eDNA metabarcoding assay yielded a total of 2,406,141 sequencing reads. After the quality control process, a total of 1,512,485 reads were retained, corresponding to an average of 68,749 reads per sample. After taxonomic annotation, a total of 2 classes, 35 orders, 64 families, and 87 genera were classified (Table 2). It was determined that all sequences from the water samples belonged to the classes Chondrichthyes and Actinopteri. The Chondrichthyes class contained three families: Carcharhinidae, Hexanchidae, and Myliobatidae. Within the Carcharhinidae family, the species Prionace glauca (blue shark) and Scoliodon laticaudus (Spadenose shark) were found. Species of the Hexanchidae family and the Myliobatidae family were unclassified. The Myliobatidae family was listed under the Convention on International Trade in Endangered Species of Wild Fauna and Flora (CITES).

3.2. Species Composition and Diversity

The relative abundance of species found at different sampling locations showed distinct differences. Statistical analysis of the non-singleton data showed that the relative abundance of genera Mugil, Scomberomorus, and Scomber was high in all samples over two years (Figure 2). Mugil accounted for 14.61% of the overall relative abundance at the genus taxonomic level, while Scomberomorus accounted for 11.66%.
Taxonomic composition analysis indicated that the genus Serrivomer showed high relative abundance at 1000 and 1200 m in 2017. The relative abundance of Scomberomorus was the highest at 3000 m, and the relative abundance of Mugil was the highest at 1000 m. The results of the species composition analysis showed that the five genera with the highest relative abundance in the whole water column of each site in 2017 were Mugil, Scomberomorus, Serrivomer, Konosirus, and Scomber. The relative abundance of the genus Scomberomorus was the highest in the whole water column of each station in 2018, and the relative abundance of Mugil, Rhynchopelates, Coryphaena, and Scomber decreased in turn. The water column was dominated by Scombriformes, followed by Lophiiformes (Figure 3). Krona diagram results showed that Scomberomorus niphonius and Mugil cephalus were the dominant fish species. Scomberomorus niphonius and Mugil cephalus are highly commercial species, according to FishBase. Lasiognathus sp. NMMBP 9030 belonged to bathypelagic fish and was the dominant species of Lophiiformes. Heatmap results (Figure 4) showed that the community of fish was different among various depths. The richness of fish at 1000 m was significantly greater than that at other sampling depths. The main distribution range of the species Scoliodon laticaudus and Prionace glauca were below 2500 m depth. Larimichthys crocea was mainly found in water depths of 700 m. The heatmap also indicated that fish richness in the seamounts was higher than that in the sea basin.

3.3. Community Diversity

The results of the alpha diversity parameters (Chao1, observed species, Shannon index, and Simpson) for each sampling depth in 2018 tended to be consistent with that of 2017, with non-significant differences (Figure 5). The variation among depths was neither significant for the Chao1 index nor for the observed species index in 2017. Shannon index and Simpson index gradually increased with sampling depths (Table 3). Alpha diversity results showed that the Chao1 index, observed species index, Shannon index, and Simpson index decreased with water depth in 2018 (Table 3). Alpha diversity parameters at 1000 m were significantly larger than those at 3000 m (p = 0.05), including for the Chao1 index and observed species index.

3.4. Species Distribution by Depth

Fish species distribution at different depths was detected using eDNA analysis. Our results indicated that the greatest number of species was detected at a depth of 2500 m, while the fewest number of species were detected at 5000 m (Figure 6). The variation at 2500 m was relatively large over three random repeated sampling efforts. Samples collected at depths of 700 m, 1200 m, and 5000 m were not subject to repeated sampling efforts, which may have contributed to accidental results. The beta diversity showed that the spatial structure based on different depths was not obvious when the fish community was ordinated by Bray–Curtis PCoA (Figure 7). The relationship between fish community and environmental factors was clarified by RDA results, and the proportion of fish community variation was explained by axe 1 (14.6%) and axe 2 (7.8%). As shown in Figure 8, Chla and DO were the main influencing factors of fish community structure in the CCZ. Chla was negatively correlated with the fish community, but DO was positively correlated with the fish community.

4. Discussion

Our study demonstrates that eDNA can be an effective method for studying fish diversity. eDNA can be collected from any type of aquatic or wild environment for monitoring fish ecology, composition, and distribution, as well as for monitoring endangered and invasive species [55,56]. eDNA metabarcoding is an efficient and versatile method that does not require extensive taxonomic expertise [31]. Compared to traditional sampling methods, eDNA methods are non-invasive and not destructive to the environment [44]. In CCZ areas, nodules support distinct species and community structures, such as sessile organisms, numerous other megafaunal, and meiofaunal and microbial taxa [4,57]. The demersal fauna has a limited supply of exotic food and is characterized by slow growth, replenishment, reproduction, and recovery after disturbance. Mining would affect benthic communities, which in turn affects fish distribution through the food chain. Removal of polymetallic nodules due to mining would lead to a loss of food-web integrity and a substantial decline in faunal biodiversity [58]. So, the determination of fish diversity by eDNA provides a valuable community assessment before mining.
Since eDNA released by different individuals within a population coexists in the aquatic environment, eDNA analysis can be extended to the assessment of diversity within populations [15]. eDNA has been widely used to detect the presence of plants and animals, and fish have become a common study subject in recent studies. The DADA2 bioinformatics pipeline uses a denoising algorithm to obtain ASVs to infer the true biological sequence, discriminating differences in sequence variants as small as one nucleotide [15,47,59]. The ASV is considered to be equivalent to the DNA sequence present in the original environmental sample and has been proposed to improve the accuracy of assessing the intraspecific diversity of fish populations [47]. In a similar study, 66 functional entities were detected using eDNA technology on Malpelo Island, a remote marine protected area, and the functional richness for eDNA was higher than that in underwater videos [60]. eDNA metabarcoding detects more fish than underwater visual census techniques [61]. Research has shown that eDNA methods are capable of gathering a spectrum of functional traits, showing the most functionally diverse and least redundant fish assemblages [62].
At the genus level, Mugil accounted for 16.03% of the total relative abundance, and Scomberomorus accounted for 10.69% of the total relative abundance during our 2017 sampling efforts. The relative abundance of Mugil was the highest at 700 m, while the relative abundance of Scomberomorus was the highest at 2500 m. In 2018, Scomberomorus accounted for 13.92% of the total relative abundance, and Mugil accounted for 11.29% of the total relative abundance. The relative abundance of Mugil was highest at 1000 m, and the relative abundance of Scomberomorus was highest at 3000 m. The following species had high relative abundance and were widespread among the study area: Mugil cephalus, Scomberomorus niphonius, Konosirus punctatus, Scomber japonicus, and Serrivomer sector.
According to data retrieved from FishBase (https://fishbase.se/search.php (24 March 2021)), we identified that four of the fish species we detected with high relative abundance were migratory fish. Due to the high mobility and widespread distribution of migratory fish, we speculated that the release of eDNA during migratory processes results in a higher detection rate than other fish. We found that the species Scoliodon laticaudus, Prionace glauca, and Harpadon nehereus are ‘Near Threatened’ fish on the IUCN Red List, according to FishBase. Other species found on the IUCN Red List included Larimichthys crocea (‘Critically Endangered’) and Epinephelus fuscoguttatus (‘Vulnerable’). Interestingly, our results showed that Chondrichthyes fish are detected in the bathyal zone. Chondrichthyan fish, such as Scoliodon laticaudus and Prionace glauca, are important consumers in most marine ecosystems that are commonly found to depths of 1000 m but are uncommon, exceedingly rare, or quite possibly absent deeper than 3000 m [63]. A survey by Priede et al. [64] illustrated that the deepest Chondrichthyes below 3000 m was a shark, Centrophorus squamosus, captured at 3280 m by baited long line. The sharks (Centroscymnus coelolepis) were reported to be deepest at 3700 m [65]. Our results showed for the first time that eDNA metabarcoding detected sharks, Scoliodon laticaudus and Prionace glauca, at depths over 1000 m in the CCZ. We investigated the diversity of fish taxa and found critically endangered species Larimichthys crocea at DY45-II-CC-S06 and DY45-III-CCW-S01, compared with other sites. However, Larimichthys crocea is a commercially important species in China and distributed in the Western Pacific regions. We hypothesize that the eDNA of Larimichthys crocea flows into the sea with domestic water on research vessels and is collected by WTS-LV Large Volume Water Transfer System.
In terms of global biodiversity in oceanic areas, Molinos et al. [66] analyzed the distribution of biodiversity under different climate change models. Their results showed that with an increase in temperature, the total number of species decreased at low latitudes, increased at middle latitudes, and remained unchanged at high latitudes. Costello and Chaudhary [20] analyzed the vertical distribution of biodiversity in a changing environment, and found that biodiversity decreased with increasing water depth (distribution law of indexing). Burrows et al. [67] further analyzed the horizontal and vertical migration rules of biodiversity under climate change, with results showing that in areas with small temperature changes, biological migration was not obvious, while in areas with large temperature changes, organisms mainly adapted to temperature changes through vertical migration. These findings were similar to our results; high Shannon diversity index results were found at a depth of 1000 m, but the largest number of fish species were detected at 2500 m. Low temperature makes eDNA gradually degrades compared with surface temperature. So, we hypothesize that the DNA degradation time, DNA sink, temperature, and light contribute to this phenomenon, and we should pay attention to this process in further study. Additionally, the sink rate in a huge depth is very low; for DNA from 2500 m to 3000 m or even deeper water, the period of degradation would be much longer, preventing the degradation that could make it so that fish are not detected. The variance explained by PCoA in different depths was not statistically significant, which may indicate the connectivity of vertical habitats through dispersal, migration, or movement of seawater that carries eDNA [33]. Low temperature has been a larger contributor to the fish community similarity at different depths. The study of Takahara et al. showed that temperature may be the main driving factor of eDNA distribution [68]. Higher temperatures directly increase DNA degradation through the denaturation of DNA molecules, and indirectly degrade eDNA by increasing enzyme kinetics and microbial metabolism [69]. The decay rates of fish eDNA in marine water appear to be between 6.9 and 71.1 h [70]. Low temperature can preserve eDNA, but eDNA gradually degrades with the increase in sedimentation time. The Spatial dynamic of the fish community was affected by environmental factors. The RDA result of this study indicated that the main environmental factors influencing fish distribution were Chla and DO. Our results are similar to the study of Diao et al., in which Chla was negatively correlated with fish assemblages and affected fish assemblages by cascade effects [29]. Studies showed that the interaction of temperature and DO drives fish to use horizontal and vertical space [71].
The eDNA metabarcoding has been widely used in research aimed at fish diversity and detecting a large number of fish species. Our results demonstrate the usefulness of eDNA metabarcoding in conservation and management purposes for marine fishes. We found the DNA signature of Near Threatened fish, Critically Endangered fish, and Vulnerable fish in the CCZ. eDNA metabarcoding in biodiversity assessments will be crucial as humans continue to balance the use and conservation of marine resources in marine ecosystems.

Author Contributions

Conceptualization, W.H. and W.L.; Methodology, L.W. and W.L.; Software, D.O., L.C. and Y.Z.; Validation, W.H. and D.O.; Formal analysis, W.H. and W.L.; Investigation, L.W.; Resources, D.O., H.H., R.O., J.Q. and L.L.; Data curation, L.W.; Writing—original draft, W.H.; Writing—review & editing, W.L. and H.H. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the Open Funding Project of the Key Laboratory of Sustainable Exploitation of Oceanic Fisheries Resources, the Ministry of Education (No. A1-2006-22-200206), the Deep Sea Habitats Discovery Project of China Deep Ocean Affairs Administration (No. DY-XZ-04), the Monitoring and Protection of Ecology and Environment in the East Pacific Ocean (No. DY135-E2-5-06), and the Monitoring and Protection of Ecology and Environment of Seamounts in the Western Pacific Ocean (No. DY135-E2-2-05).

Data Availability Statement

Data not available due to nature of this research restrictions.

Acknowledgments

The author would like to thank all the scientific researchers who provided assistance during the cruise.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Cairns, S.D. New abyssal Primnoidae (Anthozoa: Octocorallia) from the Clarion-Clipperton Fracture Zone, equatorial northeastern Pacific. Mar. Biodivers. 2016, 46, 141–150. [Google Scholar] [CrossRef]
  2. Jones, D.O.B.; Simon-Lledo, E.; Amon, D.J.; Bett, B.J.; Caulle, C.; Clement, L.; Connelly, D.P.; Dahlgren, T.G.; Durden, J.M.; Drazen, J.C.; et al. Environment, ecology, and potential effectiveness of an area protected from deep-sea mining (Clarion Clipperton Zone, abyssal Pacific). Prog. Oceanogr. 2021, 197, 102653. [Google Scholar] [CrossRef]
  3. Simon-Lledo, E.; Bett, B.J.; Huvenne, V.A.I.; Schoening, T.; Benoist, N.M.A.; Jeffreys, R.M.; Durden, J.M.; Jones, D.O.B. Megafaunal variation in the abyssal landscape of the Clarion Clipperton Zone. Prog. Oceanogr. 2019, 170, 119–133. [Google Scholar] [CrossRef]
  4. Smith, C.R.; Clark, M.R.; Goetze, E.; Glover, A.G.; Howell, K.L. Editorial: Biodiversity, Connectivity and Ecosystem Function Across the Clarion-Clipperton Zone: A Regional Synthesis for an Area Targeted for Nodule Mining. Front. Mar. Sci. 2021, 8, 797516. [Google Scholar] [CrossRef]
  5. Jakiel, A.; Palero, F.; Błażewicz, M. Deep ocean seascape and Pseudotanaidae (Crustacea: Tanaidacea) diversity at the Clarion-Clipperton Fracture Zone. Sci. Rep. 2019, 9, 17305. [Google Scholar] [CrossRef] [Green Version]
  6. Zinssmeister, C.; Wilke, T.; Hoppenrath, M. Species diversity of dinophysoid dinoflagellates in the Clarion-Clipperton Fracture Zone, eastern Pacific. Mar. Biodivers. 2017, 47, 271–287. [Google Scholar] [CrossRef]
  7. Taboada, S.; Riesgo, A.; Wiklund, H.; Paterson, G.L.J.; Koutsouveli, V.; Santodomingo, N.; Dale, A.C.; Smith, C.R.; Jones, D.O.B.; Dahlgren, T.G.; et al. Implications of population connectivity studies for the design of marine protected areas in the deep sea: An example of a demosponge from the Clarion-Clipperton Zone. Mol. Ecol. 2018, 27, 4657–4679. [Google Scholar] [CrossRef]
  8. Wang, C.-S.; Liao, L.; Xu, H.-X.; Xu, X.-W.; Wu, M.; Zhu, L.-Z. Bacterial diversity in the sediment from polymetallic nodule fields of the Clarion-Clipperton Fracture Zone. J. Microbiol. 2010, 48, 573–585. [Google Scholar] [CrossRef]
  9. Lindh, M.V.; Maillot, B.M.; Smith, C.R.; Church, M.J. Habitat filtering of bacterioplankton communities above polymetallic nodule fields and sediments in the Clarion-Clipperton zone of the Pacific Ocean. Environ. Microbiol. Rep. 2018, 10, 113–122. [Google Scholar] [CrossRef] [PubMed]
  10. Lambshead, P.J.; Brown, C.J.; Ferrero, T.J.; Hawkins, L.E.; Smith, C.R.; Mitchell, N.J. Biodiversity of nematode assemblages from the region of the Clarion-Clipperton Fracture Zone, an area of commercial mining interest. BMC Ecol. 2003, 3, 1. [Google Scholar] [CrossRef]
  11. Li, Q.; Lei, Y.; Liu, J.; Shen, Y.; Huang, H.; Wang, C.; Li, H.; Li, T. Characteristics of foraminiferal communities in the western Clarion–Clipperton Zone revealed by eDNA metabarcoding. J. Sea Res. 2022, 189, 102286. [Google Scholar] [CrossRef]
  12. Mora, C.; Tittensor, D.P.; Myers, R.A. The completeness of taxonomic inventories for describing the global diversity and distribution of marine fishes. Proc. R. Soc. B Biol. Sci. 2008, 275, 149–155. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Costello, M.J.; Wilson, S.; Houlding, B. Predicting total global species richness using rates of species description and estimates of taxonomic effort. Syst. Biol. 2012, 61, 871–883. [Google Scholar] [CrossRef] [Green Version]
  14. Stuart-Smith, R.D.; Bates, A.E.; Lefcheck, J.S.; Duffy, J.E.; Baker, S.C.; Thomson, R.J.; Stuart-Smith, J.F.; Hill, N.A.; Kininmonth, S.J.; Airoldi, L.; et al. Integrating abundance and functional traits reveals new global hotspots of fish diversity. Nature 2013, 501, 539–542. [Google Scholar] [CrossRef] [Green Version]
  15. Wang, S.P.; Yan, Z.G.; Hanfling, B.; Zheng, X.; Wang, P.Y.; Fan, J.T.; Li, J.L. Methodology of fish eDNA and its applications in ecology and environment. Sci. Total Environ. 2021, 755, 142622. [Google Scholar] [CrossRef]
  16. Myers, R.A.; Worm, B. Rapid worldwide depletion of predatory fish communities. Nature 2003, 423, 280–283. [Google Scholar] [CrossRef] [PubMed]
  17. Frank, K.T.; Petrie, B.; Choi, J.S.; Leggett, W.C. Trophic cascades in a formerly cod-dominated ecosystem. Science 2005, 308, 1621–1623. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  18. Genner, M.J.; Sims, D.W.; Wearmouth, V.J.; Southall, E.J.; Southward, A.J.; Henderson, P.A.; Hawkins, S.J. Regional climatic warming drives long-term community changes of British marine fish. Proc. R. Soc. B-Biol. Sci. 2004, 271, 655–661. [Google Scholar] [CrossRef]
  19. Thomsen, P.F.; Willerslev, E. Environmental DNA—An emerging tool in conservation for monitoring past and present biodiversity. Biol. Conserv. 2015, 183, 4–18. [Google Scholar] [CrossRef]
  20. Costello, M.J.; Chaudhary, C. Marine Biodiversity, Biogeography, Deep-Sea Gradients, and Conservation. Curr. Biol. 2017, 27, R511–R527. [Google Scholar] [CrossRef]
  21. Zou, K.; Chen, J.; Ruan, H.; Li, Z.; Guo, W.; Li, M.; Liu, L. eDNA metabarcoding as a promising conservation tool for monitoring fish diversity in a coastal wetland of the Pearl River Estuary compared to bottom trawling. Sci. Total Environ. 2020, 702, 134704. [Google Scholar] [CrossRef]
  22. Pont, D.; Rocle, M.; Valentini, A.; Civade, R.; Jean, P.; Maire, A.; Roset, N.; Schabuss, M.; Zornig, H.; Dejean, T. Environmental DNA reveals quantitative patterns of fish biodiversity in large rivers despite its downstream transportation. Sci. Rep. 2018, 8, 10361. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Burian, A.; Mauvisseau, Q.; Bulling, M.; Domisch, S.; Qian, S.; Sweet, M. Improving the reliability of eDNA data interpretation. Mol. Ecol. Resour. 2021, 21, 1422–1433. [Google Scholar] [CrossRef] [PubMed]
  24. Lacoursière-Roussel, A.; Rosabal, M.; Bernatchez, L. Estimating fish abundance and biomass from eDNA concentrations: Variability among capture methods and environmental conditions. Mol. Ecol. Resour. 2016, 16, 1401–1414. [Google Scholar] [CrossRef]
  25. Hervé, A.; Domaizon, I.; Baudoin, J.M.; Dejean, T.; Gibert, P.; Jean, P.; Peroux, T.; Raymond, J.C.; Valentini, A.; Vautier, M.; et al. Spatio-temporal variability of eDNA signal and its implication for fish monitoring in lakes. PLoS ONE 2022, 17, e0272660. [Google Scholar] [CrossRef]
  26. Takeuchi, A.; Iijima, T.; Kakuzen, W.; Watanabe, S.; Yamada, Y.; Okamura, A.; Horie, N.; Mikawa, N.; Miller, M.J.; Kojima, T.; et al. Release of eDNA by different life history stages and during spawning activities of laboratory-reared Japanese eels for interpretation of oceanic survey data. Sci. Rep. 2019, 9, 6074. [Google Scholar] [CrossRef] [Green Version]
  27. Rey, A.; Carney, K.J.; Quinones, L.E.; Pagenkopp Lohan, K.M.; Ruiz, G.M.; Basurko, O.C.; Rodriguez-Ezpeleta, N. Environmental DNA Metabarcoding: A Promising Tool for Ballast Water Monitoring. Environ. Sci. Technol. 2019, 53, 11849–11859. [Google Scholar] [CrossRef]
  28. Miya, M.; Sato, Y.; Fukunaga, T.; Sado, T.; Poulsen, J.Y.; Sato, K.; Minamoto, T.; Yamamoto, S.; Yamanaka, H.; Araki, H.; et al. MiFish, a set of universal PCR primers for metabarcoding environmental DNA from fishes: Detection of more than 230 subtropical marine species. R. Soc. Open Sci. 2015, 2, 150088. [Google Scholar] [CrossRef] [Green Version]
  29. Diao, C.; Jia, H.; Guo, S.; Hou, G.; Xian, W.; Zhang, H. Biodiversity exploration in autumn using environmental DNA in the South China sea. Environ. Res. 2022, 204, 112357. [Google Scholar] [CrossRef]
  30. Fujii, K.; Doi, H.; Matsuoka, S.; Nagano, M.; Sato, H.; Yamanaka, H. Environmental DNA metabarcoding for fish community analysis in backwater lakes: A comparison of capture methods. PLoS ONE 2019, 14, e0210357. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  31. Monuki, K.; Barber, P.H.; Gold, Z. eDNA captures depth partitioning in a kelp forest ecosystem. PLoS ONE 2021, 16, e0253104. [Google Scholar] [CrossRef]
  32. Keck, F.; Blackman, R.C.; Bossart, R.; Brantschen, J.; Couton, M.; Hürlemann, S.; Kirschner, D.; Locher, N.; Zhang, H.; Altermatt, F. Meta-analysis shows both congruence and complementarity of DNA and eDNA metabarcoding to traditional methods for biological community assessment. Mol. Ecol. 2022, 31, 1820–1835. [Google Scholar] [CrossRef]
  33. Nguyen, B.N.; Shen, E.W.; Seemann, J.; Correa, A.M.S.; O’Donnell, J.L.; Altieri, A.H.; Knowlton, N.; Crandall, K.A.; Egan, S.P.; McMillan, W.O.; et al. Environmental DNA survey captures patterns of fish and invertebrate diversity across a tropical seascape. Sci. Rep. 2020, 10, 6729. [Google Scholar] [CrossRef] [Green Version]
  34. Maruyama, A.; Sugatani, K.; Watanabe, K.; Yamanaka, H.; Imamura, A. Environmental DNA analysis as a non-invasive quantitative tool for reproductive migration of a threatened endemic fish in rivers. Ecol. Evol. 2018, 8, 11964–11974. [Google Scholar] [CrossRef] [Green Version]
  35. Sato, Y.; Miya, M.; Fukunaga, T.; Sado, T.; Iwasaki, W. MitoFish and MiFish Pipeline: A Mitochondrial Genome Database of Fish with an Analysis Pipeline for Environmental DNA Metabarcoding. Mol. Biol. Evol. 2018, 35, 1553–1555. [Google Scholar] [CrossRef] [Green Version]
  36. Bessey, C.; Jarman, S.N.; Simpson, T.; Miller, H.; Stewart, T.; Keesing, J.K.; Berry, O. Passive eDNA collection enhances aquatic biodiversity analysis. Commun. Biol. 2021, 4, 236. [Google Scholar] [CrossRef] [PubMed]
  37. Deutschmann, B.; Mueller, A.-K.; Hollert, H.; Brinkmann, M. Assessing the fate of brown trout (Salmo trutta) environmental DNA in a natural stream using a sensitive and specific dual-labelled probe. Sci. Total Environ. 2019, 655, 321–327. [Google Scholar] [CrossRef]
  38. Mondal, R.; Bhat, A. Temporal and environmental drivers of fish-community structure in tropical streams from two contrasting regions in India. PLoS ONE 2020, 15, e0227354. [Google Scholar] [CrossRef] [Green Version]
  39. Li, J.; Convertino, M. Temperature increase drives critical slowing down of fish ecosystems. PLoS ONE 2021, 16, e0246222. [Google Scholar] [CrossRef] [PubMed]
  40. Hu, C.; Harrison, D.P.; Hinton, M.G.; Siegrist, Z.C.; Kiefer, D.A.J.F.O. Habitat analysis of the commercial tuna of the Eastern Tropical Pacific Ocean. Fish. Oceanogr. 2018, 27, 417–434. [Google Scholar] [CrossRef]
  41. DiBattista, J.D.; Fowler, A.M.; Riley, I.J.; Reader, S.; Hay, A.; Parkinson, K.; Hobbs, J.A. The use of environmental DNA to monitor impacted coastal estuaries. Mar. Pollut. Bull. 2022, 181, 113860. [Google Scholar] [CrossRef] [PubMed]
  42. Kumar, G.; Reaume, A.M.; Farrell, E.; Gaither, M.R. Comparing eDNA metabarcoding primers for assessing fish communities in a biodiverse estuary. PLoS ONE 2022, 17, e0266720. [Google Scholar] [CrossRef]
  43. Shen, M.; Xiao, N.; Zhao, Z.; Guo, N.; Luo, Z.; Sun, G.; Li, J. eDNA metabarcoding as a promising conservation tool to monitor fish diversity in Beijing water systems compared with ground cages. Sci. Rep. 2022, 12, 11113. [Google Scholar] [CrossRef]
  44. Miya, M.; Gotoh, R.O.; Sado, T. MiFish metabarcoding: A high-throughput approach for simultaneous detection of multiple fish species from environmental DNA and other samples. Fish. Sci. 2020, 86, 939–970. [Google Scholar] [CrossRef]
  45. Magoč, T.; Salzberg, S.L. FLASH: Fast length adjustment of short reads to improve genome assemblies. Bioinformatics 2011, 27, 2957–2963. [Google Scholar] [CrossRef] [Green Version]
  46. Bolyen, E.; Rideout, J.R.; Dillon, M.R.; Bokulich, N.A.; Abnet, C.C.; Al-Ghalith, G.A.; Alexander, H.; Alm, E.J..; Arumugam, M.; Asnicar, F.; et al. Reproducible, interactive, scalable and extensible microbiome data science using QIIME 2. Nat. Biotechnol. 2019, 37, 852–857. [Google Scholar] [CrossRef]
  47. Callahan, B.J.; McMurdie, P.J.; Holmes, S.P. Exact sequence variants should replace operational taxonomic units in marker-gene data analysis. ISME J. 2017, 11, 2639–2643. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  48. Katoh, K.; Misawa, K.; Kuma, K.; Miyata, T. MAFFT: A novel method for rapid multiple sequence alignment based on fast Fourier transform. Nucleic Acids Res. 2002, 30, 3059–3066. [Google Scholar] [CrossRef] [Green Version]
  49. McDonald, D.; Price, M.N.; Goodrich, J.; Nawrocki, E.P.; DeSantis, T.Z.; Probst, A.; Andersen, G.L.; Knight, R.; Hugenholtz, P. An improved Greengenes taxonomy with explicit ranks for ecological and evolutionary analyses of bacteria and archaea. ISME J. 2012, 6, 610–618. [Google Scholar] [CrossRef]
  50. Ondov, B.D.; Bergman, N.H.; Phillippy, A.M. Interactive metagenomic visualization in a Web browser. BMC Bioinform. 2011, 12, 385. [Google Scholar] [CrossRef] [Green Version]
  51. Chao, A. Nonparametric-Estimation of the Number of Classes in a Population. Scand. J. Stat. 1984, 11, 265–270. [Google Scholar]
  52. Shannon, C.E. A mathematical theory of communication. Bell Syst. Tech. J. 1948, 27, 379–423. [Google Scholar] [CrossRef] [Green Version]
  53. Simpson, E.H. Measurement of Diversity. Nature 1949, 163, 688. [Google Scholar] [CrossRef]
  54. Yan, T.; He, J.; Yang, D.; Ma, Z.; Chen, H.; Zhang, Q.; Deng, F.; Ye, L.; Pu, Y.; Zhang, M.; et al. Fish Community Structure and Biomass Particle-Size Spectrum in the Upper Reaches of the Jinsha River (China). Animals 2022, 12, 3412. [Google Scholar] [CrossRef] [PubMed]
  55. Moushomi, R.; Wilgar, G.; Carvalho, G.; Creer, S.; Seymour, M. Environmental DNA size sorting and degradation experiment indicates the state of Daphnia magna mitochondrial and nuclear eDNA is subcellular. Sci. Rep. 2019, 9, 12500. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  56. Xing, Y.; Gao, W.; Shen, Z.; Zhang, Y.; Bai, J.; Cai, X.; Ouyang, J.; Zhao, Y. A Review of Environmental DNA Field and Laboratory Protocols Applied in Fish Ecology and Environmental Health. Front. Environ. Sci. 2022, 10, 725360. [Google Scholar] [CrossRef]
  57. Simon-Lledó, E.; Bett, B.J.; Huvenne, V.A.I.; Schoening, T.; Benoist, N.M.A.; Jones, D.O.B. Ecology of a polymetallic nodule occurrence gradient: Implications for deep-sea mining. Limnol. Oceanogr. 2019, 64, 1883–1894. [Google Scholar] [CrossRef] [Green Version]
  58. Stratmann, T.; Soetaert, K.; Kersken, D.; van Oevelen, D. Polymetallic nodules are essential for food-web integrity of a prospective deep-seabed mining area in Pacific abyssal plains. Sci. Rep. 2021, 11, 12238. [Google Scholar] [CrossRef]
  59. Pérez-Burillo, J.; Trobajo, R.; Leira, M.; Keck, F.; Rimet, F.; Sigró, J.; Mann, D.G. DNA metabarcoding reveals differences in distribution patterns and ecological preferences among genetic variants within some key freshwater diatom species. Sci. Total Environ. 2021, 798, 149029. [Google Scholar] [CrossRef]
  60. Marques, V.; Castagné, P.; Polanco, A.F.; Borrero-Pérez, G.H.; Hocdé, R.; Guérin, P.; Juhel, J.B.; Velez, L.; Loiseau, N.; Letessier, T.B.; et al. Use of environmental DNA in assessment of fish functional and phylogenetic diversity. Conserv. Biol. 2021, 35, 1944–1956. [Google Scholar] [CrossRef]
  61. Valdivia-Carrillo, T.; Rocha-Olivares, A.; Reyes-Bonilla, H.; Domínguez-Contreras, J.F.; Munguia-Vega, A. Integrating eDNA metabarcoding and simultaneous underwater visual surveys to describe complex fish communities in a marine biodiversity hotspot. Mol. Ecol. Resour. 2021, 21, 1558–1574. [Google Scholar] [CrossRef]
  62. Aglieri, G.; Baillie, C.; Mariani, S.; Cattano, C.; Calò, A.; Turco, G.; Spatafora, D.; Di Franco, A.; Di Lorenzo, M.; Guidetti, P.; et al. Environmental DNA effectively captures functional diversity of coastal fish communities. Mol. Ecol. 2021, 30, 3127–3139. [Google Scholar] [CrossRef] [PubMed]
  63. Treberg, J.R.; Speers-Roesch, B. Does the physiology of chondrichthyan fishes constrain their distribution in the deep sea? J. Exp. Biol. 2016, 219, 615–625. [Google Scholar] [CrossRef] [Green Version]
  64. Priede, I.G.; Froese, R.; Bailey, D.M.; Bergstad, O.A.; Collins, M.A.; Dyb, J.E.; Henriques, C.; Jones, E.G.; King, N. The absence of sharks from abyssal regions of the world’s oceans. Proc. R. Soc. B Biol. Sci. 2006, 273, 1435–1441. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  65. Forster, G.R. Line Fishing on the Continental Slope the Selective Effect of Different Hook Patterns. J. Mar. Biol. Assoc. United Kingd. 1973, 53, 749–751. [Google Scholar] [CrossRef]
  66. Molinos, J.G.; Halpern, B.S.; Schoeman, D.S.; Brown, C.J.; Kiessling, W.; Moore, P.J.; Pandolfi, J.M.; Poloczanska, E.S.; Richardson, A.J.; Burrows, M.T. Climate velocity and the future global redistribution of marine biodiversity. Nat. Clim. Chang. 2016, 28, 5849–5858. [Google Scholar] [CrossRef]
  67. Burrows, M.T.; Bates, A.E.; Costello, M.J.; Edwards, M.; Edgar, G.J.; Fox, C.J.; Halpern, B.S.; Hiddink, J.G.; Pinsky, M.L.; Batt, R.D.; et al. Ocean community warming responses explained by thermal affinities and temperature gradients. Nat. Clim. Chang. 2019, 9, 959. [Google Scholar] [CrossRef]
  68. Takahara, T.; Minamoto, T.; Yamanaka, H.; Doi, H.; Kawabata, Z. Estimation of fish biomass using environmental DNA. PLoS ONE 2012, 7, e35868. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  69. Ravanat, J.L.; Douki, T.; Cadet, J. Direct and indirect effects of UV radiation on DNA and its components. J. Photochem. Photobiol. B Biol. 2001, 63, 88–102. [Google Scholar] [CrossRef]
  70. Collins, R.A.; Wangensteen, O.S.; O’Gorman, E.J.; Mariani, S.; Sims, D.W.; Genner, M.J. Persistence of environmental DNA in marine systems. Commun. Biol. 2018, 1, 185. [Google Scholar] [CrossRef] [Green Version]
  71. Cooke, S.J.; Bergman, J.N.; Twardek, W.M.; Piczak, M.L.; Casselberry, G.A.; Lutek, K.; Dahlmo, L.S.; Birnie-Gauvin, K.; Griffin, L.P.; Brownscombe, J.W.; et al. The movement ecology of fishes. J. Fish Biol. 2022, 101, 756–779. [Google Scholar] [CrossRef] [PubMed]
Water 15 02123 g001 550
Figure 1. Location of sample sites in the Clarion–Clipperton Zone; the red points show the areas of sample collection. (PMN—Polymetallic nodules, APEIs—Areas of Particular Environmental Interest).
Figure 1. Location of sample sites in the Clarion–Clipperton Zone; the red points show the areas of sample collection. (PMN—Polymetallic nodules, APEIs—Areas of Particular Environmental Interest).
Water 15 02123 g001
Water 15 02123 g002 550
Figure 2. Relative abundance of reads classified by genus.
Figure 2. Relative abundance of reads classified by genus.
Water 15 02123 g002
Water 15 02123 g003 550
Figure 3. Krona diagram showing fish community composition.
Figure 3. Krona diagram showing fish community composition.
Water 15 02123 g003
Water 15 02123 g004 550
Figure 4. Taxonomic heatmap showing composition at different depths and sites.
Figure 4. Taxonomic heatmap showing composition at different depths and sites.
Water 15 02123 g004
Water 15 02123 g005 550
Figure 5. Alpha diversity parameters at different depths.
Figure 5. Alpha diversity parameters at different depths.
Water 15 02123 g005
Water 15 02123 g006 550
Figure 6. Species detected at various sampling depths.
Figure 6. Species detected at various sampling depths.
Water 15 02123 g006
Water 15 02123 g007 550
Figure 7. PCoA depicting similarity in community composition at different depths. The ellipses are 95% confidence interval.
Figure 7. PCoA depicting similarity in community composition at different depths. The ellipses are 95% confidence interval.
Water 15 02123 g007
Water 15 02123 g008 550
Figure 8. Redundancy analysis of the relationship between fish distribution and environmental factors.
Figure 8. Redundancy analysis of the relationship between fish distribution and environmental factors.
Water 15 02123 g008
Table
Table 1. Sampling information and samples for analysis.
Table 1. Sampling information and samples for analysis.
Sample IDDateLatitude (N)Longitude (E)Depth (m)
DY50A-KW1-S01_1000 m12 August 201810.004205.6655198
DY50A-KW1-S01_3000 m12 August 201810.004205.6655198
DY45-I-NLG-S06_2500 m6 August 201720.637161.2994941
DY45-I-NA-S05_2500 m26 August 201720.140156.6953584
DY45-II-KW1-S01_1000 m18 September 201710.990205.7415234
DY45-III-KW1-S01_4000 m10 October 201710.990205.7415236
DY45-III-KW1-S01_2500 m10 October 201710.990205.7415236
DY45-II-KW1-S05_3000 m11 September 201710.056205.6585169
DY45-II-KW1-S05_5000 m11 August 201710.056205.6585169
DY45-II-KW1-S05_1000 m11 July 201710.056205.6585169
DY45-III-CCW-S01_700 m16 October 20179.207202.0151304
DY45-III-CCW-S01_1200 m16 October 20179.207202.0151304
DY50A-KW1-S02_1000 m14 August 20189.888206.5395083
DY50A-KW1-S02_3000 m14 August 20189.888206.5395083
DY45-II-S40_2500 m21 September 201710.187205.4075147
DY45-II-S40_4000 m21 September 201710.187205.4075147
DY45-II-S40_1000 m21 September 201710.187205.4075147
DY45-II-CC-S06_1000 m22 August 201712.976206.7295524
DY45-II-CC-S06_4000 m22 August 201712.976206.7295524
DY45-II-CC-S06_2500 m22 August 201712.976206.7295524
DY50B-A8-S03_3000 m6 September 201813.323139.2954987
DY50B-A8-S03_1000 m6 September 201913.323139.2954987
Table
Table 2. Taxonomic composition of fish species.
Table 2. Taxonomic composition of fish species.
OrderFamilyGenus/Species
AlepocephaliformesAlepocephalidaeAlepocephalus sp.
Leptoderma spp.
Platytroctidae
AnguilliformesNemichthyidaeAvocettina sp
SerrivomeridaeSerrivomer beanii
Serrivomer sector
AnguillidaeAnguilla spp.
AulopiformesEvermannellidaeCoccorella sp.
SynodontidaeSaurida spp.
Harpadon nehereus
BeloniformesExocoetidaeCheilopogon sp.
BelonidaeCololabis saira
BeryciformesBerycidaeBeryx spp.
MelamphaidaePoromitra spp.
Scopelogadus sp.
CarangiformesEcheneidaeRemora remora
CoryphaenidaeCoryphaena hippurus
CarangidaeDecapterus macrosoma
CentrarchiformesTerapontidaeRhynchopelates oxyrhynchus
Kyphosidae
ClupeiformesClupeidaeKonosirus punctatus
Konosirus spp.
Sardinella spp.
Sardinops melanostictus
Sardinops sagax
EngraulidaeEncrasicholina punctifer
Engraulis japonicus
ElopiformesElopidaeElops hawaiensis
MegalopidaeMegalops spp.
GadiformesMacrouridae
GadidaeGadus spp.
MelanonidaeMelanonus zugmayeri
GobiiformesGobiidaeAcanthogobius hasta
Exyrias spp.
Luciogobius sp.
Trypauchen vagina
GonorynchiformesChanidaeChanos sp.
GonorynchidaeGonorynchus breviatus
KurtiformesApogonidaeOstorhinchus fasciatus
LampriformesLampridaeLampris guttatus
LophiiformesMelanocetidaeMelanocetus murrayi
Melanocetus johnsonii
ThaumatichthyidaeLasiognathus sp. NMMBP 9030
DiceratiidaeBufoceratias sp.
MugiliformesMugilidaePlaniliza spp.
Mugil cephalus
Mugil spp.
MyctophiformesMyctophidaeCeratoscopelus spp.
Diaphus aliciae
Diaphus spp.
Myctophum spp.
Nannobrachium sp.
Symbolophorus spp.
Lampanyctus spp.
Hygophum proximum
Protomyctophum thompsoni
Lampanyctus tenuiformis
Bolinichthys pyrsobolus
Bolinichthys spp.
OphidiiformesOphidiidaeHoplobrotula armata
PempheriformesHowellidaeHowella sp.
LateolabracidaeLateolabrax japonicus
Lateolabrax maculatus
PerciformesScorpaenidaeScorpaena spp.
SebastidaeHelicolenus spp.
SerranidaeEpinephelus fuscoguttatus
Epinephelus sp.
PleuronectiformesCynoglossidae
PomacentridaePomacentridaeAbudefduf spp.
SalmoniformesSalmonidaeSalvelinus leucomaenis
Salvelinus sp.
SciaenidaeSciaenidaeJohnius spp.
Larimichthys crocea
Pennahia argentata
Miichthys miiuy
ScombriformesBramidaeEumegistus illustris
ChiasmodontidaeDysalotus sp.
CentrolophidaePsenopsis anomala
GempylidaeGempylus serpens
Lepidocybium flavobrunneum obrunneum flavobrunneum flavobrunneum flavobrunneum
ScombridaeEuthynnus alletteratus
Euthynnus affinis
Euthynnus spp.
Auxis thazard
Scomberomorus niphonius
Scomberomorus spp.
Katsuwonus pelamis
Thunnus alalunga
Thunnus spp.
Scomber japonicus
Scomber sp.
Acanthocybium solandri
Auxis rochei
NomeidaeNomeus gronovii
Cubiceps squamiceps
SiganidaeSiganidae
SiluriformesPlotosidaePlotosus lineatus
AriidaeNetuma spp.
StomiiformesGonostomatidaeCyclothone atraria
Cyclothone pallida
Cyclothone obscura
PhosichthyidaeVinciguerria nimbaria
StomiidaeChauliodus sp.
Stomias sp.
Idiacanthus antrostomus
Thysanactis spp.
Sternoptyx obscura
SternoptychidaeArgyropelecus sladeni
Argyropelecus sp.
SpariformesSparidaeAcanthopagrus schlegelii
SphyraenidaeSphyraenidaeSphyraena spp.
SyngnathiformesSyngnathidae
TetraodontiformesMolidaeMasturus lanceolatus
UranoscopiformesAmmodytidaeAmmodytes sp.
CarcharhiniformesCarcharhinidaeScoliodon laticaudus
Prionace glauca
Hexanchidae
MyliobatiformesMyliobatidae
Table
Table 3. Alpha diversity parameters at different depths.
Table 3. Alpha diversity parameters at different depths.
DepthsChao1Observed_SpeciesShannonSimpson
DY45 < 2000 m24.1333323.933333.085990.768111
DY45 2000 m–3000 m43.24342.066673.241670.815328
DY45 4000 m–5000 m30.0354329.23.6299750.885577
DY50 1000 m36.261135.466673.784880.883085
DY50 3000 m17.6333317.566673.12690.851768
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

He, W.; Wang, L.; Ou, D.; Li, W.; Huang, H.; Ou, R.; Qiu, J.; Cai, L.; Lin, L.; Zhang, Y. Fish Diversity Monitoring Using Environmental DNA Techniques in the Clarion–Clipperton Zone of the Pacific Ocean. Water 2023, 15, 2123. https://0-doi-org.brum.beds.ac.uk/10.3390/w15112123

AMA Style

He W, Wang L, Ou D, Li W, Huang H, Ou R, Qiu J, Cai L, Lin L, Zhang Y. Fish Diversity Monitoring Using Environmental DNA Techniques in the Clarion–Clipperton Zone of the Pacific Ocean. Water. 2023; 15(11):2123. https://0-doi-org.brum.beds.ac.uk/10.3390/w15112123

Chicago/Turabian Style

He, Weiyi, Lei Wang, Danyun Ou, Weiwen Li, Hao Huang, Rimei Ou, Jinli Qiu, Lei Cai, Lianghua Lin, and Yanxu Zhang. 2023. "Fish Diversity Monitoring Using Environmental DNA Techniques in the Clarion–Clipperton Zone of the Pacific Ocean" Water 15, no. 11: 2123. https://0-doi-org.brum.beds.ac.uk/10.3390/w15112123

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop