Moderation of antidepressant response by the serotonin transporter gene

Patricia Huezo-Diaz , Rudolf Uher , Rebecca Smith , Marcella Rietschel , Neven Henigsberg , Andrej Marušič , Ole Mors , Wolfgang Maier , Joanna Hauser , Daniel Souery , Anna Placentino , Astrid Zobel , Erik Roj Larsen , Piotr M. Czerski , Bhanu Gupta , Farzana Hoda , Nader Perroud , Anne Farmer , Ian Craig , Katherine J. Aitchison , Peter McGuffin

Abstract

Background

There have been conflicting reports on whether the length polymorphism in the promoter of the serotonin transporter gene (5-HTTLPR) moderates the antidepressant effects of selective serotonin reuptake inhibitors (SSRIs). We hypothesised that the pharmacogenetic effect of 5-HTTLPR is modulated by gender, age and other variants in the serotonin transporter gene.

Aims

To test the hypothesis that the 5-HTTLPR differently influences response to escitalopram (an SSRI) compared with nortriptyline (a noradrenaline reuptake inhibitor).

Method

The 5-HTTLPR and 13 additional markers across the serotonin transporter gene were genotyped in 795 adults with moderate-to-severe depression treated with escitalopram or nortriptyline in the Genome Based Therapeutic Drugs for Depression (GENDEP) project.

Results

The 5-HTTLPR moderated the response to escitalopram, with long-allele carriers improving more than short-allele homozygotes. A significant three-way interaction between 5-HTTLPR, drug and gender indicated that the effect was concentrated in males treated with escitalopram. The single-nucleotide polymorphism rs2020933 also influenced outcome.

Conclusions

The effect of 5-HTTLPR on antidepressant response is SSRI specific conditional on gender and modulated by another polymorphism at the 5′ end of the serotonin transporter gene.

The selective serotonin reuptake inhibitor (SSRI) antidepressants are used as first-line treatment of depression, but up to 50% of people with depression do not adequately benefit from SSRI treatment.1,2 A significant proportion of the individual differences in response to SSRI is familial.3 When searching for a genetic predictor of SSRI response, the serotonin transporter gene (SLC6A4) is the natural candidate as it encodes the molecular target of SSRI antidepressants. Several functional polymorphisms have been found in the SLC6A4 gene. The common length polymorphism constituted by a deletion of 44 base pairs in the promoter region, referred to as 5-HTTLPR, has been explored in most detail and implicated in depression aetiology.4 The long allele has been associated with more efficient transcription,5 resilience to depressogenic effects of adversity4 and better response to SSRI antidepressants.6 A meta-analysis of 14 studies with 1435 participants showed that long-allele carriers had higher probability of response and remission than short-allele homozygotes.7 Therefore, SLC6A4 was the prime candidate for response to the SSRI escitalopram and the first gene to be genotyped and tested in the Genome Based Therapeutic Drugs for Depression (GENDEP) project. However, significant heterogeneity of 5-HTTLPR effect was reported in non-European populations.8,7 Whereas no effect of 5-HTTLPR on efficacy of citalopram was found in a large mixed-ethnicity sample,9 an ethnicity-sensitive re-analysis of the same sample found an effect in the expected direction among White non-Hispanic participants.10 As the 5-HTTLPR allele frequencies and its linkage disequilibrium with other polymorphisms within the SLC6A4 gene vary markedly across populations,11 it is possible that the discrepancies are a result of the influence of another polymorphism that is in population-specific linkage disequilibrium with the 5-HTTLPR.

