Mapping genomic loci implicates genes and synaptic biology in schizophrenia.

Schizophrenia has a heritability of 60-80%1, much of which is attributable to common risk alleles. Here, in a two-stage genome-wide association study of up to 76,755 individuals with schizophrenia and 243,649 control individuals, we report common variant associations at 287 distinct genomic loci. Associations were concentrated in genes that are expressed in excitatory and inhibitory neurons of the central nervous system, but not in other tissues or cell types. Using fine-mapping and functional genomic data, we identify 120 genes (106 protein-coding) that are likely to underpin associations at some of these loci, including 16 genes with credible causal non-synonymous or untranslated region variation. We also implicate fundamental processes related to neuronal function, including synaptic organization, differentiation and transmission. Fine-mapped candidates were enriched for genes associated with rare disruptive coding variants in people with schizophrenia, including the glutamate receptor subunit GRIN2A and transcription factor SP4, and were also enriched for genes implicated by such variants in neurodevelopmental disorders. We identify biological processes relevant to schizophrenia pathophysiology; show convergence of common and rare variant associations in schizophrenia and neurodevelopmental disorders; and provide a resource of prioritized genes and variants to advance mechanistic studies.

Exploring the Genetic Risk of Childhood Daytime Urinary Incontinence: A Genome-Wide Association Study.

PURPOSE: Childhood incontinence is stigmatized and underprioritized, and a basic understanding of its pathogenesis is missing. Our goal was to identify risk-conferring genetic variants in daytime urinary incontinence (DUI). MATERIALS AND METHODS: We conducted a genome-wide association study in the Danish iPSYCH2015 cohort. Cases (3024) were identified through DUI diagnosis codes and redeemed prescriptions for DUI medication in patients aged 5 to 20 years. Controls (30,240), selected from the same sample, were matched to cases on sex and psychiatric diagnoses, if any, and down-sampled to a 1:10 case:control ratio. Replication was performed in the Icelandic deCODE cohort (5475 cases/287,773 controls). Single-nucleotide polymorphism heritability was calculated using the genome-based restricted maximum likelihood method. Cross-trait genetic correlation was estimated using linkage disequilibrium score regression. Polygenic risk scores generated with LDpred2-auto and BOLT-LMM were assessed for association. RESULTS: Variants on chromosome 6 (rs12210989, odds ratio [OR] 1.24, 95% CI 1.17-1.32, P = 3.21 × 10-12) and 20 (rs4809801, OR 1.18, 95% CI 1.11-1.25, P = 3.66 × 10-8) reached genome-wide significance and implicated the PRDM13 and RIPOR3 genes. Chromosome 6 findings were replicated (P = .024, OR 1.09, 95% CI 1.01-1.16). Liability scale heritability ranged from 10.20% (95% CI 6.40%-14.00%) to 15.30% (95% CI 9.66%-20.94%). DUI and nocturnal enuresis showed positive genetic correlation (rg = 1.28 ± 0.38, P = .0007). DUI was associated with attention-deficit/hyperactivity disorder (OR 1.098, 95% CI 1.046-1.152, P < .0001) and BMI (OR 1.129, 95% CI 1.081-1.178, P < .0001) polygenic risk. CONCLUSIONS: Common genetic variants contribute to the risk of childhood DUI, and genes important in neuronal development and detrusor smooth muscle activity were implicated. These findings may help guide identification of new treatment targets.

Polygenic liabilities and treatment trajectories in early-onset depression: a Danish register-based study.

