梁釗 林栩 覃卿 杜秀日 唐志明 王晨 韋美理

【摘要】目的基于基因表達(dá)綜合數(shù)據(jù)庫(kù)（GEO）挖掘與狼瘡性腎炎（LN）相關(guān)的潛在基因。

方法從GEO中搜集LN相關(guān)的樣本，獲得GSE32591、GSE81622、GSE99967共3個(gè)數(shù)據(jù)集。利用GEO2R平臺(tái)對(duì)這3個(gè)數(shù)據(jù)集進(jìn)行分析，篩選出共同差異基因，并利用在線分析工具DAVID完成GO富集分析和KEGG通路富集分析。將篩選出的共同差異基因?qū)隨TRING在線數(shù)據(jù)庫(kù)，構(gòu)建共同差異基因的蛋白蛋白相互作用網(wǎng)絡(luò)，利用Cytoscape軟件進(jìn)行模塊分析并識(shí)別樞紐基因。

結(jié)果篩選得到110個(gè)共同差異基因。GO富集分析發(fā)現(xiàn)這些基因主要參與了防御反應(yīng)、免疫反應(yīng)、對(duì)病毒的反應(yīng)、對(duì)細(xì)菌的反應(yīng)、對(duì)傷害的反應(yīng)、炎癥反應(yīng)等生物學(xué)過程。KEGG通路富集分析主要包括了系統(tǒng)性紅斑狼瘡、抗原處理和呈遞、補(bǔ)體與凝血級(jí)聯(lián)、RIGI樣受體等信號(hào)通路。從最顯著基因模塊中識(shí)別出10個(gè)樞紐基因：RSAD2、OAS1、MX1、ISG15、DDX58、IFIH1、IFI44、IFI44L、IFIT1、IFIT3。

結(jié)論生物信息學(xué)分析顯示RSAD2、OAS1、MX1、ISG15、DDX58、IFIH1、IFI44、IFI44L、IFIT1、IFIT3可能是與LN相關(guān)的樞紐基因，為L(zhǎng)N的機(jī)制研究提供了一種全新的思路。

【關(guān)鍵詞】狼瘡性腎炎;生物信息學(xué);基因;共表達(dá)網(wǎng)絡(luò)

中圖分類號(hào)：R593.24+2文獻(xiàn)標(biāo)志碼：ADOI：10.3969/j.issn.10031383.2020.06.002

【Abstract】ObjectiveTo explore potential? the genes associated with lupus nephritis（LN） based on Gene Expression Omnibus（GEO）.

MethodsLN related samples were collected from the GEO database，and three data sets（GSE32591，GSE81622 and GSE99967） were obtained.GEO2R platform was applied to analyze these three data sets，and common differential genes were screened.In addition，online analysis tool DAVID was used to complete GO enrichment analysis and KEGG pathway enrichment analysis.The common differential genes were imported into the STRING online database，and the proteinprotein interaction network of the common differential genes was constructed.The module analysis and the pivot genes were identified by the Cytoscape software.

Results110 common differential genes were identified.The GO functional enrichment analysis showed that these genes were mainly involved in the biological processes of defense response，immune response，response to virus，response to bacteria，response to injury，and inflammatory response，etc.The KEGG signal pathway enrichment analysis mainly included systemic lupus erythematosus，antigen processing and presentation，complement and coagulation cascades，RIGIlike receptor signaling pathway，etc.Ten pivot genes were identified from the most significant modules：RSAD2，OAS1，MX1，ISG15，DDX58，IFIH1，IFI44，IFI44L，IFIT1 and IFIT3.

ConclusionBioinformatics analysis shows that RSAD2，OAS1，MX1，ISG15，DDX58，IFIH1，IFI44，IFI44L，IFIT1 and IFIT3? may be the pivot genes related to LN，which provides a new idea for the pathogenesis research of LN.

【Key words】LN;bioinformatics;gene;coexpression network

