侯永旺 楊志聰 史 麗（河北北方學(xué)院附屬第一醫(yī)院，張家口 075000）

1 型糖尿?。╰ype 1 diabetes mellitus，T1DM）是兒童最常見的慢性病之一，最常見的發(fā)病年齡為4～6歲和青春期早期（10～14歲）。T1DM的發(fā)病率和患病率一直穩(wěn)步上升，占糖尿病患者的5%～10%。在世界范圍內(nèi)，發(fā)病率也存在相當大的地域差異。據(jù)報道，發(fā)病率最高的是芬蘭和其他北歐國家，其發(fā)病率約是報告發(fā)病率最低的中國和委內(nèi)瑞拉的400倍［1-3］。在T1DM中，胰島中的β細胞在數(shù)月或數(shù)年內(nèi)被免疫破壞，導(dǎo)致胰島素絕對缺乏。目前研究人員認為該病存在遺傳傾向，與特定HLA（DR 和DQ）等位基因密切相關(guān)，且T1DM 具有家族遺傳傾向［4-8］。目前公認的T1DM 與自身免疫應(yīng)答密切相關(guān)，最近有學(xué)者發(fā)現(xiàn)miR-150 缺失通過細胞免疫應(yīng)答促進T1DM 發(fā)生，TNFSF14在STZ 誘導(dǎo)的T1DM 中發(fā)揮重要作用，其可能通過下調(diào)CD4+FoxP3+Treg 分化介導(dǎo)T1DM 發(fā)生，這些研究說明免疫應(yīng)答異常與T1DM 密切相關(guān)［9-10］。現(xiàn)有的預(yù)測T1DM 的主要血清標志物包括抗β細胞自身抗體、抗胰島素的自身抗體（IAA）、抗胰島素瘤相關(guān)抗原-2 的自身抗體（IA-2）、抗谷氨酸脫羧酶的自身抗體（GAD）、抗鋅轉(zhuǎn)運蛋白8的自身抗體（ZnT8）和胰島細胞抗體（ICA），但有些抗體的陽性率不是很高［11-13］。因此探尋新的分子標志物對于T1DM的預(yù)測和診斷治療具有重要意義。

本文主要通過生物信息學(xué)手段篩選出T1DM 患者中的71 個差異表達基因（differentially expressed genes，DEGs），最終確定免疫相關(guān)基因造血細胞信號轉(zhuǎn)導(dǎo)子（hematopoietic cell signal transducer，HCST）為T1DM 的關(guān)鍵基因。研究發(fā)現(xiàn)HCST 可以編碼跨膜信號適配器，并與C型凝集素樣受體NKG2D 形成免疫識別受體復(fù)合物，該復(fù)合物可能通過激活樹突狀細胞（dendritic cells，DCs）、自然殺傷（natural killer，NK）細胞和T 細胞在細胞的存活和增殖中發(fā)揮作用［14-15］。因此，HCST 可能是T1DM 的有效分子標志物。

1 資料與方法

1.1 資料 從GEO 數(shù)據(jù)庫搜索T1DM 和健康人群的外周血單個核細胞微陣列表達數(shù)據(jù)。通過篩選得到GPL570 平臺的GSE55098 數(shù)據(jù)集。該數(shù)據(jù)集中各包含10例健康人群和12例T1DM。

1.2 方法

1.2.1 篩選DEGs 從GEO 數(shù)據(jù)庫下載矩陣文件，借 助GEO2R 分 析DEGs（P<0.05，log2|FC| ≥1）。RStudio 軟件加載“gplots”“heatmap”“ggplot2”做熱圖和火山圖可視化分析。

