長期有機無機肥配施對東北黑土真菌群落結(jié)構(gòu)的影響

丁建莉，姜昕*，馬鳴超，關(guān)大偉，趙百鎖，魏丹，曹鳳明，李力，李俊*

（1中國農(nóng)業(yè)科學(xué)院農(nóng)業(yè)資源與農(nóng)業(yè)區(qū)劃研究所，北京100081；2農(nóng)業(yè)部微生物產(chǎn)品質(zhì)量安全風(fēng)險評估實驗室，北京100081；3黑龍江省農(nóng)業(yè)科學(xué)院土壤肥料與環(huán)境資源研究所，哈爾濱150086）

長期有機無機肥配施對東北黑土真菌群落結(jié)構(gòu)的影響

丁建莉1,2，姜昕1,2*，馬鳴超1,2，關(guān)大偉1,2，趙百鎖2，魏丹3，曹鳳明1,2，李力1,2，李俊1,2*

（1中國農(nóng)業(yè)科學(xué)院農(nóng)業(yè)資源與農(nóng)業(yè)區(qū)劃研究所，北京100081；2農(nóng)業(yè)部微生物產(chǎn)品質(zhì)量安全風(fēng)險評估實驗室，北京100081；3黑龍江省農(nóng)業(yè)科學(xué)院土壤肥料與環(huán)境資源研究所，哈爾濱150086）

【目的】分析長期有機無機肥配施對土壤真菌豐度、多樣性及其群落特征的影響，探討東北黑土真菌群落變化與施肥的相關(guān)性，為進一步調(diào)節(jié)土壤微生物結(jié)構(gòu)，改善其生態(tài)功能提供參考依據(jù)。【方法】依托黑龍江省農(nóng)業(yè)科學(xué)院36年長期定位試驗站，選取4個不同施肥處理：不施肥處理(CK)、有機肥處理(M)、氮磷鉀無機肥處理(NPK)和有機肥配施無機肥處理(MNPK)的耕作層土壤為研究對象，以真菌ITS基因為分子標(biāo)靶借助qPCR技術(shù)和Illumina Miseq高通量測序平臺，研究不同施肥處理對黑土中真菌群落豐度、多樣性和組成的影響，并與土壤化學(xué)性質(zhì)進行偶聯(lián)分析，揭示群落與施肥的相關(guān)性?！窘Y(jié)果】長期施用無機肥顯著降低土壤pH，而有機無機配施可以有效緩解土壤酸化。NPK處理的ITS基因豐度顯著高于MNPK；MNPK處理的細(xì)菌/真菌比值(26.91×104)顯著大于NPK，各處理比值由高到低為MNPK>M>CK>NPK。細(xì)菌/真菌比值與土壤pH正相關(guān)；MNPK處理的真菌α多樣性指數(shù)值略大于NPK。Ascomycota和Zygomycota為土壤中主要真菌門，不同施肥處理之間真菌組成的相對豐度存在顯著差異，對照處理Ascomycota的相對豐度為45.35%，MNPK和NPK處理分別為50.93%和56.16%。有機肥有利于降低病原真菌相對豐度，具有高度侵染性的Cochliobolus在MNPK (0.41%)和M(0.39%)中的相對豐度顯著小于CK(3.25%)和NPK(2.08%)。CCA分析表明，土壤理化性質(zhì)共解釋土壤真菌群落結(jié)構(gòu)變化的73.3%，有效磷(貢獻(xiàn)量為32.4%，P=0.002)、銨態(tài)氮(貢獻(xiàn)量為14.8%，P=0.01)和硝態(tài)氮(貢獻(xiàn)量為16.2%，P=0.048)是3個重要的影響因子?！窘Y(jié)論】不同施肥條件下土壤真菌豐度、多樣性，以及菌群組成特征不同。與無機肥相比，有機肥無機肥配施能夠有效改善真菌群落結(jié)構(gòu)，降低真菌的豐度，增加真菌多樣性，并提高土壤pH，減緩?fù)寥浪峄?。土壤有效磷、銨態(tài)氮和硝態(tài)氮含量是影響黑土土壤真菌群落結(jié)構(gòu)變化重要因素。

施肥；真菌群落；高通量測序；土壤有效磷；土壤銨態(tài)氮；土壤硝態(tài)氮

黑土有機質(zhì)含量高、土質(zhì)肥沃，是最寶貴的土壤資源。我國東北黑土作為世界四大黑土區(qū)之一，面積達(dá)103萬km2、占全國陸地總面積的10.7%，其中耕地面積3000萬hm2、占全國耕地總面積的22.2%[1]，因此東北黑土在保障國家糧食安全、區(qū)域生態(tài)環(huán)境安全和農(nóng)業(yè)可持續(xù)發(fā)展中具有不可替代的地位。然而由于長期不合理施肥導(dǎo)致土壤質(zhì)量退化，如長期施用無機肥使得黑土土壤侵蝕嚴(yán)重、土層變薄、有機質(zhì)含量下降、土壤酸化、土壤微生物活性降低等一系列問題[2–3]。而施用有機肥是提高作物產(chǎn)量、增加土壤有機質(zhì)和改善土壤生物學(xué)性狀的有效措施[4–5]。但近五十年來我國有機肥使用量持續(xù)下降，調(diào)查數(shù)據(jù)表明，在1960、1980、2000、2010年我國有機肥施用量占總施肥量的80%、60%、30%、10%[6]。面對無機肥施用量增加，有機肥施用量逐年減少的現(xiàn)狀，我國提出了有機無機肥配施的施肥指導(dǎo)策略。研究土壤微生物生物量、群落組成等對施肥、耕作等措施的快速響應(yīng)變化，可以指示土壤質(zhì)量變化[7]。

