CN105733988B - 组合物及应用 - Google Patents

组合物及应用 Download PDF

Info

Publication number
CN105733988B
CN105733988B CN201610165123.XA CN201610165123A CN105733988B CN 105733988 B CN105733988 B CN 105733988B CN 201610165123 A CN201610165123 A CN 201610165123A CN 105733988 B CN105733988 B CN 105733988B
Authority
CN
China
Prior art keywords
species
abundance
ankylosing spondylitis
composition
reads
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN201610165123.XA
Other languages
English (en)
Other versions
CN105733988A (zh
Inventor
郑智俊
张林爽
吴春燕
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Shanghai Realbio Technology Co ltd
Original Assignee
Shanghai Realbio Technology Co ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Shanghai Realbio Technology Co ltd filed Critical Shanghai Realbio Technology Co ltd
Priority to CN201610165123.XA priority Critical patent/CN105733988B/zh
Publication of CN105733988A publication Critical patent/CN105733988A/zh
Application granted granted Critical
Publication of CN105733988B publication Critical patent/CN105733988B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • C—CHEMISTRY; METALLURGY
    • C12—BIOCHEMISTRY; BEER; SPIRITS; WINE; VINEGAR; MICROBIOLOGY; ENZYMOLOGY; MUTATION OR GENETIC ENGINEERING
    • C12N—MICROORGANISMS OR ENZYMES; COMPOSITIONS THEREOF; PROPAGATING, PRESERVING, OR MAINTAINING MICROORGANISMS; MUTATION OR GENETIC ENGINEERING; CULTURE MEDIA
    • C12N1/00—Microorganisms, e.g. protozoa; Compositions thereof; Processes of propagating, maintaining or preserving microorganisms or compositions thereof; Processes of preparing or isolating a composition containing a microorganism; Culture media therefor
    • C12N1/20—Bacteria; Culture media therefor
    • C—CHEMISTRY; METALLURGY
    • C12—BIOCHEMISTRY; BEER; SPIRITS; WINE; VINEGAR; MICROBIOLOGY; ENZYMOLOGY; MUTATION OR GENETIC ENGINEERING
    • C12Q—MEASURING OR TESTING PROCESSES INVOLVING ENZYMES, NUCLEIC ACIDS OR MICROORGANISMS; COMPOSITIONS OR TEST PAPERS THEREFOR; PROCESSES OF PREPARING SUCH COMPOSITIONS; CONDITION-RESPONSIVE CONTROL IN MICROBIOLOGICAL OR ENZYMOLOGICAL PROCESSES
    • C12Q1/00—Measuring or testing processes involving enzymes, nucleic acids or microorganisms; Compositions therefor; Processes of preparing such compositions
    • C12Q1/02—Measuring or testing processes involving enzymes, nucleic acids or microorganisms; Compositions therefor; Processes of preparing such compositions involving viable microorganisms
    • C12Q1/025—Measuring or testing processes involving enzymes, nucleic acids or microorganisms; Compositions therefor; Processes of preparing such compositions involving viable microorganisms for testing or evaluating the effect of chemical or biological compounds, e.g. drugs, cosmetics
    • C—CHEMISTRY; METALLURGY
    • C12—BIOCHEMISTRY; BEER; SPIRITS; WINE; VINEGAR; MICROBIOLOGY; ENZYMOLOGY; MUTATION OR GENETIC ENGINEERING
    • C12Q—MEASURING OR TESTING PROCESSES INVOLVING ENZYMES, NUCLEIC ACIDS OR MICROORGANISMS; COMPOSITIONS OR TEST PAPERS THEREFOR; PROCESSES OF PREPARING SUCH COMPOSITIONS; CONDITION-RESPONSIVE CONTROL IN MICROBIOLOGICAL OR ENZYMOLOGICAL PROCESSES
    • C12Q1/00—Measuring or testing processes involving enzymes, nucleic acids or microorganisms; Compositions therefor; Processes of preparing such compositions
    • C12Q1/02—Measuring or testing processes involving enzymes, nucleic acids or microorganisms; Compositions therefor; Processes of preparing such compositions involving viable microorganisms
    • C12Q1/04—Determining presence or kind of microorganism; Use of selective media for testing antibiotics or bacteriocides; Compositions containing a chemical indicator therefor
    • C12Q1/06—Quantitative determination

Landscapes

  • Chemical & Material Sciences (AREA)
  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Organic Chemistry (AREA)
  • Wood Science & Technology (AREA)
  • Zoology (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • Biotechnology (AREA)
  • Proteomics, Peptides & Aminoacids (AREA)
  • Genetics & Genomics (AREA)
  • General Engineering & Computer Science (AREA)
  • General Health & Medical Sciences (AREA)
  • Microbiology (AREA)
  • Biochemistry (AREA)
  • Toxicology (AREA)
  • Physics & Mathematics (AREA)
  • Analytical Chemistry (AREA)
  • Biophysics (AREA)
  • Medicinal Chemistry (AREA)
  • Immunology (AREA)
  • Molecular Biology (AREA)
  • Biomedical Technology (AREA)
  • Virology (AREA)
  • Tropical Medicine & Parasitology (AREA)
  • Measuring Or Testing Involving Enzymes Or Micro-Organisms (AREA)

Abstract

本发明公开了一种组合物，其包括第一微生物集中的至少两种，第一微生物集由以下菌种组成：毛螺科菌Lachnospiraceae_bacterium_6_1_63FAA、死亡梭杆菌Fusobacterium_mortiferum、拟杆菌Bacteroides_pectinophilus、阴沟肠杆菌Enterobacter_cloacae、挑剔真杆菌Eubacterium_eligens、毛螺科菌Lachnospiraceae_bacterium_8_1_57FAA和梭菌Clostridium_sp._L2‑50。该组合物在健康人群和强直性脊柱炎患者中的丰度具有显著差异，可以作为检测和/或治疗强直性脊柱炎的标志物。

Description

组合物及应用
技术领域
本发明涉及生物标志物领域，具体的，本发明涉及组合物及其应用，更具体的，本发明涉及一种组合物、组合物的用途、一种利用组合物确定个体的状态的方法、一种利用组合物确定个体的状态的装置、一种利用组合物对多个个体进行分类的方法、一种治疗强直性脊柱炎的药物以及一种制备治疗强直性脊柱炎药物的方法。
背景技术
强直性脊柱炎(Ankylosing spondylitis，AS)是一种主要侵犯脊柱，并累及骶髂关节和周围关节的慢性进行性炎性疾病。起病多为15～30岁的男性，现在研究显示，该病男女之比为2:1～3:1，强直性脊柱炎的病因目前尚未完全明确，近年来研究结果表明该病与遗传素质、感染及免疫等因素有关。
强直性脊柱炎一般起病比较隐匿，早期可无任何临床症状，有些病人在早期可表现出轻度的全身症状，如乏力、消瘦、长期或间断低热、厌食、轻度贫血等。由于病情较轻，病人大多不能早期发现，致使病情延误，失去最佳治疗时机。
强直性脊柱炎与人体白细胞抗原HLA-B27密切相关，据流行病学调查，强直性脊柱炎病人HLA-B27阳性率高达90％～96％，而普通人群HLA-27阳性率仅4％～9％；HLA-B27阳性者强直性脊柱炎发病率约为10％～20％，而普通人群发病为1‰～2‰，相差约100倍。虽然HLA-B27检测对强直性脊柱炎的诊断有一定的帮助，但绝大部分的患者还要通过病史、体征和X线检查作出诊断。目前，强直性脊柱炎的早期诊断主要是用CT、放射性核扫描、磁共振等影像学进行检查，检查过程较复杂。
随着人体基因组测序完成及高通量测序技术的高速发展，基因筛查成为强直性脊柱炎诊断的方向，对于发现强直性脊柱炎潜在人群很有优势。研究表明超过70％的强直性脊柱炎患者患有肠道炎症，这些患者中5～10％肠道炎症很严重，可能会发展为临床上的炎症性肠病(IBD)或克罗恩氏病(CD)(Mielants et al.,1985)。克罗恩氏病(CD)的一些标记基因与强直性脊柱炎相关联(Parkes et al.,2013)，说明两种疾病可能有类似发病机制，可能与肠道失调相关。有研究表明与强直性脊柱炎相关的多基因在肠道免疫中起重要作用，如涉及到白介素IL-23途径的基因对肠道健康的调节起重要作用(Wellcome et al.,2007)。与健康对照相比，强直性脊柱炎患者及其直系亲属肠道通透性增加，再次说明肠道微生物在强直性脊柱炎中的重要作用(Mielants et al.,1991)。到目前为止，还没有强直性脊柱炎患者肠道微生物标记物的相关报道。
发明内容
本发明旨在至少在一定程度上解决上述技术问题之一或至少提供一种商业选择。
依据本发明的第一方面，提供一种组合物，其包括第一微生物集中的至少两种菌种，所述第一微生物集由以下菌种组成：毛螺科菌Lachnospiraceae_bacterium_6_1_63FAA、死亡梭杆菌Fusobacterium_mortiferum、拟杆菌Bacteroides_pectinophilus、阴沟肠杆菌Enterobacter_cloacae、挑剔真杆菌Eubacterium_eligens、毛螺科菌Lachnospiraceae_bacterium_8_1_57FAA和梭菌Clostridium_sp._L2-50。
根据本发明的实施例，所述组合物包括所述第一微生物集中的毛螺科菌Lachnospiraceae_bacterium_6_1_63FAA和死亡梭杆菌Fusobacterium_mortiferum；任选的，所述组合物还包括所述第一微生物集中的其它菌种中的至少一种。
根据本发明的实施例，所述组合物包括所述第一微生物集中的毛螺科菌Lachnospiraceae_bacterium_6_1_63FAA和阴沟肠杆菌Enterobacter_cloacae；任选的，所述组合物还包括所述第一微生物集中的其它菌种中的至少一种。
根据本发明的实施例，所述组合物包括所述第一微生物集中的死亡梭杆菌Fusobacterium_mortiferum和阴沟肠杆菌Enterobacter_cloacae；任选的，所述组合物还包括所述第一微生物集中的其它菌种中的至少一种。
根据本发明的实施例，所述组合物包括所述第一微生物集中的阴沟肠杆菌Enterobacter_cloacae和毛螺科菌Lachnospiraceae_bacterium_8_1_57FAA；任选的，所述组合物还包括所述第一微生物集中的其它菌种中的至少一种。
根据本发明的实施例，所述组合物还包括第二微生物集中的至少一种菌种，所述第二微生物集由以下菌种组成：放线菌Actinomyces sp.ICM47、短双歧杆菌Bifidobacterium breve、双歧杆菌Bifidobacterium sp.12_1_47BFAA、长双歧杆菌Bifidobacterium longum、产气柯林斯菌Collinsella aerofaciens、两岐双岐杆菌Bifidobacterium bifidum、假小链双歧杆菌Bifidobacterium pseudocatenulatum、迟缓埃格特菌Eggerthella lenta、链状双歧杆菌Bifidobacterium catenulatum和青春双歧杆菌Bifidobacterium adolescentis。
上述本发明一方面或者任一实施例提供的组合物是发明人通过对大量患强直性脊柱炎个体和大量健康对照个体的粪便样本中的微生物的丰度的差异比较分析、以及验证，而确定下来的，明确了肠道微生物中强直性脊柱炎相关的微生物标志物。利用所称的组合物作为强直性脊柱微生物炎标志物，能够确定个体处于患有强直性脊柱炎状态的概率高低或者处于健康状态的概率高低，能够用于非侵入性的早期发现或辅助检测强直性脊柱炎。
依据本发明的第二方面，本发明提供上述任一实施例中的组合物在检测强直性脊柱炎、监测治疗强直性脊柱炎、制备治疗强直性脊柱炎的药物和/或制备功能性食品中的用途。
上述组合物是发明人通过差异比较分析各种肠道微生物在强直性脊柱炎疾病组和健康组的粪便样本中的丰度，以及经过大量已知状态的粪便样本的验证，而确定下来的。组合物中至少包括一部分菌种，这部分菌种相较于在强直性脊柱炎患者组、在健康群体中显著富集，所称显著富集是指与在疾病组中的丰度相比，上述菌种在健康组中的丰度均具有统计意义地高于或者明显地、实质性地高于在疾病组中的丰度；能够使该部分菌种丰度增加的物质能够用于治疗强直性脊柱炎或者益于强直性脊柱炎患者服用，能够使其丰度增加的物质不限于治疗强直性脊柱炎的药物和有益肠道菌群平衡的功能性食品。所称组合物中还可能包括另一部分菌种：该部分菌种相较于在健康群体组、在强直性脊柱炎患者群体中显著富集，所称显著富集是指与在健康组中的丰度相比，上述菌种在疾病组中的丰度均具有统计意义地高于或者明显地、实质性地高于在健康组中的丰度；能够使该部分菌种丰度减小的物质能够用于治疗强直性脊柱炎或者益于强直性脊柱炎患者服用，能够使其丰度减小的物质不限于治疗强直性脊柱炎的药物和有益肠道菌群平衡的功能性食品。上述本发明一方面提供的组合物能够用于制备治疗强直性脊柱炎的药物和/或用于制备益于平衡肠道菌群的功能性食品、保健药等。
依据本发明的第三方面，提供一种利用上述任一组合物确定个体的状态的方法，该方法包括：(1)确定所述个体的粪便样本中的所述组合物中的各种菌种的丰度；(2)分别将(1)中确定的组合物中的各种菌种的丰度与其在对照组中的丰度进行比较，依据获得的比较结果确定所述个体的状态，所述对照组由一组或多组相同状态的个体的粪便样本组成，所述状态包括患有强直性脊柱炎和不患有强直性脊柱炎。所称组合物是发明人通过差异比较分析各种肠道微生物在强直性脊柱炎疾病组和健康组的粪便样本中的丰度，以及经过大量已知状态的粪便样本的验证，而确定下来的。
本发明这一方面的方法基于检测个体的粪便样本中的组合物中的各种微生物的丰度，分别将检测确定的各种微生物的丰度与其在对照组中丰度进行比较，依据获得的比较结果能够确定个体为强直性脊柱炎个体或者为健康个体的概率。为早期发现强直性脊柱炎提供一种非侵入性的辅助检测方法。
上述本发明的这一方面的利用组合物确定个体的状态的方法的全部或部分步骤，可以利用包含可拆分的相应单元功能模块的装置/系统来施行，或者将方法程序化、存储于机器可读介质，利用机器运行该可读介质来实现。
依据本发明的第四方面，提供一种利用上述任一实施例中的组合物确定个体的状态的装置，该装置用以实施上述本发明一方面的利用组合物确定个体的状态的方法的全部或部分步骤，该装置包括：数据输入单元，用于输入数据；数据输出单元，用于输出数据；处理器，用于执行可执行程序，执行所述可执行程序包括上述本发明一方面的确定个体状态的方法；存储单元，与所述数据输入单元、所述数据输出单元和所述处理器相连，用于存储数据，其中包括所述可执行程序。上述对本发明一方面的利用组合物确定个体的状态的方法的技术特征和优点的描述，同样适用本发明这一方面的装置，在此不再赘述。
依据本发明的第五方面，提供一种利用上述任一组合物确定个体状态的系统，该系统用以实施本发明一方面的利用组合物确定个体的状态的方法的全部或部分步骤，该系统包括：微生物丰度确定单元，用于确定所述个体的粪便样本中的所述组合物中的各种菌种的丰度；个体状态确定单元，用于分别将来自所述微生物丰度确定单元的所述组合物中的各种菌种的丰度与其在对照组中的丰度进行比较，依据获得的比较结果确定所述个体的状态，所述对照组由一组或多组相同状态的个体的粪便样本组成，所述状态包括患有强直性脊柱炎和不患有强直性脊柱炎。上述对本发明一方面的利用组合物确定个体的状态的方法的技术特征和优点的描述，同样适用本发明这一方面的系统，在此不再赘述。
依据本发明的第六方面，提供一种利用上述任一组合物对多个个体进行分类的方法，该方法包括：确定各个个体的状态，其中包括，(1)确定所述各个个体的粪便样本中的所述组合物中的各种菌种的丰度，(2)分别将(1)中确定的组合物中的各种菌种的丰度与其在对照组中的丰度进行比较，依据获得的比较结果确定所述个体的状态，所述对照组由一组或多组相同状态的个体的粪便样本组成，所述状态包括患有强直性脊柱炎和不患有强直性脊柱炎；依据获得的各个个体的状态对所述各个个体进行分类。该方法能够依据个体的状态的不同区分开多个个体或者区分开多个未知的粪便样本，便于归类、标记管理。另外，上述对本发明一方面的利用组合物确定个体的状态的方法的技术特征和优点的描述，同样适用本发明这一方面的方法，在此不再赘述。
依据本发明的第七方面，提供一种治疗强直性脊柱炎的药物，所述药物促使上述任一组合物中的第一微生物集中的各种菌种的丰度增加，和/或促使所述组合物中的第二微生物集中的各种菌种的丰度减少。所称组合物能够作为强直性脊柱炎微生物标志物，是发明人通过差异分析各种肠道微生物在强直性脊柱炎疾病组和健康组的粪便样本中的丰度，以及经过大量已知状态的粪便样本的验证，而确定下来的。
所称组合物中至少包括一部分菌种，即来自第一微生物集中的菌种，这部分菌种相较于在强直性脊柱炎患者组、在健康群体中显著富集，所称显著富集是指与在疾病组中的丰度相比，上述菌种在健康组中的丰度均具有统计意义地高于或者明显地、实质性地高于在疾病组中的丰度；能够使该部分菌种丰度增加的物质能够用于治疗强直性脊柱炎或者益于强直性脊柱炎患者服用，能够使其丰度增加的物质不限于治疗强直性脊柱炎的药物和有益肠道菌群平衡的功能性食品。所称组合物中还可能包括另一部分菌种，即来自第二微生物集中菌种：该部分菌种相较于在健康群体组、在强直性脊柱炎患者群体中显著富集，所称显著富集是指与在健康组中的丰度相比，上述菌种在疾病组中的丰度均具有统计意义地高于或者明显地、实质性地高于在健康组中的丰度；能够使该部分菌种丰度减小的物质能够用于治疗强直性脊柱炎或者益于强直性脊柱炎患者服用，能够使其丰度减小的物质不限于治疗强直性脊柱炎的药物和有益肠道菌群平衡的功能性食品。
利用本发明这一方面的药物或者功能性食品，合理有效地应用确定的强直性脊柱炎微生物标志物，扶持肠道有益菌的生长和/或抑制肠道潜在致病菌，对于辅助降低血内毒素水平和/或减轻强直性脊柱炎的临床症状具有重要意义。
依据本发明的第八方面，本发明提供一种制备或筛选上述本发明一方面的治疗强直性脊柱炎的药物的方法，该方法包括制备或者筛选促使上述任一组合物中的第一微生物集中的各种菌种的丰度增加、和/或促使所述组合物中的第二微生物集中的各种菌种的丰度减少的物质作为所述药物的步骤。
利用本发明这一方面的生产或筛选治疗强直性脊柱炎的药物的方法，通过合理有效地应用确定的强直性脊柱炎生物标志物进行筛选，能够获得能扶持肠道有益菌的生长和/或抑制肠道潜在致病菌的药物，对于辅助降低血内毒素水平和/或减轻强直性脊柱炎的临床症状具有重要意义。
附图说明
本发明的上述和/或附加的方面和优点从结合下面附图对实施方式的描述中将变得明显和容易理解，其中：
图1是本发明的实施例中的筛选鉴定强直性脊柱炎微生物标志物的试验分析流程示意图。
图2是本发明的实施例中的微生物物种标志物作为诊断指标的AUC评价结果示意图。
具体实施方式
下面详细描述本发明的实施例，所述实施例的示例在附图中示出，其中，自始至终相同或类似的标号表示相同或类似的元件或具有相同或类似功能的元件。下面通过参考附图描述的实施例是示例性的，仅用于解释本发明，而不能理解为对本发明的限制。需要说明的，本文中所使用的术语“第一”或者“第二”等仅为方便描述，不能理解为指示或暗示相对重要性，也不能理解为之间有先后顺序关系。
在本发明的描述中，除非另有说明，“多个”的含义是两个或两个以上。在本文中，除非另有明确的规定和限定，术语“相连”、“连接”等术语应做广义理解，例如，可以是固定连接，也可以是可拆卸连接，或一体地连接；可以是机械连接，也可以是电连接；可以是直接相连，也可以通过中间媒介间接相连，可以是两个元件内部的连通。
生物学标志物是从生物学介质中可以检测到的细胞、生物化学或分子改变。生物学介质包括各种体液、组织、细胞、粪便、头发、呼气等。
所称的某种微生物的丰度指在某一微生物群体中该种微生物的丰富程度，例如在肠道微生物群体中该种微生物的丰富程度，可表示为该种微生物在该群体中的含量。
所称序列的一致性(identity)和序列的相似度(similarity)分别指序列间相同或相似的程度。
根据本发明的一个实施方式提供的一种组合物，其包括第一微生物集中的至少两种菌种，所述第一微生物集由以下菌种组成：毛螺科菌Lachnospiraceae_bacterium_6_1_63FAA、死亡梭杆菌Fusobacterium_mortiferum、拟杆菌Bacteroides_pectinophilus、阴沟肠杆菌Enterobacter_cloacae、挑剔真杆菌Eubacterium_eligens、毛螺科菌Lachnospiraceae_bacterium_8_1_57FAA和梭菌Clostridium_sp._L2-50。
根据本发明的实施例，所述组合物包括所述第一微生物集中的毛螺科菌Lachnospiraceae_bacterium_6_1_63FAA和死亡梭杆菌Fusobacterium_mortiferum；任选的，所述组合物还包括所述第一微生物集中的其它菌种中的任意一种、两种、三种、四种或者其它全部五种。
根据本发明的实施例，所述组合物包括所述第一微生物集中的毛螺科菌Lachnospiraceae_bacterium_6_1_63FAA和阴沟肠杆菌Enterobacter_cloacae；任选的，所述组合物还包括所述第一微生物集中的其它菌种中的任意一种、两种、三种、四种或者其它全部五种。
根据本发明的实施例，所述组合物包括所述第一微生物集中的死亡梭杆菌Fusobacterium_mortiferum和阴沟肠杆菌Enterobacter_cloacae；任选的，所述组合物还包括所述第一微生物集中的其它菌种中的任意一种、两种、三种、四种或者其它全部五种。
根据本发明的实施例，所述组合物包括所述第一微生物集中的阴沟肠杆菌Enterobacter_cloacae和毛螺科菌Lachnospiraceae_bacterium_8_1_57FAA；任选的，所述组合物还包括所述第一微生物集中的其它菌种中的任意一种、两种、三种、四种或者其它全部五种。
根据本发明的实施例，所述组合物还包括第二微生物集中的任意一种、两种、三种、四种、五种、六种、七种、八种、九种或者全部十种菌种，所述第二微生物集由以下菌种组成：放线菌Actinomyces sp.ICM47、短双歧杆菌Bifidobacterium breve、双歧杆菌Bifidobacterium sp.12_1_47BFAA、长双歧杆菌Bifidobacterium longum、产气柯林斯菌Collinsella aerofaciens、两岐双岐杆菌Bifidobacterium bifidum、假小链双歧杆菌Bifidobacterium pseudocatenulatum、迟缓埃格特菌Eggerthella lenta、链状双歧杆菌Bifidobacterium catenulatum和青春双歧杆菌Bifidobacterium adolescentis。
上述本发明一方面或者任一实施例提供的组合物是发明人通过对大量患强直性脊柱炎个体和大量健康对照个体的粪便样本中的微生物的丰度的差异比较分析、以及验证，而确定下来的，明确了肠道微生物中强直性脊柱炎相关的微生物标志物。利用所称的组合物作为强直性脊柱微生物炎标志物，能够确定个体处于患有强直性脊柱炎状态的概率高低或者处于健康状态的概率高低，能够用于非侵入性的早期发现或辅助检测强直性脊柱炎。
根据本发明的另一个实施方式提供的上述任一实施例中的组合物在检测强直性脊柱炎、监测治疗强直性脊柱炎、制备治疗强直性脊柱炎的药物和/或制备功能性食品中的用途。
上述组合物是发明人通过差异比较分析各种肠道微生物在强直性脊柱炎疾病组和健康组的粪便样本中的丰度，以及经过大量已知状态的粪便样本的验证，而确定下来的。组合物中至少包括一部分菌种，即属于第一微生物集的菌种，这部分菌种相较于在强直性脊柱炎患者组、在健康群体中显著富集，所称显著富集是指与在疾病组中的丰度相比，上述菌种在健康组中的丰度均具有统计意义地高于或者明显地、实质性地高于在疾病组中的丰度；能够使该部分菌种丰度增加的物质能够用于治疗强直性脊柱炎或者益于强直性脊柱炎患者服用，能够使其丰度增加的物质不限于治疗强直性脊柱炎的药物和有益肠道菌群平衡的功能性食品。所称组合物中还可能包括另一部分菌种，即属于第二微生物集的菌种，该部分菌种相较于在健康群体组、在强直性脊柱炎患者群体中显著富集，所称显著富集是指与在健康组中的丰度相比，上述菌种在疾病组中的丰度均具有统计意义地高于或者明显地、实质性地高于在健康组中的丰度；能够使该部分菌种丰度减小的物质能够用于治疗强直性脊柱炎或者益于强直性脊柱炎患者服用，能够使其丰度减小的物质不限于治疗强直性脊柱炎的药物和有益肠道菌群平衡的功能性食品。该实施方式提供的组合物能够用于制备治疗强直性脊柱炎的药物和/或用于制备益于平衡肠道菌群的功能性食品、保健药等。
根据本发明的另一个实施方式提供的一种利用上述任一实施例中的组合物确定个体的状态的方法，该方法包括以下步骤(1)和(2)：
(1)确定待测个体中的组合物中的各种菌种的丰度。
确定待测个体的粪便样本中的组合物中的各种菌种的丰度。
根据本发明的实施例，进行以下以完成该步骤：获得所述个体的粪便样本中的核酸序列的测序数据，所述测序数据包括多个读段；将所述读段比对至所述组合物中的各种菌种的基因组，获得比对结果；依据所述比对结果，确定所述组合物中的各种菌种的丰度。
所称的测序数据通过对样本中的核酸序列进行测序得来，测序依据所选的测序平台的不同，可选择但不限于半导体测序技术平台比如PGM、Ion Proton、BGISEQ-100平台，合成边测序的技术平台比如Illumina公司的Hiseq、Miseq序列平台以及单分子实时测序平台比如PacBio序列平台。测序方式可以选择单端测序，也可以选择双末端测序，获得的下机数据是测读出来的片段，称为读段(reads)。
比对可以利用已知比对软件进行，例如SOAP、BWA和TeraMap等，在比对过程中，一般对比对参数进行设置，设置一个或者一对读段(reads)最多允许有s个碱基错配(mismatch)，例如设置s≤2，若reads中有超过s个碱基发生错配，则视为该reads无法比对到(比对上)该组装片段上。所称的获得的比对结果包含各条读段与各物种的基因组的比对情况，包括读段是否能够比对上某个或某些物种的基因组、只唯一比对到一种物种还是比对到多种物种的基因组、比对到物种的基因组的位置、比对到物种基因组的唯一位置还是多个位置等信息。根据本发明的一个实施例，利用SOAPalign 2.21进行比对，设置参数为–r2–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。
根据本发明的一个实施例，所称依据测序数据，所述依据比对结果，确定所述组合物中的各种菌种的丰度，包括利用以下公式确定菌种S的丰度：菌种S的丰度Ab(S)＝Ab(US)+Ab(MS)，其中，Ab(US)＝US/lS，US为唯一比对上该种菌S的读段数目，lS为比对上的该种菌S的基因组的长度，
Figure BDA0000947377640000091
MS为非唯一比对上该种菌S的读段的数目，i表示非唯一比对上该种菌S的读段的编号，Coi为非唯一比对上该种菌S的读段i对应的丰度系数，
Figure BDA0000947377640000092
N为非唯一比对上该种菌S的读段比对上的菌种种类数目，j为非唯一比对上该种菌S的读段比对上的菌种的编号，Uj为唯一比对上菌种j的读段数目。
对于给定的菌种S，其丰度为Ab(S)，与Unique reads和Multiple reads相关，上述公式中的Ab(U)和Ab(M)分别为该物种S的Unique reads和Multiple reads贡献的丰度。每个multiple reads，有特异的物种丰度系数Co，假设一条multiple reads比对上N个物种，可以利用下列公式计算该条multiple reads的Co：
Figure BDA0000947377640000093
即对于这类multiplereads，把其所比对上的N个物种的unique reads的丰度之和作为分母。
(2)丰度比较，以确定个体状态。
分别将(1)中确定的组合物中的各种菌种的丰度与其在对照组中的丰度进行比较，依据获得的比较结果确定所述个体的状态，所述对照组由一组或多组相同状态的个体的粪便样本组成，所述状态包括患有强直性脊柱炎和不患有强直性脊柱炎。所称组合物是发明人通过差异比较分析各种肠道微生物在强直性脊柱炎疾病组和健康组的粪便样本中的丰度，以及经过大量已知状态的粪便样本的验证，而确定下来的。
所称组合物中的各种微生物在对照组中的丰度可以与待测个体的同时测定或者先后测定。根据本发明的一个实施例，将组合物中的各种微生物在对照组中的丰度预先测定保存，用以作为设定阈值的依据。将组合物中的微生物在待测个体样本中的丰度与所称的阈值比较，基于二者的差异程度，进行待测个体的状态的确定。所称阈值可以为一数值或者数值范围，例如，基于某种微生物在对照组中的丰度均值，该微生物对应的阈值可以设为该丰度均值的95％的置信区间。
根据本发明的实施例，对照组由多个患强直性脊柱炎的个体的粪便样本组成，该步骤(2)包括：当(1)中确定的组合物中的第一微生物集中的各种菌种的丰度与其在对照组中的丰度均无显著差异，和/或当(1)中确定的组合物中的第二微生物集中的各种菌种的丰度与其在对照组中的丰度均无显著差异，确定所述个体的状态为患有强直性脊柱炎。根据本发明的实施例，所述显著高于指组合物中的第一微生物集中的各种菌种的丰度分别具有统计意义地高于其在对照组中的丰度，显著性水平α可以取常用的数值，例如设置为0.05。所称无显著差异为步骤(1)中确定的组合物中的第二微生物集中的各种微生物的丰度都落入其在对照组中的丰度的第一预定置信区间。所称的置信区间是指由样本统计量所构造的总体参数的估计区间。在统计学中，一个概率样本的置信区间(Confidence interval)是对这个样本的某个总体参数的区间估计。置信区间展现的是这个参数的真实值有一定概率落在测量结果的周围的程度。置信区间给出的是被测量参数的测量值的可信程度，即前面所要求的“一定概率”。这个概率被称为置信水平。
根据本发明的实施例，对照组由多个健康个体的粪便样本组成，该步骤(2)包括：当(1)中确定的组合物中的第一微生物集中的各种菌种的丰度均显著低于其在对照组中的丰度，和/或当(1)中确定的组合物中的第二微生物集中的各种菌种的丰度都显著高于其在对照组中的丰度，确定所述个体的状态为患有强直性脊柱炎。所称无显著差异为(1)中确定的组合物中的的第一微生物集中的各种菌种的丰度都落入其在对照组中的丰度的第二预定置信区间，任选的所述第二预定置信区间为95％置信区间。所述显著高于或低于指组合物中的第二或第一微生物集中的各种菌种的丰度分别具有统计意义地高于或低于其在对照组中的丰度，显著性水平α可以取常用的数值，例如设置为0.05。
需要说明的是，根据目的或要求不同，可能对确定个体状态结果的可信程度有不同的要求，本领域技术人员可以选择不同的显著性水平(α)，即选择不同的可能犯错误的概率，如此，确定的个体的状态的可信程度为1-α。例如，利用该实施例确定个体处于所确定出的状态95％是可靠的。
该方法基于检测个体的粪便样本中的组合物中的各种菌种的丰度，分别将检测确定的各种菌种的丰度与其在对照组中丰度进行比较，依据获得的比较结果能够确定个体为强直性脊柱炎个体或者为健康个体的概率。为早期发现强直性脊柱炎提供一种非侵入性的辅助检测或者辅助干预治疗的方法。
以上任一实施例中的利用组合物确定个体的状态的方法的全部或部分步骤，可以利用包含可拆分的相应单元功能模块的装置/系统来施行，或者将方法程序化、存储于机器可读介质，利用机器运行该可读介质来实现。
根据本发明的一个实施方式提供的一种利用上述本发明任一实施例中的的组合物确定个体的状态的装置，该装置用以实施上述本发明任一实施例中的利用组合物确定个体的状态的方法的全部或部分步骤，该装置包括：数据输入单元，用于输入数据；数据输出单元，用于输出数据；处理器，用于执行可执行程序，执行所述可执行程序包括完成上述本发明任一实施例的确定个体状态的方法；存储单元，与所述数据输入单元、所述数据输出单元和所述处理器相连，用于存储数据，其中包括所述可执行程序。上述对本发明任一实施例的利用组合物确定个体的状态的方法的技术特征和优点的描述，同样适用本发明这一方面的装置，在此不再赘述。
根据本发明的一个实施方式提供的一种利用上述本发明任一实施例中的组合物确定个体状态的系统，该系统用以实施上述本发明任一实施例中的利用组合物确定个体的状态的方法的全部或部分步骤，该系统包括：微生物丰度确定单元，用于确定所述个体的粪便样本中的所述组合物中的各种菌种的丰度；个体状态确定单元，用于分别将来自所述微生物丰度确定单元的所述组合物中的各种菌种的丰度与其在对照组中的丰度进行比较，依据获得的比较结果确定所述个体的状态，所述对照组由一组或多组相同状态的个体的粪便样本组成，所述状态包括患有强直性脊柱炎和不患有强直性脊柱炎。上述对本发明任一实施例的利用组合物确定个体的状态的方法的技术特征和优点的描述，同样适用本发明这一方面的系统，在此不再赘述。
根据本发明的实施例，所述微生物丰度确定单元用于进行以下：获得所述个体的粪便样本中的核酸序列的测序数据，所述测序数据包括多个读段；将所述读段比对至所述组合物中的各种菌种的基因组，获得比对结果；依据所述比对结果，确定所述组合物中的各种菌种的丰度。
获得所述个体的粪便样本中的核酸序列的测序数据，所述测序数据包括多个读段；将所述读段比对至所述组合物中的各种微生物的基因组，获得比对结果；依据所述比对结果，确定所述组合物中的各种微生物的丰度。
所称的测序数据通过对样本中的核酸序列进行测序得来，测序依据所选的测序平台的不同，可选择但不限于半导体测序技术平台比如PGM、Ion Proton、BGISEQ-100平台，合成边测序的技术平台比如Illumina公司的Hiseq、Miseq序列平台以及单分子实时测序平台比如PacBio序列平台。测序方式可以选择单端测序，也可以选择双末端测序，获得的下机数据是测读出来的片段，称为读段(reads)。
比对可以利用已知比对软件进行，例如SOAP、BWA和TeraMap等，在比对过程中，一般对比对参数进行设置，设置一个或者一对读段(reads)最多允许有s个碱基错配(mismatch)，例如设置s≤2，若reads中有超过s个碱基发生错配，则视为该reads无法比对到(比对上)该组装片段上。所称的获得的比对结果包含各条读段与各物种的基因组的比对情况，包括读段是否能够比对上某个或某些物种的基因组、只唯一比对到一种物种还是比对到多种物种的基因组、比对到物种的基因组的位置、比对到物种基因组的唯一位置还是多个位置等信息。根据本发明的一个实施例，利用SOAPalign 2.21进行比对，设置参数为–r2–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。
根据本发明的一个实施例，所称依据测序数据，所述依据比对结果，确定所述组合物中的各种菌种的丰度，包括利用以下公式确定菌种S的丰度：菌种S的丰度Ab(S)＝Ab(US)+Ab(MS)，其中，Ab(US)＝US/lS，US为唯一比对上该种菌S的读段数目，lS为比对上的该种菌S的基因组的长度，
Figure BDA0000947377640000121
MS为非唯一比对上该种菌S的读段的数目，i表示非唯一比对上该种菌S的读段的编号，Coi为非唯一比对上该种菌S的读段i对应的丰度系数，
Figure BDA0000947377640000122
N为非唯一比对上该种菌S的读段比对上的菌种种类数目，j为非唯一比对上该种菌S的读段比对上的菌种的编号，Uj为唯一比对上菌种j的读段数目。
对于给定的菌种S，其丰度为Ab(S)，与Unique reads和Multiple reads相关，上述公式中的Ab(U)和Ab(M)分别为该物种S的Unique reads和Multiple reads贡献的丰度。每个multiple reads，有特异的物种丰度系数Co，假设一条multiple reads比对上N个物种，可以利用下列公式计算该条multiple reads的Co：
Figure BDA0000947377640000131
即对于这类multiplereads，把其所比对上的N个物种的unique reads的丰度之和作为分母。
根据本发明的实施例，所述个体状态确定单元中的对照组包括多个患强直性脊柱炎的个体的粪便样本，所述个体状态确定单元用于进行以下：当来自微生物丰度确定单元的组合物中的第一微生物集中的各种菌种的丰度与其在对照组中的丰度均无显著差异，和/或当来自微生物丰度确定单元的组合物中的第二微生物集中的各种菌种的丰度都与其在对照组中的丰度无显著差异，确定所述个体的状态为患有强直性脊柱炎。所述无显著差异为来自所述微生物丰度确定单元的组合中的第一或第二微生物集中的各种微生物的丰度落入其在对照组中的丰度的第一预定置信区间，任选的所述第一预定置信区间为95％置信区间。所称的置信区间是指由样本统计量所构造的总体参数的估计区间。在统计学中，一个概率样本的置信区间(Confidence interval)是对这个样本的某个总体参数的区间估计。置信区间展现的是这个参数的真实值有一定概率落在测量结果的周围的程度。置信区间给出的是被测量参数的测量值的可信程度，即前面所要求的“一定概率”。这个概率被称为置信水平。
根据本发明的实施例，所述个体状态确定单元中的对照组包括多个健康个体的粪便样本，所述个体状态确定单元用于进行以下：当来自微生物丰度确定单元的组合物中的第一微生物集中的各种菌种的丰度都显著低于其在对照组中的丰度，和/或当来自微生物丰度确定单元的组合物中的第二微生物集中的各种菌种的丰度都显著高于其在对照组中的丰度，确定所述个体的状态为患有强直性脊柱炎。所称无显著差异为来自微生物丰度确定单元中的组合物中的的第一微生物集中的各种菌种的丰度都落入其在对照组中的丰度的第二预定置信区间，任选的所述第二预定置信区间为95％置信区间。所述显著低于/高于指组合物中的第一/第二微生物集中的各种菌种的丰度分别具有统计意义地低于/高于其在对照组中的丰度，显著性水平α可以取常用的数值，例如设置为0.05。
需要说明的是，根据目的或要求不同，可能对确定个体状态结果的可信程度有不同的要求，本领域技术人员可以选择不同的显著性水平(α)，即选择不同的可能犯错误的概率，如此，确定的个体的状态的可信程度为1-α。例如，利用该实施例确定个体处于所确定出的状态95％是可靠的。
根据本发明的另一个实施方式提供的一种利用上述本发明任一实施例的标志物对多个个体进行分类的方法，该方法包括：分别利用上述本发明任一实施例中的确定个体的状态的方法确定各个个体的状态；依据获得的各个个体的状态对各个个体进行分类。该方法能够依据个体的状态的不同区分开多个个体或者区分开多个未知的粪便样本，便于归类、标记管理。另外，上述对本发明实施例中的利用组合物确定个体的状态的方法的技术特征和优点的描述，同样适用本发明这一方面的方法，在此不再赘述。
根据本发明的又一个实施方式提供的一种治疗强直性脊柱炎的药物，所述药物能够促使上述本发明任一实施例中的组合物中的第一微生物集中的各种菌种的丰度增加，和/或促使所述组合物中的第二微生物集中的各种菌种的丰度减少。所称组合物能够作为强直性脊柱炎微生物标志物，是发明人通过差异分析各种肠道微生物在强直性脊柱炎疾病组和健康组的粪便样本中的丰度，以及经过大量已知状态的粪便样本的验证，而确定下来的。所称组合物中至少包括一部分菌种，即来自第一微生物集中的菌种，这部分菌种相较于在强直性脊柱炎患者组、在健康群体中显著富集，所称显著富集是指与在疾病组中的丰度相比，上述菌种在健康组中的丰度均具有统计意义地高于或者明显地、实质性地高于在疾病组中的丰度；能够使该部分菌种丰度增加的物质能够用于治疗强直性脊柱炎或者益于强直性脊柱炎患者服用，能够使其丰度增加的物质不限于治疗强直性脊柱炎的药物和有益肠道菌群平衡的功能性食品。所称组合物中还可能包括另一部分菌种，即来自第二微生物集中菌种：该部分菌种相较于在健康群体组、在强直性脊柱炎患者群体中显著富集，所称显著富集是指与在健康组中的丰度相比，上述菌种在疾病组中的丰度均具有统计意义地高于或者明显地、实质性地高于在健康组中的丰度；能够使该部分菌种丰度减小的物质能够用于治疗强直性脊柱炎或者益于强直性脊柱炎患者服用，能够使其丰度减小的物质不限于治疗强直性脊柱炎的药物和有益肠道菌群平衡的功能性食品。
利用这一实施方式的药物或者功能性食品，合理有效地应用确定的强直性脊柱炎微生物标志物，扶持肠道有益菌的生长和/或抑制肠道潜在致病菌，可以阻止肠道屏障的缺损，改善并恢复肠道微生态结构，对于辅助降低血内毒素水平和/或减轻强直性脊柱炎的临床症状具有重要意义。
根据本发明的又一个实施方式，提供一种制备或筛选上述任一实施例中的治疗强直性脊柱炎的药物的方法，该方法包括制备或者筛选能够促使上述本发明任一实施例中的组合物中的第一微生物集中的各种菌种的丰度增加、和/或促使所述组合物中的第二微生物集中的各种菌种的丰度减少的物质作为所述药物的步骤。
利用本发明这一方面的生产或筛选治疗强直性脊柱炎的药物的方法，通过合理有效地应用确定的强直性脊柱炎生物标志物进行筛选，能够获得能扶持肠道有益菌的生长和/或抑制肠道潜在致病菌的药物，对于辅助降低血内毒素水平和/或减轻强直性脊柱炎的临床症状具有重要意义。
以下结合具体实施例对本发明的方法和/或装置进行详细的描述。除另有交待，以下实施例中涉及的未特别交待的试剂、序列(接头、标签和引物)、软件及仪器，都是常规市售产品或者开源的，例如购买Illumina的转录组文库构建试剂盒。
以下实施例包括第一阶段和第二阶段，即对应发现阶段和验证阶段。发现阶段包括：基于分析比较73个强直性脊柱炎患者及83个健康对照组的肠道微生物成分及功能改变，以确定物种标志物；验证阶段包括：利用24个AS患者及31个健康对照组验证第一阶段结果的准确性。
实施例1
该示例中，发明人从73个强直性脊柱炎患者、83个健康对照的粪便样品开展整个肠道菌群微生物的关联分析研究描述粪便微生物群落及功能成分特征。总的来说，发明人下载约428.09Gb高质量的测序数据(LC健康人)以及实验测序得到的294Gb高质量测序数据构建了强直性脊柱炎参照基因集，并且和下载的LC基因集、IGC基因集构建一个更加完整的基因集。宏基因组分析显示，17个微生物物种与强直性脊柱炎疾病密切相关，其中7个菌在健康人的肠道微生物中富集，10个菌在强直性脊柱炎病人的肠道微生物中富集。
1、测序数据的获取
1.1样本收集和DNA提取
强直性脊柱炎患者来自杭州浙江中医药大学第一附属医院，实验共采集了73个中国强直性脊柱炎患者的粪便样品，其中每个个体的新鲜粪便样品分成200mg/份，共5份，立即-80℃冰箱冷冻保存。
73个中国强直性脊柱炎患者的粪便样品中提取总DNA。苯酚三氯甲烷处理提取DNA方法提取DNA。
1.2构建文库及测序，以及参照数据下载
DNA建库按仪器制造商(Illumina)的操作指南进行。对文库进行PE100bp测序。Illumina HiSeq2000(Illumina,San Diego,CA)平台对73个样品的文库进行测序。每个样本平均产生4.03Gb(sd.±0.64Gb)高质量测序结果，总计294Gb测序数据量。
从EBI下载83个健康中国人的测序数据，accession号：ERP005860。
参照图1的实验流程，鉴定强直性脊柱炎的相关生物标志物，其中省略的步骤或者细节为本领域技术人员所熟知，几个重要步骤介绍如下面步骤所述。
2、生物标志物的鉴定
2.1测序数据的基本处理
1)测序数据经过质控：获得第一阶段的73个AS样品的测序数据以后，对其进行过滤，质控按以下标准进行：a)去除低质量碱基(Q20)大于50％的reads；b)移除大于5个N碱基的reads；c)移除尾部低质量(Q20)和N碱基。丢失成对reads的序列被认为是单条reads用于组装。
2)同样处理下载的健康人数据。
3)LC基因集获自
Qin,N.et al.Alterations of the human gut microbiome in livercirrhosis.Nature 513,59-64(2014).，从ftp://climb.genomics.cn/pub/10.5524/100001_101000/100064/1.GeneCatalogs/IGC.fa.gz下载得到IGC基因集。
2.2物种丰度分析
利用SOAPalign 2.21将经过上步骤处理的paired-end clean reads比对(匹配)到冗余基因组，这里，所称的冗余基因组来自各数据库中公开的的细菌的参考基因组。比对参数为–r2–m 100–x 1000。Reads与冗余基因组比对，可能被分为两部分：a)Unique reads(U)：reads只比对上一个物种的基因组；这些reads被定义为unique reads。即，如果这些基因组来自同一物种，reads与一个以上的基因组比对，发明人将这些reads定义为uniquereads。b)Multiple reads(M)：如果reads比对上一个物种以上的基因组，定义为multiplereads。即，如果比对上的基因组来自不同物种，发明人定义这些reads为multiple reads。
对于物种S，其丰度为Ab(S)，与特有的U reads和共享的M reads相关，丰度的计算方式如下：
Ab(S)＝Ab(U)+Ab(M)
Ab(U)＝U/l
Figure BDA0000947377640000161
Ab(U)和Ab(M)分别为来自该物种S的unique reads和multiple reads的丰度，l表示基因组的长度。每个multiple reads，有特异物种系数Co；假设某一multiple reads比对上N个物种，按以下方法计算该条multiple reads的Co：
Figure BDA0000947377640000162
对于这些reads，发明人以加和的N个物种的unique reads的丰度作为标准。即对于multiple reads，发明人把其所比对上的N个物种的unique reads丰度之和作为分母。
2.3关联分析/筛选物种标记物
为了获得与强直性脊柱炎疾病密切相关的肠道微生物物种标记物，发明人利用强直性脊柱炎患者组(73例)与健康人组(83例)两组肠道微生物物种丰度数据，在物种级别做了一个与疾病相关性的研究。
基于上步骤得到的物种丰度表，发明人设置标准如下：(1)强直性脊柱炎患者组或健康人组物种丰度的中位数必须大于0.00001；(2)通过结合Benjamini Hochberg的多重检验的Wilcoxon秩和检验进行检验，得到每个物种和强直性脊柱炎疾病的相关性p值和q值；(3)使用一个相对非常严格的阈值(q_values<5.0E-05)。然后，发明人得到28个与强直性脊柱炎疾病密切相关的肠道微生物物种。其中在强直性脊柱炎(AS)患者中富集的有14个物种，在健康人(HD)中富集的有14个物种。这28个微生物物种标记物如下表1所示。
表1
Figure BDA0000947377640000171
Figure BDA0000947377640000181
实施例2
为了证实实施例1中的分析结果可以用以作为强直性脊柱炎物种标志物，进一步比较验证群体中的31个健康人及24个强直性脊柱炎患者的物种丰度，并对最终微生物物种标记物根据验证情况做出筛选。验证群体的测序数据的获得与处理，参照实施例1进行。
验证结果如下：上述富集在健康人群中的14个物种中，7个在验证集中得到高质量的验证(p_values<0.05)，健康人富集的微生物物种标志物验证的p值和q值情况如下表2所示。微生物物种分别为毛螺科菌(Lachnospiraceae bacterium 6_1_63FAA)、死亡梭杆菌(Fusobacterium mortiferum)、拟杆菌([Bacteroides]pectinophilus)、阴沟肠杆菌(Enterobacter cloacae)、挑剔真杆菌(Eubacterium eligens)、毛螺科菌(Lachnospiraceae bacterium8_1_57FAA)、梭菌(Clostridium sp.L2-50)。拟杆菌存在于人体肠道中，主要是帮助分解食物以提供人体所需的营养和能量。
表2
Figure BDA0000947377640000182
对于上述富集在强直性脊柱炎患者中的14个物种，其中10个在验证集中得到高质量的验证(p_values<0.05)，强直性脊柱炎患者富集的微生物物种标志物验证的p值和q值情况如下表3所示，分别为放线菌(Actinomyces sp.ICM47)、短双歧杆菌(Bifidobacteriumbreve)、双歧杆菌(Bifidobacterium sp.12_1_47BFAA)、长双歧杆菌(Bifidobacteriumlongum)、产气柯林斯菌(Collinsella aerofaciens)、两岐双岐杆菌(Bifidobacteriumbifidum)、假小链双歧杆菌(Bifidobacterium pseudocatenulatum)、迟缓埃格特菌(Eggerthella lenta)、链状双歧杆菌(Bifidobacterium catenulatum)、青春双歧杆菌(Bifidobacterium adolescentis)。
双歧杆菌可以预防腹泻，减少便秘，即双向调节。有研究表明超过70％的强直性脊柱炎患者患有肠道炎症，这些患者中5～10％肠道炎症很严重，可能会发展为临床上的炎症性肠病(IBD)或克罗恩氏病(CD)(Mielants et al.,1985)，研究发现IBD病人肠道中拟杆菌、梭状芽胞杆菌、厚壁菌下降，放线菌、变形菌比例升高。
表3
Figure BDA0000947377640000191
发明人认为，可以将从健康人富集的7个微生物物种标记物，作为强直性脊柱炎疾病患病的反向指标，或作为治疗强直性脊柱炎进行研发的微生物制剂药物菌群成分，或作为检测强直性脊柱炎、监测强直性脊柱炎治疗进程的恢复指标；将强直性脊柱炎患者富集的10个微生物物种标记物，作为强直性脊柱炎疾病患病的正向指标，特别用于强直性脊柱炎疾病非创伤式的检测和诊断。
发明人利用这17个微生物物种标记物，构建一个综合指标，估计ROC(Receiver-operating characteristic)曲线下面积AUC，AUC越大，表示诊断能力越高，评价综合得分对应其对强直性脊柱炎的诊断能力。通过对一期(第一阶段)的156个样品和二期(第二阶段)的55个样品进行评测，具体情况如图2所示，都表现出了很好的诊断能力，在一期得到AUC＝96.62％，如图2左边图所示，置信区间为93.93％-99.31％；二期得到AUC＝89.92％，如图2右边图所示，置信区间为81.28％-98.56％。
实施例3
利用45个粪便样本进行样本来源的个体状态的检测。
参照实施例2的方法确定各粪便样本中的表2的7种菌种中的至少两种的丰度，判断各样本中的这些菌种的丰度是否落入各自在疾病对照组或者健康对照组的丰度的95％的置信区间，判定这7种菌种中的两种、三种、四种、五种、六种或者全部七种的丰度均落入健康组的对应区间的样本所对应的个体的状态为非强直性脊柱炎患者，判定均落入疾病组的对应区间的样本所对应的个体的状态为强直性脊柱炎患者。
进一步地，使组合物包括表3中的至少一种在AS组中富集的菌种。参照实施例2的方法确定各粪便样本中的表3的放线菌(Actinomyces sp.ICM47)、短双歧杆菌(Bifidobacterium breve)、双歧杆菌(Bifidobacterium sp.12_1_47BFAA)、长双歧杆菌(Bifidobacterium longum)、产气柯林斯菌(Collinsella aerofaciens)、两岐双岐杆菌(Bifidobacterium bifidum)、假小链双歧杆菌(Bifidobacterium pseudocatenulatum)、迟缓埃格特菌(Eggerthella lenta)、链状双歧杆菌(Bifidobacterium catenulatum)和青春双歧杆菌(Bifidobacterium adolescentis)的丰度。结合上步的判断结果，进一步根据待测样本中的来自的表3的菌株的丰度是否落入各自在疾病对照组或者健康对照组的丰度的95％的置信区间来判断个体的状态，若表3中的菌种的丰度落入疾病组的对应区间，判定样本所对应的个体进一步倾向于患有强直性脊柱炎患者；若表3中的菌种的丰度落入健康组的对应区间，判定样本所对应的个体的状态进一步倾向于为非强直性脊柱炎患者。
采用不同的组合物菌种组成，结果显示：组合物只包括2种表2中的菌种和1种表3中的菌种，基本能够对全部的45个样本进行判断，但判断结果与记录的样本源自的个体的状态一致的程度达约75％；组合物由表2中的全部菌种组成，能够对其中的40个样本进行判断，判断结果与记录的个体状态的一致性达92％；组合物由表2中的全部菌种以及表3中的放线菌和产气柯林斯菌组成，能够对其中的39个样本进行状态判断，判断结果与记录的个体状态的一致性达90％。
另外，发明人发现对表2和表3中的物种联合检测，例如检测表3中的物种标志物被富集，同时表2中物种标记物不被富集，能够更准确的判断发现强直性脊柱炎患者或易感人群。
在利用标志物治疗强直性脊柱炎的方案中，发明人发现使3中的物种标志物生长得到抑制或者清除，同时使表2中物种标记物被富集，治疗效果极佳。
在本说明书的描述中，参考术语“一个实施例”、“一些实施例”、“示例”、“具体示例”、或“一些示例”等的描述意指结合该实施例或示例描述的具体特征、结构、材料或者特点包含于本发明的至少一个实施例或示例中。在本说明书中，对上述术语的示意性表述不一定指的是相同的实施例或示例。而且，描述的具体特征、结构、材料或者特点可以在任何的一个或多个实施例或示例中以合适的方式结合。
尽管已经示出和描述了本发明的实施例，本领域的普通技术人员可以理解：在不脱离本发明的原理和宗旨的情况下可以对这些实施例进行多种变化、修改、替换和变型，本发明的范围由权利要求及其等同物限定。

