QuEChERS-SPE-超高效液相色谱-串联质谱法测定调味面制品中的12种真菌毒素

徐子婷; 郝莉花; 马静; 卫润鑫; 李晓明; 范雯婷; 路书彦

doi:10.13386/j.issn1002-0306.2021070203

QuEChERS-SPE-超高效液相色谱-串联质谱法测定调味面制品中的12种真菌毒素

河南省产品质量监督检验院，河南郑州 450000

基金项目: 河南省市场监督管理局科技计划项目（2020sj19）；河南省重点研发与推广专项（科技攻关）项目（212102410005）。

详细信息

作者简介:
徐子婷（1986−），女，硕士，工程师，研究方向：食品质量与安全，E-mail：xztwin10@126.com

中图分类号: TS217.1
计量
- 文章访问数: 0183
- HTML全文浏览量: 095
- PDF下载量: 014
出版历程
- 收稿日期: 2021-07-18
- 网络出版日期: 2022-01-17
- 刊出日期: 2022-03-14

Determination of 12 Mycotoxins in Seasoning Flour Products by QuEChERS-SPE-Ultra Performance Liquid Chromatography-Tandem Mass Spectrometry

Henan Product Quality Supervision and Inspection Institute, Zhengzhou 450000, China

摘要

摘要: 建立了QuEChERS-SPE-超高效液相色谱-串联质谱法检测调味面制品中12种真菌毒素残留量的方法。样品由甲酸-乙腈-水溶液提取，加入氯化钠和无水硫酸镁盐析分层，取部分上层清液经Prime HLB小柱萃取净化，净化后的液体再加入无水硫酸镁、乙二胺-N-丙基硅烷（PSA）、C₁₈和中性氧化铝粉末（Al-N）净化；以Waters CORTECS C₁₈+色谱柱分离，梯度洗脱目标化合物，质谱分析采用电喷雾离子源（ESI）正离子模式和多反应监测方式（MRM），以阴性样品为基质空白，匹配标准曲线，外标法定量。研究结果：在浓度0.3～1000 ng/mL范围内，此12种真菌毒素线性关系良好，决定系数（R²）大于0.99，方法检出限为0.10～3.0 μg/kg，定量限为0.30～10.0 μg/kg；取低、中、高三个添加水平下，12种真菌毒素的平均加标回收率为84.2%～97.2%，相对标准偏差为3.25%～8.77%（n=6）。该方法具有回收率高、准确、快速、灵敏度高等特点，可同时对调味面制品中12种真菌毒素残留量进行测定。
- 超高效液相色谱-串联质谱 /
- QuEChERS /
- 固相萃取 /
- 调味面制品 /
- 真菌毒素
Abstract: A method was developed for the residual quantity determination of 12 mycotoxins in seasoning flour products by QuEChERS-SPE-ultra performance liquid chromatography-tandem mass spectrometry. The samples were extracted with formic acid acetonitrile aqueous solution and salted out by adding sodium chloride and anhydrous magnesium sulfate. The supernatant was taken and purified by Prime HLB column extraction. The purified liquid was further cleaned up by adding anhydrous magnesium sulfate, PSA, C₁₈ and Al-N. The target compound, which was separated with Waters CORTECS C₁₈+ column, was eluted by gradient elution. The electrospray ion source (ESI) positive ion mode and multiple reaction monitoring mode (MRM) were used for mass spectrometry. Take negative samples as ground substance whitespace along with matching standard curve. The quantification was achieved by external standard method. The results showed that within the scope of 0.3~1000 ng/mL, the linear relationships of 12 mycotoxins were very good, the determination coefficient (R²) of which were higher than 0.99. The limits of detection were 0.10~3.0 μg/kg, and the limits of quantity were 0.30~10.0 μg/kg. Under the adding levels of high and medium as well as low, the average adding standard recoveries of 12 mycotoxins were 84.2%~97.2% and the relative standard deviations (RSDs) were 3.25%~8.77% (n=6). This method with high recovery rate and high sensitivity was rapid and accurate, which would be used to simultaneously determinate the 12 mycotoxins in seasoning flour products.
- ultra performance liquid chromatography-tandem mass spectrometry /
- QuEChERS /
- solid phase extraction /
- seasoning flour products /
- mycotoxins

HTML全文

辣条是一种以小麦粉和/或其他谷物粉等为主要原料，添加或不添加辅料，经配料、熟制、挤压成型、调味、包装等步骤制作而成的调味面制品^[1]。调味面制品以独特的口感，在国内从儿童到老人不乏受众，就连国外也有很多爱好者，可谓是“史上最牛零食”^[2]。资料显示，2009年调味面制品行业的年生产总值为100多亿元，市场销售额为170多亿元^[3]，到2016年，调味面制品行业年生产总值已增至330亿元，市场销售额增至510亿元。8年间，该行业的生产总值增长230%，销售额增长200%^[4]。而到2018年底，2019年初调味面制品行业年产值规模已高达580多亿元，并且这两年调味面制品生产总值、销售额仍在持续增长，已然成为全国方便休闲食品产业的新经济增长点。

调味面制品的主要原料为粮谷粉类，粮谷在生长、贮藏及运输过程中都易滋生真菌毒素^[5]，且使用受到真菌毒素侵染的粮谷粉加工制成的调味面制品，在市场上很难被消费者所发现。污染粮谷类及粮谷类制品的真菌毒素主要有玉米赤霉烯酮（ZEN）、脱氧雪腐镰刀菌烯醇（DON）、3-乙酰脱氧雪腐镰刀菌烯醇（3-AcDON）、赭曲霉毒素A（OTA）和15-乙酰脱氧雪腐镰刀菌烯醇（15-AcDON）等。同时在调味面制品的生产过程中，需要使用食用油，而植物油和油料种子在生产和储存期间极易产生多种真菌或受到多种真菌毒素的污染。脱氧雪腐镰刀菌烯醇（DON）、伏马毒素（FB）类、黄曲霉（AFT）类、玉米赤霉烯酮（ZEN）等真菌毒素都极易造成食用油及其制品的污染^[6-7]。目前发现这些易引起粮谷、植物油污染的真菌毒素会对人体健康产生严重威胁，过量摄入受污染的食品会导致肝肾的功能性损坏、致癌、致畸，或诱发免疫抑制性疾病^[8-9]。

