梵凈山典型植被類型的根系真菌群落結構及功能類群特征研究

劉敏　李龍　李海波等

關鍵詞：植被類型；FUNGuild;根系真菌；功能類群；細根性狀

中圖分類號：S718.81 文獻標志碼：A 文章編號：1001-1498（2023）01-0047-12

梵凈山地處熱帶、亞熱帶生物區系向溫帶生物區系的生態交錯區，由于獨特的地理位置、復雜的生物和生態演化過程、多樣化的植物類群特征、眾多地方特有和瀕危動植物種類等特點，引起國內外不同領域學者的廣泛關注，相關研究涉及到土壤理化性質、物種種類及多樣性、珍稀瀕危物種種群結構特征等多個方面。關于地下微生物群落的研究則集中在叢枝菌根真菌，對真菌群落結構及功能類群特征尚無系統報道。真菌特別是某些功能類群在促進有機質分解和養分循環、增強植物對礦物質營養元素的吸收、提高植物對生物和非生物脅迫的耐受性、維持生態系統的穩定性和生產力等方面發揮著重要作用。此外，真菌群落組成不僅受土壤理化性質、氣候等非生物因素影響，其對植物種類、群落組成和多樣性等生物因素也較為敏感，而植被類型可通過影響非生物和生物因素作用于真菌群落。尤為特殊的是，不同植被類型還會塑造特定的真菌功能類群：Tedersoo等通過對全世界范圍內不同植被類型的土壤真菌群落進行分析后發現，外生菌根真菌在溫帶落葉林中較為豐富，熱帶潮濕雨林和熱帶干旱森林中植物病原菌更為多樣；Liu等發現在中國東部森林地區，溫帶森林以外生菌根真菌為主，熱帶和亞熱帶森林中具有豐富的叢枝菌根真菌和植物病原菌；Sheng等研究了神農架山地森林生態系統中的土壤真菌功能類群，發現外生菌根真菌在落葉闊葉林、針闊混交林和針葉林中尤為豐富，植物病原菌一木腐生菌主要出現在常綠闊葉林，植物病原菌和寄生性真菌在山頂灌叢中具有較高相對豐度。菌在山頂灌叢中具有較高相對豐度。

梵凈山擁有的典型植物區系過渡性及地球上同緯度保存完好的森林生態系統，為豐富且獨特真菌種類的形成創造良好生存條件。此外，由于水熱條件等差異造就了梵凈山多樣且分異明顯的植被類型，為探討植被類型與真菌群落間的關系提供了天然理想實驗室。本研究以梵凈山國家級自然保護區典型且具有代表性的常綠闊葉林、常綠落葉闊葉混交林及位于山頂附近的亞高山矮林3種森林植被類型為研究對象，分析真菌群落結構和功能類群特征與植被類型間的關系，探討土壤性質、細根性狀和葉性狀對真菌群落的影響。研究結果可為后續探討亞熱帶地區不同植被類型樹種的環境適應機制及地上和地下部分的養分交換策略提供來自于地下真菌群落方面的數據支撐。

1材料與方法

1.1研究區概況及樣本采集

研究區域位于梵凈山國家級自然保護區（27°46′50″N～28°1′30″N，108°35′55″E～108°48′30″E）內，地處貴州銅仁市江口縣、印江自治縣、松桃自治縣交界處，面積約為419 km²。亞熱帶季風濕潤氣候，年平均氣溫16℃，年平均降雨量1300 mm，平均相對濕度80%以上，土壤以山地黃壤和黃棕壤為主。植被類型豐富，具有代表性和面積較大的是常綠闊葉林和常綠落葉闊葉混交林，前者主要植物種類為栲（Castanopsisfargesii Franch.）、米櫧（C.carlesii（Hemsl.）Hayata.）、巴東櫟（Quercus engleriana Seem.）、山茶（Camellia japonica L.）等，后者植物種類包括水青岡（Fagus longipetiolata Seem.）、貴州青岡（Cyclobalanopsis argyrotricha （A.Camus）Chun et Y.T.Chang ex Y.C.Hsu et H.W. Jen）和扇葉槭（Acer flabellatum Rehd.）等。亞高山矮林是位于山頂附近的特殊植被類型，以吊鐘花（Enkianthus quinqueflorus Lour.）、耳葉杜鵑（Rhododendron auriculatum Hemsl.）和長萼杜鵑（Rhododendron longicalyx Fang f.）等杜鵑花科植物為主。

