通过多个miRNA的表达量确定雷公藤多苷片治疗类风湿性关节炎的个体有效性的系统

摘要

本公开提供了一种通过多个miRNA的表达量确定雷公藤多苷片治疗类风湿性关节炎的个体有效性的系统，所述多个miRNA包括hsa‑miR‑550b‑2‑5p、hsa‑miR‑4797‑5p、hsa‑miR‑6509‑5p和hsa‑miR‑378g。本公开还提供了一种确定雷公藤多苷片治疗类风湿性关节炎的个体有效性的分子标志物、该分子标志物在制备确定雷公藤多苷片治疗类风湿性关节炎的个体有效性的试剂盒中的用途、检测该分子标志物的表达量的试剂在制备确定雷公藤多苷片治疗类风湿性关节炎的个体有效性的试剂盒中的用途。本公开为类风湿性关节炎的个体化诊疗提供了高特异性、高敏感性的生物标志物和疗效预测模型。

说明书

技术领域

本公开涉及生物医学领域，具体地，涉及一种通过多个miRNA的表达量确定雷公藤多苷片治疗类风湿性关节炎的个体有效性的系统、一种确定雷公藤多苷片治疗类风湿性关节炎的个体有效性的分子标志物、该分子标志物在制备确定雷公藤多苷片治疗类风湿性关节炎的个体有效性的试剂盒中的用途、检测该分子标志物的表达量的试剂在制备确定雷公藤多苷片治疗类风湿性关节炎的个体有效性的试剂盒中的用途。

背景技术

类风湿性关节炎(Rheumatoid arthritis，RA)是世界难治性自身免疫性疾病，好发于青壮年，其发病率高(我国约为0.36％)、致残率高(5-10年病程时致残率约为60％)，对人类的健康危害极大，近年已被发达国家列为未来25年中影响人类健康最大、医疗消费最高的五大疾病之一。

雷公藤是我国的传统中药，最早记载于《滇南本草》，其味苦性寒，具清热解毒、祛风除湿、舒筋活血、消肿止痛功效，其成药制剂雷公藤多苷片一直被我国人民用来治疗风湿病，并已引起国际医药界的广泛关注，迄今为止，尚没有一味中药或中成药能取代雷公藤多苷片在抗类风湿性关节炎临床中的重要位置。然而，由于其存在临床疗效个体差异大、有效剂量与毒性剂量相当等问题，严重影响了临床的推广和合理应用。

本发明人近年对雷公藤及其活性成分治疗类风湿性关节炎的作用进行了较系统的研究，发现雷公藤及其活性成分具有缓解类风湿性关节炎病情、抑制关节滑膜炎症、防止或减缓关节结构破坏的改善病情抗风湿药样作用。同时，从类风湿性关节炎的关键病理环节入手，阐述了雷公藤最主要的免疫活性成分雷公藤甲素抑制关节滑膜炎症和血管新生，调节软骨基质降解和破骨细胞分化系统，抑制炎症相关信号通路活化、降低炎症反应级联和循环的相关作用机制。相关研究成果有助于了解雷公藤作用特点、推动其抗类风湿性关节炎的分子药理学发展。基于此，下一步的考虑是如何使研究成果更直接指导临床用药，怎样个体化地应用雷公藤多苷片使之达到最佳治疗效果，这也是目前国际风湿病研究领域的热点和亟待解决的难题。

因此，进一步筛选与雷公藤多苷片治疗类风湿性关节炎个体差异密切相关的新的生物标记物，对于揭示雷公藤多苷片治疗类风湿性关节炎的分子机制及提高其个体化用药水平均具有重要的意义。

发明内容

本公开的目的是提供与雷公藤多苷片治疗类风湿性关节炎个体差异密切相关的新的生物标记物，提高雷公藤多苷片治疗类风湿性关节炎的个体化用药水平，从而改善雷公藤多苷片治疗类风湿性关节炎的治疗效果。