關(guān)于施肥對土壤微生物影響的研究報道很多，有機肥與無機肥配施會提高土壤微生物生物量碳氮、土壤酶活[8]。真菌作為土壤微生物中的重要組成部分，在土壤生態(tài)中具有不可忽視的地位，與土傳病害、植物互作、分解有機物質(zhì)等密不可分[9–10]，因此探究施肥對土壤真菌的影響具有重要意義。研究表明施入氮肥會降低真菌的生物量[11]，降低多樣性和改變真菌的組成結(jié)構(gòu)[12–13]；相對細(xì)菌來說，真菌有更強的吸收氮素和磷素的能力[14]。無機肥減少菌根真菌的生長，而糞肥增加菌根真菌的生長[15]。施有機肥土壤理化性質(zhì)與常規(guī)農(nóng)田沒有顯著差異，但生物量高于常規(guī)農(nóng)田，且菌根真菌量增加[16]。在溫室試驗中發(fā)現(xiàn)，生物有機肥可以降低土壤真菌的數(shù)量，尤其是病原真菌的數(shù)量，降低作物病害發(fā)病率，且改變真菌的組成結(jié)構(gòu)[17]。高通量測序技術(shù)可進行環(huán)境中所有微生物群落的研究，獲得數(shù)據(jù)量大，能更真實地揭示微生物群落的復(fù)雜性和多樣性，極大的促進環(huán)境中不可培養(yǎng)微生物以及痕量菌的更深入研究[18–19]。本課題組的Zhou等[20]借助高通量測序技術(shù)研究不同施氮量和不同量氮、磷肥配施對東北黑土真菌群落結(jié)構(gòu)變化的影響，發(fā)現(xiàn)氮肥和氮、磷肥配施可增加真菌數(shù)量，降低真菌多樣性，改變真菌群落結(jié)構(gòu)的組成。本研究首次應(yīng)用高通量測序技術(shù)研究長期有機肥無機肥配施條件下東北黑土真菌群落結(jié)構(gòu)是如何演變的，及其與土壤理化性質(zhì)的偶聯(lián)關(guān)系。

真菌核糖rDNA內(nèi)部轉(zhuǎn)錄間隔區(qū)(internal transcribed spacer，ITS)序列分析技術(shù)已被應(yīng)用到物種屬內(nèi)、近緣屬間乃至科內(nèi)研究系統(tǒng)進化關(guān)系[21]。由于ITS的序列分析能實質(zhì)性地反映出屬間、種間的堿基對差異，且序列片段較小、易于分析，目前已被廣泛應(yīng)用于真菌屬內(nèi)不同種間或近似屬間的系統(tǒng)發(fā)育研究中[22]。本研究選取4個施肥處理，探討長期有機無機配施和單施無機肥條件下，真菌微生物豐度、多樣性和組成的演變規(guī)律，及其與土壤理化性質(zhì)的偶聯(lián)關(guān)系。借助qPCR技術(shù)和Illumina Miseq高通量測序平臺，以真菌ITS基因為分子標(biāo)靶探討以上問題。研究結(jié)果將為東北黑土建立合理施肥制度，維持土壤肥力，調(diào)節(jié)土壤微生物結(jié)構(gòu)，改善其生態(tài)功能提供參考依據(jù)。

1 材料與方法

1.1 試驗地概況

試驗在黑龍江省農(nóng)業(yè)科學(xué)院長期定位試驗站進行，試驗站位于黑龍江省哈爾濱市(45°40'N，126°35'E)，地勢平坦，屬于松花江二級階地，海拔151m，無霜期135d，年均氣溫3.5℃，年降水量533mm。該試驗站1979年成立，自次年起，按小麥–大豆–玉米順序輪作，2015年為玉米季。本文選取4個施肥處理：不施肥處理(CK)；有機肥處理(M)；施氮肥、磷肥和鉀肥的無機肥處理(NPK)；有機肥配施無機肥處理(MNPK)。每個處理3個重復(fù)。氮、磷、鉀肥分別為尿素、重過磷酸鈣、硫酸鉀。氮、磷、鉀肥小麥季和玉米季施用量為N150kg/hm2、P2O575 kg/hm2、K2O75kg/hm2；大豆季施用量為N75 kg/hm2、P2O5150kg/hm2、K2O75kg/hm2。有機肥為馬糞，施用量為18600kg/hm2?；A(chǔ)土壤理化性質(zhì)(1980年)：pH7.22、有機質(zhì)含量26.7g/kg、全氮1.47g/kg、全磷1.07g/kg、堿解氮151.1mg/kg、速效鉀200mg/kg、有效磷51.0mg/kg。

1.2 樣品采集

玉米收獲季(2015年9月)取5—25cm的耕層土壤，在每個小區(qū)隨機抽取10個點取樣，去除土壤中雜草石塊等雜質(zhì)，混合均勻作為該處理的一個平行樣品，每個處理3個重復(fù)。土壤樣本保存于封口塑料袋中，一部分于溫室風(fēng)干、研磨后過0.2mm篩，用于測定土壤理化性質(zhì)(用鮮土測定土壤銨態(tài)氮和硝態(tài)氮含量)；另一部分樣品保存于–80℃冰箱，用于微生物群落分析。

1.3 試驗方法

