刺萼龙葵（Solanum rostratum）根分泌物、根际微生物与土壤代谢物的密度依赖性变化及其“分泌物–微生物–代谢物”级联

《Plants》：Density-Dependent Changes in Root Exudates, Rhizosphere Microbes, and Soil Metabolites of Solanum rostratum

　　

种内竞争塑造了入侵植物的高密度单优群落，但密度如何通过根分泌物与根际微生物重组及土壤代谢物积累相关联仍不清楚。研究人员设置了4个种植密度（每盆1、3、5和10株），并整合了生长与光合生理、根分泌物代谢组学、根际宏基因组学和土壤代谢组学，采用Mantel检验和变量水平相关分析。密度增加降低了除株高以外的所有性状，株高表现为轻微的非显著增加后显著下降。根分泌物从低密度下以酚酸/萜类为主转向高密度下以糖/肽类为主。低密度富集的细菌属与生物量呈正相关，而高密度富集的属（如Ilumatobacter、Allosphingosinicella）呈负相关。土壤代谢物如Glutarate Semialdehyde和Sutherlandin D与生物量、根系性状和光合作用呈负相关。Mantel检验显示植物表型与根分泌物（r = 0.80）以及分泌物与土壤代谢物（r = 0.83）之间存在强相关性。变量水平分析将酚酸类分泌物与推定生长相关细菌（如Flavitalea、Blastococcus）联系起来，将糖/肽类分泌物与推定胁迫适应细菌（如Allosphingosinicella、Ilumatobacter）联系起来，后者与Glutarate Semialdehyde和Sutherlandin D积累呈正相关。这些结果表明，密度相关的分泌物变化与根际微生物重组和土壤代谢物积累同步发生，形成“分泌物–微生物–代谢物”级联。该级联可能缓冲高密度种内竞争，并有助于解释入侵范围内的单优群落形成。

研究背景方面，植物密度通过改变资源竞争格局深刻影响植物形态建成、生理性状和根际过程。在入侵生态学中，许多入侵植物常形成高密度单优群落，种内竞争可能是其适应性扩张的重要调控力量。刺萼龙葵（Solanum rostratum）入侵后同样呈现连片高密度生长，但种内竞争如何系统调控其根际过程尚不清楚。已有研究表明种内竞争可改变根分泌物代谢谱，植物密度可重塑根际化学，密度相关的分泌物重塑还可能加剧负向植物–土壤反馈；根际微生物可影响入侵植物的种内竞争，入侵植物竞争力可能受氮利用策略与根际微生物组共同调控。根分泌物是植物–土壤微生物互作的核心化学媒介，植物可通过根分泌物主动招募有益微生物并干扰养分循环。当前入侵植物根际研究多集中于种间竞争，而种内竞争加剧时根分泌物组成是否改变并进一步与根际细菌群落重组和土壤代谢物积累相关联尚缺乏研究。刺萼龙葵为原产北美的恶性一年生入侵植物，已在我国多地迅速扩散，严重威胁农业与生态系统，现有研究多集中于种子性状及化感与自毒效应，而种内竞争如何调控其幼苗生长、生理与根际微生态系统间的多层次互作仍需深入探究。

基于此，研究人员设置4个种植密度处理（每盆1、3、5和10株），系统研究种内竞争对刺萼龙葵幼苗形态、光合生理、根分泌物、根际细菌群落和土壤代谢物的影响，旨在阐明种内竞争下“分泌物–微生物–代谢物”级联的响应模式，为入侵植物的生态调控提供理论基础。研究假设：种内竞争加剧会改变根分泌物组成，与根际细菌群落从推定生长相关类群向推定胁迫适应类群的转变相关联，并伴随与生长下降相关的土壤代谢物积累，形成潜在的根际反馈。

