Regulation of Tissue-Specific Alternative Splicing: C. elegans as a Model System

Sensory neuron transcriptomes reveal complex neuron-specific function and regulation of mec-2/Stomatin splicing

The function and identity of a cell is shaped by transcription factors controlling transcriptional networks, and further shaped by RNA binding proteins controlling post-transcriptional networks. To overcome limitations inherent to analysis of sparse single-cell post-transcriptional data, we leverage the invariant Caenorhabditis elegans cell lineage, isolating thousands of identical neuron types from thousands of isogenic individuals. The resulting deep transcriptomes facilitate splicing network analysis due to increased sequencing depth and uniformity. We focus on mechanosensory touch-neuron splicing regulated by MEC-8/RBPMS. We identify a small MEC-8-regulated network, where MEC-8 establishes touch-neuron isoforms differing from default isoforms found in other cells. MEC-8 establishes the canonical long mec-2/Stomatin isoform in touch neurons, but surprisingly the non-canonical short isoform predominates in other neurons, including olfactory neurons, and mec-2 is required for olfaction. Forced endogenous isoform-specific expression reveals that the short isoform functions in olfaction but not mechanosensation. The long isoform is functional in both processes. Remarkably, restoring the long isoform completely rescues mec-8 mutant mechanosensation, indicating a single MEC-8 touch-neuron target is phenotypically relevant. Within the long isoform we identify a cassette exon further diversifying mec-2 into long/extra-long isoforms. Neither is sufficient for mechanosensation. Both are simultaneously required, likely functioning as heteromers to mediate mechanosensation.

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.

RBFOX2 alters splicing outcome in distinct binding modes with multiple protein partners

RBFOX2 controls the splicing of a large number of transcripts implicated in cell differentiation and development. Parsing RNA-binding protein datasets, we uncover that RBFOX2 can interact with hnRNPC, hnRNPM and SRSF1 to regulate splicing of a broad range of splicing events using different sequence motifs and binding modes. Using immunoprecipitation, specific RBP knockdown, RNA-seq and splice-sensitive PCR, we show that RBFOX2 can target splice sites using three binding configurations: single, multiple or secondary modes. In the single binding mode RBFOX2 is recruited to its target splice sites through a single canonical binding motif, while in the multiple binding mode RBFOX2 binding sites include the adjacent binding of at least one other RNA binding protein partner. Finally, in the secondary binding mode RBFOX2 likely does not bind the RNA directly but is recruited to splice sites lacking its canonical binding motif through the binding of one of its protein partners. These dynamic modes bind distinct sets of transcripts at different positions and distances relative to alternative splice sites explaining the heterogeneity of RBFOX2 targets and splicing outcomes.

Alternative Splicing Dynamics of the Hypothalamus-Pituitary-Ovary Axis During Pubertal Transition in Gilts

The timing of puberty in mammals marks the point at which reproduction becomes possible. Abnormalities in the timing of puberty may exert a series of negative effects on subsequent health outcomes. Alternative splicing (AS) has not only emerged as a significant factor in the transcription of genes but it is also reported to play a role in the timing of puberty. However, to date, the changes and dynamics of AS during the onset of puberty is extremely seldom explored. In the present study, we used gilts as a research model to investigated the dynamics of AS and differentially expressed AS (DEAS) events within the hypothalamus-pituitary-ovary (HPO) axis across pre-, in-, and post-puberty. We detected 3,390, 6,098, and 9,085 DEAS events in the hypothalamus, pituitary, and ovary when compared across pre-, in-, and post-pubertal stages, respectively. Within the entire HPO axis, we also identified 22,889, 22,857, and 21,055 DEAS events in the pre-, in-, and post-pubertal stages, respectively. Further analysis revealed that the differentially spliced genes (DSGs) associated with staged DEAS events were likely to be enriched in the oxytocin signaling pathway, thyroid hormone signaling pathway, GnRH signaling pathway, and oocyte meiosis signaling pathway. The DSGs associated with DEAS events across the entire HPO axis were enriched in endocytosis signaling pathway, the MAPK signaling pathway, and the Rap1 signaling pathway. Moreover. the ASs of TAC1, TACR3, CYP19A1, ESR1, ESRRA, and FSHR were likely to regulate the functions of the certain HPO tissues during the onset of puberty. Collectively, the AS dynamics and DEAS events were comprehensively profiled in hypothalamus, pituitary, and ovary across the pre-, in-, and post-pubertal stages in pigs. These findings may enhance our knowledge of how puberty is regulated by AS and shed new light on the molecular mechanisms underlying the timing of puberty in mammals.

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.

mRNA Editing, Processing and Quality Control in Caenorhabditis elegans

