礦物元素結(jié)合穩(wěn)定同位素的中華絨螯蟹產(chǎn)地溯源

張政權(quán)，黃冬梅*，蔡友瓊，孟憲菁，史永富，孔 聰，黃宣運(yùn)，湯云瑜，張 璇，楊光昕

（1.上海海洋大學(xué)食品學(xué)院，上海 201306；2.中國(guó)水產(chǎn)科學(xué)研究院東海水產(chǎn)研究所，農(nóng)業(yè)農(nóng)村部遠(yuǎn)洋與極地漁業(yè)創(chuàng)新重點(diǎn)實(shí)驗(yàn)室，上海 200090；3.賽默飛世爾科技（中國(guó)）有限公司，上海 201206）

中華絨螯蟹又稱河蟹、大閘蟹，屬于日本絨螯蟹中華亞種，因其風(fēng)味獨(dú)特且營(yíng)養(yǎng)豐富，在我國(guó)深受廣大消費(fèi)者的歡迎。自二十世紀(jì)七十年代中期以來(lái)，中華絨螯蟹在我國(guó)的養(yǎng)殖產(chǎn)量迅速增長(zhǎng)，目前已經(jīng)成為我國(guó)最主要的水產(chǎn)品之一。中華絨螯蟹的養(yǎng)殖主要集中在遼河、黃河和長(zhǎng)江三大水系中[1]，據(jù)消費(fèi)者多年的反饋，長(zhǎng)江水系養(yǎng)殖的河蟹無(wú)論在口感風(fēng)味還是個(gè)體大小上都相對(duì)較好，尤其是以陽(yáng)澄湖出產(chǎn)的大閘蟹最負(fù)盛名，每年大閘蟹上市的季節(jié)，陽(yáng)澄湖大閘蟹都是消費(fèi)者爭(zhēng)相品嘗的對(duì)象。而現(xiàn)在無(wú)論網(wǎng)絡(luò)上還是水產(chǎn)市場(chǎng)上都能夠看到“原產(chǎn)自”陽(yáng)澄湖的大閘蟹，其中大多數(shù)無(wú)疑為不良商家為了獲取更高利潤(rùn)假冒陽(yáng)澄湖大閘蟹進(jìn)行銷售，這不僅損害了消費(fèi)者和當(dāng)?shù)仞B(yǎng)殖戶的利益，同時(shí)還對(duì)地區(qū)品牌構(gòu)成了嚴(yán)重威脅。為了規(guī)范水產(chǎn)養(yǎng)殖市場(chǎng)秩序，保護(hù)消費(fèi)者和養(yǎng)殖戶的合法權(quán)益，建立完善的河蟹防偽判別體系，追溯假冒偽劣產(chǎn)品源頭至關(guān)重要。

