大豆烟酰胺介导的盐碱性耐受性的转录组学响应及推定机制

《Metabolites》：Transcriptomic Responses and Putative Mechanisms of Nicotinamide-Mediated Saline–Alkaline Tolerance in Soybean

　　

目的：盐碱性胁迫（由中性盐与碱性盐共同组成）严重制约大豆产量。尽管烟酰胺可增强植物耐盐性，但其在大豆中最有效浓度及分子基础尚不清楚。本研究旨在系统探究外源烟酰胺对大豆盐碱性耐受性的调控效应，并鉴定盐胁迫下候选基因与推定机制。方法：以大豆品种‘合丰25’幼苗为材料，施加200 mmol/L混合盐胁迫，并分别用0（CK）、10、50、200 mg/L烟酰胺处理。测定生长参数、活性氧（ROS）水平、离子含量、光合色素、渗透调节物质及抗氧化酶活性；行转录组测序与RT-qPCR鉴定差异表达基因（DEGs）并作富集分析。结果：外源烟酰胺较盐胁迫CK显著改善大豆盐碱性耐受性，其中50 mg/L（SA2）生理改良最显著，提升株高、生物量、K+含量与抗氧化酶活性，降低Na+积累与ROS水平。转录组分析确定PM1为五组非最优比较共有基因；DEGs富集于ko03110（分子伴侣与折叠催化剂）、ko00199（细胞色素P450）、ko04141（内质网蛋白加工）。RT-qPCR证实PM1、P450基因（CYP74A1、CYP83D1）、热激蛋白（HSP23、HSP18）及ROS清除基因（NAC1、GSTU43）在非最优剂量（SA1/SA3）上调、在SA2回落至胁迫CK水平；而SA2特异性恢复胁迫抑制的WRKY57与LBD15表达。SA2转录谱与盐碱CK最相似且DEGs最少。结论：50 mg/L烟酰胺最有效缓解大豆盐碱性胁迫，伴随最少转录扰动；非最优浓度引发内质网（ER）蛋白加工、解毒及氧化应激相关广泛重编程，PM1可作为剂量依赖转录调控候选标记。

研究背景与立项依据：土壤盐渍化是全球性非生物限制因子，其中盐碱性胁迫（saline–alkaline stress，由NaCl、Na₂SO₄等中性盐与NaHCO₃、Na₂CO₃等碱性盐共存）除离子毒害与渗透失衡外，还抬升环境pH、破坏胞内pH稳态，对大豆（Glycine max L.）等敏感作物造成萌发期与苗期不可逆损伤。大豆因根系排Na⁺能力弱、碱胁迫下易缺铁、共生固氮系统脆弱而耐受性偏低。烟酰胺（nicotinamide）作为维生素B3活性形式，是烟酰胺腺嘌呤二核苷酸（NAD⁺）与烟酰胺腺嘌呤二核苷酸磷酸（NADP⁺）前体，参与能量代谢、氧化还原稳态及聚ADP核糖聚合酶（PARP）/沉默信息调节蛋白（SIRT）信号，但其在大豆盐碱性胁迫下的最适浓度与转录层面机制尚未系统阐明。该研究发表于《Metabolites》，以黑龙江盐碱区主栽品种‘合丰25’为材料，试图填补最优剂量筛选、大豆专属转录景观刻画及分子通路推定的三重空白。

主要技术方法：研究人员采用完全随机设计盆栽试验，以V1期‘合丰25’幼苗接受200 mmol/L混合盐（摩尔比NaCl:NaHCO₃:Na₂CO₃:Na₂SO₄=1:1:9:9，pH≈8.6）胁迫，于胁迫第7天叶面喷施10 mg/L（SA1）、50 mg/L（SA2）、200 mg/L（SA3）烟酰胺，以等量去离子水为CK，每处理3次生物学重复；综合生理上依托火焰光度法测Na⁺/K⁺、碘化钾法测H₂O₂、羟胺氧化法测O₂^?·、乙醇萃取测光合色素、蒽酮/酸茚三酮法测渗透物、酶偶联DCPIP还原法测NAD⁺/NADH/NADP⁺/NADPH；转录层面对叶片混池总RNA行Illumina NovaSeq 6000测序（比对Glycine max Wm82.a2.v1），DESeq2筛|log₂FC|>1且p.adj<0.05的差异表达基因（DEGs），clusterProfiler行GO/KEGG富集，RT-qPCR用GAPDH归一化2^?ΔΔCt法验证。

