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

Mobile group I introns at nuclear rDNA position L2066 harbor sense and antisense homing endonuclease genes intervened by spliceosomal introns

Background

Mobile group I introns encode homing endonucleases that confer intron mobility initiated by a double-strand break in the intron-lacking allele at the site of insertion. Nuclear ribosomal DNA of some fungi and protists contain mobile group I introns harboring His-Cys homing endonuclease genes (HEGs). An intriguing question is how protein-coding genes embedded in nuclear ribosomal DNA become expressed. To address this gap of knowledge we analyzed nuclear L2066 group I introns from myxomycetes and ascomycetes.

Results

A total of 34 introns were investigated, including two identified mobile-type introns in myxomycetes with HEGs oriented in sense or antisense directions. Intriguingly, both HEGs are interrupted by spliceosomal introns. The intron in Didymium squamulosum, which harbors an antisense oriented HEG, was investigated in more detail. The group I intron RNA self-splices in vitro, thus generating ligated exons and full-length intron circles. The intron HEG is expressed in vivo in Didymium cells, which involves removal of a 47-nt spliceosomal intron (I-47) and 3′ polyadenylation of the mRNA. The D. squamulosum HEG (lacking the I-47 intron) was over-expressed in E. coli, and the corresponding protein was purified and shown to confer endonuclease activity. The homing endonuclease was shown to cleave an intron-lacking DNA and to produce a pentanucleotide 3′ overhang at the intron insertion site.

Conclusions

The L2066 family of nuclear group I introns all belong to the group IE subclass. The D. squamulosum L2066 intron contains major hallmarks of a true mobile group I intron by encoding a His-Cys homing endonuclease that generates a double-strand break at the DNA insertion site. We propose a potential model to explain how an antisense HEG becomes expressed from a nuclear ribosomal DNA locus.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13100-022-00280-4.

Structural Organization of S516 Group I Introns in Myxomycetes

Group I introns are mobile genetic elements encoding self-splicing ribozymes. Group I introns in nuclear genes are restricted to ribosomal DNA of eukaryotic microorganisms. For example, the myxomycetes, which represent a distinct protist phylum with a unique life strategy, are rich in nucleolar group I introns. We analyzed and compared 75 group I introns at position 516 in the small subunit ribosomal DNA from diverse and distantly related myxomycete taxa. A consensus secondary structure revealed a conserved group IC1 ribozyme core, but with a surprising RNA sequence complexity in the peripheral regions. Five S516 group I introns possess a twintron organization, where a His-Cys homing endonuclease gene insertion was interrupted by a small spliceosomal intron. Eleven S516 introns contained direct repeat arrays with varying lengths of the repeated motif, a varying copy number, and different structural organizations. Phylogenetic analyses of S516 introns and the corresponding host genes revealed a complex inheritance pattern, with both vertical and horizontal transfers. Finally, we reconstructed the evolutionary history of S516 nucleolar group I introns from insertion of mobile-type introns at unoccupied cognate sites, through homing endonuclease gene degradation and loss, and finally to the complete loss of introns. We conclude that myxomycete S516 introns represent a family of genetic elements with surprisingly dynamic structures despite a common function in RNA self-splicing.

Natural circularly permuted group II introns in bacteria produce RNA circles

Group II self-splicing introns are large structured RNAs that remove themselves from transcripts while simultaneously sealing the resulting gaps. Some representatives can subsequently reverse splice into DNA, accounting for their pervasive distribution in bacteria. The catalytically active tertiary structure of each group II intron is assembled from six domains that are arranged in a conserved order. Here, we report structural isomers of group II introns, called CP group II ribozymes, wherein the characteristic order of domains has been altered. Domains five and six, which normally reside at the 3' end of group II introns, instead occupy the 5' end to form circularly permuted variants. These unusual group II intron derivatives are catalytically active and generate large linear branched and small circular RNAs, reaction products that are markedly different from those generated by canonical group II introns. The biological role of CP group II ribozymes is currently unknown.

Genome Evolution from Random Ligation of RNAs of Autocatalytic Sets

