牟紅霞, 張文文, 劉秉儒

(1.寧夏大學(xué) 西北退化生態(tài)系統(tǒng)恢復(fù)與重建教育部重點(diǎn)實(shí)驗(yàn)室, 銀川 750021;2.北方民族大學(xué) 生物科學(xué)與工程學(xué)院, 銀川 750021)

紫花苜蓿(Medicagosativa)屬于多年生草本植物，且這種植物的根系比較發(fā)達(dá)、抗鹽堿、擁有較好的適應(yīng)能力、能夠防風(fēng)固沙、改良土壤等多種特性[1]。因其具有高產(chǎn)量、高品質(zhì)，可用于放牧、青貯飼料和干草的制備等多種利用方式，因此將其作為牧草栽培中的優(yōu)良品種。紫花苜蓿是我國種植面積居首位的牧草，具有十分重要的生態(tài)和經(jīng)濟(jì)效應(yīng)。關(guān)于紫花苜蓿生長年限對土壤理化性狀和肥力研究均有報(bào)道，比如，有研究已證明，生長年限對土壤理化性狀和肥力均有影響，除pH值外不同種植年限間有著一定的差異性[2-3]，土壤中有機(jī)質(zhì)和全氮的含量隨種植年限增加而增長，但其增長率會(huì)隨種植年限增加而下降，這表示其增幅具有一定的周期性[4-5]。同時(shí)苜蓿地土壤中養(yǎng)分增長年限會(huì)因?yàn)閰^(qū)域的不同而有所不同，在中國黃土高原地區(qū)，通常其土壤有機(jī)質(zhì)、全氮增長期表現(xiàn)在5～10 a[6]。但是在一些干旱地區(qū)，種植4 a左右，土壤養(yǎng)分含量相對來說比較高[7]。對“黑土灘”不同建植年限栽培草地的研究則發(fā)現(xiàn)土壤養(yǎng)分變化表現(xiàn)為2 a一個(gè)周期的變化規(guī)律[8]。土壤養(yǎng)分增加幅度和周期受到土壤微生物的調(diào)控，但是對苜蓿種植年限的土壤微生物群落結(jié)構(gòu)變化及演變規(guī)律與影響因素知之甚少。

微生物是整個(gè)土壤生態(tài)系統(tǒng)之中最為關(guān)鍵的分解者，參與到了土壤碳氮養(yǎng)分循環(huán)的過程之中，是土壤中肥力的重要活性因子[9-10]。真菌是土壤中微生物的一個(gè)主要成分，大約占微生物總量的90%[11]。而微生物量快速周轉(zhuǎn)過程中逐步釋放的供植物生長的有效營養(yǎng)，主要受真菌菌群的調(diào)控。同時(shí)在植物有機(jī)體分解的初期過程中，真菌是比較活躍的，因此真菌在生態(tài)系統(tǒng)營養(yǎng)循環(huán)中具有關(guān)鍵作用[12]。在降雨量只有200 mm左右的引黃灌區(qū)，紫花苜蓿種植有著得天獨(dú)厚的優(yōu)勢，使其產(chǎn)量更高、品質(zhì)更好，在這樣的水土條件下，微生物真菌的群落結(jié)構(gòu)和功能與其他地區(qū)有怎樣的差異尚不清楚。因此，本研究從土壤微生物中真菌多樣性特征入手，通過不同種植年限苜蓿地土壤肥力變化與真菌群落多樣性變化之間的某種規(guī)律，尋找影響土壤肥力的主要真菌菌群，為引黃灌區(qū)紫花苜蓿連種提供一個(gè)穩(wěn)定的土壤養(yǎng)分循環(huán)系統(tǒng)。