为了实现上述目的，一方面，本公开提供了一种通过多个miRNA的表达量确定雷公藤多苷片治疗类风湿性关节炎的个体有效性的系统，其中，该系统包括计算装置、用于输入类风湿性关节炎患者个体的多个miRNA的表达量的输入装置和用于输出雷公藤多苷片治疗类风湿性关节炎患者的个体有效性的输出装置；其中，所述多个miRNA包括hsa-miR-550b-2-5p、hsa-miR-4797-5p、hsa-miR-6509-5p和hsa-miR-378g；所述计算装置包括存储器和处理器，所述存储器中存储有计算机程序，所述处理器被配置为执行所述存储器中存储的计算机程序，以实现建模算法和如式(1)所示的判别函数的算法；所述建模算法为支持向量机算法和/或最小偏二乘算法；

F(c)＝sgn[f₁(c₁)+f₂(c₂)+f₃(c₃)+f₄(c₄)+b]式(1)

式(1)中，F(c)表示雷公藤多苷片治疗类风湿性关节炎患者的个体有效性，F(c)返回值为1表示有效，返回值为-1时表示无效；c₁、c₂、c₃和c₄依次分别表示hsa-miR-550b-2-5p、hsa-miR-4797-5p、hsa-miR-6509-5p和hsa-miR-378g的相对表达量；所述相对表达量是指相对于内参的表达量比值；f₁(c₁)、f₂(c₂)、f₃(c₃)和f₄(c₄)分别为依据所述建模算法训练得到的核函数，b为依据所述建模算法训练得到的临界打分值。

另一方面，本公开还提供了一种确定雷公藤多苷片治疗类风湿性关节炎的个体有效性的分子标志物，其中，该分子标志物为多个miRNA，所述多个miRNA由hsa-miR-550b-2-5p、hsa-miR-4797-5p、hsa-miR-6509-5p和hsa-miR-378g组成。

另一方面，本公开还提供了分子标志物在制备确定雷公藤多苷片治疗类风湿性关节炎的个体有效性的试剂盒中的用途，其中，该分子标志物为多个miRNA，所述多个miRNA由hsa-miR-550b-2-5p、hsa-miR-4797-5p、hsa-miR-6509-5p和hsa-miR-378g组成。

再一方面，本公开还提供了检测分子标志物的表达量的试剂在制备确定雷公藤多苷片治疗类风湿性关节炎的个体有效性的试剂盒中的用途，其中，该分子标志物为多个miRNA，所述多个miRNA由hsa-miR-550b-2-5p、hsa-miR-4797-5p、hsa-miR-6509-5p和hsa-miR-378g组成。

类风湿性关节炎作为系统性免疫性疾病，伴随着基因、蛋白质以及分子相互作用网络的动态失调，具有总体上治愈率低、致残率高等特点。随着分子生物学领域的高速发展，科学家发现多种信号分子参与类风湿性关节炎的发生和发展，信号分子的活性受到遗传和环境中多种因素的影响，调节环路复杂的相互作用和交叉干扰常常使得人们对类风湿性关节炎相关的信号输出难以作出正确的认识和预测。鉴于类风湿性关节炎的高度复杂性，一直以来沿用的针对单一高特异性靶点的治疗策略并没有取得满意的效果。

在疾病诊疗生物标志物的研究领域，通常首先通过筛选疾病组织中的差异表达分子作为标志物，并采用人工神经网络、偏最小二乘法等构建疾病预测模型。但此种分析方法的可重复性较差。可见，进行疾病诊疗生物标志物的筛选只考虑差异表达数据是远远不够的。

