橄榄油是以油橄榄油为原料经物理压榨得到的优质植物油,富含不饱和脂肪酸、多酚、角鲨烯以及维生素等营养物质[1],具有预防糖尿病、心脑血管疾病、抗衰老、防癌等功能[2],是地中海饮食的重要组成部分,在西方国家享有“液体黄金”之美誉[3]。虽然橄榄油的品质和等级主要依据其有机指标来确定[4],但其所含无机组分对于橄榄油的食用安全和保质期起重要的作用[5],同时,橄榄油中无机元素的含量不仅受油橄榄对矿物元素的吸收和预富集特性所控制,而且与油橄榄的产地和生长环境密切相关[6]。因此,通过研究橄榄油中无机元素的种类和含量,并结合数据统计分析表征橄榄油的地域信息,对橄榄油的产地溯源保护和地理标志保护具有重要意义。
目前,植物油样品的前处理技术主要有干法灰化、湿法消解、微波消解和有机溶剂稀释等方法[7-8]。干法灰化法和湿法消解法在加热处理过程会导致部分元素挥发,不能用于挥发性元素的测定;微波消解法样品消解彻底,但由于消解油脂的特殊性,长时间预消解操作费时,不适合高通量分析[9];有机溶剂稀释法避免了复杂烦琐的消解过程和分析元素的损失,但必须面对有机基质在分析过程中形成的严重干扰。在应用于植物油中微量元素的原子光谱测定技术中,电感耦合等离子体串联质谱(inductively coupled plasma tandem mass spectrometry,ICP-MS/MS)利用特有的串联质谱和碰撞反应池(collision reaction cell,CRC)系统,几乎可以消除所有干扰,适合微量以及痕量元素的测定[10-12]。因此,采用有机溶剂稀释结合ICP-MS/MS能准确测定植物油中的微量元素。
主成分分析是通过对错综复杂数据进行降维,筛选能重现的数据信息,将原变量转换为少数新变量(主成分)的多元统计分析方法,主成分分析能集中典型地表征原变量的数据特征,从而充分反映总体信息[13]。聚类分析是依据观察对象某些特征加以归类的数理统计方法[14]。主成分分析和聚类分析已广泛应用于食品分类和产地溯源的研究[15-18]。本文采用煤油稀释橄榄油后直接进样,利用ICP-MS/MS测定其中16种微量元素的含量,并根据不同产地橄榄油中微量元素的含量进行主成分分析和聚类分析,以期为橄榄油的产地溯源保护和地理标志保护提供理论依据。
18种特级初榨橄榄油样品(产地分别为意大利4种、希腊4种、西班牙4种、中国6种)购置于京东超市。
1 000 μg/g的Na、K、Ti、Cr、Mn、Fe、Co、Ni、Cu、Zn、As、Se、Sr、Cd、Hg、Pb单元素油标,美国Spex CertiPrep公司;1 000 μg/g的Sc、Y、Bi单元素内标,美国Spex CertiPrep公司;高纯煤油,国药集团化学试剂有限公司;99.99%的高纯O2、NH3、He,湖南省华中特种气体有限公司产品。
Agilent 8800型电感耦合等离子体串联质谱仪,美国Agilent公司;配置有自吸式C-flow 200 PFA雾化器,Peltier冷却双通道Scott型雾化室,窄口径(内径为1.5 mm)进样炬管,S透镜,Pt锥。Element XR型扇形磁场电感耦合等离子质谱(sector field inductuvely coupled plasma mass spectormetry,SF-ICP-MS)仪,美国Thermo Fisher Scientific公司。
1.3.1 橄榄油前处理
采用质量稀释法用煤油将橄榄油样品进行5倍稀释后待测,同时采用煤油为空白试剂。
1.3.2 ICP-MS/MS工作参数
ICP-MS/MS调谐优化后的操作条件:射频功率1 550 W;载气流速1.08 L/min;辅助气流速1.08 L/min;补偿气流速0.15 L/min;采样深度8.0 mm;O2流速0.08 L/min;操作模式MS/MS;反应气O2、NH3/He(He中混入10%的NH3);反应气流速,O20.30 mL/min、NH3/He 3.0 mL/min;八极杆偏置电压-18 V;动能歧视电压-8 V;同位素23Na、39K、48Ti、52Cr、55Mn、56Fe、59Co、60Ni、63Cu、66Zn、75As、78Se、88Sr、111Cd、202Hg、208Pb。
1.3.3 橄榄油微量元素的测定
采用煤油将1 000 μg/g的分析元素油标和内标进行稀释,配制0.00、2.00、10.、50.0、200 μg/kg系列分析元素混合油标,使用标准内标混合“T”型接头在线加入1 μg/g的Sc、Y、Bi混合内标,在优化的ICP-MS/MS操作条件下进行测定,获得各元素的校准曲线,同时对空白和样品进行测定,利用MassHunter工作站软件计算样品中分析元素的含量。
1.3.4 统计分析
利用18种橄榄油样品中16种微量元素的含量组成18×16的原始数据矩阵,应用SPSS 25软件进行统计分析。将原始数据进行标准化处理,以消除橄榄油中微量元素原始数据之间的量纲和数量级不同对分析结果的影响,使标准化后的数据具有可比性并遵循正态分布规律,从而进行主成分分析和聚类分析。
2.1.1 线性关系与检出限
按实验条件采用ICP-MS/MS测定系列分析元素混合油标,以分析元素与内标元素的信号强度比值所对应的标准溶液浓度自动进行线性回归,建立校准曲线,各元素的线性范围和相关系数见表1。采用所建立的方法对空白(煤油)连续测定10次,计算各元素标准偏差,以3倍标准偏差所对应的浓度为仪器的检出限(limit of detection,LOD)[19]。由表1可以看出,所有分析元素在具有良好的线性关系,线性相关系数≥0.999 7。各元素仪器的LOD为0.001~1.22 μg/kg。
表1 分析元素的线性关系与检出限
Table 1 Linear relationship and the limit of detection for analytes
元素线性范围/(μg·kg-1)线性相关系数(R2)LOD/(μg·kg-1)元素线性范围/(μg·kg-1)线性相关系数(R2)LOD/(μg·kg-1)Na4.10~2000.999 91.22Cu0.11~2000.999 80.034K0.48~2001.000 00.14Zn0.52~2000.999 70.16Ti0.02~2000.999 90.006As0.007~2000.999 90.002Cr0.013~2000.999 80.004Se0.018~2001.000 00.005Mn0.010~2000.999 90.003Sr0.10~2000.999 90.031Fe8.32~2001.000 00.25Cd0.034~2001.000 00.010Co0.003~2001.000 00.001Hg0.038~2000.999 90.011Ni0.04~2000.999 80.012Pb0.026~2001.000 00.008
2.1.2 方法准确性与精密度
由于缺少橄榄油标准参考物质,本实验选取实际样品中加入油标,采用所建立的分析方法对加标样品平行测定6次,计算各元素的加标回收率和相对标准偏差(relative standard deviation,RSD),同时采用SF-ICP-MS对加标样品进行分析,对2种分析方法的测定结果进行统计学分析,用来评价ICP-MS/MS的准确性和精密度,结果见表2。由表2可以看出,各元素的加标回收率为95.2%~108%,RSD≤5.3%;16种分析元素的加标测定值与采用SF-ICP-MS的加标测定值基本一致,验证了所建立的方法准确性和精密度好。
表2 加标样品的分析结果(n=6)
Table 2 Results obtained for analytes in spiked sample (n=6)
元素测定值/(μg·kg-1)加标值/(μg·kg-1)加标测定值/(μg·kg-1)加标回收率/%RSD/%SF-ICP-MS加标测定值/(μg·kg-1)Na305300621±27.81054.5610±19.0K97.2100205±7.531083.8212±7.13Ti0.281.001.26±0.0598.04.01.18±0.07Cr110100213±11.01035.2206±12.5Mn4.695.009.50±0.4896.25.19.81±0.52Fe233300519±16.195.33.1532±18.4Co0.321.001.38±0.041062.91.40±0.05Ni0.461.001.51±0.071054.61.55±0.09Cu1.835.006.93±0.2697.03.87.12±0.41Zn2.975.008.12±0.291033.68.00±0.37As2.035.006.79±0.3095.24.46.85±0.30Se0.721.001.69±0.0997.05.31.59±0.10Sr4.765.009.53±0.3495.43.69.48±0.38Cd0.601.001.56±0.0396.01.91.54±0.05Hg0.391.001.35±0.0596.03.71.27±0.04Pb0.581.001.62±0.041042.51.65±0.07
采用ICP-MS/MS法测定了来自国内外不同产地的18种橄榄油样品,其中来自意大利4种(样品编号IT1-4)、希腊4种(样品编号GR1-4)、西班牙4种(样品编号ES1-4)、中国6种(样品编号CN1-6),测定各橄榄油中的16种微量元素,每个样品重复测定6次,结果见表3。18种橄榄油中Na、K、Fe、Cr的含量相对较高,而Co、Cd、Hg、Pb的含量相对较低。对照我国橄榄油国家标准GB 23347—2009以及食品中污染物限量国家标准GB 2762—2017[20],所有橄榄油中的Fe、Cu的含量均低于GB 23347—2009中特级初榨橄榄油的限量值,重金属元素的含量远远低于GB 2762—2017中的限量标准。
表3 不同产地橄榄油中微量元素的含量(已进行了稀释校正) 单位:μg/kg
Table 3 Concentrations of trace elements in olive oil from different producing areas (Dilution correction has been carried out)
样品NaKTiCrMnFeCoNiCuZnAsSeSrCdHgPbIT13121790.2222.41.891520.132.306.383.010.334.108.590.450.470.64IT23302100.5123.15.291730.300.316.103.890.405.666.320.430.930.30IT32481890.4119.55.502140.452.185.062.140.731.370.420.600.200.17IT42272030.3824.94.362610.571.174.651.701.680.841.770.380.510.36GR130597.20.281104.692330.320.461.832.972.030.724.760.600.390.58GR22751150.1796.03.452780.410.251.461.822.450.621.920.310.280.35GR32791940.331023.563260.211.284.442.550.821.801.620.510.760.25GR42611700.2575.02.662730.181.163.972.901.362.780.350.820.340.37ES126487.20.1920.64.332000.501.434.572.821.741.850.680.260.350.23ES219790.40.2818.74.621950.471.121.821.332.801.040.550.780.460.21ES32101060.361233.561340.551.351.251.641.871.050.820.870.700.15ES41761200.1688.63.741860.601.401.521.792.100.990.720.480.360.52CN15131461.2931.14.213230.352.3112.68.500.162.0821.20.250.660.18CN25341611.9037.20.943570.513.2613.47.210.316.8222.30.710.650.48CN34352031.6427.81.362750.323.267.285.860.431.9417.60.220.540.21CN44621581.1530.21.742960.172.1411.36.920.303.1616.10.340.310.40CN54152390.4425.11.352300.352.256.505.630.573.4026.40.300.620.29CN64052251.7821.81.653480.122.209.875.282.533.9015.00.290.720.69
2.3.1 橄榄油中微量元素的适用性检验
采用SPSS 25的Bartlett's球状检验对18种橄榄油中的16种微量元素进行适用性检验,考察原始数据矩阵中各变量的相关性,当显著性值(Sig或P)<0.05时,原始数据矩阵中各变量具有相关性,适合主成分分析,反之不适合主成分分析[21]。本实验检验结果表明,原假设原始数据矩阵是单位矩阵时的显著性值为0.000,拒绝原始数据矩阵中各变量全独立的原假设,表明18种橄榄油中的16种微量元素之间具有强线性相关性,能应用主成分分析进行研究。
2.3.2 主成分的贡献率与选择
采用主成分分析法提取公因子,其特征值和总体方差描述见表4。由表4可以看出,累计方差贡献率的80.961%来源于前5个主成分,即1个5因子模型可解释80.961%的实验数据,包含了18种橄榄油中16种微量元素的大部分信息,并且这5个主成分的特征值均>1,符合累积方差贡献值大于80%且特征值>1的主成分个数选择原则[22]。因此,本实验提取前5个主成分来评判16种微量元素在18种橄榄油中的分布。
表4 主成分分析特征值和总体方差描述
Table 4 Principal component analysis characteristics and total variance explained
主成分特征值方差/%累积方差/%主成分特征值方差/%累积方差/%17.77348.58348.58390.3622.26496.13321.5699.80958.392100.2451.53097.66431.4058.78167.173110.1861.16598.82941.2057.53474.707120.0990.62199.45051.0016.25580.961130.0410.25699.70660.7844.90285.864140.0310.19199.89770.6744.21490.078150.0140.09099.98780.6073.79293.869160.0020.013100.000
2.3.3 微量元素的因子分析
根据橄榄油中微量元素的因子分析结果(见表5),分别建立5个主成分因子(F1、F2、F3、F4、F5)得分模型:
表5 橄榄油中微量元素的因子分析
Table 5 Factor analysis of trace elements in olive oil
元素主成分因子元素主成分因子F1F2F3F4F5F1F2F3F4F5Na0.945-0.0390.1530.0910.116Cu0.942-0.0990.0790.036-0.001K0.5950.037-0.524-0.2630.069Zn0.926-0.1140.1710.0810.124Ti0.8530.0280.2610.1980.079As-0.6670.3670.308-0.0210.132Cr-0.5210.3860.1550.4700.296Se0.7130.140-0.4340.235-0.252Mn-0.687-0.398-0.164-0.0680.290Sr0.909-0.0690.1010.0450.007Fe0.6230.2390.3820.0470.429Cd-0.4310.140-0.1020.697-0.365Co-0.445-0.5200.3260.290-0.126Hg0.379-0.033-0.5800.4390.404Ni0.761-0.1900.2450.048-0.407Pb0.1960.832-0.007-0.173-0.215
F1=0.945X1+0.595X2+0.853X3-0.521X4-0.687X5+0.623X6-0.445X7+0.761X8+0.942X9+0.926X10-0.667X11+0.713X12+0.909X13-0.431X14+0.379 X15+0.196X16
F2=-0.039X1+0.037X2+0.028X3+0.386X4-0.398X5+0.239X6-0.520X7-0.190X8-0.099X9-0.114X10+0.367X11+0.140X12-0.069X13+0.140X14-0.033X15+0.832X16
F3=0.153X1-0.524X2+0.261X3+0.155X4-0.164X5+0.382X6+0.326X7+0.245X8+0.079X9+0.171X10+0.308X11-0.434X12+0.101X13-0.102X14-0.580 X15-0.007X16
F4=0.091X1-0.263X2+0.198X3+0.470X4-0.068X5+0.047X6+0.290X7+0.048X8+0.036X9+0.081X10-0.021X11+0.235X12+0.045X13+0.697X14+0.439 X15-0.173X16
F5=0.116X1+0.069X2+0.079X3+0296X4+0.290X5+0.429X6-0.126X7-0.407X8-0.001X9+0.124X10+0.132X11-0.252X12+0.007X13-0.365X14+0.404X15-0.216X16
可以看出,F1在X1、X3、X9、X10、X13有较大的因子载荷值,分别为0.945、0.853、0.942、0.926、0.909,与对应的Na、Ti、Cu、Zn、Sr高度正相关,F2在X16有较大因子荷载值(0.832),与对应Pb高度正相关,由于F1和F2的累积方差贡献率(58.392%)大于50%,可确定Na、Ti、Cu、Zn、Sr、Pb是橄榄油的特征元素。
通过ICP-MS/MS测定得到橄榄油中微量元素含量的原始数据较为离散,不能直观地比较各元素的相近程度。本实验采用SPSS 25的系统聚类分析对不同产地橄榄油样品进行聚类分析。以16种微量元素含量为变量,对数据进行标准处理后,以欧氏距离平方(系统默认)为准则采用组间连接聚类,结果如图1所示。可以看出,随着聚类距离的增大,聚类结果也逐渐发生变化。当聚类距离为10时,18种橄榄油聚为两大类,国内样品聚为一类,国外样品聚为一类;当聚类距离为8时,18种橄榄油聚为四大类,国内样品聚为一类,国外样品意大利产橄榄油归为一类,希腊产橄榄油归为1类,西班牙产橄榄油归为一类。表明以橄榄油中16种微量元素为变量可以对不同产地的橄榄油进行溯源,有利于橄榄油的产地溯源保护和地理标志保护。
图1 橄榄中微量元素含量系统聚类
Fig.1 Phylogenetic clustering of olive oil based on trace elements
利用ICP-MS/MS测定了不同产地橄榄油中16种微量元素Na、K、Ti、Cr、Mn、Fe、Co、Ni、Cu、Zn、As、Se、Sr、Cd、Hg、Pb的含量,采用SPSS 25对18种橄榄油中16种微量元素组成的原始数据矩阵进行了主成分分析和聚类分析。橄榄油中的Na、K、Fe、Cr含量相对较高,Co、Cd、Hg、Pb含量相对较低,其余微量元素无明显规律性。不同产地橄榄油中微量元素的含量具有地域特征,通过聚类分析可以较好地区分橄榄油的产地来源,研究结果为橄榄油的产地判别提供了准确可靠的标记验证方法。
[1] KALUA C M, ALLEN M S, BEDGOOD D R, et al. Olive oil volatile compounds, flavour development and quality: A critical review [J]. Food Chemistry, 2007, 100(1): 273-286.
[2] 刘少敏, 何天鹏, 薛丹丹, 等. 六种特级初榨橄榄油挥发性香气成分的比较[J].食品与发酵工业, 2018, 44(6): 251-256.
[3] 邓龙, 陈雄飞, 刘贤标, 等. GC-MS-O结合电子鼻对橄榄油挥发性成分的分析与鉴别[J].中国食品学报, 2019, 19(5): 276-286.
[4] 中华人民共和国国家质量监督检验检疫总局. GB 23347—2009 橄榄油、油橄榄果渣油[S].北京: 中国标准出版社, 2009.
[5] CHOE E, MIN D B. Mechanisms and factors for edible oil oxidation [J]. Comprehensive Reviews in Food Science and Food Safety, 2006, 5(4): 169-186.
[6] ZEINER M, STEFFAN I, CINDRIC I J. Determination of trace elements in olive oil by ICP-AES and ETA-AAS: A pilot study on the geographical characterization [J]. Microchemical Journal, 2005, 81(2): 171-176.
[7] 邱会东, 赵波, 张红, 等. 食用植物油中重金属分析方法的研究进展[J]. 中国油脂, 2017, 42(1): 76-79.
[8] 张萍, 谢华林, 朱乾华, 等. 食用橄榄油中重金属元素的质谱分析[J]. 现代食品科技, 2014, 30(3): 206-209.
[9] 倪张林, 汤富彬, 屈明华. 不同前处理方法测定植物油中重金属的研究[J]. 中国油脂, 2012, 37(7): 85-87.
[10] GONG J, SOLIVIO M J, MERINO E J, et al. Developing ICP-MS/MS for the detection and determination of synthetic DNA-protein crosslink models via phosphorus and sulfur detection [J]. Analytical Bioanalytical Chemistry, 2015, 407: 2 433-2 437.
[11] NELSON J, HOPFER H, SILVA F, et al. Evaluation of GC-ICP-MS/MS as a new strategy for specific heteroatom detection of phosphorus, sulfur, and chlorine determination in foods [J]. Journal of Agricultural and Food Chemistry, 2015, 63(18): 4 478-4 483.
[12] HU X, CAO Z, SUN W, et al. Accurate determination of arsenic and selenium in plant food samples by using ICP-MS/MS [J]. Analytical Methods, 2016, 8(32): 6 150-6 157.
[13] WOLD S, ESBENSEN K, GELADI P. Principal component analysis [J]. Chemometrics and Intelligent Laboratory Systems, 1987, 2(1-3): 37-52.
[14] 倪永年. 化学计量学在分析化学中的应用[M]. 北京: 科学出版社, 2004: 228-229.
[15] LI D, MENG X, LI B. Profiling of anthocyanins from blueberries produced in China using HPLC-DAD-MS and exploratory analysis by principal component analysis [J]. Journal of Food Composition and Analysis, 2016, 47: 1-7.
[16] 郑彦婕, 胡书玉, 黎永乐, 等. 基于无机元素含量的地理标志食醋分类[J]. 食品与发酵工业, 2012, 38(9): 167-169.
[17] RANAMUKHAARACHCHI S A, PEIRIS R H, MORESOLI C. Fluorescence spectroscopy and principal component analysis of soy protein hydrolysate fractions and the potential to assess their antioxidant capacity characteristics [J]. Food Chemistry, 2017, 217: 469-475.
[18] MAIA M, NUNES F M. Authentication of beeswax (Apis mellifera) by high-temperature gas chromatography and chemometric analysis [J]. Food Chemistry, 2013, 136(2): 961-968.
[19] NIE X, LIANG Y, XIE H, et al. Trace amounts of impurities in electrolytic manganese metal by sector field inductively coupled plasma mass spectrometry [J]. Journal of Central South University, 2013, 20(12): 3 385-3 390.
[20] 中华人民共和国国家卫生部. GB 2762—2017 食品安全国家标准食品中污染物限量[S]. 北京: 中国标准出版社, 2017.
[21] PERES-NETO P R, JACKSON D A, SOMERS K M. How many principal components? Stopping rules for determining the number of non-trivial axes revisited [J]. Computational Statistics & Data Analysis, 2005, 49(4): 974-997.
[22] 刘晓燕, 李洪怡, 苏燕, 等. 臭氧和超声波对鲜切莲藕品质影响的主成分分析[J]. 食品与发酵工业, 2018, 44(5): 148-155.