LD score regression (LDSC16,17) analyses were run on summary statistics obtained from GWA analyses to estimate the proportion of phenotypic variance attributable to all SNPs in the analyses (SNP-based heritability; h2SNP). Bivariate LD Score regression (https://github.com/bulik/ldsc) on the GWAS summary statistics was used to estimate the genetic correlations between all neuroticism phenotypes. Significance of genetic correlations between all neuroticism phenotypes was determined by applying a Bonferroni corrected threshold of P = 2.22 × 10−4 (0.05/(15 × 15) = 0.05/225; Fig. 1 and Supplementary Figure 24 combined).