IF50.5｜驱动突变衰减重塑肠道转化景观

驱动突变出现的先后顺序决定克隆是被正选择还是负选择，并塑造后续肿瘤发展。

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

研究背景

结直肠癌（CRC）传统上被认为通过APC肿瘤抑制基因及其他驱动基因的逐步突变，伴随正选择克隆的扩增而发展。然而，近年研究发现许多癌前病变是多克隆的，包含表达不同突变APC蛋白的多个克隆。克隆合作可能涉及表达N端APC截短蛋白的克隆被表达更C端截短蛋白的克隆招募，提示不同Apc突变具有不等的转化景观。这一现象可能与“恰到好处”模型相关，即通过保留部分β-catenin结合能力的APC突变组合实现最佳WNT通路扰动。此前工作表明多克隆性依赖于背景：在含有KRAS激活突变的上皮背景上起始时，大多数肿瘤为单克隆。本研究旨在理解由KRAS及其他驱动突变在正常肠上皮中创造的背景如何影响Apc及其他CRC驱动基因突变的转化与选择景观。

科学问题

卡在哪里： 此前研究已发现不同Apc突变具有不等的转化能力，且多克隆性依赖于背景，但尚不清楚正常肠上皮中不同驱动突变背景如何定量影响后续Apc和Ctnnb1突变的固定与丢失，以及这些突变在无背景支持下是否经历负选择。

本文要回答： 本研究旨在确定正常小鼠肠上皮中不同驱动突变背景如何改变Apc和Ctnnb1突变的转化与克隆选择景观，并检验人类CRC中是否存在类似的启动事件证据。

技术路线

作者主要用了：条件性基因敲入/敲除小鼠模型、ENU化学诱变、靶向捕获测序、dNdScv选择分析、突变特征推断、克隆衰减模拟。

作者首先通过他莫昔芬诱导在肠上皮中建立多种促癌背景（如KrasG12D、Trp53null、Fbxw7null等），然后用ENU化学诱变剂诱导随机突变，观察肿瘤发生及驱动基因突变谱。接着，通过反转处理顺序（先ENU后他莫昔芬）进行“救援”实验，以评估在无背景支持下ENU诱导的驱动突变是否被负选择清除。最后，利用人类CRC基因组数据（GENIE和INCISE队列）分析APC突变分布与KRAS突变状态的关系，寻找启动事件的证据。

图 1　Fig. 1。图1展示了实验设计（先他莫昔芬诱导启动背景，再ENU诱变）以及不同启动背景对肿瘤负担和生存的影响。对照小鼠中位生存期251天，平均肿瘤数3.1个；而KrasG12D启动小鼠中位生存期仅37天，平均肿瘤数921个。dNdScv分析显示Apc无义突变和Ctnnb1错义突变被正选择。不同背景下Apc和Ctnnb1驱动突变的比例差异显著，且KrasG12D启动导致Ctnnb1驱动肿瘤增加596倍，Apc驱动肿瘤增加58倍。

核心亮点

亮点 1｜启动背景定量改变肿瘤易感性

不同启动背景对ENU诱变的肿瘤易感性差异巨大。KrasG12D启动小鼠中位生存期仅37天，平均肿瘤数921个；而对照小鼠中位生存期251天，平均肿瘤数3.1个。Fbxw7null和Trp53null启动小鼠中位生存期分别为64天和120天，肿瘤数分别为157和40个。肿瘤主要由Apc无义突变或Ctnnb1错义突变驱动，但比例因背景而异：KrasG12D启动小鼠中约85%的肿瘤含Ctnnb1突变，而Apchet启动小鼠中约90%的肿瘤含Apc突变。这些结果表明不同启动背景定量改变转化率，并影响驱动基因突变的选择。

图 2　Fig. 2。图2显示Ctnnb1突变几乎全部为exon 3的错义突变，集中在七个密码子。p.G34E和p.T41I占41%，但不同背景下优势突变不同。通过突变特征推断，发现p.T41I和p.G34E的高频主要由高突变易感性驱动，而p.I35S的O/E比值为4，表明受到正选择。不同背景对特定突变的O/E比值影响显著，如对照和Trp53hetN1-ICDhet背景中p.I35S的O/E比值高达14和11，而KrasG12D背景中O/E比值范围最小。

亮点 2｜Ctnnb1突变谱受背景塑造

在348个Ctnnb1驱动突变中，几乎所有都是错义突变，聚集在exon 3的七个密码子。p.G34E和p.T41I占所有突变的41%，但不同背景下优势突变不同：KrasG12D和Trp53null背景中p.G34E占优，而Fbxw7null背景中p.S37F和p.D32G更常见。通过突变特征推断，发现p.T41I和p.G34E的高频主要由高突变易感性驱动，但p.I35S的O/E比值为4，表明受到正选择。不同背景对特定突变的O/E比值影响显著，如对照和Trp53hetN1-ICDhet背景中p.I35S的O/E比值高达14和11，而KrasG12D背景中O/E比值范围最小，提示其创造了广泛容许的转化景观。

