Trans-splicing enhances translational efficiency in C. elegans

U6 snRNA m6A modification is required for accurate and efficient splicing of C. elegans and human pre-mRNAs

pre-mRNA splicing is a critical feature of eukaryotic gene expression. Both cis- and trans-splicing rely on accurately recognising splice site sequences by spliceosomal U snRNAs and associated proteins. Spliceosomal snRNAs carry multiple RNA modifications with the potential to affect different stages of pre-mRNA splicing. Here, we show that the conserved U6 snRNA m6A methyltransferase METT-10 is required for accurate and efficient cis- and trans-splicing of C. elegans pre-mRNAs. The absence of METT-10 in C. elegans and METTL16 in humans primarily leads to alternative splicing at 5' splice sites with an adenosine at +4 position. In addition, METT-10 is required for splicing of weak 3' cis- and trans-splice sites. We identified a significant overlap between METT-10 and the conserved splicing factor SNRNP27K in regulating 5' splice sites with +4A. Finally, we show that editing endogenous 5' splice site +4A positions to +4U restores splicing to wild-type positions in a mett-10 mutant background, supporting a direct role for U6 snRNA m6A modification in 5' splice site recognition. We conclude that the U6 snRNA m6A modification is important for accurate and efficient pre-mRNA splicing.

A nematode-specific ribonucleoprotein complex mediates interactions between the major nematode spliced leader snRNP and its target pre-mRNAs

Spliced leader trans-splicing of pre-mRNAs is a critical step in the gene expression of many eukaryotes. How the spliced leader RNA and its target transcripts are brought together to form the trans-spliceosome remains an important unanswered question. Using immunoprecipitation followed by protein analysis via mass spectrometry and RIP-Seq, we show that the nematode-specific proteins, SNA-3 and SUT-1, form a complex with a set of enigmatic non-coding RNAs, the SmY RNAs. Our work redefines the SmY snRNP and shows for the first time that it is essential for nematode viability and is involved in spliced leader trans-splicing. SNA-3 and SUT-1 are associated with the 5' ends of most, if not all, nascent capped RNA polymerase II transcripts, and they also interact with components of the major nematode spliced leader (SL1) snRNP. We show that depletion of SNA-3 impairs the co-immunoprecipitation between one of the SL1 snRNP components, SNA-2, and several core spliceosomal proteins. We thus propose that the SmY snRNP recruits the SL1 snRNP to the 5' ends of nascent pre-mRNAs, an instrumental step in the assembly of the trans-spliceosome.

Molecular insights and antibody response to Dr20/22 in dogs naturally infected with Dirofilaria repens

Subcutaneous dirofilariasis, caused by the parasitic nematode Dirofilaria repens, is a growing concern in Europe, affecting both dogs and humans. This study focused on D. repens Dr20/22, a protein encoded by an alt (abundant larval transcript) gene family. While well-documented in L3 larvae of other filariae species, this gene family had not been explored in dirofilariasis. The research involved cloning Dr20/22 cDNA, molecular characterization, and evaluating its potential application in the diagnosis of dirofilariasis. Although Real-Time analysis revealed mRNA expression in both adult worms and microfilariae, the native protein remained undetected in lysates from both developmental stages. This suggests the protein's specificity for L3 larvae and may be related to a process called SLTS (spliced leader trans-splicing), contributing to stage-specific gene expression. The specificity of the antigen for invasive larvae positions it as a promising early marker for dirofilariasis. However, ELISA tests using sera from infected and uninfected dogs indicated limited diagnostic utility. While further research is required, our findings contribute to a deeper understanding of the molecular and immunological aspects of host-parasite interactions and could offer insights into the parasite's strategies for evading the immune system.

Characterization of trans-spliced chimeric RNAs: insights into the mechanism of trans-splicing

Trans-splicing is a post-transcriptional processing event that joins exons from separate RNAs to produce a chimeric RNA. However, the detailed mechanism of trans-splicing remains poorly understood. Here, we characterize trans-spliced genes and provide insights into the mechanism of trans-splicing in the tunicate Ciona. Tunicates are the closest invertebrates to humans, and their genes frequently undergo trans-splicing. Our analysis revealed that, in genes that give rise to both trans-spliced and non-trans-spliced messenger RNAs, trans-splice acceptor sites were preferentially located at the first functional acceptor site, and their paired donor sites were weak in both Ciona and humans. Additionally, we found that Ciona trans-spliced genes had GU- and AU-rich 5' transcribed regions. Our data and findings not only are useful for Ciona research community, but may also aid in a better understanding of the trans-splicing mechanism, potentially advancing the development of gene therapy based on trans-splicing.

