基于Sr穩(wěn)定同位素和礦物元素的鮑魚產(chǎn)地溯源技術(shù)研究

洪赫陽 張信澤 張秀珍 韓典峰 崔艷梅 趙軍強(qiáng) 宮向紅 徐英江

（1山東省海洋資源與環(huán)境研究院，山東省海洋生態(tài)修復(fù)重點(diǎn)實(shí)驗(yàn)室/煙臺市海珍品質(zhì)量安全控制與精深加工重點(diǎn)實(shí)驗(yàn)室，山東 煙臺 264006；2上海海洋大學(xué)食品學(xué)院，上海 201306）

近年來，為保護(hù)不同產(chǎn)地資源特色，我國實(shí)施了對產(chǎn)品地理品牌的建設(shè)。以鮑魚為例，自2006年以來，遼寧大連，山東長島、榮成及福建連江鮑魚相繼被列為地理標(biāo)志產(chǎn)品［1-4］。鮑魚具有極高的營養(yǎng)和藥用價值［5］，其主產(chǎn)區(qū)分布在遼寧、山東、福建和廣州等地［6］，而不同的地域特征使得鮑魚具有不同的品質(zhì)和營養(yǎng)價值，如南方鮑的平均硬度和咀嚼度較高、北方鮑的彈性高［7］。隨著人們對鮑魚需求量增多，以次充好等情況時有發(fā)生［8］，嚴(yán)重?fù)p害了消費(fèi)者的權(quán)益，因此，分析鮑魚地域差異及研發(fā)可靠有效的產(chǎn)地溯源技術(shù)具有重要意義。

穩(wěn)定同位素技術(shù)是國際上食品產(chǎn)地溯源與保護(hù)的前沿技術(shù)［9］。產(chǎn)地溯源研究中常用的穩(wěn)定同位素包括碳（C）、氮（N）、氫（H）、氧（O）、硫（S）、鉛（Pb）、鍶（Sr）和釹（Nd）等［10］，近年來，Sr 穩(wěn)定同位素因不受環(huán)境生物代謝、周圍環(huán)境及飲食過程影響的特點(diǎn)而在農(nóng)產(chǎn)品溯源研究中得到廣泛應(yīng)用［11］。相較于輕穩(wěn)定同位素，Sr穩(wěn)定同位素在氣候條件差異較小的鄰近區(qū)域產(chǎn)品判別時能發(fā)揮關(guān)鍵作用［12］。由于礦物元素組成和含量在水、土壤、飼料和空氣環(huán)境中分布不均，導(dǎo)致不同地域礦物元素指紋特征存在差異［13］，因此，礦物元素可作為農(nóng)產(chǎn)品產(chǎn)地溯源指標(biāo)。但陸源生物中元素易受生長環(huán)境、栽培管理措施等因素的影響，水產(chǎn)品中元素易受養(yǎng)殖飼料添加、水體環(huán)境污染等因素影響［14］，使元素指紋圖譜技術(shù)的應(yīng)用存在一定局限性，因此，元素結(jié)合Sr 穩(wěn)定同位素進(jìn)行產(chǎn)地溯源的研究逐漸增多。如Li等［15］利用87Sr/86Sr 結(jié)合28 種礦物元素區(qū)分貴州藍(lán)莓產(chǎn)地；康露等［16］利用87Sr/86Sr 和8 種礦質(zhì)元素區(qū)分新疆甜瓜產(chǎn)地；還有探究葡萄酒地理起源的研究［17］。上述研究均取得了理想的判別效果，但鮮見探究鮑魚產(chǎn)地追溯可行性的報道。

