For incident depression, removing TNF-α from multi-adjusted models decreased the size of the effect for anxiety (Model 1: 2.77 (1.65, 7.89) to Model 2: 1.82 (1.02, 2.71) and increased the observed estimate for storage LUTS (Model 1: 1.45 (1.05, 2.36) to Model 2: 1.98 (1.03, 4.51) and reduced the model fit (R2change: 0.244–0.220). This evidence concerns the gene TNF and major depressive disorder.