張 琳，毛 偉，劉 芳，枉 前，呂 軍△

（1. 中國(guó)人民解放軍陸軍軍醫(yī)大學(xué)第一附屬醫(yī)院，重慶 400038； 2. 重慶市南岸區(qū)人民醫(yī)院，重慶 400060）

類風(fēng)濕關(guān)節(jié)炎（RA）屬全身性難治性自身免疫性疾?。?］，最終可導(dǎo)致關(guān)節(jié)畸形和功能喪失［2］，可并發(fā)肺部疾?。?］、心血管疾?。?］、惡性腫瘤、抑郁癥等［5］。不僅會(huì)造成患者身體機(jī)能、生活質(zhì)量和社會(huì)參與度下降，也給患者家庭和社會(huì)帶來(lái)了巨大的經(jīng)濟(jì)負(fù)擔(dān)。迄今為止，其病因及具體發(fā)病機(jī)制尚未完全明確。故RA 早期的診斷具有重要意義。部分生物標(biāo)志物已被常規(guī)用于RA 的診斷和治療［6］，但仍不能滿足臨床實(shí)踐和藥物開(kāi)發(fā)的需求。有研究發(fā)現(xiàn)，免疫功能紊亂與RA 的發(fā)病密切相關(guān)［7］。以T 淋巴細(xì)胞、B 淋巴細(xì)胞、中性粒細(xì)胞、巨噬細(xì)胞為代表的免疫細(xì)胞參與并介導(dǎo)RA 的自身免疫反應(yīng)，誘導(dǎo)細(xì)胞因子、趨化因子擾亂免疫系統(tǒng)平衡，進(jìn)而加快軟骨和骨侵蝕［8］。因此，從免疫細(xì)胞浸潤(rùn)的角度研究RA的發(fā)病機(jī)制是目前RA 研究的熱點(diǎn)和發(fā)展趨勢(shì)［9］，對(duì)于開(kāi)發(fā)新的免疫治療靶點(diǎn)具有重要價(jià)值?；诖?，本研究中利用生物信息學(xué)的方法探索RA 的生物標(biāo)志物，并預(yù)測(cè)可能用于治療RA 的中藥，為今后診斷和治療RA 提供新思路。現(xiàn)報(bào)道如下。

1 材料與方法

1.1 儀器、試藥與細(xì)胞

儀器：DS-11FX+超微量/熒光分光光度計(jì)（美國(guó)Denovix公司）。

試藥：Trizol 試劑盒（天根生化科技＜北京＞有限公司，批號(hào)為DP424）；TaKaRa PrimeScriptTM反轉(zhuǎn)錄聚合酶鏈反應(yīng)（RT - PCR）試劑盒（日本TaKaRa 株式會(huì)社，批號(hào)為RR092A）；蛋白酶抑制劑（批號(hào)為ST506）、RIPA裂解緩沖液（批號(hào)為P0013B）、甘油醛-3-磷酸脫氫酶（GAPDH，批號(hào)為A0208），均購(gòu)自美國(guó)Beyotime Biotechnology 公司；BCA 試劑盒（廣州碩譜生物科技有限公司，批號(hào)為BL521A）；Dulbecco改良Eagle培養(yǎng)基（美國(guó)Gibco公司，批號(hào)為C11995500BT）；兔抗細(xì)胞骨架相關(guān)蛋白2（CKAP2）多克隆抗體（英國(guó)Abcam 公司，批號(hào)為ab227214）。

細(xì)胞：人成纖維細(xì)胞樣滑膜細(xì)胞系（MH7A），購(gòu)自中國(guó)科學(xué)院上海細(xì)胞生物研究所。

1.2 患者（基因）收集

選擇GPL96 平臺(tái)的3 個(gè)數(shù)據(jù)集GSE55235，GSE55457，GSE55584 收集患者滑膜組織基因表達(dá)情況。共收集到33例RA患者，20例健康者，分別納入RA組和HC組。其中3個(gè)數(shù)據(jù)集分別含RA患者、健康者各10例；RA患者13例，健康者10例；RA患者10例。

1.3 基因篩選與分析

