DiGIR1 and NaGIR1: naturally occurring group I-like ribozymes with unique core organization and evolved biological role

Analysis of small and large subunit rDNA introns from several ectomycorrhizal fungi species

The small (18S) and large (28S) nuclear ribosomal DNA (rDNA) introns have been researched and sequenced in a variety of ectomycorrhizal fungal taxa in this study, it is found that both 18S and 28S rDNA would contain introns and display some degree variation in size, nucleotide sequences and insertion positions within the same fungi species (Meliniomyces). Under investigations among the tested isolates, 18S rDNA has four sites for intron insertions, 28S rDNA has two sites for intron insertions. Both 18S and 28S rDNA introns among the tested isolates belong to group I introns with a set of secondary structure elements designated P1-P10 helics and loops. We found a 12 nt nucleotide sequences TACCACAGGGAT at site 2 in the 3'-end of 28S rDNA, site 2 introns just insert the upstream or the downstream of the12 nt nucleotide sequences. Afters sequence analysis of all 18S and 28S rDNA introns from tested isolates, three high conserved regions around 30 nt nucleotides (conserved 1, conserved 2, conserved 3) and identical nucleotides can be found. Conserved 1, conserved 2 and conserved 3 regions have high GC content, GC percentage is almost more than 60%. From our results, it seems that the more convenient host sites, intron sequences and secondary structures, or isolates for 18S and 28S rDNA intron insertion and deletion, the more popular they are. No matter 18S rDNA introns or 18S rDNA introns among tested isolates, complementary base pairing at the splicing sites in P1-IGS-P10 tertiary helix around 5'-end introns and exons were weak.

Convergent evolution of twintron-like configurations: One is never enough

Introns inserted within introns are commonly referred to as twintrons, however the original definition for twintron implied that splicing of the external member of the twintron could only proceed upon splicing of the internal member. This review examines the various types of twintron-like arrangements that have been reported and assigns them to either nested or twintron categories that are subdivided further into subtypes based on differences of their mode of splicing. Twintron-like arrangements evolved independently by fortuitous events among different types of introns but once formed they offer opportunities for the evolution of new regulatory strategies and/or novel genetic elements.

Speciation of a group I intron into a lariat capping ribozyme

The lariat-capping (LC) ribozyme is a natural ribozyme isolated from eukaryotic microorganisms. Despite apparent structural similarity to group I introns, the LC ribozyme catalyzes cleavage by a 2',5' branching reaction, leaving the 3' product with a 3-nt lariat cap that functionally substitutes for a conventional mRNA cap in the downstream pre-mRNA encoding a homing endonuclease. We describe the crystal structures of the precleavage and postcleavage LC ribozymes, which suggest that structural features inherited from group I ribozymes have undergone speciation due to profound changes in molecular selection pressure, ultimately giving rise to an original branching ribozyme family. The structures elucidate the role of key elements that regulate the activity of the LC ribozyme by conformational switching and suggest a mechanism by which the signal for branching is transmitted to the catalytic core. The structures also show how conserved interactions twist residues, forming the lariat to join chemical groups involved in branching.

Predicting coaxial helical stacking in RNA junctions

RNA junctions are important structural elements that form when three or more helices come together in space in the tertiary structures of RNA molecules. Determining their structural configuration is important for predicting RNA 3D structure. We introduce a computational method to predict, at the secondary structure level, the coaxial helical stacking arrangement in junctions, as well as classify the junction topology. Our approach uses a data mining approach known as random forests, which relies on a set of decision trees trained using length, sequence and other variables specified for any given junction. The resulting protocol predicts coaxial stacking within three- and four-way junctions with an accuracy of 81% and 77%, respectively; the accuracy increases to 83% and 87%, respectively, when knowledge from the junction family type is included. Coaxial stacking predictions for the five to ten-way junctions are less accurate (60%) due to sparse data available for training. Additionally, our application predicts the junction family with an accuracy of 85% for three-way junctions and 74% for four-way junctions. Comparisons with other methods, as well applications to unsolved RNAs, are also presented. The web server Junction-Explorer to predict junction topologies is freely available at: http://bioinformatics.njit.edu/junction.

Intermolecular interaction between a branching ribozyme and associated homing endonuclease mRNA

