IF31.7｜化疗重塑造血干细胞克隆演化并加速克隆性衰老

单细胞全基因组测序揭示化疗直接诱变并重塑HSPC克隆结构，加速克隆性衰老，为t-MN起源提供新模型。

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

研究背景

造血干细胞（HSC）随年龄增长以每年约16-25个单核苷酸变异（SNV）的速率线性积累体细胞突变。部分突变赋予细胞适应性优势，驱动克隆性造血（CH），导致HSC克隆多样性随年龄下降，并增加血液恶性肿瘤风险。外源性压力如化疗可改变突变率和细胞适应性景观，但化疗对正常HSC基因组及群体动态的直接影响尚不清楚。治疗相关髓系肿瘤（t-MN）是化疗的罕见严重并发症，常携带DNA损伤反应（DDR）通路基因突变（如TP53、PPM1D），提示化疗可能选择性扩增携带这些突变的克隆。然而，化疗如何影响整体HSC群体结构、个体HSC基因组，以及这些效应如何促成t-MN发生，仍缺乏系统研究。

科学问题

卡在哪里： 既往研究多基于bulk测序，无法解析单个HSPC的基因组和克隆结构；动物模型虽提示化疗促进TP53/PPM1D突变克隆扩增，但缺乏人类体内直接证据；化疗诱导的突变是否直接驱动t-MN发生亦不明确。

本文要回答： 利用单细胞克隆全基因组测序，系统评估化疗对人类HSPC基因组和克隆动态的影响，并整合匹配的t-MN基因组，追溯t-MN的克隆起源和演化路径。

技术路线

作者主要用了：单细胞克隆全基因组测序、突变特征分析、系统发育树重建、Moran模型模拟、小鼠移植模型。

作者收集10例多发性骨髓瘤患者动员的外周血干细胞（PBSC）和6例正常供者PBSC，在甲基纤维素半固体培养基中培养单细胞来源的HSPC克隆，对1,276个克隆进行全基因组测序（WGS）。通过比较化疗暴露组与正常对照的突变负荷、突变特征和克隆结构，评估化疗的诱变和选择效应。进一步整合匹配的t-MN样本bulk WGS数据，进行系统发育分析，追溯t-MN的克隆起源。同时利用Moran模型模拟HSPC群体动态，验证化疗对克隆多样性的影响，并通过小鼠骨髓移植模型验证TP53/PPM1D突变对突变率的影响。

图 1　Fig. 1。图1展示了实验流程和队列概况。a部分为实验设计示意图；b部分总结了10例患者的临床特征、化疗暴露、克隆数和t-MN发生情况，红色框标记发生t-MN的患者；c部分为驱动突变和拷贝数变异的Oncoplot，显示PPM1D和TP53在化疗患者中富集，而正常供者中DNMT3A突变仅见于老年个体。该图奠定了研究队列和主要遗传变异基础。

核心亮点

亮点 1｜美法仑显著增加HSPC突变负荷

在10例患者中，8例的HSPC突变负荷与年龄匹配的正常对照相当（正常突变率17.2 SNV/年），但两例接受美法仑治疗的患者（PID0003和PID0004）突变负荷增加约2倍以上。突变特征分析鉴定出三个特征：SBS1、SBS5和SBS-C。SBS-C仅见于美法仑暴露患者，与已知美法仑相关特征SBS-MM1余弦相似度达90%。在PID0003和PID0004中，SBS-C分别贡献约60%和50%的突变。系统发育分析显示SBS-C突变在治疗后获得，证实美法仑直接诱导HSPC基因组突变。环磷酰胺虽同为烷化剂，但未观察到突变增加，可能因HSC通过醛脱氢酶将其代谢为无活性形式。

图 2　Fig. 2。图2聚焦突变负荷和突变特征。a部分散点图显示大多数化疗患者HSPC突变数与年龄匹配的正常突变率一致，但两例美法仑暴露患者（PID0003和PID0004）突变数显著升高；b部分比较了不同研究的正常突变率；c部分展示三个突变特征，其中SBS-C与美法仑相关特征SBS-MM1高度相似；d部分显示SBS-C在美法仑患者中贡献大部分突变；e部分系统发育树显示SBS-C突变在治疗后获得。该图证实美法仑直接诱导HSPC突变。

亮点 2｜化疗重塑HSPC克隆结构并降低多样性

