黔中山地不同種植方式茶園土壤真菌群落特征

張健，徐明，陳馳，張姣，陳進(jìn)

（1. 貴州大學(xué)生命科學(xué)學(xué)院真菌資源研究所，貴州 貴陽(yáng) 550025；2. 貴州貴茶（集團(tuán)）有限公司，貴州 貴陽(yáng) 550025）

微生物是土壤生態(tài)系統(tǒng)的重要生物組分，在土壤元素的生物地球化學(xué)循環(huán)中發(fā)揮著關(guān)鍵作用［1］。土壤微生物在保持土壤生態(tài)系統(tǒng)健康方面中發(fā)揮著不可替代的作用，它們可通過(guò)豐富多樣的生活方式、代謝途徑和復(fù)雜的生態(tài)網(wǎng)絡(luò)關(guān)系，實(shí)現(xiàn)土壤生態(tài)系統(tǒng)的物質(zhì)循環(huán)、污染物降解等生態(tài)功能［1－4］。真菌是全球生物多樣性的重要組成部分，土壤真菌通過(guò)共生－寄生的相互作用影響動(dòng)植物群落動(dòng)態(tài)［2］。土壤真菌作為監(jiān)測(cè)土壤質(zhì)量變化和健康狀況的一個(gè)重要指標(biāo)，在評(píng)價(jià)生態(tài)系統(tǒng)健康與維護(hù)生態(tài)系統(tǒng)平衡中也發(fā)揮了關(guān)鍵作用［1，5］。

我國(guó)有著悠久的植茶歷史與茶文化積淀。目前我國(guó)是茶園面積和茶產(chǎn)量、茶消費(fèi)量等位居全球最高的國(guó)家［6］。伴隨綠色可持續(xù)發(fā)展理念和茶產(chǎn)業(yè)多元途徑的開(kāi)發(fā)利用，我國(guó)茶園種植面積呈逐年增加的趨勢(shì)［7］。但我國(guó)茶產(chǎn)業(yè)發(fā)展也面臨著茶葉生產(chǎn)技術(shù)與科技創(chuàng)新、茶業(yè)專(zhuān)業(yè)技術(shù)人才與勞動(dòng)力緊缺等諸多因素的影響；鑒于我國(guó)不同茶區(qū)地理生態(tài)環(huán)境、人文、種植管理等方面差異明顯，亟待加強(qiáng)茶產(chǎn)業(yè)相關(guān)基礎(chǔ)與應(yīng)用的創(chuàng)新性研究［8］。近年來(lái)，以高原山地為主要地貌的貴州省，發(fā)揮其山地氣候條件適于茶樹(shù)生長(zhǎng)的優(yōu)勢(shì)，已發(fā)展成為我國(guó)茶園面積種植最大的省份。

土壤微生物是土壤中最活躍且多樣性最為豐富的組分之一，可作為反映茶園土壤質(zhì)量及評(píng)價(jià)土壤生態(tài)系統(tǒng)可持續(xù)性的重要生物學(xué)指標(biāo)［9］。土壤微生物在反映茶園環(huán)境與種植管理方面較為敏感，具有重要的生物指示作用［10］。目前有關(guān)茶園類(lèi)型、施肥模式、遮陰、修剪等方面對(duì)土壤細(xì)菌群落報(bào)道相對(duì)較多［11－15］，針對(duì)現(xiàn)代茶園、古茶園及不同培肥管理茶園的土壤真菌群落已有相關(guān)報(bào)道［10，16，17，22］；但有關(guān)不同茶園種植方式對(duì)土壤真菌影響的研究卻鮮有報(bào)道。本研究以黔中山地不同種植方式（新茶園、普通茶園和碾茶園）茶園土壤真菌群落結(jié)構(gòu)與功能為研究目標(biāo)，進(jìn)一步分析影響黔中山地茶園土壤真菌群落的主要環(huán)境因子，以期為該區(qū)茶園生態(tài)管理提供科學(xué)參考。

1 材料與方法

1.1 研究區(qū)域概況

