引言

同样的致癌基因变异，为什么在不同器官中会引发截然不同的后果？为什么有的组织迅速陷入癌变的深渊，而有的组织却能长期保持沉默？

一直以来，癌症基因组测序为我们描绘了不同癌症类型的突变图谱，但隐藏在这些突变背后的演化法则却始终晦暗不明。

2月25日，《Nature》的研究报道“A disease model resource reveals core principles of tissue-specific cancer evolution”，为我们揭开了这一谜题的一角。

该研究通过构建一个庞大的疾病模型库，不仅以前所未有的分辨率重构了癌症演化的时间线，更向我们展示了一个令人深思的现象：癌症的演化并非完全随机的“掷骰子”，而是受到组织原有的表观遗传状态、细胞发育层级以及突变基因剂量效应的严格制约。这场发生在细胞深处的微观博弈，其底层逻辑比我们想象的更为严谨和确定。

构建全景视角：五百九十个微观生命的拼图

要探究不同组织中癌症演化的普遍规律，仅靠单一的模型或零散的细胞系是远远不够的。研究人员开启了一项浩大的工程，构建了包含五百九十个细胞系的“小鼠癌症细胞系图谱” (Mouse Cancer Cell line Atlas, MCCA)。这个图谱涵盖了二十二个细胞谱系和四十六种疾病亚型，包含了八十一个不同的工程化等位基因 (Alleles) 或外源性致癌触发因素。

这绝不仅仅是一个简单的细胞收集库。研究人员对这些模型进行了多维度的深度测序与刻画，包括基因组学 (Genomics)、转录组学 (Transcriptomics) 以及详尽的临床元数据 (Clinical metadata)。更为巧妙的是，为了让这些细胞系能够在免疫功能健全的小鼠宿主中发挥最大价值，研究人员开发了一套基于基因组测序推断免疫表型 (Immunophenotyping) 的计算工具。

通过提取来自小鼠基因组计划中二十九个近交品系的21,923,209个单核苷酸多态性 (Single-nucleotide polymorphisms, SNPs) 数据，他们计算了品系间SNP模式的皮尔逊相关性，最终定义了十五个系统发育相关的品系群。这使得他们能够精准锁定1,097,314个“特征性SNP” (Signature SNPs)。与此同时，为了精准界定主要组织相容性复合体 (Major histocompatibility complex, MHC) 单倍型，他们将MHC基因座划分为六个基因簇（H2-K、A、E、D、Q、T），并利用375,097个MHC相关的SNP，最终确定了44,219个特征性SNP用于单倍型分类。

这些庞大的数字背后，是极高的精准度。数据表明，在进行MHC匹配的移植实验中，当供体与受体之间的SNP错配率高达42%和45%时，肿瘤细胞无法定植生长；而当错配率低于12%时，肿瘤细胞则能实现稳健的定植。

更引人深思的是，这种遗传背景的差异被巧妙地转化为探索肿瘤免疫的工具。例如，在完全匹配的受体（如C57BL/6;129-F1杂交小鼠）中，仅有肿瘤的体细胞突变会产生免疫原性，其蛋白质改变的肿瘤突变负荷 (Protein-altering tumor mutational burden, pTMB) 仅为0.4个/Mb。然而，如果将相同的细胞系移植到C57BL/6受体中，那些源自129品系的特异性生殖细胞系变异就会作为新抗原暴露出来，使得“有效”pTMB跃升至1.4个/Mb，这与人类高突变负荷的癌症特征高度吻合。这种对模型遗传背景的精准掌控，为后续的机制探索奠定了坚实的基础。

剂量的权力：突变等位基因失衡的隐秘时间线

在癌症基因组学中，致癌基因（如KRAS）的突变常常被视为一个“非黑即白”的二元事件。然而，致癌信号的真实强度，往往取决于突变基因的剂量。研究人员通过整合单核苷酸变异 (Single-nucleotide variant, SNV) 和拷贝数变异 (Copy-number variation, CNV) 数据，将KRAS突变 (KRASMUT) 的等位基因状态精确划分为三种：基因剂量降低 (dGD)、杂合状态 (HET)和基因剂量增加 (iGD)。

