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

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

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

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

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

1 材料與方法

1.1 研究區(qū)域概況

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

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

選取了新茶園（Tx）、普通茶園（Tp）和碾茶園（Tn）3種不同種植管理方式的茶園樣地，以周邊馬尾松林（Pinus massoniana，Mp）作為對照（表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ì)后，全部過2 mm土篩，均勻混合后分為2份，一份經(jīng)室內(nèi)自然風(fēng)干后用于測定土壤理化性質(zhì)；另一份經(jīng)冷凍干燥保存于-80℃冰箱，用于高通量測序的土壤總DNA提取。

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

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

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

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

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

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

經(jīng)高通量測序得到的雙端序列原始數(shù)據(jù)，根據(jù)PE reads之間的overlap關(guān)系，將成對的reads拼接（merge）成一條序列，利用FLASH軟件對reads質(zhì)量和拼接效果進(jìn)行質(zhì)控過濾，再根據(jù)序列首尾兩端的barcode和引物序列區(qū)分樣品得到有效序列。通過Uparse軟件在相似水平為97%的OTU進(jìn)行聚類分析；樣品序列按最小序列抽平，本研究中每個(gè)樣品按60540條序列進(jìn)行抽平，通過QIIME軟件獲取OTU豐度、α多樣性（包括Sobs、Shannon、Simpson、ACE和Chao等多樣性指數(shù)）、稀釋曲線和韋恩圖；采用Welch T法檢驗(yàn)不同茶園種植方式土壤真菌群落在屬水平上組間的顯著性差異，采用“fdr”法對P值進(jìn)行多重檢驗(yàn)校正，采用“Welch’s intverted”法計(jì)算95%的置信區(qū)間；基于非加權(quán)unifrac距離矩陣按照average算法進(jìn)行層級聚類；采用加權(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)方式對土壤真菌功能進(jìn)行預(yù)測分析［21］。

2 結(jié)果與分析

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

土壤真菌ITS擴(kuò)增子高通量測序，共得到優(yōu)化有效序列總數(shù)為574809，物種稀釋曲線呈現(xiàn)出隨著測得序列數(shù)量的增加逐漸趨于平坦（圖1a），這表明測試有效序列數(shù)據(jù)信息量能夠反映土壤樣本真菌群落多樣性。按＞97%相似度進(jìn)行OTU注釋，土壤真菌群落共鑒定得到13門、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 不同種植方式下茶園和馬尾松林土壤真菌的稀釋曲線（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 不同種植方式下茶園和馬尾松林土壤真菌群落的分類鑒定與多樣性指數(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)特征

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

不同種植方式茶園土壤真菌群落科水平的組成結(jié)構(gòu)間存在明顯差異（圖2b），其中未分類刺盾 炱 菌 目（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%）、未分類子囊菌門（unclassified_p__Ascomycota，1.1%～ 5.4%） 和 肉 座 菌 科（Hypocreaceae，2.1%～4.4%）等為茶園土壤真菌群落的優(yōu)勢科類群。

圖2 不同種植方式茶園和馬尾松林土壤真菌群落門和科水平上相對豐度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 不同種植方式茶園間土壤真菌群落相似性

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

圖3 不同種植方式下茶園及馬尾松林土壤真菌群落的層級聚類與非度量多維尺度分析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

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

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

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

真菌的生態(tài)功能類群是基于共生（Symbiotroph）、腐生（Saprotroph）和病理（Pathotroph）三種營養(yǎng)型來劃分的。通過FUNGuild數(shù)據(jù)庫對土壤真菌功能分類進(jìn)行統(tǒng)計(jì)，結(jié)果表明，馬尾松林與茶園土壤真菌功能類群組成結(jié)構(gòu)間存在明顯差異（圖6），主要表現(xiàn)為外生菌根真菌（34.4%）和外生菌根－蘭科菌根－根活體營養(yǎng)相關(guān)（8.8%）類群是馬尾松林土壤真菌群落優(yōu)勢的功能類型，而它們在茶園土壤真菌群落中占比則很小；這一結(jié)果與馬尾松（屬典型外生菌根真菌宿主）和茶樹（屬叢枝菌根真菌宿主）為不同菌根型植物及其群落的共存植物組成狀況是相對一致的。茶園土壤真菌優(yōu)勢功能類群主要包括未定義腐生真菌（undefined saprotrophs，21.0%～37.7%）、真菌寄生－未定義腐生生物（Fungal Parasite-Undefined Saprotroph，8.4% ～13.4%）、杜鵑花類菌根真菌（Ericoid Mycorrhizal，3.7%～7.4%）、內(nèi)生真菌－凋落物腐生營養(yǎng)－土壤腐生營養(yǎng)－未定義腐生營養(yǎng)（Endophyte-Litter Saprotroph-Soil Saprotroph-Undefined Saprotroph，1.6%～5.9%）和動(dòng)物病原體－內(nèi)生真菌－地衣寄生蟲－植物病原體－土壤腐生物－木腐生物（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%）等，它們共占茶園土壤真菌功能類群相對豐度的50.3%～58.4%；另外，未知功能真菌類群（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é)果表明：子囊菌門和擔(dān)子菌門為茶園土壤真菌群落的優(yōu)勢門類群，二者占土壤真菌群落相對豐度的88.2%～96.7%（圖2a），盡管對照馬尾松林土壤真菌群落子囊菌門和擔(dān)子菌門也是優(yōu)勢門類群，且占土壤真菌群落相對豐度的89.8%，這一比例與新茶園更為接近；但茶園與馬尾松林土壤真菌群落中子囊菌門與擔(dān)子菌門間比例（A∶B）上卻存在顯著差異；其中，不同種植方式茶園土壤真菌群落A∶B表現(xiàn)為Tx（3.5∶1）＜Tp（3.8∶1）＜Tn（4.7∶1）；而在馬尾松林A∶B比例約為0.92∶1。這一結(jié)果與云南西雙版納州報(bào)道的森林和現(xiàn)代茶園土壤真菌群落中A∶B比例較為一致［10］；可見，子囊菌門在茶園土壤真菌群落中占有絕對優(yōu)勢；而在外生菌根真菌宿主為主要林分的馬尾松林中子囊菌門與擔(dān)子菌門間比例則相對均衡；因此，針對不同植被類型土壤真菌群落子囊菌門與擔(dān)子菌門格局形成機(jī)制有待進(jìn)一步深入研究。

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

茶園土壤真菌群落的結(jié)構(gòu)與功能受到茶樹品種、植茶年限、施肥管理、種植管理、土壤環(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)度對土壤真菌多樣性有下降趨勢影響，這可能與茶園土壤酸化、茶樹修剪物的降解產(chǎn)物和茶樹本身固有的酚類物質(zhì)對茶園的真菌群落有一定的抑制作用有關(guān)［23，24］。目前有關(guān)茶園土壤真菌群落多樣性隨種植茶樹年限變化的研究仍較少、且結(jié)果不盡一致［23，24］，這可能由于受到選擇茶園地理氣候條件、種植年限梯度、管理方式及樣品檢測方法等影響；可見，針對茶園土壤真菌群落演變機(jī)制的研究有待加強(qiáng)。

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