1.3.1 土壤理化性質(zhì)測定以下土壤理化性質(zhì)參照魯如坤[23]方法，采用酸度計法(土水比為1∶1)測定土壤pH；采用半微量凱氏定氮法測定土壤全氮；Olsen法測定土壤有效磷；采用乙酸銨浸提—原子吸收分光光度法測定土壤速效鉀；通過重鉻酸鉀—容量法—外加熱法測定土壤有機質(zhì)。全自動流動分析儀(AA3)測定硝態(tài)氮和銨態(tài)氮的含量[24]。

1.3.2 土壤DNA的提取和ITS基因的高通量測序使用美國MOBIO公司的Power Max Soil DNA Isolation Kit試劑盒提取土壤總DNA，準(zhǔn)確稱取0.25g土壤樣品，依據(jù)試劑盒說明書方法提取土壤總DNA。為了使土壤總DNA具有重復(fù)性和代表性，每個處理的每個重復(fù)樣品，提取6次DNA，把所提取的6次DNA混合均勻作為該重復(fù)土壤總DNA進行后續(xù)實驗[25]。用核酸定量儀(NanoDrop ND-1000)檢測土壤總DNA濃度和純度。符合要求的DNA送至北京奧維森基因科技有限公司應(yīng)用Illumina MiSeq平臺進行測序。

測序引物為ITS1F/ITS2(ITS1F:5′-CTTGGTCATTTAGAGGAAGTAA-3′,ITS2:5′-GCTGCGTTCTTCATCGATGC-3′)[26]。Illumina Miseq平臺測序的主要步驟：首先合成帶有barcode的特異引物，并使用Qubit熒光定量系統(tǒng)對PCR產(chǎn)物進行定量檢測。再依據(jù)每個樣本的測序量要求以及定量結(jié)果進行相應(yīng)比例的混合。然后將含有標(biāo)簽的接頭與DNA片段鏈接。選擇性地富集兩端連有接頭的DNA片段，擴增DNA文庫，使用測序平臺對真菌ITS1區(qū)基因進行測序。原始序列上傳至NCBI的SRA數(shù)據(jù)庫編號為SRP085784。

1.3.3 定量PCR分析采用SYBR Green定量PCR分析ITS基因，反應(yīng)在ABI7500Real-time PCR (ABI，USA)儀器上進行。反應(yīng)體系為FastFire qPCR PreMix(TIANGEN，China)10μL，ROX Reference Dye0.4μL，1μL DNA模板，10nmol/L引物，補加ddH2O至20μL。真菌定量PCR引物為ITS4/ITS5

(ITS4:5′-TCCTCCGCTTATTGATATGC-3′,ITS5:5′-GGAAGTAAAAGTCGTAACAAGG-3′)[27]。以含有目標(biāo)基因的重組pGEMR-T載體為標(biāo)準(zhǔn)質(zhì)粒，制備質(zhì)粒和后續(xù)方法參照文獻(xiàn)[28]。細(xì)菌16S rRNA基因豐度檢測方法同真菌，其引物為515F/806R[29]。

1.3.4 高通量測序數(shù)據(jù)分析通過Illumina MiSeq平臺進行Paired-end測序，使用Trimmomatic軟件修剪數(shù)據(jù)末端堿基，以Q20≥90的標(biāo)準(zhǔn)用Readfq (vertion6.0)軟件將低質(zhì)量的成對reads過濾掉；采用FLASH(version1.2.10)軟件把成對的reads通過其overlap拼接到一條序列，拼接的錯誤匹配率為0.1。Mothur(version1.31.2)軟件去除長度小于200 bp的序列和maxhomop大于10的序列，并去除嵌合體。QIIME(v1.8.0)軟件將拼接過濾后的序列聚類為用于物種分類的OTU(Operational Taxonomic Units)，OUT相似性設(shè)置為97%。為了得到每個OTU對應(yīng)的物種分類信息，對OTU代表序列進行分類學(xué)分析，并在各個水平統(tǒng)計每個樣品的群落組成，對比UNITE數(shù)據(jù)庫。利用軟件Mothur進行α多樣分析。采用SPSS19.1軟件進行單因素ANOVA分析和Pearson相關(guān)性分析，并用Turkey顯著差異法分析處理間的差異顯著性。使用PAST軟件，基于Bray-Curtis相似距離，對黑土中真菌群落結(jié)構(gòu)進行聚類分析[30]。采用CANOCO5.0將土壤理化性質(zhì)和真菌群落組成進行典范對應(yīng)分析(CCA)。

2 結(jié)果與分析

2.1 不同施肥處理土壤理化性質(zhì)

由表1可知，與不施肥相比，單施無機肥降低了土壤pH，施有機肥提高了pH，NPK處理的土壤pH值比CK降低0.53，MNPK的pH要高出NPK處理0.40。有機無機配施和單施無機肥處理極顯著增加土壤中全氮、硝態(tài)氮、有效磷和有機質(zhì)等的含量(P<0.05)，MNPK處理的有效磷含量較CK處理提高60.14mg/kg，NPK處理也較CK處理提高57.81 mg/kg；而施肥對不同處理的銨態(tài)氮、速效鉀含量影響不顯著。總體上來說，施肥改變了土壤理化性質(zhì)。

表1 不同施肥處理對黑土土壤理化性質(zhì)的影響分析Table 1 Variations of soil physicochemical properties under different fertilizer treatments

2.2 不同施肥處理的黑土ITS基因豐度和細(xì)菌/真菌比值分析