狼瘡性腎炎（lupus nephritis，LN）是系統(tǒng)性紅斑狼瘡（systemic lupus erythematosus，SLE）累及腎臟所引起的一種免疫復(fù)合物性腎炎，該病臨床表現(xiàn)多樣，主要包括血尿、蛋白尿、輕度腎功能損傷、急進(jìn)性腎功能衰竭等，其發(fā)病機(jī)制錯(cuò)綜復(fù)雜[1]。隨著基因芯片、RNAseq等高通量測(cè)序技術(shù)的發(fā)展，利用生物信息學(xué)方法挖掘公共數(shù)據(jù)庫(kù)中與LN相關(guān)的基因，對(duì)LN的機(jī)制研究有著重要的意義。美國(guó)國(guó)立生物技術(shù)信息中心（national center for biotechnology information，NCBI）的基因表達(dá)綜合數(shù)據(jù)庫(kù)（gene expression omnibus，GEO）是一個(gè)公共基因表達(dá)數(shù)據(jù)庫(kù)，收錄和儲(chǔ)存了大量的高通量實(shí)驗(yàn)數(shù)據(jù)，研究人員可以免費(fèi)下載這些數(shù)據(jù)進(jìn)行分析研究[2]。本研究擬通過挖掘GEO數(shù)據(jù)庫(kù)，篩選在LN中有差異表達(dá)的基因，構(gòu)建共表達(dá)基因的蛋白互作網(wǎng)絡(luò)并識(shí)別與LN相關(guān)的樞紐基因，為L(zhǎng)N的機(jī)制研究提供一種全新的思路。

1對(duì)象與方法

1.1研究對(duì)象

對(duì)公共的基因芯片數(shù)據(jù)庫(kù)GEO（https：//www.ncbi.nlm.nih.gov/geo/）進(jìn)行LN相關(guān)樣本檢索，得到GEO數(shù)據(jù)庫(kù)中GSE32591、GSE81622、GSE99967共3個(gè)LN人群研究的數(shù)據(jù)集。GSE32591為來自LN患者腎小球組織和正常對(duì)照組、腎間質(zhì)組織和正常對(duì)照組進(jìn)行對(duì)比的表達(dá)譜芯片，GSE81622為將LN患者與正常對(duì)照組的外周血單個(gè)核細(xì)胞（peripheral blood mononuclear cell，PBMC）進(jìn)行比較的表達(dá)譜芯片，GSE99967為L(zhǎng)N患者與正常對(duì)照組的外周血進(jìn)行比較的表達(dá)譜芯片。

1.2研究方法

1.2.1差異基因選取

GEO2R程序是NCBI中基于R語(yǔ)言的在線分析程序，我們使用GEO2R平臺(tái)（https：//www.ncbi.nlm.nih.gov/geo/geo2r/）將樣本分組并篩選差異基因，按平臺(tái)默認(rèn)方案選出每個(gè)芯片中表達(dá)差異明顯的基因，設(shè)定P<0.05、差異倍數(shù)絕對(duì)值≥1（|logFC|≥1）作為篩選差異基因的閾值，符合條件的差異基因數(shù)據(jù)認(rèn)為其表達(dá)具有顯著性差異。

1.2.2共同差異基因選取

根據(jù)3個(gè)數(shù)據(jù)集分組的不同，將GSE32591中LN患者腎小球與正常腎小球?qū)φ盏男酒瑪?shù)據(jù)設(shè)為“腎小球組”，將GSE32591中LN患者腎間質(zhì)與正常腎間質(zhì)對(duì)照的芯片數(shù)據(jù)設(shè)為“腎間質(zhì)組”，將GSE81622的芯片數(shù)據(jù)設(shè)為“PBMC組”，將GSE99967設(shè)為“外周血組”。利用在線分析工具（https：//bioinfogp.cnb.csic.es/tools/venny/index.html）篩選出四組數(shù)據(jù)中的共同差異基因。

1.2.3共同差異基因的GO和KEGG通路分析

對(duì)共同差異基因使用在線工具DAVID 6.7（https：//davidd.ncifcrf.gov/）進(jìn)行基因本體（gene ontology，GO）富集分析和京都基因與基因組百科全書通路（kyoto encyclopedia of genes and genomes pathway，KEGG）富集分析。

1.2.4蛋白蛋白相互作用網(wǎng)絡(luò)的構(gòu)建（proteinprotein interaction network，PPI網(wǎng)絡(luò)） 及關(guān)鍵基因的識(shí)別

將篩選出的共同差異基因?qū)朐诰€數(shù)據(jù)庫(kù)STRING（https：//stringdb.org/）構(gòu)建PPI網(wǎng)絡(luò)并進(jìn)行聚類分析，綜合得分>0.4認(rèn)為蛋白間的相互作用顯著。隨后，用生物信息學(xué)繪圖軟件Cytoscape的MCODE插件篩選PPI網(wǎng)絡(luò)中關(guān)鍵節(jié)點(diǎn)基因聚類形成的顯著基因模塊。之后，用Cytoscape的cytoHubba插件從顯著基因模塊中識(shí)別與LN相關(guān)的樞紐基因。

2結(jié)果

2.1差異基因和共同差異基因的選取結(jié)果

腎小球組共篩選出361個(gè)差異表達(dá)基因，腎間質(zhì)組共篩選出133個(gè)差異表達(dá)基因，PBMC組共篩選出158個(gè)差異表達(dá)基因，外周血組共篩選出122個(gè)差異表達(dá)基因。用韋恩圖篩選出這4組數(shù)據(jù)的共同表達(dá)差異基因，共110個(gè)。見圖1。

