As shown in Table 1, MR analysis demonstrated significant causal effects of depression-related exposures on target genes: Inverse-variance weighted (IVW) model estimates for “Antidepressant”: JUN: OR = 1.029, 95% CI = 1.002–1.057, p = 0.034; PTK2: OR = 1.040, 95% CI = 1.005–1.076, p = 0.027; IVW estimates for “Depression–dysthymia” diagnosis: JUN: OR = 1.048, 95% CI = 1.021–1.076, p = 0.0004; PTK2: OR = 1.073, 95% CI = 1.039–1.109, p = 1.808 × 10−5. Here, PTK2 is linked to depressive symptom measurement.