本发明人从全局的角度出发，计算分析高通量生物分子表达谱数据，并整合生物分子网络的特征分析，筛选出了高可信度的疾病诊疗生物标志物及相应的计算公式，由此得到了本发明创造。具体地，首先，对类风湿性关节炎患者采用雷公藤多苷片治疗，通过DAS28、ACR20/50/70等疗效性指标划分雷公藤多苷片治疗达标/未达标组；接着，应用miRNA和基因表达谱芯片技术，检测两组的miRNA和基因的差异表达状况；再通过差异表达数据分析，选取差异表达miRNA和基因谱，并建立miRNA共表达调控网络和miRNA介导的基因表达调控网络；然后，筛选网络中的关键miRNAs作为评价雷公藤多苷片治疗类风湿性关节炎效果的生物标志物；最后，运用支持向量机算法和/或偏最小二乘法算法，针对上述生物标志物的表达量，构建雷公藤多苷片个体化治疗类风湿性关节炎的疗效预测模型，并对其性能进行评估。评估结果显示，本公开的技术方案为类风湿性关节炎的个体化诊疗提供了高特异性、高敏感性的生物标志物和疗效预测模型，并为临床制定类风湿性关节炎个体化治疗方案提供了一种新型、高效且无创的辅助工具，从而有效预测了类风湿性关节炎患者从临床治疗中的获益程度，从而可实现根据个体基因检测数据与药物效应的关系设计、制定个体化药物治疗方案，以充分发挥药物对机体的作用，增加首剂处方的有效性，提高治疗效果，减少药物不良反应的发生率，弥补传统或经验方案的不足，节省或合理配置医药卫生资源。

本公开的其他特征和优点将在随后的具体实施方式部分予以详细说明。

附图说明

附图是用来提供对本公开的进一步理解，并且构成说明书的一部分，与下面的具体实施方式一起用于解释本公开，但并不构成对本公开的限制。在附图中：

图1是实施例2和对比例2的ACC结果和AUC结果。

具体实施方式

以下结合附图对本公开的具体实施方式进行详细说明。应当理解的是，此处所描述的具体实施方式仅用于说明和解释本公开，并不用于限制本公开。

一方面，本公开提供了一种通过多个miRNA的表达量确定雷公藤多苷片治疗类风湿性关节炎的个体有效性的系统，其中，该系统包括计算装置、用于输入类风湿性关节炎患者个体的多个miRNA的表达量的输入装置和用于输出雷公藤多苷片治疗类风湿性关节炎患者的个体有效性的输出装置；其中，所述多个miRNA包括hsa-miR-550b-2-5p、hsa-miR-4797-5p、hsa-miR-6509-5p和hsa-miR-378g；所述计算装置包括存储器和处理器，所述存储器中存储有计算机程序，所述处理器被配置为执行所述存储器中存储的计算机程序，以实现建模算法和如式(1)所示的判别函数的算法；所述建模算法为支持向量机算法和/或最小偏二乘算法；

F(c)＝sgn[f₁(c₁)+f₂(c₂)+f₃(c₃)+f₄(c₄)+b]式(1)

式(1)中，F(c)表示雷公藤多苷片治疗类风湿性关节炎患者的个体有效性，F(c)返回值为1表示有效，返回值为-1时表示无效；c₁、c₂、c₃和c₄依次分别表示hsa-miR-550b-2-5p、hsa-miR-4797-5p、hsa-miR-6509-5p和hsa-miR-378g的相对表达量；所述相对表达量是指相对于内参的表达量比值；f₁(c₁)、f₂(c₂)、f₃(c₃)和f₄(c₄)分别为依据所述建模算法训练得到的核函数，b为依据所述建模算法训练得到的临界打分值。其中，支持向量机算法和最小偏二乘算法及其运算和训练的方式为已知的常规方式。

其中，hsa-miR-550b-2-5p的序列如SEQ ID NO.1所示，具体为AUGUGCCUGAGGGAGUAAGACA，hsa-miR-4797-5p的序列如SEQ ID NO.2所示，具体为GACAGAGUGCCACUUACUGAA，hsa-miR-6509-5p的序列如SEQ ID NO.3所示，具体为AUUAGGUAGUGGCAGUGGAAC，hsa-miR-378g的序列如SEQ ID NO.4所示具体为CACUGGGCUUGGAGUCAGAAGACCUGGCUCCAGCCCAGCUC。

可选地，该系统还包括多个miRNA的表达量的检测装置。

可选地，所述多个miRNA的表达量的检测装置包括miRNA表达量检测芯片和芯片信号读取器，miRNA表达量检测芯片包括分别检测hsa-miR-550b-2-5p、hsa-miR-4797-5p、hsa-miR-6509-5p和hsa-miR-378g的表达量的探针。或者所述多个miRNA的表达量的检测装置包括实时定量PCR仪和miRNA的实时定量PCR引物，miRNA的实时定量PCR引物包括分别检测hsa-miR-550b-2-5p、hsa-miR-4797-5p、hsa-miR-6509-5p和hsa-miR-378g的表达量的实时定量PCR引物。

