腎癌免疫細胞浸潤模式與生存預(yù)后分析*

華金駿，李 興，王 杰△

上海中醫(yī)藥大學附屬曙光醫(yī)院：1.泌尿外科；2.麻醉科，上海 201210

腎癌是泌尿系統(tǒng)常見的腫瘤，在我國發(fā)病率呈逐年升高趨勢，約25%的患者就診時已經(jīng)出現(xiàn)了轉(zhuǎn)移，轉(zhuǎn)移后的患者5年生存率不足10%[1]。雖然手術(shù)治療是腎癌的主要治療手段，但約有1/2的患者術(shù)后出現(xiàn)復發(fā)和預(yù)后不良，療效欠佳，占比最高的腎細胞癌對傳統(tǒng)放化療敏感程度也較低，總體而言，目前對于腎癌仍缺乏有效的針對性治療方法[2]。腎癌組織周圍存在大量T細胞、NK細胞、巨噬細胞等，免疫治療可能成為腎癌治療的有效手段[3]。腫瘤微環(huán)境免疫模式已經(jīng)被列入腫瘤十大特征之一，既能夠預(yù)測患者生存預(yù)后，還能預(yù)測患者放化療療效[4]。腫瘤微環(huán)境中免疫細胞分布模式與患者生存預(yù)后密切相關(guān)，免疫細胞浸潤模式可以作為治療依據(jù)來提高患者生存率與改善遠期預(yù)后[5]。目前對腎癌免疫細胞浸潤及免疫標志物的研究較少，尚缺乏精準篩選腎癌免疫治療標志物的方法。常規(guī)免疫組化和流式細胞檢測技術(shù)在評估免疫細胞模式上存在一定的局限，單次評估免疫細胞效率較低。CIBERSORT是目前生物信息學研究中新的基于基因表達大數(shù)據(jù)去卷積的復雜算法，可以根據(jù)患者的基因表達數(shù)據(jù)轉(zhuǎn)化為免疫細胞浸潤比例[6]。本研究采用生物信息學的方法，挖掘腫瘤癌癥基因圖譜(TCGA)數(shù)據(jù)庫，基于CIBERSORT復雜算法和TCGA腎癌患者樣本數(shù)據(jù)，系統(tǒng)性探討腎癌患者免疫細胞浸潤模式及其與免疫細胞的相關(guān)性，深入分析腎癌免疫細胞浸潤模式與患者生存預(yù)后、臨床分期的關(guān)系。

1 資料與方法