因其真菌的種類非常多，而且個(gè)體間的多態(tài)性不是非常明顯，因此采用傳統(tǒng)的方法對真菌進(jìn)行正確的分類存在較大的困難[13]。本研究使用高通量測序的方法對其土壤中的真菌多樣性進(jìn)行全面的分析。深度測序的性能是高通量測序技術(shù)的優(yōu)勢[14]，同時(shí)也給微生物方面的研究提供了一定技術(shù)層面的便利。本研究以寧夏賀蘭山人工紫花苜蓿農(nóng)場種植時(shí)間為1～6 a的紫花苜蓿地為研究對象，研究其土壤養(yǎng)分分布特征，并使用18SrDNA擴(kuò)增子測序技術(shù)分析土壤中真菌結(jié)構(gòu)及多樣性特征。進(jìn)而從土壤真菌群落變化和養(yǎng)分動(dòng)態(tài)變化等角度篩選探索最佳的牧草種植期，并為西北地區(qū)紫花苜蓿土壤養(yǎng)分管理和大面積栽培及高質(zhì)量生產(chǎn)提供科學(xué)指導(dǎo)[15]。

1 研究區(qū)概況

研究區(qū)位于寧夏銀川市賀蘭山的人工苜蓿農(nóng)場(38°30′N，106°06′E)，海拔高度為1 135 m，位于我國西北內(nèi)陸，是中溫帶干旱區(qū)[15]。晝夜溫差較大，該地區(qū)年平均氣溫8.4℃，年最低氣溫為-20℃，極端最高氣溫為38℃，年平均降水量為189.9～216.1 mm，年蒸發(fā)量為2 250.0 mm左右。年平均日照時(shí)間為3 075.5 h，無霜期為160 d。人工苜蓿農(nóng)場從2012年開始，每一年都有新種植的同種苜蓿。

2 材料與方法

2.1 樣品采集

本研究以種植年限為1～6 a(2012—2017年)的紫花苜蓿地為研究對象，苜蓿地的灌溉條件、灌溉制度、氣候降水及土壤類型一致，均為地形平坦的灌區(qū)土壤環(huán)境，試驗(yàn)設(shè)置的6個(gè)不同年限處理中，每個(gè)處理設(shè)置3個(gè)面積為30 m×30 m的重復(fù)，在每個(gè)重復(fù)樣地使用土壤鉆按照5點(diǎn)混合法取0—20 cm深度的土壤[15]，每個(gè)土樣采集1份，在試驗(yàn)樣地將一小部分裝入經(jīng)高壓滅菌后的凍存管，貼上標(biāo)簽，立即放入-4℃的保溫箱。剩余樣品裝入自封袋，并對其進(jìn)行相應(yīng)的標(biāo)記。隨后將土樣帶回實(shí)驗(yàn)室，凍存管轉(zhuǎn)移至-80℃冰箱保存，用于真菌群落的測定。自封袋中的樣品在實(shí)驗(yàn)室自然風(fēng)干，并將土壤中植物碎屑和細(xì)根等雜質(zhì)剔除，按要求過篩后用于測定土壤的各項(xiàng)理化指標(biāo)。

2.2 土壤理化性質(zhì)的測定

土壤的pH值、電導(dǎo)率是通過水土比2.5∶1的懸液，用PHS-3C酸度計(jì)和便攜式電導(dǎo)率儀測定。土壤全氮采用elementer元素分析儀進(jìn)行測定。采用堿解擴(kuò)散法對堿解氮進(jìn)行測定。用重鉻酸鉀氧化—外加熱法測定土壤中的有機(jī)碳。全磷采用NaOH熔融—鉬銻抗比色法。速效磷是用NaHCO3浸提—鉬銻抗比色方法測定的[16]。最后采用SPSS 19.0對土壤理化數(shù)據(jù)進(jìn)行方差分析和多重比較。

2.3 土壤基因組DNA提取及HiSeq測序

土壤在碾碎過80目篩后，第一步是先使用SDS對樣本中的DNA進(jìn)行提取，然后使用瓊脂糖凝膠電泳測定樣本中DNA的程度和濃度，取出一定量的樣品，然后將其置于離心管內(nèi)，使用無菌水稀釋樣品至1 ng/μl[17]。以稀釋后的基因組DNA為測序模板，根據(jù)測序區(qū)域的選擇，使用Barcode的特異引物528F和706R，New England Biolabs公司的Phusion?High-Fidelity PCR Master Mix with GC Buffer和高效高保真酶進(jìn)行PCR處理，確保擴(kuò)增效率和準(zhǔn)確性，之后通過HiSeq 2500 PE 250進(jìn)行上機(jī)測序[18]。

