A genomewide association study of citalopram response in major depressive disorder

Antidepressant response is likely influenced by genetic constitution, but the actual genes involved have yet to be determined. We have carried out a genomewide association study to determine whether common DNA variation influences antidepressant response.

Our sample is derived from Level 1 participants in the Sequenced Treatment Alternatives to Relieve Depression (STAR*D) study, all treated with citalopram. Association for the response phenotype included 883 responders and 608 nonresponders. For the remission phenotype, 743 subjects that achieved remission were compared with 608 nonresponders. We used a subset of single nucleotide polymorphisms (SNPs; n = 430,198) from the Affymetrix 500K and 5.0 Human SNP Arrays, and association analysis was carried out after correcting for population stratification.

We identified three SNPs associated with response with p values less than 1 x 10(-5) near the UBE3C gene (rs6966038, p = 4.65 x 10(-7)), another 100 kb away from BMP7 (rs6127921, p = 3.45 x 10(-6)), and a third that is intronic in the RORA gene (rs809736, p = 8.19 x 10(-6)). These same SNPs were also associated with remission. Thirty-nine additional SNPs are of interest with p values < or = .0001 for the response and remission phenotypes.

Although the findings reported here do not meet a genomewide threshold for significance, the regions identified from this study provide targets for independent replication and novel pathways to investigate mechanisms of antidepressant response. This study was not placebo controlled, making it possible that we are also observing associations to nonspecific aspects of drug treatment of depression.