謝建軍，曾廣豐，丁博，王璐，侯穎燁，陳文銳，王志元，李菊

（1.廣州海關(guān)技術(shù)中心，廣州 510623； 2.廣州醫(yī)科大學(xué)，廣州 511436)

大豆及其制品是我國人民膳食結(jié)構(gòu)中重要的組成部分，其富含植物蛋白質(zhì)、微量元素和多種維生素等成分，具有較高的營養(yǎng)價值［1］。大豆榨油后產(chǎn)生的豆粕也是我國養(yǎng)殖業(yè)畜禽飼料的重要來源［2］。然而當(dāng)前我國本土生產(chǎn)的大豆產(chǎn)量已經(jīng)不能滿足國內(nèi)市場的需要，每年都要進(jìn)口大量大豆以填補(bǔ)國內(nèi)市場的缺口。特別是近兩年，中國大豆進(jìn)口量穩(wěn)居世界第一，2019 年我國進(jìn)口大豆8 851 萬噸，2020 年全年累計進(jìn)口大豆10 032.82 萬噸，同比增幅為11.7%［3-5］。我國進(jìn)口大豆的主要國家是巴西、阿根廷、美國和加拿大等國［5］。2017 年，我國進(jìn)口大豆總量的53.3%來自巴西，34.4%來自美國。2018 年4 月中美爆發(fā)了貿(mào)易爭端，我國開始對進(jìn)口來自美國的大豆等農(nóng)產(chǎn)品征收25%的關(guān)稅，懲罰與制裁美國來源的大豆等農(nóng)產(chǎn)品貿(mào)易商［6］。受到中美貿(mào)易摩擦影響，我國進(jìn)口大豆來源發(fā)生劇烈變化，2018 年我國進(jìn)口大豆中巴西大豆占進(jìn)口大豆總量的75.1%，美國大豆占比下降至18.9%［7］。美國產(chǎn)大豆是否通過第三國逃避關(guān)稅制裁進(jìn)入我國，對進(jìn)入我國的大豆的產(chǎn)地溯源問題就變得尤其重要。

氫、氧、碳、氮、硫等穩(wěn)定同位素比值具有強(qiáng)大溯源、示蹤能力，近30 年來，穩(wěn)定同位素分析技術(shù)已在食品、考古、環(huán)境等領(lǐng)域取得了長足發(fā)展和應(yīng)用［8］，特別是在食品的真假、摻偽和產(chǎn)地溯源方面，如葡萄酒、飲料、蜂蜜、咖啡等成功應(yīng)用［9-16］，解決了部分棘手的食品真假及產(chǎn)地溯源問題。在大豆及其制品的摻偽、產(chǎn)地溯源等領(lǐng)域，近年來也有部分學(xué)者進(jìn)行了研究［17-19］。而利用元素分析-同位素比質(zhì)譜儀（EA-IRMS）技術(shù)，測定大豆水溶性蛋白質(zhì)的內(nèi)源性穩(wěn)定同位素δ13C、δ15N、δ18O 和δ2H 值來進(jìn)行進(jìn)口大豆的產(chǎn)地溯源尚未見文獻(xiàn)報道。

筆者利用IRMS 技術(shù)測定大豆水溶性蛋白質(zhì)的內(nèi)源性穩(wěn)定同位素δ13C、δ15N、δ18O 和δ2H 值，并結(jié)合數(shù)學(xué)統(tǒng)計學(xué)方法，研究探討了大豆產(chǎn)地溯源鑒別方法，構(gòu)建了識別準(zhǔn)確率高的大豆產(chǎn)地溯源鑒定模型?？蔀槲覈ｊP(guān)對特定貿(mào)易國實(shí)施特殊關(guān)稅措施提供有力的技術(shù)支撐，有效防止和打擊特定國家的大豆通過第三國轉(zhuǎn)運(yùn)出口至我國的非法行為，從而維護(hù)國家利益。

1 實(shí)驗部分

1.1 主要儀器與試劑

穩(wěn)定同位素比質(zhì)譜儀：Delta V 型，美國賽默飛世爾科技公司。

