基于多組學(xué)分析血紅素加氧酶1在腫瘤預(yù)后中的作用

田義恒,陳玉清

(1.江西科技師范大學(xué) 藥學(xué)院,江西 南昌 330038;2.河南省人民醫(yī)院 血液科,河南 鄭州 450000)

血紅素加氧酶1(heme oxygenase 1,HMOX1)也稱HO-1,是血紅素分解代謝過程中的限速酶,在全身組織廣泛存在,其可代謝血紅素產(chǎn)生一氧化碳、膽綠素和亞鐵,HMOX1及其酶產(chǎn)物在維持細(xì)胞穩(wěn)態(tài)和調(diào)節(jié)重要的生物過程中發(fā)揮重要作用,包括氧化應(yīng)激、炎癥、細(xì)胞凋亡、細(xì)胞增殖、纖維化和血管生成[1]。HMOX1與多種癌細(xì)胞增殖密切相關(guān),在多種腫瘤和組織中高表達(dá),與腫瘤生長、浸潤、耐藥有關(guān)[2]。正常情況下,HMOX1可以保護(hù)機(jī)體,但對腫瘤組織來說,HMOX1可以增強(qiáng)腫瘤的抗殺傷能力,抑制腫瘤細(xì)胞死亡[3]。使用公共數(shù)據(jù)庫腫瘤基因組圖譜(The Cancer Genome Atlas,TCGA)和基因表達(dá)綜合(Gene Expression Omnibus,GEO)數(shù)據(jù)庫對HMOX1從基因表達(dá)、生存狀態(tài)、基因改變、蛋白磷酸化、免疫浸潤以及細(xì)胞的相關(guān)途徑進(jìn)行基因表達(dá)、預(yù)后生存、免疫浸潤以及相關(guān)蛋白互作分析,有助于從臨床腫瘤樣本的角度理解HMOX1在腫瘤發(fā)生中的作用。

1 資料與方法

