科爾沁沙地樟子松根內真菌群落結構和功能群特征*

趙珮杉 郭米山 高廣磊 丁國棟 張 英

(1.北京林業大學水土保持學院 水土保持國家林業和草原局重點實驗室 北京 100083；2.寧夏鹽池毛烏素沙地生態系統國家定位觀測研究站 鹽池 751500)

根際微生物是連接根系與土壤的重要紐帶(Puglisietal.， 2013； Becketal.， 2015)，根內真菌定殖(或部分定殖)于植物根內，其多樣性和群落結構與宿主植物的生長、發育和更新能力存在密切關系(Philippotetal.， 2013； Baldrian， 2017)。一般來說，根內真菌可分為寄生真菌、腐生真菌和共生真菌3種類型(Bonitoetal.， 2016； Kilpelainenetal.， 2016)。其中，寄生真菌誘發植物病害，威脅森林生態系統健康； 腐生真菌是生態系統重要的分解者，參與生態系統物質轉化和能量循環(Mooreetal., 2004)； 共生真菌主要為菌根真菌，能夠促進宿主植物吸收水分和礦物質養分，協助植物提高抗逆性和適應性(Frewaetal.， 2018； Martinetal.， 2017)。因此，根內真菌組成及其功能直接影響宿主植物的生長和健康狀況，與生態系統穩定性存在密切聯系。

一般來說，根內真菌對宿主植物常表現出特定偏好，具有宿主特異性，這一現象在植物組成豐富多樣的天然次生林中尤為明顯(Tojuetal., 2013a； Picklesetal., 2015)。但在干旱荒漠地區，降水稀少，植被單一、稀疏，部分根內真菌也表現出較強的普適性，能廣泛定殖于多種宿主植物(Botnenetal.， 2014； Mergesetal., 2018)。根內真菌分布及其選擇宿主植物的影響因素較為復雜，既依賴于宿主植物特性，又與環境條件密切相關。在大空間尺度水平，氣候和生境條件是造成植物根內真菌組成差異的主要原因(Davidetal.， 2016； Barnesetal.， 2016a)； 而在小空間尺度水平，土壤理化性質則是根內真菌物種豐度和組成變化的主要影響因素(Gaoetal., 2015； Liuetal.， 2015)。林齡即是林分的重要屬性特征，又對林下土壤理化性質有重要影響(Gaoetal.， 2014)。已有研究表明，林齡是影響根內真菌定殖的重要因子(Kyaschenkoetal.， 2017)，宿主植物受病原菌侵害程度與林齡密切相關(Sangüesa-Barredaetal.， 2015)，菌根真菌的侵染率、多樣性和群落組成均隨林齡的變化而變化(陳雪冬等， 2017； Spakeetal.， 2016)，在不同演替階段根內真菌的優勢菌種存在差異(Daveyetal.， 2015)。因此，針對干旱荒漠地區不同林齡固沙植被，探究其根內真菌群落結構和功能群特征，對植被恢復和森林撫育管理具有重要意義。

樟子松(Pinussylvestrisvar.mongolica)天然分布于呼倫貝爾紅花爾基地區，具有適應性強、耐低溫干旱、生長迅速等特點(任悅等， 2018)，是我國干旱、半干旱和亞濕潤干旱區防風治沙的重要常綠喬木樹種(朱教君等， 2007)。鑒于此，本研究以科爾沁沙地樟子松人工林為研究對象，基于高通量測序和真菌功能群劃分，比較分析不同林齡樟子松根內真菌多樣性、群落組成和真菌營養類型特征，以期為樟子松人工林經營管理提供依據。

1 研究區概況

研究區位于科爾沁沙地東南部的遼寧省彰武縣章古臺沙地森林公園(42°23′N，122°22′E)，海拔226.5 m。屬溫帶亞濕潤大陸性季風氣候，光照充足，晝夜溫差大。年均氣溫8 ℃，年均降水量526.5 mm，年蒸發量約1 932.7 mm，年均日照時數為2 615.2 h。土壤類型為風沙土。研究區常見植物包括樟子松、胡枝子(Lespedezabicolor)、小葉錦雞兒(Caraganamicrophylla)、狗尾草(Setariaviridis)、羊草(Leymuschinensis)、冰草(Agropyroncristatum)等。

2 研究方法

2.1 樣品采集與處理

