喀斯特地區不同海拔對油茶根際土壤叢枝菌根真菌物種群落的影響

曾維軍 劉燕 秦維 羅文敏

摘要：本研究旨在了解喀斯特地區不同海拔油茶根際土壤叢枝菌根真菌（ａｒｂｕｓｃｕｌａｒ ｍｙｃｏｒｒｈｉｚａｌ ｆｕｎｇｉ，ＡＭＦ）的種群多樣性及物種分布特征，以提高ＡＭＦ 在喀斯特地區生態修復中的應用基礎。 以貴州喀斯特高、低海拔地區油茶根際土壤為材料，利用Ｉｌｌｕｍｉｎａ ＭｉＳｅｑ 測序平臺進行高通量測序，分析ＡＭＦ 物種多樣性及差異。 結果表明，喀斯特地區油茶根際土壤ＡＭＦ 中球囊霉屬（Ｇｌｏｍｕｓ）占主要地位，整體豐度為８７.２８％；低海拔油茶根際土壤ＡＭＦ 群落特異ＯＴＵ 數和Ｓｈａｎｎｏｎ、Ｃｈａｏ１ 指數顯著高于高海拔，物種多樣性較為豐富，且組間ＡＭＦ 物種聚類較遠；油茶根際土壤有機質和有效磷含量與球囊霉屬真菌豐度呈顯著正相關，相關系數分別為０.５８、０.６５，表明喀斯特的低磷土壤促進油茶根際球囊霉屬的生長，增加土壤碳庫，促進油茶對磷的吸收利用。其作用機理待進一步研究驗證。

關鍵詞：油茶；喀斯特地區；根際土壤；叢枝菌根真菌；多樣性

中圖分類號：Ｓ１５４.３文獻標識號：Ａ文章編號：１００１－４９４２（２０２４）０１－０１０５－０８

貴州省６２％的國土面積為喀斯特地貌，總面積達１０.９１ 萬ｋｍ２，居全國之首，是全球喀斯特發育最典型、最復雜、景觀類型最豐富的地區之一[１] 。 喀斯特地貌地表干旱缺水，容易造成石漠化等生態退化問題[２] 。 通過植被恢復緩解石漠化是該地區生態修復的關鍵，但土壤的缺水和貧瘠導致植物難以定居、生長緩慢[３] 。 油茶是具有重要經濟和生態價值的油料植物，貴州引種面積逐年增加，栽培過程中也面臨著土壤保水保肥能力弱等問題，制約著油茶果產量和質量的提升。尋找一種提高油茶抗逆性的高效方法成為解決該問題的關鍵[４－５] 。

近年來研究發現，廣泛定植于植物根部的叢枝菌根真菌（ＡＭＦ） 具有促進宿主營養吸收、抗病、抗旱和改善土壤微生態等功能[６－９] 。 世界上許多國家和地區也規定在特定環境下造林，必須采用菌根化的苗木，以增加移栽的存活率[１０－１１] 。劉燕等[１２] 前期研究也發現，耐旱且油茶果產量穩定的貴州野生油茶，其根系普遍被ＡＭＦ 定殖。 因此，對ＡＭＦ 的開發應用是未來研究的熱點和趨勢。

ＡＭＦ 的開發利用過程中存在兩個主要問題，一是不同ＡＭＦ 對環境適應性存在差異，二是其高效培養和接種技術有待突破[１３－１４] 。 貴州喀斯特高、低海拔地區都有野生油茶林分布，不同海拔的油茶可能形成了自己特有的根際ＡＭＦ 群落。 本研究通過Ｉｌｌｕｍｉｎａ ＭｉＳｅｑ 平臺對喀斯特高、低海拔地區野生油茶根際土樣進行高通量測序，探究ＡＭＦ 的物種多樣性及其在土壤營養形成中的作用，為ＡＭＦ 在貴州喀斯特地區油茶綠色種植中的應用提供理論依據。

１ 材料與方法

１.１ 研究區概況

