基于数字PCR 的不同品种鸭组织中线粒体与核DNA拷贝数差异研究

纪艺 徐晓丽 姜媛媛 汪小福 徐俊锋 李玥莹 陈笑芸

（1. 沈阳师范大学生命科学学院，沈阳 110034；2. 浙江省农业科学院农产品质量标准研究所 省部共建农产品质量安全国家重点实验室， 杭州 310021）

食品安全关乎人们的身体健康和生命安全，已经成为当前全社会关注的热点问题。近年来肉制品掺假屡禁不止，肉制品以次充好，引起人们对食品安全的高度关注［1］。肉及肉制品动物源性成分核酸检测是判断掺假肉的最常用技术之一［2］。与核DNA相比，线粒体DNA 具有其特殊的母系遗传方式和胞内多拷贝带来的高灵敏度检测的特征，因此，除被广泛用于饲料、生鲜及肉制品的来源追踪和物种成分鉴别等研究外，线粒体DNA 在国内外还被用于掺假肉类的定量检测［3-7］。

鸭肉是我国较为便宜的畜牧业肉类品种，也是最常见的被用于掺假的肉品。近年来，在国内外PCR 技术在肉类掺假以及肉类品种来源分析中以及得到广泛应用［8-9］。林霖等［10］构建了实时荧光PCR检测方法用于检测番鸭、家鸭DNA 成分，为识别掺假肉类提供了技术支持。但是，荧光定量PCR 需要依赖标准曲线，由于不同组织中线粒体的拷贝数不同［11-12］，使得基于线粒体DNA 的肉类掺假定量检测存在客观的技术局限性。微滴式数字PCR 可以通过对样品进行微滴化处理，对每个微滴逐个进行检测，从而实现绝对定量。

本研究通过数字PCR 技术测定不同品种鸭不同组织中线粒体DNA 和核DNA 的拷贝数，对线粒体和核DNA 的拷贝数比值进行分析，旨在探寻鸭肉制品掺假定量检测的最佳组织来源，并为检测机构定量鸭肉掺假比例筛选合适的靶DNA 提供一定的 借鉴。

1 材料与方法

1.1 材料

不同品种（绍兴鸭、北京鸭、缙云麻鸭、攸县麻鸭、高邮麻鸭、金定麻鸭、山麻鸭、GD-绍兴鸭和马太湖鸭）的生鲜鸭来源于浙江省农业科学院畜牧兽医研究所和中国农业科学院北京畜牧兽医研究所，其中北京鸭和绍兴鸭的取样组织包括鸭血、鸭胸、鸭肝、鸭皮、鸭心及鸭腿；缙云麻鸭、攸县麻鸭、高邮麻鸭和GD-绍兴鸭的取样组织包括鸭胸和鸭腿；山麻鸭、金定麻鸭、攸县鸭和马太湖鸭的取样组织为鸭血。

动物组织DNA 试剂盒购自杭州新景生物试剂开发有限公司；血液DNA 提取试剂盒购自浙江易思得生物科技有限公司；预混液ddPCR Supermix for Probes 和Droplet Generation Oil for Probes 购 自 美 国Bio-Rad；引物和探针由杭州尚亚赛生物科技有限公司合成。

1.2 方法

1.2.1 样品处理和DNA 提取 用无菌手术刀将试验材料切成小块，先在液氮中将样品迅速冷冻并研磨成粉末。按照试剂盒说明书的要求称取一定量的样品，采用相应的试剂盒对所有样品进行 DNA 提取，每个样品设置3 个平行重复。

1.2.2 DNA 浓度、纯度和盐残留测定 各取1 μL 的DNA 样品溶液用Qubit 2.0 进行DNA 的浓度检测，利用NanoDrop 2000 进行纯度和盐残留等指标的 检测。

1.2.3 鸭源性微滴式数字PCR 检测体系的建立 鸭线粒体基因和核单拷贝基因的PCR 引物和探针序列参考文献合成［13］，线粒体12S rRNA 基因的上游引物序列为5′-GATCAAAATGCAACTAAGCTGTCGC-3′；下游引物序列为5′-GACCTGTCTTATTAGCGGGGTGCTG-3′；探针序列为FAM-CCCTAAATCTTGATACTTACCCTACCGAA-BHQ1；IL-2 核基因的上游序列为5′-GGAGCACCTCTATCAGAGAAAGACA-3′，下 游 序列 为5′-GTGTGTAGAGCTCAAGATCAATCCC-3′，探针序列为FAM-TGGGAACAAGCATGAATGTAAGTGGATGGT-BHQ1。将上述提取的不同鸭组织DNA 稀释至同一浓度后作为模板，进行微滴式数字PCR 反应。PCR 体系为20 μL，含前述步骤所得的DNA 模板2 μL，2×ddPCR Supermix for Probes 反应液10 μL，10 μmol/L 的上下游引物各1 μL，10 μmol/L 的探针0.5 μL，其余用ddH2O 补足到20 μL。

PCR 反应程序中退火温度的高低对微滴下雨现象有显著的影响，在反应程序的退火环节，设置了温度梯度。根据PCR 反应的热图，鸭源性基因组DNA 在不同的退火温度下，微滴的荧光强度存在较显著差异（图1）。比较不同退火温度的数字PCR 扩增热图，发现在退火温度为60℃时，没有明显的下雨现象，确定数字PCR 的最佳退火温度60℃。最终，确定微滴式数字PCR 反应程序如下：95℃预变性10 min；40 个循环（94℃ 变性15 s，60℃退火延伸1 min）；98℃变性10 min；4℃保存。

图1 微滴式数字PCR 退火温度优化扩增热图

