Genome-wide association study of major depressive disorder: new results, meta-analysis, and lessons learned.

Major depressive disorder (MDD) is a common complex disorder with a partly genetic etiology. We conducted a genome-wide association study of the MDD2000+ sample (2431 cases, 3673 screened controls and >1 M imputed single-nucleotide polymorphisms (SNPs)). No SNPs achieved genome-wide significance either in the MDD2000+ study, or in meta-analysis with two other studies totaling 5763 cases and 6901 controls. These results imply that common variants of intermediate or large effect do not have main effects in the genetic architecture of MDD. Suggestive but notable results were (a) gene-based tests suggesting roles for adenylate cyclase 3 (ADCY3, 2p23.3) and galanin (GAL, 11q13.3); published functional evidence relates both of these to MDD and serotonergic signaling; (b) support for the bipolar disorder risk variant SNP rs1006737 in CACNA1C (P=0.020, odds ratio=1.10); and (c) lack of support for rs2251219, a SNP identified in a meta-analysis of affective disorder studies (P=0.51). We estimate that sample sizes 1.8- to 2.4-fold greater are needed for association studies of MDD compared with those for schizophrenia to detect variants that explain the same proportion of total variance in liability. Larger study cohorts characterized for genetic and environmental risk factors accumulated prospectively are likely to be needed to dissect more fully the etiology of MDD.

Genome-wide association for major depressive disorder: a possible role for the presynaptic protein piccolo.

Major depressive disorder (MDD) is a common complex trait with enormous public health significance. As part of the Genetic Association Information Network initiative of the US Foundation for the National Institutes of Health, we conducted a genome-wide association study of 435 291 single nucleotide polymorphisms (SNPs) genotyped in 1738 MDD cases and 1802 controls selected to be at low liability for MDD. Of the top 200, 11 signals localized to a 167 kb region overlapping the gene piccolo (PCLO, whose protein product localizes to the cytomatrix of the presynaptic active zone and is important in monoaminergic neurotransmission in the brain) with P-values of 7.7 x 10(-7) for rs2715148 and 1.2 x 10(-6) for rs2522833. We undertook replication of SNPs in this region in five independent samples (6079 MDD independent cases and 5893 controls) but no SNP exceeded the replication significance threshold when all replication samples were analyzed together. However, there was heterogeneity in the replication samples, and secondary analysis of the original sample with the sample of greatest similarity yielded P=6.4 x 10(-8) for the nonsynonymous SNP rs2522833 that gives rise to a serine to alanine substitution near a C2 calcium-binding domain of the PCLO protein. With the integrated replication effort, we present a specific hypothesis for further studies.

A common variant in the 3'UTR of the GRIK4 glutamate receptor gene affects transcript abundance and protects against bipolar disorder.

Underactivity of the glutamatergic system is an attractive model for the pathophysiology of several major mental illnesses. We previously described a chromosome abnormality disrupting the kainate class ionotropic glutamate receptor gene, GRIK4/KA1, in an individual with schizophrenia and learning disability (mental retardation). We also demonstrated in a case-control study that two physically separated haplotypes within this gene were significantly associated with increased risk of schizophrenia and decreased risk of bipolar disorder, respectively. The latter protective haplotype was located at the 3' end of the gene. We now report the identification from carriers of the protective haplotype of a deletion variant within the 3' untranslated region of the gene. The deletion allele also was found to be negatively associated with bipolar disorder in both initial (P = 0.00000019) and replication (P = 0.0107) case-control studies. Expression studies indicated that deletion-carrying mRNA transcripts were relatively more abundant. We postulate that this may be a direct consequence of the differences in the RNA secondary structures predicted for the insertion and deletion alleles. These data suggest a mechanism whereby the genetic protective effect is mediated through increased kainate receptor expression.

Whole-genome association study of bipolar disorder.

We performed a genome-wide association scan in 1461 patients with bipolar (BP) 1 disorder, 2008 controls drawn from the Systematic Treatment Enhancement Program for Bipolar Disorder and the University College London sample collections with successful genotyping for 372,193 single nucleotide polymorphisms (SNPs). Our strongest single SNP results are found in myosin5B (MYO5B; P=1.66 x 10(-7)) and tetraspanin-8 (TSPAN8; P=6.11 x 10(-7)). Haplotype analysis further supported single SNP results highlighting MYO5B, TSPAN8 and the epidermal growth factor receptor (MYO5B; P=2.04 x 10(-8), TSPAN8; P=7.57 x 10(-7) and EGFR; P=8.36 x 10(-8)). For replication, we genotyped 304 SNPs in family-based NIMH samples (n=409 trios) and University of Edinburgh case-control samples (n=365 cases, 351 controls) that did not provide independent replication after correction for multiple testing. A comparison of our strongest associations with the genome-wide scan of 1868 patients with BP disorder and 2938 controls who completed the scan as part of the Wellcome Trust Case-Control Consortium indicates concordant signals for SNPs within the voltage-dependent calcium channel, L-type, alpha 1C subunit (CACNA1C) gene. Given the heritability of BP disorder, the lack of agreement between studies emphasizes that susceptibility alleles are likely to be modest in effect size and require even larger samples for detection.

Evidence for association and epistasis at the DAOA/G30 and D-amino acid oxidase loci in an Irish schizophrenia sample.

The D-amino acid oxidase (DAO) signaling pathway has been implicated in schizophrenia pathogenesis. This may be mediated through modulation of NMDA function by DAO, which is in turn activated by DAO activator (DAOA, formerly G72). Chumakov et al. (2002); PNAS 99: 13675-13680, identifying the novel schizophrenia susceptibility gene DAOA/G30 and a number of independent studies have since reported evidence of association between the DAOA and DAO genes and schizophrenia. However, at least two studies have failed to replicate the epistatic interaction between these loci described in the original report and there have been differences in the associated alleles/haplotypes reported at each locus. In this study, we performed association and epistasis analyses of the DAOA/G30 and DAO loci in a sample of 373 cases with DSM-IV schizophrenia/schizoaffective disorder and 812 controls from the Republic of Ireland. Corrected for the number of tests performed, we found evidence for association between markers at both genes and schizophrenia: DAOA/G30 (P = 0.005, OR = 1.34 (1.09, 1.65)) and DAO (P = 0.003, OR = 1.43 (1.12, 1.84). The data suggest that evidence for association at DAO (marker rs2111902) is more consistent than previously realized, particularly in Caucasian schizophrenia populations. We identified evidence for epistatic interaction between the associated SNPs at DAOA and DAO genes in contributing to schizophrenia risk (OR = 9.3 (1.4, 60.5). Based on these data, more systematic investigation of genes involved in DAO signaling is required.