具体实施方式

为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对照附图说明本发明的具体实施方式。显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图，并获得其他的实施方式。

下面以具体实施例详细介绍本发明的技术方案。

大豆油样品数据采集过程中的实验器材包括：实验器材包括：选择美国ABsciex公司的Triple TOF5600+高分辨质谱仪；日本岛津公司的HPLC 20AD高效液相色谱仪；美国Waters公司的Xbridge BEH C18色谱柱(100mm×2.1mm,3μm)；德国IKA公司的 VORTEX 4旋涡混匀器；以及榨油机；

试剂选择：甲醇、乙腈、丙酮、乙酸乙酯(色谱纯，赛默飞世尔公司，美国)；

大豆粉样品数据采集过程中的实验器材包括：ICP 5300电感耦合等离子体发射光谱仪 (美国PE公司)；Elane-drc-e电感耦合等离子体质谱仪(美国PE公司)；ETHosA型微波消解仪(意大利Milestone公司)；Ultra WAVE超级微波消解平台(意大利Milestone公司)；DS-360石墨消解仪(广州分析测试中心)；Milli-Q超纯水系统(电阻率18.2MΩ·cm)；所有玻璃器皿均用20％HNO₃浸泡24h，用水反复冲洗，最后用去离子水冲洗干净，待用。

IV-ICP-MS-71A标准溶液(10μg/mL，inorganic Ventures)；CCS-5标准溶液：(100μg/mL，inorganic Ventures)；GNM-M292793-2013：ICP用多元素标准溶液：国家有色金属及分析测试中心；硝酸：优级纯，苏州晶瑞有限公司生产。Rh、Re、In混合内标：(100μg/mL，inorganic Ventures)；水：经超纯水机制备(18.3MΩ.Cm)；其他所用试剂均为优级纯以上；氮气：99.99％以上。

大豆水溶性蛋白质的数据采集中使用的仪器：Flash EA2000元素分析仪，DELTAVAdvantage同位素质谱仪和ConFlo IV接口(美国Thermo Scientific公司)；Milli-Q超纯水仪(美国Millipore公司)；SW 30H型超声仪(瑞士SONO SWISS公司)；SIGMA3-30KS 高速冷冻离心机(德国Sigma离心机公司)；Sartorius BT 223S千分之一天平(德国赛多利斯公司)；榨油机；真空冷冻干燥机(德国Christ公司，Alpha 1-2LD PLUS)；粉碎机(上海嘉定粮油有限公司，JL064)；银杯(德国IVAAnalysentechnik公司)；锡杯(美国赛默飞公司)；作为载气的He气(纯度大于99％)、CO₂标准参考气(纯度大于99.9995％)、 CO标准参考气(纯度大于99.9995％)、H₂标准参考气(纯度大于99.9995％)；N₂标准参考气(纯度大于99.9995％)；超滤离心管(美国密理博公司，Centrifugal Filter Units，3K/10K)。

EA-IRMS用碳、氮同位素标准物质：IAEA-600咖啡因：δ¹³C＝-27.771‰，δ¹⁵N＝+1.00‰， V-PDB为基准，购于国际原子能机构；EA-IRMS用氢氧同位素标准物质：EMA-P1聚酰胺：δ²H＝-25.3‰，δ¹⁸O＝20.99‰，以V-SMOW为基准，购于英国Elemental Microanalysis公司；实验用水为超纯水，阻抗值为18.2MΩ·cm。

大豆及大豆油样品来源：大豆标准样品来源于全国有关直属海关及国外购买共收集大豆样品807个，其中各个产地的大豆样品数量如下：其中美国样品152个，巴西样品423个，加拿大样品96个，阿根廷样品68个，乌拉圭样品14个，俄罗斯样品25个，中国样品29个；大豆油样品334个：其中美国样品36个，巴西样品226个，乌克兰样品7个，阿根廷样品46个，墨西哥样品1个，俄罗斯样品16个。

基于IRMS、LC-Q-TOF-MS和多元素分析技术相结合的大豆原产地溯源识别方法，包括如下步骤：

S1.标准样品制备：选取具有明确地区的不同产地的大豆样品，每个产地的所述大豆样品分别经压榨后得到大豆油样品、经研磨后得到标准大豆粉样品、以及提取所述标准大豆粉样品的水溶性蛋白质作为靶向目标物，将大豆油样品用稀释溶剂进行逐级稀释100～200 倍后得到标准大豆油样品，所述大豆样品来源于至少两个不同产地；

S2.标准样品的质谱数据、矿质元素含量数据和稳定同位素比值数据采集：用液相色谱- 四级杆飞行时间质谱仪分别采集不同地区所述标准大豆油样品的质谱数据信息，获得所述标准大豆油样品的甘油三酯化合物的IDA-MS高分辨质谱数据；用矿质元素含量分析仪采集不同产地的所述标准大豆粉样品的质谱数据信息，获得所述标准大豆粉样品的矿质元素含量数据；以及用同位素质谱仪采集不同产地的所述标准大豆粉样品的水溶性蛋白质的质谱数据信息，获得所述水溶性蛋白质的稳定同位素比值数据；

