Leveraging pleiotropy for the improved treatment of psychiatric disorders.

Over 90% of drug candidates fail in clinical trials, while it takes 10-15 years and one billion US dollars to develop a single successful drug. Drug development is more challenging for psychiatric disorders, where disease comorbidity and complex symptom profiles obscure the identification of causal mechanisms for therapeutic intervention. One promising approach for determining more suitable drug candidates in clinical trials is integrating human genetic data into the selection process. Genome-wide association studies have identified thousands of replicable risk loci for psychiatric disorders, and sophisticated statistical tools are increasingly effective at using these data to pinpoint likely causal genes. These studies have also uncovered shared or pleiotropic genetic risk factors underlying comorbid psychiatric disorders. In this article, we argue that leveraging pleiotropic effects will provide opportunities to discover novel drug targets and identify more effective treatments for psychiatric disorders by targeting a common mechanism rather than treating each disease separately.

The genetic architecture of substance use and its diverse correlations with mental health traits.

Although harmful substance use is common and represented by shared symptom features and high genetic correlations, the underlying genetic relationships between substance use traits have not been fully explored. We have investigated the genetic architecture of substance use traits through exploratory and confirmatory factor analyses using genomic structural equation modeling (Genomic SEM), and explored genetic correlations between different aspects of substance use and mental health-related traits. Genomic SEM was used to identify latent factors representing the relationships between 14 substance use traits (alcohol, nicotine, cannabis and opioid use), and to confirm or modify existing latent factors for 38 mental health-related traits. A bi-factor model best explained the genetic overlap between substance use traits, including a general substance use factor and two sub-factors representing genetic liability specific to alcohol use or smoking. The SNP-based heritability of these factors ranged from 2 to 7 % and each factor had 10 or more independent significant SNPs identified. Bivariate correlations revealed patterns of genetic overlap with other mental health-related factors unique to each substance use factor. Variations in the genetic overlap between psychiatric traits and different aspects of substance use can be used to further investigate the pleiotropy present between these traits, and explore commonalities in etiology.

The genetic landscape of substance use disorders.

Substance use disorders represent a significant public health concern with considerable socioeconomic implications worldwide. Twin and family-based studies have long established a heritable component underlying these disorders. In recent years, genome-wide association studies of large, broadly phenotyped samples have identified regions of the genome that harbour genetic risk variants associated with substance use disorders. These regions have enabled the discovery of putative causal genes and improved our understanding of genetic relationships among substance use disorders and other traits. Furthermore, the integration of these data with clinical information has yielded promising insights into how individuals respond to medications, allowing for the development of personalized treatment approaches based on an individual's genetic profile. This review article provides an overview of recent advances in the genetics of substance use disorders and demonstrates how genetic data may be used to reduce the burden of disease and improve public health outcomes.

Genome-wide association study identifies 30 obsessive-compulsive disorder associated loci.

Obsessive-compulsive disorder (OCD) affects ~1% of the population and exhibits a high SNP-heritability, yet previous genome-wide association studies (GWAS) have provided limited information on the genetic etiology and underlying biological mechanisms of the disorder. We conducted a GWAS meta-analysis combining 53,660 OCD cases and 2,044,417 controls from 28 European-ancestry cohorts revealing 30 independent genome-wide significant SNPs and a SNP-based heritability of 6.7%. Separate GWAS for clinical, biobank, comorbid, and self-report sub-groups found no evidence of sample ascertainment impacting our results. Functional and positional QTL gene-based approaches identified 249 significant candidate risk genes for OCD, of which 25 were identified as putatively causal, highlighting WDR6, DALRD3, CTNND1 and genes in the MHC region. Tissue and single-cell enrichment analyses highlighted hippocampal and cortical excitatory neurons, along with D1- and D2-type dopamine receptor-containing medium spiny neurons, as playing a role in OCD risk. OCD displayed significant genetic correlations with 65 out of 112 examined phenotypes. Notably, it showed positive genetic correlations with all included psychiatric phenotypes, in particular anxiety, depression, anorexia nervosa, and Tourette syndrome, and negative correlations with a subset of the included autoimmune disorders, educational attainment, and body mass index.. This study marks a significant step toward unraveling its genetic landscape and advances understanding of OCD genetics, providing a foundation for future interventions to address this debilitating disorder.