數(shù)據(jù)獲取與差異表達(dá)基因的篩選：在Gene Expression Omnibus（GEO）數(shù)據(jù)庫(kù)以“rheumatoid arthritis”“synovial”為關(guān)鍵詞檢索人類mRNA 表達(dá)譜。以前述3 個(gè)數(shù)據(jù)集聯(lián)合分析，采用GPL8300 平臺(tái)的GSE10500 數(shù)據(jù)集和GPL570 平臺(tái)的GSE77298 數(shù)據(jù)集（分別采自滑液的巨噬細(xì)胞和滑膜組織分別含RA 患者5 例、健康者3 例和RA 患者16 例、健康者7 例）進(jìn)行外部驗(yàn)證。利用limma R 包以adjustedP＜0.05，|log2FC|＞1 為條件篩選差異表達(dá)基因（DEGs）。

加權(quán)基因共表達(dá)網(wǎng)絡(luò)（WGCNA）的構(gòu)建：計(jì)算整合后的數(shù)據(jù)集中每個(gè)基因的中位數(shù)絕對(duì)偏差，選取其中最大的5 000 個(gè)基因構(gòu)成表達(dá)矩陣。利用WGCNA 包進(jìn)行層級(jí)聚類，計(jì)算軟閾值，并構(gòu)建共表達(dá)網(wǎng)絡(luò)。

核心基因與生物標(biāo)志物的篩選：將與臨床癥狀相關(guān)性最高的模塊基因?qū)隒ytoscape 3.7.1 軟件，利用MCODE 插件篩選得分最高模塊。將模塊中的基因進(jìn)行單因素邏輯回歸，篩選P＜0.05 的基因?yàn)楹诵幕?。以?shù)量7∶3隨機(jī)分為訓(xùn)練集和測(cè)試集。訓(xùn)練集使用隨機(jī)森林遞歸特征消除（RF-RFE）、支持向量機(jī)遞歸特征消除（SVM-RFE）2種算法進(jìn)行特征選擇，重疊基因即為RA的生物標(biāo)志物。使用GSE10500 和GSE77298 進(jìn)行外部驗(yàn)證。

基因集富集分析（GSEA）：使用clusterProfiler 包，分別以“c5. go. bp. v7.4. entrez. gmt”“c2. cp. kegg. v7.4.entrez. gmt”為參考基因集對(duì)RA 和生物標(biāo)志物進(jìn)行GSEA，以P＜0.05和adjustedP＜0.05篩選顯著通路。

免疫細(xì)胞浸潤(rùn)的評(píng)估：使用CiberSort算法獲取22種免疫細(xì)胞浸潤(rùn)水平，比較RA 組與HC 組免疫細(xì)胞浸潤(rùn)的差異；從InnateDB數(shù)據(jù)庫(kù)（https：//www.innatedb.ca/）中獲取先天免疫基因，進(jìn)行相關(guān)性分析。

1.4 mRNA 表達(dá)水平檢測(cè)

采用RT-PCR 法。取對(duì)數(shù)生長(zhǎng)期的MH7A 細(xì)胞，用含10%胎牛血清、100 U/mL青霉素和100 mg/mL鏈霉素的Dulbecco 改良Eagle 培養(yǎng)基，于37 ℃、5% CO2條件下培養(yǎng)，4 ℃下12 000 r/min 離心15 min，提取總RNA，并測(cè)定RNA 的濃度和純度。反轉(zhuǎn)錄得到cDNA，反應(yīng)體系20 μL，循環(huán)程序：94 ℃5 min，94 ℃30 s，57 ℃30 s，40 個(gè)循環(huán)；72 ℃30 s。引物序列：CKAP2為5'-AGGCTGTTACAACATAAGAG - 3'和5'- GAAACTGACAAGGTACGATTAG - 3'；GAPDH為5' - CACCCACTCCTCCACCTTTG-3'和5'-CCACCACCCTGTTGCTGTAG -3'。以GAPDH 為內(nèi)參，通過(guò)2-ΔΔCT分析mRNA 的相對(duì)表達(dá)水平。所有樣品重復(fù)3次。

1.5 蛋白表達(dá)水平檢測(cè)

采用Western blot 法。收集MH7A 細(xì)胞，裂解得上清液，使用BCA 試劑盒測(cè)定蛋白質(zhì)濃度。將20 μL 上清液與5 × loading buffer 混合，98 ℃煮沸，12 000 r / min離心1 min，12% SDS - PAGE 凝膠電泳，轉(zhuǎn)移至PVDF膜，5%脫脂乳（以含吐溫-20的Tris緩沖鹽溶液為溶劑）封閉1 h，4 ℃兔抗CKAP2多克隆抗體（1∶1 000，V/V）孵育過(guò)夜，二抗孵育2 h，檢測(cè)。

1.6 相關(guān)中藥預(yù)測(cè)

