一种基于单核苷酸多态性的葡萄球菌种水平鉴定和菌株分型一体化方法与流程

本发明涉及药品安全微生物检测技术领域，具体涉及一种用于葡萄球菌种水平鉴定和菌株分型的引物组合物，以及一种葡萄球菌种水平鉴定和菌株分型一体化方法。

背景技术：

随着药品生产质量规范和国家药品质量标准的提升，药品生产质量控制中对微生物检验要求和鉴定水平越来越高。新实施的《中国药典》2015年版规定：微生物鉴定水平视情况而定，包括种、属鉴定和菌株分型；无菌试验结果阳性时，对检出微生物鉴定必要时需达到菌株水平；药品微生物实验室在过程控制中应具有通过对污染菌溯源分析，发现污染源的能力。

葡萄球菌属在自然界中分布十分广泛，传统微生物检验方法通常以“凝固酶阳性”作为判定葡萄球菌具有致病性的重要标准，如凝固酶阳性的金黄色葡萄球菌能够引起化脓感染、败血症、肺炎等疾病，是药品微生物检验中重要控制菌。但近些年研究发现，凝固酶阴性葡萄球菌(如表皮葡萄球菌、溶血葡萄球菌和人葡萄球菌等)也是人类疾病重要病原体。环境微生物监测结果表明凝固酶阴性葡萄球菌是污染药品、生产环境和检测实验室环境的主要微生物，占到50％以上。传统微生物检测方法仅关注“凝固酶阳性”葡萄球菌存在较大的安全风险，亟需建立能够准确鉴定多种致病性葡萄球菌的方法。基于16srrna序列比对的方法目前被广泛应用于细菌的鉴定，但有研究表明16srrna核酸序列在葡萄球菌属中高度保守，对于亲缘关系较近的葡萄球菌种的区分能力差[speciesidentificationofcoagulase-negativestaphylococci:genotypingissuperiortophenotyping_[j].veterinarymicrobiology,2009,134(1–2):20-28]，不能够提供准确的鉴定结果，因此亟需在葡萄球菌属中寻找新的区分能力强的靶基因用于葡萄球菌种水平的鉴定。

此外，《中国药典》2015年版指出药品微生物检验实验室应具有对微生物分型溯源的能力。多位点序列分型方法(multilocussequencetyping，mlst)是一种基于核酸序列测定的细菌分型方法。这种方法通过pcr扩增多个管家基因内部片段并测定其序列，利用管家基因序列之间的变异来区分不同菌株。mlst方法操作简单，能快速得到菌株分型结果，便于不同实验室进行数据比较分析，已经用于多种细菌的流行病学监测和进化研究。目前已经针对葡萄球菌属中的金黄色葡萄球菌[multilocussequencetypingforcharacterizationofmethicillin-resistantandmethicillin-susceptibleclonesofstaphylococcusaureus[j].jclinmicrobiol38(3):1008-1015]、表皮葡萄球菌[improvedmultilocussequencetypingschemeforstaphylococcusepidermidis[j].jclinmicrobiol45(2):616-619]、人葡萄球菌[multilocussequencetypingandfurthergeneticcharacterizationoftheenigmaticpathogen,staphylococcushominis[j].plosone8(6):e66496]等建立了相应的mlst方法，但是这些葡萄球菌mlst分型方法所用的靶基因和引物都具有种特异性，并不能用于其它葡萄球菌菌株的分型。另外，mlst方法的应用需建立在菌株鉴定的基础上，极大的降低了工作效率，提高了实验成本。

因此，本领域的技术人员致力于开发一套基于核酸序列差异分析的，能够同时实现多种葡萄球菌种水平鉴定和菌株分型的一体化方法。

技术实现要素：

本发明的目的在于克服现有技术中的缺陷，提供一种用于葡萄球菌种水平鉴定和菌株分型的靶基因及其对应的pcr扩增引物，及利用该pcr扩增引物对葡萄球菌种水平鉴定和菌株分型一体化的方法，该方法能够同时实现多种葡萄球菌种水平的鉴定和分型，提高了工作效率，降低了实验成本。

为实现上述目的，本发明采用如下技术方案：

本发明的第一个目的是提供一种用于葡萄球菌种水平鉴定和菌株分型的靶基因，其为具有单核苷酸多态性特征的葡萄球菌持家基因，其包括fema(aminoacyltransferase)、ftsz(celldivisionprotein)、gap(glyceraldehyde-3-phosphatedehydrogenase)、pyrh(uridylatekinase)、rpob(rnapolymeraseb)和tufa(translationalelongationfactor)。