Genetic influences on circulating retinol and its relationship to human health.

Retinol is a fat-soluble vitamin that plays an essential role in many biological processes throughout the human lifespan. Here, we perform the largest genome-wide association study (GWAS) of retinol to date in up to 22,274 participants. We identify eight common variant loci associated with retinol, as well as a rare-variant signal. An integrative gene prioritisation pipeline supports novel retinol-associated genes outside of the main retinol transport complex (RBP4:TTR) related to lipid biology, energy homoeostasis, and endocrine signalling. Genetic proxies of circulating retinol were then used to estimate causal relationships with almost 20,000 clinical phenotypes via a phenome-wide Mendelian randomisation study (MR-pheWAS). The MR-pheWAS suggests that retinol may exert causal effects on inflammation, adiposity, ocular measures, the microbiome, and MRI-derived brain phenotypes, amongst several others. Conversely, circulating retinol may be causally influenced by factors including lipids and serum creatinine. Finally, we demonstrate how a retinol polygenic score could identify individuals more likely to fall outside of the normative range of circulating retinol for a given age. In summary, this study provides a comprehensive evaluation of the genetics of circulating retinol, as well as revealing traits which should be prioritised for further investigation with respect to retinol related therapies or nutritional intervention.

Genome-wide association study meta-analysis of dizygotic twinning illuminates genetic regulation of female fecundity.

