基于气相色谱-离子迁移谱结合多元统计方法分析辣椒品种对辣椒油理化性质和风味物质的影响

杨芳,袁海彬,贾洪锋*,邓凤琳,王珍妮

(四川旅游学院 食品学院,四川 成都,610100)

摘 要 该研究对5种不同品种辣椒(子弹头、龙脊椒、涮涮辣、七星椒、干线椒)制成的辣椒油(分别编号为PZ1、PZ2、PZ3、PZ4、PZ5)中辣椒素类物质含量、色差、过氧化值、挥发性有机化合物(volatile organic compounds,VOCs)的影响,并结合多元统计方法分析5种辣椒油之间的关键差异标志物。结果表明,PZ3的辣椒素含量、辣度、亮度、黄度均最高,过氧化值最低且与其余样品间有显著性差异(P<0.05);经气相色谱-离子迁移谱(gas chromatography-ion mobility spectrometry,GC-IMS)共分离鉴定出73种VOCs,包括醛类、酯类、酮类、醇类、杂环类、羧酸类、硫醚类,分别有19、18、12、12、6、4、2种;GC-IMS二维差异图、指纹图谱及主成分分析结果一致,PZ3风味最为独特,PZ5和PZ4次之,PZ1和PZ2风味接近;通过偏最小二乘-判别分析发现,PZ3的特征香气物质主要为酯类,包括异戊酸己酯-D、异丁酸己酯、己酸异戊酯、丁酸己酯-D等;通过变量重要性投影分析,筛选出5种辣椒油的26种差异标志物,包括乙酸乙酯、异戊酸己酯、二乙醇缩乙醛、乙酸-D、己醛-D、二甲基硫醚、3-甲基丁醛、二氢-2(3H)-呋喃酮、1-羟基-2-丙酮-D、丙酮、丁酸-M等,其中异戊酸己酯、丁酸己酯-M、糠醇-M、己酸异戊酯-M、乙酸-D、3-甲基丁醛、戊醛-D等物质在第1主成分的载荷值较大;丙酸乙酯、3-甲基-1-丁醇、乙酸乙酯、丁醛、二乙醇缩乙醛等在第2主成分上的载荷较大。

关键词 辣椒油;气相色谱-离子迁移谱;辣椒素;挥发性风味物质;偏最小二乘-判别分析;变量重要性投影

气相色谱-离子迁移谱(gas chromatography-ion mobility spectrometry, GC-IMS)是一种近年来出现的新型气相分离和检测技术,具有分辨率高、灵敏度高等特点,特别适合于挥发性有机化合物(volatile organic compounds, VOCs)的痕量检测,在果蔬[1]、香辛料[2-3]、调味油[4]和菌类[5]等VOCs的检测分析方面得到了广泛的应用[6]

辣椒(Capsicum annuum L.)属茄科一年或多年生植物。辣椒中含有丰富的辣椒素、辣椒红素、维生素C、多酚、类黄酮等物质[7-9],因而具有清除自由基、抗氧化、抗肿瘤等生物活性[8-12]。不同品种的辣椒在辣椒素含量、VOCs的种类及含量上有明显差异[13-14]。辣椒油是用干辣椒和植物油制成的有独特风味的调味油,具有香味浓郁、辣度适口的特点,在川菜的加工中具有举足轻重的地位,深受人们的喜爱[15]。张洪新等[16]采用GC-MS对14种不同品种辣椒制成的辣椒油制品的挥发性香气成分、辣椒素含量等进行了对比,结果发现样品的香气成分和辣椒素含量均存在差异。何小龙等[17]对4种辣椒和大豆油制备的辣椒油进行研究发现,辣椒品种对辣椒油中辣椒素和二氢辣椒素的含量影响显著。石自彬等[18]对10种辣椒和花生油制备的辣椒油进行研究发现,辣椒油的挥发性风味成分中,酯类化合物含量最高;辣椒油的抗氧化效果随贮藏时间延长而减弱,反复熬煮会使辣椒油的抗氧化效果减弱。然而,上述研究主要采用常规的GC-MS技术,对一些含量少但对风味至关重要的物质却无法分析鉴定[4,19]。目前,采用GC-IMS分析辣椒品种对辣椒油中VOCs的影响的研究较少。

本研究通过HPLC和GC-IMS等检测技术和方法,对5种辣椒(子弹头、龙脊椒、涮涮辣、七星椒、干线椒)制备的辣椒油中辣椒素类物质含量、色差、过氧化值、VOCs进行检测和分析,结合主成分分析(principal component analysis, PCA)、偏最小二乘-判别分析(partial least squares-discrimination analysis, PLS-DA)及热图聚类分析等多元统计方法,分析辣椒品种对辣椒油的理化性质和VOCs的影响,以对辣椒油的理化性质和风味品质进行更为科学、系统的评价,这将为不同品种辣椒制备辣椒油的品质评定和合理加工提供一定的理论支撑。

1 材料与方法

1.1 材料与试剂

辣椒品种:子弹头、龙脊椒、涮涮辣、七星椒、干线椒;金龙鱼玉米油,益海嘉里(重庆)粮油有限公司;精制食用盐,四川省盐业总公司。

冰乙酸、二氯甲烷、KI、Na2S2O3、可溶性淀粉、KIO3,均为国产分析纯。

1.2 仪器与设备

Flavor Spec®型GC-IMS仪:含CTC自动顶空进样器、Laboratory Analytical Viewer(LAV)分析软件、GC×IMS Library Search软件及软件内置的NIST数据库和IMS 数据库对物质进行定性分析,德国G.A.S公司;戴安Ultimate 3000型HPLC仪,赛默飞世尔科技有限公司;IS128型万分之一分析天平,上海西塘生物科技有限公司。

1.3 试验方法

1.3.1 辣椒油制备

参考文献[17]的方法并稍作修改。

