磷供应依赖性调控大豆叶片缺铁响应的多组学整合分析：从转录重编程到代谢缓冲的层级调控机制

《Environmental and Experimental Botany》：Multi-Omics Profiling of Soybean Leaves under Fe Deficiency with Different Phosphorus Supply

　　

磷（phosphorus, P）与铁（iron, Fe）是植物生长发育必需的两种矿质元素，二者在土壤化学行为与植物代谢过程中存在复杂互作，然而磷有效性如何在分子层面调控缺铁响应尚不清晰。研究人员以大豆品种CN15为材料，在缺铁（30 μM Fe-EDTA）条件下设置三个磷供应梯度（0 μM、150 μM、300 μM KH₂PO₄），并以正常铁营养（80 μM Fe-EDTA, 300 μM KH₂PO₄）为对照，整合表型测定、转录组测序及非靶向代谢组学分析，系统解析不同磷水平下大豆幼苗叶片对缺铁胁迫的分子响应机制。结果表明，磷供应从根本上重塑了缺铁响应的转录组格局——在三个低铁处理与对照的比较中仅共享一个差异表达基因（differentially expressed gene, DEG），而低铁处理间的比较则共享262个DEG。缺铁条件下磷供应增加促进了地上部干重（300 P处理较0 P提高44.4%），但高磷加剧了磷铁拮抗，表现为叶片铁浓度在300 P处理下较150 P下降11.7%。代谢组学分析鉴定出1,381个差异代谢物，其中360个在三个低铁处理中共享。整合分析揭示了一条由磷有效性驱动的"感知–整合–执行–缓冲"层级调控链：ABC转运蛋白（尤其是ABCB1，即ATP结合盒B亚家族成员1）可能参与营养状态与碳分配的耦合；乙烯信号通路（以ERS1，即乙烯响应传感器1为代表）整合并放大胁迫信号；戊糖磷酸途径（pentose phosphate pathway, PPP）与PRS1（磷酸核糖焦磷酸合酶1）执行碳骨架的差异化分配；VTE4（γ-生育酚甲基转移酶）介导终端氧化还原缓冲。Mantel检验进一步确认PRS1、G6PD（葡萄糖-6-磷酸脱氢酶）和ABCB1为连接基因表达与代谢及表型变异的核心枢纽。该研究为培育大豆耐磷铁复合胁迫品种提供了多个潜在遗传靶点。

**论文解读**

**一、研究背景与问题提出**

磷和铁是植物生长发育必需的两大矿质营养元素。然而，全球约30%的耕地因土壤磷的快速固定而处于缺磷状态，严重制约作物产量。植物为适应低磷胁迫进化出根系构型重塑、根系分泌物增加及与微生物共生等多种磷活化策略，且磷利用效率存在显著的基因型差异。与此同时，尽管铁在地壳中含量丰富，但在碱性或石灰性土壤中，铁主要以难溶性三价铁氧化物形式存在，生物有效性极低，常引发植物缺铁黄化。植物依赖还原策略和螯合策略进行铁吸收，而在酸性或淹水还原条件下，过量亚铁离子则通过芬顿反应引发氧化胁迫，植物需借助根表铁膜形成、胞内区隔化、铁蛋白结合及抗氧化防御系统强化等机制抵御铁毒害。

值得关注的是，磷与铁在土壤化学行为及植物体内代谢过程中并非彼此独立，而是存在复杂的相互作用。在土壤层面，磷与铁可形成难溶性复合物，导致二者有效性呈负相关；在植物体内，过量施磷会阻碍铁的吸收与向地上部的转运。这种磷铁互作在水稻、小麦等单子叶植物以及菜豆、大豆等双子叶植物中均有报道，暗示其可能是跨植物谱系的保守特征。大豆作为对铁和磷需求均较高的重要双子叶作物，是探究双子叶植物磷铁互作分子机制的适宜体系。然而，磷有效性如何在分子层面系统调控缺铁响应，其背后的转录重编程和代谢网络重塑机制尚不清晰。基于此，研究人员提出科学假设：外源磷供应将改变缺铁条件下大豆叶片的转录组和代谢组，且响应模式因磷水平而异，并据此开展整合表型、转录组和代谢组的多组学研究。

**二、研究目的与意义**

本研究旨在通过整合表型测定、转录组学与非靶向代谢组学数据，系统表征缺铁响应中磷依赖性的分子变化特征，并鉴定可能参与磷铁交互对话的调控基因与代谢通路。研究不仅有助于从系统生物学层面深入理解植物矿质营养稳态的协同调控机制，也为培育耐受磷铁复合胁迫的大豆品种提供可操作的遗传靶点和理论依据。