STUDY QUESTION: Which genetic factors regulate female propensity for giving birth to spontaneous dizygotic (DZ) twins? SUMMARY ANSWER: We identified four new loci, GNRH1, FSHR, ZFPM1, and IPO8, in addition to previously identified loci, FSHB and SMAD3. WHAT IS KNOWN ALREADY: The propensity to give birth to DZ twins runs in families. Earlier, we reported that FSHB and SMAD3 as associated with DZ twinning and female fertility measures. STUDY DESIGN, SIZE, DURATION: We conducted a genome-wide association meta-analysis (GWAMA) of mothers of spontaneous dizygotic (DZ) twins (8265 cases, 264 567 controls) and of independent DZ twin offspring (26 252 cases, 417 433 controls). PARTICIPANTS/MATERIALS, SETTING, METHODS: Over 700 000 mothers of DZ twins, twin individuals and singletons from large cohorts in Australia/New Zealand, Europe, and the USA were carefully screened to exclude twins born after use of ARTs. Genetic association analyses by cohort were followed by meta-analysis, phenome wide association studies (PheWAS), in silico and in vivo annotations, and Zebrafish functional validation. MAIN RESULTS AND THE ROLE OF CHANCE: This study enlarges the sample size considerably from previous efforts, finding four genome-wide significant loci, including two novel signals and a further two novel genes that are implicated by gene level enrichment analyses. The novel loci, GNRH1 and FSHR, have well-established roles in female reproduction whereas ZFPM1 and IPO8 have not previously been implicated in female fertility. We found significant genetic correlations with multiple aspects of female reproduction and body size as well as evidence for significant selection against DZ twinning during human evolution. The 26 top single nucleotide polymorphisms (SNPs) from our GWAMA in European-origin participants weakly predicted the crude twinning rates in 47 non-European populations (r = 0.23 between risk score and population prevalence, s.e. 0.11, 1-tail P = 0.058) indicating that genome-wide association studies (GWAS) are needed in African and Asian populations to explore the causes of their respectively high and low DZ twinning rates. In vivo functional tests in zebrafish for IPO8 validated its essential role in female, but not male, fertility. In most regions, risk SNPs linked to known expression quantitative trait loci (eQTLs). Top SNPs were associated with in vivo reproductive hormone levels with the top pathways including hormone ligand binding receptors and the ovulation cycle. LARGE SCALE DATA: The full DZT GWAS summary statistics will made available after publication through the GWAS catalog (https://www.ebi.ac.uk/gwas/). LIMITATIONS, REASONS FOR CAUTION: Our study only included European ancestry cohorts. Inclusion of data from Africa (with the highest twining rate) and Asia (with the lowest rate) would illuminate further the biology of twinning and female fertility. WIDER IMPLICATIONS OF THE FINDINGS: About one in 40 babies born in the world is a twin and there is much speculation on why twinning runs in families. We hope our results will inform investigations of ovarian response in new and existing ARTs and the causes of female infertility. STUDY FUNDING/COMPETING INTEREST(S): Support for the Netherlands Twin Register came from the Netherlands Organization for Scientific Research (NWO) and The Netherlands Organization for Health Research and Development (ZonMW) grants, 904-61-193, 480-04-004, 400-05-717, Addiction-31160008, 911-09-032, Biobanking and Biomolecular Resources Research Infrastructure (BBMRI.NL, 184.021.007), Royal Netherlands Academy of Science Professor Award (PAH/6635) to DIB, European Research Council (ERC-230374), Rutgers University Cell and DNA Repository (NIMH U24 MH068457-06), the Avera Institute, Sioux Falls, South Dakota (USA) and the National Institutes of Health (NIH R01 HD042157-01A1) and the Genetic Association Information Network (GAIN) of the Foundation for the National Institutes of Health and Grand Opportunity grants 1RC2 MH089951. The QIMR Berghofer Medical Research Institute (QIMR) study was supported by grants from the National Health and Medical Research Council (NHMRC) of Australia (241944, 339462, 389927, 389875, 389891, 389892, 389938, 443036, 442915, 442981, 496610, 496739, 552485, 552498, 1050208, 1075175). L.Y. is funded by Australian Research Council (Grant number DE200100425). The Minnesota Center for Twin and Family Research (MCTFR) was supported in part by USPHS Grants from the National Institute on Alcohol Abuse and Alcoholism (AA09367 and AA11886) and the National Institute on Drug Abuse (DA05147, DA13240, and DA024417). The Women's Genome Health Study (WGHS) was funded by the National Heart, Lung, and Blood Institute (HL043851 and HL080467) and the National Cancer Institute (CA047988 and UM1CA182913), with support for genotyping provided by Amgen. Data collection in the Finnish Twin Registry has been supported by the Wellcome Trust Sanger Institute, the Broad Institute, ENGAGE-European Network for Genetic and Genomic Epidemiology, FP7-HEALTH-F4-2007, grant agreement number 201413, National Institute of Alcohol Abuse and Alcoholism (grants AA-12502, AA-00145, AA-09203, AA15416, and K02AA018755) and the Academy of Finland (grants 100499, 205585, 118555, 141054, 264146, 308248, 312073 and 336823 to J. Kaprio). TwinsUK is funded by the Wellcome Trust, Medical Research Council, Versus Arthritis, European Union Horizon 2020, Chronic Disease Research Foundation (CDRF), Zoe Ltd and the National Institute for Health Research (NIHR) Clinical Research Network (CRN) and Biomedical Research Centre based at Guy's and St Thomas' NHS Foundation Trust in partnership with King's College London. For NESDA, funding was obtained from the Netherlands Organization for Scientific Research (Geestkracht program grant 10000-1002), the Center for Medical Systems Biology (CSMB, NVVO Genomics), Biobanking and Biomolecular Resources Research Infrastructure (BBMRI-NL), VU University's Institutes for Health and Care Research (EMGO+) and Neuroscience Campus Amsterdam, University Medical Center Groningen, Leiden University Medical Center, National Institutes of Health (NIH, ROI D0042157-01A, MH081802, Grand Opportunity grants 1 RC2 Ml-1089951 and IRC2 MH089995). Part of the genotyping and analyses were funded by the Genetic Association Information Network (GAIN) of the Foundation for the National Institutes of Health. Computing was supported by BiG Grid, the Dutch e-Science Grid, which is financially supported by NWO. Work in the Del Bene lab was supported by the Programme Investissements d'Avenir IHU FOReSIGHT (ANR-18-IAHU-01). C.R. was supported by an EU Horizon 2020 Marie Sklodowska-Curie Action fellowship (H2020-MSCA-IF-2014 #661527). H.S. and K.S. are employees of deCODE Genetics/Amgen. The other authors declare no competing financial interests. TRIAL REGISTRATION NUMBER: N/A.

Identification of drug repurposing candidates for the treatment of anxiety: A genetic approach.

Anxiety disorders are a group of prevalent and heritable neuropsychiatric diseases. We previously conducted a genome-wide association study (GWAS) which identified genomic loci associated with anxiety; however, the biological consequences underlying the genetic associations are largely unknown. Integrating GWAS and functional genomic data may improve our understanding of the genetic effects on intermediate molecular phenotypes such as gene expression. This can provide an opportunity for the discovery of drug targets for anxiety via drug repurposing. We used the GWAS summary statistics to determine putative causal genes for anxiety using MAGMA and colocalization analyses. A transcriptome-wide association study was conducted to identify genes with differential genetically regulated levels of gene expression in human brain tissue. The genes were integrated with a large drug-gene expression database (Connectivity Map), discovering compounds that are predicted to "normalise" anxiety-associated expression changes. The study identified 64 putative causal genes associated with anxiety (35 genes upregulated; 29 genes downregulated). Drug mechanisms adrenergic receptor agonists, sigma receptor agonists, and glutamate receptor agonists gene targets were enriched in anxiety-associated genetic signal and exhibited an opposing effect on the anxiety-associated gene expression signature. The significance of the project demonstrated genetic links for novel drug candidates to potentially advance anxiety therapeutics.

The association between trauma exposure, polygenic risk and individual depression symptoms.

BACKGROUND: Traumatic experiences are associated with increased risk for major depressive disorder (MDD). This study sought to determine the extent that trauma exposure, depression polygenic risk scores (PRS), and their interaction are associated with MDD and individual depression symptoms. METHODS: Data from 102,182 individuals from the large-scale UK Biobank population cohort was analysed. A series of regression analyses were conducted to estimate the association between trauma, depression PRS and 1) current depression, 2) lifetime MDD case-control status, 3) nine individual current depressive symptoms, and 4) thirteen individual symptoms experienced during a major depressive episode. Additive and multiplicative PRS-by-trauma interactions were also assessed. RESULTS: Trauma and depression PRS were significantly associated with both current depression and lifetime MDD. A positive, additive interaction effect was observed on depression, but multiplicative interactions were not significant. Trauma exposure and depression PRS were associated with specific patterns of depression symptoms; Trauma was associated with low self-esteem, suicidal ideation, and atypical (but not typical) neurovegetative symptoms. Additive interaction effects were observed on six out of nine current depressive symptoms. CONCLUSIONS: Trauma exposure and genetic predisposition to depression may lead to particular symptomatology, which may contribute to the extreme clinical heterogeneity observed in individuals with major depression.