設置高海拔和低海拔２ 個研究區，海拔高度３００～２ ４２１ ｍ。 研究區屬亞熱帶季風性和季風性濕潤氣候，年平均氣溫１１.５ ～ １９.２ ℃，年降水量９２６.２～１ ３４０.７ ｍｍ，年日照時數為１ ３２８ ～ １ ８１２ｈ，無霜期１８０～３４５ ｄ。 樣區為野生油茶林，林齡為４０～５０ 年，經度在１０４.０８° Ｅ～１０９.０７° Ｅ，緯度在２４.９３° Ｎ～２７.４７° Ｎ。 每樣地設３ 個野生油茶根際土壤采樣點，樣地詳情見表１。

１.２ 土壤樣品采集

２０２１ 年８ 月，在貴州省威寧、冊亨、從江、玉屏縣的野生油茶分布區內各設置規格為４００ ｍ２的３ 個相鄰樣方（２０ ｍ×２０ ｍ），在每個樣方內沿“Ｓ”形路徑分別選擇１０ 株胸徑和樹高相近的油茶樹，鏟除落葉層并挖開上層覆土，剪下細根分枝并輕輕抖動，用毛刷將黏在細根上的土壤刷落到自封袋中。 每個樣方內取１０ 份根際土壤，混勻后４ ℃冰箱保存，用于測定ＡＭＦ 多樣性。 樣品的采集與預處理參照?土壤微生物研究原理與方法?[１５] 。

１.３ 土壤樣品測序

根據ｓｏｉｌ ＤＮＡ ｋｉｔ（Ｏｍｅｇａ Ｂｉｏ－ｔｅｋ， Ｎｏｒｃｒｏｓｓ，ＧＡ，Ｕ.Ｓ.）使用說明進行總ＤＮＡ 提取。 １％的瓊脂糖凝膠電泳檢測基因組ＤＮＡ 質量，使用Ｎａｎｏ￣Ｄｒｏｐ２０００（Ｔｈｅｒｍｏ Ｓｃｉｅｎｔｉｆｉｃ 公司，美國）測定ＤＮＡ濃度和純度。 用兩對特異性引物ＡＭＬ１/ ＡＭＬ２（５′ － ＡＴＣＡＡＣＴＴＴＣＧＡＴＧＧＴＡＧＧＡＴＡＧＡ － ３′；５′－ＧＡＡＣＣＣＡＡＡＣＡＣＴＴＴＧＧＴＴＴＣＣ－３′）和ＡＭＶ４－５ＮＦ/ＡＭＤＧＲ（５′－ＡＧＣＴＣＧＴＡＧＴＴＧＡＡＴＴＴＣＧ－３′；５′－ＣＣＣＡＡＣＴＡＴＣＣＣＴＡＴＴＡＡＴＣＡＴ － ３′） 進行巢式ＰＣＲ 擴增[１６] 。 ＰＣＲ 反應體系為： ５ × ＴｒａｎｓＳｔａｒｔＦａｓｔＰｆｕ 緩沖液４ mＬ，２.５ ｍｍｏｌ/ Ｌ ｄＮＴＰｓ ２ mＬ，上、下游引物（５ mｍｏｌ/ Ｌ）各０.８ mＬ，ＴｒａｎｓＳｔａｒｔ ＦａｓｔＰｆｕＤＮＡ 聚合酶０.４ mＬ，ＤＮＡ 模板１０ mＬ，超純水補足至２０ mＬ。 擴增程序如下：９５ ℃ 預變性３ ｍｉｎ；９５ ℃變性３０ ｓ，５５ ℃退火３０ ｓ， ７２ ℃延伸４５ ｓ，