数据呈现出高度的一致性：在所有被分析的组织（胰腺、肺、肠道）中，无论是在小鼠细胞系还是人类癌症队列中，dGD的发生率几乎为0%。相反，KRASMUT-iGD极为普遍。这意味着，增加突变基因的剂量在演化过程中受到了强烈的正向选择。更重要的是，在人类癌症队列（如ICGC-PanCuRx胰腺癌队列、TCGA-LUAD肺癌队列和TCGA-COAD肠癌队列）中，携带KRASMUT-iGD的患者，其疾病特异性生存率均出现了显著下降（胰腺癌P=0.0207，肺癌P=0.0065，肠癌P=0.0284）。

然而，真正拉开组织特异性差异序幕的，是KRASMUT-iGD获得的时间节点。

在胰腺癌模型中，研究人员通过激光显微切割 (Laser microdissection) 技术提取了极微小的胰腺上皮内瘤变 (PanIN) 细胞，并利用TaqMan qPCR技术减去基质细胞污染的干扰，最终还原了纯净的肿瘤细胞突变等位基因频率 (Variant allele frequency, VAF)。结果显示，在极早期的PanIN病变阶段，KRASMUT剂量增加就已经发生。

相反，在肺部，针对Trp53野生型小鼠的早期肺腺瘤和晚期肺癌进行的对比分析表明，KRASMUT-iGD在早期的腺瘤中极少出现，而是在进展为癌症阶段才被频繁观察到。

肠道的情况则更为特殊。通过对从极早期增生、腺瘤到晚期癌症及转移灶的六十三个类器官 (Organoids) 进行时序分析，研究人员发现，在增生和腺瘤阶段，几乎检测不到KRASMUT-iGD；只有在演化至癌症及转移阶段时，KRASMUT-iGD的比例才大幅攀升。功能性移植实验进一步证实了这一点：当把携带不同亚克隆KRASMUT-iGD比例的肠道腺瘤类器官移植到小鼠体内时，它们的定植能力与移植时的KRASMUT-iGD水平呈强烈的正相关，而与整体的SNV或CNV负荷无关。

这引出了一个关键问题：为什么同为KRAS突变，胰腺急于在早期就增加它的剂量，而肺和肠道却显得“不紧不慢”？

逆转生命时钟：去分化风暴的引擎

为了探寻背后的机制，研究人员利用多西环素 (Doxycycline) 可调控的慢病毒载体系统，在未发生转化的正常人类胰腺、肺和肠道细胞中，引入了不同表达水平的KRASG12D。为了最真实地模拟体内环境，他们将细胞置于三维 (3D) 基质胶中进行培养。

转录组主成分分析 (Principal component analysis, PCA) 揭示了一个惊人的事实：随着多西环素浓度的增加（即KRASG12D表达量上升），所有细胞系的转录组变化都在主成分1 (PC1) 上发生了显著位移。PC1轴捕获了高达64%到89%的转录组方差，这说明基因剂量的改变是驱动细胞内部状态重塑的最核心力量。

但是，沿着PC1轴，不同组织走向了不同的命运。基因集富集分析 (Gene set enrichment analysis, GSEA) 显示，在胰腺和肺细胞中，高剂量的KRASG12D强烈激活了与侵袭（上皮间质转化 EMT、粘附斑）和胎儿程序再激活（内胚层分化/形成）相关的基因表达谱。在显微镜下观察3D球状体，也能直观地看到细胞从紧密的“粘附”状态转变为松散侵袭的“非粘附”状态。相反，在肠道细胞中，高剂量KRAS触发的主要是细胞周期和DNA复制等增殖信号，并未出现显著的组织解离现象。

对于胰腺而言，这种“去分化” (De-differentiation) 尤为关键。胰腺癌起源于高度分化的腺泡细胞或导管细胞。研究人员将人类胰腺导管上皮 (HPDE) 细胞在不同KRASG12D剂量下的转录组，与体外干细胞向导管细胞分化序列的特定时间点基因特征进行了比对。结果显示，只有在“高剂量”的KRASG12D作用下，细胞才会彻底抛弃成熟导管细胞的特征（第30、45、59天标记），转而表现出早期胰腺祖细胞的特征（第3、13、20天标记）。