根据相关文献报道，真菌毒素的检测方法主要有薄层色谱法（TLC）^[10]、酶联免疫吸附法^[11]、气相色谱法（GC）^[12]、液相色谱法（HPLC）^[13-14]及液相色谱-串联质谱法（HPLC-MS/MS）^[15-16]；其中酶联免疫法定性不稳定，易出现假阳性，仅适用于快速检测；气相色谱-质谱联用法样品检测前常需要衍生化处理，在一定程度上受到限制^[17]；液相色谱-质谱联用法拥有液相色谱高效的分离技术以及质谱高分辨率的检测能力，具有准确性好、灵敏度高、特异性和选择性强的特点，可多种组分同时检测，是现开发在复杂基体样品中同时测定多种真菌毒素的良好方法。目前我国现行的真菌毒素国标检测方法主要针对单一类或少数几种真菌毒素的检测，且前处理方法互相独立，检测效率低下。而高通量真菌毒素检测方法研究基质多为粮谷与饲料，此外关于植物油、茶叶、方便面、香辛料中真菌毒素的检测也有相应研究^[18-21]，但是基质相对复杂的调味面制品中多真菌毒素的高通量检测方法目前还没有相关的报道。

同时，高通量检测方法的广泛应用，还需要有效的样品前处理方法为基础，以保证最后的分析结果^[22]。目前真菌毒素检测常用的提取方法有加速溶剂萃取法、固相萃取法、免疫亲和柱法及超临界萃取法，这些提取方法常用于单一类真菌毒素的提取。而QuEChERS技术具有待测物回收率高、应用范围广、分析速度快、污染小和操作简单等优势，已广泛用于动植物源性食品中兽药、农药残留等的测定中。将QuEChERS和固相萃取（SPE）联合使用，能很好的减弱样品的基质效应。为此，本研究采用QuEChERS与SPE技术结合作为前处理方法，并对该方法进行改进和优化，在此基础上，建立了采用UPLC-MS/MS法同时测定调味面制品中12种真菌毒素的高通量检测和确定方法。

1. 材料与方法

1.1 材料与仪器

40份调味面制品　20份购自郑州市各大超市、小食杂店，20份购自京东、淘宝、拼多多等网站；12种真菌毒素标准品：黄曲霉毒素B₁（Aflatoxin B₁，AFB₁）、黄曲霉毒素B₂（Aflatoxin B₂，AFB₂）、黄曲霉毒素G₁（Aflatoxin G₁，AFG₁）、黄曲霉毒素G₂（Aflatoxin G₂，AFG₂）、玉米赤酶烯酮（Zearalenone，ZEN）、脱氧雪腐镰刀菌烯醇（Deoxynival fusenol，DON）、3-乙酰基脱氧雪腐镰刀菌烯醇（3-Acetyl-deoxynival fusenol，3-AcDON）、15-乙酰基脱氧雪腐镰刀菌烯醇（15-Acetyl-deoxynival fusenol，15-AcDON）、T-2毒素（T-2 toxin，T-2）、赭曲霉毒素A（Ochratoxin A，OTA）、伏马毒素B₁（Fumatoxin B₁，FB₁）、伏马毒素B₂（Fumatoxin B₂，FB₂）　均购自上海安谱实验科技股份有限公司。

PRiME HLB固相萃取柱（3 cc，150 mg）、HLB固相萃取柱（6 cc，50 mg）、I-Class-TQS超高效液相色谱-串联质谱仪　美国Waters公司；PSA、C₁₈、中性氧化铝净化填料　月旭科技（上海）股份有限公司；氯化钠、无水硫酸镁　分析纯，天津市恒兴试剂有限公司；甲酸、乙酸铵　MS级，德国赛默飞世尔公司；乙腈　MS级，美国霍尼韦尔公司；VeloCity 14R Pro离心机　英国达美柯公司；MS-TS分析天平　美国梅特勒-托利多公司；Milli-Q纯水仪　美国Millipore公司；HD-2500多管涡旋混合器　杭州佑宁仪器有限公司；5982-911012孔真空固相萃取仪　美国安捷伦公司；N-Evap氮吹仪　上海安谱实验科技股份有限公司。

1.2 实验方法

1.2.1 样品预处理

提取：在50 mL离心管中称取2 g经粉碎机粉碎的均匀试样，加入20 mL 甲酸-乙腈-水（5:45:50，v/v），涡旋振荡10 min，加入4 g MgSO₄，1.5 g NaCl，用手大力振摇1 min。离心（4000 r/min，5 min）后，取部分上清液待净化。

净化：取萃取柱PRIME HLB（3 cc，150 mg）放在12孔固相萃取仪上，无需活化，取0.5 mL上清液预过柱，弃去滤液。装入收集管，移取2 mL上清液过萃取柱并收集滤液，定量取1 mL滤液转移到2 mL离心管（提前装有150 mg MgSO₄，50 mg PSA，30 mg C₁₈，30 mg Al-N）中，涡旋1 min，离心（13000 r/min，3 min）后取500 μL上清液，氮吹后复溶于500 μL乙腈-水（90:10，v/v）中，经0.22 μm膜过滤后，供超高效液相色谱-串联质谱仪分析。

1.2.2 基质匹配标准溶液的制备

分别移取适量12种真菌毒素单标准溶液，用乙腈稀释并定容，配制成混合标准溶液，于−18 ℃储存。称取空白调味面制品样品，按1.2.1方法操作，得到空白基质溶液，用空白基质溶液稀释混合标准溶液配置成一系列混合基质匹配标准溶液，且现用现配。

1.2.3 色谱条件

UPLC色谱柱：CORTECS C₁₈+（2.1×100 mm，粒径1.6 μm）；色谱柱温：40 ℃；流动相A相：0.1%（v/v）甲酸水溶液（含5 mmol/L乙酸铵），B相：0.1%（v/v）甲酸乙腈溶液；流速：0.3 mL/min；进样体积：6 μL，UPLC色谱梯度洗脱条件见表1。

表 1 液相色谱梯度洗脱条件

Table 1. HPLC gradient elution conditions

步骤	时间（min）	流速（mL/min）	流动相体积分数（%）
步骤	时间（min）	流速（mL/min）	A	B
1	0	0.3	90.0	10.0
2	0.50	0.3	90.0	10.0
3	4.00	0.3	40.0	60.0
4	6.00	0.3	40.0	60.0
5	6.10	0.3	10.0	90.0
6	9.00	0.3	10.0	90.0
7	9.10	0.3	90.0	10.0
8	12.00	0.3	90.0	10.0