BACKGROUND: The clinical course of major depressive disorder (MDD) is heterogeneous, and early-onset MDD often has a more severe and complex clinical course. Our goal was to determine whether polygenic scores (PGSs) for psychiatric disorders are associated with treatment trajectories in early-onset MDD treated in secondary care. METHODS: Data were drawn from the iPSYCH2015 sample, which includes all individuals born in Denmark between 1981 and 2008 who were treated in secondary care for depression between 1995 and 2015. We selected unrelated individuals of European ancestry with an MDD diagnosis between ages 10-25 (N = 10577). Seven-year trajectories of hospital contacts for depression were modeled using Latent Class Growth Analysis. Associations between PGS for MDD, bipolar disorder, schizophrenia, ADHD, and anorexia and trajectories of MDD contacts were modeled using multinomial logistic regressions. RESULTS: We identified four trajectory patterns: brief contact (65%), prolonged initial contact (20%), later re-entry (8%), and persistent contact (7%). Relative to the brief contact trajectory, higher PGS for ADHD was associated with a decreased odds of membership in the prolonged initial contact (odds ratio = 1.06, 95% confidence interval = 1.01-1.11) and persistent contact (1.12, 1.03-1.21) trajectories, while PGS-AN was associated with increased odds of membership in the persistent contact trajectory (1.12, 1.03-1.21). CONCLUSIONS: We found significant associations between polygenic liabilities for psychiatric disorders and treatment trajectories in patients with secondary-treated early-onset MDD. These findings help elucidate the relationship between a patient's genetics and their clinical course; however, the effect sizes are small and therefore unlikely to have predictive value in clinical settings.

Polygenic liability, stressful life events and risk for secondary-treated depression in early life: a nationwide register-based case-cohort study.

BACKGROUND: In this study, we examined the relationship between polygenic liability for depression and number of stressful life events (SLEs) as risk factors for early-onset depression treated in inpatient, outpatient or emergency room settings at psychiatric hospitals in Denmark. METHODS: Data were drawn from the iPSYCH2012 case-cohort sample, a population-based sample of individuals born in Denmark between 1981 and 2005. The sample included 18 532 individuals who were diagnosed with depression by a psychiatrist by age 31 years, and a comparison group of 20 184 individuals. Information on SLEs was obtained from nationwide registers and operationalized as a time-varying count variable. Hazard ratios and cumulative incidence rates were estimated using Cox regressions. RESULTS: Risk for depression increased by 35% with each standard deviation increase in polygenic liability (p < 0.0001), and 36% (p < 0.0001) with each additional SLE. There was a small interaction between polygenic liability and SLEs (ß = -0.04, p = 0.0009). The probability of being diagnosed with depression in a hospital-based setting between ages 15 and 31 years ranged from 1.5% among males in the lowest quartile of polygenic liability with 0 events by age 15, to 18.8% among females in the highest quartile of polygenic liability with 4+ events by age 15. CONCLUSIONS: These findings suggest that although there is minimal interaction between polygenic liability and SLEs as risk factors for hospital-treated depression, combining information on these two important risk factors could potentially be useful for identifying high-risk individuals.

Sequencing and de novo assembly of 150 genomes from Denmark as a population reference.

Hundreds of thousands of human genomes are now being sequenced to characterize genetic variation and use this information to augment association mapping studies of complex disorders and other phenotypic traits. Genetic variation is identified mainly by mapping short reads to the reference genome or by performing local assembly. However, these approaches are biased against discovery of structural variants and variation in the more complex parts of the genome. Hence, large-scale de novo assembly is needed. Here we show that it is possible to construct excellent de novo assemblies from high-coverage sequencing with mate-pair libraries extending up to 20 kilobases. We report de novo assemblies of 150 individuals (50 trios) from the GenomeDenmark project. The quality of these assemblies is similar to those obtained using the more expensive long-read technology. We use the assemblies to identify a rich set of structural variants including many novel insertions and demonstrate how this variant catalogue enables further deciphering of known association mapping signals. We leverage the assemblies to provide 100 completely resolved major histocompatibility complex haplotypes and to resolve major parts of the Y chromosome. Our study provides a regional reference genome that we expect will improve the power of future association mapping studies and hence pave the way for precision medicine initiatives, which now are being launched in many countries including Denmark.

