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

    模擬氮沉降對三江平原小葉章濕地土壤細(xì)菌多樣性的影響

    2016-05-25 08:45:24張榮濤楊立賓柴春榮王繼豐付曉玲鐘海秀邢軍會(huì)倪紅偉
    草業(yè)科學(xué) 2016年4期
    關(guān)鍵詞:三江平原小葉菌門

    隋 心,張榮濤,楊立賓,許 楠,柴春榮,王繼豐,付曉玲,鐘海秀,邢軍會(huì),張 玉,倪紅偉

    (1.東北林業(yè)大學(xué)林學(xué)院, 黑龍江 哈爾濱 150040;2.黑龍江省科學(xué)院自然與生態(tài)研究所濕地與生態(tài)保育國家地方聯(lián)合工程實(shí)驗(yàn)室, 黑龍江 哈爾濱 150040)

    ?

    模擬氮沉降對三江平原小葉章濕地土壤細(xì)菌多樣性的影響

    隋 心1,2,張榮濤2,楊立賓2,許 楠2,柴春榮2,王繼豐2,付曉玲2,鐘海秀2,邢軍會(huì)2,張 玉2,倪紅偉2

    (1.東北林業(yè)大學(xué)林學(xué)院, 黑龍江 哈爾濱 150040;2.黑龍江省科學(xué)院自然與生態(tài)研究所濕地與生態(tài)保育國家地方聯(lián)合工程實(shí)驗(yàn)室, 黑龍江 哈爾濱 150040)

    摘要:利用黑龍江省科學(xué)院自然與生態(tài)研究所三江平原濕地生態(tài)定位研究站內(nèi)的模擬氮沉降試驗(yàn)平臺(tái),應(yīng)用Miseq測序技術(shù)對細(xì)菌16S rDNA進(jìn)行序列測定和分析,探討添加氮對三江平原小葉章濕地土壤細(xì)菌群落結(jié)構(gòu)的影響。結(jié)果表明,與對照相比,添加氮增加土壤細(xì)菌群落多樣性,低氮(4 g·hm-2·a-1)時(shí)土壤細(xì)菌多樣性最高,高氮(8 g·hm-2·a-1)次之。模擬氮沉降導(dǎo)致濕地土壤細(xì)菌群落的結(jié)構(gòu)以及豐度發(fā)生變化,低氮、高氮和對照地的優(yōu)勢種群均為酸桿菌門(Acidobacteria),變形菌門(Proteobacteria)次之,隨著氮沉降的增加,酸桿菌門的豐度逐漸下降,變形菌門的豐度則先升高后降低。外源氮的輸入使得土壤的理化性質(zhì)發(fā)生了改變,導(dǎo)致細(xì)菌群落的組成以及細(xì)菌的相對豐度發(fā)生改變。16S rDNA heatmap分析則表明,不同的細(xì)菌群落對氮沉降的響應(yīng)不同。本研究為預(yù)測未來大氣氮沉降對三江平原濕地土壤微生物的影響以及濕地生態(tài)系統(tǒng)的變化提供基礎(chǔ)數(shù)據(jù)和理論參考。

    關(guān)鍵詞:16S rDNA;小葉章濕地;Miseq;氮沉降;細(xì)菌多樣性

    近幾十年來,由于科技進(jìn)步以及工業(yè)的高速發(fā)展,燃燒大量的化石燃料導(dǎo)致了大氣氮沉降量迅速增加。根據(jù)歐洲NITREX(Nitrogen Saturation Experiments)項(xiàng)目計(jì)算,森林每年氮飽和的臨界負(fù)荷的最小值為10 kg·hm-2[1],但是到目前為止,中歐森林每年的大氣氮沉降量達(dá)到了25~60 kg·hm-2,已經(jīng)嚴(yán)重超出了森林的氮飽和的最小值[2]。從全球來看,在北美洲有些森林每年的大氣氮沉降量也超出了森林的氮飽和的最小值,達(dá)到40 kg·hm-2[3-4]。已有研究表明,我國降水中NO3-離子的濃度(3.2×10-7mol·L-1)和沉降量(2.0×10-7mol·L-1)與美國和日本接近,但是我國降水中的NH4+離子的濃度則更高,其沉降量也更大,離子濃度高達(dá)3.77×10-6mol·L-1,分別是美國和日本的4倍和3.7倍,導(dǎo)致我國氮沉降量十分巨大,分別是美國和日本的3倍和1.8倍[5]。目前,大氣氮沉降已經(jīng)成為影響全球變化的重要現(xiàn)象之一,受到了全世界的關(guān)注。近幾十年來,由于大氣氮沉降的沉降量不斷增加,給全球生態(tài)系統(tǒng)造成了嚴(yán)重的影響,例如土壤酸化的加重,生物功能和多樣性的降低等[6-8],引起了全世界的廣泛關(guān)注[9-10]。

    近年來,隨著分子生物學(xué)技術(shù)的發(fā)展,環(huán)境微生物的研究方法也從以前的實(shí)驗(yàn)室純培養(yǎng)進(jìn)入到了分子生物學(xué)方法階段。Handelsman等[11]在1998提出宏基因組(Metagenome)的概念,它包括環(huán)境中所有微生物的基因,目前主要指細(xì)菌和真菌的基因組總和。宏基因組學(xué)(Metagenomics)作為分子生態(tài)學(xué)的一種典型技術(shù)策略,能夠從環(huán)境樣品中直接獲得的基因組中進(jìn)行文庫構(gòu)建、克隆、測序[12],可以從群落水平上了解微生物完整的群落結(jié)構(gòu),并且真實(shí)地揭示環(huán)境微生物群落的多樣性和復(fù)雜性[13],在環(huán)境微生物群落結(jié)構(gòu)的研究中顯示出無比的優(yōu)越性[14],在土壤、水體、濕地以及極端環(huán)境微生物多樣性等研究方面被廣泛應(yīng)用[15-18]。

    三江平原濕地是我國面積最大、分布最為集中的濕地[19],由于人類活動(dòng)的影響,濕地遭到破壞,面積大幅減少,已經(jīng)從1949年的534.5萬 hm2減少到148.16萬hm2[9], 濕地生態(tài)系統(tǒng)結(jié)構(gòu)與功能發(fā)生了退化。目前研究主要集中在植物群落變化[10]、溫室氣體排放[19]和營養(yǎng)物質(zhì)循環(huán)[20]。國外大量的研究已經(jīng)表明,氮沉降長期增加會(huì)對土壤微生物產(chǎn)生不利影響[21-25],但是有關(guān)氮沉降對三江平原小葉章(Deyeuxiaangustifolia)濕地的土壤微生物的結(jié)構(gòu)組成和變化的影響,到目前為止還沒有相關(guān)的研究。因此,本研究通過野外原位試驗(yàn),探討模擬氮沉降對三江平原小葉章濕地土壤細(xì)菌多樣性的影響,以期為研究三江平原濕地生態(tài)系統(tǒng)在氮沉降背景下的變化提供科學(xué)依據(jù),并為三江平原濕地的保護(hù)與管理提供理論參考。

    1材料與方法

    1.1研究區(qū)域概況

    研究樣點(diǎn)設(shè)置在黑龍江省科學(xué)院自然與生態(tài)研究所三江平原濕地生態(tài)定位研究站——洪河國家級自然保護(hù)區(qū)內(nèi)(47°45′39″ N,133°37′04″ E)。研究站海拔55~65 m,多年平均氣溫1.9 ℃,≥10 ℃有效積溫2 165~2 624 ℃·d。多年平均降水量為585 mm,50%~70%集中在7-9月,多年平均蒸發(fā)量為1 166 mm。主要植被類型為草甸和沼澤,優(yōu)勢植物有小葉章、狹葉甜茅(Glyceriaspiculosa)、毛果苔草(Carexlasiocarpa)、漂筏苔草(Carexpseudocuraica)等[26]。

    1.2樣品采集

    2010 年5 月在研究區(qū)設(shè)置9塊20 m×20 m的小葉章沼澤化草甸濕地,設(shè)3個(gè)氮沉降梯度,分別為對照N1 (0 g·hm-2·a-1)、低氮N2(4 g·hm-2·a-1) 和高氮N3(8 g·hm-2·a-1),每種處理重復(fù)3次。每年5月份進(jìn)行施氮處理,將NH4NO3溶于水,用花灑均勻噴施。于2014年6月采集土壤樣品,采取隨機(jī)取樣的方式,在每個(gè)試驗(yàn)區(qū)分別用直徑為4 cm的土鉆按混合采樣法采集0-10 cm 表層土,每個(gè)土樣有5~8 個(gè)采集點(diǎn)的土壤混合,每個(gè)處理分別混合成一袋土,共計(jì)3袋,去除土樣中的植物和動(dòng)物殘?bào)w等雜質(zhì),混合均勻后用四分法取適量于自封袋中,迅速帶回實(shí)驗(yàn)室。取一部分保存在-80 ℃冰箱,用于分析土壤細(xì)菌群落結(jié)構(gòu)。另一部分風(fēng)干、磨細(xì)過0.15 mm篩、1.00 mm篩,用于測定土壤理化性質(zhì)。

    1.3土壤理化性質(zhì)檢測

    土壤有機(jī)碳利用碳氮分析儀測定(耶拿-2100S, 德國);土壤中全氮采用半微量凱氏定氮法測定;NO3--N采用酚二磺酸比色法測定;NH4+-N采用氯化鉀浸提-靛酚藍(lán)比色法測定。

    1.4土壤微生物總DNA的提取

    土壤微生物總DNA提取采用美國MO BIO的強(qiáng)力土壤DNA提取試劑盒(PowerSoil?DNA Isolation Kit),按照說明書的操作步驟進(jìn)行提取。將提取出的土壤微生物總DNA溶解在100 μL去離子水中,取5 μL的DNA用1.0%的瓊脂糖凝膠電泳檢測(0.5×TAE緩沖液),分析DNA的完整性和相對濃度。

    1.5土壤細(xì)菌16S rDNA基因測序

    由上海美吉生物醫(yī)藥科技有限公司進(jìn)行Miseq高通量測序。細(xì)菌16S rDNA擴(kuò)增引物采用通用引物(8F/533R)。

    1.6數(shù)據(jù)處理

    1.6.1有效序列和優(yōu)化序列數(shù)據(jù)統(tǒng)計(jì)由于采用多個(gè)樣品平行測序的方法,所以各樣品中的序列均引入了一段標(biāo)示其樣本來源信息的barcode標(biāo)簽序列及前引物(Forward Primer)序列。本次分析根據(jù)barcode標(biāo)簽序列和前引物序列篩選出有效序列后,將測序接頭與barcode序列去除,并對處理后的有效序列進(jìn)行數(shù)據(jù)統(tǒng)計(jì)。

    為了分析結(jié)果更精確,需要對有效序列進(jìn)行去雜。去掉小于150 bp和含有模糊堿基以及引物堿基對錯(cuò)配達(dá)到兩個(gè)以上的序列,從而得到優(yōu)化后的序列,用優(yōu)化后的序列來進(jìn)行后續(xù)生物信息分析。

    1.6.2OTU-based 分析對所有樣品進(jìn)行劃分可操作分類單元(Operationaltaxonomic unit,OTU ),并對OTU進(jìn)行生物信息統(tǒng)計(jì)分析。將優(yōu)化序列選取長度大于350 bp的序列并截齊后與silva數(shù)據(jù)庫中的微生物基因序列進(jìn)行比對后對序列進(jìn)行聚類。聚類分析使用的軟件為Mothur和Chopseq(http://www.mothur.org/wiki/Main-Page)。

    1.6.3細(xì)菌群落多樣性分析細(xì)菌群落物種的豐富度(Richness)和多樣性(Diversity)分別用Chao1和Shannon指數(shù)表示,測序深度指數(shù)用覆蓋度(Coverage)表示。計(jì)算公式如下:

    Chao1=OTUs+n1(n1-1)/2(n2+1).

    式中,OTUs=實(shí)際觀測到的OTU數(shù)目;n1=只含有一條序列OUT數(shù)目;n2=只含有兩條序列的OUT數(shù)目。

    Shannon index(H)=-∑Pi×log2Pi.

    式中,Pi=ni/N,表明第i個(gè)物種的相對多度;ni=含有i條序列的OTU數(shù)目;N=所有序列數(shù)。

    Coverage=1-n1/N.式中,n1=只含有一條序列的OTU數(shù)目,N=總序列數(shù)。

    群落多樣性分析使用的軟件為Mothur (version v.1.30.1) 。

    1.6.4群落結(jié)構(gòu)分析

    將優(yōu)化序列與silva庫中的已知序列進(jìn)行比對,對OTU序列進(jìn)行種屬鑒定。在門(Phylum)的水平上做樣品群落分布柱狀圖,比較3個(gè)土樣中的細(xì)菌群落分布情況。

    2結(jié)果

    2.1模擬氮沉降條件下土壤理化性質(zhì)的變化特征

    模擬氮沉降樣地的土壤理化性質(zhì)分析結(jié)果顯示(表1),樣地的土壤類型主要為草甸土,且不同處理間的pH值存在差異,其中高氮樣地的pH均值最低(P<0.05)。不同處理間土壤的全氮、硝態(tài)氮、銨態(tài)氮和速效磷存在顯著差異(P<0.05)。有機(jī)碳、全氮、硝態(tài)氮、全磷和銨態(tài)氮含量呈現(xiàn)隨著氮沉降的增加而升高的趨勢,而速效磷含量則呈現(xiàn)隨著氮沉降的增加而下降的趨勢。

    2.2測序數(shù)據(jù)分析

    通過Miseq高通量測序并優(yōu)化后,3個(gè)處理共獲得39 979條序列,總堿基數(shù)為15 845 087 bp,平均堿基長度為396.32 bp,其中301~400 bp長度的序列占總序列數(shù)的99.74%(表2)。

    表1 模擬氮沉降條件下的土壤理化性質(zhì)的比較

    注:同行不同小寫表示不同處理間差異顯著(P<0.05)。下同。

    Note:Different lower case letters within the same row show significant difference among different treatments at 0.05 level. The same below.

    表2 16s rDNA有效序列的數(shù)量及分布

    利用測序獲得的序列進(jìn)行隨機(jī)抽樣,將抽到的序列數(shù)所代表的OTU數(shù)目構(gòu)建稀釋性曲線。3個(gè)處理的稀釋性曲線均趨于平坦,表明測序數(shù)據(jù)接近飽和,測序的深度比較合理,增加測序數(shù)量對發(fā)現(xiàn)新的OTU貢獻(xiàn)率較小(圖1)。

    圖1 不同處理土壤細(xì)菌OTU稀釋曲線

    2.3模擬氮沉降條件下土壤細(xì)菌多樣性分析

    通過比對細(xì)菌16S rDNA的相似性,將相似性高于98%的序列歸為同一個(gè)OTU[27]。對不同氮沉降梯度下小葉章濕地土壤細(xì)菌16S rDNA多樣性指數(shù)進(jìn)行分析(表3),N1、N2和N3的Shannon指數(shù)分別為5.26、5.38和5.32,豐富度指數(shù)Chaol分別為745、698和728。與N1相比,隨著氮濃度的逐漸增加,細(xì)菌的多樣性有所增多,在N2時(shí)細(xì)菌多樣性達(dá)到最高;而細(xì)菌的豐富度指數(shù)與多樣性的變化不同,在N2時(shí)細(xì)菌豐富度最低。這表明不同氮沉降梯度下小葉章濕地土壤細(xì)菌的結(jié)構(gòu)和豐富度均存在著差異, 且不同氮沉降梯度下土壤細(xì)菌的群落結(jié)構(gòu)變化規(guī)律不同。

    表3 模擬氮沉降各梯度下小葉章濕地土壤

    2.4模擬氮沉降條件下土壤細(xì)菌16S rDNA序列分析

    將測序所得16S rDNA序列在Genbank中進(jìn)行比對,得到不同氮沉降梯度下小葉章濕地細(xì)菌的分類(表4)。N1、N2和N3小葉章濕地優(yōu)勢種群均為酸桿菌門(Acidobacteria),分別占總測序數(shù)的53.31%、45.64%和36.33%。變形菌門(Proteobacteria)分布僅次于酸桿菌門,在N1、N2和N3中也有大量的分布,分別占總序列數(shù)的21.60%、21.98%和18.20%。另外,浮霉菌門(Planctomycetes)、綠彎菌門(Chloroflexi)和硝化螺旋菌門(Nitrospirae)在3個(gè)處理中的比例均超過3%,也有大量的分布。

    表4 模擬氮沉降條件下土壤細(xì)菌16S rDNA序列分析

    不同氮沉降的濕地中土壤微生物還包括擬桿菌門(Bacteroidetes)、綠菌門(Chlorobi)、芽單胞菌門(Gemmatimonadetes)、迷蹤菌門(Elusimicrobia)、螺旋體菌門(Spirochaetae)、疣微菌門(Verrucomicrobia)、厚壁菌門(Firmicutes)、藍(lán)藻菌門(Cyanobacteria)、裝甲菌門(Armatimonadetes)、熱胞菌門(Thermotogae)和纖維桿菌門(Fibrobacteres),它們在N1、N2和N3三個(gè)處理中分布均比較少,各類細(xì)菌占總序列數(shù)的比例都小于3%。另外,N1、N2和N3三個(gè)處理中均有≥17%的細(xì)菌在Genbank中沒有被明確的分類,屬于未分類的細(xì)菌。

    與N1相比,氮沉降對酸桿菌門數(shù)量分布的影響最大,隨著氮濃度的增加,土壤中酸桿菌門呈現(xiàn)逐漸下降的趨勢。氮濃度增加后,與N1相比,N2和N3的酸桿菌門分別減少了7.67和16.98個(gè)百分點(diǎn)。變形菌門的變化規(guī)律稍有不同,與N1相比,N2增加了0.38個(gè)百分點(diǎn),N3降低了3.4個(gè)百分點(diǎn)。此外,浮霉菌門、綠彎菌門和硝化螺旋菌門受氮沉降的影響也比較大。與N1相比,土壤氮含量增加后,N2的浮霉菌門增加了0.72個(gè)百分點(diǎn),而N3下降了0.83個(gè)百分點(diǎn);而N2和N3的綠彎菌門分別增加了4.64和2.84個(gè)百分點(diǎn);硝化螺旋菌門則在N2中幾乎沒有變化,在N3中減少了0.49個(gè)百分點(diǎn)。

    2.5細(xì)菌群落結(jié)構(gòu)分析

    在屬的水平上對各處理所含菌屬進(jìn)行聚類,根據(jù)聚類后各處理中的不同OTU所含序列的豐度制作Heatmap(圖2),該圖能夠反映在菌屬水平上各樣品細(xì)菌群落結(jié)構(gòu)與豐度上的差異。不同氮沉降梯度下濕地土壤16S rDNA細(xì)菌可以劃分為5個(gè)聚類(Cluster)。菌屬的豐度受到不同氮濃度的影響,在所測得的菌屬中豐度受不同氮濃度影響比較大的是鏈霉菌(Streptomyces)、Oryzihumus和嗜酸鏈霉菌屬(Streptacidiphilus)。

    Cluster 1主要包括Bryobacter、Candidatus、硝化螺旋菌、鏈霉菌等(圖2)。氮沉降增加使得細(xì)菌結(jié)構(gòu)的分布發(fā)生了明顯的影響。例如,Bryobacter和Candidatus在N1中要高于N2和N3;鏈霉菌的豐度在N3中要高于N1和N2。

    Cluster 2主要包括綠彎菌門和硝化螺旋菌門兩大類群。其中綠彎菌門的豐度在N3中比較低,在N1和N2中較高;而硝化螺旋菌門在N3中比較高,在N1和N2中較低。

    Cluster 3主要包括螺旋體菌、放線菌、地桿菌、光合細(xì)菌、醋桿菌、Aquabacterium、Chitinophagaceae等。螺旋體菌在N3中的豐度較低,在N1和N2中豐度較高;放線菌和地桿菌則在N1中豐度較低,N2和N3中豐度較高;光合細(xì)菌和醋桿菌則在3個(gè)氮沉降樣地中豐度接近;Aquabacterium和Chitinophagaceae則在N1中豐度較高,在N2和N3中豐度較低。

    圖2 基于細(xì)菌16S rDNA序列構(gòu)建的熱圖

    Cluster 4主要包括Isosphaera、Chitinophagaceae、Syntrophorhabdaceae、黃單胞菌等。Isosphaera和Chitinophagaceae在N1中豐度較高,N2和N3中豐度較低;Syntrophorhabdaceae和黃單胞菌在N1中的豐度較低,N2和N3中的豐度較高。

    Cluster 5主要包括Singulisphaera、黏球菌、酸微菌等。Singulisphaera、黏球菌、酸微菌在N3中的豐度要大于N1和N2中。

    從總體來看,N1和N2的細(xì)菌結(jié)構(gòu)和豐度都比較接近,N3的土壤細(xì)菌群落結(jié)構(gòu)和豐度變化比較大,這說明N3對土壤細(xì)菌結(jié)構(gòu)的影響比較明顯,這從圖2可以明顯看出來。另外,對于N1和N2兩個(gè)小葉章濕地來看,二者的細(xì)菌群落結(jié)構(gòu)和豐度都比較接近,這說明在N2的條件下,土壤細(xì)菌的結(jié)構(gòu)并沒有發(fā)生顯著的改變。

    3討論

    土壤pH值是影響土壤微生物多樣性的重要因素之一,pH值的變化對土壤的有機(jī)質(zhì)分解與合成以及N、P、K營養(yǎng)元素的合成與轉(zhuǎn)化有著重要的作用[28]。在不同氮沉降條件下,N1的pH值最高,其次是N2,但二者無顯著差異(P>0.05)。氮沉降的增加會(huì)對土壤養(yǎng)分產(chǎn)生影響。將N2和N3兩種處理與N1相比,氮沉降使土壤有機(jī)碳、氮等的含量均顯著增加(表1),這可能是氮輸入能促進(jìn)碳固定,從而增加土壤有機(jī)碳含量,這與Conant等[29]的觀點(diǎn)一致。說明外源氮素的輸入增加了濕地生態(tài)系統(tǒng)可利用氮素狀況,從而對濕地土壤有機(jī)碳含量產(chǎn)生了重要影響。另外,外源氮素的增加還促進(jìn)了土壤硝態(tài)氮和銨態(tài)氮含量增加,促進(jìn)了土壤氮的礦化率,使得土壤礦質(zhì)氮的含量也相應(yīng)的增加,從而使土壤全氮的含量得到顯著增加,這和大多數(shù)學(xué)者的研究結(jié)果一致,表明土壤硝態(tài)氮和銨態(tài)氮在土壤氮飽和之前會(huì)隨著氮輸入量的增加而增加[30]。在南亞熱帶森林土壤有效氮含量對于模擬氮沉降的初期響應(yīng)的研究中,發(fā)現(xiàn)大氣的氮沉降會(huì)使土壤中的速效氮含量增加,而且氮處理水平越高、土壤速效氮含量增加越多[31]。

    多樣性指數(shù)對于分析不同土壤的微生物群落結(jié)構(gòu)是非常有效的方法,多樣性指數(shù)越高則微生物群落多樣性越高,它由種類的豐富度及種類的均勻度兩部分組成[32-33]。本研究表明,不同氮沉降條件下土壤細(xì)菌的群落結(jié)構(gòu)發(fā)生了顯著變化,Shannon指數(shù)呈現(xiàn)出N2>N3>N1的趨勢(表2),可見,氮沉降會(huì)增加速效氮含量增加,而且氮處理水平越高、土壤速效氮含量增加越多[31]。

    多樣性指數(shù)對于分析不同土壤的微生物群落結(jié)構(gòu)是非常有效的方法,多樣性指數(shù)越高則微生物群落多樣性越高,它由種類的豐富度及種類的均勻度兩部分組成[32-33]。本研究結(jié)果表明,不同氮沉降條件下土壤細(xì)菌的群落結(jié)構(gòu)發(fā)生了顯著的變化,Shannon指數(shù)呈現(xiàn)出N2>N3>N1的趨勢(表2),可見,氮沉降會(huì)增加細(xì)菌的多樣性,低氮較高,高氮略有下降,說明氮沉降增加細(xì)菌多樣性具有閾值限制性,即氮濃度增加到一定閾值時(shí),細(xì)菌多樣性則由升高轉(zhuǎn)變?yōu)橄陆?。這可能是由于氮沉降增加了土壤中的營養(yǎng)物質(zhì),為細(xì)菌的生長提供了豐富的營養(yǎng)來源,因此,一定濃度的氮素會(huì)促進(jìn)細(xì)菌的生長,這與對杉木林土壤微生物[34]的研究結(jié)果一致。但是在高氮水平下,氮增加改變了土壤中營養(yǎng)物質(zhì)的有效性,這導(dǎo)致了某些微生物定閾值時(shí),細(xì)菌多樣性則由升高轉(zhuǎn)變?yōu)橄陆?。這可能是由于氮沉降增加了土壤中的營養(yǎng)物質(zhì),為細(xì)菌的生長提供了豐富的營養(yǎng)來源,因此,一定濃度的氮素會(huì)促進(jìn)細(xì)菌的生長,這和劉彩霞等[34]的研究結(jié)果一致。但是在高氮水平下,氮增加改變了土壤中營養(yǎng)物質(zhì)的有效性,這導(dǎo)致了某些微生物種群(如偏于利用植物殘?bào)w的種群) 的生長,但是同時(shí)又抑制了其它種群,致使群落多樣性指數(shù)的下降??偟膩碚f,氮沉降會(huì)對土壤細(xì)菌的群落結(jié)構(gòu)具有一定促進(jìn)作用。由此可以推斷,目前我國一些地區(qū)的高氮沉降勢必會(huì)對土壤微生物群落產(chǎn)生影響[35-37]。

    通過16S rDNA序列分析表明,N3、N2和N1的優(yōu)勢種群都是酸桿菌門,分別占53.31%、45.64%和36.33%。很多研究都表明,濕地土壤中主要的細(xì)菌群落是酸桿菌門和變形菌門。例如,三江平原濕地土壤細(xì)菌主要是酸桿菌[38],還有研究表明,濕地土壤中主要優(yōu)勢細(xì)菌群落為變形菌[39-40],本研究得到相似結(jié)果,表明在濕地土壤中,酸桿菌門和變形菌門是最主要的兩大類群。

    氮沉降改變了土壤細(xì)菌群落組成,酸桿菌門的豐度隨著氮沉降的增加而下降,而變形菌門隨著氮沉降的增加先升高后降低,這和Zhou等[41]的研究結(jié)果一致,表明酸桿菌門與土壤氮含量呈負(fù)相關(guān)關(guān)系。對于變形菌門的豐度來說,低氮處理時(shí),變形菌門豐度最高,高氮時(shí)出現(xiàn)下降。有一些研究結(jié)果表明,變形菌門在長期施氮處理中的占主要優(yōu)勢[42-43]。而本研究發(fā)現(xiàn),在高氮處理時(shí),土壤變形菌門豐度會(huì)下降,這可能是由于在高氮濃度時(shí),土壤氮含量已經(jīng)過飽和,過飽和氮反而會(huì)抑制土壤養(yǎng)分的供給,導(dǎo)致變形菌豐度下降。

    本研究結(jié)果還顯示,氮沉降導(dǎo)致土壤礦質(zhì)氮含量大幅上升,土壤呈現(xiàn)弱酸性,而3個(gè)處理中都含有較多的酸桿菌。這可能是與弱酸性的土壤環(huán)境有關(guān),土壤的弱酸性會(huì)有利于某些酸桿菌的生長。濕地中土壤微生物在濕地土壤的物質(zhì)循環(huán)和能量流動(dòng)中起非常重要的作用。但目前對濕地的研究還相對較少,對氮沉降的研究更為薄弱。這是因?yàn)闈竦嘏c森林、草原等生態(tài)系統(tǒng)相比,由于其環(huán)境特殊,季節(jié)性或常年積水,所以在不同地點(diǎn)、不同季節(jié)的研究結(jié)果差別會(huì)較大。因此,研究濕地土壤微生物是一個(gè)十分復(fù)雜和龐大的系統(tǒng)工程,對其研究需要不斷地從各個(gè)方面積累,并應(yīng)有多個(gè)學(xué)科的滲透與加入。全球大氣氮沉降增加背景下我國東北地區(qū)的氮沉降將會(huì)增加。在此背景下,研究濕地土壤細(xì)菌群落組成和結(jié)構(gòu)對氮沉降的響應(yīng),對研究未來氣候變化對濕地生態(tài)系統(tǒng)多樣性格局和功能穩(wěn)定性的影響具有重要參考意義。

    參考文獻(xiàn)(References)

    [1]Wright R F,Roelofs J G M,Bredemeier M,Blanck K,Boxman A W,Emmett B A,Tietema A.NITREX:Responses of coniferous forest ecosystems to experimentally changed deposition of nitrogen.Forest Ecology and Management,1995,71(1):163-169.

    [2]Kazda M.Indications of unbalanced nitrogen nutrition of Norway spruce stands.Plant and Soil,1990,128(1):97-101.

    [3]Lovett G M,Reiners W A,Olson R K.Cloud droplet deposition in subalpine balsam fir forests:Hydrological and chemical inputs.Science,1982,218(4579):1303-1304.

    [4]Parker G G.Throughfall and stemflow in the forest nutrient cycle.Advances in Ecological Research,1983,13(4):57-133.

    [5]孫崇基.酸雨.北京:中國環(huán)境科學(xué)出版社,2001.

    Sun C J.Acid Rain.Beijing:China Environmental Science Press,2001.(in Chinese)

    [6]Hofhansl F,Wanek W,Drage S,Huber W,Weissenhofer A,Richter M A.Topography strongly affects atmospheric deposition and canopy exchange processes in different types of wet lowland rainforest,Southwest Costa Rica.Biogeochemistry,2011,106(3):371-396.

    [7]薛璟花,莫江明,李炯,方運(yùn)霆,李德軍.氮沉降對外生菌根真菌的影響.生態(tài)學(xué)報(bào),2004,24(8):1785-1792.

    Xue J H,Mo J M,Li J,Fan Y T,Li D J.Effects of nitrogen deposition on ectomycorrhizal fungi.Acta Ecologica Sinica,2004,24(8):1785-1792.(in Chinese)

    [8]王汝南.模擬大氣氮沉降對溫帶森林土壤溫室氣體交換通量的影響.北京:北京林業(yè)大學(xué)碩士學(xué)位論文,2012.

    Wang R N.Effects of simulated atmospheric nitrogen deposition on the exchange fluxes of greenhouse gases in the temperate forest soil.Master Thesis.Beijing:Beijing Forestry University,2012.(in Chinese)

    [9]劉興土,馬學(xué)慧.人類大面積開荒對三江平原環(huán)境影響及區(qū)域性生態(tài)環(huán)境保護(hù).地理科學(xué),2000,20(1):14-19.

    Liu X T,Ma X H.Influence of large-scale reclamation on natural environment and regional environmental protection in the Sanjiang Plain.Scientia Geographica Sinica,2000,20(1):14-19.(in Chinese)

    [10]汲玉河,呂憲國,楊青,趙魁義.三江平原小葉章群落近三十年的動(dòng)態(tài)變化.生態(tài)學(xué)雜志,2006,25(11):1328-1332.

    Ji Y H,Lyu X G,Yang Q,Zhao K Y.Recent thirty years evolvement trend ofCalamagrostisangustifoliacommunity in Sanjiang Plain.Chinese Journal of Ecology,2006,25(11):1328-1332.(in Chinese)

    [11]Handelsman J,Rondon M R,Brady S F.Molecular biological access to the chemistry of unknown soil microbes:A new frontier for natural products.Chemistry Biology,1998,5(10):245-249.

    [12]秦楠,栗東芳,楊瑞馥.高通量測序技術(shù)及其在微生物學(xué)研究中的應(yīng)用.微生物學(xué)報(bào),2011,51(4):445-457.

    Qin N,Li D F,Yang R F.Next-generation sequencing technologies and the application in microbiology——A review.Acta Microbiologica Sinica,2011,51(4):445-457.(in Chinese)

    [13]李慶崗,陶立高.高通量測序技術(shù)及其在生命科學(xué)中的應(yīng)用.畜牧與飼料科學(xué),2012,33(2):25-28.

    Li Q G,Tao L G.High-throughput sequencing technology and its application in life science.Animal Husbandry and Feed Science,2012,33(2):25-28.(in Chinese)

    [14]艾芳芳,楊樺,曲媛媛,周集體,李昂,關(guān)曉燕,茍敏.宏基因組研究及其應(yīng)用研究進(jìn)展.環(huán)境科學(xué)與技術(shù),2007,30(12):100-103.

    Ai F F,Yang H,Qu Y Y,Zhou J T,Li A,Guan X Y,Gou M.Progress on metagenomic technology and its application.Environmental Science & Technology,2007,30(12):100-103.(in Chinese)

    [15]He Z L,Xu M,Deng Y,Kang S,Kellogg L,Wu L,Van Nostrand J D,Hobbie S E,Reich P B,Zhou J.Meta genomic analysis reveals a marked divergence in the structure of belowground microbial communities at elevated CO2.Ecology Letter,2010,13(5):564-575.

    [16]Monteiro J M,Vollú R E,Coelho M R R,Fonseca A,Neto S C G,Seldin L.Bacterial communities within the rhizosphere and roots of vetiver (ChrysopogonzizanioidesL. Roberty) sampled at different growth stages.European Journal of Soil Biology,2011,47(4):236-242.

    [17]Yu Y,Wang H,Liu J,Wang Q,Shen T,Guo W,Wang R.Shifts in microbial community function and structure along the successional gradient of coastal wetlands in Yellow River Estuary.European Journal of Soil Biology,2012,49(3):12-21.

    [18]侯梅鋒,何士龍,李棟,張潔,趙云.連云港海底底泥及青海湖底泥細(xì)菌多樣性研究.環(huán)境科學(xué),2011,32(9):2681-2688.

    Hou M F,He S L,Li D,Zhang J,Zhao Y.Bacterial diversity in Lianyungang marine sediment and Qinghai Lake sediment.Chinese Journal of Environmental Science,2011,32(9):2681-2688.(in Chinese)

    [19]宋長春,王毅勇,王躍思,趙志春.人類活動(dòng)影響下淡水沼澤濕地溫室氣體排放變化.地理科學(xué),2005,26(1):82-86.

    Song C C,Wang Y Y,Wang Y S,Zhao Z C.Character of the greenhouse gas emission in the freshwater mire under human activities.Scientia Geographica Sinica,2005,26(1):82-86.(in Chinese)

    [20]孫志高,劉景雙,于君寶.三江平原小葉章不同植物群落氮素的累積與分配.應(yīng)用生態(tài)學(xué)報(bào),2009,20(2):277-284.

    Sun Z G,Liu J S,Yu J B.Nitrogen accumulation and allocation in calamagrostis angusifolia in different plant communities of Sanjiang Plain.Chinese Journal of Applied Ecology,2009,20(2):277-284.(in Chinese)

    [21]Lilleskov E A,Fahey T J,Lovett G M.Ectomycorrhizal fungal aboveground community change over an atmospheric nitrogen deposition gradient.Ecological Applications,2001(2):397-410.

    [22]Lilleskov E A,Hobbie E A,Fahey T J.Ectomycorrhizal fungal taxa differing in response to nitrogen deposition also differ in pure culture organic nitrogen use and natural abundance of nitrogen isotopes.New Phytologist,2002,154(1):219-231(13).

    [23]Waldrop M P,Zak D R,Sinsabaugh R L.Microbial community response to nitrogendeposition in northern forest ecosystems.Soil Biology & Biochemistry,2004,36(9):1443-1451.

    [24]Compton J E,Watrud L S,Porteous L A,Degrood S.Response of soil microbial biomass and community composition to chronic nitrogen additions at Harvard forest.Forest Ecology & Management,2004,196(1):143-158.

    [25]Deforest J L,Zak D R,Pregitzer K S,Burton A.Atmospheric nitrate deposition and the microbial degradation of cellobiose and vanillin in a Northern Hardwood forest.Soil Biology & Biochemistry,2004,36(6):965-971.

    [26]倪紅偉,李君.洪河自然保護(hù)區(qū)生物多樣性.哈爾濱:黑龍江科學(xué)技術(shù)出版社,1999.

    Ni H W,Li J.Biodiversity in the Honghe National Nature Reserve.Harbin:Heilongjiang Science and Technology Press,1999.(in Chinese)

    [27]Stach J E,Maldonado L A,Masson D G,Ward A C,Goodfellow M,Bull A T.Staistical approaches for estimating actinobacterial diversity in marine sediments.Applied and Environmental Microbiology,2003,69(10):6189-6200.

    [28]Ndour N Y B,Achouak W,Christen R,Heulin T,Brauman A,Chotte J L.Characteristics of microbial habitats in a tropical soil subject to different fallow management.Applied Soil Ecology,2008,38(1):51-61.

    [29]Conant R T,Paustian K,Elliot E T.Grassand management and conversion into grassland:Effects on soil carbon.Ecological Application,2001,11(2):343-355.

    [30]樊后保,劉文飛,李燕燕,廖迎春,袁穎紅,徐雷.亞熱帶杉木(Cunninghamialanceolata)人工林生長與土壤養(yǎng)分對氮沉降的響應(yīng).生態(tài)學(xué)報(bào),2007,27(11):4630-4642.

    Fan H B,Liu W F,Li Y Y,Liao Y C,Yuan Y H,Xu L.Tree growth and soil nutrients in response to nitrogen deposition in a subtropical chinese fir plantation.Acta Ecologica Sinica,2007,27(11):4630-4642.(in Chinese)

    [31]王暉,莫江明,魯顯楷,薛璟花,李炯,方運(yùn)霆.南亞熱帶森林土壤微生物量碳對氮沉降的響應(yīng).生態(tài)學(xué)報(bào),2008,28(2):470-478.

    Wang H,Mo J M,Lu X K,Xue J H,Li J,Fang Y T.Effects of elevated nitrogen deposition on soil microbial biomass carbon in the main subtropical forests of southern China.Acta Ecologica Sinica,2008,28(2):470-478.(in Chinese)

    [32]Harch B D,Meech W,Pankhurst C E,Correll R L.Using the Gini coefficient with BIOLOG substrate utilization data to provide an alternative quantitative measure for comparing bacterial soil communities.Journal of Microbiological Methods,1997,30(1):91-101.

    [33]Staddon W J,Duchesne L C,Trevors J T.Microbial diversity and community srutcture of postdisturbance forest soil as determined by soil-carbon source utilization patterns.Microbial Ecology,1997,34(2):125-130.

    [34]劉彩霞,焦如珍,董玉紅,孫啟武,李峰卿,周新華.杉木林土壤微生物區(qū)系對短期模擬氮沉降的響應(yīng).林業(yè)科學(xué)研究,2015,28(2):271-276.

    Liu C X,Jiao R Z,Dong Y H,Sun Q W,Li F Q,Zhou X H.Response of soil microbe flora to simulated nitrogen deposition at a plantation ofCunninghamialanceolata.Forest Research,2015,28(2):271-276.(in Chinese)

    [35]侯愛敏,彭少麟,周國逸.鼎湖山地區(qū)馬尾松年輪元素含量與酸雨的關(guān)系.生態(tài)學(xué)報(bào),2002,22(9):1552-1559.

    Hou A M,Peng S L,Zhou G Y.Tree-ring chemical changes and possible impacts of acid precipitation in Dinghushan,south China.Acta Ecologica Sinica,2002,22(9):1552-1559.(in Chinese)

    [36]劉菊秀,周國逸,溫達(dá)志,褚國偉.酸沉降影響下廣東陸地生態(tài)系統(tǒng)表層土壤特征.農(nóng)業(yè)環(huán)境科學(xué)學(xué)報(bào),2001,20(4):231-234.

    Liu J X,Zhou G Y,Wen D Z,Chu G W.Properties of surface soil on selected forest ecosystem affected by acid deposition in Guangdong.Agro-environmental Protection,2001,20(4):231-234.(in Chinese)

    [37]周國逸,閆俊華.鼎湖區(qū)域大氣降水特征和物質(zhì)元素輸入對森林生態(tài)系統(tǒng)存在和發(fā)育影響.生態(tài)學(xué)報(bào),2001,21(12):2002-2012.

    Zhou G Y,Yan J H.The influences of regional atmospheric precipitation characteristics and its element inputs on the existence and development of Dinghushan forest ecosystems.Acta Ecologica Sinica,2001,21(12):2002-2012.(in Chinese)

    [38]劉峰,王建波,謝立紅,許修宏.三江平原根際微生物群落結(jié)構(gòu)研究.國土與自然資源研究,2013(3):93-95.

    Liu F,Wang J B,Xie L H,Xu X H.Community structure of rhizosphere microbial in Sanjiang Plain.Territory & Natural Resources Study,2013(3):93-95.(in Chinese)

    [39]劉紹雄,王明月,王娟,楊宇明,謬???王金華,熊智.基于PCR-DGGE技術(shù)的劍湖濕地湖濱帶土壤微生物群落結(jié)構(gòu)多樣性分析.農(nóng)業(yè)環(huán)境科學(xué)學(xué)報(bào),2013,32(7):1405-1412.

    Liu S X,Wang M Y,Wang J,Yang Y M,Miao F J,Wang J H,Xiong Z.Analyzing soil microbial community structure diversity from Jianhu wetland lakeside zone using PCR-DGGE technique.Journal of Agro-environment Science,2013,32(7),1405-1412.(in Chinese)

    [40]王文雯.若爾蓋高寒濕地土壤細(xì)菌群落結(jié)構(gòu)比較.北京:中國科學(xué)院研究生院碩士學(xué)位論文,2010.

    Wang W W.Comparison of soil bacterial community structure in Ruoergai high cold Wetland.Master Thesis.Beijing:Graduate University of the Chinese Academy of Sciences,2010.(in Chinese)

    [41]Zhou J,Guan D,Zhou B,Zhao B,Ma M,Qin J.Influence of 34-years of fertilization on bacterial communities in an intensively cultivated black soil in northeast China.Soil Biology & Biochemistry,2015,90:42-51.

    [42]Campbell B J,Polson S W,HansonT E,Mack M C,Schuur E A G.The effect of nutrient deposition on bacterial communities in arctic tundra soil.Environmental Microbiology,2010,12(7):1842-1854.

    [43]Jangid K,Williams M A,Franzluebbers A J,Sanderlin J S,Reeves J H,Jenkins M B.Relative impacts of land-use,management intensity and fertilization upon soil microbial community structure in agricultural systems.Soil Biology & Biochemistry,2008,40(11):2843-2853.

    Effect of simulation nitrogen depositions on bacterial diversity ofDeyeuxiaangustifoliain wetland of Sanjiang Plain

    Sui Xin1,2, Zhang Rong-tao2, Yang Li-bin2, Xu Nan2, Chai Chun-rong2, Wang Ji-feng2,Fu Xiao-ling2, Zhong Hai-xiu2, Xing Jun-hui2, Zhang Yu2, Ni Hong-wei2

    (1.College of Forestry, Northeast Forestry University, Harbin 150040, China;2.National and Provincial Joint Engineering Laboratory of Wetlands and Ecological Conservation Institute of Nature & Ecology, Heilongjiang Academy of Sciences, Harbin 150040, China)

    Abstract:Wetland soil bacterial community structure under the simulation nitrogen depositions was studied using 16S rDNA gene sequencing. The soil samples were collected from three different nitrogen depositions in the Sanjiang Plain Field Experiment Station of the Institute of Nature & Ecology, Heilongjiang Academy Sciences. Total microbial DNA was directly extracted from soil samples and 16S rDNA genes were sequenced and analyzed by Miseq molecular technology. The results showed that both low N and high N treatment increased the bacterial diversity compared with control treatment and low N treatment had the highest diversity. The wetland bacterial community structures significantly changed under simulation N depositions. In low N treatment, high N treatment and control treatment, Acidobacteria were predominated followed by Proteobacteria. The richness of Acidobacteria decreased and the richness of Proteobacteria increased first then decreased along the deposition. Soil pH decreased with deposition, which showing acidic property and increasing the Acidobacteria distribution. Soil physicochemical properties changed by input N which resulted in the soil bacterial community composition and abundant changed. The heatmap tree based on the 16S rDNA sequences showed that different bacterial community structure differently responded on Nitrogen deposition. This research provided basic data and theory consultation for predicting the influence of future atmospheric nitrogen deposition on soil microorganism in Sanjiang Plain wetlands.

    Key words:16S rDNA; wetland; miseq; nitrogen deposition; bacterial diversity

    Corresponding author:Ni Hong-weiE-mail:nihongwei2000@163.com

    中圖分類號(hào):S154.36

    文獻(xiàn)標(biāo)識(shí)碼:A

    文章編號(hào):1001-0629(2016)4-0589-10*

    通信作者:倪紅偉(1964-),男,黑龍江雙城人,研究員,博士,研究方向?yàn)樯鷳B(tài)學(xué)。E-mail:nihongwei2000@163.com

    基金項(xiàng)目:國家自然科學(xué)基金項(xiàng)目(31470019);黑龍江省科研機(jī)構(gòu)創(chuàng)新能力提升專項(xiàng)計(jì)劃(GY2014KC0039);黑龍江省科學(xué)院院級項(xiàng)目(STJB16-01、STJB16-04、2015YZ01、2014ST05)

    收稿日期:2015-08-06接受日期:2015-12-08

    DOI:10.11829/j.issn.1001-0629.2015-0444

    隋心,張榮濤,楊立賓,許楠,柴春榮,王繼豐,付曉玲,鐘海秀,邢軍會(huì),張玉,倪紅偉.模擬氮沉降對三江平原小葉章濕地土壤細(xì)菌多樣性的影響.草業(yè)科學(xué),2016,33(4):589-598.

    Sui X,Zhang R T,Yang L B,Xu N,Chai C R,Wang J F,Fu X L,Zhong H X,Xing J H,Zhang Y,Ni H W.Effect of simulation nitrogen depositions on bacterial diversity ofDeyeuxiaangustifoliain wetland of Sanjiang Plain.Pratacultural Science,2016,33(4):589-598.

    第一作者:隋心(1982-),男,黑龍江拜泉人,研究實(shí)習(xí)員,博士,研究方向?yàn)樯鷳B(tài)學(xué)。E-mail:xinsui_cool@126.com

    猜你喜歡
    三江平原小葉菌門
    Positive unlabeled named entity recognition with multi-granularity linguistic information①
    蘋果小葉病的防治
    河北果樹(2020年4期)2020-11-26 06:05:00
    野生樹鼩與人工飼養(yǎng)樹鼩消化道不同部位微生物組成的比較研究
    饑餓與重?cái)z食對河蟹腸道菌群結(jié)構(gòu)的影響
    小葉樟樹下的遐思
    海峽姐妹(2019年1期)2019-03-23 02:42:40
    昆蟲體內(nèi)微生物多樣性的影響因素研究進(jìn)展
    妊娠期糖尿病腸道菌群變化及臨床價(jià)值分析
    小葉蓮化學(xué)成分的研究
    中成藥(2018年9期)2018-10-09 07:18:46
    三江平原氣溫變化趨勢及區(qū)域差異識(shí)別
    基于灰色絕對關(guān)聯(lián)度的三江平原降水量影響因子識(shí)別
    国产男女超爽视频在线观看| 国产乱人偷精品视频| 秋霞在线观看毛片| 午夜激情福利司机影院| 天美传媒精品一区二区| 午夜激情久久久久久久| 中文精品一卡2卡3卡4更新| 精品午夜福利在线看| 高清av免费在线| 亚洲精品色激情综合| 国产一区二区在线观看日韩| 国产精品成人在线| 亚洲精品乱码久久久久久按摩| 熟女电影av网| 人人妻人人澡人人爽人人夜夜| 国产伦理片在线播放av一区| tube8黄色片| 夜夜爽夜夜爽视频| 久久婷婷青草| 欧美97在线视频| 亚洲av中文字字幕乱码综合| 国产一区有黄有色的免费视频| 老熟女久久久| 久久精品国产鲁丝片午夜精品| 纯流量卡能插随身wifi吗| 亚洲av在线观看美女高潮| 老女人水多毛片| 男人添女人高潮全过程视频| 亚洲三级黄色毛片| av视频免费观看在线观看| 中文资源天堂在线| 亚洲av中文字字幕乱码综合| 人妻 亚洲 视频| 久久av网站| 国产在线免费精品| 久久久久久久久久人人人人人人| 久久久久久久久久久免费av| 成人亚洲欧美一区二区av| 黄片wwwwww| 菩萨蛮人人尽说江南好唐韦庄| 国产精品av视频在线免费观看| 日韩精品有码人妻一区| 美女国产视频在线观看| 色吧在线观看| 中国三级夫妇交换| 九色成人免费人妻av| 一级a做视频免费观看| 青春草亚洲视频在线观看| 午夜免费鲁丝| 男男h啪啪无遮挡| 国语对白做爰xxxⅹ性视频网站| 亚洲欧美清纯卡通| 男人爽女人下面视频在线观看| 丝袜脚勾引网站| 亚洲成人一二三区av| 九九爱精品视频在线观看| av在线老鸭窝| 99久久精品一区二区三区| 小蜜桃在线观看免费完整版高清| 小蜜桃在线观看免费完整版高清| 亚洲欧美成人精品一区二区| 十八禁网站网址无遮挡 | 中文字幕亚洲精品专区| 尤物成人国产欧美一区二区三区| 新久久久久国产一级毛片| 美女中出高潮动态图| 日日撸夜夜添| 国产亚洲欧美精品永久| 99九九线精品视频在线观看视频| 男女国产视频网站| 日本av手机在线免费观看| 久久午夜福利片| 欧美精品亚洲一区二区| 特大巨黑吊av在线直播| 国产爽快片一区二区三区| 国产日韩欧美在线精品| 蜜桃在线观看..| 国产白丝娇喘喷水9色精品| 国产成人精品久久久久久| 街头女战士在线观看网站| 日本wwww免费看| 欧美精品国产亚洲| 视频中文字幕在线观看| 好男人视频免费观看在线| 色婷婷久久久亚洲欧美| 国产精品一二三区在线看| 国产在线一区二区三区精| 亚洲三级黄色毛片| 免费大片黄手机在线观看| 日本免费在线观看一区| 在线观看免费日韩欧美大片 | 观看免费一级毛片| 大又大粗又爽又黄少妇毛片口| 久久久欧美国产精品| 一个人看的www免费观看视频| 国产精品久久久久久久电影| 男女啪啪激烈高潮av片| 26uuu在线亚洲综合色| 亚洲精品乱久久久久久| 一级二级三级毛片免费看| 好男人视频免费观看在线| 少妇 在线观看| 九九爱精品视频在线观看| 久久婷婷青草| 国产成人精品福利久久| 午夜福利高清视频| 亚洲欧洲国产日韩| 精品久久国产蜜桃| 欧美zozozo另类| 国语对白做爰xxxⅹ性视频网站| 日本一二三区视频观看| 欧美xxxx性猛交bbbb| 亚洲av综合色区一区| 精品亚洲乱码少妇综合久久| 80岁老熟妇乱子伦牲交| 亚洲经典国产精华液单| 久久久亚洲精品成人影院| 久久久久精品久久久久真实原创| 97超视频在线观看视频| 国产视频首页在线观看| 好男人视频免费观看在线| 在线免费观看不下载黄p国产| 国产精品久久久久久精品古装| 久久精品国产鲁丝片午夜精品| 久久这里有精品视频免费| 成人免费观看视频高清| 中文字幕av成人在线电影| 日韩中字成人| 免费久久久久久久精品成人欧美视频 | 日韩欧美一区视频在线观看 | 日韩,欧美,国产一区二区三区| 免费av不卡在线播放| 国产乱人视频| 日韩不卡一区二区三区视频在线| 国产日韩欧美亚洲二区| 欧美一区二区亚洲| 国产黄片视频在线免费观看| 成年人午夜在线观看视频| 狠狠精品人妻久久久久久综合| 国产高清三级在线| 亚洲欧美成人精品一区二区| 全区人妻精品视频| 亚洲av国产av综合av卡| 国产成人免费观看mmmm| 精品少妇久久久久久888优播| 高清欧美精品videossex| 亚洲精品456在线播放app| 联通29元200g的流量卡| 欧美人与善性xxx| av在线蜜桃| 亚洲成人中文字幕在线播放| 国产成人精品婷婷| 国产av国产精品国产| 狂野欧美白嫩少妇大欣赏| 国产深夜福利视频在线观看| 91精品一卡2卡3卡4卡| 国产有黄有色有爽视频| 少妇人妻一区二区三区视频| 成人国产av品久久久| 丰满人妻一区二区三区视频av| 噜噜噜噜噜久久久久久91| 性色avwww在线观看| 成年免费大片在线观看| 中文字幕人妻熟人妻熟丝袜美| 熟女av电影| 亚洲精品aⅴ在线观看| 菩萨蛮人人尽说江南好唐韦庄| 国产视频首页在线观看| 免费观看的影片在线观看| 亚洲成人一二三区av| a级一级毛片免费在线观看| 国产成人精品福利久久| 欧美亚洲 丝袜 人妻 在线| xxx大片免费视频| 大片免费播放器 马上看| 免费看日本二区| 美女中出高潮动态图| 国产熟女欧美一区二区| 日韩伦理黄色片| 久久影院123| 人妻少妇偷人精品九色| 国产在线一区二区三区精| 十八禁网站网址无遮挡 | 国产女主播在线喷水免费视频网站| 最黄视频免费看| 舔av片在线| 国产精品嫩草影院av在线观看| 日韩一区二区三区影片| 男人添女人高潮全过程视频| 国产av码专区亚洲av| 搡女人真爽免费视频火全软件| 青春草亚洲视频在线观看| 黄色视频在线播放观看不卡| 成人美女网站在线观看视频| 亚洲成色77777| 婷婷色综合大香蕉| 搡老乐熟女国产| 久久久久网色| 国产 一区 欧美 日韩| 国产精品不卡视频一区二区| 1000部很黄的大片| 日韩欧美 国产精品| 日韩制服骚丝袜av| 丰满迷人的少妇在线观看| 乱系列少妇在线播放| 丝袜喷水一区| 亚洲av二区三区四区| 国产精品国产三级专区第一集| 97超视频在线观看视频| 亚洲欧洲日产国产| 热99国产精品久久久久久7| 久久精品国产鲁丝片午夜精品| 联通29元200g的流量卡| 成年美女黄网站色视频大全免费 | 另类亚洲欧美激情| 男女免费视频国产| 黑人高潮一二区| 日韩免费高清中文字幕av| av线在线观看网站| 国产精品一区二区三区四区免费观看| 男的添女的下面高潮视频| 久久人人爽人人爽人人片va| 我要看黄色一级片免费的| 尤物成人国产欧美一区二区三区| 日日摸夜夜添夜夜添av毛片| 久久热精品热| 免费av中文字幕在线| 国产 精品1| 精品久久久噜噜| 国产熟女欧美一区二区| 亚洲av.av天堂| 啦啦啦视频在线资源免费观看| 在线免费观看不下载黄p国产| 中文字幕久久专区| 国产免费一级a男人的天堂| 久久久久久伊人网av| 亚洲国产最新在线播放| 亚洲国产精品国产精品| 女人久久www免费人成看片| 黄片wwwwww| 免费观看无遮挡的男女| 18禁动态无遮挡网站| 国产在线免费精品| 久久精品国产自在天天线| 一级毛片久久久久久久久女| 欧美丝袜亚洲另类| 尾随美女入室| 在线观看免费日韩欧美大片 | 久久精品国产亚洲av天美| 国产午夜精品一二区理论片| av网站免费在线观看视频| 国产男女内射视频| 简卡轻食公司| 精品一区二区三卡| 中文天堂在线官网| 又粗又硬又长又爽又黄的视频| videossex国产| 精品亚洲成国产av| 黄色日韩在线| av免费观看日本| 久久 成人 亚洲| 亚洲精品国产av蜜桃| 一本久久精品| 中文字幕制服av| 国产黄片美女视频| 街头女战士在线观看网站| 视频中文字幕在线观看| 欧美日韩亚洲高清精品| www.av在线官网国产| 91aial.com中文字幕在线观看| 日本猛色少妇xxxxx猛交久久| videossex国产| 一个人免费看片子| 日本黄大片高清| 国产精品一及| 好男人视频免费观看在线| 国产日韩欧美亚洲二区| 亚洲国产高清在线一区二区三| 国产精品一及| 能在线免费看毛片的网站| 国产日韩欧美在线精品| av网站免费在线观看视频| av一本久久久久| 久久99热6这里只有精品| 亚洲欧美日韩东京热| 亚洲美女黄色视频免费看| 亚洲欧美日韩卡通动漫| 蜜桃亚洲精品一区二区三区| 亚洲在久久综合| 最黄视频免费看| 最后的刺客免费高清国语| 新久久久久国产一级毛片| 制服丝袜香蕉在线| 欧美少妇被猛烈插入视频| 性色av一级| 日韩,欧美,国产一区二区三区| 一区在线观看完整版| 好男人视频免费观看在线| 免费少妇av软件| 亚洲激情五月婷婷啪啪| 青青草视频在线视频观看| 91精品国产九色| 黑丝袜美女国产一区| 高清av免费在线| 国产乱来视频区| 国产精品一区二区三区四区免费观看| 久久精品国产自在天天线| 91久久精品国产一区二区成人| 国产久久久一区二区三区| 国产精品一区二区三区四区免费观看| 少妇丰满av| 久久韩国三级中文字幕| 精品亚洲成a人片在线观看 | 久久99热这里只频精品6学生| 国产欧美日韩一区二区三区在线 | 男人爽女人下面视频在线观看| 美女中出高潮动态图| 亚洲真实伦在线观看| 久久综合国产亚洲精品| 亚洲美女搞黄在线观看| 欧美日韩一区二区视频在线观看视频在线| 国产精品一及| 久久久色成人| 赤兔流量卡办理| 少妇人妻一区二区三区视频| 久久 成人 亚洲| 欧美一级a爱片免费观看看| 午夜福利网站1000一区二区三区| 男女边吃奶边做爰视频| 久久99热6这里只有精品| 伦理电影大哥的女人| 国产精品国产三级国产av玫瑰| 99久久精品国产国产毛片| 赤兔流量卡办理| 精品国产乱码久久久久久小说| av国产久精品久网站免费入址| videos熟女内射| 日韩一本色道免费dvd| 日本爱情动作片www.在线观看| 成年av动漫网址| 两个人的视频大全免费| 国产精品爽爽va在线观看网站| 99热全是精品| 国产白丝娇喘喷水9色精品| 国产国拍精品亚洲av在线观看| 水蜜桃什么品种好| 精品少妇黑人巨大在线播放| 最近最新中文字幕免费大全7| 亚洲精品日本国产第一区| 国产伦理片在线播放av一区| 99久久精品热视频| 我要看黄色一级片免费的| 91精品一卡2卡3卡4卡| 九九久久精品国产亚洲av麻豆| 毛片女人毛片| 久久毛片免费看一区二区三区| 亚洲欧美精品专区久久| 午夜老司机福利剧场| 日韩电影二区| 五月伊人婷婷丁香| av国产精品久久久久影院| 亚洲国产毛片av蜜桃av| 亚洲欧美精品专区久久| 最后的刺客免费高清国语| 在线观看一区二区三区激情| av播播在线观看一区| 亚洲综合精品二区| 日日摸夜夜添夜夜添av毛片| 国产成人91sexporn| 最近的中文字幕免费完整| 成人免费观看视频高清| 国产男女内射视频| 欧美极品一区二区三区四区| 高清午夜精品一区二区三区| 成人漫画全彩无遮挡| 伊人久久国产一区二区| 亚洲激情五月婷婷啪啪| 嫩草影院入口| 男女无遮挡免费网站观看| 精品久久久久久久久亚洲| 热99国产精品久久久久久7| 99国产精品免费福利视频| 中文字幕亚洲精品专区| 国产男女超爽视频在线观看| 欧美少妇被猛烈插入视频| 亚洲美女搞黄在线观看| 久久热精品热| 亚洲成人手机| 午夜免费观看性视频| 日韩精品有码人妻一区| 国产成人aa在线观看| 日韩国内少妇激情av| 五月玫瑰六月丁香| 久久久久性生活片| 亚洲色图av天堂| a级毛片免费高清观看在线播放| 久久久久久久大尺度免费视频| 亚洲国产高清在线一区二区三| 我要看日韩黄色一级片| 老女人水多毛片| 亚洲精品日韩av片在线观看| 亚洲av免费高清在线观看| 熟女av电影| 中文精品一卡2卡3卡4更新| 久久久久久久国产电影| 在线观看av片永久免费下载| 亚洲,一卡二卡三卡| 日韩一区二区视频免费看| 亚洲av成人精品一二三区| 99热全是精品| 亚洲第一区二区三区不卡| 黄色怎么调成土黄色| 成人综合一区亚洲| 国产精品秋霞免费鲁丝片| 国产乱来视频区| 寂寞人妻少妇视频99o| 久久国产亚洲av麻豆专区| 人人妻人人看人人澡| 亚洲欧美成人综合另类久久久| 老司机影院成人| 国产成人freesex在线| 国产国拍精品亚洲av在线观看| 免费看av在线观看网站| 亚洲精品aⅴ在线观看| 18禁在线无遮挡免费观看视频| 联通29元200g的流量卡| 菩萨蛮人人尽说江南好唐韦庄| 韩国av在线不卡| 日本免费在线观看一区| 日日摸夜夜添夜夜爱| 欧美日本视频| 日韩中文字幕视频在线看片 | 在线观看免费日韩欧美大片 | 日本欧美视频一区| 免费观看av网站的网址| 青春草亚洲视频在线观看| 精品一区二区三卡| 日韩制服骚丝袜av| 国产精品伦人一区二区| 91午夜精品亚洲一区二区三区| 好男人视频免费观看在线| 一个人看视频在线观看www免费| 欧美精品国产亚洲| 伊人久久国产一区二区| 免费观看a级毛片全部| 51国产日韩欧美| 99精国产麻豆久久婷婷| 啦啦啦中文免费视频观看日本| 国产欧美日韩精品一区二区| 一区二区av电影网| 色视频www国产| 午夜福利网站1000一区二区三区| 国产又色又爽无遮挡免| 亚洲美女视频黄频| 99热6这里只有精品| 国产av精品麻豆| 人人妻人人添人人爽欧美一区卜 | 久久精品人妻少妇| 97热精品久久久久久| 国产毛片在线视频| 97在线视频观看| 夜夜爽夜夜爽视频| 高清日韩中文字幕在线| 欧美丝袜亚洲另类| 国产永久视频网站| 国产精品嫩草影院av在线观看| 亚洲色图av天堂| 91在线精品国自产拍蜜月| 日日摸夜夜添夜夜添av毛片| 夫妻午夜视频| 国产精品麻豆人妻色哟哟久久| 国产一区二区在线观看日韩| 久久精品久久久久久久性| 欧美精品人与动牲交sv欧美| 午夜福利影视在线免费观看| 三级经典国产精品| 国产av国产精品国产| 国产高清三级在线| 日韩精品有码人妻一区| 如何舔出高潮| 精品人妻视频免费看| 99热这里只有是精品在线观看| 制服丝袜香蕉在线| 亚洲va在线va天堂va国产| 久久久久久久久久久免费av| 最近中文字幕高清免费大全6| 亚洲av中文字字幕乱码综合| 一级毛片 在线播放| 直男gayav资源| 少妇人妻 视频| 女的被弄到高潮叫床怎么办| 欧美日韩视频精品一区| 天堂俺去俺来也www色官网| 婷婷色麻豆天堂久久| 亚洲欧美一区二区三区国产| 国产免费一级a男人的天堂| 一区二区三区免费毛片| 成人高潮视频无遮挡免费网站| 激情 狠狠 欧美| 午夜免费观看性视频| 免费久久久久久久精品成人欧美视频 | 欧美性感艳星| 欧美变态另类bdsm刘玥| 午夜激情久久久久久久| 久久国产精品男人的天堂亚洲 | 91在线精品国自产拍蜜月| 99re6热这里在线精品视频| 欧美一区二区亚洲| 免费人成在线观看视频色| 97超视频在线观看视频| 国产乱人偷精品视频| 色婷婷av一区二区三区视频| 久久精品国产鲁丝片午夜精品| 少妇被粗大猛烈的视频| 久久久久久久大尺度免费视频| 久久99蜜桃精品久久| 妹子高潮喷水视频| 免费高清在线观看视频在线观看| 国产精品久久久久久精品电影小说 | 草草在线视频免费看| 国产精品蜜桃在线观看| 久久国产乱子免费精品| 国产精品伦人一区二区| 日本午夜av视频| 中文字幕免费在线视频6| 美女主播在线视频| 国产91av在线免费观看| 草草在线视频免费看| 国产欧美日韩一区二区三区在线 | 超碰97精品在线观看| 国产精品无大码| 国产一级毛片在线| 精华霜和精华液先用哪个| 国产69精品久久久久777片| 亚洲精品日韩在线中文字幕| 久久人人爽人人片av| 午夜福利在线在线| 亚洲四区av| 寂寞人妻少妇视频99o| 日韩亚洲欧美综合| 久久人妻熟女aⅴ| 老司机影院成人| 交换朋友夫妻互换小说| 下体分泌物呈黄色| 亚洲av中文av极速乱| 建设人人有责人人尽责人人享有的 | 午夜福利在线在线| 少妇 在线观看| 国产亚洲5aaaaa淫片| 啦啦啦啦在线视频资源| 3wmmmm亚洲av在线观看| 特大巨黑吊av在线直播| 亚洲熟女精品中文字幕| 国产欧美日韩一区二区三区在线 | 亚洲人成网站高清观看| 日韩大片免费观看网站| 国产一区有黄有色的免费视频| 各种免费的搞黄视频| 国产淫语在线视频| 日韩强制内射视频| 日日摸夜夜添夜夜添av毛片| 最近中文字幕高清免费大全6| 日本黄色日本黄色录像| 毛片一级片免费看久久久久| av福利片在线观看| 性色avwww在线观看| 肉色欧美久久久久久久蜜桃| 黄色视频在线播放观看不卡| 欧美亚洲 丝袜 人妻 在线| 男人狂女人下面高潮的视频| 久久精品久久久久久噜噜老黄| 欧美成人午夜免费资源| 日韩三级伦理在线观看| 日韩欧美一区视频在线观看 | 偷拍熟女少妇极品色| 欧美老熟妇乱子伦牲交| 亚洲色图综合在线观看| 人体艺术视频欧美日本| 精品一区二区三区视频在线| 婷婷色麻豆天堂久久| 国产精品久久久久久久久免| 菩萨蛮人人尽说江南好唐韦庄| 精品亚洲成a人片在线观看 | 国产 一区 欧美 日韩| 久久久久视频综合| 人妻 亚洲 视频| 日本av手机在线免费观看| 亚洲av二区三区四区| 最后的刺客免费高清国语| 久久人妻熟女aⅴ| 狂野欧美白嫩少妇大欣赏| 久久久午夜欧美精品| 久久久久久久大尺度免费视频| 亚洲欧美成人综合另类久久久| 久久久久久久久久久丰满| 欧美精品国产亚洲| 欧美老熟妇乱子伦牲交| 亚洲aⅴ乱码一区二区在线播放| 黄色日韩在线| 春色校园在线视频观看| 男女边摸边吃奶| 免费看日本二区| 日韩不卡一区二区三区视频在线| 精品人妻一区二区三区麻豆| 国产亚洲欧美精品永久| 久久av网站| 欧美xxxx黑人xx丫x性爽| 人体艺术视频欧美日本| 国产精品一区www在线观看| 国产真实伦视频高清在线观看| 午夜精品国产一区二区电影| 免费不卡的大黄色大毛片视频在线观看| 亚洲三级黄色毛片| 99九九线精品视频在线观看视频| 成人高潮视频无遮挡免费网站|