S3.标准大豆油样品靶向化合物的确定：将步骤S2获得的IDA-MS高分辨质谱数据中的甘油三酯化合物质谱数据与脂类化合物数据库中的标准甘油三酯化合物数据进行匹配，确定所述标准大豆油样品的甘油三酯化合物靶向化合物；

S4.大豆产地溯源鉴定模型建立：确定标准大豆油样品中甘油三酯靶向化合物后，将所述IDA-MS高分辨质谱数据通过分析软件进行分析，以获取所述标准大豆油样品的甘油三酯化合物标志物观测峰值数据，将同一产地大豆的所述标准大豆油样品中的甘油三酯化合物标志物观测峰值数据、和标准大豆粉样品的矿质元素含量数据、以及标准大豆粉样品的水溶性蛋白质的稳定同位素比值数据同时导入建模软件进行主成分分析法、偏最小二乘法- 判别分析法和正交偏最小二乘法回归分析法中的一种或多种方式处理，然后依次将不同产地大豆的标准大豆油样品的甘油三酯化合物标志物观测峰值数据、和标准大豆粉样品的矿质元素含量数据、以及标准大豆粉样品的水溶性蛋白质的稳定同位素比值数据均导入建模软件中，以获得不同产地标准大豆油样品中大豆油脂特征分布规律、标准大豆粉样品矿质元素以及标准大豆粉样品的水溶性蛋白质的稳定同位素分布规律，构建基于LC-Q-TOF-MS、矿质元素含量和稳定同位素比值分析相结合的大豆产地溯源鉴定模型；所述标志物观测峰值数据即为MarkerView peaks数据，通过MasterView分析软件分析获得标准大豆油样品的甘油三酯的MarkerView peaks数据，

S5.结果预测：将待检测大豆样品分别经过压榨得到大豆油样品和经研磨后得到待检测大豆粉样品、提取待检测大豆粉样品得到待检测大豆水溶性蛋白质，将大豆油样品用稀释溶剂进行逐级稀释100～200倍后得到待检测大豆油样品，通过基液相色谱-四级杆飞行时间质谱仪采集待检测大豆油样品的IDA-MS高分辨质谱数据，用矿质元素含量分析仪采集待所述待检测大豆粉样品的矿质元素含量数据，以及用同位素质谱仪采集待检测大豆水溶性蛋白质的稳定同位素比值数据；将待检测大豆样品的IDA-MS高分辨质谱数据、矿质元素含量数据和稳定同位素比值数据同时导入步骤S4的所述大豆产地溯源鉴定模型中，进行待检测大豆的产地溯源预测。

多国OPLS-DA大豆产地溯源鉴定模型实施例1

选取具有明确地区的不同产地的大豆样品，每个产地的所述大豆样品分别经压榨后得到大豆油样品、经研磨后得到标准大豆粉样品、以及提取所述标准大豆粉样品的水溶性蛋白质作为靶向目标物，将大豆油样品用稀释溶剂进行逐级稀释200倍后得到标准大豆油样品，选取甲醇-乙酸乙酯混合液作为稀释溶剂，所述稀释溶剂中的甲醇和乙酸乙酯比例为1： 1。

采集标准大豆油样品质谱数据：

将稀释后得到的所述标准大豆油样品置于所述液相色谱-四级杆飞行时间质谱仪的进样器中，通过所述液相色谱-四级杆飞行时间质谱仪中的液相色谱仪对所述标准大豆油样品进行分离分析，然后再通过所述液相色谱-四级杆飞行时间质谱仪中的质谱仪进行所述标准大豆油样品的质谱数据采集，通过质谱仪的一级TOF-MS扫描和二级IDA-MS扫描分别得到所述标准大豆油样品的一级质谱信息和二级质谱信息，所述二级质谱信息为IDA-MS高分辨质谱数据，通过IDA-MS高分辨质谱数据确定所述标准大豆油样品的甘油三酯化合物靶向化合物，确定靶向化合物的方法包括：根据所述IDA-MS高分辨质谱数据中的目标物分子量范围，将大豆油样品中的IDA-MS高分辨质谱数据分为3个区域：分子量在800-1000 的第一区域，分子量在550-800的第二区域，以及分子量在550以下的第三区域，所述第一区域与第二区域是大豆油油脂代谢特征区域，将脂类化合物数据库中的分子量范围在 700-950之间的动植物来源甘油三酯化合物与第一区域和第二区域中的大豆油甘油三酯化合物进行匹配筛选，确定分子量范围在766-920Da之间的114个大豆油中的甘油三酯化合物作为脂类组学分析靶向化合物，应用于大豆油中甘油三酯脂质组学产地溯源的鉴定模型研究之中。

脂类化合物数据库为美国公布的脂类化合物数据库LIPID MAPS LipidomicsGateway，标准甘油三酯化合物数据为脂类化合物数据库LIPID MAPS Lipidomics Gateway中 Triradylglycerols的6899个甘油三酯化合物，通过PeakView软件定性分析确定大豆油样品常见的116个甘油三酯化合物为所述标准大豆油样品的甘油三酯化合物靶向化合物。

所述液相色谱-四级杆飞行时间质谱仪中液相色谱仪的液相色谱条件为：流速为0.5 μL/min，柱温为40℃，Xbridge BEH C18色谱柱梯度洗脱，2μL进样量；流动相中A相为异丙醇，流动相中的B相为乙腈，其中流动相在不同时间段内的B相含量如下：0min,70％ B；0-5min,70-65％B；5-8min,65％B；10-10.5min,65-70％B；10.5-15min,70％B。