調(diào)查的山地茶園位于貴州省貴陽(yáng)市花溪區(qū)久安鄉(xiāng)境內(nèi)，該區(qū)位于東經(jīng)106°35′～106°37′，北緯26°30′～26°32′；久安鄉(xiāng)地勢(shì)起伏大，屬中山山地地形，海拔1090～1400 m；年降雨量為1000～1100 mm，屬亞熱帶季風(fēng)性氣候，年平均氣溫14.8℃，無(wú)霜期260 d；土壤類(lèi)型主要為黃壤。久安鄉(xiāng)地下煤礦資源豐富，曾經(jīng)煤產(chǎn)業(yè)盛行；但由于地處阿哈湖畔上游，屬貴陽(yáng)市水源保護(hù)地水資源環(huán)境，于2010年開(kāi)始，陸續(xù)全部關(guān)停了境內(nèi)煤礦，依托久安鄉(xiāng)境內(nèi)豐富的古茶樹(shù)資源和適宜茶樹(shù)生長(zhǎng)的地理氣候條件轉(zhuǎn)型發(fā)展生態(tài)茶園［18］。

1.2 樣地設(shè)置及樣品采集

選取了新茶園（Tx）、普通茶園（Tp）和碾茶園（Tn）3種不同種植管理方式的茶園樣地，以周邊馬尾松林（Pinus massoniana，Mp）作為對(duì)照（表1），分別選取了2個(gè)樣地重復(fù)。其中，新茶園由山地自然植被轉(zhuǎn)化茶園種植，種植管理年限為7年，進(jìn)行定形修剪和有機(jī)肥培肥管理；普通茶園種植管理年限為11年，采用輕修剪和有機(jī)肥培肥管理；碾茶園種植管理年限為11年，但其在株高（相較普通茶園矮）、遮陰處理、修剪（平行或弧型修剪）和采摘以及施肥等種植管理方面差異明顯［15，19］。于2018年8月采集土壤樣品，用土鉆采集0～20 cm土層的六點(diǎn)混合土樣。土壤樣品經(jīng)去掉石頭、根系等雜質(zhì)后，全部過(guò)2 mm土篩，均勻混合后分為2份，一份經(jīng)室內(nèi)自然風(fēng)干后用于測(cè)定土壤理化性質(zhì)；另一份經(jīng)冷凍干燥保存于-80℃冰箱，用于高通量測(cè)序的土壤總DNA提取。

1.3 土壤理化性質(zhì)測(cè)定

pH值采用電位法（土水重量比為1∶2.5）測(cè)定；有機(jī)質(zhì)（SOM）采用重鉻酸鉀氧化－分光光度法測(cè)定；全氮（TN）采用半微量凱氏定氮法測(cè)定；全磷（TP）采用鉬銻抗比色法測(cè)定；有效磷（AP）采用雙酸浸提－鉬銻抗比色法測(cè)定；以上土壤指標(biāo)的測(cè)定參照魯如坤的《土壤農(nóng)業(yè)化學(xué)分析方法》［20］。不同種植茶園樣地基本信息描述見(jiàn)表1。

表1 調(diào)查樣地基本信息Table 1 Basic information on sampling sites

1.4 土壤真菌ITS擴(kuò)增子高通量測(cè)序

采用試劑盒（FastDNA?SPIN Kit for Soil）提取純化土壤總DNA。土壤真菌ITS擴(kuò)增子采用通用引物 為 ITS3F（5＇-GCATCGATGAAGAACGCAGC-3＇）和 ITS4R（5＇-TCCTCCG CTTATTGATATGC-3＇） 的核糖體RNA ITS2高變區(qū)進(jìn)行測(cè)序，樣品由上海美吉生物醫(yī)藥科技有限公司利用Illumina MiSeq PE300（2×300 bp）測(cè)序平臺(tái)完成測(cè)序。

1.5 數(shù)據(jù)處理

