|
基于氢核磁共振与偏最小二乘法对酸枣仁及其掺伪品的鉴别基于氢核磁共振与偏最小二乘法对酸枣仁及其掺伪品的鉴别 申晨曦1，杜晨晖2,*，李震宇1，李爱平1，秦雪梅1，闫 艳1,* （1.山西大学中医药现代研究中心，山西 太原 030006；2.山西中医药大学中药学院，山西 太原 030619） 摘 要：采用氢核磁共振（1H nuclear magnetic resonance，1H NMR）代谢组学技术对酸枣仁、理枣仁以及不同比例的掺伪粉末进行对比，结合多元统计分析建立一种快速鉴别酸枣仁粉末及其掺伪品的方法。主成分分析显示掺伪比例较低的样品与酸枣仁有部分的重叠。偏最小二乘法-判别分析（partial least squares-discriminant analysis，PLS-DA）显示，当掺伪比例不低于40%时，酸枣仁粉和掺伪品沿t[1]轴完全分开。在此基础上，采用PLS对酸枣仁及其不同掺伪比例粉末样品训练集数据进行回归分析，并建立掺伪判别模型。1H NMR与PLS相结合的模型预测结果准确性较高，可实现酸枣仁中理枣仁掺假的快速鉴别。该方法为酸枣仁及含酸枣仁的保健食品、复方制剂质量评价提供了一条新途径。 关键词：酸枣仁；理枣仁；掺伪；氢核磁共振；代谢组学技术；偏最小二乘法 酸枣仁（Ziziphi spinosae Semen）为鼠李科植物酸枣（Ziziphus jujuba Mill. var. spinosa (Bunge) Hu ex H. F.Chou）的干燥成熟种子[1]，具有“养心补肝，宁心安神，敛汗，生津”的功效，是中医临床养心安神的首选药物[2]。酸枣仁是国家卫生和计划生育委员会最早颁布的“既是食品又是中药材物质”的中药材之一。由于其采收成本高，临床需求量大，价格不断攀升。市场上利用同科同属植物滇刺枣（Ziziphus mauritiana Lam.）的干燥成熟种子——理枣仁（Ziziphi mauritianae Semen）冒充酸枣仁和掺伪现象越来越多。作为酸枣仁最常见的伪品，理枣仁亦“具有宁心安神、除烦敛汗之功效”，作为云南地方习用品，早在1974年已被收入《云南省药品标准》[3]。 近年来，学术界对酸枣仁和理枣仁的鉴别方法也进行了大量研究，如性状鉴别[4]、显微鉴别[5]、理化鉴别[6-8]和DNA分子生物鉴别[9]等方法。以上方法为酸枣仁与理枣仁的真伪鉴别提供了科学基础，但存在一定局限性，如性状鉴别只适用于外观未破损的样品鉴别，显微鉴别需要具有丰富鉴别经验的专业人员进行真伪品的鉴别。然而理枣仁的外观与酸枣仁极相似，不易区别，特别是打粉后或入中成药后，难以区分和鉴定。为避免将廉价的理枣仁混入酸枣仁，以次充好，亟需建立一种快速、灵敏、准确的以理枣仁粉末形式掺伪酸枣仁粉末的鉴别方法。 植物代谢组学是代谢组学的一个重要分支。氢核磁共振（1H nuclear magnetic resonance，1H NMR）是植物代谢组学中常用的研究测试工具[10]，具有样品预处理简单，对样品无损害、检测时间短、重复性好等优势。1H NMR代谢组学技术可获得复杂混合体系中高通量的化学成分信息，包括化学位移和信号响应等，但是，同时各种化学成分信号的重叠也使得图谱变得十分复杂，仅靠肉眼观察只能获得有限的数据信息。多元统计分析能够从复杂的数据中最大限度的提取信息，反映样本分组的整体差异性[11]，具有很好的应用性。常用的多元统计分析方法包括无监督识别模式的主成分分析（principal component analysis，PCA）、有监督模式的偏最小二乘法-判别分析（partial least squares-discriminant analysis，PLS-DA）和PLS回归分析[12-13]。特别是应用多元统计分析方法建立的数学模型，能从未知样品的1H NMR图谱中预测样品的性质等信息，适合于中药材及饮片掺假的判别。 课题组前期采用1H NMR代谢组学技术对黄芪、款冬等中药材[14-17]进行了质量评价。另有研究采用1H NMR代谢组学技术对茶油、老陈醋、咖啡、蜂蜜和啤酒等食品进行了质量优劣评价[18-28]。本研究利用1H NMR代谢组学技术获得酸枣仁和理枣仁的化学信息，并结合PLS法对酸枣仁粉末掺伪比例进行鉴别研究。首先表征和分析酸枣仁、理枣仁和酸枣仁粉末掺假样品的1H NMR谱图。然后利用PLS对谱图的分段积分数据矩阵进行分析，建立酸枣仁粉末掺伪定量模型，进而对模型进行检验。最后，通过PLS模型预测未知样本的掺伪比例。该方法从整体的角度出发，以1H NMR结合多元统计分析应用于酸枣仁粉末掺伪的鉴别工作中，旨在为酸枣仁的定性分析提供一种新的思路，同时也为酸枣仁及以酸枣仁粉末入药的保健食品的质量评价提供一条新途径。 1 材料与方法1.1 材料与试剂酸枣仁和理枣仁分别采购于河北安国药材市场、昆明菊花村药材市场和山西振东道地药材有限公司，经山西中医药大学杜晨晖副教授鉴定酸枣仁为Z. jujuba var. spinosa的干燥成熟种子，理枣仁为Z. mauritiana Lam.的干燥成熟种子。共收集50 批样品，酸枣仁（25 批）和理枣仁（25 批），样品保存于山西大学中医药现代研究中心。 重水 美国Norell公司；氘代甲醇、三甲基硅烷丙酸钠盐 青岛腾龙微波科技有限公司；氘代氢氧化钠瑞士Armar公司。 1.2 仪器与设备600 MHz Avance III 1H NMR 瑞士Bruker公司；Neofuge 13R高速冷冻离心机 上海力申科学仪器有限公司；RE-52A旋转蒸发仪 上海亚荣生化仪器厂；pH计 奥豪斯仪器有限公司。 1.3 方法1.3.1 样品制备 在共直流母线分布式发电系统中,与PV面板直接相连的DC/DC变换器,也被称为光伏直流模块PDM(Photovoltaic DC Module),是系统的关键装置之 25 批酸枣仁样品（酸枣仁1～25）各称取50 g，分别粉碎，过3号筛（50 目），即得25 批酸枣仁粉末样品。25 批理枣仁样品（理枣仁1～25）各称取50 g，分别粉碎，过3号筛，即得25 批理枣仁粉末样品。从25 批酸枣仁粉末样品中随机选取6 批，并从25 批理枣仁粉末样品中随机选取6 批。将所选酸枣仁粉末（n=6）中分别掺入不同质量比（5%、10%、20%、40%、60%、80%）的理枣仁粉末（n=6），共36 个掺伪样品，样品信息见表1。 表 1 实验样品信息
Table 1 Information about pure and mixed samples of ZSS and ZMS tested in this study  样品编号 样品比例 样品编号 样品比例酸枣仁1～25 100%酸枣仁 理枣仁1～25 100%理枣仁掺伪1～6 5%理枣仁 掺伪19～24 40%理枣仁掺伪7～12 10%理枣仁 掺伪25～30 60%理枣仁掺伪13～18 20%理枣仁 掺伪31～36 80%理枣仁
精密称取以上粉末样品各200 mg，置于10 mL离心管中，分别加蒸馏水及甲醇各1.5 mL，氯仿3 mL；旋涡混匀1 min；超声提取25 min；3 500 r/min离心25 min，上层为甲醇水相（水溶性部分），下层为氯仿相；将甲醇水层移至25 mL圆底烧瓶中，减压浓缩至干，氯仿层弃去。测定前用氘代甲醇400 μL与含有缓冲盐的重水（KH2PO4溶于氘代水中，以1 mol/L氘代氢氧化钠溶液调节pH值至6，含0.01%磷酸三钠（trisodium phosphate，TSP））400 μL进行溶解，移至1.5 mL离心管中，离心10 min（13 000 r/min），移取上清液600 μL于5 mm核磁管中，待测。 1.3.2 样品1H NMR测试条件 (2)FRP布加固不仅能有效减小梁的初始弹性挠度和蠕变稳态挠度,而且可以显著降低裂纹对DF梁弯曲变形的影响. 样品在25 ℃、600 MHz 1H NMR上测定，测定频率600.13 MHz，扫描范围为δ -5～15，扫描次数64 scans，采用Noesygppr1d序列压制水峰，用氘代甲醇进行锁场，内标为TSP。 1.3.3 数据处理 将测得的各样品自由衰减信号导入MestReNova（version 8.0.1, Mestrelab Research, Santiago de Compostlla,Spain）软件并将图谱以TSP（0.00）定标，进行相位调整和基线校正。然后以0.01分段对化学位移区间δ 0.72～9.15进行分段积分，其中对水峰区域δ 4.68～4.92和残余甲醇峰δ 3.34～3.38不进行积分，积分数据采用峰面积归一化处理。 在“大数据”的背景下，企业应该学会将新兴的互联网技术同人工神经网络模型的建立相结合。依据大数据的ANN模型，把获得的数据同企业内部审计预警机制相结合，并通过实际运用的结果的分析，测试数据的准确性和完整性及企业风险内部的审计过程。 1.3.4 精密度实验和重复性实验 取编号为酸枣仁1样品，按照1.3.1节方法制备样品，并按照1.3.2节方法连续测定6 次，将6 次所测图谱分别按照1.3.3节方法进行数据处理，得到6 组积分面积数据，依据夹角余弦法，采用夹角余弦公式，计算6 张谱图之间的夹角余弦值，考察仪器的精密度。取编号为酸枣仁1样品粉末，按照1.3.1节方法平行制备6 份样品，并按照1.3.2节下方法测定，将6 个所测图谱分别按照1.3.3节方法进行数据处理，得到6 组积分面积数据，计算6 份样品之间的夹角余弦值，考察制备方法的精密度[29]。夹角余弦值（cosθ）计算公式如下： 三是耳螨离开宿主后，在适宜环境温度下可以存活4～15天，在治疗螨虫病的同时需要对笼具、食盒、排泄物等进行消毒； 
式中：Xi为样品X的第i个峰的相对积分值；Yi为样品Y的第i个峰的相对积分值。1.3.5 多元统计分析 1.3.5.1 PCA和PLS-DA 由图5可以看出，当菌粉与奶果混合液1∶1000（质量浓度比）混合发酵时，酸奶凝乳最好，感官评分为86分；而在菌种添加量递减过程中，酸奶品质随之下降；当添加量为菌粉∶混合液=1∶2500时，酸奶凝乳不均匀，香味不明显，口感稀薄，综合品质最差。因而可得，菌粉最佳添加量为1∶1000。 将归一化后的数据导入SIMCA-P 13.0（Umetrics，瑞典）软件进行PCA和PLS-DA。PCA不能忽略组内误差和消除无关的随机误差，因此需要进一步采用有监督的PLS-DA，减小组内差异，扩大组间差异，以实现酸枣仁及其伪品的准确鉴别。 1.3.5.2 PLS预测模型的构建 在收集的样品中随机选取其中20 批酸枣仁（酸枣仁1、酸枣仁3～21）、21 批理枣仁（理枣仁1～21）及24 批掺伪样品（掺伪1～4、掺伪7～10、掺伪13～16、掺伪19～22、掺伪25～28、掺伪31～34）构成训练集，共65 组强度积分值数据；余下的4 批酸枣仁（酸枣仁22～25）、4 批理枣仁（理枣仁22～25）和12 批酸枣仁掺伪样品（掺伪5～6、掺伪11～12、掺伪17～18、掺伪23～24、掺伪29～30、掺伪35～36）组成测试集，共20 组强度积分值数据。 基于PLS原理[20]，按照样本类别特征赋予训练集样本分类变量值，作为因变量（Y变量）。首先将酸枣仁赋值为1，掺伪5%理枣仁赋值1.05，掺伪10%理枣仁赋值1.1，掺伪20%理枣仁赋值1.2，掺伪40%理枣仁赋值1.4，掺伪60%理枣仁赋值1.6，掺伪80%理枣仁赋值1.8，理枣仁赋值为2。然后将训练集数据导入SIMCA-P 13.0软件中，采用pareto（帕雷托）换算对训练集65 个积分强度值数据矩阵进行标准化，归一化的强度积分值作为自变量（X变量）。采用PLS法将标准化的数据矩阵X变量与类别变量（Y变量）进行线性回归，建立预测模型。 这一创新做法，从2013年起连续6年被写入安徽省委一号文件，2017年被原农业部等六部委联合发文向全国推广，2018年被写入中央一号文件。为抓好政策扶持和配套服务，宿州市委、市政府出台了《关于扶持农业产业化联合体发展的若干政策意见》等文件，市、县区两级财政每年兑现产业化联合体成员各类奖补资金3600多万元。 将测试集样品的数据代入经过检验的模型，计算得到各样品的类别变量预测值P，并将预测值与训练集得出的鉴别范围比对，对测试集样品进行鉴别。 2 结果与分析2.1 仪器精密度与方法精密度实验图1A为同一样品连续测定6 次所得图谱，按照1.3.3节方法进行数据处理，得到6 组积分面积数据，计算6 张谱图之间的夹角余弦值分别为1.000、0.997、0.999、0.998、0.998、0.997，计算值均在0.99以上，表明仪器精密度良好。图1B为相同条件下同一批次酸枣仁平行制备6 份实验所得图谱，按照1.3.3节方法进行图谱处理，得到6 组积分面积数据，计算6 份样品之间的夹角余弦值分别为1.000、0.998、0.998、0.993、0.998、0.999，计算值均在0.99以上，表明制备方法精密度良好。这为定性鉴别酸枣仁和理枣仁提供了稳定可靠的分析条件。  图 1 仪器精密度实验和重复性实验1H NMR谱图
Fig. 1 1H NMR spectra showing instrumental precision and repeatability
2.2 1H NMR图谱分析 图 2 酸枣仁、掺伪样品及理枣仁1H NMR（δ 0～3）谱图
Fig. 2 1H NMR (δ 0–3) spectra of ZSS, ZMS and their mixtures
由图1可知，酸枣仁样品的1H NMR谱图，大致可以分为脂肪区（δ 0.00～3.00）、糖类化合物区（δ 3.00～6.00）和芳香区（δ 6.00～9.00）。通过比较酸枣仁、理枣仁及酸枣仁粉掺伪样品的1H NMR谱图，脂肪区出现变化较大的信号。通过对化学位移、耦合常数J、HMDB数据库和文献数据以及1H-1H化学位移相关谱的分析，鉴定出5 个变化较为明显的物质[30]。分别为化合物1缬氨酸（δ=1.01，d，J=7.2 Hz；δ=1.07，d，J=7.2 Hz）、化合物2乳酸（δ=1.33，d，J=6.6 Hz）、化合物3丙氨酸（δ=1.49，d，J=7.2 Hz）、化合物4 γ-氨基丁酸（δ=2.31，t，J=7.2 Hz）、化合物5半胱氨酸（δ=3.05，s）。随着酸枣仁粉掺伪比例的增加，观察到5 个特征峰的信号强度逐渐减弱（图2）。表明酸枣仁粉末较伪品含有更高的氨基酸类成分。 2.3 酸枣仁粉及掺伪品的PCA在进行PCA之前，通过Hotelling’s T2 Range算法判断异常样本。当T2超过临界值时，说明该样本与其他样本间存在较大的差异。如图3所示，样品中存在一个异常值（酸枣仁2），应作为异常点剔除，剩余85 个样品用于下一步的模型建立。  图 3 酸枣仁粉及掺伪样品的Hotelling’s T2 Range图
Fig. 3 Hotelling’s T2 Range plots of pure and adulterated samples
 图 4 剔除异常点后酸枣仁粉及掺伪样品的PCA得分图