基于时钟样突变（SBS1+SBS5）重建系统发育树，发现化疗后HSPC呈现显著寡克隆性，多个大克隆平行扩增。Simpson均匀度指数显示，部分化疗患者（年龄46-65岁）的克隆多样性与75岁以上正常个体相当。细胞毒性化疗（美法仑、环磷酰胺、多柔比星、长春新碱）处理组的均匀度显著低于非细胞毒性治疗组（来那度胺、沙利度胺、硼替佐米）。Moran模型模拟显示，化疗导致克隆多样性立即且持续下降，当引入化疗耐药突变（如TP53/PPM1D）时，多样性下降更显著。这些结果表明化疗加速HSPC克隆性衰老。

图 3　Fig. 3。图3展示基于时钟样突变重建的系统发育树。a部分为接受细胞毒性化疗的患者，b部分为非细胞毒性治疗患者。树中标注了驱动突变和拷贝数变异。可见多个平行克隆携带不同TP53或PPM1D突变，显示趋同进化。该图直观呈现化疗后HSPC的寡克隆结构和驱动突变克隆的扩张。

亮点 3｜DDR通路突变克隆的趋同进化

在化疗后HSPC中，57个扩增克隆中有30个（53%）携带驱动突变，而正常对照中仅17%。TP53和PPM1D突变频繁出现趋同进化：PID0004中检测到6个平行克隆，各携带不同PPM1D突变；PID0003、PID0005和PID0008中观察到多个不同TP53突变克隆。这些克隆在第二次ASCT时仍存在，表明移植后重新扩增。值得注意的是，98%的TP53突变克隆（106/108）拷贝数正常，仅2个克隆出现17p LOH，提示TP53突变在CH阶段并不伴随基因组不稳定性。

图 4　Fig. 4。图4展示克隆多样性分析。a部分比较化疗后HSPC与正常对照的Simpson均匀度指数，显示部分患者多样性显著降低；b部分与年龄匹配正常对照比较，确认化疗组多样性下降；c部分显示细胞毒性化疗组多样性低于非细胞毒性组；d和e部分为Moran模型模拟结果，显示化疗导致多样性立即下降，且耐药突变加剧下降。该图支持化疗加速克隆性衰老的结论。

亮点 4｜TP53/PPM1D突变本身不增加突变率

比较野生型克隆与携带驱动突变的克隆，发现部分TP53/PPM1D突变克隆突变率显著升高，但并非所有驱动突变克隆均如此。一些无驱动突变的克隆也显示高突变率，且与端粒缩短相关。突变率与端粒长度呈负相关，提示高突变率可能源于加速的细胞分裂而非驱动突变直接效应。小鼠模型验证：Trp53−/−和Ppm1dmut/+ HSPC的突变率与野生型相当，顺铂处理增加突变率但突变型和野生型无差异，且拷贝数变异无差异。表明在CH阶段，TP53/PPM1D突变不直接导致基因组不稳定性。

图 5　Fig. 5。图5展示突变率与端粒长度的克隆特异性分析。a和b部分分别显示携带驱动突变的克隆在突变数和端粒长度上的分布；c部分比较不同克隆的突变率，显示部分TP53/PPM1D突变克隆突变率升高，但并非普遍；d部分显示高突变克隆端粒缩短；e部分显示突变率与端粒长度负相关；f部分展示驱动突变的突变特征归属，部分驱动突变可能由美法仑诱导。该图提示高突变率可能源于加速增殖而非驱动突变直接效应。

亮点 5｜t-MN克隆起源的两种路径

整合HSPC克隆和t-MN bulk WGS进行系统发育分析，在5/9例患者中确定了t-MN的最近共同祖先（MRCA）。例如PID0002中，t-MDS克隆起源于携带TET2、TP53突变和17p LOH的克隆；PID0005中，t-MDS起源于携带TP53和DNMT3A突变的克隆。这些病例的共同特征是继发性突变和TP53双等位基因失活是恶性转化的关键。在4例未找到MRCA的患者中，2例t-MN样本显示SBS-C特征，提示可能起源于未动员的骨髓HSPC，这些细胞在预处理时暴露于高剂量美法仑。因此提出t-MN的两种发育路径：多数起源于移植的PBSC，少数起源于未动员的骨髓HSPC。

生信可带走

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

• 方法栈： 单细胞克隆全基因组测序、突变特征分析、系统发育树重建、Moran模型模拟、小鼠移植模型
• 公开数据： 原文未给出公开组学登录号
• 代码： 原文未给出公开 GitHub/GitLab 仓库
• 谁该点开原文： 从事克隆性造血、治疗相关髓系肿瘤、单细胞基因组学和化疗基因组毒理学研究的科研人员。

带走一句

单细胞克隆WGS能区分化疗的直接诱变和选择效应，为解析t-MN起源提供高分辨率框架；分析时需注意突变特征与克隆结构解耦，避免将驱动突变与高突变率简单关联。