可选地，所述miRNA表达量检测芯片还包括内参探针，所述内参探针为检测RNU6基因的表达量的探针，所述miRNA的实时定量PCR引物，还包括内参引物，所述内参引物为检测RNU6基因的实时定量PCR引物。

可选地，其中，作为训练得到的一组参考值，式(1)中，f₁(c₁)＝0.5554×c₁，f₂(c₂)＝0.4323×c₂，f₃(c₃)＝-0.5148×c₃，f₄(c₄)＝-0.4895×c₄，b＝-0.0016。也就是说，判别函数可以简化为如式(2)所示：

F(c)＝sgn(0.5554×c₁+0.4323×c₂-0.5148×c₃-0.4895×c₄-0.0016)式(2)。

需要说明的是，f₁(c₁)、f₂(c₂)、f₃(c₃)、f₄(c₄)和b可能会随着miRNA的表达量的检测手段的偏性而发生改变，也可能会随着训练数据集的数据规模大小等因素发生改变。式(2)所示的判别函数是本公开的发明人依据实施例1中的数据以最小偏二乘的建模算法进行训练得到的，并不限制本公开的范围。当然也可以根据其他的数据，和/或，选取支持向量机算法，选择hsa-miR-550b-2-5p、hsa-miR-4797-5p、hsa-miR-6509-5p和hsa-miR-378g的相对表达量，进行训练，得到式(1)范围内的判别函数。

另一方面，本公开还提供了一种确定雷公藤多苷片治疗类风湿性关节炎的个体有效性的分子标志物，其中，该分子标志物为多个miRNA，所述多个miRNA由hsa-miR-550b-2-5p、hsa-miR-4797-5p、hsa-miR-6509-5p和hsa-miR-378g组成。

另一方面，本公开还提供了分子标志物在制备确定雷公藤多苷片治疗类风湿性关节炎的个体有效性的试剂盒中的用途，其中，该分子标志物为多个miRNA，所述多个miRNA由hsa-miR-550b-2-5p、hsa-miR-4797-5p、hsa-miR-6509-5p和hsa-miR-378g组成。

再一方面，本公开还提供了检测分子标志物的表达量的试剂在制备确定雷公藤多苷片治疗类风湿性关节炎的个体有效性的试剂盒中的用途，其中，该分子标志物为多个miRNA，所述多个miRNA由hsa-miR-550b-2-5p、hsa-miR-4797-5p、hsa-miR-6509-5p和hsa-miR-378g组成。

其中，检测分子标志物的表达量的试剂可以为探针和/或引物。

以下，通过实施例进一步详细说明本发明。

实施例1

本实施例用于说明本公开的生物标志物的发现和预测模型的建立。

病例来源及样本量：中国中医科学院广安门医院风湿病科门诊及病房就诊的RA患者40例，符合纳入标准者。同时选取来自健康体检中心的健康志愿者10例，年龄和性别与疾病组相配。

纳入标准：符合美国风湿病学会(American College of Rheumatology,ACR)1987年RA的分类标准和2010类风湿关节炎分类标准(ACR/EULAR)，从未接受过雷公藤多苷片治疗或在最近4周内没有接受过雷公藤多苷片治疗的受试对象。

1987年美国风湿病学会修订的类风湿关节炎分类标准包括：

1.晨僵至少1小时(≥6周)；

2. 3个或3个以上关节区的关节炎(≥6周)；

3.腕、掌指关节或近端指间关节炎(≥6周)；

4.对称性关节炎(≥6周)；

5.皮下结节；

6.手X线改变；

7.类风湿因子阳性。

有上述七项中四项者即可诊断为类风湿关节炎。

2010ACR/EULAR类风湿关节炎分类标准包括：

A：受累关节

1个大关节(0分)

2-10个大关节(1分)

1-3个小关节(有或没有大关节)(2分)

4-10个小关节(有或没有大关节)(3分)

