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

    基于環(huán)境DNA宏條形碼技術(shù)的大洋河下游魚類多樣性研究

    2025-08-03 00:00:00朱春月楊培民劉忠航張伯序胡宗云張健
    安徽農(nóng)業(yè)科學(xué) 2025年12期
    關(guān)鍵詞:魚類位點(diǎn)物種

    中圖分類號(hào)S932. 文獻(xiàn)標(biāo)識(shí)碼A文章編號(hào) 0517-6611(2025)12-0059-06

    doi:10.3969/j.issn.0517-6611.2025.12.015

    開放科學(xué)(資源服務(wù))標(biāo)識(shí)碼(OSID):

    Study on the Fish Diversity in the Downstream of Dayang River Based on Environmental DNA Metabarcoding Technology ZHUChun-yue,YANGPei-min,LIUZhong-hangetal(LiaoningItituteofFreshwaterFisheries/LiaoningKeyLaboratoryofquatic Animal Diseases Control,Liaoyang,Liaoning111000)

    AbstractInodetounderstandthspecicompositooftifisspecisommunitisandprotectheirdiversityintedowsteamof DayangRiverviroetalDmetabarodinehgasusdtoaletespisdvesityofsiedostreaoagi er.Throughthecolletionandhigh-troughputsequencingaalysisofDA(eDNA)samplesfromthedownstreamenvironentoftheDayang River,atotalofsiesofsectediplsfrotreaftaoitga,, 9orders.Amongtherders,Cypriniforeshdthelgestspeiesumber34speies)folowdbyGobfos(1Ospecies)llfs speciesHophalmchorinlematocesdilssdgheratiedaethotaleatiea 83.85% .Overall,thediversiyandabundanceoffshatestuarinesites(K)werelowerhilethoseinLingdongsection(S)wereighr.Comparedwithtradioalsiods,oealetabarodinhdthdaagofstigssiivitydiou beusedaasupnttodalsysouesuryTsdyouldeichtesructuraladfucioaloaoofsoucs inthlowerrachsofthengRivendproidesicfaiosuprtsforthsteatcaagentdesatoofishc in this water area.

    KeyWordsEnvironmental DNA(eDNA) technology;Fish diversity;Dayang River

    大洋河發(fā)源于遼寧省岫巖縣偏嶺鄉(xiāng)一棵樹嶺南側(cè),流經(jīng)岫巖縣,至東港市黃土坎入黃海[1]。干流全長 230km ,流域面積 6504km2 ,是遼河與鴨綠江之間最大的河流,多年平均徑流量為31億 m3[2] 。大洋河水生態(tài)系統(tǒng)保護(hù)較好,大部分河段為沙礫底質(zhì)[。該流域內(nèi)無污染,且干流沒有大型水庫[3],魚類資源較為豐富[4]。目前關(guān)于大洋河的研究多集中在單一物種生物學(xué)研究[1,5-6]、水質(zhì)資源評估與保護(hù)[7-9]、洪水分析與預(yù)報(bào)技術(shù)方面[10-I]。迄今為止關(guān)于魚類資源多樣性的研究較少,僅在1980年對大洋河進(jìn)行了魚類資源的調(diào)查。由于近年來大洋河流域內(nèi)環(huán)境的變化和支流水壩的建設(shè)[3],及時(shí)準(zhǔn)確掌握水域內(nèi)的魚類生物量和時(shí)空動(dòng)態(tài)分布有助于開展魚類資源的保護(hù)、開發(fā)和合理利用。

    環(huán)境DNA宏條形碼技術(shù)(eDNAmetabarcoding)是一種利用DNA 條形碼和高通量測序(high throughput sequencing,HTS)相結(jié)合,對環(huán)境樣本中的eDNA進(jìn)行快速檢測,從而獲取環(huán)境樣本中DNA所屬物種的分類學(xué)信息和基因功能信息的研究方法[12-13]。其原理是收集物種通過細(xì)胞脫落、尿液/糞便排泄等方式殘留在環(huán)境中的DNA,然后利用特定的引物擴(kuò)增出物種(或生物群落)特定DNA序列,并通過高通量測序確定物種種類組成和多樣性特征[14-15]。與傳統(tǒng)檢測方法相比,eDNA具有較高的靈敏度,可以檢測到少量DNA拷貝,采樣方案相對簡單,即采樣 0.5~3.0L 水即可對水域中的生物進(jìn)行檢測[16-17]。最初,eDNA 技術(shù)被用于分析土壤環(huán)境中微生物的多樣性以及物種鑒定[18-19],后來被應(yīng)用于檢測評估淡水環(huán)境中生物的多樣性和種類[20-22],現(xiàn)在已被廣泛應(yīng)用于海洋和河口等不同生態(tài)系統(tǒng)的生物多樣性檢測、評估[23-24]該研究利用環(huán)境DNA宏條形碼技術(shù)對大洋河下游開展魚類資源調(diào)查,分析了大洋河下游魚類組成、生物量和多樣性指數(shù),以客觀評估大洋河下游魚類資源現(xiàn)狀,以期為大洋河的生態(tài)修復(fù)效果評估與方案優(yōu)化提供數(shù)據(jù)支撐。

    1材料與方法

    1.1采樣時(shí)間及采樣地點(diǎn)試驗(yàn)所用環(huán)境樣本為2023年10月20—23日采自大洋河下游的水域樣本,此次調(diào)查按離海遠(yuǎn)近在河流下游選取3個(gè)地點(diǎn),即河口(K)、黃土坎段(H)和嶺東段(S)。為了檢測到整個(gè)區(qū)域,分別在每個(gè)地點(diǎn)河流的斷截面取3個(gè)點(diǎn),包括河流兩岸和中間段,共9個(gè)位點(diǎn),具體采樣位點(diǎn)信息見圖1。

    圖1大洋河下游水樣采樣位點(diǎn)示意Fig.1Sampling sites in the downstream of the Dayang River

    1.2環(huán)境DNA樣本采集、處理及提取每個(gè)采樣位點(diǎn)采集表層、中層和底層水樣各 3L ,混勻。為了避免樣品交叉感染,取樣品和手套均為一次性,采樣后丟棄。各位點(diǎn)水樣在采集后 24h 內(nèi)處理完畢,使用 0.22μm 無菌濾膜真空抽濾,抽濾前用84消毒液對抽濾裝置進(jìn)行浸泡消毒處理,得到濾膜(各點(diǎn)位5\~7個(gè)平行),放置在無菌無酶的 2mL 離心管中保存至 -80°C 冰箱中。在DNA提取前,先在無菌管中將每個(gè)濾膜剪成小塊,再利用DNeasyPowerWaterKit試劑盒(QIAGEN)對濾膜進(jìn)行基因組DNA提取。最后,利用 1% 瓊脂糖凝膠電泳對提取的基因組DNA進(jìn)行質(zhì)量檢測。

    1.3PCR擴(kuò)增及高通量測序使用已報(bào)道的魚類eDNA通用引物Tele02-F(AAACTCGTGCCAGCCACC)和Tele02-R(GGGTATCTAATCCCAGTTTG)進(jìn)行 PCR 擴(kuò)增[25]。PCR 擴(kuò)增體系為 5×Fast Pfu Buffer 4μL,2.5 mmol/L dNTPs 2μL 正反引物( 5μmol/L )各 0.8μL 、FastPfu DNA Polymerase0.4μL 、模板DNA 10ng ,補(bǔ)雙蒸水至 20μL 。PCR反應(yīng)程序?yàn)? 個(gè)循環(huán);72°C 延伸 10min 。參照電泳初步定量結(jié)果,使用QuantiFlu-orT-ST藍(lán)色熒光定量系統(tǒng)(Promega 公司)進(jìn)行檢測定量,檢測完成后利用IlluminaPE250進(jìn)行文庫構(gòu)建。

    1.4數(shù)據(jù)統(tǒng)計(jì)與分析使用USEARCH軟件(version10,ht-tp://www.drive5.com/usearch/)和 GOLD 數(shù)據(jù)庫(https://gold.jgi.doe.gov/),采用denovo和reference相結(jié)合的方式去除嵌合體,按照 97% 相似性對OTU(operational taxonomic u-nits)進(jìn)行聚類分析,得到OTU的代表序列,生成OTU表格。采用UCLUST 算法(https://www.drive5.com/usearch/manual/uclust_algo.html)對序列相似性 97% 的OTU代表序列進(jìn)行分類學(xué)分析,統(tǒng)計(jì)各樣本在各個(gè)分類水平的群落組成。使用Mothur 軟件[26-27](version v.1.30.1,http://www.mothur.org/wi-ki/Schloss_SOP#Alpha_diversity)進(jìn)行 ∝ 多樣性分析,再使用在線分析軟件(http://www.cloud.biomicroclass.com/Cloud-Platform/SoftPage/FLO)繪制OTU花瓣圖[28]?;诟鞑蓸狱c(diǎn)物種序列豐度,利用R語言vegan包繪制物種相對豐度圖[29]。

    2 結(jié)果與分析

    2.1環(huán)境DNA測序結(jié)果對來自大洋河下游的9個(gè)水樣進(jìn)行了分析,共獲得9740369個(gè)原始序列,其中9732280個(gè)高質(zhì)量序列,平均長度為 171bp 。在序列相似性 ?97% 的情況下,通過聚類共獲得963個(gè)OTU。eDNA樣品的高通量測序統(tǒng)計(jì)數(shù)據(jù)如表1所示。

    表1樣本序列統(tǒng)計(jì)Table1 Statistics of the sample sequences

    2.2魚類物種組成將 OTU與GenBank數(shù)據(jù)庫(https://www.ncbi.nlm.nih.gov/genbank/)進(jìn)行比較,對注釋序列的OTU進(jìn)行分類。共檢測到61種魚類,隸屬于9目21科52屬,其中以鯉形目(Cypriniformes)最多,有34種(占檢測物種數(shù)的 56% );其次為蝦虎魚目(Gobiiformes),有10種(占檢測物種數(shù)的 16% )。不同采樣位點(diǎn)魚類屬水平相對豐度如圖2所示。整體來看,各采樣位點(diǎn)相對豐度排名前3位的屬為鰱鰱屬 Hypophthalmichthys 鯽屬Carassius鯪屬Planiliza 棘蝦虎魚屬Acanthogobius屬Coilia 副泥鰍屬Paramisgurnus未分類Unclassified 屬 Hemibaarbus原鯉屬Procypris 小鉤屬M(fèi)icrophysogobio鲇屬Silurus 其他Other縞虎魚屬Tridentiger100Hrrraageeeiaaree 80604020】H-1 H2H-3K-1K-2K-3S-1S-2 S-3樣本編號(hào) Sample number屬(Hypophthalmichthys)鮫屬(Planiliza)和屬(Coilia),鰱屬和鮫屬在K(河口)采樣位點(diǎn)相對豐度較高,屬在H(黃土坎段)采樣位點(diǎn)相對豐度較高?;诃h(huán)境DNA檢測結(jié)果的大洋河下游優(yōu)勢魚類物種統(tǒng)計(jì)見表2。

    2.3魚類多樣性分析此次測序得到的魚類樣本稀釋曲線如圖3所示,9個(gè)采樣位點(diǎn)測序深度均達(dá)到平臺(tái)期,證實(shí)了后續(xù)多樣性分析的可行性。由表3可知,在大洋河下游魚類物種的多樣性指標(biāo)中,Chao指數(shù)的變化范圍為 2 059.667~ 3116.824,物種豐富度指數(shù)的變化范圍為 1 397~2 451 。Shannon多樣性指數(shù)的變化范圍為 3.198 1~5.302 3, Simpson多樣性指數(shù)的變化范圍為 0.1011~0.2710 。ACE指數(shù)的變化范圍為2126.576\~3234.492,Pielou均勻度指數(shù)的變化范圍為 0.301394~0.470934. 。盡管各采樣位點(diǎn)的魚類物種多樣性指標(biāo)略有差異,但K-2的Simpson多樣性指數(shù)最高,其余指標(biāo)均為S-1最高。

    表2基于環(huán)境DNA檢測結(jié)果的大洋河下游優(yōu)勢魚類物種統(tǒng)計(jì)Table 2Statistics of the dominant fish species in the downstream of Dayang River basedon eDNA detection results
    圖3大洋河下游魚類樣本稀釋曲線

    2.4魚類物種差異分析利用在線分析軟件,根據(jù)獲得的樣本OTU統(tǒng)計(jì)表,繪制9個(gè)采樣位點(diǎn)樣本的OTU花瓣圖。各樣本共有和特有OTU如圖4a所示。由圖4a可知,9個(gè)采樣位點(diǎn)樣本共有OTU占比較高,且S-1和S-2這2個(gè)采樣位點(diǎn)具有較多的特有OTU。為了比較大洋河下游不同區(qū)域魚類物種的組成,對S、H、K3個(gè)地點(diǎn)進(jìn)行了分析,各地共有和特有OTU如圖4b所示。嶺東段(S)的特有OTU數(shù)量最多,黃土坎段(H)次之,河口(K)最少。

    3討論

    3.1大洋河下游魚類物種組成及多樣性大洋河處于遼河多雨區(qū),魚類種類繁多[4],河口連接黃渤海,一些廣鹽性魚類容易進(jìn)入大洋河繁殖生存[30]。近年來,由于環(huán)境破壞和水利工程建設(shè)等原因,漁業(yè)資源呈逐年減少趨勢[31],因此了解大洋河流域魚類群落的組成和結(jié)構(gòu)有利于稀有瀕危魚類的研究、保護(hù)和可持續(xù)開發(fā)利用。目前針對大洋河魚類研究方法主要為傳統(tǒng)捕撈法,該方法對魚類種群群落具有破壞性,實(shí)施過程較為費(fèi)時(shí)、費(fèi)力,且需要專業(yè)研究人員進(jìn)行分類、鑒定[32]。該研究利用eDNA技術(shù)分析了大洋河下游3個(gè)區(qū)域共9個(gè)位點(diǎn)的魚類組成,共檢測到61種魚類,隸屬于9目21科52屬。其中,鰱屬的相對豐度較高,其次為鮫屬和屬。

    表3大洋河下游魚類群落的多樣性指標(biāo)Table3The diversity indices of the fish community in the downstream of Dayang River
    Fig.3Rarefaction curve of the fish samples in the downstream of Dayang River圖4大洋河下游各采樣位點(diǎn)魚類物種OTU花瓣圖

    Fig.4Petal map of OTU of fish species at different sampling sites in the downstream of Dayang River

    該研究采用多樣性指標(biāo)來表征大洋河下游水域的魚類物種多樣性。多樣性指標(biāo)的大小可以反映物種群落的豐度和多樣性[33],包括物種豐富度指數(shù)、Shannon 多樣性指數(shù)和Pielou均勻度指數(shù)等。其中,Chao指數(shù)和物種豐富度指數(shù)表征物種豐度;Shannon多樣性指數(shù)、Simpson多樣性指數(shù)、ACE指數(shù)和Pielou均勻度指數(shù)表征物種的豐富度。9個(gè)采樣位點(diǎn)中S-1位點(diǎn)具有較高的Chao指數(shù)、物種豐富度指數(shù)、Shan-non多樣性指數(shù)和Pielou均勻度指數(shù),表明相較于其他位點(diǎn),S-1具有較高的魚類物種豐度和多樣性。除Simpson多樣性指數(shù)外,各地點(diǎn)3個(gè)位點(diǎn)所有指標(biāo)的平均值均表現(xiàn)為 K。

    10月大洋河水溫呈上升趨勢,而溫度的升高有利于多種喜溫浮游動(dòng)植物的生長和繁殖,生物多樣性增加,初級(jí)生產(chǎn)力提高,有利于魚類的繁衍聚集,因此魚類的多樣性更高[34]。

    3.2傳統(tǒng)捕撈法與eDNA技術(shù)檢測結(jié)果的比較傳統(tǒng)捕撈法與eDNA技術(shù)在魚類資源調(diào)查中存在一定差異。凌嵐馨等[32]通過環(huán)境DNA技術(shù)與傳統(tǒng)捕撈法對崇明島內(nèi)河魚類多樣性進(jìn)行了調(diào)查研究,結(jié)果發(fā)現(xiàn)eDNA技術(shù)檢測到的魚類物種數(shù)量多于傳統(tǒng)捕獲法所獲得的物種數(shù)量,且檢測的物種多樣性明顯高于傳統(tǒng)捕獲法。根據(jù)歷史文獻(xiàn)[4記載,采用傳統(tǒng)捕撈法在大洋河下游共監(jiān)測到魚類47種,隸屬于22科。與該研究結(jié)果相比,相同種類有16種,包括花棒魚(Abbottinarivularis)鯽魚(Carassiuscuvieri)長吻魚骨(Hemibarbuslon-girostris)白鰱(Hypophthalmichthysmolitrix)等;新增了烏原鯉(Procyprismera)異育銀鯽(Carassiusgibelio)矛形蝦虎魚(Acanthogobius hasta)、副鰍(Paramisgurnus dabryanus)等 45種魚類;在歷史文獻(xiàn)中有記載而此次沒有調(diào)查到的魚類種類有31種,包括彈涂魚(Periophthalmuscantonensis)、香魚(Plecoglossusaltivelis)黃顙魚(Pelteobagrusfulvidraco姑魚(Johnius belangeri)、烏(Ophiocephalusargus)、日本鰻(Anguillajaponica)和鯰魚(Parasilurusasotus)等(表4)。利用eDNA技術(shù)檢測到的魚類種數(shù)遠(yuǎn)大于傳統(tǒng)捕撈法獲得的魚類種數(shù),造成這種差異的原因可能包含以下方面: ① eDNA技術(shù)具有較高的靈敏性,可以檢測到生物量較低的物種[33]; ② 采樣時(shí)間間隔較長,大洋河下游的優(yōu)勢魚類物種和群落組成發(fā)生了變化[35]; ③ 魚類習(xí)性和環(huán)境因素的變化可能影響了魚類多樣性[36]; ④ 由于eDNA測定選擇的引物特異性不強(qiáng),造成數(shù)據(jù)庫比對結(jié)果存在一定差異,從而導(dǎo)致結(jié)果不準(zhǔn)確[32]; ⑤ 比對數(shù)據(jù)庫不夠完整,使測得的序列無法進(jìn)行比對。該研究發(fā)現(xiàn)有7.46% 的物種未能進(jìn)行分類,這種現(xiàn)象在eDNA技術(shù)被廣泛應(yīng)用于魚類資源調(diào)查前就有出現(xiàn)[37]。此外,現(xiàn)有數(shù)據(jù)庫中還存在大量的錯(cuò)誤信息,這也影響了研究結(jié)果的準(zhǔn)確性[7]。

    表4大洋河下游魚類資源比較Table 4 Comparison of fish resources in the downstream of Dayang River

    3.3魚類生物量分析凌建忠等[38]研究表明生態(tài)系統(tǒng)中生物的eDNA豐度與生物量密切相關(guān),可用于生物量的評估。Tillotson等[39利用eDNA技術(shù)評估不同河流中鮭魚的生物量,證實(shí)了eDNA技術(shù)用于物種生物量評估的可行性。該研究eDNA檢測結(jié)果表明,大洋河下游10月魚類的優(yōu)勢種為鰱屬、鮫屬和屬,其中魴屬主要為刀(C.nasus),相對豐度為 16.08% 。作為一種中小型洄游魚類,刀原產(chǎn)于我國長江地區(qū),與河豚、魚并稱“長江三鮮”。據(jù)歷史資料記載我國野生刀資源豐富[40,但近年來由于水利建設(shè)、環(huán)境污染和過度捕撈等原因,我國野生刀資源銳減,2018年被世界自然保護(hù)聯(lián)盟(IUCN)列為瀕危魚類[4I]。遼寧是北方地區(qū)刀的主要洄游區(qū)域,主要包括遼河、大洋河和鴨綠江[42]。目前,僅在大洋河發(fā)現(xiàn)小規(guī)模刀漁汛[5.43]。刀作為我國重要的經(jīng)濟(jì)魚類之一,近年來對其開展了大規(guī)模的放流活動(dòng),用來恢復(fù)刀魚類資源。該研究在大洋河下游3個(gè)地點(diǎn)均檢測到刀魴,并發(fā)現(xiàn)其為調(diào)查區(qū)域的優(yōu)勢種,表明大洋河水域的刀資源恢復(fù)良好。

    4結(jié)論

    該研究利用eDNA宏條形碼技術(shù)分析了大洋河下游魚類物種的多樣性,獲得了較為豐富的魚類物種信息,共鑒定出61種魚類,隸屬9目21科52屬。與歷史調(diào)查數(shù)據(jù)相比,該研究發(fā)現(xiàn)通過eDNA技術(shù)獲得的魚類多樣性高于傳統(tǒng)捕撈法,證實(shí)了eDNA技術(shù)應(yīng)用于魚類資源調(diào)查的可行性。下一步將繼續(xù)擴(kuò)大采樣范圍和采樣時(shí)間,從時(shí)間和空間上開展更加全面的大洋河魚類多樣性調(diào)查,為大洋河魚類資源的保護(hù)、開發(fā)和合理利用提供科學(xué)依據(jù)。

    參考文獻(xiàn)

    [1]鄭建平,王芳.大洋河河流生態(tài)需水研究[J].河海大學(xué)學(xué)報(bào)(自然科學(xué)版),2006,34(5):502-504.

    [2]張賀.基于HYDRUS-MODFLOW的大洋河岫巖流域地表水與地下水耦 合模型[J].陜西水利,2020(3):7-9,14.

    [3]葉天舒,金鑫,蔡濤.遼寧省大洋河流域暴雨洪水特征分析[J].中國防 汛抗旱,2019,29(7):29-32.

    [4]史為良.大洋河及其毗鄰河流的魚類區(qū)系特征[J].水產(chǎn)科學(xué),1985,4 (4) :53-57,52.

    [5]張健,楊培民,姜濤,等.基于耳石微化學(xué)的大洋河刀生境履歷研究 [J].水生生物學(xué)報(bào),2024,48(1):130-137.

    [6]王慶志,常亞青.大洋河河蜆的繁殖生物學(xué)研究[J].大連水產(chǎn)學(xué)院學(xué) 報(bào),2010,25(1) :8-13.

    [7]衛(wèi)寶泉,武藝.大洋河流域河口濕地生態(tài)環(huán)境評價(jià)三維系統(tǒng)的設(shè)計(jì)與實(shí) 現(xiàn)[J].地理科學(xué)研究,2017,6(2):65-71.

    [8]胡曉松.不同分布式水文模型在大洋河流域?qū)Ρ冗\(yùn)用研究[J].黑龍江 水利科技,2015,43(7):11-13.

    [9]趙靜,趙學(xué)義.不同土地利用變化情景模式對大洋河水土流失影響研究 [J].水利天地,2015(5):20-22,43.

    [10]張靜.大洋河洪水預(yù)報(bào)技術(shù)研究[J].陜西水利,2020(11):65-67.

    [11]張蓉蓉.大洋河\"8·03\"暴雨洪水分析[C]//遼寧省水利學(xué)會(huì).水與水 技術(shù)(第8輯).北京:中國水利水電出版社,2018:19-21.

    [12]FICETOLA G F,MIAUD C,POMPANON F,et al.Species detection using environmental DNA from water samples[J].Biol Lett,2008,4(4):423- 425.

    [13]LODGE D M,TURNER CR,JERDE C L,et al.Conservation in a cup of water:Estimating biodiversity andpopulation abundance from environmental DNA[J].Mol Ecol,2012,21(11):2555-2558.

    [14] THOMSEN P F,WILLERSLEV E.Environmental DNA-An emerging tool in conservation for monitoring past and present biodiversity[J].Biol Conserv,2015,183:4-18.

    [15] LACOURSIERE-ROUSSEL A,COTE G,LECLERC V,et al.Quantifying relative fish abundance with eDNA:A promising tool for fisheries management[J].J Appl Ecol,2016,53(4) :1148-1157.

    [16]JERDECL,MAHONAR,CHADDERTONWL,etal.\"Sight-unseen\" detection of rare aquatic species using environmental DNA[J].Conserv Lett,2011,4(2):150-157.

    [17] MAHON AR,JERDE C L,GALASKA M,et al.Validation of eDNA surveillance sensitivity for detection of Asian carps incontrolledand field experiments[J].PLoS One,2013,8(3) :1-6.

    [18]OGRAMA,SAYLERGS,BARKAY T.Theextractionand purification of microbial DNA from sediments[J].J Microbiol Meth,1987,7(2/3):57- 66.

    [19]WOLCOTT R D,GONTCHAROVA V,SUN Y,et al.Evaluation of the bacterial diversity among and within individual venous leg ulcers using bacterial tag-encoded FLXand Titanium amplicon pyrosequencing and metagenomic approaches[J].BMC Microbiol,2009,9:1-11.

    [20] PILLIOD D S,GOLDBERG C S,ARKLE R S,et al.Estimating occupancy and abundance of stream amphibians using environmental DNA from filtered water samples[J].Can J Fish Aquat Sci,2013,70(8) :1123-1130.

    [21] TAKAHARA T,MINAMOTO T,YAMANAKA H,et al.Estimation of fish biomass using environmental DNA[J].PLoS One,2O12,7(4) :1-8.

    [22]THOMSENPF,KIELGASTJ,IVERSENL L,etal.Monitoring endangered freshwater biodiversity using environmental DNA[J].Mol Ecol,

    [23]YAMAMOTO S,MASUDAR,SATO Y,et al.Environmental DNAme-tabarcoding reveals local fish communities in a species-rich coastal sea[J].Sci Rep,2017,7:1-12.

    [24]THOMSENPF,KIELGASTJ,IVERSENLL,etal.Detectionofadiversemarine fish fauna using environmental DNA from seawater samples[J].PLoS One,2012,7(8) :1-9.

    [25] ZHANG S,LU Q,WANG Y Y,et al.Assessment of fish communities usingenvironmental DNA:Effect of spatial sampling design inlentic systems ofdifferent sizes[J].Mol Ecol Res,2020,20(1) :242-255.

    [26] SCHLOSS P D,GEVERS D,WESTCOTT S L.Reducing the effects of PCRamplification and sequencing artifacts on 16S rRNA-based studies[J].PLoS One,2011,6(12) :1-14.

    [27] WHITTAKER R H.Evolution and measurement of species diversity[J].Taxon,1972,21(2/3) :213-251.

    [28]FOUTS D E,SZPAKOWSKI S,PURUSHE J,et al.Next generation se-quencing to define prokaryotic and fungal diversity in the bovine rumen[J].PLoS One,2012,7(11) :1-11.

    [29]JAMI E,ISRAEL A,KOTSER A,et al.Exploring the bovine rumen bacte-rial community from birth to adulthood[J].ISME J,2013,7(6):1069-1079.

    [30]于硯民,董志剛,曹仁江.鴨綠江口濱海濕地生物多樣性與保護(hù)對策[J].資源科學(xué),1999,21(1):57-62.

    [31]楊永利.大洋河沙里寨水文站水沙變化與驅(qū)動(dòng)因素分析[J].水利技術(shù)監(jiān)督,2020(5):105-108,205.

    [32]凌嵐馨,范共,胡云,等.環(huán)境DNA技術(shù)與傳統(tǒng)捕撈揭示崇明島內(nèi)河魚類多樣性[J].上海海洋大學(xué)學(xué)報(bào),2022,31(6):1434-1444.

    [33]楊青,李宏俊,李洪波,等.海洋生物多樣性評價(jià)方法綜述[J].海洋環(huán)境科學(xué),2013,32(1):157-160.

    [34]魏洪祥,楊培民,蔣湘輝,等.大洋河刀洄游河段浮游生物群落特征研究[J].水產(chǎn)養(yǎng)殖,2021,42(1):29-36.

    [35]PLOUGH L V,OGBURN M B,F(xiàn)ITZGERALD C L,et al.EnvironmentalDNA analysis of river herring in Chesapeake Bay: A powerful tool for mo-nitoring threatened keystone species[J].PLoS One,2018,13(11) :1-25.

    [36]DE SOUZA L S,GODWIN J C,RENSHAW M A,et al.EnvironmentalDNA(eDNA) detection probability is influenced by seasonal activity oforganisms[J].PLoS One,2016,11(10):1-15.

    [37]MIYA M.Environmental DNA metabarcoding:A novel method for biodi-versity monitoring of marine fish communities[J].Ann Rev Mar Sci,2022,14:161-185.

    [38]凌建忠,姜亞洲,孫鵬,等.環(huán)境DNA技術(shù)在象山港水域魚類多樣性調(diào)查中的應(yīng)用與評估[J].中國水產(chǎn)科學(xué),2021,28(2):205-214.

    [39]TILLOTSONMD,KELLYRP,DUDAJJ,etal.Concentrations ofenvi-ronmental DNA (eDNA)reflect spawning salmon abundance at fine spa-tial and temporal scales[J].Biol Conserv,2O18,22O:1-11.

    [40]解玉浩.東北地區(qū)淡水魚類[M].沈陽:遼寧科學(xué)技術(shù)出版社,2007.

    [41]洪珍珍,梅肖樂,王苗苗,等.長江流域刀生物學(xué)特性及資源保護(hù)現(xiàn)狀[J].水產(chǎn)養(yǎng)殖,2023,44(7):36-41.

    [42]張賽賽.中國遼寧地區(qū)刀(Coilianasus)遺傳多樣性研究[D].大連:大連海洋大學(xué),2016.

    [43]劉中,于偉君,石俊艷,等.雙臺(tái)子河漁業(yè)環(huán)境調(diào)查報(bào)告[J].水產(chǎn)科學(xué),1992,11( 10) 5-9.

    猜你喜歡
    魚類位點(diǎn)物種
    基于OBE理念的魚類生理學(xué)課程教學(xué)改革實(shí)踐
    小麥條銹菌SNP分子標(biāo)記的開發(fā)與應(yīng)用
    棉花種子油分含量QTL的定位及遺傳效應(yīng)分析
    上海海洋大學(xué)解碼魚類胚胎中的新生轉(zhuǎn)錄本
    TLR4基因多態(tài)性與T2DM病人并發(fā)缺血性心腦血管病的關(guān)系探討
    基于eDNA宏條形碼技術(shù)的巢湖夏季魚類群落結(jié)構(gòu)及多樣性
    流速-水溫對草魚生殖洄游的交互作用
    維生素D受體基因多態(tài)性、維生素D水平與急性非ST段抬高型心肌梗死患者預(yù)后的相關(guān)性分析
    大醫(yī)生(2025年12期)2025-08-11 00:00:00
    吃光入侵物種真的是解決之道嗎?
    英語世界(2023年10期)2023-11-17 09:18:18
    回首2018,這些新物種值得關(guān)注
    蜜桃在线观看..| 亚洲中文av在线| 亚洲精品一区蜜桃| 欧美变态另类bdsm刘玥| 人体艺术视频欧美日本| 美女脱内裤让男人舔精品视频| 精品国产一区二区三区久久久樱花| 1024视频免费在线观看| 日韩人妻精品一区2区三区| 一区二区三区四区激情视频| 黄网站色视频无遮挡免费观看| 国产成人精品在线电影| 婷婷丁香在线五月| 丝袜喷水一区| 国产欧美日韩综合在线一区二区| 啦啦啦在线免费观看视频4| 精品亚洲成a人片在线观看| 1024香蕉在线观看| 这个男人来自地球电影免费观看| 成人影院久久| videos熟女内射| 欧美日韩视频精品一区| 超碰成人久久| 国产免费又黄又爽又色| 免费观看av网站的网址| www.自偷自拍.com| 国产亚洲精品第一综合不卡| 久久久久久免费高清国产稀缺| 美国免费a级毛片| 看免费成人av毛片| 久久影院123| 欧美+亚洲+日韩+国产| 久久久精品94久久精品| 制服诱惑二区| 亚洲九九香蕉| 亚洲精品久久成人aⅴ小说| 啦啦啦啦在线视频资源| 黄色怎么调成土黄色| 菩萨蛮人人尽说江南好唐韦庄| 伊人亚洲综合成人网| 日韩熟女老妇一区二区性免费视频| 91成人精品电影| 久久久久国产精品人妻一区二区| 99九九在线精品视频| 在线av久久热| 性色av乱码一区二区三区2| 国产又色又爽无遮挡免| 三上悠亚av全集在线观看| 国产精品久久久久久精品电影小说| 校园人妻丝袜中文字幕| 亚洲av男天堂| 欧美在线一区亚洲| 啦啦啦 在线观看视频| 考比视频在线观看| 久久精品国产a三级三级三级| 亚洲欧美日韩高清在线视频 | 国产亚洲一区二区精品| 中国国产av一级| 免费观看a级毛片全部| 久久精品国产亚洲av涩爱| 人人妻,人人澡人人爽秒播 | 午夜91福利影院| 美女扒开内裤让男人捅视频| 亚洲美女黄色视频免费看| 午夜免费观看性视频| 美国免费a级毛片| 亚洲综合色网址| 欧美日韩视频精品一区| 大型av网站在线播放| 你懂的网址亚洲精品在线观看| 97精品久久久久久久久久精品| 久久国产精品影院| 在线亚洲精品国产二区图片欧美| tube8黄色片| 国产亚洲一区二区精品| 在线 av 中文字幕| 9191精品国产免费久久| 伦理电影免费视频| 欧美人与性动交α欧美精品济南到| 免费观看人在逋| 亚洲成av片中文字幕在线观看| 欧美精品高潮呻吟av久久| 亚洲国产欧美日韩在线播放| 老司机影院成人| 亚洲国产欧美网| 黑人猛操日本美女一级片| 最黄视频免费看| 精品久久久久久电影网| 又大又黄又爽视频免费| 亚洲天堂av无毛| 欧美人与性动交α欧美精品济南到| 青草久久国产| 成年人黄色毛片网站| 亚洲熟女精品中文字幕| 最黄视频免费看| 免费人妻精品一区二区三区视频| 亚洲av电影在线进入| 午夜影院在线不卡| 国产欧美日韩一区二区三区在线| 777久久人妻少妇嫩草av网站| 亚洲 欧美一区二区三区| 久久综合国产亚洲精品| 中文精品一卡2卡3卡4更新| 国产熟女午夜一区二区三区| 国产在线一区二区三区精| 免费av中文字幕在线| 国产男人的电影天堂91| 午夜av观看不卡| 久久99一区二区三区| 国产一区二区三区综合在线观看| h视频一区二区三区| 国产无遮挡羞羞视频在线观看| 不卡av一区二区三区| 色视频在线一区二区三区| √禁漫天堂资源中文www| 国产精品人妻久久久影院| 99香蕉大伊视频| 永久免费av网站大全| 中文字幕人妻丝袜制服| 国产精品国产av在线观看| av在线app专区| 老汉色∧v一级毛片| 亚洲精品一卡2卡三卡4卡5卡 | 女人被躁到高潮嗷嗷叫费观| 亚洲国产精品成人久久小说| 免费一级毛片在线播放高清视频 | 午夜福利视频精品| 多毛熟女@视频| 久久久久久亚洲精品国产蜜桃av| 久久精品国产亚洲av高清一级| 午夜福利免费观看在线| 亚洲,欧美精品.| 久久久精品免费免费高清| 亚洲精品久久成人aⅴ小说| 亚洲 国产 在线| 纯流量卡能插随身wifi吗| 男人舔女人的私密视频| 三上悠亚av全集在线观看| 国产亚洲av片在线观看秒播厂| videos熟女内射| 国产又爽黄色视频| 午夜福利视频在线观看免费| 国产日韩欧美亚洲二区| 好男人电影高清在线观看| 晚上一个人看的免费电影| 亚洲国产欧美日韩在线播放| 婷婷成人精品国产| 飞空精品影院首页| 国产精品久久久久久精品古装| 欧美精品亚洲一区二区| 我要看黄色一级片免费的| 老司机深夜福利视频在线观看 | 好男人电影高清在线观看| 亚洲专区国产一区二区| 色视频在线一区二区三区| 制服诱惑二区| 国产精品一国产av| 色网站视频免费| 国产又色又爽无遮挡免| 赤兔流量卡办理| 在线观看免费视频网站a站| kizo精华| 国产欧美日韩精品亚洲av| 亚洲av美国av| 精品人妻一区二区三区麻豆| 天天躁夜夜躁狠狠躁躁| 狠狠精品人妻久久久久久综合| 老司机靠b影院| 日韩制服丝袜自拍偷拍| 日韩中文字幕欧美一区二区 | 涩涩av久久男人的天堂| 老熟女久久久| 天堂俺去俺来也www色官网| 亚洲欧美中文字幕日韩二区| 欧美日韩av久久| 亚洲五月色婷婷综合| 黄色 视频免费看| 欧美精品高潮呻吟av久久| 亚洲国产精品999| 90打野战视频偷拍视频| 男女免费视频国产| 国产色视频综合| 精品国产一区二区久久| 多毛熟女@视频| 精品少妇内射三级| 日韩免费高清中文字幕av| 亚洲精品成人av观看孕妇| 男人舔女人的私密视频| 精品亚洲乱码少妇综合久久| 欧美老熟妇乱子伦牲交| 性色av乱码一区二区三区2| 波多野结衣一区麻豆| 九草在线视频观看| 亚洲av电影在线进入| 国产精品 国内视频| 18在线观看网站| 桃花免费在线播放| 欧美久久黑人一区二区| 久久久久久久久免费视频了| 五月天丁香电影| 午夜福利视频在线观看免费| xxxhd国产人妻xxx| 人人澡人人妻人| 人人妻,人人澡人人爽秒播 | 成人国语在线视频| svipshipincom国产片| 亚洲综合色网址| 不卡av一区二区三区| 下体分泌物呈黄色| 十八禁网站网址无遮挡| 99久久精品国产亚洲精品| 黄色 视频免费看| 另类亚洲欧美激情| 亚洲美女黄色视频免费看| 亚洲自偷自拍图片 自拍| 日日摸夜夜添夜夜爱| 亚洲精品日本国产第一区| 五月开心婷婷网| 欧美黑人欧美精品刺激| 人成视频在线观看免费观看| 欧美日韩视频高清一区二区三区二| 精品国产国语对白av| 欧美97在线视频| 一区在线观看完整版| 欧美精品亚洲一区二区| 免费少妇av软件| 婷婷色综合大香蕉| 建设人人有责人人尽责人人享有的| 亚洲精品美女久久av网站| 中文精品一卡2卡3卡4更新| 国产精品久久久久久精品古装| 久久精品成人免费网站| 每晚都被弄得嗷嗷叫到高潮| 五月开心婷婷网| 免费高清在线观看日韩| 国产精品 欧美亚洲| 国产男女内射视频| 一区二区三区乱码不卡18| 久久国产精品男人的天堂亚洲| 国产精品一二三区在线看| 丝袜脚勾引网站| 女性被躁到高潮视频| 亚洲色图 男人天堂 中文字幕| 欧美另类一区| 一本综合久久免费| 夜夜骑夜夜射夜夜干| 国产高清视频在线播放一区 | 人人妻人人澡人人看| 久久久国产一区二区| 午夜福利视频精品| 午夜两性在线视频| 精品人妻1区二区| 午夜福利一区二区在线看| 人体艺术视频欧美日本| 青草久久国产| av一本久久久久| 91九色精品人成在线观看| 又紧又爽又黄一区二区| 日韩视频在线欧美| 99国产综合亚洲精品| 亚洲一卡2卡3卡4卡5卡精品中文| 欧美少妇被猛烈插入视频| 又黄又粗又硬又大视频| av线在线观看网站| 91精品三级在线观看| 1024视频免费在线观看| 一本—道久久a久久精品蜜桃钙片| 国产1区2区3区精品| 日韩一区二区三区影片| 青春草视频在线免费观看| videosex国产| 99国产精品一区二区蜜桃av | 波多野结衣一区麻豆| 一区二区av电影网| 男女边吃奶边做爰视频| 我要看黄色一级片免费的| 丝袜在线中文字幕| avwww免费| 国产伦人伦偷精品视频| bbb黄色大片| 亚洲欧美日韩高清在线视频 | 黄频高清免费视频| 亚洲人成电影免费在线| svipshipincom国产片| 2021少妇久久久久久久久久久| 亚洲国产成人一精品久久久| 国产在线视频一区二区| 亚洲激情五月婷婷啪啪| 爱豆传媒免费全集在线观看| 99国产精品一区二区三区| 搡老岳熟女国产| 精品久久久久久久毛片微露脸 | 国产午夜精品一二区理论片| 国产男女超爽视频在线观看| 精品国产一区二区久久| 黑人欧美特级aaaaaa片| 一二三四在线观看免费中文在| 久久国产精品男人的天堂亚洲| 天天添夜夜摸| 久9热在线精品视频| 久久ye,这里只有精品| 亚洲三区欧美一区| 成年动漫av网址| 日本a在线网址| 亚洲成国产人片在线观看| 9191精品国产免费久久| 精品人妻一区二区三区麻豆| 免费女性裸体啪啪无遮挡网站| 国产精品人妻久久久影院| 老司机亚洲免费影院| 国产免费又黄又爽又色| 我要看黄色一级片免费的| 久久这里只有精品19| 黄色a级毛片大全视频| 一边摸一边做爽爽视频免费| 成人午夜精彩视频在线观看| 婷婷色麻豆天堂久久| 久久久久久久精品精品| 性色av一级| 久久久国产精品麻豆| 久久精品久久久久久噜噜老黄| 一级毛片女人18水好多 | 精品免费久久久久久久清纯 | 国产精品欧美亚洲77777| 2021少妇久久久久久久久久久| 日韩电影二区| 在线看a的网站| 爱豆传媒免费全集在线观看| 免费av中文字幕在线| 久久精品亚洲熟妇少妇任你| 国产真人三级小视频在线观看| 大话2 男鬼变身卡| 亚洲国产欧美在线一区| 亚洲欧美一区二区三区国产| 精品一品国产午夜福利视频| 纵有疾风起免费观看全集完整版| 国产亚洲欧美精品永久| 咕卡用的链子| 欧美精品高潮呻吟av久久| 国产熟女欧美一区二区| 亚洲av男天堂| 免费在线观看视频国产中文字幕亚洲 | 男女床上黄色一级片免费看| 国产一区二区三区av在线| 欧美少妇被猛烈插入视频| 少妇 在线观看| 精品少妇黑人巨大在线播放| 国产成人免费观看mmmm| 国产精品99久久99久久久不卡| a级片在线免费高清观看视频| 大香蕉久久成人网| av国产精品久久久久影院| 欧美亚洲日本最大视频资源| 国产欧美日韩精品亚洲av| 免费在线观看黄色视频的| 亚洲欧美一区二区三区久久| 国产野战对白在线观看| 人妻 亚洲 视频| 午夜福利视频在线观看免费| 亚洲黑人精品在线| 精品一区在线观看国产| 下体分泌物呈黄色| 国产又色又爽无遮挡免| 国产熟女午夜一区二区三区| 亚洲情色 制服丝袜| 欧美日韩亚洲高清精品| 成人亚洲精品一区在线观看| 亚洲av男天堂| 久久精品亚洲熟妇少妇任你| 国产av一区二区精品久久| 婷婷色麻豆天堂久久| 久久青草综合色| 国产野战对白在线观看| 美女大奶头黄色视频| 久久综合国产亚洲精品| 欧美成人精品欧美一级黄| 免费观看av网站的网址| 亚洲国产欧美在线一区| 久久国产精品影院| 高清视频免费观看一区二区| 欧美日韩成人在线一区二区| 国产一区二区 视频在线| 国产一级毛片在线| 女人久久www免费人成看片| 日本一区二区免费在线视频| 不卡av一区二区三区| 精品一区二区三卡| 亚洲激情五月婷婷啪啪| 91老司机精品| 欧美+亚洲+日韩+国产| 丝袜美腿诱惑在线| 国产精品免费大片| 久久久久久免费高清国产稀缺| 国产免费现黄频在线看| 国产成人欧美在线观看 | 女性被躁到高潮视频| 99国产精品免费福利视频| 波多野结衣一区麻豆| 一级a爱视频在线免费观看| 国产一区二区激情短视频 | 一本—道久久a久久精品蜜桃钙片| 亚洲精品国产区一区二| 午夜福利在线免费观看网站| 久久九九热精品免费| 在线 av 中文字幕| 日韩,欧美,国产一区二区三区| 久久这里只有精品19| 秋霞在线观看毛片| 国产熟女欧美一区二区| 高潮久久久久久久久久久不卡| 男人舔女人的私密视频| 女性生殖器流出的白浆| 大码成人一级视频| 亚洲五月色婷婷综合| 少妇裸体淫交视频免费看高清 | 日韩大码丰满熟妇| 欧美久久黑人一区二区| 真人做人爱边吃奶动态| 色婷婷久久久亚洲欧美| 搡老岳熟女国产| 这个男人来自地球电影免费观看| 黄色a级毛片大全视频| 少妇粗大呻吟视频| videosex国产| 日韩中文字幕视频在线看片| 一二三四在线观看免费中文在| 一本一本久久a久久精品综合妖精| 99香蕉大伊视频| 亚洲欧美一区二区三区久久| 精品一区二区三卡| 亚洲精品久久久久久婷婷小说| 啦啦啦中文免费视频观看日本| 女人高潮潮喷娇喘18禁视频| 交换朋友夫妻互换小说| 丰满迷人的少妇在线观看| 午夜精品国产一区二区电影| 老司机影院毛片| 又粗又硬又长又爽又黄的视频| 99热全是精品| 久久久国产精品麻豆| 日韩中文字幕视频在线看片| 久久国产精品影院| 岛国毛片在线播放| 国产成人免费观看mmmm| 丝袜脚勾引网站| 人人妻人人爽人人添夜夜欢视频| 精品国产乱码久久久久久小说| 国产精品久久久久成人av| 欧美黑人精品巨大| 一区在线观看完整版| 女人爽到高潮嗷嗷叫在线视频| 国产主播在线观看一区二区 | 欧美国产精品va在线观看不卡| 亚洲成人国产一区在线观看 | 啦啦啦啦在线视频资源| 黑丝袜美女国产一区| 18禁裸乳无遮挡动漫免费视频| 超碰成人久久| 国产日韩欧美在线精品| 老司机午夜十八禁免费视频| 中文乱码字字幕精品一区二区三区| 久久人妻福利社区极品人妻图片 | 国产精品久久久久久精品古装| 国产高清不卡午夜福利| 2018国产大陆天天弄谢| 热re99久久精品国产66热6| 亚洲人成77777在线视频| 高清av免费在线| 精品国产超薄肉色丝袜足j| 国产高清不卡午夜福利| 久久久久精品国产欧美久久久 | 亚洲欧洲日产国产| 99久久人妻综合| 91国产中文字幕| 国精品久久久久久国模美| 男人添女人高潮全过程视频| 久久精品久久久久久噜噜老黄| 亚洲专区国产一区二区| 久久国产精品大桥未久av| 久久 成人 亚洲| 国产一卡二卡三卡精品| 国产成人免费观看mmmm| 麻豆乱淫一区二区| 亚洲精品美女久久久久99蜜臀 | 在线观看www视频免费| 黑丝袜美女国产一区| 精品一区在线观看国产| 久久久久久人人人人人| 美女脱内裤让男人舔精品视频| 制服人妻中文乱码| 免费黄频网站在线观看国产| 免费高清在线观看视频在线观看| 咕卡用的链子| 欧美黑人欧美精品刺激| 精品国产乱码久久久久久小说| 丝袜脚勾引网站| 国产在线免费精品| 又粗又硬又长又爽又黄的视频| 老汉色∧v一级毛片| 久久久国产精品麻豆| √禁漫天堂资源中文www| 少妇粗大呻吟视频| 男的添女的下面高潮视频| av国产精品久久久久影院| 亚洲精品久久成人aⅴ小说| 久久这里只有精品19| 国产成人欧美| 国产黄色免费在线视频| 午夜福利影视在线免费观看| 亚洲精品美女久久av网站| 免费av中文字幕在线| 日本91视频免费播放| 亚洲美女黄色视频免费看| 日本欧美国产在线视频| 看十八女毛片水多多多| 天堂8中文在线网| 免费看不卡的av| 少妇人妻久久综合中文| 午夜两性在线视频| 精品福利永久在线观看| 欧美日韩国产mv在线观看视频| 黄色一级大片看看| 亚洲精品国产色婷婷电影| 日本欧美视频一区| 成人国产av品久久久| 乱人伦中国视频| 18在线观看网站| 欧美国产精品va在线观看不卡| 日韩熟女老妇一区二区性免费视频| 色综合欧美亚洲国产小说| 精品亚洲乱码少妇综合久久| 51午夜福利影视在线观看| 国产精品三级大全| videos熟女内射| 女警被强在线播放| 日韩中文字幕欧美一区二区 | 精品卡一卡二卡四卡免费| www.999成人在线观看| 久久人人爽av亚洲精品天堂| kizo精华| 日韩av免费高清视频| 国产亚洲精品第一综合不卡| 日本wwww免费看| 免费在线观看影片大全网站 | 精品欧美一区二区三区在线| 新久久久久国产一级毛片| 亚洲欧美精品综合一区二区三区| 母亲3免费完整高清在线观看| 99国产精品一区二区三区| 1024香蕉在线观看| 男人操女人黄网站| 久久久久国产精品人妻一区二区| 日韩人妻精品一区2区三区| 亚洲欧美一区二区三区国产| 国产精品av久久久久免费| 午夜免费鲁丝| 一本一本久久a久久精品综合妖精| 精品一区在线观看国产| h视频一区二区三区| 久久久久视频综合| 嫩草影视91久久| 中文字幕人妻熟女乱码| 亚洲欧美激情在线| 母亲3免费完整高清在线观看| 男人爽女人下面视频在线观看| 免费女性裸体啪啪无遮挡网站| netflix在线观看网站| 国产在线免费精品| av在线播放精品| 久久av网站| 国产男女内射视频| 国产精品香港三级国产av潘金莲 | 黄色视频在线播放观看不卡| 少妇裸体淫交视频免费看高清 | 在线av久久热| 男女午夜视频在线观看| 午夜老司机福利片| 久久久久久人人人人人| 亚洲色图 男人天堂 中文字幕| 亚洲欧美激情在线| 18在线观看网站| 99热全是精品| 最新的欧美精品一区二区| 97精品久久久久久久久久精品| 国产女主播在线喷水免费视频网站| 中文字幕精品免费在线观看视频| 色94色欧美一区二区| 国产成人免费无遮挡视频| 深夜精品福利| 欧美久久黑人一区二区| 人妻一区二区av| av天堂久久9| 亚洲第一青青草原| 丰满迷人的少妇在线观看| 欧美老熟妇乱子伦牲交| 两个人看的免费小视频| 精品一区二区三卡| 美女主播在线视频| a 毛片基地| 波野结衣二区三区在线| 欧美大码av| 欧美av亚洲av综合av国产av| 制服诱惑二区| 9热在线视频观看99| 精品久久久久久久毛片微露脸 | 老熟女久久久| 亚洲精品日韩在线中文字幕| 欧美 亚洲 国产 日韩一| 成人亚洲精品一区在线观看| 亚洲情色 制服丝袜| 新久久久久国产一级毛片| 国产色视频综合| www日本在线高清视频| 99re6热这里在线精品视频| av又黄又爽大尺度在线免费看|