經(jīng)高通量測(cè)序得到的雙端序列原始數(shù)據(jù)，根據(jù)PE reads之間的overlap關(guān)系，將成對(duì)的reads拼接（merge）成一條序列，利用FLASH軟件對(duì)reads質(zhì)量和拼接效果進(jìn)行質(zhì)控過(guò)濾，再根據(jù)序列首尾兩端的barcode和引物序列區(qū)分樣品得到有效序列。通過(guò)Uparse軟件在相似水平為97%的OTU進(jìn)行聚類(lèi)分析；樣品序列按最小序列抽平，本研究中每個(gè)樣品按60540條序列進(jìn)行抽平，通過(guò)QIIME軟件獲取OTU豐度、α多樣性（包括Sobs、Shannon、Simpson、ACE和Chao等多樣性指數(shù)）、稀釋曲線(xiàn)和韋恩圖；采用Welch T法檢驗(yàn)不同茶園種植方式土壤真菌群落在屬水平上組間的顯著性差異，采用“fdr”法對(duì)P值進(jìn)行多重檢驗(yàn)校正，采用“Welch’s intverted”法計(jì)算95%的置信區(qū)間；基于非加權(quán)unifrac距離矩陣按照average算法進(jìn)行層級(jí)聚類(lèi)；采用加權(quán)標(biāo)準(zhǔn)化unifrac距離算法進(jìn)行非度量多維尺度（Non-metric multidimensional scaling，NMDS） 和相似性檢驗(yàn)（Analysis of similarities，ANOSIM）分析不同茶園種植方式土壤真菌群落樣本間的相似性。采用典范關(guān)聯(lián)分析（Canonical Correlation Analysis，CCA）揭示茶園土壤真菌群落與環(huán)境因子間作用關(guān)系。利用FUNGuild（http://www.funguild.org/）微生態(tài)工具按照營(yíng)養(yǎng)方式對(duì)土壤真菌功能進(jìn)行預(yù)測(cè)分析［21］。

2 結(jié)果與分析

2.1 不同種植方式茶園土壤真菌群落Alpha多樣性特征

土壤真菌ITS擴(kuò)增子高通量測(cè)序，共得到優(yōu)化有效序列總數(shù)為574809，物種稀釋曲線(xiàn)呈現(xiàn)出隨著測(cè)得序列數(shù)量的增加逐漸趨于平坦（圖1a），這表明測(cè)試有效序列數(shù)據(jù)信息量能夠反映土壤樣本真菌群落多樣性。按＞97%相似度進(jìn)行OTU注釋?zhuān)寥勒婢郝涔茶b定得到13門(mén)、54綱、137目、315科、633屬、1017種和2857 OTUs。其中，土壤真菌群落OTUs豐富度表現(xiàn)出（表2）：Tx（1636）＞Mp（1354）＞Tp（1281）＞Tn（1197），且共有OTUs和特有OTUs分別占比為47.32%和52.68%（圖1b）。不同種植方式茶園土壤真菌群落Alpha多樣性均表現(xiàn)出Tx最高和Tn最?。ū?）。

圖1 不同種植方式下茶園和馬尾松林土壤真菌的稀釋曲線(xiàn)（a）和韋恩圖（b）Fig. 1 Dilution curves （a） and Venn diagram （b） of soil fungi at different type tea plantations and surrounding P. massoniana forests

表2 不同種植方式下茶園和馬尾松林土壤真菌群落的分類(lèi)鑒定與多樣性指數(shù)Table 2 Fungal community classification and diversity indices of soils at different type tea plantations and surrounding P. massoniana forests

2.2 不同種植方式茶園土壤真菌群落組成結(jié)構(gòu)特征

茶園與對(duì)照馬尾松林的土壤真菌群落組成間存在明顯差異，且不同種植方式茶園土壤真菌群落間也存在一定差異。在門(mén)水平上，子囊菌門(mén)（Ascomycota）和擔(dān)子菌門(mén)（Basidiomycota）為不同種植方式茶園和馬尾松林土壤真菌群落的優(yōu)勢(shì)菌門(mén)（圖2a），但二者在組成結(jié)構(gòu)上存在明顯不同；其中，子囊菌門(mén)與擔(dān)子菌門(mén)間比例（A∶B）在茶園土壤真菌群落中約為（3.5～4.7）∶1，不同種植方式茶園土壤真菌群落的A∶B表現(xiàn)出：Tx（3.5∶1）＜Tp（3.8∶1）＜Tn（4.7∶1）；而在馬尾松林A∶B比例約為0.92∶1；因此，子囊菌門(mén)是茶園土壤真菌中占有絕對(duì)優(yōu)勢(shì)門(mén)類(lèi)群。不同種植方式茶園土壤真菌群落其他相對(duì)豐度占1%以上的門(mén)包括被孢霉門(mén)（Mortierellomycota，1.6%～5.9%）、毛霉門(mén)（Mucoromycota，1.1%～3.1%）和羅茲菌門(mén)（Rozellomycota，0%～1.3%），這些非優(yōu)勢(shì)門(mén)類(lèi)群可能受到不同茶園種植方式的影響而表現(xiàn)的變化幅度更加劇烈。

