丁洋 李艷艷 趙進(jìn)勇,? 彭文啟 張晶 任錦豪

北京大學(xué)學(xué)報(bào)(自然科學(xué)版) 第60卷 第1期 2024年1月

Acta Scientiarum Naturalium Universitatis Pekinensis, Vol. 60, No. 1 (Jan. 2024)

10.13209/j.0479-8023.2023.080

中國(guó)水科院“城鄉(xiāng)水系生態(tài)景觀構(gòu)建理論和技術(shù)研究創(chuàng)新團(tuán)隊(duì)”項(xiàng)目(WE0145B042021)和水利部重大科技項(xiàng)目(SKR-2022052)資助

2022–12–07;

2023–06–19

基于環(huán)境DNA宏條形碼的漢江上游黃金峽段魚類多樣性研究

丁洋1,2李艷艷3趙進(jìn)勇1,2,?彭文啟1,2張晶1,2任錦豪1

1.中國(guó)水利水電科學(xué)研究院, 北京 100038; 2.水利部京津冀水安全保障重點(diǎn)實(shí)驗(yàn)室, 北京 100038; 3.北京大學(xué)地球與空間科學(xué)學(xué)院, 北京 100871; ?通信作者, E-mail: zhaojy@iwhr.com

2020 年, 采用環(huán)境 DNA 宏條形碼對(duì)漢江上游黃金峽魚類多樣性進(jìn)行研究, 并對(duì)比歷史調(diào)查數(shù)據(jù), 構(gòu)建漢江上游黃金峽段魚類名錄。從 9 個(gè)采樣點(diǎn)中共檢測(cè)出 20 種魚類, 占名錄的 45.45%; 鯉魚()、鯽魚()、圓吻鲴()和銀鮈()為優(yōu)勢(shì)種; 黃金峽上、下游 Shannon 指數(shù)存在顯著性差異(<0.01,=9), 上游魚類 Shannon 指數(shù)明顯大于下游, 黃金峽水利樞紐是造成魚類多樣性空間差異的主要原因。環(huán)境 DNA 檢測(cè)出的魚類組成與過(guò)去傳統(tǒng)方法監(jiān)測(cè)的結(jié)果相近。環(huán)境 DNA 宏條形碼技術(shù)是一種新興的生物監(jiān)測(cè)手段, 可以輔助傳統(tǒng)魚類監(jiān)測(cè)方法, 快速檢測(cè)漢江上游魚類多樣性及其空間分布。

環(huán)境 DNA 宏條形碼; 魚類多樣性; 漢江上游黃金峽段

河流、湖泊和濕地是地球上生物多樣性最豐富的地方之一, 它們以僅占地球表面 1%的面積, 貢獻(xiàn)了全球 50%的魚類多樣性[1]。魚類多樣性是生物多樣性的重要組成部分, 保護(hù)魚類多樣性有利于維持水生態(tài)系統(tǒng)的穩(wěn)定與健康。但是, 目前針對(duì)魚類多樣性的研究缺乏足夠的調(diào)查工作[2]。2021 年, 世界自然基金會(huì)(WWF)和世界自然保護(hù)聯(lián)盟(IUCN)等 16 個(gè)組織聯(lián)合發(fā)布的報(bào)告指出, 已知有 80 種淡水魚類滅絕, 其中僅 2020 年就滅絕 16 種, 1970 年以來(lái), 淡水洄游魚類的數(shù)量下降 76%, 大型淡水魚類的數(shù)量下降 94%[3]。取水、引水以及農(nóng)業(yè)擴(kuò)張和集約化等人類活動(dòng)是對(duì)淡水供應(yīng)量和質(zhì)量的主要威脅因素, 而生境破碎化和流量變化通過(guò)改變河流和洪泛區(qū)湖泊等, 成為對(duì)淡水生物多樣性的主要威 脅[4–5]。在淡水魚類多樣性喪失日益嚴(yán)重的背景下, 掌握魚類多樣性現(xiàn)狀變得尤為重要[6]。

