A multivariate linear model of all biomarkers showed that TNF-α only had a positive correlation with BASDEC depression score (p = 0.007). This evidence concerns the gene TNF and depressive disorder.