借助適當(dāng)?shù)幕瘜W(xué)計量學(xué)方法準(zhǔn)確地篩選有效指標(biāo)并提高鑒別的正確率是鮑魚溯源研究的關(guān)鍵［11］，如Zhou 等［18］在嬰兒奶粉溯源研究中指出，應(yīng)用人工神經(jīng)網(wǎng)絡(luò)（artificial neural network，ANN）化學(xué)計量方法建立的模型判別效果較線性判別分析（linear discriminant analysis，LDA）模型更好。本試驗(yàn)以連江、長島、榮成和大連4 個產(chǎn)地皺紋盤鮑（Haliotisdiscushannai）為研究對象，利用電感耦合等離子體質(zhì)譜儀（inductively coupled plasma mass spectrometry，ICP-MS）分析不同組織中的87Sr/86Sr、88Sr/86Sr 和20 種元素含量特征及差異，通過主成分分析（principal component analysis，PCA）、正交偏最小二乘法判別分析（orthogonal partial least squares discrimination analysis，OPLS-DA）、LDA 3 種化學(xué)計量學(xué)方法建立鮑魚產(chǎn)地溯源模型，旨在為鮑魚溯源體系建立提供參考。

1 材料與方法

1.1 材料與試劑

硝酸（優(yōu)級純），德國默克公司；鍶同位素標(biāo)準(zhǔn)物質(zhì)（NIST SRM 987），美國國家標(biāo)準(zhǔn)局；扇貝標(biāo)準(zhǔn)物質(zhì)（GBW10024），地球物理地球化學(xué)勘察研究所；多元素內(nèi)標(biāo)混合物（濃度100 μg·mL-1），美國安捷倫科技有限公司。

1.2 儀器與設(shè)備

ICAP RQ 電感耦合等離子體質(zhì)譜儀，美國Thermo Scientific 公司；ASX560 自動進(jìn)樣器，美國Teledyne CETAC公司；Topex-GT400全自動微波消解儀，上海屹堯儀器科技發(fā)展有限公司；Crymill冷凍研磨機(jī)，德國萊馳公司；Alpha 1-4LDplus冷凍干燥機(jī)，德國Chirst公司。

1.3 樣品采集和制備

鮑魚樣品采集自福建省連江市、山東省（長島縣、榮成市）和遼寧省大連市3個省份的4個主要產(chǎn)區(qū)，不同產(chǎn)區(qū)地理氣候環(huán)境具有一定差異，且連江具有典型的亞熱帶氣候特征［19］。所有鮑魚樣品均于2021年9、10月采集，每個產(chǎn)區(qū)各選取大小相似的50 只鮑魚，每只鮑魚樣品濕重約57.4～63.9 g。鮑魚用超純水清洗，對整只鮑魚取殼肌、裙邊、消化腺性腺及外套膜4 個組織，置于-80 ℃冰箱預(yù)冷凍24 h，再放入干燥機(jī)干燥48 h，取出研磨過篩（100目）后放入干燥器中保存待測。

稱取樣品（0.2±0.005）g 置于聚四氟乙烯消解罐中，加入10 mL 硝酸，加蓋過夜，再放入微波消解儀中，具體條件見表1。冷卻后將消解罐放在180 ℃控溫電熱板上，趕酸至體積小于1 mL，后定容至25 mL，混勻備用，同時做空白試驗(yàn)。

1.4 測定項目與方法

1.4.1 Sr 穩(wěn)定同位素的測定 使用ICP-MS 對樣品中87Sr/86Sr、88Sr/86Sr 進(jìn)行測定。配制80 μg·L-1Sr 標(biāo)準(zhǔn)溶液，試驗(yàn)采用插入監(jiān)控物質(zhì)，即每批樣品測試時加入空白樣品、重復(fù)樣品和標(biāo)準(zhǔn)物質(zhì)［20］。Sr 穩(wěn)定同位素標(biāo)準(zhǔn)值為87Sr/86Sr=0.710 34±0.000 26、88Sr/86Sr=8.376 1±0.003 25，儀器精密度為0.35%和0.37%，表明測量結(jié)果具有良好的再現(xiàn)性。