**三、关键技术方法**

研究人员选用大豆品种CN15，在V1期（第一片三出复叶期）开始实施四个营养处理：对照（CK）为铁充足条件（300 μM KH₂PO₄, 80 μM Fe-EDTA），其余三个处理均在缺铁条件下（30 μM Fe-EDTA）分别设置0 μM、150 μM和300 μM KH₂PO₄的磷梯度。在V3期采集叶片样品进行表型测定（地上部干重，n = 5；电感耦合等离子体发射光谱测定叶片磷铁含量，n = 3）、转录组测序（Illumina HiSeq2500平台，每处理3个生物学重复，采用HISAT2比对、StringTie组装定量、DESeq2差异分析，以FDR < 0.05且|fold change| ≥ 2为差异基因筛选标准）以及非靶向代谢组分析（Waters ACQUITY UPLC系统联用Thermo Fisher Q Exactive质谱，正负离子模式采集，每处理5个生物学重复，采用OPLS-DA模型以VIP > 1且|log₂FC| ≥ 1或ANOVA P < 0.05筛选差异代谢物）。基因与代谢物的功能注释和通路富集分别基于GO、KEGG数据库（超几何检验，FDR ≤ 0.05），多组学整合采用通路共富集分析、相关性聚类及Mantel检验，并选取五个核心基因（ABCB1、ERS1、PRS1、G6PD、VTE4）进行qRT-PCR验证。

**四、主要研究结果**

**（一）表型响应：磷供应对生长和营养分配的差异性影响**

与对照相比，所有低铁处理的地上部干重均显著降低，其中0 P处理降幅最大（44.6%），表明磷铁双重缺乏严重抑制幼苗生长。缺铁条件下，地上部干重随磷供应增加而递增（300 P较0 P提高44.4%），但仍显著低于对照（降幅20.0%）。叶片磷铁浓度分析揭示了复杂的非线性磷铁互作模式：缺铁条件下，150 P处理的叶片铁浓度与对照无显著差异，而300 P处理的铁浓度较150 P下降11.7%；叶片磷浓度在300 P处理下未进一步增加，反而较对照显著降低11.3%。上述结果表明磷铁之间同时存在双向拮抗与协同的精细平衡。

**（二）磷供应重塑缺铁转录组景观**

转录组测序数据质量可靠（Q30 ≥ 93.76%，唯一比对率 ≥ 93.21%）。主成分分析显示四个处理组清晰分离，PC1与PC2共解释71.9%的方差。韦恩分析揭示了一个关键发现：在0 P vs CK、150 P vs CK及300 P vs CK三个比较中，仅共享一个差异表达基因，表明绝大多数缺铁响应基因的转录调控并非仅由铁状态决定，而是高度依赖磷供应水平。相比之下，低铁处理间的三个比较共享262个DEG，构成对磷浓度高度敏感的核心调控模块，广泛涉及磷吸收转运、信号转导、碳氮代谢及次生代谢。GO分类显示，0 P以广泛上调为主，150 P呈上下调平衡模式，300 P以下调为主导。KEGG富集分析表明，0 P突出光合作用相关通路（如光合作用-天线蛋白，21个DEG），150 P突出次生代谢（如倍半萜与三萜生物合成，6个DEG），300 P富集通路最少（如亚油酸代谢，4个DEG）。磷供应由此重定向了缺铁条件下的转录优先序：0P表现为光合补偿，150P恢复光合稳态并诱导次生代谢防御，300P则抑制整体响应但激活防御与激素信号通路，与加剧的磷铁拮抗一致。

**（三）代谢组重构：膜脂重塑与氧化还原稳态**

代谢组数据质量可靠，QC样品紧密聚类，OPLS-DA模型解释力良好（Q²Y = 0.859）。韦恩分析共鉴定出1,381个差异代谢物，其中360个在三组低铁处理中共享，表明铁有效性降低引发了广泛的代谢重编程。低铁处理间比较则鉴定出999个差异代谢物，仅57个为三组共享，表明不同磷供应在缺铁条件下引发部分特异、部分重叠的代谢响应。KEGG富集显示，0 P显著富集ABC转运蛋白（10个代谢物），150 P富集嘌呤代谢和甘油磷脂代谢，300 P响应最广，涵盖异黄酮生物合成、嘌呤代谢及泛醌/萜醌生物合成。值得注意的是，甘油磷脂代谢和α-亚麻酸代谢在150 P vs CK中特异富集，提示磷依赖性的膜脂重塑机制：光系统I和细胞色素b₆f复合体富含铁硫簇，类囊体膜对缺铁高度敏感，而磷铁互作可能通过影响膜脂脂肪酸去饱和程度来维持膜的流动性和功能。此外，戊糖磷酸途径在0 P vs 300 P中显著富集，伴随其关键限速酶G6PD转录上调，提示高磷加剧缺铁诱导的氧化胁迫，迫使碳流转向NADPH再生以维持氧化还原稳态。

