IF50.5｜组织特异性癌症演化的核心原理

KRAS突变剂量与CDKN2A染色质状态共同决定不同组织中的癌症演化轨迹。

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

研究背景

癌症基因组测序已揭示不同组织间存在显著差异，但塑造这些差异的演化原理仍未被充分理解。KRAS是最常见的突变癌基因，其致癌潜力、遗传互作和表型效应均表现出明显的组织特异性。例如，KRAS在胰腺、肺和肠道中与不同的继发改变协作，且突变KRAS的等位基因失衡模式各异。然而，这些现象背后的机制，如KRAS突变剂量增加的时间、生物学后果及与抑癌基因的互作，大多尚未阐明。小鼠作为最重要的哺乳动物模型，缺乏类似人类CCLE的泛癌种细胞系资源，限制了机制研究。

科学问题

卡在哪里： 此前研究缺乏系统性的小鼠泛癌种细胞系资源，无法在免疫健全背景下研究癌症演化。KRAS突变剂量增加在不同组织中的时间、功能后果及与抑癌基因的互作规律尚不清楚，尤其是CDKN2A等抑癌基因的组织特异性失活模式及其与KRAS剂量增加的时序关系未被揭示。

本文要回答： 建立小鼠泛癌种细胞系资源MCCA，并利用其结合人类队列和功能实验，解析KRAS起始的胰腺、肺和肠道癌症中组织特异性演化的核心原理。

技术路线

作者主要用了：MCCA多组学、免疫表型推断、KRAS剂量分析、CDKN2A染色质状态、转座子筛选。

作者首先构建了包含590个细胞系的MCCA资源，涵盖22个谱系和46种疾病类型，并整合基因组、转录组、表型及临床数据（图1）。随后开发计算方法从测序数据推断品系背景、MHC单倍型和性别，以指导免疫健全移植实验。接着聚焦KRAS起始的胰腺、肺和肠道癌症，通过分析MCCA及人类TCGA/ICGC数据，确定KRAS突变剂量状态及其演化时序。进一步利用可诱导KRASG12D的非转化人类细胞模型和类器官，研究剂量依赖的转录组和表型效应。最后通过表观基因组、单细胞测序和体内转座子筛选，揭示CDKN2A染色质状态如何决定组织特异性选择压力。

图 1　Fig. 1。图1展示MCCA的构建流程和数据集概览。a部分说明细胞系来源和表征步骤，并指向网站。b部分显示590个细胞系的多组学数据热图，按转录组聚类排序，覆盖22个谱系和46种疾病类型。该图强调MCCA的规模和多样性，为后续分析奠定基础。

核心亮点

亮点 1｜MCCA资源覆盖广泛

作者建立了包含590个小鼠癌症细胞系的MCCA资源，涵盖22个谱系和46种疾病类型。每个细胞系均提供基因组、转录组、形态学、组织病理和临床元数据，并通过cBioPortal界面公开。MCCA还包括36个常用公共细胞系。该资源填补了小鼠缺乏泛癌种细胞系集合的空白，为免疫健全移植和机制研究提供了基础。

图 2　Fig. 2。图2展示MCCA免疫表型推断。a部分显示基于SNP的品系背景、MHC单倍型和性别，并举例说明回交后129背景SNP在工程化等位基因附近持续存在。b-e部分展示遗传背景和MHC单倍型分布，以及移植实验验证。该图说明免疫表型注释对免疫健全移植的重要性。

亮点 2｜KRAS剂量增加受正选择

在胰腺、肺和肠道癌症中，KRAS突变等位基因剂量增加（iGD）普遍存在，而剂量减少（dGD）极为罕见，表明iGD受到正选择。KRASMUT-iGD与患者生存期缩短相关。在胰腺中，iGD在早期PanIN阶段即出现；在肺中，iGD在癌而非腺瘤中富集；在肠道中，iGD在腺瘤到癌的进展中被选择。功能实验证实iGD促进肠道腺瘤向癌的转变。