A polygenic resilience score moderates the genetic risk for schizophrenia.

Based on the discovery by the Resilience Project (Chen R. et al. Nat Biotechnol 34:531-538, 2016) of rare variants that confer resistance to Mendelian disease, and protective alleles for some complex diseases, we posited the existence of genetic variants that promote resilience to highly heritable polygenic disorders1,0 such as schizophrenia. Resilience has been traditionally viewed as a psychological construct, although our use of the term resilience refers to a different construct that directly relates to the Resilience Project, namely: heritable variation that promotes resistance to disease by reducing the penetrance of risk loci, wherein resilience and risk loci operate orthogonal to one another. In this study, we established a procedure to identify unaffected individuals with relatively high polygenic risk for schizophrenia, and contrasted them with risk-matched schizophrenia cases to generate the first known "polygenic resilience score" that represents the additive contributions to SZ resistance by variants that are distinct from risk loci. The resilience score was derived from data compiled by the Psychiatric Genomics Consortium, and replicated in three independent samples. This work establishes a generalizable framework for finding resilience variants for any complex, heritable disorder.

Atlantic herring (Clupea harengus) population structure in the Northeast Atlantic Ocean.

The Atlantic herring Clupea harengus L has a vast geographical distribution and a complex population structure with a few very large migratory units and many small local populations. Each population has its own spawning ground and/or time, thereby maintaining their genetic integrity. Several herring populations migrate between common feeding grounds and over-wintering areas resulting in frequent mixing of populations. Thus, many herring fisheries are based on mixed populations of different demographic status. In order to avoid over-exploitation of weak populations and to conserve biodiversity, understanding the population structure and population mixing is important for maintaining biologically sustainable herring fisheries. The aim of this study was to investigate the genetic population structure of herring in the Faroese and surrounding waters, and to develop genetic markers for distinguishing between four herring management units (often called stocks), namely the Norwegian spring-spawning herring (NSSH), Icelandic summer-spawning herring (ISSH), North Sea autumn-spawning herring (NSAH), and Faroese autumn-spawning herring (FASH). Herring from the four stocks were sequenced at low coverage, and single nucleotide polymorphisms (SNPs) were called and used for population structure analysis and individual assignment. An ancestry-informative SNP panel with 118 SNPs was developed and tested on 240 individuals. The results showed that all four stocks appeared to be genetically differentiated populations, but at lower levels of differentiation between FASH and ISSH than the other two populations. Overall assignment rate with the SNP panel was 80.7%, and agreement between the genetic and traditional visual assignment was 75.5%. The NSAH and NSSH samples had the highest assignment rate (100% and 98.3%, respectively) and highest agreement between traditional and genetic assignment methods (96.6% and 94.9%, respectively). The FASH and ISSH samples had substantially lower assignment rates (72.9% and 51.7%, respectively) and agreement between traditional and genetic methods (39.5% and 48.4%, respectively).

Whole-exome sequencing implicates DGKH as a risk gene for panic disorder in the Faroese population.

The demographic history of the isolated population of the Faroe Islands may have induced enrichment of variants rarely seen in outbred European populations, including enrichment of risk variants for panic disorder (PD). PD is a common mental disorder, characterized by recurring and unprovoked panic attacks, and genetic factors have been estimated to explain around 40% of the risk. In this study the potential enrichment of PD risk variants was explored based on whole-exome sequencing of 54 patients with PD and 211 control individuals from the Faroese population. No genome-wide significant associations were found, however several single variants and genes showed strong association with PD, where DGKH was found to be the strongest PD associated gene. Interestingly DGKH has previously demonstrated genome-wide significant association with bipolar disorder as well as evidence of association to other mental disorders. Additionally, we found an enrichment of PD risk variants in the Faroese population; variants with otherwise low frequency in more outbreed European populations. © 2016 Wiley Periodicals, Inc.

Do North Atlantic eels show parallel patterns of spatially varying selection?