2.4 生物信息學(xué)分析

下機(jī)數(shù)據(jù)通過預(yù)處理，將一些質(zhì)量較低的reads移除，然后將成對reads拼接成一條序列，為得到高質(zhì)量的reads，對tags兩端的barcode序列及引物序列進(jìn)行去除，去除掉嵌合體以及短序列之后獲得clean tags。通過對經(jīng)拼接過濾后的clean tags與物種注釋數(shù)據(jù)庫比對來檢測嵌合體序列，同時(shí)去除其中的嵌合體序列[19]，得到最終的Effective Tags，對所有樣本的全部Effective Tags進(jìn)行聚類，默認(rèn)以97%的一致性(Identity)將序列聚類成為OTU(Operational Taxonomic Units)，同時(shí)依據(jù)其算法原則篩選出OTU中出現(xiàn)頻數(shù)最高的序列作為OTUs的代表序列，對OTU代表序列進(jìn)行物種注釋，使用Qiime軟件(Version 1.9.1)之中的blast方法[20]與Unit數(shù)據(jù)庫對物種進(jìn)行注釋分析，而且對于每個(gè)分類水平：kingdom(界)，phylum(門)，class(綱)，order(目)，family(科)，genus(屬)，species(種)進(jìn)行分類統(tǒng)計(jì)，統(tǒng)計(jì)其群落組成以及各個(gè)樣本的豐度信息。

2.5 多樣性分析

Alpha多樣性分析是對單個(gè)樣品中物種多樣性分析，基于OTU的結(jié)果，同時(shí)利用Qiime軟件(Version 1.9.1)計(jì)算出菌群豐度指數(shù)(Chao1，ACE)和菌群多樣性指數(shù)的大小(Shannon，Simpson)[21]；Beta多樣性是比較和分析不同樣本的微生物群落構(gòu)成。同樣是基于OTU的結(jié)果計(jì)算Unifrac距離，之后使用OTU的豐度信息對Unifrac距離進(jìn)一步構(gòu)建Weighted Unifrac距離[22]。最后，通過多變量統(tǒng)計(jì)學(xué)方法PCA主成分分析、PCoA主坐標(biāo)分析和Beta多樣性指數(shù)組間的差異分析，從而發(fā)現(xiàn)不同樣本(組)之間的差異。

2.6 環(huán)境因子關(guān)聯(lián)分析

研究環(huán)境因子與物種之間，環(huán)境因子與物種豐富度(alpha多樣性)之間的相關(guān)性，目前常用的有Spearman相關(guān)性分析和Mantel test分析。Spearman相關(guān)性分析以Spearman相關(guān)系數(shù)作為量度，Spearman相關(guān)系數(shù)又被叫做秩相關(guān)系數(shù)。其基本原理也就是兩變量所具有的秩次大小不同進(jìn)行線性分析，對于數(shù)據(jù)的原始變量分布沒有進(jìn)行要求，這是一種非參數(shù)類型的分析辦法，其適用范圍比較廣泛；Mantel test是對兩個(gè)矩陣相關(guān)關(guān)系的檢驗(yàn)，總的來說，它是兩個(gè)不同總體之間的相關(guān)關(guān)系，多用于生態(tài)學(xué)上，Mantel test用于計(jì)算環(huán)境因子和微生物群落數(shù)據(jù)的相關(guān)性，是對兩個(gè)矩陣相關(guān)關(guān)系的檢驗(yàn)。

3 結(jié)果與分析

3.1 不同種植年限紫花苜蓿土壤基本理化指標(biāo)分析