下载: 导出CSV

| 显示表格

1.2.4 质谱条件

电离源：采用电喷雾电离源（ESI）正模式扫描；调整离子源温度：500 ℃；毛细管电压：1.50 kV；MS采用多反应监测方式（MRM）；12种化合物的质谱采集参照表2。

表 2 12种真菌毒素的串联质谱参数

Table 2. Tandem mass spectrometry parameters of 12 mycotoxins

序号	目标物	保留时间（min）	监测离子对（m·z⁻¹）	锥孔电压（V）	碰撞能量（eV）
1	黄曲霉毒素B₁ AFB₁	3.94	313.2>285.1^*	40	24
1	黄曲霉毒素B₁ AFB₁	3.94	313.2>241.1	40	36
2	黄曲霉毒素B₂ AFB₂	3.75	315.2>287.1^*	40	30
2	黄曲霉毒素B₂ AFB₂	3.75	315.2>259.1	40	26
3	黄曲霉毒素G₁ AFG₁	3.75	329.2>243.1^*	35	25
3	黄曲霉毒素G₁ AFG₁	3.75	329.2>283.1	35	25
4	黄曲霉毒素G₂ AFG₂	3.55	331.2>245.1^*	40	30
4	黄曲霉毒素G₂ AFG₂	3.55	331.2>257.1	40	30
5	玉米赤霉烯酮 ZEN	5.11	319.2>187^*	30	19
5	玉米赤霉烯酮 ZEN	5.11	319.2>185	30	23
6	脱氧雪腐镰刀菌烯醇 DON	2.10	297.1>249.1^*	30	10
6	脱氧雪腐镰刀菌烯醇 DON	2.10	297.1>231.1	30	13
7	3-乙酰基脱氧雪腐镰刀菌烯醇 3-AcDON	3.11	339.2>231.1^*	30	12
7	3-乙酰基脱氧雪腐镰刀菌烯醇 3-AcDON	3.11	339.2>137	30	15
8	15-乙酰基脱氧雪腐镰刀菌烯醇 15-AcDON	3.11	339.2>279.1^*	30	10
8	15-乙酰基脱氧雪腐镰刀菌烯醇 15-AcDON	3.11	339.2>261.1	30	12
9	T-2毒素 T-2	4.90	467.3>305.1^*	30	9
9	T-2毒素 T-2	4.90	467.3>245.1	30	9
10	赭曲霉毒素A OTA	5.06	404.2>239.1^*	30	19
10	赭曲霉毒素A OTA	5.06	404.2>358.2	30	14
11	伏马毒素B₁ FB₁	3.61	722.4>334.2^*	47	40
11	伏马毒素B₁ FB₁	3.61	722.4>352.2	47	38
12	伏马毒素B₂ FB₂	4.09	706.4>336.2^*	47	40
12	伏马毒素B₂ FB₂	4.09	706.4>318.2	47	40
注：*为定量离子对。

下载: 导出CSV

| 显示表格

2. 结果与分析

2.1 流动相的选择

本研究对测定真菌毒素的流动相进行考察，发现真菌毒素类化合物易溶于甲醇或乙腈，其中几种毒素类化合物中含有-OCOCH₃乙酰氧基团，而乙酰氧类化合物在甲醇溶剂中不稳定，因此当流动相选用乙腈:水体系时，该类化合物在色谱柱上的洗脱效果和正模式下的响应强度均好于甲醇:水体系。在乙腈:水体系中尝试加入适量甲酸，发现灵敏度得到显著提高，这应该是增加了毒素类化合物在电喷雾电离时的离子化效率。当流动相选择为0.1%（v/v）甲酸水：0.1%（v/v）甲酸乙腈以初始比例90:10进行梯度洗脱时，观察到AFG2峰形拖尾严重，加入少量乙酸铵后，这种拖尾峰的情况得到有效改善。所以本研究最终选择含5 mmol/L乙酸铵的0.1%（v/v）甲酸水溶液和0.1%（v/v）甲酸乙腈溶液作为流动相，其质谱图见图1。

下载: 全尺寸图片幻灯片

图 1 12种真菌毒素质谱图

Figure 1. Mass spectrogram of 12 mycotoxins

下载: 全尺寸图片幻灯片

2.2 样品前处理条件的优化

2.2.1 样品提取溶剂的优化

本研究探究了不同提取溶剂对真菌毒素提取效果的影响，采用50%甲醇水溶液、50%乙腈水溶液和添加5%甲酸的50%乙腈水溶液对调味面制品中的12种真菌毒素进行提取，提取后加入4 g MgSO₄，1.5 g NaCl，用手大力振摇1 min。离心（4000 r/min，5 min）后，取上层有机相过0.22 μm膜，上机进行测定。不同提取溶剂对真菌毒素的回收率影响见图2。

图 2 不同提取溶剂对回收率的影响

Figure 2. Effects of different extraction solvents on recovery rate

下载: 全尺寸图片幻灯片

由图2可知，12种真菌毒素中，AFB₂、3-ADON、OTA、FB₁、FB₂这5种真菌毒素在50%甲醇水溶液的回收率和50%乙腈水溶液的回收率差异性并不显著（P>0.05），剩余7种真菌毒素的回收率有显著性差异。这是由于甲醇、乙腈作为2种通用型提取溶剂，乙腈较甲醇具有更好的选择性，可以避免提取更多脂肪、色素等杂质，且乙腈的提取效率更好，因此本研究采用乙腈作为提取溶剂^[23]。

从图2中还可看出，添加5%甲酸50%乙腈水溶液，12种真菌毒素回收率均为最高。分析其原因，因为本研究待测真菌毒素中伏马毒素、赭曲霉毒素具有羧基集团，亲水性较强，使用高有机相体系会导致其回收率部分损失，故需要使用水含量较高的提取溶剂进行提取。但提取液中的高比例的水会大量溶出基质中的杂质，造成方法灵敏度降低，从而部分毒素的回收率也会显著下降^[24]。同时，像DON、ZEN、OTA、FB₁、FB₂对酸敏感，降低萃取液的pH能提高稳定性，增加萃取效率，因此为保证各种真菌毒素的提取效果，还需在提取液中添加辅助试剂（乙酸、甲酸等），以增强联合提取的效果^[11]，而甲酸、乙酸效果相近，本研究选用甲酸。

2.2.2 固相萃取柱的优化

