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

    黃土高原旱地麥田26年免耕覆蓋對(duì)土壤肥力及原核微生物群落多樣性的影響*

    2019-03-08 06:12:22張貴云呂貝貝張麗萍范巧蘭魏明峰袁嘉瑋柴躍進(jìn)
    關(guān)鍵詞:原核保護(hù)性耕作

    張貴云, 呂貝貝, 張麗萍, 劉 珍, 范巧蘭, 魏明峰, 姚 眾, 袁嘉瑋, 柴躍進(jìn)

    ?

    黃土高原旱地麥田26年免耕覆蓋對(duì)土壤肥力及原核微生物群落多樣性的影響*

    張貴云1, 呂貝貝1, 張麗萍1, 劉 珍1, 范巧蘭1, 魏明峰1, 姚 眾1, 袁嘉瑋1, 柴躍進(jìn)2

    (1. 山西省農(nóng)業(yè)科學(xué)院棉花研究所 運(yùn)城 044000; 2. 山西省臨汾市農(nóng)機(jī)局 臨汾 041000)

    為明確黃土高原旱作麥田長(zhǎng)期保護(hù)性耕作對(duì)土壤肥力和土壤原核微生物的效應(yīng), 以位于山西省臨汾市實(shí)施保護(hù)性耕作26年的小麥田為試驗(yàn)基地, 采用Illumina Hiseq 2500高通量測(cè)序等手段, 開展了不同耕作措施[免耕覆蓋(NTS)、深松免耕覆蓋(SNTS)和傳統(tǒng)耕作(TT1)]對(duì)土壤理化性質(zhì)和土壤原核微生物多樣性的影響分析。研究結(jié)果表明: 1)NTS和SNTS處理比TT1處理顯著提高了土壤全氮、堿解氮、速效磷及速效鉀的含量, 降低了土壤pH, 提高了土壤貯水能力和水分含量, 降低了0~10 cm土層的土壤容重, 但提高了10~20 cm土層的土壤容重; 同時(shí), SNTS處理顯著增加了土壤的有機(jī)質(zhì)含量。2)Illumina Hiseq高通量16S rRNA基因V4區(qū)測(cè)序結(jié)果表明: NTS和SNTS處理比TT1處理顯著降低了綠彎菌門的相對(duì)豐度; NTS比SNTS處理顯著降低了土壤中疣微菌門和綠彎菌門的相對(duì)豐度; NTS處理顯著增加了土壤原核微生物群落的多樣性, 但未顯著改變?cè)宋⑸锶郝涞呢S度; SNTS對(duì)原核微生物群落的多樣性和豐富度均未有顯著改變; NTS處理的顯著性差異物種(Biomarker)高于其他2個(gè)處理; 其他原核微生物門的相對(duì)豐度, 在3個(gè)處理間尚未有明顯差異。3)聚類分析可見: NTS和SNTS處理與TT1處理的微生物群落結(jié)構(gòu)差異較大; NTS處理與SNTS處理間的微生物群落結(jié)構(gòu)差異較小。4)CCA分析可知: 土壤pH、有機(jī)質(zhì)、速效氮、速效磷、速效鉀含量對(duì)土壤原核微生物群落遺傳多樣性的變化起著重要作用; 與TT1處理相比, NTS和SNTS處理在一定程度上改變了土壤原核微生物群落結(jié)構(gòu), 但仍存在結(jié)構(gòu)的相似性。綜上所述, 長(zhǎng)期進(jìn)行NTS和SNTS處理對(duì)黃土高原旱地麥田土壤微生物多樣性、豐富度以及土壤肥力因子等的正效作用明顯。

    保護(hù)性耕作; 免耕覆蓋; 深松免耕; 高通量測(cè)序; 土壤原核微生物

    保護(hù)性耕作是以秸稈覆蓋還田、農(nóng)田免耕或少耕、施肥、播種等復(fù)式作業(yè)為主要內(nèi)容的一種新型旱作機(jī)械化耕作方式[1-2]。與傳統(tǒng)耕作方式相比, 保護(hù)性耕作不僅省工、省力、減少投入, 而且可減少水土流失, 保持土層結(jié)構(gòu), 改善土壤理化性質(zhì), 促進(jìn)微生物的繁殖[3-5]。研究表明, 不同類型生境中, 微生物的群落組成存在明顯差異[6], 耕作方式不同, 土壤環(huán)境也會(huì)不同, 進(jìn)而引起土壤微生物的數(shù)量和結(jié)構(gòu)組成發(fā)生改變[7]。Yan等[8]對(duì)耕作和免耕變性土微生物群落功能多樣性的研究顯示, 免耕比傳統(tǒng)耕作更有利于增加土壤微生物多樣性。李彤等[9]研究發(fā)現(xiàn)西北地區(qū)旱作麥田采用保護(hù)性耕作, 可以影響土壤微生物群落豐度, 并顯著影響土壤理化性質(zhì), 進(jìn)而影響土壤微生物空間結(jié)構(gòu)。

    20世紀(jì)60年代保護(hù)性耕作傳入我國(guó)[10], 由于其能夠增加土壤蓄水, 增強(qiáng)土壤抗旱能力, 提高作物產(chǎn)量, 降低生產(chǎn)成本等, 在我國(guó)北方旱地得到了廣泛應(yīng)用[11]。為提高我國(guó)北方干旱地區(qū)農(nóng)業(yè)生產(chǎn)水平, 1992年中國(guó)農(nóng)業(yè)大學(xué)和山西省農(nóng)機(jī)局受澳大利亞國(guó)際農(nóng)業(yè)研究中心資助, 在山西省臨汾市開展了旱地保護(hù)性耕作體系研究。張貴云等[12]以該試驗(yàn)基地為平臺(tái), 開展了長(zhǎng)期保護(hù)性耕作對(duì)叢枝菌根真菌多樣性的影響研究, 結(jié)果表明, 保護(hù)性耕作提高了叢枝菌根真菌(AMF)多樣性, 且深松與免耕相結(jié)合的保護(hù)性耕作方式, 更有利于提高土壤AMF多樣性。土壤細(xì)菌是土壤生態(tài)系統(tǒng)中不可或缺的部分, 參與多種土壤生理生化過程, 如土壤有機(jī)質(zhì)的分解、腐殖質(zhì)的形成、養(yǎng)分的轉(zhuǎn)化和循環(huán)等, 土壤生態(tài)系統(tǒng)的質(zhì)量和健康與土壤中細(xì)菌群落結(jié)構(gòu)的變化有很大關(guān)系[13-14], 因而對(duì)土壤微生物的研究尤其是原核微生物群落結(jié)構(gòu)的研究尤為重要。近年來, 人們運(yùn)用高通量手段研究土壤微生物的技術(shù)日趨成熟, 這為深入分析微生物資源提供了可能[15]。

    高通量測(cè)序技術(shù)是對(duì)傳統(tǒng)測(cè)序技術(shù)的一次革命性創(chuàng)新[16], 其產(chǎn)生的16S rRNA基因的測(cè)序數(shù)據(jù), 覆蓋深度高, 可以用來估計(jì)微生物群落的物種組成, 能更真實(shí)地揭示環(huán)境中微生物群落的復(fù)雜度和多樣性[17]。目前為止, 國(guó)內(nèi)外有關(guān)保護(hù)性耕作田土壤微生物多樣性的研究有了一定的報(bào)道, 但是黃土高原旱作麥田長(zhǎng)期保護(hù)性耕作對(duì)土壤原核微生物群落的影響, 尚鮮見研究報(bào)道。本研究以山西臨汾旱地實(shí)施保護(hù)性耕作26年的小麥田為試驗(yàn)平臺(tái), 采用Illumina Hiseq 2500高通量測(cè)序平臺(tái)測(cè)序土壤16S rDNA, 結(jié)合土壤質(zhì)量環(huán)境因子的典型相關(guān)分析(CCA), 對(duì)不同保護(hù)性耕作措施對(duì)土壤原核微生物群落和土壤環(huán)境的影響進(jìn)行研究, 探尋保護(hù)性耕作技術(shù)應(yīng)用的新理論和新策略, 為在微生物水平上研究黃土高原旱作地區(qū)保護(hù)性耕作技術(shù)提供重要理論支撐。

    1 材料與方法

    1.1 研究區(qū)域概況

    試驗(yàn)基地位于山西省臨汾市堯都區(qū)城隍村(36°01¢53″N, 111°37¢46″E, 海拔555 m)。試驗(yàn)區(qū)域年平均氣溫12.2 ℃, 年0 ℃以上積溫4 600 ℃, 無霜期195 d, 年平均地溫13.8 ℃。年均降水量515.7 mm, 且均集中在7—9月份; 年均蒸發(fā)量1 933.1 mm。土壤類型為碳酸鹽褐土, 試驗(yàn)地種植制度為長(zhǎng)期單一種植冬小麥()。

    1.2 耕作處理方式

    2017年5月, 借助1992年在山西省臨汾市開始的保護(hù)性耕作體系長(zhǎng)期定位試驗(yàn)[18], 開展了黃土高原旱作麥田長(zhǎng)期保護(hù)性耕作對(duì)土壤原核微生物群落多樣性及土壤肥力的影響研究。本試驗(yàn)選擇26年的免耕覆蓋(NTS)、連續(xù)深松覆蓋4年后免耕22年的免耕覆蓋(SNTS)以及26年的傳統(tǒng)耕作(TT1)3種耕作處理方式。其中, NTS為全年不耕作, 播種時(shí)用免耕播種機(jī)一次性完成施肥和播種, 收獲后將各小區(qū)秸稈粉碎全覆蓋; SNTS是深松(深度約40 cm)后, 進(jìn)行長(zhǎng)期的免耕覆蓋處理; TT1是小麥?zhǔn)斋@后用鏵式犁耕翻(深度約25 cm), 無秸稈覆蓋, 播前旋耕(深度約15 cm)整地。各處理重復(fù)3次, 每小區(qū)面積為600 m2。試驗(yàn)地?zé)o灌溉條件, 作物靠自然降水生長(zhǎng)。播種時(shí)間: 2016年9月25日, 小麥品種為‘長(zhǎng)6359’, 播種量為120 kg·hm-2。試驗(yàn)田化肥使用量:尿素為225 kg·hm-2, 磷酸二銨為150 kg·hm-2。

    1.3 土壤樣品采集

    于2017年5月8日(小麥灌漿期)采取小麥根際土樣, 每小區(qū)隨機(jī)選3點(diǎn), 每個(gè)樣點(diǎn)先去掉表土, 再?gòu)?~20 cm土層取小麥根際土樣, 每點(diǎn)采集3份土樣均勻混合成一個(gè)樣品。將所采集的土樣分為兩部分, 一部分保存于滅菌的塑封袋中, 于低溫條件下帶回實(shí)驗(yàn)室, 保存于-80 ℃冰箱中用于提取DNA, 供土壤原核微生物群落結(jié)構(gòu)多樣性分析; 另一部分保存于普通塑封袋中, 帶回實(shí)驗(yàn)室自然風(fēng)干, 根據(jù)試驗(yàn)需要分別過18目、100目篩子用于土壤理化性質(zhì)的測(cè)定。土壤水分和容重指標(biāo)則分別取自同一樣點(diǎn)0~10 cm、10~20 cm土層。

    1.4 測(cè)定方法

    土壤理化性質(zhì)測(cè)定: pH采用1∶2.5土水比浸提酸度計(jì)法, 土壤水分采用烘干法, 土壤容重采用環(huán)割法, 有機(jī)質(zhì)采用重鉻酸鉀容量法—稀釋熱法, 全氮采用半微量開氏法, 堿解氮采用堿解擴(kuò)散法, 速效磷采用0.5 mol×L-1NaHCO3法, 速效鉀采用NH4OAc浸提—火焰光度法測(cè)定[19]。

    土壤總DNA的提取及16S rRNA基因的PCR擴(kuò)增: 土壤總DNA基因組的提取采用DNA提取試劑盒進(jìn)行。采用引物515F(5¢-GTGCCAGCMGCCGC-GGTAA-3¢)和806R(5¢-GGACTACHVGGGTWTCT AAT-3¢)擴(kuò)增土壤微生物16S rRNA基因V4區(qū)。30 μL PCR反應(yīng)體系: Phusion Master Mix (New England Biolabs, 2×)15 μL, Primer (2 μmol×L-1) 3 μL, Genomic DNA (1 ng×μL-1) 10 μL, dd H2O 2 μL。PCR反應(yīng)條件: 98 ℃預(yù)變性1 min; 98 ℃ 10 s, 50 ℃ 30 s, 72 ℃ 30 s, 30個(gè)循環(huán); 72 ℃, 5 min。每個(gè)樣本3次重復(fù), 根據(jù)PCR產(chǎn)物濃度進(jìn)行等濃度混樣, 充分混勻后使用1×TAE濃度2%的瓊脂糖膠電泳純化PCR產(chǎn)物, 選擇主帶大小在400~450 bp之間的序列, 割膠回收目標(biāo)條帶。產(chǎn)物純化試劑盒使用的是Thermo Scientific公司GeneJET膠回收試劑盒。

    原核微生物16S rRNA基因V4區(qū)高通量測(cè)序: 經(jīng)純化檢測(cè)之后的土樣送至北京諾禾致源科技股份有限公司, 進(jìn)行基因組DNA的提取和PCR擴(kuò)增、PCR產(chǎn)物的混樣和純化, 使用New England Biolabs公司的NEB Next?Ultra?DNA Library Prep Kit for Illumina建庫(kù)試劑盒進(jìn)行文庫(kù)的構(gòu)建, 構(gòu)建好的文庫(kù)經(jīng)過Qubit定量和文庫(kù)檢測(cè)合格后, 使用Illumina公司的HiSeq 2500測(cè)序平臺(tái), 進(jìn)行上機(jī)測(cè)序。

    運(yùn)用Uparse軟件對(duì)序列進(jìn)行聚類, 以97%的一致性(Identity)將序列聚類成為OTUs。采用軟件RDP classifier進(jìn)行物種分類。分別在門(phylum)和屬(genus)水平上統(tǒng)計(jì)各樣本的群落組成, 并對(duì)OTU進(jìn)行生態(tài)學(xué)分析, 使用Qiime軟件計(jì)算土壤物種相對(duì)豐度、Chao1、ACE、Shannon、Simpson指數(shù)等, 使用R軟件繪制稀釋曲線。對(duì)OTUs進(jìn)行多序列比對(duì), 用Qiime軟件構(gòu)建系統(tǒng)發(fā)生樹, 獲得不同樣品和分組的群落結(jié)構(gòu)差異, 通過UPGMA樣品聚類樹進(jìn)行展示。選用LEfSe統(tǒng)計(jì)分析方法對(duì)分組樣品的物種組成和群落結(jié)果進(jìn)行差異顯著性檢驗(yàn)。同時(shí), 進(jìn)行CCA分析和多樣性指數(shù)與環(huán)境因子間的相關(guān)性分析, 得到顯著影響組間群落變化的環(huán)境因子。

    采用SPSS 20.0對(duì)土壤理化數(shù)據(jù)進(jìn)行統(tǒng)計(jì)分析, 使用R語(yǔ)言分析微生物多樣性數(shù)據(jù), 運(yùn)用Canoco軟件進(jìn)行環(huán)境因子對(duì)微生物群落結(jié)構(gòu)的影響分析。

    2 結(jié)果與分析

    2.1 保護(hù)性耕作措施對(duì)土壤理化性質(zhì)的影響

    由表1可以看出, 在0~20 cm土層土壤中, NTS和SNTS 2個(gè)處理的土壤全氮、堿解氮、速效磷及速效鉀含量均顯著高于TT1處理, 而其pH顯著低于TT1, SNTS處理的土壤有機(jī)質(zhì)含量顯著高于NTS和TT1處理。TT1處理的土壤含水量和容重與其他處理間差異均顯著, NTS和SNTS處理間差異不顯著, 3個(gè)處理土壤含水量依次為NTS>SNTS>TT1; 0~10 cm土層土壤容重依次為TT1>NTS>SNTS, 10~20 cm土層土壤容重為NTS>SNTS>TT1。結(jié)果顯示, 保護(hù)性耕作與傳統(tǒng)耕作相比, 顯著提高了土壤中的全氮、堿解氮、速效磷及速效鉀的含量, 顯著降低了土壤pH, 且深松免耕處理顯著增加了土壤有機(jī)質(zhì)的含量; 同時(shí), 保護(hù)性耕作還可以增加土壤貯水能力, 提高土壤水分含量, 降低0~10 cm土層的土壤容重, 但提高10~20 cm土層的土壤容重。

    表1 長(zhǎng)期不同耕作方式下土壤的理化性質(zhì)

    同列數(shù)字后不同小寫字母表示差異顯著(<0.05)。NTS: 免耕覆蓋; SNTS: 深松免耕覆蓋; TT1: 傳統(tǒng)耕作。Different lowercase letters in the same column indicate significant differences at 0.05 level. NTS: no-tillage with stubble; SNTS: no-tillage with stubble after subsoiling; TT1: traditional tillage.

    2.2 土壤樣品測(cè)序深度及合理性分析

    為檢測(cè)樣本的數(shù)據(jù)量是否合理, 從樣品中隨機(jī)抽取一定測(cè)序量的數(shù)據(jù), 統(tǒng)計(jì)它們所代表物種數(shù)目(即OTUs數(shù)目), 以抽取的測(cè)序數(shù)據(jù)量與對(duì)應(yīng)的物種數(shù)來構(gòu)建稀釋曲線, 不同耕作方式土壤樣品原核微生物稀釋曲線見圖1。由圖1可知, 隨著測(cè)序量的增加, 原核微生物稀釋曲線逐漸趨于平坦, 但尚未達(dá)到飽和, 更多的測(cè)序量只會(huì)產(chǎn)生少量新的OTU, 表明本次試驗(yàn)的測(cè)序數(shù)據(jù)量合理, 能夠較真實(shí)地反映這3組土壤樣品的原核微生物群落, 說明測(cè)序結(jié)果包含了大多數(shù)原核微生物類群, 基本可以反映不同耕作下土壤生境中原核微生物群落結(jié)構(gòu)與組成, 但仍有少量的原核微生物種類未被發(fā)現(xiàn)。

    2.3 保護(hù)性耕作措施對(duì)土壤原核微生物物種豐度的影響

    不同處理土壤原核微生物種類OTU數(shù)量關(guān)系韋恩圖如圖2所示, 3個(gè)處理共產(chǎn)生9 235個(gè)OTU, 共同包含的OTU個(gè)數(shù)為4 466, 占總OTU數(shù)量的48.36%。其中, NTS共得到7 257個(gè)OTU, SNTS共得到6 547個(gè)OTU, TT1共得到6 384個(gè)OTU。NTS中特有的OTU為1 314個(gè), 占總OTU數(shù)量的14.23%; SNTS中特有的OTU為732個(gè), 占總OTU數(shù)量的7.93%; TT1中特有的OTU為702個(gè), 占總OTU數(shù)量的7.60%。表明3個(gè)處理原核微生物的群落結(jié)構(gòu)差異較大。分別將3組兩兩比較發(fā)現(xiàn), TT1和NTS共有的OTU為5 138個(gè), TT1和SNTS共有的OTU最少, 為5 010個(gè), NTS和SNTS共有的OTU最多, 為5 271個(gè)。

    圖1 長(zhǎng)期不同耕作方式下土壤原核微生物稀釋曲線

    NTS: 免耕覆蓋; SNTS: 深松免耕覆蓋; TT1: 傳統(tǒng)耕作。NTS: no-tillage with stubble; SNTS: no-tillage with stubble after subsoiling; TT1: traditional tillage.

    圖2 長(zhǎng)期不同耕作方式下土壤原核微生物Venn圖

    NTS: 免耕覆蓋; SNTS: 深松免耕覆蓋; TT1: 傳統(tǒng)耕作。NTS: no-tillage with stubble; SNTS: no-tillage with stubble after subsoiling; TT1: traditional tillage.

    2.4 保護(hù)性耕作措施對(duì)土壤原核微生物群落組成及相對(duì)豐度的影響

    門分類水平上, 不同耕作方式下原核微生物群落組成及相對(duì)豐度見圖3。從圖3可知, 3個(gè)處理土壤原核微生物群落在門水平上的組成基本相同, 但平均相對(duì)豐度所占比例略有差異。各處理中原核微生物類群相對(duì)豐度排名前10的有變形菌門(Proteobacteria, 占31.45%)、酸桿菌門(Acidobacteria, 占16.51%)、放線菌門(Actinobacteria, 占14.94%)、厚壁菌門(Firmicutes, 占3.69%)、擬桿菌門(Bacteroidetes, 占8.74%)、芽單胞菌門(Gemmatimonadetes, 占6.51%)、浮霉菌門(Planctomycetes, 占3.09%)、奇古菌門(Thaumarchaeota, 占2.87%)、疣微菌門(Verrucomicrobia, 占3.08%)、綠彎菌門(Chloroflexi, 占2.65%), 它們占總類群的92.43%~93.84%。其中, 變形菌門、酸桿菌門和放線菌門為明顯優(yōu)勢(shì)菌群, 分別占3個(gè)處理平均相對(duì)豐度10%以上, NTS處理3次重復(fù)平均相對(duì)豐度分別為34.24%、13.23%、14.06%, SNTS處理分別為28.02%、18.48%、17.02%, TT1處理分別為32.07%、17.83%、13.73%。統(tǒng)計(jì)分析結(jié)果顯示, 3個(gè)處理中優(yōu)勢(shì)菌門平均相對(duì)豐度均未達(dá)到顯著性差異水平(0.05)。其他7個(gè)菌門的平均相對(duì)豐度較低, 所占比例均小于10.0%。統(tǒng)計(jì)分析結(jié)果顯示, 除疣微菌門在NTS和SNTS處理中平均相對(duì)豐度達(dá)到顯著差異水平(0.05), 綠彎菌門在NTS和SNTS、NTS和TT1處理中平均相對(duì)豐度達(dá)到顯著差異水平(0.05)外, 厚壁菌門、擬桿菌門、芽單胞菌門、浮霉菌門、奇古菌門在3個(gè)處理中均未有明顯差異。此外, 還有6.47%的原核微生物類群未被分類。結(jié)果表明, 與傳統(tǒng)耕作相比, 免耕顯著降低了綠彎菌門的相對(duì)豐度; 與深松免耕相比, 免耕顯著降低了土壤中疣微菌門和綠彎菌門的相對(duì)豐度, 其他原核微生物門的相對(duì)豐度在3個(gè)處理間尚未有明顯差異。

    圖3 長(zhǎng)期不同耕作方式下土壤原核微生物群落在門水平上的組成及相對(duì)豐度

    NTS: 免耕覆蓋; SNTS: 深松免耕覆蓋; TT1: 傳統(tǒng)耕作。NTS: no-tillage with stubble; SNTS: no-tillage with stubble after subsoiling; TT1: traditional tillage.

    屬水平上, 不同耕作方式下原核微生物群落的組成及相對(duì)豐度如圖4所示。在屬水平上, 3個(gè)處理中原核微生物類群的平均相對(duì)豐度均小于5.0%, 其中RB41屬是3個(gè)處理共同的優(yōu)勢(shì)菌群(占3.26%), 在TT1處理中的相對(duì)豐度最高, 未被分類的原核微生物類群占到86.02%。在3種處理的土壤樣品中, 未分類的原核微生物成了第一大類群, 說明隨著分類的細(xì)化, 相關(guān)數(shù)據(jù)庫(kù)中可用信息量越來越少, 同時(shí)也說明土壤中含有一定數(shù)量的潛在原核微生物類群。

    圖4 長(zhǎng)期不同耕作方式下土壤原核微生物群落在屬水平上的組成及相對(duì)豐度

    TT1: 傳統(tǒng)耕作; NTS: 免耕覆蓋; SNTS: 深松免耕覆蓋。TT1: traditional tillage; NTS: no-tillage with stubble; SNTS: no-tillage with stubble after subsoiling; TT1: traditional tillage.

    使用Qiime軟件計(jì)算Unifrac距離, 并基于Weighted Unifrac方法構(gòu)建UPGMA樣品聚類樹, 對(duì)3個(gè)處理的土壤原核微生物群落構(gòu)成的相似性進(jìn)行聚類分析, 并將聚類結(jié)果與各樣品在門水平上的物種相對(duì)豐度進(jìn)行整合展示(圖5)。結(jié)果表明, 3個(gè)處理可分為2大類, 其中NTS和SNTS處理的土壤原核微生物群落組成趨于一類, 在聚類關(guān)系上趨同性較強(qiáng), TT1處理單獨(dú)為一類。表明免耕土壤原核微生物種群結(jié)構(gòu)與深松免耕覆蓋并沒有明顯差異, 親緣關(guān)系較近, 說明深松4年之后免耕覆蓋22年土壤原核微生物種群結(jié)構(gòu)變化趨于免耕。

    圖5 長(zhǎng)期不同耕作方式下土壤原核微生物群落相似性聚類樹

    NTS: 免耕覆蓋; SNTS: 深松免耕覆蓋; TT1: 傳統(tǒng)耕作。NTS: no-tillage with stubble; SNTS: no-tillage with stubble after subsoiling; TT1: traditional tillage.

    2.5 保護(hù)性耕作措施對(duì)土壤原核微生物多樣性的影響

    Chao1指數(shù)和ACE指數(shù)表征物種的豐富度, 其值越高表明原核微生物群落的物種豐富度越高; Shannon指數(shù)和Simpson指數(shù)表征微生物多樣性程度, Shannon指數(shù)越高表明原核微生物群落的多樣性越高, Simpson指數(shù)與之相反; Coverage是指各樣本文庫(kù)的覆蓋率[20-22]。

    表2是不同耕作方式下土壤原核微生物的多樣性指數(shù)。從表2可以看出, 3個(gè)處理文庫(kù)覆蓋率均在98%以上, 能夠較為全面地反映不同耕作措施下土壤原核微生物群落的種類和結(jié)構(gòu)。綜合豐富度指數(shù)Chao1和ACE以及多樣性指數(shù)Shannon和Simpson指數(shù), 3個(gè)處理土壤環(huán)境原核微生物群落多樣性由高到低排序?yàn)? NTS>SNTS>TT1。統(tǒng)計(jì)分析結(jié)果, 3個(gè)處理的Chao1指數(shù)和ACE指數(shù)差異不顯著。不同耕作處理對(duì)原核微生物群落的多樣性影響程度均不同, 與傳統(tǒng)耕作相比, 免耕顯著增加了土壤原核微生物群落的多樣性, 未顯著改變?cè)宋⑸锶郝涞呢S富度; 深松免耕處理對(duì)原核微生物群落的多樣性和豐富度均未有顯著改變。

    表2 長(zhǎng)期不同耕作方式下土壤原核微生物群落的多樣性指數(shù)

    同列數(shù)字后不同小寫字母表示差異顯著(<0.05)。NTS: 免耕覆蓋; SNTS: 深松免耕覆蓋; TT1: 傳統(tǒng)耕作。Different lowercase letters in the same column indicate significant differences at 0.05 level. NTS: no-tillage with stubble; SNTS: no-tillage with stubble after subsoiling; TT1: traditional tillage.

    2.6 不同耕作措施下土壤原核微生物群落組間差異性

    通過LEfSe分析, 可以找出組間具有顯著差異的Biomarker, 以及顯著影響組間差異性的物種或群落[23]。不同耕作方式下土壤原核微生物群落組間差異LEfSe分析柱狀圖和進(jìn)化分支圖如圖6所示, 不同處理對(duì)應(yīng)不同的Biomarker, NTS處理的顯著性差異物種(Biomarker)最多, SNTS次之, TT1最少。在門水平上, TT1樣品中有1個(gè)主要門(疣微菌門)的細(xì)菌存在顯著差異, SNTS樣品中有1個(gè)主要門(奇古菌門)的古菌存在顯著差異, TT1樣品在門水平上無顯著差異。在屬水平上, SNTS中放線菌門的、變形菌門的共2屬的細(xì)菌存在顯著差異, NTS中變形菌門的、衣原體門的、擬桿菌門的共3屬的細(xì)菌存在顯著差異。

    圖6 長(zhǎng)期不同耕作方式下土壤原核微生物群落處理間差異LEfSe分析圖(A)和進(jìn)化分支圖(B)

    NTS: 免耕覆蓋; SNTS: 深松免耕覆蓋; TT1: 傳統(tǒng)耕作。NTS: no-tillage with stubble; SNTS: no-tillage with stubble after subsoiling; TT1: traditional tillage.

    圖7 門水平下長(zhǎng)期不同耕作處理的土壤原核微生物群落與土壤理化因子的典型對(duì)應(yīng)分析

    TN、AN、AP、AK、OM分別代表土壤全氮、堿解氮、速效磷、速效鉀、有機(jī)質(zhì)含量。NTS: 免耕覆蓋; SNTS: 深松免耕覆蓋; TT1: 傳統(tǒng)耕作。TN, AN, AP, AK and OM represented total N, available N, available P, available K and organic matter, respectively. NTS: no-tillage with stubble; SNTS: no-tillage with stubble after subsoiling; TT1: traditional tillage.

    2.7 土壤原核微生物群落組成與土壤環(huán)境因子的相關(guān)分析

    土壤理化性質(zhì)的變化會(huì)直接或間接的影響土壤微生物群落結(jié)構(gòu)的組成[24-25]。為分析不同耕作方式土壤生境中原核微生物群落與土壤環(huán)境因子特性間的關(guān)系, 以原核微生物群落在門水平上的相對(duì)豐度數(shù)據(jù)為物種數(shù)據(jù), 土壤理化數(shù)據(jù)為土壤環(huán)境變量, 采用Canoco軟件進(jìn)行了CCA, 其結(jié)果排序如圖7所示, 第1排序軸對(duì)土壤原核微生物群落變化解釋最多, 達(dá)62.35%; 第2排序軸僅解釋了原核微生物群落變化的23.75%, 前兩軸共解釋了原核微生物群落總變異的86.10%, 在對(duì)原核微生物群落變異的解釋中起主導(dǎo)作用。其中, pH(=0.968 5)、速效氮(=-0.836 8)、速效磷(=-0.898 4)、速效鉀(=-0.994 1)與第1排序軸相關(guān)性較高, 全氮(=-0.719 6)、有機(jī)質(zhì)(=-0.917 7)與第2排序軸相關(guān)性較高。結(jié)果表明, 土壤pH、有機(jī)質(zhì)、速效氮、速效磷、速效鉀含量對(duì)土壤原核微生物群落遺傳多樣性的變化起著重要作用。3個(gè)處理樣本中, TT1.1、TT1.2、NTS.3相距較近, 說明這3個(gè)樣本的群落結(jié)構(gòu)相近, NTS.1、NTS.2、SNTS.2、TT1.3相距較近, 說明這4個(gè)樣本的群落結(jié)構(gòu)相近; NTS.1及NTS.2與NTS.3相距較遠(yuǎn), TT1.1及TT1.2與TT1.3相距較遠(yuǎn), SNTS.1、SNTS.2、SNTS.3相距較遠(yuǎn), 說明取樣過程中空間分布的不均勻及復(fù)雜性也會(huì)影響原核微生物的群落結(jié)構(gòu)。同時(shí), TT1的3個(gè)樣本與NTS.3、SNTS.1相距較遠(yuǎn), 說明這3組樣本之間群落結(jié)構(gòu)存在差異。結(jié)果表明, 與傳統(tǒng)耕作相比, 保護(hù)性耕作在一定程度上能夠改變土壤原核微生物群落結(jié)構(gòu), 但仍存在結(jié)構(gòu)的相似性。

    3 討論

    不同耕作方式對(duì)土壤養(yǎng)分的影響有所不同[26]。Dolan等[27]和Filho等[28]研究發(fā)現(xiàn), 長(zhǎng)期免耕提高了土壤表層有機(jī)質(zhì)的含量。本研究結(jié)果表明, 長(zhǎng)達(dá)26年的免耕和深松免耕處理, 改變了0~20 cm土層的理化性質(zhì), 其中土壤的有機(jī)質(zhì)、全氮、堿解氮、速效磷、速效鉀含量相對(duì)于傳統(tǒng)耕作均有不同程度的提高, 這一結(jié)論與楊培培等[29]、李友軍等[30]研究結(jié)果一致, 說明保護(hù)性耕作能夠提高土壤肥力, 促進(jìn)作物生長(zhǎng)。戴亮[31]在不同耕作方式對(duì)土壤微生物及土壤固碳能力影響的研究中指出, 深松免耕同時(shí)兼具翻耕和免耕的優(yōu)點(diǎn), 可以降低土壤有機(jī)質(zhì)的礦化速率, 加深土壤腐殖質(zhì)化程度, 減緩?fù)寥烙袡C(jī)碳的損失。許菁等[32]通過5年的定位試驗(yàn), 開展了長(zhǎng)期保護(hù)性耕作對(duì)小麥-玉米()兩熟農(nóng)田土壤碳儲(chǔ)量及固碳潛力的影響研究, 結(jié)果表明, 秸稈翻入土壤后在很大程度腐解, 與無秸稈還田土壤有機(jī)碳儲(chǔ)量之間的差距逐漸拉大, 土壤碳儲(chǔ)量達(dá)到峰值, 說明不同耕作方式使秸稈在土壤中的腐解環(huán)境不同, 引起秸稈腐解速率的不同, 土壤中秸稈腐解程度隨時(shí)間的推進(jìn)而提高。本試驗(yàn)結(jié)果發(fā)現(xiàn), 長(zhǎng)期保護(hù)性耕作與傳統(tǒng)耕作相比, 免耕處理和深松免耕處理的土壤有機(jī)質(zhì)含量也均有不同程度的提高, 且深松免耕優(yōu)于免耕, 前者顯著提高了土壤有機(jī)質(zhì)的含量, 較傳統(tǒng)耕作提高了61.50%, 可能因?yàn)樯钏擅飧龠M(jìn)了腐生菌增殖, 更有利于土壤有機(jī)質(zhì)的提高。

    保護(hù)性耕作可以增加土壤的貯水能力, 提高土壤水分含量, 這是因?yàn)楸Wo(hù)性耕作在土壤表面覆蓋有作物殘茬, 與傳統(tǒng)耕作相比其地表裸露面積減少, 從而土壤水分蒸發(fā)減少, 提高了土壤蓄水保墑能力[33]。本研究表明, 保護(hù)性耕作顯著提高了黃土高原旱作小麥田土壤含水量, 有較明顯的優(yōu)勢(shì)。土壤容重是衡量土壤肥力的重要指標(biāo)之一[34]。相關(guān)研究表明, 傳統(tǒng)翻耕的土壤容重大于保護(hù)性耕作[35], 本研究發(fā)現(xiàn)0~10 cm土層傳統(tǒng)耕作處理土壤容重高于免耕和深松免耕, 這與上述觀點(diǎn)相一致, 可能是由于保護(hù)性耕作措施增加了土壤的透氣性和有機(jī)質(zhì)含量, 在一定年限后有利于降低土壤容重[36]。而本研究3個(gè)處理10~20 cm和0~10 cm土層土壤容重變化規(guī)律不一致, 10~20 cm土層土壤容重變化表現(xiàn)為免耕處理高于傳統(tǒng)耕作和深松免耕, 與上述觀點(diǎn)不一致。趙洪利等[37]研究發(fā)現(xiàn), 免耕與深松、翻耕處理相比, 土壤容重增大, 和本研究10~20 cm土壤容重結(jié)論一致, 可能是因?yàn)楸驹囼?yàn)進(jìn)行的長(zhǎng)期免耕處理無土壤耕作, 在長(zhǎng)達(dá)26年的降水和自身重力影響下, 使得10~20 cm土層的土壤容重大于深松免耕和傳統(tǒng)耕作處理。CCA結(jié)果表明, 土壤原核微生物群落結(jié)構(gòu)多樣性是多個(gè)土壤環(huán)境因子共同作用的結(jié)果, 其中pH、有機(jī)質(zhì)、速效氮、速效磷、速效鉀含量對(duì)于研究區(qū)土壤原核微生物群落遺傳多樣性的變化起著重要作用, 說明土壤碳、氮等營(yíng)養(yǎng)元素對(duì)原核微生物的多樣性有顯著影響[24]。

    土壤微生物類群中, 以細(xì)菌居多, 通常能占到70%~90%, 其種類繁多、功能多樣, 具有最為豐富的遺傳多樣性[25,38]。不同耕作方式對(duì)土壤原核微生物群落結(jié)構(gòu)、種類、數(shù)量存在一定的影響, 采用高通量測(cè)序技術(shù)對(duì)土壤原核微生物16S rRNA基因的V4高變區(qū)進(jìn)行擴(kuò)增測(cè)序分析, 結(jié)果表明, 3種處理在97%的相似度水平上, 共產(chǎn)生9 878種屬分類水平的原核微生物(OTUs), 此外測(cè)序還發(fā)現(xiàn)了許多未被分類的原核微生物, 需通過深度測(cè)序或利用其他先進(jìn)手段對(duì)微生物進(jìn)行更細(xì)致地分類研究。在門分類水平上, 土壤主要優(yōu)勢(shì)菌群為變形菌門、酸桿菌門和放線菌門, 其中變形菌門的豐度最高。這與Roesch等[39]對(duì)土壤微生物的研究相一致, 可能是因?yàn)樽冃尉T是細(xì)菌中最大的一個(gè)門, 且該試驗(yàn)基地土壤呈弱堿性, 而變形菌門則是堿性土壤中的主要優(yōu)勢(shì)菌群[40]。有研究報(bào)道曾指出, 變形菌門中的許多類群可起到固氮作用, 并能適應(yīng)各種復(fù)雜環(huán)境, 環(huán)境條件的變化對(duì)其分布和相對(duì)豐度影響不大[41], 本試驗(yàn)中變形菌門在3個(gè)處理間差異并不顯著, 這一結(jié)論與上述觀點(diǎn)相符合。

    Doran[42]曾認(rèn)為, 免耕由于對(duì)土壤的機(jī)械擾動(dòng)較小, 土體基本處于一種相對(duì)厭氧的狀態(tài), 土壤微生物的數(shù)量及活性維持在相對(duì)穩(wěn)定水平, 進(jìn)而可能導(dǎo)致了土體芽孢細(xì)菌數(shù)量上升。Liu等[43]曾認(rèn)為, 土壤擾動(dòng)會(huì)降低土壤微生物的多樣性和活性。Bending等[44]曾認(rèn)為, 秸稈或根茬還田能提高土壤微生物的多樣性。本研究結(jié)果表明不同耕作措施對(duì)土壤原核微生物的多樣性影響均不同, 總體表現(xiàn)出保護(hù)性耕作高于傳統(tǒng)耕作, 且免耕覆蓋顯著提高了土壤原核微生物群落的豐富度和多樣性, 這一結(jié)論與上述觀點(diǎn)相符合??赡苁且?yàn)榻斩掃€田腐解過程中的部分產(chǎn)物, 如糖類、氨基酸、蛋白質(zhì)、維生素及有機(jī)酸、酚類等, 為微生物的生長(zhǎng)繁殖提供了豐富營(yíng)養(yǎng); 另一方面, 保護(hù)性耕作能改善土壤結(jié)構(gòu), 增加土壤透氣性, 土壤微生物種類隨之增加, 從而引起土壤原核微生物多樣性的增加[45-46]。表明長(zhǎng)期保護(hù)性耕作對(duì)提高土壤原核微生物群落豐富度和多樣性更有利, 免耕可能會(huì)更有效地為當(dāng)季作物提供養(yǎng)分[7]。

    4 結(jié)論

    黃土高原旱作區(qū)長(zhǎng)期免耕覆蓋、連續(xù)深松覆蓋免耕覆蓋以及傳統(tǒng)耕作下, 0~20 cm土層土壤中主要優(yōu)勢(shì)菌群為變形菌門、酸桿菌門和放線菌門, 均未出現(xiàn)獨(dú)有種群。與傳統(tǒng)耕作相比, 保護(hù)性耕作的原核微生物差異物種(Biomarker)明顯多于傳統(tǒng)耕作。不同耕作措施雖在一定程度上不會(huì)改變土壤原核微生物的組成和種類, 但會(huì)影響其相對(duì)分布。

    黃土高原旱作區(qū)長(zhǎng)期保護(hù)性耕作措施, 提高了土壤有機(jī)質(zhì)含量、全氮、堿解氮、速效磷、速效鉀含量, 降低了土壤pH, 改善了土壤性狀。其中, pH、有機(jī)質(zhì)、速效氮、速效磷、速效鉀含量是影響土壤原核微生物群落遺傳多樣性重要理化因子。

    本研究試驗(yàn)結(jié)果表明, 黃土高原旱作區(qū)連續(xù)多年的免耕處理是相對(duì)較優(yōu)的耕作處理措施, 適宜該地區(qū)積極推廣。

    本研究中, 3個(gè)處理土壤樣品的16S rDNA高通量測(cè)序獲取了大量的序列, 在一定程度上為后續(xù)研究提供了一定的理論支撐, 但由于稀釋曲線仍未完全飽和, 因此, 還有待進(jìn)行更深一步的研究。另外, 有關(guān)保護(hù)性耕作對(duì)土壤微生物的影響還需結(jié)合多學(xué)科、多方法進(jìn)行綜合評(píng)價(jià), 以便獲得更有效和全面的研究結(jié)論。在后續(xù)研究中, 我們將增設(shè)不同年份的取樣試材, 以便能夠更完整地反映出不同耕作方式對(duì)土壤環(huán)境的動(dòng)態(tài)影響。

    [1] 向新華, 魏巍, 張興義, 等. 保護(hù)性耕作對(duì)大豆生長(zhǎng)發(fā)育及土壤微生物多樣性影響[J]. 大豆科學(xué), 2013, 32(3): 321–327 XIANG X H, WEI W, ZHANG X Y, et al. Effect of conservation tillage systems on soybean growth and soil microbial diversity[J]. Soybean Science, 2013, 32(3): 321–327

    [2] 高旺盛. 論保護(hù)性耕作技術(shù)的基本原理與發(fā)展趨勢(shì)[J]. 中國(guó)農(nóng)業(yè)科學(xué), 2007, 40(12): 2702–2708 GAO W S. Development trends and basic principles of conservation tillage[J]. Scientia Agricultura Sinica, 2007, 40(12): 2702–2708

    [3] KITTIWORAWAT S, YOUPENSUK S, RERKASEM B. Diversity of arbuscular mycorrhizal fungi inand effect of the soil pH on the symbiosis[J]. Chiang Mai Journal of Science, 2010, 37(3): 517–527

    [4] HELGASON B L, WALLEY F L, GERMIDA J J. Fungal and bacterial abundance in long-term no-till and intensive-till soils of the Northern Great Plains[J]. Soil Science Society of America Journal, 2009, 73(1): 120–127

    [5] TREONIS A M, AUSTIN E E, BUYER J S, et al. Effects of organic amendment and tillage on soil microorganisms and microfauna[J]. Applied Soil Ecology, 2010, 46(1): 103–110

    [6] 曹鵬. 典型生態(tài)系統(tǒng)土壤微生物分布特征及生態(tài)學(xué)機(jī)制[D]. 北京: 中國(guó)科學(xué)院研究生院, 2015 CAO P. Distribution patterns and the ecological mechanisms of soil microorganism in typical ecosystems[D]. Beijing: University of Chinese Academy of Sciences, 2015

    [7] 何玉梅, 張仁陟, 張麗華, 等. 不同耕作措施對(duì)土壤真菌群落結(jié)構(gòu)與生態(tài)特征的影響[J]. 生態(tài)學(xué)報(bào), 2007, 27(1): 113–119 HE Y M, ZHANG R Z, ZHANG L H, et al. Effects of different tillage practices on fungi community structure and ecologic characteristics in loess soils[J]. Acta Ecologica Sinica, 2007, 27(1): 113–119

    [8] YAN F, MCBRATNEY A B, COPELAND L. Functional substrate biodiversity of cultivated and uncultivated A horizons of vertisols in NW New South Wales[J]. Geoderma, 2000, 96(4): 321–343

    [9] 李彤, 王梓廷, 劉露, 等. 保護(hù)性耕作對(duì)西北旱區(qū)土壤微生物空間分布及土壤理化性質(zhì)的影響[J]. 中國(guó)農(nóng)業(yè)科學(xué), 2017, 50(5): 859–870 LI T, WANG Z T, LIU L, et al. Effect of conservation tillage practices on soil microbial spatial distribution and soil physico-chemical properties of the northwest dryland[J]. Scientia Agricultura Sinica, 2017, 50(5): 859–870

    [10] 高煥文. 保護(hù)性耕作技術(shù)與機(jī)具[M]. 北京: 化學(xué)工業(yè)出版社, 2004 GAO H W. Conservation Tillage Technology and Machinery Tools[M]. Beijing: Chemical Industry Press, 2004

    [11] 王碧勝, 蔡典雄, 武雪萍, 等. 長(zhǎng)期保護(hù)性耕作對(duì)土壤有機(jī)碳和玉米產(chǎn)量及水分利用的影響[J]. 植物營(yíng)養(yǎng)與肥料學(xué)報(bào), 2015, 21(6): 1455–1464 WANG B S, CAI D X, WU X P, et al. Effects of long-term conservation tillage on soil organic carbon, maize yield and water utilization[J]. Journal of Plant Nutrition and Fertilizer, 2015, 21(6): 1455–1464

    [12] 張貴云, 張麗萍, 魏明峰, 等. 長(zhǎng)期保護(hù)性耕作對(duì)叢枝菌根真菌多樣性的影響[J]. 中國(guó)生態(tài)農(nóng)業(yè)學(xué)報(bào), 2018, 26(7): 1048–1055 ZHANG G Y, ZHANG L P, WEI M F, et al. Effect of long-term conservation tillage on arbuscular mycorrhizal fungi diversity[J]. Chinese Journal of Eco-Agriculture, 2018, 26(7): 1048–1055

    [13] 遲鳳琴, 劉晶鑫, 匡恩俊, 等. 黑土長(zhǎng)期定位試驗(yàn)原狀土搬遷對(duì)土壤細(xì)菌群落多樣性的影響[J]. 土壤通報(bào), 2015, 46(6): 1420–1427 CHI F Q, LIU J X, KUANG E J, et al. Effects of black soil long-term test undisturbed soil moved on soil bacterial community diversity[J]. Chinese Journal of Soil Science, 2015, 46(6): 1420–1427

    [14] 李橋. 基于高通量測(cè)序技術(shù)下土壤微生物群落結(jié)構(gòu)的研究[D]. 濟(jì)南: 山東師范大學(xué), 2014 LI Q. Research of soil microbial community structure based on high-throughput sequencing technology[D]. Jinan: Shandong Normal University, 2014

    [15] 韓亞飛, 伊文慧, 王文波, 等. 基于高通量測(cè)序技術(shù)的連作楊樹人工林土壤細(xì)菌多樣性研究[J]. 山東大學(xué)學(xué)報(bào): 理學(xué)版, 2014, 49(5): 1–6 HAN Y F, YI W H, WANG W B, et al. Soil bacteria diversity in continuous cropping poplar plantation by high throughput sequencing[J]. Journal of Shandong University: Natural Science, 2014, 49(5): 1–6

    [16] 岳桂東, 高強(qiáng), 羅龍海, 等. 高通量測(cè)序技術(shù)在動(dòng)植物研究領(lǐng)域中的應(yīng)用[J]. 中國(guó)科學(xué): 生命科學(xué), 2012, 42(2): 107–124 YUE G D, GAO Q, LUO L H, et al. The application of high-throughput sequencing technology in plant and animal research[J]. Scientia Sinica Vitae, 2012, 42(2): 107–124

    [17] 牛世全, 龍洋, 李海云, 等. 應(yīng)用Illumina MiSeq高通量測(cè)序技術(shù)分析河西走廊地區(qū)鹽堿土壤微生物多樣性[J]. 微生物學(xué)通報(bào), 2017, 44(9): 2067–2078 NIU S Q, LONG Y, LI H Y, et al. Microbial diversity in saline alkali soil from Hexi Corridor analyzed by Illumina MiSeq high-throughput sequencing system[J]. Microbiology China, 2017, 44(9): 2067–2078

    [18] VERE D. Research into Conservation Tillage for Dryland Cropping in Australia and China[R]. Canberra: The Australian Centre for International Agricultural Research, 2005: 8–46

    [19] 鮑士旦. 土壤農(nóng)化分析[M]. 第3版. 北京: 中國(guó)農(nóng)業(yè)出版社, 2000 BAO S D. Soil and Agricultural Chemistry Analysis[M]. 3rd ed. Beijing: China Agriculture Press, 2000

    [20] PITTA D W, PARMAR N, PATEL A K, et al. Bacterial diversity dynamics associated with different diets and different primer pairs in the rumen of Kankrej cattle[J]. PLoS One, 2014, 9(11): e111710

    [21] 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

    [22] 曹輝, 李燕歌, 周春然, 等. 炭化蘋果枝對(duì)蘋果根區(qū)土壤細(xì)菌和真菌多樣性的影響[J]. 中國(guó)農(nóng)業(yè)科學(xué), 2016, 49(17): 3413–3424 CAO H, LI Y G, ZHOU C R, et al. Effect of carbonized apple branches on bacterial and fungal diversities in apple root-zone soil[J]. Scientia Agricultura Sinica, 2016, 49(17): 3413–3424

    [23] ZHANG C H, LI S F, YANG L, et al. Structural modulation of gut microbiota in life-long calorie-restricted mice[J]. Nature Communications, 2013, 4: 2163

    [24] 戴雅婷, 閆志堅(jiān), 解繼紅, 等. 基于高通量測(cè)序的兩種植被恢復(fù)類型根際土壤細(xì)菌多樣性研究[J]. 土壤學(xué)報(bào), 2017, 54(3): 735–748 DAI Y T, YAN Z J, XIE J H, et al. Soil bacteria diversity in rhizosphere under two types of vegetation restoration based on high throughput sequencing[J]. Acta Pedologica Sinica, 2017, 54(3): 735–748

    [25] 厲桂香, 馬克明. 土壤微生物多樣性海拔格局研究進(jìn)展[J]. 生態(tài)學(xué)報(bào), 2018, 38(5): 1521–1529 LI G X, MA K M. Progress in the study of elevational patterns of soil microbial diversity[J]. Acta Ecologica Sinica, 2018, 38(5): 1521–1529

    [26] 李玉潔, 王慧, 趙建寧, 等. 耕作方式對(duì)農(nóng)田土壤理化因子和生物學(xué)特性的影響[J]. 應(yīng)用生態(tài)學(xué)報(bào), 2015, 26(3): 939–948 LI Y J, WANG H, ZHAO J N, et al. Effects of tillage methods on soil physicochemical properties and biological characteristics in farmland: A review[J]. Chinese Journal of Applied Ecology, 2015, 26(3): 939–948

    [27] DOLAN M S, CLAPP C E, ALLMARAS R R, et al. Soil organic carbon and nitrogen in a Minnesota soil as related to tillage, residue and nitrogen management[J]. Soil and Tillage Research, 2006, 89(2): 221–231

    [28] FILHO C C, LOUREN?O A, DE F GUIMAR?ES M, et al. Aggregate stability under different soil management systems in a red latosol in the state of Parana, Brazil[J]. Soil and Tillage Research, 2002, 65(1): 45–51

    [29] 楊培培, 楊明欣, 董文旭, 等. 保護(hù)性耕作對(duì)土壤養(yǎng)分分布及冬小麥吸收與分配的影響[J]. 中國(guó)生態(tài)農(nóng)業(yè)學(xué)報(bào), 2011, 19(4): 755–759 YANG P P, YANG M X, DONG W X, et al. Effect of conservation tillage on wheat and soil nutrient distribution and absorption[J]. Chinese Journal of Eco-Agriculture, 2011, 19(4): 755–759

    [30] 李友軍, 黃明, 吳金芝, 等. 不同耕作方式對(duì)豫西旱區(qū)坡耕地水肥利用與流失的影響[J]. 水土保持學(xué)報(bào), 2006, 20(2): 42–45 LI Y J, HUANG M, WU J Z, et al. Effects of Different Tillage on utilization and run-off of water and nutrient in sloping farmland of Yuxi dryland area[J]. Journal of Soil and Water Conservation, 2006, 20(2): 42–45

    [31] 戴亮. 耕作方式對(duì)土壤微生物及土壤腐殖質(zhì)影響的研究[D]. 長(zhǎng)春: 吉林農(nóng)業(yè)大學(xué), 2012 DAI L. Tillage practices on soil microbes and soil humus impact[D]. Changchun: Jilin Agricultural University, 2012

    [32] 許菁, 李曉莎, 許姣姣, 等. 長(zhǎng)期保護(hù)性耕作對(duì)麥-玉兩熟農(nóng)田土壤碳氮儲(chǔ)量及固碳固氮潛力的影響[J]. 水土保持學(xué)報(bào), 2015, 29(6): 191–196 XU J, LI X S, XU J J, et al. Effects of long-term conservation tillage on soil organic carbon and nitrogen storages, soil carbon sequestration potential and nitrogen fixation potential in wheat-maize crop system[J]. Journal of Soil and Water Conservation, 2015, 29(6): 191–196

    [33] 王克鵬, 張仁陟, 董博, 等. 長(zhǎng)期保護(hù)性耕作對(duì)黃土高原旱地土壤水分及作物葉水勢(shì)的影響[J]. 生態(tài)學(xué)報(bào), 2014, 34(13): 3752–3761 WANG K P, ZHANG R Z, DONG B, et al. Effect of long-term conservation tillage on soil water regimes and leaf water potential of crops in rainfed areas of the Loess Plateau[J]. Acta Ecologica Sinica, 2014, 34(13): 3752–3761

    [34] 郁鑫, 王旭東. 黃土高原不同耕作方式對(duì)土壤理化性質(zhì)及作物產(chǎn)量的影響[J]. 安徽農(nóng)業(yè)科學(xué), 2018, 46(5): 144–146 YU X, WANG X D. Effects of different tillage methods on soil physical-chemical properties and crop yield in Loess Plateau[J]. Journal of Anhui Agricultural Sciences, 2018, 46(5): 144–146

    [35] 嚴(yán)潔, 鄧良基, 黃劍. 保護(hù)性耕作對(duì)土壤理化性質(zhì)和作物產(chǎn)量的影響[J]. 中國(guó)農(nóng)機(jī)化, 2005, (2): 31–34 YAN J, DENG L J, HUANG J. Effect of conservation tillage on soil physicochemical properties and crop yields[J]. Chinese Agricultural Mechanization, 2005, (2): 31–34

    [36] 李娟. 渭北旱塬玉米田輪耕下土壤質(zhì)量及生產(chǎn)性能響應(yīng)研究[D]. 楊凌: 西北農(nóng)林科技大學(xué), 2017 LI J. Research on soil quality and production capability under rotation tillage on maize field in Weibei highland[D]. Yangling: Northwest A&F University, 2017

    [37] 趙洪利, 李軍, 賈志寬, 等. 不同耕作方式對(duì)黃土高原旱地麥田土壤物理性狀的影響[J]. 干旱地區(qū)農(nóng)業(yè)研究, 2009, 27(3): 17–21 ZHAO H L, LI J, JIA Z K, et al. Effect of different tillages on soil physical properties of dryland wheat field in the Loess Plateau[J]. Agricultural Research in the Arid Areas, 2009, 27(3): 17–21

    [38] DUNBAR J, BARNS S M, TICKNOR L O, et al. Empirical and theoretical bacterial diversity in four Arizona soils[J]. Applied and Environmental Microbiology, 2002, 68(6): 3035–3045

    [39] ROESCH L F W, FULTHORPE R R, RIVA A, et al. Pyrosequencing enumerates and contrasts soil microbial diversity[J]. The ISME Journal, 2007, 1(4): 283–290

    [40] 劉洋, 曾全超, 黃懿梅. 基于454高通量測(cè)序的黃土高原不同喬木林土壤細(xì)菌群落特征[J]. 中國(guó)環(huán)境科學(xué), 2016, 36(11): 3487–3494 LIU Y, ZENG Q C, HUANG Y M. Soil microbial communities by 454prosequencing under different arbor forests on the Loess Plateau[J]. China Environmental Science, 2016, 36(11): 3487–3494

    [41] LIU J J, SUI Y Y, YU Z H, et al. High throughput sequencing analysis of biogeographical distribution of bacterial communities in the black soils of northeast China[J]. Soil Biology and Biochemistry, 2014, 70: 113–122

    [42] DORAN J W. Soil microbial and biochemical changes associated with reduced tillage[J]. Soil Science Society of America Journal, 1980, 44(4): 765–771

    [43] LIU X, LINDEMANN W C, WHITFORD W G, et al. Microbial diversity and activity of disturbed soil in the northern Chihuahuan Desert[J]. Biology and Fertility of Soils, 2000, 32(3): 243–249

    [44] BENDING G D, TURNER M K, JONES J E. Interactions between crop residue and soil organic matter quality and the functional diversity of soil microbial communities[J]. Soil Biology and Biochemistry, 2002, 34(8): 1073–1082

    [45] 王伏偉, 王曉波, 李金才, 等. 施肥及秸稈還田對(duì)砂姜黑土細(xì)菌群落的影響[J]. 中國(guó)生態(tài)農(nóng)業(yè)學(xué)報(bào), 2015, 23(10): 1302–1311 WANG F W, WANG X B, LI J C, et al. Effects of fertilization and straw incorporation on bacterial communities in lime concretion black soil[J]. Chinese Journal of Eco-Agriculture, 2015, 23(10): 1302–1311

    [46] 羅培宇. 輪作條件下長(zhǎng)期施肥對(duì)棕壤微生物群落的影響[D]. 沈陽(yáng): 沈陽(yáng)農(nóng)業(yè)大學(xué), 2014 LUO P Y. Effect of long-term fertilization on microbial community in a brown soil under crop-rotation system[D]. Shenyang: Shenyang Agricultural University, 2014

    Effect of long-term no-tillage with stubble on soil fertility and diversity of prokaryotic microbiome in dryland wheat soils on the Loess Plateau, China*

    ZHANG Guiyun1, LYU Beibei1, ZHANG Liping1, LIU Zhen1, FAN Qiaolan1, WEI Mingfeng1, YAO Zhong1,YUAN Jiawei1, CHAI Yuejin2

    (1. Cotton Research Institute, Shanxi Academy of Agricultural Sciences, Yuncheng 044000, China; 2. Agricultural Machinery Bureau of Linfen City, Linfen 041000, China)

    In dryland regions, conservation tillage has the advantages of reducing soil erosion, maintaining soil structure, improving soil physical and chemical properties, and promoting microbial reproduction, compared with traditional farming methods. Up to now, the effect of long-term conservation tillage on soil prokaryotic microbial community has not been adequately reported for the China’s Loess Plateau region. To understand the effects of long-term conservation tillage on soil prokaryotic microbial community diversity and soil fertility in dryland wheat fields on the Loess Plateau, the Illumina Hiseq 2500 high-throughput sequencing and physicochemical test methods were used to analyze the effects of different tillage measures on soil physical and chemical properties and the diversity of soil prokaryotic microbial community in the wheat fields in Linfen City, Shanxi Province. Tillage measures adopted in the study area for 26 years included no-tillage with stubble (NTS), no-tillage with stubble after subsoiling (SNTS) and traditional tillage (TT1). The results showed that: 1) while soil contents of total nitrogen, available nitrogen, available phosphorus, available potassium, and soil water storage capacity and water content under NTS and SNTS treatments were significantly higher than those under TT1, soil pH under NTS and SNTS was lower than that under TT1. NTS and SNTS treatments reduced soil bulk density in the 0-10 cm soil layer, but increased it in 10-20 cm soil layer compared with TT1 treatment. Meanwhile, SNTS treatment significantly increased soil organic matter content. 2) The results of the Illumina Hiseq high-throughput 16S rRNA gene sequencing in area V4 indicated that: SNTS treatment significantly reduced the relative abundance of Chloroflexi, compared with TT1 treatment. NTS treatment significantly reduced the relative abundance of Verrucomicrobia and Chloroflexi, compared with SNTS treatment. NTS treatment significantly increased the diversity of soil prokaryote community, but not significantly changed the abundance of prokaryote community. SNTS treatment not significantly changed the diversity and richness of prokaryote community. Biomarkers under NTS treatment were higher than under other treatments. The relative abundance of other prokaryotes at the phylum division was not significantly different among three treatments. 3) UPGMA analysis showed that prokaryote community structures of NTS and SNTS were significantly different from those of TT1, and the differences between NTS and SNTS were small. 4) CCA analysis showed that soil pH, contents of organic matter, available nitrogen, available phosphorus and available potassium were important for changes in genetic diversity of soil prokaryote communities. Although NTS and SNTS treatments could change soil prokaryotic microbiome structure greatly, there was still structural similarity among the soil prokaryotic microbiomes among three treatments. To sum up, long-term conservation tillage had obvious positive effects on soil microbial diversity, richness and soil fertility in dryland wheat fields on the Loess Plateau.

    Conservation tillage; No-tillage with stubble; No-tillage with stubble after subsoiling; High-throughput sequencing; Soil prokaryotic microbiome

    , ZHANG Guiyun, E-mail: guiyunzhang@126.com

    Jun. 28, 2018;

    Oct. 25, 2018

    S154.3

    A

    2096-6237(2019)03-0358-11

    10.13930/j.cnki.cjea.180604

    張貴云, 呂貝貝, 張麗萍, 劉珍, 范巧蘭, 魏明峰, 姚眾, 袁嘉瑋, 柴躍進(jìn). 黃土高原旱地麥田26年免耕覆蓋對(duì)土壤肥力及原核微生物群落多樣性的影響[J]. 中國(guó)生態(tài)農(nóng)業(yè)學(xué)報(bào)(中英文), 2019, 27(3): 358-368

    ZHANG G Y, LYU B B, ZHANG L P, LIU Z, FAN Q L, WEI M F, YAO Z, YUAN J W, CHAI Y J. Effect of long-term no-tillage with stubble on soil fertility and diversity of prokaryotic microbiome in dryland wheat soils on the Loess Plateau, China[J]. Chinese Journal of Eco-Agriculture, 2019, 27(3): 358-368

    張貴云, 主要從事植物營(yíng)養(yǎng)學(xué)和菌根學(xué)研究。E-mail: guiyunzhang@126.com

    2018-06-28

    2018-10-25

    * This work was supported by Shanxi Key Research and Development Project (201703D321009-3), the Natural Science Foundation of Shanxi Province (2014011029-3), and the Doctor Fund Project of Shanxi Academy of Agricultural Sciences (YBSJJ1405).

    * 山西省重點(diǎn)研發(fā)計(jì)劃項(xiàng)目(201703D321009-3)、山西省自然科學(xué)基金項(xiàng)目(2014011029-3)和山西省農(nóng)業(yè)科學(xué)院博士基金項(xiàng)目(YBSJJ1405)資助

    猜你喜歡
    原核保護(hù)性耕作
    寧夏保護(hù)性耕作的必要性及技術(shù)模式探討
    保護(hù)性開發(fā)利用高鋁煤炭資源的必要性及相關(guān)建議
    耕作深度對(duì)紫色土坡地旋耕機(jī)耕作侵蝕的影響
    玉米保護(hù)性耕作的技術(shù)要領(lǐng)
    結(jié)核分枝桿菌CFP10、ESAT6蛋白的原核表達(dá)、純化及ELISPOT檢測(cè)方法的建立與應(yīng)用
    草地耕作技術(shù)在澳大利亞的應(yīng)用
    土壤與作物(2015年3期)2015-12-08 00:46:58
    癌癥標(biāo)記蛋白 AGR2的原核表達(dá)及純化
    牛分支桿菌HBHA基因的克隆及原核表達(dá)
    保護(hù)性耕作機(jī)具選型中注意事項(xiàng)
    人巨細(xì)胞病毒pp150-gp52蛋白原核可溶性表達(dá)與IgM捕獲ELISA方法建立和應(yīng)用
    亚洲人成网站在线观看播放| 免费观看a级毛片全部| av网站在线播放免费| 亚洲国产最新在线播放| 免费黄网站久久成人精品| 亚洲在久久综合| 亚洲av欧美aⅴ国产| 国产成人欧美在线观看 | 99久久99久久久精品蜜桃| 国产亚洲一区二区精品| 热99久久久久精品小说推荐| 日韩人妻精品一区2区三区| 中文字幕色久视频| 亚洲国产最新在线播放| 一个人免费看片子| 久久精品国产亚洲av高清一级| 99热网站在线观看| 日韩免费高清中文字幕av| 久久久久久人人人人人| 亚洲男人天堂网一区| 亚洲一区二区三区欧美精品| 多毛熟女@视频| 99久国产av精品国产电影| 国产日韩欧美亚洲二区| 国产成人精品久久久久久| 女人精品久久久久毛片| 亚洲综合精品二区| 亚洲欧美色中文字幕在线| 99精国产麻豆久久婷婷| 可以免费在线观看a视频的电影网站 | 人人妻,人人澡人人爽秒播 | 久久天躁狠狠躁夜夜2o2o | 五月天丁香电影| 日本91视频免费播放| 成人亚洲精品一区在线观看| 少妇人妻久久综合中文| 一本大道久久a久久精品| www.自偷自拍.com| 国产一区二区激情短视频 | av网站免费在线观看视频| 亚洲国产欧美在线一区| 夜夜骑夜夜射夜夜干| 亚洲少妇的诱惑av| 亚洲欧美精品自产自拍| 国产精品三级大全| 建设人人有责人人尽责人人享有的| 国产精品久久久人人做人人爽| 亚洲一码二码三码区别大吗| 欧美另类一区| 少妇 在线观看| 午夜福利一区二区在线看| 伊人亚洲综合成人网| 国产精品免费大片| 亚洲熟女毛片儿| 久久天堂一区二区三区四区| 国产一区二区三区av在线| 9色porny在线观看| 在线亚洲精品国产二区图片欧美| 午夜福利在线免费观看网站| 久久久久视频综合| 少妇 在线观看| 热99久久久久精品小说推荐| 日本色播在线视频| 亚洲婷婷狠狠爱综合网| 国产免费又黄又爽又色| 成年美女黄网站色视频大全免费| √禁漫天堂资源中文www| 亚洲国产精品成人久久小说| 免费在线观看完整版高清| av在线观看视频网站免费| 搡老乐熟女国产| 青春草国产在线视频| 亚洲,欧美精品.| 在线观看免费日韩欧美大片| 少妇猛男粗大的猛烈进出视频| 男人舔女人的私密视频| 看十八女毛片水多多多| 极品人妻少妇av视频| 精品卡一卡二卡四卡免费| 久久精品国产a三级三级三级| 波野结衣二区三区在线| 亚洲美女视频黄频| av电影中文网址| 国产成人系列免费观看| 亚洲国产精品一区二区三区在线| svipshipincom国产片| 丝袜脚勾引网站| 亚洲第一av免费看| 汤姆久久久久久久影院中文字幕| 国精品久久久久久国模美| 男人舔女人的私密视频| 亚洲人成77777在线视频| 操美女的视频在线观看| 欧美精品亚洲一区二区| 老汉色av国产亚洲站长工具| 国产欧美日韩综合在线一区二区| 国语对白做爰xxxⅹ性视频网站| 一个人免费看片子| 汤姆久久久久久久影院中文字幕| 黄片播放在线免费| 日韩制服丝袜自拍偷拍| 美女中出高潮动态图| 亚洲国产精品一区二区三区在线| 观看av在线不卡| 男女高潮啪啪啪动态图| 国产亚洲一区二区精品| 99精国产麻豆久久婷婷| 久久久久精品国产欧美久久久 | 精品一区二区三区av网在线观看 | 亚洲一码二码三码区别大吗| 我的亚洲天堂| 国产亚洲av片在线观看秒播厂| 自线自在国产av| 亚洲成色77777| 最新在线观看一区二区三区 | 黄色视频在线播放观看不卡| 日韩,欧美,国产一区二区三区| 国产欧美日韩综合在线一区二区| 国产片特级美女逼逼视频| 久久精品aⅴ一区二区三区四区| 午夜日本视频在线| 91精品国产国语对白视频| 人人妻,人人澡人人爽秒播 | 欧美日韩一级在线毛片| 色播在线永久视频| 两个人免费观看高清视频| 汤姆久久久久久久影院中文字幕| 亚洲国产欧美网| 欧美精品高潮呻吟av久久| 久久久久久久久久久久大奶| av在线播放精品| 十八禁高潮呻吟视频| 国产免费福利视频在线观看| 成人三级做爰电影| 女的被弄到高潮叫床怎么办| 亚洲七黄色美女视频| 国产精品 欧美亚洲| 日韩中文字幕欧美一区二区 | 欧美日韩福利视频一区二区| 嫩草影院入口| 成年动漫av网址| av在线观看视频网站免费| 啦啦啦视频在线资源免费观看| 97在线人人人人妻| 欧美日韩成人在线一区二区| 国产免费视频播放在线视频| 日韩 欧美 亚洲 中文字幕| 国产精品国产三级专区第一集| 久久99热这里只频精品6学生| 亚洲美女搞黄在线观看| 欧美精品一区二区大全| 少妇猛男粗大的猛烈进出视频| 国产精品久久久av美女十八| 亚洲精品久久久久久婷婷小说| 一区二区三区四区激情视频| 国产精品香港三级国产av潘金莲 | 久久综合国产亚洲精品| 91成人精品电影| 国产精品久久久久久人妻精品电影 | 日韩制服丝袜自拍偷拍| 亚洲av中文av极速乱| 日韩成人av中文字幕在线观看| 亚洲图色成人| 汤姆久久久久久久影院中文字幕| 亚洲精品国产av成人精品| 80岁老熟妇乱子伦牲交| 操美女的视频在线观看| 超碰97精品在线观看| 亚洲人成77777在线视频| 日本猛色少妇xxxxx猛交久久| 国产精品久久久久久精品古装| 欧美最新免费一区二区三区| 精品第一国产精品| 国产色婷婷99| 国产极品粉嫩免费观看在线| 亚洲欧美精品综合一区二区三区| 亚洲中文av在线| 成人国语在线视频| 国产成人欧美| 一级片免费观看大全| 久久精品久久精品一区二区三区| 久久婷婷青草| 精品一区二区免费观看| 亚洲精品第二区| 日韩熟女老妇一区二区性免费视频| 国产男人的电影天堂91| 九九爱精品视频在线观看| 人成视频在线观看免费观看| 久久青草综合色| 天天躁狠狠躁夜夜躁狠狠躁| 一区在线观看完整版| 无限看片的www在线观看| 免费看不卡的av| 午夜老司机福利片| 激情视频va一区二区三区| 90打野战视频偷拍视频| 9色porny在线观看| avwww免费| 哪个播放器可以免费观看大片| 最近手机中文字幕大全| 久久鲁丝午夜福利片| 成人国语在线视频| 999久久久国产精品视频| 亚洲国产精品一区二区三区在线| 国产成人a∨麻豆精品| 精品视频人人做人人爽| av片东京热男人的天堂| 亚洲国产欧美日韩在线播放| 国产福利在线免费观看视频| 久久精品久久精品一区二区三区| 久久毛片免费看一区二区三区| 又大又黄又爽视频免费| 亚洲三区欧美一区| 男女下面插进去视频免费观看| 久久青草综合色| 男人添女人高潮全过程视频| 极品人妻少妇av视频| 侵犯人妻中文字幕一二三四区| 精品酒店卫生间| 午夜激情久久久久久久| 亚洲少妇的诱惑av| 欧美日韩一级在线毛片| 天堂俺去俺来也www色官网| 欧美日韩一区二区视频在线观看视频在线| 看免费av毛片| 男女无遮挡免费网站观看| 操出白浆在线播放| 日韩中文字幕欧美一区二区 | 大码成人一级视频| 精品亚洲乱码少妇综合久久| 亚洲天堂av无毛| 五月天丁香电影| 国产一区二区激情短视频 | 一本—道久久a久久精品蜜桃钙片| 建设人人有责人人尽责人人享有的| 精品国产乱码久久久久久男人| 精品一区二区三区av网在线观看 | 国产麻豆69| 97精品久久久久久久久久精品| 国产精品一区二区在线观看99| 国产片内射在线| 精品亚洲乱码少妇综合久久| 中文字幕人妻熟女乱码| 麻豆乱淫一区二区| 午夜福利网站1000一区二区三区| 美国免费a级毛片| 国产精品二区激情视频| 男女高潮啪啪啪动态图| 日韩电影二区| 一本色道久久久久久精品综合| 亚洲精品国产区一区二| 日韩大码丰满熟妇| 国产精品99久久99久久久不卡 | 国产 精品1| 天天添夜夜摸| av国产精品久久久久影院| 久热这里只有精品99| 欧美精品一区二区大全| 又黄又粗又硬又大视频| 国产在线免费精品| 天美传媒精品一区二区| 热99久久久久精品小说推荐| 亚洲国产中文字幕在线视频| 巨乳人妻的诱惑在线观看| 日韩精品有码人妻一区| 亚洲av综合色区一区| 在线看a的网站| 91精品伊人久久大香线蕉| 尾随美女入室| 成人影院久久| 欧美变态另类bdsm刘玥| 日韩中文字幕视频在线看片| 女人高潮潮喷娇喘18禁视频| 日韩精品有码人妻一区| 久久久久久久大尺度免费视频| 欧美精品人与动牲交sv欧美| 女性被躁到高潮视频| 成人亚洲精品一区在线观看| 日日爽夜夜爽网站| 不卡视频在线观看欧美| 国产精品女同一区二区软件| 欧美变态另类bdsm刘玥| 欧美国产精品一级二级三级| 久久婷婷青草| 久久久国产一区二区| 亚洲av在线观看美女高潮| 国产一区二区三区av在线| 国产黄色视频一区二区在线观看| 欧美精品一区二区大全| a级毛片黄视频| 人人妻人人澡人人爽人人夜夜| 久久综合国产亚洲精品| 久久亚洲国产成人精品v| 一边亲一边摸免费视频| 亚洲第一av免费看| 成年女人毛片免费观看观看9 | 亚洲国产看品久久| 亚洲国产欧美在线一区| 在线观看一区二区三区激情| 中国三级夫妇交换| 97在线人人人人妻| 免费不卡黄色视频| 伊人久久国产一区二区| 国产精品秋霞免费鲁丝片| 亚洲精品国产区一区二| 香蕉国产在线看| 亚洲欧美日韩另类电影网站| 母亲3免费完整高清在线观看| 天天躁日日躁夜夜躁夜夜| 岛国毛片在线播放| 啦啦啦视频在线资源免费观看| 国产成人一区二区在线| 国产一区二区 视频在线| 欧美黑人精品巨大| 午夜久久久在线观看| 少妇精品久久久久久久| 十分钟在线观看高清视频www| 自线自在国产av| 一区在线观看完整版| 亚洲三区欧美一区| 国产精品成人在线| 国产精品一二三区在线看| 国产精品久久久久久人妻精品电影 | 黄片小视频在线播放| 国产成人啪精品午夜网站| 国产精品久久久av美女十八| 在线观看三级黄色| 国产免费又黄又爽又色| 大片电影免费在线观看免费| 久久狼人影院| 亚洲精品国产一区二区精华液| 无遮挡黄片免费观看| 精品久久久久久电影网| 伊人久久国产一区二区| 国产97色在线日韩免费| 精品亚洲成国产av| 十分钟在线观看高清视频www| 老熟女久久久| 七月丁香在线播放| 一个人免费看片子| 大码成人一级视频| 久久精品熟女亚洲av麻豆精品| 亚洲成国产人片在线观看| 女人高潮潮喷娇喘18禁视频| 成人国产av品久久久| 精品国产乱码久久久久久小说| 婷婷色av中文字幕| 最黄视频免费看| 国产一级毛片在线| 99热国产这里只有精品6| 国产一级毛片在线| 中文字幕人妻熟女乱码| 免费黄网站久久成人精品| 欧美日韩一级在线毛片| 嫩草影视91久久| 免费在线观看黄色视频的| 如何舔出高潮| 日本一区二区免费在线视频| 黑人巨大精品欧美一区二区蜜桃| 久久这里只有精品19| 日韩制服丝袜自拍偷拍| av片东京热男人的天堂| 少妇精品久久久久久久| 午夜激情久久久久久久| 成人午夜精彩视频在线观看| 国产片内射在线| 国产午夜精品一二区理论片| 国产亚洲最大av| av视频免费观看在线观看| 午夜久久久在线观看| h视频一区二区三区| www.精华液| 久热这里只有精品99| 多毛熟女@视频| 久久久精品免费免费高清| √禁漫天堂资源中文www| 亚洲精品久久成人aⅴ小说| 午夜久久久在线观看| 一区二区日韩欧美中文字幕| 中国国产av一级| 亚洲,欧美精品.| 亚洲情色 制服丝袜| 欧美人与善性xxx| 久久综合国产亚洲精品| 九草在线视频观看| av国产久精品久网站免费入址| 成人亚洲精品一区在线观看| 国语对白做爰xxxⅹ性视频网站| 香蕉丝袜av| 高清欧美精品videossex| 51午夜福利影视在线观看| 老司机影院成人| 人成视频在线观看免费观看| 99热全是精品| 亚洲一卡2卡3卡4卡5卡精品中文| 韩国精品一区二区三区| 爱豆传媒免费全集在线观看| 美女中出高潮动态图| 男女边吃奶边做爰视频| 欧美激情 高清一区二区三区| 久久综合国产亚洲精品| av有码第一页| 亚洲精品国产av成人精品| 亚洲欧美一区二区三区黑人| 国产精品一区二区精品视频观看| 91国产中文字幕| 国产成人欧美| 国产麻豆69| 在线观看人妻少妇| 国产亚洲午夜精品一区二区久久| 欧美激情 高清一区二区三区| 观看美女的网站| 午夜福利,免费看| 久久99热这里只频精品6学生| 黄色一级大片看看| 亚洲综合精品二区| 久久久久久久久久久久大奶| 久久久精品区二区三区| 视频区图区小说| 精品亚洲乱码少妇综合久久| 色综合欧美亚洲国产小说| 精品人妻在线不人妻| 久久影院123| 国产在线视频一区二区| 咕卡用的链子| 欧美在线一区亚洲| 一本大道久久a久久精品| 成年美女黄网站色视频大全免费| 亚洲欧美精品综合一区二区三区| 国产精品无大码| 亚洲国产成人一精品久久久| 最近最新中文字幕免费大全7| 免费少妇av软件| 国产亚洲午夜精品一区二区久久| 高清视频免费观看一区二区| 热re99久久国产66热| 欧美激情高清一区二区三区 | 波多野结衣一区麻豆| 久久人人97超碰香蕉20202| 一级毛片 在线播放| 美女中出高潮动态图| 欧美精品av麻豆av| 午夜福利网站1000一区二区三区| 亚洲欧美激情在线| 看非洲黑人一级黄片| 涩涩av久久男人的天堂| 男女边吃奶边做爰视频| 亚洲精品美女久久久久99蜜臀 | 精品亚洲乱码少妇综合久久| 国产亚洲欧美精品永久| 欧美黄色片欧美黄色片| 午夜日韩欧美国产| 欧美精品一区二区大全| 亚洲成人一二三区av| 日韩不卡一区二区三区视频在线| 一级爰片在线观看| 久久影院123| 国产成人精品福利久久| 日本欧美国产在线视频| 免费高清在线观看视频在线观看| 免费黄网站久久成人精品| 国产成人欧美在线观看 | 亚洲激情五月婷婷啪啪| 如何舔出高潮| 18禁裸乳无遮挡动漫免费视频| 一区二区日韩欧美中文字幕| 最近中文字幕高清免费大全6| 久久女婷五月综合色啪小说| 国产精品久久久久久精品电影小说| 国产黄色免费在线视频| 久久久久久久国产电影| 亚洲成国产人片在线观看| av视频免费观看在线观看| 久久99热这里只频精品6学生| 国产av精品麻豆| 一边摸一边抽搐一进一出视频| 在线观看免费日韩欧美大片| 乱人伦中国视频| 国产亚洲欧美精品永久| 亚洲在久久综合| 咕卡用的链子| 亚洲一区二区三区欧美精品| 熟妇人妻不卡中文字幕| 熟女少妇亚洲综合色aaa.| 男女床上黄色一级片免费看| 欧美日韩国产mv在线观看视频| 欧美在线一区亚洲| 飞空精品影院首页| 狂野欧美激情性xxxx| 国产成人午夜福利电影在线观看| 美女福利国产在线| 伊人亚洲综合成人网| 老汉色∧v一级毛片| 99热网站在线观看| 国产又色又爽无遮挡免| 中文字幕高清在线视频| 欧美精品人与动牲交sv欧美| 女人被躁到高潮嗷嗷叫费观| 777米奇影视久久| 欧美老熟妇乱子伦牲交| 99九九在线精品视频| a级毛片在线看网站| 亚洲国产欧美日韩在线播放| 大码成人一级视频| 精品少妇黑人巨大在线播放| 国语对白做爰xxxⅹ性视频网站| 麻豆av在线久日| 在线亚洲精品国产二区图片欧美| 精品午夜福利在线看| 午夜免费男女啪啪视频观看| 国产成人免费无遮挡视频| 欧美xxⅹ黑人| 激情视频va一区二区三区| av片东京热男人的天堂| 国产午夜精品一二区理论片| 国产1区2区3区精品| 亚洲av日韩在线播放| 精品国产超薄肉色丝袜足j| 日韩制服骚丝袜av| 在线观看免费视频网站a站| 午夜免费鲁丝| 大话2 男鬼变身卡| 久久久久久免费高清国产稀缺| 国产熟女午夜一区二区三区| 婷婷色av中文字幕| 亚洲色图 男人天堂 中文字幕| 妹子高潮喷水视频| 午夜av观看不卡| 国产亚洲一区二区精品| 亚洲激情五月婷婷啪啪| 久久人人爽人人片av| 男女边摸边吃奶| 水蜜桃什么品种好| 亚洲av国产av综合av卡| 狂野欧美激情性bbbbbb| 日韩av在线免费看完整版不卡| 亚洲欧美一区二区三区国产| 美女大奶头黄色视频| 黄色视频在线播放观看不卡| 老鸭窝网址在线观看| 少妇人妻精品综合一区二区| 麻豆av在线久日| 韩国av在线不卡| 亚洲色图 男人天堂 中文字幕| 韩国高清视频一区二区三区| 国产 一区精品| 亚洲精品中文字幕在线视频| 国产精品无大码| 欧美成人精品欧美一级黄| 日本午夜av视频| 狠狠精品人妻久久久久久综合| 视频区图区小说| 日日啪夜夜爽| 色吧在线观看| 一区福利在线观看| 国产av国产精品国产| 天堂中文最新版在线下载| 妹子高潮喷水视频| 国产一区亚洲一区在线观看| 久久精品国产亚洲av涩爱| 美国免费a级毛片| 18在线观看网站| 国产爽快片一区二区三区| 熟女少妇亚洲综合色aaa.| 黄色一级大片看看| 老司机在亚洲福利影院| 美女主播在线视频| 亚洲国产看品久久| 看免费av毛片| 国产成人精品久久二区二区91 | 交换朋友夫妻互换小说| 亚洲美女搞黄在线观看| 久久午夜综合久久蜜桃| 免费观看人在逋| 欧美亚洲日本最大视频资源| 一边亲一边摸免费视频| 人人妻人人澡人人爽人人夜夜| 狂野欧美激情性bbbbbb| 极品人妻少妇av视频| 亚洲欧美成人精品一区二区| 咕卡用的链子| 国产精品成人在线| av在线老鸭窝| 亚洲精品久久午夜乱码| 亚洲av日韩在线播放| 男人舔女人的私密视频| av一本久久久久| 成人亚洲欧美一区二区av| 青春草国产在线视频| 各种免费的搞黄视频| 不卡视频在线观看欧美| 男人爽女人下面视频在线观看| 操出白浆在线播放| 1024香蕉在线观看| 无限看片的www在线观看| 男人操女人黄网站| 一区二区日韩欧美中文字幕| 亚洲国产欧美网| 久久99精品国语久久久| 亚洲国产毛片av蜜桃av| 啦啦啦视频在线资源免费观看| 男女边摸边吃奶| 亚洲国产毛片av蜜桃av| 啦啦啦视频在线资源免费观看| 国产极品天堂在线| 国产xxxxx性猛交| 嫩草影院入口| 欧美精品一区二区免费开放| 精品视频人人做人人爽| 日韩欧美精品免费久久| 国产熟女午夜一区二区三区| av视频免费观看在线观看| 我要看黄色一级片免费的| 日本一区二区免费在线视频| 欧美老熟妇乱子伦牲交|