王國輝 楊雄濤 朱廣迎

肺癌已經(jīng)成為全球惡性腫瘤發(fā)病率和死亡率的首要原因。據(jù)預測，2018年有180萬肺癌患者死亡，占惡性腫瘤死亡人數(shù)的近1/5（18.4%）[1]。其中，非小細胞肺癌（non-small cell lung cancer, NSCLC）占肺癌的83%，主要包括腺癌、鱗狀細胞癌、大細胞癌。75%的NSCLC患者診斷時已處于中晚期，治療方式主要有手術(shù)、放化療、靶向治療以及免疫治療，但是其5年生存率仍然低于30%[2]。因此，需要進一步研究NSCLC發(fā)生發(fā)展的分子作用機制，以改進其現(xiàn)有的臨床診療現(xiàn)狀。

經(jīng)典的免疫學理論認為，免疫球蛋白G（immunoglobulin G, IgG）是免疫系統(tǒng)在抗原刺激下，由B淋巴細胞或記憶B細胞增殖分化成的漿細胞產(chǎn)生、可以與相應抗原發(fā)生特異性結(jié)合的物質(zhì)，是體液免疫的重要效應分子。然而我國學者邱曉彥教授在國際上首次發(fā)現(xiàn)非B細胞，特別是腫瘤細胞也可以合成和分泌IgG，稱為腫瘤來源免疫球蛋白（cancer-IgG）[3]，并發(fā)現(xiàn)cancer-IgG在腫瘤的發(fā)生發(fā)展中發(fā)揮著重要的生物學作用[4-9]。這一區(qū)別于經(jīng)典免疫學理論的觀點，得到了越來越多國內(nèi)外學者的關(guān)注。Lee等[10]利用卵巢癌細胞系OC-3-VGH的裂解產(chǎn)物去免疫小鼠，得到了一株單克隆抗體RP215，通過識別IgG重鏈恒定區(qū)的一個特殊糖基化位點，可以特異性地識別cancer-IgG。Liao等[9]發(fā)現(xiàn)，RP215識別的cancer-IgG促進了腫瘤細胞的增殖、侵襲和轉(zhuǎn)移，并且可以作為一個潛在的腫瘤干細胞標志物[4]。Tang等的研究也表明，RP215識別的cancer-IgG是通過其特殊的糖基化表位激活黏著斑通路來促進肺鱗癌的發(fā)生發(fā)展。IgG由兩條重鏈和兩條輕鏈構(gòu)成，每條重鏈和輕鏈都由可變區(qū)和恒定區(qū)構(gòu)成。IgG1重鏈恒定區(qū)的編碼基因（immunoglobulin heavy constant gamma 1, IGHG1）的表達量與cancer-IgG的表達量呈正相關(guān)[11]。已有研究[6,7,12-15]表明，IGHG1在膀胱癌、前列腺癌、腎透明細胞癌、胰腺癌、結(jié)直腸癌等處于高表達的狀態(tài)，其與誘導腫瘤細胞凋亡、腫瘤免疫逃逸相關(guān)。

本研究利用生物信息學和免疫組化的方法，分析cancer-IgG在NSCLC中的表達情況與預后生存以及臨床病例特征的關(guān)系。采用基因富集分析方法（gene set enrichment analysis, GSEA），分析IGHG1高表達樣本所富集到的信號通路，初步探討cancer-IgG參與腫瘤發(fā)生發(fā)展的機制。

1 材料與方法

1.1 生物信息學分析在基因表達匯編（gene expression omnibus, GEO） 數(shù)據(jù)庫（http://www.ncbi.nlm.nih.gov/geo）下載NSCLC相關(guān)數(shù)據(jù)集GSE30219[16]、GSE33532和GSE37745[17]的原始數(shù)據(jù)，采用穩(wěn)健多芯片平均標準化（RAM）方法對表達譜數(shù)據(jù)進行標準化處理（圖1）。以IGHG1表達值的中位數(shù)為分組依據(jù)，分為高表達組和低表達組。利用GSE33532數(shù)據(jù)集，比較NSCLC與正常組織中IGHG1表達差異；通過基因表達譜動態(tài)分析在線網(wǎng)站（Gene Expression Profiling Interactive Analysis,GEPIA）[18]分析TCGA和GTEx數(shù)據(jù)庫中IGHG1在LUAD和LUSC中的表達情況；運用GSE30219數(shù)據(jù)集中的生存信息，依據(jù)IGHG1的表達量，繪制生存曲線。利用Kaplan-MeierPlotter數(shù)據(jù)庫（http://kmplot.com/analysis）[19]中的肺癌數(shù)據(jù)集進行在線生存分析。數(shù)據(jù)集GSE37745中的組織樣本依據(jù)IGHG1表達量的中位值分為高表達組和低表達組采用GSEA3.0版本進行基因富集分析（http://software.broadinsitute.org/gsea/index.jsp）[20]。從GSEA網(wǎng)站MsigDB數(shù)據(jù)庫中獲取c2.cp.kegg.v6.2基因集為參照基因，按照默認加權(quán)富集統(tǒng)計方法進行富集分析。每次分析重復1,000次。

