[go: up one dir, main page]

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

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

Info

Publication number
CN108034708B
CN108034708B CN201711342852.9A CN201711342852A CN108034708B CN 108034708 B CN108034708 B CN 108034708B CN 201711342852 A CN201711342852 A CN 201711342852A CN 108034708 B CN108034708 B CN 108034708B
Authority
CN
China
Prior art keywords
rheumatoid arthritis
mrnas
expression
detecting
grapl
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
CN201711342852.9A
Other languages
English (en)
Other versions
CN108034708A (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.)
Institute of Materia Medica of CACMS
Original Assignee
Institute of Materia Medica of CACMS
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 Institute of Materia Medica of CACMS filed Critical Institute of Materia Medica of CACMS
Priority to CN201711342852.9A priority Critical patent/CN108034708B/zh
Publication of CN108034708A publication Critical patent/CN108034708A/zh
Application granted granted Critical
Publication of CN108034708B publication Critical patent/CN108034708B/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
    • 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/68—Measuring or testing processes involving enzymes, nucleic acids or microorganisms; Compositions therefor; Processes of preparing such compositions involving nucleic acids
    • C12Q1/6876—Nucleic acid products used in the analysis of nucleic acids, e.g. primers or probes
    • C12Q1/6883—Nucleic acid products used in the analysis of nucleic acids, e.g. primers or probes for diseases caused by alterations of genetic material
    • 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
    • C12Q2600/00—Oligonucleotides characterized by their use
    • C12Q2600/106—Pharmacogenomics, i.e. genetic variability in individual responses to drugs and drug metabolism
    • 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
    • C12Q2600/00—Oligonucleotides characterized by their use
    • C12Q2600/158—Expression markers
    • 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
    • C12Q2600/00—Oligonucleotides characterized by their use
    • C12Q2600/16—Primer sets for multiplex assays
    • 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
    • C12Q2600/00—Oligonucleotides characterized by their use
    • C12Q2600/166—Oligonucleotides used as internal standards, controls or normalisation probes

Landscapes

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

Abstract

本公开提供了一种通过多个mRNA的表达量确定雷公藤多苷片治疗类风湿性关节炎的个体有效性的系统，所述多个mRNA包括MX1、OASL、SPINK1、CRK、GRAPL和RNF2。本公开还提供了一种确定雷公藤多苷片治疗类风湿性关节炎的个体有效性的分子标志物、该分子标志物在制备确定雷公藤多苷片治疗类风湿性关节炎的个体有效性的试剂盒中的用途、检测该分子标志物的表达量的试剂在制备确定雷公藤多苷片治疗类风湿性关节炎的个体有效性的试剂盒中的用途。本公开的技术方案为类风湿性关节炎的个体化诊疗提供了高特异性、高敏感性的生物标志物和疗效预测模型，并为临床制定类风湿性关节炎个体化治疗方案提供了新型、高效且无创的辅助工具。

Description

