Page 48 - 《中国药房》2023年22期

表1 78批山楂叶样品相似度评价及GCA分析结果 直观呈现数据结果，成分含量越高则条带颜色越红，反
[12]
相对最优 相对最差 相对关 相对关联 相对最优 相对最差 相对关 相对关联 之越蓝 。本研究将4个主产地（河北省承德市、辽宁省
编号 相似度 编号 相似度
关联度 关联度 联度 度排序 关联度 关联度 联度 度排序 葫芦岛市、山东省临沂市、山西省运城市）的78批山楂叶
S1 0.950 0.668 0.539 0.553 15 S40 0.980 0.663 0.540 0.551 17
S2 0.983 0.616 0.635 0.492 39 S41 0.973 0.594 0.686 0.464 48 药材的 8 个共有峰峰面积导入 MetaboAnalyst 5.0 工具，
S3 0.871 0.816 0.449 0.645 1 S42 0.990 0.637 0.579 0.524 25 相关数据经筛选、归一化处理后，得到CA热图（图3A）。
S4 0.969 0.627 0.642 0.494 38 S43 0.993 0.606 0.636 0.488 41 由图3A可见，78批样品可分为2类，河北省承德市和辽
S5 0.974 0.627 0.642 0.494 38 S44 0.970 0.627 0.632 0.498 35
S6 0.983 0.646 0.554 0.538 22 S45 0.982 0.698 0.505 0.580 10 宁省葫芦岛市样品大致分为一类，山东省临沂市和山西
S7 0.996 0.760 0.437 0.635 2 S46 0.987 0.607 0.654 0.481 43 省运城市样品大致分为一类；结合热图的条带颜色可
S8 0.961 0.646 0.570 0.531 23 S47 0.956 0.631 0.593 0.516 27
S9 0.994 0.617 0.619 0.499 34 S48 0.951 0.544 0.874 0.384 63 知，河北省承德市和辽宁省葫芦岛市样品中峰1、3、5、6、
S10 0.998 0.706 0.482 0.595 8 S49 0.978 0.568 0.765 0.426 58 7、8 代表成分含量较高，山东省临沂市和山西省运城市
S11 0.954 0.587 0.741 0.442 54 S50 0.998 0.621 0.598 0.509 31 样品中这6个成分含量则较低。
S12 0.993 0.653 0.537 0.548 19 S51 0.988 0.617 0.636 0.492 39
S13 0.977 0.703 0.543 0.564 13 S52 0.983 0.648 0.575 0.530 24 （3）PCA：采用 SPSS 19.0 软件，以 78 批山楂叶药材
S14 0.992 0.649 0.556 0.539 21 S53 0.987 0.630 0.613 0.507 33 的 8 个共有峰峰面积为变量进行 PCA，并进行 KMO 检
S15 0.952 0.574 0.776 0.425 59 S54 0.973 0.637 0.598 0.516 27
S16 0.973 0.591 0.699 0.458 50 S55 0.978 0.568 0.766 0.426 58 验 及 Bartlett 球 形 检 验 。 结 果 显 示 ，KMO 度 量 值 为
[13]
S17 0.943 0.582 0.762 0.433 56 S56 0.974 0.587 0.718 0.450 52 0.616，Bartlett 球形检验的 P＜0.001，可进行 PCA 。将
S18 0.980 0.603 0.667 0.475 45 S57 0.939 0.652 0.558 0.539 21 8 个共有峰峰面积导入 MetaboAnalyst 5.0 工具，相关数
S19 0.988 0.630 0.606 0.510 30 S58 0.967 0.590 0.718 0.451 51
S20 0.990 0.634 0.597 0.515 28 S59 0.969 0.573 0.769 0.427 57 据经筛选、归一化处理后，得到PCA得分图（图3B）。由
S21 0.994 0.734 0.457 0.616 6 S60 0.980 0.644 0.596 0.519 26 图 3B 可见，河北省承德市和辽宁省葫芦岛市样品大致
S22 0.989 0.632 0.595 0.515 28 S61 0.941 0.566 0.791 0.417 60
S23 0.984 0.615 0.642 0.489 40 S62 0.979 0.578 0.734 0.440 55 可分为一类，山东省临沂市和山西省运城市样品大致可
S24 0.973 0.673 0.558 0.547 20 S63 0.982 0.613 0.645 0.487 42 分为一类，与CA热图结果基本一致。
S25 0.997 0.719 0.472 0.604 7 S64 0.984 0.607 0.656 0.481 43 （4）OPLS-DA：为进一步明确差异标志物，本研究将
S26 0.998 0.623 0.596 0.511 29 S65 0.983 0.562 0.789 0.416 61
S27 0.994 0.658 0.538 0.550 18 S66 0.954 0.544 0.865 0.386 62 78 批山楂叶药材的 8 个共有峰峰面积导入 SIMCA 14.1
S28 0.896 0.705 0.502 0.584 9 S67 0.943 0.544 0.880 0.382 64 软件，建立OPLS-DA模型。由模型验证参数可知，模型
S29 0.982 0.659 0.527 0.556 14 S68 0.984 0.610 0.622 0.495 37 2 2
S30 0.995 0.655 0.543 0.547 20 S69 0.978 0.630 0.610 0.508 32 的稳定性较好（R X＝0.843，R Y＝0.725）且交叉验证预测
2
2
[14]
S31 0.994 0.689 0.500 0.579 11 S70 0.981 0.617 0.619 0.499 34 能力较强（Q ＝0.532） 。置换检验的 R 拟合直线在 Y
S32 0.976 0.606 0.655 0.480 44 S71 0.982 0.616 0.636 0.492 39 坐标轴上的截距为0.049 6（＜0.3），说明所建模型可靠；
S33 0.987 0.595 0.674 0.469 47 S72 0.974 0.581 0.729 0.444 53
2
S34 0.987 0.595 0.672 0.470 46 S73 0.987 0.617 0.620 0.499 34 Q 拟合直线在 Y 坐标轴上的截距为－0.188 0（＜0.05），
S35 0.985 0.589 0.692 0.460 49 S74 0.972 0.563 0.789 0.416 61 说明所建模型不存在过度拟合，可用于分析样品的组间
S36 0.988 0.598 0.662 0.475 45 S75 0.982 0.609 0.635 0.489 40 [15]
S37 0.989 0.663 0.537 0.552 16 S76 0.935 0.747 0.460 0.619 5 差异 。提取 OPLS-DA 模型中 8 个共有峰的变量重要
S38 0.983 0.676 0.519 0.566 12 S77 0.945 0.758 0.449 0.628 3 性投影（variable importance in the projection，VIP）值，以
S39 0.973 0.625 0.632 0.497 36 S78 0.947 0.746 0.454 0.622 4 [16]
VIP 值＞1 且误差线不超过原点为标准 进行差异标志
（2）CA：CA 热图可对数据进行筛选、提取和降维， 物筛选，结果（图 4）显示，峰 1（绿原酸）、峰 2、峰 3（牡荆
通过横向聚类来反映各样品之间的关系，纵向聚类来反 素葡萄糖苷）、峰5的VIP值均大于1，其对应成分为差异
映各化学成分之间的关系，该图通过渐变的蓝红色条来 标志物，可能是影响山楂叶药材质量的主要化学成分。
河北省承德市 辽宁省葫芦岛市 山东省临沂市 山西省运城市
0.5

产地
峰1 PC2 0
峰7
峰8
峰3
峰5 －0.5 河北省承德市
峰6 辽宁省葫芦岛市
山东省临沂市
峰2 山西省运城市
峰4
S61 S67 S48 S66 S55 S72 S56 S58 S49 S65 S15 S59 S11 S17 S35 S36 S33 S34 S71 S46 S64 S18 S74 S16 S32 S41 S62 S39 S44 S4 S5 S2 S23 S13 S40 S27 S38 S69 S22 S60 S63 S68 S70 S9 S19 S20 S24 S30 S31 S73 S75 S52 S54 S50 S51 S53 S43 S42 S26 S14 S37 S45 S12 S29 S47 S1 S6 S8 S10 S25 S7 S21 S28 S57 S3 S76 S77 S78 －1 0 1
PC1
A. CA热图 B. PCA得分图
图3 78批山楂叶药材的CA热图和PCA得分图

· 2730 · China Pharmacy 2023 Vol. 34 No. 22 中国药房 2023年第34卷第22期