所述液相色谱-四级杆飞行时间质谱仪中质谱仪的四级杆飞行时间质谱条件为：质谱仪采用正离子模式采集数据，离子源为：ESI和APCI复合源；正离子扫描方式为：APCI源连接自动校正系统，一级TOF-MS扫描准确质量范围：100～2000Da，数据采集时间100ms，DP:100V，CE:10V；二级IDA-MS扫描准确质量范围：50～2000Da，DP:100V，CE:35±15V；所述质谱仪采用高灵敏模式，数据采集时间50ms，信号阈值100cps，每次循环采集6次数据，且采用动态背景减法扣除。

所述液相色谱-四级杆飞行时间质谱仪为岛津LC20AD液相色谱仪，所述质谱仪为Triple TOF 5600+质谱仪，自动校正系统为CDS系统，四级杆飞行时间质谱条件还包括：每10个样品自动校正1次，APCI正离子校正液流速0.3mL/min，气帘气压力为：40psi，离子源雾化气压力为：50psi，离子源辅助加热气压力为：50psi，离子源温度为：500℃，离子源电压5500V，质谱仪采集的所有IDA-MS高分辨质谱数据在ABSciex公司的Analyst TF 1.6软件采集，IDA-MS高分辨质谱数据在PeakView，MasterView软件上定性定量处理分析后，导入SIMCA14.0软件(瑞士Umetrics公司)中，将所述标准大豆油样品中的甘油三酯化合物标志物观测峰值数据进行正交偏最小二乘法回归分析法处理，获得不同产地大豆的甘油三酯及代谢物的分布规律。

采集标准大豆粉样品矿质元素含量：

称取0.5g标准大豆粉样品于微波消解罐的内罐中，加入4ml硝酸，盖上内罐盖，浸泡 6～12个小时；浸泡后，旋紧微波消解罐的外罐置于微波消解仪中进行消解，然后设定升温程序对消解仪中的标准大豆粉样品进行消解，消解完毕待冷却后取出内罐，于电热板上赶酸，赶酸结束后将内罐中消化液转移至50mL聚丙烯离心管中，并用少量去离子水多次洗涤微波消解罐的内罐，合并洗涤液并用去离子水将聚丙烯离心管中的消化液和去离子水的混合液定容至25mL，混匀备用，同时做空白对比试验，然后矿质元素含量分析仪采用 ICP-MS和/或ICP-AES测定分析标准大豆粉样品中的矿质元素含量。

所述空白对比试验包括两组，其中一组空白对比试验在聚丙烯离心管中放置25mL去离子水；另一组空白对比试验在聚丙烯离心管中放置25mL硝酸。

所述微波消解罐在所述微波消解仪中升温程序包括如下步骤如表1所示：

a1.在8min内将微波消解仪的温度升至120℃，然后保持120℃的温度5min；

a2.然后在10min内将微波消解仪的温度从120℃升至180℃，然后保持180℃的温度 10min。

表1微波消解罐在所述微波消解仪中升温程序

电感耦合等离子体发射光谱法ICP-AES的工作条件为：检测波长为396.152nm；功率为 1500W；等离子体气流量：10.0L/min；辅助气流量为0.55L/min；雾化气流量：0.55L/min；蠕动泵转速为50r/min，观测方式：轴向；载气：99.996％高纯氩气；

电感耦合等离子体质谱法ICP-MS的工作条件为：采用模式:标准(STD)；气体:氩气(≥99.998％,高纯)；扫描/读数:20；读数/重复次数:1；重复次数:3；雾化气流速:0.84L/min；辅助气流量:1.2L/min；等离子气流量:1500L/min；ICP射频功率:1500W；蠕动泵:20rpm。采样锥及截取锥：铂锥。

通过ICP-AES和ICP-MS标准大豆粉样品中矿质元素的含量，其中ICP-AES检测标准大豆粉样品中Si、Mg、Al、P、K、Na、Ca、Mn、Fe、Cu、Sr和Zn等13种常量元素，ICP-MS 检测分析其他39种微量元素。ICP-AES与ICP-MS检测分析获得52种矿质元素含量数据，通过描述统计分析方法，归纳整理excel列表，导入SIMCA 14.1软件(瑞士Umetrics公司) 中进行主成分分析、偏最小二乘法判别分析、正交偏最小二乘法法判别分析构建初步鉴定模型。本实施例还选取Ca、Mg、P、Si、Mn和Zn等6种常量元素，探讨大豆样品中两两元素之间的含量相关性，大豆样品中常量元素Mg和Ca相关系数是0.81，P和Si元素含量相关系数是0.72，Zn和Mn元素含量相关系数0.36，由此可知大豆样品中含量比较高的常量元素，具有内在的一致性，尤其是大豆样品中常量元素Mg和Ca，由此可以推测不同国家地区的大豆中矿质元素含量具有生物内在的一致性，因此通过矿质元素可以鉴别区分不同产地的大豆。