傳統(tǒng)的魚類調(diào)查方法(拖網(wǎng)、圍網(wǎng)、電捕魚、釣魚、水聲和視覺(jué)等)耗時(shí)費(fèi)力, 價(jià)格昂貴, 可能對(duì)魚體造成傷害, 并對(duì)調(diào)查人員具有較高的分類學(xué)知識(shí)要求[7]。目前, 基于環(huán)境 DNA 宏條形碼(environ-mental DNA metabarcoding, eDNA metabarcoding)技術(shù)的物種檢測(cè)徹底改變了生物多樣性調(diào)查的方式, 被證明是一種高效、經(jīng)濟(jì)和非侵入性的生物監(jiān)測(cè)方法[8]。環(huán)境 DNA 宏條形碼技術(shù)通過(guò)直接從環(huán)境樣本(例如沉積物、糞便、水、土壤或空氣)中提取DNA, 利用識(shí)別魚類種群的通用引物進(jìn)行 PCR 擴(kuò)增和高通量測(cè)序, 確定魚類的種類和相對(duì)豐度, 并計(jì)算不同區(qū)域的魚類多樣性[9–10]。監(jiān)測(cè)過(guò)程無(wú)需采集目標(biāo)生物, 極大地簡(jiǎn)化了魚類檢測(cè)步驟, 彌補(bǔ)了傳統(tǒng)形態(tài)學(xué)監(jiān)測(cè)的不足[11]。環(huán)境 DNA 宏條形碼在魚類物種多樣性分析、外來(lái)入侵魚類物種足跡追蹤和瀕危珍稀魚類資源調(diào)查中表現(xiàn)出明顯的優(yōu)勢(shì)和巨大的潛能。Sales 等[7]利用環(huán)境 DNA 宏條形碼, 監(jiān)測(cè)和分析巴西熱基蒂尼奧尼亞河流域魚類群落的時(shí)空動(dòng)態(tài)。Robson 等[12]利用環(huán)境 DNA 宏條形碼檢測(cè)熱帶入侵魚類物種。Song 等[13]利用環(huán)境 DNA 宏條形碼, 定量地計(jì)算水流特征對(duì)亞洲鯉魚入侵檢測(cè)概率的影響。國(guó)內(nèi)利用環(huán)境 DNA 宏條形碼開展魚類多樣性調(diào)查的研究目前較少, 但也針對(duì)魚類資源、物種多樣性及瀕危珍稀魚類資源調(diào)查開展了一些研究。Wang 等[14]利用環(huán)境 DNA 宏條形碼評(píng)價(jià)我國(guó)東海大黃魚資源。舒璐等[15]基于環(huán)境 DNA 宏條形碼研究洱海魚類多樣性。吳昀晟等[16]對(duì)比環(huán)境 DNA宏條形碼與傳統(tǒng)調(diào)查法對(duì)長(zhǎng)江江豚的檢測(cè)結(jié)果, 證明環(huán)境 DNA 宏條形碼在長(zhǎng)江江豚監(jiān)測(cè)中不僅具有較高的準(zhǔn)確性, 還具有更高的靈敏性。

漢江是長(zhǎng)江最大的支流, 上游緊依中國(guó)南北分界線(秦嶺–淮河線)的南界, 是重要的水源涵養(yǎng)生態(tài)功能保護(hù)區(qū), 也是中國(guó)南水北調(diào)中線工程和引漢濟(jì)渭工程的水源地[17–20]。典型的山地生態(tài)系統(tǒng)決定了漢江上游生態(tài)環(huán)境的敏感性和脆弱性, 極易受到人為干擾和氣候變化的影響[21]。目前針對(duì)漢江魚類的調(diào)查研究多集中在中下游地區(qū)[22–24], 對(duì)漢江上游魚類多樣性的調(diào)查研究鮮有報(bào)道。隨著長(zhǎng)江流域重點(diǎn)水域“十年禁漁”措施的全面啟動(dòng)[25], 漢江作為禁捕水域之一, 傳統(tǒng)的魚類調(diào)查方法不再適用于該區(qū)域。

基于上述背景, 本研究首次利用環(huán)境 DNA 宏條形碼, 對(duì)漢江上游黃金峽段魚類多樣性的空間分布特征進(jìn)行調(diào)查分析, 旨在探討環(huán)境 DNA 宏條形碼技術(shù)在檢測(cè)漢江上游魚類多樣性可行性, 為漢江上游魚類多樣性監(jiān)測(cè)及生物資源估算提供新的技術(shù)手段。

1 材料與方法

1.1 生物監(jiān)測(cè)點(diǎn)位選擇

