不同耕作和秸稈還田下褐土真菌群落變化特征

代紅翠，張慧，薛艷芳，高英波，錢欣，趙海軍，成浩，李宗新，劉開昌

不同耕作和秸稈還田下褐土真菌群落變化特征

代紅翠1,3，張慧2,3，薛艷芳2,3，高英波2,3，錢欣2,3，趙海軍3，成浩4，李宗新2,3，劉開昌1,3

（1山東省農業科學院作物研究所，濟南 250100；2山東省農業科學院玉米研究所，濟南 250100；3小麥玉米國家工程實驗室，濟南 250100；4青島農業大學農學院，山東青島 266109）

【】研究華北平原地區小麥-玉米周年復種模式下不同耕作和秸稈還田方式對土壤真菌群落結構及功能的影響，并探索農田土壤肥力對耕作和秸稈還田方式響應差異的生物學機制，為優化耕作與秸稈還田方式和提高農田土壤肥力提供理論依據。本研究以華北平原小麥-玉米周年復種農田土壤為研究對象，采用Miseq高通量測序技術，結合FUNGuild真菌功能預測工具，分析3種耕作方式（免耕、深耕與旋耕）與兩種秸稈還田方式（麥秸單季還田與小麥-玉米秸稈雙季還田）定位試驗條件下，小麥成熟期土壤真菌群落結構與功能的差異，結合土壤理化性質，進一步探究農田土壤真菌群落結構及功能變化的環境驅動因子。與免耕雙季還田相比，深耕秸稈雙季還田與深耕單季還田0—10 cm耕層土壤有機碳含量分別降低35.04%和44.30%；免耕秸稈單季還田10—20 cm耕層土壤中堿解氮含量顯著低于其他處理。農田土壤0—10 cm土層真菌主要包含子囊菌門(Ascomycota)、擔子菌門(Basidiomycota)和壺菌門(Chytridiomycota)，相對豐度分別為68.98%、16.96%和1.62%；10—20 cm土層真菌主要包含子囊菌門(Ascomycota)、擔子菌門（Basidiomycota）、壺菌門(Chytridiomycota)和球囊菌門(Glomeromycota)，相對豐度分別為68.44%、15.52%、1.51%和1.23%。不同處理土壤真菌群落結構存在差異，與秸稈單季還田相比，秸稈雙季還田0—10 cm和10—20 cm土層擔子菌門(Basidiomycota)相對豐度分別提高了50.07%和29.08%。進一步分析土壤群落結構發生變化的原因，結果顯示0—10 cm土層土壤真菌群落多元回歸樹第一次分割以土壤有機碳為節點，其閾值為11.17 g·kg-1，旋耕雙季還田和免耕雙季還田與其他處理分離；10—20 cm土層土壤真菌群落多元回歸樹第一次分割以堿解氮為節點，其閾值為6.52 mg·kg-1，免耕單季還田與其他處理分離。從營養類型看，各處理0—10 cm（26.84%）和10—20 cm（23.91%）土層土壤真菌均以病理營養型為主；與免耕相比，深耕和旋耕處理0—10 cm土層病理營養型真菌相對豐度分別顯著降低25.16%和16.45%，且以深耕單季還田最低。病理營養型真菌相對豐度與土壤有機碳、可溶性有機碳、全氮、堿解氮和速效鉀含量均呈顯著正相關關系。不同耕作與秸稈還田方式改變了農田土壤真菌群落結構及功能，土壤有機碳和堿解氮的改變是影響真菌群落組成的主要因素；深耕能夠降低秸稈還田后土壤病理營養型真菌相對豐度，利于保持農田土壤生態系統健康。

耕作方式；秸稈還田；土壤真菌群落；高通量測序；多元回歸樹

0 引言

【研究意義】耕地質量下降是中國農業可持續發展的主要限制性因素[1]。在以小麥-玉米周年種植制度為主的華北平原，長期的集約化種植使耕層淺薄化、土壤緊實化、養分利用率較低等一系列土壤健康問題尤為突出[2]。這可能是由傳統耕作措施（如翻耕、移除秸稈、過度施用化肥等）使表土暴露，破壞了土壤結構，引起耕地質量下降[3-4]。目前，耕作方式優化與秸稈還田被當作是現階段該地區農田土壤肥力提升的重要措施[5]。一般認為，不同耕作、施肥和管理方式等可能引起土壤微生物群落和功能的變化。土壤微生物在土壤生態系統中起著關鍵作用，80%—90%的土壤相關功能由微生物調節，包括土壤養分轉化，土壤系統的穩定性及抗干擾能力等[6]。真菌構成了大部分微生物生物量，具有促進土壤有機質循環，并且參與毒物降解和作物病害發生等重要的土壤生態過程[7]，是維持生態系統功能的基礎和生態系統健康的指示物[8]。因此，研究華北平原小麥玉米周年種植制度下，耕作措施和秸稈還田方式下真菌群落、功能變化及其驅動因素對提升農田土壤肥力具有重要意義?！厩叭搜芯窟M展】耕作方式作為重要的農田管理措施影響農業生態環境[9]，合理的耕作與秸稈還田方式可以有效地增加土壤有機質含量，為土壤微生物的生長與繁殖提供豐富的養分資源[10]。WANG等[11]在西北土中研究發現，耕作方式通過影響土壤顆粒機械組成和養分含量，進而影響土壤子囊菌門和擔子菌門的分布水平。SUN等[12]通過磷脂脂肪酸法（PLFA）研究認為在東北黑土中免耕比翻耕更有利于增加表層土壤微生物生物量，且在表層形成以真菌為優勢種群的群落，利于真菌的生長。郭梨錦等[13]在湖北小麥-水稻輪作農田中免耕顯著降低真菌和革蘭氏陰性菌的生物量以及真菌/細菌，長期免耕和秸稈還田有利于提高土壤微生物多樣性。可見，耕作與秸稈還田方式影響農田土壤真菌群落生物量與群落組成。目前，關于農田不同耕作與秸稈還田方式對土壤真菌群落結構與功能的影響研究并不完善。應用Illumina MiSeq高通量測序可以在種水平研究真菌群落組成，但很難解讀真菌的功能[14]。此外，基于PLFA法可以對主導作用的功能型真菌進行分類比較，但更側重微生物群落的總體分析，總體種類較少，數據量相對較小[15]。【本研究切入點】合理的耕作與秸稈還田方式能夠改善土壤水、肥、氣、熱等條件，促進土壤微生物生長，是改善土壤肥力狀況的重要措施。然而，土壤真菌群落及功能對不同耕作和秸稈還田方式的響應差異，尤其是不同耕作和秸稈還田方式影響了哪些關鍵養分指標，進而調控了土壤真菌群落與功能，亟需進行深入研究。本研究采用高通量測序法對土壤真菌群落多樣性進行比較，并結合FUNGuild工具預測真菌功能營養型，系統比較不同耕作與秸稈還田方式對農田土壤真菌群落和功能的影響研究?！緮M解決的關鍵問題】本文基于華北平原小麥-玉米周年復種多年定位試驗，采用Miseq高通量測序，基于FUNGuild對真菌序列進行功能分類，研究土壤真菌群落結構及功能對不同耕作和秸稈還田方式的響應差異，并結合土壤理化性狀揭示其驅動因素，以期闡明不同耕作和秸稈還田方式影響華北農田土壤肥力的生物學機制，為華北平原農田土壤地力提升提供科學依據。