２７ 個循環；７２ ℃延伸１０ ｍｉｎ；４ ℃保存。 每樣本３個重復。 將同一樣本的ＰＣＲ 產物混合后使用２％瓊脂糖凝膠回收，利用ＡｘｙＰｒｅｐ ＤＮＡ Ｇｅｌ Ｅｘｔｒａｃ￣ｔｉｏｎ Ｋｉｔ（Ａｘｙｇｅｎ Ｂｉｏｓｃｉｅｎｃｅｓ，Ｕｎｉｏｎ Ｃｉｔｙ，ＣＡ，ＵＳＡ）純化，２％瓊脂糖凝膠電泳檢測其純度，并用Ｑｕａ￣ｎｔｕｓＴＭ Ｆｌｕｏｒｏｍｅｔｅｒ（Ｐｒｏｍｅｇａ， ＵＳＡ）對回收產物進行檢測定量。 使用ＮＥＸＴＦＬＥＸ Ｒａｐｉｄ ＤＮＡ－ＳｅｑＫｉｔ 對純化后的ＰＣＲ 產物進行建庫。 委托上海美吉生物醫藥科技有限公司利用Ｉｌｌｕｍｉｎａ ＭｉＳｅｑ 平臺進行測序。

１.４ 土壤理化性質檢測

土壤全氮、全磷、全鉀、有效磷、速效鉀含量及ｐＨ 值分別參照ＮＹ/ Ｔ ５３—１９８７、ＧＢ/ Ｔ ９８３７—１９８８、ＮＹ/ Ｔ ８７—１９８８、ＮＹ/ Ｔ １１２１.７—２０１４、ＮＹ/ Ｔ１８４９—２０１０、ＮＹ/ Ｔ － １２４１—１９９９ 中方法進行測定；堿解氮采用堿解擴散法測定[１７] 。

１.５ 數據統計與分析

使用Ｆｌａｓｈ１.２.１１ 進行ｐａｉｒ－ｅｎｄ 雙端序列拼接；ＲＤＰ Ｃｌａｓｓｉｆｉｅｒ２.１３ 進行序列分類注釋；Ｕｓｅ￣ａｒｃｈ１１ 進行ＯＴＵ 統計；Ｍｏｔｈｕｒ１.３０.２ 進行Ａｌｐｈａ多樣性分析，并對樣本間的多樣性指數進行ｔ 檢驗；Ｍｅｇａ７ 和ＩＱ－ＴＲＥＥ 進行進化樹分析和構建；Ｒ語言（ｖｅｒｓｉｏｎ ３.３.１）進行物種Ｖｅｎｎ 圖分析、制作群落柱形圖以及ＮＭＤＳ 分析和作圖。 最后使用Ａｄｏｂｅ Ｉｌｌｕｓｔｒａｔｏｒ２０２２ 修飾圖形。 其中聚類方式采用ＵＳＥＡＲＣＨ７－ｕｐａｒｓｅ 算法，物種分類數據庫為ｍａａｒｊａｍ０８１/ ＡＭＯＴＵ。

２ 結果與分析

２.１ 不同海拔野生油茶根際土壤ＡＭＦ 群落ＯＴＵ組成及結構

不同海拔野生油茶根際土壤中共獲得３６７ １４８條高質量序列，樣品平均序列為２０ ３９７條，總堿基數為７９ ３３８ ２１４ ｂｐ，平均序列長度為２１６ ｂｐ，序列相似度為０.９７，ＯＴＵ 分類置信度為０.７。由稀釋曲線（圖１）可知，抽取的序列條數達到１４ ２２３條以上，曲線趨于平坦，表明不同處理的土樣樣本測序數量足夠，測序數據合理。 由Ｖｅｎｎ圖（圖２）可知，喀斯特高、低海拔地區ＡＭＦ 群落共有的ＯＴＵ 數為６１ 個，占總數的１８.３２％；高海拔地區特異ＯＴＵ 數為５７ 個，占總數的１７.１２％；低海拔地區特異ＯＴＵ 數為１５４ 個， 占總數的４６.２５％。

２.２ 不同海拔野生油茶根際土壤ＡＭＦ 群落結構及多樣性