主要关键技术方法方面，研究人员以中国新疆乌鲁木齐水磨沟区石人沟（43°49′ N，87°46′ E，海拔约950 m）未被刺萼龙葵入侵的斑块采集实验土壤并保留原生微生物群落，种子采自同一区域；设置每盆1、3、5、10株及无植物空白对照（CK）共5个处理，每处理3个生物学重复；测定生长与光合气体交换参数（A、E、Gs、Ci）；采用水培法收集根分泌物并进行LC-MS/MS非靶向代谢组学分析；采集根际土壤（0–2 mm）进行宏基因组测序（NR与KEGG注释）和土壤代谢组学分析；采用PLS-DA/OPLS-DA、LEfSe、单因素方差分析及Tukey–Kramer事后检验，并以Mantel检验和Spearman相关进行多组学整合分析。此外采用基于火山图与Venn图的交集筛选策略鉴定跨密度梯度的核心代谢物，并通过KEGG通路富集与拓扑分析确定核心代谢通路。

研究结果方面，2.1 植物密度显著抑制幼苗生长并诱导光合的气孔限制：根系构型响应最为显著，从单株组密集分枝的发达根系至十株组稀疏衰减的结构，总根长下降91.8%（由446.25降至36.79 cm）、总根表面积下降91.2%（由40.83降至3.61 cm²）、总根体积下降94.1%（由0.974降至0.057 cm³）、根尖数下降93.2%（由893降至61）。所有地上与地下生物量参数随密度增加显著下降（p < 0.05），十株组相对单株组叶面积、根鲜重、根干重、地上鲜重和地上干重分别下降78.94%、81.05%、82.84%、83.56%和78.25%，主根长下降35.66%；株高则表现为从单株组至五株组轻微非显著增加（+8.24%，p > 0.05）后在最高密度显著下降（十株组相对单株组?14.33%，p < 0.05）。效应量分析（Cohen’s d，单株组对十株组）显示所有性状均存在极大密度效应（|d| > 3.7）。光合方面，净光合速率（A）随密度增加逐渐下降，十株组相对单株组显著降低，并伴随气孔导度（Gs）和胞间CO₂浓度（Ci）同步下降，蒸腾速率（E）在十株组下降超过70%；A、Gs和Ci的同步下降表明高密度下光合抑制主要归因于气孔限制而非非气孔限制。上述结果表明高密度种内竞争对刺萼龙葵幼苗造成严重碳同化亏缺，并级联导致生物量积累和根系发育的大幅下降。

2.2 密度相关的根分泌物谱重编程：PLS-DA显示样品沿密度梯度清晰渐进分离，四个处理占据不同的排序空间（R² > 0.99，p = 0.001）；十株组相对单株组有419个代谢物显著上调、101个下调。层次聚类显示五株组与单株组先聚在一起，其次为三株组，而十株组形成独立分支，表明高种植密度引起分泌物组成的最大分异；单株组以Rhodojaponin III、Caffeine和Ferulic Acid为特征，三株组以Barrinic Acid、19Alpha-19-Hydroxy-3,11-Dioxo-12-Ursen-28-Oic Acid和16-Hydroxygypsogenic acid为特征，五株组以4-Nitrophenol、Vulgarin和Glycerol为特征，十株组以Pro-Lys-Ser、P-Toluenesulfonic Acid和Phenoxyacetic acid为特征。通过基于交集策略的核心代谢物筛选，共955个代谢物在三个两两比较中持续差异丰富，其中仅13个可依据HMDB和Metlin数据库推定注释；核心代谢物具有高质量精度（<10 ppm）。Spearman相关分析揭示两种相反响应模式：2-Benzylsuccinic Acid、2-Phenylglycine、3-Methylcrotonylglycine和5-Hydroxycalamenene与所有测定的生长和生理性状（叶面积、主根长、根与地上生物量、A、E、Gs和Ci）显著正相关（p < 0.05）；其余9个代谢物（1,2-Cyclohexanediol, 1-Methyl-4-(1-Methylethenyl)-；3-Methylgallic Acid 5-Sulfate；Acetyl-DL-Valine；Alpha-D-Glucose；Alpinumisoflavone；Arg-Met-Ser；Bakkenolide B；Chrysoobtusin；Citric Acid）与这些性状显著负相关。