超过10个关节(至少一个小关节)(5分)

B：血清学(至少需要1项结果)

-RF和CCP(抗环瓜氨酸肽抗体)阴性(0分)

-RF和CCP，至少有一项是低滴度阳性(2分)

-RF和CCP，至少有一项高滴度阳性(3分)

C：急性期反应物(至少需要1项结果)

-CRP和ESR均正常(0分)

-CRP或ESR异常(1分)

D：症状持续时间

-＜6周(0分)

-≥6周(1分)

在A-D内，取病人符合条件的最高分。例如，患者有5个小关节和4个大关节受累，评分为3分。分数相加≥6分，诊断为类风湿关节炎。

3)排除标准：由于各种原因不能坚持服药12周的RA患者，患有其他风湿病的患者，育龄期女性或有生育要求的女性。

4)知情同意：所有研究对象均自愿签署知情同意书。知情同意书的内容包括本次研究的目的、意义和方法，研究对象因参加此项研究可能获得的益处和可能出现的风险，本研究的意义以及收集到的研究对象相关信息尤其个人隐私方面的保密问题等。

5)用药方案：雷公藤多苷片20mg，口服，每日3次。每两周随访一次，随访根据问卷调查雷公藤多苷片治疗期间的各项指标变化。

6)疗程：12周；

7)观察指标：采集患者的临床资料，包括性别、年龄、民族、发病时间、抗CCP水平、RF水平、

理化检查(ESR，CRP)及X光片/CT等影像学资料，并建立相应数据库。

8)疗效判定：以是否能满足达标治疗(Tight control)来判定疗效。所谓达标治疗，就是以降低RA患者病情活动度，达到临床缓解为目标的治疗方案。采用DAS28、ACR20和ACR50等病情评价方法。具体地说，患者使用雷公藤多苷片后，每两周评价1次，每次随访观察病情是否改善在20％以上，以及12周内是否达DAS28<2.4,如果达标，则表示治疗有效，否则为无效。

DAS28(Disease Activity Score in 28 joints，28个关节的平均疾病活动性评分)。

DAS28(4)＝0.56*sqrt(t28)+0.28*sqrt(sw28)+0.70*Ln(ESR)+0.014*GH。(t28：28个关节中的疼痛个数；sw28：28关节中的肿胀关节个数；ESR：血沉；GH：患者对疾病活动的整天评估(VAS评分，单位100mm)；28个关节包括：双手近端指间关节共10个、双手掌指关节共10个、双腕关节2个、双肘关节2个、双肩关节2个、双膝关节2个。)病情高活动度：≥5.1；病情低活动度：<3.2；病情缓解：<2.6。

ACR反应标准包括：关节压痛数改善程度和关节肿胀数改善程度。并包括下列5项中3项改善程度：

患者对疼痛的评价

患者对疾病活动的总体评价

医生对疾病活动的总体评价

患者对身体功能的评价(HAQ)

急性期反应物的数值(ESR、CRP)

各项指标的改善百分率＝(治疗前值－治疗后值)/治疗前值×100％

ACR20：关节压痛数改善程度及关节肿胀数改善程度≥20％，其余5项中至少3项改善程度≥20％

ACR50：关节压痛数改善程度及关节肿胀数改善程度≥50％，其余5项中至少3项改善程度≥50％

ACR70：关节压痛数改善程度及关节肿胀数改善程度≥70％，其余5项中至少3项改善程度≥70％

9)安全性观察：定期复查肝肾功能、血常规、尿常规等，对出现严重的病人不能耐受的副作用，

或发现合并有新的药物禁忌症，经过临床判断确定病人所出现的症状是由服用本药引起的以上情况时，

停用原药或减量，换用其他药物。

