《Functional Ecology》：Below-ground drivers of rhizosphere functional gene assembly during Moso bamboo expansion
編輯推薦：
本文通過整合個(gè)體水平根系性狀測量、根際土壤分析及宏基因組測序，系統(tǒng)研究了亞熱帶毛竹在純竹林（PBF）與混交竹林（MBF）擴(kuò)張?zhí)荻认碌牡叵律鷳B(tài)過程。研究發(fā)現(xiàn)，混交林中的竹子呈現(xiàn)出更具“獲取型”的根系策略（如更高的比根長SRL和比根面積SRA），而純竹林根際則富集了多種氮（N）、磷（P）循環(huán)功能基因，并形成了連接性更強(qiáng)的微生物共現(xiàn)網(wǎng)絡(luò)。結(jié)構(gòu)方程模型表明，森林類型主要通過改變土壤性質(zhì)（而非根系性狀的直接作用）來塑造功能基因組裝。這些結(jié)果揭示了已建立的毛竹林通過土壤修飾產(chǎn)生了自我強(qiáng)化的微生物反饋，形成了功能整合、養(yǎng)分循環(huán)潛力增強(qiáng)的根際群落，為在擴(kuò)張前沿（即反饋尚未鞏固的混交林階段）進(jìn)行管理干預(yù)提供了關(guān)鍵窗口。
1 引言
毛竹（Phyllostachys edulis）擴(kuò)張是亞熱帶亞洲最顯著的本土植物入侵之一，在中國每年擴(kuò)張速率約為2%，導(dǎo)致生物多樣性顯著喪失并改變生態(tài)系統(tǒng)功能。傳統(tǒng)擴(kuò)張機(jī)制研究多關(guān)注資源競爭、化感作用等地上過程，而根-土-微生物相互作用這一關(guān)鍵的地下驅(qū)動(dòng)過程仍不清楚。特別是，種內(nèi)根系性狀變異如何影響根際微生物功能潛力，以及根系性狀與土壤性質(zhì)在功能基因組裝中的相對(duì)貢獻(xiàn)，是當(dāng)前的研究空白。本研究旨在闡明毛竹擴(kuò)張過程中，根系性狀變異與根際土壤性質(zhì)如何共同塑造微生物功能基因譜。我們整合了根際生態(tài)位構(gòu)建理論和根系經(jīng)濟(jì)學(xué)譜系概念，假設(shè)：（1）混交林（MBF）中的竹子相比純竹林（PBF）會(huì)表現(xiàn)出更“獲取型”的根系策略；（2）純竹林根際擁有功能更整合的微生物群落，其養(yǎng)分循環(huán)基因更富集，共現(xiàn)網(wǎng)絡(luò)更復(fù)雜；（3）根系形態(tài)性狀、化學(xué)組成和根際性質(zhì)通過交互作用共同驅(qū)動(dòng)微生物功能基因組成。
2 材料與方法
研究在江西井岡山國家級(jí)自然保護(hù)區(qū)進(jìn)行，涵蓋4個(gè)海拔梯度。我們比較了純竹林（PBF，毛竹蓋度>80%）和混交竹林（MBF，毛竹與鄉(xiāng)土樹種共優(yōu)），每個(gè)林型在每個(gè)海拔點(diǎn)選取6株標(biāo)準(zhǔn)竹進(jìn)行采樣。測定了竹子細(xì)根的形態(tài)性狀（如比根長SRL、比根面積SRA、根組織密度RTD）、化學(xué)性狀（如N、P、K、木質(zhì)素、纖維素等）以及根際土壤的理化性質(zhì)（如pH、SOC、TN、TP、MBC、MBN等）。通過鳥槍法宏基因組測序分析了微生物功能基因（基于KEGG、CAZy、NCycDB、PCycDB數(shù)據(jù)庫），并構(gòu)建了功能基因共現(xiàn)網(wǎng)絡(luò)。使用線性混合模型、PERMANOVA、變異分配分析（VPA）、結(jié)構(gòu)方程模型（SEM）等統(tǒng)計(jì)方法分析數(shù)據(jù)。
3 結(jié)果
3.1 林型間根系性狀與根際土壤性質(zhì)的分化
混交林（MBF）中的竹子表現(xiàn)出顯著的“獲取型”根系策略：其比根長（615.33 vs. 446.03 cm/g）和比根面積（156.06 vs. 122.86 cm2/g）顯著高于純竹林（PBF），而根組織密度（0.32 vs. 0.39 g/cm3）和根直徑顯著更低。同時(shí)，MBF中竹子根系的磷、鉀和淀粉含量也更高。在土壤性質(zhì)方面，PBF根際土壤pH值更高（5.53 vs. 5.28），但電導(dǎo)率更低；而MBF根際土壤的有效磷、速效鉀、土壤有機(jī)碳（SOC）、微生物生物量碳（MBC）和氮（MBN）則顯著更高。主成分分析顯示，根際土壤性質(zhì)在林型間的分化程度強(qiáng)于根系性狀。
3.2 毛竹根際功能基因組成與多樣性
非度量多維尺度分析表明，林型顯著區(qū)分了所有代謝本體數(shù)據(jù)庫（KEGG、CAZy、NCycDB、PCycDB）的功能基因組成，但海拔梯度解釋了更大的變異。功能基因α多樣性（如基因豐富度、香農(nóng)指數(shù)）在林型間大多無顯著差異，但普氏分析證實(shí)功能基因組成與微生物群落分類結(jié)構(gòu)（尤其是細(xì)菌群落）高度一致。
3.3 功能基因豐度與林型響應(yīng)基因
Z-score分析揭示了功能基因豐度的系統(tǒng)變化。在碳水化合物活性酶（CAZy）中，僅糖基轉(zhuǎn)移酶在PBF中豐度顯著更高。在氮循環(huán)基因（NCycDB）中，PBF根際在有機(jī)氮降解與合成、同化硝酸鹽還原、反硝化、異化硝酸鹽還原、硝化和固氮共六條通路中均顯示出顯著更高的基因豐度。在磷循環(huán)基因（PCycDB）中，PBF根際在嘌呤代謝、轉(zhuǎn)運(yùn)蛋白、嘧啶代謝、膦酸酯代謝、有機(jī)磷酸酯水解等七類功能上顯著富集。指示物種分析進(jìn)一步發(fā)現(xiàn)，PBF擁有更多與養(yǎng)分循環(huán)相關(guān)的指示基因（如反硝化基因nosZ, narG，有機(jī)磷礦化基因phoD），而MBF的指示基因則更多與碳水化合物（如半纖維素）降解相關(guān)。
3.4 林型特異的功能基因共現(xiàn)網(wǎng)絡(luò)
共現(xiàn)網(wǎng)絡(luò)分析顯示，PBF根際的微生物功能基因網(wǎng)絡(luò)通常具有更高的連接性、模塊化和更多的正相互作用。尤其是在磷循環(huán)網(wǎng)絡(luò)中，PBF的網(wǎng)絡(luò)邊數(shù)幾乎是MBF的兩倍（1520 vs. 788），網(wǎng)絡(luò)密度和聚類系數(shù)也更高，表明其微生物群落內(nèi)部代謝協(xié)作更緊密、結(jié)構(gòu)更穩(wěn)定。
3.5 根系性狀與土壤性質(zhì)共同驅(qū)動(dòng)微生物功能基因組成
距離冗余分析和變異分配分析表明，根際土壤性質(zhì)和根系性狀共同解釋了功能基因組成的大部分變異，其中土壤性質(zhì)的獨(dú)立貢獻(xiàn)（23.43%-28.64%）大于根系性狀（2.31%-8.14%），但兩者存在顯著的協(xié)同效應(yīng)（共同解釋30.29%-38.99%）。土壤有機(jī)碳（SOC）、全磷（TP）和速效鉀（AK）是驅(qū)動(dòng)功能基因組成的最主要環(huán)境因子。曼特爾檢驗(yàn)顯示，比根面積（SRA）和比根長（SRL）與所有數(shù)據(jù)庫的功能基因組成均呈最強(qiáng)正相關(guān)。結(jié)構(gòu)方程模型最終揭示，林型通過顯著改變土壤性質(zhì)來直接影響功能基因組裝，而根系性狀的直接效應(yīng)有限。
4 討論與結(jié)論
本研究表明，毛竹擴(kuò)張伴隨著顯著的地下生態(tài)位構(gòu)建過程。混交林中的竹子為應(yīng)對(duì)種間競爭壓力，采取了“獲取型”根系策略以最大化資源捕獲，這同時(shí)刺激了根際微生物活性和養(yǎng)分礦化，創(chuàng)造了動(dòng)態(tài)的根際微環(huán)境。而一旦擴(kuò)張成功形成純竹林，竹子則轉(zhuǎn)向更“保守”的策略，并通過長期土壤修飾（如提高pH），塑造出一個(gè)功能高度整合的根際微生物群落。該群落以全面富集的氮、磷循環(huán)基因?yàn)樘卣鳎⑿纬蛇B接緊密、協(xié)作性強(qiáng)的代謝網(wǎng)絡(luò)，從而建立了高效且穩(wěn)定的養(yǎng)分獲取與循環(huán)機(jī)制。這種自我強(qiáng)化的植物-土壤正反饋，是毛竹能夠成功擴(kuò)張并長期占據(jù)生態(tài)位的核心地下機(jī)制。
研究發(fā)現(xiàn)，土壤性質(zhì)（而非根系性狀）是驅(qū)動(dòng)根際功能基因組裝的主要途徑。這意味著，管理毛竹擴(kuò)張的關(guān)鍵窗口期在于其擴(kuò)張前沿（即混交林階段），此時(shí)自我強(qiáng)化的微生物反饋尚未鞏固。對(duì)于已建立的純竹林，僅移除地上部分可能不足以恢復(fù)生態(tài)系統(tǒng)，因?yàn)橥寥肋z留效應(yīng)（如改變的微生物群落和功能）會(huì)持續(xù)存在。未來的恢復(fù)策略可能需要結(jié)合土壤改良和本地微生物群落的引入，以克服這些地下遺留影響，從而有效維護(hù)亞熱帶森林的生物多樣性與生態(tài)系統(tǒng)服務(wù)。
