We found three significant biomarkers (Χ2 = 31.78, df = 3, p < 0.001; Nagelkerke = 0.303) which were positively associated with major depression, i.e., sIL-1RA, TNF-R1, and TBARS levels. The gene discussed is TNFRSF1A; the disease is major depressive disorder.