1.2.2 蛋白互作（protein-protein interaction，PPI）網(wǎng)絡(luò)和關(guān)鍵基因篩選 利用String（https：//string-db.org/）對DEGs 進行PPI 分析，設(shè)置combined sore>0.4對分析結(jié)果進行篩選，并利用Cytoscape 軟件（v3.8.0）進行可視化。采用MCODE 插件對PPI 進行加權(quán)網(wǎng)絡(luò)分析，篩選關(guān)鍵的子網(wǎng)絡(luò)和基因，并對子網(wǎng)絡(luò)做富集分析。利用cytoHubba 插件對PPI 網(wǎng)絡(luò)進行關(guān)鍵基因篩選，對4 種計算方式［MCC（maximun clique centrality）、DMNC（maximum neighborhood component）、MNC（maximun neighborhood component）和Degree］得到的前5 個關(guān)鍵基因進行Venn分析，確定最終的關(guān)鍵基因。

1.2.3 富集分析 借助DAVID 在線分析工具對PPI子網(wǎng)絡(luò)基因進行GO（gene ontology）分析［包括生物過程（biological process，BP）、細胞組分（cellular component，CC）、分子功能（molecular function，MF）］和KEGG（Kyoto Encyclopedia of Genes and Genomes）通路富集分析。富集分析的結(jié)果通過在線網(wǎng)站（http：//www.bioinformatics.com.cn）進行可視化。

1.3 統(tǒng)計學(xué)處理 使用GraphPad Prims 9 軟件對兩組關(guān)鍵基因表達量進行統(tǒng)計分析和受試者工作特征曲線（receiver operating characteristic curve，ROC）分析，兩組間比較采用t檢驗，P<0.05 表示差異具有統(tǒng)計學(xué)意義。

2 結(jié)果

2.1 篩選DEGs 通過P<0.05 和log2|FC|≥1 的條件對GSE55098 數(shù)據(jù)集進行篩選，T1DM 與對照組相比共得到71 個DEGs，其中有上調(diào)基因22 個，下調(diào)基因49個（圖1）。

圖1 DEGs的篩選Fig.1 Screen of DEGs

2.2 PPI 網(wǎng)絡(luò)構(gòu)建及子網(wǎng)絡(luò)功能富集 PPI 分析結(jié)果顯示，String 設(shè)置0.4 的PPI 信任值，且僅顯示27 個節(jié)點，58 個互作關(guān)系（圖2A）。利用MCODE 插件篩選1 個關(guān)鍵子網(wǎng)絡(luò)，得分7.5，共有9 個節(jié)點和30 個互作關(guān)系（圖2B）。通過對子網(wǎng)絡(luò)的9 個基因進行GO 和KEGG 富集分析發(fā)現(xiàn)，這9 個基因主要參與調(diào)節(jié)免疫反應(yīng)、細胞的防御反應(yīng)、免疫應(yīng)答及細胞溶解的生物過程，涉及的主要信號通路是NK 細胞介導(dǎo)的細胞毒作用的信號通路（圖2C、表1）。

表1 GO和KEGG富集分析Tab.1 GO and KEGG enrichment analysis

圖2 PPI網(wǎng)絡(luò)構(gòu)建及子網(wǎng)絡(luò)功能富集Fig.2 PPI network construction and subnetwork function enrichment

2.3 關(guān)鍵基因篩選 MCC 算法得到前5 個關(guān)鍵基因為有PRF1、CCL4、KLRD1、HCST 和GNLY；DMNC算法得到的前5 個關(guān)鍵基因為HCST、GNLY、GZMH、NKG7 和NCAM1；MNC 算法得到的前5 個關(guān)鍵基因PRF1、KLRD1、CCL4、NKG7 和HCST；Degree算法得到的前5 個關(guān)鍵基因為PRF1、CCL4、KLRD1、NKG7 和HCST（圖3A）。4 種算法得到的關(guān)鍵基因通過Venn 分析最終篩選出1 個關(guān)鍵基因，即HCST（圖3B）。

圖3 關(guān)鍵基因篩選Fig.3 Screen of hub genes

2.4 HCST 在T1DM 中的表達量及ROC 曲線 對GSE55098 數(shù)據(jù)集中10 例健康人群和12 例T1MD 患者的HCST 表達量進行對比分析，結(jié)果顯示HCST 在T1DM 組的表達量明顯低于健康對照組，差異有統(tǒng)計學(xué)意義（圖4A）。HCST 的ROC 曲線分析結(jié)果顯示曲線下面積為0.983 3，P=0.0001，差異有統(tǒng)計學(xué)意義（圖4B）。

