慢性阻塞性肺病标志微生物及其应用

摘要

本发明提出了慢性阻塞性肺病标志微生物及其应用，该慢性阻塞性肺病标志微生物包括第一微生物集，因此，进一步提出了一种试剂盒，包括适于检测第一微生物集中的至少一种菌种的试剂，所述第一微生物集由以下菌种组成：Parabacteroides_sp_CAG_2、Porphyromonas_sp_31_2、Prevotella_sp_CAG_617、Prevotella_sp_CAG_755。本发明提出的微生物在健康人群和慢性阻塞性肺病患者中的丰度具有显著差异，可以作为有效检测和/或治疗慢性阻塞性肺病的标志物。

说明书

技术领域

本发明涉及生物技术领域，具体的，本发明涉及慢性阻塞性肺病标志微生物及其用途，更具体的，本申请涉及一种试剂盒、试剂在制备试剂盒中的用途、用于预防或者治疗慢性阻塞性肺病的药物组合物或者食品组合物、确定个体是否患有慢性阻塞性肺病方法、确定个体是否患有慢性阻塞性肺病的装置、一种装置、一种筛选药物的方法。

背景技术

慢性阻塞性肺病(Chronic obstructive pulmonary dise COPDe，COPD)是一种异质性疾病，具有肺部病变，包括慢性支气管炎、气道重塑和损害肺功能的肺气肿。它有许多全身性合并症，例如心血管疾病、结肠炎和骨质疏松症。该疾病产生的主要原因是吸入香烟烟雾、空气污染或其他有害颗粒，同时，一些遗传因素也会导致该病的发生，包括遗传因素、气道反应性增高、在怀孕期、新生儿期、婴儿期或儿童期由于各种原因导致肺发育或生长不良等；另外，炎症和免疫反应等因素也会导致该疾病的发生。

当前的COPD治疗方法有限，主要治疗目的是控制症状和减少恶化。其中，高剂量吸入皮质类固醇被广泛使用，但其功效仅限于减少急性发作频率，或与支气管扩张剂联合使用以改善COPD症状。然而，许多患者对类固醇治疗没有反应，而且这些疗法无法改变引发和推动疾病进展的因素，不能逆转组织病变，并且容易发生严重的呼吸道感染和肺炎。研究发现，呼吸道中的病毒和细菌与疾病恶化有关。

随着人体基因组测序完成及高通量测序技术的高速发展，基因筛查成为慢性阻塞性肺病诊断的方向，对于发现慢性阻塞性肺病潜在人群很有优势。到目前为止，还没有慢性阻塞性肺病患者肠道微生物标记物的相关报道。因此，研究慢性阻塞性肺病患者肠道微生物群的特征对慢性阻塞性肺病的检测及治疗具有重要意义。

发明内容

本申请是基于发明人对以下事实和问题的发现和认识作出的：

本申请的申请人通过前期大量研究，意外地发现一种慢性阻塞性肺病标志微生物，为早期发现慢性阻塞性肺病提供一种非侵入性方法；合理有效地应用标志微生物，扶持肠道有益菌的生长，抑制肠道潜在致病菌，可以治疗或减轻慢性阻塞性肺病的临床症状。

为此，在本发明的第一方面，本发明提出了一种试剂盒。根据本发明的实施例，所述试剂盒包括适于检测第一微生物集中的至少一种菌种的试剂，所述第一微生物集由以下菌种组成：Parabacteroides_sp_CAG_2、Porphyromonas_sp_31_2、Prevotella_sp_CAG_617、Prevotella_sp_CAG_755。所述标志微生物是发明人对大量已知状态的粪便样本进行验证，通过差异比较分析各种肠道微生物在慢性阻塞性肺病组和健康组个体粪便样本中的丰度，而确定下来的，可以有效区分慢性阻塞性肺病组和健康组个体。因此，根据本发明具体实施例的试剂盒，可以准确的检测第一微生物集中的至少一种菌种，进行学研究或准确区分或诊断慢性阻塞性肺病患者和健康个体。

在本发明的第二方面，本发明提出了试剂在制备试剂盒中的用途，所述试剂适于检测第一微生物集中的至少一种菌种。根据本发明的实施例，所述试剂盒用于诊断慢性阻塞性肺病或检测慢性阻塞性肺病的治疗效果，所述第一微生物集由以下菌种组成：Parabacteroides_sp_CAG_2、Porphyromonas_sp_31_2、 Prevotella_sp_CAG_617、Prevotella_sp_CAG_755。根据本发明具体实施例的试剂制备的试剂盒，可以准确的检测所述第一微生物集中的至少一种菌种，极准确的区分慢性阻塞性肺病患者和健康个体，由此，可以有效的在早期进行慢性阻塞性肺病诊断，或用于检测治疗过程中慢性阻塞性肺病的变化。

