地黄根际食物网营养重组在第一季种植期间即呈现品种间差异

《Agronomy》：Early Trophic Reorganization of the Rhizosphere Food Web Differs Between Two Cultivars of Rehmannia glutinosa

　　

摘要专业翻译 连作可导致土壤退化与重茬问题，但土壤食物网的分化是否在首个种植季内即已发生，目前仍不清楚。本研究探讨了地黄（Rehmannia glutinosa）首季种植期间线虫群落组成与食物网功能的时序变化。研究人员评估了不同生长阶段、品种和土壤区室下属级线虫组成、营养结构、生态学指数和代谢足迹（metabolic footprint）。进一步整合了土壤理化性质、代谢物谱以及细菌和真菌群落，以考察它们与线虫群落变异和功能分配的关系。线虫群落在生长季内发生显著变化，采样时间是群落变异的主导来源。食细菌线虫（bacterivore）构成主要营养背景，而植物寄生线虫和食真菌线虫（fungivore）则表现出品种依赖性和阶段依赖性的增加。与敦王（DW）品种相关的土壤，尤其是根际土壤，持续表现出向植食性代谢通道的偏移。在若干处理中，植物寄生线虫尽管在总丰度中占比很小，却对总代谢足迹贡献不成比例，表明丰度与足迹之间存在明显的解耦现象。整合分析显示，线虫群落组成与土壤理化条件、细菌、真菌及代谢物谱广泛相关，而代谢足迹分配则与微生物群落组成尤其是真菌的关系更为密切。这些发现表明，线虫群落和食物网的分化可在首个种植季内出现，早于连作障碍完全形成。代谢足迹提供了丰度指标无法捕捉的功能信息，可作为地黄栽培中早期地下生态风险的敏感指示指标。

论文解读

一、研究背景与科学问题

连作和缩短轮作在集约化农业系统中被广泛采用，但同一土壤上重复种植同一种作物常导致产量下降、土壤退化和重茬问题。这些问题并非由单一因素引起，而是源于土壤理化性质、养分有效性、自毒化合物、土传病原菌和地下微生物群落的交互变化。植物–土壤反馈进一步表明，作物引起的土壤生物变化可影响后续植物的表现，在农业生态系统中形成正反馈或负反馈循环。因此，连作障碍不仅是一个作物生产问题，更是一个涉及多个地下组分相互作用的土壤生态学问题。

土壤线虫是土壤食物网的重要成员，其营养分类——包括食细菌线虫（bacterivore，Ba）、食真菌线虫（fungivore，Fu）、植物寄生线虫（plant parasite，Pp）和杂食–捕食线虫（omnivore–predator，Om）——为将土壤动物与微生物分解途径、根部取食压力和高营养级调控联系起来提供了实用框架。由于不同线虫类群在取食习性、生活史策略和环境扰动敏感性上存在差异，线虫群落结构已被广泛用作土壤条件和食物网状态的生物指示指标。线虫动物区系分析进一步整合富集指数（enrichment index，EI）和结构指数（structure index，SI），以评价土壤生态系统中的资源可利用性和食物网复杂性。除丰度和群落组成外，线虫代谢足迹可估算不同线虫功能团的生物量生产和呼吸碳利用，提供碳在营养通道间分配的功能信息。

已有研究表明，连作可强烈改变线虫群落和土壤食物网结构。在连作草莓系统中，土壤pH下降、酚酸积累、线虫多样性降低，优势分解途径从细菌通道转向真菌通道。在长期花生连作中，连作年限增加与微生物食性线虫减少、植物寄生线虫增加、线虫多样性降低以及线虫食物网内碳通量失衡相关。连续棉花连作也被报道会影响线虫代谢足迹和营养连接度。这些研究主要描述的是土壤退化或重茬问题已经形成的长期系统，而线虫食物网结构是否在首个种植季内即开始转变，仍不清楚。

地黄（Rehmannia glutinosa）是一种对连作障碍高度敏感的药用作物。已有研究表明，地黄连作会重塑根际真菌群落、促进自毒物质积累并改变根系分泌物–微生物组相互作用。然而，与酚酸、根系分泌物、土传病原菌以及细菌和真菌群落相比，土壤动物组分受到的关注较少。特别是地黄系统中线虫营养结构和代谢足迹分配的表征仍十分有限，限制了对土壤食物网在早期地下生态变化中如何响应的理解。