1.4.2 元素含量的測定 對于元素同位素的選擇，既要選擇靈敏度和豐度高的同位素，又要避開同量異位素的干擾［21］，如64Ni+和64Zn+［22］。本試驗(yàn)測定了20 種元素，其同位素分別為11B、27Al、48Ti、51V、55Mn、59Co、63Cu、66Zn、71Ga、75As、77Se、88Sr、111Cd、137Ba、208Pb、7Li、23Na、24Mg、39K、44Ca。使用多元素混合標(biāo)準(zhǔn)溶液繪制標(biāo)準(zhǔn)曲線，采用1 μg·mL-1內(nèi)標(biāo)混合溶液（Se、Ge、Rh、In）校正儀器的漂移。儀器精密度（1.2%～6.4%）良好，該方法可用于鮑魚元素的測定。

1.5 儀器工作條件

ICP-MS儀器開機(jī)預(yù)熱至穩(wěn)定后，用調(diào)諧液優(yōu)化儀器參數(shù)，獲得最佳工作條件。ICP-MS工作參數(shù)如下：發(fā)生器功率1 550 W；冷卻氣流量14.0 L·min-1；霧化氣流量1.0 mL·min-1；氦氣流量4.125 mL·min-1；霧化室溫度2 ℃；采樣深度5 mm；駐留時間0.02 s；重復(fù)次數(shù)3次；氧化物（140Ce16O/140Ce）≤2%；雙電荷（137Ba++）/（137Ba+）≤3%。在He模式下進(jìn)行樣品測定，并在線引入1 μg·mL-1的內(nèi)標(biāo)溶液。

1.6 數(shù)據(jù)處理

采用SPSS 26.0 軟件進(jìn)行單因素方差分析（One-Way ANOVA）及事后最小顯著差異法（least significant difference，LSD）多重比較、PCA和LDA，采用Origin 2022b軟件繪制雷達(dá)圖和得分散點(diǎn)圖，采用SIMCA 軟件進(jìn)行OPLS-DA及圖形繪制，用Excel 2021進(jìn)行數(shù)據(jù)處理。

2 結(jié)果與分析

2.1 鮑魚Sr穩(wěn)定同位素比值及礦物元素的特征差異分析

表2 顯示了鮑魚不同指標(biāo)指紋分布特征差異。經(jīng)過LSD多重比較分析可知，4個組織中不同產(chǎn)地的87Sr/86Sr、88Sr/86Sr比值整體存在顯著性差異（P＜0.05）。從組織來看，鮑魚殼肌中87Sr/86Sr比值分布較裙邊、消化腺性腺和外套膜偏正，而88Sr/86Sr比值分布較其他組織偏負(fù)，表明殼肌中87Sr/86Sr、88Sr/86Sr指標(biāo)展現(xiàn)了明顯的地域差異，與其他組織存在較大的分餾效應(yīng)。從產(chǎn)地來看，除殼肌指標(biāo)外，大連和榮成產(chǎn)區(qū)樣本87Sr/86Sr、88Sr/86Sr值普遍高于連江和長島，多數(shù)組織樣本呈大連＞榮成＞長島＞連江的趨勢，這與緯度變化一致。

表2 鮑魚中Sr穩(wěn)定同位素比值和20種元素含量分布情況Table 2 Distribution of Sr stable isotope ratio and 20 elementsin abalone/（mg·kg-1）

4 個組織中不同產(chǎn)地20 種礦物元素（除消化腺性腺中Li元素）含量整體存在顯著性差異（P＜0.05）。連江鮑魚樣品Sr、Pb、Cd元素含量相對偏高，長島鮑魚的4個組織中Cd元素含量均相對偏高，榮成鮑魚As元素含量和大連鮑魚B 元素含量均較其他產(chǎn)區(qū)偏高，且主要富集在殼肌中。整體來說，多數(shù)元素含量的產(chǎn)地特征存在相似性，即連江產(chǎn)地鮑魚樣品的元素含量相對偏高，而大連樣品相對偏低，且組織間特征也存在相似性，即消化腺性腺的元素含量普遍高于殼肌、裙邊和外套膜。根據(jù)表2 差異顯著性分析可知，殼肌B、Al、Se、Ba、Li、K、Ca元素，裙邊Mg元素，消化腺性腺Ti、Zn、K、Ca 元素在不同地域間展示了良好的地域差異性（P＜0.05）。