2.3 种内竞争重构根际微生物群落及分泌物–微生物关联：基于非冗余基因集的物种注释显示各处理共享4241个核心类群。α多样性分析显示，引入刺萼龙葵后细菌多样性相对未种植对照持续富集，Chao和Shannon指数在所有种植处理中显著升高，Simpson指数相应降低； richness（Chao）与evenness（Shannon）均在单株组达到峰值并随密度增加逐渐下降，但即使在十株组仍高于CK；Simpson指数随密度增加总体下降，在十株组达到最低，五株组和十株组相对CK显著降低（p < 0.05），而单株组和三株组与CK无显著差异。基于Bray–Curtis的PCA显示CK与多数种植处理沿PC1轴明显分离，表明刺萼龙葵对土壤细菌群落组成具有强根际效应；单株处理沿PC2轴与其他处理明显分离，而种植处理呈现密度依赖性聚类，高密度处理（三、五、十株组）沿PC1轴逐渐偏移。门水平上，Actinomycetota、Pseudomonadota和Acidobacteriota在所有根际群落中占主导，Pseudomonadota在十株组明显高于单株组，提示向富营养型和潜在胁迫适应类群的密度相关转变；属水平上Rubrobacter在各处理中保持高丰度，而Nitrososphaera在密度梯度间无显著差异。LEfSe分析（LDA > 3，p < 0.05）鉴定出处理特异性生物标志类群：CK组为Nitrospira，单株组为Flavitalea，三株组为Arthrobacter，五株组为Thermocoleostomius，十株组为Allosphingosinicella与Microbacterium。对优势属的单因素方差分析与Tukey–Kramer事后检验显示：Blastococcus和Arthrobacter在中等密度（三株组）达到峰值后下降；Allosphingosinicella、Ilumatobacter和Chloracidobacterium单调增加并在十株组达到峰值；Nitrososphaera随密度逐渐下降并在十株组达到最低；Nitrospira在植物引入后（单株组）急剧下降后在较高密度部分恢复；Microvirga和Geodermatophilus仅呈中度波动；共11个属显示显著的密度依赖性分化（LDA > 3，p < 0.05）。Spearman相关分析显示这些密度响应属组装为两个与根分泌物相关联的功能对比模块：Allosphingosinicella、Ilumatobacter、Microbacterium和Nitrospira与胁迫相关代谢物（如1,2-Cyclohexanediol, 1-Methyl-4-(1-Methylethenyl)-、Citric Acid、Chrysoobtusin）正相关，与推定生长相关的酚酸类（如2-Benzylsuccinic Acid、5-Hydroxycalamenene）负相关，而Flavitalea呈相反模式；Arthrobacter和Pseudarthrobacter仅与2-Benzylsuccinic Acid相关，Thermocoleostomius与3-Methylcrotonylglycine相关，Blastococcus与2-Benzylsuccinic Acid和2-Phenylglycine正相关而与Arg-Met-Ser和Alpinumisoflavone负相关，unclassified_Candidatus_Rokuibacteriota与1,2-Cyclohexanediol正相关而与3-Methylcrotonylglycine负相关，Microvirga与Chrysoobtusin正相关而与2-Phenylglycine和5-Hydroxycalamenene负相关。这些分泌物定义的微生物模块进一步转化为与植物表现的差异化关联：潜在胁迫适应模块（Allosphingosinicella、Ilumatobacter、Microbacterium和Nitrospira）与生物量和气体交换参数广泛负相关，其中Ilumatobacter除株高外与所有测定性状显著负相关；推定生长相关模块（Flavitalea和Blastococcus）与多数生长和生理性状正相关；unclassified_Candidatus_Rokuibacteriota是唯一与株高显著负相关的属，Thermocoleostomius仅与E正相关，Arthrobacter和Pseudarthrobacter未检测到显著相关。宏基因组功能分析显示，KEGG通路Level 2总体功能类别在各密度间相对稳定，提示高功能冗余；但在Pathway Level 3，与环境感知和运输相关的通路（如Two-component system和ABC transporters）在高密度下富集；CAZy分析显示核心碳水化合物活性酶如Glycoside Hydrolases（GHs）保持稳定，这些功能变化支持竞争加剧下推定胁迫适应类群的响应。