在常綠闊葉林海拔600m處、常綠落葉闊葉混交林海拔1500 m處和亞高山矮林海拔2200 m處分別垂直于海拔高度設置1條100 m的樣方線，每隔10 m設置一個采樣點，中心象限法采集樣本，即在每個采樣點設置4個方向，選擇距離采樣點最近的2個樹種（起測胸徑為2cm），去除表層石塊和落葉等雜物，從樹干開始挖掘完整根系，抖落法收集根際土壤。由于部分采樣點樹木生長在坡度較大區域或巖石縫隙中，采集根系和土壤極為困難，僅48個樣本用于后續分析（每種植被類型16個）。其中，常綠闊葉林中，栲5個，巴東櫟4個，山茶3個，黃丹木姜子（Litsea elongata （Wall.ex Nees） Benth.et Hook.f.）2個，四照花（Cornus kousa subsp.chinensis（Osborn）Q.Y.Xiang）1個，柃木（Eurya japonica Thunb.）1個；常綠落葉闊葉混交林中，貴州青岡5個，水青岡4個，多脈青岡（Cyclobalanopsis multinervis W.C.Cheng＆T.Hong）3個，云山青岡（C.sessilifolia （Blume） Schottky）3個，青岡（C.glauca （Thunberg） Oersted）1個；亞高山矮林中，耳葉杜鵑8個，吊鐘花6個，長萼杜鵑2個。同時，在每個植株上隨機摘取20片向陽、完全展開、沒有病蟲害的葉片，用于葉性狀指標分析。根系樣本裝在放有干冰的泡沫箱運回實驗室，一部分保存在4℃冰箱用于分析根性狀指標，另一部分保存在-70℃冰箱用于真菌測定。

1.2細根性狀、葉性狀和土壤理化性質分析

根系用自來水清洗掉附著土壤后，依據Pregitzer等根序分級法，按照序級進行分級，前兩個序級歸為吸收根并進行分離。EPSON Perfection V700（EPSON America Inc.）掃描細根樣本，根系圖像分析軟件WinRhizo（Regent Instruments Inc.，Quebec，Canada）分析掃描的圖片獲得平均根直徑/mm、根長/cm、根表面積／cm²和根體積／cm³。將掃描過的細根稱量質量后，放入信封袋在70℃下烘48 h，再次稱量質量得到根生物量／g。計算比根長（根長／生物量），（cm·g^-1）、根比表面積（根表面積／生物量），（cm²cm·g^-1）、根組織密度（生物量，根體積），（g·g^-3）和根干物質含量（根干質量，根鮮質量）／（g·g^-1）。

葉片放平整后利用Nikon Coolpi×5000數碼相機（Nikon Corp.，Tokyo，Japan）拍照，ImageJ（Systat Software Inc.，Richmond，CA）軟件計算葉面積／cm²。稱重后，在70℃下烘72 h，再次稱量質量得到葉生物量，g。計算比葉面積（葉面積／葉生物量），（cm²·g^-1）和葉干物質含量（葉干質量／葉鮮質量）／（g·g^-1）。

所有土壤、細根和葉樣本均重復測定3次。土壤pH測定采用pH計法（土水比例1：2.5）；土壤、細根和葉碳、氮含量測定利用CN元素分析儀（Vario EL Ⅲ，Hanau.Germany）;土壤全磷測定采用HClO₄-H₂SO₄消煮一鉬銻抗比色法，細根和葉磷含量測定采用H2S04-H202消煮一鉬銻抗比色法。

1.3根系DNA提取、PCR擴增和測序

