树种身份经由根代谢物介导反馈与亚热带森林根际微生物组组装相关

《Functional Ecology》：Tree species identity is associated with rhizosphere microbiome assembly via root metabolite-mediated feedbacks in subtropical forests

　　

根际微生物组（rhizosphere microbiome）对森林生态功能至关重要，但树种身份（tree species identity）和邻近树木相互作用如何与其组装相关联的机制仍知之甚少。
研究人员利用 Illumina MiSeq 和液相色谱–质谱联用（LC–MS），研究了亚热带森林中树种（松树 vs. 枫香）和森林组成（纯林 vs. 混交林）如何塑造根际微生物群落和根系代谢物。
根际真菌群落结构在不同根际类型间存在显著差异，而细菌群落和代谢物谱主要在纯松林与其他森林类型之间不同。预测的功能变化包括纤维素分解细菌下降（?27%）和反硝化类群增加（+54%），混交林中外生菌根真菌较纯林富集（p < 0.05）。
代谢物谱分析显示，脂质、有机酸和糖类是最丰富的化合物类别，纯松林与混交林根际之间有55种代谢物差异丰度。根际代谢物谱与微生物群落结构显著相关（p < 0.05）。研究人员鉴定了特定的微生物–代谢物关联；脂肪酸与寡营养细菌（WPS-2）正相关，糖类与反硝化菌（Reyranella）相关，二肽与有益真菌 Penicillium 相关。
这些结果表明，树种身份通过根代谢物相关反馈与根际微生物组组装相关联，而这些反馈与土壤养分循环相关。这为理解植物–微生物组相互作用并利用其增强生态系统韧性提供了机制性、功能导向的框架。

研究背景与问题
根际微生物组（rhizosphere microbiome，植物根系周围狭窄土壤区域中的微生物集合）是森林生态系统功能的重要组成部分，调控从养分循环到植物健康的过程。树种身份（tree species identity）和森林组成可通过根系分泌物、凋落物质量和土壤性质差异影响根际微生物群落。然而，特定化学机制，尤其根系代谢物在介导植物–微生物相互作用中的作用，在森林生态系统中仍知之甚少，特别是在经历组成变化的森林中。根际是植物根系与土壤和微生物相互作用的关键界面，其微生物驱动生物地球化学转化，影响养分有效性和植物生长健康；包含有益共生体和潜在病原体。根系分泌物被认为是农业系统微生物群落组装的关键驱动因素，但在森林根际的功能重要性关注较少。代谢组学可直接捕获根际化学景观，检测有机酸、糖、次生代谢物等实际小分子，这些可作为营养、信号或抗菌剂，因而适合揭示微生物招募和群落组装的化学机制。将代谢组谱与微生物群落测序整合，可识别特定代谢物–微生物关联，并生成关于森林根际微生物组结构化学驱动因素的可检验假说。本研究检查四种根际类型：纯松林松树根际（PRP）、混交松–阔叶林松树根际（PRM）、混交林枫香根际（LRM）、纯枫香林枫香根际（LRP）。研究目的为：确定树种身份对根际微生物群落的影响；研究这些森林类型相关的根际代谢物谱变化；分析根际微生物和代谢物之间潜在关联。假设为：树种身份和森林组成与不同根际微生物群落组成相关；纯林和混交林代谢物谱不同，反映宿主特异性化学特征；特定根际代谢物与特定微生物类群相关，提示微生物群落组装的潜在化学关联。研究人员利用 Illumina MiSeq 和 LC–MS，发现树种身份通过根代谢物相关反馈与根际微生物组组装相关，并与土壤养分循环相关；提供机制性、功能导向框架，以理解植物–微生物组互作并利用其增强生态系统韧性。

关键技术方法
研究在南京紫金山（Purple (Zijin) Mountain Park）亚热带季风区开展。选择纯松林（Pinus thunbergii Parl.）、松/枫香混交林（P. thunbergii + Liquidambar formosana Hance）和纯阔叶林（L. formosana）。2021年4月9日，各森林类型设3个50×80 m样地，纯林选优势树，混交林选松及邻近枫香；根际土为距根表2 mm内土壤，四类各9个，共36个。DNA提取后扩增细菌16S rRNA V3–V4区和真菌ITS1区，Illumina MiSeq测序，QIIME2、DADA2、SILVA 138.1和UNITE v.99_2022.11注释，FAPROTAX和FUNGuild预测功能。每类随机6个共24个样品进行LC–MS非靶向代谢组，Progenesis QI处理，数据库注释及KEGG映射。统计用R，包括PCA、PERMANOVA、ALDEx2、指示种分析、OPLS-DA、Procrustes、MMvec和PLS-PM。样本队列来源：南京紫金山。

