A novel protein domain in an ancestral splicing factor drove the evolution of neural microexons

SRRM2 splicing factor modulates cell fate in early development

Embryo development is an orchestrated process that relies on tight regulation of gene expression to guide cell differentiation and fate decisions. The Srrm2 splicing factor has recently been implicated in developmental disorders and diseases, but its role in early mammalian development remains unexplored. Here, we show that Srrm2 dosage is critical for maintaining embryonic stem cell pluripotency and cell identity. Srrm2 heterozygosity promotes loss of stemness, characterised by the coexistence of cells expressing naive and formative pluripotency markers, together with extensive changes in gene expression, including genes regulated by serum-response transcription factor (SRF) and differentiation-related genes. Depletion of Srrm2 by RNA interference in embryonic stem cells shows that the earliest effects of Srrm2 heterozygosity are specific alternative splicing events on a small number of genes, followed by expression changes in metabolism and differentiation-related genes. Our findings unveil molecular and cellular roles of Srrm2 in stemness and lineage commitment, shedding light on the roles of splicing regulators in early embryogenesis, developmental diseases and tumorigenesis.

Unique transcriptomes of sensory and non-sensory neurons: insights from Splicing Regulatory States

Alternative Splicing (AS) programs serve as instructive signals of cell type specificity, particularly within the brain, which comprises dozens of molecularly and functionally distinct cell types. Among them, retinal photoreceptors stand out due to their unique transcriptome, making them a particularly well-suited system for studying how AS shapes cell type-specific molecular functions. Here, we use the Splicing Regulatory State (SRS) as a novel framework to discuss the splicing factors governing the unique AS pattern of photoreceptors, and how this pattern may aid in the specification of their highly specialized sensory cilia. In addition, we discuss how other sensory cells with ciliated structures, for which data is much scarcer, also rely on specific SRSs to implement a proteome specialized in the detection of sensory stimuli. By reviewing the general rules of cell type- and tissue-specific AS programs, firstly in the brain and subsequently in specialized sensory neurons, we propose a novel paradigm on how SRSs are established and how they can diversify. Finally, we illustrate how SRSs shape the outcome of mutations in splicing factors to produce cell type-specific phenotypes that can lead to various human diseases.

SRRM2 phase separation drives assembly of nuclear speckle subcompartments

Nuclear speckles (NSs) are nuclear biomolecular condensates that are postulated to form by macromolecular phase separation, although the detailed underlying forces driving NS formation remain elusive. SRRM2 and SON are 2 non-redundant scaffold proteins for NSs. How each individual protein governs assembly of the NS protein network and the functional relationship between SRRM2 and SON are largely unknown. Here, we uncover immiscible multiphases of SRRM2 and SON within NSs. SRRM2 and SON are functionally independent, specifically regulating alternative splicing of subsets of mRNA targets, respectively. We further show that SRRM2 forms multicomponent liquid phases in cells to drive NS subcompartmentalization, which is reliant on homotypic interaction and heterotypic non-selective protein-RNA complex coacervation-driven phase separation. SRRM2 serine/arginine-rich (RS) domains form higher-order oligomers and can be replaced by oligomerizable synthetic modules. The serine residues within the RS domains, however, play an irreplaceable role in fine-tuning the liquidity of NSs.

Memory consolidation in honey bees is enhanced by down-regulation of Down syndrome cell adhesion molecule and changes its alternative splicing

Down syndrome cell adhesion molecule (Dscam) gene encodes a cell adhesion molecule required for neuronal wiring. A remarkable feature of arthropod Dscam is massive alternative splicing generating thousands of different isoforms from three variable clusters of alternative exons. Dscam expression and diversity arising from alternative splicing have been studied during development, but whether they exert functions in adult brains has not been determined. Here, using honey bees, we find that Dscam expression is critically linked to memory retention as reducing expression by RNAi enhances memory after reward learning in adult worker honey bees. Moreover, alternative splicing of Dscam is altered in all three variable clusters after learning. Since identical Dscam isoforms engage in homophilic interactions, these results suggest a mechanism to alter inclusion of variable exons during memory consolidation to modify neuronal connections for memory retention.

CCAR-1 works together with the U2AF large subunit UAF-1 to regulate alternative splicing

The Cell Division Cycle and Apoptosis Regulator (CCAR) protein family members have recently emerged as regulators of alternative splicing and transcription, as well as having other key physiological functions. For example, mammalian CCAR2/DBC1 forms a complex with the zinc factor protein ZNF326 to integrate alternative splicing with RNA polymerase II transcriptional elongation in AT-rich regions of the DNA. Additionally, Caenorhabditis elegans CCAR-1, a homolog to mammalian CCAR2, facilitates the alternative splicing of the perlecan unc-52 gene. However, much about the CCAR family's role in alternative splicing is unknown. Here, we have examined the role of CCAR-1 in genome-wide alternative splicing in Caenorhabditis elegans and have identified new alternative splicing targets of CCAR-1 using RNA sequencing. Also, we found that CCAR-1 interacts with the spliceosome factors UAF-1 and UAF-2 using mass spectrometry, and that knockdown of ccar-1 affects alternative splicing patterns, motility, and proteostasis of UAF-1 mutant worms. Collectively, we demonstrate the role of CCAR-1 in regulating global alternative splicing in C. elegans and in conjunction with UAF-1.

Regulation of alternative splicing and polyadenylation in neurons

Cell-type-specific gene expression is a fundamental feature of multicellular organisms and is achieved by combinations of regulatory strategies. Although cell-restricted transcription is perhaps the most widely studied mechanism, co-transcriptional and post-transcriptional processes are also central to the spatiotemporal control of gene functions. One general category of expression control involves the generation of multiple transcript isoforms from an individual gene, whose balance and cell specificity are frequently tightly regulated via diverse strategies. The nervous system makes particularly extensive use of cell-specific isoforms, specializing the neural function of genes that are expressed more broadly. Here, we review regulatory strategies and RNA-binding proteins that direct neural-specific isoform processing. These include various classes of alternative splicing and alternative polyadenylation events, both of which broadly diversify the neural transcriptome. Importantly, global alterations of splicing and alternative polyadenylation are characteristic of many neural pathologies, and recent genetic studies demonstrate how misregulation of individual neural isoforms can directly cause mutant phenotypes.