1 材料與方法

1.1 試驗區概況

定位試驗設于山東省農業科學院玉米研究所龍山試驗基地（117°32′E，36°43′N）。該試驗基地位于華北平原典型小麥玉米一年兩熟輪作區，土壤類型為褐土，屬于典型溫帶大陸性季風氣候，年均日照時數2 647.6 h，年均氣溫12.8℃，年均降雨量600.8 mm，無霜期209 d。試驗地土壤有機質為14.81 g·kg-1，全氮為0.85 g·kg-1，有效磷為19.26 mg kg-1，速效鉀為46.37 mg kg-1，pH為7.64。

自2012年10月，設置大田裂區試驗，主因素為小麥播前免耕（NT）、深耕（DT）和旋耕（RT）3種耕作方式（玉米均為貼茬免耕直播）；副因素為小麥-玉米秸稈雙季秸稈還田（DS）和小麥秸稈單季還田（SS），共6個處理（表1）。每個處理3次重復，共18個小區，小區面積為6 m×45 m =270 m2。小麥品種為濟麥22，播量172.5 kg·hm-2，播種方式為寬幅精播；玉米品種為魯單9066，行距60 cm，種植密度為75 000株/hm2。在小麥播種前，施用600 kg·hm-2復合肥（N﹕P2O5﹕K2O為17﹕17﹕17）作基肥，小麥拔節期追施尿素225 kg·hm-2。在玉米播種前，施用600 kg·hm-2復合肥（N﹕P2O5﹕K2O為17﹕17﹕17）作基肥，在玉米大口期追施尿素225 kg·hm-2。

1.2 樣品采集

于2017年6月小麥收獲期，小區“S”形5點采集0—10 cm 和10—20 cm 耕層土壤樣品。同小區同土層樣品混為一個土樣，用冰盒保存并迅速帶回實驗室。剔去植物根和小石塊等雜質后，部分樣品凍存于-80℃冰箱中，用于土壤微生物分析，部分鮮土樣用于測定土壤可溶性有機碳（DOC），部分樣品自然風干后過0.25 mm和1 mm篩測定土壤養分指標，剩余樣品-80℃冷藏測定土壤真菌群落。

1.3 測定方法

采用重鉻酸鉀-容量法測定土壤有機碳（TOC），凱氏定氮法測定全氮（TN），堿解氮（AN）采用堿解擴散法，有效磷（AP）采用鉬銻抗分光光度法，速效鉀（AK）采用火焰光度法。以上測定方法均參考《土壤農化分析》[16]。pH測定采用電位法，水土比為5﹕1。

表1 不同耕作與秸稈還田方式處理

土壤真菌群落的18S rRNA測序土壤DNA提取使用DNA提取試劑盒（OMEGA，美國），提取方法參照說明書。利用1%瓊脂糖凝膠電泳檢測抽提基因組DNA的純度和完整性。使用NanoDrop ND-2000分光光度計（NanoDrop，ND2000，Thermo Scientific，Wilmington）測定所提DNA的濃度與純度。選用真菌18S通用引物0817F:5'-TTAGCATGGAATAATRRAA TAGGA-3'，1196R: 5'-TCTGGACCTGGTGAGTTTCC -3'進行擴增[17]，選用TransGen AP221-02試劑盒，擴增體系為20 μL，包括4 μL的5×FastPfu Buffer，2 μL的2.5 mmol·L-1dNTPs，引物（5 μmol·L-1）各0.8 μL，0.4 μL的FastPfu Polymerase和10 ng模板DNA。 擴增條件為95℃解鏈2 min，95℃ 30 s，55℃ 30 s, 72℃ 30 s，25個循環，再72℃延伸5 min，于ABI GeneAmp? 9700型PCR儀上進行。每個樣本3次重復，將同一樣本的PCR產物混合后用2%瓊脂糖凝膠電泳檢測，回收擴增子，依照AxyPrep DNA Gel Extraction Kit說明書進行純化（Axygen Biosciences, Union City, U.S.），再用QuantiFluorTM-ST（Promega, U.S.）進行定量。采用Illumina MiSeq測序平臺標準流程進行雙端進行測序，由上海美吉生物醫藥科技有限公司提供技術支持。

1.4 數據統計

