吴茱萸品质变异的多维因素解析与层级质量评价研究

《Frontiers in Plant Science》：Dissection of quality variation and hierarchical quality evaluation of Evodiae Fructus under multidimensional factors

　　

摘要 引言： 吴茱萸（Evodiae Fructus）因多种混杂因素的影响，普遍存在显著的质量波动与异质性。为此，本研究开展多因素分析，旨在探索与其化学品质相关的潜在关键因素，为该药材的质量评价提供科学依据与理论支持。 方法： 首先，针对HPLC多组分定量数据集，研究人员采用Kruskal-Wallis H检验及随后的Dunn事后检验（经Benjamini-Hochberg FDR校正），并进一步构建以Gamma分布和对数连接函数指定的广义线性模型（GLM）以及随机森林回归模型，以量化植物来源、地理来源、商品规格和成熟度与吴茱萸品质之间的关联。其次，针对多维UPLC-MS/MS代谢组学图谱，研究人员应用主成分分析（PCA）、正交偏最小二乘判别分析（OPLS-DA）、Kruskal-Wallis H检验和差异倍数（FC）分析，在预设比较组间筛选差异代谢物，保留满足VIP > 1、|log2FC| > log2(1.5)且未校正p < 0.05的代谢物作为差异代谢物。研究人员利用基于GNPS的特征分子网络（FBMN）可视化相应的化学空间。随后，采用Pearson相关分析（|r| > 0.6，P < 0.05）和随机森林建模筛选候选核心标志物，并通过受试者工作特征（ROC）曲线分析（AUC > 0.8）进一步确认。最后，研究人员建立了比较性反向传播神经网络（BPNN）模型，其超参数通过使用10折交叉验证的网格搜索进行优化，分别以药典参考指标和筛选出的核心生物标志物作为输入，用于区分不同商品规格的吴茱萸，并在独立测试集上评估其预测性能。 结果： 研究发现商品规格与吴茱萸中活性成分的差异积累之间存在显著相关性。ROC曲线验证确认，2-羟基-4-甲氧基-3-(3'-甲基丁-2'-烯基)喹啉-4-酮和1-甲基-[(3E,6Z,9Z,12E)-十五碳-3,6,9,12-四烯-1-基]-4(1H)-喹诺酮这两种烷基喹诺酮类生物碱，可被视为吴茱萸质量评估的特异性潜在生物标志物。这两个候选标志物能够区分不同规格的吴茱萸。 结论： 综上，本研究厘清了商品规格与吴茱萸化学品质之间的独立及交互关联。研究人员获得的候选标志物使BPNN模型能够有效区分不同商品规格。这些发现为该药材的层级质量评估提供了有价值的参考。

吴茱萸品质变异的多维因素解析与层级质量评价研究——论文解读

一、研究背景与问题

吴茱萸（Evodiae Fructus）为芸香科植物吴茱萸（Evodia rutaecarpa）、石虎（Evodia rutaecarpa var. officinalis）或疏毛吴茱萸（Evodia rutaecarpa var. bodinieri）的干燥近成熟果实，具有温中止痛等功效，在中医临床广泛应用，是左金丸、华佗再造丸等经典方剂的重要组成药物。然而，受长期自然演化、环境条件和异地引种的综合影响，吴茱萸形成了丰富的变异类型，导致市售药材原植物混杂、不同产地样品质量差异显著。同时，现行标准体系中商品规格划分标准不一，市场上存在“同药不同标”的乱象，例如广西产“大花吴茱萸”果顶开裂的形态特征在现行标准中缺乏明确的分类定位，部分商家将其混入中花批次或以低价劣药出售，不仅阻碍了药材资源的合理利用，也加剧了药材质量的波动。当前《中国药典》对吴茱萸的质量控制主要以特征性生物碱和柠檬苦素含量为核心评价指标，但单一的质量控制评价体系难以全面客观地反映药材整体质量。现有质量评价研究多聚焦于地理来源或商品规格等单一因素，忽视了基于整体代谢组的药材差异特征，无法全面阐释吴茱萸的内在质量异质性。因此，建立系统全面的质量评价体系、识别特征性质量标志物，对于稳定吴茱萸质量、保障临床用药安全具有重要意义。本研究旨在基于非靶向代谢组学结合多元统计与化学计量学方法，从地理来源、成熟度、植物来源和商品规格四个维度系统解析吴茱萸化学成分的变异规律，量化各因素对化学轮廓差异的相对贡献，并筛选与不同分组相关的候选质量标志物。

二、研究设计与技术方法

研究人员共收集42批次吴茱萸果实样品，其中10个批次提供了不同成熟度的子样品，总计59个实验样品，来源于药材市场和种植基地，涵盖广西、湖南、江西、贵州等产地。所有样品经HPLC多指标定量分析（包括吴茱萸次碱、吴茱萸碱、柠檬苦素），所有近成熟期样品进一步经UPLC-Q-Exactive Orbitrap-HRMS非靶向代谢组学分析。主要关键技术方法包括：①基于HPLC定量数据，采用Kruskal-Wallis H检验结合Dunn事后检验、Gamma分布对数连接的广义线性模型（GLM）和随机森林回归模型量化各因素与品质的关联；②基于UPLC-MS/MS代谢组学数据，采用主成分分析（PCA）、正交偏最小二乘判别分析（OPLS-DA）、Kruskal-Wallis H检验和差异倍数分析筛选差异代谢物，利用GNPS特征分子网络（FBMN）可视化化学空间；③通过Pearson相关分析和随机森林建模筛选候选核心标志物，经ROC曲线分析（AUC > 0.8）确认；④建立反向传播神经网络（BPNN）模型，以药典指标和核心生物标志物为输入区分不同商品规格，超参数经网格搜索和10折交叉验证优化。