傳統(tǒng)的水產(chǎn)品產(chǎn)地溯源方法是添加物理防偽標(biāo)簽[2]，這種方式一定程度上能夠做到對(duì)不同來(lái)源的產(chǎn)品進(jìn)行區(qū)分，但也存在一些不足，這些獨(dú)立于產(chǎn)品之外的防偽標(biāo)識(shí)在加工運(yùn)輸過(guò)程中容易造成標(biāo)簽的脫落或損壞導(dǎo)致信息鏈中斷，同時(shí)虛假標(biāo)簽的使用更是會(huì)讓防偽標(biāo)簽形同虛設(shè)。因此，急需開發(fā)一種新的水產(chǎn)品產(chǎn)地溯源方式，能夠反映水產(chǎn)品產(chǎn)地的環(huán)境特征和水產(chǎn)品自身的飲食特征，而礦物元素和穩(wěn)定同位素豐度是目前最能反映這些特征的因素。目前國(guó)內(nèi)外采用礦物元素和穩(wěn)定同位素技術(shù)對(duì)茶葉[3-4]、稻米[5-8]、豆類[9-10]、牛羊肉[11-17]等農(nóng)產(chǎn)品的產(chǎn)地溯源已有一些研究報(bào)道，對(duì)中華絨螯蟹進(jìn)行產(chǎn)地溯源的研究，主要集中在對(duì)其外部形態(tài)和風(fēng)味特征方面差異性的比較[18-21]，這種研究方式存在很大的局限性，且含主觀因素較多，對(duì)外部形態(tài)相近的水產(chǎn)品很難進(jìn)行區(qū)分，也有學(xué)者通過(guò)脂肪酸含量[22]、礦物元素含量[23-24]等方面對(duì)不同地區(qū)的中華絨螯蟹進(jìn)行了區(qū)分，但僅局限于采用單一技術(shù)對(duì)中華絨螯蟹溯源的研究，暫未發(fā)現(xiàn)采用礦物元素與穩(wěn)定同位素結(jié)合的方式對(duì)長(zhǎng)江水系的中華絨螯蟹進(jìn)行產(chǎn)地溯源。本研究采用礦物元素含量和穩(wěn)定同位素豐度結(jié)合的方式，比較長(zhǎng)江水系中陽(yáng)澄湖、崇明和興化3 個(gè)地區(qū)河蟹的8 種礦物元素的含量和C、N穩(wěn)定同位素豐度，利用SPSS 19.0數(shù)據(jù)處理軟件對(duì)所得數(shù)據(jù)進(jìn)行方差分析和顯著性比較，利用聚類分析對(duì)3 個(gè)地區(qū)的樣品進(jìn)行了初步聚類，利用線性判別分析對(duì)樣品進(jìn)一步歸類比較，并建立產(chǎn)地判別模型，然后利用主成分分析篩選出能夠區(qū)分3 個(gè)地區(qū)河蟹的特征指標(biāo)，以期實(shí)現(xiàn)對(duì)中華絨螯蟹的產(chǎn)地溯源。

1 材料與方法

1.1 材料與試劑

河蟹樣品由養(yǎng)殖戶協(xié)助于2019年7月分別于陽(yáng)澄湖中、興化養(yǎng)殖塘和崇明養(yǎng)殖塘捕撈，雌蟹和雄蟹各10 只，質(zhì)量范圍在63～107 g之間，每個(gè)樣品分別用自封袋封裝編號(hào)并立即運(yùn)回實(shí)驗(yàn)室（表1）。取其腿部及腹部肌肉混合均勻，將樣品放入真空冷凍干燥機(jī)連續(xù)干燥至恒質(zhì)量后，研磨過(guò)100 目篩，放入干燥器中保存?zhèn)溆谩?/p>

表1 河蟹樣品基本信息Table 1 Information about crab samples collected for this study and their geographical origins

硝酸（分析純） 江蘇晶銳化學(xué)股份有限公司；超純水、過(guò)氧化氫（分析純） 上海安譜實(shí)驗(yàn)科技股份有限公司；1 mg/mL K、Ca、Na、Mg、Mn、Zn、Cu、Fe元素單標(biāo)準(zhǔn)液 國(guó)家標(biāo)準(zhǔn)物質(zhì)認(rèn)證中心；穩(wěn)定同位素標(biāo)準(zhǔn)品USGS40 東莞市譜標(biāo)實(shí)驗(yàn)器材科技有限公司。

1.2 儀器與設(shè)備

Milli-Q ELEMENT超純水機(jī) 美國(guó)默克密理博 公司；FreeZone真空冷凍干燥機(jī) 美國(guó)Labconco公司；JA5000C精密電子天平 常州市幸運(yùn)電子天平有限 公司；ETHOS 1微波消解儀 意大利邁爾斯通公司；HP-H35SC電加熱板 盛達(dá)杰森自動(dòng)化設(shè)備有限公司；ContrAA-700原子吸收光譜儀 德國(guó)耶拿分析儀器股份公司；EA IsoLink元素分析儀、DELTA V Advantage同位素比率質(zhì)譜儀 美國(guó)Thermo Fisher公司。

