• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

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

    2024-01-30 02:47:40丁洋李艷艷趙進勇彭文啟張晶任錦豪
    北京大學學報(自然科學版) 2024年1期
    關(guān)鍵詞:物種環(huán)境檢測

    丁洋 李艷艷 趙進勇,? 彭文啟 張晶 任錦豪

    北京大學學報(自然科學版) 第60卷 第1期 2024年1月

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

    10.13209/j.0479-8023.2023.080

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

    2022–12–07;

    2023–06–19

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

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

    1.中國水利水電科學研究院, 北京 100038; 2.水利部京津冀水安全保障重點實驗室, 北京 100038; 3.北京大學地球與空間科學學院, 北京 100871; ?通信作者, E-mail: zhaojy@iwhr.com

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

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

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

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

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

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

    1 材料與方法

    1.1 生物監(jiān)測點位選擇

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

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

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

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

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

    2 結(jié)果與分析

    2.1 物種組成

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

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

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

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

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

    2.3 物種空間分布

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

    圖2 各采樣點魚類物種相對序列豐度

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

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

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

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

    圖6 各采樣點魚類物種空間差異

    3 討論

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

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

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

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

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

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

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

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

    4 結(jié)論

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

    [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] 徐念, 熊美華, 邵科, 等. 長江中下游環(huán)境DNA宏條形碼生物多樣性檢測技術(shù)初步研究. 環(huán)境科學研究, 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宏條形碼的洱海魚類多樣性研究. 水生生物學報, 2020, 44 (5): 1080–1086

    [16] 吳昀晟, 唐永凱, 李建林, 等. 環(huán)境DNA在長江江豚監(jiān)測中的應(yīng)用. 中國水產(chǎn)科學, 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é)一, 鄒磊, 等. 氣候變化和人類活動對漢江上游徑流變化影響的定量研究. 南水北調(diào)與水利科技, 2017, 15(1): 1–6

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

    [20] 李瑛. 基于引嘉濟漢–引漢濟渭的復(fù)雜跨流域調(diào)水工程協(xié)同調(diào)度研究[D]. 西安: 西安理工大學, 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] 汪登強, 高雷, 段辛斌, 等. 漢江下游魚類早期資源及梯級聯(lián)合生態(tài)調(diào)度對魚類繁殖影響的初步分析. 長江流域資源與環(huán)境, 2019, 28(8): 1909–1917

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

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

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

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

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

    [28] 王曉臣. 漢江(陜西段)魚類種類與群落結(jié)構(gòu)組成的時空變化研究[D]. 西安: 西北大學, 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] 朱瑤. 大壩對魚類棲息地的影響及評價方法述評. 中國水利水電科學研究院學報, 2005(2): 100–103

    [31] 卜倩婷, 李獻, 朱仁, 等. 低頭壩驅(qū)動山區(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), 張建云, 莫康樂, 等. 水電工程水生態(tài)環(huán)境效應(yīng)評價方法與調(diào)控措施. 水科學進展, 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)測湖泊真核浮游植物的精準性. 環(huán)境科學, 2021, 42 (2): 796–807

    [36] 姜維, 王啟軍, 鄧捷, 等. 以川陜哲羅鮭為目標物種的水樣環(huán)境DNA分析流程的優(yōu)化. 應(yīng)用生態(tài)學報, 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

    猜你喜歡
    物種環(huán)境檢測
    吃光入侵物種真的是解決之道嗎?
    英語世界(2023年10期)2023-11-17 09:18:18
    “不等式”檢測題
    “一元一次不等式”檢測題
    “一元一次不等式組”檢測題
    長期鍛煉創(chuàng)造體內(nèi)抑癌環(huán)境
    一種用于自主學習的虛擬仿真環(huán)境
    孕期遠離容易致畸的環(huán)境
    回首2018,這些新物種值得關(guān)注
    環(huán)境
    電咖再造新物種
    汽車觀察(2018年10期)2018-11-06 07:05:26
    在线观看www视频免费| 窝窝影院91人妻| 最新美女视频免费是黄的| 亚洲激情在线av| 全区人妻精品视频| 可以在线观看的亚洲视频| 国产精品亚洲av一区麻豆| 国产黄色小视频在线观看| 极品教师在线免费播放| 国产一区二区激情短视频| 国产午夜精品久久久久久| 国产在线观看jvid| 少妇人妻一区二区三区视频| 久久九九热精品免费| 巨乳人妻的诱惑在线观看| 老汉色∧v一级毛片| 亚洲专区国产一区二区| 国产av又大| 十八禁人妻一区二区| 精品欧美国产一区二区三| 久久九九热精品免费| 亚洲人与动物交配视频| 国内精品久久久久久久电影| 又爽又黄无遮挡网站| 亚洲avbb在线观看| 久久午夜综合久久蜜桃| 精品国内亚洲2022精品成人| 日韩欧美在线乱码| 欧美 亚洲 国产 日韩一| 欧美久久黑人一区二区| 老汉色av国产亚洲站长工具| 成人国产一区最新在线观看| 免费看a级黄色片| 久久婷婷成人综合色麻豆| 国产99白浆流出| 欧美日韩一级在线毛片| 国产成人系列免费观看| 日韩欧美 国产精品| 欧美日韩精品网址| av福利片在线| 免费观看人在逋| 91字幕亚洲| 久久国产精品人妻蜜桃| 男插女下体视频免费在线播放| 岛国在线免费视频观看| 亚洲黑人精品在线| 一级片免费观看大全| 一本久久中文字幕| 在线看三级毛片| 一级毛片高清免费大全| 日本免费一区二区三区高清不卡| 国产av不卡久久| 给我免费播放毛片高清在线观看| 亚洲天堂国产精品一区在线| 久久国产精品影院| 久久99热这里只有精品18| 大型av网站在线播放| 少妇的丰满在线观看| 91九色精品人成在线观看| 精品电影一区二区在线| 视频区欧美日本亚洲| 亚洲成人久久爱视频| 国产探花在线观看一区二区| 麻豆久久精品国产亚洲av| 久久精品国产清高在天天线| 1024视频免费在线观看| 黄色丝袜av网址大全| av免费在线观看网站| 亚洲美女视频黄频| 一个人免费在线观看的高清视频| 亚洲成人国产一区在线观看| e午夜精品久久久久久久| 99在线人妻在线中文字幕| 亚洲精品久久成人aⅴ小说| 嫩草影视91久久| 亚洲乱码一区二区免费版| 搡老岳熟女国产| 国产私拍福利视频在线观看| 舔av片在线| aaaaa片日本免费| netflix在线观看网站| av中文乱码字幕在线| 国内毛片毛片毛片毛片毛片| 免费看日本二区| 欧美性长视频在线观看| 国产免费av片在线观看野外av| av有码第一页| 久久久久性生活片| 特级一级黄色大片| 中文字幕人成人乱码亚洲影| 欧美成人性av电影在线观看| 91麻豆av在线| 亚洲最大成人中文| 美女午夜性视频免费| 国产欧美日韩一区二区三| 人成视频在线观看免费观看| bbb黄色大片| 色尼玛亚洲综合影院| 精品久久久久久久人妻蜜臀av| 高清在线国产一区| 国产蜜桃级精品一区二区三区| 欧美av亚洲av综合av国产av| 丰满的人妻完整版| 99在线人妻在线中文字幕| 亚洲av日韩精品久久久久久密| 免费搜索国产男女视频| 可以在线观看的亚洲视频| 欧美中文综合在线视频| 成年女人毛片免费观看观看9| 淫秽高清视频在线观看| 精品久久久久久久毛片微露脸| 欧美日韩一级在线毛片| 免费在线观看日本一区| 精品少妇一区二区三区视频日本电影| 欧美三级亚洲精品| 日韩欧美三级三区| 青草久久国产| 99热这里只有是精品50| 全区人妻精品视频| 1024香蕉在线观看| 99国产精品一区二区蜜桃av| 狂野欧美白嫩少妇大欣赏| 91麻豆av在线| 国产探花在线观看一区二区| 国产亚洲精品综合一区在线观看 | 午夜激情av网站| 中文亚洲av片在线观看爽| 午夜影院日韩av| 亚洲国产欧洲综合997久久,| 亚洲国产精品成人综合色| 日本一二三区视频观看| av福利片在线| 日韩欧美国产在线观看| 亚洲va日本ⅴa欧美va伊人久久| 黄色 视频免费看| 久久精品aⅴ一区二区三区四区| 亚洲av熟女| 制服人妻中文乱码| 国产高清视频在线播放一区| 国产精品精品国产色婷婷| 亚洲第一欧美日韩一区二区三区| 国产精品免费视频内射| 国产av不卡久久| 国产视频内射| 亚洲一区中文字幕在线| 宅男免费午夜| 久久这里只有精品19| 熟妇人妻久久中文字幕3abv| 免费av毛片视频| av视频在线观看入口| 国产精品久久视频播放| 桃红色精品国产亚洲av| 国内精品久久久久久久电影| 黑人欧美特级aaaaaa片| 白带黄色成豆腐渣| 国产亚洲精品第一综合不卡| 天天添夜夜摸| 99热6这里只有精品| 在线播放国产精品三级| 国产精品乱码一区二三区的特点| 午夜福利高清视频| 在线免费观看的www视频| 90打野战视频偷拍视频| 男女之事视频高清在线观看| 亚洲欧美日韩高清专用| 欧美极品一区二区三区四区| 亚洲自偷自拍图片 自拍| 免费在线观看影片大全网站| 听说在线观看完整版免费高清| 少妇的丰满在线观看| 国产精品免费视频内射| 麻豆国产97在线/欧美 | 日韩中文字幕欧美一区二区| 亚洲无线在线观看| 亚洲男人的天堂狠狠| 欧美成人午夜精品| 淫秽高清视频在线观看| 欧美国产日韩亚洲一区| 嫩草影视91久久| 亚洲黑人精品在线| 国产精品野战在线观看| 国产又色又爽无遮挡免费看| 麻豆国产av国片精品| 好看av亚洲va欧美ⅴa在| 久久精品aⅴ一区二区三区四区| 久久久国产成人免费| 一本综合久久免费| 久久中文看片网| 美女 人体艺术 gogo| 亚洲国产精品合色在线| 亚洲精品色激情综合| 很黄的视频免费| 国产乱人伦免费视频| 哪里可以看免费的av片| 99久久综合精品五月天人人| 日本成人三级电影网站| 一进一出好大好爽视频| 可以在线观看毛片的网站| 欧美日韩瑟瑟在线播放| 午夜福利在线在线| www.999成人在线观看| 久久久久国内视频| 黑人欧美特级aaaaaa片| 一级黄色大片毛片| 国产欧美日韩一区二区三| 国产精品日韩av在线免费观看| 日韩大码丰满熟妇| 亚洲国产中文字幕在线视频| a在线观看视频网站| 很黄的视频免费| 一区二区三区国产精品乱码| 色在线成人网| 免费搜索国产男女视频| 亚洲欧美日韩高清专用| 99re在线观看精品视频| av有码第一页| 久久精品国产99精品国产亚洲性色| 免费观看精品视频网站| 日韩欧美 国产精品| 久久久久久国产a免费观看| 久久久国产成人免费| 国产精品野战在线观看| 欧美+亚洲+日韩+国产| 久久久久久久午夜电影| 在线看三级毛片| 国产伦在线观看视频一区| 天天躁狠狠躁夜夜躁狠狠躁| 亚洲av第一区精品v没综合| 亚洲av电影不卡..在线观看| 在线a可以看的网站| 久久精品国产亚洲av香蕉五月| 国产精品 国内视频| 一卡2卡三卡四卡精品乱码亚洲| 搡老熟女国产l中国老女人| 麻豆久久精品国产亚洲av| 麻豆国产97在线/欧美 | 成人三级做爰电影| www日本在线高清视频| www.999成人在线观看| 成熟少妇高潮喷水视频| 久久欧美精品欧美久久欧美| www.www免费av| 欧美乱色亚洲激情| 国产一区二区在线av高清观看| 欧美极品一区二区三区四区| 亚洲精品在线观看二区| 给我免费播放毛片高清在线观看| 亚洲最大成人中文| or卡值多少钱| 精品国产乱子伦一区二区三区| www日本黄色视频网| 一进一出抽搐动态| 在线a可以看的网站| 日韩欧美免费精品| 日日夜夜操网爽| 丰满的人妻完整版| 精品午夜福利视频在线观看一区| 97碰自拍视频| tocl精华| 精品少妇一区二区三区视频日本电影| 日日干狠狠操夜夜爽| 黄色毛片三级朝国网站| 丰满人妻一区二区三区视频av | 亚洲国产精品成人综合色| 色综合婷婷激情| 免费电影在线观看免费观看| 亚洲av电影不卡..在线观看| 午夜影院日韩av| 午夜精品久久久久久毛片777| 国产精品一区二区精品视频观看| 香蕉丝袜av| 中文字幕久久专区| 亚洲成人国产一区在线观看| 欧洲精品卡2卡3卡4卡5卡区| 国产成人系列免费观看| 亚洲欧美激情综合另类| 99久久久亚洲精品蜜臀av| 久久天躁狠狠躁夜夜2o2o| 女人爽到高潮嗷嗷叫在线视频| 长腿黑丝高跟| 最新美女视频免费是黄的| 国产一区二区在线观看日韩 | 国产伦人伦偷精品视频| 免费一级毛片在线播放高清视频| 99久久精品热视频| 一卡2卡三卡四卡精品乱码亚洲| 人人妻,人人澡人人爽秒播| 成年人黄色毛片网站| 免费无遮挡裸体视频| 国内毛片毛片毛片毛片毛片| 在线国产一区二区在线| 麻豆国产97在线/欧美 | 亚洲国产精品sss在线观看| 一本久久中文字幕| 欧美一级a爱片免费观看看 | 后天国语完整版免费观看| 国产探花在线观看一区二区| 亚洲一区二区三区色噜噜| 19禁男女啪啪无遮挡网站| 国产亚洲av嫩草精品影院| 怎么达到女性高潮| 婷婷精品国产亚洲av在线| 久久精品国产亚洲av高清一级| 高清在线国产一区| 欧美三级亚洲精品| 久久中文字幕人妻熟女| 婷婷丁香在线五月| 精品人妻1区二区| 一级作爱视频免费观看| 三级国产精品欧美在线观看 | 久久午夜综合久久蜜桃| 国产免费av片在线观看野外av| 欧美一级毛片孕妇| 亚洲自拍偷在线| 国产成人av教育| 色哟哟哟哟哟哟| www.www免费av| 国产乱人伦免费视频| 亚洲av第一区精品v没综合| 中文字幕久久专区| 亚洲精品中文字幕在线视频| 国产1区2区3区精品| 动漫黄色视频在线观看| 国产精品精品国产色婷婷| 日韩有码中文字幕| 日韩欧美在线乱码| 又粗又爽又猛毛片免费看| 欧美日韩乱码在线| 伦理电影免费视频| 免费在线观看影片大全网站| 久久中文字幕一级| 国产aⅴ精品一区二区三区波| 男人舔女人的私密视频| 国产欧美日韩一区二区精品| 精品一区二区三区视频在线观看免费| 丝袜人妻中文字幕| 精品不卡国产一区二区三区| 色精品久久人妻99蜜桃| 老司机在亚洲福利影院| 夜夜爽天天搞| 国产精品1区2区在线观看.| 人成视频在线观看免费观看| 精品福利观看| 久久精品91无色码中文字幕| 国语自产精品视频在线第100页| 狠狠狠狠99中文字幕| 成人永久免费在线观看视频| 国产伦人伦偷精品视频| 不卡av一区二区三区| 丝袜美腿诱惑在线| 搡老熟女国产l中国老女人| 国产一区二区三区在线臀色熟女| 久久精品国产综合久久久| 90打野战视频偷拍视频| 亚洲国产日韩欧美精品在线观看 | 少妇人妻一区二区三区视频| 欧美一级毛片孕妇| avwww免费| 国产真人三级小视频在线观看| 日本撒尿小便嘘嘘汇集6| 国内少妇人妻偷人精品xxx网站 | 99精品在免费线老司机午夜| 五月玫瑰六月丁香| 欧美久久黑人一区二区| 十八禁网站免费在线| 两个人看的免费小视频| 欧美黑人欧美精品刺激| 国产99久久九九免费精品| 波多野结衣高清作品| 国内精品久久久久精免费| 一二三四社区在线视频社区8| 男女那种视频在线观看| www.999成人在线观看| 一级毛片女人18水好多| 日韩三级视频一区二区三区| 国产成人精品久久二区二区免费| 精品久久久久久久毛片微露脸| 亚洲国产欧美网| 久久久国产精品麻豆| 日韩欧美在线乱码| 免费观看人在逋| 99国产精品一区二区三区| 国产亚洲精品久久久久5区| 香蕉久久夜色| а√天堂www在线а√下载| 欧美性长视频在线观看| 老鸭窝网址在线观看| av免费在线观看网站| 这个男人来自地球电影免费观看| 亚洲精品在线美女| 麻豆国产av国片精品| 一区二区三区高清视频在线| 亚洲人成伊人成综合网2020| 亚洲乱码一区二区免费版| 少妇的丰满在线观看| www日本黄色视频网| 精品久久久久久成人av| 在线永久观看黄色视频| 99热这里只有是精品50| 久久久久久久午夜电影| 成人永久免费在线观看视频| 人成视频在线观看免费观看| 午夜两性在线视频| 成人av在线播放网站| 午夜福利高清视频| 国产亚洲欧美98| 18禁黄网站禁片免费观看直播| 亚洲欧美精品综合久久99| 亚洲片人在线观看| 欧美久久黑人一区二区| 久久久久久国产a免费观看| 亚洲成av人片在线播放无| 国产亚洲精品久久久久5区| 欧美日韩瑟瑟在线播放| 日韩大码丰满熟妇| 久久久久久人人人人人| 国产亚洲欧美98| www日本黄色视频网| 久久伊人香网站| 香蕉av资源在线| 九色国产91popny在线| 亚洲熟妇中文字幕五十中出| 黑人操中国人逼视频| 日本免费a在线| 白带黄色成豆腐渣| 亚洲五月天丁香| 99国产精品一区二区三区| 亚洲中文av在线| 香蕉av资源在线| 久久久久九九精品影院| 一本大道久久a久久精品| 99久久精品热视频| 最新在线观看一区二区三区| 亚洲精品中文字幕在线视频| 欧美成人一区二区免费高清观看 | 深夜精品福利| 亚洲欧美激情综合另类| 国产精品国产高清国产av| 色综合亚洲欧美另类图片| 久久久久久亚洲精品国产蜜桃av| 成人特级黄色片久久久久久久| 免费人成视频x8x8入口观看| 欧美 亚洲 国产 日韩一| 欧美黑人巨大hd| 精品一区二区三区视频在线观看免费| 非洲黑人性xxxx精品又粗又长| 亚洲国产精品999在线| 日本三级黄在线观看| 深夜精品福利| 韩国av一区二区三区四区| 51午夜福利影视在线观看| 亚洲专区国产一区二区| 欧美日韩福利视频一区二区| 国产三级在线视频| svipshipincom国产片| 精品乱码久久久久久99久播| 日日摸夜夜添夜夜添小说| 俄罗斯特黄特色一大片| 精品福利观看| 别揉我奶头~嗯~啊~动态视频| 99国产极品粉嫩在线观看| 中出人妻视频一区二区| 99热只有精品国产| 日韩欧美在线乱码| 免费搜索国产男女视频| 悠悠久久av| 亚洲免费av在线视频| 婷婷六月久久综合丁香| 啦啦啦韩国在线观看视频| 最近视频中文字幕2019在线8| 色老头精品视频在线观看| 嫁个100分男人电影在线观看| 18禁观看日本| 亚洲av电影不卡..在线观看| 亚洲av美国av| 搡老熟女国产l中国老女人| 夜夜爽天天搞| 老熟妇仑乱视频hdxx| 又爽又黄无遮挡网站| 两个人免费观看高清视频| 中文字幕久久专区| 搡老妇女老女人老熟妇| 亚洲电影在线观看av| 国内精品久久久久久久电影| 久久亚洲精品不卡| 欧美最黄视频在线播放免费| 黄频高清免费视频| 动漫黄色视频在线观看| 亚洲免费av在线视频| 少妇被粗大的猛进出69影院| 99精品在免费线老司机午夜| 男女之事视频高清在线观看| 一区二区三区激情视频| 国产亚洲欧美98| 日韩欧美一区二区三区在线观看| 亚洲 欧美一区二区三区| 免费观看人在逋| 亚洲熟女毛片儿| 很黄的视频免费| 婷婷六月久久综合丁香| 真人做人爱边吃奶动态| 一个人观看的视频www高清免费观看 | 首页视频小说图片口味搜索| 男女做爰动态图高潮gif福利片| 99精品在免费线老司机午夜| 亚洲精品中文字幕一二三四区| 一区二区三区激情视频| 欧美黄色淫秽网站| 日韩欧美三级三区| 久久精品亚洲精品国产色婷小说| 麻豆av在线久日| 男人舔奶头视频| √禁漫天堂资源中文www| 国产黄a三级三级三级人| 国产亚洲精品综合一区在线观看 | 欧美一区二区精品小视频在线| 精品欧美一区二区三区在线| 深夜精品福利| 欧美性长视频在线观看| 亚洲色图 男人天堂 中文字幕| 亚洲美女黄片视频| 国产在线精品亚洲第一网站| 国产亚洲精品久久久久5区| 成人特级黄色片久久久久久久| 欧美色欧美亚洲另类二区| 香蕉av资源在线| 黑人操中国人逼视频| 在线观看午夜福利视频| 久久欧美精品欧美久久欧美| 亚洲欧洲精品一区二区精品久久久| av有码第一页| 一二三四在线观看免费中文在| 无人区码免费观看不卡| 久久精品国产亚洲av香蕉五月| 精品久久久久久久人妻蜜臀av| 母亲3免费完整高清在线观看| 18禁观看日本| 国产精品永久免费网站| or卡值多少钱| 全区人妻精品视频| 国产成人欧美在线观看| 天堂动漫精品| 亚洲一区中文字幕在线| 亚洲人与动物交配视频| 少妇被粗大的猛进出69影院| 黑人操中国人逼视频| 俺也久久电影网| 无遮挡黄片免费观看| 淫妇啪啪啪对白视频| 免费观看人在逋| 色哟哟哟哟哟哟| 成人午夜高清在线视频| 人妻夜夜爽99麻豆av| 成人亚洲精品av一区二区| 女生性感内裤真人,穿戴方法视频| 亚洲色图av天堂| 国产精品美女特级片免费视频播放器 | 久久99热这里只有精品18| 亚洲午夜理论影院| 国产亚洲av嫩草精品影院| 18禁黄网站禁片免费观看直播| 欧美最黄视频在线播放免费| 母亲3免费完整高清在线观看| 男人舔奶头视频| 成人18禁高潮啪啪吃奶动态图| 一个人免费在线观看电影 | 久久香蕉国产精品| 变态另类成人亚洲欧美熟女| 18禁黄网站禁片午夜丰满| 久久国产乱子伦精品免费另类| 精品日产1卡2卡| 亚洲精品美女久久久久99蜜臀| 婷婷亚洲欧美| 亚洲一卡2卡3卡4卡5卡精品中文| 亚洲精品一区av在线观看| 欧美成人免费av一区二区三区| 亚洲精品美女久久av网站| 免费无遮挡裸体视频| 亚洲欧美激情综合另类| 国产成年人精品一区二区| 日韩精品中文字幕看吧| 国产真实乱freesex| 国产av一区二区精品久久| 91在线观看av| 亚洲欧美日韩高清专用| 久久这里只有精品中国| 搡老熟女国产l中国老女人| 香蕉国产在线看| 国产真人三级小视频在线观看| 真人做人爱边吃奶动态| 一本久久中文字幕| 中文字幕人成人乱码亚洲影| 色综合婷婷激情| 成人手机av| 在线观看一区二区三区| 91在线观看av| 麻豆国产av国片精品| 国产蜜桃级精品一区二区三区| 国产伦一二天堂av在线观看| 成人av在线播放网站| 香蕉av资源在线| 成人高潮视频无遮挡免费网站| 无人区码免费观看不卡| 88av欧美| 777久久人妻少妇嫩草av网站| 97超级碰碰碰精品色视频在线观看| 国产欧美日韩精品亚洲av| 一二三四在线观看免费中文在| 国产成人欧美在线观看| 亚洲国产欧美人成| 日本黄色视频三级网站网址| 欧美另类亚洲清纯唯美| 亚洲性夜色夜夜综合| 成人av一区二区三区在线看| 日韩欧美国产在线观看|