研究區(qū)pH值在8.5～8.7，屬于堿性土壤，第1年和第6年存在顯著差異，但相差數(shù)值僅為0.11，所以長期種植苜蓿對土壤酸堿度影響很小。對苜蓿地土壤樣本進(jìn)一步進(jìn)行主要養(yǎng)分分析得到：電導(dǎo)率、有效磷和全磷含量均無顯著變化，有機(jī)碳、全氮和堿解氮含量均呈先降后回升的趨勢，且在第5年達(dá)到最大，說明隨著種植年限的增加，苜蓿地土壤養(yǎng)分變化可能表現(xiàn)為循環(huán)往復(fù)的波浪式變化趨勢(表1)。

表1 不同種植年限紫花苜蓿土壤理化指標(biāo)

3.2 紫花苜蓿土壤真菌OTU水平分析

為了準(zhǔn)確地認(rèn)知各樣本OTU聚類情況和注釋情況，對OTU聚類及注釋結(jié)果進(jìn)行了綜合統(tǒng)計(jì)和分析，其中Total Tags是每個(gè)樣本的總Tags數(shù)目，用于OTU聚類等后續(xù)分析的有效數(shù)據(jù)；Taxon Tags是用于構(gòu)建OTUs并且獲得注釋信息的Tags數(shù)目；Unclassified Tags指的是沒有獲得注釋信息的Tags數(shù)目；Unique Tags是頻數(shù)為1，并且無法被聚類到OTUs的Tags數(shù)目，所以不對其進(jìn)行后續(xù)分析，詳細(xì)的統(tǒng)計(jì)數(shù)據(jù)如表2所示。在樣本之中隨機(jī)性的進(jìn)行抽取測序量的數(shù)據(jù)，對他們的物種數(shù)目進(jìn)行統(tǒng)計(jì)(也就是OTUs數(shù)目)，將測序數(shù)據(jù)量和其對應(yīng)的物種數(shù)來構(gòu)建曲線，即稀釋曲線，稀釋曲線能夠直觀地反映出測序數(shù)據(jù)量是否科學(xué)合理，同時(shí)能夠間接地反映出物種的豐富度，如圖1所示，曲線傾向于平坦，表明測序數(shù)據(jù)的科學(xué)合理，而更多的數(shù)據(jù)量只會(huì)產(chǎn)生少量新的物種。

表2 不同種植年限苜蓿地OTU聚類和注釋情況統(tǒng)計(jì)

圖1 土壤真菌稀釋曲線分析

3.3 不同種植年限苜蓿地土壤真菌群落組成

圖2和圖3分別展示了6組樣品，在綱水平和屬水平下豐度前10的真菌群落組成情況。在綱水平上，鑒定的豐度前10的菌占了總數(shù)中的73%，其中包括糞殼菌綱(Sordariomycetes)(44%)、座囊菌綱(Dothideomycetes)(12%)、散囊菌綱(Eurotiomycetes)(6%)、盤菌綱(Pezizomycetes)(3%)、絲足蟲綱(Incertae)(4%)為優(yōu)勢菌，所占比例高達(dá)69%，不同菌隨種植年限的的變化呈現(xiàn)不同的變化規(guī)律(圖2)。在屬水平上，共鑒定出328個(gè)屬，豐度前10的菌占總數(shù)的27%，其中毛殼屬(Chaetomium)(8%)、鐮孢屬(Fusarium)(7%)、赤霉菌屬Gibberella(4%)、Lectera(1%)、小畫線殼屬(Monographella)(1%)，除毛殼屬隨著種植年限的不同有較大變化外，其余菌屬變化都較為穩(wěn)定(圖3)。

圖2 不同種植年限苜蓿地土壤真菌群落相對豐度前10真菌綱

圖3 不同種植年限苜蓿地土壤真菌群落相對豐度前10真菌屬

3.4 不同種植年限苜蓿地土壤真菌生物多樣性