在本发明的第三方面，本发明提出了一种用于预防或者治疗慢性阻塞性肺病的药物组合物或者食品组合物。根据本发明的实施例，含有第一微生物集中的至少一种菌种，所述第一微生物集由以下菌种组成：Parabacteroides_sp_CAG_2、Porphyromonas_sp_31_2、Prevotella_sp_CAG_617、 Prevotella_sp_CAG_755。根据本发明实施例的慢性阻塞性肺病标志微生物中第一微生物集的菌种可以非侵入性的在早期发现或辅助检测慢性阻塞性肺病，确定个体患有慢性阻塞性肺病的概率高低或者个体处于健康状态的概率高低；同时，提高慢性阻塞性肺病高风险人群或慢性阻塞性肺病患者肠道内的所述第一微生物集中的各种菌种，可以降低患慢性阻塞性肺病的概率或减缓、治愈慢性阻塞性肺病，因此，所述包含所述第一微生物集中的至少一种菌种的药物或者食品组合物能够用于平衡肠道菌群，有效预防或治疗慢性阻塞性肺病。

在本发明的第四方面，本发明提出了一种确定个体是否患有慢性阻塞性肺病的方法。根据本发明的实施例，包括：(1)确定所述个体的粪便样本中标志微生物的丰度，所述标志微生物包括第一微生物集和第二微生物集中的至少一种菌种；(2)将步骤(1)中得到的所述丰度与预定的阈值进行比较，以便确定所述个体是否患有慢性阻塞性肺病；其中，所述第一微生物集由以下菌种组成： Parabacteroides_sp_CAG_2、Porphyromonas_sp_31_2、Prevotella_sp_CAG_617、Prevotella_sp_CAG_755；所述第二微生物集由以下菌种组成：Sellimonas intestinalis、Lachnospiraceae bacterium_COE1、 Streptococcusparasanguinis、Streptococcus tigurinus、Streptococcus downei、Clostridiumsaccharolyticum、 Bifidobacterium aesculapii。所述标志微生物是发明人对大量已知状态的粪便样本进行验证，通过差异比较分析各种肠道微生物在慢性阻塞性肺病组和健康组粪便样本中的丰度，而确定下来的。根据本发明实施例的方法可以依据个体的粪便样本中的所述标志微生物的丰度确定个体是否患有慢性阻塞性肺病。

在本发明的第五方面，本发明提出了一种确定个体是否患有慢性阻塞性肺病的装置。根据本发明的实施例，包括：丰度确定单元，用于确定所述个体的粪便样本中标志微生物的丰度，所述标志微生物包括第一微生物集和第二微生物集中的至少一种菌种；比较单元，用于将所得到的所述丰度与预定的阈值进行比较，以便确定所述个体是否患有慢性阻塞性肺病；其中，所述第一微生物集由以下菌种组成： Parabacteroides_sp_CAG_2、Porphyromonas_sp_31_2、Prevotella_sp_CAG_617、Prevotella_sp_CAG_755；所述第二微生物集由以下菌种组成：Sellimonas intestinalis、Lachnospiraceae bacterium_COE1、Streptococcus parasanguinis、Streptococcus tigurinus、Streptococcus downei、Clostridium saccharolyticum、 Bifidobacterium aesculapii。所述标志微生物是发明人通过差异比较分析各种肠道微生物在慢性阻塞性肺病患者和健康人群的粪便样本中的丰度并经过分析和大量已知状态的粪便样本的验证而确定下来的，根据本发明实施例的装置可以准确确定个体是否为慢性阻塞性肺病的高风险人群或慢性阻塞性肺病患者。

在本发明的第六方面，本发明提出了一种装置。根据本发明的实施例，包括：计算机可读存储介质，其上存储有计算机程序，所述程序用于执行第四方面所述的方法；以及一个或者多个处理器，用于执行所述计算机可读存储介质中的程序。根据本发明实施例的装置可以准确确定个体是否为慢性阻塞性肺病的高风险人群或慢性阻塞性肺病患者。