通過qPCR檢測黑土中ITS基因和16S rRNA基因豐度，分析黑土土壤中真菌豐度，以及細(xì)菌/真菌比值(16S rRNA基因豐度值與ITS基因豐度值做商)，結(jié)果見圖1。各處理ITS基因拷貝數(shù)為1.66×105～3.04×105/g,soil)，NPK處理的最高(3.04×105/g,soil)，比MNPK處理(2.02×105/g,soil)的高出1.02×105/g, soil，不同處理之間存在顯著差異。通過分析細(xì)菌/真菌比值可知，MNPK處理(26.91×104)和M處理(23.61×104)均高于CK(19.58×104)，而NPK處理(9.90×104)小于CK。另外，通過土壤理化性質(zhì)和ITS基因豐度，以及細(xì)菌/真菌比值的相關(guān)性分析(表2)可知，ITS基因拷貝數(shù)同銨態(tài)氮和速效鉀含量顯著正相關(guān)；細(xì)菌/真菌比值同土壤pH值顯著正相關(guān)，與銨態(tài)氮含量顯著負(fù)相關(guān)。

2.3 真菌 α 多樣性分析

不同施肥處理之間的真菌α多樣性分析結(jié)果(表3)表明，只有Observed Species分析，CK處理的最高(735)，與其他3個處理差異達(dá)到顯著水平，3個施肥處理間雖有差異，但未達(dá)顯著水平。Chao1、PD whole tree和Shannon指數(shù)，所有處理間均沒有顯著性差異。

2.4 不同施肥處理的黑土土壤真菌組成分析

由圖2中可以看到，長期施肥的黑土土壤中檢測到6個真菌門，相對豐度大于1的有5個，分別為：Ascomycota(相對豐度35.47%～56.22%)，Zygomycota(18.37%～25.66%)，Basidiomycota (8.26%～12.14%)，Glomeromycota(1.19%～2.63%)，Chytridiomycota(1.96%～9.87%)。相對豐度大于1%的菌綱有10個和相對豐度大于0.1%的菌屬40個(圖2)。

不同施肥處理土壤5個真菌門Ascomycota、Zygomycota、Basidiomycota、Glomeromycota和Chytridiomycota的相對豐度存在差異(圖3)。CK處理Ascomycota的相對豐度為45.35%，M處理的為35.47%，而MNPK處理和NPK處理的分別提高到50.93%和56.16%。Zygomycota的相對豐度在CK處理為20.91%，在M處理比CK增加了4.75%，而在MNPK處理和NPK處理分別減少為18.37%和19.95%?？梢婇L期施肥引起了真菌優(yōu)勢菌門相對豐度的差異。

圖1 施肥處理ITS基因豐度和細(xì)菌/真菌比值的qPCR結(jié)果Fig. 1 The abundance of fungi as indicated by the number of ITS copies and bacteria-to-fungi ratio in different treatments analyzed by qPCR

表2 ITS基因拷貝數(shù)、細(xì)菌/真菌比值與土壤理化性質(zhì)Pearson相關(guān)性分析Table 2 Pearson correlation coefficients between ITS gene sequence copy numbers, bacteria-to-fungi ratio and soil physiochemical characteristics

表3 不同施肥處理黑土土壤中真菌 α 多樣性指數(shù)Table 3 α-diversity index of fungi under different fertilizer treatments

另外，不同施肥處理土壤真菌屬水平組成分析(表4)，4個處理之間優(yōu)勢菌屬Fusarium的相對豐度沒有顯著性差異；而Cochliobolus在CK處理和NPK處理中分別為3.25%和2.08%，在施入有機肥的MNPK處理和M處理中的相對豐度分別降低到0.41%和0.39%。Exophiala和Alternaria也有類似變化趨勢，它們的相對豐度在NPK處理(分別為5.38%和2.49%)均顯著高于MNPK處理(分別為3.15%和1.65%)。

2.5 不同施肥處理的聚類分析

聚類分析根據(jù)每個樣品OTU的組成情況，更直觀地展示各個樣品之間的關(guān)系，組成越相似的樣品，聚類關(guān)系越近。通過此分析，具有相似β多樣性的供試樣品聚類在一起，本研究各處理的結(jié)果(圖4)顯示CK處理的幾個重復(fù)先行相聚，再同M和MNPK聚類，最后與NPK聚類。4個處理之間彼此分開，說明處理之間的群落結(jié)構(gòu)存在差異；MNPK處理與NPK處理相比較，MNPK處理與CK處理的親緣關(guān)系更近。

圖2 長期施肥黑土土壤真菌組成Fig. 2 Fungal community composition in long-term fertilized soils

圖3 不同施肥處理的真菌門水平相對豐度Fig. 3 Relative abundances of fungal common phyla under different long-term fertilizer treatments

表4 不同施肥處理黑土土壤中真菌屬相對豐度 (%)Table 4 Relative abundances of fungal common genera under different treatments

圖4 不同施肥處理黑土土壤中真菌組成聚類分析 (OTU水平)Fig. 4 Results of clustering analysis based on fungal community composition in different fertilized soils (at OUT-level)

2.6 CCA分析

CCA分析主要反映樣品、菌群與環(huán)境因子之間關(guān)系。通過CCA分析長期不同施肥處理土壤中真菌群落與土壤理化性質(zhì)之間的相互關(guān)系(圖5)，本文選取的土壤理化性質(zhì)包括：土壤pH、有機質(zhì)、全氮、硝態(tài)氮、銨態(tài)氮、有效磷、速效鉀等。土壤理化性質(zhì)共解釋土壤真菌群落結(jié)構(gòu)變化的73.3%，前兩軸累計解釋量為56.53%，第一軸和第二軸的解釋量分別為36.54%、19.99%。分析得知，對真菌群落結(jié)構(gòu)解釋量前三位的理化指標(biāo)分別為：有效磷(貢獻(xiàn)量為32.4%，P=0.002)、銨態(tài)氮(貢獻(xiàn)量為14.8%，P=0.01)和硝態(tài)氮(貢獻(xiàn)量為16.2%，P=0.048)。