經過QIIME（v1.8.0）軟件過濾、拼接、去除嵌合體后[18-19]，聚類為用于物種分類的OTU（Operational Taxonomic Units），并將所有樣品進行抽平。采用RDP classifier貝葉斯算法對97%相似水平的OTU代表序列進行分類學分析[20]，置信度閾值為0.7，選用Silva（Release128 http://www.arb-silva.de）數據庫[21]，得到分類學信息。利用Mothur軟件（Version 1.31.2）進行α多樣性分析[22]。基于Unifrace距離進行群落結構分析[23]。

根據“1SE”準則，最小交叉驗證相對誤差值來確定分類結果，能很好地反映異質環境下群落的結構特征[24]，使用R 3.5.1軟件運行mvpart和MVPARTwrap程序包運行并繪制多元回歸樹。使用Python 3.7進行真菌功能類群FUNGuild數據庫[25]比對。不同處理方差分析（ANOVA）與多重比較（Duncan’s test）用SPSS 20.0軟件完成。

2 結果

2.1 耕作方式對土壤理化性質的影響

不同耕作方式顯著影響土壤理化性質（表2）。方差分析結果表明，耕作方式顯著影響農田0—10 cm（＜0.01）及10—20 cm土層土壤有機碳（TOC）（＜0.001），0—10 cm（＜0.01）及10—20 cm可溶性有機碳（DOC）（＜0.05）含量。土壤0—10 cm土層，與免耕雙季秸稈還田相比，深耕雙季秸稈還田TOC與DOC分別降低35.04%和23.89%（＜0.05）；深耕單季秸稈還田TOC與DOC比免耕單季秸稈還田分別降低26.09%和33.78%（＜0.05）。土壤10—20 cm土層，深耕雙季秸稈還田與旋耕雙季秸稈還田TOC含量比免耕雙季秸稈還田分別提高21.64%與38.21%（＜0.05）。秸稈還田方式顯著影響土壤0—10 cm及10—20 cm土層TOC（＜0.01；＜0.001），DOC（＜0.01；＜0.05）與TN（＜0.05；＜0.001）含量。免耕處理中，單季秸稈還田比雙季秸稈還田土壤0—10 cm土層TOC、DOC和TN分別顯著降低24.64%，21.61%和12.61%（＜0.05）。

2.2 土壤真菌多樣性

利用Illumina MiSeq平臺對真菌18S rRNA基因測序結果顯示，36個土壤樣本共獲得1 372 440條有效序列，378個OTU，各樣本的覆蓋率指數（good's coverage）顯示各樣本的覆蓋率指數顯示0—10及10—20 cm土層土壤OTUs涵蓋了土壤中99.9%的真菌（表3）。測序深度能夠比較真實地反映土壤樣本的真菌群落，代表真菌群落多樣性。將數據抽平，每個樣品得到23295條高質量序列（表3）。

OTU數目與α多樣性指數表示物種的豐富度，即一個群落或生境內物種的復雜程度，值越高表明群落內物種數目越多。方差分析結果表明，耕作方式顯著影響0—10 cm和10—20 cm土層真菌OTU數目，而秸稈還田方式顯著影響0—10 cm土層真菌香農指數。

本研究中，圖1顯示農田土壤0—10 cm土層真菌主要包含子囊菌門（Ascomycota）、擔子菌門（Basidiomycota）和壺菌門（Chytridiomycota），相對豐度分別為68.98%、16.96%和1.62%，約占87.55%，未被分類到已知真菌（Unclassified）的序列占9.45%；土壤10—20 cm土層真菌主要包含子囊菌門（Ascomycota）、擔子菌門（Basidiomycota）、壺菌門（Chytridiomycota）和球囊菌門（Glomeromycota），相對豐度分別為68.44%、15.52%、1.51%和1.23%，約占86.71%，未被分類到已知真菌的序列占10.95%。其中，雙季秸稈還田（DS）處理比單季秸稈還田（SS）處理在0—10 cm和10—20 cm土層Basidiomycota相對豐度分別顯著提高50.07%和29.08%（圖2，＜0.05）。β多樣性用來比較不同環境中真菌群落結構的差異性，UniFrac距離用于估計不同處理間的差異。在0—10 cm土層，PCoA縱坐標圖顯示PC1軸和PC2軸分別可以解釋37.29%和23.11%群落組成差異（圖3-a）。相同處理的樣本聚集在一起, 這表明來自相同樣本的相似性很高但不同處理的樣品之間存在差異。秸稈單季還田中，免耕、深耕和旋耕處理之間的簇分離，進一步說明不同處理土壤真菌群落之間存在差異。相反，在10—20 cm土層，除NTS與DTD外，不同處理樣本彼此靠近，說明處理之間真菌群落結構相似（圖3-b）?？梢? 不同耕作和秸稈還田制度對農田不同土層土壤真菌群落的影響存在差異。

表2 不同耕作和秸稈還田方式下土壤化學性質

數據為平均值±標準差（n=3），同列數值后不同小寫字母表示差異顯著（＜0.05）。*表示＜0.05，**表示＜0.01，***表示＜0.001，ns 表示差異不顯著。NTD，免耕+雙季秸稈還田；NTS，免耕+小麥單季秸稈還田；DTD，深耕+雙季秸稈還田；DTS，深耕+小麥單季秸稈還田；RTD，旋耕+雙季秸稈還田；RTS，旋耕+小麥單季秸稈還田。下同

Values are means ± standard deviation (n=3). Different letters mean significant differences (＜0.05). * indicates a significant difference between treatments at＜0.05; ** indicates a significant difference between treatments at＜0.01; *** indicates a significant difference between treatments at＜0.001. no tillage+double-season returning (NTD), no tillage+single-season returning (NTS), deep tillage+double-season returning (DTD), deep tillage+single-season returning (DTS), rotary tillage+double-season returning (RTD) and rotary tillage+single-season returning (RTS). The same as below

表3 不同耕作和秸稈還田方式下土壤真菌高通量測序及群落α多樣性指數

2.3 土壤真菌功能類群

