SL-quant: a fast and flexible pipeline to quantify spliced leader trans-splicing events from RNA-seq data

SSUP-72/PINN-1 coordinates RNA-polymerase II 3' pausing and developmental gene expression in C. elegans

During exit from Caenorhabditis elegans (C. elegans) L1 developmental arrest, a network of growth- and developmental genes is activated, many of which are organized into operons where transcriptional termination is uncoupled from mRNA 3'-end processing. CDK-12-mediated Pol II CTD S2 phosphorylation enhances SL2 trans-splicing at downstream operonic genes, preventing premature termination and ensuring proper gene expression for developmental progression. Using a genetic screen, we identified the SSUP-72/PINN-1 module as a suppressor of defects induced by CDK-12 inhibition. Loss of SSUP-72/PINN-1 bypasses the requirement for CDK-12 in post-embryonic development. Genome-wide analyses reveal that SSUP-72, a CTD S5P phosphatase, affects Pol II 3' pausing and regulates intra-operon termination. Our findings establish SSUP-72/PINN-1 as a key regulator of Pol II dynamics, coordinating operonic gene expression and growth during C. elegans post-embryonic development.

Polycistronic Expression of Multi-Subunit Complexes in the Eukaryotic Environment: A Narrative Review

Protein complexes are involved in many vital biological processes. Therefore, researchers need these protein complexes for biochemical and biophysical studies. Several methods exist for expressing multi-subunit proteins in eukaryotic cells, such as 2A sequences, IRES, or intein. Nevertheless, each of these elements has several disadvantages that limit their usage. In this article, we suggest a new system for expressing multi-subunit proteins, which have several advantages over existing methods meanwhile it, lacks most of their disadvantages. Leishmania is a unicellular eukaryote and member of the Trypanosomatidae family. In the expression system of Leishmania, pre-long RNAs that contain several protein sequences transcribe. Then these long RNAs separate into mature mRNAs in the process named trans splicing. For producing multi-subunit protein, Leishmania transformed with a vector containing the sequences of all subunits. Therefore, those subunits translate and form the complex under eukaryotic cell conditions. The sequence of each protein must separate by the spatial sequence needed for trans splicing. Based on a Leishmania expression pattern, not only is it possible to produce the complexes with the correct structures and post-translational modifications, but also it is possible to overcome previous method problems.

SLIDR and SLOPPR: flexible identification of spliced leader trans-splicing and prediction of eukaryotic operons from RNA-Seq data

BACKGROUND

Spliced leader (SL) trans-splicing replaces the 5' end of pre-mRNAs with the spliced leader, an exon derived from a specialised non-coding RNA originating from elsewhere in the genome. This process is essential for resolving polycistronic pre-mRNAs produced by eukaryotic operons into monocistronic transcripts. SL trans-splicing and operons may have independently evolved multiple times throughout Eukarya, yet our understanding of these phenomena is limited to only a few well-characterised organisms, most notably C. elegans and trypanosomes. The primary barrier to systematic discovery and characterisation of SL trans-splicing and operons is the lack of computational tools for exploiting the surge of transcriptomic and genomic resources for a wide range of eukaryotes.

RESULTS

Here we present two novel pipelines that automate the discovery of SLs and the prediction of operons in eukaryotic genomes from RNA-Seq data. SLIDR assembles putative SLs from 5' read tails present after read alignment to a reference genome or transcriptome, which are then verified by interrogating corresponding SL RNA genes for sequence motifs expected in bona fide SL RNA molecules. SLOPPR identifies RNA-Seq reads that contain a given 5' SL sequence, quantifies genome-wide SL trans-splicing events and predicts operons via distinct patterns of SL trans-splicing events across adjacent genes. We tested both pipelines with organisms known to carry out SL trans-splicing and organise their genes into operons, and demonstrate that (1) SLIDR correctly detects expected SLs and often discovers novel SL variants; (2) SLOPPR correctly identifies functionally specialised SLs, correctly predicts known operons and detects plausible novel operons.

CONCLUSIONS

SLIDR and SLOPPR are flexible tools that will accelerate research into the evolutionary dynamics of SL trans-splicing and operons throughout Eukarya and improve gene discovery and annotation for a wide range of eukaryotic genomes. Both pipelines are implemented in Bash and R and are built upon readily available software commonly installed on most bioinformatics servers. Biological insight can be gleaned even from sparse, low-coverage datasets, implying that an untapped wealth of information can be retrieved from existing RNA-Seq datasets as well as from novel full-isoform sequencing protocols as they become more widely available.

RNA polymerase II CTD S2P is dispensable for embryogenesis but mediates exit from developmental diapause in C. elegans

Serine 2 phosphorylation (S2P) within the CTD of RNA polymerase II is considered a Cdk9/Cdk12-dependent mark required for 3'-end processing. However, the relevance of CTD S2P in metazoan development is unknown. We show that cdk-12 lesions or a full-length CTD S2A substitution results in an identical phenotype in Caenorhabditis elegans Embryogenesis occurs in the complete absence of S2P, but the hatched larvae arrest development, mimicking the diapause induced when hatching occurs in the absence of food. Genome-wide analyses indicate that when CTD S2P is inhibited, only a subset of growth-related genes is not properly expressed. These genes correspond to SL2 trans-spliced mRNAs located in position 2 and over within operons. We show that CDK-12 is required for maximal occupancy of cleavage stimulatory factor necessary for SL2 trans-splicing. We propose that CTD S2P functions as a gene-specific signaling mark ensuring the nutritional control of the C. elegans developmental program.