Three polymorphisms have been proposed to modulate the functionality of 5-HTTLPR: single nucleotide polymorphisms rs25531 and rs2020933, and a 17bp variable tandem repeat in the second intron (STin2). The rs25531 is located within the repeats constituting the 5-HTTLPR. Two initial reports disagreed about the location of this putatively functional A/G single nucleotide polymorphism being either within12 or immediately upstream13 of the 5-HTTLPR insertion/deletion segment. However, sequence alignment has demonstrated that these two reports point to the same polymorphism, the rs25531 (A/G), which is located in the sixth repeat of the 5-HTTLPR outside the deleted segment and therefore can co-occur with either long or short variant of the 5-HTTLPR.14 The minor G allele of the rs23551 creates a consensus binding sequence for the transcription factor AP-2, which is involved in antidepressant action.15,16 It has been proposed that the rs25531 influences the functionality of the serotonin transporter and recoding of 5-HTTLPR long alleles based on the rs25531 has been suggested.17 The functionality of this single nucleotide polymorphism has however not been consistently replicated and it is still unclear whether the recoding of the 5-HTTLPR is warranted.4,14,18 Therefore, in the present study, the rs25531 is evaluated separately from the 5-HTTLPR. The rs2020933 single nucleotide polymorphism is located in the first intron, approximately 2.5kb 3′ to the 5-HTTLPR and contributes to the variation of serotonin transporter expression.18 The STin2 two major alleles (STin2.10 and STin2.12) correspond to 10 or 12 repeats of a 17bp sequence and the STin2.12 allele has been associated with enhanced transcription.19

It has been proposed that the moderation of antidepressant response by the 5-HTTLPR may depend on age.20 This suggestion is supported by data showing age-related variation in serotonergic function,2123 but tests of interaction between age and 5-HTTLPR have not been reported in most pharmacogenetic studies.

Gender may be another important factor influencing the relationship between 5-HTTLPR genotype and antidepressant response. Ovarian steroids have a strong influence on serotonin synthesis,24 expression of serotonergic receptors,2527 and the serotonin transporter.28,29 There are gender differences in 5-HTTLPR moderation of response to serotonergic challenges, including tryptophan depletion30 and to stressful life events.4,31 However, previous pharmacogenetic studies have not reported a test of gender×genotype interaction.

The GENDEP project aims to identify differential predictors of response to SSRI and tricyclic antidepressants. The primary aim is to test the hypothesis that the 5-HTTLPR interacts with the drug used (escitalopram v. nortriptyline) in treatment of depression. Specifically, it is expected that short-allele homozygotic status will be associated with poor response to escitalopram but will not influence response to nortriptyline. As the heterogeneity of findings across studies may be a result of the influence of a polymorphism that is in population-specific linkage disequilibrium with the 5-HTTLPR, we systematically explore DNA sequence variation across the SLC6A4 gene, including the three polymorphisms with known effects on gene expression. Additionally, the moderating effects of age and gender are tested.

Method

Study design and sample

The GENDEP project is an open-label part-randomised multi-centre pharmacogenetic study with two active pharmacological treatment arms.32 It was designed to establish clinical and genetic determinants of therapeutic response to two antidepressants with contrasting primary modes of action: nortriptyline (a tricyclic inhibitor of monoamine reuptake with strongest affinity for the noradrenaline transporter) and escitalopram (an SSRI). Eight hundred and eleven adults with ICD–1033/DSM–IV34 unipolar major depression of at least moderate severity35 were recruited in eight European countries: Belgium, Croatia, Denmark, Germany, Italy, Poland, Slovenia and the UK. To minimise confounding by population stratification, recruitment was restricted to individuals of White European parentage. Personal or family history of bipolar disorder or schizophrenia constituted exclusion criteria. The study was approved by ethics boards in all participating centres. All participants provided a written consent after the procedures were fully explained.

Participants included 296 men and 514 women between 19 and 72 years old (Table 1). The average participant was in her second episode of moderately severe depression and scored 28.7 (s.d. = 6.7) on the Montgomery–Åsberg Depression Rating Scale (MADRS),36 21.7 (s.d. = 5.3) on the 17-item Hamilton Rating Scale for Depression (HRSD–17)37,38 and 28.0 (s.d. = 9.7) on the Beck Depression Inventory (BDI)39 at baseline.

View this table:
Table 1

Sample characteristics at baseline. The demographic and clinical characteristics are given by drug and 5-HTTLPR genotype

Participants with no contraindications were randomly allocated to receive flexible dosage of nortriptyline (50–150 mg daily) or escitalopram (10–30 mg daily) for 12 weeks. Individuals with relative or absolute contraindications (medical conditions incompatible with the use of the antidepressant and history of intolerance or non-response to one study medication) for one of the drugs were allocated non-randomly to the other antidepressant. Other psychotropic medication was not allowed with the exception of occasional use of hypnotics. Adherence was monitored by weekly self-report and plasma levels of antidepressants were measured at week 8. Adverse effects of medication were measured using the UKU Side-Effect Rating Scale40 and Self-Report Antidepressant Side-Effect Checklist.41