2017年7月，在26、33和43年生樟子松人工林內各設置1個50 m×50 m的樣地(以下簡稱26a、33a和43a)。對每個樣地內所有樟子松進行樹高和胸徑的測量并進行相關調查，詳情見表1。每個樣地隨機選取間距8 m以上的5株標準木作為采樣對象。采樣時，先用鐵鍬除去林木基部半徑2 m內的地表覆蓋物，后沿其根部小心挖至根的末級，剔除雜根。將3個樣地共15株樟子松細根樣品及少量新鮮土壤(保鮮根樣)分別裝入塑封袋中并編號，并置于4 ℃便攜式保溫箱內帶回。同時，采集適量根系附近土壤樣品用于土壤理化性質分析。

2.2 根內真菌的分子鑒定

根系樣品帶回后保存于-80 ℃冰箱。試驗時，用蒸餾水輕輕沖洗根樣以去除雜物，并剪成10 cm左右的根段，洗凈后將同一樣地的根尖樣品充分混合，再隨機分成5份作為重復組。參照PowerSoil? DNA Kit(MO BIO公司)試劑盒說明書對根尖樣品進行DNA提取。采用PCR對根內真菌的rDNA ITS區進行擴增，真菌通用引物為 ITS1F(5’-CTTGGTCATTTAGAGGAAGTAA-3’)和 ITS2(5’-GCTGCGTTCTTCATCGATGC-3’)，每個樣本重復3次。對PCR擴增產物進行凝膠電泳檢測，合格后使用Illumina MiSeq平臺進行上機測序。利用NCBI數據庫中BLAST 分析對相似性≥97% OTU序列進行比對，并注釋其群落的物種信息，注釋遵循以對比相似度≥97%鑒定到種水平，相似度介于90%～97%之間則可鑒定到屬水平(Kennedyetal.， 2010)。

表1 樣地概況Tab.1 Basic characteristics of the three sampling plots

2.3 真菌功能類群分析

利用FUNGuild v1.0將樟子松根內真菌劃分為不同的營養型和功能群(表2)，結果數據僅保留置信度為“很可能”(Probable)和“極可能”(Highly probable)的OTUs及其類別； 將營養類型鑒定為復合營養型的群落全部納入“其他真菌”(Other fungi)，功能群鑒定為復合多功能方式的群落分別統一為該營養型下的“其他病理/腐生真菌”(Other pathogenic/saprotrophic fungi)(Nguyenetal., 2016)。

表2 真菌營養型和功能群的劃分Tab.2 Classification of trophic mode and guild for fungi

2.4 土壤理化性質測定

采用環刀烘干稱重法測定土壤孔隙度(total soil porosity,TSP)和土壤含水量(soil water content,SWC)，采用pH分析儀(上海雷磁，中國)測定土壤pH值，采用重鉻酸鉀外加熱法測定土壤有機碳(soil organic carbon,SOC)含量，采用全自動化學分析儀(Smartchem 200，意大利)測定土壤全氮(total nitrogen,TN)和全磷(total phosphorus,TP)含量。

2.5 數據分析

根據真菌類群在整體群落中的占比劃分優勢度，占比≥10%為優勢屬，1%～10%為常見屬，≤1%為稀有屬(譙利軍等， 2018)。采用SPSS 20進行數據處理與分析，包括最小顯著性差異法(LSD)和斯皮爾曼相關性檢驗(Spearman’s rank correlation)，運用R 3.5.1進行Shannon、Simpson和Pielou多樣性指數及物種相關性矩陣計算，采用OriginPro 2017和Cytoscape繪圖和輸出結果。

3 結果與分析

3.1 樟子松根內真菌alpha多樣性

從科爾沁沙地26a、33a和43a樣地樟子松根尖樣品中共獲得根內真菌高質量ITS序列384 016條及可劃分OTU數量832個。如表3所示，隨著林齡的增加，樟子松根內真菌Shannon、Simpson和Pielou指數逐漸增大，均在43a樣地中達到最大值，分別是2.04、0.72和0.44，但不同林齡林分之間無顯著差異(P>0.05)。

表3 科爾沁沙地樟子松根內真菌多樣性指數①Tab.3 Diversity indices of fungi in roots of P. sylvestris var. mongolica in the Horqin Sandy Land

3.2 樟子松根內真菌群落組成

