DAS28CRP到DAS28ESR换算方程的建立和阈值的推导
2022-03-28普利明马剑达林建子戴冽凌莉
普利明?马剑达?林建子?戴冽?凌莉
【摘要】目的 以基于ESR的28个关节疾病活动度评分(DAS28ESR)为标准,建立基于CRP的28个关节疾病活动度评分(DAS28CRP)到DAS28ESR的换算方程,并推导DAS28CRP的阈值。方法 选择于中山大学孙逸仙纪念医院风湿免疫科就诊且同时测量DAS28ESR和DAS28CRP的患者。患者被随机分至训练集(70%)和验证集(30%)中。采用加权Robust回归建立DAS28CRP到DAS28ESR的换算方程,利用Wilcoxon 符号秩检验比较转换值和实际值差异,采用kappa值评估疾病活动度分类一致性。结果 共纳入955例患者。换算方程为DAS28ESR1 = 1.06×DAS28CRP + 0.4275。根据换算方程,DAS28CRP相应的阈值为2.05、2.61和4.40。经过换算后,DAS28ESR转换值与实际值比较差异无统计学意义(P > 0.05),且kappa值由0.54提升至0.78。结论 成功建立从DAS28CRP到DAS28ESR的换算方程,这有助于不同研究间的比较,所推导出的DAS28CRP阈值有助于提高其与DAS28ESR在疾病活动度评估方面的一致性。
【关键词】类风湿关节炎;疾病活动度;换算方程;阈值;红细胞沉降率;C-反应蛋白
Establishment of conversion equation from DAS28CRP to DAS28ESR and derivation of threshold of DAS28CRP Pu Liming△, Ma Jianda, Lin Jianzi, Dai Lie, Ling Li. △Department of Medical Statistics, School of Public Health, Sun Yat-sen University, Guangzhou 510080, China
Corresponding author, Dai Lie, E-mail: dailie@mail.sysu.edu.cn; Ling Li, E-mail: lingli@mail.sysu.edu.cn
【Abstract】Objective To establish the conversion equation from DAS28CRP to DAS28ESR and derive the corresponding thresholds of DAS28CRP using DAS28ESR as the standard. Methods Patients who received simultaneous measurement of DAS28ESR and DAS28CRP at Department of Rheumatology of Sun Yat-sen Memorial Hospital were recruited and randomly assigned into the training set (70%) and validation set (30%). The Bisquare-weighted robust regression was applied to establish the conversion equation from DAS28CRP to DAS28ESR. The Wilcoxon signed-rank test was used to compare the difference between the converted values and actual values of DAS28ESR and the kappa values were used to assess the agreement in disease activity classification. Results A total of 955 patients were enrolled in the study. The conversion equation was: DAS28ESR1 = 1.06 × DAS28CRP + 0.4275. According to the conversion equation, the corresponding thresholds for DAS28CRP were 2.05, 2.61 and 4.40, respectively. After conversion, there was no statistical difference between the converted values and actual values of DAS28ESR (P > 0.05). The kappa values were increased from 0.54 to 0.78. Conclusions The conversion equation from DAS28CRP to DAS28ESR is successfully established, which facilitates the comparison of findings among different studies. The derived DAS28CRP thresholds contribute to improving the agreement with DAS28ESR in assessing the disease activity.
【Key words】Rheumatoid arthritis; Disease activity; Conversion equation; Threshold;ESR;CRP
類风湿关节炎(RA)是一种常见的全身性自身免疫疾病,控制疾病活动度,实现缓解是RA的治疗目标[1]。美国风湿病学会(ACR)、欧洲风湿病联盟(EULAR)以及我国《2018中国类风湿关节炎诊疗指南》均建议定期评估RA患者的疾病活动度,根据评估结果调整治疗方案[2-4]。28个关节疾病活动度评分(DAS28)是临床广泛应用的疾病活动度评估指标。1995年Prevoo等首先提出使用基于ESR的DAS28(DAS28ESR)评估疾病活动度,并根据特定阈值将患者分为缓解、低疾病活动度、中疾病活动度和高疾病活动度。但是考虑到ESR容易受年龄、性别等影响,Fransen等(2004年)和Tamhane等(2013年)建议使用CRP代替ESR进行评估,即DAS28CRP,并建议使用DAS28ESR阈值,但是未进行验证。此外,ACR和EULAR未明确区分DAS28ESR和DAS28CRP的阈值。因此,在临床实践及临床研究中,医师可能将两者互换使用[5]。也有研究使用DAS28ESR和DAS28CRP对同一患者群体进行评估,结果表明与DAS28ESR相比,DAS28CRP倾向于低估疾病活动度,建议两者不能互换使用[6-11]。这给不同研究间的比较带来挑战。考虑到DAS28ESR和DAS28CRP间有强相关性,因此可通过建立换算方程实现两者相互转换[8, 12]。本研究旨在以DAS28ESR为标准,建立从DAS28CRP到DAS28ESR的换算方程,同时基于方程以DAS28ESR为标准推导DAS28CRP相应阈值。
对象与方法
一、研究对象
选取从2010年7月至2019年3月在中山大学孙逸仙纪念医院风湿免疫科首次就诊的RA患者。纳入标准为:①符合2009年ACR/EULAR诊断标准;②同时测量DAS28ESR和DAS28CRP的RA患者[13]。本研究经中山大学孙逸仙纪念医院伦理委员会批准(批件号:SYSEC-KY-KS-011),由于为回顾性研究,且研究数据收集过程无创、不涉及敏感信息,经伦理委员会审查批准豁免知情同意。
二、临床资料的收集
收集RA患者的人口学、临床和实验室变量:①人口学变量包括性别和年龄;②临床变量包括病程、28个压痛关节计数(TJC28)、28个肿胀关节计数(SJC28)及患者病情活动度评分(PtGA);③实验室变量包括ESR、CRP、RF和抗环瓜氨酸肽抗体(ACPA)。
三、疾病活动度评估
根据指南推荐的DAS28ESR和DAS28CRP计算公式如下[14]:
DAS28ESR = 0.56 × + 0.28 ×
+ 0.70 ×ln(ESR)+ 0.014×PtGA
DAS28CRP = 0.56 × + 0.28 ×
+ 0.36×ln(CRP+1)
+ 0.014 × PtGA + 0.96
本研究中DAS28ESR1被定义为根据换算方程从DAS28CRP计算得到的DAS28ESR转换值。根据特定阈值将患者分为缓解(DAS28 < 2.6)、低疾病活动度(2.6 ≤ DAS28 < 3.2)、中疾病活动度(3.2 ≤ DAS28 < 5.1)和高疾病活动度(DAS28 ≥ 5.1)。
四、统计学处理
经正态性检验,本研究中定量资料均不符合正态分布,因此采用M(P25,P75)进行统计学描述,定性资料采用频数(频率)进行统计学描述。将患者随机分为训练集(70%)和验证集(30%),采用Wilcoxon秩和检验和Pearson χ2检验比较2个数据集变量差异是否存在統计学意义。采用加权Robust回归在训练集中建立换算方程。
换算方程的验证包括定量和定性验证,首先根据换算方程将DAS28CRP值转换为DAS28ESR转换值,即DAS28ESR1。采用Wilcoxon符号秩检验比较DAS28ESR实际值和DAS28CRP实际值、DAS28ESR1转换值差异是否存在统计学意义,采用kappa值及其95%置信区间(95%CI)评估DAS28ESR与DAS28CRP、DAS28ESR1在疾病活动度分类上的一致性。
依据换算方程,推导DAS28CRP相应缓解、低疾病活动度、中疾病活动度和高疾病活动度阈值。检验水准为双侧α= 0.05,以上统计分析在MATLAB(R2021a)和SAS 9.4中完成。
结果
一、训练集和验证集患者人口学和临床特征比较
共955例患者被纳入研究,年龄为52.00(43.00,61.00)岁,其中女性占79.39%。训练集和验证集中,一般人口学特征和临床变量比较差异均无统计学意义(P均> 0.05),见表1。
二、DAS28ESR与DAS28CRP的比较
全数据集中,DAS28ESR和DAS28CRP分别为4.54(3.11,6.08)和3.89(2.60,5.20),DAS28ESR与DAS28CRP实际值差值为0.70(0.40,0.94),Wilcoxon符号秩检验表明差值具有统计学意义(P < 0.001),DAS28CRP低于DAS28ESR。
另外,DAS28ESR缓解率为17.49% (167/955),DAS28CRP缓解率为 25.03% (239/955) 。两者不一致率为32.77% (313/955),kappa (95%CI) 值为0.54(0.50,0.58),见表2。
三、换算方程的建立及验证
DAS28ESR和DAS28CRP为线性关系(图1),建立的换算方程为
DAS28ESR1 = 1.06×DAS28CRP + 0.4275
通过换算,患者DAS28ESR1转换值与DAS28ESR实际值差值为-0.02(-0.28,-0.20),2组比较差异无统计学意义(P > 0.05),见表3。通过换算,DAS28ESR和DAS28CRP的不一致率从32.75%下降至15.33%,kappa (95%CI) 值从0.54(0.46,0.61) 上升为0.78(0.72,0.84),见表4。
四、 DAS28CRP阈值的推导
将DAS28ESR缓解、低疾病活动度、中疾病活动度的阈值2.6、3.2、5.1分别代入换算方程,推导的DAS28CRP相应阈值分别为2.05、2.61、4.40。
讨论
随着RA新治疗方法的增多,比较2种或以上治疗策略的相对疗效有助于选择最佳的治疗方案[15]。DAS28ESR和DAS28CRP是临床研究中广泛使用的结局指标,了解两者能否互换具有重要意义[16-17]。本研究显示,DAS28CRP值低于DAS28ESR值,当使用相同阈值对患者进行疾病活动度评估时,DAS28ESR与DAS28CRP不一致率高达32.77%,超过临床上接受的30%上限值[18]。因此,DAS28ESR和DAS28CRP不能直接互换使用。
其他学者也得到了类似的结果。Song等[11]分析了30 501例中国患者数据,其中80%患者为中度至高度疾病活动度,发现DAS28CRP低于DAS28ESR(4.5±1.6 vs. 5.1±1.7),袁钰荔等[9]和秦晨曼[10]的研究也得到相似的结果。Hamann等[6]、Nielung等[7]和Eissa等[8]分别在英国、丹麦和埃及人群中进行了研究,亦得到相似的结果,即 DAS28CRP值低于DAS28ESR值,DAS28CRP比DAS28ESR更容易实现缓解。Tamhane等(2013年)发现即使按照性别、年龄和病程分层,DAS28CRP值仍低于DAS28ESR值。Siemons等[12]也同样认为DAS28ESR与DAS28CRP无法互换使用。一项荟萃分析表明,DAS28CRP的缓解率高于DAS28ESR(18%~22% vs. 5%~10%)[19]。对疾病活动度分类的一致性评估中,Son等[20]分析了540例韩国RA患者数据,发现使用DAS28CRP评估的缓解率高于DAS28ESR(23.3% vs. 24.8%),两者kappa值为0.45,表现中度一致性。Greenmyer等[21]报道DAS28ESR和DAS28CRP的kappa系数仅0.307,不一致率高达33.9%。
DAS28ESR和DAS28CRP的不一致带来两方面的挑战。一方面,对于临床实践而言,医师根据患者疾病活动度评估结果调整用药方案,评估结果的不一致将导致不同的治疗决策。例如DAS28CRP低估疾病活动度可能造成患者无法得到足够的治疗。另一方面,DAS28ESR和DAS28CRP间的不一致性限制了使用不同指标的研究间的比较。
目前有2种方法可以用于提高DAS28ESR与DAS28CRP一致性,包括建立换算方程以及调整阈值。由于DAS28ESR与DAS28CRP由相似的单独指标构成,以及DAS28ESR和DAS28CRP之间的强相关性,因此本研究通过加权Robust回归建立了DAS28CRP到DAS28ESR的线性换算方程。Son等[20]和Sengul等[22]也报道了DAS28ESR和DAS28CRP之间具有线性关系。本研究将换算方程进行定量和定性两方面验证,作为定量变量时,DAS28ESR实际值与DAS28ESR1转换值的差异无统计学意义。当作为定性变量时,不一致率由32.75%降至15.33%,kappa值由0.54提升至0.78。为了便于临床实践中医师根据阈值评估疾病活动度,进而选择治疗方案,我们基于换算方程重新推导DAS28CRP阈值,结果分别为2.05、2.61、4.40。
基于其他国家患者所推导出的DAS28CRP的阈值也均低于目前指南推荐的阈值。Inoue等[23]基于日本国家大型RA观察性数据库,分析了6729例患者,以DAS28ESR为标准,建议DAS28CRP相应缓解、低疾病活动度、中疾病活动度的阈值应调整为2.3、2.7和4.1。Park等[24]基于198例韩国RA患者建议的DAS28CRP阈值为2.19、2.60、4.07。Kuriya等[25]基于995例加拿大患者数据,建议的DAS28CRP相应阈值则为2.5、2.9和4.6,高于亚洲患者。Fleischmann等[26]分析了5项随机对照试验数据,患者来自欧洲、亚洲、北美洲、拉丁美洲及澳大利亚5个国家或地区,建议DAS28CRP相应的缓解和低疾病活动度阈值分别为2.4和2.6。进一步分析表明,亚洲患者的阈值低于非亚洲患者,缓解和低疾病活动度阈值分别为2.2 vs. 2.4以及2.6 vs. 2.9。不同人群中阈值差异的原因,可能是CRP水平与基因多态性有关[22]。因此Eissa等[8]和Sengul等[22]建议在不同种族人群中建立相应阈值。本研究为单中心研究,纳入患者为中国RA患者,与其他种族相比,CRP和ESR的关系可能不同,因此无法代表更广泛的RA患者群体,方程和阈值需要进一步外部验证。但是本研究可以为其他人群中换算方程及阈值的建立提供方法学参考。
综上所述,本研究建立的换算方程,将DAS28CRP值转化为DAS28ESR值,有助于实现不同研究间的比较。同时推导的DAS28CRP阈值,有助于提高DAS28CRP与DAS28ESR评估的一致性。
参 考 文 献
[1] 陈乐锋, 马剑达, 李谦华, 等. 生物制剂依从性对类风湿关节炎治疗达标的影响.新医学, 2019, 50(5):313-318.
[2] England B R, Tiong B K, Bergman M J, et al. 2019 update of the American college of rheumatology recommended rheumatoid arthritis disease activity measures. Arthritis Care Res (Hoboken), 2019, 71(12): 1540-1555.
[3] Smolen J S, Landewé R B M, Bijlsma J W J, et al. EULAR recommendations for the management of rheumatoid arthritis with synthetic and biological disease-modifying antirheumatic drugs: 2019 update. Ann Rheum Dis, 2020, 79(6): 685-699.
[4] 中華医学会风湿病学分会. 2018中国类风湿关节炎诊疗指南. 中华内科杂志, 2018, 57(4): 242-251.
[5] Fleischmann R M, van der Heijde D, Gardiner P V, et al. DAS28-CRP and DAS28-ESR cut-offs for high disease activity in rheumatoid arthritis are not interchangeable. RMD Open, 2017, 3(1): e000382.
[6] Hamann P D H, Shaddick G, Hyrich K, et al. Gender stratified adjustment of the DAS28-CRP improves inter-score agreement with the DAS28-ESR in rheumatoid arthritis. Rheumatology (Oxford), 2019, 58(5): 831-835.
[7] Nielung L, Christensen R, Danneskiold-Sams?e B, et al. Validity and agreement between the 28-joint disease activity score based on C-reactive protein and erythrocyte sedimentation rate in patients with rheumatoid arthritis. Arthritis, 2015, 2015: 401690.
[8] Eissa M, El Shafey A, Hammad M. Comparison between different disease activity scores in rheumatoid arthritis: an Egyptian multicenter study. Clin Rheumatol, 2017, 36(10): 2217-2224.
[9] 袁鈺荔, 尹玉月, 范惠珍. DAS28-CRP与DAS28-ESR评分在类风湿关节炎中的比较. 实用中西医结合临床, 2018, 18(11): 5-6, 9.
[10] 秦晨曼. 基于C-反应蛋白和血沉的类风湿关节炎患者疾病活动度DAS28的比较. 郑州: 郑州大学, 2013.
[11] Song X, Wang Y H, Li M T, et al. Chinese registry of rheumatoid arthritis: IV. Correlation and consistency of rheumatoid arthritis disease activity indices in China. Chin Med J (Engl), 2021, 134(12): 1465-1470.
[12] Siemons L, Vonkeman H E, ten Klooster P M, et al. Interchangeability of 28-joint disease activity scores using the erythrocyte sedimentation rate or the C-reactive protein as inflammatory marker. Clin Rheumatol, 2014, 33(6): 783-789.
[13] Aletaha D, Neogi T, Silman A J, et al. 2010 Rheumatoid arthritis classification criteria: an American College of Rheumatology/European League Against Rheumatism collaborative initiative. Arthritis Rheum, 2010, 62(9): 2569-2581.
[14] Anderson J, Caplan L, Yazdany J, et al. Rheumatoid arthritis disease activity measures: American College of Rheumatology recommendations for use in clinical practice. Arthritis Care Res (Hoboken), 2012, 64(5): 640-647.
[15] Hazlewood G S, Barnabe C, Tomlinson G, et al. Methotrexate monotherapy and methotrexate combination therapy with traditional and biologic disease modifying antirheumatic drugs for rheumatoid arthritis: Abridged Cochrane systematic review and network meta-analysis. BMJ, 2016, 353: i1777.
[16] Shivacheva T K. DAS28-ESR and DAS28-CRP: are they interchangeable in measuring the activity of rheumatoid arthritis in response to treatment with biological agents? Folia Med (Plovdiv), 2020, 62(1): 46-51.
[17] Studenic P, Felson D, de Wit M, et al. Testing different thresholds for patient global assessment in defining remission for rheumatoid arthritis: are the current ACR/EULAR Boolean criteria optimal? Ann Rheum Dis, 2020, 79(4): 445-452.
[18] Critchley L A, Critchley J A. A meta-analysis of studies using bias and precision statistics to compare cardiac output measurement techniques. J Clin Monit Comput, 1999, 15(2): 85-91.
[19] Smolen J S, Aletaha D, Gruben D, et al. Brief report: remission rates with tofacitinib treatment in rheumatoid arthritis: a comparison of various remission criteria. Arthritis Rheumatol, 2017, 69(4): 728-734.
[20] Son K M, Kim S Y, Lee S H, et al. Comparison of the disease activity score using the erythrocyte sedimentation rate and C-reactive protein levels in Koreans with rheumatoid arthritis. Int J Rheum Dis, 2016, 19(12): 1278-1283.
[21] Greenmyer J R, Stacy J M, Sahmoun A E, et al. DAS28-CRP cutoffs for high disease activity and remission are lower than DAS28-ESR in rheumatoid arthritis. ACR Open Rheumatol, 2020, 2(9): 507-511.
[22] Sengul I, Akcay-Yalbuzdag S, Ince B, et al. Comparison of the DAS28-CRP and DAS28-ESR in patients with rheumatoid arthritis. Int J Rheum Dis, 2015, 18(6): 640-645.
[23] Inoue E, Yamanaka H, Hara M, et al. Comparison of Disease Activity Score (DAS)28- erythrocyte sedimentation rate and DAS28- C-reactive protein threshold values. Ann Rheum Dis, 2007, 66(3): 407-409.
[24] Park S Y, Lee H, Cho S K, et al. Evaluation of disease activity indices in Korean patients with rheumatoid arthritis. Rheumatol Int, 2012, 32(2): 545-549.
[25] Kuriya B, Schieir O, Lin D, et al. Thresholds for the 28-joint disease activity score (DAS28) using C-reactive protein are lower compared to DAS28 using erythrocyte sedimentation rate in early rheumatoid arthritis. Clin Exp Rheumatol, 2017, 35(5): 799-803.
[26] Fleischmann R, van der Heijde D, Koenig A S, et al. How much does Disease Activity Score in 28 joints ESR and CRP calculations underestimate disease activity compared with the simplified disease activity index? Ann Rheum Dis, 2015, 74(6): 1132-1137.
(收稿日期:2021-12-27)
(本文編辑:林燕薇)