Next Article in Journal
A SLAM Method Based on Multi-Robot Cooperation for Pipeline Environments Underground
Next Article in Special Issue
The Effect of the Artificial Reef on the Structure and Function of Sediment Bacterial Community
Previous Article in Journal
A Holistic Review of GM/IM Selection Methods from a Structural Performance-Based Perspective
Previous Article in Special Issue
Human Activities Aggravate VOC Pollution in the Huangshui River of the Tibetan Plateau
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

eDNA Metabarcoding Analysis of the Composition and Spatial Patterns of Fish Communities in the Sanbanxi Reservoir, China

1
School of Animal Science, Guizhou University, Guiyang 550025, China
2
Fisheries Research Institute, Academy of Agricultural Sciences of Guizhou, Guiyang 550025, China
*
Author to whom correspondence should be addressed.
Sustainability 2022, 14(20), 12966; https://doi.org/10.3390/su142012966
Submission received: 7 July 2022 / Revised: 4 October 2022 / Accepted: 8 October 2022 / Published: 11 October 2022
(This article belongs to the Special Issue Wetlands: Conservation, Management, Restoration and Policy)

Abstract

:
The construction of a reservoir dam alters the environment within its basin, including composition of the fish community, fish biodiversity, and the river ecosystem itself. This study was conducted in the Sanbanxi Reservoir and used eDNA metabarcoding technology comprising eDNA capture and extraction, PCR amplification, sequencing and database comparison analysis, and other environmental DNA metabarcoding standardized analysis processes to characterize the composition and diversity of fish communities and assess their current status. A total of 48 species of fish were detected. Previously, 68 species of fish were screened and identified in this reservoir based on the reports of Dai and Gu. The results for fish community composition showed that species of the order Cypriniformes are still the most dominant in the Reservoir with 38 species of cyprinids, accounting for 90.81% of all OTUs. Carp were no longer the dominant species, and Spinibarbus denticulatus, Homalopteridae, Cobitidae, and Sisoridae were not detected, with the exception of Misgurnus anguillicaudatus (Cobitidae). These families have the common characteristic of being adapted to survive in fast-water, sandstone substrate habitats. The results also show that two of the sampling sites, sbx03 and sbx10, significantly differed from other sampling sites due to their geographical environment. The impact of the construction of reservoirs on freshwater fish communities is extreme, since the transformation from a lotic to a lentic habitat contributes to habitat destruction and constrains fish in movement. The change in the aquatic environment before and after the storage of water in the Sanbanxi Reservoir has reduced the number of fish species found in the reservoir, and species characteristically found in fast moving, rapids habitats are virtually absent. The profound change in the aquatic environment from that of a lotic to a lentic habitat leads to changes in the composition of fish populations in the reservoir and to a certain extent a reduction in the ecological stability and species diversity within the reservoir. Therefore, the protection of fish diversity in the reservoir is of great significance to the stability of the ecosystem.

1. Introduction

Environmental DNA (eDNA) refers to the sum of DNA fragments that can be directly extracted from environmental samples (such as water, soil, air, ice cores, and others) [1,2]. The first application of eDNA analysis to water samples was in marine biology research [3]. With the continuous development of sequencing technology, such as high-throughput sequencing of tens of billions of base pairs in a relatively short period, eDNA technology has become a promising approach in monitoring, although optimization is still required [4,5]. The eDNA metabarcoding is the description of biological communities using the sequencing information resulting from the amplification of a specific locus or group of loci. Various molecular means are used in the molecular taxa identification of DNA extracted from environmental media. Compared with traditional monitoring methods, eDNA metabarcoding is non-invasive and is highly sensitive, cheap, and determinations can be made rapidly [6,7,8]. Currently, eDNA metabarcoding technology is being widely used to detect invasive [6,9,10] and endangered [11,12,13] species, as well as to assess terrestrial, aquatic, and marine biodiversity [14,15,16]. Through eDNA, biological information about species composition, such as alpha and beta diversity, can be determined more quickly than through using traditional methods. It can be seen from the research reports in recent years that there is a significant positive correlation between detected eDNA and fish population abundance [17]. In addition, eDNA can more intuitively and quickly reflect fish distribution and activity ranges [18].
Human activities, such as the construction of dams and reservoirs, adversely impact rivers. For example, the number of fish species decreased from 141 to 78 before and after the establishment of the Jinsha River Reservoir in China [19]. The impacts of dams can, therefore, potentially contribute to the decline of aquatic biodiversity globally [20,21,22]. Dams affect fish populations in a number of ways [23], including through changes in water storage and flow rate, habitat destruction, and the blocking of species movement [24]. Flow is a major determinant of the physical habitat in streams, which in turn is a major selective factor guiding biotic composition. Aquatic species have evolved life history strategies primarily in direct response to natural flow regimes [23]. Changes in discharge regimes and reduced sediment supply below dams are known to cause: an elevated base flow; channel incision, constriction, or widening; changes in bed material; and the loss of spawning habitats [25,26]. Stepped dams have resulted in a change in the composition of fish communities and the extinction of some endemic fish in the Yangtze River [24] and stingrays in other locations [25,27,28] due to a change in the river flow regime. Changes in hydrological conditions constitute a significant threat to sustaining many fish stocks [26,29,30]. Fish with floating eggs do not develop well in standing water because an adequate rafting duration is essential for their embryonic development [31]. Thus, the impact of dams on fish communities is relatively severe. As an essential part of aquatic ecosystems, species diversity and the composition of fish communities affect the health of aquatic ecosystems. Conserving fish diversity is, therefore, of great significance for the stability of these environments and the habitats within them.
Sanbanxi Reservoir, whose dam site is located in Jinping County, has an adjustable reservoir pool elevation in which the water level may be raised or lowered through the control of outflow from the dam. The reservoir was formed as a result of the construction of the Sanbanxi power station, which is the second large hydropower station on the mainstem of the Yuanjiang River. The reservoir has a total capacity of approximately 4.09 billion m3 [32]. It commenced impounding water in January 2006, and the water level rose to 475 m. The average water depth at the head of the reservoir is over 130 m, and along with rare flooding, obvious thermal stratification frequently occurs throughout the year. The outflow tends to be colder than the natural inflow in the spring and summer, which has a negative impact on the spawning of downstream fish [33]. The habitat of the reservoir has thus been affected, causing alterations to the structure and distribution of fish stocks in the basin. This study aims to explore the impact of habitat changes in the waters of the reservoir on the composition and geographical distribution of fish populations. The results of our study provide data for environmental assessment, environmental planning, governance, and restoration in this watershed.

