Convergent origins and rapid evolution of spliced leader trans-splicing in Metazoa: Insights from the Ctenophora and Hydrozoa

Replacement of mRNA 5′ UTR sequences by short sequences trans-spliced from specialized, noncoding, spliced leader (SL) RNAs is an enigmatic phenomenon, occurring in a set of distantly related animal groups including urochordates, nematodes, flatworms, and hydra, as well as in Euglenozoa and dinoflagellates. Whether SL trans-splicing has a common evolutionary origin and biological function among different organisms remains unclear. We have undertaken a systematic identification of SL exons in cDNA sequence data sets from non-bilaterian metazoan species and their closest unicellular relatives. SL exons were identified in ctenophores and in hydrozoan cnidarians, but not in other cnidarians, placozoans, or sponges, or in animal unicellular relatives. Mapping of SL absence/presence obtained from this and previous studies onto current phylogenetic trees favors an evolutionary scenario involving multiple origins for SLs during eumetazoan evolution rather than loss from a common ancestor. In both ctenophore and hydrozoan species, multiple SL sequences were identified, showing high sequence diversity. Detailed analysis of a large data set generated for the hydrozoan Clytia hemisphaerica revealed trans-splicing of given mRNAs by multiple alternative SLs. No evidence was found for a common identity of trans-spliced mRNAs between different hydrozoans. One feature found specifically to characterize SL-spliced mRNAs in hydrozoans, however, was a marked adenosine enrichment immediately 3′ of the SL acceptor splice site. Our findings of high sequence divergence and apparently indiscriminate use of SLs in hydrozoans, along with recent findings in other taxa, indicate that SL genes have evolved rapidly in parallel in diverse animal groups, with constraint on SL exon sequence evolution being apparently rare.     

Group II Introns Generate Functional Chimeric Relaxase Enzymes with Modified Specificities through Exon Shuffling at Both the RNA and DNA Level

Group II introns are large self-splicing RNA enzymes with a broad but somewhat irregular phylogenetic distribution. These ancient retromobile elements are the proposed ancestors of approximately half the human genome, including the abundant spliceosomal introns and non-long terminal repeat retrotransposons. In contrast to their eukaryotic derivatives, bacterial group II introns have largely been considered as harmful selfish mobile retroelements that parasitize the genome of their host. As a challenge to this view, we recently uncovered a new intergenic trans-splicing pathway that generates an assortment of mRNA chimeras. The ability of group II introns to combine disparate mRNA fragments was proposed to increase the genetic diversity of the bacterial host by shuffling coding sequences. Here, we show that the Ll.LtrB and Ef.PcfG group II introns from Lactococcus lactis and Enterococcus faecalis respectively can both use the intergenic trans-splicing pathway to catalyze the formation of chimeric relaxase mRNAs and functional proteins. We demonstrated that some of these compound relaxase enzymes yield gain-of-function phenotypes, being significantly more efficient than their precursor wild-type enzymes at supporting bacterial conjugation. We also found that relaxase enzymes with shuffled functional domains are produced in biologically relevant settings under natural expression levels. Finally, we uncovered examples of lactococcal chimeric relaxase genes with junctions exactly at the intron insertion site. Overall, our work demonstrates that the genetic diversity generated by group II introns, at the RNA level by intergenic trans-splicing and at the DNA level by recombination, can yield new functional enzymes with shuffled exons, which can lead to gain-of-function phenotypes.     

Characterization of RNase R-digested cellular RNA source that consists of lariat and circular RNAs from pre-mRNA splicing

Besides linear RNAs, pre-mRNA splicing generates three forms of RNAs: lariat introns, Y-structure introns from trans-splicing, and circular exons through exon skipping. To study the persistence of excised introns in total cellular RNA, we used three Escherichia coli 3′ to 5′ exoribonucleases. Ribonuclease R (RNase R) thoroughly degrades the abundant linear RNAs and the Y-structure RNA, while preserving the loop portion of a lariat RNA. Ribonuclease II (RNase II) and polynucleotide phosphorylase (PNPase) also preserve the lariat loop, but are less efficient in degrading linear RNAs. RNase R digestion of the total RNA from human skeletal muscle generates an RNA pool consisting of lariat and circular RNAs. RT–PCR across the branch sites confirmed lariat RNAs and circular RNAs in the pool generated by constitutive and alternative splicing of the dystrophin pre-mRNA. Our results indicate that RNase R treatment can be used to construct an intronic cDNA library, in which majority of the intron lariats are represented. The highly specific activity of RNase R implies its ability to screen for rare intragenic trans-splicing in any target gene with a large background of cis-splicing. Further analysis of the intronic RNA pool from a specific tissue or cell will provide insights into the global profile of alternative splicing.     

Dark-induced mRNA instability involves RNase E/G-type endoribonuclease cleavage at the AU-box and SD sequences in cyanobacteria

