Multidimensional variation and population stratification across 8000 complete human centromeres
Multidimensional variation and population stratification across 8000 complete human centromeres
Sun, Y.; Wan, S.; Nie, L.; Yu, D.; Zhou, F.; Yang, Y.; Yang, X.; Liu, A.; Chen, Q.; Fu, K.; Ni, Q.; He, Y.; Su, B.; Mao, Y.; Ye, K.; Yang, X.; Zhang, G.; Dongya, W.
AbstractHuman centromeres are indispensable for the faithful segregation of chromosomes during cell division, yet their highly repetitive nature has historically precluded comprehensive characterization, leaving fundamental questions about their sequence diversity, evolution trajectories and function dynamics unresolved. Here, we generated 6,312 complete human centromere sequences from 320 phased genome assemblies in Asian Pan-Genome project phase 1. By integrating the assemblies from the Human Pangenome Reference Consortium (HPRC) and Human Genome Structural Variation Consortium (HGSVC), we constructed a multidimensional genetic variation map encompassing over 8,000 gapless centromeres. Centromeric satellite arrays account for 4.19% to 6.01% of the whole genome, with substantial variations in size and architecture across chromosomes. Using a refined alpha satellite clustering approach that captures global diversity, we identified 195 higher-order repeat (HOR) arrays, 56.4% of which are absent from the T2T-CHM13 reference genome. Extensive structural variations across multiple dimensions exhibit population stratification, including centromeric haplotypes (CenHaps), ultra-large pericentric inversions spanning up to 36.6 Mbp, and inter-chromosomal HOR sharing that reflects sequence exchange among chromosomes. Integrating CENP-A CUT&Tag experiments and long-read-based DNA methylation profiles, we demonstrate that 16.8% of centromeres harbor multiple potential kinetochore assembly sites, and CenHap-specific local HOR homogenization is associated with kinetochore positioning. Despite global suppression of recombination at centromeres, we observed asymmetric linkage disequilibrium flanking centromeres and an ancient recombination event within the centromere of chromosome 19. Furthermore, contrary to the prevailing assumption of high mutation rates in centromeres, our estimates based on stringent orthology reveal no significantly higher single-base substitution rates for centromeres relative to flanking pericentromeric regions with substantial variations across chromosomes, despite extraordinary structural plasticity. Collectively, these multi-scale centromeric variations provide a global view of human centromere diversity and population stratification, fundamentally redefine centromere evolution through a dual-track model balancing structural innovation with mutational constraint, and establish an essential resource for investigating centromere biology and a baseline reference for diagnosing centromere-associated disorders.