從TCMSP 數(shù)據(jù)庫(kù)（https：// tcmspw. com/）中獲取所有中藥的活性成分及其對(duì)應(yīng)基因，將其與DEGs的重疊基因作為中藥治療RA的靶基因，統(tǒng)計(jì)中藥對(duì)應(yīng)的靶基因數(shù)量。篩選類藥性（DL）＞0.18的活性成分，以及靶基因數(shù)量大于上四分位數(shù)的中藥，利用Cytoscape 3.7.1 構(gòu)建中藥-活性成分-靶基因網(wǎng)絡(luò)圖，使用MCODE 插件篩選評(píng)分最高的模塊，其包含的中藥可能是用于治療RA的中藥。

1.7 數(shù)據(jù)分析

采用GraphPad Prism 5.0 軟件進(jìn)行統(tǒng)計(jì)學(xué)分析。計(jì)量資料以X±s表示，行t檢驗(yàn)。P＜0.05 為差異有統(tǒng)計(jì)學(xué)意義。

2 結(jié)果

2.1 差異表達(dá)基因分析

共有706個(gè)差異表達(dá)基因，其中包括上調(diào)基因410個(gè)，下調(diào)基因296個(gè)。詳見(jiàn)圖1（火山圖中，紅色表示上調(diào)，藍(lán)色表示下調(diào)，灰色表示不顯著）。

A. 主成分分析聚類圖 B. 差異表達(dá)基因火山圖 C. 顯著差異表達(dá)基因熱圖圖1 差異表達(dá)基因分析A.Cluster diagram of principal component analysis B.Volcano plot of differentially - expressed genes C.Heatmap of significantly differentially -expressed genesFig.1 Analysis of differentially - expressed genes

2.2 GSEA

基因本體論（GO）功能富集中的生物過(guò)程（BP）分析顯示，RA主要富集在B淋巴細(xì)胞受體信號(hào)通路、抗原受體介導(dǎo)的信號(hào)通路、免疫反應(yīng)調(diào)節(jié)細(xì)胞表面受體信號(hào)通路、免疫球蛋白的產(chǎn)生等生物學(xué)途徑；京都基因與基因組百科全書(shū)（KEGG）通路富集分析顯示，RA 主要富集在原發(fā)性免疫缺陷、免疫球蛋白A（IgA）產(chǎn)生的腸道免疫網(wǎng)絡(luò)、抗原加工和呈遞、自身免疫性甲狀腺疾病等信號(hào)通路。詳見(jiàn)圖2。

A，C. 生物學(xué)過(guò)程 B，D.KEGG通路圖2 基因集富集分析A，C.Biological processes B，D.KEGG pathwaysFig.2 Analysis of gene set enrichment

2.3 WGCNA

共得到5 個(gè)模塊。模塊與臨床癥狀相關(guān)性分析顯示，綠松石（顏色）模塊中的基因不僅與模塊相關(guān)度高，而且與RA也高度相關(guān)（r=0.87，P＜0.001）。詳見(jiàn)圖3。

A. 不同軟閾值下的擬合優(yōu)度 B. 模塊與RA相關(guān)性熱圖 C. 聚類分析樹(shù)狀圖 D. 綠松石模塊特征散點(diǎn)圖圖3 加權(quán)基因共表達(dá)網(wǎng)絡(luò)分析A.Goodness of fit under different soft thresholds B.Heatmap of correlation between module and RA C.Tree diagram of cluster analysis D.Characteristic scatter diagram of turquoise moduleFig.3 Analysis of weighted gene co - expression network

2.4 核心基因與潛在生物標(biāo)志物的篩選與驗(yàn)證

獲得32 個(gè)核心基因。通過(guò)2 種算法篩選得到重疊基因CKAP2（見(jiàn)圖4）。使用受試者操作特征（ROC）曲線檢驗(yàn)CKAP2的診斷效能。訓(xùn)練集中CKAP2的藥- 時(shí)曲線下面積（AUC）為0.994，測(cè)試集中為0.963。與HC 組比較，RA 組患者CKAP2表達(dá)水平顯著上調(diào)（見(jiàn)圖5）。外部驗(yàn)證其AUC分別為0.933 和0.866，表明CKAP2具有較高的診斷價(jià)值。BP 的GSEA 分析顯示，CKAP2主要富集于T 淋巴細(xì)胞選擇、B 淋巴細(xì)胞受體信號(hào)通路、肽抗原的抗原加工和呈遞等生物途徑；KEGG 的GSEA 分析顯示，CKAP2主要富集于同種異體排斥、原發(fā)性免疫缺陷、IgA 產(chǎn)生的腸道免疫網(wǎng)絡(luò)、抗原加工和呈遞、系統(tǒng)性紅斑狼瘡和細(xì)胞黏附分子（CAM）等信號(hào)通路。

