胡麻(Linum usitatissimum L.),又名油用亚麻,为亚麻科亚麻属草本植物,在我国主要种植于内蒙古、甘肃、宁夏、山西和新疆等地,年产量33.5万t [1-2]。胡麻油α-亚麻酸含量在45%～70%,在人体内可转换为DHA和EPA,具有调节血脂、降低低密度脂蛋白、软化血管和降血压等功效[3-4]。目前我国胡麻油市场价格是其他食用植物油的3～20倍,极易受到廉价油脂、餐厨垃圾油和地沟油的掺假或冒充。不同加工精炼方式生产的胡麻油产品级别档次不同[5]。热榨胡麻籽油烘烤和压榨温度可高达180 ℃以上,破坏不饱和脂肪酸,生成反式脂肪酸(TFA)、醛类、多环芳烃等有害物质,显著增高苯并芘[6-8]。冷榨工艺是在60 ℃以下压榨,可有效保留营养并控制有害物质生成[9]。建立判别胡麻油物种、加工方法和品级的模型具有重要意义。

指纹特征是指研究对象的一系列化学成分的种类和比例,以及仪器指纹谱数据集和矩阵中潜藏的研究对象特有的模式或特征[10]。每种植物油应该具有特定的脂肪酸(FAs)构成或模式,即其FAs指纹。化学计量学数据结构分析及可视化、模式和判别分析,为解析和挖掘数据集、数据矩阵等大数据结构和隐藏的规律模式,及利用食物大数据描绘和分类食物,提供了无限的可能[10-11]。目前基于FAs指纹进行胡麻油掺假判别的研究较少。Wang等[12]在胡麻油中掺不同比例菜籽油、玉米油、花生油、葵花籽油、大豆油和芝麻油,测定了C16:0、C18:0、C18:1ω9 、C18:2ω6和C18:3ω3这5种脂肪酸,并进行了相似性评价、判别分析和主成分分析(principal component analysis,PCA),结果表明掺假体积分数大于10%时判别准确率可达100%。但PCA对胡麻油和其他油脂的聚类分离效果并不直观。Sun等[13]用气相色谱-质谱检测胡麻油和其他5种油的FAs,用PCA和递归支持向量机进行胡麻油和非胡麻油二分类判别,鉴别掺伪体积分数为10%的胡麻油样品准确率可达95.6%。这些报道中对照油种类少,且尚未对胡麻油冷榨和热榨工艺及品级开展判别尝试,此外化学计量学解析方法方面也需要进行新的探索。

本研究拟观察胡麻油和其他包括餐厨垃圾油和地沟油样品在内的10种油脂的12种FAs聚类特征,评价以FAs指纹鉴别胡麻油物种、冷热榨工艺及品级真实性的可行性,建立胡麻油掺菜籽油样品的PLS定量预测模型。研究旨在为打击胡麻油掺假冒充,培育正常市场竞争秩序,维护消费者权益,优化榨油和精炼工艺,促进品牌化发展提供一定帮助,为胡麻油科研、生产和产品标准制定提供参考。

1 材料与方法

1.1 材料与试剂

1.1.1 油脂样品信息

在前期植物油近红外光谱和FAs指纹特征预研究的基础上,从上百份胡麻油样品中保留30份代表性品牌和小作坊胡麻油样品,进行以FAs指纹判别胡麻油真实性的可行性研究。实验所用油脂样品共116份,详细信息见表1。30份胡麻油样品分别为内蒙古呼和浩特、集宁、锡林郭勒、巴彦淖尔胡麻油加工企业采集的5份冷榨胡麻油(聚类参照样品)及家庭作坊采集的3份热榨胡麻油(聚类参照样品),其余22份市售和联系厂家邮寄的品牌胡麻油样品均来自内蒙古自治区和山西省。采集市售品牌菜籽油、葵花籽油、大豆油、花生油、玉米油、芝麻油和调和油各10份(共70份)。16份废弃油分别为12份餐厨垃圾油(食堂和饭馆废弃剩余饭菜垃圾中提取未精炼的油脂,采集自呼和浩特市某环保公司)、4份精炼地沟油(国家食品安全风险评估中心)。采样后分装于棕色样品瓶,充入氮气拧紧瓶盖用封口膜密封,4 ℃保存备用。

表1 样品信息

Tab.1 Sample information 份

样品加工方法品级产地冷榨热榨浸出水代一级二级三级四级内蒙古山西其他品牌总计胡麻油1812--2010--2010-2530菜籽油46--3-43--10810葵花籽油10---10---3-71010大豆油3-7-6-4---101010花生油46--91----101010玉米油10---10---1-91010芝麻油6--491----101010调和油----------10910废弃油餐厨垃圾油------------12地沟油------------4

-表示未采集样品。

