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

    增溫與降水變化對青藏高原高寒草甸土壤nirS反硝化菌群落豐度和群落結(jié)構(gòu)的影響

    2018-01-04 05:56:53潘曉悅郭光霞孔維棟
    生態(tài)學(xué)報(bào) 2017年23期
    關(guān)鍵詞:青藏高原硝化群落

    潘曉悅,王 曉,郭光霞,孔維棟

    1 中國礦業(yè)大學(xué)環(huán)境與測繪學(xué)院,徐州 221116 2 中國科學(xué)院青藏高原研究所,高寒生態(tài)學(xué)與生物多樣性重點(diǎn)實(shí)驗(yàn)室,北京 100101 3 中國科學(xué)院南京地理與湖泊研究所,湖泊與環(huán)境國家重點(diǎn)實(shí)驗(yàn)室,南京 210008

    增溫與降水變化對青藏高原高寒草甸土壤nirS反硝化菌群落豐度和群落結(jié)構(gòu)的影響

    潘曉悅1,2,王 曉1,郭光霞2,3,孔維棟2,*

    1 中國礦業(yè)大學(xué)環(huán)境與測繪學(xué)院,徐州 221116 2 中國科學(xué)院青藏高原研究所,高寒生態(tài)學(xué)與生物多樣性重點(diǎn)實(shí)驗(yàn)室,北京 100101 3 中國科學(xué)院南京地理與湖泊研究所,湖泊與環(huán)境國家重點(diǎn)實(shí)驗(yàn)室,南京 210008

    全球變化已成為國際研究熱點(diǎn)。青藏高原屬典型生態(tài)脆弱帶,該地區(qū)升溫幅度更加明顯,已導(dǎo)致大量冰川融化和明顯降水變化,進(jìn)而使該地區(qū)水循環(huán)和土壤水分發(fā)生巨大變化。溫度和降水的變化可能會引起土壤微生物豐度和群落結(jié)構(gòu)的改變,進(jìn)而影響生物地球化學(xué)循環(huán)。但青藏高原地區(qū)土壤微生物群落結(jié)構(gòu)和功能對全球變化響應(yīng)的研究較少。研究了模擬增溫和降水變化對青藏高原高寒草甸土壤nirS反硝化菌群落豐度和群落結(jié)構(gòu)的影響。研究表明,增溫1、2、4℃對nirS基因豐度影響不顯著;增加降水100%時(shí),增溫4℃處理顯著增加nirS基因豐度(P<0.05)。在未升溫與升溫2℃背景下增加和減少降水對nirS基因豐度的影響不顯著。增溫和增減降水均顯著影響nirS反硝化菌群落結(jié)構(gòu),且兩個因子具有一定的交互作用。CCA結(jié)果顯示,增溫和降水的共同解釋變量中,增溫對nirS反硝化菌群落結(jié)構(gòu)變化的影響達(dá)極顯著(P<0.01),解釋了其中的54.2%,降水變化解釋了45.5%(P<0.05)。

    nirS;基因豐度;群落結(jié)構(gòu);氣候變化;青藏高原

    氣候變化已成為典型的全球尺度的環(huán)境問題。根據(jù)政府間氣候變化專門委員會(IPCC)第5次評估報(bào)告,1901—2012年間,全球平均表面溫度上升0.89℃,報(bào)告預(yù)測到2100年全球氣溫將升高2—4.8℃,海平面上升26—81cm,冰川和冰原也以更快的速度融化,到2050年北極地區(qū)將成為幾近無冰區(qū)[1]。青藏高原地區(qū)平均氣溫在1961—2007年間每10年上升0.36℃[2],約為全球氣溫增幅的3倍(每10年上升溫度0.12℃,約為1951—2012年間)[1]。氣候變暖已經(jīng)導(dǎo)致青藏高原地區(qū)冰川融化速度加快,且該地區(qū)降水具有增加的趨勢[3],這些導(dǎo)致該地區(qū)水循環(huán)和土壤水分顯著變化[4-5]。

    目前,青藏高原地區(qū)功能微生物群落的研究還處于初步階段,對升溫和降水變化的響應(yīng)研究還比較少。本文以青藏高原那曲地區(qū)高寒草甸土壤為研究對象,采用微生物分子生態(tài)學(xué)研究方法研究了土壤增溫和降水變化下,青藏高原高寒草甸土壤中nirS反硝化菌群落豐度和群落結(jié)構(gòu)的變化,為研究高寒生態(tài)系統(tǒng)反硝化微生物對全球變暖的響應(yīng)以及全球變暖情況下N2O排放提供理論依據(jù)。

    1 材料與研究方法

    1.1 試驗(yàn)地點(diǎn)

    增溫試驗(yàn)平臺選擇在青藏高原那曲地區(qū)現(xiàn)代化草地畜牧業(yè)科技示范園區(qū)內(nèi)的典型高山嵩草草甸。那曲地區(qū)位于西藏自治區(qū)北部(29°55′—36°30′N,83°55′—95°05′E)。平均海拔4500m左右,多雪峰和高山,東部是高山峽谷,中西部為高原湖盆。該地區(qū)屬于青藏高寒氣候區(qū),干燥寒冷,年平均氣溫-2.9—3.4℃,最冷月為1月,月平均氣溫-14.9—7.4℃;最熱月為7月,月平均氣溫8.7—12.2℃,年平均降水量在298.6—708.4mm,年內(nèi)降雨量的80%以上都集中在5—9月[17]。那曲地區(qū)植被類型從東南向西北依次為高山疏林灌叢草甸、高寒草甸、高寒草原、高寒荒漠草原四大基帶及相鄰的過度亞帶,分布與生物氣候帶相適應(yīng)[18],該地區(qū)主要草地土壤類型有高山草甸土、亞高山草甸土、高山草原土、高山荒漠草原土等[19]。

    1.2 試驗(yàn)設(shè)計(jì)

    實(shí)驗(yàn)平臺于2014年7月份開始運(yùn)行。利用紅外照射的自動控制系統(tǒng)準(zhǔn)確控制增溫,設(shè)置4個增溫梯度:0(對照)、1.0、2.0、4.0℃,每個增溫實(shí)驗(yàn)小區(qū)為3m直徑的圓形,利用“巢式”設(shè)計(jì),在每一種增溫幅度下,利用PVC板(33cm深)將該小區(qū)隔成2個分區(qū),一個分區(qū)不增加降水(對照),另一個分區(qū)增加降水100%,共8個處理。針對學(xué)者提出的“2℃閾值”,在增溫0、2℃條件下,設(shè)置減水50%的試驗(yàn),減水小區(qū)占整個直徑為3m的圓形。每個減水小區(qū),10個透明的聚碳酸樹脂通道(投影面積占整個小區(qū)面積的50%)以15°角傾斜固定在距地面160cm處,截留的雨水通過管道流入白色聚乙烯塑料雨水采集器。每次降雨,將減水小區(qū)采集器的雨水均勻添加到增水小區(qū),達(dá)到增水100%。每個處理4個重復(fù)。在另外地方單獨(dú)放置白色聚乙烯塑料雨水采集器,采集降水,用于增水小區(qū)的增水。設(shè)有溫度探針,監(jiān)控溫度變化。

    1.3 土壤采集

    為研究nirS反硝化菌對短期增溫和降水變化的響應(yīng)。采樣時(shí)間為2014年10月份。在每個實(shí)驗(yàn)小區(qū)內(nèi),隨機(jī)選擇采樣點(diǎn),使用土鉆采集5cm深的土壤樣品,采集后去除植物根系等雜質(zhì),過篩(2mm)混勻,用于分子生物學(xué)分析的樣品放入干冰盒中保存,運(yùn)回拉薩實(shí)驗(yàn)室后轉(zhuǎn)入-80℃冰箱保存,之后放于冰盒中空運(yùn)至北京實(shí)驗(yàn)室。

    1.4 土壤DNA的提取

    DNA提取采用Powersoil? DNA Isolation Kit(MOBIO,USA),使用0.5g土壤,方法完全參照制造商說明書。提取結(jié)束后取1 μL的DNA溶液于NanoDrop2000核酸定量儀(Thermo Scientific)中測定其DNA含量,使用前貯存于-20℃。

    1.5 qPCR

    nirS基因豐度采用LightCycler 480 (Roche)實(shí)時(shí)熒光定量PCR儀進(jìn)行定量試驗(yàn),擴(kuò)增引物[20]為nirS-cd3aF(5′-GTSAACGTSAAGGARACSGG),nirS-R3cd(5′-GASTTCGGRTGSGTCTTGA)(由生工生物工程(上海)股份有限公司合成),酶使用SYBR green kit(TaKaRa),反應(yīng)體系10μL,包括:3μL緩沖液,5μL的酶,引物各0.5μL,1μL DNA模板。反應(yīng)程序?yàn)椋?4℃預(yù)變性2min,94℃變性30s,58℃復(fù)性1min,72℃延伸30s,進(jìn)行10個循環(huán),從第二次循環(huán)開始復(fù)性溫度依次降0.5℃,94℃變性30s,53℃復(fù)性1min,72℃延伸30s,進(jìn)行25個循環(huán),最后72℃延伸10min。本次定量實(shí)驗(yàn)結(jié)果擴(kuò)增效率>99%,標(biāo)準(zhǔn)曲線斜率約為-3.3。

    1.6 末端限制性片段長度多態(tài)性分析(T-RFLP)

    nirS基因PCR擴(kuò)增引物序列與定量PCR相同,引物nirS-cd3aF5′端用FAM標(biāo)記。擴(kuò)增儀器使用ABI Veriti- 96梯度PCR儀(Applied Biosystems)。PCR反應(yīng)體系25μL,包括:9.5μL緩沖液,12.5μL的2×EasyTaq PCR SuperMix(全式金),引物各1μL,1 μL DNA模板。反應(yīng)程序與定量PCR相同。由于末端限制性片段多態(tài)性分析要求PCR純化產(chǎn)物濃度較高,因此采用切膠純化的方法去除非特異性帶的影響,膠回收使用AxyPrepTMDNA Gel Extraction Kit(Axygen),然后,采用內(nèi)切酶HhaI(TaKaRa)對純化產(chǎn)物進(jìn)行酶切,PCR純化產(chǎn)物及酶混勻離心后在37℃條件下消化3h,而后65℃變性20min。純化后的酶切產(chǎn)物送至中國農(nóng)業(yè)科學(xué)院進(jìn)行基因掃描。不加水背景下增溫0、1℃的樣品,有一重復(fù)沒有收集到足量的nirS基因,只有3個重復(fù)。

    1.7 數(shù)據(jù)分析

    定量PCR數(shù)據(jù)分析采用 SPSS Statistics 23軟件中的one-way ANOVA 及LSD方法比較不同處理間的差異性,并用Excel 2010作圖。群落結(jié)構(gòu)采用T-RFLP和Canoco 5軟件中的主成分(principal component analysis, PCA)、典型相關(guān)分析(canonical correlation analysis, CCA)方法,采用PCA分析群落結(jié)構(gòu)受短期增溫或降水變化單一因素影響的分布規(guī)律,采用CCA分析短期增溫和降水變化對群落結(jié)構(gòu)綜合影響的顯著性和貢獻(xiàn)率。

    2 結(jié)果與分析

    2.1 短期增溫對nirS基因豐度的影響

    從圖1可以看出,在不加水情況下,與不升溫對照處理相比,增溫1、2、4℃對nirS基因豐度影響不大。增加降水100%情況下,增溫4℃顯著提高nirS基因豐度(P<0.05),其基因拷貝數(shù)約為其他3個升溫處理的2倍,增溫1、2℃與對照相比差異不顯著。這說明增加降水在一定程度上改變了nirS基因豐度對升溫的響應(yīng)。不同溫度梯度下,增加降水100%與不加水處理相比,基因豐度進(jìn)行顯著性差異檢驗(yàn)結(jié)果都沒有顯著差異。

    2.2 短期降水變化對nirS基因豐度的影響

    由圖2可知,不增溫時(shí),nirS基因豐度隨降水的增加有逐漸上升的趨勢,增水100%處理的基因豐度顯著高于減水50%處理(P<0.05)。增溫2℃情況下,減水50%、對照、加水100%3個處理的基因豐度差別不大。土壤在不同的降水量情況下,增溫2℃與不增溫處理之間,差異不明顯。

    圖1 nirS基因豐度隨土壤增溫的變化Fig.1 nirS gene abundance with soil temperature change

    圖2 nirS基因豐度隨降水的變化Fig.2 nirS gene abundance with the precipitation change

    2.3 短期增溫對nirS反硝化菌群落結(jié)構(gòu)的影響

    增溫對nirS反硝化菌群落結(jié)構(gòu)的影響如圖3,在土壤不加水背景下,群落結(jié)構(gòu)隨溫度的變化和在土壤增加降水100%背景下,群落結(jié)構(gòu)隨溫度的變化。

    由圖3的PCA結(jié)果可知,主成分1和主成分2(PCA1、PCA2)可解釋nirS反硝化菌群落結(jié)構(gòu)總變異的43.67%。增溫1、2、4℃的3種處理均改變了nirS反硝化菌的群落結(jié)構(gòu);增溫1℃處理的群落結(jié)構(gòu)變化與增溫2℃和增溫4℃有差異,增溫4℃與增溫2℃相比,在圖中的位置區(qū)分并不明顯。

    從圖3可以看出,PCA1、PCA2共解釋了nirS反硝化菌群落結(jié)構(gòu)總變異的63.81%。土壤增加降水100%背景下,增溫對nirS反硝化菌群的影響與降水不變背景不同。增水背景下,隨著增溫程度的增加,群落結(jié)構(gòu)持續(xù)發(fā)生變化。其中增溫1℃到增溫2℃之間差異最明顯;與對照位置相比,增溫2℃與增溫4℃的距離比增溫1℃更加明顯,而增溫2℃與增溫4℃之間差異較小。降水變化,改變了nirS反硝化菌群落結(jié)構(gòu)隨溫度變化的規(guī)律。

    圖3 增溫對nirS群落結(jié)構(gòu)影響的主成分分析Fig.3 Principal component analysis of T-RFLP profiles for nirS community structure in response to temperature increase T0, 增溫0℃increased 0℃; T1, 增溫1℃ increased 1℃; T2, 增溫2℃ increased 2℃; T4, 增溫4℃ increased 4℃; T0+W, 增溫0℃加水100% increased 0℃ and double precipitation; T1+W, 增溫1℃加水100% increased 1℃ and double precipitation; T2+W, 增溫2℃加水100% increased 2℃ and double precipitation; T4+W, 增溫4℃加水100% increased 4℃ and double precipitation

    2.4 短期降水變化對nirS反硝化菌群落結(jié)構(gòu)的影響

    降水變化對nirS反硝化菌群落結(jié)構(gòu)的影響如圖4,土壤不增溫背景下,群落結(jié)構(gòu)隨土壤增減水的變化,土壤增溫2℃背景下,群落結(jié)構(gòu)隨土壤增減水的變化。

    結(jié)果表明,土壤不增溫背景下,PCA1和PCA2共解釋nirS反硝化菌群落結(jié)構(gòu)總變異的65.04%。其中土壤減水50%在第一主軸的右側(cè)和第二主軸的下端,而不加水和增加降水100%在第一主軸的左側(cè)和第二主軸的上端。與對照相比,土壤減水50%時(shí)的群落結(jié)構(gòu)比增加降水100%變化更明顯,增加降水100%與對照相比群落結(jié)構(gòu)基本無差異。

    結(jié)果說明,土壤增溫2℃背景下,PCA1和PCA2共解釋了nirS反硝化菌群落結(jié)構(gòu)總變異的58.32%。減水50%時(shí)的群落結(jié)構(gòu)與對照相比產(chǎn)生了較明顯的變異,而增加降水100%與對照相比群落結(jié)構(gòu)基本無差異。減水50%分布在第一主軸的左側(cè)和第二主軸的上端,增加降水100%與不加水分布在第一主軸的右側(cè)和第二主軸的下端。

    由圖4可得,土壤增溫2℃和不增溫相比,降水變化時(shí),nirS反硝化菌群落結(jié)構(gòu)的在圖中的位置和變化方向發(fā)生改變,說明土壤溫度增加影響了降水變化對群落結(jié)構(gòu)影響的規(guī)律。

    圖4 降水變化對nirS群落結(jié)構(gòu)影響的主成分分析Fig.4 Principal component analysis of T-RFLP profiles for nirS community structure in response to precipitation changeT0, 增溫0℃ increased 0℃;T0-W, 增溫0℃減水50% increased 0℃ and half of the precipitation;T0+W, 增溫0℃加水100% increased 0℃ and double precipitation;T2, 增溫2℃ increased 2℃;T2-W, 增溫2℃減水50% increased 2℃ and half of the precipitation;T2+W, 增溫2℃加水100% increased 2℃ and double precipitation

    2.5 短期增溫和降水變化對nirS反硝化菌群落結(jié)構(gòu)的綜合影響

    由圖5,PCA和CCA結(jié)果可得,nirS反硝化菌群落結(jié)構(gòu)的變化是由土壤溫度和降水聯(lián)合影響產(chǎn)生。CCA分析結(jié)果顯示增溫和降水變化對整個群落結(jié)構(gòu)的變化共解釋了19.9%,溫度解釋了其中的54.2%,對nirS反硝化菌群落結(jié)構(gòu)的變化有極顯著影響(P<0.01);水分貢獻(xiàn)了其中的45.5%,對nirS反硝化菌群落結(jié)構(gòu)的變化有顯著影響(P<0.05)。增溫和降水兩個環(huán)境因子交互作用對群落結(jié)構(gòu)的影響的顯著性檢驗(yàn)結(jié)果顯示,對nirS反硝化菌影響不顯著。

    圖5 增溫和降水變化對nirS群落結(jié)構(gòu)影響的典范對應(yīng)分析Fig.5 Canonical correspondence analysis of T-RFLP profiles for nirS community structure in response to precipitation change and temperature increaseT0, 增溫0℃ increased 0℃;T0-W, 增溫0℃減水50% increased 0℃ and half of the precipitation;T0+W, 增溫0℃加水100% increased 0℃ and double precipitation;T2, 增溫2℃ increased 2℃;T2-W, 增溫2℃減水50% increased 2℃ and half of the precipitation;T2+W, 增溫2℃加水100% increased 2℃ and double precipitation

    3 討論

    本研究結(jié)果表明nirS基因豐度對短期增溫的響應(yīng)不敏感,只在增加降水100%背景下,增溫4℃,nirS基因豐度才顯著增加。nirS基因豐度對增溫響應(yīng)不敏感已有研究證明,如Jung等[21]研究表明,南極土壤增溫30d處理,nirS基因豐度變化不明顯,Zheng等[22]研究的青藏高原高寒草甸土壤結(jié)果也表明溫度的變化對nirS反硝化菌群落豐度的影響不顯著。也有可能的原因是,增溫處理時(shí)間較短,如Penton等[23]的研究結(jié)果顯示,在永久凍土增溫(1.5℃)1年后nirS基因豐度沒有明顯的改變,而高原土壤增溫(1.8—2.7℃)10年后,nirS豐度有明顯的增加。本研究出現(xiàn)的在增加降水100%和增溫4℃時(shí),nirS豐度有明顯的增加,其中的原因可能為:nirS基因微生物對于增溫2℃以內(nèi)有一定的適應(yīng)能力,如Zhang等[24]的研究,在增溫約1.79℃,持續(xù)5年后nirS基因豐度變化不明顯,而主要受到N的增加量和水分增加的影響,同時(shí),在降水減少或不變情況下,增溫后土壤濕度降低[25],不適宜nirS反硝化菌生存。

    短期降水變化對nirS基因豐度的影響規(guī)律不明顯,只有在不增溫背景下,減水50%和加水100%之間nirS基因豐度有明顯的區(qū)別。這與之前的一些研究結(jié)果不同,Ding等[26]的研究結(jié)果表明,在不同降水區(qū)域,nirS基因豐度隨區(qū)域降水量的增加而明顯增加,Zhang等[24]的研究結(jié)果也顯示,經(jīng)過5a的增水處理,nirS基因豐度有明顯的增加。出現(xiàn)不同結(jié)果的原因可能是:本次實(shí)驗(yàn)降水變化處理的時(shí)間較短,短期內(nèi)nirS反硝化菌有一定的適應(yīng)能力;研究區(qū)域土壤性質(zhì)和環(huán)境的不同,反硝化菌還受到空間位置、有機(jī)碳、N的化合物的影響[27- 29],進(jìn)而導(dǎo)致nirS微生物響應(yīng)的結(jié)果不同,上述兩個研究土壤樣品均采于內(nèi)蒙古草原地區(qū)。本研究結(jié)果顯示不同增溫幅度下,nirS基因豐度隨降水變化規(guī)律的不同,Szukics等[30]的研究,nirK基因豐度在不同的溫度處理下,隨水分變化規(guī)律也不一致,可能因?yàn)樵鰷厥雇寥罎穸扔幸欢ǔ潭鹊慕档蚚25],導(dǎo)致增溫后變化規(guī)律受到影響。

    與基因豐度不同,土壤短期增溫顯著影響青藏高原高寒草甸地區(qū)nirS反硝化菌的群落結(jié)構(gòu)。這與之前的研究結(jié)果相同,加拿大農(nóng)業(yè)土壤在-4℃,-1℃,+2℃,+5℃溫度梯度下,經(jīng)過120d后,不同的溫度處理之間nirS反硝化菌的群落結(jié)構(gòu)有明顯的區(qū)別[31]。Braker等[32]對nirS反硝化菌增溫培養(yǎng)3周后,群落多樣性和群落結(jié)構(gòu)都發(fā)生了明顯的變化。增溫對參與N循環(huán)的土壤微生物的群落結(jié)構(gòu)產(chǎn)生影響已經(jīng)被很多研究證實(shí)。Szukics等[30]研究結(jié)果顯示,參與N循環(huán)的基因(nirK、amoA)群落結(jié)構(gòu)在增溫培養(yǎng)1周之后,產(chǎn)生了明顯的改變。Yergeau等[33]利用16S rRNA基因研究南極地區(qū)微生物對增溫的響應(yīng),增溫約3年結(jié)果顯示,增溫對功能菌尤其是參與N循環(huán)微生物有明顯的影響。

    短期降水變化對青藏高原高寒草甸地區(qū)nirS反硝化菌的群落結(jié)構(gòu)產(chǎn)生了顯著的影響。這一結(jié)果與Jha等[34]對牧場土壤的研究結(jié)果相同,nirS反硝化菌與土壤含水量有顯著相關(guān)性,Hamonts等[35]對根際土壤的研究結(jié)果也表明,土壤水分明顯增加時(shí),nirS反硝化菌群落結(jié)構(gòu)也發(fā)生明顯的改變。但是本研究從PCA結(jié)果看,不增溫和增溫2℃的結(jié)果相似,均是增加降水100%和對照相比群落結(jié)構(gòu)差異不明顯,減水50%分布與前兩者有較大的區(qū)別。可能原因有:1)土壤增水造成了土壤的厭氧環(huán)境[32],反硝化菌更適應(yīng)此環(huán)境[36-37],所以經(jīng)過3個月增水100%處理沒有發(fā)現(xiàn)較明顯的改變,而減水后的環(huán)境更不利于反硝化菌的生存,進(jìn)而導(dǎo)致明顯的變化;2)降水變化改變了土壤中N的無機(jī)化合物的數(shù)量和存在形式[24,29],本研究此假設(shè)還需進(jìn)一步研究考證。

    短期增溫和降水變化對nirS反硝化菌群落結(jié)構(gòu)有顯著影響,并且兩因子之間具有一定的相互作用。其原因可能為溫度和降水的變化引起植物和土壤生物化學(xué)性質(zhì)的改變,進(jìn)而引起土壤微生物群落結(jié)構(gòu)的變異[38-39],所以不同的環(huán)境和試驗(yàn)地條件下,實(shí)驗(yàn)結(jié)果和結(jié)論會有一定的不同,與不同地點(diǎn)的植物、土壤特性、土壤營養(yǎng)元素和微生物特征等都有很大的關(guān)系。但是兩個因素對nirS基因豐度影響并不明顯,只有在降水和溫度適宜的條件下基因豐度才會有明顯的增加,這可能說明受到短期的增溫和降水變化的影響,nirS反硝化菌的優(yōu)勢菌種發(fā)生了改變[23],但所含nirS基因的量沒有明顯的變化。反硝化過程由narG或napA、nirK或nirS、norB或norZ和nosZ等多種基因編碼的酶共同完成,所以為研究溫度和降水變化對反硝化過程的影響,還需研究環(huán)境因子對其他反硝化基因的作用,以及進(jìn)一步研究長期增溫和降水變化對反硝化功能微生物的影響,進(jìn)而研究全球變暖對整個反硝化過程及N2O的排放的影響。

    4 結(jié)論

    本研究結(jié)果表明,土壤短期增溫和降水變化顯著影響青藏高原高寒草甸土壤nirS反硝化菌群落結(jié)構(gòu),降水變化調(diào)控增溫對微生物群落結(jié)構(gòu)的影響。短期增溫和降水變化對青藏高原高寒草甸土壤nirS反硝化菌群落豐度的影響較小,僅在增水100%背景下增溫4℃豐度顯著增加。

    [1] IPCC. Climate Change 2013: The Physical Science Basis. Cambridge: Cambridge University Press, 2013.

    [2] Wang B, Bao Q, Hoskins B, Wu G X, Liu Y M. Tibetan plateau warming and precipitation changes in East Asia. Geophysical Research Letters,2008,35(14): L14702.

    [3] Lutz A F, Immerzeel W W, Shrestha A B, Bierkens M F P. Consistent increase in High Asia′s runoff due to increasing glacier melt and precipitation. Nature Climate Change, 2014, 4(7): 587- 592.

    [4] Yang W, Guo X F, Yao T D, Zhu M L, Wang Y J. Recent accelerating mass loss of southeast Tibetan glaciers and the relationship with changes in macroscale atmospheric circulations.Climate Dynamics, 2016, 47(3- 4): 805- 815.

    [5] Yang W, YaoT D, Xu B Q, Ma L L, Wang Z H, Wan M. Characteristics of recent temperate glacier fluctuations in the Parlung Zangbo River basin, southeast Tibetan Plateau. Chinese Science Bulletin, 2010, 55(20): 2097- 2102.

    [6] 單文俊, 王慶貴, 閆國永, 邢亞娟. 基于土壤微生物的碳氮互作效應(yīng)綜述. 中國農(nóng)學(xué)通報(bào), 2016, 32(23): 65- 71.

    [7] 沈菊培, 賀紀(jì)正. 微生物介導(dǎo)的碳氮循環(huán)過程對全球氣候變化的響應(yīng). 生態(tài)學(xué)報(bào), 2011, 31(11): 2957- 2967.

    [8] 吳建國, 呂佳佳. 氣候變化對青藏高原高寒草甸適宜氣候分布范圍的潛在影響. 草地學(xué)報(bào), 2009, 17(6): 699- 705.

    [9] 蘆曉飛. 西藏米拉山高寒草甸土壤微生物多樣性研究[D]. 北京: 中國農(nóng)業(yè)科學(xué)院, 2009.

    [10] IPCC. Climate Change 2007: The Physical Science Basis. IPCC WGI Fourth Assessment Report, Geneva: Intergovernmental Panel on Climate Change, 2007.

    [11] 付曉青, 李勇. 土壤氧化亞氮排放時(shí)空變異性及其方法研究進(jìn)展. 生態(tài)學(xué)雜志, 2012, 31(3): 724- 730.

    [12] Cicerone R J. Analysis of sources and sinks of atmospheric nitrous oxide (N2O). Journal of Geophysical Research Atmospheres, 1989, 941(D15): 18265- 18271.

    [13] 張玉銘, 胡春勝, 張佳寶, 董文旭, 王玉英, 宋利娜.農(nóng)田土壤主要溫室氣體(CO2、CH4、N2O)的源/匯強(qiáng)度及其溫室效應(yīng)研究進(jìn)展. 中國生態(tài)農(nóng)業(yè)學(xué)報(bào), 2011, 19(4): 966- 975.

    [14] 朱永官, 王曉輝, 楊小茹, 徐會娟, 賈炎. 農(nóng)田土壤N2O產(chǎn)生的關(guān)鍵微生物過程及減排措施. 環(huán)境科學(xué), 2014, 35(2): 792- 800.

    [15] Zhu X, Burger M, Doane T A, Horwath W R. Ammonia oxidation pathways and nitrifier denitrification are significant sources of N2O and NO under low oxygen availability. Proceedings of the National Academy of Sciences of the United States of America, 2013, 110(16): 6328- 6333.

    [16] 王海濤, 鄭天凌, 楊小茹. 土壤反硝化的分子生態(tài)學(xué)研究進(jìn)展及其影響因素. 農(nóng)業(yè)環(huán)境科學(xué)學(xué)報(bào), 2013, 32(10): 1915- 1924.

    [17] 張佳華. 中國藏北地區(qū)生態(tài)環(huán)境與氣象災(zāi)害遙感. 北京: 氣象出版社, 2007: 17- 17.

    [18] 西藏自治區(qū)土地管理局、畜牧局. 西藏自治區(qū)草地資源. 北京: 科學(xué)出版社, 1994: 371- 372.

    [19] 西藏自治區(qū)土地管理局. 西藏自治區(qū)土壤資源. 北京: 科學(xué)出版社, 1994: 683- 697.

    [20] Throb?ck I N, Enwall K, Jarvis ?, Hallin S. Reassessing PCR primers targetingnirS,nirKandnosZgenes for community surveys of denitrifying bacteria with DGGE. FEMS Microbiology Ecology, 2004, 49(3): 401- 417.

    [21] Jung J, Yeom J, Kim J, Han J, Lim H S, Park H, Hyun S, Park W. Change in gene abundance in the nitrogen biogeochemical cycle with temperature and nitrogen addition in Antarctic soils. Research in Microbiology, 2011, 162(10): 1018- 1026.

    [22] Zheng Y, Yang W, Hu H W, Kim Y C, Duan J C, Luo C Y, Wang S P, Guo L D. Ammonia oxidizers and denitrifiers in response to reciprocal elevation translocation in an alpine meadow on the Tibetan Plateau. Journal of Soils and Sediments, 2014, 14(6): 1189- 1199.

    [23] Penton C R, Louis D S, Pham A, Cole J R, Wu L Y, Luo Y Q, Schuur E A G, Zhou J Z, Tiedje J M. Denitrifying and diazotrophic community responses to artificial warming in permafrost and tallgrass prairie soils. Frontiers in Microbiology, 2015, 6: 746.

    [24] Zhang X M, Liu W, Schloter M, Zhang G M, Chen Q S, Huang J H, Li L H, Elser J J, Han X G. Response of the abundance of key soil microbial nitrogen-cycling genes to multi-factorial global changes. PLoS One, 2013, 8(10): e76500.

    [25] 武丹丹, 井新, 林笠, 楊新宇, 張振華, 賀金生. 青藏高原高寒草甸土壤無機(jī)氮對增溫和降水改變的響應(yīng). 北京大學(xué)學(xué)報(bào): 自然科學(xué)版, 2016, 52(5): 959- 966.

    [26] Ding K, Zhong L, Xin X P, Xu Z H, Kang X M, Liu W J, Rui Y C, Jiang L L, Tang L, Wang Y F. Effect of grazing on the abundance of functional genes associated with N cycling in three types of grassland in Inner Mongolia. Journal of Soils and Sediments, 2015, 15(3): 683- 693.

    [27] Dandie C E, Burton D L, Zebarth B J, Henderson S L, Trevors J T, Goyer C. Changes in bacterial denitrifier community abundance over time in an agricultural field and their relationship with denitrification activity. Applied and Environmental Microbiology, 2008, 74(19): 5997- 6005.

    [28] Clark I M, Buchkina N, Jhurreea D, Goulding K W T, Hirsch P R. Impacts of nitrogen application rates on the activity and diversity of denitrifying bacteria in the Broadbalk Wheat Experiment. Philosophical Transactions of the Royal Society B: Biological Sciences, 2012, 367(1593): 1235- 1244.

    [29] Pajares S, Bohannan B J M. Ecology of nitrogen fixing, nitrifying, and denitrifying microorganisms in tropical forest soils. Frontiers in Microbiology, 2016, 7: 1045.

    [30] Szukics U, Abell G C, H?dl V, Mitter B, Sessitsch A, Hackl E, Zechmeister-Boltenstern S. Nitrifiers and denitrifiers respond rapidly to changed moisture and increasing temperature in a pristine forest soil. FEMS Microbiology Ecology, 2010, 72(3): 395- 406.

    [31] Wertz S, Goyer C, Zebarth B J, Burton D L, Tatti E, Chantigny M H, Filion M. Effects of temperatures near the freezing point on N2O emissions, denitrification and on the abundance and structure of nitrifying and denitrifying soil communities. FEMS Microbiology Ecology, 2013, 83(1): 242- 254.

    [32] Braker G, Schwarz J, Conrad R. Influence of temperature on the composition and activity of denitrifying soil communities. FEMS Microbiology Ecology, 2010, 73(1): 134- 148.

    [33] Yergeau E, Bokhorst S, Kang S, Zhou J Z, Greer C W, Aerts R, Kowalchuk G A. Shifts in soil microorganisms in response to warming are consistent across a range of Antarctic environments.The ISME Journal, 2012, 6(3): 692- 702.

    [34] Jha N, SaggarS, Giltrap D, Tillman R, Deslippe J R. Soil properties impacting denitrifier community size, structure, and activity in New Zealand dairy-grazed pasture. Biogeosciences Discuss,2016, doi:10.5194/bg- 2016- 390.

    [35] Hamonts K, Clough T J, Stewart A, Clinton P W, Richardson A E, Wakelin S A, O′Callaghan M, Condron L M. Effect of nitrogen and waterlogging on denitrifier gene abundance, community structure and activity in the rhizosphere of wheat. FEMS Microbiology Ecology, 2013, 83(3): 568- 584.

    [36] Bateman E J, Baggs E M. Contributions of nitrification and denitrification to N2O emissions from soils at different water-filled pore space. Biology and Fertility of Soils, 2005, 41(6): 379- 388.

    [37] Diba F, Shimizu M, Hatano R. Effects of soil aggregate size, moisture content and fertilizer management on nitrous oxide production in a volcanic ash soil. Soil Science and Plant Nutrition, 2011, 57(5): 733- 747.

    [38] Wang Z, Luo T X, Li R C, Tang Y H, Du M Y. Causes for the unimodal pattern of biomass and productivity in alpine grasslands along a large altitudinal gradient in semi-arid regions.Journal of Vegetation Science, 2013, 24(1): 189- 201.

    [39] Guo G X, Kong W D, Liu J B, Zhao J X, Du H D, Zhang X Z, Xia P H. Diversity and distribution of autotrophic microbial community along environmental gradients in grassland soils on the Tibetan Plateau. Applied Microbiology and Biotechnology, 2015, 99(20): 8765- 8776.

    EffectsofincreasedtemperatureandprecipitationchangeonnirSgeneabundanceandcommunitystructureinalpinemeadowsoilsontheQinghai-TibetPlateau

    PAN Xiaoyue1,2,WANG Xiao1, GUO Guangxia2,3,KONG Weidong2,*

    1SchoolofEnvironmentScienceandSpatialInformatics,ChinaUniversityofMiningandTechnology,Xuzhou221116,China2KeyLaboratoryofAlpineEcologyandBiodiversity,InstituteofTibetanPlateauResearch,ChineseAcademyofSciences,Beijing100101,China3StateKeyLaboratoryofLakeScienceandEnvironment,NanjingInstituteofGeographyandLimnology,ChineseAcademyofSciences,Nanjing210008,China

    Globally, climate change has become one of the most studied research topics. The Qinghai-Tibet Plateau covers an area of 2.5 ×106km2, with grassland being the dominant landscape. The plateau represents a typical ecologically fragile zone and has been experiencing rapid warming and precipitation change. Climate change has resulted in massive glacier melting and hydrology cycling change on the plateau. The warming and altered precipitation could profoundly alter soil microbial abundances and community structure, and subsequently influence ecosystem functions and biogeochemical cycling of carbon and nitrogen. In particular, this could influence microbe-driven greenhouse gas emissions in soils, e.g., N2O. However, responses of soil microbial community structure and functions to the warming and altered precipitation remain largely unexplored on the plateau, especially for the soil microorganisms involved in nitrogen cycling. The present study explored the effects of warming and altered precipitation on the denitrifying bacterial community in alpine meadow soils at Naqu research station on the Qinghai-Tibet Plateau. The effects of short-term (three months)warming (+1℃, +2℃,and +4 ℃)and precipitation change (50% precipitation decrease and 100% precipitation increase)on denitrifying bacteria (nirS) abundance and community structure were characterized using quantitative polymerase chain reaction (qPCR) and terminal restriction fragment length polymorphism (T-RFLP). The qPCR results showed that temperature increase in all treatments did not significantly affectnirSgene abundances, and both precipitation treatments did not significantly influencenirSgene abundances in the control and 2℃ warming treatments. In contrast, the 4℃ warming coupled with 100% precipitation increase treatments significantly elevatednirSgene abundances, suggesting that the interaction of warming and precipitation increase influenced the gene abundance at higher temperatures. Principal component analysis based on T-RFLP data demonstrated that temperature increase substantially altered thenirScommunity structure under conditions of no precipitation change and precipitation increase. ThenirScommunity structure of the control and 100% precipitation increase treatments showed similar patterns, while the 50% precipitation decrease treatment demonstrated different community patterns from the above two treatments, suggesting that thenirScommunity structure was sensitive to precipitation decrease. Canonical correspondence analysis further revealed that thenirScommunity structure was jointly driven by temperature increase and precipitation change,indicating that the interaction of temperature increase and precipitation contributed to the shift of thenirScommunity structure.The temperature increase explained thenirScommunity structure variation by 54.2% (P<0.01),precipitation change was explained by 45.5% (P<0.05), and the two jointly explained by 19.9%. Our results indicate that short-term warming and precipitation change do not influencenirSgene abundance; however, they do substantially shift thenirScommunity structure.

    nirS; gene abundance; community structure; climate change; Qinghai-Tibet Plateau

    國家自然科學(xué)基金青年科學(xué)基金項(xiàng)目(41401287);湖泊與環(huán)境國家重點(diǎn)實(shí)驗(yàn)室開放基金項(xiàng)目(2014SKL010)

    2016- 10- 18; < class="emphasis_bold">網(wǎng)絡(luò)出版日期

    日期:2017- 08- 14

    *通訊作者Corresponding author.E-mail: wdkong@itpcas.ac.cn

    10.5846/stxb201610182116

    潘曉悅,王曉,郭光霞,孔維棟.增溫與降水變化對青藏高原高寒草甸土壤nirS反硝化菌群落豐度和群落結(jié)構(gòu)的影響.生態(tài)學(xué)報(bào),2017,37(23):7938- 7946.

    Pan X Y,Wang X, Guo G X,Kong W D.Effects of increased temperature and precipitation change onnirSgene abundance and community structure in alpine meadow soils on the Qinghai-Tibet Plateau.Acta Ecologica Sinica,2017,37(23):7938- 7946.

    猜你喜歡
    青藏高原硝化群落
    青藏高原上的“含羞花”
    大學(xué)生牙齦炎齦上菌斑的微生物群落
    合成微生物群落在發(fā)酵食品中的應(yīng)用研究
    為了讓青藏高原的天更藍(lán)、水更綠、草原更美
    MBBR中進(jìn)水有機(jī)負(fù)荷對短程硝化反硝化的影響
    厭氧氨氧化與反硝化耦合脫氮除碳研究Ⅰ:
    春季和夏季巢湖浮游生物群落組成及其動態(tài)分析
    海水反硝化和厭氧氨氧化速率同步測定的15N示蹤法及其應(yīng)用
    化石見證青藏高原的隆起
    長期施用復(fù)合微生物肥對農(nóng)田土壤微生物群落的影響
    河南科技(2014年18期)2014-02-27 14:14:54
    日本成人三级电影网站| 人人妻人人看人人澡| 在线播放无遮挡| 国产精品一区二区性色av| 不卡一级毛片| 国产探花在线观看一区二区| 精品一区二区免费观看| 小说图片视频综合网站| 简卡轻食公司| 最近视频中文字幕2019在线8| 亚洲美女视频黄频| 久久久成人免费电影| 亚州av有码| 亚洲av成人av| 性插视频无遮挡在线免费观看| 日韩高清综合在线| 国产精品蜜桃在线观看 | 久久欧美精品欧美久久欧美| 最近视频中文字幕2019在线8| 一进一出抽搐动态| 日韩在线高清观看一区二区三区| 神马国产精品三级电影在线观看| 国产精品女同一区二区软件| 国产v大片淫在线免费观看| 国产精品不卡视频一区二区| 亚洲精品乱码久久久久久按摩| 中文精品一卡2卡3卡4更新| 97热精品久久久久久| 国产91av在线免费观看| 黄色日韩在线| 久久中文看片网| 亚洲精品久久国产高清桃花| 国内精品久久久久精免费| 国产熟女欧美一区二区| 五月玫瑰六月丁香| av女优亚洲男人天堂| 天堂网av新在线| 成人特级黄色片久久久久久久| 国产蜜桃级精品一区二区三区| 日本av手机在线免费观看| 啦啦啦啦在线视频资源| 老熟妇乱子伦视频在线观看| 丰满的人妻完整版| 亚洲四区av| 国产色爽女视频免费观看| 变态另类成人亚洲欧美熟女| 爱豆传媒免费全集在线观看| 乱系列少妇在线播放| 中文字幕久久专区| 午夜福利在线在线| 亚洲aⅴ乱码一区二区在线播放| 免费看日本二区| 欧美人与善性xxx| 亚洲乱码一区二区免费版| 狂野欧美激情性xxxx在线观看| 久99久视频精品免费| 大又大粗又爽又黄少妇毛片口| 精品熟女少妇av免费看| 国产一区二区亚洲精品在线观看| 在线免费十八禁| 免费观看精品视频网站| 此物有八面人人有两片| 国产黄片视频在线免费观看| 精品久久久久久久久久久久久| 欧美激情久久久久久爽电影| 久久久久久伊人网av| 一边摸一边抽搐一进一小说| 99精品在免费线老司机午夜| 22中文网久久字幕| 超碰av人人做人人爽久久| 久久精品久久久久久久性| 欧美日韩国产亚洲二区| 日韩欧美一区二区三区在线观看| 欧美色视频一区免费| 天堂网av新在线| 成人毛片a级毛片在线播放| 国产精品人妻久久久影院| 给我免费播放毛片高清在线观看| 精品久久久久久久人妻蜜臀av| 天天躁夜夜躁狠狠久久av| 色综合站精品国产| 国产视频内射| 看十八女毛片水多多多| 亚洲av.av天堂| 性色avwww在线观看| 国产在视频线在精品| 国产亚洲精品久久久久久毛片| 免费电影在线观看免费观看| 乱系列少妇在线播放| 日韩欧美在线乱码| 亚洲欧洲国产日韩| 老司机影院成人| 可以在线观看毛片的网站| 精品一区二区三区人妻视频| 亚洲最大成人av| 国产视频内射| 人人妻人人澡人人爽人人夜夜 | 在线观看66精品国产| 美女国产视频在线观看| 亚洲国产色片| 久久久久免费精品人妻一区二区| 精品人妻偷拍中文字幕| 成人亚洲欧美一区二区av| 99riav亚洲国产免费| 久久6这里有精品| 卡戴珊不雅视频在线播放| 日韩 亚洲 欧美在线| 久久这里有精品视频免费| 中文字幕av成人在线电影| 日韩强制内射视频| 国产91av在线免费观看| 悠悠久久av| 亚洲国产欧美人成| 极品教师在线视频| 欧美区成人在线视频| 亚洲精品色激情综合| 国产人妻一区二区三区在| 97超碰精品成人国产| 最好的美女福利视频网| 久久午夜福利片| 久久人妻av系列| 色哟哟哟哟哟哟| 国产视频首页在线观看| 特大巨黑吊av在线直播| 人人妻人人澡人人爽人人夜夜 | 国产免费男女视频| 国产伦精品一区二区三区视频9| 国产日本99.免费观看| 国产亚洲av片在线观看秒播厂 | 嫩草影院入口| 九九热线精品视视频播放| 色播亚洲综合网| 性插视频无遮挡在线免费观看| 国产麻豆成人av免费视频| 国产又黄又爽又无遮挡在线| 亚洲在久久综合| 亚洲激情五月婷婷啪啪| 免费不卡的大黄色大毛片视频在线观看 | 日本在线视频免费播放| 免费看光身美女| av天堂在线播放| 国产一区二区激情短视频| 国内精品宾馆在线| 国产激情偷乱视频一区二区| 自拍偷自拍亚洲精品老妇| 两个人视频免费观看高清| 国产综合懂色| 2021天堂中文幕一二区在线观| 国产高清有码在线观看视频| 国内精品久久久久精免费| av天堂在线播放| 夜夜爽天天搞| 国产在线男女| 九九久久精品国产亚洲av麻豆| 男人的好看免费观看在线视频| 国产在线男女| 美女国产视频在线观看| 国产69精品久久久久777片| 老熟妇乱子伦视频在线观看| 亚洲精品国产成人久久av| 欧美极品一区二区三区四区| 99久久久亚洲精品蜜臀av| 在线观看美女被高潮喷水网站| 免费黄网站久久成人精品| 午夜福利在线观看免费完整高清在 | 99久久人妻综合| 国产精品一区二区性色av| 免费大片18禁| 午夜精品一区二区三区免费看| 免费看a级黄色片| 18禁在线无遮挡免费观看视频| 成人综合一区亚洲| 在线观看一区二区三区| 日韩亚洲欧美综合| 精品少妇黑人巨大在线播放 | 毛片女人毛片| 不卡视频在线观看欧美| 亚洲精品粉嫩美女一区| 别揉我奶头 嗯啊视频| 一区二区三区高清视频在线| 免费搜索国产男女视频| 中文字幕av在线有码专区| 国内揄拍国产精品人妻在线| 人妻系列 视频| 免费av不卡在线播放| 变态另类丝袜制服| 亚洲久久久久久中文字幕| 高清午夜精品一区二区三区 | 国语自产精品视频在线第100页| 国产男人的电影天堂91| 亚洲精品乱码久久久v下载方式| 成人毛片60女人毛片免费| 免费看美女性在线毛片视频| 又黄又爽又刺激的免费视频.| 神马国产精品三级电影在线观看| 中国美女看黄片| 国产精品人妻久久久久久| 99久久精品热视频| 亚洲一区高清亚洲精品| 边亲边吃奶的免费视频| av卡一久久| 免费人成在线观看视频色| 黄片wwwwww| 1024手机看黄色片| 国产精品一区二区在线观看99 | 亚洲精品久久久久久婷婷小说 | 国产精品美女特级片免费视频播放器| av天堂中文字幕网| 免费观看人在逋| 成人美女网站在线观看视频| 村上凉子中文字幕在线| 亚洲人与动物交配视频| 亚洲中文字幕一区二区三区有码在线看| av在线天堂中文字幕| 天堂影院成人在线观看| 欧美日韩在线观看h| 亚洲精品国产av成人精品| 欧美日本视频| 久久精品影院6| 蜜桃亚洲精品一区二区三区| 国产乱人视频| 男插女下体视频免费在线播放| 成人永久免费在线观看视频| 免费不卡的大黄色大毛片视频在线观看 | 欧美性猛交╳xxx乱大交人| 午夜久久久久精精品| 亚洲精品久久久久久婷婷小说 | 99热这里只有是精品在线观看| 岛国毛片在线播放| 婷婷六月久久综合丁香| www.av在线官网国产| av天堂中文字幕网| 少妇丰满av| 国产精品久久久久久av不卡| 国产精品一区www在线观看| 国产色爽女视频免费观看| 婷婷精品国产亚洲av| 国产精品无大码| a级一级毛片免费在线观看| 亚洲中文字幕日韩| 婷婷色综合大香蕉| 少妇裸体淫交视频免费看高清| 麻豆av噜噜一区二区三区| 欧美激情国产日韩精品一区| 亚洲自拍偷在线| 成人午夜高清在线视频| 精品99又大又爽又粗少妇毛片| 亚洲色图av天堂| 国产成人91sexporn| 亚洲va在线va天堂va国产| 99热这里只有是精品50| 两个人的视频大全免费| 国产一区二区在线av高清观看| 日韩亚洲欧美综合| 91精品一卡2卡3卡4卡| 亚洲自拍偷在线| 夜夜爽天天搞| 欧美3d第一页| 亚洲av二区三区四区| 亚洲欧美精品自产自拍| 99久国产av精品| 麻豆一二三区av精品| 久久久久久九九精品二区国产| 色综合站精品国产| 麻豆乱淫一区二区| 91狼人影院| eeuss影院久久| 六月丁香七月| 国产精品麻豆人妻色哟哟久久 | 麻豆国产av国片精品| 国产午夜精品一二区理论片| 久久99蜜桃精品久久| 欧美性感艳星| 内地一区二区视频在线| 欧美日韩一区二区视频在线观看视频在线 | 搡老妇女老女人老熟妇| 91狼人影院| 国产一区二区亚洲精品在线观看| 联通29元200g的流量卡| 国产成人91sexporn| 91精品国产九色| 日日啪夜夜撸| 国产精品美女特级片免费视频播放器| 一进一出抽搐gif免费好疼| 性色avwww在线观看| 国产精品国产高清国产av| 日韩强制内射视频| 久久久欧美国产精品| 国产一区亚洲一区在线观看| 国产在视频线在精品| 少妇被粗大猛烈的视频| 亚洲av不卡在线观看| 自拍偷自拍亚洲精品老妇| 精品人妻熟女av久视频| 精品人妻偷拍中文字幕| 熟女人妻精品中文字幕| 99久久久亚洲精品蜜臀av| 18禁裸乳无遮挡免费网站照片| 欧美不卡视频在线免费观看| 青春草亚洲视频在线观看| 免费人成在线观看视频色| 97超视频在线观看视频| 亚洲天堂国产精品一区在线| 九九久久精品国产亚洲av麻豆| 如何舔出高潮| 国产av在哪里看| 国产真实乱freesex| 九九在线视频观看精品| 国产精品久久视频播放| 亚洲经典国产精华液单| 好男人视频免费观看在线| 久久久久国产网址| 69人妻影院| 国产成人一区二区在线| 成人漫画全彩无遮挡| 99久久精品一区二区三区| 日韩欧美三级三区| 精品久久久久久久久久免费视频| av又黄又爽大尺度在线免费看 | 欧美变态另类bdsm刘玥| 欧美性猛交黑人性爽| 免费人成在线观看视频色| 国产成年人精品一区二区| 久久精品国产亚洲av天美| 97人妻精品一区二区三区麻豆| 秋霞在线观看毛片| 午夜福利在线在线| 麻豆av噜噜一区二区三区| 波多野结衣高清无吗| 午夜爱爱视频在线播放| 亚洲激情五月婷婷啪啪| 成人性生交大片免费视频hd| 成年版毛片免费区| 久久久久久久久久黄片| 欧美激情在线99| 午夜福利在线在线| 夜夜爽天天搞| av黄色大香蕉| 亚洲精品日韩在线中文字幕 | 成年女人看的毛片在线观看| 日本一二三区视频观看| 寂寞人妻少妇视频99o| 97热精品久久久久久| 国产精品1区2区在线观看.| 免费人成视频x8x8入口观看| 白带黄色成豆腐渣| 午夜福利在线观看吧| 国产精品.久久久| 成人午夜高清在线视频| 人人妻人人澡人人爽人人夜夜 | 国产成人91sexporn| 欧洲精品卡2卡3卡4卡5卡区| 久久精品国产亚洲av香蕉五月| 精品日产1卡2卡| 精品久久久久久久久亚洲| 性欧美人与动物交配| 午夜福利在线在线| 久久精品国产99精品国产亚洲性色| 国产亚洲5aaaaa淫片| 在线免费观看不下载黄p国产| 又粗又爽又猛毛片免费看| 亚洲国产色片| 神马国产精品三级电影在线观看| 桃色一区二区三区在线观看| av免费观看日本| 中出人妻视频一区二区| 看非洲黑人一级黄片| 一卡2卡三卡四卡精品乱码亚洲| 欧美bdsm另类| 亚洲精品乱码久久久v下载方式| 亚洲欧洲国产日韩| 乱码一卡2卡4卡精品| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 最好的美女福利视频网| 国产v大片淫在线免费观看| 亚洲精品成人久久久久久| 91精品国产九色| 黄片wwwwww| 国产精品久久电影中文字幕| av免费在线看不卡| 亚洲在线观看片| 国产av不卡久久| 啦啦啦韩国在线观看视频| АⅤ资源中文在线天堂| 亚洲av中文字字幕乱码综合| 欧美性感艳星| 婷婷亚洲欧美| 男插女下体视频免费在线播放| 欧美三级亚洲精品| 亚洲成人久久爱视频| 天天躁夜夜躁狠狠久久av| 午夜激情福利司机影院| 成人三级黄色视频| 色尼玛亚洲综合影院| 久久久久久久久大av| 国产在线精品亚洲第一网站| 国产精品日韩av在线免费观看| 麻豆一二三区av精品| 久久久精品欧美日韩精品| 伦精品一区二区三区| 少妇裸体淫交视频免费看高清| 亚洲欧美日韩东京热| 日日干狠狠操夜夜爽| 亚洲色图av天堂| 最后的刺客免费高清国语| 九九在线视频观看精品| 成人av在线播放网站| 男人舔奶头视频| 欧美色欧美亚洲另类二区| 日本与韩国留学比较| 日本黄色片子视频| 搡女人真爽免费视频火全软件| 亚洲国产精品久久男人天堂| 久久精品国产清高在天天线| 欧美潮喷喷水| 免费观看人在逋| 国产精品.久久久| 国产成人影院久久av| 精品一区二区三区人妻视频| 亚洲欧美中文字幕日韩二区| 插逼视频在线观看| 两个人视频免费观看高清| 色尼玛亚洲综合影院| 国产久久久一区二区三区| 午夜福利高清视频| 亚洲七黄色美女视频| 精品国产三级普通话版| 国产淫片久久久久久久久| 插逼视频在线观看| 午夜福利在线观看吧| 夫妻性生交免费视频一级片| 亚洲av第一区精品v没综合| 波多野结衣巨乳人妻| 午夜a级毛片| 久久久久性生活片| 成年女人看的毛片在线观看| 亚洲aⅴ乱码一区二区在线播放| 国产精品人妻久久久久久| 爱豆传媒免费全集在线观看| 老司机福利观看| 久久99精品国语久久久| 99久国产av精品国产电影| 国产高清不卡午夜福利| 天美传媒精品一区二区| 赤兔流量卡办理| 爱豆传媒免费全集在线观看| 亚洲美女搞黄在线观看| 美女大奶头视频| 午夜激情欧美在线| 三级国产精品欧美在线观看| 最近手机中文字幕大全| 欧美最新免费一区二区三区| 亚洲aⅴ乱码一区二区在线播放| 欧美激情国产日韩精品一区| 乱人视频在线观看| 人人妻人人看人人澡| 精品久久久久久久久av| 国产色爽女视频免费观看| 久久婷婷人人爽人人干人人爱| 久久草成人影院| 黄色视频,在线免费观看| 亚洲18禁久久av| 久久精品国产亚洲av天美| 美女大奶头视频| 国产单亲对白刺激| 久久人人精品亚洲av| 日本撒尿小便嘘嘘汇集6| 亚洲精品自拍成人| 日韩 亚洲 欧美在线| 亚洲精品亚洲一区二区| 亚洲人与动物交配视频| 一区二区三区高清视频在线| 在线免费十八禁| 人妻夜夜爽99麻豆av| 日本一二三区视频观看| 精品熟女少妇av免费看| av女优亚洲男人天堂| 男女啪啪激烈高潮av片| 一级黄片播放器| 亚洲精品久久久久久婷婷小说 | 欧美高清性xxxxhd video| 91aial.com中文字幕在线观看| 亚洲精品色激情综合| 国产毛片a区久久久久| 亚洲成av人片在线播放无| 成年av动漫网址| 直男gayav资源| 性欧美人与动物交配| 性欧美人与动物交配| 亚洲久久久久久中文字幕| 成熟少妇高潮喷水视频| 国产av一区在线观看免费| 观看免费一级毛片| 中文字幕久久专区| 青春草视频在线免费观看| 乱码一卡2卡4卡精品| 人妻久久中文字幕网| 国产大屁股一区二区在线视频| 成熟少妇高潮喷水视频| 老师上课跳d突然被开到最大视频| 久久精品国产鲁丝片午夜精品| 小蜜桃在线观看免费完整版高清| 成人欧美大片| 精品国内亚洲2022精品成人| 国产久久久一区二区三区| 乱系列少妇在线播放| 国产精品久久久久久精品电影小说 | 一本一本综合久久| 舔av片在线| 黑人高潮一二区| 亚洲18禁久久av| 一边亲一边摸免费视频| 国产一区二区亚洲精品在线观看| 日韩强制内射视频| 亚洲av第一区精品v没综合| 日韩av在线大香蕉| 成人无遮挡网站| 在线免费观看不下载黄p国产| 久久久久性生活片| 色视频www国产| 亚洲成人久久爱视频| 成人二区视频| 麻豆av噜噜一区二区三区| 色尼玛亚洲综合影院| 99精品在免费线老司机午夜| 午夜视频国产福利| 久久久久久久久久久免费av| 亚洲中文字幕日韩| 日韩欧美 国产精品| 男女下面进入的视频免费午夜| 欧美一区二区国产精品久久精品| 91久久精品国产一区二区成人| 亚洲最大成人手机在线| 成年免费大片在线观看| 成人欧美大片| 国内精品美女久久久久久| 最新中文字幕久久久久| 草草在线视频免费看| 国产精品无大码| 国产精品日韩av在线免费观看| 亚洲精品久久国产高清桃花| 国内精品久久久久精免费| 国产乱人视频| 一级毛片电影观看 | 日韩欧美国产在线观看| 免费看日本二区| 69人妻影院| 黄色欧美视频在线观看| 日韩欧美一区二区三区在线观看| 一级二级三级毛片免费看| 亚洲精品日韩av片在线观看| 一区二区三区免费毛片| 中文字幕久久专区| 国产乱人视频| 两个人视频免费观看高清| 午夜精品国产一区二区电影 | 亚洲真实伦在线观看| 噜噜噜噜噜久久久久久91| 99热这里只有是精品在线观看| 夫妻性生交免费视频一级片| 国产男人的电影天堂91| 国产视频内射| 欧美区成人在线视频| 国产一区二区三区在线臀色熟女| 99久久成人亚洲精品观看| 22中文网久久字幕| 99国产精品一区二区蜜桃av| 国产欧美日韩精品一区二区| 最新中文字幕久久久久| 最近的中文字幕免费完整| 一级黄色大片毛片| 久久久久久大精品| 内地一区二区视频在线| 成年女人看的毛片在线观看| 欧美性感艳星| 乱人视频在线观看| 国产高清三级在线| 国产高清激情床上av| 欧美丝袜亚洲另类| 特级一级黄色大片| 国产伦精品一区二区三区视频9| 国产亚洲av嫩草精品影院| 伊人久久精品亚洲午夜| 国产成人精品久久久久久| 69人妻影院| 国产精品电影一区二区三区| 在线免费观看的www视频| 免费黄网站久久成人精品| 亚洲精品日韩在线中文字幕 | 亚洲人与动物交配视频| 岛国毛片在线播放| 免费不卡的大黄色大毛片视频在线观看 | 国产黄片视频在线免费观看| 亚洲欧美精品自产自拍| 老司机福利观看| 九色成人免费人妻av| 国产精品一区二区三区四区免费观看| 国产淫片久久久久久久久| 色噜噜av男人的天堂激情| 中文字幕av在线有码专区| 成人无遮挡网站| 国产一区二区三区av在线 | 中国国产av一级| 欧美日韩综合久久久久久| 人妻系列 视频| 欧美日韩精品成人综合77777| 亚洲成人久久爱视频| 人体艺术视频欧美日本| 不卡一级毛片| 老司机影院成人| 国产乱人视频| 最近中文字幕高清免费大全6| 免费一级毛片在线播放高清视频| 此物有八面人人有两片| 国产黄色视频一区二区在线观看 |