不同品种干红辣椒→微波(功率1 000 W)加热1 min→凉至室温→粉碎→过筛(18目)→称重(24 g)→加入食盐(1 g)→混匀辣椒面,备用。玉米油(100 g)→熬制(180 ℃)→加入至辣椒面中→搅拌(约35 s)→自然降温浸提(24 h)→编号(样品编号PZ1、PZ2、PZ3、PZ4、PZ5,分别对应子弹头、龙脊椒、涮涮辣、七星椒、干线椒5种辣椒制备的辣椒油样品)→待测

自然浸提完成后,每个样品取上层清油液,分别用于辣椒素类物质含量、色差值、过氧化值、VOCs的检测,各平行测定3次。

1.3.2 辣椒素类物质含量的测定

参照GB/T 21266—2007《辣椒及辣椒制品中辣椒素类物质测定及辣度表示方法》[20],采用高效液相色谱法对辣椒油中辣椒素含量、二氢辣椒素含量、辣椒素类物质总量、斯科维尔指数和辣度进行测定,平行测定3次。

色谱柱Zorbax SB-C18(4.6 mm×250 mm×5 μm)石英毛细管柱;流动相65%的甲醇水溶液;紫外光波长280 nm;流速1 mL/min;柱温30 ℃。

1.3.3 色差的测定

利用色差仪测定不同样品的亮度(L*)、红度(a*)、黄度(b*)值,平行测定3次。

1.3.4 过氧化值的测定

参照GB 5009.227—2016《食品安全国家标准 食品中过氧化值的测定》中的滴定法[21],对辣椒油的过氧化值进行测定,平行测定3次。

1.3.5 VOCs的检测

取1.5 g辣椒油样品于20 mL顶空瓶中,80 ℃下孵化20 min后,顶空进样,用GC-IMS仪进行测定;分析时间40 min;平行测定3 次,盲测1次。

自动进样条件:孵化温度80 ℃;孵化时间20 min;进样方式顶空进样;进样体积500 μL;进样针温度85 ℃;振荡方式加热;孵化转速500 r/min;不分流;清洗时间5 min。

GC条件:色谱柱WAX(30 m×0.53 mm×1 μm)石英毛细管柱;色谱柱温度60 ℃;载气N2(纯度≥99.999%);IMS温度45 ℃;载气流速0~10 min,2 mL/min;10~20 min,10 mL/min;20~40 min,100 mL/min。

IMS条件:漂移管长度9.8 cm,漂移管温度45 ℃,漂移气为N2(纯度≥99.999%),漂移气流速150 mL/min。

1.3.6 数据处理

利用Flavor Spec®配备的 Laboratory Analytical Viewer(LAV)分析软件及GC×IMS Library Search定性软件对辣椒油样品的VOCs进行采集和分析;利用软件内置的 NIST数据库和IMS数据库对物质进行定性分析。

采用Excel和Origin进行数据的统计和绘图;采用SPSS 26.0对样品间差异采用ANOVA检验分析,利用Duncan方法进行显著性差异分析,P<0.05,差异显著;采用SIMCA 14.1进行PLS-DA;用R语言进行热图聚类分析;结果用平均值±标准差的形式表示。

2 结果与分析

2.1 辣椒品种对辣椒油中辣椒素类物质含量的影响

辣椒素类物质是辣椒碱类化合物的总称,是辣椒油中的辛辣成分,属于酰胺类化合物[22]。辣椒素具有调味、清除自由基、抗氧化、抗癌等活性[8-12],但大量摄入辣椒素会产生毒性,局部刺激导致呼吸道感染,并可能会增加某些癌症的发生[12]。5种辣椒油的辣椒素类化合物含量如表1所示。

表1 辣椒油样品的辣椒素类物质含量
Table 1 Content of capsaicinoids in chili oil samples

编号辣椒素含量/(g/kg)二氢辣椒素含量/(g/kg)辣椒素类物质总量/(g/kg)斯科维尔指数/SHU辣度PZ10.15±0.00e0.08±0.00e0.25±0.00e3 835.67±66.15e 25.57±0.44ePZ20.28±0.00c0.13±0.00c0.46±0.01c7 098.67±92.93c47.32±0.62cPZ36.66±0.01a2.70±0.00a10.41±0.02a160 471.00±286.70a1 069.81±1.91aPZ40.51±0.01b0.34±0.00b0.95±0.01b14 597.33±79.05b97.32±0.53bPZ50.22±0.00d0.10±0.00d0.36±0.00d5 571.67±23.35d37.14±0.16d

注:同列小写字母不同表示差异显著(P<0.05)(下同)。

由表1可知,5种辣椒油的辣椒素含量、二氢辣椒素含量、辣椒素类物质总量、斯科维尔指数、辣度顺序均为PZ3>PZ4>PZ2>PZ5>PZ1;均为涮涮辣制备的辣椒油样品(PZ3)显著高于其余4种样品(P<0.05)。PZ3辣椒素类物质总量高达(10.41±0.02) g/kg、辣度高达1 069.81±1.91,这与辣椒素等辣味物质的含量多少直接影响辣度[23]的结论一致。产生这些差异的原因,可能主要是辣椒品种的不同。因此,在选择制作辣椒油的辣椒品种时,应根据产品对辣椒素类物质总量和辣度的不同需求来选择。

2.2 辣椒品种对辣椒油色差的影响

由图1可知,5种辣椒油的亮度(L*)依次为PZ3>PZ2>PZ5>PZ4>PZ1,PZ3最大,为45.57±1.49;PZ2次之,为45.11±0.01。PZ5红度(a*)最高,PZ3黄度(b*)最高且与其余样品间存在显著性差异(P<0.05)。产生这些差异的原因,可能与不同品种辣椒自身所富集的辣椒红素、β-胡萝卜素等含量不同有关[7]