不同種植方式茶園土壤真菌群落科水平的組成結(jié)構(gòu)間存在明顯差異（圖2b），其中未分類(lèi)刺盾 炱 菌 目（unclassified_o__Chaetothyriales，11.1%～25.6%）、曲霉菌科（Aspergillaceae，4.5%～21.2%）、 三 型 霉 菌 科（Trimorphomycetaceae，8.3%～ 11.9%）、 黏 毛 菌 科（Myxotrichaceae，3.7%～ 7.4%）、 被 孢 霉 科（Mortierellaceae，1.6%～ 5.9%）、 小 囊 菌 科（Microascaceae，0.9%～ 10.7%）、 叢 赤 殼 科（Nectriaceae，1.8%～5.3%）、未分類(lèi)子囊菌門(mén)（unclassified_p__Ascomycota，1.1%～ 5.4%） 和 肉 座 菌 科（Hypocreaceae，2.1%～4.4%）等為茶園土壤真菌群落的優(yōu)勢(shì)科類(lèi)群。

圖2 不同種植方式茶園和馬尾松林土壤真菌群落門(mén)和科水平上相對(duì)豐度Fig. 2 Relative abundance of fungal communities in soil of different type tea plantations and surrounding P. massoniana forests at phylum and family levels

2.3 不同種植方式茶園間土壤真菌群落相似性

對(duì)土壤真菌群落相似性進(jìn)行層級(jí)聚類(lèi)分析結(jié)果表明（圖3a），茶園與馬尾松林的土壤真菌群落聚為不同的類(lèi)群；而且不同種植方式茶園土壤真菌群落也聚為不同的二級(jí)分支，相對(duì)于新茶園和普通茶園，碾茶園土壤真菌群落可單獨(dú)聚為一個(gè)二級(jí)分支?；谕寥勒婢郝渚嚯x排序的非度量多維尺度NMDS分析結(jié)果（圖3b）具有很好的代表性（Stress＝0.018，小于0.1），進(jìn)一步表明馬尾松林與茶園土壤真菌群落間呈現(xiàn)明顯分異，而且不同種植方式茶園土壤真菌群落也存在較為明顯的分異。

圖3 不同種植方式下茶園及馬尾松林土壤真菌群落的層級(jí)聚類(lèi)與非度量多維尺度分析Fig. 3 Hierarchical clustering and non-metric multidimensional scale analysis on fungal communities in soils at different type tea plantations and surrounding P. massoniana forests

2.4 不同種植方式茶園間土壤真菌群落差異分析

不同樣地間土壤真菌群落相似性ANOSIM檢驗(yàn)分析結(jié)果表明（圖4a），馬尾松林與茶園土壤真菌群落間存在顯著性差異（P＝0.010），但不同種植方式茶園間土壤真菌群落差異未達(dá)顯著程度（P＝0.065）。Mantel檢驗(yàn)結(jié)果表明（圖4b）土壤真菌群落矩陣與環(huán)境因子矩陣間存在顯著性相關(guān)關(guān)系（P＝0.018）；經(jīng)典范關(guān)聯(lián)分析CCA排序結(jié)果表明，有機(jī)質(zhì)和有效磷是黔中山地茶園土壤真菌群落具有重要影響的土壤因子。

圖4 不同種植方式下茶園及馬尾松林土壤真菌群落間相似性差異與典范相關(guān)分析Fig. 4 Similarity difference and canonical correlation analysis on fungal communities in soils at different type tea plantations and surrounding P. massoniana forests