根系用蒸餾水沖洗10 min后，70%乙醇浸泡1min，3.3%次氯酸鈉浸泡3min，75%乙醇浸泡30 s，無菌水反復沖洗后驗證消毒效果。滅菌后根系用液氮冷凍研磨，DNeasy Plant Mini Kit（Qiagen Gmbh， Hilden， Germany）提取DNA，0.8%瓊脂糖凝膠電泳檢測DNA提取質量，紫外分光光度計（Thermo Fisher Scientific，MA，USA）準確定量DNA濃度。真菌rDNA-ITS保守區引物ITSIF12（5-GTGARTCATCGAATCTTTG-3）和ITS2 （5-TCCTCCGCTTATTGATATGC-3）進行PCR擴增。擴增體系包括12.5μL Phusion（R）Hot Start Flex 2X Master Mix（New England Biolabs，Ipswich，USA）、10 μmoKl·L^-1正向引物和反向引物各2.5 pL、Genomic DNA 50 ng，補充滅菌超純水到25 pL。擴增條件：98℃預變性30 s；98℃變性10 s，54℃退火30 s、72℃延伸45 s，共35個循環；72℃延伸10 min。擴增產物利用2%瓊脂糖凝膠電泳檢測，AxyPrepPCR Cleanup Kit（Axygen Biosciences， CA， USA）純化，Quant-iTTM PicoGreen^TMdsDNA Assay Kit（Thermo Fisher Scientific，MA，USA）定量濃度，Illumina MiSeq平臺測序。

1.4生物信息學分析和真菌功能預測

DADA2（divisive amplicon denoisingalgorithm）進行序列拼接、去除錯誤堿基序列、去掉barcode和引物序列、過濾嵌合體后，Ribosomal Database Project （RDP）Naive Bayesian Classifier得到擴增序列變異（amplicon sequence variants，ASVs）。ASVs相當于以100%相似度聚類OTUs（operational taxonomic units），但與OTUs比，具有較高數據精確度與物種分辨率。在UNITE數據庫中比對，得到序列注釋信息，FUNGuildc比對和解析真菌營養型和功能類群，置信度為“極可能”（highly probable）和“很可能”（probable）等級的功能分組被用于后續分析。

1.5數據統計分析

混合效應模型探討植被類型對土壤理化性質、細根性狀、葉性狀、和真菌α多樣性指數的影響，植被類型為固定因子，樹種為隨機因子，并檢驗模型殘差的正態性、齊性和獨立性，若不滿足則對數據進行轉換。LEfSe（http：//huttenhower.sph.harv-ard.edu/lefse/）判別不同植被類型間真菌的顯著差異類群。非度量多維標度（nonmetric multidimensional scaling，NMDS）排序和置換多元方差分析（permutational multivariate analysis of variance，PERMANOVA）對不同植被類型間真菌群落組成進行比較，Kruskal-Wallis test和Wilcoxon signed-rank tests判別植被類型間真菌Bray-Curtis相異性指數的高低。植被類型、樹種、細根性狀、葉性狀和土壤性質對真菌α多樣性指數影響的相對重要性基于ABT（aggregated boosted tree）分析，并計算其與真菌α多樣性指數的Pearson相關性。降趨勢對應分析（detrended correspondence analysis，DCA）檢驗真菌數據為單峰分布還是線性分布，DCA排序軸梯度長度最大值為2.657，小于3.0，選擇線性排序模型中的冗余分析（redundancy analysis，RDA）探討影響真菌群落的因子；執行RDA前，對所涉及解釋變量進行前向選擇。方差分解分析評估植被類型、土壤性質、細根性狀和葉性狀對真菌群落影響的大小。除特別說明外，統計分析均在R中完成。

2結果

2.1不同植被類型根系真菌群落特征概況

48個根系樣本共得到5 103個ASVs，屬于9門、193科和385個屬（不包括unidentified fungi）。子囊菌門（Ascomycota）和擔子菌門（Basidiomycota）為優勢門，占全部ASVs數的89.93%。小蔓毛殼科（Herpotrichiellaceae）、球囊霉科（Glomeraceae）和柔膜菌科（Helotiaceae）為優勢科，晶杯菌科（Hyaloscyphaceae）、Sebacinales_Group_A、被孢霉科（Mortierellaceae）、皮盤菌科（Dermateaceae）、紅菇科（Russulaceae）、小菇科（Mycenaceae）、叢赤殼科（Nectriaceae）、水盤菌科（Vibrisseaceae）、錘舌菌科（Leotiaceae）、肉座菌科（Hypocreaceae）和革菌科（Thelephoraceae）相對豐度超過1%。常綠闊葉林具有最多ASVs數（2287個），子囊菌門為優勢門，球囊霉科為優勢科；亞高山矮林共得到ASVs數為1636個，子囊菌門為優勢門，柔膜菌科和晶杯菌科為優勢科；常綠落葉闊葉混交林ASVs數最少，為1608個，子囊菌門和擔子菌門為優勢門，球囊霉科和小蔓毛殼科為優勢科。LEfSe分析（圖1A）顯示，常綠闊葉林顯著差異類群為球囊菌門（Glomeromycota）、球囊菌綱（Glomeromycetes）、糞殼菌綱（ Sordariomycetes）、Pezizomycotina_Incertae_sedis綱、Pezizomycotina_lncertae_sedis目、糞殼菌目（Sordariales）、雞油菌目（Cantharellales）、球囊霉目（Glomerales）、球囊霉科、叢赤殼科、Pezizomycotina_lncertae_sedis科和發菌科（Trichocomaceae），常綠落葉闊葉混交林中顯著差異類群為革菌目（Thelephorales）、紅菇目（Russulales）、革菌科和紅菇科，錘舌菌綱（Leotiomycetes）、柔膜菌目（Helotiales）、蠟殼耳目（Sebacinales）、Sebacinales_Group_B科、柔膜菌科和錘舌菌科是亞高山矮林的顯著差異類群。