图1 辣椒油样品的色差值
Fig.1 Chromatic aberration value of chili oil samples 注:小写字母不同表示组间差异显著(P<0.05)(下同)。

2.3 辣椒品种对辣椒油过氧化值的影响

油脂的过氧化值是评价油脂氧化程度的关键指标[24]。由图2可知,5种辣椒油的过氧化值依次为:PZ3P<0.05)。结合表1和图2可知,辣椒油中辣椒素类物质总量越高,其过氧化值越低,这与辣椒素类物质具有良好的抗氧化活性的构效特征相吻合[25]。由此可见,提高辣椒油中辣椒素类物质总量,能有效提高辣椒油的抗氧化能力,提高辣椒油的稳定性。

图2 辣椒油样品的过氧化值
Fig.2 Peroxide value of chili oil samples

2.4 辣椒品种对辣椒油VOCs的影响

2.4.1 辣椒油VOCs的定性分析

对照GC×IMS Library Search内置的NIST数据库和IMS数据库,根据辣椒油VOCs的保留指数、保留时间和迁移时间进行定性分析,结果如表2所示。部分VOCs出现了二聚体,它们具有相近的保留时间和不同的迁移时间,这是GC-MS所检测不出的[26]。5种样品中共检测出73种VOCs,94个信号峰,分别为醇类12种(二聚体3种)、醛类19种(二聚体5种)、酮类12种(二聚体3种)、羧酸类4种(二聚体3种)、酯类18种(二聚体5种)、杂环类6种(二聚体2种)、硫醚类2种。其中醛类(二乙醇缩乙醛、己醛、3-甲基丁醛、丙烯醛等)种类最多且含量高,阈值低[27],对辣椒油风味贡献大;其中,己醛主要来源于植物油的氧化[28],可为辣椒油提供青香和脂肪香;来源于异亮氨酸的Strecker降解[29]的3-甲基丁醛可为辣椒油提供果香。酯类物质主要来源于醇类和酸类化合物所发生的酯化反应[30],如乙酸乙酯、异戊酸己酯等酯类物质,可为辣椒油提供果香、青香。脂肪氧化的另一产物酮类[31]在辣椒油中含量也较高。此外,乙酸、丙酸等具有酸味的物质是辣椒油中的“异味”物质,降低“异味”物质的含量可使辣椒油风味得到提升[15]

表2 辣椒油样品已定性挥发性有机化合物
Table 2 Identified VOCs in chili oil samples

编号化合物名称CAS号保留指数保留时间/s迁移时间/ms风味描述醇类A1(Z)-3-己烯醇C9289611 359.8937.6791.232 22辛辣味、青香A2戊醇-MC714101 268.3770.8291.512 37辛辣味、杂醇香A3戊醇-DC714101 269771.8441.254 21辛辣味、杂醇香A43-甲基-1-丁醇C1235131 225.8705.8741.241 04辛辣味、果香A51-戊烯-3-醇C6162511 182.6637.8740.940 73辛辣味A61-丁醇-MC713631 166.3604.3821.378 02果香A71-丁醇-DC713631 166.8605.3971.183 09果香A82-甲基-1-丙醇C788311 114.6508.9791.172 55酒香A91-丙醇-MC712381 060.1430.4061.249 54酒精味、发酵味A101-丙醇-DC712381 061.1431.6661.111 03酒精味、发酵味A11乙醇C64175956325.2181.130 97酒香A12叔丁醇C75650930.3306.9521.319 85酒香A132-丁醇C789221 043.4409.621.147 76果香A141-己醇C1112731 373964.9361.326 26辛辣味、果香A154-甲基-1-戊醇C6268911 316.8854.0081.658 2坚果味醛类B1苯乙醛C1227811 762.92 248.6771.256 98玫瑰香、甜香B2壬醛C1241961 404.71 033.5381.476 48柑橘香、脂肪香B3(E)-2-辛烯醛C25488701 440.61 117.4161.334 25辛辣味、脂肪香B4(E)-2-庚烯醛-MC188295551 333.6885.7551.668 18脂肪香B5(E)-2-庚烯醛-DC188295551 334.3887.0861.260 05脂肪香B6辛醛C1241301 291.9809.3961.378 02脂肪香B72-己烯醛C5055771 261.9760.681.192 31青香、脂肪香B8(E)-2-己烯醛C67282631 234.8719.0681.181 77青香、脂肪香B9庚醛-MC1117171 201.6671.3671.699 41脂肪香、柑橘香B10庚醛-DC1117171 203.1673.3971.329 29脂肪香、柑橘香B11(E)-2-戊烯醛-MC15768701 155582.0541.359 58果香、辛辣味B12(E)-2-戊烯醛-DC15768701 155582.0541.105 37果香、辛辣味B13己醛-MC662511 109.1499.8451.562 43青香、脂肪香B14己醛-DC662511 107.3496.81.260 8青香、脂肪香B15戊醛-MC1106231 009369.9381.424 78发酵香、果香B16戊醛-DC1106231 008.5369.3091.173 99发酵香、果香B173-甲基丁醛C590863938.4312.621.398 55果香、奶酪味B18二乙醇缩乙醛C105577919.2299.3931.035 48甜香、脂肪香B19丁醛C123728901.1287.4261.281 02巧克力味B20丙烯醛C107028877.2272.3091.061 72果香B21丙醛C123386827.2243.3351.139 37辛辣味B222-甲基丙醛C78842839.7250.2631.279 97花香、辛辣B23苯甲醛C1005271 523.71 338.0921.458 36苦杏仁味、果香B24甲硫基丙醛C32684931 476.91 208.8251.093 56土豆味酮类C11-羟基-2-丙酮-MC1160961 316.8854.0531.230 5辛辣的甜味C21-羟基-2-丙酮-DC1160961 316.7853.8011.048 25辛辣的甜味C3乙偶姻-MC5138601 302.3827.6651.330 61苹果香、梨香C4乙偶姻-DC5138601 301.8826.651.059 27苹果香、梨香C51-辛烯-3-酮C43129961 314.6849.9931.284 51蘑菇味C6环己酮C1089411 316.2853.0381.151 47薄荷味C72-庚酮C1104301 191.1656.1431.264 75果香、辛辣味C84-甲基-3-戊烯-2-酮C1417971 133.2541.4571.111 96土豆味C92,3-戊二酮C6001461 085.3463.7891.226 46黄油味

