The contribution of alternative splicing to genetic risk for psychiatric disorders

Single-nucleus transcriptomic profiling of human orbitofrontal cortex reveals convergent effects of aging and psychiatric disease

Aging is a complex biological process and represents the largest risk factor for neurodegenerative disorders. The risk for neurodegenerative disorders is also increased in individuals with psychiatric disorders. Here, we characterized age-related transcriptomic changes in the brain by profiling ~800,000 nuclei from the orbitofrontal cortex from 87 individuals with and without psychiatric diagnoses and replicated findings in an independent cohort with 32 individuals. Aging affects all cell types, with LAMP5+LHX6+ interneurons, a cell-type abundant in primates, by far the most affected. Disrupted synaptic transmission emerged as a convergently affected pathway in aged tissue. Age-related transcriptomic changes overlapped with changes observed in Alzheimer's disease across multiple cell types. We find evidence for accelerated transcriptomic aging in individuals with psychiatric disorders and demonstrate a converging signature of aging and psychopathology across multiple cell types. Our findings shed light on cell-type-specific effects and biological pathways underlying age-related changes and their convergence with effects driven by psychiatric diagnosis.

Integrative analysis of transcriptome-wide association study and mRNA expression profile identified risk genes for bipolar disorder

OBJECTIVE

Bipolar disorder (BD) is a debilitating neuropsychiatric disorder, which is associated with genetic variation through "vast but mixed" Genome-Wide Association Studies (GWAS). Transcriptome-Wide Association Study (TWAS) is more effective in explaining genetic factors that influence complex diseases and can help identifying risk genes more reliably. So, this study aims to identify potential BD risk genes in pedigrees with TWAS.

METHODS

We conducted a TWAS analysis with expression quantitative trait loci (eQTL) analysis on extended BD pedigrees, and the BD genome-wide association study (GWAS) summary data acquired from the Psychiatric Genomics Consortium (PGC). Furthermore, the BD-associated genes identified by TWAS were validated by mRNA expression profiles from the Gene Expression Omnibus (GEO) Datasets (GSE23848 and GSE46416). Functional enrichment and annotation analysis were implemented by RStudio (version 4.2.0).

RESULTS

TWAS identified 362 genes with P value < 0.05, and 18 genes remain significant after Bonferroni correction, such as SEMA3G (PTWAS=1.07 × 10-11), ALOX5AP (PTWAS=3.12 × 10-8), and PLEC (PTWAS=1.27 × 10-7). Further 6 overlapped genes were detected in integrative analysis, such as UQCRB (PTWAS=0.0020, PmRNA=0.0000), TMPRSS9 (PTWAS=0.0405, PmRNA=0.0032), and SNX10 (PTWAS=0.0104, PmRNA=0.0015). Using genes identified by TWAS, Gene Ontology (GO) enrichment analysis identified 40 significant GO terms, such as mitochondrial ATP synthesis coupled electron transport, mitochondrial respiratory, aerobic electron transport chain, oxidative phosphorylation, mitochondrial membrane proteins, and ubiquinone activity. The Kyoto Encyclopedia of Genes and Genomes (KEGG) Pathway enrichment analysis identified significant 15 pathways for BD, such as Oxidative phosphorylation, endocannabinoids signaling, neurodegeneration, and reactive oxide species.

CONCLUSIONS

We found a set of BD-associated genes and pathways, validating the important role of neurodevelopmental abnormalities, inflammatory responses, and mitochondrial dysfunction in the pathology of BD, offering novel information for comprehending the genetic basis of BD.

Exposure to elevated glucocorticoid during development primes altered transcriptional responses to acute stress in adulthood

Early life stress (ELS) is a major risk factor for developing psychiatric disorders, with glucocorticoids (GCs) implicated in mediating its effects in shaping adult phenotypes. In this process, exposure to high levels of developmental GC (hdGC) is thought to induce molecular changes that prime differential adult responses. However, identities of molecules targeted by hdGC exposure are not completely known. Here, we describe lifelong molecular consequences of hdGC exposure using a newly developed zebrafish double-hit stress model, which shows altered behaviors and stress hypersensitivity in adulthood. We identify a set of primed genes displaying altered expression only upon acute stress in hdGC-exposed adult fish brains. Interestingly, this gene set is enriched in risk factors for psychiatric disorders in humans. Lastly, we identify altered epigenetic regulatory elements following hdGC exposure. Thus, our study provides comprehensive datasets delineating potential molecular targets mediating the impact of hdGC exposure on adult responses.

Genome-wide landscape of RNA-binding protein (RBP) networks as potential molecular regulators of psychiatric co-morbidities: a computational analysis