2.2 基于Sr穩(wěn)定同位素和礦物元素的主成分分析（PCA）

對4個產(chǎn)地鮑魚樣品進(jìn)行PCA 模型建立。結(jié)果顯示，不同組織的檢驗(yàn)統(tǒng)計量（kaiser-meyer-olkin，KMO）值均大于0.6，Bartlett 的球形度檢驗(yàn)顯著性均小于0.05，表明PCA分析可行。如表3所示，不同組織前3個主成分的累計方差貢獻(xiàn)率在66%以上，可充分反映原始數(shù)據(jù)信息。利用其載荷矩陣絕對值做雷達(dá)圖（圖1），進(jìn)一步篩選特征指標(biāo)，即載荷絕對值越大，距離中心點(diǎn)越遠(yuǎn)，貢獻(xiàn)率越大。結(jié)果表明，殼肌Ca、Sr、Mg、Mn、V、Li、Ti、Co、Se、Ba、K、87Sr/86Sr，裙邊Ti、Al、Mn、Ca、Ga、Sr、Li、Ba、V、K、Cd，消化腺性腺Ti、Ca、B、Na、Ga、Pb、Mg、Mn、Sr、Al、Cd，外套膜Sr、Mn、As、V、Ca、Se、Co、K、Ti、Mg、Zn、Pb、Cd 的載荷絕對值在0.7～1.0 之間，可基本反映主成分的信息。再繪制得分散點(diǎn)圖（圖2），可見不同組織PCA 模型可在一定程度上將4 個產(chǎn)地的鮑魚區(qū)分開來，但不同產(chǎn)地樣品間依然存在重疊現(xiàn)象，因此采用判別分析進(jìn)一步對4 個產(chǎn)地鮑魚樣品進(jìn)行鑒別。

圖1 不同組織鮑魚前3個主成分載荷絕對值雷達(dá)圖Fig.1 Radar chart of absolute values of the first three principal component loads in different tissues of abalone

圖2 不同組織鮑魚的主成分因子得分圖Fig.2 Principal component factor scores of abalone in different tissues

表3 不同組織鮑魚Sr穩(wěn)定同位素和礦物元素的PCA結(jié)果Table 3 Principal component analysis results of Sr stable isotopes and multi-element from different tissues in abalone

2.3 基于Sr穩(wěn)定同位素和礦物元素的判別分析

2.3.1 正交偏最小二乘法判別分析（OPLS-DA） 本研究以22 個指標(biāo)作為因變量，不同產(chǎn)地作為自變量，構(gòu)建溯源模型。不同組織自變量（R2x）、因變量（R2y）擬合指數(shù)和模型預(yù)測指數(shù)（Q2）均大于0.5，表示模型擬合結(jié)果可接受［23］。

如OPLS-DA 得分圖（圖3）所示，不同組織具體判別結(jié)果如下：殼肌和裙邊均可將連江、榮成和大連鮑魚明顯區(qū)分開來，消化腺性腺可顯著區(qū)分連江和大連樣品，外套膜中心點(diǎn)處得分存在重疊現(xiàn)象，判別效果略差。根據(jù)變量重要性投影（variable importance in the projection，VIP）大于1 確定的模型顯著變量，殼肌、裙邊、消化腺性腺、外套膜中顯著變量指標(biāo)分別為88Sr/86Sr、87Sr/86Sr、Cd、K、V、Zn、Cu、Se；B、Cd、As、Al、Ba、Ga、Se；88Sr/86Sr、Cd、B、Na、Ca、Se、Ba、Co、As；88Sr/86Sr、B、Li、Ga、Al、Cu、Cd、Zn。