在本发明的第七方面，本发明提出了一种筛选药物的方法。根据本发明的实施例，所述药物用于治疗或者预防慢性阻塞性肺病，所述方法包括：将候选药物施用于受试者，检测施用前后，所述受试者粪便中标志微生物的丰度，所述标志微生物包括第一微生物集和第二微生物集中的至少一种菌种，其中，满足下列条件至少之一的候选药物适于用于治疗或者预防慢性阻塞性肺病：(1)进行所述施用后，所述第一微生物集中的至少一种菌种的所述丰度升高；和(2)进行所述施用后，所述第二微生物集中的至少一种菌种的所述丰度降低；其中，所述第一微生物集由以下菌种组成：Parabacteroides_sp_CAG_2、Porphyromonas_sp_31_2、Prevotella_sp_CAG_617、Prevotella_sp_CAG_755；所述第二微生物集由以下菌种组成：Sellimonas intestinalis、Lachnospiraceae bacterium_COE1、Streptococcus parasanguinis、 Streptococcus tigurinus、Streptococcus downei、Clostridium saccharolyticum、Bifidobacterium aesculapii。根据本发明实施例的方法可以生产或筛选出促进所述标志微生物中第一微生物集中各种菌种生长，和/ 或抑制肠道标志微生物中第二微生物集中的各种菌种生长的药物，对于辅助减轻慢性阻塞性肺病的临床症状具有重要意义。

本发明的附加方面和优点将在下面的描述中部分给出，部分将从下面的描述中变得明显，或通过本发明的实践了解到。

附图说明

本发明的上述和/或附加的方面和优点从结合下面附图对实施例的描述中将变得明显和容易理解，其中：

图1是根据本发明实施例的筛选慢性阻塞性肺病标志微生物的实验分析流程示意图；以及

图2是根据本发明实施例的标志微生物综合指标AUC评价结果示意图，其中，Specificity表示特异度，即预测为阳性且实际为阳性，真阳性，纵坐标Sensitivity表示敏感度，即真阴性：

2-A为第一期37个样品数据ROC曲线下AUC值和置信区间结果图；

2-B为第二期20个样品数据ROC曲线下AUC值和置信区间结果图。

具体实施方式

下面详细描述本发明的实施例，所述实施例的示例在附图中示出。下面通过参考附图描述的实施例是示例性的，旨在用于解释本发明，而不能理解为对本发明的限制。

此外，术语“第一”、“第二”仅用于描述目的，而不能理解为指示或暗示相对重要性或者隐含指明所指示的技术特征的数量。由此，限定有“第一”、“第二”的特征可以明示或者隐含地包括至少一个该特征。在本发明的描述中，“多个”的含义是至少两个，例如两个，三个等，除非另有明确具体的限定。

术语“任选地”仅用于描述目的，而不能理解为指示或暗示相对重要性。由此，限定有“任选地”的特征可以明示或者隐含地包括或不包括该特征。

生物学标志物是从生物学介质中可以检测到的细胞/生物化学或分子改变。生物学介质包括各种体液、组织、细胞、粪便、头发、呼气等。

微生物的丰度是指在某一微生物群体中该种微生物的丰富程度，例如在肠道微生物群体中的该种微生物程度，可表示为该种微生物在该群体中的含量。

根据本发明提供的一种试剂盒，包括适于检测第一微生物集中的至少一种菌种的试剂，所述第一微生物集由以下菌种组成：Parabacteroides_sp_CAG_2、Porphyromonas_sp_31_2、Prevotella_sp_CAG_617、Prevotella_sp_CAG_755。

根据本发明的具体实施方案，所述试剂盒进一步包括适于检测第二微生物集中的至少一种菌种的试剂，所述第二微生物集由以下菌种组成：Sellimonas intestinalis、Lachnospiraceae bacterium_COE1、 Streptococcus parasanguinis、Streptococcustigurinus、Streptococcus downei、Clostridium saccharolyticum、 Bifidobacteriumaesculapii。

根据本发明的具体实施方案，所述试剂盒包括适于检测所述第一微生物集中全部所述菌种的试剂。

根据本发明的具体实施方案，所述试剂盒包括适于检测所述第二微生物集中全部所述菌种的试剂。

根据本发明的具体实施方案，所述标志微生物是发明人通过对大量患慢性阻塞性肺病个体和大量健康对照个体的粪便样本中的微生物的丰度的差异比较分析、以及验证，而确定下来的，明确了肠道微生物中慢性阻塞性肺病相关的标志微生物。利用包含检测所述标志微生物的试剂的试剂盒能够确定个体处于患有慢性阻塞性肺病状态的概率高低或者处于健康状态的概率高低，能够用于非侵入性的早期发现或辅助检测慢性阻塞性肺病。

