基于生物信息学筛选早发性乳腺癌差异表达基因

目的筛选并分析与早发性乳腺癌发生、发展相关的靶基因。 方法(1)在美国国立生物技术信息中心的公共基因芯片数据库(GEO)中检索早发性乳腺癌样本及非早发性乳腺癌样本相关基因芯片数据。对上述数据使用GEO2R、R4.1.2及Venn软件筛选出相关差异表达基因(DEGs)，并运用在线分析工具(Web Gestalt)对DEGs，进行相关功能和信号通路富集分析。(2)同时，通过String在线数据库构建DEGs编码的蛋白质-蛋白质相互作用(PPI)网络，并利用Cytohubba插件对该网络中的基因进行评分，筛选出枢纽基因。将枢纽基因导入Kaplan-Meier生存分析工具(Kaplan-Meier Plotter)，评估枢纽基因在早发性乳腺癌的预后价值。(3)将肿瘤基因组图谱(TCGA)数据库中的肿瘤组织以年龄为标准进行分组，分析枢纽基因在各年龄组中的表达，并与正常组织中的表达进行比较，对得到的枢纽基因进行验证。DEGs表达量的多组间比较使用Kruskal-Wallis H检验，使用Bonferroni法进行两两比较。 结果(1)筛选出编号为GSE89116、GSE109169、GSE36295的基因芯片数据集，共得到80个差异表达基因，其中上调差异表达基因17个，下调差异表达基因63个。富集分析显示：DEGs主要富集在脂质代谢和氧化还原过程以及PPAR信号通路、AMPK信号通路上。(2)在PPI中发现主要的关键基因为PPARG、ADIPOQ、LIPE、PCK1、PDK4、ACACB、PLIN1、CAV1、CD36、ANGPTL4。ACACB、ADIPOQ、CAV1、LIPE、PLIN1、PPARG基因的低表达与乳腺癌患者的不良OS相关(HR＝0.69、0.84、0.76、0.88、0.78、0.82；95%CI：0.59~0.80、0.76~0.93、0.67~0.83、0.79~0.97、0.70~0.86、0.73~0.90；P均<0.050)。(3)ACACB、ADIPOQ、LIPE、PLIN1、CAV1及PPARG这6个与预后相关的基因在正常组织中的表达量均远高于各年龄组肿瘤组织中的表达量(χ²＝104.03、179.57、161.85、189.87、118.56、103.62，P均<0.001)，早发性乳腺癌组(21~40岁)的LIPE、PLIN1表达量低于41~60岁、61~80岁年龄组，差异具有统计学意义(LIPE: Z＝21.07、23.12, P均<0.050; PLIN1：Z＝16.89、18.76, P均<0.050)。 结论早发性乳腺癌与非早发性乳腺癌存在差异基因表达谱，LIPE、PLIN1可能是早发性乳腺癌发生、发展的关键基因。

Screening for differentially expressed genes in early-onset breast cancer using bioinformatic analysis