Fig. 4 PCA score plot for pure and adulterated samples after elimination of outliers
从图4可知，酸枣仁粉和掺伪样品的分类情况。前5 个主成分共解释了72.4%的方差，酸枣仁和掺伪样品在t[2]轴上有较好区分。大部分酸枣仁样品在t[2]轴的下部，大部分掺伪样品在t[2]轴的上部，由于化学成分的相似性，掺伪比例较低的样品与酸枣仁有部分的重叠，随着掺伪比例的增加，样品的分布在t[1]轴上呈从右到左的分布。 2.4 酸枣仁粉及掺伪品PLS-DA无监督模式识别的PCA无法判断未知样品，且存在与研究目的无关的组内误差以及随机误差，不利于分组信息的准确性。为确定分组样品的差异，进一步采用有监督模式识别的PLS-DA。其分析过程中可以消除众多化学信息中相互重叠的部分，使得分析数据更加准确可靠[26]。 面前一杯咖啡已经冰凉了，车子还没来。上次接了她去，又还在公寓里等了快一个钟头他才到。说中国人不守时刻，到了官场才登峰造极了。再照这样等下去，去买东西店都要打烊了。 为准确识别掺伪样品，首先将25 个理枣仁样品和36 个掺伪样品归为一类，酸枣仁作为另一类，建立PLSDA模型。如图5A～C所示，当掺伪比例分别为5%、10%和20%时，酸枣仁粉和掺伪品有部分交叉。当掺伪比例不低于40%时（图5D），酸枣仁粉和掺伪品在t[1]轴上可以完全分开。随着掺伪比例的增加，酸枣仁样品所在区域和掺伪品所在区域间的距离逐渐增加（图5E～G）。为进一步验证该模型的可靠性，采用排列实验法对图5D中对应的PLS-DA模型进行验证。通过200 次随机改变类别变量Y的排列顺序，得到累计贡献率R2和预测能力Q2。其中排列实验模型的R2回归线和Q2回归线与纵轴的截距分别为0.474和-0.483，最右端鉴别模型的原始Q2值大于左边任何一个Y变量随机排列模型的Q2值（图6），说明所构建的PLS-DA模型没有出现过拟合现象，且有较好预测能力，显示模型可靠。 
 图 5 酸枣仁和不同比例掺伪品的PLS-DA图