三、主要研究结果

3.1 多维分析探索影响吴茱萸品质的关键因素

3.1.1 基于单一维度的影响程度分析

Kruskal-Wallis H检验结果显示，四个影响因素均对标志物含量产生显著影响。吴茱萸碱和吴茱萸次碱合计含量总体趋势为：广西样品高于江西和湖南；E. rutaecarpa高于E. rutaecarpa var. officinalis；大花和中花规格高于小花规格；近成熟和过成熟样品高于未成熟样品。而柠檬苦素含量在不同产地间无明显差异，其变化趋势与生物碱合计含量相反。随机森林回归模型分析表明，果实成熟度和商品规格是柠檬苦素含量的最主要影响因素，而商品规格主导吴茱萸碱-吴茱萸次碱合计含量，证实商品规格是吴茱萸样品品质变异的关键来源。

3.1.2 Gamma分布对数连接广义线性模型分析

对吴茱萸碱和吴茱萸次碱合计含量，成熟度呈现主导主效应（Wald χ2 = 119.182），各二阶交互项均无统计学意义。分层分析显示，在近成熟阶段，小花样品生物碱积累显著低于大花样品（Exp(B) = 0.250，P < 0.001），中花样品也低于大花样品（Exp(B) = 0.482，P = 0.002）。对柠檬苦素含量，商品规格是最强影响因素（Wald χ2 = 27.82，P < 0.001），且与地理来源存在显著交互作用（P = 0.011），但效应方向在三个产地中一致：小花样品柠檬苦素含量均显著高于大花样品，其中湖南群体的规格梯度效应最强。

3.2 基于多维分析的吴茱萸差异代谢物筛选

对61个样品、690个代谢特征的合并数据集进行PCA和OPLS-DA分析，PCA模型前两个主成分解释94.5%的总变异，OPLS-DA模型（R2X = 0.895，R2Y = 0.992，Q2 = 0.989）能有效区分六个比较组。基于不同成熟度的分析显示，吴茱萸酰胺I和吴茱萸碱等特征性生物碱以及喹啉类生物碱在近成熟和过成熟样品中显著上调，而多巴胺、5-甲氧基-N,N-二甲基色胺等酰胺类化合物在未成熟果实中富集。基于不同植物来源的分析显示，E. rutaecarpa样品富集吲哚类和喹啉类生物碱，而E. rutaecarpa var. officinalis样品主要富集咔啉类生物碱和苯丙素类成分。基于不同商品规格的分析显示，中花规格E. rutaecarpa var. officinalis积累丰富的喹诺酮类生物碱，小花规格柠檬苦素含量较高；大花规格E. rutaecarpa富集喹诺酮类生物碱，中花规格则显著积累酚酸和黄酮类成分。基于不同地理来源的分析显示，广西产样品富集喹诺酮类生物碱，湖南产样品以黄酮类和氨基酸代谢物为特征，江西产样品中3-羟基-1H-喹啉-4-酮显著升高。

3.3 机器学习结合相关分析筛选核心代谢标志物

Pearson相关分析显示，吲哚类生物碱（吴茱萸酰胺I、二氢吴茱萸次碱等）、喹啉类生物碱（1-甲基-2-癸基-4(1H)-喹诺酮等）和黄酮类代谢物（异鼠李素）与吴茱萸碱和吴茱萸次碱合计含量显著正相关，而5-甲氧基-N,N-二甲基色胺、柠檬苦素和部分苯丙素类代谢物与该指标显著负相关。随机森林模型筛选出变量重要性评分前15的代谢物作为候选特征标志物，ROC曲线分析综合AUC结果表明，2-羟基-4-甲氧基-3-(3'-甲基丁-2'-烯基)喹啉-4-酮和1-甲基-[(3E,6Z,9Z,12E)-十五碳-3,6,9,12-四烯-1-基]-4(1H)-喹诺酮可作为吴茱萸的潜在生物标志物，两者联合应用可进一步提高对不同商品规格样品的区分准确度。BPNN模型验证了这两个候选标志物能有效区分不同商品规格的吴茱萸样品。

四、讨论与结论

综合随机森林模型、广义线性模型和聚类热图分析结果，商品规格对吴茱萸品质变异具有显著的独立效应，应作为该药材质量控制的关键依据。吴茱萸碱和吴茱萸次碱的合计含量主要随成熟度推进而增加，并与商品规格、地理来源和植物来源存在显著独立关联；而柠檬苦素的积累主要受商品规格控制，且该规格效应的强度受地理来源显著调节。生物碱和柠檬苦素成分对商品分级的差异性响应表明，应根据吴茱萸药材的具体化学特征实施针对性的采收期优化和商品分类标准制定。

研究结论： 本研究厘清了商品规格与吴茱萸化学品质之间的独立及交互关联。筛选获得的两种烷基喹诺酮类生物碱候选标志物使BPNN模型能够有效区分不同商品规格的吴茱萸。这些发现为该药材的层级质量评估提供了有价值的参考，为吴茱萸的质量评价、商品规格分级和标准化生产提供了实验证据和数据支持。