miR-27a在原發性痛風關節炎患者中的表達水平及生物信息學分析*

茍海梅，陳 瑩，徐 磊，李九龍，郭曉蘭，，鐘曉武，△

1．川北醫學院附屬醫院檢驗科，四川南充 637000；2． 川北醫學院檢驗系，四川南充 637000； 3．川北醫學院轉化醫學研究中心，四川南充 637000；4．南充市第五人民醫院檢驗科，四川南充 637100

痛風是機體嘌呤代謝紊亂、血尿酸水平持續升高而導致單鈉尿酸鹽(MSU)晶體析出并沉積于組織或器官而引起反復發作性炎性疾病[1]。近年來，隨著生活水平不斷提高，痛風的發病率逐年呈上升趨勢，已成為臨床常見病、多發病，且痛風發病年齡有向低齡化發展趨勢[2]。MSU晶體的沉積是痛風發生的重要物質基礎，但痛風發病機制及炎癥自限的確切分子機制尚未完全闡明，仍是目前亟待解決的問題。

微小RNA(miRNA)是一類由22～24 個核苷酸組成的高度保守性內源性非編碼小RNA，它通過與mRNA的3′非翻譯區(3′UTR)特異性結合降解靶mRNA和(或)抑制靶mRNA的表達及翻譯，進行轉錄后調控[3]。多項研究表明，miRNA在痛風發病機制中有關鍵性調控作用[4-6]。miR-27a在骨性關節炎[7]及系統性紅斑狼瘡[8]中均表達異常呈下調趨勢，提示miR-27a與炎性反應負調控相關。已有研究表明，miR-27a在痛風患者中存在異常表達，其可能作為炎癥負性調控因子參與痛風的炎癥免疫反應[9]，但其與痛風發生相關的確切靶分子機制尚不清楚。本研究運用生物信息學分析方法預測miR-27a作用的靶基因，對其靶基因進行GO富集分析、KEGG信號通路富集分析，為痛風炎癥免疫反應研究提供理論依據及數據支持。

1 資料與方法

1.1一般資料 將川北醫學院附屬醫院風濕免疫科就診的20例男性原發性痛風性關節炎(PGA)患者納入PGA組，將20例體檢健康男性納入健康對照組。其中，PGA組平均年齡(44.30±13.27)歲，健康對照組平均年齡(47.70±12.65)歲，兩組年齡比較，差異無統計學意義(P>0.05)。

1.2方法

1.2.1miRNA檢測 采用人外周血淋巴細胞分離液(天津灝洋公司)對人外周血單個核細胞(PBMC)進行分離，并用Trizol試劑盒(美國 Invitrogen公司)一步法提取人外周血單個核細胞總RNA，同時用紫外分光光度儀進行濃度測定。采用TaqMan MicroRNA Assays 試劑盒(美國Life Technologies公司)將miR-27a逆轉錄為 cDNA，并采用實時熒光定量PCR進行擴增及熒光檢測，以 U6為內參，嚴格按照說明書指示完成檢測操作。采用2-ΔΔCt法分析miRNA的相對表達量。