人工制备掺假样本:选择已知背景冷榨胡麻油样本5份和5份菜籽油随机组合,将菜籽油按质量分数0、25%、50%、75%和100%掺入胡麻油中制备掺假样本,每个比例10份,共计50份。

1.1.2 实验试剂

甲醇,色谱纯,北京迈瑞达科技有限公司;正庚烷、三氟化硼-甲醇,色谱纯,上海麦克林生化科技有限公司;氢氧化钠、氯化钠、无水硫酸钠,分析纯,天津市风船化学试剂科技有限公司;十一碳酸甘油三酯标准品(纯度99%)、37 组分脂肪酸甲酯混合标准品,美国Sigma公司。

1.2 仪器与设备

FL-9720型气相色谱仪、FL-9720型氢火焰离子检测器,浙江福立分析仪器有限公司;KQ2200E型超声波清洗器,昆山市超声仪器有限公司;JG-HH-S26型数显恒温水浴锅,济南精城实验仪器有限公司;AL204型电子天平,梅特勒-托利多仪器有限公司。

1.3 实验方法

脂肪酸的测定参照GB 5009.168—2016《食品安全国家标准 食品中脂肪酸的测定》[14]进行前处理,116份样品平行测定2次,获得232套数据。

气相色谱条件:毛细管色谱柱SPTM-2560(100 m×0.25 mm×0.20 μm),美国Supplco公司。初始温度100 ℃,保持5 min,以3 ℃/min升温至240 ℃,保持35 min;分流式进样,分流比为20∶1,载气为氮气,氮气流速30 mL/min,隔垫吹扫2 mL/min,氢气流速30 mL/min,空气流速300 mL/min;注射样品量为1 μL。

1.4 数据处理

1.4.1 聚类和分类分析

116份样品平行测定获得的232套数据,用SIMCA 14.1化学计量学软件对各研究分类(油脂种类、胡麻油冷榨和热榨、胡麻油一级和二级)进行聚类分析和判别分析。对FAs数据进行center scaling(中心缩放)、unit variance scaling(单位方差缩放)和Pareto scaling(帕累托缩放)等预处理后,利用PCA观察油脂样品以FAs指纹自然聚类的特征。数据最优预处理结合正交偏最小二乘判别分析(orthogonal partial least squares discriminant analysis,OPLS-DA),进行有监督的聚类分析和研究分类的判别分析;通过根向量图找出对分类贡献大的关键FAs。

1.4.2 描述性统计和差异检验

平行测定样品脂肪酸含量并取均值,用SPSS 20.0对FAs含量按研究分类(油脂种类、胡麻油冷榨和热榨、胡麻油一级和二级)进行描述性统计和差异检验。结果用平均值±标准偏差表示。

1.4.3 相关性分析

用Excel对胡麻油5种关键FAs与菜籽油掺假量进行相关性分析。

1.4.4 偏最小二乘回归定量预测模型的建立

用SIMCA 14.1建立胡麻油掺入菜籽油的偏最小二乘回归(partial least squares,PLS)定量预测模型,估计FAs建模判别胡麻油掺菜籽油的敏感阈。利用相关系数(R2)、校正均方差(RMSEE)、交叉验证均方差(RMSEcv)评估PLS模型可靠性,R2的值越接近1,RMSEE和RMSEcv的值越小说明模型的准确性和精度越高,模型稳定性越好。

2 结果与分析

2.1 10种油脂脂肪酸指纹特征及真实性判别分析

10种油脂FAs指纹数据进行center scaling预处理后的PCA和OPLS-DA结果见图1。PCA和OPLS-DA前2个主成分方差累计贡献率为71.9%和71.0%,说明对12种FAs指标共性提取率较高。10种油脂以FAs指纹在主轴(横轴)上分开聚类,说明不同种类油脂FAs指纹差别大,且主要因素线性高。胡麻油聚类于主轴最右侧,与其他油脂距离较远,相对而言与菜籽油和大豆油聚类最靠近,与花生油最远,说明胡麻油与其他油脂FAs指纹极其不同,与菜籽油和大豆油FAs指纹最接近,也应该最难鉴别。菜籽油颜色比大豆油深,但与色浅的冷榨胡麻油很接近。并且菜籽油比胡麻油便宜得多,种植加工地域与胡麻油有重叠,我国东北、华北、西北都有菜籽种植和加工,菜籽油也最有可能作为胡麻油的掺假油。因此,本研究着重探究了FAs指纹判别胡麻油与菜籽油的可行性。

图1 10种油脂脂肪酸PCA 和 OPLS-DA 结果
Fig.1 PCA and OPLS-DA results for fatty acids of ten kinds of oil

