Hostname: page-component-7b9c58cd5d-hxdxx Total loading time: 0 Render date: 2025-03-16T09:23:39.413Z Has data issue: false hasContentIssue false

Translating GWAS findings into therapies for depression and anxiety disorders: gene-set analyses reveal enrichment of psychiatric drug classes and implications for drug repositioning

Published online by Cambridge University Press:  20 December 2018

Hon-Cheong So*
Affiliation:
School of Biomedical Sciences, The Chinese University of Hong Kong, Shatin, Hong Kong KIZ-CUHK Joint Laboratory of Bioresources and Molecular Research of Common Diseases, Kunming Zoology Institute of Zoology and The Chinese University of Hong Kong, Shatin, Hong Kong
Carlos Kwan-Long Chau
Affiliation:
School of Biomedical Sciences, The Chinese University of Hong Kong, Shatin, Hong Kong
Alexandria Lau
Affiliation:
School of Biomedical Sciences, The Chinese University of Hong Kong, Shatin, Hong Kong
Sze-Yung Wong
Affiliation:
Faculty of Medicine, The Chinese University of Hong Kong, Shatin, Hong Kong
Kai Zhao
Affiliation:
School of Biomedical Sciences, The Chinese University of Hong Kong, Shatin, Hong Kong
*
Author for correspondence: Hon-Cheong So, E-mail: hcso@cuhk.edu.hk
Rights & Permissions [Opens in a new window]

Abstract

Background

Depression and anxiety disorders (AD) are the first and sixth leading causes of disability worldwide. Despite their high prevalence and significant disability resulted, there are limited advances in new drug development. Recently, genome-wide association studies (GWAS) have greatly advanced our understanding of the genetic basis underlying psychiatric disorders.

Methods

Here we employed gene-set analyses of GWAS summary statistics for drug repositioning. We explored five related GWAS datasets, including two on major depressive disorder (MDD2018 and MDD-CONVERGE, with the latter focusing on severe melancholic depression), one on AD, and two on depressive symptoms and neuroticism in the population. We extracted gene-sets associated with each drug from DSigDB and examined their association with each GWAS phenotype. We also performed repositioning analyses on meta-analyzed GWAS data, integrating evidence from all related phenotypes.

Results

Importantly, we showed that the repositioning hits are generally enriched for known psychiatric medications or those considered in clinical trials. Enrichment was seen for antidepressants and anxiolytics but also for antipsychotics. We also revealed new candidates or drug classes for repositioning, some of which were supported by experimental or clinical studies. For example, the top repositioning hit using meta-analyzed p values was fendiline, which was shown to produce antidepressant-like effects in mouse models by inhibition of acid sphingomyelinase.

Conclusion

Taken together, our findings suggest that human genomic data such as GWAS are useful in guiding drug discoveries for depression and AD.

Type
Original Articles
Copyright
Copyright © Cambridge University Press 2018 

Introduction

Depression and anxiety disorders (AD) are among the most common psychiatric disorders. According to the latest report by the World Health Organization, depression affects more than 300 million people worldwide and is the leading cause of disability (World Health Organization, 2017). AD affect more than 260 million people and is the sixth leading cause of disability (World Health Organization, 2017). The two disorders are highly comorbid and might share common pathophysiologies (Frances et al., Reference Frances, Manning, Marin, Kocsis, Mckinney, Hall and Kline1992; Goodwin, Reference Goodwin2015). Nevertheless, pharmacological treatment for major depressive disorder (MDD) or AD has not seen much advance in the last two decades or so, with a lack of therapies having novel mechanisms of action. In addition, only about one-third of MDD patients achieve complete remission after a single antidepressant trial (Trivedi et al., Reference Trivedi, Rush, Wisniewski, Nierenberg, Warden, Ritz, Norquist, Howland, Lebowitz, Mcgrath, Shores-Wilson, Biggs, Balasubramani, Fava and Team2006) and around 10–30% of patients are treatment-resistant (Al-Harbi, Reference Al-Harbi2012).

On the other hand, with the advent of high-throughput technologies such as genome-wide association studies (GWAS) in the last decade, we have gained a much better understanding of the genetic bases of many complex diseases. It is hoped that human genomics data will accelerate drug development for psychiatric disorders, especially due to the difficulties for animal models to fully mimic human psychiatric conditions such as depression (Papassotiropoulos and de Quervain, Reference Papassotiropoulos and De Quervain2015).

We hypothesize that gene-sets associated with antidepressants or anxiolytics, or more generally with psychiatric medications, will be enriched among the GWAS results of depression and anxiety phenotypes. If the gene-set analysis (GSA) approach is able to ‘re-discover’ known treatments, it might also be able to reveal new therapies for depression and AD. Practically, if the pathway targeted by a particular drug (or chemical) is significantly enriched in GWAS but the drug is not a known treatment for the disease, it may serve as a good candidate for repositioning.

With regards to GWAS-based drug discovery, the current focus is mainly on identification of new drug targets from the top GWAS hits. In an earlier study, Sanseau et al. (Reference Sanseau, Agarwal, Barnes, Pastinen, Richards, Cardon and Mooser2012) identified the most significant GWAS hits from a range of diseases and compared them against known drug targets to find ‘mismatches’ (i.e. drug indication different from the studied disorder) as candidates for repurposing. While it is intuitive and legitimate to focus on the most significant SNPs, for many complex traits, the genetic architecture may be highly polygenic and variants of weaker effects may be ‘hidden’. Moreover, given the complex and multifactorial etiologies of many complex diseases, the development of multi-target drugs with wide-ranging biological activities (known as ‘polypharmacology’) is gaining increased attention [please refer to, e.g. Anighoro et al. (Reference Anighoro, Bajorath and Rastelli2014) for a review]. It is argued that multi-target drugs may have improved efficacy over highly selective pharmacological agents, as they tackle multiple pathogenic pathways in the system. A gene-set or pathway-based approach to drug repositioning follows this multi-target paradigm.

GSA is an established approach to gain biological insight into expression microarrays, GWAS, or other high-throughput ‘omics’ studies (Maciejewski, Reference Maciejewski2014; de Leeuw et al., Reference De Leeuw, Neale, Heskes and Posthuma2016). De Jong et al. made use of GSA to identify repurposing opportunities for schizophrenia (de Jong et al., Reference De Jong, Vidler, Mokrab, Collier and Breen2016) and several top results were supported by the literature. In another very recent study, Gaspar and Breen (Reference Gaspar and Breen2017) performed further analyses of GSA results, and reported that GWAS signals of schizophrenia are enriched for neuropsychiatric medications as sample size increases. Another related study on schizophrenia was conducted by Ruderfer et al. (Reference Ruderfer, Charney, Readhead, Kidd, Kahler, Kenny, Keiser, Moran, Hultman, Scott, Sullivan, Purcell, Dudley and Sklar2016), who collected genome-wide significant GWAS variants and exome sequencing results and compared the identified genes against known or predicted drug targets. Significant enrichment was observed for antipsychotics.

In this study, we take a different focus on depression and AD, which are highly prevalent and disabling disorders. Besides considering individual GWAS studies of phenotypes related to depression and anxiety, we also performed analyses on combined GWAS summary statistics to improve the power.

Methods

Overview of analytic approach

Our analyses can be broadly divided into two steps. Firstly, we extracted gene-sets associated with a variety of drugs (or chemicals), and tested whether the gene-set associated with each individual drug is enriched among the GWAS results. The drugs ranked among the top (i.e. those with lower p values) were considered potential candidates for repositioning. In the second step, we performed analyses on the drugs. We evaluated the prioritized drugs and tested which drug classes were enriched. As we have hypothesized above, we would specifically test whether antidepressants/anxiolytics and other psychiatric medications were enriched among the repositioning results. Figure 1 describes an overview of our analytic approach.

Fig. 1. A schematic diagram showing our analysis workflow. We first extracted gene-sets associated with a variety of drugs (and chemicals), and examined whether the gene-set associated with each drug is enriched among the GWAS results. The drugs ranked among the top were considered potential candidates for repositioning, and a systemic literature search was performed for the top 20 repositioning hits of each psychiatric trait and meta-analyzed results. We also evaluated which drug classes were enriched. We specifically tested whether antidepressants/anxiolytics and other psychiatric medications were enriched, and then also examined all ATC (level 3) drug classes for enrichment in a hypothesis-free manner (ATC, Anatomical Therapeutic Chemical Classification System).

GWAS data

In the main analysis, we considered five GWAS datasets that are associated with depression and anxiety. Two are studies of MDD, namely MDD-2018 (Wray and Sullivan, Reference Wray, Ripke, Mattheisen, Trzaskowski, Byrne, Abdellaoui, Adams, Agerbo, Air, Andlauer, Bacanu, Baekvad-Hansen, Beekman, Bigdeli, Binder, Blackwood, Bryois, Buttenschon, Bybjerg-Grauholm, Cai, Castelao, Christensen, Clarke, Coleman, Colodro-Conde, Couvy-Duchesne, Craddock, Crawford, Crowley, Dashti, Davies, Deary, Degenhardt, Derks, Direk, Dolan, Dunn, Eley, Eriksson, Escott-Price, Kiadeh, Finucane, Forstner, Frank, Gaspar, Gill, Giusti-Rodriguez, Goes, Gordon, Grove, Hall, Hannon, Hansen, Hansen, Herms, Hickie, Hoffmann, Homuth, Horn, Hottenga, Hougaard, Hu, Hyde, Ising, Jansen, Jin, Jorgenson, Knowles, Kohane, Kraft, Kretzschmar, Krogh, Kutalik, Lane, Li, Li, Lind, Liu, Lu, MacIntyre, MacKinnon, Maier, Maier, Marchini, Mbarek, McGrath, McGuffin, Medland, Mehta, Middeldorp, Mihailov, Milaneschi, Milani, Mill, Mondimore, Montgomery, Mostafavi, Mullins, Nauck, Ng, Nivard, Nyholt, O'Reilly, Oskarsson, Owen, Painter, Pedersen, Pedersen, Peterson, Pettersson, Peyrot, Pistis, Posthuma, Purcell, Quiroz, Qvist, Rice, Riley, Rivera, Saeed Mirza, Saxena, Schoevers, Schulte, Shen, Shi, Shyn, Sigurdsson, Sinnamon, Smit, Smith, Stefansson, Steinberg, Stockmeier, Streit, Strohmaier, Tansey, Teismann, Teumer, Thompson, Thomson, Thorgeirsson, Tian, Traylor, Treutlein, Trubetskoy, Uitterlinden, Umbricht, Van der Auwera, van Hemert, Viktorin, Visscher, Wang, Webb, Weinsheimer, Wellmann, Willemsen, Witt, Wu, Xi, Yang, Zhang, Arolt, Baune, Berger, Boomsma, Cichon, Dannlowski, de Geus, DePaulo, Domenici, Domschke, Esko, Grabe, Hamilton, Hayward, Heath, Hinds, Kendler, Kloiber, Lewis, Li, Lucae, Madden, Magnusson, Martin, McIntosh, Metspalu, Mors, Mortensen, Muller-Myhsok, Nordentoft, Nothen, O'Donovan, Paciga, Pedersen, Penninx, Perlis, Porteous, Potash, Preisig, Rietschel, Schaefer, Schulze, Smoller, Stefansson, Tiemeier, Uher, Volzke, Weissman, Werge, Winslow, Lewis, Levinson, Breen, Borglum and Sullivan2017) and MDD-CONVERGE (Cai et al., Reference Cai, Bigdeli, Kretzschmar, Li, Liang, Song, Hu, Li, Jin, Hu, Wang, Wang, Qian, Liu, Jiang, Lu, Zhang, Yin, Li, Xu, Gao, Reimers, Webb, Riley, Bacanu, Peterson, Chen, Zhong, Liu, Wang, Sun, Sang, Jiang, Zhou, Li, Li, Zhang, Wang, Fang, Pan, Miao, Zhang, Hu, Yu, Du, Sang, Li, Chen, Cai, Yang, Yang, Ha, Hong, Deng, Li, Li, Song, Gao, Zhang, Gan, Meng, Pan, Gao, Zhang, Sun, Li, Niu, Zhang, Liu, Hu, Zhang, Lv, Dong, Wang, Tao, Wang, Xia, Rong, He, Liu, Huang, Mei, Shen, Liu, Shen, Tian, Liu, Wu, Gu, Fu, Shi, Chen, Gan, Liu, Wang, Yang, Cong, Marchini, Yang and Wang2015). However, the samples included in these studies are different. MDD-CONVERGE (N = 10 640) is a cohort of Chinese women, and this sample mainly consists of hospital-ascertained cases affected by severe depression, of whom ~85% had melancholic symptoms (Cai et al., Reference Cai, Bigdeli, Kretzschmar, Li, Liang, Song, Hu, Li, Jin, Hu, Wang, Wang, Qian, Liu, Jiang, Lu, Zhang, Yin, Li, Xu, Gao, Reimers, Webb, Riley, Bacanu, Peterson, Chen, Zhong, Liu, Wang, Sun, Sang, Jiang, Zhou, Li, Li, Zhang, Wang, Fang, Pan, Miao, Zhang, Hu, Yu, Du, Sang, Li, Chen, Cai, Yang, Yang, Ha, Hong, Deng, Li, Li, Song, Gao, Zhang, Gan, Meng, Pan, Gao, Zhang, Sun, Li, Niu, Zhang, Liu, Hu, Zhang, Lv, Dong, Wang, Tao, Wang, Xia, Rong, He, Liu, Huang, Mei, Shen, Liu, Shen, Tian, Liu, Wu, Gu, Fu, Shi, Chen, Gan, Liu, Wang, Yang, Cong, Marchini, Yang and Wang2015). The MDD-2018 (N = 173 005, for which full summary statistics are available) samples are composed of Caucasians of both sexes, and they are more heterogeneous and not specifically enriched for any subtypes of depression (Wray et al., Reference Ripke, Wray, Lewis, Hamilton, Weissman, Breen, Byrne, Blackwood, Boomsma, Cichon, Heath, Holsboer, Lucae, Madden, Martin, McGuffin, Muglia, Noethen, Penninx, Pergadia, Potash, Rietschel, Lin, Müller-Myhsok, Shi, Steinberg, Grabe, Lichtenstein, Magnusson, Perlis, Preisig, Smoller, Stefansson, Uher, Kutalik, Tansey, Teumer, Viktorin, Barnes, Bettecken, Binder, Breuer, Castro, Churchill, Coryell, Craddock, Craig, Czamara, De Geus, Degenhardt, Farmer, Fava, Frank, Gainer, Gallagher, Gordon, Goryachev, Gross, Guipponi, Henders, Herms, Hickie, Hoefels, Hoogendijk, Hottenga, Iosifescu, Ising, Jones, Jones, Jung-Ying, Knowles, Kohane, Kohli, Korszun, Landen, Lawson, Lewis, Macintyre, Maier, Mattheisen, McGrath, McIntosh, McLean, Middeldorp, Middleton, Montgomery, Murphy, Nauck, Nolen, Nyholt, O'Donovan, Oskarsson, Pedersen, Scheftner, Schulz, Schulze, Shyn, Sigurdsson, Slager, Smit, Stefansson, Steffens, Thorgeirsson, Tozzi, Treutlein, Uhr, van den Oord, Van Grootheest, Völzke, Weilburg, Willemsen, Zitman, Neale, Daly, Levinson and Sullivan2018). The MDD-2018 study represents the largest and most recent GWAS meta-analysis on MDD, and includes all subjects from an earlier GWAS, i.e. the MDD-PGC-2012 study by Major Depressive Disorder Working Group of the Psychiatric GWAS Consortium et al. (Reference Ripke, Wray, Lewis, Hamilton, Weissman, Breen, Byrne, Blackwood, Boomsma, Cichon, Heath, Holsboer, Lucae, Madden, Martin, McGuffin, Muglia, Noethen, Penninx, Pergadia, Potash, Rietschel, Lin, Müller-Myhsok, Shi, Steinberg, Grabe, Lichtenstein, Magnusson, Perlis, Preisig, Smoller, Stefansson, Uher, Kutalik, Tansey, Teumer, Viktorin, Barnes, Bettecken, Binder, Breuer, Castro, Churchill, Coryell, Craddock, Craig, Czamara, De Geus, Degenhardt, Farmer, Fava, Frank, Gainer, Gallagher, Gordon, Goryachev, Gross, Guipponi, Henders, Herms, Hickie, Hoefels, Hoogendijk, Hottenga, Iosifescu, Ising, Jones, Jones, Jung-Ying, Knowles, Kohane, Kohli, Korszun, Landen, Lawson, Lewis, Macintyre, Maier, Mattheisen, McGrath, McIntosh, McLean, Middeldorp, Middleton, Montgomery, Murphy, Nauck, Nolen, Nyholt, O'Donovan, Oskarsson, Pedersen, Scheftner, Schulz, Schulze, Shyn, Sigurdsson, Slager, Smit, Stefansson, Steffens, Thorgeirsson, Tozzi, Treutlein, Uhr, van den Oord, Van Grootheest, Völzke, Weilburg, Willemsen, Zitman, Neale, Daly, Levinson and Sullivan2013) (N = 18 759). Note that for MDD-2018, full summary statistics (on which the current analysis is based) are only publicly available for samples excluding 23andMe subjects.