1.2 肺癌組織芯片及免疫組織化學染色肺癌組織芯片 購自上海芯超生物科技有限公司，所有患者經(jīng)病理確診為肺癌，收集性別、年齡、病理分級、T分期、淋巴結(jié)轉(zhuǎn)移、遠處轉(zhuǎn)移和美國癌癥聯(lián)合會（American Joint Committee on Cancer, AJCC）第7版TNM分期等臨床病理資料（表1）。

組織芯片置于二甲苯脫蠟20 min，更換新鮮二甲苯重復1次。將脫蠟后的芯片于100%乙醇中浸泡5 min 2次，95%乙醇、80%乙醇、蒸餾水各浸泡5 min。堿性抗原修復液（Tris-EDTA, pH=9）用高壓鍋加熱至沸騰，將芯片放入，計時2 min，自然冷卻至室溫。3%H2O2室溫下避光孵育10 min，用正常羊血清工作液室溫封閉30 min。隨后加入一抗RP215，4 ℃過夜，HRP標記的二抗室溫30 min。DAB顯色，蘇木素復染。RP215單克隆抗體由北京大學醫(yī)學部邱曉彥課題組饋贈，HRP標記的二抗購自于Cell Signaling Technology。

1.3 免疫組織化學染色判定每張切片 隨機選取10個高倍（400×）視野，由兩名病理科醫(yī)生獨立閱片。細胞質(zhì)染色評分利用四個強度等級（0：陰性，1：弱陽性，2：中度陽性，3：強陽性）以及陽性細胞百分比等級（0: 0%, 1: 1%-5%, 2: 6%-25%, 3:26%-50%, 4: 51%-100%）。最終的評分為強度等級和陽性細胞率（百分比）等級的乘積，0-3為低表達，4-9為高表達。

1.4 統(tǒng)計學分析 采用SPSS 20.0版統(tǒng)計學軟件進行數(shù)據(jù)分析。不同組織中cancer-IgG表達水平的比較采用χ2檢驗及獨立樣本t檢驗，與臨床病理特征分析采用χ2檢驗，生存分析采用Kaplan-Meier法并行Log-rank檢驗。P＜0.05定義為差異有統(tǒng)計學意義。在GSEA分析中，以P＜0.05及FDR＜0.25的基因集作為顯著富集的基因集。

2 結(jié)果

2.1 IGHG1在NSCLC中的表達與臨床預后之間的關(guān)系 在GSE33532數(shù)據(jù)集中，IGHG1在NSCLC中的表達水平為（6.635±0.125），正常肺組織的表達水平為（4.965±0.141），NSCLC中的表達量顯著高于正常肺組織，差異具有統(tǒng)計學意義（P＜0.01）（圖2A）。分析TCGA和GTEx數(shù)據(jù)庫中的LUAD和LUSC的數(shù)據(jù)，也得到了同樣的結(jié)論：IGHG1在NSCLC中的表達量顯著高于正常組織（P＜0.01）（圖2B）。分析GSE30219數(shù)據(jù)集中NSCLC患者的生存信息，IGHG1高表達組的總生存期較低表達組患者顯著縮短（P=0.008, HR=1.632, 95%CI: 1.138-2.314）（圖2C）。再利用KM-plotter（http://kmplot.com/analysis）數(shù)據(jù)庫在線分析IGHG1與NSCLC患者預后的關(guān)系，也得出了相似的結(jié)論，高表達IGHG1組的患者預后差于低表達組（P=3.4e-05, HR=1.4, 95%CI:1.19-1.64）（圖2D）。