The evolutionary origin of the genome remains elusive. Here, I hypothesize that its first iteration, the protogenome, was a multi-ribozyme RNA. It evolved, likely within liposomes (the protocells) forming in dry-wet cycling environments, through the random fusion of ribozymes by a ligase and was amplified by a polymerase. The protogenome thereby linked, in one molecule, the information required to seed the protometabolism (a combination of RNA-based autocatalytic sets) in newly forming protocells. If this combination of autocatalytic sets was evolutionarily advantageous, the protogenome would have amplified in a population of multiplying protocells. It likely was a quasispecies with redundant information, e.g., multiple copies of one ribozyme. As such, new functionalities could evolve, including a genetic code. Once one or more components of the protometabolism were templated by the protogenome (e.g., when a ribozyme was replaced by a protein enzyme), and/or addiction modules evolved, the protometabolism became dependent on the protogenome. Along with increasing fidelity of the RNA polymerase, the protogenome could grow, e.g., by incorporating additional ribozyme domains. Finally, the protogenome could have evolved into a DNA genome with increased stability and storage capacity. I will provide suggestions for experiments to test some aspects of this hypothesis, such as evaluating the ability of ribozyme RNA polymerases to generate random ligation products and testing the catalytic activity of linked ribozyme domains.

Progress toward SHAPE Constrained Computational Prediction of Tertiary Interactions in RNA Structure

As more sequencing data accumulate and novel puzzling genetic regulations are discovered, the need for accurate automated modeling of RNA structure increases. RNA structure modeling from chemical probing experiments has made tremendous progress, however accurately predicting large RNA structures is still challenging for several reasons: RNA are inherently flexible and often adopt many energetically similar structures, which are not reliably distinguished by the available, incomplete thermodynamic model. Moreover, computationally, the problem is aggravated by the relevance of pseudoknots and non-canonical base pairs, which are hardly predicted efficiently. To identify nucleotides involved in pseudoknots and non-canonical interactions, we scrutinized the SHAPE reactivity of each nucleotide of the 188 nt long lariat-capping ribozyme under multiple conditions. Reactivities analyzed in the light of the X-ray structure were shown to report accurately the nucleotide status. Those that seemed paradoxical were rationalized by the nucleotide behavior along molecular dynamic simulations. We show that valuable information on intricate interactions can be deduced from probing with different reagents, and in the presence or absence of Mg2+. Furthermore, probing at increasing temperature was remarkably efficient at pointing to non-canonical interactions and pseudoknot pairings. The possibilities of following such strategies to inform structure modeling software are discussed.

Mechanisms of catalytic RNA molecules

Ribozymes are folded catalytic RNA molecules that perform important biological functions. Since the discovery of the first RNA with catalytic activity in 1982, a large number of ribozymes have been reported. While most catalytic RNA molecules act alone, some RNA-based catalysts, such as RNase P, the ribosome, and the spliceosome, need protein components to perform their functions in the cell. In the last decades, the structure and mechanism of several ribozymes have been studied in detail. Aside from the ribosome, which catalyzes peptide bond formation during protein synthesis, the majority of known ribozymes carry out mostly phosphoryl transfer reactions, notably trans-esterification or hydrolysis reactions. In this review, we describe the main features of the mechanisms of various types of ribozymes that can function with or without the help of proteins to perform their biological functions.

Use of a Lariat Capping Ribozyme to Study Cap Function In Vivo

A lariat cap is a naturally occurring substitute of a conventional mRNA cap and is found in a particular genomic setting in a few eukaryotic microorganisms. It is installed by the lariat capping ribozyme acting in cis. In principle, any RNA molecule in any organism can be equipped with a lariat cap in vivo when expressed downstream of a lariat capping ribozyme. Lariat capping is thus a versatile tool for studying the importance of the 5' end structure of RNA molecules. In this chapter, we present protocols to validate the presence of the lariat cap and measure the efficiency of in vivo cleavage by the lariat capping ribozyme.

Conformational adaptation of UNCG loops upon crowding

If the A-form helix is the major structural motif found in RNA, the loops that cap them constitute the second most important family of motifs. Among those, two are overrepresented, GNRA and UNCG tetraloops. Recent surveys of RNA structures deposited in the PDB show that GNRA and UNCG tetraloops can adopt tertiary folds that are very different from their canonical conformations, characterized by the presence of a U-turn of a Z-turn, respectively. Crystallographic data from both a lariat-capping (LC) ribozyme and a group II intron ribozyme reveal that a given UUCG tetraloop can adopt a distinct fold depending on its structural environment. Specifically, when the crystal packing applies relaxed constraints on the loop, the canonical Z-turn conformation is observed. In contrast, a highly packed environment induces "squashing" of the tetraloop by distorting its sugar-phosphate backbone in a specific way that expels the first and fourth nucleobases out of the loop, and falls in van der Waals distance of the last base pair of the helix, taking the place of the pair formed between the first and fourth residues in Z-turn loops. The biological relevance of our observations is supported by the presence of similarly deformed loops in the highly packed environment of the ribosome and in a complex between a dsRNA and a RNase III. The finding that Z-turn loops change conformation under higher molecular packing suggests that, in addition to their demonstrated role in stabilizing RNA folding, they may contribute to the three-dimensional structure of RNA by mediating tertiary interactions with distal residues.