2.5 免疫浸潤(rùn)分析

與HC 組比較，RA 組患者M(jìn)1 巨噬細(xì)胞、漿細(xì)胞、活化的記憶T淋巴細(xì)胞CD4+、T淋巴細(xì)胞CD8+、濾泡輔助T 淋巴細(xì)胞等浸潤(rùn)較多，而在靜息的樹(shù)狀細(xì)胞、活化的自然殺傷（NK）細(xì)胞、活化的肥大細(xì)胞、靜息的記憶T 淋巴細(xì)胞CD4+等浸潤(rùn)較少。詳見(jiàn)圖6。CKAP2表達(dá)量與免疫細(xì)胞相關(guān)性分析顯示，CKAP2的表達(dá)與漿細(xì)胞、T 淋巴細(xì)胞CD8+呈顯著正相關(guān)（r= 0.84，P= 5.0 × 10-15；r= 0.70，P= 4.1 × 10-9），與靜息的記憶T 淋巴細(xì)胞CD4+、活化的肥大細(xì)胞呈顯著負(fù)相關(guān)（r= - 0.69，P=1.5×10-8；r=-0.56，P=1.2×10-5）。詳見(jiàn)圖7。

從InnateDB數(shù)據(jù)庫(kù)獲取先天免疫基因，分析CKAP2與先天免疫基因的相關(guān)性。結(jié)果顯示，CKAP2的表達(dá)與CD27，PLCG2，UCP2呈顯著正相關(guān)，與IL1R1，EGFR，PARD3呈顯著負(fù)相關(guān)（P＜0.05）。其中，CKAP2與CD27表達(dá)的相關(guān)性最高（r=0.92，P＜0.001）。

2.6 CKAP2 體外表達(dá)水平的驗(yàn)證

與不含白細(xì)胞介素1β（IL - 1β）的MH7A 細(xì)胞比較，含IL-1β的MH7A 細(xì)胞中CKAP2的mRNA、蛋白表達(dá)水平均顯著升高（P＜0.05）。詳見(jiàn)圖8。與前述分析結(jié)果一致，表明數(shù)據(jù)分析中鑒定出的RA 潛在生物標(biāo)志物CKAP2在MH7A細(xì)胞中變化顯著。

2.7 相關(guān)中藥預(yù)測(cè)

共篩選出108種中藥，構(gòu)建中藥-活性成分-靶基因網(wǎng)絡(luò)圖，使用MCODE插件篩選高得分模塊后得到11種中藥（包括蒼耳子、黃芩、黨參等）。詳見(jiàn)圖9。

A. 中藥篩選 B. 中藥- 活性成分- 靶基因網(wǎng)絡(luò) C.MCODE得分最高的模塊圖9 中藥預(yù)測(cè)A.Screening of TCMs B.Network of TCMs - active components - targets C.Modules with the highest score of MCODEFig.9 Prediction of TCMs

3 討論

RA 病因復(fù)雜，由于關(guān)節(jié)破壞的不可逆性，早期診斷對(duì)于預(yù)防RA 進(jìn)展有重要作用。目前已有的生物標(biāo)志物不能滿足臨床需求，仍有部分RA患者漏診［10］，因此，尋找新的、更有效的生物標(biāo)志物，分析RA 免疫細(xì)胞浸潤(rùn)情況，預(yù)測(cè)RA 新的治療藥物，對(duì)于RA 患者的診斷與治療具有重要意義。