调和油聚类较离散,这应该与不同品牌产品调配油种类不同有关。餐厨垃圾油和地沟油在第二轴(纵轴)方向与植物油完全分离,而两者聚类较靠近,说明两者本质上一致,都是来自剩饭菜的废弃油。OPLS-DA的聚类功能属于有监督的聚类分析,按已知分类拟合,因此聚类和判别效果优于PCA[15],不同种类油距离更大,同种油更聚集。OPLS-DA模型对不同油脂判别准确率见表2。OPLS-DA模型对胡麻油判别准确率为100.0%,说明以FAs区分胡麻油与其他油脂完全可行。对调和油误判率最高在预料之内。精炼地沟油(来自国家食品安全风险评估中心)全部被判别为餐厨垃圾油(本研究采集),这可能与2类样品采样不均匀有关;但两者都是餐厨废弃油,本应设为一类建模,因而并不算误判。

表2 10种油脂分类结果

Tab.2 Classification results of ten kinds of oil 份

样品预测结果胡麻油菜籽油葵花籽油大豆油玉米油花生油芝麻油调和油餐厨垃圾油地沟油样本数正确数正确率/%胡麻油600000000006060100.0菜籽油020000000002020100.0葵花籽油002000000002020100.0大豆油00019100000201995.0玉米油000020000002020100.0花生油000002000002020100.0芝麻油000000200002020100.0调和油002122300102000废弃油餐厨垃圾油000000002402424100.0地沟油0000000080800

10种油脂FAs质量比见表3。根向量分析[图1(b)、图1(d)]结合描述性统计结果(表3)表明:胡麻油聚类分离的关键FAs为C18:3ω3、C18:1ω9和C18:2ω6,其中C18:3ω3几乎可作为判别胡麻油真实性的单指标。但根据国家标准,亚麻籽油(胡麻油)中C18:3ω3的相对含量为45% ～ 70%[16],即掺入高达55%的低价食用油,C18:3ω3含量依然可以达到标准。因此,单指标的孤立分析不足以鉴别胡麻油与其他油脂,需要多变量指标进行分析。餐厨垃圾油和地沟油聚类分离的关键FAs为C14:0、C16:1和C18:1(9t),结合描述性统计,这3种FAs含量显著高于植物油,可能与废弃油含动物油脂有关。

表3 10种油脂FAs含量比较

Tab.3 Fatty acid content comparison of ten kinds of oil g·100g-1

脂肪酸w(油脂)胡麻油菜籽油葵花籽油大豆油玉米油花生油芝麻油调和油餐厨垃圾油地沟油胡麻/菜籽C14:00.04±0.02b0.06±0.02b0.07±0.02b0.08±0.02b0.04±0.02b0.04±0.02b0.02±0.02b0.08±0.06b0.64±1.10a0.55±0.74a0.7C16:05.04±0.96e4.18±0.98e5.77±2.68de10.00±0.98bc11.57±0.80bc9.35±1.86c8.54±0.76cd9.47±3.88c13.00±4.46b18.32±32.98a1.2C18:03.39±0.60b1.71±0.42c3.13±1.36b3.93±0.62b1.55±0.10c3.14±0.82b4.98±0.48a3.00±1.58b5.63±5.08a3.96±3.16b2.0C20:00.15±0.12f0.53±0.04b0.23±0.10ef0.36±0.10cd0.38±0.12c1.26±0.50a0.57±0.06b0.41±0.14c0.28±0.04de0.38±0.12c0.3C22:00.28±0.14e0.64±0.08d1.41±0.14b0.79±0.12cd0.25±0.06e4.77±1.34a0.29±0.08e1.01±0.78c0.58±0.16d0.32±0.50e0.4C24:00.13±0.16d0.22±0.24cd0.31±0.46bc0.29±0.06cd0.30±0.06bc2.30±0.92a0.18±0.06cd0.47±0.32b0.22±0.44cd0.17±0.20cd0.6C16:10.06±0.06c0.17±0.14c0.07±0.08c0.08±0.02c0.09±0.02c0.06±0.02c0.11±0.04c0.10±0.02c0.87±0.90a0.48±0.98b1.2C18:1(9t)0.06±0.18c0.05±0.06c0.02±0.04c0.05±0.08c0.05±0.04c0.05±0.04c0.19±0.20c0.05±0.04c0.33±0.60b0.51±1.30a1.2C18:1ω916.62±5.02f52.24±7.70a24.54±6.50d19.79±3.90e25.10±2.00d40.63±9.24b34.83±1.74c26.41±8.18d25.31±5.90d33.20±6.40c0.3C20:10.37±1.12cd0.56±1.42bc0.10±0.14d0.31±0.40cd0.22±0.26cd0.92±0.50ab0.16±0.14cd0.36±0.48cd0.37±0.20cd1.03±2.40a0.7C18:2ω614.62±4.64g19.47±3.28f57.23±8.28a49.47±1.98bc52.65±3.66ab31.14±8.14f43.14±2.82d47.36±9.50cd38.12±17.6e26.46±29.20f0.8C18:3ω350.14±9.68a7.99±2.50b0.14±0.08e7.00±3.10b0.71±0.24e0.10±0.08e0.37±0.06e3.72±4.56cd5.28±2.20bc1.89±4.20de6.3SFA9.24±1.26de7.50±1.32e11.50±0.58cd15.68±1.12b14.24±1.08bc20.83±4.28a14.70±0.80bc14.61±3.04bc20.98±11.18a24.12±34.64a1.2MUFA17.36±6.08e54.86±7.78a24.79±6.38d20.31±4.10e25.52±2.06d40.93±12.32b35.27±1.84c27.00±8.38d27.30±7.48d37.46±10.36bc0.3PUFA65.48±7.52a28.36±5.84e57.49±8.16b56.95±4.24bc53.49±3.82bc32.32±12.32e43.58±2.82d51.40±9.26c44.04±19.66d28.90±15.12e2.3TFA0.10±0.20c0.11±0.12c0.04±0.12c0.13±0.20c0.07±0.06c0.09±0.08c0.21±0.22bc0.11±0.16c0.35±0.60b0.62±1.40a0.9ω-3PUFA50.56±9.38a8.31±3.08b0.21±0.20f7.03±3.10bc0.73±0.24f0.13±0.08f0.38±0.06f3.74±4.56de5.37±2.16cd1.96±4.20ef6.1ω-6PUFA14.87±4.76h19.95±3.12g57.25±8.28a49.83±1.92bc52.73±5.68ab31.17±8.14f43.17±2.82de47.59±9.24cd38.66±17.78e26.80±29.34f0.7ω-6/ω-30.32.4272.67.172.2239.8113.612.77.213.7-

