IF50.5｜正常膀胱克隆景观的性别与吸烟偏倚

超深双链测序揭示男性膀胱对截短驱动突变的正选择更强，吸烟促进TERT启动子克隆扩增。

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

研究背景

男性膀胱癌风险约为女性的四倍，吸烟也是主要风险因素，但机制不明。肿瘤发生由体细胞突变驱动，正常组织中的克隆扩增是癌变的第一步。既往研究已利用超深测序在正常食管、子宫内膜等组织中发现了驱动突变克隆，并揭示了年龄相关的克隆扩增。在膀胱中，已有研究通过显微切割和全基因组测序发现正常尿路上皮存在携带癌症基因突变的克隆，但受限于样本量和测序深度，无法系统比较不同性别和吸烟史个体的克隆选择差异。

科学问题

卡在哪里： 此前对正常膀胱克隆的研究样本量小、测序深度有限，无法准确估计每个基因的正选择强度，因此无法回答性别和吸烟是否以及如何影响正常尿路上皮的克隆选择。

本文要回答： 利用超深双链测序系统描绘正常膀胱尿路上皮的克隆突变图谱，并定量评估性别和吸烟对驱动突变正选择的影响。

技术路线

作者主要用了：超深双链DNA测序、dN/dS正选择分析、突变特征提取、自然饱和诱变。

作者从45例尸检供体的膀胱穹顶和三角区刷取上皮细胞，获得79份多克隆样本，对16个膀胱癌相关基因（含TERT启动子）进行超深双链测序（中位深度约5,164×）。利用新开发的Omega方法计算每个基因的dN/dS值以量化正选择强度，并在样本水平估计驱动突变数量。随后通过混合效应回归模型分析年龄、性别、吸烟史等临床变量与突变密度及选择强度的关联，同时利用肿瘤基因组数据验证正常组织中的发现。

图 1　Fig. 1。图1展示研究设计和突变全景。a为膀胱刷取和双链测序示意图；b显示正常膀胱中检测到的SNV数量远超892例肿瘤；c为79份样本的突变数量和密度，按供体年龄、性别、吸烟史排列；d和e展示提取的APOBEC、年龄和化疗相关突变特征。

核心亮点

亮点 1｜正常膀胱突变负荷远超肿瘤

在79份正常尿路上皮样本中，共鉴定出64,278个体细胞突变，中位774个/样本，是892例膀胱肿瘤中相同基因突变总数的约16倍。突变密度与年龄及化疗暴露相关，提取到APOBEC、年龄相关和化疗相关三种主要突变特征，但未发现吸烟相关特征。

图 2　Fig. 2。图2展示RBM10和TP53的突变分布及选择强度。a显示截短、错义和同义突变沿编码区的分布，可见截短突变在RBM10中富集；b为dN/dS值，显示RBM10和TP53截短突变受强正选择；c显示FGFR3截短突变受负选择；d展示TERT启动子突变分布，激活突变集中于两个热点。

亮点 2｜13个基因在正常膀胱中受正选择

通过dN/dS、功能影响偏倚、3D聚类和移码插入缺失富集四种方法，确认RBM10、KDM6A、STAG2、KMT2D、ARID1A、CDKN1A、TP53、EP300、NOTCH2、CREBBP、FOXQ1、KMT2C、RB1共13个基因在正常尿路上皮中驱动克隆扩增。FGFR3截短突变则受负选择。TERT启动子中56个激活突变中至少55个为驱动突变。

图 3　Fig. 3。图3展示性别偏倚。a为13个基因在男女供体中的驱动SNV和蛋白影响插入缺失数量热图；b显示年龄与蛋白影响和非蛋白影响突变密度的关系，前者斜率更大；c显示年龄与dN/dS值的关联；d和e显示男性RBM10、CDKN1A和ARID1A截短突变dN/dS值显著高于女性。

亮点 3｜男性对截短突变的正选择更强

多变量回归显示，男性尿路上皮中RBM10、CDKN1A和ARID1A的截短突变dN/dS值显著高于女性，且独立于年龄、吸烟、饮酒、BMI和化疗暴露。STAG2也呈类似趋势。在膀胱肿瘤中，男性RBM10和CDKN1A截短突变频率也显著更高，提示正常组织中的选择偏倚可能反映在肿瘤突变谱中。

图 4　Fig. 4。图4展示TERT启动子突变与年龄和吸烟的关系。a显示55岁以下、55岁以上不吸烟和55岁以上吸烟三组供体中激活突变的数量；b为三组突变密度，吸烟组显著升高；c显示吸烟与年龄交互作用对TERT启动子突变密度的显著影响。

亮点 4｜TERT启动子突变与年龄和吸烟相关

TERT启动子激活突变几乎全部出现在55岁以上供体，且在该年龄组中，有吸烟史者的突变率显著高于不吸烟者。该关联独立于性别、BMI和化疗暴露。吸烟并未增加总体突变密度，提示烟草致癌物可能通过促进携带TERT启动子突变的克隆扩增而非直接诱变发挥作用。

图 5　Fig. 5。图5展示自然饱和诱变。a为各基因氨基酸残基突变次数分布；b为TP53、EP300和FGFR3的理论与观察饱和动力学曲线，正选择基因观察曲线高于理论曲线，FGFR3则低于；c为TP53位点选择图谱，显示DNA结合域显著富集。

亮点 5｜自然饱和诱变揭示位点选择

通过超深测序，在TP53蛋白的394个氨基酸残基中观察到230个（58%）发生突变。与肿瘤数据相比，正常组织覆盖的突变位点比例更高。利用下采样构建的观察动力学曲线显示，正选择基因的突变积累快于中性预期，而FGFR3则慢于预期。位点选择分析发现TP53的DNA结合域、EP300和CREBBP的HAT结构域等存在显著选择信号，且与实验饱和诱变结果一致。

生信可带走

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

• 方法栈： 超深双链DNA测序、dN/dS正选择分析、突变特征提取、自然饱和诱变
• 公开数据： 原文未给出公开组学登录号
• 代码： 原文未给出公开 GitHub/GitLab 仓库
• 谁该点开原文： 关注癌症基因组学、体细胞进化和正常组织克隆动态的研究者，以及从事超深测序数据分析的生信人员。

带走一句

超深双链测序结合dN/dS分析可在正常组织中定量评估驱动基因的选择强度，并揭示癌症风险因素对克隆演化的影响，为研究体细胞进化和癌症预防提供了新范式。