根据本发明的具体实施方案，所述适于检测所述第一微生物集或第二微生物集的试剂不受特别限制，任何可以检测所述微生物菌种的试剂均包含在本发明的范围内，如通过形态学特征、生理生化反应特征、生态学特征以及血清学反应、对噬菌体的敏感性、分子生物学等方面检测所述微生物菌种的试剂，具体的，如抗体、酶、核酸分子等。

本文中，所述微生物形态学特征指：在显微镜下观察微生物外形大小、形状、排列等，细胞构造，革兰氏染色反应，能否运动、鞭毛着生部位和数目，有无芽孢和荚膜、芽孢的大小和位置，放线菌和真菌的繁殖器官的形状、构造，孢子的数目、形状、大小、颜色和表面特征等。

本文中，所述微生物生理生化反应特征指：微生物利用物质的能力、代谢产物的特殊性，如是否产生H

本文中，所述微生物血清学反应指：利用抗原与抗体的高度敏感特异性反应，来鉴别相似的菌种，或对同种微生物分型，如用已知菌种、型或菌株制成的抗血清，与待鉴定的所述微生物是否发生特异性的血清学反应来鉴定微生物。

本文中，分子生物学方法检测微生物主要包括：利用PCR技术、高通量测序等方法。

根据本发明提供的试剂在制备试剂盒中的用途，所述试剂适于检测第一微生物集中的至少一种菌种，所述试剂盒用于诊断慢性阻塞性肺病或者检测慢性阻塞性肺病的治疗效果，所述第一微生物集由以下菌种组成：Parabacteroides_sp_CAG_2、Porphyromonas_sp_31_2、Prevotella_sp_CAG_617、 Prevotella_sp_CAG_755。

根据本发明的具体实施例，所述标志微生物是发明人通过对大量患慢性阻塞性肺病个体和大量健康对照个体的粪便样本中的微生物的丰度的差异比较分析、以及验证，而确定下来的，明确了肠道微生物中慢性阻塞性肺病相关的微生物标志物。利用检测所述标志微生物的试剂能够确定个体患有慢性阻塞性肺病的概率高低或者处于健康状态的概率高低，能够用于非侵入性的早期发现或辅助检测慢性阻塞性肺病。

根据本发明一些具体的实施例，所述试剂进一步适于检测第二微生物集中的至少一种菌种，所述第二微生物集由以下菌种组成：Sellimonas intestinalis、Lachnospiraceae bacterium_COE1、Streptococcus parasanguinis、Streptococcustigurinus、Streptococcus downei、Clostridium saccharolyticum、Bifidobacteriumaesculapii。

根据本发明的具体实施例，所述适于检测所述第一微生物集或第二微生物集的试剂不受特别限制，可以检测所述微生物菌种的试剂均包含在本发明的范围内，如通过形态学特征、生理生化反应特征、生态学特征以及血清学反应、对噬菌体的敏感性、分子生物学等方面检测所述微生物菌种的试剂等，具体的，如抗体、酶、核酸分子。

本文中，所述微生物形态学特征指：在显微镜下观察微生物外形大小、形状、排列等，细胞构造，革兰氏染色反应，能否运动、鞭毛着生部位和数目，有无芽孢和荚膜、芽孢的大小和位置，放线菌和真菌的繁殖器官的形状、构造，孢子的数目、形状、大小、颜色和表面特征等。

本文中，所述微生物生理生化反应特征指：微生物利用物质的能力、代谢产物的特殊性，如是否产生H

本文中，所述微生物血清学反应指：利用抗原与抗体的高度敏感特异性反应，来鉴别相似的菌种，或对同种微生物分型，如用已知菌种、型或菌株制成的抗血清，与待鉴定的所述微生物是否发生特异性的血清学反应来鉴定微生物。

本文中，分子生物学方法检测微生物主要包括：利用PCR技术、高通量测序等方法。

根据本发明提供的一种用于预防或者治疗慢性阻塞性肺病的药物组合物或者食品组合物，含有第一微生物集中的至少一种菌种，所述第一微生物集由以下菌种组成：Parabacteroides_sp_CAG_2、 Porphyromonas_sp_31_2、Prevotella_sp_CAG_617、Prevotella_sp_CAG_755。

