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

    基于eDNA宏條形碼技術(shù)的巢湖夏季魚類群落結(jié)構(gòu)及多樣性

    2025-08-12 00:00:00劉成鳳嚴(yán)云志韓勝盼張家赫謝平陳雋夏午來
    水生態(tài)學(xué)雜志 2025年4期
    關(guān)鍵詞:巢湖魚類群落

    中圖分類號(hào):S932.4 文獻(xiàn)標(biāo)志碼:A 文章編號(hào):1674-3075(2025)04-0229-09

    淡水生態(tài)系統(tǒng)在小范圍內(nèi)提供了高度的棲息地多樣性(Higginsetal,2021),它雖然占地球表面積不到 1% ,但卻承載著地球上約 12% 已知物種的繁衍生息(Dudgeon,2019)。魚類是水生生態(tài)系統(tǒng)中最多樣化的脊椎動(dòng)物群,支持著關(guān)鍵的生態(tài)系統(tǒng)服務(wù)。然而,它們的種群越來越容易受到由人類直接或間接造成的影響,如棲息地破壞、過度捕撈、富營(yíng)養(yǎng)化和污染等(Worm etal,2006;EmmettDuffy et al,2016)。魚類多樣性是生物多樣性的重要組成部分,研究魚類多樣性可以有效監(jiān)測(cè)水生生態(tài)系統(tǒng)的健康狀況,對(duì)漁業(yè)資源保護(hù)、恢復(fù)和可持續(xù)發(fā)展具有重要意義(Phangetal,2019;Zou etal,2020)。

    巢湖是中國(guó)5大淡水湖之一,位于長(zhǎng)江中下游,地處安徽省中部,承擔(dān)著供水、防洪、航運(yùn)、灌溉等多種功能,是長(zhǎng)江經(jīng)濟(jì)帶和長(zhǎng)三角經(jīng)濟(jì)圈的重要板塊。(Yuetal,2011;張方方,2014;Fangetal,2017)。近年來,巢湖周邊經(jīng)濟(jì)的迅猛發(fā)展,導(dǎo)致巢湖流域生態(tài)、湖區(qū)水域環(huán)境日益惡化以及水生生物種群結(jié)構(gòu)發(fā)生重大改變,包括魚類種類急劇減少、水土流失日益加劇、湖泊調(diào)蓄功能削弱等,這些問題的發(fā)展和蔓延導(dǎo)致巢湖生物群落結(jié)構(gòu)趨向單一化,生態(tài)系統(tǒng)穩(wěn)定性降低,嚴(yán)重影響了全流域經(jīng)濟(jì)發(fā)展和環(huán)境保護(hù)(Kong etal,2013;Yang etal,2013;Yaoamp;He,2020)。

    環(huán)境DNA(environmentDNA,簡(jiǎn)稱eDNA)是指從環(huán)境(土壤、水、空氣等)而非單個(gè)有機(jī)體中收集的DNA片段的總和(Reesetal,2014)。而eDNA宏條形碼技術(shù)則是富集空氣、水體、土壤中等各種生物在活動(dòng)場(chǎng)所中殘留于環(huán)境中的DNA,通過PCR擴(kuò)增、高通量測(cè)序等一系列操作,實(shí)現(xiàn)對(duì)物種的識(shí)別鑒定(Ruppertetal,2019)。近年來,eDNA宏條形碼技術(shù)在水生生態(tài)系統(tǒng)中的應(yīng)用涉及多種水體環(huán)境,包括池塘(Harperetal,2018)、湖泊(Tkachukamp;Dunn,2020;高婉茹,2022)、河流(Pontetal,2018)、海洋(王雪華,2021)。與傳統(tǒng)漁業(yè)資源調(diào)查(比如電魚、拖網(wǎng)、圍網(wǎng))相比,eDNA宏條形碼技術(shù)整個(gè)過程不需要采集目標(biāo)生物,以采樣方法簡(jiǎn)單、檢測(cè)靈敏度高和對(duì)生物不具破壞性的優(yōu)勢(shì)彌補(bǔ)了傳統(tǒng)物種監(jiān)測(cè)的不足,在評(píng)估生物多樣性、確定物種分布、種群動(dòng)態(tài)、生態(tài)系統(tǒng)健康等方面具有較強(qiáng)的應(yīng)用潛能(Thomsenamp;Willerslev,2015)。本研究嘗試引入新興的eDNA宏條形碼技術(shù)來探究巢湖夏季魚類群落結(jié)構(gòu)、多樣性及其與環(huán)境因子的關(guān)系,以期為巢湖水生生物監(jiān)測(cè)以及生態(tài)保護(hù)提供新的技術(shù)手段,為今后長(zhǎng)江流域開展?jié)O業(yè)資源調(diào)查提供新思路。

    1材料與方法

    1.1樣點(diǎn)設(shè)置及eDNA樣品采集

    2022年8月在巢湖設(shè)置15個(gè)采樣點(diǎn)進(jìn)行水樣采集,從左至右分為3個(gè)湖區(qū),分別為西湖區(qū)(CH1\~CH5)、中湖區(qū)(CH6\~CH10)、東湖區(qū)(CH11\~CH15),具體采樣點(diǎn)信息見圖1。水樣采集使用采水器采集表層水和底層水各 10L ,混勻后取1L水樣保存至已消毒的1L廣口瓶中。采集后在 24h 內(nèi),使用 0.7μm 玻璃纖維濾膜(Whatman公司)進(jìn)行過濾,1個(gè)采樣點(diǎn)共過濾3張濾膜,每張濾膜過濾水樣300mL ,同時(shí)過濾同等體積的純凈水作為陰性對(duì)照,將過濾的濾膜放入 2mL 滅菌離心管中,置于 -40°C 冰箱,運(yùn)回實(shí)驗(yàn)室后立即放入, -80°C 的超低溫冰箱。為盡可能減少DNA污染,采樣及過濾器具均使用 10% 的次氯酸消毒,操作人員在采樣及過濾過程中全程佩戴一次性手套,并在不同采樣點(diǎn)之間更換手套。每個(gè)濾膜DNA的提取采用OMEGA組織提取試劑盒按照說明方法進(jìn)行,并將提取到的eDNA樣品迅速置于 -80°C 保存?zhèn)溆谩?/p>

    圖1巢湖采樣點(diǎn)位置分布Fig.1 LocationofthesamplingsitesinChaohu Lake

    1.2PCR擴(kuò)增與高通量測(cè)序

    本研究使用魚類通用引物Tele02F/R進(jìn)行PCR擴(kuò)增,擴(kuò)增引物Tele02F: 5 -AAACTC-

    GTGCCAGCCACC- ?3 ;Tele02_R: 5 -GGGTATCTA-ATCCCAGTTTG- ?3 。擴(kuò)增片段長(zhǎng)度為 151~200bp 。PCR擴(kuò)增體系 20μL:4μL 5× FastPfuBuffer, 2μL dNTPs(2.5mmol/L) , 0.8μL 上下游引物 (5μmol/L) ,0.4μL FastPfu Polymerase,模 板DNA ,用 補(bǔ)至 20μL 。擴(kuò)增條件: 95°C 預(yù)變性 5min :95°C 變性 30s.55°C 退火 30s,72‰ 延伸 45s(35 個(gè)循環(huán));最后 72% 延伸 10min 。每個(gè)樣本進(jìn)行3次重復(fù)PCR擴(kuò)增,將同一樣本的PCR產(chǎn)物混合后用 2% 瓊脂糖凝膠電泳檢測(cè),使用AxyPrepDNA凝膠提取試劑盒(AXYGEN公司)進(jìn)行純化。陰性對(duì)照未擴(kuò)增出目的條帶,將檢測(cè)合格的環(huán)境DNA樣本送至上海凌恩生物科技有限公司進(jìn)行文庫(kù)構(gòu)建,用IlluminaNova-SeqPE250平臺(tái)進(jìn)行配對(duì)端測(cè)序 (2×250bp 。

    原始測(cè)序數(shù)據(jù)以FASTQ格式保存,根據(jù)barcode得到所有樣品的有效序列,然后對(duì)reads的質(zhì)量進(jìn)行質(zhì)控過濾(過濾reads尾部質(zhì)量值20以下的堿基,丟棄質(zhì)控后 50bp 以下的reads,去除含有模糊字符的、不能被組裝的reads)。使用UPARSE軟件,對(duì)優(yōu)化序列提取非重復(fù)序列,去除沒有重復(fù)的單序列;按照97% 相似性對(duì)非重復(fù)序列(不含單序列)進(jìn)行OTU聚類,在聚類過程中采用denovo和reference結(jié)合的方式去除嵌合體,得到OTU的代表序列;將所有優(yōu)化序列map至OTU代表序列,選出與OTU代表序列相似性在 97% 以上的序列,生成OTU表格。采用Uclust算法對(duì) 97% 相似水平的OTU代表序列進(jìn)行分類學(xué)分析,并分別在各個(gè)分類水平:域、界、門、綱、目、科、屬、種統(tǒng)計(jì)各樣本的群落組成。代表性O(shè)TU序列經(jīng)魚類線粒體數(shù)據(jù)庫(kù)和凌恩生物構(gòu)建的淡水魚類數(shù)據(jù)庫(kù)進(jìn)行物種注釋分析,再結(jié)合巢湖魚類生物自建庫(kù)對(duì)所注釋物種進(jìn)行人工校對(duì)并手動(dòng)去除非魚類信息。

    1.3環(huán)境因子的測(cè)定

    在采樣過程中使用便攜式水深儀現(xiàn)場(chǎng)測(cè)定水深(Depth)、多參數(shù)水質(zhì)測(cè)定儀YSI6600測(cè)定pH、水溫(WT)、鹽度(Sal)、溶解氧(DO)和濁度(Turb)等。其余理化指標(biāo)回到實(shí)驗(yàn)室后測(cè)定,其中測(cè)定溶解性總氮(DTN)使用堿性過硫酸鉀消解分光光度法;氨氮(NH4+–N) 使用納氏試劑分光光度法測(cè)定;化學(xué)需氧量(COD)使用高錳酸鉀滴定法;總有機(jī)碳(TOC)使用總有機(jī)碳分析儀(王晨等,2022)。

    1.4數(shù)據(jù)分析

    本研究基于各采樣點(diǎn)的魚類物種序列豐度進(jìn)行數(shù)據(jù)分析,所有數(shù)據(jù)預(yù)處理在Excel2019中完成,統(tǒng)計(jì)分析及作圖在R(version4.2.2)中完成,除熱圖需使用pheatmap包外,其余所有圖片均使用R的gg-plot2包制作。使用R中的vegan包計(jì)算Alpha多樣性指數(shù)(Shannon-Wiener指數(shù)、Simpson指數(shù)、Pielou均勻度指數(shù)),使用vegan包中的vegdist函數(shù)計(jì)算物種的Bray-Curtis距離;采用t檢驗(yàn)或Wilcoxon秩和檢驗(yàn)進(jìn)行Alpha多樣性組間差異分析,并用ggsigif包進(jìn)行顯著性標(biāo)記;使用R中的vegan包基于Bray-Curtis距離進(jìn)行主坐標(biāo)分析(PCoA),并結(jié)合置換多元方差分析(PERMANOVA)來檢驗(yàn)不同湖區(qū)魚類群落結(jié)構(gòu)差異是否顯著;使用vegan包進(jìn)行RDA分析,解析魚類群落結(jié)構(gòu)與環(huán)境因子的關(guān)系,并運(yùn)用蒙特卡洛置換檢驗(yàn)(MonteCarlo)分析環(huán)境因子對(duì)魚類群落影響的顯著程度。采用Origin2023繪制巢湖優(yōu)勢(shì)魚類物種組成圖。Alpha多樣性指數(shù)(Shan-non-Wiener指數(shù)、Simpson指數(shù)、Pielou均勻度指數(shù))計(jì)算公式如下:

    式中: H 為Shannon-Wiener指數(shù); Pi 為物種 i 的個(gè)體數(shù)占總個(gè)體數(shù)的比例; D 為Simpson指數(shù); J 為Pielou均勻度指數(shù);S為物種總數(shù)。

    2結(jié)果與分析

    2.1魚類物種組成

    巢湖eDNA樣品高通量測(cè)序后獲得高質(zhì)量序列9504070條,按照 97% 相似性聚類后得到7040個(gè)OTUs,平均序列長(zhǎng)度為 170.89bp ,在數(shù)據(jù)庫(kù)中匹配到魚類48種,隸屬于7目13科39屬(表1),包含5種在巢湖從未報(bào)道過的魚類物種,其中15個(gè)樣本中共有3609886reads序列數(shù)無法注釋到種的分類學(xué)水平,占比 37% 。在目分類水平中,鯉形目種類最多,共計(jì)35種,約占檢測(cè)總數(shù)的 72.92% ;其次是鱸形目6種,占 12.5% ;鲇形目和頜針魚目各2種,約占 4.17% ;而鮭形目、鯡形目和鰻目均僅有1種,約占 2.08% 。

    2.2巢湖魚類多樣性

    巢湖不同湖區(qū)魚類 a 多樣性指數(shù)如表2所示。Shannon-Wiener指數(shù)為 0.743~2.281 ,平均值為1.710±0.375 ;Simpson指數(shù)為 0.348~0.860 ,平均值為0.696±0.140 ;Pielou均勻度指數(shù)為 0.221~0.677 ,平均值為 0.504±0.112 。Shannon-Wiener、Simpson 和Pielou均勻度指數(shù)最高值均出現(xiàn)在CH15采樣點(diǎn),表明該采樣點(diǎn)的魚類物種多樣性較高。巢湖不同湖區(qū)魚類的3種多樣性指數(shù)分布趨勢(shì)基本相同,且差異均不顯著(圖2)。

    表1基于eDNA宏條形碼技術(shù)在巢湖檢測(cè)到的魚類物種Tab.1 Fish species detected in Chaohu Lake based on eDNA metabarcoding technolog:
    表2巢湖各采樣點(diǎn)魚類 ∝ 多樣性指數(shù)

    基于各采樣點(diǎn)魚類種水平的分布熱圖(圖3)可知不同采樣點(diǎn)間魚類物種組成存在差異,鳙(Aris-tichthysnobilis)和刀(Coilianasus)的相對(duì)序列豐度和分布頻率最高。由圖4可知相對(duì)豐度前10位的優(yōu)勢(shì)魚類物種分別為、刀、翹嘴(Culteralbur-nus)、鯉(Cyprinuscarpio)、鯽(Carassiusauratus)、鰱(Hypophthalmichthysmolitrix)、達(dá)氏(Chanodich-thysdabryi)、張氏(Hemicultertchangi)、團(tuán)頭魴(Megalobramaamblycephala)、似(Toxabramisswinhonis)。

    對(duì)巢湖不同湖區(qū)魚類群落組成進(jìn)行PCoA分析(圖5),結(jié)果顯示,2個(gè)主坐標(biāo)對(duì)魚類群落組成差異的解釋度為 62.35% 。東湖區(qū)樣本點(diǎn)之間距離相對(duì)較遠(yuǎn),表明東湖區(qū)各樣本的魚類群落組成差異較大;西湖區(qū)樣本點(diǎn)之間距離較近,且有聚集趨勢(shì),表明西湖區(qū)各樣本的魚類群落組成相似。PERMANOVA檢驗(yàn)表明,巢湖不同湖區(qū)魚類群落結(jié)構(gòu)差異不顯著 (R2=0.13 ,P=0.53gt;0.05 。

    2.3巢湖魚類群落與環(huán)境因子關(guān)系

    對(duì)魚類物種數(shù)據(jù)“Hellinger\"轉(zhuǎn)換后進(jìn)行去趨勢(shì)對(duì)應(yīng)分析(DCA),根據(jù)分析結(jié)果第一排序軸 lt;3 ,故選用RDA模型。在冗余分析(RDA)前對(duì)18個(gè)環(huán)境因子進(jìn)行預(yù)選,刪除共線性強(qiáng)的環(huán)境因子,最后保留的10個(gè)環(huán)境因子(Depth、WT、Turb、COD、 NH4+ -N、DO、pH、DTN、TOC、Sal)VIF值均小于10。如圖6所示,前2軸累計(jì)解釋變化率為 51.97% 。在RDA1軸上,NH4+ -N和DTN呈負(fù)相關(guān),其余因子呈正相關(guān)。在RDA2軸上,Sal、Depth、DTN呈正相關(guān),而其余因子呈負(fù)相關(guān)。蒙特卡洛置換檢驗(yàn)分析結(jié)果表明 NH4+ -N(F=4.2436,P=0.001) 對(duì)魚類群落結(jié)構(gòu)具有極顯著影響(Plt;0.01) ,WT(F=2.0486,P=0.039) 和D )(F=1.9169 .P=0.041) 對(duì)魚類群落結(jié)構(gòu)具有顯著影響 (Plt;0.05) ,其余環(huán)境因子對(duì)其影響不顯著 (Pgt;0.05 )。

    3討論

    3.1基于eDNA宏條形碼技術(shù)的巢湖魚類物種組成及多樣性

    對(duì)巢湖漁業(yè)資源的調(diào)查最早始于20世紀(jì)50年代末60年代初,此后陸續(xù)有學(xué)者對(duì)其進(jìn)行研究。1987年調(diào)查記錄了巢湖湖區(qū)魚類85種(王岐山,1987),2005年巢湖湖區(qū)魚類54種(過龍根,2005),2016年采集到巢湖流域河流魚類47種(錢紅等,2016)。以上研究都是通過傳統(tǒng)漁業(yè)資源調(diào)查方法完成的,不僅耗時(shí)費(fèi)力,還會(huì)影響漁業(yè)資源發(fā)展,破壞水生生態(tài)系統(tǒng)的穩(wěn)定性(李曉玲等,2022)。本研究利用eDNA宏條形碼技術(shù)調(diào)查分析了巢湖夏季魚類物種組成和多樣性,共檢測(cè)到48種魚類,以鯉形目魚類為主,這與梁陽陽等(2022)通過傳統(tǒng)漁業(yè)資源調(diào)查方法監(jiān)測(cè)到的結(jié)果較為吻合。研究結(jié)果表明,通過eDNA宏條形碼技術(shù)檢出的魚類物種與傳統(tǒng)調(diào)查方法監(jiān)測(cè)到的物種在很大程度上吻合,并檢測(cè)到2種在巢湖歷史記錄中存在,但是近年文獻(xiàn)報(bào)道中傳統(tǒng)調(diào)查方法均未捕獲到的物種,分別為赤眼(Squalio-barbuscurriculus)和唇(Hemibarbuslabeo)。說明eDNA宏條形碼技術(shù)可以有效實(shí)現(xiàn)魚類監(jiān)測(cè),作為未來探究巢湖魚類組成的輔助手段。此外,eDNA宏條形碼技術(shù)還檢測(cè)出了5種未報(bào)道過的魚類,分別為張氏、汪氏近紅(Ancherythroculterwangi)、斜方(Acheilognathusrhombeus)、斑鱖(Siniperca scher-zeri)、大口黑鱸(Micropterussalmoides)。這5種魚在長(zhǎng)江水系均有分布,在巢湖檢出可能是由于含有其DNA的長(zhǎng)江水體通過裕溪河流至巢湖,因?yàn)榄h(huán)境DNA可以在水體中存在幾天到幾周的時(shí)間;或者是采樣點(diǎn)附近存在住宅,居民們將含有魚類DNA的廢水直接排放到湖內(nèi);也有可能是部分采樣地點(diǎn)位于航道附近,往返漁船排出的水中存在上述魚類DNA,在采樣過程中被采集檢出;又或者是由于物種注釋有誤。這些因素都會(huì)造成對(duì)巢湖真實(shí)魚類多樣性的誤判,但其是否在巢湖真實(shí)存在,有待進(jìn)一步的調(diào)查研究。歷史記錄存在的黃鱔(Monopterusalbus)、中華螃(Rhodeussinensis)、圓尾斗魚(Macropodusocellatus)、烏(Ophiocephalusargus)等未檢測(cè)到,這可能與魚類生活習(xí)性不同有關(guān),如黃鱔(Monopterusalbus)等穴居種類不易被檢測(cè)到,也有可能是環(huán)境DNA在水體中的降解;比對(duì)數(shù)據(jù)庫(kù)不完善;采樣、提取和儲(chǔ)存過程中受到污染及通用引物存在局限性造成的。由此可見eDNA宏條形碼技術(shù)自身也存在缺點(diǎn),不能完全代替?zhèn)鹘y(tǒng)調(diào)查方法,但可以作為一種補(bǔ)充手段,減少傳統(tǒng)監(jiān)測(cè)對(duì)巢湖生態(tài)系統(tǒng)造成的不良影響,將二者結(jié)合可更完整地反映巢湖魚類多樣性及其空間分布(王汝賢等,2023)。

    圖2巢湖不同湖區(qū)的 ∝ 多樣性指數(shù)Fig.2FishαdiversityindicesforeachareasampledinChaohuLake
    圖3巢湖各采樣點(diǎn)種水平上的魚類分布
    圖4巢湖各采樣點(diǎn)優(yōu)勢(shì)魚類物種組成 Fig.4Dominantfishspeciesateach samplingpointinChaohuLake
    圖5巢湖魚類群落結(jié)構(gòu)的主坐標(biāo)分析Fig.5 Principal coordinateanalysisoffishcommunitystructureinChaohuLake
    圖6巢湖魚類群落結(jié)構(gòu)與環(huán)境因子的RDA分析* denotes Plt;0.05 *** denotes Plt;0.001 元Fig.6RedundancyanalysisoffishcommunitystructureandenvironmentalfactorsinChaohuLake

    魚類群落結(jié)構(gòu)組成差異性正在成為生態(tài)系統(tǒng)功能的關(guān)鍵驅(qū)動(dòng)因素,為生物保護(hù)提供信息(Lamyetal,2021)。本研究中, a 多樣性分析結(jié)果顯示巢湖不同湖區(qū)魚類群落結(jié)構(gòu)不存在顯著差異。 a 多樣性指數(shù)(Shannon-Wiener指數(shù)、Simpson指數(shù)、Pielou均勻度指數(shù))最高值均出現(xiàn)在CH15采樣點(diǎn),這個(gè)采樣點(diǎn)位于巢湖與柘皋河交匯處,且與裕溪河距離較近,河湖水系連通有利于河流與湖泊之間生物的交換,提高水資源優(yōu)化配置、改善河湖健康狀態(tài),這可能是該點(diǎn)位魚類多樣性較高的原因(毛志剛等,2019)。PCoA分析表明不同采樣點(diǎn)魚類群落結(jié)構(gòu)呈現(xiàn)一定的差異性,巢湖東湖區(qū)的魚類群落物種組成差異較大;西湖區(qū)群落物種組成相似;中湖區(qū)與東湖區(qū)之間存在明顯交叉,說明群落物種組成與東湖區(qū)相似性更高。PERMANOVA檢驗(yàn)結(jié)果顯示巢湖不同湖區(qū)魚類群落結(jié)構(gòu)差異不顯著 (Pgt;0.05) ,這與 a 多樣性分析結(jié)果一致。

    3.2環(huán)境因子對(duì)巢湖魚類群落結(jié)構(gòu)的影響

    魚類群落結(jié)構(gòu)受生物因素和非生物因素的雙重影響,與環(huán)境因子之間的關(guān)系在不同規(guī)模和生態(tài)系統(tǒng)類型之間差異顯著(Rotherhametal,2011;Larsonetal,2013)。相關(guān)研究表明,水溫、水深、總磷、總氮、pH、溶解氧等環(huán)境因子對(duì)湖泊魚類群落分布均存在一定影響(程琳,2011;劉燕山等,2021)。劉鵬飛等(2022)在固城湖魚類群落結(jié)構(gòu)的研究中表明總氮、高錳酸鹽指數(shù)和濁度對(duì)魚類群落結(jié)構(gòu)的時(shí)空差異有顯著性影響。Zhang等(2019)利用eDNA技術(shù)調(diào)查研究長(zhǎng)江口及其鄰近水域的魚類多樣性,結(jié)果顯示水溫、鹽度和溶氧是影響季節(jié)間魚類組成差異的主要環(huán)境因子。王銀肖等(2022)在研究中表明影響白洋淀魚類群落的主要環(huán)境因子是水溫、DO和Chl-a。本研究通過RDA分析得出,氨氮 (NH4+–N) 對(duì)巢湖魚類群落結(jié)構(gòu)具有極顯著影響 (Plt;0.01) ,水溫(WT)和溶解氧(DO)對(duì)巢湖魚類群落結(jié)構(gòu)具有顯著影響 (Plt;0.05) ,其余環(huán)境因子(DTN、TOC、pH、Depth、Sal、Turb、COD)影響不顯著 (Pgt;0.05) 。氨氮是水環(huán)境中一個(gè)重要的污染指標(biāo),是藻類必需的一種營(yíng)養(yǎng)元素(李月萍,1998)。水體中氮含量發(fā)生變化,會(huì)導(dǎo)致水質(zhì)發(fā)生改變,能直接或間接影響魚類生長(zhǎng)、繁殖。水溫會(huì)直接影響魚類的新陳代謝強(qiáng)度,從而影響到魚類生活的各個(gè)方面;對(duì)魚類性腺發(fā)育也存在一定影響,夏季水溫升高,性成熟年齡提前,會(huì)減少懷卵和排卵量。溶解氧可以促進(jìn)水體中物質(zhì)的良性循環(huán),反映棲息地的適宜性(宋超等,2022)。巢湖是典型的富營(yíng)養(yǎng)化湖泊,夏季氣溫升高,水體中藍(lán)藻數(shù)量也會(huì)增多,會(huì)導(dǎo)致溶解氧含量減少,水環(huán)境中溶解的差異會(huì)影響浮游動(dòng)物群落的大小和分類,對(duì)以浮游動(dòng)物為主要食物來源的魚類,如鰱、鳙、幼齡刀等影響較大。本研究中,巢湖以鳙、刀的OTUs序列豐度相較于其他魚類更高,這可能也是溶解氧成為影響魚類群落結(jié)構(gòu)主要環(huán)境因子的原因。

    參考文獻(xiàn)

    程琳,2011.長(zhǎng)江流域湖泊魚類群落與主要環(huán)境因子關(guān)系及 其預(yù)測(cè)模型[D].北京:中國(guó)科學(xué)院研究生院

    高婉茹,2022.基于環(huán)境DNA的湖泊型水體魚類多樣性研究 [D].上海:上海海洋大學(xué).

    過龍根,2005.大型淺水富營(yíng)養(yǎng)湖泊一巢湖的漁業(yè)生態(tài)學(xué)研 究[D].武漢:中國(guó)科學(xué)院研究生院(水生生物研究所).

    李曉玲,劉洋,王叢叢,等,2022.基于環(huán)境DNA技術(shù)的夏季 東海魚類物種多樣性研究[J].海洋學(xué)報(bào),44(4):74-84.

    LI X L,LIU Y, WANG C C, et al, 2022. Study on fish species diversity in the East China Sea in summer based on environmental DNA technology[J]. Haiyang Xuebao, 44(4): 74-84.

    梁陽陽,盧文軒,楊坤,等,2022.巢湖魚類群落結(jié)構(gòu)及多樣性 [J].上海海洋大學(xué)學(xué)報(bào),31(6):1445-1456.

    LIANG Y Y,LU W X, YANG K, et al, 2022. The community structure and diversity of fish in the ChaohuLake[J]. Journal of Shanghai Ocean University, 31(6):1445-1456.

    劉燕山,谷先坤,唐晟凱,等,2021.東太湖魚類群落結(jié)構(gòu)特征 及其與環(huán)境因子的關(guān)系[J].生態(tài)學(xué)報(bào),41(2):769-780.

    LIU Y S,GU X K, TANG S K, et al, 2021. Fish community structure and its relationship with main environmental variablesinthe East Taihu Lake[J].Acta Ecologica Sinica, 41(2):769-780.

    劉鵬飛,景麗,任瀧,等,2022.固城湖魚類群落結(jié)構(gòu)現(xiàn)狀及其 與環(huán)境因子的關(guān)系[J].大連海洋大學(xué)學(xué)報(bào),37(5):841- 849.

    LIU PF, JING L,REN L, et al, 2022. Current fish community structure and its relationship with environmental factors in Gucheng Lake,Nanjing China[J]. Journal of Dalian Ocean University,37(5):841-849.

    李月萍,1998.氨氮對(duì)魚類的影響[J].內(nèi)陸水產(chǎn),23(4):26.

    毛志剛,谷孝鴻,龔志軍,等,2019.洪澤湖魚類群落結(jié)構(gòu)及其 資源變化[J].湖泊科學(xué),31(4):1109-1119.

    MAO ZG, GU X H, GONG ZJ, et al, 2019. The structure of fish community and changes of fishery resources in Lake Hongze[J]. Journal of Lake Sciences,31(4):1109-1119.

    錢紅,嚴(yán)云志,儲(chǔ)玲,等,2016.巢湖流域河流魚類群落的時(shí)空 分布[J].長(zhǎng)江流域資源與環(huán)境,25(2):257-264.

    QIAN H, YAN Y Z, CHUL, et al, 2016. Spatial and temporal patterns of fish assemblages in the rivers of Chaohu basin [J].Resources and Environment in the Yangtze Basin, 25 (2):257-264.

    宋超,胡麗娟,趙峰,等,2022.長(zhǎng)江口海上風(fēng)電場(chǎng)水域魚類群 落結(jié)構(gòu)特征及其與環(huán)境因子的關(guān)系[J].中國(guó)水產(chǎn)科學(xué), 29(3):469-482.

    SONG C,HU L J, ZHAO F, et al, 2022. Fish community structure and its relationship with environmental factors in offshore wind farm waters of the Yangtze Estuary[J]. Journal of Fishery Sciences of China, 29(3):469-482.

    王雪華,2021.基于燈光敷網(wǎng)和環(huán)境DNA宏條形碼的西北 太平洋游泳動(dòng)物多樣性比較研究[D].上海:上海海洋 大學(xué).

    王晨,陶孟,李愛民,等,2022.基于環(huán)境DNA宏條形碼技術(shù) 的秦淮河生物多樣性研究[J].生態(tài)學(xué)報(bào),42(2):611-624.

    WANG C, TAO M,LI A M, et al,2022. Research on the biodiversity of Qinhuai River based on environmental DNA metabacroding[J]. Acta Ecologica Sinica, 42(2):611-624.

    王岐山,1987.巢湖魚類區(qū)系研究[J].安徽大學(xué)學(xué)報(bào)(自然科 學(xué)版),11(2):70-78.

    WANG Q S,1987. Study on ICHTHYOLOGlCAL fauna of Chao lake[J]. Journal of Anhui University (Natural Sciences), 11(2):70-78.

    王汝賢,楊剛,耿智,等,2023.環(huán)境DNA技術(shù)在長(zhǎng)江口魚類 多樣性分析中的應(yīng)用[J].水生生物學(xué)報(bào),47(3):365-375.

    WANG R X, YANG G, GENG Z, et al, 2023. Application of environmental DNA technology in fish diversity analysis in the Yangtze River estuary[J]. Acta Hydrobiologica Sinica, 47(3):365-375.

    王銀肖,楊慧蘭,譚慧敏,等,2022.白洋淀魚類群落結(jié)構(gòu)與環(huán) 境因子關(guān)系分析[J].上海海洋大學(xué)學(xué)報(bào),31(6):1488- 1501.

    WANG Y X, YANG H L, TAN H M, et al, 2022. Fish community structure and its relationship with environmental factors in Baiyangdian Lake[J]. Journal of Shanghai Ocean University, 31(6):1488-1501.

    張方方,2014.巢湖流域生態(tài)系統(tǒng)健康評(píng)價(jià)[D].蕪湖:安徽師 范大學(xué).

    DUDGEON D, 2019. Multiple threats imperil freshwater biodiversity in the Anthropocene[J]. Current Biology, 29(19): R960-R967.

    EMMETT DUFFY J, LEFCHECK J S, STUART-SMITH R D, et al 2n16 Rindivercitv enhanree reef fich hinmacc and resistance to climate change[J]. Proceedings of the National Academy of Sciences of the United States of America, 113(22):6230-6235.

    FANG T, LU WX,LI J,et al,2017. Levels and risk assessment of metals in sediment and fish from Chaohu Lake, Anhui Province,China[J].Environmental Science and Pollution Research International, 24(18):15390-15400.

    HARPER L R,BUXTONAS,REES HC, et al,2019.Prospects and challenges of environmental DNA (eDNA) monitoring in freshwater ponds[J].Hydrobiologia, 826(1):25-41.

    HIGGINS J, ZABLOCKI J,NEWSOCK A, et al, 2021. urable freshwater protection: a framework for establishing and maintaining long-term protection for freshwater ecosystems and the values they sustain[J]. Sustainability,13 (4):1950.

    KONG X Z, JORGENSEN S E, HE W, et al, 2013. Predicting therestoration effectsby a structural dynamic approachin Lake Chaohu, China[J]. Ecological Modelling,266:73-85.

    LARSONJH, TREBITZ A S, STEINMANAD, et al,2013. Great Lakes rivermouth ecosystems: scientific synthesis and management implications[J]. Journal of Great Lakes Research, 39(3):513-524.

    PHANG S C, COOPERMAN M, LYNCH A J, et al, 2019. Fishing for conservation of freshwater tropical fishes in the Anthropocene[J].Aquatic Conservation: Marine and Freshwater Ecosystems, 29(7):1039-1051.

    PONT D,ROCLE M, VALENTINI A, et al, 2018. Environmental DNA reveals quantitative patterns of fish biodiversity in large rivers despite its downstream transportation[J]. Scientific Reports, 8(1):10361.

    REES H C,MADDISON B C,MIDDLEDITCH D J,et al, 2014.REVIEW: the detection of aquatic animal species using environmental DNA-a review of eDNA as a survey tool in ecology[J]. Journal of Applied Ecology, 51(5): 1450-1459.

    RUPPERT K M, KLINE R J, RAHMAN M S,2019. Past, present,and future perspectives of environmental DNA (eDNA)metabarcoding: a systematic review in methods, monitoring,and applications of global eDNA[J]. Global Ecology and Conservation,17:e00547.

    ROTHERHAM D, CHAPMAN M,UNDERWOOD A, et al, 2011.Untangling spatial and temporal variation in abundances of estuarine fish sampled with multi-mesh gillnets [J].Marine Ecology Progress Series, 435:183-195.

    TKACHUK K A, DUNN D A, 2020. Analysis of sea lamprey environmental DNA during lampricide treatment in a tributary of Lake Ontario[J].Knowledge amp; Management of Aquatic Ecosystems, 421:14.

    THOMSEN PF, WILLERSLEV E, 2015. Environmental DNAan emerging tool in conservation for monitoring past and present biodiversity[J].Biological Conservation,183:4-18.

    LAMY T,PITZKJ,CHAVEZ FP, etal,2021. Environmental DNA reveals the fine-grained and hierarchical spatial structure of kelp forest fish communities[J]. Scientific Reports,11(1):14439.

    WORM B, BARBIER E B, BEAUMONT N, et al, 2006. Impacts of biodiversity loss on ocean ecosystem services[J]. Science, 314(5800):787-790.

    YU HB,XI B D, JIANG JY, et al, 2011. Environmental heterogeneity analysis, assessment of trophic state and source identification in Chaohu Lake,China[J]. Environmental Science and Pollution Research,18(8):1333-1342.

    YAO JP,HE G L,2020. Design of computational model for cyanobacterial pollutant diffusion change in Chaohu Lake of China based on large data[J]. Discrete Dynamics in Nature and Society, (1):2036349.

    YANG L B, LEI K, MENG W, et al, 2013. Temporal and spatial changes in nutrients and chlorophyll- a in a shallow lake,Lake Chaohu, China:an 11-year investigation[J]. Journal of Environmental Sciences,25(6):1117-1123.

    ZOU K S, CHENJW,RUANHT, et al,2020. eDNA metabarcoding as a promising conservation tool for monitoring fish diversity in a coastal wetland of the Pearl River Estuary compared to bottom trawling[J]. Science of the Total Environment, 702:134704.

    ZHANG H, YOSHIZAWA S, IWASAKI W, et al, 2019. Seasonal fish assemblage structure using environmental DNA in the Yangtze estuary and its adjacent waters[J]. Frontiers in Marine Science, 11(37):33931-33940.

    (責(zé)任編輯 熊美華 崔莎莎)

    Summer Fish Community Structure and Diversity in Chaohu Lake Based on eDNA Metabarcoding Technology

    LIU Chengfeng1.2, YAN Yunzhi1, HAN Shengpan2, ZHANG Jiahe23, XIE Ping2, CHEN Jun2, XIA Wulai2

    (1. Collaborative Innovation Center of Recovery and Reconstruction of Degraded Ecosystem in Wanjiang Basin Co-founded by Anhui Province and Ministry of Education, School of Ecology and Environment, Anhui Normal University, Wuhu 2410oo, P.R. China ; 2. Donghu Experimental Station of Lake Ecosystems, Institute of Hydrobiology, Chinese Academy of Sciences,Wuhan 430072, P.R. China; 3. College of Ocean Science,Agricultural University of Hebei, Qinhuangdao O66000,P.R. China)

    Abstract: In this study, we explored fish community structure, diversity and its relationship with environmental factors during summer in Chaohu Lake (Hefei, Anhui Province) using eDNA metabarcoding and redundancy analysis (RDA).We aimed to provide a new technical means of monitoring and protecting fish resources in Chaohu Lake,and provide new insights for fish resource investigation in the Yangtze River basin. In August 2O22,water samples at 15 sampling points of three different areas in Chaohu Lake were collected for eDNA high-throughput sequencing,as wellas for the determination of water environment parameters.A total of 7 040 OTUs were obtained, with an average sequence length of 170.89 bp, whichidentified48 fish species inthedatabase,belongingto39 genera,13families,7ordersandalso including two species (Squaliobarbus curriculus and Hemibarbus labeo), which were not collected by traditional survey methods reported in the literature over the past 10 years.Cypriniformes (35 species) were the dominant group, accounting for 72.92% of the total number of detected fishes.Aristichys nobilis and Coilia nasus were the dominant species. The average values of the Shannon-Wiener, Simpson and Pielou evenness indices of the 15 sampling points were, respectively, 1.710±0.375 , 0.696±0.140 and 0.504±0.112 The spatial differences in fish community diversity was not significant in Chaohu Lake, with similar distributions in all three lake areas. Principal coordinate analysis (PCoA)and permutation multivariate analysis of variance (PERMANOVA) showed no significant diferences in fish community structure among the three areas of Chaohu Lake ( R2=0.13 , P=0.53 ). Redundancy analysis shows that ammonia nitrogen, water temperature and dissolved oxygen significantly affected fish community structure in Chaohu Lake (Plt;0.050) ). Our research demonstrates that eDNA metabarcoding technology is feasible and very useful for monitoring fish resources in Chaohu Lake.

    Key words: eDNA metabarcoding technology; fish diversity; community structure; environmental factors; ChaohuLake

    猜你喜歡
    巢湖魚類群落
    上海海洋大學(xué)解碼魚類胚胎中的新生轉(zhuǎn)錄本
    鐵對(duì)短程反硝化的性能影響研究進(jìn)展
    遼寧化工(2025年7期)2025-08-18 00:00:00
    不同基肥及施用量對(duì)設(shè)施大棚土壤養(yǎng)分和微生物的影響
    不同海拔文冠果根系性狀、土壤理化性質(zhì)與叢枝菌根真菌群落特征的相關(guān)性
    基于高通量測(cè)序?qū)Σ煌鑸@土壤真菌群落結(jié)構(gòu)特征的分析
    遷徙大軍的定位密碼
    基于計(jì)算機(jī)視覺的魚類智能投喂方法研究進(jìn)展
    流速-水溫對(duì)草魚生殖洄游的交互作用
    基于線粒體CO/基因序列的中國(guó)蜆屬遺傳多樣性和結(jié)構(gòu)分析
    亚洲aⅴ乱码一区二区在线播放 | 国产成人精品在线电影| 一进一出抽搐gif免费好疼 | 两性午夜刺激爽爽歪歪视频在线观看 | 中文字幕高清在线视频| 两个人看的免费小视频| 久久久国产一区二区| 大码成人一级视频| 成年女人毛片免费观看观看9 | 国产激情欧美一区二区| 久久人人爽av亚洲精品天堂| 亚洲精品国产色婷婷电影| 可以免费在线观看a视频的电影网站| 老司机午夜十八禁免费视频| 日韩有码中文字幕| 中国美女看黄片| 18禁裸乳无遮挡免费网站照片 | 亚洲人成77777在线视频| 国产精品九九99| 捣出白浆h1v1| 岛国毛片在线播放| 啦啦啦 在线观看视频| videosex国产| 国产一区二区三区在线臀色熟女 | 国产精品二区激情视频| 99精品欧美一区二区三区四区| 亚洲男人天堂网一区| 亚洲熟妇熟女久久| 日韩三级视频一区二区三区| 精品久久久久久久毛片微露脸| 国产精品av久久久久免费| 精品免费久久久久久久清纯 | 久久久久久久久免费视频了| 黑人欧美特级aaaaaa片| 免费观看精品视频网站| 亚洲欧美精品综合一区二区三区| 脱女人内裤的视频| 99久久99久久久精品蜜桃| 中文字幕精品免费在线观看视频| 亚洲自偷自拍图片 自拍| 老司机亚洲免费影院| 一边摸一边抽搐一进一出视频| 水蜜桃什么品种好| 97人妻天天添夜夜摸| 一级a爱视频在线免费观看| 国产亚洲一区二区精品| 亚洲全国av大片| 天天躁夜夜躁狠狠躁躁| 黑人巨大精品欧美一区二区蜜桃| 国产男靠女视频免费网站| 女人高潮潮喷娇喘18禁视频| 亚洲全国av大片| 999久久久国产精品视频| 午夜视频精品福利| 中文字幕制服av| 99久久人妻综合| 国产高清视频在线播放一区| 国产一区二区三区视频了| 久久久久久免费高清国产稀缺| 免费在线观看影片大全网站| 日本五十路高清| 久久性视频一级片| 别揉我奶头~嗯~啊~动态视频| 一进一出抽搐动态| 色婷婷av一区二区三区视频| 久久久国产成人精品二区 | 黄色片一级片一级黄色片| 欧美丝袜亚洲另类 | 纯流量卡能插随身wifi吗| 欧美亚洲日本最大视频资源| 国产无遮挡羞羞视频在线观看| 成人国产一区最新在线观看| 欧美日韩亚洲高清精品| 精品一区二区三区视频在线观看免费 | 久久久久久久精品吃奶| 久久久久国内视频| 欧美乱妇无乱码| 热99re8久久精品国产| 久久久国产成人精品二区 | 最新在线观看一区二区三区| 亚洲第一欧美日韩一区二区三区| 欧美成人免费av一区二区三区 | 热99国产精品久久久久久7| 老司机在亚洲福利影院| 欧美亚洲日本最大视频资源| 欧美+亚洲+日韩+国产| 搡老岳熟女国产| 国产高清激情床上av| 99热只有精品国产| 亚洲性夜色夜夜综合| 亚洲一码二码三码区别大吗| 成人黄色视频免费在线看| 很黄的视频免费| 欧美另类亚洲清纯唯美| 黄片播放在线免费| 免费一级毛片在线播放高清视频 | 精品国产一区二区三区久久久樱花| 亚洲一码二码三码区别大吗| 18在线观看网站| 欧美成狂野欧美在线观看| 99热网站在线观看| 欧美一级毛片孕妇| 久久热在线av| 色综合婷婷激情| 一二三四在线观看免费中文在| 免费在线观看影片大全网站| 女人高潮潮喷娇喘18禁视频| 一a级毛片在线观看| 亚洲熟妇熟女久久| 亚洲欧美一区二区三区黑人| 国产精华一区二区三区| 免费在线观看亚洲国产| 大片电影免费在线观看免费| 制服人妻中文乱码| 在线国产一区二区在线| 两性午夜刺激爽爽歪歪视频在线观看 | 欧美激情 高清一区二区三区| 亚洲精品美女久久av网站| 香蕉丝袜av| 国产精品 国内视频| 无人区码免费观看不卡| 人妻久久中文字幕网| svipshipincom国产片| 啦啦啦在线免费观看视频4| 国产成+人综合+亚洲专区| 欧美不卡视频在线免费观看 | 99香蕉大伊视频| 欧美黄色淫秽网站| 欧美大码av| 亚洲国产欧美网| 黄色成人免费大全| av天堂在线播放| 99国产极品粉嫩在线观看| 黄片小视频在线播放| 免费高清在线观看日韩| 丝袜人妻中文字幕| 久久久国产成人免费| 中出人妻视频一区二区| 亚洲情色 制服丝袜| 真人做人爱边吃奶动态| 一边摸一边抽搐一进一出视频| 女人被躁到高潮嗷嗷叫费观| 美女午夜性视频免费| 精品一品国产午夜福利视频| 黄色视频不卡| 亚洲欧美一区二区三区黑人| 久久性视频一级片| 久久久久国产一级毛片高清牌| 最近最新中文字幕大全电影3 | 久久人人爽av亚洲精品天堂| 一本大道久久a久久精品| 亚洲精品在线美女| videos熟女内射| 在线免费观看的www视频| 国产男女超爽视频在线观看| 18在线观看网站| 首页视频小说图片口味搜索| 搡老乐熟女国产| netflix在线观看网站| 中文亚洲av片在线观看爽 | x7x7x7水蜜桃| √禁漫天堂资源中文www| 一级毛片高清免费大全| 不卡一级毛片| 丁香欧美五月| 捣出白浆h1v1| 一级片免费观看大全| 欧美丝袜亚洲另类 | 男女之事视频高清在线观看| 麻豆av在线久日| 亚洲五月婷婷丁香| 久久久国产成人精品二区 | 国产精品永久免费网站| 国产1区2区3区精品| 一边摸一边抽搐一进一出视频| 韩国精品一区二区三区| 在线永久观看黄色视频| 中文字幕人妻丝袜一区二区| 少妇裸体淫交视频免费看高清 | 亚洲色图综合在线观看| 中文字幕精品免费在线观看视频| 美女国产高潮福利片在线看| 999精品在线视频| 成年动漫av网址| 一个人免费在线观看的高清视频| 人成视频在线观看免费观看| 操出白浆在线播放| 夫妻午夜视频| 丝袜美足系列| 国产精品免费视频内射| 人人妻,人人澡人人爽秒播| 久久天堂一区二区三区四区| 久久中文字幕一级| 国产精华一区二区三区| 国产成人欧美在线观看 | 亚洲国产精品sss在线观看 | 老司机在亚洲福利影院| 欧美日韩黄片免| 欧美精品高潮呻吟av久久| 久久久国产精品麻豆| 久久久国产成人免费| 黑人猛操日本美女一级片| av电影中文网址| 亚洲成人免费电影在线观看| 国产亚洲精品第一综合不卡| 人人妻人人爽人人添夜夜欢视频| 俄罗斯特黄特色一大片| 一进一出抽搐动态| 久久午夜综合久久蜜桃| 中文字幕精品免费在线观看视频| 黄网站色视频无遮挡免费观看| 午夜影院日韩av| 亚洲av成人av| 久久性视频一级片| 老司机午夜福利在线观看视频| 亚洲av片天天在线观看| 午夜两性在线视频| 老熟女久久久| 午夜福利一区二区在线看| 香蕉丝袜av| 99热国产这里只有精品6| 久久久久精品国产欧美久久久| netflix在线观看网站| 亚洲中文字幕日韩| 欧美 日韩 精品 国产| 他把我摸到了高潮在线观看| 欧美日韩亚洲综合一区二区三区_| 免费观看a级毛片全部| 国内久久婷婷六月综合欲色啪| 看黄色毛片网站| 久久久久久免费高清国产稀缺| 久久精品aⅴ一区二区三区四区| 亚洲av日韩精品久久久久久密| 成人手机av| 久久人妻av系列| 国产精品免费视频内射| 国产在线观看jvid| 露出奶头的视频| 久99久视频精品免费| 亚洲在线自拍视频| 麻豆成人av在线观看| 精品国产美女av久久久久小说| 亚洲九九香蕉| 男人的好看免费观看在线视频 | 亚洲精品av麻豆狂野| 亚洲熟女毛片儿| 黄片小视频在线播放| 亚洲精品自拍成人| 夫妻午夜视频| 欧美不卡视频在线免费观看 | aaaaa片日本免费| 欧美最黄视频在线播放免费 | 在线播放国产精品三级| 下体分泌物呈黄色| 亚洲片人在线观看| a级毛片在线看网站| 色综合欧美亚洲国产小说| av电影中文网址| 波多野结衣av一区二区av| 亚洲熟妇熟女久久| 久久久国产成人免费| 久久热在线av| 亚洲三区欧美一区| 丰满的人妻完整版| 一边摸一边抽搐一进一小说 | 新久久久久国产一级毛片| 色94色欧美一区二区| 性少妇av在线| 国产成+人综合+亚洲专区| av线在线观看网站| 91精品三级在线观看| 国产av又大| 国精品久久久久久国模美| 嫩草影视91久久| 久久国产精品大桥未久av| 极品人妻少妇av视频| 韩国av一区二区三区四区| 欧美av亚洲av综合av国产av| 国产亚洲精品久久久久5区| 18禁裸乳无遮挡动漫免费视频| 91九色精品人成在线观看| 久久久国产欧美日韩av| 国产av精品麻豆| 久久香蕉激情| 国产99久久九九免费精品| 国产欧美日韩一区二区三区在线| 久久久久久久午夜电影 | 无人区码免费观看不卡| 国产1区2区3区精品| 欧美日韩亚洲高清精品| 黄色视频不卡| 99久久国产精品久久久| 久久 成人 亚洲| 18禁观看日本| 五月开心婷婷网| 99热国产这里只有精品6| 99精品在免费线老司机午夜| 亚洲av成人不卡在线观看播放网| 人人妻人人澡人人爽人人夜夜| 国产麻豆69| 黑人操中国人逼视频| 村上凉子中文字幕在线| 大型av网站在线播放| 成人影院久久| 一进一出抽搐gif免费好疼 | 欧美不卡视频在线免费观看 | 亚洲第一欧美日韩一区二区三区| 久久久久国内视频| 在线av久久热| 欧美乱码精品一区二区三区| 午夜福利视频在线观看免费| 51午夜福利影视在线观看| 久久久国产一区二区| 国产乱人伦免费视频| 宅男免费午夜| 亚洲五月婷婷丁香| 免费在线观看日本一区| 成人国产一区最新在线观看| 大陆偷拍与自拍| 精品欧美一区二区三区在线| 亚洲va日本ⅴa欧美va伊人久久| 久久国产精品男人的天堂亚洲| 夜夜爽天天搞| www.熟女人妻精品国产| 在线观看免费视频日本深夜| 亚洲,欧美精品.| 一级作爱视频免费观看| 丝袜美足系列| 国产精品久久久av美女十八| 精品福利永久在线观看| 日韩欧美在线二视频 | 午夜亚洲福利在线播放| 高清av免费在线| 一二三四在线观看免费中文在| 国产精品免费一区二区三区在线 | 两个人看的免费小视频| 久99久视频精品免费| 久久 成人 亚洲| 精品亚洲成a人片在线观看| 狠狠狠狠99中文字幕| 美女扒开内裤让男人捅视频| 国产亚洲av高清不卡| 国产一区二区三区综合在线观看| 午夜老司机福利片| 色播在线永久视频| 欧美av亚洲av综合av国产av| 男女之事视频高清在线观看| 国产激情久久老熟女| 一级a爱视频在线免费观看| 99热只有精品国产| 久99久视频精品免费| 9191精品国产免费久久| 免费看a级黄色片| 交换朋友夫妻互换小说| 夫妻午夜视频| 女同久久另类99精品国产91| 亚洲精品一二三| 久久久国产成人免费| 亚洲午夜理论影院| 啦啦啦免费观看视频1| 欧美亚洲 丝袜 人妻 在线| 制服人妻中文乱码| 美女视频免费永久观看网站| 亚洲性夜色夜夜综合| 熟女少妇亚洲综合色aaa.| 国产99白浆流出| 视频区欧美日本亚洲| 人妻久久中文字幕网| 在线观看免费视频日本深夜| www.999成人在线观看| 午夜精品在线福利| 午夜日韩欧美国产| 国产精品久久久人人做人人爽| 在线观看免费日韩欧美大片| 一级毛片高清免费大全| 国产又爽黄色视频| 一边摸一边抽搐一进一小说 | 久久中文看片网| 操出白浆在线播放| 香蕉久久夜色| 午夜福利一区二区在线看| 在线观看日韩欧美| 欧美国产精品va在线观看不卡| 80岁老熟妇乱子伦牲交| 在线观看免费午夜福利视频| 中文字幕精品免费在线观看视频| 精品国产美女av久久久久小说| 免费在线观看黄色视频的| 欧美日韩瑟瑟在线播放| 亚洲男人天堂网一区| 日本wwww免费看| 九色亚洲精品在线播放| 亚洲精品国产精品久久久不卡| 成人国语在线视频| 国产精品免费视频内射| 女人久久www免费人成看片| 高清在线国产一区| 亚洲精品在线观看二区| 亚洲国产精品一区二区三区在线| av天堂久久9| 91成人精品电影| 最新的欧美精品一区二区| 在线观看免费高清a一片| 人妻 亚洲 视频| 国产熟女午夜一区二区三区| 一级毛片高清免费大全| 久久狼人影院| 免费在线观看黄色视频的| 美国免费a级毛片| 欧美激情 高清一区二区三区| 亚洲人成伊人成综合网2020| 久热爱精品视频在线9| 欧美黄色片欧美黄色片| 热99re8久久精品国产| 又黄又粗又硬又大视频| 在线免费观看的www视频| 两人在一起打扑克的视频| 亚洲欧美日韩另类电影网站| 成人免费观看视频高清| svipshipincom国产片| 久久久久久久国产电影| 精品福利永久在线观看| 中文字幕av电影在线播放| 午夜福利欧美成人| 成年人黄色毛片网站| 不卡av一区二区三区| 岛国毛片在线播放| 人妻久久中文字幕网| 欧美亚洲日本最大视频资源| 香蕉丝袜av| 99精品在免费线老司机午夜| avwww免费| 日本撒尿小便嘘嘘汇集6| 国产无遮挡羞羞视频在线观看| 国产精品自产拍在线观看55亚洲 | 国产99白浆流出| 国产成人精品久久二区二区91| 99精品在免费线老司机午夜| 80岁老熟妇乱子伦牲交| 亚洲人成电影免费在线| 王馨瑶露胸无遮挡在线观看| 亚洲视频免费观看视频| 操美女的视频在线观看| 超碰97精品在线观看| 啪啪无遮挡十八禁网站| 精品国产国语对白av| 久久久久精品人妻al黑| 日韩欧美在线二视频 | 在线观看免费午夜福利视频| 一本大道久久a久久精品| 国产成人系列免费观看| 亚洲黑人精品在线| 欧美久久黑人一区二区| 一本大道久久a久久精品| 成人特级黄色片久久久久久久| 中文字幕另类日韩欧美亚洲嫩草| 脱女人内裤的视频| 国产蜜桃级精品一区二区三区 | 国产高清视频在线播放一区| 中亚洲国语对白在线视频| 亚洲va日本ⅴa欧美va伊人久久| 亚洲精品国产精品久久久不卡| 国产片内射在线| 日韩成人在线观看一区二区三区| 欧美人与性动交α欧美精品济南到| 97人妻天天添夜夜摸| 免费女性裸体啪啪无遮挡网站| 日日摸夜夜添夜夜添小说| 热99久久久久精品小说推荐| 精品久久蜜臀av无| 黄色 视频免费看| 最近最新中文字幕大全免费视频| 在线观看免费高清a一片| 欧美黑人欧美精品刺激| 免费av中文字幕在线| 日韩成人在线观看一区二区三区| 日本一区二区免费在线视频| 建设人人有责人人尽责人人享有的| 天天操日日干夜夜撸| 9191精品国产免费久久| 国产欧美日韩精品亚洲av| 亚洲欧洲精品一区二区精品久久久| 欧美精品高潮呻吟av久久| 黄色片一级片一级黄色片| 身体一侧抽搐| 日韩大码丰满熟妇| 欧美精品啪啪一区二区三区| 亚洲av电影在线进入| 久久 成人 亚洲| 日本精品一区二区三区蜜桃| 黑人猛操日本美女一级片| 涩涩av久久男人的天堂| 日本wwww免费看| 老汉色∧v一级毛片| 超色免费av| 老司机在亚洲福利影院| 日本a在线网址| 亚洲午夜理论影院| 成人亚洲精品一区在线观看| 亚洲精品久久午夜乱码| 欧美黄色片欧美黄色片| 成人18禁在线播放| 免费少妇av软件| 狂野欧美激情性xxxx| 午夜免费成人在线视频| bbb黄色大片| 亚洲 国产 在线| 麻豆国产av国片精品| 婷婷成人精品国产| 高清毛片免费观看视频网站 | 高清在线国产一区| 国产主播在线观看一区二区| 久久久久国产一级毛片高清牌| videos熟女内射| 夜夜夜夜夜久久久久| 午夜免费成人在线视频| 国产成人一区二区三区免费视频网站| 国产精品综合久久久久久久免费 | 美女国产高潮福利片在线看| 亚洲欧美一区二区三区黑人| 老司机深夜福利视频在线观看| 日韩视频一区二区在线观看| 成人精品一区二区免费| 欧美日韩中文字幕国产精品一区二区三区 | 日韩欧美在线二视频 | 国产激情久久老熟女| 国产精品香港三级国产av潘金莲| 久久久久国产一级毛片高清牌| 亚洲国产毛片av蜜桃av| 十分钟在线观看高清视频www| 村上凉子中文字幕在线| 亚洲人成伊人成综合网2020| 18禁观看日本| av福利片在线| 中文字幕制服av| 国产极品粉嫩免费观看在线| 中文字幕制服av| 视频区图区小说| 久久人人97超碰香蕉20202| 操出白浆在线播放| 国产精品影院久久| 99久久国产精品久久久| 中文字幕制服av| 一级作爱视频免费观看| 午夜精品国产一区二区电影| 一区在线观看完整版| 国产精品久久久久久精品古装| 超色免费av| 日本黄色视频三级网站网址 | 黑人巨大精品欧美一区二区mp4| 亚洲精品自拍成人| 亚洲专区国产一区二区| 国产黄色免费在线视频| 亚洲五月色婷婷综合| 男人的好看免费观看在线视频 | 男女高潮啪啪啪动态图| 久热这里只有精品99| 欧美中文综合在线视频| 久久精品aⅴ一区二区三区四区| √禁漫天堂资源中文www| 国产av一区二区精品久久| 久久影院123| 一本一本久久a久久精品综合妖精| 91老司机精品| 日本黄色日本黄色录像| 国产1区2区3区精品| 日韩熟女老妇一区二区性免费视频| 在线观看免费高清a一片| 精品久久久久久久毛片微露脸| 咕卡用的链子| 精品国产乱码久久久久久男人| 国产亚洲精品一区二区www | 美国免费a级毛片| 国产亚洲精品久久久久5区| 在线十欧美十亚洲十日本专区| 在线看a的网站| 国产精品98久久久久久宅男小说| 免费看十八禁软件| 亚洲精品国产区一区二| 免费观看a级毛片全部| 精品电影一区二区在线| 精品一区二区三区四区五区乱码| 亚洲伊人色综图| 久久精品国产99精品国产亚洲性色 | 日韩 欧美 亚洲 中文字幕| 大码成人一级视频| 国产主播在线观看一区二区| 国产精品欧美亚洲77777| 日韩欧美三级三区| av国产精品久久久久影院| 成在线人永久免费视频| 又黄又粗又硬又大视频| 交换朋友夫妻互换小说| 在线观看66精品国产| 久久九九热精品免费| 一进一出抽搐gif免费好疼 | 老司机午夜十八禁免费视频| 午夜影院日韩av| 老司机午夜福利在线观看视频| 精品一品国产午夜福利视频| 欧美精品av麻豆av| 中文字幕人妻丝袜一区二区| 黄色怎么调成土黄色| 亚洲中文av在线| 好看av亚洲va欧美ⅴa在| 一夜夜www| 色老头精品视频在线观看| 久久久久国内视频| 亚洲精品国产区一区二| 国产精品秋霞免费鲁丝片| 天天躁狠狠躁夜夜躁狠狠躁| 欧美日韩亚洲国产一区二区在线观看 | 亚洲第一av免费看| 黑人猛操日本美女一级片| a级片在线免费高清观看视频|