圖3 不同組織鮑魚OPLS-DA得分圖（左）和VIP貢獻(xiàn)圖（右）Fig.3 Discriminatory score graph （left） and VIP contribution graph （right） of OPLS-DA of different tissues in abalone

對4個產(chǎn)地鮑魚模型進(jìn)行判別分析，結(jié)果如表4所示。殼肌、裙邊、消化腺性腺、外套膜總體正確判別率分別為99%、97%、99%和95%，其中連江和榮成消化腺性腺指標(biāo)均可100%正確識別樣品產(chǎn)地，長島和大連正確判別率均為98%，溯源效果最優(yōu)。因此，OPLSDA 可將4 個產(chǎn)地鮑魚樣品準(zhǔn)確地區(qū)分開來，展現(xiàn)了OPLS-DA較好的判別能力。

表4 不同組織鮑魚OPLS-DA判別結(jié)果Table 4 Discriminant results of OPLS-DA of different tissues in abalone

2.3.2 線性判別分析（LDA） 利用Fisher 判別函數(shù)對鮑魚不同產(chǎn)地樣品進(jìn)行線性判別分析，并結(jié)合留一交叉驗(yàn)證法（Leave-One-Out Cross Validation）對判別結(jié)果進(jìn)行準(zhǔn)確性評估。結(jié)果顯示，單獨(dú)以87Sr/86Sr和88Sr/86Sr 為分析指標(biāo)，鮑魚4 種組織正確判別率較低，均小于90.0%。故針對此問題，進(jìn)一步探究Sr 穩(wěn)定同位素與礦物元素指標(biāo)結(jié)合的溯源能力。

為了提高實(shí)際檢測效率，采用逐步判別分析法進(jìn)行分析。通過威爾克Lambda 分析可知，4 種組織均得出了3 個典型判別函數(shù)（α=0.01），其中，前2 個判別函數(shù)的累積貢獻(xiàn)率在76.8%以上，以此做得分散點(diǎn)圖，如圖4-A～D 所示，殼肌和裙邊可將連江、長島和榮成鮑魚樣品有效區(qū)分開來；消化腺性腺可將大連鮑魚與長島榮成兩地鮑魚區(qū)分開，但榮成部分樣品與長島分?jǐn)?shù)重疊；外套膜中不同產(chǎn)地樣品均出現(xiàn)重疊現(xiàn)象。

對4 個產(chǎn)地的鮑魚樣品進(jìn)行歸類（表5），結(jié)果如下：殼肌B、Al、Ti、As、Se、Sr、Cd、Ba、Li、Na、Mg、K、87Sr/86Sr，裙邊B、Ti、Co、Zn、As、Se、Sr、Cd、Ba、Na、Mg、Ca、88Sr/86Sr，消化腺性腺V、Mn、Cu、Zn、Ga、As、Se、Sr、Cd、Pb、Mg、K、Ca、87Sr/86Sr，外套膜B、Ti、V、Cu、Ga、Se、Sr、Li、Na、Mg、Ca、87Sr/86Sr進(jìn)入了判別模型，其交叉驗(yàn)證正確判別率分別為97.0%、97.5%、99.0%和96.0%，正確率均在96%以上，效果良好。與單獨(dú)使用Sr 穩(wěn)定同位素相比，結(jié)合判別使4 種組織正確率均明顯提高，尤其是殼肌指標(biāo)的正確率由56%提升至97%，說明結(jié)合判別溯源效果優(yōu)于單一判別。此外，如圖4-A～C 和表5 所示，肌肉組織（殼肌、裙邊）可區(qū)分長島和榮成樣本，消化腺性腺可區(qū)分大連與長島樣本，二者結(jié)合建模可成功將所有產(chǎn)地樣品區(qū)分開（圖4-E、F），100%精準(zhǔn)識別4個產(chǎn)地的樣本。