3.4.1 單個(gè)樣品Alpha多樣性分析 對指示群落結(jié)構(gòu)改變的4個(gè)Alpha多樣性指數(shù)Shannon，Simpson，Chao1以及ACE指數(shù)展開組間差異性分析。表3結(jié)果表明，Simpson指數(shù)隨種植年限的不同并沒有發(fā)生顯著性變化，Shannon，Chao1和ACE指數(shù)均在種植時(shí)間為5 a時(shí)達(dá)到最大，2 a時(shí)達(dá)到最小，且差異達(dá)顯著水平，由此可知在種植時(shí)間為5 a時(shí)，土壤真菌菌群豐度和多樣性均達(dá)到最高，變化規(guī)律從大到小依次為5 a>6 a>1 a>4 a>3 a>2 a。

表3 不同種植年限苜蓿地土壤真菌Alpha多樣性指數(shù)

3.4.2 復(fù)雜樣品的真菌群落Beta 多樣性分析 Beta多樣性是指利用Weighted Unifrac距離對兩個(gè)樣品之間的相異系數(shù)進(jìn)行衡量，其數(shù)值越小，表示兩個(gè)樣本在物種多樣性方面差異性越小。Beta多樣性分析表明，種植時(shí)間為1 a和2 a，3 a和4 a，2 a和4 a之間的相異系數(shù)從大到小分別為0.554，0.489，0.482，這3組樣本物種多樣性差異達(dá)到最大；而2 a和3 a，4 a和5 a，2 a和5 a之間相異系數(shù)從小到大分別為0.26，0.298，0.364，因此可以得到種植時(shí)間為2 a和3 a時(shí)，其物種多樣性組成相似性最高(圖4)。

圖4 加權(quán)Unifrac的熱圖

3.5 土壤環(huán)境因子與真菌群落組成相關(guān)性分析

對屬水平下真菌群落組成與土壤環(huán)境因子展開Spearman相關(guān)性分析，結(jié)果表明：pH值、AN和SOC與一部分真菌間存在顯著或極顯著相關(guān)性。pH與Gibellulopsis、鏈格孢屬(Alternaria)具有極顯著正相關(guān)關(guān)系，與Lectera、小畫線殼屬(Monographella)具有顯著負(fù)相關(guān)關(guān)系，與金孢子菌屬(Chrysosporium)呈顯著正相關(guān)關(guān)系；AN與Lectera具有極顯著負(fù)相關(guān)關(guān)系，與木霉屬(Trichoderma)呈顯著正相關(guān)關(guān)系，與金孢子菌屬(Chrysosporium)呈顯著負(fù)相關(guān)關(guān)系；SOC與綠僵菌屬(Metarhizium)、小畫線殼屬(Monographella)、漆斑菌屬(Myrothecium)具有顯著負(fù)相關(guān)關(guān)系。EC、TP和TN與個(gè)別屬類存在顯著相關(guān)性(圖5)。

圖5 環(huán)境因子相關(guān)性分析熱圖

4 討 論

土壤真菌數(shù)量大，種類多，是檢測土壤質(zhì)量變化的重要指標(biāo)，在評價(jià)土壤生態(tài)系統(tǒng)等功能方面具有重要作用[23]。本研究通過對不同種植年限苜蓿地土壤在綱水平和屬水平下真菌群落組成中豐度占比前10的菌類分析發(fā)現(xiàn)，糞殼菌綱、座囊菌綱、散囊菌綱、盤菌綱、絲足蟲綱為優(yōu)勢菌綱，毛殼屬、鐮孢屬、赤霉菌屬、Lectera、小畫線殼屬為優(yōu)勢菌屬。其中糞殼菌綱占總綱類真菌的44%，而且隨著種植年限的增加，沒有大的波動(dòng)，糞殼菌是一種寄生在動(dòng)物糞便上的真菌，因此糞殼菌成為最具優(yōu)勢菌的原因可能是農(nóng)田管理過程中大量農(nóng)家肥施入的結(jié)果，其他優(yōu)勢菌類隨著種植年限的增加，表現(xiàn)為先減小后增大的變化趨勢，這與土壤理化指標(biāo)變化基本保持一致。在屬水平下，土壤真菌群落由于受pH值，AN和SOC等理化因子的影響，隨著種植年限的增加，各菌屬變化不再呈現(xiàn)規(guī)律性變化。相關(guān)的研究顯示，鐮刀屬、毛殼屬和木霉屬在土壤中是常見的纖維素分解者[24]，而本研究發(fā)現(xiàn)毛殼屬為6個(gè)樣本中最具優(yōu)勢的真菌，占總菌屬的8%，因此毛殼屬可能是影響苜蓿地土壤碳循環(huán)的主要真菌菌屬[25]。