While DNA serves as the blueprint of life, the distinct functions of each cell are determined by the dynamic expression of genes from the static genome. The amount and specific sequences of RNAs expressed in a given cell involves a number of regulated processes including RNA synthesis (transcription), processing, splicing, modification, polyadenylation, stability, translation, and degradation. As errors during mRNA production can create gene products that are deleterious to the organism, quality control mechanisms exist to survey and remove errors in mRNA expression and processing. Here, we will provide an overview of mRNA processing and quality control mechanisms that occur in Caenorhabditis elegans, with a focus on those that occur on protein-coding genes after transcription initiation. In addition, we will describe the genetic and technical approaches that have allowed studies in C. elegans to reveal important mechanistic insight into these processes.

Splicing in a single neuron is coordinately controlled by RNA binding proteins and transcription factors

Single-cell transcriptomes are established by transcription factors (TFs), which determine a cell's gene-expression complement. Post-transcriptional regulation of single-cell transcriptomes, and the RNA binding proteins (RBPs) responsible, are more technically challenging to determine, and combinatorial TF-RBP coordination of single-cell transcriptomes remains unexplored. We used fluorescent reporters to visualize alternative splicing in single Caenorhabditis elegans neurons, identifying complex splicing patterns in the neuronal kinase sad-1. Most neurons express both isoforms, but the ALM mechanosensory neuron expresses only the exon-included isoform, while its developmental sister cell the BDU neuron expresses only the exon-skipped isoform. A cascade of three cell-specific TFs and two RBPs are combinatorially required for sad-1 exon inclusion. Mechanistically, TFs combinatorially ensure expression of RBPs, which interact with sad-1 pre-mRNA. Thus a combinatorial TF-RBP code controls single-neuron sad-1 splicing. Additionally, we find ‘phenotypic convergence,’ previously observed for TFs, also applies to RBPs: different RBP combinations generate similar splicing outcomes in different neurons.

All the cells in the human nervous system contain the same genetic information, and yet there are many kinds of neurons, each with different features and roles in the body. Proteins known as transcription factors help to establish this diversity by switching on different genes in different types of cells.

A mechanism known as RNA splicing, which is regulated by RNA binding proteins, can also provide another layer of regulation. When a gene is switched on, a faithful copy of its sequence is produced in the form of an RNA molecule, which will then be ‘read’ to create a protein. However, the RNA molecules may first be processed to create templates that can differ between cell types: this means that a single gene can code for slightly different proteins, some of them specific to a given cell type. Yet, very little is known about how RNA splicing can generate more diversity in the nervous system.

To investigate, Thompson et al. developed a fluorescent reporter system that helped them track how the RNA of a gene called sad-1 is spliced in individual neurons of the worm Caenorhabditis elegans. This showed that sad-1 was turned on in all neurons, but the particular spliced versions varied widely between different types of nerve cells.

Additional experiments combined old school and cutting-edge genetics technics such as CRISPR/Cas9 to identify the proteins that control the splicing of sad-1 in different kinds of neurons. Despite not directly participating in RNA splicing, a number of transcription factors were shown to be involved. These molecular switches were turning on genes that code for RNA binding proteins differently between types of neurons, which in turn led sad-1 to be spliced according to neuron-specific patterns.

The findings by Thompson et al. could provide some insight into how mammals can establish many types of neurons; however, a technical hurdle stands in the way of this line of research, as it is still difficult to detect splicing in single neurons in these species.

Synthetic Genetic Interaction (CRISPR-SGI) Profiling in Caenorhabditis elegans

Genetic interaction screens are a powerful methodology to establish novel roles for genes and elucidate functional connections between genes. Such studies have been performed to great effect in single-cell organisms such as yeast and E. coli (Schuldiner et al., 2005; Butland et al., 2008; Costanzo et al., 2010), but similar large-scale interaction studies using targeted reverse-genetic deletions in multi-cellular organisms have not been feasible. We developed a CRISPR/Cas9-based method for deleting genes in C. elegans and replacing them with a heterologous fluorescent reporter (Norris et al., 2015). Recently we took advantage of that system to perform a large-scale, reverse genetic screen using null alleles in animals for the first time, focusing on RNA binding protein genes (Norris et al., 2017). This type of approach should be similarly applicable to many other gene classes in C. elegans. Here we detail the protocols involved in generating a library of double mutants and performing medium-throughput competitive fitness assays to test for genetic interactions resulting in fitness changes.

The combinatorial control of alternative splicing in C. elegans