科爾沁沙地樟子松根內真菌隸屬于5門16綱54目84科182屬。從門水平上看(圖1)，不同林齡均包含子囊菌門(Ascomycota)、擔子菌門(Basidiomycota)、接合菌門(Zygomycota)、壺菌門(Chytridiomycota)和球囊菌門(Glomeromycota)，且每種菌門在不同林齡中的相對豐度均不存在顯著差異(P>0.05)。子囊菌門的相對豐度隨林齡增加逐漸降低，由68.11%下降至57.79%。擔子菌門相對豐度隨林齡增加而增加，變化范圍為29.32%～38.41%。接合菌門相對豐度均值為2.99%，壺菌門和球囊菌門占所有根內真菌的0.05%。

從屬水平上看(圖2)，各林齡沙地樟子松根內真菌優勢屬和常見屬共包括22屬。其中，26a樣地優勢屬為鞘孢屬(Chalara)和肉座菌屬(Hypocrea)，33a樣地優勢屬為Phialocephala、鞘孢屬和肉座菌屬，43a樣地優勢屬為Phialocephal。26a、33a和43a樣地根內真菌優勢屬占所有真菌屬的比例分別為42.89%、47.58%和17.83%。26a、33a和43a樣地根內真菌常見屬分別有10、9和18屬，占比分別為 46.19%、42.13%和72.34%。稀有屬個數和占比在不同林齡間相近，分別是128、130和126屬以及10.92%、10.30%和9.83%。隨著林齡增長，優勢屬比例呈下降趨勢，常見屬呈增加趨勢，而稀有屬則保持在一個較為穩定的水平。結合它們的相對豐度進行聚類分析(圖2)，26a和33a樣地聚合性較好，表明其2組根內真菌組成結構具有較高的相似度。

圖2 科爾沁沙地樟子松根內真菌聚類熱圖Fig.2 Cluster heatmap of fungi in roots of P. sylvestris var. mongolica in the Horqin Sandy Land

圖1 科爾沁沙地樟子松根內真菌門的相對豐度Fig.1 The relative abundance of fungi in roots of P. sylvestris var. mongolica in the Horqin Sandy Land

3.3 沙地樟子松根內真菌生態功能群

科爾沁沙地樟子松不同營養型真菌相對豐度隨著林齡的增加呈現出不同的變化趨勢(圖3)。病理營養型真菌相對豐度呈現下降趨勢，由17.48%下降到7.73%，腐生營養型較為穩定(7.84%～10.06%)，共生營養型則表現出較大波動，在26a、33a和43a樣地中的相對豐度分別為22.05%、38.90%和48.79%。其中，不同林齡中共生營養型根內真菌主要由外生菌根真菌和內生真菌組成，且均存在少量叢枝菌根真菌。26a和33a樣地的共有優勢屬(鞘孢屬和肉座菌屬)均為復合營養型真菌，33a和 43a樣地的共有優勢屬Phialocephala為內生真菌。43a和33a樣地中相對豐度較高的乳菇屬(Lactarius)、絲蓋傘屬(Inocybe)、糙緣腺革菌屬(Amphinema)和須腹菌屬(Rhizopogon)均為共生營養型中的外生菌根真菌。

圖3 科爾沁沙地樟子松根內真菌功能群的網絡分布Fig.3 Function network of fungi in roots of P. sylvestrisvar. Mongolica in the Horqin Sandy Land不同真菌營養型在各樣地中的相對豐度用扇形圖表示,真菌屬在不同林齡的相對豐度由邊厚度表示,節點的面積表示所有樣品中屬的絕對豐度,不同林齡根內真菌共16個優勢屬和常見屬被可靠鑒定。The pie diagram indicates the relative abundance of trophic mode in different age groups, the thickness of links indicates the relative abundance of genera, the square measure of nodes roughly represents the absolute abundance of genera in all samples, 16 dominant genera and common genera were identified reliably.

3.4 土壤理化性質對根內真菌多樣性和營養型的影響

不同林齡樟子松根圍土壤的理化性質如圖4所示。 26a樣地中土壤含水量顯著高于33a和43a樣地(P<0.05)，土壤有機碳、全氮和全磷含量均呈隨林齡增加而均出現顯著升高(P<0.05)，表明土壤養分含量上升。根系周圍土壤理化性質與根內真菌多樣性指數和根內真菌營養型相對豐度的相關性如表4所示，根內真菌多樣性與土壤理化性質均無顯著相關關系(P>0.05)，不同真菌營養型中僅共生營養型真菌相對豐度與土壤含水量顯著負相關(P<0.05)。