漢江發(fā)源于陜西省寧強(qiáng)縣, 自西向東流經(jīng)陜西省內(nèi)漢中、安康兩市, 并于安康市白河縣流出陜西省。漢江上游流域位于亞熱帶季風(fēng)氣候區(qū), 年均氣溫為 13.2℃, 多年平均降雨量為 895mm[26]。本次調(diào)查區(qū)域位于引漢濟(jì)渭水源地——黃金峽水利樞紐上游及下游, 涉及漢中市洋縣和西鄉(xiāng)縣。于 2020 年11 月 13 日至 16 日進(jìn)行采樣, 分別選取 9 個(gè)采樣點(diǎn), 樣點(diǎn)編號(hào)及位置見(jiàn)圖 1。其中, 在黃金峽水利樞紐壩上布設(shè) 5 個(gè)采樣點(diǎn)(S1～S5), 壩下布設(shè) 4 個(gè)采樣點(diǎn)(S6～S9); 漢江干流全長(zhǎng)約 100km, 設(shè)置 6 個(gè)樣點(diǎn), 漢江支流黨水河、金水河及子午河各設(shè)一個(gè)樣點(diǎn)。

1.2 環(huán)境DNA樣品采集及環(huán)境DNA測(cè)定

使用 2.5L 的采水器, 在每個(gè)樣點(diǎn)采集表層水和底層水各 1L, 混合后保存至已消毒的 2L 廣口瓶中。所有樣品均在 24 小時(shí)內(nèi)使用 0.45μm 的混合纖維素濾膜(MCE; Whatman公司)進(jìn)行真空抽濾, 玻璃抽濾漏斗在每次抽濾前用 10%的次氯酸鈉消毒液浸泡30 分鐘, 并用純凈水沖洗干凈, 防止樣品間的交叉污染。為評(píng)估是否存在外源 DNA 污染, 每次過(guò)濾設(shè)置一個(gè)陰性對(duì)照。樣品過(guò)濾后, 富集環(huán)境DNA 的濾膜分別裝入 5mL 的離心管中, 于?20℃保存, 并盡快送往實(shí)驗(yàn)室進(jìn)行 DNA 提取。濾膜送往上海美吉生物醫(yī)藥科技有限公司進(jìn)行環(huán)境 DNA測(cè)定、PCR 擴(kuò)增及高通量測(cè)序。

1.3 歷史數(shù)據(jù)調(diào)查

通過(guò)對(duì)照洋縣水利局魚類資源數(shù)據(jù), 整理得到漢江上游黃金峽段魚類名錄。

2 結(jié)果與分析

2.1 物種組成

通過(guò)環(huán)境 DNA 宏條形碼檢測(cè), 從 9 個(gè)采樣點(diǎn)的樣品中共得到有效拼接序列 542400 條, 注釋出魚類20 種, 隸屬 3 目 7 科 18 屬(表 1)。各采樣點(diǎn)魚類物種及其相對(duì)序列豐度見(jiàn)圖 2, 采樣點(diǎn)與魚類物種關(guān)系如圖 3 所示。各采樣點(diǎn)的魚類物種組成存在明顯的差異, 其中鯉科(如鯽魚、鯉魚和花魚骨)在各采樣點(diǎn)均顯示出相對(duì)較高的序列豐度。

2.2 歷史調(diào)查數(shù)據(jù)

