黃冠又,劉家駿,杜永貴,蘇 恒,葛學(xué)成,張 欣,甘鴻川

(貴陽(yáng)市第二人民醫(yī)院神經(jīng)外科,貴陽(yáng) 550081)

神經(jīng)膠質(zhì)瘤(glioma)起源于神經(jīng)上皮組織,是一種位于成人原發(fā)性中樞神經(jīng)系統(tǒng)中的惡性腫瘤,具有發(fā)病率高、惡性程度高、致殘致死率和復(fù)發(fā)率高等特點(diǎn),嚴(yán)重威脅患者的生命健康[1-2]。目前膠質(zhì)瘤的治療方案是在手術(shù)切除腫瘤的基礎(chǔ)上,結(jié)合病理和分子病理結(jié)果,輔助放療和(或)化療。但由于腫瘤異質(zhì)性、耐藥性及血腦屏障等因素,導(dǎo)致患者對(duì)各種治療的敏感性差異較大。據(jù)報(bào)道,低級(jí)別膠質(zhì)瘤患者5年生存率為60%～74%[3],而高級(jí)別的膠質(zhì)瘤患者的中位生存時(shí)間僅為15個(gè)月左右,5年生存率低于6%[4]。因此,深入了解膠質(zhì)瘤發(fā)生、發(fā)展的特點(diǎn)和分子機(jī)制,發(fā)現(xiàn)新的分子標(biāo)記物,對(duì)于開(kāi)發(fā)更好的膠質(zhì)瘤治療方法至關(guān)重要。

甲基化CpG結(jié)合蛋白2(methyl-CpG-binding protein 2,MeCP2)是甲基化結(jié)合蛋白家族中的主要成員,是維系大腦功能所必需的,并參與神經(jīng)元的成熟和突觸發(fā)生,在基因轉(zhuǎn)錄調(diào)控中發(fā)揮重要作用[5]。當(dāng)MeCP2基因功能異常時(shí),不能對(duì)相應(yīng)的下游靶基因產(chǎn)生抑制作用,會(huì)造成DNA甲基化改變,進(jìn)而影響DNA的轉(zhuǎn)錄和翻譯,從而使基因的表達(dá)異常,與腫瘤的發(fā)生密切相關(guān),并促進(jìn)腫瘤細(xì)胞的增殖和遷移,如肺癌、結(jié)腸癌、卵巢癌等[6-7]。目前關(guān)于MeCP2與神經(jīng)膠質(zhì)瘤的相關(guān)研究較少,其在膠質(zhì)瘤中的具體作用和機(jī)制仍不明確。

本研究基于對(duì)癌癥基因組圖譜(the cancer genome atlas,TCGA)數(shù)據(jù)庫(kù)進(jìn)行挖掘分析,同時(shí)結(jié)合GTEx數(shù)據(jù)集,分析mRNA在膠質(zhì)瘤(低級(jí)別膠質(zhì)瘤LGG和高級(jí)別膠質(zhì)瘤GBM)和正常腦組織的表達(dá)差異,目的在于分析MeCP2在膠質(zhì)瘤患者中的表達(dá)與臨床病理特征和預(yù)后的關(guān)系,通過(guò)對(duì)樣本數(shù)據(jù)的挖掘,探討MeCP2在膠質(zhì)瘤中的潛在作用,為膠質(zhì)瘤的靶向治療提供理論依據(jù)。

1 材料與方法

1.1 數(shù)據(jù)來(lái)源及臨床資料

