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

    丹江口水庫上下游古菌優(yōu)勢菌群落結(jié)構(gòu)特征分析

    2020-06-07 10:09:22林亞萱黨晨原鐘思寧王佳文鄭彤倪晉仁
    關(guān)鍵詞:古菌丹江口水庫產(chǎn)甲烷

    林亞萱 黨晨原 鐘思寧 王佳文 鄭彤 倪晉仁

    丹江口水庫上下游古菌優(yōu)勢菌群落結(jié)構(gòu)特征分析

    林亞萱 黨晨原 鐘思寧 王佳文 鄭彤 倪晉仁?

    教育部水沙科學(xué)重點(diǎn)實(shí)驗(yàn)室, 北京大學(xué)環(huán)境工程系, 北京 100871; ?通信作者, E-mail:nijinren@iee.pku.edu.cn

    采用 16S rRNA Illumina Miseq 高通量測序技術(shù), 分析丹江口水庫庫區(qū)及漢江下游古菌物種組成, 并對大壩上、下游水體與沉積物中占優(yōu)勢的氨氧化古菌(Ammonia-oxidizing archaea, AOA)和產(chǎn)甲烷古菌(Methanogenic archaea)群落結(jié)構(gòu)進(jìn)行分析。結(jié)果表明, 優(yōu)勢種群落結(jié)構(gòu)組成受到水體與沉積物樣本差異的影響, 可由氨氧化古菌的好氧特性與產(chǎn)甲烷古菌的厭氧特性合理地解釋。網(wǎng)絡(luò)圖分析表明, 丹江口水庫上游氨氧化古菌與產(chǎn)甲烷古菌具有顯著的相關(guān)關(guān)系。受丹江口水庫運(yùn)行的影響, 大壩下游水體及沉積物中氨氧化古菌豐度皆比大壩上游少, 而沉積物中產(chǎn)甲烷古菌豐度較高, 二者間相關(guān)性不明顯。

    大壩; 氨氧化古菌; 產(chǎn)甲烷古菌; 高通量測序

    河流是重要的生態(tài)系統(tǒng), 在生態(tài)平衡和人類生活中起著至關(guān)重要的作用[1]。河流是地球陸地表面因流水作用而形成的典型地貌類型[2], 具有泄洪、供水、自凈、航運(yùn)、發(fā)電和景觀等多種功能[3]。河流中含有豐富的微生物, 它們是生物地球化學(xué)過程(如氮、硫和碳循環(huán))的關(guān)鍵參與者, 是各生態(tài)系統(tǒng)的主要生產(chǎn)者和消費(fèi)者[4], 對水質(zhì)有重要的影響[5]。為了實(shí)現(xiàn)水資源優(yōu)化配置, 利用水庫可以實(shí)現(xiàn)蓄水、防洪和發(fā)電等諸多功能, 但也會(huì)導(dǎo)致河流生境和生物的變化[6]。

    漢江是長江最大的支流, 水利資源豐富[7]。漢江上游山谷具有天然的地勢落差, 丹江口水庫的建設(shè)使庫區(qū)水位上升, 由于周邊土壤被淹沒, 導(dǎo)致土壤中養(yǎng)分進(jìn)入水體, 水中氮和磷的含量增加[8]。前人的研究表明, 大壩對生境以及細(xì)菌、魚類等群落結(jié)構(gòu)和多樣性造成影響[9–10]。Reyes-Gavilán 等[11]指出, 水壩建設(shè)會(huì)對魚類群落的遷徙造成影響, 從而影響整體多樣性。Wang 等[12]指出, 大壩會(huì)對浮游細(xì)菌的生物量和多樣性造成影響。另有研究表明, 大壩會(huì)導(dǎo)致沉積物物理化學(xué)性質(zhì)改變, 使得微生物群落結(jié)構(gòu)發(fā)生變化[13]。

    古菌中的氨氧化古菌(ammonia-oxidizing archa-ea, AOA)在河流生態(tài)系統(tǒng)的氮循環(huán)中發(fā)揮著重要作用[14]。例如, Zhou等[15]的研究發(fā)現(xiàn), 氨氧化古菌可以適應(yīng)極端環(huán)境, 并在極端環(huán)境中對全球氮循環(huán)發(fā)揮不可或缺的作用。產(chǎn)甲烷古菌(methanogenic ar-chaea)則是介導(dǎo)土壤和河流底泥中甲烷循環(huán)的重要功能物種[16], 廣泛分布在海洋沉積物[17]、淡水沉積物[18]、水稻田[19]、濕地[20]、人類和動(dòng)物的腸道[21]以及厭氧消化器等極端環(huán)境。有研究估計(jì), 河流對全球甲烷的貢獻(xiàn)將從每年 1.5Tg CH4升高到 26.8 Tg CH4, 相當(dāng)于濕地與湖泊甲烷排放量的 15%和40%[22–23]。鑒于氨氧化古菌及產(chǎn)甲烷古菌存在的廣泛性和重要性, 本文選擇漢江代表性河段進(jìn)行采樣, 重點(diǎn)了解漢江受丹江口水庫建設(shè)影響區(qū)域的古菌群落與優(yōu)勢種, 揭示大壩對研究區(qū)域古菌優(yōu)勢種的影響。

    1 研究區(qū)域與研究方法

    1.1 研究區(qū)域及樣點(diǎn)布設(shè)

    漢江全長 1532km, 流域面積有 15.9 萬 km2[24]。根據(jù)漢江的國家水文站點(diǎn)分布情況, 設(shè)置 6 個(gè)采樣點(diǎn), 從上游到下游分別是白河、陶岔、丹江口壩下、襄陽、仙桃和集家嘴(圖 1), 于 2014 年春季(3月)和秋季(10月)同步采集水樣與沉積物樣本。

    1.2 樣品采集與處理

    參照中華人民共和國環(huán)境保護(hù)行業(yè)標(biāo)準(zhǔn) HJ/T- 2002《地表水和污水檢測技術(shù)規(guī)范》中的等比例混合水樣采集法[25], 將采集的水樣在 24 小時(shí)內(nèi)用孔徑為 0.22μm 的聚碳酸酯膜(Millipore, 美國)進(jìn)行過濾。過濾水樣的濾膜與沉積物樣本都置于干冰中, 運(yùn)輸至實(shí)驗(yàn)室, 貯藏在?80oC冰箱中備用。

    水樣 DNA 提取使用液氮反復(fù)凍融法來裂解細(xì)胞, 后續(xù)與沉積物樣本一致, 依據(jù) FastDNA SPIN Kit for Soil (MP Biomedicals, 美國)試劑盒協(xié)議步驟進(jìn)行。

    1.3 數(shù)據(jù)處理與分析

    采用 524F10extF/Arch958RmodR[26]引物對 16S rRNA 進(jìn)行 3 次 PCR 重復(fù)擴(kuò)增。本研究采用 MiSeq高通量測序技術(shù), 以 97% 的相似性閾值進(jìn)行 OUT 聚類, 后與 Silva 數(shù)據(jù)庫比對, 獲得 OTU 物種表格。

    使用 R 語言軟件(R version 3.5.1)的 pheatmap, vegan, ggplot2 和 igraph程序包, 分別對數(shù)據(jù)進(jìn)行熱圖(Heatmap)分析、非度量多維標(biāo)定 (non-metric multi-dimensional scaling, NMDS)分析、ANOSIM分析、方差分析和網(wǎng)絡(luò)關(guān)系(Network)圖分析。其中, NMDS 使用 Bray-Curtis 距離, 方差分析以< 0.05 作為篩選標(biāo)準(zhǔn)。Network 以= 0.85 作為閾值, 做 Spearman 相似性分析, 使用Gephi 作圖。

    2 結(jié)果與討論

    2.1 古菌及優(yōu)勢種的群落組成

    通過 16S rRNA 高通量測序技術(shù), 分析漢江生態(tài)系統(tǒng)中古菌群落組成。此次調(diào)查共得到 2908 個(gè)OTU, reads 數(shù) 966898 條。所有樣品共檢測出 21 個(gè)門, 相對豐度大于 0.01%的共有 12 個(gè)門以及一個(gè)未分類微生物, 分別為(奇古菌門, 53.9%)、(廣古菌門, 26%)、(深古菌門, 9.5%)、(0.7%)、(洛基古菌, 0.8%)、(1.1%)、(0.6%)、(0.4%)、(0.3%)、(0.3%)、(0.3%)、(0.3%)及unclassified (5.3%), 如圖 2 所示。奇古菌門和廣古菌門是漢江流域的主要優(yōu)勢菌門, 數(shù)量占有絕對優(yōu)勢, 但兩者在水樣和泥樣中的相對豐度有顯著區(qū)別: 奇古菌門在水體樣品中的相對豐度(超過 50%)都比在沉積物中高(約占 20%); 廣古菌門是漢江流域沉積物中古菌豐度最高的門(43.5% (春泥)和 42.8%(秋泥)), 也是浮游古菌豐度第二高的優(yōu)勢門(23.5% (春水)和4.7% (秋水))。

    奇古菌門中以氨氧化古菌(87.90%)為主, 廣古菌門內(nèi)物種豐富, 包括產(chǎn)甲烷菌(65.34%)、甲烷氧化菌(7.64%)與嗜鹽菌(0.02%)等, 由此得出氨氧化古菌是奇古菌門的優(yōu)勢種, 產(chǎn)甲烷菌是廣古菌門的優(yōu)勢種的結(jié)論。漢江中的優(yōu)勢功能物種是氨氧化古菌和產(chǎn)甲烷古菌, 二者分別在漢江的水體與沉積物中起重要作用。發(fā)現(xiàn)的氨氧化古菌共有 4 種屬, 分別是(),(),()和(), 其中在氨氧化古菌中相對豐度最高。產(chǎn)甲烷古菌中是屬屬和屬的相對豐度最高。屬和屬是已知產(chǎn)甲烷古菌中僅有的兩種可以使用乙酸作為底物的屬[27]。研究表明, 在淡水環(huán)境中是乙酸型的產(chǎn)甲烷古菌占主導(dǎo), 其中乙酸型和氫營養(yǎng)產(chǎn)甲烷古菌分別占 CH4產(chǎn)量的 70%和 30%[28–29]。

    漢江上游的丹江口水壩是南水北調(diào)中線工程的調(diào)水源頭。圖 3 比較古菌優(yōu)勢種在丹江口水壩上、下游的相對豐度差異, 其中 4 個(gè)區(qū)域分別是丹江口水壩上游水體、上游沉積物、下游水體和下游沉積物中古菌的相對豐度, 顏色的深淺代表相對豐度的高低。4 個(gè)分組中, 氨氧化古菌的屬、屬綱和綱以及產(chǎn)甲烷古菌的屬、屬和屬都是豐度較高的物種。屬和綱在水體中相對豐度高, 體現(xiàn)氨氧化古菌對好氧環(huán)境的偏好。屬與屬則在水體中相對豐度較低, 在沉積物中相對豐度較高, 顯示產(chǎn)甲烷古菌的厭氧特性。

    本研究中, 浮游古菌以好氧的氨氧化古菌為主, 沉積物古菌則以厭氧的產(chǎn)甲烷古菌為主, 這種現(xiàn)象可由氨氧化古菌的好氧特性與產(chǎn)甲烷古菌的厭氧特性來合理的解釋。雖然在水體與沉積物中也分別發(fā)現(xiàn)產(chǎn)甲烷古菌與氨氧化古菌, 但這是由于氨氧化古菌可以在低氧環(huán)境中生存[15]。另一方面, 產(chǎn)甲烷古菌也在好氧環(huán)境中被發(fā)現(xiàn)[30–32], 但有氧環(huán)境對產(chǎn)甲烷有抑制作用。

    分別篩選出壩上與壩下水體和沉積物中具有顯著性差異的優(yōu)勢物種(<0.05)。圖 4 顯示, 水庫下游浮游氨氧化古菌的豐度略少于上游, 而產(chǎn)甲烷古菌的豐度在上、下游的差別不顯著(尤其是屬和屬)。圖 5 顯示, 沉積物中氨氧化古菌的豐度在丹江口水壩下游減少, 表明沉積物中氨氧化古菌受到水壩運(yùn)行的影響, 而產(chǎn)甲烷古菌的豐度在下游比上游有明顯的提升, 可能與上、下游沉積物中有機(jī)質(zhì)含量的差異有關(guān)。

    2.2 古菌優(yōu)勢種群落結(jié)構(gòu)空間分布

    采用 NMDS 分析法, 分別分析大壩上、下游兩個(gè)優(yōu)勢菌群落結(jié)構(gòu)的時(shí)空格局(圖 6)。壩上與壩下的氨氧化古菌群落組成在沉積物與水體間有差異, 沉積物中的物種形成一個(gè)明顯的區(qū)塊, 浮游物種相對分散, 但也形成一個(gè)區(qū)塊, 表明漢江流域的氨氧化古菌群落結(jié)構(gòu)組成受到水體和沉積物樣本差異的影響。產(chǎn)甲烷古菌與氨氧化古菌的分析結(jié)果一致, 產(chǎn)甲烷古菌的群落結(jié)構(gòu)組成在壩上與壩下都受到樣品差異的影響。

    為驗(yàn)證結(jié)果的可靠性, 對 NMDS 分析結(jié)果進(jìn)行相似性分析(圖 7)。壩上氨氧化古菌的值為 0.5,值為 0.01; 壩下值為 0.5643,值為 0.001。壩上產(chǎn)甲烷古菌的值為 0.3056,值為 0.04; 壩下值為 0.3047,值為 0.007。值表示結(jié)果具有統(tǒng)計(jì)學(xué)意義, 表明優(yōu)勢種在丹江口水庫上、下游的樣品差異分組具有意義, 顯示氨氧化古菌和產(chǎn)甲烷古菌都受到水體與沉積物樣品差異的影響。

    對古菌優(yōu)勢種的時(shí)空格局分析結(jié)果表明, 雖然丹江口水庫在一定程度上影響古菌優(yōu)勢種在上、下游的相對豐度, 但并未顯著地影響古菌群落結(jié)構(gòu), 表明古菌對不同環(huán)境具有更廣泛的適應(yīng)性。溫度、pH和氧氣濃度等的變化是導(dǎo)致水體與沉積物中優(yōu)勢種群落結(jié)構(gòu)變化的關(guān)鍵因素, 例如,屬和屬常在海洋與河口沉積物中被發(fā)現(xiàn)[33],屬偏好酸性土壤,屬則主要在陸地土壤中被發(fā)現(xiàn)[34]。產(chǎn)甲烷古菌內(nèi)部物種同樣對環(huán)境有不同的需求, 如甲烷屬要求乙酸濃度高的生存環(huán)境,屬則可在乙酸濃度偏低的環(huán)境中生存[35]。產(chǎn)甲烷菌根據(jù)對溫度的要求不同, 可分為嗜熱菌(55~80°C 左右)、嗜溫菌(35°C 左右)以及嗜冷菌(低于 25°C)。例如, 研究表明低溫條件下產(chǎn)甲烷古菌以乙酸型的屬占主導(dǎo)[36], 嗜冷菌的屬物種不能在超過 20°C的地方生長[37]。

    圖5 沉積物中優(yōu)勢物種在大壩上游與下游的豐度

    2.3 古菌優(yōu)勢種之間的關(guān)系

    通過網(wǎng)絡(luò)關(guān)系圖(圖 8), 分析氨氧化古菌與產(chǎn)甲烷古菌在大壩上、下游物種之間的關(guān)系。丹江口水壩上游優(yōu)勢物種共生網(wǎng)絡(luò)的節(jié)點(diǎn)中, 氨氧化古菌和產(chǎn)甲烷古菌分別占 44.44%和 55.56%, 兩者之間的關(guān)系以正相關(guān)為主。上游優(yōu)勢物種的共生網(wǎng)絡(luò)(圖 8(a))可以分為 3 個(gè)部分: 氨氧化古菌的獨(dú)立部分、產(chǎn)甲烷古菌獨(dú)立部分以及兩者相互關(guān)聯(lián)的部分。兩者相互關(guān)聯(lián)的部分主要包括 OTU2577, OTU 4187, OTU1436, OTU4358, OTU2748, OTU2799, OTU2545, OTU2543, OTU3255, OTU3089, OTU 3469 和 OTU4227, 分別屬于氨氧化古菌的綱與, 以及產(chǎn)甲烷古菌的,,,,和。

    丹江口水壩下游氨氧化古菌與產(chǎn)甲烷古菌的比例分別是 47.22%和 52.78%。兩個(gè)優(yōu)勢種之間的關(guān)系以正相關(guān)為主, 負(fù)相關(guān)僅占總體的不到 1%。下游兩個(gè)優(yōu)勢種各自形成獨(dú)立的區(qū)塊(圖 8(b)), 彼此之間基本上沒有關(guān)聯(lián)。既屬于上游相互聯(lián)系的又在下游相關(guān)的 OTU 有 OTU1436, OTU2577, OUT2545, OTU3255, OTU3089 和 OTU4227, 分別屬于氨氧化古菌的綱以及產(chǎn)甲烷古菌的和。

    對比丹江口水庫上、下游物種間的關(guān)系, 可知丹江口水庫上游氨氧化古菌與產(chǎn)甲烷古菌具有顯著的相關(guān)關(guān)系, 而下游的相關(guān)性并不明顯。黃柳琴等[38]的研究表明, 雖然水庫運(yùn)行對氨氧化古菌群落結(jié)構(gòu)的影響并不明顯, 但會(huì)改變其生存方式, 如在三峽大壩上游的氨氧化古菌以附著為主, 而下游以浮游為主。這與我們的研究結(jié)果十分類似, 即壩上游氨氧化古菌與產(chǎn)甲烷古菌關(guān)聯(lián)密切, 但壩下游因河床的劇烈沖刷大大削弱了二者之間的聯(lián)系。

    3 結(jié)論

    本文基于 2014 年春季(4 月)和秋季(10 月)在受丹江口水庫影響的漢江河段監(jiān)測結(jié)果, 對古菌物種組成、群落結(jié)構(gòu)及優(yōu)勢種的相互關(guān)系進(jìn)行分析, 得到以下主要結(jié)論。

    1)漢江流域古菌物種的組成由奇古菌門和廣古菌門主導(dǎo), 水體中以奇古菌門中的氨氧化古菌為主, 沉積物中則以廣古菌門的產(chǎn)甲烷古菌為主, 這可用氨氧化古菌和產(chǎn)甲烷古菌分別適宜好氧和厭氧環(huán)境進(jìn)行很好的解釋。

    2)丹江口水庫對氨氧化古菌豐度的影響更主要地表現(xiàn)為對沉積物中古菌的影響。沉積物中產(chǎn)甲烷古菌的豐度在壩下游明顯增高, 可能與上、下游沉積物中有機(jī)質(zhì)含量的差異有關(guān)。大壩上、下游氨氧化古菌與產(chǎn)甲烷古菌的群落結(jié)構(gòu)組成受水體與沉積物樣本類型差異的影響更大。

    3)網(wǎng)絡(luò)圖分析表明, 丹江口水庫上游氨氧化古菌與產(chǎn)甲烷古菌具有顯著的相關(guān)性, 但是在大壩下游二者之間未呈現(xiàn)明顯的相關(guān)性, 推測與壩下河床沖刷造成物種生存方式的變化有關(guān)。

    [1] Vos W M D, Eekert M H A V, Smidt H, et al. Structure and function of microbial communities in river ecosystems in response to global environmental changes and pollution. ACM Communications in Computer Algebrahelp, 2004, 42(1/2): 55–56

    [2] 欒建國, 陳文祥. 河流生態(tài)系統(tǒng)的典型特征和服務(wù)功能. 人民長江, 2004, 35(9): 41–43

    [3] 倪晉仁, 劉元元. 河流健康診斷與生態(tài)修復(fù). 中國水利, 2006(13): 4–10

    [4] Fan L, Song C, Meng S, et al. Spatial distribution of planktonic bacterial and archaeal communities in the upper section of the tidal reach in Yangtze River. Scientific Reports, 2016, 6(1): 39147

    [5] Hahn M W. The microbial diversity of inland waters. Current Opinion in Biotechnology, 2006, 17(3): 256–261

    [6] 蔡玉鵬. 大型水利工程對長江中下游關(guān)鍵生態(tài)功能區(qū)影響研究[D]. 南京: 河海大學(xué), 2007

    [7] 邢雯慧, 王堅(jiān)紅, 張方偉, 等. 秋汛期影響漢江流域降雨的水汽通道特征研究. 人民長江, 2019, 50 (2): 1–6

    [8] 李雯, 逯建維, 王蜜. 南水北調(diào)工程對丹江口生態(tài)環(huán)境影響. 居舍, 2018(32): 151–151

    [9] Jansson R. The effect of dams on biodiversity [D]. Ume?: Ume? University, 2006

    [10] Liu T, Zhang A N, Wang J, et al. Integrated bio-geography of planktonic and sedimentary bacterial communities in the Yangtze River. Microbiome, 2018, 6(1): 16

    [11] Reyes-Gavilán F G, Garrido R, Nicieza A G, et al. Fish community variation along physical gradients in short streams of northern Spain and the disruptive effect of dams. Hydrobiologia, 1996, 321(2): 155–163

    [12] Wang P, Wang X, Wang C, et al. Shift in bacterio-plankton diversity and structure: influence of anthro-pogenic disturbances along the Yarlung Tsangpo River on the Tibetan Plateau, China. Scientific Reports, 2017, 7(1): 12529

    [13] Li J P, Dong S K, Yang Z F, et al. Effects of cascade hydropower dams on the structure and distribution of riparian and upland vegetation along the middle-lower Lancang-Mekong River. Forest Ecology & Manage-ment, 2012, 284(3): 251–259

    [14] 劉正輝. 氨氧化古菌及其對氮循環(huán)貢獻(xiàn)的研究進(jìn)展. 微生物學(xué)通報(bào), 2015, 42(4): 774–782

    [15] Zhou L, Wang S, Zou Y, et al. Species, abundance and function of ammonia-oxidizing Archaea in Inland Wa-ters across China. Scientific Reports, 2015, 5: 15969

    [16] 佘晨興, 仝川. 閩江口蘆葦沼澤濕地土壤產(chǎn)甲烷菌群落結(jié)構(gòu)的垂直分布. 生態(tài)學(xué)報(bào), 2012, 32(17): 5299–5308

    [17] Whiticar M J, Faber E, Schoell M. Biogenic methane formation in marine and freshwater environments: CO2reduction vs. acetate fermentation — isotope evidence. Geochimica et Cosmochimica Acta, 1986, 50(5): 693–709

    [18] Chan O C, Wolf M, Hepperle D, et al. Methanogenic archaeal community in the sediment of an artificially partitioned acidic bog lake. FEMS Microbiology Ecology, 2002, 42(1): 119–129

    [19] Sakai S, Imachi H, Hanada S, et al.gen. nov. sp. nov. a methane-producing archaeon, the first isolate of the lineage “Rice Cluster I”, and proposal of the new archaeal orderord. nov. International Journal of Systematic and Evolutionary Microbiology, 2008, 58(4): 929–936

    [20] Angle J C, Morin T H, Solden L M, et al. Metha-nogenesis in oxygenated soils is a substantial fraction of wetland methane emissions. Nature Commu-nications, 2017, 8(1): 1567

    [21] Pester M, Brune A. Hydrogen is the central free intermediate during lignocellulose degradation by termite gut symbionts. ISME Journal, 2007, 1(6): 551–565

    [22] Bastviken D, Tranvik L J, Downing J A, et al. Fresh-water methane emissions offset the continental car-bon sink. Science, 2011, 331: 50–50

    [23] Stanley E H, Casson N J, Christel S T, et al. The ecology of methane in streams and rivers: patterns, controls, and global significance. Ecological Mono-graphs, 2016, 86(2): 146–171

    [24] 況琪軍, 譚渝云, 萬登榜, 等. 漢江中下游江段藻類現(xiàn)狀調(diào)查及“水華”成因分析. 長江流域資源與環(huán)境, 2000, 9(1): 63–70

    [25] 國家環(huán)境保護(hù)局. 地表水和污水監(jiān)測技術(shù)規(guī)范. HJ/T91-2002. 北京: 中國環(huán)境科學(xué)出版社, 2003

    [26] Pires A C C, Cleary D F R, Almeida A, et al. Denaturing gradient gel electrophoresis and barcoded pyrosequencing reveal unprecedented archaeal diver-sity in mangrove sediment and rhizosphere samples. Applied and Environmental Microbiology, 2012, 78 (16): 5520–5528

    [27] Liu Y, Whitman W B. Metabolic, phylogenetic, and ecological diversity of the methanogenic archaea. Annals of the New York Academy of Sciences, 2010, 1125(1): 171–189

    [28] Whiticar M J. Carbon and hydrogen isotope syste-matics of bacterial formation and oxidation of me-thane. Chemical Geology, 1999, 161(1/2/3): 291–314

    [29] Conrad R. Contribution of hydrogen to methane pro-duction and control of hydrogen concentrations in methanogenic soils and sediments. Fems Microbiology Ecology, 1999, 28(3): 193–202

    [30] Grossart H P, Frindte K, Dziallas C, et al. Microbial methane production in oxygenated water column of an oligotrophic lake. Proceedings of the National Aca-demy of Sciences, 2011, 108(49): 19657–19661

    [31] Bogard M J, Del Giorgio P A, Boutet L, et al. Oxic water column methanogenesis as a major component of aquatic CH4fluxes. Nature Communications, 2014, 5: 5350

    [32] Tang K W, Mcginnis D F, Ionescu D, et al. Methane production in oxic Lake waters potentially increases aquatic methane flux to air. Environmental Science & Technology Letters, 2016, 3(6): 227–233

    [33] Mosier A C, Allen E E, Kim M, et al. Genome sequence of “Nitrosoarchaeum limnia” BG20, a low-salinity ammonia-oxidizing archaeon from the San Francisco Bay estuary. American So-ciety for Microbiology, 2012, 194(8): 2119–2120

    [34] Tourna M, Stieglmeier M, Spang A, et al. Nitroso-sphaera viennensis, an ammonia oxidizing archaeon from soil. Proceedings of the National Academy of Sciences of the United States of America, 2011, 108 (20): 8420–8425

    [35] Conrad R. Microbial ecology of methanogens and methanotrophs. Advances in Agronomy, 2007, 96(7): 1–63

    [36] Grosskopf R. Diversity and structure of the metha-nogenic community in anoxic rice paddy soil mic-rocosms as examined by cultivation and direct 16S rRNA gene sequence retrieval. Applied and Environ-mental Microbiology, 1998, 64(3): 960–969

    [37] Franzmann P D, Liu Y, Balkwill D L, et al.sp. nov. a psychrophilic, H2-using methanogen from Ace Lake, Antarctica. International Journal of Systematic Bacteriology, 1997, 47(4): 1068–1072

    [38] 黃柳琴, 蔣宏忱, 王尚, 等. 長江三峽大壩兩側(cè)水體中氨氧化微生物種群結(jié)構(gòu)分析. 微生物學(xué)報(bào), 2011, 51(1): 83–90

    Community Characteristics of Dominant Archaea before and after the Danjiangkou Dam

    LIN Yaxuan,DANG Chenyuan,ZHONG Sining,WANG Jiawen,ZHENG Tong,NI Jinren?

    Key Laboratory of Water and Sediment Sciences (MOE), Department of Environmental Engineering, Peking University, Beijing 100871;? Corresponding author, E-mail: nijinren@iee.pku.edu.cn

    The 16S rRNA Illumina Miseq high-throughput sequencing technology was used to analyze the species composition of archaea in the Danjiangkou reservoir area and the downstream Hanjiang River. The structure of ammonia-oxidizing archaea (AOA) and methanogenic archaeathat dominated the water and sediments in the upstream and downstream dam was analyzed. The results showed that the structure of the dominant species was affected by the samples difference, which could be explained by the aerobic preference of the ammonia-oxidizing archaea and the anaerobic properties of the methanogenic archaea. According to the network analysis, there was a significant correlation between the ammonia-oxidizing archaea and the methanogenic archaea in the upstream of the Danjiangkou reservoir. Because of the operation of the Danjiangkou reservoir, the abundance of ammonia-oxidizing archaea in the water and sediments of downstream of the dam was reduced, and the abundance of methanogenic archaea in the sediments was relatively high. The relation between them was not obvious.

    dam; ammonia-oxidizing archaea; methane-oxidizing archaea; high-throughput sequencing

    10.13209/j.0479-8023.2020.034

    國家自然科學(xué)基金(51539001)資助

    2019–04–18;

    2019–06–25

    猜你喜歡
    古菌丹江口水庫產(chǎn)甲烷
    不同pH和氧氣條件下土壤古菌與海洋古菌的競爭適應(yīng)機(jī)制*
    變油為氣,“榨干”廢棄油田
    大自然探索(2022年5期)2022-07-11 03:10:33
    海洋古菌
    丹江口水庫的水資源管理與保護(hù)概述
    水庫移民農(nóng)地流轉(zhuǎn)與生計(jì)恢復(fù)問題及對策研究
    基于ETM+數(shù)據(jù)的水體信息提取
    綠色科技(2017年10期)2017-07-05 12:18:01
    丹江口水庫河南外遷農(nóng)村移民安置效果研究
    科技資訊(2017年7期)2017-05-06 11:42:44
    零級(jí)動(dòng)力學(xué)模型產(chǎn)甲烷量估算值與實(shí)測值的比較研究
    普通一級(jí)模型產(chǎn)甲烷量估算
    東太平洋海隆深海熱液區(qū)沉積物古菌多樣性分析
    久久人人爽av亚洲精品天堂| 国产亚洲一区二区精品| 热99国产精品久久久久久7| 我的亚洲天堂| 在线国产一区二区在线| 啦啦啦视频在线资源免费观看| 久久午夜综合久久蜜桃| 女警被强在线播放| 国产av又大| 国产在线一区二区三区精| 国产精品美女特级片免费视频播放器 | 久久久久久久国产电影| 黄色丝袜av网址大全| 制服诱惑二区| 超碰97精品在线观看| 黄色视频,在线免费观看| 校园春色视频在线观看| 久久天堂一区二区三区四区| 国产深夜福利视频在线观看| 国产免费男女视频| 丁香六月欧美| 亚洲五月天丁香| 久久久久久久久免费视频了| 欧美成狂野欧美在线观看| 一夜夜www| 亚洲国产看品久久| 国产一区在线观看成人免费| 免费在线观看日本一区| 亚洲视频免费观看视频| 国产深夜福利视频在线观看| 老熟女久久久| 女人高潮潮喷娇喘18禁视频| 天天躁狠狠躁夜夜躁狠狠躁| av天堂在线播放| 国产野战对白在线观看| 国产欧美日韩一区二区三| 一区二区日韩欧美中文字幕| 久久热在线av| 黑人巨大精品欧美一区二区蜜桃| 人人妻人人添人人爽欧美一区卜| 亚洲av片天天在线观看| 看黄色毛片网站| 老司机深夜福利视频在线观看| 精品福利观看| 亚洲专区国产一区二区| 欧美不卡视频在线免费观看 | 最新的欧美精品一区二区| 中国美女看黄片| 午夜免费观看网址| 韩国av一区二区三区四区| 1024香蕉在线观看| av天堂在线播放| 一区二区三区激情视频| 国产成人欧美| 桃红色精品国产亚洲av| 久久这里只有精品19| 一区二区三区精品91| 欧美日韩一级在线毛片| 精品久久蜜臀av无| 久久青草综合色| 国产乱人伦免费视频| 韩国精品一区二区三区| videos熟女内射| 好男人电影高清在线观看| 亚洲av成人不卡在线观看播放网| 巨乳人妻的诱惑在线观看| 欧美黑人欧美精品刺激| 久久久久久久午夜电影 | 久久国产精品大桥未久av| 久久天躁狠狠躁夜夜2o2o| 亚洲av熟女| 五月开心婷婷网| 亚洲精品成人av观看孕妇| 国产精品久久久久久人妻精品电影| 久久精品国产亚洲av高清一级| 亚洲九九香蕉| 电影成人av| 麻豆国产av国片精品| 久久久国产成人免费| 亚洲人成电影观看| 久久中文字幕人妻熟女| 岛国毛片在线播放| www.999成人在线观看| 露出奶头的视频| cao死你这个sao货| 99精品久久久久人妻精品| 久久草成人影院| 欧美最黄视频在线播放免费 | 国产不卡一卡二| 欧美日韩中文字幕国产精品一区二区三区 | 中文亚洲av片在线观看爽 | 国产精品香港三级国产av潘金莲| 搡老岳熟女国产| 精品一品国产午夜福利视频| 女人高潮潮喷娇喘18禁视频| 黑人巨大精品欧美一区二区蜜桃| www.自偷自拍.com| 侵犯人妻中文字幕一二三四区| 亚洲成人免费av在线播放| 欧美黄色片欧美黄色片| 69av精品久久久久久| 一边摸一边抽搐一进一小说 | 国产精品一区二区精品视频观看| 在线观看舔阴道视频| 亚洲专区国产一区二区| 久久久国产欧美日韩av| 精品人妻1区二区| 黄色视频,在线免费观看| 亚洲欧美日韩另类电影网站| 国产99久久九九免费精品| 免费在线观看视频国产中文字幕亚洲| 亚洲五月色婷婷综合| 欧美+亚洲+日韩+国产| 国产一区在线观看成人免费| 热99国产精品久久久久久7| 老司机影院毛片| 黄频高清免费视频| 国产黄色免费在线视频| 久久九九热精品免费| 国产成人av教育| 亚洲综合色网址| 免费不卡黄色视频| 老司机靠b影院| 国产成人精品久久二区二区免费| av不卡在线播放| 欧美亚洲日本最大视频资源| 久久久国产欧美日韩av| 人妻一区二区av| 十八禁人妻一区二区| av片东京热男人的天堂| 亚洲黑人精品在线| 在线观看免费午夜福利视频| xxx96com| 欧美 亚洲 国产 日韩一| 又紧又爽又黄一区二区| 日韩 欧美 亚洲 中文字幕| 亚洲熟妇熟女久久| 久久久久精品人妻al黑| 少妇裸体淫交视频免费看高清 | av福利片在线| 大型av网站在线播放| 久久久久久久久久久久大奶| 亚洲色图 男人天堂 中文字幕| av国产精品久久久久影院| 自拍欧美九色日韩亚洲蝌蚪91| 黑人猛操日本美女一级片| 我的亚洲天堂| 窝窝影院91人妻| 久久久久国产精品人妻aⅴ院 | 一夜夜www| av福利片在线| 欧美一级毛片孕妇| 免费av中文字幕在线| 男女午夜视频在线观看| av在线播放免费不卡| a级片在线免费高清观看视频| 欧美午夜高清在线| 久久久久视频综合| 丝袜美足系列| 91精品国产国语对白视频| 亚洲成人免费电影在线观看| 热99久久久久精品小说推荐| 夜夜爽天天搞| 色婷婷av一区二区三区视频| 美女福利国产在线| 国产亚洲av高清不卡| 亚洲专区国产一区二区| 精品国产亚洲在线| 亚洲伊人色综图| 老汉色av国产亚洲站长工具| 十八禁人妻一区二区| 亚洲国产中文字幕在线视频| 精品国内亚洲2022精品成人 | 成人影院久久| 日本a在线网址| 日韩大码丰满熟妇| 人人妻人人添人人爽欧美一区卜| 国产无遮挡羞羞视频在线观看| 80岁老熟妇乱子伦牲交| 久久性视频一级片| 色尼玛亚洲综合影院| 国产亚洲精品一区二区www | 国产亚洲欧美在线一区二区| 黄色片一级片一级黄色片| 日日爽夜夜爽网站| 一级毛片女人18水好多| 精品久久久精品久久久| 侵犯人妻中文字幕一二三四区| 一区二区日韩欧美中文字幕| 久久精品人人爽人人爽视色| 欧美色视频一区免费| 精品国内亚洲2022精品成人 | 80岁老熟妇乱子伦牲交| 亚洲五月天丁香| 日韩欧美在线二视频 | 大香蕉久久网| 在线观看舔阴道视频| 精品一区二区三区av网在线观看| 欧美乱色亚洲激情| 女性被躁到高潮视频| 国产无遮挡羞羞视频在线观看| 国产成人精品在线电影| 桃红色精品国产亚洲av| 欧美中文综合在线视频| 黄色a级毛片大全视频| 亚洲欧美精品综合一区二区三区| 新久久久久国产一级毛片| 国产亚洲欧美在线一区二区| 无人区码免费观看不卡| 精品国内亚洲2022精品成人 | 日韩成人在线观看一区二区三区| 国产有黄有色有爽视频| 美女扒开内裤让男人捅视频| 免费在线观看视频国产中文字幕亚洲| 精品人妻1区二区| 黄色视频不卡| 亚洲黑人精品在线| svipshipincom国产片| 美女国产高潮福利片在线看| av在线播放免费不卡| 久久久久久久精品吃奶| 黄色片一级片一级黄色片| 成人精品一区二区免费| 久久久久久人人人人人| 亚洲av成人av| 51午夜福利影视在线观看| 亚洲精品久久成人aⅴ小说| 伊人久久大香线蕉亚洲五| 国产一区二区三区在线臀色熟女 | 少妇被粗大的猛进出69影院| 啪啪无遮挡十八禁网站| xxx96com| 99国产极品粉嫩在线观看| 亚洲五月天丁香| 人妻丰满熟妇av一区二区三区 | 亚洲精品中文字幕一二三四区| 亚洲国产中文字幕在线视频| 欧美成人午夜精品| 无遮挡黄片免费观看| 久久久久久久精品吃奶| 亚洲情色 制服丝袜| 午夜视频精品福利| 叶爱在线成人免费视频播放| 亚洲成a人片在线一区二区| 国产精品98久久久久久宅男小说| 在线观看免费视频网站a站| 国产精品 国内视频| 国产97色在线日韩免费| 韩国av一区二区三区四区| 精品久久久久久,| 免费女性裸体啪啪无遮挡网站| 婷婷丁香在线五月| 黄网站色视频无遮挡免费观看| 美女高潮到喷水免费观看| 18禁裸乳无遮挡动漫免费视频| 久热爱精品视频在线9| 国产一卡二卡三卡精品| 久久久国产成人免费| 99久久99久久久精品蜜桃| 国产在线观看jvid| 交换朋友夫妻互换小说| 两个人看的免费小视频| 亚洲 国产 在线| 国内毛片毛片毛片毛片毛片| 欧美丝袜亚洲另类 | 999精品在线视频| 女人爽到高潮嗷嗷叫在线视频| 涩涩av久久男人的天堂| 天天操日日干夜夜撸| 人人妻人人添人人爽欧美一区卜| 嫩草影视91久久| 69av精品久久久久久| 正在播放国产对白刺激| 午夜视频精品福利| 欧美在线一区亚洲| 怎么达到女性高潮| 自线自在国产av| 亚洲中文av在线| 精品无人区乱码1区二区| 免费日韩欧美在线观看| 国产xxxxx性猛交| 男女下面插进去视频免费观看| 午夜成年电影在线免费观看| 黄网站色视频无遮挡免费观看| 欧美日韩中文字幕国产精品一区二区三区 | 啦啦啦 在线观看视频| 免费不卡黄色视频| 99久久99久久久精品蜜桃| 国产不卡一卡二| 久久精品亚洲av国产电影网| 美女福利国产在线| 一a级毛片在线观看| 波多野结衣一区麻豆| 久久精品亚洲精品国产色婷小说| 日韩欧美一区二区三区在线观看 | 人妻丰满熟妇av一区二区三区 | 久久久久久久精品吃奶| 黄色怎么调成土黄色| 国产男靠女视频免费网站| 一级,二级,三级黄色视频| 丁香六月欧美| 亚洲精品久久成人aⅴ小说| 欧美色视频一区免费| 久久精品亚洲熟妇少妇任你| 午夜免费观看网址| 视频区图区小说| 99久久99久久久精品蜜桃| 国产高清激情床上av| 亚洲一卡2卡3卡4卡5卡精品中文| 91在线观看av| 他把我摸到了高潮在线观看| 日本撒尿小便嘘嘘汇集6| 在线十欧美十亚洲十日本专区| 天天躁夜夜躁狠狠躁躁| 欧美在线一区亚洲| 国产1区2区3区精品| 亚洲少妇的诱惑av| 国产午夜精品久久久久久| 国产精品免费大片| 建设人人有责人人尽责人人享有的| 一个人免费在线观看的高清视频| 亚洲avbb在线观看| 亚洲欧美日韩另类电影网站| 亚洲av美国av| 欧美黄色淫秽网站| 亚洲五月色婷婷综合| 91成年电影在线观看| x7x7x7水蜜桃| a级毛片黄视频| 久久香蕉国产精品| 中出人妻视频一区二区| 老司机靠b影院| 满18在线观看网站| 国内久久婷婷六月综合欲色啪| 悠悠久久av| 亚洲一码二码三码区别大吗| 成人影院久久| 国产成人系列免费观看| 男人的好看免费观看在线视频 | 国产精品 欧美亚洲| 18禁观看日本| 狠狠狠狠99中文字幕| 久久精品国产99精品国产亚洲性色 | 亚洲国产欧美日韩在线播放| 欧美激情高清一区二区三区| 黄片播放在线免费| 精品亚洲成国产av| 91国产中文字幕| 手机成人av网站| 国产三级黄色录像| 免费久久久久久久精品成人欧美视频| 日本黄色视频三级网站网址 | 丝袜美足系列| 欧美激情 高清一区二区三区| 日韩一卡2卡3卡4卡2021年| 亚洲一区二区三区不卡视频| 老熟女久久久| 亚洲一区二区三区欧美精品| netflix在线观看网站| 一区福利在线观看| 亚洲熟女精品中文字幕| www日本在线高清视频| a级片在线免费高清观看视频| 首页视频小说图片口味搜索| 亚洲精品国产区一区二| 久久99一区二区三区| 久久青草综合色| 欧美日韩成人在线一区二区| 国精品久久久久久国模美| 精品一品国产午夜福利视频| 视频区图区小说| 亚洲三区欧美一区| 中文字幕另类日韩欧美亚洲嫩草| 夜夜躁狠狠躁天天躁| 久久狼人影院| 日韩有码中文字幕| 少妇裸体淫交视频免费看高清 | 国产精品乱码一区二三区的特点 | 窝窝影院91人妻| 日韩欧美一区视频在线观看| 久久精品成人免费网站| 99久久综合精品五月天人人| 国产免费男女视频| 久久精品aⅴ一区二区三区四区| 大型av网站在线播放| 亚洲精品av麻豆狂野| 热99久久久久精品小说推荐| 女人精品久久久久毛片| 999精品在线视频| 久久久久久久久久久久大奶| 热re99久久国产66热| 少妇被粗大的猛进出69影院| 首页视频小说图片口味搜索| 又黄又爽又免费观看的视频| 天堂√8在线中文| 大片电影免费在线观看免费| 欧美+亚洲+日韩+国产| 中亚洲国语对白在线视频| netflix在线观看网站| 日本精品一区二区三区蜜桃| 久久久国产一区二区| 日韩欧美国产一区二区入口| 一区二区三区精品91| 99国产精品一区二区三区| 亚洲精品粉嫩美女一区| 久久午夜综合久久蜜桃| 国产又色又爽无遮挡免费看| 久久这里只有精品19| 亚洲精品在线美女| av网站免费在线观看视频| 91在线观看av| 捣出白浆h1v1| av一本久久久久| 国产男靠女视频免费网站| 欧美人与性动交α欧美精品济南到| 亚洲一区中文字幕在线| 天天添夜夜摸| 欧美最黄视频在线播放免费 | 男人舔女人的私密视频| 在线观看免费日韩欧美大片| 又紧又爽又黄一区二区| 亚洲av日韩在线播放| 啦啦啦视频在线资源免费观看| 在线免费观看的www视频| 成人18禁在线播放| 国产成人系列免费观看| 免费在线观看黄色视频的| 精品国产超薄肉色丝袜足j| 交换朋友夫妻互换小说| 中文字幕精品免费在线观看视频| 亚洲中文av在线| 精品一区二区三区av网在线观看| 两个人免费观看高清视频| 欧美成人免费av一区二区三区 | 久久人妻熟女aⅴ| 9色porny在线观看| 精品一品国产午夜福利视频| 老司机午夜福利在线观看视频| 婷婷精品国产亚洲av在线 | 自线自在国产av| 男人操女人黄网站| 9热在线视频观看99| 中文字幕人妻丝袜制服| 美女福利国产在线| 老汉色av国产亚洲站长工具| 精品久久久久久久久久免费视频 | 国产精品影院久久| 久9热在线精品视频| 亚洲中文字幕日韩| 国产成人精品无人区| 久久精品亚洲熟妇少妇任你| 校园春色视频在线观看| 黑人巨大精品欧美一区二区mp4| 日本a在线网址| 国产成+人综合+亚洲专区| 午夜福利在线免费观看网站| 99精品在免费线老司机午夜| 欧美黑人欧美精品刺激| 丰满迷人的少妇在线观看| 在线观看免费视频网站a站| 狂野欧美激情性xxxx| 精品亚洲成a人片在线观看| 美女 人体艺术 gogo| 久久久久久人人人人人| 一级作爱视频免费观看| 999久久久精品免费观看国产| 免费看十八禁软件| 欧美 亚洲 国产 日韩一| 成人免费观看视频高清| 国精品久久久久久国模美| 亚洲熟妇中文字幕五十中出 | 精品欧美一区二区三区在线| 久久人人97超碰香蕉20202| 国产精品偷伦视频观看了| 免费在线观看视频国产中文字幕亚洲| 久久国产精品大桥未久av| 亚洲av片天天在线观看| 欧美日本中文国产一区发布| 丰满迷人的少妇在线观看| 国产不卡一卡二| 久久这里只有精品19| 999精品在线视频| 国产黄色免费在线视频| 老司机深夜福利视频在线观看| 90打野战视频偷拍视频| 法律面前人人平等表现在哪些方面| videos熟女内射| 99久久国产精品久久久| 建设人人有责人人尽责人人享有的| 亚洲欧美一区二区三区黑人| 日韩人妻精品一区2区三区| 亚洲一区二区三区不卡视频| 人人妻人人澡人人爽人人夜夜| 50天的宝宝边吃奶边哭怎么回事| 一进一出好大好爽视频| 欧美中文综合在线视频| 久久久国产欧美日韩av| 久久久精品国产亚洲av高清涩受| 免费在线观看影片大全网站| 香蕉国产在线看| 欧美精品亚洲一区二区| 久久99一区二区三区| 亚洲av成人一区二区三| 中文字幕人妻丝袜一区二区| 国产精品永久免费网站| 国产成人系列免费观看| 日本vs欧美在线观看视频| 国产一区二区三区视频了| 老鸭窝网址在线观看| 美女国产高潮福利片在线看| 国内毛片毛片毛片毛片毛片| 国产成人精品久久二区二区免费| 19禁男女啪啪无遮挡网站| 成人18禁高潮啪啪吃奶动态图| 午夜福利乱码中文字幕| 两人在一起打扑克的视频| 国产欧美日韩一区二区三区在线| 久久午夜亚洲精品久久| 欧美成人午夜精品| 99国产精品一区二区蜜桃av | 成年人免费黄色播放视频| 淫妇啪啪啪对白视频| 国产精品免费一区二区三区在线 | 在线十欧美十亚洲十日本专区| 国产av精品麻豆| 精品国产乱子伦一区二区三区| 久久天堂一区二区三区四区| 在线播放国产精品三级| 精品一品国产午夜福利视频| 欧美黄色淫秽网站| av线在线观看网站| av福利片在线| 欧美日韩亚洲高清精品| 男人的好看免费观看在线视频 | 免费高清在线观看日韩| 国产成人精品久久二区二区免费| 欧美日韩中文字幕国产精品一区二区三区 | 国产成人啪精品午夜网站| av免费在线观看网站| 亚洲人成77777在线视频| 中文字幕av电影在线播放| 亚洲精品粉嫩美女一区| 亚洲色图综合在线观看| 十八禁网站免费在线| 精品少妇久久久久久888优播| 在线天堂中文资源库| 日日夜夜操网爽| 大香蕉久久网| 日日爽夜夜爽网站| 亚洲人成77777在线视频| 王馨瑶露胸无遮挡在线观看| 天天添夜夜摸| 精品高清国产在线一区| 午夜福利免费观看在线| 搡老熟女国产l中国老女人| 中文字幕另类日韩欧美亚洲嫩草| av不卡在线播放| 国产伦人伦偷精品视频| 色播在线永久视频| 精品国产一区二区三区久久久樱花| 成人永久免费在线观看视频| 丝袜美腿诱惑在线| netflix在线观看网站| 国产有黄有色有爽视频| 国产精品国产高清国产av | 91在线观看av| 久久国产精品大桥未久av| 亚洲午夜理论影院| 国产av精品麻豆| 搡老熟女国产l中国老女人| 亚洲七黄色美女视频| 精品福利观看| 黄片大片在线免费观看| 黄网站色视频无遮挡免费观看| 波多野结衣一区麻豆| 久久久久久久国产电影| 老司机午夜十八禁免费视频| 每晚都被弄得嗷嗷叫到高潮| 久久久精品国产亚洲av高清涩受| 午夜免费成人在线视频| 极品人妻少妇av视频| 999久久久国产精品视频| x7x7x7水蜜桃| av天堂在线播放| 高清av免费在线| 一区二区日韩欧美中文字幕| xxx96com| 91精品三级在线观看| 啦啦啦视频在线资源免费观看| 亚洲中文av在线| 99精品欧美一区二区三区四区| 777米奇影视久久| 久久久精品免费免费高清| 国产片内射在线| 人人妻人人爽人人添夜夜欢视频| 黄色毛片三级朝国网站| 久久国产精品大桥未久av| 两性午夜刺激爽爽歪歪视频在线观看 | 极品人妻少妇av视频| 国产精品久久久av美女十八| 国产精品久久电影中文字幕 | 久久久久久久精品吃奶| 免费日韩欧美在线观看| 国产黄色免费在线视频| 久久婷婷成人综合色麻豆| 精品久久久久久电影网| 乱人伦中国视频| 狂野欧美激情性xxxx| av在线播放免费不卡| 午夜成年电影在线免费观看| 18禁国产床啪视频网站| 久久人人爽av亚洲精品天堂| 国产一区二区激情短视频| 一本综合久久免费| 久久久久久久久久久久大奶|