Claims (1)

1.一种用于判断个体是否患有强直性脊柱炎的组合物，其包括第一微生物集，所述第一微生物集由以下菌种组成：毛螺科菌Lachnospiraceae_bacterium_6_1_63FAA、死亡梭杆菌Fusobacterium_mortiferum、拟杆菌Bacteroides_pectinophilus、阴沟肠杆菌Enterobacter_cloacae、挑剔真杆菌Eubacterium_eligens、毛螺科菌Lachnospiraceae_bacterium_8_1_57FAA和梭菌Clostridium_sp._L2-50，
所述组合物还包括第二微生物集中的至少一种菌种，所述第二微生物集由以下菌种组成：放线菌Actinomyces sp.ICM47、短双歧杆菌Bifidobacterium breve、双歧杆菌Bifidobacterium sp.12_1_47BFAA、长双歧杆菌Bifidobacterium longum、产气柯林斯菌Collinsella aerofaciens、两岐双岐杆菌Bifidobacterium bifidum、假小链双歧杆菌Bifidobacterium pseudocatenulatum、迟缓埃格特菌Eggerthella lenta、链状双歧杆菌Bifidobacterium catenulatum和青春双歧杆菌Bifidobacterium adolescentis。
CN201610165123.XA 2016-03-22 2016-03-22 组合物及应用 Active CN105733988B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610165123.XA CN105733988B (zh) 2016-03-22 2016-03-22 组合物及应用

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610165123.XA CN105733988B (zh) 2016-03-22 2016-03-22 组合物及应用

