基于礦物元素指紋分析的藥食兼用百合產地溯源

李瑞琴 白濱 于安芬 徐瑞 郭斌

摘 要：利用iCE-3500原子吸收光譜儀等儀器測定來自甘肅省、湖南省、江西省和安徽省4個省8個市11個縣共119份藥食兼用百合樣品中 10 種礦物元素（Ca、Mg、Zn、Cu、Fe、Se、K、P、Cd、Pb）的含量，對數據進行方差分析、主成分分析和判別分析。結果表明：不同地域來源藥食兼用百合樣品中礦物元素含量不但與其產地的土壤質地、土壤類型密切相關，此外還受品種和栽培措施等因素的影響;確定Ca、Mg、Zn、Cu、Fe為百合的特征元素;建立的判別模型可以明顯區分各產地及不同品種的百合，為名特優新農產品產地溯源及鑒定提供了科學依據。

關鍵詞：藥食兼用百合;品種;產地;礦物元素

百合含有蛋白質、碳水化合物、脂肪、礦物元素等營養物質及多糖、甾體皂苷、酚類化合物、黃酮類化合物、生物堿等生物活性物質，其以干燥的肉質鱗葉入藥，常用于陰虛久咳，痰中帶血，虛煩驚悸等，還具有抗抑郁、抗癌、降血糖、增強機體免疫力等功效[1-7]。不同產地和品種百合中礦物元素的含量差異能反映百合的品質及其產地特征，結合數據統計學方法對其進行分析比較，對百合品質評價、產地判別和品牌驗證具有重要意義。目前我國規模化種植的百合主要有4 種，即蘭州百合、川百合、卷丹百合和龍牙百合，研究以栽培技術、繁育和主要營養成分含量及健康功效等為主[8-11]。礦物元素溯源技術結合統計學方法已經廣泛應用于小麥、茶葉、水果、蜂蜜、葡萄酒產地溯源的研究[12-21]。因此，研究百合礦物元素的種類及其含量，探尋能表征百合地域信息的特異性指標，為產地保護和產品溯源提供技術支撐。

1 材料與方法

1.1 試驗材料

從甘肅省、湖南省、江西省和安徽省采集2016—2017 年度百合樣品。每個省選1～3個主產市，每個市選1～5個主產縣，每個縣選3個主產鄉（鎮），每個鄉（鎮）選種植面積最大的藥食兼用百合基地（表1）。

1.2 試驗方法

1.2.1 采樣方法 在百合收獲期，選擇代表性農戶，直接從田間收獲藥食兼用百合，均勻取樣3 kg，儲藏運輸過程均保持冷鏈狀態。

1.2.2 樣品預處理方法 樣品運回實驗室先除去殘根和有瑕疵的鱗瓣，剝瓣、用水沖洗干凈，隨即瀝干水分，放入烘箱中105℃殺青，隨即降至60℃烘干。再用 1093旋風磨（丹麥 Foss Tecator 公司）磨制得到百合粉，供試驗分析使用。

1.2.3 測試方法 稱取0.05～0.1 g 百合粉，樣品采用HNO3-HClO4法消煮。用iCE-3500原子吸收光譜儀（美國Thermo公司）測定樣品中Ca、Mg、Zn、Cu、Fe、Cd、Pb元素的含量，P，酸溶-釩鉬黃比色法，用CARY50紫外可見分光光度計;K，酸溶-火焰光度法，用M410型火焰光度計;Se，氫化物原子熒光光譜法，用RGF-8740原子熒光光度計，其中Cd、Pb用石墨爐原子吸收法測定。所有樣品的測定均在農業農村部農產品質量安全風險評估實驗室（蘭州）完成。

1.2.4 數據處理方法 用SPSS軟件對數據進行方差分析（Duncan多重比較分析，顯著性水平設為α=0.05）、相關分析、主成分分析和判別分析。

2 結果與分析

2.1 藥食兼用百合中礦物元素含量差異分析

