IF50.5｜结直肠癌突变过程的地理与年龄差异

细菌源诱变剂colibactin在早发性结直肠癌中富集，提示生命早期暴露可能驱动年轻化趋势。

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

研究背景

结直肠癌发病率存在显著地理差异，且近二十年来早发性结直肠癌（<50岁）在许多国家翻倍增长，原因不明。传统流行病学难以捕捉外源性诱变暴露，而突变特征分析可揭示基因组中遗留的诱变过程指纹。既往结直肠癌全基因组研究多集中于北美和欧洲，早发性病例数量有限，且未系统比较不同地理区域。已知colibactin（由携带pks岛的大肠杆菌产生）可诱导特征性突变特征SBS88和ID18，并在正常结直肠上皮中早期印迹，但其与早发性结直肠癌的关联尚缺乏大规模证据。

科学问题

卡在哪里： 此前缺乏跨多国、大样本的结直肠癌全基因组数据，无法评估突变过程的地理差异及其与早发性结直肠癌的关联；colibactin暴露与年龄的关系及其在肿瘤发生中的作用时间点亦未明确。

本文要回答： 利用11国981例结直肠癌全基因组，系统比较突变特征的地理和年龄差异，并检验colibactin相关突变特征是否在早发性病例中富集及其在肿瘤演化中的时序。

技术路线

作者主要用了：全基因组测序、突变特征提取、进化分析、驱动基因检测、微生物组分析。

研究纳入11国981例原发结直肠癌及配对正常组织，进行全基因组测序（中位覆盖度53×/27×）。首先按分子亚型分为MSS和MSI，聚焦802例MSS进行突变特征提取。利用SigProfilerExtractor和mSigHdp两种方法提取SBS、ID、DBS、CN和SV特征，并与COSMIC参考特征比对。随后通过逻辑回归评估特征在不同国家和年龄组中的富集，利用进化分析（MutationTimeR）将突变分为早期克隆、晚期克隆和亚克隆，以确定colibactin特征的时序。最后结合驱动基因分析和微生物组测序数据，评估colibactin暴露与APC驱动突变的关系。

图 1　Fig. 1。图1展示981例结直肠癌的地理分布、亚部位和分子亚型。a显示11个国家病例数及早发性比例，颜色代表ASR。b显示肿瘤亚部位分布，早发性病例在远端结肠和直肠富集1.88倍。c显示MSS和MSI的突变负担差异：MSI具有高SBS和ID负担但低基因组畸变，MSS相反。d显示早发性MSS肿瘤SBS和ID负担略低。e-g显示早发与晚发肿瘤的平均突变谱高度相似（余弦相似度>0.97）。

核心亮点

亮点 1｜多国队列揭示特征地理差异

在802例MSS结直肠癌中，阿根廷富集SBS89（OR=28.0）、DBS8（OR=8.9）和新型ID_J（OR=9.6），三者显著共现；哥伦比亚富集SBS94（OR=19.7）、SBS_F（OR=10.7）和DBS6（OR=12.5）；俄罗斯富集SBS2和SBS_H，巴西富集CN_F，泰国DBS2减少。这些特征大多病因未知，提示不同国家存在特异性诱变暴露。此外，colibactin相关特征SBS88和ID18与结直肠癌ASR正相关，尤其与直肠癌ASR相关，表明微生物组来源的诱变剂可能影响发病率。

图 2　Fig. 2。图2展示突变特征的国家特异性富集。a为各特征在不同国家与所有其他国家的OR值，点大小代表显著性。b显示阿根廷富集SBS89、DBS8和ID_J。c显示哥伦比亚富集SBS94、SBS_F和DBS6。d-e显示SBS88、ID18、SBS1和SBS_H与ASR正相关，CN_F负相关。f显示SBS88和ID18与直肠癌ASR相关。

亮点 2｜colibactin特征在早发性癌症中富集

SBS88和ID18在早发性（<50岁）结直肠癌中分别富集2.5倍和4倍（q=0.006和q=3.7×10^-7），且随诊断年龄增加患病率显著下降（P trend=1.3×10^-4和2.0×10^-7）。以任一特征定义colibactin暴露，其在<40岁患者中比>70岁患者常见3.3倍。该效应在远端结肠和直肠更明显。此外，未知病因特征SBS_M和ID14也在早发性病例中富集。

图 3　Fig. 3。图3展示特征与诊断年龄的关联。a为早发与晚发组各特征的OR值，SBS88和ID18显著富集于早发组。b显示时钟样特征SBS1、SBS5、ID1、ID2在晚发组富集。c显示SBS88和ID18患病率随年龄增加而下降。d显示colibactin暴露者诊断年龄中位数更低（62 vs 67岁）。e显示该效应在远端结肠和直肠更显著。

亮点 3｜colibactin诱变发生在肿瘤早期

通过进化分析将突变分为早期克隆、晚期克隆和亚克隆，发现SBS88和ID18在早期克隆突变中显著富集（q=4.2×10^-4和q=6.1×10^-5），与时钟样特征SBS1、SBS5和ID1相似，表明colibactin诱变发生在肿瘤发生的早期阶段，可能早于大多数驱动事件。该模式在早发性和晚发性病例中均存在，提示colibactin暴露可能为早期事件，为后续肿瘤发生奠定基础。

图 4　Fig. 4。图4展示colibactin特征的时序和微生物组关联。a显示SBS88和ID18在早期克隆突变中富集，与时钟样特征相似。b显示colibactin暴露状态与肿瘤中pks+细菌检测无关联。c-d显示携带特征但未检测到pks+细菌的患者诊断年龄更早。

亮点 4｜colibactin驱动APC基因失活

在colibactin阳性肿瘤中，ID18贡献了约25%的APC驱动插入缺失（indels），而SBS88仅贡献3.9%的APC驱动替换。由于APC失活通常是结直肠癌发生的起始事件，colibactin诱导的APC移码突变可能在生命早期赋予细胞生长优势，加速肿瘤发生。此外，SBS88和ID18还贡献了染色质修饰基因的突变，进一步支持其致癌作用。

图 5　Fig. 5。图5展示驱动基因突变与colibactin的关联。a显示48个驱动基因的突变频率，APC、TP53、KRAS等最常见。b显示晚发性肿瘤总驱动突变数高于早发性。c-d显示APC驱动突变在晚发性病例中更常见。e-f显示在colibactin阳性肿瘤中，ID18贡献25.3%的APC驱动indels，而SBS88仅贡献3.9%的APC驱动替换。

亮点 5｜微生物组检测与特征不一致

在肿瘤样本中检测pks+细菌时，未发现SBS88或ID18阳性肿瘤与pks+细菌存在关联。然而，携带特征但未检测到pks+细菌的患者诊断年龄更早（P=1.3×10^-7），提示colibactin暴露可能发生在生命早期，随后微生物组发生改变。这一发现强调了突变特征作为历史暴露记录的价值，而非当前微生物状态的反映。

生信可带走

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

• 方法栈： 全基因组测序、突变特征提取、进化分析、驱动基因检测、微生物组分析
• 公开数据： 原文未给出公开组学登录号
• 代码： 原文未给出公开 GitHub/GitLab 仓库
• 谁该点开原文： 关注癌症基因组学、突变特征分析、微生物组与肿瘤发生交叉领域的研究者。

带走一句

突变特征可作为历史诱变暴露的分子档案，揭示地理和年龄特异的致癌过程；分析时应结合进化时序和驱动基因，以阐明诱变剂的功能影响。