Using MAGMA23, we tested the association between all 13 neuroticism phenotypes and 7,297 gene-sets. Specifically, we tested for association with 7,244 gene ontology gene-sets (MsigDB version 6.0; http://software.broadinstitute.org/gsea/msigdb/collections.jsp) and 53 tissue expression profiles (https://www.gtexportal.org/home/). For all gene-sets competitive P values were computed, which result from the test whether the combined effect of genes in a gene-set is significantly larger than the combined effect of a same number of randomly selected genes (in contrast, self-contained P values result from testing against the null hypothesis of no effect). We only report competitive P values, which are more conservative compared to self-contained P values. Competitive P values were Bonferroni corrected (α = 0.05/7,297 = 6.85 × 10−6).