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

    秸稈還田與氮肥運(yùn)籌對(duì)水稻不同生育期土壤細(xì)菌群落結(jié)構(gòu)的影響

    2021-01-29 10:44:14王娟娟,胡珈瑋,狄霖,劉玲玲,王桂良,錢曉晴,張振華
    關(guān)鍵詞:秸稈還田氮肥水稻

    王娟娟,胡珈瑋,狄霖,劉玲玲,王桂良,錢曉晴,張振華

    摘要:秸稈還田是當(dāng)前資源循環(huán)利用的重要方式。氮肥運(yùn)籌對(duì)于秸稈還田效率及土壤細(xì)菌群落結(jié)構(gòu)的影響較大。本研究采用盆栽試驗(yàn)的方法,探究秸稈還田和不同氮肥運(yùn)籌條件下水稻主要生育期的土壤細(xì)菌群落結(jié)構(gòu)變化。試驗(yàn)設(shè)置了常規(guī)施氮(T1)、秸稈還田+氮肥后移(T2)、秸稈還田+常規(guī)施氮(T3)、秸稈還田+氮肥前移(T4)4個(gè)處理,分別測(cè)定了水稻苗期、分蘗期、成熟期土壤部分理化性質(zhì)及細(xì)菌群落組成,并于成熟期測(cè)產(chǎn)。結(jié)果表明,各處理中優(yōu)勢(shì)細(xì)菌包括Chloroflexi、Proteobacteria、Acidobacteria、Firmicutes、Actinobacteria,總占比>78%,Chloroflexi在水稻苗期和成熟期占比最高,而Proteobacteria在水稻分蘗期最多。秸稈還田+氮肥后移處理提高了水稻苗期土壤細(xì)菌的豐富度,但對(duì)其均勻度影響較小。在水稻分蘗期和成熟期,秸稈還田與否及氮肥運(yùn)籌對(duì)土壤細(xì)菌多樣性影響較小。秸稈還田提高了水稻產(chǎn)量,以秸稈還田+氮肥前移處理的增產(chǎn)效果最大,且該處理增加了水稻分蘗期土壤硝態(tài)氮、速效磷含量。在水稻成熟期,土壤中與產(chǎn)量呈顯著正相關(guān)的細(xì)菌群落包括Acidobacteria和Gemmatinomadetes等。綜上所述,秸稈還田并增加氮肥基施比例可以增加水稻產(chǎn)量,并對(duì)改善水稻生長(zhǎng)后期土壤細(xì)菌群落結(jié)構(gòu)產(chǎn)生積極作用。

    關(guān)鍵詞:秸稈還田;氮肥;土壤;細(xì)菌群落結(jié)構(gòu);水稻

    中圖分類號(hào):S154.3;S511文獻(xiàn)標(biāo)識(shí)碼:A文章編號(hào):1000-4440(2021)06-1460-11

    Effects of straw returning and nitrogen management on soil microbial community structure at different rice growth stages

    WANG Juan-juan1,HU Jia-wei1,DI Lin2,LIU Ling-ling1,WANG Gui-liang1,QIAN Xiao-qing1,3,ZHANG Zhen-hua3,4

    (1.College of Environmental Science and Engineering, Yangzhou University, Yangzhou 225127, China;2.Zhenjiang Agricultural and Rural Bureau, Zhenjiang 212009, China;3.Institute of Soil Health, Yangzhou University, Yangzhou 225127, China;4.Institute of Agricultural Resources and Environment, Jiangsu Academy of Agricultural Sciences, Nanjing 210014, China)

    Abstract:Straw returning is an important way of resource recycling. Nitrogen (N) application has great influence on straw returning efficiency and soil bacterial community structure. In this study, a pot experiment was conducted to explore the changes of soil bacterial community structure at different rice growth stages under the treatments of straw returning combined with N management. There were four treatments in the experiment, including conventional N fertilization (T1), straw returning with postponed N fertilization (T2), straw returning with conventional N fertilization (T3) and straw returning with moved forward N fertilization (T4). Soil physiochemical properties and bacterial community composition were determined at seedling, tillering and mature stage, and yield was measured at harvest stage. The results showed that the dominant bacteria phyla in all treatments were Chloroflexi, Proteobacteria, Acidobacteria, Firmicutes and Actinobacteria (together accounted for over 78% of the total sequences). Chloroflexi represented the most dominant phylum at seedling and mature stages, while Proteobacteria was the most abundant phylum at tillering stage. Straw returning with postponed N fertilization increased the richness of soil bacteria, but had little influence on bacterial evenness at seedling stage. Straw returning and nitrogen management had little influence on soil bacterial diversity at tillering and maturity stages of rice. Straw returning increased rice yield, and the highest yield was obtained in the treatment of straw returning with moved forward N fertilization (T4). Moreover, the contents of soil available phosphorus (P) and nitrate nitrogen were increased at tillering stage of rice in T4 treatment. The relative abundance of bacterial communities, including Acidobacteria and Gemmatinomadetes was positively correlated with the yield at ripen stage. In conclusion, straw returning combined with moved forward N fertilization can increase rice yield and improve soil microbial community structure at late growth stage.

    Key words:straw returning;nitrogen fertilization;soil;bacterial community structure;rice

    農(nóng)作物秸稈是農(nóng)業(yè)生產(chǎn)過程中主要的副產(chǎn)品,數(shù)量巨大,中國(guó)是農(nóng)業(yè)大國(guó),秸稈資源也非常豐富,但傳統(tǒng)的不合理利用方式不僅會(huì)導(dǎo)致資源嚴(yán)重浪費(fèi),還會(huì)造成嚴(yán)重的環(huán)境問題。秸稈中富含有機(jī)質(zhì)、氮、磷、鉀以及多種中、微量元素,能較好地平衡土壤有機(jī)碳與土壤養(yǎng)分結(jié)構(gòu),因而被廣泛用于可持續(xù)農(nóng)業(yè)生產(chǎn)中[1],秸稈還田利用已成為當(dāng)前研究的熱點(diǎn)[2]。秸稈腐解后可有效提高土壤養(yǎng)分,同時(shí)也可以改善土壤結(jié)構(gòu)和土壤水熱狀況[3]。土壤性質(zhì)的改變則會(huì)進(jìn)一步影響土壤中微生物群落結(jié)構(gòu)以及相關(guān)功能的表達(dá)。

    土壤微生物是農(nóng)田系統(tǒng)重要的組成部分,對(duì)調(diào)控土壤中各類養(yǎng)分循環(huán)、土壤團(tuán)粒結(jié)構(gòu)形成、土壤碳庫動(dòng)態(tài)變化以及作物產(chǎn)量形成等均有深刻影響[4-6]。秸稈還田可能通過直接提供養(yǎng)分或間接改變土壤理化性質(zhì)而影響土壤細(xì)菌組成,而土壤微生物則決定著秸稈腐解過程,并進(jìn)一步影響土壤養(yǎng)分和植物生長(zhǎng)。大量研究結(jié)果表明,秸稈還田后土壤微生物組成及功能均有一定變化[5,7-8]。例如,秸稈還田引起的土壤pH變化,可能對(duì)土壤細(xì)菌群落改變起到關(guān)鍵作用[9]。然而,土壤細(xì)菌群落結(jié)構(gòu)、代謝以及活動(dòng)對(duì)秸稈還田的反饋影響因秸稈、土壤類型、作物品種、輪作方式、肥料運(yùn)籌等不同而異[10-11],如在土壤有機(jī)碳充足的條件下,大部分細(xì)菌的生長(zhǎng)可能受秸稈的影響較小[12]。此外,秸稈還田對(duì)土壤不同營(yíng)養(yǎng)型細(xì)菌群落豐度的影響也有差異[13]。

    秸稈還田對(duì)水稻產(chǎn)量的影響受還田量、耕作方式、水稻品種及肥料運(yùn)籌等多方面因素影響[14]。不合理的秸稈施用不僅不會(huì)提高水稻產(chǎn)量,還可能帶來負(fù)面效應(yīng)[15]。過多的秸稈施用可能會(huì)影響作物根系生長(zhǎng),降低土壤養(yǎng)分,如氮素水平,甚至可能引起雜草、病蟲為害[16],其主要原因在于秸稈施用后土壤的碳氮比并不一定能滿足水稻植株和土壤微生物對(duì)該類營(yíng)養(yǎng)的需要[17]。秸稈還田帶來的大量有機(jī)碳促進(jìn)土壤微生物繁殖,而要求更多的氮源,導(dǎo)致土壤微生物與植物的爭(zhēng)氮現(xiàn)象[18],尤其是秸稈剛腐解時(shí)需要消耗大量的氮營(yíng)養(yǎng),因此,秸稈還田時(shí)通常會(huì)配合施用一定量的氮肥,用以減輕微生物固持無機(jī)氮對(duì)水稻生長(zhǎng)的不良影響[19]。有研究結(jié)果表明,秸稈還田配施氮肥可以有效改善土壤生態(tài)環(huán)境,從而增加作物產(chǎn)量[20]。

    水稻氮肥施用包括基肥和追肥,氮肥的運(yùn)籌方式對(duì)水稻生長(zhǎng)及肥料利用效率至關(guān)重要,同時(shí),施肥方式也影響到農(nóng)業(yè)生產(chǎn)效率[21-22]。有研究發(fā)現(xiàn),氮肥的合理運(yùn)籌可直接決定秸稈和氮肥的利用率,以及作物產(chǎn)量、品質(zhì)[23]。氮肥對(duì)土壤細(xì)菌群落結(jié)構(gòu)組成、優(yōu)勢(shì)細(xì)菌群落相對(duì)豐度的影響甚至要大于秸稈的作用[24]。秸稈還田條件下,氮肥的不同運(yùn)籌方式對(duì)土壤細(xì)菌群落結(jié)構(gòu)的影響較為復(fù)雜。本研究擬探究秸稈還田配合氮肥不同運(yùn)籌方式對(duì)水稻產(chǎn)量以及不同水稻生育時(shí)期土壤細(xì)菌群落組成的影響,明確秸稈還田條件下氮肥運(yùn)籌是否在一定程度上通過影響土壤微生物而最終影響水稻產(chǎn)量,以期為水稻生產(chǎn)中秸稈還田與否以及氮肥合理運(yùn)籌方式的選擇提供科學(xué)依據(jù)和技術(shù)支撐。

    1材料與方法

    1.1土壤采集

    盆栽試驗(yàn)所用土壤采自江蘇省鎮(zhèn)江市丹徒區(qū)榮炳鎮(zhèn)(31°53′36.70″N,119°24′2.81″E),該區(qū)域長(zhǎng)期用作鎮(zhèn)江市農(nóng)業(yè)農(nóng)村局現(xiàn)代農(nóng)業(yè)(稻麥)科技綜合示范項(xiàng)目。土壤為潴育型水稻土,質(zhì)地為黏質(zhì)壤土,試驗(yàn)所用土樣取自稻麥輪作的小麥成熟期。采集0~20 cm耕層土壤,隨機(jī)選擇10個(gè)樣點(diǎn),樣點(diǎn)間距離不少于5 m,每點(diǎn)取長(zhǎng)、寬、高分別為40 cm、40 cm、20 cm的土塊,破碎混勻,自然風(fēng)干后過孔徑30 mm篩,以用于盆栽試驗(yàn)。另從上述混合土樣中隨機(jī)取20點(diǎn),每點(diǎn)約50 g土,磨碎過20目篩和100目篩,測(cè)定土壤的基本理化性質(zhì)(表1)。

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

    秸稈還田盆栽試驗(yàn)于2019年6月在揚(yáng)州大學(xué)環(huán)境科學(xué)與工程學(xué)院的試驗(yàn)大棚中進(jìn)行。試驗(yàn)共設(shè)4個(gè)處理:T1處理,常規(guī)施氮(基蘗肥、穗肥、粒肥的尿素施用量分別為1盆0.87 g、1.01 g、1.01 g);T2處理,秸稈還田+氮肥后移(秸稈的施用量為1 kg土壤2 g秸稈,基蘗肥、穗肥、粒肥的尿素施用量分別為1盆0 g、1.45 g、1.45 g);T3處理,秸稈還田+常規(guī)施氮(秸稈的施用量為1 kg土壤2 g秸稈,基蘗肥、穗肥、粒肥的尿素施用量分別為1盆0.87 g、1.01 g、1.01 g);T4處理,秸稈還田+氮肥前移(秸稈的施用量為1 kg土壤2 g秸稈,基蘗肥、穗肥、粒肥的尿素施用量為1盆1.45 g、0.73 g、0.73 g)。秸稈為小麥秸稈,其總C、N、P2O5、K2O含量分別為460.00 g/kg、3.17 g/kg、0.38 g/kg和14.60 g/kg。

    盆栽使用圓柱形塑料盆缽,直徑為20 cm,高度為30 cm。每缽用土總量為10 kg。裝盆前,按各處理方案將風(fēng)干過篩后的土壤與不同用量的秸稈及氮肥(尿素)混勻。除此之外,每盆施用的其他營(yíng)養(yǎng)元素肥料包括磷肥0.39 g(P2O5)、鉀肥0.52 g(K2O)、硫酸鎂1.30 g、硫酸鋅0.13 g、硫酸錳0.13 g、氯化鈣1.30 g、三氯化鐵0.67 g、硼酸0.07 g、鉬酸銨0.01 g。采用逐步擴(kuò)大法將上述成分與土壤混均后裝入塑料盆中,淹水平衡3 d后移栽水稻幼苗。水稻品種為南粳9108,選擇長(zhǎng)勢(shì)一致的三葉期水稻幼苗移栽至各處理土壤中,每盆3穴,每穴4株。每個(gè)處理3盆,隨機(jī)擺放,且每隔21 d調(diào)整一次盆缽位置。水稻生長(zhǎng)期視蒸發(fā)情況及時(shí)補(bǔ)充水分,觀測(cè)水稻生長(zhǎng)情況,并于收獲期測(cè)定產(chǎn)量。

    分別于水稻苗期(移栽后7 d)、分蘗期(移栽后30 d)及成熟期(移栽后112 d)采集土壤樣品,用于分析測(cè)定。采樣選用直徑為1.5 cm的自制土鉆采集各處理表層(0~20.0 cm)土壤樣品,每盆隨機(jī)取3個(gè)樣點(diǎn)的土壤樣品進(jìn)行混勻?;靹蚝蟮耐寥罉悠贩殖?份,一份自然風(fēng)干用于測(cè)定樣品理化性質(zhì),另一份于-80 ℃冰箱中冷凍用于土壤微生物高通量測(cè)序分析。

    1.3土壤養(yǎng)分測(cè)定

    不同處理下水稻苗期、分蘗期、成熟期的土壤樣品理化性質(zhì)測(cè)定參照文獻(xiàn)[25]的方法。其中,土壤電導(dǎo)率和pH采用純水浸提法(水∶土為5∶1,體積質(zhì)量比)測(cè)定;土壤速效磷含量采用鉬藍(lán)比色分光光度法(0.5 mol/L NaHCO3浸提)測(cè)定;土壤速效鉀含量采用火焰光度法(NH4OAc浸提)測(cè)定;土壤有機(jī)質(zhì)含量采用重鉻酸鉀外加熱法測(cè)定;土壤全氮含量采用半微量開氏法測(cè)定;土壤銨態(tài)氮、硝態(tài)氮含量分別采用KCl浸提靛酚藍(lán)比色法、紫外分光光度法測(cè)定;土壤堿解氮含量采用擴(kuò)散法測(cè)定。

    1.4細(xì)菌多樣性測(cè)定與數(shù)據(jù)分析

    不同處理土壤樣品的細(xì)菌DNA提取使用PowerSoil DNA提取試劑盒(Qiagen公司產(chǎn)品)操作,詳細(xì)步驟參照說明書。所用引物為通用引物338F(5′-ACTCCTACGGGAGGCAGCAG-3′)和806R(5′-GGACTACHVGGGTWTCTAAT-3′),擴(kuò)增細(xì)菌16S rRNA基因的V3~V4區(qū)域。所得PCR擴(kuò)增產(chǎn)物的回收、純化采用AxyPrepDNA Gel Extraction試劑盒(Axygen公司產(chǎn)品)進(jìn)行,并用QuantiFluorTM-ST試劑盒(Promega公司產(chǎn)品)依據(jù)公司提供的標(biāo)準(zhǔn)程序?qū)兓腜CR產(chǎn)物進(jìn)行定量,最后于Illumina MiSeq平臺(tái)進(jìn)行高通量測(cè)序,測(cè)序工作由上海美吉生物醫(yī)藥科技有限公司完成。

    高通量數(shù)據(jù)的生物信息學(xué)分析在上海美吉生物醫(yī)藥科技有限公司提供的生物云平臺(tái)(http://www.i-sanger.com)上進(jìn)行。根據(jù)97%的序列相似度,將測(cè)定結(jié)果中有效序列聚類為不同的分類操作單元(OTUs)。采用RDP classifier貝葉斯算法對(duì)97%相似水平的OTU代表序列進(jìn)行分類,參照細(xì)菌16S rRNA對(duì)比庫Silva(http://www.arb-silva.de)對(duì)各OTU進(jìn)行物種注釋,并獲得各樣本的細(xì)菌分類學(xué)信息,以及不同分類水平上的群落物種組成。按照各樣品細(xì)菌物種注釋及豐富度繪制群落柱狀圖,并計(jì)算細(xì)菌多樣性指數(shù)。采用Student’s t檢驗(yàn)方法比較各土壤樣品的細(xì)菌群落多樣性及物種相對(duì)豐度的差異。采用冗余分析(RDA)對(duì)主要環(huán)境因子和樣品群落組成之間的相關(guān)性進(jìn)行分析。最后采用相關(guān)性Heatmap圖分析各土壤樣品中不同物種與環(huán)境變量以及產(chǎn)量各因素之間的關(guān)系。

    1.5統(tǒng)計(jì)分析

    土壤理化性質(zhì)和水稻產(chǎn)量等的統(tǒng)計(jì)分析采用SPSS 26軟件的單因素方差(ANOVA)分析,差異性比較采用最小顯著性差異(LSD)法。

    2結(jié)果與分析

    2.1不同處理下水稻不同生育期的土壤性質(zhì)及水稻產(chǎn)量

    土壤基本理化性質(zhì)和水稻產(chǎn)量的分析結(jié)果(表2)顯示,土壤pH在苗期以T4處理為最低,在分蘗期,秸稈還田3個(gè)處理(T2、T3、T4)的pH均顯著低于常規(guī)施氮(T1)處理,秸稈還田配合基施氮肥處理(T3、T4處理)則使土壤pH進(jìn)一步下降。秸稈還田+氮肥前移(T4)處理下,水稻苗期、分蘗期的土壤硝態(tài)氮含量和速效磷含量在4個(gè)處理中均最高。在水稻分蘗期,T4處理下的土壤硝態(tài)氮含量顯著高于T1處理、T2處理。此外,各處理間土壤速效鉀含量和有機(jī)質(zhì)含量的差異未達(dá)到顯著水平。在水稻成熟期,各處理的土壤養(yǎng)分含量之間差異不顯著。秸稈還田總體上增加了水稻產(chǎn)量,T3、T4處理下水稻產(chǎn)量分別為1盆54.2 g和58.4 g,均顯著高于T1處理(1盆35.5 g)。

    2.2不同處理下土壤細(xì)菌群落動(dòng)態(tài)變化

    2.2.1細(xì)菌多樣性本試驗(yàn)共采集水稻苗期、分蘗期、成熟期不同處理的盆栽土壤樣品36個(gè),高通量測(cè)序獲得通過質(zhì)量控制的優(yōu)質(zhì)16S rRNA基因序列1 361 200個(gè),平均每個(gè)土壤樣品獲得約37 811個(gè)優(yōu)質(zhì)序列,其平均長(zhǎng)度為463個(gè)堿基對(duì)。按97%的相似度對(duì)序列進(jìn)行聚類分析,共獲得4 435個(gè)OTU。Chao指數(shù)表示一個(gè)樣品中物種的豐富度,而Shannon指數(shù)則表示物種的均勻度。對(duì)各土壤樣品的細(xì)菌多樣性進(jìn)行計(jì)算,結(jié)果(表3)顯示,水稻苗期T2處理的土壤細(xì)菌豐富度顯著高于T1處理和T3處理,分別增加了8.3%和4.6%。然而在分蘗期和成熟期,各處理土樣的細(xì)菌豐富度之間差異均不顯著。土壤細(xì)菌均勻度的分析結(jié)果表明,在水稻苗期,秸稈還田處理的細(xì)菌Shannon指數(shù)與T1處理相比均有所增加;在水稻分蘗期,秸稈還田處理的細(xì)菌Shannon指數(shù)與T1處理相比均有所下降,但各處理之間的差異沒有達(dá)到顯著水平。

    2.2.2細(xì)菌群落組成及變化對(duì)所獲得的OTU進(jìn)行分類,共得到細(xì)菌種類49個(gè)門,124個(gè)綱,233個(gè)目,404個(gè)科,735個(gè)屬。圖1顯示,相對(duì)豐度(即序列占比)較高的優(yōu)勢(shì)細(xì)菌門有Chloroflexi(22.5%)、Proteobacteria(22.1%)、Acidobacteria(17.8%)、 Firmicutes(9.4%)和Actinobacteria(6.5%)。其中,在水稻生長(zhǎng)苗期、成熟期,平均豐度最高的為Chloroflexi;在分蘗期,土壤中細(xì)菌則以Proteobacteria最多。

    同一生育期不同處理下土壤細(xì)菌群落的組間差異顯著性測(cè)定結(jié)果(圖1)表明,在水稻苗期,T4處理下土壤中Bacteroidetes的相對(duì)豐度顯著高于其他處理(P<0.05),而T3處理下Saccharibacteria的相對(duì)豐度則顯著低于其他處理(P<0.001)。在水稻分蘗期,T2處理下Firmicutes和Actinobacteria的相對(duì)豐度分別顯著高于T4處理和T1處理。在水稻成熟期,T1處理下Armatimonadetes的相對(duì)豐度顯著高于其他處理(P<0.001),T2處理下Nitrospirae的相對(duì)豐度顯著高于T1處理。

    秸稈還田后同一處理下土壤細(xì)菌門水平的組成在水稻苗期、分蘗期、成熟期有較大變化。其中,T2處理下土壤中Nitrospirae的相對(duì)豐度在分蘗期低于苗期、成熟期(P<0.05);T2處理下土壤中Saccharibacteria的相對(duì)豐度在分蘗期低于苗期、成熟期(P<0.01)。T3處理下土壤中Firmicutes、Actinobacteria的相對(duì)豐度在苗期顯著高于分蘗期和成熟期(P<0.05)。T4處理下水稻苗期土壤中Chloroflexi的相對(duì)豐度顯著低于成熟期(P<0.05),而水稻苗期土壤中Firmicutes的相對(duì)豐度則顯著高于其他2個(gè)生育時(shí)期(P<0.05)。T1處理下土壤細(xì)菌在門水平上水稻不同生育期間沒有顯著變化。

    在屬水平上,所有處理中相對(duì)豐度最高的5個(gè)優(yōu)勢(shì)細(xì)菌屬為norank_c__SBR2076、norank_c__Acidobacteria、norank_o__Subgroup_7、Anaeromyxobacter和norank_f__Anaerolineaceae(圖2)??傮w而言,屬于Chloroflexi的未分類的norank_c__SBR2076在水稻不同生育期的平均占比均為最高,占5.4%~8.0%,其他細(xì)菌屬的相對(duì)豐度隨水稻生育期變化而波動(dòng)較大。在水稻苗期,T2處理下norank_p__Saccharibacteria的相對(duì)豐度高于T3處理(P<0.001),其他處理各細(xì)菌屬相對(duì)豐度差異不顯著(圖2A)。在水稻分蘗期,T4處理下Bryobacter的相對(duì)豐度高于T2處理(P<0.05),其他處理各細(xì)菌屬相對(duì)豐度差異均不顯著(圖2B)。在水稻成熟期,T1處理下norank_f__Anaerolineaceae的相對(duì)豐度高于T2、T3、T4處理(P<0.001),而Nitrospira、Bryobacter和norank_f__Acidobacteriaceae__Subgroup_1這3個(gè)細(xì)菌屬的相對(duì)豐度則低于T2、T3、T4處理(P<0.05)(圖2C)。

    2.3不同處理下土壤中細(xì)菌群落組成與環(huán)境因子的相關(guān)性

    冗余分析結(jié)果(圖3)表明,總體而言,在水稻苗期、分蘗期和成熟期,秸稈還田處理的土壤細(xì)菌群落結(jié)構(gòu)與秸稈不還田處理差異均較小。氮肥調(diào)控對(duì)各土壤樣品的細(xì)菌群落組成有一些影響,具體表現(xiàn)為,在水稻苗期,T2處理(秸稈還田+氮肥后移)與其他處理土壤微生物群落結(jié)構(gòu)能區(qū)分開來;在水稻分蘗期,T2處理與T4處理,即秸稈還田+氮肥后移處理與秸稈還田+氮肥前移處理土壤微生物群落結(jié)構(gòu)之間可以區(qū)分。秸稈不還田處理的土壤微生物群落結(jié)構(gòu)在水稻成熟期相對(duì)于苗期更集中一些。在水稻苗期,土壤銨態(tài)氮含量對(duì)細(xì)菌群落組成影響顯著(P<0.05),影響最大的處理為T1和T4;在水稻分蘗期,土壤電導(dǎo)率對(duì)土壤細(xì)菌群落組成的影響達(dá)到顯著水平(P<0.05),影響最大的是T3處理;在水稻成熟期,土壤銨態(tài)氮含量對(duì)細(xì)菌群落組成的影響達(dá)到顯著水平(P<0.05),其他土壤理化因子則對(duì)細(xì)菌群落組成沒有顯著影響。

    各處理下土壤細(xì)菌類群相對(duì)豐度與環(huán)境因子之間的相關(guān)性可采用Spearman秩相關(guān)來進(jìn)行分析,表4只列舉了與環(huán)境因子顯著相關(guān)的細(xì)菌門或細(xì)菌屬。表4顯示,土壤中有2大類細(xì)菌群容易受到環(huán)境因子影響,Actinobacteria、Firmicultes、Saccharibacteria與速效磷含量、速效鉀含量、土壤電導(dǎo)率、硝態(tài)氮含量均呈顯著正相關(guān),與pH均呈極顯著負(fù)相關(guān);而Nitrospirae、Chloroflexi和Latescibacteria則與速效鉀含量、土壤電導(dǎo)率和銨態(tài)氮含量均呈顯著負(fù)相關(guān)。土壤有機(jī)質(zhì)含量與細(xì)菌門水平上各細(xì)菌類群之間沒有顯著相關(guān)性。

    各土壤細(xì)菌屬相對(duì)豐度與環(huán)境因子的相關(guān)性分析結(jié)果(表4)表明,相對(duì)豐度排名前15的細(xì)菌屬中只有norank_c__SBR2076、norank_f__Anaerolineaceae、Nitrospira、Pseudarthrobacter、Massilia等受環(huán)境因子影響顯著。其中,Pseudarthrobacter的相對(duì)豐度與各環(huán)境因子主要呈正相關(guān),而其他則主要為負(fù)相關(guān)。

    2.4水稻產(chǎn)量及其構(gòu)成因素與各細(xì)菌類群的關(guān)系

    對(duì)水稻成熟期各處理下土壤不同細(xì)菌門相對(duì)豐度與產(chǎn)量及其構(gòu)成因素的相關(guān)性進(jìn)行分析,表5只列出具有顯著相關(guān)性的細(xì)菌類群,Acidobacteria、Chlorflexi、Gemmatimonadetes的相對(duì)豐度與水稻每穗粒數(shù)呈顯著正相關(guān),Acidobacteria和Gemmatimonadetes的相對(duì)豐度與產(chǎn)量呈正相關(guān),而Nitrospirae的相對(duì)豐度則與千粒質(zhì)量呈顯著正相關(guān)。在屬水平上,相對(duì)豐度與每穗粒數(shù)呈正相關(guān)的細(xì)菌屬最多,包括Gemmatimonas、norank_c_Acidobacteria和4個(gè)未知屬(norank_c__SBR2076、norank_o__Subgroup_7、norank_c__KD4-96和norank_o__SC-I-84)。水稻產(chǎn)量與Gemmatimonas、norank_c_Acidobacteria、norank_c_SBR2076和nonrank_o_Subgroup_7的相對(duì)豐度呈正相關(guān)。其中,主要細(xì)菌屬norank_c_SBR2076、norank_c_Acidobacteria和nonrank_o_Subgroup_7的相對(duì)豐度均與水稻產(chǎn)量及每穗粒數(shù)呈顯著正相關(guān)。

    3討論

    秸稈還田是提高土壤養(yǎng)分、減少資源浪費(fèi)最有效的途徑之一。有研究結(jié)果表明,秸稈還田能增加土壤中有機(jī)質(zhì)含量,并且在一定程度上提高土壤中有效養(yǎng)分(如氮、磷、鉀等)的含量[26-27]。本研究中,秸稈還田對(duì)苗期、分蘗期、成熟期土壤養(yǎng)分含量的影響比較復(fù)雜。就氮素而言,秸稈還田配合常規(guī)施氮(T3)處理下土壤硝態(tài)氮含量在苗期顯著低于常規(guī)施氮(T1)處理,在分蘗期顯著高于T1處理,成熟期則差異不顯著,這可能與水稻生長(zhǎng)對(duì)有效養(yǎng)分的吸收利用有關(guān)。此外,秸稈還田還可能引起土壤中氮的反硝化損失[28]。在秸稈還田配合不同氮肥運(yùn)籌的幾個(gè)處理中,氮肥前移處理提高了水稻苗期和分蘗期的土壤硝態(tài)氮含量,表明增加氮肥基施有利于水稻生育早期土壤養(yǎng)分的供應(yīng)和積累。與常規(guī)施氮處理相比,秸稈還田配施氮肥處理降低了水稻分蘗期的土壤pH。前人研究發(fā)現(xiàn)秸稈分解過程可能產(chǎn)生大量有機(jī)酸和酰胺態(tài)氮或銨態(tài)氮,進(jìn)入土壤后發(fā)生硝化作用釋放氫離子,導(dǎo)致土壤酸化[28-29]。

    土壤微生物多樣性對(duì)維持土壤肥力與生態(tài)系統(tǒng)平衡起著重要作用,但秸稈還田處理對(duì)不同土壤微生物的影響各異。Yan等[30]發(fā)現(xiàn),秸稈還田可以提高土壤細(xì)菌豐富度,但對(duì)其均勻度沒有影響。同樣,本試驗(yàn)中水稻苗期各秸稈還田處理的土壤細(xì)菌豐富度與T1處理相比有所提高,秸稈還田+氮肥后移處理的細(xì)菌豐富度最高,這可能與初期的秸稈投入導(dǎo)致一些異養(yǎng)的細(xì)菌大量生長(zhǎng)有關(guān),而氮肥添加可能更有利于優(yōu)勢(shì)菌的生長(zhǎng),這種差異在其他生育期有所降低。Zhao等[29]研究發(fā)現(xiàn),秸稈還田降低了土壤細(xì)菌OTU、豐度和均勻度,可能是因?yàn)橥寥兰?xì)菌群落受土壤原來的基本理化性質(zhì)影響更大。秸稈還田+氮肥前移處理的細(xì)菌豐富度和均勻度在成熟期均高于其他處理,在一定程度上說明氮肥前移可能更有利于水稻生長(zhǎng)后期土壤細(xì)菌多樣性的提高。

    秸稈還田后,土壤主要細(xì)菌類群組成有所變化,并且在水稻不同生育期表現(xiàn)不同。在所有處理中,Chloroflexi、 Proteobacteria、Acidobacteria、Firmicutes和Actinobacteria為主要優(yōu)勢(shì)菌門。其中,Chloroflexi、Proteobacteria和Firmicutes均為典型的富營(yíng)養(yǎng)細(xì)菌類群,而Acidobacteria則為貧營(yíng)養(yǎng)細(xì)菌類群。秸稈腐解需要各種細(xì)菌類群的參與,Chloroflexi、Proteobacteria和Firmicutes等直接參與秸稈的降解,而放線菌和真菌等則利用降解產(chǎn)物,在后期生長(zhǎng)更快[27]。

    比較各處理間細(xì)菌群落組成的差異,在苗期,T4處理下Bacteroidetes的相對(duì)豐度顯著高于其他處理,表明秸稈還田+氮肥前移處理可以提高該細(xì)菌門的相對(duì)豐度。擬桿菌門是農(nóng)田利用的一個(gè)重要生物學(xué)指標(biāo)[31],前人研究發(fā)現(xiàn)肥沃農(nóng)田的擬桿菌門相對(duì)豐度要高于貧瘠草地,且該細(xì)菌門的相對(duì)豐度與土壤pH呈正相關(guān)[30,32]。本試驗(yàn)同樣證明了在養(yǎng)分供應(yīng)充足的條件下Bacteroidetes的相對(duì)豐度較高,但與前人不同的是,其豐度并不受pH影響。在水稻分蘗期,秸稈還田+氮肥后移處理的Firmicutes和Actinobacteria相對(duì)豐度均高于常規(guī)施氮處理和秸稈還田+氮肥前移處理。厚壁菌門被認(rèn)為與有機(jī)碳降解關(guān)系較大[33],可能T2處理的氮肥后移減緩了秸稈降解速率,到了分蘗期仍有大量可利用于碳源,而T4處理的秸稈因氮肥前移而秸稈在較短時(shí)間內(nèi)得到較充分的降解,因此在水稻分蘗期厚壁菌門相對(duì)豐度較其他細(xì)菌群落低,這與前人研究結(jié)果相近[34]。此外,F(xiàn)irmicutes和Actinobacteria的相對(duì)豐度均與土壤pH呈負(fù)相關(guān),這也能進(jìn)一步解釋兩者在T2處理中相對(duì)豐度遠(yuǎn)高于T1處理。在水稻成熟期,T2處理下Nitrospirae的相對(duì)豐度顯著高于T1處理。硝化螺旋菌門的活動(dòng)與土壤氮素相關(guān),有研究發(fā)現(xiàn)秸稈還田處理會(huì)增加硝化螺旋菌門的相對(duì)豐度[35],同時(shí)配合氮肥后移可能使得該細(xì)菌門的細(xì)菌在水稻生長(zhǎng)后期更加活躍。而在本研究秸稈還田條件下,氮肥運(yùn)籌對(duì)于硝化螺旋菌門的影響較小。

    秸稈還田后,各處理下水稻不同生育期的土壤細(xì)菌群落組成均有一定變化,而不加秸稈的常規(guī)施氮處理下土壤細(xì)菌群落組成在整個(gè)生育期比較穩(wěn)定,可見秸稈還田對(duì)土壤微生物還有季節(jié)性的影響,這與Liu等[27]提出的秸稈降低過程通常伴隨著土壤細(xì)菌組成的季節(jié)性變化的結(jié)論相一致。冗余分析結(jié)果表明,土壤理化因子對(duì)細(xì)菌群落分布的影響在不同生育期表現(xiàn)不同,主要是銨態(tài)氮含量和土壤電導(dǎo)率這2個(gè)因子,銨態(tài)氮含量在苗期、成熟期均與細(xì)菌群落組成顯著相關(guān)。有研究結(jié)果表明,土壤微生物群落受氮肥影響較大,且因環(huán)境條件而異[36]。尿素的添加會(huì)改變土壤中銨態(tài)氮含量,而水稻土壤中與氮循環(huán)有關(guān)的微生物往往占優(yōu)勢(shì)地位[37],一些相關(guān)細(xì)菌則很容易受到影響[38],從而進(jìn)一步影響細(xì)菌群落組成。土壤電導(dǎo)率受秸稈分解產(chǎn)生的一些有機(jī)酸等化合物影響[39],對(duì)水稻分蘗期的細(xì)菌群落組成有顯著影響。

    有研究結(jié)果表明,秸稈還田可以改善土壤養(yǎng)分、水稻產(chǎn)量和品質(zhì)[6]。水稻成熟期,土壤中Acidobacteria和Gemmatimonadetes與產(chǎn)量呈顯著正相關(guān);Chloroflexi、Acidobacteria、Gemmatimonadetes與每穗粒數(shù)呈正顯著或極顯著相關(guān),Nitrospirae與千粒質(zhì)量呈顯著正相關(guān)。

    4結(jié)論

    秸稈還田且氮肥后移處理提高了水稻苗期土壤細(xì)菌的豐富度,但對(duì)其均勻度影響較小。秸稈還田與否以及氮肥運(yùn)籌對(duì)水稻分蘗期、成熟期土壤細(xì)菌豐富度和均勻度影響較小;水稻各生育期,對(duì)細(xì)菌群落結(jié)構(gòu)組成有顯著影響的土壤理化因子不同,苗期、成熟期以銨態(tài)氮影響最大,而分蘗期以土壤電導(dǎo)率的影響最為顯著。相對(duì)豐度(即序列占比)較高的優(yōu)勢(shì)細(xì)菌門有Chloroflexi、Proteobacteria、Acidobacteria、Firmicutes和Actinobacteria。其中,在水稻生長(zhǎng)苗期、成熟期,平均豐度最高的為Chloroflexi;在分蘗期,土壤中細(xì)菌則以Proteobacteria最多;Chloroflexi的未分類屬norank_c__SBR2076在水稻苗期、分蘗期、成熟期的相對(duì)豐度均最高。秸稈還田增加了水稻產(chǎn)量,各處理中以T4處理的增產(chǎn)效果最大。相對(duì)豐度與水稻產(chǎn)量顯著正相關(guān)的細(xì)菌包括Acidobacteria和Gemmatimonadetes,與千粒質(zhì)量正相關(guān)的為Nitrospirae,norank_o_Subgroup_7的相對(duì)豐度與水稻千粒質(zhì)量、每穗粒數(shù)、產(chǎn)量均呈顯著正相關(guān)。因此,秸稈還田并增加氮肥基施比例可以明顯增加水稻產(chǎn)量,并對(duì)改善水稻生長(zhǎng)后期土壤細(xì)菌群落結(jié)構(gòu)產(chǎn)生有益的影響。

    參考文獻(xiàn):

    [1]ZHAO S C, LI K J, ZHOU W, et al. Changes in soil microbial community, enzyme activities and organic matter fractions under long-term straw return in north-central China [J]. Agriculture, Ecosystems and Environment, 2016, 216: 82-88.

    [2]SU Y, KWONG R W M, TANG W L, et al. Straw return enhances the risks of metals in soil? [J]. Ecotoxicology and Environmental Safety, 2021, 207: 111201.

    [3]YANG H S, ZHAI S L, LI Y F, et al. Waterlogging reduction and wheat yield increase through long-term ditch-buried straw return in a rice-wheat rotation system [J]. Field Crops Research, 2017, 209: 189-197.

    [4]王青霞,李美霖,陳喜靖,等. 秸稈還田下氮肥運(yùn)籌對(duì)水稻各生育期土壤微生物群落結(jié)構(gòu)的影響[J]. 應(yīng)用生態(tài)學(xué)報(bào), 2020, 31(3):935-944.

    [5]HARRIS J. Soil microbial communities and restoration ecology: facilitators or followers? [J]. Science, 2009, 325: 573-574.

    [6]CHEN Z, WANG H, LIU X, et al. Changes in soil microbial community and organic carbon fractions under short-term straw return in a rice-wheat cropping system [J]. Soil and Tillage Research, 2017, 165: 121-127.

    [7]DING L J, SU J Q, SUN G X, et al. Increased microbial functional diversity under long-term organic and integrated fertilization in a paddy soil [J]. Applied Microbiology and Biotechnology, 2018, 102(4): 1969-1982.

    [8]BEI S K, ZHANG Y L, Li T T, et al. Response of the soil microbial community to different fertilizer inputs in a wheat-maize rotation on a calcareous soil [J]. Agriculture, Ecosystems and Environment, 2018, 260: 58-69.

    [9]WANG N, YU J G, ZHAO Y H, et al. Straw enhanced CO2 and CH4 but decreased N2O emissions from flooded paddy soils: changes in microbial community compositions [J]. Atmospheric Environment, 2017, 174: 171-179.

    [10]YANG H S, MENG Y, FENG J X, et al. Direct and indirect effects of long-term ditch-buried straw return on soil bacterial community in a rice-wheat rotation system [J]. Land Degradation & Development, 2019, 31: 851-867.

    [11]ZENG X N, SONG Q L, YU S, et al. Effects of fertilization and straw incorporation on bacterial communities in black soil, northeastern China [J]. Applied Ecology and Environmental Research, 2020, 18(6): 7535-7550.

    [12]ZHU J, PENG H, JI X H, et al. Effects of reduced inorganic fertilization and rice straw recovery on soil enzyme activities and bacterial community in double-rice paddy soils [J]. European Journal of Soil Biology, 2019, 94: 103116.

    [13]ZHAO S C, QIU S J, XU X P, et al. Change in straw decomposition rate and soil microbial community composition after straw addition in different long-term fertilization soils [J]. Applied Soil Ecology, 2019, 138: 123-133.

    [14]王子陽. 秸稈還田與氮肥運(yùn)籌對(duì)雙季稻產(chǎn)量、品質(zhì)與土壤養(yǎng)分的影響 [D].揚(yáng)州:揚(yáng)州大學(xué),2021

    [15]HUANG S, ZENG Y J, WU J F, et al. Effect of crop residue retention on rice yield in China: a meta-analysis [J]. Field Crops Research, 2013, 154: 188-194.

    [16]SU Y, YU M, XI H, et al. Soil microbial community shifts with long-term of different straw return in wheat-corn rotation system [J]. Scientific Reports, 2020, 10(1): 6360.

    [17]LI T, GAO J S, BAI L Y, et al. Influence of green manure and rice straw management on soil organic carbon, enzyme activities, and rice yield in red paddy soil [J]. Soil & Tillage Research, 2019, 195: 104428.

    [18]孔麗麗,侯云鵬,尹彩俠,等.秸稈還田下寒地水稻實(shí)現(xiàn)高產(chǎn)高氮肥利用率的氮肥運(yùn)籌模式[J].植物營(yíng)養(yǎng)與肥料學(xué)報(bào),2021,27(7):1282-1293.

    [19]何虎,吳建富,曾研華,等.稻草全量還田下氮肥運(yùn)籌對(duì)雙季晚稻產(chǎn)量及其氮素吸收利用的影響[J].植物營(yíng)養(yǎng)與肥料學(xué)報(bào),2014,20(4):811-820.

    [20]侯賢清,吳鵬年,王艷麗,等. 秸稈還田配施氮肥對(duì)土壤水肥狀況和玉米產(chǎn)量的影響 [J]. 應(yīng)用生態(tài)學(xué)報(bào), 2018, 29(6): 1928-1934.

    [21]鄒露陽,陳震,段福義, 等. 噴灌水肥與種植密度互作對(duì)水氮利用效率和冬小麥產(chǎn)量的影響[J]. 排灌機(jī)械工程學(xué)報(bào), 2020, 38(6):632-636.

    [22]李晶晶,臧文靜,黎耀軍,等.噴灌施氮管理對(duì)春玉米產(chǎn)量及水氮利用的影響[J]. 排灌機(jī)械工程學(xué)報(bào),2020,38(12):1277-1283.

    [23]嚴(yán)奉君,孫永健,馬均,等. 不同土壤肥力條件下麥稈還田與氮肥運(yùn)籌對(duì)雜交稻氮素利用、產(chǎn)量及米質(zhì)的影響[J]. 中國(guó)水稻科學(xué), 2015, 29(1):56-64.

    [24]SUN R B, ZHANG X X, GUO X S, et al. Bacterial diversity in soils subjected to long-term chemical fertilization can be more stably maintained with the addition of livestock manure than wheat straw [J]. Soil Biology and Biochemistry, 2015, 88: 9-18.

    [25]鮑士旦.土壤農(nóng)化分析[M]. 3版. 北京:中國(guó)農(nóng)業(yè)出版社, 2000:30, 44, 50, 53, 56, 81, 106, 185.

    [26]LIU C, LU M, CUI J, et al. Effects of straw carbon input on carbon dynamics in agricultural soils: a meta-analysis [J]. Global Change Biology, 2014, 20(5): 1366-1381.

    [27]LIU L, DING M J, ZHOU L K, et al. Effects of different rice straw on soil microbial community structure [J]. Agronomy Journal, 2021, 113(2): 794-805.

    [28]HU X J, LIU J J, WEI D, et al. Soil bacterial communities under different long-term fertilization regimes in three locations across Effects of different rice straw agronomy journal the black soil region of Northeast China[J]. Pedosphere, 2018, 28(5): 751-763.

    [29]ZHAO S C, QIU S J, XU X P, et al. Change in straw decomposition rate and soil microbial community composition after straw addition in different long-term fertilization soils [J]. Applied Soil Ecology, 2019, 138: 123-133.

    [30]YAN S S, SONG J M, FAN J S, et al. Changes in soil organic carbon fractions and microbial community under rice straw return in Northeast China [J]. Global Ecology and Conservation, 2020, 22:e00962.

    [31]WOLINSKA A, KUZNIAR A, ZIELENKIEWICZ U, et al. Bacteroidetes as a sensitive biological indicator of agricultural soil usage revealed by a culture-independent approach [J]. Applied Soil Ecology, 2017, 119: 128-137.

    [32]LAUBER C L, HAMADY M, KNIGHT R, 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.

    [33]BONANOMI G, FILIPPIS F D, CESARANO G, et al. Organic farming induces changes in soil microbiota that affect agro-ecosystem functions [J]. Soil Biology & Biochemistry, 2016, 103: 327-336.

    [34]GUO H H, NASIR M, LV J L, et al. Understanding the variation of microbial community in heavy metals contaminated soil using high throughput sequencing [J]. Ecotoxicology and Environmental Safety, 2017, 144: 300-306.

    [35]劉驍蒨. 秸稈還田方式與施肥對(duì)水稻土壤微生物學(xué)特性的影響 [D].成都:四川農(nóng)業(yè)大學(xué), 2013.

    [36]CHU H, LIN X, FUJII T, et al. Soil microbial biomass, dehydrogenase activity, bacterial community structure in response to longterm fertilizer management[J]. Soil Biology and Biochemistry, 2007, 39:2971-2976.

    [37]LIGI T, TRUU M, TRUU J, et al. Effects of soil chemical characteristics and water regime on denitrification genes (nirS, nirK, and nosZ) abundances in a created riverine wetland complex[J]. Ecological Engineering, 2014, 72:47-55.

    [38]WANG N, DING L J, XU H J, et al. Variability in responses of bacterial communities and nitrogen oxide emission to urea fertilization among various flooded paddy soils[J]. FEMS Microbiology Ecology, 2015, 91(3):1-11

    [39]SU Y, YU M, XI H, et al. Soil microbial community shifts with long-term of different straw return in wheat-corn rotation system [J]. Scientific Reports, 2020, 10(1): 6360.

    (責(zé)任編輯:王妮)

    收稿日期:2021-08-24

    基金項(xiàng)目:國(guó)家自然科學(xué)青年科學(xué)基金項(xiàng)目(41701093);中國(guó)博士后基金項(xiàng)目(2017M611928);國(guó)家重點(diǎn)研發(fā)計(jì)劃課題項(xiàng)目(2017YFD0200107)

    作者簡(jiǎn)介:王娟娟(1979-),女,江蘇泗洪人,博士,副教授,主要從事環(huán)境微生物學(xué)研究。(E-mail)wangjuanjuan@yzu.edu.cn

    通訊作者:張振華,(E-mail) zhenhuaz70@hotmail.com

    猜你喜歡
    秸稈還田氮肥水稻
    什么是海水稻
    水稻種植60天就能收獲啦
    軍事文摘(2021年22期)2021-11-26 00:43:51
    氮肥供應(yīng)充足 春耕生產(chǎn)有保障
    一季水稻
    文苑(2020年6期)2020-06-22 08:41:52
    水稻花
    文苑(2019年22期)2019-12-07 05:29:00
    江淮小氮肥 耕耘六十年——紀(jì)念安徽小氮肥誕生六十周年
    抓住機(jī)遇 主動(dòng)作為 努力推進(jìn)我國(guó)氮肥市場(chǎng)穩(wěn)步前行
    2017春季各地氮肥市場(chǎng)掠影
    玉米高產(chǎn)栽培秸稈還田保護(hù)性耕作技術(shù)措施
    赤壁市秸稈還田替代鉀肥在中稻上的應(yīng)用效果研究
    亚洲国产精品合色在线| 91av网一区二区| 久久精品影院6| 久久久久久大精品| 舔av片在线| a级毛片a级免费在线| 丰满的人妻完整版| 少妇丰满av| 国产精品无大码| 自拍偷自拍亚洲精品老妇| 动漫黄色视频在线观看| 欧美区成人在线视频| 国产白丝娇喘喷水9色精品| a级一级毛片免费在线观看| 精品人妻一区二区三区麻豆 | 亚洲内射少妇av| 天堂√8在线中文| 琪琪午夜伦伦电影理论片6080| 亚洲精品色激情综合| or卡值多少钱| 欧美黑人巨大hd| 亚洲人成伊人成综合网2020| 亚洲无线观看免费| 欧美+亚洲+日韩+国产| 五月玫瑰六月丁香| 18禁裸乳无遮挡免费网站照片| 婷婷丁香在线五月| 久久久久免费精品人妻一区二区| ponron亚洲| 91午夜精品亚洲一区二区三区 | 99热精品在线国产| 九九爱精品视频在线观看| 一进一出抽搐动态| 午夜福利视频1000在线观看| 日韩强制内射视频| 亚洲18禁久久av| 午夜a级毛片| 美女大奶头视频| 国产三级中文精品| 精品一区二区三区视频在线观看免费| 欧美日韩国产亚洲二区| 亚洲在线自拍视频| 国产精华一区二区三区| 偷拍熟女少妇极品色| eeuss影院久久| 日韩一本色道免费dvd| 99久久九九国产精品国产免费| 最近最新中文字幕大全电影3| 久久6这里有精品| 欧美bdsm另类| 久久这里只有精品中国| 亚州av有码| 成年女人永久免费观看视频| 日日干狠狠操夜夜爽| 国产亚洲av嫩草精品影院| 51国产日韩欧美| 国产v大片淫在线免费观看| 久久人妻av系列| 麻豆一二三区av精品| 久久久久国产精品人妻aⅴ院| 国产极品精品免费视频能看的| 国产一区二区三区av在线 | 欧美一区二区精品小视频在线| 两性午夜刺激爽爽歪歪视频在线观看| 国产精品人妻久久久久久| 免费一级毛片在线播放高清视频| 成人美女网站在线观看视频| 日日摸夜夜添夜夜添小说| 亚洲午夜理论影院| 免费看光身美女| 一个人看视频在线观看www免费| 全区人妻精品视频| 欧美日本亚洲视频在线播放| 色噜噜av男人的天堂激情| 日本 av在线| 亚洲七黄色美女视频| 日本爱情动作片www.在线观看 | 国产一区二区亚洲精品在线观看| 最近视频中文字幕2019在线8| 97热精品久久久久久| 国产日本99.免费观看| 成年免费大片在线观看| 色尼玛亚洲综合影院| 亚洲成a人片在线一区二区| 熟女人妻精品中文字幕| 很黄的视频免费| 一级av片app| 亚洲人成网站高清观看| 少妇的逼水好多| 国产精品乱码一区二三区的特点| ponron亚洲| 国产私拍福利视频在线观看| a级毛片a级免费在线| 亚洲精品影视一区二区三区av| 51国产日韩欧美| 99久久成人亚洲精品观看| 看片在线看免费视频| 亚洲av五月六月丁香网| 久久久久久九九精品二区国产| 久久精品夜夜夜夜夜久久蜜豆| 欧美一区二区国产精品久久精品| 一进一出抽搐动态| 在线观看舔阴道视频| 国产午夜精品久久久久久一区二区三区 | 99久久久亚洲精品蜜臀av| 麻豆国产97在线/欧美| 简卡轻食公司| 亚洲欧美日韩东京热| 色噜噜av男人的天堂激情| 中亚洲国语对白在线视频| 国内久久婷婷六月综合欲色啪| 国产亚洲精品久久久久久毛片| 亚洲av美国av| bbb黄色大片| 免费av观看视频| 亚洲成人精品中文字幕电影| 久久午夜福利片| 九色国产91popny在线| 悠悠久久av| 久久精品综合一区二区三区| 国产v大片淫在线免费观看| 看十八女毛片水多多多| 免费电影在线观看免费观看| 亚洲人与动物交配视频| 欧美人与善性xxx| 黄片wwwwww| 嫩草影院新地址| 九九爱精品视频在线观看| 国产熟女欧美一区二区| 97超视频在线观看视频| 国产美女午夜福利| 欧美人与善性xxx| 国产主播在线观看一区二区| 久久久国产成人免费| 欧美高清成人免费视频www| 精品久久久久久,| 少妇猛男粗大的猛烈进出视频 | 干丝袜人妻中文字幕| 日韩精品中文字幕看吧| 香蕉av资源在线| av女优亚洲男人天堂| 免费一级毛片在线播放高清视频| 午夜a级毛片| 日本一二三区视频观看| 精品久久久久久成人av| 亚洲无线观看免费| 国产男靠女视频免费网站| 天堂动漫精品| 欧美成人a在线观看| 中出人妻视频一区二区| 美女被艹到高潮喷水动态| 日本黄大片高清| 精品一区二区免费观看| 精品欧美国产一区二区三| 波多野结衣巨乳人妻| 日本色播在线视频| 亚洲aⅴ乱码一区二区在线播放| 三级毛片av免费| 少妇的逼水好多| 美女免费视频网站| 最后的刺客免费高清国语| 99精品在免费线老司机午夜| 国产一区二区三区av在线 | 欧美性猛交黑人性爽| 国产精品久久久久久久久免| 国产一级毛片七仙女欲春2| a级毛片免费高清观看在线播放| 免费看日本二区| 亚洲精品一卡2卡三卡4卡5卡| 亚洲国产欧洲综合997久久,| 欧美极品一区二区三区四区| 亚洲色图av天堂| 亚洲专区中文字幕在线| 美女大奶头视频| 黄色女人牲交| 亚洲狠狠婷婷综合久久图片| 午夜福利在线观看吧| 国产日本99.免费观看| 日日夜夜操网爽| 亚洲人成网站在线播放欧美日韩| 美女大奶头视频| 听说在线观看完整版免费高清| 成年人黄色毛片网站| 精品久久久久久久久久久久久| 国产免费男女视频| 波多野结衣高清无吗| 国产精品久久久久久精品电影| 久久精品国产亚洲av香蕉五月| 99久久精品热视频| 亚州av有码| 精品午夜福利在线看| 国产黄a三级三级三级人| 亚洲国产欧美人成| 日本 欧美在线| 干丝袜人妻中文字幕| 午夜视频国产福利| 亚洲成人久久性| 亚洲av.av天堂| 亚洲欧美日韩卡通动漫| 成年免费大片在线观看| 看黄色毛片网站| 1000部很黄的大片| 免费不卡的大黄色大毛片视频在线观看 | 中文资源天堂在线| 日韩人妻高清精品专区| 国产精品爽爽va在线观看网站| 直男gayav资源| 亚洲精品亚洲一区二区| 久久亚洲真实| 内射极品少妇av片p| 亚洲自偷自拍三级| 欧美3d第一页| 久久热精品热| 午夜免费激情av| 中文资源天堂在线| 88av欧美| 亚洲第一电影网av| 日日啪夜夜撸| 麻豆国产av国片精品| 国产精品乱码一区二三区的特点| 淫秽高清视频在线观看| 精品久久久久久久久久久久久| 99国产极品粉嫩在线观看| 99riav亚洲国产免费| 中国美白少妇内射xxxbb| 国产精品一及| 国产69精品久久久久777片| 中文亚洲av片在线观看爽| 床上黄色一级片| 欧洲精品卡2卡3卡4卡5卡区| 日韩欧美精品v在线| 亚洲av免费高清在线观看| 欧美精品啪啪一区二区三区| 听说在线观看完整版免费高清| 观看免费一级毛片| 美女大奶头视频| 成人性生交大片免费视频hd| 国产精品久久久久久久电影| 91久久精品国产一区二区成人| 赤兔流量卡办理| 亚洲av.av天堂| 成人鲁丝片一二三区免费| 国内精品美女久久久久久| 欧美zozozo另类| 国产伦精品一区二区三区四那| 国产蜜桃级精品一区二区三区| 免费在线观看成人毛片| 免费观看在线日韩| 亚洲av熟女| 国产激情偷乱视频一区二区| 亚州av有码| 老司机福利观看| 国产三级中文精品| 禁无遮挡网站| 在线观看免费视频日本深夜| 男女视频在线观看网站免费| 一个人看视频在线观看www免费| 麻豆国产97在线/欧美| 色尼玛亚洲综合影院| 国产精品福利在线免费观看| 亚洲成人久久爱视频| 99视频精品全部免费 在线| 日韩强制内射视频| 校园春色视频在线观看| 欧美日韩综合久久久久久 | 亚洲欧美激情综合另类| 亚洲中文字幕一区二区三区有码在线看| 两个人视频免费观看高清| 国产色婷婷99| 亚洲av.av天堂| 精品日产1卡2卡| 午夜福利在线在线| 欧美日韩乱码在线| 欧洲精品卡2卡3卡4卡5卡区| 国产高清视频在线播放一区| 我要看日韩黄色一级片| 老司机午夜福利在线观看视频| 天天一区二区日本电影三级| 99久久中文字幕三级久久日本| 日本爱情动作片www.在线观看 | 日本a在线网址| 99九九线精品视频在线观看视频| 欧美激情在线99| 夜夜爽天天搞| 一本一本综合久久| 欧美3d第一页| 又黄又爽又刺激的免费视频.| 韩国av一区二区三区四区| 又紧又爽又黄一区二区| 搡老岳熟女国产| 美女xxoo啪啪120秒动态图| 性插视频无遮挡在线免费观看| 久久精品国产亚洲网站| 国产精品久久电影中文字幕| 久久久久久久亚洲中文字幕| 老司机深夜福利视频在线观看| 又黄又爽又免费观看的视频| 我的老师免费观看完整版| 琪琪午夜伦伦电影理论片6080| 91精品国产九色| 熟女电影av网| 十八禁网站免费在线| 欧美丝袜亚洲另类 | 自拍偷自拍亚洲精品老妇| 少妇猛男粗大的猛烈进出视频 | 深夜a级毛片| 在线观看午夜福利视频| 国产av麻豆久久久久久久| 婷婷色综合大香蕉| 久久国内精品自在自线图片| 欧美高清成人免费视频www| 国产不卡一卡二| 嫁个100分男人电影在线观看| 午夜视频国产福利| 一a级毛片在线观看| 男人的好看免费观看在线视频| 99精品在免费线老司机午夜| 国产真实伦视频高清在线观看 | 婷婷精品国产亚洲av| 精品午夜福利在线看| 日本撒尿小便嘘嘘汇集6| 亚洲国产精品合色在线| 亚洲av免费在线观看| 成年免费大片在线观看| 精品不卡国产一区二区三区| 国产精品日韩av在线免费观看| 欧洲精品卡2卡3卡4卡5卡区| 国产男靠女视频免费网站| 久久九九热精品免费| 一个人看视频在线观看www免费| 久久亚洲精品不卡| 男人狂女人下面高潮的视频| 久久国产乱子免费精品| 成年免费大片在线观看| av天堂在线播放| 熟女人妻精品中文字幕| 国产一区二区三区视频了| 欧美成人一区二区免费高清观看| 黄色一级大片看看| 在线观看美女被高潮喷水网站| 99久久精品一区二区三区| 国产蜜桃级精品一区二区三区| bbb黄色大片| 亚洲熟妇熟女久久| 性色avwww在线观看| 联通29元200g的流量卡| 久久久午夜欧美精品| 色播亚洲综合网| 欧美一级a爱片免费观看看| 午夜免费成人在线视频| 夜夜爽天天搞| 赤兔流量卡办理| 国产精品av视频在线免费观看| 麻豆国产av国片精品| 搞女人的毛片| 国产高清视频在线观看网站| 性欧美人与动物交配| 午夜福利在线观看免费完整高清在 | 精品久久久久久成人av| 两个人的视频大全免费| 亚洲av熟女| 日日摸夜夜添夜夜添小说| 男女下面进入的视频免费午夜| 日韩 亚洲 欧美在线| 动漫黄色视频在线观看| 丝袜美腿在线中文| 真人做人爱边吃奶动态| 欧美日韩精品成人综合77777| 真人一进一出gif抽搐免费| 熟妇人妻久久中文字幕3abv| 欧美又色又爽又黄视频| 亚洲最大成人手机在线| 日韩国内少妇激情av| 国产成年人精品一区二区| 国产日本99.免费观看| 中文亚洲av片在线观看爽| 最近中文字幕高清免费大全6 | 熟女电影av网| 久久久久久大精品| 亚洲av电影不卡..在线观看| 一个人看的www免费观看视频| 久久久久久久久大av| 99久久精品国产国产毛片| 免费一级毛片在线播放高清视频| 国产真实伦视频高清在线观看 | 变态另类成人亚洲欧美熟女| av黄色大香蕉| 黄色日韩在线| 日韩精品有码人妻一区| 国产成人影院久久av| 久久久久免费精品人妻一区二区| 精华霜和精华液先用哪个| 国产人妻一区二区三区在| 国产黄片美女视频| 舔av片在线| 国产高清三级在线| 国产一区二区在线观看日韩| 乱系列少妇在线播放| 国产精品电影一区二区三区| 日本 av在线| 亚洲人与动物交配视频| 国产精品亚洲一级av第二区| 精品免费久久久久久久清纯| 亚洲人成网站高清观看| 日本色播在线视频| bbb黄色大片| 亚洲成人免费电影在线观看| 超碰av人人做人人爽久久| 亚洲在线观看片| 3wmmmm亚洲av在线观看| 免费在线观看成人毛片| av在线蜜桃| 在线观看午夜福利视频| 天天一区二区日本电影三级| 国产一区二区激情短视频| 天堂影院成人在线观看| 少妇人妻一区二区三区视频| 少妇人妻精品综合一区二区 | 99在线人妻在线中文字幕| 别揉我奶头 嗯啊视频| 国产三级在线视频| 一本一本综合久久| 又紧又爽又黄一区二区| 久久6这里有精品| 在线播放无遮挡| 国产精品精品国产色婷婷| 搞女人的毛片| videossex国产| 熟女电影av网| 国产精品美女特级片免费视频播放器| 村上凉子中文字幕在线| 人妻久久中文字幕网| 日本 av在线| 国产aⅴ精品一区二区三区波| 97超视频在线观看视频| 精品久久久久久久人妻蜜臀av| 最近视频中文字幕2019在线8| 内射极品少妇av片p| 免费搜索国产男女视频| 亚洲第一区二区三区不卡| 国产精品一区www在线观看 | 久久久久精品国产欧美久久久| 亚洲一级一片aⅴ在线观看| 99热精品在线国产| 欧美一区二区亚洲| 国产黄片美女视频| 日韩亚洲欧美综合| 神马国产精品三级电影在线观看| 欧美成人免费av一区二区三区| 99久久精品热视频| 精品国内亚洲2022精品成人| 俄罗斯特黄特色一大片| 亚洲经典国产精华液单| 好男人在线观看高清免费视频| 97碰自拍视频| 一卡2卡三卡四卡精品乱码亚洲| 日韩欧美一区二区三区在线观看| 黄色欧美视频在线观看| 中文字幕熟女人妻在线| 少妇熟女aⅴ在线视频| www.www免费av| 成人av一区二区三区在线看| 干丝袜人妻中文字幕| 搡女人真爽免费视频火全软件 | 男女啪啪激烈高潮av片| 亚洲成a人片在线一区二区| 亚洲精品在线观看二区| 无人区码免费观看不卡| 变态另类成人亚洲欧美熟女| 午夜免费成人在线视频| 久久九九热精品免费| 免费看光身美女| 99在线人妻在线中文字幕| 国产精品伦人一区二区| 亚洲成人免费电影在线观看| 日韩强制内射视频| 精品人妻一区二区三区麻豆 | 性插视频无遮挡在线免费观看| 亚洲在线观看片| 国产av在哪里看| 欧美激情在线99| 精品一区二区三区av网在线观看| 日韩欧美国产一区二区入口| 一级毛片久久久久久久久女| 成年免费大片在线观看| 国产精品98久久久久久宅男小说| 天堂动漫精品| 久久精品国产99精品国产亚洲性色| 国产精品1区2区在线观看.| 国产精品人妻久久久影院| 国内精品久久久久久久电影| av天堂中文字幕网| 欧美日韩黄片免| av在线天堂中文字幕| 国产精品永久免费网站| 国产真实伦视频高清在线观看 | 国产极品精品免费视频能看的| 春色校园在线视频观看| 精品不卡国产一区二区三区| 性插视频无遮挡在线免费观看| 一个人看的www免费观看视频| 少妇猛男粗大的猛烈进出视频 | 成人特级av手机在线观看| 精品久久久久久久久久免费视频| 久久中文看片网| 亚洲人成网站高清观看| 亚州av有码| 亚洲在线观看片| 国产亚洲91精品色在线| 免费人成在线观看视频色| 少妇的逼好多水| 天堂av国产一区二区熟女人妻| 床上黄色一级片| 亚洲av中文字字幕乱码综合| 欧美日韩中文字幕国产精品一区二区三区| 午夜视频国产福利| 狂野欧美激情性xxxx在线观看| 赤兔流量卡办理| 中文字幕熟女人妻在线| 色5月婷婷丁香| 日本爱情动作片www.在线观看 | 看十八女毛片水多多多| 欧美xxxx黑人xx丫x性爽| 三级男女做爰猛烈吃奶摸视频| www.www免费av| 精品一区二区三区视频在线观看免费| 午夜精品久久久久久毛片777| 国产爱豆传媒在线观看| 中出人妻视频一区二区| 亚洲av成人av| 又黄又爽又免费观看的视频| 午夜爱爱视频在线播放| 国产精品福利在线免费观看| 成年女人毛片免费观看观看9| eeuss影院久久| 精品一区二区免费观看| 亚洲人成网站在线播| 国产不卡一卡二| 欧美高清性xxxxhd video| 亚洲美女视频黄频| 别揉我奶头 嗯啊视频| 久久精品国产亚洲av香蕉五月| 久久精品国产鲁丝片午夜精品 | 亚洲欧美日韩高清专用| 欧美一区二区精品小视频在线| 可以在线观看的亚洲视频| 嫩草影视91久久| 两人在一起打扑克的视频| 亚洲美女黄片视频| 亚洲精品456在线播放app | 蜜桃久久精品国产亚洲av| 一进一出好大好爽视频| 国产aⅴ精品一区二区三区波| 天堂动漫精品| 中文字幕av在线有码专区| 最近视频中文字幕2019在线8| 中国美白少妇内射xxxbb| 国产高清有码在线观看视频| 一个人看的www免费观看视频| 午夜福利视频1000在线观看| 老熟妇仑乱视频hdxx| 99热这里只有是精品在线观看| 在线观看一区二区三区| 一卡2卡三卡四卡精品乱码亚洲| 日韩欧美三级三区| 伊人久久精品亚洲午夜| 男人和女人高潮做爰伦理| 97碰自拍视频| 女的被弄到高潮叫床怎么办 | 最近最新中文字幕大全电影3| 床上黄色一级片| 亚洲三级黄色毛片| 91麻豆精品激情在线观看国产| 在线免费十八禁| 99热精品在线国产| 日韩中字成人| 国产日本99.免费观看| 国产一区二区三区视频了| 丝袜美腿在线中文| 日日摸夜夜添夜夜添av毛片 | 男女视频在线观看网站免费| 亚洲国产日韩欧美精品在线观看| 波多野结衣高清作品| 国产久久久一区二区三区| 久久欧美精品欧美久久欧美| 久久九九热精品免费| av天堂中文字幕网| 久久久久久久精品吃奶| 国产伦一二天堂av在线观看| 久久久久久国产a免费观看| 乱码一卡2卡4卡精品| avwww免费| 女人被狂操c到高潮| 欧美又色又爽又黄视频| 99在线视频只有这里精品首页| 国产老妇女一区| 国产黄色小视频在线观看| 少妇的逼水好多| 有码 亚洲区| 欧美最黄视频在线播放免费| 亚洲av免费在线观看| 中文字幕人妻熟人妻熟丝袜美| 久久精品国产自在天天线| 国产单亲对白刺激| 亚洲精品一区av在线观看| 国产视频内射| 国产精品,欧美在线| 老司机深夜福利视频在线观看| 亚洲成人久久爱视频| www.色视频.com| 免费看光身美女| x7x7x7水蜜桃| 欧美精品国产亚洲| 亚洲专区国产一区二区| 久久午夜亚洲精品久久|