采用FUNGuild預測農田土壤真菌群落的營養型，本研究結果鑒定為病理營養型（Pathotroph）、腐生營養型（Saprotroph）和共生營養型（Symbiotroph）以及其他無法鑒定營養型的種群。從營養類型看，華北平原小麥玉米農田土壤真菌以病理營養型為主（圖4），與NT相比，DT和RT處理分別顯著降低25.16%和16.45%（圖4-a，＜0.05），且以DTS處理最低。RT處理腐生營養型真菌豐度顯著高于DT（0—10 cm）和NT（10—20 cm）（圖4-b，＜0.05），說明秸稈還田后旋耕更有利于腐生營養型真菌的生長。共生營養型真菌在0—10 cm土層SS處理中含量顯著高于DS（圖4-c，＜0.05），說明連續多年集約化的雙季秸稈還田不利于共生營養型真菌在農田耕層土壤中的富集。

2.4 土壤真菌群落結構與環境因子的關系

多元回歸樹（MRT）是一種研究物種分布與環境因子之間關系的數量分類方法。將微生物群落進行多元回歸分析，以明確驅動土壤真菌群落的主要因子。圖5顯示不同處理0—10 cm土層土壤真菌群落多元回歸樹第一次分割以有機碳為節點，將所有處理分為兩組，將RTD和NTD處理與其他處理分解，該次分組能解釋群落結構變化的40.58%，其閾值為11.17；10—20 cm土層土壤真菌群落多元回歸樹第一次分割以堿解氮為節點，將NTS處理與其他處理分解，該次分組能解釋群落結構變化的65.84%，其閾值為6.52。

圖1 不同處理0—10 cm（a）和10—20 cm（b）土層真菌優勢門組成（平均相對豐度＞1%）

NT：免耕；DT：深耕；RT：旋耕；DS：雙季秸稈還田；SS：單季秸稈還田。下同

土壤0—10 cm（a）和10—20 cm（b）土層真菌群落主坐標分析，采用UniFrac距離算法比較不同處理間群落的差異

數據為平均值±標準差，不同小寫字母表示差異顯著（＜0.05）

