食管癌甲基化驱动基因的筛选、功能及调控通路分析

摘要

cqvip:目的筛选食管癌甲基化驱动基因,并进行功能及调控通路分析。方法搜集TCGA数据库中基因表达数据、DNA甲基化数据、临床信息完整的食管癌患者资料162例,记为肿瘤组;其中16例有癌旁对照资料,记为对照组。使用R语言Methyl Mix包综合分析两组患者的基因表达数据和DNA甲基化数据。R语言Methyl Mix包定义肿瘤组和对照组甲基化平均值的差值倍数为差异甲基化值(DM值),Spearman相关分析法分析基因表达与甲基化相关性,以|DM值|>0、|相关系数(Cor)|>0. 3为截断值筛选甲基化驱动基因。使用在线分析软件DAVID 6. 8(https://david. ncifcrf. gov/)对食管癌甲基化驱动基因进行基因本体(GO)功能富集分析。使用在线数据库KOBAS 3. 0(http://kobas. cbi. pku. edu. cn/)对食管癌甲基化驱动基因进行京都基因与基因组百科全书(KEGG)通路分析。结果筛选得到88个食管癌甲基化驱动基因,其中74(84. 1%)个为高甲基化驱动基因,14个为低甲基化驱动基因。最可能发生甲基化的食管癌甲基化驱动基因有Makorin环指蛋白3(MKRN3)基因、人Fermitin家族同系物3(FERMT3)基因、含V-Set免疫球蛋白域蛋白2(VSIG2)基因等。在分子功能层面,食管癌甲基化驱动基因主要影响了金属离子结合、DNA结合、转录因子活性、序列特异性、核酸结合;在生物学过程层面,食管癌甲基化驱动基因影响了转录调控、DNA模板化;在细胞成分层面,食管癌甲基化驱动基因影响了核组分、胞内组分。锌指蛋白家族如ZNF582、ZNF69、ZKSCAN7、ZNF583、ZNF568、ZNF454、ZNF790、ZNF880等和SOX1基因被富集在了多个GO中。食管癌甲基化驱动基因主要参与了乙醛酸和二羧酸代谢通路,甘氨酸、丝氨酸和苏氨酸代谢通路,碳代谢通路及谷胱甘肽代谢通路的调控。结论食管癌甲基化驱动基因有88个,以高甲基化驱动基因为主;最可能发生甲基化调控的食管癌甲基化驱动基因有MKRN3基因、FERMT3基因、VSIG2基因等;食管癌甲基化驱动基因主要影响了金属离子结合、DNA结合等分子功能,转录调控、DNA模板化等生物学过程,核组分、胞内组分等细胞成分;食管癌甲基化驱动基因主要调控代谢相关的多个通路。