2.2 Cancer-IgG在NSCLC中的表達與臨床病例特征及預后的相關(guān)性 如圖3（A-D）所示，Cancer-IgG在胞膜、胞質(zhì)和核膜上均有表達，陽性細胞呈巢狀排列，在NSCLC中的陽性率為81.9%（59/72）。如表1所示，依據(jù)cancer-Ig G表達量的高低，進一步對該7 2例NSCLC患者的臨床病例特征進行分析，結(jié)果顯示：cancer-IgG與NSCLC患者的臨床分期（P=0.042）、T分期（P=0.044）和轉(zhuǎn)移（P=0.007）相關(guān)。生存分析顯示，高表達cancer-IgG組中的總生存期顯著低于低表達組（P=0.008,HR=1.746, 95%CI: 1.39-1.94）。

2.3 IGHG1的功能基因富集選用GSE37745數(shù)據(jù)集，運用GSEA分析方法分析IGHG1表達水平對調(diào)控基因集富集的影響。結(jié)果顯示IGHG1高表達的腫瘤樣本可以富集在細胞黏附（NES=1.654、NOMP=0.007、F D Rq=0.0 8 5）、細胞因子-細胞因子相互作用（NES=1.690、NOMP=0.001、FDRq=0.075）和趨化因子信號通路（NES=1.564、NOMP=0.019、FDRq=0.125）（圖4）。

3 討論

目前，經(jīng)典的免疫學理論認為IgG僅由B淋巴細胞合成和分泌，是體液免疫的重要效應分子。IgG是由兩個重鏈和兩個輕鏈構(gòu)成的Y形分子，輕鏈和重鏈之間以二硫鍵連接。每條重鏈和輕鏈都是由恒定區(qū)和可變區(qū)構(gòu)成，可變區(qū)是由若干個基因片段所編碼，包括可變段（V區(qū)）、多樣段（D區(qū)）和連接段（J區(qū)）[21]。B淋巴細胞通過重排V、D和J基因片段，產(chǎn)生IgG的可變區(qū)。由于每段可變區(qū)基因都有不同的拷貝，各段之間還存在不同的組合方式，因此IgG的可變區(qū)可以發(fā)生數(shù)量巨大的變化，與抗原發(fā)生特異性結(jié)合，發(fā)揮抗體的功能[22]。我國學者邱曉彥教授在國際上率先發(fā)現(xiàn)腫瘤細胞也可以產(chǎn)生IgG，其結(jié)構(gòu)與B淋巴細胞分泌的IgG最大的不同就是可變區(qū)呈現(xiàn)相對固定的重排模式。λ鏈的可變區(qū)主要呈現(xiàn)VH5-51/D3-9/JH4、VH3-30/D6-19/JH4兩種重排模式，μ鏈的可變區(qū)主要呈現(xiàn)VH3-15/D3-10/JH4、VH6-1/D6-13/JH4和VH4-30-2/D3-22/JH4三種重排模式[23]。

Cancer-IgG在功能方面與B細胞來源的IgG也有很大的區(qū)別。Liao[4]發(fā)現(xiàn)在乳腺癌細胞中cancer-IgG與腫瘤細胞的侵襲和轉(zhuǎn)移密切相關(guān)，高表達cancer-IgG的細胞具有腫瘤干細胞的特征。Liu[13]發(fā)現(xiàn)IgG1重鏈編碼基因IGHG1在前列腺癌組織中的表達顯著高于正常前列腺組織，并且與前列腺癌的組織學分級相關(guān)。顧江[11]等的研究證明了IGHG1與cancer-IgG的表達量呈正相關(guān)。本研究在GEO、TCGA和GTEx數(shù)據(jù)庫中發(fā)現(xiàn)IGHG1在NSCLC中的表達量顯著高于正常肺組織，并且高表達IGHG1的NSCLC患者預后較差。目前商業(yè)化的抗IgG1抗體是用人外周血提取的IgG做為免疫原獲得的，不能區(qū)分B細胞來源的及腫瘤來源的IgG，我們及前期的報道已表明，用商品化抗人IgG進行免疫組化染色時，腫瘤細胞及間質(zhì)細胞都有陽性反應。RP215抗體是用腫瘤細胞裂解液作為免疫原獲得的，可特異性識別腫瘤來源免疫球蛋白重鏈恒定區(qū)的一個特殊唾液酸化表位，而不識別B細胞來源的IgG[9]。所以RP215對腫瘤來源的IgG更具有特異性。故本研究進一步利用RP215免疫組織化學染色的方法發(fā)現(xiàn)cancer-IgG在NSCLC中的表達水平顯著高于正常肺組織，低表達cancer-IgG患者總生存期顯著長于高表達的患者，分析患者臨床特征發(fā)現(xiàn)cancer-IgG與NSCLC患者的臨床分期、T分期和轉(zhuǎn)移相關(guān)。