Data are means±standard deviation; Different small letters mean significant differences (＜0.05）

圖4 不同耕作和秸稈還田方式下土壤0—10 cm和10—20 cm土層真菌功能營養型相對豐度

Fig. 4 The relative abundance of soil fungal functional groups in 0-10 cm and 10-20 cm layers with different tillage and straw returning methods

TOC：有機碳；DOC：可溶性有機碳；TN：全氮；AN：堿解氮；AP：有效磷；AK：速效鉀；Error：相對誤差；CV error：交叉驗證誤差；SE：交叉驗證標準誤

2.5 土壤真菌功能群與環境因子的關系

相關性分析結果顯示，Pathotroph豐度與土壤TOC（=0.600，＜0.001）、DOC（=0.563，＜0.001）、TN（=0.428，＜0.01）、AN（=0.394，＜0.05）和速效鉀（=0.517，＜0.01）均呈顯著正相關關系。除土壤pH外（= 0.570，＜0.01），Symbiotroph豐度與土壤TOC（= -0.535，＜0.01）、TN（= -0.597，＜0.001）、AN（= -0.525，＜0.01）、AP（= -0.698，＜0.001）和AK（= -0.503，＜0.01）均呈顯著負相關關系。

3 討論

3.1 耕作和秸稈還田方式對土壤真菌群落組成的影響

土壤真菌是土壤的重要組成成分，與土壤肥力密切相關[26]。不同的農田管理措施，如種植方式[27-28]，耕作措施[29]及施肥[15,30]等，對土壤真菌群落結構及多樣性的變化均有較大的影響。本研究結果顯示，不同耕作方式與秸稈還田方式下，子囊菌門、擔子菌門、壺菌門和接合菌門，均為農田耕層土壤耕層的優勢真菌。這與前人研究一致，如WANG等[31]江蘇稻麥輪作農田中發現子囊菌門（68.9%）與擔子菌門（19.9%）占統治地位。與滇南茶園[30]，黃土高原春小麥-豌豆輪作[32]和美國小麥-休閑不同耕作方式農田[33]對土壤真菌群落結構的研究一致。然而，肖禮等[27]發現黃土高原3種典型種植類型梯田玉米地、蘋果、蘋果與土豆間作土壤真菌優勢菌群為子囊菌門與接合菌門。另外，也有研究報道波蘭小麥-玉米農田土壤真菌優勢菌群為接合菌門與擔子菌門[34]。

表4 不同耕作與施肥制度土壤真菌功能營養型豐度與土壤理化性質的關系

*表示＜0.05，**表示＜0.01，***表示＜0.001；樣本數=36

*Indicates a significant difference between treatments at＜0.05; **Indicates a significant difference between treatments at＜0.01; ***Indicates a significant difference between treatments at＜0.001; Sample size (=36)

雖然不同耕作方式之間土壤真菌群落物種組成相似，但不同秸稈還田方式土壤主要真菌門相對豐度存在差異。秸稈雙季還田擔子菌相對豐度高于單季還田。這說明擔子菌門的相對豐度受到秸稈還田量的影響。這可能是因為秸稈富含木質素[35]，擔子菌門中的某些種群，如白腐真菌，分解碳氮比高的木質素[36]。雙季秸稈還田為擔子菌提供了一個較好的生長環境，更多地利用降解作物殘留物，從而促進其快速增長。不同處理因作物殘體量的差異及土壤微環境的不同，進而改變了真菌物種組成。

3.2 真菌功能分類對耕作和秸稈還田方式的響應

選擇適宜的耕作措施，增加真菌生物多樣性，可以預防或減少根系疾病的損害，在保持土壤質量和健康方面發揮至關重要的作用。土壤真菌的多樣性決定了生態系統的多樣性和植物生產力[26]。根據營養方式可將真菌分為病理營養型、腐生營養型和共生營養型3種類型。本研究結果顯示，秸稈單季還田共生營養真菌豐度高于秸稈雙季還田（圖4-c），其在作物健康、營養和品質方面具有重要作用[15]。因此，我們認為秸稈單季還田比雙季還田更有益于土壤肥力和作物生產。病理營養型真菌從宿主細胞獲取營養來源，土壤中的病理營養型真菌對植物生長具有一定負面影響[31,37]。本研究中，病理營養型真菌相對豐度為12.4%—37.2%，高于聶三安等[15]黃泥土不同施肥試驗中的相對含量。這說明連續的秸稈還田促進了土壤中的病理營養型真菌生長。不過，深耕一定程度上降低了病理營養型真菌豐度，這說明深耕能夠降低秸稈還田后作物生長過程中潛在的負面影響。不同的耕作方式及秸稈還田方式帶來的營養條件與土壤微環境的改變對病理營養型真菌存在顯著影響，土壤TOC與病理營養型真菌顯著正相關（＜0.001），這說明通過大量秸稈還田，將土壤有機碳富集在耕層土壤，存在一定的負面影響。秸稈還田增加農田病蟲害是我國華北地區麥玉輪作兩季秸稈還田農業生態系統中普遍存在的問題[38]，本研究結果表明深耕能夠減少土壤病原真菌，降低農田病蟲害風險，促進作物生長，提高產量與品質。腐生真菌汲取生長所需營養及氨基酸類物質[39]，但研究結果顯示腐生營養型真菌豐度與土壤TOC和DOC并無相關關系。這可能是因為FUNGuild功能比對是基于已有文獻和數據，僅在一定程度上解析了真菌的功能[31]。本研究中62.00%的土壤真菌功能未被解析出來，復雜的土壤真菌群落功能仍有待深入研究。

3.3 土壤真菌群落組成及功能的相關因子分析

本研究發現土壤TOC（0—10 cm）是影響褐土真菌群落的主要因素，這與LIU等[40]報道東北黑土區與北方旱地區[41]土壤真菌群落的主要影響因子為土壤TOC相一致。對于免耕秸稈雙季還田與旋耕雙季秸稈還田來說，大量秸稈還田提高了表層土壤有機質。有機質作為微生物的能量來源和代謝底物，土壤有機質的含量及其成分顯著影響著微生物的群落組成[42]。秸稈富含木質素，通過秸稈還田土壤中增加的有機碳成分多為復雜且難被分解的木質素等，更易被K-策略的貧營養型真菌群落分解利用[43]。故免耕秸稈雙季還田與旋耕雙季秸稈還田真菌群落區別于其他處理。與此同時，研究發現土壤堿解氮是影響10—20 cm土層真菌群落的主要因素，與裴振等[44]在黃河三角洲鹽堿地土壤中研究結果相似。然而，這與北極圈群島[45]、瑞典北部森林[46]、我國毛竹林[47]和英國農田土壤中[48]，pH是影響真菌群落結構的主要理化性質的研究結果不同。這可能是因為本研究中，不同處理pH變化范圍不大，不足以對真菌群落結構和功能產生影響。以上分析表明，不同秸稈還田及耕作方式對土壤TOC及AN的改變是影響了真菌群落的主要原因。

4 結論

本文通過對土壤真菌Illumina高通量測序結合FUNGuild功能比對，研究華北平原多年定位試驗條件下小麥-玉米周年復種農田土壤真菌群落多樣性和功能對不同耕作與秸稈還田方式的響應。本研究中，Ascomycota和Basidiomycota是華北平原小麥-玉米農田土壤中的優勢真菌，以子囊菌為主。此外，雙季秸稈還田比單季秸稈還田提高了Basidiomycota相對豐度，土壤有機碳和堿解氮是影響土壤真菌群落結構發生變化的主要因素。與免耕和旋耕相比，深耕能夠降低秸稈還田后病理營養型真菌相對豐度，降低作物生長存在潛在的負面影響，有利于保持農田土壤生態系統健康。

[1] FAN M, SHEN J, YUAN L, JIANG R, CHEN X, DAVIES W J, ZHANG F. Improving crop productivity and resource use efficiency to ensure food security and environmental quality in China., 2012, 63(1): 13.

[2] 田慎重, 王瑜, 張玉鳳, 邊文范, 董亮, 羅加法, 郭洪海. 旋耕轉深松和秸稈還田增加農田土壤團聚體碳庫. 農業工程學報, 2017, 24(33): 133-140.

TIAN S Z, WANG Y, ZHANG Y F, BIAN W F, DONG L, LUO J F, GUO H H. Residue returning with subsoiling replacing rotary tillage improving aggregate and associated carbon., 2017, 24(33): 133-140.(in Chinese)

[3] 宋明偉, 李愛宗, 蔡立群, 張仁陟. 耕作方式對土壤有機碳庫的影響. 農業環境科學學報,2008, 27(2): 622-626.

SONG M W, LI A Z, CAI L Q, ZHANG R Z. Effects of different tillage methods on soil organic carbon pool., 2008, 27(2): 622- 626. (in Chinese)

[4] 張建軍, 王勇, 樊廷錄, 郭天文, 趙剛, 黨翼, 王磊, 李尚中. 耕作方式與施肥對隴東旱塬冬小麥-春玉米輪作農田土壤理化性質及產量的影響. 應用生態學報, 2013, 24(4): 1001-1008.

ZHANG J J, WANG Y, FAN T L, GUO T W, ZHAO G, DANG Y, WANG L, LI S Z. Effects of different tillage and fertilization modes on the soil physical and chemical properties and crop yield under winter wheat/spring corn rotation on dryland of east Gansu, Northwest China., 2013, 24(4): 1001-1008.(in Chinese)

[5] 趙亞麗, 薛志偉, 郭海斌, 穆心愿, 李潮海. 耕作方式與秸稈還田對冬小麥-夏玉米耗水特性和水分利用效率的影響. 中國農業科學, 2014, 47(17): 3359-3371.

ZHAO Y L, XUE Z W, GUO H B, MU X Y, LI C H. Effects of tillage and straw returning on water consumption characteristics and water use efficiency in the winter wheat and summer maize rotation system.2014, 47(17): 3359-3371. (in Chinese)

[6] 楊學明, 張曉平, 方華軍, 梁愛珍, 齊曉寧, 王洋. 北美保護性耕作及對中國的意義. 應用生態學報, 2004, 2(15): 335-340.

YANG X M, ZHANG X P, FANG H J, LIANG A Z, QI X N, WANG Y. Conservation tillage systems in North America and their significance for China., 2004, 2(15): 335-340.(in Chinese)

[7] 陳丹梅, 袁玲, 黃建國, 冀建華, 侯紅乾, 劉益仁. 長期施肥對南方典型水稻土養分含量及真菌群落的影響. 作物學報, 2017, 43(2): 286-295.

CHEN D M, YUAN L, HUANG J G, JI J H, HOU H Q, LIU Y R. Influence of long-term fertilizations on nutrients and fungal communities in typical paddy soil of South China., 2017, 43(2): 286-295. (in Chinese)

[8] 孫冰潔, 賈淑霞, 張曉平, 梁愛珍, 陳學文, 張士秀, 劉四義, 陳升龍. 耕作方式對黑土表層土壤微生物生物量碳的影響. 應用生態學報, 2015, 26(1): 101-107.

SUN B J, JIA S X, ZHANG X P, LIANG A Z, CHEN X W, ZHANG S X, LIU S Y, CHEN S L. Impact of tillage practices on microbial biomass carbon in top layer of black soils., 2015, 26(1): 101-107.(in Chinese)

[9] AZOOZ R H, ARSHAD M A, FRANZLUEBBERS A J. Pore size distribution and hydraulic conductivity affected by tillage in northwestern Canada．, 1996, 60: 1197-1201.

[10] 郭梨錦. 免耕與稻稈還田對稻麥種植系統土壤有機碳庫與微生物多樣性的影響[D]. 武漢: 華中農業大學, 2018.

GUO L J. effects of no-tillage and straw return on soil organic carbon pool and microbial diversity in rice-wheat systems[D]. Wuhan: Huazhong Agricultural University, 2018. (in Chinese)

[11] WANG Z T, LIU L, CHEN Q, WEN X X, LIAO Y C. Conservation tillage increases soil bacterial diversity in the dryland of northern China., 2016, 36(2): 28-36.

[12] SUN B J, JIA S X, ZHANG S X, MCLAUGHLIN N B, ZHANG X P, LIANG A Z, CHEN X W, WEI S C, LIU S Y. Tillage, seasonal and depths effects on soil microbial properties in black soil of Northeast China., 2016, 155: 421-428.

[13] 郭梨錦, 曹湊貴, 張枝盛, 劉天奇, 李成芳. 耕作方式和秸稈還田對稻田表層土壤微生物群落的短期影響. 農業環境科學學報, 2013, 32(8): 1577-1584.

GUO L J, CAO C G, ZHANG Z S, LIU T Q, LI C F.Short-term effects of tillage practices and wheat-straw returned to rice fields on topsoil microbial community structure and microbial diversity in central China., 2013, 32(8): 1577-1584. (in Chinese)

[14] 姚曉東, 王娓, 曾輝. 磷脂脂肪酸法在土壤微生物群落分析中的應用. 微生物學通報, 2016, 43(9): 2086-2095.

YAO X D, WANG W, ZENG H. Application of phospholipid fatty acid method in analyzing soil microbial community composition., 2016, 43(9): 2086-2095. (in Chinese)

[15] 聶三安, 王祎, 雷秀美, 趙麗霞, 林瑞余, 王飛, 邢世和. 黃泥田土壤真菌群落結構和功能類群組成對施肥的響應. 應用生態學報, 2018, 29(8): 2721-2729.

NIE S A, WANG W, LEI X M, ZHAO L X, LIN R Y, WANG F, XING S H. response of fungal community structure and functional group to fertilization in yellow clayey soil., 2018, 29(8): 2721-2729. (in Chinese)

[16] 鮑士旦. 土壤農化分析. 北京: 中國農業出版社, 2000: 34-35, 44-49, 56-58, 81, 106-107.

BAO S D.. Beijing: China Agriculture Press, 2000: 34-35, 44-49, 56-58, 81, 106-107. (in Chinese)

[17] ROUSK J, B??TH E, BROOKES P C, LAUBER C L, LOZUPONE C, CAPORASO J G, KNIGHT R, FIERER N. Soil bacterial and fungal communities across a pH gradient in an arable soil., 2010, 4(10): 1340-1351.

[18] EDGAR R C. UPARSE: highly accurate OTU sequences from microbial amplicon reads., 2013, 10(10): 996-998.

[19] CAPORASO J G, KUCZYNSKI J, STOMBAUGH J, BITTINGER K, BUSHMAN F D, COSTELLO E K, FIERER N, PEA A G, GOODRICH J K, GORDON J I, HUTTLEY G A, KELLEY S T, KNIGHTS D, KOENIG J E, LEY R E, LOZUPONE C A, MCDONALD D, MUEGGE B D, PIRRUNG M, REEDER J, SEVINSKY J R, TURNBAUGH P J, WALTERS W A, WIDMANN J, YATSUNENKO T, ZANEVELD J, KNIGHT R. QIIME allows analysis of high-throughput community sequencing data., 2010, 7: 335-336.

[20] WANG Q, GARRITY G M, TIEDJE J M, COLE J R. Naive bayesian classifier for rapid assignment of rrna sequences into the new bacterial taxonomy., 2007, 73(16): 5261-5267.

[21] QUAST C, PRUESSE E, YILMAZ P, GERKEN J, SCHWEER T, YARZA P, PEPLIES J, GL?CKNER F O. The silva ribosomal RNA gene database project: improved data processing and web-based tools., 2013, 41(1): 590-596.

[22] SCHLOSS P D, GEVERS D, WESTCOTT S L. Reducing the effects of PCR amplification and sequencing artifacts on 16S rRNA-based studies., 2011, 6(12): 1-14.

[23] CATHERINE L, ROB K. UniFrac: a new phylogenetic method for comparing microbial communities., 2005, 71: 8228-8235.

[24] DE'ATH G. Multivariate Regression trees: a new technique for modeling species-environment relationships., 2002, 83(4): 1105-1117.

[25] NGUYEN N H, SONG Z, BATES S T, BRANCO S, TEDERSOO L, MENKE J, SCHILLING, JONATHAN S, KENNEDY P G. Funguild: an open annotation tool for parsing fungal community datasets by ecological guild., 2015, 20: 241-248.

[26] FR?C M, HANNULA S E, BE?KA M, J?DRYCZKA M. Fungal biodiversity and their role in soil health., 2018, 9: 1-9.

[27] 肖禮, 黃懿梅, 趙俊峰, 周俊英, 郭澤慧, 劉洋. 土壤真菌組成對黃土高原梯田種植類型的響應. 中國環境科學, 2017, 37(8): 3151-3158.

XIAO L, HUANG Y M, ZHAO J F, ZHOU J Y, GUO Z H, LIU Y. High-throughput sequencing revealed soil fungal communities under three terrace agrotypes on the loess plateau., 2017, 37(8): 3151-3158. (in Chinese)

[28] 李銳, 劉瑜, 褚貴新. 不同種植方式對綠洲農田土壤酶活性與微生物多樣性的影響. 應用生態學報, 2015, 26(2): 490-496.

LI R, LIU Y, CHU G X. Effects of different cropping patterns on soil enzyme activities and soil microbial community diversity in oasis farmland.,2015, 26(2): 490-496. (in Chinese)

[29] 何玉梅, 張仁陟, 張麗華, 解開治. 不同耕作措施對土壤真菌群落結構與生態特征的影響. 生態學報, 2007, 27(1): 113-119.

HE Y M, ZHANG R Z, ZHANG L H, XIE K Z. Effects of different tillage practices on fungi community structure and ecologic characteristics in loess soils., 2007, 27(1): 113-119. (in Chinese)

[30] 季凌飛, 倪康, 馬立鋒, 陳兆杰, 趙遠艷, 阮建云, 郭世偉. 不同施肥方式對酸性茶園土壤真菌群落的影響. 生態學報, 2018, 38(22): 1-8.

JI L F, NI K, MA L F, CHEN Z J, ZHAO Y Y, RUAN J Y, GUO S W. Effect of different fertilizer regimes on the fungal community of acidic tea-garden soil.,2018, 38(22): 1-8. (in Chinese)

[31] WANG J, RHODES G, HUANG Q, SHEN Q R. Plant growth stages and fertilization regimes drive soil fungal community compositions in a wheat-rice rotation system., 2018, 54(6): 731-742.

[32] ESSEL E, LI L, DENG C, XIE J, ZHANG R, LUO Z, CAI L. Evaluation of bacterial and fungal diversity in a long-term spring wheat-field pea rotation field under different tillage practices., 2018, 98: 1-19.

[33] SCHLATTER D C, SCHILLINGER W F, BARY A I, SHARRATT B, PAULITZ T C. Biosolids and conservation tillage: Impacts on soil fungal communities in dryland wheat-fallow cropping systems., 2017, 115: 556-567.

[34] ANNA G, JAROS?AW G. Fungal genetics and functional diversity of microbial communities in the soil under long-term monoculture of maize using different cultivation techniques., 2018, 9: 1-9.

[35] 張經廷, 張麗華, 呂麗華, 董志強, 姚艷榮, 金欣欣, 姚海坡, 賈秀領. 還田作物秸稈腐解及其養分釋放特征概述. 核農學報, 2018, 32(11): 2274-2280.

ZHANG J T, ZHANG L H, LYU L H, DONG Z Q, YAO Y R, JIN X X, YAO H P, JIA X L. Overview of the characteristics of crop straw decomposition and nutrients release of returned field crops., 2018, 32(11): 2274-2280. (in Chinese)

[36] LAUBER C L, STRICKLAND M S, BRADFORD M A, FIERER N. The influence of soil properties on the structure of bacterial and fungal communities across land-use types., 2008, 40(9): 2407-2415.

[37] NTHONY M A, FREY S D, STINSON K A. Fungal community homogenization, shift in dominant trophic guild, and appearance of novel taxa with biotic invasion., 2017, 8(9): 1-17.

[38] 董印麗, 李振峰, 王若倫,卜學平, 付建敏, 董秀秀. 華北地區小麥、玉米兩季秸稈還田存在問題及對策研究. 中國土壤與肥料, 2018(1): 159-163.

DONG Y L, LI Z F, WANG R L, BU X P, FU J M，DONG X X. Study on the problems and countermeasures of returning wheat and corn stalks into the soil in north China., 2018(1): 159-163. (in Chinese)

[39] 李秀璋. 醉馬草內生真菌與宿主種帶真菌、根際微生物的互作及其進化研究[D]. 蘭州: 蘭州大學, 2017.

LI X Z. Study on the evolution and interactions ofwith host seed-borne fungi and rhizospheric microorganism[D]. Lanzhou: Lanzhou University, 2017. (in Chinese)

[40] LIU J, SUI Y, YU Z, SHI Y, CHU, H, JIN J, LIU X, WANG G. Soil carbon content drives the biogeographical distribution of fungal communities in the black soil zone of northeast China., 2015, 83: 29-39.

[41] WANG Z, CHEN Q, LIU L, WEN X, LIAO Y. Responses of soil fungi to 5-year conservation tillage treatments in the drylands of northern China., 2016, 101: 132-140.

[42] 王慧穎, 徐明崗, 周寶庫, 段英華. 黑土細菌及真菌群落對長期施肥響應的差異及其驅動因素. 中國農業科學, 2018, 51(5): 914-925.

WANG H Y, XU M G, ZHOU B K, MA X, DUAN Y H. Response and driving factors of bacterial and fungal community to long-term fertilization in black soil., 2018, 51(5): 914-925. (in Chinese)

[43] STRICKLAND M S, ROUSK J. Considering fungal:bacterial dominance in soils–Methods, controls, and ecosystem implications., 2010, 42(9): 1385-1395.

[44] 裴振, 孔強, 郭篤發. 鹽生植被演替對土壤微生物碳源代謝活性的影響. 中國環境科學, 2017, 37(1): 373-380.

PEI Z, KONG Q, GUO D F. Effect of succession of halophytic vegetation on soil microbial carbon metabolic activity., 2017, 37(1): 373-380. (in Chinese)

[45] ZHANG T, WANG N, LIU H, ZHANG Y Q, YU L Y. Soil pH is a key determinant of soil fungal community composition in the Ny-?lesund region, Svalbard (High Arctic)., 2016, 7: 1-10.

[46] H?GBERG M N, YARWOOD S A, MYROLD D D. Fungal but not bacterial soil communities recover after termination of decadal nitrogen additions to boreal forest., 2014, 72: 35-43.

[47] 李永春, 梁雪, 李永夫, 王祈, 陳俊輝, 徐秋芳. 毛竹入侵闊葉林對土壤真菌群落的影響. 應用生態學報, 2016, 27(2): 585-592.

LI Y C, LIANG X, LI Y F, WANG Q, CHEN J H, XU Q F. Effects ofinvasion of native broadleaf forest on soil fungal community., 2016, 27(2): 585-592. (in Chinese)

[48] ROUSK J, BROOKES P C, BAATH E. Contrasting soil pH effects on fungal and bacterial growth suggest functional redundancy in carbon mineralization., 2009, 75(6): 1589-1596.

Response of Fungal Community and Function to Different Tillage and Straw Returning Methods

DAI HongCui1,3, ZHANG Hui2,3, XUE YanFang2,3, GAO YingBo2,3, QIAN Xin2,3，ZHAO HaiJun3, CHENG Hao4, LI ZongXin2,3, LIU KaiChang1,3

(1Crop Research Institute, Shandong Academy of Agricultural Sciences, Jinan 250100;2Maize Research Institute, Shandong Academy of Agricultural Sciences, Jinan 250100;3National Engineering Laboratory for Wheat and Maize, Jinan 250100;4College of Agronomy, Qingdao Agricultural University, Qingdao 266109, Shandong)

【】This study was conducted to explore the change of fungal community structure and function response to different tillage and straw returning methods in wheat-maize rotation system in the North China Plain. It aimed to clarify the biological mechanism of soil fertility improvement, which provided a theoretical support for sustainable development for agricultural production.【】A six-year field study with split plot design was conducted to investigate the effects of different soil tillage methods (no tillage, CT; deep tillage, DT; rotation tillage, ST) and straw returning methods (wheat and maize straws were returned to the field, DS; only wheat straw was returned to the field, SS) on changes of fungal community structure and function in soils from wheat-maize rotation system in the North China Plain. In combination with soil properties, multiple regression trees and correlation analysis was carried out to investigate driving factors of fungal community structure and function in soil．【】The results showed that, compared with NT, soil organic carbon content under DS and SS were reduced by 35.04% and 44.30% in 0-10 cm layer, respectively. The available nitrogen of NT under SS treatment was significantly lower than that under other treatments in 10-20 cm layer. Ascomycetes (68.98%), Basidiomycetes (16.96%) and Chytridiomycetes (1.62%) were the dominant fungus in 0-10 cm layer, while Ascomycetes (68.44%), Basidiomycetes (15.52%), Chytridiomycetes (1.51%) and Coccidiomycetes (1.23%) were the dominant fungus in 10-20 cm layer. The different tillage and straw returning methods changed soil fungal community structure. Specifically, the relative abundance ofin DS increased 50.07% and 29.08% respectively in 0-10 cm and 10-20 cm layers than that of SS. The multiple regression trees showed that soil fungal communities were divided into soil organic carbon nodes with a threshold of 11.17 mg·kg-1in 0-10 cm layer, additionally, the soil fungi community were divided into available nitrogen nodes with a threshold of 6.52 mg·kg-1in 10-20 cm layer. In this study, Pathotroph was mainly function type of soil fungi in 0-10 cm (26.84%) and 10-20 cm (23.91%) layers in wheat-maize rotation in the North China Plain. Compared with NT, Pathotroph relative abundance under DT and RT treatments were reduced by 25.16% and 16.45%, respectively. The results of correlation analysis showed that Pathotroph relative abundance were positively correlated with soil total organic carbon, dissolved organic carbon, total nitrogen, available nitrogen and available potassium.【】In general, our results indicated that different tillage and straw returning methods changed soil fungal community structure and relative abundance of functional groups. The content of soil organic carbon and available nitrogen were the driving factors shaped the fungal community structure. Besides, DT could reduce Pathotroph relative abundance, which was conducive to maintaining the soil ecosystem health.

tillage practice; straw returning; soil fungal community; high-throughput sequencing; multiple regression trees

10.3864/j.issn.0578-1752.2019.13.008

2018-12-24；

2019-02-26

國家重點研發計劃（2018YFD0300606）、山東省現代農業產業技術體系玉米創新團隊項目（SDAIT-02-07，SDAIT-02-11）、山東省農業科學院農業科技創新工程（CXGC2016A05）

代紅翠，E-mail：daihongcui2013@163.com。

劉開昌，Tel：0531-66659091；E-mail：liukc1971@126.com。通信作者李宗新，Tel：0531-66659402；E-mail：sdaucliff@sina.com

（責任編輯 李云霞）