由于调味面制品中有大量亲脂性杂质、中性盐的存在，会降低色谱柱的寿命，同时增大基质效应，干扰目标化合物的分析；若选择专一性较强的萃取小柱，虽然对测定某种或某类化合物的效果较好，但对多组分中性质差异较大的样品，会对某些组分的吸附性较差而影响回收率^[25]。

HLB为反相固相萃取柱，对目标物进行吸附，样品中的杂质将在淋洗步骤中除去，反相萃取柱的上样液一般为大比例的水相，本研究需要把乙腈蒸发并置换成水相，才能使目标物在柱子上保留，增加了前处理时间；而Prime-HLB柱选择吸附杂质，对目标物没有吸附，作为聚合物柱不需要活化，可直接萃取，简化了过柱子流程^[26]；通过考察3种加标水平下多种真菌毒素的回收率发现，Prime-HLB柱均优于HLB柱，特别是ZEN经HLB净化后回收率稳定性较差，无法准确定量。造成回收率损失的原因是在HLB这类柱子的操作中，目标物的吸附和选择淋洗液进行杂质淋洗都容易丢失化合物，且ZEN在酸性溶剂净化中才能获得较高的回收率^[27]，详见图3~图5。

图 3 低浓度水平（0.6～20 μg/kg）加标回收率

Figure 3. Recovery rate at low concentration level (0.6~20 μg/kg)

下载: 全尺寸图片幻灯片

图 4 中浓度水平（1.5～50 μg/kg）加标回收率

Figure 4. Recovery rate at medium concentration level (1.5~50 μg/kg)

下载: 全尺寸图片幻灯片

图 5 高浓度水平（3.0～100 μg/kg）加标回收率

Figure 5. Recovery rate at high concentration level (3.0~100 μg/kg)

下载: 全尺寸图片幻灯片

2.2.3 净化条件的优化

调味面制品中存在大量干扰物如色素、脂质等，因此净化步骤对定性、定量结果，尤其是定量的准确性非常重要。QuEChERS法常用净化填料有：石墨化碳黑（GCB）、十八烷基硅烷键合硅胶（C₁₈）、中性氧化铝（Al-N）和乙二胺-N-丙基硅烷（PSA）等。据文献报道，GCB对具有平面分子结构的化合物有显著的吸附效果^[28]，可以脱除某些色素，同时它也能脱除具有平面结构的毒素，如AFT和ZEN等^[29]；C₁₈作为非极性吸附剂能够吸附基质中非极性物质，在去除样品中的脂质、糖类等有机化合物起到有效作用，但也会使得脂溶性物质的检测结果变差，且单独使用净化能力有限；PSA是一种碱性吸附剂，在结构上有两个氨基，可以通过氢键结合去除碳水化合物、有机酸和少量色素，但它也会导致酸性化合物的检测结果减低^[22]；中性氧化铝Al-N可以吸附芳香烃和脂肪胺等富电化合物。因此本研究在选取QuEChERS净化试剂时，分别选取了A组：50 mg PSA，30 mg C₁₈，30 mg Al-N，150 mg MgSO₄；B组：50 mg PSA，150 mg MgSO₄；C组：30 mg C₁₈，150 mg MgSO₄；D组：30 mg Al-N，150 mg MgSO₄。通过比较经PSA、C₁₈和Al-N分别与MgSO₄组合净化后的样品提取液净化结果发现，对于12种真菌毒素，B组和D组回收率较差。A组和C组对AFB₁、AFB₂、DON、3-ADON和FB₂这5种毒素回收率差异性并不显著（P>0.05），对AFG₁、AFG₂、ZEN、15-ADON、T-2、OTA和FB₁这7种毒素回收率差异性明显，且A组回收率普遍高于C组，表明A组不仅能起到良好的净化作用，还能得到较高的回收率，最终确定50 mg PSA、30 mg C₁₈、30 mg Al-N、150 mg MgSO₄混合填料作为QuEChERS净化吸附剂，详见图6。

图 6 不同净化剂对回收率的影响

Figure 6. Effects of diffent purifiers on recovery rate

下载: 全尺寸图片幻灯片

2.3 方法验证

2.3.1 标准曲线与检出限

根据12种化合物的响应强弱显示，在空白调味面制品基质溶液中添加系列浓度的真菌毒素，配制标准溶液。用超高效液相色谱-串联质谱仪进行检测，纵坐标选为峰面积（Y），横坐标为对应的质量浓度（X，ng/mL），绘制标准曲线，参见表3。

表 3 12种真菌毒素的线性范围、线性方程、决定系数、检出限和定量限

Table 3. Linear ranges, linear equations, determination coefficients (R²), LODs and LOQs of 12 mycotoxins

序号	真菌毒素	线性方程	线性范围（ng/mL）	决定系数R²	检出限（μg/kg）	定量限（μg/kg）
1	AFB₁	Y=33887X−13260.8	0.3~100	0.9989	0.10	0.30
2	AFB₂	Y=18634.9X−5523.12	0.3~100	0.9989	0.10	0.30
3	AFG₁	Y=24469.6X−6523.41	0.3~100	0.9976	0.10	0.30
4	AFG₂	Y=6539.12X−9929.47	0.3~100	0.9953	0.10	0.30
5	ZEN	Y=5515.82X−23299.1	0.5~100	0.9973	0.50	1.0
6	DON	Y=383.675X−4922.25	1.0~100	0.9957	0.50	1.0
7	3-ADON	Y=1994.44X−8549.09	2.0~200	0.9938	1.0	2.0
8	15-ADON	Y=707.06X−4070.51	5.0~500	0.9961	2.0	5.0
9	T-2	Y=40.8983X−176.602	1.0~100	0.9902	1.0	2.0
10	OTA	Y=59309.3X−107169	1.0~100	0.9957	0.50	1.0
11	FB₁	Y=2691.25X−64092.1	10~100	0.9926	3.0	10.0
12	FB₂	Y=5724.1X−162588	10~1000	0.9906	3.0	10.0

下载: 导出CSV

| 显示表格

由表3可知，这12种化合物在0.3～1000 ng/mL质量浓度范围内能够呈现良好的线性关系，其决定系数（R²）均不低于0.99。根据添加回收试验，以3倍信噪比（S/N）确定该12种真菌毒素的检出限（LOD）为0.10～3.0 μg/kg，以10倍信噪比（S/N）确定该12种真菌毒素的定量限（LOQ）为0.30～10.0 μg/kg。

2.3.2 回收率与精密度