研究结果
3.1 根际土壤微生物多样性和群落组成
通过Illumina MiSeq测序和α多样性分析，研究人员发现细菌丰富度（observed ASVs）和Shannon指数在森林类型间无差异，但Pielou均匀度指数在混交林（PRM和LRM）高于纯林（PRP和LRP）（Dunn's test，p=0.001）；纯枫香林（LRP）真菌多样性和均匀度高于混交林（LRM）（p=0.003和p=0.002）。PCA和PERMANOVA显示森林类型显著影响细菌和真菌群落结构（p=0.001）。成对检验表明，纯松林根际细菌群落结构与混交林和纯阔叶林显著不同；真菌群落结构在所有三类森林间均不同；混交林内松（PRM）和枫香（LRM）根际细菌和真菌群落无差异。测序获得1,532,141条细菌序列（5573个ASVs），代表26门，Proteobacteria最丰富，其次是Acidobacteriota和Actinobacteriota；属水平以Acidothermus、Candidatus_Solibacter和Bradyrhizobium优势。差异丰度分析发现42个细菌属差异；PRP中Acidothermus、Subgroup_2、Mycobacterium、Roseiarcus、Occallatibacter、KF-JG30-B3和WPS-2富集，混交林中Reyranella和TK10较多；LRP中Bradyrhizobium和Mycobacterium富集。真菌获得4,209,004条序列（2786个ASVs），分属12门，Basidiomycota和Ascomycota最丰富，优势属包括Archaeorhizomyces、Cladophialophora和Russula。20个真菌属差异；PRP中Penicillium、Podila、Geminibasidium和Saitozyma富集，混交林中Clavulina较高；LRP中Saitozyma和Solicoccozyma高于混交林；混交林中Clavulina、Lactarius和Geminibasidium较LRP更普遍。

3.2 微生物群落功能谱
通过FAPROTAX和FUNGuild预测，细菌和真菌功能谱在森林类型间显著不同（PERMANOVA，p<0.05）。PRP细菌功能谱与混交林和纯阔叶林不同；真菌功能谱在纯林与混交林间不同，而混交林内松和枫香相似。FAPROTAX注释16.7% ASVs，主要功能为chemoheterotrophy（33.98%）、aerobic chemoheterotrophy（33.76%）和cellulolysis（13.47%）。PRP中chemoheterotrophy、cellulolysis、nitrogen fixation和ureolysis预测富集；混交林中denitrification相关功能更丰富。cellulolysis从PRP的11.6%降至PRM的8.43%（?27%），denitrification从0.57%升至0.88%（+54%）。FUNGuild注释27.8% ASVs，外生菌根真菌（ectomycorrhizal fungi）在混交林多于两种纯林，undefined saprotrophs在PRP富集，animal pathogens在LRP富集。原文强调这些为生物信息学推断，应视为假设生成而非确认功能变化。

3.3 根际代谢物谱
LC–MS非靶向代谢组检测到7149种代谢物，分17类；脂质和类脂分子最多（2068），其次有机酸及衍生物（556）、organoheterocyclic compounds（475）、organic oxygen compounds（427）、benzenoids（417）、phenylpropanoids and polyketides（369），2647种未分类。PCA显示PRP代谢物谱与混交林和纯阔叶林明显分离（PRP vs PRM：p=0.006；PRP vs LRP：p=0.006）；混交林内松和枫香，以及LRP与LRM之间代谢谱无显著差异。PRP与PRM之间有55种差异代谢物（VIP≥2，p<0.01），LRP与LRM之间3种，PRM与LRM无。PRP富集脂质（如eicosanoids、diterpenoids、fatty acids）、有机酸（包括氨基酸）、benzenoids和organosulfur compounds；organic oxygen compounds、organoheterocyclic compounds和phenylpropanoids在PRP和PRM中均存在；PRM富集coumarins和flavonoids等特定phenylpropanoids，而PRP中stilbenes更多。LRP与LRM的3种差异代谢物为脂质，属fatty alcohols和monoradylglycerols，在LRP中更丰富。

3.4 根际微生物组与代谢物谱的相关性
Procrustes分析显示根际代谢物谱与细菌群落结构（M²=0.401，p=0.015）和真菌群落结构（M²=0.404，p=0.020）显著关联。MMvec分析进一步识别特定微生物–代谢物共现模式：脂肪酸（如9,12,15-Eicosatriynoic acid）与寡营养细菌门WPS-2（logCP=2.31）和Occallatibacter（logCP=1.89）正共现；糖stachyose与反硝化属Reyranella（logCP=2.05）和TK10（logCP=1.76）共现；PRP富集真菌Saitozyma、Penicillium、Podila和Geminibasidium与fatty acid（8,11,14,18-Eicosatetraynoic acid）共现，Penicillium、Podila和Geminibasidium还与dipeptide（Phenylalanyl-Isoleucine）共现。PLS-PM解释根际类型间细菌和真菌群落结构52.7%的方差；森林类型与代谢物谱关联最强（路径系数=?0.679）；代谢物对细菌群落结构的正关联（0.682）大于对真菌群落结构的关联（0.538）。原文说明MMvec产生的是log条件概率而非标准相关系数，PLS-PM和MMvec为相关性分析。