High frequency of alternative splicing variants of the oncogene Focal Adhesion Kinase in neuroendocrine tumors of the pancreas and breast

The characteristic genetic abnormality of neuroendocrine neoplasms (NENs), a heterogeneous group of tumors found in various organs, remains to be identified. Here, based on the analysis of the splicing variants of an oncogene Focal Adhesion Kinase (FAK) in The Cancer Genome Atlas datasets that contain 9193 patients of 33 cancer subtypes, we found that Box 6/Box 7-containing FAK variants (FAK6/7) were observed in 7 (87.5%) of 8 pancreatic neuroendocrine carcinomas and 20 (11.76%) of 170 pancreatic ductal adenocarcinomas (PDACs). We tested FAK variants in 157 tumor samples collected from Chinese patients with pancreatic tumors, and found that FAK6/7 was positive in 34 (75.6%) of 45 pancreatic NENs, 19 (47.5%) of 40 pancreatic solid pseudopapillary neoplasms, and 2 (2.9%) of 69 PDACs. We further tested FAK splicing variants in breast neuroendocrine carcinoma (BrNECs), and found that FAK6/7 was positive in 14 (93.3%) of 15 BrNECs but 0 in 23 non-NEC breast cancers. We explored the underlying mechanisms and found that a splicing factor serine/arginine repetitive matrix protein 4 (SRRM4) was overexpressed in FAK6/7-positive pancreatic tumors and breast tumors, which promoted the formation of FAK6/7 in cells. These results suggested that FAK6/7 could be a biomarker of NENs and represent a potential therapeutic target for these orphan diseases.

NvPrdm14d-expressing neural progenitor cells contribute to non-ectodermal neurogenesis in Nematostella vectensis

Neurogenesis has been studied extensively in the ectoderm, from which most animals generate the majority of their neurons. Neurogenesis from non-ectodermal tissue is, in contrast, poorly understood. Here we use the cnidarian Nematostella vectensis as a model to provide new insights into the molecular regulation of non-ectodermal neurogenesis. We show that the transcription factor NvPrdm14d is expressed in a subpopulation of NvSoxB(2)-expressing endodermal progenitor cells and their NvPOU4-expressing progeny. Using a new transgenic reporter line, we show that NvPrdm14d-expressing cells give rise to neurons in the body wall and in close vicinity of the longitudinal retractor muscles. RNA-sequencing of NvPrdm14d::GFP-expressing cells and gene knockdown experiments provide candidate genes for the development and function of these neurons. Together, the identification of a population of endoderm-specific neural progenitor cells and of previously undescribed putative motoneurons in Nematostella provide new insights into the regulation of non-ectodermal neurogenesis.

Insplico: effective computational tool for studying splicing order of adjacent introns genome-wide with short and long RNA-seq reads

Although splicing occurs largely co-transcriptionally, the order by which introns are removed does not necessarily follow the order in which they are transcribed. Whereas several genomic features are known to influence whether or not an intron is spliced before its downstream neighbor, multiple questions related to adjacent introns' splicing order (AISO) remain unanswered. Here, we present Insplico, the first standalone software for quantifying AISO that works with both short and long read sequencing technologies. We first demonstrate its applicability and effectiveness using simulated reads and by recapitulating previously reported AISO patterns, which unveiled overlooked biases associated with long read sequencing. We next show that AISO around individual exons is remarkably constant across cell and tissue types and even upon major spliceosomal disruption, and it is evolutionarily conserved between human and mouse brains. We also establish a set of universal features associated with AISO patterns across various animal and plant species. Finally, we used Insplico to investigate AISO in the context of tissue-specific exons, particularly focusing on SRRM4-dependent microexons. We found that the majority of such microexons have non-canonical AISO, in which the downstream intron is spliced first, and we suggest two potential modes of SRRM4 regulation of microexons related to their AISO and various splicing-related features. Insplico is available on gitlab.com/aghr/insplico.

Small but mighty: microexons in glucose homeostasis

Many molecular mechanisms underlying blood glucose homeostasis remain elusive. Juan-Mateu et al. find that pancreatic islet cells utilize a regulatory program, originally identified in neurons, that involves alternative splicing of microexons in genes important for insulin secretion or diabetes risk.

Pancreatic microexons regulate islet function and glucose homeostasis

Pancreatic islets control glucose homeostasis by the balanced secretion of insulin and other hormones, and their abnormal function causes diabetes or hypoglycaemia. Here we uncover a conserved programme of alternative microexons included in mRNAs of islet cells, particularly in genes involved in vesicle transport and exocytosis. Islet microexons (IsletMICs) are regulated by the RNA binding protein SRRM3 and represent a subset of the larger neural programme that are particularly sensitive to SRRM3 levels. Both SRRM3 and IsletMICs are induced by elevated glucose levels, and depletion of SRRM3 in human and rat beta cell lines and mouse islets, or repression of particular IsletMICs using antisense oligonucleotides, leads to inappropriate insulin secretion. Consistently, mice harbouring mutations in Srrm3 display defects in islet cell identity and function, leading to hyperinsulinaemic hypoglycaemia. Importantly, human genetic variants that influence SRRM3 expression and IsletMIC inclusion in islets are associated with fasting glucose variation and type 2 diabetes risk. Taken together, our data identify a conserved microexon programme that regulates glucose homeostasis.

Human EWS-FLI protein recapitulates in Drosophila the neomorphic functions that induce Ewing sarcoma tumorigenesis