從TCGA數(shù)據(jù)庫(kù)(https://portal.gdc.cancer.gov)中下載神經(jīng)膠質(zhì)瘤患者基因表達(dá)譜數(shù)據(jù)(TCGA-GBMLGG、level3、RNA-seq)及臨床和生存信息,共獲得696例膠質(zhì)瘤患者組織和5例癌旁組織mRNA轉(zhuǎn)錄組數(shù)據(jù)。納入標(biāo)準(zhǔn):(1)有完整的患者臨床和基因表達(dá)數(shù)據(jù);(2)病理結(jié)果為L(zhǎng)GG或GBM。對(duì)TCGA數(shù)據(jù)庫(kù)中基因表達(dá)水平的標(biāo)準(zhǔn)化是通過(guò)對(duì)FPKM值以2為底(log2)取對(duì)數(shù)實(shí)現(xiàn)的。本研究獲貴陽(yáng)市第二人民醫(yī)院倫理委員會(huì)批準(zhǔn)(審批號(hào):202153)。

1.2 MeCP2的表達(dá)情況與膠質(zhì)瘤臨床病理和預(yù)后關(guān)系

從UCSC XENA數(shù)據(jù)庫(kù)(https://xenabrowser.net/datapages/s)經(jīng)Toil流程[8]統(tǒng)一處理的TCGA和GTEx的TPM格式的RNAseq數(shù)據(jù),提取TCGA中GBM、LGG(696例)和GTEx中對(duì)應(yīng)的正常組織(1 157例)數(shù)據(jù)。對(duì)數(shù)據(jù)進(jìn)行l(wèi)og2轉(zhuǎn)化后進(jìn)行膠質(zhì)瘤和正常腦組織樣本間的表達(dá)比較。

696例膠質(zhì)瘤臨床數(shù)據(jù)中,男398例,女498例;≤60歲553例,>60歲143例;WHO中樞神經(jīng)系統(tǒng)分類(lèi)中,WHO Ⅱ級(jí)224例,WHO Ⅲ級(jí)243例,WHO Ⅳ級(jí)168例;IDH突變者440例,野生型246例;1p/19q聯(lián)合缺失171例,非聯(lián)合缺失518例。經(jīng)合并后分析GBM和LGG兩組數(shù)據(jù)中MeCP2表達(dá)量的中位值為4.114,將>4.114設(shè)置為高表達(dá)組,≤4.114為低表達(dá)組,繪制Kaplan-Meier生存曲線(xiàn),對(duì)MeCP2不同表達(dá)組進(jìn)行總生存期(overall survival,OS)和疾病特異生存期(disease specific survival,DSS)分析。

1.3 差異表達(dá)基因的篩選和富集分析

應(yīng)用R軟件(4.1.3)limma包進(jìn)行數(shù)據(jù)歸一化和標(biāo)準(zhǔn)化進(jìn)行差異表達(dá)基因篩選,按基因表達(dá)差異倍數(shù)的絕對(duì)值(logFC>2,P<0.05)的標(biāo)準(zhǔn),得到有統(tǒng)計(jì)學(xué)意義的差異基因?；赗軟件使用ggplot2包對(duì)于差異基因進(jìn)行GO功能富集和KEGG信號(hào)通路富集分析。

1.4 MeCP2表達(dá)與免疫浸潤(rùn)和腫瘤微環(huán)境的相關(guān)性

根據(jù)TCGA數(shù)據(jù)庫(kù)中膠質(zhì)瘤的免疫浸潤(rùn)情況進(jìn)行評(píng)估,通過(guò)CIBERSORT工具分析免疫浸潤(rùn)狀況。通過(guò)Wilcoxon秩和檢驗(yàn)比較和分析MeCP2低表達(dá)組和高表達(dá)組之間免疫細(xì)胞浸潤(rùn)差異。每個(gè)樣本的腫瘤免疫微環(huán)境(TME)評(píng)分,包括基質(zhì)評(píng)分、免疫評(píng)分和ESTIMATE評(píng)分,通過(guò)Restimate包來(lái)確定。

1.5 統(tǒng)計(jì)學(xué)處理

使用R軟件(版本4.1.2)進(jìn)行統(tǒng)計(jì)分析。兩組間比較采用Mann-WhitneyU檢驗(yàn),多組見(jiàn)比較采用Kruskal-Wallis秩和檢驗(yàn);log-rank檢驗(yàn)MeCP2與膠質(zhì)瘤患者預(yù)后的關(guān)系,并繪制Kaplan-Meier生存曲線(xiàn);采用單因素和多因素Cox回歸模型分析相關(guān)臨床病理因素和MeCP2表達(dá)與膠質(zhì)瘤患者預(yù)后的關(guān)系;相關(guān)性分析采用Spearman相關(guān)系數(shù)。以P<0.05為差異有統(tǒng)計(jì)學(xué)意義。

2 結(jié) 果

2.1 MeCP2在膠質(zhì)瘤組織中表達(dá)情況

TCGA數(shù)據(jù)庫(kù)及GTEx數(shù)據(jù)集聯(lián)合分析結(jié)果表明,與正常腦組織相比,MeCP2基因在膠質(zhì)瘤組織中明顯上調(diào),兩組間差異有統(tǒng)計(jì)學(xué)意義(P<0.001),見(jiàn)圖1。

***:P<0.001。圖1 MeCP2在膠質(zhì)瘤和正常腦組織中的表達(dá)

2.2 MeCP2的表達(dá)情況與膠質(zhì)瘤臨床病理特征的關(guān)系

對(duì)TCGA數(shù)據(jù)庫(kù)中與膠質(zhì)瘤患者相關(guān)的臨床病理特征數(shù)據(jù)進(jìn)行分析,結(jié)果表明MeCP2表達(dá)與性別無(wú)關(guān),但與年齡、病理分級(jí)和分子分型相關(guān)。MeCP2高表達(dá)與低齡(≤60歲)、低級(jí)別膠質(zhì)瘤(WHO Ⅱ、Ⅲ級(jí))、異檸檬酸脫氫酶(IDH)突變型和1p/19q共缺失相關(guān)(P<0.01或P<0.001),見(jiàn)圖2。

ns:差異無(wú)統(tǒng)計(jì)學(xué)意義。圖2 TCGA數(shù)據(jù)庫(kù)中MeCP2基因表達(dá)與膠質(zhì)瘤患者臨床病理特征關(guān)系

2.3 MeCp2表達(dá)水平對(duì)膠質(zhì)瘤患者生存期的影響

對(duì)696例膠質(zhì)瘤患者進(jìn)行Kaplan-Meier生存分析顯示,MeCP2高表達(dá)組(紅線(xiàn))患者的生存曲線(xiàn)在低表達(dá)組(藍(lán)線(xiàn))之上,即MeCP2低表達(dá)組OS較MeCP2高表達(dá)組明顯縮短,差異有統(tǒng)計(jì)學(xué)意義(P<0.001),見(jiàn)圖3。

圖3 MeCP2在TCGA數(shù)據(jù)庫(kù)中的生存分析

2.4 MeCP2表達(dá)水平與患者預(yù)后的關(guān)系

將患者的MeCP2表達(dá)水平及相關(guān)臨床病理因素納入Cox單因素、多因素風(fēng)險(xiǎn)比例回歸模型分析,結(jié)果顯示高齡(>60歲)、WHO分級(jí)(Ⅲ、Ⅳ級(jí))和IDH突變型是膠質(zhì)瘤患者生存的獨(dú)立影響因素。單因素分析顯示MeCP2的危險(xiǎn)度(HR)=0.412(95%CI:0.319～0.533),見(jiàn)表1。

表1 MeCP2表達(dá)水平與患者預(yù)后的單因素和多因素分析

2.5 MeCP2的差異表達(dá)基因和富集分析

根據(jù)TCGA數(shù)據(jù)庫(kù)中696個(gè)膠質(zhì)瘤樣本,比較樣本中MeCP2高表達(dá)和低表達(dá)組中基因轉(zhuǎn)錄水平差異,進(jìn)行差異分析,篩選條件為基因差異倍數(shù)的絕對(duì)值(logFC>1,P<0.05),共篩選到1 025個(gè)差異表達(dá)基因,表達(dá)水平升高的有266個(gè),表達(dá)水平降低的有759個(gè)。GO分析顯示,這些差異表達(dá)基因參與生物學(xué)過(guò)程包含白細(xì)胞遷移、細(xì)胞外結(jié)構(gòu)/基質(zhì)組織、細(xì)胞趨化等。KEGG信號(hào)通路分析顯示,差異表達(dá)基因參與的分子信號(hào)通路包括細(xì)胞因子受體相互作用、神經(jīng)活性配體受體相互作用、趨化因子信號(hào)通路等。

2.6 MeCP2表達(dá)與膠質(zhì)瘤中免疫細(xì)胞浸潤(rùn)和腫瘤微環(huán)境相關(guān)性

CIBERSORT分析結(jié)果表明,MeCP2基因高表達(dá)組中靜止樹(shù)突狀細(xì)胞、中性粒細(xì)胞和M2巨噬細(xì)胞較多(P<0.05或P<0.001),而低表達(dá)組中記憶B細(xì)胞、初始B細(xì)胞和漿細(xì)胞較多(P<0.01)。與MeCP2高表達(dá)組相比,低表達(dá)組的免疫評(píng)分、基質(zhì)評(píng)分和ESTIMATE評(píng)分都較高(P<0.001),見(jiàn)圖4。

圖4 膠質(zhì)瘤中MeCP2高表達(dá)和低表達(dá)組中免疫相關(guān)評(píng)分

3 討 論

膠質(zhì)瘤的發(fā)生發(fā)展是多基因參與的異常復(fù)雜過(guò)程,一些重要基因的缺失或改變能導(dǎo)致細(xì)胞異常增殖,該過(guò)程可能涉及癌基因EGFR激活和抑癌基因如p53、p16等失活改變或異常表達(dá)。此外,與膠質(zhì)瘤風(fēng)險(xiǎn)增加相關(guān)的基因,如表皮生長(zhǎng)因子受體EGFR基因擴(kuò)增/突變、1p/19q缺失、IDH突變、端粒酶逆轉(zhuǎn)錄酶(TERT)啟動(dòng)子突變等,這些基因突變也是膠質(zhì)瘤發(fā)病機(jī)制的重要因素和判斷預(yù)后的重要標(biāo)志物[9-10]。由于膠質(zhì)瘤高度異質(zhì)性,其致病基因突變復(fù)雜,針對(duì)膠質(zhì)瘤的特異性靶點(diǎn)仍需更多研究,為其個(gè)性化治療提供依據(jù)。

MeCP2基因是一種轉(zhuǎn)錄調(diào)節(jié)因子,能夠調(diào)控基因表達(dá),具有轉(zhuǎn)錄激活或抑制雙重作用。神經(jīng)元中MeCP2缺失會(huì)導(dǎo)致神經(jīng)發(fā)育障礙、Rett綜合征及主要的神經(jīng)精神疾病(例如自閉癥和精神分裂癥)[11]。MeCP2可以在DNA中結(jié)合甲基化的CpG二核苷酸,有研究認(rèn)為MeCP2通過(guò)招募轉(zhuǎn)錄抑制因子Sin3A和組蛋白脫乙?；?HDAC)作為轉(zhuǎn)錄抑制啟動(dòng)子[12]。除了在神經(jīng)元中的作用外,MeCP2還在膠質(zhì)細(xì)胞(星形膠質(zhì)細(xì)胞、少突膠質(zhì)細(xì)胞和微膠質(zhì)細(xì)胞)中發(fā)揮調(diào)節(jié)作用,MeCP2基因調(diào)節(jié)失調(diào)會(huì)損害膠質(zhì)細(xì)胞和神經(jīng)元功能[13]。KIM等[14]研究發(fā)現(xiàn)丙戊酸可以使C6膠質(zhì)瘤細(xì)胞中褪黑激素MT1受體、HDAC 1、HDAC 2、HDAC 3和MeCP2的mRNA表達(dá)出現(xiàn)明顯的時(shí)間依賴(lài)性變化。BIAN等[15]研究揭示MeCP2是腦膠質(zhì)瘤上皮間質(zhì)轉(zhuǎn)化(EMT)的1個(gè)新的調(diào)節(jié)因子,提示抑制MeCP2可能是抑制腦膠質(zhì)瘤EMT的治療選擇。從以上研究報(bào)道可知,MeCP2可能是膠質(zhì)瘤中治療的一個(gè)新的靶點(diǎn),但其在膠質(zhì)瘤中的作用目前鮮有相關(guān)研究。

本研究基于TCGA數(shù)據(jù)庫(kù)和GTEx數(shù)據(jù)集分析MeCP2在膠質(zhì)瘤和正常腦組織間表達(dá)的差異性,結(jié)果顯示MeCP2在膠質(zhì)瘤中表達(dá)水平高于正常腦組織,且表達(dá)水平與患者年齡、病理分級(jí)和分型相關(guān),在低齡(<60歲)、低級(jí)別膠質(zhì)瘤、IDH突變型和1p/19q共缺失中表達(dá)水平明顯增高。通過(guò)Kaplan-Meier生存分析發(fā)現(xiàn),MeCP2的表達(dá)水平可以反映預(yù)后,表現(xiàn)為MeCP2高表達(dá)水平的患者預(yù)后總生存時(shí)間明顯高于MeCP2低水平患者,差異有統(tǒng)計(jì)學(xué)意義?？紤]在上述因素中膠質(zhì)瘤患者預(yù)后較好,該結(jié)果同MeCP2表達(dá)上調(diào)預(yù)示患者預(yù)后較好相吻合。研究表明表達(dá)抑制或缺失是MeCP2在惡性腫瘤中的主要表達(dá)失調(diào)形式[16],本研究中,MeCP2在GBM患者中表達(dá)水平低于LGG,表明其表達(dá)水平隨著膠質(zhì)瘤的發(fā)展有逐步下降趨勢(shì),提示MeCP2的表達(dá)抑制可能參與膠質(zhì)瘤的發(fā)生發(fā)展。但通過(guò)Cox風(fēng)險(xiǎn)比例回歸模型分析,并未證實(shí)MeCP2是影響膠質(zhì)瘤患者生存的獨(dú)立影響因子,提示MeCP2并不能簡(jiǎn)單地歸結(jié)為腫瘤抑制因子,MeCP2高表達(dá)改善膠質(zhì)瘤患者預(yù)后的機(jī)制尚需進(jìn)一步研究證實(shí)。