2.2GO富集分析和KEGG通路富集分析

GO富集分析包括生物學(xué)功能、細(xì)胞組分和分子功能三個(gè)方面。將110個(gè)共同表達(dá)差異基因?qū)隓AVID數(shù)據(jù)庫(kù)進(jìn)行GO富集分析和KEGG通路富集分析，發(fā)現(xiàn)富集到的生物學(xué)過程主要包括防御反應(yīng)、免疫反應(yīng)、對(duì)病毒的反應(yīng)、對(duì)細(xì)菌的反應(yīng)、對(duì)傷害的反應(yīng)、炎癥反應(yīng)等過程;參與的細(xì)胞組分包括細(xì)胞外間隙、細(xì)胞外區(qū)域、細(xì)胞外區(qū)域部分、血紅蛋白復(fù)合體、補(bǔ)體成分C1復(fù)合物等;發(fā)揮了血紅蛋白結(jié)合、核糖核酸內(nèi)切酶活性、核酸內(nèi)切酶活性、氧結(jié)合、抗原結(jié)合等分子功能;KEGG通路主要包括了系統(tǒng)性紅斑狼瘡、抗原處理和呈遞、補(bǔ)體與凝血級(jí)聯(lián)、RIGI樣受體等信號(hào)通路。見圖2。

2.3PPI網(wǎng)絡(luò)

共同差異基因編碼蛋白之間的相互作用網(wǎng)絡(luò)及聚類分析，發(fā)現(xiàn)有5個(gè)不同類別的聚類，其中不同顏色代表蛋白的不同類別。見圖3。

將PPI網(wǎng)絡(luò)導(dǎo)入Cytoscape，使用MCODE插件篩選顯著基因模塊，用cytoHubba插件從最顯著基因模塊中篩選出10個(gè)hub基因：RSAD2、OAS1、MX1、ISG15、DDX58、IFIH1、IFI44、IFI44L、IFIT1、IFIT3。見圖4。

3討論

LN是SLE累及腎臟的自身免疫性疾病，遺傳、環(huán)境、病毒感染、性激素等多種因素均可導(dǎo)致，其發(fā)病機(jī)制尚未完全明確[3]。GEO數(shù)據(jù)庫(kù)收錄了大量與疾病相關(guān)的高通量數(shù)據(jù)，對(duì)LN的機(jī)制研究具有十分重要的意義。本研究納入GEO數(shù)據(jù)庫(kù)中LN患者腎組織、PBMC和外周血的基因表達(dá)數(shù)據(jù)，篩選出共同差異基因110個(gè)。

GO富集分析發(fā)現(xiàn)，這些共同差異基因在防御反應(yīng)、免疫反應(yīng)、對(duì)病毒的反應(yīng)、對(duì)細(xì)菌的反應(yīng)、對(duì)傷害的反應(yīng)、炎癥反應(yīng)、類固醇激素刺激的反應(yīng)、對(duì)糖皮質(zhì)激素刺激的反應(yīng)、對(duì)細(xì)胞因子刺激的反應(yīng)和固有免疫應(yīng)答等生物學(xué)過程發(fā)揮作用。有研究顯示，LN的發(fā)病是免疫復(fù)合物沉積于腎組織引起腎臟損傷的過程，這個(gè)過程包括炎癥細(xì)胞募集、細(xì)胞因子產(chǎn)生、氧化損傷、補(bǔ)體吸收和成纖維細(xì)胞活化增殖等，其中，炎癥和纖維化是LN發(fā)生發(fā)展的關(guān)鍵過程，它們涉及先天性免疫和適應(yīng)性免疫的免疫細(xì)胞與常駐腎細(xì)胞的相互作用[4～5]。SLE的發(fā)病可能與多種病毒感染相關(guān)，如EB病毒、人類皰疹病毒、巨細(xì)胞病毒、水痘帶狀皰疹病毒等[6]。YU等[7]研究發(fā)現(xiàn)，在LN患者的腎組織標(biāo)本中EB病毒陽(yáng)性率顯著高于對(duì)照組腎組織標(biāo)本，且腎組織EB病毒標(biāo)記陽(yáng)性組抗Sm陽(yáng)性率也高于腎組織EB病毒標(biāo)記陰性組，這些發(fā)現(xiàn)提示抗Sm的產(chǎn)生可能與EB病毒感染有關(guān)，腎EB病毒感染可能通過誘導(dǎo)產(chǎn)生抗Sm抗體而參與LN的發(fā)病機(jī)制。