Ewing sarcoma (EwS) is a human malignant tumor typically driven by the Ewing sarcoma-Friend leukemia integration (EWS-FLI) fusion protein. A paucity of genetically modified animal models, partially owed to the high toxicity of EWS-FLI, hinders research on EwS. Here, we report a spontaneous mutant variant, EWS-FLI_1FS, that circumvents the toxicity issue in Drosophila. Through proteomic and genomic analyses, we show that human EWS-FLI_1FS interacts with the Drosophila homologues of EWS-FLI human protein partners, including core subunits of chromatin remodeling complexes, the transcription machinery, and the spliceosome; brings about a massive dysregulation of transcription that affects a significant fraction of known targets of EWS-FLI in human cells; and modulates splicing. We also show that EWS-FLI_1FS performs in Drosophila the two major neomorphic activities that it is known to have in human cells: activation of transcription from GGAA microsatellites and out competition of ETS transcription factors. We conclude that EWS-FLI_1FS reproduces in Drosophila the known oncogenic activities of EWS-FLI that drive EwS tumorigenesis in humans. These results open up an unprecedented opportunity to investigate EWS-FLI's oncogenic pathways in vivo in a genetically tractable organism.

SRRM2 organizes splicing condensates to regulate alternative splicing

SRRM2 is a nuclear-speckle marker containing multiple disordered domains, whose dysfunction is associated with several human diseases. Using mainly EGFP-SRRM2 knock-in HEK293T cells, we show that SRRM2 forms biomolecular condensates satisfying most hallmarks of liquid-liquid phase separation, including spherical shape, dynamic rearrangement, coalescence and concentration dependence supported by in vitro experiments. Live-cell imaging shows that SRRM2 organizes nuclear speckles along the cell cycle. As bona-fide splicing factor present in spliceosome structures, SRRM2 deficiency induces skipping of cassette exons with short introns and weak splice sites, tending to change large protein domains. In THP-1 myeloid-like cells, SRRM2 depletion compromises cell viability, upregulates differentiation markers, and sensitizes cells to anti-leukemia drugs. SRRM2 induces a FES splice isoform that attenuates innate inflammatory responses, and MUC1 isoforms that undergo shedding with oncogenic properties. We conclude that SRRM2 acts as a scaffold to organize nuclear speckles, regulating alternative splicing in innate immunity and cell homeostasis.

Specialization of the photoreceptor transcriptome by Srrm3-dependent microexons is required for outer segment maintenance and vision

Retinal photoreceptors have a distinct transcriptomic profile compared to other neuronal subtypes, likely reflecting their unique cellular morphology and function in the detection of light stimuli by way of the ciliary outer segment. We discovered a layer of this molecular specialization by revealing that the vertebrate retina expresses the largest number of tissue-enriched microexons of all tissue types. A subset of these microexons is included exclusively in photoreceptor transcripts, particularly in genes involved in cilia biogenesis and vesicle-mediated transport. This microexon program is regulated by Srrm3, a paralog of the neural microexon regulator Srrm4. Despite the fact that both proteins positively regulate retina microexons in vitro, only Srrm3 is highly expressed in mature photoreceptors. Its deletion in zebrafish results in widespread down-regulation of microexon inclusion from early developmental stages, followed by other transcriptomic alterations, severe photoreceptor defects, and blindness. These results shed light on the transcriptomic specialization and functionality of photoreceptors, uncovering unique cell type-specific roles for Srrm3 and microexons with implications for retinal diseases.

Alternative splicing as a source of phenotypic diversity

A major goal of evolutionary genetics is to understand the genetic processes that give rise to phenotypic diversity in multicellular organisms. Alternative splicing generates multiple transcripts from a single gene, enriching the diversity of proteins and phenotypic traits. It is well established that alternative splicing contributes to key innovations over long evolutionary timescales, such as brain development in bilaterians. However, recent developments in long-read sequencing and the generation of high-quality genome assemblies for diverse organisms has facilitated comparisons of splicing profiles between closely related species, providing insights into how alternative splicing evolves over shorter timescales. Although most splicing variants are probably non-functional, alternative splicing is nonetheless emerging as a dynamic, evolutionarily labile process that can facilitate adaptation and contribute to species divergence.

Hnrnpul1 controls transcription, splicing, and modulates skeletal and limb development in vivo

Mutations in RNA-binding proteins can lead to pleiotropic phenotypes including craniofacial, skeletal, limb, and neurological symptoms. Heterogeneous nuclear ribonucleoproteins (hnRNPs) are involved in nucleic acid binding, transcription, and splicing through direct binding to DNA and RNA, or through interaction with other proteins in the spliceosome. We show a developmental role for Hnrnpul1 in zebrafish, resulting in reduced body and fin growth and missing bones. Defects in craniofacial tendon growth and adult-onset caudal scoliosis are also seen. We demonstrate a role for Hnrnpul1 in alternative splicing and transcriptional regulation using RNA-sequencing, particularly of genes involved in translation, ubiquitination, and DNA damage. Given its cross-species conservation and role in splicing, it would not be surprising if it had a role in human development. Whole-exome sequencing detected a homozygous frameshift variant in HNRNPUL1 in 2 siblings with congenital limb malformations, which is a candidate gene for their limb malformations. Zebrafish Hnrnpul1 mutants suggest an important developmental role of hnRNPUL1 and provide motivation for exploring the potential conservation of ancient regulatory circuits involving hnRNPUL1 in human development.

A developmentally programmed splicing failure contributes to DNA damage response attenuation during mammalian zygotic genome activation