为填补这一知识空白，研究人员考察了地黄首季种植期间线虫群落组成和食物网功能的时序变化，旨在确定线虫指示指标能否在连作障碍完全形成之前揭示早期的地下生态分化。

二、研究设计与主要技术方法

田间试验于2023年在中国河南省焦作市武陟县乔庙乡（113°35′E，35°03′N）进行，该地为地黄主产区，土壤类型为褐土。两个地黄品种——金九（JJ）和敦王（DW）——于2023年4月同时种植，处于地黄首年种植状态。试验包含两个品种、两个土壤区室（根际土和 bulk 土）和四个采样月份（2023年7月、8月、10月和11月），每个品种×区室×月份组合设置三个独立生物学重复，共48份土壤样品。

主要技术方法包括：采用改良贝尔曼漏斗法提取土壤线虫并在属级水平进行形态学鉴定；基于营养类群和c–p值（colonizer–persister value）进行功能团分类；计算香农多样性指数、成熟度指数（maturity index）、富集指数和结构指数；按照Ferris的方法估算各类代谢足迹（bacterivore footprint，BFP；fungivore footprint，FFP；herbivore footprint，HFP；omnivore–predator footprint，PFP；enrichment footprint，EFP；structure footprint，SFP；total metabolic footprint，MFP）；采用三因素方差分析、基于Bray–Curtis距离的主坐标分析（PCoA）和置换多元方差分析（PERMANOVA）检验各因素效应；对匹配的7月和11月样本整合土壤理化数据、代谢组学数据（UHPLC-Q-Exactive HF-X非靶向代谢组学）以及细菌16S rRNA和真菌ITS扩增子测序数据，开展Mantel检验、偏距离冗余分析（db-RDA）、跨层Spearman相关网络和偏最小二乘路径模型（PLS-PM）等多层整合分析。

三、主要研究结果

3.1 土壤理化性质的时序变化

土壤理化性质在首季种植期间表现出显著的时序变化和处理特异性差异。土壤pH在各处理中保持弱碱性范围，根际土壤pH略低于 bulk 土。总酚酸浓度主要受品种和采样月份影响，JJ相关土壤普遍高于DW相关土壤，并在10月出现较高值。微生物生物量碳（MBC）在8月下降、10月显著升高、11月再次下降。根际土壤的铵态氮（NH₄⁺-N）浓度普遍高于对应bulk土，而硝态氮（NO₃^?-N）则呈现相反模式，在bulk土中相对富集。

3.2 线虫群落组成的时序与品种依赖性变化

共鉴定出44个线虫属，包括13个植物寄生属、13个食细菌属、6个食真菌属和12个杂食–捕食属。食细菌线虫在大多数样本中构成最大营养类群，形成首季种植期间以细菌通道为背景的食物网结构。在此背景之上，植物寄生线虫表现出处理特异性和阶段特异性的增加。主坐标分析表明采样月份是属级线虫群落组成分化的主要来源，PERMANOVA确认月份解释了33.3%的总变异，月份×品种交互作用解释了14.1%，品种单独解释了11.3%，而土壤区室仅解释了2.0%。在属级水平上，Acrobeloides平均相对丰度最高（23.8%），其次为Aphelenchus（9.2%）、Cephalobus（8.7%）、Dolichodorus（8.2%）和Eucephalobus（7.0%）。Dolichodorus在JJ根际土中显著增加（7月至10月从4.0增至33.7个/100g干土），而Heterodera和Meloidogyne则更特征性地出现在DW相关土壤中。

3.3 线虫食物网状况与代谢足迹分配

EI–SI动物区系图谱显示大多数处理均值位于低至中等富集和结构范围内，表明首季种植期间线虫食物网整体结构化程度不高。DW相关土壤在中后期表现出更大的总代谢足迹，且主要由植食性足迹（HFP）驱动。DWR从8月起HFP成为主导组分，在10月和11月占总代谢足迹的大部分。植食性偏移指数（herbivory-shift index，HFP/(BFP+FFP)）表现出显著的品种×区室×月份交互作用（F=16.93，p<0.001），DWR在8月、10月和11月分别达到10.94、10.60和16.64，DWN在8月出现33.59的瞬时高峰，而JJ处理全年维持在较低水平。丰度–足迹解耦分析显示，在DWR 10月处理中，植物寄生线虫仅占总丰度的24.7%，却贡献了总代谢足迹的89.5%；DWR 11月为32.0%对92.0%；DWN 8月为42.6%对95.8%。