不同種植方式茶園土壤真菌群落差異顯著的屬類(lèi)群間存在明顯差異（圖5）。新茶園與普通茶園間（Tx vs. Tp）土壤真菌群落間有4個(gè)屬表現(xiàn)出顯著性差異（圖5a），分別為綠僵菌屬（Metarhizium）、未 分 類(lèi) 煤 炱 目（unclassified_o__Capnodiales）、樹(shù)枝擔(dān)菌屬（Basidiodendron）和未分類(lèi)壺菌門(mén)（unclassified_p__Chytridiomycota）。新茶園與碾茶園間（Tx vs. Tn）土壤真菌群落有9個(gè)屬呈顯著性差異（圖5b），分別為未分類(lèi)刺盾炱菌目（unclassified_o__Chaetothyriales）、綠僵菌屬（Metarhizium）、頭寬盤(pán)菌屬（Phialocephala）、未分類(lèi)微球黑粉菌綱（unclassified_c__Microbotryomycetes）、膠柄殼屬（Collophora）、未分類(lèi)座囊菌綱（unclassified_c__Dothideomycetes）、土生空?qǐng)F(tuán)菌屬（Cenococcum）、鬼傘屬（Coprinellus）和未分類(lèi)壺菌門(mén)（unclassified_p__Chytridiomycota）。普通茶園與碾茶園間（Tp vs. Tn）土壤真菌群落有6個(gè)屬表現(xiàn)為顯著性差異（圖5c），分別為刺盾炱菌目（unclassified_o__Chaetothyriales）、擬棘殼孢屬（Pyrenochaetopsis）、柱孢屬（Cylindrocarpon）、林地小片耳屬（Holtermanniella）、樹(shù)枝擔(dān)菌屬（Basidiodendron）和須散囊菌屬（Connersia）。(接下一頁(yè)）

圖5 不同種植方式下茶園土壤真菌群落屬水平組間差異Fig. 5 Differences in fungal communities in soils of different type tea plantations at genus level

2.5 土壤真菌群落功能預(yù)測(cè)分析

真菌的生態(tài)功能類(lèi)群是基于共生（Symbiotroph）、腐生（Saprotroph）和病理（Pathotroph）三種營(yíng)養(yǎng)型來(lái)劃分的。通過(guò)FUNGuild數(shù)據(jù)庫(kù)對(duì)土壤真菌功能分類(lèi)進(jìn)行統(tǒng)計(jì)，結(jié)果表明，馬尾松林與茶園土壤真菌功能類(lèi)群組成結(jié)構(gòu)間存在明顯差異（圖6），主要表現(xiàn)為外生菌根真菌（34.4%）和外生菌根－蘭科菌根－根活體營(yíng)養(yǎng)相關(guān)（8.8%）類(lèi)群是馬尾松林土壤真菌群落優(yōu)勢(shì)的功能類(lèi)型，而它們?cè)诓鑸@土壤真菌群落中占比則很小；這一結(jié)果與馬尾松（屬典型外生菌根真菌宿主）和茶樹(shù)（屬叢枝菌根真菌宿主）為不同菌根型植物及其群落的共存植物組成狀況是相對(duì)一致的。茶園土壤真菌優(yōu)勢(shì)功能類(lèi)群主要包括未定義腐生真菌（undefined saprotrophs，21.0%～37.7%）、真菌寄生－未定義腐生生物（Fungal Parasite-Undefined Saprotroph，8.4% ～13.4%）、杜鵑花類(lèi)菌根真菌（Ericoid Mycorrhizal，3.7%～7.4%）、內(nèi)生真菌－凋落物腐生營(yíng)養(yǎng)－土壤腐生營(yíng)養(yǎng)－未定義腐生營(yíng)養(yǎng)（Endophyte-Litter Saprotroph-Soil Saprotroph-Undefined Saprotroph，1.6%～5.9%）和動(dòng)物病原體－內(nèi)生真菌－地衣寄生蟲(chóng)－植物病原體－土壤腐生物－木腐生物（Animal Pathogen-Endophyte-Lichen Parasite-Plant Pathogen-Soil Saprotroph-Wood Saprotroph，1.5%～4.8%）、動(dòng)物病原體－內(nèi)生真菌－植物病原體（Animal Pathogen-Endophyte-Plant Pathogen-Undefined Saprotroph，0.3%～3.9%）等，它們共占茶園土壤真菌功能類(lèi)群相對(duì)豐度的50.3%～58.4%；另外，未知功能真菌類(lèi)群（Unknown）占31.4%～34.6%，這些相關(guān)土壤真菌的生理生態(tài)功能有待進(jìn)一步深入研究。