圖4 科爾沁沙地樟子松根圍土壤理化性質Fig.4 Soil properties of rhizospheric soil under P. sylvestris var. mongolica in the Horqin Sandy LandSWC： 土壤含水量Soil water content; TSP： 土壤孔隙度Total soil porosity; SOC： 有機碳Soil organic carbon; TN： 全氮Total nitrogen; TP： 全磷Total phosphorus. 下同The same below.

表4 沙地樟子松根內真菌多樣性和營養型相對豐度與土壤理化性質的相關性①Tab.4 Correlation coefficients between diversity and relative abundance of the trophic mode with soil properties

4 討論

4.1 科爾沁沙地樟子松根內真菌多樣性和群落組成

土壤真菌多樣性指數的高低代表著群落穩定性和抵抗力強弱(高尚坤等， 2018)，43a樣地中樟子松根內真菌群落抗擾能力強于26a和33a樣地，這對樟子松的健康生長具有積極效應。眾所周知，植物根內真菌豐富度和多樣性常受土壤pH值和有機質含量等因素影響(Toljanderetal., 2006； Barnesetal.， 2016b)，然而沙地樟子松根內真菌多樣性與環境干擾因素的相關性較差，這種現象在先前研究中同樣存在(Jumpponenetal., 2014； Entwistleetal., 2013)，這是因為真菌多樣性與環境因子的相關性取決于所處生態系統的范圍和類型，在某些生態系統中真菌對氣候變化的承受力具有較大彈性(Zhengetal., 2017)。

從門水平上看，科爾沁沙地樟子松根內真菌主要由子囊菌門和擔子菌門組成，其余為少量結合菌門及極少量壺菌門和球囊菌門，這與櫟樹(Quercus)和極地柳(Salixpolaris)等外生菌根依賴型樹種根內真菌群落組成結果相似(Tojuetal., 2013b； Kernaghanetal., 2011)。其中，不同林齡中子囊菌門相對豐度均高于擔子菌門，一是因為子囊菌的進化速度和物種多樣性均高于擔子菌(王海英等, 2010)，在干旱地區的適應性更強； 二是因為沙地土壤中養分含量較低，而子囊菌門包含的許多腐生真菌能將難降解性有機質分解、促進營養元素礦化過程，使土壤中養分含量增加并提高可利用率，從而間接影響著植物的生長(Barnesetal., 2018)。從屬水平上看，沙地樟子松根內真菌優勢屬為鞘孢屬、肉座菌屬和Phialocephala，其中： 鞘孢屬被認為是造成植物枯梢的植物病原菌(Tuliketal., 2018)，肉座菌屬是通過產生抗生素抑制植物病原真菌的腐生真菌(Jaklitsch, 2009)，而FUNGuild卻將鞘孢屬鑒定為復合營養型，肉座菌屬鑒定為病理營養型，這可能是因為內生真菌在發育階段的不穩定性導致它們不停地改變生態策略，于是生活史復雜的病原菌在不同的宿主中表現出不同的生活狀態(Arnoldetal., 2007)； 也可能是伴隨著根內真菌所在植物根系組織的死亡，其為維持正常的生理活動而改變自身的營養方式(Paulin-Mahadyetal., 2002)。Phialocephala為常出現于健康植物中的深色有隔內生真菌(dark septate endophytes,DSE)，在針葉樹種DSE中占主要地位(Suronoetal., 2018)，它能激發植物對非生物脅迫的抗性，加速根系內物質運輸和礦化作用，抑制根系病原菌(Mandyametal., 2005； Tellenbachetal., 2013)。總體上看，不同林齡沙地樟子松根內真菌優勢屬和常見屬的種類和占比波動較大，出現真菌組成異質性的原因可能是林齡影響了宿主植物根系分泌物的組成成分(李嬌等， 2014)，而不同根系分泌物將直接影響根際微生物的數量和群落結構(Lladóetal., 2018)。

4.2 科爾沁沙地樟子松根內真菌生態功能特征

