Identifying drug targets for schizophrenia through gene prioritization.

Background: Schizophrenia genome-wide association studies (GWASes) have identified >250 significant loci and prioritized >100 disease-related genes. However, gene prioritization efforts have mostly been restricted to locus-based methods that ignore information from the rest of the genome. Methods: To more accurately characterize genes involved in schizophrenia etiology, we applied a combination of highly-predictive tools to a published GWAS of 67,390 schizophrenia cases and 94,015 controls. We combined both locus-based methods (fine-mapped coding variants, distance to GWAS signals) and genome-wide methods (PoPS, MAGMA, ultra-rare coding variant burden tests). To validate our findings, we compared them with previous prioritization efforts, known neurodevelopmental genes, and results from the PsyOPS tool. Results: We prioritized 62 schizophrenia genes, 41 of which were also highlighted by our validation methods. In addition to DRD2 , the principal target of antipsychotics, we prioritized 9 genes that are targeted by approved or investigational drugs. These included drugs targeting glutamatergic receptors ( GRIN2A and GRM3 ), calcium channels ( CACNA1C and CACNB2 ), and GABA B receptor ( GABBR2 ). These also included genes in loci that are shared with an addiction GWAS ( e.g. PDE4B and VRK2 ). Conclusions: We curated a high-quality list of 62 genes that likely play a role in the development of schizophrenia. Developing or repurposing drugs that target these genes may lead to a new generation of schizophrenia therapies. Rodent models of addiction more closely resemble the human disorder than rodent models of schizophrenia. As such, genes prioritized for both disorders could be explored in rodent addiction models, potentially facilitating drug development.

Fine-mapping genomic loci refines bipolar disorder risk genes.