對(duì)MeCP2高表達(dá)和低表達(dá)組篩選出的差異表達(dá)基因進(jìn)行GO和KEGG富集分析,結(jié)果顯示,差異表達(dá)基因明顯富集在白細(xì)胞遷移和趨化的生物學(xué)過(guò)程及相關(guān)分子信號(hào)通路,如細(xì)胞因子和趨化因子信號(hào)通路等。研究證實(shí)多種趨化因子及受體與膠質(zhì)瘤細(xì)胞的增殖、遷移和侵襲密切相關(guān),并在膠質(zhì)瘤中表達(dá)升高,抑制趨化因子及受體的表達(dá)有助于膠質(zhì)瘤的抗血管生成治療[17-18]。并且一些細(xì)胞因子和趨化因子具有免疫調(diào)節(jié)作用,能影響GBM的腫瘤微環(huán)境,有望成為GBM免疫治療靶向的候選分子[19]。因此,MeCP2可能參與調(diào)控細(xì)胞因子或趨化因子信號(hào)通路并在膠質(zhì)瘤中起作用,為今后膠質(zhì)瘤的免疫靶向治療提供新的靶點(diǎn)。

目前已知,腫瘤微環(huán)境對(duì)免疫治療有重要影響,因此,本研究最后還分析了MeCP2表達(dá)與膠質(zhì)瘤免疫浸潤(rùn)和腫瘤微環(huán)境的關(guān)系。結(jié)果表明,低表達(dá)組有更多的免疫細(xì)胞浸潤(rùn)和更高的免疫評(píng)分,表明低表達(dá)組患者可能對(duì)免疫療法有更好的治療反應(yīng),提示MeCP2可能是膠質(zhì)瘤的一個(gè)免疫治療靶點(diǎn)。

綜上所述,高級(jí)別膠質(zhì)瘤中,MeCP2表達(dá)水平低于低級(jí)別膠質(zhì)瘤,其低表達(dá)水平預(yù)示膠質(zhì)瘤惡性程度高,患者生存期短。MeCP2參與調(diào)控細(xì)胞因子或趨化因子信號(hào)通路,并與免疫細(xì)胞浸潤(rùn)密切相關(guān),MeCP2可能與膠質(zhì)瘤的侵襲和發(fā)展有關(guān),可能是一個(gè)潛在的膠質(zhì)瘤基因治療靶點(diǎn)。但MeCP2在膠質(zhì)瘤中的具體生物學(xué)功能和作用機(jī)制還需進(jìn)一步研究證實(shí)。