Outcome measurement

The response to antidepressant medication is a complex phenotype that includes changes in a number of symptoms occurring over a period of up to 12 weeks following the initiation of an antidepressant. In GENDEP, the responses to escitalopram and nortriptyline were assessed by weekly administration of three established measures of depression severity: the clinician-rated 10-item MADRS, the HRSD–17 and the self-report 21-item BDI were administered at week 0 (baseline) and then weekly for 12 weeks. All raters were trained and achieved high interrater reliability, with intraclass correlation of at least 0.9.42 A psychometric analysis found that the MADRS was the most internally consistent and informative of the three scales and that depressive symptoms could be described in more detail by three symptom dimensions derived by categorical item factor analysis: observed mood, cognitive symptoms and neurovegetative symptoms.42 In the pharmacogenetic analyses we used the MADRS score as a primary continuous outcome variable. Findings of interest were further explored using the three symptom dimensions.

Genotyping

Adequate DNA samples were available for 795 of the 811 participants. A sample of 5–8 ml of blood was collected in ethylenediamine tetra-acetic acid and frozen. DNA was extracted using a standard procedure.43

A two-stage method was used to type the 5HTTLPR and rs25531. In the first stage, short and long alleles were determined by polymerase chain reaction (PCR). In the second stage, the A/G alleles of rs25531 single nucleotide polymorphism within the 5HTTLPR repeats were identified using restriction fragment length polymorphism. The PCR product was digested with the Msp1 enzyme, which cuts the amplified PCR product at two or three places, depending on the rs25531 genotype. The combinations of products that determine the genotypes are shown in Appendix 1.