對喀斯特野生油茶根際土壤ＡＭＦ 群落進行物種注釋得到１ 綱５ 目８ 科１４ 屬７５ 種。 由圖３可知，在屬水平上，相對豐度排前４ 位的屬分別為球囊霉屬（Ｇｌｏｍｕｓ，８７.２８％）、無梗囊霉屬（Ａｃａｕ￣ｌｏｓｐｏｒａ，５. ９７％）、盾巨孢囊霉屬（ Ｓｃｕｔｅｌｌｏｓｐｏｒａ，２.７８％）、球囊菌綱中不能注釋的菌屬（ｕｎｃｌａｓｓｉｆｉｅｄ＿ｃ＿Ｇｌｏｍｅｒｏｍｙｃｅｔｅｓ，２.１２％），球囊霉屬為優勢屬。喀斯特高、低海拔地區相比，球囊霉屬相對豐度差異不顯著，無梗囊霉屬和球囊菌綱中不能注釋的菌屬在高海拔地區相對豐度較高。 高海拔地區未發現盾巨孢囊霉屬。不同海拔地區野生油茶根際土壤ＡＭＦ 群落α 多樣性指數如表２ 所示，高海拔地區Ｓｈａｎｎｏｎ指數和Ｃｈａｏ１ 指數均顯著低于低海拔地區，ＳＨＥＩ指數和ＰＤ 指數與低海拔地區相比差異不顯著。表明喀斯特高海拔地區的野生油茶根際ＡＭＦ 群落的多樣性和物種總數均顯著低于低海拔地區。

２.３ 不同海拔野生油茶根際土壤ＡＭＦ 群落ＰＣｏＡ 分析與樣本聚類分析

ＰＣｏＡ 分析可通過降維找出影響野生油茶根際土壤ＡＭＦ 樣本群落組成差異的潛在主成分。由圖４ 可知，在種水平上，主成分１（ＰＣ１）與主成分２（ＰＣ２）解釋方差分別為５１.０５％和３２.５９％；在ＯＴＵ 水平上，主成分１（ＰＣ１）與主成分２（ＰＣ２）解釋方差分別為２４.６０％和２２.８６％，影響野生油茶根際土壤ＡＭＦ 群落結構的主導因子明顯（Ｐ <０.０５）。在種水平上，同一海拔、同一地點根際土壤３ 個重復比較聚集，不同地點之間分化明顯，高海拔樣點之間的物種分化較低海拔明顯。 在ＯＴＵ水平上聚類情況與種水平上相反。

非度量多維尺度分析（ＮＭＤＳ）可將多維空間的研究對象簡化到低維空間進行定位、分析和歸類，揭示復雜背景下的簡單規律。 由圖５ 可知，在種水平上，高海拔野生油茶根際土壤樣本組的ＡＭＦ 群落包含了低海拔土壤樣本組，且Ｓｔｒｅｓｓ<０.０５，具有很好的代表性；在ＯＴＵ 水平上，不同海拔樣本組間存在差異，且Ｓｔｒｅｓｓ<０.２。 ＰＣｏＡ 與ＮＭＤＳ 分析共同驗證了貴州喀斯特高、低海拔間野生油茶根際土壤ＡＭＦ 群落組成差異明顯，高海拔樣本組涵蓋了低海拔的ＡＭＦ 分子種。

基于ｗｅｉｇｈｔｅｄ ｎｏｒｍａｌｉｚｅｄ ｕｎｉｆｒａｃ 距離算法對貴州喀斯特高、低海拔野生油茶根際土壤樣本進行層級聚類分析（圖６）可知，高海拔地區Ａ３ 和Ａ６ 樣本組的ＡＭＦ 種群相似，聚為一類。 高、低海拔的ＡＭＦ 物種組成豐度上有較大差異，說明海拔高低能改變ＡＭＦ 種群結構。 但是低海拔樣本組Ｅ２ 與高海拔Ａ２ 組的ＡＭＦ 種群相似，聚為一類。威寧縣草海鎮（Ａ３）與從江縣拱孑村（Ｄ２）兩地海拔差為２ ０５２ ｍ，其ＡＭＦ 群落結構相似性最低。

２.４ ＡＭＦ 群落與土壤養分因子的相關性分析