Publications (2)

Publication Number Publication Date
CN105733988A CN105733988A (zh) 2016-07-06
CN105733988B true CN105733988B (zh) 2020-09-01

Family

ID=56251135

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610165123.XA Active CN105733988B (zh) 2016-03-22 2016-03-22 组合物及应用

Country Status (1)

Country Link
CN (1) CN105733988B (zh)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112094904B (zh) * 2020-09-30 2022-08-16 中国人民解放军军事科学院军事医学研究院 急性电离辐射早期损伤检测的微生物种水平标志物及应用
CN114317784B (zh) * 2021-12-08 2024-02-13 上海锐翌医学检验实验室有限公司 白塞病标志微生物及其应用

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105296590A (zh) * 2015-09-30 2016-02-03 上海锐翌生物科技有限公司 大肠癌标志物及其应用
CN105378739A (zh) * 2012-12-13 2016-03-02 麦太宝根有限公司 鉴定有发生2型糖尿病风险的人

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105378739A (zh) * 2012-12-13 2016-03-02 麦太宝根有限公司 鉴定有发生2型糖尿病风险的人
CN105296590A (zh) * 2015-09-30 2016-02-03 上海锐翌生物科技有限公司 大肠癌标志物及其应用

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
Intestinal Dysbiosis in Ankylosing Spondylitis;Mary-Ellen Costello等;《ARTHRITIS & RHEUMATOLOGY》;20150331;686–691 *

