周新朋 金曄華 張潤潤 何東儀

1.上海中醫(yī)藥大學，上海 201203；2.上海中醫(yī)藥大學附屬光華醫(yī)院關節(jié)內(nèi)科，上海 200052

類風濕關節(jié)炎（rheumatoid arthritis，RA）是一種慢性炎癥性疾病，以關節(jié)受累為主，引起關節(jié)軟骨和骨的破壞，并最終導致關節(jié)的畸形[1]。其發(fā)病與遺傳、性別、年齡、吸煙等密切相關[2]。成纖維樣滑膜細胞（fibroblast-like synoviocytes，F(xiàn)LS）是間充質(zhì)來源的、具有成纖維細胞特征的滑膜細胞，是RA 發(fā)病過程中的主要的效應細胞[3]。

人類基因組絕大部分都可以被轉(zhuǎn)錄，但只有一小部分（大約2%）轉(zhuǎn)錄為mRNA，剩余部分則轉(zhuǎn)錄為不同種類的非編碼RNA[4]。長鏈非編碼RNA（lncRNA）是指長度超過200 nt，但不具有蛋白編碼能力的RNA[5]。lncRNA 可通過競爭內(nèi)源性RNA（competitive endogenous RNA，ceRNA）機制，通過微小RNA（miRNA）反應元件（microRNA response elements，MREs）競爭性結合miRNA，對基因的表達進行轉(zhuǎn)錄后調(diào)節(jié)[6]。

本課題組在美國國立生物技術信息中心的基因表達數(shù)據(jù)庫（Gene Expression Omnibus，GEO）中獲得RA-FLS 細胞的lncRNA 表達芯片數(shù)據(jù)集，將多個數(shù)據(jù)集進行合并，構建lncRNA-miRNA-mRNA網(wǎng)絡并預測關鍵的lncRNA。

1 材料與方法

1.1 基因表達芯片數(shù)據(jù)集的獲取

在GEO 中獲取相關的數(shù)據(jù)集GSE103578、GSE8 3147、GSE128813，每個數(shù)據(jù)集中對照組含有3 個健康人滑膜分離的FLS，RA 組含有3 個RA 患者滑膜分離的FLS。使用limma 包[7]進行標準化，通過idmap3 包[8]對探針進行注釋，根據(jù)基因名稱將數(shù)據(jù)集合并，并去除批次效應。

1.2 差異基因的篩查

根據(jù)基因類型篩選mRNA 及l(fā)ncRNA，通過limma包[7]分別進行基因差異表達分析，以logFC >1.0 并且P ＜0.05 的基因為差異表達基因。

1.3 構建lncRNA-miRNA-mRNA 網(wǎng)絡

在lncbase[9]、starBase[10]、miRcode[11]數(shù)據(jù)庫中，獲得與差異表達的lncRNA 相關miRNA，在miRWalk[12]、TargetScan[13]數(shù)據(jù)庫中尋找與mRNA 相關的miRNA，二者取交集，通過Cytoscape 軟件[14]，構建lncRNAmiRNA-mRNA 網(wǎng)絡，將網(wǎng)絡中的lncRNA 作為關鍵lncRNA。

2 結果

2.1 芯片數(shù)據(jù)預處理

對獲得的數(shù)據(jù)集進行注釋及合并，得到對照組、RA 組各9 個樣本，共有20 050 條基因的數(shù)據(jù)集，其中mRNA 有14 760 條、lncRNA 有5290 條。去除批次效應后各樣本間數(shù)據(jù)達到一致性要求（批次效應處理前后數(shù)據(jù)集均一性變化如圖1 所示），可以進行差異表達分析。

2.2 基因差異表達分析

分別對數(shù)據(jù)集中的mRNA 和lncRNA 進行差異表達分析。結果顯示，與對照組比較，RA 組mRNA 上調(diào)31 條，下調(diào)56 條，見圖2A；RA 組lncRNA上調(diào)9 條，下調(diào)14 條，見圖2B。

圖1 去除批次效應前后基因的表達情況

