枳椇Hovenia acerba是鼠李科Rhamnaceae枳椇属Hovenia高大落叶乔木，俗称拐枣。枳椇果实利用价值高，其果梗质量占果实总质量的90%。枳椇果梗肥大、肉质多汁、营养丰富，具很高的食用价值^[1−6]。光叶毛果枳椇H. trichocarpa var. robusta为枳椇属的变种^[7]。叶下面无毛或沿脉被疏柔毛，萼片和果实被锈色密绒毛。肥大果梗含丰富糖，可生食、酿酒、制醋和熬糖等^{[2, 8]}。枳椇果梗是主要食用部位，其滋味主要由非挥发性代谢物主导，不同种枳椇果梗常呈不同风味^[1]。卢孟飞等^[1]通过比较主栽枳椇品种果梗的营养成分差异发现：不同枳椇品种果梗营养成分和内在品质存在显著差异。因此，为更好地利用枳椇不同种及品种，需对其开展代谢物分析和评价。

广泛靶向代谢组学是一种高通量、高精度和广覆盖的代谢产物检测方法，在茶树Camellia sinensis^[9−12]、草莓Fragaria × ananassa^[13]、玉米Zea mays^[14]和黄纹早竹Phyllostachys violascens var. flavistriatus^[15]等植物的化学成分和品质特征研究中应用广泛。PENG等^[16]采用LC-ESI-QTOF-MS/MS技术对枳椇果梗进行测定，发现枳椇果梗的抗氧化活性与酚类、单宁酸总质量分数呈正相关，表明酚类物质可能是枳椇果梗抗氧化活性的主要成分，且首次鉴定出枳椇果梗提取物中含有低聚糖、有机酸、脂肪酸和黄酮类等影响风味的化合物。YANG等^[17]通过定量分析发现：枳椇果肉的主要营养成分包括糖类、蛋白质、脂质和黄酮等化合物，且糖类十分丰富。

目前，对枳椇研究集中于果梗的代谢物提取和功能食品的开发利用等^[17−18]，但对枳椇不同种的全面分析报道较少。诸多研究表明：代谢组学可对植物各器官内含物质进行准确解析，并筛选不同种或品种的特征化合物。因此，本研究以枳椇和光叶毛果枳椇果梗为试验材料，通过广泛靶向代谢组学鉴定2种枳椇果梗的代谢物差异，为枳椇种特征化合物的提取与利用及优良品种的选育提供参考。

1.   材料与方法

1.1   材料

材料选自浙江省丽水市白云山的枳椇和光叶毛果枳椇。枳椇浆果状核果，近球形，无毛(图1A)；光叶毛果枳椇浆果状核果，近球形，密被锈色或棕色绒毛(图1B)。2023年11月13日，在2种枳椇树的东、西、南、北4个方位，分别采集完全成熟的果梗，均匀混合后液氮速冻，−80 ℃保存。在液氮环境下用植物粉碎机粉碎，真空冷冻干燥后用于广泛靶向代谢组学测定。

图 1  2种枳椇属植物的果实表型

Figure 1  Fruit phenotypes of two Hovenia species

下载: 全尺寸图片 幻灯片

1.2   代谢组测定

称取50 mg样品粉末，溶解于1200 μL提取液(−20 ℃预冷的体积分数为70%的甲醇水内标提取液)；涡旋6次，每次间隔30 min，每次30 s；离心3 min (12 000 r·min⁻¹，4 ℃)，吸取上清液，用0.22 μm微孔滤膜过滤样品，保存于进样瓶，进行超高效液相色谱(UPLC，ExionLC™ AD)和串联质谱(MS/MS，QTRAP® 6500+)分析。

色谱条件：Agilent SB-C18 1.8 µm，2.1 mm×100.0 mm，流速0.35 mL·min⁻¹；柱温40 ℃；进样量2 μL。流动相条件：流动相A为超纯水(加入体积分数为0.1%的甲酸)，流动相B为乙腈(加入体积分数为0.1%的甲酸)。

