Hyperactivity of the hypothalamic-pituitary-adrenal (HPA) axis is among the most consistent neuroendocrine abnormalities in major depressive disorder (MDD). The peptide adrenocorticotropin hormone (ACTH) mediates HPA axis function during stress and is encoded by the proopiomelanocortin (POMC) gene polycistronically. After screening 39 POMC polymorphisms, we evaluated the association of polymorphisms with susceptibility to MDD in 145 MDD patients and 193 normal subjects; in patients, we also evaluated the response to treatment with antidepressants. Additionally, we investigated the role of gene-environment interaction between POMC haplotypes and stressful life events (SLE) in the treatment response. Although genotypes and haplotypes were not significantly associated with the risk of MDD, non-remitters were more likely to carry haplotype 1 (ht1) and to have no ht2 than were remitters (corrected P = 0.010-0.035). Although observations were limited in patients without SLE, a significant haplotype-SLE interaction was observed (P = 0.020). Additionally, at 1, 2, and 8 weeks of treatment, the 21-item Hamilton Depression Rating scores of MDD subjects with POMC ht2 were significantly (P = 0.003-0.044) lower than those of patients with ht1 in subjects those did not experience SLE. MDD subjects possessing POMC ht2 achieved remission significantly (P = 0.013; survival analysis) faster than patients with ht1. This study suggests that POMC haplotypes, via an interaction with SLE, are associated with antidepressant treatment outcomes in MDD patients. Regarding SLE, haplotypes of the POMC gene could be useful markers for predicting the response to antidepressant treatment in MDD patients.