续表2

编号化合物名称CAS号保留指数保留时间/s迁移时间/ms风味描述C101-戊烯-3-酮-MC16295891 046.5413.3991.311 45辛辣味C111-戊烯-3-酮-DC16295891 048415.2891.076 41辛辣味C122-戊酮C1078791 004.4364.91.347 13果香C132-丁酮C78933926.6304.4321.247 44果香C14丙酮C67641852.9257.8221.114 18苹果香、梨香C15丁二酮C431038976.1340.3351.183 44奶油味羧酸类D1丁酸-DC1079261 715.72 029.7331.160 14乳酪味D2丁酸-MC1079261 715.72 029.7331.363 14乳酪味D3丙酸C790941 640.61 724.4761.109 86酸味D4异丁酸-MC793121 631.51 690.7921.368 72酸味D5异丁酸-DC793121 631.51 690.7921.147 1酸味D6乙酸-MC641971 504.91 284.4841.154 55酸味D7乙酸-DC641971 504.91 284.4841.052 12酸味酯类E1乙酸丁酯C1238641 096.6479.5361.235 9果香E2丁酸乙酯C1055441 069.8443.0031.199 18果香E3乙酸异丁酯C1101901 035.5400.1721.232 75果香E4乙酸乙酯C141786907.9291.8351.338 74果香E5乙酸甲酯C79209844.1252.7831.200 23果香E6乙酸丙酯C1096041 010.2371.1981.504 53果香E7己酸异戊酯-MC21986101 460.51 166.7122.156 88果香E8己酸异戊酯-DC21986101 461.81 169.8911.523 46果香E9异戊酸己酯-MC100321521 411.61 049.1192.164 28青香、果香E10异戊酸己酯-DC100321521 411.61 049.1191.525 92青香、果香E11丁酸己酯-MC26396361 389.4999.8562.142 09青香、果香E12丁酸己酯-DC26396361 388996.6781.501 27青香、果香E13异丁酸己酯C23490771 355.1928.1681.455 78青香、果香E14乙酸叶醇酯C36817181 331.6881.9121.322 21青香、果香E15异戊酸异戊酯-MC6597011 305.5833.3682.042 3果香E16异戊酸异戊酯-DC6597011 306.5835.2451.453 22果香E172-甲基丁酸丁酯-MC157067371 202.1672.0011.898 26果香E182-甲基丁酸丁酯-DC157067371 202.7672.9391.379 35果香E19丙酸异戊酯C1056801 187.1647.6081.353 5果香E20异戊酸甲酯C5562411 039404.2771.528 84果香E21丙酸乙酯C105373997.6357.6291.477 92果香E22γ-丁内酯C964801 6701 838.1581.094 96奶油香E23二氢-2(3H)-呋喃酮C964801 715.22 027.6281.087 51奶油香杂环类F1糠醇-MC980001 695.91 944.1991.371 38面包味F2糠醇-DC980001 694.91 940.1321.148 88面包味F32-糠基硫醇C980221 430.21 092.2681.104 74烘焙香F4糠醛-MC980111 496.31 260.8271.335 67木香、甜香F5糠醛-DC980111 493.71 253.6541.091 7木香、甜香F62,5-二甲基吡嗪C1233201 338.2894.651.097 47坚果香、烘烤香F7甲基吡嗪C1090801 280.9791.1281.082 98坚果香F8噻吩C1100211 037.6402.6921.035 48大蒜味硫醚类G1二乙基二硫醚C1108161 201.6671.3671.135 67大蒜、洋葱香G2二甲基硫醚C75-18-3808.4233.2570.957 83洋葱香

注:“-M”表示单体,“-D”表示二聚体;编号为已定性的物质(下同)。

2.4.2 辣椒油VOCs的GC-IMS二维差异图

图3是以PZ1为参照,其余谱图扣除PZ1中的信号峰,得到的GC-IMS二维平面差异谱图。

图3 辣椒油样品的GC-IMS二维差异图
Fig.3 GC-IMS two-dimensional qualitative topographic diffence plots

整个二维谱图背景为蓝色,横坐标1.0处为反应离子(reactive ion peak,RIP)峰,RIP峰两侧的每一个点均代表一种VOC,白色表示VOC浓度较低,红色表示浓度较高,颜色越深表示浓度越大。由图3可知,矩形区域红色VOCs的含量相对于PZ1高,椭圆形区域蓝色VOCs的含量相对于PZ1低;PZ3、PZ4和PZ5的VOCs种类及含量与PZ1差异较大,PZ1和PZ2的VOCs种类和含量相近。

2.4.3 辣椒油VOCs的GC-IMS指纹图谱

指纹图谱(图4)可以直观地呈现5种辣椒油VOCs之间的差异,整个图背景为蓝色,Y轴为样品编号(每1行为1个样品的指纹图谱),X轴为VOCs的名称或编号,图中点的颜色深浅和点的面积表示VOCs含量的高低,颜色越深面积越大则含量越高,即白色点表示VOCs含量低,红色点表示含量高[1]

图4 辣椒油样品的指纹谱图
Fig.4 Fingerprint spectra of chili oil samples
注:21种未定性的VOCs以阿拉伯数字进行编号(下同)。