不同字母表示同行数据差异显著(P<0.05)或差异极显著(P<0.01)。

由10种油脂FAs质量比结果可知,胡麻油以不饱和脂肪酸为主,占总FAs质量的90%,多不饱和脂肪酸(polyunsaturated fatty acid,PUFA)占总FAs的71%。在测得的12种FAs中,以C18:3ω3、C18:1ω9、C18:2ω6、C16:0和C18:0为主,分别占总FAs质量的54%、18%、16%、5%和4%。胡麻油C18:3ω3质量比极显著高于其他油脂(P<0.01),分别是菜籽油、葵花籽油、大豆油、玉米油、花生油、芝麻油、调和油、餐厨垃圾油和地沟油的6.3、358.1、7.2、70.6、501.4、135.5、13.5、9.5和26.5倍;C18:1ω9和C18:2ω6质量比显著低于其他油脂(P<0.05)。菜籽油C18:1ω9,葵花籽油C18:2ω6及花生油C20:0、C22:0、C24:0质量比显著高于其他油脂(P<0.05)。餐厨垃圾油和地沟油C14:0、C16:1和C18:1(9t)质量比显著高于其他油脂(P<0.05)。大豆油、玉米油、芝麻油和调和油没有明显的特征FAs。

从营养学角度出发,胡麻油、菜籽油、葵花籽油、大豆油、玉米油、花生油、芝麻油、调和油、餐厨垃圾油和地沟油ω-6/ω-3为0.3、2.4、272.6、7.1、72.2、239.8、113.6、12.7、7.2和13.7,胡麻油和菜籽油比值最接近中国营养学会建议的1或2[17]。

2.2 冷榨和热榨胡麻油脂肪酸指纹特征及真实性判别分析

冷榨和热榨胡麻油FAs指纹数据进行center scaling预处理后的PCA和OPLS-DA结果见图2。PCA和OPLS-DA前2个主成分方差累计贡献率为92.3%和78.7%,说明对12种FAs指标共性提取率较高。OPLS-DA分离效果优于PCA,冷榨和热榨样品在OPLS-DA分析中可以完全分离[图2(c)],说明不同加工方法的胡麻油FAs存在差异,可能是高温热榨破坏部分FAs,导致聚类分离。此外,同种油脂样本点分布比较分散,可能是不同厂家或品牌的生产工艺参数存在差异造成的。OPLS-DA模型对冷榨和热榨胡麻油准确率见表4,冷热榨胡麻油整体判别准确率为96.7%,说明以FAs区分冷榨和热榨胡麻油可行。

图2 冷榨和热榨胡麻油脂肪酸PCA和OPLS-DA结果
Fig.2 PCA and OPLS-DA results for fatty acids of cold and hot pressed flaxseed oil

表4 冷榨和热榨胡麻油分类结果

Tab.4 Classification results of cold and hot pressed flaxseed oil 份

样品预测结果冷榨热榨样本数正确数正确率/%冷榨3603636100.0热榨222242291.7合计3822605896.7

冷榨和热榨胡麻油FAs含量比较见表5。根向量分析[图2(b)、图2(d)]结合描述性统计结果(表5)表明:冷榨和热榨胡麻油分类的关键FAs 为C18:3ω3、C18:2ω6和C18:1ω9。

