In a fully adjusted model, depression was significantly associated with higher levels of IFN-γ and IL-1ra; IL-33, CCL-2, CCL-4 and TNF- α remained significant below the Bonferroni threshold. The gene discussed is IL33; the disease is major depressive disorder.