肌痛性脑脊髓炎标志微生物及其应用的制作方法

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


背景技术：

2.肌痛性脑脊髓炎(myalgic encephalomyelitis,me)又称慢性疲劳综合症(chronic fatigue syndrome， cfs)，也称慢性疲劳免疫功能紊乱综合征(chronic fatigue immune dysfunction syndrome,cfids)。是一种累及机体多个系统的疾病，该病病程长，病人感到极度疲倦，肌痛，并且会引起认知功能障碍。各年龄段均存在患肌痛性脑脊髓炎的风险。肌痛性脑脊髓炎的发病机理尚不清楚，但其自身因素和外界因素已被确认。自身因素包括年龄、遗传易感性、家庭集聚性、激素状况等，外界因素包括传染、免疫接种、情绪紧张、机体损伤(如意外事故、手术)和环境中的毒素等。
3.肌痛性脑脊髓炎通常起病比较隐匿，症状表现为持续性疲倦、活动不适(类似流感症状)、肌痛、关节痛、神经痛(针刺样痛)、头痛、认知和睡眠障碍、其他神经系统症状、起立眩晕、平衡障碍对光和声音过敏，且病人可能会出现反复性喉痛和消化紊乱(如恶心、无食欲、消化不良、腹胀、腹泻和便秘)等症状。临床医生常常将其误诊为神经衰弱、更年期综合症、内分泌失调、神经官能症等，从而失去最佳治疗时机。
4.目前尚无确定的化验诊断标准，确诊取决于对典型症状的识别以及排除其他疾病。临床评估包括全面采集病史、睡眠评估、机体检查、心理健康评估(包括精神状态和心理学评估)。基本检查包括全血化验、尿常规、c-反应蛋白、肌酐、电解质、肝功能、甲状腺功能、肌酸磷酸激酶和类风湿因子检测。应排除的疾病有甲状腺病、慢性传染病、炎症、代谢性疾病、抑郁症。在体检发现异常结果时，应针对特异性症状进行特殊化验检查以排除其他疾病。一般在6个月内做一次鉴别诊断，再满3个月再次检查并调整治疗方案。
5.随着人体基因组测序完成及高通量测序技术的高速发展，基因筛查成为肌痛性脑脊髓炎诊断的方向，对于发现肌痛性脑脊髓炎潜在人群很有优势。此外，有研究表明肌痛性脑脊髓炎患者肠道微生物的水平出现异常。然而，到目前为止，还没有肌痛性脑脊髓炎患者肠道微生物标记物的相关报道。


技术实现要素：

