五生长年限三叶青块根相对代谢物丰度的UHPLC–Orbitrap-MS非靶向代谢组学解析

《Metabolites》：Untargeted Metabolomics Profiling of Growth-Year-Associated Changes in Relative Metabolite Abundance in Tetrastigma hemsleyanum Tuberous Roots

　　

背景/目的：三叶青（Tetrastigma hemsleyanum）是一种传统药用植物，其块根富含黄酮类和酚酸类成分，但这些成分的相对丰度如何随生长年限变化仍研究不足。方法：采用基于UHPLC–Orbitrap-MS的非靶向代谢组学分析五个生长年限（S1–S5；每年三株独立生长植株）的块根。注释依据代谢组学标准倡议（Metabolomics Standards Initiative，MSI）标准分级：349个整理条目中，252个为Level 2，56个为Level 3，因未分析真实标准品，无条目达到Level 1。结果：主成分分析（PCA）显示非线性分化，第一年根最独特。结合变量投影重要性（VIP）、倍数变化和t检验标准，鉴定出相邻生长年限间264个差异条目，经Benjamini–Hochberg（BH）校正后均保持显著（q < 0.01）。黄酮醇相关特征（如槲皮素相关信号、杨梅素和芦丁）在S1相对丰富；黄酮C-糖苷、黄烷-3-醇和酚酸特征（如牡荆素、儿茶素和阿魏酸）在S3；氨基酸、有机酸和核苷酸特征在S4。46个差异注释的京都基因与基因组百科全书（KEGG）映射提示黄酮类和氨基酸相关通路。结论：鉴于样本量小、相对峰面积数据和推定注释，这些结果具有假设生成性；它们识别出用于靶向定量的候选阶段和成分，而非最佳采收年限。

研究背景与问题：三叶青（Tetrastigma hemsleyanum Diels et Gilg）为葡萄科多年生攀缘藤本，块根是常用药用部位，传统用于清热解毒、活血祛风化痰等，现代药理研究显示其提取物及活性部位具有免疫调节、抗病毒、抗肿瘤、抗炎、抗氧化、解热和保肝等作用，黄酮类和酚酸类等次生代谢物被认为是重要物质基础。生长年限可显著影响药用植物次生代谢物积累，但三叶青大规模栽培历史较短、基础研究有限，其块根特征活性代谢物随发育年份的积累模式仍不清楚；已有代谢组学研究多集中于地理起源判别和机器学习辅助发现潜在标志物，发育研究虽开始描述块根形成相关分子变化，但缺乏对已注释代谢物类别跨连续栽培年限变化的系统信息，且数据库匹配可能产生假阳性，需明确报告注释置信度。因此，研究人员开展五生长年限块根非靶向代谢组学研究，以描述主要注释成分类相对丰度的年限相关变化、筛选差异特征并进行通路层面解读。该研究结论具有假设生成性，识别候选发育阶段和成分供后续真实标准品确认、靶向定量及与生物量、产量和药理质量指标整合，而非确定最佳采收年限。

关键技术方法（不超过250字）：样本采自中国浙江省丽水市毛竹林下半野生栽培基地的绿藤型三叶青块根，按块根特征划分S1–S5五个连续生长年限，每年限3株独立生长植株为3个生物学重复（n = 3，共15份）。采用UHPLC–Orbitrap-MS非靶向代谢组学，正负离子模式分别采集；Compound Discoverer 3.3进行峰检测、对齐及MS1/MS2数据库匹配，并按MSI标准分级。统计采用PCA、正交偏最小二乘判别分析（OPLS-DA）及单组分偏最小二乘判别分析（PLS-DA）敏感性分析，以VIP > 1、|log₂FC| > 1和p < 0.05联合筛选差异注释条目，并进行BH校正和KEGG过代表分析。

研究结果：
数据质量与方法可靠性：通过QC样本总离子流图重叠及PCA聚类评价，正负模式保留时间和峰面积一致，QC紧密聚类，满足验收标准；各生长年限内三株间399候选行中位变异系数（CV）分别为5.4%、5.1%、3.9%、4.6%、4.8%，1.5–4.3%的行CV > 20%，这些为生物与分析合并变异而非QC精密度。结论为分析方法稳定，但样本来自单一栽培群体。

五个生长年限间的整体代谢分化：通过PCA分析349行统计矩阵，五类呈分离趋势，组内三重复紧密；正负合并PC1/PC2解释42.1%/33.2%，正模式40.6%/33.1%，负模式45.2%/33.2%；S1最独特，分布非单调。结论为检测化学谱随年限呈非线性变化，仅描述所采样植株。

监督判别分析与模型验证：通过相邻年限OPLS-DA及200次置换，并另做单组分PLS-DA留一法及穷举平衡标签精确置换；OPLS-DA得分图相邻组分离，R²X为0.942、0.921、0.934、0.856，R²Y均为1.000，Q²为1.000、0.999、0.997、0.969；PLS-DA留一分类准确率1.00，交叉验证Q²为0.999、0.998、0.995、0.974，但精确置换p = 0.10。结论为近完美统计可能过拟合，不能作为预测泛化验证，VIP仅作筛选标准之一。

代谢物谱与分类：通过Compound Discoverer导出399候选行，经人工置信度分级、合理性审查和模式内去重得349整理条目（215 POS、106 NEG、28 NEG + POS合并）；252个Level 2（189高置信和63 probable）、56个Level 3、40个Level 4/未鉴定、1个待复核，无Level 1；303个可分类，包括氨基酸及衍生物33、黄酮类31、生物碱22、糖和醇22、酚酸22、脂质及类脂17、有机酸15、核苷酸及衍生物12、维生素4。保留时间检查发现山柰酚、异鼠李素、槲皮素苷元相关信号与糖苷共洗脱或正负模式保留时间差2.34 min，可能为源内碎片或异构体，称为苷元相关信号。结论为注释属推定性质。