KEGG通路富集分析顯示，這些共同差異基因主要富集在系統(tǒng)性紅斑狼瘡、抗原處理和呈遞、補(bǔ)體與凝血級(jí)聯(lián)、RIGI樣受體等信號(hào)通路。當(dāng)SLE累及腎臟時(shí)，免疫復(fù)合物沉積于腎小球內(nèi)皮細(xì)胞和系膜區(qū)引起炎癥細(xì)胞募集、補(bǔ)體激活等過程，最終導(dǎo)致不同程度的腎臟損傷[4]。有研究顯示，RIGI在LN患者的腎臟組織標(biāo)本中有顯著的表達(dá)[8]，此外，還有研究發(fā)現(xiàn)LN患者的尿沉渣中RIGI mRNA的表達(dá)水平高于IgA腎病患者[9]，這些結(jié)果表明RIGI可能參與了LN的發(fā)生過程。

通過構(gòu)建PPI網(wǎng)絡(luò)篩選顯著基因模塊，我們從最顯著基因模塊中識(shí)別出了10個(gè)可能與LN發(fā)生發(fā)展高度相關(guān)的樞紐基因：RSAD2、OAS1、MX1、ISG15、DDX58、IFIH1、IFI44、IFI44L、IFIT1、IFIT3，這些基因均參與Ⅰ型干擾素信號(hào)通路。環(huán)境因素與SLE的發(fā)病密切相關(guān)，其中EB病毒感染被認(rèn)為是SLE的環(huán)境危險(xiǎn)因素之一。在感染EB病毒后，機(jī)體通過多種途徑誘導(dǎo)Ⅰ型干擾素信號(hào)，并激活免疫系統(tǒng)產(chǎn)生炎癥介質(zhì)和自身抗體參與SLE的發(fā)病過程，自身抗體沉積于腎臟引起腎臟損害便可導(dǎo)致LN[10～13]。最近的研究表明，與正常人相比，SLE患者的Ⅰ型干擾素和γ干擾素表達(dá)水平均明顯升高，且Ⅰ型干擾素和γ干擾素的激活密切相關(guān)[14～15]。DER等[16]研究表明，與健康對(duì)照組相比，LN患者腎小管細(xì)胞中Ⅰ型干擾素的表達(dá)水平顯著升高，同時(shí)Ⅰ型干擾素在LN患者的腎小球系膜細(xì)胞和內(nèi)皮細(xì)胞中也有不同程度的升高。這些研究說明，Ⅰ型干擾素與LN的發(fā)生發(fā)展密不可分，對(duì)Ⅰ型干擾素在LN中的發(fā)病機(jī)制進(jìn)行深入研究有重要的意義。

RSAD2是干擾素誘導(dǎo)的鐵硫簇結(jié)合抗病毒蛋白，有研究發(fā)現(xiàn)，在SLE患者幼稚CD4+T細(xì)胞中甲基化的RSAD2表達(dá)明顯升高，RSAD2基因CpG位點(diǎn)甲基化與SLE相關(guān)自身抗體的產(chǎn)生相關(guān)[17～18]。JOSEPH等[19]最近的研究證實(shí)，在SLE患者中CpG位點(diǎn)甲基化的RSAD2基因與SLE疾病活動(dòng)指數(shù)相關(guān)。而在一項(xiàng)LN全基因組DNA甲基化的研究中，RSAD2基因CpG位點(diǎn)甲基化被證明與LN有很高的相關(guān)性[20]。

OAS1是2'5'寡腺苷酸合成酶（OAS）的一種亞型，已知其功能是抗病毒。有研究表明，與對(duì)照組相比，SLE患者外周血中的OAS1和OAS1 mRNA表達(dá)水平顯著升高[21～22]。LANDOLTMARTICORENA等[23]研究顯示，B細(xì)胞活化因子和OAS1的表達(dá)水平在SLE中均升高，且兩者的表達(dá)水平緊密相關(guān)，與對(duì)照組相比，SLE患者外周血中的B細(xì)胞增生異常，表明OAS1的高表達(dá)可能與B細(xì)胞的不同表型相關(guān)。FENG等[24]對(duì)具有不同疾病活動(dòng)指數(shù)的SLE患者進(jìn)行了研究，并根據(jù)病情的變化情況對(duì)患者進(jìn)行3～12個(gè)月的隨訪，發(fā)現(xiàn)與正常對(duì)照組相比，SLE患者外周血中OAS1的表達(dá)水平顯著升高，且OAS1的高表達(dá)與SLE疾病活動(dòng)性、LN均密切相關(guān)。