在矿质元素含量的测定过程中，为了数据的准确性，要对大豆样品进行优化，初步构建的基于矿质元素含量产地溯源鉴定模型中一些建模样品数据存在显著离群现象，同时，产地溯源模型参数指标Hotelling’s T2值，也显示一些建模数据是异常值数据，根据数学建模技术路线原理，需要剔除异常值数据，为此本实施例构建大豆多元素分析产地溯源鉴定模型删除异常值数据，进行产地溯源鉴定模型的优化。

标准大豆粉样品的水溶性蛋白质提取过程如下：将大豆样品研磨粉碎制备30g标准大豆粉样品，称取1g标准大豆粉样品，加入10mL超纯水，经过超声波超声30min，10000 转的转动条件下离心5min，提取1mL上清蛋白质提取液加入3K超滤离心管，离心分别获得分子量3k Da以上的蛋白质，将离心后的蛋白质转移至进样瓶中，冷冻干燥8h以上至全部干燥，获得水溶性蛋白质的固态蛋白质粉末。

所述同位素质谱仪为在线气体分析-稳定同位素质谱仪GasBench-IRMS、元素分析-同位素质谱仪EA-IRMS、气相色谱-同位素质谱仪GC-IRMS、液相色谱-同位素质谱仪LC-IRMS 中的一种。本实施例选取元素分析-同位素质谱仪EA-IRMS来分析大豆水溶性蛋白质的稳定同位素比值。大豆水溶性蛋白质的稳定同位素比值分析通过元素分析仪串联同位素质谱仪实现。EA-IRMS可以高精度测定有机或无机样品中氢、碳、氧、氮和硫等的稳定同位素的平均组成。EA-IRMS方法是将样品包裹于锡纸中，称量后直接放入EA自动进样器，此后样品在燃烧管中经过高温燃烧转变为NO_x，CO₂，SO₂或H₂O等气体。根据检测需求，混合气体经气相色谱柱等分离后除去干扰气体，目标气体最终进入同位素质谱仪进行检测。在碳同位素质谱测定中，样品燃烧生成的混合气体随着氦气流进入还原管，在还原管中，NO_x转变为N₂，过量的O₂被除去，随后混合气体经干燥管除去H₂O，再经气相色谱柱实现CO₂和N₂的分离获得单一的CO₂气体。

标准大豆粉样品的水溶性蛋白质的稳定同位素比值数据采集方法为：

a1.选取元素分析-同位素质谱仪EA-IRMS，在锡杯和银杯中均放置称量后的水溶性蛋白质的固态蛋白质粉末，分别检测水溶性蛋白质中稳定碳和氮的同位素比值，以及稳定氢和氧的同位素比值；

a2.调谐校正元素分析-同位素质谱仪EA-IRMS；

a3.将锡杯包裹的水溶性蛋白质放入EA自动进样器中，然后引入980℃高温燃烧炉，获得CO₂和N₂气体，将CO₂和N₂混合氦气导入同位素质谱仪主机中进行水溶性蛋白质中的δ¹⁵N和δ¹³C检测；将银杯包裹的水溶性蛋白质引入1380℃高温裂解炉中，获得CO和H₂气体，将CO和H₂混合氦气导入同位素质谱仪主机中进行水溶性蛋白质的δ²H和δ¹⁸O检测。

所述步骤a2中，调谐校正元素分析-同位素质谱仪EA-IRMS的方法如下：在所述元素分析-同位素质谱仪EA-IRMS中通入高纯度CO、H2、N2和CO2作为参考气体，并通入高纯度He气作为引流气体，然后用已知δ²H、δ¹⁸O、δ¹⁵N和δ¹³C值的标准物质标定所述参考气体的稳定同位素δ值，再以参考气体标定后的稳定同位素δ值作为检测标准，用于校正步骤a3中检测得到的标准大豆粉样品的水溶性蛋白质中的稳定同位素δ值。校正过程中的数据的处理以及分析均通过美国Thermo Scientific公司的Isodat3.0软件实现。已知δ²H、δ¹⁸O、δ¹⁵N和δ¹³C值的标准物质为：碳、氮同位素标准物质IAEA-600咖啡因：δ¹³C＝-27.771‰，δ¹⁵N＝+1.00‰，V-PDB为基准，购于国际原子能机构；氢氧同位素标准物质EMA-P1聚酰胺：δ²H＝-25.3‰，δ¹⁸O＝20.99‰，以V-SMOW为基准。

元素分析-同位素质谱仪EA-IRMS同位素质谱不同于常规测定样品中目标物分子量的有机质谱，而通过CO₂、H₂、CO、N₂小分子化合物的间接检测，确定碳、氢、氧、氮的同位素丰度比。通过CO₂目标气体检测标准大豆粉样品的水溶性蛋白质的δ¹³C值，H₂目标气体检测标准大豆粉样品的水溶性蛋白质的δ²H，CO目标气体检测标准大豆粉样品的水溶性蛋白质的δ¹⁸O，N₂目标气体测定标准大豆粉样品的水溶性蛋白质的δ¹⁵N。

