Investigation of possible epistatic interactions between GRIA2 and GRIA4 variants on clinical outcomes in patients with major depressive disorder

Abstract

Objectives

To investigate the effects of glutamate receptor, ionotropic, alpha-amino-3-hydroxy-5-methyl-4-isoxazole propionate (AMPA) 2 (GRIA2) rs4260586 and glutamate receptor, ionotropic, AMPA 4 (GRIA4) rs10736648 single nucleotide polymorphisms (SNPs) on response to antidepressants in Korean patients with major depressive disorder (MDD), and to ascertain whether epistatic interactions might exist between these SNPs.

Methods

In this retrospective analysis, patients were assessed at hospital admission and discharge using the Montgomery–Åsberg depression rating scale (MADRS). A multiple regression model was employed to investigate the effects of the two SNP variants on clinical/sociodemographic outcomes relating to MDD.

Results

Out of 145 Korean patients, the presence of both GRIA2 rs4260586 and GRIA4 rs10736648 polymorphisms had no significant association with MADRS improvement scores or other clinical/sociodemographic variables.

Conclusions

These data potentially suggest a lack of epistatic interaction between GRIA2 and GRIA4 variants, regarding clinical outcomes in patients with MDD. The study was limited by small sample size, use of different antidepressants and incomplete coverage of genes under investigation. Future research should include larger patient samples treated with different antidepressants, analysis of different SNPs and/or investigation of different gene–gene interactions within the glutamatergic system.

Keywords

Glutamate receptor ionotropic glutamate receptor epistasis major depression glutamatergic single nucleotide polymorphism

Introduction

Major depressive disorder (MDD) is a common mental disorder, with a lifetime US incidence of 12% in men and 20% in women.¹ Evidence suggests that MDD has a strong genetic aetiology.² Several genetic variants – including single nucleotide polymorphisms (SNPs) within genes encoding the serotonin transporter and several serotonin receptors – significantly modulate the response to antidepressants.³

Pharmacogenetic studies of psychotropic drug response that focus on one only SNP have shown inconsistent results or relatively small effects.⁴ Consequently, there is increasing emphasis on investigating SNPs within the genes encoding transporters and receptors that are unrelated to the serotonergic system. These include investigations of the glutamatergic system,⁵ or of epistatic interactions between genetic polymorphisms in candidate genes that may influence treatment outcomes in complex diseases, such as psychiatric disorders.^6,7

Research investigating possible associations between the glutamatergic system and antidepressant response has gained attention, with Maeng and Zarate hypothesizing that the efficacy of antidepressant drugs might be explained (at least in part) by direct or indirect actions of these drugs on the glutamatergic system.⁵ Consistent with this hypothesis, there is evidence that antidepressants can modulate glutamatergic transmission,⁸ and that variations within some glutamatergic genes can modulate antidepressant response.⁹

Epistasis is defined as a functional interaction between genes, encompassing various events including promoter activity control, epigenetic control, chromatin remodelling, and many other molecular reactions.¹⁰ These events can impact cell life-cycles and complex traits, and are orchestrated through genetically driven complex, yet flexible, activities.¹⁰ Thus, research into epistatic interactions between different candidate genes could be relevant in understanding the biological diversity that could influence antidepressant response.

Interactions between trace amine-associated receptor 6 and heat shock protein-70 gene polymorphisms have been found to affect the development of schizophrenia and mood disorders.^11,12 The present retrospective analysis of a published dataset¹³ investigated the relationship between genetic variants of glutamate receptor, ionotropic, alpha-amino-3-hydroxy-5-methyl-4-isoxazole propionate (AMPA) 2 (GRIA2) and glutamate receptor, ionotropic, AMPA 4 (GRIA4), in a sample of 145 Korean MDD patients treated with different antidepressants; such variants may influence antidepressant response. These two genes, which code for two different subunits of the glutamatergic AMPA-sensitive receptors, have been investigated individually;¹³ no association between the GRIA2 and GRIA4 polymorphisms investigated and clinical outcomes was found.

Patients and methods

Study population