通過(guò)收集洋縣水利局魚類資源數(shù)據(jù), 統(tǒng)計(jì)得到魚類 37 種, 屬于 4 目 7 科 30 屬。其中, 鯉形目物種數(shù)最多, 有 30 種, 占比為 81.08%; 其次是鱸形目, 有 5 種, 占比為 13.51%; 鲇形目有兩種, 占比為5.41%; 鯰形目有一種, 占比為 2.70%。結(jié)合環(huán)境DNA 宏條碼的檢測(cè)結(jié)果, 最終統(tǒng)計(jì)得到魚類 44 種, 屬于 4 目 8 科, 整理得到漢江上游黃金峽段魚類名錄, 見(jiàn)附錄(訪問(wèn) http://xbna.pku.edu.cn 查看附錄)。

從圖 4 可以看出, 由歷史調(diào)查數(shù)據(jù)和環(huán)境 DNA宏條形碼共同檢測(cè)出的魚類共 13 種, 占總種數(shù)的29.55%; 僅由歷史調(diào)查數(shù)據(jù)統(tǒng)計(jì)得到的魚類共 24種, 占總種數(shù)的 54.54%; 僅由環(huán)境 DNA 宏條形碼檢測(cè)出的魚類共 7 種, 占總種數(shù)的 15.91%。

僅由環(huán)境 DNA 宏條形碼檢測(cè)出的魚類為激流馬口鱲(_sp._ZF-11309)、福建小鰾鮈(_)、紅尾副鰍(_)、圓尾擬鲿(_)、大鰭鳠(_)、神農(nóng)吻蝦虎魚(-)和真吻鰕虎魚(_)。

2.3 物種空間分布

根據(jù)研究區(qū)位置, 將采樣點(diǎn)分為黃金峽水利樞紐上游區(qū)域和下游區(qū)域(簡(jiǎn)稱黃金峽上游、黃金峽下游)。t 檢驗(yàn)方法顯示(圖 5), 黃金峽上下游 Shannon指數(shù)存在顯著性差異(<0.01,=9), 黃金峽上游魚類 Shannon 指數(shù)明顯大于下游?；诓蓸狱c(diǎn)中魚類物種序列豐度數(shù)據(jù), 利用 Wilcoxon 秩和檢驗(yàn)方法, 檢驗(yàn)黃金峽水利樞紐上下游魚類物種差異的顯著性。由圖 6 可以看出, 黃金峽上游及下游魚類物種差異明顯。黃金峽下游魚類監(jiān)測(cè)物種均為鯉科, 以圓吻鲴(41.56%)、鯉魚(32.65%)和銀鮈(25%)為主; 黃金峽上游魚類監(jiān)測(cè)物種以鯽魚(34.32%)、鯉魚(13.92%)和波氏吻鰕虎魚(11.09%)為主, 且鲇魚、斷線麥穗魚、激流馬口鱲、馬口魚、花魚骨、細(xì)紋頜須鮈和紅尾副鰍等魚類物種均有涉及。

圖2 各采樣點(diǎn)魚類物種相對(duì)序列豐度

圖3 各采樣點(diǎn)與魚類物種關(guān)系Circos圖

表1 漢江上游黃金峽段9個(gè)采樣點(diǎn)環(huán)境DNA宏條形碼注釋魚類物種

灰色區(qū)域代表本研究環(huán)境 DNA 宏條碼所檢測(cè)物種, 白色區(qū)域代表歷史調(diào)查數(shù)據(jù)

圖5 黃金峽上下游 Shannon 指數(shù)組間差異

圖6 各采樣點(diǎn)魚類物種空間差異

3 討論

3.1 基于環(huán)境DNA宏條形碼的漢江上游黃金峽段魚類多樣性檢測(cè)

本研究運(yùn)用環(huán)境 DNA 宏條形碼技術(shù)分析漢江上游黃金峽段魚類多樣性, 共檢測(cè)出 20 種魚類, 魚類物種以鯉科居多, 共 12 種, 占魚類物種總數(shù)的 60%, 其中鯉魚、鯽魚、圓吻鲴和銀鮈為優(yōu)勢(shì)種。對(duì)比漢江上游黃金峽段魚類名錄, 在 9 個(gè)采樣點(diǎn)位利用環(huán)境 DNA 宏條碼檢測(cè)出 20 種魚類, 占名錄的 45.45%, 表明環(huán)境 DNA 宏條碼具有較高的準(zhǔn)確率。此外, 有 7 種魚類未出現(xiàn)在歷史統(tǒng)計(jì)數(shù)據(jù)中, 如激流馬口鱲(spZF-11309)、福建小鰾鮈()和紅尾副鰍()等小型魚類?？赡苁怯捎跐h江上游黃金峽段魚類調(diào)查數(shù)據(jù)十分稀少, 導(dǎo)致本研究魚類目錄收集不完整, 因此部分檢測(cè)物種未出現(xiàn)在歷史調(diào)查數(shù)據(jù)中。環(huán)境 DNA 宏條形碼技術(shù)能克服傳統(tǒng)調(diào)查方法的一些局限性, 對(duì)傳統(tǒng)捕撈方法進(jìn)行補(bǔ)充觀測(cè), 從而減少傳統(tǒng)調(diào)查方法觀測(cè)不到位導(dǎo)致的生物多樣性下降情況。

1980—1982 年, 陜西省動(dòng)物研究所等機(jī)構(gòu)的研究人員在石泉大壩以上干支流共采集到 61 種魚類, 種類組成同樣以鯉科魚類為主, 占物種總數(shù)的63.93%[27]。王曉臣[28]在 2010—2012 年開展整個(gè)漢江陜西段的魚類調(diào)查, 在石泉大壩以上干支流共采集 63 種魚類, 種類組成以鯉科魚類為主, 占物種總數(shù)的 64.52%, 以棒花魚、麥穗魚、寬鰭眠、馬口魚及鲇魚為優(yōu)勢(shì)物種。采樣點(diǎn)涉及的河段為歷史調(diào)查河段的一部分, 可能是本研究檢測(cè)出的魚類種數(shù)不如傳統(tǒng)漁獲物方法監(jiān)測(cè)的魚類種數(shù)多的主要原因, 但檢測(cè)出的魚類組成與過(guò)去傳統(tǒng)方法監(jiān)測(cè)的結(jié)果相近, 并均顯示魚類以鯉科為主, 占物種總數(shù)的 60%左右。盡管低豐度魚類物種的 DNA 難以采集, 但可以通過(guò)增加采樣點(diǎn)和采樣次數(shù)等方式, 提高檢測(cè)率[29]。

3.2 漢江上游黃金峽段魚類空間分布分析

黃金峽下游檢測(cè)的魚類物種僅有鯉科, 其中圓吻鲴、鯉魚和銀鮈等靜水性魚類為優(yōu)勢(shì)物種?，F(xiàn)場(chǎng)調(diào)研發(fā)現(xiàn), 黃金峽下游取樣點(diǎn)為石門水庫(kù)庫(kù)區(qū)回水區(qū), 屬靜水區(qū)域, 水面寬, 流速小且穩(wěn)定, 水位較高, 水流對(duì)河岸坡棲息地的沖蝕能力較弱。受水庫(kù)回水的影響, 原有的河道失去急流、淺灘和較大的蜿蜒度[30], 破壞了喜急流性魚類的棲息環(huán)境, 因而黃金峽下游中急流性魚類數(shù)目減少, 靜水性魚類數(shù)目相應(yīng)地增加, 以圓吻鲴、鯉魚和銀鮈為主的靜水魚類在黃金峽下游區(qū)域占主導(dǎo)地位, 魚類多樣性較低。黃金峽上游檢測(cè)魚類物種包括鯉科、鰍科、鲿科、鲇科、鱧科和鰕虎魚科, 涵蓋所有檢測(cè)出的魚類科目, 魚類多樣性高。其中, S1 和 S5 采樣點(diǎn)分別位于漢江支流黨水河和金水河, 為山區(qū)溪流河道, 受地勢(shì)起伏及河流蜿蜒等因素影響, 水流流態(tài)多樣, 生境復(fù)雜多樣, 底質(zhì)條件顯著差異, 地貌單元復(fù)雜多樣[31], 因而這兩處采樣點(diǎn)檢測(cè)出的魚類多樣性高, 且種類也包括急流性魚類及靜水性魚類。黃金峽水利樞紐的建設(shè)及截流, 在很大程度上阻斷了上下游的魚類活動(dòng), 影響魚類壩上與壩下的遷移交流, 同時(shí)大壩極大地影響魚類的生存環(huán)境(改變棲息地的物理和生態(tài)特征, 如改變水流、營(yíng)養(yǎng)動(dòng)態(tài)、水質(zhì)和溫度)[32–34]。水利工程的建設(shè)是造成魚類多樣性及物種的差異的主要原因。

3.3 漢江上游環(huán)境下段環(huán)境DNA宏條形碼技術(shù)對(duì)比原始調(diào)查方法的優(yōu)劣性探討

與傳統(tǒng)的魚類監(jiān)測(cè)方法相比, 環(huán)境 DNA 宏條形碼技術(shù)在魚類物種鑒定及多樣性評(píng)估等方面具有明顯的優(yōu)勢(shì)。1)非破壞性采樣: 少量的水樣采集不會(huì)危害當(dāng)?shù)貫l危魚類, 也不會(huì)破壞當(dāng)?shù)厣鷳B(tài)環(huán)境; 2)無(wú)采樣限制: 不受禁漁期和禁漁區(qū)的限制, 可以根據(jù)研究需求設(shè)置監(jiān)測(cè)點(diǎn)位; 3)高效便捷, 節(jié)約成本。環(huán)境 DNA 宏條形碼技術(shù)使用 2 L 左右的水樣便可實(shí)現(xiàn)多個(gè)分類單元的快速鑒定, 可減少船只租賃或購(gòu)買漁獲物的費(fèi)用, 無(wú)需依賴形態(tài)學(xué)和分類學(xué)的專業(yè)知識(shí)[35–36]。