Resolution of polycistronic RNA by SL2 trans-splicing is a widely conserved nematode trait

Spliced leader trans-splicing is essential for the processing and translation of polycistronic RNAs generated by eukaryotic operons. In C. elegans, a specialized spliced leader, SL2, provides the 5' end for uncapped pre-mRNAs derived from polycistronic RNAs. Studies of other nematodes suggested that SL2-type trans-splicing is a relatively recent innovation, confined to Rhabditina, the clade containing C. elegans and its close relatives. Here we conduct a survey of transcriptome-wide spliced leader trans-splicing in Trichinella spiralis, a distant relative of C. elegans with a particularly diverse repertoire of 15 spliced leaders. By systematically comparing the genomic context of trans-splicing events for each spliced leader, we identified a subset of T. spiralis spliced leaders that are specifically used to process polycistronic RNAs-the first examples of SL2-type spliced leaders outside of Rhabditina. These T. spiralis spliced leader RNAs possess a perfectly conserved stem-loop motif previously shown to be essential for SL2-type trans-splicing in C. elegans We show that genes trans-spliced to these SL2-type spliced leaders are organized in operonic fashion, with short intercistronic distances. A subset of T. spiralis operons show conservation of synteny with C. elegans operons. Our work substantially revises our understanding of nematode spliced leader trans-splicing, showing that SL2 trans-splicing is a major mechanism for nematode polycistronic RNA processing, which may have evolved prior to the radiation of the Nematoda. This work has important implications for the improvement of genome annotation pipelines in nematodes and other eukaryotes with operonic gene organization.

RNA structures in alternative splicing and back-splicing

Alternative splicing greatly expands the transcriptomic and proteomic diversities related to physiological and developmental processes in higher eukaryotes. Splicing of long noncoding RNAs, and back- and trans- splicing further expanded the regulatory repertoire of alternative splicing. RNA structures were shown to play an important role in regulating alternative splicing and back-splicing. Application of novel sequencing technologies made it possible to identify genome-wide RNA structures and interaction networks, which might provide new insights into RNA splicing regulation in vitro to in vivo. The emerging transcription-folding-splicing paradigm is changing our understanding of RNA alternative splicing regulation. Here, we review the insights into the roles and mechanisms of RNA structures in alternative splicing and back-splicing, as well as how disruption of these structures affects alternative splicing and then leads to human diseases. This article is categorized under: RNA Processing > Splicing Regulation/Alternative Splicing RNA Structure and Dynamics > Influence of RNA Structure in Biological Systems.

SLFinder, a pipeline for the novel identification of splice-leader sequences: a good enough solution for a complex problem

BACKGROUND

Spliced Leader trans-splicing is an important mechanism for the maturation of mRNAs in several lineages of eukaryotes, including several groups of parasites of great medical and economic importance. Nevertheless, its study across the tree of life is severely hindered by the problem of identifying the SL sequences that are being trans-spliced.

RESULTS

In this paper we present SLFinder, a four-step pipeline meant to identify de novo candidate SL sequences making very few assumptions regarding the SL sequence properties. The pipeline takes transcriptomic de novo assemblies and a reference genome as input and allows the user intervention on several points to account for unexpected features of the dataset. The strategy and its implementation were tested on real RNAseq data from species with and without SL Trans-Splicing.

CONCLUSIONS

SLFinder is capable to identify SL candidates with good precision in a reasonable amount of time. It is especially suitable for species with unknown SL sequences, generating candidate sequences for further refining and experimental validation.

Heterodera glycines utilizes promiscuous spliced leaders and demonstrates a unique preference for a species-specific spliced leader over C. elegans SL1

Spliced leader trans-splicing (SLTS) plays a part in the maturation of pre-mRNAs in select species across multiple phyla but is particularly prevalent in Nematoda. The role of spliced leaders (SL) within the cell is unclear and an accurate assessment of SL occurrence within an organism is possible only after extensive sequencing data are available, which is not currently the case for many nematode species. SL discovery is further complicated by an absence of SL sequences from high-throughput sequencing results due to incomplete sequencing of the 5'-ends of transcripts during RNA-seq library preparation, known as 5'-bias. Existing datasets and novel methodology were used to identify both conserved SLs and unique hypervariable SLs within Heterodera glycines, the soybean cyst nematode. In H. glycines, twenty-one distinct SL sequences were found on 2,532 unique H. glycines transcripts. The SL sequences identified on the H. glycines transcripts demonstrated a high level of promiscuity, meaning that some transcripts produced as many as nine different individual SL-transcript combinations. Most uniquely, transcriptome analysis revealed that H. glycines is the first nematode to demonstrate a higher SL trans-splicing rate using a species-specific SL over well-conserved Caenorhabditis elegans SL-like sequences.