圖6 不同種植方式茶園與馬尾松林土壤真菌群落功能組成Fig. 6 Composition of soil fungal communities at different type tea plantations and surrounding P. massoniana forests

3 結(jié)論與討論

土壤真菌結(jié)構(gòu)與功能特征在反映茶園土壤環(huán)境質(zhì)量與健康方面提供著重要的依據(jù)。本文探討了黔中山地不同種植方式茶園土壤真菌群落組成、結(jié)構(gòu)與功能特征。結(jié)果表明：子囊菌門(mén)和擔(dān)子菌門(mén)為茶園土壤真菌群落的優(yōu)勢(shì)門(mén)類(lèi)群，二者占土壤真菌群落相對(duì)豐度的88.2%～96.7%（圖2a），盡管對(duì)照馬尾松林土壤真菌群落子囊菌門(mén)和擔(dān)子菌門(mén)也是優(yōu)勢(shì)門(mén)類(lèi)群，且占土壤真菌群落相對(duì)豐度的89.8%，這一比例與新茶園更為接近；但茶園與馬尾松林土壤真菌群落中子囊菌門(mén)與擔(dān)子菌門(mén)間比例（A∶B）上卻存在顯著差異；其中，不同種植方式茶園土壤真菌群落A∶B表現(xiàn)為T(mén)x（3.5∶1）＜Tp（3.8∶1）＜Tn（4.7∶1）；而在馬尾松林A∶B比例約為0.92∶1。這一結(jié)果與云南西雙版納州報(bào)道的森林和現(xiàn)代茶園土壤真菌群落中A∶B比例較為一致［10］；可見(jiàn)，子囊菌門(mén)在茶園土壤真菌群落中占有絕對(duì)優(yōu)勢(shì)；而在外生菌根真菌宿主為主要林分的馬尾松林中子囊菌門(mén)與擔(dān)子菌門(mén)間比例則相對(duì)均衡；因此，針對(duì)不同植被類(lèi)型土壤真菌群落子囊菌門(mén)與擔(dān)子菌門(mén)格局形成機(jī)制有待進(jìn)一步深入研究。

土壤真菌群落不同科類(lèi)群隨茶園的不同種植方式表現(xiàn)出不同的變化方式。本研究中茶園土壤真菌群落中未分類(lèi)刺盾炱菌目（unclassified_o__Chaetothyriales）為科水平相對(duì)豐度最高的類(lèi)群，在碾茶園土壤真菌群落的相對(duì)豐度（25.6%）明顯高于新茶園（11.8%）和普通茶園（11.1%）；另外，本研究中該類(lèi)群共包含了130個(gè)OTUs，目前這些OTUs功能均未得到注釋?zhuān)@也意味著黔中山地土壤真菌該類(lèi)群真菌物種多樣性有待進(jìn)一步深入探索。曲霉菌科（主要為腐生真菌）在碾茶園土壤真菌群落的相對(duì)豐度（21.2%）也明顯高于新茶園（4.5%）和普通茶園（5.7%），這與未分類(lèi)刺盾炱菌目的變化趨勢(shì)相似；而被孢霉科、叢赤殼科和未分類(lèi)子囊菌門(mén)在新茶園和普通茶園土壤真菌群落的相對(duì)豐度（5.9%和5.0%、3.7%和5.3%、3.9%和5.4%）明顯高于碾茶園（1.6%、1.8%和1.1%）；本研究中茶園土壤真菌群落被孢霉科僅包含的被孢霉屬（Mortierella）1個(gè)屬，其一般作為有益菌，可誘發(fā)系統(tǒng)性抵抗和促進(jìn)植物生長(zhǎng)，加強(qiáng)對(duì)植物病原體抵抗；而赤殼科主要包括了世界廣泛地分布于土壤和有機(jī)體內(nèi)的鐮刀菌屬（Fusarium）和其它少量的病原真菌和腐生真菌，如麗赤殼屬（Calonectria）、柱孢屬（Cylindrocarpon）和泥赤叢殼菌屬（Ilyonectria）等；這預(yù)示著變化趨勢(shì)相似的土壤真菌類(lèi)群，但它們?cè)谏鷳B(tài)功能并不趨同，不同種植方式茶園土壤真菌群落結(jié)構(gòu)上的變化可能是真菌群落構(gòu)建機(jī)制中不同物種間相互作用的動(dòng)態(tài)權(quán)衡。

