葉洵，劉子博，張婷，張萌萌，何林，王丹，彭偉，吳純潔

(成都中醫(yī)藥大學(xué) 藥學(xué)院，成都 611137)

花椒為蕓香科植物花椒(ZanthoxylumbungeanumMaxim.)的干燥成熟果皮[1]?；ń肥羌彝コＳ谜{(diào)味品，深受人們的喜愛，最早在春秋戰(zhàn)國時期記載使用?；ń泛胸S富的揮發(fā)油、酰胺等活性物質(zhì)，在醫(yī)療領(lǐng)域也有一定的應(yīng)用價值，具有溫中行氣、驅(qū)寒、止痛、殺蟲等功效，可用于緩解胃痛、嘔吐以及腹瀉等癥狀。到目前為止，從科學(xué)研究到生產(chǎn)，研究人員和操作人員對花椒的分類依據(jù)都不同，叫法各異，這導(dǎo)致了國內(nèi)花椒的品種和名稱的混亂，給生產(chǎn)和銷售帶來諸多不便，產(chǎn)品質(zhì)量不一[2-3]。花椒商品依據(jù)顏色分類，分為“紅花椒”和“青花椒”。其中，“紅花椒”的知名產(chǎn)區(qū)為四川漢源和茂汶、陜西韓城、甘肅武都等地區(qū)，而“青花椒”的知名產(chǎn)區(qū)則是四川金陽和洪雅、重慶江津等地區(qū)。根據(jù)栽培品種分類，主要有漢源貢椒、正路花椒、大紅袍、大紅椒、小紅椒、金陽青花椒、洪雅藤椒等品種[4]?；ń放c中藥材性質(zhì)類似，其質(zhì)量受種植地的經(jīng)緯度、氣候(日照時間、降雨量、干燥度等)以及土壤等因素的影響較大，各地種植的花椒品質(zhì)差異較大，導(dǎo)致花椒揮發(fā)油成分含量也有所不同[5-6]，揮發(fā)性香味成分是構(gòu)成和影響其加工產(chǎn)品質(zhì)量的重要因素。

指紋圖譜技術(shù)多用于鑒別道地藥材，因此可采用該技術(shù)建立花椒揮發(fā)油指紋圖譜，作為花椒產(chǎn)地鑒別的依據(jù)[7]。本研究采用水蒸氣蒸餾的方法提取花椒揮發(fā)油，建立花椒揮發(fā)油GC-MS指紋圖譜，分析揮發(fā)油化學(xué)成分差異，結(jié)合保留指數(shù)法、聚類分析和主成分分析對不同產(chǎn)地的花椒揮發(fā)油成分進行研究。

1 材料與方法

1.1 材料與試劑

10批花椒樣品具體信息見表1，于2019年8月-10月采集，經(jīng)成都中醫(yī)藥大學(xué)藥學(xué)院李敏教授鑒定為蕓香科植物花椒(ZanthoxylumbungeanumMaxim.)的干燥成熟果皮。

表1 10批花椒藥材名稱及來源Table 1 The names and sources of 10 batches of Z. bungeanum

正己烷、無水硫酸鈉(分析純)：購自成都市科隆化學(xué)品有限公司；正構(gòu)烷烴(C10～C25)：購自成都埃法生物科技有限公司。

1.2 儀器與設(shè)備

揮發(fā)油提取器 四川蜀牛玻璃儀器有限公司；QE-100二兩裝高速粉碎機 浙江省蘭溪市偉能達電器有限公司；Agilent 7890A-5975C氣相色譜-質(zhì)譜聯(lián)用儀、Agilent GC-MSD工作站 美國安捷倫公司；Sartorius BP211D電子天平 德國賽多利斯公司；MH-1000調(diào)溫型電熱套 北京科偉永興儀器有限公司。

2 實驗方法與結(jié)果

2.1 指紋圖譜研究

2.1.1 色譜-質(zhì)譜條件

氣相色譜柱為HP-5MS毛細管質(zhì)譜柱(30 m×250 μm×0.25 μm)；載氣為氦氣(He)，流速為1 mL/min；程序升溫：初始溫度為50 ℃，保持2 min，以5 ℃/min升至150 ℃，最后以25 ℃/min升至240 ℃；進樣量為1 μL，采用分流模式，分流比為30∶1；溶劑延遲5 min；離子源溫度為230 ℃；EI電離，電子能量為70 eV；掃描模式為全掃描模式；質(zhì)譜數(shù)據(jù)標(biāo)準(zhǔn)譜庫為NIST 14.L。在相同實驗條件下檢測所有供試品，記錄樣品GC-MS總離子流圖。

2.1.2 花椒揮發(fā)油提取

按照2020年版《中國藥典》四部(通則2204)中揮發(fā)油測定法，采用揮發(fā)油蒸餾法提取花椒揮發(fā)油。具體操作：花椒樣品粉碎通過2號、3號篩混合均勻。取粉末適量20 g，稱定重量(準(zhǔn)確至0.01 g)，置于500 mL圓底燒瓶中，加10倍量水與玻璃珠數(shù)粒，連接揮發(fā)油測定器與回流冷凝管，放入電熱套中緩緩加熱至沸騰，并保持微沸約5 h，至測定器中油量不再增加，停止加熱，冷卻后，讀取揮發(fā)油體積，并計算供試品中揮發(fā)油的含量。加適量無水硫酸鈉放于-20 ℃保存過夜。

2.1.3 樣品溶液制備

花椒樣品溶液：精密吸取揮發(fā)油溶液50 μL，加正己烷定容至5 mL，0.22 μm微孔濾膜過濾，即得供試品溶液。

正構(gòu)烷烴溶液：精密吸取C10～C25正構(gòu)烷烴100 μL，加正己烷定容至1 mL，搖勻，0.22 μm微孔濾膜過濾，即得0.1 mg/mL正構(gòu)烷烴溶液。

2.1.4 保留指數(shù)(RI)的計算

按照2.1.1項實驗條件，吸取花椒樣品溶液與正構(gòu)烷烴溶液進行檢測，獲得總離子色譜圖。根據(jù)二者的保留時間，計算待鑒定化合物的保留指數(shù)RI，將獲得的RI值與文獻RI值進行比較。

保留指數(shù)的計算公式為：

(1)

式中：tR為保留時間，min；x為供試品化合物；Z、Z+1為待分析物前后兩個正構(gòu)烷烴的碳原子數(shù)，即tR(Z)

2.1.5 精密度考察

取1份供試品溶液，按照2.1.1的色譜條件，連續(xù)進樣6次，進樣量1 μL，記錄總離子流圖。以芳樟醇為參照峰，計算共有峰保留時間RSD值均小于3%，峰面積RSD值均小于3%，表明儀器具有良好的精密度。

2.1.6 穩(wěn)定性考察

取2.1.3項的供試品溶液，按照2.1.1的色譜條件，于0，2，4，8，12，24 h測定峰面積，以芳樟醇為參照峰，計算共有峰保留時間RSD值均小于3%，峰面積RSD值均小于3%，表明花椒揮發(fā)油在24 h內(nèi)穩(wěn)定。

2.1.7 重復(fù)性考察

取6份同一批次花椒藥材粉末，按照2.1.2和2.1.3中供試品溶液的制備方法制備，按照2.1.1的色譜條件，進樣1 μL，記錄總離子流色譜圖。以芳樟醇為參照峰，計算其他峰保留時間RSD值均小于3%，峰面積RSD值小于3%，表明該方法重復(fù)性較好，適合花椒指紋圖譜的研究。

2.1.8 指紋圖譜的建立

2.1.8.1 花椒揮發(fā)油成分共有峰確認

GC-MS總離子流圖顯示：花椒揮發(fā)油確定了28～37個峰，共有峰23個，共有峰面積占總峰面積86%以上，以D-Limonene、Linalool、Geranyl acetate、Caryophyllene、γ-Elemene、Humulene、Germacrene D、β-Myrcene、α-Terpineol等為主，具體結(jié)果見表2。

表2 10批花椒揮發(fā)油GC-MS分析結(jié)果Table 2 The analysis results of 10 batches of Z. bungeanum volatile oils by GC-MS

2.1.8.2 指紋圖譜的建立

將其導(dǎo)入“中藥色譜指紋圖譜相似度評價系統(tǒng)”(2012版)，以表1中M2為參照圖譜，時間寬度設(shè)為0.10 min，多點校正后，生成花椒揮發(fā)油共有模式對照指紋圖譜，標(biāo)定23個峰，10批花椒揮發(fā)油樣品指紋圖譜疊加圖見圖1，對照圖譜見圖2。

圖1 10批花椒揮發(fā)油樣品指紋圖譜疊加對照圖Fig.1 The fingerprint spectra overlay reference diagram of 10 batches of Z. bungeanum volatile oil samples

圖2 花椒揮發(fā)油對照指紋圖譜Fig.2 The reference fingerprint spectra of Z. bungeanum volatile oils

各批次花椒揮發(fā)油色譜圖與對照圖譜相似度結(jié)果見表3。

表3 10批花椒揮發(fā)油對照圖譜相似度Table 3 The similarity of reference fingerprint spectra of 10 batches of Z. bungeanum volatile oils

由表3可知，10批花椒揮發(fā)油樣品的相似度較好，均大于0.9，表明各產(chǎn)地的花椒質(zhì)量具有較高的一致性。

2.2 不同產(chǎn)地花椒揮發(fā)油成分差異研究

揮發(fā)油是花椒香氣的主要成分，花椒揮發(fā)油含量的多少可直接反映花椒香味程度?；ń窊]發(fā)油含量是反映花椒風(fēng)味物質(zhì)的一項重要指標(biāo)。根據(jù)式(2)計算樣品中揮發(fā)油含量。

