CN112614596B - 一种肠道菌群移植治疗溃疡性结肠炎的供受体配型方法 - Google Patents
一种肠道菌群移植治疗溃疡性结肠炎的供受体配型方法 Download PDFInfo
- Publication number
- CN112614596B CN112614596B CN202011527189.1A CN202011527189A CN112614596B CN 112614596 B CN112614596 B CN 112614596B CN 202011527189 A CN202011527189 A CN 202011527189A CN 112614596 B CN112614596 B CN 112614596B
- Authority
- CN
- China
- Prior art keywords
- donor
- pwy
- ulcerative colitis
- recipient
- matching
- 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
Links
- 206010009900 Colitis ulcerative Diseases 0.000 title claims abstract description 58
- 201000006704 Ulcerative Colitis Diseases 0.000 title claims abstract description 58
- 238000000034 method Methods 0.000 title claims abstract description 49
- 238000002054 transplantation Methods 0.000 title claims abstract description 28
- 230000000968 intestinal effect Effects 0.000 title claims abstract description 25
- 239000000370 acceptor Substances 0.000 claims abstract description 25
- 239000000203 mixture Substances 0.000 claims abstract description 21
- 238000012545 processing Methods 0.000 claims abstract description 4
- 230000037361 pathway Effects 0.000 claims description 26
- 239000011159 matrix material Substances 0.000 claims description 22
- 230000009286 beneficial effect Effects 0.000 claims description 21
- 241000894007 species Species 0.000 claims description 13
- 241000192031 Ruminococcus Species 0.000 claims description 9
- 238000002560 therapeutic procedure Methods 0.000 claims description 9
- 230000002550 fecal effect Effects 0.000 claims description 7
- 238000011160 research Methods 0.000 claims description 7
- 241000186394 Eubacterium Species 0.000 claims description 5
- 238000004364 calculation method Methods 0.000 claims description 5
- 210000001035 gastrointestinal tract Anatomy 0.000 claims description 5
- 230000002939 deleterious effect Effects 0.000 claims description 4
- 241000702460 Akkermansia Species 0.000 claims description 3
- 241000193830 Bacillus <bacterium> Species 0.000 claims description 3
- 241000588914 Enterobacter Species 0.000 claims description 3
- 241000194033 Enterococcus Species 0.000 claims description 3
- 241001531192 Eubacterium ventriosum Species 0.000 claims description 3
- 241001112693 Lachnospiraceae Species 0.000 claims description 3
- 241000122116 Parvimonas Species 0.000 claims description 3
- 241000351207 Peptoniphilus Species 0.000 claims description 3
- 238000012300 Sequence Analysis Methods 0.000 claims description 3
- 241000194017 Streptococcus Species 0.000 claims description 3
- 238000012163 sequencing technique Methods 0.000 claims description 3
- 241000662772 Flavonifractor Species 0.000 claims description 2
- 241000785902 Odoribacter Species 0.000 claims description 2
- 241000093136 Sellimonas Species 0.000 claims description 2
- 241000607768 Shigella Species 0.000 claims description 2
- 210000000936 intestine Anatomy 0.000 claims description 2
- 150000001875 compounds Chemical class 0.000 claims 2
- 241000605314 Ehrlichia Species 0.000 claims 1
- 239000000126 substance Substances 0.000 claims 1
- 230000007613 environmental effect Effects 0.000 description 4
- 241000566145 Otus Species 0.000 description 2
- 241000191992 Peptostreptococcus Species 0.000 description 2
- 241000607142 Salmonella Species 0.000 description 2
- 230000001580 bacterial effect Effects 0.000 description 2
- 230000002950 deficient Effects 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 201000010099 disease Diseases 0.000 description 2
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 description 2
- 238000011156 evaluation Methods 0.000 description 2
- 239000002245 particle Substances 0.000 description 2
- NBGAYCYFNGPNPV-UHFFFAOYSA-N 2-aminooxybenzoic acid Chemical class NOC1=CC=CC=C1C(O)=O NBGAYCYFNGPNPV-UHFFFAOYSA-N 0.000 description 1
- 208000004998 Abdominal Pain Diseases 0.000 description 1
- 206010012735 Diarrhoea Diseases 0.000 description 1
- 201000000297 Erysipelas Diseases 0.000 description 1
- 241001531182 Eubacterium xylanophilum Species 0.000 description 1
- 208000035126 Facies Diseases 0.000 description 1
- 208000012671 Gastrointestinal haemorrhages Diseases 0.000 description 1
- 208000022559 Inflammatory bowel disease Diseases 0.000 description 1
- 241000736262 Microbiota Species 0.000 description 1
- 108060008682 Tumor Necrosis Factor Proteins 0.000 description 1
- 102100040247 Tumor necrosis factor Human genes 0.000 description 1
- 206010054094 Tumour necrosis Diseases 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 230000000259 anti-tumor effect Effects 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 239000003246 corticosteroid Substances 0.000 description 1
- 229960001334 corticosteroids Drugs 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000012217 deletion Methods 0.000 description 1
- 230000037430 deletion Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 210000003608 fece Anatomy 0.000 description 1
- 208000030304 gastrointestinal bleeding Diseases 0.000 description 1
- 230000002068 genetic effect Effects 0.000 description 1
- 244000005709 gut microbiome Species 0.000 description 1
- 230000001900 immune effect Effects 0.000 description 1
- 239000002955 immunomodulating agent Substances 0.000 description 1
- 229940121354 immunomodulator Drugs 0.000 description 1
- 208000027866 inflammatory disease Diseases 0.000 description 1
- 238000002360 preparation method Methods 0.000 description 1
- 108090000623 proteins and genes Proteins 0.000 description 1
- 230000000306 recurrent effect Effects 0.000 description 1
- 230000035945 sensitivity Effects 0.000 description 1
- 238000001356 surgical procedure Methods 0.000 description 1
- 208000024891 symptom Diseases 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 230000004580 weight loss Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16H—HEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
- G16H50/00—ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics
- G16H50/70—ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics for mining of medical data, e.g. analysing previous cases of other patients
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16H—HEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
- G16H50/00—ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics
- G16H50/50—ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics for simulation or modelling of medical disorders
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16H—HEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
- G16H70/00—ICT specially adapted for the handling or processing of medical references
- G16H70/60—ICT specially adapted for the handling or processing of medical references relating to pathologies
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/10—Information and communication technologies [ICT] supporting adaptation to climate change, e.g. for weather forecasting or climate simulation
Landscapes
- Engineering & Computer Science (AREA)
- Medical Informatics (AREA)
- Health & Medical Sciences (AREA)
- Public Health (AREA)
- General Health & Medical Sciences (AREA)
- Data Mining & Analysis (AREA)
- Epidemiology (AREA)
- Primary Health Care (AREA)
- Databases & Information Systems (AREA)
- Pathology (AREA)
- Biomedical Technology (AREA)
- Investigating Or Analysing Biological Materials (AREA)
- Acyclic And Carbocyclic Compounds In Medicinal Compositions (AREA)
Abstract
本发明公开了一种肠道菌群移植治疗溃疡性结肠炎的供受体配型方法，属于临床医疗技术领域，包括如下步骤：(S1)：收集临床溃疡性结肠炎患者受体以及健康人员供体的肠道菌群数据，建立数据库；(S2)：通过分析处理数据库数据，提取供体独有多样性、菌属分类单元、通路分类单元以及用以比较供体与受体之间组成差异的供受体距离数据组成特征，构建溃疡性结肠炎供受体配型模型；(S3)：配型模型通过层次分析法分解特征层次结构，对供体库中供体进行打分排序，选择最优值供体，使患者获得最优供体配型。本发明公开的肠道菌群移植治疗溃疡性结肠炎的供受体配型方法，便于患者匹配合适的菌群供体，以保证菌群移植治疗的有效性。
Description
技术领域
本发明涉及临床医疗技术领域，尤其涉及一种肠道菌群移植治疗溃疡性结肠炎的供受体配型方法。
背景技术
溃疡性结肠炎(UC)是炎症性肠病的一种亚型，是一种缓解性和复发性炎症性疾病，症状包括腹痛、腹泻、粪便急迫性、胃肠出血和体重减轻。UC已成为21世纪的一种全球性疾病，在西方国家的合并发病率为每10万人中6-15人，在新兴工业化国家发病率在加速上升。尽管有大量可用的治疗方法，包括皮质类固醇、抗肿瘤坏死因子(TNF-a)制剂、氨水杨酸盐、免疫调节剂和手术，但许多患者对这些治疗没有反应或在治疗期间出现继发性失效。因此，开发新疗法和研究替代策略是必要的。虽然UC的病因尚不清楚，但它被认为是复杂的遗传、免疫和环境因素。
肠道菌群作为一种重要的环境因子，在UC发育过程中起着重要的作用。通过粪便微生物群移植(FMT)操纵肠道微生物群已经在临床试验中证明了对UC缓解以及其他疾病的有效性。然而，FMT的有效性在不同的研究中存在差异，在这些研究中，特定的供体可能发挥关键作用。例如，使用5名供体对UC进行的FMT随机对照试验显示，78％的患者在接受了来自单一供体的粪便后达到了缓解。此外，接受来自同一供体粪便的FMT批次的患者表现出更高的缓解率(37％)，而接受FMT批次不包括该供体的患者(18％)。因此，如何准确选取溃疡性结肠炎医疗研究的特征参数，从而为溃疡性结肠炎患者匹配合适的菌群供体以保证菌群移植治疗的有效性，是菌群移植医疗技术急需解决的技术难题。
发明内容
为了克服现有技术的缺陷，本发明所要解决的技术问题在于提出一种肠道菌群移植治疗溃疡性结肠炎的供受体配型方法，通过本方法，便于溃疡性结肠炎患者匹配合适的菌群供体，以保证菌群移植治疗的有效性。
为达此目的，本发明采用以下技术方案：
本发明提供的一种肠道菌群移植治疗溃疡性结肠炎的供受体配型方法，包括如下步骤：(S1)：收集临床溃疡性结肠炎患者受体以及健康人员供体的肠道菌群数据，建立数据库；(S2)：通过分析处理数据库数据，提取供体独有多样性、菌属分类单元、通路分类单元以及用以比较供体与受体之间组成差异的供受体距离数据组成特征，构建溃疡性结肠炎供受体配型模型，并对各指标特征进行计算分析；(S3)：配型模型通过层次分析法分解特征层次结构，对供体库中供体进行打分排序，计算全局权重，选择最优值供体，使患者获得最优供体配型。
本发明优选地技术方案在于，在步骤(S2)中，采用求同存异法，根据供体和受体肠道内多样性的种类进行比对，找出供体独有多样性特征；
求同存异法计算供体独有多样性的公式为：
Obsalpha-diver＝ObsDonor-ObsDonor∩patient；
本发明优选地技术方案在于，在步骤(S2)中，在溃疡性结肠炎患者中选择至少一种富集或缺失的菌属分类单元，组成菌属分类特征；
利用增益减害原则计算菌属分类单元，其公式为：
本发明优选地技术方案在于，在步骤(S2)中，在溃疡性结肠炎患者中选择至少一种富集或缺失的通路分类单元，组成通路分类特征；
利用增益减害原则计算通路分类单元，其公式为：
本发明优选地技术方案在于，供体和受体肠道内多样性的种类是基于数据库中供体和受体的粪便样本进行16S序列分析得到的。
本发明优选地技术方案在于，在步骤(S3)中，层次分析法进行供受体配型时包括如下步骤：
(S31)：将溃疡性结肠炎患者供受体配型涉及的因素划分为目标层、准则层以及方案层，建立供受体配型层次结构模型，其中，准则层包括供体独有多样性、有益菌属分类单元、有害菌属分类单元、有益通路分类单元、有害通路分类单元以及距离；
(S32)：构造两两比较判断矩阵，判断矩阵用数值表示下层元素对上层元素的相对重要性；
(S33)：进行层次单排序，计算单一准则下元素的相对权值，通过随机一致性比率CR来检验判断矩阵的一致性，当CR<0.1时，认为判断矩阵具有满意的一致性，当CR≥0.1时，则需重新调整判断矩阵，使其具有一致性；
(S34)：计算各层元素对目标层的总排序权重，检验一致性；
(S35)：根据综合权重得分进行排序，选择权重高的供体作为最优供体。
本发明优选地技术方案在于，供受体距离计算公式为：
其中，Sp，i、SD，i分别代表受体和供体第i个多样性种类的个数。
本发明优选地技术方案在于，肠道菌群移植治疗溃疡性结肠炎的供受体配型方法，还包括如下步骤：(S4)：将根据步骤(S1)、步骤(S2)以及步骤(S3)确定的供受体配型模型运用于临床研究，与临床缓解率相结合评定其准确率。
本发明优选地技术方案在于，溃疡性结肠炎患者中富集的菌属分类单元为Sellimonas、Parvimonas、消化链球菌属、埃希志贺氏杆菌属、Flavonifractor、Erysipelatoclostridium、链球菌属、Peptoniphilus、肠球菌属中的至少一种；溃疡性结肠炎患者中缺失的菌属分类单元为考拉杆菌属、瘤胃菌科_UCG-002、瘤胃球菌属_2、粪球菌属_2、Akkermansia、Odoribacter、[Eubacterium]_ventriosum_group、毛螺菌科_NK4A136_group、瘤胃菌科_UCG-014、瘤胃菌科_UCG-013、un_f_瘤胃菌科、Fusicatenibacter、Ruminiclostridium_6、Parabacteroides、[Eubacterium]_xylanophilum_group、氏菌属、Parasutterella、Christensenellaceae_R-7_group、Alloprevotella、毛螺菌科_UCG-008、[Eubacterium]_hallii_group、Enterorhabdus、Senegalimassilia、un_o_Gastranaerophilales、Tyzzerella_3中的至少一种。
本发明优选地技术方案在于，溃疡性结肠炎患者中富集的通路分类单元为PWY-6470、HEMESYN2-PWY、GALACTARDEG-PWY、METHGLYUT-PWY、PWY-6891、PWY-5910、ORNDEG-PWY、3-HYDPHENACETATE-PWY、KDO-NAGLIPASYN-PWY、NAD-BIOSYNTHESIS-II、PPGPPMET-PWY、RUMP-PWY、PWY-7377中的至少一种；溃疡性结肠炎患者中缺失的通路分类单元为PWY-7374、PWY-6151、PWY-7090、PWY-6572、PWY-7111、COBALSYN-PWY、PWY-6478、PWY-5005、GLUTORN-PWY、PWY4FS-7、PWY-5103、PWY-5101、PANTOSYN-PWY、PWY-6507、ILEUSYN-PWY、PWY-3001、PWY-6385、PWY-5659、PWY-5676的至少一种。
本发明的有益效果为：
本发明提供的肠道菌群移植治疗溃疡性结肠炎的供受体配型方法，结合供体和受体的信息，根据相生相克的原则采用相关方法计算其供体独有的多样性种类、菌属分类特征、通路分类特征、距离，通过层次分析法建立层次结构模型计算相应的权重来确定配型模型，并应用于预先临床试验进行性能比较，评定模型的高灵敏性，从而验证本发明策略方法的可行性。通过上述过程，便于患者匹配合适的菌群供体，以保证菌群移植治疗的有效性。
附图说明
图1是本发明具体实施方式中提供的肠道菌群移植治疗溃疡性结肠炎的供受体配型方法的原理流程示意图。
具体实施方式
下面结合附图并通过具体实施方式来进一步说明本发明的技术方案。本领域技术人员应当理解，下面所具体描述的内容是说明性的而非限制性的，不应以此限制本发明的保护范围。
如图1所示，为了便于患者匹配合适的菌群供体，以保证菌群移植治疗的有效性。进一步地，本实施例中提供的一种肠道菌群移植治疗溃疡性结肠炎的供受体配型方法，包括如下步骤：(S1)：收集临床溃疡性结肠炎患者受体以及健康人员供体的肠道菌群数据，建立数据库；(S2)：通过分析处理数据库数据，提取供体独有多样性、菌属分类单元、通路分类单元以及用以比较供体与受体之间组成差异的供受体距离数据组成特征，构建溃疡性结肠炎供受体配型模型，并对各指标特征进行计算分析；(S3)：配型模型通过层次分析法分解特征层次结构，对供体库中供体进行打分排序，计算全局权重，选择最优值供体，使患者获得最优供体配型。其中，供体独有多样性特征指供体肠道中不同于受体肠道内的独有多样性菌种，菌属分类特征指在溃疡性结肠炎患者中富集或缺失的菌属分类单元，通路分类单元特征指在溃疡性结肠炎患者中富集或缺失的通路分类单元，供受体距离特征指Bray-Curtis距离，主要基于OTUs的计数统计，比较供体群落和受体群落之间的组成差异，DBray-Curtis值越小表示二者组成差异小。通过提取上述多角度多方位的供受体数据特征来建立配型模型，能提供更多的信息，使得配型模型特征多样化，配型匹配度更高，提高菌群移植治疗效果。在层次分析法供体与受体配型中，层次分析法是将为患者匹配合适的菌群供体的决策问题按总目标、各层子目标、评价准则直至具体的供体配型备投方案的顺序分解为不同的层次结构，然后用求解判断矩阵特征向量的办法，求得每一层次的各元素对上一层次某元素的优先权重，最后再加权和的方法递阶归并各供体配型备择方案对为患者匹配合适的菌群供体的总目标的最终权重，此最终权重最大者即为最优方案。通过上述过程，便于患者匹配合适的菌群供体，以保证菌群移植治疗的有效性。
优选地，在步骤(S2)中，采用求同存异法，根据供体和受体肠道内多样性的种类进行比对，找出供体独有多样性特征；
求同存异法计算供体独有多样性的公式为：
Obsalpha-diver＝ObsDonor-ObsDonor∩patient；
式中，Obsalpha-diver代表供体独有的多样性种类个数；ObsDonor代表供体实际的多样性种类个数；ObsDonor∩patient代表供体和受体共有的多样性种类个数；
代表供体独有多样性种类的相对丰度。通过上述公式，可计算环境样本内供体独有的物种多样性大小，并对供体各个独有物种多样性的具体数据进行打分排序。
优选地，在步骤(S2)中，在溃疡性结肠炎患者中选择至少一种富集或缺失的菌属分类单元，组成菌属分类特征；
利用增益减害原则计算菌属分类单元，其公式为：
式中，Ratio.btDonor代表受体缺少的有益菌属分类单元但在供体中出现的有益菌属分类单元个数在受体缺少的有益菌属分类单元总数的比例；Ratio.htDonor代表在供体和受体中均出现的有害菌属分类单元的个数在受体出现的有害菌属分类单元总数的比例；Abu.btDonor和Abu.htDonor分别代表对应的相对丰度。通过上述公式，可计算环境样本内受体在富集或缺失的有益菌属或者有害菌属分类单元的比例构成以及相对丰度，并对不同组成的菌属分类单元的具体数据进行打分排序，以便构建更准确可靠的配型模型。
优选地，在步骤(S2)中，在溃疡性结肠炎患者中选择至少一种富集或缺失的通路分类单元，组成通路分类特征；
利用增益减害原则计算通路分类单元，其公式为：
式中，Ratio.bpDonor代表受体缺少的有益通路分类单元但在供体中出现的有益通路分类单元个数在受体缺少的有益通路分类单元总数的比例；Ratio.hpDonor代表在供体和受体中均出现的有害通路分类单元的个数在受体出现的有害通路分类单元总数的比例；Abu.bpDonor和Abu.hpDonor分别代表对应的相对丰度。通过上述公式，可计算环境样本内受体在富集或缺失的有益通路或者有害通路分类单元的比例构成以及相对丰度，并对不同组成的通路分类单元的具体数据进行打分排序，以便构建更准确可靠的配型模型。
优选地，供体和受体肠道内多样性的种类是基于数据库中供体和受体的粪便样本进行16S序列分析得到的。通过16S序列分析法对肠道菌群进行高通量测序，通过分析测序序列的构成了解特定环境中肠道菌群的构成情况或基因的组成以及功能，核心是研究菌群样本中的物种分类、物种丰度以及系统进化，通过OTU分析，可以知道研究样本中的菌群多样性和不同菌群的丰度。
优选地，在步骤(S3)中，层次分析法进行供受体配型时包括如下步骤：
(S31)：将溃疡性结肠炎患者供受体配型涉及的因素划分为目标层、准则层以及方案层，建立供受体配型层次结构模型，其中，准则层包括供体独有多样性、有益菌属分类单元、有害菌属分类单元、有益通路分类单元、有害通路分类单元以及距离；
(S32)：构造两两比较判断矩阵，判断矩阵用数值表示下层元素对上层元素的相对重要性；
(S33)：进行层次单排序，计算单一准则下元素的相对权值，通过随机一致性比率CR来检验判断矩阵的一致性，当CR<0.1时，认为判断矩阵具有满意的一致性，当CR≥0.1时，则需重新调整判断矩阵，使其具有一致性；
(S34)：计算各层元素对目标层的总排序权重，检验一致性；
(S35)：根据综合权重得分进行排序，选择权重高的供体作为最优供体。
在步骤(S32)中，构造两两比较判断矩阵，所述判断矩阵用数值表示下层元素对上层元素的相对重要性，将比较的结果用数值表示。一般使用数字1～9作为标度(见表1)。
表1：因素比值规则表
例如，对准则层的六个特征因素(供体独有多样性、有益菌属分类单元、有害菌属分类单元、有益通路分类单元、有害通路分类单元以及距离)进行专家打分(见表2)，
表2：六个特征因素的打分表
因此构造的判断矩阵M为：
由判断矩阵M计算出6个特征指标的全局权重分别为0.072，0.433，0.210，0.037，0.210，0.037。
在步骤(S33)中，进行层次单排序，计算单一准则下元素的相对权值，采用随机一致性比率CR来检验判断矩阵的一致性，其计算公式为：
CR＝CI/RI (1)，
公式(1)中RI为随机一致性指标(见表3)，CI为一致性指标，公式(2)中n为本层指标的个数，λmax为判断矩阵的最大特征根；当CR<0.1时，认为判断矩阵具有满意的一致性，反之，则需重新调整判断矩阵，使其具有一致性。
表3：RI数值表
在步骤(S34)中，进行层次综合排序，检验一致性，即计算各层元素对目标层的总排序权重，计算公式为：
公式(3)中ckj表示C层包含的n个元素相对于上一层B层因素的层次单排序，bj为B层单排序权值。
通过上述过程，根据综合权重得分进行排序，选择权重高的作为最优供体，便于患者匹配最合适的菌群供体，以保证菌群移植治疗的有效性。
优选地，供受体距离计算公式为：
其中，Sp，i、SD，i分别代表受体和供体第i个多样性种类的个数。Bray-Curtis距离主要基于OTUs的计数统计，比较供体群落和受体群落之间的组成差异，值越小表示二者组成差异小。通过上述过程，以便根据构成准则层的距离因素的具体数据进行打分排序。
优选地，肠道菌群移植治疗溃疡性结肠炎的供受体配型方法，还包括如下步骤：(S4)：将根据步骤(S1)、步骤(S2)以及步骤(S3)确定的供受体配型模型运用于临床研究，与临床缓解率相结合评定其准确率。通过本方法构建的供受体配型模型的准确率处于50％以上，保证菌群移植治疗的有效性。
优选地，溃疡性结肠炎患者中富集的菌属分类单元为Sellimonas、Parvimonas、消化链球菌属、埃希志贺氏杆菌属、Flavonifractor、Erysipelatoclostridium、链球菌属、Peptoniphilus、肠球菌属中的至少一种；溃疡性结肠炎患者中缺失的菌属分类单元为考拉杆菌属、瘤胃菌科_UCG-002、瘤胃球菌属_2、粪球菌属_2、Akkermansia、Odoribacter、[Eubacterium]_ventriosum_group、毛螺菌科_NK4A136_group、瘤胃菌科_UCG-014、瘤胃菌科_UCG-013、un_f_瘤胃菌科、Fusicatenibacter、Ruminiclostridium_6、Parabacteroides、[Eubacterium]_xylanophilum_group、氏菌属、Parasutterella、Christensenellaceae_R-7_group、Alloprevotella、毛螺菌科_UCG-008、[Eubacterium]_hallii_group、Enterorhabdus、Senegalimassilia、un_o_Gastranaerophilales、Tyzzerella_3中的至少一种。通过上述数据，使得菌属分类单元数据信息构成更多样化，保证供受体配型模型的构成更合理可靠。
优选地，溃疡性结肠炎患者中富集的通路分类单元为PWY-6470、HEMESYN2-PWY、GALACTARDEG-PWY、METHGLYUT-PWY、PWY-6891、PWY-5910、ORNDEG-PWY、3-HYDPHENACETATE-PWY、KDO-NAGLIPASYN-PWY、NAD-BIOSYNTHESIS-II、PPGPPMET-PWY、RUMP-PWY、PWY-7377中的至少一种；溃疡性结肠炎患者中缺失的通路分类单元为PWY-7374、PWY-6151、PWY-7090、PWY-6572、PWY-7111、COBALSYN-PWY、PWY-6478、PWY-5005、GLUTORN-PWY、PWY4FS-7、PWY-5103、PWY-5101、PANTOSYN-PWY、PWY-6507、ILEUSYN-PWY、PWY-3001、PWY-6385、PWY-5659、PWY-5676的至少一种。通过上述数据，使得通路分类单元数据信息构成更多样化，保证供受体配型模型的构成更合理可靠。
本发明提供如下具体的实施例以阐述本方法的应用原理，有如下步骤组成：
1)全局权重设计：
根据专家打分建立判断矩阵，从而计算全局权重；
因此由判断矩阵M计算出6个特征指标的全局权重分别为0.072，0.433，0.210，0.037，0.210，0.037；
2)模型构建：
利用AHP进行模型构建，采用上述适合的理论方法与结论，优化参数；
3)指标计算：
对某个特定的受体运用求同存异的方法计算其供体独有的多样性指标特征权重为：
Walpha-diver＝(0.043，0.055，0.052，0.060，0.070，0.055，0.040，0.040，0.071，0.067，0.060，0.081，0.044，0.045，0.052，0.068，0.096)T
利用增益减害原则计算其菌属分类单元的特征权重分别为：Wbt＝(0.067，0.040，0.040，0.053，0.054，0.040，0.054，0.041，0.079，0.053，0.026，0.082，0.070，0.053，0.076，0.096，0.077)T
Wht＝(0.026，0.039，0.036，0.051，0.089，0.073，0.090，0.019，0.046，0.083，0.089，0.039，0.068，0.055，0.064，0.056，0.078)T
利用增益减害原则计算其通路分类单元的特征权重分别为：
Wbp＝(0.052，0.055，0.055，0.058，0.064，0.053，0.063，0.056，0.061.0.061，0.059，0.059，0.061，0.059，0.061，0.061，0.062)T
Whp＝(0.058，0.059，0.059，0.058，0.059，0.059，0.059，0.058，0.059，0.059，0.059，0.059，0.059，0.059，0.059，0.059，0.059)T
最后运用bray-curtis方法计算其距离的特征权重为：
Wdist＝(0.050，0.056，0.057，0.057，0.062，0.062，0.061，0.052，0.060，0.062，0.063，0.058，0.061，0.059，0.060，0.058，0.061)T
4)性能评定：
将根据1)、2)和3)确定的模型计算得分，公式为：
ScoreDonor＝GW[1]*Wdist+GW[2]*Walpha-diver+GW[3]*Wbt+GW[4]*Wht+GW[5]*Wbp+GW[6]*Whp
将上述供受体配型模型运用于已公开的临床研究的21个受体和17个供体，与临床缓解率相结合评定其准确率达到62％。
本发明是通过优选实施例进行描述的，本领域技术人员知悉，在不脱离本发明的精神和范围的情况下，可以对这些特征和实施例进行各种改变或等效替换。本发明不受此处所公开的具体实施例的限制，其他落入本申请的权利要求内的实施例都属于本发明保护的范围。
Claims (6)
1.一种肠道菌群移植治疗溃疡性结肠炎的供受体配型方法，其特征在于，包括如下步骤：
（S1）：收集临床溃疡性结肠炎患者受体以及健康人员供体的肠道菌群数据，建立数据库；
（S2）：通过分析处理数据库数据，提取供体独有多样性、菌属分类单元、通路分类单元以及用以比较供体与受体之间组成差异的供受体距离数据组成特征，构建溃疡性结肠炎供受体配型模型，并对各指标特征进行计算分析；
所述步骤（S2）中，采用求同存异法，根据供体和受体肠道内多样性的种类进行比对，找出供体独有多样性特征；
所述求同存异法计算供体独有多样性的公式为：
式中：
代表供体独有的多样性种类个数；
代表供体实际的多样性种类个数；
代表供体和受体共有的多样性种类个数；
代表供体独有多样性种类的相对丰度；
所述步骤（S2）中，在溃疡性结肠炎患者中选择至少一种富集或缺失的菌属分类单元，组成菌属分类特征；
利用增益减害原则计算菌属分类单元，其公式为：
式中，
代表受体缺少的有益菌属分类单元但在供体中出现的有益菌属分类单元个数在受体缺少的有益菌属分类单元总数的比例；
代表在供体和受体中均出现的有害菌属分类单元的个数在受体出现的有害菌属分类单元总数的比例；
和
分别代表对应的相对丰度；
所述步骤（S2）中，在溃疡性结肠炎患者中选择至少一种富集或缺失的通路分类单元，组成通路分类特征；
利用增益减害原则计算通路分类单元，其公式为：
式中，
代表受体缺少的有益通路分类单元但在供体中出现的有益通路分类单元个数在受体缺少的有益通路分类单元总数的比例；
代表在供体和受体中均出现的有害通路分类单元的个数在受体出现的有害通路分类单元总数的比例；
和
分别代表对应的相对丰度；
所述供受体距离计算公式为：
其中，
、
分别代表受体和供体第i个多样性种类的个数；
（S3）：所述配型模型通过层次分析法分解特征层次结构，对供体库中供体进行打分排序，计算全局权重，选择最优值供体，使患者获得最优供体配型。
2.根据权利要求1所述的肠道菌群移植治疗溃疡性结肠炎的供受体配型方法，其特征在于：
所述供体和受体肠道内多样性的种类是基于所述数据库中供体和受体的粪便样本进行16S序列分析得到的。
3.根据权利要求1所述的肠道菌群移植治疗溃疡性结肠炎的供受体配型方法，其特征在于：
在所述步骤（S3）中，所述层次分析法进行供受体配型时包括如下步骤：
（S31）：将溃疡性结肠炎患者供受体配型涉及的因素划分为目标层、准则层以及方案层，建立供受体配型层次结构模型，其中，所述准则层包括供体独有多样性、有益菌属分类单元、有害菌属分类单元、有益通路分类单元、有害通路分类单元以及距离；
（S32）：构造两两比较判断矩阵，所述判断矩阵用数值表示下层元素对上层元素的相对重要性；
（S33）：进行层次单排序，计算单一准则下元素的相对权值，通过随机一致性比率CR来检验判断矩阵的一致性，当CR＜0.1时，认为判断矩阵具有满意的一致性，当CR≥0.1时，则需重新调整判断矩阵，使其具有一致性；
（S34）：计算各层元素对目标层的总排序权重，检验一致性；
（S35）：根据综合权重得分进行排序，选择权重高的供体作为最优供体。
4.根据权利要求1所述的肠道菌群移植治疗溃疡性结肠炎的供受体配型方法，其特征在于，还包括如下步骤：
（S4）：将根据所述步骤（S1）、所述步骤（S2）以及所述步骤（S3）确定的供受体配型模型运用于临床研究，与临床缓解率相结合评定其准确率。
5.根据权利要求1所述的肠道菌群移植治疗溃疡性结肠炎的供受体配型方法，其特征在于：
所述溃疡性结肠炎患者中富集的菌属分类单元为Sellimonas、Parvimonas、消化链球菌属、埃希志贺氏杆菌属、Flavonifractor、Erysipelatoclostridium、链球菌属、Peptoniphilus、肠球菌属中的至少一种；
所述溃疡性结肠炎患者中缺失的菌属分类单元为考拉杆菌属、瘤胃菌科_UCG-002、瘤胃球菌属_2、粪球菌属_2、Akkermansia、Odoribacter、[Eubacterium]_ventriosum_group、毛螺菌科_NK4A136_group、瘤胃菌科_UCG-014、瘤胃菌科_UCG-013、un_f_瘤胃菌科、Fusicatenibacter、Ruminiclostridium_6、Parabacteroides、[Eubacterium]_xylanophilum_group、氏菌属、Parasutterella、Christensenellaceae_R-7_group、Alloprevotella、毛螺菌科_UCG-008、[Eubacterium]_hallii_group、Enterorhabdus、Senegalimassilia、un_o_Gastranaerophilales、Tyzzerella_3中的至少一种。
6.根据权利要求1所述的肠道菌群移植治疗溃疡性结肠炎的供受体配型方法，其特征在于：
所述溃疡性结肠炎患者中富集的通路分类单元为PWY-6470、HEMESYN2-PWY、GALACTARDEG-PWY、METHGLYUT-PWY、PWY-6891、PWY-5910、ORNDEG-PWY、3-HYDPHENACETATE-PWY、KDO-NAGLIPASYN-PWY、NAD-BIOSYNTHESIS-II、PPGPPMET-PWY、RUMP-PWY、PWY-7377中的至少一种；
所述溃疡性结肠炎患者中缺失的通路分类单元为PWY-7374、PWY-6151、PWY-7090、PWY-6572、PWY-7111、COBALSYN-PWY、PWY-6478、PWY-5005、GLUTORN-PWY、PWY4FS-7、PWY-5103、PWY-5101、PANTOSYN-PWY、PWY-6507、ILEUSYN-PWY、PWY-3001、PWY-6385、PWY-5659、PWY-5676的至少一种。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011527189.1A CN112614596B (zh) | 2020-12-22 | 2020-12-22 | 一种肠道菌群移植治疗溃疡性结肠炎的供受体配型方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011527189.1A CN112614596B (zh) | 2020-12-22 | 2020-12-22 | 一种肠道菌群移植治疗溃疡性结肠炎的供受体配型方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN112614596A CN112614596A (zh) | 2021-04-06 |
| CN112614596B true CN112614596B (zh) | 2023-01-10 |
Family
ID=75244057
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202011527189.1A Active CN112614596B (zh) | 2020-12-22 | 2020-12-22 | 一种肠道菌群移植治疗溃疡性结肠炎的供受体配型方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN112614596B (zh) |
Families Citing this family (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN114446396A (zh) * | 2021-12-17 | 2022-05-06 | 广州保量医疗科技有限公司 | 肠道菌群移植的群组配型方法、系统、设备及存储介质 |
| CN114496094A (zh) * | 2021-12-17 | 2022-05-13 | 广州保量医疗科技有限公司 | 肠道菌群移植配型的度量学习方法、系统、设备及介质 |
| CN114496279B (zh) * | 2022-01-12 | 2022-08-30 | 广州保量医疗科技有限公司 | 菌群移植配型的排序方法、系统、计算机设备及存储介质 |
| CN114496277B (zh) * | 2022-01-12 | 2022-07-26 | 广州保量医疗科技有限公司 | 肠道菌群配型的多组学数据优化方法、系统、设备及介质 |
| CN114496278B (zh) * | 2022-01-12 | 2022-09-27 | 广州保量医疗科技有限公司 | 肠道菌群移植配型的数据处理方法、系统、设备及介质 |
| CN114794016B (zh) * | 2022-05-09 | 2023-07-18 | 承葛生物技术(广州)有限公司 | 一种构建肠道菌群分布紊乱和抗肿瘤免疫能力失调模型的方法 |
| CN117312893B (zh) * | 2023-03-28 | 2024-04-09 | 广东弘元普康医疗科技有限公司 | 一种菌群匹配度的评估方法及相关装置 |
| CN117037914B (zh) * | 2023-07-25 | 2024-04-02 | 中科国康(浙江)生命科学有限公司 | 一种基于个体特征定向制备的肠道菌群的系统 |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106661109A (zh) * | 2014-04-30 | 2017-05-10 | 马克思-德布鲁克-分子医学中心亥姆霍兹联合会 | 抗cd269(bcma)人源化抗体 |
| CN107708811A (zh) * | 2015-04-21 | 2018-02-16 | 恩立夫克治疗有限责任公司 | 治疗性汇集的血液凋亡细胞制剂与其用途 |
| CN110176305A (zh) * | 2019-05-27 | 2019-08-27 | 天益健康科学研究院（镇江）有限公司 | 一种利用高通量基因测序评估肠道菌群健康的方法 |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20200265940A1 (en) * | 2019-02-15 | 2020-08-20 | The Regents Of The University Of California | Predictive treatment controller for vedolizumab |
-
2020
- 2020-12-22 CN CN202011527189.1A patent/CN112614596B/zh active Active
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106661109A (zh) * | 2014-04-30 | 2017-05-10 | 马克思-德布鲁克-分子医学中心亥姆霍兹联合会 | 抗cd269(bcma)人源化抗体 |
| CN107708811A (zh) * | 2015-04-21 | 2018-02-16 | 恩立夫克治疗有限责任公司 | 治疗性汇集的血液凋亡细胞制剂与其用途 |
| CN110176305A (zh) * | 2019-05-27 | 2019-08-27 | 天益健康科学研究院（镇江）有限公司 | 一种利用高通量基因测序评估肠道菌群健康的方法 |
Non-Patent Citations (1)
| Title |
|---|
| "肠菌移植规范化管理中的伦理思考";刘贇鹏，许鸿志，肖传兴，杨晓宁，任建林;《中国医学伦理学》;20170831(第8期);948-952 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN112614596A (zh) | 2021-04-06 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN112614596B (zh) | 一种肠道菌群移植治疗溃疡性结肠炎的供受体配型方法 | |
| Clark et al. | Strongyle infection and gut microbiota: profiling of resistant and susceptible horses over a grazing season | |
| Dettman et al. | Reproductive isolation and phylogenetic divergence in Neurospora: comparing methods of species recognition in a model eukaryote | |
| Gandelman et al. | Machine learning reveals chronic graft-versus-host disease phenotypes and stratifies survival after stem cell transplant for hematologic malignancies | |
| CN105296590B (zh) | 大肠癌标志物及其应用 | |
| Pulkstenis et al. | Two goodness‐of‐fit tests for logistic regression models with continuous covariates | |
| Sahri et al. | Towards a comprehensive characterization of durum wheat landraces in Moroccan traditional agrosystems: analysing genetic diversity in the light of geography, farmers’ taxonomy and tetraploid wheat domestication history | |
| CN110097928B (zh) | 一种基于肠道菌群预测组织微量元素含量的预测方法和预测模型 | |
| CN108064272B (zh) | 用于类风湿性关节炎的生物标记物及其用途 | |
| De Luca et al. | Conservation and genetic characterisation of common bean landraces from Cilento region (southern Italy): high differentiation in spite of low genetic diversity | |
| CN112435756B (zh) | 基于多数据集差异互证的肠道菌群关联疾病风险预测系统 | |
| CN107194137A (zh) | 一种基于医疗数据建模的坏死性小肠结肠炎分类预测方法 | |
| Wang et al. | The nature of interspecific interactions and co‐diversification patterns, as illustrated by the fig microcosm | |
| CN113380396A (zh) | 一种基于粪便微生物标志物和人dna含量的多种肠道疾病风险评估的方法及应用 | |
| CN106156541A (zh) | 分析个体两类状态的免疫差异的方法和装置 | |
| CN112951423B (zh) | 培美曲塞化疗后不良反应发生风险预测模型及其构建方法 | |
| CN109949866A (zh) | 病原体操作组的检测方法、装置、计算机设备和存储介质 | |
| Roquet et al. | Evolutionary origins and species delineation of the two Pyrenean endemics Campanula jaubertiana and C. andorrana (Campanulaceae): evidence for transverse alpine speciation | |
| Cloete et al. | Genetic analysis of faecal worm egg count in South African Merinos under natural challenge | |
| CN115831375A (zh) | 一种抗乳腺癌候选药物的效果预测模型构造方法 | |
| Trinajstic et al. | Phenotypic plasticity, polymorphism and phylogeny within placoderms | |
| CN116364179A (zh) | 结直肠癌预后标志物筛选系统及方法、结直肠癌预后风险评估系统 | |
| Pfohl et al. | Association between post-dexamethasone cortisol level and blood pressure in depressed inpatients | |
| Hall et al. | The effects of coffee consumption and antibiotic use on gut microbial community structure of Parkinson’s disease patients | |
| Bergmans et al. | Competition between the sibling species Tisbe holothuriae and T. battagliai (Harpacticoida) |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant | ||
| EE01 | Entry into force of recordation of patent licensing contract |
Application publication date: 20210406 Assignee: Ningde Chengge Pharmaceutical Co.,Ltd. Assignor: Xiamen chengge Biotechnology Co.,Ltd. Contract record no.: X2024980004686 Denomination of invention: A donor receptor matching method for intestinal microbiota transplantation in the treatment of ulcerative colitis Granted publication date: 20230110 License type: Common License Record date: 20240424 |
|
| EE01 | Entry into force of recordation of patent licensing contract |