Molecular interpretation of genome-wide association studies using multiomics analysis

Genome-wide association studies have found thousands of single-nucleotide polymorphisms associated with various human traits. Recently, a powerful statistical approach called MetaXcan has been proposed for interpreting genome-wide associations at the gene level. We extended MetaXcan to a multiomics application, using a brain cortex reference dataset that includes gene expression, DNA methylation, and histone acetylation data from approximately 400 individuals. Our approach, Multi-MetaXcan, consists of three steps. In the first step, we use regularized regression to build models that predict gene expression and variation in epigenomic modifications from single-nucleotide polymorphisms. We call these models genotype-based imputation models. In the second step, we apply these models to map genome-wide associations to gene-level and epigenomic-level associations. Finally, in the third step, our model summarizes all molecular-level associations at the gene level by building epigenome-based imputation models that predict gene expression levels from nearby epigenomic marks like CpG sites and transcriptionally active regions. In summary, Multi-MetaXcan identifies trait-associated genes whose expression levels are impacted by single-nucleotide polymorphisms and their influence on intermediate molecular traits such as DNA methylation and histone acetylation. We applied Multi-MetaXcan to a major depressive disorder genome-wide association study. As the result, we discovered 12 genes, 25 transcriptionally active regions, and 163 CpG sites associated with major depressive disorder corresponding to 74 genes in total. 26 of these genes fall within or close to previously identified major depressive disorder-associated genomic regions. Importantly, the inclusion of epigenomic data resulted in an additional 62 genes that were not identified by gene expression imputation model alone.Science, Faculty ofGraduat

Genome-wide interaction study of a proxy for stress-sensitivity and its prediction of major depressive disorder

Individual response to stress is correlated with neuroticism and is an important predictor of both neuroticism and the onset of major depressive disorder (MDD). Identification of the genetics underpinning individual differences in response to negative events (stress-sensitivity) may improve our understanding of the molecular pathways involved, and its association with stress-related illnesses. We sought to generate a proxy for stress-sensitivity through modelling the interaction between SNP allele and MDD status on neuroticism score in order to identify genetic variants that contribute to the higher neuroticism seen in individuals with a lifetime diagnosis of depression compared to unaffected individuals. Meta-analysis of genome-wide interaction studies (GWIS) in UK Biobank (N = 23,092) and Generation Scotland: Scottish Family Health Study (N = 7,155) identified no genome-wide significance SNP interactions. However, gene-based tests identified a genome-wide significant gene, ZNF366, a negative regulator of glucocorticoid receptor function implicated in alcohol dependence (p = 1.48x10-7; Bonferroni-corrected significance threshold p < 2.79x10-6). Using summary statistics from the stress-sensitivity term of the GWIS, SNP heritability for stress-sensitivity was estimated at 5.0%. In models fitting polygenic risk scores of both MDD and neuroticism derived from independent GWAS, we show that polygenic risk scores derived from the UK Biobank stress-sensitivity GWIS significantly improved the prediction of MDD in Generation Scotland. This study may improve interpretation of larger genome-wide association studies of MDD and other stress-related illnesses, and the understanding of the etiological mechanisms underpinning stress-sensitivity

Dissecting the Shared Genetic Architecture of Suicide Attempt, Psychiatric Disorders, and Known Risk Factors

Background Suicide is a leading cause of death worldwide, and nonfatal suicide attempts, which occur far more frequently, are a major source of disability and social and economic burden. Both have substantial genetic etiology, which is partially shared and partially distinct from that of related psychiatric disorders. Methods We conducted a genome-wide association study (GWAS) of 29,782 suicide attempt (SA) cases and 519,961 controls in the International Suicide Genetics Consortium (ISGC). The GWAS of SA was conditioned on psychiatric disorders using GWAS summary statistics via multitrait-based conditional and joint analysis, to remove genetic effects on SA mediated by psychiatric disorders. We investigated the shared and divergent genetic architectures of SA, psychiatric disorders, and other known risk factors. Results Two loci reached genome-wide significance for SA: the major histocompatibility complex and an intergenic locus on chromosome 7, the latter of which remained associated with SA after conditioning on psychiatric disorders and replicated in an independent cohort from the Million Veteran Program. This locus has been implicated in risk-taking behavior, smoking, and insomnia. SA showed strong genetic correlation with psychiatric disorders, particularly major depression, and also with smoking, pain, risk-taking behavior, sleep disturbances, lower educational attainment, reproductive traits, lower socioeconomic status, and poorer general health. After conditioning on psychiatric disorders, the genetic correlations between SA and psychiatric disorders decreased, whereas those with nonpsychiatric traits remained largely unchanged. Conclusions Our results identify a risk locus that contributes more strongly to SA than other phenotypes and suggest a shared underlying biology between SA and known risk factors that is not mediated by psychiatric disorders.Peer reviewe

Evidence for Increased Genetic Risk Load for Major Depression in Patients Assigned to Electroconvulsive Therapy