Another two GWAS studies were meta-analyses on depressive symptoms and neuroticism conducted by the Social Science Genetics Association Consortium (SSGAC) (Okbay et al., Reference Okbay, Baselmans, De Neve, Turley, Nivard, Fontana, Meddens, Linner, Rietveld, Derringer, Gratten, Lee, Liu, de Vlaming, Ahluwalia, Buchwald, Cavadino, Frazier-Wood, Furlotte, Garfield, Geisel, Gonzalez, Haitjema, Karlsson, van der Laan, Ladwig, Lahti, van der Lee, Lind, Liu, Matteson, Mihailov, Miller, Minica, Nolte, Mook-Kanamori, van der Most, Oldmeadow, Qian, Raitakari, Rawal, Realo, Rueedi, Schmidt, Smith, Stergiakouli, Tanaka, Taylor, Thorleifsson, Wedenoja, Wellmann, Westra, Willems, Zhao, LifeLines Cohort, Amin, Bakshi, Bergmann, Bjornsdottir, Boyle, Cherney, Cox, Davies, Davis, Ding, Direk, Eibich, Emeny, Fatemifar, Faul, Ferrucci, Forstner, Gieger, Gupta, Harris, Harris, Holliday, Hottenga, De Jager, Kaakinen, Kajantie, Karhunen, Kolcic, Kumari, Launer, Franke, Li-Gao, Liewald, Koini, Loukola, Marques-Vidal, Montgomery, Mosing, Paternoster, Pattie, Petrovic, Pulkki-Raback, Quaye, Raikkonen, Rudan, Scott, Smith, Sutin, Trzaskowski, Vinkhuyzen, Yu, Zabaneh, Attia, Bennett, Berger, Bertram, Boomsma, Snieder, Chang, Cucca, Deary, van Duijn, Eriksson, Bultmann, de Geus, Groenen, Gudnason, Hansen, Hartman, Haworth, Hayward, Heath, Hinds, Hypponen, Iacono, Jarvelin, Jockel, Kaprio, Kardia, Keltikangas-Jarvinen, Kraft, Kubzansky, Lehtimaki, Magnusson, Martin, McGue, Metspalu, Mills, de Mutsert, Oldehinkel, Pasterkamp, Pedersen, Plomin, Polasek, Power, Rich, Rosendaal, den Ruijter, Schlessinger, Schmidt, Svento, Schmidt, Alizadeh, Sorensen, Spector, Starr, Stefansson, Steptoe, Terracciano, Thorsteinsdottir, Thurik, Timpson, Tiemeier, Uitterlinden, Vollenweider, Wagner, Weir, Yang, Conley, Smith, Hofman, Johannesson, Laibson, Medland, Meyer, Pickrell, Esko, Krueger, Beauchamp, Koellinger, Benjamin, Bartels and Cesarini2016). The meta-analysis on depressive symptoms (SSGAC-DS) included the MDD-PGC-2012 study (N = 18 759) and a case–control sample from the Genetic Epidemiology Research on Aging (GERA) Cohort (N = 56 368), and it also comprised a UK BioBank sample made up of general population (N = 105 739). Depressive symptoms were measured by a self-reported questionnaire. We also included another study on neuroticism (SSGAC-NEU) (N = 170 906), as this personality trait is known to be closely associated with depression and AD (Lahey, Reference Lahey2009). In addition, antidepressants may affect personality traits, including a reduction in neuroticism, independent of their effects on depressive symptoms (Tang et al., Reference Tang, Derubeis, Hollon, Amsterdam, Shelton and Schalet2009). The aim of including SSGAC-DS and SSGAC-NEU is to investigate whether the study of depression/anxiety-related phenotypes in the general population, other than clinical depression or AD, may affect the potential of guiding drug discovery. An analogous example in cardiovascular medicine is that some genetic variants associated with raised fasting glucose may not be linked to diabetes (Florez et al., Reference Florez, Jablonski, Mcateer, Franks, Mason, Mather, Horton, Goldberg, Dabelea, Kahn, Arakaki, Shuldiner and Knowler2012; Merino et al., Reference Merino, Leong, Posner, Porneala, Masana, Dupuis and Florez2017); hence the results from our gene-set and drug repositioning analysis might also differ when different phenotypes are studied.

The last dataset is a GWAS meta-analysis of AD, including generalized AD, panic disorder, social phobia, agoraphobia, and specific phobias (Otowa et al., Reference Otowa, Hek, Lee, Byrne, Mirza, Nivard, Bigdeli, Aggen, Adkins, Wolen, Fanous, Keller, Castelao, Kutalik, Der Auwera, Homuth, Nauck, Teumer, Milaneschi, Hottenga, Direk, Hofman, Uitterlinden, Mulder, Henders, Medland, Gordon, Heath, Madden, Pergadia, Van Der Most, Nolte, Van Oort, Hartman, Oldehinkel, Preisig, Grabe, Middeldorp, Penninx, Boomsma, Martin, Montgomery, Maher, Van Den Oord, Wray, Tiemeier and Hettema2016b). We extracted the GWAS results of the case–control analyses (N = 17 310).

GWAS summary results were downloaded from https://www.med.unc.edu/pgc/results-and-downloads and https://www.thessgac.org/data.

Extracting gene-sets associated with each drug

We made use of the DSigDB database (Yoo et al., Reference Yoo, Shin, Kim, Ryall, Lee, Lee, Jeon, Kang and Tan2015) to extract gene-sets related to each drug. DSigDB holds gene-sets for a total of 17 839 unique compounds. The gene-sets were compiled according to multiple sources: (1) bioassay results from PubChem (Kim et al., Reference Kim, Thiessen, Bolton, Chen, Fu, Gindulyte, Han, He, He, Shoemaker, Wang, Yu, Zhang and Bryant2016) and ChEMBL (Gaulton et al., Reference Gaulton, Bellis, Bento, Chambers, Davies, Hersey, Light, Mcglinchey, Michalovich, Al-Lazikani and Overington2012); (2) kinase profiling assay from the literature and two kinase databases (Medical Research Council Kinase Inhibitor database and Harvard Medical School Library of Integrated Network-based Cellular Signatures database); (3) differentially expressed genes after drug treatment (with >2 fold-change compared with controls), as derived from the Connectivity Map (Lamb et al., Reference Lamb, Crawford, Peck, Modell, Blat, Wrobel, Lerner, Brunet, Subramanian, Ross, Reich, Hieronymus, Wei, Armstrong, Haggarty, Clemons, Wei, Carr, Lander and Golub2006); and (4) manually curated and text mined drug targets from the Therapeutics Targets Database (Li et al., Reference Li, Yu, Li, Zhang, Tang, Yang, Fu, Zhang, Cui, Tu, Zhang, Li, Yang, Sun, Qin, Zeng, Chen, Chen and Zhu2018) and the Comparative Toxicogenomics Database (Mattingly et al., Reference Mattingly, Colby, Forrest and Boyer2003). We downloaded the entire database from http://tanlab.ucdenver.edu/DSigDB. The above sources captured different aspects of drug properties, for example, differentially expressed genes from perturbation experiments might be different from drug target genes defined in bioassay studies. Besides performing analyses based on the whole database which incorporates a broad definition of drug-related genes, we also performed a separate enrichment analysis for genes derived from PubChem and ChEMBL only, as they represent conventionally defined ‘drug targets’ that are more well-studied and perhaps more directly associated with drug actions.

It should be noted that although the focus is on drug ‘repositioning’, the analytic framework is general and can apply to any drugs with some known associated genes. Indeed DSigDB contains a substantial number of drugs which do not have an approved indication yet, which were still included in our analyses.

GSA approach

We first converted the SNP-based test results to gene-based test results. We employed fastBAT (Bakshi et al., Reference Bakshi, Zhu, Vinkhuyzen, Hill, Mcrae, Visscher and Yang2016) (included in the software package GCTA) for gene-based analyses. FastBAT computes the sum of χ2 statistics over all SNPs within a gene and uses an analytic approach to compute the p value. Gene size and linkage disequilibrium patterns are taken into account when computing the p values. The same statistical approach for gene-based tests is also used by two other popular programs, VEGAS (Liu et al., Reference Liu, Mcrae, Nyholt, Medland, Wray, Brown, Hayward, Montgomery, Visscher, Martin, Macgregor and Investigators2010) and PLINK (Purcell et al., Reference Purcell, Neale, Todd-Brown, Thomas, Ferreira, Bender, Maller, Sklar, De Bakker, Daly and Sham2007), although they computed p values by simulations or permutations. FastBAT has been shown to be equivalent to VEGAS and PLINK at higher p values (>1E-06) and more accurate than these two programs for smaller p (Bakshi et al., Reference Bakshi, Zhu, Vinkhuyzen, Hill, Mcrae, Visscher and Yang2016). We ran fastBAT with the default settings and used the 1000 Genomes genotype data as the reference panel.

We then performed a standard GSA by comparing gene-based test statistics within and outside the gene-set. We adopted the same approach as implemented in MAGMA (de Leeuw et al., Reference De Leeuw, Mooij, Heskes and Posthuma2015), which is also reviewed in de Leeuw et al. (Reference De Leeuw, Neale, Heskes and Posthuma2016). Briefly, gene-based p values are first converted to z-statistics by z = Φ−1(p), where Φ−1 is the probit function (more negative z-values represent stronger statistical associations). We then employed a single-sided two-sample t test to see if the mean z-statistics of genes within the gene-set is lower than that outside the gene-set. To avoid results driven by only a few genes, we only considered drugs with at least five genes in their gene-sets. A total of 5232 drugs were included for final analyses.

Combining p values across datasets

Besides analyzing each GWAS dataset in turn for repositioning opportunities, we also considered the aggregate contribution of all datasets, as depression, anxiety, and neuroticism are closely connected to each other. Depression and AD are highly clinically comorbid (Lamers et al., Reference Lamers, Van Oppen, Comijs, Smit, Spinhoven, Van Balkom, Nolen, Zitman, Beekman and Penninx2011; Kessler et al., Reference Kessler, Sampson, Berglund, Gruber, Al-Hamzawi, Andrade, Bunting, Demyttenaere, Florescu, De Girolamo, Gureje, He, Hu, Huang, Karam, Kovess-Masfety, Lee, Levinson, Medina Mora, Moskalewicz, Nakamura, Navarro-Mateu, Browne, Piazza, Posada-Villa, Slade, Ten Have, Torres, Vilagut, Xavier, Zarkov, Shahly and Wilcox2015), demonstrate significant genetic correlations (Otowa et al., Reference Otowa, Hek, Lee, Byrne, Mirza, Nivard, Bigdeli, Aggen, Adkins, Wolen, Fanous, Keller, Castelao, Kutalik, Der Auwera, Homuth, Nauck, Teumer, Milaneschi, Hottenga, Direk, Hofman, Uitterlinden, Mulder, Henders, Medland, Gordon, Heath, Madden, Pergadia, Van Der Most, Nolte, Van Oort, Hartman, Oldehinkel, Preisig, Grabe, Middeldorp, Penninx, Boomsma, Martin, Montgomery, Maher, Van Den Oord, Wray, Tiemeier and Hettema2016a), and share similar pharmacological treatments (Ballenger, Reference Ballenger2000). As mentioned above, neuroticism is associated with both depression and anxiety, and levels of neuroticism may be affected by antidepressant treatment. Studies also suggest a shared genetic basis between neuroticism and depression or AD (Jardine et al., Reference Jardine, Martin, Henderson and Rao1984; Gale et al., Reference Gale, Hagenaars, Davies, Hill, Liewald, Cullen, Penninx, Boomsma, Pell, Mcintosh, Smith, Deary, Harris and Gwas2016).

The combined analysis is complementary to the study of individual phenotypes. Based on previous studies and as discussed earlier, depression, anxiety, and neuroticism may have some common biological basis, although most likely their etiologies do not completely overlap. A combined analysis may improve the power to detect drugs which target the shared pathophysiological processes, by increasing the total sample size. On the other hand, study of individual phenotypes may reveal drugs or drug classes that are specifically useful for particular disease symptoms or subgroups, for example, melancholic depression.

We performed meta-analysis of p values based on two methods, the Simes’ method (Simes, Reference Simes1986) and the Brown's approach (Brown, Reference Brown1975; Poole et al., Reference Poole, Gibbs, Shmulevich, Bernard and Knijnenburg2016). The Simes’ method is valid under positive regression dependencies (Sarkar and Chang, Reference Sarkar and Chang1997). Brown's method is similar to Fisher's method but also accounts for dependencies in p values. Briefly, assuming k p values, Fisher showed that the statistic $T = \sum {-2\log P_i} $ should follow a χ2 distribution with 2k degrees of freedom if the p values are independent. Brown's method is an extension of Fisher's approach by estimating the statistic T with a re-scaled χ2 distribution $T\sim c{\rm \chi} _{2f}^2 $ where

$$f = \displaystyle{{E{(T)}^2} \over {{\mathop{\rm var}} (T)}}\; \quad {\rm and} \;c\, = \,k/f.$$

The expectation of the statistic T can be estimated by E(T) = 2k and the variance by ${\mathop{\rm var}} (T) = 4k + 2\sum\limits_{i \lt j} {{\mathop{\rm cov}} (-2\log P_i,-2\log P_j)} $. We estimated the covariance empirically from the observed p value vectors of different phenotypes. The Simes’ method is more powerful when there are few false null hypotheses (H0) (i.e. few true associations) whereas Fisher's or Brown's method are more powerful when there are multiple false H0.