BACKGROUND: The two North Atlantic eel species, the European and the American eel, represent an ideal system in which to study parallel selection patterns due to their sister species status and the presence of ongoing gene flow. A panel of 80 coding-gene SNPs previously analyzed in American eel was used to genotype European eel individuals (glass eels) from 8 sampling locations across the species distribution. We tested for single-generation signatures of spatially varying selection in European eel by searching for elevated genetic differentiation using FST-based outlier tests and by testing for significant associations between allele frequencies and environmental variables. RESULTS: We found signatures of possible selection at a total of 11 coding-gene SNPs. Candidate genes for local selection constituted mainly genes with a major role in metabolism as well as defense genes. Contrary to what has been found for American eel, only 2 SNPs in our study correlated with differences in temperature, which suggests that other explanatory variables may play a role. None of the genes found to be associated with explanatory variables in European eel showed any correlations with environmental factors in the previous study in American eel. CONCLUSIONS: The different signatures of selection between species could be due to distinct selective pressures associated with the much longer larval migration for European eel relative to American eel. The lack of parallel selection in North Atlantic eels could also be due to most phenotypic traits being polygenic, thus reducing the likelihood of selection acting on the same genes in both species.

Sex differences in shared genetic determinants between severe mental disorders and metabolic traits.

High rates of metabolic risk factors contribute to premature mortality in patients with severe mental disorders, but the molecular underpinnings of this association are largely unknown. We performed the first analysis on shared genetic factors between severe mental disorders and metabolic traits considering the effect of sex. We applied an integrated analytical pipeline on the largest sex-stratified genome-wide association datasets available for bipolar disorder (BD), major depressive disorder (MDD), schizophrenia (SZ), and for body mass index (BMI) and waist-to-hip ratio (WHR) (all including participants of European origin). We observed extensive genetic overlap between all severe mental disorders and variants associated with BMI in women or men and identified several genetic loci shared between BD, or SZ and BMI in women (24 and 91, respectively) or men (13 and 208, respectively), with mixed directions of effect. A large part of the identified genetic variants showed sex differences in terms of location, genes modulated in adipose tissue and/or brain regions, and druggable targets. By providing a complete picture of disorder specific and cross-disorder shared genetic determinants, our results highlight potential sex differences in the genetic liability to metabolic comorbidities in patients with severe mental disorders.

Distinct genomic signatures and modifiable risk factors underly the comorbidity between major depressive disorder and cardiovascular disease.

Major depressive disorder (MDD) and cardiovascular disease (CVD) are often comorbid, resulting in excess morbidity and mortality. Using genomic data, this study elucidates biological mechanisms, key risk factors, and causal pathways underlying their comorbidity. We show that CVDs share a large proportion of their genetic risk factors with MDD. Multivariate genome-wide association analysis of the shared genetic liability between MDD and atherosclerotic CVD (ASCVD) revealed seven novel loci and distinct patterns of tissue and brain cell-type enrichments, suggesting a role for the thalamus. Part of the genetic overlap was explained by shared inflammatory, metabolic, and psychosocial/lifestyle risk factors. Finally, we found support for causal effects of genetic liability to MDD on CVD risk, but not from most CVDs to MDD, and demonstrated that the causal effects were partly explained by metabolic and psychosocial/lifestyle factors. The distinct signature of MDD-ASCVD comorbidity aligns with the idea of an immunometabolic sub-type of MDD more strongly associated with CVD than overall MDD. In summary, we identify plausible biological mechanisms underlying MDD-CVD comorbidity, as well as key modifiable risk factors for prevention of CVD in individuals with MDD.

Distinct biological signature and modifiable risk factors underlie the comorbidity between major depressive disorder and cardiovascular disease.