将元素分析-同位素质谱仪EA-IRMS采集标准大豆粉样品的水溶性蛋白质的δ²H、δ¹⁸O、δ¹⁵N和δ¹³C值，通过描述统计分析方法，归纳整理excel列表，导入SIMCA 14.1软件中，进行正交偏最小二乘法回归分析法的方式处理，获得不同产地大豆同位素分布规律，从而构建大豆产地溯源的同位素鉴别模型。如图7所示，选择标准大豆粉样品的水溶性蛋白质的δ²H、δ¹⁸O、δ¹⁵N和δ¹³C构建的OPLS-DA多国鉴定模型，结果显示巴西、美国、阿根廷、乌拉圭及加拿大等国家大豆样品产地区分度相对较好。通过对大豆中水溶性蛋白质的4个变量δ²H、δ¹⁸O、δ¹⁵N和δ¹³C对产地区分的贡献度进行分析，δ¹⁵N对产地溯源的贡献最大，其次是δ²H和δ¹⁸O，δ¹³C的贡献最低。

本实施例选取具有明确地区的巴西大豆样品、美国大豆样品、中国大豆样品、阿根廷大豆样品、加拿大大豆样品、乌拉圭大豆样品、俄罗斯大豆样品来制作标准样品，并将该标准样品中不同地区的大豆分别通过液相色谱-四级杆飞行时间质谱仪采集IDA-MS高分辨质谱数据，并经过PeakView，MasterView软件上定性定量处理分析后获取标准样品中不同产地大豆样品的甘油三酯化合物MarkerView peaks数据；以及用ICP-MS和ICP-AES采集优化后的标准大豆粉样品的矿质元素含量，以及用元素分析-同位素质谱仪EA-IRMS采集标准大豆粉样品的水溶性蛋白质的δ²H、δ¹⁸O、δ¹⁵N和δ¹³C值，将处理后的甘油三酯化合物标志物观测峰值数据、矿质元素含量数据和水溶性蛋白质的稳定同位素比值数据归纳整理 excel列表，导入SIMCA 14.1软件进行正交偏最小二乘法回归分析法处理，从而构建基于 IRMS、LC-Q-TOF-MS和多元素分析技术相结合的大豆产地溯源鉴定模型。

如图1可以看出，单独建立基于LC-Q-TOF-MS多国OPLS-DA大豆产地溯源鉴定模型，可显著区分巴西、俄罗斯、阿根廷及美国产地的大豆样品，单由于美国、加拿大及阿根廷产地样品交叉分布在一起，影响模型的预测鉴定准确度，因此需要对大豆和大豆油产地溯源鉴定模型进行进一步优化，为了避免部分国家或地区的样品在鉴定模型中交叉分布在一起，影响模型的预测鉴定准确度，因此需要对大豆产地溯源鉴定模型进行进一步优化，因此需要对大豆和大豆油产地溯源鉴定模型进行进一步优化，所述优化步骤包括：通过大豆产地溯源鉴定模型的VIP值确定所述靶向化合物中贡献值度大的部分甘油三酯化合物，并根据hotelling’s和DModx指标，删除超出99％置信区间的异常值样品，以及删除在鉴定模型中出现数目比较少的所在地区的全部大豆油样品，以优化大豆和大豆油产地溯源鉴定模型。如图2所示，通过鉴定模型的VIP值确定哪些变量贡献度大，最终确定分子量为873.6967、 875.7123、851.7123、877.7280、853.7280、865.7280、913.7280、829.7280、915.7436、 767.6184、835.6810、879.7436、917.7593、859.7749、855.7436、895.779、825.6967、 921.7906、887.8062及869.7593等甘油三酯化合物对多国OPLS-DA大豆和大豆油产地溯源鉴定模型贡献度大，研究根据hotelling’s和DModx指标，并删除数目比较少的中国和乌拉圭样品，优化OPLS-DA大豆和大豆油产地溯源鉴定模型。如图3所示，通过优化后的基于脂质组学的多国OPLS-DA大豆产地溯源鉴定模型中能够对不同国家的样品较为显著的区分，尤其是显著区分巴西、俄罗斯、美国等产地样品；

如图5所示，大豆样品中Sr、As、Ni、Cd、Mo、Se、Cs、Cr、Mn、Tl、Zr、Co、Be、 U、V、Al、Eu、Na、Ce、Cu、Nd、P、Gd、Sn、Sc、I和Sm元素是贡献度较大的元素，选取这27种元素作为优化后的特征贡献度元素，如4和6所示，经优化后的基于矿质元素含量的多国OPLS-DA大豆产地溯源鉴定模型中能够对不同国家的样品较为显著的区分，尤其是显著区分巴西和非巴西大豆样品；以及巴西、俄罗斯、美国、阿根廷和加拿大产地的样品。

因此，最终的基于LC-Q-TOF-MS、矿质元素含量和稳定同位素比值分析相结合的大豆产地溯源鉴定模型是由优选后的甘油三酯靶向化合物、特征贡献度矿质元素数据、和稳定同位素δ¹⁵N、δ²H、δ¹⁸O和δ¹³C构建而成的。其能够明显区分巴西、俄罗斯、美国、阿根廷和加拿大产地的大豆样品。