1.3 方法

1.3.1 元素測(cè)定

準(zhǔn)確稱取0.5 g樣品于50 mL尖底離心管，加入5 mL硝酸和1 mL過(guò)氧化氫，放置于通風(fēng)櫥酸解過(guò)夜，第2天將酸解后的樣品轉(zhuǎn)移至微波消解罐中放入微波消解儀中消解。將消解后的樣品放置加熱板上趕酸至溶液呈澄清并有白色煙霧冒出，取出放置至室溫，定容稀釋上機(jī)。樣品液經(jīng)進(jìn)樣管進(jìn)入霧化器，經(jīng)霧化成為水蒸氣進(jìn)入原子化器，水蒸氣在原子化器內(nèi)被空氣-乙炔混合氣體燃燒產(chǎn)生的高溫原子化，同時(shí)吸收光源產(chǎn)生的電磁輻射，電磁輻射能的吸收與蒸汽中被測(cè)元素的含量呈正比，通過(guò)最終測(cè)定的吸光度可以計(jì)算出待測(cè)元素的含量。本研究測(cè)定8 種礦物元素的標(biāo)準(zhǔn)曲線線性相關(guān)系數(shù)均在0.999以上，方法的回收率為96.0%～106.9%，精密度為1.1%～5.4%，各元素的方法檢出限分別為K 118.0 mg/kg、Ca 19.0 mg/kg、 Na 101.0 mg/kg、Mg 89.0 mg/kg、Mn 0.06mg/kg、Zn 0.70 mg/kg、Fe 0.44 mg/kg、Cu 0.10 mg/kg。

1.3.2 穩(wěn)定同位素測(cè)定

用百萬(wàn)分之一天平準(zhǔn)確稱取400～450 μg蟹肉干粉包裹在錫紙杯里，依次放入進(jìn)樣盤中，樣品通過(guò)自動(dòng)進(jìn)樣器進(jìn)入EA IsoLink元素分析單元，在元素分析單元中燃燒并被還原成純凈的CO2和N2，經(jīng)過(guò)ConFlo IV連續(xù)流接口裝置經(jīng)氦氣稀釋后進(jìn)入DELTA V Advantage穩(wěn)定同位素質(zhì)譜儀進(jìn)行檢測(cè)。儀器條件如下：載氣流速180 mL/min，進(jìn)樣口吹掃速率70 mL/min，氧氣流速175 mL/min，燃燒爐溫度980 ℃，還原爐溫度55 ℃，每測(cè)定20 個(gè)樣品插入2 個(gè)USGS40（咖啡因）標(biāo)準(zhǔn)品進(jìn)行質(zhì)量控制。CN穩(wěn)定同位素的測(cè)試精度均優(yōu)于0.1‰。

穩(wěn)定性C、N同位素豐度分別表示為δ13C、δ15N，其中δ13C的相對(duì)標(biāo)準(zhǔn)為V-PDB（維也納-PeeDee箭石標(biāo)準(zhǔn)），δ15N的相對(duì)標(biāo)準(zhǔn)為空氣中的N2。計(jì)算公式如下：

式中：R為重同位素與輕同位素的豐度比，即13C/12C、15N/14N。

1.4 數(shù)據(jù)處理軟件

采用SPSS 19.0數(shù)據(jù)處理軟件對(duì)實(shí)驗(yàn)結(jié)果進(jìn)行方差分析、顯著性比較、系統(tǒng)聚類分析、線性判別分析（linear discriminant analysis，LDA）和主成分分析（principal component analysis，PCA），采用SPSS 19.0作樣品分布散點(diǎn)圖。

2 結(jié)果與分析

2.1 元素含量和穩(wěn)定同位素豐度產(chǎn)地差異分析

