Abelson Helper Integration Site-1 Gene Variants on Major Depressive Disorder and Bipolar Disorder

Article information

Psychiatry Investig. 2014;11(4):481-486
Publication date (electronic) : 2014 October 20
doi : https://doi.org/10.4306/pi.2014.11.4.481
1Institute of Psychiatry, Department of Biomedical and NeuroMotor Sciences, University of Bologna, Bologna, Italy.
2Department of Psychiatry, The Catholic University of Korea College of Medicine, Seoul, Republic of Korea.
3Department of Psychiatry and Behavioural Sciences, Duke University Medical Center, Durham, NC, USA.
4Department of Psychiatry, Korea University, College of Medicine, Seoul, Republic of Korea.
5Global Medical Education, New York, NY, USA.
Correspondence: Chi-Un Pae, MD. Department of Psychiatry, Bucheon St. Mary's Hospital, The Catholic University of Korea College of Medicine, 327 Sosa-ro, Wonmi-gu, Bucheon 420-717, Republic of Korea. Tel: +82-32-340-7067, Fax: +82-32-340-2255, pae@catholic.ac.kr
Received 2013 October 29; Revised 2014 January 18; Accepted 2014 January 19.

Abstract

Objective

The present study aimed to explore whether 4 single nucleotide polymorphisms (SNPs) within the AHI1 gene could be associated with major depressive disorder (MD) and bipolar disorder (BD), and whether they could predict clinical outcomes in mood disorders.

Methods

One hundred and eighty-four (184) patients with MD, 170 patients with BD and 170 healthy controls were genotyped for 4 AHI1 SNPs (rs11154801, rs7750586, rs9647635 and rs9321501). Baseline and final clinical measures for MD patients were assessed through the Hamilton Rating Scale for Depression (HAM-D). Allelic and genotypic frequencies in MD and BD subjects were compared with those of each disorder and healthy group using the χ2 statistics. Repeated measures ANOVA was used to test possible influences of SNPs on treatment efficacy.

Results

The rs9647635 A/A was more represented in subjects with BD as compared with MD and healthy subjects together. The rs9647635 A/A was also more presented in patients with MD than in healthy subjects. With regard to the allelic analysis, rs9647635 A allele was more represented in subjects with BD compared with healthy subjects, while it was not observed between patients with MD and healthy subjects.

Conclusion

Our findings provide potential evidence of an association between some variants of AHI1 and mood disorders susceptibility but not with clinical outcomes. However, we will need to do more adequately-powered and advanced association studies to draw any conclusion due to clear limitations.

INTRODUCTION

Mood disorders are severe, chronic psychiatric disorders that represent a major public health concern.1,2 Family, twin and adoption studies show evidence for a strong genetic component in mood disorders, with a relative contribution of genetic factors as high as 30-40% for major depressive disorder (MD)3 and 60-85% for bipolar disorder (BD).4

In the present study we focused on 4 single nucleotide polymorphisms (SNPs) within the Abelson helper integration site-1 gene (AHI1), which could be potentially involved in the etiology of MD and BD as well as in their treatment response.

The AHI1 gene locus, which is located on chromosome 6q23 (positive linkage with mood disorders in genome wide association study),5 encodes the protein Jouberin and is widely expressed in the brain.

The mouse orthologue of Jouberin, AHI1, binds to huntingtin-associated protein 1 (Hap1) to form a stable protein complex in the brain that is important for maintaining the level of tyrosine kinase receptor B (TrkB), which is critical for neuronal differentiation and brain development.6 Interestingly, the endogenous TrkB ligand in humans [i.e., the brain derived neurotrophic factor (BDNF)] is a survival factor for parvalbumine-positive interneurons, which have been thought to be involved in the physiopathology of mood disorders.2,7,8 Further, in the mouse brain AHI1 is abundant in the hypothalamus and amygdale,6,9 two important brain regions whose dysfunction can lead to affective disturbances.10

