基于高通量測序對不同茶園土壤真菌群落結構特征的分析

中圖分類號：S182 文獻標志碼：A 文章編號：1002-1302（2025）11-0220-13

土壤是茶樹生長發育的基質，也是茶樹根系獲取養分、抵抗疾病、健康成長的關鍵所在。它是一個復雜的動態生態系統，據估計植物根部周圍的根際土壤中每克含有數十億微生物，種類數以方計[]。而土壤微生物在生態系統功能中發揮著重要作用，尤其是對土壤中的某些過程很重要，如礦質養分循環、有機質周轉和土壤結構等[2]。同時，微生物群落結構組成及土壤理化因子直接或間接影響植物產量和品質[3]。湖北省恩施州是中國重要的優質茶葉產區，種植規模位居湖北省第一，全國地級市產茶區第四。截至2023年年底，恩施州茶園面積為 1 200km² ，茶產業是恩施州的第一特色產業和優勢富民民產業。因此，分析恩施州茶園土攘微生物群落組成和土壤養分指標之間的相關性對培育高效優質茶樹樹體和明晰的茶樹地下生態模式，形成茶樹健康生長地下土壤模式服務具有十分重要的指導意義。土壤微生物是土壤中最活躍且多樣性最為豐富的組分之一，但受外界因素影響也大。萬人源調查了云南茶園土壤真菌群落結構在門分類水平上，以子囊菌、擔子菌、接合菌、壺菌和球囊菌5個門為主[4]；王峰研究指出有機管理方式對茶園土壤真菌群落結構能產生顯著影響[5]；王華等研究指出，茶園間作綠肥可有效改善土壤微生物群落結構[；劉威調查了豫南地區信陽不同樹齡的茶園，發現在一定植茶年限內，種植茶樹有利于提高土壤細菌多樣性和相對豐度，其中，酸桿菌綱（Acidobacteriia）和擬桿菌綱（Bacteroidia）相對豐度變化最為顯著[7];胡磊探究了施化肥會顯著增加土層中真菌的數量，但可能會使細菌型土壤向真菌型土壤轉化[8]；馬震珠等發現蘋果園覆蓋白三葉后，土層中芽單胞菌的相對豐度增高[9]；有研究發現土壤微生物群通過礦物顆粒與有機質結合，能促進團聚體的形成[10]，對土壤結構的發育具有積極作用。茶樹作為一種葉用經濟植物，種植年限長，種植模式單一。且茶樹在長期的生存進化中形成了喜溫喜濕、喜酸嫌鈣、喜銨厭硝、聚鋁富錳等生態環境特性[11]，其微生物群落結構也獨具特色。恩施州地處30^°N “名茶緯度”之上，高海拔、低緯度、多云霧，茶園種植面積大、分布廣、品種多，各個縣市管理方式不一，土壤微生物群落結構特性的研究尚未見報道。本研究采用高通量測序技術對恩施州5個主產區的茶園土壤真菌群落組成進行分析，探明其土壤真菌群落結構與土壤理化因子間的相關性，以期為揭示恩施州不同茶園土壤真菌群落結構、真菌多樣性及其與茶園土壤環境質量關系提供理論依據。

1材料與方法

1.1 取樣位點概況

恩施州5個縣市，每個產區采樣10個，采樣點品種、海拔信息見表1。

1.2樣本采集與處理

取樣時間為2023年10月。根據采樣茶園的地形條件，按“S\"形在 05：00-08：00 時取 0～20cm 的土樣，每個采樣點沿著茶樹樹冠邊緣垂直投影取樣，適當清除地表的枯枝落葉，采用土鉆采集土壤，混合均勻，剔除雜質。一部分土樣（約 20g 裝人無菌袋及時采用液氮冷凍后，于 -80°C 保存，供分子生物學研究，另一部分土樣（約 500g 室內自然風干，過 1mm 篩網，用于土壤 ΔpH 值和有效養分的測定，然后再取部分土壤樣品研磨并過 0.1mm 篩網，用于土壤有機質和全氮含量的測定。

1.3土壤養分指標檢測方法

土壤 pH 值采用電位法測定；全氮和土壤有機質含量采用VarioMaxCN分析儀測定；土壤中有效養分含量用Mehlich-3聯合浸提劑提取（Mehlich1984），按照土：液比 1mg：10mL 混合， 180r/min 振蕩 30min ，濾液過 0.45μm 濾膜，最后采用電感耦合等離子體原子發射光譜（ICP－AES；ThermoJarrellAshLtd.，USA）測定。

1.4DNA的提取及高通量測序

根據E. Z.N.A.@ soil DNA kit（Omega Bio -tek，Norcross，GA，U.S.）說明書進行茶園土壤微生物群落總基因組DNA抽提，使用 1% 瓊脂糖凝膠電泳檢測抽提的基因組DNA的質量，使用NanoDrop2000（美國ThermoScientific公司）測定DNA濃度和純度。以該DNA為模板，選取引物ITS1F/ITS2R進行PCR（PCR 儀：ABI GeneAmp °ledast 9700 型）擴增，PCR反應體系為： 10× Buffer緩沖液 2. 0μL 2.5mmol/L dNTPs 2. 0μL ，5 μmol/L 上游引物0.8μL ，5 μmol/L 下游引物 ，TransStartFastPfuDNA聚合酶 0.4μL ，模板DNA 10ng ，補足至 20.0μL 。擴增程序如下： 95°C 預變性 3min ：95^°C 變性 30s，55‰ 退火 30s，72°C 延伸 ^45s，35 個循環； 72^°C 穩定延伸 10min，10^qC 保存。使用2% 瓊脂糖凝膠回收PCR產物，利用DNA凝膠回收純化試劑盒（PCRClean-UpKit，中國逾華）進行回收產物純化，并用Qubit4.O（ThermoFisherScientific，USA對回收產物進行檢測定量。IlluminaPE300/PE250平臺進行測序，測序后對原始序列進行拼接、質控和過濾，最后采用UPARSEv7.1軟件，設置 97% 相似度分類值對優化序列進行聚類。得到代表性序列與Silva數據庫對比注釋用于后續分析。

1.5 數據分析

利用SPSS27對數據整理與方差分析，Qiime1.9繪制真菌群落結構圖與豐度，Mothur1.30.1對序列進行 ∝ 多樣性分析，RDPClassifier2.13軟件對序列進行分類注釋，基于QIIME2（2020.2）分析平臺和R語言（3.3.1）進行物種組成、不同樣品組間的相似性和環境因子關聯分析。

2 結果與分析

2.1茶園土壤養分指標

由表2可知，恩施州5縣市茶園樣地土壤pH值在 4.61～5.48 之間，利川市茶園樣地土壤 pH 值最高，而咸豐縣茶園樣地土壤pH值最低。結合國家綠色食品產地環境質量標準（NY/T391—2021）和《茶葉產地環境技術條件（NY/T853—2004）》中茶園土壤的肥力分級標準，本研究中恩施州茶園樣地土壤pH值均分布在I級（ ΔpH 值 4.5～5.5） 茶園范圍內。有效鈣、有效鎂、有效磷含量在5縣市茶園樣地土壤中無顯著差異，有效銅、有效鐵、有效鉀、有效錳、有效鋅含量存在顯著差異（ Plt;0.05， 。有機質、全氮指標方面，恩施市茶園樣地土壤含量最高，分別達到 20.18.2.18g/kg ，全氮含量方面，咸豐縣、宣恩縣、恩施市3個樣地的茶園土壤存在顯著性差異（ Plt;0.05） 。有機質含量方面，4個縣市之間存在顯著性差異，且咸豐縣最低。

2.2茶園土壤真菌OTU及多樣性分析

50個土壤樣品通過高通量測序，共獲得8484個真菌OTU，共注釋物種1個界、17門、68綱、165目、397科、944屬、1611種。

ACE、Chaol、Coverage、Shannon、simpson等多樣性指數信息見表3。由表3可知，ACE和Chao1指數中，Xe和Hf樣本數值最大，與Es、Lc、Xf等3個縣市存在顯著差異;BDQT和EC樣本與其他品種存在顯著性差異；E1樣本與E3存在顯著差異；Shannon指數排序為 Hfgt;Lcgt;Xegt;Xfgt;Es ，Simpson指數排序依次為 LcE1gt; E2，Simpson指數海拔分組排序依次為 E3

PLS-DA分析結果（圖1）表明，不同分組樣本可明顯區分，并聚成其對應類群，表示不同縣市、不同品種、不同海拔茶園土壤真菌菌群組成具有顯著差異。從不同縣市對比結果可看出，Hf和Xe樣本聚集程度較高（圖1－A），不同品種中BDQT樣本聚集程度最高（圖1-B），不同海拔中E1與E2聚集程度較高（圖1-C）。以上結果表明，恩施州茶園土壤真菌群落結構受不同縣市管理水平、海拔、品種的影響，群落結構的構建受隨機性過程的影響較小。

2.3茶園土壤真菌群落結構分析

由圖2可知，門分類水平下，不同茶樹品種、不同海拔、不同縣市共注釋到8個門，分別是子囊菌門、擔子菌門、被孢霉門、未分類真菌、羅茲菌門、壺菌門、球囊菌門和其他。2.3.1不同縣市茶園土壤真菌群落結構分析以茶園土壤真菌群落為分析對象，設置相對豐度 gt;1% 的篩選條件，由圖2－A可知，在茶園土壤真菌豐度表現上，不同縣市樣本結果存在差異。其中，Hf共注釋到5個門，子囊菌門（55. 56% ）、擔子菌門（24.28% ）、被孢霉門（ 10. 96% ）、未分類真菌0 5.38% ）、羅茲菌門 （1.95% ）。Xf共注釋到6個門，子囊菌門 56.52% ）擔子菌門 （24.37% ）被孢霉門 （7.61% ）、未分類真菌 （5.69% ）、羅茲菌門（ 3.59% ）、壺菌門 （1.35% ）； Xe 共注釋到子囊菌門（ 65.48% ）、擔子菌門（21. 26% ）、被孢霉門（ 4.40% ）、未分類真菌（4. 27% ）、羅茲菌門（ 2.64% ）壺菌門 （1.41% ）； Lc 共注釋到6個門，子囊菌門（ 46.41% ）擔子菌門（ 20.00% ）、被孢霉門（12. 78% ）、未分類真菌（13. 06% ）、羅茲菌門（5.03% ）球囊菌門 （1.87% ）；Es共注釋到5個門，子囊菌門 33.80% ）擔子菌門 32.18% ）被孢霉門8.40% ）、未分類真菌（11. 33% ）、羅茲菌門中 13.07% ）。茶園土壤樣本相對豐度 gt;10% 水平方面，Hf分組樣本共注釋到3個門類真菌， Xf，Xe 注釋得到2個門類真菌； Lc 、Es注釋得到4個門類真菌。不同分組樣本排序表現依次為Hf（90. 81% ） gt; Xe（86.74%）gt;Xf（80.88%）gt;Lc（79.47%）gt;Es（79.05%） 。2.3.2不同品種茶園土壤真菌群落分析由圖2-B可知，在相對豐度 gt;1% 的條件下，不同品種中茶園土壤真菌豐度占比不一樣。TY共注釋到6個門，子囊菌門 （54.04% ）擔子菌門 （26.73% ）被孢霉門 （8.12% ）、羅茲菌門 （4.15% ）、未分類真菌（ 4.10% ）壺菌門 （2.45% ）。BDQT共注釋到5個門，為子囊菌門（ 49.26% ）、擔子菌門（ 23.58% ）被孢霉門 （10.45% ）、未分類真菌 （8.98% ）、羅茲菌門L 5.93% ）；EC共注釋到7個門，依次為子囊菌門（60. 45% ）擔子菌門（ 18.45% ）、未分類真菌（ 8.61% ）、被孢霉門 （6.16% ）、羅茲菌門 （2.86% ）、球囊菌門（ 1.73% ）壺菌門（ 1.44% ）;FY共注釋到5個門，依次為子囊菌門（ 46.76% ）、擔子菌門（ 27.77% ）、被孢霉門（9. 97% ）、羅茲菌門（ 7.82% ）、未分類真菌（ 6.34% ）； ΔZC 注釋到5個門，子囊菌門 47.13% ）擔子菌門 29.56% ）、未分類真菌（ 10.13% ）、被孢霉門 （7.23% ）、羅茲菌門（ 4.56% ）;ZYQ共注釋到6個門，依次是子囊菌門（51. 69% ）擔子菌門（25. 84% ）、未分類真菌（ 8.90% ）、被孢霉門（7. 43% ）、羅茲菌門（ 3.83% ）球囊菌門 （1.35% ）。不同品種排名前3的3個門累計相對豐度依次為 TY（88.89%）gt;EC （87.52%）gt;ZC（86.83%）gt;ZYQ（86.43%）gt;FY （204號 （84.50%）gt;BDQT（83.28%） ，表示優勢種群在不同品種中相對豐度占比有差異。

2.3.3不同海拔茶園土壤真菌群落結構分析由圖2-C可知，在相對豐度 gt;1% 的條件下，不同海拔茶園土壤真菌群落組成有差異，E1（ 500～800m） 共注釋到6個門，能準確分類的有5個門，分別是子囊菌門 （53.95% ）擔子菌門（ 24.45% ）、被孢霉門1 8.01% ）、未分類真菌（ 7.37% ）、羅茲菌門 （4.42% ）、壺菌門 （1.00% ）;E2（ 800～1000m 共注釋到6個門，能準確分類的有5個門，分別是子囊菌門（44.49% ）、擔子菌門（23. 38% ）、被孢霉門（ 19.86% ）、未分類真菌 （7.67% ）、羅茲菌門中 2.26% ）。優勢種群子囊菌門與擔子菌門隨海拔升高而有降低趨勢，但被孢霉門在E3，即高海拔區域的茶園相對豐度占比較高，表示不同海拔區域的茶園土壤真菌種群組成是不一樣的。

2.4恩施州茶園土壤真菌核心菌群的構建

核心微生物菌群是微生物中的關鍵組成部分，在進化過程中通過自然選擇得以保留、富集和遺傳，在不同環境樣品中發揮重要的生態作用[12]，本研究參考許國琪的方法[12]，將不同分組樣本中共存的微生物菌群定義為核心微生物菌群。不同樣本經過篩選處理后，構建的核心微生物菌群結果見圖3。由圖3可知，不同縣市茶園土壤樣本之間有核心真菌屬294屬，恩施州茶園土壤真菌核心菌群的分布表現不同。不同縣市樣本之間核心土壤微生物種群存在相同真菌，且核心真菌種類表現相似。在獨有差異真菌屬方面，5縣市結果依次表現為Xf（ 79）gt;Xe（74）gt;Hf（69）gt;Es（39）gt;Lc（31） 。不同品種方面共檢測得到核心微生物群落190個屬，占總數 40.59% ；每個品種樣本真菌屬差異明顯，其中，BDQT土壤樣本中獨有真菌125屬，占總數的26.71% ，最低的ZYQ土壤樣本中有9屬。不同海拔條件的茶園，共有核心真菌屬382屬，占總數的40.60% ，E1與E2共有真菌154屬，E1與E3共有真菌47屬，E2與E3共有16屬，E1獨有235屬，E2獨有73屬，E3僅有34屬。綜合所有樣本發現，本研究核心種群相對豐度排前5的依次為原隱球菌屬、被孢霉屬、g_unclassified_k_Fung、鐮刀菌屬等。

（ 50.63% ）、擔子菌門（23. 41% ）、被孢霉門中 （6.41% ）、未分類真菌（9. 27% ）、羅茲菌門0 8.01% ）球囊菌門（ 1.13% ）。E3（ gt;1 000m 共注釋到5個門，能準確分類的有4個門，子囊菌門

將構建出的核心微生物菌群與茶園土壤養分指標進行RDA/CCA關聯分析。由圖4可知，不同品種核心真菌群落與AFe含量表現呈顯著正相關關系，與 ΔpH 值、AAl、 ACa 、 ACu 、 AMn 、 AMg 、 AZn 、SOM含量等表現為正相關。不同縣市核心真菌群落與AK、AMg、AP、TN含量表現為負相關。不同海拔核心真菌群落與 ACu 、AFe、SOM含量呈顯著正相關，與pH值、 含量呈正相關。AK和TN含量與不同分組土壤真菌核心群落呈負相關。

2.5茶園土壤真菌群落差異性分析

2.5.1不同縣市茶園土壤真菌LEfSe分析由圖5可知，將LDA閾值設置為3進行LEfSe分析，得到恩施州主產茶5縣市茶園土壤中差異顯著的真菌群落（圖5-A），得到從門到屬水平差異豐度分類群共59個（圖5-B）。其中，Hf茶園土壤樣本中16種菌群顯著高于其他縣市，屬水平有Lectera、粉紅螺旋聚孢霉、Plectosphaerella、漢納酵母屬、Devriesia、Corynascella糞殼科未分類屬、Simplicillium、刺果蠅屬、擔子曲霉屬等。Xf茶園土壤樣本中，11種菌群顯著高于其他縣市，屬水平有擔子菌門未分類屬、擬盤多毛孢屬、Strelitziana、Cyphellophora。Xe茶園土壤樣本中有20 種菌群顯著高于其他縣市，屬水平上有枝孢屬、GS13綱未分類屬、星核衣屬、Stachybotryaceae未分類屬、烏氏霉屬、Titaea、彎孢菌屬。Lc茶園土壤樣本中有5種菌群顯著高于其他縣市，屬水平上有Neonectria、Sordariales未分類屬。Es茶園土壤樣本中有7種菌群顯著高于其他縣市。2.5.2不同品種茶園土壤真菌LEfSe分析通過LEfSe分析，不同品種分組均有分析出顯著性差異的真菌類群，由圖6可知，TY茶園土壤樣本中有9種菌群顯著高于其他品種，主要是隸屬擔子菌門未分類綱、目、科、屬4個分類群，和GS13綱及其他屬于其相關的目、科、屬3個分類群，還有1個分類群為假擬盤多毛孢屬;BDQT茶園土壤樣本中有16種菌群顯著高于其他品種，其中， LDAgt;4 的有8個分類群，分別是Plectosphaerellaceae、Lectera、古根菌綱及其相關目與科共3個分類群，線黑粉菌目、Solicoccozyma、Piskurozymaceae；剩余8個分類群分別是木耳目及其相關科、屬共3個分類群；油壺菌門及其相關綱、目、科、屬共5個分類群；EC茶園土壤樣本中有6種菌群顯著高于其他品種，古根菌綱及隸屬于其相關目、科、屬4個分類群，且LDA值 gt;4 ：還有紅酵母屬與Robillarda;FY茶園樣本中僅2種菌群顯著高于其他品種，隸屬于Trechisporales的相關的科、屬2個分類群；同樣ZYQ茶園土壤樣本中也僅有2種菌群顯著高于其他品種，小菇科及隸屬其屬的2個分類群；ZC茶園土壤樣本中有Albifimbria隸屬于Platygloeales的科與屬、Bannoa共4種菌群顯著高于其他品種。

2.5.3不同海拔茶園土壤真菌LEfSe分析E1 與E3這2個海拔梯度LEfSe分析出了顯著性差異真菌類群，由圖7可知，E1茶園土壤樣本中有3種菌群顯著高于其他海拔梯度，分別為假擬盤多毛孢屬、隸屬于肉座菌目未分類的科和屬2個分類群。但E3茶園土壤樣本中共分析出具有顯著性差異的26種菌群，分別為被孢霉門及隸屬于其相關的綱、目、科、屬9個分類群，隸屬于糞殼菌目未命名的科、屬2個分類群，白冬孢酵母目、白冬孢酵母屬、白冬孢酵母科，其他科水平上為Saccharomycetales_fam_Incertae _ sedis、Tetragoniomycetaceae、Platygloeaceae、Microsporidea_ord_Incertae_sedis、Niessliaceae5個分類群，其他屬水平上的為粒毛盤菌屬、Platygloeaceae

A為進化分支圖。圈從外到里的分類水平依次為門、綱、目、科、屬。各顏色的點表示在其組別中重要的真菌類群，黃色的點表示在各組中豐度差異不顯著的真菌類群。B為LDA值分布柱狀圖，展示的是LDA值 gt; 3的顯著差異菌群圖5不同縣市茶園土壤真菌LEfSe分析

未分類屬、鬼傘屬、菌刺孢屬、Paramicrosporidium、CercophoraMonocillium共7個分類群。精準分類的被孢霉門、被孢霉綱、被孢霉目、被孢霉科、被孢霉屬4個類群 LDAgt;4.5 。

2.6茶園土壤微生物菌群與土壤理化因子相關性分析

基于Spearman算法對茶園土壤理化因子與不同分組相對豐度前20的微生物菌群進行相關性分析。由圖8可知，小光殼屬（Leptosphaerulina）真菌與ACa含量成顯著性正相關；假散囊菌（Pseudeurotium）相對豐度在不同分組中與ACu含量呈極顯著正相關，pH值呈顯著性正相關;unclassified_o_Trechisporales相對豐度在不同分組中僅與ACu含量呈顯著性正相關；unclassified_o_GS11相對豐度僅在不同分組中與AP含量呈顯著性正相關；青霉菌屬（Penicillium）相對豐度在不同分組中與 、AFe含量呈顯著性正相關，但與AMn含量呈顯著性負相關，在不同海拔分組中還與ACa含量呈顯著性負相關;unclassified_c_Sordariomycetes相對豐度在不同縣市、不同品種中與ACu、AFe含量成顯著性正相關，在不同海拔分組中僅與ACu含量成顯著性正相關；unclassified Sordariomycetes、unclassified_p_Basidiomycota、踝節菌屬（Talaromyces）3種菌群屬的相對豐度在不同縣市分組與ACu、AFe含量呈顯著性正相關，在不同品種分組中與AFe含量成顯著性正相關，在不同海拔分組中與ACu含量呈顯著性正相關;unclassified_o_Helotiales相對豐度在不同分組中均與AAl含量呈顯著性正相關;unclassified_p_Rozellomycota相對豐度在不同分組中均與SOM含量呈顯著性正相關；但

AK含量在不同分組中對茶園真菌群落相對豐度均無影響，但其他環境因子與微生物的相對豐度均密切相關。

3討論

3.1不同縣市茶園土壤真菌群落結構分析

本研究中茶園土壤真菌多樣性結果表明，ACE和Chaol指數存在顯著性差異，Simpson、Shannon指數不存在顯著性差異，表明不同縣市茶園土壤真菌豐度存在顯著差異，而微生物種群群落差異性不大。真菌群落結構結果反映，不同縣市分組在門水平組成上均注釋到8個門，均以子囊菌門和擔子菌門為主要優勢菌群，這與其他植被類型一致。有研究發現，這2類真菌門是土壤分解的重要參與者，其中，子囊菌門的真菌多為腐生菌，在降解復雜有機質中發揮關鍵作用[13-15]，而擔子菌則是對于植物殘

A為進化分支圖。圈從外到里的分類水平依次為門、綱、目、科、屬。各顏色的點表示在其組別中重要的真菌類群，黃色的點表示在各組中豐度差異不顯著的真菌類群。B為LDA值分布柱狀圖，展示的是LDA值gt;3的顯著差異菌群體中的木質纖維素有較強的分解作用[16]。Xe子囊菌門相對豐度最高，但其有機質含量較低，這可能與土壤通透性有關，有研究表明，土壤腐生菌與土壤孔隙度呈正相關，Xe茶園建設年代較為久遠，土壤較為板結，不利于子囊菌門發揮功能。而ES子囊菌門相對豐度較低，但有機質含量最高，有研究表明茶園間作可以顯著降低土壤容重，改善土壤結構，提高有機質含量以及土壤微生物種類和數量[17]，ES習慣間作紅薯、大豆等，可能提高了子囊菌門發揮功能的潛力。茶園土壤微生物數量和組成本就受土壤類型、土壤深度、理化性質、植茶年齡、人為管理措施等因素影響比較顯著[18]，5 縣市在茶園管理方面也有各自特色，也就造就了茶園土壤真菌群落結構的差異。

3.2不同品種茶園土壤真菌群落組成分析

PLS-DA分析中，BDQT、FY、EC樣本可以明顯區分并聚成3個類群，說明3個樣本的土壤真菌菌群組成具有顯著差異，而且FY、EC樣本點分布的離散情況也表明，同一品種內真菌菌群組成也具有一定差異，這可能與植茶年限有關[19]。BDQT是當地種植最老的品種，且在LEfSe分析中BDQT有16 種菌群顯著高于其他品種， ∝ 多樣性指數也表明BDQT真菌群落豐富度和多樣性高。后來大面積發展茶園面積的時候開始大量引進FY，種植年限相對久遠，ZYQ、TY、ZC是最近20年來新進的新品種，而基于屬水平分類，不同品種具有核心物種191屬，各品種之前真菌群落不再重疊，而是單獨擁有，因此表示不同品種中茶園土壤真菌群落組成結構不黃色的點表示在各組中豐度差異不顯著的真菌類群。B為LDA值分布柱狀圖，展示的是LDA值 gt; 3的顯著差異菌群一樣。

3.3不同海拔茶園土壤真菌群落組成分析恩施州茶園主要分布在海拔為 500～800m 范圍內，占本次茶園土壤樣本調查的 56% 。本研究結果表明，真菌優勢菌門中子囊菌門、擔子菌門這2個門類真菌相對豐度隨海拔的升高而降低，而被孢霉門相對豐度在E3占比較高，LEfSe分析表明，E3中精準分類的有被孢霉門、被孢霉綱、被孢霉目、被孢霉科、被孢霉屬4個類群，且 LDAgt;4.5 。同時，從環境因子相關性分析中可發現，被孢霉屬相對豐度是受土壤養分指標影響較少，而其為好氣性真菌[20]可能高海拔的茶園土壤存在凍融現象，保持了良好的土壤通氣性。

4結論

本研究初步探索了恩施州不同縣市、不同品種、不同海拔茶園土壤真菌群落組成及多樣性特征，明確了優勢種群、核心微生物類群、顯著差異分類群；探明了除AK含量之外，其他土壤養分指標的波動將會對真菌群落結構組成有顯著性影響。本研究結果將為恩施州主產茶區各品種茶園科學管理，以及茶園土壤真菌群落結構改善提供一定的科學依據。

參考文獻：

[1]MendesR，GarbevaP，RaaijmakersJM.The rhizosphere microbiome：significanceofplant beneficial，plantpathogenic，and human pathogenic microorganisms[J]. FEMS Microbiology Reviews， 2013，37（5） ：634-663.

[2]Panke-BuisseK，PooleAC，Goodrich JK，etal．Selection on soil

microbiomesrevealsreproducibleimpactsonplant function[J].The ISMEJournal，2015，9（4）：980-989.

[3]張二豪，達瓦秋宗，徐雨婷，等.西藏濕生扁蕾根際土壤微生物群落結構特征分析[J]．西南農業學報，2023，36（12）：2741-2753.

[4]萬人源，馬會杰，蔣賓，等.茶園土壤真菌群落組成及影響因素研究［J].中國農學通報，2021，37（33）：88-97.

