Hamilton Depression Scale scores showed a positive correlation with IL-10 (β = 0.512, 95% CI: 0.128–0.896, p = 0.012), whereas SCOPA-AUT scores were inversely linked to IFN-γ (β = −0.588, 95% CI: −0.962 to −0.214, p = 0.003), with adjusted R2 values of 41.2 and 49.5%, respectively. This evidence concerns the gene IL10 and major depressive disorder.