Pangenome alignment reveals global diversity and evolution of human centromeric regions
Authors:
- Jordan EizengaJEJordan Eizenga
- Mira MastorasMMMira Mastoras
- Julian K. LucasJLJulian K. Lucas
- Julian MenendezJMJulian Menendez
- Faith OkamotoFOFaith Okamoto
- Glenn HickeyGHGlenn Hickey
- Sasha A. LangleySLSasha A. Langley
- Hailey LoucksHLHailey Loucks
- Fedor RyabovFRFedor Ryabov
- Yulia ZybinaYZYulia Zybina
- Mobin AsriMAMobin Asri
Preprint
Abstract
Centromeres play essential roles in chromosome segregation and genome stability, yet they remain among the least characterized regions of the human genome. Despite advances in long-read sequencing and complete genome assembly, the extreme repetitiveness and structural complexity of these regions still challenge population-scale analysis, obscuring their mutational dynamics. The Human Pangenome Reference Consortium has now accurately assembled over 6,000 centromeres, providing an opportunity to catalog global centromere variation. However, centromeric regions have been systematically excluded from pangenome alignments due to the technical challenge of aligning their highly repetitive tandem arrays and extreme structural variability. Here we introduce Centrolign, a graph-based multiple sequence alignment tool that combines a uniqueness-driven objective function with partial-order partial-order alignment to accurately align alpha satellite higher-order repeats. By prioritizing rare matches within tandem arrays and leveraging extended centromere-spanning haplotypes formed by suppressed recombination, Centrolign produces progressive multiple sequence alignments that preserve ancestral repeat organization. Applied across human centromeres, these alignments reveal the phylogenetic structure of similar satellite array haplotypes and enable precise estimation of variation rates, structural variant frequencies, and spatial patterns of mutation within satellite arrays. Integrating Centrolign graphs with repeat annotation tools and pangenome mapping algorithms allows accurate variant calling and genotyping from long reads without prior assembly. Moreover, we show that centromere haplotypes can be accurately subtyped with k-mers alone. Together, these advances establish a robust framework for incorporating centromeres into broader pangenomes, and population genomics in general, advancing our understanding of human genome evolution and diversity.Keywords
Funding Information
- NIH/NHGRI (R01HG011274)
- NIH/NHGRI (UM1HG010971)
- NIH/NHGRI (R01HG014490)
- NIH/NHGRI (U01HG010961)
- NIH/NHGRI (U24HG010262)
- NIH/NHGRI (OT2OD026682)
This publication has 82 references indexed in Scilit:
- SNV/indel hypermutator phenotype in biallelic RAD51C variant: Fanconi anemiaHuman Genetics, 2023
- Differential Relationship of DNA Replication Timing to Different Forms of Human Mutation and VariationThe American Journal of Human Genetics, 2012
- Alpha-satellite DNA of primates: old and new familiesChromosoma, 2001
- Linear-Time Longest-Common-Prefix Computation in Suffix Arrays and Its Applications2001
- Organization and evolution of an alpha satellite DNA subset shared by human chromosomes 13 and 21Journal of Molecular Evolution, 1993
- Color Set Size problem with applications to string matching1992
- The phylogeny of human chromosome specific alpha satellitesChromosoma, 1988
- Hierarchical order in chromosome-specific human alpha satellite DNATrends in Genetics, 1987
- An improved algorithm for matching biological sequencesJournal of Molecular Biology, 1982
- Comparison of phylogenetic treesMathematical Biosciences, 1981