[5]王峰，陳玉真，吳志丹，等.有機管理模式對茶園土壤真菌群落結構及功能的影響[J].茶葉科學，2022，42（5）：672-688.

[6]王華，傅海平，周品謙，等.茶園間作綠肥對土壤微生物群落代謝功能差異的影響[J].茶葉通訊，2023，50（3）：327-332.

[7]劉威，夏雨，呂自敏，等.豫南不同樹齡茶園土壤微生物群落結構及多樣性分析[J].南方農業學報，2023，54（5）：1455-1464.

[8]胡磊.套種圓葉決明和施肥對茶園土壤固氮微生物群落的影響[D].福州：福建農林大學，2010：48-49.

[9]馬震珠，魏倩倩，李尚瑋，等.覆蓋和埋置白三葉對黃土高原蘋果園土壤細菌多樣性的影響[J]．水土保持通報，2019，39（5）：143-151.

[10]OrtizC，Fernández-AlonsoMJ，KitzlerB，etal.Variationsin soilaggregation，microbial community structureand soil organicmattercyclingassociated tolong-term afforestation andwoodyencroachment in a Mediterranean alpine ecotone[J].Geoderma，2022，405：115450.

[11]楊廣容，王秀青，謝瑾，等.云南古茶園和現代茶園土壤養分與茶葉品質成分關系的研究[J]．茶葉科學，2015，35（6）：574-582.

