• <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| 亚洲欧美精品专区久久| 3wmmmm亚洲av在线观看| 嫩草影院精品99| 卡戴珊不雅视频在线播放| 秋霞在线观看毛片| 国产色婷婷99| av专区在线播放| 九草在线视频观看| 久久国产乱子免费精品| 国产精品爽爽va在线观看网站| 国产一区二区三区av在线 | 日本欧美国产在线视频| а√天堂www在线а√下载| 色噜噜av男人的天堂激情| 亚州av有码| 免费大片18禁| 日本免费a在线| 校园人妻丝袜中文字幕| 又粗又硬又长又爽又黄的视频 | 国产精品精品国产色婷婷| 精品久久国产蜜桃| 在线免费观看的www视频| 亚洲欧洲国产日韩| 国产黄片美女视频| 国产亚洲精品久久久久久毛片| 亚洲性久久影院| ponron亚洲| 久久精品夜夜夜夜夜久久蜜豆| 国产一区二区三区在线臀色熟女| 51国产日韩欧美| 精品国内亚洲2022精品成人| 身体一侧抽搐| 99在线人妻在线中文字幕| 最近中文字幕高清免费大全6| 一个人免费在线观看电影| 能在线免费观看的黄片| 天堂影院成人在线观看| .国产精品久久| 久久国产乱子免费精品| 亚洲av熟女| 久久久国产成人免费| 久久久精品大字幕| 日本-黄色视频高清免费观看| 国内揄拍国产精品人妻在线| 丝袜美腿在线中文| 在线播放国产精品三级| 午夜视频国产福利| 一级黄色大片毛片| 婷婷六月久久综合丁香| 97超视频在线观看视频| 99久国产av精品国产电影| 国产又黄又爽又无遮挡在线| av在线蜜桃| 亚洲av成人av| 能在线免费观看的黄片| 亚洲性久久影院| 亚洲欧美精品专区久久| 男女视频在线观看网站免费| 国产乱人偷精品视频| 小说图片视频综合网站| 九九久久精品国产亚洲av麻豆| 观看美女的网站| 日本撒尿小便嘘嘘汇集6| 人人妻人人看人人澡| 精品人妻偷拍中文字幕| 欧美xxxx黑人xx丫x性爽| 国产精品久久电影中文字幕| 不卡视频在线观看欧美| 婷婷色综合大香蕉| 免费大片18禁| 日本三级黄在线观看| 少妇的逼水好多| 亚洲欧美成人精品一区二区| 五月玫瑰六月丁香| 国产成人福利小说| 成人漫画全彩无遮挡| 亚洲av第一区精品v没综合| 淫秽高清视频在线观看| 黑人高潮一二区| 国产精品人妻久久久影院| 听说在线观看完整版免费高清| 美女cb高潮喷水在线观看| 国产又黄又爽又无遮挡在线| 亚洲精品日韩在线中文字幕 | av在线老鸭窝| 亚洲精品日韩在线中文字幕 | 九草在线视频观看| 亚洲欧洲国产日韩| 高清毛片免费观看视频网站| 国产av一区在线观看免费| 91狼人影院| 国产午夜精品久久久久久一区二区三区| 国国产精品蜜臀av免费| 啦啦啦啦在线视频资源| 亚洲最大成人中文| avwww免费| 欧美成人a在线观看| 老师上课跳d突然被开到最大视频| 三级经典国产精品| 91麻豆精品激情在线观看国产| 18禁裸乳无遮挡免费网站照片| 国产久久久一区二区三区| 国产精品1区2区在线观看.| 热99在线观看视频| 亚洲av免费高清在线观看| av在线老鸭窝| 91午夜精品亚洲一区二区三区| 99热网站在线观看| 高清毛片免费观看视频网站| 国产探花在线观看一区二区| 村上凉子中文字幕在线| 欧美zozozo另类| 久久精品国产鲁丝片午夜精品| 看十八女毛片水多多多| 亚洲国产欧美在线一区| 久久精品久久久久久久性| 精品少妇黑人巨大在线播放 | 99热只有精品国产| 精品免费久久久久久久清纯| 成人性生交大片免费视频hd| 十八禁国产超污无遮挡网站| 日日啪夜夜撸| 一区二区三区四区激情视频 | 草草在线视频免费看| 亚洲va在线va天堂va国产| 国产成人a区在线观看| 一进一出抽搐gif免费好疼| 有码 亚洲区| 国产中年淑女户外野战色| 尾随美女入室| 日韩中字成人| 午夜福利视频1000在线观看| 波多野结衣巨乳人妻| 国内少妇人妻偷人精品xxx网站| 精品国产三级普通话版| 天美传媒精品一区二区| 非洲黑人性xxxx精品又粗又长| 国产久久久一区二区三区| 麻豆成人av视频| 最好的美女福利视频网| 我要搜黄色片| 亚洲不卡免费看| 久久精品91蜜桃| 久久中文看片网| 精品人妻熟女av久视频| 成年版毛片免费区| 五月玫瑰六月丁香| 欧美成人免费av一区二区三区| 欧美不卡视频在线免费观看| 国产一区二区亚洲精品在线观看| 老司机福利观看| 欧美三级亚洲精品| 在线观看美女被高潮喷水网站| 99热6这里只有精品| 不卡视频在线观看欧美| АⅤ资源中文在线天堂| 美女高潮的动态| 国产黄片美女视频| 色哟哟·www| 美女高潮的动态| 在线观看66精品国产| 亚州av有码| 51国产日韩欧美| 国产精品乱码一区二三区的特点| 美女国产视频在线观看| 欧美色欧美亚洲另类二区| 大香蕉久久网| 成年女人永久免费观看视频| 亚洲美女搞黄在线观看| 国产精品综合久久久久久久免费| 日韩强制内射视频| 一本精品99久久精品77| 亚洲欧美日韩无卡精品| 美女黄网站色视频| 免费搜索国产男女视频| 日韩欧美精品免费久久| 亚洲图色成人| 亚洲av二区三区四区| 亚洲精品自拍成人| 亚洲av成人av| 亚洲精品久久国产高清桃花| 成人亚洲精品av一区二区| 中文资源天堂在线| 久久久午夜欧美精品| 国产一区二区三区在线臀色熟女| 一卡2卡三卡四卡精品乱码亚洲| 国产精品一区二区在线观看99 | 人妻制服诱惑在线中文字幕| 精品久久国产蜜桃| 中文字幕人妻熟人妻熟丝袜美| 免费观看在线日韩| 丝袜美腿在线中文| 少妇的逼水好多| av免费观看日本| 亚洲av一区综合| 又爽又黄无遮挡网站| 一个人免费在线观看电影| 精品久久久久久成人av| 国产精品野战在线观看| 欧美三级亚洲精品| av在线播放精品| 十八禁国产超污无遮挡网站| 亚洲乱码一区二区免费版| 亚洲精品影视一区二区三区av| 一本精品99久久精品77| 少妇被粗大猛烈的视频| 少妇熟女欧美另类| 国产成人精品久久久久久| 欧美性感艳星| 免费观看精品视频网站| 国产精品福利在线免费观看| 久久久国产成人免费| 欧美激情久久久久久爽电影| 九九热线精品视视频播放| 日本在线视频免费播放| 色哟哟哟哟哟哟| 久久综合国产亚洲精品| 日韩大尺度精品在线看网址| 色噜噜av男人的天堂激情| 国产老妇伦熟女老妇高清| 日韩成人伦理影院| 国内揄拍国产精品人妻在线| 欧美xxxx性猛交bbbb| 国产精品.久久久| 亚洲国产欧美人成| 乱码一卡2卡4卡精品| 亚洲成人精品中文字幕电影| 97热精品久久久久久| 亚洲欧美成人精品一区二区| 长腿黑丝高跟| 中文字幕制服av| 人妻夜夜爽99麻豆av| 久久6这里有精品| 成人欧美大片| 日韩制服骚丝袜av| 桃色一区二区三区在线观看| 精品不卡国产一区二区三区| 天美传媒精品一区二区| 91精品一卡2卡3卡4卡| 亚洲av免费在线观看| 国产精品蜜桃在线观看 | 久久综合国产亚洲精品| 赤兔流量卡办理| 菩萨蛮人人尽说江南好唐韦庄 | 欧美另类亚洲清纯唯美| 国产乱人视频| 欧美色视频一区免费| 欧美人与善性xxx| 亚洲国产欧美在线一区| 高清毛片免费看| 亚洲精品亚洲一区二区| 国产精品久久久久久精品电影小说 | 在线观看午夜福利视频| 国内久久婷婷六月综合欲色啪| 欧美激情久久久久久爽电影| 深夜a级毛片| 男插女下体视频免费在线播放| 亚洲国产日韩欧美精品在线观看| 亚洲国产欧洲综合997久久,| 桃色一区二区三区在线观看| 亚洲av中文字字幕乱码综合| 寂寞人妻少妇视频99o| 色综合站精品国产| 久久午夜亚洲精品久久| 亚洲天堂国产精品一区在线| 日韩人妻高清精品专区| 毛片女人毛片| 精品不卡国产一区二区三区| 国产伦精品一区二区三区视频9| 免费搜索国产男女视频| 婷婷色av中文字幕| 伦理电影大哥的女人| 91麻豆精品激情在线观看国产| 久久九九热精品免费| 免费av观看视频| 五月伊人婷婷丁香| 亚洲欧美日韩高清专用| 小说图片视频综合网站| 最近中文字幕高清免费大全6| 日韩大尺度精品在线看网址| 精品一区二区免费观看| 欧美高清性xxxxhd video| 欧美日本亚洲视频在线播放| 成人av在线播放网站| 日本成人三级电影网站| 一区二区三区免费毛片| 国产白丝娇喘喷水9色精品| 美女xxoo啪啪120秒动态图| 国产色婷婷99| 中文欧美无线码| av卡一久久| 尤物成人国产欧美一区二区三区| 此物有八面人人有两片| 99在线视频只有这里精品首页| 国产伦精品一区二区三区四那| 国产av在哪里看| 日本免费a在线| 国产色爽女视频免费观看| www日本黄色视频网| 麻豆精品久久久久久蜜桃| 在线观看免费视频日本深夜| 亚洲自偷自拍三级| 国产综合懂色| 美女国产视频在线观看| 久久草成人影院| 寂寞人妻少妇视频99o| 久久午夜福利片| 日本熟妇午夜| 日本免费a在线| 亚洲三级黄色毛片| 国产成人精品婷婷| 亚洲天堂国产精品一区在线| 国产精品av视频在线免费观看| 在线免费观看不下载黄p国产| 你懂的网址亚洲精品在线观看 | 色视频www国产| 亚洲激情五月婷婷啪啪| 99在线视频只有这里精品首页| 亚洲经典国产精华液单| 欧美激情在线99| 精品少妇黑人巨大在线播放 | 精品人妻熟女av久视频| 久久久久久九九精品二区国产| 99国产精品一区二区蜜桃av| av在线亚洲专区| 国产高清视频在线观看网站| 国产一级毛片在线| 免费av不卡在线播放| 舔av片在线| 99久久精品一区二区三区| 国产一区亚洲一区在线观看| 老女人水多毛片| www.色视频.com| 一卡2卡三卡四卡精品乱码亚洲| 日韩av不卡免费在线播放| 男人的好看免费观看在线视频| a级一级毛片免费在线观看| 欧美日韩在线观看h| 日韩高清综合在线| 精品久久久噜噜| 国产精品一及| 日本一本二区三区精品| 日韩欧美三级三区| 午夜精品国产一区二区电影 | 一级毛片电影观看 | 国产日韩欧美在线精品| 国产伦精品一区二区三区视频9| 性欧美人与动物交配| 1000部很黄的大片| 久久九九热精品免费| 欧美日韩一区二区视频在线观看视频在线 | 亚洲欧美中文字幕日韩二区| 午夜福利高清视频| 插阴视频在线观看视频| 全区人妻精品视频| 禁无遮挡网站| av天堂中文字幕网| 综合色丁香网| 精品久久久噜噜| 国产 一区精品| 国产精品嫩草影院av在线观看| 成人午夜精彩视频在线观看| 97人妻精品一区二区三区麻豆| 亚洲国产精品成人综合色| 最近最新中文字幕大全电影3| 日韩av不卡免费在线播放| 六月丁香七月| 国产伦在线观看视频一区| 久久久久久久久大av| a级毛色黄片| 国产一级毛片七仙女欲春2| 日韩欧美精品免费久久| 永久网站在线| 久久国内精品自在自线图片| 边亲边吃奶的免费视频| 最近2019中文字幕mv第一页| 亚洲精品久久久久久婷婷小说 | 成人性生交大片免费视频hd| 婷婷亚洲欧美| 日韩强制内射视频| 真实男女啪啪啪动态图| 美女国产视频在线观看| 精品久久久噜噜| 又粗又硬又长又爽又黄的视频 | 亚洲美女视频黄频| 国产亚洲av片在线观看秒播厂 | 变态另类成人亚洲欧美熟女| 美女cb高潮喷水在线观看| 你懂的网址亚洲精品在线观看 | 欧美成人免费av一区二区三区| 久久99精品国语久久久| 干丝袜人妻中文字幕| 亚洲精品亚洲一区二区| 亚洲一区高清亚洲精品| 人人妻人人澡欧美一区二区| 特大巨黑吊av在线直播| 国产一区二区三区av在线 | 在线a可以看的网站| 国产一区二区亚洲精品在线观看| 男人狂女人下面高潮的视频| 国产精品一区二区性色av| 99国产精品一区二区蜜桃av| 久久这里有精品视频免费| 久久99热这里只有精品18| 美女高潮的动态| 国内少妇人妻偷人精品xxx网站| 99在线视频只有这里精品首页| 日本成人三级电影网站| 亚洲av第一区精品v没综合| 成年av动漫网址| 波多野结衣高清无吗| 最近手机中文字幕大全| 色综合色国产| 欧美zozozo另类| 国产精品1区2区在线观看.| av卡一久久| 国产熟女欧美一区二区| 简卡轻食公司| 久久人人精品亚洲av| 99热全是精品| 色哟哟哟哟哟哟| 六月丁香七月| av.在线天堂| 1024手机看黄色片| 男人舔奶头视频| 99热这里只有是精品在线观看| av在线观看视频网站免费| 国产日本99.免费观看| 最好的美女福利视频网| 极品教师在线视频| 91麻豆精品激情在线观看国产| 亚洲天堂国产精品一区在线| 色综合亚洲欧美另类图片| 中出人妻视频一区二区| 成人国产麻豆网| 久久人妻av系列| 亚洲综合色惰| 国产av麻豆久久久久久久| 舔av片在线| 成人高潮视频无遮挡免费网站| 亚洲av第一区精品v没综合| 人妻制服诱惑在线中文字幕| 18禁裸乳无遮挡免费网站照片| 精品不卡国产一区二区三区| 午夜精品在线福利| 一进一出抽搐gif免费好疼| 熟女人妻精品中文字幕| 一区二区三区高清视频在线| 欧美成人一区二区免费高清观看| 伦精品一区二区三区| 一级毛片电影观看 | 免费人成在线观看视频色| a级毛片免费高清观看在线播放| 国产日韩欧美在线精品| 中国国产av一级| 岛国毛片在线播放| 亚洲高清免费不卡视频| 熟妇人妻久久中文字幕3abv| 精品久久国产蜜桃| 国产探花在线观看一区二区| 国产成年人精品一区二区| 亚洲久久久久久中文字幕| 国产不卡一卡二| 日韩一区二区三区影片| 亚洲国产高清在线一区二区三| 两个人视频免费观看高清| 国产精品伦人一区二区| 国内精品美女久久久久久| 国产成人aa在线观看| 国产真实伦视频高清在线观看| 国产一级毛片在线| 色吧在线观看| kizo精华| 亚洲婷婷狠狠爱综合网| 欧美一区二区精品小视频在线| 熟女人妻精品中文字幕| 精品久久久久久久久亚洲| 国产又黄又爽又无遮挡在线| 淫秽高清视频在线观看| 熟妇人妻久久中文字幕3abv| av卡一久久| 成人毛片60女人毛片免费| 桃色一区二区三区在线观看| 日韩精品有码人妻一区| 国内精品久久久久精免费| 深夜精品福利| 色5月婷婷丁香| 亚洲欧美日韩卡通动漫| 亚洲一区二区三区色噜噜| 久久久国产成人免费| 久久99精品国语久久久| 国产乱人偷精品视频| 国产精品免费一区二区三区在线| 久久99热6这里只有精品| 国产精品综合久久久久久久免费| 一本久久精品| 插逼视频在线观看| 日韩 亚洲 欧美在线| 亚洲国产精品sss在线观看| 亚洲图色成人| 亚洲国产欧美人成| 国产精品.久久久| 亚洲国产色片| 久久精品国产自在天天线| 国产成人a区在线观看| 欧美性猛交╳xxx乱大交人| 波野结衣二区三区在线| 国产麻豆成人av免费视频| 国产 一区精品| 国产高清三级在线| 久久精品影院6| 乱码一卡2卡4卡精品| 久久这里只有精品中国| 女人被狂操c到高潮| 91久久精品电影网| 在线免费十八禁| 身体一侧抽搐| videossex国产| 久久久久免费精品人妻一区二区| 岛国毛片在线播放| 99九九线精品视频在线观看视频| 亚洲丝袜综合中文字幕| 又粗又硬又长又爽又黄的视频 | 欧美激情国产日韩精品一区| 黄色视频,在线免费观看| 看黄色毛片网站| 亚洲自偷自拍三级| 亚洲av成人精品一区久久| 免费电影在线观看免费观看| 亚洲国产精品国产精品| 色综合亚洲欧美另类图片| 麻豆av噜噜一区二区三区| 人人妻人人澡欧美一区二区| 简卡轻食公司| 麻豆久久精品国产亚洲av| 一区二区三区免费毛片| 国产精品久久久久久久电影| 午夜老司机福利剧场| .国产精品久久| 亚洲精品国产av成人精品| 欧美日韩综合久久久久久| 男人的好看免费观看在线视频| 欧美区成人在线视频| 国产亚洲av嫩草精品影院| 在线观看66精品国产| 美女内射精品一级片tv| 特级一级黄色大片| 日韩成人av中文字幕在线观看| 成人永久免费在线观看视频| 免费大片18禁| 麻豆一二三区av精品| 深夜精品福利| 国产探花在线观看一区二区| 久久精品综合一区二区三区| 在线播放无遮挡| 国国产精品蜜臀av免费| kizo精华| 中文字幕熟女人妻在线| 联通29元200g的流量卡| 国产精品一区二区三区四区久久| 亚洲国产精品久久男人天堂| 在线观看一区二区三区| 国产日本99.免费观看| 亚洲精品乱码久久久久久按摩| 国内精品宾馆在线| 欧美日韩综合久久久久久| 好男人视频免费观看在线| 看十八女毛片水多多多| 卡戴珊不雅视频在线播放| 波野结衣二区三区在线| 岛国在线免费视频观看| 老熟妇乱子伦视频在线观看| 久久韩国三级中文字幕| 在线国产一区二区在线| 草草在线视频免费看| 又爽又黄无遮挡网站| 国产精品无大码| 熟女人妻精品中文字幕| 我的老师免费观看完整版| 老司机影院成人| 天堂中文最新版在线下载 | 男人舔奶头视频| 国产中年淑女户外野战色| 男女啪啪激烈高潮av片| 亚洲成人久久爱视频| 日日啪夜夜撸| 亚洲精品456在线播放app| av免费观看日本| 精品日产1卡2卡| 亚洲欧美精品专区久久| 男人狂女人下面高潮的视频| 国产不卡一卡二| 少妇熟女欧美另类| av天堂在线播放| 成年免费大片在线观看| 亚洲乱码一区二区免费版| 日本撒尿小便嘘嘘汇集6| 国产av在哪里看| 亚洲av免费高清在线观看| 欧美日本视频| 亚洲国产欧美人成| 中文字幕av在线有码专区| 国产精品爽爽va在线观看网站| 精品国内亚洲2022精品成人| 久久久久久久久久黄片| 国产综合懂色| 亚洲内射少妇av| 人人妻人人看人人澡| 日本免费a在线| 日韩国内少妇激情av| 91久久精品国产一区二区成人| 免费无遮挡裸体视频| 日韩av在线大香蕉| 久久久久国产网址|