The MDD-CONVERGE sample includes only Chinese subjects and does not overlap with other datasets, but we accounted for overlapping samples (with Brown's or Simes’ method) in the remaining GWAS studies.

Testing for enrichment of psychiatric and other drug classes

We considered three sources when defining psychiatric drug-sets in our analyses. The first set came from drugs listed in the Anatomical Therapeutic Classification (ATC) system downloaded from KEGG. We extracted three groups of drugs: (1) all psychiatric drugs (coded ‘N05’ or ‘N06’); (2) antipsychotics (coded ‘N05A’); (3) antidepressants and anxiolytics (coded ‘N05B’ or ‘N06A’). The second source was from MEDication Indication resource (MEDI) (Wei et al., Reference Wei, Cronin, Xu, Lasko, Bastarache and Denny2013) derived from four public medication resources, namely RxNorm, Side Effect Resource 2 (SIDER2), Wikipedia, and MedlinePlus. A random subset of the extracted indications was checked by physicians. The MEDI high-precision subset (MEDI-HPS), with an estimated curation precision of 92%, was used in our analyses (Wei et al., Reference Wei, Cronin, Xu, Lasko, Bastarache and Denny2013). Since only known drug indications are included in ATC or MEDI-HPS, we also included an expanded set of drugs that are considered for clinical trials (as listed on https://clinicaltrials.gov). These drugs are usually promising candidates supported by preclinical or clinical studies. A precompiled list of these drugs was obtained from https://doi.org/10.15363/thinklab.d212. We also examined enrichment for closely related disorders in combinations, including schizophrenia with bipolar disorder (BD), as well as depression with anxiety.

The above represents a hypothesis-driven analysis of psychiatric drug classes. To explore whether other drug groups may be repositioned for disease treatment, we also performed a comprehensive enrichment analysis of all ATC level 3 drug classes. To avoid the results being driven by too few drugs in a class, we limited the analyses to drug classes with at least five members. Note that in the previous step we have found individual drugs as repositioning candidates, but we also hope to find out which drug classes may be promising for repositioning, which can also provide insights into potentially new mechanisms of actions in future development. A total of 129 and 79 ATC drug classes (for gene-sets derived from the entire DSigDB and PubChem/ChEMBL, respectively) were included in the final analyses.

We performed enrichment tests of repositioning hits for known drug classes, in a manner similar to the GSA described above. The p values are first converted to z-statistics, and the mean z-score within each drug class is compared against the theoretical null of zero (self-contained test) and against other drugs outside the designated drug class (competitive test) with one-sided tests.

In the context of GSA, the self-contained test examines the null hypothesis that none of the genes in the gene-set are associated with the phenotype, while the competitive test examines the null hypothesis that genes in the set are no more strongly associated with the phenotype than those outside the set. If only a small number of genes contribute to the phenotype, the difference between the two tests is small. However, if the trait is highly polygenic (which may be the case for many complex traits), the competitive test is more appropriate and usually preferred for revealing the biological processes involved. Using clinical trials as an analogy, self-contained analysis is akin to a treatment-only design, whereas the competitive analysis is akin to a treatment versus control design.

Here we are testing for enrichment of sets of drugs, and there is one complication to be noted. It is reasonable to believe that the current antidepressants or anxiolytics are not the only drugs that have therapeutic effects; in other words, a certain proportion of drugs in the ‘competing set’ might also have therapeutic potential against depression or anxiety. Therefore, results of the competitive tests should be interpreted with this potential limitation in mind. In this paper, we presented the drug-set enrichment results of both self-contained and competitive tests. Ideally, the results will be significant regardless of the type of analysis (i.e. large ‘treatment’ effects as well as an edge over the ‘control’ group, or in our case, significant enrichment for each psychiatric drug gene-set as well as stronger enrichment over non-psychiatric drug gene-sets).

Literature support of results

We extracted the top 20 repositioning hits of each psychiatric trait and meta-analyzed results, with drug-related gene-sets derived from either the entire database or drug targets defined by PubChem or ChEMBL. We performed a systematic search in PubMed and Google scholar using the following terms: Drug_name AND (depression OR depressive OR antidepressant OR anxiety OR panic OR phobia OR anxiolytic). References therein were looked up as necessary. The search was performed in June to August 2017 and March 2018 (for MDD-2018 results).

A limitation of manual search is that assessment of all drugs is almost impossible due to the time involved. As will be discussed later, we revealed literature support for the top repositioning candidates, but it could be argued that evidence of support may also be found for lower ranking drugs. One may also wish to evaluate whether drugs with less (or no) support by the literature are associated with weaker significance (i.e. higher p values) in our analysis. These drugs are roughly analogous to ‘negative controls’ in an experiment.

As a complementary strategy to manual literature search of the top results, we conducted an ‘automated’ literature search on all drugs in PubMed in an unbiased manner. We extracted the number of research articles supporting each drug's association with depression or AD using a Python script. To enhance specificity of the search, we employed MeSH terms for the search using the following terms: [‘drug_name’(Title/Abstract)] AND [anxiety(MeSH Terms) OR anxiety disorders(MeSH Terms) OR depressive disorder major(MeSH Terms) OR depressive disorder(MeSH Terms)] AND [therapeutics(MeSH Terms) OR antidepressive agents(MeSH Terms) OR anti-anxiety agents(MeSH Terms)]. We then examined the correlation between the number of research articles of support and p value for each drug. We hypothesize that drugs with lower p values (i.e. stronger candidates from our analysis) will correspond to a greater number of articles retrieved from PubMed, leading to negative correlations. Similar approaches for validating repositioning candidates based on literature mining have also been used in other studies (e.g. Huang et al., Reference Huang, Li, Sheng, Xia, Ma, Zhan and Wong2014). As the number of articles is typically skewed and not normally distributed, we employed Spearman and Kendall correlation measures, which are based on ranks of observations. In addition, we compared the p values of drugs with no article support versus those with at least one article support, and hypothesized that the former group would have higher p values. We used the Wilcoxon rank-sum test for such comparison. Statistical tests were carried out in R and tests were one-tailed.

We note that the number of research articles is only a crude proxy of the level of literature support. This is considered a secondary analysis and the results should be interpreted with these limitations in mind.

Correction for multiple testing

We employed the false discovery rate (FDR) approach to account for multiple testing (Benjamini and Hochberg, Reference Benjamini and Hochberg1995). The FDR approach controls the expected proportion of false positives among those declared to be significant (with consideration of the multiple testing performed). It has been widely adopted in genomic studies (So and Sham, Reference So and Sham2011). As a simple example, if there are 100 hypotheses with ‘FDR-adjusted p value’ ⩽0.05, we would expect that among these 100 hypotheses, there are ⩽5% (i.e. five) false positives. ‘FDR-adjusted p values’ (also known as ‘q-values’) were computed by the R function p.adjust with the Benjamini–Hochberg procedure (Benjamini and Hochberg, Reference Benjamini and Hochberg1995). The primary threshold for FDR-adjusted p values was set at 0.05, while results with FDR-adjusted p between 0.05 and 0.1 were regarded as suggestive associations.

Results

Enrichment of psychiatric drug classes among the drugs repositioned from GSA

Tables 1–3 show the FDR-adjusted p values for enrichment of major psychiatric drug classes among the drugs repositioned from GSA. Gene-sets were derived from the entire DSigDB. We observed that the drugs repositioned from most GWAS of anxiety and depressive traits are generally enriched for known psychiatric medications.

Table 1. Enrichment of repositioning hits derived from GWAS of major depressive disorder (MDD) and depressive symptoms for psychiatric drug categories (FDR-adjusted p values presented; gene-sets derived from entire DSigDB)

Note that we present the FDR-adjusted p values, which are already corrected for multiple testing. Gene-sets were derived from the entire DSigDB. We used three databases to define psychiatric drug categories. ATC and MEDI-HPS records known psychiatric drugs while clinicalTrial.gov records the drugs that were tested in clinical trials. Test results with FDR < 0.05 are in bold. Results with FDR between 0.05 and 0.1 are in italics.

Self: self-contained test; Compet, competitive test. In the context of gene-set analysis, the self-contained test examines the null hypothesis that none of the genes in the gene-set are associated with the phenotype, while the competitive test examines the null hypothesis that genes in the set are no more strongly associated with the phenotype than those outside the set. MDD-CONVERGE, MDD with GWAS data from the CONVERGE Consortium; MDD with GWAS data from the Psychiatric Genomics Consortium; DepSym, GWAS of depressive symptoms from the Social Science Genetics Association Consortium (SSGAC).

Table 2. Enrichment of repositioning hits derived from GWAS of anxiety disorders and neuroticism for psychiatric drug categories (FDR-adjusted p values presented; gene-sets derived from entire DSigDB)

The FDR-adjusted p values are presented, which are already corrected for multiple testing. Gene-sets were derived from the entire DSigDB.

Test results with FDR < 0.05 are in bold. Results with FDR between 0.05 and 0.1 are in italics. Self: self-contained test; Compet, competitive test. Anxiety CC, GWAS of anxiety disorders case–control sample; neurotic, GWAS of neuroticism in general population.

Table 3. Enrichment of repositioning hits derived from meta-analyzed GWAS p values for psychiatric drug categories (FDR-adjusted p values presented; gene-sets derived from entire DSigDB)

Self, self-contained test; Compet, competitive test.

The presented Brown and Simes p values are derived from meta-analyzing all five datasets, i.e. MDD-CONVERGE, MDD2018, DepSym, AnxietyCC, and neuroticism.

The FDR-adjusted p values are presented, which are already corrected for multiple testing. Gene-sets were derived from the entire DSigDB.

Test results with FDR < 0.05 are in bold. Results with FDR between 0.05 and 0.1 are in italics.

First we consider the datasets (MDD-CONVERGE, MDD2018, SSGAC-DS) which focus on depression traits (Table 1). We observed significant enrichment for various psychiatric drug classes in all the above datasets. The MDD-CONVERGE sample showed significant enrichment in both self-contained and competitive analyses. Significant results were seen for antipsychotics and antidepressants/anxiolytics within ATC and MEDI-HPS categories. For SSGAC-DS, we observed enrichment of drugs for schizophrenia and BD, and suggestive associations with medications for anxiety and depression listed in clinicalTrial.gov. For the MDD2018 sample, we also observed enrichment for antipsychotics and most psychiatric drug groups listed in clinicalTrial.gov. As a subsidiary analysis, we also examined whether there is any enrichment for psychiatric drugs based on the much smaller MDD-PGC-2012 sample, but no significant results were observed (online Supplementary Table S4)

As for GWAS studies on neuroticism (SSGAC-NEU) and AD (Table 2), there was evidence of enrichment for most psychiatric drug classes under study. Interestingly, for neuroticism, the strongest enrichment was for antipsychotics (lowest FDR-adjusted p = 2.28E-09) instead of antidepressants.

For analyses involving meta-analyzed GWAS data across all datasets (Table 3), enrichment was observed for all psychiatric drug classes, with generally stronger or at least comparable statistical associations when compared with enrichment tests of individual GWAS. The results of Brown's and Simes’ tests were largely consistent with each other.

We also performed an analysis based on drug target genes derived from PubChem or ChEBML only. The results (FDR-adjusted p values) were given in online Supplementary Tables S1–S3 and were largely similar to the above findings.

Table 4 and online Supplementary Table S5 show the results of enrichment tests across all ATC level 3 drug classes based on gene-sets derived from the whole DSigDB database. Table 5 and online Supplementary Table S6 show the findings from the same analysis but with gene-sets derived from PubChem or ChEBML only. We consider results from both of the analyses here. (Note that we mainly focus on the top 5 ATC drug classes, in contrast to top 20 repositioning candidates in the subsequent section, as the number of drug classes considered was smaller than the number of drugs.) For MDD-CONVERGE, drug classes pertaining to corticosteroids were ranked highly. Notably, antidepressants and antipsychotics were ranked within the top 5 in our two sets of analyses. For MDD-PGC, anti-infective and anti-inflammatory agents were significantly enriched. Interestingly, hypnotics/sedatives and anxiolytics, which were commonly used as short-term therapies especially at the acute phase of illness (Kanba, Reference Kanba2004), are also ranked among the top. Antidepressants or other psychiatric drugs however were not enriched. For MDD2018, interestingly, antiepileptics was ranked the highest. Other top-ranked drug classes that may be clinically relevant include antimigraine preparations, hypnotics, and dopaminergic agents.

Table 4. Top 5 enriched ATC drug classes (gene-sets derived from the whole DSigDB database)

The FDR-adjusted p values are presented, which are already corrected for multiple testing. Test results with FDR < 0.05 are in bold. Results with FDR between 0.05 and 0.1 are in italics. Self, self-contained test; Compet, competitive test. The definitions of self-contained and competitive tests are very similar to those in Tables 1–3. In the context, the above ‘drug-group’ analysis, the self-contained test examines the null hypothesis that none of the drugs in the group are associated with repositioning potential, while the competitive test examines the null hypothesis that drugs inside the group are no more significantly associated with repositioning potential than those outside the group.

Table 5. Top 5 enriched ATC drug classes with gene-sets derived from PubChem and ChEMBL only

The FDR-adjusted p values are presented, which are already corrected for multiple testing. Test results with FDR <0.05 are in bold. Results with FDR between 0.05 and 0.1 are in italics. Self, self-contained test; Compet, competitive test. The definitions of self-contained and competitive tests are very similar to those in Tables 1–3. In the context the above ‘drug-group’ analysis, the self-contained test examines the null hypothesis that none of the drugs in the group are associated with repositioning potential, while the competitive test examines the null hypothesis that drugs inside the group are no more significantly associated with repositioning potential than those outside the group.

For AD, antipsychotics and antidepressants were the two most strongly enriched medication classes. The β-blocking agents were ranked among the top for AD. The β-adrenergic blockers, especially propranolol, have been clinically used in AD for a long time, although the efficacy still remains uncertain (Steenen et al., Reference Steenen, Van Wijk, Van Der Heijden, Van Westrhenen, De Lange and De Jongh2016). As for depressive symptoms, dopaminergic and antiepileptic agents were listed among the top, and interestingly lipid-lowering agent was also on the top list. As for neuroticism, antipsychotics were the most strongly enriched drug class, and antidepressants were also ranked within the top 5. Other drug classes listed in top 5 also included hypnotics/sedatives, anxiolytics, and dopaminergic agents. The full results are shown in online Supplementary Tables S5 and S6.

We also examined the presence and extent of overlap between the top-listed ATC drug categories. The results are presented in online Supplementary Table S7. In general, there is a little overlap across most categories and almost zero overlap for categories starting with N (neuropsychiatric drugs). Overlap was concentrated in categories of anti-inflammatory agents, anti-infective agents, and corticosteroids.

Correlation of level of drug p values with degree of literature support

Online Supplementary Table S8 shows the correlation between p value of each drug (derived from GSA) with the number of research articles retrieved (gene-sets derived from the whole DSigDB database). All results were statistically significant and passed FDR correction for multiple testing (FDR<0.05), suggesting that drugs with lower p values from our GSA may have stronger literature support, and those with higher p values have weaker support. The strongest correlations were observed for combined analysis using Simes’ or Brown's method. When restricted to gene-sets from PubChem and ChEMBL (online Supplementary Table S9), the correlations were weaker, with significant results (FDR<0.05) observed for depressive symptoms, MDD2018, and combined p values using Brown's method.

However, we observed that the levels of correlation were only modest (<0.1). This analysis alone may not be sufficient to provide evidence for the validity of repositioning candidates; however, it may be considered a complementary validation strategy for our repositioning approach as a whole (as compared with manual literature review of the top candidates), as this analysis is more quantitative, unbiased, and comprehensive (covering all drugs).

Top repositioning hits and literature support from previous studies

Some of the top hits are mentioned below. However, we emphasize that these findings are tentative and require validations in further experimental and clinical studies. In addition, although some of the drugs received support from the literature, a certain degree of publication bias is likely present. Also, our literature search is not meant to be a systematic review due to the number of drugs involved.

Based on gene-sets derived from the whole DSigDB database

We found several interesting repositioning hits that were supported by previous studies (Table 6; online Supplementary Table S10 gives a fully annotated list of the top 20 drugs). The top-ranked repositioning hit identified in the meta-analysis was fendiline (Brown's p = 1.06E-11, q = 5.55E-8), a non-selective calcium channel blocker (CCB). Fendiline was recently shown to exert antidepressant-like effects in a mouse model by inhibition of acid sphingomyelinase activity and reduction of ceramide concentrations in the hippocampus (Gulbins et al., Reference Gulbins, Palmada, Reichel, Luth, Bohmer, Amato, Muller, Tischbirek, Groemer, Tabatabai, Becker, Tripal, Staedtler, Ackermann, Van Brederode, Alzheimer, Weller, Lang, Kleuser, Grassme and Kornhuber2013). A drop in ceramide concentrations might lead to increased neurogenesis and improved neuronal maturation and survival (Gulbins et al., Reference Gulbins, Palmada, Reichel, Luth, Bohmer, Amato, Muller, Tischbirek, Groemer, Tabatabai, Becker, Tripal, Staedtler, Ackermann, Van Brederode, Alzheimer, Weller, Lang, Kleuser, Grassme and Kornhuber2013). The drug was ranked among the top for MDD-PGC-2012, MDD2018, SSGAC-DS and neuroticism.

Table 6. Selected repositioning hits (ranked within top 20 for a phenotype) with literature support

FDR-adj p, FDR-adjusted p values.

Here the repositioning hits are listed regardless of whether the underlying gene-sets were derived from the entire DsigDB or from PubChem/ChEBML only. Results from meta-analysis (Simes’ or Brown's method) are not specifically listed as they are included in the results of individual phenotypes.

Please also refer to the legends of Table 1 and the main text for the references and more detailed discussions.

As for individual phenotypes, some interesting candidates for MDD-CONVERGE include DHEA, a neurosteroid with evidence of antidepressive effects in a double-blind randomized controlled trial (RCT); naringenin chalcone, a citrus bioflavonoid shown to be effective in mouse models (Yi et al., Reference Yi, Liu, Li, Luo, Liu, Geng, Tang, Xia and Wu2014); amoxapine, a tetracyclic antidepressant (Jue et al., Reference Jue, Dawson and Brogden1982). Other candidates listed for MDD-PGC or MDD2018 include ibuprofen, an NSAID shown to be associated with lower depressive symptoms in osteoarthritis patients; piperlongumine, a constituent of the fruit of Piper longum, which was shown to confer resistance against stress in a mouse model (Yadav et al., Reference Yadav, Chatterjee, Majeed and Kumar2015); and sanguinarine, a selective mitogen-activated protein kinase phosphatase-1 (Mkp-1) inhibitor which produced antidepressant-like effect in rats (Chen et al., Reference Chen, Wang, Zhang, Wang, Peng, Sun and Tan2012).

For AD, candidates included a known antidepressant trazodone and a serotonin agonist quipazine that may increase brain serotonin levels (Fuller et al., Reference Fuller, Snoddy, Perry, Roush, Molloy, Bymaster and Wong1976); for depressive symptoms, alsterpaullone, a glycogen synthase kinase-3β (GSK-3β) inhibitor (Leost et al., Reference Leost, Schultz, Link, Wu, Biernat, Mandelkow, Bibb, Snyder, Greengard, Zaharevitz, Gussio, Senderowicz, Sausville, Kunick and Meijer2000). Increased activation of GSK-3β was also associated with depression-like behavior in mouse models, which could be alleviated by GSK-3β inhibitors (Beaulieu et al., Reference Beaulieu, Zhang, Rodriguiz, Sotnikova, Cools, Wetsel, Gainetdinov and Caron2008). In addition, inhibition of GSK3 has been postulated as a major mechanism of action by the mood stabilizer lithium (Freland and Beaulieu, Reference Freland and Beaulieu2012). A NMDA receptor antagonist H-7 was also top-listed (Amador and Dani, Reference Amador and Dani1991).

Among the top repositioning hits from the individual and meta-analysis results, a number of them are CCB. These include fendiline, perhexiline, prenylamine, and felodipine [prenylamine was withdrawn from the market due to risk of QT prolongation and torsades de pointes (Fung, Reference Fung2001)]. Although with the exception of fendiline, no direct experimental or clinical studies have shown antidepressant or anxiolytic properties of the above drugs, CCB as a whole has been proposed as treatment for various psychiatric disorders. CCB has been mostly studied for the treatment of mania, recently reviewed in Cipriani et al. (Reference Cipriani, Saunders, Attenburrow, Stefaniak, Panchal, Stockton, Lane, Tunbridge, Geddes and Harrison2016). However, the number of quality double-blind RCTs was small, and there is yet no sufficient evidence to suggest the use of CCB in treating manic symptoms. As for depression, a recent pilot (patient-only) study of isradipine on bipolar depression showed positive results (Ostacher et al., Reference Ostacher, Iosifescu, Hay, Blumenthal, Sklar and Perlis2014). Another CCB, nicardipine, was reported to enhance the antidepressant action of electroconvulsive therapy (Dubovsky et al., Reference Dubovsky, Buzan, Thomas, Kassner and Cullum2001). Notwithstanding the mixed evidence, CCB are probably still worthy of further investigation for depression and AD, given the biological relevance of calcium signaling, preliminary support from clinical studies, and the wide availability as well as known safety profiles of this drug class.

Based on gene-sets derived from PubChem and ChEMBL

A manually annotated list of the top 20 prioritized drugs is given in online Supplementary Table S11. Some highlighted candidates are also shown in Table 6. Some repositioning hits for MDD-CONVERGE included oxotremorine-M, a muscarinic acetylcholine receptor agonist that may ameliorate depressive symptoms and restore hippocampal neurogenesis in an animal model (Veena et al., Reference Veena, Srikumar, Mahati, Raju and Shankaranarayana Rao2011) and capeserod, a 5-HT4 receptor agonist shown effective in a mouse model (Tamburella et al., Reference Tamburella, Micale, Navarria and Drago2009). For MDD-PGC-2012, top-listed drugs include three NSAIDs, namely ibuprofen, meclofenamate, and flufenamic acid, which possess anti-inflammatory actions that may also be beneficial for depression (Iyengar et al., Reference Iyengar, Gandhi, Aneja, Thorpe, Razzouk, Greenberg, Mosovich and Farkouh2013; Zaminelli et al., Reference Zaminelli, Gradowski, Bassani, Barbiero, Santiago, Maria-Ferreira, Baggio and Vital2014). Remarkably, our analysis recovered a known augmentation therapy for depression, namely Cytomel [or synthetic triiodothyronine (T3)] (Touma et al., Reference Touma, Zoucha and Scarff2017). Another candidate was guanidinonaltrindole di-trifluoroacetate, a κ-opioid receptor inhibitor with antidepressant-like and anxiolytic-like efficacy in rat models (Peters et al., Reference Peters, Zacco, Gordon, Maciag, Litwin, Thompson, Schroeder, Sygowski, Piser and Brugel2011). For MDD2018, resveratrol was reported to have antidepressant-like effects (de Oliveira et al., Reference De Oliveira, Chenet, Duarte, Scaini and Quevedo2017) and is currently under study in an RCT (https://clinicaltrials.gov/ct2/show/NCT03384329); valproate was also on the list which is effective against bipolar depression (Smith et al., Reference Smith, Cornelius, Azorin, Perugi, Vieta, Young and Bowden2010). Preliminary evidence from a small clinical study showed possible benefits for treatment-resistant depression (Ghabrash et al., Reference Ghabrash, Comai, Tabaka, Saint-Laurent, Booij and Gobbi2016).

For depressive symptoms, the top-listed candidates include prochlorperazine and raclopride, which are antipsychotics; protriptyline, a tricyclic antidepressant; oleamide, a cannabinoid receptor type 1 (CB1) receptor agonist of potential antidepressant effects in animal models (Fedorova et al., Reference Fedorova, Hashimoto, Fecik, Hedrick, Hanus, Boger, Rice and Basile2001; Wei et al., Reference Wei, Yang, Dong and Wu2007); ellagic acid and chebulinic acid (an ellagitannin), natural phenol antioxidants with some evidence of improving depressive traits again in animal models (Dhingra and Chhillar, Reference Dhingra and Chhillar2012; Girish et al., Reference Girish, Raj, Arya and Balakrishnan2013; Onasanwo et al., Reference Onasanwo, Faborode, Agrawal, Ijiwola, Jaiyesimi and Narender2014). Also of note is the nicotinic agonist DMPP, as nicotine has been shown to have antidepressant properties in pre-clinical and clinical studies (Philip et al., Reference Philip, Carpenter, Tyrka and Price2012); the antidepressant effects may be due to initial activation of the nicotinic receptor followed by desensitization leading to long-term antagonism (Philip et al., Reference Philip, Carpenter, Tyrka and Price2012). For neuroticism, again quite a number of top hits are CCB, and the CCB fendiline was the top-ranked candidate. The repositioning hits derived from meta-analyzed p values were covered above.

There were a number of noteworthy candidates from the analysis results of AD. Notably, the top-ranked drugs contain pregabalin and gabapentin, two drugs with similar mechanisms which have been shown anxiolytic effects in clinical trials (especially pregabalin) (Pande et al., Reference Pande, Crockatt, Feltner, Janney, Smith, Weisler, Londborg, Bielski, Zimbroff, Davidson and Liu-Dumaw2003; Rickels et al., Reference Rickels, Pollack, Feltner, Lydiard, Zimbroff, Bielski, Tobias, Brock, Zornberg and Pande2005); clomipramine, a known tricyclic antidepressant; and trazodone, another known antidepressant belonging to the serotonin antagonist and reuptake inhibitor class. Some other noteworthy candidates include ritanserin, a selective 5-HT2A and 5-HT2C receptor antagonist (Ceulemans et al., Reference Ceulemans, Hoppenbrouwers, Gelders and Reyntjens1985; Nappi et al., Reference Nappi, Sandrini, Granella, Ruiz, Cerutti, Facchinetti, Blandini and Manzoni1990; Strauss and Klieser, Reference Strauss and Klieser1991); pindolol, a β-blocker with partial β-adrenergic receptor agonist activity and a possible augmentation therapy for depression (Blier and Bergeron, Reference Blier and Bergeron1995; Artigas et al., Reference Artigas, Celada, Laruelle and Adell2001); atomoxetine, a norepinephrine reuptake inhibitor.

Discussion

In this study, we leveraged large-scale GWAS summary data and analyzed gene-sets associated with drugs to uncover repositioning opportunities for depression and AD. It is encouraging that we observed significant enrichment for known psychiatric medications or drugs considered in clinical trials. Remarkably, antipsychotics and antidepressants were the two most significantly enriched drug classes in our combined analysis (with gene-sets derived from the whole DSigDB). Our findings provide support for the validity of GSA in drug repurposing. In addition, we revealed a number of interesting candidates for repurposing, some of which were supported by prior studies. Although relatively few susceptibility variants of genome-wide significance have been found for depression and AD, our findings suggest that leveraging variants with weaker associations, for example, by GSA, might still contribute valuable information to the prioritization of new drug candidates.

While we have included several datasets (SSGAC-DS, MDD2018, MDD-CONVERGE) directly related to depression, the enrichment results are not entirely the same. MDD2018 and MDD-CONVERGE are case–control GWAS studies on MDD; however, the drugs and drug classes enriched were quite different. Another noteworthy finding is that despite similar sample sizes, MDD-CONVERGE demonstrated much stronger enrichment for psychiatric drugs than MDD-PGC-2012, suggesting its better value in guiding drug discoveries. The discrepancy in results between MDD-CONVERGE and other samples might be due to different characteristics of the recruited subjects. The MDD-CONVERGE sample differs from other MDD GWAS samples by gender, ethnicity, and other characteristics of the case and control subjects (Flint et al., Reference Flint, Chen, Shi, Kendler and Consortium2012; Cai et al., Reference Cai, Bigdeli, Kretzschmar, Li, Liang, Song, Hu, Li, Jin, Hu, Wang, Wang, Qian, Liu, Jiang, Lu, Zhang, Yin, Li, Xu, Gao, Reimers, Webb, Riley, Bacanu, Peterson, Chen, Zhong, Liu, Wang, Sun, Sang, Jiang, Zhou, Li, Li, Zhang, Wang, Fang, Pan, Miao, Zhang, Hu, Yu, Du, Sang, Li, Chen, Cai, Yang, Yang, Ha, Hong, Deng, Li, Li, Song, Gao, Zhang, Gan, Meng, Pan, Gao, Zhang, Sun, Li, Niu, Zhang, Liu, Hu, Zhang, Lv, Dong, Wang, Tao, Wang, Xia, Rong, He, Liu, Huang, Mei, Shen, Liu, Shen, Tian, Liu, Wu, Gu, Fu, Shi, Chen, Gan, Liu, Wang, Yang, Cong, Marchini, Yang and Wang2015). The CONVERGE sample only included women, and cases were required to have two or more depressive episodes. The proportion of patients with melancholia is ~85%, which is higher than the prevalence in a general MDD sample (~23.5%) (McGrath et al., Reference Mcgrath, Khan, Trivedi, Stewart, Morris, Wisniewski, Miyahara, Nierenberg, Fava and Rush2008). To reduce the probability that some patients initially diagnosed with MDD may go on to develop BD, cases were required to be at least 30 years of age. Regarding control subjects, all controls had been screened by personal interview to exclude history of depressive episode(s). There was also an age limit for controls (>40), to reduce the chance that some of them may develop MDD at later ages. Finally, all four grandparents of cases and controls were required to be of Chinese ancestry, in order to maximize genetic homogeneity.

It is widely accepted that MDD is a heterogeneous disorder, with a variety of clinical presentations and possibly divergent pathophysiologies (Goldberg, Reference Goldberg2011). By recruiting a more homogeneous group of patients, the CONVERGE study might have better power in detecting susceptibility genes despite a lower sample size. Indeed, MDD-CONVERGE revealed two genome-wide significant loci while none was found in the MDD-PGC-2012 study although the sample sizes were similar. It is also worth mentioning that previous meta-analyses showed that the response to antidepressant depends on the baseline severity of depression (Kirsch et al., Reference Kirsch, Deacon, Huedo-Medina, Scoboria, Moore and Johnson2008; Fournier et al., Reference Fournier, Derubeis, Hollon, Dimidjian, Amsterdam, Shelton and Fawcett2010). The studies reported that effects of antidepressants were largest for the most severely depressed group, but smaller for mild-to-moderate depression. Our findings, although based on a different study paradigm, are broadly in line with this clinical observation.

Given that the prioritized candidates from MDD-CONVERGE and other samples of European ancestry are different, an interesting question is whether there may exist ‘ethnicity-specific drugs’, or at least drug efficacy differing by ethnicities. This is an interesting direction for future work, but our present understanding of MDD and the data at hand do not allow a concrete conclusion to be drawn yet. While the CONVERGE sample differs from other samples in ethnicity (Chinese versus European ancestry), there are also marked differences in patient characteristics and other ascertainment criteria as detailed above. It is therefore difficult to tease out whether the differences in results are due to ethnic-specific factors alone. Our current understanding of MDD in different ethnic groups is still relatively limited. Also, few studies have explored differences in genetic basis of MDD across ethnicities. A recent work reported statistically significant but only partial genetic correlations between MDD in Chinese and European populations (0.33 for lifetime MDD; 0.40 and 0.41 for female-only and recurrent MDD, respectively; the Chinese sample was based on MDD-CONVERGE) (Bigdeli et al., Reference Bigdeli, Ripke, Peterson, Trzaskowski, Bacanu, Abdellaoui, Andlauer, Beekman, Berger, Blackwood, Boomsma, Breen, Buttenschon, Byrne, Cichon, Clarke, Couvy-Duchesne, Craddock, De Geus, Degenhardt, Dunn, Edwards, Fanous, Forstner, Frank, Gill, Gordon, Grabe, Hamilton, Hardiman, Hayward, Heath, Henders, Herms, Hickie, Hoffmann, Homuth, Hottenga, Ising, Jansen, Kloiber, Knowles, Lang, Li, Lucae, Macintyre, Madden, Martin, Mcgrath, Mcguffin, Mcintosh, Medland, Mehta, Middeldorp, Milaneschi, Montgomery, Mors, Muller-Myhsok, Nauck, Nyholt, Nothen, Owen, Penninx, Pergadia, Perlis, Peyrot, Porteous, Potash, Rice, Rietschel, Riley, Rivera, Schoevers, Schulze, Shi, Shyn, Smit, Smoller, Streit, Strohmaier, Teumer, Treutlein, Van Der Auwera, Van Grootheest, Van Hemert, Volzke, Webb, Weissman, Wellmann, Willemsen, Witt, Levinson, Lewis, Wray, Flint, Sullivan and Kendler2017). With regards to drug response, Murphy et al. (Reference Murphy, Hou, Maher, Woldehawariat, Kassem, Akula, Laje and Mcmahon2013) reported that patients of African ancestry had poorer response to treatment based on the STAR*D trial; genetic ancestry explained a significant proportion of residual variance even after accounting for socioeconomic and baseline clinical differences. However, some studies reported no racial group differences in treatment response (Lesser et al., Reference Lesser, Myers, Lin, Bingham Mira, Joseph, Olmos, Schettino and Poland2010; Lesser et al., Reference Lesser, Zisook, Gaynes, Wisniewski, Luther, Fava, Khan, Mcgrath, Warden, Rush and Trivedi2011). However, to our knowledge, there are no similar studies comparing Chinese with other ethnicities. Overall, we believe that further studies are required to elucidate the differences and similarities in clinical/genetic determinants of depression and drug response across ethnic groups.

It is noteworthy that the repositioning hits are not only enriched for antidepressants or anxiolytics but also antipsychotics. A meta-analysis by Spielmans et al. revealed that atypical antipsychotics are effective as adjunctive treatment for treatment-resistant depression (Spielmans et al., Reference Spielmans, Berman, Linardatos, Rosenlicht, Perry and Tsai2013). Zhou et al. also reached a similar conclusion in a recent network meta-analysis (Zhou et al., Reference Zhou, Keitner, Qin, Ravindran, Bauer, Del Giovane, Zhao, Liu, Fang, Zhang and Xie2015). Atypical antipsychotics may also be useful for AD and symptoms (Gao et al., Reference Gao, Muzina, Gajwani and Calabrese2006; LaLonde and Van Lieshout, Reference Lalonde and Van Lieshout2011; Samuel et al., Reference Samuel, Zimovetz, Gabriel and Beard2011), although further studies are required and that the benefits need to be balanced against the side-effects. Furthermore, a shared genetic basis between schizophrenia and depression is well-established (Cross-Disorder Group of the Psychiatric Genomics et al., Reference Lee, Ripke, Neale, Faraone, Purcell, Perlis, Mowry, Thapar, Goddard, Witte, Absher, Agartz, Akil, Amin, Andreassen, Anjorin, Anney, Anttila, Arking, Asherson, Azevedo, Backlund, Badner, Bailey, Banaschewski, Barchas, Barnes, Barrett, Bass, Battaglia, Bauer, Bayes, Bellivier, Bergen, Berrettini, Betancur, Bettecken, Biederman, Binder, Black, Blackwood, Bloss, Boehnke, Boomsma, Breen, Breuer, Bruggeman, Cormican, Buccola, Buitelaar, Bunney, Buxbaum, Byerley, Byrne, Caesar, Cahn, Cantor, Casas, Chakravarti, Chambert, Choudhury, Cichon, Cloninger, Collier, Cook, Coon, Cormand, Corvin, Coryell, Craig, Craig, Crosbie, Cuccaro, Curtis, Czamara, Datta, Dawson, Day, De Geus, Degenhardt, Djurovic, Donohoe, Doyle, Duan, Dudbridge, Duketis, Ebstein, Edenberg, Elia, Ennis, Etain, Fanous, Farmer, Ferrier, Flickinger, Fombonne, Foroud, Frank and Franke2013), and a recent study also found significant genetic correlation between neuroticism and schizophrenia (Smith et al., Reference Smith, Escott-Price, Davies, Bailey, Colodro-Conde, Ward, Vedernikov, Marioni, Cullen, Lyall, Hagenaars, Liewald, Luciano, Gale, Ritchie, Hayward, Nicholl, Bulik-Sullivan, Adams, Couvy-Duchesne, Graham, Mackay, Evans, Smith, Porteous, Medland, Martin, Holmans, Mcintosh, Pell, Deary and O'Donovan2016). Epidemiology studies also demonstrated associations of neuroticism with schizophrenia (Malouff et al., Reference Malouff, Thorsteinsson and Schutte2005).

Antidepressants and antipsychotics were the most consistently associated categories based on our analysis. Here we also briefly discuss a few other drug categories listed among the top. However, we caution that the other drug categories are less well-studied for depression and anxiety, with considerably weaker support of therapeutic effect from the literature. They are also less consistently implicated in our analysis. While discussed below, we stress that our results are tentative and further validation studies are necessary to verify the findings.

For MDD-CONVERGE, which is composed of mainly severe melancholic depressive patients, we found that corticosteroids were highly ranked. A number of studies (Gold and Chrousos, Reference Gold and Chrousos2002; Stetler and Miller, Reference Stetler and Miller2011; Lamers et al., Reference Lamers, Vogelzangs, Merikangas, De Jonge, Beekman and Penninx2013) have shown that hyperactivity of the hypothalamic–pituitary axis is an important feature of melancholic depression. Nevertheless, as with other kinds of GSA, the current approach of GSA does not delineate the exact direction of drug effects. It may be concluded that corticosteroids or drugs targeting relevant pathways might be of clinical significance to melancholic depression; however, whether corticosteroids themselves or drugs blocking their actions would be useful remains to be investigated. It is also worth noting that currently there is no reliable approach to subtype depression yet. While some researchers propose melancholic features may be used to define subtypes of depression (Parker et al., Reference Parker, Fink, Shorter, Taylor, Akiskal, Berrios, Bolwig, Brown, Carroll, Healy, Klein, Koukopoulos, Michels, Paris, Rubin, Spitzer and Swartz2010), the biological and clinical validity of such classification remains to be established (Melartin et al., Reference Melartin, Leskela, Rytsala, Sokero, Lestela-Mielonen and Isometsa2004; Harald and Gordon, Reference Harald and Gordon2012).

We briefly highlight a few other drug classes here. Dopaminergic agents were ranked among the top for depressive symptoms; stimulants including dopaminergic agents have been tested in clinical trials for unipolar and bipolar depression (Nierenberg et al., Reference Nierenberg, Dougherty and Rosenbaum1998; Dell'Osso et al., Reference Dell'Osso, Ketter, Cremaschi, Spagnolin and Altamura2013; Corp et al., Reference Corp, Gitlin and Altshuler2014; Szmulewicz et al., Reference Szmulewicz, Angriman, Samame, Ferraris, Vigo and Strejilevich2017). However, further randomized controlled trials of large sample sizes are required to confirm the efficacy and clarify possible side-effects of such treatment. We also observed that anti-migraine medications were highly ranked in the list. An increased rate of depression among migraine patients is well-established (Antonaci et al., Reference Antonaci, Nappi, Galli, Manzoni, Calabresi and Costa2011), and antidepressants (mainly tricyclics) have been used for migraine prophylaxis (Koch and Jurgens, Reference Koch and Jurgens2009). The effects of anti-migraine medications on depressive symptoms however remain to be elucidated. Lipid-lowering agents was top-listed for depressive symptoms (based on all genes in DSigDB) and AD (based on drug target genes from PubChem and ChEMBL), and studies have shown potential benefits of statins in depression (Salagre et al., Reference Salagre, Fernandes, Dodd, Brownstein and Berk2016), and when used with concomitant SSRI (Kohler et al., Reference Kohler, Gasse, Petersen, Ingstrup, Nierenberg, Mors and Ostergaard2016). However, some studies have also reported increased depressive symptoms with statins (You et al., Reference You, Lu, Zhao, Hu and Zhang2013), hence the exact relationship may be complex and may differ by patient characteristics. Whether other types of lipid-lowering drugs may be beneficial is another topic worthy of further investigations.

In this study, we employed the GSA approach to drug repositioning. The current study is complementary to our recent repositioning attempt using a new framework in which the drug-induced transcriptome is compared against GWAS-imputed expression profiles (So et al., Reference So, Chau, Chiu, Ho, Lo, Yim and Sham2017). Each of these two methods has their own advantages and disadvantages. The methodology of finding reversed expression patterns [as detailed in So et al. (Reference So, Chau, Chiu, Ho, Lo, Yim and Sham2017)] has a unique advantage of accounting for the directions of associations. It also takes into account the functional impact of variants on expression and is intuitive from a biological point of view. While differentially expressed genes can be included in gene-sets, the actual (quantitative) expression changes are not considered which results in a loss of information. GSA also does not delineate the directions of effects. Nevertheless, GSA can directly make use of knowledge concerning known drug targets and other information on drug-related genes, for which more databases are available. Moreover, in the ‘reversed transcriptome’ approach (So et al., Reference So, Chau, Chiu, Ho, Lo, Yim and Sham2017), one only considers functional effects of the SNPs on expression, but genetic variants may also contribute to disease pathogenesis via other mechanisms, such as splicing or changes in protein function. The gene-based test in the current analysis is based on combining statistical evidence of individual SNPs in a general manner, and may capture a wider range of effects on disease risk. Also, the transcriptome comparison approach involves ‘imputing’ expression levels; since the major reference transcriptome dataset (GTEx) is mainly composed of Caucasians (84.6%) with greater proportion of males (65.6%) (https://www.gtexportal.org/home/tissueSummaryPage, accessed 7 September 2017), the quality of imputation for other ethnicities and females might be less reliable, for example, when applied to the MDD-CONVERGE dataset.

Just as medications acting on different pathways might have synergistic therapeutic effects, we believe that it is beneficial to have different approaches to drug repositioning to complement each other. Of course, computational methods leveraging human genomic data are not the only means to drug discovery. We believe that a combination of a variety of approaches, including experimental and computational ones, is required to speed up drug repurposing and discoveries.

Our enrichment analyses support the application of GSA in drug repositioning in depression and anxiety. However, we stress that our repositioning results should be validated in further pre-clinical and clinical studies before translation to practice. All the findings should be regarded tentative until further validations. In addition, GSA analyses do not provide information on the direction of effects, as discussed previously. Measures of statistical significance also do not provide definitive evidence for the actual therapeutic effects of the repositioned drugs.

The current study provides a proof-of-concept of how human genomic data may assist drug discovery or repositioning in depression/AD. However, each candidate has to be individually and carefully assessed (e.g. for their side-effects, efficacy, interaction with other agents, ability to cross blood–brain barrier, etc.) before being tested in clinical trials or applied in the clinic. For the above reasons, we do not believe every top ‘hit’ will be useful clinically, although in the future one may develop a more ‘intelligent’ repositioning framework to factor in the above elements. Our described methodology may be considered as a general way to prioritize promising candidates for further research. Given the huge cost [up to ~US$2.6 billion (105)] in developing a new drug (the cost is lower but still considerable for repositioning), even a minute reduction in failure rate will represent very substantial savings in absolute terms.

In summary, we have performed a drug repositioning analyses on depression and AD, using a GSA approach considering five related GWAS studies. We showed that the repositioned drugs are in general enriched for known psychiatric medications or those considered in clinical trials. Remarkably, antipsychotics and antidepressants were ranked among the top even if we considered all level 3 ATC drug classes. Our findings lend further support to the usefulness of human genomic data in guiding drug development in psychiatry.

Supplementary material

The supplementary material for this article can be found at https://doi.org/10.1017/S0033291718003641.

Author ORCIDs

Hon-Cheong So 0000-0002-7102-833X.

Acknowledgements

This work is partially supported by the Lo Kwee Seong Biomedical Research Fund and a Direct Grant from the Chinese University of Hong Kong. We also thank Professor Stephen K.W. Tsui and the Hong Kong Bioinformatics Centre for computing support. We would also like to acknowledge the Psychiatric Genomics Consortium, the CONVERGE Consortium, the Social Science Genetics Association Consortium and Otawa et al. for providing open access to full GWAS summary results.

Author contributions

H.-C.S. conceived and designed the study. H.-C.S. performed data analyses with input from C.K.L.C and K.Z. H.-C.S. interpreted the data. A.L. and S-Y.W. performed drug annotations. H.-C.S. wrote the manuscript and supervised the study.

Conflict of interest

None.

References

Al-Harbi, KS (2012) Treatment-resistant depression: therapeutic trends, challenges, and future directions. Patient Preference and Adherence 6, 369388.Google Scholar
Amador, M and Dani, JA (1991) Protein kinase inhibitor, H-7, directly affects N-methyl-D-aspartate receptor channels. Neuroscience Letters 124, 251255.Google Scholar
Anighoro, A, Bajorath, J and Rastelli, G (2014) Polypharmacology: challenges and opportunities in drug discovery. Journal of Medicinal Chemistry 57, 78747887.Google Scholar
Antonaci, F, Nappi, G, Galli, F, Manzoni, GC, Calabresi, P and Costa, A (2011) Migraine and psychiatric comorbidity: a review of clinical findings. Journal of Headache and Pain 12, 115125.Google Scholar
Artigas, F, Celada, P, Laruelle, M and Adell, A (2001) How does pindolol improve antidepressant action? Trends in Pharmacological Sciences 22, 224228.Google Scholar
Bakshi, A, Zhu, ZH, Vinkhuyzen, AAE, Hill, WD, Mcrae, AF, Visscher, PM and Yang, J (2016) Fast set-based association analysis using summary data from GWAS identifies novel gene loci for human complex traits. Scientific Reports 6, 32894.Google Scholar
Ballenger, JC (2000) Anxiety and depression: optimizing treatments. Primary Care Companion Journal Clinical Psychiatry 2, 7179.Google Scholar
Beaulieu, JM, Zhang, X, Rodriguiz, RM, Sotnikova, TD, Cools, MJ, Wetsel, WC, Gainetdinov, RR and Caron, MG (2008) Role of GSK3 beta in behavioral abnormalities induced by serotonin deficiency. Proceedings of the National Academy of Sciences of the USA 105, 13331338.Google Scholar
Benjamini, Y and Hochberg, Y (1995) Controlling the false discovery rate – a practical and powerful approach to multiple testing. Journal of the Royal Statistical Society Series B-Methodological 57, 289300.Google Scholar
Bigdeli, TB, Ripke, S, Peterson, RE, Trzaskowski, M, Bacanu, SA, Abdellaoui, A, Andlauer, TF, Beekman, AT, Berger, K, Blackwood, DH, Boomsma, DI, Breen, G, Buttenschon, HN, Byrne, EM, Cichon, S, Clarke, TK, Couvy-Duchesne, B, Craddock, N, De Geus, EJ, Degenhardt, F, Dunn, EC, Edwards, AC, Fanous, AH, Forstner, AJ, Frank, J, Gill, M, Gordon, SD, Grabe, HJ, Hamilton, SP, Hardiman, O, Hayward, C, Heath, AC, Henders, AK, Herms, S, Hickie, IB, Hoffmann, P, Homuth, G, Hottenga, JJ, Ising, M, Jansen, R, Kloiber, S, Knowles, JA, Lang, M, Li, QS, Lucae, S, Macintyre, DJ, Madden, PA, Martin, NG, Mcgrath, PJ, Mcguffin, P, Mcintosh, AM, Medland, SE, Mehta, D, Middeldorp, CM, Milaneschi, Y, Montgomery, GW, Mors, O, Muller-Myhsok, B, Nauck, M, Nyholt, DR, Nothen, MM, Owen, MJ, Penninx, BW, Pergadia, ML, Perlis, RH, Peyrot, WJ, Porteous, DJ, Potash, JB, Rice, JP, Rietschel, M, Riley, BP, Rivera, M, Schoevers, R, Schulze, TG, Shi, J, Shyn, SI, Smit, JH, Smoller, JW, Streit, F, Strohmaier, J, Teumer, A, Treutlein, J, Van Der Auwera, S, Van Grootheest, G, Van Hemert, AM, Volzke, H, Webb, BT, Weissman, MM, Wellmann, J, Willemsen, G, Witt, SH, Levinson, DF, Lewis, CM, Wray, NR, Flint, J, Sullivan, PF and Kendler, KS (2017) Genetic effects influencing risk for major depressive disorder in China and Europe. Translational Psychiatry 7, e1074.Google Scholar
Blier, P and Bergeron, R (1995) Effectiveness of pindolol with selected antidepressant drugs in the treatment of major depression. Journal of Clinical Psychopharmacology 15, 217222.Google Scholar
Brown, MB (1975) 400: a method for combining non-independent, one-sided tests of significance. Biometrics 31, 987992.Google Scholar
Cai, N, Bigdeli, TB, Kretzschmar, W, Li, YH, Liang, JQ, Song, L, Hu, JC, Li, QB, Jin, W, Hu, ZF, Wang, GB, Wang, LM, Qian, PY, Liu, Y, Jiang, T, Lu, Y, Zhang, XQ, Yin, Y, Li, YR, Xu, X, Gao, JF, Reimers, M, Webb, T, Riley, B, Bacanu, S, Peterson, RE, Chen, YP, Zhong, H, Liu, ZR, Wang, G, Sun, J, Sang, H, Jiang, GQ, Zhou, XY, Li, Y, Li, Y, Zhang, W, Wang, XY, Fang, X, Pan, RD, Miao, GD, Zhang, QW, Hu, J, Yu, FY, Du, B, Sang, WH, Li, KQ, Chen, GB, Cai, M, Yang, LJ, Yang, DL, Ha, BW, Hong, XH, Deng, H, Li, GY, Li, K, Song, Y, Gao, SG, Zhang, JB, Gan, ZY, Meng, HQ, Pan, JY, Gao, CG, Zhang, KR, Sun, N, Li, YH, Niu, QH, Zhang, YT, Liu, TQ, Hu, CM, Zhang, Z, Lv, LX, Dong, JC, Wang, XP, Tao, M, Wang, XM, Xia, J, Rong, H, He, Q, Liu, TB, Huang, GP, Mei, QY, Shen, ZM, Liu, Y, Shen, JH, Tian, T, Liu, XJ, Wu, WY, Gu, DH, Fu, GY, Shi, JG, Chen, YC, Gan, XC, Liu, LF, Wang, LN, Yang, FZ, Cong, EZ, Marchini, J, Yang, HM, Wang, J, et al. (2015) Sparse whole-genome sequencing identifies two loci for major depressive disorder. Nature 523, 588– +.Google Scholar
Ceulemans, DL, Hoppenbrouwers, ML, Gelders, YG and Reyntjens, AJ (1985) The influence of ritanserin, a serotonin antagonist, in anxiety disorders: a double-blind placebo-controlled study versus lorazepam. Pharmacopsychiatry 18, 303305.Google Scholar
Chen, YC, Wang, HN, Zhang, RG, Wang, HH, Peng, ZW, Sun, RZ and Tan, QR (2012) Microinjection of sanguinarine into the ventrolateral orbital cortex inhibits Mkp-1 and exerts an antidepressant-like effect in rats. Neuroscience Letters 506, 327331.Google Scholar
Cipriani, A, Saunders, K, Attenburrow, MJ, Stefaniak, J, Panchal, P, Stockton, S, Lane, TA, Tunbridge, EM, Geddes, JR and Harrison, PJ (2016) A systematic review of calcium channel antagonists in bipolar disorder and some considerations for their future development. Molecular Psychiatry 21, 13241332.Google Scholar
Corp, SA, Gitlin, MJ and Altshuler, LL (2014) A review of the use of stimulants and stimulant alternatives in treating bipolar depression and major depressive disorder. Journal of Clinical Psychiatry 75, 10101018.Google Scholar
Cross-Disorder Group of the Psychiatric Genomics C, Lee, SH, Ripke, S, Neale, BM, Faraone, SV, Purcell, SM, Perlis, RH, Mowry, BJ, Thapar, A, Goddard, ME, Witte, JS, Absher, D, Agartz, I, Akil, H, Amin, F, Andreassen, OA, Anjorin, A, Anney, R, Anttila, V, Arking, DE, Asherson, P, Azevedo, MH, Backlund, L, Badner, JA, Bailey, AJ, Banaschewski, T, Barchas, JD, Barnes, MR, Barrett, TB, Bass, N, Battaglia, A, Bauer, M, Bayes, M, Bellivier, F, Bergen, SE, Berrettini, W, Betancur, C, Bettecken, T, Biederman, J, Binder, EB, Black, DW, Blackwood, DH, Bloss, CS, Boehnke, M, Boomsma, DI, Breen, G, Breuer, R, Bruggeman, R, Cormican, P, Buccola, NG, Buitelaar, JK, Bunney, WE, Buxbaum, JD, Byerley, WF, Byrne, EM, Caesar, S, Cahn, W, Cantor, RM, Casas, M, Chakravarti, A, Chambert, K, Choudhury, K, Cichon, S, Cloninger, CR, Collier, DA, Cook, EH, Coon, H, Cormand, B, Corvin, A, Coryell, WH, Craig, DW, Craig, IW, Crosbie, J, Cuccaro, ML, Curtis, D, Czamara, D, Datta, S, Dawson, G, Day, R, De Geus, EJ, Degenhardt, F, Djurovic, S, Donohoe, GJ, Doyle, AE, Duan, J, Dudbridge, F, Duketis, E, Ebstein, RP, Edenberg, HJ, Elia, J, Ennis, S, Etain, B, Fanous, A, Farmer, AE, Ferrier, IN, Flickinger, M, Fombonne, E, Foroud, T, Frank, J, Franke, B, et al. (2013) Genetic relationship between five psychiatric disorders estimated from genome-wide SNPs. Nature Genetics 45, 984994.Google Scholar
De Jong, S, Vidler, LR, Mokrab, Y, Collier, DA and Breen, G (2016) Gene-set analysis based on the pharmacological profiles of drugs to identify repurposing opportunities in schizophrenia. Journal of Psychopharmacology 30, 826830.Google Scholar
De Leeuw, CA, Mooij, JM, Heskes, T and Posthuma, D (2015) MAGMA: generalized gene-set analysis of GWAS data. PLoS Computational Biology 11, e1004219.Google Scholar
De Leeuw, CA, Neale, BM, Heskes, T and Posthuma, D (2016) The statistical properties of gene-set analysis. Nature Reviews Genetics 17, 353364.Google Scholar
De Oliveira, MR, Chenet, AL, Duarte, AR, Scaini, G and Quevedo, J (2017) Molecular mechanisms underlying the anti-depressant effects of resveratrol: a review. Molecular Neurobiology 55, 45434559.Google Scholar
Dell'Osso, B, Ketter, TA, Cremaschi, L, Spagnolin, G and Altamura, AC (2013) Assessing the roles of stimulants/stimulant-like drugs and dopamine-agonists in the treatment of bipolar depression. Current Psychiatry Reports 15, 378.Google Scholar
Dhingra, D and Chhillar, R (2012) Antidepressant-like activity of ellagic acid in unstressed and acute immobilization-induced stressed mice. Pharmacological Reports 64, 796807.Google Scholar
Dubovsky, SL, Buzan, R, Thomas, M, Kassner, C and Cullum, CM (2001) Nicardipine improves the antidepressant action of ECT but does not improve cognition. Journal of ECT 17, 310.Google Scholar
Fedorova, I, Hashimoto, A, Fecik, RA, Hedrick, MP, Hanus, LO, Boger, DL, Rice, KC and Basile, AS (2001) Behavioral evidence for the interaction of oleamide with multiple neurotransmitter systems. Jouranl of Pharmacology and Expertimental Therapeutics 299, 332342.Google Scholar
Flint, J, Chen, Y, Shi, S, Kendler, KS and Consortium, C (2012) Epilogue: lessons from the CONVERGE study of major depressive disorder in China. Journal of Affective Disorders 140, 15.Google Scholar
Florez, JC, Jablonski, KA, Mcateer, JB, Franks, PW, Mason, CC, Mather, K, Horton, E, Goldberg, R, Dabelea, D, Kahn, SE, Arakaki, RF, Shuldiner, AR, Knowler, WC and Diabetes Prevention Program Research G (2012) Effects of genetic variants previously associated with fasting glucose and insulin in the Diabetes Prevention Program. PLoS ONE, 7, e44424.Google Scholar
Fournier, JC, Derubeis, RJ, Hollon, SD, Dimidjian, S, Amsterdam, JD, Shelton, RC and Fawcett, J (2010) Antidepressant drug effects and depression severity: a patient-level meta-analysis. JAMA 303, 4753.Google Scholar
Frances, A, Manning, D, Marin, D, Kocsis, J, Mckinney, K, Hall, W and Kline, M (1992) Relationship of anxiety and depression. Psychopharmacology 106, S82S86.Google Scholar
Freland, L and Beaulieu, JM (2012) Inhibition of GSK3 by lithium, from single molecules to signaling networks. Frontiers in Molecular Neuroscience 5, 14.Google Scholar
Fuller, RW, Snoddy, HD, Perry, KW, Roush, BW, Molloy, BB, Bymaster, FP and Wong, DT (1976) The effects of quipazine on serotonin metabolism in rat brain. Life Science 18, 925933.Google Scholar
Fung, M (2001) Evaluation of the characteristics of safety withdrawal of prescription drugs from worldwide pharmaceutical markets-1960 to 1999. Drug Information Journal 35, 293317.Google Scholar
Gale, CR, Hagenaars, SP, Davies, G, Hill, WD, Liewald, DCM, Cullen, B, Penninx, BW, Boomsma, DI, Pell, J, Mcintosh, AM, Smith, DJ, Deary, IJ, Harris, SE and Gwas, ICBP (2016) Pleiotropy between neuroticism and physical and mental health: findings from 108038 men and women in UK Biobank. Translational Psychiatry 6, e791.Google Scholar
Gao, K, Muzina, D, Gajwani, P and Calabrese, JR (2006) Efficacy of typical and atypical antipsychotics for primary and comorbid anxiety symptoms or disorders: a review. Journal of Clinical Psychiatry 67, 13271340.Google Scholar
Gaspar, HA and Breen, G (2017) Drug enrichment and discovery from schizophrenia genome-wide association results: an analysis and visualisation approach. Scientific Reports 7, 12460.Google Scholar
Gaulton, A, Bellis, LJ, Bento, AP, Chambers, J, Davies, M, Hersey, A, Light, Y, Mcglinchey, S, Michalovich, D, Al-Lazikani, B and Overington, JP (2012) ChEMBL: a large-scale bioactivity database for drug discovery. Nucleic Acids Research 40, D1100D1107.Google Scholar
Ghabrash, MF, Comai, S, Tabaka, J, Saint-Laurent, M, Booij, L and Gobbi, G (2016) Valproate augmentation in a subgroup of patients with treatment-resistant unipolar depression. World of Journal Biological Psychiatry 17, 165170.Google Scholar
Girish, C, Raj, V, Arya, J and Balakrishnan, S (2013) Involvement of the GABAergic system in the anxiolytic-like effect of the flavonoid ellagic acid in mice. European Journal of Pharmacology 710, 4958.Google Scholar
Gold, PW and Chrousos, GP (2002) Organization of the stress system and its dysregulation in melancholic and atypical depression: high vs low CRH/NE states. Molecular Psychiatry 7, 254275.Google Scholar
Goldberg, D (2011) The heterogeneity of ‘major depression’. World Psychiatry 10, 226228.Google Scholar
Goodwin, GM (2015) The overlap between anxiety, depression, and obsessive-compulsive disorder. Dialogues in Clinical Neuroscience 17, 249260.Google Scholar
Gulbins, E, Palmada, M, Reichel, M, Luth, A, Bohmer, C, Amato, D, Muller, CP, Tischbirek, CH, Groemer, TW, Tabatabai, G, Becker, KA, Tripal, P, Staedtler, S, Ackermann, TF, Van Brederode, J, Alzheimer, C, Weller, M, Lang, UE, Kleuser, B, Grassme, H and Kornhuber, J (2013) Acid sphingomyelinase-ceramide system mediates effects of antidepressant drugs. Nature Medicine 19, 934–+.Google Scholar
Harald, B and Gordon, P (2012) Meta-review of depressive subtyping models. Journal of Affective Disorders 139, 126140.Google Scholar
Huang, L, Li, F, Sheng, J, Xia, X, Ma, J, Zhan, M and Wong, ST (2014) Drugcomboranker: drug combination discovery based on target network analysis. Bioinformatics (Oxford, England) 30, i228i236.Google Scholar
Iyengar, RL, Gandhi, S, Aneja, A, Thorpe, K, Razzouk, L, Greenberg, J, Mosovich, S and Farkouh, ME (2013) NSAIDs are associated with lower depression scores in patients with osteoarthritis. American Journal of Medicine 126, 1017. e11–8.Google Scholar
Jardine, R, Martin, N, Henderson, A and Rao, D (1984) Genetic covariation between neuroticism and the symptoms of anxiety and depression. Genetic Epidemiology 1, 89107.Google Scholar
Jue, SG, Dawson, GW and Brogden, RN (1982) Amoxapine – a review of its pharmacology and efficacy in depressed states. Drugs 24, 123.Google Scholar
Kanba, S (2004) Although antidepressants and anxiolytics are frequently used together to treat depression in the acute phase, how effective is the concomitant use of these drugs? Journal of Psychiatry and Neuroscience 29, 485486.Google Scholar
Kessler, RC, Sampson, NA, Berglund, P, Gruber, MJ, Al-Hamzawi, A, Andrade, L, Bunting, B, Demyttenaere, K, Florescu, S, De Girolamo, G, Gureje, O, He, Y, Hu, C, Huang, Y, Karam, E, Kovess-Masfety, V, Lee, S, Levinson, D, Medina Mora, ME, Moskalewicz, J, Nakamura, Y, Navarro-Mateu, F, Browne, MA, Piazza, M, Posada-Villa, J, Slade, T, Ten Have, M, Torres, Y, Vilagut, G, Xavier, M, Zarkov, Z, Shahly, V and Wilcox, MA (2015) Anxious and non-anxious major depressive disorder in the World Health Organization World Mental Health Surveys. Epidemiology and Psychiatry Sciences 24, 210226.Google Scholar
Kim, S, Thiessen, PA, Bolton, EE, Chen, J, Fu, G, Gindulyte, A, Han, LY, He, JE, He, SQ, Shoemaker, BA, Wang, JY, Yu, B, Zhang, J and Bryant, SH (2016) Pubchem substance and compound databases. Nucleic Acids Research 44, D1202D1213.Google Scholar
Kirsch, I, Deacon, BJ, Huedo-Medina, TB, Scoboria, A, Moore, TJ and Johnson, BT (2008) Initial severity and antidepressant benefits: a meta-analysis of data submitted to the Food and Drug Administration. PLoS Medicine 5, e45.Google Scholar
Koch, HJ and Jurgens, TP (2009) Antidepressants in long-term migraine prevention. Drugs 69, 119.Google Scholar
Kohler, O, Gasse, C, Petersen, L, Ingstrup, KG, Nierenberg, AA, Mors, O and Ostergaard, SD (2016) The effect of concomitant treatment with SSRIs and statins: a population-based study. American Journal of Psychiatry 173, 807815.Google Scholar
Lahey, BB (2009) Public health significance of neuroticism. American Psychologist 64, 241256.Google Scholar
Lalonde, CD and Van Lieshout, RJ (2011) Treating generalized anxiety disorder with second generation antipsychotics: a systematic review and meta-analysis. Journal of Clinical Psychopharmacology 31, 326333.Google Scholar
Lamb, J, Crawford, ED, Peck, D, Modell, JW, Blat, IC, Wrobel, MJ, Lerner, J, Brunet, JP, Subramanian, A, Ross, KN, Reich, M, Hieronymus, H, Wei, G, Armstrong, SA, Haggarty, SJ, Clemons, PA, Wei, R, Carr, SA, Lander, ES and Golub, TR (2006) The connectivity map: using gene-expression signatures to connect small molecules, genes, and disease. Science 313, 19291935.Google Scholar
Lamers, F, Van Oppen, P, Comijs, HC, Smit, JH, Spinhoven, P, Van Balkom, AJ, Nolen, WA, Zitman, FG, Beekman, AT and Penninx, BW (2011) Comorbidity patterns of anxiety and depressive disorders in a large cohort study: the Netherlands Study of Depression and Anxiety (NESDA). Journal of Clinical Psychiatry 72, 341348.Google Scholar
Lamers, F, Vogelzangs, N, Merikangas, KR, De Jonge, P, Beekman, ATF and Penninx, BWJH (2013) Evidence for a differential role of HPA-axis function, inflammation and metabolic syndrome in melancholic versus atypical depression. Molecular Psychiatry 18, 692699.Google Scholar
Leost, M, Schultz, C, Link, A, Wu, YZ, Biernat, J, Mandelkow, EM, Bibb, JA, Snyder, GL, Greengard, P, Zaharevitz, DW, Gussio, R, Senderowicz, AM, Sausville, EA, Kunick, C and Meijer, L (2000) Paullones are potent inhibitors of glycogen synthase kinase-3 beta and cyclin-dependent kinase 5/p25. European Journal of Biochemistry 267, 59835994.Google Scholar
Lesser, IM, Myers, HF, Lin, KM, Bingham Mira, C, Joseph, NT, Olmos, NT, Schettino, J and Poland, RE (2010) Ethnic differences in antidepressant response: a prospective multi-site clinical trial. Depression and Anxiety 27, 5662.Google Scholar
Lesser, IM, Zisook, S, Gaynes, BN, Wisniewski, SR, Luther, JF, Fava, M, Khan, A, Mcgrath, P, Warden, D, Rush, AJ and Trivedi, M (2011) Effects of race and ethnicity on depression treatment outcomes: the CO-MED trial. Psychiatry Service 62, 11671179.Google Scholar
Li, YH, Yu, CY, Li, XX, Zhang, P, Tang, J, Yang, Q, Fu, T, Zhang, X, Cui, X, Tu, G, Zhang, Y, Li, S, Yang, F, Sun, Q, Qin, C, Zeng, X, Chen, Z, Chen, YZ and Zhu, F (2018) Therapeutic target database update 2018: enriched resource for facilitating bench-to-clinic research of targeted therapeutics. Nucleic Acids Research 46, D1121D1127.Google Scholar
Liu, JZ, Mcrae, AF, Nyholt, DR, Medland, SE, Wray, NR, Brown, KM, Hayward, NK, Montgomery, GW, Visscher, PM, Martin, NG, Macgregor, S and Investigators, A (2010) A versatile gene-based test for genome-wide association studies. American Journal of Human Genetics 87, 139145.Google Scholar
Maciejewski, H (2014) Gene set analysis methods: statistical models and methodological differences. Briefings in Bioinformatics 15, 504518.Google Scholar
Major Depressive Disorder Working Group of the Psychiatric GWAS Consortium, Ripke, S, Wray, NR, Lewis, CM, Hamilton, SP, Weissman, MM, Breen, G, Byrne, EM, Blackwood, DH, Boomsma, DI, Cichon, S, Heath, AC, Holsboer, F, Lucae, S, Madden, PA, Martin, NG, McGuffin, P, Muglia, P, Noethen, MM, Penninx, BP, Pergadia, ML, Potash, JB, Rietschel, M, Lin, D, Müller-Myhsok, B, Shi, J, Steinberg, S, Grabe, HJ, Lichtenstein, P, Magnusson, P, Perlis, RH, Preisig, M, Smoller, JW, Stefansson, K, Uher, R, Kutalik, Z, Tansey, KE, Teumer, A, Viktorin, A, Barnes, MR, Bettecken, T, Binder, EB, Breuer, R, Castro, VM, Churchill, SE, Coryell, WH, Craddock, N, Craig, IW, Czamara, D, De Geus, EJ, Degenhardt, F, Farmer, AE, Fava, M, Frank, J, Gainer, VS, Gallagher, PJ, Gordon, SD, Goryachev, S, Gross, M, Guipponi, M, Henders, AK, Herms, S, Hickie, IB, Hoefels, S, Hoogendijk, W, Hottenga, JJ, Iosifescu, DV, Ising, M, Jones, I, Jones, L, Jung-Ying, T, Knowles, JA, Kohane, IS, Kohli, MA, Korszun, A, Landen, M, Lawson, WB, Lewis, G, Macintyre, D, Maier, W, Mattheisen, M, McGrath, PJ, McIntosh, A, McLean, A, Middeldorp, CM, Middleton, L, Montgomery, GM, Murphy, SN, Nauck, M, Nolen, WA, Nyholt, DR, O'Donovan, M, Oskarsson, H, Pedersen, N, Scheftner, WA, Schulz, A, Schulze, TG, Shyn, SI, Sigurdsson, E, Slager, SL, Smit, JH, Stefansson, H, Steffens, M, Thorgeirsson, T, Tozzi, F, Treutlein, J, Uhr, M, van den Oord, EJ, Van Grootheest, G, Völzke, H, Weilburg, JB, Willemsen, G, Zitman, FG, Neale, B, Daly, M, Levinson, DF and Sullivan, PF (2012) A mega-analysis of genome-wide association studies for major depressive disorder. Molecular Psychiatry 18, 497511.Google Scholar
Malouff, JM, Thorsteinsson, EB and Schutte, NS (2005) The relationship between the five-factor model of personality and symptoms of clinical disorders: a meta-analysis. Journal of Psychopathology and Behavioral Assessment 27, 101114.Google Scholar
Mattingly, CJ, Colby, GT, Forrest, JN and Boyer, JL (2003) The Comparative Toxicogenomics Database (CTD). Environmental Health Perspectives 111, 793795.Google Scholar
Mcgrath, PJ, Khan, AY, Trivedi, MH, Stewart, JW, Morris, DW, Wisniewski, SR, Miyahara, S, Nierenberg, AA, Fava, M and Rush, AJ (2008) Response to a selective serotonin reuptake inhibitor (citalopram) in major depressive disorder with melancholic features: a STAR*D report. Journal of Clinical Psychiatry 69, 18471855.Google Scholar
Melartin, T, Leskela, U, Rytsala, H, Sokero, P, Lestela-Mielonen, P and Isometsa, E (2004) Co-morbidity and stability of melancholic features in DSM-IV major depressive disorder. Psychological Medicine 34, 14431452.Google Scholar
Merino, J, Leong, A, Posner, DC, Porneala, B, Masana, L, Dupuis, J and Florez, JC (2017) Genetically driven hyperglycemia increases risk of coronary artery disease separately from type 2 diabetes. Diabetes Care 40, 687693.Google Scholar
Murphy, E, Hou, L, Maher, BS, Woldehawariat, G, Kassem, L, Akula, N, Laje, G and Mcmahon, FJ (2013) Race, genetic ancestry and response to antidepressant treatment for major depression. Neuropsychopharmacology 38, 25982606.Google Scholar
Nappi, G, Sandrini, G, Granella, F, Ruiz, L, Cerutti, G, Facchinetti, F, Blandini, F and Manzoni, GC (1990) A New 5-Ht2 antagonist (Ritanserin) in the treatment of chronic headache with depression – a double-blind-study vs amitriptyline. Headache 30, 439444.Google Scholar
Nierenberg, AA, Dougherty, D and Rosenbaum, JF (1998) Dopaminergic agents and stimulants as antidepressant augmentation strategies. Journal of Clinical Psychiatry 59, 6064.Google Scholar
Okbay, A, Baselmans, BM, De Neve, JE, Turley, P, Nivard, MG, Fontana, MA, Meddens, SF, Linner, RK, Rietveld, CA, Derringer, J, Gratten, J, Lee, JJ, Liu, JZ, de Vlaming, R, Ahluwalia, TS, Buchwald, J, Cavadino, A, Frazier-Wood, AC, Furlotte, NA, Garfield, V, Geisel, MH, Gonzalez, JR, Haitjema, S, Karlsson, R, van der Laan, SW, Ladwig, KH, Lahti, J, van der Lee, SJ, Lind, PA, Liu, T, Matteson, L, Mihailov, E, Miller, MB, Minica, CC, Nolte, IM, Mook-Kanamori, D, van der Most, PJ, Oldmeadow, C, Qian, Y, Raitakari, O, Rawal, R, Realo, A, Rueedi, R, Schmidt, B, Smith, AV, Stergiakouli, E, Tanaka, T, Taylor, K, Thorleifsson, G, Wedenoja, J, Wellmann, J, Westra, HJ, Willems, SM, Zhao, W, LifeLines Cohort, S, Amin, N, Bakshi, A, Bergmann, S, Bjornsdottir, G, Boyle, PA, Cherney, S, Cox, SR, Davies, G, Davis, OS, Ding, J, Direk, N, Eibich, P, Emeny, RT, Fatemifar, G, Faul, JD, Ferrucci, L, Forstner, AJ, Gieger, C, Gupta, R, Harris, TB, Harris, JM, Holliday, EG, Hottenga, JJ, De Jager, PL, Kaakinen, MA, Kajantie, E, Karhunen, V, Kolcic, I, Kumari, M, Launer, LJ, Franke, L, Li-Gao, R, Liewald, DC, Koini, M, Loukola, A, Marques-Vidal, P, Montgomery, GW, Mosing, MA, Paternoster, L, Pattie, A, Petrovic, KE, Pulkki-Raback, L, Quaye, L, Raikkonen, K, Rudan, I, Scott, RJ, Smith, JA, Sutin, AR, Trzaskowski, M, Vinkhuyzen, AE, Yu, L, Zabaneh, D, Attia, JR, Bennett, DA, Berger, K, Bertram, L, Boomsma, DI, Snieder, H, Chang, SC, Cucca, F, Deary, IJ, van Duijn, CM, Eriksson, JG, Bultmann, U, de Geus, EJ, Groenen, PJ, Gudnason, V, Hansen, T, Hartman, CA, Haworth, CM, Hayward, C, Heath, AC, Hinds, DA, Hypponen, E, Iacono, WG, Jarvelin, MR, Jockel, KH, Kaprio, J, Kardia, SL, Keltikangas-Jarvinen, L, Kraft, P, Kubzansky, LD, Lehtimaki, T, Magnusson, PK, Martin, NG, McGue, M, Metspalu, A, Mills, M, de Mutsert, R, Oldehinkel, AJ, Pasterkamp, G, Pedersen, NL, Plomin, R, Polasek, O, Power, C, Rich, SS, Rosendaal, FR, den Ruijter, HM, Schlessinger, D, Schmidt, H, Svento, R, Schmidt, R, Alizadeh, BZ, Sorensen, TI, Spector, TD, Starr, JM, Stefansson, K, Steptoe, A, Terracciano, A, Thorsteinsdottir, U, Thurik, AR, Timpson, NJ, Tiemeier, H, Uitterlinden, AG, Vollenweider, P, Wagner, GG, Weir, DR, Yang, J, Conley, DC, Smith, GD, Hofman, A, Johannesson, M, Laibson, DI, Medland, SE, Meyer, MN, Pickrell, JK, Esko, T, Krueger, RF, Beauchamp, JP, Koellinger, PD, Benjamin, DJ, Bartels, M and Cesarini, D (2016) Genetic variants associated with subjective well-being, depressive symptoms, and neuroticism identified through genome-wide analyses. Nature Genetics 48, 624633.Google Scholar
Onasanwo, S, Faborode, S, Agrawal, M, Ijiwola, O, Jaiyesimi, B and Narender, T (2014) Antidepressant and anxiolytic potentials of chebulinic acid in laboratory rodent. Annals of Depression and Anxiety 1, 1032.Google Scholar
Ostacher, MJ, Iosifescu, DV, Hay, A, Blumenthal, SR, Sklar, P and Perlis, RH (2014) Pilot investigation of isradipine in the treatment of bipolar depression motivated by genome-wide association. Bipolar Disorders 16, 199203.Google Scholar
Otowa, T, Hek, K, Lee, M, Byrne, EM, Mirza, SS, Nivard, MG, Bigdeli, T, Aggen, SH, Adkins, D, Wolen, A, Fanous, A, Keller, MC, Castelao, E, Kutalik, Z, Der Auwera, SV, Homuth, G, Nauck, M, Teumer, A, Milaneschi, Y, Hottenga, JJ, Direk, N, Hofman, A, Uitterlinden, A, Mulder, CL, Henders, AK, Medland, SE, Gordon, S, Heath, AC, Madden, PA, Pergadia, ML, Van Der Most, PJ, Nolte, IM, Van Oort, FV, Hartman, CA, Oldehinkel, AJ, Preisig, M, Grabe, HJ, Middeldorp, CM, Penninx, BW, Boomsma, D, Martin, NG, Montgomery, G, Maher, BS, Van Den Oord, EJ, Wray, NR, Tiemeier, H and Hettema, JM (2016 a) Meta-analysis of genome-wide association studies of anxiety disorders. Molecular Psychiatry 21, 1485.Google Scholar
Otowa, T, Hek, K, Lee, M, Byrne, EM, Mirza, SS, Nivard, MG, Bigdeli, T, Aggen, SH, Adkins, D, Wolen, A, Fanous, A, Keller, MC, Castelao, E, Kutalik, Z, Der Auwera, SV, Homuth, G, Nauck, M, Teumer, A, Milaneschi, Y, Hottenga, JJ, Direk, N, Hofman, A, Uitterlinden, A, Mulder, CL, Henders, AK, Medland, SE, Gordon, S, Heath, AC, Madden, PAF, Pergadia, ML, Van Der Most, PJ, Nolte, IM, Van Oort, FVA, Hartman, CA, Oldehinkel, AJ, Preisig, M, Grabe, HJ, Middeldorp, CM, Penninx, BWJH, Boomsma, D, Martin, NG, Montgomery, G, Maher, BS, Van Den Oord, EJ, Wray, NR, Tiemeier, H and Hettema, JM (2016 b) Meta-analysis of genome-wide association studies of anxiety disorders. Molecular Psychiatry 21, 13911399.Google Scholar
Pande, AC, Crockatt, JG, Feltner, DE, Janney, CA, Smith, WT, Weisler, R, Londborg, PD, Bielski, RJ, Zimbroff, DL, Davidson, JRT and Liu-Dumaw, M (2003) Pregabalin in generalized anxiety disorder: a placebo-controlled trial. American Journal of Psychiatry 160, 533540.Google Scholar
Papassotiropoulos, A and De Quervain, DJF (2015) Failed drug discovery in psychiatry: time for human genome-guided solutions. Trends in Cognitive Sciences 19, 183187.Google Scholar
Parker, G, Fink, M, Shorter, E, Taylor, MA, Akiskal, H, Berrios, G, Bolwig, T, Brown, WA, Carroll, B, Healy, D, Klein, DF, Koukopoulos, A, Michels, R, Paris, J, Rubin, RT, Spitzer, R and Swartz, C (2010) Issues for DSM-5: whither melancholia? The case for its classification as a distinct mood disorder. American Journal of Psychiatry 167, 745747.Google Scholar
Peters, MF, Zacco, A, Gordon, J, Maciag, CM, Litwin, LC, Thompson, C, Schroeder, P, Sygowski, LA, Piser, TM and Brugel, TA (2011) Identification of short-acting kappa-opioid receptor antagonists with anxiolytic-like activity. European Journal of Pharmacology 661, 2734.Google Scholar
Philip, NS, Carpenter, LL, Tyrka, AR and Price, LH (2012) The nicotinic acetylcholine receptor as a target for antidepressant drug development. Scientific World Journal 2012, 104105.Google Scholar
Poole, W, Gibbs, DL, Shmulevich, I, Bernard, B and Knijnenburg, TA (2016) Combining dependent p values with an empirical adaptation of Brown's method. Bioinformatics (Oxford, England) 32, 430436.Google Scholar
Purcell, S, Neale, B, Todd-Brown, K, Thomas, L, Ferreira, MAR, Bender, D, Maller, J, Sklar, P, De Bakker, PIW, Daly, MJ and Sham, PC (2007) PLINK: a tool set for whole-genome association and population-based linkage analyses. American Journal of Human Genetics 81, 559575.Google Scholar
Rickels, K, Pollack, MH, Feltner, DE, Lydiard, RB, Zimbroff, DL, Bielski, RJ, Tobias, K, Brock, JD, Zornberg, GL and Pande, AC (2005) Pregabalin for treatment of generalized anxiety disorder – a 4-week, multicenter, double-blind, placebo-controlled trial of pregabalin and alprazolam. Archives of General Psychiatry 62, 10221030.Google Scholar
Ruderfer, DM, Charney, AW, Readhead, B, Kidd, BA, Kahler, AK, Kenny, PJ, Keiser, MJ, Moran, JL, Hultman, CM, Scott, SA, Sullivan, PF, Purcell, SM, Dudley, JT and Sklar, P (2016) Polygenic overlap between schizophrenia risk and antipsychotic response: a genomic medicine approach. The Lancet. Psychiatry 3, 350357.Google Scholar
Salagre, E, Fernandes, BS, Dodd, S, Brownstein, DJ and Berk, M (2016) Statins for the treatment of depression: a meta-analysis of randomized, double-blind, placebo-controlled trials. Journal of Affective Disorders 200, 235242.Google Scholar
Samuel, M, Zimovetz, EA, Gabriel, Z and Beard, SM (2011) Efficacy and safety of treatments for refractory generalized anxiety disorder: a systematic review. International Clinical Psychopharmacology 26, 6368.Google Scholar
Sanseau, P, Agarwal, P, Barnes, MR, Pastinen, T, Richards, JB, Cardon, LR and Mooser, V (2012) Use of genome-wide association studies for drug repositioning. Nature Biotechnology 30, 317320.Google Scholar
Sarkar, SK and Chang, C-K (1997) The Simes method for multiple hypothesis testing with positively dependent test statistics. Journal of the American Statistical Association 92, 16011608.Google Scholar
Simes, RJ (1986) An improved Bonferroni procedure for multiple tests of significance. Biometrika 73, 751754.Google Scholar
Smith, LA, Cornelius, VR, Azorin, JM, Perugi, G, Vieta, E, Young, AH and Bowden, CL (2010) Valproate for the treatment of acute bipolar depression: systematic review and meta-analysis. Journal of Affective Disorders 122, 19.Google Scholar
Smith, DJ, Escott-Price, V, Davies, G, Bailey, MES, Colodro-Conde, L, Ward, J, Vedernikov, A, Marioni, R, Cullen, B, Lyall, D, Hagenaars, SP, Liewald, DCM, Luciano, M, Gale, CR, Ritchie, SJ, Hayward, C, Nicholl, B, Bulik-Sullivan, B, Adams, M, Couvy-Duchesne, B, Graham, N, Mackay, D, Evans, J, Smith, BH, Porteous, DJ, Medland, SE, Martin, NG, Holmans, P, Mcintosh, AM, Pell, JP, Deary, IJ and O'Donovan, MC (2016) Genome-wide analysis of over 106 000 individuals identifies 9 neuroticism-associated loci. Molecular Psychiatry 21, 749757.Google Scholar
So, HC and Sham, PC (2011) Multiple testing and power calculations in genetic association studies. Cold Spring Harbor Protocols 2011, pdb.top95.Google Scholar
So, HC, Chau, CK, Chiu, WT, Ho, KS, Lo, CP, Yim, SH and Sham, PC (2017) Analysis of genome-wide association data highlights candidates for drug repositioning in psychiatry. Nature Neuroscience 20, 13421349.Google Scholar
Spielmans, GI, Berman, MI, Linardatos, E, Rosenlicht, NZ, Perry, A and Tsai, AC (2013) Adjunctive atypical antipsychotic treatment for major depressive disorder: a meta-analysis of depression, quality of life, and safety outcomes. PLoS Medicine 10, e1001403.Google Scholar
Steenen, SA, Van Wijk, AJ, Van Der Heijden, GJ, Van Westrhenen, R, De Lange, J and De Jongh, A (2016) Propranolol for the treatment of anxiety disorders: systematic review and meta-analysis. Journal of Psychopharmacology 30, 128139.Google Scholar
Stetler, C and Miller, GE (2011) Depression and hypothalamic-pituitary-adrenal activation: a quantitative summary of four decades of research. Psychosomatic Medicine 73, 114126.Google Scholar
Strauss, WH and Klieser, E (1991) Psychotropic effects of ritanserin, a selective S2 antagonist: an open study. European Neuropsychopharmacology 1, 101105.Google Scholar
Szmulewicz, AG, Angriman, F, Samame, C, Ferraris, A, Vigo, D and Strejilevich, SA (2017) Dopaminergic agents in the treatment of bipolar depression: a systematic review and meta-analysis. Acta Psychiatrica Scandinavica 135, 527538.Google Scholar
Tamburella, A, Micale, V, Navarria, A and Drago, F (2009) Antidepressant properties of the 5-HT4 receptor partial agonist, SL65.0155: behavioral and neurochemical studies in rats. Progress in Neuro-Psychopharmacology & Biological Psychiatry 33, 12051210.Google Scholar
Tang, TZ, Derubeis, RJ, Hollon, SD, Amsterdam, J, Shelton, R and Schalet, B (2009) Personality change during depression treatment a placebo-controlled trial. Archives of General Psychiatry 66, 13221330.Google Scholar
Touma, KTB, Zoucha, AM and Scarff, JR (2017) Liothyronine for depression: a review and guidance for safety monitoring. Innovations in Clinical Neuroscience 14, 2429.Google Scholar
Trivedi, MH, Rush, AJ, Wisniewski, SR, Nierenberg, AA, Warden, D, Ritz, L, Norquist, G, Howland, RH, Lebowitz, B, Mcgrath, PJ, Shores-Wilson, K, Biggs, MM, Balasubramani, GK, Fava, M and Team, SS (2006) Evaluation of outcomes with citalopram for depression using measurement-based care in STAR*D: implications for clinical practice. American Journal of Psychiatry 163, 2840.Google Scholar
Veena, J, Srikumar, BN, Mahati, K, Raju, TR and Shankaranarayana Rao, BS (2011) Oxotremorine treatment restores hippocampal neurogenesis and ameliorates depression-like behaviour in chronically stressed rats. Psychopharmacology (Berl) 217, 239253.Google Scholar
Wei, XY, Yang, JY, Dong, YX and Wu, CF (2007) Anxiolytic-like effects of oleamide in group-housed and socially isolated mice. Progress in Neuro-Psychopharmacology & Biological Psychiatry 31, 11891195.Google Scholar
Wei, WQ, Cronin, RM, Xu, H, Lasko, TA, Bastarache, L and Denny, JC (2013) Development and evaluation of an ensemble resource linking medications to their indications. Journal of American Medical Informatic Association 20, 954961.Google Scholar
World Health Organization (2017) Depression and Other Common Mental Disorders Global Health Estimates. Available at http://apps.who.int/iris/bitstream/10665/254610/1/WHO-MSD-MER-2017.2-eng.pdf?ua=1.Google Scholar
Wray, NR, Ripke, S, Mattheisen, M, Trzaskowski, M, Byrne, EM, Abdellaoui, A, Adams, MJ, Agerbo, E, Air, TM, Andlauer, TMF, Bacanu, SA, Baekvad-Hansen, M, Beekman, AFT, Bigdeli, TB, Binder, EB, Blackwood, DRH, Bryois, J, Buttenschon, HN, Bybjerg-Grauholm, J, Cai, N, Castelao, E, Christensen, JH, Clarke, TK, Coleman, JIR, Colodro-Conde, L, Couvy-Duchesne, B, Craddock, N, Crawford, GE, Crowley, CA, Dashti, HS, Davies, G, Deary, IJ, Degenhardt, F, Derks, EM, Direk, N, Dolan, CV, Dunn, EC, Eley, TC, Eriksson, N, Escott-Price, V, Kiadeh, FHF, Finucane, HK, Forstner, AJ, Frank, J, Gaspar, HA, Gill, M, Giusti-Rodriguez, P, Goes, FS, Gordon, SD, Grove, J, Hall, LS, Hannon, E, Hansen, CS, Hansen, TF, Herms, S, Hickie, IB, Hoffmann, P, Homuth, G, Horn, C, Hottenga, JJ, Hougaard, DM, Hu, M, Hyde, CL, Ising, M, Jansen, R, Jin, F, Jorgenson, E, Knowles, JA, Kohane, IS, Kraft, J, Kretzschmar, WW, Krogh, J, Kutalik, Z, Lane, JM, Li, Y, Li, Y, Lind, PA, Liu, X, Lu, L, MacIntyre, DJ, MacKinnon, DF, Maier, RM, Maier, W, Marchini, J, Mbarek, H, McGrath, P, McGuffin, P, Medland, SE, Mehta, D, Middeldorp, CM, Mihailov, E, Milaneschi, Y, Milani, L, Mill, J, Mondimore, FM, Montgomery, GW, Mostafavi, S, Mullins, N, Nauck, M, Ng, B, Nivard, MG, Nyholt, DR, O'Reilly, PF, Oskarsson, H, Owen, MJ, Painter, JN, Pedersen, CB, Pedersen, MG, Peterson, RE, Pettersson, E, Peyrot, WJ, Pistis, G, Posthuma, D, Purcell, SM, Quiroz, JA, Qvist, P, Rice, JP, Riley, BP, Rivera, M, Saeed Mirza, S, Saxena, R, Schoevers, R, Schulte, EC, Shen, L, Shi, J, Shyn, SI, Sigurdsson, E, Sinnamon, GBC, Smit, JH, Smith, DJ, Stefansson, H, Steinberg, S, Stockmeier, CA, Streit, F, Strohmaier, J, Tansey, KE, Teismann, H, Teumer, A, Thompson, W, Thomson, PA, Thorgeirsson, TE, Tian, C, Traylor, M, Treutlein, J, Trubetskoy, V, Uitterlinden, AG, Umbricht, D, Van der Auwera, S, van Hemert, AM, Viktorin, A, Visscher, PM, Wang, Y, Webb, BT, Weinsheimer, SM, Wellmann, J, Willemsen, G, Witt, SH, Wu, Y, Xi, HS, Yang, J, Zhang, F, eQtlgen, andMe, Arolt, V, Baune, BT, Berger, K, Boomsma, DI, Cichon, S, Dannlowski, U, de Geus, ECJ, DePaulo, JR, Domenici, E, Domschke, K, Esko, T, Grabe, HJ, Hamilton, SP, Hayward, C, Heath, AC, Hinds, DA, Kendler, KS, Kloiber, S, Lewis, G, Li, QS, Lucae, S, Madden, PFA, Magnusson, PK, Martin, NG, McIntosh, AM, Metspalu, A, Mors, O, Mortensen, PB, Muller-Myhsok, B, Nordentoft, M, Nothen, MM, O'Donovan, MC, Paciga, SA, Pedersen, NL, Penninx, B, Perlis, RH, Porteous, DJ, Potash, JB, Preisig, M, Rietschel, M, Schaefer, C, Schulze, TG, Smoller, JW, Stefansson, K, Tiemeier, H, Uher, R, Volzke, H, Weissman, MM, Werge, T, Winslow, AR, Lewis, CM, Levinson, DF, Breen, G, Borglum, AD, Sullivan, PF and Major Depressive Disorder Working Group of the Psychiatric Genomics C (2018) Genome-wide association analyses identify 44 risk variants and refine the genetic architecture of major depression. Nature Genetics 50, 668681.Google Scholar
Yadav, V, Chatterjee, SS, Majeed, M and Kumar, V (2015) Long lasting preventive effects of piperlongumine and a Piper longum extract against stress triggered pathologies in mice. Journal of Intercultural Ethnopharmacology 4, 277283.Google Scholar
Yi, LT, Liu, BB, Li, J, Luo, L, Liu, Q, Geng, D, Tang, Y, Xia, Y and Wu, D (2014) BDNF signaling is necessary for the antidepressant-like effect of naringenin. Progress in Neuro-Psychopharmacology & Biological Psychiatry 48, 135141.Google Scholar
Yoo, M, Shin, J, Kim, J, Ryall, KA, Lee, K, Lee, S, Jeon, M, Kang, J and Tan, AC (2015) DSigDB: drug signatures database for gene set analysis. Bioinformatics (Oxford, England) 31, 30693071.Google Scholar
You, H, Lu, W, Zhao, SP, Hu, ZP and Zhang, JN (2013) The relationship between statins and depression: a review of the literature. Expert Opinion on Pharmacotherapy 14, 14671476.Google Scholar
Zaminelli, T, Gradowski, RW, Bassani, TB, Barbiero, JK, Santiago, RM, Maria-Ferreira, D, Baggio, CH and Vital, MABF (2014) Antidepressant and antioxidative effect of ibuprofen in the rotenone model of Parkinson's disease. Neurotoxicity Research 26, 351362.Google Scholar
Zhou, XY, Keitner, GI, Qin, B, Ravindran, AV, Bauer, M, Del Giovane, C, Zhao, JP, Liu, YY, Fang, YR, Zhang, YQ and Xie, P (2015) Atypical antipsychotic augmentation for treatment-resistant depression: a systematic review and network meta-analysis. International Journal of Neuropsychopharmacology 18, pyv060.Google Scholar
Figure 0

Fig. 1. A schematic diagram showing our analysis workflow. We first extracted gene-sets associated with a variety of drugs (and chemicals), and examined whether the gene-set associated with each drug is enriched among the GWAS results. The drugs ranked among the top were considered potential candidates for repositioning, and a systemic literature search was performed for the top 20 repositioning hits of each psychiatric trait and meta-analyzed results. We also evaluated which drug classes were enriched. We specifically tested whether antidepressants/anxiolytics and other psychiatric medications were enriched, and then also examined all ATC (level 3) drug classes for enrichment in a hypothesis-free manner (ATC, Anatomical Therapeutic Chemical Classification System).

Figure 1

Table 1. Enrichment of repositioning hits derived from GWAS of major depressive disorder (MDD) and depressive symptoms for psychiatric drug categories (FDR-adjusted p values presented; gene-sets derived from entire DSigDB)

Figure 2

Table 2. Enrichment of repositioning hits derived from GWAS of anxiety disorders and neuroticism for psychiatric drug categories (FDR-adjusted p values presented; gene-sets derived from entire DSigDB)

Figure 3

Table 3. Enrichment of repositioning hits derived from meta-analyzed GWAS p values for psychiatric drug categories (FDR-adjusted p values presented; gene-sets derived from entire DSigDB)

Figure 4

Table 4. Top 5 enriched ATC drug classes (gene-sets derived from the whole DSigDB database)

Figure 5

Table 5. Top 5 enriched ATC drug classes with gene-sets derived from PubChem and ChEMBL only

Figure 6

Table 6. Selected repositioning hits (ranked within top 20 for a phenotype) with literature support

Supplementary material: File

So et al. supplementary material

So et al. supplementary material 1

Download So et al. supplementary material(File)
File 6.7 MB