本发明之所以选择fema、ftsz、gap、pyrh、rpob和tufa等6个持家基因作为葡萄球菌种水平鉴定和菌株分型靶标是通过以下方法筛选到的：

从ncbi网站上下载57株代表性葡萄球菌全基因组序列，涵盖28个葡萄球菌种。然后进行泛基因组学分析，通过以下筛选条件的设置：1)核心基因(coregenes)，2)靶基因编码的氨基酸序列在葡萄球菌之间的同源性大于50％，3)靶基因在葡萄球菌基因组中必须是单拷贝基因，4)氨基酸序列长度大于200个氨基酸。进行泛基因组学数据整理分析，对筛选出的靶基因序列逐一进行序列分析，剔除有明显片段缺失和没有明显保守区段的靶基因序列，最后筛选出具有单核苷酸多态性特征的葡萄球菌持家基因。

优选的，上述各靶基因的标准分析序列为seqidno：1～seqidno：6所示序列。

本发明的第二个目的是提供一种用于葡萄球菌种水平鉴定和菌株分型的引物组合物，其包括：

fema引物对：

fema-ftaayaaartcaccaacataytcseqidno：7

fema-rtymgntcatttatggaagatacseqidno：8

ftsz引物对：

ftsz-fgcbacdytdaargtyatyggseqidno：9

ftsz-rttrtcwtcraadccwgthgcseqidno：10

gap引物对：

gap-fatggttttggtagaattggtcgtttaseqidno：11

gap-rgacatttcgttatcataccaagctgseqidno：12

pyrh引物对：

pyrh-faytdagygghgaagcdytmgseqidno：13

pyrh-ragaargangangchgtwgaseqidno：14

rpob引物对：

rpob-fcaattcatggaccaagcseqidno：15

rpob-rccgtcccaagtcatgaaacseqidno：16

tufa引物对：

tufa-fccaatgccacaaactcgtseqidno：17

tufa-rcctgaaccaacagtacgtseqidno：18。

上述6个靶基因的pcr扩增引物是通过对ncbi数据库中下载的葡萄球菌靶基因序列，序列比对分析后，利用primer5.0软件在靶基因保守序列区域设计引物，然后对引物的有效性和通用性进行评价，最终筛选出6对性能较优的引物。所述6对pcr扩增引物序列，在葡萄球菌中具有较好的通用性，能够实现葡萄球菌属多种葡萄球菌的有效扩增。

本发明第三个目的是提供一种葡萄球菌种水平鉴定和菌株分型一体化的方法，其包括以下步骤：

步骤1、提取待检测菌株dna；

步骤2、将步骤1中提取的待检测菌株dna作为模板，利用上述引物组合物分别对6种靶基因：fema、ftsz、gap、pyrh、rpob和tufa进行pcr扩增；

步骤3、通过凝胶电泳，检测步骤2中所得的扩增产物，定性判断待检测菌株是否为葡萄球菌；

步骤4、将步骤3中定性判断为葡萄球菌的产物进行核酸序列测定，

步骤5、将步骤4中的测定的6个靶基因的序列与对应靶基因的标准序列进行序列比对分析，截取片段大小一致的序列；

步骤6、将步骤5中截取的待检测菌株靶基因序列按fema-ftsz-gap-pyrh-rpob-tufa顺序进行序列串联；

步骤7、将步骤6中所得的待检测菌株串联序列与葡萄球菌菌种数据库中的标准菌株串联序列进行聚类分析，构建邻接法(neighbor-joiningmethod)进化树，根据邻接法进化树中菌株之间的亲缘聚类关系，进行葡萄球菌种水平鉴定和菌株分型。

为了进一步优化上述技术方案，本发明所采取的技术措施还包括：

优选地，在步骤1和步骤2之间，包括筛选具有单核苷酸多态性的葡萄球菌持家基因fema、ftsz、gap、pyrh、rpob和tufa，作为葡萄球菌种水平鉴定和菌株分型的靶基因的步骤。

优选地，所述筛选方法为葡萄球菌泛基因组学分析方法，其筛选条件如下：1)核心基因(coregenes)，2)靶基因编码的氨基酸序列在葡萄球菌之间的同源性大于50％，3)靶基因在葡萄球菌基因组中必须是单拷贝基因，4)氨基酸序列长度大于200个氨基酸。

优选地，所述葡萄球菌泛基因组学分析方法中包括57株代表性葡萄球菌全基因组序列，涵盖28个葡萄球菌种。

优选地，所述葡萄球菌泛基因组学分析方法包括：

i)对氨基酸序列的比对分析，按同源性从高到低逐一进行序列分析，剔除有明显片段缺失和没有明显保守区段的氨基酸靶序列；