结合表2和图4可知,5种样品共有且含量较高的VOCs集中在A(A1~A4)区域,包括3-甲基丁醛、2-甲基丙醛、庚醛、2,3-戊二酮、环己酮、1-辛烯-3-酮、2,5-二甲基吡嗪、甲基吡嗪、(Z)-3-己烯醇、二氢-2(3H)-呋喃酮等,为辣椒油贡献辛辣味、脂肪香、青香、果香、花香等总体风味特征。B(B1~B2)区域,包括异戊酸己酯、异丁酸己酯、己酸异戊酯、丁酸己酯、异戊酸异戊酯、2-甲基丁酸丁酯、异戊酸甲酯、乙酸丙酯、乙酸叶醇酯、丙酸异戊酯、苯甲醛、丙烯醛、糠醛、壬醛、糠醇、1-己醇等,多数为酯类,在PZ3中含量高且与其余4种样品间存在显著性差异(P<0.05),所以PZ3果香味浓,风味最为特别。C区域的乙酸丁酯、丙酸乙酯、乙酸乙酯在PZ4中含量高且与其余样品间存在显著性差异(P<0.05),可能为PZ4与其余样品香气的差异所在。D区域的(E)-2-戊烯醛、乙酸甲酯、丁酸乙酯、辛醛在PZ5中含量高,可能为PZ5的主要差异VOCs。PZ1和PZ2的VOCs接近。

2.4.4 辣椒油样品VOCs的PCA

利用Flavor Spec®配备的LAV软件对5种辣椒油样品所有的VOCs进行PCA,结果如图5所示。由图5可知,第1主成分(58%)和第2主成分(19%)占方差的77%,能反映样品的总体特征[32]。PZ3与其余样品的距离非常远,可以根据第1主成分进行区分;其次,PZ5、PZ4分别与PZ1和PZ2的距离较远,可以根据第2主成分将PZ5、PZ4与PZ1和PZ2进行区分;PZ1和PZ2在第1与第2主成分上均接近,但各自成组;盲样(图中5个黑点)均准确地落在相应的组别中。因此,使用GC-IMS结合PCA,可将5种辣椒油样品的VOCs进行较好地判别和区分。

图5 辣椒油样品的PCA结果
Fig.5 PCA results of chili oil samples

2.4.5 辣椒油样品VOCs的PLS-DA结果

PLS-DA是一种基于偏最小二乘回归算法的有监督的分析方法,可实现复杂数据的可视化、判别分析和预测[33]。在PLS-DA中使用7倍交叉验证和200次随机重分类对辣椒油VOCs进行建模、因子载荷分析和变量重要性投影(variable importance in projection,VIP)计算。利用SIMCA 14.1软件对5种辣椒油所有VOCs进行PLS-DA,结果如图6所示。

a-PLS-DA置换图;b-因子载荷图

图6 辣椒油样品VOCs的PLS-DA置换图和因子载荷图
Fig.6 Permutation test plots and scores plot of PLS-DA based on the VOCs of chili oil samples

由图6-a可知,RX2=0.977,RY2=0.992,Q2=0.972,说明该模型较可靠,能对5种辣椒油样品的风味进行较好的预测。根据PLS-DA因子载荷图对5种辣椒油重要的香气成分进行分析,如图6-b所示,PZ3的香气物质主要为酯类,包括E10(异戊酸己酯-D)、E13(异丁酸己酯)、E7(己酸异戊酯)、E12(丁酸己酯-D)等,因此PZ3果香味浓郁。PZ4的主要香气成分有E1(乙酸丁酯)、A4(3-甲基-1-丁醇)、E21(丙酸乙酯)等;PZ5的主要香气成分有B11[(E)-2-戊烯醛-M]、F6(2,5-二甲基吡嗪)、E5(乙酸甲酯)等;C11(1-戊烯-3-酮-D)、C7(2-庚酮)等是PZ1的主要香气成分,PZ1和PZ2风味相似。

VIP可用于筛选对香气轮廓有重要影响的关键差异标志物[33]。如图7-a所示,26种VOCs的VIP>1.0(共32个峰),分别为E4(乙酸乙酯)、E10(异戊酸己酯)、B18(二乙醇缩乙醛)、D7(乙酸-D)、B14(己醛-D)、G2(二甲基硫醚)、B17(3-甲基丁醛)、E23(二氢-2(3H)-呋喃酮)、C2(1-羟基-2-丙酮)、C14(丙酮)、1、D2(丁酸-M)、B20(丙烯醛)、A4(3-甲基-1-丁醇)、D4(异丁酸-M)、E21(丙酸乙酯)、F1(糠醇-M)、B21(丙醛)、C11(1-戊烯-3-酮)、B19(丁醛)、A11(乙醇)、E7(己酸异戊酯-M)、C15(丁二酮)、E11(丁酸己酯-M)、B16(戊醛-D)、F4(糠醛-M),VIP值越大,则差异越显著。由图7-b可知,E9(异戊酸己酯-M)、E10(异戊酸己酯-D)、E11(丁酸己酯-M)、F1(糠醇-M)、E7(己酸异戊酯-M)、D7(乙酸-D)、B17(3-甲基丁醛)、B16(戊醛-D)等物质对第1主成分贡献较大;E21(丙酸乙酯)、A4(3-甲基-1-丁醇)、E4(乙酸乙酯)、B19(丁醛)、B18(二乙醇缩乙醛)对第2主成分贡献较大。

a-VIP值图;b-载荷图;c-聚类热图

图7 辣椒油重要风味物质(VIP>1.0)的PLS-DA的VIP值图、载荷图和聚类热图
Fig.7 VIP scores, scores plot, and cluster heat map of PLS-DA based on the VOCs of chili oil samples