通過方差膨脹因子分析（ＶＩＦ）對７ 個影響樣本ＡＭＦ 菌群組成的土壤養分因子進行篩選，保留多重共線性較小的因子，進行后續研究。 由表３可知，篩選后的土壤養分因子分別為全磷、全鉀、有效磷、速效鉀和有機質，ＶＩＦ 值均小于１０。

通過計算篩選后的土壤養分因子與ＡＭＦ 物種之間的Ｐｅａｒｓｏｎ 相關系數，獲得相關性熱圖（圖７）。 結果表明，球囊霉屬與土壤有機質和有效磷呈顯著正相關，相關系數分別為０.５８、０.６５；多樣孢囊霉屬（Ｄｉｖｅｒｓｉｓｐｏｒａ） 和原囊霉屬（Ａｒｃｈａｅｏｓｐｏ￣ｒａ）與全磷呈極顯著正相關，相關系數分別為０.９１、０.９０；盾巨孢囊霉屬與速效鉀呈顯著負相關，相關系數為－０.６１；無梗囊霉屬（Ａｃａｕｌｏｓｐｏｒａ）與全鉀呈極顯著負相關，相關系數為－０.７８。

３ 討論

除宿主特異性外，ＡＭＦ 的分布還受氣候、土壤條件、緯度、空間隔離、干擾等多種非生物因子的影響[１８－２２] 。 研究發現，隨著海拔的升高，青藏高原草甸ＡＭＦ 的物種豐度逐漸降低[２２] ；在宿主植物和土壤條件相似的情況下，ＡＭＦ 多樣性隨氣溫升高而增加[２３－２４] 。 本研究發現，低海拔地區油茶根際ＡＭＦ 群落特異ＯＴＵ 數、Ｓｈａｎｎｏｎ 指數、Ｃｈａｏ１ 指數顯著高于高海拔地區，表明ＡＭＦ 物種總數和多樣性隨海拔的升高而降低；同時ＡＭＦ 的群落結構也發生了改變，與低海拔地區相比，高海拔地區油茶根際土壤中未發現盾巨孢囊霉屬真菌。 本研究中，海拔高度的改變直接導致兩地區間較大的年平均氣溫差，達７.７ ℃，說明野生油茶根際土壤ＡＭＦ 的物種結構和多樣性差異也可能是由緯度、海拔差造成的環境溫差導致。 另外，高、低海拔地區的油茶根際ＡＭＦ 物種可以通過聚類分析分別聚在一起，同時高、低海拔組間又有一定的物種隔離，說明在喀斯特地區海拔差異越大，油茶根際ＡＭＦ 物種的親緣關系就越遠。

石灰土是喀斯特地區分布最廣的非地帶性土壤，表現出高鈣鎂、低磷鉀和土壤偏堿性的地球化學特征[２５] ，而在一定范圍內，低磷土壤會促進ＡＭＦ 的生長[２６－２７] 。 研究表明，球囊霉屬真菌能分泌球囊霉素，作為一種非專性蛋白，其主要作用是穩定和增加土壤有機碳庫、改善土壤團聚體，還可以促進植物對磷的吸收，并能與鉛結合在土壤重金屬形態轉化中具有潛在的重要意義[２８－３０] 。 本研究中，球囊霉屬在ＡＭＦ 中占據主要地位，并與油茶根際土壤有機質和有效磷含量呈顯著正相關，說明喀斯特的低磷土壤能提高油茶根際球囊霉屬豐度，使土壤中的球囊霉素含量理論上得到提高，從而改善土壤質量，增加土壤碳庫，促進油茶對磷的吸收。

４ 結論

貴州喀斯特高、低海拔地區的野生油茶根際土壤ＡＭＦ 多樣性、群落結構、種水平上的聚類關系均具有顯著差異，其中低海拔地區的ＡＭＦ 多樣性較豐富，并與高海拔地區存在較為明顯的種群親緣隔離。 另外，喀斯特地區特有的石灰土中缺乏磷、鉀，促進了油茶根際ＡＭＦ 的生長，尤其是球囊霉屬，其分泌較多的球囊霉素擴撒到土壤中，可改善土壤質量、增加土壤碳庫、促進油茶對磷的吸收利用，從而使油茶根際土壤有機質和有效磷含量增加。 其相關作用機理有待進一步研究驗證。