MX1是一種干擾素誘導(dǎo)基因。有學(xué)者發(fā)現(xiàn)，在LN患者腎小球系膜區(qū)中MX1染色呈顯著陽(yáng)性，而在IgA腎病、紫癜性腎炎、直立性蛋白尿患者的腎組織中未能檢測(cè)到MX1的表達(dá)[25]。SHIMIZU等[26]研究發(fā)現(xiàn)，在LN患者PBMC中MX1蛋白的表達(dá)顯著高于IgA腎病患者、ANCA相關(guān)性血管炎患者和健康對(duì)照組，在腎組織活檢中發(fā)現(xiàn)，LN患者的腎小球和腎小管區(qū)域均可檢測(cè)到MX1，而在IgA腎病和ANCA相關(guān)性血管炎患者中僅觀察到非常微弱的MX1染色。有研究表明miR155可能通過靶向STAT1信號(hào)通路調(diào)節(jié)Ⅰ型干擾素的應(yīng)答[27]，而LEISS等[28]在使用miR155敲除小鼠研究SLE發(fā)病機(jī)制時(shí)發(fā)現(xiàn)，與誘導(dǎo)的野生型狼瘡小鼠相比，miR155缺陷型狼瘡小鼠自身抗體的血清水平和腎臟損傷程度明顯減少，MX1的表達(dá)也減少。這些發(fā)現(xiàn)提示MX1可能是LN發(fā)生發(fā)展中的一個(gè)重要的調(diào)節(jié)因子。

ISG15也是一種干擾素誘導(dǎo)基因。YUAN等[29]研究顯示，SLE患者血液中ISG15的表達(dá)與正常對(duì)照組相比明顯升高，ISG15的表達(dá)水平與SLE疾病活動(dòng)相關(guān)，同時(shí)該研究還發(fā)現(xiàn)SLE患者體內(nèi)淋巴細(xì)胞的減少與ISG15的表達(dá)水平密切相關(guān)，其機(jī)制可能是ISG15參與淋巴細(xì)胞凋亡的過程。HAN等[14]在研究SLE患者EB病毒感染和ISG15的關(guān)系時(shí)，發(fā)現(xiàn)ISG15與EB病毒基因LMP1的表達(dá)呈正相關(guān)，表明可能是EB病毒的LMP1基因通過激活Ⅰ型干擾素通路參與SLE的發(fā)病機(jī)制。此前的研究報(bào)道顯示，使用姥鮫烷（Pristane）誘導(dǎo)的狼瘡小鼠ISG15的表達(dá)顯著上調(diào)，而ISG15在用Pristane誘導(dǎo)的miR155敲除小鼠中表達(dá)減少，miR155缺陷型狼瘡小鼠的自身抗體血清水平和腎臟損傷程度也明顯減少，提示ISG15與LN的腎臟損傷密切相關(guān)[28]。

DDX58是RIGI的編碼基因。RIGI是誘導(dǎo)Ⅰ型干擾素信號(hào)的最常見受體，作為病毒傳感器，RIGI可以檢測(cè)到病毒的入侵并觸發(fā)抗病毒應(yīng)答，其持續(xù)異常的激活會(huì)引起自身免疫性疾病[30～31]。SUZUKI等[8]在活動(dòng)性LN患者的腎臟標(biāo)本中發(fā)現(xiàn)RIGI有顯著的表達(dá)，且RIGI的表達(dá)水平與腎小球病變的嚴(yán)重程度相關(guān)。之前有學(xué)者發(fā)現(xiàn)，地塞米松能夠抑制星形膠質(zhì)細(xì)胞的RIGI表達(dá)水平[32]，IMAIZUMI等[33]的研究發(fā)現(xiàn)，γ干擾素能夠誘導(dǎo)腎小球系膜細(xì)胞RIGI的表達(dá)，而使用地塞米松干預(yù)不能夠影響γ干擾素誘導(dǎo)的RIGI的表達(dá)，表明RIGI可能與腎小球的炎癥相關(guān)，腎小球系膜細(xì)胞產(chǎn)生的RIGI可能參與了LN的發(fā)病機(jī)制。