2. Materials and Methods

2.1. Sampling Scheme Design

The experimental samples were collected in November 2021 at the Sanbanxi Reservoir in Jinping County, Guizhou Province. The sampling range was within 30 km above the dam site. Originally, each sample point was set to be spaced 3 km apart, but due to the irregularity of the riverbanks and slower water flow, the sample points were fine-tuned to allow collection of more eDNA residue (see Figure 1). At each sampling point, we took 5 L of water from the upper (surface), middle (2.5 m), and lower (5 m) levels of the water column for a total of 15 L of water, and then concentrated this to a final volume of 1 L. The water samples were stored in a refrigerator at 0 °C after collection. The samples were pumped and filtered in the laboratory on the same day (filter membrane with 0.45 μm pore size), followed by storage in a freezer at −80 °C.

2.2. DNA Extraction and PCR Amplification

Total genomic DNA was extracted from the filter membrane using PowerWater DNA Isolation Kits. To avoid contamination, extraction was carried out independently for each sample. The mitochondrial gene 12S rRNA, the most widely used molecular marker in fish diversity monitoring, was used as the target in amplification. The fish-specific primer “Tele02” (Tele02-F: 5′-AAACTCGTGCCAGCCACC-3′, Tele02-R: 5′-GGGTATCTAATCCCAGTTTG-3′), which was developed by Taberlet [34], was used in amplification with addition of the sample-specific barcode sequence. The PCR reactions were prepared in total volumes of 20 μL as follows: 4 μL 5× FastPfu Buffer, 2 μL dNTPs, 0.4 μL FastPfu Polymerase, 0.8 μL of each 10 μM primer, 2–5 μL of genomic DNA (10 ng/μL). Finally, ddH2O was added to a total system volume of 20 μL. The PCR conditions included an initial denaturation step at 95 °C for 5 min, followed by 27 cycles at 95 °C for 30 s, 55 °C for 30 s, and 72 °C for 45 s, with a final duration of 10 min at 72 °C. The products of PCR were sent to the Biozeron Company for sequencing.

2.3. Data Analysis

2.3.1. DNA Analysis and Species Identification

Samples were sequenced and 250 bp paired-end readings were obtained using a MiSeq sequencer (Illumina, Inc., San Diego, CA, USA). Sequence data were sorted for each sample using the indices. The quality of the data from the readings was controlled and the results filtered using Trimmomatic software (parameters—leading: 0; trailing: 20; sliding window: 10/20; and minlen: 75). Finally, the high-quality sequences of each sample were retained by splitting the barcode and primer sequences using cutadapt with default parameters, and the sequence direction was corrected according to the positive and negative barcode and primer directions during the process. Paired reads were merged. OTU (operational taxonomic unit) clustering analysis was performed according to a sequence similarity threshold of ≥97%, and representative OTU sequences were compared with sequences in the MitoFish “http://mitofish.aori.u-tokyo.ac.jp (accessed on 10 March 2022)” and NCBI “https://www.ncbi.nlm.nih.gov (accessed on 4 March 2022)” databases for classification annotation and generation of the corresponding OTU abundance table. Sequences of unmatched species were culled from the samples at each sampling point and the sequences of matched species were averaged.

2.3.2. Statistical Analysis

After removing the sequence data showing high identity to non-fish organisms (such as bacteria, birds, amphibians, mammals, and so forth.), the remaining filtered data were compared with fish sequences and showed ≥97% identity and E-values ≤ 10−5, and OTUs corresponding to the same species were merged. If any OTU did not show adequate similarity for identification at the species level, statistical analyses were then carried out for a higher taxonomic level such as genus and family. The sequence number proportion of each species in each sample was calculated in Excel. The taxonomic information of fish was improved by referring to the Fishbase database “https://www.fishbase.in/home.htm (accessed on 18 April 2022)”. Finally, a bar chart of fish composition at each sampling point was drawn using R software (Version 3.1.3; Packets: ggplot2, picante, scatterplot3d, ellipse, maptools, vegan and ape).

Alpha Diversity Analysis

