In the final regression model, only oxytocin level was a predictor (B = − 0.090, p < 0.0001), the model explaining 65.2% of the depression variation. The gene discussed is OXT; the disease is major depressive disorder.