1.3   统计分析

使用Analyst 1.6.3软件分析质谱数据中的代谢物，得到代谢物列表，包括代谢物名称和丰度。利用KEGG数据库(https://www.metaboanalyst.ca)对鉴定到的代谢物进行功能和分类注释。使用SIMCA 14.1软件对代谢物进行主成分分析(PCA)和正交偏最小二乘判别分析(OPLS-DA)。其中OPLS-DA 的S-plot图用于识别和解释模型中重要变量的关系，位于右上角和左下角的变量表示其差异显著，红色点通常表示分组贡献值(V_IP)≥1的变量，而绿色点表示V_IP＜1的变量。根据V_IP和组间变化的显著性筛选得到差异代谢物。再应用排列检验(permutation test)验证OPLS-DA模型可靠性，通过随机打乱样本的分组标签，重新建立和评估模型，以此来确定模型可靠性，R²为模型的解释率，代表了模型的可解释性，即模型中有多少信息能被解释。其中：R²X表示模型对预测变量(X矩阵)的解释程度；R²Y表示模型对响应变量(Y矩阵)的解释程度；Q²表示模型的预测能力，Q²＞0.5被认为是好模型。使用Excel 2019软件对数据进行标准化处理和t检验。接着通过t检验计算得到2组样品间的P值，筛选显著性差异代谢物。使用KEGG对差异代谢物进行通路分析，得到差异代谢通路的富集结果。在本研究中，V_IP为变量重要性投影值，表示代谢物对分组的贡献。log₂F_C表征差异倍数，为每个代谢物在样本组与对照组中3个重复均值的比值，正值表示样本组中代谢物相对含量高于对照组，负值则相反。P值用于判断实验结果的显著性。

数据标准化计算公式为：

其中：x 是原始数据，μ 是数据的均值，σ 是数据的标准差。

2.   结果与分析

2.1   2种枳椇果梗的代谢组成分总体分析

对2种枳椇果梗样品进行UPLC-MS/MS分析，共鉴定出13类1 810种代谢物(图2)，包括黄酮类307种(16.96%)、生物碱195种(10.77%)、脂类221种(12.21%)、有机酸65种(3.59%)、萜类191种(10.55%)、单宁13种(0.72%)、甾体1种(0.06%)、醌类18种(0.99%)、酚酸类207种(11.44%)、氨基酸及其衍生物146种(8.07%)、核苷酸及其衍生物54种(2.98%)、木脂素和香豆素130种(7.18%)以及其他代谢物262种(14.48%)。

图 2  2种枳椇属植物果梗代谢物分类

Figure 2  Metabolite classification of fruit pedicels in two Hovenia species　　　　　　　

下载: 全尺寸图片 幻灯片

2.2   2种枳椇果梗代谢组的主成分分析

为揭示2种枳椇果梗的代谢物组成差异，采用多维模式识别方法对代谢组数据进行多元统计分析。主成分分析得分表明：2组样品的3个重复良好，代谢组学数据质量可靠(图3A)。PC1和PC2分别解释了总变量的52.9%和14.9%，累计贡献率达67.8%。来自2组的样品表现出明显分层，并在第一主成分上分离显著，表明2组样品的代谢物在组间差异显著。载荷分析确定各代谢物与2个种之间存在相关性(图3B)。

图 3  2种枳椇属植物果梗PCA得分图(A)和载荷图(B)

Figure 3  PCA score diagram (A) and loading diagram (B) of fruit pedicels in two Hovenia species

下载: 全尺寸图片 幻灯片

OPLS-DA分析表明：第一主成分和第二主成分的贡献率分别为57.7%和13.7%，累计贡献率为71.4%，可更好地解释总体变量(图4A)。R²X=0.775，R²Y=0.998，Q²=0.972，其中Q²＞0.9表明模型可靠(图4B)。S-plot图表示主成分与代谢物的相关系数，可通过变量重要性投影值(V_IP)筛选差异代谢物(图4C)。

图 4  2种枳椇属植物果梗OPLS-DA图(A)、排序检验图(B)和S-plot图(C)

Figure 4  OPLS-DA score plot (A), ranked test plot (B) and S-plot (C) of fruit pedicels in two Hovenia species

下载: 全尺寸图片 幻灯片

2.3   2种枳椇果梗代谢物的鉴定与分析

为研究枳椇和光叶毛果枳椇果梗的代谢物差异，以V_IP＞1、|log₂F_C|＞1和P＜0.05为标准筛选出782种差异代谢物。与枳椇果梗相比，378种代谢物在光叶毛果枳椇果梗中显著上调，404种代谢物显著下调(P＜0.05，图5A)。为更好地了解差异代谢物种类，对筛选的差异代谢物进行标准化处理，并根据代谢物种类进行聚类分析(图5B)。结果表明：782种差异代谢物中，黄酮和酚酸类种类最多，分别为140和102种，其次为萜类95种、生物碱89种、木脂素和香豆素75种、脂类70种、氨基酸及其衍生物56种、核苷酸及其衍生物27种、有机酸22种、单宁4种、醌类4种和其他类98种。与枳椇果梗相比，单宁、醌类代谢物相对含量在光叶毛果枳椇果梗中显著下调(P＜0.05，图5C)。

图 5  2种枳椇属植物果梗差异代谢物火山图(A)、差异代谢物聚类热图(B)和差异代谢物散点图(C)

Figure 5  Differential metabolite volcano plot (A), differential metabolite clustering heat map (B) and differential metabolite scatter plot (C) of fruit pedicels in two Hovenia species

下载: 全尺寸图片 幻灯片

2.4   2种枳椇果梗差异代谢物通路分析

为进一步预测差异代谢物的主要生物学功能，将枳椇和光叶毛果枳椇果梗的差异代谢物注释到KEGG数据库，差异代谢物的代谢通路富集分析筛选到3条差异显著的代谢通路(图6)，每条通路中富集的代谢物如表1所示。差异显著的代谢物主要富集在苯丙烷类化合物、黄酮类化合物和异喹啉生物碱的生物合成途径(图7)。在这3条代谢通路中，黄酮类化合物生物合成途径富集的差异代谢物在枳椇果梗中的相对含量显著高于光叶毛果枳椇果梗(P＜0.05)；在苯丙烷类化合物生物合成途径中，除了芥子醛外，其他富集的差异代谢物在光叶毛果枳椇果梗中的相对含量均显著高于枳椇果梗(P＜0.05)；在异喹啉生物碱生物合成途径中，除了番荔枝碱外，其他富集的差异代谢物在光叶毛果枳椇果梗中的相对含量均显著高于枳椇果梗(P＜0.05)。

图 6  2种枳椇属植物果梗差异代谢物的富集通路气泡图

Figure 6  Enrichment pathway bubble diagram of fruit pedicels in two Hovenia species

下载: 全尺寸图片 幻灯片

表 1  2种枳椇属植物果梗差异代谢物显著富集通路

Table 1  Metabolites significantly enriched pathways of fruit pedicels in two Hovenia species

通路	通路名称	富集的代谢物名称
Map00940	苯丙烷类化合物生物合成	咖啡酰奎宁酸、对香豆酸、芥子醛、芥子酸、咖啡酸、阿魏酸
Map00941	黄酮类化合物生物合成	咖啡酰奎宁酸、柚皮素查耳酮、柚皮素、3, 4, 2’, 4’, 6’-五羟基查耳酮-4’-O-葡萄糖苷、樱桃苷、 　3, 4, 2’, 4’, 6’-五羟基查耳酮、3’, 4, 4’, 5, 7-五羟基黄烷、短叶松素、圣草酚、根皮苷、 　白矢车菊素
Map00950	异喹啉生物碱生物合成	甲基衡州乌药碱、N-甲基乌药碱、异波尔定碱、华青藤碱、乌药碱、原儿茶醛、 　对香豆酸、番荔枝碱

下载: 导出CSV 

| 显示表格

图 7  2种枳椇属植物果梗代谢差异中苯丙烷类化合物、黄酮类化合物和异喹啉生物碱的生物合成通路及相对含量差异

Figure 7  Biosynthesis pathways and contents of phenylpropanoid, flavonoid and isoquinoline alkaloid in fruit pedicels from two Hovenia species

下载: 全尺寸图片 幻灯片

3.   讨论

本研究对枳椇果梗和光叶毛果枳椇果梗进行广泛靶向代谢组学分析，以探究2种枳椇属植物果梗的代谢物差异，共鉴定出1 810种代谢物，其中枳椇果梗中化合物数量最多的三类分别是黄酮类、脂类和酚酸类化合物。KEGG代谢通路分析发现：2种枳椇果梗差异代谢物苯丙烷类化合物、黄酮类化合物和异喹啉生物碱的生物合成通路富集程度最高。前人研究表明：乌药碱、甲基衡州乌药碱和N-甲基乌药碱等生物碱化合物具有显著镇痛效果^[19−20]，在临床上被用于缓解疼痛，且具有抗炎、抗高血脂和抗肿瘤等多种药理作用^[21]。柚皮素、查尔酮等黄酮类化合物亦具有多种生物活性和药理作用，如较强抗氧化性和抗菌性^[22−25]，且查尔酮等黄酮类化合物还具有高甜度和低热量^[26]，可作为食品甜味剂。本研究发现：乌药碱、甲基衡州乌药碱和N-甲基乌药碱等生物碱化合物在光叶毛果枳椇果梗中显著上调，而查耳酮、柚皮素和2’, 3, 4, 4’, 6’-五羟基查耳酮等黄酮化合物在枳椇果梗中显著上调。因此，2种枳椇属植物果梗在潜在应用领域可能存在差异，需进行分类加工利用。

目前，枳椇的应用不再局限于果汁、果酒类产品^[27]，更多聚焦于果梗中多糖、黄酮和酚酸类化合物的提取和功能性研究^[28−31]。例如，黄酮类化合物具有抗菌性^[27]，枳椇花粉提取物和茎、叶提取物可抑制粪肠球菌Enterococcus faecalis介导的炎症和丙酸杆菌 Cutibacterium acnes引起的痤疮^[32−33]。本研究发现：枳椇果梗黄酮类化合物相对含量比光叶毛果枳椇果梗更高，更适合用于提取黄酮类化合物。然而，目前对于枳椇果梗黄酮类化合物抗菌性的研究及其提取物在抗菌上的应用尚欠缺。单宁、醌类具多样生物活性，可用于抗菌和抗病毒作用^[34−36]，但也存在潜在危害，例如，较高水平的单宁可与消化酶结合，影响消化功能^[37]，且可与蛋白质形成不溶性复合物，可能会干扰蛋白质和其他营养元素的吸收^[38]。部分醌类物质可能对肾脏、肝脏有损害^[39]，还被认为具有致癌作用^[40−41]。因此，在开发利用枳椇产品时，需进行充分的安全性评价。

本研究表明，苯丙烷类化合物、黄酮类化合物和异喹啉生物碱的生物合成途径是2种枳椇果梗中差异最显著的3条代谢通路，且多种中间产物也存在显著差异。这些差异代谢物与果梗感官差异的关系还需更多试验解析。

4.   结论

同一生态条件下，2种枳椇果梗的代谢谱存在显著差异，特别是具有药用价值和高甜度的生物碱和黄酮类代谢产物相对含量。这可能与2种枳椇果梗感官质量差异的形成有一定关系。本研究结果可为2种枳椇属植物的特征化合物提取与利用及优良品种的选育提供理论支撑。