Copepods are one of the most abundant metazoans in the marine ecosystem, constituting a critical link in aquatic food webs and contributing significantly to the global carbon budget, yet molecular mechanisms of their gene expression are not well understood. Here we report the detection of spliced leader (SL) trans-splicing in calanoid copepods. We have examined nine species of wild-caught copepods from Jiaozhou Bay, China that represent the major families of the calanoids. All these species contained a common 46-nt SL (CopepodSL). We further determined the size of CopepodSL precursor RNA (slRNA; 108-158 nt) through genomic analysis and 3'-RACE technique, which was confirmed by RNA blot analysis. Structure modeling showed that the copepod slRNA folded into typical slRNA secondary structures. Using a CopepodSL-based primer set, we selectively enriched and sequenced copepod full-length cDNAs, which led to the characterization of copepod transcripts and the cataloging of the complete set of 79 eukaryotic cytoplasmic ribosomal proteins (cRPs) for a single copepod species. We uncovered the SL trans-splicing in copepod natural populations, and demonstrated that CopepodSL was a sensitive and specific tool for copepod transcriptomic studies at both the individual and population levels and that it would be useful for metatranscriptomic analysis of copepods.

High throughput sequencing techniques are poorly adapted for in vivo studies of parasites, which require prior in vitro culturing and purification. Trypanosomatids, a group of kinetoplastid protozoans, possess a distinctive feature in their transcriptional mechanism whereby a specific Spliced Leader (SL) sequence is added to the 5'end of each mRNA by trans-splicing. This allows to discriminate Trypansomatid RNA from mammalian RNA and forms the basis of our new multiplexed protocol for high-throughput, selective RNA-sequencing called SL-seq. We provided a proof-of-concept of SL-seq in Leishmania donovani, the main causative agent of visceral leishmaniasis in humans, and successfully applied the method to sequence Leishmania mRNA directly from infected macrophages and from highly diluted mixes with human RNA. mRNA profiles obtained with SL-seq corresponded largely to those obtained from conventional poly-A tail purification methods, indicating both enumerate the same mRNA pool. However, SL-seq offers additional advantages, including lower sequencing depth requirements, fast and simple library prep and high resolution splice site detection. SL-seq is therefore ideal for fast and massive parallel sequencing of parasite transcriptomes directly from host tissues. Since SLs are also present in Nematodes, Cnidaria and primitive chordates, this method could also have high potential for transcriptomics studies in other organisms.

The lack of general class II transcription factors was a hallmark of the genomic sequences of the human parasites Trypanosoma brucei, Trypanosoma cruzi and Leishmania major. However, the recent identification of TFIIA as part of a protein complex essential for RNA polymerase II-mediated transcription of SLRNA genes, which encode the trans splicing-specific spliced leader RNA, suggests that trypanosomatids assemble a highly divergent set of these factors at the SLRNA promoter. Here we report the identification of a trypanosomatid TFIIB-like (TFIIB(like)) protein which has limited overall sequence homology to eukaryotic TFIIB and archaeal TFB but harbors conserved residues within the N-terminal zinc ribbon domain, the B finger and cyclin repeat I. In accordance with the function of TFIIB, T.brucei TFIIB(like) is encoded by an essential gene, localizes to the nucleus, specifically binds to the SLRNA promoter, interacts with RNA polymerase II, and is absolutely required for SLRNA transcription.

