人参叶片衰老过程中人参皂苷积累的时间动态与调控：基于UPLC-HRMS代谢组与转录组的整合分析

《Plants》：Temporal Dynamics and Regulatory Control of Ginsenoside Accumulation During Leaf Senescence in Panax ginseng

　　

人参皂苷（ginsenosides）是人参（Panax ginseng C.A. Meyer）的主要生物活性三萜皂苷，在叶片中大量积累。其构成一个包含300多种结构不同化合物的大次生代谢物家族，但叶片衰老如何重塑其积累仍不清楚。研究人员通过超高效液相色谱–高分辨质谱（UPLC-HRMS）分析了跨越自然衰老进程（T₁–T₃）和霜终止阶段（T₄）的四个采样阶段的376个皂苷代谢物特征，采用广义估计方程（GEE）对自然衰老梯度（T₁–T₃）的时间效应建模，将霜终止阶段T₄作为独立终点分析，经组织学验证衰老，并以基于重复的统计检验分析T₁–T₃转录组。叶片衰老全局重塑皂苷代谢组，376个皂苷特征中326个（86.7%）在自然衰老梯度（T₁–T₃）显示显著时间效应。四个经典苷元类别——原人参二醇（PPD）、原人参三醇（PPT）、奥克梯隆（OT）和齐墩果酸（OA）——在自然衰老进程中遵循分化积累程序：持续积累（PPT、OT和PPD）对U形曲线（OA）。在霜终止阶段T₄，各类群急剧分化：PPD型丰度回落至T₁水平（丰度加权类水平T₄/T₃ = 0.63），而其他类别进一步升高——这是一种霜关联重塑模式，而非衰老程序的延续。霜关联PPD下降依赖糖基化：低糖基化化合物比高糖基化化合物下降更剧烈（逐化合物比值T₄/T₃中位数0.51对1.04；单侧置换p = 1.2 × 10^?4），并与T₃时β-葡萄糖苷酶（β-glucosidase）基因诱导（探索性和假设水平，基于两个合格T₃重复，需更大重复数验证）及组织学验证的液泡塌陷一致，符合去糖基化相关分解代谢机制。衰老信号与差异激素响应（水杨酸（SA）激活，茉莉酸（JA）和乙烯减弱）、Ca²⁺信号组分广泛抑制及WRKY/NAC中心转录控制一致，联合分析提示二分皂苷–基因相关模式。总之，这些结果建立了人参叶片衰老过程中人参皂苷积累的时间图谱，刻画了皂苷次生代谢与衰老程序的类特异性耦合，并提供新视角：人参皂苷是嵌入人参生命周期生理活动的代谢物。

研究背景与问题。人参（Panax ginseng C.A. Meyer）主要活性成分人参皂苷（ginsenosides）为三萜皂苷，具抗肿瘤、免疫调节、神经保护和抗衰老等活性。按苷元分为原人参二醇（PPD）、原人参三醇（PPT）、奥克梯隆（OT）和齐墩果酸（OA）型，其生物合成经甲羟戊酸（MVA）和甲基赤藓醇磷酸（MEP）途径、鲨烯环氧化环化、细胞色素P450氧化及UDP-糖基转移酶（UGT）糖基化。皂苷稳态还受降解、转运和区室储存影响。叶片富含人参皂苷，是栽培副产物。叶片衰老受控，涉及基因重编程、细胞结构有序降解和营养再动员。拟南芥中衰老重构次生代谢已有系统记录，但药用植物叶片衰老与专用代谢物积累的耦合尚不清楚，人参中尚无研究追踪叶片衰老过程人参皂苷动态。前期分析显示类特异性：PPT和OT单调升高，PPD在晚衰峰后于枯萎期急降。因此研究人员提出问题：PPD终末期下降是否由衰老程序驱动？衰老信号如何差异调控不同皂苷类？

研究人员整合靶向代谢组、GEE时间建模、石蜡组织学、转录组及代谢组–转录组联合分析。结论为：自然衰老全局重塑皂苷代谢组，326/376特征出现显著时间效应；四类经典苷元分化积累，PPT、OT、PPD持续积累，OA呈U形；霜终止T₄时PPD回落至T₁水平而其他类上升；PPD下降具糖基化依赖性，与β-葡萄糖苷酶诱导和液泡塌陷一致，提示去糖基化分解，但为假设。衰老信号涉及SA激活、JA和乙烯减弱、Ca²⁺信号抑制及WRKY/NAC转录控制，联合分析显示二分相关模式。意义在于建立时间图谱，提出衰老–皂苷调控框架，提名WRKY和NAC087候选，并提示人参皂苷嵌入生命周期生理活动。

关键技术方法（≤250字）。三年生田间人参（Panax ginseng cv. Damaya）种植于吉林抚松，按表型采样：T₁红果后全绿、T₂叶尖早黄、T₃约50%黄化、T₄霜后萎蔫干枯。代谢组覆盖四阶段，每阶段3个生物学重复；转录组覆盖T₁–T₃，原计划每阶段6个重复，因RNA降解最终T₁=6、T₂=4、T₃=2。UPLC-HRMS靶向代谢组，MS/MS诊断碎片注释；石蜡切片400×；Illumina NovaSeq 6000 RNA-seq，fastp、STAR、RSEM，TPM定量，Welch t检验和FDR控制；GEE建模T₁–T₃时间效应，T₄作独立终点；K-means聚类；Spearman相关和Fisher富集。