(2)

花椒揮發(fā)油含量測定結(jié)果見表4。

表4 10批花椒揮發(fā)油含量的分析結(jié)果Table 4 The analysis results of 10 batches of Z. bungeanum volatile oil content

由表4可知，不同產(chǎn)地花椒揮發(fā)油含量在2.8%～6.7%之間。不同產(chǎn)地的花椒揮發(fā)油含量有差異，且同一產(chǎn)地的花椒揮發(fā)油含量也有所不同。其中10批花椒樣品中，揮發(fā)油含量最多的是茂縣花椒，為6.7%，含量最少的是甘肅武都，為2.8%，但均滿足2020年版《中國藥典》規(guī)定的揮發(fā)油含量不低于1.5%。

2.3 不同產(chǎn)地花椒揮發(fā)油聚類分析

系統(tǒng)聚類分析可根據(jù)樣品之間的親疏遠近關(guān)系進行分類，樣品的相似度越大，兩者之間的距離越近。適合樣品數(shù)量較少時的聚類分析，最終輸出樹狀圖的聚類結(jié)果。根據(jù)本實驗所得花椒總離子流色譜圖，經(jīng)過化學(xué)成分差異分析，分析不同產(chǎn)地花椒樣品揮發(fā)油的共有峰，并將共有峰峰面積作為花椒揮發(fā)油聚類分析的分析對象，將數(shù)據(jù)導(dǎo)入SPSS 25.0軟件，采用瓦爾德法，并用歐氏距離的平方為測度對其進行聚類分析。