2.2不同植被類型根系真菌功能類群分析

3種植被類型樹種根系真菌在“極可能”和“很可能”水平上進行功能分組的ASVs數為2231個，共被分為病理營養型（pathotroph）、共生營養型（symbiotroph）、腐生營養型（saprotroph）、病理-共生營養型（pathotroph-symbiotroph）和病理-腐生營養型（pathotroph-saprotroph）5個營養型。其中，41.86%的ASVs被歸為腐生營養型，其次是共生營養型（28.73%）和病理營養型（15.37%），其余2個營養型所占比例較低。在功能亞類水平上，undefined saprotroph相對豐度最高，其次是叢枝菌根（arbuscular mycorrhizal）真菌、外生菌根（ectomycorrhizal）真菌、植物病原菌（plant pathogen）、內生真菌（endophyte）、寄生性真菌（fungal parasite）和杜鵑花類菌根（ericoid mycorrhizal）真菌，其余功能亞類相對豐度較低。此外，叢枝菌根真菌僅出現在常綠闊葉林和常綠落葉闊葉混交林，且在常綠闊葉林（7.77%）的相對豐度高于常綠落葉闊葉混交林（7.30%）；外生菌根真菌在常綠落葉闊葉混交林（10.20%）的相對豐度明顯高于常綠闊葉林（3.28%）和亞高山矮林（0.98%），內生真菌和杜鵑花類菌根真菌只出現在亞高山矮林。LEfSe分析（圖1B）還發現，叢枝菌根真菌是常綠闊葉林的顯著差異功能類群，外生菌根真菌是常綠落葉闊葉混交林的顯著差異功能類群，亞高山矮林的顯著差異功能類群為杜鵑花類菌根真菌。

2.3不同植被類型根系真菌群落結構分析

混合效應模型結果表明，植被類型能顯著影響真菌豐富度和Shannon指數（P<0.05），對Simpson指數和Pielou指數無顯著影響（P>0.05），常綠落葉林中真菌豐富度顯著高于常綠落葉闊葉混交林和亞高山矮林，Shannon指數顯著高于亞高山矮林（圖2）。

NMDS分析結果顯示真菌群落組成在植被類型間存在差異（圖3A），PERMANOVA揭示了相同的結果（F=2.856，P=0.001）。此外，Bray-Curtis相異性指數（圖3B）在3種植被類型間存在極顯著差異（P<0.001），且在常綠落葉闊葉混交林中最高，常綠闊葉林次之，亞高山矮林最低，表明真菌群落組成在亞高山矮林樣本間具有最高相似性，其次是常綠闊葉林，常綠落葉闊葉混交林的相似性最低，NMDS分析也印證了該結果。

2.4植被類型、樹種、土壤性質、細根性狀和葉性狀與根系真菌群落的關系分析

以植被類型為固定因子，樹種為隨機因子構建的混合效應模型發現，植被類型可顯著影響土壤pH、土壤全碳、土壤全氮、土壤全磷、比根長、根組織密度、平均根直徑、根比表面積、根碳含量、根氮含量、葉面積和葉碳含量（P<0.05）。此外，常綠闊葉林和常綠落葉闊葉混交林土壤pH顯著高于亞高山矮林，亞高山矮林土壤具有最高的全碳和全氮含量，常綠闊葉林具有最高的土壤全磷含量；比根長、根比表面積和根碳含量在亞高山矮林顯著高于常綠闊葉林和常綠落葉闊葉混交林，常綠闊葉林具有最高的根組織密度、平均根直徑和根氮含量；葉面積在常綠闊葉林和常綠落葉闊葉混交林顯著高于亞高山矮林，葉碳含量在亞高山矮林顯著高于常綠闊葉林和常綠落葉闊葉混交林（表1）。