在實(shí)際分析中, 需要考慮采樣點(diǎn)上游環(huán)境 DNA輸移的影響。受環(huán)境 DNA 輸移的影響, 某個(gè)采樣點(diǎn)的魚類檢測(cè)結(jié)果不一定意味著在采集時(shí)魚類個(gè)體就在此處, 可能導(dǎo)致對(duì)采樣點(diǎn)魚類豐富度的高估。因此, 本研究位于河口的采樣點(diǎn)(S1 和 S5)檢測(cè)出的高魚類多樣性需要進(jìn)一步驗(yàn)證。河口處受到支流上游來(lái)水和干流頂托來(lái)水的雙重影響, 流速減緩, 泥沙沉積明顯, 導(dǎo)致環(huán)境 DNA 可能存在于水體中, 也可能與表層沉積物結(jié)合, 同時(shí)受環(huán)境 DNA 的輸移和積累影響, 較高濃度的環(huán)境 DNA 可能有助于沉積物中 DNA 分子的再懸浮, 從而影響在空間和時(shí)間尺度上魚類物種及多樣性的推斷。然而, 由于在漢江上游水環(huán)境中還沒(méi)有進(jìn)行過(guò)有關(guān)環(huán)境 DNA 輸移和擴(kuò)散的研究, 很難就這一問(wèn)題得出可靠的結(jié)論。因此, 開展環(huán)境 DNA 的時(shí)空動(dòng)力學(xué)研究有助于更精確地使用環(huán)境 DNA 宏條形碼技術(shù)來(lái)檢測(cè)魚類多樣性。

