Genome-wide Burden of Rare Short Deletions is Enriched in Major Depressive Disorder in Four Cohorts

Major depressive disorder (MDD) is moderately heritable, with a high prevalence and a presumed high heterogeneity. Copy number variants (CNVs) could contribute to the heritable component of risk, but the two previous genome-wide association studies of rare CNVs did not report significant findings. In this meta-analysis of four cohorts (5780 patients and 6626 control subjects), we analyzed the association of MDD to 1) genome-wide burden of rare deletions and duplications, partitioned by length (100 kb) and other characteristics, and 2) individual rare exonic CNVs and CNV regions. Patients with MDD carried significantly more short deletions than control subjects (p = .0059) but not long deletions or short or long duplications. The confidence interval for long deletions overlapped with that for short deletions, but long deletions were 70% less frequent genome-wide, reducing the power to detect increased burden. The increased burden of short deletions was primarily in intergenic regions. Short deletions in cases were also modestly enriched for high-confidence enhancer regions. No individual CNV achieved thresholds for suggestive or significant association after genome-wide correction. p values < .01 were observed for 15q11.2 duplications (TUBGCP5, CYFIP1, NIPA1, and NIPA2), deletions in or near PRKN or MSR1, and exonic duplications of ATG5. The increased burden of short deletions in patients with MDD suggests that rare CNVs increase the risk of MDD by disrupting regulatory regions. Results for longer deletions were less clear, but no large effects were observed for long multigenic CNVs (as seen in schizophrenia and autism). Further studies with larger sample sizes are warranted. [Abstract copyright: Copyright © 2019 Society of Biological Psychiatry. All rights reserved.

Identification of common genetic risk variants for autism spectrum disorder

Autism spectrum disorder (ASD) is a highly heritable and heterogeneous group of neurodevelopmental phenotypes diagnosed in more than 1% of children. Common genetic variants contribute substantially to ASD susceptibility, but to date no individual variants have been robustly associated with ASD. With a marked sample-size increase from a unique Danish population resource, we report a genome-wide association meta-analysis of 18,381 individuals with ASD and 27,969 controls that identified five genome-wide-significant loci. Leveraging GWAS results from three phenotypes with significantly overlapping genetic architectures (schizophrenia, major depression, and educational attainment), we identified seven additional loci shared with other traits at equally strict significance levels. Dissecting the polygenic architecture, we found both quantitative and qualitative polygenic heterogeneity across ASD subtypes. These results highlight biological insights, particularly relating to neuronal function and corticogenesis, and establish that GWAS performed at scale will be much more productive in the near term in ASD.Peer reviewe

Genome-wide association study identifies 30 Loci Associated with Bipolar Disorder

This paper is dedicated to the memory of Psychiatric Genomics Consortium (PGC) founding member and Bipolar disorder working group co-chair Pamela Sklar. We thank the participants who donated their time, experiences and DNA to this research, and to the clinical and scientific teams that worked with them. We are deeply indebted to the investigators who comprise the PGC. The views expressed are those of the authors and not necessarily those of any funding or regulatory body. Analyses were carried out on the NL Genetic Cluster Computer (http://www.geneticcluster.org ) hosted by SURFsara, and the Mount Sinai high performance computing cluster (http://hpc.mssm.edu).Bipolar disorder is a highly heritable psychiatric disorder. We performed a genome-wide association study including 20,352 cases and 31,358 controls of European descent, with follow-up analysis of 822 variants with P<1x10-4 in an additional 9,412 cases and 137,760 controls. Eight of the 19 variants that were genome-wide significant (GWS, p < 5x10-8) in the discovery GWAS were not GWS in the combined analysis, consistent with small effect sizes and limited power but also with genetic heterogeneity. In the combined analysis 30 loci were GWS including 20 novel loci. The significant loci contain genes encoding ion channels, neurotransmitter transporters and synaptic components. Pathway analysis revealed nine significantly enriched gene-sets including regulation of insulin secretion and endocannabinoid signaling. BDI is strongly genetically correlated with schizophrenia, driven by psychosis, whereas BDII is more strongly correlated with major depressive disorder. These findings address key clinical questions and provide potential new biological mechanisms for BD.This work was funded in part by the Brain and Behavior Research Foundation, Stanley Medical Research Institute, University of Michigan, Pritzker Neuropsychiatric Disorders Research Fund L.L.C., Marriot Foundation and the Mayo Clinic Center for Individualized Medicine, the NIMH Intramural Research Program; Canadian Institutes of Health Research; the UK Maudsley NHS Foundation Trust, NIHR, NRS, MRC, Wellcome Trust; European Research Council; German Ministry for Education and Research, German Research Foundation IZKF of Münster, Deutsche Forschungsgemeinschaft, ImmunoSensation, the Dr. Lisa-Oehler Foundation, University of Bonn; the Swiss National Science Foundation; French Foundation FondaMental and ANR; Spanish Ministerio de Economía, CIBERSAM, Industria y Competitividad, European Regional Development Fund (ERDF), Generalitat de Catalunya, EU Horizon 2020 Research and Innovation Programme; BBMRI-NL; South-East Norway Regional Health Authority and Mrs. Throne-Holst; Swedish Research Council, Stockholm County Council, Söderström Foundation; Lundbeck Foundation, Aarhus University; Australia NHMRC, NSW Ministry of Health, Janette M O'Neil and Betty C Lynch

Improving Genetic Prediction by Leveraging Genetic Correlations Among Human Diseases and Traits

Genomic prediction has the potential to contribute to precision medicine. However, to date, the utility of such predictors is limited due to low accuracy for most traits. Here theory and simulation study are used to demonstrate that widespread pleiotropy among phenotypes can be utilised to improve genomic risk prediction. We show how a genetic predictor can be created as a weighted index that combines published genome-wide association study (GWAS) summary statistics across many different traits. We apply this framework to predict risk of schizophrenia and bipolar disorder in the Psychiatric Genomics consortium data, finding substantial heterogeneity in prediction accuracy increases across cohorts. For six additional phenotypes in the UK Biobank data, we find increases in prediction accuracy ranging from 0.7 for height to 47 for type 2 diabetes, when using a multi-trait predictor that combines published summary statistics from multiple traits, as compared to a predictor based only on one trait. © 2018 The Author(s)

Classical Human Leukocyte Antigen Alleles and C4 Haplotypes Are Not Significantly Associated With Depression

Background: The prevalence of depression is higher in individuals with autoimmune diseases, but the mechanisms underlying the observed comorbidities are unknown. Shared genetic etiology is a plausible explanation for the overlap, and in this study we tested whether genetic variation in the major histocompatibility complex (MHC), which is associated with risk for autoimmune diseases, is also associated with risk for depression. Methods: We fine-mapped the classical MHC (chr6: 29.6–33.1 Mb), imputing 216 human leukocyte antigen (HLA) alleles and 4 complement component 4 (C4) haplotypes in studies from the Psychiatric Genomics Consortium Major Depressive Disorder Working Group and the UK Biobank. The total sample size was 45,149 depression cases and 86,698 controls. We tested for association between depression status and imputed MHC variants, applying both a region-wide significance threshold (3.9 × 10−6) and a candidate threshold (1.6 × 10−4). Results: No HLA alleles or C4 haplotypes were associated with depression at the region-wide threshold. HLA-B*08:01 was associated with modest protection for depression at the candidate threshold for testing in HLA genes in the meta-analysis (odds ratio = 0.98, 95% confidence interval = 0.97–0.99). Conclusions: We found no evidence that an increased risk for depression was conferred by HLA alleles, which play a major role in the genetic susceptibility to autoimmune diseases, or C4 haplotypes, which are strongly associated with schizophrenia. These results suggest that any HLA or C4 variants associated with depression either are rare or have very modest effect sizes

Investigating rare pathogenic/likely pathogenic exonic variation in bipolar disorder

Bipolar disorder (BD) is a serious mental illness with substantial common variant heritability. However, the role of rare coding variation in BD is not well established. We examined the protein-coding (exonic) sequences of 3,987 unrelated individuals with BD and 5,322 controls of predominantly European ancestry across four cohorts from the Bipolar Sequencing Consortium (BSC). We assessed the burden of rare, protein-altering, single nucleotide variants classified as pathogenic or likely pathogenic (P-LP) both exome-wide and within several groups of genes with phenotypic or biologic plausibility in BD. While we observed an increased burden of rare coding P-LP variants within 165 genes identified as BD GWAS regions in 3,987 BD cases (meta-analysis OR = 1.9, 95% CI = 1.3–2.8, one-sided p = 6.0 × 10−4), this enrichment did not replicate in an additional 9,929 BD cases and 14,018 controls (OR = 0.9, one-side p = 0.70). Although BD shares common variant heritability with schizophrenia, in the BSC sample we did not observe a significant enrichment of P-LP variants in SCZ GWAS genes, in two classes of neuronal synaptic genes (RBFOX2 and FMRP) associated with SCZ or in loss-of-function intolerant genes. In this study, the largest analysis of exonic variation in BD, individuals with BD do not carry a replicable enrichment of rare P-LP variants across the exome or in any of several groups of genes with biologic plausibility. Moreover, despite a strong shared susceptibility between BD and SCZ through common genetic variation, we do not observe an association between BD risk and rare P-LP coding variants in genes known to modulate risk for SCZ

Genome-wide by Environment Interaction Studies of Depressive Symptoms and Psychosocial Stress in UK Biobank and Generation Scotland

Stress is associated with poorer physical and mental health. To improve our understanding of this link, we performed genome-wide association studies (GWAS) of depressive symptoms and genome-wide by environment interaction studies (GWEIS) of depressive symptoms and stressful life events (SLE) in two UK population-based cohorts (Generation Scotland and UK Biobank). No SNP was individually significant in either GWAS, but gene-based tests identified six genes associated with depressive symptoms in UK Biobank (DCC, ACSS3, DRD2, STAG1, FOXP2 and KYNU; p < 2.77 x 10(-6)). Two SNPs with genome-wide significant GxE effects were identified by GWEIS in Generation Scotland: rs12789145 (53-kb downstream PIWIL4; p = 4.95 x 10(-9); total SLE) and rs17070072 (intronic to ZCCHC2; p = 1.46 x 10(-8); dependent SLE). A third locus upstream CYLC2 (rs12000047 and rs12005200, p < 2.00 x 10(-8); dependent SLE) when the joint effect of the SNP main and GxE effects was considered. GWEIS gene-based tests identified: MTNR1B with GxE effect with dependent SLE in Generation Scotland; and PHF2 with the joint effect in UK Biobank (p < 2.77 x 10(-6)). Polygenic risk scores (PRSs) analyses incorporating GxE effects improved the prediction of depressive symptom scores, when using weights derived from either the UK Biobank GWAS of depressive symptoms (p = 0.01) or the PGC GWAS of major depressive disorder (p = 5.91 x 10(-3)). Using an independent sample, PRS derived using GWEIS GxE effects provided evidence of shared aetiologies between depressive symptoms and schizotypal personality, heart disease and COPD. Further such studies are required and may result in improved treatments for depression and other stress-related conditions

Integrated analysis of environmental and genetic influences on cord blood DNA methylation in new-borns

Epigenetic processes, including DNA methylation (DNAm), are among the mechanisms allowing integration of genetic and environmental factors to shape cellular function. While many studies have investigated either environmental or genetic contributions to DNAm, few have assessed their integrated effects. Here we examine the relative contributions of prenatal environmental factors and genotype on DNA methylation in neonatal blood at variably methylated regions (VMRs) in 4 independent cohorts (overall n = 2365). We use Akaike’s information criterion to test which factors best explain variability of methylation in the cohort-specific VMRs: several prenatal environmental factors (E), genotypes in cis (G), or their additive (G + E) or interaction (GxE) effects. Genetic and environmental factors in combination best explain DNAm at the majority of VMRs. The CpGs best explained by either G, G + E or GxE are functionally distinct. The enrichment of genetic variants from GxE models in GWAS for complex disorders supports their importance for disease risk

Genome-wide association of mood-incongruent psychotic bipolar disorder

Mood-incongruent psychotic features (MICP) are familial symptoms of bipolar disorder (BP) that also occur in schizophrenia (SZ), and may represent manifestations of shared etiology between the major psychoses. In this study we have analyzed three large samples of BP with imputed genome-wide association data and have performed a meta-analysis of 2196 cases with MICP and 8148 controls. We found several regions with suggestive evidence of association (P<10(-6)), although no marker met genome-wide significance criteria. The top associations were on chromosomes: 6q14.2 within the PRSS35/SNAP91 gene complex (rs1171113, P=9.67 × 10(-8)); 3p22.2 downstream of TRANK/LBA1 (rs9834970, P=9.71 × 10(-8)); and 14q24.2 in an intron of NUMB (rs2333194, P=7.03 × 10(-7)). These associations were present in all three samples, and both rs1171113 and rs2333194 were found to be overrepresented in an analysis of MICP cases compared with all other BP cases. To test the relationship of MICP with SZ, we performed polygenic analysis using the Psychiatric GWAS Consortium SZ results and found evidence of association between SZ polygenes and the presence of MICP in BP cases (meta-analysis P=0.003). In summary, our analysis of the MICP phenotype in BP has provided suggestive evidence for association of common variants in several genes expressed in the nervous system. The results of our polygenic analysis provides support for a modest degree of genetic overlap between BP with MICP and SZ, highlighting that phenotypic correlations across syndromes may be due to the influence of polygenic risk factors