ABT分析用來揭示植被類型、樹種、土壤理化性質、細根性狀和葉性狀對真菌豐富度和多樣性指數影響的相對重要性（圖4）。結果表明：根碳含量是影響真菌豐富度的首要因子，相對重要性可達到24.47%（圖4A）；影響Shannon指數、Simpson指數和Pielou指數的首要因子均為葉干物質含量，相對重要性分別為14.63%、18.08%和15.55%（圖4B、圖4C、圖4D）。Pearson相關性分析結果指出，根碳含量與真菌豐富度（r=-0.626，P<0.001）、Shannon指數（r=-0.407，P=0.004）呈顯著負相關，葉干物質含量與Shannon指數（r=-0.293，P=0.043）、Simpson指數（r=-0.591，P=0.002）、Pielou指數（r=-0.531，P=0.0259）呈顯著負相關。

對包括植被類型和樹種在內的22個解釋變量進行前向選擇后剩余7個，分別是植被類型、土壤pH、土壤全磷、根干物質含量、根碳含量、葉干物質含量和葉碳含量，與真菌群落構建RDA模型（圖5A）。結果發現，植被類型、土壤pH、根碳含量和葉碳含量是影響真菌群落的顯著因子（P<0.05）。方差分解分析進一步揭示植被類型、土壤性質（土壤pH、土壤全磷）、細根性狀（根干物質含量、根碳含量）和葉性狀（葉干物質含量、葉碳含量）對真菌群落的影響（圖5B），發現植被類型（26.96%）是影響真菌群落的首要因子，其次是細根性狀（20.10%）和土壤性質（13.53%），葉性狀（3.92%）影響最小。

3討論

3.1不同植被類型根系真菌和功能類群

本研究基于高通量測序技術探討世界自然遺產地梵凈山地區3種植被類型根系真菌群落，發現子囊菌門和擔子菌門為優勢真菌門，森林生態系統中枯枝落葉和枯死木是植物和微生物主要養分來源，子囊菌門和擔子菌門的大多數真菌能促進纖維素和木質素等的分解和礦化從而調控生物地球化學循環。此外，叢枝菌根真菌為常綠闊葉林的顯著差異功能類群，外生菌根真菌是常綠落葉闊葉混交林的顯著差異功能類群，可能因為常綠闊葉林養分周轉速率一般高于落葉闊葉林，不同菌根真菌采取的養分循環策略存在差異，叢枝菌根真菌通常在養分周轉速率快的地區占優勢，外生菌根真菌普遍存在于養分循環慢的區域。同時，本研究還發現，杜鵑花類菌根真菌主要出現在亞高山矮林，一方面梵凈山地區亞高山矮林的優勢植物主要以杜鵑花科為主，這為杜鵑花類菌根真菌的生存提供大量可供選擇的宿主；另一方面亞高山矮林生長在山頂土層淺薄、土壤貧瘠、晝夜溫差大、溫度低的惡劣生境中，豐富和多樣化的杜鵑花類菌根真菌可通過促進杜鵑花科植物對礦物質營養元素的吸收及產生具有抗氧化和水解作用的酶來幫助宿主植物獲得復合有機物中的養分，從而增強杜鵑花科植物對山頂惡劣環境的適應能力。

3.2不同植被類型真菌群落結構的差異