系統(tǒng)聚類分析結(jié)果見圖3。

圖3 10批花椒揮發(fā)油聚類分析結(jié)果Fig.3 The cluster analysis results of 10 batches of Z. bungeanum volatile oils

由圖3可知，當(dāng)類間距離為20時，10批樣品分為兩大類；H1、H2、H3為一類，其他聚為一類。當(dāng)類間距離為10時，10批樣品被分為三大類：M4聚為一類，H1、H2、H3聚為一類，其余為一類。同一產(chǎn)地的花椒樣品揮發(fā)油相似度較好，但是亦被聚類分析劃分為不同的類別，說明同一產(chǎn)地花椒既具有一致性又存在一定的差異。

2.4 不同產(chǎn)地花椒揮發(fā)油主成分分析

主成分分析是一種降維的數(shù)據(jù)分析方法，根據(jù)樣品之間的相似度將樣品分為若干組。為了找出不同產(chǎn)地花椒揮發(fā)油質(zhì)量差異性，利用SIMCA-P 14.1軟件對10批花椒揮發(fā)油共有峰峰面積進行主成分分析，以23個共有峰峰面積為變量進行PCA運算，共生成4個主成分，累積貢獻率達到97.36%，其中前2個主成分累積貢獻率達到83.23%，證明模型有效。以PCA得分和峰面積變量得到Score圖(見圖4)和Loading圖(見圖5)。漢源貢椒與茂縣大紅袍、武都大紅袍沿著t[1]軸能明顯分開，分為兩大類，各組內(nèi)差異較小，無異常離群樣本，說明不同產(chǎn)地花椒揮發(fā)油化學(xué)成分具有較大差異。根據(jù)PCA分析結(jié)果，將漢源、茂縣和武都3個產(chǎn)地的樣品大致分為2組，以23個共有指紋峰峰面積積分值為變量進行PLS-DA分析，得到23個變量的VIP圖(見圖6)(R2Y=0.927，Q2=0.703，模型有效)，VIP可用來描述PLS-DA分析中變量的變異程度，通常選VIP>1的變量為特征變量。以VIP>1.0為顯著影響，共找到10個差異標(biāo)志物，對其影響顯著性排序，分別為峰22>峰23>峰19>峰18>峰10>峰16>峰11>峰6>峰14>峰20>峰1，說明不同產(chǎn)區(qū)的花椒揮發(fā)油中Germacrene D、(+)-delta-Cadinene、Caryophyllene、Geranyl acetate、Linalool、α-Terpinyl acetate、(E,Z)-2,6-Dimethylocta-2,4,6-triene、trans-β-Ocimene、Nerol和γ-Elemene的含量存在顯著差異，可作為不同產(chǎn)地花椒揮發(fā)油質(zhì)量差異的主要共有峰。