Dual coding potential of a 2',5'-branched ribonucleotide in DNA

Branchpoints in RNA templates are highly mutagenic, but it is not known yet whether this also applies to branchpoints in DNA templates. Here, we report how nucleic acid polymerases replicate a 2',5'-branched DNA (bDNA) molecule. We constructed long-chained bDNA templates containing a branch guanosine and T7 promoters at both arms by splinted ligation. Quantitative real-time PCR analysis was used to investigate whether a branchpoint blocks DNA synthesis from the two arms in the same manner. We find that the blocking effect of a branchpoint is arm-specific. DNA synthesis from the 2'-arm is more than 20,000-fold decreased, whereas from the 3'-arm only 15-fold. Our sequence analysis of full-length nucleic acid generated by Taq DNA polymerase, Moloney murine leukemia virus reverse transcriptase, and T7 RNA polymerase from the 2'-arm of bDNA shows that the branched guanine has a dual coding potential and can base-pair with cytosine and guanine. We find that branchpoint templating is influenced by the type of the surrounding nucleic acid and is probably modulated by polymerase and RNase H active sites. We show that the branchpoint bypass by the polymerases from the 3'-arm of bDNA is predominantly error-free, indicating that bDNA is not as highly mutagenic as 2',5'-branched RNA.

Lariat capping as a tool to manipulate the 5' end of individual yeast mRNA species in vivo

The 5' cap structure of eukaryotic mRNA is critical for its processing, transport, translation, and stability. The many functions of the cap and the fact that most, if not all, mRNA carries the same type of cap makes it difficult to analyze cap function in vivo at individual steps of gene expression. We have used the lariat capping ribozyme (LCrz) from the myxomycete Didymium to replace the mRNA m⁷G cap of a single reporter mRNA species with a tiny lariat in which the first and the third nucleotide are joined by a 2', 5' phosphodiester bond. We show that the ribozyme functions in vivo in the budding yeast Saccharomyces cerevisiae presumably without cofactors and that lariat capping occurs cotranscriptionally. The lariat-capped reporter mRNA is efficiently exported to the cytoplasm where it is found to be oligoadenylated and evenly distributed. Both the oligoadenylated form and a lariat-capped mRNA with a templated poly(A) tail translates poorly, underlining the critical importance of the m⁷G cap in translation. Finally, the lariat-capped RNA exhibits a threefold longer half-life compared to its m⁷G-capped counterpart, consistent with a key role for the m⁷G cap in mRNA turnover. Our study emphasizes important activities of the m⁷G cap and suggests new utilities of lariat capping as a molecular tool in vivo.

Accumulation of Stable Full-Length Circular Group I Intron RNAs during Heat-Shock

Group I introns in nuclear ribosomal RNA of eukaryotic microorganisms are processed by splicing or circularization. The latter results in formation of full-length circular introns without ligation of the exons and has been proposed to be active in intron mobility. We applied qRT-PCR to estimate the copy number of circular intron RNA from the myxomycete Didymium iridis. In exponentially growing amoebae, the circular introns are nuclear and found in 70 copies per cell. During heat-shock, the circular form is up-regulated to more than 500 copies per cell. The intron harbours two ribozymes that have the potential to linearize the circle. To understand the structural features that maintain circle integrity, we performed chemical and enzymatic probing of the splicing ribozyme combined with molecular modeling to arrive at models of the inactive circular form and its active linear counterpart. We show that the two forms have the same overall structure but differ in key parts, including the catalytic core element P7 and the junctions at which reactions take place. These differences explain the relative stability of the circular species, demonstrate how it is prone to react with a target molecule for circle integration and thus supports the notion that the circular form is a biologically significant molecule possibly with a role in intron mobility.

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.

Survey of chimeric IStron elements in bacterial genomes: multiple molecular symbioses between group I intron ribozymes and DNA transposons

IStrons are chimeric genetic elements composed of a group I intron associated with an insertion sequence (IS). The group I intron is a catalytic RNA providing the IStron with self-splicing ability, which renders IStron insertions harmless to the host genome. The IS element is a DNA transposon conferring mobility, and thus allowing the IStron to spread in genomes. IStrons are therefore a striking example of a molecular symbiosis between unrelated genetic elements endowed with different functions. In this study, we have conducted the first comprehensive survey of IStrons in sequenced genomes that provides insights into the distribution, diversity, origin and evolution of IStrons. We show that IStrons have a restricted phylogenetic distribution limited to two bacterial phyla, the Firmicutes and the Fusobacteria. Nevertheless, diverse IStrons representing two major groups targeting different insertion site motifs were identified. This taken with the finding that while the intron components of all IStrons belong to the same structural class, they are fused to different IS families, indicates that multiple intron–IS symbioses have occurred during evolution. In addition, introns and IS elements related to those that were at the origin of IStrons were also identified.

