Abstract
Despite strong evidence supporting the heritability of major depressive disorder (MDD), previous genome-wide studies were unable to identify risk loci among individuals of European descent. We used self-report data from 75,607 individuals reporting clinical diagnosis of depression and 231,747 individuals reporting no history of depression through 23andMe and carried out meta-analysis of these results with published MDD genome-wide association study results. We identified five independent variants from four regions associated with self-report of clinical diagnosis or treatment for depression. Loci with a P value <1.0 × 10−5 in the meta-analysis were further analyzed in a replication data set (45,773 cases and 106,354 controls) from 23andMe. A total of 17 independent SNPs from 15 regions reached genome-wide significance after joint analysis over all three data sets. Some of these loci were also implicated in genome-wide association studies of related psychiatric traits. These studies provide evidence for large-scale consumer genomic data as a powerful and efficient complement to data collected from traditional means of ascertainment for neuropsychiatric disease genomics.
This is a preview of subscription content, access via your institution
Access options
Subscribe to this journal
Receive 12 print issues and online access
£139.00 per year
only £11.58 per issue
Buy this article
- Purchase on SpringerLink
- Instant access to full article PDF
Prices may be subject to local taxes which are calculated during checkout
Similar content being viewed by others
References
Angst, F., Stassen, H.H., Clayton, P.J. & Angst, J. Mortality of patients with mood disorders: follow-up over 34-38 years. J. Affect. Disord. 68, 167–181 (2002).
Lopez, A.D., Mathers, C.D., Ezzati, M., Jamison, D.T. & Murray, C.J. Global and regional burden of disease and risk factors, 2001: systematic analysis of population health data. Lancet 367, 1747–1757 (2006).
Wittchen, H.U. et al. The size and burden of mental disorders and other disorders of the brain in Europe 2010. Eur. Neuropsychopharmacol. 21, 655–679 (2011).
Lichtenstein, P. et al. Recurrence risks for schizophrenia in a Swedish national cohort. Psychol. Med. 36, 1417–1425 (2006).
Sullivan, P.F., Kendler, K.S. & Neale, M.C. Schizophrenia as a complex trait: evidence from a meta-analysis of twin studies. Arch. Gen. Psychiatry 60, 1187–1192 (2003).
CONVERGE Consortium. Sparse whole-genome sequencing identifies two loci for major depressive disorder. Nature 523, 588–591 (2015).
Hawrylycz, M.J. et al. An anatomically comprehensive atlas of the adult human brain transcriptome. Nature 489, 391–399 (2012).
Le Meur, N. et al. MEF2C haploinsufficiency caused by either microdeletion of the 5q14.3 region or mutation is responsible for severe mental retardation with stereotypic movements, epilepsy and/or cerebral malformations. J. Med. Genet. 47, 22–29 (2010).
Paciorkowski, A.R. et al. MEF2C haploinsufficiency features consistent hyperkinesis, variable epilepsy, and has a role in dorsal and ventral neuronal developmental pathways. Neurogenetics 14, 99–111 (2013).
Barbosa, A.C. et al. MEF2C, a transcription factor that facilitates learning and memory by negative regulation of synapse numbers and function. Proc. Natl. Acad. Sci. USA 105, 9391–9396 (2008).
Wilkinson, M.B. et al. Imipramine treatment and resiliency exhibit similar chromatin regulation in the mouse nucleus accumbens in depression models. J. Neurosci. 29, 7820–7832 (2009).
Schizophrenia Working Group of the Psychiatric Genomics Consortium. Biological insights from 108 schizophrenia-associated genetic loci. Nature 511, 421–427 (2014).
Cross-Disorder Group of the Psychiatric Genomics Consortium. Identification of risk loci with shared effects on five major psychiatric disorders: a genome-wide analysis. Lancet 381, 1371–1379 (2013).
Major Depressive Disorder Working Group of the Psychiatric GWAS Consortium. A mega-analysis of genome-wide association studies for major depressive disorder. Mol. Psychiatry 18, 497–511 (2013).
Whitehead, A. Meta-Analysis of Controlled Clinical Trials (John Wiley & Sons, 2002).
Agoston, Z. et al. Meis2 is a Pax6 co-factor in neurogenesis and dopaminergic periglomerular fate specification in the adult olfactory bulb. Development 141, 28–38 (2014).
Agoston, Z., Li, N., Haslinger, A., Wizenmann, A. & Schulte, D. Genetic and physical interaction of Meis2, Pax3 and Pax7 during dorsal midbrain development. BMC Dev. Biol. 12, 10 (2012).
Pers, T.H. et al. Biological interpretation of genome-wide association studies using predicted gene functions. Nat. Commun. 6, 5890 (2015).
Purcell, S., Cherny, S.S. & Sham, P.C. Genetic Power Calculator: design of linkage and association genetic mapping studies of complex traits. Bioinformatics 19, 149–150 (2003).
Bulik-Sullivan, B.K. et al. LD Score regression distinguishes confounding from polygenicity in genome-wide association studies. Nat. Genet. 47, 291–295 (2015).
Okbay, A. et al. Genetic variants associated with subjective well-being, depressive symptoms, and neuroticism identified through genome-wide analyses. Nat. Genet. 48, 624–633 (2016).
Psychiatric GWAS Consortium Bipolar Disorder Working Group. Large-scale genome-wide association analysis of bipolar disorder identifies a new susceptibility locus near ODZ4. Nat. Genet. 43, 977–983 (2011).
Keller, M.B. et al. Results of the DSM-IV mood disorders field trial. Am. J. Psychiatry 152, 843–849 (1995).
Regier, D.A. et al. DSM-5 field trials in the United States and Canada, Part II: test–retest reliability of selected categorical diagnoses. Am. J. Psychiatry 170, 59–70 (2013).
Wisniewski, S.R. et al. Can phase III trial results of antidepressant medications be generalized to clinical practice? A STAR*D report. Am. J. Psychiatry 166, 599–607 (2009).
Castro, V.M. et al. Validation of electronic health record phenotyping of bipolar disorder cases and controls. Am. J. Psychiatry 172, 363–372 (2015).
Eriksson, N. et al. Web-based, participant-driven studies yield novel genetic associations for common traits. PLoS Genet. 6, e1000993 (2010).
Tung, J.Y. et al. Efficient replication of over 180 genetic associations with self-reported medical data. PLoS One 6, e23473 (2011).
1000 Genomes Project Consortium. An integrated map of genetic variation from 1,092 human genomes. Nature 491, 56–65 (2012).
Browning, S.R. & Browning, B.L. Rapid and accurate haplotype phasing and missing-data inference for whole-genome association studies by use of localized haplotype clustering. Am. J. Hum. Genet. 81, 1084–1097 (2007).
Fuchsberger, C., Abecasis, G.R. & Hinds, D.A. minimac2: faster genotype imputation. Bioinformatics 31, 782–784 (2015).
Durand, E.Y., Do, C.B., Mountain, J.L. & Macpherson, J.M. Ancestry Composition: a novel, efficient pipeline for ancestry deconvolution. Preprint at bioRxiv http://dx.doi.org/10.1101/010512 (2014).
Henn, B.M. et al. Cryptic distant relatives are common in both isolated and cosmopolitan genetic samples. PLoS One 7, e34267 (2012).
Lehoucq, R.B., Sorensen, D.C. & Yang, C. ARPACK Users' Guide: Solution of Large-Scale Eigenvalue Problems with Implicitly Restarted Arnoldi Methods (SIAM, 1998).
Tucker, G., Price, A.L. & Berger, B. Improving the power of GWAS and avoiding confounding from population stratification with PC-Select. Genetics 197, 1045–1049 (2014).
Fehrmann, R.S. et al. Gene expression analysis identifies global gene dosage sensitivity in cancer. Nat. Genet. 47, 115–125 (2015).
Acknowledgements
We would like to thank the research participants and employees of 23andMe for making this work possible. The authors thank the investigators and patient participants of the Psychiatric Genomic Consortium Major Depressive Disorder study for making the PGC MDD phase 1 results available for download. This study was supported in part by Pfizer, Inc. R.H.P. is supported in part by the National Institute of Mental Health and the National Human Genome Research Institute (P50 MH106933). We also thank the Social Science Genetics Association Consortium (SSGAC) for sharing results for subjective well-being, depressive symptoms, and neuroticism.
Author information
Authors and Affiliations
Contributions
A.R.W., C.L.H., and J.R.W. conceived the meta-analysis and statistical analysis. A.R.W., C.L.H., and R.H.P. oversaw data set analysis and primary data interpretation. C.L.H. designed and performed meta-analysis and further statistical analysis of the three data sets. X.C. provided statistical support and data visualization for the meta-analysis. M.W.N. provided DEPICT functional annotation and LD score regression analyses. R.H.P., A.R.W., and C.L.H. wrote the manuscript. A.R.W., R.H.P., C.L.H., D.A.H., S.A.P., and M.W.N. provided data interpretation and revised the manuscript. J.Y.T. and D.A.H. conceived and designed the 23andMe MDD GWAS. D.A.H. and C.T. performed GWAS for 23andMe data sets and statistical support.
Corresponding authors
Ethics declarations
Competing interests
R.H.P. has served on scientific advisory boards for or consulted to Genomind, Healthrageous, Perfect Health, Proteus Biomedical, Psybrain, and RID Ventures and receives royalties through Massachusetts General Hospital from Concordant Rater Systems (now Bracket). C.L.H., X.C., M.W.N., and S.A.P. are all employees and stockholders of Pfizer, Inc. C.T., D.A.H., and J.Y.T. are employees of and own stock or stock options in 23andMe, Inc. A.R.W. is a former employee and stockholder of Pfizer, Inc., and a current employee of the Perelman School of Medicine at the University of Pennsylvania Orphan Disease Center in partnership with the Loulou Foundation. J.R.W. is a former employee and stockholder of Pfizer, Inc., and a current employee and stockholder of Nestlé Health Science.
Integrated supplementary information
Supplementary Figure 1 Discovery-phase GWAS of 23andMe self-report ascertainment of major depression (75,607 cases and 231,747 controls).
Manhattan plot of 23andMe self-report ascertainment of major depression (75,607 cases and 231,747 controls). The threshold for genome-wide significance (P < 1 × 10−8) is indicated by the horizontal line. Red dots represent SNPs with P values smaller than the genome-wide significance threshold. Peaks with P values less than 5 × 10−8 (standard GWAS significance threshold) are labelled in black. (b) Quantile–quantile plot for the 23andMe MDD GWAS.
Supplementary Figure 2 The proportion of variance explained by each principal component in the subset of 23andMe participants with European ancestry.
(a,b) The proportion of variance explained by each principal component (a) and the proportion of each component’s variance that is explained by country of ancestry, for a set of individuals reporting four grandparents from a single country (b). The first five principal components are largely explained by geographical ancestry, whereas higher-order principal components are not.
Supplementary information
Supplementary Text and Figures
Supplementary Figures 1 and 2. (PDF 238 kb)
Supplementary Tables 1–11, 13 and 14
Supplementary Tables 1–11, 13 and 14. (XLSX 80 kb)
Supplementary Table 12
Summary statistics from the most significant 10,000 SNPs from the 23andMe discovery data set. (XLSX 1012 kb)
Rights and permissions
About this article
Cite this article
Hyde, C., Nagle, M., Tian, C. et al. Identification of 15 genetic loci associated with risk of major depression in individuals of European descent. Nat Genet 48, 1031–1036 (2016). https://doi.org/10.1038/ng.3623
Received:
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1038/ng.3623