For differentially methylated region analysis, we applied a linear model on VMRs (see ‘VMR analysis’) as a function of: (1) global genetic ancestry; (2) local genetic ancestry; (3) sex; (4) age; and (5) top five principal components of DNA methylation derived from the top one million variable CpGs. We corrected both statistical tests for multiple testing using the Benjamini–Hochberg method.