阶段特异性差异注释条目筛选：通过VIP > 1、|log₂FC| > 1、p < 0.05联合筛选并BH校正；S2 vs S1为171个（63高/108低），S3 vs S2为173个（143高/30低），S4 vs S3为197个（43高/154低），S5 vs S4为123个（105高/18低）；合并去重264个；664个比较层面调用BH后q < 0.01（最大q = 5.5 × 10^?3，最大nominal p = 4.0 × 10^?3），四比较联合BH结果相同；敏感性分析中243/264（92.0%）仍满足q < 0.05且|log₂FC| > 1。结论为S3→S4低丰度条目最多，FDR不弥补n = 3，差异条目作探索候选。

主要注释成分类的生长年限相关相对丰度模式：通过主要类差异注释层次聚类热图（组均值z-score，Ward linkage）分析；黄酮类呈双峰，S1槲皮素糖苷、杨梅素、异鼠李素、芦丁相关最高；S3牡荆素、apigetrin等黄酮C-糖苷及儿茶素、表儿茶素、表没食子儿茶素等黄烷-3-醇富集；S5木犀草素、山柰酚、山柰酚-7-O-葡萄糖苷较高；酚酸阿魏酸、丁香酸、2,6-二甲氧基苯酚S3最高，水杨酸S1较高。氨基酸中犬尿喹啉酸和色氨酸衍生物S3，甲硫氨酸、酪氨酸、苏氨酸、脯氨酸S4，5-羟色胺/5-羟色氨酸S5，γ-氨基丁酸（GABA）/组氨酸/苯丙氨酸S2。有机酸琥珀酸、葡萄糖酸、顺乌头酸S4。脂质N-酰基乙醇胺S2，茉莉酸相关S3，13-HOTrE和12-氧代植物二烯酸S5。糖醇含环烯醚萜苷S3，蔗糖S4；核苷酸S4。结论为S1黄酮醇相关、S3多类次生代谢物共富集、S4氨基酸/有机酸/核苷酸增加；均为相对峰面积和推定注释。

KEGG通路富集分析：通过264差异条目中46个可映射KEGG、38个落入候选通路，对19条通路进行过代表分析（单侧超几何检验和BH校正，背景945个化合物）；氨酰-tRNA生物合成9个命中（p = 1.1 × 10^?4，FDR = 2.1 × 10^?3）、黄酮和黄酮醇生物合成8个命中（p = 5.8 × 10^?4，FDR = 5.5 × 10^?3）、黄酮类生物合成8个命中（p = 6.9 × 10^?3，FDR = 4.4 × 10^?2）达FDR < 0.05；丙氨酸、天冬氨酸和谷氨酸代谢、苯丙烷生物合成、维生素B6代谢呈富集趋势。Reactome靶向交叉检查21个关键注释仅蔗糖精确匹配。结论为提供有限通路背景，与聚类一致但不全面；背景非可检测化合物可能夸大显著性。

讨论部分总结：讨论指出PCA非线性分化且S1最独特；OPLS-DA近完美统计须结合n = 3及单独PLS-DA精确置换p = 0.10谨慎解释，仅支持探索性可视化。黄酮类呈双峰模式，S1黄酮醇相关相对富集，S3黄酮C-糖苷、黄烷-3-醇和酚酸相关富集，S4氨基酸、有机酸和核苷酸相关增加；S2 GABA和N-酰基乙醇胺相关较高，S3/S5茉莉酸和氧脂素相关注释可能与发育防御及次生代谢信号有关，但转录、酶、激素和通量未测，机制联系为假设。方法学上，研究以MSI置信度分级和合理性审查重新审视候选注释，无Level 1；跨模式去重按名称，保留时间检查发现苷元相关信号可能来自糖苷源内碎片或异构体，相关化合物需真实标准品和靶向MS/MS确认。局限包括单基地取样、未考察季节、生长年限按块根特征而非种植记录且未跨基地重复；峰面积为相对半定量非绝对浓度且未评估 exhaustive extraction efficiency；n = 3在高维设置下易过拟合且无外部验证；注释不确定；KEGG仅基于38个映射标识和通路背景；未测块根重量、产量、绝对成分含量和生物活性，不能直接推荐采收年限。

研究结论翻译：采用UHPLC–Orbitrap-MS非靶向代谢组学并结合MSI置信度分级，该研究描述了五个生长年限三叶青块根中阶段相关注释模式。整理数据集包含349个注释条目：252个Level 2、56个Level 3、40个Level 4/未鉴定和1个需进一步复核；无MSI Level 1确认。264个报告差异注释条目经BH校正后均显著，提示S1黄酮醇相关信号相对较高，S3黄酮C-糖苷、黄烷-3-醇和酚酸相关信号较高，S4氨基酸、有机酸和核苷酸相关信号较高。46个差异注释的KEGG映射提供有限通路层面背景，涉及黄酮类和氨基酸代谢。鉴于n = 3、单独敏感性分析精确置换结果不显著及通路覆盖不完整，这些结果应视为假设生成性；它们识别未来真实标准品确认、靶向定量及与生物量、产量和药理质量指标整合的候选发育阶段和成分类别，而非确定最佳采收年限。