来源：BMC Genomics | 编译：DNA Lab Space
沙棘果实发育早期，哪些基因在悄悄塑造果实的质地？一项针对中国沙棘（Hippophae rhamnoides subsp. sinensis Rousi）的全基因组分析给出了线索。研究者从沙棘中鉴定出71个Ⅲ类过氧化物酶（POD）基因，通过表达谱、共表达网络和胁迫处理实验，锁定5个与木质素沉积密切相关的候选成员。这项工作发表在BMC Genomics上，为理解沙棘果实硬度形成和广谱抗逆性提供了分子层面的新视角。
研究背景
Ⅲ类过氧化物酶是一类含血红素的氧化还原酶，在植物中普遍参与木质素聚合、活性氧清除、细胞壁加固和胁迫响应等过程。木质素是植物次生细胞壁的关键组分，它的含量和分布直接决定果实质地——木质素过多，果实偏硬；过少，则果实软烂易损。沙棘作为西北干旱半干旱地区的重要经济林树种，果实富含维生素和黄酮类化合物，但果实软化快、货架期短的问题长期困扰产业发展。
然而，沙棘POD基因家族的系统鉴定和功能研究此前一直空白。哪些POD成员参与沙棘果实木质素合成？它们是否响应干旱、盐碱等逆境信号？这些问题的答案，关系到沙棘果实品质改良和抗逆品种选育的分子靶点选择。
研究方法
植物材料与取样
研究团队从青海省祁连县沙棘天然杂交区（38°15′N，100°16′E）采集中国沙棘果实样本。依据此前发表的采样标准，选取花后30天（Hrh_S2）、70天（Hrh_S4）和120天（Hrh_S6）三个发育阶段的果实，分别对应幼果期、膨大期和成熟期。凭证标本保存于西北师范大学植物标本馆（NWTC），编号NWTC-HR-2209，经孙坤教授鉴定确认物种身份。果实样品经液氮速冻后转入超低温冰箱保存，用于转录组测序、木质素含量测定、过氧化物酶活性检测和实时荧光定量PCR分析。
基因家族鉴定与系统发育分析
利用已知植物POD蛋白序列作为查询序列，对沙棘全基因组进行BLAST搜索和HMMER结构域筛选，共获得71个非冗余HrPOD基因。将这些基因的蛋白序列与拟南芥、水稻等物种的POD家族成员进行多序列比对，采用邻接法构建系统发育树，根据进化关系将71个HrPODs划分为7个亚组。随后用MEME工具和CDD数据库分析蛋白保守基序和结构域组成。
表达模式分析
利用RNA-seq数据计算71个HrPODs在三个果实发育阶段的FPKM值，绘制热图展示表达动态。同时，对果实样品进行木质素含量测定（采用Klason法或硫代乙酰基法）和过氧化物酶活性检测（愈创木酚氧化法），获得与转录组数据对应的生理指标。
功能富集与共表达分析
对71个HrPODs进行Gene Ontology（GO）和Kyoto Encyclopedia of Genes and Genomes（KEGG）功能富集分析，明确这些基因参与的生物学过程和代谢通路。将HrPODs表达量与木质素含量、POD活性进行Pearson相关性分析，以|r|>0.8且P<0.05为阈值筛选候选基因。
非生物胁迫处理
鉴于野生沙棘常分布于干旱和高盐碱环境，且胁迫从种子萌发到果实成熟全程存在，研究团队重点考察萌发期的胁迫响应。用不同浓度的PEG6000（0、5、10、15、20%）、NaCl（0、0.05、0.10、0.15、0.20 mol/L）和NaHCO₃（0、0.005、0.010、0.015、0.020 mol/L）溶液分别处理沙棘种子，模拟干旱、盐和碱胁迫，检测处理后的基因表达变化和木质素含量。
研究结果
71个HrPODs的系统鉴定与分类
研究共鉴定出71个非冗余HrPOD基因，这一数量与拟南芥（73个）和水稻（138个）的POD家族规模相当，表明沙棘POD家族在进化过程中经历了复杂的扩张与分化。系统发育分析将这71个基因清晰划分为7个亚组，各亚组内成员在保守基序组成上高度相似，进一步验证了分类的可靠性。蛋白结构分析显示，所有HrPODs均含有典型的植物过氧化物酶结构域，包含血红素结合位点和底物结合口袋的关键残基，这些高度保守的氨基酸位点是酶催化活性所必需的。
表达谱揭示发育期特异性调控
转录组数据显示，71个HrPODs在果实三个发育阶段的表达呈现明显差异。大部分基因（超过60%）在花后30天的幼果期（Hrh_S2）表达量最高，随着果实发育逐渐下降。这一表达趋势与木质素在果实发育早期的快速积累相吻合——幼果期细胞壁正在发生次生加厚，需要大量POD参与木质素单体聚合。少数基因在成熟期（Hrh_S6）出现表达峰值，可能与果实成熟过程中的乙烯响应或衰老相关氧化过程有关。
功能富集指向苯丙烷代谢通路
GO富集分析显示，HrPODs显著富集在“氧化应激响应”和“过氧化物酶活性”等条目，这与Ⅲ类POD清除H₂O₂、缓解氧化损伤的已知功能一致。KEGG富集分析则将这些基因主要定位到“苯丙烷生物合成”通路——木质素正是苯丙烷代谢途径的末端产物之一。两条富集结果相互印证：HrPODs既参与木质素聚合这一合成方向，又参与氧化还原平衡这一调控方向，暗示它们在果实发育中承担双重角色。
共表达分析筛选出5个关键候选基因
将71个HrPODs的表达量与果实木质素含量、POD活性进行Pearson相关分析，发现5个基因的表达模式与木质素沉积曲线显著正相关（相关系数r>0.9，P<0.01）。这5个候选基因分布在不同的系统发育亚组中，说明木质素合成的调控并非由单一进化分支负责，而是多个亚组成员协同参与。进一步分析发现，这5个基因在幼果期的表达量是成熟期的3-8倍，与木质素含量的时间变化趋势高度同步。
胁迫处理影响基因表达与木质素积累
萌发期胁迫实验显示，PEG6000、NaCl和NaHCO₃处理均能改变HrPODs的表达水平，且响应模式随胁迫类型和浓度而异。以20% PEG6000处理为例，部分候选基因表达量较对照上调2-5倍；而0.20 mol/L NaCl处理则使另一些基因表达下调。木质素含量检测同样呈现浓度依赖性变化：适度的干旱或盐胁迫（如10% PEG6000、0.10 mol/L NaCl）促进木质素积累，而高浓度胁迫（20% PEG6000、0.20 mol/L NaCl）则抑制木质素合成。这种非线性响应关系提示，POD介导的木质素沉积在轻度胁迫下可能作为防御机制增强细胞壁强度，而重度胁迫则导致代谢紊乱。
讨论与展望
这项研究首次系统完成了沙棘POD基因家族的全基因组鉴定，71个成员的规模、7个亚组的分类格局以及保守结构域特征，为后续功能研究奠定了框架。5个候选HrPODs的筛选尤其值得关注——它们同时满足“与木质素含量显著相关”和“响应多种非生物胁迫”两个条件，暗示这些基因可能是连接果实质地形成与逆境适应的分子节点。
研究选择萌发期进行胁迫处理，这一设计颇具巧思。种子萌发阶段是植物生活史中对环境最敏感的窗口期，沙棘在干旱盐碱地自然分布的特性，决定了其萌发期必然演化出高效的胁迫响应机制。该研究选择在萌发期而非营养生长期或生殖生长期开展胁迫实验，可能正是基于对这一生态适应性的考量。
当然，目前的工作仍以生物信息学分析和表达相关性为主，5个候选基因是否真正参与木质素合成、如何被上游转录因子调控，尚需通过异源表达、基因编辑和酶动力学实验加以验证。木质素含量的测定方法、POD活性的具体数值，原文摘要未详述具体参数。未来若能将候选基因在沙棘中进行稳定转化，观察果实硬度和木质素沉积的实际变化，将是对当前结论最直接的检验。
沙棘果实品质改良和抗逆育种，正需要这样从基因家族全景到候选靶点的系统性工作。71个POD基因的图谱已经绘就，接下来要做的，是让这些分子开关在田间地头真正转动起来。
参考来源
Zhao J, Li K, Zhao M, Jiang X, Liu X. Identification and characterization of the class III peroxidase gene family in sea buckthorn (Hippophae rhamnoides subsp. sinensis Rousi) and its potential role in lignin biosynthesis. BMC Genomics. PMID: 41572154.
PMID: 41572154