Molecular characterization of a new member of the lariat capping twin-ribozyme introns

BACKGROUND

Twin-ribozyme introns represent a complex class of mobile group I introns that harbour a lariat capping (LC) ribozyme and a homing endonuclease gene embedded in a conventional self-splicing group I ribozyme (GIR2). Twin-ribozyme introns have so far been confined to nucleolar DNA in Naegleria amoeboflagellates and the myxomycete Didymium iridis.

RESULTS

We characterize structural organization, catalytic properties and molecular evolution of a new twin-ribozyme intron in Allovahlkampfia (Heterolobosea). The intron contains two ribozyme domains with different functions in ribosomal RNA splicing and homing endonuclease mRNA maturation. We found Allovahlkampfia GIR2 to be a typical group IC1 splicing ribozyme responsible for addition of the exogenous guanosine cofactor (exoG), exon ligation and circularization of intron RNA. The Allovahlkampfia LC ribozyme, by contrast, represents an efficient self-cleaving ribozyme that generates a small 2',5' lariat cap at the 5' end of the homing endonuclease mRNA, and thus contributes to intron mobility.

CONCLUSIONS

The discovery of a twin-ribozyme intron in a member of Heterolobosea expands the distribution pattern of LC ribozymes. We identify a putative regulatory RNA element (AP2.1) in the Allovahlkampfia LC ribozyme that involves homing endonuclease mRNA coding sequences as an important structural component.

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.

cis-Acting 5' hammerhead ribozyme optimization for in vitro transcription of highly structured RNAs

RNA-mediated biological processes usually require precise definition of 5' and 3' ends. RNA ends obtained by in vitro transcription using T7 RNA polymerase are often heterogeneous in length and sequence. An efficient strategy to overcome these drawbacks consists of inserting an RNA with known boundaries in between two ribozymes, usually a 5' hammerhead and a 3' hepatitis delta virus ribozymes, that cleave off the desired RNA. In practice, folding of the three RNAs challenges each other, potentially preventing thorough processing. Folding and cleavage of the 5' hammerhead ribozyme relies on a sequence of nucleotides belonging to the central RNA making it more sensitive than the usual 3' hepatitis delta virus ribozyme. The intrinsic stability of the central RNA may thus prevent correct processing of the full transcript. Here, we present a method in which incorporation of a full-length hammerhead ribozyme with a specific tertiary interaction prevents alternative folding with the lariat capping GIR1 ribozyme and enables complete cleavage in the course of the transcription. This strategy may be transposable for in vitro transcription of any highly structured RNA.

Spliceosome twin introns in fungal nuclear transcripts

The spliceosome is an RNA/protein complex, responsible for intron excision from eukaryotic nuclear transcripts. In bacteria, mitochondria and plastids, intron excision does not involve the spliceosome, but occurs through mechanisms dependent on intron RNA secondary and tertiary structure. For group II/III chloroplast introns, "twintrons" (introns within introns) have been described. The excision of the external intron, and thus proper RNA maturation, necessitates prior removal of the internal intron, which interrupts crucial sequences of the former. We have here predicted analogous instances of spliceosomal twintrons ("stwintrons") in filamentous fungi. In two specific cases, where the internal intron interrupts the donor of the external intron after the first or after the second nucleotide, respectively, we show that intermediates with the sequence predicted by the "stwintron" hypothesis, are produced in the splicing process. This implies that two successive rounds of RNA scanning by the spliceosome are necessary to produce the mature mRNA. The phylogenetic distributions of the stwintrons we have identified suggest that they derive from "late" events, subsequent to the appearance of the host intron. They may well not be limited to fungal nuclear transcripts, and their generation and eventual disappearance in the evolutionary process are relevant to hypotheses of intron origin and alternative splicing.

Nuclear group I introns in self-splicing and beyond

Group I introns are a distinct class of RNA self-splicing introns with an ancient origin. All known group I introns present in eukaryote nuclei interrupt functional ribosomal RNA genes located in ribosomal DNA loci. The discovery of the Tetrahymena intron more than 30 years ago has been essential to our understanding of group I intron catalysis, higher-order RNA structure, and RNA folding, but other intron models have provided information about the biological role. Nuclear group I introns appear widespread among eukaryotic microorganisms, and the plasmodial slime molds (myxomycetes) contain an abundance of self-splicing introns. Here, we summarize the main conclusions from previous work on the Tetrahymena intron on RNA self-splicing catalysis as well as more recent work on myxomycete intron biology. Group I introns in myxomycetes that represent different evolutionary stages, biological roles, and functional settings are discussed.