RNA tertiary interactions involving docking of GNRA (N; any base; R; purine) hairpin loops into helical stem structures on other regions of the same RNA are one of the most common RNA tertiary interactions. In this study, we investigated a tertiary association between a GAAA hairpin tetraloop in a small branching ribozyme (DiGIR1) and a receptor motif (HEG P1 motif) present in a hairpin structure on a separate mRNA molecule. DiGIR1 generates a 2', 5' lariat cap at the 5' end of its downstream homing endonuclease mRNA by catalysing a self-cleavage branching reaction at an internal processing site. Upon release, the 5' end of the mRNA forms a distinct hairpin structure termed HEG P1. Our biochemical data, in concert with molecular 3D modelling, provide experimental support for an intermolecular tetraloop receptor interaction between the L9 GAAA in DiGIR1 and a GNRA tetraloop receptor-like motif (UCUAAG-CAAGA) found within the HEG P1. The biological role of this interaction appears to be linked to the homing endonuclease expression by promoting post-cleavage release of the lariat capped mRNA. These findings add to our understanding of how protein-coding genes embedded in nuclear ribosomal DNA are expressed in eukaryotes and controlled by ribozymes.

Learning to live together: mutualism between self-splicing introns and their hosts

Group I and II introns can be considered as molecular parasites that interrupt protein-coding and structural RNA genes in all domains of life. They function as self-splicing ribozymes and thereby limit the phenotypic costs associated with disruption of a host gene while they act as mobile DNA elements to promote their spread within and between genomes. Once considered purely selfish DNA elements, they now seem, in the light of recent work on the molecular mechanisms regulating bacterial and phage group I and II intron dynamics, to show evidence of co-evolution with their hosts. These previously underappreciated relationships serve the co-evolving entities particularly well in times of environmental stress.

The Naegleria genome: a free-living microbial eukaryote lends unique insights into core eukaryotic cell biology

Naegleria gruberi, a free-living protist, has long been treasured as a model for basal body and flagellar assembly due to its ability to differentiate from crawling amoebae into swimming flagellates. The full genome sequence of Naegleria gruberi has recently been used to estimate gene families ancestral to all eukaryotes and to identify novel aspects of Naegleria biology, including likely facultative anaerobic metabolism, extensive signaling cascades, and evidence for sexuality. Distinctive features of the Naegleria genome and nuclear biology provide unique perspectives for comparative cell biology, including cell division, RNA processing and nucleolar assembly. We highlight here exciting new and novel aspects of Naegleria biology identified through genomic analysis.

RNA reprogramming and repair based on trans-splicing group I ribozymes

While many traditional gene therapy strategies attempt to deliver new copies of wild-type genes back to cells harboring the defective genes, RNA-directed strategies offer a range of novel therapeutic applications. Revision or reprogramming of mRNA is a form of gene therapy that modifies mRNA without directly changing the transcriptional regulation or the genomic gene sequence. Group I ribozymes can be engineered to act in trans by recognizing a separate RNA molecule in a sequence-specific manner, and to covalently link a new RNA sequence to this separate RNA molecule. Group I ribozymes have been shown to repair defective transcripts that cause human genetic or malignant diseases, as well as to replace transcript sequences by foreign RNA resulting in new cellular functions. This review provides an overview of current strategies using trans-splicing group I ribozymes in RNA repair and reprogramming.

Natural and unnatural ribozymes: back to the primordial RNA world

We review natural and in vitro selected ribozymes, for which combined studies could provide us with both insight into the functions performed by ancient RNA molecules in a primitive RNA world and a hypothesis about evolutionary steps that led to the contemporary world.

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.

Molecular modelling of the GIR1 branching ribozyme gives new insight into evolution of structurally related ribozymes

Twin-ribozyme introns contain a branching ribozyme (GIR1) followed by a homing endonuclease (HE) encoding sequence embedded in a peripheral domain of a group I splicing ribozyme (GIR2). GIR1 catalyses the formation of a lariat with 3 nt in the loop, which caps the HE mRNA. GIR1 is structurally related to group I ribozymes raising the question about how two closely related ribozymes can carry out very different reactions. Modelling of GIR1 based on new biochemical and mutational data shows an extended substrate domain containing a GoU pair distinct from the nucleophilic residue that dock onto a catalytic core showing a different topology from that of group I ribozymes. The differences include a core J8/7 region that has been reduced and is complemented by residues from the pre-lariat fold. These findings provide the basis for an evolutionary mechanism that accounts for the change from group I splicing ribozyme to the branching GIR1 architecture. Such an evolutionary mechanism can be applied to other large RNAs such as the ribonuclease P.

Expression of protein-coding genes embedded in ribosomal DNA