Background

Psychiatric disorders are a major burden on global health. These illnesses manifest as co-morbid conditions, further complicating the treatment. There is a limited understanding of the molecular and regulatory basis of psychiatric co-morbidities. The existing research in this regard has largely focused on epigenetic modulators, non-coding RNAs, and transcription factors. RNA-binding proteins (RBPs) functioning as multi-protein complexes are now known to be predominant controllers of multiple gene regulatory processes. However, their involvement in gene expression dysregulation in psychiatric co-morbidities is yet to be understood.

Results

Ten RBPs (QKI, ELAVL2, EIF2S1, SRSF3, IGF2BP2, EIF4B, SNRNP70, FMR1, DAZAP1, and MBNL1) were identified to be associated with psychiatric disorders such as schizophrenia, major depression, and bipolar disorders. Analysis of transcriptomic changes in response to individual depletion of these RBPs showed the potential influence of a large number of RBPs driving differential gene expression, suggesting functional cross-talk giving rise to multi-protein networks. Subsequent transcriptome analysis of post-mortem human brain samples from diseased and control individuals also suggested the involvement of ~ 100 RBPs influencing gene expression changes. These RBPs were found to regulate various processes including transcript splicing, mRNA transport, localization, stability, and translation. They were also found to form an extensive interactive network. Further, hnRNP, SRSF, and PCBP family RBPs, Matrin3, U2AF2, KHDRBS1, PTBP1, and also PABPN1 were found to be the hub proteins of the RBP network.

Conclusions

Extensive RBP networks involving a few hub proteins could result in transcriptome-wide dysregulation of post-transcriptional modifications, potentially driving multiple psychiatric disorders. Understanding the functional involvement of RBP networks in psychiatric disorders would provide insights into the molecular basis of psychiatric co-morbidities.

From bugs to bedside: functional annotation of human genetic variation for neurological disorders using invertebrate models

The exponential accumulation of DNA sequencing data has opened new avenues for discovering the causative roles of single-nucleotide polymorphisms (SNPs) in neurological diseases. The opportunities emerging from this are staggering, yet only as good as our abilities to glean insights from this surplus of information. Whereas computational biology continues to improve with respect to predictions and molecular modeling, the differences between in silico and in vivo analysis remain substantial. Invertebrate in vivo model systems represent technically advanced, experimentally mature, high-throughput, efficient and cost-effective resources for investigating a disease. With a decades-long track record of enabling investigators to discern function from DNA, fly (Drosophila) and worm (Caenorhabditis elegans) models have never been better poised to serve as living engines of discovery. Both of these animals have already proven useful in the classification of genetic variants as either pathogenic or benign across a range of neurodevelopmental and neurodegenerative disorders-including autism spectrum disorders, ciliopathies, amyotrophic lateral sclerosis, Alzheimer's and Parkinson's disease. Pathogenic SNPs typically display distinctive phenotypes in functional assays when compared with null alleles and frequently lead to protein products with gain-of-function or partial loss-of-function properties that contribute to neurological disease pathogenesis. The utility of invertebrates is logically limited by overt differences in anatomical and physiological characteristics, and also the evolutionary distance in genome structure. Nevertheless, functional annotation of disease-SNPs using invertebrate models can expedite the process of assigning cellular and organismal consequences to mutations, ascertain insights into mechanisms of action, and accelerate therapeutic target discovery and drug development for neurological conditions.

The Long Non-Coding RNA GOMAFU in Schizophrenia: Function, Disease Risk, and Beyond

Neuropsychiatric diseases are among the most common brain developmental disorders, represented by schizophrenia (SZ). The complex multifactorial etiology of SZ remains poorly understood, which reflects genetic vulnerabilities and environmental risks that affect numerous genes and biological pathways. Besides the dysregulation of protein-coding genes, recent discoveries demonstrate that abnormalities associated with non-coding RNAs, including microRNAs and long non-coding RNAs (lncRNAs), also contribute to the pathogenesis of SZ. lncRNAs are an actively evolving family of non-coding RNAs that harbor greater than 200 nucleotides but do not encode for proteins. In general, lncRNA genes are poorly conserved. The large number of lncRNAs specifically expressed in the human brain, together with the genetic alterations and dysregulation of lncRNA genes in the SZ brain, suggests a critical role in normal cognitive function and the pathogenesis of neuropsychiatric diseases. A particular lncRNA of interest is GOMAFU, also known as MIAT and RNCR2. Growing evidence suggests the function of GOMAFU in governing neuronal development and its potential roles as a risk factor and biomarker for SZ, which will be reviewed in this article. Moreover, we discuss the potential mechanisms through which GOMAFU regulates molecular pathways, including its subcellular localization and interaction with RNA-binding proteins, and how interruption to GOMAFU pathways may contribute to the pathogenesis of SZ.

