类风湿关节炎的诊断及预后标志物的制作方法

1.本技术涉及类风湿关节炎诊断技术领域，尤其是涉及类风湿关节炎的诊断及预后标志物。


背景技术：

2.类风湿关节炎(rheumatoid arthritis，ra)是一种以周身多关节疼痛、肿胀及进行性侵蚀病变为特点的全身自身免疫性疾病，治疗不及时可致关节畸形、功能丧失，最终会导致残疾、生活质量降低甚至死亡，时常还伴有包括心血管和肺等组织器官受累病变。类风湿关节炎的治疗措施包括一般性治疗、药物治疗和外科手术等。但是，目前的治疗措施通常仅能控制病情而无法治愈，因此需要患者进行定期随访，及时评估以根据实际情况考虑是否改变治疗方案。然而，现有的评估手段是基于观察和问卷评分的acr标准(如acr70)，相对于生理学改变，无法及时准确反映患者的状态，所以有必要从ra的发病机制出发，寻找合适的诊断或预后的标志物。
3.关于ra的发病机制目前尚不明确，目前的主流观点认为是遗传、环境与免疫之间相互作用的结果。ra在遗传因素方面的研究已经取得了一定的进展，但仅仅依靠遗传因素仍不足以解释ra的发病机制。部分研究结果显示，ra患者在关节炎症状发生前就能在血液中检测到类风湿因子和抗环胍氨酸肽抗体等自身免疫性抗体。该结果表明自身免疫的发生部位可能起源于关节外，也提示研究人员有必要关注遗传学指标以外的生物标志物。
4.随着微生物dna测序技术的不断发展，微生物与自身免疫性疾病的研究得到了广泛的关注和重视，肠道菌群的稳定是维持机体正常新陈代谢和免疫防御体系的重要保证，研究发现肠道菌群失调与疾病发生有着密切的关系。肠道菌群失调会改变上皮和粘膜通透性，使机体丧失免疫耐受功能，致敏免疫细胞和抗原在关节处聚集，导致类风湿关节炎等疾病的发生。但肠道菌群失调与ra疾病发生之间的因果关系仍不清楚，在疾病发生发展过程中起关键作用的肠道菌群的筛选、鉴定及功能作用有待进一步研究。因此，有必要从微生物组学出发，寻找具有良好诊断价值的ra代谢标志物。


技术实现要素：