2.4 高密度竞争重塑根际代谢组并富集潜在与生长下降相关的代谢物：土壤代谢物谱的PLS-DA显示清晰的密度依赖性分离（R = 0.84，p = 0.001）；与微生物群落结构的均质化模式一致，三、五、十株组在排序空间中紧密聚类，而CK和单株组占据不同位置，表明种内竞争加剧与土壤化学的趋同变化相关联；单株组与十株组比较鉴定出36个显著上调和142个下调代谢物。层次聚类热图显示：未种植CK组以芳香烃和药物类似物（如Chlorobenzene和Lamivudine Sulfoxide）为主；单株组代谢组富集肽类和长链脂肪酸衍生物（如Phe-Ser-Arg和(S)-10,16-Dihydroxyhexadecanoic Acid）；随密度增加，三株组富集聚酮抗生素和长链酰基肉碱（如Filipin II和O-(17-Carboxyheptadecanoyl)carnitine），五株组以渗透调节碳水化合物（如Trehalose和Raffinose）为主，十株组以鞘脂和甾体糖苷（如SM(14:3_2O/18:4)和Furostane Base-2H+O-Hex）为特征；这种从肽和脂肪酸主导谱向渗透调节物和膜脂主导谱的渐进转变反映了种内竞争加剧下根际累积的生物学胁迫化学特征。KEGG通路富集分析（p < 0.05）鉴定出8条显著富集通路，包括ABC transporters、Galactose metabolism、Arachidonic acid metabolism、Cutin, suberine and wax biosynthesis、Biotin metabolism、Nucleotide metabolism、Tropane, piperidine and pyridine alkaloid biosynthesis和Tryptophan metabolism。KEGG拓扑分析（p < 0.05且通路impact > 0.1）进一步缩小至3条核心通路：Biotin metabolism、Tropane, piperidine, and Pyridine alkaloid biosynthesis和Pantothenate and CoA biosynthesis；某些通路（如Pantothenate and CoA biosynthesis）因高Impact值在拓扑分析中被确定为核心通路，尽管未在常规富集分析中达到统计学显著。这些通路中4个代谢物成为核心贡献者（VIP > 1，FDR < 0.1）：N-Feruloyloctopamine、N-Carboxyethyl-Gamma-Aminobutyric Acid、Glutarate Semialdehyde和Sutherlandin D。Spearman相关分析显示4个核心代谢物均与幼苗生长和光合表现广泛显著负相关：N-Feruloyloctopamine、N-Carboxyethyl-Gamma-Aminobutyric Acid和Glutarate Semialdehyde与所有测定性状（叶面积、主根长、根与地上生物量、A、E、Gs和Ci）显著负相关，Sutherlandin D呈类似但略弱的模式，与株高关联不显著。上述结果表明高密度竞争不仅耗竭根际促生物代谢物，还与具有潜在自毒或化感特性的生物活性化合物积累相关联，可能在土壤中形成化学遗产并强化拥挤条件下刺萼龙葵的生长抑制。