ii)对筛选出的氨基酸靶序列的核酸序列进行比对分析，剔除没有明显保守区段的靶基因序列，筛选出具有单核苷酸多态性特征的葡萄球菌持家基因。

pcr的反应过程是本领域技术人员所能掌握的，一般为五步法，预变性、变性、退火、延伸、再延伸。根据本发明人研究，优选地，步骤2中的pcr扩增反应体系包括：takara2×premix预混液12.5μl，引物对(10μmol/l)各1.0μl，菌株dna模板2.0μl，使用灭菌ddh2o水补足20.0μl。

优选地，步骤2中的pcr扩增程序如下：

fema：95℃预变性5min，95℃预变性30s，48℃退火30s，72℃延伸30s，进行35个循环，最后72℃延伸10min；

ftsz：95℃预变性5min，95℃预变性30s，50℃退火30s，72℃延伸45s，进行35个循环，最后72℃延伸10min；

gap：95℃预变性5min，95℃预变性30s，50℃退火30s，72℃延伸45s，进行35个循环，最后72℃延伸10min；

pyrh：95℃预变性5min，95℃预变性30s，50℃退火30s，72℃延伸45s，进行35个循环，最后72℃延伸10min；

rpob：95℃预变性5min，95℃预变性30s，50℃退火30s，72℃延伸45s，进行35个循环，最后72℃延伸10min；

tufa：95℃预变性5min，95℃预变性30s，50℃退火30s，72℃延伸45s，进行35个循环，最后72℃延伸10min；

优选地，所述步骤3中定性判断待检测菌株的标准为：如果凝胶电泳结果中出现6个靶基因的明显亮色条带，则表明该条带对应的菌株为葡萄球菌菌株，否则为非葡萄球菌菌株。

优选地，所述步骤3中凝胶电泳可采用本领域常用的凝胶电泳类型，更优选为琼脂糖凝胶电泳，琼脂糖凝胶电泳操作简单，电泳速度快，样品不需事先处理就可以进行电泳，且对样品吸附极微，因此电泳图谱清晰，分辨率高，重复性好。

本发明收集整理国家权威微生物菌种保藏中心葡萄球菌标准菌株，及经过明确鉴定的微生物检验实验室常见葡萄球菌菌株。进行靶基因pcr扩增和核酸测序分析后，建立不同葡萄球菌种靶基因fema、ftsz、gap、pyrh、rpob和tufa的标准序列及串联序列数据库。将待检测菌株的6个靶基因串联序列与数据库中不同葡萄球菌标准菌株进行比对和聚类分析，进行葡萄球菌菌种水平鉴定和菌株分型。

本发明分别在上述6个靶基因的snp多态性丰富的区域定义每个靶基因的标准分析序列(约500bp)，其序列如下所示：

fema标准分析序列(seqidno：1)：

tcttcggcatcttcactaaagtcaccactaataccatagaaattataccgatcaataccatgttcaattgcatagttaatcatcttccattgaaccgcatagctccctgcaaaatggcgataacgatttgaagttccaccagcgtagtaaactacttcaaacggattgattataaagaagccagcagagatgggtaattcattgccatgttcttgttttaagtttttggcttcattaattttttgctggtttgcatcaagttgttgttctaaattttcttttttgttatgtgcttttttattctctggacgcttctcaatgtcttttaaagctttattataatctttattaagcacatttctttcattatttagttcctctatgtactcatcgaagttaatataggctagtggtactaaaacacggtctttatagtgtttgaatctgttgtaataaaaactatcatctctatctgcaaaatctttagtttcagaggta；

ftsz标准分析序列(seqidno：2)：

acacgtggattaggtgctggtgctaatcctgaaattggaaagaaagcagcagaagaatcaagagaacaaattgaagatgctatccaaggtgctgatatggtattcgtaactgctggtatgggtggcggtactggaacaggtgctgcgccggttgttgctaaaattgcaaaagaaatgggtgctttaactgtaggtgttgttacgcgtccattcggtttcgaaggtcgtaagcgccaaacacaagcggcagctggtgtagaatctatgaaagcagcagtggatacattaattgttattccaaatgatcgcttattagatatcgttgacaaatctacgccgatgatggaagcatttaaagaagcggataatgtattacgtcaaggtgtacaaggtatttcagatttaattgcagtatcaggtgaagtgaacttagactttgcagacgttaaaacaattatgtctaatcaaggttctgcgttaatgggtatcggtgtgtcatctggtgaaaacagagcagtcgaa；

gap标准分析序列(seqidno：3)：

