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

    基于eDNA技術(shù)的長(zhǎng)江下游典型江段浮游植物生物完整性評(píng)價(jià)

    2025-04-07 00:00:00周宇恒趙旭王玉清陳麗麗章守宇王凱

    DOI:10.3969/j.issn.2095-1191.2025.01.024

    摘要:【目的】基于環(huán)境DNA(eDNA)技術(shù)評(píng)價(jià)長(zhǎng)江下游典型江段浮游植物生物完整性,為評(píng)估長(zhǎng)江大保護(hù)實(shí)施后的生態(tài)恢復(fù)效果提供理論依據(jù)。【方法】分別在長(zhǎng)江下游安慶(AQ)、蕪湖(EH)、靖江(JJ)3個(gè)典型江段豐水期(2021年9月)和枯水期(2021年4月)采集eDNA樣本。在3個(gè)典型江段選擇參考點(diǎn)和受損點(diǎn),通過(guò)篩選物種豐富度、群落物種豐度、群落結(jié)構(gòu)組成和群落營(yíng)養(yǎng)結(jié)構(gòu),選擇5個(gè)候選指標(biāo),然后采用0~5評(píng)分法構(gòu)建長(zhǎng)江下游典型江段浮游植物生物完整性指數(shù)(P-IBI)評(píng)價(jià)體系?!窘Y(jié)果】在安慶、蕪湖、靖江3個(gè)典型江段共鑒定出8個(gè)門(mén)96個(gè)屬浮游植物,包括硅藻門(mén)、綠藻門(mén)、藍(lán)藻門(mén)、金藻門(mén)、甲藻門(mén)、褐藻門(mén)、紅藻門(mén)、定鞭藻門(mén)和舟形藻屬、根管藻屬、衣藻屬、色球藻屬、微囊藻屬、團(tuán)藻屬和骨條藻屬等。在豐水期,安慶、蕪湖和靖江江段在屬分類(lèi)水平上浮游植物物種總數(shù)分別為4849、5746和7617種,枯水期分別下降至3933、4243和5683種。豐水期,安慶、蕪湖江段的優(yōu)勢(shì)種為硅藻門(mén)、藍(lán)藻門(mén)和綠藻門(mén);靖江江段硅藻門(mén)占主導(dǎo)地位;枯水期,3個(gè)江段的浮游植物群落均以硅藻門(mén)和綠藻門(mén)為主。在豐水期,安慶、蕪湖、靖江江段的得分分別為17、16和15,浮游植物生物完整性等級(jí)分別為好、好和一般。在枯水期,安慶、蕪湖、靖江江段的得分分別為13、16和17,浮游植物生物完整性等級(jí)分別為一般、好和好?!窘Y(jié)論】eDNA技術(shù)在物種分類(lèi)數(shù)量上展現(xiàn)出高豐富性、高物種多樣性和詳細(xì)的群落組成,對(duì)浮游植物物種檢測(cè)具有高效且準(zhǔn)確的特點(diǎn),在生物參數(shù)指標(biāo)的選取中具有良好的篩選作用,對(duì)P-IBI評(píng)價(jià)體系構(gòu)建具有良好的支撐作用。

    關(guān)鍵詞:浮游植物;eDNA技術(shù);浮游植物生物完整性指數(shù)(P-IBI);水文期;長(zhǎng)江下游

    中圖分類(lèi)號(hào):S917.3文獻(xiàn)標(biāo)志碼:A文章編號(hào):2095-1191(2025)01-0275-12

    Evaluation of biointegrity of phytoplankton in typical sections of the lower reaches of Yangtze River using eDNA technology

    ZHOU Yu-heng1,ZHAO Xu1,2,WANG Yu-qing1,2,CHEN Li-li1,2,ZHANG Shou-yu1,2,WANG Kai1,2*

    (1College of Oceanography and Ecological Science,Shanghai Ocean University,Shanghai 201306,China;2Engineering Research Center of Marine Pasture,Shanghai Ocean University,Shanghai 201306,China)

    Abstract:【Objective】To assess the phytoplankton biointegrity in typical sections of the lower reaches of Yangtze River based on environmental DNA(eDNA)technology,which could provide theoretical basis for evaluating the ecologi-cal restoration effects following the implementation of Yangtze River conservation initiative.【Method】eDNA samples were collected during wet season(September 2021)and dry season(April 2021)from 3 typical sections of the lowerreaches of Yangtze River:Anqing(AQ),Wuhu(EH)and Jingjiang(JJ).Reference and damaged points were selected inthese 3 sections.Five candidate indicators were chosen by evaluating species richness,community species abundance,community structure composition and community trophic structure.A scoring system ranging from 0 to 5 was used to con-struct an evaluation system for phytoplankton biointegrity in typical sections of the lower reaches of Yangtze River.【Result】Atotal of 96 genera from 8 phyla of phytoplankton were identified in the 3 typical sections of Anqing,Wuhu andJingjiang including Bacillariophyta,Chlorophyta,Cyanobacteria,Chrysophyta,Pyrrophyta,Phaeophyta,Rhodophyta,Haptophyta and Navicula,Rhizosolenia,Scenedesmus,Chroococcus,Microcystis,Volvox,Skeletonema.The total num-bers of phytoplankton species in Anqing,Wuhu and Jingjiang sections at genus level were 4849,5746 and 7617 respec-tively in wet season,and decreased to 3933,4243 and 5683 respectively in dry season.During wet season,the dominant species in Anqing and Wuhu sections were Bacillariophyta,Cyanobacteria and Chlorophyta.In Jingjiang section,besides the dominant Bacillariophyta was the olominant.During dry season,the phytoplankton communities in all 3 sections were primarily composed of Bacillariophyta and Chlorophyta.In wet season,the scores for Anqing,Wuhu and Jingjiang sec-tions were 17,16 and 15 respectively,with phytoplankton biointegrity grades of good,good and general.In dry season,the scores for Anqing,Wuhu and Jingjiang sections were 13,16 and 17 respectively,with phytoplankton biointegrity grades of genera,good and good.【Conclusion】eDNA technology demonstrates high richness,high species diversity and detailed community composition in terms of species classification.It proves to be efficient and accurate in detecting phyto-plankton species,plays a significant role in the selection of biological parameter indexes,and provides strong support for the construction of P-IBI evaluation system.

    Key words:phytoplankton;eDNA technology;phytoplankton biointegrity index(P-IBI);hydrological period;lower reaches of Yangtze River

    Foundation items:National Key Research and Development Program of China(2018YFD0900904)

    0引言

    【研究意義】河流生態(tài)系統(tǒng)是以河流為主體,涵蓋水體、陸地、河岸帶及周邊濕地與沼澤等的復(fù)合生態(tài)系統(tǒng)(徐宗學(xué)和馬欣洋,2023)。豐水期,河流降雨量增加,水位上升,水體鹽度降低,水溫升高,營(yíng)養(yǎng)鹽含量增加,有利于浮游植物生長(zhǎng)(Guo et al.,2015);而枯水期,河流降雨量減少,陸源有機(jī)物數(shù)量減少,水體滯留時(shí)間延長(zhǎng),不利于浮游植物生長(zhǎng)(邱陽(yáng)凌等,2018)。浮游植物能反映水體生態(tài)系統(tǒng)的健康狀況(房冬梅,2015;Stamou et al.,2019;田豐歌等,2023),適宜作為評(píng)估河流生態(tài)健康的指示物種(錢(qián)爍統(tǒng)等,2018;肖善勢(shì)等,2021;李琦,2023;盛祥銳等,2024),其群落結(jié)構(gòu)特征常用于構(gòu)建生物完整性指數(shù)(Index of biotic integrity,IBI)以評(píng)估河流生態(tài)系統(tǒng)的健康狀況。【前人研究進(jìn)展】IBI是指經(jīng)演化后生物群落在結(jié)構(gòu)和功能方面與其所處的區(qū)域環(huán)境相適應(yīng)的特性,反映了環(huán)境因子與生物參數(shù)的定量關(guān)系,能通過(guò)比較參數(shù)值與系統(tǒng)標(biāo)準(zhǔn)值來(lái)評(píng)價(jià)生態(tài)系統(tǒng)的健康狀況(沈強(qiáng)等,2012)。浮游植物生物完整性指數(shù)(Phytoplanktonic index of biotic integrity,P-IBI)具有對(duì)水生態(tài)變化敏感的特點(diǎn),已廣泛用于評(píng)價(jià)水域生態(tài)健康狀況(譚巧等,2017;黎明民等,2018;紀(jì)瑩璐等,2023)。王芳等(2022)對(duì)長(zhǎng)江下游城市湖泊的生態(tài)健康狀況進(jìn)行評(píng)估,通過(guò)冗余分析(RDA)識(shí)別影響浮游植物生物量的主要環(huán)境因素,并選取生物量、密度及Margalef指數(shù)作為核心指標(biāo)構(gòu)建P-IBI評(píng)價(jià)體系,為城市湖泊生態(tài)健康評(píng)估提供了參考;譚巧等(2023)利用P-IBI評(píng)價(jià)了長(zhǎng)江江蘇段在豐水期(7月)和枯水期(12月)的生態(tài)健康狀況,發(fā)現(xiàn)長(zhǎng)江江蘇段水環(huán)境較好;廖延科等(2024)在豐水期(8月)和枯水期(12月)構(gòu)建了修河中下游P-IBI并評(píng)價(jià)其生態(tài)健康狀況,整體評(píng)價(jià)結(jié)果為亞健康。環(huán)境DNA(eDNA)技術(shù)結(jié)合了高通量測(cè)序技術(shù)和DNA條形碼技術(shù),能快速識(shí)別群落物種組成(陳宇飛等,2022),且與傳統(tǒng)技術(shù)相比,其成本更低、檢測(cè)速度更快、檢測(cè)結(jié)果更準(zhǔn)確(Quastet al.,2013),因而在水生生態(tài)系統(tǒng)的生物群落研究和健康狀況評(píng)價(jià)中得到廣泛應(yīng)用(DeSantis et al.,2006;Turak et al.,2017)。呂嘉誠(chéng)等(2022)在豐水期(7月)通過(guò)eDNA宏條形碼技術(shù)調(diào)查了浮游藻類(lèi)的生物多樣性和生態(tài)狀況,發(fā)現(xiàn)浮游藻類(lèi)污染程度為中等;姜山等(2023)采用eDNA技術(shù)監(jiān)測(cè)了沙潁河豐水期和枯水期硅藻群落生態(tài)健康狀況,發(fā)現(xiàn)豐水期硅藻群落生態(tài)健康狀況優(yōu)于枯水期?!颈狙芯壳腥朦c(diǎn)】長(zhǎng)江流域水生生物資源豐富,是我國(guó)生物多樣性最具代表性的地區(qū)之一。長(zhǎng)江下游浮游植物群落門(mén)類(lèi)多樣,涵蓋硅藻門(mén)、藍(lán)藻門(mén)、綠藻門(mén)、甲藻門(mén)等多個(gè)類(lèi)群,既包括淡水種類(lèi),也有海水種及半咸水種(張瀟月,2018)。近年來(lái),人類(lèi)日益增加的涉水活動(dòng)對(duì)長(zhǎng)江水文節(jié)律和生態(tài)環(huán)境產(chǎn)生嚴(yán)重影響,導(dǎo)致浮游植物物種數(shù)量和多樣性急劇下降(楊桂山和徐昔保,2020;王迪等,2024)。雖然農(nóng)業(yè)農(nóng)村部自2021年1月1日實(shí)施為期10年的禁捕政策,但棲息環(huán)境的改善對(duì)浮游植物群落結(jié)構(gòu)及其多樣性變化的具體影響仍需進(jìn)一步研究。【擬解決的關(guān)鍵問(wèn)題】利用eDNA技術(shù)構(gòu)建P-IBI評(píng)價(jià)體系,對(duì)長(zhǎng)江下游安慶(AQ)、蕪湖(WH)和靖江(JJ)3個(gè)典型江段在枯水期(2021年4月)與豐水期(2021年9月)的浮游植物生物完整性進(jìn)行評(píng)價(jià),為評(píng)估長(zhǎng)江大保護(hù)實(shí)施后的生態(tài)恢復(fù)效果提供理論依據(jù)。

    1材料與方法

    1.1 eDNA分析與處理

    1.1.1采樣點(diǎn)位及eDNA樣本采集在長(zhǎng)江下游安慶(30°30′0″N,117°3′0″E)、蕪湖(31°20′0″N,118°21′30″E)、靖江(32°1′0″N,120°16′30″E)3個(gè)典型江段各設(shè)6個(gè)站點(diǎn)(圖1),于長(zhǎng)江下游枯水期和豐水期采集水樣。從各站點(diǎn)表、中、底3個(gè)水層各采集500 mL水樣,每500 mL水樣加入5 mL魯格氏液,靜置48 h,采用虹吸法濃縮至30 mL。每次取100μL濃縮水樣滴入浮游生物計(jì)數(shù)框中,在光學(xué)顯微鏡下進(jìn)行計(jì)數(shù),每個(gè)樣品重復(fù)3次,計(jì)數(shù)視野≥100個(gè),最終計(jì)數(shù)平均值間的差異≤15%。

    1.1.2 PCR擴(kuò)增使用引物Euk-F(5'-AACCTGG TTGATCCTGCCAGT-3')和Euk-R(5'-GATCCTTCT GCAGGTTCACCTAC-3')(王汝賢等,2023)對(duì)eDNA樣本進(jìn)行擴(kuò)增,擴(kuò)增區(qū)域?yàn)閂4區(qū)。PCR反應(yīng)體系24.0μL:2×PCR Mix 10.0μL,dNTP 4.0μL,上、下游引物各1.0μL,DNA模板5.0μL,ddH2O 3.0μL。擴(kuò)增程序:95℃預(yù)變性5 min;95℃30 s,50℃30 s,72℃3 min,進(jìn)行25個(gè)循環(huán);72℃延伸10 min。每個(gè)樣品設(shè)3個(gè)生物學(xué)重復(fù),使用無(wú)菌水作為陰性對(duì)照。擴(kuò)增產(chǎn)物使用AxyPrep DNA凝膠回收試劑盒[康寧生命科學(xué)(吳江)有限公司]進(jìn)行純化回收。

    1.1.3高通量測(cè)序純化回收產(chǎn)物中添加Phix(Illumina,美國(guó))構(gòu)建測(cè)序文庫(kù),采用Illumina HiSeqTM 2500平臺(tái)進(jìn)行高通量測(cè)序。測(cè)序獲得的原始讀序(Raw reads)經(jīng)FastQC和Trimmomatic進(jìn)行質(zhì)量評(píng)估和過(guò)濾(王晨等,2022),去除接頭、標(biāo)簽序列、低質(zhì)量及含N堿基較多的Reads,從而得到高質(zhì)量讀序(Clean reads)。利用FLASH對(duì)雙端測(cè)序的Reads進(jìn)行拼接,再以USEARCH對(duì)拼接后的序列進(jìn)行聚類(lèi),并篩除引物錯(cuò)配超過(guò)2個(gè)堿基、長(zhǎng)度短于100bp或總堿基錯(cuò)誤率超過(guò)2%的Reads,以獲得優(yōu)化后的Clean reads。通過(guò)18S rRNA高通量測(cè)序,獲得平均長(zhǎng)度為450 bp的Clean reads。在此基礎(chǔ)上,根據(jù)最小樣本序列數(shù)進(jìn)行標(biāo)準(zhǔn)化(抽平),并將序列歸類(lèi)為操作分類(lèi)單元(Operational taxonomic units,OTU),同時(shí)提取代表序列。

    1.1.4生物信息學(xué)分析按97%的相似水平對(duì)Clean reads進(jìn)行OTU聚類(lèi)分析,并基于OTU計(jì)算群落多樣性指標(biāo),包括Shannon-wiener多樣性指數(shù)(H′)、Margalef豐富度指數(shù)(D)和Pielou均勻度指數(shù)(J)。為獲取每個(gè)OTU對(duì)應(yīng)的分類(lèi)信息,采用RDP classifier貝葉斯算法對(duì)97%相似水平的OTU代表序列進(jìn)行分析,并輸入Silva和Protist Ribosomal Reference(PR2)數(shù)據(jù)庫(kù)進(jìn)行比對(duì),以獲得OTU的物種注釋信息(胡愈炘等,2021)。使用H′、D、J來(lái)表征浮游植物群落多樣性變化,并以eDNA豐度衡量浮游植物在屬分類(lèi)水平上的數(shù)量分布情況。計(jì)算公式如下:

    式中,Pi為第i種浮游植物屬水平數(shù)量占總屬水平數(shù)量(N)的比例,S為樣點(diǎn)浮游植物屬水平數(shù)量。

    1.2長(zhǎng)江下游典型江段P-IBI構(gòu)建

    1.2.1參考點(diǎn)與受損點(diǎn)選取根據(jù)人類(lèi)干擾程度,將監(jiān)測(cè)站點(diǎn)劃分為參考點(diǎn)和受損點(diǎn)。物種多樣性是評(píng)估水體健康狀況的重要生態(tài)指標(biāo),多樣性指數(shù)越低,表明人類(lèi)干擾程度越高(鄭丙輝等,2007)。以H′為分類(lèi)標(biāo)準(zhǔn),H'≥2的站點(diǎn)定義為參考點(diǎn),反之則歸為受損點(diǎn)(王備新等,2006;劉凌等,2020)。

    1.2.2候選指標(biāo)篩選從浮游植物物種豐富度、群落物種豐度、群落結(jié)構(gòu)組成及群落營(yíng)養(yǎng)結(jié)構(gòu)4個(gè)方面進(jìn)行候選指標(biāo)的篩選,篩選過(guò)程包括:(1)分布范圍分析。若某指標(biāo)在超過(guò)90%的采樣點(diǎn)中未被檢測(cè)到(數(shù)值為0),則認(rèn)為其分布范圍過(guò)于局限,予以剔除。(2)判別能力分析。采用箱線(xiàn)圖分析各指標(biāo)在參考點(diǎn)與受損點(diǎn)25%~75%分位數(shù)范圍重疊情況,以此評(píng)估判別能力,依據(jù)重疊程度賦值如下:完全無(wú)重疊,賦值3;存在重疊但中位線(xiàn)未進(jìn)入對(duì)方箱體,賦值2;存在重疊且僅有一方中位線(xiàn)進(jìn)入對(duì)方箱體,賦值1;雙方中位線(xiàn)均位于對(duì)方箱體,賦值0;僅保留判別能力指數(shù)(IQ)≥2的指標(biāo)用于后續(xù)分析(蔡琨等,2016)。(3)相關(guān)分析。對(duì)篩選后的生物指標(biāo)進(jìn)行Pearson相關(guān)分析,若2個(gè)指標(biāo)的相關(guān)系數(shù)|r|gt;0.75,表示2個(gè)指標(biāo)高度相關(guān),可能提供相似的生態(tài)信息,在此情況下,僅保留其中1個(gè)代表性指標(biāo),以減少冗余。

    1.2.3 P-IBI評(píng)分計(jì)算依據(jù)各指標(biāo)值的10%、25%、50%、75%和90%分位數(shù)設(shè)定閾值,采用0~5評(píng)分法進(jìn)行量化,并對(duì)不同指標(biāo)進(jìn)行標(biāo)準(zhǔn)化處理(蔣為等,2023),累加各指標(biāo)得分,以計(jì)算長(zhǎng)江下游典型江段的P-IBI??偟梅郑╒)計(jì)算公式如下:

    式中,Mi為第i個(gè)指標(biāo)的評(píng)價(jià)得分,n為納入計(jì)算的指標(biāo)總數(shù)。

    1.2.4 P-IBI評(píng)價(jià)標(biāo)準(zhǔn)根據(jù)P-IBI評(píng)分結(jié)果,將浮游植物生物完整性劃分為五個(gè)等級(jí):0~5分為極差,6~10分為差,11~15分為一般,16~20分為好,21~25分為極好。

    2結(jié)果與分析

    2.1長(zhǎng)江下游典型江段浮游植物群落結(jié)構(gòu)

    在長(zhǎng)江下游安慶、蕪湖、靖江3個(gè)典型江段檢測(cè)到的浮游植物隸屬于8個(gè)門(mén)類(lèi),包括硅藻門(mén)(Bacil-lariophyta)、綠藻門(mén)(Chlorophyta)、藍(lán)藻門(mén)(Cyanobac-teria)、金藻門(mén)(Chrysophyta)、甲藻門(mén)(Pyrrophyta)、褐藻門(mén)(Phaeophyta)、紅藻門(mén)(Rhodophyta)和定鞭藻門(mén)(Haptophyta)。在屬分類(lèi)水平上,共鑒定出96個(gè)屬,主要包括舟形藻屬(Navicula)、根管藻屬(Rhizo-solenia)、衣藻屬(Scenedesmus)、色球藻屬(Chroo-coccus)、微囊藻屬(Microcystis)、團(tuán)藻屬(Volvox)和骨條藻屬(Skeletonema)等。浮游植物物種組成及豐度在不同水文時(shí)期存在明顯差異。豐水期,安慶、蕪湖和靖江江段在屬分類(lèi)水平上浮游植物種總數(shù)分別為4849、5746和7617種,枯水期則分別下降至3933、4243和5683種。

    豐水期(表1和圖2),安慶江段的優(yōu)勢(shì)屬包括硅藻門(mén)舟形藻屬、藍(lán)藻門(mén)微囊藻屬、綠藻門(mén)衣藻屬,其中硅藻門(mén)舟形藻屬占比最高(20.35%);蕪湖江段的硅藻門(mén)舟形藻屬仍為主要優(yōu)勢(shì)屬,占比21.71%,其次為藍(lán)藻門(mén)微囊藻屬(10.42%);靖江江段的優(yōu)勢(shì)屬涵蓋多個(gè)門(mén)類(lèi),包括硅藻門(mén)舟形藻屬和骨條藻屬、褐藻門(mén)海帶屬(Laminaria)等,其中硅藻門(mén)骨條藻屬占比最高(23.56%)??菟冢ū?和圖3),安慶江段以硅藻門(mén)舟形藻屬(28.08%)和綠藻門(mén)衣藻屬(20.55%)占據(jù)主導(dǎo)地位;蕪湖江段以綠藻門(mén)衣藻屬(30.32%)為絕對(duì)優(yōu)勢(shì)屬,其次為硅藻門(mén)根管藻屬(25.38%);靖江江段以綠藻門(mén)衣藻屬(19.19%)和硅藻門(mén)骨條藻屬(17.57%)為主要優(yōu)勢(shì)屬。

    2.2長(zhǎng)江下游典型江段P-IBI評(píng)價(jià)體系建立

    2.2.1參考點(diǎn)和受損點(diǎn)選取結(jié)果在豐水期,安慶、蕪湖和靖江江段浮游植物的H′分別為2.01、2.02和1.98,在枯水期分別為2.03、2.10和1.95。選定安慶和蕪湖江段作為參考點(diǎn),靖江江段作為受損點(diǎn)。

    2.2.2候選指標(biāo)確定在安慶、蕪湖和靖江3個(gè)典型江段的豐水期和枯水期,選取17個(gè)對(duì)水生態(tài)環(huán)境變化較敏感的候選指標(biāo)(表2)構(gòu)建長(zhǎng)江下游典型江段P-IBI評(píng)價(jià)體系。

    2.2.3候選指標(biāo)分布范圍和判別能力分析結(jié)果

    17個(gè)候選指標(biāo)分布范圍分析結(jié)果顯示,豐水期M7、M15、M16等3個(gè)候選指標(biāo)的變化幅度較大,其中M7僅提供數(shù)量信息,不適合參與構(gòu)建P-IBI評(píng)價(jià)體系,予以剔除。判別能力分析結(jié)果(圖4)顯示,M2、M3、M4、M5、M6、M8、M9、M10的IQlt;2,不滿(mǎn)足入選條件,最終選定M1、M11、M12、M13、M14和M17等6個(gè)指標(biāo)作為候選指標(biāo)??菟冢琈2、M7和M13等3個(gè)候選指標(biāo)的變化幅度較大,其中M2、M7也僅提供數(shù)量信息,予以剔除。判別能力分析結(jié)果(圖5)顯示,M3、M4、M5、M6、M8、M9、M10的IQlt;2,因此僅有M1、M11、M12、M14、M15、M17等6個(gè)指標(biāo)符合要求。

    2.2.4候選指標(biāo)間相關(guān)分析結(jié)果對(duì)豐水期和枯水期6個(gè)候選指標(biāo)進(jìn)行Pearson相關(guān)分析,結(jié)果(表3和表4)顯示,在豐水期,M11與M12、M14和M17極顯著相關(guān)(Plt;0.01,|r|gt;0.75),表明M11所表達(dá)的生態(tài)信息在6個(gè)候選指標(biāo)中具有較高的可替代性,因此剔除M11。最終確定M1、M12、M13、M14和M17共5個(gè)候選指標(biāo)用于構(gòu)建長(zhǎng)江下游典型江段豐水期P-IBI評(píng)價(jià)體系。在枯水期,M17與M11和M15極顯著相關(guān)(Plt;0.01,|r|gt;0.75),因此剔除M17。最終確定M1、M11、M12、M14和M15共5個(gè)指標(biāo)用于構(gòu)建長(zhǎng)江下游典型江段枯水期P-IBI評(píng)價(jià)體系。

    2.3 P-IBI評(píng)價(jià)結(jié)果

    根據(jù)長(zhǎng)江下游典型江段候選指標(biāo)評(píng)價(jià)標(biāo)準(zhǔn)(表5和表6),在豐水期,安慶、蕪湖和靖江江段的得分分別為17、16和15,浮游植物生物完整性等級(jí)分別為好、好和一般(表7)。在枯水期,安慶、蕪湖和靖江江段的得分分別為13、16和17,浮游植物生物完整性等級(jí)分別為一般、好和好(表8)。

    3討論

    3.1長(zhǎng)江下游典型江段浮游植物組成與差異

    eDNA技術(shù)在物種檢測(cè)中具有高度敏感性,不僅能揭示更高的物種多樣性,還能提供更詳細(xì)的群落結(jié)構(gòu)數(shù)據(jù)。吳泳江等(2023)研究發(fā)現(xiàn),eDNA技術(shù)有助于提高浮游植物的物種多樣性檢出結(jié)果;谷思雨等(2024)研究表明,eDNA技術(shù)在水生生物研究中能以非入侵且高效的方式探測(cè)物種的存在和分布。本研究應(yīng)用eDNA技術(shù)分析長(zhǎng)江下游安慶、蕪湖和靖江3個(gè)典型江段在豐水期(2021年9月)和枯水期(2021年4月)的水樣,共檢測(cè)到浮游植物8門(mén)96屬。與張瀟月(2018)對(duì)長(zhǎng)江下游干流豐水期和枯水期的研究結(jié)果、張靜等(2023)對(duì)長(zhǎng)江流域豐水期和枯水期的研究結(jié)果相比,本研究所檢測(cè)到的浮游植物門(mén)類(lèi)數(shù)和屬數(shù)均有所增加,可能與采用的浮游植物分析方法有關(guān)。

    在各江段浮游植物物種組成上,豐水期安慶和蕪湖江段以硅藻門(mén)、綠藻門(mén)和藍(lán)藻門(mén)為主,與謝滿(mǎn)華和任青松(2016)、劉明典等(2017)對(duì)長(zhǎng)江干流安慶江段和蕪湖江段的研究結(jié)果一致;而靖江江段由于受潮汐和徑流影響,生態(tài)環(huán)境不同,浮游植物物種組成有所差異,以硅藻門(mén)占主導(dǎo)地位。在枯水期,3個(gè)江段的浮游植物群落均以硅藻門(mén)和綠藻門(mén)為主??傮w來(lái)看,長(zhǎng)江下游浮游植物群落結(jié)構(gòu)在不同水文期表現(xiàn)出一定的一致性,硅藻門(mén)始終占據(jù)主導(dǎo)地位,安慶和蕪湖江段的群落組成相似,靖江江段則表現(xiàn)出不同的特征,與鐘可兒等(2022)、魏念等(2023)的研究結(jié)果相似,展示了長(zhǎng)江下游不同江段浮游植物群落結(jié)構(gòu)在水文季節(jié)性變化中的動(dòng)態(tài)特征。

    3.2長(zhǎng)江下游典型江段浮游植物生物完整性評(píng)價(jià)

    在河流生態(tài)系統(tǒng)健康評(píng)估中,候選指標(biāo)的選擇至關(guān)重要,常見(jiàn)的候選指標(biāo)包括水文水資源(蔡琨等,2016)、水質(zhì)狀況(古小超等,2023)及水生生物多樣性(劉惠秋等,2023)等。本研究利用eDNA技術(shù),選擇多個(gè)生物參數(shù)作為候選指標(biāo)檢測(cè)浮游植物的物種信息和群落結(jié)構(gòu),與傳統(tǒng)的顯微鏡觀(guān)察和分類(lèi)學(xué)鑒定方法相比,具有便捷和高效的優(yōu)勢(shì),能克服微型物種觀(guān)察上的困難。參考點(diǎn)和受損點(diǎn)的選擇是生物完整性評(píng)價(jià)體系構(gòu)建中的關(guān)鍵步驟,由于人類(lèi)活動(dòng)對(duì)環(huán)境的廣泛干擾,選擇適當(dāng)?shù)膮⒖键c(diǎn)至關(guān)重要(余景等,2017)。馬芊芊(2015)、熊梅君等(2019)通過(guò)浮游植物Shannon-wiener多樣性指數(shù),劃定了貴州高原百花水庫(kù)和長(zhǎng)江上游河段的參考點(diǎn)與受損點(diǎn)。鑒于河流生態(tài)系統(tǒng)的縱向廊道特征(胡鳳啟等,2017),本研究選擇長(zhǎng)江下游的安慶、蕪湖和靖江3個(gè)典型江段作為研究區(qū)域,并利用Shannon-wiener多樣性指數(shù)作為劃分參考點(diǎn)和受損點(diǎn)的標(biāo)準(zhǔn),評(píng)估長(zhǎng)江下游江段的浮游植物生物完整性,其研究結(jié)果能代表不同空間特征的浮游植物群落結(jié)構(gòu),也能反映水體健康狀態(tài)的多樣性變化。

    本研究結(jié)果顯示,豐水期和枯水期的浮游植物生物完整性有所不同,可能與水文期的變化有關(guān),豐水期降水量增加,水量上升,陸源物質(zhì)流向下游,營(yíng)養(yǎng)物質(zhì)增多,為浮游植物生長(zhǎng)提供了有利條件(李照等,2017)。在枯水期,由于降水量減少、蒸發(fā)作用和水量降低,安慶江段浮游植物生物完整性受到負(fù)面影響;而蕪湖、靖江江段由于位于近河口區(qū)域,受海水與淡水交互影響,浮游植物生物完整性保持較好(邵海燕等,2024)。

    利用eDNA技術(shù)檢測(cè)浮游植物群落結(jié)構(gòu)和物種多樣性不僅具有高效、精準(zhǔn)的特點(diǎn),還能捕捉到生態(tài)系統(tǒng)中隱性或稀有物種的信息,在浮游植物生物完整性評(píng)價(jià)體系中占據(jù)優(yōu)勢(shì)。與傳統(tǒng)定量方法相比,eDNA技術(shù)仍存在一定局限性。呂嘉誠(chéng)等(2022)研究表明,eDNA技術(shù)無(wú)法準(zhǔn)確反映種群數(shù)量,且對(duì)物種數(shù)據(jù)庫(kù)的完整性和準(zhǔn)確性有較高依賴(lài);陳立國(guó)等(2024)指出,通過(guò)結(jié)合傳統(tǒng)鏡檢方法,優(yōu)化實(shí)驗(yàn)設(shè)計(jì)和數(shù)據(jù)處理方式,并加強(qiáng)物種數(shù)據(jù)庫(kù)的建設(shè)與維護(hù),能提升eDNA技術(shù)在浮游植物生物完整性評(píng)價(jià)中的應(yīng)用效果;程云山等(2024)對(duì)太湖梅梁灣的研究發(fā)現(xiàn),eDNA技術(shù)在物種豐度的定量上存在一定偏差,結(jié)果可能受DNA濃度和PCR效率等因素的影響,導(dǎo)致定量結(jié)果不準(zhǔn)確。

    4結(jié)論

    eDNA技術(shù)在物種分類(lèi)數(shù)量上展現(xiàn)出高豐富性、高物種多樣性和詳細(xì)的群落組成,對(duì)浮游植物物種檢測(cè)具有高效且準(zhǔn)確的特點(diǎn),在生物參數(shù)指標(biāo)的選取中具有良好的篩選作用,對(duì)P-IBI體系構(gòu)建具有良好的支撐作用。

    參考文獻(xiàn)(References):

    蔡琨,秦春燕,李繼影,張?jiān)?,牛志春,李旭?2016.基于浮游植物生物完整性指數(shù)的湖泊生態(tài)系統(tǒng)評(píng)價(jià):以2012年冬季太湖為例[J].生態(tài)學(xué)報(bào),36(5):1431-1441.[Cai K,Qin CY,Li JY,Zhang Y,Niu Z C,Li X W.2016.Prelimi-nary study on phytoplanktonic index of biotic integrity(P-IBI)assessment for lake ecosystem health:A case of Taihu Lake in winter,2012[J].Acta Ecologica Sinica,36(5):1431-1441.]doi:10.5846/stxb201407211483.

    陳立國(guó),萬(wàn)年新,劉本文,劉國(guó)祥,張燁.2024.eDNA揭示調(diào)水對(duì)棘洪灘水庫(kù)浮游植物的影響[J].環(huán)境科學(xué)與技術(shù),47(2):80-88.[Chen L G,Wan N X,Liu B W,Liu G X,Zhang Y.2024.Effect of water diversion on phytoplankton in the Jihongtan reservoir revealed by eDNA[J].Environ-mental scienceamp;Technology,47(2):80-88.]doi:10.19672/j.cnki.1003-6504.1978.23.338.

    陳宇飛,嚴(yán)航,夏霆,姚春雨,王浩,邱堯,劉雷.2022.基于浮游動(dòng)物生物完整性指數(shù)的太湖流域生態(tài)系統(tǒng)評(píng)價(jià)[J].南京工業(yè)大學(xué)學(xué)報(bào)(自然科學(xué)版),44(3):335-343.[Chen Y F,Yan H,Xia T,Yao C Y,Wang H,Qiu Y,Liu L.2022.Ecosystem assessment of Taihu Lake Basin based on zoo-plankton index of biological integrity[J].Journal of Nan-jing University of Technology(Natural Science Edition),44(3):335-343.]doi:10.3969/j.issn.1671-7627.2022.03.013.

    程云山,任藝晨,席貽龍,劉燕茹,何虎,范方威,陳東宜,李化炳.2024.基于環(huán)境DNA技術(shù)和形態(tài)學(xué)鑒定的浮游植物多樣性比較[J].湖泊科學(xué),36(5):1336-1352.[Cheng Y S,Ren Y C,Xi Y L,Liu Y R,He H,F(xiàn)an F W,Chen D Y,Li H B.2024.Comparative study on the biodiversity of phytoplankton based on environmental DNA technology and morphological identification[J].Journal of Lake Scien-ces,36(5):1336-1352.]doi:10.18307/2024.0511.

    房冬梅.2015.長(zhǎng)江中游宜昌至武穴江段浮游生物群落結(jié)構(gòu)研究[D].南京:南京農(nóng)業(yè)大學(xué).[Fang D M.2015.Annual characteristics of planktoncommunities at yichang to wuxue section inthe middle reach of the Yangtze River[D].Nanjing:Nanjing Agricultural University.]doi:10.7666/d.Y3185099.

    古小超,王子璐,趙興華,高鍇,羅娜,梅鵬蔚,張震.2023.河流生態(tài)環(huán)境健康評(píng)價(jià)技術(shù)體系構(gòu)建及應(yīng)用[J].中國(guó)環(huán)境監(jiān)測(cè),39(3):87-98.[Gu X C,Wang Z L,Zhao X H,Gao K,Luo N,Mei P W,Zhang Z.2023.Establishment and application of the index system for river ecosystem healthassessment[J].Environmental Monitoring in China,39(3):87-98.]doi:10.19316/j.issn.1002-6002.2023.03.11.

    谷思雨,陳凱,金小偉,李文攀,陳曉飛,熊晶,湯敏喆,姜傳奇,熊杰,李濤,張琪,崔永德,曾宏輝,何舜平,王業(yè)耀,繆煒.2024.水生生物環(huán)境DNA監(jiān)測(cè)技術(shù)的發(fā)展、應(yīng)用與標(biāo)準(zhǔn)化[J].水生生物學(xué)報(bào),48(8):1443-1458.[Gu S Y,Chen K,Jin X W,Li W P,Chen X F,Xiong J,Tang M Z,Jiang C Q,Xiong J,Li T,Zhang Q,Cui Y D,Zeng H H,He S P,Wang Y Y,Miao W.2024.Development,applica-tions,and standardization of environmental DNA monito-ring technology for aquatic organisms[J].Acta Hydrobio-logica Sinica,48(8):1443-1458.]doi:10.7541/2024.2023.0409.

    胡鳳啟,宋平原,李濤.2017.河流特征與河流生態(tài)系統(tǒng)的功能探討[J].黃河水利職業(yè)技術(shù)學(xué)院學(xué)報(bào),29(1):14-16.[Hu F Q,Song PY,Li T.2017.Discussion on river charac-teristics and function of river ecosystem[J].Journal of Yellow River Conservancy Technical Institute,29(1):14-16.]doi:10.13681/j.cnki.cn41-1282/tv.2017.01.004.

    胡愈炘,彭玉,李瑞雯,黃杰,周正,胡圣,王英才,邱光勝.2021.基于環(huán)境DNA宏條形碼的丹江口水庫(kù)浮游生物多樣性及群落特征[J].湖泊科學(xué),33(6):1650-1659.[Hu Y X,Peng Y,Li R W,Huang J,Zhou Z,Hu S,Wang Y C,Qiu G S.2021.Plankton diversity and community charac-teristics in Danjiangkou Reservoir based on environmental DNA metabarcoding[J].Journal of Lake Sciences,33(6):1650-1659.]doi:10.18307/2021.0604.

    紀(jì)瑩璐,王盡文,宿凱,張海舟,趙傳庭,張亮.2023.膠州灣口大型底棲動(dòng)物群落結(jié)構(gòu)及底棲生境健康狀況[J].廣東海洋大學(xué)學(xué)報(bào),43(5):1-10.[Ji Y L,Wang J W,Su K,Zhang H Z,Zhao C T,Zhang L.2023.Community struc-ture of macrobenthos and benthic habitat health assess-ment at the mouth of Jiaozhou Bay[J].Journal of Guang-dong Ocean University,43(5):1-10.]doi:10.3969/j.issn.1673-9159.2023.05.001.

    姜山,張顏,李飛龍,張效偉.2023.基于eDNA的硅藻群落時(shí)空異質(zhì)性及生態(tài)健康評(píng)價(jià)[J].環(huán)境科學(xué),44(1):272-281.[Jiang S,Zhang Y,Li F L,Zhang X W.2023.Spatial and temporal heterogeneity and ecological health evalua-tion of diatom community based on eDNA[J].Environ-mental Science,44(1):272-281.]doi:10.13227/j.hjkx.202204019.

    蔣為,李杰,譚志衛(wèi).2023.基于浮游植物生物完整性指數(shù)的洱海水生態(tài)健康評(píng)價(jià)[J].環(huán)境科學(xué)與技術(shù),46(S1):224-230.[Jiang W,Li J,Tan Z W.2023.The ecosystem health assessment based on phytoplanktonic indexof biotic integ-rity,Erhai Lake[J].Environmental Science and Techno-logy,46(S1):224-230.]doi:10.19672/j.cnki.1003-6504.0434.22.338.

    黎明民,駱鑫,付家想,藍(lán)文陸.2018.基于浮游植物生物完整性的北部灣生態(tài)健康評(píng)價(jià)[J].中國(guó)環(huán)境監(jiān)測(cè),34(6):113-121.[Li M M,Luo X,F(xiàn)u J X,Lan W L.2018.Study on phytoplanktonic index of biotic integrity(P-IBI)assessment for marine ecosystem health:A case of Beibu Gulf[J].Environmental Monitoring in China,34(6):113-121.]doi:10.19316/j.issn.1002-6002.2018.06.12.

    李琦.2023.基于浮游藻類(lèi)生物完整性指數(shù)的伊洛河生態(tài)健康評(píng)價(jià)[D].鄭州:華北水利水電大學(xué).[Li Q.2023.Eco-logical health assessment of Yiluo River based on the bio-logical integrity indexof planktonic algae[D].Zhengzhou:North China University of Water Resources and Electric Power.]doi:10.27144/d.cnki.ghbsc.2023.000870.

    李照,宋書(shū)群,李才文.2017.豐、枯水期長(zhǎng)江口鄰近海域浮游植物群落結(jié)構(gòu)特征及其環(huán)境影響初探[J].海洋學(xué)報(bào),39(10):124-144.[Li Z,Song S Q,Li C W.2017.Prelimi-nary discussion on the phytoplankton assemblages and itsresponse to the environmental changes in the Changjiang(Yangtze)River Estuary and its adjacent waters during the dry season and the wet season[J].Acta Oceanologica Sinica,39(10):124-144.]doi:103969/i.issn.0253-41932017.10011.

    廖延科,張萌,馮兵,張成,夏威,劉足根,田翠翠.2024.基于P-IBI指數(shù)的修河中下游水生態(tài)健康評(píng)價(jià)[J].環(huán)境科學(xué)與技術(shù),47(2):104-115.[Liao Y K,Zhang M,F(xiàn)eng B,Zhang C,Xia W,Liu Z G,Tian C C.2024.Evaluation of aquatic eco-health in the middle and lower reaches of Xiuhe River basin based on the P-IBI index[J].Environ-mental Scienceamp;Technology,47(2):104-115.]doi:10.19672/j.cnki.1003-6504.2127.23.338.

    劉惠秋,李曉東,楊清,巢欣,楊勝嫻,閆冰潔,巴桑.2023.基于浮游植物完整性指數(shù)的雅魯藏布江中上游河流水生態(tài)健康評(píng)價(jià)[J].干旱區(qū)資源與環(huán)境,37(9):109-117.[Liu H Q,Li X D,Yang Q,Chao X,Yang S X,Yan B J,Ba S.2023.Evaluation on water ecological health in middle and upper reaches of the Yarlung Zangbo River based on phytoplankton-index of biotic integrity[J].Jour-nal of Arid Land Resources and Environment,37(9):109-117.]doi:10.13448/j.cnki.jalre.2023.215.

    劉凌,朱燕,李博韜,陳翔.2020.基于MBFG分類(lèi)法的長(zhǎng)江江蘇段浮游植物生物完整性評(píng)價(jià)[J].水資源保護(hù),36(4):13-20.[Liu L,Zhu Y,Li B T,Chen X.2020.Assessment of phytoplankton biological integrity in Jiangsu section of Yangtze River based on MBFG[J].Water Resources Pro-tection,36(4):13-20.]doi:10.3880/j.issn.1004-6933.2020.04.003.

    劉明典,李鵬飛,曾澤國(guó),黃翠,劉紹平.2017.長(zhǎng)江干流安慶段浮游植物群落結(jié)構(gòu)特征[J].淡水漁業(yè),47(4):29-36.[Liu M D,Li P F,Zeng Z G,Huang C,Liu S P.2017.The characteristic of phytoplankton community structure in Anqing section of the Yangtze River[J].Freshwater Fishe-ries,47(4):29-36.]doi:10.3969/j.issn.1000-6907.2017.04.006.

    呂嘉誠(chéng),林淵源,李?lèi)?ài)軍,趙崢.2022.基于環(huán)境DNA宏條形碼的柴河浮游藻類(lèi)研究[J].環(huán)境科學(xué)導(dǎo)刊,41(6):1-7.[LüJ C,Lin Y Y,Li A J,Zhao Z.2022.Study on plankticalgae in Chai River based on environmental DNA macro-barcodes[J].Environmental Science Survey,41(6):1-7.]doi:10.13623/j.cnki.hkdk.2022.06.003.

    馬芊芊.2015.以浮游生物完整性指數(shù)評(píng)價(jià)長(zhǎng)江上游干流宜賓至江津段河流健康度[D].重慶:西南大學(xué).[Ma Q Q.2015.Assessment of ecosystem health of upper Yangtze River in Yibin to Jiangjin using plankton-index of biotic integrity[D].Chongqing:Southwest University.]

    錢(qián)爍統(tǒng),劉凌,陳寧,燕文明.2018.河生藻類(lèi)形態(tài)功能群分布及其與環(huán)境因子關(guān)系的分析[J].水電能源科學(xué),36(6):146-149.[Qian S T,Liu L,Chen N,Yan W M.2018.Distribution of fiuvial planktonic algae based on functional groups and its correlations with environmental factors[J].Water Resources and Power,36(6):146-149.]

    邱陽(yáng)凌,林育青,劉俊杰,唐磊,關(guān)鐵生,陳求穩(wěn),陳凱,王麗.2018.淮河干流及主要支流夏季浮游植物群落生物多樣性評(píng)價(jià)[J].環(huán)境科學(xué)學(xué)報(bào),38(4):1665-1672.[Qiu Y L,Lin Y Q,Liu J J,Yang L,Guan T S,Chen Q W,Chen K,Wang L.2018.The biodiversity assessment of phytoplank-ton community in summer within main stream and tribu-taryof Huaihe River[J].Acta Scientiae Circumstantiae,38(4):1665-1672.]doi:10.13671/j.hjkxxb.2017.0408.

    邵海燕,王卿,高春霞,趙靜.2024.長(zhǎng)江口浮游植物群落特征及影響因素分析[J].大連海洋大學(xué)學(xué)報(bào),39(1):124-133.[Shao H Y,Wang Q,Gao C X,Zhao J.2024.Analy-sis of phytoplankton community characteristics and influen-cing factors in the Yangtze River Estuary[J].Journal of Dalian Fisheries University,39(1):124-133.]doi:10.16535/j.cnki.dlhyxb.2023-171.

    沈強(qiáng),俞建軍,陳暉,胡菊香,陳明秀,鄭金秀,李聃,張俊芳.2012.浮游生物完整性指數(shù)在浙江水源地水質(zhì)評(píng)價(jià)中的應(yīng)用[J].水生態(tài)學(xué)雜志,33(2):26-31.[Shen Q,Yu J J,Chen H,Hu J X,Chen M X,Zheng J X,Li D,Zhang J F.2012.Application of plankton integrity index to water qua-lity assessment of water source in Zhejiang Province[J].Journal of Hydroecology,33(2):26-31.]doi:10.15928/j.1674-3075.2012.02.013.

    盛祥銳,羅遵蘭,孫光,肖能文,康玲玲,李慧清,韓靜.2024.基于浮游植物的生物完整性指數(shù)開(kāi)發(fā)與驗(yàn)證——以永定河北京段為例[J].環(huán)境科學(xué)學(xué)報(bào),44(2):489-500.[Sheng X R,Luo Z L,Sun G,Xiao N W,Kang L L,Li H Q,Han J.2024.Development and verification of the phyto-plankton biological integrity index—A case study in Yong-ding River in Beijing[J].Acta Scientiae Circumstantiae,44(2):489-500.]doi:10.13671/j.hjkxxb.2023.0295.

    譚巧,劉凌,黃艷芬,張順婷.2023.長(zhǎng)江江蘇段著生藻類(lèi)群落結(jié)構(gòu)與生物完整性分析[J].人民長(zhǎng)江,54(12):74-81.[Tan Q,Liu L,Huang Y F,Zhang S T.2023.Analysis ofcommunity structure and biotic integrity of periphyton inJiangsu section of Changjiang River[J].Yangtze River,54(12):74-81.]doi:10.16232/j.cnki.1001-4179.2023.12.010.譚巧,馬芊芊,李斌斌,呂紅健,付梅,姚維志.2017.應(yīng)用浮游

    植物生物完整性指數(shù)評(píng)價(jià)長(zhǎng)江上游河流健康[J].淡水漁業(yè),47(3):97-104.[Tan Q,Ma Q Q,Li B B,LüH J,F(xiàn)u M,Yao W Z.2017.Ecological health assessment of the upper reaches of the Yangtze River,based on biotic integ-rity index of phytoplankton[J].Freshwater Fisheries,47(3):97-104.]doi:10.3969/j.issn.1000-6907.2017.03.015.

    田豐歌,何靜,王符菁,黃穎華,莊澤彬,黃彬彬.2023.夏季珠江口兩種網(wǎng)型網(wǎng)采浮游動(dòng)物的群落特征比較[J].廣東海洋大學(xué)學(xué)報(bào),43(2):51-59.[Tian F G,He J,Wang F J,Huang Y H,Zhuang Z B,Huang B B.2023.Comparison of community characteristics of zooplankton sampled by two mesh nets in the Pearl River Estuary in summer[J].Journal of Guangdong Ocean University,43(2):51-59.]doi:10.3969/j.issn.1673-9159.2023.02.007.

    王備新,楊蓮,劉正文.2006.生物完整性指數(shù)與水生態(tài)系統(tǒng)健康評(píng)價(jià)[J].生態(tài)學(xué)雜志,25(6):707-710.[Wang B X,Yang L,Liu Z W.2006.Bindex of biological integrity andits application in health assessment of aquatic ecosystem[J].Chinese Journal of Ecology,25(6):707-710.]

    王晨,陶孟,李?lèi)?ài)民,施鵬,楊江華,王志浩,張效偉.2022.基于環(huán)境DNA宏條形碼技術(shù)的秦淮河生物多樣性研究[J].生態(tài)學(xué)報(bào),42(2):611-624.[Wang C,Tao M,LiAM,Shi P,Yang J H,Wang Z H,Zhang X W.2022.Research on the biodiversity of Qinhuai River based on environmen-tal DNAmetabacroding[J].Acta Ecologica Sinica,42(2):611-624.]doi:10.5846/stxb202007221920.

    王迪,梁皓瑞,張軍曉,陳際雨,梁宇釗,曾芳,馬媛,張際標(biāo),張鵬,何嘉鏵,張偉杰.2024.湛江灣浮游植物優(yōu)勢(shì)種空間生態(tài)位[J].廣東海洋大學(xué)學(xué)報(bào),44(5):40-49.[Wang D,Liang H R,Zhang J X,Chen J Y,Liang Y Z,Zeng F,Ma Y,Zhang J B,Zhang P,He J H,Zhang W J.2024.Spa-tial ecological niche of dominant phytoplankton species in Zhan-jiang Bay[J].Journal of Guangdong Ocean Univer-sity,44(5):40-49.]doi:10.3969/j.issn.1673-9159.2024.05.005.

    王芳,李永吉,馬廷婷,陳輝輝,王曉龍,李寬意,吳召仕.2022.基于浮游植物的城市湖泊生態(tài)健康評(píng)價(jià)—以長(zhǎng)江下游銅陵市西湖為例[J].湖泊科學(xué),34(6):1890-1900.[Wang F,Li Y J,Ma T T,Chen H H,Wang X L,Li K Y,Wu Z S.2022.Ecological health assessment of urban lake based on phytoplankton—A case study of Lake Xihu,Tongling,lower reaches of the Yangtze River[J].Journal of Lake Science,34(6):1890-1900.]doi:10.18307/2022.0628.

    王汝賢,楊剛,耿智,趙峰,馮雪松,張濤.2023.環(huán)境DNA技術(shù)在長(zhǎng)江口魚(yú)類(lèi)多樣性分析中的應(yīng)用[J].水生生物學(xué)報(bào),47(3):365-375.[Wang R X,Yang G,Geng Z,Zhao F,F(xiàn)eng X S,Zhang T.2023.Application of environmental DNA technology in fish diversity analysis in the Yangtze River Estuary[J].Acta Hydrobiologica Sinica,47(3):365-375.]doi:10.7541/2023.2022.0272.

    魏念,李云峰,吳凡,茹輝軍,吳湘香,楊傳順,倪朝輝.2023.長(zhǎng)江流域重點(diǎn)水域浮游植物群落結(jié)構(gòu)特征及影響因素[J].長(zhǎng)江流域資源與環(huán)境,32(12):2519-2527.[Wei N,Li Y F,Wu F,Ru H J,Wu X X,Yang C S,Ni Z H.2023.Phytoplankton community structures and impact factors in key water areas of Yangtze River Basin[J].Resources and Environment in the Yangtze Basin,32(12):2519-2527.]doi:10.11870/cjlyzyyhj202312006.

    吳泳江,蔡明成,孫翰昌,胡慧蝶,鄧雅心,龍應(yīng)根.2023.環(huán)境DNA宏條形碼在浮游動(dòng)植物多樣性研究中的應(yīng)用[J].漁業(yè)研究,45(4):399-407.[Wu Y J,Cai M C,Sun H C,Hu H D,Deng Y X,Long Y G.2023.Application of envi-ronmental DNA metabarcoding in studying biodiversity of zooplankton and phytoplankton[J].Journal of Fisheries Research,45(4):399-407.]doi:10.14012/j.cnki.fjsc.2023.04.010.

    肖善勢(shì),郝雅賓,劉金殿,張愛(ài)菊,王俊,羅偉,何海生,周志明.2021.應(yīng)用生物完整性指數(shù)評(píng)價(jià)錢(qián)塘江流域—浙江段水生態(tài)系統(tǒng)健康[J].水產(chǎn)科學(xué),40(5):740-749.[Xiao S S,Hao Y B,Liu J D,Zhang A J,Wang J,Luo W,He H S,Zhou Z M.2021.Evaluation of water ecosystem health in Qiantang River Basin and Zhejiang section by biologi-cal integrity index[J].Fisheries Sciences,40(5):740-749.]doi:10.16378/j.cnki.1003-1111.20031.

    謝滿(mǎn)華,任青松.2016.長(zhǎng)江蕪湖段水域浮游生物群落年度結(jié)構(gòu)變化研究[J].現(xiàn)代農(nóng)業(yè)科技,(11):239-241.[Xie M H,Ren Q S.2016.Annua1 changes of plankton commu-nity structure in Wuhu mainstream of Yangtze River[J].Modern Agricultural Science and Technology,(11):239-241.]doi:10.3969/j.issn.1007-5739.2016.11.150.

    熊梅君,李秋華,陳倩,何應(yīng),馬欣洋,韓孟書(shū),唐重麗,張華俊.2019.基于GIS應(yīng)用P-IBI評(píng)價(jià)貴州高原百花水庫(kù)水生態(tài)系統(tǒng)健康[J].生態(tài)學(xué)雜志,38(10):3093-3101.[Xiong M J,Li Q H,Chen Q,He Y,Ma X Y,Han M S,Tang Z L,Zhang H J.2019.Evaluation of water ecosystem health of Baihua Reservoir in Guizhou based on GIS and P-IBI[J].Chinese Journal of Ecology,38(10):3093-3101.]doi:10.13292/j.1000-4890.201910.010.

    徐宗學(xué),馬欣洋.2023.河流生態(tài)系統(tǒng)健康評(píng)價(jià)——以銀川市典農(nóng)河為例[J].水利發(fā)展研究,23(9):1-7.[Xu Z X,Ma X Y.2023.[Assessment of river ecosystem health:A case study of Diannong River in Yinchuan City[J].Water Re-sources Development Research,23(9):1-7.]doi:10.13928/j.cnki.wrdr.2023.09.001.]

    楊桂山,徐昔保.2020.長(zhǎng)江經(jīng)濟(jì)帶“共抓大保護(hù),不搞大開(kāi)發(fā)”的基礎(chǔ)與策略[J].中國(guó)科學(xué)院院刊,35(8):940-950.[Yang G S,Xu X B.2020.Foundation and strategy of well-coordinated environmental conservation and avoiding excessive development in the Yangtze River economic belt[J].Bulletin of Chinese Academy of Sciences,35(8):940-950.]doi:10.16418/j.issn.1000-3045.20200511002.

    余景,趙漫,胡啟偉,陳丕茂.2017.基于魚(yú)類(lèi)生物完整性指數(shù)的深圳鵝公灣漁業(yè)水域健康評(píng)價(jià)[J].南方農(nóng)業(yè)學(xué)報(bào),48(3):524-531.[Yu J,Zhao M,Hu Q W,Chen P M.2017.Ecosystem health assessment of Egong Bay fisheries areain Shenzhen based on index of biological integrity of fish[J].Journal of Southern Agriculture,48(3):524-531.]doi:10.3969/j:issn.2095-1191.2017.03.024.

    張靜,胡愈炘,胡圣,黃杰.2023.長(zhǎng)江流域浮游植物群落的環(huán)境驅(qū)動(dòng)及生態(tài)評(píng)價(jià)[J].環(huán)境科學(xué),44(4):2072-2082.[Zhang J,Hu Y X,Hu S,Huang J.2023.Environmental driving factors and assessment on the aquatic ecosystem of periphytic algae of six inflow rivers in Yangtze River basin[J].Environmental Science,44(4):2072-2082.]doi:10.13227/j.hjkx.202206095.

    張瀟月.2018.長(zhǎng)江下游干流浮游植物多樣性研究[D].上海:上海師范大學(xué).[Zhang XY.2018.A study on the phy-toplankton diversity in the lower reaches of the Yangtze River[D].Shanghai:Shanghai Normal University.]doi:10.7666/d.D01425867.

    鄭丙輝,田自強(qiáng),張雷,鄭凡東.2007.太湖西岸湖濱帶水生生物分布特征及水質(zhì)營(yíng)養(yǎng)狀況[J].生態(tài)學(xué)報(bào),(10):4214-4223.[Zheng B H,Tian Z Q,Zhang L,Zheng F D.2007.The characteristics of the hydrobios'distribution and the analysis of water quality along the West Shore of Taihu Lake[J].Acta Ecologica Sinica,(10):4214-4223.]doi:10.3321/j.issn:1000-0933.2007.10.032.

    鐘可兒,楊瀟,馬吉順,周瓊.2022.長(zhǎng)江干流浮游植物群落結(jié)構(gòu)空間異質(zhì)性及其影響因素[J].長(zhǎng)江流域資源與環(huán)境,31(11):2458-2472.[Zhong K E,Yang X,Ma J S,Zhou Q.2022.Spatial Heterogeneity and driving factors of phy-toplankton community structure in the main stem of the Yangtze River[J].Resources and Environment in the Yang-tze Basin,31(11):2458-2472.]doi:10.11870/cjlyzyyhj 202211011.

    DeSantis T Z,Hugenholtz P,Larsen N,Rojas M,Brodie E L,Keller K,Huber T,Dalevi D,Hu P,Andersen G L.2006.Greengenes,a chimera-checked 16S rRNA gene database and workbench compatible with ARB[J].Applied and Environmental Microbiology,72(7):5069-5072.doi:10.1128/AEM.03006-05.

    Guo W,Ye F,Xu S D,Jia G D.2015.Seasonal variation in sources and processing of particulate organic carbon in the Pearl River Estuary,South China[J].Estuarine,Coastal Shelf Science,167:540-548.doi:10.1016/j.ecss.2015.11.004.

    Quast C,Pruesse E,Yilmaz P,Gerken J,Schweer T,Yarza P,Peplies J,Gl?ckner F O.2013.The SILVA ribosomal RNA gene database project:Improved data processing and web-based tools[J].Nucleic Acids Research,41(Dl):D590-D596.doi:10.1093/nar/gks 1219.

    Stamou G,Katsiapi M,Moustaka-Gouni M,Michaloudi E.2019.Trophic state assessment based on zooplankton com-munities in Mediterranean lakes[J].Hydrobiologia,844:83-103.doi:10.1007/s 10750-018-3880-9.

    Turak E,Harrison I,Dudgeon D,Abell R,Bush A,Darwall W,F(xiàn)inlayson M,F(xiàn)errie S,F(xiàn)reyhof J,Hermoso V,Juffe-Bignoli D,Linke Simon,Nel J,C.Patricio H,Pittock J,Raghavan R,Revenga C,P.Simaika J,De Wever A.2017.Essential biodiversity variables for measuring change in global freshwater biodiversity[J].Biological Conserva-tion,213:272-279.doi:10.1016/j.biocon.2016.09.005.

    (責(zé)任編輯:蘭宗寶,張博)

    女生性感内裤真人,穿戴方法视频| 国产淫片久久久久久久久| 欧美高清性xxxxhd video| 欧美日韩国产亚洲二区| 悠悠久久av| 99久久精品热视频| 久久久精品欧美日韩精品| 午夜福利在线在线| 最近视频中文字幕2019在线8| 中文字幕av在线有码专区| 国产成人a区在线观看| 中出人妻视频一区二区| 精品久久久久久成人av| 国产精品一区二区免费欧美| 91麻豆精品激情在线观看国产| 精品一区二区三区人妻视频| 久久久久国产精品人妻aⅴ院| 天堂av国产一区二区熟女人妻| 成人精品一区二区免费| 色综合婷婷激情| 亚洲欧美日韩东京热| 网址你懂的国产日韩在线| 97人妻精品一区二区三区麻豆| 我要看日韩黄色一级片| 日本撒尿小便嘘嘘汇集6| 成人精品一区二区免费| 色在线成人网| 中文字幕熟女人妻在线| 免费看日本二区| 色5月婷婷丁香| 精品人妻1区二区| 欧美性猛交黑人性爽| 国产精品一区www在线观看 | 亚洲av中文av极速乱 | 99热这里只有是精品在线观看| 国产aⅴ精品一区二区三区波| 99热这里只有是精品50| 国产av麻豆久久久久久久| 国产探花极品一区二区| 成人国产综合亚洲| 小说图片视频综合网站| 身体一侧抽搐| 国产精品久久久久久亚洲av鲁大| 一区二区三区高清视频在线| 无遮挡黄片免费观看| 国产欧美日韩精品亚洲av| 99久久精品国产国产毛片| 老女人水多毛片| 日韩中文字幕欧美一区二区| 亚洲av免费在线观看| 一本久久中文字幕| 老女人水多毛片| 亚洲国产高清在线一区二区三| 91麻豆精品激情在线观看国产| 午夜老司机福利剧场| 久久精品国产鲁丝片午夜精品 | 男人和女人高潮做爰伦理| 亚洲天堂国产精品一区在线| 久久久久久九九精品二区国产| 欧美精品国产亚洲| 午夜激情欧美在线| 国产精品亚洲美女久久久| 亚洲乱码一区二区免费版| a级毛片免费高清观看在线播放| 国产成人av教育| 波多野结衣巨乳人妻| 日韩一本色道免费dvd| 琪琪午夜伦伦电影理论片6080| 能在线免费观看的黄片| 搞女人的毛片| 精品久久久久久久久亚洲 | 三级国产精品欧美在线观看| 狂野欧美激情性xxxx在线观看| 99热6这里只有精品| 亚洲中文日韩欧美视频| 日韩欧美精品v在线| 观看美女的网站| 亚洲av电影不卡..在线观看| 成人特级黄色片久久久久久久| 97超视频在线观看视频| 国产成人av教育| 亚洲在线观看片| 免费av不卡在线播放| 亚洲精品一区av在线观看| 国产欧美日韩精品亚洲av| 精品免费久久久久久久清纯| 18禁裸乳无遮挡免费网站照片| 精品无人区乱码1区二区| 老司机午夜福利在线观看视频| 国产成人a区在线观看| 免费看光身美女| 我的老师免费观看完整版| 午夜福利18| av视频在线观看入口| 内地一区二区视频在线| 波野结衣二区三区在线| 天堂√8在线中文| 好男人在线观看高清免费视频| 欧美潮喷喷水| 亚洲欧美日韩无卡精品| 亚洲精华国产精华精| 中文在线观看免费www的网站| 免费观看在线日韩| 91久久精品国产一区二区三区| 韩国av在线不卡| 91狼人影院| 国产真实伦视频高清在线观看 | 欧美区成人在线视频| 国产私拍福利视频在线观看| 国产麻豆成人av免费视频| 亚洲成av人片在线播放无| 日本色播在线视频| 午夜福利视频1000在线观看| 日韩在线高清观看一区二区三区 | 麻豆国产av国片精品| 国产精品野战在线观看| 国产一区二区在线观看日韩| 久久久色成人| 欧美高清成人免费视频www| 无遮挡黄片免费观看| 久久精品久久久久久噜噜老黄 | 小蜜桃在线观看免费完整版高清| 亚洲精品成人久久久久久| 免费搜索国产男女视频| 久久国产精品人妻蜜桃| 亚洲成人精品中文字幕电影| 两个人视频免费观看高清| 在线播放国产精品三级| 国产精品综合久久久久久久免费| 最近视频中文字幕2019在线8| 91麻豆av在线| 亚洲一区二区三区色噜噜| 国产精品日韩av在线免费观看| 国产视频内射| 日韩中文字幕欧美一区二区| av在线天堂中文字幕| 91麻豆av在线| 国产女主播在线喷水免费视频网站 | 亚洲经典国产精华液单| 三级国产精品欧美在线观看| 99精品在免费线老司机午夜| 亚洲国产欧洲综合997久久,| 欧美黑人巨大hd| 成人二区视频| 少妇的逼水好多| 亚洲一区二区三区色噜噜| 欧美黑人欧美精品刺激| 亚洲欧美日韩高清专用| 十八禁国产超污无遮挡网站| 999久久久精品免费观看国产| 香蕉av资源在线| 欧美在线一区亚洲| 美女高潮的动态| 中文字幕人妻熟人妻熟丝袜美| 日韩在线高清观看一区二区三区 | 无遮挡黄片免费观看| 久久久久久国产a免费观看| 久久久久国内视频| 日本在线视频免费播放| 我的老师免费观看完整版| 偷拍熟女少妇极品色| 国内久久婷婷六月综合欲色啪| 麻豆久久精品国产亚洲av| 91久久精品国产一区二区三区| 特大巨黑吊av在线直播| 亚洲国产高清在线一区二区三| a级毛片a级免费在线| 老司机福利观看| 中文亚洲av片在线观看爽| 国产精品三级大全| 高清在线国产一区| 精品久久久久久,| 色精品久久人妻99蜜桃| av在线亚洲专区| 久久精品国产亚洲av香蕉五月| 国产真实乱freesex| 日韩国内少妇激情av| 日日撸夜夜添| 国产v大片淫在线免费观看| 国产亚洲精品久久久久久毛片| 国产精品三级大全| 亚洲天堂国产精品一区在线| 国产一区二区三区av在线 | 国产黄色小视频在线观看| 国内久久婷婷六月综合欲色啪| 成人毛片a级毛片在线播放| 少妇裸体淫交视频免费看高清| 级片在线观看| 五月伊人婷婷丁香| 亚洲人与动物交配视频| 国产视频一区二区在线看| 国产一区二区亚洲精品在线观看| 直男gayav资源| 国产乱人视频| 禁无遮挡网站| 成人国产综合亚洲| 极品教师在线免费播放| 啦啦啦韩国在线观看视频| 亚洲精品日韩av片在线观看| 桃红色精品国产亚洲av| 少妇熟女aⅴ在线视频| 亚洲熟妇中文字幕五十中出| 色尼玛亚洲综合影院| 日韩一本色道免费dvd| 毛片一级片免费看久久久久 | 亚洲av免费在线观看| 欧美丝袜亚洲另类 | 熟女人妻精品中文字幕| 小说图片视频综合网站| 嫩草影院新地址| 淫妇啪啪啪对白视频| 精品久久久久久,| 禁无遮挡网站| 黄色日韩在线| 69人妻影院| 成人综合一区亚洲| 午夜a级毛片| 人妻少妇偷人精品九色| 久久香蕉精品热| 一级毛片久久久久久久久女| 99热这里只有是精品50| 久久精品久久久久久噜噜老黄 | 亚洲成人中文字幕在线播放| 午夜激情欧美在线| 两个人的视频大全免费| www.www免费av| 成人精品一区二区免费| 亚洲综合色惰| 99久国产av精品| 久久久久久九九精品二区国产| 久久亚洲真实| 欧美+日韩+精品| 国产精品人妻久久久影院| 成年女人毛片免费观看观看9| 精品一区二区三区人妻视频| 日本 欧美在线| 欧美高清成人免费视频www| 高清日韩中文字幕在线| 中出人妻视频一区二区| 人人妻人人澡欧美一区二区| 能在线免费观看的黄片| 国产视频一区二区在线看| 一级毛片久久久久久久久女| 久久久久久久久中文| 国产精品久久久久久精品电影| 精品无人区乱码1区二区| 全区人妻精品视频| 精品久久久久久久末码| 国产精品一区二区三区四区免费观看 | 亚洲熟妇熟女久久| 日本一本二区三区精品| 黄色女人牲交| 国产日本99.免费观看| 在线天堂最新版资源| 日本免费a在线| 大又大粗又爽又黄少妇毛片口| 国内揄拍国产精品人妻在线| 国产v大片淫在线免费观看| 国产爱豆传媒在线观看| 最近在线观看免费完整版| 女同久久另类99精品国产91| 国产黄色小视频在线观看| 成人欧美大片| 国产男人的电影天堂91| 又黄又爽又刺激的免费视频.| 精品久久久久久久久亚洲 | 最新中文字幕久久久久| 蜜桃亚洲精品一区二区三区| 日本五十路高清| 亚洲最大成人手机在线| 全区人妻精品视频| 少妇猛男粗大的猛烈进出视频 | 久久精品国产鲁丝片午夜精品 | 久久国产精品人妻蜜桃| 亚洲综合色惰| 国产真实伦视频高清在线观看 | 国产黄片美女视频| 我要搜黄色片| 色5月婷婷丁香| 有码 亚洲区| 中文字幕人妻熟人妻熟丝袜美| 美女大奶头视频| 热99re8久久精品国产| 动漫黄色视频在线观看| 一区二区三区免费毛片| 婷婷精品国产亚洲av| 欧美性感艳星| 久久久久免费精品人妻一区二区| 欧美日韩中文字幕国产精品一区二区三区| 免费观看精品视频网站| 十八禁网站免费在线| 99久久成人亚洲精品观看| 91av网一区二区| 一级黄片播放器| 午夜福利在线观看免费完整高清在 | 国产午夜精品久久久久久一区二区三区 | 日本爱情动作片www.在线观看 | 国产男靠女视频免费网站| 国产淫片久久久久久久久| 99视频精品全部免费 在线| 久久6这里有精品| 欧美日本视频| 亚洲av免费在线观看| 别揉我奶头 嗯啊视频| 亚洲欧美日韩无卡精品| 桃色一区二区三区在线观看| 国产三级在线视频| 夜夜夜夜夜久久久久| 日本免费一区二区三区高清不卡| 简卡轻食公司| 日本黄色视频三级网站网址| 高清毛片免费观看视频网站| 女同久久另类99精品国产91| 欧美黑人巨大hd| 亚洲最大成人中文| 两个人视频免费观看高清| 中国美白少妇内射xxxbb| 亚洲精品粉嫩美女一区| 99久久精品一区二区三区| 国产成年人精品一区二区| 精品久久久久久久久久久久久| 婷婷色综合大香蕉| 老司机午夜福利在线观看视频| 亚洲精品一区av在线观看| 国产成人av教育| 日韩欧美 国产精品| a在线观看视频网站| 亚洲av五月六月丁香网| 嫩草影院入口| av天堂中文字幕网| 亚洲精品久久国产高清桃花| 国产老妇女一区| 亚洲精华国产精华精| 成人国产综合亚洲| 天天躁日日操中文字幕| 麻豆成人av在线观看| 日韩av在线大香蕉| 欧美不卡视频在线免费观看| 久久亚洲真实| 国产精品98久久久久久宅男小说| 欧美激情在线99| 亚洲欧美日韩东京热| 午夜亚洲福利在线播放| 国内精品一区二区在线观看| 99国产精品一区二区蜜桃av| 成年版毛片免费区| 最近最新免费中文字幕在线| 精品不卡国产一区二区三区| 亚洲五月天丁香| 亚洲av五月六月丁香网| 亚洲精华国产精华精| 国产乱人伦免费视频| 欧美在线一区亚洲| av福利片在线观看| 亚洲欧美清纯卡通| 动漫黄色视频在线观看| 久久久国产成人免费| av在线蜜桃| 好男人在线观看高清免费视频| 精品一区二区三区视频在线观看免费| 亚洲欧美日韩高清专用| ponron亚洲| 国产精品嫩草影院av在线观看 | 少妇猛男粗大的猛烈进出视频 | av在线亚洲专区| 欧美性感艳星| 日韩,欧美,国产一区二区三区 | 99视频精品全部免费 在线| 国模一区二区三区四区视频| 亚洲精品一卡2卡三卡4卡5卡| 性插视频无遮挡在线免费观看| 国产乱人伦免费视频| 欧美黑人巨大hd| 亚洲av二区三区四区| 色精品久久人妻99蜜桃| 免费观看在线日韩| 人妻少妇偷人精品九色| 亚洲精品在线观看二区| 成人二区视频| 男女下面进入的视频免费午夜| 亚洲人成网站在线播| avwww免费| 女的被弄到高潮叫床怎么办 | 亚洲在线自拍视频| 男女边吃奶边做爰视频| 日韩欧美在线乱码| 男女之事视频高清在线观看| 看免费成人av毛片| 欧美成人一区二区免费高清观看| 成人特级黄色片久久久久久久| 春色校园在线视频观看| 国产精品一区二区性色av| 麻豆国产av国片精品| 一进一出抽搐gif免费好疼| av女优亚洲男人天堂| 人妻夜夜爽99麻豆av| 午夜爱爱视频在线播放| 久久欧美精品欧美久久欧美| 国产精品伦人一区二区| 国内揄拍国产精品人妻在线| 国产大屁股一区二区在线视频| 亚洲精品色激情综合| 国内精品久久久久精免费| 少妇的逼好多水| 成人欧美大片| 深爱激情五月婷婷| 亚洲国产精品合色在线| 成人无遮挡网站| 亚洲第一区二区三区不卡| 国内精品宾馆在线| 国产三级中文精品| 免费看av在线观看网站| 久久精品夜夜夜夜夜久久蜜豆| 久久久久精品国产欧美久久久| 在线免费观看的www视频| 十八禁国产超污无遮挡网站| 精品久久国产蜜桃| 成人永久免费在线观看视频| 日本爱情动作片www.在线观看 | 日韩欧美一区二区三区在线观看| 天堂网av新在线| 久久精品综合一区二区三区| 国产免费男女视频| 久久午夜福利片| 精品一区二区三区视频在线观看免费| 51国产日韩欧美| 亚洲美女黄片视频| 久久久久久久久久久丰满 | 欧美极品一区二区三区四区| www.色视频.com| 国产 一区 欧美 日韩| 亚洲18禁久久av| 午夜视频国产福利| 日本爱情动作片www.在线观看 | www.色视频.com| 精品久久久久久久久久免费视频| a级毛片a级免费在线| 国产精品99久久久久久久久| 中亚洲国语对白在线视频| 成人特级黄色片久久久久久久| 国产真实乱freesex| 午夜免费成人在线视频| 欧美另类亚洲清纯唯美| 美女高潮喷水抽搐中文字幕| 日韩欧美在线乱码| 久久久久久久精品吃奶| 日本黄色视频三级网站网址| 精品无人区乱码1区二区| 一级a爱片免费观看的视频| 亚洲不卡免费看| 亚洲一级一片aⅴ在线观看| 最近最新中文字幕大全电影3| 99在线人妻在线中文字幕| 老女人水多毛片| 久久这里只有精品中国| 久久久久久久亚洲中文字幕| 露出奶头的视频| 日韩国内少妇激情av| 欧美xxxx黑人xx丫x性爽| 免费观看在线日韩| 韩国av在线不卡| 国产一区二区在线观看日韩| 精品久久久久久久久久免费视频| 久久久精品欧美日韩精品| 黄色配什么色好看| 禁无遮挡网站| 婷婷精品国产亚洲av在线| 精品日产1卡2卡| 久久久久久久久大av| 国产av麻豆久久久久久久| 免费看美女性在线毛片视频| 三级男女做爰猛烈吃奶摸视频| 久久午夜亚洲精品久久| 国产精品福利在线免费观看| 别揉我奶头 嗯啊视频| 久久精品影院6| netflix在线观看网站| 日本a在线网址| 国产精品亚洲一级av第二区| 精品人妻熟女av久视频| 国产一区二区在线观看日韩| 午夜精品一区二区三区免费看| 亚洲狠狠婷婷综合久久图片| 白带黄色成豆腐渣| 午夜精品久久久久久毛片777| 国产精品福利在线免费观看| 级片在线观看| 露出奶头的视频| 伦理电影大哥的女人| 在现免费观看毛片| 少妇人妻一区二区三区视频| 十八禁国产超污无遮挡网站| 欧美性猛交黑人性爽| 午夜日韩欧美国产| 亚洲无线在线观看| 亚洲在线自拍视频| 真人做人爱边吃奶动态| 一区二区三区高清视频在线| 国产精品美女特级片免费视频播放器| 神马国产精品三级电影在线观看| 国产中年淑女户外野战色| 男女啪啪激烈高潮av片| 可以在线观看的亚洲视频| 国产精品久久久久久亚洲av鲁大| 久久精品国产亚洲网站| 日韩欧美国产一区二区入口| 亚洲成人久久爱视频| 欧美+亚洲+日韩+国产| 美女 人体艺术 gogo| 欧美+亚洲+日韩+国产| 久久久久性生活片| 最好的美女福利视频网| 国产精品美女特级片免费视频播放器| 亚洲成人免费电影在线观看| 毛片女人毛片| 看免费成人av毛片| 免费观看人在逋| 国产午夜精品论理片| 丰满乱子伦码专区| 真人做人爱边吃奶动态| 精品久久久久久久久av| 欧美中文日本在线观看视频| 日韩欧美国产一区二区入口| 美女高潮的动态| 草草在线视频免费看| 日本一本二区三区精品| 黄片wwwwww| 欧美另类亚洲清纯唯美| 亚洲精品成人久久久久久| 尤物成人国产欧美一区二区三区| 黄色女人牲交| 又黄又爽又刺激的免费视频.| 国产一区二区亚洲精品在线观看| 免费av观看视频| 国产一区二区三区在线臀色熟女| 人人妻人人看人人澡| a级毛片免费高清观看在线播放| 欧美日本视频| 亚洲欧美日韩卡通动漫| 看十八女毛片水多多多| 嫁个100分男人电影在线观看| 亚洲中文字幕日韩| 欧美成人性av电影在线观看| 99视频精品全部免费 在线| 狂野欧美激情性xxxx在线观看| 淫妇啪啪啪对白视频| 国产黄a三级三级三级人| 亚洲午夜理论影院| 69av精品久久久久久| av.在线天堂| 国产人妻一区二区三区在| av中文乱码字幕在线| 国产69精品久久久久777片| 色在线成人网| 国产v大片淫在线免费观看| 免费电影在线观看免费观看| 国产伦精品一区二区三区四那| 久久精品综合一区二区三区| 不卡视频在线观看欧美| 毛片一级片免费看久久久久 | 国产精品av视频在线免费观看| 亚洲美女搞黄在线观看 | 亚洲精品国产成人久久av| 人妻制服诱惑在线中文字幕| 亚洲精品久久国产高清桃花| 成人毛片a级毛片在线播放| 99热6这里只有精品| 日韩中字成人| 国产 一区 欧美 日韩| 少妇的逼好多水| 国产免费男女视频| 午夜免费激情av| 成人av一区二区三区在线看| 午夜激情福利司机影院| 久久精品久久久久久噜噜老黄 | 国产一区二区三区视频了| 国产成人av教育| 国产 一区精品| 国产av一区在线观看免费| 国产一区二区在线av高清观看| 亚洲久久久久久中文字幕| 精品久久久久久久久亚洲 | 国产午夜精品久久久久久一区二区三区 | 亚洲精品成人久久久久久| 欧美成人免费av一区二区三区| 日韩欧美三级三区| 九色国产91popny在线| 欧美日韩中文字幕国产精品一区二区三区| 欧美激情久久久久久爽电影| 欧美黑人巨大hd| 三级国产精品欧美在线观看| 欧美色欧美亚洲另类二区| 久久精品夜夜夜夜夜久久蜜豆| 国产aⅴ精品一区二区三区波| 国产精品爽爽va在线观看网站| 免费无遮挡裸体视频| av女优亚洲男人天堂| 精品久久久久久久久亚洲 | 有码 亚洲区| 亚洲性久久影院| 午夜福利在线观看吧| 一区二区三区高清视频在线| 国产精品久久久久久久久免| 国产精品久久久久久精品电影| 赤兔流量卡办理| 免费搜索国产男女视频| 亚洲国产高清在线一区二区三| 狂野欧美白嫩少妇大欣赏| av国产免费在线观看| 日本一本二区三区精品| 亚洲av一区综合| 一区二区三区激情视频| 亚洲四区av| 精品人妻1区二区| 亚洲精品日韩av片在线观看|