The efficacy of antidepressant drug treatment in depression is unsatisfactory; 1 in 3 patients does not fully recover even after several treatment trials. Genetic factors and clinical characteristics contribute to the failure of a favorable treatment outcome.To identify genetic and clinical determinants of antidepressant drug treatment outcome in depression.Genomewide pharmacogenetic association study with 2 independent replication samples.We performed a genomewide association study in patients from the Munich Antidepressant Response Signature (MARS) project and in pooled DNA from an independent German replication sample. A set of 328 single-nucleotide polymorphisms highly related to outcome in both genomewide association studies was genotyped in a sample of the Sequenced Treatment Alternatives to Relieve Depression (STAR*D) study.A total of 339 inpatients with a depressive episode (MARS sample), a further 361 inpatients with depression (German replication sample), and 832 outpatients with major depression (STAR*D sample).We generated a multilocus genetic variable that described the individual number of alleles of the selected single nucleotide polymorphisms associated with beneficial treatment outcome in the MARS sample ("response" alleles) to evaluate additive genetic effects on antidepressant drug treatment outcome.Multilocus analysis revealed a significant contribution of a binary variable that categorized patients as carriers of a high vs low number of response alleles in the prediction of antidepressant drug treatment outcome in both samples (MARS and STAR*D). In addition, we observed that patients with a comorbid anxiety disorder combined with a low number of response alleles showed the least favorable outcome.These results demonstrate the importance of multiple genetic factors combined with clinical features in the prediction of antidepressant drug treatment outcome, which underscores the multifactorial nature of this trait.