aggtcgtttcactggagaagttgaagttatcgaaggtggtttccgtgtaaacggtaaagaaattaaatcattcgatgaaccagatgctggtaaattaccatggggcgatttagatatcgacgtagtattagaatgtactggtttctatactgataaagaaaaagcacaagctcacatcgatgcaggtgctaaaaaagtattaatctccgctccagctaaaggtgatgtaaaaacaatcgtattcaacactaaccatgatacattagatggttcagaaacagttgtttcaggtgcttcttgtactactaactcattagcaccagttgcaaaagttttaagtgacgaattcggtttagttgaaggtttcatgactacaattcacgcttacactggtgaccaaaatacacaagacgcacctcacagaaaaggtgacaaacgtcgtgcacgtgcagcagcagaaaatattattcctaactcaacaggtgctgctaaagctatcggtaaagttattccagaaatcgatggtaaa；

pyrh标准分析序列(seqidno：4)：

ttgaatatgtgtcaaatgttcatatttaattgcatttgcatcaacttttggatcagctgaataaacaccatctacattatttttacccattaaaattacatctgcttctacttcagccgcacgtaaagcagcagttgtgtcagttgagaagtaaggattaccaattccagcggcaaatataacaacacgttttttctccaaatgtctaattgcacgtcgacgaatatatggttcagcaacttgtttcatttcaatagaagttaaaacacgtgtatcacaatctaattgttctaaactgtcttgtaatgctaaagcattcattacagtagcaagcatccccatgtaatcagccgtaccacgatccatacctaggtcgctaccagttttgcctctccaaatatttccaccgccgacaatgacagcaatttcacaatccattttggcaacttcagcaacttgctgtgcaacgctttt；

rpob标准分析序列(seqidno：5)：

taacaagtaccagagtttgaacgcttgaatttcgctaaaggataacgatctaattcaccttcatgctcagttccattttcttcaactaaacgacgaactaaaatttcattagattcaacatgttcaacgcgccctctatgcttagcagtgattgcagcaccagagtctcttgcggctacatgttccatacctgtacccacaaatggagcttcaggattcattaaaggcactgcttgacgttgcatgtttgctcccattaacgcacggttagagtcgtcattttctaagaatggaatacatgctgttgctgctgaaacaacttgttttggtgaaacgtccatgtaatccattttctctttagccataacagtgttattaccacggaaacgacaaacaacttcatcatctaagaaacgtccattttcatcaagtctagaattagcctgtgcaacaacgtagctatcctcttcatcagctgtcaaataatctatttgatcggtgattgagtttgtatctaaatccactttacgatatgg；

tufa标准分析序列(seqidno：6)：

cccaggtgacgatgtacctgtaatcgctggttctgcattaaaagcattagaaggcgatgctgaatacgaacaaaaaatcttagacttaatgcaagcagttgatgattacattccaactccagaacgtgattctgacaaaccattcatgatgccagttgaggacgtattctcaatcactggtcgtggtactgttgctacaggccgtgttgaacgtggtcaaatcaaagttggtgaagaagttgaaatcatcggtatgcatgaaacttc---taaaacaactgttactggtgtagaaatgttccgtaagttattagactacgctgaagctggtgacaacatcggtgctttattacgtggtgttgcacgtgaagacgtacaacgtggtcaagtattagctgctcctggttctattacaccacacacaaaattcaaagctgaagtatacgtattatctaaagatgaaggtggacgtcacactccattcttcactaactatcgccc。

与现有技术相比，本发明具有以下有益效果：

本发明采用以上优化方案，优选6个靶基因，基于靶基因的单核苷酸多态性建立了葡萄球菌种水平鉴定和菌株分型一体化方法，该方法可对国内常见的葡萄球菌种能够快速有效的鉴定，鉴定同时能够进行菌株分型。本发明所述的葡萄球菌鉴定及分型的一体化方法，具有通用性高、鉴定准确和分辩力强等优点，且通量高、操作简单、快速有效。解决了现有分型方法需在菌株鉴定的基础上才能实现菌株分型的弊端，实现了葡萄球菌种水平鉴定和分型的一体化，降低了试验成本，提高了工作效率，结果判定方便准确，更加具有实用性。

本发明所述的葡萄球菌鉴定分型一体检测鉴定方法提供了新的技术平台，可用于药品质量监管部门、药品生产企业，对药品、生产和检验环境中污染葡萄球菌的种水平鉴定和菌株分型溯源，具有广阔市场前景和较大的经济、社会效益。

附图说明

图1是靶基因pcr扩增产物琼脂糖凝胶电泳检测图谱；

图2是172株葡萄球菌6个靶基因串联标准序列聚类分析图。