表5 冷榨和热榨胡麻油FAs含量比较

Tab.5 Comparison of fatty acid content between cold and hot pressed flaxseed oil g·100g-1

脂肪酸w(冷榨)w(热榨)w(冷榨)/w(热榨)C14:00.04±0.020.04±0.021.0C16:04.91±0.62b5.25±1.22a0.9C18:03.49±0.54a3.24±0.56b1.1C20:00.12±0.08b0.20±0.12a0.6C22:00.26±0.12b0.30±0.16a0.9C24:00.17±0.04a0.07±0.18b2.4C16:10.06±0.040.06±0.101.0C18:1(9t)0.09±0.22a0.02±0.04b4.5C18:1ω915.60±2.62b18.15±6.14a0.9C20:10.41±1.280.30±0.861.4C18:2ω613.73±1.68b15.96±6.20a0.9C18:3ω353.21±5.32a45.54±7.00b1.2SFA9.21±1.189.30±1.381.0MUFA16.23±2.66b19.05±8.00a0.9PUFA67.23±5.30a62.85±7.34b1.1TFA0.13±0.26a0.05±0.08b2.6ω-3PUFA53.29±5.34a46.48±8.16b1.1ω-6PUFA13.91±1.86b16.31±6.18a0.9ω-6/ω-30.30.4-

不同字母表示同行数据差异显著(P<0.05)或差异极显著(P<0.01)。

由表5可知,冷榨与热榨胡麻油饱和脂肪酸(saturated fatty acid,SFA)差异不显著,冷榨PUFA显著高于热榨(P<0.05),热榨单不饱和脂肪酸(monounsaturated fatty acid,MUFA)显著高于冷榨(P<0.05),这是由于高温破坏部分PUFA,从而使MUFA相对含量增加,与糟帆等[18]研究结果一致。12种FAs中有9种FAs存在显著差异(P<0.05)或极显著差异(P<0.01),其中C18:3ω3、C18:1ω9和C18:2ω6的绝对差异较大。冷榨C18:3ω3质量比显著高于热榨(P<0.05),C18:1ω9和C18:2ω6显著低于热榨(P<0.05),这是由于C18:3ω3对热处理敏感,高温热榨时被破坏,从而使C18:1ω9和C18:2ω6相对含量增加[19]。

2.3 一级和二级胡麻油脂肪酸指纹特征及真实性判别分析

以我国油料品种、质量以及相应加工方法为指标,确定植物油的质量等级。我国食用油质量标准体系规定,市场上的一般食用油根据精炼程度(橄榄油和特种油脂除外)分为4个等级。国家标准将传统压榨胡麻油成品依据色泽、气味、透明度、水分及挥发物、不溶性杂质和酸价分为一级和二级[16]。

一级和二级胡麻油FAs指纹数据进行unit variance scaling预处理后的PCA和OPLS-DA结果见图3。PCA和OPLS-DA前2个主成分累计方差贡献率为81.9%和77.6%,说明对12种FAs指标共性提取率较高。OPLS-DA分离效果优于PCA,一级和二级胡麻油在OPLS-DA分析中可以完全分离[图3(c)],说明不同品级胡麻油FAs存在差异。不同品级胡麻油有不同加工方法,因此不同品级胡麻油聚类分离结果不能排除是否有加工方法的交互作用。OPLS-DA模型对一级和二级胡麻油判别准确率见表6,一级、二级胡麻油整体判别准确率为90.0%,说明以FAs区分一级和二级胡麻油可行。

图3 一级和二级胡麻油脂肪酸PCA和OPLS-DA结果
Fig.3 PCA and OPLS-DA results for fatty acids of first and second grade flaxseed oil

表6 一级和二级胡麻油分类结果

Tab.6 Classification results of first and second grade flaxseed oil 份

样品预测结果一级二级样本数正确数正确率/%一级382403895.0二级416201680.0合计4218605490.0

一级和二级胡麻油FAs比较见表7。根向量分析[图3(b)、图3(d)]结合描述性统计结果(表7)表明:一级和二级胡麻油分类的关键FAs为C18:3ω3、C18:1ω9和C18:2ω6。

表7 一级和二级胡麻油FAs含量比较

Tab.7 Comparison of fatty acid content between first and second grade flaxseed oil g·100g-1

脂肪酸w(一级)w(二级)w(一级)/w(二级)C14:00.04±0.020.04±0.021.0C16:04.97±0.545.19±1.880.9C18:03.48±0.56a3.21±0.52b1.1C20:00.12±0.10b0.21±0.08a0.6C22:00.25±0.12b0.32±0.10a0.8C24:00.15±0.12a0.08±0.18b1.9C16:10.06±0.040.06±0.101.0C18:1(9t)0.08±0.22a0.01±0.04b8.0C18:1ω915.73±2.64b18.38±6.64a0.9C20:10.37±1.240.35±0.901.1C18:2ω613.93±2.52b16.01±6.46a0.9C18:3ω352.53±6.52a45.38±7.66b1.2SFA9.24±1.149.25±1.501.0MUFA16.32±2.60b19.43±8.58a0.8PUFA66.75±6.02a62.93±7.74b1.1TFA0.13±0.24a0.04±0.06b3.3ω-3PUFA52.60±6.52a46.49±8.96b1.1ω-6PUFA14.11±2.62b16.39±6.46a0.9ω-6/ω-30.30.4-

