IF50.5｜异基因造血干细胞移植后克隆动态

移植后干细胞克隆多样性加速丢失，相当于额外衰老10-15年。

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

研究背景

异基因造血干细胞移植（HCT）是治疗血液恶性肿瘤和骨髓衰竭性疾病的重要手段，其成功依赖于供体造血干/祖细胞（HSPC）在受者体内长期植入并重建多系造血。既往研究通过逆转录病毒整合位点、慢病毒条形码等克隆示踪技术，在模式生物或自体基因治疗中估计了植入干细胞数量，但无法解析移植前已存在的克隆结构对植入的影响。体细胞突变作为内源性谱系标记，可回溯至胚胎发育期，已用于描绘健康人一生中造血系统的克隆动态，但在异基因移植背景下的长期克隆演化尚缺乏系统研究。

科学问题

卡在哪里： 此前对人类异基因HCT后长期植入的干细胞数量、克隆多样性变化及选择压力来源缺乏直接测量。动物模型和自体基因治疗中的克隆示踪方法无法反映异体移植中免疫排斥、化疗损伤骨髓微环境等独特挑战，且无法评估供体原有克隆结构对植入的影响。

本文要回答： 利用全基因组测序重建供受者对造血系统系统发育树，定量长期植入干细胞数量，比较供受者克隆多样性差异，并解析移植特异性选择压力的模式。

技术路线

作者主要用了：单细胞来源造血集落全基因组测序、系统发育重建、深度靶向测序、近似贝叶斯计算。

研究纳入10对HLA全相合同胞供受者，移植后9-31年采集外周血，分离CD34+ HSPC进行单细胞来源集落培养，对每个集落进行全基因组测序（WGS），利用体细胞突变重建供者和受者各自的造血系统系统发育树。同时，对分选的粒细胞、单核细胞、B细胞和T细胞进行定制靶向深度测序，以量化各克隆在不同成熟血细胞谱系中的贡献。通过近似贝叶斯计算（ABC）模拟造血衰老和移植瓶颈，估计长期植入干细胞数量，并比较供受者克隆多样性及选择模式。

图 1　Fig. 1。图1展示了研究设计和系统发育树构建流程。a部分显示从10对供受者采集外周血，分离CD34+细胞进行单细胞集落培养，对每个集落进行全基因组测序，利用体细胞突变重建系统发育树，并对成熟血细胞亚群进行靶向测序。b部分展示了一对供受者各自的系统发育树，显示供者来源的造血系统在受者体内重建后的克隆结构。

核心亮点

亮点 1｜长期植入干细胞数量估计

通过ABC模拟，估计每个受者长期植入的造血干细胞数量中位数在700至25,000之间。年轻供者（移植时18-47岁）的植入数量显著高于年长供者（50-66岁），后者约低10倍。植入数量与输注的CD34+细胞剂量无显著相关性，但与供者年龄呈负相关。这一结果提示，供者年龄是决定长期植入干细胞池大小的关键因素，而非输注细胞数量。

图 2　Fig. 2。图2展示了长期植入造血干细胞数量的估计结果。a部分显示每对供受者通过ABC估计的植入干细胞数量后验分布，中位数从700到25,000不等。b部分显示植入干细胞数量与输注CD34+细胞剂量无显著相关性。c部分显示植入干细胞数量与供者年龄呈负相关，年轻供者植入数量更高。

亮点 2｜受者克隆多样性加速丢失

与匹配供者相比，受者造血系统克隆多样性显著降低，相当于额外衰老10-15年。受者中扩增克隆（≥2%造血贡献）的比例平均增加23个百分点。系统发育多样性指标（平均成对距离、平均最近分类单元距离）均显示受者多样性下降。这种加速的寡克隆化不能用简单的移植瓶颈解释，提示存在移植特异性的选择压力。

图 3　Fig. 3。图3展示了受者克隆多样性丢失和谱系偏倚。a部分显示各供受者对中扩增克隆（≥2%造血贡献）的数量和大小，受者中扩增克隆比例更高。b部分显示受者系统发育多样性低于供者。c部分显示受者系统发育年龄比供者平均大12年。d-f部分显示各成熟血细胞亚群中克隆多样性和谱系偏倚，T细胞中扩增克隆贡献最低。

亮点 3｜两种移植特异性选择模式

系统发育树分析揭示两种不同的选择模式：修剪选择（pruning selection）和生长选择（growth selection）。修剪选择中，克隆扩增的细胞分裂发生在供者体内，移植前已存在，其优势来自动员、采集、体外存活或初始归巢等环节。生长选择中，克隆扩增的细胞分裂发生在受者骨髓内，移植后发生，尤其在携带多个驱动突变的克隆中更为显著。这两种模式解释了受者克隆多样性丢失的机制。

图 4　Fig. 4。图4展示了两种选择模式的系统发育树证据。a部分显示pair 3中受者存在一个在供者中不明显的巨大克隆扩增，其合并事件发生在移植前，支持修剪选择。b部分显示pair 9中受者存在一个携带Y染色体丢失的克隆扩增，其合并事件发生在移植时，支持生长选择。c部分显示pair 7中类似生长选择的模式。

亮点 4｜驱动突变克隆的移植后动态

比较受者与供者单核细胞中52个已知驱动突变的变异等位基因频率（VAF），发现大多数克隆（33/52）在受者与供者间无显著差异，14个克隆在受者中扩增，5个克隆缩小。CHEK2突变克隆一致在受者中扩增。携带两个或以上驱动突变的克隆在受者中显著更大，而单驱动突变克隆则无一致趋势。dN/dS分析显示全基因组编码区存在净正选择，约1/11的非同义突变受正选择，提示存在大量未知的移植特异性驱动基因。

图 5　Fig. 5。图5展示了驱动突变克隆的移植后动态。a部分显示各驱动基因突变在受者与供者单核细胞中VAF的log2倍数变化，CHEK2突变一致在受者中扩增。b部分显示携带多个驱动突变的克隆在受者中显著更大，而单驱动突变克隆无一致趋势。

亮点 5｜APOBEC突变在受者中富集

在2,824个造血集落中，63个（2.2%）携带APOBEC特征突变，其中61个来自受者，且几乎全部发生在移植后。APOBEC突变与驱动突变、预处理方案、供者性别、干细胞来源或克隆扩增状态无关。其基因组分布特征（靠近十字形反向重复、高Alu密度、低GC含量）与APOBEC酶活性一致。受者免疫抑制状态下病毒感染可能诱导APOBEC活性，导致移植后突变负荷增加。

生信可带走

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

• 方法栈： 单细胞来源造血集落全基因组测序、系统发育重建、深度靶向测序、近似贝叶斯计算
• 公开数据： 原文未给出公开组学登录号
• 代码： 公开仓库 http://www.github.com/nickwilliamssanger/Rsimpop、https://github.com/MathijsSanders/SangerLCMFiltering、https://github.com/cancerit/vafCorrect；本地已克隆：Rsimpop、hdp、SangerLCMFiltering、vafCorrect、Clonal_dynamics_of_HSCT
• 谁该点开原文： 从事造血干细胞移植、克隆造血和体细胞进化研究的科研人员及临床医生。

带走一句

移植后克隆多样性丢失由移植特异性选择驱动，而非简单瓶颈；分析克隆动态时需区分修剪选择和生长选择，并关注未知驱动基因。