40例RA患者服用雷公藤多苷片12周，期间定期随访检测达标治疗的指标，根据临床调查量表将这些RA患者分为两组：达标治疗组和未达标治疗组。同时选取10名健康志愿者作为对照组；分离患者PBMC，分别使用Affymetrix miRNA 4.0芯片和Affymetrix EG1.0芯片进行miRNA和基因表达谱检测；通过Affymetrix GeneChip Scanner 3000 7G扫描仪对杂交完毕的芯片进行扫描，获得芯片原始数据。对原始数据进行标准化，内参为RNU6基因，筛选差异表达的miRNA，筛选标准为P值及倍数改变(Fold Change)。采用生物信息学差异表达数据分析方法，筛选雷公藤多苷片治疗RA疗效个体差异相关的miRNA和基因列表；基于上述差异miRNA的共表达相关性，建立miRNA共表达网络，筛选与其他差异miRNA具有显著共表达相关性的关键miRNAs；基于上述关键miRNAs靶基因与上述差异基因的相互作用信息，建立差异miRNA介导的基因表达调控网络，并通过对上述网络的拓扑特征进行分析，筛选网络中的关键基因对应的上游关键miRNAs作为评价雷公藤多苷片治疗RA效果的生物标志物。详细地说，从现有的生物分子相互作用数据库，如HAPPI、HPRD、Reactome等，提取雷公藤多苷片治疗RA达标组和未达标组差异表达基因与差异miRNA靶基因的相互作用信息，并进行去冗余处理。采用Cytoscape和NaviGator软件进行基因相互作用网络的构建。对上述相互作用网络中节点的拓扑特征进行分析，计算其相应的网络连接度(degree)、紧密度(closeness)、介度(betweenness)的值，从中挑取网络中发挥重要功能并且具有紧密相互作用关系的hub基因，选择标准是：①hub基因在网络中的连接度、紧密度、介度均大于相应计算结果的中位数(median)；②hub基因在网络中的隐藏边数(hidden-edge)占总连接边数的比例<50％；③hub基因间要有直接相互作用关系，能够形成子网。将筛选出来的hub基因的上游调控miRNAs作为评价雷公藤多苷片治疗RA疗效的生物标志物。所述生物标志物具体为hsa-miR-550b-2-5p、hsa-miR-4797-5p、hsa-miR-6509-5p和hsa-miR-378g。

基于上述生物标志物的表达量，采用偏最小二乘法算法，构建雷公藤多苷片个体化治疗RA的疗效预测模型，详细地为：将通过miRNA表达谱芯片检测技术得到的miRNA表达谱数据集，随机分为训练集和独立测试集，重复100次；用训练集训练模型中每个组分(miRNA标志物)的权重值(weight)和卡分阈值(cutoff)；用独立测试集进行性能评估；重复验证100次，计算预测准确性和ROC曲线下面积的均值及标准差；进行五倍交叉验证，评估模型的稳定性；该模型如式(1)所示。

F(c)＝sgn[f₁(c₁)+f₂(c₂)+f₃(c₃)+f₄(c₄)+b]式(1)

式(1)中，F(c)表示雷公藤多苷片治疗类风湿性关节炎患者的个体有效性，F(c)返回值为1表示有效，返回值为-1时表示无效；c₁、c₂、c₃和c₄依次分别表示hsa-miR-550b-2-5p、hsa-miR-4797-5p、hsa-miR-6509-5p和hsa-miR-378g的相对表达量；所述相对表达量是指相对于内参的表达量比值；f₁(c₁)、f₂(c₂)、f₃(c₃)和f₄(c₄)分别为依据所述建模算法训练得到的核函数，b为依据所述建模算法训练得到的临界打分值；具体地，式(1)中，f₁(c₁)＝0.5554×c₁，f₂(c₂)＝0.4323×c₂，f₃(c₃)＝-0.5148×c₃，f₄(c₄)＝-0.4895×c₄，b＝-0.0016，也就是说判别函数具体为：

F(c)＝sgn(0.5554×c₁+0.4323×c₂-0.5148×c₃-0.4895×c₄-0.0016)式(2)。

实施例2

本实施例用于说明本公开的生物标志物和预测模型的验证。