具体实施方式

本发明提供了一种用于葡萄球菌种水平鉴定和菌株分型的靶基因及其对应的引物序列，以及葡萄球菌种水平鉴定和菌株分型一体化方法。

下面结合附图和实施例，对本发明的具体实施方式作进一步描述。以下实施例仅用于更加清楚地说明本发明的技术方案，而不能以此来限制本发明的保护范围。

本发明使用的葡萄球菌菌株包含国内常见的18种葡萄球菌，具体菌株信息如表1所示。共172株葡萄球菌，分别来自于国家权威微生物菌种保藏中心，及经过明确鉴定的微生物检验实验室。为保证菌株信息的准确性，我们分别使用16srrna序列比对和全自动微生物生化鉴定系统vitek2.0compact，对菌株信息进行了再确认。

实施例1

本实施例为用于葡萄球菌种水平鉴定和菌株分型的靶基因pcr扩增方法的建立，其包括如下步骤：

步骤一：筛选葡萄球菌属具有单核苷酸多态性的持家基因；

从ncbigenome网站下载57株高质量的代表性葡萄球菌全基因组核酸序列，涵盖28个葡萄球菌种，进行葡萄球菌泛基因组学分析。其中，以氨基酸序列作为分析对象，选择orthomcl软件进行泛基因组学分析，利用mysql软件进行数据管理，采用mcl软件算法进行序列之间的比对分析。通过以下筛选条件的设置：1)核心基因(coregenes)，在所有葡萄球菌菌株中都具有的基因，2)靶基因编码的氨基酸序列在葡萄球菌之间的同源性大于50％，提供一定的保守性和足够的变异性，以区分不同的葡萄球菌菌株，3)靶基因在葡萄球菌基因组中必须是单拷贝基因，保证基因的唯一性，4)氨基酸序列长度大于200个氨基酸，便于设计引物并提供足够的序列多样性。通过对葡萄球菌泛基因组学数据进行整理，筛选出具有单核苷酸多态性的葡萄球菌持家基因序列。利用软件mega7.0进行氨基酸序列的比对分析，按同源性从高到低逐一进行序列分析。首先剔除有明显片段缺失和没有明显保守区段的氨基酸靶序列。然后再针对筛选出的氨基酸靶序列的核酸序列进行比对分析，剔除没有明显保守区段的靶基因序列，最后筛选出具有单核苷酸多态性特征的葡萄球菌持家基因。

步骤二：靶基因pcr扩增方法的评价；

靶基因pcr扩增引物设计：靶基因核酸序列比对分析后，尽量避开单核苷酸多态性位点，选择保守核酸序列区域，利用引物设计软件primer5.0，针对筛选出的靶基因设计多对pcr扩增引物。

靶基因pcr扩增有效性评价：选择表1中的172株葡萄球菌，提取dna后，在优化的退火温度下，采用seqidno：1-seqidno：12的引物对分别对候选靶基因的pcr扩增，验证靶基因pcr扩增引物的有效性和通用性。琼脂糖凝胶电泳检测pcr扩增产物，阳性pcr扩增产物进行核酸测序。对测序后获得的靶基因序列进行分析，进一步判定pcr扩增的准确性。

表1.172株葡萄球菌菌株信息表

如图1所示，结果表明候选靶基因fema，ftsz，gap，pyrh，rpob和tufa的pcr扩增引物，在172株葡萄球菌中都能够进行有效pcr扩增(图1)，且核酸测序结果表明6个靶基因的pcr扩增序列与目标序列一致。以上结果表明，针对具有单核苷酸多态性的靶基因fema，ftsz，gap，pyrh，rpob和tufa，成功建立了pcr扩增方法，且靶基因pcr扩增引物在葡萄球菌中具有非常好的扩增效率和适用性。

上述pcr反应体系如下：采用20.0μl的反应体系，其中包括：takara2×premix预混液12.5μl，上下游引物(10μmol/l)各1.0μl，菌株dna模板2.0μl，使用灭菌ddh2o水补足20.0μl。

上述pcr反应循环参数如下：

fema:95℃预变性5min，95℃预变性30s，48℃退火30s，72℃延伸30s，进行35个循环，最后72℃延伸10min；

ftsz:95℃预变性5min，95℃预变性30s，50℃退火30s，72℃延伸45s，进行35个循环，最后72℃延伸10min；

gap:95℃预变性5min，95℃预变性30s，50℃退火30s，72℃延伸45s，进行35个循环，最后72℃延伸10min；

pyrh:95℃预变性5min，95℃预变性30s，50℃退火30s，72℃延伸45s，进行35个循环，最后72℃延伸10min；