Ten challenges for clinical translation in psychiatric genetics.

Genome-wide association studies have identified hundreds of robust genetic associations underlying psychiatric disorders and provided important biological insights into disease onset and progression. There is optimism that genetic findings will pave the way to precision psychiatry by facilitating the development of more effective treatments and the identification of groups of patients that these treatments should be targeted toward. However, there are several challenges that must be addressed before genetic findings can be translated into the clinic. In this Perspective, we highlight ten challenges for the field of psychiatric genetics, focused on the robust and generalizable detection of genetic risk factors, improved definition and assessment of psychopathology and achieving better clinical indicators. We discuss recent advancements in the field that will improve the explanatory and predictive power of genetic data and ultimately contribute to improving the management and treatment of patients with a psychiatric disorder.

Genetic evidence that the causal association of educational attainment with reduced risk of Alzheimer's disease is driven by intelligence.

Alzheimer's disease (AD) is predicted to affect 132 million people by 2050. Targeting modifiable lifestyle risk factors that are associated with an increased risk of AD could prevent a large proportion of dementia cases, allowing people to reach the end of their life dementia free. However, evidence obtained from the observational studies does not take into account how risk factors are correlated with one another, and whether they causally contribute to increased AD risk. In this study, we determine whether the relationship between previously speculated AD risk factors and AD susceptibility is consistent with causality using large-scale genetic data. We focus on educational attainment (EA), intelligence and household income which have been previously shown to be causally associated with AD. Using GWAS-by-subtraction and Multivariable Mendelian Randomization we show that of these, only the cognitive component of EA (intelligence) is independently causally associated with AD. This work has ramifications for the modifiability of lifestyle risk factors for AD.