6.本技术是基于发明人对以下事实和问题的发现和认识作出的：
7.本技术的申请人通过前期大量研究，意外地发现一些微生物可以作为检测肌痛性脑脊髓炎的标志微生物，合理有效地应用标志微生物，扶持肠道有益菌的生长，抑制肠道潜在致病菌，可以治疗或减轻肌痛性脑脊髓炎的临床症状。
8.为此，在本发明的第一方面，本发明提出了一种试剂盒。根据本发明的实施例，包
括适于检测第一微生物集中的至少一种菌种的试剂，所述第一微生物集由以下菌种组成：dorea_longicatena、厚壁菌 cag56(firmicutes_bacterium_)cag_56、栖粪杆菌(faecalibacterium_sp_)cag_82和霍氏真杆菌 (eubacterium_hallii_)cag_12。根据本发明具体实施例的试剂盒，可以准确的区分肌痛性脑脊髓炎患者和健康个体。
9.在本发明的第二方面，本发明提出了试剂在制备试剂盒中的用途，所述试剂适于检测第一微生物集中的至少一种菌种。根据本发明的实施例，所述试剂盒用于诊断肌痛性脑脊髓炎或检测肌痛性脑脊髓炎的治疗效果，所述第一微生物集由以下菌种组成：dorea_longicatena、厚壁菌cag56 (firmicutes_bacterium_)cag_56、栖粪杆菌(faecalibacterium_sp_)cag_82和霍氏真杆菌 (eubacterium_hallii_)cag_12。根据本发明具体实施例的试剂制备的试剂盒，可以准确的检测所述第一微生物集中的至少一种菌种，极准确的区分肌痛性脑脊髓炎患者和健康个体，由此，可以有效的在早期进行肌痛性脑脊髓炎诊断，或用于检测治疗过程中肌痛性脑脊髓炎的变化。
10.在本发明的第三方面，本发明提出了一种用于预防或者治疗肌痛性脑脊髓炎的药物组合物或者食品组合物。根据本发明的实施例，含有第一微生物集中的至少一种菌种，所述第一微生物集由以下菌种组成：dorea_longicatena、厚壁菌cag56(firmicutes_bacterium_)cag_56、栖粪杆菌(faecalibacterium_sp_) cag_82和霍氏真杆菌(eubacterium_hallii_)cag_12。根据本发明实施例的肌痛性脑脊髓炎标志微生物中第一微生物集的菌种可以非侵入性的在早期发现或辅助检测肌痛性脑脊髓炎，确定个体患有肌痛性脑脊髓炎的概率高低或者个体处于健康状态的概率高低；同时，提高肌痛性脑脊髓炎高风险人群或已肌痛性脑脊髓炎患者肠道内的所述第一微生物集中的各种菌种，可以降低患肌痛性脑脊髓炎的概率或减缓、治愈肌痛性脑脊髓炎，因此，所述包含所述第一微生物集中的至少一种菌种的药物或者食品组合物能够用于平衡肠道菌群，有效预防或治疗肌痛性脑脊髓炎。
11.在本发明的第四方面，本发明提出了一种确定个体是否患有肌痛性脑脊髓炎的方法。根据本发明的实施例，包括：(1)确定所述个体的粪便样本中标志微生物的丰度，所述标志微生物包括第一微生物集和第二微生物集中的至少一种菌种；(2)将步骤(1)中得到的所述丰度与预定的阈值进行比较，以便确定所述个体是否患有肌痛性脑脊髓炎；其中，所述第一微生物集由以下菌种组成：dorea_longicatena、厚壁菌cag56(firmicutes_bacterium_)cag_56、栖粪杆菌(faecalibacterium_sp_)cag_82和霍氏真杆菌(eubacterium_hallii_)cag_12；所述第二微生物集由以下菌种组成：毛螺菌(anaerostipes_sp_) cag_276、梭状芽孢杆菌(clostridium_clostridioforme)、粪厌氧棒杆菌(anaerostipes_caccae)和梭状芽孢杆菌(clostridium_clostridioforme_)cag_511。根据本发明实施例的方法可以依据个体的粪便样本中的所述标志微生物的丰度确定个体是否患有肌痛性脑脊髓炎，所述标志微生物是发明人对大量已知状态的粪便样本进行验证，通过差异比较分析各种肠道微生物在肌痛性脑脊髓炎组和健康组粪便样本中的丰度，而确定下来的。
12.在本发明的第五方面，本发明提出了一种确定个体是否患有肌痛性脑脊髓炎的装置。根据本发明的实施例，包括：丰度确定单元，用于确定所述个体的粪便样本中标志微生物的丰度，所述标志微生物包括第一微生物集和第二微生物集中的至少一种菌种；比较单元，用于将所得到的所述丰度与预定的阈值进行比较，以便确定所述个体是否患有肌痛性
脑脊髓炎；其中，所述第一微生物集由以下菌种组成： dorea_longicatena、厚壁菌cag56(firmicutes_bacterium_)cag_56、栖粪杆菌(faecalibacterium_sp_) cag_82和霍氏真杆菌(eubacterium_hallii_)cag_12；所述第二微生物集由以下菌种组成：毛螺菌 (anaerostipes_sp_)cag_276、梭状芽孢杆菌(clostridium_clostridioforme)、粪厌氧棒杆菌 (anaerostipes_caccae)和梭状芽孢杆菌(clostridium_clostridioforme_)cag_511。所述标志微生物是发明人通过差异比较分析各种肠道微生物在肌痛性脑脊髓炎患者和健康人群的粪便样本中的丰度并经过分析，以及将大量已知状态的粪便样本的验证而确定下来的，根据本发明实施例的装置可以准确确定个体是否为肌痛性脑脊髓炎的高风险人群或肌痛性脑脊髓炎患者。
13.在本发明的第六方面，本发明提出了一种装置。根据本发明的实施例，包括：计算机可读存储介质，其上存储有计算机程序，所述程序用于执行第四方面所述的方法；以及一个或者多个处理器，用于执行所述计算机可读存储介质中的程序。根据本发明实施例的装置可以准确确定个体是否为肌痛性脑脊髓炎的高风险人群或肌痛性脑脊髓炎患者。
14.在本发明的第七方面，本发明提出了一种筛选药物的方法。根据本发明的实施例，所述药物用于治疗或者预防肌痛性脑脊髓炎，所述方法包括：将候选药物施用于受试者，检测施用前后，所述受试者粪便中标志微生物的丰度，所述标志微生物包括第一微生物集和第二微生物集中的至少一种菌种，其中，满足下列条件至少之一的候选药物适于用于治疗或者预防肌痛性脑脊髓炎：(1)进行所述施用后，所述第一微生物集中的至少一种菌种的所述丰度升高；和(2)进行所述施用后，所述第二微生物集中的至少一种菌种的所述丰度降低；其中，所述第一微生物集由以下菌种组成：dorea_longicatena、厚壁菌 cag56(firmicutes_bacterium_)cag_56、栖粪杆菌(faecalibacterium_sp_)cag_82和霍氏真杆菌 (eubacterium_hallii_)cag_12；所述第二微生物集由以下菌种组成：毛螺菌(anaerostipes_sp_) cag_276、梭状芽孢杆菌(clostridium_clostridioforme)、粪厌氧棒杆菌(anaerostipes_caccae)和梭状芽孢杆菌(clostridium_clostridioforme_)cag_511。根据本发明实施例的方法可以生产或筛选出促进所述标志微生物中第一微生物集中各种菌种生长和/或抑制肠道标志微生物中第二微生物集中的各种菌种生长的药物，对于辅助减轻肌痛性脑脊髓炎的临床症状具有重要意义。
15.本发明的附加方面和优点将在下面的描述中部分给出，部分将从下面的描述中变得明显，或通过本发明的实践了解到。
附图说明
16.本发明的上述和/或附加的方面和优点从结合下面附图对实施例的描述中将变得明显和容易理解，其中：
17.图1是根据本发明实施例的筛选肌痛性脑脊髓炎标志微生物的实验分析流程示意图；以及
18.图2是根据本发明实施例的标志微生物综合指标auc评价结果示意图，其中，横坐标specificity 表示特异度，即预测为阳性且实际为阳性，真阳性，纵坐标sensitivity表示敏感度，即真阴性：
19.2-a为第一期74个样品数据roc曲线下auc值和置信区间结果图；
20.2-b为第二期26个样品数据roc曲线下auc值和置信区间结果图。
具体实施方式
21.下面详细描述本发明的实施例，所述实施例的示例在附图中示出。下面通过参考附图描述的实施例是示例性的，旨在用于解释本发明，而不能理解为对本发明的限制。
22.此外，术语“第一”、“第二”仅用于描述目的，而不能理解为指示或暗示相对重要性或者隐含指明所指示的技术特征的数量。由此，限定有“第一”、“第二”的特征可以明示或者隐含地包括至少一个该特征。在本发明的描述中，“多个”的含义是至少两个，例如两个，三个等，除非另有明确具体的限定。
23.术语“任选地”仅用于描述目的，而不能理解为指示或暗示相对重要性。由此，限定有“任选地”的特征可以明示或者隐含地包括或不包括该特征。
24.生物学标志物是从生物学介质中可以检测到的细胞/生物化学或分子改变。生物学介质包括各种体液、组织、细胞、粪便、头发、呼气等。
25.微生物的丰度是指在某一微生物群体中该种微生物的丰富程度，例如在肠道微生物群体中的该种微生物程度，可表示为该种微生物在该群体中的含量。
26.根据本发明的一个实施方案提供的一种试剂盒，包括适于检测第一微生物集中的至少一种菌种的试剂，所述第一微生物集由以下菌种组成：dorea_longicatena、厚壁菌cag56(firmicutes_bacterium_) cag_56、栖粪杆菌(faecalibacterium_sp_)cag_82，霍氏真杆菌(eubacterium_hallii_)cag_12。
27.根据本发明的具体实施方案，所述试剂盒进一步包括适于检测第二微生物集中的至少一种菌种的试剂，所述第二微生物集由以下菌种组成：毛螺菌(anaerostipes_sp_)cag_276、梭状芽孢杆菌 (clostridium_clostridioforme)、粪厌氧棒杆菌(anaerostipes_caccae)、梭状芽孢杆菌 (clostridium_clostridioforme_)cag_511。
28.根据本发明的具体实施方案，所述试剂盒包括适于检测所述第一微生物集中全部所述菌种的试剂。
29.根据本发明的具体实施方案，所述试剂盒包括适于检测所述第二微生物集中全部所述菌种的试剂。
30.根据本发明的具体实施方案，所述标志微生物是发明人通过对大量患肌痛性脑脊髓炎个体和大量健康对照个体的粪便样本中的微生物的丰度的差异比较分析、以及验证，而确定下来的，明确了肠道微生物中肌痛性脑脊髓炎相关的标志微生物。利用包含检测所述标志微生物的试剂的试剂盒能够确定个体处于患有肌痛性脑脊髓炎状态的概率高低或者处于健康状态的概率高低，能够用于非侵入性的早期发现或辅助检测肌痛性脑脊髓炎。
31.根据本发明的具体实施方案，所述适于检测所述第一微生物集或第二微生物集的试剂不受特别限制，任何可以检测所述微生物菌种的试剂均包含在本发明的范围内，如通过形态学特征、生理生化反应特征、生态学特征以及血清学反应、对噬菌体的敏感性、分子生物学等方面检测所述微生物菌种的试剂，具体的，如抗体、酶、核酸分子等。
32.本文中，所述微生物形态学特征指：在显微镜下观察微生物外形大小、形状、排列等，细胞构造，革兰氏染色反应，能否运动、鞭毛着生部位和数目，有无芽孢和荚膜、芽孢的大小和位置，放线菌和真菌的繁殖器官的形状、构造，孢子的数目、形状、大小、颜色和表面
特征等。
33.本文中，所述微生物生理生化反应特征指：微生物利用物质的能力、代谢产物的特殊性，如是否产生h2s、吲哚、co2、醇、有机酸，能否还原硝酸盐，能否使牛奶凝固、冻化等、生长环境(温度、湿度、氧气及二氧化碳等气体的浓度、ph、是否耐高渗、是否有嗜盐性等)、与其它生物间的关系(如：共生、寄生、宿主范围及致病情况)等。
34.本文中，所述微生物血清学反应指：利用抗原与抗体的高度敏感特异性反应，来鉴别相似的菌种，或对同种微生物分型，如用已知菌种、型或菌株制成的抗血清，与待鉴定的所述微生物是否发生特异性的血清学反应来鉴定微生物。
35.本文中，分子生物学方法检测微生物主要包括：利用pcr技术、高通量测序等方法。
36.根据本发明提供的试剂在制备试剂盒中的用途，所述试剂适于检测第一微生物集中的至少一种菌种，所述试剂盒用于诊断肌痛性脑脊髓炎或者检测肌痛性脑脊髓炎的治疗效果，所述第一微生物集由以下菌种组成：dorea_longicatena、厚壁菌cag56(firmicutes_bacterium_)cag_56、栖粪杆菌 (faecalibacterium_sp_)cag_82，霍氏真杆菌(eubacterium_hallii_)cag_12。
37.根据本发明的具体实施方案，所述标志微生物是发明人通过对大量患肌痛性脑脊髓炎个体和大量健康对照个体的粪便样本中的微生物的丰度的差异比较分析、以及验证，而确定下来的，明确了肠道微生物中肌痛性脑脊髓炎相关的微生物标志物。利用检测所述标志微生物的试剂能够确定个体患有肌痛性脑脊髓炎的概率高低或者处于健康状态的概率高低，能够用于非侵入性的早期发现或辅助检测肌痛性脑脊髓炎。
38.根据本发明一些具体的实施方案，所述试剂进一步适于检测第二微生物集中的至少一种菌种，所述第二微生物集由以下菌种组成：毛螺菌(anaerostipes_sp_)cag_276、梭状芽孢杆菌 (clostridium_clostridioforme)、粪厌氧棒杆菌(anaerostipes_caccae)、梭状芽孢杆菌 (clostridium_clostridioforme_)cag_511。
39.根据本发明的具体实施方案，所述适于检测所述第一微生物集或第二微生物集的试剂不受特别限制，可以检测所述微生物菌种的试剂均包含在本发明的范围内，如通过形态学特征、生理生化反应特征、生态学特征以及血清学反应、对噬菌体的敏感性、分子生物学等方面检测所述微生物菌种的试剂等，具体的，如抗体、酶、核酸分子。
40.本文中，所述微生物形态学特征指：在显微镜下观察微生物外形大小、形状、排列等，细胞构造，革兰氏染色反应，能否运动、鞭毛着生部位和数目，有无芽孢和荚膜、芽孢的大小和位置，放线菌和真菌的繁殖器官的形状、构造，孢子的数目、形状、大小、颜色和表面特征等。
41.本文中，所述微生物生理生化反应特征指：微生物利用物质的能力、代谢产物的特殊性，如是否产生h2s、吲哚、co2、醇、有机酸，能否还原硝酸盐，能否使牛奶凝固、冻化等、生长环境(适宜生长的温度、湿度、氧气及二氧化碳等气体的浓度、ph、是否耐高渗、是否有嗜盐性等)、与其它生物间的关系(如：共生、寄生、宿主范围及致病情况)等。
42.本文中，所述微生物血清学反应指：利用抗原与抗体的高度敏感特异性反应，来鉴别相似的菌种，或对同种微生物分型，如用已知菌种、型或菌株制成的抗血清，与待鉴定的所述微生物是否发生特异性的血清学反应来鉴定微生物。
43.本文中，分子生物学方法检测微生物主要包括：利用pcr技术、高通量测序等方法。
44.根据本发明提供的一种用于预防或者治疗肌痛性脑脊髓炎的药物组合物或者食品组合物，含有第一微生物集中的至少一种菌种，所述第一微生物集由以下菌种组成：dorea_longicatena、厚壁菌cag56 (firmicutes_bacterium_)cag_56、栖粪杆菌(faecalibacterium_sp_)cag_82，霍氏真杆菌 (eubacterium_hallii_)cag_12。
45.上述标志微生物是发明人通过差异比较分析各种肠道微生物在肌痛性脑脊髓炎疾病组和健康组的粪便样本中的丰度，以及经过大量已知状态的粪便样本的验证，而确定下来的。所述标志微生物中第一微生物集中的菌种相较于肌痛性脑脊髓炎患者群体，在健康群体组中显著富集，所述显著富集是指与在肌痛性脑脊髓炎患者组中的丰度相比，上述菌种在健康组中的丰度均具有统计意义地高于或者明显地、实质性地高于在肌痛性脑脊髓炎患者组中的丰度；能够使该部分菌种丰度提高的物质能够用于治疗肌痛性脑脊髓炎或者益于肌痛性脑脊髓炎患者服用，能够使其丰度提高的物质不限于治疗肌痛性脑脊髓炎的药物和有益肠道菌群平衡的功能性食品。因此，该实施方式提供的标志微生物能够用于制备治疗肌痛性脑脊髓炎的药物和/或用于制备益于平衡肠道菌群的功能性食品、保健药等，所述药物或食品可有效治疗或缓解肌痛性脑脊髓炎。
46.根据本发明提供的一种确定个体是否患有肌痛性脑脊髓炎的方法，包括步骤(1)和(2)。
47.(1)确定所述个体的粪便样本中的标志微生物的丰度。
48.所述标志微生物包括第一微生物集和第二微生物集中的至少一种菌种。其中，所述第一微生物集由以下菌种组成：dorea_longicatena、厚壁菌cag56(firmicutes_bacterium_)cag_56、栖粪杆菌 (faecalibacterium_sp_)cag_82，霍氏真杆菌(eubacterium_hallii_)cag_12；所述第二微生物集由以下菌种组成：毛螺菌(anaerostipes_sp_)cag_276、梭状芽孢杆菌(clostridium_clostridioforme)、粪厌氧棒杆菌(anaerostipes_caccae)、梭状芽孢杆菌(clostridium_clostridioforme_)cag_511。
49.根据本发明的一些具体的实施方案，步骤(1)进一步包括：获得所述个体的粪便样本中的核酸测序数据；将所述测序数据与参考基因组进行比对；基于所述比对的结果，确定所述标志微生物的丰度。
50.根据本发明的具体实施方案，在步骤(1)中，按照下列公式确定所述标志微生物的丰度： 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能够比对的微生物的编号。
51.比对可以利用已知比对软件进行，例如soap、bwa和teramap等，在比对过程中，一般对比对参数进行设置，设置一个或者一对读段(reads)最多允许有s个碱基错配
(mismatch)，例如设置s≤2，若reads中有超过s个碱基发生错配，则视为该reads无法比对到(比对上)该组装片段上。所述的获得的比对结果包含各条读段与各物种的参考基因组的比对情况，包括读段是否能够比对上某个或某些物种的参考基因组、只唯一比对到一种物种还是比对到多种物种的参考基因组、比对到物种的参考基因组的位置、比对到物种参考基因组的唯一位置还是多个位置等信息。
52.所述菌种/微生物的参考基因组指预先确定的该微生物物种的序列，可以是预先获得的待测样本所属或者所包含的生物类别的任意参考模板，例如，目标是待测样本中的微生物，参考序列可选择ncbi数据库中的各种微生物的参考基因组和/或hmp、metahit项目公开的dacc肠道参考基因组，进一步地，也可以预先配置包含更多参考序列的资源库，例如依据待测样本来源的个体的状态、地域等因素选择或是测定组装出更接近的序列作为参考序列。根据本发明的一个实施例，各种微生物的参考基因组从公开的数据库中获得，通常的，一种微生物的参考基因组有多个版本，即一种微生物有多个公开的参考基因组。
53.reads与物种的参考基因组比对，比对上的可以被分为两部分：a)unique reads(u)：唯一比对上一个物种的参考基因组；称这些reads为unique reads。即，如果reads比对上的参考基因组均来自同一物种，定义这些reads为unique reads；b)multiple reads(m)：比对上一个以上物种的参考基因组，定义为multiple reads。即，如果reads比对上的参考基因组来自至少两种物种，定义这些reads为multiple reads。
54.(2)丰度比较，以确定个体是否患有肌痛性脑脊髓炎。
55.根据本发明的一个实施方案，将步骤(1)中得到的所述丰度与预定的阈值进行比较，以便确定所述个体是否患有肌痛性脑脊髓炎。
56.根据本发明的一些具体实施方案，所述阈值为预先设定的。将标志微生物在待测个体样本中的丰度与所述的阈值比较，进行待测个体的状态的确定。所述阈值可以为一数值或者数值范围，阈值设定方式不受特别限制，能够对所述个体的是否患炎症性肠病进行判断的阈值均可使用，例如，基于已知患病或健康状态个体中的标志微生物的丰度均值，该微生物对应的阈值可以设为该丰度均值的95％的置信区间 (confidence interval)。
57.所述的置信区间是指由样本统计量所构造的总体参数的估计区间。在统计学中，一个概率样本的置信区间是对这个样本的某个总体参数的区间估计。置信区间展现的是这个参数的真实值有一定概率落在测量结果的周围的程度。置信区间给出的是被测量参数的测量值的可信程度，即前面所要求的“一定概率”，这个概率被称为置信水平。
58.根据本发明的一些具体实施方案，当步骤(1)中确定的标志微生物的丰度达到所述患肌痛性脑脊髓炎丰度阈值，未达到所述不患肌痛性脑脊髓炎丰度阈值时，确定所述个体患有肌痛性脑脊髓炎，当(1) 中确定的标志微生物的丰度达到不患肌痛性脑脊髓炎丰度阈值，未达到患肌痛性脑脊髓炎丰度阈值时，确定所述个体不患肌痛性脑脊髓炎。
59.需要说明的是，根据目的或要求不同，可能对确定个体状态结果的可信程度、目的有不同的要求，本领域技术人员可以选择不同的显著性水平或阈值。
60.该方法基于检测个体的粪便样本中的标志微生物中的各种菌种的丰度，分别将检测确定的各种菌种的丰度与其阈值进行比较，依据获得的比较结果能够确定个体为肌痛性脑脊髓炎个体或者为健康个体的概率。为早期发现肌痛性脑脊髓炎提供一种非侵入性的辅助检测或者辅助干预治疗的方法。
61.以上任一实施例中的利用标志微生物确定个体是否患有肌痛性脑脊髓炎的方法的全部或部分步骤，可以利用包含可拆分的相应单元功能模块的装置/系统来施行，或者将方法程序化、存储于机器可读介质，利用机器运行该可读介质来实现。
62.根据本发明提供的一种确定个体是否患有肌痛性脑脊髓炎的装置，该装置包括：丰度确定单元，用于确定所述个体的粪便样本中标志微生物的丰度；比较单元，用于将所得到的所述丰度与预定的阈值进行比较，以便确定所述个体是否患有肌痛性脑脊髓炎；其中，所述第一微生物集由以下菌种组成：dorea_longicatena、厚壁菌cag56(firmicutes_bacterium_)cag_56、栖粪杆菌(faecalibacterium_sp_) cag_82，霍氏真杆菌(eubacterium_hallii_)cag_12；所述第二微生物集由以下菌种组成：毛螺菌 (anaerostipes_sp_)cag_276、梭状芽孢杆菌(clostridium_clostridioforme)、粪厌氧棒杆菌 (anaerostipes_caccae)、梭状芽孢杆菌(clostridium_clostridioforme_)cag_511。上述对本发明任一实施例的利用标志微生物确定个体是否患有肌痛性脑脊髓炎的方法的技术特征和优点的描述，同样适用本发明这一方面的装置，在此不再赘述。
63.根据本发明的实施方案，所述丰度确定单元适于通过下列步骤确定所述丰度：获得所述个体的粪便样本中的核酸测序数据；将所述测序数据与参考基因组进行比对；基于所述比对的结果，确定所述标志微生物的丰度。
64.比对可以利用已知比对软件进行，例如soap、bwa和teramap等，在比对过程中，一般对比对参数进行设置，设置一个或者一对读段(reads)最多允许有s个碱基错配(mismatch)，例如设置s≤2，若reads中有超过s个碱基发生错配，则视为该reads无法比对到(比对上)该组装片段上。所述的获得的比对结果包含各条读段与各物种的参考基因组的比对情况，包括读段是否能够比对上某个或某些物种的参考基因组、只唯一比对到一种物种还是比对到多种物种的参考基因组、比对到物种的参考基因组的位置、比对到物种参考基因组的唯一位置还是多个位置等信息。根据本发明的一个实施例，利用 soapalign 2.21进行比对，设置参数为
–
r 2
–
m 100
–
x 1000。
65.所述微生物的参考基因组指预先确定的该微生物物种的序列，可以是预先获得的待测样本所属或者所包含的生物类别的任意参考模板，例如，目标是待测样本中的微生物，参考序列可选择ncbi数据库中的各种微生物的参考基因组和/或hmp、metahit项目公开的dacc肠道基因组，进一步地，也可以预先配置包含更多参考序列的资源库，例如依据待测样本来源的个体的状态、地域等因素选择或是测定组装出更接近的序列作为参考序列。根据本发明的一个实施例，各种微生物的参考基因组从公开的数据库中获得，通常的，一种微生物的参考基因组有多个版本，即一种微生物有多个公开的参考基因组。
66.reads与物种的参考基因组比对，比对上的可以被分为两部分：a)unique reads(u)：唯一比对上一个物种的基因组；称这些reads为unique reads。即，如果reads比对上的参考基因组均来自同一物种，定义这些reads为unique reads。b)multiple reads(m)：比对上一个以上物种的参考基因组，定义为multiplereads。即，如果reads比对上的参考基因组来自至少两种物种，定义这些reads为multiple reads。
67.根据本发明的一个实施方案，按照下列公式确定所述标志微生物的丰度：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能够比对的微生物的编号。上述对本发明任一实施例的利用标志微生物确定个体是否患有肌痛性脑脊髓炎的方法的技术特征和优点的描述，同样适用本发明这一方面的装置，在此不再赘述。
68.根据本发明提供的一种装置，包括：计算机可读存储介质，其上存储有计算机程序，所述程序用于执行前面所述的一种确定个体是否患有肌痛性脑脊髓炎方法；以及一个或者多个处理器，用于执行所述计算机可读存储介质中的程序。
69.根据本发明提供的一种筛选药物的方法，所述药物用于治疗或者预防肌痛性脑脊髓炎，所述方法包括：将候选药物施用于受试者，检测施用前后，所述受试者粪便中标志微生物的丰度，所述标志微生物包括第一微生物集和第二微生物集中的至少一种菌种，其中，满足下列条件至少之一的候选药物适于用于治疗或者预防肌痛性脑脊髓炎：(1)进行所述施用后，所述第一微生物集中的至少一种菌种的所述丰度升高；和(2)进行所述施用后，所述第二微生物集中的至少一种菌种的所述丰度降低；其中，所述第一微生物集由以下菌种组成：dorea_longicatena、厚壁菌cag56(firmicutes_bacterium_)cag_56、栖粪杆菌(faecalibacterium_sp_)cag_82，霍氏真杆菌(eubacterium_hallii_)cag_12；所述第二微生物集由以下菌种组成：毛螺菌(anaerostipes_sp_)cag_276、梭状芽孢杆菌(clostridium_clostridioforme)、粪厌氧棒杆菌(anaerostipes_caccae)、梭状芽孢杆菌(clostridium_clostridioforme_)cag_511。
70.利用本发明这一方面的生产或筛选治疗肌痛性脑脊髓炎的药物的方法，通过合理有效地应用确定的肌痛性脑脊髓炎生物标志物进行筛选，能够获得能扶持肠道有益菌的生长和/或抑制肠道潜在致病菌的药物，对于辅助减轻肌痛性脑脊髓炎的临床症状具有重要意义。
71.下面将对实施例作具体介绍。以下实施例中涉及的未特别交待的试剂、序列、软件及仪器，都是常规市售产品。
72.实施例1生物标志物的鉴定
73.本实施例中，发明人通过对37个肌痛性脑脊髓炎患者和37个健康对照的粪便样品进行研究，从而获得肠道菌群微生物微生物群落及功能成分特征。总的来说，发明人通过参考文献得到约116.6gb高质量数据构建了肌痛性脑脊髓炎患者(as)参照基因集，以及115.2gb高通量数据构建健康人(lc)基因集，且igc基因集构建一个更加完整的基因集。宏基因组分析显示，18个微生物物种与肌痛性脑脊髓炎病密切相关，其中12个菌在健康人的肠道微生物中富集，6个菌在肌痛性脑脊髓炎病患者的肠道微生物中富集。
74.1、测序数据的获取：
75.肌痛性脑脊髓炎患者来自美国四个地点，分别为纽约州纽约市、犹他州盐湖城、内
华达州因克莱恩村；和佛罗里达州迈阿密招募的37名me/cfs病例和37名健康对照者，数据下载的ncbi的项目编号是srp102150。
76.参照图1所示的实验流程，鉴定肌痛性脑脊髓炎的相关生物标志物，其中省略的步骤或者细节为本领域技术人员所熟知，几个重要步骤介绍如下面所述。
77.2、微生物物种丰度分析
78.2.1序列优化统计
79.1)首先进行第一期测序，本期测序采集74个样品的数据，获得第一期的74个样品的测序数据以后，对其进行过滤，质控按以下标准进行：a)移除大于5个n碱基的reads；b)去除低质量碱基(q20) 大于50％的reads；c)移除尾部低质量(q20)和n碱基。丢失成对的reads序列被认为是单条reads用于组装。
80.2)采用1)所述方法同样处理下载的健康人数据。
81.3)从ftp.cngb.org/pub/sciraid/microbiome/humangut_9.9m/genecatalog/igc.fa.gz下载得到igc基因集。
82.2.2物种丰度分析
83.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。
84.对于物种s，其丰度为ab(s)，与特有的u reads和共享的m reads相关，丰度的计算方式如下：
85.ab(s)＝ab(us) ab(ms)，
86.其中，s表示所述标志微生物的编号，
87.ab(s)表示所述标志微生物s的丰度，
88.ab(us)＝us/ls，
89.us为所述测序数据中与所述标志微生物s的参考基因组唯一比对的读段数目，
90.ls为所述标志微生物s的参考基因组的总长度，
[0091][0092]ms
为所述测序数据中与所述标志微生物s的参考基因组非唯一比对的读段的数目，
[0093]
coi为所述第i读段对应的丰度系数，
[0094][0095]
i表示所述非唯一比对读段的编号，
[0096]
co
i,s
表示针对所述标志微生物s，
[0097]
所述非唯一比对的读段i的丰度系数，
[0098]
n为所述非唯一比对的读段i能够比对的微生物的总数，
[0099]
j表示所述非唯一比对的读段i能够比对的微生物的编号。
[0100]
每个样品中计算得到的物种丰度值，全部除以每个样品的丰度总数后，得到归一化后的物种丰度表。
[0101]
2.3筛选微生物物种标记物
[0102]
为了获得与肌痛性脑脊髓炎疾病密切相关的肠道微生物物种标记物，发明人利用肌痛性脑脊髓炎患者(as)组(37例)与正常人(hd)组(37例)两组肠道微生物物种丰度数据，在物种级别做了一个与疾病相关性的研究。基于步骤2.2所得到的物种丰度表，发明人设置标准如下：(1) 肌痛性脑脊髓炎患者组或健康人组物种丰度的中位数必须大于0.000001；(2)通过结合benjaminihochberg的多重检验的wilcoxon秩和检验进行检验，得到每个物种和肌痛性脑脊髓炎疾病的相关性p值；(3)筛选的阈值p_values《0.005，利用上述参数进行筛选。发明人得到18个与肌痛性脑脊髓炎疾病密切相关的肠道微生物物种，其中，在肌痛性脑脊髓炎患者肠道内中富集的微生物有12个物种，在健康人中富集的有6个物种，这18个微生物物种标记物如表1所示。
[0103]
表1：
[0104][0105][0106]
实施例2微生物物种标记物的验证
[0107]
为了验证实施例1中的发现，发明人进一步分析验证群体中的13个健康人及13个肌痛性脑脊髓炎患者的粪便样本中的所述18种菌属的丰度，并根据各物种在健康组及疾病组的富集情况对所述18个微生物物种标记物做出删选，物种丰度的分析参照实施例1进行。
[0108]
验证结果如下：上述富集在健康人群中的6个物种，4个在验证集中得到高质量的验证 (p_values《0.05)，健康人富集的微生物物种标志物验证的均值和p值结果如表2所示，其中，克罗恩病人术后缓解患者与dorea longicaten菌群呈现强相关性。栖粪杆菌与特
发性钙结石的患病呈现负相关性霍氏真杆菌是另一种产丁酸的细菌，可以提高胰岛素敏感性并增强糖尿病小鼠中的能量消耗，这表明它可能在疾病中具有保护作用。
[0109]
表2：
[0110][0111]
对于上述富集在肌痛性脑脊髓炎患者中的12个物种中，其中4个在验证集中得到高质量的验证(p值《0.05)，验证富集在肌痛性脑脊髓炎患者中的4个微生物标志物验证的p值和q值情况在验证集数据中的情况如表3所示。
[0112]
已有研究报道，食欲较差的老年人的肠道菌群中毛螺菌属的丰度显著降低，而在老年女性中，食欲较差者的肌肉力量显著降低，提示毛罗菌属的丰度与肌肉力量具有一定的联系。
[0113]
表3：
[0114][0115]
发明人认为，可以将从健康人富集的4个微生物物种标记物，作为肌痛性脑脊髓炎疾病患病的反向指标，或作为治疗肌痛性脑脊髓炎进行研发的微生物制剂药物菌群成分，或作为检测肌痛性脑脊髓炎、监测肌痛性脑脊髓炎治疗进程的恢复指标；将肌痛性脑脊髓炎患者富集的4个微生物物种标记物，作为肌痛性脑脊髓炎疾病患病的正向指标，特别用于肌痛性脑脊髓炎疾病非创伤式的检测和诊断。
[0116]
发明人利用这8个微生物物种标记物，构建一个综合指标，估计roc(receiver-operatingcharacteristic)曲线下面积auc，auc越大，表示诊断能力越高，评价综合得分对应其对肌痛性脑脊髓炎的诊断能力。通过对一期(第一阶段)的74个样品和二期(第二阶段)的26个样品进行评测，具体情况如图2所示，都表现出了很好的诊断能力，在一期得到auc＝84.8％，如图2-a 所示，置信区间为68.8％-100.0％；二期得到auc＝80.0％，如图2-b所示，置信区间为47.2％-100.0％。
[0117]
实施例3个体状态的检测
[0118]
本实施例中，发明人利用45个粪便样本进行样本来源的个体状态的检测。
[0119]
参照实施例2的方法确定各粪便样本中表3所示的毛螺菌(anaerostipes_sp_)cag_276、梭状芽孢杆菌(clostridium_clostridioforme)、粪厌氧棒杆菌(anaerostipes_caccae)、梭状芽孢杆菌 (clostridium_clostridioforme_)cag_511的丰度，判断各样本中的这4种菌株丰度是否落入各自在疾病对照组或者健康对照组的丰度的95％的置信区间，判定这4种菌种的丰度均落入疾病组的对应区间的样本所对应的个体的状态为肌痛性脑脊髓炎患者，判定4种菌种的丰度均落入健康组的对应区间的样本所对应的个体的状态为非
肌痛性脑脊髓炎患者。
[0120]
结果显示，利用本实施例所述方法能够对其中的45个样本进行个体状态判断，而且对这45 个样本中的31个样本对应个体的状态的判断，与记录的该样本来源个体的状态一致。
[0121]
另外，发明人发现对表2和表3中的物种联合检测，例如检测表3中的物种标志物被富集，同时表2中物种标记物不被富集，能够更准确的判断发现肌痛性脑脊髓炎患者或易感人群。
[0122]
在利用标志物治疗肌痛性脑脊髓炎的方案中，发明人发现使表3中的物种标志物生长得到抑制或者清除，同时使表2中物种标记物被富集，治疗效果极佳。
[0123]
在本说明书的描述中，参考术语“一个实施例”、“一些实施例”、“示例”、“具体示例”、或“一些示例”等的描述意指结合该实施例或示例描述的具体特征、结构、材料或者特点包含于本发明的至少一个实施例或示例中。在本说明书中，对上述术语的示意性表述不必须针对的是相同的实施例或示例。而且，描述的具体特征、结构、材料或者特点可以在任一个或多个实施例或示例中以合适的方式结合。此外，在不相互矛盾的情况下，本领域的技术人员可以将本说明书中描述的不同实施例或示例以及不同实施例或示例的特征进行结合和组合。
[0124]
尽管上面已经示出和描述了本发明的实施例，可以理解的是，上述实施例是示例性的，不能理解为对本发明的限制，本领域的普通技术人员在本发明的范围内可以对上述实施例进行变化、修改、替换和变型。