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.     

Spliced leader trans-splicing in the nematode Trichinella spiralis uses highly polymorphic, noncanonical spliced leaders

The trans-splicing of short spliced leader (SL) RNAs onto the 5' ends of mRNAs occurs in a diverse range of taxa. In nematodes, all species so far characterized utilize a characteristic, conserved spliced leader, SL1, as well as variants that are employed in the resolution of operons. Here we report the identification of spliced leader trans-splicing in the basal nematode Trichinella spiralis, and show that this nematode does not possess a canonical SL1, but rather has at least 15 distinct spliced leaders, encoded by at least 19 SL RNA genes. The individual spliced leaders vary in both size and primary sequence, showing a much higher degree of diversity compared to other known trans-spliced leaders. In a survey of T. spiralis mRNAs, individual mRNAs were found to be trans-spliced to a number of different spliced leader sequences. These data provide the first indication that the last common ancestor of the phylum Nematoda utilized spliced leader trans-splicing and that the canonical spliced leader, SL1, found in Caenorhabditis elegans, evolved after the divergence of the major nematode clades. This discovery sheds important light on the nature and evolution of mRNA processing in the Nematoda.     

Increased efficiency of evolved group I intron spliceozymes by decreased side product formation

The group I intron ribozyme from Tetrahymena was recently reengineered into a trans-splicing variant that is able to remove 100-nt introns from pre-mRNA, analogous to the spliceosome. These spliceozymes were improved in this study by 10 rounds of evolution in Escherichia coli cells. One clone with increased activity in E. coli cells was analyzed in detail. Three of its 10 necessary mutations extended the substrate binding duplexes, which led to increased product formation and reduced cleavage at the 5′-splice site. One mutation in the conserved core of the spliceozyme led to a further reduction of cleavage at the 5′-splice site but an increase in cleavage side products at the 3′-splice site. The latter was partially reduced by six additional mutations. Together, the mutations increased product formation while reducing activity at the 5′-splice site and increasing activity at the 3′-splice site. These results show the adaptation of a ribozyme that evolved in nature for cis-splicing to trans-splicing, and they highlight the interdependent function of nucleotides within group I intron ribozymes. Implications for the possible use of spliceozymes as tools in research and therapy, and as a model for the evolution of the spliceosome, are discussed.     

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.     

A conformational switch in the DiGIR1 ribozyme involved in release and folding of the downstream I-DirI mRNA

DiGIR1 is a group I-like cleavage ribozyme found as a structural domain within a nuclear twin-ribozyme group I intron. DiGIR1 catalyzes cleavage by branching at an Internal Processing Site (IPS) leading to formation of a lariat cap at the 5′-end of the 3′-cleavage product. The 3′-cleavage product is subsequently processed into an mRNA encoding a homing endonuclease. By analysis of combinations of 5′- and 3′-deletions, we identify a hairpin in the 5′-UTR of the mRNA (HEG P1) that is formed by conformational switching following cleavage. The formation of HEG P1 inhibits the reversal of the branching reaction, thus giving it directionality. Furthermore, the release of the mRNA is a consequence of branching rather than hydrolytic cleavage. A model is put forward that explains the release of the I-DirI mRNA with a lariat cap and a structured 5′-UTR as a direct consequence of the DiGIR1 branching reaction. The role of HEG P1 in GIR1 branching is reminiscent of that of hairpin P-1 in splicing of the Tetrahymena rRNA group I intron and illustrates a general principle in RNA-directed RNA processing.     

Protein trans-splicing and characterization of a split family B-type DNA polymerase from the hyperthermophilic archaeal parasite Nanoarchaeum equitans

Nanoarchaeum equitans family B-type DNA polymerase (Neq DNA polymerase) is encoded by two separate genes, the large gene coding for the N-terminal part (Neq L) of Neq DNA polymerase and the small gene coding for the C-terminal part (Neq S), including a split mini-intein sequence. The two Neq DNA polymerase genes were cloned and expressed in Escherichia coli individually, together (for the Neq C), and as a genetically protein splicing-processed form (Neq P). The protein trans-spliced Neq C was obtained using the heating step at 80 degrees C after the co-expression of the two genes. The protein trans-splicing of the N-terminal and C-terminal parts of Neq DNA polymerase was examined in vitro using the purified Neq L and Neq S. The trans-splicing was influenced mainly by temperature, and occurred only at temperatures above 50 degrees C. The trans-splicing reaction was inhibited in the presence of zinc. Neq S has no catalytic activity and Neq L has lower 3'-->5' exonuclease activity; whereas Neq C and Neq P have polymerase and 3'-->5' exonuclease activities, indicating that both Neq L and Neq S are needed to form the active DNA polymerase that possesses higher proofreading activity. The genetically protein splicing-processed Neq P showed the same properties as the protein trans-spliced Neq C. Our results are the first evidence to show experimentally that natural protein trans-splicing occurs in an archaeal protein, a thermostable protein, and a family B-type DNA polymerase.     

Circular RNAs are abundant,conserved, and associated with ALU repeats