3 討論

3.1 長期施肥對黑土理化性質(zhì)的影響

有機肥配施無機肥和單施無機肥處理極顯著增加土壤中全氮、硝態(tài)氮、有效磷和有機質(zhì)等的含量?？傮w上來說，施肥改變了土壤理化性質(zhì)。對比4個處理1980年的土壤基礎(chǔ)pH值(7.22)和施肥36年后的pH值，可知多年的施肥引起了土壤的嚴(yán)重酸化，應(yīng)采取優(yōu)化施肥措施緩解酸化。NPK處理的pH (5.89)明顯低于CK、M和MNPK三個處理，說明施無機肥導(dǎo)致了土壤的酸化，而有機肥或者有機無機配施可以減緩?fù)寥纏H值的降低。酸性土壤pH因為施入糞肥而增加土壤pH，這可能是由于在糞肥中存在著大量的碳酸鹽和碳酸氫鹽[29,31]，另外，帶有羥基和酚羥基的有機酸類在土壤中可以起到緩沖作用，而增加土壤pH[32]。

3.2 長期施肥對黑土中ITS基因豐度的影響

通過qPCR分析結(jié)果顯示，不施肥的ITS基因豐度明顯小于無機肥處理，顯著大于單施有機肥，略高于有機無機配施。這些結(jié)果說明施用無機肥會提高土壤中真菌的數(shù)量，施用有機肥減少了土壤中真菌數(shù)量，而有機無機配施的真菌數(shù)量與不施肥處理無顯著差異。因為真菌適應(yīng)生存的pH范圍較寬(pH 5～9)[33]，且較低的pH相對更有利于真菌的繁殖生長[34]，通過土壤pH分析得知，施用無機肥明顯的降低了土壤pH，4個處理中施無機肥的土壤pH最低，所以施無機肥的真菌數(shù)最高。有機無機配施相對單施無機肥來說，可減緩?fù)寥纏H降低，所以土壤真菌含量相對較低。通過ITS基因拷貝數(shù)與土壤理化性質(zhì)Pearson相關(guān)性分析，可知ITS基因拷貝數(shù)與土壤中的銨態(tài)氮含量顯著正相關(guān)，這個結(jié)果與Zhou等[20]在研究氮肥對真菌豐度影響時得到的結(jié)論一致。

圖5 不同施肥處理黑土土壤中真菌群落與土壤理化性質(zhì)的CCA分析Fig. 5 Relationship between dominant soil fungal communities and soil chemical characteristics presented by CCA

3.3 長期施肥對黑土中細(xì)菌與真菌比值的影響

本文研究發(fā)現(xiàn)長期施用有機肥或者有機無機配施會提高黑土土壤中細(xì)菌/真菌比值(MNPK處理的為26.91×104)，施無機肥(NPK處理的為9.90×104)會降低細(xì)菌/真菌比值。這種比值變化是因為土壤pH與細(xì)菌數(shù)量有相關(guān)性[35]，pH值較低時，細(xì)菌數(shù)量較??；無機肥處理的pH明顯小于有機無機配施處理，分析可知長期施用無機肥會降低黑土土壤中細(xì)菌數(shù)量，而有機無機配施會提高土壤中細(xì)菌數(shù)量(MNPK處理的16S rRNA基因拷貝數(shù)比NPK的高出2.40×1010/g,soil)。另外一個原因是pH值越低越有利于真菌的生長[34]。綜合長期施肥而使得細(xì)菌數(shù)量和真菌數(shù)量的變化，發(fā)現(xiàn)長期施用無機肥導(dǎo)致黑土由“細(xì)菌型”向“真菌型”轉(zhuǎn)變，而有機無機配施可以避免“細(xì)菌型”向“真菌型”的轉(zhuǎn)變。通過細(xì)菌/真菌比值和土壤理化性質(zhì)Pearson相關(guān)性分析，可知其與土壤pH有正相關(guān)性。

3.4 長期施肥對黑土中真菌 α 多樣性的影響

本研究發(fā)現(xiàn)雖然不施肥處理、無機肥處理、有機無機配施處理的真菌α多樣性總體上沒有顯著性差異，但是有機無機配施相對于無機肥有略提高真菌的多樣性，無機肥有略降低土壤真菌的多樣性的趨勢，這與Zhou等[20]發(fā)現(xiàn)長期施用無機氮肥會降低土壤真菌的α多樣性相一致。而Lou等[36]通過PCRDGGE方法研究報道無機肥和有機肥均會顯著增加真菌多樣性，尤其是有機肥，結(jié)果的不同很可能是因為試驗方法和土壤類型的差別。

3.5 長期施肥對黑土中真菌組成的影響

Ascomycota的相對豐度在不施肥處理(45.35%)，有機無機配施處理(50.93%)和無機肥(56.16%)處理中大于單施有機肥處理(35.47%)的相對豐度，說明Ascomycota對無論有機肥和無機肥均不穩(wěn)定。有機無機配施處理和無機肥處理中Ascomycota的相對豐度增加，是因為土壤中施入含氮的肥料使土壤中氮含量相對高，而土壤中施入適量氮能增加Ascomycota的百分比[13]。Zygomycota的相對豐度在單有機肥處理中要高于不施肥、有機無機配施和單施無機肥處理，有機糞肥提高了土壤中Zygomycota的含量，Neher等[37]研究發(fā)現(xiàn)Zygomycota的相對豐度與糞肥具有相關(guān)性。Cochliobolus在不施肥處理和單施無機肥處理均有較高的相對豐度，而在施入有機肥的MNPK處理和M處理中均大幅度降低；Exophiala和Alternaria有類似變化趨勢。通過門和屬兩個水平的相對豐度分析得知，長期施用有機肥和無機肥會改變土壤中真菌組成的相對豐度。