IFIH1是一種干擾素誘導(dǎo)基因，可以編碼黑色素瘤分化相關(guān)蛋白（MDA5）。既往的研究顯示，MDA5參與病毒感染后干擾素應(yīng)答調(diào)節(jié)的過程，且SLE的疾病活性與MDA5的水平呈負(fù)相關(guān)[34]。FUNABIKI等[35]的研究顯示IFIH1 G821S錯(cuò)義突變的小鼠可出現(xiàn)狼瘡樣癥狀，該研究在小鼠血清和腎臟中分別檢測(cè)到抗核抗體、抗dsDNA抗體的表達(dá)和免疫球蛋白、補(bǔ)體的沉積，同時(shí)，腎臟中干擾素、IL6、IP10等炎性介質(zhì)的表達(dá)均顯著上調(diào)，說明錯(cuò)義突變的IFIH1可能通過激活干擾素信號(hào)觸發(fā)機(jī)體的自身免疫反應(yīng)。為了驗(yàn)證FUNABIKI等人的研究結(jié)果，MUNROE等[36]對(duì)SLE患者血液中IFIH1與炎癥介質(zhì)、自身抗體的相關(guān)性進(jìn)行評(píng)估，發(fā)現(xiàn)IFIH1基因與炎性介質(zhì)IL6、IP10、自身抗體的表達(dá)密切相關(guān)。這些研究表明，IFIH1可能與LN的腎臟損傷有關(guān)。

IFI44、IFI44L屬于干擾素誘導(dǎo)蛋白家族，是干擾素誘導(dǎo)基因，參與了Ⅰ型干擾素信號(hào)通路。KIROU等[37]研究顯示，在SLE的PBMC中IFI44的表達(dá)顯著增高，提示IFI44參與了α干擾素信號(hào)通路介導(dǎo)的SLE發(fā)生發(fā)展過程，但他們的研究并未揭示IFI44與SLE腎臟損傷的關(guān)系。ZHAO等[38]研究發(fā)現(xiàn)，SLE患者IFI44L啟動(dòng)子DNA甲基化水平顯著低于正常對(duì)照組，同時(shí)，有腎臟受累的SLE患者IFI44L啟動(dòng)子DNA甲基化水平也顯著低于沒有腎臟受累的SLE患者，提示IF144L啟動(dòng)子的DNA甲基化水平可能與SLE患者的腎臟損害有關(guān)。

IFIT1、IFIT3屬于干擾素誘導(dǎo)的四肽重復(fù)蛋白家族，它們同樣是干擾素誘導(dǎo)基因。有研究發(fā)現(xiàn)，IFIT1在SLE患者的外周血中表達(dá)升高，而在LN患者的腎臟組織中同樣也觀察到IFIT1的表達(dá)升高[39]。HU等[40]使用MRL/lpr狼瘡小鼠模型研究IFIT1的表達(dá)在LN腎臟病理變化中的作用，發(fā)現(xiàn)IFIT1與MRL/lpr小鼠腎組織中Factin、Nephrin和Podocin三種足細(xì)胞蛋白的表達(dá)呈負(fù)相關(guān)，進(jìn)一步的病理分析提示隨著IFIT1表達(dá)的升高，MRL/lpr小鼠的腎組織中足細(xì)胞大量丟失，說明上調(diào)IFIT1的表達(dá)可能通過破壞足細(xì)胞的細(xì)胞結(jié)構(gòu)引起腎臟的損傷。WANG等[41]的研究結(jié)果表明，與正常對(duì)照組相比，IFIT3在SLE患者PBMC表達(dá)明顯升高，且IFIT3可以激活cGASsting信號(hào)通路。XIAO等[42]利用CRISPR/Cas9技術(shù)，構(gòu)建了可產(chǎn)生自身免疫癥狀的基因突變Trex1D18N/D18N小鼠，當(dāng)刪除Trex1D18N/D18N小鼠cGAS的一個(gè)等位基因時(shí)發(fā)現(xiàn)腎臟的炎癥顯著減輕，進(jìn)一步的研究發(fā)現(xiàn)，是cGAS介導(dǎo)了Trex1D18N/D18N自身抗體的產(chǎn)生，而自身抗體的產(chǎn)生是狼瘡樣自身免疫的標(biāo)志之一。這些發(fā)現(xiàn)表明IFIT3和cGAS信號(hào)通路可能參與了狼瘡樣Trex1D18N/D18N小鼠的腎臟損害。

綜上所述，本研究運(yùn)用生物信息學(xué)的方法挖掘了與LN相關(guān)的樞紐基因，發(fā)現(xiàn)這些樞紐基因在LN的發(fā)生發(fā)展中起到關(guān)鍵的作用，這為L(zhǎng)N的機(jī)制研究提供了一種全新的思路。本研究的不足之處在于，我們未能通過實(shí)驗(yàn)進(jìn)一步驗(yàn)證這些樞紐基因的表達(dá)水平。因此，未來我們可以通過后續(xù)的實(shí)驗(yàn)研究，驗(yàn)證這些樞紐基因的功能，這可能是未來研究的熱點(diǎn)。

聲明：本研究不存在任何利益沖突。

參考文獻(xiàn)

[1] KRONBICHLER A，BREZINA B，GAUCKLER P，et al.Refractory lupus nephritis：When，why and how to treat[J].Autoimmun Rev，2019，18（5）：510518.

