These genes included Asxl1, Zfp292, Wrap73, Col5a3, Abcc5, Fscn1, Wdfy3, Pou6f2, Svil, Prlr, Gnat2, Slc30a7, Kdm5a, Slco1a2, Nlrp1a, Crlf3, Tpcn1, Pigr, Pou6f2 and Reep3. Among these genes, GNAT2 has a variant in human (rs6537837) that was reported to be associated with unipolar depression with a genome wide significance of p = 1e−646, while Prlr and Nlrp1a were implicated in depression-like behavior in animal models47,48. Here, WDFY3 is linked to major depressive disorder.