SNPs from the summary statistics, PRS were coded for every individual by multiplying an individual’s number of effect alleles at a particular SNP by that SNP’s effect size (beta) from the discovery GWAS, then averaging across SNPs to create one score per person. Clumping was done using the European subset of the 1000 Genomes Phase 3 sample (Consortium, 2015) as an external linkage disequilibrium (LD) reference panel, using a 500 kb physical distance and an LD threshold of r2>=0.25. Scores representing effect sizes with increasingly lenient thresholds of statistical significance in the discovery GWAS were constructed (pT <0.0001, pT <0.001, pT <0.01, pT <0.05, pT <0.10, pT <0.20, pT <0.30, pT <0.40, pT <0.50). Scores were standardized before statistical analysis.