Group I intron ribozymes

Group I intron ribozymes constitute one of the main classes of ribozymes and have been a particularly important model in the discovery of key concepts in RNA biology as well as in the development of new methods. Compared to other ribozyme classes, group I intron ribozymes display considerable variation both in their structure and the reactions they catalyze. The best described pathway is the splicing pathway that results in a spliced out intron and ligated exons. This is paralleled by the circularization pathway that leads to full-length circular intron and un-ligated exons. In addition, the intronic products of these pathways have the potential to integrate into targets and to form various types of circular RNA molecules. Thus, group I intron ribozymes and associated elements found within group I introns is a rich source of biological phenomena. This chapter provides a strategy and protocols for initial characterization of new group I intron ribozymes.

Recurrent insertion of 5'-terminal nucleotides and loss of the branchpoint motif in lineages of group II introns inserted in mitochondrial preribosomal RNAs

A survey of sequence databases revealed 10 instances of subgroup IIB1 mitochondrial ribosomal introns with 1 to 33 additional nucleotides inserted between the 5' exon and the consensus sequence at the intron 5' end. These 10 introns depart further from the IIB1 consensus in their predicted domain VI structure: In contrast to its basal helix and distal GNRA terminal loop, the middle part of domain VI is highly variable and lacks the bulging A that serves as the branchpoint in lariat formation. In vitro experiments using two closely related IIB1 members inserted at the same ribosomal RNA site in the basidiomycete fungi Grifola frondosa and Pycnoporellus fulgens revealed that both ribozymes are capable of efficient self-splicing. However, whereas the Grifola intron was excised predominantly as a lariat, the Pycnoporellus intron, which possesses six additional nucleotides at the 5' end, yielded only linear products, consistent with its predicted domain VI structure. Strikingly, all of the introns with 5' terminal insertions lack the EBS2 exon-binding site. Moreover, several of them are part of the small subset of group II introns that encode potentially functional homing endonucleases of the LAGLIDADG family rather than reverse transcriptases. Such coincidences suggest causal relationships between the shift to DNA-based mobility, the loss of one of the two ribozyme sites for binding the 5' exon, and the exclusive use of hydrolysis to initiate splicing.

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.

Structure-based search reveals hammerhead ribozymes in the human microbiome

Deep sequencing of viral or bacterial nucleic acids monitors the presence and diversity of microbes in select populations and locations. Metagenomic study of mammalian viromes can help trace paths of viral transmissions within or between species. High throughput sequencing of patient and untreated sewage microbiomes showed many sequences with no similarity to genomic sequences of known function or origin. To estimate the distribution of functional RNAs in these microbiomes, we used the hammerhead ribozyme (HHR) motif to search for sequences capable of assuming its three-way junction fold. Although only two of the three possible natural HHR topologies had been known, our analysis revealed highly active ribozymes that terminated in any of the three stems. The most abundant of these are type II HHRs, one of which is the fastest natural cis-acting HHR yet discovered. Altogether, 13 ribozymes were confirmed in vitro, but only one showed sequence similarity to previously described HHRs. Sequences surrounding the ribozymes do not generally show similarity to known genes, except in one case, where a ribozyme is immediately preceded by a bacterial RadC gene. We demonstrate that a structure-based search for a known functional RNA is a powerful tool for analysis of metagenomic datasets, complementing sequence alignments.

DNA as a versatile chemical component for catalysis, encoding, and stereocontrol

DNA (deoxyribonucleic acid) is the genetic material common to all of Earth's organisms. Our biological understanding of DNA is extensive and well-exploited. In recent years, chemists have begun to develop DNA for nonbiological applications in catalysis, encoding, and stereochemical control. This Review summarizes key advances in these three exciting research areas, each of which takes advantage of a different subset of DNA's useful chemical properties.

The tertiary structure of group II introns: implications for biological function and evolution

