• <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ū)沉積物古菌多樣性分析
    最近2019中文字幕mv第一页| 日韩中文字幕视频在线看片 | 亚洲精品第二区| 欧美xxxx性猛交bbbb| 熟女人妻精品中文字幕| 成人18禁高潮啪啪吃奶动态图 | 久久青草综合色| 一区二区av电影网| 免费少妇av软件| av一本久久久久| 国产一区有黄有色的免费视频| 日韩不卡一区二区三区视频在线| 伦理电影大哥的女人| 久久国内精品自在自线图片| 欧美bdsm另类| av在线app专区| 一本久久精品| 少妇猛男粗大的猛烈进出视频| 亚洲av欧美aⅴ国产| 成人亚洲欧美一区二区av| 18+在线观看网站| 精品久久久噜噜| 亚洲国产av新网站| 一级毛片黄色毛片免费观看视频| 亚洲一级一片aⅴ在线观看| 免费av不卡在线播放| 人妻制服诱惑在线中文字幕| 欧美激情国产日韩精品一区| 国产白丝娇喘喷水9色精品| 精品国产一区二区三区久久久樱花 | 少妇人妻一区二区三区视频| 99精国产麻豆久久婷婷| 九九在线视频观看精品| 亚洲精华国产精华液的使用体验| 黄色日韩在线| 久久国产精品大桥未久av | 成人免费观看视频高清| 亚洲国产精品999| 久久久久国产网址| 国产黄频视频在线观看| 欧美日韩视频高清一区二区三区二| h视频一区二区三区| 精品一品国产午夜福利视频| 亚洲人成网站高清观看| 成年av动漫网址| 亚洲精品色激情综合| av免费观看日本| 久久久久精品久久久久真实原创| 自拍偷自拍亚洲精品老妇| 日本午夜av视频| 亚洲综合色惰| 高清av免费在线| 一级毛片久久久久久久久女| 夜夜骑夜夜射夜夜干| 午夜激情福利司机影院| 久久久色成人| 国产国拍精品亚洲av在线观看| 少妇被粗大猛烈的视频| 女人久久www免费人成看片| 多毛熟女@视频| 精品一区二区免费观看| 精品国产露脸久久av麻豆| a级毛片免费高清观看在线播放| 免费看不卡的av| 大又大粗又爽又黄少妇毛片口| 久久精品国产a三级三级三级| 日本欧美视频一区| 久久精品熟女亚洲av麻豆精品| 中文字幕亚洲精品专区| 少妇精品久久久久久久| 久久婷婷青草| 一个人看视频在线观看www免费| 99国产精品免费福利视频| 国产成人精品久久久久久| 亚洲综合色惰| 久久99热这里只频精品6学生| 精品久久国产蜜桃| 高清不卡的av网站| 观看美女的网站| 久久综合国产亚洲精品| 好男人视频免费观看在线| 亚洲欧美中文字幕日韩二区| 亚洲自偷自拍三级| 在线天堂最新版资源| 中国美白少妇内射xxxbb| 99re6热这里在线精品视频| 日本黄色日本黄色录像| 人人妻人人看人人澡| 欧美精品一区二区大全| 国产爱豆传媒在线观看| 最新中文字幕久久久久| 少妇人妻 视频| 精品少妇久久久久久888优播| 最近手机中文字幕大全| 男人爽女人下面视频在线观看| 一本一本综合久久| 国产精品久久久久久av不卡| 亚洲精品国产av成人精品| 狂野欧美激情性bbbbbb| 亚洲国产精品国产精品| 国产午夜精品一二区理论片| 日韩,欧美,国产一区二区三区| 亚洲欧美日韩东京热| 人人妻人人澡人人爽人人夜夜| 肉色欧美久久久久久久蜜桃| 2022亚洲国产成人精品| www.色视频.com| 777米奇影视久久| 日本av手机在线免费观看| 男人和女人高潮做爰伦理| 久久久色成人| 晚上一个人看的免费电影| videos熟女内射| 亚洲av二区三区四区| 国产午夜精品久久久久久一区二区三区| 夫妻午夜视频| 亚洲人与动物交配视频| 欧美日韩亚洲高清精品| 国产伦精品一区二区三区四那| 91精品一卡2卡3卡4卡| 亚洲成色77777| 精品一区二区免费观看| 亚洲精品日本国产第一区| 大又大粗又爽又黄少妇毛片口| 熟女电影av网| 国产欧美另类精品又又久久亚洲欧美| 能在线免费看毛片的网站| 亚洲国产毛片av蜜桃av| 亚洲激情五月婷婷啪啪| 嫩草影院新地址| 国产久久久一区二区三区| 日韩成人伦理影院| 深夜a级毛片| 精品久久久久久久末码| 亚洲av在线观看美女高潮| 一级毛片黄色毛片免费观看视频| 啦啦啦在线观看免费高清www| 毛片一级片免费看久久久久| 成人二区视频| 一个人看的www免费观看视频| 中文欧美无线码| 亚洲欧美日韩另类电影网站 | 能在线免费看毛片的网站| 日韩一本色道免费dvd| 99热6这里只有精品| 爱豆传媒免费全集在线观看| 亚洲熟女精品中文字幕| 高清毛片免费看| 精品人妻一区二区三区麻豆| 亚洲人成网站在线观看播放| 欧美区成人在线视频| 亚洲av日韩在线播放| 最新中文字幕久久久久| 伦理电影免费视频| 男人和女人高潮做爰伦理| 在线亚洲精品国产二区图片欧美 | 日韩av在线免费看完整版不卡| 亚洲国产av新网站| 久久久久久久久大av| 一本色道久久久久久精品综合| 日本vs欧美在线观看视频 | 纵有疾风起免费观看全集完整版| 91久久精品电影网| 自拍偷自拍亚洲精品老妇| 大片电影免费在线观看免费| 亚洲婷婷狠狠爱综合网| 久热久热在线精品观看| 久久婷婷青草| 身体一侧抽搐| 日本wwww免费看| 我要看日韩黄色一级片| 亚洲四区av| 丝瓜视频免费看黄片| 欧美激情国产日韩精品一区| 91午夜精品亚洲一区二区三区| 亚洲欧美成人综合另类久久久| 一级毛片 在线播放| 久久久久网色| 人妻一区二区av| 成人国产麻豆网| 亚洲精品日本国产第一区| 色5月婷婷丁香| 日韩不卡一区二区三区视频在线| 日本猛色少妇xxxxx猛交久久| 久久99热这里只有精品18| 在线观看国产h片| 国产精品国产三级国产av玫瑰| 成人国产av品久久久| 多毛熟女@视频| 国产精品熟女久久久久浪| 久久久久久久久久久丰满| 国产一区二区三区综合在线观看 | 又大又黄又爽视频免费| 国产精品精品国产色婷婷| 国产精品久久久久久精品电影小说 | 日韩电影二区| 国产精品99久久久久久久久| 国产黄色免费在线视频| 青春草亚洲视频在线观看| 日韩亚洲欧美综合| 少妇精品久久久久久久| 亚洲美女黄色视频免费看| 欧美成人a在线观看| 丝瓜视频免费看黄片| 人妻夜夜爽99麻豆av| 国产亚洲91精品色在线| 少妇 在线观看| 性高湖久久久久久久久免费观看| 国产成人91sexporn| 熟女人妻精品中文字幕| 亚洲国产日韩一区二区| 黑人高潮一二区| 成年免费大片在线观看| 黄色配什么色好看| 大片电影免费在线观看免费| 晚上一个人看的免费电影| 制服丝袜香蕉在线| 久久6这里有精品| videos熟女内射| 亚洲av中文av极速乱| 卡戴珊不雅视频在线播放| 国产色爽女视频免费观看| 你懂的网址亚洲精品在线观看| xxx大片免费视频| 天堂俺去俺来也www色官网| av天堂中文字幕网| 七月丁香在线播放| 青春草国产在线视频| 老司机影院毛片| 国产亚洲最大av| 亚洲精品国产av蜜桃| 中文天堂在线官网| 丰满少妇做爰视频| 国产成人一区二区在线| 日韩中文字幕视频在线看片 | 久久国产乱子免费精品| 亚洲欧美清纯卡通| 欧美日本视频| 国产精品偷伦视频观看了| 少妇裸体淫交视频免费看高清| 亚洲精品456在线播放app| 国产成人精品久久久久久| 日韩av免费高清视频| 一级毛片我不卡| 蜜桃在线观看..| 国产爱豆传媒在线观看| 亚洲内射少妇av| 狂野欧美激情性bbbbbb| 大陆偷拍与自拍| 丰满人妻一区二区三区视频av| 99热全是精品| 永久网站在线| 国产成人免费无遮挡视频| 国产色婷婷99| 女的被弄到高潮叫床怎么办| 国产成人精品婷婷| 91在线精品国自产拍蜜月| 黑人猛操日本美女一级片| 91久久精品国产一区二区成人| 国产成人精品福利久久| 中文字幕av成人在线电影| 国产精品一区二区在线观看99| 97超碰精品成人国产| 日本av免费视频播放| 国产av国产精品国产| 狠狠精品人妻久久久久久综合| 国产色婷婷99| 99精国产麻豆久久婷婷| 国产亚洲最大av| 亚洲国产成人一精品久久久| 身体一侧抽搐| 深爱激情五月婷婷| av女优亚洲男人天堂| 黄色欧美视频在线观看| 国产精品久久久久久久电影| 十分钟在线观看高清视频www | 久久国产精品男人的天堂亚洲 | 日本与韩国留学比较| 大片电影免费在线观看免费| 最近2019中文字幕mv第一页| 国产久久久一区二区三区| 国产精品一区二区在线观看99| av播播在线观看一区| 国产免费视频播放在线视频| 涩涩av久久男人的天堂| 大话2 男鬼变身卡| 男女边吃奶边做爰视频| 国产精品国产三级国产专区5o| 国产av一区二区精品久久 | 大码成人一级视频| 老司机影院成人| 国产伦精品一区二区三区四那| 男人舔奶头视频| 欧美成人一区二区免费高清观看| 久久亚洲国产成人精品v| 一本—道久久a久久精品蜜桃钙片| 日韩 亚洲 欧美在线| 熟女电影av网| 校园人妻丝袜中文字幕| 在线播放无遮挡| 亚洲av成人精品一二三区| 久久精品国产亚洲网站| 中文字幕av成人在线电影| 少妇猛男粗大的猛烈进出视频| 亚洲av中文字字幕乱码综合| 五月伊人婷婷丁香| 建设人人有责人人尽责人人享有的 | 欧美精品一区二区免费开放| 亚洲精品第二区| 伦理电影大哥的女人| 丰满迷人的少妇在线观看| 日本欧美视频一区| 丰满人妻一区二区三区视频av| 一级黄片播放器| 久久久久久久久久人人人人人人| 在线观看国产h片| 国产免费又黄又爽又色| 成人国产麻豆网| av女优亚洲男人天堂| 国产在线视频一区二区| 成人漫画全彩无遮挡| 青春草视频在线免费观看| 亚洲av.av天堂| 晚上一个人看的免费电影| 大话2 男鬼变身卡| 一本久久精品| 一个人看的www免费观看视频| 日日撸夜夜添| 亚洲不卡免费看| 国产永久视频网站| 成年人午夜在线观看视频| 三级经典国产精品| 亚洲aⅴ乱码一区二区在线播放| 成年女人在线观看亚洲视频| 亚洲av福利一区| av在线老鸭窝| 国产精品国产三级国产av玫瑰| 亚洲成人一二三区av| 在线观看av片永久免费下载| 日韩免费高清中文字幕av| 欧美激情极品国产一区二区三区 | 99视频精品全部免费 在线| 交换朋友夫妻互换小说| 婷婷色综合大香蕉| 视频中文字幕在线观看| 女的被弄到高潮叫床怎么办| 黄色日韩在线| 国产精品秋霞免费鲁丝片| 能在线免费看毛片的网站| 美女中出高潮动态图| 黄片无遮挡物在线观看| 九九久久精品国产亚洲av麻豆| 五月开心婷婷网| 国产精品av视频在线免费观看| 亚洲欧美精品专区久久| 18禁在线播放成人免费| 精品久久久久久久久av| 国产爽快片一区二区三区| 欧美区成人在线视频| 99热网站在线观看| 狂野欧美白嫩少妇大欣赏| 日韩电影二区| 欧美日韩国产mv在线观看视频 | 99久久精品国产国产毛片| 91狼人影院| 色哟哟·www| av线在线观看网站| 日本与韩国留学比较| 亚洲精品色激情综合| 国产精品久久久久久精品古装| 九九在线视频观看精品| 国产精品一及| 在线观看免费高清a一片| 久久久a久久爽久久v久久| 欧美区成人在线视频| 精品午夜福利在线看| 欧美极品一区二区三区四区| 国产伦在线观看视频一区| 国产精品久久久久久久久免| 亚洲精品456在线播放app| 欧美老熟妇乱子伦牲交| 日本黄大片高清| 亚洲精品国产成人久久av| 亚洲国产成人一精品久久久| 久久久久久九九精品二区国产| 亚洲综合色惰| 久久精品国产亚洲网站| 一个人看的www免费观看视频| 男女免费视频国产| 99热这里只有是精品在线观看| 欧美变态另类bdsm刘玥| 亚洲精品自拍成人| 久久综合国产亚洲精品| 丝袜脚勾引网站| 日韩欧美精品免费久久| 精品一区二区三卡| av播播在线观看一区| 校园人妻丝袜中文字幕| 在线观看一区二区三区激情| 男人舔奶头视频| 国产精品爽爽va在线观看网站| 啦啦啦啦在线视频资源| 九九在线视频观看精品| 久久ye,这里只有精品| 美女高潮的动态| 超碰97精品在线观看| 欧美一级a爱片免费观看看| 国产爱豆传媒在线观看| 国产精品精品国产色婷婷| 亚洲av综合色区一区| 身体一侧抽搐| 色视频在线一区二区三区| 久久国产精品男人的天堂亚洲 | 国产欧美亚洲国产| 色吧在线观看| 国产 精品1| 成人免费观看视频高清| xxx大片免费视频| 99热这里只有精品一区| 亚洲经典国产精华液单| 国产男人的电影天堂91| 99热全是精品| 99视频精品全部免费 在线| 色吧在线观看| 中文字幕精品免费在线观看视频 | 亚洲色图av天堂| 夫妻午夜视频| 亚洲精品一区蜜桃| av在线观看视频网站免费| 一本一本综合久久| 尾随美女入室| 在线观看三级黄色| 国产免费又黄又爽又色| 一级毛片电影观看| 亚洲最大成人中文| av国产久精品久网站免费入址| 免费av不卡在线播放| 日韩欧美精品免费久久| 国产永久视频网站| 欧美3d第一页| 精品久久久久久久久亚洲| 精品亚洲乱码少妇综合久久| 国产亚洲91精品色在线| 亚洲国产精品999| 久久6这里有精品| 亚洲经典国产精华液单| 亚洲精品日本国产第一区| 1000部很黄的大片| 国产美女午夜福利| 免费高清在线观看视频在线观看| 美女高潮的动态| 交换朋友夫妻互换小说| 成年人午夜在线观看视频| 亚洲精品乱码久久久久久按摩| 免费黄色在线免费观看| 在线观看一区二区三区| 亚洲国产毛片av蜜桃av| 一边亲一边摸免费视频| 久久精品国产亚洲av涩爱| 精品一区在线观看国产| 午夜免费鲁丝| 一本一本综合久久| 久久久久国产精品人妻一区二区| 永久网站在线| 高清欧美精品videossex| 你懂的网址亚洲精品在线观看| 中国国产av一级| 国产伦理片在线播放av一区| 国产亚洲最大av| 人妻 亚洲 视频| 丰满乱子伦码专区| 久久国产精品男人的天堂亚洲 | 舔av片在线| 大又大粗又爽又黄少妇毛片口| 成人国产麻豆网| 中文字幕制服av| 免费人妻精品一区二区三区视频| 激情 狠狠 欧美| 午夜老司机福利剧场| 国产在线男女| 亚洲一级一片aⅴ在线观看| 免费观看无遮挡的男女| 成人午夜精彩视频在线观看| 各种免费的搞黄视频| 国产精品久久久久久久久免| av在线蜜桃| 亚洲内射少妇av| 色视频在线一区二区三区| 国产精品蜜桃在线观看| 久久国内精品自在自线图片| 嫩草影院新地址| 狂野欧美白嫩少妇大欣赏| 国产成人a区在线观看| xxx大片免费视频| 久久国产亚洲av麻豆专区| 国产精品麻豆人妻色哟哟久久| 毛片女人毛片| 女人十人毛片免费观看3o分钟| av女优亚洲男人天堂| av在线老鸭窝| 国产伦在线观看视频一区| 91精品一卡2卡3卡4卡| 国产淫片久久久久久久久| 最近最新中文字幕大全电影3| 大又大粗又爽又黄少妇毛片口| 中文字幕亚洲精品专区| 91精品国产国语对白视频| 亚洲精品乱码久久久v下载方式| 国产久久久一区二区三区| 看十八女毛片水多多多| 哪个播放器可以免费观看大片| a级毛片免费高清观看在线播放| 老师上课跳d突然被开到最大视频| 久久久久久久久久成人| freevideosex欧美| 亚洲国产欧美在线一区| 午夜福利影视在线免费观看| 日韩在线高清观看一区二区三区| 亚洲国产精品专区欧美| 99热全是精品| 国产真实伦视频高清在线观看| 91精品国产国语对白视频| 国产成人精品福利久久| 亚洲国产精品一区三区| 精品亚洲乱码少妇综合久久| 各种免费的搞黄视频| 国产淫语在线视频| 伦理电影大哥的女人| 嫩草影院入口| 人人妻人人爽人人添夜夜欢视频 | 女人十人毛片免费观看3o分钟| 国产又色又爽无遮挡免| 欧美日韩精品成人综合77777| 天堂俺去俺来也www色官网| 国产伦在线观看视频一区| 自拍欧美九色日韩亚洲蝌蚪91 | 最后的刺客免费高清国语| 嘟嘟电影网在线观看| 亚洲欧美日韩无卡精品| 男男h啪啪无遮挡| 黄色日韩在线| 欧美另类一区| 人妻一区二区av| 高清日韩中文字幕在线| 热re99久久精品国产66热6| 日本黄色日本黄色录像| 国产精品av视频在线免费观看| 97超视频在线观看视频| 人妻系列 视频| 日本免费在线观看一区| 日韩免费高清中文字幕av| 久久毛片免费看一区二区三区| 大码成人一级视频| 干丝袜人妻中文字幕| 久久精品人妻少妇| 久久6这里有精品| 久久亚洲国产成人精品v| 一本—道久久a久久精品蜜桃钙片| 国产人妻一区二区三区在| 精品人妻偷拍中文字幕| 日本爱情动作片www.在线观看| av不卡在线播放| 联通29元200g的流量卡| 大片电影免费在线观看免费| 在线亚洲精品国产二区图片欧美 | 免费人成在线观看视频色| 日韩av在线免费看完整版不卡| 高清欧美精品videossex| 伊人久久精品亚洲午夜| 爱豆传媒免费全集在线观看| 色视频www国产| 人妻夜夜爽99麻豆av| 噜噜噜噜噜久久久久久91| av在线老鸭窝| 亚洲欧美成人综合另类久久久| 99热6这里只有精品| 国产男女超爽视频在线观看| 国模一区二区三区四区视频| 欧美高清成人免费视频www| 国产亚洲午夜精品一区二区久久| 国产爽快片一区二区三区| 最新中文字幕久久久久| 精品少妇久久久久久888优播| 亚洲,欧美,日韩| 又粗又硬又长又爽又黄的视频| 大又大粗又爽又黄少妇毛片口| 丰满人妻一区二区三区视频av| 中国美白少妇内射xxxbb| 爱豆传媒免费全集在线观看| 九草在线视频观看| 下体分泌物呈黄色| 国产乱人视频| 欧美丝袜亚洲另类| 久久久久久久久大av| 高清毛片免费看| 国产亚洲最大av| 国产成人一区二区在线| 日本欧美国产在线视频| 精品国产三级普通话版| 亚洲国产成人一精品久久久| 亚洲精品国产成人久久av| 色5月婷婷丁香| 免费观看无遮挡的男女| 高清视频免费观看一区二区| 亚洲高清免费不卡视频| 免费高清在线观看视频在线观看| 深爱激情五月婷婷| 国产熟女欧美一区二区| 妹子高潮喷水视频| 久久国产乱子免费精品| 久久精品国产鲁丝片午夜精品| 色哟哟·www| 国产有黄有色有爽视频| 色视频www国产| 深夜a级毛片|