The genetic architecture of human DNA methylome variation

Population DNA methylome analyses remain largely confined to reference-defined genomic space despite extensive human genetic diversity. A central challenge is therefore to define comparable methylation phenotypes across diverse genomes and map the genetic architecture of methylome variation. Here, using 460 haplotypes from the Human Pangenome Reference Consortium Release 2, we develop two complementary pangenome-based frameworks: one compares methylation within homologous regions while retaining haplotype-specific structural configuration and CpG content; the other represents CpG island (CGI) variation at homologous loci as sequence-resolved alleles. We identified 47,039 CGI alleles, including 8,040 absent from CHM13; 62% were attributable to structural variants and 31% to small variants. Through local genealogies, we reconstructed CGI gain, loss and turnover, revealing how present-day methylation differences align with CGI evolutionary trajectories. We found that genome-wide, methylation variability was buffered against genetic diversity across more than 90% of the genome, with exceptions at extreme nucleotide diversity and locally around structural variants. Our methylation framework enabled population-scale allele-aware methylation quantitative trait locus (mQTL) mapping into non-reference and structurally complex sequence, where mQTL-associated intervals were enriched. mQTLs linked methylation to transcriptomic phenotypes, population differentiation and complex traits, supporting 8% of population differentially methylated regions and colocalizing with more than half of genome-wide association study loci. Together, our study redefines human methylome analysis and our findings reveal a genome-wide architecture of genetic-epigenetic coupling, reshaping how methylation variation is interpreted across human development and disease.