rpob:95℃预变性5min，95℃预变性30s，50℃退火30s，72℃延伸45s，进行35个循环，最后72℃延伸10min；

tufa:95℃预变性5min，95℃预变性30s，50℃退火30s，72℃延伸45s，进行35个循环，最后72℃延伸10min。

实施例2

本实施例为葡萄球菌种水平鉴定和菌株分型方法的建立，其包括如下步骤：

步骤一：靶基因序列标准化及多样性分析；

对实施例1中获取得172株葡萄球菌中6个靶基因的核酸序列，利用软件mega7.0进行比对分析，去除两端准确度较低的核酸序列，在snp多态性丰富的区域定义每个靶基因的标准分析序列(约500bp)，其序列为seqidno：1-seqidno：6所示序列。

利用软件dnasp5.0分析每个靶基因标准序列的多样性，分别从等位基因(alleles)数目、snp位点数、nucleotidediversity(π)和allelicdiversity(hd)等方面，评价靶基因在不同葡萄球菌之间的序列多样性(如表2所示)。

6个靶基因在172株菌中的等位基因数目分别为fema(53)，ftsz(50)，gap(51)，pyrh(48)，rpob(56)和tufa(41)，snp位点在每个靶基因中的比率都达到了30％以上，其中靶基因fema的snp位点数最多，占到了60％以上。综合以上分析结果，本发明筛选的靶基因标准序列在不同葡萄球菌之间具有较好的多样性，提供了区分不同葡萄球菌菌株的能力。

步骤二：建立葡萄球菌鉴定和分型方法；

本发明选取每个菌株6个靶基因的标准序列，按照fema-ftsz-gap-pyrh-rpob-tufa顺序进行标准序列的串联，串联后利用mega7.0进行比对分析。将定性判断为葡萄球菌的靶基因pcr产物进行序列测定，并将其与对应靶基因的标准序列进行比对分析，截取片段大小一致的序列；利用上述截取的序列按fema-ftsz-gap-pyrh-rpob-tufa顺序进行序列串联。

基于序列比对结果，对所得的待检测菌株的串联序列与葡萄球菌菌种数据库中的标准菌株串联序列，选用邻接法(neighbor-joiningmethod)和bootstrapmethod(1000)方法进行聚类分析，其结果如图2和表2所示。

表2.172株葡萄球菌菌种鉴定与分型结果

由如图2所示，对属于同一个葡萄球菌种的菌株标记为同一种颜色符号。聚类结果表明，本发明建立的基于单核苷酸多态性分析的葡萄球菌种水平鉴定和分型一体化方法，可以准确的把172株葡萄球菌鉴定为18个种，区分为95个序列型，其结果如表2所示。

由上述结果分析可知本发明依据172株葡萄球菌的鉴定和分型结果，建立葡萄球菌靶基因标准序列及串联序列数据库，用于鉴定待检葡萄球菌，并对葡萄球菌菌株进行分型溯源。

由上述实施例可知，解决了现有分型方法需在菌株鉴定的基础上才能实现菌株分型的弊端，实现了葡萄球菌种水平鉴定和菌株分型的一体化，降低了试验成本，提高了工作效率，结果判定方便准确，更加具有实用性

以上对本发明的具体实施例进行了详细描述，但其只作为范例，本发明并不限制于以上描述的具体实施例。对于本领域技术人员而言，任何对本发明进行的等同修改和替代也都在本发明的范畴之中。因此，在不脱离本发明的精神和范围下所作的均等变换和修改，都应涵盖在本发明的范围内。

sequencelisting

<110>上海市食品药品检验所

<120>一种基于单核苷酸多态性的葡萄球菌种鉴定和菌株分型一体化方法

<160>18

<170>patentinversion3.5

<210>1

<211>498

<212>dna

<213>fema标准分析序列

<400>1

tcttcggcatcttcactaaagtcaccactaataccatagaaattataccgatcaatacca60

tgttcaattgcatagttaatcatcttccattgaaccgcatagctccctgcaaaatggcga120

taacgatttgaagttccaccagcgtagtaaactacttcaaacggattgattataaagaag180

ccagcagagatgggtaattcattgccatgttcttgttttaagtttttggcttcattaatt240

ttttgctggtttgcatcaagttgttgttctaaattttcttttttgttatgtgctttttta300

ttctctggacgcttctcaatgtcttttaaagctttattataatctttattaagcacattt360

ctttcattatttagttcctctatgtactcatcgaagttaatataggctagtggtactaaa420

acacggtctttatagtgtttgaatctgttgtaataaaaactatcatctctatctgcaaaa480

tctttagtttcagaggta498