De novo truncating NOVA2 variants affect alternative splicing and lead to heterogeneous neurodevelopmental phenotypes

Alternative splicing (AS) is crucial for cell-type specific gene transcription and plays a critical role in neuronal differentiation and synaptic plasticity. De novo frameshift variants in NOVA2, encoding a neuron-specific key splicing factor, have been recently associated with a new neurodevelopmental disorder (NDD) with hypotonia, neurological features, and brain abnormalities. We investigated eight unrelated individuals by exome sequencing (ES) and identified seven novel pathogenic NOVA2 variants, including two with a novel localization at the KH1 and KH3 domains. In addition to a severe NDD phenotype, novel clinical features included psychomotor regression, attention deficit-hyperactivity disorder (ADHD), dyspraxia, and urogenital and endocrinological manifestations. To test the effect of the variants on splicing regulation, we transfected HeLa cells with wildtype and mutant NOVA2 cDNA. The novel variants NM_002516.4: c.754_756delCTGinsTT p.(Leu252Phefs*144) and c.1329dup p.(Lys444Glnfs*82) all negatively affected AS events. The distal p.(Lys444Glnfs*82) variant, causing a partial removal of the KH3 domain, had a milder functional effect leading to an intermediate phenotype. Our findings expand the molecular and phenotypic spectrum of NOVA2-related NDD, supporting the pathogenic role of AS disruption by truncating variants and suggesting that this is a heterogeneous condition with variable clinical course. This article is protected by copyright. All rights reserved.

An alternative splicing hypothesis for neuropathology of schizophrenia: evidence from studies on historical candidate genes and multi-omics data

Alternative splicing of schizophrenia risk genes, such as DRD2, GRM3, and DISC1, has been extensively described. Nevertheless, the alternative splicing characteristics of the growing number of schizophrenia risk genes identified through genetic analyses remain relatively opaque. Recently, transcriptomic analyses in human brains based on short-read RNA-sequencing have discovered many "local splicing" events (e.g., exon skipping junctions) associated with genetic risk of schizophrenia, and further molecular characterizations have identified novel spliced isoforms, such as AS3MT^d2d3 and ZNF804A^E3E4. In addition, long-read sequencing analyses of schizophrenia risk genes (e.g., CACNA1C and NRXN1) have revealed multiple previously unannotated brain-abundant isoforms with therapeutic potentials, and functional analyses of KCNH2-3.1 and Ube3a1 have provided examples for investigating such spliced isoforms in vitro and in vivo. These findings suggest that alternative splicing may be an essential molecular mechanism underlying genetic risk of schizophrenia, however, the incomplete annotations of human brain transcriptomes might have limited our understanding of schizophrenia pathogenesis, and further efforts to elucidate these transcriptional characteristics are urgently needed to gain insights into the illness-correlated brain physiology and pathology as well as to translate genetic discoveries into novel therapeutic targets.

The impact of psychosocial defeat stress on the bed nucleus of the stria terminalis transcriptome in adult male mice

The bed nucleus of the stria terminalis (BNST) is a focal point for the convergence of inputs from canonical stress-sensitive structures to fine-tune the response to stress. However, its role in mediating phenotypes of stress resilience or susceptibility is yet to be fully defined. In this study, we carried out unbiased RNA-sequencing to analyse the BNST transcriptomes of adult male mice, which were classified as resilient or susceptible following a 10-day chronic psychosocial defeat stress paradigm. Pairwise comparisons revealed 20 differentially expressed genes in resilience (6) and susceptible (14) mice compared with controls. An in silico validation of our data against an earlier study revealed significant concordance in gene expression profiles associated with resilience to chronic stress. Enrichment analysis revealed that resilience is linked to functions including retinoic acid hydrolase activity, phospholipase inhibitor and tumour necrosis factor (TNF)-receptor activities, whereas susceptibility is linked to alterations in amino acid transporter activity. We also identified differential usage of 134 exons across 103 genes associated with resilience and susceptibility; enrichment analysis for genes with differential exon usage in resilient mice was linked to functions including adrenergic receptor binding mice and oxysterol binding in susceptible mice. Our findings highlight the important and underappreciated role of the BNST in stress resilience and susceptibility and reveal research avenues for follow-up investigations.

