小柳-原田综合症标志微生物及其应用的制作方法

1.本发明涉及生物技术领域，具体的，本发明涉及小柳-原田综合症标志微生物及其用途，更具体的，本发明涉及一种试剂盒、试剂在制备试剂盒中的用途、用于预防或者治疗小柳-原田综合症的药物组合物或者食品组合物、确定个体是否患有小柳-原田综合症方法、确定个体是否患有小柳-原田综合症的装置、一种装置、一种筛选药物的方法。


背景技术：

2.小柳-原田综合症(vogt-koyanagi-harada，vkh)又称葡萄膜大脑炎(uveoencephalitis)综合征，也称色素膜脑膜炎或眼-脑-耳-皮综合征。是一种自身免疫反应引起的攻击黑色素细胞抗原的致盲眼病，往往还会侵袭耳、脑(脊)膜、毛发以及皮肤等器官。起病多发于青壮年，发病率与性别无明显的关系，黄种人种较为常见。小柳-原田综合症的病因前尚未完全明确，文献报告有内分泌学说、病毒学说和免疫学说等三种学说。
3.小柳-原田综合症一般起病比较隐匿，临床上一般分为三期，分别是前驱期、眼病期、恢复期。前驱期有感冒样的症状，如发热、头痛、耳鸣、听力异常、颈项僵硬、头皮过敏等。由于病情较轻，病人大多会误，致使病情延误，失去最佳治疗时机。
4.小柳-原田综合症的诊断标准定为：
①
病前无眼球外伤或手术史；
②
临床与实验室排除了其他病因的葡萄膜炎；
③
双眼累及具有小柳-原田综合症的急性或慢性期，眼前节或后节的表现；
④
有神经、耳部症状或脑脊液细胞增多；
⑤
皮肤、毛发改变。虽然临床表现加上器械检查可以对病做出正确诊断，但是大部分需要患者还要通过病史、眼科常规检查、造影、光学成像等作出诊断。目前，小柳-原田综合症的早期诊断主要是用血常规、超声显微镜、眼底血管造影、x射线检查等技术进行检查，检查过程较复杂。
5.随着人体基因组测序完成及高通量测序技术的高速发展，基因筛查成为小柳-原田综合症诊断的方向，对于发现小柳-原田综合症潜在人群很有优势。有研究表明与小柳-原田综合症相关的多基因在肠道免疫中起重要作用，如tlr受体能够识别病原微生物配体，随即引发自身免疫反应响应。途径的基因对肠道健康的调节起重要作用(kawai t et al.,2010)。到目前为止，还没有小柳-原田综合症患者肠道微生物标记物的相关报道。因此，研究小柳-原田综合症患者肠道微生物群的特征对小柳-原田综合症的检测及治疗具有重要意义。


技术实现要素：