Major depressive disorder (MDD) and cardiovascular disease (CVD) are often comorbid, resulting in excess morbidity and mortality. Here we show that CVDs share most of their genetic risk factors with MDD. Multivariate genome-wide association analysis of shared genetic liability between MDD and atherosclerotic CVD revealed seven loci and distinct patterns of tissue and brain cell-type enrichments, suggesting the involvement of the thalamus. Part of the genetic overlap was explained by shared inflammatory, metabolic and psychosocial or lifestyle risk factors. Our data indicated causal effects of genetic liability to MDD on CVD risk, but not from most CVDs to MDD, and showed that the causal effects were partly explained by metabolic and psychosocial or lifestyle factors. The distinct signature of MDD-atherosclerotic CVD comorbidity suggests an immunometabolic subtype of MDD that is more strongly associated with CVD than overall MDD. In summary, we identified biological mechanisms underlying MDD-CVD comorbidity and modifiable risk factors for prevention of CVD in individuals with MDD.

Distinct biological signature and modifiable risk factors underlie the comorbidity between major depressive disorder and cardiovascular disease.

Major depressive disorder (MDD) and cardiovascular disease (CVD) are often comorbid, resulting in excess morbidity and mortality. Here we show that CVDs share most of their genetic risk factors with MDD. Multivariate genome-wide association analysis of shared genetic liability between MDD and atherosclerotic CVD revealed seven loci and distinct patterns of tissue and brain cell-type enrichments, suggesting the involvement of the thalamus. Part of the genetic overlap was explained by shared inflammatory, metabolic and psychosocial or lifestyle risk factors. Our data indicated causal effects of genetic liability to MDD on CVD risk, but not from most CVDs to MDD, and showed that the causal effects were partly explained by metabolic and psychosocial or lifestyle factors. The distinct signature of MDD-atherosclerotic CVD comorbidity suggests an immunometabolic subtype of MDD that is more strongly associated with CVD than overall MDD. In summary, we identified biological mechanisms underlying MDD-CVD comorbidity and modifiable risk factors for prevention of CVD in individuals with MDD.

Genome-wide association study meta-analysis of 9,619 cases with tic disorders.

BACKGROUND: Despite the significant personal and societal burden of tic disorders (TD), treatment outcomes remain modest, necessitating a deeper understanding of their etiology. Family history is the biggest known risk factor and identifying risk genes could accelerate progress in the field. METHODS: Expanding upon previous sample size limitations, we added 4,800 new TD cases and 971,560 controls, conducting a GWAS meta-analysis with 9,619 cases and 981,048 controls of European ancestry. We attempted to replicate the results in an independent deCODE Genetics GWAS (885 TD cases and 310,367 controls). To characterize GWAS findings, we conducted several post-GWAS gene-based and enrichment analyses. RESULTS: A genome-wide significant hit (rs79244681, p=2.27x10-08) within MCHR2-AS1 was identified, though it was not replicated. Post-GWAS analyses revealed a 13.8% SNP-heritability and three significant genes: BCL11B, NDFIP2, and RBM26. Common variant risk for TD was enriched within genes preferentially expressed in the cortico-striato-thalamo-cortical circuit (including the putamen, caudate, nucleus accumbens, and Brodman area 9) and five brain cell types (excitatory and inhibitory telencephalon-, inhibitory- di- and mesencephalon, hindbrain-, and medium spiny neurons). TD polygenic risk was enriched within loss-of-function intolerant genes (p=0.0017) and high-confidence neurodevelopmental disorder genes (p=0.0108). Of 112 genetic correlations, 43 were statistically significant, showing high positive correlations with most psychiatric disorders. Of the two SNPs previously associated with TD, one (rs2453763) replicated in an independent sub-sample of our GWAS (p=0.00018). CONCLUSIONS: This GWAS was still underpowered to identify high-confidence, replicable loci, but the results suggest imminent discovery of common genetic variants for TD.

Genome-wide association study identifies new loci associated with OCD.