1.2.4 DNA 拷贝数定量分析 微滴式数字PCR 结束后，将96 孔板置入微滴读取仪中读取信号，使用软件QuantaSoft Version 1.6.6.0320 分析实验数据，获得DNA 模板拷贝数的绝对定量结果。

2 结果

2.1 DNA质量测定

采用荧光计Qubit 2.0 测定DNA 浓度，利用NanoDrop 2000 核酸蛋白测定仪来测定DNA 纯度和盐残留情况，提取DNA 的浓度及纯度等参数见表1。结果显示，提取的DNA 浓度在20 ng/μL-200 ng/μL 之间，说明在不同组织中都能提取到较理想的DNA；OD260/OD280比值基本都在1.8-2.0 之间，说明蛋白质和RNA 等杂质残留少；OD260/OD230比值普遍大于2.0，说明核酸中盐分等杂质残留少，提取的DNA 质量较好。鉴于所提取组织的DNA 浓度不同，因此将工作浓度统一调整为8.89 ng/μL，并用于后续实验。

2.2 鸭源性微滴式数字PCR检测体系的建立

微滴式数字PCR 可实现对基因拷贝数的绝对定量，如在反应体系中加入GD-绍兴鸭胸和马太湖麻鸭血的DNA，利用微滴式数字PCR 分别读取IL-2核基因和12S rRNA 线粒体基因的拷贝数。从图2 可知，灰色为阴性微滴，蓝色为阳性微滴，在FAM 荧光通道上均扩增出鸭源性的信号，且得到良好的扩增，满足分析要求。从实验结果来看，不同鸭品种间的DNA 其IL-2 核基因的拷贝数较为稳定，而12S rRNA 线粒体基因的拷贝数变化较大，表明线粒体的丰度存在较大差异。

2.3 DNA拷贝数分析

结果如表2 所示，微滴式数字PCR 反应测得8.89 ng/μL 的DNA 中含有的核DNA 拷贝数平均值为13 470，与理论值［14］接近，所有组织中核DNA 拷贝数的变异系数为11.16%，小于线粒体DNA 拷贝数的变异系数（26.23%），表明核DNA 是开展鸭肉制品掺假定量检测的最适DNA 来源。线粒体DNA的拷贝数在不同鸭品种乃至同一品种不同组织间存在较大差异。整体而言，鸭肝中线粒体拷贝数最高；鸭心、鸭胸肉、鸭腿肉和鸭皮中的线粒体拷贝数次之；鸭血中的线粒体拷贝数最低。

表1 不同品种鸭组织DNA 的提取质量

2.4 线粒体和核DNA拷贝数稳定性分析

线粒体DNA 的拷贝数和核DNA 拷贝数的比值能反映提取DNA 中线粒体的相对丰度。如表2 所示，鸭肝中线粒体/核DNA 拷贝数比值最高；鸭心及鸭胸/腿肉中的拷贝数比值次之；鸭皮的拷贝数比值再次之；鸭血的拷贝数比值最低。分析6 个不同组织中线粒体/核DNA 拷贝数比值后可知，鸭肝中线粒体/核DNA 拷贝数比值的变异系数最大；鸭血和鸭心次之；而鸭腿肉的变异系数最小（4.10%）。

图2 微滴式数字PCR 检测不同品种鸭组织DNA 核和线粒体DNA 的拷贝数

3 讨论

以DNA 为基础的PCR 技术被广泛用于食品中动物源性成分的鉴别，目前国内外学者常利用线粒体基因来开展肉类掺假的定性和定量检测［15-18］。众所周知，线粒体是产生ATP 的细胞器，在能量代谢旺盛的细胞中，线粒体的数量相对丰富［19］。然而，由于不同组织的能量代谢水平不同，其线粒体的丰度也不尽相同，这势必影响基于线粒体DNA 的肉类掺假定量检测的准确性。本研究提取了9 个鸭品种不同组织的总DNA，并测定了线粒体DNA 和核DNA 的拷贝数及其比例情况。

从检测结果可知，不同品种鸭源性核DNA 的拷贝数稳定性好，而线粒体DNA 拷贝数在不同组织中存在较大的差异，核DNA 是开展鸭肉掺假定量检测的最优选择。在检测的6 个组织中，鸭肝中线粒体/核DNA 拷贝数比值及其最大，这与肝脏是重要的代谢器官密切相关［20］；鸭胸肉、鸭心、鸭腿肉和鸭皮次之，这与心肌细胞参与心脏舒缩活动和骨骼肌参与机体运动密不可分［21］；鸭血中拷贝数比值最低，说明在开展基于线粒体DNA 的掺假比例定量检测时须充分考虑掺假物的组织来源，以免影响检测的准确性。本研究发现鸭腿肉的变异系数最小，这可能与腿部骨骼肌在鸭子日常运动中发挥动力作用紧密相关［22］。因此，在开展基于核DNA 的鸭肉掺假比例定量检测时，不同品种以及组织均可进行测定；线粒体DNA 作为靶基因的鸭肉掺假比例定量检测时，鸭腿肉来源的肉制品是最佳选择。本文为鸭肉制品甚至其他肉制品开展掺假定量检测选择合适的靶DNA 提供参考。

表2 不同品种鸭组织中线粒体和核DNA 拷贝数及拷贝数比值情况

4 结论

核DNA 的拷贝数在不同品种鸭组织中具有良好的稳定性，是开展鸭肉制品掺假定量检测的最适DNA 来源。在开展利用线粒体DNA 为参比的掺假定量检测时，需充分考虑不同组织中线粒体的丰度问题，以免影响检测的准确性。不同品种鸭腿肉中线粒体/核DNA 拷贝数比值的变异系数最小。

致谢：

感谢浙江省农业科学院畜牧兽医研究所和中国农业科学院北京畜牧兽医研究所提供实验材料。