2.5 多组学整合揭示根分泌物–微生物–代谢物级联与种内竞争下根际重塑相关联：研究人员整合植物表型、根分泌物、根际细菌群落和土壤代谢物四层数据集，采用Mantel检验评估四层数据在整体群落水平的关联，采用变量水平Spearman相关分析解析具体两两关系，并构建概念模型总结级联途径。Mantel检验显示植物表型与根分泌物强相关（r = 0.80，p < 0.05），根分泌物与土壤代谢物高度相关（r = 0.83，p < 0.05）；植物表型与根际微生物和土壤代谢物的关联为中等（r = 0.47和0.43，p < 0.05），根分泌物与根际微生物的关联显著但相对较弱（r = 0.27，p < 0.05）；微生物群落与土壤代谢物在整体群落水平无显著相关（r = 0.06，p > 0.05），这种弱关联可能源于整体群落距离矩阵被众多稳定、非响应类群稀释，因此进一步通过变量水平相关分析解析特定类群与代谢物的两两关系。变量水平Spearman相关分析显示密度响应细菌属可分为两种相反相关模式的组：以Flavitalea和Blastococcus为代表的推定生长相关类群与植物生物量正相关，以Allosphingosinicella、Ilumatobacter和Microbacterium为代表的潜在胁迫适应类群与Glutarate Semialdehyde和Sutherlandin D正相关。整合上述多组学证据，研究人员总结了与种内竞争相关的基于相关的根际重塑级联途径：种内竞争加剧与植物表型重塑和根分泌物重编程相关联；分泌物组成变化与根际细菌群落从推定生长相关类群向推定胁迫适应类群的转变相关联；该微生物群落重组进一步与两个生长下降相关土壤代谢物（Glutarate Semialdehyde和Sutherlandin D）的积累相关联，可能形成土壤化学遗产并对后续幼苗生长产生潜在反馈。

讨论部分总结：3.1 种内竞争对幼苗表型和根分泌物的影响表现为显著不对称性，竞争加剧时根系构型显著简化、生物量和光合速率大幅下降，而株高仅轻微增加后下降，表明资源有限条件下植物优先维持地上光竞争能力而减少对根系的碳投入；表型变化改变碳分配模式并最直接表现为根分泌物组成的显著变化，分泌物从酚酸主导转向糖、肽和有机酸主导，Mantel检验显示表型与分泌物强相关（r = 0.80，p < 0.05），表明二者同步发生；根分泌物并非光合碳分配的被动产物而是植物响应环境变化的主动化学策略，酚酸减少与糖和有机酸增加意味着从防御性分泌向资源获取和胁迫适应型分泌的转变。3.2 尽管Mantel检验在整体群落水平显示根分泌物与根际微生物仅弱相关（r = 0.27，p < 0.05），变量水平相关分析揭示清晰的成对关系：弱竞争下占优的酚酸类分泌物（如2-Benzylsuccinic Acid）与Flavitalea和Blastococcus等属正相关，强竞争下富集的糖、肽和有机酸类分泌物（如Citric Acid、Arg-Met-Ser）与Allosphingosinicella、Ilumatobacter和Microbacterium等属正相关；整体水平弱相关可能源于距离矩阵被众多稳定非响应类群稀释，而聚焦特定类群时关联变得清晰，这与根分泌物化学组成可作为影响根际微生物底物利用和定殖能力的选择压力的近期研究一致；研究人员推测竞争加剧时分泌物从酚酸向糖、肽和有机酸转变可能改变根际化学环境从而选择性影响不同微生物类群的生长和定殖，推定生长相关细菌可能更依赖酚酸底物而潜在胁迫适应细菌更适应糖和有机酸主导的根际环境；环境感知和运输相关通路（如Two-component system和ABC transporters）的富集为该转变提供了直接功能支持。3.3 微生物群落与土壤代谢物在整体群落水平无显著相关（r = 0.06，p > 0.05），但变量水平分析显示强竞争下富集的Allosphingosinicella、Ilumatobacter和Microbacterium与Glutarate Semialdehyde和Sutherlandin D显著正相关，而推定生长相关的Flavitalea和Blastococcus与这两个代谢物显著负相关，同时这两个代谢物与植物生物量、根系性状和光合参数显著负相关，表明微生物群落重组与特定土壤代谢物积累同步发生；研究人员提出可检验假设：竞争加剧下富集的潜在胁迫适应细菌可能通过代谢转化或与植物互作促进Glutarate Semialdehyde和Sutherlandin D在土壤中积累，这些代谢物可能进一步反馈抑制后续幼苗生长，但该推断基于相关分析无法建立因果关系，仍需通过外源代谢物添加和关键菌株接种实验验证；本研究将刺萼龙葵种内化感互作框架扩展至活体根系，揭示了种内竞争下从分泌物到微生物重组再到代谢物积累的关联链。3.4 对入侵植物种内竞争策略的意义方面，传统观点认为高密度导致的资源耗竭和自毒效应限制种群扩张，但本研究表明刺萼龙葵并非被动忍受密度胁迫，而是通过多层次根际重塑缓冲该胁迫；需要指出自然高密度种群中密度效应与单株资源可利用性降低本身耦合，盆栽固定土壤体积使高密度必然导致单株资源可利用性降低，因此本研究并非旨在将资源限制与竞争信号分离，而是评估密度胁迫对根际微生态系统的综合生态后果，观察到的从酚酸向糖/肽分泌物的转变及后续微生物重组代表植物应对该复合胁迫的主动策略而非单纯的养分耗竭被动结果；本研究观察到竞争加剧时植物优先维持株高，研究人员推测这种株高维持可能是争夺光的策略，但需遮荫实验区分地上与地下竞争；该“分泌物–微生物–代谢物”级联途径可能帮助植物在高密度种群中维持根际功能的相对稳定；该根际重塑策略为入侵生态学中的“密度悖论”提供了潜在解释，即入侵植物如何在高密度下维持高种群增长；经典假说如“天敌释放”和“生物抗性”主要关注新入侵范围内专性天敌的缺失，却很大程度上忽视了植物自身主动调节根际微环境以缓冲种内竞争的能力；本研究提示入侵植物可能通过分泌物相关的微生物重组和土壤化学遗产在高密度下构建对自身有利或至少不损害自身种群建立的根际条件；当然本研究基于相关分析无法建立因果关系，该级联途径的普遍性及其在野外多物种群落中的实际作用仍需通过外源添加、菌株接种和野外密度梯度实验进一步验证。3.5 局限性方面：破坏性拔取和24 h水培收集可能诱导伤口胁迫和微生物降解，未来应采用非破坏性原位收集方法（如吸附树脂陷阱或根系观测窗）或无菌水培系统；固定土壤体积的盆栽设计无法将竞争信号与资源限制分离，观察到的根际效应应解释为对二者的综合响应，未来需采用分根系统或持续养分补充解耦这些机制；非靶向代谢组学注释率有限，鉴定的代谢物应视为推定鉴定并需验证；小样本量（n = 3）增加过拟合风险；尽管在KEGG Pathway Level 3鉴定了特定功能变化，分析仅限于通路水平丰度，未广泛开展特定基因水平注释等更深层次探索；此外本研究主要基于相关分析无法建立直接因果关系，提出的级联途径为概念模型，未来需外源代谢物添加和菌株接种实验验证因果关系。