通过多个mRNA的表达量确定雷公藤多苷片治疗类风湿性关节 炎的个体有效性的系统
技术领域
本公开涉及生物医学领域，具体地，涉及一种通过多个mRNA的表达量确定雷公藤多苷片治疗类风湿性关节炎的个体有效性的系统、一种确定雷公藤多苷片治疗类风湿性关节炎的个体有效性的分子标志物、该分子标志物在制备确定雷公藤多苷片治疗类风湿性关节炎的个体有效性的试剂盒中的用途、检测该分子标志物的表达量的试剂在制备确定雷公藤多苷片治疗类风湿性关节炎的个体有效性的试剂盒中的用途。
背景技术
类风湿性关节炎(Rheumatoid arthritis，RA)是世界难治性自身免疫性疾病，好发于青壮年，其发病率高(我国约为0.36％)、致残率高(5-10年病程时致残率约为60％)，对人类的健康危害极大，近年已被发达国家列为未来25年中影响人类健康最大、医疗消费最高的五大疾病之一。
雷公藤是我国的传统中药，最早记载于《滇南本草》，其味苦性寒，具清热解毒、祛风除湿、舒筋活血、消肿止痛功效，其成药制剂雷公藤多苷片一直被我国人民用来治疗风湿病，并已引起国际医药界的广泛关注，迄今为止，尚没有一味中药或中成药能取代雷公藤多苷片在抗类风湿性关节炎临床中的重要位置。然而，由于其存在临床疗效个体差异大、有效剂量与毒性剂量相当等问题，严重影响了临床的推广和合理应用。
本发明人近年对雷公藤及其活性成分治疗类风湿性关节炎的作用进行了较系统的研究，发现雷公藤及其活性成分具有缓解类风湿性关节炎病情、抑制关节滑膜炎症、防止或减缓关节结构破坏的改善病情抗风湿药样作用。同时，从类风湿性关节炎的关键病理环节入手，阐述了雷公藤最主要的免疫活性成分雷公藤甲素抑制关节滑膜炎症和血管新生，调节软骨基质降解和破骨细胞分化系统，抑制炎症相关信号通路活化、降低炎症反应级联和循环的相关作用机制。相关研究成果有助于了解雷公藤作用特点、推动其抗类风湿性关节炎的分子药理学发展。基于此，下一步的考虑是如何使研究成果更直接指导临床用药，怎样个体化地应用雷公藤多苷片使之达到最佳治疗效果，这也是目前国际风湿病研究领域的热点和亟待解决的难题。
因此，进一步筛选与雷公藤多苷片治疗类风湿性关节炎个体差异密切相关的新的生物标记物，对于揭示雷公藤多苷片治疗类风湿性关节炎的分子机制及提高其个体化用药水平均具有重要的意义。
发明内容
本公开的目的是提供与雷公藤多苷片治疗类风湿性关节炎个体差异密切相关的新的生物标记物，提高雷公藤多苷片治疗类风湿性关节炎的个体化用药水平，从而改善雷公藤多苷片治疗类风湿性关节炎的治疗效果。
为了实现上述目的，一方面，本公开提供了一种通过多个mRNA的表达量确定雷公藤多苷片治疗类风湿性关节炎的个体有效性的系统，其中，该系统包括计算装置、用于输入类风湿性关节炎患者个体的多个mRNA的表达量的输入装置和用于输出雷公藤多苷片治疗类风湿性关节炎患者的个体有效性的输出装置；其中，所述多个mRNA包括MX1、OASL、SPINK1、CRK、GRAPL和RNF2；所述计算装置包括存储器和处理器，所述存储器中存储有计算机程序，所述处理器被配置为执行所述存储器中存储的计算机程序，以实现建模算法和如式(1)所示的判别函数的算法；所述建模算法为支持向量机算法和/或最小偏二乘算法；
F(c)＝sgn[f1(c1)+f2(c2)+f3(c3)+f4(c4)+f5(c5)+f6(c6)+b] 式(1)
式(1)中，F(c)表示雷公藤多苷片治疗类风湿性关节炎患者的个体有效性，F(c)返回值为1表示有效，返回值为-1时表示无效；c1、c2、c3、c4、c5和c6依次分别表示MX1、OASL、SPINK1、CRK、GRAPL和RNF2的相对表达量；所述相对表达量是指相对于内参的表达量比值；f1(c1)、f2(c2)、f3(c3)、f4(c4)、f5(c5)和f6(c6)分别为依据所述建模算法训练得到的核函数，b为依据所述建模算法训练得到的临界打分值。
另一方面，本公开还提供了一种确定雷公藤多苷片治疗类风湿性关节炎的个体有效性的分子标志物，其中，该分子标志物为多个mRNA，所述多个mRNA由MX1、OASL、SPINK1、CRK、GRAPL和RNF2组成。
另一方面，本公开还提供了分子标志物在制备确定雷公藤多苷片治疗类风湿性关节炎的个体有效性的试剂盒中的用途，其中，该分子标志物为多个mRNA，所述多个mRNA由MX1、OASL、SPINK1、CRK、GRAPL和RNF2组成。
再一方面，本公开还提供了检测分子标志物的表达量的试剂在制备确定雷公藤多苷片治疗类风湿性关节炎的个体有效性的试剂盒中的用途，其中，该分子标志物为多个mRNA，所述多个mRNA由MX1、OASL、SPINK1、CRK、GRAPL和RNF2组成。
类风湿性关节炎作为系统性免疫性疾病，伴随着基因、蛋白质以及分子相互作用网络的动态失调，具有总体上治愈率低、致残率高等特点。随着分子生物学领域的高速发展，科学家发现多种信号分子参与类风湿性关节炎的发生和发展，信号分子的活性受到遗传和环境中多种因素的影响，调节环路复杂的相互作用和交叉干扰常常使得人们对类风湿性关节炎相关的信号输出难以作出正确的认识和预测。鉴于类风湿性关节炎的高度复杂性，一直以来沿用的针对单一高特异性靶点的治疗策略并没有取得满意的效果。
在疾病诊疗生物标志物的研究领域，通常首先通过筛选疾病组织中的差异表达分子作为标志物，并采用人工神经网络、偏最小二乘法等构建疾病预测模型。但此种分析方法的可重复性较差。可见，进行疾病诊疗生物标志物的筛选只考虑差异表达数据是远远不够的。
本发明人从全局的角度出发，计算分析高通量生物分子表达谱数据，并整合生物分子网络的特征分析，筛选出了高可信度的疾病诊疗生物标志物及相应的计算公式，由此得到了本发明创造。具体地，首先，对类风湿性关节炎患者采用雷公藤多苷片治疗，通过DAS28、ACR20/50/70等疗效性指标划分雷公藤多苷片治疗达标/未达标组；接着，应用基因表达谱芯片技术，检测两组的基因的差异表达状况；再通过差异表达数据分析，选取差异表达基因谱，并建立基因表达调控网络；然后，筛选网络中的关键mRNAs作为评价雷公藤多苷片治疗类风湿性关节炎效果的生物标志物；最后，运用支持向量机算法和/或偏最小二乘法算法，针对上述生物标志物的表达量，构建雷公藤多苷片个体化治疗类风湿性关节炎的疗效预测模型，并对其性能进行评估。评估结果显示，本公开的技术方案为类风湿性关节炎的个体化诊疗提供了高特异性、高敏感性的生物标志物和疗效预测模型，并为临床制定类风湿性关节炎个体化治疗方案提供了一种新型、高效且无创的辅助工具，从而有效预测了类风湿性关节炎患者从临床治疗中的获益程度，从而可实现根据个体基因检测数据与药物效应的关系设计、制定个体化药物治疗方案，以充分发挥药物对机体的作用，增加首剂处方的有效性，提高治疗效果，减少药物不良反应的发生率，弥补传统或经验方案的不足，节省或合理配置医药卫生资源。
本公开的其他特征和优点将在随后的具体实施方式部分予以详细说明。
附图说明
附图是用来提供对本公开的进一步理解，并且构成说明书的一部分，与下面的具体实施方式一起用于解释本公开，但并不构成对本公开的限制。在附图中：
图1是实施例2和对比例2的ACC结果和AUC结果。
具体实施方式
以下结合附图对本公开的具体实施方式进行详细说明。应当理解的是，此处所描述的具体实施方式仅用于说明和解释本公开，并不用于限制本公开。
一方面，本公开提供了一种通过多个mRNA的表达量确定雷公藤多苷片治疗类风湿性关节炎的个体有效性的系统，其中，该系统包括计算装置、用于输入类风湿性关节炎患者个体的多个mRNA的表达量的输入装置和用于输出雷公藤多苷片治疗类风湿性关节炎患者的个体有效性的输出装置；其中，所述多个mRNA包括MX1、OASL、SPINK1、CRK、GRAPL和RNF2；所述计算装置包括存储器和处理器，所述存储器中存储有计算机程序，所述处理器被配置为执行所述存储器中存储的计算机程序，以实现建模算法和如式(1)所示的判别函数的算法；所述建模算法为支持向量机算法和/或最小偏二乘算法；
F(c)＝sgn[f1(c1)+f2(c2)+f3(c3)+f4(c4)+f5(c5)+f6(c6)+b] 式(1)
式(1)中，F(c)表示雷公藤多苷片治疗类风湿性关节炎患者的个体有效性，F(c)返回值为1表示有效，返回值为-1时表示无效；c1、c2、c3、c4、c5和c6依次分别表示MX1、OASL、SPINK1、CRK、GRAPL和RNF2的相对表达量；所述相对表达量是指相对于内参的表达量比值；f1(c1)、f2(c2)、f3(c3)、f4(c4)、f5(c5)和f6(c6)分别为依据所述建模算法训练得到的核函数，b为依据所述建模算法训练得到的临界打分值。其中，支持向量机算法和最小偏二乘算法及其运算和训练的方式为已知的常规方式。
其中，MX1、OASL、SPINK1、CRK、GRAPL和RNF2是指NCBI数据库中的基因名称(GeneSymbol)，具体注释如表1所示，表1中的参考序列仅用于举例说明而不用于限制。
表1
Figure BDA0001508790540000061
可选地，该系统还包括多个mRNA的表达量的检测装置。
可选地，所述多个mRNA的表达量的检测装置包括mRNA表达量检测芯片和芯片信号读取器，mRNA表达量检测芯片包括分别检测MX1、OASL、SPINK1、CRK、GRAPL和RNF2的表达量的探针。或者所述多个mRNA的表达量的检测装置包括实时定量PCR仪和mRNA的实时定量PCR引物，mRNA的实时定量PCR引物包括分别检测MX1、OASL、SPINK1、CRK、GRAPL和RNF2的表达量的实时定量PCR引物。
可选地，所述mRNA表达量检测芯片还包括内参探针，所述内参探针为检测RPS18基因和/或GAPDH基因的表达量的探针，所述mRNA的实时定量PCR引物，还包括内参引物，所述内参引物为检测RPS18基因和/或GAPDH基因的实时定量PCR引物。
可选地，其中，作为以GAPDH基因为内参训练得到的一组参考值，式(1)中，f1(c1)＝-0.4694×c1，f2(c2)＝-0.2494×c2，f3(c3)＝-0.5592×c3，f4(c4)＝0.3429×c4，f5(c5)＝0.4054×c5，f6(c6)＝0.3504×c6，b＝-0.03。也就是说，判别函数可以简化为如式(2)所示：
F(c)＝sgn(-0.4694×c1-0.2494×c2-0.5592×c3+0.3429×c4+0.4054×c5+0.3504×c6-0.03)
式(2)。
可选地，其中，作为以RPS18基因训练得到的一组参考值，式(1)中，f1(c1)＝-0.3031×c1，f2(c2)＝-0.3644×c2，f3(c3)＝-0.6511×c3，f4(c4)＝0.3852×c4，f5(c5)＝0.2624×c5，f6(c6)＝0.3663×c6，b＝-0.0036。也就是说，判别函数可以简化为如式(3)所示：
F(c)＝sgn(-0.3031×c1-0.3644×c2-0.6511×c3+0.3852×c4+0.2624×c5+0.3663×c6-0.0036)
式(3)。
需要说明的是，f1(c1)、f2(c2)、f3(c3)、f4(c4)、f5(c5)、f6(c6)和b可能会随着mRNA的表达量的检测手段的偏性而发生改变，也可能会随着训练数据集的数据规模大小等因素发生改变。式(2)和式(3)所示的判别函数是本公开的发明人依据实施例1中的数据以最小偏二乘的建模算法进行训练得到的，并不限制本公开的范围。当然也可以根据其他的数据，和/或，选取支持向量机算法，选择MX1、OASL、SPINK1、CRK、GRAPL和RNF2的相对表达量，进行训练，得到式(1)范围内的判别函数。
另一方面，本公开还提供了一种确定雷公藤多苷片治疗类风湿性关节炎的个体有效性的分子标志物，其中，该分子标志物为多个mRNA，所述多个mRNA由MX1、OASL、SPINK1、CRK、GRAPL和RNF2组成。
另一方面，本公开还提供了分子标志物在制备确定雷公藤多苷片治疗类风湿性关节炎的个体有效性的试剂盒中的用途，其中，该分子标志物为多个mRNA，所述多个mRNA由MX1、OASL、SPINK1、CRK、GRAPL和RNF2组成。
再一方面，本公开还提供了检测分子标志物的表达量的试剂在制备确定雷公藤多苷片治疗类风湿性关节炎的个体有效性的试剂盒中的用途，其中，该分子标志物为多个mRNA，所述多个mRNA由MX1、OASL、SPINK1、CRK、GRAPL和RNF2组成。
其中，检测分子标志物的表达量的试剂可以为探针和/或引物。
以下，通过实施例进一步详细说明本发明。
实施例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 28joints，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 EG1.0芯片进行基因表达谱检测；通过Affymetrix GeneChip Scanner 3000 7G扫描仪对杂交完毕的芯片进行扫描，获得芯片原始数据。对原始数据进行标准化，内参分别为RPS18基因和GAPDH基因，筛选差异表达的mRNA，筛选标准为P值及倍数改变(Fold Change)。采用生物信息学差异表达数据分析方法，筛选雷公藤多苷片治疗RA疗效个体差异相关的mRNA列表；基于上述差异mRNA的共表达相关性，建立mRNA共表达网络，筛选与其他差异mRNA具有显著共表达相关性的关键mRNAs；基于上述关键mRNAs靶基因与上述差异基因的相互作用信息，建立差异mRNA介导的基因表达调控网络，并通过对上述网络的拓扑特征进行分析，筛选网络中的关键基因对应的上游关键mRNAs作为评价雷公藤多苷片治疗RA效果的生物标志物。详细地说，从现有的生物分子相互作用数据库，如HAPPI、HPRD、Reactome等，提取雷公藤多苷片治疗RA达标组和未达标组差异表达基因与差异mRNA靶基因的相互作用信息，并进行去冗余处理。采用Cytoscape和NaviGator软件进行基因相互作用网络的构建。对上述相互作用网络中节点的拓扑特征进行分析，计算其相应的网络连接度(degree)、紧密度(closeness)、介度(betweenness)的值，从中挑取网络中发挥重要功能并且具有紧密相互作用关系的hub基因，选择标准是：①hub基因在网络中的连接度、紧密度、介度均大于相应计算结果的中位数(median)；②hub基因在网络中的隐藏边数(hidden-edge)占总连接边数的比例<50％；③hub基因间要有直接相互作用关系，能够形成子网。将筛选出来的hub基因的上游调控mRNAs作为评价雷公藤多苷片治疗RA疗效的生物标志物。所述生物标志物具体为MX1、OASL、SPINK1、CRK、GRAPL和RNF2。
基于上述生物标志物的表达量，采用偏最小二乘法算法，构建雷公藤多苷片个体化治疗RA的疗效预测模型，详细地为：将通过mRNA表达谱芯片检测技术得到的mRNA表达谱数据集，随机分为训练集和独立测试集，重复100次；用训练集训练模型中每个组分(mRNA标志物)的权重值(weight)和卡分阈值(cutoff)；用独立测试集进行性能评估；重复验证100次，计算预测准确性和ROC曲线下面积的均值及标准差；进行五倍交叉验证，评估模型的稳定性；该模型如式(1)所示。
F(c)＝sgn[f1(c1)+f2(c2)+f3(c3)+f4(c4)+f5(c5)+f6(c6)+b] 式(1)
式(1)中，F(c)表示雷公藤多苷片治疗类风湿性关节炎患者的个体有效性，F(c)返回值为1表示有效，返回值为-1时表示无效；c1、c2、c3、c4、c5和c6依次分别表示MX1、OASL、SPINK1、CRK、GRAPL和RNF2的相对表达量；所述相对表达量是指相对于内参的表达量比值；f1(c1)、f2(c2)、f3(c3)、f4(c4)、f5(c5)和f6(c6)分别为依据所述建模算法训练得到的核函数，b为依据所述建模算法训练得到的临界打分值；具体地，以GAPDH基因为内参时，式(1)中，f1(c1)＝-0.4694×c1，f2(c2)＝-0.2494×c2，f3(c3)＝-0.5592×c3，f4(c4)＝0.3429×c4，f5(c5)＝0.4054×c5，f6(c6)＝0.3504×c6，b＝-0.03，也就是说判别函数具体为：
F(c)＝sgn(-0.4694×c1-0.2494×c2-0.5592×c3+0.3429×c4+0.4054×c5+0.3504×c6-0.03)
式(2)。
具体地，以RPS18基因为内参时，式(1)中，f1(c1)＝-0.3031×c1，f2(c2)＝-0.3644×c2，f3(c3)＝-0.6511×c3，f4(c4)＝0.3852×c4，f5(c5)＝0.2624×c5，f6(c6)＝0.3663×c6，b＝-0.0036。也就是说，判别函数具体为：
F(c)＝sgn(-0.3031×c1-0.3644×c2-0.6511×c3+0.3852×c4+0.2624×c5+0.3663×c6-0.0036)
式(3)。
实施例2
本实施例用于说明本公开的生物标志物和预测模型的验证。
重新选取60例RA患者作为独立测试集，纳入标准、排除标准同实施例1；服用雷公藤多苷片12周，期间定期随访检测达标治疗的指标，根据临床调查量表将这些RA患者分为两组：达标治疗组和未达标治疗组；利用实时定量PCR(qRT-PCR)，检测上述独立测试样本中MX1、OASL、SPINK1、CRK、GRAPL和RNF2的表达量，进一步验证雷公藤多苷片个体化治疗RA的疗效预测模型的性能，评价指标包括预测准确性(accuracy,ACC)接受者操作特性曲线(receiver operating characteristic curve，简称ROC曲线)下面积(area under thecurve,AUC)。以GAPDH基因为内参时，式(2)判别式的ACC结果为80.20％，AUC结果为0.825。以RPS18基因为内参时，式(3)判别式的ACC结果为83.10％，AUC结果为0.802。
实时定量PCR(qRT-PCR)所用的引物如表2所示：
表2
Figure BDA0001508790540000141
对比例1
在进行实施例1的同时，在筛选出来的差异表达mRNA中选取不同于MX1、OASL、SPINK1、CRK、GRAPL和RNF2的组合的mRNA组合，作为评价雷公藤多苷片治疗的生物标志物，具体为ACTL6B,IGF1,RNF8,RPL23,GHR,ZNF384。
在进行实施例2的同时，采用与实施例2相同的方法进行评估，结果显示以GAPDH基因为内参时，ACC结果为44.44％，AUC结果为0.650。以RPS18基因为内参时，ACC结果为58.12％，AUC结果为0.625。
对比例2
本对比例对比了将筛选出来的六个mRNA:MX1、OASL、SPINK1、CRK、GRAPL和RNF2单独表达量或取其表达量平均值用于预测雷公藤多苷片疗效的预测性能，对比结果如表3和图1所示。图1中PLS模型是以RPS18为内参。
表3
Figure BDA0001508790540000151
对比例3
本对比例对比了常用临床特征值对雷公藤多苷片疗效的预测性能与PLS模型的预测性能，结果见表4所示：
表4
Figure BDA0001508790540000152
根据实施例2和对比例1-3的比较可以看出，本公开的技术方案为类风湿性关节炎的个体化诊疗提供了高特异性、高敏感性的生物标志物和疗效预测模型。
以上结合附图详细描述了本公开的优选实施方式，但是，本公开并不限于上述实施方式中的具体细节，在本公开的技术构思范围内，可以对本公开的技术方案进行多种简单变型，这些简单变型均属于本公开的保护范围。
另外需要说明的是，在上述具体实施方式中所描述的各个具体技术特征，在不矛盾的情况下，可以通过任何合适的方式进行组合，为了避免不必要的重复，本公开对各种可能的组合方式不再另行说明。
此外，本公开的各种不同的实施方式之间也可以进行任意组合，只要其不违背本公开的思想，其同样应当视为本公开所公开的内容。
序列表
<110> 中国中医科学院中药研究所
<120> 通过多个mRNA的表达量确定雷公藤多苷片治疗类风湿性关节炎的个体有效性的系统
<130> 8333ICMM
<160> 12
<170> SIPOSequenceListing 1.0
<210> 1
<211> 20
<212> DNA
<213> 人(Homo sapiens)
<400> 1
ctccgacacg agttccacaa 20
<210> 2
<211> 20
<212> DNA
<213> 人(Homo sapiens)
<400> 2
ggctcttcca gtgccttgat 20
<210> 3
<211> 20
<212> DNA
<213> 人(Homo sapiens)
<400> 3
ttcagcgagc tgcagagaaa 20
<210> 4
<211> 20
<212> DNA
<213> 人(Homo sapiens)
<400> 4
ccctctgctc cactgtcaag 20
<210> 5
<211> 20
<212> DNA
<213> 人(Homo sapiens)
<400> 5
tgaccctgtc tgtgggactg 20
<210> 6
<211> 21
<212> DNA
<213> 人(Homo sapiens)
<400> 6
tcagcaaggc ccagattttt g 21
<210> 7
<211> 20
<212> DNA
<213> 人(Homo sapiens)
<400> 7
aggcagggta gtggagtgat 20
<210> 8
<211> 20
<212> DNA
<213> 人(Homo sapiens)
<400> 8
cttcaggctt gtcccggatt 20
<210> 9
<211> 20
<212> DNA
<213> 人(Homo sapiens)
<400> 9
ctacagagag cgacgagctg 20
<210> 10
<211> 20
<212> DNA
<213> 人(Homo sapiens)
<400> 10
gaaggctccc agatggttcc 20
<210> 11
<211> 20
<212> DNA
<213> 人(Homo sapiens)
<400> 11
cgaacacctc aggaggcaat 20
<210> 12
<211> 20
<212> DNA
<213> 人(Homo sapiens)
<400> 12
tgccacttct aagggctgtg 20

