Surprising conservation of schizophrenia risk genes in lower organisms reflects their essential function and the evolution of genetic liability

Evolutionary conservation of putative suicidality-related risk genes that produce diminished motivation corrected by clozapine, lithium and antidepressants

Background

Genome wide association studies (GWAS) and candidate gene analyses have identified genetic variants and genes that may increase the risk for suicidal thoughts and behaviors (STBs). Important unresolved issues surround these tentative risk variants such as the characteristics of the associated genes and how they might elicit STBs.

Methods

Putative suicidality-related risk genes (PSRGs) were identified by comprehensive literature search and were characterized with respect to evolutionary conservation, participation in gene interaction networks and associated phenotypes. Evolutionary conservation was established with database searches and BLASTP queries, whereas gene-gene interactions were ascertained with GeneMANIA. We then examined whether mutations in risk-gene counterparts in C. elegans produced a diminished motivation phenotype previously connected to suicide risk factors.

Results and conclusions

From the analysis, 105 risk-gene candidates were identified and found to be: 1) highly conserved during evolution, 2) enriched for essential genes, 3) involved in significant gene-gene interactions, and 4) associated with psychiatric disorders, metabolic disturbances and asthma/allergy. Evaluation of 17 mutant strains with loss-of-function/deletion mutations in PSRG orthologs revealed that 11 mutants showed significant evidence of diminished motivation that manifested as immobility in a foraging assay. Immobility was corrected in some or all of the mutants with clozapine, lithium and tricyclic antidepressant drugs. In addition, 5-HT2 receptor and muscarinic receptor antagonists restored goal-directed behavior in most or all of the mutants. These studies increase confidence in the validity of the PSRGs and provide initial clues about possible mechanisms that mediate STBs.

Identifying enhancer properties associated with genetic risk for complex traits using regulome-wide association studies

Genetic risk for complex traits is strongly enriched in non-coding genomic regions involved in gene regulation, especially enhancers. However, we lack adequate tools to connect the characteristics of these disruptions to genetic risk. Here, we propose RWAS (Regulome Wide Association Study), a new application of the MAGMA software package to identify the characteristics of enhancers that contribute to genetic risk for disease. RWAS involves three steps: (i) assign genotyped SNPs to cell type- or tissue-specific regulatory features (e.g., enhancers); (ii) test associations of each regulatory feature with a trait of interest for which genome-wide association study (GWAS) summary statistics are available; (iii) perform enhancer-set enrichment analyses to identify quantitative or categorical features of regulatory elements that are associated with the trait. These steps are implemented as a novel application of MAGMA, a tool originally developed for gene-based GWAS analyses. Applying RWAS to interrogate genetic risk for schizophrenia, we discovered a class of risk-associated AT-rich enhancers that are active in the developing brain and harbor binding sites for multiple transcription factors with neurodevelopmental functions. RWAS utilizes open-source software, and we provide a comprehensive collection of annotations for tissue-specific enhancer locations and features, including their evolutionary conservation, AT content, and co-localization with binding sites for hundreds of TFs. RWAS will enable researchers to characterize properties of regulatory elements associated with any trait of interest for which GWAS summary statistics are available.

Enhancers are regulatory regions that influence gene expression via the binding of transcription factors. Risk for many heritable diseases is enriched in regulatory regions, including enhancers. In this study, we introduce a novel application of the MAGMA software tool that enables testing for associations between enhancer attributes and risk, and we use this method to determine the enhancer characteristics that are associated with risk for schizophrenia. We found that enhancers associated with schizophrenia risk are both evolutionarily conserved and in physical contact with mutation-intolerant genes, many of which have neurodevelopmental functions. Risk-associated enhancers are also AT-rich and contain binding sites for neurodevelopmental transcription factors.

How Variation in Risk Allele Output and Gene Interactions Shape the Genetic Architecture of Schizophrenia

Schizophrenia is a highly heritable polygenic psychiatric disorder. Characterization of its genetic architecture may lead to a better understanding of the overall burden of risk variants and how they determine susceptibility to disease. A major goal of this project is to develop a modeling approach to compare and quantify the relative effects of single nucleotide polymorphisms (SNPs), copy number variants (CNVs) and other factors. We derived a mathematical model for the various genetic contributions based on the probability of expressing a combination of risk variants at a frequency that matched disease prevalence. The model included estimated risk variant allele outputs (VAOs) adjusted for population allele frequency. We hypothesized that schizophrenia risk genes would be more interactive than random genes and we confirmed this relationship. Gene–gene interactions may cause network ripple effects that spread and amplify small individual effects of risk variants. The modeling revealed that the number of risk alleles required to achieve the threshold for susceptibility will be determined by the average functional locus output (FLO) associated with a risk allele, the risk allele frequency (RAF), the number of protective variants present and the extent of gene interactions within and between risk loci. The model can account for the quantitative impact of protective variants as well as CNVs on disease susceptibility. The fact that non-affected individuals must carry a non-trivial burden of risk alleles suggests that genetic susceptibility will inevitably reach the threshold for schizophrenia at a recurring frequency in the population.