研究结论部分翻译：本研究整合植物表型、根分泌物、根际细菌群落和土壤代谢物四层数据，揭示了刺萼龙葵在种内竞争下的根际重塑策略。种内竞争对幼苗的影响呈不对称性：根系构型简化，生物量和光合速率显著下降，而株高仅轻微增加后下降，表明植物优先维持地上光竞争能力。随着种内竞争加剧，根分泌物从酚酸主导转向糖、肽和有机酸主导，分泌物组成变化与根际微生物重组同步发生：弱竞争下占优的酚酸类分泌物与Flavitalea和Blastococcus等推定生长相关细菌正相关，而强竞争下富集的糖/肽类分泌物与Allosphingosinicella、Ilumatobacter和Microbacterium等潜在胁迫适应细菌正相关。微生物重组进一步与土壤中Glutarate Semialdehyde和Sutherlandin D的积累相关联；这两个代谢物与植物生物量、根系性状和光合参数显著负相关，可能构成对后续幼苗生长产生潜在反馈的土壤化学遗产。总之，刺萼龙葵在种内竞争下通过“分泌物–微生物–代谢物”级联重塑根际，该策略可能帮助其耐受高密度种群胁迫并部分解释其入侵成功。由于本研究基于相关分析，上述级联途径的因果关系仍需通过外源添加和菌株接种实验进一步验证。