The cis-eQTL data for depression GWAS SNPs contained 428 significant SNP–gene interactions, with the top effects on GSDME (rs2721811, p = 6.07 × 10−123), PCDHA8 (rs3806843, p = 7.23 × 10−88), PCDHA13 (rs3806843, p = 1.45 × 10−71), PCDHA10 (rs3806843, p = 8.73 × 10−61) and ZSCAN32 (rs853679, p = 4.47 × 10−58) (Table 1). Here, PCDHA8 is linked to depressive disorder.