[2] CLOUGH E，BARRETT T.The gene expression omnibus database[J].Methods Mol Biol，2016，1418：93110.

[3] ZANDMANGODDARD G，SOLOMON M，ROSMAN Z，et al.Environment and lupusrelated diseases[J].Lupus，2012，21（3）：241250.

[4] KOMOLAFE O O.Rapidly progressive glomerulonephritis：a wild card manifestation of lupus nephritis[J].Saudi J Kidney Dis Transpl，2018，29（2）：443451.

[5] YUNG S，YAP D Y，CHAN T M.Recent advances in the understanding of renal inflammation and fibrosis in lupus nephritis[J].F1000Res，2017，6：874.

[6] DRABORG A，IZARZUGAZA J M，HOUEN G.How compelling are the data for EpsteinBarr virus being a trigger for systemic lupus and other autoimmune diseases？[J].Curr Opin Rheumatol，2016，28（4）：398404.

[7] YU X X，YAO C W，TAO J L，et al.The expression of renal EpsteinBarr virus markers in patients with lupus nephritis[J].Exp Ther Med，2014，7（5）：11351140.

[8] SUZUKI K，IMAIZUMI T，TSUGAWA K，et al.Expression of retinoic acidinducible geneI in lupus nephritis[J].Nephrol Dial Transplant，2007，22（8）：24072409.

[9] TSUGAWA K，OKI E，SUZUKI K，et al.Expression of mRNA for functional molecules in urinary sediment in glomerulonephritis[J].Pediatr Nephrol，2008，23（3）：395401.

[10] HARLEY J B，HARLEY I T，GUTHRIDGE J M，et al.The curiously suspicious：a role for EpsteinBarr virus in lupus[J].Lupus，2006，15（11）：768777.

[11] CROW M K，OLFERIEV M，KIROU K A.Type I interferons in autoimmune disease[J].Annu Rev PatholMech Dis，2019，14（1）：369393.

[12] RNNBLOM L，ALM G V，ELORANTA M L.The type I interferon system in the development of lupus[J].Semin Immunol，2011，23（2）：113121.

[13] MOHAN C，PUTTERMAN C.Genetics and pathogenesis of systemic lupus erythematosus and lupus nephritis[J].Nat Rev Nephrol，2015，11（6）：329341.

[14] HAN L，ZHANG Y，WANG Q，et al.EpsteinBarr virus infection and type I interferon signature in patients with systemic lupus erythematosus[J].Lupus，2018：961203317753069.

[15] LIU M M，LIU J L，HAO S M，et al.Higher activation of the interferonGamma signaling pathway in systemic lupus erythematosus patients with a high type I IFN score：relation to disease activity[J].Clin Rheumatol，2018，37（10）：26752684.

[16] DER E，SURYAWANSHI H，MOROZOV P，et al.Tubular cell and keratinocyte singlecell transcriptomics applied to lupus nephritis reveal type I IFN and fibrosis relevant pathways[J].Nat Immunol，2019，20（7）：915927.

[17] COIT P，JEFFRIES M，ALTOROK N，et al.Genomewide DNA methylation study suggests epigenetic accessibility and transcriptional poising of interferonregulated genes in naive CD4+T cells from lupus patients[J].J Autoimmun，2013，43：7884.

[18] CHUNG S A，NITITHAM J，ELBOUDWAREJ E，et al.Genomewide assessment of differential DNA methylation associated with autoantibody production in systemic lupus erythematosus[J].PLoS One，2015，10（7）：e0129813.

[19] JOSEPH S，GEORGE N I，GREENKNOX B，et al.Epigenomewide association study of peripheral blood mononuclear cells in systemic lupus erythematosus：Identifying DNA methylation signatures associated with interferonrelated genes based on ethnicity and SLEDAI[J].J Autoimmun，2019，96：147157.

[20] MOK A，SOLOMON O，NAYAK R R，et al.Genomewide profiling identifies associations between lupus nephritis and differential methylation of genes regulating tissue hypoxia and type 1 interferon responses[J].Lupus Sci Med，2016，3（1）：e000183.

[21] YE S，GUO Q，TANG J P，et al.Could 2'5'oligoadenylate synthetase isoforms be biomarkers to differentiate between disease flare and infection in lupus patients？A pilot study[J].Clin Rheumatol，2007，26（2）：186190.

[22] FENG X B，WU H，GROSSMAN J M，et al.Association of increased interferoninducible gene expression with disease activity and lupus nephritis in patients with systemic lupus erythematosus[J].Arthritis Rheum，2006，54（9）：29512962.