Elucidating the relationship between migraine risk and brain structure using genetic data.

Migraine is a highly common and debilitating disorder that often affects individuals in their most productive years of life. Previous studies have identified both genetic variants and brain morphometry differences associated with migraine risk. However, the relationship between migraine and brain morphometry has not been examined on a genetic level, and the causal nature of the association between brain structure and migraine risk has not been determined. Using the largest available genome-wide association studies to date, we examined the genome-wide genetic overlap between migraine and intracranial volume, as well as the regional volumes of nine subcortical brain structures. We further focused the identification and biological annotation of genetic overlap between migraine and each brain structure on specific regions of the genome shared between migraine and brain structure. Finally, we examined whether the size of any of the examined brain regions causally increased migraine risk using a Mendelian randomization approach. We observed a significant genome-wide negative genetic correlation between migraine risk and intracranial volume (rG = -0.11, P = 1 × 10-3) but not with any subcortical region. However, we identified jointly associated regional genomic overlap between migraine and every brain structure. Gene enrichment in these shared genomic regions pointed to possible links with neuronal signalling and vascular regulation. Finally, we provide evidence of a possible causal relationship between smaller total brain, hippocampal and ventral diencephalon volume and increased migraine risk, as well as a causal relationship between increased risk of migraine and a larger volume of the amygdala. We leveraged the power of large genome-wide association studies to show evidence of shared genetic pathways that jointly influence migraine risk and several brain structures, suggesting that altered brain morphometry in individuals with high migraine risk may be genetically mediated. Further interrogation of these results showed support for the neurovascular hypothesis of migraine aetiology and shed light on potentially viable therapeutic targets.

A Local Genetic Correlation Analysis Provides Biological Insights Into the Shared Genetic Architecture of Psychiatric and Substance Use Phenotypes.

BACKGROUND: Global genetic correlation analysis has provided valuable insight into the shared genetic basis between psychiatric and substance use disorders. However, little is known about which regions disproportionately contribute to the global correlation. METHODS: We used Local Analysis of [co]Variant Annotation to calculate bivariate local genetic correlations across 2495 approximately equal-sized, semi-independent genomic regions for 20 psychiatric and substance use phenotypes. We performed a transcriptome-wide association study using expression weights from the prefrontal cortex to identify risk genes for each phenotype, followed by probabilistic fine-mapping to prioritize credible causal genes within each bivariate locus. RESULTS: We detected 80 significant (p < 2.08 × 10-6) bivariate local genetic correlations across 61 loci. The expression effect directions for risk genes within each bivariate locus were largely consistent with the local correlation coefficients, suggesting that genetically regulated gene expression may be used in the functional interpretation of local genetic correlations. Probabilistic fine-mapping identified several genes that may drive pleiotropic mechanisms for genetically correlated phenotypes. For example, we confirmed a local genetic correlation between schizophrenia and smoking behavior at 15q25 and prioritized PSMA4 as the most credible gene candidate underlying both phenotypes. CONCLUSIONS: Our study reveals previously unreported local bivariate genetic correlations between psychiatric and substance use phenotypes, which we fine-mapped to identify shared credible causal genes underlying genetically correlated phenotypes.

An analysis of genetically regulated gene expression and the role of co-expression networks across 16 psychiatric and substance use phenotypes.