为了直观区分特征差异标志物(VIP>1.0)在5种辣椒油中的变化,根据筛选出的26种特征差异标志物的峰强度绘制聚类热图(图7-c)。由图7-c可知,PZ3最先分组,风味最为特别,E10(异戊酸己酯)、E11(丁酸己酯-M)、E7(己酸异戊酯-M)、B20(丙烯醛)、B16(戊醛-D)等物质在PZ3中含量高;PZ5中 C11(1-戊烯-3-酮)、B21(丙醛)、C15(丁二酮)含量高而再分组;PZ4中E4(乙酸乙酯)、A4(3-甲基-1-丁醇)、E21(丙酸乙酯)含量高;PZ1和PZ2聚类而风味相似。该结果与GC-IMS二维差异图、指纹图谱和PCA结果一致。

3 结论

本研究采用HPLC和GC-IMS等检测技术,对5种辣椒(子弹头、龙脊椒、涮涮辣、七星椒、干线椒)制备的辣椒油进行了辣椒素类物质含量、色差、过氧化值、VOCs的检测和分析,辣椒的品种对制成的辣椒油中辣椒素类物质含量、色差、过氧化值、VOCs均有影响。结果表明,5种辣椒油辣度有显著性差异(P<0.05),其中PZ3的辣椒素类物质总量、辣度显著高于其他样品,PZ4其次,PZ1最低;PZ3亮度L*和黄度b*最高;5种辣椒油的过氧化值呈现出随着辣椒油中辣椒素类物质总量的增加而降低的趋势,PZ3过氧化值最低,辣椒素类物质总量最高;经GC-IMS共分离鉴定出73种VOCs,包括醛类19种、酯类18种、酮类12种、醇类12种、杂环类6种、羧酸类4种和硫醚类2种;其中,GC-IMS二维差异图和指纹图谱显示,PZ3样品VOCs种类多含量高,风味最为独特,PZ1和PZ2风味接近;通过PLS-DA进一步分析发现,PZ3与其余样品的特征差异香气物质主要为酯类,包括异戊酸己酯-D、异丁酸己酯、己酸异戊酯、丁酸己酯-D等;经变量重要性投影分析筛选出5种辣椒油的26种关键差异标志物(VIP>1.0),包括E4(乙酸乙酯)、E10(异戊酸己酯)、B18(二乙醇缩乙醛)、D7(乙酸-D)、B14(己醛-D)、G2(二甲基硫醚)、B17(3-甲基丁醛)等,其中E9(异戊酸己酯-M)、E10(异戊酸己酯-D)、E11(丁酸己酯-M)、F1(糠醇-M)、E7(己酸异戊酯-M)、D7(乙酸-D)、B17(3-甲基丁醛)、B16(戊醛-D)等物质,对第1主成分贡献较大;E21(丙酸乙酯)、A4(3-甲基-1-丁醇)、E4(乙酸乙酯)、B19(丁醛)、B18(二乙醇缩乙醛)对第2主成分贡献较大。PCA结果、聚类热图分析结果与GC-IMS二维差异图和指纹图谱结果一致。研究结果将为合理选用辣椒品种,工业化制备辣椒油提供一定的理论支撑。

参考文献

[1] YANG L Z, LIU J, WANG X Y, et al.Characterization of volatile component changes in jujube fruits during cold storage by using headspace-gas chromatography-ion mobility spectrometry[J].Molecules, 2019, 24(21):3904.

[2] GE X Z, LIANG Q H, LONG Y, et al.Assessment of fresh Alpinia galanga (A.galanga) drying techniques for the chemical composition of essential oil and its antioxidant and biological activity[J].Food Chemistry, 2022, 392:133314.

[3] LI W D, LI J J, QIN Z, et al.Insights into the composition and antibacterial activity of Amomum tsao-ko essential oils from different regions based on GC-MS and GC-IMS[J].Foods, 2022, 11(10):1402.

[4] 倪瑞洁, 詹萍, 田洪磊.基于GC-IMS结合多元统计方法分析炸制时间对花椒调味油挥发性物质的影响[J].食品科学, 2022, 43(6):279-286.

NI R J, ZHAN P, TIAN H L.Effects of frying time on volatile flavor compounds in fried pepper (Zanthoxylum bungeanum) oil as analyzed by gas chromatography-ion mobility spectrometry and multivariate statistical analysis[J].Food Science, 2022, 43(6):279-286.

[5] 杨芳, 范成梦, 贾洪锋, 等.基于气相色谱-离子迁移谱对不同产地羊肚菌风味化合物的分析[J].食品与发酵工业, 2021, 47(10):207-213.

YANG F, FAN C M, JIA H F, et al.Analysis of flavor compounds of Morchella spp.from different habitats based on gas chromatography-ion mobility spectrum[J].Food and Fermentation Industries, 2021, 47(10):207-213.

[6] WANG S Q, CHEN H T, SUN B G.Recent progress in food flavor analysis using gas chromatography-ion mobility spectrometry (GC-IMS)[J].Food Chemistry, 2020, 315:126158.

[7] DE S MENDES N, BRANCO DE ANDRADE GONÇALVES É C.The role of bioactive components found in peppers[J].Trends in Food Science &Technology, 2020, 99:229-243.

[8] SINISGALLI C, VEZZA T, DIEZ-ECHAVE P, et al.The beneficial effects of red sun-dried Capsicum annuum L.cv senise extract with antioxidant properties in experimental obesity are associated with modulation of the intestinal microbiota[J].Molecular Nutrition &Food Research, 2021, 65(3):e2000812.

[9] IC E.Quantitative viscosity determination in irradiated major spices (black pepper, cardamom, cinnamon, ginger, and turmeric) by using a vibro viscometer for long-term storage[J].Food Control, 2022, 133:108679.

[10] 杜坚, 暴天超, 李垠青, 等.大蒜油和辣椒油对童鸽生长性能、免疫器官指数、血清生化与抗氧化指标的影响[J].中国家禽, 2021, 43(4):54-59.

