完整人类着丝粒的变异与进化

The variation and evolution of complete human centromeres

作者信息Glennis A Logsdon, Allison N Rozanski, Fedor Ryabov, Tamara Potapova, Valery A Shepelev, Claudia R Catacchio, David Porubsky, Yafei Mao, DongAhn Yoo, Mikko Rautiainen, Sergey Koren, Sergey Nurk, Julian K Lucas, Kendra Hoekzema, Katherine M Munson, Jennifer L Gerton, Adam M Phillippy, Mario Ventura, Ivan A Alexandrov, Evan E Eichler
PMID38570684
期刊Nature
发布时间2024-05
DOI10.1038/s41586-024-07278-3

实验完整度

包含两个独立证据层级:完整着丝粒组装与比较基因组分析,以及CENP-A ChIP-seq和DNA甲基化的功能验证,并涉及多物种比较分析。

主要模型

CHM1细胞系 CHM13细胞系 黑猩猩、猩猩、猕猴细胞系 56个人类基因组样本

重点核对

CHM1和CHM13细胞系的着丝粒组装完整性 CENP-A ChIP-seq和DNA甲基化数据 黑猩猩、猩猩、猕猴的着丝粒组装 56个人类基因组的着丝粒比对 系统发育分析中的比对策略和参数

摘要

人类着丝粒因其重复性和大规模的特性,传统上极难进行测序和组装。因此,尽管着丝粒是突变最快速的区域之一,人类着丝粒的变异模式和进化模型仍不完整。本研究利用长读长测序技术,对第二个人类基因组的所有着丝粒进行了完整测序和组装,并与已完成的参考基因组进行了比较。我们发现,相较于其独特的侧翼序列,这两组着丝粒的单核苷酸变异至少增加了4.1倍,且大小差异可达3倍。此外,由于新的α-卫星高阶重复序列的出现,45.8%的着丝粒序列无法使用标准方法进行可靠比对。DNA甲基化和CENP-A染色质免疫沉淀实验表明,26%的着丝粒的动粒位置差异超过500 kb。为了理解进化变化,我们选择了六条染色体,并从黑猩猩、猩猩和猕猴基因组中测序和组装了31个直系同源着丝粒。比较分析揭示了α-卫星HOR几乎完全的更替,每个物种都有特征性的α-卫星HOR特异性变化。人类单倍型的系统发育重建支持着丝粒的短臂和长臂之间几乎没有重组,并揭示了新的α-卫星HOR具有单系起源,为估计人类着丝粒DNA的跳跃式扩增和突变速率提供了策略。

实验结论

提炼研究问题、关键发现与证据,快速把握文章的核心贡献。

研究问题
人类着丝粒的遗传变异和进化机制是什么?
核心机制
着丝粒的快速变异和进化主要由α-卫星HOR的跳跃式扩增驱动,伴随着新的HOR结构的出现,且动粒位置可变。
主要证据
基于CHM1和CHM13完整着丝粒的比较,结合56个人类基因组和多种灵长类动物直系同源着丝粒的分析,以及CENP-A ChIP-seq和DNA甲基化实验。
研究意义
提供了人类着丝粒变异和进化的基线,并为理解着丝粒突变速率和动粒位置可塑性提供了框架。

研究路径

按研究推进顺序梳理实验设计、验证步骤与关键观察。

1

完成CHM1着丝粒的组装与验证

获得第二个人类基因组的完整着丝粒序列,并验证其准确性。

利用PacBio HiFi和Oxford Nanopore长读长测序数据,通过hifiasm组装和SUNK辅助scaffolding,最终生成完整的着丝粒组装,并通过多种方法验证其完整性。

2

比较CHM1与CHM13着丝粒的遗传变异

量化两个人类基因组之间着丝粒的序列和结构差异。

通过三种比对策略,比较CHM1和CHM13着丝粒的序列同一性、长度和结构变异。

3

分析着丝粒的表观遗传差异