表5 不同組織鮑魚LDA結(jié)果aTable 5 Fisher linear discriminant analysis results of different tissues in abalone

3 討論

3.1 不同產(chǎn)地對Sr 穩(wěn)定同位素比值及礦物元素的指紋分布影響

Sr有4種天然存在的同位素，即84Sr、86Sr、87Sr和88Sr，一般認(rèn)為87Sr 是可變的，84Sr、86Sr 和88Sr 之間的比值為常數(shù)［24］。河流的87Sr/86Sr 取決于平均基巖年齡和成分流域［25］，在食品溯源中的優(yōu)勢在于，在加工之后直至最終產(chǎn)物，Sr 始終保持其原始同位素比率不變。本研究中，我國不同地區(qū)的87Sr/86Sr 范圍較Gregor?i? 等［26］對斯洛文尼亞地域牛奶溯源研究的比值范圍偏正，表明87Sr/86Sr 在不同國家有較大的地質(zhì)差異。前人研究表明，地質(zhì)體在化學(xué)風(fēng)化和碳酸鹽結(jié)核等過程中，其88Sr/86Sr 也可能發(fā)生顯著變化［27-28］。本研究同樣發(fā)現(xiàn)，88Sr/86Sr 在不同產(chǎn)地分布范圍有較顯著差異，說明88Sr/86Sr是可變的，兩種Sr比值特征分布均表現(xiàn)為大連樣品比值相對偏高，連江樣品比值相對偏低，究其原因，一是風(fēng)化期間的分餾效應(yīng)會導(dǎo)致同位素Sr 的損失［27］，一般情況下，在濕熱地區(qū)化學(xué)風(fēng)化作用占主導(dǎo)地位［29］，使連江樣品比值偏低；二是Sr 穩(wěn)定同位素與地質(zhì)特征和地理氣候環(huán)境相關(guān)。本試驗(yàn)中4 個產(chǎn)地緯度值與Sr 穩(wěn)定同位素特征分布一致。此外，本研究發(fā)現(xiàn)，殼肌中87Sr/86Sr 存在較大的地域差異，由于肌肉的細(xì)胞更新速度較慢，穩(wěn)定同位素特征達(dá)到穩(wěn)定狀態(tài)的平衡時間較長，更容易反映時間尺度上的變化情況［30-31］，因此，殼肌指標(biāo)更能反映長期鮑魚產(chǎn)地差異。

生物體不能自行合成礦物質(zhì)元素［32］，因此其元素特征的差異反映了周圍環(huán)境（即水和沉積物）和飲食的化學(xué)成分［33］。由于海水中重金屬會在食物網(wǎng)中累積，導(dǎo)致不同地區(qū)生物體中重金屬含量產(chǎn)生較大差異，重金屬含量偏高的現(xiàn)象與人類活動息息相關(guān)［34］。福建連江養(yǎng)殖密度過大、生活垃圾、工業(yè)排污等水體污染都會造成水體富營養(yǎng)化［35］。本研究連江鮑魚中尤其以Sr、Pb 和Cd元素含量較高，且富集在不同組織中，這是因?yàn)樵氐拇x在不同組織和器官中存在差異，攝入的元素首先經(jīng)過消化道，再到達(dá)二次積累器官［36］。工業(yè)廢水是海水主要污染途徑，其排放會使近海海水和浮游生物體內(nèi)的Cd含量不斷增高［37］。由此推測，本研究長島鮑魚4個組織中Cd元素含量均相對偏高與工業(yè)廢水的排放相關(guān)。有研究表明，野生捕撈魚類較養(yǎng)殖的As、Se、K濃度較高［38］，這與本研究中榮成鮑魚指紋特征一致，推測三種指標(biāo)可在一定程度上反映鮑魚原產(chǎn)地信息。