土壤微生物群落多樣性通常用多樣性指數(shù)的大小來分析，多樣性指數(shù)越高則微生物群落多樣性越高[26]。本試驗(yàn)通過對指示群落結(jié)構(gòu)變化的Alpha多樣性指數(shù)進(jìn)行組間差異性分析，最終研究結(jié)果表明，Simpson指數(shù)不受種植時(shí)間的影響，而Shannon指數(shù)、Chao1指數(shù)和ACE指數(shù)表現(xiàn)出相同的變化規(guī)律，即均在種植2 a時(shí)達(dá)到最小，5 a時(shí)達(dá)到最大，且差異達(dá)顯著水平，因此連種苜蓿會(huì)使土壤中真菌菌群的多樣性處于一個(gè)較高的水平，這與Wal等[27]發(fā)現(xiàn)休耕地土壤真菌多樣性指數(shù)也處于較高水平的研究結(jié)果保持一致，畢明麗[28]、孫倩[29]等研究結(jié)果也支持這一結(jié)論。Beta多樣性分析結(jié)果顯示，隨著種植年限的增加，土壤真菌多樣性變化表現(xiàn)為由變幅較大到逐漸趨于平穩(wěn)的趨勢，各種植年限間相異系數(shù)可反映這一變化規(guī)律。

已有研究表明，土壤真菌多樣性和群落組成受多種因素影響，土壤環(huán)境是其中一個(gè)重要因素[30-31]，因此對屬水平下真菌群落組成與土壤環(huán)境因子進(jìn)行Spearman相關(guān)性分析，有助于探索影響不同種植年限苜蓿地土壤真菌群落結(jié)構(gòu)及多樣性的主要影響因子。結(jié)果表明土壤pH值，AN和SOC與一部分真菌間存在顯著或極顯著相關(guān)性，EC、TP和TN與個(gè)別真菌存在顯著相關(guān)性，AP與所有真菌不具有相關(guān)性。Liu等[32]研究發(fā)現(xiàn)，速效磷對土壤真菌群落結(jié)構(gòu)具有較大影響，且施磷肥會(huì)提高土壤真菌多樣性[33]，而本研究中，速效磷變化范圍為7.89～14.15，隨著種植年限的增加未發(fā)生顯著性變化，因此未能體現(xiàn)影響程度。Beare等[34]提出免耕少耕能增加真菌菌群的多樣性，增加不穩(wěn)定碳的固定積累，進(jìn)而影響土壤有機(jī)碳含量，這與本研究Spearman相關(guān)性分析結(jié)果一致。本研究中pH值隨著種植年限的變化第1年與第6年存在顯著性變化，且pH值與真菌屬水平下存在顯著或極顯著相關(guān)性，這與韓世忠等[35]研究發(fā)現(xiàn)土壤pH值會(huì)對真菌群落結(jié)構(gòu)組成造成影響的結(jié)果一致。

5 結(jié) 論

毛殼屬、鐮孢屬、赤霉菌屬、Lectera、小畫線殼屬這是在不同的種植年限苜蓿地土壤中的優(yōu)勢菌群，隨著種植年限的增加并無顯著性變化；多樣性指數(shù)分析結(jié)果顯示，土壤中的養(yǎng)分以及真菌群落組成和豐度都在種植5 a的時(shí)間達(dá)到峰值；Simpson指數(shù)隨種植年限的增加無顯著性變化；種植時(shí)間為1～6 a的苜蓿地土壤真菌的多樣性組成相似度較低，土壤pH值、速效氮和有機(jī)碳是影響真菌群落組成及多樣性的主要環(huán)境因子。