eIF4E as a molecular wildcard in metazoans RNA metabolism

The evolutionary origin of eukaryotes spurred the transition from prokaryotic-like translation to a more sophisticated, eukaryotic translation. During this process, successive gene duplication of a single, primordial eIF4E gene encoding the mRNA cap-binding protein eukaryotic translation initiation factor 4E (eIF4E) gave rise to a plethora of paralog genes across eukaryotes that underwent further functional diversification in RNA metabolism. The ability to take different roles is due to eIF4E promiscuity in binding many partner proteins, rendering eIF4E a highly versatile and multifunctional player that functions as a molecular wildcard. Thus, in metazoans, eIF4E paralogs are involved in various processes, including messenger RNA (mRNA) processing, export, translation, storage, and decay. Moreover, some paralogs display differential expression in tissues and developmental stages and show variable biochemical properties. In this review, we discuss recent advances shedding light on the functional diversification of eIF4E in metazoans. We emphasise humans and two phylogenetically distant species which have become paradigms for studies on development, namely the fruit fly Drosophila melanogaster and the roundworm Caenorhabditis elegans.

Quantitative analysis of C. elegans transcripts by Nanopore direct-cDNA sequencing reveals terminal hairpins in non trans-spliced mRNAs

In nematodes and kinetoplastids, mRNA processing involves a trans-splicing step through which a short sequence from a snRNP replaces the original 5' end of the primary transcript. It has long been held that 70% of C. elegans mRNAs are submitted to trans-splicing. Our recent work suggested that the mechanism is more pervasive but not fully captured by mainstream transcriptome sequencing methods. Here we use Oxford Nanopore's long-read amplification-free sequencing technology to perform a comprehensive analysis of trans-splicing in worms. We demonstrate that spliced leader (SL) sequences at the 5' end of the mRNAs affect library preparation and generate sequencing artefacts due to their self-complementarity. Consistent with our previous observations, we find evidence of trans-splicing for most genes. However, a subset of genes appears to be only marginally trans-spliced. These mRNAs all share the capacity to generate a 5' terminal hairpin structure mimicking the SL structure and offering a mechanistic explanation for their non conformity. Altogether, our data provide a comprehensive quantitative analysis of SL usage in C. elegans.

Alternative Splicing in the Regulatory Circuit of Plant Temperature Response

As sessile organisms, plants have evolved complex mechanisms to rapidly respond to ever-changing ambient temperatures. Temperature response in plants is modulated by a multilayer regulatory network, including transcriptional and post-transcriptional regulations. Alternative splicing (AS) is an essential post-transcriptional regulatory mechanism. Extensive studies have confirmed its key role in plant temperature response, from adjustment to diurnal and seasonal temperature changes to response to extreme temperatures, which has been well documented by previous reviews. As a key node in the temperature response regulatory network, AS can be modulated by various upstream regulations, such as chromatin modification, transcription rate, RNA binding proteins, RNA structure and RNA modifications. Meanwhile, a number of downstream mechanisms are affected by AS, such as nonsense-mediated mRNA decay (NMD) pathway, translation efficiency and production of different protein variants. In this review, we focus on the links between splicing regulation and other mechanisms in plant temperature response. Recent advances regarding how AS is regulated and the following consequences in gene functional modulation in plant temperature response will be discussed. Substantial evidence suggests that a multilayer regulatory network integrating AS in plant temperature response has been unveiled.

slORFfinder: a tool to detect open reading frames resulting from trans-splicing of spliced leader sequences