上述第一微生物集中的菌种是发明人通过差异比较分析各种肠道微生物在慢性阻塞性肺病疾病组和健康组的粪便样本中的丰度，以及经过大量已知状态的粪便样本的验证，而确定下来的。所述标志微生物中的第一微生物集菌种相较于慢性阻塞性肺病患者群体，在健康群体组中显著富集，所述显著富集是指与在慢性阻塞性肺病患者组中的丰度相比，上述菌种在健康组中的丰度均具有统计意义地高于或者明显地、实质性地高于在慢性阻塞性肺病患者组中的丰度；能够使该部分菌种丰度提高的物质能够用于治疗慢性阻塞性肺病或者益于慢性阻塞性肺病患者服用，能够使其丰度提高的物质不限于治疗慢性阻塞性肺病的药物和有益肠道菌群平衡的功能性食品。该实施方式提供的标志微生物能够用于制备治疗慢性阻塞性肺病的药物和/或用于制备益于平衡肠道菌群的功能性食品、保健药等，所述药物或食品可有效治疗或缓解慢性阻塞性肺病。

根据本发明提供的一种确定个体是否患有慢性阻塞性肺病的方法，包括步骤(1)和(2)。

所述标志微生物包括第一微生物集和第二微生物集中的至少一种菌种。其中，所述第一微生物集由以下菌种组成：Parabacteroides_sp_CAG_2、Porphyromonas_sp_31_2、Prevotella_sp_CAG_617、 Prevotella_sp_CAG_755；所述第二微生物集由以下菌种组成：Sellimonas intestinalis、Lachnospiraceae bacterium_COE1、Streptococcusparasanguinis、Streptococcus tigurinus、Streptococcus downei、Clostridiumsaccharolyticum、Bifidobacterium aesculapii。

根据本发明的一些具体的实施例，步骤(1)进一步包括：获得所述个体的粪便样本中的核酸测序数据；将所述测序数据与参考基因组进行比对；基于所述比对的结果，确定所述标志微生物的丰度。

根据本发明的具体实施例，在步骤(1)中，按照下列公式确定所述标志微生物的丰度： Ab(S)＝Ab(U

比对可以利用已知比对软件进行，例如SOAP、BWA和TeraMap等，在比对过程中，一般对比对参数进行设置，设置一个或者一对读段(reads)最多允许有s个碱基错配(mismatch)，例如设置s≤2，若reads中有超过s个碱基发生错配，则视为该reads无法比对到(比对上)该组装片段上。所述的获得的比对结果包含各条读段与各物种的参考基因组的比对情况，包括读段是否能够比对上某个或某些物种的参考基因组、只唯一比对到一种物种还是比对到多种物种的参考基因组、比对到物种的参考基因组的位置、比对到物种参考基因组的唯一位置还是多个位置等信息。

所述菌种/微生物的参考基因组指预先确定的该微生物物种的序列，可以是预先获得的待测样本所属或者所包含的生物类别的任意参考模板，例如，目标是待测样本中的微生物，参考序列可选择NCBI数据库中的各种微生物的参考基因组和/或HMP、MetaHIT项目公开的DACC肠道参考基因组，进一步地，也可以预先配置包含更多参考序列的资源库，例如依据待测样本来源的个体的状态、地域等因素选择或是测定组装出更接近的序列作为参考序列。根据本发明的一个实施例，各种微生物的参考基因组从公开的数据库中获得，通常的，一种微生物的参考基因组有多个版本，即一种微生物有多个公开的参考基因组。

reads与物种的参考基因组比对，比对上的可以被分为两部分：a)Unique reads(U)：唯一比对上一个物种的参考基因组；称这些reads为unique reads。即，如果reads比对上的参考基因组均来自同一物种，定义这些reads为unique reads；b)Multiple reads(M)：比对上一个以上物种的参考基因组，定义为multiple reads。即，如果reads比对上的参考基因组来自至少两种物种，定义这些reads为multiple reads。

根据本发明的一个实施例，将步骤(1)中得到的所述丰度与预定的阈值进行比较，以便确定所述个体是否患有慢性阻塞性肺病。

根据本发明的一些具体实施例，所述阈值为预先设定的。将标志微生物在待测个体样本中的丰度与所述的阈值比较，进行待测个体的状态的确定。所述阈值可以为一数值或者数值范围，阈值设定方式不受特别限制，能够对所述个体的是否患慢性阻塞性肺病进行判断的阈值均可使用，例如，基于已知患病或健康状态个体中的标志微生物的丰度均值，该微生物对应的阈值可以设为该丰度均值的95％的置信区间。

