Genome-wide approaches to antidepressant treatment: working towards understanding and predicting response
© BioMed Central Ltd. 2012
Published: 27 June 2012
Skip to main content
We’re sorry, something doesn't seem to be working properly.
Please try refreshing the page. If that doesn't work, please contact us so we can address the problem.
© BioMed Central Ltd. 2012
Published: 27 June 2012
We’re sorry, something doesn't seem to be working properly.
Please try refreshing the page. If that doesn't work, please contact us so we can address the problem.
Antidepressants are among the most commonly prescribed drugs, and a range of medications are available. However, treatment response to a particular drug varies greatly between patients, with only 30% of patients responding well to the first treatment administered. Given evidence that antidepressant treatment response is a heritable trait, together with technological advances in genetic research, three recently published genome-wide investigations into antidepressant responses have examined the determinants of variability in treatment outcomes between depressed patients. Here, we review these studies within the context of wider research efforts to identify treatment response predictors. Some interesting genes have been implicated, but no variants have yet been robustly and reliably linked to response. This may suggest that genetic effect sizes are smaller than originally anticipated. Candidate gene approaches in these samples have lent support to the involvement of serotonergic, glutamatergic and stress-response systems in treatment response, although corroborative evidence from genome-wide analyses indicates these results should be interpreted cautiously. Closer examination of antidepressant response, considering it as a complex trait, has indicated that multiple genes of small effect are likely to be involved. Furthermore, there is some evidence that genetic influence on response to treatment may vary between patients with different symptom profiles or environmental exposures. This has implications for the translation of pharmacogenetic findings into clinical practice: genotypic information from multiple loci and data on non-genetic factors are likely to be needed to tailor antidepressant treatment to the individual patient.
Depression is a serious and prevalent psychiatric disorder, and, while there are a range of treatment options available, there is a high degree of variability between patients in terms of their response to a particular treatment. Genes are likely to play an important role in this variability, and with the rapid pace of technological development in the field of genetics there is a growing interest in using pharmacogenetic approaches to identify predictors of antidepressant response.
This review will focus on the three large genome-wide analyses of antidepressant response that have recently been published, and consider the findings within the context of wider research efforts to identify treatment response predictors. While genetic effect sizes appear to be smaller than originally anticipated, analyses considering possible interactions between both genetic and environmental factors, as well as methods that attempt to address the symptomatic heterogeneity of depression, may point the way to fruitful new research avenues for identifying clinically valuable predictors of treatment response.
Depression is a common and disabling illness with a lifetime prevalence of up to 17% . The World Health Organization projects that by 2020 depression will be the second leading contributor to the global burden of disease . The disorder is characterized by low mood, loss of interest and reduced energy. Depression is also associated with cognitive symptoms such as reduced concentration, low self-esteem and suicidal ideations, as well as somatic symptoms such as early morning wakening, and loss of appetite and libido. There is a relatively high degree of symptomatic heterogeneity between depressed patients, with some showing 'atypical' features such as increased sleep and appetite.
Symptoms and classification of depression
DSM-IV classification of depressive episode
ICD-10 classification of depressive episode
(1) Depressed mood for at least 2 weeks
Five or more symptoms, including (1) or (2)
Mild: four or more symptoms, including two of (1), (2) or (3)Moderate: six or more symptoms, including two of (1), (2) or (3)Severe: eight or more symptoms, including (1), (2) and (3)
(2) Loss of interest and enjoyment
(3) Increased fatigability
(4) Loss of confidence/self-esteema
(6) Suicidal thoughts or intent
(7) Reduced concentration/indecisiveness
(9) Sleep disturbance
(10) Altered appetite
Single episode or recurrent
Single episode or recurrent
With/without psychotic featuresb
With/without psychotic featuresb (severe depression only)
With/without catatonic features
With/without somatic symptoms
With/without atypical features
With/without postpartum onset
Nevertheless, definitions of depression should not be considered absolute or immutable; there is continued debate over the best way to understand and define the illness. Indeed there is an ongoing research effort to confirm whether depression is best considered as a homogeneous clinical entity, given the variability seen between patients in terms of symptoms, course of illness and treatment response.
In order to measure symptoms and establish treatment response over time, numerous assessment tools are available, including clinician-rated scales such as the Hamilton Depression Rating Scale (HAM-D) , as well as self-report alternatives such as the Beck Depression Inventory (BDI) . However, there are differences in the items included in each of these scale, as well as variability in the weightings given to certain symptoms, reflecting the lack of consensus regarding the core features of depression.
Commonly used antidepressant medications
Mode of action
Tricyclic antidepressants (TCAs)
Nortriptyline, amitriptyline, imipramine, dothiepin
Block norepinephrine and serotonin reuptake transporters; also antagonize serotoninergic, histaminergic, adrenergic and cholinergic receptors
Common side effects include dry mouth and constipation; cardiovascular issues can occur; fatal in overdose
Selective serotonin reuptake inhibitors (SSRIs)
Fluoxetine, sertraline, paroxetine, citalopram, escitalopram
Block serotonin reuptake transporters; more selective than TCAs
Common side effects include nausea, sexual dysfunction; lower cardiotoxicity than TCAs (therefore safer in overdose)
Serotonin-norepinephrine reuptake inhibitors (SNRIs)
Block norepinephrine and serotonin reuptake transporters; more selective than TCAs
Common side effects similar to SSRIs; also associated with increased blood pressure
TCAs act by blocking both norepinephrine and serotonin transporters. This inhibits reuptake of these two neurotransmitters into the presynaptic neuron, so increasing available levels in the synaptic cleft to bind to postsynaptic receptors. TCAs also act as antagonists at serotonergic, histaminergic, cholinergic and adrenergic receptors. However, these actions are associated with commonly reported side effects (such as dry mouth and constipation). Furthermore, in some cases TCAs can cause cardiovascular issues and be fatal in overdose.
Given these adverse effects, SSRIs were developed to target the serotonin transporter much more selectively. SSRIs are widely used in depression, and have replaced TCAs as the first line of treatment for depression due to their more limited side-effect profile and improved safety. Nevertheless, TCAs remain useful for patients who do not respond to newer antidepressants. The dual-action SNRIs have similar selectivity and side-effect profiles to SSRIs, but inhibit the reuptake of both serotonin and norepinephrine.
The antidepressants described above have comparable average efficacy . However, individuals vary widely in their response to treatments, with only an estimated 30% of patients responding sufficiently to their first treatment [10, 11]. Antidepressants have a delayed onset of action, and it can take up to 8 weeks to assess treatment response  and, prior to commencing treatment, it is unclear which drug will produce an optimal response for a particular patient. Given the time and costs required to identify the most effective medication, together with the negative impact of multiple treatment steps (which are associated with increased risk of relapse ), there is considerable clinical value in identifying predictors of treatment response.
While it is logistically difficult to collect the samples required to investigate whether antidepressant treatment response is a heritable trait, family studies looking at response to antidepressants indicate that variability in response to these medications is likely to be at least partially genetic in nature [12, 13]. Given this evidence, a growing number of studies are investigating which genes might determine antidepressant treatment response. This research has the dual aims of better understanding how these drugs exert their therapeutic action, and enabling individualized treatment recommendations.
The US-based Sequenced Treatment Alternatives to Relieve Depression (STAR*D) project  enrolled a total of 4,041 treatment-seeking adults. All patients included in the study had a diagnosis of non-psychotic major depressive disorder (as confirmed by DSM-IV criteria), scoring greater than 14 on the HAM-D (that is, suffered from moderate to severe depression). Among the sample, there were high rates of chronic or recurrent depression. STAR*D had a multistep design; all patients received citalopram (an SSRI) in the first treatment level, then those without sufficient symptomatic benefit moved onto further treatment levels where they received either alternative or augmented treatment options.
STAR*D is the largest and longest study to ever be conducted evaluating antidepressant treatments. It was found that with first-level treatment with citalopram, approximately 30% of patients achieved remission. While the study was not expressly designed for pharmacogenetic purposes, DNA was collected from 1,953 patients, allowing a GWAS of treatment outcomes .
Looking at those patients who were enrolled on the trial for at least 6 weeks, the investigators elected to explore two outcome phenotypes: response and remission to citalopram treatment. Treatment response was defined as a greater than 50% reduction in scores on the Quick Inventory of Depressive Symptomatology Self-Report version (QIDS-SR)  at the patient's final visit, giving 883 responders and 608 non-responders. Remission was defined as scoring 5 or less on QIDS-SR at follow-up, and patients achieving remission (n = 743) were compared against those classified as non-responders. Genotyping was performed on either the Affymetrix 500K or 5.0 Human SNP arrays, and 430,198 SNPs were used in the genome-wide analysis. The study was designed to be representative of the ancestral diversity within the US population and so stringent statistical controls for population stratification were required.
Three SNPs were associated with treatment response at a suggested level of significance (P <1 × 10-5), and these same SNPs were also among the top 25 hits associated with the remission phenotype (however, it should be noted that remitters represent a subgroup of responding patients). Rs6966038, which is located near the ubiquitin protein ligase E3C gene (UBE3C; involved in protein ubiquitination to promote degradation ), was associated with response at P = 4.65 × 10-7 and remission at P = 3.63 × 10-7. Rs6127921 is closest in location to the bone morphogenic protein 7 gene (BMP7; linked to bone development), and was associated with response at P = 3.45 × 10-6 and remission at P = 1.07 × 10-6. Finally, rs809736 is an intronic SNP in the RAR-related orphan receptor alpha gene (RORA; a nuclear receptor ) and was associated with response at P = 8.19 × 10-6 and remission at P = 7.64 × 10-5. There was no evidence that treatment side effects might confound the above findings.
The Munich Antidepressant Response Signatures (MARS) study is a German pharmacogenomic project, with genome-wide data available for 339 inpatients with either major depression or bipolar disorder (11.2%) . The majority of patients were German (85.1%), and all reported white European ancestry. The study was naturalistic in design; patients were treated with a range of antidepressants. Depression symptoms were measured using the HAM-D  and three outcome phenotypes were considered: early partial response (HAM-D score improvement greater than 25% after 2 weeks), response (greater than 50% improvement in scores after 5 weeks) and remission (a score of <10 on the HAM-D after 5 weeks or before hospital discharge). Genotyping was performed using the Sentrix Human-1 and HumanHap300 arrays. While MARS was a smaller study than STAR*D, the genome-wide analysis not only used a replication sample of 361 German inpatients, but also a subsample of 832 white outpatients from the STAR*D study to further validate findings.
The most significant association (P = 7.6 × 10-7) was seen between early partial response and rs6989467, a SNP in the 5' region of the cadherin 17 gene (CDH17; encoding a calcium-dependent membrane-associated glycoprotein). Looking across all three phenotypes, the SNP most consistently associated with treatment outcome was rs1502174 (P = 8.5 × 10-5), located downstream of EPHB1 (encoding ephrin type-B receptor 1 protein, a transmembrane protein involved in developmental processes in the nervous system).
The Genome Based Therapeutic Drugs for Depression (GENDEP) study is a European, multicenter, 12-week study with two active pharmacological treatment arms (for further details see ). A total of 811 outpatients with moderate to severe depression received one of two antidepressants, nortriptyline (a TCA) or escitalopram (an SSRI), in a partially randomized design. The primary outcome measure was the clinician-rated Montgomery-Asberg Depression Rating scale (MADRS) , on which each patient was measured weekly. For the pharmacogenomic analysis , treatment response was defined using percentage change in the MADRS from baseline to week 12 of the study and genome-wide data were available (using the Illumnia Human610 Quad chip) for 706 patients, 394 of whom were taking escitalopram and 312 nortriptyline. In a whole-sample analysis, associations at a suggestive level of significance were seen in intergenic regions on chromosome 1 (rs2136093, rs6701608, rs2136094) and chromosome 10 (rs16920624, rs11598854, rs7081156). However, when drug-specific analyses were performed, a genome-wide significant association with nortriptyline response was seen with an intronic SNP (rs 2500535) contained within the UST gene. UST encodes uronyl-2-sulphotransferase, an enzyme with an important link to neurotrophic mechanisms . This SNP also reached suggested levels of significance when testing gene-drug interactions (that is, the SNP differentially predicted treatment outcome with escitalopram and nortriptyline). The best performing SNP within the escitalopram-specific analysis was rs1126757 (P = 2.83 × 10-6), a synonymous SNP within the coding region of the IL11 (interleukin 11) gene. Interleukins are neuropoietic cytokines involved in inflammatory responses, and inflammation processes have been associated with depression .
Genome-wide association studies of antidepressant treatment response
GWAS sample size (N)
Primary outcome measure
Response phenotype definition
SNPs reaching suggestive levels of significance (nearest gene)
Affymetrix 500K/5.0 human SNP arrays
Response: >50% reduction in score
rs6966038 (UBE3C), rs6127921 (BMP7), rs809736 (RORA)
Remission: score <5
339 (two replication samples, N=361 and N=831)
Sentrix Human-1 and HumanHap300 BeadChip
Early partial response: >25% reduction in score after 2 weeks of treatment
rs6989467 (CDH17)a, rs1502174 (EPHB1)
Response: >50% reduction in score
Remission: score <10
Illumina 610 array
Percentage change from baseline to week 12
rs2136093, rs6701608 and rs2136094 (intergenic, on chromosome 1)rs16920624, rs11598854 and rs7081156 (intergenic, on chromosome 10)
Candidate gene studies of antidepressant treatment response undertaken in large cohorts
SLC6A4 (solute carrier family 6 (neurotransmitter transporter, serotonin), member 4)
No (but see )
HTR2A (5-hydroxytryptamine (serotonin) receptor 2A, G-protein-coupled)
rs7997012 (also see )
451 (plus replication sample)
rs7997012: reverse allelic association to 
FKBP5 (FK506 binding protein 5)
294 (plus replication sample)
rs1360780, rs3800373, rs4713916
Escitalopram or nortriptyline
GRIK4 (glutamate receptor, ionotropic, kainate 4)
The serotonin transporter is a key target for many antidepressants (including SSRIs, SNRIs and TCAs), and so the gene encoding this protein (SLC6A4) represents a prominent candidate gene for pharmacogenetics. In the promoter region of this gene there is a repeat polymorphism comprising a 44 bp deletion , known as 5HTTLPR (serotonin transport-linked polymorphic region). This variant was first associated with treatment response by Smeraldi et al. , a finding supported by a meta-analysis in 2007 . However, when 5HTTLPR was examined within large cohort samples, the results were mixed. Within the GENDEP sample, the polymorphism did indeed predict treatment outcome, but only for those patients on escitalopram . Sex-specific analysis suggested that this effect was concentrated in males. Furthermore, an additional variant (rs2020933) within the SLC6A4 gene also demonstrated a significant effect on treatment outcome.
However, other large studies report conflicting results. The Genetic and Clinical Predictors of Treatment Response in Depression (GenPod) project is a pharmacogenetic study based in the UK, where genotypic and outcome data are available for 520 patients, randomized to receive citalopram or reboxetine . Reports from GenPod as well as those from STAR*D  found no association between 5HTTLPR and treatment outcome. Inconsistent findings could result from ethnic variability  or differences in outcome response definitions, but in the recent meta-analysis by Taylor et al. , including both GENDEP and STAR*D data (although not GenPod), a weak effect of 5HTTLPR genotype on remission did not reach significance once potential publication bias was accounted for.
Another serotonergic gene that has been implicated in antidepressant response is HTR2A. This encodes the serotonin receptor 2A, which is downregulated by antidepressant treatment in parallel with clinical improvement . A significant association between a marker in the HTR2A gene and treatment outcomes was observed when looking at 68 candidate genes (768 SNPs) in the STAR*D cohort (n = 1,297) . A significant association was seen between rs7997012 and treatment response (the association with remission was not significant), and other markers within this gene also reached suggested levels of significance. Given differences in allele frequencies, the authors suggested that HTR2A may play a role in the racial differences in treatment response found in their mixed-ancestry sample. The top two SNPs implicated in STAR*D were also followed up in the MARS cohort . In a combined sample of 637 inpatients (as in the genome-wide analysis, a second independent German sample was used in addition to the 451 patients recruited as part of the MARS project) a significant association was seen between rs7997012 and both remission and response phenotypes. In that report, however, the allele association was the inverse of that observed in STAR*D; in MARS the AA genotype was associated with impaired treatment response whereas in STAR*D the G allele was associated with impaired treatment response.
When HTR2A was examined in the GENDEP cohort, as part of an analysis looking at 10 different candidate genes (n = 760), once again an association was seen with treatment outcomes, but there was no replication of the specific allelic relationship seen in STAR*D . In an escitalopram-specific analysis, no association was observed with the rs7997012 SNP. Instead, two intronic markers in the HTR2A gene (rs931623 and rs2224721, which are in linkage disequilibrium) were both significantly associated with response to the SSRI. Therefore, while several studies indicated a relationship between genetic variation in the HTR2A gene and treatment response, there has been no replicated evidence that specific variants confer increased likelihood of favorable outcomes. Further attempts at replication in large cohorts are needed to clarify this issue, and it may also be valuable to perform sequencing in this region in an attempt to identify the causal variant that might underlie these findings. Interestingly, Murphy et al.  have implicated genetic variation in HTR2A in side effects to the SSRI paroxetine. Thus, the potentially confounding effects of treatment discontinuation should also be considered.
In addition to serotonergic pathways, the hypothalamic-pituitary-adrenal (HPA) axis has also been implicated in antidepressant treatment response . When looking across eight genes in the HPA axis in the MARS sample , significant associations were observed between treatment outcomes and three SNPs within the glucocorticoid-receptor-associated co-chaperone gene FKBP5 (rs1360780, rs3800373 and rs4713916). Follow-up work demonstrated that these SNPs were linked to increased intracellular expression of FKBP5, which leads to regulation of the HPA axis via adaptive changes in the glucocorticoid receptor. This evidence supports the proposal that effective antidepressant treatment may involve modulation of the HPA axis. Examination of these three SNPs in STAR*D replicated the association between rs4713916 and remission , and it was noted that the linkage disequilibrium patterns observed in the mixed-ancestry sample of STAR*D could suggest that the causal variant underlying this relationship may lie in the promoter region of the gene. However, the association with FKBP5 was not replicated in GENDEP . As part of a ten-candidate gene analysis, the HPA axis genes FKBP5 and NR3C1 (which encodes the glucocorticoid receptor) were both examined. In a whole-sample analysis, no associations between treatment response and any of the markers tested in the FKBP5 gene (including the rs1360780 SNP implicated by Binder et al. ) were identified, although two SNPs within NR3C1 were associated at suggested levels of significance (rs852977 and rs10482633).
The final gene implicated in these large cohort studies of antidepressant treatment response is GRIK4, a gene encoding glutamatergic receptor, kainate 4. Paddock et al.  examined the 68 candidate genes previously considered in the STAR*D sample , but were able to incorporate genetic data from a second wave of genotyping, taking the sample size from 1,297 to 1,816. The authors not only replicated the HTR2A signal outlined above , but also reported evidence implicating rs1954787 within GRIK4 in this enlarged sample. The glutamatergic system has previously been considered within treatment response pathways  and this association could lend weight to this hypothesis, but, as highlighted, further replication is crucial. It was also observed that the combined effect of GRIK4 and HTR2A on treatment outcomes was more robust than when each gene was considered separately.
Despite the evidence above linking candidate genes and treatment outcome, if these findings were true, it would be expected that there would be a consistent pattern of associations seen across GWASs. Nevertheless, in both the STAR*D and GENDEP genome-wide analyses [15, 16], a priori candidate genes did not show associations with treatment outcomes above the levels expected by chance (although the STAR*D analysis does draw attention to 'nominally significant' associations with several SNPs of interest in the HTR2A gene). Although meta-analyses are needed, the results so far indicate that the current genes of interest do not exert a sizable direct effect on antidepressant treatment response.
It has been calculated that STAR*D and GENDEP are powered to detect single SNP main effects on general treatment response that are clinically significant in size . However, it now seems likely that treatment response is a complex trait that is not determined by a single genetic variant, but results from multiple interacting genetic effects. In which case, the failure of pharmacogenetic and pharmacogenomic approaches to identify any replicated genetic variants associated with clinically meaningful effects on treatment response may be because studies are underpowered to detect these effects.
Effects at one gene locus may not only be moderated by other genetic variants, but may also be dependent on environmental exposures. Environments such as stressful life events (SLEs) have been associated with depression , and investigations into this relationship have indicated that 5HTTLPR moderates the effect of SLEs on the emergence of depressive symptoms [48, 49]. Mandelli et al.  reported initial evidence that this same interaction may also predict antidepressant treatment response, in a sample of 159 mood-disorder patients treated with an SSRI. This interaction between SLEs and 5HTTLPR was also replicated in GENDEP, although only among patients taking escitalopram, not among those taking nortriptyline  (Figure 2). This could suggest that the gene-environment interaction is drug-specific, and may be why no statistically significant interactions were found with 5HTTLPR in a sample of 290 patients who were receiving a range of different antidepressants . Thus, it is likely that the accurate prediction of antidepressant treatment response will require the consideration of multiple factors, both genetic and environmental, and the interaction between them.
Beyond SLEs, other non-genetic factors that may influence antidepressant response, such as symptomatic differences between patients, have also been explored as potential indicators of antidepressant response. This reflects the ongoing debate as to whether depression is a single, homogeneous clinical entity. Three commonly considered subgroups, defined on their symptomatic presentation, are: melancholic depression (akin to 'endogenous' depression, with pervasive anhedonia, sleep and appetite disturbance plus diurnal variation of mood); atypical depression (with features including mood reactivity, hyperphagia and hypersomnia); and anxious depression (reflecting the high rates of anxious symptoms and co-morbid anxiety disorders seen in depression [53, 54]). In the STAR*D sample, differences in response to citalopram seen between melancholic and non-melancholic and between atypical and non-atypical depression were accounted for by demographic and clinical factors at baseline [55, 56], but anxious depression was linked to poorer treatment response . However, this pattern with anxious depression is not consistently replicated [58, 59]. Studies comparing the efficacy of different antidepressant drugs in these subtypes have also been inconsistent in their findings [60–63].
As noted by Uher et al. , it is important to consider not only whether associations between depression subtype and treatment response can be replicated, but also whether the observed effect sizes are sufficiently large to be of clinical value.
Interestingly, more robust predictors of treatment response were identified when the symptomatic variation in depression was explored in terms of dimensional traits (rather than the above categorical subtypes). Using detailed symptomatic data available in GENDEP, six dimensional facets of depressive symptoms were described: mood, anxiety, pessimism, interest activity (describing symptoms of indecisiveness and reduced activity), sleep and appetite . The interest-activity facet was the most robust predictor of response in GENDEP; symptoms of low interest and decreased activity predicted significantly worse treatment outcomes . Importantly, the association between interest-activity symptoms and treatment outcome was also replicated in the STAR*D sample.
The predictive value of the interest-activity facet is substantial enough to be of clinical importance (one additional remission could be accurately predicted for every three or four patients where interest activity was measured), and obtaining the required information is relatively inexpensive. Thus, interest-activity scores appear to be valuable in translational terms. Nonetheless, low levels of interest and activity predicted poorer outcomes for patients, regardless of which antidepressant was used. As Simon and Perlis  highlight, the greatest clinical utility is obtained from factors that differentially predict response to alternative treatments, as these may be used to guide treatment recommendations. Consequently, assessment of whether other treatment options might be effective in patients displaying loss of interest and decreased activity would be of translational value. Therefore, in addition to between-drug comparisons, work in 'therapygenetics'  (where predictors of response to talking therapies are considered) points to a potentially fruitful direction for exploring differential predictors of treatment response. If studies allowing direction comparisons between psychological and pharmacological treatments could be undertaken, response predictors with high translational value might be revealed. For example, there is preliminary evidence that the 5HTTLPR variant linked with poorer antidepressant response is associated with better response to cognitive behavioral therapy .
The observed heterogeneity in treatment response associated with symptomatic differences also points to interesting new avenues for pharmacogenetic research, as it may indicate that genetic moderation of response to antidepressants could be symptom-specific. For example, among patients taking nortriptyline in GENDEP, genetic variation in the beta-3 subunit of the G-protein complex (GNB3) was associated specifically with differences in improvement of neurovegetative symptoms (that is, symptoms relating to sleep and appetite), but was not associated with improvement in other symptom dimensions (such as observed mood or cognitive symptoms) . Therefore, consideration of the symptomatic diversity encompassed by the diagnostic category of depression (possibly through dimensional rather than categorical models) may allow the identification of genetic effects on specific components of treatment response that are not otherwise seen.
Genetic research into antidepressant response not only aims to identify outcome predictors, but also provides valuable insight into the physiological processes underlying successful treatment of depression. The lack of strong genetic signals in either genome-wide or candidate gene studies indicates that antidepressant treatment response is not a phenotype determined by a small number of highly influential genes. Instead, the genetic evidence suggests that antidepressant treatment response is best considered a complex phenotype, determined by multiple genes of small effect. While candidate gene findings require further replication and confirmation, the evidence suggests that these genes lie in multiple, interacting biological systems.
The serotonergic system, as the canonical target for antidepressants, is a key candidate for investigators, and has been implicated with associations observed between treatment response and variation in SLC6A4  and HTR2A [36–38], although effect sizes are small when considered alone. The effects of successful antidepressant treatment on the HPA axis suggest that this system also plays an important role in the physiological processes that lead to remission . This proposal is supported by the observed association between FKBP5 (which acts as a regulator within this system) and treatment response [41, 42]. Indeed, there has been a large research effort in attempting to develop novel antidepressants targeting the HPA system , although this is yet to achieve success . Similarly, evidence of glutamatergic involvement in both depression and treatment response  has led to investigations of ketamine as another potential new treatment for depression , and the association between variability in the glutamatergic gene GRIK4 and treatment response strengthens the hypothesis that this system also plays an important role in antidepressant action.
Arguably, what is most interesting about the genetic findings is that they provide preliminary indications that, rather than purely additive effects of each of the different pathways contributing to treatment response, the serotonergic, glutamatergic and HPA systems may have interactive effects. This might indicate a physiological interaction, as experimental observations in humans seem to suggest (for example, ).
The non-genetic predictors of response outlined here also may give an insight into the mechanisms underlying antidepressant response. For example, the finding that interest-activity symptoms (of low interest, reduced activity, indecisiveness and lack of enjoyment) predict poorer treatment response  could be interpreted to suggest that the pathophysiological processes producing these symptoms are, to some extent, distinct, and less well targeted by antidepressants.
Nevertheless, the identification of treatment response predictors is only the first step in understanding the actions of antidepressants; further exploration of the biological impact of these genetic variants and the mechanism by which suggested predictors may be associated with treatment response is needed, as well as a greater understanding of how these predictors may interact.
Genome-wide analyses of antidepressant treatment response have proved disappointing to date, with no strong signals emerging from the three studies conducted thus far. Given the small genetic effect sizes observed in candidate gene studies, it seems that existing studies are likely to be underpowered on a genome-wide scale, and strong conclusions cannot be drawn without increasing sample sizes and undertaking formal meta-analyses.
Nonetheless, there is some evidence from candidate gene approaches within these large samples indicating that genetic variability in the serotonergic, glutamatergic and stress-response pathways may influence antidepressant treatment response. These associations remain small in effect size; no genes of interest have yet been reliably linked with clinically significant effects. This suggests that, rather than single gene effects, there may instead be multiple genetic influences that determine treatment outcomes in combination with environmental factors. The interactions between each of these various elements must be considered if we wish to better provide treatment recommendations that are optimized for each individual.
Indeed, the best-performing predictor identified in the three large GWAS samples has been the symptom dimension of interest activity, and there is some evidence for symptom-specific improvements. This indicates that the future for pharmacogenomic studies lies not only in larger samples, but also in the collection of detailed and accurate patient data, to allow the identification of any subgroup-specific relationships.
The translational value of pharmacogenetic work with antidepressants has yet to be fully realized; indeed, for real clinical utility, studies must focus on comparing different treatments for depression in order to identify differential predictors of response. Furthermore, it is likely that translation into clinics will only be achieved by combining information from multiple response predictors. Nonetheless, the candidate genes identified to date provide an insight into the mechanisms of antidepressant action, and may enable the development of novel treatments for depression.
Beck Depression Inventory
Diagnostic and Statistical Manual of Mental Disorders: 4th edition
Genome Based Therapeutic Drugs for Depression
Genetic and Clinical Predictors of Treatment response in Depression
genome-wide association study
Hamilton Depression Rating Scale
International Classification of Diseases: 10th revision
Montgomery-Asberg Depression Rating scale
Munich Antidepressant Response Signatures
Quick Inventory of Depressive Symptomatology: self-report version
stressful life event
single nucleotide polymorphism
serotonin-norepinephrine reuptake inhibitor
selective serotonin reuptake inhibitor
Sequenced Treatment Alternatives to Relieve Depression
serotonin transport-linked polymorphic region.
Karen Hodgson is supported by an MRC PhD scholarship. This work was in part supported by the South London and Maudsley NHS Foundation Trust NIHR Biomedical Research Centre.