影響因子：8.7
關(guān)于非腫瘤生信总珠，我們也解讀過很多办绝，主要有以下類型
1 單個疾病分析
2 單個疾病結(jié)合免疫浸潤，鐵死亡姚淆，自噬等相關(guān)分析
3 兩種相關(guān)疾病聯(lián)合分析孕蝉，包括非腫瘤結(jié)合非腫瘤，非腫瘤結(jié)合腫瘤或者非腫瘤結(jié)合泛癌分析
目前非腫瘤生信發(fā)文的門檻較低腌逢，有需要的朋友歡迎交流降淮！
研究概述
骨關(guān)節(jié)炎（OA）是一種對老年人影響較大的退行性關(guān)節(jié)疾病，以關(guān)節(jié)軟骨的退變搏讶、變性為典型表現(xiàn)佳鳖，常引起關(guān)節(jié)軟骨丟失及軟骨下骨改變，患病率隨著年齡的增長而上升媒惕。目前其病因及發(fā)病機(jī)制仍不十分清楚系吩。
在這項研究中，作者使用生物信息學(xué)方法識別了與OA免疫相關(guān)的生物標(biāo)志物妒蔚，確定了兩種不同的分子亞型（ClusterA和ClusterB）穿挨，并對不同分子亞型的預(yù)后月弛、免疫浸潤和功能途徑進(jìn)行了分析】剖ⅲ總之帽衙，本研究對于OA免疫微環(huán)境的早期診斷、治療和預(yù)后具有重要的臨床意義贞绵。
流程圖
研究結(jié)果
一厉萝、識別OA免疫相關(guān)差異表達(dá)基因
1、從GEO數(shù)據(jù)庫下載微陣列數(shù)據(jù)集榨崩，通過主成分分析（PCA）處理后谴垫，篩選得到69個OA樣本和34個對照樣本（下圖）。
2母蛛、通過差異表達(dá)分析確定OA樣本和對照樣本之間的差異基因（DEGs）弹渔，其中83個基因表達(dá)上調(diào)，22個基因表達(dá)下調(diào)（下圖A-B）溯祸。通過交叉DEGs和免疫基因識別差異表達(dá)的免疫基因（下圖C）肢专。這些交叉點產(chǎn)生了28個Imm-DEG（下圖D），其中26個基因表達(dá)上調(diào)（下圖E）焦辅，2個基因表達(dá)下調(diào)（下圖F）博杖。
3、繪制OA樣本和對照樣本中Imm-DEGs表達(dá)差異的熱圖（下圖A）筷登。結(jié)果顯示剃根，大多數(shù)基因在OA樣本中的表達(dá)水平高于正常樣本（下圖C）。
二前方、構(gòu)建風(fēng)險預(yù)后特征
1狈醉、利用LASSO回歸從28個Imm-DEG中識別出16個對OA有較大影響的特征基因（下圖A-B）。
2惠险、根據(jù)這16個特征基因的系數(shù)得出風(fēng)險評分（下圖C-D）苗傅，ROC曲線結(jié)果顯示這16基因特征的預(yù)測效能較好（下圖E）。
3班巩、根據(jù)風(fēng)險評分和16特征基因創(chuàng)建列線圖（下圖A-B）渣慕。校準(zhǔn)曲線和決策曲線分析顯示，該列線圖的臨床預(yù)測性能較為可靠（下圖C-D）抱慌。
4逊桦、分析16基因特征與基因表達(dá)水平之間的相關(guān)性。結(jié)果顯示抑进，在所有樣本中强经，IGKV1-17與IGLC1的相關(guān)系數(shù)為0.95；ADIPOQ與LEP呈正相關(guān)寺渗；LEP與除ADIPOQ外的所有其他基因呈負(fù)相關(guān)（下圖A）匿情；其中IGKV1-17兰迫、IGLC1和其他基因之間的功能相關(guān)系數(shù)達(dá)到0.9（下圖B）。
5码秉、在OA樣本中，IGKV1-17與IGLC1的相關(guān)系數(shù)為0.94鸡号；ADIPOQ與LEP和SPP1呈正相關(guān)转砖；LEP與除ADIPOQ之外的所有其他基因呈負(fù)相關(guān)（下圖C）；OGN與大多數(shù)其他特征基因呈負(fù)相關(guān)鲸伴；其中IGKV1-17和IGLC1之間的相關(guān)系數(shù)為0.89（下圖D）府蔗。
三、16特征基因的不同免疫學(xué)亞型
1汞窗、基于16特征基因的不同表達(dá)模式建立了兩種免疫亞型（clusterA和clusterB）姓赤。其中clusterA中有34個樣本，clusterB中有35個樣本仲吏。
2不铆、構(gòu)建所有差異表達(dá)免疫基因的熱圖（下圖A）。結(jié)果顯示裹唆，ADIPOQ誓斥、LEP、OGN和TAC1在clusterA中的表達(dá)水平顯著低于clusterB许帐；而HLA-DMA劳坑、HLA-DMB、HLA-DPA1成畦、HLA-DQA1距芬、HLA-DQB1、HLA-DRA循帐、TLR7框仔、CCL18、CXCL10拄养、MMP9存和、BLNK、IGHM衷旅、IGKV1-17捐腿、IGLC1、IGLV1-44柿顶、CXCL13茄袖、CXCL6、CXCL9嘁锯、SPP1宪祥、IL10RA聂薪、SDC1、TNFRSF17和TRBC1在clusterA中的表達(dá)水平顯著高于clusterB（下圖B)蝗羊。
3藏澳、評估這兩種免疫亞型的ROC曲線。結(jié)果表明耀找，這16個基因都可以獨立區(qū)分clusterA和clusterB（下圖C）翔悠。
四、構(gòu)建PPI免疫基因網(wǎng)絡(luò)
1野芒、構(gòu)建DEG蓄愁、Imm-DEG和16基因特征的PPI網(wǎng)絡(luò)以探索差異表達(dá)的免疫基因之間的關(guān)系。結(jié)果顯示狞悲，DEG的PPI網(wǎng)絡(luò)包含211個相互作用對和75個基因撮抓，其中MMP9與19個DEG密切相關(guān)，而CXCL10與16個DEG相關(guān)（下圖A）摇锋；Imm-DEGs的PPI網(wǎng)絡(luò)包含58個相互作用對和24個基因丹拯，其中CXCL10和MMP9與10個差異表達(dá)的免疫基因密切相關(guān)，而CXCL13和HLA-DRA與8個差異表達(dá)的免疫基因相關(guān)（下圖C）荸恕；基因特征的PPI網(wǎng)絡(luò)包含25個相互作用對和15個基因咽笼，其中CXCL10和MMP9與10個Imm-DEG密切相關(guān)，而CXCL13和HLA-DRA與8個Imm-DEG相關(guān)（下圖E）戚炫。
2剑刑、利用富集分析驗證這3個PPI網(wǎng)絡(luò)中基因的功能。結(jié)果顯示双肤，DEG主要富集在脂肪細(xì)胞因子信號施掏、糖酵解/糖異生、Toll樣受體茅糜、細(xì)胞因子受體IL-17信號以及類風(fēng)濕性關(guān)節(jié)炎等通路中（下圖B）七芭。差異表達(dá)的免疫基因主要富集在腸道免疫網(wǎng)絡(luò)中IgA的產(chǎn)生、病毒蛋白與細(xì)胞因子和細(xì)胞因子受體的相互作用等途徑中（下圖D）蔑赘。
五狸驳、兩種免疫亞型的差異分析
1、在兩種免疫亞型中獲得了168個差異表達(dá)基因缩赛。其中cluster A中有93個上調(diào)基因和75個下調(diào)基因（下圖A）耙箍。熱圖顯示這些DEG可以區(qū)分兩種不同的免疫亞型（下圖B）。
2酥馍、分析兩種免疫亞型之間的DEG的功能（下圖A）辩昆。結(jié)果顯示，這些DEG富集在細(xì)胞因子-細(xì)胞因子受體相互作用旨袒、趨化因子信號傳導(dǎo)汁针、瘧疾和TNF信號傳導(dǎo)以及NOD樣受體信號通路中（下圖B）术辐。
3、GSEA分析結(jié)果顯示施无，骨關(guān)節(jié)炎主要與免疫反應(yīng)激活和脂肪酸代謝密切相關(guān)辉词。
六、兩種免疫亞型之間的免疫特征差異
1猾骡、利用CIBERSORT算法評估兩種不同免疫亞型之間的免疫細(xì)胞浸潤水平瑞躺。分析顯示，cluster A中M2巨噬細(xì)胞卓练、靜息狀態(tài)下樹突細(xì)胞和肥大細(xì)胞隘蝎、活化的NK細(xì)胞和調(diào)節(jié)性T細(xì)胞的水平顯著低于cluster B（下圖A）购啄。然而襟企，M1巨噬細(xì)胞、活化肥大細(xì)胞狮含、漿細(xì)胞和T濾泡輔助細(xì)胞的水平顯著高于cluster B（下圖B）顽悼。
2、計算cluster A和cluster B中的免疫細(xì)胞含量几迄。分析顯示蔚龙，cluster A中的T濾泡輔助細(xì)胞與活化的記憶性CD4+T細(xì)胞呈顯著正相關(guān)，而與活化和靜息狀態(tài)下的肥大細(xì)胞映胁、M2巨噬細(xì)胞和B細(xì)胞呈顯著負(fù)相關(guān)（下圖C）木羹；在cluster B中，γδT細(xì)胞解孙、M1巨噬細(xì)胞坑填、T濾泡輔助細(xì)胞和活化的肥大細(xì)胞與調(diào)節(jié)性T細(xì)胞呈正相關(guān)，而調(diào)節(jié)性T細(xì)胞與M1巨噬細(xì)胞呈負(fù)相關(guān)（下圖D）弛姜。
3脐瑰、計算所有OA樣本（下圖E）、clusterA（下圖F）和cluster B（下圖G）樣本中的16種基因特征和免疫細(xì)胞類型之間的相關(guān)性廷臼。結(jié)果顯示苍在，TAC1與cluster B患者的免疫細(xì)胞含量無關(guān)，記憶性CD4+T細(xì)胞和活化的NK細(xì)胞與基因特征無關(guān)荠商。
七寂恬、qRT-PCR驗證hub基因的表達(dá)水平
1、實驗結(jié)果顯示：IL-1β誘導(dǎo)后莱没，TCA1掠剑、TLR7、MMP9郊愧、CXCL10朴译、CXCL13井佑、HLA-DRA、ADIPOQ和SPP1的mRNA表達(dá)水平在正常人（下圖A）和骨關(guān)節(jié)炎患者（下圖B）中均顯著升高眠寿，其中CXCL10的上升趨勢最明顯躬翁。
總結(jié)
在這項研究中，作者鑒定出骨關(guān)節(jié)炎的關(guān)鍵生物標(biāo)志物盯拱，確定了兩種不同的免疫亞型（ClusterA和ClusterB）盒发，并使用CIBERSORT等方法分析了不同免疫亞型之間的免疫細(xì)胞浸潤水平的差異。最后實驗結(jié)果顯示狡逢，TCA1宁舰、TLR7、MMP9奢浑、CXCL10蛮艰、CXCL13、HLA-DRA和ADIPOQSPP1在IL-1β誘導(dǎo)組中的表達(dá)水平顯著高于骨關(guān)節(jié)炎組雀彼，還表明骨關(guān)節(jié)炎可能與免疫反應(yīng)有關(guān)壤蚜。