所述的置信区间是指由样本统计量所构造的总体参数的估计区间。在统计学中，一个概率样本的置信区间(Confidence interval)是对这个样本的某个总体参数的区间估计。置信区间展现的是这个参数的真实值有一定概率落在测量结果的周围的程度。置信区间给出的是被测量参数的测量值的可信程度，即前面所要求的“一定概率”，这个概率被称为置信水平。

根据本发明的一些具体实施例，当步骤(1)中确定的标志微生物的丰度达到所述患慢性阻塞性肺病丰度阈值，未达到所述不患慢性阻塞性肺病丰度阈值时，确定所述个体患有慢性阻塞性肺病，当(1)中确定的标志微生物的丰度达到不患慢性阻塞性肺病丰度阈值，未达到患慢性阻塞性肺病丰度阈值时，确定所述个体不患慢性阻塞性肺病。

需要说明的是，根据目的或要求不同，可能对确定个体状态结果的可信程度有不同的要求，本领域技术人员可以选择不同的显著性水平或阈值。

该方法基于检测个体的粪便样本中的标志微生物中的各种菌种的丰度，分别将检测确定的各种菌种的丰度与其阈值进行比较，依据获得的比较结果能够确定个体为慢性阻塞性肺病个体或者为健康个体的概率。为早期发现慢性阻塞性肺病提供一种非侵入性的辅助检测或者辅助干预治疗的方法。

以上任一实施例中的利用标志微生物确定个体是否患有慢性阻塞性肺病的方法的全部或部分步骤，可以利用包含可拆分的相应单元功能模块的装置/系统来施行，或者将方法程序化、存储于机器可读介质，利用机器运行该可读介质来实现。

根据本发明提供的一种确定个体是否患有慢性阻塞性肺病的装置，该装置包括：丰度确定单元，用于确定所述个体的粪便样本中标志微生物的丰度；比较单元，用于将所得到的所述丰度与预定的阈值进行比较，以便确定所述个体是否患有慢性阻塞性肺病；其中，所述第一微生物集由以下菌种组成： Parabacteroides_sp_CAG_2、Porphyromonas_sp_31_2、Prevotella_sp_CAG_617、Prevotella_sp_CAG_755；所述第二微生物集由以下菌种组成：Sellimonas intestinalis、Lachnospiraceae bacterium_COE1、 Streptococcusparasanguinis、Streptococcus tigurinus、Streptococcus downei、Clostridiumsaccharolyticum、 Bifidobacterium aesculapii。上述对本发明任一实施例的利用标志微生物确定个体是否患有慢性阻塞性肺病的方法的技术特征和优点的描述，同样适用本发明这一方面的装置，在此不再赘述。

根据本发明的实施例，所述丰度确定单元适于通过下列步骤确定所述丰度：获得所述个体的粪便样本中的核酸测序数据；将所述测序数据与参考基因组进行比对；基于所述比对的结果，确定所述标志微生物的丰度。

比对可以利用已知比对软件进行，例如SOAP、BWA和TeraMap等，在比对过程中，一般对比对参数进行设置，设置一个或者一对读段(reads)最多允许有s个碱基错配(mismatch)，例如设置s≤2，若reads中有超过s个碱基发生错配，则视为该reads无法比对到(比对上)该组装片段上。所述的获得的比对结果包含各条读段与各物种的参考基因组的比对情况，包括读段是否能够比对上某个或某些物种的参考基因组、只唯一比对到一种物种还是比对到多种物种的参考基因组、比对到物种的参考基因组的位置、比对到物种参考基因组的唯一位置还是多个位置等信息。根据本发明的一个实施例，利用 SOAPalign 2.21进行比对，设置参数为–r 2–m 100–x 1000。

所述微生物的参考基因组指预先确定的该微生物物种的序列，可以是预先获得的待测样本所属或者所包含的生物类别的任意参考模板，例如，目标是待测样本中的微生物，参考序列可选择NCBI数据库中的各种微生物的参考基因组和/或HMP、MetaHIT项目公开的DACC肠道基因组，进一步地，也可以预先配置包含更多参考序列的资源库，例如依据待测样本来源的个体的状态、地域等因素选择或是测定组装出更接近的序列作为参考序列。根据本发明的一个实施例，各种微生物的参考基因组从公开的数据库中获得，通常的，一种微生物的参考基因组有多个版本，即一种微生物有多个公开的参考基因组。