Trans-splicing of a spliced leader (SL) to the 5' ends of mRNAs is used to produce mature mRNAs in several phyla of great importance to human health and the marine ecosystem. One of the consequences of the addition of SL sequences is the change or disruption of the open reading frames (ORFs) in the recipient transcripts. Given that most SL sequences have one or more of the trinucleotide NUG, including AUG in flatworms, trans-splicing of SL sequences can potentially supply a start codon to create new ORFs, which we refer to as slORFs, in the recipient mRNAs. Due to the lack of a tool to precisely detect them, slORFs were usually neglected in previous studies. In this work, we present the tool slORFfinder, which automatically links the SL sequences to the recipient mRNAs at the trans-splicing sites identified from SL-containing reads of RNA-Seq and predicts slORFs according to the distribution of ribosome-protected footprints (RPFs) on the trans-spliced transcripts. By applying this tool to the analyses of nematodes, ascidians and euglena, whose RPFs are publicly available, we find wide existence of slORFs in these taxa. Furthermore, we find that slORFs are generally translated at higher levels than the annotated ORFs in the genomes, suggesting they might have important functions. Overall, this study provides a tool, slORFfinder (https://github.com/songbo446/slORFfinder), to identify slORFs, which can enhance our understanding of ORFs in taxa with SL machinery.

Distance-dependent inhibition of translation initiation by downstream out-of-frame AUGs is consistent with a Brownian ratchet process of ribosome scanning

BACKGROUND

Eukaryotic ribosomes are widely presumed to scan mRNA for the AUG codon to initiate translation in a strictly 5'-3' movement (i.e., strictly unidirectional scanning model), so that ribosomes initiate translation exclusively at the 5' proximal AUG codon (i.e., the first-AUG rule).

RESULTS

We generate 13,437 yeast variants, each with an ATG triplet placed downstream (dATGs) of the annotated ATG (aATG) codon of a green fluorescent protein. We find that out-of-frame dATGs can inhibit translation at the aATG, but with diminishing strength over increasing distance between aATG and dATG, undetectable beyond ~17 nt. This phenomenon is best explained by a Brownian ratchet mechanism of ribosome scanning, in which the ribosome uses small-amplitude 5'-3' and 3'-5' oscillations with a net 5'-3' movement to scan the AUG codon, thereby leading to competition for translation initiation between aAUG and a proximal dAUG. This scanning model further predicts that the inhibitory effect induced by an out-of-frame upstream AUG triplet (uAUG) will diminish as uAUG approaches aAUG, which is indeed observed among the 15,586 uATG variants generated in this study. Computational simulations suggest that each triplet is scanned back and forth approximately ten times until the ribosome eventually migrates to downstream regions. Moreover, this scanning process could constrain the evolution of sequences downstream of the aATG to minimize proximal out-of-frame dATG triplets in yeast and humans.

CONCLUSIONS

Collectively, our findings uncover the basic process by which eukaryotic ribosomes scan for initiation codons, and how this process could shape eukaryotic genome evolution.

Transgenic Expression of Haemonchus contortus Cytochrome P450 Hco-cyp-13A11 Decreases Susceptibility to Particular but Not All Macrocyclic Lactones in the Model Organism Caenorhabditis elegans

The number of reported macrocyclic lactones (ML) resistance cases across all livestock hosts is steadily increasing. Different studies in the parasitic nematode Haemonchus contortus assume the participation of cytochrome P450s (Cyps) enzymes in ML resistance. Still, functional data about their individual contribution to resistance or substrate specificity is missing. Via microinjection, transgenic Caenorhabditis elegans expressing HCON_00141052 (transgene-Hco-cyp-13A11) from extrachromosomal arrays were generated. After 24 h of exposure to different concentrations of ivermectin (IVM), ivermectin aglycone (IVMa), selamectin (SEL), doramectin (DRM), eprinomectin (EPR), and moxidectin (MOX), motility assays were performed to determine the impact of the H. contortus Cyp to the susceptibility of the worms against each ML. While transgene-Hco-cyp-13A11 significantly decreased susceptibility to IVM (four-fold), IVMa (2-fold), and SEL (3-fold), a slight effect for DRM and no effect for MOX, and EPR was observed. This substrate specificity of Hco-cyp-13A11 could not be explained by molecular modeling and docking studies. Hco-Cyp-13A11 molecular models were obtained for alleles from isolates with different resistance statuses. Although 14 amino acid polymorphisms were detected, none was resistance specific. In conclusion, Hco-cyp-13A11 decreased IVM, IVMa, and SEL susceptibility to a different extent, but its potential impact on ML resistance is not driven by polymorphisms.

A novel, essential trans-splicing protein connects the nematode SL1 snRNP to the CBC-ARS2 complex