1.1獲取數(shù)據(jù) 進入在線TCGA網(wǎng)站(https://portal.gdc.cancer.gov/)選擇疾病腎癌，數(shù)據(jù)來源方式選擇TCGA，然后選擇轉(zhuǎn)錄本數(shù)據(jù)，工作流數(shù)據(jù)選擇矯正過的FPKM數(shù)據(jù)，下載組織樣本基因表達數(shù)據(jù)[7]，共1 023例，包含正常組織128例，腎癌組織895例。在TCGA網(wǎng)站腎癌中選擇臨床樣本數(shù)據(jù)，選擇TCGA數(shù)據(jù)來源，選擇數(shù)據(jù)格式bcr xml，方便后期程序化處理，下載關(guān)于腎癌的臨床數(shù)據(jù)，包括年齡、性別、生存時間、生存狀態(tài)、腫瘤分期信息，共943例，運用Excel對臨床信息進行篩選，剔除數(shù)據(jù)缺失、分期狀態(tài)未知、生存狀態(tài)未知的患者，保留完整數(shù)據(jù)的臨床信息患者，最終獲得526例有效腎癌患者臨床特征信息。

1.2評估腎癌免疫細胞浸潤 運用Excel對526例腎癌患者臨床信息中的年齡、性別、臨床分期繪制直觀三線表。利用Perl軟件(5.26.1)提取腎癌相關(guān)轉(zhuǎn)錄本中mRNA表達譜數(shù)據(jù)，生成轉(zhuǎn)錄本mRNA的矩陣文件，并進行基因名稱轉(zhuǎn)換標準化注釋[8]。運用R軟件調(diào)用CIBERSORT數(shù)據(jù)包(https://cibersort.stanford.edu/)檢查標準化后的矩陣文件中22種腫瘤免疫細胞浸潤情況[9]。這22種免疫細胞包括M0巨噬細胞、M1巨噬細胞、M2巨噬細胞、濾泡輔助T細胞等，其中幼稚CD4+T細胞在全部樣本中水平為0，故針對性分析其余21種免疫細胞。運用Perl數(shù)據(jù)包，以P<0.05為差異有統(tǒng)計學意義，最終篩選得到10例正常組織和558例腎癌組織免疫細胞浸潤信息。運用Excel分別繪制正常和腎癌組織免疫細胞占比餅圖。運用R軟件包中vioplot 包繪制可視化腎癌組織與正常組織免疫細胞構(gòu)成比小提琴圖，利用R軟件包中的corroplot包分析免疫細胞間的相關(guān)性。將免疫細胞水平矩陣與其對應(yīng)的TCGA臨床數(shù)據(jù)合并，生成包含臨床數(shù)據(jù)和免疫細胞浸潤信息的TXT文本，運用R軟件包中的survival包完成生存分析，以及免疫細胞浸潤與腎癌臨床分期相關(guān)性評價。

1.3統(tǒng)計學處理 總生存期定義為從診斷日期開始到死亡日期之間的間隔，應(yīng)用R軟件包批量繪制K-M生存分析曲線，并進行假設(shè)檢驗(Log-Rank test)分析差異。腎癌分期中兩組樣本的比較采用Wilcoxon秩和檢驗，兩組以上樣本的比較則采用Kruskal-Wallis秩和檢驗，以P<0.05為差異有統(tǒng)計學意義[10]。

2 結(jié) 果

2.1腎癌患者臨床特征 從TCGA下載腎癌臨床相關(guān)數(shù)據(jù)，剔除不完整的數(shù)據(jù)，得到526例腎癌患者臨床信息，包含年齡、性別、臨床分期信息，見表1。

表1 腎癌患者臨床病理特征

2.2腎癌免疫細胞浸潤分布 從TCGA下載到1 023例免疫細胞浸潤信息，經(jīng)過篩選后獲得10例正常組織，558例腎癌組織免疫細胞浸潤信息，繪制正常組織和腎癌組織免疫細胞構(gòu)成比圖及腎癌免疫細胞相關(guān)性熱圖。結(jié)果顯示，腎癌組織與正常組織相比，M2巨噬細胞水平較高(P<0.05)；幼稚B細胞、漿細胞、未活化的CD4+記憶T細胞、靜息態(tài)樹突狀細胞水平較低(P<0.05)。同時，在腎癌組織中，CD8+T細胞與濾泡輔助性T細胞(r=0.59)、γδT細胞水平呈正相關(guān)(r=0.60)；CD8+T細胞與M2巨噬細胞水平(r=-0.55)、未活化的CD4+記憶T細胞水平(r=-0.62)呈負相關(guān)。

2.3腎癌浸潤性免疫細胞水平與臨床相關(guān)性評價

2.3.1與臨床分期的相關(guān)性 將獲得的臨床分期信息與其對應(yīng)的免疫細胞浸潤信息進行合并，繪制腎癌浸潤性免疫細胞與臨床分期箱式圖，見圖1。隨著臨床分期的進展，未活化的NK細胞、單核細胞、未活化的肥大細胞、未活化的CD4+記憶T細胞水平逐漸降低，提示這些細胞抑制腎癌進展；M0巨噬細胞、濾泡輔助性T細胞、CD8+T細胞、活化的CD4+記憶T細胞、漿細胞、調(diào)節(jié)性T細胞水平隨著腎癌進展逐漸遞增，提示這些細胞促進腎癌進展。

圖1 腎癌浸潤性免疫細胞水平與臨床分期箱式圖

2.3.2與生存預(yù)后的相關(guān)性 根據(jù)腎癌免疫細胞構(gòu)成比高和低，結(jié)合假設(shè)檢驗(Log-Rank test)計算其與生存的相關(guān)性，批量繪制腎癌組織免疫細胞K-M生存曲線，見圖2。可以發(fā)現(xiàn)，低水平調(diào)節(jié)性T細胞、濾泡輔助性T細胞的患者5年生存率更高；高水平單核細胞、靜息態(tài)樹突狀細胞、未活化的肥大細胞的患者5年生存率更高。

注：A、B、C、D、E分別代表低水平和高水平調(diào)節(jié)性T細胞、濾泡輔助性T細胞、單核細胞、靜息態(tài)樹突狀細胞、未活化的肥大細胞患者生存率的比較。

3 討 論

目前，相關(guān)學者對腫瘤免疫細胞浸潤的研究和認識逐步提升，發(fā)現(xiàn)某些特定的免疫細胞亞群在預(yù)測個體化特異性治療效果上起著主要作用[11]。通常實驗室采用流式細胞術(shù)和免疫組化技術(shù)分析腫瘤免疫細胞，但流程相對較復雜，周期較長。CIBERSORT是斯坦福大學研究人員寫出的分析包，創(chuàng)新性采用反卷積算法，可以準確、快速地分析基因表達數(shù)據(jù)，目前得到廣泛應(yīng)用[12]。本研究基于TCGA數(shù)據(jù)庫分析腎癌患者免疫細胞浸潤和臨床相關(guān)信息，對腎癌免疫細胞浸潤模式進行了分析評價，發(fā)現(xiàn)對比正常組織，腎癌組織M2巨噬細胞水平較高；幼稚B細胞、漿細胞、未活化的CD4+記憶T細胞、靜息態(tài)樹突狀細胞水平較低。已有研究表明，M2巨噬細胞與腫瘤增殖、侵襲、轉(zhuǎn)移密切相關(guān)[13]，高水平的M2巨噬細胞是腎乳頭狀癌的核心特征之一，高水平的靜息態(tài)樹突狀細胞與透明腎細胞癌的預(yù)后呈正相關(guān)[14]。在腎癌組織中，CD8+T細胞與濾泡輔助性T細胞、γδT細胞水平呈正相關(guān)(P<0.05)；CD8+T細胞與M2巨噬細胞、未活化的CD4+記憶T細胞水平呈負相關(guān)(P<0.05)，說明CD8+T細胞改變與腎癌免疫微環(huán)境的變化密切相關(guān)，針對243例腎透明細胞癌患者的隊列研究也證實CD8+T細胞對誘導腫瘤免疫和細胞分布具有重要意義[15]，這提示CD8+T細胞可能在免疫微環(huán)境平衡調(diào)控中發(fā)揮重要作用。

通過將腎癌免疫細胞浸潤模式與臨床分期特征結(jié)合分析，發(fā)現(xiàn)隨著腎癌進展，未活化的NK細胞、單核細胞、未活化的肥大細胞、未活化的CD4+記憶T細胞水平逐漸降低；M0巨噬細胞、濾泡輔助性T細胞、CD8+T細胞、活化的CD4+記憶T細胞、漿細胞、調(diào)節(jié)性T細胞水平逐漸增加，提示免疫細胞浸潤模式的改變與腎癌的進展密切相關(guān)。有研究表明，相對于腎透明細胞癌低危患者，高危患者的濾泡輔助性T細胞、調(diào)節(jié)性T細胞、M0巨噬細胞水平明顯升高[16]。M0巨噬細胞主要作用是吞噬清除凋亡細胞及致病物質(zhì)，能夠調(diào)控免疫應(yīng)答，促進組織重構(gòu)、纖維化等[17-18]。單核細胞可以通過吞噬和產(chǎn)生抗體消滅外來入侵病原體，隨著腫瘤進展，機體防御系統(tǒng)功能下降，失去對腫瘤的控制。當效應(yīng)T細胞初次與腫瘤作用后，腫瘤抗原得到清除，T細胞逐步凋亡，但其中有少量的效應(yīng)T細胞轉(zhuǎn)為記憶T細胞，長期存活發(fā)揮效應(yīng)，在遭遇同種抗原攻擊時，快速介導發(fā)生效應(yīng)。伴隨腫瘤的進展，T細胞逐步變?yōu)樾?yīng)T細胞，發(fā)揮腫瘤抗原清除作用，活化的CD4+記憶T細胞增多。調(diào)節(jié)性T細胞在腫瘤浸潤淋巴細胞、外周淋巴細胞和局部淋巴結(jié)中占比較高，與包括腎癌在內(nèi)的多種腫瘤進展密切相關(guān)，可以增強腫瘤抗原識別，加速抗原清除[19]。這些免疫細胞可能參與了腎癌的發(fā)生和進展，為腎癌的診斷和轉(zhuǎn)歸提供了參考。