本研究中發(fā)現(xiàn)，CKAP2 可能為RA 的潛在生物標(biāo)志物，且篩選出的CKAP2在內(nèi)外部驗(yàn)證中均可靠，表明本研究的整合策略是可行的。CKAP2 為微管結(jié)合蛋白，具有穩(wěn)定微管的作用，可調(diào)節(jié)非整倍體、細(xì)胞周期和細(xì)胞凋亡［11］。既往研究表明，CKAP2 在乳腺癌［12］、卵巢癌［13］、胃癌［14］及其他惡性腫瘤［15］中可通過(guò)各種機(jī)制縮短患者的總生存期。然而，其在RA中的表達(dá)水平與作用尚無(wú)研究報(bào)道。有文獻(xiàn)報(bào)道，肝癌細(xì)胞中CKAP2基因沉默顯著降低了其細(xì)胞內(nèi)JAK2 和STAT3 的磷酸化水平，抑制了JAK/STAT 信號(hào)通路的激活［16］。動(dòng)物實(shí)驗(yàn)表明，磷酸化JAK2 和STAT3 的表達(dá)與RA 模型大鼠損傷程度高度相關(guān)［17］。此外，RA中JAK/STAT信號(hào)通路持續(xù)激活可導(dǎo)致基質(zhì)金屬蛋白酶（MMP）基因表達(dá)水平升高，而MMP 可降解關(guān)節(jié)細(xì)胞外基質(zhì)成分，導(dǎo)致關(guān)節(jié)軟骨病變［18］。且激活JAK/STAT信號(hào)通路可阻礙成纖維樣滑膜細(xì)胞（FLS）生長(zhǎng)而促進(jìn)其凋亡，導(dǎo)致滑膜細(xì)胞增生［19］。

本研究中發(fā)現(xiàn)，CKAP2在RA 患者中的表達(dá)上調(diào)，推測(cè)其可能參與了RA 的病理過(guò)程。GSEA 分析顯示，RA 的發(fā)生及CKAP2高表達(dá)主要與T 淋巴細(xì)胞、樹(shù)狀細(xì)胞遷移與趨化、B淋巴細(xì)胞、抗原加工、CAM及腸道分泌型IgA 免疫作用有關(guān)。既往研究表明，樹(shù)突狀細(xì)胞、B 淋巴細(xì)胞可介導(dǎo)T 淋巴細(xì)胞活化，活化的T 淋巴細(xì)胞分泌促炎細(xì)胞因子引起滑膜炎癥，導(dǎo)致關(guān)節(jié)病理性損傷，促進(jìn)RA 疾病的進(jìn)展［20］。CAM 在RA 的發(fā)病中起關(guān)鍵作用，其隨白細(xì)胞浸潤(rùn)到滑膜室中，滑膜微血管中的內(nèi)皮細(xì)胞激活，CAM表達(dá)量增加［21］。這些研究結(jié)果進(jìn)一步表明CKAP2 與RA 進(jìn)展有關(guān)，且可能作為其潛在的生物標(biāo)志物。

通過(guò)分析CKAP2與免疫細(xì)胞和先天免疫基因的相關(guān)性，發(fā)現(xiàn)CKAP2的表達(dá)與漿細(xì)胞、T 淋巴細(xì)胞CD8+和CD27呈顯著正相關(guān)，而與靜息記憶T 淋巴細(xì)胞CD4+、活化的肥大細(xì)胞呈顯著負(fù)相關(guān)。既往研究表明，漿細(xì)胞、T 淋巴細(xì)胞、肥大細(xì)胞和自然殺傷細(xì)胞在RA 中起重要作用［20，22］。此外，CD27已被證明可防止T 淋巴細(xì)胞凋亡，并在RA 的發(fā)生和發(fā)展中發(fā)揮重要作用［23］。本研究結(jié)果表明，CKAP2 與CD27呈強(qiáng)相關(guān)，因此推測(cè)CKAP2通過(guò)上調(diào)漿細(xì)胞、T 淋巴細(xì)胞CD8+、CD27或下調(diào)靜息的記憶T 淋巴細(xì)胞CD4+、活化的肥大細(xì)胞來(lái)參與RA 的發(fā)生、發(fā)展。并驗(yàn)證了CKAP2 在MH7A 細(xì)胞中的mRNA 和蛋白表達(dá)水平與生物信息學(xué)分析結(jié)果一致，進(jìn)一步證明了本次分析的可靠性和CKAP2對(duì)RA發(fā)展的重要性。此外，為尋找治療RA 的潛在中藥，本研究中利用TCMSP 數(shù)據(jù)庫(kù)獲得了DEGs 相關(guān)的中藥11 種。這些中藥未來(lái)均可能被開(kāi)發(fā)為治療RA的新藥。

但本研究也存在一定不足，如初步篩選的所有數(shù)據(jù)均來(lái)源于公共數(shù)據(jù)庫(kù)，屬回溯性研究，缺乏對(duì)其在RA中的作用及其他基因間詳細(xì)關(guān)系的具體研究。

綜上所述，CKAP2 可作為RA 早期診斷的候選生物標(biāo)志物，蒼耳子、黃芩、黨參等11 種中藥對(duì)RA 具有潛在治療作用。本研究可為治療RA 中藥的研發(fā)提供新思路。