6.本技术是基于发明人对以下事实和问题的发现和认识作出的：
7.本技术的申请人通过前期大量研究，意外地发现一些微生物可以作为检测小柳-原田综合症的标志微生物，为早期发现小柳-原田综合症提供一种非侵入性方法；合理有效地应用标志微生物，扶持肠道有益菌的生长，抑制肠道潜在致病菌，可以治疗或减轻小柳-原田综合症的临床症状。
8.为此，在本发明的第一方面，本发明提出了一种试剂盒。根据本发明的实施例，所
述试剂盒包括适于检测第一微生物集中的至少一种菌种的试剂，所述第一微生物集由以下菌种组成：固氮螺菌(azospirillum_sp.)_cag:260和梭菌(clostridium_sp.)_cag:349。根据本发明实施例的试剂盒可以准确的区分小柳-原田综合症患者和健康个体。
9.在本发明的第二方面，本发明提出了试剂在制备试剂盒中的用途，所述试剂适于检测第一微生物集中的至少一种菌种。根据本发明的实施例，所述试剂盒用于诊断小柳-原田综合症或检测小柳-原田综合症的治疗效果，所述第一微生物集由以下菌种组成：固氮螺菌(azospirillum_sp.)_cag:260和梭菌(clostridium_sp.)_cag:349。根据本发明具体实施例的试剂制备的试剂盒，可以准确的检测所述第一微生物集中的至少一种菌种，极准确的区分小柳-原田综合症患者和健康个体，由此，可以有效的在早期进行小柳-原田综合症诊断，或用于检测治疗过程中小柳-原田综合症的变化。
10.在本发明的第三方面，本发明提出了一种用于预防或者治疗小柳-原田综合症的药物组合物或者食品组合物。根据本发明的实施例，含有第一微生物集中的至少一种菌种，所述第一微生物集由以下菌种组成：固氮螺菌(azospirillum_sp.)_cag:260和梭菌(clostridium_sp.)_cag:349。根据本发明实施例的小柳-原田综合症标志微生物中第一微生物集的菌种可以非侵入性的在早期发现或辅助检测小柳-原田综合症，确定个体患有小柳-原田综合症的概率高低或者个体处于健康状态的概率高低；同时，提高小柳-原田综合症高风险人群或小柳-原田综合症患者肠道内的所述第一微生物集中的各种菌种，可以降低患小柳-原田综合症的概率或减缓、治愈小柳-原田综合症，因此，所述包含所述第一微生物集中的至少一种菌种的药物或者食品组合物能够用于平衡肠道菌群，有效预防或治疗小柳-原田综合症。
11.在本发明的第四方面，本发明提出了一种确定个体是否患有小柳-原田综合症的方法。根据本发明的实施例，包括：(1)确定所述个体的粪便样本中标志微生物的丰度，所述标志微生物包括第一微生物集和第二微生物集中的至少一种菌种；(2)将步骤(1)中得到的所述丰度与预定的阈值进行比较，以便确定所述个体是否患有小柳-原田综合症；其中，所述第一微生物集由以下菌种组成：固氮螺菌(azospirillum_sp.)_cag:260和梭菌(clostridium_sp.)_cag:349；所述第二微生物集由以下菌种组成：副普氏菌(paraprevotella_clara)、厚壁菌(firmicutes_bacterium)_cag:194和酪黄肠球菌(enterococcus_casseliflavus)。根据本发明实施例的方法可以依据个体的粪便样本中的所述标志微生物的丰度确定个体是否患有小柳-原田综合症，所述标志微生物是发明人对大量已知状态的粪便样本进行验证，通过差异比较分析各种肠道微生物在小柳-原田综合症组和健康组粪便样本中的丰度，而确定下来的。
12.在本发明的第五方面，本发明提出了一种确定个体是否患有小柳-原田综合症的装置。根据本发明的实施例，包括：丰度确定单元，用于确定所述个体的粪便样本中标志微生物的丰度，所述标志微生物包括第一微生物集和第二微生物集中的至少一种菌种；比较单元，用于将所得到的所述丰度与预定的阈值进行比较，以便确定所述个体是否患有小柳-原田综合症；其中，所述第一微生物集由以下菌种组成：固氮螺菌(azospirillum_sp.)_cag:260和梭菌(clostridium_sp.)_cag:349；所述第二微生物集由以下菌种组成：副普氏菌(paraprevotella_clara)、厚壁菌(firmicutes_bacterium)_cag:194和酪黄肠球菌(enterococcus_casseliflavus)。所述标志微生物是发明人通过差异比较分析各种肠道微
生物在小柳-原田综合症患者和健康人群的粪便样本中的丰度并经过分析和大量已知状态的粪便样本的验证而确定下来的，根据本发明实施例的装置可以准确确定个体是否为小柳-原田综合症的高风险人群或小柳-原田综合症患者。
13.在本发明的第六方面，本发明提出了一种装置。根据本发明的实施例，包括：计算机可读存储介质，其上存储有计算机程序，所述程序用于执行第四方面所述的方法；以及一个或者多个处理器，用于执行所述计算机可读存储介质中的程序。根据本发明实施例的装置可以准确确定个体是否为小柳-原田综合症的高风险人群或小柳-原田综合症患者。
14.在本发明的第七方面，本发明提出了一种筛选药物的方法。根据本发明的实施例，所述药物用于治疗或者预防小柳-原田综合症，所述方法包括：将候选药物施用于受试者，检测施用前后，所述受试者粪便中标志微生物的丰度，所述标志微生物包括第一微生物集和第二微生物集中的至少一种菌种，其中，满足下列条件至少之一的候选药物适于用于治疗或者预防小柳-原田综合症：(1)进行所述施用后，所述第一微生物集中的至少一种菌种的所述丰度升高；和(2)进行所述施用后，所述第二微生物集中的至少一种菌种的所述丰度降低；其中，所述第一微生物集由以下菌种组成：固氮螺菌(azospirillum_sp.)_cag:260和梭菌(clostridium_sp.)_cag:349；所述第二微生物集由以下菌种组成：副普氏菌(paraprevotella_clara)、厚壁菌(firmicutes_bacterium)_cag:194和酪黄肠球菌(enterococcus_casseliflavus)。根据本发明实施例的方法可以生产或筛选出促进所述标志微生物中第一微生物集中各种菌种生长和/或抑制肠道标志微生物中第二微生物集中的各种菌种生长的药物，对于辅助减轻小柳-原田综合症的临床症状具有重要意义。
15.本发明的附加方面和优点将在下面的描述中部分给出，部分将从下面的描述中变得明显，或通过本发明的实践了解到。
附图说明
16.本发明的上述和/或附加的方面和优点从结合下面附图对实施例的描述中将变得明显和容易理解，其中：
17.图1是根据本发明实施例的筛选小柳-原田综合症标志微生物的实验分析流程示意图；以及
18.图2是根据本发明实施例的标志微生物综合指标auc评价结果示意图，其中，specificity表示特异度，即预测为阳性且实际为阳性，真阳性，纵坐标sensitivity表示敏感度，即真阴性，confidence interval表示置信区间：
19.2-a为第一期74个样品数据roc曲线下auc值和置信区间结果图；
20.2-b为第二期33个样品数据roc曲线下auc值和置信区间结果图。
具体实施方式
21.下面详细描述本发明的实施例，所述实施例的示例在附图中示出。下面通过参考附图描述的实施例是示例性的，旨在用于解释本发明，而不能理解为对本发明的限制。
22.此外，术语“第一”、“第二”仅用于描述目的，而不能理解为指示或暗示相对重要性或者隐含指明所指示的技术特征的数量。由此，限定有“第一”、“第二”的特征可以明示或者隐含地包括至少一个该特征。在本发明的描述中，“多个”的含义是至少两个，例如两个，三
个等，除非另有明确具体的限定。
23.术语“任选地”仅用于描述目的，而不能理解为指示或暗示相对重要性。由此，限定有“任选地”的特征可以明示或者隐含地包括或不包括该特征。
24.生物学标志物是从生物学介质中可以检测到的细胞/生物化学或分子改变。生物学介质包括各种体液、组织、细胞、粪便、头发、呼气等。
25.微生物的丰度是指在某一微生物群体中该种微生物的丰富程度，例如在肠道微生物群体中的该种微生物程度，可表示为该种微生物在该群体中的含量。
26.根据本发明的一个实施方式提供的一种试剂盒，包括适于检测第一微生物集中的至少一种菌种的试剂，所述第一微生物集由以下菌种组成：固氮螺菌(azospirillum_sp.)_cag:260和梭菌(clostridium_sp.)_cag:349。
27.根据本发明的具体实施方案，所述试剂盒进一步包括适于检测第二微生物集中的至少一种菌种的试剂，所述第二微生物集由以下菌种组成：副普氏菌(paraprevotella_clara)、厚壁菌(firmicutes_bacterium)_cag:194和酪黄肠球菌(enterococcus_casseliflavus)。
28.根据本发明的具体实施方案，所述试剂盒包括适于检测所述第一微生物集中全部所述菌种的试剂。
29.根据本发明的具体实施方案，所述试剂盒包括适于检测所述第二微生物集中全部所述菌种的试剂。
30.根据本发明的具体实施方案，所述标志微生物是发明人通过对大量患小柳-原田综合症个体和大量健康对照个体的粪便样本中的微生物的丰度的差异比较分析、以及验证，而确定下来的，明确了肠道微生物中小柳-原田综合症相关的标志微生物。利用包含检测所述标志微生物的试剂的试剂盒能够确定个体处于患有小柳-原田综合症状态的概率高低或者处于健康状态的概率高低，能够用于非侵入性的早期发现或辅助检测小柳-原田综合症。
31.根据本发明的具体实施方案，所述适于检测所述第一微生物集或第二微生物集的试剂不受特别限制，任何可以检测所述微生物菌种的试剂均包含在本发明的范围内，如通过形态学特征、生理生化反应特征、生态学特征以及血清学反应、对噬菌体的敏感性、分子生物学等方面检测所述微生物菌种的试剂，具体的，如抗体、酶、核酸分子等。
32.本文中，所述微生物形态学特征指：在显微镜下观察微生物外形大小、形状、排列等，细胞构造，革兰氏染色反应，能否运动、鞭毛着生部位和数目，有无芽孢和荚膜、芽孢的大小和位置，放线菌和真菌的繁殖器官的形状、构造，孢子的数目、形状、大小、颜色和表面特征等。
33.本文中，所述微生物生理生化反应特征指：微生物利用物质的能力、代谢产物的特殊性，如是否产生h2s、吲哚、co2、醇、有机酸，能否还原硝酸盐，能否使牛奶凝固、冻化等、生长环境(温度、湿度、氧气及二氧化碳等气体的浓度、ph、是否耐高渗、是否有嗜盐性等)、与其它生物间的关系(如：共生、寄生、宿主范围及致病情况)等。
34.本文中，所述微生物血清学反应指：利用抗原与抗体的高度敏感特异性反应，来鉴别相似的菌种，或对同种微生物分型，如用已知菌种、型或菌株制成的抗血清，与待鉴定的所述微生物是否发生特异性的血清学反应来鉴定微生物。
35.本文中，分子生物学方法检测微生物主要包括：利用pcr技术、高通量测序等方法。
36.根据本发明提供的试剂在制备试剂盒中的用途，所述试剂适于检测第一微生物集中的至少一种菌种，所述试剂盒用于诊断小柳-原田综合症或者检测小柳-原田综合症的治疗效果，所述第一微生物集由以下菌种组成：固氮螺菌(azospirillum_sp.)_cag:260和梭菌(clostridium_sp.)_cag:349。
37.根据本发明的具体实施例，所述标志微生物是发明人通过对大量患小柳-原田综合症个体和大量健康对照个体的粪便样本中的微生物的丰度的差异比较分析、以及验证，而确定下来的，明确了肠道微生物中小柳-原田综合症相关的微生物标志物。利用检测所述标志微生物的试剂能够确定个体患有小柳-原田综合症的概率高低或者处于健康状态的概率高低，能够用于非侵入性的早期发现或辅助检测小柳-原田综合症。
38.根据本发明一些具体的实施例，所述试剂进一步适于检测第二微生物集中的至少一种菌种，所述第二微生物集由以下菌种组成：副普氏菌(paraprevotella_clara)、厚壁菌(firmicutes_bacterium)_cag:194和酪黄肠球菌(enterococcus_casseliflavus)。
39.根据本发明的具体实施例，所述适于检测所述第一微生物集或第二微生物集的试剂不受特别限制，可以检测所述微生物菌种的试剂均包含在本发明的范围内，如通过形态学特征、生理生化反应特征、生态学特征以及血清学反应、对噬菌体的敏感性、分子生物学等方面检测所述微生物菌种的试剂等，具体的，如抗体、酶、核酸分子。
40.本文中，所述微生物形态学特征指：在显微镜下观察微生物外形大小、形状、排列等，细胞构造，革兰氏染色反应，能否运动、鞭毛着生部位和数目，有无芽孢和荚膜、芽孢的大小和位置，放线菌和真菌的繁殖器官的形状、构造，孢子的数目、形状、大小、颜色和表面特征等。
41.本文中，所述微生物生理生化反应特征指：微生物利用物质的能力、代谢产物的特殊性，如是否产生h2s、吲哚、co2、醇、有机酸，能否还原硝酸盐，能否使牛奶凝固、冻化等、生长环境(适宜生长的温度、湿度、氧气及二氧化碳等气体的浓度、ph、是否耐高渗、是否有嗜盐性等)、与其它生物间的关系(如：共生、寄生、宿主范围及致病情况)等。
42.本文中，所述微生物血清学反应指：利用抗原与抗体的高度敏感特异性反应，来鉴别相似的菌种，或对同种微生物分型，如用已知菌种、型或菌株制成的抗血清，与待鉴定的所述微生物是否发生特异性的血清学反应来鉴定微生物。
43.本文中，分子生物学方法检测微生物主要包括：利用pcr技术、高通量测序等方法。
44.根据本发明提供的一种用于预防或者治疗小柳-原田综合症的药物组合物或者食品组合物，含有第一微生物集中的至少一种菌种，所述第一微生物集由以下菌种组成：固氮螺菌(azospirillum_sp.)_cag:260和梭菌(clostridium_sp.)_cag:349。
45.上述标志微生物是发明人通过差异比较分析各种肠道微生物在小柳-原田综合症疾病组和健康组的粪便样本中的丰度，以及经过大量已知状态的粪便样本的验证，而确定下来的。所述标志微生物中的菌种相较于小柳-原田综合症患者群体，在健康群体组中显著富集，所述显著富集是指与在小柳-原田综合症患者组中的丰度相比，上述菌种在健康组中的丰度均具有统计意义地高于或者明显地、实质性地高于在小柳-原田综合症患者组中的丰度；能够使该部分菌种丰度提高的物质能够用于治疗小柳-原田综合症或者益于小柳-原田综合症患者服用，能够使其丰度提高的物质不限于治疗小柳-原田综合症的药物和有益
肠道菌群平衡的功能性食品。该实施方式提供的标志微生物能够用于制备治疗小柳-原田综合症的药物和/或用于制备益于平衡肠道菌群的功能性食品、保健药等，所述药物或食品可有效治疗或缓解小柳-原田综合症。
46.根据本发明提供的一种确定个体是否患有小柳-原田综合症的方法，包括步骤(1)和(2)。
47.(1)确定所述个体的粪便样本中的标志微生物的丰度。
48.所述标志微生物包括第一微生物集和第二微生物集中的至少一种菌种。其中，所述第一微生物集由以下菌种组成：固氮螺菌(azospirillum_sp.)_cag:260和梭菌(clostridium_sp.)_cag:349；所述第二微生物集由以下菌种组成：副普氏菌(paraprevotella_clara)、厚壁菌(firmicutes_bacterium)_cag:194和酪黄肠球菌(enterococcus_casseliflavus)。
49.根据本发明的一些具体的实施例，步骤(1)进一步包括：获得所述个体的粪便样本中的核酸测序数据；将所述测序数据与参考基因组进行比对；基于所述比对的结果，确定所述标志微生物的丰度。
50.根据本发明的具体实施例，在步骤(1)中，按照下列公式确定所述标志微生物的丰度：
51.ab(s)＝ab(us) ab(ms)，其中，s表示所述标志微生物的编号，ab(s)表示所述标志微生物s的丰度，ab(us)＝us/ls，us为所述测序数据中与所述标志微生物s的参考基因组唯一比对的读段数目，ls为所述标志微生物s的参考基因组的总长度，ms为所述测序数据中与所述标志微生物s的参考基因组非唯一比对的读段的数目，i表示所述非唯一比对读段的编号，coi为所述第i读段对应的丰度系数，co
i,s
表示针对所述标志微生物s，所述非唯一比对的读段i的丰度系数，n为所述非唯一比对的读段i能够比对的微生物的总数，j表示所述非唯一比对的读段i能够比对的微生物的编号。
52.比对可以利用已知比对软件进行，例如soap、bwa和teramap等，在比对过程中，一般对比对参数进行设置，设置一个或者一对读段(reads)最多允许有s个碱基错配(mismatch)，例如设置s≤2，若reads中有超过s个碱基发生错配，则视为该reads无法比对到(比对上)该组装片段上。所述的获得的比对结果包含各条读段与各物种的参考基因组的比对情况，包括读段是否能够比对上某个或某些物种的参考基因组、只唯一比对到一种物种还是比对到多种物种的参考基因组、比对到物种的参考基因组的位置、比对到物种参考基因组的唯一位置还是多个位置等信息。
53.所述菌种/微生物的参考基因组指预先确定的该微生物物种的序列，可以是预先获得的待测样本所属或者所包含的生物类别的任意参考模板，例如，目标是待测样本中的微生物，参考序列可选择ncbi数据库中的各种微生物的参考基因组和/或hmp、metahit项目公开的dacc肠道参考基因组，进一步地，也可以预先配置包含更多参考序列的资源库，例如
依据待测样本来源的个体的状态、地域等因素选择或是测定组装出更接近的序列作为参考序列。根据本发明的一个实施例，各种微生物的参考基因组从公开的数据库中获得，通常的，一种微生物的参考基因组有多个版本，即一种微生物有多个公开的参考基因组。
54.reads与物种的参考基因组比对，比对上的可以被分为两部分：a)unique reads(u)：唯一比对上一个物种的参考基因组；称这些reads为unique reads。即，如果reads比对上的参考基因组均来自同一物种，定义这些reads为unique reads；b)multiple reads(m)：比对上一个以上物种的参考基因组，定义为multiple reads。即，如果reads比对上的参考基因组来自至少两种物种，定义这些reads为multiple reads。
55.(2)丰度比较，以确定个体是否患有小柳-原田综合症。
56.根据本发明的一个实施例，将步骤(1)中得到的所述丰度与预定的阈值进行比较，以便确定所述个体是否患有小柳-原田综合症。
57.根据本发明的一些具体实施例，所述阈值为预先设定的。将标志微生物在待测个体样本中的丰度与所述的阈值比较，进行待测个体的状态的确定。所述阈值可以为一数值或者数值范围，具体的设定方式不受任何限制，例如，基于已知患病或健康状态个体中的标志微生物的丰度均值，该微生物对应的阈值可以设为该丰度均值的95％的置信区间(confidence interval)。
58.所述的置信区间是指由样本统计量所构造的总体参数的估计区间。在统计学中，一个概率样本的置信区间是对这个样本的某个总体参数的区间估计。置信区间展现的是这个参数的真实值有一定概率落在测量结果的周围的程度。置信区间给出的是被测量参数的测量值的可信程度，即前面所要求的“一定概率”，这个概率被称为置信水平。
59.根据本发明的一些具体实施例，当步骤(1)中确定的标志微生物的丰度达到所述患小柳-原田综合症丰度阈值，未达到所述不患小柳-原田综合症丰度阈值时，确定所述个体患有小柳-原田综合症，当(1)中确定的标志微生物的丰度达到不患小柳-原田综合症丰度阈值，未达到患小柳-原田综合症丰度阈值时，确定所述个体不患小柳-原田综合症。
60.需要说明的是，根据目的或要求不同，可能对确定个体状态结果的可信程度有不同的要求，本领域技术人员可以选择不同的显著性水平或阈值。
61.该方法基于检测个体的粪便样本中的标志微生物中的各种菌种的丰度，分别将检测确定的各种菌种的丰度与其阈值进行比较，依据获得的比较结果能够确定个体为小柳-原田综合症个体或者为健康个体的概率。为早期发现小柳-原田综合症提供一种非侵入性的辅助检测或者辅助干预治疗的方法。
62.以上任一实施例中的利用标志微生物确定个体是否患有小柳-原田综合症的方法的全部或部分步骤，可以利用包含可拆分的相应单元功能模块的装置/系统来施行，或者将方法程序化、存储于机器可读介质，利用机器运行该可读介质来实现。
63.根据本发明提供的一种确定个体是否患有小柳-原田综合症的装置，该装置包括：丰度确定单元，用于确定所述个体的粪便样本中标志微生物的丰度；比较单元，用于将所得到的所述丰度与预定的阈值进行比较，以便确定所述个体是否患有小柳-原田综合症；其中，所述第一微生物集由以下菌种组成：固氮螺菌(azospirillum_sp.)_cag:260和梭菌(clostridium_sp.)_cag:349；所述第二微生物集由以下菌种组成：副普氏菌(paraprevotella_clara)、厚壁菌(firmicutes_bacterium)_cag:194和酪黄肠球菌
(enterococcus_casseliflavus)。上述对本发明任一实施例的利用标志微生物确定个体是否患有小柳-原田综合症的方法的技术特征和优点的描述，同样适用本发明这一方面的装置，在此不再赘述。
64.根据本发明的实施方案，所述丰度确定单元适于通过下列步骤确定所述丰度：获得所述个体的粪便样本中的核酸测序数据；将所述测序数据与参考基因组进行比对；基于所述比对的结果，确定所述标志微生物的丰度。
65.比对可以利用已知比对软件进行，例如soap、bwa和teramap等，在比对过程中，一般对比对参数进行设置，设置一个或者一对读段(reads)最多允许有s个碱基错配(mismatch)，例如设置s≤2，若reads中有超过s个碱基发生错配，则视为该reads无法比对到(比对上)该组装片段上。所述的获得的比对结果包含各条读段与各物种的参考基因组的比对情况，包括读段是否能够比对上某个或某些物种的参考基因组、只唯一比对到一种物种还是比对到多种物种的参考基因组、比对到物种的参考基因组的位置、比对到物种参考基因组的唯一位置还是多个位置等信息。根据本发明的一个实施例，利用soapalign 2.21进行比对，设置参数为
–
r 2
–
m 100
–
x 1000。
66.所述微生物的参考基因组指预先确定的该微生物物种的序列，可以是预先获得的待测样本所属或者所包含的生物类别的任意参考模板，例如，目标是待测样本中的微生物，参考序列可选择ncbi数据库中的各种微生物的参考基因组和/或hmp、metahit项目公开的dacc肠道基因组，进一步地，也可以预先配置包含更多参考序列的资源库，例如依据待测样本来源的个体的状态、地域等因素选择或是测定组装出更接近的序列作为参考序列。根据本发明的一个实施例，各种微生物的参考基因组从公开的数据库中获得，通常的，一种微生物的参考基因组有多个版本，即一种微生物有多个公开的参考基因组。
67.reads与物种的参考基因组比对，比对上的可以被分为两部分：a)unique reads(u)：唯一比对上一个物种的基因组；称这些reads为unique reads。即，如果reads比对上的参考基因组均来自同一物种，定义这些reads为unique reads。b)multiple reads(m)：比对上一个以上物种的参考基因组，定义为multiple reads。即，如果reads比对上的参考基因组来自至少两种物种，定义这些reads为multiple reads。
68.根据本发明的一个实施例，按照下列公式确定所述标志微生物的丰度：
69.ab(s)＝ab(us) ab(ms)，其中其中，s表示所述标志微生物的编号，ab(s)表示所述标志微生物s的丰度，ab(us)＝us/ls，us为所述测序数据中与所述标志微生物s的参考基因组唯一比对的读段数目，ls为所述标志微生物s的参考基因组的总长度，ms为所述测序数据中与所述标志微生物s的参考基因组非唯一比对的读段的数目，i表示所述非唯一比对读段的编号，coi为所述第i读段对应的丰度系数，co
i,s
表示针对所述标志微生物s，所述非唯一比对的读段i的丰度系数，n为所述非唯一比对的读段i能够比对的微生物的总数，j表示所述非唯一比对的读段i能够比对的微生物的编号。上述对本发明任一实施例的利用标志微生物确定个体是否患有
小柳-原田综合症的方法的技术特征和优点的描述，同样适用本发明这一方面的装置，在此不再赘述。
70.根据本发明的又一个实施方式提供的一种装置，包括：计算机可读存储介质，其上存储有计算机程序，所述程序用于执行前面所述的一种确定个体是否患有小柳-原田综合症方法；以及一个或者多个处理器，用于执行所述计算机可读存储介质中的程序。
71.根据本发明的又一个实施方式提供的一种筛选药物的方法，所述药物用于治疗或者预防小柳-原田综合症，所述方法包括：将候选药物施用于受试者，检测施用前后，所述受试者粪便中标志微生物的丰度，所述标志微生物包括第一微生物集和第二微生物集中的至少一种菌种，其中，满足下列条件至少之一的候选药物适于用于治疗或者预防小柳-原田综合症：(1)进行所述施用后，所述第一微生物集中的至少一种菌种的所述丰度升高；和(2)进行所述施用后，所述第二微生物集中的至少一种菌种的所述丰度降低；其中，所述第一微生物集由以下菌种组成：固氮螺菌(azospirillum_sp.)_cag:260和梭菌(clostridium_sp.)_cag:349；所述第二微生物集由以下菌种组成：副普氏菌(paraprevotella_clara)、厚壁菌(firmicutes_bacterium)_cag:194和酪黄肠球菌(enterococcus_casseliflavus)。
72.利用本发明这一方面的生产或筛选治疗小柳-原田综合症的药物的方法，通过合理有效地应用确定的小柳-原田综合症生物标志物进行筛选，能够获得能扶持肠道有益菌的生长和/或抑制肠道潜在致病菌的药物，对于辅助减轻小柳-原田综合症的临床症状具有重要意义。
73.下面将对实施例作具体介绍。以下实施例中涉及的未特别交待的试剂、序列、软件及仪器，都是常规市售产品。
74.实施例1生物标志物的鉴定
75.本实施例中，发明人通过对41个小柳-原田综合症患者和38个健康对照的粪便样品进行研究，从而获得肠道菌群微生物微生物群落及功能成分特征。总的来说，发明人通过实验测序得到的216.6gb高质量测序数据构建了小柳-原田综合症患者(as)参照基因集及205.2g高通量数据构建健康人(lc)基因集，并且igc基因集构建一个更加完整的基因集。宏基因组分析显示，52个微生物物种与小柳-原田综合症疾病密切相关，其中28个菌在健康人的肠道微生物中富集，24个菌在小柳-原田综合症病人的肠道微生物中富集。
76.1、测序数据的获取：
77.小柳-原田综合症患者来自杭州重庆医科大学附属第一医院，实验共采集了41个中国小柳-原田综合症患者的粪便样品和38个健康对照的粪便样品，其中每个个体的新鲜粪便样品分成200mg/份，共5份，立即-80℃冰箱冷冻保存。
78.41个中国小柳-原田综合症患者的粪便样品和38个健康对照的粪便样品提取总dna。苯酚三氯甲烷处理提取dna方法提取dna。
79.2、构建dna文库及测序
80.dna建库按仪器制造商(illumina)的操作指南进行。对文库进行pe150 bp测序，illumina hiseq2000(illumina,san diego,ca)平台对79个样品的文库进行测序。每个样本平均产生5.7gb(sd.
±
0.74gb)高质量测序结果，总计421gb测序数据量。
81.参照图1所示的实验流程，鉴定小柳-原田综合症的相关生物标志物，其中省略的步骤或者细节为本领域技术人员所熟知，几个重要步骤介绍如下面所述。
82.3、微生物物种丰度分析
83.3.1序列优化统计
84.1)首先进行第一期测序，本期测序采集79个样品的数据，获得第一期的52个样品的测序数据以后，对其进行过滤，质控按以下标准进行：a)移除大于5个n碱基的reads；b)去除低质量碱基(q20)大于50％的reads；c)移除尾部低质量(q20)和n碱基。丢失成对的reads序列被认为是单条reads用于组装。
85.2)采用1)所述方法同样处理下载的健康人数据。
86.3)ftp://climb.genomics.cn/pub/10.5524/100001_101000/100064/1.genecatalogs/igc.fa.gz链接下载得到igc基因集。
87.3.2物种丰度分析
88.soapalign 2.21用于匹配针对冗余基因组的paired-end clean reads，这里，所称的冗余基因组来自各数据库中公开的的细菌的参考基因组，比对参数为
–
r 2
–
m 200
–
x 1000。reads与冗余基因组的比对结果，可分为两部分：a)unique reads(u)：reads只比对上一个物种的基因组；这些reads被定义为unique reads。即，如果这些基因组来自同一物种，发明人将这些reads定义为unique reads。b)multiple reads(m)：如果reads比对上两个及两个以上物种的基因组，定义为multiple reads。即，如果比对上的基因组来自不同物种，发明人定义这些reads为multiple reads。
89.对于物种s，其丰度为ab(s)，与特有的u reads和共享的m reads相关，丰度的计算方式如下：
90.ab(s)＝ab(us) ab(ms)，
91.其中，s表示所述标志微生物的编号，
92.ab(s)表示所述标志微生物s的丰度，
93.ab(us)＝us/ls，
94.us为所述测序数据中与所述标志微生物s的参考基因组唯一比对的读段数目，
95.ls为所述标志微生物s的参考基因组的总长度，
[0096][0097]ms
为所述测序数据中与所述标志微生物s的参考基因组非唯一比对的读段的数目，
[0098]
coi为所述第i读段对应的丰度系数，
[0099][0100]
i表示所述非唯一比对读段的编号，
[0101]
co
i,s
表示针对所述标志微生物s，
[0102]
所述非唯一比对的读段i的丰度系数，
[0103]
n为所述非唯一比对的读段i能够比对的微生物的总数，
[0104]
j表示所述非唯一比对的读段i能够比对的微生物的编号。
[0105]
每个样品中计算得到的物种丰度值，全部除以每个样品的丰度总数后，得到归一化后的物种丰度表。
[0106]
3.3筛选微生物物种标记物
[0107]
为了获得与小柳-原田综合症疾病密切相关的肠道微生物物种标记物，发明人利用小柳-原田综合症患者(as)组(41例)与正常人(hd)组(38例)两组肠道微生物物种丰度数据，在物种级别做了一个与疾病相关性的研究。基于步骤3.2所得到的物种丰度表，发明人设置标准如下：(1)小柳-原田综合症患者组或健康人组物种丰度的中位数必须大于0.00001；(2)通过结合benjamini hochberg的多重检验的wilcoxon秩和检验进行检验，得到每个物种在小柳-原田综合症人和健康个体样本间的相关性p值；(3)筛选阈值为p_values《0.05，利用上述参数进行筛选。发明人得到52个与小柳-原田综合症疾病密切相关的肠道微生物物种，其中，在小柳-原田综合症患者肠道内中富集的微生物有24个物种，在健康人中富集的有28个物种，这52个微生物物种标记物如表1所示。
[0108]
表1：
[0109]
[0110][0111]
实施例2微生物物种标记物的验证
[0112]
为了验证实施例1中的发现，发明人进一步分析验证群体中的14个健康人及14个小柳-原田综合症患者的粪便样本中的所述52种菌属的丰度，并根据各物种在健康组及疾病组的富集情况对所述52个微生物物种标记物做出删选，验证群体的dna提取、测序以及物种丰度的分析参照实施例1进行。
[0113]
验证结果如下：上述富集在健康人群中的28个物种，其中，2个在验证集中得到高质量的验证(p_values《0.05)，健康人富集的微生物物种标志物验证的均值和p值结果如表2所示，分别为固氮螺菌(azospirillum_sp.)_cag:260、梭菌(clostridium_sp.)_cag:349。固氮螺菌多于热带植物根际固氮作用。梭菌是肠道内一大类正常的兼性厌氧细菌,属于厚壁菌门,能产生芽孢.目前已经报道的梭菌有十几种,其中有一些主要用于疾病的预防、诊断或治疗,还有部分被证实与人类某些疾病的发生、发展、预后转归息息相关,比如抗生素相关性腹泻、伪膜性结肠炎等疾病.其主要包括有益的梭菌和有害的梭菌两大类,主要来源于肠道和粪便。
[0114]
表2：
[0115][0116]
对于上述富集在小柳-原田综合症患者中的24个物种，其中，3个在验证集中得到高质量的验证(p值《0.05)，小柳-原田综合症患者富集的微生物物种标志物验证的均值和p值具体结果如表3所示，分别为副普氏菌(paraprevotella_clara)、厚壁菌(firmicutes_bacterium)_cag:194和酪黄肠球菌(enterococcus_casseliflavus)。
[0117]
副普氏菌在小柳原田病患者的肠道富集(yang et al.,2020)。肠球菌通过逃避免
疫系统、黏附于宿主细胞、形成生物膜等方法引发疾病，毒力因子包括溶细胞素、明胶酶等。其基因组可塑性高，可通过水平基因转移获得新的毒力因子和耐药基因。
[0118]
表3：
[0119][0120]
发明人认为，可以将从健康人富集的2个微生物物种标记物，作为小柳-原田综合症疾病患病的反向指标，或作为治疗小柳-原田综合症进行研发的微生物制剂药物菌群成分，或作为检测小柳-原田综合症、监测小柳-原田综合症治疗进程的恢复指标；将小柳-原田综合症患者富集的3个微生物物种标记物，作为小柳-原田综合症疾病患病的正向指标，特别用于小柳-原田综合症疾病非创伤式的检测和诊断。
[0121]
发明人利用这5个微生物物种标记物，构建一个综合指标，估计roc(receiver-operating characteristic)曲线下面积auc，auc越大，表示诊断能力越高，评价综合得分对应其对小柳-原田综合症的诊断能力。通过对一期(第一阶段)的74个样品和二期(第二阶段)的33个样品进行评测，具体情况如图2所示，都表现出了很好的诊断能力，在一期得到auc＝78.6％，如图2a所示，置信区间为59.5％-97.6％；二期得到auc＝100.0％，如图2b所示，置信区间为100.0％-100.0％。
[0122]
实施例3个体状态的检测
[0123]
本实施例中，发明人利用45个粪便样本进行样本来源的个体状态的检测。
[0124]
参照实施例2的方法确定各粪便样本中表3所示的副普氏菌(paraprevotella_clara)、厚壁菌(firmicutes_bacterium)_cag:194和酪黄肠球菌(enterococcus_casseliflavus)的丰度，判断各样本中的这3种菌株丰度是否落入各自在疾病对照组或者健康对照组的丰度的95％的置信区间，判定这3种菌种的丰度均落入疾病组的对应区间的样本所对应的个体的状态为小柳-原田综合症患者，判定3种菌种的丰度均落入健康组的对应区间的样本所对应的个体的状态为非小柳-原田综合症患者。
[0125]
结果显示，利用本实施例所述方法能够对其中的45个样本进行个体状态判断，而且对这45个样本中的29个样本对应个体的状态的判断，与记录的该样本来源个体的状态一致。
[0126]
在本说明书的描述中，参考术语“一个实施例”、“一些实施例”、“示例”、“具体示例”、或“一些示例”等的描述意指结合该实施例或示例描述的具体特征、结构、材料或者特点包含于本发明的至少一个实施例或示例中。在本说明书中，对上述术语的示意性表述不必须针对的是相同的实施例或示例。而且，描述的具体特征、结构、材料或者特点可以在任一个或多个实施例或示例中以合适的方式结合。此外，在不相互矛盾的情况下，本领域的技术人员可以将本说明书中描述的不同实施例或示例以及不同实施例或示例的特征进行结合和组合。
[0127]
尽管上面已经示出和描述了本发明的实施例，可以理解的是，上述实施例是示例性的，不能理解为对本发明的限制，本领域的普通技术人员在本发明的范围内可以对上述实施例进行变化、修改、替换和变型。