如表2、3所示，8 種元素和2 種穩(wěn)定同位素在3 個(gè)地區(qū)的樣品中呈現(xiàn)不同程度的差異。單因素方差分析和Duncan多重比較分析結(jié)果顯示K、Mg、Cu三種元素在3 個(gè)地區(qū)間均存在顯著差異（P＜0.05），且K和Cu在崇明地區(qū)的樣品中含量最高，Mg在陽(yáng)澄湖地區(qū)的樣品中含量最高。Na在崇明和興化地區(qū)的樣品中含量差異不顯著 （P＞0.05），但Na在兩地區(qū)樣品中的含量均明顯高于陽(yáng)澄湖地區(qū)（P＜0.05）。Ca和Zn在崇明和興化2 個(gè)地區(qū)樣品中的含量差異不顯著（P＞0.05），且均顯著低于陽(yáng)澄湖地區(qū)（P＜0.05）。而Fe、Mn和δ15N在崇明和陽(yáng)澄湖2 個(gè)地區(qū)樣品中的含量差異不顯著（P＞0.05），但均顯著高于興化地區(qū)（P＜0.05）。δ13C在陽(yáng)澄湖與興化兩地區(qū)樣品中差異不顯著（P＞0.05），且都顯著高于崇明地區(qū)（P＜0.05）。鑒于三地區(qū)樣品中8 種礦物元素含量和C、N穩(wěn)定同位素豐度的差異，表明利用這種差異區(qū)別三地區(qū)的樣品具有一定的可行性。

表2 礦物元素含量（干質(zhì)量）測(cè)定結(jié)果Table 2 Mineral element contents of crab samples mg/kg

表3 穩(wěn)定同位素豐度測(cè)定結(jié)果Table 3 Stable isotope abundance in crab samples‰

2.2 聚類分析

圖1 3 個(gè)地區(qū)樣品的系統(tǒng)聚類結(jié)果Fig. 1 Systematic clustering resultsof three regional samples

為觀察礦物元素與穩(wěn)定同位素結(jié)合在對(duì)三地區(qū)樣品歸類方面的能力，將這8 種元素和C、N穩(wěn)定同位素豐度作為共同變量采用系統(tǒng)聚類方式對(duì)三地區(qū)的60 個(gè)樣品進(jìn)行聚類，聚類結(jié)果（圖1）顯示，在歐式距離為6處將樹狀圖切斷，可將樣品大致分為兩大類，第1大類22 個(gè)樣品中有16 個(gè)為興化地區(qū)樣品，6 個(gè)為陽(yáng)澄湖地區(qū)樣品，興化地區(qū)樣品占總數(shù)的72.7%，占興化地區(qū)樣品總量的80.0%，第2大類30 個(gè)樣品中有16 個(gè)為崇明地區(qū)樣品，占樣品總量的53.0%，占崇明地區(qū)樣品總量的80%，11 個(gè)為陽(yáng)澄湖地區(qū)樣品，占樣品總量的36.7%，占陽(yáng)澄湖樣品總量的55.0%，3 個(gè)為興化地區(qū)樣品，另外有1 個(gè)興化樣品、3 個(gè)陽(yáng)澄湖地區(qū)樣品和4 個(gè)崇明地區(qū)樣品被錯(cuò)誤分成了3 小類。聚類分析結(jié)果總體上將興化地區(qū)樣品單獨(dú)聚為一類，將陽(yáng)澄湖和崇明地區(qū)樣品共同聚為一類，顯示綜合礦物元素和穩(wěn)定同位素在對(duì)3 個(gè)地區(qū)的樣品進(jìn)行探索性聚類方面表現(xiàn)出一定的能力，但是聚類效果還不夠理想，存在較多的誤判樣品。

2.3 LDA結(jié)果