reads与物种的参考基因组比对，比对上的可以被分为两部分：a)Unique reads(U)：唯一比对上一个物种的基因组；称这些reads为unique reads。即，如果reads比对上的参考基因组均来自同一物种，定义这些reads为unique reads。b)Multiple reads(M)：比对上一个以上物种的参考基因组，定义为multiple reads。即，如果reads比对上的参考基因组来自至少两种物种，定义这些reads为multiple reads。

根据本发明的一个实施例，按照下列公式确定所述标志微生物的丰度： Ab(S)＝Ab(U

根据本发明的又一个实施方式提供的一种装置，包括：计算机可读存储介质，其上存储有计算机程序，所述程序用于执行前面所述的一种确定个体是否患有慢性阻塞性肺病方法；以及一个或者多个处理器，用于执行所述计算机可读存储介质中的程序。

根据本发明的又一个实施方式提供的一种筛选药物的方法，所述药物用于治疗或者预防慢性阻塞性肺病，所述方法包括：将候选药物施用于受试者，检测施用前后，所述受试者粪便中标志微生物的丰度，所述标志微生物包括第一微生物集和第二微生物集中的至少一种菌种，其中，满足下列条件至少之一的候选药物适于用于治疗或者预防慢性阻塞性肺病：(1)进行所述施用后，所述第一微生物集中的至少一种菌种的所述丰度升高；和(2)进行所述施用后，所述第二微生物集中的至少一种菌种的所述丰度降低；其中，所述第一微生物集由以下菌种组成：Parabacteroides_sp_CAG_2、Porphyromonas_sp_31_2、Prevotella_sp_CAG_617、Prevotella_sp_CAG_755；所述第二微生物集由以下菌种组成：Sellimonas intestinalis、Lachnospiraceae bacterium_COE1、Streptococcusparasanguinis、Streptococcus tigurinus、 Streptococcus downei、Clostridiumsaccharolyticum、Bifidobacterium aesculapii。

利用本发明这一方面的生产或筛选治疗慢性阻塞性肺病的药物的方法，通过合理有效地应用确定的慢性阻塞性肺病生物标志物进行筛选，能够获得能扶持肠道有益菌的生长和/或抑制肠道潜在致病菌的药物，对于辅助减轻慢性阻塞性肺病的临床症状具有重要意义。

下面将对实施例作具体介绍。以下实施例中涉及的未特别交待的试剂、序列、软件及仪器，都是常规市售产品。

实施例1标志微生物的鉴定

本实施例中，发明人共获得28个慢性阻塞性肺病患者和29个健康对照的粪便样品，通过对其中18 个慢性阻塞性肺病患者和19个健康对照的粪便样品进行研究，从而获得肠道菌群微生物微生物群落及功能成分特征。总的来说，发明人下载高质量的测序数据(control健康人)以及实验测序得到的高质量慢性阻塞性肺病测序数据。宏基因组分析显示，21个微生物物种与慢性阻塞性肺病疾病密切相关，其中 9个菌在健康人的肠道微生物中富集，12个菌在慢性阻塞性肺病病人的肠道微生物中富集。

1、测序数据的获取：

1.1参照数据下载

从EBI下载28个疾病组数据和29个健康的测序数据，下载地址为PRJNA562766。

参照图1所示的实验流程，鉴定慢性阻塞性肺病的相关生物标志物，其中省略的步骤或者细节为本领域技术人员所熟知，几个重要步骤介绍如下面所述。

2、微生物物种丰度分析

2.1序列优化统计

1)首先进行第一期测序，本期测序采集所述28个患病个体的数据后，对其进行过滤，质控按以下标准进行：a)移除大于5个N碱基的reads；b)去除低质量碱基(Q20)大于50％的reads；c)移除尾部低质量(Q20)和N碱基。丢失成对的reads序列被认为是单条reads用于组装。

2)采用1)所述方法同样处理下载的健康人数据。

2.2物种丰度分析

SOAPalign 2.21用于匹配针对冗余基因组的paired-end clean reads，这里，所称的冗余基因组来自各数据库中公开的的细菌的参考基因组，比对参数为–r 2–m 200–x1000。Reads与冗余基因组的比对结果，可分为两部分：a)Unique reads(U)：reads只比对上一个物种的基因组；这些reads被定义为unique reads。即，如果这些基因组来自同一物种，发明人将这些reads定义为unique reads。b)Multiple reads(M)：如果reads比对上两个及两个以上物种的基因组，定义为multiple reads。即，如果比对上的基因组来自不同物种，发明人定义这些reads为multiple reads。

对于物种S，其丰度为Ab(S)，与特有的U reads和共享的M reads相关，丰度的计算方式如下：

