For these analyses there was no need for a stringent genetic relatedness cut-off as in the heritability estimation described earlier, but we did exclude one of each pair of individuals with a genetic relatedness larger than 0.3 so that twin and sibling pairs did not bias the p-values. This resulted in a sample of 10,247 individuals. Population structure (first 20 principle components) was corrected for before analysis.