不同字母表示同行数据差异显著(P<0.05)或差异极显著(P<0.01)。

由表7可知,一级胡麻油PUFA显著高于二级(P<0.05),二级MUFA显著高于一级(P<0.05),两者SFA差异不显著(P>0.05)。12种FAs中8种FAs存在差异显著(P<0.05)或差异极显著(P<0.01),其中C18:3ω3、C18:1ω9和C18:2ω6的绝对差异较大。一级C18:3ω3质量比显著高于二级(P<0.05),C18:1ω9和C18:2ω6质量比显著低于二级(P<0.05)。由于C18:3ω3不稳定极易被氧化,加工等级越高其含量应越少,而本结果则相反,可能与加工方法、胡麻籽品种和品质等多方面因素有关[19-20],还有待进一步扩大样本量深入研究。

2.4 胡麻油掺菜籽油判别分析

掺不同比例菜籽油的胡麻油FAs指纹数据进行unit variance scaling预处理后的OPLS-DA结果见图4。OPLS-DA前2个主成分方差累计贡献率为90.1%,说明对12种FAs指标共性提取率高。随着菜籽油质量分数的增加,样本呈现明显梯度分离,掺质量分数25%菜籽油的样本最靠近胡麻油,50%和75%的样品依次远离[图4(a)],说明随着菜籽油质量分数的增加,FAs含量有明显差异。OPLS-DA模型对纯胡麻油和掺入25%～100%菜籽油的胡麻油样品的二分类判别结果见表8,判别准确率平均为100%。

图4 纯胡麻油和掺菜籽油的胡麻油脂肪酸OPLS-DA结果
Fig.4 OPLS-DA results for fatty acids of pure flaxseed oil and flaxseed oil mixed with rapeseed oil

表8 掺不同比例菜籽油的胡麻油分类结果

Tab.8 Classification results of flaxseed oil mixed with different proportions of rapeseed oil 份

样品预测结果胡麻油胡麻油掺菜籽油样本数正确数正确率/%胡麻油1001010100胡麻油掺菜籽油0404040100合计10405050100

由图4(b)可知,变量X中C18:3ω3、C18:0和C16:0与变量Y中胡麻油的位置最接近,与掺25%、50%、75%和100%菜籽油的位置依次远离,说明这3种FAs含量随着菜籽油质量比增加而减少,其余FAs则相反。胡麻油中5种主要FAs质量比与掺菜籽油质量分数进行相关性分析,见图5。C18:3ω3 质量比与菜籽油质量分数呈极显著负相关,C16:0和C18:0与菜籽油质量分数呈负相关但不显著,C18:1ω9和C18:2ω6呈极显著正相关,与根向量分析结果一致。

图5 胡麻油5种关键FAs与掺菜籽油质量分数相关性分析
Fig.5 Correlation analysis of five key FAs of flaxseed oil and rapeseed oil mass fraction

胡麻油掺菜籽油PLS定量预测模型见图6。PLS模型的R2为0.986,RMSEE为4.263 5,RMSEcv为4.264 8,说明模型具有很好的预测能力。依据样品分布范围可确定模型鉴别掺假能力,掺假量预测误差不超过10%,即当掺假质量分数达10%以上,可准确判别纯胡麻油和掺假胡麻油。本研究建立模型的样本量少,有很大的优化空间,完全有望提高模型的精确度和准确度。

图6 胡麻油掺菜籽油的PLS定量预测模型
Fig.6 PLS quantitative prediction model of flaxseed oil mixed with rapeseed oil

3 结 论

对10种油脂(包括植物油和餐厨废弃油)的FAs指纹进行聚类分析,结果表明,油脂种类能被有效区分,尤其是胡麻油整体上与其他油脂聚类最远,与菜籽油和大豆油距离最近。OPLS-DA模型对胡麻油物种判别准确率为100.0%,对胡麻油冷热榨工艺判别准确率96.7%,对一级和二级胡麻油判别准确率90.0%。胡麻油物种、冷热榨工艺和品级鉴别的关键FAs为C18:3ω3、C18:1ω9和C18:2ω6。PLS定量分析表明,胡麻油掺入菜籽油的预测误差不超过10%。基于FAs指纹特征建立胡麻油真实性判别模型较为可行,并有继续优化的空间。

参考文献：