Ab(S)＝Ab(U

其中，S表示所述标志微生物的编号，

Ab(S)表示所述标志微生物S的丰度，

Ab(U

U

l

M

Co

i表示所述非唯一比对读段的编号，

Co

所述非唯一比对的读段i的丰度系数，

N为所述非唯一比对的读段i能够比对的微生物的总数，

j表示所述非唯一比对的读段i能够比对的微生物的编号。

每个样品中计算得到的物种丰度值，全部除以每个样品的丰度总数后，得到归一化后的物种丰度表。

3.3筛选微生物物种标记物

为了获得与慢性阻塞性肺病疾病密切相关的肠道微生物物种标记物，发明人利用慢性阻塞性肺病患者(COPD)组(18例)与正常人(H)组(19例)两组肠道微生物物种丰度数据，在物种级别做了一个与疾病相关性的研究。

基于步骤2.2所得到的物种丰度表，发明人设置标准如下：(1)慢性阻塞性肺病患者组或健康人组物种丰度的中位数必须大于0.00001；(2)通过结合Benjamini Hochberg的多重检验的 Wilcoxon秩和检验进行检验，得到每个物种和慢性阻塞性肺病疾病的的相关性p值和fdr值；(3)使用一个严格的阈值p_values<0.05，利用上述参数进行筛选。发明人得到21个与慢性阻塞性肺病疾病密切相关的肠道微生物物种，其中，在慢性阻塞性肺病患者肠道内中富集的微生物有12个物种，在健康人中富集的有9个物种，这21个微生物物种标记物如表1所示。

表1：

实施例2微生物物种标记物的验证

为了验证实施例1中的发现，发明人进一步分析验证群体中的10个健康人及10个慢性阻塞性肺病患者的粪便样本中的所述21种菌属的丰度，并根据各物种在健康组及疾病组的富集情况对所述21个微生物物种标记物做出删选，验证群体的DNA提取、测序以及物种丰度的分析已在进行实施例1执行，可以理解为此处的实验操作同实施例1。

验证结果如下：上述富集在健康人群中的9个物种，4个在验证集中得到高质量的验证(p_values<0.05)，健康人富集的微生物物种标志物验证的均值和P值结果如表2所示。

表2：

对于上述富集在慢性阻塞性肺病患者中的12个物种中，其中7个在验证集中得到高质量的验证(p values<0.05)，慢性阻塞性肺病患者中的富集的微生物标志物验证的p值和q值情况如表3所示。

表3：

发明人认为，可以将从健康人富集的4个微生物物种标记物，作为慢性阻塞性肺病疾病患病的反向指标，或作为治疗慢性阻塞性肺病进行研发的微生物制剂药物菌群成分，或作为检测慢性阻塞性肺病、监测慢性阻塞性肺病治疗进程的恢复指标；将慢性阻塞性肺病患者富集的个微生物物种标记物，作为慢性阻塞性肺病疾病患病的正向指标，特别用于慢性阻塞性肺病疾病非创伤式的检测和诊断。

发明人利用这11个微生物物种标记物，构建一个综合指标，估计ROC(Receiver-operating characteristic)曲线下面积AUC，AUC越大，表示诊断能力越高，评价综合得分对应其对慢性阻塞性肺病的诊断能力。通过对一期(第一阶段)的37个样品和二期(第二阶段)的20个样品进行评测，具体情况如图2所示，都表现出了很好的诊断能力，在一期得到AUC＝88.9％，如图2-A所示，置信区间为58.1％-100％；二期得到AUC＝85.7％，如图2-B所示，置信区间为63.1％-100％。

在本说明书的描述中，参考术语“一个实施例”、“一些实施例”、“示例”、“具体示例”、或“一些示例”等的描述意指结合该实施例或示例描述的具体特征、结构、材料或者特点包含于本发明的至少一个实施例或示例中。在本说明书中，对上述术语的示意性表述不必须针对的是相同的实施例或示例。而且，描述的具体特征、结构、材料或者特点可以在任一个或多个实施例或示例中以合适的方式结合。此外，在不相互矛盾的情况下，本领域的技术人员可以将本说明书中描述的不同实施例或示例以及不同实施例或示例的特征进行结合和组合。

尽管上面已经示出和描述了本发明的实施例，可以理解的是，上述实施例是示例性的，不能理解为对本发明的限制，本领域的普通技术人员在本发明的范围内可以对上述实施例进行变化、修改、替换和变型。