Do Fragile X Syndrome and Other Intellectual Disorders Converge at Aberrant Pre-mRNA Splicing?

Fragile X Syndrome is a neuro-developmental disorder caused by the silencing of the FMR1 gene, resulting in the loss of its protein product, FMRP. FMRP binds mRNA and represses general translation in the brain. Transcriptome analysis of the Fmr1-deficient mouse hippocampus reveals widespread dysregulation of alternative splicing of pre-mRNAs. Many of these aberrant splicing changes coincide with those found in post-mortem brain tissue from individuals with autism spectrum disorders (ASDs) as well as in mouse models of intellectual disability such as PTEN hamartoma syndrome (PHTS) and Rett Syndrome (RTT). These splicing changes could result from chromatin modifications (e.g., in FXS, RTT) and/or splicing factor alterations (e.g., PTEN, autism). Based on the identities of the RNAs that are mis-spliced in these disorders, it may be that they are at least partly responsible for some shared pathophysiological conditions. The convergence of splicing aberrations among these autism spectrum disorders might be crucial to understanding their underlying cognitive impairments.

Sex-dependent complex association of TPH2 with multiple dimensions of ADHD

BACKGROUND

Tryptophan hydroxylase 2 (TPH2) is a key enzyme in the biosynthesis of serotonin in the brain. This study aims to investigate the role of a functional variant in TPH2 (rs17110747) in the pathophysiology of ADHD. This variant has been implicated in mood disorders in recent meta-analysis. This study uses a comprehensive approach that combines association testing and pharmaco-dynamic evaluation of behaviour, in a large sample of children with ADHD (n = 570).

METHODS

The association between various ADHD relevant traits and rs17110747 was analyzed using family-based association tests (FBAT). Children were assessed by parents, teachers and research staff under three experimental conditions (EC): baseline, placebo, and methylphenidate using a double-blind placebo-controlled crossover trial.

OUTCOMES

FBAT analysis conducted in a sample stratified based on sex of the proband, showed that there was a highly significant overtransmission of the G allele from parents to affected girls. In addition, significant association with several behavioral and cognitive dimensions of ADHD was observed only when the proband was female. Further, girls with the G/G genotype (rs17110747) had greater response to placebo when evaluated by parents.

INTERPRETATION

These results suggest that there may be a complex association of TPH2 in the etiology of ADHD, with a sex-specific effect.

Single molecule in situ hybridization reveals distinct localizations of schizophrenia risk-related transcripts SNX19 and AS3MT in human brain

Genome-wide association studies have identified single nucleotide polymorphisms (SNPs) associated with schizophrenia risk. Integration of RNA-sequencing data from postmortem human brains with these risk SNPs identified transcripts associated with increased schizophrenia susceptibility, including a class of exon 9-spliced isoforms of Sorting nexin-19 (SNX19d9) and an isoform of Arsenic methyltransferase (AS3MT) splicing out exons 2 and 3 (AS3MTd2d3). However, the biological function of these transcript variants is unclear. Defining the cell types where these risk transcripts are dominantly expressed is an important step to understand function, in prioritizing specific cell types and/or neural pathways in subsequent studies. To identify the cell type-specific localization of SNX19 and AS3MT in the human dorsolateral prefrontal cortex (DLPFC), we used single-molecule in situ hybridization techniques combined with automated quantification and machine learning approaches to analyze 10 postmortem brains of neurotypical individuals. These analyses revealed that both pan-SNX19 and pan-AS3MT were more highly expressed in neurons than non-neurons in layers II/III and VI of DLPFC. Furthermore, pan-SNX19 was preferentially expressed in glutamatergic neurons, while pan-AS3MT was preferentially expressed in GABAergic neurons. Finally, we utilized duplex BaseScope technology, to delineate the localization of SNX19d9 and AS3MTd2d3 splice variants, revealing consistent trends in spatial gene expression among pan-transcripts and schizophrenia risk-related transcript variants. These findings demonstrate that schizophrenia risk transcripts have distinct localization patterns in the healthy human brains, and suggest that SNX19 transcripts might disrupt the normal function of glutamatergic neurons, while AS3MT may lead to disturbances in the GABAergic system in the pathophysiology of schizophrenia.

The RNA-Binding Protein HuD Regulates Alternative Splicing and Alternative Polyadenylation in the Mouse Neocortex