真空冷凍干燥機(jī)：Alpha 1-2LD PLUS 型，德國Christ 公司。

粉碎機(jī)：JL064 型，上海嘉定糧油有限公司。

超純水儀：Milli-Q Elix-advantage A10 型，美國密理博公司。超聲儀：SW 30H 型，瑞士SONO SWISS 公司。高速冷凍離心機(jī)：SIGMA 3-30KS 型，德國西格瑪離心機(jī)公司。

電子分析天平：Sartorius BT 223S 型，分度值為0.001 g，德國賽多利斯公司。

He 氣（載氣）、CO2標(biāo)準(zhǔn)參考?xì)?、CO 標(biāo)準(zhǔn)參考?xì)狻2標(biāo)準(zhǔn)參考?xì)?、N2標(biāo)準(zhǔn)參考?xì)猓杭兌龋w積分?jǐn)?shù)）均大于99.999 5%，空氣化工產(chǎn)品(中國)投資有限公司。

超濾離心管：3K，美國密理博公司。

EA-IRMS 用碳、氮同位素標(biāo)準(zhǔn)物質(zhì)：IAEA-600咖啡因（δ13C= -27.771‰，δ15N=1.00‰，分別以V-PDB 和空氣中的N 為基準(zhǔn)），天津歐捷科技有限公司。

EA-IRMS 用氫、氧同位素標(biāo)準(zhǔn)物質(zhì)：EMA-P1聚酰胺（δ2H= -25.3‰，δ18O=20.99‰，以V-SMOW為基準(zhǔn)），英國Elemental Microanalysis 公司。

大豆樣品：100 份，分別產(chǎn)自阿根廷（ARG）、美國（USA）、加拿大（CAN）和巴西（BRA）。

實(shí)驗用水為超純水，電阻率為18.2 MΩ·cm。

1.2 樣品處理

將樣品直接粉碎制備30 g 大豆粉，稱取1 g 大豆粉，加入10 mL 超純水，超聲30 min 提取水溶性蛋白質(zhì)，以10 000 r／min 轉(zhuǎn)速離心5 min，取1 mL上清液加入3K 超濾離心管中，離心分別獲得相對分子質(zhì)量大于3 kDa 的蛋白質(zhì)，轉(zhuǎn)移至進(jìn)樣瓶中，冷凍干燥8 h 以上至全干，獲得固態(tài)蛋白質(zhì)粉末［20-21］。

1.3 EA-IRMS 檢測方法

用錫杯和銀杯分別包裹1.2 中的固態(tài)蛋白質(zhì)粉末，校正EA-IRMS，分開檢測樣品中穩(wěn)定碳和氮同位素比，以及穩(wěn)定氫和氧同位素比。錫杯包裹的樣品引入元素分析儀的980 ℃高溫燃燒爐，獲得CO2和N2氣體，進(jìn)入穩(wěn)定同位素比質(zhì)譜儀主機(jī)檢測樣品中δ15N 和δ13C 值；銀杯包裹的樣品引入元素分析儀的1 380 ℃高溫裂解爐中，獲得CO 和H2氣體，導(dǎo)入穩(wěn)定同位素比質(zhì)譜主機(jī)進(jìn)行檢測樣品中δ2H和δ18O 值。

1.4 同位素比數(shù)值校正

樣品測定時，每個樣品的分析起始階段通入高純CO、H2、N2和CO2作為參考?xì)怏w，然后用已知δ13C、δ15N、δ18O 和δ2H 值的標(biāo)準(zhǔn)物質(zhì)標(biāo)定參考?xì)怏w的δ值，再以參考?xì)怏w的標(biāo)定δ值為標(biāo)準(zhǔn)，校正樣品中的δ值。數(shù)據(jù)處理及分析均通過Isodat3.0軟件（美國賽默飛世爾公司)進(jìn)行。δ13C、δ15N、δ18O 和δ2H 的計算以δ13C 為例，大豆水溶性蛋白中δ13C 的計算是基于國際標(biāo)準(zhǔn)物質(zhì)美國白堊紀(jì)皮狄組層位中的擬箭石化石（V-PDB），按式(1)計算［22］：