4 結(jié)論

本研究首次使用環(huán)境 DNA 宏條形碼技術(shù)評(píng)估漢江上游黃金峽段魚類多樣性, 證明了該技術(shù)監(jiān)測(cè)漢江上游魚類物種組成和魚類空間分布情況的可行性。環(huán)境 DNA 方法可作為評(píng)估研究區(qū)魚類多樣性可靠的補(bǔ)充方法, 為漁業(yè)資源的調(diào)查和保護(hù)提供新的視角。本文結(jié)論如下: 1)通過(guò)環(huán)境 DNA 宏條形碼技術(shù)和歷史調(diào)查數(shù)據(jù)構(gòu)建漢江上游黃金峽段魚類名錄, 環(huán)境 DNA 宏條形碼技術(shù)檢測(cè)出 20 種魚類, 占名錄的 45.45%, 魚類物種以鯉科居多, 占魚類物種總數(shù)的 60%, 其中鯉魚、鯽魚、圓吻鲴和銀鮈為優(yōu)勢(shì)種; 2)與黃金峽下游相比, 黃金峽上游魚類多樣性較高, 黃金峽水利樞紐是造成魚類多樣性空間差異的主要原因; 3)在今后魚類監(jiān)測(cè)過(guò)程中, 應(yīng)開展環(huán)境 DNA 的時(shí)空動(dòng)力學(xué)研究, 可減少環(huán)境 DNA 宏條形碼技術(shù)檢測(cè)的不確定性。

[1] Barbarossa V, Bosmans J, Wanders N, et al. Threats of global warming to the world’s freshwater fishes. Nature Communications, 2021, 12: 1–10

[2] 徐念, 熊美華, 邵科, 等. 長(zhǎng)江中下游環(huán)境DNA宏條形碼生物多樣性檢測(cè)技術(shù)初步研究. 環(huán)境科學(xué)研究, 2020, 33(5): 1187–1196

[3] World Wide Fund for Nature, International Union for Conservation of Nature, Alliance for Freshwater Life, et al. The World’s Forgotten Fishes. Gran: WWF, 2021

[4] Su G H, Logez M, Xu J, et al. Human impacts on global freshwater fish biodiversity. Science, 2021, 371: 835–838

[5] Albert J S, Destouni G, Sylvester S M, et al. Scientists’ warning to humanity on the freshwater biodiversity crisis. Ambio, 2021, 50(1): 85–94