研究结果。2.1 皂苷成分清单、检测稳定性与衰老驱动样本分离：通过UPLC-HRMS检测385个代谢特征，其中376个皂苷（139 PPT、87 PPD、32 OT、40 OA、78变异苷元其他）和9个非皂苷小分子；363/376在四阶段稳定检出；PCA显示衰老阶段主导代谢变异，样本按T₁–T₄分离。2.2 四类经典皂苷类群分化时间积累程序：GEE限制T₁–T₃，374/376收敛，326个（86.7%）显著时间效应；PPT/OT单调升高，其他皂苷T₁–T₃升高并在T₄激增，OA呈U形，PPD持续至T₃后T₄降37%至T₁；K-means分为慢升和陡升两组。2.3 霜关联PPD下降具糖基化依赖性：在86个完整PPD化合物中，60个T₄仍高于T₁，但丰度加权类总量下降；低糖基化PPD（1–2糖，n=19）T₄/T₃中位0.51，高糖基化（≥3糖，n=67）中位1.04，单侧置换p=1.2×10^?4；T₃β-葡萄糖苷酶基因诱导、UGT抑制，符合去糖基化分解假设，需酶活和蛋白验证。2.4 有序衰老降解的组织学验证：石蜡切片400×显示T₁健康，T₂细胞肿胀、叶绿体部分解体，T₃–T₄细胞失形、叶绿体广泛破坏、液泡边界不清；霜害叶呈皱缩、壁内陷和浓缩，区别于自然衰老，支持PCD样降解。2.5 晚衰全局转录重编程：T₃ vs T₁，56,128可测转录本中26,632个（47.4%）下调、4,740个（8.4%）上调，抑制偏倚5.6倍；因T₃仅两个合格重复，解释为指示性。2.6 衰老核心执行者顺序激活：metacaspase异构体差异表达，bVPE下调，SAG12和ENDO1强诱导，自噬基因ATG5、ATG10、ATG16抑制；与T₃进入不可逆降解阶段一致，但为假设水平。2.7 上游衰老信号通路差异响应：SA信号激活（PR1B1强诱导），JA信号抑制（MYC2沉默、AOS1下调），乙烯通路多数下调，Ca²⁺信号广泛抑制，MAPK中MPK3/MPK10下降而MKK9上升，鞘脂代谢重塑；T₃重复少，为假设水平。2.8 转录因子家族相反调控偏倚：WRKY、ERF、NAC上偏，ARF、HD-ZIP、bHLH抑制；WRKY22、WRKY9、WRKY72A、WRKY6及NAC087、NAC073诱导，MYC2沉默；探索性。2.9 联合分析揭示二分调控模式：在四类经典苷元中，PPT/PPD/OT积累与约80%相关基因负相关，OA积累与80.1%相关基因正相关；WRKY富集于PPT/PPD/OT正相关和OA负相关基因集（OR 3.61–4.04），ARF呈镜像富集；分析为探索性。2.10 衰老–皂苷调控耦合工作模型：整合提出五条信号通道经MAPK/Ca²⁺枢纽重塑WRKY/NAC诱导和ARF/bHLH/MYC2抑制，驱动衰老执行子和皂苷生物合成；T₃–T₄液泡塌陷暴露PPD给糖苷酶，优先降解低糖基化PPD，其他类积累；去糖基化-暴露为假设。

讨论部分总结：中心发现是PPD在自然衰老持续积累与衰老执行耦合，霜终点下降；T₄无转录组，液泡完整性和T₃分子起始提供假设框架。低糖基化PPD易完全去糖基化成不稳定苷元，高糖基化保留溶解性或经韧皮部输出至库，霜终点下降可能混合分解与源–库再动员；未测生物量、总产量和收获成本，实际采收窗口需农艺试验。人参衰老程序与拟南芥相似，SAG12、核酸酶、SA/JA、WRKY/NAC保守，NAC087和WRKY候选需验证。联合分析显示负相关为被动积累，正相关为主动共调控，WRKY可能具双重角色但需直接证据。变异苷元皂苷为苷元骨架多样化产物，因叶片液氮速冻可排除加工起源。源–库模型与降解假设互补；冷害、源–库再动员、氧化降解等不能完全排除。局限包括T₄无转录组、T₃重复少、代谢组与转录组时间分辨率不匹配、关联非因果、未测亚细胞定位和酶活、注释多为类别水平。

研究结论翻译：通过整合376个皂苷特征靶向代谢组、GEE时间建模、石蜡组织学和转录组，本研究建立人参叶片衰老中人参皂苷积累时间图谱并提出假设水平调控逻辑。靶向分析检测385个代谢特征，其中376个皂苷（139 PPT、87 PPD、32 OT、40 OA、78变异苷元其他）和9个非皂苷小分子；皂苷含1–5个糖基，中位2–3个，葡萄糖为主，其次芹菜糖/戊糖和鼠李糖；包括259个未修饰糖苷、81个乙酰化和47个丙二酰化衍生物。自然衰老中326/376（86.7%）有显著时间效应，四类经典苷元呈持续积累（PPT、OT、PPD）对U形（OA），T₄时PPD下降37%至T₁，其他类上升。PPD下降具糖基化依赖性（低糖基化中位T₄/T₃=0.51 vs 高糖基化1.04），与SAG12/ENDO1诱导和液泡塌陷耦合，符合液泡塌陷/糖苷酶暴露降解假设，仍需酶活和蛋白验证。衰老信号围绕差异激素响应、WRKY/NAC控制及二分相关模式，WRKY富集OR 3.61–4.04；因T₃两个重复，推断为探索性。这些发现首次提供药用植物叶片衰老重塑人参皂苷积累的整合图景，提名WRKY和NAC087为候选调节因子；针对BGLU、UGT和衰老执行基因的T₁–T₄ qPCR验证正在进行。