多国OPLS-DA大豆产地溯源鉴定模型建立好之后，需要对该模型进行盲样验证，所述盲样验证步骤包括：选取同一地区的多个大豆验证样品，将具有明确地区的大豆验证样品分别经过压榨得到待验证大豆油样品和经研磨、提取后得到待验证大豆粉样品的水溶性蛋白，通过液相色谱-四级杆飞行时间质谱仪采集待验证大豆油样品的IDA-MS高分辨质谱数据，用矿质元素含量分析仪采集待所述待验证大豆粉样品的矿质元素含量数据，以及用同位素质谱仪采集待所述待验证大豆粉样品的水溶性蛋白稳定同位素比值数据；将大豆验证样品的IDA-MS高分辨质谱数据、矿质元素含量数据和水溶性蛋白稳定同位素比值数据同时导入优化后的的所述大豆产地溯源鉴定模型中，验证所述大豆验证样品的产地溯源准确度。

将待检测大豆样品分别经过压榨得到大豆油样品和经研磨后得到待检测大豆粉样品、提取待检测大豆粉样品得到待检测大豆水溶性蛋白质，将大豆油样品用稀释溶剂进行逐级稀释100～200倍后得到待检测大豆油样品，通过基液相色谱-四级杆飞行时间质谱仪采集待检测大豆油样品的IDA-MS高分辨质谱数据，用矿质元素含量分析仪采集待所述待检测大豆粉样品的矿质元素含量数据，以及用同位素质谱仪采集待检测大豆水溶性蛋白质的稳定同位素比值数据；将待检测大豆样品的IDA-MS高分辨质谱数据、矿质元素含量数据和稳定同位素比值数据同时导入优化后的所述大豆产地溯源鉴定模型中，进行待检测大豆的产地溯源预测。

多国大豆产地溯源鉴定模型实施例2

该实施例仅描述与上述实施例的不同之处，在本实施例中，标准大豆粉样品的矿质元素含量数据采集方法为：步骤S2中标准大豆粉样品的矿质元素含量数据采集方法为：称取 0.5g标准大豆粉样品置于聚乙烯消化管中，称取0.5g的标准大豆粉样品的精确为0.001g，加入4ml硝酸，盖上密封盖，浸泡6～12个小时；浸泡后，将密封后的聚乙烯消化管置于消解仪中，并在消解仪的内罐体中预加惰性气体，惰性气体压力4000KPa，然后设定升温程序对消解仪中的标准大豆粉样品进行消解，消解完毕待冷却后取出内罐，于电热板上赶酸，所述电热板设定的温度为150℃。赶酸结束后将内罐中消化液转移至50mL聚丙烯离心管中，并用少量去离子水多次洗涤微波消解罐的内罐，合并洗涤液并用去离子水将聚丙烯离心管中的消化液和去离子水的混合液定容至25mL，混匀备用，同时做空白对比试验，矿质元素含量分析仪采用ICP-MS和/或ICP-AES测定分析标准大豆粉样品中的矿质元素含量。

所述惰性气体为氮气，消解仪的外腔温度≤40℃。

所述聚乙烯消化管在所述消解仪中升温程序包括如下步骤如表2所示：

a1.在15000KPa压力和1500W功率的条件下，在8min内将微波消解仪的温度升至120℃；

a2.然后在15000KPa压力和1500W功率的条件下，在8min内将微波消解仪的温度从120℃升至240℃；

a3.进而在15000KPa压力和1500W功率的条件下，保持240℃的温度持续20min。

表2超级微波升温程序

两国大豆产地溯源鉴定模型实施例1

该实施例仅描述与上述实施例的不同之处，在本实施例中，所述步骤S1中标准大豆油样品制备的过程中，设所述大豆样品来源于n个不同的地区，将所述大豆样品分成n×(n-1)/2 组，每组大豆样品中均由两个不同产地的大豆样品组成，将每组大豆样品均按照步骤S2、 S3、S4建立两国大豆产地溯源鉴定模型，用于鉴定两国大豆产地溯源鉴定模型中所对应的大豆国家或地区。

该实施例中，选取具有明确地区的巴西大豆样品美国大豆样品来制作标准大豆油样品、标准大豆粉样品、及标准大豆粉样品的水溶性蛋白质，并将这两个不同地区的大豆油样品分别通过液相色谱-四级杆飞行时间质谱仪采集IDA-MS高分辨质谱数据，IDA-MS高分辨质谱数据在PeakView，MasterView软件上定性定量处理分析后，获得所述标准大豆油样品中的甘油三酯化合物标志物观测峰值数据，并将这两个不同地区的大豆粉样品分别通过ICP-MS和/或ICP-AES采集矿质元素含量数据，以及将这两个不同地区的大豆粉样品的水溶性蛋白质分别通过EA-IRMS采集稳定同位素比值数据，将甘油三酯化合物标志物观测峰值数据、矿质元素含量数据和稳定同位素比值数据导入SIMCA 14.1软件中进行正交偏最小二乘法回归分析法处理，从而建立基于IRMS、LC-Q-TOF-MS和多元素分析技术相结合的大豆产地溯源模型。从图8中可以看出巴西-美国两国OPLS-DA大豆产地溯源鉴定模型中可显著区分美国和巴西大豆样品，为了进一步验证两国压榨大豆油产地溯源鉴定模型的判定准确度，选择美国大豆样品和巴西大豆样品分别进行压榨和磨粉提取水溶性蛋白质后进行模型盲样验证，验证结果表明，巴西来源的大豆油样品判定准确度为95％，美国来源的样品判定准确度为95％。

两国大豆产地溯源鉴定模型实施例2

