柯 帥,王文波,廖紅展,彭志柱,邱縣生,唐慧敏,李清華,夏學(xué)巍
(桂林醫(yī)學(xué)院 附屬醫(yī)院 神經(jīng)外科,廣西 桂林 541000)
多型性膠質(zhì)母細(xì)胞瘤(GBM)又稱(chēng)為膠質(zhì)母細(xì)胞瘤,是最常見(jiàn)也最具侵襲性的原發(fā)惡性腦腫瘤。經(jīng)過(guò)標(biāo)準(zhǔn)化治療后的膠質(zhì)母細(xì)胞瘤患者,中位總生存期為14.6個(gè)月,診斷為膠質(zhì)母細(xì)胞瘤就預(yù)示了患者的不良預(yù)后[1]。由于其低存活率、短生存期、高復(fù)發(fā)的特點(diǎn),讓膠質(zhì)母細(xì)胞瘤已經(jīng)成為令人聞風(fēng)喪膽的可怕疾病[2]。傳統(tǒng)的手術(shù)聯(lián)合術(shù)后大劑量放化療的治療方案無(wú)法有效提高膠質(zhì)母細(xì)胞瘤患者的生存時(shí)間,保障患者預(yù)后[3]。因此,針對(duì)膠質(zhì)母細(xì)胞瘤藥物治療靶點(diǎn)的研究成為全世界研究者們臻待解決的問(wèn)題。
1986年, Renato Dulbecco首次提出全基因組測(cè)序的概念,為全世界癌癥研究點(diǎn)亮了一束光[4]。通過(guò)生物信息手段,分析不同個(gè)體基因組間的結(jié)構(gòu)差異,同時(shí)完成SNP及基因組結(jié)構(gòu)注釋?zhuān)梢杂珊暧^到微觀來(lái)認(rèn)識(shí)惡性腫瘤的發(fā)生機(jī)制[5]。Cameron W. Brennan等提出在一個(gè)統(tǒng)計(jì)支持的隊(duì)列中,系統(tǒng)的基因組分析可以定義核心生物學(xué)通路的概念,以此來(lái)促進(jìn)在分子層面上對(duì)膠質(zhì)母細(xì)胞瘤的了解[6]。全面闡明膠質(zhì)母細(xì)胞瘤中的基因變化可以提供新的靶點(diǎn),達(dá)到診斷、預(yù)后預(yù)測(cè)和治療的目的,基因組學(xué)亞型的確定,更可以將膠質(zhì)母細(xì)胞瘤類(lèi)型細(xì)化,進(jìn)行分類(lèi)診療。從最開(kāi)始的基因組學(xué)亞分型研究,到藥物靶點(diǎn)的確定,再到目前藥物通過(guò)血腦屏障、藥物載體的研發(fā),針對(duì)膠質(zhì)母細(xì)胞瘤的研究一直沒(méi)有間斷過(guò)[7]。但是,即便如此,影響膠質(zhì)母細(xì)胞瘤疾病進(jìn)展的因素和作用機(jī)制,仍然是困擾研究者們的難題。
在轉(zhuǎn)錄水平分析膠質(zhì)母細(xì)胞瘤的基因表達(dá)水平,試圖尋找一個(gè)新的標(biāo)記物,有利于從分子層面了解膠質(zhì)母細(xì)胞瘤。數(shù)據(jù)均來(lái)自TCGA(the Cancer Genome Atlas,https://cancergenome.nih.gov/)數(shù)據(jù)庫(kù)。在TCGA數(shù)據(jù)庫(kù)中搜索膠質(zhì)母細(xì)胞瘤mRNA樣本,經(jīng)過(guò)仔細(xì)篩查,從中選取符合研究標(biāo)準(zhǔn)的樣本集,其中共有膠質(zhì)母細(xì)胞瘤169例和正常腦組織5例,隨之進(jìn)行差異基因的篩選。將篩選出來(lái)的差異基因利用GO(Gene Ontology)和KEGG(Kyoto Encyclopedia of Gene and Genomes pathway enrichment analysis)富集分析進(jìn)行功能注解。然后,通過(guò)構(gòu)建差異基因的蛋白互作網(wǎng)絡(luò)篩選出連通性最高的核心基因。最后,利用Kaplan-Meier法對(duì)各個(gè)核心基因進(jìn)行生存分析,從而獲得與膠質(zhì)母細(xì)胞瘤生存時(shí)間相關(guān)的核心基因KCNAB2。
膠質(zhì)母細(xì)胞瘤樣本均來(lái)自TCGA,下載腫瘤樣本和正常對(duì)照樣本的mRNA表達(dá)數(shù)據(jù)。將基因的表達(dá)量轉(zhuǎn)換為數(shù)值,對(duì)重復(fù)的基因取平均值并過(guò)濾表達(dá)量低的基因,然后構(gòu)建設(shè)計(jì)矩陣和對(duì)比矩陣,計(jì)算標(biāo)準(zhǔn)化因子,估計(jì)散度,隨后對(duì)兩組數(shù)據(jù)間進(jìn)行雙馬尾檢驗(yàn),矯正p值,導(dǎo)出基因。根據(jù)所獲得的各個(gè)基因的p值及l(fā)ogFC值篩選差異基因。以上計(jì)算均依靠R軟件(R 3.5.2)進(jìn)行。
GO富集分析是一種常用的基因功能分析方法,通過(guò)將數(shù)量龐大的基因按照生物學(xué)進(jìn)程、分子功能、細(xì)胞組分進(jìn)行聚類(lèi)富集調(diào)查的方法,對(duì)基因功能進(jìn)行注解。KEGG富集分析是一種現(xiàn)在廣泛使用的系統(tǒng)分析基因產(chǎn)物在細(xì)胞中代謝途徑及基因產(chǎn)物功能的分析方法。KEGG可以將基因按照所在的生物學(xué)通路進(jìn)行富集調(diào)查,從而對(duì)目標(biāo)基因做進(jìn)一步注解[8]。利用DAVID (https://david.ncifcrf.gov/tools.jsp)工具[9]對(duì)差異基因進(jìn)行注解和探索。
蛋白互作網(wǎng)絡(luò)(Protein-protein Interaction,PPI)是由單個(gè)蛋白間的相互作用關(guān)系連接而成,蛋白互作網(wǎng)絡(luò)的構(gòu)建可以幫助我們了解疾病狀態(tài)下生物信號(hào)、能量代謝以及蛋白之間的功能聯(lián)系[13]。利用STRING(The search tool for the retrieval of interacting genes,https://string-db.org/)平臺(tái)分析差異基因的蛋白互作信息。將蛋白互作網(wǎng)絡(luò)導(dǎo)入cytoscape軟件,制作蛋白互作網(wǎng)絡(luò)鳥(niǎo)瞰圖。利用cytohubba計(jì)算各差異基因之間連通性得分,得分高的選定為核心基因。鑒于BottleNeck在大數(shù)量基因數(shù)據(jù)計(jì)算時(shí)具有更好的表現(xiàn),并且,BottleNeck算法得到的核心基因(Hub-bottlenecks)更傾向于神經(jīng)信號(hào)轉(zhuǎn)導(dǎo)通路,BottleNeck得分被選用來(lái)選取核心基因[10]。
ONCOMINE(https://www.oncomine.org/)被用來(lái)對(duì)所獲得的核心基因進(jìn)行進(jìn)一步差異分析。ONCOMINE癌基因芯片數(shù)據(jù)挖掘平臺(tái)是目前世界上最大的癌基因芯片數(shù)據(jù)庫(kù)和整合數(shù)據(jù)挖掘平臺(tái),可用于比較主要癌癥類(lèi)型和各自正常組織的差異表達(dá)分析,進(jìn)行一個(gè)或多個(gè)基因的查詢(xún)[11]。 借助已有的研究,綜合回顧、分析核心基因在各個(gè)探針中的表達(dá)情況。
為探索核心基因與膠質(zhì)母細(xì)胞瘤患者生存時(shí)間的聯(lián)系,評(píng)估核心基因?qū)δz熱質(zhì)母細(xì)胞瘤患者生存時(shí)間的影響,利用患者的生存數(shù)據(jù)進(jìn)行生存分析。將核心基因表達(dá)量的改變與膠質(zhì)母細(xì)胞瘤患者的生存時(shí)間和生存幾率的關(guān)系繪制為生存曲線,計(jì)算p值,p值<0.05為具有統(tǒng)計(jì)學(xué)的生存曲線。
從TCGA中下載獲取膠質(zhì)母細(xì)胞瘤轉(zhuǎn)錄組數(shù)據(jù)169例,正常對(duì)照腦組織標(biāo)本5例,基因17 848個(gè)。共獲得差異基因3 183個(gè),上調(diào)基因1 582個(gè),下調(diào)基因1 601個(gè),所獲結(jié)果可視化為熱圖(見(jiàn)圖1)。
表1 GO富集分析結(jié)果Table 1 Results of GO enrichment analysis
*注: BP: Biological process生物學(xué)進(jìn)程; CC:Cellular component細(xì)胞組分; MF:Molecular function分子功能。
表2 KEGG富集分析結(jié)果Table 2 Results of KEGG enrichment analysis
利用STRING得到的蛋白互作網(wǎng)絡(luò),共有1 951個(gè)節(jié)點(diǎn),6 967條邊(見(jiàn)圖2)。將結(jié)果導(dǎo)入cytoscape。利用cytohubba對(duì)差異基因連通性進(jìn)行計(jì)算,將BottleNeck得分前十的差異基因選定為核心基因,得分第一為HIST1H2BH(Histone Cluster 1 H2B Family Member H),BottleNeck得分64分,得分第二為AR(Androgen Receptor),BottleNeck得分54分。KCNAB2(Potassium Voltage-Gated Channel Subfamily A Regulatory Beta Subunit 2)排在第九位,BottleNeck得分26分。在后續(xù)的生存分析中,發(fā)現(xiàn)只有AR與KCNAB2與膠質(zhì)母細(xì)胞瘤患者生存時(shí)間存在統(tǒng)計(jì)學(xué)意義的相關(guān)性,因此將研究重點(diǎn)放在尚未被驗(yàn)證過(guò)的KCNAB2基因上。
圖1 差異基因分布熱圖Fig.1 DEGs distribution heatmap
注:X軸代表TCGA樣本號(hào),Y軸代表基因名。黑色代表非差異基因,紅色代表上調(diào)差異基因,黃色代表下調(diào)差異基因。由于基因數(shù)量巨大,只在Y軸顯示部分不同表達(dá)量的基因。
圖2 差異基因蛋白互作網(wǎng)絡(luò)鳥(niǎo)瞰圖 Fig.2 Bird’s eye view of protein-protein interaction network
注:圖中共有1 951個(gè)節(jié)點(diǎn),6 967條邊。
將KCNAB2基因?qū)隣NCOMINE,各研究中KCNAB2表達(dá)量箱線圖輸出結(jié)果見(jiàn)圖3,結(jié)果提示KCNAB2在各個(gè)研究的不同探針結(jié)果里均存在不同程度的低表達(dá)。綜合三個(gè)已有的研究進(jìn)行meta分析。meta分析的結(jié)果提示,在Bredel Brain2的研究中KCNAB2的低表達(dá)程度最顯著,綜合分析的p值為2.67X10-10。
圖3 ONCOMINE聯(lián)合分析結(jié)果Fig.3 Result of ONCOMINE analysis
Fig.3 Result of ONCOMINE analysis
注:KCNAB2在不同探針下的表達(dá)量,x軸的1代表KCNAB2在正常腦組織中的表達(dá)量,2代表KCNAB2在GBM中的表達(dá)量。(a)中p值1.28X10-10,F(xiàn)C值-7.957;(b)中p值0.016,F(xiàn)C值-2.264;(c)中p值2.67X10-10,F(xiàn)C值-9.866;(d)中p值2.72X10-7,F(xiàn)C值-3.100;(e)中p值7.6 X10-14,F(xiàn)C值-6.271;(f)中p值8.56X10-10,F(xiàn)C值-7.526。
利用Rstudio將168例膠質(zhì)母細(xì)胞瘤樣本的基因表達(dá)信息與生存信息進(jìn)行匹配,繪制核心基因生存曲線,其中只有KCNAB2和AR具有統(tǒng)計(jì)學(xué)意義,結(jié)果可視化輸出見(jiàn)圖4。
圖4 AR和KCNAB2的存活曲線Fig.4 Survival curves of AR and KCNAB2
注:x軸為存活時(shí)間(d),y軸為存活率,紅色曲線代表基因的高表達(dá),綠色曲線代表基因的低表達(dá)。
與膠質(zhì)母細(xì)胞瘤生存相關(guān)的核心基因主要為AR和KCNAB2。AR基因最早在前列腺癌的研究中被發(fā)現(xiàn),在Zalcman,N等的文獻(xiàn)里已經(jīng)證實(shí)其可作為膠質(zhì)母細(xì)胞瘤的一個(gè)重要化學(xué)藥物治療靶點(diǎn),AR的表達(dá)沉默可以導(dǎo)致膠質(zhì)母細(xì)胞瘤的細(xì)胞活力降低;在細(xì)胞系中,AR的沉默會(huì)導(dǎo)致膠質(zhì)母細(xì)胞瘤細(xì)胞在體內(nèi)和體外的凋亡[12-13]。KCNAB2膠質(zhì)母細(xì)胞瘤的表達(dá)差異關(guān)系及對(duì)患者生存時(shí)間的影響均為首次提出,通過(guò)分析KCNAB2的生存曲線可以發(fā)現(xiàn),KCNAB2過(guò)表達(dá)時(shí)膠質(zhì)母細(xì)胞瘤患者表現(xiàn)出更低的生存率。這提示我們,在膠質(zhì)母細(xì)胞瘤中存在著一條KCNAB2相關(guān)的作用機(jī)制,可以影響膠質(zhì)母細(xì)胞瘤患者的生存率。
腦膠質(zhì)瘤來(lái)源于神經(jīng)上皮細(xì)胞,占顱腦腫瘤的40%~50%,是最常見(jiàn)的顱內(nèi)原發(fā)腫瘤。多型性膠質(zhì)母細(xì)胞瘤又稱(chēng)為膠質(zhì)母細(xì)胞瘤(GBM),由星形膠質(zhì)瘤惡變而來(lái),占全部顱腦腫瘤的10.2%[14]。是成人最常見(jiàn)也最具侵襲性的惡性顱腦腫瘤,grade IV級(jí)[15]。膠質(zhì)母細(xì)胞瘤患者中位生存期僅為12~16個(gè)月,5年生存率不足5%[16],診斷為膠質(zhì)母細(xì)胞瘤就預(yù)示了患者的不良預(yù)后[17]。目前,針對(duì)膠質(zhì)母細(xì)胞瘤的方案,主要是手術(shù)聯(lián)合術(shù)后大劑量放療和化療,盡管近年來(lái),手術(shù)、放療和化療的綜合治療手段得到了長(zhǎng)足的發(fā)展,但是膠質(zhì)母細(xì)胞瘤的預(yù)后仍取決于首次手術(shù)切除程度。膠質(zhì)母細(xì)胞瘤具有浸潤(rùn)性生長(zhǎng)的特性,這使得手術(shù)完全切除腫瘤組織成為幾乎不可能的事情[18]。隨著分子生物學(xué)的飛速發(fā)展,針對(duì)膠質(zhì)母細(xì)胞瘤的基因治療被全世界的研究者寄予厚望,基因治療有望在膠質(zhì)母細(xì)胞瘤患者中取得良好的治療效果,甚至最終治愈這一可怕疾病。
篩選得到的十個(gè)核心基因分別是:HIST1H2BH、AR、SNAP25、VAMP2、GNAI1、TP53、CDK1、DYNC1I1,KCNAB2 andSLC04C1。相較于正常腦組織,SNAP25、VAMP2、GNAI1、DYNC1I1、KCNAB2,SLC04C1在膠質(zhì)母細(xì)胞瘤中低表達(dá),而HIST1H2BH、AR、TP53、CDK1在膠質(zhì)母細(xì)胞瘤中高表達(dá)。
離子通道基因KCNAB2,別名AKR6A5、HKvbeta2、HKvbeta2.1、HKvbeta2.2、KCNA2B,KV-BETA-2。編碼電壓門(mén)控K+通道β亞單位蛋白Kvb2,位于1號(hào)染色體1p36區(qū)域。1p36缺失綜合征的患者,在受精卵形成時(shí),由于該染色體片段的缺失,導(dǎo)致部分患者KCNAB2基因的缺失,這一部分缺失KCNAB2基因的患者出現(xiàn)了聯(lián)想學(xué)習(xí)和聯(lián)想記憶功能為主要表現(xiàn)的認(rèn)知功能障礙,通過(guò)研究發(fā)現(xiàn),KCNAB2基因的缺失可以導(dǎo)致杏仁核外側(cè)核投射神經(jīng)元的神經(jīng)生理學(xué)改變,包括動(dòng)作電位的慢后超極化(SAHP)的降低和神經(jīng)元興奮性的增加,從而導(dǎo)致患者聯(lián)想記憶和聯(lián)想學(xué)習(xí)能力受[19-20]。另外,KCNAB2基因缺失的患者中絕大多數(shù)出現(xiàn)了嚴(yán)重癲癇表型,由于離子通道的異??梢栽黾影d癇的易感性,可以認(rèn)為KCNAB2基因作為癲癇表型基因家族的一個(gè)候選基因。但是仍有很大一部分1p36缺失綜合征的患者出現(xiàn)了癲癇表型,卻未發(fā)現(xiàn)存在KCNAB2基因的缺失,提示:KCNAB2基因的缺失是患者發(fā)生嚴(yán)重癲癇的一個(gè)重要危險(xiǎn)因素[21]。除了KCNAB2基因,在膠質(zhì)母細(xì)胞瘤患者中一定還存在一個(gè)潛在的機(jī)制,或其他責(zé)任基因的共同作用,從而導(dǎo)致患者嚴(yán)重癲癇表型的出現(xiàn)。在膠質(zhì)母細(xì)胞瘤的患者中,KCNAB2基因的低表達(dá)是否也會(huì)導(dǎo)致患者出現(xiàn)嚴(yán)重癲癇和認(rèn)知功能障礙有待我們進(jìn)一步的實(shí)驗(yàn)去尋求實(shí)踐支持。
1)KCNAB2基因在膠質(zhì)母細(xì)胞瘤患者中存在顯著的低表達(dá);
2)KCNAB2基因的高表達(dá)與膠質(zhì)母細(xì)胞瘤患者生存時(shí)間的縮短存在相關(guān)性,可能作為膠質(zhì)母細(xì)胞瘤藥物治療的有效靶標(biāo);
3)Kvb2為鉀離子通道蛋白,KCNAB2基因的低表達(dá)可能是導(dǎo)致膠質(zhì)母細(xì)胞瘤細(xì)胞膜電位改變的因素,甚至是導(dǎo)致膠質(zhì)母細(xì)胞瘤癲癇表型的因素之一。