Information about replication cohorts and their respective settings for replication analyses is outlined in the Supplementary Note. If applicable, summary statistics from different replication datasets for the same cell type or tissue were meta-analyzed using a weighted Z-score method61. Benjamini-Hochberg FDR17 was used to adjust replication analysis P-values for multiple testing. We required FDR<0.05 and the same effect direction with discovery to declare effect replicating. R package pwr (https://cran.r-project.org/web/packages/pwr/index.html) was used to conduct power analyses for replication datasets.