Spliced leader trans-splicing is essential for gene expression in many eukaryotes. To elucidate the molecular mechanism of this process, we characterise the molecules associated with the Caenorhabditis elegans major spliced leader snRNP (SL1 snRNP), which donates the spliced leader that replaces the 5' untranslated region of most pre-mRNAs. Using a GFP-tagged version of the SL1 snRNP protein SNA-1 created by CRISPR-mediated genome engineering, we immunoprecipitate and identify RNAs and protein components by RIP-Seq and mass spectrometry. This reveals the composition of the SL1 snRNP and identifies associations with spliceosome components PRP-8 and PRP-19. Significantly, we identify a novel, nematode-specific protein required for SL1 trans-splicing, which we designate SNA-3. SNA-3 is an essential, nuclear protein with three NADAR domains whose function is unknown. Mutation of key residues in NADAR domains inactivates the protein, indicating that domain function is required for activity. SNA-3 interacts with the CBC-ARS2 complex and other factors involved in RNA metabolism, including SUT-1 protein, through RNA or protein-mediated contacts revealed by yeast two-hybrid assays, localisation studies and immunoprecipitations. Our data are compatible with a role for SNA-3 in coordinating trans-splicing with target pre-mRNA transcription or in the processing of the Y-branch product of the trans-splicing reaction.

Efficient visual screening of CRISPR/Cas9 genome editing in the nematode Pristionchus pacificus

CRISPR/Cas9 genome editing has been applied to a wide variety of organisms, including nematodes such as Caenorhabditis elegans and Pristionchus pacificus. In these nematodes, genome editing is achieved by microinjection of Cas9 protein and guide RNA into the hermaphrodite gonads. However, P. pacificus is less efficient in CRISPR/Cas9 genome editing and exogenous gene expression. Therefore, it takes considerable time and effort to screen for target mutants if there are no visual markers that indicate successful injection. To overcome this problem, co-injection markers (gRNA for Ppa-prl-1, which induces the roller phenotype, and Ppa-egl-20p::turboRFP, a plasmid expressing a fluorescent protein) have been developed in P. pacificus. By selecting worms with the roller phenotype or turboRFP expression, screening efficiency is substantially increased to obtain worms with desired mutations. Here, we describe a step-by-step protocol for the visual screening system for CRISPR/Cas9 genome editing in P. pacificus. We also describe technical tips for microinjection, which is difficult for beginners. This protocol will facilitate genome editing in P. pacificus and may be applied to other nematode species.

Optimizing expression of a single copy transgene in C. elegans

The utility of single copy transgenic insertions in C. elegans is often limited by low expression. We examined the effects of modifying the trans-splicing signal, the Kozak ribosome binding site, the N-terminal amino acid of the reporter and the 3′ UTR sequences on the expression level of a mec-4 promoter GFP transgene. The trans-splicing signal and the 3′ UTR had most dramatic effects on expression while modifying the Kozak signal or the N-terminal amino acid had less influence on expression.

Determinants of genome-wide distribution and evolution of uORFs in eukaryotes

Upstream open reading frames (uORFs) play widespread regulatory functions in modulating mRNA translation in eukaryotes, but the principles underlying the genomic distribution and evolution of uORFs remain poorly understood. Here, we analyze ~17 million putative canonical uORFs in 478 eukaryotic species that span most of the extant taxa of eukaryotes. We demonstrate how positive and purifying selection, coupled with differences in effective population size (Ne), has shaped the contents of uORFs in eukaryotes. Besides, gene expression level is important in influencing uORF occurrences across genes in a species. Our analyses suggest that most uORFs might play regulatory roles rather than encode functional peptides. We also show that the Kozak sequence context of uORFs has evolved across eukaryotic clades, and that noncanonical uORFs tend to have weaker suppressive effects than canonical uORFs in translation regulation. This study provides insights into the driving forces underlying uORF evolution in eukaryotes.

Engineering rules that minimize germline silencing of transgenes in simple extrachromosomal arrays in C. elegans