The present retrospective study was conducted at Bucheon St Mary’s Hospital, The Catholic University of Korea College of Medicine, Bucheon, Kyounggi-Do, Republic of Korea, in October 2012. Details of the study population have been reported.¹³ Briefly, consecutive inpatients diagnosed with MDD according to the Diagnostic and Statistical Manual of Mental Disorders fourth edition (DSM-IV) criteria¹⁴ were included and underwent clinical interview using the Mini International Neuropsychiatric Interview (MINI).¹⁵ There was no restriction regarding the type of antidepressant treatment administered. Patients were excluded if they had: current severe (or unstable) medical or neurological conditions; current treatment with long-acting antipsychotic drugs; concomitant alcohol and substance abuse disorders. Patients were also excluded if they were not of Korean ethnicity.¹³ All patients admitted to the hospital – including those in the present study – are routinely assessed for severity of illness at admission and at discharge, by means of the Montgomery–Åsberg Depression Rating Scale (MADRS).¹⁶ Response was defined a priori as a ≥50% reduction in symptoms between admission and discharge; remission was defined as a MADRS score ≤7 at discharge, both as described previously.¹⁷ The following additional clinical and sociodemographic variables were recorded: total MADRS scores at admission and discharge; MDD subtypes; age; sex; age at MDD onset; family history of MDD; past history of suicide attempt; type of antidepressant and concomitant medication used.

The study protocol was approved by the Institutional Review Board of Bucheon St Mary’s Hospital, The Catholic University of Korea College of Medicine (approval number, HC10TISI0031). All patients included in the study provided written informed consent.

Outcome measures

The main outcome measure was the relationship between GRIA2 rs4260586 and GRIA4 rs10736648 SNPs and improvements in depressive symptoms, as measured by the change in total MADRS scores between admission and discharge. Secondary outcome measures included: presence of both GRIA2 rs4260586 and GRIA4 rs10736648 polymorphisms and total MADRS scores at admission and discharge; any response to treatment; remission; MDD subtypes; age; sex; age at MDD onset; family history of MDD; past history of suicide attempt; duration of admission; type of antidepressant used; concomitant anxiolytic therapy.

Genotyping of GRIA2 rs4260586 and GRIA4 rs10736648 SNPs

The protocol used to genotype the patient population has been described previously.¹³ To cover the highest possible proportion of the gene, the GRIA2 rs4260586 and GRIA4 rs10736648 SNPs were selected from SNPs within the original dataset.¹³ One SNP was selected for each gene that tagged the largest number of SNPs within the same gene. This approach ensured that ∼14% of GRIA2 variations and ∼1% of GRIA4 variations were covered.

Statistical analyses

Statistical analyses were performed using the STATISTICA software package, version 5 (StatSoft Italia, Vigonza, Padua, Italy) for Windows®. A multiple regression model was used to investigate the relationship between the two genotypes and clinical and sociodemographic variables. All P-values were 2-tailed and statistical significance was set at the P ≤ 0.004 level, after correction of multiple comparisons. These parameters gave sufficient power (0.80) to detect a medium–large (d = 0.29) effect size for patients carrying the TT rs4260586 genotype, compared with those carrying the AT genotype.¹⁸ Such an effect size corresponded with the possibility of detecting four-point differences in final MADRS scores.

Results

A total of 145 Korean inpatients with MDD, diagnosed according to DSM-IV criteria, were included. No significant correlation was observed between GRIA2 rs4260586 and GRIA4 rs10736648 and improvement of depressive symptoms, as measured by the change in total MADRS scores from admission to discharge (β = 0.17). In addition, no significant correlation was observed between rs4260586 within GRIA2 and rs10736648 within GRIA4, and secondary outcome measures (Table 1).

Table 1.

Influence of GRIA2 rs4260586 and GRIA4 rs10736648 single nucleotide polymorphisms on clinical and sociodemographic outcomes, in Korean patients with major depressive disorders (MDD).

Clinical or demographic characteristic	Patients n = 14	Influence of rs4260586	Influence of rs10736648	Influence of both rs4260586 and rs10736648
Clinical or demographic characteristic	Patients n = 14	β	β	β
Male sex	75 (52)	0.05	0.07	0.03
Age, years	41.37 ± 14.07	0.13	−0.02	0.05
Total MADRS score
Baseline	34.35 ± 8.95	−0.36	−0.43	0.48
Discharge	17.12 ± 9.88	−0.05	0.18	−0.23
Positive responders	78 (54)	−0.07	0.13	−0.13
Remitters	22 (15)	0.13	0.18	−0.23
Clinical subtypes^a
MDD without PF	105 (72)	0.14	0.28	−0.32
MDD with PF	11 (8)
Dysthymia	4 (3)
MDD NOS	5 (3)
Age at onset, years	38.08 ± 13.29	−0.04	−0.23	0.40
Family history of psychiatric disorders, yes	30 (21)	−0.07	0.11	−0.13
Suicide attempts, yes	36 (25)	0.01	−0.02	0.17
Duration of admission, days	32.31 ± 20.55	−0.37	0.02	0.18
Antipsychotic drug^b
Paroxetine	40 (28)	−0.02	0.20	−0.95
Venlafaxine	35 (24)
Fluoxetine	23 (16)
Mirtazapine	21 (14)
Other	3 (2)
Concomitant anxiolytics
Alprazolam	30 (21)	0.08	0.16	−0.24
Lorazepam	73 (50)
Clonazepam	3 (2)
Buspirone	4 (3)
None	35 (24)