參 考 文 獻：

[１] Ｔａｎｇ Ｊ， Ｔａｎｇ Ｘ Ｘ， Ｑｉｎ Ｙ Ｍ， ｅｔ ａｌ. Ｋａｒｓｔ ｒｏｃｋｙ ｄｅｓｅｒｔｉｆｉｃａｔｉｏｎｐｒｏｇｒｅｓｓ： ｓｏｉｌ ｃａｌｃｉｕｍ ａｓ ａ ｐｏｓｓｉｂｌｅ ｄｒｉｖｉｎｇ ｆｏｒｃｅ[Ｊ]. Ｓｃｉｅｎｃｅｏｆ ｔｈｅ Ｔｏｔａｌ Ｅｎｖｉｒｏｎｍｅｎｔ， ２０１９， ６４９：１２５０－１２５９.

[２] 王世杰. 喀斯特石漠化———中國西南最嚴重的生態地質環境問題[Ｊ]. 礦物巖石地球化學通報，２００３，２２（２）：１２０－１２６.

[３] 魏源，王世杰，劉秀明，等. 叢枝菌根真菌及在石漠化治理中的應用探討[Ｊ]. 地球與環境，２０１２，４０（１）：８４－９２.

[４] 黃安香，柏文戀，鄔能英，等. 貴州油茶主產區土壤重金屬環境質量狀況及風險評價[Ｊ]. 環境化學， ２０２３，４２（４）：１３２５－１３３６.

[５] 柏文戀，黃安香，楊守祿，等. 貴州油茶主產區土壤養分特征及綜合肥力評價[Ｊ]. 西南農業學報，２０２２，３５（９）：２０９６－２１０６.

[６] 劉代軍，涂波，施松梅，等. 石漠化地區的生態危機及菌根桑生物修復潛力研究進展[Ｊ]. 中國巖溶，２０１２，３１（２）：１８５－１９０.

[７] Ｄｕｒｕ Ｍ， Ｃｒｕｚ Ｐ Ｐ， Ｒａｏｕｄａ Ａ Ｈ Ｋ， ｅｔ ａｌ. Ｒｅｌｅｖａｎｃｅ ｏｆ ｐｌａｎｔｆｕｎｃｔｉｏｎａｌ ｔｙｐｅｓ ｂａｓｅｄ ｏｎ ｌｅａｆ ｄｒｙ ｍａｔｔｅｒ ｃｏｎｔｅｎｔ ｆｏｒ ａｓｓｅｓｓｉｎｇｄｉｇｅｓｔｉｂｉｌｉｔｙ ｏｆ ｎａｔｉｖｅ ｇｒａｓｓ ｓｐｅｃｉｅｓ ａｎｄ ｓｐｅｃｉｅｓ￣ｒｉｃｈ ｇｒａｓｓｌａｎｄｃｏｍｍｕｎｉｔｉｅｓ ｉｎ ｓｐｒｉｎｇ[Ｊ]. Ａｇｒｏｎｏｍｙ Ｊｏｕｒｎａｌ， ２００８， １００（６）：１６２２－１６３０.

[８] 姜海燕，李敏. 西瓜－辣椒輪作體系ＡＭＦ 和ＰＧＰＲ 防控線蟲的效果[Ｊ]. 山東農業科學，２０２１，５３（４）：１３１－１３４.

[９] 康佳，梁夕金，楊文龍，等. 叢枝菌根真菌（ＡＭＦ）對鹽堿地花生根系土壤微生物的影響[Ｊ]. 山東農業科學，２０２２，５４（７）：７７－８４.