Electroconvulsive therapy (ECT) is the treatment of choice for severe and treatment-resistant depression; disorder severity and unfavorable treatment outcomes are shown to be influenced by an increased genetic burden for major depression (MD). Here, we tested whether ECT assignment and response/nonresponse are associated with an increased genetic burden for major depression (MD) using polygenic risk score (PRS), which summarize the contribution of diseaserelated common risk variants. Fifty-one psychiatric inpatients suffering from a major depressive episode underwent ECT. MD-PRS were calculated for these inpatients and a separate population-based sample (n = 3,547 healthy; n = 426 self-reported depression) based on summary statistics from the Psychiatric Genomics Consortium MDD-working group (Cases: n = 59,851; Controls: n = 113,154). MD-PRS explained a significant proportion of disease status between ECT patients and healthy controls (p = .022, R2 = 1.173%); patients showed higher MD-PRS. MD-PRS in population-based depression self-reporters were intermediate between ECT patients and controls (n.s.). Significant associations between MD-PRS and ECT response (50% reduction in Hamilton depression rating scale scores) were not observed. Our findings indicate that ECT cohorts show an increased genetic burden for MD and are consistent with the hypothesis that treatment-resistant MD patients represent a subgroup with an increased genetic risk for MD. Larger samples are needed to better substantiate these findings

A phenome-wide association and Mendelian Randomisation study of polygenic risk for depression in UK Biobank.

Depression is a leading cause of worldwide disability but there remains considerable uncertainty regarding its neural and behavioural associations. Here, using non-overlapping Psychiatric Genomics Consortium (PGC) datasets as a reference, we estimate polygenic risk scores for depression (depression-PRS) in a discovery (N = 10,674) and replication (N = 11,214) imaging sample from UK Biobank. We report 77 traits that are significantly associated with depression-PRS, in both discovery and replication analyses. Mendelian Randomisation analysis supports a potential causal effect of liability to depression on brain white matter microstructure (β: 0.125 to 0.868, pFDR < 0.043). Several behavioural traits are also associated with depression-PRS (β: 0.014 to 0.180, pFDR: 0.049 to 1.28 × 10-14) and we find a significant and positive interaction between depression-PRS and adverse environmental exposures on mental health outcomes. This study reveals replicable associations between depression-PRS and white matter microstructure. Our results indicate that white matter microstructure differences may be a causal consequence of liability to depression

Association of whole-genome and NETRIN1 signaling pathway-derived polygenic risk scores for Major Depressive Disorder and white matter microstructure in UK Biobank

Background: Major depressive disorder is a clinically heterogeneous psychiatric disorder with a polygenic architecture. Genome-wide association studies have identified a number of risk-associated variants across the genome and have reported growing evidence of NETRIN1 pathway involvement. Stratifying disease risk by genetic variation within the NETRIN1 pathway may provide important routes for identification of disease mechanisms by focusing on a specific process, excluding heterogeneous risk-associated variation in other pathways. Here, we sought to investigate whether major depressive disorder polygenic risk scores derived from the NETRIN1 signaling pathway (NETRIN1-PRSs) and the whole genome, excluding NETRIN1 pathway genes (genomic-PRSs), were associated with white matter microstructure. Methods: We used two diffusion tensor imaging measures, fractional anisotropy (FA) and mean diffusivity (MD), in the most up-to-date UK Biobank neuroimaging data release (FA: n = 6401; MD: n = 6390). Results: We found significantly lower FA in the superior longitudinal fasciculus (β = −.035, p =.029) and significantly higher MD in a global measure of thalamic radiations (β =.029, p =.021), as well as higher MD in the superior (β =.034, p =.039) and inferior (β =.029, p =.043) longitudinal fasciculus and in the anterior (β =.025, p =.046) and superior (β =.027, p =.043) thalamic radiation associated with NETRIN1-PRS. Genomic-PRS was also associated with lower FA and higher MD in several tracts. Conclusions: Our findings indicate that variation in the NETRIN1 signaling pathway may confer risk for major depressive disorder through effects on a number of white matter tracts

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.Peer reviewe

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

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

The genetics of the mood disorder spectrum:genome-wide association analyses of over 185,000 cases and 439,000 controls

Background Mood disorders (including major depressive disorder and bipolar disorder) affect 10-20% of the population. They range from brief, mild episodes to severe, incapacitating conditions that markedly impact lives. Despite their diagnostic distinction, multiple approaches have shown considerable sharing of risk factors across the mood disorders. Methods To clarify their shared molecular genetic basis, and to highlight disorder-specific associations, we meta-analysed data from the latest Psychiatric Genomics Consortium (PGC) genome-wide association studies of major depression (including data from 23andMe) and bipolar disorder, and an additional major depressive disorder cohort from UK Biobank (total: 185,285 cases, 439,741 controls; non-overlapping N = 609,424). Results Seventy-three loci reached genome-wide significance in the meta-analysis, including 15 that are novel for mood disorders. More genome-wide significant loci from the PGC analysis of major depression than bipolar disorder reached genome-wide significance. Genetic correlations revealed that type 2 bipolar disorder correlates strongly with recurrent and single episode major depressive disorder. Systems biology analyses highlight both similarities and differences between the mood disorders, particularly in the mouse brain cell-types implicated by the expression patterns of associated genes. The mood disorders also differ in their genetic correlation with educational attainment – positive in bipolar disorder but negative in major depressive disorder. Conclusions The mood disorders share several genetic associations, and can be combined effectively to increase variant discovery. However, we demonstrate several differences between these disorders. Analysing subtypes of major depressive disorder and bipolar disorder provides evidence for a genetic mood disorders spectrum