The individuals were genotyped using the Illumina OmniExpress 700K chip, imputed to the 1000 genome imputation panel. After imputation we selected SNPs from the HapMap3 panel since this SNP set was optimized to capture common genetic variation,21 (link) which is required for the GREML analysis. For quality control, SNPs with a minor allele frequency less than 1%, a higher missing rate than 0.03 and that failed the Hardy-Weinberg equilibrium for a threshold of 10-6 were removed. The first 20 principal components are included as covariates to adjust for population stratification.22 (link)