1.2.2miR-27a在人類基因組中的位置及保守性分析 應用UCSC基因組在線軟件工具(http://genome.ucsc.edu/)分析miR-27a在人類基因組中的位置，查找不同物種間的成熟序列并分析其保守性。

1.2.3miR-27a靶基因篩選 應用TargetScan (http://www.targetscanorg/)數據庫，miRDB(http://mirdb.org/)數據庫和RNA22(https://cm.jefferson.edu/rna22/precomputed/)數據庫分別預測 miR-27a 靶基因，同時運用GeneCards數據庫(http://www.genecards.org/)查詢與痛風發病相關的疾病基因。為了增加預測結果的可靠性，將4種數據庫的預測結果取交集作為與痛風發病相關的miR-27a靶基因。

1.2.4miR-27a靶基因的功能注釋 應用DAVID數據庫(https://david.ncifcrf.gov/)對miR-27a預測靶基因集合進行GO富集分析和KEGG通路富集分析。GO富集分析包括生物過程、細胞組分和分子功能。

2 結 果

2.1miR-27a在人類基因組中的位置及保守性分析 應用UCSC基因組在線軟件獲取miR-27a在人類基因組中的位置(圖1)，miR-27a定位于人類19q13.12染色體上chr19:13 836 440-13 836 517位置，長度共78 bp，其核苷酸序列在人、黑猩猩、大猩猩、獼猴、狒狒、小鼠、狗、大象物種中呈高度保守。

圖1 miR-27a 在基因組中的位置和保守性分析

2.2miR-27a 在PGA患者中的表達 應用實時熒光定量PCR法檢測PGA組和健康對照組PBMC中miR-27a的相對表達水平，結果顯示，miR-27a在PGA組患者中的相對表達水平明顯低于健康對照組，差異有統計學意義(P<0.05)。見圖2。

圖2 miR-27a 在PGA組和健康對照組中的表達水平

2.3篩選與痛風發病相關的miR-27a靶基因 應用TargetScan、miRDB和RNA22數據庫分別預測miR-27a的靶基因，然后對miR-27a靶基因結果取交集，共得到 585個潛在調控靶基因。此外，在GeneCards數據庫中查詢與痛風發病相關基因共980個基因。將GeneCards數據庫所得基因與585個潛在調控靶基因再次取交集，最終得到MAPK14、RFX3、ABCA1、HS2ST1、ENPEP和CACNA2D3等共33個與痛風發病相關的miR-27a靶基因。見圖3。

2.4預測靶基因的GO富集分析 對與痛風發病相關的miR-27a靶基因進一步進行GO富集分析，GO富集分析結果顯示：miR-27a的靶基因顯著富集到生物過程調節、細胞過程調節、有機物代謝過程、多細胞生物過程、生物過程正調控、對外界刺激的反應、主要代謝調節等多個生物學過程中；富集到質膜的外側、細胞表面、血影蛋白細胞骨架及細胞膜的細胞組分；富集有機環化合物結合、蛋白結合、催化活性、離子結合等分子功能。見表1。

圖3 miR-27a靶基因預測

2.5預測靶基因的KEGG信號通路分析 對與痛風發病相關的miR-27a靶基因進一步進行KEGG功能富集分析，KEGG信號通路分析顯示，miR-27a的靶基因顯著富集于MAPK信號通路、癌癥相關通路、Rap1信號通路、細胞因子-細胞因子受體相互作用、干細胞多能性調節信號通路、青少年成人起病型糖尿病通路和醛固酮控制的鈉鹽重吸收通路。見表2。

表1 miR-27a 預測靶基因 GO 富集分析結果

續表1 miR-27a 預測靶基因 GO 富集分析結果

表2 miR-27a預測靶基因 KEGG 信號通路分析

3 討 論

miRNA作為調控分子，可同時調節多個靶基因，參與生物體細胞信號轉導通路調節，從而對細胞特定功能和表型進行有效調控，發揮生物學功能[3]。現有研究表明，miRNA與免疫功能失調、多種慢性炎癥或自身免疫病等相關[10-11]。通過UCSC基因組在線軟件分析，miR-27a核苷酸序列在多個物種中呈高度保守，提示其發生突變的概率低，有潛在的重要生物學功能。本研究對PGA組和健康對照組中miR-27a基因表達進行了檢測，發現在PGA組中miR-27a基因表達水平較健康對照組明顯降低，與李寧寧等[9]研究結果一致，進一步提示miR-27a可能作為炎癥負性調控因子在痛風的炎癥免疫反應過程中起到重要作用。

為進一步分析miR-27a參與痛風炎癥免疫反應相關的靶分子機制，準確預測其靶基因，本研究運用TargetScan、RNA22和miRDB 3個使用較多的靶基因預測軟件，取3個軟件共同預測的靶基因作為結果，利用不同軟件、不同算法從多角度分析，以增強預測結果的可靠性。本研究再將生物信息學方法預測的靶基因和 GeneCards數據庫中與痛風相關的基因聯合分析，取交集最終得到MAPK14、RFX3、ABCA1、HS2ST1、ENPEP、CACNA2D3等33個潛在的靶基因。MAPK14(即p38α)是p38MAPK最重要的亞型之一，p38MAPK信號通路與炎性反應調控密切相關。有研究顯示，p38MAPK信號通路參與調控MSU誘導的單核細胞IL-8表達，進而參與痛風性關節炎的炎性反應[12]。在中國人群全基因組關聯分析(GWAS)研究中發現，RFX3是中國人群特有的原發性痛風易感基因[13]。RFX3由胰島β細胞分泌，可結合到葡萄糖激酶基因的啟動子區域，調控其表達，從而影響胰島β細胞的功能[14]。胰島β細胞主要參與血糖的調節，糖代謝的紊亂不僅直接影響嘌呤的代謝及尿酸的合成，而且影響尿酸的排泄[15]。因此，RFX3突變可能引起胰島β細胞功能缺陷，從而直接或間接引起高尿酸血癥和痛風。

GO富集分析表明，miR-27a靶基因存在細胞代謝調節及對外界刺激的反應等多種生物學功能。痛風屬于代謝性風濕病范疇，與嘌呤代謝紊亂及尿酸排泄減少所致的高尿酸血癥相關，高尿酸血癥所致的MSU沉積可作為外界刺激反復引起炎性反應。KEGG通路富集分析顯示，miR-27a可能通過調控MAPK信號通路、細胞因子-細胞因子受體相互作用和醛固酮控制的鈉鹽重吸收通路等發揮生物學功能。MAPK信號通路是一類蘇氨酸蛋白激酶，可以被細胞外信號或刺激等激活，包括ERK1/2 MAPK、JNK MAPK、p38MAPK等活化途徑。目前研究發現，ERK1/2 MAPK、JNK MAPK和p38MAPK通路可誘導前炎癥細胞因子IL-8的產生，且參與調控MSU誘導的單核細胞IL-8表達[14，16-17]。在體外MSU誘導的單核細胞模型中發現，中性粒細胞趨化活性90%來自IL-8，而大量中性粒細胞浸潤到關節腔并介導炎性反應是痛風性關節炎的核心病理機制[14]。MAPK14是miR-27a預測靶基因，故筆者推測痛風性關節炎經MSU刺激后，miR-27a基因表達水平降低，且作為炎癥負性調控因子靶向作用于MAPK14，MAPK14在痛風性關節炎呈高表達，MAPK14經p38MAPK通路誘導單核細胞IL-8表達，IL-8激活中性粒細胞并浸潤到關節腔介導痛風炎性反應。近年來的研究發現，IL-1、IL-8、IL-10等細胞因子在痛風性關節炎的發生、發展過程中發揮重要作用[14]，miR-27a參與細胞因子-細胞因子受體相互作用，因此miR-27a可能通過細胞因子-細胞因子受體相互作用參與痛風性關節炎的發病機制。醛固酮控制的鈉鹽重吸收系統失衡會引起腎臟的損傷，進而影響嘌呤的代謝及尿酸的合成，而且影響尿酸的排泄，miR-27a也可能通過參與醛固酮控制的鈉鹽重吸收直接或間接引起高尿酸血癥和痛風。

本研究通過生物信息學分析挖掘miR-27a在痛風炎癥免疫反應相關的靶分子機制，結果表明miR-27a可能通過調控多個與痛風發病相關的靶基因，影響多條信號通路的網絡調節，從而參與痛風炎癥免疫反應機制，進一步為miR-27a調控痛風炎癥免疫反應機制及痛風治療靶點提供了新的依據。