取空白调味面制品样品，以低、中、高（2倍定量限、5倍定量限、10倍定量限）三种浓度，每个加标水平选取6个平行样，添加12种真菌毒素标准品后按照1.2.1进行试验，12种真菌毒素的平均回收率为84.2%～97.2%，相对标准偏差为3.25%～8.77%，符合GB/T 27417-2017《合格评定化学分析方法确认和验证指南》^[30]标准中关于定性定量要求，同时也符合痕量分析标准含量低于百万分之一的要求。回收率与精密度数据见表4。

表 4 空白样品中12种真菌毒素的加标回收率与相对标准偏差

Table 4. Recoveries and RSDs for 12 mycotoxins

真菌毒素	添加水平（μg/kg）	平均回收率（n=6，%）	相对标准偏差（n=6，%）
AFB₁	0.6	92.2	6.71
	1.5	95.4	4.28
	3	91.8	5.20
AFB₂	0.6	92.6	8.68
	1.5	94.8	5.03
	3	92.8	6.72
AFG₁	0.6	90.8	7.54
	1.5	90.5	5.26
	3	89.8	6.83
AFG₂	0.6	91.8	8.25
	1.5	92.3	5.63
	3	90.6	6.45
ZEN	2	96.6	8.77
	5	97.2	4.96
	10	94.3	5.02
DON	2	85.6	8.70
	5	90.3	6.52
	10	88.2	7.34
3-ADON	4	84.2	6.17
	10	86.7	6.45
	20	84.6	7.03
15-ADON	10	93.6	6.23
	25	95.1	7.54
	50	94.2	7.46
T-2	4	89.9	8.25
	10	92.7	7.12
	20	90.6	4.23
OTA	2	88.1	5.21
	5	92.5	7.54
	10	90.6	6.83
FB₁	20	86.6	3.25
	50	84.8	8.75
	100	86.8	6.34
FB₂	20	88.1	8.42
	50	90.3	6.56
	100	86.7	8.52

下载: 导出CSV

| 显示表格

2.4 实际样品的检测

采用建立的QuEChERS-SPE-超高效液相色谱-串联质谱法对40份调味面制品进行分析，其中3份样品检出黄曲霉毒素B₁，含量分别为0.32、0.62、0.66 μg/kg；2份样品检出脱氧雪腐镰刀菌烯醇，含量分别为2.4、3.1 μg/kg。现行的国家标准GB 2761-2017《食品安全国家标准食品中真菌毒素限量》^[31]并未对调味面制品中的黄曲霉毒素B₁和脱氧雪腐镰刀菌烯醇做出限量规定，可以参考GB 2761-2017植物油脂（花生油、玉米油除外）中黄曲霉毒素B₁限量10 μg/kg、花生油、玉米油黄曲霉毒素B₁限量20 μg/kg；谷物及其制品中脱氧雪腐镰刀菌烯醇限量1000 μg/kg。根据以上标准，本研究所取的40份样品中，虽然部分样品黄曲霉毒素B₁和脱氧雪腐镰刀菌烯醇检出结果呈阳性，但其含量均低于真菌毒素限量标准。

3. 结论

本研究结合QuEChERS-固相萃取SPE前处理方法的优点，建立了UPLC超高效液相色谱-串联质谱法同时检测调味面制品中12种真菌毒素的检测方法。在对实际样品进行测定时，电喷雾电离源ESI⁺模式下，MRM多反应监测扫描，6 min内12种化合物均能出峰，同时本方法还具有较低的检出限和定量限。因此本方法具有快速、准确、回收率高，灵敏度高等特点，适用于调味面制品中12种真菌毒素残留的定性分析和定量分析，补充了调味面制品中多种真菌毒素的定量确证检测方法。

下载: 全尺寸图片幻灯片

图 1 12种真菌毒素质谱图

Figure 1. Mass spectrogram of 12 mycotoxins

下载: 全尺寸图片幻灯片

图 2 不同提取溶剂对回收率的影响

Figure 2. Effects of different extraction solvents on recovery rate

下载: 全尺寸图片幻灯片

图 3 低浓度水平（0.6～20 μg/kg）加标回收率

Figure 3. Recovery rate at low concentration level (0.6~20 μg/kg)

下载: 全尺寸图片幻灯片

图 4 中浓度水平（1.5～50 μg/kg）加标回收率

Figure 4. Recovery rate at medium concentration level (1.5~50 μg/kg)

下载: 全尺寸图片幻灯片

图 5 高浓度水平（3.0～100 μg/kg）加标回收率

Figure 5. Recovery rate at high concentration level (3.0~100 μg/kg)

下载: 全尺寸图片幻灯片

图 6 不同净化剂对回收率的影响

Figure 6. Effects of diffent purifiers on recovery rate

下载: 全尺寸图片幻灯片

表 1 液相色谱梯度洗脱条件

Table 1 HPLC gradient elution conditions

步骤	时间（min）	流速（mL/min）	流动相体积分数（%）
步骤	时间（min）	流速（mL/min）	A	B
1	0	0.3	90.0	10.0
2	0.50	0.3	90.0	10.0
3	4.00	0.3	40.0	60.0
4	6.00	0.3	40.0	60.0
5	6.10	0.3	10.0	90.0
6	9.00	0.3	10.0	90.0
7	9.10	0.3	90.0	10.0
8	12.00	0.3	90.0	10.0

下载: 导出CSV

表 2 12种真菌毒素的串联质谱参数

Table 2 Tandem mass spectrometry parameters of 12 mycotoxins