Fig. 5 Score plots of PLS-DA for pure ZSS and adulterated samples
 图 6 排列实验对PLS-DA模型的可靠性验证
Fig. 6 Validation of PLS-DA model by permutation test
2.5 基于PLS的酸枣仁粉及掺伪品的定量分析为更准确地检测实际掺伪比例，进一步采用PLS法，建立掺假量预测回归模型。训练集均方根误差（root mean square error of estimation，RMSEE）、预测集均方根误差（root mean square error of prediction，RMSEP）和训练集相关系数（R2）可以评价PLS模型，R2越大，RMSEE和RMSEP越小，说明模型的准确性和精确度越高。本研究训练集由65 个样品构成（酸枣仁样品20 个、理枣仁样品21 个，掺伪样品24 个），共获得802 个强度积分变量，构成65×802的数据矩阵。利用PLS对其进行分析，建立模型。以训练集的真实值为横坐标，预测值为纵坐标，绘制的散点图（图7）。训练集模型中R2为0.986 5，大于0.98，RMSEE为0.077 0。利用PLS预测模型对训练集的类别变量进行计算。由表2可知，酸枣仁粉的计算值范围为0.86≤P≤1.15；掺伪5%理枣仁的预测值范围为1.00≤P≤1.18；掺伪10%理枣仁的预测值范围为1.06≤P≤1.18；掺伪20%理枣仁的预测值范围为1.15≤P≤1.33；掺伪40%理枣仁的预测值为1.39≤P≤1.47；掺伪60%理枣仁的预测值为1.51≤P≤1.69；掺伪80%理枣仁的预测值范围为1.68≤P≤1.75，理枣仁粉的计算值范围为1.93≤P≤2.07。  图 7 训练集和验证集真实值和预测值关系图
Fig. 7 Plot of true values versus estimated values for training set and validation set
表 2 训练集的真实值和预测值
Table 2 True values and estimated values for training set samples  样品编号 真实值 预测值 样品编号 真实值 预测值 样品编号 真实值 预测值酸枣仁1 1 0.88 理枣仁3 2 1.97 掺伪4 1.05 1.09酸枣仁3 1 1.07 理枣仁4 2 2.00 掺伪7 1.1 1.12酸枣仁4 1 0.96 理枣仁5 2 1.97 掺伪8 1.1 1.18酸枣仁5 1 0.97 理枣仁6 2 1.94 掺伪9 1.1 1.06酸枣仁6 1 0.98 理枣仁7 2 1.97 掺伪10 1.1 1.11酸枣仁7 1 0.86 理枣仁8 2 1.96 掺伪13 1.2 1.23酸枣仁8 1 0.92 理枣仁9 2 1.93 掺伪14 1.2 1.33酸枣仁9 1 1.08 理枣仁10 2 2.02 掺伪15 1.2 1.15酸枣仁10 1 0.89 理枣仁11 2 2.00 掺伪16 1.2 1.22酸枣仁11 1 1.01 理枣仁12 2 2.07 掺伪19 1.4 1.39酸枣仁12 1 1.15 理枣仁13 2 2.06 掺伪20 1.4 1.47酸枣仁13 1 1.01 理枣仁14 2 1.99 掺伪21 1.4 1.40酸枣仁14 1 1.04 理枣仁15 2 2.05 掺伪22 1.4 1.40酸枣仁15 1 1.00 理枣仁16 2 2.03 掺伪25 1.6 1.58酸枣仁16 1 1.15 理枣仁17 2 2.02 掺伪26 1.6 1.69酸枣仁17 1 1.12 理枣仁18 2 2.02 掺伪27 1.6 1.51酸枣仁18 1 1.09 理枣仁19 2 2.01 掺伪28 1.6 1.58酸枣仁19 1 0.98 理枣仁20 2 2.03 掺伪31 1.8 1.73酸枣仁20 1 0.96 理枣仁21 2 2.02 掺伪32 1.8 1.75酸枣仁21 1 0.98 掺伪1 1.05 1.08 掺伪33 1.8 1.68理枣仁1 2 1.99 掺伪2 1.05 1.18 掺伪34 1.8 1.72理枣仁2 2 1.96 掺伪3 1.05 1.00
表 3 测试集的真实值和预测值
Table 3 True values and estimated values for validation set samples  样品编号 预测值 鉴别结果 样品编号 预测值 鉴别结果酸枣仁22 1.05 正确 掺伪11 1.07 错误酸枣仁23 1.14 正确 掺伪12 1.16 错误酸枣仁24 1.12 正确 掺伪17 1.20 正确酸枣仁25 1.12 正确 掺伪18 1.25 正确理枣仁22 2.05 正确 掺伪23 1.43 正确理枣仁23 2.03 正确 掺伪24 1.45 正确理枣仁24 2.01 正确 掺伪29 1.62 正确理枣仁25 2.04 正确 掺伪30 1.57 正确掺伪5 1.05 错误 掺伪35 1.74 正确掺伪6 1.07 错误 掺伪36 1.78 错误
为进一步检验PLS模型，将酸枣仁粉及不同比例掺伪品的测试集样品代入训练集创建的PLS模型中（图7），预测结果见表3。测试集RMSEP为0.058 200。RMSEE和RMSEP值均接近于0，说明模型可靠，可准确预测掺假理枣仁的量。测试集的预测结果P值表明，酸枣仁粉掺伪5%和10%时，与酸枣仁粉预测值有交叉，无法准确预测其含量。当掺伪比例达到20%时，样品计算所得P值除掺伪36号样品稍有偏差其余样本的P值均在在所建立模型的P值范围之内。说明该模型具有较好的预测和评价未知样品掺伪量的能力。掺伪36号样品为掺伪80%，与模型中的P值较为接近。 算法：设g:H→H是单值映射，Ni:H→H(i=1,2)是单值映射；A，B:H→C(H)是集值映射，设泛函 φ:H→是真凸下半连续，次可微的，且 g(u)∈∂φ其中∂φ表示φ的次微分。 3 结 论本研究利用1H NMR代谢组学技术获得了酸枣仁粉及不同比例掺伪品的化学信息。直观分析显示缬氨酸、乳酸、丙氨酸、γ-氨基丁酸、半胱氨酸5 个氨基酸成分随着掺伪比例的增加表现出信号强度逐渐减弱的趋势。PCA表明掺伪比例较低的样品与酸枣仁有部分的重叠。PLS-DA显示，当掺伪比例不低于40%时，酸枣仁粉和掺伪品沿t[1]轴完全分开。在此基础上，通过PLS模型可准确预测掺假酸枣仁粉中理枣仁的掺入量。1H NMR结合PLS-DA和PLS可以预测酸枣仁粉掺假量，具有分析时间短、可控性强、准确性高等优点，可实现未知样本的准确判别，对于酸枣仁和理枣仁的鉴别及掺伪区分有明显优势。该方法特别适合于性状鉴别和薄层色谱等方法难以实现的酸枣仁药材真伪鉴别。在后续研究中，应增加样本量，对本鉴别模型数据扩充，使鉴别结果具有代表性，为酸枣仁及含酸枣仁的复方制剂质量评价提供一种新方法和新途径。 参考文献： [1] 国家药典委员会. 中国药典: 一部[M]. 2015版. 北京: 中国医药科技出版社, 2015: 366. [2] SHERGIS J L, NI X, SARRIS J, et al. Ziziphus spinosa seeds for insomnia: a review of chemistry and psychopharmacology[J].Phytomedicine, 2017, 34: 38-43. DOI:10.1016/j.phymed.2017.07.004. [3] 云南省卫生局. 云南省药品标准[M]. 昆明: 云南省卫生厅, 1974. [4] 黄德耀. 从外观性状识别酸枣仁伪品[J]. 中国实用医药, 2015, 10(1):221-222. DOI:10.14163/j.cnki.11-5547/r.2015.01.165. [5] 王雪, 杨光, 周修腾, 等. 酸枣仁、理枣仁及枳椇子的显微鉴别[J].现代中药研究与实践, 2018, 32(3): 9-12. DOI:10.13728/j.1673-6427.2018.03.003. [6] SUN S, LIU H L, XU S J, et al. Quality analysis of commercial samples of Ziziphi spinosae Semen (suanzaoren) by means of chromatographic fingerprinting assisted by principal component analysis[J]. Journal of Pharmaceutical Analysis, 2014, 4(3): 217-222.DOI:10.1016/j.jpha.2014.01.003. [7] GUO S, DUAN J N, LI Y Q, et al. Comparison of the bioactive components in two seeds of Ziziphus species by different analytical approaches combined with chemometrics[J]. Frontiers in Pharmacology, 2017, 8: 609. DOI:10.3389/fphar.2017.00609. [8] ZHANG F X, LI M, QIAO L R, et al. Rapid characterization of Ziziphi spinosae Semen by UPLC/Qtof MS with novel informatics platform and its application in evaluation of two seeds from Ziziphus species[J].Journal of Pharmaceutical & Biomedical Analysis, 2016, 122: 59-80.DOI:10.1016/j.jpba.2016.01.047. [9] 刘萍, 古今, 周文斌, 等. 酸枣仁与其混淆品及伪品的RAPD分析[J].解放军药学学报, 2006(1): 67-70. [10] KORTESNIEMI M, SINKKONEN J, YANG B, et al. NMR metabolomics demonstrates phenotypic plasticity of sea buckthorn(Hippophae rhamnoides) berries with respect to growth conditions in Finland and Canada[J]. Food Chemistry, 2017, 219: 139-147.DOI:10.1016/j.foodchem.2016.09.125. [11] 李震宇, 段亚辉, 秦雪梅, 等. 中药质量差异性研究的思考[J]. 药学学报, 2017, 52(12): 1820-1826. DOI:10.16438/j.0513-4870.2017-0622. [12] 阿基业, 何骏, 孙润彬. 代谢组学数据处理: 主成分分析十个要点问题[J]. 药学学报, 2018, 53(6): 929-937. DOI:10.16438/j.0513-4870.2017-1288. [13] BYLESJÖ M, RANTALAINEN M, CLOAREC O, et al. OPLS discriminant analysis: combining the strengths of PLS-DA and SIMCA classification[J]. Journal of Chemometrics, 2010, 20: 341-351.DOI:10.1002/cem.1006. [14] LI A P, LI Z Y, SUN H F, et al. Comparison of two different Astragali Radix by a ¹H NMR-based metabolomic approach[J]. Journal of Proteome Research, 2015, 14(5): 2005-2016. DOI:10.1021/pr501167u. [15] LI Z Y, ZHANG Z Z, DU G H, et al. Comparative analysis of Danggui and European Danggui using nuclear magnetic resonance-based metabolic fingerprinting[J]. Journal of Pharmaceutical & Biomedical Analysis, 2015, 103: 44-51. DOI:10.1016/j.jpba.2014.10.028. [16] LI Z Y, SUN H M, XING J, et al. Chemical and biological comparison of raw and vinegar-baked Radix Bupleuri[J]. Journal of Ethnopharmacology, 2015, 165: 20-28. DOI:10.1016/j.jep.2015.02.024. [17] LI Z Y, LI J, ZHANG Z Z, et al. NMR-based metabolomic analyses for the componential differences and the corresponding metabolic responses of three batches of Farfarae Flos[J]. Chemometrics &Intelligent Laboratory Systems, 2017, 165: 1-10. DOI:10.1016/j.chemolab.2017.03.010. [18] 王晓玲, 吴晶, 谭明乾. 低场核磁共振结合化学计量学方法快速检测掺假核桃油[J]. 分析测试学报, 2015, 34(7): 789-794.DOI:10.3969/j.issn.1004-4957.2015.07.006. [19] ZHU M T, SHI T, CHEN Y, et al. Prediction of fatty acid composition in camellia oil by 1H NMR combined with PLS regression[J]. Food Chemistry, 2019, 279: 339-346. DOI:10.1016/j.foodchem.2018.12.025. [20] 范玛莉, 李震宇, 李爱平, 等. 偏最小二乘方法在不同食醋判别分析中的应用[J]. 食品工业科技, 2014, 35(21): 324-328. DOI:10.13386/j.issn1002-0306.2014.21.061. [21] CONSONNI R, POLLA D, CAGLIANI L R. Organic and conventional coffee differentiation by NMR spectroscopy[J]. Food Control, 2018,94: 284-288. DOI:10.1016/j.foodcont.2018.07.013. [22] CONSONNI R, POLLA D, CAGLIANI L R. NMR-based metabolomic approach to differentiate organic and conventional Italian honey[J]. Food Control, 2019, 98: 133-140. DOI:10.1016/j.foodcont.2018.11.007. [23] 陈雷, 刘红兵, 罗立廷, 等. 氢核磁共振结合正交偏最小二乘法对油菜蜜中果葡糖浆掺假的判别分析[J]. 食品科学, 2017, 38(4): 275-282. DOI:10.7506/spkx1002-6630-201704045. [24] AUGUSTO D S L, LUIZ F D, GOMES T A, et al. Discrimination of Brazilian lager beer by 1H NMR spectroscopy combined with chemometrics[J]. Food Chemistry, 2019, 272: 488-493. DOI:10.1016/j.foodchem.2018.08.077. [25] 佘僧, 李熠, 宋洪波, 等. 低聚糖和多酚结合偏最小二乘判别分析鉴别油菜蜜产地[J]. 食品科学, 2019, 40(12): 290-295. DOI:10.7506/spkx1002-6630-20180424-308. [26] 石婷, 陈倩, 闫小丽, 等. 核磁共振氢谱结合化学计量学快速检测掺假茶油[J]. 食品科学, 2018, 39(22): 241-248. DOI:10.7506/spkx1002-6630-201822037. [27] SHI T, ZHU M T, ZHOU X Y, et al. 1H NMR combined with PLS for the rapid determination of squalene and sterols in vegetable oils[J]. Food Chemistry, 2019, 287: 46-54. DOI:10.1016/j.foodchem.2019.02.072. [28] 何瑶, 郑彦婕, 邓伶莉, 等. 五常稻花香米的1H-NMR波谱分析及掺假鉴别[J]. 食品工业科技, 2016, 37(12): 80-84; 171. DOI:10.13386/j.issn1002-0306.2016.12.007. [29] 王亚男, 孙伟, 贺文义, 等. 基于1H-NMR指纹图谱的四批市售桑白皮药材质量分析靳香菊[J]. 沈阳药科大学学报, 2017, 34(1): 62-70.DOI:10.14066/ j.cnki.cn21-1349/r.2017.01.012. [30] 李静, 魏玉海, 秦雪梅, 等. 基于NMR代谢组学技术的当归不同部位化学成分比较[J]. 中草药, 2017, 48(7): 1409-1415. DOI:10.7501/j.issn.0253-2670.2017.07.025.
Differentiation between Authentic and Adulterated Ziziphi spinosae Semen by 1H NMR Spectroscopy Combined with Partial Least Squares SHEN Chenxi1, DU Chenhui2,*, LI Zhenyu1, LI Aiping1, QIN Xuemei1, YAN Yan1,*
(1. Modern Research Center for Traditional Chinese Medicine, Shanxi University, Taiyuan 030006, China;2. School of Traditional Chinese Materia Medica, Shanxi University of Traditional Chinese Medicine, Taiyuan 030619, China) Abstract: In the present study, the chemical composition of pure powders of Ziziphi spinosae Semen (ZSS) and Ziziphi mauritianae Semen (ZMS) and their mixtures in different proportions was compared by 1H nuclear magnetic resonance(1H NMR)-based metabolomic fingerprinting. Then the obtained spectral data were analyzed by multivariate analysis to develop a new method for rapid detection of ZSS adulteration. Principal component analysis (PCA) performed on the 1H NMR data revealed that only a few samples at low adulteration ratios were overlapped with ZSS. In addition, the data were also analyzed by partial least squares-discriminant analysis (PLS-DA). A clear separation was seen between ZSS and 40% adulterated samples. Then a PLS model was established based on the 1H NMR data and was applied to detect adulteration in ZSS falsified by intentional addition of different concentrations of ZMS, and a quantitative calibration model was established. The PLS regression approach provided a highly reliable and accurate prediction model. The results demonstrated that 1H NMR spectroscopy coupled with PLS is a fast and effective technique for checking adulteration in ZSS. Meanwhile, the method provides a potential tool for quality control of ZSS. Keywords: Ziziphi spinosae Semen; Ziziphi mauritianae Semen; adulterants; 1H nuclear magnetic resonance; metabolomics;partial least squares
收稿日期：2019-03-14 基金项目：国家自然科学基金青年科学基金项目（81603289；81603251）；山西省研究生联合培养基地人才培养项目（2017JD01）；山西省科技攻关计划-振东专项（2016ZD0105）；山西省科学技术厅-地产中药功效物质研究与利用山西省重点实验室项目（201605D111004）；山西大学引进人才事业发展经费项目（226545006）；山西省研究生教育创新项目 第一作者简介：申晨曦（1993—）（ORCID: 0000-0003-2629-7343），女，硕士研究生，研究方向为中药质量评价。E-mail: 939948694@qq.com*通信作者简介： 目前，我国全面提倡实施素质教育，同时，提倡减轻学生的学业负担。“变题”的方法和技术成为教学中常用的教学方法。通过“变题”练习把规律性的问题联系在一起，不仅对此类题型做了归纳总结，对解题方法有了进一步的认识，避免了学生的“题海”战术，减轻学生的课业负担，提高课堂的教学效率，对知识的掌握，思维和能力培养都起着非常重要的作用。 DOI:10.7506/spkx1002-6630-20190314-180 中图分类号：TS207.7 文献标志码：A 文章编号：1002-6630（2020）08-0275-07 引文格式： 申晨曦, 杜晨晖, 李震宇, 等. 基于氢核磁共振与偏最小二乘法对酸枣仁粉末及其掺伪品的鉴别[J]. 食品科学, 2020,41(8): 275-281. DOI:10.7506/spkx1002-6630-20190314-180. http://www.spkx.net.cn例如，教师在给学生讲解《鲸》的课文内容时，教师可以为学生制作多媒体动画，让学生欣赏鲸鱼的进化过程，了解更多的科学知识。教师再引导学生学习鲸鱼不是鱼的原因，启发学生进行学习和思考。然后，教师可以让学生简单概括鲸鱼的特点，然后再为学生拟定标题，让学生练习写科普类的作文，从而发展学生的科学修养，锻炼学生的写作能力。 SHEN Chenxi, DU Chenhui, LI Zhenyu, et al. Differentiation between authentic and adulterated Ziziphi spinosae Semen by 1H NMR spectroscopy combined with partial least squares[J]. Food Science, 2020, 41(8): 275-281. (in Chinese with English abstract) DOI:10.7506/spkx1002-6630-20190314-180. http://www.spkx.net.cn
| 
|Archiver|手机版|小黑屋|www.openhelp100.com
( 冀ICP备19026749号-1 )
发表于 2021-2-3 22:06:24