Ribosomal DNA (rDNA) is a specialised chromosomal location that is dedicated to high-level transcription of ribosomal RNA genes. Interestingly, rDNAs are frequently interrupted by parasitic elements, some of which carry protein genes. These are non-LTR retrotransposons and group II introns that encode reverse transcriptase-like genes, and group I introns and archaeal introns that encode homing endonuclease genes (HEGs). Although rDNA-embedded protein genes are widespread in nuclei, organelles and bacteria, there is surprisingly little information available on how these genes are expressed. Exceptions include a handful of HEGs from group I introns. Recent studies have revealed unusual and essential roles of group I and group I-like ribozymes in the endogenous expression of HEGs. Here we discuss general aspects of rDNA-embedded protein genes and focus on HEG expression from group I introns in the nucleolus.

Short-term sequence evolution and vertical inheritance of the Naegleria twin-ribozyme group I intron

BACKGROUND

Ribosomal DNA of several species of the free-living Naegleria amoeba harbors an optional group I intron within the nuclear small subunit ribosomal RNA gene. The intron (Nae.S516) has a complex organization of two ribozyme domains (NaGIR1 and NaGIR2) and a homing endonuclease gene (NaHEG). NaGIR2 is responsible for intron excision, exon ligation, and full-length intron RNA circularization, reactions typical for nuclear group I intron ribozymes. NaGIR1, however, is essential for NaHEG expression by generating the 5' end of the homing endonuclease messenger RNA. Interestingly, this unusual class of ribozyme adds a lariat-cap at the mRNA.

RESULTS

To elucidate the evolutionary history of the Nae.S516 twin-ribozyme introns we have analyzed 13 natural variants present in distinct Naegleria isolates. Structural variabilities were noted within both the ribozyme domains and provide strong comparative support to the intron secondary structure. One of the introns, present in N. martinezi NG872, contains hallmarks of a degenerated NaHEG. Phylogenetic analyses performed on separate data sets representing NaGIR1, NaGIR2, NaHEG, and ITS1-5.8S-ITS2 ribosomal DNA are consistent with an overall vertical inheritance pattern of the intron within the Naegleria genus.

CONCLUSION

The Nae.S516 twin-ribozyme intron was gained early in the Naegleria evolution with subsequent vertical inheritance. The intron was lost in the majority of isolates (70%), leaving a widespread but scattered distribution pattern. Why the apparent asexual Naegleria amoebae harbors active intron homing endonucleases, dependent on sexual reproduction for its function, remains a puzzle.

An mRNA is capped by a 2', 5' lariat catalyzed by a group I-like ribozyme

Twin-ribozyme introns are formed by two ribozymes belonging to the group I family and occur in some ribosomal RNA transcripts. The group I-like ribozyme, GIR1, liberates the 5' end of a homing endonuclease messenger RNA in the slime mold Didymium iridis. We demonstrate that this cleavage occurs by a transesterification reaction with the joining of the first and the third nucleotide of the messenger by a 2',5'-phosphodiester linkage. Thus, a group I-like ribozyme catalyzes an RNA branching reaction similar to the first step of splicing in group II introns and spliceosomal introns. The resulting short lariat, by forming a protective 5' cap, might have been useful in a primitive RNA world.

Reverse splicing of a mobile twin-ribozyme group I intron into the natural small subunit rRNA insertion site

A mobile group I intron containing two ribozyme domains and a homing endonuclease gene (twin-ribozyme intron organization) can integrate by reverse splicing into the small subunit rRNA of bacteria and yeast. The integration is sequence-specific and corresponds to the natural insertion site (homing site) of the intron. The reverse splicing is independent of the homing endonuclease gene, but is dependent on the group I splicing ribozyme domain. The observed distribution of group I introns in nature can be explained by horizontal transfer between natural homing sites by reverse splicing and subsequent spread in populations by endonuclease-dependent homing.

The recent transfer of a homing endonuclease gene

The myxomycete Didymium iridis (isolate Panama 2) contains a mobile group I intron named Dir.S956-1 after position 956 in the nuclear small subunit (SSU) rRNA gene. The intron is efficiently spread through homing by the intron-encoded homing endonuclease I-DirI. Homing endonuclease genes (HEGs) usually spread with their associated introns as a unit, but infrequently also spread independent of introns (or inteins). Clear examples of HEG mobility are however sparse. Here, we provide evidence for the transfer of a HEG into a group I intron named Dir.S956-2 that is inserted into the SSU rDNA of the Costa Rica 8 isolate of D.iridis. Similarities between intron sequences that flank the HEG and rDNA sequences that flank the intron (the homing endonuclease recognition sequence) suggest that the HEG invaded the intron during the recent evolution in a homing-like event. Dir.S956-2 is inserted into the same SSU site as Dir.S956-1. Remarkably, the two group I introns encode distantly related splicing ribozymes with phylogenetically related HEGs inserted on the opposite strands of different peripheral loop regions. The HEGs are both interrupted by small spliceosomal introns that must be removed during RNA maturation.

