• <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一区| 欧美一区二区亚洲| 欧美激情极品国产一区二区三区 | 久久久成人免费电影| 色5月婷婷丁香| 欧美另类一区| 联通29元200g的流量卡| 王馨瑶露胸无遮挡在线观看| 汤姆久久久久久久影院中文字幕| 身体一侧抽搐| 久久久精品免费免费高清| 久久久久国产精品人妻一区二区| 高清av免费在线| 日韩不卡一区二区三区视频在线| 插逼视频在线观看| 91精品国产九色| 少妇人妻久久综合中文| 久久久精品免费免费高清| 91精品国产国语对白视频| 三级国产精品欧美在线观看| 一本一本综合久久| 久久久久久久久久成人| 少妇熟女欧美另类| 国内揄拍国产精品人妻在线| 国产色爽女视频免费观看| 99久久综合免费| 国产精品免费大片| 欧美日韩视频精品一区| 99re6热这里在线精品视频| 日韩一区二区三区影片| 青青草视频在线视频观看| 国产成人精品久久久久久| 久热久热在线精品观看| 亚洲国产欧美在线一区| 午夜福利在线观看免费完整高清在| 亚洲精品一区蜜桃| 免费在线观看成人毛片| 国产在视频线精品| 国产高清国产精品国产三级 | 精品久久久久久电影网| h日本视频在线播放| 国产一级毛片在线| 亚洲国产最新在线播放| 99久久精品热视频| 成人毛片a级毛片在线播放| 狂野欧美白嫩少妇大欣赏| 婷婷色综合大香蕉| 亚洲人成网站高清观看| 国产片特级美女逼逼视频| 久久精品国产自在天天线| 免费av不卡在线播放| a级一级毛片免费在线观看| 妹子高潮喷水视频| 好男人视频免费观看在线| 26uuu在线亚洲综合色| 一个人免费看片子| 国产精品女同一区二区软件| 日本午夜av视频| av专区在线播放| 熟女人妻精品中文字幕| 在线观看国产h片| 日本wwww免费看| 99精国产麻豆久久婷婷| 美女高潮的动态| 男女边吃奶边做爰视频| 国产黄色免费在线视频| 另类亚洲欧美激情| 在线免费观看不下载黄p国产| 国产精品一区二区在线不卡| 亚洲综合色惰| 日韩av不卡免费在线播放| 性色av一级| 日韩av免费高清视频| 中国美白少妇内射xxxbb| 欧美日韩国产mv在线观看视频 | 久久 成人 亚洲| 亚洲精品aⅴ在线观看| 欧美激情极品国产一区二区三区 | 国产男人的电影天堂91| 身体一侧抽搐| 看免费成人av毛片| 一级毛片aaaaaa免费看小| 久久久久久久大尺度免费视频| 一区二区三区乱码不卡18| 啦啦啦视频在线资源免费观看| 在线观看一区二区三区| 寂寞人妻少妇视频99o| 国产精品一区二区三区四区免费观看| 欧美人与善性xxx| 日本欧美视频一区| 人妻系列 视频| 亚洲国产最新在线播放| 麻豆国产97在线/欧美| 观看av在线不卡| 成年美女黄网站色视频大全免费 | 国产精品偷伦视频观看了| 久久久久久久久久人人人人人人| 精品熟女少妇av免费看| 国产精品一区二区在线观看99| 男女无遮挡免费网站观看| 中文字幕亚洲精品专区| 我要看黄色一级片免费的| 青春草亚洲视频在线观看| 爱豆传媒免费全集在线观看| 五月伊人婷婷丁香| 91久久精品电影网| 久久久久久久久久人人人人人人| 国产伦在线观看视频一区| 成人美女网站在线观看视频| 高清在线视频一区二区三区| 成年免费大片在线观看| 国产精品一区www在线观看| 99久国产av精品国产电影| 在线观看免费高清a一片| 91久久精品国产一区二区三区| 国产亚洲91精品色在线| 五月伊人婷婷丁香| 观看av在线不卡| 又大又黄又爽视频免费| 亚洲一区二区三区欧美精品| 国产白丝娇喘喷水9色精品| 亚洲久久久国产精品| 另类亚洲欧美激情| 99久久人妻综合| 极品教师在线视频| 午夜激情福利司机影院| 春色校园在线视频观看| 在线观看美女被高潮喷水网站| 99九九线精品视频在线观看视频| 国产av国产精品国产| 麻豆国产97在线/欧美| 国内揄拍国产精品人妻在线| 亚洲精品亚洲一区二区| 日本-黄色视频高清免费观看| 国产一区二区三区综合在线观看 | 好男人视频免费观看在线| 美女xxoo啪啪120秒动态图| 亚洲一级一片aⅴ在线观看| 免费观看a级毛片全部| 国产免费福利视频在线观看| 亚洲欧美日韩东京热| 91精品一卡2卡3卡4卡| 免费少妇av软件| 国产91av在线免费观看| 国产乱来视频区| 老司机影院成人| 97热精品久久久久久| 久久久久久久精品精品| 啦啦啦中文免费视频观看日本| 欧美国产精品一级二级三级 | 精品久久久久久久久亚洲| 精品熟女少妇av免费看| 人人妻人人澡人人爽人人夜夜| 波野结衣二区三区在线| 久久久欧美国产精品| 国产毛片在线视频| 人妻制服诱惑在线中文字幕| 最近中文字幕2019免费版| 高清日韩中文字幕在线| 97在线视频观看| xxx大片免费视频| 日韩制服骚丝袜av| 99热国产这里只有精品6| 日本与韩国留学比较| 日产精品乱码卡一卡2卡三| freevideosex欧美| 亚洲av日韩在线播放| 国产精品.久久久| 成人高潮视频无遮挡免费网站| 欧美日本视频| 伊人久久国产一区二区| 亚洲图色成人| 日韩一区二区三区影片| 三级国产精品片| 少妇熟女欧美另类| 秋霞在线观看毛片| 亚洲精品第二区| 精品人妻熟女av久视频| 麻豆精品久久久久久蜜桃| 欧美精品亚洲一区二区| 狂野欧美激情性bbbbbb| 狂野欧美激情性bbbbbb| 蜜桃亚洲精品一区二区三区| 国产精品蜜桃在线观看| 老师上课跳d突然被开到最大视频| 内地一区二区视频在线| 免费黄频网站在线观看国产| 少妇精品久久久久久久| 精品一区二区三卡| 黄色一级大片看看| 国产高清不卡午夜福利| 久久精品人妻少妇| 日本午夜av视频| 菩萨蛮人人尽说江南好唐韦庄| 婷婷色综合大香蕉| 最新中文字幕久久久久| 久久久成人免费电影| 黑人高潮一二区| 夫妻午夜视频| 国产精品精品国产色婷婷| 久久97久久精品| 婷婷色综合www| 欧美精品亚洲一区二区| 嘟嘟电影网在线观看| 日日啪夜夜爽| 精品一区二区三卡| 人妻制服诱惑在线中文字幕| 国产精品一区二区性色av| 我的老师免费观看完整版| 中文欧美无线码| 少妇丰满av| 一本色道久久久久久精品综合| 国产亚洲91精品色在线| 舔av片在线| 91精品国产国语对白视频| 一本一本综合久久| 99热全是精品| 免费黄频网站在线观看国产| 另类亚洲欧美激情| 黄色配什么色好看| 高清午夜精品一区二区三区| 亚洲真实伦在线观看| 欧美bdsm另类| 日本vs欧美在线观看视频 | 91午夜精品亚洲一区二区三区| 深夜a级毛片| 精华霜和精华液先用哪个| 精品国产三级普通话版| 国产成人免费观看mmmm| 一本—道久久a久久精品蜜桃钙片| 欧美精品国产亚洲| 亚洲精品色激情综合| 激情 狠狠 欧美| 亚洲国产精品一区三区| videossex国产| 国产亚洲精品久久久com| 亚洲精品456在线播放app| 亚洲美女搞黄在线观看| 日韩电影二区| 尾随美女入室| www.av在线官网国产| 亚洲av二区三区四区| 三级经典国产精品| 大陆偷拍与自拍| 欧美最新免费一区二区三区| www.av在线官网国产| 成年美女黄网站色视频大全免费 | 18禁裸乳无遮挡免费网站照片| 亚洲国产色片| 美女视频免费永久观看网站| 在线观看一区二区三区激情| 97超碰精品成人国产| 国产成人freesex在线| 插逼视频在线观看| 欧美精品人与动牲交sv欧美| 久久精品国产自在天天线| 亚洲国产欧美在线一区| 超碰av人人做人人爽久久| 一区二区三区四区激情视频| 色网站视频免费| 国产熟女欧美一区二区| 精品人妻偷拍中文字幕| 亚洲熟女精品中文字幕| 亚洲精品乱久久久久久| 91久久精品国产一区二区成人| 国产精品一区www在线观看| 亚洲成人一二三区av| tube8黄色片| 人妻一区二区av| 成年美女黄网站色视频大全免费 | 国产精品嫩草影院av在线观看| 国产乱来视频区| 欧美 日韩 精品 国产| av专区在线播放| 女人久久www免费人成看片| 边亲边吃奶的免费视频| 最近最新中文字幕免费大全7| 偷拍熟女少妇极品色| 一级毛片我不卡| 纯流量卡能插随身wifi吗| kizo精华| 大话2 男鬼变身卡| 全区人妻精品视频| 国产男女超爽视频在线观看| 欧美精品一区二区大全| 26uuu在线亚洲综合色| 亚洲精品视频女| 一区二区av电影网| 深夜a级毛片| 亚洲av.av天堂| 成人漫画全彩无遮挡| 亚洲国产最新在线播放| 丝袜喷水一区| av福利片在线观看| 91午夜精品亚洲一区二区三区| 日韩伦理黄色片| 亚洲国产欧美人成| 丰满少妇做爰视频| 蜜桃亚洲精品一区二区三区| 人妻 亚洲 视频| 亚洲经典国产精华液单| 亚洲美女视频黄频| 伦理电影免费视频| 2021少妇久久久久久久久久久| 韩国高清视频一区二区三区| 观看免费一级毛片| 99re6热这里在线精品视频| av在线蜜桃| 草草在线视频免费看| 国产综合精华液| 国产精品熟女久久久久浪| 国产在线一区二区三区精| 国产探花极品一区二区| 日本av免费视频播放| 国产亚洲91精品色在线| 天天躁日日操中文字幕| 欧美日韩亚洲高清精品| 97超碰精品成人国产| 亚洲国产高清在线一区二区三| 3wmmmm亚洲av在线观看| 日本爱情动作片www.在线观看| 亚洲欧美日韩另类电影网站 | 久久久久久人妻| 日韩在线高清观看一区二区三区| 国产熟女欧美一区二区| 观看av在线不卡| 久久影院123| 五月天丁香电影| av黄色大香蕉| 国产在线视频一区二区| 日日啪夜夜撸| 国产精品免费大片| 国产在视频线精品| 中文字幕免费在线视频6| 日韩在线高清观看一区二区三区| 精品熟女少妇av免费看| 免费观看a级毛片全部| 99九九线精品视频在线观看视频| 国产亚洲精品久久久com| 尤物成人国产欧美一区二区三区| 不卡视频在线观看欧美| 国产极品天堂在线| 爱豆传媒免费全集在线观看| 卡戴珊不雅视频在线播放| 国产精品一二三区在线看| 日韩欧美一区视频在线观看 | 日韩成人伦理影院| 午夜福利网站1000一区二区三区| 少妇人妻 视频| 久久人人爽av亚洲精品天堂 | 人人妻人人澡人人爽人人夜夜| 国产精品福利在线免费观看| 午夜激情久久久久久久| 美女cb高潮喷水在线观看| 亚洲精品乱码久久久v下载方式| 人妻一区二区av| 日日啪夜夜撸| 男女边吃奶边做爰视频| 日本欧美国产在线视频| 高清在线视频一区二区三区| 亚洲最大成人中文| 在线观看免费日韩欧美大片 | 国产深夜福利视频在线观看| 亚洲欧美日韩东京热| 日本av免费视频播放| 搡老乐熟女国产| 国产精品99久久99久久久不卡 | 免费不卡的大黄色大毛片视频在线观看| 亚洲伊人久久精品综合| 亚洲av成人精品一二三区| 九九久久精品国产亚洲av麻豆| 七月丁香在线播放| 国产精品人妻久久久影院| av又黄又爽大尺度在线免费看| 亚洲av福利一区| 久久久精品94久久精品| 一级二级三级毛片免费看| 91在线精品国自产拍蜜月| 美女内射精品一级片tv| 国产精品熟女久久久久浪| 水蜜桃什么品种好| 久久人人爽人人片av| 中文字幕亚洲精品专区| 老师上课跳d突然被开到最大视频| 久久这里有精品视频免费| 色5月婷婷丁香| 久久国产精品大桥未久av | 欧美+日韩+精品| 久久国产乱子免费精品| 国产高清有码在线观看视频| 欧美日韩亚洲高清精品| 五月开心婷婷网| 精品久久久久久电影网| 午夜老司机福利剧场| 日日啪夜夜爽| 亚洲第一区二区三区不卡| 在线观看三级黄色| 色婷婷久久久亚洲欧美| 在线观看一区二区三区激情| 我的老师免费观看完整版| 男女边吃奶边做爰视频| 日韩一区二区三区影片| 夫妻性生交免费视频一级片| 99久久精品热视频| 97热精品久久久久久| 国产日韩欧美亚洲二区| 夜夜看夜夜爽夜夜摸| 啦啦啦在线观看免费高清www| 黄色视频在线播放观看不卡| 久热这里只有精品99| 亚洲精品,欧美精品| 91在线精品国自产拍蜜月| 内地一区二区视频在线| 成人亚洲欧美一区二区av| 我的女老师完整版在线观看| 欧美xxxx性猛交bbbb| 国产毛片在线视频| 亚洲美女视频黄频| 毛片女人毛片| 国产欧美日韩精品一区二区| 国产真实伦视频高清在线观看| 卡戴珊不雅视频在线播放| videos熟女内射| 亚洲激情五月婷婷啪啪| 91在线精品国自产拍蜜月| 国产在线一区二区三区精| 我的老师免费观看完整版| 91久久精品国产一区二区三区| 热re99久久精品国产66热6| 亚洲国产日韩一区二区| 搡女人真爽免费视频火全软件| 大片免费播放器 马上看| 新久久久久国产一级毛片| 人体艺术视频欧美日本| 国产亚洲欧美精品永久| 性色av一级| 99国产精品免费福利视频| 一区二区av电影网| 丝袜脚勾引网站| 欧美老熟妇乱子伦牲交| 日韩制服骚丝袜av| 久久精品夜色国产| 九草在线视频观看| 97超视频在线观看视频| 国产 精品1| 久久鲁丝午夜福利片| 又粗又硬又长又爽又黄的视频| 美女高潮的动态| av天堂中文字幕网| 啦啦啦在线观看免费高清www| 久久人人爽人人片av| 亚洲在久久综合| 国产欧美另类精品又又久久亚洲欧美| a 毛片基地| 欧美老熟妇乱子伦牲交| 精品人妻视频免费看| 丰满少妇做爰视频| 国产欧美日韩精品一区二区| 日本与韩国留学比较| 美女福利国产在线 | 少妇的逼水好多| 在线观看一区二区三区| 国产精品熟女久久久久浪| 青春草视频在线免费观看| 亚洲精品久久久久久婷婷小说| 嫩草影院入口| 99久久精品一区二区三区| 草草在线视频免费看| av专区在线播放| 久久精品熟女亚洲av麻豆精品| 午夜日本视频在线| 午夜福利视频精品| 亚洲无线观看免费| 五月玫瑰六月丁香| 欧美3d第一页| 天堂中文最新版在线下载| 久久综合国产亚洲精品| 精品亚洲成国产av| 久久国产乱子免费精品| 国内精品宾馆在线| 成人高潮视频无遮挡免费网站| 夫妻午夜视频| 国产乱人偷精品视频| 国产亚洲5aaaaa淫片| 成人特级av手机在线观看| 久久人人爽人人片av| 成年美女黄网站色视频大全免费 | 久久久欧美国产精品| 日韩伦理黄色片| 51国产日韩欧美| 欧美最新免费一区二区三区| 青青草视频在线视频观看| 少妇精品久久久久久久| av免费观看日本| 国产精品一区www在线观看| 久久久精品免费免费高清| 国产精品不卡视频一区二区| 午夜福利高清视频| 欧美成人a在线观看| 十八禁网站网址无遮挡 | 国产一区有黄有色的免费视频| 在线观看av片永久免费下载| 一级二级三级毛片免费看| 亚洲人成网站在线播| 高清午夜精品一区二区三区| 国国产精品蜜臀av免费| 嘟嘟电影网在线观看| 国产在线免费精品| 国产国拍精品亚洲av在线观看| 热re99久久精品国产66热6| 色网站视频免费| 免费av中文字幕在线| 亚洲精品成人av观看孕妇| 欧美日韩综合久久久久久| 成人漫画全彩无遮挡| 99热这里只有精品一区| freevideosex欧美| 好男人视频免费观看在线| 美女中出高潮动态图| av在线观看视频网站免费| 一区二区三区四区激情视频| 尤物成人国产欧美一区二区三区| 中文字幕av成人在线电影| 啦啦啦在线观看免费高清www| 国产精品久久久久久久久免| 黑丝袜美女国产一区| 国产精品一区二区在线不卡| 精品亚洲成国产av| 欧美日韩精品成人综合77777| 国产免费一区二区三区四区乱码| 亚洲三级黄色毛片| 亚洲美女黄色视频免费看| 99久国产av精品国产电影| 在线播放无遮挡| 五月伊人婷婷丁香| 男人和女人高潮做爰伦理| 亚洲av中文字字幕乱码综合| 亚洲精品日韩av片在线观看| 国产精品嫩草影院av在线观看| 亚洲无线观看免费| 日产精品乱码卡一卡2卡三| 国内精品宾馆在线| 成人高潮视频无遮挡免费网站| 国产精品爽爽va在线观看网站| 国产大屁股一区二区在线视频| 亚洲电影在线观看av| 久久青草综合色| 久久精品熟女亚洲av麻豆精品| 婷婷色综合大香蕉| 国产精品熟女久久久久浪| 亚洲三级黄色毛片| 九九久久精品国产亚洲av麻豆| 久久久久国产精品人妻一区二区| 精品国产一区二区三区久久久樱花 | 亚洲精品第二区| 日日摸夜夜添夜夜添av毛片| 亚洲美女视频黄频| 国产中年淑女户外野战色| 日本-黄色视频高清免费观看| 91精品国产国语对白视频| a级一级毛片免费在线观看| 草草在线视频免费看| 亚洲成人手机| 内射极品少妇av片p| 国产成人精品婷婷| 国产v大片淫在线免费观看| 亚洲av二区三区四区| 激情五月婷婷亚洲| 久久久久久久亚洲中文字幕| 26uuu在线亚洲综合色| a 毛片基地| 国产免费福利视频在线观看| 日本av手机在线免费观看| 成人美女网站在线观看视频| 国产伦精品一区二区三区视频9| 特大巨黑吊av在线直播| 91久久精品电影网| 欧美xxxx黑人xx丫x性爽| 亚洲久久久国产精品| 亚州av有码| 九九爱精品视频在线观看| 亚洲精品亚洲一区二区| 大香蕉97超碰在线| 日韩国内少妇激情av| 亚洲欧洲日产国产| 国产伦精品一区二区三区四那| 在线观看三级黄色| 久久国产精品男人的天堂亚洲 | 国产永久视频网站| 精品人妻熟女av久视频| 午夜福利高清视频| 乱系列少妇在线播放| 久久鲁丝午夜福利片| 色哟哟·www| av视频免费观看在线观看| 97在线视频观看| a级一级毛片免费在线观看| av免费在线看不卡| 亚洲国产毛片av蜜桃av| 精品午夜福利在线看| 亚洲国产毛片av蜜桃av| 一本—道久久a久久精品蜜桃钙片| 涩涩av久久男人的天堂| 国产亚洲5aaaaa淫片| 高清欧美精品videossex| 国产熟女欧美一区二区| 国产综合精华液| 久久久久国产网址| 免费久久久久久久精品成人欧美视频 | av不卡在线播放| 91狼人影院| 成人特级av手机在线观看|