Also Published As

Publication number Publication date
CN105733988A (zh) 2016-07-06

Similar Documents

Publication Publication Date Title
CN107217089B (zh) 确定个体状态的方法及装置
CN107217088B (zh) 强直性脊柱炎微生物标志物
CN107034279A (zh) 结核病微生物标志物在制备诊断结核病的试剂中的应用
CN105132518A (zh) 大肠癌标志物及其应用
CN110838365A (zh) 肠易激综合症相关菌群标志物及其试剂盒
CN113913490B (zh) 非酒精性脂肪肝标志微生物及其应用
CN114182007A (zh) 白塞病标志基因及其应用
WO2017156739A1 (zh) 分离的核酸及应用
CN105733988B (zh) 组合物及应用
CN107217086B (zh) 疾病标志物及应用
CN105671177B (zh) 强直性脊柱炎标志物及应用
WO2017156764A1 (zh) 分离的核酸及应用
CN113999922A (zh) 急性腹泻标志微生物及其应用
CN114836508A (zh) 慢性阻塞性肺病标志微生物及其应用
CN113930479B (zh) 系统性红斑狼疮标志微生物及其应用
CN114317674B (zh) 类风湿性关节炎标志微生物及其应用
CN114317784B (zh) 白塞病标志微生物及其应用
CN113999923B (zh) 小柳-原田综合症标志微生物及其应用
CN114317718A (zh) 类风湿性关节炎标志ko及其应用
CN114262743B (zh) 中风标志微生物及其应用
CN114517235A (zh) 肌痛性脑脊髓炎标志基因及其应用
CN114292932A (zh) 慢性传输便秘标志ko及其应用
CN114317717A (zh) 儿童克罗恩病标志基因及其应用
CN114317716A (zh) 痛风症标志基因及其应用
CN114292933A (zh) 格雷夫斯病标志基因及其应用

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
CB02 Change of applicant information

Address after: Room 202-3 and 302, Building No. 138, Xinjunhuan Road, Minhang District, Shanghai, 20114

Applicant after: SHANGHAI RUIYI BIOTECHNOLOGY CO., LTD.

Address before: Room 119, 1st floor, 3058 Pusan Road, Pudong New Area, Shanghai 200050

Applicant before: SHANGHAI RUIYI BIOTECHNOLOGY CO., LTD.

CB02 Change of applicant information
GR01 Patent grant
GR01 Patent grant