Site-specific reverse splicing of a HEG-containing group I intron in ribosomal RNA

The wide, but scattered distribution of group I introns in nature is a result of two processes; the vertical inheritance of introns with or without losses, and the occasional transfer of introns across species barriers. Reversal of the group I intron self-splicing reaction, termed reverse splicing, coupled with reverse transcription and genomic integration potentially mediate an RNA-based intron mobility pathway. Compared to the well characterized endonuclease-mediated intron homing, reverse splicing is less specific and represents a likely explanation for many intron transpositions into new genomic sites. However, the frequency and general role of an RNA-based mobility pathway in the spread of natural group I introns is still unclear. We have used the twin-ribozyme intron (Dir.S956-1) from the myxomycete Didymium iridis to test how a mobile group I intron containing a homing endonuclease gene (HEG) selects between potential insertion sites in the small subunit (SSU) rRNA in vitro, in Escherichia coli and in yeast. Surprisingly, the results show a site-specific RNA-based targeting of Dir.S956-1 into its natural (S956) SSU rRNA site. Our results suggest that reverse splicing, in addition to the established endonuclease-mediated homing mechanism, potentially accounts for group I intron spread into the homologous sites of different strains and species.

The molecular evolution and structural organization of self-splicing group I introns at position 516 in nuclear SSU rDNA of myxomycetes

Group I introns are relatively common within nuclear ribosomal DNA of eukaryotic microorganisms, especially in myxomycetes. Introns at position S516 in the small subunit ribosomal RNA gene are particularly common, but have a sporadic occurrence in myxomycetes. Fuligo septica, Badhamia gracilis, and Physarum flavicomum, all members of the family Physaraceae, contain related group IC1 introns at this site. The F. septica intron was studied at the molecular level and found to self-splice as naked RNA and to generate full-length intron RNA circles during incubation. Group I introns at position S516 appear to have a particularly widespread distribution among protists and fungi. Secondary structural analysis of more than 140 S516 group I introns available in the database revealed five different types of organization, including IC1 introns with and without His-Cys homing endonuclease genes, complex twin-ribozyme introns, IE introns, and degenerate group I-like introns. Both intron structural and phylogenetic analyses indicate a multiple origin of the S516 introns during evolution. The myxomycete introns are related to S516 introns in the more distantly related brown algae and Acanthamoeba species. Possible mechanisms of intron transfer both at the RNA- and DNA-levels are discussed in order to explain the observed widespread, but scattered, phylogenetic distribution.

Hydrolytic cleavage by a group I intron ribozyme is dependent on RNA structures not important for splicing

DiGIR2 is the group I splicing-ribozyme of the mobile twin-ribozyme intron Dir.S956-1, present in Didymium nuclear ribosomal DNA. DiGIR2 is responsible for intron excision, exon ligation, 3'-splice site hydrolysis, and full-length intron RNA circle formation. We recently reported that DiGIR2 splicing (intron excision and exon ligation) competes with hydrolysis and subsequent full-length intron circularization. Here we present experimental evidence that hydrolysis at the 3'-splice site in DiGIR2 is dependent on structural elements within the P9 subdomain not involved in splicing. Whereas the GCGA tetra-loop in P9b was found to be important in hydrolytic cleavage, probably due to tertiary RNA-RNA interactions, the P9.2 hairpin structure was found to be essential for hydrolysis. The most important positions in P9.2 include three adenosines in the terminal loop (L9.2) and a consensus kink-turn motif in the proximal stem. We suggest that the L9.2 adenosines and the kink-motif represent key regulatory elements in the splicing and hydrolytic reaction pathways.

The ability to form full-length intron RNA circles is a general property of nuclear group I introns

In addition to splicing, group I intron RNA is capable of an alternative two-step processing pathway that results in the formation of full-length intron circular RNA. The circularization pathway is initiated by hydrolytic cleavage at the 3' splice site and followed by a transesterification reaction in which the intron terminal guanosine attacks the 5' splice site presented in a structure analogous to that of the first step of splicing. The products of the reactions are full-length circular intron and unligated exons. For this reason, the circularization reaction is to the benefit of the intron at the expense of the host. The circularization pathway has distinct structural requirements that differ from those of splicing and appears to be specifically suppressed in vivo. The ability to form full-length circles is found in all types of nuclear group I introns, including those from the Tetrahymena ribosomal DNA. The biological function of the full-length circles is not known, but the fact that the circles contain the entire genetic information of the intron suggests a role in intron mobility.