Genome-wide association studies (GWASs) have identified thousands of risk loci for psychiatric and substance use phenotypes, however the biological consequences of these loci remain largely unknown. We performed a transcriptome-wide association study of 10 psychiatric disorders and 6 substance use phenotypes (GWAS sample size range, N = 9725-807,553) using expression quantitative trait loci data from 532 prefrontal cortex samples. We estimated the correlation of genetically regulated expression between phenotype pairs, and compared the results with the genetic correlations. We identified 393 genes with at least one significant phenotype association, comprising 458 significant associations across 16 phenotypes. Overall, the transcriptomic correlations for phenotype pairs were significantly higher than the respective genetic correlations. For example, attention deficit hyperactivity disorder and autism spectrum disorder, both childhood developmental disorders, had significantly higher transcriptomic correlation (r = 0.84) than genetic correlation (r = 0.35). Finally, we tested the enrichment of phenotype-associated genes in gene co-expression networks built from human prefrontal cortex samples. Phenotype-associated genes were enriched in multiple gene co-expression modules and the implicated modules contained genes involved in mRNA splicing and glutamatergic receptors, among others. Together, our results highlight the utility of gene expression data in the understanding of functional gene mechanisms underlying psychiatric disorders and substance use phenotypes.

Integrative Network-Based Analysis Reveals Gene Networks and Novel Drug Repositioning Candidates for Alzheimer Disease.

BACKGROUND AND OBJECTIVES: To integrate genome-wide association study data with tissue-specific gene expression information to identify coexpression networks, biological pathways, and drug repositioning candidates for Alzheimer disease. METHODS: We integrated genome-wide association summary statistics for Alzheimer disease with tissue-specific gene coexpression networks from brain tissue samples in the Genotype-Tissue Expression study. We identified gene coexpression networks enriched with genetic signals for Alzheimer disease and characterized the associated networks using biological pathway analysis. The disease-implicated modules were subsequently used as a molecular substrate for a computational drug repositioning analysis, in which we (1) imputed genetically regulated gene expression within Alzheimer disease implicated modules; (2) integrated the imputed gene expression levels with drug-gene signatures from the connectivity map to identify compounds that normalize dysregulated gene expression underlying Alzheimer disease; and (3) prioritized drug compounds and mechanisms of action based on the extent to which they normalize dysregulated expression signatures. RESULTS: Genetic factors for Alzheimer disease are enriched in brain gene coexpression networks involved in the immune response. Computational drug repositioning analyses of expression changes within the disease-associated networks retrieved known Alzheimer disease drugs (e.g., memantine) as well as biologically meaningful drug categories (e.g., glutamate receptor antagonists). DISCUSSION: Our results improve the biological interpretation of genetic data for Alzheimer disease and provide a list of potential antidementia drug repositioning candidates for which the efficacy should be investigated in functional validation studies.

Symptom-level modelling unravels the shared genetic architecture of anxiety and depression.

Depression and anxiety are highly prevalent and comorbid psychiatric traits that cause considerable burden worldwide. Here we use factor analysis and genomic structural equation modelling to investigate the genetic factor structure underlying 28 items assessing depression, anxiety and neuroticism, a closely related personality trait. Symptoms of depression and anxiety loaded on two distinct, although highly genetically correlated factors, and neuroticism items were partitioned between them. We used this factor structure to conduct genome-wide association analyses on latent factors of depressive symptoms (89 independent variants, 61 genomic loci) and anxiety symptoms (102 variants, 73 loci) in the UK Biobank. Of these associated variants, 72% and 78%, respectively, replicated in an independent cohort of approximately 1.9 million individuals with self-reported diagnosis of depression and anxiety. We use these results to characterize shared and trait-specific genetic associations. Our findings provide insight into the genetic architecture of depression and anxiety and comorbidity between them.

E-MAGMA: an eQTL-informed method to identify risk genes using genome-wide association study summary statistics.

