一种检测多种细菌的方法、检测多种细菌的系统及其应用与流程

1.本发明属于生物技术领域，具体涉及一种检测多种细菌的方法、检测多种细菌的系统及其应用，所述方法是一种专门针对多种细菌的基因组序列进行核酸扩增的方法。


背景技术：

2.核酸扩增技术(nucleic acid amplificationtechnology，naat)是一类分子生物学技术的总称，这类技术通过引物、dna聚合酶和其它试剂在特定温度下进行反应，实现微量核酸的快速、特异性扩增，可广泛应用于疾病诊断、病原微生物检测、食品安全检测、动植物检疫以及涉及分子克隆的各项应用，如测序、基因克隆、基因操作、等位基因分析、突变检测等。其中聚合酶链式反应(polymerase chain reaction，pcr)是最早实现应用的核酸扩增技术，它模仿dna在体内复制的过程，通过一对特异的寡核苷酸引物与待扩增dna片段互补，再经过若干个循环的变性、退火、延伸过程，实现dna片段的指数级扩增。pcr技术由cetus公司的mullis于1985年发明，随后saiki等人将热稳定dna聚合酶引入pcr反应体系，解决了早期pcr技术因热变性循环过程中dna聚合酶钝化导致需要不断人工添加聚合酶的缺点，大大提高了核酸扩增效率并使该技术可以实现自动化，自此以后，以pcr为代表的核酸扩增技术得到广泛的推广和应用，也极大的推动了分子生物学的发展。
3.众所周知，现有核酸扩增技术通常会避开富含gc或at的区域，其原因在于，这些区域易通过自身互补配对形成发夹环二级结构，从而阻碍引物与模板的结合；即使引物与模板能够勉强结合，也容易影响dna聚合酶延模板链的延伸，从而导致dna聚合酶沿模板扩增时“卡顿”并干扰dna合成。因此，为了提高反应成功率，传统的核酸扩增技术对扩增区域的碱基组成具有选择偏好性，其引物设计往往集中在gc含量40％～60％的区域，以45％～55％为最适宜，不同扩增技术略有差别，如pcr的引物设计建议g+c含量为40％～60％(见分子克隆实验指南第四版，表7-1引物设计)，lamp引物设计建议g+c含量为40％～65％(见a guide to lamp primer designing(primerexplorer v3))。
4.然而，自然界中生物体的核酸碱基组成差别巨大，如疟原虫基因组at碱基含量约为82％，周慧琦等研究了2670株细菌和古细菌，发现其基因组gc含量变化范围从14％到75％(周慧琦,2014)。从检测角度讲，现有核酸扩增技术的碱基组成选择偏好性忽略了大量潜在的靶序列，放大了目标对象被成功检出的难度；从基因操作角度讲，该选择偏好性导致相当多的目标序列难以被顺利扩增，从而难以开展进一步分子操作。如能针对传统核酸扩增方法所回避的gc碱基含量、at碱基含量非均衡域设计适当的引物，避免形成二级结构，同时设置适宜的反应条件，顺利开展核酸扩增实验，将显著增加目标对象被检出的成功率，同时使得更多序列能够被扩增进而方便地实施下游分子操作。
5.针对gc碱基含量、at碱基含量非均衡区域开展核酸扩增反应，其关键步骤是引物设计环节。鉴于以上分析，领域内公开发布的种属特异基因(或序列)通常来自于gc含量40％～60％甚至45％～55％的区域，因而无法应用于针对gc、at碱基含量非均衡区域的引物设计。另外，此类序列中连续存在的gc-rich或at-rich区，其碱基组成多样性相比其它区
域偏弱，又为特异性引物的设计带来额外的困难。因此首先需要提出一套能够针对gc、at碱基含量非均衡序列进行高通量、自动化、高效率的引物设计算法流程。
6.除了扩增区域选择和引物设计外，扩增区域和引物的tm值计算也是扩增过程能否顺利完成的重要影响因素。扩增区域和引物的tm值通常按照最近邻二态模型进行计算，但是不同研究者及引物生产公司采用的具体计算公式各有不同，如分子克隆实验指南第四版中建议的引物tm计算公式为tm＝4
×
(g,c数)+2
×
(a,t数)，takara公司建议的引物tm计算公式为tm＝4
×
(g,c数)+2
×
(a,t数)+32-2
×
(总碱基数)；生工生物的引物tm计算公式为tm(0.05mna
+
)＝59.94+1
×
(gc的百分含量)-(675/引物序列长度)。对于gc、at碱基含量非均衡区域而言，应当如何设计模型进行引物tm值计算，成为确保实验成功的重要因素。
7.综上所述，领域内亟需发展一种针对传统核酸扩增技术难以涉及的区域的核酸扩增技术，能够高通量、自动化地设计特异性引物，计算其tm值，并设置适宜的变性温度、退火温度等反应条件，拓展核酸扩增技术的应用范围，提升扩增的特异性，满足核酸检测、分子遗传学研究等多方面的需求。


技术实现要素：