序号	目标物	保留时间（min）	监测离子对（m·z⁻¹）	锥孔电压（V）	碰撞能量（eV）
1	黄曲霉毒素B₁ AFB₁	3.94	313.2>285.1^*	40	24
1	黄曲霉毒素B₁ AFB₁	3.94	313.2>241.1	40	36
2	黄曲霉毒素B₂ AFB₂	3.75	315.2>287.1^*	40	30
2	黄曲霉毒素B₂ AFB₂	3.75	315.2>259.1	40	26
3	黄曲霉毒素G₁ AFG₁	3.75	329.2>243.1^*	35	25
3	黄曲霉毒素G₁ AFG₁	3.75	329.2>283.1	35	25
4	黄曲霉毒素G₂ AFG₂	3.55	331.2>245.1^*	40	30
4	黄曲霉毒素G₂ AFG₂	3.55	331.2>257.1	40	30
5	玉米赤霉烯酮 ZEN	5.11	319.2>187^*	30	19
5	玉米赤霉烯酮 ZEN	5.11	319.2>185	30	23
6	脱氧雪腐镰刀菌烯醇 DON	2.10	297.1>249.1^*	30	10
6	脱氧雪腐镰刀菌烯醇 DON	2.10	297.1>231.1	30	13
7	3-乙酰基脱氧雪腐镰刀菌烯醇 3-AcDON	3.11	339.2>231.1^*	30	12
7	3-乙酰基脱氧雪腐镰刀菌烯醇 3-AcDON	3.11	339.2>137	30	15
8	15-乙酰基脱氧雪腐镰刀菌烯醇 15-AcDON	3.11	339.2>279.1^*	30	10
8	15-乙酰基脱氧雪腐镰刀菌烯醇 15-AcDON	3.11	339.2>261.1	30	12
9	T-2毒素 T-2	4.90	467.3>305.1^*	30	9
9	T-2毒素 T-2	4.90	467.3>245.1	30	9
10	赭曲霉毒素A OTA	5.06	404.2>239.1^*	30	19
10	赭曲霉毒素A OTA	5.06	404.2>358.2	30	14
11	伏马毒素B₁ FB₁	3.61	722.4>334.2^*	47	40
11	伏马毒素B₁ FB₁	3.61	722.4>352.2	47	38
12	伏马毒素B₂ FB₂	4.09	706.4>336.2^*	47	40
12	伏马毒素B₂ FB₂	4.09	706.4>318.2	47	40
注：*为定量离子对。

下载: 导出CSV

表 3 12种真菌毒素的线性范围、线性方程、决定系数、检出限和定量限

Table 3 Linear ranges, linear equations, determination coefficients (R²), LODs and LOQs of 12 mycotoxins

序号	真菌毒素	线性方程	线性范围（ng/mL）	决定系数R²	检出限（μg/kg）	定量限（μg/kg）
1	AFB₁	Y=33887X−13260.8	0.3~100	0.9989	0.10	0.30
2	AFB₂	Y=18634.9X−5523.12	0.3~100	0.9989	0.10	0.30
3	AFG₁	Y=24469.6X−6523.41	0.3~100	0.9976	0.10	0.30
4	AFG₂	Y=6539.12X−9929.47	0.3~100	0.9953	0.10	0.30
5	ZEN	Y=5515.82X−23299.1	0.5~100	0.9973	0.50	1.0
6	DON	Y=383.675X−4922.25	1.0~100	0.9957	0.50	1.0
7	3-ADON	Y=1994.44X−8549.09	2.0~200	0.9938	1.0	2.0
8	15-ADON	Y=707.06X−4070.51	5.0~500	0.9961	2.0	5.0
9	T-2	Y=40.8983X−176.602	1.0~100	0.9902	1.0	2.0
10	OTA	Y=59309.3X−107169	1.0~100	0.9957	0.50	1.0
11	FB₁	Y=2691.25X−64092.1	10~100	0.9926	3.0	10.0
12	FB₂	Y=5724.1X−162588	10~1000	0.9906	3.0	10.0

下载: 导出CSV

表 4 空白样品中12种真菌毒素的加标回收率与相对标准偏差

Table 4 Recoveries and RSDs for 12 mycotoxins

真菌毒素	添加水平（μg/kg）	平均回收率（n=6，%）	相对标准偏差（n=6，%）
AFB₁	0.6	92.2	6.71
	1.5	95.4	4.28
	3	91.8	5.20
AFB₂	0.6	92.6	8.68
	1.5	94.8	5.03
	3	92.8	6.72
AFG₁	0.6	90.8	7.54
	1.5	90.5	5.26
	3	89.8	6.83
AFG₂	0.6	91.8	8.25
	1.5	92.3	5.63
	3	90.6	6.45
ZEN	2	96.6	8.77
	5	97.2	4.96
	10	94.3	5.02
DON	2	85.6	8.70
	5	90.3	6.52
	10	88.2	7.34
3-ADON	4	84.2	6.17
	10	86.7	6.45
	20	84.6	7.03
15-ADON	10	93.6	6.23
	25	95.1	7.54
	50	94.2	7.46
T-2	4	89.9	8.25
	10	92.7	7.12
	20	90.6	4.23
OTA	2	88.1	5.21
	5	92.5	7.54
	10	90.6	6.83
FB₁	20	86.6	3.25
	50	84.8	8.75
	100	86.8	6.34
FB₂	20	88.1	8.42
	50	90.3	6.56
	100	86.7	8.52

下载: 导出CSV

参考文献(31)

