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.

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.

Spliced-Leader RNA as a Dynamic Marker for Monitoring Viable Leishmania Parasites During and After Treatment.

Accurate detection of viable Leishmania parasites is critical for evaluating visceral leishmaniasis (VL) treatment response at an early timepoint. We compared the decay of kinetoplast DNA (kDNA) and spliced-leader RNA (SL-RNA) in vitro, in vivo, and in a VL patient cohort. An optimized combination of blood preservation and nucleic acid extraction improved efficiency for both targets. SL-RNA degraded more rapidly during treatment than kDNA, and correlated better with microscopic examination. SL-RNA quantitative polymerase chain reaction emerges as a superior method for dynamic monitoring of viable Leishmania parasites. It enables individualized treatment monitoring for improved prognoses and has potential as an early surrogate endpoint in clinical trials.

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.

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.

Hypothesis: Trans-splicing Generates Evolutionary Novelty in the Photosynthetic Amoeba Paulinella.

Plastid primary endosymbiosis has occurred twice, once in the Archaeplastida ancestor and once in the Paulinella (Rhizaria) lineage. Both events precipitated massive evolutionary changes, including the recruitment and activation of genes that are horizontally acquired (HGT) and the redeployment of existing genes and pathways in novel contexts. Here we address the latter aspect in Paulinella micropora KR01 (hereafter, KR01) that has independently evolved spliced leader (SL) trans-splicing (SLTS) of nuclear-derived transcripts. We investigated the role of this process in gene regulation, novel gene origination, and endosymbiont integration. Our analysis shows that 20% of KR01 genes give rise to transcripts with at least one (but in some cases, multiple) sites of SL addition. This process, which often occurs at canonical cis-splicing acceptor sites (internal introns), results in shorter transcripts that may produce 5'-truncated proteins with novel functions. SL-truncated transcripts fall into four categories that may show: (i) altered protein localization, (ii) altered protein function, structure, or regulation, (iii) loss of valid alternative start codons, preventing translation, or (iv) multiple SL addition sites at the 5'-terminus. The SL RNA genes required for SLTS are putatively absent in the heterotrophic sister lineage of photosynthetic Paulinella species. Moreover, a high proportion of transcripts derived from genes of endosymbiotic gene transfer (EGT) and HGT origin contain SL sequences. We hypothesize that truncation of transcripts by SL addition may facilitate the generation and expression of novel gene variants and that SLTS may have enhanced the activation and fixation of foreign genes in the host genome of the photosynthetic lineages, playing a key role in primary endosymbiont integration.

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.

An RNA polymerase II-associated TFIIF-like complex is indispensable for SL RNA gene transcription in Trypanosoma brucei.

Trypanosomes are protistan parasites that diverged early in evolution from most eukaryotes. Their streamlined genomes are packed with arrays of tandemly linked genes that are transcribed polycistronically by RNA polymerase (pol) II. Individual mRNAs are processed from pre-mRNA by spliced leader (SL) trans splicing and polyadenylation. While there is no strong evidence that general transcription factors are needed for transcription initiation at these gene arrays, a RNA pol II transcription pre-initiation complex (PIC) is formed on promoters of SLRNA genes, which encode the small nuclear SL RNA, the SL donor in trans splicing. The factors that form the PIC are extremely divergent orthologues of the small nuclear RNA-activating complex, TBP, TFIIA, TFIIB, TFIIH, TFIIE and Mediator. Here, we functionally characterized a heterodimeric complex of unannotated, nuclear proteins that interacts with RNA pol II and is essential for PIC formation, SL RNA synthesis in vivo, SLRNA transcription in vitro, and parasite viability. These functional attributes suggest that the factor represents TFIIF although the amino acid sequences are too divergent to firmly make this conclusion. This work strongly indicates that early-diverged trypanosomes have orthologues of each and every general transcription factor, requiring them for the synthesis of SL RNA.

An in vivo genetic screen for genes involved in spliced leader trans-splicing indicates a crucial role for continuous de novo spliced leader RNP assembly.

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.

Cyclin-Dependent Kinase CRK9, Required for Spliced Leader trans Splicing of Pre-mRNA in Trypanosomes, Functions in a Complex with a New L-Type Cyclin and a Kinetoplastid-Specific Protein.