LDA是一種有監(jiān)督的數(shù)據(jù)歸類方式，可以在假設(shè)樣品分類的基礎(chǔ)上把已知樣品作為訓(xùn)練集，利用各類樣品之間的相關(guān)性建立判別模型，從而可以利用判別模型對(duì)未知樣品進(jìn)行判別歸類，還可以采用“留一法”交叉驗(yàn)證把訓(xùn)練集中的除某一個(gè)樣品之外的其余樣品作為訓(xùn)練集，建立相關(guān)的判別函數(shù)對(duì)該樣品進(jìn)行判別，在這個(gè)過(guò)程中，總樣品中的每一個(gè)個(gè)體都被除該個(gè)體之外的樣品建立的判別函數(shù)進(jìn)行驗(yàn)證，從而極大地檢驗(yàn)了判別模型的穩(wěn)定性。本研究采用步進(jìn)式DA方法，按照每個(gè)變量在判別分析過(guò)程中的重要性由大到小逐步選為判別變量，建立分類判別函數(shù)。如表4所示，崇明和興化地區(qū)的樣品實(shí)現(xiàn)了全部100%的正確分類，陽(yáng)澄湖地區(qū)的樣品有1 個(gè)被誤判為崇明地區(qū)，正確判別率為95.0%，三地區(qū)樣品的總體正確判別率為98.3%，交叉驗(yàn)證結(jié)果中陽(yáng)澄湖地區(qū)的1 個(gè)樣品被誤判到崇明地區(qū)，雖然3 個(gè)地區(qū)中有個(gè)別樣品出現(xiàn)了誤判，其總體的交叉驗(yàn)證判別正確率仍然高達(dá)98.3 %。

表4 LDA及交叉驗(yàn)證結(jié)果Table 4 LDA and cross-validation results

本次判別分析采用礦物元素結(jié)合穩(wěn)定同位素建立2 個(gè)典型判別式函數(shù)，其中函數(shù)1的方差貢獻(xiàn)率為88.9%，函數(shù)2的方差貢獻(xiàn)率為11.1%，利用這2 個(gè)判別函數(shù)分別作為橫縱坐標(biāo)建立坐標(biāo)系，得到3 個(gè)組的散點(diǎn)圖，由圖2可以直觀地看出3 個(gè)地區(qū)的所有60 個(gè)樣品均被明顯分為3 個(gè)區(qū)域。說(shuō)明本次判別分析建立的判別函數(shù)對(duì)來(lái)自3 個(gè)地區(qū)的樣品具有很好的分類效果。

圖2 典型判別式函數(shù)分類圖Fig. 2 Good classification of three regional samples with typical discriminant functions

2.4 PCA結(jié)果

PCA運(yùn)用數(shù)學(xué)降維的分析方式，能夠?qū)捉M具有相關(guān)性或相似性質(zhì)的變量合并為一組，選取主要成分，減少分析變量，最終實(shí)現(xiàn)簡(jiǎn)化判別模型的目的。本研究通過(guò)設(shè)置特征值大于1的成分作為PC，最終篩選出4 個(gè)PC，累計(jì)方差貢獻(xiàn)率為74.8%，其中PC1主要綜合了K、δ15N的信息，PC2中主要綜合了Mg、Ca的信息，PC3中的主要貢獻(xiàn)元素為Mn、Fe，PC4中的主要貢獻(xiàn)元素為Na、Zn。通過(guò)比較各個(gè)主成分組合對(duì)樣品的分類效果，選取PC1、PC2、PC4分別作為Z軸、Y軸和X軸作散點(diǎn)圖，得到三地區(qū)樣品在三維坐標(biāo)系中的空間分布圖，由圖3可以看出，3 個(gè)地區(qū)的樣品分別分布在空間坐標(biāo)系的3 個(gè)不同的區(qū)域，PCA結(jié)果表明，采用PC1、PC2、PC4可對(duì)三地區(qū)的樣品實(shí)現(xiàn)正確分類，說(shuō)明K、Mg、Ca、Na、Zn和δ15N是區(qū)分3 個(gè)地區(qū)河蟹的特征指標(biāo)。

圖3 三地區(qū)樣品在PC坐標(biāo)系中的分布Fig. 3 Distribution of three regional samples in the PC coordinate system

3 討 論