DU J, BAO T C, LI Y Q, et al.Effects of garlic oil and chili oil on growth performance, immune organ index, serum biochemical and antioxidant indexes of squeaker[J].China Poultry, 2021, 43(4):54-59.

[11] 萨仁高娃, 胡文忠, 姜爱丽.辣椒营养保健功能及辣椒食品的研究进展[J].食品工业科技, 2012, 33(15):371-375.

SA R G W, HU W Z, JIANG A L.Research process in pharmacological function and products of chilli[J].Science and Technology of Food Industry, 2012, 33(15):371-375.

[12] DE S MENDES N, SANTOS M C P, SANTOS M C B, et al.Characterization of pepper (Capsicum baccatum)—A potential functional ingredient[J].LWT, 2019, 112:108209.

[13] 王兴波, 饶雷, 王永涛, 等.9个品种干辣椒的品质分析及评价[J].食品工业科技, 2022, 43(18):300-310.

WANG X B, RAO L, WANG Y T, et al.Quality analysis and evaluation of nine varieties of dried peppers[J].Science and Technology of Food Industry, 2022, 43(18):300-310.

[14] 杜勃峰, 李达, 肖仕芸, 等.基于HS-SPME-GC-MS及主成分分析综合评价贵州典型辣椒品种香气品质[J].食品研究与开发, 2019, 40(7):149-155.

DU B F, LI D, XIAO S Y, et al.Comprehensive evaluation of aroma quality of different pepper based on HS-SPME-GC-MS and principal component analysis[J].Food Research and Development, 2019, 40(7):149-155.

[15] 杨慧, 黄绿红, 张帆, 等.基于感官和气相色谱-离子迁移谱分析油温对辣椒油风味的影响[J].中国食品学报, 2021, 21(9):328-335.

YANG H, HUANG L H, ZHANG F, et al.Effect of oil temperature on the flavor of chilli oil by GC-IMS and sensory evaluation[J].Journal of Chinese Institute of Food Science and Technology, 2021, 21(9):328-335.

[16] 张洪新, 刘辉, 陈光静, 等.不同品种辣椒油制品的品质研究[J].食品与发酵工业, 2019, 45(5):190-198.

ZHANG H X, LIU H, CHEN G J, et al.Study on the quality of chili oil products from different chili varieties[J].Food and Fermentation Industries, 2019, 45(5):190-198.

[17] 何小龙, 周晓燕, 李辉, 等.辣椒油制作过程中的品质变化研究[J].食品研究与开发, 2014, 35(11):91-93.

HE X L, ZHOU X Y, LI H, et al.Study on different varieties pepper effect on quality of chili oil[J].Food Research and Development, 2014, 35(11):91-93.

[18] 石自彬, 周世中.辣椒品种对辣椒油品质影响研究[J].粮食与油脂, 2022, 35(2):69-73.

SHI Z B, ZHOU S Z.Study on the effect of chili varieties on the quality of chili oil[J].Cereals &Oils, 2022, 35(2):69-73.

[19] 郭家刚, 杨松, 丁思年, 等.基于气相离子迁移谱的不同产地生姜挥发性有机物指纹图谱分析[J].食品科学, 2021, 42(24):236-241.

GUO J G, YANG S, DING S N, et al.Fingerprint analysis of volatile organic compounds in ginger rhizomes from different geographical origins by gas chromatography-ion mobility spectrometry[J].Food Science, 2021, 42(24):236-241.

[20] 国家质量监督检验检疫总局, 中国国家标准化管理委员会.辣椒及辣椒制品中辣椒素类物质测定及辣度表示方法:GB/T 21266—2007[S].北京:中国标准出版社, 2008.

General Administration of Quality Supervision, Inspection and Quarantine of the People’s Republic of China, Standardization Administration of China.Determination of total capsaicinoid content and representation of pungency degree in capsicum and its products:GB/T 21266—2007[S].Beijing:Standards Press of China,2008.

[21] 中华人民共和国国家卫生和计划生育委员会. 食品安全国家标准 食品中过氧化值的测定:GB 5009.227—2016[S].北京:中国标准出版社, 2016.

National health and family planning commission of the People’s Republic of China.National food safety standard determination of peroxide value in foods:GB 5009.227—2016[S].Beijing:Standards Press of China,2016.

[22] 王鑫, 王兴国, 金青哲, 等.HPLC法测定红辣椒油树脂中辣椒素的含量[J].食品科学, 2005, 26(11):207-209.

WANG X, WANG X G, JIN Q Z, et al.Determination of capsaicinoids in paprika oleoresin by HPLC[J].Food Science, 2005, 26(11):207-209.

[23] 杨莉, 张淼, 贾洪锋, 等.基于高效液相色谱法的麻辣火锅底料辣度分级模型建立及其在熬煮过程中的变化[J].食品科学, 2021, 42(4):233-239.

YANG L, ZHANG M, JIA H F, et al.Modelling for pungency grading of spicy hot pot seasonings based on capsaicinoid content determined by HPLC and analysis of its changes during boiling[J].Food Science, 2021, 42(4):233-239.

[24] 易志, 吴雪辉, 沈冰, 等.温度及光照对亚麻籽油贮藏稳定性影响研究[J].粮食与油脂, 2016, 29(6):17-21.

YI Z, WU X H, SHEN B, et al.Effects of temperature and light on the storage stability of flaxseed oil[J].Cereals &Oils, 2016, 29(6):17-21.

[25] 王梦, 赵佩霞, 张鹏, 等.辣椒碱单体、二氢辣椒碱单体和降二氢辣椒碱单体抗氧化效果研究[J].化学与生物工程, 2019, 36(2):8-11.

WANG M, ZHAO P X, ZHANG P, et al.Antioxidant effect of capsaicin monomer, dihydrocapsaicin monomer, and nordihydrocapsaicin monomer[J].Chemistry &Bioengineering, 2019, 36(2):8-11.