讨论部分总结
讨论指出，森林类型与根际微生物组结构差异显著相关，但可能部分由与森林组成共变的土壤理化性质介导；采样设计不能完全区分树种身份效应与松材线虫病造成的历史背景。细菌群落中PRP与其他森林类型不同，而真菌群落在三类森林间均不同，提示真菌对树种组成变化更敏感；PRP与LRP分离确认树种身份是根际定殖的关键决定因素，但土壤pH、有机质、氮和磷等随森林类型变化且未在根际测定，不能完全分解方差。与假设相反，混交林中松和枫香根际无差异，提示局部土壤条件或共享根系分泌物可能均质化微生物栖息地，与 edaphic factors 可覆盖植物身份的结构作用一致。PRP富集耐酸细菌Roseiarcus、Occallatibacter和候选门WPS-2，与松林较低pH和较低养分有效性一致；WPS-2寡营养特性与PRP低养分一致。混交林Reyranella反硝化菌增加与反硝化功能富集一致。真菌中Penicillium解磷菌在PRP较多，Saitozyma纤维素分解属在纯林富集，可能反映凋落物质量和碳循环差异；病原拮抗菌Solicoccozyma在LRP较多。细菌丰富度和多样性无差异但纯林均匀度较低，机制未明；一种可能是纯松林中萜类和脂肪酸等代谢物有利于某些优势类群而抑制其他类群，混交林资源异质性更大，降低竞争排斥。LRP真菌多样性较高与脂质如fatty alcohols较高一致，可能刺激真菌生长和多样性。根际代谢物组成在松树和阔叶林间显著不同，强调树种塑造根际化学；混交林中松与枫香代谢物趋同，原文提出资源均质化、化感作用和共享微生物代谢根分泌物三个非互斥且未经检验的假设。微生物–代谢物关联中，WPS-2与脂肪酸9,12,15-Eicosatriynoic acid正相关，WPS-2作为好氧有机异养菌可清除有机氮化合物，可能在营养限制下利用该代谢物；Reyranella与stachyose相关，stachyose为已知趋化剂，可能反映混交林根际氮循环增强。PLS-PM表明森林类型强关联代谢物谱，代谢物对细菌群落结构影响大于真菌；差异可能因细菌生长更快、世代更短，对代谢物可用性响应更快，更依赖不稳定碳源，且有机酸可改变土壤pH和微环境选择细菌；真菌代谢多样性更高，能降解复杂聚合物，并可通过菌丝网络获取根际外养分，缓冲局部化学变化。但原文强调这些共现必须谨慎解释：PLS-PM和MMvec固有相关性，不能建立因果；共现可能来自直接代谢物效应、间接经其他类群、共享环境偏好或植物介导效应；需要代谢物添加或根分泌物操控等实验验证。预测功能模式与观察到的群落和环境变化一致：PRP纤维素分解功能富集与松树难降解、富纤维素凋落物一致；混交林和阔叶林反硝化增加与较高凋落物输入、土壤碳和氮升高及厌氧微位点相关；混交林外生菌根真菌富集可能反映养分获取和种间植物相互作用增强。但FAPROTAX和FUNGuild依赖分类注释，土壤环境数据库覆盖有限，预测代表潜在生态角色而非确认活动，需直接酶活或宏转录组验证。

研究结论部分翻译
总之，本研究表明，树种身份和森林组成与不同的根际微生物群落结构和代谢物谱相关。真菌对树种组成变化更敏感，而细菌和代谢物变化在纯松林与其他森林类型之间最为明显。预测的功能变化提示纤维素分解下降和反硝化增加，同时混交林中外生菌根真菌富集。在混交林内，松树和枫香根际在微生物群落和代谢物谱上表现趋同，提示邻近树木效应，如共享根系分泌物或均质化土壤条件，可能覆盖宿主特异性信号。代谢物谱与微生物群落结构强相关，表明根系分泌物可能影响微生物组装。然而，研究人员承认，研究设计不能完全区分树种身份效应与病害遗留效应，因为森林组成梯度由松材线虫病引发。此外，研究结果为相关性，未建立因果关系。这些发现强调了植物–代谢物–微生物关联在跨森林类型地下生态变化中的重要性，同时强调未来需要合适对照的研究以建立因果关系并区分物种身份与干扰遗留机制。