基于宏基因组测序技术的非靶向筛查方法检测肉类食品中的动物源性成分

基金项目：国家市场监督管理总局科技计划项目（2019MK057）；广东省市场监督管理局科技计划项目（2021ZS03）。

作者简介：丁清龙（1990—），男，河北沧州人，硕士，工程师。研究方向：食品安全检测。

通信作者：周露（1982—），女，河南洛阳人，博士，正高级工程师。研究方向：食品安全检测。E-mail：zhoulu1982@sohu.com。

摘 要：目的：建立肉类食品中动物源性成分非靶向筛查方法。方法：基于宏基因组测序技术建立肉类食品中动物源性成分非靶向筛查方法，将该方法用于模拟样品和实际样品检测，并用现有标准检测方法对实际样品检测结果进行确认。结果：模拟样品检测结果与模拟情况一致；在实际样品中检出与样品名称不一致或样品标签未标识的动物源性成分，且非靶向筛查方法检测结果与现有标准方法检测结果一致。结论：该方法可以快速锁定样品中未知的动物源性成分，检测结果准确可靠，可为政府打击肉类食品掺假提供更为有力的技术支撑。

关键词：宏基因组测序；动物源性成分；非靶向；掺假

Determination of Animal-Derived Ingredients in Meat Products by Non-Targeted Screening Method Based on Metagenomic Sequencing Technology

DING Qinglong， YANG Danting， XIE Aihua， WEI Yun， CHEN Xiufen， ZHOU Lu*

（Guangdong Institute of Food Inspection， Guangzhou 510435， China）

Abstract： Objective： To establish the non-targeted screening method for animal-derived ingredients in meat products. Method： The non-targeted screening method for animal-derived ingredients in meat products was established based on metagenomic sequencing technology. The method was applied to the determination of simulated samples and actual samples， and the results of actual samples were confirmed by existing standard determination methods. Result： The determination results of simulated samples were consistent with design of simulated samples. The animal-derived ingredients inconsistent with sample name or label identification were detected in actual samples， and the results of non-targeted screening method and existing standard methods were consistent. Conclusion： This method could quickly test the unknown animal-derived ingredients of samples， and the determination results were accurate and reliable. It could provide more powerful technical support to combat meat adulteration for the government.

Keywords： metagenomic sequencing technology; animal-derived ingredients; non-targeted; adulteration

肉类食品是人们餐桌必备的主要食品之一，常见的高经济价值的肉类食品有牛肉、羊肉及其制品等。但受经济利益驱动，肉类食品中掺杂使假、以次充好事件频繁发生。例如，欧洲市场的“马肉风波”，以鸭肉代替牛羊肉的“假羊肉串”“假牛肉干”“假肥牛片”等[1]。这种掺假行为不仅侵害了消费者的权益，扰乱市场，更降低了消费者对相关部门肉类食品监管的信任度。

我国现行有效的肉类食品动物源性成分检测标准仅能进行单一成分的定向检测。但随着监管力度的加大，肉类食品掺假物种的种类向多元化方向发展，样品掺假的物种及物种数目存在很大的不确定性。单一成分的定性检测极易造成漏检导致假阴性，给肉类食品市场监管带来巨大挑战。因此，急需建立肉类食品中动物源性成分非靶向筛查方法。

目前，可能实现动物源性成分非靶向筛查的技术分为两类：一类是基于质谱、光谱等产生的“指纹”图谱进行鉴别[2-3]；一类是基于分子生物学的方法。“指纹”图谱方法对基质要求高，仅适用于未加工的肉类，同时该方法需要复杂的数学计算。在分子生物学方法中，多重PCR法和微流控芯片法可以实现对样品多个靶标的同时检测，但若样品中含有实验用引物探针检测靶标之外的成分，仍会存在漏检的情况[4-7]。将DNA条形码技术与Sanger测序技术相结合，挑取多个细菌克隆的方式具有很强的随机性，该方法漏检的可能性仍然很大[8]。而基于二代测序的宏基因组技术在鉴定土壤、粪便等含有混合微生物的样品中菌种方面已得到广泛应用，该技术可以实现对扩增产物中多种扩增条带的同时检测，为基于宏基因组测序技术建立肉类食品中动物源性成分非靶向筛查方法奠定基础[9]。

本研究基于宏基因组测序技术建立肉类食品中动物源性成分非靶向筛查方法，并将该方法用于模拟样品和实际样品检测，以期为相关部门打击肉类食品掺假提供可靠的技术支撑。

1 材料与方法

1.1 材料、试剂与仪器

1.1.1 实验样品

实验中用到的猪、水牛、黄牛、山羊、绵羊、鸡、鸭、鹅、马和驴等肉样于农贸市场购买；牦牛、鼠、雀、鸵鸟、貂、骆驼、乳鸽、猫、雉鸡、鸽子和貉肉源从中国检验检疫科学研究院获得；香辣味牛肉粒、猪血豆腐、牛排等样品于超市购买；餐馆牛肉和羊肉串在餐饮店获得。

1.1.2 仪器与试剂