In eukaryotes, cyclin-dependent kinases (CDKs) control the cell cycle and critical steps in gene expression. The lethal parasite Trypanosoma brucei, member of the phylogenetic order Kinetoplastida, possesses eleven CDKs which, due to high sequence divergence, were generically termed CDC2-related kinases (CRKs). While several CRKs have been implied in the cell cycle, CRK9 was the first trypanosome CDK shown to control the unusual mode of gene expression found in kinetoplastids. In these organisms, protein-coding genes are arranged in tandem arrays which are transcribed polycistronically. Individual mRNAs are processed from precursor RNA by spliced leader (SL) trans splicing and polyadenylation. CRK9 ablation was lethal in cultured trypanosomes, causing a block of trans splicing before the first transesterification step. Additionally, CRK9 silencing led to dephosphorylation of RNA polymerase II and to hypomethylation of the SL cap structure. Here, we tandem affinity-purified CRK9 and, among potential CRK9 substrates and modifying enzymes, discovered an unusual tripartite complex comprising CRK9, a new L-type cyclin (CYC12) and a protein, termed CRK9-associated protein (CRK9AP), that is only conserved among kinetoplastids. Silencing of either CYC12 or CRK9AP reproduced the effects of depleting CRK9, identifying these proteins as functional partners of CRK9 in vivo. While mammalian cyclin L binds to CDK11, the CRK9 complex deviates substantially from that of CDK11, requiring CRK9AP for efficient CRK9 complex formation and autophosphorylation in vitro. Interference with this unusual CDK rescued mice from lethal trypanosome infections, validating CRK9 as a potential chemotherapeutic target.

Minimum free energy predicted base pairing in the 39 nt spliced leader and 5' UTR of calmodulin mRNA from Trypanosoma cruzi: influence of the multiple trans-splicing sites.

We analyzed the compositional changes and the stable base pairs in the predicted secondary structure of the 5' UTR calmodulin mRNA in T. cruzi. The three copies of calmodulin in T. cruzi genome display variable position of the trans splicing sites and give rise to several mRNA that differs slightly on 5' UTR composition in the epimastigote stage. We show that the pattern of high probability base pairs in the minimum free energy predicted secondary structures of the calmodulin 5' UTR remains unchanged despite the nucleotide composition variation. However, the 39 nt spliced leader (mini-exon, the 5' exon sequence transferred to trypanosome mRNAs by the mechanism of trans splicing) shows a variable pattern of high and low probability base pairing as consequence of the altered composition of the 5' UTR.

Polycistronic pre-mRNA processing in vitro: snRNP and pre-mRNA role reversal in trans-splicing.

Spliced leader (SL) trans-splicing in Caenorhabditis elegans attaches a 22-nucleotide (nt) exon onto the 5' end of many mRNAs. A particular class of SL, SL2, splices mRNAs of downstream operon genes. Here we use an embryonic extract-based in vitro splicing system to show that SL2 specificity information is encoded within the polycistronic pre-mRNA, and that trans-splicing specificity is recapitulated in vitro. We define an RNA sequence required for SL2 trans-splicing, the U-rich (Ur) element, through mutational analysis and bioinformatics as a short stem-loop followed by a sequence motif, UAYYUU, located approximately 50 nt upstream of the trans-splice site. Furthermore, this element is predicted in intercistronic regions of numerous operons of C. elegans and other species that use SL2 trans-splicing. We propose that the UAYYUU motif hybridizes with the 5' splice site on the SL2 RNA to recruit the SL to the pre-mRNA. In this way, the UAYYUU motif in the pre-mRNA would serve an analogous function to the similar sequence in the U1 snRNA, which binds to the 5' splice site of introns, effectively reversing the roles of snRNP and pre-mRNA in trans-splicing.

On the Possibility of an Early Evolutionary Origin for the Spliced Leader Trans-Splicing.

Trans-splicing is a process by which 5'- and 3'-ends of two pre-RNA molecules transcribed from different sites of the genome can be joined together to form a single RNA molecule. The spliced leader (SL) trans-splicing is mediated by the spliceosome and it allows the replacement of 5'-end of pre-mRNA by 5'(SL)-end of SL-RNA. This form of splicing has been observed in many phylogenetically unrelated eukaryotes. Either the SL trans-splicing (SLTS) originated in the last eukaryotic common ancestor (LECA) (or even earlier) and it was lost in most eukaryotic lineages, or this mechanism of RNA processing evolved several times independently in various unrelated eukaryotic taxa. The bioinformatic comparisons of SL-RNAs from various eukaryotic taxonomic groups have revealed the similarities of secondary structures of most SL-RNAs and a relative conservation of their splice sites (SSs) and Sm-binding sites (SmBSs). We propose that such structural and functional similarities of SL-RNAs are unlikely to have evolved repeatedly many times. Hence, we favor the scenario of an early evolutionary origin for the SLTS and multiple losses of SL-RNAs in various eukaryotic lineages.

Symbiodinium transcriptome and global responses of cells to immediate changes in light intensity when grown under autotrophic or mixotrophic conditions.

