孫曉瑩,鄧雪蓉,李光韜,謝文慧,樊勇,張慧娟,楊新蕾,張卓莉
隨著達標治療(treat to target, T2T)理念的提出,類風(fēng)濕關(guān)節(jié)炎(rheumatoid arthritis, RA)患者的治療目標為疾病緩解或低疾病活動度。在臨床緩解的關(guān)節(jié)中,通過關(guān)節(jié)超聲檢查仍可以發(fā)現(xiàn)亞臨床滑膜炎癥[1]。亞臨床滑膜炎的存在可預(yù)測RA復(fù)發(fā)及關(guān)節(jié)結(jié)構(gòu)破壞進展[2]。
很少有研究關(guān)注RA亞臨床滑膜炎相關(guān)的血清學(xué)生物標志物。先前研究發(fā)現(xiàn),鈣衛(wèi)蛋白與人堿性成纖維細胞生長因子(basic fibroblast growth factor, bFGF)可作為在RA臨床緩解患者中識別超聲定義的活動性滑膜炎的生物標志物[3-4]。然而,上述標志物仍需進一步驗證,目前尚未見關(guān)于RA亞臨床滑膜炎轉(zhuǎn)錄組學(xué)變化的研究。在既往研究中,RA的轉(zhuǎn)錄組分析通常集中在外周血單個核細胞(peripheral blood mononuclear cell, PBMC)[5]。本研究旨在探索RA亞臨床滑膜炎患者PBMC轉(zhuǎn)錄組學(xué)特點,研究亞臨床滑膜炎發(fā)生機制。
本研究前瞻性入選2019年5月至2019年7月就診于北京大學(xué)第一醫(yī)院風(fēng)濕免疫科臨床緩解的RA患者共10例,其中亞臨床滑膜炎組及超聲下緩解組患者各5例。定義所有關(guān)節(jié)超聲下能量多普勒(power doppler, PD)分級均為0級的患者為超聲下緩解組,任意一個關(guān)節(jié)PD≥1級的患者為亞臨床滑膜炎組。
RA患者納入標準:(1)符合1987年美國風(fēng)濕病學(xué)會(American College of Rheumatology, ACR)或2010年ACR/歐洲抗風(fēng)濕病聯(lián)盟(European League Against Rheumatism, EULAR) RA分類標準;(2)年齡>18歲;(3)基于C反應(yīng)蛋白的28個關(guān)節(jié)疾病活動度評分(disease activity score based on 28 joint count and CRP, DAS28-CRP)<2.6。排除標準:(1)合并其他關(guān)節(jié)炎性疾病;(2)存在關(guān)節(jié)明顯畸形或殘毀,或有缺如及外傷等。
本研究方案經(jīng)本院倫理委員會審批,獲得所有入組患者的口頭或書面知情同意。
所有患者在同一天完成臨床評估、實驗室及超聲檢查,并用EDTA抗凝采血管采集所有患者外周靜脈血10 mL。記錄患者的臨床資料,包括姓名、性別、年齡、28個關(guān)節(jié)(包括雙側(cè)肩、肘、腕、膝、掌指關(guān)節(jié)[metacarpophalangeal, MCP]1-5和近端指間關(guān)節(jié)[proximal interphalangeal, PIP]1-5關(guān)節(jié))的腫脹關(guān)節(jié)計數(shù)(swollen joint counts, SJC)及壓痛關(guān)節(jié)計數(shù)(tender joint counts, TJC)、患者整體評估(patient’s global assessment, PGA)、評估者整體評估(evaluator’s global assessment,EGA)以及實驗室資料,包括紅細胞沉降率(erythrocyte sedimentation rate, ESR)及C反應(yīng)蛋白(C-reactive protein, CRP)。通過DAS28-ESR、DAS28-CRP、簡化臨床疾病活動指數(shù)(simplified disease activity index, SDAI)、臨床疾病活動指數(shù)(clinical disease activity index, CDAI)等指標評估患者的疾病活動度。應(yīng)用美國GE LOGIQ-E9型彩色超聲波診斷儀和ML探頭(6-15MHz),分別從掌側(cè)及背側(cè)對28個關(guān)節(jié)及雙踝、雙跖趾關(guān)節(jié)(metatarsophalangeal, MTP)1-5關(guān)節(jié)共40個關(guān)節(jié)進行超聲檢查評估滑膜炎?;译A(grey scale, GS)及PD超聲半定量分級標準均按照Szkudlarek等[6]提出的0~3級半定量評分系統(tǒng)。
從患者新鮮的外周靜脈血中分離PBMC,用Trizol法進行細胞總RNA提取。應(yīng)用Agilent 2100生物分析儀檢測總RNA的濃度、RIN/RQN值、28S/18S及片斷大小。
使用BGISEQ平臺測序,測序所得數(shù)據(jù)為原始測序數(shù)據(jù),隨后對其進行質(zhì)控,以確定測序數(shù)據(jù)是否適合后續(xù)分析。質(zhì)控后,經(jīng)過濾得到過濾后的數(shù)據(jù)比對到參數(shù)序列。比對完,通過統(tǒng)計比對率、讀長在參考序列上的分布情況,判斷比對結(jié)果是否通過第二次質(zhì)控。若通過,進行基因定量分析、基于基因表達水平的各項分析,篩選出兩組間差異表達基因(differentially expressed genes,DEGs),進行GO功能顯著性富集分析、KEGG通路顯著性富集分析、聚類、蛋白網(wǎng)絡(luò)互作(protein-protein interaction,PPI)分析及轉(zhuǎn)錄因子深入挖掘分析。本研究使用DEseq2算法進行差異基因檢測,采集差異倍數(shù)在兩倍以上且校正P值≤0.05篩選差異基因。應(yīng)用DIAMOND將DEGs對比至STRING數(shù)據(jù)庫,利用與已知蛋白的同源性獲得DEG編碼蛋白的互作關(guān)系。
計量資料先進行正態(tài)分布檢驗,符合正態(tài)分布計量資料用均值±標準差表示,非正態(tài)分布的計量資料用中位數(shù)(四分位間距)表示;計數(shù)資料以率表示。組間計量資料均數(shù)比較采用兩個獨立樣本的t檢驗(正態(tài)分布)。非正態(tài)分布計量資料的比較使用Wilcoxon符號秩和檢驗,計數(shù)資料兩組間比較用卡方檢驗或Fisher精確檢驗。用 SPSS 21.0 軟件進行統(tǒng)計分析。以P<0.05為差異有統(tǒng)計學(xué)意義。
臨床緩解RA超聲下緩解組與亞臨床滑膜炎組患者人口學(xué)、實驗室指標及疾病活動度間無顯著差異(表1)。40個關(guān)節(jié)GS總評分在兩組患者之間也無顯著差異。亞臨床滑膜炎組患者40個關(guān)節(jié)總PD評分為2(1,4.5),超聲下緩解組患者40個關(guān)節(jié)總PD評分為0(0,0),兩組間差異有統(tǒng)計學(xué)意義(P<0.05)(表1)。
表1 臨床緩解RA患者亞臨床滑膜炎組及超聲下緩解組人口學(xué)及臨床特征Table 1 Demographic and clinical characteristics of RA patients with ultrasound-defined subclinical synovitis and patientswith ultrasound determined remission
本項目使用BGISEQ平臺一共測了10個樣本,亞臨床滑膜炎組及超聲下緩解組各5個,樣本比對基因組的平均比對率為 96.43%,比對基因集的平均比對率為67.29%;一共檢測到18 397個基因。按照|log 2 FC|>1.0及P<0.05的標準,與超聲緩解組患者對比,亞臨床滑膜炎組患者DEGs共 304個,其中106個基因上調(diào),198 個基因下調(diào)。
從篩選出的DEGs 中,利用 STRING 平臺構(gòu)建PPI網(wǎng)絡(luò)圖(圖1)。其中篩選連接節(jié)點數(shù)排序前20個樞紐基因(表2), 包括2′-5′-寡腺苷酸合成酶(2′-5′-oligoadenylate synthetase, OAS)1、OAS2、OAS3、OASL、干擾素誘導(dǎo)蛋白(interferon-induced protein, IFI)44、IFI44L、黏液病毒/流感病毒抗性1(myxovirus resistance 1, MX1)、MMP9、干擾素誘導(dǎo)蛋白四肽重復(fù)序列(interferon-induced protein with tetratricopeptide repeats, IFIT)1、IFIT2、IFIT3、ISG15泛素樣修飾劑(ISG15 ubiquitin like modifier, ISG15)、XAF1、干擾素調(diào)節(jié)因子7(interferon regulatory factor 7, IRF7)、含2個自由基的S-腺苷甲硫氨酸域(radical S-adenosyl methionine domain containing 2, RSAD2)、泛素特異性肽酶18(Ubl carboxyl-terminal hydrolase 18, USP18)、受體轉(zhuǎn)運蛋白4(receptor-transporting protein 4, RTP4)、鳥苷酸結(jié)合蛋白1(guanylate-binding protein 1, GBP1)、CXCL10、DDX60。在RA亞臨床滑膜炎組患者的PBMC中,除MMP9為上調(diào)基因外,其余均為下調(diào)基因。
圖 1 PPI 網(wǎng)絡(luò)分析圖
表2 PPI網(wǎng)絡(luò)圖中連接點數(shù)排序前20個差異表達基因Table 2 Top 20 differentially expressed genes (DEGs)in PPI network
基于生物過程(biological process, BP)、細胞組成(cellular component, CC)、分子功能(molecular function, MF)的GO富集分析如圖2所示,與超聲下緩解組相比,亞臨床滑膜炎組PBMC顯示出多種富集,如BP類別中的免疫功能、Ⅰ型干擾素(interferon,IFN)組信號通路,CC類別中的細胞外間隙,MF類別中的絲氨酸型內(nèi)肽酶活性及OAS活性。KEGG富集結(jié)果顯示這些靶基因在抗原處理和遞呈、IL-17信號通路及自然殺傷細胞介導(dǎo)的細胞毒性等通路高度富集(圖3),且這些富集通路中的DEGs 主要影響運輸和分解代謝過程以及細胞生長和死亡過程。
圖 2 GO生物過程、細胞組成及分子功能分析
圖 3 KEGG通路富集分析差異表達基因
在RA患者的日常臨床實踐中,主要通過患者的體格檢查及急性期反應(yīng)物等臨床指標評估疾病活動性[7-10]。然而,部分處于臨床緩解狀態(tài)的患者仍會出現(xiàn)疾病復(fù)發(fā)及放射學(xué)進展[1,11]。近年來,隨著超聲技術(shù)的發(fā)展,處于臨床緩解的RA患者使用敏感的超聲檢查可發(fā)現(xiàn)臨床檢查及常規(guī)實驗室檢查不容易檢測到的殘留炎癥,即亞臨床滑膜炎[12]。多項研究證實,超聲下亞臨床滑膜炎與疾病復(fù)發(fā)及影像學(xué)進展密切相關(guān)[2,12]。與大多數(shù)高脂血癥或糖尿病等慢性疾病不同,RA沒有金標準的血清生物標志物。因此,鑒定RA亞臨床滑膜炎的敏感血清生物標志物具有很大的挑戰(zhàn)性。
本研究是第一個關(guān)于RA亞臨床滑膜炎患者PBMC轉(zhuǎn)錄組學(xué)特點的研究。在達臨床緩解的RA患者中初步篩選出可能與亞臨床滑膜炎發(fā)生相關(guān)的20個DEGs,其中MMP9基因在RA亞臨床滑膜炎組患者的PBMC中上調(diào),而下調(diào)的DEGs主要為Ⅰ型IFN信號通路相關(guān)的基因。篩選出的上述DEGs主要與RA發(fā)病機制和藥物代謝有關(guān),可能是闡明RA亞臨床滑膜炎分子發(fā)病機制的關(guān)鍵基因標志。
MMP9屬于MMPs家族的明膠醇類,在滑膜被覆細胞以及滑膜下組織中浸潤的炎性細胞中表達,可以溶解明膠、Ⅰ型、Ⅴ型和Ⅹ型膠原以及粘合素,而這些成分也是關(guān)節(jié)軟骨的重要組成成分。此外,MMP9能夠激活和調(diào)節(jié)與血管新生相關(guān)的生長因子。Stojanovic等[13]發(fā)現(xiàn)MMP9水平在RA患者血漿及滑液中均高于骨關(guān)節(jié)炎患者。Hakim等[14]研究顯示MMP9基因的表達水平隨DMARDs藥物治療而下降。本研究顯示相對超聲下緩解組患者,MMP9基因在RA亞臨床滑膜炎組患者的PBMC中上調(diào),提示MMP9表達上調(diào)可能參與亞臨床滑膜炎的發(fā)生,是一個潛在的鑒別亞臨床滑膜炎的血清學(xué)標志物。MMP9參與亞臨床滑膜炎的發(fā)生可能與其激活和調(diào)節(jié)血管新生相關(guān)的生長因子等機制有關(guān),目前尚不清楚,需要在未來的研究中進一步闡明。
本研究發(fā)現(xiàn)在亞臨床滑膜炎組患者中,下調(diào)的DEGs大部分為Ⅰ型IFN相關(guān)的基因,主要參與Ⅰ型IFN介導(dǎo)的信號通路。IFN是細胞因子超家族,分為Ⅰ型和Ⅱ型,Ⅰ型包括IFNα和IFNβ。IFN均通過Janus激酶(JAK)信號轉(zhuǎn)導(dǎo)子和STAT通路觸發(fā)細胞內(nèi)信號轉(zhuǎn)導(dǎo)級聯(lián)反應(yīng),從而在生物學(xué)和病理學(xué)條件下發(fā)揮多種免疫功能。Ⅰ型IFN已被證明起著將先天免疫和適應(yīng)性免疫聯(lián)系在一起的信號作用,其對免疫系統(tǒng)具有雙重作用,表現(xiàn)為免疫刺激和免疫抑制功能[15]。Ⅰ型IFN在RA發(fā)病機制中的作用目前尚不清楚。van der Pouw Kraan等[16]學(xué)者曾發(fā)現(xiàn)RA患者外周血細胞的基因表達分析顯示了大量與Ⅰ型IFN活性相關(guān)的基因過表達,Ⅰ型IFN介導(dǎo)的免疫也明顯上調(diào)。而He等[17]提出Ⅰ型IFN相關(guān)基因可能是RA的特征基因,可能通過影響免疫反應(yīng)而與RA發(fā)生相關(guān)。Ⅰ型IFN相關(guān)基因還與雷公藤、TNF抑制劑、IL-6抑制劑及CD20單抗等藥物治療反應(yīng)相關(guān),其表達水平的變化對不同藥物治療有所區(qū)別[5,18-20]。另外,Ⅰ型IFN與TNF-α之間存在交叉調(diào)節(jié)。IFNβ在體外以及關(guān)節(jié)炎動物模型體內(nèi)均可下調(diào)TNF-α生成[21-22]。相反,TNF-α在體外抑制漿細胞樣樹突狀細胞產(chǎn)生的IFN-β,在系統(tǒng)性幼年關(guān)節(jié)炎患者及干燥綜合征患者體內(nèi)TNF-α的阻斷可導(dǎo)致Ⅰ型IFN通路激活增加[23-24]??傊?,先前研究的結(jié)果表明Ⅰ型IFN相關(guān)基因可能是RA特征基因,可能通過影響免疫反應(yīng)而與RA的發(fā)展相關(guān)。根據(jù)其免疫雙重作用,它們在RA的作用可能有弊或有利。既往研究顯示在RA患者的滑膜組織中可以檢測到IFNβ蛋白,其中FLS是導(dǎo)致RA滑膜中IFNβ水平升高的原因[25-26]。另外,有學(xué)者提出RA滑膜中炎性細胞可能是從預(yù)先激活的外周細胞中募集的,而不是或不僅是通過局部產(chǎn)生[27]。因此,外周血PBMC中Ⅰ型IFN表達水平可反應(yīng)局部滑膜組織中的表達水平。本研究Ⅰ型IFN相關(guān)的基因在超聲下緩解組PBMC中相對上調(diào)可能暗示Ⅰ型IFN在RA中主要起到抗炎作用,亦或者由于Ⅰ型IFN高水平患者滑膜增殖主要以炎癥活動為主,使其對藥物治療反應(yīng)好,滑膜炎可很快消退,這值得進一步探討。
此外,亞臨床滑膜炎組患者下調(diào)的DEGs還包括CXCL10。CXCL10由IFN-γ刺激后分泌產(chǎn)生,具有強大的招募中性粒細胞、促進細胞因子分泌的功能,可介導(dǎo)炎癥反應(yīng),通常在RA患者外周血及滑膜組織表達增高[28]。然而,CXCL10在RA整個病程中是否持續(xù)發(fā)揮作用及其與RA疾病活動度及骨破壞等之間的關(guān)系目前尚未闡明。有研究發(fā)現(xiàn),CXCL10在部分RA患者的滑液和滑膜組織中水平很高,但血清中幾乎檢測不到CXCL10[29]。Kraan等[30]研究顯示,RA患者滑膜組織中CXCL10在臨床受累及臨床未受累關(guān)節(jié)中的表達無明顯差異,表明CXCL10水平與RA臨床活動性可能無明顯相關(guān)。Strieter等[31]體外和體內(nèi)研究均顯示,CXCL10是一種血管生成抑制劑, 其作用機制可能與清除一些其他的生長因子, 如IL-8和bFGF等有關(guān)。提示在本研究中CXCL10基因水平下調(diào)與亞臨床滑膜炎發(fā)生有關(guān)可能是由于其對血管形成相關(guān)因子的抑制作用下降有關(guān),值得進一步深入研究。
本研究存在一些局限性:第一,本研究納入的樣本量小,兩組人群性別、年齡等基本信息存在一定差異,可能導(dǎo)致偏倚。第二,兩組人群均為RA達臨床緩解的患者,目前疾病均處于相對穩(wěn)定狀態(tài)。部分超聲下緩解組患者仍有輕度滑膜增生,使得兩組患者之間差異較小。第三,本研究采用全外周血單個核細胞,未區(qū)分細胞群。白細胞亞群的動態(tài)變化也可能會混淆結(jié)果[32]。第四,本研究為初步探索性研究,僅研究了亞臨床滑膜炎相關(guān)的轉(zhuǎn)錄組學(xué),而mRNA水平未必能代表蛋白水平,本研究結(jié)果有待在基因及蛋白水平上進一步驗證。另外,由于RA亞臨床滑膜炎患者的滑膜組織及滑液不容易獲取,本研究僅對患者PBMC進行轉(zhuǎn)錄組學(xué)分析。未來將收集患者滑膜組織及滑液進行轉(zhuǎn)錄組學(xué)檢測。第五,本研究僅納入臨床及超聲均緩解的RA患者作為對照,缺乏健康人群及活動期RA患者作為對照,未來研究將納入這兩組人群作為對照組進行驗證,探索特異于RA亞臨床滑膜炎患者的血清學(xué)標志物。
綜上所述,在RA達臨床緩解的患者,外周血MMP9升高以及Ⅰ型干擾素相關(guān)通路、CXCL10下調(diào)可能提示亞臨床滑膜炎存在的潛在生物學(xué)標志物,或許可以為揭示亞臨床滑膜炎的發(fā)生提供新的視角和潛在的新靶點。但其分子發(fā)病機制目前尚不清楚,仍需在基因及蛋白水平進一步驗證,并深入探索。