[23] LANDOLTMARTICORENA C，WITHER R，REICH H，et al.Increased expression of B cell activation factor supports the abnormal expansion of transitional B cells in systemic lupus erythematosus[J].J Rheumatol，2011，38（4）：642651.

[24] FENG X B，HUANG J，LIU Y，et al.Identification of interferoninducible genes as diagnostic biomarker for systemic lupus erythematosus[J].Clin Rheumatol，2015，34（1）：7179.

[25] WATANABE S，IMAIZUMI T，TSURUGA K，et al.Glomerular expression of myxovirus resistance protein 1 in human mesangial cells：possible activation of innate immunity in the pathogenesis of lupus nephritis[J].Nephrology （Carlton），2013，18（12）：833837.

[26] SHIMIZU Y，YASUDA S，KIMURA T，et al.Interferoninducible Mx1 protein is highly expressed in renal tissues from treatmentnaive lupus nephritis，but not in those under immunosuppressive treatment[J].Mod Rheumatol，2018，28（4）：661669.

[27] GRACIAS D T，STELEKATI E，HOPE J L，et al.The microRNA miR155 controls CD8（+） T cell responses by regulating interferon signaling[J].Nat Immunol，2013，14（6）：593602.

[28] LEISS H，SALZBERGER W，JACOBS B，et al.MicroRNA 155deficiency leads to decreased autoantibody levels and reduced severity of nephritis and pneumonitis in pristaneinduced lupus[J].PLoS One，2017，12（7）：e0181015.

[29] YUAN Y，MA H，YE Z，et al.Interferonstimulated gene 15 expression in systemic lupus erythematosus：Diagnostic value and association with lymphocytopenia[J].Z Rheumatol，2018，77（3）：256262.

[30] KATO H，F(xiàn)UJITA T.RIGIlike receptors and autoimmune diseases[J].Curr Opin Immunol，2015，37：4045.

[31] BUERS I，NITSCHKE Y，RUTSCH F.Novel interferonopathies associated with mutations in RIGI like receptors[J].Cytokine Growth Factor Rev，2016，29：101107.

[32] YOSHIDA H，IMAIZUMI T，LEE S J，et al.Retinoic acidinducible geneI mediates RANTES/CCL5 expression in U373MG human astrocytoma cells stimulated with doublestranded RNA[J].Neurosci Res，2007，58（2）：199206.

[33] IMAIZUMI T，TANAKA H，TAJIMA A，et al.Retinoic acidinducible geneI （RIGI） is induced by IFN{Gamma} in human mesangial cells in culture：possible involvement of RIGI in the inflammation in lupus nephritis[J].Lupus，2010，19（7）：830836.

[34] SU Y J，CHIU W C，KUO H C.Inverse association between antiviral immunity and lupus disease activity[J].Viral Immunol，2018，31（10）：689694.

[35] FUNABIKI M，KATO H，MIYACHI Y，et al.Autoimmune disorders associated with gain of function of the intracellular sensor MDA5[J].Immunity，2014，40（2）：199212.

[36] MUNROE M E，PEZANT N，BROWN M A，et al.Association of IFIH1 and proinflammatory mediators：Potential new clues in SLEassociated pathogenesis[J].PLoS One，2017，12（2）：e0171193.

[37] KIROU K A，LEE C，GEORGE S，et al.Coordinate overexpression of interferonalphainduced genes in systemic lupus erythematosus[J].Arthritis Rheum，2004，50（12）：39583967.

[38] ZHAO M，ZHOU Y，ZHU B C，et al.IFI44L promoter methylation as a blood biomarker for systemic lupus erythematosus[J].Ann Rheum Dis，2016，75（11）：19982006.

[39] WECKERLE C E，F(xiàn)RANEK B S，KELLY J A，et al.Network analysis of associations between serum interferonα activity，autoantibodies，and clinical features in systemic lupus erythematosus[J].Arthritis Rheum，2011，63（4）：10441053.

[40] HU W P，NIU G D，LI H B，et al.The association between expression of IFIT1 in podocytes of MRL/lpr mice and the renal pathological changes it causes：an animal study[J].Oncotarget，2016，7（47）：7646476470.

[41] WANG J H，DAI M，CUI Y G，et al.Association of abnormal elevations in IFIT3 with overactive cyclic GMPAMP synthase/Stimulator of interferon genes signaling in human systemic lupus erythematosus monocytes[J].Arthritis & Rheumatology（Hoboken，N.J.），2018，70（12）：20362045.

[42] XIAO N Y，WEI J J，XU S，et al.cGAS activation causes lupuslike autoimmune disorders in a TREX1 mutant mouse model[J].J Autoimmun，2019，100：8494.

（收稿日期：2020-04-10修回日期：2020-05-08）

（編輯：潘明志）