图 3　Fig. 3。图3展示KRAS突变等位基因剂量状态及其演化时序。a-c部分显示小鼠和人类胰腺、肺、肠道癌症中KRASMUT-iGD普遍存在而dGD罕见。d部分显示iGD与生存期缩短相关。e和f部分显示肺和肠道中iGD在癌阶段富集。该图支持KRAS剂量增加受正选择且时序组织特异。

亮点 3｜KRAS效应剂量与组织依赖

利用多西环素滴定表达KRASG12D的非转化人类细胞模型，发现KRASG12D以剂量依赖方式驱动转录组变化，但诱导的分子过程因组织而异：胰腺和肺富集侵袭和发育程序，肠道则富集增殖特征。表型上，胰腺和肺细胞出现失黏附生长，而肠道细胞无此现象。在胰腺中，高剂量KRASG12D诱导去分化为早期祖细胞状态，且Tgfbr2失活可绕过对KRAS剂量增加的需求。

图 4　Fig. 4。图4展示KRASG12D剂量依赖的转录组和表型效应。a部分PCA显示剂量驱动转录组分离。b部分基因集富集显示组织特异性通路。c部分显示胰腺和肺出现失黏附生长而肠道无。d和e部分显示胰腺去分化剂量敏感。f和g部分显示肠道腺瘤中WNT通路激活。该图揭示KRAS效应的组织和剂量依赖性。

亮点 4｜肠道需WNT阻断分化

在KRAS起始的肠道肿瘤中，KRASG12D不能阻断分化，因此KRASG12D-iGD在增生阶段无法克隆扩增。作者发现Apc和Ctnnb1突变在腺瘤阶段出现，且所有腺瘤均显示WNT通路上调。WNT激活阻断分化后，KRASG12D-iGD才获得克隆生长优势。转座子筛选证实WNT通路基因插入在肠道癌中排名靠前，而在胰腺癌中则否。

图 5　Fig. 5。图5展示CDKN2A失活模式及其机制。a和b部分显示胰腺癌中CDKN2A纯合缺失频繁，肺和肠道罕见。c部分显示CDKN2A在胰腺活性染色质、肠道Polycomb抑制、肺二价状态。d部分单细胞数据验证表达差异。e-i部分显示KRAS诱导Cdkn2a表达和衰老的组织差异及转座子筛选结果。j和k部分显示基因改变时序和人类数据验证。该图建立CDKN2A染色质状态与选择压力的因果链。

亮点 5｜CDKN2A染色质状态决定选择压力

CDKN2A失活在KRAS突变胰腺癌中频繁（82%小鼠，64%人），但在肺和肠道中罕见。机制上，CDKN2A在胰腺中处于活性染色质状态并高表达，在肠道中受Polycomb抑制，在肺中呈二价状态。KRASG12D在胰腺中强烈诱导Cdkn2a表达和衰老，但在肠道中仅轻微诱导。体内转座子筛选显示Cdkn2a插入在胰腺癌中排名靠前（44%），而在肠道中仅9%。

亮点 6｜基因改变时序组织特异

在胰腺中，Cdkn2a纯合缺失先于KrasG12D-iGD出现，且KRASMUT-iGD依赖于CDKN2A或TP53失活。在肺和肠道中，KRASMUT-iGD可在CDKN2A完整时发生，Cdkn2a缺失可在其后获得。人类数据验证了这种依赖性：CDKN2APROFTP53PROF胰腺癌中无KRASMUT-iGD，而肺和肠道癌中KRASMUT-iGD与TSG状态无关。

生信可带走

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

• 方法栈： MCCA多组学、免疫表型推断、KRAS剂量分析、CDKN2A染色质状态、转座子筛选
• 公开数据： 原文未给出公开组学登录号
• 代码： 原文未给出公开 GitHub/GitLab 仓库
• 谁该点开原文： 从事癌症基因组学、演化生物学和临床转化的研究人员。

带走一句

分析癌症基因组时，需整合基因剂量、染色质状态和演化时序，才能理解组织特异性改变模式。