3 討論

T1MD 是由胰腺β 細胞與免疫系統(tǒng)間的復(fù)雜相互作用所引起。但是機體中針對β細胞的免疫反應(yīng)是隨機的還是存在誘發(fā)因素尚存在爭議［16-17］。一種長期存在的假說認為，病毒感染是T1DM 自身免疫反應(yīng)的潛在觸發(fā)因素。在病毒中，腸道病毒是研究最多的。糖尿病預(yù)測與預(yù)防研究證實了腸道病毒感染與第1個自身抗體出現(xiàn)間的關(guān)系［18］。但目前大多數(shù)對于T1DM 免疫功能障礙的研究集中在下游T細胞和B 細胞的適應(yīng)性反應(yīng)，而忽視了上游先天參與者如NK細胞的作用。研究發(fā)現(xiàn)，與T1DM 有關(guān)的許多病毒是慢性潛伏或溶原性感染，它們通過類似于妊娠耐受的機制降低NK細胞的反應(yīng)和數(shù)量。NK細胞對感染的抑制反應(yīng)導(dǎo)致不正確的信號傳導(dǎo)，這種不適當?shù)目乖蔬f導(dǎo)致了CD8+淋巴細胞的增殖和細胞毒性［19-20］。

NK 細胞對病原體做出反應(yīng)并直接將其殺死，此外，NK 細胞還參與了適應(yīng)性免疫反應(yīng)的強度和精細控制。最近的研究表明，NK 細胞控制下游反應(yīng)的效應(yīng)和抑制活性，包括激活或殺死抗原提呈細胞和Treg、細胞毒性T 淋巴細胞（CTLs）、T 輔助細胞（Th）和B 細胞［21-23］。NK 細胞的作用是調(diào)節(jié)細胞毒性CD8+T 淋巴細胞反應(yīng)，從而控制異常或慢性炎癥反應(yīng)，避免未識別的細胞和組織遭到破壞［24］。EB病毒被證明能有效抑制NK 細胞，從而抑制NK 細胞對細胞毒性CD8+T 淋巴細胞反應(yīng)的調(diào)節(jié)［25-26］。研究發(fā)現(xiàn)，在NOD 模型中，NK 細胞介導(dǎo)CD8+細胞毒性T 淋巴細胞的破壞可改善自身免疫性糖尿?。?7］。

本文通過篩選GSE55098 數(shù)據(jù)集中T1DM 的差異基因并構(gòu)建PPI網(wǎng)絡(luò)，借助MCODE 插件篩選關(guān)鍵子網(wǎng)絡(luò)，通過對子網(wǎng)絡(luò)基因功能富集分析發(fā)現(xiàn)HCST及相關(guān)基因主要富集在NK 細胞導(dǎo)致的細胞毒信號途徑，并參與調(diào)節(jié)免疫反應(yīng)、細胞的防御反應(yīng)、免疫應(yīng)答及細胞溶解的生物過程。同時，本研究結(jié)果表明HCST 在T1DM 中表達量顯著降低。因此課題組認為HCST 基因的降低引起了NK 細胞參與的固有免疫及相關(guān)信號通路的功能發(fā)生紊亂，這使得在發(fā)生病毒感染或存在其他誘導(dǎo)因素的情況下，NK 細胞的錯誤遞呈作用導(dǎo)致了異常的細胞毒作用，最終胰島β 細胞遭到慢性損傷并導(dǎo)致T1DM 發(fā)生。這些結(jié)果與BAYER 等［19］的觀點相似，進一步證實了本研究結(jié)果。ROC 曲線結(jié)果表明，HCST 可以成為診斷T1DM的潛在分子標志物。

綜上所述，HCST 在T1DM 中表達量降低，或可成為診斷T1DM 的潛在分子標志物。此外，HCST 引起了NK 細胞參與的固有免疫及相關(guān)信號通路的功能紊亂，導(dǎo)致T1DM 發(fā)生，因此HCST 可能成為治療T1DM的潛在靶點。