重新选取60例RA患者作为独立测试集，纳入标准、排除标准同实施例1；服用雷公藤多苷片12周，期间定期随访检测达标治疗的指标，根据临床调查量表将这些RA患者分为两组：达标治疗组和未达标治疗组；利用实时定量PCR(qRT-PCR)，检测上述独立测试样本中hsa-miR-550b-2-5p、hsa-miR-4797-5p、hsa-miR-6509-5p和hsa-miR-378g的表达量，进一步验证实施例1得到的雷公藤多苷片个体化治疗RA的疗效预测模型的性能，评价指标包括预测准确性(accuracy,ACC)接受者操作特性曲线(receiver operating characteristiccurve，简称ROC曲线)下面积(area under the curve,AUC)。ACC结果为87.09％。AUC结果为1.0。

实时定量PCR(qRT-PCR)所用的引物如表1所示：

表1

对比例1

在进行实施例1的同时，在筛选出来的差异表达miRNA中选取不同于hsa-miR-550b-2-5p、hsa-miR-4797-5p、hsa-miR-6509-5p和hsa-miR-378g的组合的miRNA组合，作为评价雷公藤多苷片治疗的生物标志物，具体为hsa-miR-3160-3p、hsa-miR-5196-5p、hsa-miR-6129、hsa-miR-6779-5p。

在进行实施例2的同时，采用与实施例2相同的方法进行评估，结果显示ACC结果为58.33％。AUC结果为0.516。

对比例2

本对比例对比了将筛选出来的四个miRNA:hsa-miR-550b-2-5p、hsa-miR-4797-5p、hsa-miR-6509-5p和hsa-miR-378g单独表达量或取其表达量平均值用于预测雷公藤多苷片疗效的预测性能，对比结果如表2和图1所示：

表2

对比例3

本对比例对比了常用临床特征值对雷公藤多苷片疗效的预测性能与PLS模型的预测性能，结果见表3所示：

表3

根据实施例2和对比例1-3的比较可以看出，本公开的技术方案为类风湿性关节炎的个体化诊疗提供了高特异性、高敏感性的生物标志物和疗效预测模型。

以上结合附图详细描述了本公开的优选实施方式，但是，本公开并不限于上述实施方式中的具体细节，在本公开的技术构思范围内，可以对本公开的技术方案进行多种简单变型，这些简单变型均属于本公开的保护范围。

另外需要说明的是，在上述具体实施方式中所描述的各个具体技术特征，在不矛盾的情况下，可以通过任何合适的方式进行组合，为了避免不必要的重复，本公开对各种可能的组合方式不再另行说明。

此外，本公开的各种不同的实施方式之间也可以进行任意组合，只要其不违背本公开的思想，其同样应当视为本公开所公开的内容。

序列表

<110> 中国中医科学院中药研究所

<120> 通过多个miRNA的表达量确定雷公藤多苷片治疗类风湿性关节炎的个体有效性的系统

<130> 8332ICMM

<160> 9

<170> SIPOSequenceListing 1.0

<210> 1

<211> 22

<212> RNA

<213> 人(Homo sapiens)

<400> 1

augugccuga gggaguaaga ca 22

<210> 2

<211> 21

<212> RNA

<213> 人(Homo sapiens)

<400> 2

gacagagugc cacuuacuga a 21

<210> 3

<211> 21

<212> RNA

<213> 人(Homo sapiens)

<400> 3

auuagguagu ggcaguggaa c 21

<210> 4

<211> 41

<212> RNA

<213> 人(Homo sapiens)

<400> 4

cacugggcuu ggagucagaa gaccuggcuc cagcccagcu c 41

<210> 5

<211> 24

<212> DNA

<213> 人工序列(Artificial Sequence)

<400> 5

cgatgtgcct gagggagtaa gaca 24

<210> 6

<211> 23

<212> DNA

<213> 人工序列(Artificial Sequence)

<400> 6

cggacagagt gccacttact gaa 23

<210> 7

<211> 24

<212> DNA

<213> 人工序列(Artificial Sequence)

<400> 7

gcgattaggt agtggcagtg gaac 24

<210> 8

<211> 22

<212> DNA

<213> 人工序列(Artificial Sequence)

<400> 8

cgactgggct tggagtcaga ag 22

<210> 9

<211> 23

<212> DNA

<213> 人工序列(Artificial Sequence)

<400> 9

cgagtgcagg gtccgaggta ttc 23