Symbiosis between unicellular dinoflagellates (genus Symbiodinium) and their cnidarian hosts (e.g. corals, sea anemones) is the foundation of coral reef ecosystems. Dysfunction of this symbiosis under changing environmental conditions has led to global reef decline. Little information is known about Symbiodinium gene expression and mechanisms by which light impacts host-symbiont associations. To address these issues, we generated a transcriptome from axenic Symbiodinium strain SSB01. Here we report features of the transcriptome, including occurrence and length distribution of spliced leader sequences, the functional landscape of encoded proteins and the impact of light on gene expression. Expression of many Symbiodinium genes appears to be significantly impacted by light. Transcript encoding cryptochrome 2 declined in high light while some transcripts for Regulators of Chromatin Condensation (RCC1) declined in the dark. We also identified a transcript encoding a light harvesting AcpPC protein with homology to Chlamydomonas LHCSR2. The level of this transcript increased in high light autotrophic conditions, suggesting that it is involved in photo-protection and the dissipation of excess absorbed light energy. The most extensive changes in transcript abundances occurred when the algae were transferred from low light to darkness. Interestingly, transcripts encoding several cell adhesion proteins rapidly declined following movement of cultures to the dark, which correlated with a dramatic change in cell surface morphology, likely reflecting the complexity of the extracellular matrix. Thus, light-sensitive cell adhesion proteins may play a role in establishing surface architecture, which may in turn alter interactions between the endosymbiont and its host.

Trypanosoma brucei gambiense Spliced Leader RNA Is a More Specific Marker for Cure of Human African Trypanosomiasis Than T. b. gambiense DNA.

To assess the efficacy of treatment for human African trypanosomiasis, accurate tests that can discriminate relapse from cure are needed. We report the first data that the spliced leader (SL) RNA is a more specific marker for cure of human African trypanosomiasis than parasite DNA. In blood samples obtained from 61 patients in whom human African trypanosomiasis was cured, SL RNA detection had specificities of 98.4%-100%, while DNA detection had a specificity of only 77%. Data from our proof-of-concept study show that SL RNA detection has high potential as a test of cure.

Single-cell transcriptomics using spliced leader PCR: Evidence for multiple losses of photosynthesis in polykrikoid dinoflagellates.

BACKGROUND: 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. RESULTS: We applied spliced-leader targeted PCR (SLPCR) to obtain dinoflagellate-specific transcriptomes on single-cell isolates of P. lebouriae from marine sediments. Polykrikos lebouriae expressed nuclear-encoded photosynthetic genes that were characteristic of the peridinin-plastids of dinoflagellates, rather than those from a diatom of haptophyte. We confirmed these findings at the genomic level using multiple displacement amplification (MDA) to obtain a partial plastome of P. lebouriae. CONCLUSION: From these data, we infer that P. lebouriae has retained the peridinin plastids ancestral for dinoflagellates as a whole, while its closest relatives have lost photosynthesis multiple times independently. We discuss these losses with reference to mixotrophy in polykrikoid dinoflagellates. Our findings demonstrate new levels of variation associated with the peridinin plastids of dinoflagellates and the usefulness of SLPCR approaches on single cell isolates. Unlike other transcriptomic methods, SLPCR has taxonomic specificity, and can in principle be adapted to different splice-leader bearing groups.

RNA-seq analysis of the C. briggsae transcriptome.

Curation of a high-quality gene set is the critical first step in genome research, enabling subsequent analyses such as ortholog assignment, cis-regulatory element finding, and synteny detection. In this project, we have reannotated the genome of Caenorhabditis briggsae, the best studied sister species of the model organism Caenorhabditis elegans. First, we applied a homology-based gene predictor genBlastG to annotate the C. briggsae genome. We then validated and further improved the C. briggsae gene annotation through RNA-seq analysis of the C. briggsae transcriptome, which resulted in the first validated C. briggsae gene set (23,159 genes), among which 7347 genes (33.9% of all genes with introns) have all of their introns confirmed. Most genes (14,812, or 68.3%) have at least one intron validated, compared with only 3.9% in the most recent WormBase release (WS228). Of all introns in the revised gene set (103,083), 61,503 (60.1%) have been confirmed. Additionally, we have identified numerous trans-splicing leaders (SL1 and SL2 variants) in C. briggsae, leading to the first genome-wide annotation of operons in C. briggsae (1105 operons). The majority of the annotated operons (564, or 51.0%) are perfectly conserved in C. elegans, with an additional 345 operons (or 31.2%) somewhat divergent. Additionally, RNA-seq analysis revealed over 10 thousand small-size assembly errors in the current C. briggsae reference genome that can be readily corrected. The revised C. briggsae genome annotation represents a solid platform for comparative genomics analysis and evolutionary studies of Caenorhabditis species.

Direct evidence for RNA-RNA interactions at the 3' end of the Hepatitis C virus genome using surface plasmon resonance.

