For the analysis, we selected SNPs significantly associated to depression at p < 1 × 10−5. After clumping those SNPs using 500 kb flanking regions and an LD cutoff threshold r2 > 0.1, 63 independent SNP signals were identified from 816 top variants. These 63 top SNPs were further merged into 59 non-overlapping loci containing 157 genes, which were then assessed using the DEPICT algorithm for gene set and tissue enrichment18. Results shown in Supplementary Table 3 are not corrected for multiple testing.