由于生物體本身不能合成礦物元素，其體內(nèi)的礦物元素含量均來(lái)自日常飲食和其生長(zhǎng)的環(huán)境中[25-26]，因此兩者之間存在顯著相關(guān)性，這就為利用礦物元素區(qū)分不同地區(qū)的水產(chǎn)品提供了理論基礎(chǔ)，郭利攀等[27]利用電感耦合等離子體質(zhì)譜儀和原子吸收光譜測(cè)定東海3 個(gè)地區(qū)經(jīng)濟(jì)魚類的25 種金屬元素，采用偏最小二乘-判別分析對(duì)3 個(gè)地區(qū)的樣品實(shí)現(xiàn)了97.92%的正確分類，并且通過(guò)主成分分析篩選出對(duì)魚類產(chǎn)地區(qū)分影響最為顯著的礦物元素，在自然狀態(tài)下，生物體內(nèi)的穩(wěn)定同位素組成由于氣候、環(huán)境和生物代謝類型等影響因素發(fā)生自然分餾效應(yīng)而存在一定的差異[28-29]，這種差異是生物體一種自然屬性，可為區(qū)分不同來(lái)源的生物提供判別依據(jù)。郭捷敏等[30]比較了野生和人工養(yǎng)殖環(huán)境中的紅樹林青蟹的C、N穩(wěn)定同位素組成，發(fā)現(xiàn)野生環(huán)境下的青蟹與人工養(yǎng)殖環(huán)境下的青蟹具有不同的C、N穩(wěn)定同位素組成，可作為青蟹產(chǎn)地區(qū)分的主要依據(jù)。馬冬紅等[31]利用氫穩(wěn)定同位素對(duì)我國(guó)廣東、海南、廣西和福建4 個(gè)地區(qū)的羅非魚實(shí)現(xiàn)了85%的判別正確率，Sant’Ana等[32]發(fā)現(xiàn)C、N穩(wěn)定同位素指標(biāo)的結(jié)合在無(wú)論雨季還是旱季都能夠?qū)Π臀鞯T魚提供很好的溯源效果，在本研究中，N穩(wěn)定同位素在區(qū)分3 個(gè)地區(qū)的樣品時(shí)表現(xiàn)出一定的優(yōu)勢(shì)，δ15N在3 個(gè)地區(qū)蟹肉中的含量由大到小為陽(yáng)澄湖＞崇明＞興化，通過(guò)對(duì)3 個(gè)地區(qū)河蟹人工餌料和池塘底泥中N穩(wěn)定同位素的測(cè)定，發(fā)現(xiàn)3 個(gè)地區(qū)餌料中δ15N豐度比由大到小為陽(yáng)澄湖（3.78）＞崇明（2.95）＞興化（2.20），底泥中δ15N豐度比由大到小為陽(yáng)澄湖（9.46）＞崇明（3.35）＞興化（1.60），雖然河蟹在長(zhǎng)期的飼養(yǎng)過(guò)程中攝入食物的種類較為復(fù)雜，除了人工餌料還可能攝入養(yǎng)殖水域中其他小型魚類和植物，但本研究證實(shí)了δ15N在3 個(gè)地區(qū)的蟹體、餌料和底泥中呈現(xiàn)顯著相關(guān)性，表明環(huán)境和飲食方面的N攝入對(duì)蟹體中δ15N具有積極的影響，從而為利用δ15N對(duì)不同地區(qū)河蟹的區(qū)分提供了有力的證據(jù)，通過(guò)將穩(wěn)定同位素豐度與礦物元素含量結(jié)合對(duì)3 個(gè)地區(qū)的樣品實(shí)現(xiàn)了98.3%的正確分類。這些結(jié)果表明，不管是采用多種礦物元素還是采用礦物元素與穩(wěn)定同位素結(jié)合的方法都能夠產(chǎn)生很好的判別效果，并且采用礦物元素和穩(wěn)定同位素相結(jié)合的方法能更好地提高判別模型的穩(wěn)定性。