[１０] Ｓｋｕｊｉ?? Ｊ， Ａｌｌｅｎ Ｍ Ｆ. Ｕｓｅ ｏｆ ｍｙｃｏｒｒｈｉｚａｅ ｆｏｒ ｌａｎｄ ｒｅｈａｂｉｌｉｔａｔｉｏｎ[Ｊ]. ＭＩＲＣＥＮ Ｊｏｕｒｎａｌ ｏｆ Ａｐｐｌｉｅｄ Ｍｉｃｒｏｂｉｏｌｏｇｙ ａｎｄ Ｂｉｏｔｅｃｈｎｏｌ￣ｏｇｙ， １９８６， ２：１６１－１７６.

[１１] 魏源，王世杰，劉秀明，等. 喀斯特地區叢枝菌根真菌遺傳多樣性[Ｊ]. 生態學雜志，２０１１，３０（１０）：２２２０－２２２６.

[１２] 劉燕，曾維軍，秦維，等. 威寧短柱油茶叢枝菌根真菌和深色有隔內生真菌定殖特征及其與果實經濟性狀的相關性[Ｊ]. 江蘇農業科學，２０２３，５１（２２）：１４５－１５３.

[１３] 楊佳絨，戴冬，陳俊芳，等. 叢枝菌根真菌多樣性對植物群落構建和稀有種維持的研究進展[ Ｊ]. 植物生態學報，２０２３，４７（６）：７４５－７５５.

[１４] 儲薇，郭信來，張晨，等. 叢枝菌根真菌－植物－根際微生物互作研究進展與展望[Ｊ]. 中國生態農業學報：中英文，２０２２，３０（１１）：１７０９－１７２１.

[１５] 林先貴. 土壤微生物研究原理與方法[Ｍ]. 北京：高等教育出版社，２０１０：１１２－１３４.

[１６] Ｖａｎ Ｇｅｅｌ Ｍ， Ｂｕｓｓｃｈａｅｒｔ Ｐ， Ｈｏｎｎａｙ Ｏ， ｅｔ ａｌ. Ｅｖａｌｕａｔｉｏｎ ｏｆ ｓｉｘｐｒｉｍｅｒ ｐａｉｒｓ ｔａｒｇｅｔｉｎｇ ｔｈｅ ｎｕｃｌｅａｒ ｒＲＮＡ ｏｐｅｒｏｎ ｆｏｒ ｃｈａｒａｃｔｅｒｉｚａ￣ｔｉｏｎ ｏｆ ａｒｂｕｓｃｕｌａｒ ｍｙｃｏｒｒｈｉｚａｌ ｆｕｎｇａｌ （ＡＭＦ） ｃｏｍｍｕｎｉｔｉｅｓ ｕ￣ｓｉｎｇ ４５４ ｐｙｒｏｓｅｑｕｅｎｃｉｎｇ[Ｊ]. Ｊｏｕｒｎａｌ ｏｆ Ｍｉｃｒｏｂｉｏｌｏｇｉｃａｌ Ｍｅｔｈ￣ｏｄｓ， ２０１４， １０６：９３－１００.

[１７] 馬獻發. 土壤學實驗[Ｍ]. 北京：中國林業出版社，２０２０：２１８－２３４.[１８] Ｂｅｖｅｒ Ｊ Ｄ. Ｈｏｓｔ￣ｓｐｅｃｉｆｉｃｉｔｙ ｏｆ ＡＭ ｆｕｎｇａｌ ｐｏｐｕｌａｔｉｏｎ ｇｒｏｗｔｈｒａｔｅｓ ｃａｎ ｇｅｎｅｒａｔｅ ｆｅｅｄｂａｃｋ ｏｎ ｐｌａｎｔ ｇｒｏｗｔｈ [ Ｊ]. Ｐｌａｎｔ ａｎｄＳｏｉｌ， ２００２， ２４４：２８１－２９０.

[１９] 楊海水. 宿主植物對叢枝菌根真菌的影響———共生功能、地理分布及多樣性[Ｄ]. 杭州：浙江大學，２０１３.

[２０] 楊蕊毓，鄧錡璋，田麗平，等. 不同生境下川麥冬根圍土壤叢枝菌根真菌多樣性[Ｊ]. 西北植物學報，２０２２，４２（１）：１４５－１５３.