In a recent animal study,11 mice with neuronal AHI1 deficiency show reduced TrkB level in the brain and depressive phenotypes, which can be alleviated by antidepressant drugs or by overexpression of TrkB in the amygdala. It can be possibly assumed that overexpressed TrkB in the mouse amygdala could ameliorate the depressive phenotype of AHI1 mutant mice since the TrkB signaling plays a critical role in the depression associated with neuropsychiatric disorders.11,12,13 In addition, the loss of human AHI1 or mouse AHI1 in distinct types of cells seems to influence different cellular signaling pathways or function, which is also crucial pathway to be involved in the development and treatment response in mood disorders.14,15 Recently, animal research was taken to investigate the mechanism whereby alteration in AHI1 expression may be implicated in the pathogenesis of neuropsychiatric disorders. Lotan et al.16 conducted AHI1 heterozygous knockout (AHI1+/-) mice study to see the relationship between AHI1 gene and stress resilience, where they evidenced that under-expression of the AHI1 gene during neurodevelopment brings about relative resilience to various stressors during adulthood.16 Such resilience has been manifested by an anxiolytic-like behavioral phenotype, accompanied by a blunted response of the autonomic nervous system and the hypothalamus-pituitary-adrenal axis (HPA axis).16 That study proposes a potential link of AHI1 under-expression with a defect in the process of threat detection.16

Based on aforementioned, the present study aimed to investigate whether a set of SNPs within the AHI1 gene could be associated with MD and BD, and to explore whether such variants could be involved in clinical outcomes in patients with MD and BD.

METHODS

Subjects

The sample under investigation in the present study comprised 184 in-patients with MD and 170 in-patients with bipolar disorder (BD). All patients are native Koreans and were consecutively recruited at the Department of Psychiatry of The Catholic University of Korea. Patients were eligible for inclusion if they had a documented clinical diagnosis of MD and BD according to the DSM-IV criteria, as assessed by the Mini-International Neuropsychiatric Interview (MINI).17 The treatments were based on naturalistic treatment settings without any particular restrictions (i.e., type of pharmacological treatments and presence of concomitant co-morbidities). However, patients were excluded if they currently had severe or unstable medical and neurological conditions, were being treated with a long-acting antipsychotic, and had concomitant alcohol and substance abuse disorders. The choice not to adopt restrictive inclusion and exclusion criteria was motivated by the decision to include a sample of subjects that could be representative of the usual psychiatric inpatients in Korea. A further sample of 170 Korean psychiatrically healthy subjects, who underwent the same assessment as psychiatric patients to exclude possible psychiatric disorders and who were asked for the presence of any known psychiatric disorder in first and second degree relatives, and came from the same location as the psychiatric patients in the present study, was also included to compare genotype and allelic frequencies among three groups under investigation. All MD and BD patients admitted to the hospital were assessed for the severity of illness both at baseline and discharge by the administration of psychometric rating scales. In particular, MD severity was assessed by the administration of the Hamilton Depression Rating Scale-17 items (HAM-D).18 Raters were trained with the specific instruments with good inter-rater reliability (k>0.8). On the other hand, no psychometric evaluation was available for BD patients. Additionally, the following clinical and demographic variables were recorded for all the patients recruited: gender, age, age at onset, familiar history of psychiatric disorders, lifetime suicide attempts, duration of admission, medications at discharge and concomitant anxiolytics. The study protocol was approved by the institutional review board (approval number HC10TISI0031).

Statistical analyses