1.5 數(shù)據(jù)統(tǒng)計分析

將樣品中水溶性蛋白質(zhì)中的δ13C、δ15N、δ18O和δ2H 值歸納整理成excel 表格，導(dǎo)入SIMCA 14.1軟件(瑞典Umetrics 公司)中，進(jìn)行主成分分析（PCA）與正交偏最小二乘法判別分析（OPLS-DA），獲得不同產(chǎn)地大豆同位素分布規(guī)律，從而構(gòu)建大豆產(chǎn)地溯源的同位素鑒別模型。

2 結(jié)果與討論

2.1 EA-IRMS 色譜圖

同位素質(zhì)譜不同于常規(guī)測定樣品中目標(biāo)物相對分子質(zhì)量的有機(jī)質(zhì)譜，它是通過間接檢測CO2、H2、CO、N2小分子化合物來確定碳、氫、氧、氮的同位素豐度比。通過CO2目標(biāo)氣體檢測樣品的δ13C值，以H2目標(biāo)氣體檢測樣品的δ2H，以CO 目標(biāo)氣體檢測樣品的δ18O，以N2目標(biāo)氣體測定樣品的δ15N。樣品同位素比色譜圖如圖1～圖4 所示。

圖1 大豆水溶性蛋白質(zhì)δ 15N 的EA-IRMS 色譜圖

圖2 大豆水溶性蛋白質(zhì)δ 13C 的EA-IRMS 色譜圖

圖3 大豆水溶性蛋白質(zhì)δ 2H 的EA-IRMS 色譜圖

圖4 大豆水溶性蛋白質(zhì)δ 18O 的EA-IRMS 色譜圖

2.2 大豆水溶性蛋白質(zhì)的穩(wěn)定同位素比值

采用EA-IRMS 同位素比質(zhì)譜儀分析測定大豆水溶性蛋白質(zhì)的δ13C、δ15N、δ18O 和δ2H 值見表1。大豆屬于C3植物，根據(jù)文獻(xiàn)［23］，C3植物的δ13C 值為-20‰～-35‰(平均值為-26‰)。本研究大豆水溶性蛋白質(zhì)測定的δ13C 值為-24.675‰～-29.490‰，利用t分布統(tǒng)計學(xué)分析60 種大豆樣品水溶性蛋白質(zhì)δ13C 平均值為-26.931‰，標(biāo)準(zhǔn)偏差為0.995%。由此可知，本方法檢測大豆水溶性蛋白質(zhì)δ13C 值與文獻(xiàn)報道數(shù)據(jù)基本一致。同時，由表1 可知，測定的大豆水溶性蛋白質(zhì)δ15N值為-1.595‰～2.537‰，平均值為0.320‰，標(biāo)準(zhǔn)偏差為1.046%；δ18O 值為15.749‰～25.127‰，平均值為19.503‰，標(biāo)準(zhǔn)偏差為2.250%；δ2H 值為-33.504‰～-118.697‰，平均值為-83.534‰，標(biāo)準(zhǔn)偏差為21.365%。

表1 大豆水溶性蛋白質(zhì)的δ 13C、δ 15N、δ 18O 和δ 2H 值

續(xù)表1

2.3 統(tǒng)計分析方法比較

利用SIMCA14.1 軟件對測定的大豆水溶性蛋白質(zhì)δ13C、δ15N、δ18O 和δ2H 值進(jìn)行多元變量統(tǒng)計分析，分別建立PCA 和OPLS-DA 兩種統(tǒng)計分析方法的大豆產(chǎn)地溯源鑒別模型，如圖5 和圖6 所示（橢圓形區(qū)域為Hotelling T2 檢驗，95%置信區(qū)間）。

圖5 大豆水溶性蛋白質(zhì)δ 13C、δ 15N、δ 18O 和δ 2H 的PCA 鑒別模型