[２１] 朱亮，郭可馨，藍麗英，等. 亞高山森林類型轉換對土壤叢枝菌根真菌多樣性的影響[Ｊ]. 生態學雜志，２０２０，３９（１２）：３９４３－３９５１.

[２２] Ｇａｉ Ｊ Ｐ， Ｔｉａｎ Ｈ， Ｙａｎｇ Ｆ Ｙ， ｅｔ ａｌ. Ａｒｂｕｓｃｕｌａｒ ｍｙｃｏｒｒｈｉｚａｌ ｆｕｎ￣ｇａｌ ｄｉｖｅｒｓｉｔｙ ａｌｏｎｇ ａ Ｔｉｂｅｔａｎ ｅｌｅｖａｔｉｏｎ ｇｒａｄｉｅｎｔ[Ｊ]. Ｐｅｄｏｂｉｏｌｏ￣ｇｉａ， ２０１２， ５５（３）：１４５－１５１.

[２３] Ｋｏｓｋｅ Ｒ Ｅ. Ｄｉｓｔｒｉｂｕｔｉｏｎ ｏｆ ＶＡ ｍｙｃｏｒｒｈｉｚａｌ ｆｕｎｇｉ ａｌｏｎｇ ａ ｌａｔｉｔｕ￣ｄｉｎａｌ ｔｅｍｐｅｒａｔｕｒｅ ｇｒａｄｉｅｎｔ[Ｊ]. Ｍｙｃｏｌｏｇｉａ， １９８７， ７９：５５－６８.

[２４] 石國璽，王芳萍，馬麗，等. 長期、短期增溫對高寒草甸ＡＭ真菌群落結構的影響[Ｊ]. 草地學報，２０２１，２９（增刊１）：１７９－１８９.

[２５] 羅緒強，王程媛，楊鴻雁，等. 喀斯特優勢植物種干旱和高鈣適應性機制研究進展[Ｊ]. 中國農學通報，２０１２，２８（１６）：１－５.

[２６] Ｔａｗａｒａｙａ Ｋ， Ｓａｉｔｏ Ｍ， Ｍｏｒｉｏｋａ Ｍ， ｅｔ ａｌ. Ｅｆｆｅｃｔ ｏｆ ｐｈｏｓｐｈａｔｅａｐｐｌｉｃａｔｉｏｎ ｔｏ ａｒｂｕｓｃｕｌａｒ ｍｙｃｏｒｒｈｉｚａｌ ｏｎｉｏｎ ｏｎ ｔｈｅ ｄｅｖｅｌｏｐｍｅｎｔａｎｄ ｓｕｃｃｉｎａｔｅ ｄｅｈｙ￣ｄｒｏｇｅｎａｓｅ ａｃｔｉｖｉｔｙ ｏｆ ｉｎｔｅｒｎａｌ ｈｙｐｈａｅ[Ｊ].Ｓｏｉｌ Ｓｃｉｅｎｃｅ ａｎｄ Ｐｌａｎｔ Ｎｕｔｒｉｔｉｏｎ， １９９４， ４０（４）：６６７－６７３.

[２７] 吳強盛，夏仁學，鄒英寧. 柑橘叢枝菌根真菌生長與根際有效磷和磷酸酶活性的相關性[Ｊ]. 應用生態學報，２００６，１７（４）：４６８５－４６８９.

[２８] 張健，徐明，吳麗豐，等. 干擾對黔中地區馬尾松林土壤球囊霉素的影響[Ｊ]. 森林與環境學報，２０２２，４２（１）：６３－７０.

[２９] 張健，徐明，王陽，等. 黔中地區不同馬尾松群叢土壤球囊霉素分布特征[Ｊ]. 生態環境學報，２０２１，３０（１２）：２３０３－２３０８.

[３０] 甘佳偉，韓曉增，鄒文秀. 球囊霉素及其在土壤生態系統中的作用[Ｊ]. 土壤與作物，２０２２，１１（１）：４１－５３.