鄭文新 劉 紅 汪 偉

摘? ?要：目的? 分析銀屑病易感基因的功能，篩選關(guān)鍵基因。方法? 結(jié)合OMIM數(shù)據(jù)庫以及近年來發(fā)表的文章收集銀屑病易感基因，利用Gather在線服務(wù)對(duì)易感基因進(jìn)行富集分析，蛋白質(zhì)相互作用網(wǎng)絡(luò)由String數(shù)據(jù)庫構(gòu)建，計(jì)算節(jié)點(diǎn)的度數(shù)，據(jù)此篩選關(guān)鍵易感基因。結(jié)果? 銀屑病的易感基因功能主要富集在免疫應(yīng)答，防御反應(yīng)，生物應(yīng)激反應(yīng)，有機(jī)生理過程，應(yīng)激反應(yīng)，對(duì)害蟲寄生蟲病原體的反應(yīng)，對(duì)外部生物刺激的反應(yīng)，對(duì)外傷的反應(yīng)，抗原表達(dá)，抗原呈遞及抗原加工等類別。我們通過計(jì)算易感基因所編碼的蛋白質(zhì)相互作用網(wǎng)絡(luò)中節(jié)點(diǎn)的度數(shù)，按度數(shù)的高低，選出排在前十的10個(gè)基因，作為關(guān)鍵易感基因，功能分析和通路分析提示關(guān)鍵基因和很多功能明確與銀屑病相關(guān)的基因相互作用。結(jié)論? 易感基因分析提示，銀屑病是與環(huán)境因素密切相關(guān)，具有遺傳性的免疫性皮膚病。蛋白質(zhì)相互作用網(wǎng)絡(luò)分析提示，關(guān)鍵易感基因可以沒有明確的功能，但是會(huì)和很多與銀屑病相關(guān)功能明確的基因相互作用。

關(guān)鍵詞：銀屑病;易感基因;蛋白質(zhì)相互作用網(wǎng)絡(luò);富集分析

中圖分類號(hào)：R758.63? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? 文獻(xiàn)標(biāo)識(shí)碼：B? ? ? ? ? ? ? ? ? ? ? ? ? ? ? DOI：10.3969/j.issn.1006-1959.2018.21.001文章編號(hào)：1006-1959（2018）21-0001-03

Bioinformatics Analysis of the Susceptible Genes of Psoriasis

ZHENG Wen-Xin1，LIU Hong1，WANG Wei2

（1.Capital Medical University School of Biomedical Engineering，Beijing 100069，China;

2.Central Laboratory of Capital Medical University，Beijing 100069，China）

Abstract：Objective? To analyze the functional enrichment of psoriasis susceptible genes and screen the key susceptible genes.Methods? According to the OMIM database and the published articles in recent years，the susceptibility genes of psoriasis were collected，and the susceptibility genes were enriched and analyzed by Gather online service. The protein interaction network was constructed by using the String database， and the key susceptible genes were screened. Results? The function of susceptibility gene of psoriasis is mainly enriched in immune response， defense response， biological stress response， organic physiological process， stress response， response to insect parasite pathogens， response to external biological stimuli， response to trauma， antigen expression， antigen presentation and antigen processing. By calculating the degree of the nodes in the protein interaction network encoded by the susceptibility genes， we selected the first 10 genes， as the key susceptible genes. Both the functional analysis and pathway analysis suggested that the key genes interact with the many genes related to psoriasis. Conclusion? Susceptibility gene analysis indicates that psoriasis is a hereditary immune skin disease， closely related to environmental factors. Protein interaction network analysis suggests that though some key susceptible genes do not have definite functions， but they interact with a lot of genes with specific functions associated with psoriasis.

Key words：Psoriasis;Susceptible genes;Protein interaction network;Enrichment analysis

在自然人群中銀屑?。╬soriasis）的發(fā)病率為0.1%～3.0%，是一種與環(huán)境因素相關(guān)的、遺傳性免疫性皮膚病[1]。隨著全基因組關(guān)聯(lián)分析方法（genome-wide association study，GWAS）的發(fā)展，越來越多的銀屑病易感基因被發(fā)現(xiàn)，OMIM數(shù)據(jù)庫收錄了14個(gè)銀屑病易感位點(diǎn)：PSORS1-PSORS14。對(duì)于這些位點(diǎn)的研究有些已達(dá)到了蛋白質(zhì)水平，確定了易感基因，有些位點(diǎn)還沒有確定易感基因。因大部分候選位點(diǎn)中的基因相互連鎖，很難確定真正的易感基因，而且大部分作用機(jī)制還不明確[2，3]。我們進(jìn)一步分析易感基因的功能，找出關(guān)鍵易感基因。

1材料與方法

1.1候選易感基因? 銀屑病14個(gè)易感位點(diǎn)收錄在OMIM數(shù)據(jù)庫：PSORS1-PSORS14，查詢文獻(xiàn)中的相關(guān)研究內(nèi)容，收集銀屑病候選易感基因如下[3]。

