《Agronomy》：Integrated Transcriptomic and Metabolomic Analyses Reveal the Auxin-Mediated Regulatory Network Governing Alfalfa Responses to Phosphorus Deficiency Stress
编辑推荐：
生长素在植物响应低磷胁迫中发挥正向作用。然而，吲哚-3-乙酸（IAA）介导紫花苜蓿（Medicago sativa L.）缺磷响应的分子机制仍不清楚。研究人员以紫花苜蓿品种‘Zhongmu No. 3’为材料，设正常磷（NP，1000 μM KH2PO4）和低
生长素在植物响应低磷胁迫中发挥正向作用。然而，吲哚-3-乙酸（IAA）介导紫花苜蓿（Medicago sativa L.）缺磷响应的分子机制仍不清楚。研究人员以紫花苜蓿品种‘Zhongmu No. 3’为材料，设正常磷（NP，1000 μM KH2PO4）和低磷补充1 μM IAA（LP + IAA，10 μM KH2PO4）两个处理，在处理后48 h和10 d测定形态性状、生理参数、转录组谱式和代谢物积累。与NP组相比，LP + IAA组生长相关性状发生显著变化：48 h时株高增加34%，总根长为NP组1.24倍，侧根数增加23.48%；10 d时根鲜重增加28%，总根长为NP组1.25倍，侧根数增加36.58%，根体积增加22.2%，而株高、茎粗和地上部鲜重与NP组相近；48 h时根酸性磷酸酶活性比NP组高114.8%。转录组分析在48 h和10 d分别鉴定到1274和2277个差异表达基因（DEGs），48 h上调基因包括氨基磷脂ATP酶9（amino phospholipid ATPase 9）和氨基醇磷酸转移酶1（amino alcohol phosphotransferase 1），10 d上调基因包括磷酸转运蛋白1（phosphate transporter 1）和紫色酸性磷酸酶12（purple acid phosphatase 12）。代谢组分析在48 h和10 d分别鉴定到308和1296个差异积累代谢物（DAMs），早期响应富集于嘌呤代谢并涉及（5′-phosphoribosyl）-5-formamido-4-imidazolecarboxamide（FAICAR），长期处理涉及L-天冬氨酸、腺嘌呤和cAMP。整合分析确定色氨酸代谢、半胱氨酸与甲硫氨酸代谢、甘油磷脂代谢和植物激素信号转导为主要调控途径；CDP-胆碱积累及脂质重塑基因变化进一步表明膜磷脂重塑增强。上述结果揭示了外源IAA存在下紫花苜蓿对低磷的形态、生理、转录和代谢响应。
研究背景方面，紫花苜蓿（Medicago sativa L.）为优质豆科牧草，磷（P）是植物生长发育大量元素，约全球半数耕地存在有效磷不足，土壤有机磷占比25%–56%且难被直接吸收，磷肥当季回收率仅10%–25%，过量施磷又引发水体富营养化，缺磷会抑制苜蓿生长与产量。植物生长受生长素（auxin，IAA）、乙烯、赤霉素等激素与缺磷信号互作调控，其中生长素通过调控根细胞分裂、伸长、侧根形成参与根构型重塑；模式植物中缺磷可改变局部IAA分布促进侧根发育，苜蓿中ARF等生长素信号基因已有研究，但IAA如何协同转录与代谢网络调控苜蓿缺磷适应仍不清楚，因此研究人员开展转录组与代谢组整合研究。
该研究发表于《Agronomy》，研究人员以紫花苜蓿品种‘Zhongmu No. 3’为对象，设正常磷（NP）与低磷加1 μM IAA（LP + IAA）处理，在48 h和10 d取样，证明外源IAA在低磷下显著促进根系构型改变、激活磷动员相关酶活、诱导转录与代谢重编程，并揭示色氨酸代谢、半胱氨酸与甲硫氨酸代谢、甘油磷脂代谢及植物激素信号转导共同构成IAA介导的缺磷适应调控网络，对减少磷肥依赖、改良牧草抗逆育种具有理论意义。
关键技术方法上，样本来自紫花苜蓿‘Zhongmu No. 3’水培苗，分NP与LP + IAA两组、48 h与10 d两个时间点；采用RNA-seq进行链特异性转录组测序与DESeq2差异基因分析，UPLC-QTOF-MS进行非靶向代谢组检测与OPLS-DA差异代谢物筛选，qRT-PCR验证转录本，O2PLS与KEGG/Pathview做转录-代谢整合分析，SPSS与R完成统计与通路制图。
研究结果如下。
3.1 外源IAA对低磷下苜蓿形态性状的影响：通过LA-S根析系统测根构型，研究人员发现IAA组48 h株高增34%、总根长与侧根数显著增加；10 d根鲜重、总根长、侧根数、根表面积和根体积均高于NP组，而茎粗与地上部鲜重无大变化，说明IAA主要重塑地下部吸收器官。
3.2 外源IAA对低磷下苜蓿生理性状的影响：通过酸磷酸酶/植酸酶试剂盒与元素测定，研究人员发现48 h叶、根酸磷酸酶和植酸酶活性升高，根酸磷酸酶升114.8%；10 d叶绿素与叶氮高于NP组，总磷低于NP组，总硫高于NP组，提示磷动员、硫代谢与膜脂替代过程被激活。
3.3 转录组分析：经Illumina NovaSeq测序比对至Medicago sativa V3基因组，48 h得1274个DEGs、10 d得2277个DEGs；48 h上调ALA9、AAPT1，下调PLC6、LPCAT1、PEAMT；10 d上调PHT1、PAP12、GH3.1、PLDζ1、GDPD、PLA1。GO/KEGG显示短期富集MAPK与β-丙氨酸代谢，长期富集植物激素信号、糖酵解、甘油磷脂代谢等。TF分析见AP2/ERF-ERF、MYB、WRKY、bHLH、NAC、bZIP、GRAS大量响应。qRT-PCR验证16个基因中81.25%趋势一致。
3.4 代谢组分析：UPLC-QTOF-MS显示48 h有308个DAMs、10 d有1296个DAMs，共有105个；48 h富集嘌呤代谢与FAICAR，10 d富集α-亚麻酸代谢、戊糖与葡萄糖醛酸转化、抗坏血酸与醛糖酸代谢、淀粉蔗糖代谢。
3.5 转录组与代谢组整合分析：O2PLS与相关网络显示，48 h甘氨酸-丝氨酸-苏氨酸代谢与嘌呤代谢仅单点关联；10 d L-天冬氨酸连8个基因，次黄苷、腺嘌呤、cAMP、鸟苷形成嘌呤辐射网。整合KEGG中48 h代谢端嘌呤代谢显著、转录端β-丙氨酸代谢显著；10 d转录端富集甘油酯、核黄素、碳代谢，代谢端富集α-亚麻酸、玉米素生物合成等。关键途径中，色氨酸代谢里YUC与TAA1/TAR上调、ACO差异、ASMT下调；半胱氨酸与甲硫氨酸代谢里PHGDH上调、OAS-TL下调、SAMDC与ACO诱导；甘油磷脂代谢里AGPAT、LPCAT、LPIAT、PLC、GDPD等上调且CDP-胆碱积累；激素信号里AUX1受抑、TIR1上调、ARF下调、GH3强诱导，CK受体CRE1下调、GA受体GID1下调、DELLA上调，说明IAA统筹多种激素信号。
讨论部分总结：研究人员指出根比地上部更敏感，因根直接感知无机磷酸（Pi）且IAA直接调控侧根起始与伸长；酸磷酸酶与植酸酶增强支持有机磷动员与胞内Pi回收，但未直接测磷再利用率；总磷低源于外供磷少，总硫高与SQDG等硫脂替代膜磷脂有关；PLDζ1与GDPD长期上调印证膜脂降解与磷回收被激活；ALA9属氨基磷脂转运家族，参与早期膜脂运输；PEAMT短期下调与拟南芥结果不同，可能缘于物种、时长、组织及外源IAA；AP2/ERF、MYB、WRKY、bHLH等TF家族整合缺磷与其他非生物胁迫；TIR1/AFB–AUX/IAA–ARF模块中长期TIR1上调、AUX1下调、GH3诱导，平衡根生长与缺磷适应；YUC与TAA1/TAR上调表明色氨酸依赖的IAA生物合成被激活，与水稻OsTAR2/OsYUC8机制类似。
结论部分翻译：与NP组相比，IAA组根构型发生显著变化，总根长、侧根数、根表面积和根体积增加，长期处理下地上部生长相对稳定；早期酸磷酸酶和植酸酶活性升高，表明磷动员相关响应增强。转录组与代谢组分析揭示基因表达和代谢物积累的时间依赖性变化，包括PHT1、PAP12、GH3.1、PLDζ1、GDPD和PLA1的改变。色氨酸代谢、半胱氨酸与甲硫氨酸代谢、甘油磷脂代谢和植物激素信号转导被确定为主要响应途径。总之，这些发现表明根系发育、磷动员相关过程、膜脂代谢、硫代谢和激素相关信号与外源IAA存在下苜蓿对低磷的响应有关。