Transgenes are prone to progressive silencing due to their structure, copy number, and genomic location. In C. elegans, repressive mechanisms are particularly strong in the germline with almost fully penetrant transgene silencing in simple extrachromosomal arrays and frequent silencing of single-copy transgene insertions. A class of non-coding DNA, Periodic An/Tn Clusters (PATCs) can prevent transgene-silencing in repressive chromatin or from small interfering RNAs (piRNAs). Here, we describe design rules (codon-optimization, intron and PATC inclusion, elevated temperature (25 °C), and vector backbone removal) for efficient germline expression from arrays in wildtype animals. We generate web-based tools to analyze PATCs and reagents for the convenient assembly of PATC-rich transgenes. An extensive collection of silencing resistant fluorescent proteins (e.g., gfp, mCherry, and tagBFP) can be used for dissecting germline regulatory elements and a set of enhanced enzymes (Mos1 transposase, Cas9, Cre, and Flp recombinases) enable efficient genetic engineering in C. elegans.

Improving Transgenesis Efficiency and CRISPR-Associated Tools Through Codon Optimization and Native Intron Addition in Pristionchus Nematodes

A lack of appropriate molecular tools is one obstacle that prevents in-depth mechanistic studies in many organisms. Transgenesis, clustered regularly interspaced short palindromic repeats (CRISPR)-associated engineering, and related tools are fundamental in the modern life sciences, but their applications are still limited to a few model organisms. In the phylum Nematoda, transgenesis can only be performed in a handful of species other than Caenorhabditis elegans, and additionally, other species suffer from significantly lower transgenesis efficiencies. We hypothesized that this may in part be due to incompatibilities of transgenes in the recipient organisms. Therefore, we investigated the genomic features of 10 nematode species from three of the major clades representing all different lifestyles. We found that these species show drastically different codon usage bias and intron composition. With these findings, we used the species Pristionchus pacificus as a proof of concept for codon optimization and native intron addition. Indeed, we were able to significantly improve transgenesis efficiency, a principle that may be usable in other nematode species. In addition, with the improved transgenes, we developed a fluorescent co-injection marker in P. pacificus for the detection of CRISPR-edited individuals, which helps considerably to reduce associated time and costs.

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.

Tissue-Specific Transcription Footprinting Using RNA PoI DamID (RAPID) in Caenorhabditis elegans

Differential gene expression across cell types underlies development and cell physiology in multicellular organisms. Caenorhabditis elegans is a powerful, extensively used model to address these biological questions. A remaining bottleneck relates to the difficulty to obtain comprehensive tissue-specific gene transcription data, since available methods are still challenging to execute and/or require large worm populations. Here, we introduce the RNA Polymerase DamID (RAPID) approach, in which the Dam methyltransferase is fused to a ubiquitous RNA polymerase subunit to create transcriptional footprints via methyl marks on the DNA of transcribed genes. To validate the method, we determined the polymerase footprints in whole animals, in sorted embryonic blastomeres and in different tissues from intact young adults by driving tissue-specific Dam fusion expression. We obtained meaningful transcriptional footprints in line with RNA-sequencing (RNA-seq) studies in whole animals or specific tissues. To challenge the sensitivity of RAPID and demonstrate its utility to determine novel tissue-specific transcriptional profiles, we determined the transcriptional footprints of the pair of XXX neuroendocrine cells, representing 0.2% of the somatic cell content of the animals. We identified 3901 candidate genes with putatively active transcription in XXX cells, including the few previously known markers for these cells. Using transcriptional reporters for a subset of new hits, we confirmed that the majority of them were expressed in XXX cells and identified novel XXX-specific markers. Taken together, our work establishes RAPID as a valid method for the determination of RNA polymerase footprints in specific tissues of C. elegans without the need for cell sorting or RNA tagging.

Alternative ORFs and small ORFs: shedding light on the dark proteome

Traditional annotation of protein-encoding genes relied on assumptions, such as one open reading frame (ORF) encodes one protein and minimal lengths for translated proteins. With the serendipitous discoveries of translated ORFs encoded upstream and downstream of annotated ORFs, from alternative start sites nested within annotated ORFs and from RNAs previously considered noncoding, it is becoming clear that these initial assumptions are incorrect. The findings have led to the realization that genetic information is more densely coded and that the proteome is more complex than previously anticipated. As such, interest in the identification and characterization of the previously ignored 'dark proteome' is increasing, though we note that research in eukaryotes and bacteria has largely progressed in isolation. To bridge this gap and illustrate exciting findings emerging from studies of the dark proteome, we highlight recent advances in both eukaryotic and bacterial cells. We discuss progress in the detection of alternative ORFs as well as in the understanding of functions and the regulation of their expression and posit questions for future work.