圖6 大豆水溶性蛋白質(zhì)δ 13C、δ 15N、δ 18O 和δ 2H 的OPLS-DA 鑒別模型

PCA 溯源鑒定模型結(jié)果顯示，利用大豆水溶性蛋白質(zhì)中的δ13C、δ15N、δ18O 和δ2H 值可區(qū)分四個國家的大豆樣品。構(gòu)建OPLS-DA 大豆產(chǎn)地溯源鑒定模型獲得大豆中水溶性蛋白質(zhì)的4 個變量δ13C、δ15N、δ18O 和δ2H 對產(chǎn)地區(qū)分模型的貢獻(xiàn)度用VIP值評定。結(jié)果表明，δ15N 在產(chǎn)地溯源模型中對產(chǎn)地溯源的貢獻(xiàn)最大，其次是δ2H 和δ18O，δ13C 的貢獻(xiàn)最低，其VIP 值依次為1.069，1.045，1.035，0.834。對比分析大豆產(chǎn)地溯源的兩種統(tǒng)計分析方法的模型鑒定準(zhǔn)確率，結(jié)果表明，OPLS-DA 模型分析結(jié)果顯著，建立的識別模型對美國、巴西和阿根廷產(chǎn)大豆的識別準(zhǔn)確率均為100%，對加拿大產(chǎn)大豆的識別準(zhǔn)確率為90%；PCA 識別模型對美國和阿根廷產(chǎn)大豆的識別準(zhǔn)確率均為80%，對巴西和加拿大產(chǎn)大豆的識別準(zhǔn)確率均為90%，因此最終選擇OPLS-DA統(tǒng)計方法建立的模型進(jìn)行大豆產(chǎn)地溯源識別分析。

2.4 大豆產(chǎn)地鑒定模型建立

分別對美國與其它國家大豆水溶性蛋白質(zhì)δ13C、δ15N、δ18O 和δ2H 值采用OPLS-DA 建立兩兩國家同位素產(chǎn)地溯源鑒定模型，結(jié)果如圖7～圖9 所示（橢圓形區(qū)域為Hotelling T2檢驗，95%置信區(qū)間）。

圖7 美國-加拿大大豆產(chǎn)地溯源OPLS-DA 鑒別模型

圖8 美國-巴西大豆產(chǎn)地溯源OPLS-DA 鑒別模型

圖9 美國-阿根廷大豆產(chǎn)地溯源OPLS-DA 鑒別模型

圖7～圖9 結(jié)果表明，所建立的同位素產(chǎn)地溯源鑒定模型識別準(zhǔn)確率達(dá)到100%，應(yīng)用該方法可準(zhǔn)確鑒別美國與其他國家大豆產(chǎn)地。

2.5 大豆產(chǎn)地鑒定模型驗證

為了進(jìn)一步驗證上述模型的可靠性，選用四個國家各10 份非建模大豆樣品，提取其中的水溶性蛋白質(zhì)，測定δ13C、δ15N、δ18O 和δ2H 值，分別帶入已建立的OPLS-DA 兩兩國家識別模型進(jìn)行驗證，結(jié)果見表2。盲樣驗證結(jié)果表明，除了加拿大產(chǎn)大豆中有1 份被誤判為美國產(chǎn)外，其他3 個國家盲樣均可準(zhǔn)確識別，識別準(zhǔn)確率為100%。

表2 OPLS-DA 識別模型的盲樣驗證

3 結(jié)語

建立了基于EA-IRMS 測定大豆水溶性蛋白質(zhì)δ13C、δ15N、δ18O 和δ2H 的OPLS-DA 大豆產(chǎn)地溯源模型。結(jié)果表明，美國與加拿大、阿根廷、巴西產(chǎn)的大豆在建立的兩兩識別模型中可準(zhǔn)確鑒別，識別準(zhǔn)確率均可達(dá)100%。該方法和所構(gòu)建的模型準(zhǔn)確可靠，可有效防止和打擊特定國家的大豆通過第三國轉(zhuǎn)運(yùn)出口至我國的非法行為，維護(hù)國家利益。