PSORS1（6p21.3）： HLA-C、HCG9、C6orf10、HLA-A、HLA-B、MICA、CDSN、HCR。

PSORS2（17q25）：CMRF35、RUNX1和RAPTOR 基因的第三內(nèi)含子，CARD14 基因等。

PSORS3（4q）：IRF2、IL-21。

PSORS4（1q21上）：LCE3B、LCE3C、FLG。

PSORS5（3q21上）：CSTA、SLC12A8。

PSORS6（19P13）：JUNB、AP-1轉(zhuǎn)錄因子、MUC16、BSG。

PSORS7（1p）：IL-23R、IL-12B。

PSORS8（16q）：CARD15。

PSORS9（4q31）：MGST2、IL15、VEGFC、TLR2、TLR3。

PSORS10（18p11）：候選基因目前未見相關(guān)報(bào)道。

PSORS11（5q31.1）：IL12B 及IL12R。

PSORS12（20q13）：RNF114/ZNF313。

PSORS13（6q21）：TRAF3IP2/CIKS/ACT1。

PSORS14（2q13）：IL-36RN。

1.2富集分析? 利用Gather（http：//changlab.uth.tmc.edu/gather/gather.py）在線工具，對(duì)銀屑病易感基因進(jìn)行功能和通路分析。

1.3蛋白質(zhì)相互作用網(wǎng)絡(luò)? 銀屑病易感基因所表達(dá)的產(chǎn)物所構(gòu)成的蛋白質(zhì)相互作用網(wǎng)絡(luò)，由String（https：//string-db.org/）數(shù)據(jù)庫的在線工具構(gòu)建，這些相互作用即包括了生理上的直接作用也包括了功能上的間接作用。然后把網(wǎng)絡(luò)數(shù)據(jù)導(dǎo)入Cytoscape軟件，計(jì)算節(jié)點(diǎn)度數(shù)，選擇節(jié)點(diǎn)度數(shù)高的基因作為關(guān)鍵的易感基因。

2結(jié)果

2.1 GO分類結(jié)果? 對(duì)收集到的銀屑病易感基因進(jìn)行GO分類，共有97個(gè)類別（Gather在線工具），其中P<0.0001的共有12個(gè)類。在P<0.0001的類別中，有13個(gè)基因（CCRL2，CD300C，HLA-A，HLA-B，HLA-C，IL12B，IL15，IL17C，IRF2，MGST2，MICA，TLR2，TLR3）富集在以下分類：免疫應(yīng)答，防御反應(yīng)，生物應(yīng)激反應(yīng)，有機(jī)生理過程，應(yīng)激反應(yīng)。有8個(gè)基因（CCRL2，CD300C，IL12B，IL17C，MGST2，MICA，TLR2，TLR3）富集在以下分類：對(duì)害蟲寄生蟲病原體的反應(yīng)，對(duì)外部生物刺激的反應(yīng)。有6個(gè)基因（CD300C，IL12B，IL17C，MICA，TLR2，TLR3）富集在對(duì)外傷的反應(yīng)。有4個(gè)基因（HLA-A，HLA-B，HLA-C，MICA）富集在抗原表達(dá)。有3個(gè)基因（HLA-A，HLA-B，HLA-C）富集在抗原呈遞及內(nèi)源性抗原，抗原加工及通過MHCⅠ類的內(nèi)源性抗原。

2.2 KEGG通路分析? KEGG通路富集分析提示，本文所研究的易感基因顯著富集在四條代謝通路：Toll樣受體信號(hào)通路，細(xì)胞因子-細(xì)胞因子受體相互作用，JAK-STAT信號(hào)轉(zhuǎn)導(dǎo)途徑，谷胱甘肽代謝。

2.3蛋白質(zhì)相互作用網(wǎng)絡(luò)? 由String數(shù)據(jù)庫得到銀屑病易感基因的蛋白質(zhì)相互作用網(wǎng)絡(luò)，見圖1。該網(wǎng)絡(luò)有34個(gè)節(jié)點(diǎn)對(duì)應(yīng)34個(gè)蛋白質(zhì)，有42條邊對(duì)應(yīng)蛋白質(zhì)的相互作用，節(jié)點(diǎn)的平均度數(shù)為2.47。與一個(gè)基因相互作用的基因數(shù)目越多節(jié)點(diǎn)的度數(shù)越高，下載string相互作用網(wǎng)絡(luò)數(shù)據(jù)，計(jì)算出節(jié)點(diǎn)的度數(shù)，選擇度數(shù)高的基因作為關(guān)鍵易感基因。按照度數(shù)由高到低排列排在前十的分別是IL15（9），HLA-C（9），IL12B（7），NOD2（6），CARD14（6），HLA-A（5），CDSN（4），HLA-B（4），IRF2（4），IL23R（4）。

3討論