Statistical analyses were performed using "R cran" environment (http://cran.r-project.org/). The main outcome measures of the present study were: 1) differences between genotypic and allelic frequencies in patients with BD and MD compared with each group, 2) possible influence of the 4 SNPs under investigation on clinical improvement as measured by HAM-D in MD patients.

Differences in the allelic and genetic frequencies between healthy subjects and patients were calculated using the χ2 statistics. Repeated measures ANOVA was used to investigate the association among genotypes and the variation over time of HAM-D scores in MD patients. In the case of positive findings, the following clinical variables were added as covariates in order to investigate possible confounders: gender, age, age at onset, familiar history of psychiatric disorders, lifetime suicide attempts, duration of admission, medications at discharge and concomitant anxiolytics.

Haploview 3.2 was used to generate a linkage disequilibrium (LD) map and to test for Hardy-Weinberg equilibrium (HWE).19 Tests for associations using multi-marker haplotypes were performed using the statistics environment "R" package "haplo.score". Gender, age, age at onset, familiar history of psychiatric disorders, lifetime suicide attempts, duration of admission, medications at discharge and concomitant anxiolytics were added as covariates in the case of positive findings. Permutations (n=10.000) were performed to estimate the global significance of the results for all haplotype analysis. Only haplotypes with >1% prevalence were included in the analysis. All p-values were 2-tailed. In order to reduce the likelihood of false positive findings, statistical significance was set at the level of 0.01, as calculated by means of the false discovery rate (FDR) which allows for a correction of multiple testing without being as conservative as the Bonferroni correction.20

G Power (http://www.psycho.uni-duesseldorf.de/aap/projects/gpower/) was employed to perform the power analysis. With these parameters (p=0.01) we had a sufficient power (0.80) to detect an effect size of 0.2, that as an example, corresponded to an odds ratio of 1.9 between the 3 groups of patients and the group of controls and to detect small-medium effect sizes (d=0.3 for each diagnosis) for patients with BD and MD, carrying the CC genotype of rs11154801 as compared with those carrying the AC genotype.

SNP selection

SNP relative information was retrieved by the Database of SNPs. A total of four SNPs (rs11154801, rs7750586, rs9647635, rs9321501) for AHI1 were selected to be investigated. The main criteria for this selection concern the report of their genetic frequency, their role in translation, their potential functional significance as well as bibliographic references were strongly taken under consideration.

Genotyping

Genomic DNA was extracted from blood by standard methods and quantified. High-throughput genotyping using a pyrosequencer (Biotage AB, Sweden) was used for genotyping the four SNPs (rs11154801, rs7750586, rs9647635, rs9321501) of AHI1 under investigation. PCR primers (Bioneer, Daejeon, Korea) and sequencing primers (Bioneer) used for the pyrosequencing assay were designed by using the Pyrosequencing Assay Design Software v.1 (Biotage), and 1 primer of each primer set was biotinylated. Two independent investigators blind to the clinical status of samples checked all genotypes independently. Samples showing ambiguous alleles were discarded if they showed the same features on repeated genotyping. The final call rate was more than 99.6% for each SNP.

RESULTS

Socio-demographic features of the samples such as gender, age and further clinical and socio-demographical variables are reported in Table 1. For control subjects, only data on gender and age were collected. Patients and controls did not differ with regard to gender and age. All the considered SNPs were in HWE in the whole sample (Table 2). Strong linkage disequilibrium was observed between all SNPs, particularly between rs7750586 and rs9647635, rs11154801 and rs9647635, rs11154801 and rs7750586 (Figure 1).

Table 1

Clinical and demographic characteristics of the sample

Table 2

Single Nucleotide Polymorphisms considered in this study

Figure 1

Linkage disequilibrium (LD) plot.

Differences between Genotype and Allele Frequencies in MD and BD Patients and Healthy Controls

There were significant differences in the genotype frequencies of rs9647635 across different groups of subjects (Table 3). The rs9647635 A/A was more represented in subjects with BD as compared with patients with MD (χ2=20.345; p<0.001) and healthy subjects (χ2=34.537; p=0.008). In addition, the rs9647635 A/A was also more presented in patients with MD than in healthy subjects (χ2=9.723; p=0.008). These results were partially confirmed in the allelic analysis: The rs9647635 A allele was more represented in subjects with BD (Table 3).

Table 3

Genotype and allele frequency of the SNPs under investigation in the present study

However, there were no further significant differences between genotype and allele frequencies in patients with MD, BD and in healthy controls. Further, there were no differences in haplotype frequencies across samples.

AHI1 variants and clinical improvement in MD subjects

The haplotype analysis did not reveal any significant association with clinical improvement. No further genotype, allele and haplotype under investigation were significantly associated with improvement on HAM-D total scores.

DISCUSSION

The present study aimed to investigate whether 4 SNPs within the AHI1 gene could be associated with MD and BD, and whether the same variants could predict clinical outcomes in MD patients treated with antidepressants.

To the best of our knowledge, this is the first study investigating a relationship between AHI1 and mood disorders. Our results suggest a possible role of rs9647635 on BD and MD susceptibility: a significant higher proportion of subjects with BD carried the A/A genotype and the A allele when compared with other groups. Although such significant genotype difference of rs9647635 was also seen between patients with MD and healthy subjects, other SNP markers failed completely to do so. Despite scarce information about AHI1 gene for mood disorders, our exploratory results suggest that rs 9647635 should be possibly involved in the development of MD and BD, respectively, although such variants were not associated with the treatment response in patients with MD.

Considering that AHI1's extensive developmental role, from essential cellular signaling organelles, such as the primary cilium, to neuronal networks and complex organ systems, it should not be excessive if we think that Ahi1 gene should be a promising candidate gene for neuropsychiatric disorders including mood disorders.16 In fact, a number of research groups have shown significant association of the AHI1 gene with schizophrenia21 and autism spectrum disorders,22 which are replicated in a large-scale association studies.23,24 Furthermore, according to the recent AHI1 expression study in patients with BD and schizophrenia,25 the patients with early age of onset showed higher AHI1 expression than controls and later onset patients with schizophrenia. In such study, although no difference in brain expression of AHI1 in schizophrenia or BD patients compared to controls was found, there was a genotypic difference in AHI1 expression for SNP rs9321501. Genotypes that included the under-transmitted C allele (CC/AC) showed lower expression than the homozygous AA genotype. It might be possibly interpreted that the higher expression in AHI1 mRNA in early-onset and late-onset schizophrenia patients may be a putative endophenotypes of schi-zophrenia, such as neurocognitive performance.

According to the recent animal model,16 the under-expression of the Ahi1 gene during neurodevelopment was found to cause a relative resilience to various stressors during adulthood. Such resilience was exerted by an anxiolytic-like behavioral phenotype, accompanied by a blunted response of the autonomic nervous system and the HPA axis. Another animal model11 also proposed that neuronal AHI1 deficiency may cause a reduced TrkB level in the brain and should bring about depressive phenotypes, which can be alleviated by antidepressant drugs or by overexpression of TrkB in the amygdala. AHI1 deficiency also prompts the degradation of endocytic-TrkB and reduces TrkB signaling in neuronal cells, supporting that increased degradation of TrkB can induce depression and that such impaired pathway may serve as one of putative therapeutic targets for depression.11 Hence, our preliminary findings may add valuable human genetic data for the investigation of the role of AHI1 gene in the development and treatment response for mood disorders in combination with those form previous innovative animal model studies.

However, there are many limitations of the present work. Unfortunately, it was not possible to investigate the role of rs9647635 in response to mood stabilizers due to the lack of psychometric assessment in bipolar patients. Candidate gene studies, such as the present one, are associated with a high likelihood of false positive findings.26 Hence, further replications in independent samples are needed to confirm our results. A further concern is related to the use of several drugs with different mechanisms of action for each cohort of patients that do not allow one to draw definitive conclusions with regard to the influence of the SNPs under investigation upon specific drugs or classes of drugs. However, our decision to include patients treated with different drugs could have the advantage of being closer to "real world" clinical practice. Additionally, our categorization of psychiatric disorder was based on current DSM-IV criteria, as assessed by the MINI. However, with this instrument there is neither the ability to rule out the possibility that patients with MD could switch in the future to a diagnosis of BD, nor to exclude a possible switch to schizoaffective disorder. Such issues are important for genetic studies, as shown by recent findings suggesting that BD and schizophrenia could share many important risk genes.27 Furthermore, we asked healthy controls to report only known psychiatric disorders among first and second-degree relatives, thus limiting the possibility to detect whether sub-threshold or untreated psychiatric disorders among family members of healthy control subjects could exist. A further possible limitation of the present study could be imputed to the incomplete coverage of genes under investigation. Also, we did not control for population heterogeneity, though all subjects belonged to Korean ethnicity, which is considered to be genetically homogenous.28 Finally the small sample size does not allow our data to have generalizability and thus we have to in-crease the sample power in future studies.

In conclusion, our study preliminary suggests a possible role of AHI1 variants on mood disorders susceptibility, with rs9647635 being the most promising candidate polymorphism. However, further research is needed to confirm and extend our results, in particular future studies will need to include a higher number of SNPs covering larger portions of genes and more homogeneous groups of drugs.

Acknowledgments

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

References

1. Alda M. Genetic factors and treatment of mood disorders. Bipolar Disord 2001;3:318–324. 11843781.
2. Pae CU, Chiesa A, Porcelli S, Han C, Patkar AA, Lee SJ, et al. Influence of BDNF variants on diagnosis and response to treatment in patients with major depression, bipolar disorder and schizophrenia. Neuropsychobiology 2012;65:1–11. 22094229.
3. Kendler KS, Gatz M, Gardner CO, Pedersen NL. A Swedish national twin study of lifetime major depression. Am J Psychiatry 2006;163:109–114. 16390897.
4. Smoller JW, Finn CT. Family, twin, and adoption studies of bipolar disorder. Am J Med Genet C Semin Med Genet 2003;123C:48–58. 14601036.
5. Zubenko GS, Maher BS, Hughes HB 3rd, Zubenko WN, Scott Stiffler J, Marazita ML. Genome-wide linkage survey for genetic loci that affect the risk of suicide attempts in families with recurrent, early-onset, major depression. Am J Med Genet B Neuropsychiatr Genet 2004;129B:47–54. 15274040.
6. Sheng G, Xu X, Lin YF, Wang CE, Rong J, Cheng D, et al. Huntingtin-associated protein 1 interacts with Ahi1 to regulate cerebellar and brainstem development in mice. J Clin Invest 2008;118:2785–2795. 18636121.
7. Wang AY, Lohmann KM, Yang CK, Zimmerman EI, Pantazopoulos H, Herring N, et al. Bipolar disorder type 1 and schizophrenia are accompanied by decreased density of parvalbumin- and somatostatin-positive interneurons in the parahippocampal region. Acta Neuropathol 2011;122:615–626. 21968533.
8. Lisman JE, Coyle JT, Green RW, Javitt DC, Benes FM, Heckers S, et al. Circuit-based framework for understanding neurotransmitter and risk gene interactions in schizophrenia. Trends Neurosci 2008;31:234–242. 18395805.
9. Doering JE, Kane K, Hsiao YC, Yao C, Shi B, Slowik AD, et al. Species differences in the expression of Ahi1, a protein implicated in the neurodevelopmental disorder Joubert syndrome, with preferential accumulation to stigmoid bodies. J Comp Neurol 2008;511:238–256. 18785627.
10. David DJ, Samuels BA, Rainer Q, Wang JW, Marsteller D, Mendez I, et al. Neurogenesis-dependent and -independent effects of fluoxetine in an animal model of anxiety/depression. Neuron 2009;62:479–493. 19477151.
11. Xu X, Yang H, Lin YF, Li X, Cape A, Ressler KJ, et al. Neuronal Abelson helper integration site-1 (Ahi1) deficiency in mice alters TrkB signaling with a depressive phenotype. Proc Natl Acad Sci U S A 2010;107:19126–19131. 20956301.
12. Martinowich K, Manji H, Lu B. New insights into BDNF function in depression and anxiety. Nat Neurosci 2007;10:1089–1093. 17726474.
13. Angelucci F, Brene S, Mathe AA. BDNF in schizophrenia, depression and corresponding animal models. Mol Psychiatry 2005;10:345–352. 15655562.
14. Sperlagh B, Csolle C, Ando RD, Goloncser F, Kittel A, Baranyi M. The role of purinergic signaling in depressive disorders. Neuropsychopharmacol Hung 2012;14:231–238. 23269209.
15. Chiu CT, Wang Z, Hunsberger JG, Chuang DM. Therapeutic potential of mood stabilizers lithium and valproic acid: beyond bipolar disorder. Pharmacol Rev 2013;65:105–142. 23300133.
16. Lotan A, Lifschytz T, Slonimsky A, Broner EC, Greenbaum L, Abedat S, et al. Neural mechanisms underlying stress resilience in Ahi1 knockout mice: relevance to neuropsychiatric disorders. Mol Psychiatry 2014;19:243–252. 24042478.
17. Sheehan DV, Lecrubier Y, Sheehan KH, Amorim P, Janavs J, Weiller E, et al. 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. 9881538.
18. Montgomery SA, Asberg M. A new depression scale designed to be sensitive to change. Br J Psychiatry 1979;134:382–389. 444788.
19. Barrett JC, Fry B, Maller J, Daly MJ. Haploview: analysis and visualization of LD and haplotype maps. Bioinformatics 2005;21:263–265. 15297300.
20. Benjamini Y, Drai D, Elmer G, Kafkafi N, Golani I. Controlling the false discovery rate in behavior genetics research. Behav Brain Res 2001;125:279–284. 11682119.
21. Amann-Zalcenstein D, Avidan N, Kanyas K, Ebstein RP, Kohn Y, Hamdan A, et al. AHI1, a pivotal neurodevelopmental gene, and C6orf217 are associated with susceptibility to schizophrenia. Eur J Hum Genet 2006;14:1111–1119. 16773125.
22. Alvarez Retuerto AI, Cantor RM, Gleeson JG, Ustaszewska A, Schackwitz WS, Pennacchio LA, et al. Association of common variants in the Joubert syndrome gene (AHI1) with autism. Hum Mol Genet 2008;17:3887–3896. 18782849.
23. Ingason A, Giegling I, Cichon S, Hansen T, Rasmussen HB, Nielsen J, et al. A large replication study and meta-analysis in European samples provides further support for association of AHI1 markers with schizophrenia. Hum Mol Genet 2010;19:1379–1386. 20071346.
24. Ingason A, Sigmundsson T, Steinberg S, Sigurdsson E, Haraldsson M, Magnusdottir BB, et al. Support for involvement of the AHI1 locus in schizophrenia. Eur J Hum Genet 2007;15:988–991. 17473831.
25. Slonimsky A, Levy I, Kohn Y, Rigbi A, Ben-Asher E, Lancet D, et al. Lymphoblast and brain expression of AHI1 and the novel primate-specific gene, C6orf217, in schizophrenia and bipolar disorder. Schizophr Res 2010;120:159–166. 20452750.
26. Sullivan PF. Spurious genetic associations. Biol Psychiatry 2007;61:1121–1126. 17346679.
27. Lichtenstein P, Yip BH, Bjork C, Pawitan Y, Cannon TD, Sullivan PF, et al. Common genetic determinants of schizophrenia and bipolar disorder in Swedish families: a population-based study. Lancet 2009;373:234–239. 19150704.
28. Cavalli-Sforza L. The History and Geography of Human Genes New Jersey: Princeton University Press; 1994.
29. Chiesa A, Marsano A, Han C, Lee SJ, Patkar AA, Pae CU, et al. Epistatic Interactions between CREB and CREM Variants in Affective Disorder. Psychiatry Investig 2014;11:200–203.

Article information Continued

Funded by : Ministry of Health and Welfare
Award ID : HI12C0003

Figure 1

Linkage disequilibrium (LD) plot.

Table 1

Clinical and demographic characteristics of the sample

Table 1

Data represent number (%) or mean±standard deviation. BD: bipolar disorder, MD: major depressive disorder

Table 3

Genotype and allele frequency of the SNPs under investigation in the present study

Table 3

*statistical significance. BD: bipolar disorder, MD: major depressive disorder, SNPs: single nucleotide polymorphisms