研究结果

3.1 烟酰胺处理改善大豆盐碱耐受性：生长参数（株高、根长、鲜干重）呈先升后降趋势，SA2最优；叶中NAD⁺、NADH、NADP⁺、NADPH含量及NAD⁺/NADH比显著升高，NADP⁺/NADPH比下降；K⁺上升、Na⁺与Na⁺/K⁺比下降，SA1与SA3多数指标无显著差异，证明促进效应不随浓度单调增强而存在最适态。

3.2 烟酰胺降低ROS积累：SA1–SA3均显著降低O₂^?·与H₂O₂含量，与离子稳态改善协同。

3.3 光合色素提升：叶绿素a/b及类胡萝卜素在三类烟酰胺处理下均高于CK，但SA1与SA3间无差异，呈饱和而非连续剂量效应；综合生理仍以SA2最优。

3.4 渗透调节物与抗氧化酶增强：可溶性蛋白、可溶性糖、脯氨酸及超氧化物歧化酶（SOD）、过氧化氢酶（CAT）、过氧化物酶（POD）、抗坏血酸过氧化物酶（APX）活性均以SA2最高，SA3优于或持平SA1/CK。

3.5 DEGs鉴定：SA2 vs CK仅5个DEGs，SA1 vs CK为445个、SA3 vs CK为553个；排除SA2 vs CK后在五组非最优比较交集中获唯一共享基因PM1（LEA_1亚族），RT-qPCR示其在SA1/SA3显著诱导、SA2回落至胁迫CK水平；XTH1在SA1上调、SA2压回CK水平；WRKY57与LBD15胁迫下被抑，仅SA2相对CK去抑制。

3.6 GO与KEGG富集：非最优组DEGs富集ko03110、ko04141、ko00199及响应刺激、氧化还原酶活性等GO项；SA2 vs CK几乎无通路富集。CYP74A1、CYP83D1、HSP23、HSP18、NAC1、GSTU43在SA1/SA3上调、SA2与CK持平，印证非最优浓度触发ER蛋白加工、P450解毒与氧化应激转录重编程，SA2转录扰动最小。

讨论部分总结：研究人员指出DEGs数量反映转录扰动程度而非缓解效力本身，SA2少DEGs贴近胁迫CK谱，与玉米/大豆干旱研究中适度胁迫少DEGs现象一致；PM1作为剂量偏离候选标记需功能验证，其LEA蛋白属性支持离子螯合与膜稳定推测但非直接证据；NAD⁺/NADH与NADP⁺/NADPH比值偏移间接提示谷胱甘肽-抗坏血酸循环与ER redox稳态可能受益，但缺乏谷胱甘肽、UPR标记直接测量；细胞色素P450富集仅示解毒通路可能涉及。作者明确承认缺无胁迫对照与烟酰胺单施对照，无法断言“恢复至常态”，结论严格限定于盐碱CK比较框架内；单内参GAPDH、单品种、温室条件限制外推。

结论翻译：本研究证明外源烟酰胺增强大豆盐碱性耐受性，50 mg/L在全部测试浓度中生理改良最显著。转录组显示50 mg/L仅诱发有限转录扰动（5个DEGs），而亚最适与超最适浓度引发数百DEGs，提示非最优浓度触发应激通路广泛重编程。分层Venn分析确定PM1为烟酰胺剂量偏离的候选转录标记，其表达仅在50 mg/L时与盐碱CK持平。非最优比较DEGs的KEGG富集指向ko03110、ko00199、ko04141，推定性表明非最优浓度下ER蛋白加工、解毒及氧化应激相关转录响应，需实验验证。研究为受控盐碱条件下大豆烟酰胺剂量依赖生理与转录响应提供基础，田间适用性待评。