• <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ū)沉積物古菌多樣性分析
    18禁黄网站禁片午夜丰满| 欧美极品一区二区三区四区| videossex国产| 亚洲精品在线观看二区| 亚洲av中文av极速乱 | 超碰av人人做人人爽久久| 嫩草影院新地址| 亚洲精品色激情综合| 好男人在线观看高清免费视频| 日本 欧美在线| 久久天躁狠狠躁夜夜2o2o| 观看美女的网站| 两个人视频免费观看高清| 午夜福利视频1000在线观看| 亚洲最大成人手机在线| 欧美成人一区二区免费高清观看| 国内精品宾馆在线| 能在线免费观看的黄片| 非洲黑人性xxxx精品又粗又长| 97热精品久久久久久| 69av精品久久久久久| 亚洲国产精品成人综合色| 色精品久久人妻99蜜桃| 嫁个100分男人电影在线观看| 毛片女人毛片| 欧美色视频一区免费| 国产真实乱freesex| 狂野欧美白嫩少妇大欣赏| 最新在线观看一区二区三区| 一卡2卡三卡四卡精品乱码亚洲| 91麻豆av在线| 国产亚洲精品综合一区在线观看| av天堂在线播放| 99在线人妻在线中文字幕| 亚洲精品久久国产高清桃花| 国产大屁股一区二区在线视频| 久久久久国产精品人妻aⅴ院| 免费看a级黄色片| 变态另类丝袜制服| 在线免费十八禁| 国产综合懂色| 免费观看在线日韩| 日本精品一区二区三区蜜桃| 精品久久久久久成人av| 欧美激情久久久久久爽电影| 国产 一区精品| 九九热线精品视视频播放| 我要看日韩黄色一级片| 国产精品美女特级片免费视频播放器| 亚洲图色成人| 免费av观看视频| 欧美一区二区精品小视频在线| 日韩欧美三级三区| 成年女人永久免费观看视频| 国产精品人妻久久久影院| 色吧在线观看| 国产国拍精品亚洲av在线观看| 丰满乱子伦码专区| 99国产精品一区二区蜜桃av| 亚洲人成网站在线播| 不卡视频在线观看欧美| 国产男靠女视频免费网站| 亚洲熟妇中文字幕五十中出| 精品一区二区三区视频在线| 色av中文字幕| a级一级毛片免费在线观看| 国产成人一区二区在线| 中文字幕精品亚洲无线码一区| 在线观看66精品国产| 国产色婷婷99| 一卡2卡三卡四卡精品乱码亚洲| 18禁黄网站禁片午夜丰满| 亚洲久久久久久中文字幕| 国产色婷婷99| 亚洲色图av天堂| 全区人妻精品视频| 色在线成人网| 久久久久久伊人网av| 国产国拍精品亚洲av在线观看| 97超视频在线观看视频| 国内毛片毛片毛片毛片毛片| 欧美成人a在线观看| 少妇裸体淫交视频免费看高清| 日本 欧美在线| 中文字幕精品亚洲无线码一区| 老司机福利观看| 又黄又爽又刺激的免费视频.| 日本一本二区三区精品| 亚洲国产高清在线一区二区三| 亚洲一级一片aⅴ在线观看| 毛片女人毛片| 麻豆一二三区av精品| 亚洲专区中文字幕在线| 亚洲av第一区精品v没综合| 人妻久久中文字幕网| 亚洲av五月六月丁香网| 五月玫瑰六月丁香| xxxwww97欧美| av福利片在线观看| 色哟哟·www| 精品无人区乱码1区二区| netflix在线观看网站| 午夜福利18| 91在线精品国自产拍蜜月| 精品久久久久久成人av| 成人av一区二区三区在线看| 一个人免费在线观看电影| 欧美成人免费av一区二区三区| 国语自产精品视频在线第100页| 国产精品久久久久久精品电影| 99九九线精品视频在线观看视频| 亚洲欧美日韩高清专用| 一级av片app| 亚洲在线自拍视频| 久久久久久久久久成人| 亚洲电影在线观看av| 国产精品无大码| 精品久久国产蜜桃| 18禁裸乳无遮挡免费网站照片| 干丝袜人妻中文字幕| 成人综合一区亚洲| .国产精品久久| 亚洲美女黄片视频| 成年女人永久免费观看视频| 久久草成人影院| 成年版毛片免费区| 亚洲精品粉嫩美女一区| 蜜桃亚洲精品一区二区三区| 亚洲第一区二区三区不卡| 亚洲av成人av| 非洲黑人性xxxx精品又粗又长| 欧美绝顶高潮抽搐喷水| 国产精品免费一区二区三区在线| 成人性生交大片免费视频hd| 欧美日韩国产亚洲二区| 美女大奶头视频| 69人妻影院| 麻豆成人午夜福利视频| 国产亚洲精品av在线| 国产不卡一卡二| 欧美一区二区亚洲| 久久精品国产99精品国产亚洲性色| 色综合婷婷激情| 97人妻精品一区二区三区麻豆| or卡值多少钱| 亚洲av不卡在线观看| 两性午夜刺激爽爽歪歪视频在线观看| eeuss影院久久| 性色avwww在线观看| 国产av一区在线观看免费| 亚洲经典国产精华液单| 久久精品影院6| 国产亚洲精品久久久com| 在线播放国产精品三级| 欧美三级亚洲精品| 久久久久久久久久黄片| 欧美丝袜亚洲另类 | 大型黄色视频在线免费观看| 国产成人福利小说| 午夜爱爱视频在线播放| 不卡视频在线观看欧美| 国产不卡一卡二| 日本爱情动作片www.在线观看 | 国产精品一区二区三区四区久久| 国产精品久久久久久av不卡| 欧洲精品卡2卡3卡4卡5卡区| 欧美+亚洲+日韩+国产| 国产一区二区在线av高清观看| 国产精品嫩草影院av在线观看 | www日本黄色视频网| 欧美一区二区亚洲| 免费黄网站久久成人精品| 黄色配什么色好看| 午夜福利欧美成人| 男女之事视频高清在线观看| 亚洲不卡免费看| 一进一出抽搐gif免费好疼| 国产精品人妻久久久影院| 午夜免费激情av| 乱码一卡2卡4卡精品| 在线天堂最新版资源| 色综合亚洲欧美另类图片| 麻豆国产97在线/欧美| 精品欧美国产一区二区三| 十八禁国产超污无遮挡网站| 国产高潮美女av| 精品久久久久久久久av| 国产视频内射| 国产免费一级a男人的天堂| 国产亚洲91精品色在线| 国内精品久久久久久久电影| 欧美绝顶高潮抽搐喷水| 88av欧美| 国产中年淑女户外野战色| 亚洲欧美日韩东京热| 亚洲成人中文字幕在线播放| 欧美激情在线99| 99热这里只有精品一区| 国产极品精品免费视频能看的| 欧美精品国产亚洲| 别揉我奶头~嗯~啊~动态视频| 亚洲最大成人中文| 久久中文看片网| 成年人黄色毛片网站| 欧美一级a爱片免费观看看| 三级国产精品欧美在线观看| 国产一区二区激情短视频| 国内精品宾馆在线| 亚洲欧美日韩东京热| 日本 av在线| 成人av在线播放网站| 国产毛片a区久久久久| 亚洲国产欧洲综合997久久,| 在线看三级毛片| 在线观看66精品国产| 亚洲av免费在线观看| 亚洲一区高清亚洲精品| 伦理电影大哥的女人| 久久久久久久午夜电影| 又爽又黄无遮挡网站| 国产一级毛片七仙女欲春2| АⅤ资源中文在线天堂| 久久99热6这里只有精品| 亚洲成人免费电影在线观看| 亚洲国产欧美人成| 成年女人毛片免费观看观看9| 成人国产综合亚洲| 亚洲无线在线观看| 校园春色视频在线观看| 亚洲电影在线观看av| 一个人免费在线观看电影| 蜜桃久久精品国产亚洲av| 夜夜看夜夜爽夜夜摸| 国产精品女同一区二区软件 | 此物有八面人人有两片| 日本欧美国产在线视频| 久久久久免费精品人妻一区二区| 日韩欧美一区二区三区在线观看| 在线免费观看不下载黄p国产 | 欧美黑人巨大hd| 毛片女人毛片| 亚洲久久久久久中文字幕| 赤兔流量卡办理| 人人妻人人澡欧美一区二区| 久久久久国内视频| 亚洲最大成人中文| 91av网一区二区| 久久久久久久午夜电影| 日本a在线网址| 日韩亚洲欧美综合| 99热这里只有是精品50| 成人永久免费在线观看视频| 999久久久精品免费观看国产| 日本五十路高清| 国产亚洲91精品色在线| av视频在线观看入口| 亚洲av电影不卡..在线观看| 亚洲精华国产精华液的使用体验 | 亚洲无线观看免费| 久久国产乱子免费精品| 国产一区二区激情短视频| 亚洲欧美日韩高清专用| 欧美最黄视频在线播放免费| 日本与韩国留学比较| 丰满人妻一区二区三区视频av| 日韩亚洲欧美综合| 国产视频一区二区在线看| 久久精品影院6| 高清在线国产一区| 久久久午夜欧美精品| 成年人黄色毛片网站| 免费av观看视频| 一区福利在线观看| 麻豆av噜噜一区二区三区| 国语自产精品视频在线第100页| 天堂动漫精品| 色视频www国产| 日本黄大片高清| 日韩强制内射视频| 成人美女网站在线观看视频| 久久精品国产自在天天线| 禁无遮挡网站| 夜夜看夜夜爽夜夜摸| 俄罗斯特黄特色一大片| 亚洲图色成人| 国产真实伦视频高清在线观看 | 中文字幕av在线有码专区| 丰满的人妻完整版| 亚洲欧美激情综合另类| 亚洲经典国产精华液单| 亚洲熟妇中文字幕五十中出| 国产精品嫩草影院av在线观看 | 特级一级黄色大片| 大型黄色视频在线免费观看| 麻豆国产97在线/欧美| xxxwww97欧美| 午夜日韩欧美国产| 亚洲 国产 在线| 国产av不卡久久| 国产人妻一区二区三区在| 欧美不卡视频在线免费观看| 搡老岳熟女国产| 高清毛片免费观看视频网站| 两性午夜刺激爽爽歪歪视频在线观看| 天堂动漫精品| 老女人水多毛片| 国产三级在线视频| 麻豆av噜噜一区二区三区| 女同久久另类99精品国产91| 成人午夜高清在线视频| 国产精品野战在线观看| 波多野结衣高清作品| 国产精品久久久久久久久免| 99久久久亚洲精品蜜臀av| 成人无遮挡网站| 成人精品一区二区免费| 精品乱码久久久久久99久播| av黄色大香蕉| 美女大奶头视频| 亚洲成a人片在线一区二区| 国产黄片美女视频| 听说在线观看完整版免费高清| 精品国内亚洲2022精品成人| 草草在线视频免费看| 色av中文字幕| 日日啪夜夜撸| 97超级碰碰碰精品色视频在线观看| 国产高清不卡午夜福利| or卡值多少钱| av在线亚洲专区| 国产三级在线视频| 久久亚洲真实| 中文在线观看免费www的网站| 国产精品久久久久久久电影| 久久6这里有精品| 国产高清三级在线| 亚洲自偷自拍三级| 校园春色视频在线观看| 国产不卡一卡二| 国产成人aa在线观看| 亚洲性久久影院| 波多野结衣高清作品| 亚洲av不卡在线观看| 精华霜和精华液先用哪个| 成人av在线播放网站| 欧美激情国产日韩精品一区| 嫩草影院新地址| 欧美一级a爱片免费观看看| 少妇人妻一区二区三区视频| 51国产日韩欧美| 精品人妻偷拍中文字幕| 久久国内精品自在自线图片| 黄色一级大片看看| АⅤ资源中文在线天堂| 真人做人爱边吃奶动态| 丰满人妻一区二区三区视频av| av黄色大香蕉| 亚洲精品亚洲一区二区| 51国产日韩欧美| 欧美+亚洲+日韩+国产| 日韩高清综合在线| 午夜福利在线观看免费完整高清在 | 日韩欧美精品免费久久| 看免费成人av毛片| 性插视频无遮挡在线免费观看| 中文字幕高清在线视频| 国产精品精品国产色婷婷| 国产高清三级在线| 国产成年人精品一区二区| 午夜福利成人在线免费观看| 国产高潮美女av| 久99久视频精品免费| 又紧又爽又黄一区二区| 亚洲va在线va天堂va国产| 熟妇人妻久久中文字幕3abv| 欧美在线一区亚洲| 国产亚洲精品综合一区在线观看| АⅤ资源中文在线天堂| 窝窝影院91人妻| 最新中文字幕久久久久| 午夜激情欧美在线| 99久久久亚洲精品蜜臀av| 国产精品久久久久久av不卡| 特级一级黄色大片| 久久亚洲真实| 搞女人的毛片| 狂野欧美激情性xxxx在线观看| 亚洲图色成人| 男女做爰动态图高潮gif福利片| 久久精品国产亚洲网站| 国产高潮美女av| 一进一出抽搐动态| 日本黄大片高清| 成人性生交大片免费视频hd| 亚洲第一区二区三区不卡| 日韩亚洲欧美综合| 一区二区三区激情视频| 好男人在线观看高清免费视频| 看十八女毛片水多多多| 99riav亚洲国产免费| 黄色一级大片看看| 超碰av人人做人人爽久久| 日韩国内少妇激情av| 国产精品爽爽va在线观看网站| 国产精品人妻久久久久久| 亚洲18禁久久av| 免费人成在线观看视频色| 香蕉av资源在线| 国产成人a区在线观看| 国产一区二区在线观看日韩| 直男gayav资源| 亚洲成人中文字幕在线播放| 麻豆av噜噜一区二区三区| 亚洲美女黄片视频| 午夜激情欧美在线| 亚洲成人中文字幕在线播放| 高清毛片免费观看视频网站| 欧美精品啪啪一区二区三区| 久久99热这里只有精品18| 亚洲五月天丁香| 人妻久久中文字幕网| 国产午夜福利久久久久久| 精品久久久久久久久av| 国产中年淑女户外野战色| 91狼人影院| 国产高清激情床上av| 春色校园在线视频观看| 精品一区二区三区av网在线观看| 久久热精品热| 亚洲性夜色夜夜综合| 国产v大片淫在线免费观看| 国语自产精品视频在线第100页| 免费不卡的大黄色大毛片视频在线观看 | 又爽又黄无遮挡网站| 久久欧美精品欧美久久欧美| 午夜精品久久久久久毛片777| 精品不卡国产一区二区三区| 国产精品一区www在线观看 | 一本久久中文字幕| 不卡视频在线观看欧美| 色精品久久人妻99蜜桃| 少妇高潮的动态图| 精品免费久久久久久久清纯| 午夜日韩欧美国产| 麻豆国产97在线/欧美| 最近最新中文字幕大全电影3| 一a级毛片在线观看| 春色校园在线视频观看| 国产精品综合久久久久久久免费| 亚洲四区av| 日本黄色视频三级网站网址| 成人性生交大片免费视频hd| 国产精品亚洲一级av第二区| 岛国在线免费视频观看| 欧美+日韩+精品| 国产精品av视频在线免费观看| 国产综合懂色| 色综合亚洲欧美另类图片| 亚洲avbb在线观看| 亚洲成人精品中文字幕电影| 亚洲avbb在线观看| 性插视频无遮挡在线免费观看| 一本精品99久久精品77| 精品久久久久久,| 成人av在线播放网站| 少妇高潮的动态图| 日本精品一区二区三区蜜桃| 中文字幕免费在线视频6| 又粗又爽又猛毛片免费看| 亚洲精品影视一区二区三区av| 日日摸夜夜添夜夜添av毛片 | 久久久久久久久中文| www.色视频.com| 精品99又大又爽又粗少妇毛片 | 国产aⅴ精品一区二区三区波| 一a级毛片在线观看| 国产精品久久久久久久久免| 97超视频在线观看视频| 国产av一区在线观看免费| 日本色播在线视频| 成人二区视频| 国产精品国产三级国产av玫瑰| 高清在线国产一区| 黄色女人牲交| 1024手机看黄色片| 国产精品无大码| 日本黄色片子视频| 不卡一级毛片| 97超级碰碰碰精品色视频在线观看| 悠悠久久av| 99久久精品国产国产毛片| 欧美最新免费一区二区三区| 色播亚洲综合网| 搞女人的毛片| 嫩草影院精品99| h日本视频在线播放| 女人十人毛片免费观看3o分钟| 国产伦在线观看视频一区| 国产精品日韩av在线免费观看| 国产精品亚洲美女久久久| 久久久精品欧美日韩精品| 日韩一区二区视频免费看| 国产一区二区激情短视频| 黄色配什么色好看| 夜夜爽天天搞| 欧美最新免费一区二区三区| 亚洲成人久久爱视频| 国产爱豆传媒在线观看| 色av中文字幕| 久久久成人免费电影| av在线老鸭窝| 国产欧美日韩精品亚洲av| 欧美丝袜亚洲另类 | 色av中文字幕| 欧美bdsm另类| 中文字幕久久专区| 天堂√8在线中文| 久久久久精品国产欧美久久久| 精品久久久久久,| 日本一本二区三区精品| 久久人人爽人人爽人人片va| 中文字幕久久专区| 亚洲欧美日韩卡通动漫| 少妇高潮的动态图| 啪啪无遮挡十八禁网站| 国产精品综合久久久久久久免费| 国产中年淑女户外野战色| 欧美日本亚洲视频在线播放| 国产高清三级在线| 精品久久久久久久久av| 色综合婷婷激情| 97热精品久久久久久| 国产精品久久久久久久电影| 日韩亚洲欧美综合| 国产一区二区在线观看日韩| 黄色丝袜av网址大全| 亚洲av美国av| 国产免费男女视频| 久久6这里有精品| 亚洲色图av天堂| 国产麻豆成人av免费视频| 99国产极品粉嫩在线观看| 成人鲁丝片一二三区免费| 18禁黄网站禁片免费观看直播| 亚洲在线自拍视频| 男女之事视频高清在线观看| 国产精品一区二区性色av| 欧美成人免费av一区二区三区| 亚洲综合色惰| 久久国产精品人妻蜜桃| 中文字幕av在线有码专区| 老熟妇仑乱视频hdxx| 少妇猛男粗大的猛烈进出视频 | 一本精品99久久精品77| 淫秽高清视频在线观看| 久久精品国产亚洲av涩爱 | 亚洲国产精品sss在线观看| av在线老鸭窝| 国产精品人妻久久久影院| 日本熟妇午夜| 成人无遮挡网站| 亚洲乱码一区二区免费版| 免费看日本二区| 久久精品久久久久久噜噜老黄 | 婷婷六月久久综合丁香| 精品人妻一区二区三区麻豆 | 哪里可以看免费的av片| 黄色一级大片看看| 性色avwww在线观看| 亚洲欧美日韩无卡精品| 日本在线视频免费播放| 欧美性感艳星| 亚洲色图av天堂| 日韩亚洲欧美综合| 人人妻,人人澡人人爽秒播| 欧美另类亚洲清纯唯美| 天堂√8在线中文| 亚洲经典国产精华液单| 午夜福利视频1000在线观看| 3wmmmm亚洲av在线观看| 中文字幕人妻熟人妻熟丝袜美| 国产成年人精品一区二区| 国产av一区在线观看免费| 免费看a级黄色片| 99精品久久久久人妻精品| 亚洲人成网站在线播放欧美日韩| 听说在线观看完整版免费高清| 内射极品少妇av片p| av在线天堂中文字幕| 中文字幕高清在线视频| 超碰av人人做人人爽久久| 在线国产一区二区在线| 九九热线精品视视频播放| 久久久久久久久久成人| 日韩欧美免费精品| 白带黄色成豆腐渣| 男女之事视频高清在线观看| 午夜免费男女啪啪视频观看 | 日本黄大片高清| 在线免费十八禁| 日韩中文字幕欧美一区二区| 亚洲精品影视一区二区三区av| 成年女人看的毛片在线观看| 最后的刺客免费高清国语| 伦精品一区二区三区| av视频在线观看入口| 日韩精品中文字幕看吧| 97热精品久久久久久| 国产乱人伦免费视频| 亚洲av熟女| 麻豆成人av在线观看| 亚洲五月天丁香| 国产伦在线观看视频一区| 亚洲最大成人av| 国产精品人妻久久久久久|