Group II introns are some of the largest ribozymes in nature, and they are a major source of information about RNA assembly and tertiary structural organization. These introns are of biological significance because they are self-splicing mobile elements that have migrated into diverse genomes and played a major role in the genomic organization and metabolism of most life forms. The tertiary structure of group II introns has been the subject of many phylogenetic, genetic, biochemical and biophysical investigations, all of which are consistent with the recent crystal structure of an intact group IIC intron from the alkaliphilic eubacterium Oceanobacillus iheyensis. The crystal structure reveals that catalytic intron domain V is enfolded within the other intronic domains through an elaborate network of diverse tertiary interactions. Within the folded core, DV adopts an activated conformation that readily binds catalytic metal ions and positions them in a manner appropriate for reaction with nucleic acid targets. The tertiary structure of the group II intron reveals new information on motifs for RNA architectural organization, mechanisms of group II intron catalysis, and the evolutionary relationships among RNA processing systems. Guided by the structure and the wealth of previous genetic and biochemical work, it is now possible to deduce the probable location of DVI and the site of additional domains that contribute to the function of the highly derived group IIB and IIA introns.

Spliceosome discards intermediates via the DEAH box ATPase Prp43p

To promote fidelity in nuclear pre-mRNA splicing, the spliceosome rejects and discards suboptimal substrates that have engaged the spliceosome. Whereas DExD/H box ATPases have been implicated in rejecting suboptimal substrates, the mechanism for discarding suboptimal substrates has remained obscure. Corroborating evidence that suboptimal, mutated lariat intermediates can be exported to the cytoplasm for turnover, we have found that the ribosome can translate mutated lariat intermediates. By glycerol gradient analysis, we have found that the spliceosome can dissociate mutated lariat intermediates in vivo in a manner that requires the DEAH box ATPase Prp43p. Through an in vitro assay, we demonstrate that Prp43p promotes the discard of suboptimal and optimal 5' exon and lariat intermediates indiscriminately. Finally, we demonstrate a requirement for Prp43p in repressing splicing at a cryptic splice site. We propose a model for the fidelity of exon ligation in which the DEAH box ATPase Prp22p slows the flow of suboptimal intermediates through exon ligation and Prp43p generally promotes discard of intermediates, thereby establishing a pathway for turnover of stalled intermediates. Because Prp43p also promotes spliceosome disassembly after exon ligation, this work establishes a parallel between the discard of suboptimal intermediates and the dissociation of a genuine excised intron product.

Exploring RNA structure by integrative molecular modelling

RNA molecular modelling is adequate to rapidly tackle the structure of RNA molecules. With new structured RNAs constituting a central class of cellular regulators discovered every year, the need for swift and reliable modelling methods is more crucial than ever. The pragmatic method based on interactive all-atom molecular modelling relies on the observation that specific structural motifs are recurrently found in RNA sequences. Once identified by a combination of comparative sequence analysis and biochemical data, the motifs composing the secondary structure of a given RNA can be extruded in three dimensions (3D) and used as building blocks assembled manually during a bioinformatic interactive process. Comparing the models to the corresponding crystal structures has validated the method as being powerful to predict the RNA topology and architecture while being less accurate regarding the prediction of base-base interactions. These aspects as well as the necessary steps towards automation will be discussed.

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.

Long noncoding RNAs: implications for antigen receptor diversification

Noncoding RNAs (ncRNAs), both small and large, have recently risen to prominence as surprisingly versatile regulators of gene expression. In fact, eukaryotic transcriptomes are rife with RNAs that do not code for protein, though the majority of these species remains wholly uncharacterized. The functional diversity among the mere handful of validated ncRNAs hints at the vast regulatory potential of these silent biomolecules. Though the act of noncoding transcription and the resultant ncRNAs do not directly produce proteins, they represent powerful means of gene control. Here we survey the accumulating literature on the myriad functions of long ncRNAs and emphasize one curious case of noncoding transcription at antigen receptor loci in lymphocytes.

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.

A structural basis for the recognition of 2'-deoxyguanosine by the purine riboswitch

Riboswitches are noncoding RNA elements that are commonly found in the 5'-untranslated region of bacterial mRNA. Binding of a small-molecule metabolite to the riboswitch aptamer domain guides the folding of the downstream sequence into one of two mutually exclusive secondary structures that directs gene expression. The purine riboswitch family, which regulates aspects of purine biosynthesis and transport, contains three distinct classes that specifically recognize guanine/hypoxanthine, adenine, or 2'-deoxyguanosine (dG). Structural analysis of the guanine and adenine classes revealed a binding pocket that almost completely buries the nucleobase within the core of the folded RNA. Thus, it is somewhat surprising that this family of RNA elements also recognizes dG. We have used a combination of structural and biochemical techniques to understand how the guanine riboswitch could be converted into a dG binder and the structural basis for dG recognition. These studies reveal that a limited number of sequence changes to a guanine-sensing RNA are required to cause a specificity switch from guanine to 2'-deoxyguanosine, and to impart an altered structure for accommodating the additional deoxyribose sugar moiety.

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.

