We constructed a GRS for each of the three BP traits for each individual by summing the additive coding of each set of SNPs associated with the particular BP trait weighted by the previously-described effect size from ICBP (phs000585.v1.p1), and then standardized the distribution of all groups simultaneously by the mean and standard deviation (i.e., to a standard normal distribution) for interpretability. We used the leading SNP from each locus.