常綠闊葉林中植物根系真菌具有最高的α多樣性，可能原因如下：第一，常綠闊葉林分布區域的溫度高于常綠落葉闊葉混交林和亞高山矮林，而較高的溫度能夠增強初級生產力及枯枝落葉物等的分解速率，為異養真菌生存提供更多資源；第二，亞熱帶森林中大部分樹種為菌根依賴型樹種，菌根真菌的存在能幫助宿主植物獲取更多可利用土壤養分。雖然常綠闊葉林土壤磷含量高于其它兩種植被類型，但在整個亞熱帶森林，磷是植物生長的重要限制元素，再加上常綠闊葉林具有最低的土壤氮含量。一方面某些菌根真菌的存在有助于溶解土壤中不能被植物吸收利用的磷酸鹽和固定生物氮，另一方面宿主植物會傾向于選擇更多的菌根真菌搭檔來獲得生長所必須的氮和磷；使得常綠闊葉林具有豐富且多樣的真菌種類。此外，真菌群落組成在3種植被類型間存在差異，RDA結果也發現植被類型是影響真菌群落組成的顯著因子。同在中國亞熱帶地區，Han等的研究表明真菌群落組成在常綠闊葉林、落葉闊葉林、針闊混交林、針葉林和次生灌木林存在顯著差異。植被類型間不僅植物種類及群落組成、植物功能性狀等生物因素存在差異，凋落物厚度及成分、土壤理化性質等非生物因素也有較大差別，從而可能直接或間接影響真菌群落組成。亞高山矮林真菌具有最低的Bray-Curtis相異性指數，是由于亞高山矮林所處區域植物種類（以杜鵑花科植物為主）及土壤、氣候條件的獨特性而孕育出高度相似的真菌群落。

3.3驅動根系真菌群落組成的主要影響因素

植物功能性狀是植物長期適應外部環境變化的結果，近年來越來越多的研究關注根和葉功能性狀與真菌群落的關系：Prada-Salcedo等指出根氮含量和根組織密度與真菌群落組成密切相關；Sweeney等發現根直徑、根氮含量和比根長在驅動真菌群落組成中發揮著重要作用；而本研究發現根碳含量和葉碳含量能顯著影響真菌群落組成，出現這些不同結果可能與宿主植物種類及采樣部位的差異有關。當土壤養分受到限制時，宿主植物通常會將自身大量碳交易給菌根共生體，讓更多菌根真菌搭檔去競爭碳水化合物從而提高養分吸收效益，可能是根碳含量和葉碳含量顯著影響真菌群落組成且根碳含量與真菌豐富度和Shannon指數呈顯著負相關的原因。本研究還發現土壤pH是影響真菌群落組成的顯著因子，先前同在森林生態系統中的研究也指出，土壤pH可以影響腐生真菌、叢枝菌根真菌以及外生菌根真菌，這是因為土壤pH會改變養分可利用性。此外，也有研究進一步探討了土壤性質、根性狀和葉性狀對真菌群落的影響大小：Sweeney等指出與葉性狀相比，根性狀能決定叢枝菌根真菌和腐生真菌的豐富度和多樣性；Dong等發現根系真菌多樣性和群落組成主要受土壤性質（土壤pH、土壤碳磷比和氮磷比）的顯著影響，而與葉性狀無顯著相關性。在本研究中，方差分解分析結果表明，細根性狀對真菌群落的影響大于土壤性質和葉性狀，可能原因：一方面植物細根（吸收根）是大多數真菌的棲息場所，較高比根長和根比表面積能為其生存提供更多空間；另一方面細根是某些真菌與宿主植物進行物質交換區域，植物可能會更加依賴于共生真菌來應對亞熱帶地區缺磷的土壤環境，使某些真菌功能類群（如叢枝菌根真菌、外生菌根真菌和杜鵑花類菌根真菌等）對細根性狀變化較為敏感。對上述原因分析只是基于當前生態學知識的推測，后續研究會從植物種類與真菌群落結構的關系、不同菌根類型與宿主植物的養分獲取途徑等方面深入解析梵凈山地區地上植物群落和地下真菌群落的關聯機制。

4結論

世界自然遺產地梵凈山3種典型植被類型的植物根系中，子囊菌門和擔子菌門為優勢真菌門，小蔓毛殼科、球囊霉科和柔膜菌科為優勢真菌科，腐生營養型真菌所占比例較高。叢枝菌根真菌是常綠闊葉林的顯著差異功能類群，外生菌根真菌是常綠落葉闊葉混交林的顯著差異功能類群，而亞高山矮林的顯著差異功能類群為杜鵑花類菌根真菌。植被類型、土壤pH、根碳含量和葉碳含量顯著影響根系真菌群落組成，且細根性狀對真菌群落的影響大于土壤性質和葉性狀。研究結果不僅為豐富世界自然遺產地的物種多樣性提供了數據支撐，還能夠為深入探討亞熱帶地區土壤·真菌·植物養分交易策略奠定基礎。