植物根系中不同營養型真菌之間相互促進與制約，根內真菌功能群結構組成與植物生長發育狀況密不可分。43a和33a樣地中，相對豐度高于26a樣地的根內真菌為乳菇屬、絲蓋傘屬、糙緣腺革菌屬和須腹菌屬，它們均為樟子松林常見的外生菌根真菌(郭米山等, 2018)，這證實了外生菌根真菌多樣性在森林發育的前30～40年增長較快(Wuetal., 2013； Wallanderetal., 2010)。共生營養型真菌的增加同時對病理營養型和腐生營養型真菌產生影響： 一方面，共生營養型真菌可對植物致病菌產生拮抗作用，減少病原菌的定殖以保護植物免受病害(Nagyetal.， 2013； Frewaetal., 2018)，導致病理營養型真菌不斷減少； 另一方面，菌根真菌(尤其是外生菌根真菌)對土壤中營養物質的大量攝取，腐生真菌和其他微生物分解者受到抑制(Grauetal., 2017)。因此，即使土壤養分含量隨林齡增加，科爾沁沙地樟子松根內腐生營養型真菌相對豐度一直處于較為平穩的水平。沙地樟子松根內真菌營養型結構隨林齡的變化表現出一定的波動，共生營養型真菌占優勢時宿主的抗擾能力增加(Chuetal., 2016)，反之，若病理營養型和腐生營養型真菌比例高于共生營養型真菌，則極易導致植物病害發生(Millbergetal., 2015)。

土壤理化性質和根內真菌功能類型組成之間的因果關系難以確定，真菌自身生理活動直接影響土壤生態系統的健康與穩定，同時土壤水分和養分又是真菌的主要物質來源(Lietal., 2018)。干旱區沙地樟子松的生長發育過程需要大量水分，增加共生真菌有助于植物根系從土壤中吸收水分(Groveretal., 2011)，從而導致土壤含水率降低。土壤pH值在土壤真菌群落組成和功能類型上的影響普遍存在(Barnesetal., 2018)，研究區沙地樟子松人工林根圍土壤pH值在6～7之間，并與根內真菌功能群結構無顯著相關性，這是因為真菌對pH值變化的敏感程度較低，具有較廣的土壤酸堿耐受范圍，且無明顯生長抑制區(許美玲等， 2004； Rousketal., 2010)。本研究中大部分土壤理化性質并未與根內真菌營養型組成出現顯著相關關系，這是因為宿主根系與根內真菌保持著相對穩定的供養和反饋相互作用(陳智裕等， 2017)，此時根系微環境已經過強烈的環境過濾，改變土壤性質對根內真菌群落的變化影響程度較低(Becketal., 2015)。

值得注意的是，每個樣地的樟子松根尖樣品均鑒定出極少量的叢枝菌根真菌，同樣作為具有外生菌根真菌侵染偏好型的北美紅櫟(Quercusrubra)和赤松(Pinusdensiflora)根系中也發現了此類現象(Dickieetal., 2001； Tojuetal.， 2018a)。在26a和43a樣地根尖樣品中鑒定出的杜鵑花類菌根真菌粉枝孢屬(Oidiodendron)，在紙樺(Betulapapyrifera)、香脂冷杉(Abiesbalsamea)和白云杉(Piceaglauca)等非杜鵑花植物根系中也被鑒定出(Kernaghanetal., 2011)。這跨越了根據宿主植物種類對菌根真菌進行劃分的傳統分類模式，更進一步證實宿主與真菌共生體存在較強的環境依賴性(Tojuetal., 2018b)。

5 結論

1)科爾沁沙地樟子松根尖樣品共獲得832個OTUs，根內真菌多樣性指數隨林齡的增加而不斷增大，但不同林齡樟子松人工林不存在顯著差異，且根內真菌多樣性和土壤理化性質無顯著相關關系。

2)不同林齡樟子松根內真菌隸屬于5門16綱54目84科165屬，子囊菌門和擔子菌門為主要組成部分，根內真菌優勢屬包括鞘孢屬、肉座菌屬和Phialocephala。不同林齡根內真菌群落組成存在異質性，中優勢屬和常見屬菌種類型和相對豐度變動較大。

3)科爾沁沙地樟子松根內真菌功能群組成豐富，隨林齡的增加科爾沁沙地樟子松根內共生營養型真菌急劇增加，同時病理營養型呈現下降趨勢，腐生營養型維持在一個穩定水平。土壤含水量和共生營養型真菌相對豐度存在顯著相關關系。