圖2 基因差異表達分析

2.3 lncRNA-mRNA 共表達分析

對差異表達的mRNA 和lncRNA 進行Pearson 相關性分析顯示，具有明顯相關性的lncRNA 有18 條，mRNA有51 條（R ＞0.85 且P ＜0.05）。構造的lncRNAmRNA共表達網(wǎng)絡見圖3。

2.4 通過數(shù)據(jù)庫構建lncRNA-miRNA-mRNA 網(wǎng)絡

通過獲取多個數(shù)據(jù)庫的RNA 信息，得到與共表達的lncRNA、mRNA 有結合位點的miRNA，并取交集，獲得相關miRNA，構建lncRNA-miRNA-mRNA 相關網(wǎng)絡，其中l(wèi)ncRNA 有13 條，包括lncRNA XIST、NR2F2-AS1、LINC01018 等，miRNA 有26 條，mRNA 有46 條。lncRNA-miRNA-mRNA 共表達網(wǎng)絡見圖4。

3 討論

圖3 lncRNA-mRNA 共表達網(wǎng)絡

FLS 是RA 血管翳中最常見的一種細胞類型，通過產(chǎn)生細胞因子、趨化因子、基質(zhì)金屬蛋白酶及侵襲關節(jié)軟骨導致關節(jié)的破壞。RA-FLS 與健康人滑膜FLS 比較，在細胞因子的表達、細胞信號通路的激活、原癌基因的表達方面存在明顯不同[15]，這些不同可能與RA-FLS 的表觀遺傳學改變有關[16-17]。

lncRNA 的作用機制多樣，可在轉(zhuǎn)錄和轉(zhuǎn)錄后改變基因的表達，調(diào)節(jié)各種生物學過程。最近發(fā)現(xiàn)可通過與miRNA 相互作用，改變miRNA 的表達水平，在轉(zhuǎn)錄后調(diào)節(jié)基因表達。lncRNA 與mRNA 競爭性結合MRE，可先通過計算機預測lncRNA-miRNA-mRNA 之間的相互關系[18]，為進一步實驗驗證提供指導。

有研究顯示XIST 是多種腫瘤的原癌基因[19]，本研究發(fā)現(xiàn)XIST 在基因表達芯片中差異表達，并且與多個miRNA、mRNA 相互作用，Wang 等[20]發(fā)現(xiàn)XIST 在骨性關節(jié)炎（OA）樣本中表達升高，XIST 及MMP-13、ADAMTS5 為miR-1277-5p 的靶基因，通過ceRNA 機制競爭性結合miR-1277-5p，導致MMP-13、ADAMTS5升高，促進細胞外基質(zhì)的降解，RA 與OA 發(fā)病存在共有機制[21]，XIST 的研究對RA 骨破壞機制的研究具有參考意義。NR2F2-AS1 在多種腫瘤中具有重要作用，NR2F2-AS1 在非小細胞肺癌組織及細胞中表達升高，通過與BMI1 競爭性結合miR-320b 而與腫瘤的TNM 分期、患者的淋巴結轉(zhuǎn)移相關[22]。LINC01018 在肝臟內(nèi)可以通過RNA 結合蛋白HuR 調(diào)節(jié)脂肪酸氧化相關基因的表達[23]。肝癌細胞中LINC01018 的DNA甲基化水平改變，導致在肝細胞癌中表達水平下降，抑制肝癌細胞增生減弱[24]。

圖4 lncRNA-miRNA-mRNA 共表達網(wǎng)絡

lncRNA 在RA 中的作用已在多項研究中被證實，但是主要集中在HOTAIR、LincRNA-p21 等被廣泛研究的lncRNA 中，其他lncRNA 在RA 中的作用還需要進一步尋找并驗證[25]。本研究發(fā)現(xiàn)關鍵lncRNA 有13 條，可以構建復雜的lncRNA-miRNA-mRNA 網(wǎng)絡，進一步揭示了lncRNA 在RA-FLS 細胞中的功能機制。