图 3　Fig. 3。图3展示了Apc无义突变在N端区域的分布及不同背景下的域偏好。总体而言，bin A和bin E的突变被负选择，而bin B和D被正选择。Ptennull和Apchet背景中bin A突变显著减少，而Apchet背景中bin B突变富集。双等位基因突变组合分析显示，对照和Trp53相关背景倾向于BB组合，而Ptennull背景中突变向远端偏移，保留完整Armadillo域但丢失bin E。

亮点 3｜Apc突变域偏好因背景而异

Apc无义突变广泛分布于N端区域，但不同背景对突变域的选择压力不同。总体而言，bin A（最N端）和bin E（含20氨基酸重复）的突变被负选择，而bin B和D被正选择。Ptennull和Apchet背景中bin A突变显著减少，而Apchet背景中bin B（Armadillo域）突变富集。双等位基因突变组合分析显示，对照和Trp53相关背景倾向于BB组合，而Ptennull背景中突变向远端偏移，保留完整Armadillo域但丢失bin E。这表明不同背景对Apc突变的选择压力不同，Ptennull背景独特地需要完整的Armadillo域。

图 4　Fig. 4。图4展示了救援实验的结果。KrasG12D救援30天后，小肠肿瘤数从388个降至11个，减少97%。不同Ctnnb1突变丢失效率不同：p.S33、p.G34、p.D32、p.S37、p.T41和p.I35的丢失概率依次降低。Apc突变也发生衰减，但分层分析显示bin A和B的突变在KrasG12D、Trp53null和Fbxw7null救援小鼠中呈现显著负选择，而bin C-E突变衰减接近中性。在Apchet救援中，bin B-D突变强负选择，而bin A和E突变中性衰减。

亮点 4｜救援实验揭示驱动突变衰减

通过反转处理顺序（先ENU后他莫昔芬）的救援实验，作者发现ENU诱导的驱动突变在无背景支持下经历显著负选择。KrasG12D救援30天后，小肠肿瘤数从388个降至11个，减少97%。不同Ctnnb1突变丢失效率不同：p.S33、p.G34、p.D32、p.S37、p.T41和p.I35的丢失概率依次降低。Apc突变也发生衰减，但分层分析显示bin A和B的突变在KrasG12D、Trp53null和Fbxw7null救援小鼠中呈现显著负选择，而bin C-E突变衰减接近中性。在Apchet救援中，bin B-D突变强负选择，而bin A和E突变中性衰减。这些结果表明不同Apc突变的转化潜力与负选择强度相关。

图 5　Fig. 5。图5展示了人类CRC中APC突变分布与KRAS突变状态的关系。在GENIE队列中，bin E的突变显著富集（log2[FC]=2.1），而bin B显著减少（log2[FC]=-1.7）。KRAS突变型CRC倾向于保留两个20氨基酸重复，而KRAS野生型保留一个。在INCISE息肉队列中，KRAS/APC VAF比值较高的息肉更可能保留更多20氨基酸重复，提示启动事件可能影响息肉进展风险。

亮点 5｜突变顺序决定选择结果

在Fbxw7null救援实验中观察到上位效应：先有Fbxw7缺失再获得Ctnnb1 p.I35S突变时，肿瘤形成受抑制；但先有Ctnnb1 p.I35S突变再丢失Fbxw7时，肿瘤形成被促进。类似地，先有Apc bin A截短突变再丢失Fbxw7也促进肿瘤形成。这些结果表明突变发生的顺序决定克隆是被正选择还是负选择，与之前报道的Fbxw7和Apc之间的上位关系一致。

亮点 6｜人类CRC中存在启动证据

在GENIE队列的9,278例CRC中，APC突变分布与预期突变特征比较显示，bin E（含20氨基酸重复）的突变显著富集（log2[FC]=2.1），而bin B（Armadillo域）显著减少（log2[FC]=-1.7），与小鼠救援实验中观察到的负选择一致。进一步分析发现，KRAS突变型CRC倾向于保留两个20氨基酸重复，而KRAS野生型保留一个。这一趋势在INCISE息肉队列中也存在，且KRAS/APC VAF比值较高的息肉（提示KRAS突变较早）更可能保留更多20氨基酸重复，提示启动事件可能影响息肉进展风险。

生信可带走

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

• 方法栈： 条件性基因敲入/敲除小鼠模型、ENU化学诱变、靶向捕获测序、dNdScv选择分析、突变特征推断、克隆衰减模拟
• 公开数据： 原文未给出公开组学登录号
• 代码： 原文未给出公开 GitHub/GitLab 仓库
• 谁该点开原文： 研究结直肠癌演化、克隆选择、驱动突变功能的生信和实验研究者。

带走一句

驱动突变的命运取决于它们出现的顺序和背景。分析肿瘤基因组时，不仅要看突变是否存在，还要考虑突变发生的先后和共存背景，才能理解克隆选择的方向。