Transition from maternal to embryonic transcriptional control is crucial for embryogenesis. However, alternative splicing regulation during this process remains understudied. Using transcriptomic data from human, mouse, and cow preimplantation development, we show that the stage of zygotic genome activation (ZGA) exhibits the highest levels of exon skipping diversity reported for any cell or tissue type. Much of this exon skipping is temporary, leads to disruptive noncanonical isoforms, and occurs in genes enriched for DNA damage response in the three species. Two core spliceosomal components, Snrpb and Snrpd2, regulate these patterns. These genes have low maternal expression at ZGA and increase sharply thereafter. Microinjection of Snrpb/d2 messenger RNA into mouse zygotes reduces the levels of exon skipping at ZGA and leads to increased p53-mediated DNA damage response. We propose that mammalian embryos undergo an evolutionarily conserved, developmentally programmed splicing failure at ZGA that contributes to the attenuation of cellular responses to DNA damage.

Non-Canonical Splicing and Its Implications in Brain Physiology and Cancer

The advance of experimental and computational techniques has allowed us to highlight the existence of numerous different mechanisms of RNA maturation, which have been so far unknown. Besides canonical splicing, consisting of the removal of introns from pre-mRNA molecules, non-canonical splicing events may occur to further increase the regulatory and coding potential of the human genome. Among these, splicing of microexons, recursive splicing and biogenesis of circular and chimeric RNAs through back-splicing and trans-splicing processes, respectively, all contribute to expanding the repertoire of RNA transcripts with newly acquired regulatory functions. Interestingly, these non-canonical splicing events seem to occur more frequently in the central nervous system, affecting neuronal development and differentiation programs with important implications on brain physiology. Coherently, dysregulation of non-canonical RNA processing events is associated with brain disorders, including brain tumours. Herein, we summarize the current knowledge on molecular and regulatory mechanisms underlying canonical and non-canonical splicing events with particular emphasis on cis-acting elements and trans-acting factors that all together orchestrate splicing catalysis reactions and decisions. Lastly, we review the impact of non-canonical splicing on brain physiology and pathology and how unconventional splicing mechanisms may be targeted or exploited for novel therapeutic strategies in cancer.

Parallel evolution of a splicing program controlling neuronal excitability in flies and mammals

Alternative splicing increases neuronal transcriptomic complexity throughout animal phylogeny. To delve into the mechanisms controlling the assembly and evolution of this regulatory layer, we characterized the neuronal microexon program in Drosophila and compared it with that of mammals. In nonvertebrate bilaterians, this splicing program is restricted to neurons by the posttranscriptional processing of the enhancer of microexons (eMIC) domain in Srrm234. In Drosophila, this processing is dependent on regulation by Elav/Fne. eMIC deficiency or misexpression leads to widespread neurological alterations largely emerging from impaired neuronal activity, as revealed by a combination of neuronal imaging experiments and cell type-specific rescues. These defects are associated with the genome-wide skipping of short neural exons, which are strongly enriched in ion channels. We found no overlap of eMIC-regulated exons between flies and mice, illustrating how ancient posttranscriptional programs can evolve independently in different phyla to affect distinct cellular modules while maintaining cell-type specificity.

Histone demethylase Lsd1 is required for the differentiation of neural cells in Nematostella vectensis

Chromatin regulation is a key process in development but its contribution to the evolution of animals is largely unexplored. Chromatin is regulated by a diverse set of proteins, which themselves are tightly regulated in a cell/tissue-specific manner. Using the cnidarian Nematostella vectensis as a basal metazoan model, we explore the function of one such chromatin regulator, Lysine specific demethylase 1 (Lsd1). We generated an endogenously tagged allele and show that NvLsd1 expression is developmentally regulated and higher in differentiated neural cells than their progenitors. We further show, using a CRISPR/Cas9 generated mutant that loss of NvLsd1 leads to developmental abnormalities. This includes the almost complete loss of differentiated cnidocytes, cnidarian-specific neural cells, as a result of a cell-autonomous requirement for NvLsd1. Together this suggests that the integration of chromatin modifying proteins into developmental regulation predates the split of the cnidarian and bilaterian lineages and constitutes an ancient feature of animal development.

The evolutionary point where chromatin modifier function integrated into regulation of specific cell types is unclear. In the cnidarian Nematostella vectensis, the authors here show that lysine specific demethylase Lsd1 is developmentally regulated and required for normal development including cnidocyte differentiation.

Computational Analysis of Alternative Splicing Using VAST-TOOLS and the VastDB Framework

Alternative splicing (AS) can vastly expand animal transcriptomes and proteomes. Two main open questions in the field are how AS is regulated across cell/tissue types and disease, and what roles different AS events play. To facilitate AS research, we have created the computational VastDB framework, which comprises a series of complementary software and resources that we describe in this chapter. The VastDB framework is especially designed to aid biomedical researchers without a strong computational background. It offers tools and resources to: (a) quantify AS and identify differentially spliced AS events using RNA-seq data (vast-tools), (b) perform multiple genomic and sequence analyses for investigating AS events (Matt), (c) identify AS events with genomic and regulatory conservation among species (ExOrthist), and (d) help with the biological interpretation of the results, and, ultimately, with the identification of interesting AS events to design wet-lab experiments (VastDB and PastDB).

Dynamically expressed single ELAV/Hu orthologue elavl2 of bees is required for learning and memory

Changes in gene expression are a hallmark of learning and memory consolidation. Little is known about how alternative mRNA processing, particularly abundant in neuron-specific genes, contributes to these processes. Prototype RNA binding proteins of the neuronally expressed ELAV/Hu family are candidates for roles in learning and memory, but their capacity to cross-regulate and take over each other's functions complicate substantiation of such links. Honey bees Apis mellifera have only one elav/Hu family gene elavl2, that has functionally diversified by increasing alternative splicing including an evolutionary conserved microexon. RNAi knockdown demonstrates that ELAVL2 is required for learning and memory in bees. ELAVL2 is dynamically expressed with altered alternative splicing and subcellular localization in mushroom bodies, but not in other brain regions. Expression and alternative splicing of elavl2 change during memory consolidation illustrating an alternative mRNA processing program as part of a local gene expression response underlying memory consolidation.

