发明内容
本发明的第一个目的在于提供一种特异性强、灵敏度高、扩增效果好的用于检测总霍乱弧菌与致病性霍乱弧菌的ddPCR引物和探针;
本发明的第二个目的在于提供一种操作简单、精确度高的绝对定量检测总霍乱弧菌与致病性霍乱弧菌的ddPCR试剂盒;
本发明的第三个目的在于提供一种绝对定量检测总霍乱弧菌与致病性霍乱弧菌的ddPCR方法。
为达到上述目的,本发明采用下述技术方案:
本发明通过选取单拷贝具有种属特异性且高度保守的N-乙酰葡糖胺(GlcNAc)结合蛋白A的编码基因gbpA(GenBank no.EU072441.1)及毒力协同调节菌毛编码基因tcpA(GenBank no.KP187623.1)作为靶序列,各设计出3对引物和探针组合,具体序列见表1,经过筛选最终获得特异性强、灵敏度高、扩增效果好的用于检测总霍乱弧菌与致病性霍乱弧菌的ddPCR引物和探针(gbpA-F2/R2/P2和tcpA-F1/R1/P1),其中,gbpA基因探针5’端标记VIC,tcpA基因探针5’端标记FAM,两个探针的3’端均标记BHQ。
表1ddPCR引物/探针序列
a)引物/探针位置分别参考GenBank no.EU072441.1和KP187623.1序列信息。
本发明绝对定量检测总霍乱弧菌与致病性霍乱弧菌的ddPCR试剂盒,该试剂盒包含上述ddPCR引物和探序列。进一步,试剂盒还包括完成ddPCR反应所需材料和试剂,如阳性对照DNA、ddPCR master mix、微滴生成油、微滴生成卡、Twin Tec Semi-Skirted 96孔板、铝箔热封膜。
本发明还提供了绝对定量检测总霍乱弧菌与致病性霍乱弧菌的ddPCR方法,包括:
(1)提取待测样品基因组的DNA;
(2)配置ddPCR反应体系:具体见表2,包括序列表中SEQ ID NO.1至SEQ ID NO.6所示的引物和探针;
表2ddPCR反应体系
(3)微滴生成:将ddPCR反应体系和微滴生成油加入到微滴生成卡中,置于微滴生成仪中生成微滴;
(4)扩增反应:将生成的微滴转移到96孔板中,封膜后置于PCR仪上进行扩增反应,具体条件见表3;
表3ddPCR反应条件
注:升降温速度应≤2.5℃/s
(5)结果判定
将扩增后的96孔板置于微滴读取仪(QX100,Bio-Rad,Pleasanton,CA)中,利用QuantaSoft软件进行结果读取和分析。含有扩增产物的阳性微滴与不含扩增产物的阴性微滴会呈现出荧光信号强度的差异,以阴性微滴簇荧光幅度的最高点为界来设定阈值线。按照泊松分布原理计算获得待测样品中扩增靶基因gbpA和tcpA的拷贝数,分别代表总霍乱弧菌和致病性霍乱弧菌的含量。
本发明的有益效果如下:
本发明首次开发出能同时区分致病和非致病性霍乱弧菌的双重ddPCR定量检测方法和试剂盒,其不仅适用于粪便等临床样本,更适于对水产等食品样品的检测,为水产等高风险食品提供了一种新的霍乱弧菌定量检测方法,为保障食品安全、开展风险评估等工作提供了有力的武器。具体而言,该ddPCR定量检测技术较传统培养法和qPCR法具有以下优势:1)ddPCR针对霍乱弧菌单拷贝种属特异性基因gbpA和致病性基因tcpA的双重检测体系,通过单管反应即可鉴定出样品中是否含有产毒或非产毒霍乱弧菌,并对其实现定量检测;2)对gbpA和tcpA基因的最低定量限均可达100拷贝(CFU)/g,满足国内外对食品中霍乱弧菌污染最低控制水平的要求;3)无需进行增菌培养及后续的特征性菌落分离、生化分析、血清学鉴定等步骤,大幅缩短检测时间,简化操作流程;4)直接提取水产等样品匀浆DNA,无需对样品进行复杂的稀释和检测重复,减少了稀释偏差对定值结果的影响;5)ddPCR针对细菌基因组的检测,能有效检测到活的不可培养状态(VBNC)的细菌,从而弥补培养法在VBNC细菌检测中的欠缺;6)对靶基因的直接绝对定量,无需依赖任何核酸标准品,不需要绘制标准曲线,避免了由PCR扩增效率差异所导致的偏差;7)由于ddPCR是终点检测,对样品中抑制剂的耐受度更高;8)特别适于低污染水平样品的定量分析,具有更高的精确度和重复性,易于标准化。
具体实施方式
为了更清楚地说明本发明,下面结合优选实施例和附图对本发明做进一步的说明。附图中相似的部件以相同的附图标记进行表示。本领域技术人员应当理解,下面所具体描述的内容是说明性的而非限制性的,不应以此限制本发明的保护范围。
实施例1试验菌株及培养条件
共使用178株菌(表4),其中包括1株O1群霍乱弧菌、1株O139群霍乱弧菌、100株非O1/非O139群霍乱弧菌自分离菌株、其他弧菌标准/自分离菌株56株,以及其他常见的食源性致病菌标准菌株20株。所有弧菌和非弧菌菌株分别在T1N1琼脂平板和BHI琼脂平板上活化3代后,挑取单菌落分别接种于5mL碱性蛋白胨水(APW)和5mL脑心浸液肉汤(BHI)液体培养基中,37℃110r/min震荡培养18h。
表4实验用菌株列表
注:ATCC:美国典型培养物保藏中心;CMCC:医学微生物菌种保藏管理中心;CGMCC:中国一般微生物典藏中心;BJCIQ:北京出入境检验检疫局;GDCIQ:广东出入境检验检疫局。
实施例2ddPCR引物、探针的设计
为实现对总霍乱弧菌与致病性霍乱弧菌的特异性检测和绝对定量分析,我们分别选取单拷贝具有种属特异性且高度保守的N-乙酰葡糖胺(GlcNAc)结合蛋白A的编码基因gbpA(GenBank no.EU072441.1)及毒力协同调节菌毛编码基因tcpA(GenBankno.KP187623.1)作为靶序列,通过NCBI在线工具进行序列分析和比对,利用Prime Express软件V3.0(ABI,Foster City,CA,USA)各设计出3对引物和探针,序列见表1。
实施例3ddPCR检测
1、反应预混液的配制
按照表2配制20μL的ddPCR反应预混液。
2、微滴生成
将20μL的ddPCR反应预混液和70μL微滴生成油分别加入到8孔微滴生成卡中,置于微滴生成仪(QX100,Bio-Rad,Pleasanton,CA)中生成微滴。
3、扩增反应
将生成的油包水的微滴(40μL)缓慢转移至96孔板中,封膜后置于PCR仪(GeneAmp9700,Applied BioSystems,Foster City,CA)上进行扩增反应,扩增条件如表3所示。
4、结果判定
将扩增后的96孔板置于微滴读取仪(QX100,Bio-Rad,Pleasanton,CA)中,利用QuantaSoft软件进行结果读取和分析。含有扩增产物的阳性微滴与不含扩增产物的阴性微滴会呈现出荧光信号强度的差异,以阴性微滴簇荧光幅度的最高点为界来设定阈值线。按照泊松分布原理计算获得待测样品中扩增靶基因gbpA和tcpA的拷贝数,分别代表总霍乱弧菌和致病性霍乱弧菌的含量。
实施例4qPCR检测
为便于同ddPCR方法进行比较,选用同一套引物和探针进行qPCR检测。25μL的双重qPCR反应体系包括:2×Master Mix 12.5μL,10μM的gbpA-F、gbpA-R、tcpA-F和tcpA-R各1μL,10μM的gbpA-P和tcpA-P各0.5μL,待测DNA模板2μL。于ABI7900荧光PCR仪(美国,ABI公司)或Light cycler 480II(瑞士,罗氏公司)上进行扩增,50℃2min,95℃10min,进而95℃15s、60℃1min,45个循环。反应完成后利用SDS3.2或SW1.5.1软件进行分析,根据标准曲线计算靶基因的拷贝数。
实施例5ddPCR引物、探针的筛选
首先,采用qPCR方法对已设计合成的gbpA与tcpA基因引物探针进行筛选。通过对O139群霍乱弧菌基因组DNA的检测,发现所有引物探针均能有效扩增出靶基因,但引物探针组gbpA-F2/R2/P2和tcpA-F1/R1/P1的荧光信号出现时间最早,扩增效果最好(图1A、B)。为保证扩增反应的特异性,本实验选取gbpA-F2/R2/P2和tcpA-F1/R1/P1两组引物探针,对表4中所列常见弧菌和食源性致病菌DNA进行检测,结果发现O1群霍乱弧菌、O139群霍乱弧菌和非O1/非O139群霍乱弧菌均出现gbpA基因阳性扩增信号,而只有O1群霍乱弧菌和O139群霍乱弧菌出现tcpA基因阳性扩增信号,非O1/非O139群霍乱弧菌呈tcpA基因阴性,此外,其他常见弧菌及食源性致病菌gbpA和tcpA的检测均呈阴性(图1C、D)。表明这两组引物、探针均具有良好的特异性。
继而采用ddPCR对O139群霍乱弧菌基因组DNA进行检测,进一步验证上述引物、探针的有效性。结果发现所有引物、探针扩增后均能出现明显的阳性微滴簇和阴性微滴簇,且引物探针组gbpA-F2/R2/P2和tcpA-F1/R1/P1的荧光信号强,阴性和阳性微滴簇更加集中(图2A、B),扩增效果最好。进一步利用常见弧菌对ddPCR检测体系中引物探针的特异性进行验证,发现同qPCR检测结果一样,O1群、O139群和非O1/非O139群霍乱弧菌均出现gbpA基因阳性微滴簇,副溶血性弧菌、溶藻弧菌、创伤弧菌、拟态弧菌与河流弧菌gbpA基因ddPCR检测均呈阴性(图2C);此外,只有O1群和O139群霍乱弧菌出现tcpA基因阳性微滴簇,非O1/非O139群霍乱弧菌和其他5株弧菌均无tcpA基因阳性微滴出现(图2D)。综上可见,gbpA-F2/R2/P2和tcpA-F1/R1/P1灵敏、特异,且同时满足qPCR和ddPCR的扩增条件,将其用于后续的定量分析试验和试剂盒的组装。gbpA基因探针5’端标记VIC,tcpA基因探针5’端标记FAM,两个探针的3’端均标记BHQ,从而可在同一ddPCR反应体系中完成对两种基因的同步扩增检测,实现对总霍乱弧菌和致病性霍乱弧菌的分别鉴定。
实施例6检测总霍乱弧菌和致病性霍乱弧菌的ddPCR试剂盒
试剂盒主要包括:
(1)用于检测总霍乱弧菌的引物和探针(gbpA-2F、gbpA-2R、gbpA-2P);
(2)用于检测致病性霍乱弧菌的引物和探针(tcpA-1F、tcpA-1R、tcpA-1P);
(3)阳性对照:O139群霍乱弧菌基因组DNA,1000拷贝/μL;
(4)ddPCR master mix:购自美国BioRad公司;
(5)微滴生成油:购自美国BioRad公司;
(6)微滴生成卡:购自美国BioRad公司;
(7)Twin Tec Semi-Skirted 96孔板:购自德国Eppendorf公司;
(8)铝箔热封膜:购自美国BioRad公司。
实施例7霍乱弧菌纯菌液的检测
1、基因组DNA提取与核酸蛋白分析仪测定
利用氯化钠碱性蛋白胨水对O139群霍乱弧菌参照菌株的过夜培养物进行10倍梯度稀释。采用细菌基因组DNA提取试剂盒(天根生化科技有限公司,货号:DP302)提取各稀释度菌液DNA,提取步骤按照试剂盒说明书进行,最终将核酸溶解在50μL TE缓冲液中。对提取的核酸用核酸蛋白分析仪(Beckman Coulter,DU800)进行浓度和纯度的测定,确保提取核酸的A260/A280在1.8-2.0之间,并按下列公式计算提取霍乱弧菌基因组DNA的拷贝浓度:
其中m代表核酸蛋白分析仪测得的核酸浓度(ng/μL);n代表细菌基因组的长度(bp),根据NCBI上已经发布的霍乱弧菌基因组的测序数据,其平均长度为4.01×106bp。
O139群霍乱弧菌参照菌株在氯化钠碱性蛋白胨水中的过夜培养物1mL,提取的基因组DNA(共50μL)浓度为250ng/μL,代入上述公式算得提取到的基因组DNA浓度为5.7×107拷贝/μL,按照1拷贝基因组对应1个CFU菌落,则原始菌液浓度为2.85×109CFU/mL,其对数值为9.45。
2、平板计数
同时对10倍梯度稀释的O139群霍乱弧菌参照菌株纯菌液,分别吸取1mL菌悬液于无菌平皿内,每个稀释度做两个平皿。及时将15~20mL冷却至46℃的平板计数T1N1琼脂培养基倾注平皿,并转动平皿使其混合均匀。待琼脂凝固后,36℃倒置培养24h±2h。选取菌落数在30~300CFU之间、无蔓延菌落生长的平板计数菌落总数,每个稀释度菌落数应取两个平板的平均数。结果测得原始菌液浓度为1.1×109CFU/mL,其对数值为9.04。
3、qPCR和ddPCR检测
各稀释度菌液提取到的DNA同时进行双重qPCR和双重ddPCR检测,结果发现两者对gbpA与tcpA两个基因均可检测到10-7稀释度(图3,图4),具有相同的检测灵敏度。双重ddPCR在10-2稀释度达到检测上限,阳性微滴的比率达100%,自10-3至10-7稀释度菌液测得的gbpA与tcpA基因浓度均呈现出良好的线性关系(图4,表5)。在定值范围内ddPCR测定结果具有高度一致性,由不同稀释度菌液计算出原菌液中霍乱弧菌gbpA与tcpA基因含量的对数平均值值分别为9.22与9.25,同核酸蛋白分析仪和平板计数法测得值仅相差0.18~0.23个LOG值(图5)。
表5梯度稀释的O139群霍乱弧菌纯菌液ddPCR定量检测结果
实施例8人工污染虾皮样品的检测
选择经传统方法验证无霍乱弧菌的虾皮样品作为添加基质。对O139群霍乱弧菌参照菌株的过夜培养菌悬液进行10倍梯度稀释并平板计数,同时取10-4、10-5、10-6、10-7稀释度菌液各25mL分别添加到盛有25g虾皮样品的均质袋中,各添加浓度分别设5个样品重复,另取1份25g虾皮样品加入APW代替添加菌液作为阴性对照,上述样品用拍击式均质器拍击2min,制成系列1:1含有不同浓度O139群霍乱弧菌的污染样品。取各样品匀液2mL分别利用试剂盒(天根生化科技有限公司,货号:DP302)法提取细菌基因组DNA,最终将核酸溶解在50μL TE缓冲液中,进行双重ddPCR和双重qPCR检测,比较两者的定值效果。
在进行双重qPCR时,以10倍梯度稀释的O139群霍乱弧菌基因组DNA(经核酸蛋白分析仪测定其原液浓度为4.66×107拷贝/μL)作为标准品,同人工污染样品中提取的DNA一起进行扩增检测。以标准品检测Ct值作为纵坐标,其对应核酸浓度的对数值为横坐标,分别绘制得到gbpA基因和tcpA基因的标准曲线(图6A、B)。同时qPCR检测结果显示4个不同污染水平的虾皮样品均呈阳性扩增,其扩增曲线如图6C、D所示,根据检测Ct值和得到的标准曲线计算污染样品中gbpA基因平均含量分别是368500、32500、5500、300拷贝/g(表6);tcpA基因的平均含量则分别是238500、14050、3100、250拷贝/g(表6)。
而双重ddPCR检测结果显示,随着污染水平的下降,gbpA基因和tcpA基因出现的阳性微滴数量均逐渐减少(图7),根据阳性微滴数量和泊松分布原理,由QuantaSoft软件计算直接测定的gbpA基因含量分别为99250、9530、830、98拷贝/g(表6);tcpA基因的平均含量则分别是105600、9800、1020、118拷贝/g(表6)。
根据添加菌液的平板计数结果,预计10-7至10-4菌液污染样品中O139群霍乱弧菌含量为110~110000CFU/g。因此,按照1个拷贝靶基因对应1个CFU菌来计算,无论是gbpA还是tcpA基因,ddPCR的定值结果更接近平板计数预计值,线性分析发现ddPCR定值结果与平板计数结果存在0.90~0.96倍的关系(图8B);而qPCR的定值结果却明显偏高,约为平板计数预计值的2.16~3.35倍(图8A)。其次,ddPCR针对不同污染水平虾皮样品中gbpA与tcpA两个基因测定结果的拷贝数比值在1.03~1.20范围内(图9),比较接近理论比值1;而qPCR测得的gbpA与tcpA基因含量相差较大,其拷贝数比值在1.20~2.31范围内(图9),这与两个基因的qPCR扩增效率存在的差异密切相关(根据qPCR检测标准曲线的斜率计算得到gbpA与tcpA基因扩增效率分别是87.3%与95.3%)。而ddPCR是终点检测,定值结果不会受到扩增效率差异的影响,因此根据不同基因测得的目标菌含量具有更好的一致性,定值结果更加准确。此外,对于低污染水平样品(10-7菌液添加),qPCR的检测Ct值多在37~39之间,有1份平行样品的gbpA基因检测甚至为阴性,而ddPCR对于低污染水平样品的测定结果更加准确,在不经增菌培养的条件下,其对gbpA与tcpA两个基因的最低定量限可分别达98拷贝/g和118拷贝/g(表6),满足国内外对食品中致病性弧菌污染最低控制水平的要求。
表6ddPCR与qPCR对人工污染虾皮样品中O139群霍乱弧菌的定量检测结果
显然,本发明的上述实施例仅仅是为清楚地说明本发明所作的举例,而并非是对本发明的实施方式的限定,对于所属领域的普通技术人员来说,在上述说明的基础上还可以做出其它不同形式的变化或变动,这里无法对所有的实施方式予以穷举,凡是属于本发明的技术方案所引伸出的显而易见的变化或变动仍处于本发明的保护范围之列。
序列表
<110> 北京出入境检验检疫局检验检疫技术中心
<120> 一种绝对定量检测总霍乱弧菌与致病性霍乱弧菌的试剂盒及检测方法
<130> JLC16I0608E
<160> 18
<170> PatentIn version 3.3
<210> 1
<211> 21
<212> DNA
<213> 人工合成引物gbpA-F2
<400> 1
caaaccaaac tggaacccaa a 21
<210> 2
<211> 21
<212> DNA
<213> 人工合成引物gbpA-R2
<400> 2
tcaacgacac agaacggatt g 21
<210> 3
<211> 25
<212> DNA
<213> 人工合成引物gbpA-P2
<400> 3
ccattgtcgc gtgatgcatt tgacc 25
<210> 4
<211> 25
<212> DNA
<213> 人工合成引物tcpA-F1
<400> 4
gggatatgtt tccatttatc aacgt 25
<210> 5
<211> 21
<212> DNA
<213> 人工合成引物tcpA-R1
<400> 5
gcgacactcg tttcgaaatc a 21
<210> 6
<211> 25
<212> DNA
<213> 人工合成引物tcpA-P1
<400> 6
tgctttcgct gctgtcgctg atctt 25
<210> 7
<211> 24
<212> DNA
<213> 人工合成引物gbpA-F1
<400> 7
ggacgtgtca ccttgtgtaa attt 24
<210> 8
<211> 24
<212> DNA
<213> 人工合成引物gbpA-R1
<400> 8
cactttgtgg ttcgtattga atcg 24
<210> 9
<211> 28
<212> DNA
<213> 人工合成引物gbpA-P1
<400> 9
ccgctaacgg cactggagar aaaaacac 28
<210> 10
<211> 17
<212> DNA
<213> 人工合成引物gbpA-F3
<400> 10
cgtggctggc accaaag 17
<210> 11
<211> 21
<212> DNA
<213> 人工合成引物gbpA-R3
<400> 11
ccagagtatg gccatggctt a 21
<210> 12
<211> 27
<212> DNA
<213> 人工合成引物gbpA-P3
<400> 12
ctcgctagtg atggcgcaat ctaccaa 27
<210> 13
<211> 23
<212> DNA
<213> 人工合成引物tcpA-F2
<400> 13
tgatgaggca aagaatcctt tca 23
<210> 14
<211> 21
<212> DNA
<213> 人工合成引物tcpA-R2
<400> 14
ccaccgactg taattgcgaa t 21
<210> 15
<211> 30
<212> DNA
<213> 人工合成引物tcpA-P2
<400> 15
attttctcat ttccacgaaa ctctgcagcg 30
<210> 16
<211> 25
<212> DNA
<213> 人工合成引物tcpA-F3
<400> 16
cgattgattc gcagaatatg actaa 25
<210> 17
<211> 22
<212> DNA
<213> 人工合成引物tcpA-R3
<400> 17
tgcggtagct ggataattac ca 22
<210> 18
<211> 26
<212> DNA
<213> 人工合成引物tcpA-P3
<400> 18
ctgcgcaaaa tctaaacagc gtgcaa 26