圖4 10批花椒揮發(fā)油Score圖Fig.4 The score diagram of 10 batches of Z. bungeanum volatile oils

圖5 10批花椒揮發(fā)油Loading圖Fig.5 The loading diagram of 10 batches of Z. bungeanum volatile oils

圖6 10批花椒揮發(fā)油VIP圖Fig.6 The VIP diagram of 10 batches of Z. bungeanum volatile oils

3 討論

揮發(fā)油作為花椒風(fēng)味物質(zhì)的主要成分之一，本研究采用GC-MS技術(shù)，測定花椒揮發(fā)油總離子流圖，并建立指紋圖譜。結(jié)果表明花椒揮發(fā)油相似度大于0.90，證明指紋圖譜建立的方法有效。說明不同產(chǎn)地的花椒揮發(fā)油既有一定相似性，又存在明顯差異。茂縣與武都大紅袍相似度在0.929～0.992，揮發(fā)油含量在2.8%～6.7%，兩者跨度較大，可能是由于采收季節(jié)時間、保存方式等的不一樣導(dǎo)致。聚類分析、主成分分析結(jié)果表明不同產(chǎn)地花椒揮發(fā)油成分之間存在差異，該結(jié)論與課凈璇等的研究結(jié)果一致[8]。

本研究以當(dāng)前為主流市場品種的漢源花椒、武都花椒以及茂縣花椒的10批樣品進行GC-MS指紋圖譜研究及相似度評價，共標(biāo)定了23個共有峰。聚類分析與主成分分析結(jié)果發(fā)現(xiàn)，漢源貢椒歸為一類；茂縣與武都花椒二者聚為一類，可能是由于兩地的花椒品種均為大紅袍，在地理上具有一定的相似性。漢源花椒被列為皇室貢品1000余年，成為中國進貢史長之最，漢源花椒依然是我國品質(zhì)最好的花椒之一。漢源花椒土壤條件以黃棕壤為主，海拔在1100～1800 m。隴南大紅袍花椒主要分布于武都、文縣、舟曲等地，隴南武都區(qū)光熱資源豐富，氣候十分適合花椒生長，被冠以“千年椒鄉(xiāng)”[9-10]，花椒果實迅速膨大，椒皮鮮紅，油腺多而密，麻味重，粒大飽滿，香氣濃郁，品質(zhì)特異，被稱為“大紅袍”。而茂縣干旱河谷地區(qū)晝夜溫差大的獨特氣候條件，塑造了茂縣特有的“大紅袍”花椒鮮紅透亮、麻味素含量高的特征。武都、茂縣花椒產(chǎn)地范圍內(nèi)海拔900～2600 m，土壤為砂壤土，土壤pH值6.5～7.5。茂縣與武都地理環(huán)境均適宜“大紅袍”品種的栽種，故在本實驗中采集的大紅袍樣品出現(xiàn)了聚為一類的現(xiàn)象。

