Identification and validation of key molecules associated with humoral immune modulation in Parkinson's disease based on bioinformatics
基于生物信息學(xué)的帕金森病體液免疫調(diào)節(jié)相關(guān)關(guān)鍵分子的鑒定和驗(yàn)證
發(fā)表期刊：Front Immunol
發(fā)表日期：2022 Sep 15
影響因子：8.786
DOI:? 10.3389/fimmu.2022.948615
一匙隔、研究背景
????????隨著醫(yī)學(xué)發(fā)展的進(jìn)步，人們的預(yù)期壽命大大增加熏版；特別是隨著老齡化程度的增加纷责，與年齡相關(guān)的神經(jīng)退行性疾病的發(fā)病率也隨之提高。帕金森埠扯獭（PD）是最常見的神經(jīng)退行性運(yùn)動(dòng)障礙再膳，其發(fā)病年齡通常在65-70歲之間。
????????神經(jīng)免疫性炎癥可以保護(hù)神經(jīng)元免受損害曲横，但其毒性作用也會(huì)加劇神經(jīng)元的紊亂喂柒，這主要是由細(xì)胞因子和趨化因子調(diào)節(jié)的中樞神經(jīng)系統(tǒng)（CNS）的先天免疫細(xì)胞和外周侵入性免疫細(xì)胞引起的。神經(jīng)退行性疾病中血腦屏障（BBB）的破壞使血液中的神經(jīng)毒性物質(zhì)禾嫉、細(xì)胞和微生物病原體進(jìn)入大腦灾杰，這些與炎癥和免疫反應(yīng)有關(guān)，可以啟動(dòng)各種神經(jīng)退行性疾病的途徑熙参。炎癥性中樞神經(jīng)系統(tǒng)細(xì)胞表達(dá)某些粘性分子艳吠，吸引更多的外周免疫細(xì)胞和抗體通過BBB進(jìn)入中樞神經(jīng)系統(tǒng)。
二孽椰、材料與方法
1昭娩、數(shù)據(jù)來源
1） 基因表達(dá)數(shù)據(jù)GSE100054（10個(gè)PD和9個(gè)HC）和GSE49126（30個(gè)PD和20個(gè)HC）
2） GSE22491（10個(gè)PD和8個(gè)HC）被用來作為驗(yàn)證數(shù)據(jù)集
3） 從三位神經(jīng)退行性疾病專家那里招募了69名PD患者凛篙，以及36名年齡和性別匹配的健康對照者
2、分析流程
1） 差異表達(dá)基因的鑒定：limma栏渺，采用線性模型來評(píng)估PD患者和健康對照組之間的差異表達(dá)呛梆；所有的DEGs與疾病數(shù)據(jù)庫DisGeNET中記錄的與帕金森病相關(guān)的基因相交，得到感興趣的DEGs
2） PPI網(wǎng)絡(luò)的構(gòu)建和模塊分析：蛋白質(zhì)-蛋白質(zhì)相互作用網(wǎng)絡(luò)（PPI）是用STRING構(gòu)建的
3） DEGs和重要模塊的KEGG和GO富集度分析：Fisher's精確檢驗(yàn)對GO分類進(jìn)行評(píng)估迈嘹，使用關(guān)鍵詞 "免疫"削彬，確定了39條參與免疫相關(guān)生物過程的途徑，共獲得13個(gè)分子
4） 對PD分子特征的驗(yàn)證：在數(shù)據(jù)集GSE22491中秀仲，比較了PD患者和健康對照組中13個(gè)分子的表達(dá)水平融痛；選擇PPBP、PROS1和LCN2三個(gè)分子進(jìn)行實(shí)驗(yàn)驗(yàn)證
5） PBMC分離神僵、RNA提取和反轉(zhuǎn)錄雁刷、定量實(shí)時(shí)聚合酶鏈反應(yīng)、酶聯(lián)免疫吸附試驗(yàn)
三保礼、實(shí)驗(yàn)結(jié)果
01 - 識(shí)別差異表達(dá)的基因
????????作者從GEO數(shù)據(jù)庫下載了微陣列表達(dá)數(shù)據(jù)集沛励。根據(jù)上述選擇標(biāo)準(zhǔn)進(jìn)行差異表達(dá)分析后，在GSE100054中得到373個(gè)DEGs炮障，包括337個(gè)上調(diào)基因和36個(gè)下調(diào)基因目派，在GSE49126中得到384個(gè)DEGs，包括200個(gè)上調(diào)基因和184個(gè)下調(diào)基因胁赢。進(jìn)行主成分分析（PCA）并以散點(diǎn)圖顯示企蹭，其中每個(gè)點(diǎn)代表一個(gè)樣本（圖2A，C）智末，在火山圖中以不同顏色顯示篩選出的組間差異明顯的基因（圖2B谅摄，D）。用GSE100054和GSE49126來識(shí)別不同表達(dá)的基因（DEGs）系馆，然后分別與DisGeNET數(shù)據(jù)庫記錄的PD相關(guān)基因相交（圖3A）送漠。
????????為了描述差異表達(dá)基因所編碼的蛋白質(zhì)之間的相互作用，從STRING數(shù)據(jù)庫中構(gòu)建了一個(gè)PPI網(wǎng)絡(luò)（圖3B）由蘑。該網(wǎng)絡(luò)由92個(gè)節(jié)點(diǎn)和293條邊組成闽寡；節(jié)點(diǎn)代表感興趣的DEGs，邊代表基因之間的相互作用尼酿。然后將PPI網(wǎng)絡(luò)導(dǎo)入cytoscape進(jìn)行可視化爷狈，并使用MCODE進(jìn)行聚類相關(guān)分析。最后谓媒，得到三個(gè)高分的聚類，即聚類1（17個(gè)節(jié)點(diǎn)何乎，得分8.7句惯，圖3C）土辩、聚類2（7個(gè)節(jié)點(diǎn)，得分4.3抢野，圖3D）和聚類3（3個(gè)節(jié)點(diǎn)拷淘，得分3.0，圖3E）指孤。
02 - DEGs的功能富集和重要模塊
????????基于智人背景的clusterProfiler軟件包對感興趣的DEGs（94個(gè)基因）和Hubgenes（27個(gè)基因）進(jìn)行了富集分析启涯，以獲得富集信息。通過GO分類和富集分析恃轩，如生物過程结洼、細(xì)胞成分、分子功能等叉跛，揭示了與PD相關(guān)的功能注釋詞松忍，并探討了其生物學(xué)意義。為了揭示與感興趣的DEGs和Hubgenes相關(guān)的生物學(xué)途徑筷厘，利用KEGG數(shù)據(jù)庫獲得了重要的信號(hào)通路鸣峭。
????????DEGs富集的生物學(xué)過程主要包括炎癥反應(yīng)的調(diào)節(jié)、血小板脫顆粒酥艳、免疫效應(yīng)過程的調(diào)節(jié)摊溶、免疫反應(yīng)分子介質(zhì)的產(chǎn)生、白細(xì)胞增殖充石、參與免疫反應(yīng)的細(xì)胞因子產(chǎn)生莫换、神經(jīng)炎癥反應(yīng)等（圖4A）。在 "細(xì)胞成分"方面赫冬，DEGs主要富集在囊泡腔浓镜、血小板α顆粒腔、質(zhì)膜外側(cè)劲厌、運(yùn)輸囊泡等區(qū)域（圖4B）膛薛；在 "分子功能"方面，DEGs主要與細(xì)胞因子活性补鼻、受體配體活性哄啄、NAD+核苷酶活性、白介素-1受體結(jié)合等有關(guān)（圖4C）风范。在KEGG的富集分析中咨跌，DEGs主要涉及脂質(zhì)和動(dòng)脈硬化、MAPK信號(hào)通路硼婿、細(xì)胞因子-細(xì)胞因子受體相互作用锌半、HIF-1信號(hào)通路、炎癥性腸病和膽固醇代謝等（圖4D）寇漫。
????????Hubgenes富集的生物過程主要包括參與免疫反應(yīng)的細(xì)胞因子產(chǎn)生刊殉、炎癥反應(yīng)的調(diào)節(jié)殉摔、血小板脫顆粒山憨、膠質(zhì)細(xì)胞激活顶燕、細(xì)胞因子產(chǎn)生的正向調(diào)節(jié)瓣窄、免疫效應(yīng)過程的調(diào)節(jié)等（圖5A）狞贱。在 "細(xì)胞成分"方面勒虾，Hubgenes主要富集在囊泡腔捷雕、分泌顆粒腔吝秕、膜區(qū)和其他區(qū)域（圖5B）绵咱。在 "分子功能"方面瓢颅，表明Hubgenes主要與細(xì)胞因子活性恩尾、細(xì)胞因子受體結(jié)合、淀粉樣蛋白-β結(jié)合惜索、C-C趨化因子受體活性等有關(guān)（圖5C）特笋。在KEGG的富集分析中，Hubgenes主要涉及類風(fēng)濕性關(guān)節(jié)炎巾兆、細(xì)胞因子-細(xì)胞因子受體相互作用猎物、癌癥中的蛋白多糖、PI3K-Akt信號(hào)通路角塑、壞死蔫磨、糖尿病并發(fā)癥的AGE信號(hào)通路等（圖5D）。
03 - 鑒定和驗(yàn)證PD的分子特征
????????作者基于27個(gè)Hubgenes（以 "免疫 "為關(guān)鍵詞）進(jìn)行了基因本體論（GO）的生物過程分析圃伶，篩選出39條免疫相關(guān)的通路堤如，并從上述Hubgenes中選出13個(gè)靶分子，如CCR2窒朋、CD 44搀罢、TLR2、TLR4侥猩、IL4榔至、IL-1B、IL-33欺劳、TGFB2唧取、PPBP、PROS1划提、CLU枫弟、LCN2和NLRP3。GSE22491被選為驗(yàn)證數(shù)據(jù)集以比較這些分子在PD和HC中的表達(dá)水平鹏往。IL-1B淡诗、TGFB2、CLU、IL33和NLRP3的表達(dá)在兩組之間沒有明顯差異韩容；其他基因被認(rèn)為有明顯差異（圖6）绪爸。根據(jù)以往的研究和實(shí)驗(yàn)實(shí)施的客觀條件，選擇了PPBP宙攻、LCN2和PROS1進(jìn)行以下實(shí)驗(yàn)。與HC相比介褥，這三個(gè)分子在PD中都有所下降（圖6）座掘；但在GSE100054中，PPBP柔滔、PROS1和LCN2的表達(dá)水平高于HC（圖7A-C）溢陪。在免疫相關(guān)的生物過程中，PPBP和LCN2主要參與體液免疫反應(yīng)和中性粒細(xì)胞激活參與免疫反應(yīng)睛廊；然而形真，除了體液免疫反應(yīng)，PROS1也參與調(diào)節(jié)體液免疫反應(yīng)和免疫效應(yīng)器過程超全。
04 - 在轉(zhuǎn)錄水平上對選定的差異分子進(jìn)行驗(yàn)證
????????作者招募了69名PD患者（39名男性和30名女性）和36名健康對照（19名男性和17名女性）。從這些參與者中嘶朱，首先蛾坯，選擇了一些進(jìn)行mRNA表達(dá)水平分析。觀察到疏遏，與健康對照組相比脉课，PD患者的PBMC中LCN2的表達(dá)明顯增加；但是财异，PPBP的mRNA水平下降（圖8B倘零，C）。PROS1的表達(dá)水平相當(dāng)?shù)痛链纾虼藷o法評(píng)估呈驶。此外，還從全血中提取了RNA庆揩，PROS1的表達(dá)在PD中明顯增加（圖8A）俐东，而LCN2和PPBP在組間沒有明顯差異。
05 - 在翻譯水平上對選定的差異分子進(jìn)行驗(yàn)證
????????對于ELISA虏辫，上述三個(gè)分子的血漿水平與mRNA表達(dá)一致。研究結(jié)果顯示锈拨，PD患者和健康對照組的外周血樣本中PPBP水平分別為1.770微克/毫升和2.269微克/毫升（圖9A）砌庄。與健康對照組相比，檢測結(jié)果顯示，PD患者血漿中PROS1和LCN2的濃度均有所提高（圖9B, C）娄昆。應(yīng)用接收操作特征（ROC）曲線來確定上述三種分子的診斷效果佩微，發(fā)現(xiàn)PPBP、PROS1和LCN2的曲線下面積（AUC）分別為0.663萌焰、0.674和0.885（圖10A-C）哺眯。
四奶卓、結(jié)論
????????在這項(xiàng)研究中，作者利用基因表達(dá)數(shù)據(jù)集來識(shí)別DEGs撼玄，并通過構(gòu)建PPI網(wǎng)絡(luò)和識(shí)別關(guān)鍵模塊發(fā)現(xiàn)了13個(gè)感興趣的基因夺姑。基于GO生物過程分析掌猛，PPBP盏浙、PROS1和LCN2參與了體液免疫反應(yīng)。盡管用Q-PCR和ELISA方法證明了PPBP荔茬、PROS1和LCN2在帕金森病患者和健康對照組中的統(tǒng)計(jì)學(xué)差異废膘，但還需要進(jìn)一步的調(diào)查，以便更深入地了解涉及這三種分子的外周體液免疫反應(yīng)對帕金森病潛在機(jī)制的影響慕蔚。