Circular RNAs composed of exonic sequence have been described in a small number of genes. Thought to result from splicing errors, circular RNA species possess no known function. To delineate the universe of endogenous circular RNAs, we performed high-throughput sequencing (RNA-seq) of libraries prepared from ribosome-depleted RNA with or without digestion with the RNA exonuclease, RNase R. We identified >25,000 distinct RNA species in human fibroblasts that contained non-colinear exons (a “backsplice”) and were reproducibly enriched by exonuclease degradation of linear RNA. These RNAs were validated as circular RNA (ecircRNA), rather than linear RNA, and were more stable than associated linear mRNAs in vivo. In some cases, the abundance of circular molecules exceeded that of associated linear mRNA by >10-fold. By conservative estimate, we identified ecircRNAs from 14.4% of actively transcribed genes in human fibroblasts. Application of this method to murine testis RNA identified 69 ecircRNAs in precisely orthologous locations to human circular RNAs. Of note, paralogous kinases HIPK2 and HIPK3 produce abundant ecircRNA from their second exon in both humans and mice. Though HIPK3 circular RNAs contain an AUG translation start, it and other ecircRNAs were not bound to ribosomes. Circular RNAs could be degraded by siRNAs and, therefore, may act as competing endogenous RNAs. Bioinformatic analysis revealed shared features of circularized exons, including long bordering introns that contained complementary ALU repeats. These data show that ecircRNAs are abundant, stable, conserved and nonrandom products of RNA splicing that could be involved in control of gene expression.     

Relationship between the self-splicing activity and the solidity of the master domain of the Tetrahymena group I ribozyme

The highly conserved P3-P7 domain of the Group I intron ribozymes is known to contain essential elements, such as the binding site for the cofactor guanosine, required for conducting the splicing reaction. We investigated the domain of the Tetrahymena intron ribozyme and its variants in order to clarify the relationship between its stability and function. We found that the destabilization of the P3-P7 domain facilitates the active structure formation at high magnesium ion concentrations where the formation is retarded for the wild type. The destabilized domain also increases K(GTP)(m) although this can be compensated by increasing the concentration of Mg(2+), indicating that the stable domain is required for establishing a tight guanosine binding site. The results suggest that the stability of the domain affects the rate-limiting step in the RNA folding pathway and also regulates the efficiency of the splicing reaction.     

A novel spliceosome-mediated trans-splicing can change our view on genome complexity of the divergent eukaryote Giardia intestinalis

Although spliceosomal introns are an abundant landmark in eukaryotic genomes, the nuclear genome of the divergent eukaryote Giardia intestinalis, the causative agent of giardiasis, has been considered as “intron-poor” with only five canonical (cis-spliced) introns. However, three research groups (including ours) have independently reported a novel class of spliceosomal introns in the G. intestinalis genome. Three protein-coding genes are split into pieces in the G. intestinalis genome, and each of the partial coding regions was independently transcribed into polyadenylated premature mRNAs (pre-mRNAs). The two pre-mRNAs directly interact with each other by an intermolecular-stem structure formed between their non-coding portions, and are then processed into mature mRNAs by spliceosome-mediated trans-splicing. Here, we summarize the recently published works on split introns (“splintrons”) in the G. intestinalis genome, and then provide our speculation on the functional property of the Giardia spliceosomes based on the putative ratio of splintrons to canonical introns. Finally, we discuss a scenario for the transition from typical GT-AG boundaries to non-typical AT-AC boundaries in a particular splintron of Giardia.     

The linear form of a group II intron catalyzes efficient autocatalytic reverse splicing, establishing a potential for mobility

Self-splicing group II introns catalyze their own excision from pre-RNAs, thereby joining the flanking exons. The introns can be released in a lariat or linear form. Lariat introns have been shown to reverse the splicing reaction; in contrast, linear introns are generally believed to perform no or only poor reverse splicing. Here, we show that a linear group II intron derived from ai5γ can reverse the second step of splicing with unexpectedly high efficiency and precision. Moreover, the linear intron generates dramatically more reverse-splicing product than its lariat equivalent. The finding that linear group II introns can readily undergo the critical first step of mobility by catalyzing efficient reverse splicing into complementary target molecules demonstrates their innate potential for mobility and transposition and raises the possibility that reverse splicing by linear group II introns may have played a significant role in certain forms of intron mobility and lateral gene transfer during evolution.     

Metal ion coordination to 2′ functionality of guanosine mediates substrate-guanosine coupling in group I ribozymes: implications for conserved role of metal ions and for variability in RNA folding in ribozyme catalysis

Divalent metal ions are necessary in the self splicing reaction of group I introns, and we report that metal interaction to the 2′ position of guanosine for the Azoarcus ribozyme is required for catalysis. Moreover, this metal coordination promotes the guanosine-substrate coupled binding to the ribozyme, which is another conserved feature seen across phylogenetic boundaries. Typically there is a 4-9-fold difference in binding of G to E_free versus E · S. In the Tetrahymena ribozyme’s case this substrate-guanosine communication was attributed to conformational change(s) that lead to cooperative binding of the two cofactors which is almost nonexistent at low temperatures (4 °C). In the prokaryotic Azoarcus ribozyme we also see a 4-5-fold difference in binding of the guanosine/substrate to E_free versus E · G or E · S at 10 °C that is attributed to guanosine-substrate coupling. This coupling is diminished when the metal (Mg²⁺) coordination to the 2′ is disrupted with use of 2′-amino-2′-deoxyguanosine. The coupling is restored when softer Mn²⁺ ions are added to the buffer. This evidence generalizes a model for group I ribozyme catalysis that involves metal coordination to the 2′ position of guanosine. However, we see one striking difference in that the guanosine-substrate coupling is reversed. In the Azoarcus system (10 °C) the guanosine/substrate binds 5-fold more tightly to E_free than to E · S or E · G, which is the opposite for Tetrahymena even when the later is run at 4 °C. One implication for this difference in coupling is that the Azoarcus is in a folded state well accommodated for guanosine or substrate binding. This initial binding actually causes a conformational change that retards the subsequent binding of the second cofactor, which contrasts what was found for the Tetrahymena ribozyme. These results indicate that while the role for the metal ions in the chemical catalysis is conserved across phylogenetic boundaries, there is variability in the folding pattern of the ribozyme that leads to phosphoryl transfer.