For this study, the Chao1 index, ACE index, Shannon index, Simpson index, and Coverage index were selected to reflect community richness, community diversity, and community coverage.
Chao1: is an index used to estimate the number of OTU contained in a sample using the chao1 algorithm. The calculation formula used in this analysis is as follows:
S chao1   =   S obs + n 1 n 1 1 2 n 2 + 1
ACE: is an index that estimates species diversity using rare species and used to estimate the number of OTUs contained in a community.
ACE = S abund   + S rare C ace + F 1 C ace + γ 2 ace
Shannon: One of the biodiversity indices used to estimate the environmental sample. It similarly to the Simpson Index.
H shannon   = i = 1 S obs n i N ln n i N
Simpson Index: One of the biodiversity indices used to estimate environmental samples, proposed by Edward Hugh Simpson (1949). In ecology it is commonly used to quantitatively describe the biodiversity of an area.
D simpson   = i = 1 S obs n i ( n i 1 ) N ( N 1 )
Coverage: refers to the coverage of various data libraries, and the higher the value, the higher the probability that the sequence in the sample will be measured, and the lower the probability that it will not be measured. The index reflects whether the sequencing results represent the true condition of the organisms in the environmental sample.
C   =   1 n 1 N
The alpha diversity index analysis utilized in this experiment was the mothur software (version v.1.30.1 http://www.mothur.org/wiki/Schloss_SOP#Alpha_diversity (accessed on 30 April 2022)) and the OTU similarity level for exponential evaluation was 97% (0.97) (Sobs: the number of OTU actually observed; n1 = number of OTUs with only one sequence; n2 = number of OTUs with only two sequences; ni= the number of sequences contained in the OTU; Sabund = number of species enriched (abundance threshold greater than n); Srare = number of rare (abundance threshold less than or equal to n) species; F1 = number of species containing only 1 individual; γ2ace = estimate of the coefficient of variation for rare species).

Beta Diversity Analysis

In this study, to identify potential principal components that affect differences in community composition in the sample by reducing dimensions. PCoA analysis of community compositions at different sample points was based on the Bray–Curtis distance algorithm to explore differences or similarities in community composition among different groups of samples (PCoA uses the ellipse package and the maptools package).

3. Results

3.1. Species Composition

The eDNA metabarcoding analysis resulted in a total of 648 OTUs from 10 sampling sites. The Venn diagram of the OTU distribution at each sampling point is shown in Figure 2. Sampling point sbx05 had 537 had the largest number of OTUs. The lowest number, 446 OTUs, occurred at sbx07. Sbx03 and sbx10 had 7 OTUs. A total of 648 OTUs were designated for 48 species of fish belonging to 38 genera and 12 families in 4 orders (Table 1). Among them, species of the order Cypriniformes were dominant, with 32 species belonging to 26 genera and 3 families, accounting for 93.05% of the OTU sequence counts. Siluriformes comprised 3.50% of OTUs and included 9 species belonging to 6 genera and 3 families. Perciformes comprised 2.78% of OTUs and included 6 species belonging to 5 genera and 5 families. Cyprinodontiformes accounted for 0.67% OTUs with only 1 species. Economically important fish such as Onychostoma lini, Onychostoma rara, and Siniperca kneri were present, as well as ornamental fish such as Myxocyprinus asiaticus. In addition, a total of 10 exotic species were found, including Megalobrama terminalis, Megalobrama amblycephala, Culter mongolicus, Squaliobarbus curriculus, Ictalurus punctatus, Cirrhinus mrigala, Myxocyprinus asiaticus, Micropterus salmoides, Oreochromis niloticus, and Procypris rabaudi.

3.2. Spatial Patterns of the Community Structure Composition and Species Diversity Analysis

The composition of the fish communities at sbx03 and sbx10 was significantly different from that of the other eight sample points. The proportions of the fish species also varied greatly between sbx3 and sbx10 (Figure 3), contributing to apparent differences in the composition of community structure. Alpha diversity analysis was performed on all sampling sites, and the data are shown in Table 2. The values ranged from 681 to 824 in the Chao Index, from 685 to 811 in the ACE Index, from 0.06 to 0.24 for the Simpson Index, and from 4.4 to 5.3 for the Shannon Index. Based on the Alpha Diversity Index, there was no difference in fish species at the 10 sampling points. Coverage values of each sample were 0.99.

3.3. Beta Diversity Analysis

The Bray–Curtis distance matrix was used for PCoA analysis. It showed that sampling points 3 and 10 were distant from the remaining points that were clustered (Figure 4). The divergence of these points from those remaining is also reflected in the heat map (Figure S1). The divergence value between sbx03, sbx10 and other sampling points ranges between 0.6 and 0.8 (Figure S1), indicating the fish composition of sampling sites 3 and 10 were significantly different from those of the other sampling sites. The other eight sampling sites had similar fish diversity.

4. Discussion

4.1. Changes in Fish Composition

The fish composition results showed that Cypriniformes are still the dominant fish species in the reservoir. There were 48 species of fish detected in the study, whereas 68 species of fish were earlier documented in the basin [35]. Changes in fish diversity and community structure are inevitable during reservoir filling. Many reports have documented the elimination of fish species and the reduction in fishery productivity following the construction of reservoirs [36,37]. The reservoir waters were inhabited primarily by 38 species of cyprinids, accounting for 90.81% of OTUs. The dominant species at present is Spinibarbus denticulatus. Cobitidae, Homalopteridae, and Sisoridae species were not found in the present study, although Sioridae were considered relatively abundant in this reservoir according to a study conducted several decades ago [35]. The Three Gorges Dam showed a similar phenomenon of certain fishes adapting to the rapids in lotic conditions prior to water storage, and other fish adapting to lentic, still water, after the reservoir was filled for water storage [38]. The results indicate that water storage in the reservoir does clearly influence the structure of its fish community.
Both the Chao1 and ACE Indexes were used to estimate the species richness. The higher the value, the greater the species richness of the sample. The Simpson Index is used to assess all organisms’ diversity index in the sample, and the higher the Shannon value, the higher the community diversity. The Simpson Index is also used to evaluate the microbial diversity index in the sample, while the larger the value of the Simpson Index, the lower the community diversity. Based on the alpha diversity results of this experiment, the large values of the Chao1 and ACE indices indicate the high abundance of unique OTUs in the 10 sampling points. The Simpson and Shannon Indices provide the comparative species-richness and community diversity of each sampling site. The overall ecological environment of the reservoir is relatively good. The chao1 and ACE indices of the sbx03 sample point were lower relative to other sample points, while the Simpson and Shannon Indices reflect the higher species diversity of the sample point, as can be seen in (Figure 3). This may be due to the geographical location of the site and human factors. The sbx03 sampling point is located at a bend in the river which may form an area of food sedimentation and fish accumulation. In addition, there are some human activities around that point that may have a certain impact on its richness and species diversity. Similar results were found in the PCoA analysis. Coverage values of each sample were 0.99, indicating the high coverage of the OTU sequencing data, and that they reflect the actual situation.

4.2. The Associated Effects of eDNA Monitoring

The eDNA detection is a new method that has emerged in studying fish diversity and species monitoring [39], but it lacks technical standardization and generalization. In this study, various factors affected the detection of eDNA metabarcoding. The first was the temperature. Higher temperatures exacerbate eDNA degradation, with temperatures close to 0 °C allowing better preservation of eDNA [40,41]. The samples were placed in a 0 °C car refrigerator after collection, although percolation in open air may have led to partial eDNA degradation. This may have caused the sbx03 sample point data to appear abnormal. Temperature also affects the degradation rate of eDNA in the water sample, so there were differences between different locations, but because the sampling was conducted in the winter, the temperature differences are not large. The second factor was that of transportation and storage time. The eDNA analysis is time-sensitive, as eDNA degrades within a certain period of time, so the preservation and transportation of eDNA water samples affect the detection rate, as concluded in Hinlo’s study [42,43]. We completed processing within 24 h from the termination of sampling to the extraction of filtration and preservation, so preservation of the detection rate of eDNA in our study was accurate, given that the rate of eDNA degradation is relatively low within 24 h at low temperatures. In addition, false negatives and false positives can occur in eDNA detection, possibly due to the use of PCR primers and probes of insufficient specificity [44]. The primer we used was the commonly used “Tele02”. Some false negatives and false positives occurred, and sequences corresponding to some deep-sea marine life that would not be found in this basin were detected and we manually screened them out. We tried to reduce the impact of these factors on this study as much as possible during sample preparation.
The turbidity of the water and the distribution of eDNA substances in the water body also had an impact on this study [45,46]. However, these uncontrollable factors were relatively mild in this experiment, and the water samples we collected were clear, the filtration time was fast, and the sampling was divided to reduce these effects as much as possible.

4.3. Changes in Fish Community Composition and Structure Caused by the Construction of the Reservoir

Dams block about 63% of the world’s rivers and are one of the most important factors contributing to the decline in global freshwater fish diversity [47]. In this study, we verified 48 species of fish, a decrease from the previously report 68 species for this reservoir [35]. The construction of large hydropower dams affects biodiversity by changes at the community, species, and even at genetic levels. In the Jinsha River in China, a total of 141 fish species were present before the establishment of the reservoir, but the number decreased to 78 species after the construction of the dam and filling of the reservoir [19]. The construction of the Three Gorges Dam reservoir impacted the composition and abundance of fish in its waters (Perera et al., 2014) and seriously affected the spawning of Chinese sturgeon, resulting in a sharp decline in the Chinese sturgeon population [48].
The construction of a reservoir widens the original water surface, increasing the reservoir capacity. Impacts of dams may, therefore, contribute to the decline of aquatic biodiversity throughout the reservoir [20,21,22]. Dams affect fish populations in several ways [23], such as through changes in water storage and flow rate, habitat destruction, and by restricting species movement [24]. These factors cause ecological disturbances resulting in changes in the aquatic animals present and in the biodiversity [37]. Based on the results of this study, some fish adapted to a lotic, rapid-flow sand and gravel environment were eliminated, while fish of the genus Spinibarbus and Megalobrama survived. Homalopteridae, Cobitidae, and Sisoridae were not detected, all having the common shared characteristic of preferring fast-water and sandstone environments. It is suggested that the construction of dams or reservoirs has a great impact on rivers by changing the flow state from lotic to lentic conditions, resulting in changes in the aquatic animal community and the disturbance of the overall biodiversity, which is similar to earlier findings [49]. With the completion of the dam, the resulting impoundment for water storage causes the river to change from turbulent flowing water to a static reservoir [24]. In addition, the completion of the dam results in a large temperature difference between the dam site and the rest of the reservoir. Our results showed there were apparent differences in the composition and abundance of fish between sbx03, sbx10 and the other sampling sites. The sbx10 sample site was located at the dam site of the reservoir, with low temperature and dissolved oxygen [50]. The optimal temperature for Spinibarbus is between 19 and 31 °C, whereas the temperature at the dam site was only 6 °C [50,51].
After the Sanbanxi Reservoir was impounded, reservoir water level fluctuation has become the major factor affecting part of the spawning grounds. Changes in hydrological conditions constitute a significant threat to fish stocks [26,29,30]. Fishes with floating eggs cannot develop adequately in standing water because sufficient rafting duration is essential for the embryonic development of drifting eggs [31]. Thus, the impact of dams on fish communities is relatively severe. As an essential part of aquatic ecosystems, species diversity and the composition of fish communities determine the health of these ecosystems. Conserving fish diversity is therefore of great significance for their stability. Based on our findings, we propose the following recommendations: first, strengthen the conservation of species diversity, increase stability, conduct appropriate release of stocked and planted fish into the reservoir, increase the diversity of species in reservoirs, and control fishing and fish harvest methods; secondly, establish spawning grounds around reservoirs to provide a stable place for fish to reproduce; thirdly, conduct regular surveys to tailor conservation measures in real time to the fish diversity of the reservoir.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/su142012966/s1, Figure S1. Heatmap analysis for Distance matrix.

Author Contributions

F.H. and X.M. were involved in the conception and design. X.M. and H.Y. wrote the draft. H.Y. analyzed and interpreted the data. X.Z. (Xue Zhong) and P.Z. ran the experiments. X.Z. (Xianjun Zhou) and X.D. completed the review and editing. W.M., S.Z. completed sample collection. All authors have read and agreed to the published version of the manuscript.

Funding

This study was funded by the Science and Technology Support Project of Guizhou Province (2020-1Y099), Young scientist fund of Guizhou Academy of Agricultural Sciences (2021-16) and the Introduced Talent Research Fund of Guizhou University (2014-11).

Institutional Review Board Statement

The animal study protocol was conducted following the guidelines for the care and use of experimental animals of China (GB/T35892 2018) and approved by the ethics committee of Guizhou University.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript; or in the decision to publish the results.

References

  1. Levy-Booth, D.J.; Campbell, R.G.; Gulden, R.H.; Hart, M.M.; Powell, J.R. Cycling of extracellular DNA in the soil environment. Soil Biol. Biochem. 2007, 39, 2977–2991. [Google Scholar] [CrossRef]
  2. Pietramellara, G.; Ascher, J.; Borgogni, F.; Borgogni, F.; Ceccherini, M.T. Extracellular DNA in soil and sediment: Fate and ecological relevance. Biol. Fertil. Soils 2009, 45, 219–235. [Google Scholar] [CrossRef]
  3. Ficetola, G.F.; Miaud, C.; Pompanon, F.; Taberlet, P. Species detection using environmental DNA from water samples. Biol. Lett. 2008, 4, 423–425. [Google Scholar] [CrossRef] [Green Version]
  4. Shendure, J.; Ji, H. Next-generation DNA sequencing. Nat. Biotechnol. 2008, 26, 1135–1145. [Google Scholar] [CrossRef]
  5. Baird, D.J.; Hajibabaei, M. Biomonitoring 2.0: A new paradigm in ecosystem assessment made possible by next-generation DNA sequencing. Mol. Ecol. 2012, 21, 2039–2044. [Google Scholar] [CrossRef] [PubMed]
  6. Jerde, C.L.; Mahon, A.R.; Chadderton, W.L.; Lodge, D.M. “Sight-unseen” detection of rare aquatic species using environmental DNA. Conserv. Lett. 2011, 4, 150–157. [Google Scholar] [CrossRef]
  7. Biggs, J.; Ewald, N.; Valentini, A.; Gaboriaud, C. Using eDNA to develop a national citizen science-based monitoring programme for the great crested newt (Trituruscristatus). Biol. Conserv. 2015, 183, 19–28. [Google Scholar] [CrossRef]
  8. Sigsgaard, E.E.; Carl, H.; Møller, P.R.; Thomsen, P.F. Monitoring the near-extinct European weather loach in Denmark based on environmental DNA from water samples. Biol. Conserv. 2015, 183, 46–52. [Google Scholar] [CrossRef] [Green Version]
  9. Takahara, T.; Minamoto, T.; Doi, H. Using environmental DNA to estimate the distribution of an invasivefish species in ponds. PLoS ONE 2013, 8, e56584. [Google Scholar] [CrossRef] [Green Version]
  10. Fujiwara, A.; Matsuhashi, S.; Doi, H.; Yamamoto, S.; Minamoto, T. Use of en-vironmental DNA to survey the distribution of an invasive submerged plant in ponds. Freshw. Sci. 2016, 35, 748–754. [Google Scholar] [CrossRef]
  11. Laramie, M.B.; Pilliod, D.S.; Goldberg, C.S. Characterizing the distribution of an endangered salmonid using environmental DNA analysis. Biol. Conserv. 2015, 183, 29–37. [Google Scholar] [CrossRef] [Green Version]
  12. Thomsen, P.F.; Kielgast, J.; Iversen, L.L.; Wiuf, C.; Rasmussen, M.; Gilbert, M.T.P.; Orlando, L.; Willerslev, E. Monitoring endangered freshwater biodiversity using environmental DNA. Mol. Ecol. 2012, 21, 2565–2573. [Google Scholar] [CrossRef] [PubMed]
  13. Wilcox, T.M.; McKelvey, K.S.; Young, M.K.; Jane, S.F.; Lowe, W.H.; Whiteley, A.R.; Schwartz, M.K. Robust detection of rare species using environmental DNA: The importance of primer specificity. PLoS ONE 2013, 8, e59520. [Google Scholar] [CrossRef] [Green Version]
  14. Lodge, D.M.; Turner, C.R.; Jerde, C.L.; Barnes, M.A.; Chadderton, L.; Egan, S.P.; Feder, J.L.; Mahon, A.R.; Pfrender, M.E. Conservation in a cup of water: Estimating biodiversity and population abundance from environmental DNA. Mol. Ecol. 2012, 21, 2555–2558. [Google Scholar] [CrossRef] [PubMed]
  15. Kelly, R.P.; O’Donnell, J.L.; Lowell, N.C.; Shelton, A.O.; Samhouri, J.F.; Hennessey, S.M.; Feist, B.E.; Williams, G.D. Genetic signatures of ecological diversity along an urbanization gradient. PeerJ 2016, 4, e2444. [Google Scholar] [CrossRef] [Green Version]
  16. Port, J.A.; O’Donnell, J.L.; Romero-Maraccini, O.C.; Leary, P.R.; Litvin, S.Y.; Nickols, K.J.; Yamahara, K.M.; Kelly, R.P. Assessing vertebrate biodiversity in a kelp forest ecosystem using environmental DNA. Mol. Ecol. 2016, 25, 527–541. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  17. Tillotson, M.D.; Kelly, R.P.; Duda, J.J.; Hoy, M.; Kralj, J.; Quinn, T.P. Concentrations of environmental DNA (eDNA) reflect spawning salmon abundance at fine spatial and temporal scales. Biol. Conserv. 2018, 220, 1–11. [Google Scholar] [CrossRef]
  18. Pfleger, M.O.; Rider, S.J.; Johnston, C.E.; Janosik, A.M. Saving the doomed: Using eDNA to aid in detection of rare sturgeon for conservation (Acipenseridae). Glob. Ecol. Conserv. 2016, 8, 99–107. [Google Scholar] [CrossRef] [Green Version]
  19. Gao, S.B.; Tang, H.Y.; Qiao, H.; Yang, Z.; Chen, J.S. A study on the current status of fish resources in the lower reaches of the Jinsha River. J. Water Ecol. 2013, 34, 44–49. [Google Scholar]
  20. Allan, J.D.; Flecker, A.S. Biodiversity conservation in running waters. BioScience 1993, 43, 32–43. [Google Scholar] [CrossRef]
  21. Marchant, R.; Hehir, G. The use of AUSRIVAS predictive models to assess the response of lotic macroinvertebrates to dams in south-east Australia. Freshw. Biol. 2002, 47, 1033–1050. [Google Scholar] [CrossRef]
  22. March, J.G.; Benstead, J.P.; Pringle, C.M.; Scatena, F.N. Damming tropical island streams: Problems, solutions, and alternatives. BioScience 2003, 53, 1069–1078. [Google Scholar] [CrossRef] [Green Version]
  23. Bunn, S.E.; Arthington, A.H. Basic principles and ecological consequences of altered flow regimes for aquatic biodiversity. Environ. Manag. 2002, 30, 492–507. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Cheng, F.; Li, W.; Castello, L.; Murphy, B.R.; Xie, S.G. Potential effects of dam cascade on fish: Lessons from the Yangtze River. Rev. Fish Biol. Fish. 2015, 25, 569–585. [Google Scholar] [CrossRef]
  25. Glowacki, L.B.; Penczak, T. Drivers of fish diversity, homogenization/differentiation and species range expansions at the watershed scale. Divers. Distrib. 2013, 19, 907–918. [Google Scholar] [CrossRef]
  26. Cowx, I.G.; Noble, R.A.; Nunn, A.D.; Bolland, J.; Walton, S.; Peirson, G.; Harvey, J.P. Flow requirements of non-salmonids. Fish. Manag. Ecol. 2012, 19, 548–556. [Google Scholar] [CrossRef]
  27. Rahel, F.J. Homogenization of freshwater faunas. Annu. Rev. Ecol. Syst. 2002, 33, 291–315. [Google Scholar] [CrossRef] [Green Version]
  28. Clavero, M.; Blanco-Garrido, F.; Prenda, J. Fish fauna in Iberian Mediterranean River basins: Biodiversity, introduced species and damming impacts. Aquat. Conserv. Mar. Freshw. Ecosyst. 2004, 14, 575–585. [Google Scholar] [CrossRef] [Green Version]
  29. Nunn, A.D.; Harvey, J.P.; Britton, J.R.; Frear, P.A.; Cowx, I.G. Fish, climate and the Gulf Stream: The influence of abiotic factors on the recruitment success of cyprinid fishes in lowland rivers. Freshw. Biol. 2007, 52, 1576–1586. [Google Scholar] [CrossRef]
  30. Young, P.S.; Cech, J.J.; Thompson, L.C. Hydropower-related pulsed-flow impacts on stream fishes: A brief review, conceptual model, knowledge gaps, and research needs. Rev. Fish Biol. Fish. 2011, 21, 713–731. [Google Scholar] [CrossRef]
  31. Agostinho, A.A.; Gomes, L.C.; Veríssimo, S.; Okada, E.K. Flood regime, dam regulation and fish in the Upper Paraná River: Effects on assemblage attributes, reproduction and recruitment. Rev. Fish Biol. Fish. 2004, 14, 11–19. [Google Scholar] [CrossRef]
  32. He, W.; Lian, J.; Zhang, J.; Yu, X.; Chen, S. Impact of intra-annual runoff uniformity and global warming on the thermal regime of a large reservoir. Sci. Total Environ. 2019, 658, 1085–1097. [Google Scholar] [CrossRef] [PubMed]
  33. He, W.; Lian, J.; Yao, Y.; Wu, M.; Ma, C. Modeling the effect of temperature-control curtain on the thermal structure in a deep stratified reservoir. J. Environ. Manag. 2017, 202, 106–116. [Google Scholar] [CrossRef] [PubMed]
  34. Taberlet, P.; Bonin, A.; Zinger, L.; Coissac, E. Environmental DNA: For Biodiversity Research and Monitoring; Oxford University Press: Oxford, UK, 2018. [Google Scholar]
  35. Gu, Y.H.; Li, D.H.; Wang, C.L. Analysis of fish resources in the reservoir area of Sanbanxi Power Station of Qingshui River. J. Guizhou Norm. Univ. (Nat. Sci. Ed.) 1998, 3, 5–12. [Google Scholar]
  36. Gao, X.; Zeng, Y.; Wang, J.; Liu, H. Immediate impacts of the second impoundment on fish communities in the Three Gorges Reservoir. Environ. Biol. Fishes 2010, 87, 163–173. [Google Scholar] [CrossRef]
  37. Santos, R.M.; Fernandes, L.F.; Cortes, R.M.; Varandas, S.G. Integrative assessment of river damming impacts on aquatic fauna in a Portuguese reservoir. Sci. Total Environ. 2017, 601–602, 1108–1118. [Google Scholar] [CrossRef]
  38. Perera, H.A.C.C.; Li, Z.J.; Silva, S.S.; Zhang, T.L.; Yuan, J.; Ye, S.W.; Xia, Y.G.; Liu, J.S. Effect of the distance from the dam on river fish community structure and compositional trends, with reference to the Three Gorges Dam, Yangtze River, China. Acta Hydrobiol. Sin. 2014, 38, 438–445. [Google Scholar]
  39. Deiner, K.; Walser, J.; Machler, E.; Altermatt, F. Choice of capture and extraction methods affect detection of freshwater biodiversity from environmental DNA. Biol. Conserv. 2015, 183, 53–63. [Google Scholar] [CrossRef]
  40. Strickler, K.M.; Fremier, A.K.; Goldberg, C.S. Quantifying effects of UV-B, temperature, and pH on eDNA degrada tion in aquatic microcosms. Biol. Conserv. 2015, 183, 85–92. [Google Scholar] [CrossRef]
  41. Turner, C.R.; Barnes, M.A.; Xu, C.C.Y.; Jones, S.E.; Christopher, L.J.; Lodge, D.M. Particle size distribution and optimal capture of aqueous macrobial eDNA. Methods Ecol. Evol. 2014, 5, 676–684. [Google Scholar] [CrossRef] [Green Version]
  42. Hinlo, R.; Gleeson, D.; Lintermans, M.; Furlan, E. Methods to maximise recovery of environmental DNA from water samples. PLoS ONE 2017, 12, e0179251. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  43. Dejean, T.; Valentini, A.; Miquel, C.; Taberlet, P.; Bellemain, E.; Miaud, C. Improved detection of an alien invasive species through environmental DNA barcoding: The example of the American bullfrogLithobates catesbeianus. J. Appl. Ecol. 2012, 49, 953–959. [Google Scholar] [CrossRef]
  44. Bohmann, K.; Evans, A.; Gilbert, M.T.P.; Carvalho, G.R.; Creer, S.; Knapp, M.; Yu, D.W.; Bruyn, M. Environmental DNA for wildlife biology and biodiversity monitoring. Trends Ecol. Evol. 2014, 29, 358–367. [Google Scholar] [CrossRef] [PubMed]
  45. Eichmiller, J.J.; Bajer, P.G.; Sorensen, P.W. The relationship between the distribution of common carp and their environ mental DNA in a small lake. PLoS ONE 2014, 9, e112611. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  46. Robson, H.L.; Noble, T.H.; Saunders, R.J.; Robson, S.K.; Burrows, D.W.; Jerry, D.R. Fine-tuning for the tropics: Application of eDNA technology for invasive fish detection in tropical freshwater ecosystems. Mol. Ecol. Resour. 2016, 16, 922–932. [Google Scholar] [CrossRef]
  47. Grill, G.; Lehner, B.; Thieme, M.; Geenen, B.; Tickner, D. Mapping the world’s free-flowing rivers. Nature 2019, 569, 215–221. [Google Scholar] [CrossRef] [PubMed]
  48. Tao, Y.; Wang, Y.; Wang, D.; Wu, J.; Ni, L. Assessing water temperature variations and impacts on fish spawning downstream of Three Gorges dam. J. Hydroelectr. Eng. 2018, 37, 48–55. [Google Scholar]
  49. Dudgeon, D. The ecology of tropical Asian rivers and streams in relation to biodiversity conservation. Ann. Rev. Ecol. Syst. 2000, 31, 239–263. [Google Scholar] [CrossRef]
  50. Dai, Y.G.; Chen, Y.F. Fish fauna and its ecological types of Qingshui River in GuizhouProvince, Southwest China. Chin. J. Ecol. 2007, 26, 682–687. [Google Scholar]
  51. Xie, G.; Chen, K.C.; Hu, Y.C.; Deng, G.C.; Li, C.W. Relationship between embryonic development of Spinibarbus denticulatus and water temperature and salinity. J. Dalian Fish. Univ. 2003, 18, 95–98. [Google Scholar]
Figure 1. Geographic location of sample sites using ArcMap 9.1 software (ESRI Inc., Redlands, CA, USA).
Figure 1. Geographic location of sample sites using ArcMap 9.1 software (ESRI Inc., Redlands, CA, USA).
Sustainability 14 12966 g001
Figure 2. The Venn diagram of OTU distribution (OTU distribution across the 10 sampling points).
Figure 2. The Venn diagram of OTU distribution (OTU distribution across the 10 sampling points).
Sustainability 14 12966 g002
Figure 3. Fish community structure distribution heatmap, different color reflecting the relative abundance of community (%).
Figure 3. Fish community structure distribution heatmap, different color reflecting the relative abundance of community (%).
Sustainability 14 12966 g003
Figure 4. Multiple samples PCA analysis. The difference in the x-axis can explain 74.77% of the full analysis results.
Figure 4. Multiple samples PCA analysis. The difference in the x-axis can explain 74.77% of the full analysis results.
Sustainability 14 12966 g004
Table 1. List of fish species and habitat detected based on environmental DNA metabarcoding and documentary records. Observations are marked with “+” sign.
Table 1. List of fish species and habitat detected based on environmental DNA metabarcoding and documentary records. Observations are marked with “+” sign.
TaxonHabitat TypeDocumentary Records (before Impounding)Detected on eDNA
Cypriniformes
Cobitidae
Leptobotia guilinensisFlowing and benthic fish+
Leptobotia tchangiFlowing and benthic fish+
Parabotia banarescuiFlowing and benthic fish+
Parabotia maculosusFlowing and benthic fish+
Schistura fasciolataFlowing and benthic fish+
Cobitis sinensisFlowing and benthic fish+
Misgurnus anguillicaudatusstill water++
Cyprinidae
Zacco platypusFlowing fish+
Opsariichthys bidensstreaming flow+
Pseudolaubuca sinensisstreaming flow+
Rectoris luxiensisstreaming flow++
Sinilabeo tungtingrapid flow+
Acheilognathus meridianustreaming flow++
Rhodeus ocellatusstreaming flow or still water++
Microphysogobio kiatingensisstreaming flow++
Microphysogobio fukiensisstreaming flow++
Gobiobotia meridionalisstreaming flow+
Pseudorasbora parvastreaming flow+
Sarcocheilichthys kiansiensisstreaming flow+
Cyprinus carpiostreaming flow or still water++
Carassius auratusstreaming flow or still water++
Ctenopharyngodon idellusstreaming flow or still water++
Abbottina rivularisstreaming flow+
Squalidus argentatusstreaming flow or still water++
Squalidus wolterstorffistreaming flow +
Hemibarbusmaculatusstreaming flow or still water++
Hemibarbus labeorapid flow+
Saurogobio dabryistreaming flow+
Platysmacheilus exiguusstreaming flow+
Sinibrama macropsstreaming flow or still water++
Sinibrama taeniatusstreaming flow or still water++
Hemicculter Leuciclusstreaming flow or still water++
Pseudohemiculter hainanensisstreaming flow or still water++
Spinibarbus denticulatusstreaming flow++
Acrossocheilus monticolastreaming flow++
Onychostoma linistreaming flow++
Onychostoma rarastreaming flow++
Folifer brevifilisrapid flow+
Distoechodon tumirostrisstreaming flow+
Xenocypris argenteastreaming flow++
Megalobrama terminalisstreaming flow or still water +
Megalobrama amblycephalastreaming flow or still water +
Hypophthalmichthys molitrixstreaming flow or still water++
Hypophthalmichthys nobilisstreaming flow or still water++
Procypris rabaudistreaming flow +
Cirrhinus mrigala +
Mylopharyngodon piceusstreaming flow++
Chanodichthys mongolicusstreaming flow +
Culter alburnusstreaming flow or still water++
Squaliobarbus curriculusstreaming flow +
Homalopteridae
Sinogastromuzon hisashiensisrapid flow+
Lepturichthys fimbriatarapid flow+
Hemimyzon macropterarapid flow+
Vanmanenia pingchowensisrapid flow+
Catostomidae
Myxocyprinus asiaticusstreaming flow +
Siluriformes +
Siluridae +
Silurus asotusstreaming flow or still water++
Silurus meriordinalisstreaming flow or still water++
Bagridae
Pelteobagrus fulvidracostreaming flow++
Pseudobagrus adiposalisstreaming flow++
Pseudobagrus truncatsustreaming flow++
Leiocassis crassilabrisstreaming flow+
Leiocassis longirostrisstreaming flow++
Hemibagrusmacropterusstreaming flow++
Hemibagrus guttatusstreaming flow +
Sisoridae
Glyptothorax fukiensisrapid flow+
Ictaluridae
Ictalurus punctatus +
Perciformes
Serranidae
Siniperca scherzeristreaming flow++
Siniperca kneristreaming flow++
Siniperca undulatarapid flow+
Siniperca obscurarapid flow+
Coreosiniperca rouleirapid flow+
Channidae +
Channa asiaticastreaming flow or still water++
Eleotridae
Odontobutis obscurusstreaming flow or still water+
Gobiidae
Rhinogobius giurinusstreaming flow or still water++
Cichlidae
Oreochromis niloticus +
Centrarchidae
Micropterus salmoides +
Synbranchiformes
Synbranchidae
Monopterus albuscavernicolous+
Cyprinodontiformes
Poeciliidae
Gambusia affinis +
Table 2. Alpha Diversity Synthesis Analysis.
Table 2. Alpha Diversity Synthesis Analysis.
SampleChao1SimpsonACEShannonCoverage
sbx01779.65150.1204765.81615.31540.999286599
sbx02773.01750.1102747.18494.86420.999147635
sbx03681.64930.0635685.42795.23860.998604079
sbx04760.80.2033763.47244.55170.999425456
sbx05769.51510.1315753.78765.3080.998973452
sbx06764.73910.202753.80154.58410.999160983
sbx07756.15870.2403755.1174.58890.999725285
sbx08736.47050.1865725.09264.73340.999213728
sbx09824.08570.172811.76445.02110.999317919
sbx10729.24190.212723.46294.42050.999445499
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Ma, X.; Yang, H.; Zhong, X.; Zeng, P.; Zhou, X.; Zeng, S.; Dong, X.; Min, W.; Huang, F. eDNA Metabarcoding Analysis of the Composition and Spatial Patterns of Fish Communities in the Sanbanxi Reservoir, China. Sustainability 2022, 14, 12966. https://doi.org/10.3390/su142012966

AMA Style

Ma X, Yang H, Zhong X, Zeng P, Zhou X, Zeng S, Dong X, Min W, Huang F. eDNA Metabarcoding Analysis of the Composition and Spatial Patterns of Fish Communities in the Sanbanxi Reservoir, China. Sustainability. 2022; 14(20):12966. https://doi.org/10.3390/su142012966

Chicago/Turabian Style

Ma, Xiuhui, Hanwei Yang, Xue Zhong, Peng Zeng, Xianjun Zhou, Sheng Zeng, Xianghong Dong, Wenwu Min, and Fujiang Huang. 2022. "eDNA Metabarcoding Analysis of the Composition and Spatial Patterns of Fish Communities in the Sanbanxi Reservoir, China" Sustainability 14, no. 20: 12966. https://doi.org/10.3390/su142012966

APA Style

Ma, X., Yang, H., Zhong, X., Zeng, P., Zhou, X., Zeng, S., Dong, X., Min, W., & Huang, F. (2022). eDNA Metabarcoding Analysis of the Composition and Spatial Patterns of Fish Communities in the Sanbanxi Reservoir, China. Sustainability, 14(20), 12966. https://doi.org/10.3390/su142012966

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