8.本发明提出了一种创新的核酸扩增方法，并应用于检测多种细菌。该方法以富含at的核酸序列为扩增目标，进行目标序列识别、引物设计、变性温度及退火温度计算、模板局部解链的核酸扩增等。首先，采用自动化的引物设计流程，充分利用公共数据资源中丰富的基因组序列信息，高通量地设计针对富含at的目标序列的特异性引物对。其次，拟合了符合富含at序列特点的tm值计算公式，可根据理论扩增产物序列及引物的序列长度、碱基组成等因素计算变性温度和退火温度；考虑到富含at的核酸序列变性温度显著低于gc含量40％～60％以及高gc的序列，因此将计算所得的理论扩增产物序列tm值作为最低变性温度，在该温度下，能够在不添加任何化学变性剂的前提下，对双链dna的富含at区域进行局部解链，进而启动核酸扩增反应。
9.本发明一方面在引物设计环节引入公共数据资源中的海量基因组数据，进行序列特异性排查；另一方面，在变性温度设置环节有别于传统核酸扩增方法在93～95℃变性温度下打开所有双链结构，而是通过降低变性温度仅打开富含at序列，从而在源头上极大地减少了非靶向区域发生非特异扩增的可能性。这两方面的设计策略显著降低了核酸扩增反应中常见的非特异性扩增的可能性，在此基础上形成了针对富含at序列的核酸扩增方法。该方法的引物设计步骤采用c语言和perl语言编程实现，具有高通量、自动化的特点，整个扩增方法特异性强、成功率高，有效地拓展了传统核酸扩增技术的应用范围，同时在引物设计通量、反应特异性等方面显著提升了反应性能。
10.本发明中，所述核酸扩增方法包括识别待扩增核酸序列中富含at碱基的序列作为目标序列；针对目标序列自动化、高通量地设计特异性引物；设计特定的tm值计算公式计算反应变性及退火温度，并筛选合适的引物组合；设置核酸扩增反应条件，在模板局部解链的条件下进行核酸扩增反应。
11.具体包括如下步骤：
12.(1)筛选待扩增核酸序列中富含at碱基的目标序列；
13.(2)针对所述目标序列自动化、高通量地设计兼具通用性和特异性的引物；
14.(3)通过特定公式计算反应变性温度及退火温度，并设置核酸扩增反应条件；即，基于gc的百分含量、引物序列长度、引物对的理论扩增产物序列长度计算反应变性温度及退火温度，并设置核酸扩增反应条件；
15.(4)在模板局部解链的条件下进行核酸扩增反应，得到扩增产物。
16.本发明中，所述检测多种细菌的方法，具体步骤为：
17.(1)筛选不同细菌基因组序列中富含at碱基的目标序列；
18.(2)针对单一细菌基因组序列设计兼具通用性和特异性的引物；
19.(3)基于gc的百分含量、引物序列长度、引物对的理论扩增产物序列长度计算反应变性温度及退火温度；
20.(4)根据序列扩增长度、反应变性温度及退火温度筛选合适的单一细菌引物引物对并组成检测多种细菌的引物组；
21.(5)设置核酸扩增反应条件，在模板局部解链的条件下进行核酸扩增反应，得到扩增产物。
22.本发明步骤(1)中，对于待扩增核酸序列，从第一个碱基开始滑动宽度为1000bp的窗口，步长为5～100bp。针对每次窗口所在位置包含的序列计算at碱基含量，保留序列at碱基含量大于60％的区域作为目标序列。优选地，保留序列at碱基含量在60-80％的区域作为目标序列。
23.本发明步骤(1)中，对于多种细菌基因组序列，从第一个碱基开始滑动宽度为1000bp的窗口，步长为5～100bp。针对每次窗口所在位置包含的序列计算at碱基含量，保留序列at碱基含量大于60％的区域作为目标序列。优选地，保留序列at碱基含量在60-80％的区域作为目标序列。
24.本发明步骤(2)中，所述引物的设计方法，包括：(2.1)针对一条目标序列进行单条引物设计，得到候选引物；(2.2)对所述候选引物进行理化性质判定，筛选符合条件的单条引物；(2.3)将步骤(2.2)筛选得到的单条引物组合成引物对；(2.4)对所述引物对进行通用性、特异性的判定；(2.5)输出符合条件的引物对，得到所述特异性引物。
25.本发明中，可以采用任何能实现高通量的引物设计的编程语言来实现，如可操作性强、速度较快的c、perl等编程语言实现。
26.步骤(2.1)中，针对目标序列设计候选引物，候选引物需满足以下条件：a)引物序列长度在20bp～36bp之间；b)at碱基含量在55％～80％，c)连续gc数≤5；同时记录所述候选引物匹配到目标序列上的位置信息、正负链信息。
27.其中，优选地，所述at碱基含量为60％～75％。
28.其中，所述“连续gc数”是指引物序列中连续碱基g或者连续碱基c的数量，如引物序列aagggggttccaggcatta(seq id no.15)中连续gc数为5/2/2。
29.步骤(2.2)中，对满足上述(2.1)条件的单条引物进行理化性质判定，包含但不限于3’端稳定性、5’端稳定性和/或二级结构稳定性等物理化学性质，保留符合要求的单条引物。
30.其中，所述“符合要求的单条引物”是指具备3’端稳定性、5’端稳定性和二级结构稳定性的单条引物。
31.步骤(2.3)中，对上述步骤(2.2)筛选得到的单条引物，根据其匹配到目标序列上
的位置信息、正负链信息，将单条引物组合成引物对；引物对的理论扩增产物序列的长度应在200bp～600bp之间。根据公式0.466
×
(gc的百分含量)
×
100+66.04-(450/引物序列长度)计算单条引物的tm值；以引物对tm差值≤3℃，以及引物之间不能发生相互作用为条件，对引物对进行筛选，得到候选引物对。其中，gc的百分含量是指引物中的碱基g和c的个数占引物总碱基个数的百分比；引物序列长度是指引物的碱基个数。
32.步骤(2.4)中，对步骤(3)得到的候选引物对中的每一条引物做通用性、特异性判定。其中通用性判定指检查引物对是否严格匹配到所有的目标序列(如沙门氏菌多个菌株基因组)；
33.特异性判定指检查单条引物是否不能与所有目标序列之外的非目标序列发生特异性匹配，所述非目标序列指除待扩增核酸序列之外的其它核酸序列，所述特异性匹配指不大于2个错配的匹配。
34.通过通用性和特异性判定的候选引物可进入下一步骤。
35.步骤(2.5)中，输出符合条件的引物对，即得到所述核酸扩增反应的引物对。
36.本发明所述步骤(3)中，所述反应变性温度的计算公式为0.357
×
(gc的百分含量)
×
100+70.582-(990/引物对的理论扩增产物序列长度)，记为tma；其中，gc的百分含量是指引物对的理论扩增产物序列中的碱基g和c的个数占引物对的理论扩增产物序列总碱基个数的百分比。
37.本发明所述步骤(3)中，所述反应退火温度为上述步骤(2.3)中两条引物的tm值的平均值，记为tmb。即，根据公式0.466
×
(gc的百分含量)
×
100+66.04-(450/引物序列长度)计算单条引物的tm值，将所述tm值取平均值即为所述反应退火温度，记为tmb；其中，gc的百分含量是指引物中的碱基g和c的个数占引物总碱基个数的百分比。
38.本发明所述步骤(4)中，所述引物组合的筛选条件包括不同细菌扩增序列长度至少相差50bp，反应变性温度最多相差2℃，反应退火温度最多相差2℃。
39.本发明所述步骤(5)中，所述核酸扩增反应条件包括：在变性温度下反应5秒，在退火温度下反应5秒，在延伸温度下反应20秒，以上过程重复30-40次。其中变性温度根据需要可在所有引物对扩增片段平均变性温度
±
5℃之间调整；退火温度根据需要可在所有引物对平均退火温度
±
2℃之间调整；延伸温度为72℃。
40.本发明还提供了所述的方法在用于富含at区的目标核酸序列进行核酸扩增中的应用。
41.本发明还提供了如上所述设计方法获得的引物对。
42.本发明中，所述引物对为：
43.引物对a：
44.sal-f1：5
’‑
tcagacatccgttcagaaaat-3’(seq id no.1)
45.sal-r1：5
’‑
gttcaactgtcgacaagattaa-3’(seq id no.2)
46.引物对b：
47.sta-f1：5
’‑
gtaggtatggtaaatagttacac-3’(seq id no.3)
48.sta-r1:5
’‑
cactaatgccaaatttacttaaaatcg-3’(seq id no.4)
49.引物对c:
50.cro-f1：5
’‑
cgccataactgcataatcat-3’(seq id no.5)
51.cro-r1:5
’‑
ataacgagttaccgtgcaga-3’(seq id no.6)
52.引物对d:
53.yer-f1：5
’‑
tgtgcggtggatgtaaataattc-3’(seq id no.7)
54.yer-r1:5
’‑
gctttgaaactcaaggactg-3’(seq id no.8)
55.引物对e:
56.sta-f2：5
’‑
cctttcatctaaaaacctcca-3’(seq id no.9)
57.sta-r2:5
’‑
gaaatggatgttttaaaagaagg-3’(seq id no.10)
58.引物对f:
59.sal-f2：5
’‑
tgggttgaaatagcccatta-3’(seq id no.11)
60.sal-r2：5
’‑
gacgtgacacacttcgtttt-3’(seq id no.12)
61.引物对g：
62.yer-f2：5
’‑
atggaaaataacataatttctattaccgg-3’(seq id no.13)
63.yer-r2:5
’‑
tctctgcgaataaccttgtg-3’(seq id no.14)
64.本发明还提供了如上所述的方法或所述的引物对在多种细菌相关基因或区域的扩增和/或检测中的应用。
65.本发明还提供了如上所述的方法或所述的引物对在制备扩增和/或检测多种细菌基因组相关基因或区域的产品中的应用。
66.本发明提供了一种用于诊断多种细菌的试剂及其在检测多种细菌中的应用，所述诊断试剂包括如上所述的引物对。
67.本发明还提供了一种诊断试剂及其在检测多种细菌中的应用，所述诊断试剂包括如上所述的引物对。
68.所述多种细菌包括阪崎克罗诺杆菌、沙门氏菌、金黄色葡萄球菌、小肠结炎耶尔森菌等。
69.本发明还提供了一种同时检测多种细菌的系统，所述系统包括筛选模块，用于筛选不同细菌基因组序列中富含at碱基的目标序列；引物设计模块，其用于针对单一细菌基因组序列设计兼具通用性和特异性的引物；计算模块，用于基于gc的百分含量、引物序列长度、引物对的理论扩增产物序列长度计算反应变性温度及退火温度；组成模块，用于根据序列扩增长度、反应变性温度及退火温度筛选合适的单一细菌引物引物对并组成检测多种细菌的引物组；扩增模块，用于在模板局部解链的条件下进行核酸扩增反应，得到扩增产物。
70.本发明还提供了一种同时检测多种细菌的设备，所述设备包括：存储器、处理器和扩增模块；所述存储器上存储有计算机程序，当所述计算机程序被所述处理器执行时，实现上述所述方法的步骤1-4。
71.本发明的有益效果在于，本发明靶向传统核酸扩增方法通常需要规避的富含at碱基的序列，极大地拓展了候选靶标区域；通过基于海量基因组数据的特异性排查，以及针对富含at区域的tm值的精确计算，从而对双链dna进行局部解链，大大降低了核酸扩增技术中普遍存在的非特异性扩增的可能性。该方法有效地拓展了传统核酸扩增技术的应用范围，同时在引物设计通量、反应特异性等方面显著提升了反应性能，可满足核酸检测、分子遗传学研究等多方面的需求。
附图说明
72.图1为本发明针对肠炎沙门氏菌亚种、金黄色葡萄球菌亚种、小肠结肠炎耶尔森菌、阪崎克罗诺杆菌四种细菌核酸扩增反应的荧光染料染色结果。
73.图2为本发明针对肠炎沙门氏菌亚种、金黄色葡萄球菌亚种、小肠结肠炎耶尔森菌、阪崎克罗诺杆菌四种细菌核酸扩增反应的琼脂糖凝胶电泳结果。
74.图3为本发明针对肠炎沙门氏菌亚种、金黄色葡萄球菌亚种、小肠结肠炎耶尔森菌、阪崎克罗诺杆菌四种细菌的核酸扩增反应在81℃变性温度下的特异性检测的sybrgreen i染料显色结果。
75.图4为基于本发明的阪崎克罗诺杆菌(cronobactersakazakii)不同变性温度下核酸扩增特异性检测的染料显色结果。
76.图5为基于本发明的阪崎克罗诺杆菌(cronobactersakazakii)94℃变性温度下核酸扩增特异性检测的扩增产物电泳结果。
77.图6为基于本发明的阪崎克罗诺杆菌(cronobactersakazakii)82℃变性温度下核酸扩增特异性检测的扩增产物电泳结果。
78.图7为本发明针对肠炎沙门氏菌亚种、金黄色葡萄球菌亚种两种细菌核酸扩增反应的琼脂糖凝胶电泳结果。
79.图8为本发明针对肠炎沙门氏菌亚种、小肠结肠炎耶尔森菌两种细菌核酸扩增反应的琼脂糖凝胶电泳结果。
80.图9为本发明针对肠炎沙门氏菌亚种、金黄色葡萄球菌亚种和小肠结炎耶尔森菌三种细菌核酸扩增反应的琼脂糖凝胶电泳结果。
具体实施方式
81.结合以下具体实施例和附图，对本发明作进一步的详细说明。实施本发明的过程、条件、实验方法等，除以下专门提及的内容之外，均为本领域的普遍知识和公知常识，本发明没有特别限制内容。
82.本发明还提供了一种同时检测多种细菌的系统，所述系统包括筛选模块，用于筛选不同细菌基因组序列中富含at碱基的目标序列；引物设计模块，其用于针对单一细菌基因组序列设计兼具通用性和特异性的引物；计算模块，用于基于gc的百分含量、引物序列长度、引物对的理论扩增产物序列长度计算反应变性温度及退火温度；组成模块，用于根据序列扩增长度、反应变性温度及退火温度筛选合适的单一细菌引物引物对并组成检测多种细菌的引物组；扩增模块，用于在模板局部解链的条件下进行核酸扩增反应，得到扩增产物。
83.实施例1同时检测沙门氏菌、金黄色葡萄球菌、小肠结炎耶尔森菌和阪崎克罗诺杆菌
84.本发明针对沙门氏菌、金黄色葡萄球菌、小肠结炎耶尔森菌和阪崎克罗诺杆菌基因组筛选富含at的序列并设计特异性引物，根据计算的引物tm和目标序列tm筛选引物组合，设置反应条件进行核酸扩增，通过判断反应结果是否为阳性，确定待测样品中是否存在目标序列，进而确定待测样品中是否存在沙门氏菌、金黄色葡萄球菌、小肠结炎耶尔森菌和阪崎克罗诺杆菌中的一种或多种。具体步骤如下：
85.(1)富含at碱基的目标序列的筛选：
86.使用2019年8月5号从ncbi的ftp上下载的有完整基因组序列的细菌、古菌和病毒的数据，共2896个全基因组序列。设定集合a，其包含所有沙门氏菌基因组序列；设定集合b，其包含所有非沙门氏菌基因组序列。沙门氏菌基因组序列作为参考基因组，从基因组第一个碱基开始滑动宽度为1000bp的窗口，步长为50bp；每次滑动前计算一次序列at碱基含量，保留序列at碱基含量大于60％的区域作为候选目标序列。以上过程采用perl脚本实现。
87.(2)特异性引物的设计：
88.根据pcr引物的特性及本发明的需求，设定引物的at碱基含量为55％～75％、3’端稳定性δg《4、5’端稳定性δg《3及引物序列长度(20～36bp)等特性参数，同时设定单条引物不能产生发卡结构、自身不能发生相互作用等条件，以步骤(1)中的候选目标序列作为单条引物设计的候选序列，计算得出符合上述设定条件的单条引物。通过公式0.466
×
(gc的百分含量)
×
100+66.04-(450/引物序列长度)计算单条引物的tm值，并记录每条引物在目标序列上的位置、正负链信息(即来源于正链还是负链)以及引物的长度等信息。
89.根据单条引物的位置信息进行引物配对，保留同时满足引物对tm差值《3℃以及引物对扩增区域在200～600bp两个条件的引物对作为候选引物对。
90.使用比对软件bowtie，将上一步设计出的候选引物对中的每一条引物分别与集合a中的目标基因组序列和集合b中的非目标基因组进行序列比对。为了保证引物的通用性，当该单条引物与集合a中的目标序列进行比对时，使用参数设置
“‑aꢀ‑
n 0”，即要求该单条引物与目标序列完全匹配；为了保证引物的特异性，当该单条引物与集合b中的非目标序列进行比对时，使用参数设置
“‑aꢀ‑
n 3”，即允许该单条引物与非目标序列有不超过3个错配。系统输出满足条件的引物对，其中引物对的数量可预先设置，在本实施例中引物对的数量预先设置为100。以上高通量、自动化引物设计流程采用c、perl脚本实现。程序运行后在at富集区设计出沙门氏菌扩增引物对100个。
91.通过同样的方法分别设计出金黄色葡萄球菌、小肠结炎耶尔森菌和阪崎克罗诺杆菌扩增引物对100个。使用公式0.466
×
(gc的百分含量)
×
100+66.04-(450/引物序列长度)计算每条引物的tm并计算该引物对平均tm值作为反应退火温度；使用公式0.357
×
(gc的百分含量)
×
100+70.582-(990/扩增产物序列长度)计算每个引物对扩增区域的变性温度。
92.以不同细菌扩增序列长度至少相差50bp、反应变性温度最多相差2℃、反应退火温度最多相差2℃为条件筛选引物组合。系统输出满足条件的引物组合，数量可预先设置，在本实施例中引物组合的数量预先设置为10。随机选择一个引物组进行有效性验证。所述引物组合的序列为：
93.引物对a：
94.sal-f1：5
’‑
tcagacatccgttcagaaaat-3’(seq id no.1)
95.sal-r1：5
’‑
gttcaactgtcgacaagattaa-3’(seq id no.2)
96.引物对b：
97.sta-f1：5
’‑
gtaggtatggtaaatagttacac-3’(seq id no.3)
98.sta-r1:5
’‑
cactaatgccaaatttacttaaaatcg-3’(seq id no.4)引物对c:
99.cro-f1：5
’‑
cgccataactgcataatcat-3’(seq id no.5)
100.cro-r1:5
’‑
ataacgagttaccgtgcaga-3’(seq id no.6)
101.引物对d:
102.yer-f1：5
’‑
tgtgcggtggatgtaaataattc-3’(seq id no.7)
103.yer-r1:5
’‑
gctttgaaactcaaggactg-3’(seq id no.8)
104.所述引物对a的at碱基含量分别为60％和64％，理论平均退火温度为62℃；所述引物对的理论扩增产物序列片段长度为355bp，at的百分含量为71％，理论变性温度为78.14℃。
105.所述引物对b的at碱基含量分别为65.22％和70.37％，理论平均退火温度为62.5℃；所述引物对的理论扩增产物序列片段长度为465bp，at的百分含量为73.76％，理论变性温度为77.82℃。
106.所述引物对c的at碱基含量分别为60％和55％，理论平均退火温度为63℃；所述引物对的理论扩增产物序列片段长度为291bp，at的百分含量为68.73％，理论变性温度为78.34℃。
107.所述引物对d的at碱基含量分别为60.87％和55％，理论平均退火温度为64℃；所述引物对的理论扩增产物序列片段长度为235bp，at的百分含量为68.09％，理论变性温度为77.76℃。
108.(4)核酸扩增反应及结果检测：
109.核酸扩增反应体系配置如下表1所示。根据计算结果，采用81℃变性5秒，62℃退火5秒，72℃延伸20秒，推荐以上过程重复35次。检测对象清单见表2。
110.在使用普通梯度pcr仪反应完成后可通过两种方式判定扩增结果，一是加入终浓度25x的sybrgreen i染料，通过颜色判定扩增结果是否为阳性，即待测样品中是否存在目标序列；二是对扩增产物进行琼脂糖凝胶电泳，根据电泳条带判断扩增结果是否为阳性，即待测样品中是否存在目标序列。
111.表1.四种细菌核酸扩增反应体系
[0112][0113][0114]
表2.四种细菌核酸扩增反应特异性检测对象清单
[0115]
1.金黄色葡萄球菌2160016.福氏志贺氏菌1.18682.金黄色葡萄球菌亚种1.246517.大肠埃希氏菌107383.表皮葡萄球菌1.426018.致病性大肠埃希氏菌103724.马红球菌1.426219.致泻大肠埃希氏菌104115.蜡样芽胞杆菌1.376020.产肠毒素大肠埃希氏菌104156.蕈样芽胞杆菌2147321.肠产毒性大肠埃希氏菌106657.单增李斯特菌2163522.出血性大肠埃希氏菌215308.英诺克李斯特菌1041723.阪崎克罗诺杆菌215609.伊氏李斯特菌2166324.小肠结肠炎耶尔森菌2166910.肠炎沙门氏菌亚种1.185925.假结核耶尔森菌5350411.肠炎沙门氏菌2148226.创伤弧菌21615
12.鼠伤寒沙门氏菌1042027.副溶血弧菌1.199713.乙型副伤寒沙门菌1043728.弗式弧菌1.1613
ꢀꢀ
15.鲍氏志贺氏菌1.1061830.宋氏志贺氏菌
[0116]
图1和2展示了本发明针对肠炎沙门氏菌亚种、金黄色葡萄球菌亚种、小肠结肠炎耶尔森菌、阪崎克罗诺杆菌四种细菌核酸扩增反应的荧光染料染色结果和琼脂糖凝胶电泳结果。其中“n”为其对应的阴性模板扩增结果，“m1”指marker dl2000，“m2”指marker b(100～600bp)，“混合”指四种细菌基因组dna的混合模板。图1各单菌、混合模板及阴性模板扩增显色结果情况符合预期。图2各单菌扩增片段均符合理论计算值且混合模板种各菌扩增片段能彼此分开。
[0117]
图3展示了本发明针金黄色葡萄球菌、对沙门氏菌、阪崎克罗诺杆菌和小肠结炎耶尔森菌的核酸扩增反应在81℃变性温度下的特异性检测的sybrgreen i染料显色结果。若扩增产物呈亮绿色，则为阳性；若扩增产物呈橙色，则为阴性。图3中，1～2号金黄色葡萄球菌、金黄色葡萄球菌亚种，10～13号肠炎沙门氏菌亚种、肠炎沙门氏菌、鼠伤寒沙门氏菌和乙型副伤寒沙门氏菌，23号阪崎克罗诺杆菌，24号小肠结炎耶尔森菌呈现亮绿色，为阳性结果，符合预期；同时，其它细菌基因组dna模板检测结果均呈阴性，如第3～9、14～22号、25～30号管所示。所有结果符合预期，提示在81℃变性温度下未产生非特异性扩增或即使产生微量非特异扩增但不足以影响染料法结果判定。
[0118]
实施例2阪崎克罗诺杆菌(cronobactersakazakii)的检测
[0119]
本发明针对阪崎克罗诺杆菌(cronobactersakazakii)基因组筛选富含at的序列并设计特异性引物，根据计算的引物tm和目标序列tm设置反应条件进行核酸扩增，通过判断反应结果是否为阳性，确定待测样品中是否存在目标序列，进而确定待测样品中是否存在阪崎克罗诺杆菌。具体步骤如下：
[0120]
(1)富含at碱基的目标序列的筛选：
[0121]
使用2019年8月5号从ncbi的ftp上下载的有完整基因组序列的细菌、古菌和病毒的数据，共2896个全基因组序列。设定集合a，其包含所有阪崎克罗诺杆菌基因组序列；设定集合b，其包含所有非阪崎克罗诺杆菌基因组序列。以gi号为156932229的阪崎克罗诺杆菌基因组序列作为参考基因组，从基因组第一个碱基开始滑动宽度为1000bp的窗口，步长为50bp；每次滑动前计算一次序列at碱基含量，保留序列at碱基含量大于60％的区域作为候选目标序列。以上过程采用perl脚本实现。
[0122]
(2)特异性引物的设计：
[0123]
根据pcr引物的特性及本发明的需求，设定引物的at碱基含量为55％～75％、3’端稳定性δg《4、5’端稳定性δg《3及引物序列长度(20～36bp)等特性参数，同时设定单条引物不能产生发卡结构、自身不能发生相互作用等条件，以步骤(1)中的候选目标序列作为单条引物设计的候选序列，计算得出符合上述设定条件的单条引物。通过公式0.466
×
(gc的百分含量)
×
100+66.04-(450/引物序列长度)计算单条引物的tm值，并记录每条引物在目标序列上的位置、正负链信息(即来源于正链还是负链)以及引物的长度等信息。
[0124]
根据单条引物的位置信息进行引物配对，保留同时满足引物对tm差值《3℃以及引物对扩增区域在200～600bp两个条件的引物对作为候选引物对。
[0125]
使用比对软件bowtie，将上一步设计出的候选引物对中的每一条引物分别与集合a中的目标基因组序列和集合b中的非目标基因组进行序列比对。为了保证引物的通用性，当该单条引物与集合a中的目标序列进行比对时，使用参数设置
“‑
a-n 0”，即要求该单条引物与目标序列完全匹配；为了保证引物的特异性，当该单条引物与集合b中的非目标序列进行比对时，使用参数设置
“‑
a-n 3”，即允许该单条引物与非目标序列有不超过3个错配。系统输出满足条件的引物对，其中引物对的数量可预先设置，在本实施例中引物对的数量预先设置为20。以上高通量、自动化引物设计流程采用c、perl脚本实现。
[0126]
程序运行后在at富集区设计出pcr反应引物对20个，随机选择一个引物组进行有效性验证。所述引物对的序列为：
[0127]
cro-f1：5
’‑
cgccataactgcataatcat-3’(seq id no：5)
[0128]
cro-r1：5
’‑
ataacgagttaccgtgcaga-3’(seq id no：6)
[0129]
所述引物的at碱基含量分别为60％和55％，所述引物对的理论扩增产物序列at的百分含量为69％。
[0130]
(3)pcr反应变性温度和退火温度的计算：
[0131]
使用公式0.466
×
(gc的百分含量)
×
100+66.04-(450/引物序列长度)计算可知引物cro-f的tm＝62.2℃，cro-r的tm＝64.5℃，该引物对平均tm值为63.4℃；使用公式0.357
×
(gc的百分含量)
×
100+70.582-(990/扩增产物序列长度)计算扩增区域变性温度为78.3℃。
[0132]
(4)核酸扩增反应及结果检测：
[0133]
核酸扩增反应体系配置如下表3所示。根据计算结果，实验共设3组，分别测试在94℃/90℃/82℃变性温度下引物针对不同检测对象的特异性，核酸扩增反应条件为94℃/90℃/82℃变性5秒，63℃退火5秒，72℃延伸20秒，以上过程重复35次。特异性检测对象清单见表4。
[0134]
反应完成后通过两种方式判定扩增结果，一是加入终浓度25x的sybrgreen i染料，通过颜色判定扩增结果是否为阳性，即待测样品中是否存在目标序列；二是对扩增产物进行琼脂糖凝胶电泳，根据电泳条带判断扩增结果是否为阳性，即待测样品中是否存在目标序列。
[0135]
表3.阪崎克罗诺杆菌(cronobactersakazakii)核酸扩增反应体系
[0136]
[0137][0138]
表4.阪崎克罗诺杆菌(cronobactersakazakii)核酸扩增反应特异性检测对象清单
[0139]
1.金黄色葡萄球菌2160016.福氏志贺氏菌1.18682.金黄色葡萄球菌亚种1.246517.大肠埃希氏菌107383.表皮葡萄球菌1.426018.致病性大肠埃希氏菌103724.马红球菌1.426219.致泻大肠埃希氏菌104115.蜡样芽胞杆菌1.376020.产肠毒素大肠埃希氏菌104156.蕈样芽胞杆菌2147321.肠产毒性大肠埃希氏菌106657.单增李斯特菌2163522.出血性大肠埃希氏菌215308.英诺克李斯特菌1041723.阪崎克罗诺杆菌215609.伊氏李斯特菌2166324.小肠结肠炎耶尔森菌2166910.肠炎沙门氏菌亚种1.185925.假结核耶尔森菌5350411.肠炎沙门氏菌2148226.创伤弧菌2161512.鼠伤寒沙门氏菌1042027.副溶血弧菌1.199713.乙型副伤寒沙门菌1043728.弗式弧菌1.161314.痢疾志贺氏菌1.186929.霍乱弧菌1.867615.鲍氏志贺氏菌1.1061830.福式志贺氏菌
[0140]
实验结果见图4、5、6，其中1-22分别为金黄色葡萄球菌、金黄色葡萄球菌金黄亚种、表皮葡萄球菌、马红球菌、蜡样芽孢杆菌、蕈样芽孢杆菌、单核增生李斯特氏菌、英诺克李斯特氏菌、伊氏李斯特氏菌、肠炎沙门氏菌亚种、肠炎沙门氏菌、鼠伤寒沙门氏菌、乙型副伤寒沙门氏菌、痢疾志贺氏菌、鲍氏志贺氏菌、福氏志贺氏菌、大肠埃希氏菌(含肉毒梭菌a型基因)、致病性大肠埃希氏菌、致泻大肠埃希氏菌、产肠毒素大肠埃希氏菌、肠产毒性大肠埃希氏菌、出血性大肠埃希氏菌，24-30分别为小肠结肠炎耶尔森氏菌、假结核耶尔森氏菌、创伤弧菌、副溶血弧菌、弗氏弧菌、霍乱弧菌和福式志贺氏菌，n：阴性对照，p：阳性对照(含目的序列的质粒)；23为阪崎克罗诺杆菌。
[0141]
图4展示了本发明针对阪崎克罗诺杆菌(cronobactersakazakii)的核酸扩增反应在不同变性温度下的特异性检测的sybrgreen i染料显色结果。若扩增产物呈亮绿色，则为阳性；若扩增产物呈橙色，则为阴性。图4中，当变性温度为94℃/90℃，阴性对照(n)呈现橙色，为阴性结果，符合预期；阳性对照(p)和23号阪崎克罗诺杆菌呈现亮绿色，为阳性结果，符合预期。说明整个反应体系能正常工作。然而，其它细菌基因组dna模板扩增产物出现大量阳性结果，如变性温度94℃条件下的3、4、10～18、20～22、24、26、27、29、30号，变性温度
90℃条件下的3、7、10～18、20～22、24、26、27、29、30号管所示，与预期不符，提示在94℃/90℃变性温度下易产生非特异性扩增，导致结果呈假阳性。当变性温度为82℃，阴性对照(n)呈现橙色，为阴性结果，符合预期；阳性对照(p)和23号阪崎克罗诺杆菌呈现亮绿色，为阳性结果，符合预期。说明整个反应体系能正常工作。同时，其它细菌基因组dna模板检测结果均呈阴性，如第1～22、24～30号管所示。所有结果符合预期，提示在82℃变性温度下未产生非特异性扩增或即使产生微量非特异扩增但不足以影响染料法结果判定。
[0142]
图5展示了本发明针对阪崎克罗诺杆菌(cronobactersakazakii)的核酸扩增反应在94℃变性温度下的特异性检测的扩增产物电泳结果。若扩增产物电泳后在291bp处有一条单一的条带，则为阳性；若扩增产物电泳后无条带，则为阴性；若扩增反应产物电泳后在291bp之外出现一条或多条条带，则为非特异扩增所导致的假阳性。图5中，阴性对照(n)无条带，符合预期；阳性对照(p)和23号阪崎克罗诺杆菌在291bp处均有一条清晰条带，符合预期。说明整个反应体系能正常工作。然而，其它细菌基因组dna模板扩增产物在500bp～2000bp范围内出现多条条带，与预期不符，提示在94℃变性温度下反应体系产生大量非特异性扩增，该结果与图4中染料显色结果一致。
[0143]
图6展示了本发明针对阪崎克罗诺杆菌(cronobactersakazakii)的核酸扩增反应在82℃变性温度下的特异性检测的扩增产物电泳结果。若扩增产物电泳后在291bp处有一条单一的条带，则为阳性；若扩增反应产物电泳后无条带，则为阴性；若扩增反应产物电泳后在291bp之外出现一条或多条条带，则为非特异扩增所导致的假阳性。图6中，阴性对照(n)无条带，符合预期；阳性对照(p)和23号阪崎克罗诺杆菌在291bp处均有一条清晰条带，符合预期。说明整个反应体系能正常工作。同时，其它细菌基因组dna模板扩增反应产物大部分无条带，呈阴性结果，如2～6、10、12、13、15～22、26～30号泳道所示；少数模板扩增反应产物电泳后出现极微弱条带，易于与阳性结果区分，如1、7、8、9、11、14、24、25管所示，提示在82℃变性温度下不产生或产生在电泳结果中易于区分的微量非特异性扩增。
[0144]
综合图4～6及表4可看出，本发明核酸扩增方法针对阪崎克罗诺杆菌的应用具有良好的菌株特异性，即，在本发明提出的反应体系和反应条件下，仅阪崎克罗诺杆菌呈现阳性结果，而非阪崎克罗诺杆菌为阴性结果。值得注意的是，当使用传统pcr方法的变性温度(90℃或者94℃)进行反应时，尽管电泳结果(图2)显示目标检测对象(阪崎克罗诺杆菌，第23泳道)有单一清晰的条带，然而非目标检测对象普遍出现大量非特异性扩增，如第3泳道的表皮葡萄球菌；若使用染料显色法判定检测结果，如图4所示，在90℃和94℃变性条件下，呈现大量假阳性结果，从而无法区分检测对象和非检测对象。在同样的反应体系下，当使用本发明计算生成的变性温度82℃进行反应时，电泳结果(图6)显示目标检测对象(阪崎克罗诺杆菌，第23泳道)有单一清晰的条带，而非目标检测对象少有非特异性扩增发生，即使有，非特异性扩增条带也很微弱，极易与阳性结果区分；若使用染料显色法判定检测结果，如图4所示，在82℃变性条件下，仅目标检测对象(阪崎克罗诺杆菌，23号管)呈阳性结果，非目标检测对象均呈阴性结果，表现出良好的区分性。
[0145]
实施例3同时检测沙门氏菌和金黄色葡萄球菌
[0146]
本发明针对沙门氏菌、金黄色葡萄球菌基因组筛选富含at的序列并设计特异性引物，根据计算的引物tm和目标序列tm筛选引物组合，设置反应条件进行核酸扩增，通过判断反应结果是否为阳性，确定待测样品中是否存在目标序列，进而确定待测样品中是否存在
沙门氏菌、金黄色葡萄球菌中的一种或多种。具体步骤如下：
[0147]
(1)富含at碱基的目标序列的筛选：
[0148]
使用2019年8月5号从ncbi的ftp上下载的有完整基因组序列的细菌、古菌和病毒的数据，共2896个全基因组序列。设定集合a，其包含所有沙门氏菌基因组序列；设定集合b，其包含所有非沙门氏菌基因组序列。沙门氏菌基因组序列作为参考基因组，从基因组第一个碱基开始滑动宽度为1000bp的窗口，步长为50bp；每次滑动前计算一次序列at碱基含量，保留序列at碱基含量大于60％的区域作为候选目标序列。以上过程采用perl脚本实现。
[0149]
(2)特异性引物的设计：
[0150]
根据pcr引物的特性及本发明的需求，设定引物的at碱基含量为55％～75％、3’端稳定性δg《4、5’端稳定性δg《3及引物序列长度(20～36bp)等特性参数，同时设定单条引物不能产生发卡结构、自身不能发生相互作用等条件，以步骤(1)中的候选目标序列作为单条引物设计的候选序列，计算得出符合上述设定条件的单条引物。通过公式0.466
×
(gc的百分含量)
×
100+66.04-(450/引物序列长度)计算单条引物的tm值，并记录每条引物在目标序列上的位置、正负链信息(即来源于正链还是负链)以及引物的长度等信息。
[0151]
根据单条引物的位置信息进行引物配对，保留同时满足引物对tm差值《3℃以及引物对扩增区域在200～600bp两个条件的引物对作为候选引物对。
[0152]
使用比对软件bowtie，将上一步设计出的候选引物对中的每一条引物分别与集合a中的目标基因组序列和集合b中的非目标基因组进行序列比对。为了保证引物的通用性，当该单条引物与集合a中的目标序列进行比对时，使用参数设置
“‑
a-n 0”，即要求该单条引物与目标序列完全匹配；为了保证引物的特异性，当该单条引物与集合b中的非目标序列进行比对时，使用参数设置
“‑
a-n 3”，即允许该单条引物与非目标序列有不超过3个错配。系统输出满足条件的引物对，其中引物对的数量可预先设置，在本实施例中引物对的数量预先设置为100。以上高通量、自动化引物设计流程采用c、perl脚本实现。程序运行后在at富集区设计出沙门氏菌扩增引物对100个。
[0153]
通过同样的方法分别设计出金黄色葡萄球菌扩增引物对100个。使用公式0.466
×
(gc的百分含量)
×
100+66.04-(450/引物序列长度)计算每条引物的tm并计算该引物对平均tm值作为反应退火温度；使用公式0.357
×
(gc的百分含量)
×
100+70.582-(990/扩增产物序列长度)计算每个引物对扩增区域的变性温度。
[0154]
以不同细菌扩增序列长度至少相差50bp、反应变性温度最多相差2℃、反应退火温度最多相差2℃为条件筛选引物组合。系统输出满足条件的引物组合，数量可预先设置，在本实施例中引物组合的数量预先设置为10。随机选择一个引物组进行有效性验证。所述引物组合的序列为：
[0155]
引物对a：
[0156]
sal-f1：5
’‑
tcagacatccgttcagaaaat-3’(seq id no.1)
[0157]
sal-r1：5
’‑
gttcaactgtcgacaagattaa-3’(seq id no.2)
[0158]
引物对e：
[0159]
sta-f2：5
’‑
cctttcatctaaaaacctcca-3’(seq id no.9)
[0160]
sta-r2:5
’‑
gaaatggatgttttaaaagaagg-3’(seq id no.10)
[0161]
所述引物对a的at碱基含量分别为60％和64％，理论平均退火温度为62℃；所述引
物对的理论扩增产物序列片段长度为355bp，at的百分含量为71％，理论变性温度为78.14℃。
[0162]
所述引物对e的at碱基含量分别为61.90％和69.57％，理论平均退火温度为61℃；所述引物对的理论扩增产物序列片段长度为595bp，at的百分含量为76.30％，理论变性温度为77.38℃。
[0163]
(4)核酸扩增反应及结果检测：
[0164]
核酸扩增反应体系配置如下表5所示。根据计算结果，采用81℃变性5秒，62℃退火5秒，72℃延伸20秒，推荐以上过程重复35次。检测对象清单见表2。
[0165]
在使用普通梯度pcr仪反应完成后通过对扩增产物进行琼脂糖凝胶电泳，根据电泳条带判断扩增结果是否为阳性，即待测样品中是否存在目标序列。
[0166]
表5.两种细菌核酸扩增反应体系
[0167]
体系体积(μl)终浓度taqr300mix151xsal-f1/r1(50μm)0.25+0.250.5μmsta-f2/r2(50μm)0.25+0.250.5μm100％dmso0.52％模板(10ng)0/1/ddh2oupto25/
[0168]
图7展示了本发明针对肠炎沙门氏菌亚种和金黄色葡萄球菌亚种两种细菌核酸扩增反应的琼脂糖凝胶电泳结果。其中“n”为其对应的阴性模板扩增结果，“m”指marker dl2000，“沙门”指肠炎沙门氏菌亚种，“金葡”指金黄色葡萄球菌亚种，“混合”指两种种细菌基因组dna的混合模板。图7各单菌扩增片段均符合理论计算值且混合模板种各菌扩增片段能彼此分开。
[0169]
实施例4同时检测沙门氏菌和小肠结炎耶尔森菌
[0170]
本发明针对沙门氏菌、小肠结炎耶尔森菌基因组筛选富含at的序列并设计特异性引物，根据计算的引物tm和目标序列tm筛选引物组合，设置反应条件进行核酸扩增，通过判断反应结果是否为阳性，确定待测样品中是否存在目标序列，进而确定待测样品中是否存在沙门氏菌、小肠结炎耶尔森菌菌中的一种或多种。具体步骤如下：
[0171]
(1)富含at碱基的目标序列的筛选：
[0172]
使用2019年8月5号从ncbi的ftp上下载的有完整基因组序列的细菌、古菌和病毒的数据，共2896个全基因组序列。设定集合a，其包含所有沙门氏菌基因组序列；设定集合b，其包含所有非沙门氏菌基因组序列。沙门氏菌基因组序列作为参考基因组，从基因组第一个碱基开始滑动宽度为1000bp的窗口，步长为50bp；每次滑动前计算一次序列at碱基含量，保留序列at碱基含量大于60％的区域作为候选目标序列。以上过程采用perl脚本实现。
[0173]
(2)特异性引物的设计：
[0174]
根据pcr引物的特性及本发明的需求，设定引物的at碱基含量为55％～75％、3’端稳定性δg《4、5’端稳定性δg《3及引物序列长度(20～36bp)等特性参数，同时设定单条引物不能产生发卡结构、自身不能发生相互作用等条件，以步骤(1)中的候选目标序列作为单条引物设计的候选序列，计算得出符合上述设定条件的单条引物。通过公式0.466
×
(gc的
百分含量)
×
100+66.04-(450/引物序列长度)计算单条引物的tm值，并记录每条引物在目标序列上的位置、正负链信息(即来源于正链还是负链)以及引物的长度等信息。
[0175]
根据单条引物的位置信息进行引物配对，保留同时满足引物对tm差值《3℃以及引物对扩增区域在200～600bp两个条件的引物对作为候选引物对。
[0176]
使用比对软件bowtie，将上一步设计出的候选引物对中的每一条引物分别与集合a中的目标基因组序列和集合b中的非目标基因组进行序列比对。为了保证引物的通用性，当该单条引物与集合a中的目标序列进行比对时，使用参数设置
“‑
a-n 0”，即要求该单条引物与目标序列完全匹配；为了保证引物的特异性，当该单条引物与集合b中的非目标序列进行比对时，使用参数设置
“‑
a-n 3”，即允许该单条引物与非目标序列有不超过3个错配。系统输出满足条件的引物对，其中引物对的数量可预先设置，在本实施例中引物对的数量预先设置为100。以上高通量、自动化引物设计流程采用c、perl脚本实现。程序运行后在at富集区设计出沙门氏菌扩增引物对100个。
[0177]
通过同样的方法分别设计出小肠结炎耶尔森菌扩增引物对100个。使用公式0.466
×
(gc的百分含量)
×
100+66.04-(450/引物序列长度)计算每条引物的tm并计算该引物对平均tm值作为反应退火温度；使用公式0.357
×
(gc的百分含量)
×
100+70.582-(990/扩增产物序列长度)计算每个引物对扩增区域的变性温度。
[0178]
以不同细菌扩增序列长度至少相差50bp、反应变性温度最多相差2℃、反应退火温度最多相差2℃为条件筛选引物组合。系统输出满足条件的引物组合，数量可预先设置，在本实施例中引物组合的数量预先设置为10。随机选择一个引物组进行有效性验证。所述引物组合的序列为：
[0179]
引物对f：
[0180]
sal-f2：5
’‑
tgggttgaaatagcccatta-3’(seq id no.11)
[0181]
sal-r2：5
’‑
gacgtgacacacttcgtttt-3’(seq id no.12)
[0182]
引物对g：
[0183]
yer-f2：5
’‑
atggaaaataacataatttctattaccgg-3’(seq id no.13)
[0184]
yer-r2:5
’‑
tctctgcgaataaccttgtg-3’(seq id no.14)
[0185]
所述引物对f的at碱基含量分别为60％和55％，理论平均退火温度为63℃；所述引物对的理论扩增产物序列片段长度为352bp，at的百分含量为74.43％，理论变性温度为76.90℃。
[0186]
所述引物对g的at碱基含量分别为72.41％和55％，理论平均退火温度为63.5℃；所述引物对的理论扩增产物序列片段长度为242bp，at的百分含量为72.02％，理论变性温度为76.48℃。
[0187]
(4)核酸扩增反应及结果检测：
[0188]
核酸扩增反应体系配置如下表6所示。根据计算结果，采用81℃变性5秒，62℃退火5秒，72℃延伸20秒，推荐以上过程重复35次。检测对象清单见表2。
[0189]
在使用普通梯度pcr仪反应完成后通过对扩增产物进行琼脂糖凝胶电泳，根据电泳条带判断扩增结果是否为阳性，即待测样品中是否存在目标序列。
[0190]
表6.两种细菌核酸扩增反应体系
[0191]
体系体积(μl)终浓度
taqr300mix151xsal-f2/r2(50μm)0.25+0.250.5μmyer-f2/r2(50μm)0.25+0.250.5μm100％dmso0.52％模板(10ng)0/1/ddh2oupto25/
[0192]
图8展示了本发明针对肠炎沙门氏菌亚种和小肠结炎耶尔森菌两种细菌核酸扩增反应的琼脂糖凝胶电泳结果。其中“n”为其对应的阴性模板扩增结果，“m”指marker dl2000，“沙门”指肠炎沙门氏菌亚种，“小耶”指小肠结炎耶尔森菌亚种，“混合”指两种种细菌基因组dna的混合模板。图8各单菌扩增片段均符合理论计算值且混合模板种各菌扩增片段能彼此分开。
[0193]
实施例5同时检测沙门氏菌、金黄色葡萄球菌和小肠结炎耶尔森菌
[0194]
本发明针对沙门氏菌、金黄色葡萄球菌和小肠结炎耶尔森菌基因组筛选富含at的序列并设计特异性引物，根据计算的引物tm和目标序列tm筛选引物组合，设置反应条件进行核酸扩增，通过判断反应结果是否为阳性，确定待测样品中是否存在目标序列，进而确定待测样品中是否存在沙门氏菌、金黄色葡萄球菌和小肠结炎耶尔森菌菌中的一种或多种。具体步骤如下：
[0195]
(1)富含at碱基的目标序列的筛选：
[0196]
使用2019年8月5号从ncbi的ftp上下载的有完整基因组序列的细菌、古菌和病毒的数据，共2896个全基因组序列。设定集合a，其包含所有沙门氏菌基因组序列；设定集合b，其包含所有非沙门氏菌基因组序列。沙门氏菌基因组序列作为参考基因组，从基因组第一个碱基开始滑动宽度为1000bp的窗口，步长为50bp；每次滑动前计算一次序列at碱基含量，保留序列at碱基含量大于60％的区域作为候选目标序列。以上过程采用perl脚本实现。
[0197]
(2)特异性引物的设计：
[0198]
根据pcr引物的特性及本发明的需求，设定引物的at碱基含量为55％～75％、3’端稳定性δg《4、5’端稳定性δg《3及引物序列长度(20～36bp)等特性参数，同时设定单条引物不能产生发卡结构、自身不能发生相互作用等条件，以步骤(1)中的候选目标序列作为单条引物设计的候选序列，计算得出符合上述设定条件的单条引物。通过公式0.466
×
(gc的百分含量)
×
100+66.04-(450/引物序列长度)计算单条引物的tm值，并记录每条引物在目标序列上的位置、正负链信息(即来源于正链还是负链)以及引物的长度等信息。
[0199]
根据单条引物的位置信息进行引物配对，保留同时满足引物对tm差值《3℃以及引物对扩增区域在200～600bp两个条件的引物对作为候选引物对。
[0200]
使用比对软件bowtie，将上一步设计出的候选引物对中的每一条引物分别与集合a中的目标基因组序列和集合b中的非目标基因组进行序列比对。为了保证引物的通用性，当该单条引物与集合a中的目标序列进行比对时，使用参数设置
“‑
a-n 0”，即要求该单条引物与目标序列完全匹配；为了保证引物的特异性，当该单条引物与集合b中的非目标序列进行比对时，使用参数设置
“‑
a-n 3”，即允许该单条引物与非目标序列有不超过3个错配。系统输出满足条件的引物对，其中引物对的数量可预先设置，在本实施例中引物对的数量预先设置为100。以上高通量、自动化引物设计流程采用c、perl脚本实现。程序运行后在at富
集区设计出沙门氏菌扩增引物对100个。
[0201]
通过同样的方法分别设计出小肠结炎耶尔森菌扩增引物对100个。使用公式0.466
×
(gc的百分含量)
×
100+66.04-(450/引物序列长度)计算每条引物的tm并计算该引物对平均tm值作为反应退火温度；使用公式0.357
×
(gc的百分含量)
×
100+70.582-(990/扩增产物序列长度)计算每个引物对扩增区域的变性温度。
[0202]
以不同细菌扩增序列长度至少相差50bp、反应变性温度最多相差2℃、反应退火温度最多相差2℃为条件筛选引物组合。系统输出满足条件的引物组合，数量可预先设置，在本实施例中引物组合的数量预先设置为10。随机选择一个引物组进行有效性验证。所述引物组合的序列为：
[0203]
引物对a：
[0204]
sal-f1：5
’‑
tgggttgaaatagcccatta-3’(seq id no.1)
[0205]
sal-r1：5
’‑
gacgtgacacacttcgtttt-3’(seq id no.2)
[0206]
引物对b：
[0207]
yer-f1：5
’‑
atggaaaataacataatttctattaccgg-3’(seq id no.7)
[0208]
yer-r1:5
’‑
tctctgcgaataaccttgtg-3’(seq id no.8)
[0209]
引物对e：
[0210]
sta-f2：5
’‑
cctttcatctaaaaacctcca-3’(seq id no.9)
[0211]
sta-r2:5
’‑
gaaatggatgttttaaaagaagg-3’(seq id no.10)
[0212]
所述引物对a的at碱基含量分别为60％和55％，理论平均退火温度为63℃；所述引物对的理论扩增产物序列片段长度为352bp，at的百分含量为74.43％，理论变性温度为76.90℃。
[0213]
所述引物对b的at碱基含量分别为72.41％和55％，理论平均退火温度为63.5℃；所述引物对的理论扩增产物序列片段长度为242bp，at的百分含量为72.02％，理论变性温度为76.48℃。
[0214]
所述引物对e的at碱基含量分别为61.90％和69.57％，理论平均退火温度为61℃；所述引物对的理论扩增产物序列片段长度为595bp，at的百分含量为76.30％，理论变性温度为77.38℃。
[0215]
(4)核酸扩增反应及结果检测：
[0216]
核酸扩增反应体系配置如下表7所示。根据计算结果，采用81℃变性5秒，62℃退火5秒，72℃延伸20秒，推荐以上过程重复35次。检测对象清单见表2。
[0217]
在使用普通梯度pcr仪反应完成后通过对扩增产物进行琼脂糖凝胶电泳，根据电泳条带判断扩增结果是否为阳性，即待测样品中是否存在目标序列。
[0218]
表7.三种细菌核酸扩增反应体系
[0219]
体系体积(μl)终浓度taqr300mix151xsal-f1/r1(50μm)0.25+0.250.5μmyer-f1/r1(50μm)0.25+0.250.5μmsta-f2/r2(50μm)0.25+0.250.5μm100％dmso0.52％
模板(10ng)0/1/ddh2oupto25/
[0220]
图9展示了本发明针对肠炎沙门氏菌亚种、金黄色葡萄球菌亚种和小肠结炎耶尔森菌三种细菌核酸扩增反应的琼脂糖凝胶电泳结果。其中“n”为其对应的阴性模板扩增结果，“m”指marker dl2000，“肠亚”指肠炎沙门氏菌亚种，“小耶”指小肠结炎耶尔森菌亚种，“金亚”指金黄色葡萄球菌亚种，“混合”指三种细菌基因组dna的混合模板。图9各单菌扩增片段均符合理论计算值且混合模板种各菌扩增片段能彼此分开。
[0221]
本发明还提供了一种同时检测多种细菌的设备，所述设备包括：存储器、处理器和扩增模块；所述存储器上存储有计算机程序，当所述计算机程序被所述处理器执行时，实现本发明实施例1-4所述方法的步骤。
[0222]
参考文献
[0223]
周慧琦.(2014).基因组gc含量与碱基、密码子和氨基酸使用偏好的关系.(master master),电子科技大学,(67)
[0224]
本发明的保护内容不局限于以上实施例。在不背离发明构思的精神和范围下，本领域技术人员能够想到的变化和优点都被包括在本发明中，并且以所附的权利要求书为保护范围。