在小鼠腺泡细胞模型 (266-6) 中，同样观察到了剂量依赖性的去分化：随着突变剂量的增加，腺泡细胞成熟标志物（如Cpa2、Prss1、Ptf1a）表达骤降，而导管前体标志物（如Sox9）则急剧攀升。

更令人惊叹的是，细胞并非仅仅依靠被动的基因组扩增来增加KRAS剂量。在腺泡细胞向导管细胞化生 (Acinar-to-ductal metaplasia, ADM) 的离体实验中，当正常的腺泡细胞在体外自发转分化为导管样细胞的24小时内，即使没有外源性突变，其内源性Kras的表达量也飙升了5.7倍。这意味着，在胰腺演化的极早期，细胞状态的转换（去分化）本身就会通过非遗传机制极大地放大内源性KRAS的信号。

这就是为什么胰腺癌在早期阶段就需要KRASMUT-iGD的原因——它是启动去分化风暴，打破终末分化细胞身份锁链的必要引擎。如果在成年小鼠的腺泡细胞中，利用scAAV8载体递送sgRNA敲除维持腺泡细胞特征的Tgfbr2基因，人为诱导去分化，那么细胞对早期获得KRASMUT-iGD的选择压力就会大幅降低。

增殖与分化的博弈：肠道突围的特殊入场券

与胰腺的宁静不同，肠道上皮是哺乳动物体内更新最快的组织之一。这里的细胞处于持续的增殖和剥落之中。在这样的环境中，单凭一个KRAS突变，虽然能引起细胞的过度增殖（形成增生），但突变细胞很快就会随着上皮的快速脱落而被冲刷出肠道。

既然KRAS剂量增加无法在肠道中触发类似胰腺的“去分化”来锁定细胞状态，那么肠道肿瘤是如何获得立足之地的？

答案是：阻断分化 (Block of differentiation)。

研究人员对处于不同演化阶段的肠道类器官进行了基因组测序，结果发现，在从增生走向腺瘤的关键节点，有48%的样本突然获得了WNT通路关键基因（Apc或Ctnnb1）的体细胞突变。不仅如此，转录组数据显示，无论是否携带可检测到的突变，所有的腺瘤组织都表现出WNT信号通路的强烈上调。

这一发现揭示了肠道癌症演化中一个极其重要的顺序法则：在KRAS突变引发的锯齿状肿瘤发生路径中，必须首先通过WNT通路的异常激活来阻断肠道干细胞/前体细胞的正常分化路径。只有当分化被阻断，细胞免于被推向肠道绒毛顶端并脱落的命运后，KRASMUT-iGD所赋予的竞争性生长优势，才能真正在癌变阶段转化为不可逆的克隆扩张。

这解释了为什么在肠道中，KRASMUT-iGD总是姗姗来迟——在没有获得这特殊的“入场券”（WNT通路异常）之前，任何试图提前增加KRAS剂量的细胞，都不过是肠道快速更迭中的匆匆过客。

染色质的记忆：抑癌基因的枷锁与反击

在致癌基因疯狂踩下油门的同时，细胞内的抑癌基因 (Tumor suppressor genes, TSGs) 绝不会坐以待毙。其中，CDKN2A作为应对异常增殖信号的核心刹车系统，其与KRAS的博弈成为了决定演化走向的关键。

研究人员在微切割的小鼠肿瘤组织中发现了一个极其悬殊的比例：在KRAS突变的胰腺癌中，Cdkn2a的双等位基因纯合缺失 (Homozygous loss) 发生率高达82%；而在肺癌和肠癌中，这一比例仅为5%和11%。人类癌症队列的数据也印证了这一鸿沟：人类胰腺癌中CDKN2A失活率为64%，肺癌为15%，肠癌仅为2%。

为什么同样是面对KRAS的失控，不同组织对CDKN2A的依赖程度和抛弃意愿差异如此巨大？