MOTIVATION: Genome-wide association studies have successfully identified multiple independent genetic loci that harbour variants associated with human traits and diseases, but the exact causal genes are largely unknown. Common genetic risk variants are enriched in non-protein-coding regions of the genome and often affect gene expression (expression quantitative trait loci, eQTL) in a tissue-specific manner. To address this challenge, we developed a methodological framework, E-MAGMA, which converts genome-wide association summary statistics into gene-level statistics by assigning risk variants to their putative genes based on tissue-specific eQTL information. RESULTS: We compared E-MAGMA to three eQTL informed gene-based approaches using simulated phenotype data. Phenotypes were simulated based on eQTL reference data using GCTA for all genes with at least one eQTL at chromosome 1. We performed 10 simulations per gene. The eQTL-h2 (i.e. the proportion of variation explained by the eQTLs) was set at 1%, 2% and 5%. We found E-MAGMA outperforms other gene-based approaches across a range of simulated parameters (e.g. the number of identified causal genes). When applied to genome-wide association summary statistics for five neuropsychiatric disorders, E-MAGMA identified more putative candidate causal genes compared to other eQTL-based approaches. By integrating tissue-specific eQTL information, these results show E-MAGMA will help to identify novel candidate causal genes from genome-wide association summary statistics and thereby improve the understanding of the biological basis of complex disorders. AVAILABILITY AND IMPLEMENTATION: A tutorial and input files are made available in a github repository: https://github.com/eskederks/eMAGMA-tutorial. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

An integrative systems-based analysis of substance use: eQTL-informed gene-based tests, gene networks, and biological mechanisms.

Genome-wide association studies have identified multiple genetic risk factors underlying susceptibility to substance use, however, the functional genes and biological mechanisms remain poorly understood. The discovery and characterization of risk genes can be facilitated by the integration of genome-wide association data and gene expression data across biologically relevant tissues and/or cell types to identify genes whose expression is altered by DNA sequence variation (expression quantitative trait loci; eQTLs). The integration of gene expression data can be extended to the study of genetic co-expression, under the biologically valid assumption that genes form co-expression networks to influence the manifestation of a disease or trait. Here, we integrate genome-wide association data with gene expression data from 13 brain tissues to identify candidate risk genes for 8 substance use phenotypes. We then test for the enrichment of candidate risk genes within tissue-specific gene co-expression networks to identify modules (or groups) of functionally related genes whose dysregulation is associated with variation in substance use. We identified eight gene modules in brain that were enriched with gene-based association signals for substance use phenotypes. For example, a single module of 40 co-expressed genes was enriched with gene-based associations for drinks per week and biological pathways involved in GABA synthesis, release, reuptake and degradation. Our study demonstrates the utility of eQTL and gene co-expression analysis to uncover novel biological mechanisms for substance use traits.

An analysis of genetically regulated gene expression across multiple tissues implicates novel gene candidates in Alzheimer's disease.

INTRODUCTION: Genome-wide association studies (GWAS) have successfully identified multiple independent genetic loci that harbour variants associated with Alzheimer's disease, but the exact causal genes and biological pathways are largely unknown. METHODS: To prioritise likely causal genes associated with Alzheimer's disease, we used S-PrediXcan to integrate expression quantitative trait loci (eQTL) from the Genotype-Tissue Expression (GTEx) study and CommonMind Consortium (CMC) with Alzheimer's disease GWAS summary statistics. We meta-analysed the GTEx results using S-MultiXcan, prioritised disease-implicated loci using a computational fine-mapping approach, and performed a biological pathway analysis on the gene-based results. RESULTS: We identified 126 tissue-specific gene-based associations across 48 GTEx tissues, targeting 50 unique genes. Meta-analysis of the tissue-specific associations identified 73 genes whose expression was associated with Alzheimer's disease. Additional analyses in the dorsolateral prefrontal cortex from the CMC identified 12 significant associations, 8 of which also had a significant association in GTEx tissues. Fine-mapping of causal gene sets prioritised gene candidates in 10 Alzheimer's disease loci with strong evidence for causality. Biological pathway analyses of the meta-analysed GTEx data and CMC data identified a significant enrichment of Alzheimer's disease association signals in plasma lipoprotein clearance, in addition to multiple immune-related pathways. CONCLUSIONS: Gene expression data from brain and peripheral tissues can improve power to detect regulatory variation underlying Alzheimer's disease. However, the associations in peripheral tissues may reflect tissue-shared regulatory variation for a gene. Therefore, future functional studies should be performed to validate the biological meaning of these associations and whether they represent new pathogenic tissues.