T100 PCR仪（美国Bio-Rad公司），Biospec-nano核酸蛋白测定仪（日本岛津公司），1300系列A2级别生物安全柜（美国Thermofisher公司），Sorvall™ ST 8离心机（美国Thermofisher公司），Sub-cell GT核酸电泳仪（美国Bio-Rad公司），PowerPac Basic电泳仪电泳（美国Bio-Rad公司），GelDoc XR+凝胶成像系统（美国Bio-Rad公司），QuantStudio 6 Flex荧光定量PCR仪（美国Thermofisher公司），Qubit 4.0荧光计（美国Thermofisher公司），MiSeq高通量测序仪（美国illumina公司）。

DNA提取试剂，DNeasy mericon Food Kit（QIAGEN，货号69514）；DNA浓度检测试剂盒，QubitTM 1×dsDNA HS Assay Kit（Invitrogen，货号Q33231）；PCR扩增试剂，PrimeSTAR® HS DNA Polymerase（Takara，货号R010A）；基因文库构建试剂，NexteraR XT Index Kit 24 Indexes-96 samples（illumina，货号15055293）；高通量测序试剂，PhiX Control Kit v3（illumina，货号FC-110-3001）和MiSeq Reagent Micro Kit v2（300-cycles）（illumina，货号FC-103-1002）；荧光PCR试剂，Takara Premix Ex Taq （Probe qPCR）（Takara，货号RR390A）。

用于宏基因组测序的扩增子上游引物COI-F序列（5’-3’）为TCGTCGGCAGCGTCAGATGTGTATAAGAGACAGCTCAACYAATCAYAAAGATATYGGCAC，下游引物COF-R序列（5’-3’）为GTCTCGTGGGCTCGGAGATGTGTATAAGAGACAGGGYATNACTATRAAGAAAATTATTAC。该对引物由Illumina宏基因组测序引物和推荐性国家标准《动物制品中动物源性检测基因条码技术 Sanger测序法》（GB/T 35918—2018）中的微型基因条码引物合并而成。相关引物、探针合成及Sanger测序服务均委托生工生物（上海）股份有限公司完成。

1.2 实验方法

1.2.1 动物源性成分非靶向筛查方法的技术路线

动物线粒体DNA上编码细胞色素酶亚基I（COI）基因序列相对保守，不同种动物的COI基因序列存在差异。本方法使用的扩增子引物为通用引物，可以扩增出不同种动物的COI基因序列，技术路线如图1所示。使用该扩增子引物对样品中的所有动物源性成分进行扩增，利用宏基因组测序技术对扩增产物中不同种动物的COI基因序列进行分析检测，再将获得的序列在基因数据库（美国国家生物信息中心）中比对，从而实现样品中动物源性成分的非靶向检测。

1.2.2 样品DNA提取

按照商业化DNA提取试剂盒说明书的要求提取DNA。提取过程中设置提取对照，以灭菌去离子水代替样品进行DNA提取。

1.2.3 DNA浓度测定

取3.5 μL已提取的DNA样品，采用核酸蛋白分析仪测定DNA浓度。

1.2.4 PCR扩增及扩增产物电泳检测

采用25 μL PCR反应体系，其中10×PCR反应Buffer 2.5 μL、dNTP 2 μL、上下游引物各1 μL（引物浓度为10 μmol·L-1）及DNA模板25～100 ng，补充无菌双蒸水至总体积为25 μL。

PCR反应程序：94 ℃预变性5 min；94 ℃变性30 s，51 ℃退火30 s，72 ℃延伸30 s，40个循环；

72 ℃，10 min。

将5 μL扩增产物与适量上样缓冲液混合，用2%的琼脂糖胶进行电泳检测。

1.2.5 基因文库构建

采用COI-F和COI-R为引物，以样品DNA为模板，进行第1轮PCR。25 μL反应体系：模板

5 μL，5×Buffer 5 μL，dNTP 2 μL，上游引物1 μL，下游引物1 μL，DNA聚合酶（高保真）0.2 μL，超纯水10.8 μL。反应程序：95 ℃，3 min；95 ℃ 30 s，55 ℃ 30 s，72 ℃ 30 s，30个循环；72 ℃，5 min；4 ℃保存。采用磁珠法对第1轮扩增产物进行纯化，以灭菌超纯水重悬纯化后的DNA。

再以纯化后的DNA为模板，以Nextera○R XT Index Kit 24 Indexes-96 samples中的Index primer为引物，进行第2轮PCR（接头Index PCR）。25 μL反应体系：模板2 μL，5×Buffer 5 μL，dNTP 2 μL，Index 1引物2.5 μL，Index 2引物2.5 μL，DNA聚合酶（高保真）0.2 μL，超纯水12.8 μL。反应程序：95 ℃，3 min；95 ℃ 30 s，55 ℃ 30 s，72 ℃ 30 s，7个循环；72 ℃，5 min；4 ℃保存。采用磁珠法对第2轮扩增产物进行纯化，以灭菌超纯水重悬纯化后的DNA，即为构建好的基因文库。