Light-responsive gene expression is crucial to photosynthesizing organisms. Here, we studied functions of cis-elements (AU-box and SD sequences) and a trans-acting factor (ribonuclease, RNase) in light-responsive expression in cyanobacteria. The results indicated that AU-rich nucleotides with an AU-box, UAAAUAAA, just upstream from an SD confer instability on the mRNA under darkness. An RNase E/G homologue, Slr1129, of the cyanobacterium Synechocystis sp. strain PCC 6803 was purified and confirmed capable of endoribonucleolytic cleavage at the AU- (or AG)-rich sequences in vitro. The cleavage depends on the primary target sequence and secondary structure of the mRNA. Complementation tests using Escherichia coli rne/rng mutants showed that Slr1129 fulfilled the functions of both the RNase E and RNase G. An analysis of systematic mutations in the AU-box and SD sequences showed that the cis-elements also affect significantly mRNA stability in light-responsive genes. These results strongly suggested that dark-induced mRNA instability involves RNase E/G-type cleavage at the AU-box and SD sequences in cyanobacteria. The mechanical impact and a possible common mechanism with RNases for light-responsive gene expression are discussed. Electronic supplementary material The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Exonic Sequences in the 5′ Untranslated Region of α-Tubulin mRNA Modulate trans Splicing in Trypanosoma brucei

Previous studies have identified a conserved AG dinucleotide at the 3′ splice site (3′SS) and a polypyrimidine (pPy) tract that are required for trans splicing of polycistronic pre-mRNAs in trypanosomatids. Furthermore, the pPy tract of the Trypanosoma brucei α-tubulin 3′SS region is required to specify accurate 3′-end formation of the upstream β-tubulin gene and trans splicing of the downstream α-tubulin gene. Here, we employed an in vivo cis competition assay to determine whether sequences other than those of the AG dinucleotide and the pPy tract were required for 3′SS identification. Our results indicate that a minimal α-tubulin 3′SS, from the putative branch site region to the AG dinucleotide, is not sufficient for recognition by the trans-splicing machinery and that polyadenylation is strictly dependent on downstream trans splicing. We show that efficient use of the α-tubulin 3′SS is dependent upon the presence of exon sequences. Furthermore, β-tubulin, but not actin exon sequences or unrelated plasmid sequences, can replace α-tubulin exon sequences for accurate trans-splice-site selection. Taken together, these results support a model in which the informational content required for efficient trans splicing of the α-tubulin pre-mRNA includes exon sequences which are involved in modulation of trans-splicing efficiency. Sequences that positively regulate trans splicing might be similar to cis-splicing enhancers described in other systems.     

Operons Are a Conserved Feature of Nematode Genomes

The organization of genes into operons, clusters of genes that are co-transcribed to produce polycistronic pre-mRNAs, is a trait found in a wide range of eukaryotic groups, including multiple animal phyla. Operons are present in the class Chromadorea, one of the two main nematode classes, but their distribution in the other class, the Enoplea, is not known. We have surveyed the genomes of Trichinella spiralis, Trichuris muris, and Romanomermis culicivorax and identified the first putative operons in members of the Enoplea. Consistent with the mechanism of polycistronic RNA resolution in other nematodes, the mRNAs produced by genes downstream of the first gene in the T. spiralis and T. muris operons are trans-spliced to spliced leader RNAs, and we are able to detect polycistronic RNAs derived from these operons. Importantly, a putative intercistronic region from one of these potential enoplean operons confers polycistronic processing activity when expressed as part of a chimeric operon in Caenorhabditis elegans. We find that T. spiralis genes located in operons have an increased likelihood of having operonic C. elegans homologs. However, operon structure in terms of synteny and gene content is not tightly conserved between the two taxa, consistent with models of operon evolution. We have nevertheless identified putative operons conserved between Enoplea and Chromadorea. Our data suggest that operons and “spliced leader” (SL) trans-splicing predate the radiation of the nematode phylum, an inference which is supported by the phylogenetic profile of proteins known to be involved in nematode SL trans-splicing.     

An in vivo selection method to optimize trans-splicing ribozymes

Group I intron ribozymes can repair mutated mRNAs by replacing the 3′-terminal portion of the mRNA with their own 3′-exon. This trans-splicing reaction has the potential to treat genetic disorders and to selectively kill cancer cells or virus-infected cells. However, these ribozymes have not yet been used in therapy, partially due to a low in vivo trans-splicing efficiency. Previous strategies to improve the trans-splicing efficiencies focused on designing and testing individual ribozyme constructs. Here we describe a method that selects the most efficient ribozymes from millions of ribozyme variants. This method uses an in vivo rescue assay where the mRNA of an inactivated antibiotic resistance gene is repaired by trans-splicing group I intron ribozymes. Bacterial cells that express efficient trans-splicing ribozymes are able to grow on medium containing the antibiotic chloramphenicol. We randomized a 5′-terminal sequence of the Tetrahymena thermophila group I intron and screened a library with 9 × 10⁶ ribozyme variants for the best trans-splicing activity. The resulting ribozymes showed increased trans-splicing efficiency and help the design of efficient trans-splicing ribozymes for different sequence contexts. This in vivo selection method can now be used to optimize any sequence in trans-splicing ribozymes.