To date, four genome-wide association studies (GWAS) of obsessive-compulsive disorder (OCD) have been published, reporting a high single-nucleotide polymorphism (SNP)-heritability of 28% but finding only one significant SNP. A substantial increase in sample size will likely lead to further identification of SNPs, genes, and biological pathways mediating the susceptibility to OCD. We conducted a GWAS meta-analysis with a 2-3-fold increase in case sample size (OCD cases: N = 37,015, controls: N = 948,616) compared to the last OCD GWAS, including six previously published cohorts (OCGAS, IOCDF-GC, IOCDF-GC-trio, NORDiC-nor, NORDiC-swe, and iPSYCH) and unpublished self-report data from 23andMe Inc. We explored the genetic architecture of OCD by conducting gene-based tests, tissue and celltype enrichment analyses, and estimating heritability and genetic correlations with 74 phenotypes. To examine a potential heterogeneity in our data, we conducted multivariable GWASs with MTAG. We found support for 15 independent genome-wide significant loci (14 new) and 79 protein-coding genes. Tissue enrichment analyses implicate multiple cortical regions, the amygdala, and hypothalamus, while cell type analyses yielded 12 cell types linked to OCD (all neurons). The SNP-based heritability of OCD was estimated to be 0.08. Using MTAG we found evidence for specific genetic underpinnings characteristic of different cohort-ascertainment and identified additional significant SNPs. OCD was genetically correlated with 40 disorders or traits-positively with all psychiatric disorders and negatively with BMI, age at first birth and multiple autoimmune diseases. The GWAS meta-analysis identified several biologically informative genes as important contributors to the aetiology of OCD. Overall, we have begun laying the groundwork through which the biology of OCD will be understood and described.

Microevolution in time and space: SNP analysis of historical DNA reveals dynamic signatures of selection in Atlantic cod.

Little is known about how quickly natural populations adapt to changes in their environment and how temporal and spatial variation in selection pressures interact to shape patterns of genetic diversity. We here address these issues with a series of genome scans in four overfished populations of Atlantic cod (Gadus morhua) studied over an 80-year period. Screening of >1000 gene-associated single-nucleotide polymorphisms (SNPs) identified 77 loci that showed highly elevated levels of differentiation, likely as an effect of directional selection, in either time, space or both. Exploratory analysis suggested that temporal allele frequency shifts at certain loci may correlate with local temperature variation and with life history changes suggested to be fisheries induced. Interestingly, however, largely nonoverlapping sets of loci were temporal outliers in the different populations and outliers from the 1928 to 1960 period showed almost complete stability during later decades. The contrasting microevolutionary trajectories among populations resulted in sequential shifts in spatial outliers, with no locus maintaining elevated spatial differentiation throughout the study period. Simulations of migration coupled with observations of temporally stable spatial structure at neutral loci suggest that population replacement or gene flow alone could not explain all the observed allele frequency variation. Thus, the genetic changes are likely to at least partly be driven by highly dynamic temporally and spatially varying selection. These findings have important implications for our understanding of local adaptation and evolutionary potential in high gene flow organisms and underscore the need to carefully consider all dimensions of biocomplexity for evolutionarily sustainable management.

FarGen: Elucidating the distribution of coding variants in the isolated population of the Faroe Islands.

Here we present results from FarGen Phase I exomes. This dataset is based on the FarGen cohort, which consists of 1,541 individuals from the isolated population of the Faroe Islands. The purpose of this cohort is to serve as a reference catalog of coding variants, and to conduct population genetic studies to better understand the genetic contribution to various diseases in the Faroese population. The first whole-exome data set comprise 465 individuals and a total of 148,267 genetic variants were discovered. Principle Component Analysis indicates that the population is isolated and weakly structured. The distribution of variants in various functional classes was compared with populations in the gnomAD dataset; the results indicated that the proportions were consistent across the cohorts, but probably due to a small sample size, the FarGen dataset contained relatively few rare variants. We identified 19 variants that are classified as pathogenic or likely pathogenic in ClinVar; several of these variants are associated with monogenetic diseases with increased prevalence in the Faroe Islands. The results support previous studies, which indicate that the Faroe Islands is an isolated and weakly structured population. Future studies may elucidate the significance of the 19 pathogenic variants that were identified. The FarGen Phase I dataset is an important step for genetic research in the Faroese population, and the next phase of FarGen will increase the sample size and broaden the scope.

