Multivariate analysis identified PD-1 (OR = 3.32) and MMP-9 (OR = 2.18) as independent risk factors for depression, while IL-10 (OR = 0.62) and BDNF (OR = 0.12) emerged as protective factors. The gene discussed is BDNF; the disease is major depressive disorder.