Data presented as n (%) of patients or mean ± SD.

20 patients were excluded from this part of the analysis due to insufficient information; ^b23 patients were excluded from this part of the analysis due to insufficient information.

A multiple regression model was used to investigate the existence of possible interactions between the two genotypes, and clinical and sociodemographic variables.

No statistically significant differences were observed (P ≥ 0.004).

MADRS, Montgomery–Asberg Depression Rating Scale; PF, psychotic feature; NOS, not otherwise specified.

Discussion

In an attempt to understand whether potential epistatic interactions exist, the present study explored whether the presence of both the GRIA2 rs4260586 and GRIA4 rs10736648 SNPs was associated with improved MADRS scores and other clinical variables, in a sample of patients with MDD who were being treated with various antidepressants. There was no significant correlation between the presence of the two SNPs and improvement in MADRS scores, which potentially suggests that there is no epistatic interaction between the GRIA2 rs4260586 and GRIA4 rs10736648 polymorphisms, in terms of MDD. A possible explanation for this negative finding could be that the sample size was too small to detect subtle differences that are usually associated with single genes or gene–gene interactions, in complex disorders.¹⁹ Another possibility is that the results reflect a lack of epistatic interaction between the two SNPs investigated and clinical improvement, in terms of MADRS scores. In another study that focused on single SNPs within GRIA2 and GRIA4, no significant association with MADRS scores was observed.¹³ There is, however, some evidence to suggest that epistatic interactions can be observed, even when single SNPs considered on an individual basis do not lead to positive findings.^20,21 Such a possibility raises substantial computational and theoretical difficulties, when considering all possible gene–gene interactions and other forms of gene-expression control.²²

No statistically significant effects were observed regarding any of the secondary outcome measures in the presence of both GRIA2 rs4260586 and GRIA4 rs10736648 polymorphisms, further suggesting the possibility that no relevant epistatic interactions exist between the two SNPs under investigation. These data support results from a study that investigated the presence of these polymorphisms on an individual basis.¹³

Several limitations of the present study should be considered. First, the lack of associations observed could be simply due to the lack of statistical power that, in turn, could obscure small effects exerted by single SNPs. This issue is of particular concern if one considers that, even among pharmacogenetic studies with large sample sizes, results tend to remain conflicting.^23,24 A further limitation is the use of antidepressants with different mechanisms of action; this prevents conclusions from being drawn, with regard to the influence of these SNPs on specific drugs or classes of drugs. The decision to include patients who were being treated with different drugs could, however, have the advantage of being closer to ‘real-world’ clinical practice. The duration of hospitalization in the present study may be considered insufficient to ascertain a lack of response and remission; however, this time frame is consistent with common clinical practice.^25,26 Finally, this study was limited by the incomplete coverage of the genes under investigation, due to the tagging approach used.

The analysis of epistasis between different genes is particularly relevant to genetic-association studies, as it increases knowledge of the mechanism of action of the so-called reactome.^11,12 The analysis of both gene epistasis and the impact of single genotypes and haplotypes should be standard practice in future association studies.^11,12

The preliminary findings that are reported as a consequence of the present study suggest that it is likely that no epistatic interactions exist between GRIA2 and GRIA4, in terms of clinical outcomes, in patients with MDD treated with different antidepressants. Further research should include a larger sample size of patients with MDD and treated with different antidepressants, or include different SNPs and/or different gene–gene interactions within other genes that are known to be involved in the glutamatergic system.

Declaration of conflicting interest

The authors declare that there are no conflicts of interest.

Footnotes

Funding

This study was supported by a grant of the Korean Health Technology R&D Project, Ministry of Health & Welfare, Republic of Korea (A120004).

References

Kessler

Berglund

Demler

. The epidemiology of major depressive disorder: results from the National Comorbidiy Survey Replication (NCS-R). JAMA 2003; 289: 3095–3105.

Belmaker

Agam

. Major depressive disorder. N Engl J Med 2008; 358: 55–68.

Kato

. Possibility of a pharmacogenetic approach for prediction and personalized medication in major depressive disorder treatment. Nihon Shinkei Seishin Yakurigaku Zasshi 2010; 30: 83–92. [in Japanese, English abstract].

