极端盐胁迫下马齿苋抗氧化网络的系统发育基因组学、时空转录动态及顺式调控架构解析

《Plant Science》：Phylogenomics, Spatiotemporal Transcriptional Dynamics and Cis-Regulatory Architecture Underlying the Portulaca oleracea Antioxidant Network Under Extreme Salinity

　　

• 在马齿苋（Portulaca oleracea）中绘制了跨13个基因类别的436个抗氧化基因座。• 根执行广泛的活性氧（reactive oxygen species, ROS）清除；成熟叶片聚焦于质外体（apoplastic）的III类过氧化物酶（class III peroxidase, PRX）屏障。• 急性胁迫下的幼苗依赖谷胱甘肽过氧化物酶（glutathione peroxidase, GPX）/谷胱甘肽还原酶（glutathione reductase, GR）途径。• 慢性盐胁迫下的成年植株部署独特且精简的抗氧化谱。• 预测到7个与系统性氧化还原协调相关的胁迫专属候选枢纽调控因子。• 盐诱导型过氧化氢酶（catalase, CAT）旁系同源基因具有近端富含GCC的AP2/ERF结合区块。

研究背景与问题：全球土壤盐渍化对农业可持续性构成持续威胁，盐胁迫通过先渗透冲击、后细胞内离子毒害的两阶段过程影响植物，Na⁺与Cl^-积累可抑制酶促反应、破坏光合反应中心并造成电子传递瓶颈，进而诱发活性氧（reactive oxygen species, ROS）爆发。耐盐植物依赖多基因调控的抗氧化防御维持氧化还原稳态，但马齿苋中协调该防御的分子架构仍不清楚。已有生化研究记录了超氧化物歧化酶（superoxide dismutase, SOD）、过氧化氢酶（catalase, CAT）、抗坏血酸过氧化物酶（ascorbate peroxidase, APX）、III类过氧化物酶（class III peroxidase, PRX）、谷胱甘肽还原酶（glutathione reductase, GR）等活性变化，单组学和多组学也勾勒了马齿苋胁迫响应轮廓，但尚缺乏针对完整抗氧化机器、跨器官和发育转变的靶向转录组分析，酶清除剂的时空表达动力学及调控网络架构亦未表征。为此，研究人员提出马齿苋的氧化还原韧性由器官特异和阶段依赖的转录策略介导：幼叶以快速抗氧化屏蔽应对急性渗透冲击，成年根则通过网络级酶促响应维持长期盐胁迫下的代谢稳态。该研究的意义在于将系统发育基因组学、时空转录动态、基因调控网络（gene regulatory network, GRN）与顺式调控架构整合，为解析极端盐胁迫下抗氧化网络提供高分辨率框架，并给出候选位点和顺式元件。
主要关键技术方法：研究人员使用马齿苋B1材料（Embrapa Agroenergia种质），采用两个队列：Silva et al.的幼年急性盐胁迫队列（1和4 DAT）与Rodrigues Neto et al.的成年慢性盐胁迫队列（12 DAT），分析叶和根组织。方法包括高通量RNA-Seq与qRT-PCR验证；基于参考基因组ASM2958216v1和PoDE数据库进行注释与表达分析；以Pfam/HMM、催化基序和OrthoFinder跨154个植物物种开展系统发育基因组学；用GENIE3构建条件特异GRN；用TargetP-2.0与DeepLoc-2.1进行亚细胞定位预测；并对CAT家族进行基因结构、SWISS-MODEL三维建模、启动子MAFFT比对及FIMO/JASPAR 2026顺式基序分析。
研究结果：3.1 抗氧化防御基因座的挖掘与丰度分析：通过关键词与EC号检索、结构域和催化基序筛选，研究人员在13个类别中鉴定436个抗氧化基因座，包括6个CAT、11个SOD、26个APX、196个PRX、4个DHAR、9个MDHAR、11个GPX、4个GR、67个GST、40个TRX和37个GRX；盐胁迫下不同家族呈现器官和发育阶段差异表达。3.2 成年植株组织中的差异性高幅抗氧化响应：以FC≥5.0且FDR≤0.05筛选12天盐胁迫后的高幅响应，根中18个基因达阈值，叶中仅7个基因高度上调；根呈现多通路响应，GRX3转录本升高50.24倍，叶中单个PRX转录本升高110.24倍并预测具N端信号肽，提示质外体H₂O₂清除。3.3 核心抗氧化候选基因的qRT-PCR验证：对6个候选基因在成年叶和根中进行qRT-PCR，12个组织—处理对比与RNA-Seq方向一致；叶中PRX-1降至0.02倍、CAT降至0.17倍，根中PRX-1和PRX-2分别升高2.51倍和2.18倍。3.4 抗氧化防御基因的系统发育基因组学分析：OrthoFinder分析显示PRX分布于43个直系同源组，GST为20个、TRX为19个、GRX为13个；CAT仅2个、SOD1为4个、SOD2为3个、CCS为1个，APX为11个；若干多基因直系同源组在成年根中协同激活。3.5 抗氧化基因座的调控网络建模：以140个差异表达抗氧化基因为终端节点构建条件特异GRN，UpSet分析显示179个转录因子（transcription factor, TF）为成年根基线调控因子、42个为成年叶基线调控因子；根中27个胁迫专属调控因子以M型MADS和LBD为主，叶中10个以B3为主；7个胁迫专属候选枢纽调控因子分属B3、LBD、C3H、ERF和bHLH五个家族。3.6 CAT基因家族的基因组结构、三维建模、进化一致性与顺式调控分化：6个CAT位于6条染色体并形成3个转录偶联对；CDS长度严格为1479–1488 bp，而总基因长度4,682–8,797 bp，差异主要来自内含子；三维建模显示89.02%–93.50%序列同一性及保守Fe³⁺-heme同源四聚体架构；启动子邻接树恢复同样3个姐妹对并具高支持率；FIMO扫描在盐诱导偶联对1近端启动子发现富含GCC区块，分别位于?147至?125 bp和?134至?112 bp，重叠14个AP2/ERF矩阵，而偶联对2和3近端区域无此完整特征。
讨论总结：作者将盐胁迫下ROS来源与马齿苋抗氧化响应置于器官分工框架中：根作为盐进入的主要位点，部署SOD2、PRX、GST、TRX/GRX等多层级缓冲以维持组织活力；成熟叶则以单个预测定位于质外体的PRX和稳定CAT蛋白池保护光合装置。幼年组织在急性胁迫下优先激活GPX/GR启动，成年组织在慢性盐胁迫下转向精简稳态网络。转录组与既往蛋白/代谢数据总体一致，但CAT呈现叶中转录本下降而蛋白稳定的转录—蛋白不对称。进化分析显示PRX、GST、TRX、GRX为扩张家族，而CAT、SOD、CCS等核心酶受较强约束；GRN预测出7个胁迫专属候选枢纽调控因子。CAT偶联对1具盐诱导性并拥有近端GCC-rich AP2/ERF结合块，偶联对3不响应，提示重复后启动子亚功能化。研究结论翻译：本研究建立了盐胁迫下马齿苋抗氧化防御系统的高分辨率空间、时间和系统发育基因组图谱。对13个基因类别436个基因座的全基因组分析揭示了显著调控区室化。在长期高盐下，成年植株执行独特的器官特异分工：根部署多层级酶和硫氧还蛋白/谷氧还蛋白缓冲以管理原初离子毒性，而成熟叶将转录资源集中于靶向质外体的III类过氧化物酶屏障，同时维持既有过氧化氢酶稳定性。跨发育阶段和胁迫动力学，幼嫩组织在急性渗透转变中优先谷胱甘肽依赖（GR）启动，而慢性盐胁迫下成年叶通过由稳定过氧化氢酶池支持的精简稳态框架运作。此外，条件特异网络重建鉴定出7个胁迫专属候选枢纽转录因子，预测参与系统氧化还原调控，并得到盐诱导型过氧化氢酶旁系同源基因中AP2/ERF相关富含GCC近端启动子亚功能化的支持。这些优先位点和顺式调控元件为旨在增强敏感作物非生物胁迫抗性和氧化保护的定向遗传工程提供了具体分子候选。