[26] LI M Q, YANG R W, ZHANG H, et al.Development of a flavor fingerprint by HS-GC-IMS with PCA for volatile compounds of Tricholoma matsutake Singer[J].Food Chemistry, 2019, 290:32-39.

[27] ADEBOWALE K O, LAWAL O S.Microstructure, physicochemical properties and retrogradation behaviour of Mucuna bean (Mucuna pruriens) starch on heat moisture treatments[J].Food Hydrocolloids, 2003, 17(3):265-272.

[28] XU L R, MEI X, CHANG J R, et al.Comparative characterization of key odorants of French fries and oils at the break-in, optimum, and degrading frying stages[J].Food Chemistry, 2022, 368:130581.

[29] COUNET C, CALLEMIEN D, OUWERX C, et al.Use of gas chromatography-olfactometry to identify key odorant compounds in dark chocolate comparison of samples before and after conching[J].Journal of Agricultural and Food Chemistry, 2002, 50(8):2385-2391.

[30] MATEO J, ZUMALACRREGUI J M.Volatile compounds in chorizo and their changes during ripening[J].Meat Science, 1996, 44(4):255-273.

[31] 陈海涛, 张宁, 徐晓兰, 等.SPME和SDE-GC-MS分析贾永信腊羊肉挥发性风味成分[J].食品科学, 2013, 34(14):187-191.

CHEN H T, ZHANG N, XU X L, et al.Comparison of simultaneous distillation extraction (SDE) and solid-phase microextraction (SPME) for the analysis of volatile flavor compounds in jiayongxin preserved mutton by GC-MS[J].Food Science, 2013, 34(14):187-191.

[32] WU Z B, CHEN L Z, WU L M, et al.Classification of Chinese honeys according to their floral origins using elemental and stable isotopic compositions[J].Journal of Agricultural and Food Chemistry, 2015, 63(22):5388-5394.

[33] EVA.GC-MS analysis and PLS-DA validation of the trimethyl silyl-derivatization techniques[J].American Journal of Applied Sciences, 2012, 9(7):1124-1136.

Effects of chili varieties on physicochemical properties and flavor compounds of chili oil based on GC-IMS combined with multivariate statistical methods

YANG Fang, YUAN Haibin, JIA Hongfeng*, DENG Fenglin, WANG Zhenni

(College of Food Science and Technology, Sichuan Tourism University, Chengdu 610100, China)

ABSTRACT In this research, the content of capsaicinoids, chromatic aberration value, peroxide value, and volatile organic compounds (VOCs) in chili oil prepared from 5 different varieties of chili (Capsicum frutescens var. breviconoideum Haz, Longji chili, Capsicum frutescens L. var. shuanlaense, Qixing chili, and Trunk chili) were investigated. Combined with multivariate statistical analysis, the markers of chili oil samples were explored. Results showed that the capsaicin content, hotness, brightness, and yellowness of PZ3 were the highest, and the peroxide value was the lowest, which was a significant difference from the rest of the samples (P<0.05). A total of 73 VOCs were separated and identified by gas chromatography-ion mobility spectrometry (GC-IMS, including aldehydes, esters, ketones, alcohols, heterocycles, carboxylic acids, and thioethers, corresponding to 19, 18, 12, 12, 6, 4, and 2 respectively. The results of the principal component analysis were consistent with GC-IMS two-dimensional qualitative topographic difference plots and fingerprint spectra, that, PZ3 had the most unique flavor, followed by PZ5 and PZ4, and the flavors of PZ1 and PZ2 were very close. The results of the partial least squares-discriminant analysis showed that the characteristic aroma substances of PZ3 were mainly esters, including hexyl 2-methylbutyrate-D, hexyl isobutyrate, isoamyl hexanoate, hexyl butyrate-D, etc. Through variable importance in projection analysis, 26 marker differences of 5 kinds of chili oil were screened out, including ethyl acetate, hexyl 2-methylbutyrate, diethyl acetal, acetic acid-D, hexanal-D, dimethyl sulfide, 3-methylbutanal, dihydro-2(3H)-furanone, 1-hydroxy-2-acetone-D, acetone, butanoic acid-M, et al. Among the above-mentioned compounds, hexyl 2-methylbutyrate, hexyl butyrate-M, furfuryl alcohol-M, isoamyl hexanoate-M, acetic acid-D, 3-methylbutanal, and pentanal-D had larger load values in the first main component, while ethyl propanoate, 3-methyl-1-butanol, ethyl acetate, butanal, and diethyl acetal had larger loads values in the second main component.

Key words chili oil; gas chromatography-ion mobility spectrometry; capsaicin; volatile flavor components; partial least squares-discriminant analysis; variable importance in projection

DOI:10.13995/j.cnki.11-1802/ts.033362

引用格式:杨芳,袁海彬,贾洪锋,等.基于气相色谱-离子迁移谱结合多元统计方法分析辣椒品种对辣椒油理化性质和风味物质的影响[J].食品与发酵工业,2023,49(19):319-328.YANG Fang, YUAN Haibin, JIA Hongfeng, et al.Effects of chili varieties on physicochemical properties and flavor compounds of chili oil based on GC-IMS combined with multivariate statistical methods[J].Food and Fermentation Industries,2023,49(19):319-328.

第一作者:硕士,副教授(贾洪锋教授为通信作者,E-mail:jiahongfeng_cq@163.com)

基金项目:四川省自然科学基金项目(2022NSFSC1750);川菜工业化四川省高等学校工程研究中心项目(GCZX22-23);四川旅游学院餐饮食品感官品质智能评价科研创新团队项目(19SCTUTY04);四川旅游学院大学生创新创业训练计划项目(202111552051)

收稿日期:2022-08-18,改回日期:2022-10-08