3.2 影響鮑魚不同產(chǎn)地判別率的因素

不同海域根據(jù)地理條件不同，使Sr 穩(wěn)定同位素均具有產(chǎn)地指紋特征。連江和榮成分別位于臺灣海峽和黃海海域，其地理氣候特征具有顯著差異［19］，緯度也相差較大，這使得連江和榮成多數(shù)指標(biāo)產(chǎn)地判別良好。渤海海域因三面環(huán)陸的特殊地理構(gòu)造，使得海水交換能力受到限制，海洋生態(tài)系統(tǒng)易受到沿岸人類活動的干擾［39］。大連和長島同處于渤海海域，地域較近，導(dǎo)致本研究單獨(dú)使用一種組織建模難以將4個產(chǎn)地鮑魚區(qū)分開，小空間尺度難以對樣品100%判別［16］。OPLSDA 和LDA 結(jié)果表明，消化腺性腺單獨(dú)判別時，均體現(xiàn)了較好的產(chǎn)地區(qū)分效果，但消化腺性腺只能反映鮑魚短期內(nèi)溯源指標(biāo)信息［29］，因此可將肌肉組織和消化腺性腺結(jié)合溯源，肌肉組織代表了長期元素富集情況，而具有高代謝活性的消化腺性腺主要傾向于反映近期的元素富集情況［40］，二者結(jié)合能100%區(qū)分4個產(chǎn)地鮑魚樣品。兩種組織的聯(lián)合鑒別為距離較小的產(chǎn)地溯源研究提供了方法參考。

雖然單獨(dú)使用Sr 穩(wěn)定同位素對鮑魚溯源的表現(xiàn)并不樂觀，但將其作為輔助指標(biāo)仍有加成的效果，如Won等［41］認(rèn)為多重穩(wěn)定同位素分析可有效應(yīng)用于菲律賓蛤仔的產(chǎn)地追溯。李啟航［20］在藍(lán)莓產(chǎn)地溯源研究中指出，Sr穩(wěn)定同位素結(jié)合元素建模后判別率由81.3%提升至93.8%。本試驗(yàn)將Sr穩(wěn)定同位素與20種元素結(jié)合對鮑魚產(chǎn)地溯源，其判別效果最優(yōu)，高于單獨(dú)使用Sr穩(wěn)定同位素溯源，表明Sr 穩(wěn)定同位素和元素結(jié)合判別是鮑魚產(chǎn)地溯源的有效方法。

4 結(jié)論

本研究采用Sr 穩(wěn)定同位素與礦物元素聯(lián)合的分析方法對鮑魚的產(chǎn)地進(jìn)行追溯。結(jié)果表明，不同產(chǎn)地指標(biāo)（除消化腺性腺中Li 元素）整體差異顯著，具有較明顯的地理特征。采用PCA、OPLS-DA和LDA對不同產(chǎn)地鮑魚建模，結(jié)果表明，殼肌、裙邊、消化腺性腺和外套膜可基本實(shí)現(xiàn)鮑魚4 個產(chǎn)地的判別，其中無監(jiān)督的PCA方法無法完全區(qū)分4個產(chǎn)地鮑魚樣品，但證明了指紋圖譜技術(shù)溯源的可行性；OPLS-DA模型溯源的正確判別率分別為99.0%、97.0%、99.0%和95.5%，LDA模型溯源交叉驗(yàn)證正確率分別為97%、97.5%、99.0%和96.0%，兩種化學(xué)計量學(xué)方法均可有效鑒別鮑魚產(chǎn)地，但OPLS-DA模型在肌肉組織的判別中效果較好。此外，肌肉組織（殼肌或裙邊）和消化腺性腺的溯源效果互補(bǔ)，二者指標(biāo)結(jié)合溯源可對4個產(chǎn)地鮑魚100%精準(zhǔn)判別。