Protein Receptors Evolved from Homologous Cohesion Modules That Self-Associated and Are Encoded by Interactive Networked Genes

Previously, it was proposed that protein receptors evolved from self-binding peptides that were encoded by self-interacting gene segments (inverted repeats) widely dispersed in the genome. In addition, self-association of the peptides was thought to be mediated by regions of amino acid sequence similarity. To extend these ideas, special features of receptors have been explored, such as their degree of homology to other proteins, and the arrangement of their genes for clues about their evolutionary origins and dynamics in the genome. As predicted, BLASTP searches for homologous proteins detected a greater number of unique hits for queries with receptor sequences than for sequences of randomly-selected, non-receptor proteins. This suggested that the building blocks (cohesion modules) for receptors were duplicated, dispersed, and maintained in the genome, due to structure/function relationships discussed here. Furthermore, the genes coding for a representative panel of receptors participated in a larger number of gene-gene interactions than for randomly-selected genes. This could conceivably reflect a greater evolutionary conservation of the receptor genes, with their more extensive integration into networks along with inherent properties of the genes themselves. In support of the latter possibility, some receptor genes were located in active areas of adaptive gene relocation/amalgamation to form functional blocks of related genes. It is suggested that adaptive relocation might allow for their joint regulation by common promoters and enhancers, and affect local chromatin structural domains to facilitate or repress gene expression. Speculation is included about the nature of the coordinated communication between receptors and the genes that encode them.

Analysis of Major Depression Risk Genes Reveals Evolutionary Conservation, Shared Phenotypes, and Extensive Genetic Interactions

Major depressive disorder (MDD) affects around 15% of the population at some stage in their lifetime. It can be gravely disabling and it is associated with increased risk of suicide. Genetics play an important role; however, there are additional environmental contributions to the pathogenesis. A number of possible risk genes that increase liability for developing symptoms of MDD have been identified in genome-wide association studies (GWAS). The goal of this study was to characterize the MDD risk genes with respect to the degree of evolutionary conservation in simpler model organisms such as Caenorhabditis elegans and zebrafish, the phenotypes associated with variation in these genes and the extent of network connectivity. The MDD risk genes showed higher conservation in C. elegans and zebrafish than genome-to-genome comparisons. In addition, there were recurring themes among the phenotypes associated with variation of these risk genes in C. elegans. The phenotype analysis revealed enrichment for essential genes with pleiotropic effects. Moreover, the MDD risk genes participated in more interactions with each other than did randomly-selected genes from similar-sized gene sets. Syntenic blocks of risk genes with common functional activities were also identified. By characterizing evolutionarily-conserved counterparts to the MDD risk genes, we have gained new insights into pathogenetic processes relevant to the emergence of depressive symptoms in man.

Candidate risk genes for bipolar disorder are highly conserved during evolution and highly interconnected

OBJECTIVES

Bipolar disorder (BPD) is a highly heritable psychiatric disorder whose genetic complexity and pathogenetic mechanisms are still being unraveled. The main goal of this work was to characterize BPD risk-gene candidates (identified by Nurnberger et al., JAMA Psychiatry 71:657, 2014, and Stahl et al., Nat. Genet. 51:793, 2019) with respect to their evolutionary conservation, associated phenotypes, and extent of gene-gene interactions.

METHODS

Database searches and BLAST were used to identify homologous counterparts of human BPD risk genes in C. elegans, zebrafish, and Drosophila. Phenotypes associated with the C. elegans genes were annotated and searched. With GeneMANIA, we characterized and quantified gene-gene interactions among members of the BPD gene set in comparison to randomly chosen gene sets of the same size.

RESULTS

BPD risk genes are highly conserved across species and are enriched for essential genes and genes associated with lethality and altered life span. They are significantly more interactive with each other in comparison to random genes. We identified syntenic blocks of risk genes, which provided potential insights into molecular pathways and co-morbidities associated with BPD including coronary disease, obesity, and decreased life expectancy.

CONCLUSIONS

BPD risk genes appear to be special in terms of their degree of conservation, interconnectedness, and pleiotropic effects that extend beyond a role in brain function. Key hub genes or pleiotropic regulatory components may represent attractive targets for future drug discovery.

The Drosophila ortholog of the schizophrenia-associated CACNA1A and CACNA1B voltage-gated calcium channels regulate memory, sleep and circadian rhythms