絲狀子囊菌屬Cochliobolus具有高度的侵染性，其病原菌種對宿主植物具有特異性[38]，引起作物病害，在M處理(0.39%)和MNPK處理(0.41%)含量均較低，而NPK處理(2.08%)比MNPK處理高1.67%，說明有機無機配施可以減少Cochliobolus的相對豐度，降低作物病害發(fā)生率。具有腐生性和病原性的真菌Exophiala[39]在無機肥處理中含量顯著高于有機無機配施處理；雖然我們沒有監(jiān)測作物的發(fā)病情況，但從真菌的組成也可以初步分析得到有機肥無機肥配施具有潛在減少作物土傳病害發(fā)生的作用。另外，通過土壤理化性質(zhì)與真菌組成的CCA分析，可以知道土壤理化性質(zhì)影響著土壤真菌群落結(jié)構(gòu)，尤其是土壤的有效磷和銨態(tài)氮濃度是兩個重要因素。影響真菌群落結(jié)構(gòu)的主要理化性質(zhì)均為速效養(yǎng)分，這可能是和土壤養(yǎng)分物質(zhì)存在的形式有關(guān)系，土壤中養(yǎng)分物質(zhì)存在形式又與地理環(huán)境、氣候溫度有相關(guān)性，例如，寒溫帶土壤中氮素化合物33%～42%的氮是以氨基酸形式存在，大約有10%的氨基酸水解成銨態(tài)氮[40]。

4 結(jié)論

本研究通過對36年連續(xù)施用有機肥和無機肥土壤中真菌群落特征的研究，發(fā)現(xiàn)連續(xù)施用無機肥會導(dǎo)致土壤酸化，改變真菌群落結(jié)構(gòu)，降低土壤真菌的多樣性；而有機肥無機肥配施減緩?fù)寥浪峄?，提高土壤真菌的多樣性。從長期有機無機配施對真菌組成影響的分析可知，長期有機無機配施的真菌組成結(jié)構(gòu)更有利于降低土傳病害發(fā)生率。長期施用無機肥已導(dǎo)致黑土中真菌的數(shù)量增加，致使土壤由“細(xì)菌型”轉(zhuǎn)化為“真菌型”，而有機無機配施可以改變這種變化趨勢。CCA分析明確有效磷、銨態(tài)氮和硝態(tài)氮是影響土壤真菌群落結(jié)構(gòu)變化重要土壤理化因素，同時土壤細(xì)菌/真菌比值與土壤pH正相關(guān)。

[1]國家統(tǒng)計局黑龍江調(diào)查總隊.2014年黑龍江統(tǒng)計年鑒[M].北京:中國統(tǒng)計出版社,2014. Heilongjiang survery office of national bureau of statistics of China. Heilongjiang statistical year book[M].Beijing:China Statistica Press,2014.

[2]Wei D,Yang Q,Zhang JZ,et al.Bacterial communities structure and diversity in ablack soil as affected by long-term fertilization[J]. Pedosphere,2008,18(5):582–592.

[3]Zhou J,Guan D,Zhou B,et al.Influence of34-years of fertilization on bacterial communities in an intensively cultivated black soil in northeast China[J].Soil Biology and Biochemistry,2015,90:42–51.

[4]陶磊,褚貴新,劉濤,等.有機肥替代部分化肥對長期連作棉田產(chǎn)量、土壤微生物數(shù)量及酶活性的影響[J].生態(tài)學(xué)報,2014,34(21): 6137–6146. Tao L,Chu GX,Liu T,et al.Impacts of organic manure partial substitution for chemical fertilizer on cotton yield,soil microbial community and enzyme activities in mono-cropping system in drip irrigation condition[J].Acta Ecologica Sinica,2014,34(21): 6137–6146.

[5]Zhao J,Ni T,Li J,et al.Effects of organic-inorganic compound fertilizer with reduced chemical fertilizer application on crop yields, soil biological activity and bacterial community structure in a rice–wheat cropping system[J].Applied Soil Ecology,2016,99: 1–12.

[6]Shen QR,Tan JF,Qian XQ.An introduction to soil and fertilizer [J].Beijing:Higher Education Press,2001:4–6,270–271.

[7]Krashevska V,Klarner B,Widyastuti R,et al.Impact of tropical lowland rainforest conversion into rubber and oil palm plantations on soil microbial communities[J].Biology and Fertility of Soils,2015, 51:697–705.

[8]劉益仁,郁潔,李想,等.有機無機肥配施對麥-稻輪作系統(tǒng)土壤微生物學(xué)特性的影響[J].農(nóng)業(yè)環(huán)境科學(xué)學(xué)報,2012,31(5):989–994. Liu YR,Yu J,Li X,et al.Effects of combined application of organic and inorganic fertilizers on soil microbiological characteristics in a wheat-rice rotation system[J].Journal of Agro-Environment Science, 2012,31(5):989–994.

[9]Van der Putten WH,Klironomos JN,Wardle DA.Microbial ecology of biological invasions[J].International Society for Microbial Ecology,2007,1:28–37.

[10]Per?oh D.Plant-associated fungal communities in the light of meta’omics[J].Fungal Diversity,2015,75(1):1–25.