[1]	国家食品药品监督管理总局. 总局办公厅关于印发国家食品安全监督抽检实施细则(2018年版)的通知 (食药监办食监三〔2018〕14号) [Z]. 2018-01-18. State Food and Drug Administration. The general office of the State Food and Drug Administration issued the detailed rules for the implementation of the state food safety supervision and sampling inspection (2018 edition) [Z]. 2018-01-18.
[2]	安虹, 谭炜. 调味面制品产品质量现状调研及分析[J]. 安徽农业科学,2016,44(3):88−89,122. [AN H, TAN W. The quality status investigation and analysis of seasoning flour products[J]. J Anhui Agric Sci,2016,44(3):88−89,122. doi: 10.3969/j.issn.0517-6611.2016.03.033
[3]	湘式挤压糕点等休闲食品质量安全监管工作座谈会召开[N]. 中国食品安全报, 2015-04-14(1). Convene xiang style extrusion cakes and other leisure food quality and safety supervision symposium [N]. China Food Safety News, 2015-04-14(1).
[4]	黄河新闻网. 调味面制品行业加速升级“五毛”辣条渐退出市场[EB/OL]. [2017-07-11]. https://www.sohu.com/a/156265439_183391. Yellow River News Net. Seasoning flour products to accelerate upgrading “50 cent” seasoning flour products gradually withdraw from the market [EB/OL]. [2017-07-11]. https://www.sohu.com/a/156265439_183391.
[5]	马皎洁, 邵兵, 林肖惠,等. 我国部分地区2010年产谷物及其制品中多组分真菌毒素污染状况研究[J]. 中国食品卫生杂志,2011,23(6):481−488. [MA J J, SHAO B, LIN X H, et al. Study on the natural occurrence of multi-mycotoxin in cereal and cereal-based product samples collected from parts of China in 2010[J]. Chinese Journal of Food Hygiene,2011,23(6):481−488.
[6]	马治良, 徐同城, 刘丽娜, 等. 食用油及油料作物中真菌毒素研究进展[J]. 农产品加工,2014(2):45−48. [MA Z L, XU T C, LIU L N, et al. Research progress of mycotoxins in edible oils and oil crops[J]. Farm Products Processing,2014(2):45−48.
[7]	何景, 杨丹. 北京市地区小包装食用油中真菌毒素污染状况调查[J]. 中国油脂,2019,44(6):79−82. [HE J, YANG D. Investigation on contamination of mycotoxins in small package edible oil in Beijing[J]. China Oils and Fats,2019,44(6):79−82. doi: 10.3969/j.issn.1003-7969.2019.06.017
[8]	ALASSANE-KPEMBI I, SCHATZMAYR G, TARANU I, et al. Mycotoxins co-contamination: Methodological aspects and biological relevance of combined toxicity studies[J]. Crit Rev Food Sci Nutr,2017,57(16):3489. doi: 10.1080/10408398.2016.1140632
[9]	MANIZAN A I, OPLATOWSKA-STACHOWIAK M, PIRO-METAYER I, et al. Multi-mycotoxin determination in rice, maize and peanut products most consumed in Côte d’Ivoire by UHPLC-MS/MS[J]. Food Control,2018,87:22. doi: 10.1016/j.foodcont.2017.11.032
[10]	AIZZO I, VEDOYA G, MAURUTTO S, et al. Assessment of toxigenic fungi on Argentinean medicinal herbs[J]. Microbiol Res,2004,159(2):113. doi: 10.1016/j.micres.2004.01.013
[11]	肖娟, 辛苑, 冯莎, 等. 酶联免疫吸附法测定小麦中脱氧雪腐镰刀菌烯醇[J]. 食品安全质量检测学报,2013,4(1):224−228. [XIAO J, XIN Y, FENG S, et al. Determination of deoxynivalenol in wheat by Elisa[J]. Journal of Food Safety & Quality,2013,4(1):224−228.
[12]	MANY, LIANG G, LI A, et al. Recent advances in mycotoxin determination for food monitoring via microchip[J]. Toxins,2017,9(10):324. doi: 10.3390/toxins9100324
[13]	邵珍美. 免疫亲和柱-高效液相色谱法测定牛奶中的呕吐毒物[J]. 食品研究与开发,2016,37(19):172−175. [SHAO Z M. Determination of vomitoxin in milk by high perfonmance liquid chromatography with immunoaffinity column[J]. Food Research and Development,2016,37(19):172−175. doi: 10.3969/j.issn.1005-6521.2016.19.040
[14]	李丽, 叶金, 辛媛媛,等. 高效液相色谱法测定粮食及其制品中赭曲霉毒素A提取条件的优化[J]. 食品安全质量检测学报,2018,9(21):5602−5607. [LI L, YE J, XIN Y Y, et al. Determination of vomitoxin in milk by high perfonmance liquid chromatography with immunoaffinity column[J]. Journal of Food Safety & Quality,2018,9(21):5602−5607. doi: 10.3969/j.issn.2095-0381.2018.21.012
[15]	孟娟, 张娟, 张楠,等. 固相萃取-超高效液相色谱串联质谱法检测粮食及其制品中的玉米赤霉烯酮类真菌毒素[J]. 色谱,2010,28(28):601−607. [MENG J, ZHANG J, ZHANG N, et al. Determination of zearalenone and related mycotoxins in grains and its products by solid-phase extraction coupled with ultra perfonmance liquid chromatography-tandem mass spectrometry[J]. Chinese Journal of Chromatography,2010,28(28):601−607.
[16]	和佳鸳, 蒲彦利, 徐虹,等. UPLC-MS/MS测定婴幼儿谷物辅食中脱氧雪腐镰刀菌烯醇及其衍生物[J]. 中国卫生检验杂志,2015,25(17):2860−2862. [HE J Y, PU Y L, XU H, et al. Determination of dexynivalenol and its derivants in cereal-based complementary foods by UPLC-MS/MS[J]. Chinese Journal of Health Laboratory Technology,2015,25(17):2860−2862.
[17]	肖全伟, 吴文林, 刘玲利. QuEChERS-高效液相色谱-串联质谱法测定小麦中的8种真菌毒素[J]. 四川农业大学学报,2018,36(5):584−591. [XIAO Q W, WU W L, LIU L L. Determination of eight mycotoxins in wheat by QuEChERS-high performance liquid chromatography-tandem mass spectrometry[J]. Journal of Sichuan Agricultural University,2018,36(5):584−591.
[18]	吴宇, 叶金, 张冰,等. 稳定同位素稀释-超高效液相色谱-串联质谱法快速测定植物油中16种真菌毒素[J]. 分析化学,2018(6):975−984. [WU Y, YE J, ZHANG B, et al. Rapid determination of 16 mycotoxins in vegetable oils by stable isotope-High- performance liquid chromatography[J]. Chinese Journal of Analytical Chemistry,2018(6):975−984. doi: 10.11895/j.issn.0253-3820.171315
[19]	胡琳, 师真, 赵丽,等. 液相色谱-串联质谱法同时测定普洱茶中16种真菌毒素[J]. 浙江农业学报,2019,31(10):1700−1708. [HU L, SHI Z, ZHAO L, et al. Simultaneous detection and analysis of 16 kinds of mycotoxins in Pu-erh tea[J]. Acta Agriculturae Zhejiangensis,2019,31(10):1700−1708. doi: 10.3969/j.issn.1004-1524.2019.10.16
[20]	蒋孟圆, 张瑞雨, 赵丽, 等. 方便面6种真菌毒素的监测分析[J]. 环境化学,2020,39(3):828−830. [JIANG M Y, ZHANG R Y, ZHAO L, et al. Monitoring and analysis of 6 mycotoxins in instant noodles[J]. Environmental Chemistry,2020,39(3):828−830.
[21]	黄思瑜, 董宪兵, 周纯洁,等. 超高效液相色谱串联-质谱法同时测定香辛料中7种真菌毒素[J]. 食品工业科技,2019,40(14):270−274,279. [HUANG S Y, DONG X B, ZHOU C J, et al. Determination of 7 mycotoxin contaminants in spices by ultra-high performance liquid chromatography-tandem mass spectrometry[J]. Science and Technology of Food Industry,2019,40(14):270−274,279.
[22]	赵英莲, 张梓琪, 赵鑫, 等. QuEChERS技术在食品真菌毒素检测中的研究进展[J]. 中国酿造,2020,39(1):1−5. [ZHAO Y L, ZHANG Z Q, ZHAO X, et al. Research progress on QuEChERS technology in the detection of mycotoxins in food[J]. China Brewing,2020,39(1):1−5. doi: 10.11882/j.issn.0254-5071.2020.01.001
[23]	王少敏, 杜春晓, 刘贤贤,等. QuEChERS-超高效液相色谱-串联质谱法同时测定三七中26种真菌毒素[J]. 世界中医药,2019,14(4):798−804. [WANG S M, DU C X, LIU X X, et al. Simultaneous determination of 26 mycotoxins in notoginseng radix et rhizoma by QuEChERS-ultra-high-performance liquid chromatography-tandem mass spectrometry[J]. World Chinese Medicine,2019,14(4):798−804. doi: 10.3969/j.issn.1673-7202.2019.04.003
[24]	杨帅, 穆蕾, 杨悠悠, 等. 高效液相色谱-串联质谱法同时测定小麦粉中28种真菌毒素[J]. 中国粮油学报,2020,35(6):159−164. [YANG S, MU L, YANG Y Y, et al. Simultaneous determination of 28 mycotoxins in wheat flour by high-performance liquid chromatography-tandem mass spectrometry[J]. Journal of the Chinese Cereals and Oils Association,2020,35(6):159−164. doi: 10.3969/j.issn.1003-0174.2020.06.024
[25]	钟世欢, 叶佳明, 叶磊海, 等. 超高压液相色谱-串联质谱法测定粮食中多种真菌毒素残留[J]. 农产品加工,2018(12):59−62. [ZHONG S H, YE J M, YE L H, et al. Determination of mycotoxins in cereals by UPLC-MS/MS[J]. Farm Products Processing,2018(12):59−62.
[26]	徐潇颖, 刘柱, 梁晶晶,等. 基于2种固相萃取方法净化的液相色谱-串联质谱法测定粮食制品中多种真菌毒素[J]. 粮食与油脂,2018,31(11):76−80. [XU X Y, LIU Z, LIANG J J, et al. Determination of mycotoxins in grain products by liquid chromatography-tandem mass spectrometry based on two solid phase extraction purification methods[J]. Cereals & Oils,2018,31(11):76−80. doi: 10.3969/j.issn.1008-9578.2018.11.020
[27]	尹青冈. 玉米中玉米赤霉烯酮辐照讲解技术研究[D]. 重庆: 西南大学, 2009. YIN Q G. Study on irradiation technology of zearalenone in maize[D]. Chongqing: Southwest University, 2009.
[28]	Casiane Salete Tibola, Cunha Fernandes. Distribution of Fusarium mycotoxins in wheat milling process[J]. Food Control,2015:91−95.
[29]	GONG X M, WANG H B, ZHANG Y B, et al. Determination of 15 mycotoxins in foods and feeds using high performance liquid chromatography-tandem mass spectrometry with gel permeation chromatography combined QuEChERS purification[J]. J Chromatogr Sep T,2012,3(2):2−8.
[30]	GB/T27417-2017合格评定化学分析方法确认和验证指南[S]. 北京: 中国标准出版社, 2017. GB/T27417-2017 Comformity assessment-Guidance on validation and verification of chemical analytical methods[S]. Beijing: Standards Press of China, 2017.
[31]	GB2761-2017食品安全国家标准食品中真菌毒素限量[S]. 北京: 中国标准出版社, 2017. GB 2761-2017 National food safety standard-Limit of mycotoxins in food[S]. Beijing: Standards Press of China, 2017.

施引文献(14)

期刊类型引用(12)

1.	李敏晶，李丹阳，范梓伊. 响应面优化海星肽-锌螯合物的制备. 农产品加工. 2025(04): 11-14+19 . 百度学术
2.	张迎阳，董亚云，邹平，郑文武. 橄榄仁抗氧化肽的分离鉴定及分子对接. 食品与发酵工业. 2024(05): 189-198 . 百度学术
3.	王扬铎，苏永昌，王晓燕，施文正，刘智禹. 仿刺参不同部位ACE抑制活性分析及活性肽制备工艺优化. 食品工业科技. 2024(10): 187-197 . 本站查看
4.	张薇，李爱琳，李书国. 超声辅助酶解制备大米蛋白黄嘌呤氧化酶活性抑制肽的工艺条件优化. 中国粮油学报. 2024(06): 51-58 . 百度学术
5.	叶灏铎，管晓盛，马凤，毛远辉，孙世利，曹庸，苗建银. 英红九号茶蛋白降尿酸肽的酶解制备及不同分子量组分的活性对比. 现代食品科技. 2023(03): 147-155 . 百度学术
6.	沈畅华，杨娟，张远红，曾晓房. 双酶酶解鸽胸肉工艺优化及抗氧化性评价. 食品与机械. 2023(04): 163-169 . 百度学术
7.	杨宝强，刘生杰. 响应面法优化猴头菇蛋白抗氧化肽工艺及抗运动疲劳活性研究. 粮食与油脂. 2023(05): 142-147 . 百度学术
8.	牛金鸽，吴海玥，马世科，闫忠心，王学江，李婧，胡蓉，祁全青. 响应面优化藏羊皮胶原蛋白肽超声辅助提取工艺及其体内抗氧化活性分析. 食品工业科技. 2023(11): 163-170 . 本站查看
9.	马凤，叶灏铎，夏珍，徐燕，孙世利，曹庸，苗建银. 英红九号茶蛋白ACE抑制肽的制备、氨基酸组成及不同超滤组分的活性评价. 现代食品科技. 2023(07): 237-245 . 百度学术
10.	张艳，江波涛，李莉，王佳辉，王爽，李平，徐婷婷，刘雪松，薛沾枚，吴宪，邱景会. 驴血血红蛋白肽的制备工艺研究. 黑龙江畜牧兽医. 2023(17): 134-138+142 . 百度学术
11.	张恒慧，张志军，陈士国，叶兴乾，张国华. 紫苏粕蛋白抗氧化活性肽的制备、分离纯化及序列鉴定. 中国食品学报. 2023(09): 347-355 . 百度学术
12.	孙雪芹，廖彭莹，刘华宇，韦金锐，刘舒凌，陈俊，林燕媚，李锦清. 中华鳖肉水溶蛋白酶解工艺优化. 食品科技. 2022(11): 121-128 . 百度学术