Schizophrenia exhibits up to 80% heritability. A number of genome wide association studies (GWAS) have repeatedly shown common variants in voltage-gated calcium (Ca_v) channel genes CACNA1C, CACNA1I and CACNA1G have a major contribution to the risk of the disease. More recently, studies using whole exome sequencing have also found that CACNA1B (Ca_v2.2 N-type) deletions and rare disruptive variants in CACNA1A (Ca_v2.1 P/Q-type) are associated with schizophrenia. The negative symptoms of schizophrenia include behavioural defects such as impaired memory, sleep and circadian rhythms. It is not known how variants in schizophrenia-associated genes contribute to cognitive and behavioural symptoms, thus hampering the development of treatment for schizophrenia symptoms. In order to address this knowledge gap, we studied behavioural phenotypes in a number of loss of function mutants for the Drosophila ortholog of the Ca_v2 gene family called cacophony (cac). cac mutants showed several behavioural features including decreased night-time sleep and hyperactivity similar to those reported in human patients. The change in timing of sleep-wake cycles suggested disrupted circadian rhythms, with the loss of night-time sleep being caused by loss of cac just in the circadian clock neurons. These animals also showed a reduction in rhythmic circadian behaviour a phenotype that also could be mapped to the central clock. Furthermore, reduction of cac just in the clock resulted in a lengthening of the 24 h period. In order to understand how loss of Ca_v2 function may lead to cognitive deficits and underlying cellular pathophysiology we targeted loss of function of cac to the memory centre of the fly, called the mushroom bodies (MB). This manipulation was sufficient to cause reduction in both short- and intermediate-term associative memory. Memory impairment was accompanied by a decrease in Ca²⁺ transients in response to a depolarizing stimulus, imaged in the MB presynaptic terminals. This work shows loss of cac Ca_v2 channel function alone causes a number of cognitive and behavioural deficits and underlying reduced neuronal Ca²⁺ transients, establishing Drosophila as a high-throughput in vivo genetic model to study the Ca_v channel pathophysiology related to schizophrenia.

Genomic Chaos Begets Psychiatric Disorder

The processes that created the primordial genome are inextricably linked to current day vulnerability to developing a psychiatric disorder as summarized in this review article. Chaos and dynamic forces including duplication, transposition, and recombination generated the protogenome. To survive early stages of genome evolution, self-organization emerged to curb chaos. Eventually, the human genome evolved through a delicate balance of chaos/instability and organization/stability. However, recombination coldspots, silencing of transposable elements, and other measures to limit chaos also led to retention of variants that increase risk for disease. Moreover, ongoing dynamics in the genome creates various new mutations that determine liability for psychiatric disorders. Homologous recombination, long-range gene regulation, and gene interactions were all guided by spooky action-at-a-distance, which increased variability in the system. A probabilistic system of life was required to deal with a changing environment. This ensured the generation of outliers in the population, which enhanced the probability that some members would survive unfavorable environmental impacts. Some of the outliers produced through this process in man are ill suited to cope with the complex demands of modern life. Genomic chaos and mental distress from the psychological challenges of modern living will inevitably converge to produce psychiatric disorders in man.

Per- and Polyfluoroalkyl Substances (PFAS) Neurotoxicity in Sentinel and Non-Traditional Laboratory Model Systems: Potential Utility in Predicting Adverse Outcomes in Human Health

Per- and polyfluoroalkyl substances (PFAS) are a class of chemicals that were widely used in manufacturing and are now present in the environment throughout the world. It is known that various PFAS are quantifiable in human in blood, but potential adverse health outcomes remain unclear. Sentinel and non-traditional model species are useful to study potential toxicity of PFAS in order to understand the relationship between environmental and human health. Here, we present a critical review of studies on the neurotoxicity of PFAS in sentinel and non-traditional laboratory model systems, including Caenorhabditis elegans (nematode), Dugesia japonica (planarian), Rana pipiens (frogs), Danio rerio and Oryzias melastigma (fish), and Ursus maritimus (polar bears). PFAS have been implicated in developmental neurotoxicity in non-traditional and traditional model systems as well as sentinel species, including effects on neurotransmitter levels, especially acetylcholine and its metabolism. However, further research on the mechanisms of toxicity needs to be conducted to determine if these chemicals are affecting organisms in a similar manner. Overall, findings tend to be similar among the various species, but bioaccumulation may vary, which needs to be taken into account in future studies by quantifying target organ concentrations of PFAS to better compare different species. Furthermore, data on the majority of PFAS is lacking in neurotoxicity testing, and additional studies are needed to corroborate findings thus far.

A critical review of zebrafish schizophrenia models: Time for validation?

Schizophrenia is a mental disorder that affects 1% of the population worldwide and is manifested as a broad spectrum of symptoms, from hallucinations to memory impairment. It is believed that genetic and/or environmental factors may contribute to the occurrence of this disease. Recently, the zebrafish has emerged as a valuable and attractive model for various neurological disorders including schizophrenia. In this review, we describe current pharmacological models of schizophrenia with special emphasis on providing insights into the pros and cons of using zebrafish as a behavioural model of this disease. Moreover, we highlight the advantages and utility of using zebrafish for elucidating the genetic mechanisms underlying this psychiatric disorder. We believe that the zebrafish has high potential also in the area of precision medicine and may complement the development of therapeutics, especially for pharmacoresistant patients.