《Journal of Future Foods》：Analysis of metabolomics and transcriptomics of
Camellia drupifera seeds during maturation provides new insights into bioactive compound biosynthesis
编辑推荐：
本研究针对油茶种子成熟过程中生物活性化合物动态积累机制不清的问题，通过整合代谢组学、转录组学和生理指标分析，系统阐明了黄酮类、维生素E等关键活性成分的合成规律及其与抗氧化活性的关联，发现了果糖/甘露糖代谢在碳分配中的阶段性转换作用，并识别了跨通路调控枢纽基因，为油茶品质改良提供了重要理论依据。
在中国南方郁郁葱葱的山林中，生长着一种具有2300多年应用历史的珍贵植物——油茶。其中，石果油茶（Camellia drupifera）作为茶油 Camellia 物种的重要成员，以其种子中富含的多种生物活性化合物而备受关注。这些化合物不仅赋予茶油独特的营养价值，还具有抗氧化、抗炎等多种健康功效。然而，与广泛研究的油茶（C. oleifera）等物种相比，石果油茶种子发育过程中的生物活性成分积累规律及其调控机制仍是一个未被充分探索的领域。
目前的研究多集中于单一化合物的提取和评价，对种子发育过程中代谢物动态变化的系统性研究明显不足。特别是在不同成熟阶段，各类生物活性化合物如何协调积累，以及背后的转录调控网络如何运作，这些关键科学问题亟待解答。这种知识空白严重限制了我们对石果油茶种子品质形成机制的理解，也阻碍了其遗传改良和栽培优化。
为了揭开这一谜团，海南省农业科学院热带园艺研究所的研究团队在《Journal of Future Foods》上发表了最新研究成果。他们采用多组学整合分析策略，对石果油茶种子五个发育阶段进行了系统研究，旨在阐明生物活性化合物积累的动态规律及其调控机制，特别是揭示不同代谢通路之间的协同调控关系。
关键技术方法概述
研究团队收集了石果油茶五个发育阶段（2023年5月5日至10月22日）的种子样本，每个阶段设三个生物学重复。通过超高效液相色谱-电喷雾电离串联质谱（UPLC-ESI-MS/MS）进行广泛靶向代谢组学分析，共鉴定1559种代谢物。同时使用Illumina NovaSeq 6000平台进行转录组测序，以油茶"南荣油茶"品种基因组为参考，识别5971个差异表达基因。生理指标包括黄酮含量（FC）、维生素E（VE）、类胡萝卜素含量（CC）和多糖含量（PC2）等，并测定ABTS、DPPH、FRAP三种抗氧化活性。
发育过程中生物活性成分与抗氧化活性变化
研究发现，石果油茶种子在不同发育阶段表现出独特的生物活性成分谱和抗氧化能力。黄酮含量、维生素E和类胡萝卜素等在S2阶段最低，而在S2到S3阶段显著增加，这一时期种子处于液态胚乳阶段，代谢活跃。从S3到S5阶段，这些成分含量下降并趋于稳定，这与茶油种子8月至10月的主要油脂积累期相符。
相关性分析显示，维生素E与黄酮含量呈显著正相关（|r|=0.74），而类胡萝卜素含量、黄酮含量、维生素E和多糖含量这四个代表性指标与抗氧化能力密切相关，被选作后续分析的关键指标。
代谢组学揭示发育阶段特异性代谢特征
代谢组学分析共鉴定1559种代谢物，其中1382个为差异积累代谢物（DAMs）。黄酮类化合物占比最大（24.9%，344种），其次是酚酸（15%，208种）和氨基酸及其衍生物（12.8%，177种）。主成分分析（PCA）显示，不同发育阶段的样本明显分离，S3-S5阶段聚在一起，而S1和S2阶段形成独立簇群，反映了早期和晚期发育阶段之间清晰的生化转变。
KEGG富集分析表明，差异积累代谢物显著富集于槲皮素苷元、山奈酚苷元、黄酮苷元等生物合成途径，其中黄酮类生物合成途径在所有比较组中均被富集，突出了其在石果油茶种子代谢中的主要作用。
转录组动态揭示关键代谢通路调控机制
转录组分析发现，不同发育阶段间存在大量差异表达基因（DEGs），其中5971个基因为四个比较组所共有。GO富集分析显示这些基因与多种生物过程、分子功能和细胞组分相关。KEGG通路分析进一步表明，差异表达基因显著富集于淀粉和蔗糖代谢、黄酮生物合成、丙酮酸代谢等关键代谢途径，这与代谢物的通路富集模式相一致。
研究还鉴定出1707个转录因子（TF）来自87个转录因子家族，其中bHLH和MYB家族各包含90多个成员。这些转录因子表现出阶段特异性表达模式，可能在碳分配调控中协调初级和次级代谢途径。
黄酮类生物合成通路的关键作用
通过整合代谢组和转录组数据，研究发现黄酮类生物合成在种子发育过程中起着核心作用。32个差异积累的黄酮类化合物被定位到已鉴定的生物合成通路中，其中17种化合物在早期阶段达到峰值，13种主要在后期积累，表明成熟过程中黄酮组成发生动态变化。
研究鉴定了26个参与黄酮生物合成及黄酮和黄酮醇生物合成通路的差异表达基因，包括细胞色素P450（CYP450）、查尔酮合成酶（CHS）、查尔酮异构酶（CHI）、黄烷酮3-双加氧酶等编码基因。这些基因在早期阶段（S1-S2）高表达，有助于黄酮衍生物的强劲积累，增强抗氧化能力并保护发育中的组织。
碳分配的两阶段模型：从保护到储存
研究提出了一个描述石果油茶种子碳重新分配的发育模型。在早期阶段（S1-S2），甘露糖-1-磷酸鸟苷酰转移酶高表达，同时GDP-甘露糖、D-甘露糖和D-山梨醇水平升高，这些糖可能作为黄酮糖基化的供体。与此同时，CHS、CHI、F3H和其他黄酮生物合成基因的强表达促进了黄酮衍生物的强劲积累。
从S3到S5，GDP-甘露糖和相关糖的水平下降，而参与糖酵解和脂肪酸生物合成的基因上调，表明碳流向乙酰辅酶A生产和油脂积累重定向。共表达分析进一步揭示了四个关键基因（augustus_masked-HiC_scaffold_10-processed-gene-1914.34等）与多个通路的代谢物相关，表明了这个代谢转变背后潜在的跨通路调控节点。
特别值得注意的是，果糖和甘露糖代谢通路似乎经历了从为黄酮糖基化提供糖供体到可能支持脂质生物合成的转变。在类胡萝卜素和α-亚麻酸代谢中也观察到阶段特异性变化，脱落酸（ABA）含量在S1阶段积累到高水平，随后急剧下降，然后在S4和S5阶段适度回升。
跨通路调控网络的发现
研究首次揭示了石果油茶种子发育中存在高度整合的转录调控景观，其中初级和次级代谢过程紧密相连，而非作为离散模块运作。研究发现，单个代谢通路内的编码酶基因不仅与该通路的代谢物强相关，还与其他通路的关键代谢物相关。
相反，单个代谢物，如柚皮素和山奈酚-3-O-桑布双苷，与多个代谢模块的基因显著相关，包括类胡萝卜素、脂质和糖代谢。这些跨通路关联表明存在中央调控枢纽，在种子成熟过程中协调多个通路间的代谢流。
研究结论与意义
这项研究系统揭示了石果油茶种子成熟遵循两阶段策略：早期阶段（S1-S2）糖类优先导向黄酮糖基化及相关保护通路；后期阶段（S3-S5）碳流主要转向脂质生物合成和储存。这一转变由跨通路调控枢纽协调，包括参与糖、脂质和类胡萝卜素代谢的关键基因。
研究结果不仅阐明了石果油茶种子发育的协调代谢逻辑，还识别了 orchestrate 代谢重编程的候选调控枢纽。这些发现为理解种子发育过程中的碳分配提供了机制性见解，为提高茶油 Camellia 种子的营养品质和利用潜力提供了候选分子靶点和管理策略。
该研究的创新之处在于首次通过多组学整合分析，系统揭示了石果油茶种子发育过程中生物活性成分积累的动态规律及其转录调控网络，特别是发现了代谢通路间的协同调控关系，提出了碳分配的两阶段模型，为同类研究提供了重要参考。