Claims (9)

1.一种通过多个mRNA的表达量确定雷公藤多苷片治疗类风湿性关节炎的个体有效性的系统，其特征在于，该系统包括计算装置、用于输入类风湿性关节炎患者个体的多个mRNA的表达量的输入装置和用于输出雷公藤多苷片治疗类风湿性关节炎患者的个体有效性的输出装置；其中，所述多个mRNA包括MX1、OASL、SPINK1、CRK、GRAPL和RNF2；所述计算装置包括存储器和处理器，所述存储器中存储有计算机程序，所述处理器被配置为执行所述存储器中存储的计算机程序，以实现建模算法和如式(1)所示的判别函数的算法；所述建模算法为支持向量机算法和/或最小偏二乘算法；
F(c)＝sgn[f1(c1)+f2(c2)+f3(c3)+f4(c4)+f5(c5)+f6(c6)+b] 式(1)
式(1)中，F(c)表示雷公藤多苷片治疗类风湿性关节炎患者的个体有效性，F(c)返回值为1表示有效，返回值为-1时表示无效；c1、c2、c3、c4、c5和c6依次分别表示MX1、OASL、SPINK1、CRK、GRAPL和RNF2的相对表达量；所述相对表达量是指相对于内参的表达量比值；f1(c1)、f2(c2)、f3(c3)、f4(c4)、f5(c5)和f6(c6)分别为依据所述建模算法训练得到的核函数，b为依据所述建模算法训练得到的临界打分值。
2.根据权利要求1所述的系统，其中，该系统还包括多个mRNA的表达量的检测装置。
3.根据权利要求2所述的系统，其中，所述多个mRNA的表达量的检测装置包括mRNA表达量检测芯片和芯片信号读取器，mRNA表达量检测芯片包括分别检测MX1、OASL、SPINK1、CRK、GRAPL和RNF2的表达量的探针。
4.根据权利要求3所述的系统，其中，所述mRNA表达量检测芯片还包括内参探针，所述内参探针为检测RPS18基因和/或GAPDH基因的表达量的探针。
5.根据权利要求2所述的系统，其中，所述多个mRNA的表达量的检测装置包括实时定量PCR仪和mRNA的实时定量PCR引物，mRNA的实时定量PCR引物包括分别检测MX1、OASL、SPINK1、CRK、GRAPL和RNF2的表达量的实时定量PCR引物。
6.根据权利要求5所述的系统，其中，所述mRNA的实时定量PCR引物，还包括内参引物，所述内参引物为检测RPS18基因和/或GAPDH基因的实时定量PCR引物。
7.根据权利要求1所述的系统，其中，式(1)中，f1(c1)＝-0.4694×c1，f2(c2)＝-0.2494×c2，f3(c3)＝-0.5592×c3，f4(c4)＝0.3429×c4，f5(c5)＝0.4054×c5，f6(c6)＝0.3504×c6，b＝-0.03；或者，式(1)中，f1(c1)＝-0.3031×c1，f2(c2)＝-0.3644×c2，f3(c3)＝-0.6511×c3，f4(c4)＝0.3852×c4，f5(c5)＝0.2624×c5，f6(c6)＝0.3663×c6，b＝-0.0036。
8.一种确定雷公藤多苷片治疗类风湿性关节炎的个体有效性的分子标志物，其特征在于，该分子标志物为多个mRNA，所述多个mRNA由MX1、OASL、SPINK1、CRK、GRAPL和RNF2组成。
9.检测分子标志物的表达量的试剂在制备确定雷公藤多苷片治疗类风湿性关节炎的个体有效性的试剂盒中的用途，其特征在于，该分子标志物为多个mRNA，所述多个mRNA由MX1、OASL、SPINK1、CRK、GRAPL和RNF2组成。
CN201711342852.9A 2017-12-14 2017-12-14 通过多个mRNA的表达量确定雷公藤多苷片治疗类风湿性关节炎的个体有效性的系统 Active CN108034708B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201711342852.9A CN108034708B (zh) 2017-12-14 2017-12-14 通过多个mRNA的表达量确定雷公藤多苷片治疗类风湿性关节炎的个体有效性的系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201711342852.9A CN108034708B (zh) 2017-12-14 2017-12-14 通过多个mRNA的表达量确定雷公藤多苷片治疗类风湿性关节炎的个体有效性的系统