**（四）多组学整合鉴定核心调控枢纽**

整合分析聚焦四个代表性通路：ABC转运蛋白（物质运输与离子稳态核心）、戊糖磷酸途径（抗氧化与碳代谢重排关键节点）、植物激素信号转导（整合胁迫信号的中心调控网络）及泛醌/萜醌生物合成（线粒体能量与氧化还原平衡代表）。通路共富集显示，ABC转运蛋白在缺铁条件下普遍激活，0 P vs 300 P比较中基因富集数最高（18个）；戊糖磷酸途径在高磷加剧缺铁时最为显著（0 P vs 300 P中37个基因）；植物激素信号转导在0 P vs 300 P中富集基因数高达204个，远超其他比较。相关性聚类分析鉴定出五对关键转录本–代谢物配对关系：ABCB1（Glyma.10G055000）与蔗糖、G6PD（Glyma.08G199000）与α-D-葡萄糖-6-磷酸、ERS1（Glyma.03G216700）与反式玉米素、VTE4（Glyma.12G014300）与植基二磷酸、以及PRS1（Glyma.05G070500）与α-D-葡萄糖-6-磷酸。qRT-PCR验证了这五个基因的表达趋势与RNA-Seq数据一致。

**（五）Mantel检验确认核心枢纽**

Mantel检验整合基因表达距离、代谢物丰度距离与表型变异距离，发现PRS1、G6PD和ABCB1的表达距离与代谢组及综合表型（地上部干重、铁浓度、磷浓度）变异均显著相关，其中PRS1相关性最强。基因间表达相关性分析揭示两个正相关模块——ERS1-VTE4（r = 0.884）与G6PD-ABCB1（r = 0.806），二者呈拮抗连接，而PRS1与G6PD（r = -0.764）和ABCB1（r = -0.715）均呈负相关，充当两模块间的桥梁。由此确认PRS1、G6PD和ABCB1为连接基因表达与代谢及表型变化的核心枢纽。

**五、讨论与结论**

在讨论部分，研究人员基于多组学整合结果提出一个"感知–整合–执行–缓冲"的层级调控工作模型：ABC转运蛋白（以ABCB1为代表）可能参与营养状态感知与碳分配耦合——ABCB1与蔗糖的强正相关暗示其可能在碳代谢状态与激素信号间发挥协调功能；乙烯信号通路（以ERS1为代表）承担整合与放大胁迫信号的核心角色——ERS1与反式玉米素的显著正相关及0 P vs 300 P中多达204个激素信号转导相关DEG的支持了该推论；戊糖磷酸途径和PRS1分别执行碳骨架向NADPH再生（氧化防御）和核苷酸合成的差异化分配——PPP在300 P下的富集与G6PD上调提示高磷加剧氧化胁迫下的碳流重定向，而PRS1在低磷下的持续上调及与葡萄糖-6-磷酸的负相关则指示另一条碳利用分支；VTE4作为终端氧化还原缓冲节点——其在0 P下上调、300 P下下调的双向表达模式与膜系统保护优先级的转变相呼应，且VTE4既是PPP产生NADPH的下游功能出口，也是连接代谢重编程与膜系统完整性保护的最终环节。

需指出的是，该模型完全基于相关性数据构建，文中以工作假设形式提出，尚待直接的遗传学和生物化学实验验证。未来研究应增加转录组生物学重复以提高对细微变化的检出能力，将分析扩展至根系组织以解析地上部系统性响应之前的早期信号事件，并补充叶绿素含量、光合参数及抗氧化酶活性等生理指标，以加强分子特征与生理状态之间的关联。此外，鉴于已鉴定的核心基因–代谢物–表型关联均为相关性关系，后续功能验证实验对于确认这些候选因子的调控角色和完善所提出的工作模型至关重要。

**研究结论**：本研究表明磷有效性是重塑大豆缺铁分子响应的主控开关。通过整合生理表型、转录组和代谢组数据，研究人员解析了一个层级组织且相互关联的调控网络。基于整合数据，鉴定了四个核心基因——ERS1、ABCB1、PRS1和VTE4，其表达模式与特定代谢物和表型性状紧密关联。研究人员提出一个工作模型，其中这些基因及其关联模块可能构成一条"感知–整合–执行–缓冲"链。该模型为未来机制研究提供了概念框架，但仍待实验验证。