该实施例中，选取具有明确地区的阿根廷大豆样品美国大豆样品来制作标准大豆油样品、标准大豆粉样品、及标准大豆粉样品的水溶性蛋白质，并将这两个不同地区的大豆油样品分别通过液相色谱-四级杆飞行时间质谱仪采集IDA-MS高分辨质谱数据，IDA-MS高分辨质谱数据在PeakView，MasterView软件上定性定量处理分析后，获得所述标准大豆油样品中的甘油三酯化合物标志物观测峰值数据，并将这两个不同地区的大豆粉样品分别通过ICP-MS和/或ICP-AES采集矿质元素含量数据，以及将这两个不同地区的大豆粉样品的水溶性蛋白质分别通过EA-IRMS采集稳定同位素比值数据，将甘油三酯化合物标志物观测峰值数据、矿质元素含量数据和稳定同位素比值数据导入SIMCA 14.1软件中进行正交偏最小二乘法回归分析法处理，从而建立基于IRMS、LC-Q-TOF-MS和多元素分析技术相结合的大豆产地溯源模型。从图9中可以看出阿根廷-美国两国OPLS-DA大豆产地溯源鉴定模型中可显著区分美国和阿根廷大豆样品，为了进一步验证两国压榨大豆油产地溯源鉴定模型的判定准确度，选择美国大豆样品和阿根廷大豆样品分别进行压榨和磨粉提取水溶性蛋白质后进行模型盲样验证，验证结果表明，阿根廷来源的大豆油样品判定准确度为95％，美国来源的样品判定准确度为95％。

两国大豆产地溯源鉴定模型实施例3

该实施例中，选取具有明确地区的加拿大大豆样品美国大豆样品来制作标准大豆油样品、标准大豆粉样品、及标准大豆粉样品的水溶性蛋白质，并将这两个不同地区的大豆油样品分别通过液相色谱-四级杆飞行时间质谱仪采集IDA-MS高分辨质谱数据，IDA-MS高分辨质谱数据在PeakView，MasterView软件上定性定量处理分析后，获得所述标准大豆油样品中的甘油三酯化合物标志物观测峰值数据，并将这两个不同地区的大豆粉样品分别通过ICP-MS和/或ICP-AES采集矿质元素含量数据，以及将这两个不同地区的大豆粉样品的水溶性蛋白质分别通过EA-IRMS采集稳定同位素比值数据，将甘油三酯化合物标志物观测峰值数据、矿质元素含量数据和稳定同位素比值数据导入SIMCA 14.1软件中进行正交偏最小二乘法回归分析法处理，从而建立基于IRMS、LC-Q-TOF-MS和多元素分析技术相结合的大豆产地溯源模型。从图10中可以看出加拿大-美国两国OPLS-DA大豆产地溯源鉴定模型中可显著区分美国和加拿大大豆样品，为了进一步验证两国压榨大豆油产地溯源鉴定模型的判定准确度，选择美国大豆样品和加拿大大豆样品分别进行压榨和磨粉提取水溶性蛋白质后进行模型盲样验证，验证结果表明，加拿大来源的大豆油样品判定准确度为90％，美国来源的样品判定准确度为90％。

两国大豆产地溯源鉴定模型实施例4

该实施例中，选取具有明确地区的俄罗斯、中国大豆样品美国大豆样品来制作标准大豆油样品、标准大豆粉样品、及标准大豆粉样品的水溶性蛋白质，并将两个不同地区的大豆油样品分别通过液相色谱-四级杆飞行时间质谱仪采集IDA-MS高分辨质谱数据，IDA-MS 高分辨质谱数据在PeakView，MasterView软件上定性定量处理分析后，获得所述标准大豆油样品中的甘油三酯化合物标志物观测峰值数据，并将两个不同地区的大豆粉样品分别通过ICP-MS和/或ICP-AES采集矿质元素含量数据，以及将两个不同地区的大豆粉样品的水溶性蛋白质分别通过EA-IRMS采集稳定同位素比值数据，将甘油三酯化合物标志物观测峰值数据、矿质元素含量数据和稳定同位素比值数据导入SIMCA 14.1软件中进行正交偏最小二乘法回归分析法处理，从而建立基于IRMS、LC-Q-TOF-MS和多元素分析技术相结合的俄罗斯-美国、中国-美国这两种大豆产地溯源模型。从图11和12中可以看出俄罗斯-美国、中国-美国两国OPLS-DA大豆产地溯源鉴定模型中可显著区分美国和俄罗斯，美国和中国大豆样品，为了进一步验证两国压榨大豆油产地溯源鉴定模型的判定准确度，选择美国大豆样品、俄罗斯大豆样品和中国大豆样品分别进行压榨和磨粉及提取水溶性蛋白质后进行模型盲样验证，验证结果表明，俄罗斯、中国来源的大豆油样品判定准确度为100％，美国来源的样品判定准确度为95％以上。

由上述四个两国大豆产地溯源鉴定模型实施例和多国鉴定模型可以看出，两国大豆产地溯源鉴定模型的溯源准确度更高。

对比实施例1

该对比实施例中，为验证大豆不同类型的靶向目标物对建立基于稳定同位素比值分析的OPLS-DA鉴定模型的影响，本对比实施例选取具有明确产地的大豆样品分别经过研磨后得到大豆粉样品作为靶向目标物，用EA-IRMSS采集大豆粉样品中δ²H、δ¹⁸O、δ¹⁵N和δ¹³C数据，利用正交偏最小二乘法回归分析法建立多国OPLS-DA大豆产地溯源鉴定模型，建立的模型如图13所示。

