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

    氮添加和水分改變對亞熱帶人工林紅壤微生物多樣性及群落組成的影響

    2024-04-28 05:46:54孫浩釗劉桂珍朱玉璘林雪萍郭豪錢笑杰方熊
    生態(tài)科學(xué) 2024年1期
    關(guān)鍵詞:施氮土壤水分含水量

    孫浩釗, 劉桂珍, 朱玉璘, 林雪萍, 郭豪, 錢笑杰, 方熊

    氮添加和水分改變對亞熱帶人工林紅壤微生物多樣性及群落組成的影響

    孫浩釗1, 劉桂珍1, 朱玉璘1, 林雪萍1, 郭豪2, 錢笑杰1, 方熊3,*

    1. 福建農(nóng)林大學(xué)資源與環(huán)境學(xué)院, 土壤環(huán)境健康與調(diào)控福建省重點實驗室, 福州 350000 2. 福建農(nóng)林大學(xué)機電工程學(xué)院, 福州 350000 3. 江西農(nóng)業(yè)大學(xué)國土資源與環(huán)境學(xué)院,南昌 330045

    氮沉降加劇與全球降雨格局改變是全球變化的重要特征, 外源氮素的輸入及土壤水分條件的變化勢必會影響土壤微生物的多樣性和群落組成, 進而影響地球化學(xué)過程。運用Illumina Miseq高通量測序技術(shù), 研究氮添加、水分改變及二者交互作用對亞熱帶人工林紅壤細(xì)菌、真菌多樣性、群落組成的影響及其與土壤理化因子間的關(guān)系。結(jié)果表明: 土壤細(xì)菌多樣性隨水分含量的降低而減少, 且氮添加增強了這一抑制作用。施氮后中水分處理土壤真菌Shannon指數(shù)顯著高于低水分處理。氮添加顯著降低Acidobacteria的相對豐度, 增加Actinobacteria、Bacteroidetes的相對豐度。土壤Acidobacteria、Bacteroidetes的相對豐度隨土壤含水量的升高而增加, Actinobacteria、Firmicutes則相反。干旱加劇了氮添加對Acidobacteria相對豐度的抑制作用, 氮水交互處理顯著影響Cyanobacteria的相對豐度。而氮添加、土壤水分及二者交互作用對土壤真菌群落組成均無顯著影響。冗余分析表明土壤細(xì)菌群落組成的變化與土壤含水量、pH、NO3–-N含量、TN含量顯著相關(guān), 分別解釋了土壤細(xì)菌群落組成變化的65.6%、26.1%、21.5%、20.7%。研究表明氮添加會改變干旱對土壤細(xì)菌多樣性和群落組成的影響的強度和方向。

    氮沉降; 土壤含水量; 微生物多樣性; 微生物群落組成

    0 前言

    氮沉降加劇及降水格局變化是全球變化的兩個重要特征, 當(dāng)前我國華南地區(qū)的氮沉降量已經(jīng)達到32—49 kg·N·ha–1·a[2]。與此同時, 全球降水格局也在發(fā)生著變化, 氣候變化預(yù)測模型進一步表明, 未來全球范圍內(nèi)將經(jīng)歷更極端、異常的干旱—降雨事件[3], 我國華南地區(qū)也將面臨季節(jié)性干旱的威脅[4]。

    土壤微生物是土壤生態(tài)系統(tǒng)的重要組成部分, 在土壤生物地球化學(xué)循環(huán)等過程中發(fā)揮著重要作用。氮沉降通過改變氮素有效性、pH值、土壤C/N, 改變植物與土壤微生物之間的養(yǎng)分分配等方式, 直接或間接地影響微生物多樣性及群落組成[5]。如Barbara等長期施氮肥的實驗發(fā)現(xiàn)土壤細(xì)菌多樣性隨施氮量的增加而降低[6]。但Xu等的研究結(jié)果表明土壤細(xì)菌多樣性隨施氮量(0、50、150、250 kg·N·ha–1·a–1)的增加而增加[7]。通過Meta分析的方法研究氮添加對土壤微生物多樣性影響, Wang等的研究結(jié)果發(fā)現(xiàn)低水平的施氮量對土壤細(xì)菌和真菌的Shannon指數(shù)均有負(fù)面影響, 且氮添加對細(xì)菌Shannon指數(shù)的抑制作用隨施氮量增加而增強; 真菌則不然, 當(dāng)施氮量大于10 kg·N·ha–1·a–1時, 氮添加反而增加了真菌Shannon指數(shù)[8]。由于不同的土壤微生物類群可能在土壤氮循環(huán)過程(如固氮、硝化、反硝化等)中發(fā)揮的功能不同, 因而不同的微生物類群對外源氮素輸入的響應(yīng)各異[9]。Moorhead等認(rèn)為: 施氮會增加土壤無機氮含量, 提高土壤碳的有效性, 進而有助于增加富營養(yǎng)型菌種的相對豐度, 減少貧營養(yǎng)型菌種的相對豐度, 因為隨著氮素限制的減弱, 微生物轉(zhuǎn)向分解不穩(wěn)定碳源, 從而降低了這類偏向分解難降解碳源的貧營養(yǎng)型菌群相對豐度[10]。例如Acidobacteria是一種貧營養(yǎng)型土壤細(xì)菌, 其相對豐度通常在氮添加后降低[11], 而Proteobacteria是一種富營養(yǎng)型土壤細(xì)菌, 其相對豐度通常在氮添加后增加[12]。施氮可能通過改變氮有效性直接影響微生物相對豐度, 同時也可以通過降低土壤pH間接影響土壤微生物的相對豐度[13–14]。

    土壤水分在土壤微生物生命活動中同樣起著重要作用。土壤水分可以通過直接影響微生物細(xì)胞內(nèi)滲透壓, 影響微生物的活性[15]; 或通過影響營養(yǎng)物質(zhì)在土壤界面的遷移, 間接影響微生物獲取養(yǎng)分[16]; 此外土壤水分的變化還會影響植物生理活動, 從而通過植物根系分泌物間接影響土壤微生物[17]。如Preece等利用溫室盆栽實驗研究土壤水分下降對土壤微生物的影響, 結(jié)果表明, 細(xì)菌Shannon指數(shù)隨干旱程度的加劇而降低[18]。同時Nicholas等在研究干旱對熱帶雨林土壤微生物的影響時也有類似的發(fā)現(xiàn), 干旱組的細(xì)菌多樣性顯著低于對照組[19]。但Acosta—Martínez等在干旱和熱浪對土壤微生物影響的研究中發(fā)現(xiàn), 干旱處理中土壤細(xì)菌、真菌Shannon指數(shù)均有顯著提升[20]。岳丹丹等通過盆栽控水研究干旱對棉花根系土壤真菌多樣性的影響, 發(fā)現(xiàn)與正常澆水處理相比, 干旱處理增加了真菌的Shannon指數(shù)和Simpson指數(shù)[21]。土壤微生物群落中, 具有不同生理和生態(tài)策略的微生物類群對土壤水分有效性變化的響應(yīng)不同, 有研究表明Proteobacteria與土壤含水量顯著正相關(guān), 土壤含水量的增加會增加其相對豐度, 而Acidobacteria相對豐度會隨著土壤干旱程度的加劇持續(xù)減少, Gemmatimonadetes則與Acidobacteria相反[22–23]。

    我國的亞熱帶地區(qū), 氮沉降量高, 降水季節(jié)性變化大, 生態(tài)系統(tǒng)服務(wù)功能強[24]。在氮沉降增加和干旱加劇背景下, 氮素和土壤水分的交互作用將如何影響土壤微生物多樣性和群落組成, 氮素與土壤水分對土壤微生物多樣性和群落組成的影響是否存在協(xié)同或拮抗作用還不清楚[25]?;诖? 本實驗以亞熱帶人工林紅壤為研究對象, 模擬研究氮添加、土壤水分及二者交互作用對土壤細(xì)菌和真菌多樣性及群落組成的影響。本研究結(jié)果有助于加深區(qū)域氮、水格局變化對森林土壤微生物群落影響的理解。

    1 材料與方法

    1.1 研究方法

    實驗地點位于福州市福建農(nóng)林大學(xué)實驗溫室(26°08′N, 119°24′E), 于2017年10月在溫室內(nèi)建設(shè)18個生長池(1.4 m×1.4 m×0.6 m)。同年11月將從城郊人工林運來的紅壤混勻并填入池中, 填土深度為0.5 m。靜置穩(wěn)定1個月后將馬尾松()、杉木()、海南紅豆()和木荷()1年生幼苗各一株移入池中, 其中杉木和馬尾松為針葉植物, 木荷和紅豆為闊葉植物, 此外海南紅豆還具有固氮作用。緩苗4個月, 于2018年4月進行氮水交互處理, 根據(jù)實驗地區(qū)氮沉降水平, 設(shè)置兩組氮處理, 分別為不施氮(CK, 0 kg·N·ha–1·a–1)和施氮(AN, 80kg·N·ha–1·a–1); 每個氮處理下設(shè)置3個水分梯度, 分別為低水分(LW, 約40%田間最大持水量)、中水分(MW, 約60%田間最大持水量)和高水分(HW, 約80%田間最大持水量)處理, 在實際實驗處理中分別相當(dāng)于約15%、22%、28%土壤體積含水率。實驗共設(shè)計6個處理, 各處理設(shè)置3個重復(fù), 各處理和重復(fù)隨機分布。以硝酸銨為氮源, 每月定期將1.905 g硝酸銨溶于5 L水中噴于施氮處理中, 未施氮組噴水5 L。定期測定土壤含水率并補水分至實驗控制水平。處理8個月后, 每個樣方隨機采集0—10 cm土壤樣品4份并充分混勻,過2 mm篩后, 一份用于土壤理化性質(zhì)的分析, 一份保存于–80 ℃冰箱用于分析土壤微生物。

    1.2 樣品分析

    1.2.1 土壤基本理化性質(zhì)分析

    土壤pH使用電位法, 水土比為2.5︰1。使用元素分析儀(Vario CUBE, 艾力蒙塔, 德國)分析土壤總碳(TC)、總氮(TN)含量。土壤銨態(tài)氮(NH4+-N)、硝態(tài)氮(NO3–-N)含量用氯化鉀溶液浸提后使用間斷性流動分析儀(Smart Chem200, MS-Allianc, 意大利)測得。土壤總磷(TP)使用硫酸—高氯酸消煮—鉬銻抗比色法測得, 使用鹽酸—硫酸浸提法測定土壤有效磷含量(AP)。速效鉀(AK)使用乙酸銨浸提—火焰光度法測定, 具體實驗步驟參考《土壤理化分析與剖面描述》[26]。

    1.2.2 土壤微生物群落分析(PCR-AGE)

    本研究采用土壤DNA提取試劑盒E.Z.N.A.? Soil DNA Kit(OMEGA, D5625-01)提取土壤微生物DNA。PCR采用15 μL 2×Phanta Master Mix、引物各1 μL、10—20 ng Genomic DNA、30 μL ddH2O。細(xì)菌采用引物為341F(5’—CCTACGGGNGGCWGC AG—3’)和806R(5’—GGACTACHVGGGTWTCTA AT—3’), 真菌采用引物為ITS1F(5’—CTTGGTCA TTTAGAGGAAGTAA—3’)和ITS1R(5’—GCTGCG TTCTTCATCGATGC—3’)。PCR反應(yīng)程序為: 95 ℃預(yù)變性(5 min), 95 ℃變性(30 s), 55 ℃退火(30 s), 72 ℃延伸(45 s)(以上其中細(xì)菌循環(huán)27次、真菌循環(huán)29次); 72 ℃終延伸(10 min)。PCR結(jié)束后電泳驗證并引入Illumina橋式PCR兼容引物進行第二輪擴增, 以16S rRNA基因產(chǎn)物為細(xì)菌第二輪擴增模板, 以ITS rRNA基因產(chǎn)物為真菌第二輪擴增模板, 選擇細(xì)菌特異性較高的V3—V4區(qū), 真菌特異性較高的ITS1區(qū)作為擴增片段, 引物為F(AATGATACGG CGACCACCGAGATCTACACTATAGCCTACACTC TTTCCCTACACGACGCTCTTCCGATCT)和R(CAA GCAGAAGACGGCATACGAGATCGAGTAATGTG ACTGGAGTTCAGACGTGTGCTCTTCCGATCT), 反應(yīng)過程為95 ℃預(yù)變性(30 s), 95 ℃變性(15 s), 55℃退火(15 s), 72 ℃延伸(30 s), 5個循環(huán); 72 ℃終延伸(5 min)。PCR結(jié)束后, PCR產(chǎn)物進行瓊脂糖電泳檢測, 瓊脂糖濃度: 2%, 電壓80 V, 電泳時間40 min。選擇EqualbitTM dsDNA HSAssayK-it (Vazyme, EQ111-01)試劑盒切膠回收PCR產(chǎn)物, 由生物工程公司進行測序。

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

    使用USEARCH在0.97相似度下進行聚類, 對聚類后的序列進行嵌合體過濾后, 得到用于物種分類的OUT序列, 通過比對silva rRNA數(shù)據(jù)庫, 得到每個OTU對應(yīng)的物種分類信息。采用香農(nóng)指數(shù)(Shannon index)和辛普森指數(shù)(Simpson index)計算土壤細(xì)菌、真菌群落的均勻性及豐富度, 表征土壤微生物群落多樣性, 計算方法參考郭萍萍等[27]。使用R語言的vegan包基于Weigted Unifrac距離矩陣對所有OUT進行NMDS分析。使用SPSS26.0雙因素方差分析法檢驗氮添加、土壤水分及二者交互作用對土壤理化性質(zhì)、微生物多樣性、微生物門水平優(yōu)勢種群相對豐度的影響。采用單因素方差分析檢驗相同氮水平不同水分處理之間或相同水分水平不同氮處理之間的土壤理化性質(zhì)、微生物多樣性、微生物優(yōu)勢種群相對豐度的顯著性, 用LSD法進行多重比較。采用雙變量相關(guān)性分析對土壤相對豐度占比超過1%的優(yōu)勢菌群與土壤環(huán)境因子之間的相關(guān)性進行分析, 顯著性水平設(shè)置為0.05。使用Canoco 4.5對相對豐度占比大于1%的微生物做DCA (Detrended corresp-ondence analysis)分析, 前四軸最大值小于3, 選擇RDA(Redundancy analysis)分析其與環(huán)境因子之間關(guān)系。使用Origin軟件繪制柱狀圖。

    2 結(jié)果與分析

    2.1 土壤理化性質(zhì)

    總體上看, 氮添加顯著增加土壤TN、NH4+-N、NO3–-N的含量(表1)。不同水分條件下, 氮添加對土壤TN含量無顯著影響, 均顯著增加土壤NO3–-N含量, 但僅顯著增加了低水分處理組土壤NH4+-N的含量。水分處理顯著影響土壤TN和NH4+-N含量, 同時對土壤pH(=0.051)和NO3–-N含量(=0.067)的影響也幾乎達到顯著水平。土壤TN含量表現(xiàn)為低水分處理組高于(<0.05)中水分處理組。隨著土壤水分降低, 土壤pH有降低趨勢。氮添加與土壤水分交互作用顯著影響土壤NH4+-N、NO3–-N的含量(<0.01)。在不加氮組, 土壤NO3–-N含量隨土壤含水量的降低而降低; 但是加氮后NH4+-N、NO3–-N含量隨土壤含水量的降低而增加。在不加氮組, 水分處理對土壤pH無顯著影響; 但在加氮組, 低水分處理顯著降低土壤pH。

    2.2 土壤微生物多樣性

    如圖1所示: 整體上看, 氮添加對土壤細(xì)菌多樣性無顯著影響, 而水分處理顯著影響細(xì)菌的多樣性, 表現(xiàn)為低水分處理細(xì)菌多樣性顯著低于中、高水分處理(<0.01)。在不加氮組, 水分對土壤細(xì)菌多樣的無顯著影響; 但是加氮后, 低水分處理顯著降低土壤細(xì)菌多樣性(<0.05)。氮、水及二者的交互處理對土壤真菌多樣性無顯著影響, 但中水分處理土壤真菌Shannon指數(shù)顯著大于低水分處理。

    2.3 土壤微生物群落組成

    基于Weighted Unifrac distance的微生物群落組成非度量多維尺度分析(NMDS)分析氮添加、土壤水分及二者交互作用對土壤微生物群落組成的影響(圖2), 結(jié)果顯示, 施氮使低水分與中水分處理的細(xì)菌群落分開, 但真菌存在明顯的重疊。不同水分處理的土壤細(xì)菌群落在第一軸上分開(圖2a), 低水分處理和高水分處理的土壤真菌群落在第二軸上分開(圖2b); 總體而言, 土壤微生物群落對土壤含水量的響應(yīng)大于施氮, 而氮添加對細(xì)菌群落組成的影響主要體現(xiàn)在低水分處理。

    表1 不同處理的土壤理化性質(zhì)

    注: 表中土壤理化性質(zhì)數(shù)值為平均值±標(biāo)準(zhǔn)差,=3。不同小寫字母表示相同施氮處理不同土壤水分間差異顯著(<0.05), 不同大寫字母表示相同土壤水分不同施氮處理間差異顯著(<0.05)。N, W, N*W的值為雙因素方差分析的值, 分別表示施氮、土壤水分及二者交互作用對土壤理化性質(zhì)的影響, **表示差異極顯著(<0.01)。

    注: 平均值±標(biāo)準(zhǔn)差, n=3, 不同小寫字母表示相同施氮處理不同土壤水分間差異顯著(P<0.05); 圖中CK表示未施氮處理, AN表示氮添加處理, LW、MW、HW分別表示低、中、高水分處理。N、W、N*W分別表示施氮、土壤水分及二者交互作用對細(xì)菌, 真菌群落shannon指數(shù)、Simpson指數(shù)的影響, P<0.05時影響顯著, 下同。

    Figure 1 Shannon index (lift) and Simpson index (right) of soil bacteria (a, b) and fungi (c, d) under different treatments

    注: CL表示不施氮低水分處理, NL表示施氮低水分處理, CM表示不施氮中水分處理, NM表示施氮中水分處理, CH表示不施氮高水分處理, NH表示施氮高水分處理, 下同。

    Figure 2 Non-metric multidimensional scaling (NMDS) of bacterial (a) and fungal (b) community structure based on Weighted Unifrac distance

    選擇相對豐度在1%以上的門水平優(yōu)勢微生物種群分析其相對豐度的變化, 細(xì)菌分別為變形菌門(Proteobacteria)、酸桿菌門(Acidobacteria)、放線菌門(Actinobacteria)、綠彎菌門(Chloroflexi)、厚壁菌門(Firmicutes)、芽單胞菌門(Gemmatimonadetes)、藍(lán)藻菌門(Cyanobacteria)和擬桿菌門(Acidobacteria), 其中變形菌門相對豐度最高(29.54%—37.90%), 其次是酸桿菌(14.90%—21.69%)、放線菌(5.60%— 19.15%) (圖3)。真菌分別為子囊菌門(Ascomycota)、擔(dān)子菌門(Basidiomycota)、球囊菌門(Mortierellomycota)和羅茲菌門(Rozellomycota)。其中子囊菌門相對豐度最高(59.74%—82.76%), 其次是擔(dān)子菌門(10.12%— 29.07%)、球囊菌門(1.25%—8.80%)(圖4)。如圖3所示: 氮添加顯著增加Actinobacteria、Bacteroidetes的相對豐度(<0.05), 降低Acidobacteria的相對豐度(<0.05)。土壤水分條件的改變顯著影響除Chloroflexi、Cyanobacteria外其余6種優(yōu)勢細(xì)菌相對豐度, 具體表現(xiàn)為: Acidobacteria、Bacteroidetes的相對豐度隨土壤含水量的升高而升高(<0.05); 但Actinobacteria (<0.01)、Firmicutes (<0.05)、Gemmatimonadetes (<0.05)的相對豐度隨土壤含水量的降低而降低。氮添加使低水分處理組Cyanobacteria的相對豐度顯著下降(<0.05), 但隨著土壤水分的增加, Cyanobacteria的相對豐度對氮添加的響應(yīng)由降低變?yōu)樯摺3酥? 不施氮處理組, 中、高水分處理下Proteobacteria、Bacteroidetes的相對豐度顯著高于低水分處理(<0.05), 但其顯著性在氮添加后消失。土壤水分的增加對Firmicutes相對豐度的抑制作用也在氮添加后達到了顯著水平。施氮、土壤水分及二者的交互作用對土壤門水平優(yōu)勢真菌的影響均未表現(xiàn)出顯著性(圖4), 但干旱有增加Ascomycota和降低Basidiomycota相對豐度的趨勢。

    注: 不同大寫字母表示相同土壤水分不同施氮處理間差異顯著(P<0.05)。

    Figure 3 Relative abundance of dominant soil bacteria in phylum level under different treatments

    2.4 土壤微生物群落與環(huán)境因子的相關(guān)性分析

    土壤細(xì)菌RDA前兩個排序軸分別解釋73.0%、10.1%的土壤細(xì)菌優(yōu)勢種群與環(huán)境因子之間的關(guān)系(圖5a), 土壤含水量(=0.002)、pH(=0.004)、NO3–-N (=0.032)和TN(=0.022)與細(xì)菌群落組成顯著相關(guān), 分別解釋了土壤細(xì)菌群落組成變化的65.6%、26.1%、21.5%和20.7%(圖5a)。土壤真菌冗余分析前兩個排序軸分別解釋38.2%、2.2%的土壤真菌優(yōu)勢種群與環(huán)境因子之間的關(guān)系(圖5b), 土壤理化性質(zhì)對真菌優(yōu)勢菌群的影響不顯著。

    3 討論

    3.1 氮添加和水分改變對土壤細(xì)菌、真菌多樣性的影響

    整體上看, 氮添加及氮水交互作用對土壤細(xì)菌的多樣性無顯著影響, 僅土壤含水量對土壤細(xì)菌多樣性有顯著影響。已有許多研究發(fā)現(xiàn)施氮對土壤細(xì)菌多樣性有明顯的降低作用[8], 但也有研究表明: 在氮限制環(huán)境下, 短期施氮會增加細(xì)菌活性, 少量的氮能夠促進某些土壤細(xì)菌的生長[28]。在本研究中, 施氮有增加細(xì)菌多樣性的趨勢, 但是沒有達到顯著水平(圖1)。Christian等人認(rèn)為在微酸性土壤中, 土壤細(xì)菌多樣性與土壤pH值之間存在很強的相關(guān)關(guān)系[29]。因此在本研究中可能是由于氮添加量及實驗處理時間的原因, 施氮并沒有對土壤pH產(chǎn)生顯著影響(表1), 因此也未顯著影響細(xì)菌多樣性。土壤水分在土壤細(xì)菌生命活動中起重要作用, 一般來說, 土壤水分降低會影響細(xì)菌體內(nèi)滲透壓及土壤養(yǎng)分有效性因而對土壤細(xì)菌的多樣性有負(fù)面影響[16, 30]。這與本實驗的研究結(jié)果低水分處理土壤細(xì)菌的Shannon指數(shù)顯著低于中、高水分處理一致。本研究還發(fā)現(xiàn), 增加土壤含水量后, 氮添加對土壤細(xì)菌多樣性的影響由降低變?yōu)樯? 且施氮只顯著降低了低水分處理土壤pH值。Zhang等的實驗也發(fā)現(xiàn)了類似的結(jié)果, 增加土壤含水量增加了低施氮水平處理的細(xì)菌多樣性[31]。一方面這可能與土壤pH對細(xì)菌多樣性的影響有關(guān), 土壤水分含量增加后, 增強了土壤硝態(tài)氮的淋失, 此外土壤含水量的增加也會增強土壤反硝化作用, 從而導(dǎo)致酸性土壤pH升高[32];另一方面, 還可能是因為土壤水分的增加對細(xì)菌多樣性的促進作用。土壤水分增加可能會增強短期內(nèi)土壤氮素對土壤細(xì)菌多樣性的促進作用。

    圖4 不同處理土壤真菌門水平優(yōu)勢種群的相對豐度

    Figure 4 Relative abundance of dominant soil fungi in phylum level under different treatments

    注: 各種圖形含義同上, 理化性質(zhì)實線表示影響顯著, 虛線表示不顯著。

    Figure 5 Redundancy analysis results for bacterial community and soil environmental parameters (Bacteria a, Fungi b)

    本研究發(fā)現(xiàn)土壤真菌多樣性對土壤氮水條件變化的耐受性強于細(xì)菌, 施氮、土壤水分及二者的交互作用對土壤真菌多樣性均無顯著作用。施氮增加了真菌Shannon指數(shù)與Simpson指數(shù), 但未達到顯著水平。Zhou等在中國東北黑土上連續(xù)施用了34年氮肥, 研究其對土壤真菌的影響, 結(jié)果發(fā)現(xiàn)真菌多樣性顯著降低[33]。這與本實驗研究結(jié)果不符, 可能與實驗處理時間長短與實驗地區(qū)土壤類型有關(guān)[8, 34–35]。Zhang等人在研究施氮與降雨增加對內(nèi)蒙古草原土壤真菌影響的實驗中發(fā)現(xiàn): 正常降水條件真菌多樣性在施氮量為100 kg·N·ha–1·a–1時顯著升高, 但增加降水后其顯著性消失, 并推測施氮與增雨相結(jié)合有利于維持土壤真菌的多樣性[36]。這與本研究的結(jié)果相似, 本實驗發(fā)現(xiàn): 施氮后, 中水分處理的真菌Shannon指數(shù)顯著高于低水分處理, 但高水分處理與低水分處理之間并無顯著性。這可能是由于中土壤水分條件下更適宜土壤真菌的生命活動, 也可能由于細(xì)菌與真菌之間對有機碳的競爭所導(dǎo)致[37], 本結(jié)果也進一步說明氣候變化對土壤真菌多樣性的影響可能要從多因素交互作用預(yù)測。

    3.2 氮添加和水分改變對土壤細(xì)菌、真菌群落組成的影響

    氮素是影響土壤細(xì)菌群落組成的重要因素, 不同類群土壤微生物對施氮的響應(yīng)不同[38]。畢捷等研究了半干旱草原土壤微生物群落組成對氮水條件變化的響應(yīng), 發(fā)現(xiàn)其群落組成受到氮添加的顯著影響[39]。本研究也發(fā)現(xiàn)類似的結(jié)果, 這可能是因為外源氮素的輸入改變了土壤微生物的底物供應(yīng), 本研究中氮添加使土壤NH4+-N、NO3–-N含量顯著增加, 降低低水分處理組土壤pH(表1), 且RDA的結(jié)果表明土壤pH、NO3–-N和TN含量與細(xì)菌群落組成顯著相關(guān)(圖5a)。另一方面, 本實驗結(jié)果表明氮添加顯著降低Acidobacteria的相對豐度, 顯著增加Actinobacteria的相對豐度。這與Zeng等的研究結(jié)果一致, 他們的結(jié)果表明, Acidobacteria的相對豐度隨施氮量的升高而下降、Actinobacteria隨施氮量的升高而升高, 這可能是因為氮添加提高了土壤的養(yǎng)分有效性, 直接影響微生物的生理活動, 也可能由于氮添加降低了土壤pH或通過植物根系等作用間接影響其相對豐度[40]。此外, 本研究還發(fā)現(xiàn)土壤水分的變化改變了Cyanobacteria、Actinobacteria的相對豐度對氮添加的響應(yīng), 低水分處理中, 施氮顯著降低了Cyanobacteria的相對豐度, 增加土壤含水量后, 施氮反而增加了Cyanobacteria的相對豐度, 有研究表明: 水分是影響土壤微藻分布的重要限制因素, 在一些水分充足的土壤環(huán)境中, 藍(lán)藻的豐富度最高, 并且大多數(shù)藍(lán)藻都具有固氮功能, 能將無機氮轉(zhuǎn)化為有機氮[41], 氮源豐富的環(huán)境可能會促進其生長。張海芳等的研究結(jié)果表明, 增加土壤水分導(dǎo)致Actinobacteria相對豐度在施氮后由顯著升高變?yōu)轱@著降低[42], 這與本實驗結(jié)果相似, 本研究中Actinobacteria的相對豐度與土壤含水量極顯著負(fù)相關(guān), 與土壤硝態(tài)氮的含量顯著正相關(guān)(圖5a), 高土壤水分抑制了施氮對Actinobacteria生長的促進作用, 這可能與Actinobacteria在不同水分條件下的形態(tài)特征有關(guān), 在不同的水分條件下, Actinobacteria對氮添加的響應(yīng)存在差異。

    四種土壤優(yōu)勢真菌的相對豐度對施氮、土壤水分及二者的交互作用的響應(yīng)均不顯著, 可能是因為真菌具有細(xì)胞壁結(jié)構(gòu)這一生理結(jié)構(gòu)的優(yōu)勢, 在短期內(nèi)使其更容易抵抗外界環(huán)境因素的脅迫[43–44]。Zhou等的研究結(jié)果表明, 在300 kg·N·ha–1·a–1時Ascomycota的相對豐度顯著降低, Basidiomycota的相對豐度顯著升高[33]。這與本實驗研究結(jié)果不一致, 本研究發(fā)現(xiàn), Ascomycota與土壤NH4+-N、NO3–-N含量呈正相關(guān),則相反(圖5b), 說明實驗的施氮量及實驗處理時間的長短也可能是本實驗土壤真菌相對豐度無顯著變化的原因之一。綜上所述, 土壤微生物多樣性和群落組成與土壤氮水條件變化之間的關(guān)系較為復(fù)雜, 研究需綜合考慮土壤類型、地上部植物、實驗條件對土壤氮水交互效應(yīng)的影響。

    4 結(jié)論

    (1)氮添加對土壤細(xì)菌和真菌多樣性無顯著影響, 土壤水分的變化顯著影響土壤細(xì)菌多樣性, 表現(xiàn)為高水分處理>中水分處理>低水分處理, 且氮添加增強了低水分對細(xì)菌多樣性的抑制作用。氮、水及二者的交互處理對土壤真菌多樣性無顯著影響

    (2)不同類群細(xì)菌對氮添加、土壤水分變化及二者交互作用響應(yīng)不同, 氮添加顯著降低Acidobacteria的相對豐度, 增加Actinobacteria和Bacteroidetes的相對豐度; Acidobacteria、Bacteroidetes的相對豐度隨土壤含水量的增加增加, Actinobacteria、Firmicutes和Gemmatimonadetes則相反; 干旱使氮添加對Acidobacteria相對豐度的抑制作用達到顯著水平; 土壤水分的增加改變了Cyanobacteria的相對豐度對氮添加的響應(yīng)。

    (3)氮水單獨或交互作用對土壤優(yōu)勢真菌相對豐度均無顯著影響。但干旱有增加Ascomycota和降低Basidiomycota相對豐度的趨勢。

    [1] ZHAO Y H, ZHANG L, CHEN Y F, et al. Atmospheric nitrogen deposition to China: A model analysis on nitrogen budget and critical load exceedance[J]. Atmospheric Environment, 2017, 153: 32–40.

    [2] 遆超普, 顏曉元. 基于氮排放數(shù)據(jù)的中國大陸大氣氮素濕沉降量估算[J]. 農(nóng)業(yè)環(huán)境科學(xué)學(xué)報, 2010, 29(8): 1606–1611.

    [3] 朱義族, 李雅穎, 韓繼剛, 等. 水分條件變化對土壤微生物的影響及其響應(yīng)機制研究進展[J]. 應(yīng)用生態(tài)學(xué)報, 2019, 30(12): 4323–4332.

    [4] 吳華清, 陳小梅, 林媚珍, 等. 降水處理對南亞熱帶季風(fēng)林土壤微生物群落結(jié)構(gòu)的影響[J]. 生態(tài)環(huán)境學(xué)報, 2016, 25(4): 583–590.

    [5] 劉紅梅, 張愛林, 皇甫超河, 等. 氮沉降增加對貝加爾針茅草原土壤微生物群落結(jié)構(gòu)的影響[J]. 生態(tài)環(huán)境學(xué)報, 2017, 26(7): 1100–1106.

    [6] CAMPBELL B J, POLSON S W, HANSON T E, et al. The effect of nutrient deposition on bacterial communities in Arctic tundra soil[J]. Environmental Microbiology, 2010, 12(7): 1842–1854.

    [7] XU Y G, YU W T, MA Q, et al. Responses of bacterial and archaeal ammonia oxidisers of an acidic luvisols soil to different nitrogen fertilization rates after 9 years[J]. Biology and Fertility of Soils, 2012, 48(7): 827–837.

    [8] WANG C, LIU D, BAI E. Decreasing soil microbial diversity is associated with decreasing microbial biomass under nitrogen addition[J]. Soil Biology and Biochemistry, 2018, 120: 126–133.

    [9] 鄭勇, 賀紀(jì)正. 森林土壤微生物對干旱和氮沉降的響應(yīng)[J]. 應(yīng)用生態(tài)學(xué)報, 2020, 31(7): 2464–2472.

    [10] MOORHEAD D L, SINSABAUGH R L. A theoretical model of litter decay and microbial interaction[J]. Ecological Monographs, 2006, 76(2): 151–174.

    [11] 王光華, 劉俊杰, 于鎮(zhèn)華, 等. 土壤酸桿菌門細(xì)菌生態(tài)學(xué)研究進展[J]. 生物技術(shù)通報, 2016, 32(2): 14–20.

    [12] 郝亞群, 謝麟, 陳岳民, 等. 中亞熱帶地區(qū)氮沉降對杉木幼林土壤細(xì)菌群落多樣性及組成的影響[J]. 應(yīng)用生態(tài)學(xué)報, 2018, 29(1): 53–58.

    [13] RAMIREZ K S, CRAINE J M, FIERER N. Nitrogen fertilization inhibits soil microbial respiration regardless of the form of nitrogen applied[J]. Soil Biology & Biochemistry, 2010, 42(12): 2336–2338.

    [14] WESSEN E, HALLIN S, PHILIPPOT L. Differential responses of bacterial and archaeal groups at high taxonomical ranks to soil management[J]. Soil Biology & Biochemistry, 2010, 42(10): 1759–1765.

    [15] ALSTER C J, GERMAN D P, LU Y, et al. Microbial enzymatic responses to drought and to nitrogen addition in a southern California grassland[J]. Soil Biology and Biochemistry, 2013, 64: 68–79.

    [16] DAN N, DEVIN C D. Drought Stress and Root-Associated Bacterial Communities[J]. Frontiers in plant science, 2017, 8.

    [17] HAICHAR F Z, MAROL C, BERGE O, et al. Plant host habitat and root exudates shape soil bacterial community structure[J]. The ISME Journal: Multidisciplinary Journal of Microbial Ecology, 2008, 2(12): 1221–1230.

    [18] PREECE C, VERBRUGGEN E, LIU L, et al. Effects of past and current drought on the composition and diversity of soil microbial communities[J]. Soil Biology and Biochemistry, 2019, 131: 28–39.

    [19] NICHOLAS B J, LIM H C, BORGLIN S, et al. Pre-exposure to drought increases the resistance of tropical forest soil bacterial communities to extended drought[J]. The ISME Journal: Multidisciplinary Journal of Microbial Ecology, 2013, 7(2): 384–394.

    [20] ACOSTA-MARTINEZ V, COTTON J, GARDNER T, et al. Predominant bacterial and fungal assemblages in agricultural soils during a record drought/heat wave and linkages to enzyme activities of biogeochemical cycling[J]. Applied Soil Ecology, 2014, 84: 69–82.

    [21] 岳丹丹, 韓貝, ULLAH A, 等. 干旱條件下棉花根際真菌多樣性分析[J]. 作物學(xué)報2021, 47(9): 1806–1815.

    [22] CHODAK M, GOLEBIEWSKI M, MORAWSKA- PLOSKONKA J, et al. Soil chemical properties affect the reaction of forest soil bacteria to drought and rewetting stress[J]. Annals of microbiology, 2015, 65(3): 1627–1637.

    [23] ZHAO Q, JIAN S G, NUNAN N, et al. Altered precipitation seasonality impacts the dominant fungal but rare bacterial taxa in subtropical forest soils[J]. Biology and Fertility of Soils, 2017, 53(2): 231–245.

    x

    [24] 謝高地, 張彩霞, 張昌順, 等. 中國生態(tài)系統(tǒng)服務(wù)的價值[J]. 資源科學(xué), 2015, 37(9): 1740–1746.

    [25] DIJKSTRA F A, AUGUSTINE D J, BREWER P, et al. Nitrogen cycling and water pulses in semiarid grasslands: are microbial and plant processes temporally asynchronous?[J]. Oecologia, 2012, 170(3): 799–808.

    [26] 劉光崧. 土壤理化分析與剖面描述[M]. 北京: 中國標(biāo)準(zhǔn)出版社, 1996.

    [27] 郭萍萍, 黃幸然, 吳旺旺, 等. 不同施氮方式和施氮量對馬尾松和木荷幼苗根系土壤細(xì)菌群落的影響[J]. 生態(tài)學(xué)報, 2021, 41(1): 149–161.

    [28] SCHAEFFER S M, EVANS R D. Pulse additions of soil carbon and nitrogen affect soil nitrogen dynamics in an arid Colorado Plateau shrubland[J]. Oecologia, 2005, 145(3): 425–433.

    [29] CHRISTIAN L, MICAH H, ROB K, et al. Pyrosequencing- based assessment of soil pH as a predictor of soil bacterial community structure at the continental scale. [J]. Applied and environmental microbiology, 2009, 75(15): 5111–5120.

    [30] HUESO S, GARCIA C, HERNANDEZ T. Severe drought conditions modify the microbial community structure, size and activity in amended and unamended soils[J]. Soil Biology and Biochemistry, 2012, 50: 167–173.

    [31] ZHANG H, LIU H, ZHAO J, et al. Elevated precipitation modifies the relationship between plant diversity and soil bacterial diversity under nitrogen deposition in Stipa baicalensis steppe[J]. Applied Soil Ecology, 2017, 119: 345–353.

    [32] ZHANG X, LIU W, SCHLOTER M, et al. Response of the Abundance of Key Soil Microbial Nitrogen-Cycling Genes to Multi-Factorial Global Changes[J]. Plos One, 2013, 8(10).

    [33] ZHOU J, JIANG X, ZHOU B, et al. Thirty four years of nitrogen fertilization decreases fungal diversity and alters fungal community composition in black soil in northeast China[J]. Soil Biology and Biochemistry, 2016, 95: 135–143.

    [34] FERY S D, KNORR M, PARRENT J L, et al. Chronic nitrogen enrichment affects the structure and function of the soil microbial community in temperate hardwood and pine forests[J]. Forest Ecology and Management, 2004, 196(1): 159–171.

    [35] ZHOU Z H, WANG C K, LUO Y Q. Meta-analysis of the impacts of global change factors on soil microbial diversity and functionality[J]. Nature Communications, 2020, 11(1): 3027.

    [36] ZHANG H, WANG L, LIU H, et al. Nitrogen deposition combined with elevated precipitation is conducive to maintaining the stability of the soil fungal diversity on the Stipa baicalensis steppe[J]. Soil Biology and Biochemistry, 2018, 117: 135–138.

    [37] BAHRAM M, HILDEBRAND F, FORSLUND S K, et al. Structure and function of the global topsoil microbiome[J]. Nature, 2018, 560: 233–237.

    [38] CEDERLUND H, WESSEN E, ENWALL K, et al. Soil carbon quality and nitrogen fertilization structure bacterial communities with predictable responses of major bacterial phyla[J]. Applied Soil Ecology, 2014, 84: 62–68.

    [39] 畢捷, 張乃莉, 梁宇, 等. 氮水添加對半干旱地區(qū)草原土壤微生物群落結(jié)構(gòu)的交互影響[J]. 廣東農(nóng)業(yè)科學(xué), 2012, 39(19): 166–170.

    [40] ZENG J, LIU X J, SONG L, et al. Nitrogen fertilization directly affects soil bacterial diversity and indirectly affects bacterial community composition[J]. Soil Biology and Biochemistry, 2016, 92: 41–49.

    [41] 杜瓊, 孔維寶, 韓銳, 等. 土壤微藻的種類及其功能研究進展[J]. 生物學(xué)通報, 2015, 50(1): 1–5.

    [42] 張海芳, 劉紅梅, 趙建寧, 等. 模擬氮沉降和降雨變化對貝加爾針茅草原土壤細(xì)菌群落結(jié)構(gòu)的影響[J]. 生態(tài)學(xué)報, 2018, 38(1): 244–253.

    [43] ROMAIN B L, CATHERINE O A, MARY F K. Changing precipitation pattern alters soil microbial community response to wet-up under a Mediterranean-type climate[J]. The ISME journal, 2015, 9(4): 946–957.

    [44] FRANCISKA T V, ROB I G, MARK B, et al. Soil bacterial networks are less stable under drought than fungal networks[J]. Nature communications, 2018, 9(3303).

    Effects of nitrogen addition and moisture change on microbial diversity and community composition in subtropical plantation red soil

    SUN Haozhao1, LIU Guizhen1, ZHU Yulin1, LIN Xueping1, GUO Hao2, QIAN Xiaojie1, FANG Xiong3,*

    1. Fujian Provincial Key Laboratory of Soil Environmental Health and Regulation, College of Resources and Environment, Fujian Agriculture and Forestry University, Fuzhou 350000, China 2. College of Mechanical and Electrical Engineering, Fujian Agriculture and Forestry University, Fuzhou 350000, China 3. College of Land Resources and Environment, Jiangxi Agricultural University, Nanchang 330045, China

    The intensification of nitrogen deposition and the change of global rainfall pattern are important features of global change. The input of exogenous nitrogen and the change of soil water conditions will inevitably affect soil microbial diversity and community composition. Illumina Miseq high-throughput sequencing technology was used to study the effects of nitrogen addition, water change and their interaction on bacterial and fungal diversity and community composition of red soil in a modeling subtropical plantation, as well as their relationships with soil properties. The results showed that: soil bacterial diversity decreased with the decrease of water content, and the inhibition effect was enhanced by the addition of nitrogen. The Shannon index of soil fungi under medium water treatment was significantly higher than that under low water treatment. Nitrogen addition significantly decreased the relative abundance of Acidobacteria, while increased the relative abundance of Actinobacteriaand Bacteroidetes. In the soil, the relative abundance of Acidobacteria and Bacteroidetes increased with the increase of soil water content, while Actinobacteria and Firmicutes were on the contrary. Drought intensified the inhibition effect of nitrogen addition on the relative abundance of Acidobacteria. Nitrogen addition and soil water interaction significantly affected the relative abundance of Cyanobacteria. But nitrogen addition, soil moisture and their interactions had no significant effects on soil fungal community composition. The results of redundancy analysis showed that the changes of soil bacterial community composition were significantly related to soil water content, pH, NO3–-N content and TN, which accounted for 65.6%, 26.1%, 21.5% and 20.7% of the changes of soil bacterial community composition, respectively. These results suggest that nitrogen addition may alter the intensity and direction of the effects of drought on soil bacterial diversity and community composition.

    nitrogen deposition; soil water content; microbial diversity; microbial community composition

    10.14108/j.cnki.1008-8873.2024.01.012

    S714.3

    A

    1008-8873(2024)01-094-11

    2021-08-18;

    2021-11-09

    國家自然科學(xué)基金項目(42267034); 江西省“千人計劃”創(chuàng)新領(lǐng)軍人才長期青年項目(jxsq2023101103); 福建省松林改造提升行動重點科技攻關(guān)項目(2021FKJ01)

    孫浩釗(1997—), 男, 河南許昌人, 碩士, 研究方向為土壤生態(tài)學(xué), E-mail: 18238670521@163.com

    通信作者:方熊, 博士, 副教授, 主要從事土壤生態(tài)學(xué)研究, E-mail: fangxiong_1103@126.com

    孫浩釗, 劉桂珍, 朱玉璘, 等. 氮添加和水分改變對亞熱帶人工林紅壤微生物多樣性及群落組成的影響[J]. 生態(tài)科學(xué), 2024, 43(1): 94–104.

    SUN Haozhao, LIU Guizhen, ZHU Yulin, et al. Effects of nitrogen addition and moisture change on microbial diversity and community composition in subtropical plantation red soil[J]. Ecological Science, 2024, 43(1): 94–104.

    猜你喜歡
    施氮土壤水分含水量
    不同施氮水平對春玉米光合參數(shù)及產(chǎn)量的影響
    結(jié)合Sentinel- 1B和Landsat8數(shù)據(jù)的針葉林葉片含水量反演研究
    森林工程(2018年4期)2018-08-04 03:23:16
    一次冰雹過程成雹機理的數(shù)值模擬
    施氮水平對冬小麥冠層氨揮發(fā)的影響
    西藏高原土壤水分遙感監(jiān)測方法研究
    基于兩相混合流理論PEMFC含水量特性分析
    均勻施氮利于玉米根系生長及產(chǎn)量形成
    不同覆蓋措施對棗園土壤水分和溫度的影響
    液氯中含水量測定注意事項及改進
    中國氯堿(2014年10期)2014-02-28 01:05:00
    植被覆蓋區(qū)土壤水分反演研究——以北京市為例
    极品人妻少妇av视频| 中文字幕av电影在线播放| 国产成人aa在线观看| 丝瓜视频免费看黄片| 亚洲欧美成人精品一区二区| 2022亚洲国产成人精品| 伦理电影免费视频| 国产 精品1| 亚洲国产毛片av蜜桃av| 制服丝袜香蕉在线| 免费不卡的大黄色大毛片视频在线观看| 国产成人a∨麻豆精品| 麻豆乱淫一区二区| 另类亚洲欧美激情| 国产极品天堂在线| 伦理电影免费视频| 国产日韩欧美视频二区| 亚洲第一青青草原| av在线观看视频网站免费| 亚洲av福利一区| 欧美日韩国产mv在线观看视频| 女的被弄到高潮叫床怎么办| 日本黄色日本黄色录像| 精品国产一区二区久久| 久久精品国产综合久久久| 久久久久久人妻| 亚洲,欧美,日韩| 婷婷色麻豆天堂久久| 777米奇影视久久| 一本大道久久a久久精品| 欧美成人午夜免费资源| 亚洲经典国产精华液单| 亚洲国产av新网站| 我要看黄色一级片免费的| 极品少妇高潮喷水抽搐| 欧美激情高清一区二区三区 | 国产熟女午夜一区二区三区| 日本vs欧美在线观看视频| 国产精品亚洲av一区麻豆 | 国产精品人妻久久久影院| 亚洲内射少妇av| 中文字幕最新亚洲高清| 国产av国产精品国产| 久久精品国产综合久久久| 曰老女人黄片| 婷婷色av中文字幕| 飞空精品影院首页| 精品亚洲成国产av| 日韩欧美一区视频在线观看| 日韩熟女老妇一区二区性免费视频| 久久久久久久久久久免费av| 日本黄色日本黄色录像| 久久久久久久久免费视频了| 亚洲综合色惰| 青草久久国产| 免费高清在线观看日韩| 国产av码专区亚洲av| av网站免费在线观看视频| 欧美精品人与动牲交sv欧美| 电影成人av| 成人黄色视频免费在线看| 亚洲综合色惰| 熟女少妇亚洲综合色aaa.| 黄片播放在线免费| 嫩草影院入口| 欧美日韩精品网址| 蜜桃国产av成人99| 亚洲少妇的诱惑av| 午夜福利一区二区在线看| 欧美日本中文国产一区发布| 人体艺术视频欧美日本| 日韩av在线免费看完整版不卡| 美女国产视频在线观看| 一二三四中文在线观看免费高清| 在线观看人妻少妇| 精品福利永久在线观看| 国产亚洲欧美精品永久| 日韩中文字幕欧美一区二区 | 精品国产国语对白av| 亚洲欧美一区二区三区国产| 久久毛片免费看一区二区三区| 日韩一本色道免费dvd| 亚洲欧洲精品一区二区精品久久久 | 免费观看a级毛片全部| 久久国产精品男人的天堂亚洲| 在线观看免费视频网站a站| 欧美成人精品欧美一级黄| 韩国精品一区二区三区| 国产午夜精品一二区理论片| 黄色配什么色好看| 精品国产一区二区三区久久久樱花| 久久免费观看电影| 国产精品秋霞免费鲁丝片| 高清av免费在线| 黄色毛片三级朝国网站| 国产精品99久久99久久久不卡 | 欧美亚洲日本最大视频资源| 欧美变态另类bdsm刘玥| 成人毛片a级毛片在线播放| 少妇人妻 视频| videossex国产| 亚洲欧美色中文字幕在线| 男的添女的下面高潮视频| av一本久久久久| 亚洲国产色片| 久久亚洲国产成人精品v| 日产精品乱码卡一卡2卡三| 日韩成人av中文字幕在线观看| 欧美日韩亚洲高清精品| 女性生殖器流出的白浆| 欧美日韩亚洲国产一区二区在线观看 | 免费人妻精品一区二区三区视频| 国产成人免费无遮挡视频| 国产深夜福利视频在线观看| 老鸭窝网址在线观看| 成人影院久久| 国产在线一区二区三区精| 欧美日韩一区二区视频在线观看视频在线| 伦理电影免费视频| 99久久人妻综合| 欧美日韩一区二区视频在线观看视频在线| 在线观看国产h片| 超碰97精品在线观看| 黄色配什么色好看| 久久久国产一区二区| 麻豆乱淫一区二区| 久久国产精品大桥未久av| 爱豆传媒免费全集在线观看| 曰老女人黄片| 国产一区二区激情短视频 | 国产不卡av网站在线观看| 制服诱惑二区| 春色校园在线视频观看| 国产av一区二区精品久久| 在线观看免费日韩欧美大片| 亚洲国产最新在线播放| av视频免费观看在线观看| 日本av免费视频播放| 丝袜美足系列| 成人午夜精彩视频在线观看| 日韩伦理黄色片| av电影中文网址| 人妻 亚洲 视频| 欧美亚洲 丝袜 人妻 在线| 成人国语在线视频| 久久久久视频综合| 伦精品一区二区三区| 日韩,欧美,国产一区二区三区| 中文字幕人妻丝袜制服| 国产精品无大码| 日本猛色少妇xxxxx猛交久久| 亚洲欧美一区二区三区久久| 久久久久久免费高清国产稀缺| 久久久a久久爽久久v久久| 亚洲精品美女久久av网站| 亚洲精品久久久久久婷婷小说| 成年女人毛片免费观看观看9 | 超碰97精品在线观看| 少妇的逼水好多| 国产精品无大码| 亚洲欧美一区二区三区久久| 尾随美女入室| 日韩一卡2卡3卡4卡2021年| 曰老女人黄片| 人人妻人人澡人人爽人人夜夜| 大码成人一级视频| 久久久久久久大尺度免费视频| 精品国产一区二区三区久久久樱花| 国产精品三级大全| 久热久热在线精品观看| 永久免费av网站大全| 久久综合国产亚洲精品| 波多野结衣av一区二区av| 18在线观看网站| 精品少妇内射三级| 热99久久久久精品小说推荐| 18禁裸乳无遮挡动漫免费视频| 一边摸一边做爽爽视频免费| 丝袜美腿诱惑在线| 人人澡人人妻人| 中文字幕制服av| 看免费成人av毛片| 一级a爱视频在线免费观看| av电影中文网址| 亚洲成人一二三区av| 亚洲精品久久成人aⅴ小说| 女性生殖器流出的白浆| 成人毛片a级毛片在线播放| 免费观看a级毛片全部| 亚洲国产精品一区三区| 一级a爱视频在线免费观看| 美女中出高潮动态图| 中文字幕另类日韩欧美亚洲嫩草| 搡女人真爽免费视频火全软件| 黄色怎么调成土黄色| 国产亚洲一区二区精品| 午夜免费观看性视频| 欧美日韩一区二区视频在线观看视频在线| 人人妻人人澡人人看| av国产精品久久久久影院| 日韩一区二区三区影片| a级片在线免费高清观看视频| 这个男人来自地球电影免费观看 | 大片电影免费在线观看免费| 欧美成人午夜免费资源| 中文欧美无线码| 天天躁日日躁夜夜躁夜夜| 国产午夜精品一二区理论片| 日韩中文字幕视频在线看片| 国产爽快片一区二区三区| 永久免费av网站大全| 啦啦啦视频在线资源免费观看| 久久精品熟女亚洲av麻豆精品| 男人操女人黄网站| 久久精品夜色国产| 久久女婷五月综合色啪小说| 青春草亚洲视频在线观看| 乱人伦中国视频| 一级,二级,三级黄色视频| 好男人视频免费观看在线| 久久毛片免费看一区二区三区| 少妇人妻久久综合中文| 欧美 日韩 精品 国产| 男女下面插进去视频免费观看| 黄色毛片三级朝国网站| 深夜精品福利| www.熟女人妻精品国产| 久久久国产欧美日韩av| 一区二区av电影网| 精品少妇一区二区三区视频日本电影 | 免费看不卡的av| 国产麻豆69| 人妻系列 视频| 国产精品av久久久久免费| 亚洲欧洲精品一区二区精品久久久 | 建设人人有责人人尽责人人享有的| 香蕉丝袜av| 精品少妇久久久久久888优播| 天天影视国产精品| 最黄视频免费看| 欧美日韩av久久| 夫妻性生交免费视频一级片| 日本爱情动作片www.在线观看| 巨乳人妻的诱惑在线观看| 欧美日韩精品成人综合77777| 亚洲欧美中文字幕日韩二区| 叶爱在线成人免费视频播放| 黄色毛片三级朝国网站| 深夜精品福利| 一边亲一边摸免费视频| 18禁观看日本| 青草久久国产| 麻豆乱淫一区二区| 日韩人妻精品一区2区三区| 午夜老司机福利剧场| 熟女av电影| 亚洲综合色网址| 午夜精品国产一区二区电影| 亚洲av中文av极速乱| 国产精品久久久久成人av| 91精品三级在线观看| 午夜久久久在线观看| 女人高潮潮喷娇喘18禁视频| 亚洲美女视频黄频| 久久久久久久大尺度免费视频| 亚洲精品自拍成人| 水蜜桃什么品种好| 十八禁网站网址无遮挡| 又大又黄又爽视频免费| 国产精品女同一区二区软件| 亚洲精品在线美女| 日韩视频在线欧美| 亚洲精品日韩在线中文字幕| 啦啦啦在线观看免费高清www| 一边摸一边做爽爽视频免费| 黄色配什么色好看| 日本vs欧美在线观看视频| 国产精品99久久99久久久不卡 | 久久久久久久大尺度免费视频| 性高湖久久久久久久久免费观看| 一区二区av电影网| 国产精品熟女久久久久浪| 赤兔流量卡办理| 哪个播放器可以免费观看大片| av卡一久久| av在线老鸭窝| 日韩一本色道免费dvd| 看免费av毛片| 高清在线视频一区二区三区| 午夜福利在线观看免费完整高清在| 高清av免费在线| 成人二区视频| 午夜福利乱码中文字幕| a级毛片黄视频| 国产亚洲精品第一综合不卡| 国产在线免费精品| 少妇人妻久久综合中文| 曰老女人黄片| av有码第一页| 桃花免费在线播放| 国产 精品1| 久久精品国产亚洲av涩爱| 欧美日韩精品成人综合77777| av电影中文网址| 成人午夜精彩视频在线观看| 日韩中字成人| 亚洲av免费高清在线观看| 久久婷婷青草| 十分钟在线观看高清视频www| 久久久久久久大尺度免费视频| 91国产中文字幕| 国产黄色视频一区二区在线观看| 久久精品国产亚洲av高清一级| 久久精品国产综合久久久| 国产成人精品福利久久| 国精品久久久久久国模美| 18禁国产床啪视频网站| 国产片内射在线| 日韩视频在线欧美| 亚洲内射少妇av| 国产精品久久久av美女十八| 桃花免费在线播放| 丝袜喷水一区| 大码成人一级视频| 美女中出高潮动态图| 一本色道久久久久久精品综合| 丝瓜视频免费看黄片| 最黄视频免费看| 蜜桃国产av成人99| 美女主播在线视频| 纯流量卡能插随身wifi吗| 91aial.com中文字幕在线观看| 99热全是精品| 久久精品亚洲av国产电影网| 婷婷色av中文字幕| 婷婷色综合大香蕉| 三上悠亚av全集在线观看| 欧美亚洲 丝袜 人妻 在线| 少妇人妻精品综合一区二区| 婷婷色麻豆天堂久久| 国产成人aa在线观看| 天天操日日干夜夜撸| 春色校园在线视频观看| 大话2 男鬼变身卡| 国产老妇伦熟女老妇高清| av.在线天堂| 午夜福利影视在线免费观看| 永久网站在线| 国产精品熟女久久久久浪| 免费观看在线日韩| 日本欧美国产在线视频| 纵有疾风起免费观看全集完整版| 国产精品99久久99久久久不卡 | 国产极品粉嫩免费观看在线| 丝瓜视频免费看黄片| 人人妻人人澡人人看| 精品99又大又爽又粗少妇毛片| av免费观看日本| 美女视频免费永久观看网站| 日本色播在线视频| 亚洲精品国产av蜜桃| 99热国产这里只有精品6| 日韩中文字幕欧美一区二区 | 交换朋友夫妻互换小说| 一级a爱视频在线免费观看| 国产精品 国内视频| 黄色视频在线播放观看不卡| 亚洲av国产av综合av卡| 亚洲三区欧美一区| 亚洲熟女精品中文字幕| 中文字幕亚洲精品专区| 亚洲伊人久久精品综合| 久久久久久久亚洲中文字幕| 国产精品免费大片| 亚洲人成网站在线观看播放| 美女xxoo啪啪120秒动态图| kizo精华| 综合色丁香网| 伦理电影大哥的女人| 久久 成人 亚洲| 国产成人精品福利久久| 久久午夜福利片| 亚洲人成77777在线视频| 丰满少妇做爰视频| 欧美日韩一区二区视频在线观看视频在线| 久久久久网色| 一级,二级,三级黄色视频| 秋霞在线观看毛片| 亚洲熟女精品中文字幕| 黑人巨大精品欧美一区二区蜜桃| av免费在线看不卡| 久久狼人影院| 亚洲精品国产色婷婷电影| 少妇的逼水好多| 制服人妻中文乱码| 亚洲久久久国产精品| 黄色怎么调成土黄色| 九草在线视频观看| 咕卡用的链子| 国产女主播在线喷水免费视频网站| 人妻少妇偷人精品九色| 日产精品乱码卡一卡2卡三| 久久精品久久久久久久性| 日本av免费视频播放| 国产精品蜜桃在线观看| 久久精品国产自在天天线| 自线自在国产av| 日韩制服骚丝袜av| 国产免费又黄又爽又色| 久久ye,这里只有精品| 亚洲欧美成人精品一区二区| 亚洲色图综合在线观看| 久久女婷五月综合色啪小说| 免费日韩欧美在线观看| 婷婷色麻豆天堂久久| 精品国产超薄肉色丝袜足j| 国产 一区精品| 国产精品蜜桃在线观看| 七月丁香在线播放| 美女中出高潮动态图| 一级毛片电影观看| 欧美+日韩+精品| 日日摸夜夜添夜夜爱| 欧美日韩综合久久久久久| 热99国产精品久久久久久7| 一级片免费观看大全| 麻豆精品久久久久久蜜桃| 亚洲天堂av无毛| 黑人巨大精品欧美一区二区蜜桃| xxx大片免费视频| 成人午夜精彩视频在线观看| xxxhd国产人妻xxx| 精品国产露脸久久av麻豆| 五月天丁香电影| 国产黄色免费在线视频| 精品99又大又爽又粗少妇毛片| 亚洲av.av天堂| 天天影视国产精品| 在线观看国产h片| 在线天堂中文资源库| 久久精品国产综合久久久| 国产福利在线免费观看视频| 秋霞在线观看毛片| 最新中文字幕久久久久| 午夜福利视频在线观看免费| 国产精品 欧美亚洲| 国产毛片在线视频| 卡戴珊不雅视频在线播放| 一级毛片黄色毛片免费观看视频| 最近的中文字幕免费完整| 两性夫妻黄色片| 黑人巨大精品欧美一区二区蜜桃| 国产成人精品婷婷| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 国产精品久久久久久精品电影小说| 国产淫语在线视频| 欧美精品一区二区免费开放| 少妇人妻精品综合一区二区| 大话2 男鬼变身卡| 一本大道久久a久久精品| 欧美最新免费一区二区三区| 亚洲精品av麻豆狂野| 国产成人午夜福利电影在线观看| 成人国产av品久久久| 下体分泌物呈黄色| 国产黄频视频在线观看| 在线观看免费日韩欧美大片| 中文字幕亚洲精品专区| 成人国产av品久久久| 一级毛片黄色毛片免费观看视频| 欧美精品亚洲一区二区| 亚洲美女黄色视频免费看| 999久久久国产精品视频| 亚洲一级一片aⅴ在线观看| 亚洲国产看品久久| 在线观看美女被高潮喷水网站| 久久久久久久久久久免费av| 国产精品久久久久成人av| 一级,二级,三级黄色视频| 成人二区视频| 99久久中文字幕三级久久日本| 精品国产一区二区三区四区第35| 久久综合国产亚洲精品| 亚洲人成77777在线视频| 成年人午夜在线观看视频| 最新的欧美精品一区二区| 晚上一个人看的免费电影| 亚洲伊人久久精品综合| 午夜免费观看性视频| 观看美女的网站| 黑丝袜美女国产一区| a 毛片基地| 国产精品偷伦视频观看了| 亚洲美女搞黄在线观看| 久久久国产一区二区| 国产毛片在线视频| 欧美日韩国产mv在线观看视频| 成人毛片a级毛片在线播放| 另类精品久久| 十分钟在线观看高清视频www| 两个人免费观看高清视频| 天美传媒精品一区二区| 你懂的网址亚洲精品在线观看| 一级黄片播放器| 亚洲美女视频黄频| 人妻 亚洲 视频| 高清不卡的av网站| 人成视频在线观看免费观看| 黄网站色视频无遮挡免费观看| 亚洲国产av影院在线观看| 久热这里只有精品99| 日本爱情动作片www.在线观看| 狠狠婷婷综合久久久久久88av| 麻豆乱淫一区二区| 久久久久久伊人网av| 亚洲视频免费观看视频| 9191精品国产免费久久| 最新的欧美精品一区二区| 亚洲精品一区蜜桃| 中文字幕最新亚洲高清| 精品国产一区二区三区久久久樱花| 少妇人妻精品综合一区二区| 日韩制服骚丝袜av| 18+在线观看网站| 久久午夜福利片| 亚洲激情五月婷婷啪啪| 亚洲伊人色综图| 97在线视频观看| 黄色一级大片看看| av女优亚洲男人天堂| 观看av在线不卡| 日韩电影二区| 少妇猛男粗大的猛烈进出视频| 欧美精品人与动牲交sv欧美| 伊人久久大香线蕉亚洲五| 有码 亚洲区| 看十八女毛片水多多多| 侵犯人妻中文字幕一二三四区| 亚洲国产av新网站| 蜜桃在线观看..| 午夜老司机福利剧场| 超色免费av| 精品亚洲乱码少妇综合久久| 国产日韩一区二区三区精品不卡| 亚洲av电影在线进入| 一边摸一边做爽爽视频免费| av免费在线看不卡| 国产精品二区激情视频| 亚洲精品日本国产第一区| av福利片在线| 日韩电影二区| 热re99久久国产66热| 亚洲av中文av极速乱| 男女午夜视频在线观看| 国产一区二区 视频在线| 亚洲第一区二区三区不卡| 婷婷色av中文字幕| 大片电影免费在线观看免费| 女性生殖器流出的白浆| 少妇 在线观看| 桃花免费在线播放| 亚洲久久久国产精品| 女性被躁到高潮视频| 久久久久精品久久久久真实原创| 亚洲人成网站在线观看播放| 性色avwww在线观看| 久久精品亚洲av国产电影网| 巨乳人妻的诱惑在线观看| 日日摸夜夜添夜夜爱| 日韩av在线免费看完整版不卡| 亚洲少妇的诱惑av| 在线观看美女被高潮喷水网站| 最近最新中文字幕大全免费视频 | 国产亚洲一区二区精品| 日韩中文字幕视频在线看片| 国产欧美亚洲国产| 最近的中文字幕免费完整| 国产成人免费观看mmmm| 精品少妇黑人巨大在线播放| 老司机影院成人| 国产国语露脸激情在线看| 蜜桃在线观看..| av不卡在线播放| 99热国产这里只有精品6| 午夜免费男女啪啪视频观看| 性色avwww在线观看| 波多野结衣一区麻豆| 久久这里有精品视频免费| 欧美日韩国产mv在线观看视频| 国产精品无大码| 在线观看免费视频网站a站| av国产久精品久网站免费入址| 国产精品偷伦视频观看了| 国产色婷婷99| 国产精品久久久久久av不卡| 一区二区三区精品91| 男人添女人高潮全过程视频| 国产成人a∨麻豆精品| 在线观看免费日韩欧美大片| 青春草视频在线免费观看| 毛片一级片免费看久久久久| 国产精品成人在线| 一区福利在线观看| 午夜日韩欧美国产| 中文欧美无线码| 日韩一区二区视频免费看| 99九九在线精品视频| av免费观看日本| 国产午夜精品一二区理论片| av不卡在线播放| 精品国产一区二区久久| 高清视频免费观看一区二区| a 毛片基地| 99久久中文字幕三级久久日本| 在线天堂中文资源库|