[6] He D K, Sui X Y, Sun H Y, et al. Diversity, pattern and ecological drivers of freshwater fish in China and adjacent areas. Reviews in Fish Biology and Fisheries, 2020, 30: 387–404

[7] Sales N G, Wangensteen O S, Carvalho D C, et al. Space-time dynamics in monitoring neotropical fish communities using eDNA metabarcoding. Science of The Total Environment, 2020, 754: 142096

[8] Zhang S, Zhao J D, Yao M. A comprehensive and comparative evaluation of primers for metabarcoding eDNA from fish. Methods in Ecology and Evolution, 2020, 11: 1609–1625

[9] Wang S P, Yan Z G, H?nfling B, et al. Methodology of fish eDNA and its applications in ecology and envi-ronment. Science of The Total Environment, 2021, 755 (2): 142622

[10] Zhang X W. Environmental DNA shaping a new era of ecotoxicological research. Environmental Science & Technology, 2019, 53(10): 5605–5612

[11] Thomsen P F, Willerslev E. Environmental DNA — an emerging tool in conservation for monitoring past and present biodiversity. Biological Conservation, 2015, 183: 4–18

[12] Robson H, Noble T H, Saunders R J, et al. Fine-tuning for the tropics: application of eDNA technology for invasive fish detection in tropical freshwater ecosy-stems. Molecular Ecology Resources, 2016, 16(4): 922–932

[13] Song J W, Small M J, Casman E A. Making sense of the noise: the effect of hydrology on silver carp eDNA detection in the Chicago area waterway system. Science of the Total Environment, 2017, 605/606: 713–720

[14] Wang X Y, Lu G Q, Zhao L L, et al. Assessment of fishery resources using environmental DNA: the large yellow croaker () in the East China Sea. Fisheries Research, 2021, 235: 105813

[15] 舒璐, 林佳艷, 徐源, 等. 基于環(huán)境DNA宏條形碼的洱海魚類多樣性研究. 水生生物學(xué)報(bào), 2020, 44 (5): 1080–1086

[16] 吳昀晟, 唐永凱, 李建林, 等. 環(huán)境DNA在長(zhǎng)江江豚監(jiān)測(cè)中的應(yīng)用. 中國(guó)水產(chǎn)科學(xué), 2019, 26(1): 124–132

[17] Liu X, Zhu X, Pan Y, et al. Vegetation dynamics in Qinling-Daba Mountains in relation to climate factors between 2000 and 2014. Journal of Geographical Sciences, 2016, 26(1): 45–58

[18] 夏軍, 馬協(xié)一, 鄒磊, 等. 氣候變化和人類活動(dòng)對(duì)漢江上游徑流變化影響的定量研究. 南水北調(diào)與水利科技, 2017, 15(1): 1–6

[19] 王鵬濤, 張立偉, 李英杰, 等. 漢江上游生態(tài)系統(tǒng)服務(wù)權(quán)衡與協(xié)同關(guān)系時(shí)空特征. 地理學(xué)報(bào), 2017, 72(11): 2064–2078

[20] 李瑛. 基于引嘉濟(jì)漢–引漢濟(jì)渭的復(fù)雜跨流域調(diào)水工程協(xié)同調(diào)度研究[D]. 西安: 西安理工大學(xué), 2020

[21] Chang Y, Hou K, Wu Y, et al. A conceptual framework for establishing the index system of ecological envi-ronment evaluation: a case study of the upper Han-jiang River, China. Ecological Indicators, 2019, 107: 105568

[22] 汪登強(qiáng), 高雷, 段辛斌, 等. 漢江下游魚類早期資源及梯級(jí)聯(lián)合生態(tài)調(diào)度對(duì)魚類繁殖影響的初步分析. 長(zhǎng)江流域資源與環(huán)境, 2019, 28(8): 1909–1917

[23] 秦?zé)@, 陳君, 向芳. 漢江中下游梯級(jí)開發(fā)對(duì)產(chǎn)漂流性卵魚類繁殖的影響. 環(huán)境科學(xué)與技術(shù), 2014, 37 (增刊 2): 501–506

[24] 謝文星, 黃道明, 謝山, 等. 丹江口水利樞紐興建后漢江中下游四大家魚等早期資源及其演變. 水生態(tài)學(xué)雜志, 2009, 30(2): 44–49