花椒作為高級調(diào)味香料和其深加工產(chǎn)品，其揮發(fā)性香味成分不僅是作為香原料的重要考察指標(biāo)，也是構(gòu)成和影響其加工產(chǎn)品質(zhì)量的重要因素[11]。目前，花椒質(zhì)量等級按照林業(yè)行業(yè)標(biāo)準(zhǔn)(LY/T 1652-2005)、國標(biāo)GB/T 30391-2019等分級，主要包括感官評價和理化指標(biāo)評價，其中感官指標(biāo)包括色澤、氣味、味道和果形特征等，理化指標(biāo)則以揮發(fā)油含量為主[12]。本研究建立了基礎(chǔ)指紋圖譜數(shù)據(jù)的花椒整體質(zhì)量評價的方法，可找出質(zhì)量影響較大的差異性成分，為后續(xù)選擇指標(biāo)性成分進行含量測定來評價花椒質(zhì)量優(yōu)劣與揮發(fā)性成分相關(guān)性研究奠定了理論與實驗基礎(chǔ)。同時不同產(chǎn)地商品的花椒揮發(fā)油中Germacrene D、(+)-delta-Cadinene、Caryophyllene、Geranyl acetate、Linalool、α-Terpinyl acetate、(E,Z)-2,6-Dimethylocta-2,4,6-triene、trans-β-Ocimene、Nerol和γ-Elemene的含量存在顯著差異，可作為不同產(chǎn)地花椒揮發(fā)油質(zhì)量差異的主要共有峰。

花椒揮發(fā)油含有大量同分異構(gòu)體，在分析過程中僅利用標(biāo)準(zhǔn)質(zhì)譜庫檢索進行定性，容易出現(xiàn)誤判。對于一些常規(guī)GC-MS定性方法難以鑒別的化合物，采用保留指數(shù)法可以得到更準(zhǔn)確的定性結(jié)果[13]，如：保留時間為10.147 min的色譜峰，通過NIST 14.L標(biāo)準(zhǔn)質(zhì)譜庫檢索分析為γ-Terpinene、3-Carene，匹配度均為94，經(jīng)查閱文獻[14]，γ-Terpinene的保留指數(shù)為1060，3-Carene的保留指數(shù)為1007，與實測值進行比較，可知該化合物為γ-Terpinene。通過質(zhì)譜庫檢索結(jié)果結(jié)合保留指數(shù)法，可準(zhǔn)確性判斷花椒揮發(fā)油中的化學(xué)成分包括順反式以及α，β，γ等各種同分異構(gòu)體。

由于花椒產(chǎn)地較多，本研究僅收集到部分產(chǎn)地的花椒品種，研究結(jié)果也僅能說明漢源、武都以及茂縣花椒揮發(fā)性成分的初步特征及差異，后期還需收集足夠多代表性的樣品，進一步驗證該實驗結(jié)果。本研究建立的花椒揮發(fā)油指紋圖譜對花椒的質(zhì)量控制有積極意義，所含揮發(fā)性共有成分含量存在明顯差異，可為花椒品種選育和其在醫(yī)藥、食品、調(diào)味品等行業(yè)中的應(yīng)用提供科學(xué)依據(jù)。