Searching genomes for ribozymes and riboswitches

A discussion of experimental approaches and theoretical difficulties in the identification of ribozymes with novel catalytic functions.

New regulatory RNAs with complex structures have recently been discovered, among them the first catalytic riboswitch, a gene-regulatory RNA sequence with catalytic activity. Here we discuss some of the experimental approaches and theoretical difficulties attached to the identification of new ribozymes in genomes.

The molecular evolution and structural organization of group I introns at position 1389 in nuclear small subunit rDNA of myxomycetes

The number of nuclear group I introns from myxomycetes is rapidly increasing in GenBank as more rDNA sequences from these organisms are being sequenced. They represent an interesting and complex group of intervening sequences because several introns are mobile (or inferred to be mobile) and many contain large and unusual insertions in peripheral loops. Here we describe related group I introns at position 1389 in the small subunit rDNA of representatives from the myxomycete family Didymiaceae. Phylogenetic analyses support a common origin and mainly vertical inheritance of the intron. All S1389 introns from the Didymiaceae belong to the IC1 subclass of nuclear group I introns. The central catalytic core region of about 100 nt appears divergent in sequence composition even though the introns reside in closely related species. Furthermore, unlike the majority of group I introns from myxomycetes the S1389 introns do not self-splice as naked RNA in vitro under standard conditions, consistent with a dependence on host factors for folding or activity. Finally, the myxomycete S1389 introns are exclusively found within the family Didymiaceae, which suggests that this group I intron was acquired after the split between the families Didymiaceae and Physaraceae.

Ribozymes, riboswitches and beyond: regulation of gene expression without proteins

Although various functions of RNA are carried out in conjunction with proteins, some catalytic RNAs, or ribozymes, which contribute to a range of cellular processes, require little or no assistance from proteins. Furthermore, the discovery of metabolite-sensing riboswitches and other types of RNA sensors has revealed RNA-based mechanisms that cells use to regulate gene expression in response to internal and external changes. Structural studies have shown how these RNAs can carry out a range of functions. In addition, the contribution of ribozymes and riboswitches to gene expression is being revealed as far more widespread than was previously appreciated. These findings have implications for understanding how cellular functions might have evolved from RNA-based origins.

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.

RNABC: forward kinematics to reduce all-atom steric clashes in RNA backbone

Although accurate details in RNA structure are of great importance for understanding RNA function, the backbone conformation is difficult to determine, and most existing RNA structures show serious steric clashes (>or= 0.4 A overlap) when hydrogen atoms are taken into account. We have developed a program called RNABC (RNA Backbone Correction) that performs local perturbations to search for alternative conformations that avoid those steric clashes or other local geometry problems. Its input is an all-atom coordinate file for an RNA crystal structure (usually from the MolProbity web service), with problem areas specified. RNABC rebuilds a suite (the unit from sugar to sugar) by anchoring the phosphorus and base positions, which are clearest in crystallographic electron density, and reconstructing the other atoms using forward kinematics. Geometric parameters are constrained within user-specified tolerance of canonical or original values, and torsion angles are constrained to ranges defined through empirical database analyses. Several optimizations reduce the time required to search the many possible conformations. The output results are clustered and presented to the user, who can choose whether to accept one of the alternative conformations. Two test evaluations show the effectiveness of RNABC, first on the S-motifs from 42 RNA structures, and second on the worst problem suites (clusters of bad clashes, or serious sugar pucker outliers) in 25 unrelated RNA structures. Among the 101 S-motifs, 88 had diagnosed problems, and RNABC produced clash-free conformations with acceptable geometry for 71 of those (about 80%). For the 154 worst problem suites, RNABC proposed alternative conformations for 72. All but 8 of those were judged acceptable after examining electron density (where available) and local conformation. Thus, even for these worst cases, nearly half the time RNABC suggested corrections suitable to initiate further crystallographic refinement. The program is available from http://kinemage.biochem.duke.edu .

Mitochondrial introns: a critical view

Although group I and group II introns were discovered more than 25 years ago, they are still difficult to identify. Modeling their RNA structure also remains particularly challenging for organelle sequences, owing to their great diversity. In fact, accelerated evolution in organelles often results in a reduced RNA structure and a loss of autocatalytic splicing and intron mobility. We set out to identify all mitochondrial group I and II introns in published sequences, and, to this end, we developed and applied a new search approach: RNAweasel. On the basis of the results, we focus here on building a comprehensive picture of mitochondrial group I introns, including a modified (reduced) consensus RNA secondary structure and a concise phylogeny-based subclassification.

Group II introns: structure, folding and splicing mechanism