[1] YANG Y N, DENG Q C, JIA X, et al. Characterization of key odorants in peeled and unpeeled flaxseed powders using solvent-assisted flavor evaporation and odor activity value calculation[J]. LWT-Food Science and Techno-logy, 2021, 138: 110724.

[2] 中华人民共和国统计局.年度数据-胡麻籽产量[EB/OL].(2018-12-31)[2022-12-30]. https://data.stats.gov.cn/easyquery.htm?cn=C01&zb=A0D0F&sj=2021.

[3] AKRAMI A, NIKAEIN F, BABAJAFARI S, et al. Comparison of the effects of flaxseed oil and sunflower seed oil consumption on serum glucose, lipid profile, blood pressure, and lipid peroxidation in patients with metabolic syndrome[J]. Journal of Clinical Lipidology, 2018, 12(1): 70-77.

[4] TANG Z X, YING R F, LV B F, et al. Flaxseed oil:extraction, health benefits and products[J]. Quality Assurance and Safety of Crops &Foods, 2021, 13(1): 1-19.

[5] 王淑珍, 李应霞, 王兴瑞, 等. 不同制油方法对青海亚麻籽油品质及货架期的影响[J]. 食品科学, 2022, 43(11):92-98.

WANG S Z, LI Y X, WANG X R, et al. Effects of different preparation methods on quality and shelf life of Qinghai flaxseed oil[J]. Food Science, 2022, 43(11):92-98.

[6] 赵丹, 汪学德, 张润阳, 等. 制油工艺对油脂品质的影响研究[J]. 中国油脂, 2018, 43(6):11-15.

ZHAO D, WANG X D, ZHANG R Y, et al. Effect of oil preparation process on oil quality[J]. China Oils and Fats, 2018, 43(6):11-15.

[7] SUN X L, ZHANG B, HAN J J, et al. Effect of roasting temperature and time on volatile compounds, total tocopherols, and fatty acids of flaxseed oil[J]. Journal of Food Science, 2022, 87(4): 1624-1638.

[8] 郝梦. 胡麻油苯并芘污染水平的ELISA和HPLC法检测评价[D]. 呼和浩特: 内蒙古农业大学, 2018: 1-48.

HAO M. Detection and evaluation of benzopyrene pollution level of linseed oil by ELISA and HPLC[D]. Hohhot: Inner Mongolia Agricultural University, 2018: 1-48.

[9] ZENG J P, XIAO T, NI X G, et al. The comparative analysis of different oil extraction methods based on the quality of flaxseed oil[J]. Journal of Food Composition and Analysis, 2022, 107: 104373.

[10] 郭军, 杨月欣, 王竹, 等. 原奶或液态奶真实性鉴别方法和系统: CN105021736B[P]. 2017-11-07.

GUO J, YANG Y X, WANG Z, et al.Authenticity identifying method and system for raw milk or liquid milk: CN105021736B[P]. 2017-11-07.

[11] 杨惠茹, 郭军, 古斯愣图, 等. 六种家畜乳氨基酸特征聚类分析[J]. 食品与发酵工业, 2022, 48(6):233-240.

YANG H R, GUO J,GUSILENGTU, et al. Cluster analysis on amino acid properties of six livestock milk[J]. Food and Fermentation Industries, 2022, 48(6):233-240.

[12] WANG X R, HAN Y Z, LI Y X, et al. Detection of Qinghai-Tibet plateau flaxseed oil adulteration based on fatty acid profiles and chemometrics[J]. Food Control, 2021, 130: 108332.

[13] SUN X M, ZHANG L X, LI P W, et al. Fatty acid profiles based adulteration detection for flaxseed oil by gas chromatography mass spectrometry[J]. LWT-Food Science and Technology, 2015, 63(1): 430-436.

[14] 国家卫生和计划生育委员会, 国家食品药品监督管理总局. 食品安全国家标准 食品中脂肪酸的测定: GB 5009.168—2016[S]. 北京: 中国标准出版社, 2017.

[15] DOS SANTOS F A, SOUSA I P, FURTADO N A J C, et al. Combined OPLS-DA and decision tree as a stra-tegy to identify antimicrobial biomarkers of volatile oils analyzed by gas chromatography-mass spectrometry[J]. Revista Brasileira De Farmacognosia, 2018, 28(6): 647-653.

[16] 国家市场监督管理总局, 国家标准化管理委员会. 亚麻籽油: GB/T 8235—2019[S]. 北京: 中国标准出版社, 2019.

[17] 杨月欣, 葛可佑. 中国营养科学全书[M]. 2版. 北京: 人民卫生出版社, 2019: 58-75.

YANG Y X, GE K Y. Encyclopedia of nutrition science[M]. 2nd ed. Beijing: People’s Medical Publishing House, 2019: 58-75.

[18] 糟帆, 丁彩云, 马玉婷, 等. 不同制油工艺对亚麻籽油品质及抗氧化活性的影响[J]. 中国油脂, 2022, 47(9):13-18.

