花椒是（Zanthoxylum bungeanum Maxim.）芸香科植物青椒或花椒的干燥成熟果皮[1-2]，全世界约有250 种，我国约有45 种和14 个变种[3]，按照果实颜色分为红花椒和青花椒。花椒具有特殊的辛、麻、香味[4]，是一种典型的药食同源的植物[5]，位列卫健委公布的既是食品又是药品的中药名单，有治疗肝损伤[6]、抑菌杀虫[7-10]、缓解疼痛、消炎[11-13]、抗氧化、抗癌[14-19]等作用。花椒的活性化学物质丰富，富含挥发油[20]、酰胺类物质[21]、生物碱[22]、香豆素[23]、木脂素[24]和脂肪酸[25]等。目前，脂肪酸分析技术主要有气相色谱法、高效液相色谱法、红外光谱法、质谱法以及薄层层析法[26]。
目前国内外对于花椒的研究主要集中在栽培技术[27]、采后干制、加工技术[28]以及花椒挥发油的含量、化学组成等方面。候丽秀[28]对不同地区青花椒籽和红花椒籽进行油脂萃取，再利用气相色谱-质谱联用仪（gas chromatography-mass spectrometer，GC-MS）对籽油进行脂肪酸测定。曹蕊等[3]采用索氏抽提法从花椒中提取花椒油，用GC-MS 进行分析并对两种花椒油的脂肪酸组成和相对含量进行比较。张仁凤等[29]采用GC-MS 对6 个不同产地红花椒籽油的脂肪酸组成进行了分析。已有的研究分析了不同产地花椒的脂肪酸种类，表明花椒样品的脂肪酸组成与产地存在一定的相关性。
在我国花椒主要产自四川、陕西等地，产地不同造成市场价格差异。本文选取陕西、山东、四川、甘肃4个省份共32 份花椒，并对脂肪酸成分进行提取，结合气相色谱法对其进行分析。利用37 种脂肪酸甲酯混合标准品确定脂肪酸成分并采用外标法确定各成分的含量。进一步采用方差分析和多重比较分析不同省份花椒脂肪酸组分的显著性差异，进一步探究花椒的脂肪酸含量及其组成差异与产地的相关性，旨在为花椒在食品、医药等领域的进一步应用提供参考。
1 材料与方法
1.1 材料与试剂
花椒样品（干制品）：来自陕西韩城8 份，山东泰安、枣庄各4 份，四川汉源2 份、茂汶6 份，甘肃天水3份、陇南4 份、临夏1 份，4 个省份共计32 份样品。所有样品均自行采摘，在未检测分析前均在-20 ℃避光密封保存。
盐酸（分析纯）、焦性没食子酸（C6H6O3、分析纯）：西陇化工股份有限公司；石油醚（沸程30～60 ℃、色谱纯）：无锡市展望化工试剂有限公司；正己烷（C6H14、色谱纯）、氢氧化钠（分析纯）、无水硫酸钠（分析纯）、氯化钠（分析纯）：国药集团化学试剂有限公司；乙醇（纯度95%，色谱纯）：宿州众望健康科技有限公司；甲醇（色谱纯）：安徽高地高纯溶剂有限公司；正庚烷（色谱纯）：上海跃胜贸易有限公司；三氟化硼甲醇溶液（纯度15%，色谱纯）、37 种混合脂肪酸甲酯标准品：上海安谱实验科技股份有限公司。
1.2 仪器与设备
气相色谱仪（Agilent 7890B）：美国安捷伦科技公司；研磨机：浙江苏泊尔有限公司；PL403 电子天平（精确度0.1 mg）：梅特勒-托利多仪器上海有限公司；恒温水浴锅（40～100 ℃，控温±1 ℃）：上海申贤恒温设备厂；RV-10digitalV 数显型立式旋转蒸发仪：德国艾卡集团。
1.3 试验方法
1.3.1 标准溶液的配制
取出适量脂肪酸甲酯混合标准品移至10mL 容量瓶中，用正庚烷稀释定容成50、100、500、1 000、2 500 mg/L和5 000 mg/L 6 个浓度，贮存于-20 ℃冰箱，过滤膜后装瓶待测。
1.3.2 脂肪酸的提取
在GB 5009.168—2016《食品安全国家标准食品中脂肪酸的测定》[30]的基础上改进了溶剂，以正己烷代替乙醚[31]，对花椒的脂肪酸成分进行提取，减少了测定样品中杂质的产生且在一定程度上提高了试验的安全性和稳定性。
先将花椒样品使用研磨机粉碎，称取均匀试样1.0 g移入到250 mL 平底烧瓶中，加入100 mg 焦性没食子酸，再加入2 mL 95%乙醇和4 mL 水，混匀；加入盐酸溶液10 mL，混匀。将烧瓶放入75 ℃水浴中水解40 min。每隔10 min 振荡一下烧瓶，使黏附在烧瓶壁上的颗粒物混入溶液中。
水解完成后，取出烧瓶冷却至室温。水解后的试样，加入10 mL 95%乙醇，混匀。将烧瓶中的水解液转移到分液漏斗中，用50 mL 正己烷石油醚混合液冲洗烧瓶和塞子，冲洗液并入分液漏斗中，加盖。振摇5 min，静置10 min。将醚层提取液收集到250 mL 烧瓶中。按照以上步骤重复提取水解液3 次，最后用正己烷石油醚混合液冲洗分液漏斗，并收集250 mL 烧瓶中。
旋转蒸发浓缩至干，残留物为脂肪提取物。在脂肪提取物中加入2%氢氧化钠甲醇溶液8 mL，连接回流冷凝器，（80±1）℃水浴上回流，直至油滴消失。从回流冷凝器上端加入7 mL 15%三氟化硼甲醇溶液，在（80±1）℃水浴中继续回流2 min。用少量水冲洗回流冷凝器。停止加热，从水浴上取下烧瓶，迅速冷却至室温。准确加入20 mL 正庚烷，振摇2 min，再加入饱和氯化钠水溶液，静置分层。吸取上层正庚烷提取溶液5 mL至25 mL 试管中，加入大约4 g 无水硫酸钠，振摇1 min静置5 min，吸取上层溶液到进样瓶中待测定。
1.3.3 色谱柱及程序升温条件
HP-88 毛细管柱（100 m×0.25 mm×0.20 μm）；进样口温度：270 ℃；FID 检测器温度：280 ℃；程序升温：起始温度100 ℃，保持13 min，以10 ℃/min 升至180 ℃，保持6 min，再以1 ℃/min 升至200 ℃，保持20 min，再以4 ℃/min 升至230 ℃，保持10.5 min；载气为高纯度氮气（N2 99.999%）；分流比20∶1；进样体积1.0 μL。
1.4 数据分析
采用Excel 2010 进行花椒脂肪酸数据处理分析并制作37 种混合脂肪酸甲酯溶液的标准曲线；脂肪酸含量柱形图采用GraphPad Prism 7.00 软件完成；单因素方差分析采用GraphPad Prism 7.00 和IBM SPSS Statistics 22.0 软件完成；层次聚类分析（hierarchical cluster analysis，HCA）和相关性分析采用Matlab 软件完成；热图采用HemI 软件完成。
2 结果与分析
2.1 脂肪酸的GC 检测
按照1.3.3 中的方法对6 个浓度梯度的37 种脂肪酸甲酯混合标准品进行GC 检测，以5 000 mg/L 浓度的色谱图为代表，如图1a 所示。按照1.3.2 中的方法提取32 个不同产地的花椒脂肪酸成分并进行GC 检测，每个样品重复3 次，随机以茂汶、陇南、泰安、韩城4 处的色谱图分别作为四川、甘肃、山东、陕西4 省花椒样品色谱图，如图1b～图1e 所示。
图1 37 种混合脂肪酸甲酯标准品和花椒中脂肪酸甲酯组分的气相色谱图
Fig.1 Gas chromatography profiles of mixed methyl fatty acid for 37 standard and Zanthoxylum bungeanum
a.37 种混合脂肪酸甲酯标准品气相色谱图；b～e.不同省份花椒中脂肪酸甲酯组分气相色谱图。1.丁酸甲酯；2.己酸甲酯；3.辛酸甲酯；4.癸酸甲酯；5.十一碳酸甲酯；6.十二碳酸甲酯；7.十三碳酸甲酯；8.十四碳酸甲酯；9.顺-9-十四碳一烯酸甲酯；10.十五碳酸甲酯；11.顺-10-十五碳一烯酸甲酯；12.十六碳酸甲酯；13.顺-9-十六碳一烯酸甲酯；14.十七碳酸甲酯；15.顺-10-十七碳一烯酸甲酯；16.十八碳酸甲酯；17.反-9-十八碳一烯酸甲酯；18.顺-9-十八碳一烯酸甲酯；19.反，反-9，12-十八碳二烯酸甲酯；20.顺，顺-9，12-十八碳二烯酸甲酯；21.二十碳酸甲酯；22.顺-顺-顺-6，9，12-十八碳三烯酸甲酯；23.顺-11-二十碳一烯酸甲酯；24.顺，顺，顺-9，12，15-十八碳三烯酸甲酯；25.二十一碳酸甲酯；26.顺，顺-11，14-二十碳二烯酸甲酯；27.二十二碳酸甲酯；28.顺，顺，顺-8，11，14-二十碳三烯酸甲酯；29.顺-13-二十二碳一烯酸甲酯；30.顺11，14，17-二十碳三烯酸甲酯；31.顺-5，8，11，14-二十碳四烯酸甲酯；32.二十三碳酸甲酯；33.顺13，16-二十二碳二烯酸甲酯；34.二十四碳酸甲酯；35.顺-5，8，11，14，17-二十碳五烯酸甲酯；36.顺-15-二十四碳一烯酸甲酯；37.顺-4，7，10，13，16，19-二十二碳六烯酸甲酯。
从色谱图中可以看出，4 省花椒样品所含主要脂肪酸组分出峰时间大致相同，但峰高有所不同。
2.2 脂肪酸的定性定量结果
不同省份花椒脂肪酸含量见表1。
表1 不同省份花椒脂肪酸含量
Table 1 Fatty acid content of Zanthoxylum bungeanum from different provinces
序号脂肪酸名称保留时间/min标准曲线R2含量/（mg/g）陕西山东四川甘肃1己酸17.995 Y=0.007 2x+0.782 0 0.995 1 2.73±2.261.72±2.181.51±2.942.37±3.81 2十六碳酸35.556 Y=0.045 7x+3.792 8 0.996 3 25.61±6.44 31.46±15.85 27.91±17.91 16.09±7.90 3顺-9-十六碳一烯酸36.442 Y=0.022 4x+2.093 9 0.995 7 20.98±7.05 26.00±29.50 17.16±24.61 0.38±0.71 4十七碳酸38.198 Y=0.023 1x+2.251 5 0.992 5 0.02±0.0701.90±2.400 5顺-10-十七碳一烯酸39.174 Y=0.022 8x+2.184 2 0.995 7 0.06±0.110.18±0.380.35±0.640 6十八碳酸40.303 Y=0.024 4x+2.017 6 0.995 4 2.20±1.014.28±3.631.57±1.340.86±1.18 7顺-9-十八碳一烯酸42.222 Y=0.026 2x+0.961 7 0.997 7 52.93±15.92 77.25±49.27 33.84±10.93 35.21±24.06 8顺-顺-9，12-十八碳二烯酸45.581 Y=0.036 7x-4.262 9 0.987 8 38.78±10.52 69.95±45.40 35.46±16.91 25.31±9.54 9二十碳酸45.728 Y=0.025 8x+0.982 7 0.996 4 0.19±0.130.61±0.230.44±0.480.21±0.33 10顺-顺-顺-6，9，12-十八碳三烯酸47.089 Y=0.022 5x+2.093 4 0.994 5 4.59±1.154.92±2.042.09±1.390.81±0.43 11顺-11-二十碳一烯酸47.637 Y=0.023 1x+2.127 1 0.995 5 53.00±17.90 86.77±51.78 52.74±19.65 34.91±8.09 12顺-顺-顺-9，12，15-十八碳三烯酸49.291 Y=0.023 4x+2.197 0 0.994 600.03±0.0800 13二十一碳酸51.106 Y=0.026 4x+1.194 5 0.996 6 0.08±0.070.69±0.310.56±0.310.52±0.40 14二十二碳酸53.875 Y=0.028 2x+0.065 4 0.998 900.45±0.500.59±0.720.17±0.20 15二十三碳酸58.879 Y=0.028 9x+0.484 4 0.998 3 0.78±0.612.25±0.761.61±0.770.91±0.59 16顺-5，8，11，14，17-二十碳五烯酸65.236 Y=0.030 0x-0.255 8 0.999 5 0.24±0.000.39±0.220.42±0.210.27±0.07
采用外标法，通过6 个浓度梯度标准溶液的浓度C与标准品峰面积A 之间的关系，绘制标准曲线A=kC＋B[32-33]。37 种脂肪酸在50～5 000 mg/L 线性范围内，相关系数R2 在0.99～0.999 9 之间，线性关系良好。通过GC检测与定性分析，共鉴定出16 种脂肪酸，其中饱和脂肪酸和不饱和脂肪酸均有8 种。4 省份共有的脂肪酸有12 种，其中饱和脂肪酸和不饱和脂肪酸均有6 种。
四川和山东花椒样品均检测出15 种脂肪酸，在4 省中种类最多。四川花椒样品饱和脂肪酸有8 种，不饱和脂肪酸有7 种，顺-顺-顺-9，12，15-十八碳三烯酸未检出；山东花椒样品饱和脂肪酸有7 种，不饱和脂肪酸有8 种，十七碳酸未检出；陕西花椒样品共检测出14 种脂肪酸，其中饱和脂肪酸和不饱和脂肪酸均有7 种，顺-顺-顺-9，12，15-十八碳三烯酸和二十二碳酸未检出；甘肃花椒样品共检测出13 种脂肪酸，在4 省中种类最少，其中饱和脂肪酸有7 种，不饱和脂肪酸有6 种，十七碳酸、顺-10-十七碳一烯酸和顺-顺-顺-9，12，15-十八碳三烯酸未检出。顺-11-二十碳一烯酸（花生烯酸）、顺-9-十八碳一烯酸（油酸）、顺-顺-9，12-十八碳二烯酸（亚油酸）、十六碳酸（棕榈酸）和顺-9-十六碳一烯酸（棕榈油酸）5 种脂肪酸的含量明显高于其他脂肪酸，而甘肃省样品顺-9-十六碳一烯酸（棕榈油酸）的含量明显低于其他3 省份。陕西、山东和四川花椒样品所含脂肪酸含量最高的均为顺-11-二十碳一烯酸（花生烯酸），而甘肃花椒样品含量最高的为顺-9-十八碳一烯酸（油酸）。
在表1 中共有8 组脂肪酸含量数据的标准差略大于均值，说明这些数据较离散并存在异常值。不同产地虽在行政区域上隶属于同一省份，但地理位置有时相差甚远，甚至毗邻其他省份，推测可能是造成数据存在异常值的原因，但大部分的结果在统计学上是合理的。
2.3 不同省份花椒脂肪酸的显著性差异分析
在统计学中，差异显著性检验是用于检测试验组之间是否有差异及差异是否显著的方法，方差分析用于两个及两个以上样本均数差别的显著性检验。将4 个省份32 份花椒样品每份做3 个平行，共96 个样本通过GC 检测与定性定量分析，利用脂肪酸的种类和含量数据进行单因素方差分析，以产地为因子，以脂肪酸含量为因变量，结果见图2。多重比较检验见表2。
图2 4 个省份花椒脂肪酸含量
Fig.2 Fatty acid content of Zanthoxylum bungeanum from four provinces
A～P 依次对应己酸、十六碳酸、顺-9-十六碳一烯酸、十七碳酸、顺-10-十七碳一烯酸、十八碳酸、顺-9-十八碳一烯酸、顺-顺-9，12-十八碳二烯酸、二十碳酸、顺-顺-顺-6，9，12-十八碳三烯酸、顺-11-二十碳一烯酸、顺-顺-顺-9，12，15-十八碳三烯酸、二十一碳酸、二十二碳酸、二十三碳酸、顺-5，8，11，14，17-二十碳五烯酸16 种脂肪酸；不同小写字母表示具有显著性差异，P＜0.05。
表2 多重比较检验
Table 2 Multiple comparisons test
种类己酸十六碳酸顺-9-十六碳一烯酸十七碳酸顺-10-十七碳一烯酸产地陕西和山东陕西和四川陕西和甘肃山东和四川山东和甘肃四川和甘肃陕西和山东陕西和四川陕西和甘肃山东和四川山东和甘肃四川和甘肃陕西和山东陕西和四川陕西和甘肃山东和四川山东和甘肃四川和甘肃陕西和山东陕西和四川陕西和甘肃山东和四川山东和甘肃四川和甘肃陕西和山东陕西和四川陕西和甘肃山东和四川山东和甘肃四川和甘肃平均秩差显著性 P 值10.58否 >0.999 9 18.4否0.098 9 10.44否 >0.999 9 7.813否 >0.999 9-0.145 8 否 >0.999 9-7.958否 >0.999 9-4.167否 >0.999 9 11.98否0.817 9 31.02*** 0.000 7 16.15否0.268 0 35.19**** <0.000 1 19.04否0.107 3 3.292否 >0.999 9 28.38** 0.002 2 46.67**** <0.000 1 25.08** 0.009 6 43.38**** <0.000 1 18.29否0.128 6 3.75否 >0.999 9-14.73** 0.008 6 1.979否 >0.999 9-18.48*** 0.000 4-1.771否 >0.999 9 16.71** 0.001 8-2.583否 >0.999 9-2.042否 >0.999 9 5.792否 >0.999 9 0.541 7否 >0.999 9 8.375否0.358 2 7.833否0.469 2
续表2 多重比较检验
Continue table 2 Multiple comparisons test
种类产地平均秩差显著性 P 值十八碳酸陕西和山东-9.646否 >0.999 9顺-9-十八碳一烯酸顺-顺-9，12-十八碳二烯酸二十碳酸陕西和四川陕西和甘肃山东和四川山东和甘肃四川和甘肃陕西和山东陕西和四川陕西和甘肃山东和四川山东和甘肃四川和甘肃陕西和山东陕西和四川陕西和甘肃山东和四川山东和甘肃四川和甘肃陕西和山东陕西和四川17.81否 0.160 0 29.25** 0.001 6 27.46** 0.003 8 38.9**** <0.000 1 11.44否0.928 2-4.458否 >0.999 9 28.58** 0.002 3 29.21** 0.001 7 33.04*** 0.000 2 33.67*** 0.000 2 0.625否 >0.999 9-16.38否0.250 3 16.33否0.253 5 28.21** 0.002 7 32.71*** 0.000 3 44.58**** <0.000 1 11.88否0.838 5-31.65*** 0.000 4-13.33否0.542 7陕西和甘肃山东和四川 18.31否 0.120 3山东和甘肃 36.54**** <0.000 1四川和甘肃 18.23否 0.123 8顺-顺-顺-6，9，12-十八碳 陕西和山东-0.333 3 否 >0.999 9三烯酸陕西和四川 33.25*** 0.000 2陕西和甘肃 51.92**** <0.000 1山东和四川 33.58*** 0.000 2山东和甘肃 52.25**** <0.000 1四川和甘肃 18.67否 0.121 6顺-11-二十碳一烯酸陕西和山东 -20.21否 0.071 8陕西和四川 2.917否 >0.999 9陕西和甘肃 26.29** 0.006 5山东和四川 23.13*0.024 2山东和甘肃 46.5**** <0.000 1四川和甘肃 23.38*0.021 9顺-顺-顺-9，12，15-十八碳陕西和山东-4否 0.266 3三烯酸陕西和四川0否 >0.999 9陕西和甘肃0否 >0.999 9山东和四川4否 0.266 3山东和甘肃4否 0.266 3 4.896否 >0.999 9四川和甘肃0否 >0.999 9二十一碳酸陕西和山东 -38.04 **** <0.000 1陕西和四川 -31.54*** 0.000 5陕西和甘肃 -30.33*** 0.000 8山东和四川6.5否 >0.999 9山东和甘肃 7.708否 >0.999 9四川和甘肃 1.208否 >0.999 9
续表2 多重比较检验
Continue table 2 Multiple comparisons test
注：*表示差异显著，P<0.05，**表示差异极显著，P<0.01，***表示差异高度显著，P<0.001，****表示差异极高度显著P<0.000 1。
种类二十二碳酸二十三碳酸顺-5，8，11，14，17-二十碳五烯酸产地陕西和山东陕西和四川陕西和甘肃山东和四川山东和甘肃四川和甘肃陕西和山东陕西和四川陕西和甘肃山东和四川山东和甘肃四川和甘肃陕西和山东陕西和四川陕西和甘肃山东和四川山东和甘肃四川和甘肃平均秩差显著性 P 值-21.38**0.003 4-20.71**0.005 1-7.917否 >0.999 9 0.666 7否 >0.999 9 13.46否0.180 6 12.79否0.235 6-40.98 **** <0.000 1-23.27*0.022 3-4.75否 >0.999 9 17.71否0.163 8 36.23**** <0.000 1 18.52否0.125 8-6.042否0.573 5-6.083否0.559 9-1.875否 >0.999 9-0.041 67 否 >0.999 9 4.167否 >0.999 9 4.208否 >0.999 9
由图2 可知，在方差分析多重比较结果中，己酸、顺-10-十七碳一烯酸、顺-顺-顺-9，12，15-十八碳三烯酸和顺-5，8，11，14，17-二十碳五烯酸共4 种脂肪酸在陕西、山东、四川和甘肃4 省份花椒样品中无显著性差异，其余12 种脂肪酸均有不同程度的显著性差异（见表2），可以作为花椒产地鉴别的有效指标。通过使用IBM SPSS Statistics 22.0 分析来自陕西、山东、四川和甘肃省花椒样品之间脂肪酸数据的显著性差异，进行进一步的验证，统计显著性结果与GraphPad Prism 7.00 一致。
对4 省份共32 份花椒样品脂肪酸进行层次聚类分析，结果见图3。
图3 4 个省份不同花椒聚类图
Fig.3 Dendrogram of Zanthoxylum bungeanum from four provinces
0.陕西；1.山东；2.四川；3.甘肃。
层次聚类分析图以更加简单、直观的方式反映各省份花椒样品脂肪酸种类和含量的差异。结果表明，陕西和山东、四川和甘肃的部分样品有交叉。
基于花椒样品的脂肪酸组分建立热图，结果见图4a。
图4 4 个省份热图和相关性热图分析
Fig.4 Heat map and correlation heat map analysis of Zanthoxylum bungeanum from four provinces
a.热图；b.相关性热图分析。2．己酸甲酯；13.顺-9-十六碳一烯酸甲酯；14.十七碳酸甲酯；15.顺-10-十七碳一烯酸甲酯；16.十八碳酸甲酯；18．顺-9-十八碳一烯酸甲酯；20.顺，顺-9，12-十八碳二烯酸甲酯；21.二十碳酸甲酯；22.顺-顺-顺-6，9，12-十八碳三烯酸甲酯；23．顺-11-二十碳一烯酸甲酯；24.顺，顺，顺-9，12，15-十八碳三烯酸甲酯；25.二十一碳酸甲酯；27.二十二碳酸甲酯；32.二十三碳35.顺-5，8,11，14，17-二十碳五烯酸甲酯。
在相关性热图分析中，相关系数r 是研究变量之间线性相关程度的指标，当公式中X（产地）和Y（脂肪酸含量）的相关程度最高时，r 接近1，颜色为红色。相关性热图分析见图4b。
由图4a 可知，陕西、山东、四川和甘肃4 个省份的样品中顺-9-十六碳一烯酸、十八碳酸和顺-顺-顺-6，9，12-十八碳三烯酸3 种脂肪酸的变化趋势相近；顺-10-十七碳一烯酸、顺-9-十八碳一烯酸、顺-顺-9，12-十八碳二烯酸、顺-顺-顺-9，12，15-十八碳三烯酸、顺-11-二十碳一烯酸和顺-5，8，11，14，17-二十碳五烯酸6 种脂肪酸的变化趋势相近；二十一碳酸、二十二碳酸、二十三碳酸3 种脂肪酸的变化趋势相近。
由图4b 可知，4 个省份的花椒样品中顺-9-十八碳一烯酸、顺-顺-9，12-十八碳二烯酸、顺-顺-顺-9，12，15-十八碳三烯酸、顺-11-二十碳一烯酸、顺-5，8，11，14，17-二十碳五烯酸共5 种脂肪酸的含量呈正相关，相关性程度较高。
3 结论
陕西、山东、四川和甘肃4 个省份的花椒样品共鉴定出16 种脂肪酸，其中共有脂肪酸12 种。四川和山东样品均检测出15 种脂肪酸；陕西共检测出14 种脂肪酸；甘肃共检测出13 种脂肪酸，在4 省中种类最少。方差分析和多重比较结果表明，不同省份花椒的脂肪酸种类大致相同，但其中12 种脂肪酸含量有不同程度的显著性差异。
植物源性食品中脂肪酸组成和含量与其生长地域密切相关，但各项脂肪酸指标易受气候、土壤等因素的影响，尽管产地不同，但纬度、气候等因素相近可能造成脂肪酸指标的相似性。基于不同省份花椒脂肪酸种类和含量进行比较和分析，层级聚类分析表明陕西和山东省、四川和甘肃省的部分花椒样品在脂肪酸种类和含量方面具有相似性；进一步建立热图并进行相关性分析，结果表明陕西、山东、四川和甘肃4 个省份的花椒样品中共有5 种脂肪酸的含量相关性程度较高，而变化趋势相近的脂肪酸成分有所差异。
[1] 袁海梅,邱露,谢贞建,等.花椒属植物生物碱类成分及其药理活性研究进展[J].中国中药杂志,2015,40(23):4573-4584.YUAN Haimei,QIU Lu,XIE Zhenjian,et al.Research progress on alkaloids constituents from Zanthoxylum and their pharmacological activities[J].China Journal of Chinese Materia Medica,2015,40(23):4573-4584.
[2] 张国琳.基于风味成分的花椒品质评价研究[D].南京:南京农业大学,2014.ZHANG Guolin.Quality evaluation of Zanthoxylum bungeanum based on the analysis of flavor components[D].Nanjing:Nanjing Agricultural University,2014.
[3] 曹蕊,杨潇,蒋珍菊,等.红花椒油和青花椒油中脂肪酸组成的GC-MS 的对比研究[J].中国调味品,2011,36(8):102-105.CAO Rui,YANG Xiao,JIANG Zhenju,et al.The GC-MS contrast study of fatty acids in essential oil of Zanthoxylum bungeanum and Zanthoxylum schinifolium[J].China Condiment,2011,36(8): 102-105.
[4] 凌智群,程宝宏,魏居国,等.花椒功效发展的历史沿革[J].云南中医学院学报,2008,31(4):49-50.LING Zhiqun,CHENG Baohong,WEI Juguo,et al.Historical evolution of the efficacy development of Zanthoxylum bungeanum[J].Journal of Yunnan University of Traditional Chinese Medicine,2008,31(4):49-50.
[5] 曾京京.我国花椒的栽培起源和地理分布[J].中国农史,2000,19(4):68-75.ZENG Jingjing.Cultivation origin and geographical distribution of Zanthoxylum bungeanum in China[J].Agricultural History of China,2000,19(4):68-75.
[6] 张继,白贞芳,杨永利,等.甘肃省花椒属药用植物资源[J].中国药学杂志,2002,37(4):255-257.ZHANG Ji,BAI Zhenfang,YANG Yongli,et al.Medicinal plant resources of Zanthoxylum in Gansu Province[J].Chinese Pharmaceutical Journal,2002,37(4):255-257.
[7] PHUYAL N,JHA P K,RATURI P P,et al.In vitro antibacterial activities of methanolic extracts of fruits,seeds,and bark of Zanthoxylum armatum DC[J].Journal of Tropical Medicine,2020,2020:2803063.
[8] 寇芸芸.花椒成分分析和杀虫活性研究[D].武汉: 武汉轻工大学,2015.KOU Yunyun.The extraction process and insecticidal activity research of Zanthoxylum[D].Wuhan: Wuhan Polytechnic University,2015.
[9] LU Q,MA R F,YANG Y,et al.Zanthoxylum nitidum (roxb.) DC:Traditional uses,phytochemistry,pharmacological activities and toxicology[J].Journal of Ethnopharmacology,2020,260:112946.
[10] 李晴,许腾腾,张佳慧,等.花椒提取物对嗜水气单胞菌群体感应的抑制作用[J].生物加工过程,2020,18(2):263-268.LI Qing,XU Tengteng,ZHANG Jiahui,et al.Inhibition of Aeromonas hydrophila quorum sensing by extract from Zanthoxylum bungeanum[J].Chinese Journal of Bioprocess Engineering,2020,18(2):263-268.
[11] 宋荣,曹亮,周佳民,等.花椒种质资源及其功能成分和生物学效应研究进展[J].湖南农业科学,2014(17):23-26,29.SONG Rong,CAO Liang,ZHOU Jiamin,et al.Research progress in Zanthoxylum germplasm resources and chemical composition and biological effects[J].Hunan Agricultural Sciences,2014(17):23-26,29.
[12] WANG Z Y,HAN Q Q,DENG M Y,et al.Lemairamin,isolated from the Zanthoxylum plants,alleviates pain hypersensitivity via spinal α7 nicotinic acetylcholine receptors[J].Biochemical and Biophysical Research Communications,2020,525(4):1087-1094.
[13] HWANG K A,KWON J E,NOH Y,et al.Effects of Zanthoxylum piperitum ethanol extract on osteoarthritis inflammation and pain[J].Biomedicine&Pharmacotherapy,2018,105:481-490.
[14] 袁太宁,王艳林,汪鋆植.花椒挥发油抗宫颈癌Hela 细胞作用研究[J].湖北民族学院学报(医学版),2008,25(3):26-27.YUAN Taining,WANG Yanlin,WANG Junzhi.The study of the anti-cervix cancer hela cells effect of zanthoxylunm maxim[J].Journal of Hubei Institute for Nationalities (Medical Edition),2008,25(3):26-27.
[15] PANG W W,LIU S,HE F T,et al.Anticancer activities of Zanthoxylum bungeanum seed oil on malignant melanoma[J].Journal of Ethnopharmacology,2019,229:180-189.
[16] TIAN Y Q,ZHANG C Y,GUO M Q.Comparative study on alkaloids and their anti -proliferative activities from three Zanthoxylum species[J].BMC Complementary and Alternative Medicine,2017,17(1):460.
[17] ALAM F,NAJUM US SAQIB Q,WAHEED A.Cytotoxic activity of extracts and crude saponins from Zanthoxylum armatum DC.against human breast (MCF-7,MDA-MB-468) and colorectal (Caco-2)cancer cell lines[J].BMC Complementary and Alternative Medicine,2017,17(1):368.
[18] 游玉明,周敏,王倩倩,等.花椒麻素的抗氧化活性[J].食品科学,2015,36(13):27-31.YOU Yuming,ZHOU Min,WANG Qianqian,et al.Antioxidant activity in vitro of sanshool from Sichuan peppers (Zanthoxylum bungeanum)[J].Food Science,2015,36(13):27-31.
[19] 王孟,吕倩,董林娟.花椒醇提物抗氧化性分析及应用[J].盐科学与化工,2020,49(3):13-16.WANG Meng,LV Qian,DONG Linjuan.Antioxidant analysis and application of alcohol extract from Zanthoxylum bungeanum [J].Journal of Salt Science and Chemical Industry,2020,49(3):13-16.
[20] ROUT P K,NAIK S N,RAO Y R,et al.Extraction and composition of volatiles from Zanthoxylum rhesta:Comparison of subcritical CO2 and traditional processes[J].The Journal of Supercritical Fluids,2007,42(3):334-341.
[21] OCHIENG C O,NYONGESA D W,YAMO K O,et al.α-Amylase and α-glucosidase inhibitors from Zanthoxylum chalybeum Engl.root bark[J].Fitoterapia,2020,146:104719.
[22] 韩丽,高燕,公维娟,等.花椒中生物碱的提取与含量测定[J].化学工程与装备,2019(12):9-11.HAN Li,GAO Yan,GONG Weijuan,et al.Extraction and content determination of alkaloids in Zanthoxylum bungeanum[J].Chemical Engineering&Equipment,2019(12):9-11.
[23] 张东博,王征,梁艳妮,等.青花椒根皮中1 个新的香豆素类化合物[J].中草药,2019,50(12):2798-2801.ZHANG Dongbo,WANG Zheng,LIANG Yanni,et al.A new coumarin isolated from root barks of Zanthoxylum schinifolium[J].Chinese Traditional and Herbal Drugs,2019,50(12):2798-2801.
[24] 蔡沁玥.竹叶花椒果皮木脂素分离、鉴定、苦度分析及其体外抗氧化活性研究[D].雅安:四川农业大学,2019.CAI Qinyue.Study on separation,identification,bitterness analysis and antioxidant activity of lignan from Zanthoxylum armatum DC[D].Ya'an:Sichuan Agricultural University,2019.
[25] 林云,邓娟,周玉碧.花椒籽油中脂肪酸成分的GC-MS 分析[J].甘肃农业大学学报,2015,50(6):132-135.LIN Yun,DENG Juan,ZHOU Yubi.Analysis of fatty acids in the prickly ash seed oil by GC-MS[J].Journal of Gansu Agricultural University,2015,50(6):132-135.
[26] 张翔.花椒挥发油和脂肪酸的GC 指纹图谱研究[D].杨凌:西北农林科技大学,2011.ZHANG Xiang.Study on the GC fingerprint of essential oil and fatty acids in Zanthoxylum[D].Yangling: Northwest A & F University,2011.
[27] 唐玉莲,林红卫,陈迪钊,等.武陵青花椒壳挥发油成分及脂肪酸分析[J].食品工业科技,2012,33(10):74-78.TANG Yulian,LIN Hongwei,CHEN Dizhao,et al.Analysis of chemical components of volatile oil and fatty acid from pericarps of Wuling Zanthoxylum schinifolium Sieb.Et Zucc[J].Science and Technology of Food Industry,2012,33(10):74-78.
[28] 候丽秀.花椒籽油的化学计量学分类及脂肪酸产地差异研究[D].杨凌:西北农林科技大学,2019.HOU Lixiu.Chemometric classification and habitat variations in the fatty acids of Zanthoxylum seed oils[D].Yangling:Northwest A&F University,2019.
[29] 张仁凤,陈光静,屈立武,等.不同产地红花椒籽油脂肪酸组成的比较研究[J].中国调味品,2017,42(6):1-4.ZHANG Renfeng,CHEN Guangjing,QU Liwu,et al.Comparative analysis of fatty acids composition of seed oil in Zanthoxylum bungeanum maxim.from different orgins[J].China Condiment,2017,42(6):1-4.
[30] 国家卫生和计划生育委员会,国家食品药品监督管理总局.食品安全国家标准食品中脂肪酸的测定:GB 5009.168—2016[S].北京:中国标准出版社,2017.National Health and Family Planning Commission,State Food and Drug Administration,National standard for food safety Determination of eicosapentaenoic acid and docosahexaenoic acid in foods:GB 5009.168—2016[S].Beijing:Standards Press of China,2017.
[31] 范胜栩,李斌,孙君明,等.气相色谱方法定量检测大豆5 种脂肪酸[J].中国油料作物学报,2015,37(4):548-553.FAN Shengxu,LI Bin,SUN Junming,et al.A quantitative gas chromatographic method for determination of soybean seed fatty acid components[J].Chinese Journal of Oil Crop Sciences,2015,37(4):548-553.
[32] 刘皓月,李萌,朱庆珍,等.3 种萃取方法炸蒜油特征风味的比较分析[J].食品科学,2020,41(12):180-187.LIU Haoyue,LI Meng,ZHU Qingzhen,et al.Comparison of three extraction methods used to determine the characteristic flavor of garlic frying oil[J].Food Science,2020,41(12):180-187.
[33] 张伟,尹震花.HS-SPME-GC/MS 分析鸡矢藤挥发性成分[J].中国实验方剂学杂志,2015,21(23):55-57.ZHANG Wei,YIN Zhenhua.Analysis of volatile constituents in Paederia scandens by HS-SPME-GC/MS[J].Chinese Journal of Experimental Traditional Medical Formulae,2015,21(23):55-57.