1.1 利用TIMER2.0和UALCAN分析HMOX1在各癌癥中的表達(dá)在腫瘤免疫數(shù)據(jù)庫TIMER2.0(http://timer.cistrome.org/)“Gene-DE模塊”中輸入HMOX1,觀察不同腫瘤或特定腫瘤亞型下HMOX1在腫瘤與相鄰正常組織中的表達(dá)差異[4];對于某些無正常組織或者正常組織非常有限的腫瘤,使用GEPIA2(http://gepia.cancer-pku.cn/)數(shù)據(jù)庫,選擇Expression DIY的“Box plot板塊”,設(shè)置P值為0.01,Log2FC為1,以及“匹配TCGA正常和GTEx數(shù)據(jù)”[5]。獲取這些腫瘤組織與基因型和基因表達(dá)量關(guān)聯(lián)(Genotype Tissue Expression,GTEx)數(shù)據(jù)庫中相應(yīng)正常組織表達(dá)差異的箱線圖。UALCAN(http://ualcan.path.uab.edu/)是用于分析癌癥組學(xué)數(shù)據(jù)的交互式網(wǎng)絡(luò)資源,通過使用臨床蛋白質(zhì)組學(xué)腫瘤分析聯(lián)盟(Clinical Proteomic Tumor Analysis Consortium,CPTAC)數(shù)據(jù)集進(jìn)行蛋白表達(dá)分析[6]。選擇乳腺癌、卵巢癌、腎細(xì)胞癌、肺腺癌、子宮內(nèi)膜癌、結(jié)腸癌、頭頸部鱗狀細(xì)胞癌、胰腺癌、多形膠質(zhì)母細(xì)胞瘤、肝細(xì)胞癌等,分析HMOX1在10種原發(fā)腫瘤組織和正常組織蛋白中的表達(dá)情況。

1.2 利用基因表達(dá)譜在線分析平臺(gene expression profiling interactive analysis,GEPIA)2數(shù)據(jù)庫分析HMOX1在各癌癥中的生存預(yù)后使用GEPIA2中“Survival Map”模塊,獲取所有TCGA數(shù)據(jù)庫中表達(dá)HMOX1的腫瘤的總生存期(overall survival,OS)和無病生存期(disease-free survival,DFS)顯著性圖,設(shè)定臨界值-高、臨界值-低分別為50%、50%,初步探索HMOX1表達(dá)與腫瘤患者OS和DFS的關(guān)系。

1.3 利用腫瘤免疫浸潤評估資源(Tumor Immune Estimation Resource,TIMER)2.0數(shù)據(jù)庫分析HMOX1在各癌癥中的免疫浸潤使用TIMER2.0數(shù)據(jù)庫中“Immune”模塊,輸入HMOX1并選擇癌癥相關(guān)成纖維細(xì)胞,研究TCGA不同癌癥類型中HMOX1表達(dá)與癌癥免疫浸潤的關(guān)系。應(yīng)用TIDE、MCPCOUNTER、EPIC、EXCELL算法進(jìn)行免疫浸潤評估。

1.4 利用STRING進(jìn)行HMOX1相關(guān)蛋白互作分析利用STRING網(wǎng)站(https://cn.string-db.org/),輸入蛋白名稱HMOX1,物種選擇“Homo sapiens”,得到HMOX1的蛋白互作用網(wǎng)絡(luò)[7];使用GEPIA2網(wǎng)站獲取與HMOX1表達(dá)相關(guān)的前100個基因,并選取5個基因,使用GEPIA2工具的“Correlation Analysis”模塊分析這5個基因與HMOX1在各癌癥中表達(dá)的相關(guān)性,并利用TIMER2.0數(shù)據(jù)庫的“Gene_Corr”模塊將這 5個靶基因在各個癌癥中的表達(dá)以熱圖形式展示出來。此外,將前100個基因和蛋白互作的基因結(jié)合起來,并進(jìn)行京都基因與基因組百科全書(kyoto Encyclopedia of Genes and Genomes,KEGG)通路分析。將基因列表上傳至DAVID,選擇“OFFICIAL_GENE_SYMBOL”和“Homo sapiens”,獲取功能注釋圖的數(shù)據(jù),使用OmicShare(https://www.omicshare.com/)平臺繪制氣泡圖使KEGG富集通路可視化,使用微生信(http://www.bioinformatics.com.cn/)進(jìn)行基因本體(Gene Ontology,GO)富集分析,使生物學(xué)過程(biological process,BP)、細(xì)胞成分(cell components ,CC)和分子功能(molecular function,MF)數(shù)據(jù)可視化。

2 結(jié)果

2.1 基因表達(dá)分析HMOX1在TCGA數(shù)據(jù)庫多種癌癥中的表達(dá):HMOX1在膀胱尿路上皮癌(bladder cancer,BLCA)、乳腺癌(breast carcinoma,BRCA)、食道癌(esophageal cancer,ESCA)、多形膠質(zhì)母細(xì)胞瘤(glioblastoma multiforme,GBM)、腎透明細(xì)胞癌(clear cell carcinoma of the kidney,KIRC)、腎乳頭狀細(xì)胞癌(papillary cell carcinoma of the kidney,KIRP)、甲狀腺癌(thyroid cancer,THCA)和子宮內(nèi)膜癌(endometrial cancer,UCEC)等癌癥組織中的表達(dá)均高于相應(yīng)正常組織(P<0.05)(圖1A)。

使用GEPIA2數(shù)據(jù)庫將GTEx數(shù)據(jù)集的正常組織作為對照,進(jìn)一步評估HMOX1在正常組織和各癌癥組織中的表達(dá)。HMOX1在GBM、KIRC、急性髓系細(xì)胞白血病(acute myeloid leukemia,AML)、腦低級別膠質(zhì)瘤(lower grade glioma,LGG)、胃癌(stomach cancer,STAD)、睪丸生殖細(xì)胞瘤(testicular cancer,TGCT)和子宮癌肉瘤(uterine carcinosarcoma,UCS)等癌癥組織中的表達(dá)高于相應(yīng)正常組織,并且HMOX1在GBM和KIRC中呈高表達(dá)(P<0.05),與TIMER2.0分析結(jié)果一致(圖1B)。

運(yùn)用UALCAN平臺中CPTAC數(shù)據(jù)集提供的蛋白質(zhì)組數(shù)據(jù)評估HMOX1在腫瘤組織中的表達(dá),結(jié)果顯示與正常組織相比,HMOX1在腎細(xì)胞癌(renal cell cancer,RCC)、UCEC、頭頸鱗狀細(xì)胞癌(head and neck squamous cell carcinoma,HNSC)、胰腺癌(pancreatic cancer,PAAD)和GBM中的總蛋白表達(dá)高于正常組織(圖1C)。

2.2 生存分析使用GEPIA2網(wǎng)站,結(jié)合TCGA和GEO數(shù)據(jù)庫進(jìn)行生存預(yù)后分析,結(jié)果顯示HMOX1高表達(dá)與LGG(P=0.025)和葡萄膜黑色素瘤(uveal melanoma,UVM)(P<0.05)的OS預(yù)后不良相關(guān);DFS分析結(jié)果顯示,高表達(dá)HMOX1與前列腺癌(prostate cancer,PRAD)(P=0.013)、皮膚黑色素瘤(skin cutaneous melanoma,SKCM)(P=0.001 1)和肺腺癌(lung adenocarcinoma,LUAD)(P=0.048)預(yù)后不良相關(guān)。此外,HMOX1低表達(dá)的KIRC(P=0.000 62)與OS預(yù)后不良相關(guān),低表達(dá)HMOX1與膽管癌(bile duct cancer,CHOL)(P=0.007 9)和KIRC(P=0.001 1)的DFS預(yù)后不良相關(guān)。見圖2。

2.3 免疫浸潤分析應(yīng)用TIDE、MCPCOUNTER、EPIC、EXCELL算法探索不同癌癥免疫浸潤水平與HMOX1表達(dá)的相關(guān)性,發(fā)現(xiàn)HMOX1表達(dá)與BLCA、結(jié)腸癌(colon cancer,COAD)、HNSC、LGG、PAAD、PRAD、STAD、TGCT、UCS的癌癥相關(guān)纖維細(xì)胞浸潤水平呈正相關(guān)(P<0.05)。見圖3。

A為HMOX1表達(dá)與BLCA、COAD、HNSC、LGG、PAAD、PRAD、STAD、TGCT、UCS的關(guān)系;B為HMOX1表達(dá)癌癥相關(guān)纖維細(xì)胞浸潤水平呈正相關(guān)(P<0.05)。

2.4 蛋白互作分析為進(jìn)一步研究HMOX1在腫瘤發(fā)生和發(fā)展中的分子機(jī)制,篩選出HMOX1相互作用蛋白以及與HMOX1蛋白表達(dá)的相關(guān)基因,并對HMOX1表達(dá)相關(guān)基因進(jìn)行富集分析。使用STRING數(shù)據(jù)庫獲取10個HMOX1結(jié)合蛋白,使用GEPIA2工具獲取HMOX1表達(dá)相關(guān)前100個基因,結(jié)果顯示HMOX1與LIPA(r=0.41)、DAB2(r=0.52)、ABI3(r=0.39)、CD4(r=0.32)和CD163(r=0.40)基因表達(dá)均呈正相關(guān)。相應(yīng)熱圖數(shù)據(jù)顯示,HMOX1與上述5個基因在大部分癌癥中存在正相關(guān)。將兩組數(shù)據(jù)合并,通過DAVID網(wǎng)站進(jìn)行KEGG和GO富集分析。KEGG數(shù)據(jù)表明,HMOX1在癌癥中發(fā)病機(jī)制可能與破骨細(xì)胞分化、卟啉和葉綠素代謝、金黃色葡萄球菌感染、TNF信號通路有關(guān)。GO富集分析進(jìn)一步表明,這些基因與炎癥反應(yīng)、細(xì)胞黏附有關(guān),并且大多數(shù)基因編碼定位在質(zhì)膜上,促使三磷酸鳥苷和蛋白結(jié)合。見圖4。

A為使用STRING數(shù)據(jù)庫獲取10個HMOX1結(jié)合蛋白;B為熱圖,數(shù)據(jù)顯示HMOX1與LIPA、DAB2、ABI3、CD4和CD163基因在大部分癌癥中呈正相關(guān);C顯示HMOX1在癌癥中發(fā)病機(jī)制可能與破骨細(xì)胞分化、卟啉和葉綠素代謝、金黃色葡萄球菌感染、TNF信號通路有關(guān);D為HMOX1與LIPA(r=0.41)、DAB2(r=0.52)、ABI3(r=0.39)、CD4(r=0.32)和CD163(r=0.40)基因表達(dá)均呈正相關(guān)。

3 討論

腫瘤浸潤性免疫細(xì)胞在腫瘤微環(huán)境中與腫瘤起始、發(fā)展和轉(zhuǎn)移密切相關(guān)[8-9]。研究報道成纖維細(xì)胞促進(jìn)腫瘤細(xì)胞生長、侵襲及免疫抑制,而且基因組穩(wěn)定不易耐藥,有望成為腫瘤免疫治療的新靶標(biāo)[10]。HMOX1是一種催化血紅素在體內(nèi)氧化降解的限速酶,廣泛存在于人體各器官中,其酶活性降低可以導(dǎo)致氧化應(yīng)激降低、炎癥反應(yīng)減弱和細(xì)胞凋亡率降低[11]。另外,HMOX1是血管內(nèi)皮生長因子的有效誘導(dǎo)劑,而血管內(nèi)皮生長因子是引起腫瘤血管新生的一種重要活化因子,可以促進(jìn)內(nèi)皮細(xì)胞增殖,從而使血管通透性增加,使細(xì)胞外基質(zhì)對血管形成的誘導(dǎo)性生物學(xué)作用發(fā)生改變,并且高表達(dá)的HMOX1可以刺激癌細(xì)胞浸潤,導(dǎo)致血管生成、腫瘤轉(zhuǎn)移速度加快、血管通透性增加。在很多癌癥中,HMOX1被認(rèn)為是一種重要的生存因子,研究表明HMOX1的過度表達(dá)可促使腫瘤細(xì)胞增殖、抵抗氧化應(yīng)激、血管生成及病灶轉(zhuǎn)移[12],例如抑制卵巢癌細(xì)胞中HMOX1,可以抑制癌細(xì)胞增殖[13]。雖然有報道表明了腫瘤與HMOX1之間有聯(lián)系,但是HMOX1在不同腫瘤中的發(fā)病機(jī)制是否起到關(guān)鍵作用,通過怎樣的分子機(jī)制起作用還有待回答。

因此,本研究使用TCGA和GEO數(shù)據(jù)庫分析了數(shù)據(jù)庫中不同腫瘤類型,發(fā)現(xiàn)HMOX1在大多數(shù)腫瘤中高表達(dá)。研究表明HMOX1在膀胱癌、卵巢癌、食管鱗狀上皮細(xì)胞癌中的表達(dá)均高于相應(yīng)的正常組織,發(fā)揮促進(jìn)腫瘤增殖和遷移的作用[9,14-15]。但是預(yù)后生存分析表明HMOX1表達(dá)的高低在不同腫瘤中的預(yù)后有一定的差異性。DFS分析結(jié)果顯示,高表達(dá)HMOX1與PRAD、SKCM和LUAD預(yù)后不良相關(guān)。低表達(dá)KIRC與OS預(yù)后不良相關(guān),低表達(dá)HMOX1與CHOL和KIRC的DFS預(yù)后不良相關(guān)。

腫瘤微環(huán)境主要由免疫細(xì)胞、腫瘤相關(guān)成纖維細(xì)胞、炎癥因子、細(xì)胞外基質(zhì)以及癌細(xì)胞自身等共同組成,癌細(xì)胞的存活、生長、侵襲、遷移都受到周圍腫瘤微環(huán)境的影響。癌癥相關(guān)成纖維細(xì)胞作為腫瘤微環(huán)境中的主要成分之一,被認(rèn)為是影響腫瘤微環(huán)境免疫調(diào)節(jié)的主要參與者,可通過多種途徑促進(jìn)腫瘤的免疫逃逸[16],并且癌癥相關(guān)成纖維細(xì)胞在腫瘤的發(fā)展中發(fā)揮著多種作用,例如可以通過分泌促進(jìn)癌細(xì)胞增殖、免疫排斥的生長因子和細(xì)胞外基質(zhì)蛋白發(fā)揮促進(jìn)腫瘤細(xì)胞增殖、轉(zhuǎn)移和侵襲作用,促進(jìn)腫瘤耐藥,調(diào)節(jié)免疫細(xì)胞活性并抑制抗腫瘤免疫應(yīng)答[17]。在原發(fā)性和轉(zhuǎn)移性腫瘤中發(fā)現(xiàn)的癌癥相關(guān)成纖維細(xì)胞是高度通用的、可塑的和有彈性的細(xì)胞,它們通過與腫瘤微環(huán)境中其他細(xì)胞類型的復(fù)雜相互作用積極參與癌癥進(jìn)展[18]。本研究主要分析了HMOX1表達(dá)與癌癥相關(guān)成纖維細(xì)胞浸潤水平之間的相關(guān)性,發(fā)現(xiàn)HMOX1表達(dá)與BLCA、COAD、HNSC、HPV-HNSC、LGG、PAAD、PRAD、STAD、TGCT、UCS的癌癥相關(guān)纖維細(xì)胞浸潤值呈正相關(guān)。

本研究對HMOX1靶向結(jié)合蛋白和HMOX1蛋白表達(dá)相關(guān)基因進(jìn)行了一系列富集分析,結(jié)果顯示HMOX1在癌癥中發(fā)病機(jī)制可能和腫瘤壞死因子(tumor necrosis factor,TNF)信號通路、炎癥反應(yīng)、細(xì)胞黏附有關(guān)。有研究表明TNF信號與塑造腫瘤微環(huán)境有關(guān),并且在癌癥發(fā)展過程中,TNF有一定的抗癌作用,其途徑有誘導(dǎo)細(xì)胞凋亡、影響腫瘤血管系統(tǒng)和增強(qiáng)宿主免疫力3個方面。但是TNF對腫瘤的作用具備雙向性,既有抗瘤作用也具有促進(jìn)腫瘤作用,其濃度非常關(guān)鍵,高濃度TNF-α可選擇性破壞腫瘤血管系統(tǒng)發(fā)揮抗腫瘤效應(yīng),低濃度的TNF通過激活NF-κB信號通路,發(fā)揮抑制細(xì)胞凋亡和促進(jìn)腫瘤增殖作用[19]。通過建立小鼠結(jié)腸癌肝轉(zhuǎn)移模型,研究TNF與HO-1與結(jié)腸癌發(fā)展之間關(guān)系,發(fā)現(xiàn)增加TNF水平,降低HO-1水平可以抑制腫瘤生長[20]。

綜上所述,本研究通過對HMOX1進(jìn)行了全面的泛癌分析,發(fā)現(xiàn)HMOX1表達(dá)與多種人類癌癥的臨床預(yù)后、免疫浸潤等存在相關(guān)性,有助于從臨床腫瘤樣本的角度闡述HMOX1在腫瘤發(fā)生中的作用,為臨床診斷和治療提供依據(jù)與幫助。