The neuronal Hu/ELAV-like proteins HuB, HuC and HuD are a class of RNA-binding proteins that are crucial for proper development and maintenance of the nervous system. These proteins bind to AU-rich elements (AREs) in the untranslated regions (3'-UTRs) of target mRNAs regulating mRNA stability, transport and translation. In addition to these cytoplasmic functions, Hu proteins have been implicated in alternative splicing and alternative polyadenylation in the nucleus. The purpose of this study was to identify transcriptome-wide effects of HuD deletion on both of these nuclear events using RNA sequencing data obtained from the neocortex of Elavl4-/- (HuD KO) mice. HuD KO affected alternative splicing of 310 genes, including 17 validated HuD targets such as Cbx3, Cspp1, Snap25 and Gria2. In addition, deletion of HuD affected polyadenylation of 53 genes, with the majority of significantly altered mRNAs shifting towards usage of proximal polyadenylation signals (PAS), resulting in shorter 3'-UTRs. None of these genes overlapped with those showing alternative splicing events. Overall, HuD KO had a greater effect on alternative splicing than polyadenylation, with many of the affected genes implicated in several neuronal functions and neuropsychiatric disorders.

Meta-Analyses of Splicing and Expression Quantitative Trait Loci Identified Susceptibility Genes of Glioma

Background

The functions of most glioma risk alleles are unknown. Very few studies had evaluated expression quantitative trait loci (eQTL), and insights of susceptibility genes were limited due to scarcity of available brain tissues. Moreover, no prior study had examined the effect of glioma risk alleles on alternative RNA splicing.

Objective

This study explored splicing quantitative trait loci (sQTL) as molecular QTL and improved the power of QTL mapping through meta-analyses of both cis eQTL and sQTL.

Methods

We first evaluated eQTLs and sQTLs of the CommonMind Consortium (CMC) and Genotype-Tissue Expression Project (GTEx) using genotyping, or whole-genome sequencing and RNA-seq data. Alternative splicing events were characterized using an annotation-free method that detected intron excision events. Then, we conducted meta-analyses by pooling the eQTL and sQTL results of CMC and GTEx using the inverse variance-weighted model. Afterward, we integrated QTL meta-analysis results (Q < 0.05) with the Glioma International Case Control Study (GICC) GWAS meta-analysis (case:12,496, control:18,190), using a summary statistics-based mendelian randomization (SMR) method.

Results

Between CMC and GTEx, we combined the QTL data of 354 unique individuals of European ancestry. SMR analyses revealed 15 eQTLs in 11 loci and 32 sQTLs in 9 loci relevant to glioma risk. Two loci only harbored sQTLs (1q44 and 16p13.3). In seven loci, both eQTL and sQTL coexisted (2q33.3, 7p11.2, 11q23.3 15q24.2, 16p12.1, 20q13.33, and 22q13.1), but the target genes were different for five of these seven loci. Three eQTL loci (9p21.3, 20q13.33, and 22q13.1) and 4 sQTL loci (11q23.3, 16p13.3, 16q12.1, and 20q13.33) harbored multiple target genes. Eight target genes of sQTLs (C2orf80, SEC61G, TMEM25, PHLDB1, RP11-161M6.2, HEATR3, RTEL1-TNFRSF6B, and LIME1) had multiple alternatively spliced transcripts.

Conclusion

Our study revealed that the regulation of transcriptome by glioma risk alleles is complex, with the potential for eQTL and sQTL jointly affecting gliomagenesis in risk loci. QTLs of many loci involved multiple target genes, some of which were specific to alternative splicing. Therefore, quantitative trait loci that evaluate only total gene expression will miss many important target genes.

Cannabis Use Induces Distinctive Proteomic Alterations in Olfactory Neuroepithelial Cells of Schizophrenia Patients

A close epidemiological link has been reported between cannabis use and schizophrenia (SCZ). However, biochemical markers in living humans related to the impact of cannabis in this disease are still missing. Olfactory neuroepithelium (ON) cells express neural features and offer a unique advantage to study biomarkers of psychiatric diseases. The aim of our study was to find exclusively deregulated proteins in ON cells of SCZ patients with and without a history of cannabis use. Thus, we compared the proteomic profiles of SCZ non-cannabis users (SCZ/nc) and SCZ cannabis users (SCZ/c) with control subjects non-cannabis users (C/nc) and control cannabis users (C/c). The results revealed that the main cascades affected in SCZ/nc were cell cycle, DNA replication, signal transduction and protein localization. Conversely, cannabis use in SCZ patients induced specific alterations in metabolism of RNA and metabolism of proteins. The levels of targeted proteins in each population were then correlated with cognitive performance and clinical scores. In SCZ/c, the expression levels of 2 proteins involved in the metabolism of RNA (MTREX and ZNF326) correlated with several cognitive markers and clinical signs. Moreover, use duration of cannabis negatively correlated with ZNF326 expression. These findings indicate that RNA-related proteins might be relevant to understand the influence of cannabis use on SCZ.

Aberrant alternative splicing of HTR2A exon II in peripheral blood lymphocytes of drug-naïve schizophrenic patients

The aim of the study was to characterize the pattern of transcript isoforms of HTR2A exon II in lymphocytes of the blood as peripheral biomarkers of schizophrenia development and the effectiveness of antipsychotic therapy. We primarily observed an increase in mRNA levels and elevation of alternative variants in a sample of drug-naïve schizophrenic patients compared to the control group. There was no association of the expression level of HTR2A transcript isoforms with the effectiveness of the antipsychotic therapy. Antipsychotic-induced akathisia was associated with a significant reduction in the mRNA levels of the studied isoforms. In summary, our results suggest that levels of HTR2A exon II transcript isoforms are upregulated in patients with schizophrenia, but at the same time, elevated expression level of the studied HTR2A transcripts is associated with fewer side effects of the therapy.

The Secret Garden of Neuronal circRNAs

High-throughput transcriptomic profiling approaches have revealed that circular RNAs (circRNAs) are important transcriptional gene products, identified across a broad range of organisms throughout the eukaryotic tree of life. In the nervous system, they are particularly abundant, developmentally regulated, region-specific, and enriched in genes for neuronal proteins and synaptic factors. These features suggested that circRNAs are key components of an important layer of neuronal gene expression regulation, with known and anticipated functions. Here, we review major recognized aspects of circRNA biogenesis, metabolism and biological activities, examining potential new functions in the context of the nervous system.

Ubiquitin-proteasome system, lipid metabolism and DNA damage repair are triggered by antipsychotic medication in human oligodendrocytes: implications in schizophrenia

Schizophrenia is a chronic, severe and disabling psychiatric disorder, whose treatment is based on psychosocial interventions and the use of antipsychotic drugs. While the effects of these drugs are well elucidated in neuronal cells, they are still not so clear in oligodendrocytes, which play a vital role in schizophrenia. Thus, we aimed to characterize biochemical profiles by proteomic analyses of human oligodendrocytes (MO3.13) which were matured using a protocol we developed and treated with either haloperidol (a typical antipsychotic), clozapine (an atypical antipsychotic) or a clozapine + D-serine co-treatment, which has emerged lately as an alternative type of treatment. This was accomplished by employing shotgun proteomics, using nanoESI-LC-MS/MS label-free quantitation. Proteomic analysis revealed biochemical pathways commonly affected by all tested antipsychotics were mainly associated to ubiquitination, proteasome degradation, lipid metabolism and DNA damage repair. Clozapine and haloperidol treatments also affected proteins involved with the actin cytoskeleton and with EIF2 signaling. In turn, metabolic processes, especially the metabolism of nitrogenous compounds, were a predominant target of modulation of clozapine + D-serine treatment. In this context, we seek to contribute to the understanding of the biochemical and molecular mechanisms involved in the action of antipsychotics on oligodendrocytes, along with their possible implications in schizophrenia.

Predicting the effect of variants on splicing using Convolutional Neural Networks

Mutations that cause an error in the splicing of a messenger RNA (mRNA) can lead to diseases in humans. Various computational models have been developed to recognize the sequence pattern of the splice sites. In recent studies, Convolutional Neural Network (CNN) architectures were shown to outperform other existing models in predicting the splice sites. However, an insufficient effort has been put into extending the CNN model to predict the effect of the genomic variants on the splicing of mRNAs. This study proposes a framework to elaborate on the utility of CNNs to assess the effect of splice variants on the identification of potential disease-causing variants that disrupt the RNA splicing process. Five models, including three CNN-based and two non-CNN machine learning based, were trained and compared using two existing splice site datasets, Genome Wide Human splice sites (GWH) and a dataset provided at the Deep Learning and Artificial Intelligence winter school 2018 (DLAI). The donor sites were also used to test on the HSplice tool to evaluate the predictive models. To improve the effectiveness of predictive models, two datasets were combined. The CNN model with four convolutional layers showed the best splice site prediction performance with an AUPRC of 93.4% and 88.8% for donor and acceptor sites, respectively. The effects of variants on splicing were estimated by applying the best model on variant data from the ClinVar database. Based on the estimation, the framework could effectively differentiate pathogenic variants from the benign variants (p = 5.9 × 10⁻⁷). These promising results support that the proposed framework could be applied in future genetic studies to identify disease causing loci involving the splicing mechanism. The datasets and Python scripts used in this study are available on the GitHub repository at https://github.com/smiile8888/rna-splice-sites-recognition.

Long non-coding RNA Neat1 regulates adaptive behavioural response to stress in mice

NEAT1 is a highly and ubiquitously expressed long non-coding RNA (lncRNA) which serves as an important regulator of cellular stress response. However, the physiological role of NEAT1 in the central nervous system (CNS) is still poorly understood. In the current study, we addressed this by characterising the CNS function of the Neat1 knockout mouse model (Neat1-/- mice), using a combination of behavioural phenotyping, electrophysiology and expression analysis. RNAscope® in situ hybridisation revealed that in wild-type mice, Neat1 is expressed across the CNS regions, with high expression in glial cells and low expression in neurons. Loss of Neat1 in mice results in an inadequate reaction to physiological stress manifested as hyperlocomotion and panic escape response. In addition, Neat1-/- mice display deficits in social interaction and rhythmic patterns of activity but retain normal motor function and memory. Neat1-/- mice do not present with neuronal loss, overt neuroinflammation or gross synaptic dysfunction in the brain. However, cultured Neat1-/- neurons are characterised by hyperexcitability and dysregulated calcium homoeostasis, and stress-induced neuronal activity is also augmented in Neat1-/- mice in vivo. Gene expression analysis showed that Neat1 may act as a weak positive regulator of multiple genes in the brain. Furthermore, loss of Neat1 affects alternative splicing of genes important for the CNS function and implicated in neurological diseases. Overall, our data suggest that Neat1 is involved in stress signalling in the brain and fine-tunes the CNS functions to enable adaptive behaviour in response to physiological stress.

Transcription of PIK3CD in human brain and schizophrenia: regulation by proinflammatory cytokines

PIK3CD encodes the phosphoinositide 3-kinase (PI3K) catalytic subunit, p110δ, a lipid kinase linked to neurodevelopmental disorders, including schizophrenia (SZ). PIK3CD is regulated at the transcript level through alternate use of 5' untranslated exons (UTRs), promoters, and proinflammatory cytokines. Increases in global PIK3CD expression and downregulation by neuroleptics are observed in SZ, and preclinical efficacy of a p110δ-selective inhibitor is seen in rodent models of risk. Here, we cloned PIK3CD alternative transcripts in human brain and evaluated temporal- and tissue-specific expression. We quantified PIK3CD transcripts in B-lymphoblastoid cells from patients with SZ and examined 5' UTR transcriptional regulation by tumor necrosis factor α (TNFα) and interleukin-1β (IL1β) in patient-derived fibroblasts. We report that PIK3CD transcripts are differentially expressed in human brain in a developmental-specific manner. Transcripts encoding 5' UTRs -2A and alternative exon -1 (Alt1), P37 and AS1 and AS2 were increased in SZ. Alt1, P37, and AS2 were also preferentially expressed in fetal brain, and all transcripts were regulated by TNFα and IL1β. Our findings provide novel insight into the complexity of PIK3CD regulation in human brain, implicate PIK3CD in human neurodevelopment, and identify isoform-specific disruption in SZ.

Pre-mRNA splicing defects and RNA binding protein involvement in Niemann Pick type C disease

Niemann-Pick type C (NPC) is an autosomal recessive lysosomal storage disorder due to mutations in NPC1 (95 % cases) or NPC2 genes, encoding NPC1 and NPC2 proteins, respectively. Both NPC1 and NPC2 proteins are involved in transport of intracellular cholesterol and their alteration leads to the accumulation of unesterified cholesterol and other lipids within the lysosomes. The disease is characterized by visceral, neurological and psychiatric symptoms. However, the pathogenic mechanisms that lead to the fatal neurodegeneration are still unclear. To date, several mutations leading to the generation of aberrant splicing variants or mRNA degradation in NPC1 and NPC2 genes have been reported. In addition, different lines of experimental evidence have highlighted the possible role of RNA-binding proteins and RNA-metabolism, in the onset and progression of many neurodegenerative disorders, that could explain NPC neurological features and in general, the disease pathogenesis. In this review, we will provide an overview of the impact of mRNA processing and metabolism on NPC disease pathology.

Autism-associated missense genetic variants impact locomotion and neurodevelopment in Caenorhabditis elegans

Autism spectrum disorder (ASD) involves thousands of alleles in over 850 genes, but the current functional inference tools are not sufficient to predict phenotypic changes. As a result, the causal relationship of most of these genetic variants in the pathogenesis of ASD has not yet been demonstrated and an experimental method prioritizing missense alleles for further intensive analysis is crucial. For this purpose, we have designed a pipeline that uses Caenorhabditis elegans as a genetic model to screen for phenotype-changing missense alleles inferred from human ASD studies. We identified highly conserved human ASD-associated missense variants in their C. elegans orthologs, used a CRISPR/Cas9-mediated homology-directed knock-in strategy to generate missense mutants and analyzed their impact on behaviors and development via several broad-spectrum assays. All tested missense alleles were predicted to perturb protein function, but we found only 70% of them showed detectable phenotypic changes in morphology, locomotion or fecundity. Our findings indicate that certain missense variants in the C. elegans orthologs of human CACNA1D, CHD7, CHD8, CUL3, DLG4, GLRA2, NAA15, PTEN, SYNGAP1 and TPH2 impact neurodevelopment and movement functions, elevating these genes as candidates for future study into ASD. Our approach will help prioritize functionally important missense variants for detailed studies in vertebrate models and human cells.

Accelerated evolution of oligodendrocytes in the human brain

Recent discussions of human brain evolution have largely focused on increased neuron numbers and changes in their connectivity and expression. However, it is increasingly appreciated that oligodendrocytes play important roles in cognitive function and disease. Whether both cell types follow similar or distinctive evolutionary trajectories is not known. We examined the transcriptomes of neurons and oligodendrocytes in the frontal cortex of humans, chimpanzees, and rhesus macaques. We identified human-specific trajectories of gene expression in neurons and oligodendrocytes and show that both cell types exhibit human-specific up-regulation. Moreover, oligodendrocytes have undergone more pronounced accelerated gene expression evolution in the human lineage compared to neurons. We highlighted human-specific coexpression networks with specific functions. Our data suggest that oligodendrocyte human-specific networks are enriched for alternative splicing and transcriptional regulation. Oligodendrocyte networks are also enriched for variants associated with schizophrenia and other neuropsychiatric disorders. Such enrichments were not found in neuronal networks. These results offer a glimpse into the molecular mechanisms of oligodendrocytes during evolution and how such mechanisms are associated with neuropsychiatric disorders.

DNA Variant in the RPGRIP1L Gene Influences Alternative Splicing

The retinitis pigmentosa GTPase regulator interacting protein 1-like (RPGRIP1L) gene encodes a ciliary protein that is critical for processes related to brain development, including development of left-right asymmetry, sonic hedgehog signaling, and neural tube formation. RPGRIP1L is a risk factor for retinal degeneration, and rare, deleterious variants in the RPGRIP1L gene cause Joubert syndrome and Meckel syndrome, both autosomal recessive disorders. These syndromes are characterized by dysfunctional primary cilia that result in abnormal development - and even lethality in the case of Meckel syndrome. Genetic studies have also implicated RPGRIP1L in psychiatric disorders by suggestive findings from genome-wide association studies and findings from rare-variant exome analyses for bipolar disorder and de novo mutations in autism. In this study we identify a common variant in RPGRIP1L, rs7203525, that influences alternative splicing, increasing the inclusion of exon 20 of RPGRIP1L. We detected this alternative splicing association in human postmortem brain tissue samples and, using a minigene assay combined with in vitro mutagenesis, confirmed that the alternative splicing is attributable to the alleles of this variant. The predominate RPGRIP1L isoform expressed in adult brains does not contain exon 20; thus, a shift to include this exon may impact brain function.

Genetic and Epigenetic Alterations Underlie Oligodendroglia Susceptibility and White Matter Etiology in Psychiatric Disorders

Numerous genetic risk loci are found to associate with major neuropsychiatric disorders represented by schizophrenia. The pathogenic roles of genetic risk loci in psychiatric diseases are further complicated by the association with cell lineage- and/or developmental stage-specific epigenetic alterations. Besides aberrant assembly and malfunction of neuronal circuitry, an increasing volume of discoveries clearly demonstrate impairment of oligodendroglia and disruption of white matter integrity in psychiatric diseases. Nonetheless, whether and how genetic risk factors and epigenetic dysregulations for neuronal susceptibility may affect oligodendroglia is largely unknown. In this mini-review, we will discuss emerging evidence regarding the functional interplay between genetic risk loci and epigenetic factors, which may underlie compromised oligodendroglia and myelin development in neuropsychiatric disorders. Transcriptional and epigenetic factors are the major aspects affected in oligodendroglia. Moreover, multiple disease susceptibility genes are connected by epigenetically modulated transcriptional and post-transcriptional mechanisms. Oligodendroglia specific complex molecular orchestra may explain how distinct risk factors lead to the common clinical expression of white matter pathology of neuropsychiatric disorders.