Genome-wide association study of school grades identifies genetic overlap between language ability, psychopathology and creativity.

Cognitive functions of individuals with psychiatric disorders differ from that of the general population. Such cognitive differences often manifest early in life as differential school performance and have a strong genetic basis. Here we measured genetic predictors of school performance in 30,982 individuals in English, Danish and mathematics via a genome-wide association study (GWAS) and studied their relationship with risk for six major psychiatric disorders. When decomposing the school performance into math and language-specific performances, we observed phenotypically and genetically a strong negative correlation between math performance and risk for most psychiatric disorders. But language performance correlated positively with risk for certain disorders, especially schizophrenia, which we replicate in an independent sample (n = 4547). We also found that the genetic variants relating to increased risk for schizophrenia and better language performance are overrepresented in individuals involved in creative professions (n = 2953) compared to the general population (n = 164,622). The findings together suggest that language ability, creativity and psychopathology might stem from overlapping genetic roots.

Genome-wide analyses of ADHD identify 27 risk loci, refine the genetic architecture and implicate several cognitive domains.

Attention-deficit hyperactivity disorder (ADHD) is a prevalent neurodevelopmental disorder with a major genetic component. Here, we present a genome-wide association study meta-analysis of ADHD comprising 38,691 individuals with ADHD and 186,843 controls. We identified 27 genome-wide significant loci, highlighting 76 potential risk genes enriched among genes expressed particularly in early brain development. Overall, ADHD genetic risk was associated with several brain-specific neuronal subtypes and midbrain dopaminergic neurons. In exome-sequencing data from 17,896 individuals, we identified an increased load of rare protein-truncating variants in ADHD for a set of risk genes enriched with probable causal common variants, potentially implicating SORCS3 in ADHD by both common and rare variants. Bivariate Gaussian mixture modeling estimated that 84-98% of ADHD-influencing variants are shared with other psychiatric disorders. In addition, common-variant ADHD risk was associated with impaired complex cognition such as verbal reasoning and a range of executive functions, including attention.

Depression pathophysiology, risk prediction of recurrence and comorbid psychiatric disorders using genome-wide analyses.

Depression is a common psychiatric disorder and a leading cause of disability worldwide. Here we conducted a genome-wide association study meta-analysis of six datasets, including >1.3 million individuals (371,184 with depression) and identified 243 risk loci. Overall, 64 loci were new, including genes encoding glutamate and GABA receptors, which are targets for antidepressant drugs. Intersection with functional genomics data prioritized likely causal genes and revealed new enrichment of prenatal GABAergic neurons, astrocytes and oligodendrocyte lineages. We found depression to be highly polygenic, with ~11,700 variants explaining 90% of the single-nucleotide polymorphism heritability, estimating that >95% of risk variants for other psychiatric disorders (anxiety, schizophrenia, bipolar disorder and attention deficit hyperactivity disorder) were influencing depression risk when both concordant and discordant variants were considered, and nearly all depression risk variants influenced educational attainment. Additionally, depression genetic risk was associated with impaired complex cognition domains. We dissected the genetic and clinical heterogeneity, revealing distinct polygenic architectures across subgroups of depression and demonstrating significantly increased absolute risks for recurrence and psychiatric comorbidity among cases of depression with the highest polygenic burden, with considerable sex differences. The risks were up to 5- and 32-fold higher than cases with the lowest polygenic burden and the background population, respectively. These results deepen the understanding of the biology underlying depression, its disease progression and inform precision medicine approaches to treatment.