线索深藏在表观基因组 (Epigenome) 的折叠方式中。通过挖掘ROADMAP和ENCODE数据库中的表观遗传数据，研究人员比较了健康人类胰腺、肺和肠道组织中CDKN2A基因座位的染色质状态。

在胰腺中，CDKN2A的转录起始位点被高水平的组蛋白H3K4me3修饰占据。这是一种代表染色质处于“开放”和“活跃转录”状态的表观遗传标记。此时的CDKN2A就像一位全副武装、时刻保持警惕的守卫。

而在肠道中，情况则完全相反。肠道细胞的CDKN2A基因座位上，覆盖着致密的H3K27me3修饰。这是一种由多梳抑制复合物 (Polycomb repressive complexes, PRCs) 催化的抑制性标记，它将CDKN2A紧紧锁死在沉默状态。肺部组织则处于一种中间地带，表现为活跃与抑制标记并存的“二价” (Bivalent) 状态。

通过分析单细胞RNA测序 (scRNA-seq) 数据，研究人员证实，作为癌症假定细胞来源的胰腺腺泡细胞和导管细胞中，CDKN2A本底表达明显；而在肺肺泡细胞和肠道干细胞中，其表达极低甚至缺乏。

这种原生的染色质记忆，直接决定了细胞在遭遇KRAS突变时的反应机制。

在体外实验中，当利用小分子抑制剂（A-395和UNC1999）同时抑制PRC2复合物的活性，阻断H3K27me3的重塑时，肠道类器官中被锁死的Cdkn2a被瞬间唤醒，表达量激增，导致细胞陷入严重的生长停滞。如果提前在肠道细胞中敲除Cdkn2a基因，这种由PRC2抑制剂引起的生长停滞就能被完全挽救。而在胰腺细胞中，由于Cdkn2a本来就处于开放状态，PRC2抑制剂并没有引发明显的额外效应。

回到体内的演化现场，表观遗传的差异带来了截然不同的表型。在携带KRASG12D的小鼠模型中，研究人员进行了衰老相关的β-半乳糖苷酶 (SA-βGal) 染色。结果显示，胰腺的早期病变（ADM和PanIN）区域呈现出大面积的强阳性染色，这表明开放的CDKN2A敏锐地捕捉到了异常的KRAS信号，并成功触发了强大的细胞衰老（一种强力的抑癌机制）。相反，在肺部和肠道的早期病变中，衰老细胞几乎难觅踪影，取而代之的是高度活跃的细胞增殖（Ki-67染色强阳性）。

为了量化这种不同组织中失活Cdkn2a所带来的“演化收益”，研究人员利用piggyBac转座子系统在小鼠体内进行了正向遗传学筛选。转座子可以在基因组中随机插入，如果它破坏了某个抑癌基因，获得了生存优势的细胞就会疯狂扩增。通过对肿瘤样本进行插入位点深度测序并排序，数据显示，在胰腺癌中，高达44%的Cdkn2a转座子插入在所有突变中排名前十，这意味着破坏Cdkn2a赋予了胰腺细胞极大的生存竞争优势；而在肠癌中，这个数字仅为9%（取而代之的是WNT通路基因的插入频频上榜）。

胰腺细胞因为具备活跃的CDKN2A，能够成功对早期致癌信号进行拦截；但这同时也意味着，肿瘤如果想要继续发展，就面临着极高的选择压力去破坏这个防线。而肠道细胞本就因频繁分裂的生理需求而关闭了CDKN2A，因此对失去它并不敏感。

演化的定数：基因改变的先来后到

当我们将基因剂量效应（KRASMUT-iGD）与抑癌基因的拦截机制（CDKN2AHOM，即双等位基因纯合丢失）结合在一起考察时，一个关于癌症演化时序的清晰法则浮出水面。

研究人员对小鼠的原发灶和转移灶进行了细致的系统发育分析，通过比较配对样本中Cdkn2a和Kras的拷贝数变异 (CNV) 和杂合性缺失 (LOH) 模式，重构了时间线。