<210>2

<211>522

<212>dna

<213>ftsz标准分析序列

<400>2

acacgtggattaggtgctggtgctaatcctgaaattggaaagaaagcagcagaagaatca60

agagaacaaattgaagatgctatccaaggtgctgatatggtattcgtaactgctggtatg120

ggtggcggtactggaacaggtgctgcgccggttgttgctaaaattgcaaaagaaatgggt180

gctttaactgtaggtgttgttacgcgtccattcggtttcgaaggtcgtaagcgccaaaca240

caagcggcagctggtgtagaatctatgaaagcagcagtggatacattaattgttattcca300

aatgatcgcttattagatatcgttgacaaatctacgccgatgatggaagcatttaaagaa360

gcggataatgtattacgtcaaggtgtacaaggtatttcagatttaattgcagtatcaggt420

gaagtgaacttagactttgcagacgttaaaacaattatgtctaatcaaggttctgcgtta480

atgggtatcggtgtgtcatctggtgaaaacagagcagtcgaa522

<210>3

<211>529

<212>dna

<213>gap标准分析序列

<400>3

aggtcgtttcactggagaagttgaagttatcgaaggtggtttccgtgtaaacggtaaaga60

aattaaatcattcgatgaaccagatgctggtaaattaccatggggcgatttagatatcga120

cgtagtattagaatgtactggtttctatactgataaagaaaaagcacaagctcacatcga180

tgcaggtgctaaaaaagtattaatctccgctccagctaaaggtgatgtaaaaacaatcgt240

attcaacactaaccatgatacattagatggttcagaaacagttgtttcaggtgcttcttg300

tactactaactcattagcaccagttgcaaaagttttaagtgacgaattcggtttagttga360

aggtttcatgactacaattcacgcttacactggtgaccaaaatacacaagacgcacctca420

cagaaaaggtgacaaacgtcgtgcacgtgcagcagcagaaaatattattcctaactcaac480

aggtgctgctaaagctatcggtaaagttattccagaaatcgatggtaaa529

<210>4

<211>474

<212>dna

<213>pyrh标准分析序列

<400>4

ttgaatatgtgtcaaatgttcatatttaattgcatttgcatcaacttttggatcagctga60

ataaacaccatctacattatttttacccattaaaattacatctgcttctacttcagccgc120

acgtaaagcagcagttgtgtcagttgagaagtaaggattaccaattccagcggcaaatat180

aacaacacgttttttctccaaatgtctaattgcacgtcgacgaatatatggttcagcaac240

ttgtttcatttcaatagaagttaaaacacgtgtatcacaatctaattgttctaaactgtc300

ttgtaatgctaaagcattcattacagtagcaagcatccccatgtaatcagccgtaccacg360

atccatacctaggtcgctaccagttttgcctctccaaatatttccaccgccgacaatgac420

agcaatttcacaatccattttggcaacttcagcaacttgctgtgcaacgctttt474

<210>5

<211>536

<212>dna

<213>rpob标准分析序列

<400>5

taacaagtaccagagtttgaacgcttgaatttcgctaaaggataacgatctaattcacct60

tcatgctcagttccattttcttcaactaaacgacgaactaaaatttcattagattcaaca120

tgttcaacgcgccctctatgcttagcagtgattgcagcaccagagtctcttgcggctaca180

tgttccatacctgtacccacaaatggagcttcaggattcattaaaggcactgcttgacgt240

tgcatgtttgctcccattaacgcacggttagagtcgtcattttctaagaatggaatacat300

gctgttgctgctgaaacaacttgttttggtgaaacgtccatgtaatccattttctcttta360

gccataacagtgttattaccacggaaacgacaaacaacttcatcatctaagaaacgtcca420

ttttcatcaagtctagaattagcctgtgcaacaacgtagctatcctcttcatcagctgtc480

aaataatctatttgatcggtgattgagtttgtatctaaatccactttacgatatgg536

<210>6

<211>498

<212>dna

<213>tufa标准分析序列

<400>6

cccaggtgacgatgtacctgtaatcgctggttctgcattaaaagcattagaaggcgatgc60

tgaatacgaacaaaaaatcttagacttaatgcaagcagttgatgattacattccaactcc120

agaacgtgattctgacaaaccattcatgatgccagttgaggacgtattctcaatcactgg180

tcgtggtactgttgctacaggccgtgttgaacgtggtcaaatcaaagttggtgaagaagt240

tgaaatcatcggtatgcatgaaacttctaaaacaactgttactggtgtagaaatgttccg300

taagttattagactacgctgaagctggtgacaacatcggtgctttattacgtggtgttgc360