表 1 RP215在NSCLC中的表達及臨床特征分析Tab 1 Association between RP215 expression and clinicopathological features of lung cancer patients

GSEA結(jié)果顯示，高表達IGHG1樣本富集到了細胞黏附、細胞因子-細胞因子相互作用和趨化因子信號通路。細胞黏附分子（cell adhesion molecules,CAM）是細胞與細胞、細胞與細胞外基質(zhì)之間黏附作用和信息傳遞的膜蛋白受體，存在于正常細胞中維持組織正常形態(tài)和功能。CAM表達異常，可導致腫瘤細胞間黏附能力降低，侵襲轉(zhuǎn)移能力增強[24]。也有研究[25]表明，CAM通過調(diào)節(jié)血管內(nèi)皮細胞的功能，來促進腫瘤血管的形成，導致腫瘤的復發(fā)和轉(zhuǎn)移。細胞因子-細胞因子相互作用和趨化因子信號通路在調(diào)節(jié)細胞生長分化、細胞間相互作用、免疫調(diào)節(jié)和腫瘤轉(zhuǎn)移方面發(fā)揮重要的作用[26,27]。高表達IGHG1樣本富集到的基因集都參與了腫瘤的侵襲和轉(zhuǎn)移，提示IGHG1可能通過增加腫瘤細胞的侵襲轉(zhuǎn)移能力，來影響NSCLC患者的疾病進程。但是，該分析結(jié)果基于mRNA表達水平，并不能完全代表信號通路蛋白水平的改變，應進一步結(jié)合基礎實驗，驗證該假設。

圖 1 對GEO數(shù)據(jù)進行RAM標準化。A：GSE30219標準化前后對比圖；B：GSE33532標準化前后對比圖；C：GSE37745標準化前后對比圖。Fig 1 Standardization of gene expression.A:The standardization of GSE30219; B: The standardization of GSE33532;C: The standardization of GSE37745.

圖 2 IGHG1在NSCLC中的表達情況及其與預后的關(guān)系。在GEO數(shù)據(jù)集GSE33532（A）和TCGA+GTEx數(shù)據(jù)庫（B）中IGHG1在NSCLC中的表達顯著高于正常肺組織（P＜0.01）。在GSE30219數(shù)據(jù)集（C）和Kaplan-Meier Plotter（D）中的NSCLC數(shù)據(jù)進行生存曲線分析，IGHG1高表達的患者預后顯著差于低表達的患者。Fig 2 IGHG1 expression in NSCLC and its correlation with prognosis.The expression of IGHG1 in NSCLC was significantly higher than that in normal lung tissues (P＜0.01) in GSE33532 (A) and TCGA+GTEx datasets (B); In the GSE30219 (C) and Kaplan-Meier Plotter (D) NSCLC with high expression of IGHG1 had significantly worse prognosis than with low expression.

圖 3 Cancer-IgG免疫組化示意圖及與預后的相關(guān)性（×100）。Cancer-IgG在NSCLC中陰性表達（A）、弱陽性（B）、中度陽性（C）和強陽性（D）表達情況；E：高表達cancer-IgG的NSCLC患者生存期顯著低于低表達的患者。Fig 3 The immunohistochemistry staining results of cancer-IgG in NSCLC and Kaplan-Meier survival analysis.The expression of cancer-IgG in NSCLC with negative (A)、weak positive (B), moderate positive (C) and strong positive (D); E: The survival of patients with high expression of cancer-IgG in NSCLC was significantly lower than patients with low expression.

圖 4 基于GSE37745樣本的基因富集分析。GSEA結(jié)果顯示IGHG1的高表達與細胞黏附分子（A）、細胞因子-細胞因子相互作用（B）和趨化因子通路（C）有關(guān)。Fig 4 Enrichment plots from gene set enrichment analysis.GSEA results showed that high expression of IGHG1 was associated with cell adhesion molecule (A), cytokine-cytokine interaction (B) and chemokine pathway (C).

綜上所述，cancer-IgG可能通過細胞黏附、細胞因子-細胞因子相互作用和趨化因子信號通路影響肺癌細胞的侵襲和轉(zhuǎn)移，并且可以作為評價NSCLC患者預后的指標。