在胰腺中，所有的证据都指向了一个不可逾越的次序：Cdkn2a的纯合失活，必然发生在KRASMUT-iGD之前。这就像是必须先拆除刹车系统，才能猛踩油门。事实上，在保留了Cdkn2a功能的小鼠PanIN病变中，KRASMUT-iGD的发生率极低（3.8%）；而一旦敲除Cdkn2a，PanIN中KRASMUT-iGD的发生率迅速飙升至30.8%。在分析人类胰腺癌队列（结合COMPASS和TCGA-PAAD）时，这一铁律同样适用：在那些CDKN2A和TP53均保持完好的胰腺癌样本中，KRASMUT-iGD的发生率为0（0/16）。没有拆除路障，KRAS的扩增就不会被允许。

然而，在肺和肠道中，时序被颠倒了。系统发育树显示，肺癌和肠癌可以在保留部分或全部Cdkn2a功能的情况下，率先获得KRASG12D-iGD。Cdkn2a的彻底丢失（如果发生的话），往往出现在演化的晚期。为什么肺癌晚期仍然会选择丢失CDKN2A？对人类肺腺癌 (LUCA) 的分析给出了答案：携带CDKN2AHOM的肺癌患者，其转录组显示出更强的KRAS信号激活，且患者生存期更差。这意味着，虽然在肺部前期无需拆除刹车也能加速，但在后期，彻底摘除刹车能够换来更极致的增殖优势。

为了验证这种“表观遗传状态决定突变压力”的逻辑是否具有普适性，研究人员将目光投向了涵盖一万多个人类样本、三十三种癌症类型的TCGA泛癌数据集。他们挑选了膀胱癌、直肠癌、胃癌和子宫癌等常发生KRAS突变的癌症。

数据完美契合了模型预测。通过单细胞测序数据分析健康组织的基线状态，研究人员发现：膀胱上皮细胞具有极高的基线CDKN2A转录活性，与此对应，在KRAS突变的膀胱癌中，出现CDKN2AHOM的比例高达50%，并且频繁伴随KRASMUT-iGD。相反，在基线CDKN2A活性极低的直肠癌和子宫癌中，CDKN2AHOM的发生率均为0%。胃癌则处于中间状态，发生率为12%。

这就解释了为什么同样的突变会在不同器官中造就不同的疾病进程。小鼠肺部的致癌过程通常是迅速且多灶性的（因为刹车机制本来就不强）；而在小鼠胰腺中，由于强大的CDKN2A防线，通常需要漫长的时间才能突围，且往往只形成单发性肿瘤。

这份基于庞大疾病模型资源的研究，为我们描绘了一幅令人震撼的微观生态图景。它告诉我们，癌症基因组图谱上那些看似纷繁复杂、甚至充满随机性的突变频率、纯合或杂合状态，其实都是受到底层规则严格约束的必然结果。

这些底层规则，根植于细胞的谱系身份之中。胰腺等增殖静止器官，需要致癌基因以高剂量强制拉动“去分化”，这必然会唤醒活跃的抑癌基因，从而形成了一套必须按特定顺序破解的“密码”。而像肠道这样高度增殖、不断更新的器官，其首要任务不是去分化，而是通过WNT通路阻断分化并避免脱落；由于其特殊的生理功能，它自带对部分抑癌基因的免疫（通过多梳复合物压缩染色质），从而演化出了截然不同的突变轨迹。

这种将抑癌基因视为一个具有“连续谱”功能的连续统模型 (Continuum model)，以及在不同组织背景下基因剂量效应展现出的巨大差异，不仅解答了基础生物学中的深刻疑问，更为未来的癌症防治提供了极具启发的思路。当我们在临床上遇到特定的共突变模式或基因拷贝数改变时，我们或许能够通过解读其背后的演化法则，更精准地预测肿瘤的耐药机制，并找到干预的黄金窗口。在这场隐秘的微观演化史中，理解了生命的底层规则，我们才能真正把握住对抗疾病的先机。

参考文献

Mueller, S., de Andrade Krätzig, N., Tschurtschenthaler, M. et al. A disease model resource reveals core principles of tissue-specific cancer evolution. Nature (2026). https://doi.org/10.1038/s41586-026-10187-2

声明：本文仅用于分享，不代表平台立场，如涉及版权等问题，请尽快联系我们，我们第一时间更正，谢谢！