[12]許國琪，劉怡萱，曹鵬熙，等.青藏高原冰川棘豆（Oxytropisglacialis）內生菌核心微生物組的界定及其互作網絡分析［J].微生物學通報，2020，47（9）：2746-2758.

[13]阿拉坦存布爾，曹烏吉斯古楞，包金花，等.內蒙古地區手參根際土壤真菌群落組成及多樣性[J]．西南農業學報，2023，36（12）：2754-2762.

[14]劉煒璇，李依蒙，江紅星，等.吉林莫莫格國家級自然保護區四種典型植物群落下土壤微生物組成的對比分析［J]．生態學雜志，2024，43（10） ：2988-2998.

[15]BeimfordeC，FeldbergK，NylinderS，etal.EstimatingthePhanerozoic history of the ascomycota lineages：combining fossil andmoleculardata[J].MolecularPhylogeneticsandEvolution，2014，78：386-398.

[16]FreySD，KnorrM，ParrentJL，etal.Chronicnitrogenenrichmentaffects the structure and function of the soil microbial communityintemperate hardwood and pine forests[J]．Forest EcologyandManagement，2004，196（1） ：159-171.

[17]de Costa WAJM，SurenthranP.Tree-crop interactions inhedgerow intercroppingwith different treespeciesand tea in SriLanka：1.production and resource competition[J].AgroforestrySystems，2005，63（3）：199-209.

[18]薛冬．茶園土壤微生物群落多樣性及硝化作用研究[D]．杭州：浙江大學，2007：2-8.

[19]姚澤秀，李永春，李永夫，等.植茶年限對土壤微生物群落結構及多樣性的影響[J].應用生態學報，2020，31（8）：2749-2758.

[20]羅伶書，王一佩，杜盛.模擬降雨減少對刺槐人工林土壤養分和微生物群落的影響[J]．生態學報，2023，43（14）：5916-5925.