Five of the microsatellite markers listed on the UCSC bioinformatics website (http://genome.ucsc.edu/ accessed in September 2004) were selected for investigation and genotyped on 20 individuals, using PCR and electrophoresis, to check whether common polymorphisms at these loci are present. Only STin2 (SLC6A4, intron2) and STin4 (SLC6A4, intron4) proved to be polymorphic, and thus were carried forward for full analysis on the GENDEP sample. Fluorescent forward primers were used for amplification of these microsatellites. Genotyping of the PCR products was performed using the ABI3130 sequencer.

Single nucleotide polymorphisms were selected that tag DNA sequence variation within the SLC6A4 gene and its 1000 bp margins using the SNPTagger program (www.broad.mit.edu/mpg/tagger) run in Windows XP and HAPMAP data on the CEPH CEU population with European ancestry (CEPH NCBI Build 35/UCSC hg17/May 2004 coordinates).44 Criteria for single nucleotide polymorphism selection included a minimal allele frequency of 5% in the White population and a pairwise r2 = 0.8. Ten single nucleotide polymorphisms met these criteria (Table 2). According to SNPTagger, these 10 single nucleotide polymorphisms provided 92% coverage of the DNA sequence variation in the SLC6A4 gene. The single nucleotide polymorphisms were genotyped using the SNPlexTM method. This included oligonucleotide ligation assay/PCR technology for allelic discrimination and ligation product amplification. The assay products were run on the ABI3130 sequencer and data were exported and analysed using GENEMAPPER software version 3.7 for Windows XP (Applied Biosystems). Eleven samples were genotyped in blind duplicates and the agreement was 100% for all markers.

View this table:
Table 2

Linkage disequilibrium, allele frequency and Hardy–Weinberg equilibriuma

In addition to the SLC6A4 markers, a panel of 35 single nucleotide polymorphism markers selected as being informative for population stratification in a European population45 were genotyped in all individuals to establish the presence of and allow effective control for any population stratification.

The linkage disequilibrium measures (D′ and R2) were calculated using the Haploview program version 4 for Windows XP.46 Markers were included irrespective of genotype distribution as deviations from the Hardy–Weinberg equilibrium are expected in a cases-only sample.47 However, Hardy–Weinberg equilibrium was tested and results are given in Table 2.

Statistical analysis

To use all available weekly data on response to antidepressants and provide unbiased estimates in the presence of missing values, effects of genotype markers were tested using linear mixed models with individual random intercept and slope and fitted with maximum likelihood. All models included the fixed effects of time (linear and quadratic), baseline depression severity, history of using antidepressants, drug, age and gender, and random effects of the individual and the recruitment centre. To account for correlations between repeated observations, the intercept and slope of time were allowed to vary randomly between participants. The results of regression are presented as standardised regression coefficients (β) with 95% confidence intervals. Standardised regression coefficients denote the number of standard deviations on the outcome measure per unit of the predictor variable. As the standard deviation of the MADRS was 7, a β of 0.14 corresponds to a difference of one point on the MADRS.

As there was strong prior evidence for the influence of 5-HTTLPR on response to serotonergic antidepressants, the 5-HTTLPR was first tested as a single marker without correcting for multiple testing and the model was further developed using a stepwise approach. The effect of genotype and its interactions with drug, gender and age were tested as fixed explanatory variables in the mixed linear models. Competing nested models were compared using likelihood ratio tests. The previously reported recessive short allele model7,6 wascomparedwithafull genotype model. To test the hypothesis that a polymorphism in linkage disequilibrium with the 5-HTLPR modulates its effect, other markers within the SLC6A4 gene were then added in order of linkage disequilibrium with the 5-HTTLPR and retained if model fit improved significantly according to likelihood ratio tests (P<0.05).

Haplotypes for markers within the SLC6A4 gene that were in linkage disequilibrium (D′>0.7) were estimated 12 times using a Bayesian procedure in Phase, version 2.1 for Windows XP48,49 and used as multiply imputed data-sets in the mixed linear models.50 This method has been demonstrated to provide unbiased estimates and to have an adequate confidence interval coverage even when haplotypes are estimated based on cases only.50

The ancestry-informative markers were analysed using 10 000 burn-ins and 10 000 repetitions under an admixture model with correlated allele frequencies in STRUCTURE version 2.2 for Windows XP (http://pritch.bsd.uchicago.edu/software.html).51 This analysis found no clustering within the sample. Stratification was further estimated using the genomic control approach52 implemented in PLINK for Windows XP (http://pngu.mgh.harvard.edu/~purcell/plink).53 Based on the same set of ancestry-informative single nucleotide polymorphisms and all candidate gene markers, PLINK estimated a genomic inflation factor of 1, confirming absence of significant stratification in the GENDEP sample.

Power analysis

The sample size required for detecting an interaction between antidepressant and the 5-HTTLPR genotype was estimated based on data from published European studies that provided continuous outcome measures by genotype group.6,54 These studies showed that 5-HTTLPR influenced outcome with an effect size of 0.65–1.00 s.d. on a continuous outcome scale. Assuming a recessive model with minor allele frequency of 0.4, an effect size of 0.65 in participants treated with escitalopram and no effect of genotype on outcome in participants treated with nortriptyline, 270 participants per treatment arm were needed to detect a drug×genotype interaction at α = 0.05 with a power of 0.8. The GENDEP sample provided a power of more than 0.9.

Results

Genotyping

The 5-HTTLPR, intron 2 and intron 4 tandem repeats and 11 single nucleotide polymorphisms were typed in 795 (98%) of the 811 participants. Sample demographic and clinical characteristics at baseline were unrelated to the 5-HTTLPR genotype (all P>0.05, Table 1). There was also no significant relationship between the 5-HTTLPR genotype and the mode of allocation (random v. non-random, χ2 = 1.34, P = 0.246) or the drug allocated among non-randomised participants (χ2 = 0.0001, P = 0.994). The 16 participants with missing or inadequate DNA samples did not differ from the genotyped participants on depression severity, age or gender (all P>0.05). Table 2 shows the linkage disequilibrium among the 14 genotyped markers and minor allele frequencies. Frequencies of the 5-HTTLPR, s25531, rs2020933 and STin2 genotypes are given in Table 3.

View this table:
Table 3

Distribution of genotypes at the 5-HTTLPR, rs25531, rs2020933 and STin2 locia

The 5-HTTLPR

The 5-HTTLPR genotype moderated response to the serotonergic antidepressant in the expected direction. A significant three-way interaction between the 5-HTTLPR, drug and gender (P = 0.026; Table 4) indicated that this moderating effect is conditional on the type of antidepressant drug and on gender. Specifically, male carriers of the 5-HTTLPR long allele had a better response to escitalopram than male short-allele homozygotes (Fig. 1). The best-fitting model included a recessive effect of the 5-HTTLPR short allele, its two-way interactions with gender and drug, as well as the three-way interaction. Age did not significantly interact with genotype.

View this table:
Table 4

Results of the final best-fitting models for overall depression symptoms (Montgomery–Åsberg Depression Rating Scale (MADRS))a

Fig. 1

Moderation of antidepressant response by 5-HTTLPR genotype and gender: (a) men treated with escitalopram, (b) women treated with escitalopram, (c) men treated with nortriptyline, (d) women treated with nortriptyline. Depression severity measured by the Montgomery–Åsberg Depression Rating Scale (MADRS) in observed cases over the 12 weeks of treatment. Details of genotype are given in key; l, long allele; s, short allele. Note that values for weeks 9–12 are less reliable because of missing values resulting from drop-out and switching.

The interaction effects were qualified by significant main effects in participants treated with escitalopram and especially men treated with escitalopram. Among participants treated with escitalopram those homozygous for the short allele had significantly worse outcome than long-allele carriers (β = 0.154, 95% CI 0.003–0.305, P = 0.045). The effect of the 5-HTTLPR genotype was strong among males treated with escitalopram (β = 0.333, 95% CI 0.085–0.580, P = 0.008, Fig. 1a) but absent among females treated with escitalopram (β = 0.036, 95% CI –0.155 to 0.228, P = 0.709, Fig. 1b). The 5-HTTLPR genotype did not affect outcome among the participants treated with nortriptyline (β = –0.063, 95% CI – 0.231 to 0.106, P = 0.466, Fig. 1c and 1d). The 5-HTTLPR genotype explained 1.1% of variance in the repeated measurement of depression severity among men treated with escitalopram, but only 0.3% of variance in antidepressant response across the whole sample.

The additive or full genotype model provided no advantage over the recessive short-allele model (all likelihood ratio tests had P>0.1). The 5-HTTLPR genotype was unrelated to baseline severity, antidepressant dosage, self-reported adherence, plasma levels of either nortriptyline or escitalopram and adverse effects.

rs25531

The single nucleotide polymorphism rs25531 (A/G; minor allele frequency 0.075) was in linkage disequilibrium with the 5-HTTLPR (D′ = 0.82, R2 = 0.18).

Given the controversies about the moderation of 5-HTTLPR functionality by the rs25531,17,18 we tested whether the recoding of the 5-HTTLPR based on the rs25531 presumed functionality,17 or its inclusion as a separate marker improved the model fit. This was not the case: the recoding impaired model fit and its inclusion as a separate marker added no information to the prediction of antidepressant response (likelihood ratio test P>0.1). Therefore, the rs25531 was dropped from further analyses. The rs25531-recoded 5-HTTLPR genotype was also not related to baseline severity, antidepressant dosage, self-reported adherence, plasma levels of either nortriptyline or escitalopram and adverse effects.

rs2020933

The single nucleotide polymorphism rs2020933 (A/T; minor allele frequency 0.06) was in linkage disequilibrium with the 5-HTTLPR (D′ = 0.90, R2 = 0.04). The addition of rs2020933 (additive or dominant coding) to the model as a fixed first-order effect provided a significant improvement in the prediction of antidepressant response. Interaction between rs2020933 and 5-HTTLPR, drug, gender and age did not significantly improve the model. The final best-fitting model included both 5-HTTLPR with its interacting effects and the rs2020933 and led to a significant improvement in the prediction of antidepressant response compared with the model with no genotype data (likelihood ratio test χ 2(5) = 13.57, P = 0.0219, Table 4).

The intron 2 simple tandem repeat

The STin2 was in weak linkage disequilibrium with the 5-HTTLPR (D′ = 0.40, R2 = 0.07). The STin2 did not contribute to the prediction of treatment response (likelihood ratio test P>0.1).

Haplotype analysis

As there was significant linkage disequilibrium between the 5-HTTLPR and the rs2020933, the influence of estimated haplotypes comprising these two markers on antidepressant response was investigated. Although the haplotypes containing the 5-HTTLPR long allele were associated with better response to escitalopram, the inclusion of haplotypes provided no advantage over single marker analysis in predicting drug response.

Other markers within the SLC6A4 gene

Other markers within the SLC6A4 gene had low linkage disequilibrium with 5-HTTLPR and their inclusion as single markers or haplotypes conferred no improvement in the prediction model.

Effects on symptom dimensions

Our previous analyses have indicated that separable dimensions of depressive symptoms may differ in their response to antidepressants.32 We therefore investigated the genotype effect on change on the three symptoms dimensions.42 The 5-HTTLPR and rs2020933 significantly moderated response on the observed mood and cognitive symptom dimensions. The change on the neurovegetative symptom dimension was relatively independent of genotype (Table 5). The effect on the observed mood dimension was stronger than on the MADRS.

View this table:
Table 5

Results of the final model for specific dimensions of observed mood, cognitive and neurovegetative symptomsa

Discussion

The GENDEP data support the previously reported moderation of antidepressant response by the length polymorphism in the promoter region of the serotonin transporter gene. This effect is specific to the serotonergic mode of antidepressant action, appears to be concentrated in men and may be moderated by another polymorphism at the 5′ end of the serotonin transporter gene.

Gender effects

Within the GENDEP sample, the effect of the 5-HTTLPR on response to escitalopram was marked in men but absent in women. Among men treated with escitalopram, the 5-HTTLPR polymorphism is associated with a five-point difference on the MADRS between short-allele homozygotes and long-allele carriers. This finding is consistent with previous reports of gender-specific genetic influences on serotonergic function.55 It may reflect a biological interaction between the serotonergic system and ovarian hormones. Through the oestrogen alpha receptor, oestrogens stimulate the production of the 5-HT1A receptor that is involved in the regulation of serotonin release and is down-regulated in response to serotonin reuptake inhibitors.27 Oestrogens also increase the expression of the serotonin transporter.29 The impact of experimental manipulations of serotonergic function also depends on gender and hormonal status. Decrease in serotonergic function during tryptophan depletion has different effects on men and women30 and oestrogen supplementation may enhance the action of antidepressants in perimenopause.57,58 It is possible that the impact of the less functional short 5-HTTLPR allele is moderated by the oestrogen-induced stimulatory effect on serotonin transporter expression in hormonally active women. Although this accumulated evidence adds biological plausibility to the observed 5-HTTLPR×gender interaction, this finding requires replication in an independent sample.

Age effects

Although gender appears to be an important factor, age did not significantly interact with genotype in its effect on outcome. Joyce et al20 reported interaction between age and 5-HTTLPR with age dichotomised at 25. As the GENDEP sample is older, separate analysis of individuals below 25 years of age was not feasible. Therefore, the current negative results should not be interpreted as a refutation of the findings by Joyce and colleagues.20

The effect of rs2020933

The response to antidepressants in the GENDEP study was also associated with the single nucleotide polymorphism rs2020933, which is located in the first intron of the SLC6A4 gene. The major T allele was associated with a better response to both nortriptyline and escitalopram and there was no significant interaction with gender. The rs2020933 single nucleotide polymorphism has been reported to be functional in an allelic expression imbalance study.18 It is in significant linkage disequilibrium with the 5-HTTLPR, but because of a large difference in minor allele frequency, only a relatively small proportion of variance in the 5-HTTLPR is shared with the rs2020933. As this genomic region varies significantly between populations,11 the interplay between these two functional loci has the potential to explain discordant findings in samples of different ethnic origins.7,8 This remains to be tested in non-European samples. However, as the effect of rs2020933 was relatively weak and was not specific to the serotonergic mode of action, it is unlikely to fully explain the influence of the 5-HTTLPR on response to serotonergic antidepressants.

Drug effects on symptom dimensions

Depression is a heterogeneous group of disorders and the genetic mediation of response to antidepressants may differ between separable dimensions of depressive symptoms. For example, it has been demonstrated that change in sleep and appetite-related symptoms is relatively independent of change in core depressive symptoms such as mood and anhedonia.59 In a psychometric analysis of the GENDEP data, we have identified partially separable dimensions of observed mood, cognitive symptoms and neurovegetative symptoms that differ in their response to the two antidepressants.32,38 It has been previously reported that the moderation by the 5-HTTLPR of response to antidepressants is relatively specific to core mood symptoms and does not affect changes in sleep-related symptoms.60 The present analysis supports this distinction: the 5-HTTLPR and rs2020933 genotypes had the strongest influence on the observed mood dimension, comprising the core depressive symptoms and anxiety. These genotypes also moderated response on the cognitive symptom dimension. However, changes on the neurovegetative symptom dimension which comprises sleep, appetite and libido, were relatively independent of the 5-HTTLPR and rs2020933 genotypes. These findings add to the accumulating evidence indicating that mood and neurovegetative symptoms have distinct pathogenesis.

Effects of genotype on tolerability

It has been reported that tolerability rather than efficacy of antidepressants may be influenced by the 5-HTTLPR genotypes.61,62 We have therefore explored the relationship between the 5-HTTLPR genotype and adverse effects and related variables. We found that the 5-HTTLPR genotype was not related to escitalopram or nortriptyline dosage, plasma levels, adverse effects, drop-out or self-reported adherence in the GENDEP sample. These results remained negative after recoding according to the rs25531 marker. Thus we can conclude that, in the European population, the 5-HTTLPR genotype is directly related to antidepressant efficacy and there is no significant effect on tolerability.

Strengths and limitations

The GENDEP project is the largest study to date to compare a tricyclic antidepressant with an SSRI and is the second largest sample with pharmacogenetic data on antidepressant treatment outcome. It is however only powered to detect interacting effects of moderate size. The associations between polymorphisms in the SLC6A4 gene and response to antidepressants are of small to moderate size. Subsequently, the reported findings are detected with a modest level of certainty that satisfies the nominal statistical significance threshold but would not survive corrections for the number of polymorphisms that have been reported or suggested to moderate antidepressant response. Therefore, the results of the present study must be interpreted against the background of other published findings. The association of 5-HTTLPR with response to serotonergic antidepressants appears to be a robust and replicable finding in populations of European origin. However, the concentration of effect among males and the additional moderation by the rs2020933 polymorphism are novel findings that have to be tested in other large samples.

In conclusion, the available data suggest that variations in the 5′ end of the serotonin transporter gene have an effect on antidepressant efficacy, which is of modest effect size and depends on type of antidepressant and gender. Future studies should evaluate the population-specificity of this effect and the role of other polymorphisms including the rs2020933.

Appendix

5-HTTLPR and rs25531 genotyping

Restriction fragment length polymorphism analysis of PCR products carrying either the long (L, 16 repeat) or short (S, 14 repeat) alleles of the 5HTTLPR (L/S) and rs25531 (a/g).

Eight combinations were identified within the GENDEP sample: LaLa (324bp, 62bp and 33bp fragment), LaLg (324bp, 174bp, 150bp), LaSa (324bp, 281bp), SaLg (281bp, 174bp, 150bp), SaSa (281bp), LgLg (174bp, 150bp) and SaSg (281bp, 150bp, 131bp).

Acknowledgments

We would like to thank the following collaborators for their contribution: Helen Dean, Joanna Gray, Cerisse Gunasinghe, Amanda Elkin, Desmond Campbell, Richard J. Williamson, Julien Mendlewicz, Mara Barreto, Thomas Schulze, Christine Schmael, Susanne Höfels, Anna Schuhmacher, Ute Pfeiffer, Sandra Weber, Lisbeth Jorgensen, Anne Schinkel Stamp, Caterina Giovannini, Dejan Kozel, Moica Z. Dernovsek, Alenka Tancic, Jerneja Sveticic, Zrnka Kovacic, Pawe Kapelski, Maria Skibiñska, Aleksandra Rajewska, Monika Dmitrzak-Weglarz, Aleksandra Szczepankiewicz, and Elzbieta Cegielska. We would like to specially acknowledge the contribution of Jorge Perez, who was the principal investigator at Brescia, Italy, and who passed away in October 2007, and the important contribution made by Dr Andrej Marušič, the principal investigator at Ljubljana, Slovenia, who passed away in June 2008.

  • Received December 4, 2008.
  • Revision received January 12, 2009.
  • Accepted February 5, 2009.

References

View Abstract