IF31.7｜墨西哥城与英国生物库比较揭示克隆性造血的祖先特异性效应

最大规模非欧洲人群CH研究，发现CH频率与欧洲血统比例正相关，并鉴定出TCL1B位点祖先特异性变异。

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

研究背景

克隆性造血（CH）是造血干细胞随年龄积累体细胞突变并发生克隆扩增的现象，与血液肿瘤及心血管、肾脏等非肿瘤疾病风险升高相关。已知年龄、性别、吸烟、遗传变异等是CH的风险因素，但绝大多数CH研究集中在欧洲或欧洲血统为主的人群，遗传祖先对CH的影响尚不清楚。墨西哥城前瞻性研究（MCPS）是一个大型拉丁美洲人群队列，具有印第安、欧洲和非洲混合血统，为研究祖先背景与CH的关系提供了独特机会。

科学问题

卡在哪里： 此前CH研究主要基于欧洲血统人群，非欧洲人群数据缺乏，无法评估遗传祖先对CH频率和遗传易感性的影响。不同人群间CH频率差异及其遗传基础尚不明确，且非欧洲人群中CH相关的种系变异发现有限。

本文要回答： 比较MCPS与UKB的CH频率差异，分析MCPS内部欧洲血统比例与CH的关系，并通过全基因组和全外显子关联分析鉴定祖先特异性的CH遗传变异。

技术路线

作者主要用了：WES、WGS、GWAS、ExWAS、基因负荷分析、跨祖先荟萃分析。

作者首先对MCPS和UKB的WES数据采用统一流程检测CH，比较两队列CH频率。然后利用MCPS的混合血统，分析欧洲血统比例与CH频率的相关性。接着在MCPS中开展LTL GWAS和CH GWAS/ExWAS，鉴定祖先特异性变异。最后进行MCPS与UKB的跨祖先荟萃分析，发现新的CH易感位点，并通过基因负荷分析验证CHEK2关联。

图 1　Fig. 1。图1展示MCPS中CH驱动基因的突变人数排序（DNMT3A最多），CH频率随年龄增长，以及UKB与MCPS总体CH和基因特异性CH的频率比较，显示UKB在多数基因上频率更高。

核心亮点

亮点 1｜CH频率在MCPS中显著低于UKB

在136,401名MCPS和416,118名UKB参与者中，MCPS总体CH频率为3.12%，显著低于UKB的4.92%（P<2.2×10^-16）。年龄和性别匹配后差异依然存在（2.87% vs 4.55%）。逻辑回归调整年龄、性别和吸烟后，UKB相对MCPS的CH风险OR为1.69（95% CI 1.63-1.75，P=7.31×10^-185）。基因特异性分析显示，DNMT3A、TET2、ASXL1、PPM1D、TP53、JAK2和SRSF2的CH在UKB中均更常见，其中DNMT3A-CH的OR高达1.99。

图 2　Fig. 2。图2显示MCPS中欧洲血统比例>50%的个体CH频率高于印第安血统>50%的个体，且CH频率随欧洲基因组比例增加而升高；吸烟频率也随欧洲血统比例增加而升高。

亮点 2｜欧洲血统比例与CH频率正相关

在MCPS内部，欧洲血统比例>50%的个体CH频率为4.65%，而印第安血统>50%的个体为2.83%（P<2.2×10^-16）。以欧洲基因组比例每增加10%为连续变量，调整年龄、性别和吸烟后，CH风险随欧洲血统比例增加而升高（beta=0.84，95% CI 0.66-1.03，P=7.35×10^-19）。基因特异性分析显示，DNMT3A-CH、ASXL1-CH和SRSF2-CH与欧洲血统比例显著正相关。吸烟频率在欧洲血统比例高的个体中更高，但调整吸烟后关联仍存在。

图 3　Fig. 3。图3展示MCPS中LTL GWAS的曼哈顿图，标注了6个全基因组显著位点；LTL PRS随欧洲血统比例增加而升高；遗传预测LTL与总体CH、DNMT3A-CH、TET2-CH和SRSF2-CH正相关，与PPM1D-CH负相关。

亮点 3｜端粒长度部分解释祖先差异

利用9,598名有WGS数据的MCPS参与者推断白细胞端粒长度（LTL），发现LTL随年龄缩短。LTL GWAS鉴定出6个全基因组显著位点（TERC、NAF1、TERT、TERF1、STN1、TINF2），均为已知位点。基于MCPS构建的LTL多基因风险评分（PRS）显示，欧洲血统比例高的个体遗传预测LTL更长。将LTL PRS加入回归模型后，欧洲血统与CH的关联仅轻度减弱（总体CH beta从0.80降至0.73），提示端粒长度仅部分解释祖先差异。

图 4　Fig. 4。图4展示MCPS中CH GWAS的曼哈顿图，突出TERT和TCL1B位点；rs187319135与多种CH亚型的关联方向（增加TET2-CH和SF3B1+SRSF2-CH风险，降低DNMT3A-CH风险）；以及TCL1B启动子罕见变异rs774615666与TET2-CH和SF3B1+SRSF2-CH的强关联。

亮点 4｜TCL1B位点祖先特异性变异

MCPS的CH GWAS发现TERT位点（rs2853677）与总体CH显著关联，与欧洲研究一致。此外，在TCL1B上游82.9kb和115kb处发现两个新的低频变异rs968294563和rs187319135（MAF约1%），与总体CH、TET2-CH和SF3B1+SRSF2-CH风险增加相关，但与DNMT3A-CH风险降低相关。rs187319135在印第安血统中富集（MAF=1.62%），在欧洲血统中罕见（MAF=0.06%），且与已知欧洲信号rs10131341和rs2887399连锁不平衡低，提示独立信号。

图 5　Fig. 5。图5展示MCPS与UKB跨祖先荟萃分析的曼哈顿图，分别对应总体CH、DNMT3A-CH、TET2-CH和ASXL1-CH，标注了5个新位点（MYCN、MEIS1、PARP11-CCND2、UBE2G1-SPNS3、PPIA-H2AZ2-PURB）。

亮点 5｜跨祖先荟萃分析发现五个新位点

MCPS与UKB的跨祖先GWAS荟萃分析鉴定出5个新的CH相关位点：总体CH的MYCN（rs12471506）、DNMT3A-CH的MEIS1（rs2280334）和PARP11-CCND2（rs582975）、TET2-CH的UBE2G1-SPNS3（rs73332852）、ASXL1-CH的PPIA-H2AZ2-PURB（rs13245012）。此外，UKB单独分析发现DNMT3A-CH的新位点PDXDC1。基因负荷分析在荟萃分析中复制了CHEK2与总体CH和DNMT3A-CH的关联，且CHEK2罕见变异组成在MCPS和UKB间存在差异。

生信可带走

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

• 方法栈： WES、WGS、GWAS、ExWAS、基因负荷分析、跨祖先荟萃分析
• 公开数据： 原文未给出公开组学登录号
• 代码： 原文未给出公开 GitHub/GitLab 仓库
• 谁该点开原文： 从事血液肿瘤、衰老或群体遗传学研究的生信人员，尤其是关注克隆性造血遗传基础的研究者。

带走一句

跨祖先比较能揭示CH的遗传异质性，非欧洲人群中的低频变异可能具有较大效应，分析时应纳入多祖先数据并关注祖先特异性信号。