從表2可以看出，同種元素在不同品種百合之間差異較大，含量最高的元素是K、Mg、P、Ca。不同地區元素含量都表現為K、Mg、P、Ca、Fe、Zn、Cu、Pb、Cd、Se依次降低。甘肅蘭州百合的Ca含量最高，與其他產地百合差異顯著;江西龍牙百合的Ca、Mg、Se、P含量最低;湖南龍牙百合的Fe、Cd、Pb含量最低;安徽卷丹百合樣品的Mg、Zn、Cu、Fe、Se、K、P、Cd、Pb共9 種元素含量顯著高于其他地區。但是Cd和Pb重金屬含量均未超過《GB 2762—2019食品安全國家標準食品中污染物限量》中關于多年生蔬菜的限量要求。此結果說明了樣品的來源地域、溯源范圍與礦物元素溯源指標也是密切相關的。

2.2 藥食兼用百合中礦物元素含量相關性分析

如表3所示，百合中有諸多對元素呈極顯著正相關（P<0.01），揭示這些元素間具有相互協同、促進吸收的關系，同時說明百合在富集以上礦物元素時具有較強的協同作用。

2.3 藥食兼用百合中礦物元素含量的主成分分析

主成分分析結果表明，第1～5個主成分的貢獻率分別為41.84%、14.23%、11.85%、8.41%和6.64%，其累計方差貢獻率為82.973%。根據主成分1的貢獻率可以判定Mg、Cu、Fe、Se、K、P、Cd、Pb這8種元素的含量對于百合產地的分類的貢獻較其他元素大（表4）。

2.4 百合中礦物元素含量的判別分析

從不同省份百合樣品各元素指標的方差分析、主成分分析結果中可以對10種差異顯著的元素進行逐步分析，篩選出有效變量，建立如下判別模型。

類別1（湖南卷丹）=-23.436+0.002X1（Ca）+0.014X2（Mg）+0.236X3（Zn）+4.589X4（Cu）-0.082X5（Fe）-23.091 X6（Se）+0.002X7（K）+0.014 X8（P）+66.479X9（Cd）+56.933 X10（Pb）

類別2（安徽卷丹）=-35.03+0.006X1（Ca）+0.016X2（Mg）-0.551X3（Zn）+8.406X4（Cu）-0.082X5（Fe）-18.069 X6（Se）+0.003X7（K）+0.010 X8（P）+98.438X9（Cd）+64.909 X10（Pb）

類別3（江西龍牙）=-45.658+0.002X1（Ca）+0.019X2（Mg）+0.378X3（Zn）+5.207X4（Cu）-0.117X5（Fe）-46.237 X6（Se）+0.002X7（K）+0.032 X8（P）+132.875X9（Cd）+83.755 X10（Pb）

類別4（湖南龍牙）=-109.547+0.003X1（Ca）+0.017X2（Mg）-1.059X3（Zn）+13.829X4（Cu）+0.196X5（Fe）+87.477 X6（Se）+0.004X7（K）+0.036 X8（P）+232.938X9（Cd）+118.601 X10（Pb）

類別5（蘭州百合）=-43.803+0.014X1（Ca）+0.006X2（Mg）-818X3（Zn）+5.797X4（Cu）+0.361X5（Fe）-9.838X6（Se）+0.002X7（K）+0.026 X8（P）+101.636X9（Cd）-58.823 1 X10（Pb）

利用此判別模型對百合樣品進行回代檢驗和交叉檢驗分析，回代檢驗判別率為92.2%，交叉檢驗判別率為87.8%（表5）。不同地區的樣品位于不同的區域，說明礦物元素指紋分析技術是可以用于藥食兼用百合產地溯源的一種有效方法。

3 結論

（1）4個主要產地的藥食兼用百合中的礦物元素含量有其自身的特點。安徽卷丹百合樣品的Mg、Zn、Cu、Fe、Se、K、P、Cd、Pb共9種元素含量顯著高于其他地區，原因可能是因為安徽省的樣品雖然不是同一個產地，但是均為同一個企業所租用的生產基地，施用肥料和栽培技術基本一致;甘肅蘭州百合的Ca、Fe含量最高;江西龍牙百合的Ca、Mg、Se、P含量最低;湖南龍牙百合的Fe、Cd、Pb含量最低。由此可見，不同地域來源藥食兼用百合樣品中礦物元素含量不但與其產地的土壤質地、土壤類型密切相關，此外還受品種和栽培措施等因素的影響。（2）本研究通過百合產地與其礦物元素含量之間的相關性研究，確定Ca、Mg、Zn、Cu、Fe為百合的特征元素。（3）由主成分分析和逐步判別分析得到的判別模型，可以明顯區分各產地的百合，并且可以區分百合的不同品種;同時，建立的方法對各產地及品種的百合可以進行溯源。

