IF50.5｜墨西哥肢端黑色素瘤基因组与转录组图谱

欧洲血统比例越高，BRAF突变越多，且BRAF突变型更像皮肤黑色素瘤。

这篇按生信解读号的读法拆：先看背景和问题，再看作者怎么做、做出了哪几条硬发现，最后把数据和代码收走。

研究背景

肢端黑色素瘤（acral melanoma）是黑色素瘤的一种亚型，好发于手掌、足底和甲下等无毛皮肤，与紫外线暴露无关。在欧美人群中，该亚型占比较低，但在拉丁美洲、非洲和亚洲部分人群中却是最常见的黑色素瘤类型。既往研究主要基于欧洲和亚洲人群，发现肢端黑色素瘤的突变负荷较低，BRAF、NRAS、NF1等经典驱动基因突变频率低于皮肤黑色素瘤，而结构变异和拷贝数改变更为突出。然而，拉丁美洲人群在癌症基因组研究中严重不足，其遗传背景（美洲原住民、欧洲和非洲混合血统）对肿瘤基因组的影响尚不清楚。

科学问题

卡在哪里： 此前肢端黑色素瘤基因组研究多局限于欧洲或亚洲血统人群，拉丁美洲样本几乎空白。由于遗传祖先背景可能影响癌症风险和基因组特征，缺乏该人群数据导致无法评估血统与体细胞突变谱的关联，也限制了对不同祖先背景下肢端黑色素瘤分子机制的全面理解。

本文要回答： 本研究旨在解析墨西哥肢端黑色素瘤患者的基因组和转录组特征，揭示遗传祖先与体细胞突变的关系，并寻找与预后相关的分子分型。

技术路线

作者主要用了：全外显子测序、基因分型、拷贝数分析、RNA-seq、共识聚类、生存分析。

研究团队收集了92例墨西哥肢端黑色素瘤患者的123份肿瘤样本（含原发、转移和复发），对肿瘤和配对正常组织进行全外显子测序，并用基因分型芯片推断患者的遗传祖先成分。随后，通过多种算法识别体细胞单核苷酸变异、插入缺失和拷贝数改变，并利用dNdScv鉴定显著突变基因。对77例原发肿瘤进行RNA测序，通过共识聚类划分转录组亚型，并分析各亚型与临床特征及生存的关联。此外，利用已发表的肢端与皮肤黑色素瘤基因表达数据构建分类器，评估BRAF突变型肢端黑色素瘤的转录组相似性。

图 1　Fig. 1。图1展示了92例患者中突变频率最高的基因（NRAS、KIT、BRAF、NF1等）的突变谱，以及临床特征（分期、溃疡、部位等）的分布。右侧显示BRAF突变与欧洲血统比例的正相关关系，点图直观呈现随着欧洲血统增加，BRAF突变患者比例上升。

核心亮点

亮点 1｜经典驱动基因突变频率低

在92例患者中，NRAS、KIT、BRAF和NF1是最常见的显著突变基因，突变率分别为14%、14%、13%和9%。约40%的肿瘤携带这些经典驱动基因突变，其余60%为“四重野生型”（QWT）。与欧洲血统为主的队列相比，BRAF突变率较低，与亚洲队列相近。此外，女性患者携带驱动基因突变的几率显著高于男性（OR=3.83，P=0.013）。

图 2　Fig. 2。图2呈现了47例样本的拷贝数改变图谱，红色代表扩增，蓝色代表缺失。TERT、CRKL、GAB2、CCND1等区域频繁扩增，CDKN2A/B、ATM、TP53等区域频繁缺失。热图显示不同基因组亚组间的拷贝数模式差异，箱线图比较了不同突变状态和部位的基因组不稳定性评分（GCS）。

亮点 2｜欧洲血统与BRAF突变正相关

通过ADMIXTURE推断遗传祖先，发现患者欧洲血统比例与BRAF突变存在显著正相关（P=0.02），而KIT突变患者倾向于更高的美洲原住民血统。这一结果在控制年龄、性别和总突变负荷后仍然成立，提示遗传背景可能影响肢端黑色素瘤的分子亚型。

图 3　Fig. 3。图3展示了肢端与皮肤黑色素瘤分类器的构建和验证。左图PCA显示两类样本分离，右图显示BRAF突变型肢端黑色素瘤的A:C评分显著低于野生型，更接近皮肤黑色素瘤。下方实验表明BRAF下游信号并不直接激活皮肤黑色素瘤特征基因。

亮点 3｜BRAF突变型转录组更像皮肤黑色素瘤

利用20个基因构建的肢端/皮肤黑色素瘤分类器，发现BRAF激活突变（V600E或L597R）的肢端黑色素瘤转录组评分更接近皮肤黑色素瘤（P=0.045），并在独立队列中得到验证（P=0.039）。进一步分析显示，该差异并非由BRAF下游信号直接驱动，提示这些肿瘤可能起源于不同的黑色素细胞亚群。

图 4　Fig. 4。图4展示了44例原发肿瘤的转录组共识聚类结果。热图显示三个聚类分别高表达表皮/免疫、增殖/色素和氧化磷酸化相关基因模块。箱线图显示Cluster 1有丝分裂指数最低，Cluster 2最高；免疫细胞反卷积显示Cluster 1富集CD4+ T细胞和CAF，Cluster 2富集B细胞。

亮点 4｜三个转录组亚型预后不同

共识聚类将44例原发肿瘤分为三个亚型：Cluster 1高表达表皮/免疫相关基因，Cluster 2高表达增殖和色素相关基因，Cluster 3高表达氧化磷酸化相关基因。Cluster 1患者复发率最低（35.7%），Cluster 2最高（81.2%），且Cluster 2死亡风险显著升高（log-rank P=0.011）。Cluster 1富集CD4+ T细胞和癌症相关成纤维细胞，Cluster 2富集B细胞。

图 5　Fig. 5。图5展示了驱动基因突变和转录组聚类与生存的关联。Kaplan-Meier曲线显示携带驱动基因突变的患者复发率更高，Cluster 2患者的总生存最差。森林图总结了多变量调整后的风险比，表明驱动基因突变和Cluster 2是独立的预后不良因素。

亮点 5｜驱动基因突变与复发和死亡相关

携带任何驱动基因突变的患者复发率显著高于QWT患者（66.7% vs 37.2%，P=0.007），调整混杂因素后OR=5.31（P=0.008）。死亡风险也呈升高趋势（HR=3.19，P=0.1），其中NF1突变患者全部复发。这些结果提示驱动基因状态可能作为肢端黑色素瘤的预后指标。

生信可带走

这一块是给做分析的人用的，不是科普点缀。

• 方法栈： 全外显子测序、基因分型、拷贝数分析、RNA-seq、共识聚类、生存分析
• 公开数据： 原文未给出公开组学登录号
• 代码： 原文未给出公开 GitHub/GitLab 仓库
• 谁该点开原文： 关注黑色素瘤或癌症基因组学、遗传流行病学、肿瘤异质性和精准医学的研究者。

带走一句

做癌症基因组研究时，务必纳入不同遗传祖先背景的样本，否则可能漏掉血统相关的分子亚型，影响精准治疗策略的制定。