确定着丝粒功能区域(动粒)的位置和大小差异。

通过CENP-A ChIP-seq和DNA甲基化分析,定位CHM1和CHM13的动粒位置,并比较其差异。

4

跨物种比较着丝粒进化

重建着丝粒在灵长类中的进化历史。

对黑猩猩、猩猩和猕猴的直系同源着丝粒进行测序和组装,比较其结构和序列,并构建系统发育树。

5

估算着丝粒突变速率并分析新HOR起源

估计着丝粒的突变速率和新HOR的起源时间。

基于人类和NHP的比对,估算侧翼区域突变速率;利用HPRC/HGSVC数据,构建系统发育树并估计新HOR的出现时间。

研究方法

按研究目的归类文中使用的方法,便于定位所需技术。

产品清单

实验环节名称品牌货号
AmnioMax C-100基础培养基Thermo Fisher Scientific17001082
AmnioMax C-100补充物Thermo Fisher Scientific12556015
青霉素-链霉素Thermo Fisher Scientific15140122
RPMI-1650培养基Sigma-AldrichR8758
胎牛血清Thermo Fisher Scientific16000-044
MEM α培养基Thermo Fisher Scientific12571063
Qiagen Gentra Puregene细胞试剂盒Qiagen--
核糖核酸酶AQiagen19101
蛋白酶KQiagen19131
Triton X-100----
模板制备试剂盒v1PacBio100-259-100
SMRTbell Express模板制备试剂盒v2PacBio100-938-900
SMRTbell酶清洁试剂盒PacBio101-746-400
SageELFSage ScienceELF001
FEMTO PulseAgilentM5330AA
Sequel II测序化学1.0PacBio101-717-200
Sequel II测序化学2.0PacBio101-820-200
SMRT Cell 8MPacBio101-389-001
超长DNA测序试剂盒ONTSQK-ULK001
Nanobind UL文库制备试剂盒CirculomicsNB-900-601-01
FLO-PRO002 R9.4.1流动池ONT--
连接测序试剂盒ONTSQK-LSK110
KaryoMAX秋水仙胺溶液Thermo Fisher Scientific15212012
含DAPI的Vectashield封片剂Vector LaboratoriesH-1200-10
抗人CENP-A抗体EnzoADI-KAM-CC006-E
Halt蛋白酶抑制剂混合物Thermo Fisher Scientific78429
微球菌核酸酶Thermo Fisher Scientific88216
蛋白G琼脂糖珠Abcamab193259
TruSeq ChIP文库制备试剂盒,Set AIlluminaIP-202-1012
抗CENP-C抗体----
山羊抗兔IgG-FITCSigma-AldrichF0382
Cy3-dUTPEnzo42501

关键环节

汇总复现实验时建议重点确认的条件及原文阅读提示。

环节核对要点
细胞培养
培养基成分(如AmnioMax C-100、RPMI-1650、MEM α),FBS比例,青霉素-链霉素浓度,培养温度(37°C)
阅读提示:Methods - Cell culture
基因组组装与验证
PacBio HiFi和ONT测序覆盖度,组装软件(hifiasm、SUNK、HiCanu)及其版本,验证方法(VerityMap、GAVISUNK、Merqury)
阅读提示:Methods - Targeted sequence assembly and validation
表观遗传分析
CENP-A ChIP-seq的抗体、细胞数、MNase浓度,甲基化分析中ONT读长阈值(>30 kb)和Nanopolish参数
阅读提示:Methods - Native CENP-A ChIP–seq and analysis; CpG methylation analysis
免疫荧光染色
CENP-C抗体稀释度(1:40),免疫FISH的杂交温度和时间
阅读提示:Methods - Immuno-FISH on stretched metaphase chromosome spreads
系统发育分析
比对方法(minimap2参数)、系统发育软件(IQ-TREE)和模型选择,分子钟设定(人-黑猩猩分歧时间)
阅读提示:Methods - Human-specific phylogenetic analysis; Mutation rate estimation