acgtgaagacgtacaacgtggtcaagtattagctgctcctggttctattacaccacacac420

aaaattcaaagctgaagtatacgtattatctaaagatgaaggtggacgtcacactccatt480

cttcactaactatcgccc498

<210>7

<211>22

<212>dna

<213>artificialsequence

<220>

<223>fema-f

<220>

<221>misc_feature

<222>(4)..(4)

<223>yiscort

<220>

<221>misc_feature

<222>(8)..(8)

<223>risaorg

<220>

<221>misc_feature

<222>(20)..(20)

<223>yiscort

<400>7

taayaaartcaccaacataytc22

<210>8

<211>22

<212>dna

<213>artificialsequence

<220>

<223>fema-r

<220>

<221>misc_feature

<222>(2)..(2)

<223>yiscort

<220>

<221>misc_feature

<222>(3)..(3)

<223>misaorc

<220>

<221>misc_feature

<222>(5)..(5)

<223>nisa,c,g,ort

<400>8

tymgntcatttatggaagatac22

<210>9

<211>20

<212>dna

<213>artificialsequence

<220>

<223>ftsz-f

<220>

<221>misc_feature

<222>(3)..(3)

<223>bisg,t,orc

<220>

<221>misc_feature

<222>(6)..(6)

<223>disg,a,ort

<220>

<221>misc_feature

<222>(7)..(7)

<223>yiscort

<220>

<221>misc_feature

<222>(9)..(9)

<223>disg,a,ort

<220>

<221>misc_feature

<222>(12)..(12)

<223>risaorg

<220>

<221>misc_feature

<222>(15)..(15)

<223>yiscort

<220>

<221>misc_feature

<222>(18)..(18)

<223>yiscort

<400>9

gcbacdytdaargtyatygg20

<210>10

<211>20

<212>dna

<213>artificialsequence

<220>

<223>ftsz-r

<220>

<221>misc_feature

<222>(3)..(3)

<223>risaorg

<220>

<221>misc_feature

<222>(6)..(6)

<223>wisaort

<220>

<221>misc_feature

<222>(9)..(9)

<223>risaorg

<220>

<221>misc_feature

<222>(12)..(12)

<223>disg,aort

<220>

<221>misc_feature

<222>(15)..(15)

<223>wisaort

<220>

<221>misc_feature

<222>(18)..(18)

<223>hisa,torc

<400>10

ttrtcwtcraadccwgthgc20

<210>11

<211>26

<212>dna

<213>artificialsequence

<220>

<223>gap-f

<400>11

atggttttggtagaattggtcgttta26

<210>12

<211>25

<212>dna

<213>artificialsequence

<220>

<223>gap-r

<400>12

gacatttcgttatcataccaagctg25

<210>13

<211>20

<212>dna

<213>artificialsequence

<220>

<223>pyrh-f

<220>

<221>misc_feature

<222>(2)..(2)

<223>yiscort

<220>

<221>misc_feature

<222>(4)..(4)

<223>disg,aort

<220>

<221>misc_feature

<222>(7)..(7)

<223>yiscort

<220>

<221>misc_feature

<222>(10)..(10)

<223>hisa,torc

<220>

<221>misc_feature

<222>(16)..(16)

<223>disg,aort

<220>

<221>misc_feature

<222>(17)..(17)

<223>yiscort

<220>

<221>misc_feature

<222>(19)..(19)

<223>misaorc

<400>13

aytdagygghgaagcdytmg20

<210>14

<211>19

<212>dna

<213>artificialsequence

<220>

<223>pyrh-r

<220>

<221>misc_feature

<222>(5)..(5)

<223>risaorg

<220>

<221>misc_feature

<222>(8)..(8)

<223>nisa,c,g,ort

<220>

<221>misc_feature

<222>(11)..(11)

<223>nisa,c,g,ort

<220>

<221>misc_feature

<222>(14)..(14)

<223>hisa,torc

<220>

<221>misc_feature

<222>(17)..(17)

<223>wisaort

<400>14

agaargangangchgtwga19

<210>15

<211>17

<212>dna

<213>artificialsequence

<220>

<223>rpob-f

<400>15

caattcatggaccaagc17

<210>16

<211>19

<212>dna

<213>artificialsequence

<220>

<223>rpob-r

<400>16

ccgtcccaagtcatgaaac19

<210>17

<211>18

<212>dna

<213>artificialsequence

<220>

<223>tufa-f

<400>17

ccaatgccacaaactcgt18

<210>18

<211>18

<212>dna

<213>artificialsequence

<220>

<223>tufa-r

<400>18

cctgaaccaacagtacgt18