Normal development requires the right splice variants to be made in the right tissues at the right time. The core splicing machinery is engaged in all splicing events, but which precise splice variant is made requires the choice between alternative splice sites—for this to occur, a set of splicing factors (SFs) must recognize and bind to short RNA motifs in the pre-mRNA. In C. elegans, there is known to be extensive variation in splicing patterns across development, but little is known about the targets of each SF or how multiple SFs combine to regulate splicing. Here we combine RNA-seq with in vitro binding assays to study how 4 different C. elegans SFs, ASD-1, FOX-1, MEC-8, and EXC-7, regulate splicing. The 4 SFs chosen all have well-characterised biology and well-studied loss-of-function genetic alleles, and all contain RRM domains. Intriguingly, while the SFs we examined have varied roles in C. elegans development, they show an unexpectedly high overlap in their targets. We also find that binding sites for these SFs occur on the same pre-mRNAs more frequently than expected suggesting extensive combinatorial control of splicing. We confirm that regulation of splicing by multiple SFs is often combinatorial and show that this is functionally significant. We also find that SFs appear to combine to affect splicing in two modes—they either bind in close proximity within the same intron or they appear to bind to separate regions of the intron in a conserved order. Finally, we find that the genes whose splicing are regulated by multiple SFs are highly enriched for genes involved in the cytoskeleton and in ion channels that are key for neurotransmission. Together, this shows that specific classes of genes have complex combinatorial regulation of splicing and that this combinatorial regulation is critical for normal development to occur.

Alternative splicing (AS) is a highly regulated process that is crucial for normal development. It requires the core splicing machinery, but the specific choice of splice site during AS is controlled by splicing factors (SFs) such as ELAV or RBFOX proteins that bind to specific sequences in pre-mRNAs to regulate usage of different splice sites. AS varies across the C. elegans life cycle and here we study how diverse SFs combine to regulate AS during C. elegans development. We selected 4 RRM-containing SFs that are all well studied and that have well-characterised loss-of-function genetic alleles. We find that these SFs regulate many of the same targets, and that combinatorial interactions between these SFs affect both individual splicing events and organism-level phenotypes including specific effects on the neuromuscular system. We further show that SFs combine to regulate splicing of an individual pre-mRNA in two distinct modes—either by binding in close proximity or by binding in a defined order on the pre-mRNA. Finally, we find that the genes whose splicing are most likely to be regulated by multiple SFs are genes that are required for the proper function of the neuromuscular system. These genes are also most likely to have changing AS patterns across development, suggesting that their splicing regulation is highly complex and developmentally regulated. Taken together, our data show that the precise splice variant expressed at any point in development is often the outcome of regulation by multiple SFs.

An emerging model organism Caenorhabditis elegans for alternative pre-mRNA processing in vivo

A nematode Caenorhabditis elegans is an intron-rich organism and up to 25% of its pre-mRNAs are estimated to be alternatively processed. Its compact genomic organization enables construction of fluorescence splicing reporters with intact genomic sequences and visualization of alternative processing patterns of interest in the transparent living animals with single-cell resolution. Genetic analysis with the reporter worms facilitated identification of trans-acting factors and cis-acting elements, which are highly conserved in mammals. Analysis of unspliced and partially spliced pre-mRNAs in vivo raised models for alternative splicing regulation relying on specific order of intron excision. RNA-seq analysis of splicing factor mutants and CLIP-seq analysis of the factors allow global search for target genes in the whole animal. An mRNA surveillance system is not essential for its viability or fertility, allowing analysis of unproductively spliced noncoding mRNAs. These features offer C. elegans as an ideal model organism for elucidating alternative pre-mRNA processing mechanisms in vivo. Examples of isoform-specific functions of alternatively processed genes are summarized. WIREs RNA 2017, 8:e1428. doi: 10.1002/wrna.1428 For further resources related to this article, please visit the WIREs website.

Caenorhabditis elegans PRMT-7 and PRMT-9 Are Evolutionarily Conserved Protein Arginine Methyltransferases with Distinct Substrate Specificities

Caenorhabditis elegans protein arginine methyltransferases PRMT-7 and PRMT-9 are two evolutionarily conserved enzymes, with distinct orthologs in plants, invertebrates, and vertebrates. Biochemical characterization of these two enzymes reveals that they share much in common with their mammalian orthologs. C. elegans PRMT-7 produces only monomethylarginine (MMA) and preferentially methylates R-X-R motifs in a broad collection of substrates, including human histone peptides and RG-rich peptides. In addition, the activity of the PRMT-7 enzyme is dependent on temperature, the presence of metal ions, and the reducing agent dithiothreitol. C. elegans PRMT-7 has a substrate specificity and a substrate preference different from those of mammalian PRMT7, and the available X-ray crystal structures of the PRMT7 orthologs show differences in active site architecture. C. elegans PRMT-9, on the other hand, produces symmetric dimethylarginine and MMA on SFTB-2, the conserved C. elegans ortholog of human RNA splicing factor SF3B2, indicating a possible role in the regulation of nematode splicing. In contrast to PRMT-7, C. elegans PRMT-9 appears to be biochemically indistinguishable from its human ortholog.