本发明属于生物鉴别和/或区分领域，具体涉及琥珀酸弧菌属在鉴别和/或区分不同种族个体中的应用。
背景技术：
：：人体肠道中存在大量的共生菌群，其数量超过1000万亿个，大约是人体细胞总数的10倍。同时肠道内微生物基因的数量约为300万个，大约是人类基因组基因数量的100多倍，如此海量的基因能够帮助微生物适应多变的环境，形成了与人体密不可分的共生关系。不同遗传背景和生活习惯的居民其肠道菌群组成存在较大差异，例如拟杆菌属是西方居民肠道中相对含量最多且个体间差异最大的细菌属，韩国居民是栖粪杆菌属。目前我国不同种族肠道菌群结构差异对比的研究数据尚不足。琥珀酸弧菌属属变形杆菌门，琥珀酸弧菌科，可代谢产生乙酸和琥珀酸，参与维生素等抗炎物质的合成，有潜在的抗炎作用。迄今为止，尚未见粪便中琥珀酸弧菌属相对丰度检测应用于鉴别和/或区分人群种族的报道。专利cn107746874a中公开了一种精神分裂症的标志物，所述的标志物包括琥珀酸弧菌属等。现有技术中并没有公开琥珀酸弧菌属在鉴别和/或区分人群种族中的应用。技术实现要素：第一方面，本发明提供一种微生物标志物在鉴别和/或区分不同种族个体中的应用，所述的微生物标志物为琥珀酸弧菌属。所述的琥珀酸弧菌属包括溶糊精琥珀酸弧菌(succinivibriodextrinosolvens)。所述不同种族为藏族和汉族，优选为高原上的藏族和汉族。所述的高原为海拔2000m以上，具有低压、缺氧的条件。进一步优选的，所述的高原环境为海拔2700m以上，具有低压、缺氧的条件下。在本发明的一个具体实施方式中，所述的高原环境为海拔3500m以上，具有低压、缺氧的条件下。在本发明的一个具体实施方式中，所述的高原环境为海拔5500m以上，具有低压、缺氧的条件下。第二方面，本发明提供一种鉴别和/或区分不同种族个体的方法，所述的方法的步骤包括：(1)确定所述的个体的粪便或肠道内容物中的微生物标志物的丰度；(2)将微生物标志物的丰度值与阈值比较，得知个体的种族。优选的，所述的步骤(1)中利用测序法测定微生物标志物的丰度。更优选的，所述的步骤(1)中包括：1)收集待检测个体的粪便或肠道内容物，并提取粪便或肠道内容物中微生物的dna；2)目的片段的扩增；3)测序；4)计算微生物标志物的丰度。所述的步骤1)中粪便或肠道内容物中微生物的dna的提取方法选自：ctab法、gitc法或利用商业化的试剂盒进行。优选的，利用ctab法提取粪便或肠道内容物中微生物的dna。具体的，将收集到的粪便或肠道内容物用ctab、溶菌酶和裂解液处理后，用苯酚、氯仿和异戊醇处理使dna游离及去除杂蛋白、利用异丙醇使dna析出，用乙醇纯化dna后得到粪便或肠道内容物中微生物的dna。更具体的，收集待检测个体的粪便，吸取500-1000ulctab裂解液至ep管里，加入溶菌酶，将200-500ul左右的样品加入裂解液中，50-70℃水浴条件下反应，离心取上清，加苯酚(ph8.0)：氯仿：异戊醇(20-30：18-28：1)，颠倒混匀，10000-12000rpm离心10-15min。取上清，加氯仿：异戊醇(20-30:1)，颠倒混匀，10000-12000rpm离心10-15min。取上清，加入异丙醇，上下摇晃，-20℃沉淀。10000-12000rpm离心10-15min，将沉淀用乙醇洗涤2-5次，将沉淀晾干，加入水溶解沉淀，加rnasea消化rna后既得粪便或肠道内容物中微生物的dna。所述的步骤2)中使用带barcode的特异引物515f和806r对步骤(1)中获得的dna进行扩增，所述的515f的序列为5’-gtgccagcmgccgcggtaa-3’(seqidno：1)，所述的806r的序列为5’-ggactachvgggtwtctaat-3’(seqidno：2)。所述的步骤3)中是利用高通量测序法对步骤(2)中获得的产物进行测序。所述的步骤4)中将测序得到的数据聚类成为otus(operationaltaxonomicunits)，对物种注释分析后，计算微生物标志物的丰度。具体的，所述的步骤4)中将步骤3)中测序得到的数据进行处理，利用uparse软件(uparsev7.0.1001，http://www.drive5.com/uparse/)对有效数据进行聚类成为otus(operationaltaxonomicunits)，用mothur方法与silva(http://www.arb-silva.de/)的ssurrna数据库对otus序列进行物种注释分析，以(微生物标志物otu值)/(检测样本检测到的所有菌属otu值)表示微生物标志物相对丰度。所述的步骤(2)中将测得的微生物标志物的丰度与阈值比较，大于阈值为一个种族，小于阈值为另一个种族。所述的阈值是通过对已知种族的人群进行微生物标志物的丰度的检测后所建立的受试者工作特征曲线的约登指数为最大值时得到的微生物标志物的丰度。优选的，本发明提供一种鉴别和/或区分个体为藏族或汉族的方法，所述的方法的步骤包括：(1)确定所述的个体的粪便或肠道内容物中的琥珀酸弧菌属的丰度；(2)将琥珀酸弧菌属的丰度与阈值比较，所述的阈值为5.62×10-5，当个体的粪便或肠道内容物中的琥珀酸弧菌属的丰度大于5.62×10-5，则该个体为藏族，当个体的粪便或肠道内容物中的琥珀酸弧菌属的丰度小于5.62×10-5，则该个体为汉族。所述的步骤(1)中包括：1)收集待检测个体的粪便或肠道内容物，并提取粪便或肠道内容物中微生物的dna；2)目的片段的扩增；3)测序；4)计算琥珀酸弧菌属的丰度。所述的步骤1)中粪便或肠道内容物中微生物的dna的提取方法选自：ctab法、gitc法或利用商业化的试剂盒进行。优选的，利用ctab法提取粪便或肠道内容物中微生物的dna。具体的，将收集到的粪便或肠道内容物用ctab、溶菌酶和裂解液处理后，用苯酚、氯仿和异戊醇处理使dna游离及去除杂蛋白、利用异丙醇使dna析出，用乙醇纯化dna后得到粪便或肠道内容物中微生物的dna。更具体的，收集待检测个体的粪便，吸取500-1000ulctab裂解液至ep管里，加入溶菌酶，将200-500ul左右的样品加入裂解液中，50-70℃水浴条件下反应，离心取上清，加苯酚(ph8.0)：氯仿：异戊醇(20-30：18-28：1)，颠倒混匀，10000-12000rpm离心10-15min。取上清，加氯仿：异戊醇(20-30:1)，颠倒混匀，10000-12000rpm离心10-15min。取上清，加入异丙醇，上下摇晃，-20℃沉淀。10000-12000rpm离心10-15min，将沉淀用乙醇洗涤2-5次，将沉淀晾干，加入水溶解沉淀，加rnasea消化rna后既得粪便或肠道内容物中微生物的dna。所述的步骤2)中使用带barcode的特异引物515f和806r对步骤(1)中获得的dna进行扩增，所述的515f的序列为5’-gtgccagcmgccgcggtaa-3’(seqidno：1)，所述的806r的序列为5’-ggactachvgggtwtctaat-3’(seqidno：2)。所述的步骤3)中是利用高通量测序法对步骤(2)中获得的产物进行测序。所述的步骤4)中将测序得到的数据聚类成为otus(operationaltaxonomicunits)，对物种注释分析后，计算琥珀酸弧菌属的丰度。具体的，所述的步骤4)中将步骤3)中测序得到的数据进行处理，利用uparse软件(uparsev7.0.1001，http://www.drive5.com/uparse/)对有效数据进行聚类成为otus(operationaltaxonomicunits)，用mothur方法与silva(http://www.arb-silva.de/)的ssurrna数据库对otus序列进行物种注释分析，以(琥珀酸弧菌属otu值)/(检测样本检测到的所有菌属otu值)表示琥珀酸弧菌属相对丰度。所述的琥珀酸弧菌属包括溶糊精琥珀酸弧菌(succinivibriodextrinosolvens)。第三方面，本发明提供检测微生物标志物丰度的产品在鉴别和/或区分不同种族个体中的应用。所述的检测微生物标志物丰度的产品包括用于检测微生物标志物丰度的试剂、仪器或装置。优选的，所述的检测微生物标志物丰度的产品为检测琥珀酸弧菌属丰度的产品。所述的琥珀酸弧菌属包括溶糊精琥珀酸弧菌(succinivibriodextrinosolvens)。优选的，所述的鉴别和/或区分不同种族个体是指鉴别和/或区分个体为藏族或汉族。第四方面，本发明提供一种用于鉴别和/或区分不同种族个体的产品，所述的产品中含有用于检测微生物标志物丰度的试剂、仪器或装置。优选的，所述的鉴别和/或区分不同种族个体的产品为鉴别和/或区分个体为藏族或汉族的产品。优选的，所述的产品中含有检测琥珀酸弧菌属丰度的试剂、仪器或装置。所述的琥珀酸弧菌属包括溶糊精琥珀酸弧菌(succinivibriodextrinosolvens)。本发明中所述的微生物丰度指在某一微生物群体中该种微生物的丰度程度，例如，在肠道微生物群体中该种微生物的丰富程度，可表示为该种微生物在该群体中的含量。附图说明图1所示为箱式图显示的汉藏族肠道群中琥珀酸弧菌属的相对丰度的差异；图2所示为roc曲线显示琥珀酸弧菌属相对丰度鉴别诊断高原人群汉藏来源的敏感性和特异性。图3为藏族士兵组(zang)粪便中群落组成。图4汉族士兵(han)组粪便中群落组成。具体实施方式实施例1利用个体粪便中琥珀酸弧菌属的相对丰度鉴别受试者为藏族或汉族(1)确定所述的个体的粪便样本中的所述的微生物标志物的丰度：1)收集待检测个体的粪便，并提取粪便中微生物的dna：收集待检测个体的粪便，吸取1000ulctab裂解液至2.0mlep管里，加入溶菌酶，将500ul左右的样品加入裂解液中，65℃水浴，期间颠倒混匀数次，以使样品充分裂解。离心取上清，加苯酚(ph8.0)：氯仿：异戊醇(25：24：1)，颠倒混匀，12000rpm离心10min。取上清，加氯仿：异戊醇(24:1)，颠倒混匀，12000rpm离心10min。吸取上清至1.5ml离心管里，加入异丙醇，上下摇晃，-20℃沉淀。12000rpm离心10分钟，倒出液体，注意不要倒出沉淀。用1ml75％乙醇洗涤2次，剩余的少量液体可再次离心收集，然后用枪头吸出。超净工作台吹干或者室温晾干。加入ddh2o溶解dna样品，加rnasea1ul消化rna，37℃放置15min。之后利用琼脂糖凝胶电泳检测dna的纯度和浓度，取适量的样本dna于离心管中，使用无菌水稀释样本至1ng/μl。2)目的片段的扩增：以稀释后的基因组dna为模板，根据测序区域16sv3-v4，选择使用带barcode的特异引物：515f：5’-gtgccagcmgccgcggtaa-3’(seqidno：1)，806r：5’-ggactachvgggtwtctaat-3’(seqidno：2)，newenglandbiolabs公司的high-fidelitypcrmastermixwithgcbuffer，和高效高保真酶进行pcr，确保扩增效率和准确性。pcr产物使用2％浓度的琼脂糖凝胶进行电泳检测；根据pcr产物浓度进行等量混样,充分混匀后使用1×tae浓度2％的琼脂糖胶电泳纯化pcr产物，剪切回收目标条带。产物纯化试剂盒使用的是thermoscientific公司genejet胶回收试剂盒回收产物。3)测序：使用thermofisher公司的ionplusfragmentlibrarykit48rxns建库试剂盒进行文库的构建，构建好的文库经过qubit定量和文库检测合格后，使用thermofisher的ions5tmxl进行上机测序。4)计算琥珀酸弧菌属的丰度：①测序数据处理：使用cutadapt(v1.9.1,http://cutadapt.readthedocs.io/en/stable/)先对reads进行低质量部分剪切，再根据barcode从得到的reads中拆分出各样品数据，截去barcode和引物序列初步质控得到原始数据，经过以上处理后得到的reads需要进行去除嵌合体序列的处理，reads序列通过与物种注释数据库进行比对检测嵌合体序列，并最终去除其中的嵌合体序列，得到最终的有效数据。②otu聚类和物种注释：利用uparse软件(uparsev7.0.1001，http://www.drive5.com/uparse/)对所有样品的全部有效数据进行聚类，默认以97％的一致性将序列聚类成为otus(operationaltaxonomicunits)，同时会选取otus的代表性序列，依据其算法原则，筛选的是otus中出现频数最高的序列作为otus的代表序列。对otus序列进行物种注释，用mothur方法与silva(http://www.arb-silva.de/)的ssurrna数据库进行物种注释分析(设定阈值为0.8～1)，获得分类学信息并分别在属水平统计汉族士兵组、藏族士兵组的群落组成(见图3、4)。琥珀酸弧菌属相对丰度用(琥珀酸弧菌属otu值)/(检测样本检测到的所有菌属otu值)表示。(2)将待检测个体粪便中琥珀酸弧菌属的丰度与5.62×10-5比较，当个体的粪便或肠道内容物中的琥珀酸弧菌属的丰度大于5.62×10-5，则该个体为藏族，当个体的粪便或肠道内容物中的琥珀酸弧菌属的丰度小于5.62×10-5，则该个体为汉族。实施例2利用琥珀酸弧菌属的相对丰度鉴别受试者的种族的可行性鉴定研究人群：受试者为中国服役男性士兵，实验组为128个习服海拔3500米高原以上的藏族男性士兵，对照组为与藏族士兵年龄匹配128个平原汉族士兵。严格控制两组人群的进行相同的中国部队混合饮食，保持相同的海拔3500米训练环境和训练强度，戒烟、酒，持续3个月。排除有慢性炎性疾病、口服抗生素、急性感染和胃肠道疾病的士兵。方法：(1)确定所述的个体的粪便样本中的所述的微生物标志物的丰度：1)收集待检测个体的粪便，并提取粪便中微生物的dna：收集待检测个体的粪便，吸取1000ulctab裂解液至2.0mlep管里，加入溶菌酶，将500ul左右的样品加入裂解液中，65℃水浴，期间颠倒混匀数次，以使样品充分裂解。离心取上清，加苯酚(ph8.0)：氯仿：异戊醇(25：24：1)，颠倒混匀，12000rpm离心10min。取上清，加氯仿：异戊醇(24:1)，颠倒混匀，12000rpm离心10min。吸取上清至1.5ml离心管里，加入异丙醇，上下摇晃，-20℃沉淀。12000rpm离心10分钟，倒出液体，注意不要倒出沉淀。用1ml75％乙醇洗涤2次，剩余的少量液体可再次离心收集，然后用枪头吸出。超净工作台吹干或者室温晾干。加入ddh2o溶解dna样品，加rnasea1ul消化rna，37℃放置15min。之后利用琼脂糖凝胶电泳检测dna的纯度和浓度，取适量的样本dna于离心管中，使用无菌水稀释样本至1ng/μl。2)目的片段的扩增：以稀释后的基因组dna为模板，根据测序区域16sv3-v4，选择使用带barcode的特异引物：515f：5’-gtgccagcmgccgcggtaa-3’(seqidno：1)，806r：5’-ggactachvgggtwtctaat-3’(seqidno：2)，newenglandbiolabs公司的high-fidelitypcrmastermixwithgcbuffer，和高效高保真酶(taqdnapolymerase)进行pcr，确保扩增效率和准确性。pcr产物使用2％浓度的琼脂糖凝胶进行电泳检测；根据pcr产物浓度进行等量混样,充分混匀后使用1×tae浓度2％的琼脂糖胶电泳纯化pcr产物，剪切回收目标条带。产物纯化试剂盒使用的是thermoscientific公司genejet胶回收试剂盒回收产物。3)测序：使用thermofisher公司的ionplusfragmentlibrarykit48rxns建库试剂盒进行文库的构建，构建好的文库经过qubit定量和文库检测合格后，使用thermofisher的ions5tmxl进行上机测序。4)计算琥珀酸弧菌属的丰度：①测序数据处理：使用cutadapt(v1.9.1,http://cutadapt.readthedocs.io/en/stable/)先对reads进行低质量部分剪切，再根据barcode从得到的reads中拆分出各样品数据，截去barcode和引物序列初步质控得到原始数据，经过以上处理后得到的reads需要进行去除嵌合体序列的处理，reads序列通过与物种注释数据库进行比对检测嵌合体序列，并最终去除其中的嵌合体序列，得到最终的有效数据。②otu聚类和物种注释：利用uparse软件(uparsev7.0.1001，http://www.drive5.com/uparse/)对所有样品的全部有效数据进行聚类，默认以97％的一致性将序列聚类成为otus(operationaltaxonomicunits)，同时会选取otus的代表性序列，依据其算法原则，筛选的是otus中出现频数最高的序列作为otus的代表序列。对otus序列进行物种注释，用mothur方法与silva(http://www.arb-silva.de/)的ssurrna数据库进行物种注释分析(设定阈值为0.8～1)，获得分类学信息并分别在属水平统计汉族士兵组、藏族士兵组的群落组成(见图3、4)。琥珀酸弧菌属相对丰度用(琥珀酸弧菌属otu值)/(检测样本检测到的所有菌属otu值)表示。结果：图1所示为汉藏族肠道群中琥珀酸弧菌属的相对丰度的差异的箱式图，由图1所示，高原汉族人群中琥珀酸弧菌属相对丰度为4.17×10-5，高原藏族人群中琥珀酸弧菌属的相对丰度为1.49×10-4，藏族人群琥珀酸弧菌属相对丰度显著高于汉族人群。通过上述验证试验说明可以利用本发明提供的利用个体粪便中的琥珀酸弧菌属相对丰度鉴别个体的种族。实施例3构建roc曲线比较评估菌属相对丰度鉴别高原汉族人群和藏族人群的能力采用受试者工作曲线(roc)法进行验证，通过受试者粪便中的琥珀酸弧菌属相对丰度判断其用于鉴别高原汉族人群和藏族人群的能力。差异比较因为样本非正态分布采用mann-whitney检验。绘制roc曲线评估菌属相对丰度对高原汉族人群和藏族人群的鉴别诊断敏感性和特异性，计算界值，p＜0.05为差异有统计学意义，采用spss22.0软件包进行统计学处理。由图2的结果表明，roc曲线显示琥珀酸弧菌属相对丰度鉴别诊断高原人群汉藏来源的auc(roc曲线下面积)值为0.933，说明琥珀酸弧菌属相对丰度判断其用于鉴别高原汉族人群和藏族人群，在菌属相对丰度为5.62×10-5时，可得到最大的约登指数0.863，敏感度和特异性分别为0.873和0.956，正确率为0.914，f1分数为0.912，说明利用16srrna测序检测粪便标本琥珀酸弧菌属相对丰度，可以鉴别高原人群来源是藏族还是汉族，当琥珀酸弧菌属相对丰度大于5.62×10-5可判定为藏族。序列表<110>中国人民解放军总医院<120>琥珀酸弧菌属在鉴别和/或区分不同种族个体中的应用<130>1<160>2<170>patentinversion3.5<210>1<211>19<212>dna/rna<213>人工序列(artificialsequence)<400>1gtgccagcmgccgcggtaa19<210>2<211>20<212>dna/rna<213>人工序列(artificialsequence)<400>2ggactachvgggtwtctaat20当前第1页12当前第1页12