对比实施例2

该对比实施例中，为验证大豆不同类型的靶向目标物对建立基于稳定同位素比值分析的OPLS-DA鉴定模型的影响，本对比实施例选取具有明确产地的大豆样品分别经过压榨后得到大豆油样品作为靶向目标物，用EA-IRMSS采集大豆油样品中δ²H、δ¹⁸O、δ¹⁵N和δ¹³C数据，利用正交偏最小二乘法回归分析法建立多国OPLS-DA大豆产地溯源鉴定模型，建立的模型如图14所示。

如图7、13、14所示，对比分析多国OPLS-DA大豆产地溯源鉴定模型实施例1中基于标准大豆粉样品的水溶性蛋白质的δ²H、δ¹⁸O、δ¹⁵N和δ¹³C构建的OPLS-DA多国鉴定模型，和对比实施例1和对比实施例2的对比，来区分不同靶向目标物对稳定同位素溯源鉴定模型的影响。对比结果显示，相同样品，大豆粉样品、大豆油样品及水溶性蛋白质的同位素鉴定模型可以基本区分美国和巴西两地大豆及大豆油样品，但是对其它国家样品与美国样品的区分上，选取大豆水溶性蛋白作为靶向目标物明显优于大豆粉样品和大豆油样品。为此，单独选择大豆粉样品的水溶性蛋白质的δ²H、δ¹⁸O、δ¹⁵N和δ¹³C来建立基于 MALDI-TOF/TOF和稳定同位素技术相结合的大豆产地溯源，能够对巴西、美国、阿根廷、乌拉圭及加拿大等国家大豆样品产地明显区分开。

对比实施例3

该对比实施例中，为验证不同稀释溶剂对大豆油样品的基于甘油三酯脂质组学产地溯源的鉴定模型聚类效应的影响，本对比实施例选取具有明确产地的大豆样品分别经过压榨得到大豆油样品，用丙酮作为稀释溶剂稀释大豆油样品，然后用液相色谱-四级杆飞行时间质谱仪LC-Q-TOF-MS采集大豆油样品中甘油三酯化合物的IDA-MS高分辨质谱数据，利用主成分分析PCA、正交偏最小二乘法回归分析法这两种多元变量统计分析模型处理甘油三酯高分辨质谱数据。

对比实施例4

该对比实施例中，选取具有明确产地的大豆样品分别经过压榨得到大豆油样品，用异丙醇三酯作为稀释溶剂稀释大豆油样品，然后用液相色谱-四级杆飞行时间质谱仪LC-Q-TOF-MS采集大豆油样品中甘油三酯化合物的IDA-MS高分辨质谱数据，利用主成分分析PCA、正交偏最小二乘法回归分析法OPLS-DA这两种多元变量统计分析模型处理甘油三酯高分辨质谱数据。

对比分析多国OPLS-DA大豆产地溯源鉴定模型实施例1、对比实施例3和对比实施例 4不同稀释溶剂对大豆油甘油三酯的产地溯源鉴定的影响。如图15所示，由构建丙酮(acetone)、甲醇与乙酸乙酯混合溶剂(EA-MEOH)、异丙醇(Isopropanol)等溶剂稀释大豆油样品后的PCA模型与OPLS-DA模型可知，丙酮、甲醇与乙酸乙酯溶剂稀释大豆油样品，样品聚集效应具有一定同步性，而异丙醇稀释大豆油样品与其他两种溶剂稀释方法的聚集效应显著不同。尤其是OPLS-DA模型分析结果显示，相同样品，异丙醇稀释大豆油样品聚集区域显著区分甲醇与乙酸乙酯混合溶液及异丙醇稀释大豆油样品。为了进一步考证不同溶剂对大豆油中甘油三酯脂质组学分析的影响，如图11所示构建PCA-class模型，判别三种溶剂稀释的大豆油样品的产地溯源鉴定结果优劣问题。

异丙醇稀释大豆油溯源鉴定模型仅仅可区分压榨大豆油与成品大豆油差异，不能显著鉴定不同国家来源的大豆油样品，甲醇与乙酸乙酯混合溶液稀释大豆油样品可无差别鉴定压榨大豆油及成品大豆油的样品来源，丙酮稀释大豆油样品也可鉴定大豆油样品的产地来源，但是样品聚集的分散度比较高，以及考虑到丙酮对液相色谱柱的影响比较大，不适合用做大豆油样品的稀释溶剂，为此选择甲醇与乙酸乙酯混合溶液作为大豆油样品的稀释剂明显优于其余两种稀释溶剂。

以上对本发明所提供的一种基于IRMS、LC-Q-TOF-MS和多元素分析技术相结合的大豆产地溯源识别方法的实施例进行了详细阐述。本文中应用了具体个例对本发明的原理及实施方式进行了阐述，以上实施例的说明只是用于帮助理解本发明的核心思想。应当指出，对于本技术领域的普通技术人员来说，在不脱离本发明的原理的前提下，还可以本发明进行若干改进和修饰，这些改进和修饰也落入本发明权利要求的保护范围内。