Group II introns are large autocatalytic RNAs found in organellar genomes of plants and lower eukaryotes, as well as in some bacterial genomes. Interestingly, these ribozymes share characteristic traits with both spliceosomal introns and non-LTR retrotransposons and may have a common evolutionary ancestor. Furthermore, group II intron features such as structure, folding and catalytic mechanism differ considerably from those of other large ribozymes, making group II introns an attractive model system to gain novel insights into RNA biology and biochemistry. This review explores recent advances in the structural and mechanistic characterization of group II intron architecture and self-splicing.

Eukaryotic operon-like transcription of functionally related genes in Drosophila

Complex biological processes require coordinated function of many genes. One evolutionary solution to the problem of coordinately expressing functionally related genes in bacteria and nematodes is organization of genes in operons. Surprisingly, eukaryotic operons are considered rare outside the nematode lineage. In Drosophila melanogaster, we found lounge lizard (llz), which encodes a degenerin/ENaC cation channel, cotranscribed with CheB42a, a nonhomologous gene of unknown function residing <100 bp upstream. These two genes were transcribed from a single promoter as one primary transcript and were processed posttranscriptionally to generate individual mRNAs. The mechanism did not involve alternative splicing, and it differed from the trans splicing used in nematode operons. Both genes were expressed in the same tissues, and previous work suggested that both may be involved in courtship behavior. A bioinformatic approach identified numerous additional loci as potential Drosophila operons. These data reveal eukaryotic operon-like transcription of functionally related genes in Drosophila. The results also suggest that operon-based transcription may be more common in eukaryotes than previously appreciated.

A new RNA branching activity: the GIR1 ribozyme

The formation of lariat intermediates during the first step of splicing of group II introns and spliceosomal introns is a well-studied fundamental reaction in molecular biology. Apart from this prominent example, there are surprisingly few occurrences of branched nucleotides or even 2',5'-phosphodiester bonds in biology. We recently described a new ribozyme, the GIR1 branching ribozyme, which catalyzes the formation of a tiny lariat that caps an mRNA. This new example together with work on artificial branching ribozymes and deoxyribozymes shows that branching is facile and points to the possibility that branching reactions could be more prevalent than previously recognized.

Identification of aberrant 2'-5' RNA linkage isomers by pellicular anion exchange chromatography

During chemical RNA synthesis, many undesired products may be formed. In addition to the "n-x" sequences, depurination products, and incompletely deprotected oligonucleotides, linkage isomers may form during condensation and/or deprotection of the synthetic products. Under acidic conditions, bond migration may alter normal 3'-5' diesters to aberrant 2'-5' diesters. This results in isomers that are difficult to identify by MS and LC-MS techniques because the isomers have identical masses. HPLC methods for identification of these isomers have not advanced because the isomers are not expected to exhibit differences in hydrophobicity that allow resolution by reversed-phase columns. Neither are changes in ionic interactions anticipated for these isomers that would allow resolution by ion exchange methods. We observed that chromatography on pellicular anion exchange phases, but not on porous anion exchange phases, completely resolves oligonucleotides with very slight conformation differences (e.g., DNA vs. RNA of identical sequence). Because incorporation of 2'-5' linkages in RNA will alter solution conformation slightly, we considered that this pellicular ion exchanger might also allow resolution of identical RNA sequences harboring aberrant 2'-5' linkages from those lacking aberrant 2'-5' linkages. Using the nonporous DNAPac PA200 column, we demonstrated a chromatographic procedure for resolving synthetic RNA with aberrant linkages from their normally linked counterparts. Under certain conditions, aberrant isomers are not completely resolved from those containing only normal linkages. Therefore, we also developed an independent linkage-confirming method using a 5'-3' exonuclease. This enzyme produces incomplete digestion products during digestion of synthetic RNA containing aberrant 2'-5' linkages, and these are readily resolved by DNAPac PA200 chromatography.

The receptor for branch-site docking within a group II intron active site

The distinguishing feature of group II introns, and the property that links them with spliceosomal catalysis, is their ability to undergo splicing through branching. In this reaction, the 2'-hydroxyl group of a specific adenosine within intron domain 6 serves as the nucleophile for attack on the 5' splice site. We know less about branching than any other feature of group II intron catalysis, largely because the receptor structure for activating the branch site is unknown. Here, we identify the intronic region that binds the branch site of a group IIB intron. Located in domain 1, close to receptors for intron domain 5 and both splice sites, we demonstrate that the branch-site receptor is a functional element required for transesterification. Furthermore, we show that crosslinked branch sites can carry out both steps of splicing, suggesting that the conformational state of the intron core is set early and that it persists throughout the entire splicing process.