茶園土壤真菌群落的結(jié)構(gòu)與功能受到茶樹(shù)品種、植茶年限、施肥管理、種植管理、土壤環(huán)境等方面的影響［10，16，22，23］。本研究中調(diào)查的不同種植方式茶園土壤真菌群落在組成、結(jié)構(gòu)和功能上存在著明顯的差異（圖3～6）。其中，不同種植方式茶園土壤真菌群落OTUs豐富度表現(xiàn)出（表2）：Tx（1636）＞Tp（1281）＞Tn（1197），這表明新茶園隨種植年限和管理強(qiáng)度對(duì)土壤真菌多樣性有下降趨勢(shì)影響，這可能與茶園土壤酸化、茶樹(shù)修剪物的降解產(chǎn)物和茶樹(shù)本身固有的酚類(lèi)物質(zhì)對(duì)茶園的真菌群落有一定的抑制作用有關(guān)［23，24］。目前有關(guān)茶園土壤真菌群落多樣性隨種植茶樹(shù)年限變化的研究仍較少、且結(jié)果不盡一致［23，24］，這可能由于受到選擇茶園地理氣候條件、種植年限梯度、管理方式及樣品檢測(cè)方法等影響；可見(jiàn)，針對(duì)茶園土壤真菌群落演變機(jī)制的研究有待加強(qiáng)。

層級(jí)聚類(lèi)分析和NMDS分析結(jié)果表明，不同種植方式茶園土壤真菌群落存在明顯的分異（圖3），碾茶園與新茶園和普通茶園土壤真菌群落間的相似性較小?？梢?jiàn)，碾茶園由于茶樹(shù)高度控制、遮陰處理和機(jī)械收割等種植管理差異因素，可對(duì)土壤真菌群落產(chǎn)生了更加深刻的影響。茶園土壤真菌群落通過(guò)FUNGuild功能分類(lèi)結(jié)果表明，馬尾松林與茶園土壤真菌功能類(lèi)群組成結(jié)構(gòu)存在明顯差異（圖6），這可能是由于不同植被類(lèi)型會(huì)通過(guò)植物－真菌－土壤間的協(xié)同演變與生境的環(huán)境質(zhì)量改變，進(jìn)而直接和間接地影響土壤真菌群落多樣性與功能特征［25，26］。不同種植方式茶園土壤真菌功能類(lèi)群組成結(jié)構(gòu)也存在一定差異，這可能主要由于不同人為種植干擾方式的差異影響著茶園的土壤環(huán)境和微生境，以及茶樹(shù)的生長(zhǎng)狀況，進(jìn)而間接的影響著土壤真菌群落多樣性及功能特征。本研究茶園生境中杜鵑花類(lèi)菌根真菌占有一定比例，針對(duì)該類(lèi)菌根真菌特征也有待進(jìn)一步探索。另外，土壤的非生物因子往往對(duì)土壤真菌群落產(chǎn)生重要影響，如pH、C/N和SOC等［25，26］。本研究通過(guò)CCA分析結(jié)果表明，有機(jī)質(zhì)和有效磷是茶園土壤真菌群落重要的土壤環(huán)境影響因子（圖4），該結(jié)果與已報(bào)道的茶園土壤真菌重要影響的土壤因子相對(duì)一致［10，24］；這也預(yù)示著在茶園種植管理過(guò)程中應(yīng)注意茶園土壤C和P養(yǎng)分的科學(xué)管理，可作為對(duì)茶園土壤真菌群落的一種潛在的生態(tài)調(diào)控方式。