3.4 线虫群落分化与植物寄生风险的整合驱动因素

Mantel检验显示，属级线虫群落组成与土壤理化性质（r=0.585，p=0.0005）、真菌群落组成（r=0.522，p=0.0005）、细菌群落组成（r=0.392，p=0.0005）和全代谢物谱（r=0.295，p=0.0030）显著相关。相比之下，代谢足迹谱与真菌群落组成的关联最强（r=0.379，p=0.0005），其次为细菌群落组成（r=0.235，p=0.0115），与土壤理化性质（r=0.049，p=0.3065）和全代谢物谱（r=0.216，p=0.0520）的关联较弱。偏db-RDA在控制月份效应后，保留的土壤变量、代谢物PC1、细菌PCoA1和真菌PCoA1共同解释了线虫群落组成17%的变异（调整R²）。跨层关联网络显示Dactylonectria和Heterophoma与植物寄生或HFP相关节点呈正相关，而Chaetomium和Lysobacter与HFP相关指标呈负相关。月份校正后的PLS-PM显示，土壤性质到代谢物谱、细菌群落和植物寄生压力代理的路径系数分别为0.41、?0.43和0.39；代谢物谱到细菌群落、真菌群落和植物寄生压力的系数分别为0.68、?0.56和?0.32；细菌和真菌群落连接到植物寄生压力代理的系数分别为0.25和?0.64。

四、讨论与结论总结

讨论部分指出，线虫群落分化可在首季种植期间出现，采样时间解释了属级群落变异的最大比例，且品种间呈现不同轨迹。与长期草莓和花生连作系统中报道的线虫多样性广泛丧失和植物寄生线虫强烈优势不同，本研究中观察到的变化主要表现为异步群落更替和营养再分配，因此更适合视为早期生态分化而非已完全形成的连作障碍。植物基因型可通过根系结构、组织质量、化学防御、寄主适宜性和根系分泌物组成差异来改变地下群落。土壤区室效应相对有限，可能因为线虫是运动的消费者而非根际的空间固定组分。

丰度–足迹解耦是本研究的关键发现。丰度指标描述类群的数值代表性，而代谢足迹根据体重、生产量、呼吸和生活史特征对类群加权。DW相关土壤中不成比例的植食性足迹表明，线虫群落估计代谢活性中更大份额集中在植食通道，可能反映了持续的寄主关联资源利用。代谢活性在单一通道中的集中可产生更不对称的能量结构，可能代表营养平衡的早期降低。

整合分析揭示了线虫群落组成和代谢足迹分配与土壤–代谢物–微生物组环境的不同关联模式。群落组成与土壤理化条件以及细菌、真菌和代谢物谱广泛对应，而足迹分配与微生物群落组成的关系更密切。真菌群落组成与线虫群落和足迹分配均密切相关，提示真菌可能构成根际变化与线虫食物网功能之间的重要界面。代谢物对线虫的影响不能简化为对单一化合物的直接响应，植物代谢物可能通过代谢物–微生物组–线虫的顺序途径影响植物寄生线虫。

研究结论指出，地黄首季种植期间可出现清晰的线虫群落和食物网功能的时序与品种依赖性分化。季节进程是群落变异的主导来源，而土壤区室效应相对较弱且依赖于品种和生长阶段。DW相关土壤——尤其是根际——发展出向植食通道的代谢足迹不成比例分配，这种功能转变远强于植物寄生线虫相对丰度所暗示的程度。线虫群落组成和代谢足迹分配反映了地下环境的不同组分，二者为检测早期地下生态分化提供了互补信息。然而，这些模式应被视为潜在前兆而非已建立连作障碍的直接证据，需要多年试验将线虫动态与根系侵染、植物表现和产量联系起来，以确定本研究中观察到的早期营养失衡是否发展为具有农艺后果的重茬问题。