Surface plasmon resonance was used to investigate two previously described interactions analyzed by reverse genetics and complementation mutation experiments, involving 5BSL3.2, a stem-loop located in the NS5B coding region of HCV. 5BSL3.2 was immobilized on a sensor chip by streptavidin-biotin coupling, and its interaction either with the SL2 stem-loop of the 3' end or with an upstream sequence centered on nucleotide 9110 (referred to as Seq9110) was monitored in real-time. In contrast with previous results obtained by NMR assays with the same short RNA sequences that we used or SHAPE analysis with longer RNAs, we demonstrate that recognition between 5BSL3.2 and SL2 can occur in solution through a kissing-loop interaction. We show that recognition between Seq9110 and the internal loop of 5BSL3.2 does not prevent binding of SL2 on the apical loop of 5BSL3.2 and does not influence the rate constants of the SL2-5BSL3.2 complex. Therefore, the two binding sites of 5BSL3.2, the apical and internal loops, are structurally independent and both interactions can coexist. We finally show that the stem-loop SL2 is a highly dynamic RNA motif that fluctuates between at least two conformations: One is able to hybridize with 5BSL3.2 through loop-loop interaction, and the other one is capable of self-associating in the absence of protein, reinforcing the hypothesis of SL2 being a dimerization sequence. This result suggests also that the conformational dynamics of SL2 could play a crucial role for controlling the destiny of the genomic RNA.

Affinity purification of T7 RNA transcripts with homogeneous ends using ARiBo and CRISPR tags.

Affinity purification of RNA using the ARiBo tag technology currently provides an ideal approach to quickly prepare RNA with 3' homogeneity. Here, we explored strategies to also ensure 5' homogeneity of affinity-purified RNAs. First, we systematically investigated the effect of starting nucleotides on the 5' heterogeneity of a small SLI RNA substrate from the Neurospora VS ribozyme purified from an SLI-ARiBo precursor. A series of 32 SLI RNA sequences with variations in the +1 to +3 region was produced from two T7 promoters (class III consensus and class II 2.5) using either the wild-type T7 RNA polymerase or the P266L mutant. Although the P266L mutant helps decrease the levels of 5'-sequence heterogeneity in several cases, significant levels of 5' heterogeneity (≥1.5%) remain for transcripts starting with GGG, GAG, GCG, GGC, AGG, AGA, AAA, ACA, AUA, AAC, ACC, AUC, and AAU. To provide a more general approach to purifying RNA with 5' homogeneity, we tested the suitability of using a small CRISPR RNA stem-loop at the 5' end of the SLI-ARiBo RNA. Interestingly, we found that complete cleavage of the 5'-CRISPR tag with the Cse3 endoribonuclease can be achieved quickly from CRISPR-SLI-ARiBo transcripts. With this procedure, it is possible to generate SLI-ARiBo RNAs starting with any of the four standard nucleotides (G, C, A, or U) involved in either a single- or a double-stranded structure. Moreover, the 5'-CRISPR-based strategy can be combined with affinity purification using the 3'-ARiBo tag for quick purification of RNA with both 5' and 3' homogeneity.

RNA folding pathways and kinetics using 2D energy landscapes.

RNA folding pathways play an important role in various biological processes, such as (i) the hok/sok (host-killing/suppression of killing) system in E. coli to check for sufficient plasmid copy number, (ii) the conformational switch in spliced leader (SL) RNA from Leptomonas collosoma, which controls trans splicing of a portion of the '5 exon, and (iii) riboswitches--portions of the 5' untranslated region of messenger RNA that regulate genes by allostery. Since RNA folding pathways are determined by the energy landscape, we describe a novel algorithm, FFTbor2D, which computes the 2D projection of the energy landscape for a given RNA sequence. Given two metastable secondary structures A, B for a given RNA sequence, FFTbor2D computes the Boltzmann probability p(x, y) = Z(x,y)/Z that a secondary structure has base pair distance x from A and distance y from B. Using polynomial interpolationwith the fast Fourier transform,we compute p(x, y) in O(n(5)) time and O(n(2)) space, which is an improvement over an earlier method, which runs in O(n(7)) time and O(n(4)) space. FFTbor2D has potential applications in synthetic biology, where one might wish to design bistable switches having target metastable structures A, B with favorable pathway kinetics. By inverting the transition probability matrix determined from FFTbor2D output, we show that L. collosoma spliced leader RNA has larger mean first passage time from A to B on the 2D energy landscape, than 97.145% of 20,000 sequences, each having metastable structures A, B. Source code and binaries are freely available for download at http://bioinformatics.bc.edu/clotelab/FFTbor2D. The program FFTbor2D is implemented in C++, with optional OpenMP parallelization primitives.