Trans-splicing of mRNAs links gene transcription to translational control regulated by mTOR

BACKGROUND

In phylogenetically diverse organisms, the 5' ends of a subset of mRNAs are trans-spliced with a spliced leader (SL) RNA. The functions of SL trans-splicing, however, remain largely enigmatic.

RESULTS

We quantified translation genome-wide in the marine chordate, Oikopleura dioica, under inhibition of mTOR, a central growth regulator. Translation of trans-spliced TOP mRNAs was suppressed, consistent with a role of the SL sequence in nutrient-dependent translational control of growth-related mRNAs. Under crowded, nutrient-limiting conditions, O. dioica continued to filter-feed, but arrested growth until favorable conditions returned. Upon release from unfavorable conditions, initial recovery was independent of nutrient-responsive, trans-spliced genes, suggesting animal density sensing as a first trigger for resumption of development.

CONCLUSION

Our results are consistent with a proposed role of trans-splicing in the coordinated translational down-regulation of nutrient-responsive genes under growth-limiting conditions.

Dietary restriction induces posttranscriptional regulation of longevity genes

Dietary restriction (DR) increases life span through adaptive changes in gene expression. To understand more about these changes, we analyzed the transcriptome and translatome of Caenorhabditis elegans subjected to DR. Transcription of muscle regulatory and structural genes increased, whereas increased expression of amino acid metabolism and neuropeptide signaling genes was controlled at the level of translation. Evaluation of posttranscriptional regulation identified putative roles for RNA-binding proteins, RNA editing, miRNA, alternative splicing, and nonsense-mediated decay in response to nutrient limitation. Using RNA interference, we discovered several differentially expressed genes that regulate life span. We also found a compensatory role for translational regulation, which offsets dampened expression of a large subset of transcriptionally down-regulated genes. Furthermore, 3' UTR editing and intron retention increase under DR and correlate with diminished translation, whereas trans-spliced genes are refractory to reduced translation efficiency compared with messages with the native 5' UTR. Finally, we find that smg-6 and smg-7, which are genes governing selection and turnover of nonsense-mediated decay targets, are required for increased life span under DR.

Function and Evolution of Upstream ORFs in Eukaryotes

There is growing interest in the role of translational regulation in cellular homeostasis during organismal development. Translation initiation is the rate-limiting step in mRNA translation and is central to translational regulation. Upstream open reading frames (uORFs) are regulatory elements that are prevalent in eukaryotic mRNAs. uORFs modulate the translation initiation rate of downstream coding sequences (CDSs) by sequestering ribosomes. Over the past several years, genome-wide studies have revealed the widespread regulatory functions of uORFs in different species in different biological contexts. Here, we review the current understanding of uORF-mediated translational regulation from the perspective of functional and evolutionary genomics and address remaining gaps that deserve further study.

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.

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

Background

The spliceosomal transfer of a short spliced leader (SL) RNA to an independent pre-mRNA molecule is called SL trans-splicing and is widespread in the nematode Caenorhabditis elegans. While RNA-sequencing (RNA-seq) data contain information on such events, properly documented methods to extract them are lacking.

Findings

To address this, we developed SL-quant, a fast and flexible pipeline that adapts to paired-end and single-end RNA-seq data and accurately quantifies SL trans-splicing events. It is designed to work downstream of read mapping and uses the reads left unmapped as primary input. Briefly, the SL sequences are identified with high specificity and are trimmed from the input reads, which are then remapped on the reference genome and quantified at the nucleotide position level (SL trans-splice sites) or at the gene level.

Conclusions

SL-quant completes within 10 minutes on a basic desktop computer for typical C. elegans RNA-seq datasets and can be applied to other species as well. Validating the method, the SL trans-splice sites identified display the expected consensus sequence, and the results of the gene-level quantification are predictive of the gene position within operons. We also compared SL-quant to a recently published SL-containing read identification strategy that was found to be more sensitive but less specific than SL-quant. Both methods are implemented as a bash script available under the MIT license [1]. Full instructions for its installation, usage, and adaptation to other organisms are provided.