As results suggested that there was a significant interaction between the DVL3 (rs1709642, rs1969253) and GSK3β (rs334558, rs6438552, and rs2199503) polymorphisms on MDD susceptibility (P < 0.05), the interaction between DVL3 (rs1709642) and GSK3β (rs334558, rs6438552) showed a CV consistency of 10/10 and a testing accuracy of 55.31% which was considered as the best multilocus model. Here, GSK3B is linked to major depressive disorder.