ExOrthist: a tool to infer exon orthologies at any evolutionary distance

Several bioinformatic tools have been developed for genome-wide identification of orthologous and paralogous genes. However, no corresponding tool allows the detection of exon homology relationships. Here, we present ExOrthist, a fully reproducible Nextflow-based software enabling inference of exon homologs and orthogroups, visualization of evolution of exon-intron structures, and assessment of conservation of alternative splicing patterns. ExOrthist evaluates exon sequence conservation and considers the surrounding exon-intron context to derive genome-wide multi-species exon homologies at any evolutionary distance. We demonstrate its use in different evolutionary scenarios: whole genome duplication in frogs and convergence of Nova-regulated splicing networks (https://github.com/biocorecrg/ExOrthist).

Spliceosomal component PRP-40 is a central regulator of microexon splicing

Microexons (≤27 nt) play critical roles in nervous system development and function but create unique challenges for the splicing machinery. The mechanisms of microexon regulation are therefore of great interest. We performed a genetic screen for alternative splicing regulators in the C. elegans nervous system and identify PRP-40, a core component of the U1 snRNP. RNA-seq reveals that PRP-40 is required for inclusion of alternatively spliced, but not constitutively spliced, exons. PRP-40 is particularly required for inclusion of neuronal microexons, and our data indicate that PRP-40 is a central regulator of microexon splicing. Microexons can be relieved from PRP-40 dependence by artificially increasing exon size or reducing flanking intron size, indicating that PRP-40 is specifically required for microexons surrounded by conventionally sized introns. Knockdown of the orthologous PRPF40A in mouse neuroblastoma cells causes widespread dysregulation of microexons but not conventionally sized exons. PRP-40 regulation of neuronal microexons is therefore a widely conserved phenomenon.

Matrix-screening reveals a vast potential for direct protein-protein interactions among RNA binding proteins

RNA-binding proteins (RBPs) are crucial factors of post-transcriptional gene regulation and their modes of action are intensely investigated. At the center of attention are RNA motifs that guide where RBPs bind. However, sequence motifs are often poor predictors of RBP-RNA interactions in vivo. It is hence believed that many RBPs recognize RNAs as complexes, to increase specificity and regulatory possibilities. To probe the potential for complex formation among RBPs, we assembled a library of 978 mammalian RBPs and used rec-Y2H matrix screening to detect direct interactions between RBPs, sampling > 600 K interactions. We discovered 1994 new interactions and demonstrate that interacting RBPs bind RNAs adjacently in vivo. We further find that the mRNA binding region and motif preferences of RBPs deviate, depending on their adjacently binding interaction partners. Finally, we reveal novel RBP interaction networks among major RNA processing steps and show that splicing impairing RBP mutations observed in cancer rewire spliceosomal interaction networks. The dataset we provide will be a valuable resource for understanding the combinatorial interactions of RBPs with RNAs and the resulting regulatory outcomes.

RNA splicing factors SRRM3 and SRRM4 distinguish molecular phenotypes of castration-resistant neuroendocrine prostate cancer

Neuroendocrine (NE) differentiation in metastatic castration-resistant prostate cancer (mCRPC) is an increasingly common clinical feature arising from cellular plasticity. We recently characterized two mCRPC phenotypes with NE features: androgen receptor (AR)-positive NE-positive amphicrine prostate cancer (AMPC) and AR-negative small cell or neuroendocrine prostate cancer (SCNPC). Here, we interrogated the regulation of RE1-silencing transcription factor (REST), a transcriptional repressor of neuronal genes, and elucidated molecular programs driving AMPC and SCNPC biology. Analysis of prostate cancer (PC) cell lines, mCRPC specimens, and LuCaP patient-derived xenograft models detected alternative splicing of REST to REST4 and attenuated REST repressor activity in AMPC and SCNPC. The REST locus was also hypermethylated and REST expression was reduced in SCNPC. While serine/arginine repetitive matrix protein 4 (SRRM4) was previously implicated in alternative splicing of REST in mCRPC, we detected SRRM3 expression in REST4-positive, SRRM4-negative AMPC and SCNPC. In CRPC cell lines, SRRM3 induced alternative splicing of REST to REST4 and exacerbated the expression of REST-repressed genes. Furthermore, SRRM3 and SRRM4 expression defined molecular subsets of AMPC and SCNPC across species and tumor types. Two AMPC phenotypes and three SCNPC phenotypes were characterized, denoted either by REST attenuation and ASCL1 activity or by progressive activation of neuronal transcription factor programs, respectively. These results nominate SRRM3 as the principal REST splicing factor expressed in early NE differentiation and provide a framework to molecularly classify diverse NE phenotypes in mCRPC.

Microexon alternative splicing of small GTPase regulators: Implication in central nervous system diseases

Microexons are small sized (≤51 bp) exons which undergo extensive alternative splicing in neurons, microglia, embryonic stem cells, and cancer cells, giving rise to cell type specific protein isoforms. Due to their small sizes, microexons provide a unique challenge for the splicing machinery. They frequently lack exon splicer enhancers/repressors and require specialized neighboring trans-regulatory and cis-regulatory elements bound by RNA binding proteins (RBPs) for their inclusion. The functional consequences of including microexons within mRNAs have been extensively documented in the central nervous system (CNS) and aberrations in their inclusion have been observed to lead to abnormal processes. Despite the increasing evidence for microexons impacting cellular physiology within CNS, mechanistic details illustrating their functional importance in diseases of the CNS is still limited. In this review, we discuss the unique characteristics of microexons, and how RBPs participate in regulating their inclusion and exclusion during splicing. We consider recent findings of microexon alternative splicing and their implication for regulating the function of small GTPases in the context of the microglia, and we extrapolate these findings to what is known in neurons. We further discuss the emerging evidence for dysregulation of the Rho GTPase pathway in CNS diseases and the consequences contributed by the mis-splicing of microexons. This article is categorized under: RNA Processing > Splicing Mechanisms RNA Processing > Splicing Regulation/Alternative Splicing RNA in Disease and Development > RNA in Disease.

FINDER: an automated software package to annotate eukaryotic genes from RNA-Seq data and associated protein sequences

BACKGROUND

Gene annotation in eukaryotes is a non-trivial task that requires meticulous analysis of accumulated transcript data. Challenges include transcriptionally active regions of the genome that contain overlapping genes, genes that produce numerous transcripts, transposable elements and numerous diverse sequence repeats. Currently available gene annotation software applications depend on pre-constructed full-length gene sequence assemblies which are not guaranteed to be error-free. The origins of these sequences are often uncertain, making it difficult to identify and rectify errors in them. This hinders the creation of an accurate and holistic representation of the transcriptomic landscape across multiple tissue types and experimental conditions. Therefore, to gauge the extent of diversity in gene structures, a comprehensive analysis of genome-wide expression data is imperative.

RESULTS

We present FINDER, a fully automated computational tool that optimizes the entire process of annotating genes and transcript structures. Unlike current state-of-the-art pipelines, FINDER automates the RNA-Seq pre-processing step by working directly with raw sequence reads and optimizes gene prediction from BRAKER2 by supplementing these reads with associated proteins. The FINDER pipeline (1) reports transcripts and recognizes genes that are expressed under specific conditions, (2) generates all possible alternatively spliced transcripts from expressed RNA-Seq data, (3) analyzes read coverage patterns to modify existing transcript models and create new ones, and (4) scores genes as high- or low-confidence based on the available evidence across multiple datasets. We demonstrate the ability of FINDER to automatically annotate a diverse pool of genomes from eight species.

CONCLUSIONS

FINDER takes a completely automated approach to annotate genes directly from raw expression data. It is capable of processing eukaryotic genomes of all sizes and requires no manual supervision-ideal for bench researchers with limited experience in handling computational tools.

NvPOU4/Brain3 Functions as a Terminal Selector Gene in the Nervous System of the Cnidarian Nematostella vectensis

Terminal selectors are transcription factors that control the morphological, physiological, and molecular features that characterize distinct cell types. Here, we show that, in the sea anemone Nematostella vectensis, NvPOU4 is expressed in post-mitotic cells that give rise to a diverse set of neural cell types, including cnidocytes and NvElav1-expressing neurons. Morphological analyses of NvPOU4 mutants crossed to transgenic reporter lines show that the loss of NvPOU4 does not affect the initial specification of neural cells. Transcriptomes derived from the mutants and from different neural cell populations reveal that NvPOU4 is required for the execution of the terminal differentiation program of these neural cells. These findings suggest that POU4 genes have ancient functions as terminal selectors for morphologically and functionally disparate types of neurons and they provide experimental support for the relevance of terminal selectors for understanding the evolution of cell types.

The emergence of the brain non-CpG methylation system in vertebrates

Mammalian brains feature exceptionally high levels of non-CpG DNA methylation alongside the canonical form of CpG methylation. Non-CpG methylation plays a critical regulatory role in cognitive function, which is mediated by the binding of MeCP2, the transcriptional regulator that when mutated causes Rett syndrome. However, it is unclear whether the non-CpG neural methylation system is restricted to mammalian species with complex cognitive abilities or has deeper evolutionary origins. To test this, we investigated brain DNA methylation across 12 distantly related animal lineages, revealing that non-CpG methylation is restricted to vertebrates. We discovered that in vertebrates, non-CpG methylation is enriched within a highly conserved set of developmental genes transcriptionally repressed in adult brains, indicating that it demarcates a deeply conserved regulatory program. We also found that the writer of non-CpG methylation, DNMT3A, and the reader, MeCP2, originated at the onset of vertebrates as a result of the ancestral vertebrate whole-genome duplication. Together, we demonstrate how this novel layer of epigenetic information assembled at the root of vertebrates and gained new regulatory roles independent of the ancestral form of the canonical CpG methylation. This suggests that the emergence of non-CpG methylation may have fostered the evolution of sophisticated cognitive abilities found in the vertebrate lineage.

Silencing of SRRM4 suppresses microexon inclusion and promotes tumor growth across cancers

RNA splicing is widely dysregulated in cancer, frequently due to altered expression or activity of splicing factors (SFs). Microexons are extremely small exons (3–27 nucleotides long) that are highly evolutionarily conserved and play critical roles in promoting neuronal differentiation and development. Inclusion of microexons in mRNA transcripts is mediated by the SF Serine/Arginine Repetitive Matrix 4 (SRRM4), whose expression is largely restricted to neural tissues. However, microexons have been largely overlooked in prior analyses of splicing in cancer, as their small size necessitates specialized computational approaches for their detection. Here, we demonstrate that despite having low expression in normal nonneural tissues, SRRM4 is further silenced in tumors, resulting in the suppression of normal microexon inclusion. Remarkably, SRRM4 is the most consistently silenced SF across all tumor types analyzed, implying a general advantage of microexon down-regulation in cancer independent of its tissue of origin. We show that this silencing is favorable for tumor growth, as decreased SRRM4 expression in tumors is correlated with an increase in mitotic gene expression, and up-regulation of SRRM4 in cancer cell lines dose-dependently inhibits proliferation in vitro and in a mouse xenograft model. Further, this proliferation inhibition is accompanied by induction of neural-like expression and splicing patterns in cancer cells, suggesting that SRRM4 expression shifts the cell state away from proliferation and toward differentiation. We therefore conclude that SRRM4 acts as a proliferation brake, and tumors gain a selective advantage by cutting off this brake.

Using data from The Cancer Genome Atlas, this study shows that the splicing factor SRRM4 and its program of differentiation-promoting microexons are downregulated across tumor types with remarkable consistency, providing tumors with a proliferative advantage.

Alternative splicing landscapes in Arabidopsis thaliana across tissues and stress conditions highlight major functional differences with animals

BACKGROUND

Alternative splicing (AS) is a widespread regulatory mechanism in multicellular organisms. Numerous transcriptomic and single-gene studies in plants have investigated AS in response to specific conditions, especially environmental stress, unveiling substantial amounts of intron retention that modulate gene expression. However, a comprehensive study contrasting stress-response and tissue-specific AS patterns and directly comparing them with those of animal models is still missing.

RESULTS

We generate a massive resource for Arabidopsis thaliana, PastDB, comprising AS and gene expression quantifications across tissues, development and environmental conditions, including abiotic and biotic stresses. Harmonized analysis of these datasets reveals that A. thaliana shows high levels of AS, similar to fruitflies, and that, compared to animals, disproportionately uses AS for stress responses. We identify core sets of genes regulated specifically by either AS or transcription upon stresses or among tissues, a regulatory specialization that is tightly mirrored by the genomic features of these genes. Unexpectedly, non-intron retention events, including exon skipping, are overrepresented across regulated AS sets in A. thaliana, being also largely involved in modulating gene expression through NMD and uORF inclusion.

CONCLUSIONS

Non-intron retention events have likely been functionally underrated in plants. AS constitutes a distinct regulatory layer controlling gene expression upon internal and external stimuli whose target genes and master regulators are hardwired at the genomic level to specifically undergo post-transcriptional regulation. Given the higher relevance of AS in the response to different stresses when compared to animals, this molecular hardwiring is likely required for a proper environmental response in A. thaliana.

NineTeen Complex-subunit Salsa is required for efficient splicing of a subset of introns and dorsal-ventral patterning

The NineTeen Complex (NTC), also known as pre-mRNA-processing factor 19 (Prp19) complex, regulates distinct spliceosome conformational changes necessary for splicing. During Drosophila midblastula transition, splicing is particularly sensitive to mutations in NTC-subunit Fandango, which suggests differential requirements of NTC during development. We show that NTC-subunit Salsa, the Drosophila ortholog of human RNA helicase Aquarius, is rate-limiting for splicing of a subset of small first introns during oogenesis, including the first intron of gurken Germline depletion of Salsa and splice site mutations within gurken first intron impair both adult female fertility and oocyte dorsal-ventral patterning, due to an abnormal expression of Gurken. Supporting causality, the fertility and dorsal-ventral patterning defects observed after Salsa depletion could be suppressed by the expression of a gurken construct without its first intron. Altogether, our results suggest that one of the key rate-limiting functions of Salsa during oogenesis is to ensure the correct expression and efficient splicing of the first intron of gurken mRNA. Retention of gurken first intron compromises the function of this gene most likely because it undermines the correct structure and function of the transcript 5'UTR.

SRRM4 Expands the Repertoire of Circular RNAs by Regulating Microexon Inclusion

High-throughput RNA sequencing (RNA-seq) and dedicated bioinformatics pipelines have synergized to identify an expansive repertoire of unique circular RNAs (circRNAs), exceeding 100,000 variants. While the vast majority of these circRNAs comprise canonical exonic and intronic sequences, microexons (MEs)-which occur in 30% of functional mRNA transcripts-have been entirely overlooked. CircRNAs which contain these known MEs (ME-circRNAs) could be identified with commonly utilized circRNA prediction pipelines, CIRCexplorer2 and CIRI2, but were not previously recognized as ME-circRNAs. In addition, when employing a bespoke bioinformatics pipeline for identifying RNA chimeras, called Hyb, we could also identify over 2000 ME-circRNAs which contain novel MEs at their backsplice junctions, that are uncalled by either CIRCexplorer2 or CIRI2. Analysis of circRNA-seq datasets from gliomas of varying clinical grades compared with matched control tissue has shown circRNAs have potential as prognostic markers for stratifying tumor from healthy tissue. Furthermore, the abundance of microexon-containing circRNAs (ME-circRNAs) between tumor and normal tissues is correlated with the expression of a splicing associated factor, Serine/arginine repetitive matrix 4 (SRRM4). Overexpressing SRRM4, known for regulating ME inclusion in mRNAs critical for neural differentiation, in human HEK293 cells resulted in the biogenesis of over 2000 novel ME-circRNAs, including ME-circEIF4G3, and changes in the abundance of many canonical circRNAs, including circSETDB2 and circLBRA. This shows SRRM4, in which its expression is correlated with poor prognosis in gliomas, acts as a bona fide circRNA biogenesis factor. Given the known roles of MEs and circRNAs in oncogenesis, the identification of these previously unrecognized ME-circRNAs further increases the complexity and functional purview of this non-coding RNA family.

Global regulatory features of alternative splicing across tissues and within the nervous system of C. elegans

Alternative splicing plays a major role in shaping tissue-specific transcriptomes. Among the broad tissue types present in metazoans, the central nervous system contains some of the highest levels of alternative splicing. Although many documented examples of splicing differences between broad tissue types exist, there remains much to be understood about the splicing factors and the cis sequence elements controlling tissue and neuron subtype-specific splicing patterns. By using translating ribosome affinity purification coupled with deep-sequencing (TRAP-seq) in Caenorhabditis elegans, we have obtained high coverage profiles of ribosome-associated mRNA for three broad tissue classes (nervous system, muscle, and intestine) and two neuronal subtypes (dopaminergic and serotonergic neurons). We have identified hundreds of splice junctions that exhibit distinct splicing patterns between tissue types or within the nervous system. Alternative splicing events differentially regulated between tissues are more often frame-preserving, are more highly conserved across Caenorhabditis species, and are enriched in specific cis regulatory motifs, when compared with other types of exons. By using this information, we have identified a likely mechanism of splicing repression by the RNA-binding protein UNC-75/CELF via interactions with cis elements that overlap a 5′ splice site. Alternatively spliced exons also overlap more frequently with intrinsically disordered peptide regions than constitutive exons. Moreover, regulated exons are often shorter than constitutive exons but are flanked by longer intron sequences. Among these tissue-regulated exons are several highly conserved microexons <27 nt in length. Collectively, our results indicate a rich layer of tissue-specific gene regulation at the level of alternative splicing in C. elegans that parallels the evolutionary forces and constraints observed across metazoa.

Microexons: at the nexus of nervous system development, behaviour and autism spectrum disorder

The discovery and characterization of a network of highly conserved neuronal microexons has provided fundamental new insight into mechanisms underlying nervous system development and function, as well as an important basis for pathway convergence in autism spectrum disorder. In the past few years, considerable progress has been made in comprehensively determining the repertoires of factors that control neuronal microexons. These results have illuminated molecular mechanisms that activate the splicing of microexons, including those that control gene expression programs critical for neurogenesis, as well as synaptic protein translation and neuronal activity. Remarkably, individual disruption of specific microexons in these pathways results in autism-like phenotypes and cognitive impairment in mice. This review discusses these findings and their implications for delivering new therapeutic strategies for neurological disorders.

Alternative Splicing Regulatory Networks: Functions, Mechanisms, and Evolution

High-throughput sequencing-based methods and their applications in the study of transcriptomes have revolutionized our understanding of alternative splicing. Networks of functionally coordinated and biologically important alternative splicing events continue to be discovered in an ever-increasing diversity of cell types in the context of physiologically normal and disease states. These studies have been complemented by efforts directed at defining sequence codes governing splicing and their cognate trans-acting factors, which have illuminated important combinatorial principles of regulation. Additional studies have revealed critical roles of position-dependent, multivalent protein-RNA interactions that direct splicing outcomes. Investigations of evolutionary changes in RNA binding proteins, splice variants, and associated cis elements have further shed light on the emergence, mechanisms, and functions of splicing networks. Progress in these areas has emphasized the need for a coordinated, community-based effort to systematically address the functions of individual splice variants associated with normal and disease biology.

De novo identification of essential protein domains from CRISPR-Cas9 tiling-sgRNA knockout screens

High-throughput CRISPR-Cas9 knockout screens using a tiling-sgRNA design permit in situ evaluation of protein domain function. Here, to facilitate de novo identification of essential protein domains from such screens, we propose ProTiler, a computational method for the robust mapping of CRISPR knockout hyper-sensitive (CKHS) regions, which refer to the protein regions associated with a strong sgRNA dropout effect in the screens. Applied to a published CRISPR tiling screen dataset, ProTiler identifies 175 CKHS regions in 83 proteins. Of these CKHS regions, more than 80% overlap with annotated Pfam domains, including all of the 15 known drug targets in the dataset. ProTiler also reveals unannotated essential domains, including the N-terminus of the SWI/SNF subunit SMARCB1, which is validated experimentally. Surprisingly, the CKHS regions are negatively correlated with phosphorylation and acetylation sites, suggesting that protein domains and post-translational modification sites have distinct sensitivities to CRISPR-Cas9 mediated amino acids loss.

Srrm234, but not canonical SR and hnRNP proteins, drive inclusion of Dscam exon 9 variable exons

Alternative splicing of pre-mRNA is a major mechanism to diversify protein functionality in metazoans from a limited number of genes. The Drosophila melanogaster Down syndrome cell adhesion molecule (Dscam) gene, which is important for neuronal wiring and phagocytosis of bacteria, can generate up to 38,016 isoforms by mutually exclusive alternative splicing in four clusters of variable exons. However, it is not understood how a specific exon is chosen from the many variables and how variable exons are prevented from being spliced together. A main role in the regulation of Dscam alternative splicing has been attributed to RNA binding proteins (RBPs), but how they impact on exon selection is not well understood. Serine-arginine rich (SR) proteins and hnRNP proteins are the two main types of RBPs with major roles in exon definition and splice site selection. Here, we analyzed the role of SR and hnRNP proteins in Dscam exon 9 alternative splicing in mutant Drosophila melanogaster embryos because of their essential function for development. Strikingly, loss or overexpression of canonical SR and hnRNP proteins even when multiple proteins are depleted together, does not affect Dscam alternative exon selection very dramatically. Conversely, noncanonical SR protein Serine-arginine repetitive matrix 2/3/4 (Srrm234) is a main determinant of exon inclusion in the Dscam exon 9 cluster. Since long-range base-pairings are absent in the exon 9 cluster, our data argue for a small complement of regulatory factors as main determinants of exon inclusion in the Dscam exon 9 cluster.

The Protein-Coding Human Genome: Annotating High-Hanging Fruits

The major transcript variants of human protein-coding genes are annotated to a certain degree of accuracy combining manual curation, transcript data, and proteomics evidence. However, there is considerable disagreement on the annotation of about 2000 genes-they can be protein-coding, noncoding, or pseudogenes-and on the annotation of most of the predicted alternative transcripts. Pure transcriptome mapping approaches seem to be limited in discriminating functional expression from noise. These limitations have partially been overcome by dedicated algorithms to detect alternative spliced micro-exons and wobble splice variants. Recently, knowledge about splice mechanism and protein structure are incorporated into an algorithm to predict neighboring homologous exons, often spliced in a mutually exclusive manner. Predicted exons are evaluated by transcript data, structural compatibility, and evolutionary conservation, revealing hundreds of novel coding exons and splice mechanism re-assignments. The emerging human pan-genome is necessitating distinctive annotations incorporating differences between individuals and between populations.

Exciting times to study the identity and evolution of cell types

The EMBO/EMBL Symposium on ‘The Identity and Evolution of Cell Types’ took place in Heidelberg, Germany, on 15-19 May 2019. The symposium, which brought together a diverse group of speakers addressing a wide range of questions in multiple model systems, provided a platform to discuss how the concept of a cell type should be considered in the era of single cell omics techniques and how cell type evolution can be studied.