參考文獻

[1]趙祥云，等.百合[M].北京：中國農業出版社，2000.

[2]國家藥典委員會.中華人民共和國藥典：第一部[M].北京：中國醫學科學出版社，2010：11.

[3]楊云光，鄧成忠.食用百合品種介紹[J].中國果菜，2002（5）：32.

[4]李瑞琴，等.食用百合—土壤體系中鎘、鉛和汞的潛在生態和健康風險[J].食品科學，2016，37（5）：186-191.

[5]吳曉斌，等.百合皂苷的提取、純化及其對自由基的清除作用[J].天然產物研究與開發，2005，17（6）：777-780.

[6]周中流，石任兵，劉斌，等.卷丹甾體皂苷和酚類成分及其抗氧化活性研究[J].中草藥，2011，42（1）：21-24.

[7]李章念.兩種百合中黃酮類物質研究[D].陜西楊凌：西北農林科技大學，2007.

[8]付艷華，劉湘林，魏先運.食用百合引種鑒定和新品系選育研究[J].浙江農業科學，2010 （3）：468.

[9]李玉帆，明軍，王良桂，等.百合基本營養成分和活性物質研究進展[J].中國蔬菜，2012 （24）：7.

[10]孫紅梅，李天來，王化玲，等.不同溫度下貯藏的食用百合種球發根和完好率以及嗣后的植株生長[J].植物生理學通訊，2006，42（1）：47.

[11]王潤豐，等.野生卷丹百合黃酮類化合物抗氧化能力的研究[J].西北農業學報，2011，20（11）：152.

[12]黃淑珍.湖南香稻產地土壤特性與稻米品質的關系[J].湖南農業科學，1990（4）：37-40.

[13]吳繼紅，廖小軍，汪厚銀，等.綠茶礦物元素特征分析及產地判別研究[J].食品科學，2008，29（11）：494-497.

[14]童成英，何守陽，丁虎.茶葉產地與品質的元素、同位素鑒別技術研究進展[J].生態學雜志，2018，37（5）：1574-1583.

[15]何忠俊，等.三七主根穩定碳同位素組成與生態因子的關系[J].生態環境學報，2015，24（4）：561-568.

[16]李艷敏，等.主成分和判別分析在清香型白酒產地溯源中的應用[J].中國釀造，2018，37（1）：145-148.

[17]張玥，等.基于主成分分析和判別分析的大米產地溯源[J].中國糧油學報，2016，31（4）：1-5.

[18]興麗，趙鳳敏，曹有福，等.不同產地亞麻籽礦物質元素及脂肪酸組成的主成分分析[J].光譜學與光譜分析，2014，34（9）：2538-2543.

[19]程鐵轅，等.主成分分析法在濃香型白酒酒質評價中的應用研究[J].中國釀造，2011，30（1）：95-96.

[20]李昕，張明明.SPSS 22.0統計分析從入門到精通[M]. 北京：電子工業出版社，2015：280-312.

[21]謝漢龍，尚濤.SPSS統計分析與數據挖掘[M]. 北京：電子工業出版社，2012：484-488.

Abstract：The preliminary investigation was carried out to examine the feasibility of multi-element analysis in identifying the geographical origin of medicine and edible Lily，and to screen out the effective indicators in wheat origin assessment. The concentrations of 10 chemical elements including Ca、Mg、Zn、Cu、Fe、Se、K、P、Cd and Pb have been determined using atomic absorption spectrometer（iCE-3500）in 119 medicine and edible Lily samples from Gansu，Hunan，Jiangxi and Anhui，provinces of China，and analysis of variance （ANOVA），principal component analysis （PCA）and discriminant analysis （DA）were applied in data analysis.The results showed that the content of mineral elements in Lily samples from different geographical sources was not only closely related to soil texture and soil type，but also affected by varieties and cultivation measures.Ca，Mg，Zn，Cu and Fe were identified as characteristic elements of Lily.The established discriminant model could clearly distinguish lilies from different producing areas and varieties，and it provides scientific evidence for the origin traceability and identification of famous，special and new product.

Keywords：medicine and edible Lily;variety;geographical origin;multi-element

（責任編輯 唐建敏）