[25] 李儲(chǔ)信, 朱堯虎. 長(zhǎng)江流域即將全面禁漁 10 年. 生態(tài)經(jīng)濟(jì), 2019, 35(12): 9–12

[26] 侯易明, 潘保柱, 蔣小明, 等. 漢江上游干流和秦嶺南麓典型支流的底棲動(dòng)物群落特征及水質(zhì)生物評(píng)價(jià). 湖泊科學(xué), 2020, 32(4): 1140–1153

[27] 陜西省動(dòng)物研究所, 陜西省科學(xué)院水生生物研究所, 蘭州大學(xué). 秦嶺魚類志. 北京: 科學(xué)出版社, 1987

[28] 王曉臣. 漢江(陜西段)魚類種類與群落結(jié)構(gòu)組成的時(shí)空變化研究[D]. 西安: 西北大學(xué), 2013

[29] Shu L, Ludwig A, Peng Z. Standards for methods utilizing environmental DNA for detection of fish species. Genes, 2020, 11(3): 296

[30] 朱瑤. 大壩對(duì)魚類棲息地的影響及評(píng)價(jià)方法述評(píng). 中國(guó)水利水電科學(xué)研究院學(xué)報(bào), 2005(2): 100–103

[31] 卜倩婷, 李獻(xiàn), 朱仁, 等. 低頭壩驅(qū)動(dòng)山區(qū)溪流局域棲息地和魚類群落的同質(zhì)化. 生物多樣性, 2017, 25(8): 830–839

[32] Shao X J, Fang Y, Jawitz J W, et al. River network connectivity and fish diversity. Science of the Total Environment, 2019, 689: 21–30

[33] 陳求穩(wěn), 張建云, 莫康樂(lè), 等. 水電工程水生態(tài)環(huán)境效應(yīng)評(píng)價(jià)方法與調(diào)控措施. 水科學(xué)進(jìn)展, 2020, 31(5): 793–810

[34] Covino T. Hydrologic connectivity as a framework for understanding biogeochemical flux through water-sheds and along fluvial networks. Geomorphology, 2017, 227: 133–144

[35] 張麗娟, 徐杉, 趙崢, 等. 環(huán)境DNA宏條形碼監(jiān)測(cè)湖泊真核浮游植物的精準(zhǔn)性. 環(huán)境科學(xué), 2021, 42 (2): 796–807

[36] 姜維, 王啟軍, 鄧捷, 等. 以川陜哲羅鮭為目標(biāo)物種的水樣環(huán)境DNA分析流程的優(yōu)化. 應(yīng)用生態(tài)學(xué)報(bào), 2016, 27(7): 2372–2378

Investigating Fish Diversity in Huangjinxia Section of the Upper Hanjiang River Based on Environmental DNA Metabarcoding

DING Yang1,2, LI Yanyan3, ZHAO Jinyong1,2,?, PENG Wenqi1,2, ZHANG Jing1,2, REN Jinhao1

1. China Institute of Water Resources and Hydropower Research, Beijing 100038; 2. Key Laboratory of Water Safety for Beijing-Tianjin-Hebei Region of Ministry of Water Resources, Beijing 100038; 3. School of Earth and Space Sciences, Peking University, Beijing 100871; ? Corresponding author, E-mail: zhaojy@iwhr.com

In 2020, Environmental DNA Metabarcoding was used to study the fish diversity in the Huangjinxia section of the upper Hanjiang River, and the historical survey data were compared to construct a fish list in the Huangjinxia section of the upper Hanjiang River. A total of 20 fish species were detected from 9 sampling sites, accounting for 45.45% of the list;,,andare the dominant species. The Shannon index of fish in the upper and lower reaches of the Huangjinxia was significantly different (< 0.01,=9), and the Shannon index of fish in the upper reaches was significantly larger than that in the lower reaches. The main reason for the spatial differences in fish diversity is the Huangjinxia hydraulic hub. The fish composition detected by environmental DNA was similar to that obtained by traditional methods. As an emerging biomonitoring tool, environmental DNA macrobarcoding technology can complement traditional fish monitoring methods to rapidly detect fish diversity and its spatial distribution in the upper reaches of the Hanjiang River.

environmental DNA metabarcoding; fish diversity; Huangjinxia section of the upper Hanjiang River