[11]Wallenstein MD,Mcnulty S,Fernandez IJ,et al.Nitrogen fertilization decreases forest soil fungal and bacterial biomass in three long-term experiments[J].Forest Ecology and Management,2006, 222(1–3):459–468.

[12]Edwards IP,Zak DR,Kellner H,et al.Simulated atmospheric N deposition alters fungal community composition and suppresses ligninolytic gene expression in anorthern hard wood forest[J].PLoS One,2011,6:e20421.

[13]Paungfoo-Lonhienne C,Yeoh YK,Kasinadhuni NR,et al.Nitrogen fertilizer dose alters fungal communities in sugar cane soil and rhizosphere[J].Scientific Report,2015,5:8678.

[14]Wu T.Can ectomycorrhizal fungi circumvent the nitrogen mineralization for plant nutrition in temperate forest ecosystems[J]. Soil Biology and Biochemistry,2011,43:1109–1117.

[15]Gryndler M,Larsen J,HrselováH,et al.Organic and mineral fertilization,respectively,increase and decrease the development of external mycelium of arbuscular mycorrhizal fungi in along-term field experiment[J].Mycorrhiza,2006,16(3):159–166.

[16]M?der P,Fliessbach A,Dubois D,et al.Soil fertility and biodiversity in organic farming[J].Science,2002,296(5573):1694–1697.

[17]Lang J,Hu J,Ran W,et al.Control of cotton Verticillium wilt and fungal diversity of rhizosphere soils by bio-organic fertilizer[J]. Biology and Fertility of Soils,2011,48(2):191–203.

[18]Kiely PD,Call DF,Yates MD,et al.Anodic biofilms in microbial fuel cells harbor low numbers of higher-power-producingbacteria than abundant genera[J].Applied Microbiology and Biotechnology, 2010,88(1):371–380.

[19]Shokralla S,Spall JL,Gibson JF,et al.Next-generation sequencing

technologies for environmental DNA research[J].Molecular Ecology,2012,21(8):1794–1805.

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

[21]陳劍山,鄭服叢.ITS序列分析在真菌分類鑒定中的應(yīng)用[J].安徽農(nóng)業(yè)科學(xué),2007,35(13):3785–3786. Chen JS,Zheng FC.Application of ITS sequences in fungi classification and identification[J].Journal of Anhui Agricultural Sciences,2007,35(13):3785–3786.

[22]燕勇,李衛(wèi)平,高雯潔,等.rDNA-ITS序列分析在真菌鑒定中的應(yīng)用[J].中國衛(wèi)生檢驗雜志,2008,18(10):1958–1961. Yan Y,Li WP,Gao WJ,et al.Application of rDNA ITS sequence analysis in fungus identification[J].Chinese Journal of Health Laboratory Technology,2008,18(10):1958–1961.

[23]魯如坤.土壤農(nóng)業(yè)化學(xué)分析方法[M].北京:中國農(nóng)業(yè)科技出版社, 1999:431–472. Lu RK.Soil and agricultural chemistry analysis[M].Beijing:China Agricultural Science and Technology Press,1999:431–472.

[24]馬琳.玉米施氮量與土壤硝態(tài)氮、銨態(tài)氮的動態(tài)變化[D].哈爾濱:東北農(nóng)業(yè)大學(xué)碩士學(xué)位論文,2015. Ma L.Nitrogen fertilizer rates on maize and the dynamic changes of nitrate and ammonium in soil[D].Harbin:Master Dissertation of Northeast Agricultural University,2015.

[25]Fierer N,Lauber CL,Ramirez KS,et al.Comparative metagenomic, phylogenetic and physiological analyses of soil microbial communities across nitrogen gradients[J].The ISME Journal,2012, 6:1007–1017.

[26]Orgiazzi A,Lumini E,Nilsson RH,et al.Unravelling soil fungal communities from different Mediterranean land-use back grounds[J]. PloS One,2012,7:e34847.

[27]Schoch CL,Seifert KA,Huhndorf S,et al.Nuclear ribosomal internal transcribed spacer(ITS)region as auniversal DNA barcode marker for fungi[J].Proceedings of the National Academy of Sciences,2012,109:6241–6246.

[28]Chinnadurai C,Gopalaswamy G,Balachandar D.Long term effects of nutrient management regimes on abundance of bacterial genes and soil biochemical processes for fertility sustainability in asemi-arid tropical Alfisol[J].Geoderma,2014,232:563–572.

[29]Ding JL,Jiang X,Ma MC,et al.Effect of35years inorganic fertilizer and manure amendment on structure of bacterial and archaeal communities in black soil of northeast China[J].Applied Soil Ecology,2016,105:187–195.

[30]周晶,姜昕,周寶庫,等.長期施用尿素對東北黑土中氨氧化古菌群落的影響[J].中國農(nóng)業(yè)科學(xué),2016,49(2):294–304. Zhou J,Jiang X,Zhou BK,et al.Effects of long term application of urea on ammonia oxidizing archaea community in black soil in northeast China[J].Scientia Agricultura Sinica,2016,49(2): 294–304.

[31]Whalen JK,Chang C,Clayton GW,et al.Cattle manure amendments can increase the pH of acid soils[J].Soil Science Society of America,2000,64(3):962–966.

[32]Garc?á-Gil JC,Ceppi SB,Velasco MI,et al.Long-term effects of amendment with municipal solid waste compost on the elemental and acidic functional group composition and pH-buffer capacity of soil humic acids[J].Geoderma,2004,121(1–2):135–142.