銀屑病目前雖不能完全確定病因，但本病是與環(huán)境因素相關(guān)的遺傳性免疫皮膚病[1]，受潮、感染、吸煙、飲酒、藥物等都影響著本病的發(fā)病和病程[2]。我們對(duì)易感基因的功能富集分析發(fā)現(xiàn)，銀屑病的易感基因功能主要富集在免疫應(yīng)答，防御反應(yīng)，生物應(yīng)激反應(yīng)，對(duì)害蟲寄生蟲病原體的反應(yīng)，對(duì)外傷外部生物刺激的反應(yīng)，有機(jī)生理過程，以及抗原和抗原呈遞等類別。功能富集在這些類的易感基因?qū)е裸y屑病是和吸煙、嗜酒、受潮、感染和藥物等環(huán)境因素密切相關(guān)。和免疫應(yīng)答相關(guān)的基因表明銀屑病是具有遺傳性的免疫性皮膚病。我們對(duì)基因功能的分析和目前對(duì)銀屑病發(fā)病原因的理解一致，是對(duì)發(fā)病原因的一次理論驗(yàn)證。

我們根據(jù)銀屑病易感基因所編碼的蛋白質(zhì)的相互作用網(wǎng)絡(luò)中節(jié)點(diǎn)的度數(shù)高低，篩選出10個(gè)關(guān)鍵基因。度數(shù)最高的兩個(gè)基因是與漢族人群銀屑病顯著相關(guān)的基因IL-15位于PSORS9（4q31），和位于最主要的易感位點(diǎn)PSORS1（6p21.3）的基因HLA-C，經(jīng)實(shí)驗(yàn)證實(shí)該基因是歐美及其他人群銀屑病易感位點(diǎn)[4]，也是銀屑病性關(guān)節(jié)炎的易感位點(diǎn)[5]，并經(jīng)過Meta分析證實(shí)[6]。在波蘭患者中HLA-Cw*06基因頻率顯著升高，證實(shí)與早發(fā)型銀屑病顯著相關(guān)[7]。由于HLA-C基因沒有功能性作用，所以不是易感基因本身只一個(gè)遺傳標(biāo)記[8]。根據(jù)富集分析可知，這兩個(gè)基因富集在免疫應(yīng)答，防御反應(yīng)，生物應(yīng)激反應(yīng)，有機(jī)生理過程，應(yīng)激反應(yīng)，可見這兩個(gè)關(guān)鍵基因也佐證了銀屑病是和環(huán)境因素相關(guān)的免疫性皮膚病這一判斷。HLA-C這個(gè)基因是關(guān)鍵的銀屑病易位點(diǎn)，但卻缺少明確功能。我們所篩選出的關(guān)鍵易感基因，也是許多實(shí)驗(yàn)驗(yàn)證的關(guān)鍵易感基因，這是對(duì)實(shí)驗(yàn)結(jié)果的一次理論驗(yàn)證。實(shí)驗(yàn)證明該關(guān)鍵基因無具體功能，只是易感位點(diǎn)，通過互作網(wǎng)絡(luò)發(fā)現(xiàn)，該基因和很多功能明確與銀屑病相關(guān)的基因相互作用，所以該基因雖無具體功能，但理論和實(shí)驗(yàn)都提示該基因是關(guān)鍵易感基因。通過互作網(wǎng)絡(luò)分析，解釋了關(guān)鍵易感基因缺少實(shí)際功能的現(xiàn)象，生物信息分析提供了一個(gè)新的視角。

參考文獻(xiàn)：

[1]Nickoloff BJ.The immunologic and genetic basis of psoriasis[J].Arch Dermatol，1999，135（9）：1104-1110.

[2]張亞南，吳閩楓，金性秀，等.泛發(fā)性膿皰型銀屑病易感基因研究進(jìn)展[J].中國皮膚性病學(xué)雜志，2016，27（7）：751-752.

[3]趙恬，羅權(quán)，林玲，等.銀屑病易感基因的研究進(jìn)展[J].中國中西醫(yī)結(jié)合皮膚性病學(xué)雜志，2015，14（5）：333-338.

[4]Ellinghaus E，Ellinghaus D，Stuart PE，et al.Genome-wide association study identifies a psoriasis susceptibility locus at TRAF3IP2[J].Nat Genet，2010，42（11）：991-995.

[5]Hüffmeier U，Uebe S，Ekici AB，et al.Common variants at TRAF3IP2 are associated with susceptibility to psoriatic arthritis and psoriasis[J].Nat Genet，2010，42（11）：996-999.

[6]Ellinghaus E，Stuart PE，Ellinghaus D，et al.Genome-wide metaanalysis of psoriatic arthritis identifies susceptibility locus at REL[J].J Invest Dermatol，2011，132（4）：1133-1140.

[7]Rebala K，Szczerkowska-Dobosz A，Niespodziana K，et al.Simple and rapid screening for HLA-CW*06 in Polish patients with psoriasis[J].Clinical EXP Dermatol，2010，35（4）：431-436.

[8]Jordan CT，Cao L，Roberson ED，et al.PSORS2 Is Due to Mutations in CARD14[J].Am J Hum Genet，2012，90（4）：784-795.