• <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í)別
    一级a做视频免费观看| 久久久久网色| 亚洲精品一区蜜桃| 三级经典国产精品| 草草在线视频免费看| 亚洲精品国产成人久久av| 一级毛片我不卡| 亚洲久久久久久中文字幕| 国产精品人妻久久久久久| 午夜福利网站1000一区二区三区| 亚洲av免费在线观看| 国产免费一区二区三区四区乱码| 国产大屁股一区二区在线视频| 精品酒店卫生间| 只有这里有精品99| 欧美少妇被猛烈插入视频| 国产精品久久久久久av不卡| 欧美极品一区二区三区四区| h日本视频在线播放| 日本色播在线视频| 啦啦啦中文免费视频观看日本| 建设人人有责人人尽责人人享有的 | 亚洲自拍偷在线| 国产亚洲91精品色在线| 男人和女人高潮做爰伦理| 97在线人人人人妻| 亚洲欧美一区二区三区黑人 | 男女啪啪激烈高潮av片| 亚洲av电影在线观看一区二区三区 | 女人久久www免费人成看片| 欧美日韩视频精品一区| 免费观看无遮挡的男女| 涩涩av久久男人的天堂| 亚洲国产高清在线一区二区三| 菩萨蛮人人尽说江南好唐韦庄| 国产成人免费无遮挡视频| 国产精品伦人一区二区| 伊人久久精品亚洲午夜| av播播在线观看一区| 日韩强制内射视频| 国产精品一二三区在线看| 欧美高清成人免费视频www| 亚洲精品乱码久久久久久按摩| 日韩av免费高清视频| 亚洲av免费高清在线观看| 日本一本二区三区精品| 国产免费一级a男人的天堂| 亚洲国产欧美在线一区| 人妻 亚洲 视频| 白带黄色成豆腐渣| 亚洲av日韩在线播放| 国产永久视频网站| 亚洲精品日韩在线中文字幕| 精华霜和精华液先用哪个| 免费不卡的大黄色大毛片视频在线观看| 精品久久国产蜜桃| 王馨瑶露胸无遮挡在线观看| 国产精品99久久99久久久不卡 | 国产精品久久久久久精品古装| 国产综合懂色| 亚洲性久久影院| 天堂中文最新版在线下载 | 好男人视频免费观看在线| 亚洲不卡免费看| h日本视频在线播放| 亚洲国产av新网站| 色婷婷久久久亚洲欧美| 日韩在线高清观看一区二区三区| 午夜福利在线在线| 国产精品爽爽va在线观看网站| 最近2019中文字幕mv第一页| 中国美白少妇内射xxxbb| 自拍偷自拍亚洲精品老妇| 欧美性感艳星| 欧美丝袜亚洲另类| 亚洲国产日韩一区二区| 国产 一区精品| 校园人妻丝袜中文字幕| 亚洲国产成人一精品久久久| 熟女电影av网| 国产亚洲91精品色在线| 啦啦啦啦在线视频资源| 九草在线视频观看| 一级毛片 在线播放| 美女cb高潮喷水在线观看| 国产成人a区在线观看| 永久免费av网站大全| 麻豆国产97在线/欧美| 少妇人妻精品综合一区二区| 久久久久国产网址| 国产精品一及| 狂野欧美激情性bbbbbb| 国产精品人妻久久久久久| 国产精品偷伦视频观看了| 午夜亚洲福利在线播放| 乱码一卡2卡4卡精品| 欧美最新免费一区二区三区| 国产成人免费观看mmmm| 国产精品国产三级国产av玫瑰| 最近中文字幕2019免费版| 99视频精品全部免费 在线| 亚洲精品aⅴ在线观看| 成人国产麻豆网| 日韩强制内射视频| 日本爱情动作片www.在线观看| 成人国产麻豆网| 日本色播在线视频| 我要看日韩黄色一级片| av女优亚洲男人天堂| 国产乱人偷精品视频| 一级毛片黄色毛片免费观看视频| 少妇的逼好多水| 亚洲欧美一区二区三区国产| 九草在线视频观看| 黑人高潮一二区| 国产淫片久久久久久久久| 亚洲人成网站在线观看播放| eeuss影院久久| 国模一区二区三区四区视频| 国产在线男女| 久久精品久久久久久噜噜老黄| 一区二区三区免费毛片| 免费看不卡的av| 国产精品成人在线| 一本久久精品| 日日摸夜夜添夜夜添av毛片| 搡女人真爽免费视频火全软件| 午夜老司机福利剧场| 青春草视频在线免费观看| 男女啪啪激烈高潮av片| 一区二区三区四区激情视频| 欧美一级a爱片免费观看看| 亚洲激情五月婷婷啪啪| 毛片一级片免费看久久久久| 欧美潮喷喷水| 国产一区二区三区av在线| 色播亚洲综合网| 又爽又黄a免费视频| av福利片在线观看| 亚洲精品国产色婷婷电影| 久热久热在线精品观看| 一级毛片黄色毛片免费观看视频| 深爱激情五月婷婷| 中文在线观看免费www的网站| 麻豆久久精品国产亚洲av| 噜噜噜噜噜久久久久久91| 国产69精品久久久久777片| 五月开心婷婷网| 亚洲精品乱久久久久久| 国产69精品久久久久777片| 网址你懂的国产日韩在线| 亚洲精品乱久久久久久| 国产一级毛片在线| 我的老师免费观看完整版| 青春草亚洲视频在线观看| 如何舔出高潮| 国产免费一级a男人的天堂| 成人特级av手机在线观看| 亚洲精品456在线播放app| 人人妻人人澡人人爽人人夜夜| 亚洲自偷自拍三级| 国产精品蜜桃在线观看| 91久久精品电影网| 嫩草影院入口| 美女被艹到高潮喷水动态| 91aial.com中文字幕在线观看| av卡一久久| tube8黄色片| 久久久久久九九精品二区国产| 欧美一区二区亚洲| 51国产日韩欧美| 久热这里只有精品99| 精品人妻偷拍中文字幕| 真实男女啪啪啪动态图| 亚洲在线观看片| 少妇高潮的动态图| 亚洲av国产av综合av卡| 午夜激情久久久久久久| 亚洲欧美清纯卡通| 国产亚洲av片在线观看秒播厂| 一个人看视频在线观看www免费| 91午夜精品亚洲一区二区三区| 国产欧美另类精品又又久久亚洲欧美| 国产精品久久久久久av不卡| 亚洲精品国产av蜜桃| 中国美白少妇内射xxxbb| 国产永久视频网站| 久久人人爽人人爽人人片va| 性色avwww在线观看| 国产免费又黄又爽又色| 97人妻精品一区二区三区麻豆| 日本一二三区视频观看| 91精品一卡2卡3卡4卡| 男女下面进入的视频免费午夜| 日韩亚洲欧美综合| 国产精品.久久久| 啦啦啦在线观看免费高清www| 天堂俺去俺来也www色官网| 九九在线视频观看精品| 日韩成人av中文字幕在线观看| av在线老鸭窝| 日韩不卡一区二区三区视频在线| 三级国产精品片| 在线a可以看的网站| 国产视频首页在线观看| 国产成人精品一,二区| 国产黄频视频在线观看| 成人无遮挡网站| 久久久欧美国产精品| 亚洲欧美成人综合另类久久久| 久久久久精品性色| 精品国产露脸久久av麻豆| 亚洲av电影在线观看一区二区三区 | 久久久久久久久久久免费av| 嫩草影院精品99| 尤物成人国产欧美一区二区三区| 亚洲真实伦在线观看| 18禁裸乳无遮挡免费网站照片| 午夜视频国产福利| 免费不卡的大黄色大毛片视频在线观看| 久久精品久久久久久久性| 亚洲国产高清在线一区二区三| 午夜免费鲁丝| 人人妻人人看人人澡| 亚洲怡红院男人天堂| 亚洲成人久久爱视频| 网址你懂的国产日韩在线| 亚洲精品国产成人久久av| 国产乱人偷精品视频| 国产免费一级a男人的天堂| 99热网站在线观看| 亚洲av二区三区四区| 亚洲天堂av无毛| 日韩强制内射视频| 国产亚洲av嫩草精品影院| 国产成人精品一,二区| 制服丝袜香蕉在线| av又黄又爽大尺度在线免费看| 麻豆成人午夜福利视频| 波多野结衣巨乳人妻| 亚洲av中文av极速乱| 美女xxoo啪啪120秒动态图| 日本黄色片子视频| 亚洲欧美日韩东京热| 男女下面进入的视频免费午夜| 好男人视频免费观看在线| 日本一本二区三区精品| 欧美精品一区二区大全| 久久影院123| 青春草国产在线视频| 亚洲国产日韩一区二区| 26uuu在线亚洲综合色| 下体分泌物呈黄色| 夫妻性生交免费视频一级片| 免费观看无遮挡的男女| 亚洲欧美成人精品一区二区| 亚洲国产欧美在线一区| 久久久a久久爽久久v久久| 美女主播在线视频| 啦啦啦啦在线视频资源| 国产亚洲5aaaaa淫片| 国内精品宾馆在线| 全区人妻精品视频| 美女内射精品一级片tv| 黄色配什么色好看| 亚洲精品国产色婷婷电影| 国产色爽女视频免费观看| 深夜a级毛片| 简卡轻食公司| 免费观看在线日韩| 精品久久久精品久久久| 欧美激情久久久久久爽电影| 一级毛片久久久久久久久女| 国产又色又爽无遮挡免| 精品久久国产蜜桃| 欧美高清性xxxxhd video| 国产黄频视频在线观看| 久久久久性生活片| 久久久亚洲精品成人影院| 高清欧美精品videossex| 色综合色国产| 亚洲色图综合在线观看| 国产精品99久久久久久久久| 听说在线观看完整版免费高清| 久久久a久久爽久久v久久| 男女下面进入的视频免费午夜| 99热这里只有是精品在线观看| 国内精品美女久久久久久| 六月丁香七月| 五月开心婷婷网| 亚洲av一区综合| 美女被艹到高潮喷水动态| 成人国产av品久久久| 寂寞人妻少妇视频99o| 亚洲国产精品成人综合色| 91精品一卡2卡3卡4卡| 国产黄a三级三级三级人| av黄色大香蕉| 免费播放大片免费观看视频在线观看| 女的被弄到高潮叫床怎么办| 国产真实伦视频高清在线观看| 亚洲欧美日韩另类电影网站 | 亚洲熟女精品中文字幕| 国产精品一二三区在线看| 国产永久视频网站| 水蜜桃什么品种好| 午夜福利高清视频| 高清欧美精品videossex| 一个人看的www免费观看视频| 校园人妻丝袜中文字幕| 一个人看视频在线观看www免费| 激情 狠狠 欧美| 欧美潮喷喷水| 男女边吃奶边做爰视频| 久久精品熟女亚洲av麻豆精品| 久久久久久伊人网av| 99热这里只有精品一区| 亚洲欧美成人综合另类久久久| 久热久热在线精品观看| 亚洲天堂av无毛| 大香蕉久久网| 国产精品一区二区在线观看99| 国产男人的电影天堂91| 天堂俺去俺来也www色官网| av天堂中文字幕网| 夫妻午夜视频| 熟妇人妻不卡中文字幕| 身体一侧抽搐| 观看美女的网站| videos熟女内射| 亚洲精品,欧美精品| 亚洲国产最新在线播放| 亚洲无线观看免费| 日本wwww免费看| 国产精品久久久久久精品电影小说 | 搞女人的毛片| 黄色怎么调成土黄色| 尾随美女入室| 偷拍熟女少妇极品色| 午夜视频国产福利| 九色成人免费人妻av| 国产高清三级在线| 国产成人精品婷婷| 国产探花在线观看一区二区| 亚洲av福利一区| 极品少妇高潮喷水抽搐| 97超视频在线观看视频| 日本一本二区三区精品| 国产成人精品婷婷| 人妻 亚洲 视频| 麻豆国产97在线/欧美| 联通29元200g的流量卡| 国产成人福利小说| 另类亚洲欧美激情| 日韩强制内射视频| 国产伦在线观看视频一区| 综合色av麻豆| 日日摸夜夜添夜夜添av毛片| 综合色丁香网| 少妇高潮的动态图| videossex国产| 三级男女做爰猛烈吃奶摸视频| 亚洲高清免费不卡视频| 一级毛片黄色毛片免费观看视频| 亚洲av免费高清在线观看| 午夜福利在线观看免费完整高清在| 日韩免费高清中文字幕av| 熟女电影av网| 在线 av 中文字幕| 91午夜精品亚洲一区二区三区| 亚洲成人久久爱视频| 日韩人妻高清精品专区| 亚洲精品第二区| freevideosex欧美| 国产精品久久久久久精品古装| 亚洲自拍偷在线| 婷婷色av中文字幕| 日韩精品有码人妻一区| 又爽又黄a免费视频| 国产精品一及| 国产色婷婷99| 国产男人的电影天堂91| 如何舔出高潮| 亚洲精品自拍成人| 人体艺术视频欧美日本| 亚洲伊人久久精品综合| 国产精品熟女久久久久浪| 尾随美女入室| 五月玫瑰六月丁香| 欧美 日韩 精品 国产| 国产爱豆传媒在线观看| 91精品伊人久久大香线蕉| 熟女av电影| 久久久精品94久久精品| 男女那种视频在线观看| eeuss影院久久| 99热6这里只有精品| 高清午夜精品一区二区三区| 久久99热这里只频精品6学生| 精品视频人人做人人爽| 亚洲欧美成人综合另类久久久| 亚洲欧美一区二区三区国产| 国产老妇女一区| 少妇裸体淫交视频免费看高清| 亚洲三级黄色毛片| 日日撸夜夜添| 香蕉精品网在线| 婷婷色综合大香蕉| 亚洲天堂国产精品一区在线| 高清av免费在线| 看黄色毛片网站| 久久精品综合一区二区三区| 欧美高清成人免费视频www| 国产亚洲91精品色在线| 欧美高清成人免费视频www| 超碰av人人做人人爽久久| www.av在线官网国产| 亚洲婷婷狠狠爱综合网| 成年女人看的毛片在线观看| 久久精品夜色国产| 五月伊人婷婷丁香| 亚洲天堂国产精品一区在线| 有码 亚洲区| 亚洲欧美成人精品一区二区| 欧美zozozo另类| 欧美最新免费一区二区三区| 国产老妇伦熟女老妇高清| 日韩av不卡免费在线播放| 国产精品一二三区在线看| 亚洲成人中文字幕在线播放| av线在线观看网站| 国产亚洲91精品色在线| 中文资源天堂在线| 国产免费视频播放在线视频| 国产午夜福利久久久久久| 国产成人aa在线观看| av播播在线观看一区| 成年av动漫网址| 丰满人妻一区二区三区视频av| 精品一区二区三卡| 国产精品精品国产色婷婷| 日韩,欧美,国产一区二区三区| 王馨瑶露胸无遮挡在线观看| 边亲边吃奶的免费视频| 成人鲁丝片一二三区免费| 99热国产这里只有精品6| 国产大屁股一区二区在线视频| 国产av不卡久久| 亚洲三级黄色毛片| 在线观看美女被高潮喷水网站| 看十八女毛片水多多多| 大香蕉97超碰在线| 天天躁日日操中文字幕| 久久久久久久久久成人| 成人无遮挡网站| 午夜福利网站1000一区二区三区| 性色av一级| 久久国产乱子免费精品| 在线观看人妻少妇| 国产视频首页在线观看| 国产精品.久久久| 大香蕉97超碰在线| 只有这里有精品99| 精品亚洲乱码少妇综合久久| 少妇的逼好多水| 国产成人一区二区在线| 好男人视频免费观看在线| av在线蜜桃| 国产亚洲一区二区精品| 久久99热这里只有精品18| 青春草视频在线免费观看| 美女xxoo啪啪120秒动态图| 成年版毛片免费区| 亚洲最大成人手机在线| 日本与韩国留学比较| 蜜桃久久精品国产亚洲av| 在线观看三级黄色| 亚洲性久久影院| 国产免费一区二区三区四区乱码| 精品人妻偷拍中文字幕| 亚洲av不卡在线观看| av免费在线看不卡| 97超碰精品成人国产| 国产成人免费无遮挡视频| 制服丝袜香蕉在线| 又爽又黄无遮挡网站| 99久久精品热视频| 18+在线观看网站| 校园人妻丝袜中文字幕| 欧美三级亚洲精品| 色视频www国产| 国产成人91sexporn| 涩涩av久久男人的天堂| 久久久久久久久大av| 少妇人妻精品综合一区二区| 熟妇人妻不卡中文字幕| 亚洲天堂av无毛| 超碰97精品在线观看| 国产亚洲5aaaaa淫片| 亚洲欧美日韩另类电影网站 | 免费人成在线观看视频色| 五月玫瑰六月丁香| 亚洲av成人精品一区久久| 大码成人一级视频| 99热这里只有是精品50| 日韩大片免费观看网站| tube8黄色片| 精品午夜福利在线看| 国产免费又黄又爽又色| 高清av免费在线| 可以在线观看毛片的网站| 国产亚洲91精品色在线| 男女无遮挡免费网站观看| 国产男女内射视频| 久久久精品免费免费高清| 精品一区在线观看国产| 国产精品不卡视频一区二区| 99九九线精品视频在线观看视频| 国产成人freesex在线| 国产亚洲精品久久久com| 蜜臀久久99精品久久宅男| 高清在线视频一区二区三区| 欧美少妇被猛烈插入视频| 看十八女毛片水多多多| 欧美xxxx性猛交bbbb| 又大又黄又爽视频免费| 人妻 亚洲 视频| 天天躁日日操中文字幕| 亚洲精品日韩在线中文字幕| 日日啪夜夜爽| 99九九线精品视频在线观看视频| 国产精品久久久久久久电影| 成人鲁丝片一二三区免费| 亚州av有码| 一个人看视频在线观看www免费| 亚洲四区av| 欧美日韩视频精品一区| 精品一区二区三卡| 一二三四中文在线观看免费高清| 久久热精品热| 亚洲国产精品国产精品| 久久ye,这里只有精品| 三级男女做爰猛烈吃奶摸视频| 久久人人爽人人片av| 亚洲伊人久久精品综合| 国产伦理片在线播放av一区| 高清欧美精品videossex| 深夜a级毛片| 成人综合一区亚洲| 亚洲aⅴ乱码一区二区在线播放| 成人毛片a级毛片在线播放| 熟妇人妻不卡中文字幕| av一本久久久久| 老司机影院毛片| 国产视频内射| 五月玫瑰六月丁香| 国产毛片a区久久久久| 国产精品一区www在线观看| 亚洲美女视频黄频| 在线播放无遮挡| 亚洲综合精品二区| 亚洲精品日韩av片在线观看| 成年女人在线观看亚洲视频 | 纵有疾风起免费观看全集完整版| 99热这里只有精品一区| 日韩强制内射视频| 精品熟女少妇av免费看| 女的被弄到高潮叫床怎么办| 久久精品综合一区二区三区| 国产老妇女一区| 又爽又黄无遮挡网站| 高清av免费在线| 亚洲av二区三区四区| 麻豆精品久久久久久蜜桃| 亚洲精品色激情综合| 九草在线视频观看| 在线免费观看不下载黄p国产| 七月丁香在线播放| 搞女人的毛片| 精品人妻熟女av久视频| 日韩制服骚丝袜av| 高清视频免费观看一区二区| 亚洲精品日本国产第一区| 搡女人真爽免费视频火全软件| 91精品一卡2卡3卡4卡| 日本欧美国产在线视频| 在线观看一区二区三区| 亚洲欧美成人精品一区二区| 日本欧美国产在线视频| 精品国产三级普通话版| 老师上课跳d突然被开到最大视频| 在线观看三级黄色| 久久久久久久亚洲中文字幕| 久久女婷五月综合色啪小说 | 色婷婷久久久亚洲欧美| 日韩人妻高清精品专区| 韩国高清视频一区二区三区| 日韩伦理黄色片| 男女那种视频在线观看| 狂野欧美激情性xxxx在线观看| 亚洲国产日韩一区二区| 一本久久精品| 亚洲精品亚洲一区二区| 中文字幕亚洲精品专区| 日韩大片免费观看网站| 小蜜桃在线观看免费完整版高清| 超碰av人人做人人爽久久| 三级男女做爰猛烈吃奶摸视频| 麻豆成人午夜福利视频| 亚洲精品aⅴ在线观看| 国产精品嫩草影院av在线观看| 日本黄大片高清| 一级黄片播放器| 亚洲真实伦在线观看| 99热全是精品|