Trans-splicing generates the mature 5' ends of certain mRNAs through the addition of a small spliced leader (SL) exon to pre-mRNAs. To search for novel flatworm spliced leaders, degenerate oligonucleotides and 5' RACE [corrected was used to isolate and characterize the 5' terminal sequences of enolase mRNAs in diverse flatworms. Several new spliced leaders and their SL RNA genes were identified, characterized, and compared. All parasitic trematodes examined trans-splice enolase. A primitive polyclad turbellarian, Stylochus zebra, also contains a trans-spliced enolase mRNA. The S. zebra SL is the longest SL yet identified, 51 nucleotides. Comparison of flatworm SLs indicates that they vary significantly in sequence and length. This suggests that neither spliced leader exon sequence nor size is likely to be essential for trans-splicing in flatworms. Flatworm SL RNAs have unusual Sm binding sites with characteristics distinct from other known flatworm snRNA Sm binding sites. Predicted flatworm SL RNA secondary structures show variation exhibiting 2-4 stem loops. Although limited in sequence similarity, phylogenetically conserved regions within the diverse flatworm SL RNAs suggest that they are likely to be derived from a common ancestor and provide information on potentially important SL RNA elements. The identification of a SL in a primitive flatworm suggests that trans-splicing may have been an ancestral feature in the phylum. Representative species of other early and more recent clades within the phylum, however, do not trans-splice enolase, nor do they or representatives of several other flatworm groups, have an SL RNA with a phylogenetically conserved region identified in the current study.

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.

Spliced leader (SL) trans-splicing is a critical element of gene expression in a number of eukaryotic groups. This process is arguably best understood in nematodes, where biochemical and molecular studies in Caenorhabditis elegans and Ascaris suum have identified key steps and factors involved. Despite this, the precise details of SL trans-splicing have yet to be elucidated. In part, this is because the systematic identification of the molecules involved has not previously been possible due to the lack of a specific phenotype associated with defects in this process. We present here a novel GFP-based reporter assay that can monitor SL1 trans-splicing in living C. elegans. Using this assay, we have identified mutants in sna-1 that are defective in SL trans-splicing, and demonstrate that reducing function of SNA-1, SNA-2 and SUT-1, proteins that associate with SL1 RNA and related SmY RNAs, impairs SL trans-splicing. We further demonstrate that the Sm proteins and pICln, SMN and Gemin5, which are involved in small nuclear ribonucleoprotein assembly, have an important role in SL trans-splicing. Taken together these results provide the first in vivo evidence for proteins involved in SL trans-splicing, and indicate that continuous replacement of SL ribonucleoproteins consumed during trans-splicing reactions is essential for effective trans-splicing.

Diverse mRNAs of Trypanosoma brucei possess the same 5' terminal 35 nucleotides, termed the spliced leader (SL), which appears to be derived from a separate 135 nucleotide transcript. This small SL RNA is encoded within a 1.4 kb unit of DNA which is tandemly reiterated in the genome. In addition, there are at least 4 orphon elements containing SL sequences dispersed from the tandem array. Here we show that during the trypanosome life cycle one of the SL orphons undergoes a stage-specific modification that prevents cleavage of an EcoRV site and we further demonstrate that although only one orphon is modified, three of the SL orphons are flanked by very similar sequences. Each of these contains SL reiteration units including the non-transcribed spacer DNA, suggesting that they did not originate through an RNA intermediate. In addition no evidence of direct repeats at the junction of 1.4 kb and non-1.4 kb DNA was observed. Finally, a phylogenetic survey indicates that while many trypanosomatid species possess similarly organized SL-like sequences, only the SL orphons of closely related subspecies of the T. brucei - T. evansi complex share similar flanking regions.

We have assessed the potential of using the spliced leader (SL) or mini-exon gene as a marker for molecular phylogenetic analysis of genus Trypanosoma. A total of 27 trypanosome sequences were compared, 18 of these being newly reported. In contrast to genus Leishmania, we found the non-transcribed spacer region of the SL locus in trypanosomes to be far too variable for informative comparison of all but the most closely related species. At the other extreme, the short (39 nt) SL exon was usually completely conserved and hence uninformative. The SL RNA showed variation in both length (97-152 nt) and sequence among different trypanosome species, with most variation occurring in stem-loop II. Consequently, this region could not be aligned with confidence in multiple sequence alignment, severely reducing the number of phylogenetically informative nucleotide positions. In computer simulation, most of the SL RNAs readily folded into the 3 stem-loop secondary structure predicted previously, but again stem-loop II was highly variable. No obvious correlation could be seen between the length of this stem-loop and trypanosome biology. We conclude that the SL repeat is not an informative phylogenetic marker for long range evolutionary studies of genus Trypanosoma.

Spliced leader dependent trans-splicing (SLTS) has been described as an important RNA regulatory process that occurs in different organisms, including the trematode Schistosoma mansoni. We identified more than seven thousand putative SLTS sites in the parasite, comprising genes with a wide spectrum of functional classes, which underlines the SLTS as a ubiquitous mechanism in the parasite. Also, SLTS gene expression levels span several orders of magnitude, showing that SLTS frequency is not determined by the expression level of the target gene, but by the presence of particular gene features facilitating or hindering the trans-splicing mechanism. Our in-depth investigation of SLTS events demonstrates widespread alternative trans-splicing (ATS) acceptor sites occurring in different regions along the entire gene body, highlighting another important role of SLTS generating alternative RNA isoforms in the parasite, besides the polycistron resolution. Particularly for introns where SLTS directly competes for the same acceptor substrate with cis-splicing, we identified for the first time additional and important features that might determine the type of splicing. Our study substantially extends the current knowledge of RNA processing by SLTS in S. mansoni, and provide basis for future studies on the trans-splicing mechanism in other eukaryotes.

Spliced leader trans-splicing (SLTS) is a poorly understood mechanism that is found in a diversity of eukaryotic lineages. In SLTS, a short RNA sequence is added near the 5' ends of the transcripts of protein-coding genes by a modified spliceosomal reaction. Available data suggest that SLTS has evolved many times, and might be more likely to evolve in animals. That SLTS might be more likely to evolve in the context of the generally complex transcriptomes characteristic of animals suggests the possibility that SLTS functions in gene regulation or transcriptome diversification, however no general novel function for SLTS is known. Here, I report SLTS in a lineage of cellularly complex unicellular eukaryotes. Cryptomonads are a group of eukaryotic algae that acquired photosynthetic capacity by secondary endosymbiosis of a red alga, and that retain a reduced copy of the nucleus of the engulfed alga. I estimate that at least one-fifth of genes in the model cryptomonad Guillardia theta and its relative Hanusia phi undergo SLTS. I show that hundreds of genes in G. theta generate alternative transcripts by SLTS at alternative sites, however I find little evidence for alternative protein production by alternative SLTS splicing. Interestingly, I find no evidence for substantial operon structure in the G. theta genome, in contrast to previous findings in other lineages with SLTS. These results extend SLTS to another major group of eukaryotes, and heighten the mystery of the evolution of SLTS and its association with cellular and transcriptomic complexity.

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.

Since the initial publication of the trypanosomatid genomes, curation has been ongoing. Here we make use of existing Trypanosoma brucei ribosome profiling data to provide evidence of ribosome occupancy (and likely translation) of mRNAs from 225 currently unannotated coding sequences (CDSs). A small number of these putative genes correspond to extra copies of previously annotated genes, but 85% are novel. The median size of these novels CDSs is small (81 aa), indicating that past annotation work has excelled at detecting large CDSs. Of the unique CDSs confirmed here, over half have candidate orthologues in other trypanosomatid genomes, most of which were not yet annotated as protein-coding genes. Nonetheless, approximately one-third of the new CDSs were found only in T. brucei subspecies. Using ribosome footprints, RNA-Seq and spliced leader mapping data, we updated previous work to definitively revise the start sites for 414 CDSs as compared to the current gene models. The data pointed to several regions of the genome that had sequence errors that altered coding region boundaries. Finally, we consolidated this data with our previous work to propose elimination of 683 putative genes as protein-coding and arrive at a view of the translatome of slender bloodstream and procyclic culture form T. brucei.

Spliced-leader trans-splicing (SLTS) has been described in distantly related eukaryotes and acts to mark mRNAs with a short 5' exon, giving different mRNAs identical 5' sequence-signatures. The function of these systems is obscure. Perkinsozoa encompasses a diversity of parasitic protists that infect bivalves, toxic-tide dinoflagellates, fish and frog tadpoles. Here, we report considerable sequence variation in the SLTS-system across the Perkinsozoa and find that multiple variant SLTS-systems are encoded in parallel in the ecologically important Perkinsozoa parasite Parvilucifera sinerae. These results demonstrate that the transcriptome of P. sinerae is segregated based on the addition of different spliced-leader (SL) exons. This segregation marks different gene categories, suggesting that SL-segregation relates to functional differentiation of the transcriptome. By contrast, both sets of gene categories are present in the single SL-transcript type sampled from Maranthos, implying that the SL-segregation of the Parvilucifera transcriptome is a recent evolutionary innovation. Furthermore, we show that the SLTS-system marks a subsection of the transcriptome with increased mRNA abundance and includes genes that encode the spliceosome system necessary for SLTS-function. Collectively, these data provide a picture of how the SLTS-systems can vary within a major evolutionary group and identify how additional transcriptional-complexity can be achieved through SL-segregation.

We have previously identified a major proximal sequence element (PSE) responsible for transcription of the spliced leader (SL) gene from Trypanosoma cruzi strain CL, and showed that the sequence encompassing this PSE exhibits approximately 30% divergence between two major groups of T. cruzi isolates, but strong conservation within the groups. In this report, we show that the SL RNA gene promoter from the CL strain (group I) is efficiently expressed only in T. cruzi isolates from group I. Similarly, the sequence of the approximately 643 bp promoter region of the T. cruzi rRNA is strongly conserved within, but diverged approximately 20% between, the two groups. Reporter constructs driven by the rRNA promoter sequences from group I strains are strongly expressed after electroporation into other group I strains, but not expressed in group II strains. In contrast, constructs bearing rRNA promoter sequences from group II strains are active in strains from both groups. Phylogenetic analyses performed with both the rRNA and the SL RNA gene promoter sequences yielded similar trees, and these trees strongly reinforce the partitioning of known T. cruzi into two major groups that parallel the observed functional specificity of the promoters. Given the well-documented species specific pattern of both rRNA promoters and PSEs in higher eukaryotes, these results suggest an ancient evolutionary divergence among organisms currently classified as T. cruzi.

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.

THE SPLICED LEADER (SL) IS A GENE THAT GENERATES A FUNCTIONAL NCRNA THAT IS COMPOSED OF TWO REGIONS: an intronic region of unknown function (SLi) and an exonic region (SLe), which is transferred to the 5' end of independent transcripts yielding mature mRNAs, in a process known as spliced leader trans-splicing (SLTS). The best described function for SLTS is to solve polycistronic transcripts into monocistronic units, specifically in Trypanosomatids. In other metazoans, it is speculated that the SLe addition could lead to increased mRNA stability, differential recruitment of the translational machinery, modification of the 5' region or a combination of these effects. Although important aspects of this mechanism have been revealed, several features remain to be elucidated. We have analyzed 157 SLe sequences from 148 species from seven phyla and found a high degree of conservation among the sequences of species from the same phylum, although no considerable similarity seems to exist between sequences of species from different phyla. When analyzing case studies, we found evidence that a given SLe will always be related to a given set of transcripts in different species from the same phylum, and therefore, different SLe sequences from the same species would regulate different sets of transcripts. In addition, we have observed distinct transcript categories to be preferential targets for the SLe addition in different phyla. This work sheds light into crucial and controversial aspects of the SLTS mechanism. It represents a comprehensive study concerning various species and different characteristics of this important post-transcriptional regulatory mechanism.

Most microbial eukaryotes are uncultivated and thus poorly suited to standard genomic techniques. This is the case for Polykrikos lebouriae, a dinoflagellate with ultrastructurally aberrant plastids. It has been suggested that these plastids stem from a novel symbiosis with either a diatom or haptophyte, but this hypothesis has been difficult to test as P. lebouriae dwells in marine sand rife with potential genetic contaminants.

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.

Infections with parasitic nematodes are among the most significant of the neglected tropical diseases affecting about a billion people living mainly in tropical regions with low economic activity. The most effective current strategy to control nematode infections involves large scale treatment programs with anthelmintic drugs. This strategy is at risk from the emergence of drug resistant parasites. Parasitic nematodes also affect livestock, which are treated with the same limited group of anthelmintic drugs. Livestock parasites resistant to single drugs, and even multi-drug resistant parasites, are appearing in many areas. There is therefore a pressing need for new anthelmintic drugs. Here we use the nematode Caenorhabditis elegans as a model for parasitic nematodes and demonstrate that sinefungin, a competitive inhibitor of methyltransferases, causes a delay in development and reduced fecundity, and inhibits spliced leader trans-splicing. Spliced leader trans-splicing is an essential step in gene expression that does not occur in the hosts of parasitic nematodes, and is therefore a potential target for new anthelmintic drugs. We have exploited the ability of sinefungin to inhibit spliced leader trans-splicing to adapt a green fluorescent protein based reporter gene assay that monitors spliced leader trans-splicing for high-throughput screening for new anthelmintic compounds. We have established a protocol for robust high-throughput screening, combining mechanical dispensing of living C. elegans into 384- or 1536- well plates with addition of compounds using an acoustic liquid dispenser, and the detection of the inhibition of SL trans-splicing using a microplate reader. We have tested this protocol in a first pilot screen and envisage that this assay will be a valuable tool in the search for new anthelmintic drugs.

In the parasite Trypanosoma brucei, the causative agent of human African sleeping sickness, all mRNAs are trans-spliced to generate a common 5' exon derived from the spliced leader (SL) RNA. Perturbations of protein translocation across the endoplasmic reticulum (ER) induce the spliced leader RNA silencing (SLS) pathway. SLS activation is mediated by a serine-threonine kinase, PK3, which translocates from the cytosolic face of the ER to the nucleus, where it phosphorylates the TATA-binding protein TRF4, leading to the shutoff of SL RNA transcription, followed by induction of programmed cell death. Here, we demonstrate that SLS is also induced by depletion of the essential ER-resident chaperones BiP and calreticulin, ER oxidoreductin 1 (ERO1), and the Golgi complex-localized quiescin sulfhydryl oxidase (QSOX). Most strikingly, silencing of Rhomboid-like 1 (TIMRHOM1), involved in mitochondrial protein import, also induces SLS. The PK3 kinase, which integrates SLS signals, is modified by phosphorylation on multiple sites. To determine which of the phosphorylation events activate PK3, several individual mutations or their combination were generated. These mutations failed to completely eliminate the phosphorylation or translocation of the kinase to the nucleus. The structures of PK3 kinase and its ATP binding domain were therefore modeled. A conserved phenylalanine at position 771 was proposed to interact with ATP, and the PK3F771L mutation completely eliminated phosphorylation under SLS, suggesting that the activation involves most if not all of the phosphorylation sites. The study suggests that the SLS occurs broadly in response to failures in protein sorting, folding, or modification across multiple compartments. IMPORTANCE In this study, we found that SLS is induced by depletion of the essential ER-resident chaperones BiP and calreticulin, ER oxidoreductin 1 (ERO1), and the Golgi complex-localized quiescin sulfhydryl oxidase (QSOX). Most strikingly, silencing of Rhomboid-like 1 (TIMRHOM1), involved in mitochondrial protein import, also induces SLS. We also report on the autophosphorylation of PK3 during SLS induction. This study has implications for our understanding of how trypanosomes keep the homeostasis between the ER and the mitochondria and suggests that PK3 may participate in the connection between these two organelles. The pathway, when induced, leads to the suicide of these parasites, and its induction offers a potential novel drug target against these parasites.