[33]Rousk J,B??th E,Brookes PC,et al.Soil bacterial and fungal communities across apH gradient in an arable soil[J].The ISME Journal,2010,4:1340–1351.

[34]Rousk J,Brookes PC,B??th E.Contrasting soil pH effects on fungal and bacterial growth suggest functional redundancy in carbon mineralization[J].Applied and Environmental Microbiology,2009, 75(6):1589–1596.

[35]Ahn JH,Song J,Kim BY,et al.Characterization of the bacterial and archaeal communities in rice field soils subjected to long-term fertilization practices[J].Journal of Microbiology,2012,50(5): 754–765.

[36]Luo P,Han X,Wang Y,et al.Influence of long-term fertilization on soil microbial biomass,dehydrogenase activity,and bacterial and fungal community structure in abrown soil of northeast China[J]. Annals of Microbiology,2015,65:533–542.

[37]Neher DA,Weicht TR,Bates ST,et al.Changes in bacterial and fungal communities across compost recipes,preparation methods, and composting times[J].PLoS One,2013,8(11):e79512.

[38]Condon BJ,Leng Y,Wu D,et al.Comparative genome structure, secondary metabolite,and effector coding capacity across cochliobolus pathogens[J].PLoS Genetics,2013,9(1):e1003233.

[39]Yu L,Nicolaisen M,Larsen J,et al.Organic fertilization alters the community composition of root associated fungi in Pisum sativum[J]. Soil Biology and Biochemistry,2013,58(8):36–41.

[40]Sowden FJ,Chen Y,Schnitzer M.The nitrogen distribution in soils formed under widely differing climatic conditions[J].Geochim Cosmochim Acta,1977,41(10):1524–1526.

Structure of soil fungal communities under long-term inorganic and organic fertilization in black soil of Northeast China

DING Jian-li1,2,JIANG Xin1,2*,MA Ming-chao1,2,GUAN Da-wei1,2,ZHAO Bai-suo2,WEI Dan3,CAO Feng-ming1,2,LI Li1,2,LI Jun1,2*
[1 Institute of Agricultural Resources and Regional Planning, Chinese Academy of Agricultural Sciences, Beijing 100081, China; 2 Laboratory of Quality &Safety Risk Assessment for Microbial Products (Beijing), Ministry of Agriculture, Beijing 100081, China; 3 Institute of Soil Fertility and Environmental Sources, Heilongjiang Academy of Agricultural Sciences, Harbin 150086, China]

【Objectives】The abundance,diversity and composition of fungal community were identified,and the main relationship between soil fungi community variation and soil fertilizer would be revealed in this paper.【Methods】Based on36-years’fertilization experiment in Heilongjiang Academy of Agricultural Sciences,four treatments were selected to study:no fertilizer(CK);manure(M);chemical N,P and Kfertilizer(NPK); and manure with NPK(MNPK).The ITS gene was applied as target gene,using Illumina Miseq sequencing and quantitative PCR method.The fungal abundance,diversity and composition in black soil were analyzed.The relationship of fertilization and fungi community was calculated.【Results】Long term application of chemical fertilizers reduced soil pH and led to soil acidification,while manure did not show same impact.ITS gene abundance in ratio in MNPK(26.91×104)was significantly higher than that in NPK,and the ratio was in order of MNPK>M>CK>NPK.NPK treatment was more than in MNPK.Bacteria-to-fungi ratio was positively correlated with soil pH.Alpha diversity of fungi community of MNPK was slightly greater than that of NPK. Ascomycota and Zygomycota were the top two majority fungi found in black soil.The relative abundance of the fungi community was significantly different among the four fertilization treatments.The relative abundanceof Ascomycota in CK was45.35%,in MNPK and NPK were50.93%and56.16%,respectively.Manure reduced the relative abundance of pathogenic fungi.The relative abundances of highly aggressive Cochliobolus in MNPK (0.41%)and M(0.39%)were significantly lower than those of CK(3.25%)and NPK(2.08%).Canonical correlation analysis indicated that the selected soil properties could explain73.3%of the variation,and fungal communities were closely positively correlated with soil available P(Contribution=32.4%,P=0.002),NH4+-N concentration(Contribution=14.8,P=0.01)and NO3–-N concentration(Contribution=16.2%,P=0.048).【Conclusions】Long-term fertilization affects fungal abundance,diversity and community characteristics of black soil.Compared with chemical fertilizers,manure alone and manure combined with chemical NPK fertilizers will improve the fungal communities structures,decrease fungal abundance and increase diversity activity, enhance the soil pH,slowed down soil acidification.Among the soil fertility factors,the soil available P,NH4+-N and NO3–-N concentrations are the three dominant factors determining the variation of fungal communities structures.

fertilization;fungal communities;high-throughput sequencing;soil available P;soil NH4+-N; soil NO3–-N

2016–09–18接受日期：2017–03–01

國家自然科學(xué)基金項目（41573066）；國家“973”計劃（2015CB150506）；農(nóng)業(yè)微生物產(chǎn)品質(zhì)量安全風(fēng)險評估專項任務(wù)（GJFP201601202，GJFP201601203）；國家“863”計劃（2013AA102802-04）；國家現(xiàn)代農(nóng)業(yè)產(chǎn)業(yè)技術(shù)體系建設(shè)專項（CARS-04）資助。

丁建莉（1981—），女，遼寧凌源人，博士研究生，主要從事環(huán)境微生物與農(nóng)業(yè)生態(tài)方面的研究。E-mail：lndjl@126.com *通信作者E-mail：jiangxin@caas.cn；E-mail：lijun01@caas.cn