《HUMAN MUTATION》:Molecular Subtyping Based on EGFR Mutation-Associated Genes and the Prognostic Role of TRAF2 in Lung Adenocarcinoma
编辑推荐:
本研究旨在通过整合生物信息学方法系统鉴定与EGFR(表皮生长因子受体)突变显著相关的关键基因,并在肺腺癌(LUAD)中建立分子分型模型。研究人员利用来自cBioPortal和癌症基因组图谱(TCGA)LUAD队列的多组学数据集,筛选与EGFR突变状态相关的基因
本研究旨在通过整合生物信息学方法系统鉴定与EGFR(表皮生长因子受体)突变显著相关的关键基因,并在肺腺癌(LUAD)中建立分子分型模型。研究人员利用来自cBioPortal和癌症基因组图谱(TCGA)LUAD队列的多组学数据集,筛选与EGFR突变状态相关的基因,鉴定出18个与EGFR突变频率及患者预后均显著相关的核心基因。基于该基因特征建立了两簇分子亚型分层,这些亚型在总生存期、免疫细胞浸润谱以及免疫治疗干预预测应答方面存在统计学显著分歧。进一步分析包括机器学习算法、多变量Cox回归和分子对接,研究人员确定TRAF2(TNF受体相关因子2)为与EGFR密切相关的重要预后基因。体外实验表明TRAF2促进LUAD细胞增殖、迁移和侵袭。附加分析提示TRAF2可能通过表观遗传调控及相关信号通路参与肿瘤进展。综上,这些发现为EGFR突变LUAD的分子异质性提供了新见解,并为精准预后评估和联合治疗策略提供了潜在靶点。
该研究针对肺腺癌(LUAD)分子异质性强、EGFR(表皮生长因子受体)突变背景下精准预后与免疫治疗获益人群界定困难等问题,在EGFR突变不仅作为驱动事件且伴随复杂共突变与免疫微环境多样性的背景下,开展基于EGFR突变关联基因的系统鉴定与分子分型研究,论文发表于《HUMAN MUTATION》。研究人员整合cBioPortal的“Lung Adenocarcinoma Met Organotropism(MSK, Cancer Cell 2023)”数据集与TCGA-LUAD队列共539例LUAD样本及59例正常样本,通过生物信息学筛选、非负矩阵分解(NMF)共识聚类、免疫去卷积(CIBERSORT、immunedeconv)、TIDE(Tumor Immune Dysfunction and Exclusion)算法、分子对接(CB-Dock2)、机器学习(XGBoost)、多变量Cox回归、GOsemsim功能相似性、GSEA(基因集富集分析)、SMART数据库甲基化分析以及体外细胞功能实验(qRT-PCR、划痕愈合、克隆形成、Transwell),鉴定出18个与EGFR突变及预后相关的核心基因,建立两簇分子亚型,并确定TRAF2(TNF receptor-associated factor 2)为关键预后基因,体外证实TRAF2敲低抑制LUAD细胞增殖、迁移和侵袭;讨论部分总结认为TRAF2通过低甲基化上调,与DNA复制、ATR信号等相关,两簇分型可解析EGFR突变LUAD免疫异质性以指导联合治疗;结论部分指出基于整合生物信息学、分子对接和体外实验确定TRAF2为LUAD中与EGFR相关的重要预后基因并促进恶性表型,18基因分子分型可有效区分具不同免疫特征与免疫治疗潜力的患者亚组,拓展了EGFR中枢调控网络认知并为表观遗传与靶向联合策略提供新路径。
主要关键技术方法:研究人员使用cBioPortal中“Lung Adenocarcinoma Met Organotropism(MSK, Cancer Cell 2023)”数据集分析EGFR突变频率及预后,并纳入TCGA-LUAD队列539例LUAD样本与59例正常样本进行表达与甲基化分析;采用非负矩阵分解(NMF)共识聚类基于18个核心基因表达谱进行分子分型并依据CDF曲线确定最佳簇数;利用R包immunedeconv及CIBERSORT算法量化22种免疫细胞浸润水平,TIDE算法预测免疫检查点阻断应答;借助CB-Dock2平台对EGFR与相关蛋白及药物(吉西他滨、奥希替尼)做分子对接;应用XGBoost、多变量Cox回归、GOsemsim进行基因预后重要性排序;采用GSEA进行功能富集;通过SMART数据库分析TRAF2甲基化探针与表达相关性;体外使用LUAD细胞系A549、NCI-H1299,以qRT-PCR验证siRNA敲低效率,开展划痕愈合、克隆形成、Transwell实验评估增殖、迁移、侵袭表型。
- 3.
Result(结果)
3.1. Screening of EGFR Mutation-Related Genes(EGFR突变相关基因筛选):研究人员在cBioPortal数据集中发现EGFR突变率约35%,EGFR突变组总生存期显著优于非突变组;差异分析筛选阈值|Log2Ratio|>1且p<0.05得到74个突变频率差异基因,最终确定18个在LUAD中同时具差异表达与预后价值的基因,热图与森林图展示其表达与HR值,突变频率在EGFR突变与非突变组间存在差异。
3.2. Molecular Subtyping Based on Key Mutant Genes(基于关键突变基因的分子分型):研究人员对TCGA-LUAD基于18基因表达进行NMF共识聚类,CDF曲线显示在k=5时下降最明显但易致过度分割,综合考虑稳定性与生存差异选定k=2为最优;两簇(Cluster 1、Cluster 2)生存有显著差异(Cluster 1较差);两簇间18基因中除NKX3-1和KRAS外均差异表达显著,而k=3时多数基因无显著簇间差异,故采用两簇分型后续分析。
3.3. Correlation Analysis Between EGFR Mutation-Related Genes and Immunotherapy(EGFR突变相关基因与免疫治疗相关性分析):研究人员用CIBERSORT发现Cluster 2记忆B细胞、静息记忆CD4+T细胞、M2巨噬细胞、静息髓样树突细胞、单核细胞、活化肥大细胞浸润更高,Cluster 1浆B细胞、CD8+T细胞、活化记忆CD4+T细胞、滤泡辅助T细胞、活化NK细胞、M0/M1巨噬细胞更高;Cluster 2高表达CD274、CTLA4、HAVCR2、PDCD1、PDCD1LG2、SIGLEC15、ITPRIPL1、TIGIT等免疫检查点基因;TIDE评分Cluster 1显著高于Cluster 2,提示Cluster 2可能对免疫检查点阻断疗效更好;KEGG富集显示Cluster 1富集细胞周期、DNA复制、谷胱甘肽代谢、氨基酸生物合成、p53信号通路,Cluster 2富集造血细胞谱系、细胞黏附分子、细胞因子–受体相互作用、IgA产生肠免疫网络、Th17分化;临床分期分布显示两簇在T分期、N分期有显著差异。
3.4. Correlation Analysis Between EGFR Mutation-Related Genes and EGFR(EGFR突变相关基因与EGFR相关性分析):研究人员在TCGA-LUAD中发现EGFR表达与KIT、TRAF2负相关,与其余10个基因正相关;STRING互作网络显示HGF、TRAF2、KRAS、PIK3CG、PAK1、BRCA1与EGFR有功能/物理互作;分子对接证实这六个蛋白均与EGFR具强结合潜力。
3.5. TRAF2 Is Identified as a Key Prognostic Gene in LUAD(TRAF2被确定为LUAD关键预后基因):研究人员通过GOsemsim功能相似性、XGBoost预后重要性排序、多变量Cox回归综合判定TRAF2为关键预后基因;虽GOsemsim与XGBoost未排第一,但Cox与分子对接(图4E最高对接分数)支持TRAF2与EGFR互作且预后显著;TRAF2高表达关联年龄较大、男性、较差总生存期、吸烟史;TRAF2表达与17种免疫细胞浸润显著相关;功能分析显示TRAF2参与凋亡、肿瘤增殖、DNA修复、p53信号通路。
3.6. Expression and Functional Analysis of TRAF2(TRAF2表达与功能分析):研究人员证实TRAF2在LUAD组织(配对及非配对)显著上调,吸烟患者高于非吸烟;ROC曲线显示良好诊断价值;HPA数据库蛋白水平LUAD高于正常;GSEA显示TRAF2关联DNA甲基化、ATR信号通路、DNA复制;分子对接表明TRAF2与吉西他滨、奥希替尼具强亲和力。
3.7. Analysis of TRAF2 Methylation Levels(TRAF2甲基化水平分析):研究人员发现LUAD中TRAF2 DNA甲基化显著低于正常,且男女之间、吸烟与非吸烟之间有差异;SMART数据库定位TRAF2相关甲基化探针染色体分布与基因组信息,跨肿瘤类型显示LUAD甲基化较低;16个甲基化探针与TRAF2表达负相关,提示低甲基化可能上调TRAF2参与致癌。
3.8. Knockdown of TRAF2 Inhibits LUAD Cell Proliferation and Metastasis(TRAF2敲低抑制LUAD细胞增殖与转移):研究人员用siTRAF2#1、siTRAF2#2高效敲低TRAF2(qRT-PCR验证),在A549与NCI-H1299中划痕愈合实验显示迁移能力下降;克隆形成实验显示增殖受抑;Transwell(有/无Matrigel)显示迁移与侵袭能力均显著降低。
讨论部分总结:研究人员指出LUAD分子异质性强致治疗耐药与预后差,单一EGFR突变检测不足;整合多组学鉴定18个EGFR突变相关核心基因并建立两簇分型,Cluster 2高免疫检查点表达、低TIDE评分提示更好免疫阻断应答潜力,挑战了EGFR突变LUAD均为“免疫荒漠”的观念;TRAF2经多变量Cox、功能相似性、XGBoost及分子对接确定为与EGFR互作的关键预后基因,表达与EGFR负相关,体外敲低抑制恶性表型;GSEA显示TRAF2关联DNA甲基化、ATR信号、DNA复制,提示其促基因组不稳定;甲基化分析显示LUAD中TRAF2低甲基化、探针负相关,与环境及宿主因素(性别、吸烟)相关,表观上调连接外部暴露与EGFR驱动网络;分子对接显示TRAF2与吉西他滨、奥希替尼强结合,或为药效学生物标志物;局限包括仅用A549(KRAS突变)、NCI-H1299(无EGFR蛋白)两株细胞,需在EGFR经典突变细胞(PC9、HCC827)及类器官、动物模型验证TRAF2对EGFR-TKI敏感性影响;TRAF2与EGFR互作机制(直接泛素化、支架、转录调控)待阐明;18基因分型与两簇需独立免疫治疗队列前瞻验证;未来应开发选择性TRAF2抑制剂并在小鼠模型评估,探索TRAF2预测奥希替尼或化疗耐药的价值。结论部分翻译:综上所述,通过整合生物信息学、分子对接和体外实验,研究人员确定TRAF2为LUAD中与EGFR相关的重要预后基因,可促进肿瘤增殖、迁移和侵袭;建立的18基因分子分型有效区分具不同免疫谱与免疫治疗应答潜力的患者亚组;这些发现拓展了LUAD中EGFR中枢调控网络的认知,为表观遗传与靶向联合治疗提供了新途径。