5.本技术旨在至少解决现有技术中存在的技术问题之一。为此，本技术提出一种具有良好的诊断价值的类风湿关节炎的代谢相关的标志物，定量检测这些标志物的试剂能够用于类风湿关节炎的诊断或预后判断。
6.本技术的第一方面，提供定量检测肠道菌群的标志物的试剂在制备用于类风湿关节炎的诊断和/或预后的产品中的应用，标志物选自以下a和b中四种中的至少一种：
7.a.代谢物：9,12-十八碳二炔酸、10z壬烯酸；
8.b.微生物：疣微菌门(verrucomicrobia)、阿克曼菌属(akkermansia)。
9.根据本技术实施例的应用，至少具有如下有益效果：
10.申请人在实验中发现，类风湿关节炎患者肠道中9,12-十八碳二炔酸、10z壬烯酸
的代谢水平以及疣微菌门(verrucomicrobia)和阿克曼菌属(akkermansia)的微生物含量相比于正常人相比存在显著的上调，通过对这些代谢物或微生物中的至少一种的水平的定量检测，可以有效地区分类风湿关节炎患者和健康的正常人群，具有良好的诊断和预后价值。
11.在本技术的一些实施方式中，标志物来源于粪便样本。
12.在本技术的一些实施方式中，试剂盒还包括标志物的提取试剂。通过提取试剂将肠道菌群代谢物或微生物从粪便样本中富集出来，提取试剂可以是本领域任选的能够有效提取标志物的试剂。可以理解的是，为了不影响后续检测试剂的测定，对提取试剂的成分做一定调节。优选的，代谢物的提取试剂包括乙腈、甲醇和水，进一步三者体积比例优选为(1～5)：(1～5)：(1～3)，优选为2：2：1。微生物的提取试剂包括微生物核酸的抽提试剂。
13.在本技术的一些实施方式中，试剂通过核磁共振法、色谱法、光谱法、质谱法、pcr、测序、抗体-抗原反应中的至少一种对标志物进行定量检测。
14.在本技术的一些实施方式中，试剂通过核磁共振法、色谱法、光谱法、质谱法、抗体-抗原反应中的至少一种对代谢物进行定量检测。
15.在本技术的一些实施方式中，试剂通过pcr、测序、抗体-抗原反应中的至少一种对微生物进行定量检测。
16.其中，核磁共振法包括但不限于一维核磁共振法、二维核磁共振法等，色谱法包括但不限于气相色谱法、液相色谱法、高效液相色谱法、超高效液相色谱法等，光谱法包括但不限于紫外可见吸收光谱法、红外吸收光谱法、近红外光谱法等，质谱法包括但不限于串联质谱法、飞行时间质谱法等。而上述多种方法的组合包括但不限于液相色谱-质谱联用法。
17.pcr或测序的试剂包括扩增肠道微生物群16s rrna的试剂，优选的，扩增其v3～v4区。优选的，试剂包括疣微菌门和/或阿克曼菌属细菌特异性保守区序列的特异性扩增引物、特异性识别探针等。
18.抗体-抗原反应的试剂通过特异性抗体检测样本中代谢物或微生物的含量。
19.本技术的第二方面，提供类风湿关节炎的诊断和/或预后的试剂盒，该试剂盒包括定量检测肠道菌群的标志物的试剂，标志物选自以下a和b中四种中的至少一种：
20.a.代谢物：9,12-十八碳二炔酸、10z壬烯酸；
21.b.微生物：疣微菌门(verrucomicrobia)、阿克曼菌属(akkermansia)。
22.在本技术的一些实施方式中，试剂盒为粪便检测试剂盒。
23.在本技术的一些实施方式中，试剂盒还包括标志物的提取试剂。
24.本技术的第三方面，提供计算机可读存储介质，计算机可读存储介质存储有计算机可执行指令，计算机可执行指令用于使计算机执行以下操作：
25.步骤1：获取来自受试者的样本中的标志物的水平，标志物选自以下a和b中四种中的至少一种：
26.a.代谢物：9,12-十八碳二炔酸、10z壬烯酸；
27.b.微生物：疣微菌门(verrucomicrobia)、阿克曼菌属(akkermansia)；
28.步骤2：将标志物的水平与设定值比较，根据比较结果对受试者进行类风湿关节炎的诊断或预后判断。
29.在本技术的一些实施方式中，步骤1还包括将标志物的水平进行标准化。通过标准
化处理以进一步避免可能引起的诊断结果误差。
30.在本技术的一些实施方式中，可以获取至少两种标志物的水平，将这些水平进行数学关联以获得评分，将评分与设定的阈值评分进行比较，根据比较结果对受试者进行类风湿关节炎的诊断或预后判断。数学关联的方法可以是对每一种标志物赋予一个参数向量，得到评分与标志物水平的线性回归模型。可以理解的是，线性回归模型中各个参数向量可以通过训练以获得更好的取值以达到更准确的诊断或预后作用。
31.在本技术的一些实施方式中，设定值为用作对照的健康人的对应标志物的水平。
32.在本技术的一些实施方式中，当受试者的9,12-十八碳二炔酸和/或10z壬烯酸和/或疣微菌门(verrucomicrobia)和/或阿克曼菌属(akkermansia)的水平为健康人的相应水平的2倍以上时，受试者可以被诊断为患有类风湿关节炎，或预后不良；优选的，在2.5倍、3倍、4倍、5倍以上时，受试者可以被诊断为患有类风湿关节炎，或预后不良。
33.本技术的第四方面，提供电子设备，该电子设备包括处理器和存储器，存储器上存储有可在处理器上运行的计算机程序，处理器在运行计算机程序时实现以下操作：
34.步骤1：获取来自受试者的样本中的标志物的水平，标志物选自以下a和b中四种中的至少一种：
35.a.代谢物：9,12-十八碳二炔酸、10z壬烯酸；
36.b.微生物：疣微菌门(verrucomicrobia)、阿克曼菌属(akkermansia)；
37.步骤2：将标志物的水平与设定值比较，根据比较结果对受试者进行类风湿关节炎的诊断或预后判断。
38.其中，存储器作为一种非暂态计算机可读存储介质，可用于存储非暂态软件程序以及非暂态性计算机可执行程序，如前述的对受试者进行类风湿关节炎的诊断或预后判断。处理器通过运行存储在存储器中的非暂态软件程序以及指令，从而实现上述对受试者进行类风湿关节炎的诊断或预后判断。
39.存储器可以包括存储程序区和存储数据区，其中，存储程序区可存储操作系统、至少一个功能所需要的应用程序；存储数据区可存储执行上述步骤。此外，存储器可以包括高速随机存取存储器，还可以包括非暂态存储器，比如至少一个磁盘存储器件、闪存器件、或其他非暂态固态存储器件。例如，存储器可选包括相对于处理器远程设置的存储器，这些远程存储器可以通过网络连接至该处理器。上述网络的实例包括但不限于互联网、企业内部网、局域网、移动通信网及其组合。
40.实现上述步骤所需的非暂态软件程序以及指令存储在存储器中，当被一个或者多个处理器执行时，执行上述步骤。
41.以上所描述的装置实施例仅仅是示意性的，其中作为分离部件说明的单元可以是或者也可以不是物理上分开的，即可以位于一个地方，或者也可以分布到多个网络单元上。可以根据实际的需要选择其中的部分或者全部模块来实现本实施例方案的目的。
42.另外，上文中所公开的系统可以被实施为软件、固件、硬件及其适当的组合。某些物理组件或所有物理组件可以被实施为由处理器，如中央处理器、数字信号处理器或微处理器执行的软件，或者被实施为硬件，或者被实施为集成电路，如专用集成电路。这样的软件可以分布在计算机可读介质上，计算机可读介质可以包括计算机存储介质(或非暂时性介质)和通信介质(或暂时性介质)。如本领域普通技术人员公知的，术语计算机存储介质包
括在用于存储信息(诸如计算机可读指令、数据结构、程序模块或其他数据)的任何方法或技术中实施的易失性和非易失性、可移除和不可移除介质。计算机存储介质包括但不限于ram、rom、eeprom、闪存或其他存储器技术、cd-rom、数字多功能盘(dvd)或其他光盘存储、磁盒、磁带、磁盘存储或其他磁存储装置、或者可以用于存储期望的信息并且可以被计算机访问的任何其他的介质。此外，通信介质通常包含计算机可读指令、数据结构、程序模块或者诸如载波或其他传输机制之类的调制数据信号中的其他数据，并且可包括任何信息递送介质。
43.申请人在实验中发现，类风湿关节炎患者的肠道中9,12-十八碳二炔酸、10z壬烯酸的代谢水平以及疣微菌门(verrucomicrobia)、阿克曼菌属(akkermansia)的微生物的含量相比于正常人相比存在显著的上调，通过对这些标志物中的至少一种的水平的定量检测，可以有效地区分类风湿关节炎患者和健康的正常人群，具有良好的诊断和预后价值。
44.本技术的附加方面和优点将在下面的描述中部分给出，部分将从下面的描述中变得明显，或通过本技术的实践了解到。
附图说明
45.图1是本技术的实施例1中16s rrna测序在门分类水平下的比较结果。其中，a是hf组和ra组的不同个体所检出的不同肠道微生物的堆积柱形图；b是a中各个菌门的肠道微生物在两组样本中检出的比例；c是b中各个菌门在两组样本中检出的比例的堆积柱形图；d是在门分类水平下，hf组和ra组样本中存在显著差异的四个菌门的检验结果。
46.图2是本技术的实施例1中16s rrna测序在属分类水平下的比较结果。其中，a是hf组和ra组的不同个体所检出的不同肠道微生物的堆积柱形图；b是a中各个菌门的肠道微生物在两组样本中检出的比例；c是b中各个菌门在两组样本中检出的比例的堆积柱形图。
47.图3是本技术的实施例1中采用lefse方法筛选得到的不同分类水平下的差异表达的生物标志物，其中，a是筛选结果的分支图，从内侧到外侧的圆圈表示从门到种的分类水平，每个小圆圈的直径与肠道微生物群的相对丰度成正比；b是筛选出的显著差异的肠道菌群成分的柱状图。
48.图4是本技术的实施例1中采用lc-ms方法筛选得到的在hf组和ra组中存在显著差异的代谢物。
49.图5是本技术的实施例1中差异代谢通路富集以及差异代谢物诊断价值的结果。其中，a是代谢通路富集的气泡图；b是差异代谢物与不同分类水平的差异微生物的spearman相关分析结果；c和d分别是9,12-十八碳二炔酸和10z壬烯酸的roc曲线；e和f分别是hf组和ra组样本中9,12-十八碳二炔酸和10z壬烯酸的差异表达情况。
具体实施方式
50.以下将结合实施例对本技术的构思及产生的技术效果进行清楚、完整地描述，以充分地理解本技术的目的、特征和效果。显然，所描述的实施例只是本技术的一部分实施例，而不是全部实施例，基于本技术的实施例，本领域的技术人员在不付出创造性劳动的前提下所获得的其他实施例，均属于本技术保护的范围。
51.下面详细描述本技术的实施例，描述的实施例是示例性的，仅用于解释本技术，而
不能理解为对本技术的限制。
52.在本技术的描述中，若干的含义是一个以上，多个的含义是两个以上，大于、小于、超过等理解为不包括本数，以上、以下、以内等理解为包括本数。如果有描述到第一、第二只是用于区分技术特征为目的，而不能理解为指示或暗示相对重要性或者隐含指明所指示的技术特征的数量或者隐含指明所指示的技术特征的先后关系。
53.本技术的描述中，参考术语“一个实施例”、“一些实施例”、“示意性实施例”、“示例”、“具体示例”、或“一些示例”等的描述意指结合该实施例或示例描述的具体特征、结构、材料或者特点包含于本技术的至少一个实施例或示例中。在本说明书中，对上述术语的示意性表述不一定指的是相同的实施例或示例。而且，描述的具体特征、结构、材料或者特点可以在任何的一个或多个实施例或示例中以合适的方式结合。
54.实施例1
55.实验方法
56.1.肠道微生物的16s rrna测序
57.分别收集30例健康人和29例初发未治疗的ra患者的粪便标本分别作为对照组(hf)和疾病组(ra)，使用dna isolation kit提取试剂盒提取肠道菌群的总基因组d na。用正向引物5'-actcctacggaggcagca-3'(seq id no.1)和反向引物5'-ggacta chggggtwtctaat-3'(seq id no.2)进行pcr扩增肠道微生物群16s rrna基因的v3-v4区。用vahtstm dna磁珠纯化pcr产物，然后将得到的pcr产物作为dna模板进行solexa pcr扩增。纯化后，用nanodrop 2000测定dna浓度，用1.8％琼脂糖凝胶电泳(120v，40min)回收dna片段。基于illumina-hiseq-2500测序平台，采用双端法构建小片段文库并测序。
58.使用flash软件对每个样本的reads进行拼接，得到原始tags数据。然后用trimmomatic软件对原始tags数据进行过滤，得到高质量的tags数据。利用uchime软件去除嵌合体序列，得到有效的数据。然后用qiime软件对数据进行聚类，相似度大于97％的otu被认为是同一个otu。然后，将每个otu序列的物种与微生物参考数据库获取每个otu分类的物种信息。α多样性分析用于分析肠道微生物群落的物种丰富度(chao1和ace)和多样性(shannon和simpson指数)，采用β多样性分析法，包括主坐标分析法(pcoa)，用r语言分析不同样品肠道菌群组成的差异。用lefse(linear discriminant analysis effect size)法分析两组之间存在显著差异的生物标志物。
59.2.非靶标代谢组学检测
60.称取50mg粪便样品于ep管中，加入1000μl提取溶剂(乙腈-甲醇-水，2:2:1，含内标物1μg/ml)后，涡旋混匀30s，加入钢珠，在45hz下研磨处理4min，在冰水浴中超声5min。将匀浆和超声循环重复3次，然后在-20℃下静置1小时，并在12000rpm和4℃条件下离心15分钟。所得上清液转移到进样瓶中进行检测。质控品由所有样品中混合等量的上清液制得。
61.lc-ms/ms分析使用uhplc系统(1290，agilent technologies)，uplc hss t3柱(2.1mm
×
100mm，1.8μm)与q exactive(orbitrap ms，thermo)耦合。流动相a为0.1％甲酸 5mmol/l醋酸铵的水溶液，流动相b为乙腈。洗脱梯度为：0min，1％b；1分钟，1％b；8min，99％b；10min，99％b；10.1分钟，1％b；12分钟，1％b。流速为0.5ml/min，进样量为2μl。在lc/ms实验中，qe质谱仪能够在信息依赖的基础上获得ms/ms光谱。在这种模式下，采集软件(xcalibur 4.0.27，thermo)在根据预选标准采集和触发ms/ms光谱采集时，不断评估全扫
描测量ms数据。esi源条件设定为：鞘层气体流量为45arb，辅助气体流量为15arb，毛细管温度为400℃，全ms分辨率为70000，ms/ms分辨率为17500，碰撞能量为20/40/60ev，喷射电压为4.0kv(正)或-3.6kv(负)。
62.使用proteowizard软件将质谱原始数据转成mzxml格式，再使用maps软件(xcms内核)做保留时间矫正、峰识别、峰提取、峰积分、峰对齐等工作。使用r程序包和二级质谱数据库对峰进行物质鉴定。利用simca软件对隐结构判别分析(opls-da)进行正交投影。对照组和疾病组之间的差异代谢产物根据vip得分的变量重要性和student’t检验分析的p值进行选择。通过内部ms2数据库和以及hmdb和pubchem在线数据库的比较，根据准确的分子量和片段模式鉴定代谢物。
63.实验结果
64.下面关于肠道菌群拉丁名的命名中，部分命名前端出现的单一字母用以区分其分类水平，其中，p、c、o、f、g和s分别代表门、纲、目、科、属和种的不同分类水平。
65.16s rrna测序结果如图1所示，其中，a是在门分类水平下，hf组和ra组的不同个体所检出的肠道微生物的具体成分和比例结果；b是a中各个菌门的肠道微生物在两组样本中检出的比例；c是b中各个菌门在两组样本中检出的比例的堆积柱形图；d是在门分类水平下，hf组和ra组样本中存在显著差异的四个菌门的检验结果。从a中可以看出，对照组和实验组中一共检测出10个门分类水平的肠道微生物菌群，其中，厚壁菌门(firmicutes)、拟杆菌门(bacteroidetes)、变形菌门(proteobacteria)和放线菌门(actinobacteria)是hf和ra组的主要成分。从b和c中可以看出，firmicutes是最主要的门分类水平的菌，分别占hf组和ra组物种的59.4％和59.1％。hf组中bacteroidetes的含量高于ra组(分别为23.0％和15.2％)，而proteobacteria(hf组和ra组的含量分别为8.6％和12.4％)和verrucomicrobia的含量低于ra组(8.6％/12.4％及6.7％/7.2％)。为了评估hf组和ra组的微生物群之间的差异，对其进行wilcoxon检验，结果如d所示，发现hf组和ra组的bacteroidetes(p＝0.0132)、cyanobacteria(p＝2.63e-08)、synergistetes(p＝0.0102)和verrucomicrobia(p＝0.0075)存在显著差异。
66.而在属分类水平下对hf组和ra组的肠道菌群进行筛选，共得到159个细菌属，参考图2，其中，a是在属分类水平下，hf组和ra组的不同个体所检出的不同肠道微生物的堆积柱形图；b是a中各个菌属的肠道微生物在两组样本中检出的比例；c是b中各个菌属在两组样本中检出的比例的堆积柱形图。如图2所示，faecalibacterium、bacteroides、escheri chia-shigella、blautia、prevotella_9、bifidobacterium、subdoligranulum、agathobacter、akke rmansia、streptococcus是hf组和ra组的主要成分。在159个细菌属中，其中有62个细菌属在ra组中含量减少，包括g_uncultured_bacterium_f_prevotellaceae、g_prevotella、g_faec alibacterium、g_bifidobacterium、g_bacteroides、g_ruminococcus_1、g_ruminococcaceae_uc g-002、g_lachnospiraceae_ucg-001和g_fusobacterium等；而97个细菌属在ra组的含量增加，包括g_lactobacillus、g_streptococcus、g_akkermansia、g_blautia和g_eggerthella等。另外，筛选到43个细菌属在hf组与ra组的含量存在显著差异(p＜0.05)。
67.采用lefse方法，以lda评分》4为筛选标准，筛选出hf组与ra组间差异显著的生物标志物，结果如图3所示。其中，a是分支图，从内侧到外侧的圆圈代表了从门到种的分类水
平，每个小圆圈的直径与肠道微生物群的相对丰度成正比；b是从中筛选的显著差异的肠道菌群成分丰度的柱状图。从图中可以看出，从门到种分类水平，hf组和ra组共检测出20种显著不同的肠道菌群成分。在这些成分中，p_verrucomicrobia、c_verrucomicrobiae、c_bacilli、o_verrucomicrobiales、o_lactobacillales、f_streptococcaceae、f_akkermansiaceae、g_blautia、g_akkermansia、g_streptococcus、s_uncultured_bacterium_g_blautia、s_uncultured_bacterium_g_akkermansia和s_uncultured_bacterium_g_streptococcus在ra组的含量较丰富，而p_bacteroidetes、c_bacteroidia、o_bacteroidales、f_ruminococcaceae、f_prevotellaceae、g_faecalibacterium和s_uncultured_bacterium_g_faecalibacterium在hf组中的含量较丰富。因此，以上筛选出的组间的差异菌群有可能成为ra疾病潜在的生物标志物。
68.16s rrna测序结果表明hf组和ra组的肠道菌群存在显著差异，推测肠道微生物的改变可能引起粪便代谢物的改变。进一步采用lc-ms分析，在负离子模式和正离子模式下分别获得279个和1045个代谢物。根据vip(正交偏最小二乘判别分析opls-da的第一主成分的变量权重)》1和p《0.05筛选条件，负离子模式中获得29个hf组和ra组的差异代谢物，其中3种代谢产物含量减少，26种代谢产物含量增加。在正离子模式下，hf组和ra组共获得45个差异代谢物，其中12种代谢产物含量减少，33种代谢产物含量增加。另外，ra组与hf组相比，部分肽类物质(如val-thr-ile、ile-gly-gly-ile、val-asn-ile等)、氨基酸类物质(如赖氨酸lysine、蛋氨酸methionine)和核苷酸(如胸腺嘧啶核苷thymidine、脱氧尿苷deoxyuridine、脱氧核苷deoxyinosine、脱氧鸟苷deoxyguanosine等)含量增加；而一些脂质(如lpa(0:0/16:0)、16,16-dimethyl-pga、1,5-hete、12(r)-hepe等)含量减少，如图4所示。
69.通过metaboanalyst网站分析上述差异代谢物富集的途径，结果如图5中的a所示。从图中可以看出，有两条代谢产物显著富集的途径：嘌呤代谢(purine metabolism，p＝0.01，impact value＝0.02)和组氨酸代谢(histidine metabolism，p＝0.02，impact value＝0.12)途径。结合图4，代谢产物黄嘌呤(xanthine)、脱氧核苷(deoxyinosine)、脱氧鸟苷(deoxyguanosine)和尿酸(uric acid)在嘌呤代谢通路中被富集出来，而代谢产物尿酸盐(urocanate)和4-咪唑乙酸(4-imidazolaetic acid)在组氨酸代谢通路中被富集出来。该结果表明，肠道微生物的改变可能引起了代谢途径的改变，尤其是核酸代谢和氨基酸代谢。
70.根据vip》1.5和fc(ra/hf)》2的筛选标准，得到38种丰度在hf组和ra组之间存在显著差异的代谢物。对显著差异的代谢物与显著差异的肠道菌群进行spearman相关分析，结果如图5的b所示，从图中可以看出，疣微菌门(p_verrucomicrobia)和阿克曼菌属(g_akkermansia)与9,12-十八碳二炔酸(r＝0.54，corrected p-value＝0.01)和10z壬烯酸(r＝0.54，corrected p-value＝0.03)显著正相关。值得注意的是，其中的两种差异代谢产物都是长链脂肪酸，该结果表明肠道微生物群的改变可能导致长链脂肪酸代谢紊乱，影响ra疾病的发生和发展。
71.根据上述两种代谢物的水平分别绘制roc曲线，结果如图5的c和d所示，从图中可以看出，9,12-十八碳二炔酸的roc曲线的下面积auc为0.718(p＝0.0043)，10z壬烯酸的roc曲线的下面积auc为0.739(p＝0.0064)，该结果表明采用这两个代谢物预测受试者是否患
有类风湿关节炎具有良好的特异性和和灵敏度。另外，hf组和ra组中这两个代谢物的相对丰度的结果见图5的e和f，从图中可以看出，这两个代谢物在ra组的含量显著高于h f组。另外，由于疣微菌门(verrucomicrobia)和阿克曼菌属(akkermansia)的微生物于上述两种代谢物存在显著的正相关，因而也可以通过对这两种微生物的含量来间接对受试者的类风湿关节炎的诊断或预后做出有效的判断。
72.实施例2
73.本实施例提供一种对受试者风湿性关节炎进行诊断的设备，该设备包括处理器和存储器，存储器上存储有可被处理器运行的计算机程序。运用该设备对受试者进行风湿性关节炎诊断的方法如下：
74.步骤1：获取受试者的粪便样本中9,12-十八碳二炔酸、10z壬烯酸、疣微菌门、阿克曼菌属中至少一种标志物的含量。
75.步骤2：将受试者的含量水平与健康人的设定值比较，如果高出3倍，判定受试者患有类风湿关节炎。
76.上面结合实施例对本技术作了详细说明，但是本技术不限于上述实施例，在所属技术领域普通技术人员所具备的知识范围内，还可以在不脱离本技术宗旨的前提下作出各种变化。此外，在不冲突的情况下，本技术的实施例及实施例中的特征可以相互组合。