Malhotra

Murphy

GM Jr

Kennedy

. Pharmacogenetics of psychotropic drug response. Am J Psychiatry 2004; 161: 780–796.

Maeng

Zarate

CA Jr

. The role of glutamate in mood disorders: results from the ketamine in major depression study and the presumed cellular mechanism underlying its antidepressant effects. Curr Psychiatry Rep 2007; 9: 467–474.

Carlborg

Haley

. Epistasis: too often neglected in complex trait studies? Nat Rev Genet 2004; 5: 618–625.

Segurado

Bellgrove

Manconi

. Epistasis between neurochemical gene polymorphisms and risk for ADHD. Eur J Hum Genet 2011; 19: 577–582.

Kokras

Antoniou

Polissidis

. Antidepressants induce regionally discrete, sex-dependent changes in brain's glutamate content. Neurosci Lett 2009; 464: 98–102.

Fabbri C, Drago A and Serretti A. Early antidepressant efficacy modulation by glutamatergic gene variants in the STAR^*D. Eur Neuropsychopharmacol 2012; epub ahead of print 9 August 2012.

10.

Lane

Tsai

Lin

. Assessing gene-gene interactions in pharmacogenomics. Mol Diagn Ther 2012; 16: 15–27.

11.

Pae

Drago

Patkar

. Epistasis between a set of variations located in the TAAR6 and HSP-70 genes toward schizophrenia and response to antipsychotic treatment. Eur Neuropsychopharmacol 2009; 19: 806–811.

12.

Pae

Drago

Forlani

. Investigation of an epistastic effect between a set of TAAR6 and HSP-70 gene variations and major mood disorders. Am J Med Genet B Neuropsychiatr Genet 2010; 153B: 680–683.

13.

Chiesa

Crisafulli

Porcelli

. Influence of GRIA1, GRIA2 and GRIA4 polymorphisms on diagnosis and response to treatment in patients with major depressive disorder. Eur Arch Psychiatry Clin Neurosci 2012; 262: 305–311.

14.

American Psychiatric Association. Diagnostic and Statistical Manual of Mental Disorders, 4th edition (DSM-IV). Washington DC: American Psychiatric Press, 2004.

15.

Sheehan

Lecrubier

Sheehan

. The Mini-International Neuropsychiatric Interview (M.I.N.I.): the development and validation of a structured diagnostic psychiatric interview for DSM-IV and ICD-10. J Clin Psychiatry 1998; 59(suppl 20): 22–33. quiz 34–57.

16.

Montgomery

Asberg

. A new depression scale designed to be sensitive to change. Br J Psychiatry 1979; 134: 382–389.

17.

Riedel

Möller

Obermeier

. Response and remission criteria in major depression – a validation of current practice. J Psychiatr Res 2010; 44: 1063–1068.

18.

Cohen

. Statistical Power Analysis for the Behavioral Sciences, 2nd ed. Hillsdale, NJ: Lawrence Erlbaum Associates, 1988.

19.

Risch

Herrell

Lehner

. Interaction between the serotonin transporter gene (5-HTTLPR), stressful life events, and risk of depression: a meta-analysis. JAMA 2009; 301: 2462–2471.

20.

Arias

Beaumont

López

. Association of the peroxisome proliferator-activated receptor α gene L162V polymorphism with stage C heart failure. J Hypertens 2011; 29: 876–883.

21.

Arias

Freire

De La Roche

. Identification, inheritance, and variation of microsatellite markers in the black scallop Mimachlamys varia. Biochem Genet 2011; 49: 139–152.

22.

Serretti

Gibiino

Olgiati

. Pharmacogenetics of antidepressants and mood stabilizers. Handb Clin Neurol 2012; 106: 715–744.

23.

Abou Jamra

Becker

Georgi

. Genetic variation of the FAT gene at 4q35 is associated with bipolar affective disorder. Mol Psychiatry 2008; 13: 277–284.

24.

Shi

Badner

Gershon

. Allelic association of G72/G30 with schizophrenia and bipolar disorder: a comprehensive meta-analysis. Schizophr Res 2008; 98: 89–97.

25.

Zimmerman

Chelminski

McDermut

. Major depressive disorder and axis I diagnostic comorbidity. J Clin Psychiatry 2002; 63: 187–193.

26.

Zimmerman

Mattia

Posternak

. Are subjects in pharmacological treatment trials of depression representative of patients in routine clinical practice? Am J Psychiatry 2002; 159: 469–473.