Bipolar disorder (BD) is a heritable mental illness with complex etiology. While the largest published genome-wide association study identified 64 BD risk loci, the causal SNPs and genes within these loci remain unknown. We applied a suite of statistical and functional fine-mapping methods to these loci, and prioritized 22 likely causal SNPs for BD. We mapped these SNPs to genes, and investigated their likely functional consequences by integrating variant annotations, brain cell-type epigenomic annotations, brain quantitative trait loci, and results from rare variant exome sequencing in BD. Convergent lines of evidence supported the roles of SCN2A, TRANK1, DCLK3, INSYN2B, SYNE1, THSD7A, CACNA1B, TUBBP5, PLCB3, PRDX5, KCNK4, AP001453.3, TRPT1, FKBP2, DNAJC4, RASGRP1, FURIN, FES, YWHAE, DPH1, GSDMB, MED24, THRA, EEF1A2, and KCNQ2 in BD. These represent promising candidates for functional experiments to understand biological mechanisms and therapeutic potential. Additionally, we demonstrated that fine-mapping effect sizes can improve performance and transferability of BD polygenic risk scores across ancestrally diverse populations, and present a high-throughput fine-mapping pipeline (https://github.com/mkoromina/SAFFARI).

Identifying genetic differences between bipolar disorder and major depression through multiple GWAS.

Background: Accurate diagnosis of bipolar disorder (BD) is difficult in clinical practice, with an average delay between symptom onset and diagnosis of about 7 years. A key reason is that the first manic episode is often preceded by a depressive one, making it difficult to distinguish BD from unipolar major depressive disorder (MDD). Aims: Here, we use genome-wide association analyses (GWAS) to identify differential genetic factors and to develop predictors based on polygenic risk scores that may aid early differential diagnosis. Methods: Based on individual genotypes from case-control cohorts of BD and MDD shared through the Psychiatric Genomics Consortium, we compile case-case-control cohorts, applying a careful merging and quality control procedure. In a resulting cohort of 51,149 individuals (15,532 BD cases, 12,920 MDD cases and 22,697 controls), we perform a variety of GWAS and polygenic risk scores (PRS) analyses. Results: While our GWAS is not well-powered to identify genome-wide significant loci, we find significant SNP-heritability and demonstrate the ability of the resulting PRS to distinguish BD from MDD, including BD cases with depressive onset. We replicate our PRS findings, but not signals of individual loci in an independent Danish cohort (iPSYCH 2015 case-cohort study, N=25,966). We observe strong genetic correlation between our case-case GWAS and that of case-control BD. Conclusions: We find that MDD and BD, including BD with a depressive onset, are genetically distinct. Further, our findings support the hypothesis that Controls - MDD - BD primarily lie on a continuum of genetic risk. Future studies with larger and richer samples will likely yield a better understanding of these findings and enable the development of better genetic predictors distinguishing BD and, importantly, BD with depressive onset from MDD.

Distinguishing different psychiatric disorders using DDx-PRS.

Despite great progress on methods for case-control polygenic prediction (e.g. schizophrenia vs. control), there remains an unmet need for a method that genetically distinguishes clinically related disorders (e.g. schizophrenia (SCZ) vs. bipolar disorder (BIP) vs. depression (MDD) vs. control); such a method could have important clinical value, especially at disorder onset when differential diagnosis can be challenging. Here, we introduce a method, Differential Diagnosis-Polygenic Risk Score (DDx-PRS), that jointly estimates posterior probabilities of each possible diagnostic category (e.g. SCZ=50%, BIP=25%, MDD=15%, control=10%) by modeling variance/covariance structure across disorders, leveraging case-control polygenic risk scores (PRS) for each disorder (computed using existing methods) and prior clinical probabilities for each diagnostic category. DDx-PRS uses only summary-level training data and does not use tuning data, facilitating implementation in clinical settings. In simulations, DDx-PRS was well-calibrated (whereas a simpler approach that analyzes each disorder marginally was poorly calibrated), and effective in distinguishing each diagnostic category vs. the rest. We then applied DDx-PRS to Psychiatric Genomics Consortium SCZ/BIP/MDD/control data, including summary-level training data from 3 case-control GWAS ( N =41,917-173,140 cases; total N =1,048,683) and held-out test data from different cohorts with equal numbers of each diagnostic category (total N =11,460). DDx-PRS was well-calibrated and well-powered relative to these training sample sizes, attaining AUCs of 0.66 for SCZ vs. rest, 0.64 for BIP vs. rest, 0.59 for MDD vs. rest, and 0.68 for control vs. rest. DDx-PRS produced comparable results to methods that leverage tuning data, confirming that DDx-PRS is an effective method. True diagnosis probabilities in top deciles of predicted diagnosis probabilities were considerably larger than prior baseline probabilities, particularly in projections to larger training sample sizes, implying considerable potential for clinical utility under certain circumstances. In conclusion, DDx-PRS is an effective method for distinguishing clinically related disorders.

Mapping genomic loci implicates genes and synaptic biology in schizophrenia.

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

Brain-Heart Link in Schizophrenia: Cognitive Inhibitory Control Deficit in Patients Is Specifically Related to Parasympathetic Dysregulation.

BACKGROUND: This study examined the connection between two prominent deficits in schizophrenia: the deficit in parasympathetic regulation and the deficit in cognitive inhibitory control, within the framework of the Neurovisceral Integration Model (NIM). STUDY DESIGN: Thirty healthy controls and 30 patients with schizophrenia performed the internationally standardized antisaccade protocol while their electrocardiographic data were recorded. The interaction between the group, the cognitive inhibitory control as measured with error rate (ER) in the antisaccade task and parasympathetic activity as measured with the High Frequency power component of Heart Rate Variability (HF-HRV) was tested. STUDY RESULTS: Findings confirmed that decreased HF-HRV was specifically related to increased ER in patients with schizophrenia. In contrast, patient deficits in other oculomotor function measures such as reaction time and reaction time variability related to volitional movement control and cognitive stability respectively were not linked to the deficit in parasympathetic regulation. CONCLUSIONS: Our study validates the theory behind NIM proposing that cognitive inhibition has common physiological substrate with parasympathetic regulation. Future research could test this brain-heart link in other mental disorders especially those with a prominent deficit in inhibitory cognitive function.

Genome-wide association study of more than 40,000 bipolar disorder cases provides new insights into the underlying biology.

Bipolar disorder is a heritable mental illness with complex etiology. We performed a genome-wide association study of 41,917 bipolar disorder cases and 371,549 controls of European ancestry, which identified 64 associated genomic loci. Bipolar disorder risk alleles were enriched in genes in synaptic signaling pathways and brain-expressed genes, particularly those with high specificity of expression in neurons of the prefrontal cortex and hippocampus. Significant signal enrichment was found in genes encoding targets of antipsychotics, calcium channel blockers, antiepileptics and anesthetics. Integrating expression quantitative trait locus data implicated 15 genes robustly linked to bipolar disorder via gene expression, encoding druggable targets such as HTR6, MCHR1, DCLK3 and FURIN. Analyses of bipolar disorder subtypes indicated high but imperfect genetic correlation between bipolar disorder type I and II and identified additional associated loci. Together, these results advance our understanding of the biological etiology of bipolar disorder, identify novel therapeutic leads and prioritize genes for functional follow-up studies.

Characterisation of age and polarity at onset in bipolar disorder.

BACKGROUND: Studying phenotypic and genetic characteristics of age at onset (AAO) and polarity at onset (PAO) in bipolar disorder can provide new insights into disease pathology and facilitate the development of screening tools. AIMS: To examine the genetic architecture of AAO and PAO and their association with bipolar disorder disease characteristics. METHOD: Genome-wide association studies (GWASs) and polygenic score (PGS) analyses of AAO (n = 12 977) and PAO (n = 6773) were conducted in patients with bipolar disorder from 34 cohorts and a replication sample (n = 2237). The association of onset with disease characteristics was investigated in two of these cohorts. RESULTS: Earlier AAO was associated with a higher probability of psychotic symptoms, suicidality, lower educational attainment, not living together and fewer episodes. Depressive onset correlated with suicidality and manic onset correlated with delusions and manic episodes. Systematic differences in AAO between cohorts and continents of origin were observed. This was also reflected in single-nucleotide variant-based heritability estimates, with higher heritabilities for stricter onset definitions. Increased PGS for autism spectrum disorder (ß = -0.34 years, s.e. = 0.08), major depression (ß = -0.34 years, s.e. = 0.08), schizophrenia (ß = -0.39 years, s.e. = 0.08), and educational attainment (ß = -0.31 years, s.e. = 0.08) were associated with an earlier AAO. The AAO GWAS identified one significant locus, but this finding did not replicate. Neither GWAS nor PGS analyses yielded significant associations with PAO. CONCLUSIONS: AAO and PAO are associated with indicators of bipolar disorder severity. Individuals with an earlier onset show an increased polygenic liability for a broad spectrum of psychiatric traits. Systematic differences in AAO across cohorts, continents and phenotype definitions introduce significant heterogeneity, affecting analyses.

RICOPILI: Rapid Imputation for COnsortias PIpeLIne.

SUMMARY: Genome-wide association study (GWAS) analyses, at sufficient sample sizes and power, have successfully revealed biological insights for several complex traits. RICOPILI, an open-sourced Perl-based pipeline was developed to address the challenges of rapidly processing large-scale multi-cohort GWAS studies including quality control (QC), imputation and downstream analyses. The pipeline is computationally efficient with portability to a wide range of high-performance computing environments. RICOPILI was created as the Psychiatric Genomics Consortium pipeline for GWAS and adopted by other users. The pipeline features (i) technical and genomic QC in case-control and trio cohorts, (ii) genome-wide phasing and imputation, (iv) association analysis, (v) meta-analysis, (vi) polygenic risk scoring and (vii) replication analysis. Notably, a major differentiator from other GWAS pipelines, RICOPILI leverages on automated parallelization and cluster job management approaches for rapid production of imputed genome-wide data. A comprehensive meta-analysis of simulated GWAS data has been incorporated demonstrating each step of the pipeline. This includes all the associated visualization plots, to allow ease of data interpretation and manuscript preparation. Simulated GWAS datasets are also packaged with the pipeline for user training tutorials and developer work. AVAILABILITY AND IMPLEMENTATION: RICOPILI has a flexible architecture to allow for ongoing development and incorporation of newer available algorithms and is adaptable to various HPC environments (QSUB, BSUB, SLURM and others). Specific links for genomic resources are either directly provided in this paper or via tutorials and external links. The central location hosting scripts and tutorials is found at this URL: https://sites.google.com/a/broadinstitute.org/RICOPILI/home. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Evaluating cognitive models of visual word recognition using fMRI: Effects of lexical and sublexical variables.

In this study predictions of the dual-route cascaded (DRC) model of word reading were tested using fMRI. Specifically, patterns of co-localization were investigated: (a) between pseudoword length effects and a pseudowords vs. fixation contrast, to reveal the sublexical grapho-phonemic conversion (GPC) system; and (b) between word frequency effects and a words vs. pseudowords contrast, to reveal the orthographic and phonological lexicon. Forty four native speakers of Greek were scanned at 3T in an event-related lexical decision task with three event types: (a) 150 words in which frequency, length, bigram and syllable frequency, neighborhood, and orthographic consistency were decorrelated; (b) 150 matched pseudowords; and (c) fixation. Whole-brain analysis failed to reveal the predicted co-localizations. Further analysis with participant-specific regions of interest defined within masks from the group contrasts revealed length effects in left inferior parietal cortex and frequency effects in the left middle temporal gyrus. These findings could be interpreted as partially consistent with the existence of the GPC system and phonological lexicon of the model, respectively. However, there was no evidence in support of an orthographic lexicon, weakening overall support for the model. The results are discussed with respect to the prospect of using neuroimaging in cognitive model evaluation.

Judging the judges' performance in rhythmic gymnastics.

INTRODUCTION: Rhythmic gymnastics (RG) is an aesthetic event balancing between art and sport that also has a performance rating system (Code of Points) given by the International Gymnastics Federation. It is one of the sports in which competition results greatly depend on the judges' evaluation. In the current study, we explored the judges' performance in a five-gymnast ensemble routine. METHODS: An expert-novice paradigm (10 international-level, 10 national-level, and 10 novice-level judges) was implemented under a fully simulated procedure of judgment in a five-gymnast ensemble routine of RG using two videos of routines performed by the Greek national team of RG. Simultaneous recordings of two-dimensional eye movements were taken during the judgment procedure to assess the percentage of time spent by each judge viewing the videos and fixation performance of each judge when an error in gymnast performance had occurred. RESULTS: All judge level groups had very modest performance of error recognition on gymnasts' routines, and the best international judges reported approximately 40% of true errors. Novice judges spent significantly more time viewing the videos compared with national and international judges and spent significantly more time fixating detected errors than the other two groups. National judges were the only group that made efficient use of fixation to detect errors. CONCLUSIONS: The fact that international-level judges outperformed both other groups, while not relying on visual fixation to detect errors, suggests that these experienced judges probably make use of other cognitive strategies, increasing their overall error detection efficiency, which was, however, still far below optimum.