多元統(tǒng)計(jì)分析能夠在多個(gè)對(duì)象和多個(gè)指標(biāo)相互關(guān)聯(lián)的情況下分析它們的統(tǒng)計(jì)規(guī)律，其主要內(nèi)容包括多元線性回歸分析、PCA與因子分析、DA與聚類分析[33]。近年來(lái)隨著科研工作的深入，多元統(tǒng)計(jì)分析在食品的產(chǎn)地溯源領(lǐng)域得到了越來(lái)越廣泛的應(yīng)用[34-35]。本研究采用8 種礦物元素和2 種穩(wěn)定同位素作為變量，對(duì)崇明、興化和陽(yáng)澄湖3 個(gè)地區(qū)的河蟹進(jìn)行方差分析和顯著性比較，發(fā)現(xiàn)各變量在3 個(gè)地區(qū)的樣品都存在不同程度的差異，同時(shí)聚類分析和LDA在對(duì)3 個(gè)地區(qū)的樣品分類方面也取得了良好的效果，聚類分析大體上將興化地區(qū)的樣品分為一類，將陽(yáng)澄湖和崇明的樣品共同歸為一類，這可能是因?yàn)殛?yáng)澄湖和崇明兩個(gè)地區(qū)在地理位置上距離較近，環(huán)境條件相似所致，LDA對(duì)3 個(gè)地區(qū)樣品的交叉驗(yàn)證判別正確率為98.3%，證明線性判別分析對(duì)三地區(qū)河蟹的判別分類具有很好的效果。PCA最終篩選出4 個(gè)主成分，利用不同元素在各個(gè)PC中的載荷及各元素在3 個(gè)地區(qū)樣品中含量的差異性，最終選取了PC1、PC2和PC4成功將3 個(gè)地區(qū)的樣品分成了3 個(gè)不同的區(qū)域，從而提取出K、Mg、Ca、Na、Zn和δ15N 6個(gè)特征變量，為判別模型的優(yōu)化提供了方向。

4 結(jié) 論

本研究利用原子吸收光譜儀和穩(wěn)定同位素比率質(zhì)譜儀分別測(cè)定了崇明、陽(yáng)澄湖和興化3 個(gè)地區(qū)河蟹肌肉組織中的8 種礦物元素含量和C、N穩(wěn)定同位素豐度，結(jié)果顯示，這8 種元素和2 種穩(wěn)定同位素豐度在3 個(gè)地區(qū)間存在不同程度的差異，其中K、Mg和Cu這3 種元素在3 個(gè)地區(qū)間都存在顯著差異（P＜0.05），系統(tǒng)聚類結(jié)果整體上將興化地區(qū)的樣品單獨(dú)歸為一類，將陽(yáng)澄湖和崇明樣品共同歸為一類，采用LDA對(duì)3 個(gè)地區(qū)的樣品建立判別模型，得到崇明、陽(yáng)澄湖、興化三地區(qū)樣品的初始判別正確率分別為100%、95%和100%，整體初始判別正確率為98.3%，交叉驗(yàn)證結(jié)果為98.3%，表明采用這8 種礦物元素建立判別模型可以對(duì)3 個(gè)地區(qū)的樣品實(shí)現(xiàn)正確區(qū)分。利用PCA篩選出4 個(gè)PC，主要貢獻(xiàn)元素為K、Mg、Ca、Na、Zn和δ15N，這4 個(gè)PC的累計(jì)方差貢獻(xiàn)率為74.8%，選用PC1、PC2、PC4建立坐標(biāo)系，繪制散點(diǎn)圖，對(duì)三地區(qū)的樣品實(shí)現(xiàn)了正確區(qū)分。本研究結(jié)果表明，礦物元素結(jié)合穩(wěn)定同位素的分析技術(shù)對(duì)中華絨螯蟹的產(chǎn)地溯源具有一定的優(yōu)越性，可以為全國(guó)范圍內(nèi)的河蟹產(chǎn)地溯源鑒別提供理論基礎(chǔ)和方法參考。