Univariate analysis using the Kruskal–Wallis one-way ANOVA on ranks yielded one conventionally statistically significant variable (MCP-1; p = 0.027) for self-reported depression, one conventionally significant variable (TNF-α; p = 0.013) for clinical evidence of depression, and one conventionally significant variable (G-CSF; p = 0.030) for history of WNND. Here, CSF3 is linked to major depressive disorder.