ZAO F, DING C Y, MA Y T, et al. Effects of different extraction processes on quality and antioxidant activity of linseed oil[J]. China Oils and Fats, 2022, 47(9):13-18.

[19] 余盖文, 史训旺, 洪梦佳, 等. 炒籽温度对压榨亚麻籽油品质的影响[J]. 中国油脂, 2019, 44(6):29-33.

YU G W, SHI X W, HONG M J, et al. Effects of frying temperature of flaxseed on quality of pressed flaxseed oil[J]. China Oils and Fats, 2019, 44(6):29-33.

[20] 缪纯庆, 刘秦, 范惠玲, 等. 张掖市不同海拔地区5个胡麻品种品质特性比较分析[J]. 生物技术进展, 2021, 11(3): 330-337.

MIAO C Q, LIU Q, FAN H L, et al. Analysis on quality character of five flaxseed varieties in different altitude regions of Zhangye[J]. Current Biotechnology, 2021, 11(3): 330-337.

Establishment of Authenticity Discrimination Model of Flaxseed Oil Based on Fatty Acid Fingerprint

YE Le, LI Ru, GUO Jun*

(College of Food Science and Engineering, Inner Mongolia Agricultural University, Hohhot 010018, China)

Abstract：Flaxseed oil is rich in α-linolenic acid, and its price is 3 to 20 times higher than that of other edible oils, and it is very vulnerable to be adulterated and counterfeited. The authenticity discrimination model of flaxseed oil based on fatty acid (FAs) fingerprints was established. Eight kinds of edible oil samples (flaxseed oil, rapeseed oil, sunflower oil, soybean oil, peanut oil, corn oil, sesame oil, blended oil), two kinds of waste oil samples (kitchen waste oil and gutter oil) were collected, and 12 FAs were determined by gas chromatography. Principal component analysis (PCA) was conducted to observe the clustering properties of these 10 kinds of oils according to their FAs fingerprint. Orthogonal partial least squares discriminant analysis (OPLS-DA) was applied to evaluate the feasibility of FAs fingerprint on authenticity of species origin, processing methods (cold or hot pressing) and grade of flaxseed oil, and partial least squares (PLS) quantitative prediction model was established to estimate the adulteration dosage of rapeseed oil in flaxseed oil. The FAs fingerprint of flaxseed oil was extremely different from other oils. Content of C18:3ω3[(50.14±9.68)g/100g] was significantly higher than that of other oils, while content of C18:1ω9[(16.62±5.02) g/100g] and C18:2ω6[(14.62±4.64)g/100g] were significantly lower than that of other oils. PCA and OPLS-DA both could distinguish species origin, processing methods and grade of flaxseed oil. OPLS-DA had better discrimination formance. PCA and OPLS-DA showed that flaxseed oil was the closest cluster to rapeseed and soybean oil. Accuracy of OPLS-DA for flaxseed oil species discrimination was 100%, accuracy of cold and hot pressing discrimination was 96.7%, and accuracy of grade discrimination was 90.0%. The key FAs for discrimination of flaxseed oil species, processing methods and grade were C18:3ω3, C18:1ω9 and C18:2ω6. Correlation analysis of five key FAs contents in flaxseed oil and adulteration dosage of rapeseed oil showed that C18:3ω3 content was negatively correlated with rapeseed oil adulteration dosage, C18:1ω9 and C18:2ω6 content were positively correlated with rapeseed oil adulteration dosage. R2 of PLS quantitative discrimination model for flaxseed oil adulterated with rapeseed oil was 0.986, RMSEE was 4.263 5, while RMSEcv was 4.264 8, and the prediction error of adulteration dosage was no more than 10%. The results showed that it was feasible to establish the authenticity discrimination model of flaxseed oil by FAs fingerprint, and there was optimization space of this mold.

Keywords：Linum usitatissimum L. oil; fatty acids; fingerprint profile; discrimination; authenticity

doi：10.12301/spxb202300049

文章编号：2095-6002(2023)06-0139-11

引用格式：叶乐,李茹,郭军.基于脂肪酸指纹的胡麻油真实性判别模型的建立[J]. 食品科学技术学报,2023,41(6):139-149.

YE Le, LI Ru, GUO Jun. Establishment of authenticity discrimination model of flaxseed oil based on fatty acid fingerprint[J]. Journal of Food Science and Technology, 2023,41(6):139-149.

中图分类号：TS225.1

文献标志码：A

收稿日期：2023-01-31

基金项目：国家自然科学基金资助项目(31760489)。

Foundation: National Natural Science Foundation of China(31760489).

第一作者：叶 乐,女,硕士研究生,研究方向为食品营养与安全。

*通信作者：郭 军,男,教授,博士,主要从事食品营养与安全方面的研究。

(责任编辑：张逸群)