通過結(jié)合腎癌患者生存狀態(tài)進一步分析，發(fā)現(xiàn)低水平調(diào)節(jié)性T細胞、濾泡輔助性T細胞的腎癌患者5年生存率更高，高水平單核細胞、靜息態(tài)樹突狀細胞、未活化的肥大細胞的腎癌患者5年生存率更高。調(diào)節(jié)性T細胞在維持免疫穩(wěn)態(tài)，介導外周信號上發(fā)揮重要作用，臨床研究表明，外周血調(diào)節(jié)性T細胞的高水平與腎細胞癌患者的不良生存預(yù)后密切相關(guān)[20]。濾泡輔助性T細胞是近些年來發(fā)現(xiàn)的一群獨立的CD4+T細胞亞群，但目前有關(guān)腎癌中濾泡輔助性T細胞的研究較少，仍有待進一步探索[19,21-22]。樹突狀細胞是最強的專職抗原提呈細胞，是CD4+和CD8+T初始細胞的主要激活細胞，然而腫瘤微環(huán)境中的因素會減弱樹突狀細胞的抗原提呈功能，從而限制T細胞活性，促進腫瘤生長[23]。靜息態(tài)樹突狀細胞水平與腎癌5年生存率相關(guān)，對比發(fā)現(xiàn)，其在正常組織中水平較高，在腎癌組織中水平較低，但是正常組織和腫瘤組織中活化的樹突狀細胞水平均不高，可能是腎癌腫瘤微環(huán)境抑制了樹突狀細胞的功能，減少抗原提呈，發(fā)生了免疫逃逸。已有學者研究表明，肥大細胞能夠加速腎癌細胞侵襲和遷徙，促進腎癌惡化[24]。

目前有關(guān)腎癌組織免疫細胞浸潤模式的研究較少，本文通過CIBERSORT算法結(jié)合生物信息學對腎癌免疫細胞浸潤模式進行了系統(tǒng)性分析，腎癌患者發(fā)病年齡整體較高，發(fā)病相對隱匿，如何早期鑒別腎癌、及時預(yù)測腎癌轉(zhuǎn)歸是今后重要的研究方向。相對正常組織，腎癌組織免疫細胞構(gòu)成明顯不同，M2巨噬細胞明顯增多，同時伴隨腫瘤的惡化，腎癌免疫細胞浸潤模式逐步改變，單核細胞、NK細胞、樹突狀細胞等免疫細胞的抗原清除能力逐步降低，腫瘤微環(huán)境的免疫抑制逐步提升。臨床上可以將腎癌免疫細胞浸潤模式作為參照，預(yù)測腎癌患者的生存期及腫瘤進展情況，開展腎癌的個體化免疫治療，改善患者預(yù)后，提高患者生存率，使更多患者獲益。