Publications (2)

Publication Number Publication Date
CN108034708A CN108034708A (zh) 2018-05-15
CN108034708B true CN108034708B (zh) 2020-05-22

Family

ID=62103018

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201711342852.9A Active CN108034708B (zh) 2017-12-14 2017-12-14 通过多个mRNA的表达量确定雷公藤多苷片治疗类风湿性关节炎的个体有效性的系统

Country Status (1)

Country Link
CN (1) CN108034708B (zh)

Families Citing this family (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111518882B (zh) * 2019-02-03 2022-06-28 陈卫衡 通过分子标志物辨识激素性股骨头坏死的中医证型的系统
CN110199921A (zh) * 2019-06-28 2019-09-06 广东食品药品职业学院 一种雷公藤多苷肝损伤动物模型的构建方法及其应用
CN110675956A (zh) * 2019-08-26 2020-01-10 南京医渡云医学技术有限公司 一种面瘫治疗方案的确定方法、装置、可读介质及电子设备
CN111383774A (zh) * 2020-03-13 2020-07-07 北京市神经外科研究所 用于筛选脑胶质瘤治疗方案的系统
CN111968704B (zh) * 2020-07-20 2023-05-16 中国中医科学院中药研究所 通过ceRNA调控组合的表达量确定雷公藤多苷片治疗类风湿关节炎的个体有效性的系统
CN111951883A (zh) * 2020-08-04 2020-11-17 广东省第二人民医院（广东省卫生应急医院） 一种特征mRNA表达谱组合及结肠癌早期预测方法

Family Cites Families (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5580562A (en) * 1990-03-14 1996-12-03 Board Of Regents The University Of Texas System Preparations and uses thereof for immunosuppression
CN102634561A (zh) * 2012-03-14 2012-08-15 厦门大学 雷公藤甲素及衍生物和类似物在制备抗肿瘤药物中的应用
WO2016050111A1 (en) * 2014-09-30 2016-04-07 Bgi Shenzhen Biomarkers for rheumatoid arthritis and usage thereof

Also Published As

Publication number Publication date
CN108034708A (zh) 2018-05-15

Similar Documents

Publication Publication Date Title
CN108034708B (zh) 通过多个mRNA的表达量确定雷公藤多苷片治疗类风湿性关节炎的个体有效性的系统
CN108085382B (zh) 通过多个miRNA的表达量确定雷公藤多苷片治疗类风湿性关节炎的个体有效性的系统
Tebani et al. Integration of molecular profiles in a longitudinal wellness profiling cohort
JP6873950B2 (ja) 免疫署名：早期診断および健康モニタリングの方法
Shi et al. Alterations in the intestinal microbiota of patients with severe and active Graves’ orbitopathy: a cross-sectional study
Reif et al. Integrated analysis of genetic and proteomic data identifies biomarkers associated with adverse events following smallpox vaccination
Heard et al. A computational method to differentiate normal individuals, osteoarthritis and rheumatoid arthritis patients using serum biomarkers
CN108064272B (zh) 用于类风湿性关节炎的生物标记物及其用途
Xiang et al. Associations of geriatric nutrition risk index and other nutritional risk-related indexes with sarcopenia presence and their value in sarcopenia diagnosis
Lou et al. Performance and operational feasibility of Epstein-Barr virus–based screening for detection of nasopharyngeal carcinoma: Direct comparison of two alternative approaches
CN115856309A (zh) 蛋白标志物在制备诊断帕金森或预测帕金森病的产品中的应用
Montagnese et al. 268th ENMC workshop-Genetic diagnosis, clinical classification, outcome measures, and biomarkers in Facioscapulohumeral Muscular Dystrophy (FSHD): Relevance for clinical trials
Lin et al. Development and verification of a combined diagnostic model for sarcopenia with random forest and artificial neural network
Ha et al. Blood transcriptomics to facilitate diagnosis and stratification in pediatric rheumatic diseases–a proof of concept study
Jayaseelan et al. A computational data mining strategy to identify the common genetic markers of temporomandibular joint disorders and osteoarthritis
Hu et al. Identification and construction of a disulfidptosis-mediated diagnostic model and associated immune microenvironment of osteoarthritis from the perspective of PPPM
Yang et al. Transcriptome analysis of cold syndrome using microarray
CN117457063A (zh) 一种同时推断抗乳腺癌和covid-19药物的网络药理学方法
WO2022210606A1 (ja) 認知症の将来の発症リスクの評価方法
US8738338B2 (en) Methods and compositions for biomarkers of fatigue, fitness and physical performance capacity
CN111968704B (zh) 通过ceRNA调控组合的表达量确定雷公藤多苷片治疗类风湿关节炎的个体有效性的系统
van Vliet et al. Circulating microRNAs and isomiRs as biomarkers for the initial insult and epileptogenesis in four experimental epilepsy models: The EPITARGET study
Zhou Evaluating the impacts of biological treatment in inflammatory arthritis patients with quasi-experimental study designs
Sadioglu Cagdas et al. Monitoring disease activity and treatment response in ankylosing spondylitis: a retrospective study of hematologic inflammatory markers
Liu et al. Artificial intelligence in rheumatoid arthritis research: A bibliometric analysis from 2004 to 2023

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